Essentially non-hourglass and non-tensile-instability SPH elastic dynamics

Abstract

Since the tension instability was discovered in updated Lagrangian smoothed particle hydrodynamics (ULSPH) at the end of the 20th century, researchers have made considerable efforts to suppress its occurrence. However, up to the present day, this problem has not been fundamentally resolved. In this paper, the concept of hourglass modes is firstly introduced into ULSPH, and the inherent causes of tension instability in elastic dynamics are clarified based on this brand-new perspective. Specifically, we present an essentially non-hourglass formulation by decomposing the shear acceleration with the Laplacian operator, and a comprehensive set of challenging benchmark cases for elastic dynamics is used to showcase that our method can completely eliminate tensile instability by resolving hourglass modes. The present results reveal the true origin of tension instability and challenge the traditional understanding of its sources, i.e., hourglass modes are the real culprit behind inducing this instability in tension zones rather that the tension itself. Furthermore, a time integration scheme known as dual-criteria time stepping is adopted into the simulation of solids for the first time, to significantly enhance computational efficiency.

arXiv comment43 pages 21 figures

References in corpus (7)

• Smoothed Particle Hydrodynamics and Magnetohydrodynamics
• SPHinXsys: an open-source multi-physics and multi-resolution library based on smoothed particle hydrodynamics
• An Hourglass Control Algorithm for Lagrangian Smooth Particle Hydrodynamics
• Dual-criteria time stepping for weakly compressible smoothed particle hydrodynamics
• An integrative smoothed particle hydrodynamics framework for modeling cardiac function
• An hourglass-free formulation for total Lagrangian smoothed particle hydrodynamics
• Level-set based pre-processing algorithm for particle-based methods