/MAT/LAW101

#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/LAW101/1/1
Talc filled Polypropylene
#              RHO_I
            9.05E-3
#               EREF                  E1                  Nu                 VE1
              2700.0                                  0.2600              0.0E-1
#                VE2            EDOT_REF        GAMA_DOT_REF              ALPHAP
              1.0E-2               1.0E3           3.464E+18            8.313E-2
#            delta_H                   V                   m                  C3
            109000.0           1.325E-27               5.000            
#                 C4             ALPHAK1             ALPHAK2                  H0
                  12              1.0E-2                 0.0                  40
#            ZETA1_i                  C5                  C6                  C7
               0.000             -9.0E-3              0.6600                   0  
#                 C8                  C9                 C10                  h1
               0.200              -0.300               6.700               0.000
#            ZETA2_i                 C11                 C12                 C13
               0.000                 0.0              12.000             -7.0E-3
#                C14                  C1                  C2            LAMBDA_L
               0.600             -0.1000               8.000               4.800
#        RHO_theta_0          CV_theta_0              THETA0            ALPHA_TH
            9.05E-10               2.0E9             298.000             7.70E-5
#        THETA_GLASS         TEMP_FACTOR          THETA_FLAG              THETAi
             373.000               0.000               0.000             298.000
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#enddata

1. This material is only compatible with solid elements with Lagrange type total strain. The strain formulation flag is automatically set to I_smstr =10 in /PROP/SOLID.
2. The computational cost of this material law is high. However, this model gives very good results for polypropylene materials.
3. This model is proposed by Bouvard. ¹ It is based on a thermodynamic framework in which physically based Internal State Variables (ISVs) are selected to accurately represent the underlying physics of the polymer chain deformation. These ISVs describe the current energetic state of the polymer network and are included in the Helmholtz free energy. The model is based on three ISVs:

   ${\overline{z}}_1$

   Internal strain-like scalar induced by entanglement points. These points represent a sort of obstacle to the chain motion and, during the chain entanglement slippage at a certain amount of stress leads to a strain softening. (observed experimentally on Figure 1)

   ${\overline{z}}_2$

   Internal strain-like scalar associated to with chain alignment and coiling at large strains, that causes material hardening.

   $\beta$

   Internal strain-like tensor associated to with chain orientation and chain stretching at large strains, causing material hardening.

   
   
   Figure 1. Decomposition of the total deformation gradient

   The description of the kinematics of the problem is based on the decomposition of the deformation gradient into elastic, viscoplastic, and isotropic thermal components.(1)

   $$F=F_e{F}_p{F}_{\theta }$$
4. The model considers the incompressible plastic flow, that brings $J^p=1$ . Considering that $F_e=R_e{U}_e$ and that the elastic logarithmic strain tensor $E_e=\ln U_e$ , the mandel stress tensor $M$ can be calculated as:(2)
   $$M=2\mu \left(\theta \right)E_e+\left(K\left(\theta \right)-\frac{2}{3}\mu \left(\theta \right)\right)tr\left(E_e\right)I$$
   With,(3)

   $$\mu \left(\theta \right)=\frac{E\left(\theta \right)}{2+2\nu }$$
   And,(4)

   $$K\left(\theta \right)=2\mu \left(\theta \right)\frac{1+\nu }{3\left(1-2\nu \right)}$$
   Where,

   $E\left(\theta \right)$

   Young's modulus

   $\nu$

   Poisson coefficient

   $K\left(\theta \right)$

   Bulk moduli

   $\mu \left(\theta \right)$

   Shear moduli

