Computational & Technology Resources
an online resource for computational,
engineering & technology publications 

CivilComp Proceedings
ISSN 17593433 CCP: 88
PROCEEDINGS OF THE NINTH INTERNATIONAL CONFERENCE ON COMPUTATIONAL STRUCTURES TECHNOLOGY Edited by: B.H.V. Topping and M. Papadrakakis
Paper 150
On the Spurious Mechanisms of an EightNode Mindlin Plate Finite Element Model J.E. Abdalla Filho^{1}^{,}^{2}, I.M. Belo^{1} and R.D. Machado^{1}
^{1}Mechanical Engineering, Pontificia Catholic University of Paraná (PUCPR), Curitiba, Brazil
J.E. Abdalla Filho, I.M. Belo, R.D. Machado, "On the Spurious Mechanisms of an EightNode Mindlin Plate Finite Element Model", in B.H.V. Topping, M. Papadrakakis, (Editors), "Proceedings of the Ninth International Conference on Computational Structures Technology", CivilComp Press, Stirlingshire, UK, Paper 150, 2008. doi:10.4203/ccp.88.150
Keywords: Mindlin plate, finite elements, strain gradient notation, shear locking, spurious zero energy modes, parasitic shear.
Summary
An eightnode serendipity plate element based on Mindlin's plate theory is
formulated using strain gradient notation. Strain gradient notation is a physically
interpretable notation which allows for an apriori evaluation of the finite element
modeling capabilities [1]. It is possible to precisely identify spurious terms which
are responsible for locking and afterwards remove those terms to render an
errorfree element.
In the present case, the shear strains polynomial expansions of the element are inspected for the identification of spurious terms associated with shear locking. Several spurious terms are identified and removed. In addition, compatible mode terms, which might be mistaken as spurious, are clearly identified by recognizing that they comprise compatibility equations. Those terms are retained in order to avoid the introduction of spurious zero energy modes. Three example problems are solved to validate the procedure and to show the effectiveness of the corrected version of the finite element; namely, the pergola roof [2], a corner supported square plate [3], and a symmetric crossply laminated composite plate [4]. In each case, solutions provided by the model containing the spurious terms are compared to solutions provided by the corrected model. In the pergola roof problem, only transverse displacements are calculated. In the corner supported plate problem, transverse displacements and inplane normal stresses are calculated. In the laminated composite plate problem, inplane normal and transverse shear stresses are calculated. The solutions are displayed via curves which clearly demonstrate the differences in stiffness of the model with and without the spurious terms. Those numerical analyses demonstrate that the identified spurious terms are the cause of shear locking. Those analyses also demonstrate the effectiveness of the procedure employed to eliminate the spurious terms as, in general, solutions provided by the corrected model converge monotonically and faster to analytical solutions. A limitation of the current work is that the element has not been designed to assume arbitrary shapes as its formulation does not embed a geometric mapping procedure such as that used in the isoparametric formulation. Nevertheless, it is shown that the element behaves quite well as a rectangular element after elimination of the spurious terms. The use of strain gradient notation may be viewed as advantageous as it allows spurious terms responsible for shear locking to be identified and eliminated apriori. The authors claim to have demonstrated theoretically the sources of shear locking and spurious zero energy modes in the eightnode serendipity plate element, and to have built an efficient firstorder shear deformation theory element to analyze rectangular plates. References
purchase the fulltext of this paper (price £20)
go to the previous paper 
