| In this paper,the dynamics problem of the incompressible hyperelastic beam is modeled and analyzed based on the absolute nodal coordinate formulation.The existing small deformation materials and linear elastic constitutive models are no longer suitable for the study of incompressible hyperelastic soft structures.The absolute nodal coordinate method uses slope vector instead of the coordinates of the rotation angle to model the dynamics of nonlinear materials with large rotation and large deformation.The soft structures made of nonlinear incompressible hyperelastic materials have good flexibility and strong environmental adaptability and play an important role in medical surgery,geological exploration,disaster rescue and so on.The absolute node coordinate formulation is used to model the dynamics of the low-order beam element,the improved beam element by selective reduced integration method,the highorder beam element and the solid element.The combined use of the beam element with NeoHookean model,Yeoh model and Arruda-Boyce model is proposed,the incompressibility of materials is guaranteed by penalty function method,to implement dynamic modeling and simulation analysis on the silicone beam.The improved beam element,the high-order beam element and the solid element can restrain the volumetric locking effectively.Numerical simulation results demonstrate effectiveness and accuracy of the dynamic model.The dynamics modeling and simulation of contact collision of silicone beam are carried out.The improved beam element and the high-order beam element are combined with the Arruda-Boyce constitutive model respectively to establish the non-contact dynamic equation.Then,combining the contact model based on penalty function method with the regular Coulomb friction model effectively,the corresponding contact forces matrix and friction matrix are derived,so the dynamic model of contact collision is established.Finally,an example of oblique collision between a cantilever beam and a rigid ball is analyzed,the system dynamics equation is iteratively solved by gauss integral method and generalized α method,and the dynamic characteristics of silicone beam in the contact process were observed. |
