A Quasi Energy and Momentum Conservative Algorithm Implemented With a Co-Rotational Quadrilateral Shell Element Formulation Using Vectorial Rotational Variables
Zhongxue Li, Xunda Lin, Loc Vu-Quoc, Bassam A. Izzuddin, Haoyan Wei, Jin Xu, Hongtao Qian, Xin Zhuo
{"title":"A Quasi Energy and Momentum Conservative Algorithm Implemented With a Co-Rotational Quadrilateral Shell Element Formulation Using Vectorial Rotational Variables","authors":"Zhongxue Li, Xunda Lin, Loc Vu-Quoc, Bassam A. Izzuddin, Haoyan Wei, Jin Xu, Hongtao Qian, Xin Zhuo","doi":"10.1002/nme.70128","DOIUrl":null,"url":null,"abstract":"<div>\n \n <p>This paper proposes a flexible multi-body dynamics approach for elastic smooth and non-smooth shells undergoing large deformations and large overall motions. The formulation is based on a co-rotational curved quadrilateral shell element employing vectorial rotational variables and a quasi energy and momentum conservation algorithm. Hamilton's principle is adopted to derive the system's dynamic differential equations. When differentiating the kinetic energy functional with respect to time, the part involving the first-order differentiation of vectorial rotational variables with respect to time is integrated into an equivalent load vector, yielding a symmetric equivalent mass matrix. Accelerations, velocities, displacements, body, and surface loads of the generalized midpoints are generated by convex functions. The Newmark scheme is applied to transform the dynamic differential equations of the system into a set of nonlinear equations. Instead of using strains and their first/second derivatives with respect to local nodal variables at the midpoint configuration, the formulation employs time-averaged assumed strains (computed via the MITC method) and their corresponding derivatives evaluated at both ends of the time step for calculating the internal force vector and element tangent stiffness matrix in the local coordinate system. The transformation matrix from local to global coordinates, however, remains computed at the midpoint configuration. This approach ensures near-exact conservation of total energy and exact conservation of linear and angular momenta once the external loads vanish, while also yielding symmetric tangent stiffness matrices in both local and global coordinate systems. Finally, four examples of three smooth shells and one non-smooth shell problems subjected to impulse loads are solved to verify the proposed formulation for flexible multi-body dynamics of shells. It is shown that the results exhibit excellent agreement with those from other references, demonstrating the reliability, accuracy, and long-term stability of the proposed quasi energy and momentum conserving algorithm.</p>\n </div>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 17","pages":""},"PeriodicalIF":2.9000,"publicationDate":"2025-09-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"International Journal for Numerical Methods in Engineering","FirstCategoryId":"5","ListUrlMain":"https://onlinelibrary.wiley.com/doi/10.1002/nme.70128","RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"ENGINEERING, MULTIDISCIPLINARY","Score":null,"Total":0}
引用次数: 0
Abstract
This paper proposes a flexible multi-body dynamics approach for elastic smooth and non-smooth shells undergoing large deformations and large overall motions. The formulation is based on a co-rotational curved quadrilateral shell element employing vectorial rotational variables and a quasi energy and momentum conservation algorithm. Hamilton's principle is adopted to derive the system's dynamic differential equations. When differentiating the kinetic energy functional with respect to time, the part involving the first-order differentiation of vectorial rotational variables with respect to time is integrated into an equivalent load vector, yielding a symmetric equivalent mass matrix. Accelerations, velocities, displacements, body, and surface loads of the generalized midpoints are generated by convex functions. The Newmark scheme is applied to transform the dynamic differential equations of the system into a set of nonlinear equations. Instead of using strains and their first/second derivatives with respect to local nodal variables at the midpoint configuration, the formulation employs time-averaged assumed strains (computed via the MITC method) and their corresponding derivatives evaluated at both ends of the time step for calculating the internal force vector and element tangent stiffness matrix in the local coordinate system. The transformation matrix from local to global coordinates, however, remains computed at the midpoint configuration. This approach ensures near-exact conservation of total energy and exact conservation of linear and angular momenta once the external loads vanish, while also yielding symmetric tangent stiffness matrices in both local and global coordinate systems. Finally, four examples of three smooth shells and one non-smooth shell problems subjected to impulse loads are solved to verify the proposed formulation for flexible multi-body dynamics of shells. It is shown that the results exhibit excellent agreement with those from other references, demonstrating the reliability, accuracy, and long-term stability of the proposed quasi energy and momentum conserving algorithm.
期刊介绍:
The International Journal for Numerical Methods in Engineering publishes original papers describing significant, novel developments in numerical methods that are applicable to engineering problems.
The Journal is known for welcoming contributions in a wide range of areas in computational engineering, including computational issues in model reduction, uncertainty quantification, verification and validation, inverse analysis and stochastic methods, optimisation, element technology, solution techniques and parallel computing, damage and fracture, mechanics at micro and nano-scales, low-speed fluid dynamics, fluid-structure interaction, electromagnetics, coupled diffusion phenomena, and error estimation and mesh generation. It is emphasized that this is by no means an exhaustive list, and particularly papers on multi-scale, multi-physics or multi-disciplinary problems, and on new, emerging topics are welcome.