/MAT/LAW110 (VEGTER)

Block Format Keyword Elasto-plastic constitutive law using the interpolated yield criterion of Corus-Vegter and the hardening law of Vegter, accounting for strain rate dependency and thermal effect.

Format

(1)	(2)	(3)	(5)	(7)	(8)	(9)
/MAT/LAW110/`mat_ID`/`unit_ID` or /MAT/VEGTER/`mat_ID`/`unit_ID`
`mat_title`
$ρ_{i}$
`E`		`v`	`Ires`
$I_{c r i t}$	`TAB_YLD`	`MAT_Xscale`	`MAT_Yscale`	$ƒ_{b i}$		$ρ_{b i}^{0}$
$σ_{y l d}^{0}$		$Δ σ_{m}$	$β$	$Ω$		$η$
$ε_{0}$		$σ_{0}^{⋆}$	$Δ G_{0}$	${\dot{ε}}_{0}$		$m$
`T`_ini		`C`_hard	`F`_cut	`VP`	`I`_smooth	`TAB_TEMP`

I_{c r i t}

= 1: Read

N_{a n g l e}

(Number of experimental angles, at least 1) cards

(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)	(10)
$f_{u n}$		`R`		$f_{p s}^{1}$		$f_{p s}^{2}$		$f_{s h}$

I_{c r i t}

= 2: Read

N_{a n g l e}

(Number of experimental angles, at least 1) cards

(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)	(10)
$f_{u n}$		`R`		$f_{p s}^{1}$		$α_{p s}$		$f_{s h}$

I_{c r i t}

= 3: Read the parameters (

N_{a n g l e}

is locked to 3 in this case, 0, 45, and 90)

(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}$

I_{c r i t}

= 4: Read

N_{a n g l e}

cards

(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)	(10)
$f_{u n}$		`R`		$w_{p s}$		$w_{s h}$

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)
$ρ_{i}$	Initial density. (Real)	$[\frac{kg}{m^{3}}]$
`E`	Young‘s modulus. (Real)	$[Pa]$
`v`	Poisson's ratio. (Real)
`Ires`	Return mapping algorithm flag. = 1 Nice explicit method. = 2 (Default) Newton iterative semi-implicit method – Cutting Plane. (Integer)
$I_{c r i t}$	Vegter formulation selection: = 1 Classic Vegter formulation. = 2 Standard Vegter formulation. = 3 Vegter 2017 formulation. = 4 Simplified Vegter-Lite formulation. (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)
$ƒ_{b i}$	Biaxial scale factor. (Real > 0.0)
$ρ_{b i}^{0}$	Biaxial strain rate ratio in the 0 degree direction with respect to RD. (Real > 0.0)
$σ_{y l d}^{0}$	Initial yield stress. (Real)	$[Pa]$
$Δ σ_{m}$	Hardening stress increment. (Real)	$[Pa]$
$β$	Large strain hardening parameter. (Real)
$Ω$	Small strain hardening parameter. (Real)
$η$	Hardening exponent. (Real)
$ε_{0}$	Initial plastic strain. (Real)
$σ_{0}^{⋆}$	Limit dynamic flow stress. (Real)	$[Pa]$
$Δ G_{0}$	Maximum activation enthalpy. (Real)	$[eV]$
${\dot{ε}}_{0}$	Limit strain rate for thermal activated movement. (Real)	$[Hz]$
$m$	Strain rate behavior exponent. (Real)
`T`_ini	Initial temperature. (Real)	$[K]$
`C`_hard	Hardening coefficient. = 0 Hardening is a full isotropic model. = 1 Hardening uses the kinematic Prager-Ziegler model. = value between 0 and 1 Hardening is interpolated between the two models. (Real)
`F`_cut	Cutoff frequency for strain rate filtering. Default = 1.0 x 10²⁰ (Real)	$[Hz]$
`VP`	Strain rate choice flag. = 1 Strain rate effect on yield stress depends on the plastic strain rate. = 2 (Default) Strain rate effect on yield depends on the total strain rate. = 3 Strain rate effect on yield depends on the deviatoric strain rate. (Integer)
`I`_smooth	Interpolation type (in case of tabulated yield function). = 1 Linear interpolation. = 2 Logarithmic interpolation base 10. = 3 Logarithmic interpolation base n. (Integer)
`TAB_TEMP`	Tabulated yield stress – plastic strain – temperature identifier. (Integer)
$f_{u n}$	Uniaxial scale factor. (Real > 0.0)
`R`	Lankford coefficient. Default = 1.0 (Real > 0.0)
$f_{p s}^{1}$	First component of the plane strain scale factor. (Real > 0.0)
$f_{p s}^{2}$	Second component of the plane strain scale factor ( $I_{c r i t}$ = 1). (Real > 0.0)
$α_{p s}$	Average factor to compute the second component of the plane strain scale factor ( $I_{c r i t}$ = 2). Default = 0.5 (Real > 0.0)
$f_{s h}$	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_{c r i t}$ = 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_{c r i t}$ = 3). (Real > 0.0)
$A_{g}^{i}$	Lankford coefficient for the direction at $i$ degree with respect to the rolling direction (RD) ( $I_{c r i t}$ = 3). Default = 1.0 (Real > 0.0)
$w_{p s}$	Plane strain weight factor ( $I_{c r i t}$ = 4). (Real > 0.0)
$w_{s h}$	Shear weight factor ( $I_{c r i t}$ = 4). (Real > 0.0)

