Numerical simulation of impact and entrainment behaviors of debris flow by using SPH–DEM–FEM coupling method

Increasing rain levels can easily destabilize and destroy particulate matter in mountainous areas, which can cause natural disasters, such as debris flow and landslides. Constitutive equations and numerical simulation are the theoretical bases for understanding the behavior of these disasters. Thus, this study aimed to investigate the impact of the debris flow and its entrainment behavior on gully bed sediments. We adopted a coupled analysis method based on elastic–plastic constitutive equations by considering the elasto-plasticity of slurry and the elastic characteristics of debris materials. The coupled method consisted of smooth particle hydrodynamic (SPH), discrete element method (DEM), and finite element method (FEM) (SPH–DEM–FEM). SPH particles represented fluid, DEM particles denoted solid immersed in fluid, and FEM elements represented the terrain and structures. The coupling analysis model was used to simulate the coupling contact of solid, liquid, and structures and to describe the entrainment behavior between solid and liquid phases. The model feasibility was verified by comparing the basic simulation results with experimental values of the dam break model and the rotating cylindrical tank model. The coupled model was then combined with the data management and modeling of geographic information system to simulate the 2010 Yohutagawa debris flow event. Finally, we explored the influence of debris shape-related parameters on the debris flow erosion entrainment process.