See All

Coupling Multibody systems and granular theory

The study of some industrial process requires to model a machine interacting with granular media (ie. the tamping process of railway track). Granular simulations model the contact between huge amount of grains. Some contact laws introduce introduce non-penetration constraint between the surfaces of the grains. Such constraints, called unilateral constraints, requires specific numerical process. This specific numerical solution have to guarantee the respect of the kinematical constraints inside the multibody systems.

In Robotran we developed and implemented a strong coupling methodology for unifying our multibody formalism using relative coordinates and a discrete element method based on non-smooth contact dynamics (NSCD). Simulations done with this approach highlight the huge impact of neglecting the tamping machine dynamics when evaluating the tamping process performances.

Keywords: Contact, Granular, Non-Smooth Dynamics, Lmgc90, DEM


An efficient approach to model contact in granular media is to treat them as constraints. However this approach requires specific numerical scheme to solve the contact problem and the non-smooth dynamics due to impacts.

For application with lot of grains, it has been observed that the resolution of the contact constraints at the velocity is sufficient. The small penetration between grains does not prevent to get the correct bulk behavior of the granular material.

On the other side multibody system contains bilateral constraints. Such constraints must be solved exactly at position, velocity and acceleration levels. Small deviations in bilateral constraints quickly lead to incorrect system dynamic.


To summary the solution consist in introducing the bilateral constraints of the multibody system in the contact solver, as detailed in [1]. This guarantee that the contact solution respect the constraints of the multibody system.

The implementation of this solution has been done by deep-coupling Robotran and LMGC90. The latter dedicated to physics simulation of discrete material.

Study case: tamping process of railway tracks

The tamping process is a maintenance operation performed on railway tracks. This operation restores a smooth track geometry. While the rail is lift to the desired position, the ballast under the sleeper is compacted by the tamping machine.

Steps of the tamping process.
The tamping of ballast consists of three main phases: penetration, squeezing and lifting

This process is permitted due to the vibration of the tamping machine. The vibrations induces a viscous-like state in the ballast allowing easy diving and squeezing.

Our simulations shows that neglecting the effect of the granular media on the machine leads to erroneous results. Indeed imposing the kinematic lead to the conclusion that a vibration of 15Hz is sufficient to dive. On the opposite taking into account for the dynamics lead to a required vibration frequency of 30Hz with is in the range of existing machine.

Imposing the kinematic underestimate the required vibration frequency


[1] Docquier, N., Lantsoght, O., Dubois, F. et al. Modelling and simulation of coupled multibody systems and granular media using the non-smooth contact dynamics approach. Multibody Syst Dyn49, 181–202 (2020).