5. The Cauchy stress tensor $\sigma$ can be obtained from the Mandel stress:(5)
   $$\sigma =J_e{}^{-1}{R}_{e}M{R}_e{}^T$$
   The viscoplastic flow rule, which describes the evolution equation for $F^p$ is:(6)

   $${\dot{F}}_p=D_p{F}_p$$
   With,(7)

   $$D_p=\frac{1}{\sqrt{2}}{\dot{\gamma }}_p{N}_p$$
   Where,

   $N_p$

   Viscous flow direction.

   $dev$

   Deviatoric part of the tensor.

   (8)

   $$N_p=\frac{dev\left(M-\alpha \right)}{\Vert dev\left(M-\alpha \right)\Vert }$$
   The $\alpha$ tensor is the corresponding stress-like tensor to the strain-like tensor $\beta$ , representing the chain stretching at large deformation. The viscous flow is given by:(9)

   $${\dot{\gamma }}_p={\dot{\gamma }}_{p0}\times e^{-\frac{\text{Δ}H}{R\theta }}{\left[\sinh \left(\frac{{\tau }_{eq}V}{k_B\theta }\right)\right]}^m$$
   Where,

   ${\dot{\gamma }}_p$

   Viscous shear strain rate.

   ${\tau }_{eq}$

   Amount of plastic shear stress.

   Finally, the internal state variables evolution rules are:(10)

   $${\dot{\overline{z}}}_1=h_0\left(1-\frac{{\overline{z}}_1}{z*}\right){\dot{\gamma }}^p$$

   (11)

   $${\dot{z}}^{*}=\left({\overline{z}}_{sat}^{*}-{\mathrm{g}}_0\left(\theta \right)z^{*}\right){\dot{\gamma }}^p$$

   (12)

   $${\dot{\overline{z}}}_2=h_1\left({\lambda }^p-1\right)\left(1-\frac{{\overline{z}}_2}{z_{2sat}\left(\theta \right)}\right){\dot{\gamma }}^p$$
   And(13)

   $$\dot{\beta }=R_s\left(\theta \right)\left(D_p\beta +\beta D_p\right)$$
6. Young's modulus as a function of temperature is:(14)
   $$\mathrm{E}\left(\theta \right)=\left[E_{ref}+{\mathrm{E}}_1\left(\theta -{\theta }_{ref}\right)\right]\left[1+\frac{V{E}_1}{1+\exp \left(-\frac{\log \dot{\epsilon }-\log {\dot{\epsilon }}_{ref}}{V{E}_2}\right)}\right]$$
7. Internal State Variables (ISVs):

   ISV 1	ISV 2	ISV 3
   $k_1={\alpha }_{k1}\mu \left(\theta \right){\overline{z}}_1$	$k_2={\alpha }_{k2}\mu \left(\theta \right){\overline{z}}_2$	$\alpha ={\mu }_B\left(\theta \right)\beta$
   ${\dot{\overline{z}}}_1=h_0\left(1-\frac{{\overline{z}}_1}{z*}\right){\dot{\gamma }}^p$	${\dot{\overline{z}}}_2=h_1\left({\lambda }^p-1\right)\left(1-\frac{{\overline{z}}_2}{z_{2sat}\left(\theta \right)}\right){\dot{\gamma }}^p$	$\dot{\beta }=R_s\left(\theta \right)\left(D_p\beta +\beta D_p\right)$
   ${\dot{z}}^{*}=\left({\overline{z}}_{sat}^{*}-{\mathrm{g}}_0\left(\theta \right)z^{*}\right){\dot{\gamma }}^p$	${\overline{z}}_{2sat}\left(\theta \right)=C_{11}\left(\theta -{\theta }_{ref}\right)+C_{12}$	${\mu }_B\left(\theta \right)={\mu }_R\left(\theta \right){\left(1-\frac{tr\left(\beta \right)-3}{{\lambda }_L}\right)}^{-1}$
   ${\overline{z}}_{sat}^{*}\left(\theta \right)=C_7\left(\theta -{\theta }_{ref}\right)+C_8$		${\mu }_R\left(\theta \right)=C_1\left(\theta -{\theta }_{ref}\right)+C_2$
   $g_0\left(\theta \right)=C_9\left(\theta -{\theta }_{ref}\right)+C_{10}$		$R_s\left(\theta \right)=C_{13}\left(\theta -{\theta }_{ref}\right)+C_{14}$
   $z_0^{*}\left(\theta \right)=C_5\left(\theta -{\theta }_{ref}\right)+C_6$

