Frictional contact algorithms in SPH for the simulation of soil–structure interaction

Abstract

SUMMARY Simulation of frictional contact between soils and rigid or deformable structure in the framework of smoothed particle hydrodynamics (SPH) is presented in this study. Two algorithms are implemented into the SPH code to describe contact behavior, where the contact forces are calculated using the law of conservation of momentum based on ideal plastic collision or using the criteria of partial penetrating. In both algorithms, the problem of boundary deficiency inherited from SPH is properly handled so that the particles located at contact boundary can have precise acceleration, which is critical for contact detection. And the movement and rotation of the rigid structure are taken into account so that it is easy to simulate the process of pile driving or movement of a retaining wall in geotechnical engineering analysis. Furthermore, the capability of modeling deformability of a structure during frictional contact simulations broadens the fields of SPH application. In contrast to previous work dealing with contact in SPH, which usually use particle‐to‐particle contact or ignoring sliding between particles and solid structure, the method proposed here is more efficient and accurate, and it is suitable to simulate interaction between soft materials and rigid or deformable structures, which are very common in geotechnical engineering. A number of numerical tests are carried out to verify the accuracy and stability of the proposed algorithms, and their results are compared with analytical solutions or results from finite element method analysis. Good agreement obtained from these comparisons suggests that the proposed algorithms are robust and can be applied to extend the capability of SPH in solving geotechnical problems. Copyright © 2013 John Wiley & Sons, Ltd.

Human swallowing simulation based on videofluorography images using Hamiltonian MPS method

“…A higher value for the penalty coefficient is better as long as the analysis is stable because it means reduced penetration and increased accuracy for the analysis [51]. We used the following formula for the penalty coefficient:…”

A generic approach to modelling flexible confined boundary conditions in SPH and its application

“…Consequently, as a truly meshless method, SPH is a good alternative to MPM for modelling large deformation problems. SPH was originally proposed for astrophysical applications 13,14 ; it has been then advanced to numbers of applications in geomechanics such as large deformation and failure of geomaterials [17][18][19][20][21][22][23][24] , slope failures and landslides, [25][26][27] ) soil-structure interaction [28][29][30][31] , coupled soil-water problems 18,32,33 , and most recently the scale-dependent rock fracture. 34,35 In SPH, velocity and displacement boundary conditions can be enforced through non-slip and free-slip solid boundary conditions by using ghost 36 and virtual 19,37,38 particles.…”

On Modeling Drilling Load in Lunar Regolith Simulant

“…It may suffer from mesh distortion, resulting in an inaccurate solution and failure of computational convergence [1,2]. Although the DEM does not have that limitation, the high computational cost restricts DEM application to small-scale or shortduration simulations [13].…”