/FAIL/HASHIN
ブロックフォーマットキーワード Hashin破壊モデルを記述します。この破壊モデルはシェルおよびソリッドで使用できます。
フォーマット
(1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | (10) |
---|---|---|---|---|---|---|---|---|---|
/FAIL/HASHIN/mat_ID/unit_ID | |||||||||
Iform | Ifail_sh | Ifail_so | ratio | I_Dam | Imod | I_frwave | |||
Sdel | Tcut |
(1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | (10) |
---|---|---|---|---|---|---|---|---|---|
Soft |
(1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | (10) |
---|---|---|---|---|---|---|---|---|---|
fail_ID |
定義
フィールド | 内容 | SI 単位の例 |
---|---|---|
mat_ID | 材料識別子 (整数、最大10桁) |
|
unit_ID | 単位識別子 (整数、最大10桁) |
|
Iform | 定式化フラグ
(整数) |
|
Ifail_sh | シェル破壊フラグ。
(整数) |
|
Ifail_so | ソリッド破壊フラグ
(整数) |
|
ratio | Isolid=2またはIfail_sh=2の場合、層(または積分点)の比よりも多くが破断すると要素が削除されます。 デフォルト = 1.0(実数) |
|
I_Dam | 損傷計算フラグ 6
(整数) |
|
Imod | 緩和時間計算。
(整数) |
|
I_frwave | 近傍の要素間の破壊伝播フラグ。
(整数) |
|
ひずみ率下限。 デフォルト = 0.0(実数) |
||
縦方向引張強度(繊維方向)。 デフォルト = 1020(実数) |
||
横方向引張強度(繊維方向に垂直)。 デフォルト = 1020(実数) |
||
板厚方向の引張り強度 デフォルト = 1020(実数) |
||
縦方向圧縮強度(繊維方向)。 デフォルト = 1020(実数) |
||
横方向圧縮強度(繊維方向に垂直)。 デフォルト = 1020(実数) |
||
破砕強度 デフォルト = 1020(実数) |
||
繊維せん断強度 デフォルト = 1020(実数) |
||
マトリックスせん断強度12 デフォルト = 1020(実数) |
||
マトリックスせん断強度23 デフォルト = 1020(実数) |
||
マトリックスせん断強度13 デフォルト = 1020(実数) |
||
マトリックスおよび剥離のCoulomb摩擦角 < 90° デフォルト = 0(実数) |
||
Sdel | 剥離基準スケールファクター デフォルト = 1.0(実数) |
|
動的緩和 5 デフォルト = 1020(実数) |
||
参照ひずみ速度 デフォルト = 10-20(実数) |
||
Tcut | ひずみ速度カットオフ期間 デフォルト = (実数) |
|
Soft | 近傍の要素のうちの1つがすでに破断している際は、減少係数が破壊基準に適用されます。 I_frwave、=2の場合のみ使用。 0.0.≤ Soft ≤ 1.0 デフォルト = 0.0(実数) |
|
fail_ID | (オプション)破壊基準識別子。 4 (整数、最大10桁) |
例(複合材)
#RADIOSS STARTER
/UNIT/1
unit for mat and failure
# MUNIT LUNIT TUNIT
kg mm ms
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/MAT/COMPSH/1/1
composite material
# RHO_I
1.5E-6
# E11 E22 NU12 Iform E33
42 40 .05 1 .5
# G12 G23 G31 EPS_f1 EPS_f2
3.4 3 3 0 0
# EPS_t1 EPS_m1 EPS_t2 EPS_m2 dmax
0 0 0 0 .9999
# Wpmax Wpref Ioff IFLAWP ratio
0 0 5 0 0
# c EPS_rate_0 alpha ICC_global
0 2E-4 0 1
# sig_1yt b_1t n_1t sig_1maxt c_1t
.1 25 .1 0 0
# EPS_1t1 EPS_2t1 SIGMA_rst1 Wpmax_t1
0 0 0 0
# sig_2yt b_2t n_2t sig_2maxt c_2t
.1 20 .1 0 0
# EPS_1t2 EPS_2t2 sig_rst2 Wpmax_t2
0 0 0 0
# sig_1yc b_1c n_1c sig_1maxc c_1c
.005 800 .5 0 0
# EPS_1c1 EPS_2c1 sig_rsc1 Wpmax_c1
.08 .15 .1 0
# sig_2yc b_2c n_2c sig_2maxc c_2c
.005 2000 .5 0 0
# EPS_1c2 EPS_2c2 sig_rsc2 Wpmax_c2
0 0 0 0
# sig_12yt b_12t n_12t sig_12maxt c_12t
.004 83 .31 0 0
# EPS_1t12 EPS_2t12 sig_rst12 Wpmax_t12
.075 .085 .05 0
# GAMMA_ini GAMMA_max d3max
1E31 1E31 .9999
# Fsmooth Fcut
0 0
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/FAIL/HASHIN/1/1
# Iform Ifail_sh Ifail_so Ratio I_Dam Imod Ifrwave EPS_DOT_MIN
2 1 0 0 1
# Sigma1_T Sigma2_T Sigma3_T Sigma1_C Sigma2_C
2 .525 1E30 1.7 1.7
# Sigma_C SigmaF_12 SigmaM_12 SigmaM_23 SigmaM_13
1E30 1E30 .075 1E30 1E30
# Phi Sdelam Tau_max EPS_DOT_0 Tcut
0 1 .01
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#enddata
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
コメント
- ratioの例: ratio=0.5、Ifail_sh=2(またはIfail_so=2)の場合、層(または積分点)の半分以上が破断すると要素が削除されます。
- 3次元材料破壊モデル:
- 一方向薄層モデル:引張 / せん断繊維モード:
(1) 圧縮繊維モード:(2) ここで、
クラッシュモード(3) ここで、
破壊マトリックスモード:(4) 剥離モード:(5) ここで、
注:(6) - 繊維薄層モデル:引張 / せん断繊維モード:
(7) (8) ここで、
圧縮繊維モード:(9) ここで、(10) ここで、
クラッシュモード(11) ここで、
せん断破壊マトリックスモード:(12) マトリックス破壊モード:(13) ここで、
損傷パラメータがFi ≥ 1.0の場合、数値の不安定性を回避するため、指数関数を使用して応力を減少させます。緩和手法は、応力を徐々に減少させて使用します:(14) ここで、(15) および
ここで、- 時間
- 損傷基準が推定される場合における緩和の開始時間
- 動的緩和の時間
- 損傷開始時の応力
- 一方向薄層モデル:
- 損傷値Dは、0 ≤ D ≤ 1です。破壊の状態は次のとおりです。
- 0 ≤ D > 1の場合: 破壊なし
- D = 1の場合: 破壊
ここで、一方向薄層モデルの場合は 、繊維薄層モデルの場合は 。この損傷値は/ANIM/BRICK/DAMAまたは/ANIM/SHELL/DAMAで示されます。
- fail_IDは、/STATE/BRICK/FAILおよび/INIBRI/FAILで使用されます。デフォルト値はありません。この行が空白の場合、/INIBRI/FAIL内の破壊モデル変数のために出力される値はありません(/STATE/BRICK/FAILオプションで.staファイルに書き込まれます)。
- 破壊基準に達した後、 の値は破断要素での応力が徐々に減少して0になるまでの期間を決定します。応力が破壊の開始時における応力値の1%に達すると、要素は削除されます。これは、要素が突然削除され、隣接要素』でa破断の“連鎖反応“が起こることによる不安定性を回避するために必要となります。破壊基準に達した場合でも、 のデフォルト値であれば要素は削除されません。したがって、 をシミュレーションの時間ステップの10倍大きく定義することが推奨されます。
- I_Damオプションは、損傷計算と損傷計算時の安定性を改善します。