Example (Steel)

#RADIOSS STARTER
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/UNIT/25
Local unit system
                  Mg                  mm                   s
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#-  2. MATERIALS:
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/MAT/LAW110/1/25
Steel: Icrit = 1, example with 3 angles (0Â°, 45Â° and 90Â° to the RD)
#        Init. dens.
             7.85E-9
#                  E                  nu
            194200.0                 0.3
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#    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.61e-3                 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 
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/MAT/LAW110/2/25
Steel :Icrit = 2, example with 3 angles (0Â°, 45Â° and 90Â° to the RD)
#        Init. dens.
             7.85E-9
#                  E                  nu
            194200.0                 0.3
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#    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.61e-3                 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 
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/MAT/LAW110/3/25
Steel
#        Init. dens.
             7.85E-9
#                  E                  nu
            194200.0                 0.3
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#    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.61e-3                 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    
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/MAT/LAW110/4/25
Steel :example with 3 angles (0Â°, 45Â° and 90Â° to the RD)
#        Init. dens.
             7.85E-9
#                  E                  nu
            194200.0                 0.3
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#    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.61e-3                 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 
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#enddata

Comments

If $I_{c r i t}$ = 1, the law uses the classic Vegter yield criterion which is defined by:
(1)
$ϕ = \bar{σ} - σ_{Y}$

Where,

$ϕ$

Yield function.

$σ_{Y}$

Yield stress.

$\bar{σ}$

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).

