Slope stability analysis using smoothed particle hydrodynamics (SPH) method

Abstract

In conventional deformation analysis of geomaterials, the infinitesimal and the finite deformation theories have been widely used. These theories have been successfully implemented in several numerical methods, such as finite element method (FEM). As a result, it is now possible to predict a wide variety of deformation behaviors of geomaterials. However, when dealing with large deformation problems using the framework of the FEM, excess distortion of the FEM mesh may lead to instability of the calculation.In this study, in order to solve large deformation problem of geomaterials, the smoothed particle hydrodynamics (SPH) method is used. The method is a kind of particle method based on the mesh-free Lagrangian scheme, and is one of the promising numerical methods in the field of geotechnical engineering. The method can solve large deformation problems without mesh distortion. Moreover, it can handle the governing equations and existing constitutive models for geomaterials based on a continuum mechanics. Therefore, this method can represent the entire deformation process of a geomaterial from the small strain region to the large deformation region.In this paper, first, basic theory and formulation of the SPH method based on solid mechanics are summarized. Then, the result of a simple calculation is shown to verify the accuracy of the spatial derivatives based on the theory of the SPH method. Also, simulations of simple shear tests of both an elastic and elasto-plastic material are carried out and the obtained results are compared with theoretical solutions. Based on the obtained results, calculation accuracy of the method is discussed. Finally, a series of slope stability analyses are carried out. The numerical results obtained from the SPH method and the safety factors obtained from the Fellenius method are compared. The results indicated that the SPH method is able to express the same tendencies of safety factor obtained from the conventional circular slippage calculations. Moreover, the SPH method can evaluate both the deformation and stability simultaneously. Based on the series of validations and simulations, the effectiveness of the SPH method is discussed from the point of view of geotechnical engineering.

The particle finite element method for transient granular material flow: modelling and validation

“…In recent versions of SPH, many of the inherent drawbacks of the original version have been addressed and solved. Today, SPH is used for modelling in a wide range of engineering applications, including simulation of slope stability analysis and failure [6,58], and for the modelling of granular material flow [5,30,[37][38][39]48,57,63].…”

Free-Surface Flow Simulations with Smoothed Particle Hydrodynamics Method using High-Performance Computing

“…Eqs. (16) and (18). A CUDA kernel, however, is defined as a CUDA function that is set up and executed N times in parallel by N different CUDA threads.…”

“…Since then, it has been used in several research areas, e.g. coastal engineering [3][4][5][6][7], flooding forecast [8][9][10][11], solid body transport [12][13][14][15], soil mechanics [16][17][18][19][20], sediment erosion or entrainment processes [21][22][23][24], fastmoving non-Newtonian flows [25][26][27][28][29][30][31][32][33], flows in porous media [34][35][36], solute transport [37][38][39], turbulent flows [40][41][42] and multiphase flows [43][44][45][46][47], not to mention manifold industrial applications (see, for instance [48][49][50]…”

WCSPH with Limiting Viscosity for Modeling Landslide Hazard at the Slopes of Artificial Reservoir

“…Even if this field is rather explored with SPH method [37][38][39], this work represented the first attempt to apply the Kinetic Theory of Granular Flow to simulate, at the slope of an artificial basin, the fast dynamics of a complex landslide undergoing large deformation and involving a sediment volume comparable to the stored water volume.…”