/MAT/LAW110 (VEGTER)
Block Format Keyword Elastoplastic constitutive law using the interpolated yield criterion of CorusVegter and the hardening law of Vegter, accounting for strain rate dependency and thermal effect.
Format
(1)  (2)  (3)  (4)  (5)  (6)  (7)  (8)  (9)  (10) 

/MAT/LAW110/mat_ID/unit_ID or /MAT/VEGTER/mat_ID/unit_ID  
mat_title  
${\rho}_{i}$  
E  v  Ires  
${I}_{crit}$  TAB_YLD  MAT_Xscale  MAT_Yscale  ${\u0192}_{bi}$  ${\rho}_{bi}^{0}$  
${\sigma}_{yld}^{0}$  $\text{\Delta}{\sigma}_{m}$  $\beta $  $\text{\Omega}$  $\eta $  
${\epsilon}_{0}$  ${\sigma}_{0}^{\star}$  $\text{\Delta}{G}_{0}$  ${\dot{\epsilon}}_{0}$  $m$  
T_{ini}  C_{hard}  F_{cut}  VP  I_{smooth}  TAB_TEMP 
(1)  (2)  (3)  (4)  (5)  (6)  (7)  (8)  (9)  (10) 

${f}_{un}$  R  ${f}_{ps}^{1}$  ${f}_{ps}^{2}$  ${f}_{sh}$ 
(1)  (2)  (3)  (4)  (5)  (6)  (7)  (8)  (9)  (10) 

${f}_{un}$  R  ${f}_{ps}^{1}$  ${\alpha}_{ps}$  ${f}_{sh}$ 
(1)  (2)  (3)  (4)  (5)  (6)  (7)  (8)  (9)  (10) 

${R}_{m}^{0}$  ${R}_{m}^{45}$  ${R}_{m}^{90}$  ${A}_{g}^{0}$  ${A}_{g}^{45}$  
${A}_{g}^{90}$  ${R}^{0}$  ${R}^{45}$  ${R}^{90}$ 
(1)  (2)  (3)  (4)  (5)  (6)  (7)  (8)  (9)  (10) 

${f}_{un}$  R  ${w}_{ps}$  ${w}_{sh}$ 
Definition
Field  Contents  SI Unit Example 

mat_ID  Material identifier. (Integer, maximum 10 digits) 

unit_ID  (Optional) Unit Identifier. (Integer, maximum 10 digits) 

mat_title  Material title. (Character, maximum 100 characters) 

${\rho}_{i}$  Initial density. (Real) 
$\left[\frac{\text{kg}}{{\text{m}}^{\text{3}}}\right]$ 
E  Young‘s modulus. (Real) 
$\left[\text{Pa}\right]$ 
v  Poisson's ratio. (Real) 

Ires  Return mapping algorithm flag.
(Integer) 

${I}_{crit}$  Vegter formulation selection:
(Integer) 

TAB_YLD  Tabulated yield stress – plastic strain  strain rate
function identifier. (Integer) 

MAT_Xscale  X scale factor of the tabulated yield – plastic strain 
strain rate function. Default = 10 (Real) 

MAT_Yscale  Y scale factor of the tabulated yield – plastic strain 
strain rate function. Default = 1.0 (Real) 

${\u0192}_{bi}$  Biaxial scale factor. (Real > 0.0) 

${\rho}_{bi}^{0}$  Biaxial strain rate ratio in the 0 degree direction with
respect to RD. (Real > 0.0) 

${\sigma}_{yld}^{0}$  Initial yield stress. (Real) 
$\left[\text{Pa}\right]$ 
$\text{\Delta}{\sigma}_{m}$  Hardening stress increment. (Real) 
$\left[\text{Pa}\right]$ 
$\beta $  Large strain hardening parameter. (Real) 

$\text{\Omega}$  Small strain hardening parameter. (Real) 

$\eta $  Hardening exponent. (Real) 

${\epsilon}_{0}$  Initial plastic strain. (Real) 

