Iform = 2

ブロックフォーマットキーワードこの境界材料を使用すると、グローバル材料状態を計算するためにも使用される副材料状態（密度、エネルギー、および体積比率）を適用できます。副材料EOSパラメータは、（ドメインの）隣接要素のパラメータで構成される必要があります。

フォーマット

(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)	(10)
/MAT/LAW51/`mat_ID`/`unit_ID`
`mat_title`
空白のフォーマット
`I`_form

#グローバルパラメータ

(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)	(10)
`Scale`_time		`P`_EXT				$V E L_{i n}$			$f c t_I D_{V E L}$

#材料1パラメータ

(1)	(3)	(5)	(7)	(8)	(9)
$α_{0}^{mat_1}$	$ρ_{0}^{mat_1}$	$E_{0}^{m a t_ 1}$	`fct_ID`_{$α$ 1}	`fct_ID`_{$ρ$ 1}	$f c t_I D_{E 1}$
$C_{1}^{m a t_1}$	$C_{2}^{m a t_1}$	$C_{3}^{m a t_1}$	$C_{4}^{m a t_1}$		$C_{5}^{m a t_1}$
$Δ P_{min}^{mat_1}$	$C_{0}^{m a t_1}$

#材料2パラメータ

(1)	(3)	(5)	(7)	(8)	(9)
$α_{0}^{mat_2}$	$ρ_{0}^{mat_2}$	$E_{0}^{m a t_ 2}$	`fct_ID`_{$α$ 2}	`fct_ID`_{$ρ$ 2}	`fct_ID`_E2
$C_{1}^{m a t_2}$	$C_{2}^{m a t_2}$	$C_{3}^{m a t_2}$	$C_{4}^{m a t_2}$		$C_{5}^{m a t_2}$
$Δ P_{min}^{mat_2}$	$C_{0}^{m a t_2}$

#材料3パラメータ

(1)	(3)	(5)	(7)	(8)	(9)
$α_{0}^{mat_3}$	$ρ_{0}^{mat_3}$	$E_{0}^{m a t_ 3}$	`fct_ID`_{$α$ 3}	`fct_ID`_{$ρ$ 3}	`fct_ID`_E3
$C_{1}^{m a t_3}$	$C_{2}^{m a t_3}$	$C_{3}^{m a t_3}$	$C_{4}^{m a t_3}$		$C_{5}^{m a t_3}$
$Δ P_{min}^{mat_3}$	$C_{0}^{m a t_3}$

定義

フィールド	内容	SI単位の例
`mat_ID`	材料識別子（整数、最大10桁）
`unit_ID`	単位識別子。（整数、最大10桁）
`mat_title`	材料のタイトル（文字、最大100文字）
`I`_form	定式化フラグ（整数） =2：強制状態
`Scale`_time	入力の関数の横軸のスケールファクタデフォルト = 1（実数）
$V E L_{i n}$	流入速度のスケールファクタ 5 （実数）	$[\frac{m}{s}]$
`fct_ID`_VEL	（オプション）右記の速度関数の識別子； $f_{V E L} (t)$ 5 = 0： $f_{ρ i} (t)$ > 0： $v (t) = V E L_{i n} \cdot f_{V E L} (t)$ （整数）
$α_{0}^{mat_ i}$	初期強制体積比率 2 （実数）
$ρ_{0}^{mat_ i}$	初期強制密度（実数）	$[\frac{kg}{m^{3}}]$
$E_{0}^{m a t_ i}$	単位体積あたりの初期強制エネルギー（実数）	$[\frac{J}{m^{3}}]$
`fct_ID`_{$α$ i}	（オプション）右記の体積比率スケーリング関数の識別子； $f_{α i} (t)$ 3 = 0： $α^{m a t_i} (t) = α_{0}^{m a t_i}$ > 0： $α^{m a t_i} (t) = α_{0}^{m a t_i} \cdot f_{α_{i}} (t)$ （整数）
`fct_ID`_{$ρ_{i}$}	（オプション）右記の密度スケーリング関数の識別子； $f_{ρ i} (t)$ 3 = 0： $ρ^{m a t_i} (t) = ρ_{0}^{m a t_i}$ > 0： $ρ^{m a t_i} (t) = ρ_{0}^{m a t_i} \cdot f_{ρ_{i}} (t)$ （整数）
`fct_ID`_Ei	（オプション）右記の密度エネルギースケーリング関数の識別子； $f_{E i} (t)$ 3 = 0： $E^{{mat}_{i}} (t) = E_{0}^{mat_i}$ > 0： $E^{m a t_i} (t) = E_{0}^{m a t_i} \cdot f_{E_{i}} (t)$ （整数）
$C_{1}^{m a t_i}$	多項式EOSの係数（実数）	$[Pa]$
$C_{2}^{m a t_i}$	多項式EOSの係数（実数）	$[Pa]$
$C_{3}^{m a t_i}$	多項式EOSの係数（実数）	$[Pa]$
$C_{4}^{m a t_i}$	多項式EOSの係数（実数）
$C_{5}^{m a t_i}$	多項式EOSの係数（実数）
$Δ P_{\min}^{mat_ i}$	カットオフ圧力 4 デフォルト = -10^-30（実数）	$[Pa]$
$C_{0}^{m a t_i}$	多項式EOSの係数（実数）	$[Pa]$

