TY - JOUR
T1 - A locking-free meshless local Petrov-Galerkin formulation for thick and thin plates
AU - Li, Qiuzhan
AU - Soric, J.
AU - Jarak, T.
AU - Atluri, S. N.
N1 - Funding Information:
This work was supported by ARL and ONR. The Authors thanks Dr. Raju Namburu, and Y.D.S. Rajapakse for helpful comments.
PY - 2005/9/1
Y1 - 2005/9/1
N2 - In this paper, a locking-free meshless local Petrov-Galerkin formulation is presented for shear flexible thick plates, which remains theoretically valid in the thin-plate limit. The kinematics of a three-dimensional solid is used, instead of the conventional plate assumption. The local symmetric weak form is derived for cylindrical shaped local sub-domains. The numerical characteristics of the local symmetric weak form, in the thin plate limit, are discussed. Based on this discussion, the shear locking is theoretically eliminated by changing the two dependent variables in the governing equations. The moving least square interpolation is utilized in the in-plane numerical discretization for all the three displacement components. In the thickness direction, on the other hand, a linear interpolation is used for in-plane displacements, while a hierarchical quadratic interpolation is utilized for the transverse displacement, in order to eliminate the thickness locking. Numerical examples in both the thin plate limit and the thick plate limit are presented, and the results are compared with available analytical solutions.
AB - In this paper, a locking-free meshless local Petrov-Galerkin formulation is presented for shear flexible thick plates, which remains theoretically valid in the thin-plate limit. The kinematics of a three-dimensional solid is used, instead of the conventional plate assumption. The local symmetric weak form is derived for cylindrical shaped local sub-domains. The numerical characteristics of the local symmetric weak form, in the thin plate limit, are discussed. Based on this discussion, the shear locking is theoretically eliminated by changing the two dependent variables in the governing equations. The moving least square interpolation is utilized in the in-plane numerical discretization for all the three displacement components. In the thickness direction, on the other hand, a linear interpolation is used for in-plane displacements, while a hierarchical quadratic interpolation is utilized for the transverse displacement, in order to eliminate the thickness locking. Numerical examples in both the thin plate limit and the thick plate limit are presented, and the results are compared with available analytical solutions.
KW - Local symmetric weak form
KW - Meshless local Petrov-Galerkin methods
KW - Meshless method
KW - Moving least squares interpolation
KW - Shear locking
KW - Solid plate
KW - Thickness locking
UR - http://www.scopus.com/inward/record.url?scp=19044381983&partnerID=8YFLogxK
U2 - 10.1016/j.jcp.2005.02.008
DO - 10.1016/j.jcp.2005.02.008
M3 - Article
AN - SCOPUS:19044381983
SN - 0021-9991
VL - 208
SP - 116
EP - 133
JO - Journal of Computational Physics
JF - Journal of Computational Physics
IS - 1
ER -