/MAT/LAW87 (BARLAT2000)

Yield stresses can be defined either by user-defined functions (plastic strain versus stress) or analytically by a combination of Swift-Voce model. The model is based on Barlat YLD2000 criterion. ¹

Format

Definitions

Example 1 (with Barlat parameters input Iflag=0 and I_fit=0)

#RADIOSS STARTER
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/UNIT/1
unit for mat
                  kg                  mm                  ms
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#-  2. MATERIALS:
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/MAT/LAW87/1/1
Steel 
#              RHO_I
              7.8E-6                   0
#                  E                  Nu     IFlag        VP             coeff_c               exp_p              
                 210                 0.3         0         1             4.15401                3.57             
#                 a1                  a2                  a3                  a4     I_fit
                 1.0                 1.0                 1.0                 1.0         0
#                 a5                  a6                  a7                  a8
                 1.0                 1.0                 1.0                 1.0
#              Chard
                   0
#              exp_a               ALPHA                NEXP                Fcut   Fsmooth     NRATE
                   2                   0                   0                   0         1         1
# Blank

#  func_id                        YSCALE         strain rate
         4                           1.5                   1
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/FUNCT/4
Steel
#                  X                   Y
                   0                  .3
               0.007                  .5
                0.05                  .7
                 0.1                 .75
                 0.3                  .9
                   1                 1.2				 
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#ENDDATA
/END
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

Example 2 (with experiment data input I_fit=1)

#RADIOSS STARTER
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/UNIT/1
unit for mat
                   g                  mm                  ms
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/MAT/LAW87/1/1
Aluminum
#              RHO_I
              2.7E-3                   0
#                  E                  Nu     IFlag        VP             coeff_c               exp_p 
               70000                 0.3         1         0                   0                   0
#              sig00               sig45               sig90                sigb     I_fit
          133.179899          133.102756          132.330693          162.330301         1
#                r00                 r45                 r90                  rb
         0.703242569         0.486264221         0.865336191         0.546807587
#              Chard
                   0
#              exp_a               ALPHA                NEXP                Fcut   Fsmooth
                   8                0.55                0.21                   0         1
#             ASwift                Eps0               Qvoce                Beta                  KO
                415.             0.00220               174.7               11.19               132.4
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#ENDDATA
/END
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

Example 3 (with Hansel yield model (Iflag=2) and kinematic hardening model (Chard=1))

#RADIOSS STARTER
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/UNIT/1
unit for mat
                  kg                  mm                  ms
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/MAT/BARLAT2000/2/1
Steel
#              RHO_I
            7.800E-6                   0
#                  E                  Nu     IFlag        VP                   c                   P
                 210                  .3         2         0                   0                   0
#                 a1                  a2                  a3                  a4     I_fit
              0.4865              1.3783              0.7536              1.0246         0
#                 a5                  a6                  a7                  a8
              1.0363              0.9036              1.2321              1.4858
#              Chard  
                   1
#              exp_a
                   8                                                                                
#                 AM                  BM                  CM                  DM                  PM
               0.578               0.185               -6.78                0.02                7.54
#                 QM              E0MART                 VM0
              1379.0                0.01              0.1690
#                AHS                 BHS                 MHS                 NHS              EPS0HS
              -0.261               9.170               0.118               0.401              0.0988
#              HMART                  K1                  K2
              0.5490                3.95            -0.00681
#              TEMP0                TREF                  CP                 ETA
                300.                293.                460.                 0.1
#               CRC1                CRA1                CRC2                CRA2
                  80               0.052                   0                  0. 
#               CRC3                CRA3                CRC4                CRA4
                   0                 0.0                   0                  0.  
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#ENDDATA
/END
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

Comments

