/FAIL/ORTHBIQUAD

図 1. 双2次破壊基準形状

この2つの放物線破壊ひずみ曲線vs.応力の状態（応力の軸性）のパラメータは、入力c1-c5値を用いた初期化フェーズ中にRadiossによってインタラクティブに計算されます。

計算された放物線破壊ひずみ曲線は負の破壊ひずみ値を有します。これらの負の値は1E-6の破壊ひずみによって置き換えられ、これは、非常に大きな損傷累積と脆性の挙動をもたらします。

この破壊基準は、シェルおよびソリッドを含むすべての弾塑性材料にも使用できます。

ここで、 $Q_m^i$ はパラメータ $i$ の補間係数です。これらの補間係数は、Radioss Starterで自動的に計算されます。

たとえば、0度と90度の方向について2セットのパラメータが指定された場合、フーリエ級数補間によって45度の曲線を次のように求めることができます：

図 3. 直交異方性双2次基準形状の例

入力で指定されていない荷重方向は、フーリエ級数を使用して補間されます。2軸引張において、すべての方向の破壊時の塑性ひずみは同じです。

• 材料破壊の挙動が未知である場合、c1からc5までは0.0に設定され、軟鋼挙動（M-Flag=1）が使用されます。
• 引張破壊値のみが既知である場合、c3が定義されます（c1=c2=c4=c5 = 0.0）。軽度な鋼鉄挙動が使用され、ユーザー定義のc3値によってスケーリングされます。
• 材料挙動が既知である場合、M-Flagが定義され、c3は、期待される引張破壊に従った破壊モデルの調整に使用できます。選択された材料挙動は、ユーザー定義のc3値によってスケーリングされます。
• それ以外の場合、c1～c5が定義されることが望ましいが、定義されていない場合は0.0が使用されます。
  注: c5 = 0.0の場合、これはM-Flagの値に応じて自動的に計算されます。この場合、すべての入力方向で計算された中の最小値が保持されます。

• M-Flag > 0の場合、入力されたc1、c2、c4、c5の値は使用されず、テーブルからの事前定義済み速度値を使用してこれらの値が計算されます。

  $c1=r1\cdot c3$

  $c2=r2\cdot c3$

  $c3=c3$

  $c4=r4\cdot c3$

  $c5=r5\cdot c3$

• c3=0の場合は、テーブルに表示されているデフォルト値が使用されます。

  M-Flag	材料にほぼ対応	c3 （デフォルト）	r1	r2	r4	r5
  1	軟鋼	0.60	3.5	1.6	0.6	1.5
  2	HSS鋼	0.50	4.3	1.4	0.6	1.6
  3	UHSS鋼	0.12	5.2	3.1	0.8	3.5
  4	アルミニウムAA5182	0.30	5.0	1.0	0.4	0.8
  5	アルミニウムAA6082-T6	0.17	7.8	3.5	0.6	2.8
  6	塑性PA6GF30	0.10	3.6	0.6	0.5	0.6
  7	塑性PP T40	0.11	10.0	2.7	0.6	0.7
  99	ユーザー定義の値（オプション行）	0.30	オプション入力	オプション入力	オプション入力	オプション入力

  重要: Altairおよび執筆者のいずれも、これらの値から得られる結果、有効性、精度について一切の責任を負いません。ユーザーは、妥当な試験結果により自身の値を検証する必要があります。使用は、初期設計調査にのみ推奨されます。
• c3 > 0の場合、選択された材料の挙動は、c3およびr1からr5までの予め定義された速度値によってスケーリングされます。
• c5 = 0.0の場合は、すべての入力方向のうちの最小値が保持されます。
• M-Flag=99の場合、破壊ひずみ速度r1、r2、r4、r5は、その後に追加される行に入力される必要があります。

シェル要素については、積分点が $D=1$ に達すると、積分点の応力テンソルがゼロに設定されます。シェル要素は、P_thick_failの値に基づいて削除されます。

