/FAIL/MULLINS_OR

#RADIOSS STARTER
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/UNIT/1
unit for material and failure
                  Mg                  mm                   s
/MAT/LAW100/1/1 
Neo Hookean material
#
 1.000000000000E-09
#N_NETWORK   FLAG_HE   FLAG_CR
         0         3          
#                C10                  D1 
              0.5000
/FAIL/MULLINS_OR/1/1
#                  R                BETA                   m
                 2.0                0.02                 0.2
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#ENDDATA
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

1. This failure model can only be used with materials /MAT/LAW92, /MAT/LAW95, and /MAT/LAW100.
2. The stress during the first loading process is equal to the undamaged stress. Upon unloading and reloading the stress is multiplied by a positive softening factor function:(1)
   $$\sigma =\eta dev\left(\sigma \right)-p\text{Ι}$$
   Where,

   $dev\left(\sigma \right)$

   The deviatoric part of the stress.

   $p$

   Hydrostatic pressure.

   $\eta$

   Damage factor that is a function of the strain energy of the hyperelastic model which respects:

   (2)

   $$\eta =\{\begin{array}{c}=1\text{, } \text{if} W=W_{max}\\ <1\text{, if} W<W_{max}\end{array}$$
   Where, $W_{max}$ is the maximum strain energy that the material has been subjected to during its loading history:(3)

   $$\eta =1-\frac{1}{R}erf\left(\frac{W_{max}-W}{m+\beta W_{max}}\right)$$

   Where, $erf$ is the Gauss error function.

3. The larger the parameter $R$ , the less $\eta$ can depart from unity and hence less damage can occur. For small values of m, there is more damage at smaller strains. For higher values of m, there is less damage in the small strain region during the initial loading, but during reloading there is will be more damage at smaller strains. Smaller values of $\beta$ result in increased damage.
4. There is no failure or element deletion with this failure model.