/FAIL/CHANG

フォーマット

定義

コメント

1. この破壊モデルはシェル専用です。
2. ここでは、方向1が繊維方向です。繊維破損の破壊基準は下記のように記述されます。
   引張繊維モード： ${\sigma }_{11}>0$

   • シェルの場合：
     $$e_f{}^2={\left(\frac{{\sigma }_{11}}{{\sigma }_1^t}\right)}^2+\beta {\left(\frac{{\sigma }_{12}}{{\overline{\sigma }}_{12}}\right)}^2$$
   • ソリッドの場合：
     $$e_f{}^2={\left(\frac{{\sigma }_{11}}{{\sigma }_1^t}\right)}^2+\beta {\left(\frac{{\sigma }_{12}}{{\overline{\sigma }}_{12}}\right)}^2+{\left(\frac{{\sigma }_{13}}{{\overline{\sigma }}_{12}}\right)}^2$$

   圧縮繊維モード： ${\sigma }_{11}<0$

   $${e_c}^2={\left(\frac{{\sigma }_{11}}{{\sigma }_1^c}\right)}^2$$

3. マトリックス亀裂の破壊基準は：
   引張マトリックスモード： ${\sigma }_{22}>0$

   • シェルの場合：
     $$e_m{}^2={\left(\frac{{\sigma }_{22}}{{\sigma }_2^t}\right)}^2+{\left(\frac{{\sigma }_{12}}{{\overline{\sigma }}_{12}}\right)}^2$$
   • ソリッドの場合：
     $$e_{m1}{}^2={\left(\frac{{\sigma }_{22}}{{\sigma }_2^t}\right)}^2+{\left(\frac{{\sigma }_{12}}{{\overline{\sigma }}_{12}}\right)}^2$$
     $$e_{m2}{}^2={\left(\frac{\sigma 33}{{\sigma }_2^t}\right)}^2+{\left(\frac{{\sigma }_{13}}{{\overline{\sigma }}_{12}}\right)}^2$$
   圧縮マトリックスモード： ${\sigma }_{22}<0$

   • シェルの場合：
     $${e_d}^2={\left(\frac{{\sigma }_{22}}{2{\overline{\sigma }}_{12}}\right)}^2+\left[{\left(\frac{{\sigma }_2^c}{2{\overline{\sigma }}_{12}}\right)}^2-1\right]\frac{{\sigma }_{22}}{{\sigma }_2^c}+{\left(\frac{{\sigma }_{12}}{{\overline{\sigma }}_{12}}\right)}^2$$
   • ソリッドの場合：
     $${e_{d1}}^2={\left(\frac{{\sigma }_{22}}{2{\overline{\sigma }}_{12}}\right)}^2+[{\left(\frac{{\sigma }_2^c}{2{\overline{\sigma }}_{12}}\right)}^2-1]\frac{{\sigma }_{22}}{{\sigma }_2^c}+{\left(\frac{{\sigma }_{12}}{{\overline{\sigma }}_{12}}\right)}^2$$
     $${e_{d2}}^2={\left(\frac{\sigma 33}{2{\overline{\sigma }}_{12}}\right)}^2+[{\left(\frac{{\sigma }_2^c}{2{\overline{\sigma }}_{12}}\right)}^2-1]\frac{\sigma 33}{{\sigma }_2^c}+{\left(\frac{{\sigma }_{13}}{{\overline{\sigma }}_{12}}\right)}^2$$
4. 損傷パラメータ $e_f{}^2,e_c{}^2,e_m{}^2$ 、または ${e_d}^2\ge 1.0$ の場合、数値の不安定性を回避するため、指数関数を使用して応力を減少させます。緩和手法は、応力を徐々に減少させて使用します：
   $$\sigma (t)=\mathrm{f}(t)\cdot {\sigma }_d(t_r)$$

   $\mathrm{f}(t)=\exp \left(-\frac{t-t_r}{{\tau }_{\max }}\right)$ と $t\ge t_r$ 。

   ここで、

   $t$

   時間

   $t_r$

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

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

   動的緩和の時間

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

   損傷開始時の応力

5. 損傷値Dは、 $0\le D\le 1$ です。破壊の状態は次のとおりです。
   • 次の場合、破壊なし： $0\le D<1$
   • 次の場合、破壊： $D=1$

   ここで、 $D=Max\left(e_f{}^2,e_c{}^2,e_m{}^2,e_d{}^2\right)$ 。この損傷値は/ANIM/SHELL/DAMAで示されます。

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

   MODEオプションなしで/H3D/ELEM/FAILURE/ID=FAILIDを使用すると、グローバルの破壊指数（全モード間の最大値に対応）をプロットできます。