/MAT/LAW112 (PAPER or XIA)

#RADIOSS STARTER
/UNIT/1
unit for mat
                  Mg                  mm                   s
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/MAT/LAW112/1/1
Xia
#              RHO_I
            7.83E-10
#                 E1                  E2                  E3      Ires      Itab   Ismooth       
                4193                1554                1554         2         0         0
#               nu21                 G12                 G23                 G13  
              0.1011                 988                  76                  76
#                  K                 E3C                  CC                 
                 2.0                47.2               24.46                 
#               nu1p                nu2p                nu4p                nu5p
               0.555              0.1537                0.18               0.145
#                S01                 A01                 B01                 C01
                12.0                19.0               260.0               800.0
#                S02                 A02                 B02                 C02
                 6.5                40.0               160.0               250.0
#                S03                 A03                 B03                 C03       
                 6.0                11.0               100.0               125.0
#                S04                 A04                 B04                 C04
                 7.3                 6.0               160.0               300.0
#                S05                 A05                 B05                 C05
                 6.3                 9.0               310.0               225.0
#               ASIG                BSIG                CSIG
               16.55               16.55                3.16
#               TAU0                ATAU                BTAU
                 2.1                 9.0                 2.0    
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#ENDDATA

#RADIOSS STARTER
/UNIT/1
unit for mat
                  Mg                  mm                   s
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|          
/MAT/LAW112/1/1
Xia_tab
#              RHO_I
            7.83E-10
#                 E1                  E2                  E3      Ires      Itab   Ismooth       
                4193                1554                1554         1         1         1
#               nu21                 G12                 G23                 G13 
              0.1011                 988                  76                  76
#                  K                 E3C                  CC                 
                 2.0                47.2               24.46                 
#               nu1p                nu2p                nu4p                nu5p
               0.555              0.1537                0.18               0.145
#           TAB_YLD1         MAT_Xscale1         MAT_Yscale1
                  25                 1.0                 1.0
#           TAB_YLD2         MAT_Xscale2         MAT_Yscale2
                  25                 1.0                0.35
#           TAB_YLD3         MAT_Xscale3         MAT_Yscale3
                  25                 1.0                0.75
#           TAB_YLD4         MAT_Xscale4         MAT_Yscale4
                  25                 1.0              0.6341
#           TAB_YLD5         MAT_Xscale5         MAT_Yscale5
                  25                 1.0                 0.5
#           TAB_YLDC         MAT_XscaleC         MAT_YscaleC
                  25                 1.0                 0.5
#           TAB_YLDS         MAT_XscaleS         MAT_YscaleS
                  25                 1.0                 0.5 
/FUNCT/46
ecoulement2   
#     plastic strain              stress                                                                                
                 0.0	           12.00
               0.012	 32.979020979021
               0.025	50.4615384615385
                0.05	            74.5
               0.075	90.9473684210526
                 0.1	102.909090909091
               0.125	          112.00
                0.15	119.142857142857
               0.175	124.903225806452
                 0.2	129.647058823529
                0.25	          137.00
                 0.3	142.434782608696
                 0.4	149.931034482759
                 0.5	154.857142857143
                 1.0	165.846153846154   
/TABLE/1/25
Yld Functions : plastic strain + strain rate dependency
#DIMENSION
         2
#   FCT_ID                   strain rate                                                     Scale_y
        46                           0.0                                                        1.00
        46                           1.0                                                        1.10
        46                           5.0                                                        1.15
        46                          10.0                                                        1.20
        46                         100.0                                                        1.25
        46                      100000.0                                                        1.35  
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#enddata

1. To describe the behavior of the paperboard material law, the following orthotropic direction is considered.

   
   
   Figure 1.

