VM-LSDYNA-BENCHMARK-007

VM-LSDYNA-BENCHMARK-007
Indentation of a Semi-Rigid Long Cylinder in an Elastic, Perfectly Plastic Half-Space

Overview

Reference: Jackson, R. L. (2018). A solution of rigid perfectly plastic cylindrical indentation in plane strain and comparison to elastic-plastic finite element predictions with hardening. Journal of Applied Mechanics, 85 (2).
Analysis Type(s): Explicit Dynamics with Workbench LS-DYNA
Element Type(s): Solid
Input FilesLink to Input Files Download Page

Test Case

A semi-rigid, long, cylinder is pressed into a half-space (approximated by a parallelogram) having isotropic, elastic, perfectly plastic, mechanical properties. The indentation causes plastic deformation of the half-space. The half-space was modeled with Young's modulus (E) of 200 GPa, Poisson's ratio of 0.3, and a yield strength of 200 MPa. The semi-rigid indenter was modeled with Young's modulus (E) of 2.0E+6 GPa and Poisson's ratio of zero.

Figure 243: Problem Sketch

Problem Sketch

Figure 244: Ansys Workbench Model

Ansys Workbench Model

Material PropertiesGeometric PropertiesLoading

E1 = 200 GPa

ν = 0.3

L = 50 in

B = 2 in

H = 2 in

Disp Y = 3.0E-4 m

Analysis Assumptions and Modeling Notes

The closed form solution using slip line theory is discussed in the reference. Given an indentation depth, the length of the indentation mark on the top free surface of the half-space (along the x-direction) and the force on the indenter can be calculated using a closed form solution as shown in Figure 245: Slip Line Theory and the equation below:

(41)

In the above equation, F is the force on the indenter, R is the radius of the indenter, Sy is the yield strength, L is the thickness of the simulation domain which is also the length of the indenter, and b is the size of the indent along the X direction on the top free surface of the indenter.

Figure 245: Slip Line Theory

Slip Line Theory

For LS-DYNA simulations (Explicit), AUTOMATIC_SURFACE_TO_SURFACE contacts were used with the MORTAR option.

The half-space was discretized in two separate regions. The region away from the indenter was meshed by prescribing an element size of 0.25 mm which was not changed. The indenter and the small 2 x 1 x 0.2 mm region under the indenter were meshed according to the prescribed convergence values.

The load application was defined over 0.01 seconds.

The element used was LS-DYNA Explicit Hex 1 with a mesh size of 0.1 mm.

Results Comparison values are reported for the Windows 11 platform.

Results Comparison

ResultsTargetWorkbench LS-DYNAError (%)

Equivalent Stress Max

eq, MPa)

200.0200.00.0

Normal Stress Max

norm, MPa)

26.9626.71–0.91

Directional Def X1 Max

(u, mm)

0.1690.1700.225

Directional Def X2 Max

(u, mm)

0.007440.007440.00

Shear XY Max

xy, MPa)

757.13746.164–1.448

Total strain intensity

(EPTOINT Max)

0.133560.133560.0

Deformation X Max

(u, mm)

0.168510.16910.2314