${\sigma}_{0}^{\star}$  Limit dynamic flow stress. (Real) 
$\left[\text{Pa}\right]$ 
$\text{\Delta}{G}_{0}$  Maximum activation enthalpy. (Real) 
$\left[\text{eV}\right]$ 
${\dot{\epsilon}}_{0}$  Limit strain rate for thermal activated
movement. (Real) 
$\text{[Hz]}$ 
$m$  Strain rate behavior exponent. (Real) 

T_{ini}  Initial temperature. (Real) 
$\left[\text{K}\right]$ 
C_{hard}  Hardening coefficient.
(Real) 

F_{cut}  Cutoff frequency for strain rate filtering. Default = 1.0 x 10^{20} (Real) 
$\text{[Hz]}$ 
VP  Strain rate choice flag.
(Integer) 

I_{smooth}  Interpolation type (in case of tabulated yield function).
(Integer) 

TAB_TEMP  Tabulated yield stress – plastic strain – temperature
identifier. (Integer) 

${f}_{un}$  Uniaxial scale factor. (Real > 0.0) 

R  Lankford coefficient. Default = 1.0 (Real > 0.0) 

${f}_{ps}^{1}$  First component of the plane strain scale factor. (Real > 0.0) 

${f}_{ps}^{2}$  Second component of the plane strain scale factor (
${I}_{crit}$
= 1). (Real > 0.0) 

${\alpha}_{ps}$  Average factor to compute the second component of the plane
strain scale factor (
${I}_{crit}$
= 2). Default = 0.5 (Real > 0.0) 

${f}_{sh}$  Shear scale factor. (Real > 0.0) 

${R}_{m}^{i}$  Maximum uniaxial engineering stress for the direction at
$i$
degree with respect to the
rolling direction (RD) (
${I}_{crit}$
= 3). (Real > 0.0) 

${A}_{g}^{i}$  Maximum uniaxial uniform elongation in % for the direction at
$i$
degree with respect to the
rolling direction (RD) (
${I}_{crit}$
= 3). (Real > 0.0) 

${A}_{g}^{i}$  Lankford coefficient for the direction at
$i$
degree with respect to the
rolling direction (RD) (
${I}_{crit}$
= 3). Default = 1.0 (Real > 0.0) 

${w}_{ps}$  Plane strain weight factor (
${I}_{crit}$
= 4). (Real > 0.0) 

