/EOS/COMPACTION

ブロックフォーマットキーワード圧縮されて除荷されることが可能な土やフォームなどの多孔質媒体について圧力を記述します。

フォーマット

(1)	(3)	(5)	(7)	(10)
/EOS/COMPACTION/`mat_ID`/`unit_ID`
`eos_title`
`C`₀	`C`₁	`C`₂	`C`₃	`I`_form
$μ_{\min}$	$μ_{\max}$	`B`
`P`_sh

定義

フィールド	内容	SI単位の例
`mat_ID`	材料識別子（整数、最大10桁）
`unit_ID`	単位識別子。（整数、最大10桁）
`eos_title`	状態方程式のタイトル（文字、最大100文字）
`C`₀	状態方程式のパラメータ（実数）
`C`₁	状態方程式のパラメータ（実数）
`C`₂	状態方程式のパラメータ（実数）
`C`₃	状態方程式のパラメータ（実数）
`I`_form	除荷挙動のフラグ定義。 5 0 デフォルトは1 1（デフォルト）一定の除荷勾配 2 線形の除荷勾配（整数）
$μ_{\min}$	弾性限界。 $μ < μ_{\min}$ 以下の場合圧縮なし。 $P_{e} = C_{0} + C_{1} μ_{\min} + C_{2} μ_{\min}^{2} + C_{3} μ_{\min}^{3}$ デフォルト = 0.0（実数）
$μ_{\max}$	最大圧縮。デフォルト = 10²⁰（実数）
`B`	除荷体積弾性係数（実数）	$[Pa]$
`P`_sh	圧力シフト（実数）	$[Pa]$

例

#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/UNIT/1
unit for mat
                   g                  mm                  ms
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/MAT/LAW06/6/1
law06
#              RHO_I        
             2.33e-3          
#                 NU                PMIN
                0.22               -0.02
/EOS/COMPACTION/6/1
law06-compaction
#                 C0                  C1                  C2                  C3               Iform
                1E-2               0.256               0.256                   1                   1
#              MUMIN               MUMAX                BUNL
                 0.0               0.115                1.44
#                PSH               
                   0                            
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#enddata

この状態方程式は、土やフォームといった多孔質媒体の圧力をモデル化するために使用されます。圧縮は、弾性限界 $μ = μ_{\min}$ より小である際は、弾性応答を擁します。

図 1.
塑性圧縮は、下記の曲線で定義されるパス方向です；(1)
$P = C_{0} + C_{1} μ + C_{2} μ^{2} + C_{3} μ^{3}$

ここで、

$P$

材料内の圧力

$μ$

右記の式で示される体積ひずみ； $μ = \frac{ρ}{ρ_{0}} - 1$

除荷と再載荷は、勾配がBである直線に沿って為されます。

P_sh

相対圧力定式化をモデル化するために用いられる圧力シフト値

図 2. 圧力 vs. 体積ひずみ
2019.1より前のバージョンでは、この状態方程式は/MAT/LAW10に組み込まれていました（定弾性限界 $μ_{\min} = 0$ ）。
Radiossにより流体力学的圧力の計算に用いられ、右記の材料則と適合性のある状態方程式。
- /MAT/LAW3 (HYDPLA)
- /MAT/LAW4 (HYD_JCOOK)
- /MAT/LAW6 (HYDROまたはHYD_VISC)
- /MAT/LAW10 (DPRAG1)
- /MAT/LAW12 (3D_COMP)
- /MAT/LAW49 (STEINB)
- /MAT/LAW102 (DPRAG2)
- /MAT/LAW103 (HENSEL-SPITTEL)
除荷経路は、デフォルトでは一定の勾配 $B$ に沿って伸びています。I_form=2を使用して、圧縮に伴って増大する可変勾配を定義できます。この勾配は、 $µ_{\min}$ におけるC₁から $μ_{\max}$ における $B$ に進展します。
図 3. I_form=1：一定勾配に沿った除荷 $B$

図 4. I_form=2：除荷勾配が次の範囲においてC₁から $B$ に進展： $[µ_{\min}, µ_{\max}]$

/EOS/COMPACTION

フォーマット

定義

例

コメント