Examples

All benchmarks, examples and applications cases to be run by Kratos. Note that unit tests are in Kratos repository and NOT here

View the Project on GitHub KratosMultiphysics/Examples

Cylinder on inclined plane 2D - comparison between analytical and numerical solution with MPM

Author: Philip Franz

Kratos version: 9.5

Source files: cylinder_on_inclined_plane_2D

Case Specification

This is a 2D simulation of a cylinder on an inclined plane. A rotating as well as a frictionless sliding behaviour of the cylinder are regarded subsequently. The simulation is set up according to section 4.5.2 of (Iaconeta, 2019). Linear, unstructured, triangular elements with a size of 0.01m are used to initialize the MPs. Three MPs per cell are considered. For the backgroundmesh linear, unstructured, triangular elements with a size of 0.02m are used. However, in contrast to section 4.5.2 of (Iaconeta, 2019), the inclined plane is modelled by a line with unstructured elements with size 0.01m. On that line a non conforming Dirichlet boundary condition is imposed by using the penalty method based on (Chandra et al., 2021).

The following applications of Kratos are used:

The problem geometry as well as the boundary conditions are sketched below. The non conforming boundary condition is respresented by the copper coloured line.

Initial geometry and boundary conditions.

A hyper elastic Neo Hookean Plane strain (2D) constitutive law with unit thickness is considered with the following material parameters:

The time step is 0.001 seconds; the total simulation time is 1.0 seconds. The angle (α) of the inclined plane is 60°. The penalty-factor is 1e13.

The contact between cylinder and inclined plane is modelled with the option “contact” (see line 53, file ProjectParameters_contact.json) in the first and with “slip” in the second case, based on (Chandra et al., 2021). Choosing “contact” leads to a rolling behaviour of the cylinder; “slip” to frictionless sliding.

Results

The analytical and numerical solution for the displacement function of the respective case of the above stated problem are compared afterwards:

Initial geometry and boundary conditions.

The left image displays the rolling cylinder - modelled with option “contact”. The right one shows the sliding cylinder (frictionless) - modelled with option “slip”.

Light         Dark

References