A new five‐node locking‐free quadrilateral element based on smoothed FEM for near‐incompressible linear elasticity. (11th August 2014)

• Record Type: Journal Article
• Authors: Zhang, Yiping; Wang, Sizhao; Chan, Dave
• Is Part Of: International journal for numerical methods in engineering. Volume 100:Number 9(2014)
• Page Start: 633
• Page End: 668
• Publication Date: 2014-08-11
• Subjects: Numerical analysis -- Periodicals; Engineering mathematics -- Periodicals; 620.001518
• DOI: 10.1002/nme.4754
• Languages: English
• ISSNs: 0029-5981

Abstract

<abstract abstract-type="main" id="nme4754-abs-0001"> <title>SUMMARY</title> <p id="nme4754-para-0001">The volumetric locking issue is critical in finite element analysis of nearly incompressible problems. In this paper, a new five‐node quadrilateral element (Q5) is proposed by enriching the four‐node quadrilateral element (Q4) with a centroid node to solve the volumetric locking problem in FEM. The cell‐based smoothed FEM is employed with Q5 element (referred as Q5‐CS‐SC4) to soften the stiffness in order to obtain a better solution. To eliminate pressure oscillation in near‐incompressible problems, an edge‐based area‐weighted smoothing scheme incorporated with the cell‐wise divergence‐free Q5 element is carried out (referred as Q5‐EAW), and an adjustable area‐weighted strain smoothing scheme using a parameter <italic>p</italic> is proposed to improve the performance of the Q5 element in dealing with incompressible media (referred as Q5‐<italic>p</italic>EAW). The formulation of Q5‐<italic>p</italic>EAW is a simple combination of Q5‐CS‐SC4 and Q5‐EAW by an adjustable area weight. It can search the exact strain energy of the problem in near‐incompressible cases. We also introduce another node‐based strain smoothing technique (Q5‐NAW) into the domain‐based selective scheme to obtain a Q5‐<italic>p</italic>EAW/NAW model to solve the pressure oscillation, which gives a much smoother pressure solution than Q5‐EAW. Finally, the Q5 element is extended into hexahedral element to develop a nine‐node hexahedral (H9) element shape function by enriching the eight‐node hexahedral element (Q8) with a centroid node. An H9‐Gi/NAW model similar to Q5‐<italic>p</italic>EAW/NAW is proposed by using H9 element to solve the 3D volumetric locking. A series of benchmark problems are provided to demonstrate that the proposed Q5‐<italic>p</italic>EAW/NAW for 2D plane strain problems, and H9‐Gi/NAW model for 3D cases are locking‐free for nearly incompressible problems. Copyright © 2014 John Wiley & Sons, Ltd.</p> </abstract>