${w}_{sh}$  Shear weight factor (
${I}_{crit}$
= 4). (Real > 0.0) 
Example (Steel)
#RADIOSS STARTER
#12345678910
/UNIT/25
Local unit system
Mg mm s
#12345678910
# 2. MATERIALS:
#12345678910
/MAT/LAW110/1/25
Steel: Icrit = 1, example with 3 angles (0Â°, 45Â° and 90Â° to the RD)
# Init. dens.
7.85E9
# E nu
194200.0 0.3
#12345678910
# Icrit TAB_YLD MAT_Xscale MAT_Yscale fBI rhoBI
1 0 0.0 0.0 1.004 0.889
# YLD0 DSIGM BETA OMEGA n
107.1 179.6 0.25 8.07 1.0
# EPS0 SIGS DG0 Deps0 m
0.0 20.0 800 3.61e3 1.0
# TINI C_HARD F_CUT VP Ismooth TAB_TEMP
293.0 0.0 10000.0 1 1 0
# fUN_THETA R_THETA fPS1_THETA fPS2_THETA fSH_THETA
1.021 0.64 1.061 0.5305 0.560
0.987 0.48 1.037 0.5185 0.640
1.009 0.76 1.048 0.5240 0.560
#12345678910
#12345678910
/MAT/LAW110/2/25
Steel :Icrit = 2, example with 3 angles (0Â°, 45Â° and 90Â° to the RD)
# Init. dens.
7.85E9
# E nu
194200.0 0.3
#12345678910
# Icrit TAB_YLD MAT_Xscale MAT_Yscale fBI rhoBI
2 0 0.0 0.0 1.004 0.889
# YLD0 DSIGM BETA OMEGA n
107.1 179.6 0.25 8.07 1.0
# EPS0 SIGS DG0 Deps0 m
0.0 20.0 800 3.61e3 1.0
# TINI C_HARD F_CUT VP Ismooth TAB_TEMP
293.0 0.0 10000.0 1 1 0
# fUN_THETA R_THETA fPS1_THETA ALPS_THETA fSH_THETA
1.021 0.64 1.061 0.5 0.560
0.987 0.48 1.037 0.5 0.640
1.009 0.76 1.048 0.5 0.560
#12345678910
#12345678910
/MAT/LAW110/3/25
Steel
# Init. dens.
7.85E9
# E nu
194200.0 0.3
#12345678910
# Icrit TAB_YLD MAT_Xscale MAT_Yscale fBI rhoBI
3 0 0.0 0.0 1.004 0.889
# YLD0 DSIGM BETA OMEGA n
107.1 179.6 0.25 8.07 1.0
# EPS0 SIGS DG0 Deps0 m
0.0 20.0 800 3.61e3 1.0
# TINI C_HARD F_CUT VP Ismooth TAB_TEMP
293.0 0.0 10000.0 1 1 0
# RM_0 RM_45 RM_90 AG_0 AG_45
408.4 408.4 408.4 20.0 20.0
# AG_90 R_0 R_45 R_90
20.0 0.64 0.48 0.76
#12345678910
#12345678910
/MAT/LAW110/4/25
Steel :example with 3 angles (0Â°, 45Â° and 90Â° to the RD)
# Init. dens.
7.85E9
# E nu
194200.0 0.3
#12345678910
# Icrit TAB_YLD MAT_Xscale MAT_Yscale fBI rhoBI
4 0 0.0 0.0 1.004 0.889
# YLD0 DSIGM BETA OMEGA n
107.1 179.6 0.25 8.07 1.0
# EPS0 SIGS DG0 Deps0 m
0.0 20.0 800 3.61e3 1.0
# TINI C_HARD F_CUT VP Ismooth TAB_TEMP
293.0 0.0 10000.0 1 1 0
# fUN_THETA R_THETA W_PS W_SH
1.021 0.64 0.4125 0.75
0.987 0.48 0.4125 0.75
1.009 0.76 0.4125 0.75
#12345678910
#enddata
Comments
 If
${I}_{crit}$
= 1, the law uses the
classic Vegter yield criterion which is defined by:$$\varphi =\overline{\sigma}{\sigma}_{Y}$$Where,
 $\varphi $
 Yield function.
 ${\sigma}_{Y}$
 Yield stress.
 $\overline{\sigma}$
 Interpolated Vegter equivalent stress.
The Vegter yield locus, under plane stress conditions, is obtained by defining three second order Bezier interpolation curves in the principal stresses space (Figure 1).These curves are used to link four reference points which are measured experimentally at different loading conditions: shear, uniaxial tension, plane strain, and equibiaxial tension. Between two reference points, the yield locus is then defined in:
$$\sigma =\left(\begin{array}{c}{\sigma}_{1}\\ {\sigma}_{2}\end{array}\right)=\overline{\sigma}\left(\begin{array}{c}{f}_{1}\\ {f}_{2}\end{array}\right)=\left(\begin{array}{c}{\sigma}_{1}^{r1}\\ {\sigma}_{2}^{r1}\end{array}\right)+2\mu \left(\begin{array}{c}{\sigma}_{1}^{r1}{\sigma}_{1}^{h}\\ {\sigma}_{2}^{r1}{\sigma}_{2}^{h}\end{array}\right)+{\mu}^{2}\left(\begin{array}{c}{\sigma}_{1}^{r2}+{\sigma}_{1}^{r1}2{\sigma}_{1}^{h}\\ {\sigma}_{2}^{r2}+{\sigma}_{2}^{r1}2{\sigma}_{2}^{h}\end{array}\right)$$Where, ${\sigma}^{r1}$ , ${\sigma}^{r2}$
 Two reference points.
 ${\sigma}^{h}$
 Hinge point, which is computed automatically.
 $\mu $
 Parameter computed at each time step to determine the position on the yield locus.
