/MAT/LAW110 (VEGTER)
ブロックフォーマットのキーワード ひずみ速度依存性と熱効果を考慮し、Corus-Vegterの補間された降伏基準とVegterの硬化則を使用した弾塑性構成則。
フォーマット
(1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | (10) |
---|---|---|---|---|---|---|---|---|---|
/MAT/LAW110/mat_ID/unit_IDまたは/MAT/VEGTER/mat_ID/unit_ID | |||||||||
mat_title | |||||||||
E | v | Ires | |||||||
TAB_YLD | MAT_Xscale | MAT_Yscale | |||||||
Tini | Chard | Fcut | VP | Ismooth | TAB_TEMP |
(1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | (10) |
---|---|---|---|---|---|---|---|---|---|
R |
(1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | (10) |
---|---|---|---|---|---|---|---|---|---|
R |
(1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | (10) |
---|---|---|---|---|---|---|---|---|---|
(1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | (10) |
---|---|---|---|---|---|---|---|---|---|
R |
定義
フィールド | 内容 | SI単位の例 |
---|---|---|
mat_ID | 材料識別子 (整数、最大10桁) |
|
unit_ID | 単位の識別子(オプション) (整数、最大10桁) |
|
mat_title | 材料のタイトル (文字、最大100文字) |
|
初期密度。 (実数) |
||
E | ヤング率。 (実数) |
|
v | ポアソン比。 (実数) |
|
Ires | マッピングアルゴリズムフラグを返します。
(整数) |
|
Vegter定式化の選択:
(整数) |
||
TAB_YLD | 表形式の降伏応力 – 塑性ひずみ - ひずみ速度関数の識別子。 (整数) |
|
MAT_Xscale | 表形式の降伏 – 塑性ひずみ - ひずみ速度関数のXスケールファクター。 デフォルト = 10(実数) |
|
MAT_Yscale | 表形式の降伏 – 塑性ひずみ - ひずみ速度関数のYスケールファクター。 デフォルト = 1.0(実数) |
|
2軸スケールファクター。 (実数 > 0.0) |
||
RDに対して0度の方向の2軸ひずみ速度比率。 (実数 > 0.0) |
||
初期降伏応力 (実数) |
||
硬化応力増分。 (実数) |
||
大ひずみ硬化パラメータ。 (実数) |
||
微小ひずみ硬化パラメータ。 (実数) |
||
硬化指数 (実数) |
||
初期塑性ひずみ (実数) |
||
限界動的流れ応力。 (実数) |
||
最大活性化エンタルピー。 (実数) |
||
熱活性化運動の限界ひずみ速度。 (実数) |
||
ひずみ速度挙動指数。 (実数) |
||
Tini | 初期温度。 (実数) |
|
Chard | 硬化係数。
(実数) |
|
Fcut | ひずみ速度フィルタリングのカットオフ周波数。 デフォルト = 1.0 x 1020(実数) |
|
VP | ひずみ速度選択フラグ
(整数) |
|
Ismooth | 補間タイプ(表形式降伏関数の場合)。
(整数) |
|
TAB_TEMP | 表形式の降伏応力 – 塑性ひずみ - 温度の識別子。 (整数) |
|
単軸スケールファクター。 (実数 > 0.0) |
||
R | Lankford係数。 デフォルト = 1.0(実数 > 0.0) |
|
平面ひずみスケールファクターの第1成分。 (実数 > 0.0) |
||
平面ひずみスケールファクターの第2成分(
= 1)。 (実数 > 0.0) |
||
平面ひずみスケールファクターの第2成分を計算するための平均ファクター(
= 2)。 デフォルト = 0.5(実数 > 0.0) |
||
せん断倍率 (実数 > 0.0) |
||
回転方向(RD)を基準にして
度の方向の最大単軸工学応力(
= 3)。 (実数 > 0.0) |
||
回転方向(RD)を基準にして
度の方向の最大単軸均一伸張量(%単位)(
= 3)。 (実数 > 0.0) |
||
回転方向(RD)を基準にして
度の方向のLankford係数(
= 3)。 デフォルト = 1.0(実数 > 0.0) |
||
平面ひずみ重み係数(
= 4)。 (実数 > 0.0) |
||
せん断重み係数(
= 4)。 (実数 > 0.0) |
例
#RADIOSS STARTER
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/UNIT/25
Local unit system
Mg mm s
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#- 2. MATERIALS:
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/MAT/LAW110/1/25
Steel: Icrit = 1, example with 3 angles (0°, 45° and 90° to the RD)
# Init. dens.
7.85E-9
# E nu
194200.0 0.3
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
# Icrit TAB_YLD MAT_Xscale MAT_Yscale fBI rhoBI
1 0 0.0 0.0 1.004 0.889
# YLD0 DSIGM BETA OMEGA n
107.1 179.6 0.25 8.07 1.0
# EPS0 SIGS DG0 Deps0 m
0.0 20.0 800 3.61e-3 1.0
# TINI C_HARD F_CUT VP Ismooth TAB_TEMP
293.0 0.0 10000.0 1 1 0
# fUN_THETA R_THETA fPS1_THETA fPS2_THETA fSH_THETA
1.021 0.64 1.061 0.5305 0.560
0.987 0.48 1.037 0.5185 0.640
1.009 0.76 1.048 0.5240 0.560
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/MAT/LAW110/2/25
Steel :Icrit = 2, example with 3 angles (0°, 45° and 90° to the RD)
# Init. dens.
7.85E-9
# E nu
194200.0 0.3
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
# Icrit TAB_YLD MAT_Xscale MAT_Yscale fBI rhoBI
2 0 0.0 0.0 1.004 0.889
# YLD0 DSIGM BETA OMEGA n
107.1 179.6 0.25 8.07 1.0
# EPS0 SIGS DG0 Deps0 m
0.0 20.0 800 3.61e-3 1.0
# TINI C_HARD F_CUT VP Ismooth TAB_TEMP
293.0 0.0 10000.0 1 1 0
# fUN_THETA R_THETA fPS1_THETA ALPS_THETA fSH_THETA
1.021 0.64 1.061 0.5 0.560
0.987 0.48 1.037 0.5 0.640
1.009 0.76 1.048 0.5 0.560
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/MAT/LAW110/3/25
Steel
# Init. dens.
7.85E-9
# E nu
194200.0 0.3
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
# Icrit TAB_YLD MAT_Xscale MAT_Yscale fBI rhoBI
3 0 0.0 0.0 1.004 0.889
# YLD0 DSIGM BETA OMEGA n
107.1 179.6 0.25 8.07 1.0
# EPS0 SIGS DG0 Deps0 m
0.0 20.0 800 3.61e-3 1.0
# TINI C_HARD F_CUT VP Ismooth TAB_TEMP
293.0 0.0 10000.0 1 1 0
# RM_0 RM_45 RM_90 AG_0 AG_45
408.4 408.4 408.4 20.0 20.0
# AG_90 R_0 R_45 R_90
20.0 0.64 0.48 0.76
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/MAT/LAW110/4/25
Steel :example with 3 angles (0°, 45° and 90° to the RD)
# Init. dens.
7.85E-9
# E nu
194200.0 0.3
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
# Icrit TAB_YLD MAT_Xscale MAT_Yscale fBI rhoBI
4 0 0.0 0.0 1.004 0.889
# YLD0 DSIGM BETA OMEGA n
107.1 179.6 0.25 8.07 1.0
# EPS0 SIGS DG0 Deps0 m
0.0 20.0 800 3.61e-3 1.0
# TINI C_HARD F_CUT VP Ismooth TAB_TEMP
293.0 0.0 10000.0 1 1 0
# fUN_THETA R_THETA W_PS W_SH
1.021 0.64 0.4125 0.75
0.987 0.48 0.4125 0.75
1.009 0.76 0.4125 0.75
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#enddata
コメント
-
= 1の場合、この材料則では、次のように定義される従来のVegter降伏基準が使用されます:
(1) ここで、- 降伏関数。
- 降伏応力。
- 補間されたVegter相当応力。
平面応力条件下のVegter 降伏局面は、主応力空間で3つの2次ベジェ補間曲線を定義することで得られます。(図 1)。これらの曲線を使用して、異なる荷重条件(せん断、単軸引張、平面ひずみ、等2軸引張)で実験的に測定された4つの参照ポイントがリンクされます。2つの参照ポイント間で、降伏局面が次のように定義されます:(2) ここで、- 、
- 2つの参照ポイント。
- ヒンジポイント(自動的に計算されます)。
- 各時間ステップで計算される、降伏局面上の位置を特定するパラメータ。
測定された4つの参照ポイントは次のように表されます: 、 、 、 。これらの参照ポイントの座標はユーザーが指定します。
ここで、 = 、 、 です。回転方向を基準にして指定された方向について、せん断 、単軸引張 、等2軸引張の参照ポイント には、1つの成分のみを設定する必要があります。平面ひずみ参照ポイントについては、2つの成分 と をユーザーが指定する必要があります。
2つ目の座標 は自由に選択できます。この座標が設定されていない場合は、2つの隣接するヒンジポイントの2つ目の座標の平均が、この座標の値として使用されます。
2つの参照ポイントの間に配置されるヒンジポイントは、それらの参照ポイントの座標に加え、各参照ポイントの降伏曲面の法線も使用して計算されます。降伏曲面の法線は、法線則により、ひずみ速度テンソル成分を使用して次のように表すことができます:(3) ここで、 は指定された参照ポイントのひずみ速度比率です。せん断の参照ポイントについては、 です。平面ひずみの参照ポイントについては、 です。単軸引張の参照ポイントについては、この比率は、Lankford係数( )を使用して次のように計算できます:(4) Lankford係数はユーザーが指定する必要があります。最後に、等2軸引張 のひずみ速度比率も設定する必要があります。
-
= 2の場合は、標準のVegter基準が使用されます。これは従来のVegterと同じ基準ですが、入力カードが異なります。このカードで、平面ひずみポイントの2つ目の座標は、2つの隣接するヒンジポイントの2つ目の座標
および
の加重平均として次のように計算されます:
(5) - = 3の場合は、Vegter 2017降伏基準が使用されます。これは従来のVegterと同じ基準ですが、入力カードが異なります。このカード内のすべてのパラメータは、最大単軸工学応力 、最大均一伸張量 、およびLankford係数 によって決定されます。これら3つのパラメータは、回転方向(RD)を基準にして角度 度の3つの方向(0度、45度、90度)について指定する必要があります。この基準では、角度の数は3に固定されています。
-
= 4の場合、この材料則では、次のように定義される簡易化されたVegter Lite降伏基準が使用されます:
(6) ここで、- 降伏関数。
- 降伏応力。
- 補間されたVegter相当応力。
平面応力条件下のVegter 降伏局面は、主応力空間内で2つの2次Nurbs補間曲線を定義することで得られます。これらの曲線を使用して、異なる荷重条件(単軸圧縮、単軸引張、等2軸引張)で実験的に測定された3つの参照ポイントがリンクされます。2つの参照ポイント間で、降伏局面が次のように定義されます:(7) (8) ここで、- 、
- 2つの参照ポイント。
- ヒンジポイント(自動的に計算されます)。
- 各時間ステップで計算される、降伏局面上の位置を特定するパラメータ。
- ヒンジポイントに関連付けられた重み係数。
図 2では、降伏基準は、平面応力における従来のVegter定式化によって定義されています( = 2)。
従来のVegter定式化の場合、パラメータ 、 、 、 を設定する必要があります。この簡易化された定式化では、2つのヒンジポイントは純せん断状態と平面ひずみ状態にあります。2つの関連付けられた重み係数 と も定義する必要があります。
- 異方性を考慮するには、回転方向RDを基準にして異なる角度
のいくつかの
方向について一連のパラメータを定義します。これらの方向は、0と
の間に均等に分散されている必要があります。指定された複数の方向の間に配置されたすべての荷重方向について、Vegterモデルではフーリエ級数補間の使用が提示されています。各方向に対してすべてのパラメータを定義する必要はありません。パラメータ は、すべての方向について同じです。方向0(回転方向)のひずみ速度比率 だけで、他のすべての方向のこの比率を特定できます。
(9) 他のすべてのパラメータでは、次のフーリエ補間が使用されます:(10) 従来、標準、および2017のVegter定式化( = 1、2、3)では、同様の処理が行われます。(11) 簡易化されたVegter Lite定式化( = 4)では、同様の補間が 、 、 に対して行われます。この場合、次のようになります:(12) - パラメータによって凸状の降伏局面が定義される必要があります。定義されない場合は、シミュレーションが不安定になる可能性があります。このため、Radioss Starterには降伏基準の凸性を確認するためのデータチェックが用意されています。
- 1つの方向のみが定義されている場合( = 1)、この材料則は等方性であるため、/PROP/TYPE1(等方性シェル)と組み合わせて使用される必要があります。複数の方向が定義されている場合( > 1または = 3)、この材料則は異方性になるため、/PROP/TYPE9(直交異方性シェル)と組み合わせて使用される必要があります。
- 降伏応力
は次の式によって定義されます:
(13) ここで、- 初期降伏応力
- 硬化応力増分。
- 大ひずみ硬化パラメータ。
- 初期塑性ひずみ
- 微小ひずみ硬化パラメータ。
- 硬化指数
- 限界動的流れ応力。
- ボルツマン定数
- 温度
- 最大活性化エンタルピー。
- 非粘性極限ひずみ速度。
- ひずみ速度依存性指数。
- Tiniが定義されている場合、シミュレーション全体にわたって、温度はTiniの値で一定です。これにより、ひずみ速度依存性が一定になります(図 4a)。熱活性化効果(温度上昇後のひずみ速度依存性の微増)を考慮するために温度変化を計算する場合は、/HEAT/MATオプションを定義する必要があります(図 4a)。
- ひずみ速度
の計算は、次のようにフラグ
の値に依存します:
- = 1の場合: 塑性ひずみ速度が使用されます。
- = 2の場合: 全ひずみ速度が使用されます。
- = 3の場合: 偏差ひずみ速度が使用されます。
どの場合でも、ひずみ速度の計算には、カットオフ周波数 を使用したフィルタリングが含まれます。このカットオフ周波数は、次の関係を満たすようにユーザーによって定義されます:(14) ここで、 。
- 表形式の硬化降伏応力を使用する場合は、表形式関数TAB_YLDのIDを定義する必要があります。このテーブルを使用して、複数のひずみ速度で、塑性ひずみに応じたいくつかの降伏応力変化を定義できます。X方向とY方向に2つのスケールファクターを定義することもできます。この場合は、硬化パラメータ
、
、
、
、
は無視され、降伏応力は次のようになります:
(15) さらに、断熱条件での熱軟化を考慮するために、塑性ひずみと温度TAB_TEMPを使用して表形式の降伏応力を定義することもできます。この場合、降伏応力は次のようになります:(16) ここで、- Tini
- 参照温度。
- /HEAT/MATオプションを使用して計算される実際の温度。
- 係数Chardを設定することで、移動硬化を使用することもできます:
- Chard = 0の場合: 等方性硬化が使用されます。
- Chard = 1の場合: Prager-Ziegler移動硬化が使用されます。
- 0 ≤ Chard ≤ 1の場合: 等方性硬化と移動硬化の間でモデルが補間されます。