/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. この破壊モデルは、/MAT/LAW92、/MAT/LAW95、/MAT/LAW100の各材料でのみ使用できます。
2. 最初の載荷プロセス時の応力は、無損傷応力に等しい値にします。除荷時および再載荷時には、次のように正の軟化係数関数を乗算した応力とします。(1)
   $$\sigma =\eta dev\left(\sigma \right)-p\text{Ι}$$
   ここで、

   $dev\left(\sigma \right)$

   応力の偏差部分。

   $p$

   静水圧。

   $\eta$

   以下を考慮した超弾性モデルのひずみエネルギーの関数である損傷係数

   (2)

   $$\eta =\{\begin{array}{c}=1\text{, } \text{if} W=W_{max}\\ <1\text{, if} W<W_{max}\end{array}$$
   ここで、 $W_{max}$ は、載荷の過程でこの材料に作用した最大ひずみエネルギーです。(3)

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

   ここで、 $erf$ はガウスの誤差関数です。

3. パラメータ $R$ の値が大きいほど、1からの $\eta$ の偏差が小さくなり、損傷が発生する可能性も低くなります。mの値が小さいと、ひずみが少なくても損傷が大きくなります。mの値が大きいと、初期載荷時ではひずみが軽微であれば損傷が小さくすみますが、再載荷時ではひずみが少なくても損傷が大きくなります。 $\beta$ の値が小さいほど、損傷が大きくなります。
4. この破壊モデルでは、破壊や要素の除去は発生しません。