P_thick_failが空白または0に設定されている場合は、シェルプロパティのP_thick_failの値が使用されます。P_thick_fail > 0の場合は、シェルプロパティで定義されたP_thick_failの値は無視され、この破壊モデルに入力された値が使用されます。

P_thick_failが0より大きい値に設定されている場合は、板厚方向の破断した積分点の比率がP_thick_fail以上になると、要素が破断して、削除されます。

ソリッド要素では、積分点が $D=1$ に達すると、要素は削除されます。

• S-Flag = 1： 破壊曲線が作成されます（コメント1を参照）。この場合、これらの曲線は同じ応力3軸度の最小値に達しない可能性があります（図 4）。

  
  
  

  図 4. S-Flag =1の場合の直交異方性二重2次破壊基準形状

• S-Flag = 2： デフォルトで設定されます。これにより、c4の値がすべての方向の全体での最小値となります。すべての曲線についてこれを得るには、2つ目の式を、曲線の最小値がc4で生じる2つの2次サブ関数に分割します。ここで、 ${\sigma }^{*}=\frac{1}{\sqrt{3}}$ です（図 5）。

  
  
  

  図 5. S-Flag =2の場合の直交異方性二重2次破壊基準形状

• S-Flag = 3： S-Flag = 2と同じで、それに簡易化されたネッキング基準が追加されます（シェルの場合のみ）。ネッキング基準は、Marciniak-Kuczynski解析に基づきます。これは、応力の軸性 $\frac{1}{3}$ および $\frac{2}{3}$ の間のネッキングの開始を表す曲線を定義する2つの追加2次関数を使用します。この曲線の最小値は、Inst_startフィールドでユーザーによって定義され、平面ひずみ引張 ${\sigma }^{*}=\frac{1}{\sqrt{3}}$ で生じます（図 6）。

  この曲線を用いて、2番目のネッキング損傷値が計算され、破断は、すべての積分点が $D=1$ に達した際にのみ起こります。

  Inst_startの値は、単軸引張試験から、最大応力における（真に塑性の）ひずみとして推定できます。

  
  
  

  図 6. S-Flag =3の場合の直交異方性二重2次破壊基準形状

• ${\dot{\epsilon }}_0\ne 0$ かつ $C_{JC}\ne 0$ の場合は、Johnson-Cookのひずみ速度依存性が使用されます。この場合、破壊値における塑性ひずみは、ひずみ速度依存係数によって次のように乗算されます：(3)
  $$D={\displaystyle \sum _{t=0}^{\infty }\frac{\text{Δ}{\epsilon }_p}{{\epsilon }_f^{\theta }(\eta )\left(1+C_{JC}{〈\ln \frac{\dot{\epsilon }}{{\dot{\epsilon }}_0}〉}_{+}\right)}}$$
  ここで、

  ${\dot{\epsilon }}_0$

  非粘性極限ひずみ速度

  $C_{JC}$

  ひずみ速度依存パラメータ

  ${\langle \rangle }_{+}$

  Macaulayの括弧（正の値のみが考慮されます）

  注: /FAIL/ORTHBIQUAD基準と組み合わせたJohnson-Cook材料則を使用している場合は、構成則に使用されるJohnson-Cookパラメータは、破壊基準の場合と同じではない可能性があります。
• fct_ID_rate ≠ 0の場合は、ひずみ速度依存係数の表形式関数が使用されます。この場合、ひずみ速度に応じたひずみ速度係数（ $f_{rate}$ で表される）の進展を表す関数（/FUNCT）を定義する必要があります。この関数の横座標のひずみ速度スケールファクターを入力することもできます。Xscale_rate（デフォルトでは、このスケールファクターは1.0に設定されています）。表形式のひずみ速度依存性を使用して、損傷変数の計算は次のようになります：(4)
  $$D={\displaystyle \sum _{t=0}^{\infty }\frac{\text{Δ}{\epsilon }_p}{{\epsilon }_f^{\theta }(\eta )f_{rate}\left(\dot{\epsilon }\right)}}$$