While the Johnson-Cook model predicted the behavior of most materials in the Taylor tests, the model's prediction and test results for OFHC (oxygen free high conductivity) copper did not agree well.
In an approach seeking to improve on Johnson-Cook, Zerilli and Armstrong proposed a more sophisticated constitutive relation obtained through the use of dislocation dynamics.
The effects of strain hardening, strain-rate hardening and thermal softening (based on thermal activation analysis) have been incorporated into the formulation. The effect of grain size has also been included.
The relation has a relatively simple expression and should be applicable to a wide range of fcc (face centered cubic) materials.
A relation for iron has also been developed and is also applicable to other bcc (body centered cubic) materials.
An important point made by Zerilli and Armstrong is that each material structure type (fcc, bcc, hcp) will have its own constitutive behavior, dependent on the dislocation characteristics for that particular structure. For example, a stronger dependence of the plastic yield stress on temperature and strain rate is known to result for bcc metals as compared with fcc metals.
With this model, the yield stress varies depending on strain, strain rate and temperature.
The yield stress is given by:
For fcc metals
For bcc metals:
where
ε = effective plastic strain |
= normalized effective plastic strain rate |
T = temperature (degrees K) |
The parameters Y0, C1, C2, C3, C4, C5 and n are material constants.
The plastic flow algorithm used in this model has an implicit strain rate correction option to reduce high frequency oscillations that are sometimes observed in the yield surface under high strain rates. The strain rate correction algorithm will be at the expense of increased CPU usage.
Note: A specific heat capacity property should be defined to enable the calculation of temperature hence the melting effect.
Table 11.3: Input Data
Name | Symbol | Units | Notes |
---|---|---|---|
Initial Yield Stress | Y0 | Stress | |
Hardening Constant #1 | C1 | Stress | |
Hardening Constant #2 | C2 | Stress | |
Hardening Constant #3 | C3 | None | |
Hardening Constant #4 | C4 | None | |
Hardening Constant #5 | C5 | Stress | |
Hardening Constant n | n | None | |
Reference Strain Rate | None | Units fixed at 1/sec Default = 1.0 | |
Strain Rate Correction | None | Option List: None (Default) Implicit |
Custom results variables available for this model:
Name | Description | Solids | Shells | Beams |
---|---|---|---|---|
EFF_PL_STN | Effective Plastic Strain | Yes | Yes* | Yes* |
EFF_PL_STN_RATE | Effective Plastic Strain Rate | Yes | Yes* | Yes* |
TEMP | Temperature** | Yes | Yes* | Yes* |
SUBL_EPS | Effective sublayer plastic strain | No | Yes | No |
*Resultant value over shell/beam section.
**Temperature will be non-zero only if a specific heat capacity is defined.