The four measured reference points are referred to as: ${f}_{sh}$ , ${f}_{un}$ , ${f}_{ps}$ , and ${\u0192}_{bi}$ . The coordinates of these reference points are provided by you.
Where, ${f}_{sh}^{2}$ = ${f}_{sh}^{1}={f}_{sh}$ , ${f}_{un}^{2}=0$ , ${f}_{bi}^{2}={f}_{bi}^{1}={f}_{bi}$ . For a given direction with respect to the rolling direction, only one component must be set for the shear ${f}_{sh}$ , uniaxial tension ${f}_{un}$ , and equibiaxial tension reference point ${\u0192}_{bi}$ . For the planestrain reference point, the two components ${f}_{ps}^{1}$ and ${f}_{ps}^{2}$ must be provided by you.
There is a freedom on the choice of the second coordinate ${f}_{ps}^{2}$ . If it is not set, its value will be taken as an average of the second coordinates of the two neighboring hinge points.
The hinge point located between two reference points is computed using the reference points coordinates, but also the normal to the yield surface in each reference point. Using the normality rule, the normal to the yield surface can be expressed using the strain rate tensor component:
$$n=\left(\begin{array}{c}{n}_{1}\\ {n}_{2}\end{array}\right)=\left(\begin{array}{c}{\dot{\epsilon}}_{1}\\ {\dot{\epsilon}}_{2}\end{array}\right)=\left(\begin{array}{c}1\\ \frac{{\dot{\epsilon}}_{2}}{{\dot{\epsilon}}_{1}}\end{array}\right)=\left(\begin{array}{c}1\\ {\rho}_{r}\end{array}\right)$$Where, ${\rho}^{r}$ is the strain rate ratio in a given reference point. For the shear reference point, ${\rho}_{sh}=1$ . For the planestrain reference point, ${\rho}_{ps}=0$ . For the uniaxial tension reference point, this ratio can be computed from the Lankford coefficient denoted $R$ :
$${\rho}_{un}=\frac{R}{R+1}$$The Lankford coefficient must be provided by you. Finally, the strain rate ratio in the equibiaxial tension ${\u0192}_{bi}$ must also be set.
 If
${I}_{crit}$
= 2, the standard Vegter
criterion is used. It is the same criterion as the classical Vegter, with a
different input card. In this card, the second coordinate of the plane
strain point is computed as a weighted average of the second coordinate of
the two neighboring hinge points
${\sigma}_{h1}$
and
${\sigma}_{h2}$
: $${f}_{ps}^{2}={\sigma}_{h1}^{2}+{\alpha}_{ps}\left({\sigma}_{h2}^{2}{\sigma}_{h1}^{2}\right)$$
 If ${I}_{crit}$ = 3, the Vegter 2017 yield criterion is used. It is the same criterion as the classical Vegter, with a different input card. In this card, all parameters are determined from the maximum uniaxial engineering stress ${R}_{m}^{i}$ , the maximum uniform elongation ${A}_{g}^{i}$ , and the Lankford coefficient ${R}^{i}$ . These three parameters must be given for three directions at an angle $i$ degrees with respect to the rolling direction RD: 0 degrees, 45 degrees, and 90 degrees. For this criterion, the number of angles is locked at 3.
 If
${I}_{crit}$
= 4, the law uses the
simplified Vegter Lite yield criterion which is defined by:$$\varphi =\overline{\sigma}{\sigma}_{Y}$$Where,
 $\varphi $
 Yield function.
 ${\sigma}_{Y}$
 Yield stress.
 $\overline{\sigma}$
 Interpolated Vegter equivalent stress.
