/FAIL/WIERZBICKI

フォーマット

(1)	(3)	(4)	(5)	(7)	(9)
/FAIL/WIERZBICKI/`mat_ID`/`unit_ID`
`C`₁	`C`₂		`C`₃	`C`₄	`m`
`n`	`I`_{fail_sh}	`I`_{fail_so}	`I`_moy

オプションの行

(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)	(10)
`fail_ID`

定義

フィールド	内容	SI単位の例
`mat_ID`	材料識別子（整数、最大10桁）
`unit_ID`	単位識別子。（整数、最大10桁）
`C`₁	第1パラメータ（実数）
`C`₂	第2パラメータ（実数）
`C`₃	第3パラメータ（実数）
`C`₄	第4パラメータ（実数）
`m`	第5パラメータ（実数）
`n`	硬化指数（実数）
`I`_{fail_sh}	シェル破壊モデルフラグ。 = 1 シェルが削除されるのは、 $D = \sum \frac{Δ ε_{p}}{{\bar{ε}}_{f}} \geq 1$ に１つの積分点または層で該当する場合です。 = 2 各積分点において、応力テンソルがゼロに設定されるのは、 $D = \sum \frac{Δ ε_{p}}{{\bar{ε}}_{f}} \geq 1$ に該当する場合です。シェルが削除されるのは、 $D = \sum \frac{Δ ε_{p}}{{\bar{ε}}_{f}} \geq 1$ にすべての積分点または層で該当する場合です。（整数）
`I`_{fail_so}	3次元ソリッド破壊モデルフラグ。 = 1 ソリッド要素が削除されるのは、 $D = \sum \frac{Δ ε_{p}}{{\bar{ε}}_{f}} \geq 1$ に、1つの積分点で該当する場合です。 = 2 積分点ごとに、 $D = \sum \frac{Δ ε_{p}}{{\bar{ε}}_{f}} \geq 1$ の場合は偏差応力テンソルが消失します。（整数）
`I`_moy	破壊3次元モデル（3次元ソリッド）フラグ。 =0 ${\begin{matrix} \bar{η} = \frac{σ_{m}}{σ_{V M}} \\ \bar{ξ} = \frac{27 J_{3}}{2 σ_{V M}^{3}} \end{matrix}$ =1 ${\begin{matrix} \bar{η} = \frac{\int_{0}^{ε_{p}} \frac{σ_{m}}{σ_{V M}} d ε_{p}}{ε_{p}} \\ \bar{ξ} = \frac{\int_{0}^{ε_{p}} \frac{27 J_{3}}{2 σ_{V M}^{3}} d ε_{p}}{ε_{p}} \end{matrix}$ （整数）
`fail_ID`	破壊基準識別子。 2 （整数、最大10桁）

例（金属）

#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
Metal
#              RHO_I
              2.8E-9                   0
#                  E                  Nu
               60000                  .3
#                  a                   b                   n           EPS_p_max            SIG_max0
                 370                 190                  .4                   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/WIERZBICKI/1/1
#                 C1                  C2                  C3                  C4                   m
                0.87                1.77                0.21                0.56                   2
#                  n  Ifail_sh  Ifail_so      Imoy
                 0.2         0         1         0
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#enddata
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

破壊フラグで使用される ${\bar{ε}}_{f}$ の場合は、次のように定義されます：
(1)
${\bar{ε}}_{f} = {{\bar{ε}}_{\max}^{n} - [{\bar{ε}}_{\max}^{n} - {\bar{ε}}_{\min}^{n}] {(1 - {\bar{ξ}}^{m})}^{\frac{1}{m}}}^{\frac{1}{n}}$

図 1.

ここで、

${\bar{ε}}_{\max} = C_{1} e^{- C_{2} \bar{η}}$

${\bar{ε}}_{\min} = C_{3} e^{- C_{4} \bar{η}}$
3次元ソリッドの場合：
- I_moy = 0であれば、
- $\bar{η} = \frac{σ_{m}}{σ_{V M}}$ $\bar{ξ} = \frac{27 J_{3}}{2 σ_{V M}^{3}}$
- I_moy = 1であれば、
- $\bar{η} = \frac{\int_{0}^{ε_{p}} \frac{σ_{m}}{σ_{V M}} d ε_{p}}{ε_{p}}$ 、 $\bar{ξ} = \frac{\int_{0}^{ε_{p}} \frac{27 J_{3}}{2 σ_{V M}^{3}} d ε_{p}}{ε_{p}}$
シェルの場合：
- $\bar{η} = \frac{σ_{m}}{σ_{V M}}$ 、 $\bar{ξ} = - \frac{27}{2} \bar{η} ({\bar{η}}^{2} - \frac{1}{3})$
ここで、

$σ_{m}$

静水圧応力

$σ_{V M}$

フォンミーゼス応力

$J_{3}$

偏差応力の第3不変量
fail_IDは、/STATE/BRICK/FAILおよび/INIBRI/FAILで使用されます。デフォルト値はありません。この行が空白の場合、/INIBRI/FAIL内の破壊モデル変数のために出力される値はありません（/STATE/BRICK/FAILオプションで.staファイルに書き込まれます）。

/FAIL/WIERZBICKI

フォーマット

定義

例（金属）

コメント