/FAIL/LAD_DAMA

この破壊モデルは直交異方性ソリッドと厚肉シェルに対して使用可能で、Plyxfemと共にシェルプロパティ/PROP/TYPE17で積層間材料破壊モデルとしても使用できます。この破壊モデルは/MAT/LAW12 (3D_COMP)、/MAT/LAW14 (COMPSO)および/MAT/LAW25 (COMPSH)および/MAT/LAW1 (ELAST)（Plyxfemと使用される場合のみ）と適合します。

フォーマット

定義

例（複合材）

#RADIOSS STARTER
/UNIT/1
unit for mat and failure
#              MUNIT               LUNIT               TUNIT
                   g                  mm                  ms
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/MAT/COMPSH/1/1
composite example
#              RHO_I
             .001506
#                E11                 E22                NU12     Iform                           E33
              144000               10000                 .25         0                         20000
#                G12                 G23                 G31              EPS_f1              EPS_f2
                4200                4200                4200                   0                   0
#             EPS_t1              EPS_m1              EPS_t2              EPS_m2                dmax
                   0                   0                   0                   0                   0
#              Wpmax               Wpref      Ioff                         ratio
             1000000                   0         0                             0
#                  b                   n                fmax
                   0                   0             1000000
#            sig_1yt             sig_2yt             sig_1yc             sig_2yc               alpha
               10100               10100               10100               10100                   0
#           sig_12yc            sig_12yt                c_12          Eps_rate_0       ICC
               10068               10068                   0                   0         0
#          GAMMA_ini           GAMMA_max               d3max
                   0                   0                   0
#  Fsmooth                Fcut
         0                   0
/FAIL/LAD_DAMA/1/1
#                 K1                  K2                  K3              Gamma1              Gamma2
                2000                2000                2000               1E-20               1E-20
#                 Y0                  YC                   K                   A             Tau_max
                  40                 160              100000                   1                 .01
# Ifail_sh  Ifail_so                                                                               
         1         3
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#enddata
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

コメント

1. 剥離用のLadeveze破壊損傷モデル：(1)
   $${\text{Y}}_{{\text{d}}_3}={\frac{\partial {\text{E}}_{\text{D}}}{\partial {\text{d}}_3}|}_{\sigma =\text{cst}}=\frac{1}{2}\frac{{\langle {\sigma }_{33}\rangle }^2}{{\text{K}}_3{\left(1-{\text{d}}_3\right)}^2}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\text{Mode}\hspace{0.17em}I$$
   (2)
   $${\text{Y}}_{{\text{d}}_2}={\frac{\partial {\text{E}}_{\text{D}}}{\partial {\text{d}}_2}|}_{\sigma =\text{cst}}=\frac{1}{2}\frac{{{\sigma }_{32}}^2}{{\text{K}}_2{\left(1-{\text{d}}_2\right)}^2}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\text{Mode}\hspace{0.17em}II$$
   (3)
   $${\text{Y}}_{{\text{d}}_1}={\frac{\partial {\text{E}}_{\text{D}}}{\partial {\text{d}}_1}|}_{\sigma =\text{cst}}=\frac{1}{2}\frac{{{\sigma }_{31}}^2}{{\text{K}}_1{\left(1-{\text{d}}_1\right)}^2}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\text{Mode}\hspace{0.17em}\mathrm{III}$$

   
   
   図 1.

   ²

   2次元ソリッド要素の場合は、モードIIとモードIIIのみを使用可能です。

   ここで、 $d_i$ はその破壊モードに伴う内部損傷パラメータです。

   ダメージ進展則は等価なダメージエネルギー解放比によってコントロールされます。

   $Y=Y_{d_3}+{\gamma }_1{Y}_{d_1}+{\gamma }_2{Y}_{d_2}$ ここで、 $Y_{d_i}|{}_t=\sup Y_{d_i}|{}_{\tau \le t}$

   損傷パラメータの進展は連成係数 ${\gamma }_1$ と ${\gamma }_2$ に強く連成します。これらの2つのパラメータは剥離試験から来ます。

   現在の破壊モデルについて、 $d_1=d_2=d_3=d$ であるとみなします。

   損傷値dは特定の速度で増大します：(4)

   $$\dot{d}=\frac{k}{a}\left[1-\exp \left(-a\langle w\left(Y\right)-d\rangle \right)\right]$$

   $d<1$ の場合。

   それ以外の場合は、 $d=1$

   このとき、

   $a$ は、延性破壊の尺度です。値が低いほど、破壊の延性が強くなります。

   $\frac{a}{k}$ は、最小破壊期間です。 $Y_0$ と $Y_c$ の間のエネルギーの期間は少なくとも $\frac{a}{k}$ に等しくなければなりません。(5)

   $$\langle x\rangle =\{\begin{array}{c}x\mathit{if}x>0\\ 0\mathit{if}x<0\end{array}$$
   関数 $w(Y)$ は次のように計算されます：(6)

   $$w(Y)=\frac{\langle \sqrt{Y}-\sqrt{Y_0}\rangle }{\sqrt{Y_c}-\sqrt{Y_0}}$$

   損傷パラメータ $d\le 1.0$ ,の場合は、応力 ${\sigma }_{33}$ 、 ${\sigma }_{13}$ および ${\sigma }_{23}$ は次の関数に従って減少します：

   緩和のテクニックが応力を徐々に減少させるために使用されます：(7)

   $$\sigma (t)=\mathrm{f}(t)\cdot {\sigma }_d(t_r)$$
   ここで、(8)

   $$\mathrm{f}(t)=\exp \left(-\frac{t-t_r}{{\tau }_{\max }}\right)$$
   ここで、

   $t$

   時間

   $t_r$

   損傷基準が推定される場合における緩和の開始時間

   ${\tau }_{\text{max}}$

   動的緩和の時間

   ${\sigma }_d\left(t_r\right)$

   損傷開始時の応力

2. fail_IDは、/STATE/BRICK/FAILおよび/INIBRI/FAILで使用されます。デフォルト値はありません。この行が空白の場合、/INIBRI/FAIL内の破壊モデル変数のために出力される値はありません（/STATE/BRICK/FAILオプションで.staファイルに書き込まれます）。
3. 破壊基準に達した後、 ${\tau }_{\max }$ の値は破断要素での応力が徐々に減少して0になるまでの期間を決定します。応力が破壊の開始時における応力値の1%に達すると、要素は削除されます。これは、要素が突然削除され、隣接要素で破断の“連鎖反応“が起こることによる不安定性を回避するために必要となります。破壊基準に達した場合でも、 ${\tau }_{\max }=1.0E30$ のデフォルト値であれば要素は削除されません。したがって、 ${\tau }_{\max }$ をシミュレーションの時間ステップの10倍大きく定義することが推奨されます。