The Vegter yield locus, under plane stress conditions, is obtained by defining two second order Nurbs interpolation curves in the principal stresses space. These curves are used to link three reference points which are measured experimentally at different loading conditions: uniaxial compression, uniaxial tension, and equibiaxial tension. Between two reference points, the yield locus is then defined in:
$$\sigma =\left(\begin{array}{c}{\sigma}_{1}\\ {\sigma}_{2}\end{array}\right)=\overline{\sigma}\left(\begin{array}{c}{f}_{1}\\ {f}_{2}\end{array}\right)=\frac{{\left(1\mu \right)}^{2}\left(\begin{array}{c}{\sigma}_{1}^{r1}\\ {\sigma}_{2}^{r1}\end{array}\right)+2\mu \left(1\mu \right){w}_{h}\left(\begin{array}{c}{\sigma}_{1}^{h}\\ {\sigma}_{2}^{h}\end{array}\right)+{\mu}^{2}\left(\begin{array}{c}{\sigma}_{1}^{r2}\\ {\sigma}_{2}^{r2}\end{array}\right)}{{\left(1\mu \right)}^{2}+2\mu \left(1\mu \right){w}_{h}+{\mu}^{2}}$$$$0\le \mu \le 1,\text{}{w}_{h}\ge \text{0}$$Where, ${\sigma}^{r1}$ , ${\sigma}^{r2}$
 Two reference points.
 ${\sigma}^{h}$
 Hinge point, which is computed automatically.
 $\mu $
 Parameter computed at each time step to determine the position on the yield locus.
 ${w}_{h}$
 Weight factor associated with the hinge point.
In Figure 2, yield criterion is defined by the classic Vegter formulation in planestress ( ${I}_{crit}$ = 2).
As for the classical Vegter formulation, you must set the parameters ${f}_{un}$ , ${\u0192}_{bi}$ , $R$ , and ${\rho}_{bi}$ . In this simplified formulation, the two hinge points are in pure shear and plane strain condition. The two associated weight factors ${w}_{sh}$ and ${w}_{ps}$ must also be defined.
 To take into account the anisotropy, you can define a set of parameters for
several
${N}_{angle}$
directions at different angle $\varphi $
with respect to the rolling direction RD. These directions must be
equally distributed between 0 and
$\frac{\pi}{2}$
. For all loading directions that are located
between the given directions, the Vegter model proposes to use a Fourier
series interpolation.
Not all parameters must be defined for each direction. The parameter ${\u0192}_{bi}$ is uniform among all the directions. The strain rate ratio ${\rho}_{bi}$ in the direction 0 (rolling direction) is sufficient to determine this ratio for all other directions:
$${\rho}_{bi}\left(\theta \right)=\frac{\left({\rho}_{bi}+1\right)+\left({\rho}_{bi}1\right)\mathrm{cos}2\theta}{\left({\rho}_{bi}+1\right)\left({\rho}_{bi}1\right)\mathrm{cos}2\theta}$$For all other parameters, the following Fourier interpolation is used:
$${f}_{un}\left(\theta \right)={\displaystyle \sum}_{m=0}^{{N}_{angle}1}{\text{\psi}}_{\text{un}}^{m}\text{cos}\left(2m\theta \right)$$In the case of the classic, standard, and 2017 Vegter formulation ( ${I}_{crit}$ = 1, 2, 3), a similar operation is used.
$${f}_{sh}^{2}\left(\theta \right)={f}_{sh}^{1}\left(\frac{\pi}{2}\theta \right)$$In the case of the simplified Vegter Lite formulation ( ${I}_{crit}$ = 4), a similar interpolation is used on $R$ , ${w}_{sh}$ , and ${w}_{ps}$ . In this case:
$${w}_{sh}\left(\theta \right)={w}_{sh}\left(\frac{\pi}{2}\theta \right)$$  Parameters must define a convex yield locus, otherwise the simulation might be unstable. That is why a data check is present in the Radioss Starter to check the yield criterion convexity.
 When only one direction is defined, ${N}_{angle}$ = 1, the material law is isotropic and must be used with /PROP/TYPE1 (isotropic shell). If more than 1 direction is defined, ${N}_{angle}$ > 1 or ${I}_{crit}$ = 3, the material law becomes anisotropic and must be used with /PROP/TYPE9 (orthotropic shells).
 The yield stress
