/FAIL/TBUTCHER

フォーマット

定義

例（鋼材）

#RADIOSS STARTER
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/UNIT/1
unit for mat
#              MUNIT               LUNIT               TUNIT
                  Mg                  mm                   s
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#-  1. MATERIALS:
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/MAT/PLAS_JOHNS/1/1
Steel + ductile failure model
#              RHO_I
              7.8E-9                   0
#                  E                  Nu
              210000                  .3
#                  a                   b                   n           EPS_p_max            SIG_max0
                 270                 450                  .6                   0                   0
#                  c           EPS_DOT_0       ICC   Fsmooth               F_cut               Chard
                   0                   0         0         0                   0                   0
#                  m              T_melt              rhoC_p                 T_r
                   0                   0                   0                   0
/FAIL/TBUTCHER/1/1
#              LAMDA                   K             SIGMA_R  Ifail_sh  Ifail_so    I_DUCT     Ixfem
                   2              0.3e+8                 150         0         1         0         0
#                  A                   B                                    Dadv
                   0                   0                                       0
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#enddata
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

コメント

1. Tuler-Butcher損傷累積パラメータが、指定された臨界損傷値Kを超えると、要素は破断します。損傷パラメータの累積は、延性材料と脆性材料では異なる方法で計算されます（I_duct = 1または2）。ソリッドの場合は、延性破断のみ使用でき、I_ductフラグは無効となります。
   延性材料の場合、損傷累積パラメータは下記のとおりです：(1)

   $$D={\displaystyle {\int }_0^t\max \left(0,{\left(\sigma -{\sigma }_r\right)}^{\lambda }\right)dt}>K$$
   ここで、

   ${\sigma }_r$

   初期破壊応力

   $\sigma$

   最大主応力

   $\lambda$

   材料定数（正の実数値）

   t

   要素に亀裂が入る時間

   D

   損傷積分

   K

   損傷積分の限界値

   脆性材料（シェルのみ）の場合、損傷パラメータは下記のとおりです。

   $\dot{D}=\frac{1}{K}{\left(\sigma -{\sigma }_r\right)}^a$

   ${\sigma }_r={\sigma }_0{\left(1-D\right)}^b$

   $D=D+\dot{D}\text{Δ}t$

2. XFEM定式化（Ixfem=1）は、Belytchko（I_shell=1または2）、I_shell=3または4、およびQEPH（I_shell=24）シェル要素とのみ適合性があります。XFEMフラグがアクティブ（Ixfem=1）な場合は、この破壊基準によって、要素または層が削除されるのではなく、要素に亀裂が入ります。
   単層と多層の2つのXFEMオプションが利用可能です。XFEMオプションは、材料識別子に適用された破壊基準に関連付けられているプロパティタイプに依存します。

   1. /PROP/SHELL（TYPE1）を使用する場合は、単層XFEMが適用されます。

      この場合、要素全体の厚さが単一層とみなされます。破壊基準は各積分点で計算されますが、この要素には亀裂が１つだけ現れます。このアプローチは、シェルフラグのすべての値（I_{fail_sh}=1または2）と適合性があります。亀裂方向は、最後に破断した積分点での主制約によって決定されます。

   2. PROP/SH_SANDW（TYPE11）を使用する場合は、多層XFEMが適用されます。

      この場合、板厚全体の各積分点は個別の層とみなされます。破壊基準は別々に計算されるため、亀裂方向は各層で異なる可能性があります。各層における亀裂方向は、ある要素から別の要素へ個別に伝わります。多層XFEMは、I_{fail_sh}=1とは適合性がありません。この場合、その値は自動でI_{fail_sh}=2に設定されます。

   3. /PROP/TYPE51を使用する場合は、多層XFEMが適用され、各層に別々の亀裂が現れ、要素間を個別に伝わります。このため、亀裂の方向とパターンは各層で異なります。破壊基準は各積分点で個別に計算され、層内の全積分点が破断した時点で亀裂が伝えられます。多層XFEMは、I_{fail_sh}=1とは適合性がありません。その値は自動でI_{fail_sh}=2に設定されます。

   警告: 単層および複数層XFEM定式化は、今のところ、同じモデル内では混在できません。それらのどちらかをモデル全体に対して選択する必要があります。
3. fail_IDは、/STATE/BRICK/FAILおよび/INIBRI/FAILで使用されます。デフォルト値はありません。この行が空白の場合、/INIBRI/FAIL内の破壊モデル変数のために出力される値はありません（/STATE/BRICK/FAILオプションで.staファイルに書き込まれます）。
4. 破壊基準は次のように計算されます。

   損傷累積パラメータDは、Ixfem=0の場合に使用される唯一の破断基準です。

   Ixfem=1の場合は、DとD_advの両方を、それぞれ亀裂の初期化および拡大のために使用できます。要素に破断した隣接要素がなく、基準を満たしている場合は、新しい亀裂を初期化するためにDが使用されます。そうでない場合は、既存の亀裂が要素の境界に発生すると、亀裂拡大基準であるD_advがDの代わりに使用されます。最終的に、2つ目の亀裂が同じ要素に達すると、この基準は削除されます。