1. The yield function is expressed as:(1)
   $$f=\overline{\sigma }-{\sigma }_y$$
   (2)
   $$\overline{\sigma }= \frac{1}{2^{\raisebox{1ex}{$1$}\!\left/ \!\raisebox{-1ex}{$a$}\right.}}\left({\phi }^{\prime }\left(X^{\prime }\right)+{\phi }^{″}\left(X^{″}\right)\right){}^{\raisebox{1ex}{$1$}\!\left/ \!\raisebox{-1ex}{$a$}\right.}$$
   (3)
   $${\phi }^{\prime }\left(X^{\prime }\right)={\left|{X^{\prime }}_1-{X^{\prime }}_2\right|}^a$$
   (4)
   $${\phi }^{″}\left(X^{″}\right)={\left|2{X^{″}}_2+{X^{″}}_1\right|}^a+{\left|2{X^{″}}_1+{X^{″}}_2\right|}^a$$
   $X\text{'}$ and $X"$ denote the principal values of the tensors $X\text{'}$ and $X"$ which are a linear transformation of the stress deviator, which leads to:(5)

   $${\phi }^{\prime }\left(X^{\prime }\right)= {\left[{\left({X^{\prime }}_{xx}-{X^{\prime }}_{yy}\right)}^2+4{\left({X^{\prime }}_{xy}\right)}^2\right]}^{\frac{a}{2}}$$

   (6)

   $$\begin{array}{l} {\phi }^{″}\left(X^{″}\right)={\left[\frac{3}{2}\left({X^{″}}_{xx}-{X^{″}}_{yy}\right)+\frac{1}{2}\sqrt{{\left({X^{″}}_{xx}-{X^{″}}_{yy}\right)}^2+4{\left({X^{″}}_{xy}\right)}^2}\right]}^a+\\ {\left[\frac{3}{2}\left({X^{″}}_{xx}-{X^{″}}_{yy}\right)-\frac{1}{2}\sqrt{{\left({X^{″}}_{xx}-{X^{″}}_{yy}\right)}^2+4{\left({X^{″}}_{xy}\right)}^2}\right]}^a\end{array}$$

   The tensors $X\text{'}$ and $X"$ are linear transformations of the stress tensor:

   $X^{\prime }=L^{\prime }\sigma and X^{″}=L^{″}\sigma$ (7)

   $$L^{\prime }=\frac{1}{3}\left[\begin{array}{ccc}2{\alpha }_1& -{\alpha }_1& 0\\ -{\alpha }_2& 2{\alpha }_2& 0\\ 0& 0& 3{\alpha }_7\end{array}\right]$$

   (8)

   $$L^{″}=\frac{1}{9}\left[\begin{array}{ccc}-2{\alpha }_3+2{\alpha }_4+8{\alpha }_5-2{\alpha }_6& {\alpha }_3-4{\alpha }_4-4{\alpha }_5+4{\alpha }_6& 0\\ 4{\alpha }_3-4{\alpha }_4-4{\alpha }_5+{\alpha }_6& -2{\alpha }_3+8{\alpha }_4+2{\alpha }_5-2{\alpha }_6& 0\\ 0& 0& 9{\alpha }_8\end{array}\right]$$
2. The yield stress could be defined either by tabulated input or using the analytic Swift-Voce model.
   • Iflag=0: Tabulated.
     • It is possible to add total strain rate dependency by defining a number N_rate of functions.
   • Iflag=1: The analytic Swift-Voce model is expressed as:(9)
     $${\sigma }_y={\alpha }_{sv}\left[A{\left({\overline{\epsilon }}_p+{\epsilon }_0\right)}^n\right]+\left(1-{\alpha }_{sv}\right)\left[K_0+Q\left(1-\exp \left(-B{\overline{\epsilon }}_p\right)\right)\right]$$
     Where,

     ${\overline{\epsilon }}_p$

     Equivalent plastic strain.

   • Iflag=2: Hansel hardening model is considered.(10)
     $${\sigma }_y=\left\{B_{HS}-\left(B_{HS}-A_{HS}\right)e^{\left(-m{\left[{\overline{\epsilon }}^p+{\epsilon }_0\right]}^n\right)}\right\}\left(K_1-K_{2}T\right)+\text{Δ}{H}_{\gamma \alpha }{V}_m$$
     Temperature is updated in the law when adiabatic conditions:(11)

     $$\text{Δ}T=\eta \frac{\overline{\sigma }d{\overline{\epsilon }}_p}{\rho C_p}$$
     The martensite rate equation is computed as follows:(12)

     $$\frac{\partial V_m}{\partial \epsilon }=\left\{\begin{array}{cc}0& if\epsilon <{\epsilon }_0\\ \frac{B}{A}\cdot e^{\left(\frac{Q}{T}\right)}\cdot {\left(\frac{1-V_m}{V_m}\right)}^{\left(\frac{B+1}{B}\right)}\cdot \frac{{\left(V_m\right)}^p}{2}\cdot \left(1-\tanh \left[C+DT\right]\right)& f\epsilon \ge {\epsilon }_0\end{array}\right.$$
3. If Chard>0, a kinematic hardening model of Chaboche Rousselier is used:
   • The back stress is calculated as follows:(13)
     $$a={\displaystyle \sum _{i=1}^4{a}_i}$$
     With,(14)
     $$a_i=A_i{C}_{i}d{\epsilon }_p-C_i{a}_i\text{Δ}{\overline{\epsilon }}_p$$
   • The yield stress is computed as follows when combined isotropic kinematic hardening is chosen:(15)
     $${\sigma }_y=(1-Chard).{\sigma }_{iso\_hard}+Chard.{\sigma }_{kin\_hard}$$
4. The strain rate filtering is available to smooth strain rates when tabulated input is chosen.
   List of Animation output (in /ANIM/SHELL/USRII/JJ):

   • USR 1= plastic strain
   • USR 2= effective stress
   • USR 3= increment of plastic strain
5. When Iflag=1 (analytic Swift-Voce formulation is used) strain rates effect is taken into account using Cowper-Symonds expression:(16)
   $${\sigma }_y={\sigma }_y\left(1+{\left(\frac{\dot{\epsilon }}{c}\right)}^{\raisebox{1ex}{$1$}\!\left/ \!\raisebox{-1ex}{$p$}\right.}\right)$$

   If VP=0: $\dot{\epsilon }$ is the total strain rate.

   If VP=1: $\dot{\epsilon }$ is the plastic strain rate.

   If c=0 or p=0, the strain rate effects are not taken into account.

6. When Iflag=0 (tabulated formulation) then:

   If VP=0: ${\dot{\epsilon }}_i$ is the total strain rate.

   If VP =1: ${\dot{\epsilon }}_i$ is the plastic strain rate.

7. Strain rate filtering:

   If VP=0 (dependency on strain rate), the default value of F_cut = 10KHz.

   If VP=1 (dependency on plastic strain rate), F_smooth and F_cut are ignored.

8. If I_fit=1, the coefficients ${\alpha }_i$ will be automatically fit in the Radioss Starter. The tensile yield strengths ${\sigma }_{00},{\sigma }_{45},{\sigma }_{90}$ and Lankford ratios $r_{00},r_{45},r_{90}$ must be determined from uniaxial tension experiments along the rolling, diagonal and transverse directions at an amount of plastic work corresponding to a plastic strain equal to 0.2%. ${\sigma }_b$ and $r_b$ should be determined from biaxial test, for the same amount of plastic strain.