/FAIL/RTCL

この理論は、ボイド成長モデリングに基づいています。この破壊モデルは、シェル要素とソリッド要素に使用できます。

フォーマット

(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)	(10)
/FAIL/RTCL/`mat_ID`/`unit_ID`
`EPScal`		`Inst`	`n`

オプションの行

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

定義

フィールド	内容	SI単位の例
`mat_ID`	材料識別子（整数、最大10桁）
`unit_ID`	（オプション）単位の識別子。（整数、最大10桁）
`EPScal`	キャリブレートされた単純引張破壊ひずみ $ε_{c a l}$ （シェルに対して正則化がアクティブ化されている場合は、 $L_{e}$ の参照メッシュサイズについて）。（実数）
`Inst`	シェルの損傷正則化をアクティブ化するフラグ。 = 0 デフォルト値。 = 1 ネッキング不安定性正則化はアクティブ化されません。 = 2（デフォルト）ネッキング不安定性正則化はアクティブ化されます。（整数）
`n`	シェル損傷正則化の硬化指数。（実数）
`fail_ID`	（オプション）破壊基準識別子。（整数、最大10桁）

例（アルミニウム）

#RADIOSS STARTER
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/UNIT/1
unit for mat and failure
#              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
Aluminum
#              RHO_I
              2.7E-9                   0
#                  E                  Nu     Iflag
               70000                  .3         0
#                  a                   b                   n           EPS_p_max            SIG_max0
                  90                 200                  .3                   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/RTCL/1/1
#             EPScal      Inst                   n         
                  .2         0                 .67                   
#  fail_ID
         1
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#enddata
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

この係数は、次のように応力軸性に従って計算されます：(1)
$D = \frac{1}{ε_{c r}^{f}} \int_{0}^{\infty} f_{R T C L} (η) d ε_{p}$

ここで、

$η$

次のように定義された応力軸性。 $\frac{σ_{m}}{σ_{V M}}$

$σ_{m}$

平均応力

$σ_{V M}$

フォンミーゼス相当応力。

$ε_{p}$

累積塑性ひずみ。

$ε_{c r}^{f}$

単純引張での破壊時の塑性ひずみ。

$f_{R T C L}$

計算が次のように定義される係数。
この係数は、次のように応力軸性に従って計算されます：
$f_{R T C L} = {\begin{matrix} 0 & if & η < - \frac{1}{3} \\ 2 \frac{1 + η \sqrt{12 - 27 η^{2}}}{3 η + \sqrt{12 - 27 η^{2}}} & if & - \frac{1}{3} \leq η < \frac{1}{3} \\ e^{- \frac{1}{2}} e^{\frac{3}{2} η} & if & η \geq \frac{1}{3} \end{matrix}$
ソリッド要素の破壊時の塑性ひずみ $ε_{c r}^{f} = ε_{c a l}$ 。ただし、シェル要素の場合は、次の2つのケースが発生する可能性があります：
- $I n s t = 1$ の場合、 $ε_{c r}^{f} = ε_{c a l}$
- $I n s t = 2$ の場合、 $ε_{c r}^{f} = n + (ε_{c a l} - n) \frac{t_{e}}{L_{e}}$
  ここで、
  
  $n$
  
  硬化指数（べき乗タイプ硬化を想定： $A + B ε_{p}^{n}$ ）
  
  $t_{e}$
  
  初期シェル厚。
  
  $L_{e}$
  
  シェル面積の平方根。
  
  この最後の式では、シェルのネッキング不安定性を考慮し、結果を正則化することができます。
  注: $t_{e} = L_{e}$ の場合は、キャリブレートされた値 $ε_{c a l}$ が発生します。
損傷は、出力リクエストDAMAを使用してアニメーションファイル内でポスト処理できます。シェル要素については、いずれかの積分点がD=1に達すると、積分点の応力テンソルがゼロに設定されます。板厚方向の破壊積分点の比率が、シェルプロパティで定義されたP_thick_failと等しくなると、要素が破壊されて削除されます。ソリッド要素では積分点がD=1に達すると、要素は削除されます。
fail_IDは、キーワード/INISHE/FAIL、/INISH3/FAIL、または/INIBRI/FAILを使用して、要素内の破壊初期化に使用されます。これらの値は、/STATE/SHELL/FAILまたは/STATE/BRICK/FAILオプションを使用して、.staファイルに書き込むことができます。

/FAIL/RTCL

フォーマット

定義

例（アルミニウム）

コメント