“…In the model tests, the pipes exhibited only slight deformation compared to their displacement, meeting the requirement for implementation of the oneway transient DEM-FEM coupling method [48,49]. As in previous studies [33,46,50,51], this one-way transient coupling algorithm involved computing boundary data (i.e., particle forces on the surface of the structure) using the DEM and subsequently integrating it into the FEM simulation as a condition at each step, without providing feedback from the FEM to the DEM. During the DEM-FEM coupling, the particle forces on the boundary of the pipe (triangular geometry) were collected first in the DEM part, and then the vertex force per boundary triangle was computed as the average of the values corresponding to the surrounding triangles.…”