A nine‐node corotational curved quadrilateral shell element for smooth, folded, and multishell structures
Open Access
- 24 August 2018
- journal article
- research article
- Published by Wiley in International Journal for Numerical Methods in Engineering
- Vol. 116 (8) , 570-600
Abstract
Summary: A nine‐node corotational curved quadrilateral shell element with novel treatment for rotation at intersection of folded and multishell structures is presented. The element's corotational reference frame is defined by the two bisectors of the diagonal vectors generated using the four corner nodes and their cross product. In reference frame, the element rigid‐body rotations are excluded in calculating the local nodal variables from the global nodal variables. Rotations are not represented by axial (pseudo) vectors but by components of polar (proper) vectors, of which additivity and commutativity lead to symmetry of the tangent stiffness matrix. In the global coordinate system, the two smallest components of the midsurface normal vector at each node of a smooth shell or at nodes away from the intersection of nonsmooth shells are defined as rotational variables. In addition, of the two nodal orientation vectors at intersections of folded and multishell structures, two smallest components of one vector, together with one smaller component of another vector, are employed as rotational variables, leading to the desired additive property for all nodal variables in a nonlinear incremental solution procedure. In the local coordinate system, the two smallest components of the midsurface normal vector(s) at any node of a smooth shell or in each smooth patch of nonsmooth shell are defined as rotational variables. Different from other existing corotational finite‐element formulations, the resulting element tangent stiffness matrix is symmetric owing to the commutativity of the local nodal variables in calculating the second derivative of strain energy with respect to these nodal variables. To alleviate membrane and shear locking phenomena, the membrane strains and the out‐of‐plane shear strains are replaced with assumed strains, using the Mixed Interpolation of Tensorial Components approach, for obtaining the element tangent stiffness matrices and the internal force vector. Finally, a series of benchmark, challenging smooth, folded, and multishell structures undergoing large displacements and large rotations are analyzed to demonstrate the reliability and computational accuracy of the proposed formulation.Keywords
Funding Information
- National Natural Science Foundation of China (11672266, 11672266)
This publication has 61 references indexed in Scilit:
- Modeling of Composite Shells in 6–Parameter Nonlinear Theory with Drilling Degree of FreedomMechanics of Advanced Materials and Structures, 2011
- 4-Node combined shell element with semi-EAS-ANS strain interpolations in 6-parameter shell theories with drilling degrees of freedomComputational Mechanics, 2008
- A 9-node co-rotational quadrilateral shell elementComputational Mechanics, 2008
- A survey of recent shell finite elementsInternational Journal for Numerical Methods in Engineering, 2000
- Computational aspects of vector‐like parametrization of three‐dimensional finite rotationsInternational Journal for Numerical Methods in Engineering, 1995
- Variational principles and membrane finite elements with drilling rotations for geometrically non-linear elasticityInternational Journal for Numerical Methods in Engineering, 1995
- A unified formulation for triangular and quadrilateral flat shell finite elements with six nodal degrees of freedomCommunications in Applied Numerical Methods, 1991
- Non-conforming modes stabilization of a nine-node stress-resultant degenerated shell element with drilling freedomComputers & Structures, 1991
- Solution method for nonlinear problems with multiple critical pointsAIAA Journal, 1990
- A compatible triangular element including vertex rotations for plane elasticity analysisComputers & Structures, 1984