TY - JOUR
T1 - Procedure for non-smooth contact for planar flexible beams with cone complementarity problem
AU - Yu, Xinxin
AU - Matikainen, Marko K.
AU - Harish, Ajay B.
AU - Mikkola, Aki
N1 - Funding Information:
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the Academy of Finland (Application No. 299033 for funding of Academy Research Fellow) for supporting Marko K. Matikainen.
Publisher Copyright:
© IMechE 2020.
PY - 2021/6
Y1 - 2021/6
N2 - Contact description plays an important role in modeling of applications involving flexible multibody dynamics. Example of such applications include contact between a belt and pulley, crash-worthiness analysis in aerospace and automotive engineering. Approaches such as the linear complementarity problem (LCP), nonlinear complementarity problem (NCP) and penalty method have been proposed for contact detection and imposition of contact constraints. Contact description within multibody dynamics, however, continues to be a challenging topic, particularly in the case of flexible bodies. This paper describes and compares the use of two contact descriptions in the framework of flexible multibody dynamics; (1) the use of nonlinear cone complementarity approach (CCP) and (2) the penalty method. Both contact models are presented together with a master-slave detection algorithm. The modified form of node-to-node approach presented facilitates creation of pseudo-nodes where gap function can be calculated. This reduces the cumbersome effort of contact search. Since large deformations can be an important phenomenon in flexible multibody applications, beam elements based on the absolute nodal coordinate formulation (ANCF) are implemented in this study. To make a comparison of two approaches, the damping component is included in the penalty method by using the continuous contact model introduced by Hunt and Crossley. Numerical results are based on the simulation of ANCF beam contact with rigid ground, rigid body with an arbitrary shape and pendulum contact. Although kinematic results show a good agreement between both approaches when the coefficient of restitution is zero, the unphysical interpenetration appears in the penalty method. Nonlinear minimization problem solved by CCP approach helps to prevent the penetration during contact event. Furthermore, the proposed contact detection algorithm is proved to be capable of being used for multiple contact between beam and arbitrary shape rigid body; different contact types, such as side-by-side and corner-by-side, can be performed without prediction.
AB - Contact description plays an important role in modeling of applications involving flexible multibody dynamics. Example of such applications include contact between a belt and pulley, crash-worthiness analysis in aerospace and automotive engineering. Approaches such as the linear complementarity problem (LCP), nonlinear complementarity problem (NCP) and penalty method have been proposed for contact detection and imposition of contact constraints. Contact description within multibody dynamics, however, continues to be a challenging topic, particularly in the case of flexible bodies. This paper describes and compares the use of two contact descriptions in the framework of flexible multibody dynamics; (1) the use of nonlinear cone complementarity approach (CCP) and (2) the penalty method. Both contact models are presented together with a master-slave detection algorithm. The modified form of node-to-node approach presented facilitates creation of pseudo-nodes where gap function can be calculated. This reduces the cumbersome effort of contact search. Since large deformations can be an important phenomenon in flexible multibody applications, beam elements based on the absolute nodal coordinate formulation (ANCF) are implemented in this study. To make a comparison of two approaches, the damping component is included in the penalty method by using the continuous contact model introduced by Hunt and Crossley. Numerical results are based on the simulation of ANCF beam contact with rigid ground, rigid body with an arbitrary shape and pendulum contact. Although kinematic results show a good agreement between both approaches when the coefficient of restitution is zero, the unphysical interpenetration appears in the penalty method. Nonlinear minimization problem solved by CCP approach helps to prevent the penetration during contact event. Furthermore, the proposed contact detection algorithm is proved to be capable of being used for multiple contact between beam and arbitrary shape rigid body; different contact types, such as side-by-side and corner-by-side, can be performed without prediction.
KW - ANCF
KW - cone complementarity problem
KW - contact mechanics
KW - penalty method
U2 - 10.1177/1464419320957450
DO - 10.1177/1464419320957450
M3 - Article
SN - 1464-4193
VL - 235
SP - 179
EP - 196
JO - Institution of Mechanical Engineers. Proceedings. Part K: Journal of Multi-Body Dynamics
JF - Institution of Mechanical Engineers. Proceedings. Part K: Journal of Multi-Body Dynamics
IS - 2
ER -