Figure 1. Yield criterion defined by the classic Vegter formulation in plane-stress (I_crit=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:(2)
$σ = (\begin{matrix} σ_{1} \\ σ_{2} \end{matrix}) = \bar{σ} (\begin{matrix} f_{1} \\ f_{2} \end{matrix}) = (\begin{matrix} σ_{1}^{r 1} \\ σ_{2}^{r 1} \end{matrix}) + 2 μ (\begin{matrix} σ_{1}^{r 1} - σ_{1}^{h} \\ σ_{2}^{r 1} - σ_{2}^{h} \end{matrix}) + μ^{2} (\begin{matrix} σ_{1}^{r 2} + σ_{1}^{r 1} - 2 σ_{1}^{h} \\ σ_{2}^{r 2} + σ_{2}^{r 1} - 2 σ_{2}^{h} \end{matrix})$

Where,

$σ^{r 1}$ , $σ^{r 2}$

Two reference points.

$σ^{h}$

Hinge point, which is computed automatically.

$μ$

Parameter computed at each time step to determine the position on the yield locus.

The four measured reference points are referred to as: $f_{s h}$ , $f_{u n}$ , $f_{p s}$ , and $ƒ_{b i}$ . The coordinates of these reference points are provided by you.

Where, $f_{s h}^{2}$ = $- f_{s h}^{1} = - f_{s h}$ , $f_{u n}^{2} = 0$ , $f_{b i}^{2} = f_{b i}^{1} = f_{b i}$ . For a given direction with respect to the rolling direction, only one component must be set for the shear $f_{s h}$ , uniaxial tension $f_{u n}$ , and equibiaxial tension reference point $ƒ_{b i}$ . For the plane-strain reference point, the two components $f_{p s}^{1}$ and $f_{p s}^{2}$ must be provided by you.

There is a freedom on the choice of the second coordinate $f_{p s}^{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:(3)
$n = (\begin{matrix} n_{1} \\ n_{2} \end{matrix}) = (\begin{matrix} {\dot{ε}}_{1} \\ {\dot{ε}}_{2} \end{matrix}) = (\begin{matrix} 1 \\ \frac{{\dot{ε}}_{2}}{{\dot{ε}}_{1}} \end{matrix}) = (\begin{matrix} 1 \\ ρ_{r} \end{matrix})$

Where, $ρ^{r}$ is the strain rate ratio in a given reference point. For the shear reference point, $ρ_{s h} = - 1$ . For the plane-strain reference point, $ρ_{p s} = 0$ . For the uniaxial tension reference point, this ratio can be computed from the Lankford coefficient denoted $R$ :(4)
$ρ_{u n} = \frac{- R}{R + 1}$

The Lankford coefficient must be provided by you. Finally, the strain rate ratio in the equibiaxial tension $ƒ_{b i}$ must also be set.
If $I_{c r i t}$ = 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 $σ_{h 1}$ and $σ_{h 2}$ :
(5)
$f_{p s}^{2} = σ_{h 1}^{2} + α_{p s} (σ_{h 2}^{2} - σ_{h 1}^{2})$
If $I_{c r i t}$ = 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_{c r i t}$ = 4, the law uses the simplified Vegter Lite yield criterion which is defined by:
(6)
$ϕ = \bar{σ} - σ_{Y}$

Where,

$ϕ$

Yield function.

$σ_{Y}$

Yield stress.

$\bar{σ}$

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:(7)
$σ = (\begin{matrix} σ_{1} \\ σ_{2} \end{matrix}) = \bar{σ} (\begin{matrix} f_{1} \\ f_{2} \end{matrix}) = \frac{{(1 - μ)}^{2} (\begin{matrix} σ_{1}^{r 1} \\ σ_{2}^{r 1} \end{matrix}) + 2 μ (1 - μ) w_{h} (\begin{matrix} σ_{1}^{h} \\ σ_{2}^{h} \end{matrix}) + μ^{2} (\begin{matrix} σ_{1}^{r 2} \\ σ_{2}^{r 2} \end{matrix})}{{(1 - μ)}^{2} + 2 μ (1 - μ) w_{h} + μ^{2}}$
(8)
$0 \leq μ \leq 1, w_{h} \geq 0$

Where,

$σ^{r 1}$ , $σ^{r 2}$

Two reference points.

$σ^{h}$

Hinge point, which is computed automatically.

$μ$

Parameter computed at each time step to determine the position on the yield locus.

$w_{h}$

Weight factor associated with the hinge point.

Figure 2. Yield criterion defined by the classic Vegter formulation in plane-stress (I_crit=2)

In Figure 2, yield criterion is defined by the classic Vegter formulation in plane-stress ( $I_{c r i t}$ = 2).

As for the classical Vegter formulation, you must set the parameters $f_{u n}$ , $ƒ_{b i}$ , $R$ , and $ρ_{b i}$ . In this simplified formulation, the two hinge points are in pure shear and plane strain condition. The two associated weight factors $w_{s h}$ and $w_{p s}$ must also be defined.
To take into account the anisotropy, you can define a set of parameters for several $N_{a n g l e}$ directions at different angle $ϕ$ with respect to the rolling direction RD. These directions must be equally distributed between 0 and $\frac{π}{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 $ƒ_{b i}$ is uniform among all the directions. The strain rate ratio $ρ_{b i}$ in the direction 0 (rolling direction) is sufficient to determine this ratio for all other directions:(9)
$ρ_{b i} (θ) = \frac{(ρ_{b i} + 1) + (ρ_{b i} - 1) \cos 2 θ}{(ρ_{b i} + 1) - (ρ_{b i} - 1) \cos 2 θ}$

For all other parameters, the following Fourier interpolation is used:(10)
$f_{u n} (θ) = \sum_{m = 0}^{N_{a n g l e} - 1} ψ_{un}^{m} cos (2 m θ)$

In the case of the classic, standard, and 2017 Vegter formulation ( $I_{c r i t}$ = 1, 2, 3), a similar operation is used. (11)
$f_{s h}^{2} (θ) = - f_{s h}^{1} (\frac{π}{2} - θ)$

In the case of the simplified Vegter Lite formulation ( $I_{c r i t}$ = 4), a similar interpolation is used on $R$ , $w_{s h}$ , and $w_{p s}$ . In this case:(12)
$w_{s h} (θ) = w_{s h} (\frac{π}{2} - θ)$
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_{a n g l e}$ = 1, the material law is isotropic and must be used with /PROP/TYPE1 (isotropic shell). If more than 1 direction is defined, $N_{a n g l e}$ > 1 or $I_{c r i t}$ = 3, the material law becomes anisotropic and must be used with /PROP/TYPE9 (orthotropic shells).
The yield stress $σ_{Y}$ is defined by the following expression:
(13)
$σ_{Y} = σ_{0} + Δ σ_{m} [β (ε_{p} + ε_{0}) + {(1 - e^{- Ω (ε_{p} + ε_{0})})}^{n}] + σ_{0}^{⋆} {[1 + \frac{k T}{Δ G_{0}} \ln (\frac{\dot{ε}}{\dot{ε_{0}}})]}^{m}$

Where,

$σ_{0}$

Initial yield stress.

$Δ σ_{m}$

Hardening stress increment.

$β$

Large strain hardening parameter.

$ε_{0}$

Initial plastic strain.

$Ω$

Small strain hardening parameter.

$η$

Hardening exponent.

$σ_{0}^{⋆}$

Limit dynamic flow stress.

$k$

Boltzmann constant.

$T$

Temperature.

$Δ G_{0}$

Maximum activation enthalpy.

${\dot{ε}}_{0}$

Inviscid limit strain rate.

$m$

Strain rate dependence exponent.

The parameter, $Ω$ is called small strain hardening parameter, as its effect has a strong influence on the hardening curve at small strain (Figure 3a). Similarly, the large strain hardening parameter has a bigger influence at large strain (Figure 3b).

Figure 3. Effect of hardening parameter on the yield stress evaluation
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).

Figure 4. Effect of temperature setting on the strain rate dependency
The strain rate $\dot{ε}$ computation depends on the value of the flag $V P$ :
- If $V P$ = 1: the plastic strain rate is used.
- If $V P$ = 2: the total strain rate is used.
- If $V P$ = 3: the deviatoric strain rate is used.
In all cases the strain rate computation includes a filtering, using the cutting off frequency $F_{c u t}$ , which is user-defined, such that:(14)
${\dot{ε}}_{f} = α {\dot{ε}}^{n} + (1 - α) {\dot{ε}}^{n - 1}$

With $α = 2 π F_{c u t} Δ 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 $σ_{0}$ , $Δ σ_{m}$ , $β$ , $Ω$ , and $η$ are ignored, and the yield stress becomes:(15)
$σ_{Y} = f_{Y} (ε_{p}, \dot{ε})$

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:(16)
$σ_{Y} = f_{Y} (ε_{p}, \dot{ε}) \frac{f_{temp} (ε_{p}, T)}{f_{temp} (ε_{p}, T_{ini})}$

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 Prager-Ziegler kinematic hardening is used.
- If 0 ≤ C_hard ≤ 1: the model is interpolated between the isotropic and the kinematic hardening.
  
  Figure 5.