/MAT/LAW42 (OGDEN)

ブロックフォーマットキーワードこのキーワードは、Ogden、Mooney-Rivlin材料モデルを使用して指定された超弾性、粘性、および非圧縮性の材料を定義します。

この法則は、通常、非圧縮性のゴム、ポリマー、フォーム、およびエラストマーのモデル化に使用されます。この材料は、シェル要素とソリッド要素に使用できます。

フォーマット

(1)	(3)	(5)	(6)	(7)	(9)	(10)
/MAT/LAW42/`mat_ID`/`unit_ID`または/MAT/OGDEN/`mat_ID`/`unit_ID`
`mat_title`
$ρ_{i}$
$ν$	$σ_{cut}$		`fct_ID`_blk	`Fscale`_blk	`M`	`I`_form
$μ_{1}$	$μ_{2}$	$μ_{3}$		$μ_{4}$	$μ_{5}$
空白のフォーマット
$α_{1}$	$α_{2}$	$α_{3}$		$α_{4}$	$α_{5}$
空白のフォーマット

M > 0の場合

(1)	(3)	(5)	(7)	(9)
`G`₁	`G`₂	`G`₃	`G`₄	`G`₅
...`M`個のGの値（1行あたり5つ）
$τ_{1}$	$τ_{2}$	$τ_{3}$	$τ_{4}$	$τ_{5}$
...`M`個の $τ$ の値（1行あたり5つ）

定義

フィールド	内容	SI単位の例
`mat_ID`	材料識別子（整数、最大10桁）
`unit_ID`	単位識別子（整数、最大10桁）
`mat_title`	材料のタイトル（文字、最大100文字）
$ρ_{i}$	初期密度（実数）	$[\frac{kg}{m^{3}}]$
$ν$	ポアソン比 3 デフォルト = 0.495（実数）
$σ_{cut}$	引張におけるカットオフ応力。デフォルト = 10³⁰（実数）	$[Pa]$
`fct_ID`_blk	体積係数を相対体積の関数としてスケーリングする関数識別子 6 （整数）
`Fscale`_blk	体積弾性関数のスケールファクターデフォルト = 1.0（実数）
`M`	Prony級数内の粘性項の数（Maxwellモデルの順序）（整数）
`I`_form	シェル要素に対する非圧縮性の定式化フラグ。このフラグは、ソリッド要素に対しては効果がありません。 = 0（デフォルト）シェル要素の法線方向のストレッチを計算して、要素の非圧縮性を提供します。 = 1 ソリッド要素に対して使用するのと同様の定式化を使用します。
$μ_{1}$	せん断超弾性係数の第1パラメータ（実数）	$[Pa]$
$μ_{2}$	せん断超弾性係数の第2パラメータ（実数）	$[Pa]$
$μ_{3}$	せん断超弾性係数の第3パラメータ（実数）	$[Pa]$
$μ_{4}$	せん断超弾性係数の第4パラメータ（実数）	$[Pa]$
$μ_{5}$	せん断超弾性係数の第5パラメータ（実数）	$[Pa]$
$α_{1}$	第1材料指数（実数）
$α_{2}$	第2材料指数（実数）
$α_{3}$	第3材料指数（実数）
$α_{4}$	第4材料指数（実数）
$α_{5}$	第5材料指数（実数）
`G`_i	`ith`番目の乗数（Prony粘性項の） 7 （実数）	$[Pa]$
$τ_{i}$	`ith`番目の時間緩和（Prony粘性項の）（実数）	$[s]$

例（ゴム材）

#RADIOSS STARTER
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/UNIT/1
unit for mat
                  kg                  mm                  ms
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#-  2. MATERIALS:
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/MAT/OGDEN/1/1
rubber
#              RHO_I
                1E-6
#                 Nu           sigma_cut           funIDbulk         Fscale_bulk         M     Iform
                .495                   0                   0                   0         0         0
#               Mu_1                Mu_2                Mu_3                Mu_4                Mu_5
                2e-3               -1e-3                   0                   0                   0
# blank card

#            alpha_1             alpha_2             alpha_3             alpha_4             alpha_5
                   2                  -2                   0                   0                   0
# blank card

#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#enddata
/END
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

