/FAIL/JOHNSON

フォーマット

定義

例

3つの一般的な荷重（例： 純引張、純せん断、または純圧縮）では、材料パラメータD₁、D₂、D₃は、上記の簡素化された応力-ひずみ関係を使用して決定できます。

図 1.

例1（鋼材）

#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
#              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/JOHNSON/1/1
#                 D1                  D2                  D3                  D4                  D5
                0.11                0.08                -1.5                   0                   0
#              EPS_0  Ifail_sh  Ifail_so                                    Dadv               Ixfem
                   1         1         1                                       0                   0
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#enddata
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

例2（鋼材）

#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
#              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/JOHNSON/1/1
#                 D1                  D2                  D3                  D4                  D5
                0.11                0.08                -1.5                   0                   0
#              EPS_0  Ifail_sh  Ifail_so                                    Dadv               Ixfem
                   1         1         1                                     0.5                   1
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#enddata
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

コメント

1. パラメータは、次のような応力とひずみとの関係で使用されます：(2)
   $${\epsilon }_f=\left[D_1+D_2\text{exp}\left(D_3{\sigma }^{*}\right)\right]\hspace{0.17em}\left[1+D_4\mathit{ln}\left({\dot{\epsilon }}^{*}\right)\right]\hspace{0.17em}\left[1+D_5{T}^{*}\right]$$
   ここで、 ${\sigma }^{*}=\frac{{\sigma }_m}{{\sigma }_{VM}}$ （ ${\sigma }^{*}$ は応力軸性）(3)

   $${\dot{\epsilon }}^{*}=\frac{\dot{\epsilon }}{{\dot{\epsilon }}_0}$$
   $T*$ はすべての材料について次のように計算されます：(4)

   $$T^{*}=\frac{T-T_r}{T_{melt}-T_r}$$
   ここで、

   T_r

   初期温度

   T_melt

   材料LAW2およびLAW4の溶融温度

   /HEAT/MAT（I_form =1）がこの材料モデルを参照している場合、このカード内で定義されたT_rとT_meltの値が、/HEAT/MATで定義された対応するT₀とT_meltで上書きされます。

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

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

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

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

   • /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. この破壊モデルには、2つの異なる破壊（破断または亀裂）が導入されています。破壊基準は次のように計算されます：
   • 要素の破断Ixfem=0）：

     D > 1の場合、要素の破断（削除される）

   • 要素の亀裂（Ixfem=1）：

     次の場合に要素に亀裂が入ります：

     その要素に破断した隣接要素がない場合、新しい亀裂の初期化。

     その要素に破断した隣接要素があり、D > D_advの場合、破断は進展し、破断の進展にD_advが使用されます。_adv を使用（既存の亀裂が要素の境界に達した場合）。

     2つ目の亀裂が同じ要素に達した場合は、その要素が削除されます。

     ここで、(5)

     $$D=\sum \frac{\text{Δ}{\epsilon }_p}{{\epsilon }_f}\ge 1$$

   Dadvは常に1より小さくなくてはなりません。