A study of moderately thick quadrilateral plate elements based on the absolute nodal coordinate formulation

Finite element analysis using plate elements based on the absolute nodal coordinate formulation (ANCF) can predict the behaviors of moderately thick plates subject to large deformation. However, the formulation is subject to numerical locking, which compromises results. This study was designed to investigate and develop techniques to prevent or mitigate numerical locking phenomena. Three different ANCF plate element types were examined. The first is the original fully parameterized quadrilateral ANCF plate element. The second is an update to this element that linearly interpolates transverse shear strains to overcome slow convergence due to transverse shear locking. Finally, the third is based on a new higher order ANCF plate element that is being introduced here. The higher order plate element makes it possible to describe a higher than first-order transverse displacement field to prevent Poisson thickness locking. The term “higher order” is used, because some nodal coordinates of the new plate element are defined by higher order derivatives.The performance of each plate element type was tested by (1) solving a comprehensive set of small deformation static problems, (2) carrying out eigenfrequency analyses, and (3) analyzing a typical dynamic scenario. The numerical calculations were made using MATLAB. The results of the static and eigenfrequency analyses were benchmarked to reference solutions provided by the commercially available finite element software ANSYS.The results show that shear locking is strongly dependent on material thickness. Poisson thickness locking is independent of thickness, but strongly depends on the Poisson effect. Poisson thickness locking becomes a problem for both of the fully parameterized element types implemented with full 3-D elasticity. Their converged results differ by about 18 % from the ANSYS results. Corresponding results for the new higher order ANCF plate element agree with the benchmark. ANCF plate elements can describe the trapezoidal mode; therefore, they do not suffer from Poisson locking, a reported problem for fully parameterized ANCF beam elements. For cases with shear deformation loading, shear locking slows solution convergence for models based on either the original fully parameterized plate element or the newly introduced higher order plate element.

[1]  Aki Mikkola,et al.  Development of elastic forces for a large deformation plate element based on the absolute nodal coordinate formulation , 2006 .

[2]  Theodore E. Simos,et al.  ICNAAM 2010: International Conference of Numerical Analysis and Applied Mathematics 2010 , 2010 .

[3]  E. Stein,et al.  A 4-node finite shell element for the implementation of general hyperelastic 3D-elasticity at finite strains , 1996 .

[4]  R. Y. Yakoub,et al.  Three Dimensional Absolute Nodal Coordinate Formulation for Beam Elements: Theory , 2001 .

[5]  S. Timoshenko,et al.  THEORY OF PLATES AND SHELLS , 1959 .

[6]  Jaewook Rhim,et al.  A vectorial approach to computational modelling of beams undergoing finite rotations , 1998 .

[7]  R. D. Wood,et al.  Nonlinear Continuum Mechanics for Finite Element Analysis , 1997 .

[8]  Rita G. Toscano,et al.  A shell element for finite strain analyses: hyperelastic material models , 2007 .

[9]  Erasmo Carrera,et al.  Analysis of thickness locking in classical, refined and mixed multilayered plate theories , 2008 .

[10]  E. Ramm,et al.  Shear deformable shell elements for large strains and rotations , 1997 .

[11]  Ahmed A. Shabana,et al.  Poisson modes and general nonlinear constitutive models in the large displacement analysis of beams , 2007 .

[12]  J. N. Reddy,et al.  Shear deformable beams and plates: relationships with classical solutions , 2000 .

[13]  K. Bathe,et al.  A continuum mechanics based four‐node shell element for general non‐linear analysis , 1984 .

[14]  J. Domínguez,et al.  An Internal Damping Model for the Absolute Nodal Coordinate Formulation , 2005 .

[15]  R. Hauptmann,et al.  `Solid-shell' elements with linear and quadratic shape functions at large deformations with nearly incompressible materials , 2001 .

[16]  A. Shabana,et al.  Three-dimensional absolute nodal co-ordinate formulation: Plate problem , 1997 .

[17]  O. C. Zienkiewicz,et al.  Analysis of thick and thin shell structures by curved finite elements , 1970 .

[18]  Jan B. Jonker,et al.  SPACAR — Computer Program for Dynamic Analysis of Flexible Spatial Mechanisms and Manipulators , 1990 .

[19]  Werner Schiehlen,et al.  Multibody Systems Handbook , 2012 .

[20]  A. Mikkola,et al.  Elimination of high frequencies in the absolute nodal coordinate formulation , 2010 .

[21]  Gennady M. Kulikov,et al.  Equivalent Single-Layer and Layerwise Shell Theories and Rigid-Body Motions—Part I: Foundations , 2005 .

[22]  Aki Mikkola,et al.  Comparison between ANCF and B-spline surfaces , 2013 .

[23]  A. Mikkola,et al.  A geometrically exact beam element based on the absolute nodal coordinate formulation , 2008 .

[24]  Ahmed A. Shabana,et al.  Analysis of Thin Plate Structures Using the Absolute Nodal Coordinate Formulation , 2005 .

[25]  Ahmed A. Shabana,et al.  On the integration of computer aided design and analysis using the finite element absolute nodal coordinate formulation , 2009 .

[26]  Oleg Dmitrochenko,et al.  Generalization of Plate Finite Elements for Absolute Nodal Coordinate Formulation , 2003 .

[27]  Aki Mikkola,et al.  A Non-Incremental Finite Element Procedure for the Analysis of Large Deformation of Plates and Shells in Mechanical System Applications , 2003 .

[28]  A. Schwab,et al.  COMPARISON OF TWO MODERATELY THICK PLATE ELEMENTS BASED ON THE ABSOLUTE NODAL COORDINATE FORMULATION , 2009 .

[29]  K. Bathe.  Finite Element Procedures , 1995 .

[30]  Arend L. Schwab,et al.  Comparison of Three-Dimensional Flexible Beam Elements for Dynamic Analysis: Classical Finite Element Formulation and Absolute Nodal Coordinate Formulation , 2010 .

[31]  Arend L. Schwab,et al.  COMPARISON OF THREE-DIMENSIONAL FLEXIBLE THIN PLATE ELEMENTS FOR MULTIBODY DYNAMIC ANALYSIS: FINITE ELEMENT FORMULATION AND ABSOLUTE NODAL COORDINATE FORMULATION , 2007 .

[32]  J. N. Reddy,et al.  Tensor-based finite element formulation for geometrically nonlinear analysis of shell structures , 2007 .

[33]  A. Shabana,et al.  EFFICIENT INTEGRATION OF THE ELASTIC FORCES AND THIN THREE-DIMENSIONAL BEAM ELEMENTS IN THE ABSOLUTE NODAL COORDINATE FORMULATION , 2005 .

[34]  Aki Mikkola,et al.  The Simplest 3- and 4-Noded Fully-Parameterized ANCF Plate Elements , 2012 .

[35]  Hiroyuki Sugiyama,et al.  Formulation of Three-Dimensional Joint Constraints Using the Absolute Nodal Coordinates , 2003 .

[36]  Aki Mikkola,et al.  Beam Elements with Trapezoidal Cross Section Deformation Modes Based on the Absolute Nodal Coordinate Formulation , 2010 .

[37]  A. Shabana.  Definition of the Slopes and the Finite Element Absolute Nodal Coordinate Formulation , 1997 .