Skip to content

Basic DEM Benchmarking

Ferran Arrufat edited this page Sep 8, 2017 · 18 revisions

Benchmark definition

The DEM Benchmarks consist of a set of 9 simple tests which are run every night and whose object is to make sure both that the application performs correctly and that the code did not break after the daily changes. They are the following:

Benchmark 1: Elastic normal impact of two identical spheres

The diagram shows the time evolution of the elastic contact force between two identical elements.

The velocity after the impact depends on the coefficient of restitution value. For COR=1 the module of the initial and final velocities should be the same. Furthermore, velocities should be equal for both spheres.

Results

bench1

Benchmark 2: Elastic normal impact of a sphere against a rigid wall

The diagram shows the time evolution of the elastic contact force between a discrete element and a finite element.

The velocity after the impact depends on the coefficient of restitution value. For COR=1 the module of the initial and final velocities should be the same.

Benchmark 3: Impact of a sphere against a rigid wall with different coefficients of restitution

The diagram shows the effect of different restitution coefficients on the damping ratio.

If total energy of the system is conserved, the restitution coefficient and the damping ratio values should be identical.

Benchmark 4: Oblique impact of a sphere with a rigid plane with constant velocity module and variable incident angles

The diagram shows the tangential restitution coefficient, final angular velocity and rebound angle of the sphere.

Benchmark 5: Oblique impact of a sphere with a rigid plane with constant normal velocity and different angular velocities

The diagram shows the final linear and angular velocities of the sphere.

Benchmark 6: Oblique impact of a sphere with a rigid plane with constant normal velocity and different angular velocities

The diagram shows the final linear and angular velocities of the sphere.

Benchmark 7: Impact of two identical spheres with a constant normal velocity and different angular velocities

The diagram shows the final linear and angular velocities of both spheres.

By symmetry, the tangential final velocity of both spheres should be zero. Additionally, for a coefficient of restitution of 1, there should be no changes in the modules of both linear and angular velocities and their values should conserve symmetry.

Benchmark 8: Impact of two differently sized spheres with a constant normal velocity and variable angular velocities

The diagram shows the final linear and angular velocities of both spheres.

In this case, it is interesting to note that, the bigger and/or denser the second sphere is, the more this test resembles the sphere versus wall simulation.

Benchmark 9: Impact of two identical spheres with a constant normal velocity and different coefficients of restitution

The diagram shows the effect of different restitution coefficients on the damping ratio.

If total energy of the system is conserved, the restitution coefficient and the damping ratio values should be identical.

References Y.C.Chung, J.Y.Ooi. Benchmark tests for verifying discrete element modelling codes at particle impact level (2011).

Project information

Getting Started

Tutorials

Developers

Kratos structure

Conventions

Solvers

Debugging, profiling and testing

HOW TOs

Utilities

Kratos API

Kratos Structural Mechanics API

Clone this wiki locally