/MAT/LAW21 (DPRAG)
ブロックフォーマットキーワード この材料則は、Drucker-Prager降伏基準に基づくもので、岩石-コンクリートのように内部摩擦を伴う材料のモデル化に使用されます。これらの材料の塑性挙動は、材料内の圧力に依存します。
フォーマット
(1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | (10) |
---|---|---|---|---|---|---|---|---|---|
/MAT/LAW21/mat_ID/unit_IDまたは/MAT/DPRAG/mat_ID/unit_ID | |||||||||
mat_title | |||||||||
ρi | |||||||||
E | ν | ||||||||
A0 | A1 | A2 | Amax | ||||||
fct_IDf | Kt | FscaleP | |||||||
ΔPmin | Pext | ||||||||
B | μmax |
定義
フィールド | 内容 | SI単位の例 |
---|---|---|
mat_ID | 材料識別子 (整数、最大10桁) |
|
unit_ID | 単位識別子 (整数、最大10桁) |
|
mat_title | 材料のタイトル (文字、最大100文字) |
|
ρi | 初期密度 (実数) |
[kgm3] |
E | ヤング率 (実数) |
[Pa] |
ν | ポアソン比 (実数) |
|
A0 | 材料塑性係数 (実数) |
[Pa2] |
A1 | 材料塑性係数 (実数) |
[Pa] |
A2 | 材料塑性係数 (実数) |
|
Amax | 制限的なフォンミーゼス応力。 デフォルトは1030に設定されます(実数) |
[Pa2] |
fct_IDf |
P(μ)
を記述する関数識別子 (整数) |
|
Kt | 引張体積弾性係数 3 (実数) |
[Pa] |
FscaleP | 圧力関数スケールファクター デフォルト = 1.0(実数) |
[Pa] |
ΔPmin | 最小圧力 デフォルト = -1030(実数) |
[Pa] |
Pext | 外部圧力 4 デフォルト = 0(実数) |
|
B | 除荷体積弾性係数 3 (実数) |
[Pa] |
μmax | 圧縮における最大体積ひずみ 5 (実数) |
▸例(Sand)
コメント
- 流体力学的挙動はユーザー定義の関数
P=f(μ)
で定義されます。ここで、
- P
- 材料内の圧力
- μ
- 右記の式で示される体積ひずみ; μ=ρρ0−1
図 1. - Drucker-Prager降伏基準では、修正フォンミーゼス降伏基準を使用して塊状構造の圧力の影響を取り込みます:
(1) F=J2−(A0+A1P+A2P2)図 2.ここで、- J2
- 偏差応力の第2不変量で、 σVM=√3J2
- P
- 圧力で、 P=−I13 で与えられる圧力( I1 は第1応力不変量)
- A0、A1、およびA2
- 材料塑性係数
- A1=A2=0
- 降伏基準はフォンミーゼス( σVM=√3A0 )
- 除荷体積弾性率Bは P(μ) を記述する関数の初期勾配に、引張体積弾性率Ktは除荷体積弾性率 Bの1/100に、Ktは正の値に、それぞれ設定することが推奨されます。
- 相対圧力定式化の場合は、外部圧力が必要です。この特定のケースでは、降伏基準とエネルギー統合に全圧力の値が必要とされます。Radiossは、
Pext
に相対する圧力を出力します。全圧力の値は以下から算出できます:
(2) P=Pext+ΔP全圧力の限界は以下から算出できます:(3) Pmin=Pext+ΔPminPext=0 の場合、出力結果が全圧力です:
P=ΔP および Pmin=ΔPmin
- Bは除荷体積弾性率Bが定義されている場合、
[0,μmax]
における勾配
dPdμ
よりも大きくする必要があります。
- B=0 で μmax=0 の場合、除荷経路と荷重経路は同じです。
- B=0 または μmax≠0 の場合、Bのデフォルト値は B=dPdμ|μmax です。
- B≠0 または μmax=0 の場合、 μmax のデフォルト値は B=dPdμ|μmax です。
図 3.