/FAIL/CHANG

フォーマット

定義

コメント

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

   $$e_f{}^2={\left(\frac{{\sigma }_{11}}{{\sigma }_1^t}\right)}^2+\beta {\left(\frac{{\sigma }_{12}}{{\overline{\sigma }}_{12}}\right)}^2$$
   圧縮繊維モード： ${\sigma }_{11}<0$ (2)

   $${e_c}^2={\left(\frac{{\sigma }_{11}}{{\sigma }_1^c}\right)}^2$$
3. マトリックス亀裂の破壊基準は：
   引張マトリックスモード： ${\sigma }_{22}>0$ (3)

   $$e_m{}^2={\left(\frac{{\sigma }_{22}}{{\sigma }_2^t}\right)}^2+{\left(\frac{{\sigma }_{12}}{{\overline{\sigma }}_{12}}\right)}^2$$
   圧縮マトリックスモード： ${\sigma }_{22}<0$ (4)

   $${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$$
4. 損傷パラメータ $e_f{}^2,e_c{}^2,e_m{}^2$ 、または ${e_d}^2\ge 1.0$ の場合、数値の不安定性を回避するため、指数関数を使用して応力を減少させます。緩和手法は、応力を徐々に減少させて使用します：(5)
   $$\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倍大きく定義することが推奨されます。