A penalty approach for contact in smoothed particle hydrodynamics

✍ Scribed by R. Vignjevic; J. Campbell

• Publisher: Elsevier Science
• Year: 1999
• Tongue: English
• Weight: 962 KB
• Volume: 23
• Category: Article
• ISSN: 0734-743X
• DOI: 10.1016/s0734-743x(99)00137-2

✦ Synopsis

This paper describes the development and testing of a contact algorithm for SPH. The treatment of contact boundary conditions in SPH has not been adequately addressed, and the development of the normalised smoothing function approach has highlighted the need for correct treatment of boundary conditions. A particle to particle contact algorithm was developed for 2D. The penalty formulation was used to enforce the contact condition, and several equations for the penalty force calculation were considered. The contact algorithm was tested for one and two dimensional problems for the velocity range between 0.2 and 4.0 km/s to determine the best penalty force equation and the best approach for applying the contact force.

The tests showed that the zero-energy mode problem in the SPH method had to be addressed, as contact excited a zero-energy mode that caused non-physical motion of particles. The test results were compared to the DYNA3D results for the same problems. Background to boundary condition treatment in SPH is given in Section 1. In Section 2 the development of a contact algorithm for SPH in 1 and 2D is described.

Testing of the contact algorithm in 1 and 2D is described in Section 3. Conclusions are presented in Section 4.