2. The elastic behavior of this material law is orthotropic.
   The in-plane behavior should be fully uncoupled with the out-of-plane behavior, computed as:(1)

   $$\{\begin{array}{l}{\sigma }_{xx}=C_{11}{\epsilon }_{xx}+C_{12}{\epsilon }_{yy}\\ {\sigma }_{yy}=C_{21}{\epsilon }_{xx}+C_{22}{\epsilon }_{yy}\\ {\sigma }_{xy}=G_{12}{\gamma }_{xy}\end{array}$$

   With $C=\frac{1}{1-{\nu }_{12}{\nu }_{21}}\left[\begin{array}{cc}{E}_1& {\nu }_{12}{E}_2\\ {\nu }_{21}{E}_1& E_2\end{array}\right]$

   The transverse shear components are computed as:(2)

   $$\{\begin{array}{l}{\sigma }_{yz}=G_{23}{\epsilon }_{yz}\\ {\sigma }_{zx}=G_{21}{\epsilon }_{zx}\end{array}$$
   The out-of-plane elastic behavior (for solids only) is treated as a uniaxial equivalent problem. However, the computation of the stress may differ between tension and compression. The elasticity becomes nonlinear for compressive loadings.(3)

   $$\begin{array}{ccc}{\sigma }_{zz}=E_3{\epsilon }_{zz}^e& \text{if}& {\epsilon }_{zz}^e\ge 0\\ {\sigma }_{zz}=E_{3C}(1-{\mathrm{e}}^{-C_c{\epsilon }_{zz}^e})& \text{if}& {\epsilon }_{zz}^e<0\end{array}$$
3. In the Xia 2002 formulation, the in-plane yield criterion, denoted as $f$ , is defined as:(4)
   $$f={{\displaystyle \sum _{I=1}^6{\chi }_I\left(\frac{\sigma :N_I}{{\sigma }_Y^I}\right)}}^{2k}-1$$

   Where,

   ${\chi }_I=\{\begin{array}{ccc}1& \text{if}& \sigma :N_I>0\\ 0& \text{otherwise}& \end{array}$

   ${\chi }_I$

   Switching parameters.

   $\sigma$

   Cauchy stress tensor.

   $N_I$

   Normal direction of the yield planes.

   ${\sigma }_Y^I$

   Yield stresses.

   $k$

   Positive integer.

   Each direction is associated to a given loading direction following the order defined below:

   1

   Tension in orthotropic direction 1.

   2

   Tension in orthotropic direction 2.

   3

   Positive in-plane shear.

   4

   Compression in orthotropic direction 1.

   5

   Compression in orthotropic direction 2.

   6

   Negative in-plane shear (same input as positive in-plane shear ${\sigma }_Y^6={\sigma }_Y^3$ ).

   The normal direction vector to the yield planes are:

   $\begin{array}{l}{N}_1=\left[\begin{array}{cccccc}\frac{1}{\sqrt{1+{\nu }_{1p}^2}}& -\frac{{\nu }_{1p}}{\sqrt{1+{\nu }_{1p}^2}}& 0& 0& 0& 0\end{array}\right]\\ N_2=\left[\begin{array}{cccccc}-\frac{{\nu }_{2p}}{\sqrt{1+{\nu }_{2p}^2}}& \frac{1}{\sqrt{1+{\nu }_{2p}^2}}& 0& 0& 0& 0\end{array}\right]\\ N_3=\left[\begin{array}{cccccc}0& 0& 0& 1& 0& 0\end{array}\right]\\ N_4=\left[\begin{array}{cccccc}-\frac{1}{\sqrt{1+{\nu }_{4p}^2}}& \frac{{\nu }_{4p}}{\sqrt{1+{\nu }_{4p}^2}}& 0& 0& 0& 0\end{array}\right]\\ N_5=\left[\begin{array}{cccccc}\frac{{\nu }_{5p}}{\sqrt{1+{\nu }_{5p}^2}}& -\frac{1}{\sqrt{1+{\nu }_{5p}^2}}& 0& 0& 0& 0\end{array}\right]\\ N_6=\left[\begin{array}{cccccc}0& 0& 0& -1& 0& 0\end{array}\right]\end{array}$

   Each direction $I$ is then associated to a specific yield stress whose expression is:

   $\begin{array}{cc}{\sigma }_Y^I=S_I^0+A_I^0\tanh \left(B_I^0{\epsilon }_p^f\right)+C_I^0{\epsilon }_p^f& \begin{array}{cc}\text{with}& I\in \left[1,6\right]\end{array}\end{array}$

   Where, ${\epsilon }_p^f$ is the in-plane equivalent plastic strain (associated to the yield function $f$ ).

   The out-of-plane yield function is denoted as $g$ is defined as:

   $g=-{\sigma }_{zz}-{\sigma }_Y^C$ with ${\sigma }_Y^C=A_{\sigma }+B_{\sigma }\exp (C_{\sigma }{\epsilon }_p^g)$

   Where, ${\epsilon }_p^g$ is the out-of-plane equivalent plastic strain (associated to the yield function $g$ ).

   The transverse shear yield function is:(5)

   $$h=\frac{\sqrt{{\sigma }_{yz}^2+{\sigma }_{zx}^2}}{{\sigma }_Y^S}-1$$
   Where,

   ${\sigma }_Y^S={\tau }_0+\left[A_{\tau }-\min (0,{\sigma }_{zz})B_{\tau }\right]{\epsilon }_p^h$

   ${\epsilon }_p^h$

   Out-of-plane equivalent plastic strain (associated to the yield function $h$ ).

   If the tabulated yield stress option is selected (Itab = 1), each yield stress is associated to a table (TAB_YLDi) to define the stress evolution with the plastic strain, at several plastic strain-rate. Two scale factors can be also defined in the X and Y direction for each table. In this case, the hardening parameters $S0i$ , $A0i$ , $B0i$ , $C0i$ , $A_{\sigma }$ , $B_{\sigma }$ , $C_{\sigma }$ , ${\tau }_0$ , $A_{\tau }$ , and $B_{\tau }$ are ignored, and the yield stress becomes:

   $\begin{array}{l}\begin{array}{cc}{\sigma }_Y^I=f_Y^{tab\_YLDI}({\epsilon }_p^f,{\dot{\epsilon }}_p^f)& I\in \left[1,6\right]\end{array}\\ {\sigma }_Y^C=f_Y^{tab\_YLDC}({\epsilon }_p^g,{\dot{\epsilon }}_p^g)\\ {\sigma }_Y^S=f_Y^{tab\_YLDS}({\epsilon }_p^h,{\dot{\epsilon }}_p^h)\end{array}$

   For output field, an equivalent “global” plastic strain is computed as:(6)

   $${\epsilon }_p=\sqrt{{({\epsilon }_p^f)}^2+{({\epsilon }_p^g)}^2+{({\epsilon }_p^h)}^2}$$