この定式化は、ユーザーデータからの副材料状態を適用します。
体積比率： (1)
$α^{m a t_i} (t)$

密度： (2)
$ρ^{m a t_i} (t)$

密度エネルギー： (3)
$E^{m a t_i} (t)$

これを使用して、特定の多項式EOSから圧力を計算できます：(4)
$Δ P (t) = \min {Δ P_{\min} , C_{0} + C_{1} μ + C_{2}^{'} μ^{2} + C_{3}^{'} μ^{3} + (C_{4} + C_{5} μ) E (t)}$

ここで、 $P = Δ P + P_{E X T}$ かつ $μ (t) = \frac{ρ (t)}{ρ_{0}} = 1, E (t) = E_{int} (t) / V_{0}, C_{2}^{'} = C_{2} δ_{μ \geq 0}, and C_{3}^{'} = C_{3} δ_{μ \geq 0}$ 。これは、EOSが膨張に対して線形であり、圧縮に対して3次式であることを意味します。

また、グローバル材料状態は以下によって求められます：

圧力

$Δ P = \sum_{i} α^{m a t_i} (t) \cdot Δ P^{m a t_i}$

密度

$ρ = \sum_{i} α^{m a t_i} (t) \cdot ρ^{m a t_i}$

エネルギー

$\frac{E_{i n t}}{V} = \sum_{i} α^{m a t_i} (t) \cdot E^{m a t_i}$
体積比率によって、要素体積を3つの異なる材料で分け合うことができます。
材料ごとに、 $α_{0}^{mat_ i}$ を0と1の間に定義する必要があります。

初期体積比率の合計 $\sum_{i = 1}^{3} α_{0}^{mat _ i}$ は1に等しい必要があります。

体積の自動初期比率については、/INIVOLをご参照ください。
関数が定義されていない場合は、関連する量が一定になり、初期値に設定されます。ただし、入力量は、指定された関数識別子を使用して時間依存関数として定義できます。横軸関数は、 $f (t)$ ではなく $f (S c a l e_{t i m e}, t)$ の使用につながるFscale_tパラメータを使用してスケーリングすることもできます。
$Δ P_{\min}^{mat_ i}$ フラグは、計算される圧力の最小値です。
$P = Δ P + P_{E X T}$ のため、 $P_{E X T} = 0$ の定義は $Δ P = P$ および $Δ P_{\min} = P_{\min}$ を意味します。

液体材料圧力を正のままにして、引張り強度を避ける必要があります。これにより、

$P_{\min} = 0$ または $Δ P_{\min} = - P_{E X T}$

ソリッド材料については、 $Δ P_{\min}^{mat_ i}$ = 10³⁰ のデフォルト値が適切です。
速度が定義されていない場合、ユーザーは/IMPVELを節点と用いて定義する必要があります。もしくは、法線速度を入力します。法線速度は、各サブ材料に使用される同じ速度でグローバル材料に適用されます。停滞点における状態のみが既知の場合、代わりに、/MAT/LAW51、I_form=4,5を使用します。

Iform = 2

フォーマット

定義

コメント