/MAT/PLAS_ZERIL
Block Format Keyword This law defines an isotropic elasto-plastic material using the Zerilli-Armstrong plasticity model.
Format
(1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | (10) |
---|---|---|---|---|---|---|---|---|---|
/MAT/PLAS_ZERIL/mat_ID/unit_ID | |||||||||
mat_title | |||||||||
ρi | |||||||||
E | v | VP | |||||||
C0 | C5 | n | εmaxp | σmax0 | |||||
C1 | ˙ε0 | ICC | Fsmooth | Fcut | |||||
C3 | C4 | ρCp | Tr |
Definition
Field | Contents | SI Unit Example |
---|---|---|
mat_ID | Material
identifier. (Integer, maximum 10 digits) |
|
unit_ID | Unit identifier. (Integer, maximum 10 digits) |
|
mat_title | Material
title. (Character, maximum 100 characters) |
|
ρi | Initial
density. (Real) |
[kgm3] |
E | Young's
modulus. (Real) |
[Pa] |
v | Poisson's
ratio. (Real) |
|
VP | Formulation for rate effects.
(Integer) |
|
C0 | Plasticity yield
stress. (Real) |
[Pa] |
C5 | Plasticity hardening
parameter. (Real) |
[Pa] |
n | Plasticity hardening
exponent. 5 Default = 1.0 (Real) |
|
εmaxp | Failure plastic
strain. Default = 1030 (Real) |
|
σmax0 | Plasticity maximum
stress. Default = 1030 (Real) |
[Pa] |
C1 | Strain rate formulation
coefficient. (Real) |
[Pa] |
˙ε0 | Reference strain rate
(must be 1 s-1 converted into user's
units). (Real) |
[1s] |
ICC | Strain rate computation
flag. 7
(Integer) |
|
Fsmooth | Smooth strain rate option flag.
(Integer) |
|
Fcut | Cutoff frequency for
strain rate filtering. 8 Default = 1030 (Real) |
[Hz] |
C3 | Temperature effect
coefficient. (Real) |
[1K] |
C4 | Temperature effect coefficient.
(Real) |
[1K] |
ρCp | Specific heat per unit of volume.
(Real) |
[Jm3⋅K] |
Tr | Reference
temperature. Default = 298 K (Real) |
[K] |
Comments
- The Zerilli-Armstrong law is applicable only to shells and solids.
- The equation that describes
stress during plastic deformation is: σ=C0+(C1exp((−C3T+C4Tln(˙ε˙ε0))))+C5εpnWhere,
- εp
- Plastic strain
- ˙ε
- Strain rate
- T
- Temperature
- Yield stress should be strictly positive.
- When
ˉεp
reaches
εmaxp
in one integration point, then based on the element
type:
- Shell elements:
The corresponding shell element is deleted.
- Solid elements:
The deviatoric stress of the corresponding integral point is permanently set to 0; however, the solid element is not deleted.
- Shell elements:
- n must be less than 1.
- If ˙ε0 is 0, there is no strain rate effect.
- ICC is a flag of the strain rate
effect on material maximum stress
σmax
:
σ=σy(1+cln(˙ε˙εo)) σ=σy(1+cln(˙ε˙εo)) σmax=σmax0(1+cln(˙ε˙εo)) σmax=σmax0 Figure 1. - Strain rate filtering input (Fcut) is only available for shell and solid elements.
- The strain rate filtering is used to smooth strain rates.
- Temperature is computed
assuming adiabatic conditions:Τ=Τr+EintρCρ(Volume)
Where, Eint is the internal energy computed by Radioss.
- When the temperature is not initialized using /HEAT/MAT or /INITEMP, the reference temperature (Tr) is also the initial temperature.