8. Flow rule:
   • ${\dot{F}}_p=D_p{F}_p$
   • $D_p=\frac{1}{\sqrt{2}}{\dot{\gamma }}_p{N}_p$
   • $N^p=\frac{dev\left(M-\alpha \right)}{\Vert dev\left(M-\alpha \right)\Vert }$
   • ${\dot{\gamma }}_p={\dot{\gamma }}_{p0}{e}^{-\frac{\Delta H}{k\theta }}{\left[\sinh \left(\frac{{\tau }_{eq}V}{k_B\theta }\right)\right]}^m$
   • ${\tau }_{eq}=\tau -\left(\mathrm{Y}\left(\theta \right)+{\kappa }_1+{\kappa }_2+{\alpha }_p\overline{\pi }\right)$
   • $\mathrm{Y}\left(\theta \right)=C_3\left(\theta -{\theta }_{ref}\right)+C_4$
   • $\tau =\frac{1}{\sqrt{2}}\Vert dev\left(M-\alpha \right)\Vert$
9. Heat generation (in adiabatic conditions):
   • $c_v\dot{\theta }=\omega M:D_p$
   • $\rho \left(\theta \right)=\rho \left({\theta }_{ref}\right)\frac{1.42{\theta }_g+44.7}{1.42{\theta }_g+0.15\theta }$
   • $c_v=c_v\left({\theta }_{ref}\right)\left(0.106+3\times {10}^{-3}\theta \right)$

   Where,

   Notation	Description	Notation	Description
   ${\alpha }_{k1}$	Material parameter	$D_p$	Plastic component of the velocity gradient
   ${\alpha }_{k2}$	Material parameter	$E$	Logarithmic strain tensor
   $\alpha$	Internal shear stress	$E_{ref}$	Young's modulus at reference temperature
   ${\alpha }_p$	Pressure sensitivity parameter	$E_1$	Material parameter
   $\beta$	ISV 3 tensor	$F$	Deformation gradient
   ${\dot{\gamma }}_p$	Viscous flow	$F_e$	Elastic component of $F$
   ${\epsilon }^{nom}$	Nominal strain	$F_p$	Plastic component of $F$
   ${\epsilon }^{true}$	Hencky (true) strain	$F_{\theta }$	Thermal component of $F$
   $\eta$	Viscosity	$g_0$	Material parameter
   $\theta$	Temperature	$\text{Δ}H$	Activation energy
   ${\kappa }_1\text{, }{\kappa }_2$	Internal stress fields induced by entanglement points	$h_0$	Hardening modulus
   $\lambda$	Stretch	$h_1$	Hardening modulus
   ${\lambda }_L$	Network locking stretch	$J$	Determinant of
   ${\lambda }_p$	Equivalent plastic stretch	$K$	Elastic bulk modulus
   $\mu$	Elastic shear modulus	$k_B$	Boltzmann constant
   ${\mu }_B$	Internal shear stress modulus	$M$	Mandel stress
   ${\mu }_R$	Rubbery modulus	$m$	Viscous flow exponent
   $\nu$	Poisson coefficient	$N_p$	Viscous flow direction
   $\overline{\pi }$	Effective pressure	$R$	Gas constant
   $\sigma$	Cauchy (true) stress	$R_s$	Material parameter
   $\tau$	Equivalent shear stress	$R\text{, }U$	Antisymmetric, symmetric part of $F$
   $\psi$	Free energy	$V$	Activation volume
   $\omega$	Conversion factor	$Y$	Yield surface
   $C$	Cauchy-Green tensor	${\overline{z}}_1$	ISV 1
   $c_v$	Heat capacity	${\overline{z}}_2$	ISV 2
   $C_1,C_2,\mathrm{...},C_{14}$	Material parameters	${\overline{z}}^{*}$	Strain criteria for chain slip
   $D$	Velocity gradient	${\overline{z}}_{sat}^{*}$	Saturation value for ${\overline{z}}^{*}$