/MAT/LAW57 (BARLAT3)

This is an elasto-plastic orthotropic law for modeling anisotropic materials in forming processes especially aluminum alloys. This material law must be used with property set type /PROP/TYPE9 (SH_ORTH) or /PROP/TYPE10 (SH_COMP).

Format

Definitions

Example (Steel)

#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----|
#-  2. MATERIALS:
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/MAT/LAW57/1/1
Steel
#              RHO_I
                .008
#                  E                  NU
              206000          .300000012
#  fct_IDE                         E_INF                  CE
         0                             0                   0
#                r00                 r45                 r90              C_hard                   m
                1.79                1.51                2.27                   0                   0
#           EPSP_max               EPS_T               EPS_M
                   0                   0                   0
#   fct_ID                      Fscale_i               EPS_i
         5                             0                   0
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#-  3. FUNCTIONS:
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/FUNCT/5
function_5
#                  X                   Y
                   0                 157
                  .1                 320
                  .5                 480
                 1.2                 600
#---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 anisotopic yield criteria F for plane stress is defined by:(1)
   $$F=a{\left|K_1+K_2\right|}^m+a{\left|K_1-K_2\right|}^m+c{\left|2K_2\right|}^m-2{{\sigma }_y}^m=0$$

   Where,

   ${\sigma }_y$ is the yield stress

   $K_1=\frac{{\sigma }_{xx}+h{\sigma }_{yy}}{2}$ and $K_2=\sqrt{{\left(\frac{{\sigma }_{xx}-h{\sigma }_{yy}}{2}\right)}^2+p^2{\sigma }_{xy}^2}$

2. Angles for Lankford parameters are defined with respect to orthotropic direction 1. The material constants a, c, h, and p are obtained from the three Lankford parameters:(2)
   $$\begin{array}{l}a=2-2\sqrt{\frac{r_{00}}{1+r_{00}}}\sqrt{\frac{r_{90}}{1+r_{90}}}\\ c=2-a\\ h=\sqrt{\frac{r_{00}}{1+r_{00}}}\sqrt{\frac{1+r_{90}}{r_{90}}}\end{array}$$
   Material constant p is calculated by solving:(3)

   $$\frac{2m\cdot {\sigma }_y{}^m}{\left(\frac{\partial F}{\partial {\sigma }_{xx}}+\frac{\partial F}{\partial {\sigma }_{yy}}\right){\sigma }_{45}}-1-r_{45}=0$$
3. If the last point of the first (static) function equals 0 in stress, the default value of ${\epsilon }_p^{max}$ is set to the corresponding value of ${\epsilon }_p$ .
4. If ${\epsilon }_p$ (plastic strain) reaches ${\epsilon }_p^{\max }$ , in one integration point, the corresponding shell element is deleted.
5. If the largest principal strain ${\epsilon }_1>{\epsilon }_t$ , the stress is reduced using the following relation:(4)
   $$\sigma =\sigma \left(\frac{{\epsilon }_m-{\epsilon }_1}{{\epsilon }_m-{\epsilon }_t}\right)$$
6. If ${\epsilon }_1>{\epsilon }_m$ , the stress is reduced to 0 (but the element is not deleted).
7. The maximum number of curves is 10.
8. If $\dot{\epsilon }\le {\dot{\epsilon }}_n$ , the yield is interpolated between f_n and f_n-1.
9. If $\dot{\epsilon }\le {\dot{\epsilon }}_1$ , function f₁ is used.
10. Above ${\dot{\epsilon }}_{\max }$ , yield is extrapolated.

    
    
    Figure 1.

11. The evolution of Young's modulus:
    • If fct_ID_E > 0, the curve defines a scale factor for Young's modulus evolution with equivalent plastic strain, which means the Young's modulus is scaled by the function $\mathrm{f}\left({\overline{\epsilon }}_p\right)$ :(5)
      $$E\left(t\right)=E\cdot \mathrm{f}\left({\overline{\epsilon }}_p\right)$$

      The initial value of the scale factor should be equal to 1 and it decreases.

    • If fct_ID_E = 0, the Young's modulus is calculated as:(6)
      $$E\left(t\right)=E-\left(E-E_{inf}\right)\left[1-\exp \left(-C_E{\overline{\epsilon }}_p\right)\right]$$

      Where, E and E_inf are respectively the initial and asymptotic value of Young's modulus, and ${\overline{\epsilon }}_p$ is the accumulated equivalent plastic strain.

      Note: If fct_ID_E = 0 and C_E = 0, Young's modulus, E is kept constant.