# Calibrate Abaqus Johnson-Cook Model using MCalibration

## Introduction

In this tutorial I will show how to use MCalibration to calibrate the rate-dependent Abaqus Johnson-Cook plasticity model with damage initiation and damage evolution to experimental data. In the YouTube video below, I will also show how you can use the calibrated material model in an Abaqus Explicit FE simulation with element deletion.

The yield stress for the Johnson-Cook model is given by the equation:

$$\sigma_y = \displaystyle \left[ A + B \varepsilon_p^n \right] \cdot \left[ 1 + C \cdot \text{ln} \left( \frac{\dot{\varepsilon}_p} {\dot{\varepsilon}_0} \right) \right] \cdot \left[ 1 – \theta^m \right],$$

where A, B, C, $$\dot{\varepsilon}_0$$, n, m are material parameters, $$\varepsilon_p$$ is the plastic strain, and $$\theta$$ the temperature.

In this example I will use the experimental data shown in the figure to the right. I studied this set of experimental data in detail in the following post.

The first step is to calibrate the stress-strain part of the Abaqus Johnson-Cook material model.

The results from the final stress-strain calibration are shown in the following figure.

The next step is to setup the calibration of the damage initiation and damage evolution parameters. This is a multi-step procedure, the first of which is to export the initially calibrated version of the Johnson-Cook material model.

After that we need to setup an material model template in MCalibration. This figure shows the template definition.

Before  running the final material model calibration it is useful to deactivate all load cases that were not run to failure, and to turn on Failure Time from the load case dialog box. This will add additional zero-stress failure data at the end of the experimental data file. This makes it easier to find the material model parameters.

Here are the final calibration results. The model can tune a bit more, but the results already capture the strain-rate dependence of the failure response.

The final material model is listed below.

				
*Material, name=Mat
** Units: [length]=millimeter, [force]=Newton, [time]=seconds, [temperature]=Kelvin
*Density
1e-09
*Elastic
** E, nu
53.1556584454, 0.4
*Plastic, hardening=Johnson Cook
** A, B, n, m, ThetaMelt, ThetaTrans
1.10531318429, 23.203093989, 0.576311807976, 1, 2000, 1000
*Rate Dependent, type=Johnson Cook
** C, epsDot0
0.0749390218307, 1.06412678755
*Damage initiation, criterion=ductile
**[plast strain], [triax], [strain rate]
0.6, ,        0.001
0.805, ,          0.1
*Damage Evolution, type=displacement, softening=linear
** [plastic displacement]
0.02



### Anisotropic Viscoplastic Material Model Calibration

In this article I show how to quickly calibrate an anisotropic viscoplastic material model to uniaxial tension data.

### Abaqus 2022 – A Great Update for Linux Users

Abaqus 2022 contains some really nice improvements for Linux users. This article summarizes the key improvements.

### All about the Mullins Effect

Detailed tutorial on how to: (1) experimentally test the Mullins damage of rubbers and elastomers; (2) calibrate the Ogden-Roxburgh Mullins model; (3) how to distinguish between Mullins damage and viscoelastic effects.