この材料モデルは、Ogden、Neo-Hookean、またはMooney-Rivlin材料モデルを使用して指定された超弾性、粘性、および非圧縮性の材料を定義します。この法則は、通常、非圧縮性のゴム、ポリマー、フォーム、およびエラストマーのモデル化に使用されます。この材料は、シェル要素とソリッド要素に使用できます。
LAW42は、以下のOgden材料モデルのひずみエネルギー密度表現を使用します。(1)
$W (λ_{1}, λ_{2}, λ_{3}) = \sum_{p = 1}^{5} \frac{μ_{p}}{α_{p}} ({\bar{λ}}_{1}^{α_{p}} + {\bar{λ}}_{2}^{α_{p}} + {\bar{λ}}_{3}^{α_{p}} - 3) + \frac{K}{2} {(J - 1)}^{2}$

ここで、

$W$

ひずみエネルギー密度

$λ_{i}$

ith主工学ストレッチ

$J$

相対体積は次のように定義されます： $J = λ_{1} \cdot λ_{2} \cdot λ_{3} = \frac{ρ_{0}}{ρ}$

${\bar{λ}}_{i} = J^{- \frac{1}{3}} λ_{i}$

偏差ストレッチ

$α_{p}$ および $μ_{p}$

材料定数係数ペア。

５つまでの材料定数ペアを定義できます。

初期せん断弾性率 $μ$ および体積弾性率（ $K$ ）は、次のように与えられます：(2)
$μ = \frac{\sum_{p = 1}^{5} μ_{p} \cdot α_{p}}{2}$

および(3)
$K = μ \cdot \frac{2 (1 + ν)}{3 (1 - 2 ν)}$

ここで、 $ν$ はポアソン比で、体積弾性率の計算のみに使用されます。
パラメータ $α_{p}$ および $μ_{p}$ は、初期せん断弾性率 $μ > 0$ となるように選択されなくてはなりません。
材料の安定性のためには、各材料定数ペアが $μ_{p} \cdot α_{p} > 0$ である必要があります。
ポアソン比 $ν$ は、体積弾性率の計算のみに使用されます式 3。
純粋な非圧縮性材料の場合は、 $ν = 0.5$ です。この値は、体積弾性率（ $K$ ）の無限大値を意味します。そのため、非圧縮性材料に推奨される最大ポアソン比は、 $ν = 0.495$ です。

陰解法および陽解法シミュレーションの場合は、ポアソン比を高くすると、時間ステップ値の減少や発散につながる可能性があります。
Ogden材料モデルの特殊なケースが、ひずみエネルギー密度関数の次の式を使用して表現可能な非圧縮性Mooney-Rivlinモデルです：(4)
$W = C_{10} (I_{1} - 3) + C_{01} (I_{2} - 3)$

ここで、 $I_{1}$ と $I_{2}$ は、右Cauchy-Greenテンソルの第1不変量と第2不変量です。

この表現は、次の場合に、Ogdenひずみエネルギー密度関数から求められます：(5)
$μ_{1} = 2 \cdot C_{10}$
(6)
$μ_{2} = - 2 \cdot C_{01}$
(7)
$α_{1} = 2$
(8)
$α_{2} = - 2$
Ogden材料モデルのシンプルなケースが、ひずみエネルギー密度関数の次の式を使用して表現可能なNeo-Hookeanモデルです：(9)
$W = C_{10} (I_{1} - 3)$

ここで、

$I_{1}$

右Cauchy-Greenテンソルの第1不変量

$C_{10}$

材料の定数

この表現は、次の場合に、LAW42 Ogdenひずみエネルギー密度関数から求められます：(10)
$μ_{1} = 2 \cdot C_{10}$
(11)
$α_{1} = 2$

および(12)
$μ_{2} = α_{2} = 0$
材料の体積弾性率が圧縮を防ぐのに十分大きくない場合、LAW42では関数（fct_ID_blk）の入力が可能で、それにより体積弾性率を大きくして非圧縮性を保持できるよう、体積弾性率を相対体積の関数としてスケーリングします。
LAW42ではMaxwellモデル（剛性 $G_{i}$ とダンパ $η_{i}$ を持つn個のスプリングの系として簡単に記述することができる）を使用した粘性（速度）効果がモデル化されます：

図 1. Maxwellモデル

Maxwellモデルは、Prony級数の入力（ $G_{i}, τ_{i}$ ）を使って表現されます。超弾性初期せん断弾性率 $μ$ は、Maxwellモデルの長期せん断弾性率 $G_{\infty}$ と同じであり、 $τ_{i}$ は緩和時間です：(13)
$τ_{i} = \frac{η_{i}}{G_{i}}$

$G_{i}$ および $τ_{i}$ の値は正の値である必要があります。
粘性効果を含めるには、この材料則と共に/VISC/PRONYを使用する必要があります。

/MAT/LAW42 (OGDEN)

フォーマット

定義

例（ゴム材）

コメント