/MAT/LAW93 (ORTH_HILL) or (CONVERSE)

It could be used with property set /PROP/TYPE11, /PROP/TYPE17, /PROP/TYPE51, /PROP/PCOMPP for shell and /PROP/TYPE6 for solid.

Format

Definitions

Example (Curve Input)

#RADIOSS STARTER
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/UNIT/1
unit for mat
                  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/LAW93/1/1
plastic
#              RHO_I
           2.730E-09
#                E11                 E22                 E33                 G12                Nu12
              225654              195400              178526            75187.97                0.30
#                G13                 G23                Nu13                Nu23
            75187.97            75187.97                0.28                0.32
#     Nrate       VP                Fcut
         2         1                 0.0
#   Ifunct                        Yscale              Epsdot
         5                           1.0                0.01	     
         5                           1.5               100.0	       
#               SigY                 QR1                 CR1                 QR2                 CR2
                   0                   0                   0                 0.0                 0.0
#                R11                 R22                 R12
                 1.0             1.05626             0.96425
#                R33                 R13                 R23
              0.9337                 1.0                 1.0
/FUNCT/5
plastic
#                  X                   Y           
                   0         165.6362749
               0.002         173.8123558
               0.005         180.2967164
                0.01         186.5926709
                0.02         193.8182168
                0.05         204.4407991
                0.07         208.5903797
                 0.1         213.1182051
                0.12         215.4817557
                0.15         218.4183864
                0.17         220.0863912
                 0.2         222.2743041
                0.22         223.5689486
                0.25         225.3186882
                0.27         226.3794409
                 0.3         227.840544
                0.32         228.7406278
                0.35         229.996802
                0.37         230.7795124
                 0.4         231.8824363
                 0.5         235.0704031
                 0.6         237.7095003
                 0.7         239.9650034
                 0.8         241.9367878
                 0.9         243.689935
                   1         245.2692715
                 1.5         251.4456403
                   2         255.9237789
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#ENDDATA
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

Comments

1. The yield stress is compared to an equivalent stress in the orthotropic frame. For solid elements, this equivalent stress is defined as:(1)
   $${\sigma }_{eq}=\sqrt{F\left({\sigma }_{22}^2-{\sigma }_{33}^2\right)+G\left({\sigma }_{33}^2-{\sigma }_{11}^2\right)+H\left({\sigma }_{11}^2-{\sigma }_{22}^2\right)+2L{\sigma }_{23}^2+2M{\sigma }_{31}^2+2N{\sigma }_{12}^2}$$
   Where,

   $F=\frac{1}{2}(\frac{1}{R_{22}^2}+\frac{1}{R_{33}^2}-\frac{1}{R_{11}^2})$

   $G=\frac{1}{2}(\frac{1}{R_{33}^2}+\frac{1}{R_{11}^2}-\frac{1}{R_{22}^2})$

   $H=\frac{1}{2}(\frac{1}{R_{22}^2}+\frac{1}{R_{11}^2}-\frac{1}{R_{33}^2})$

   $L=\frac{3}{2R_{23}^2}$

   $M=\frac{3}{2R_{31}^2}$

   $N=\frac{3}{2R_{12}^2}$

   $R_{ii}^{}=\frac{{\sigma }_F^{ii}}{{\sigma }_F}$

   Normal directions

   $R_{ij}^{}=\frac{\sqrt{3}{\sigma }_F^{ij}}{{\sigma }_F}$

   Shear directions

   ${\sigma }_F^{ij}$

   Yield stress in direction $ij$

   ${\sigma }_F^{}$

   Global flow stress that can either be defined with a sum of Voce hardening, or can be tabulated (see below).

   Under plane-stress conditions, for shell elements, the equivalent yield stress becomes:(2)

   $${\sigma }_{eq}=\sqrt{\left(G+H\right){\sigma }_{11}^2+\left(F+H\right){\sigma }_{22}^2-2H{\sigma }_{11}{\sigma }_{22}+2N{\sigma }_{12}^2}$$

2. The yield function $\Phi$ will compare the Hill’s equivalent stress ${\sigma }_{eq}$ to the flow stress ${\sigma }_F^{}$ as:(3)
   $$\Phi ={\sigma }_{eq}-{\sigma }_F$$
   The two different ways to define the flow stress ${\sigma }_F^{}$ are: parameter input or curve input

   • For parameter input, the flow stress is defined with an initial yield stress and a double Voce hardening as:(4)
     $${\sigma }_F={\sigma }_Y^0+R({\epsilon }_p)$$

     With $R({\epsilon }_p)={\displaystyle \sum _i^{2}Q{R}_i\cdot \left(1-e^{-C{R}_i\cdot {\epsilon }_p}\right)}$ .

   • For curve input, the parameters input values will be ignored.
     The yield can be defined with using stress versus plastic strain curve taking in account the strain rate effect. When the stress versus strain curves are defined, this is the default method for defining the hardening.

     1. If $\dot{\epsilon }\le {\dot{\epsilon }}_n$ , the yield is interpolated between ${\mathrm{f}}_n$ and ${\mathrm{f}}_{n-1}$ .
     2. If $\dot{\epsilon }\le {\dot{\epsilon }}_1$ function, ${\mathrm{f}}_1$ is used.
     3. Above ${\dot{\epsilon }}_{\max }$ , yield is extrapolated.

        
        
        Figure 1.

3. For tabulated flow stress, 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 total strain rate is used

   In all cases the strain-rate computation includes a filtering. The cutoff frequency is automatically set for VP = 1. However, for VP = 1 or 3, you can input a cutoff frequency F_cut; otherwise, a default value will be set.