TY - JOUR
T1 - An enrichment of Q4γs plate finite element using incomplete quadratic functions, an assumed energy orthogonality of Bergan's free formulation, and mixed transverse shear strains
AU - Makarim Katili, Andi
AU - Bletzinger, Kai Uwe
AU - Katili, Irwan
N1 - Publisher Copyright:
© 2024 Elsevier Ltd
PY - 2025/1/15
Y1 - 2025/1/15
N2 - This paper introduces a new quadrilateral plate element named DSPM4, which improves upon the previous Q4γs element. The DSPM4 element has twelve DOFs and four temporary DOFs at the mid-sides of the element. The rotation functions βs are modified by adding an incomplete quadratic function to improve the bending performance. An orthogonality condition between the lower and higher-order bending energy is assumed to satisfy the constant bending patch test. The element development involves kinematic and mixed transverse shear strains. In the local normal-tangential (n-s) coordinate system of each side, the mixed transverse shear strains are obtained by combining the constitutive law for shear forces, equilibrium equations, the constitutive law for bending, and curvatures. The relationship between the temporary DOFs and the final DOFs at the corner nodes is established by equating the constant kinematic and mixed transverse shear strains in the direction cosines of each side of the element. As in the Q4γs element, the discrete shear projection method is applied to avoid the shear-locking effect. Numerical analyses are conducted for thin and thick isotropic plates, and the results are compared with those of the DKMQ and MITC4 elements. Numerical results prove that the DSPM4 element is free of shear locking in thin and thick plates, has a proper rank, and provides better convergence behaviour than the Q4γs element. Its performance in static and free vibration analyses of isotropic plates is highly comparable to the DKMQ element.
AB - This paper introduces a new quadrilateral plate element named DSPM4, which improves upon the previous Q4γs element. The DSPM4 element has twelve DOFs and four temporary DOFs at the mid-sides of the element. The rotation functions βs are modified by adding an incomplete quadratic function to improve the bending performance. An orthogonality condition between the lower and higher-order bending energy is assumed to satisfy the constant bending patch test. The element development involves kinematic and mixed transverse shear strains. In the local normal-tangential (n-s) coordinate system of each side, the mixed transverse shear strains are obtained by combining the constitutive law for shear forces, equilibrium equations, the constitutive law for bending, and curvatures. The relationship between the temporary DOFs and the final DOFs at the corner nodes is established by equating the constant kinematic and mixed transverse shear strains in the direction cosines of each side of the element. As in the Q4γs element, the discrete shear projection method is applied to avoid the shear-locking effect. Numerical analyses are conducted for thin and thick isotropic plates, and the results are compared with those of the DKMQ and MITC4 elements. Numerical results prove that the DSPM4 element is free of shear locking in thin and thick plates, has a proper rank, and provides better convergence behaviour than the Q4γs element. Its performance in static and free vibration analyses of isotropic plates is highly comparable to the DKMQ element.
KW - Discrete shear constraints
KW - Discrete shear projection method
KW - Energy-orthogonal free formulation
KW - Kinematic and mixed transverse shear strains
KW - Reissner-Mindlin plate theory
KW - Static and free vibration analysis
UR - http://www.scopus.com/inward/record.url?scp=85211478183&partnerID=8YFLogxK
U2 - 10.1016/j.compstruc.2024.107619
DO - 10.1016/j.compstruc.2024.107619
M3 - Article
AN - SCOPUS:85211478183
SN - 0045-7949
VL - 307
JO - Computers and Structures
JF - Computers and Structures
M1 - 107619
ER -