${\sigma}_{Y}$
is defined by the following expression:$${\sigma}_{Y}={\sigma}_{0}+\text{\Delta}{\sigma}_{m}\left[\beta \left({\epsilon}_{p}+{\epsilon}_{0}\right)+{\left(1{e}^{\text{\Omega}\left({\epsilon}_{p}+{\epsilon}_{0}\right)}\right)}^{n}\right]+{\sigma}_{0}^{\star}{\left[1+\frac{kT}{\text{\Delta}{G}_{0}}\mathrm{ln}\left(\frac{\dot{\epsilon}}{\dot{{\epsilon}_{0}}}\right)\right]}^{m}$$Where,
 ${\sigma}_{0}$
 Initial yield stress.
 $\text{\Delta}{\sigma}_{m}$
 Hardening stress increment.
 $\beta $
 Large strain hardening parameter.
 ${\epsilon}_{0}$
 Initial plastic strain.
 $\text{\Omega}$
 Small strain hardening parameter.
 $\eta $
 Hardening exponent.
 ${\sigma}_{0}^{\star}$
 Limit dynamic flow stress.
 $k$
 Boltzmann constant.
 $T$
 Temperature.
 $\text{\Delta}{G}_{0}$
 Maximum activation enthalpy.
 ${\dot{\epsilon}}_{0}$
 Inviscid limit strain rate.
 $m$
 Strain rate dependence exponent.
 If T_{ini} is defined, the temperature will be constant at the T_{ini} value over the simulation. This will enable a constant strain rate
dependency (Figure 4a).If you want to compute a temperature evolution to take into account the thermal activation effect (slight increase of strain rate dependency following the increase of the temperature), the option /HEAT/MAT must be defined (Figure 4a).
 The strain rate
$\dot{\epsilon}$
computation depends on the value of the flag
$VP$
:
 If $VP$ = 1: the plastic strain rate is used.
 If $VP$ = 2: the total strain rate is used.
 If $VP$ = 3: the deviatoric strain rate is used.
In all cases the strain rate computation includes a filtering, using the cutting off frequency ${F}_{cut}$ , which is userdefined, such that:
$${\dot{\epsilon}}_{f}=\alpha {\dot{\epsilon}}^{n}+\left(1\alpha \right){\dot{\epsilon}}^{n1}$$With $\alpha =2\pi {F}_{cut}\text{\Delta}t$ .
 If you want to use a tabulated hardening yield stress, the ID of the
tabulated function TAB_YLD must be defined. This table
can be used to define several yield stress evolutions with the plastic
strain, at several strain rates. Two scale factor can be also defined in the
X and Y directions. In this case, the hardening parameters
${\sigma}_{0}$
,
$\text{\Delta}{\sigma}_{m}$
,
$\beta $
,
$\text{\Omega}$
, and
$\eta $
are ignored, and the yield stress
becomes:$${\sigma}_{Y}={\text{f}}_{\text{Y}}\left({\epsilon}_{\text{p}},\dot{\epsilon}\right)$$
Furthermore, you can also define a tabulated yield stress with plastic strain and temperature TAB_TEMP, to consider the thermal softening in adiabatic condition. In this case, the yield stress becomes:
$${\sigma}_{Y}={\text{f}}_{\text{Y}}\left({\epsilon}_{\text{p}},\dot{\epsilon}\right)\frac{{\text{f}}_{\text{temp}}\left({\epsilon}_{\text{p}},\text{T}\right)}{{\text{f}}_{\text{temp}}\left({\epsilon}_{\text{p}},{\text{T}}_{\text{ini}}\right)}$$Where, T_{ini}
 Reference temperature.
 $T$
 Actual temperature computed using the /HEAT/MAT option.
 You might also want to use a kinematic hardening by setting coefficient C_{hard}:
 If C_{hard} = 0: the isotropic hardening is used.
 If C_{hard} = 1: the PragerZiegler kinematic hardening is used.
 If 0 ≤ C_{hard} ≤ 1: the model is interpolated between the isotropic and the kinematic hardening.