Factors Influencing Quasistatic Modeling of Deformation and Failure in Rock-Like Solids by the Smoothed Particle Hydrodynamics Method

As a Lagrangian mesh-free numerical method, the Smoothed Particle Hydrodynamics (SPH) method has been traditionally applied for modeling astrophysics, fluid flows and thermal problems, and there has been a growing interest in applying SPH to solid deformation problems. However, the potential of this method for quasistatic analysis of rock-like brittle materials has not been clearly explored. The major aim of this paper is to investigate the effects of key factors in SPH on the load-deformation response of rock-like solids, including variations in the particle approximation theory, the magnitude of the smoothing length and its variable method. Simple uniaxial compression (UC) loading conditions were chosen, and a series of numerical studies were carried out sequentially on an idealized elastic case and an actual test of marble material. Typical results of the axial stress-strain response from infinitesimal to finite deformation as well as the progressive failure process for the marble tests are given and the influences of various factors are discussed. It is found that only provided proper choices of particle momentum equation and the smoothing length parameter, the SPH method is capable for favorably reproducing the deformation and progressive failure evolution in rock-like materials under quasistatic compression loads.