MATFVE

バルクデータエントリ周波数依存性粘弾性材料の材料プロパティを定義します。

フォーマット1（`TYPE`=FORMULA）

(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)	(10)
MATFVE	`MID`	`FORMULA`
	`R1`	`I1`	`a`	`R1k`	`I1k`	`b`

フォーマット2（`TYPE`=TABLE）

(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)	(10)
MATFVE	`MID`	`TABLE`
	`TID1`	`TID2`	`TID3`	`TID4`

フォーマット3（`TYPE`=PRONY）

(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)
MATFVE	`MID`	`PRONY`			`gD1`	`tD1`	`gB1`	`tB1`
	`gD2`	`tD2`	`gD3`	`tD3`	`gD4`	`tD4`	`gD5`	`tD5`
	`gB2`	`tB2`	`gB3`	`tB2`	`gB4`	`tB4`	`gB5`	`tB5`

フォーマット4（`TYPE`=PRELOAD）

(1)	(2)	(3)	(4)	(5)
MATFVE	`MID`	`PRELOAD`
	`EUNIAXI`
	`El`	`Es`	`f`	$ε$
	etc.
	`EVOLUME`
	`Kl`	`Ks`	`f`	`J`
	etc.

例

(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)	(10)
MATVE	2	PRONY			0.25	5e-2	0.25	5e-2

定義

フィールド	内容	SI単位の例
`MID`	固有の材料識別番号。デフォルトなし（整数 > 0）
`TYPE`	周波数依存性粘弾性材料のモデルタイプ。 FORMULA 4 TABLE 5 PRONY 6 PRELOAD 7 デフォルト値はありません。
`R1`	FORMULAタイプの偏差部分の実数部。デフォルト = 空白（実数）
`I1`	FORMULAタイプの偏差部分の虚数部。デフォルト = 空白（実数）
`a`	FORMULAタイプの偏差部分の実数。デフォルト = 空白（実数）
`R1k`	FORMULAタイプのバルク部分の実数部。デフォルト = 空白（実数）
`I1k`	FORMULAタイプのバルク部分の虚数部。デフォルト = 空白（実数）
`b`	FORMULAタイプのバルク部分の実数。デフォルト = 空白（実数）
`TID1`	偏差部分の実数部を指定するための表ID（`TABLEDi`のみ）。デフォルト = 空白（整数 > 0）
`TID2`	偏差部分の虚数部を指定するための表ID（`TABLEDi`のみ）。デフォルト = 空白（整数 > 0）
`TID3`	バルク部分の実数部を指定するための表ID（`TABLEDi`のみ）。デフォルト = 空白（整数 > 0）
`TID4`	バルク部分の虚数部を指定するための表ID（`TABLEDi`のみ）。デフォルト = 空白（整数 > 0）
`gDi`	$i$ 番目の偏差プロニー級数の弾性係数比。デフォルトなし（実数 > 0.0）
`tDi`	$i$ 番目の偏差プロニー級数の緩和時間。デフォルトなし（実数 > 0.0）
`gBi`	$i$ 番目のバルクプロニー級数の弾性係数比。デフォルトなし（実数 > 0.0）
`tBi`	$i$ 番目のバルクプロニー級数の緩和時間。デフォルトなし（実数 > 0.0）
`EUNIAXI`	`TYPE` = PRELOADの場合に、単軸試験からの貯蔵弾性係数と損失弾性係数を示すための継続行。
`EVOLUME`	`TYPE` = PRELOADの場合に、体積試験からの貯蔵弾性係数と損失弾性係数を示すための継続行。
`El`	単軸損失弾性係数。デフォルトなし（実数）
`Es`	単軸貯蔵弾性係数。デフォルトなし（実数）
`f`	周波数。デフォルトなし（実数 ≥ 0.0）
$ε$	単軸公称ひずみ。デフォルトなし（実数）
`Kl`	バルク損失弾性係数。デフォルトなし（実数）
`Ks`	バルク貯蔵弾性係数。デフォルトなし（実数）
`J`	現在の体積と元の体積の間の体積比。デフォルトなし（実数 > 0.0）

長期応答は、弾性材料の記述で指定され、これは同じMIDを持っている必要があります。MATHE/MAT1/MAT9がサポートされています。
CHEXA、CTETRA、CPENTA、およびCPYRA要素がサポートされています。
時間領域内の偏差緩和の関数 $g * (t)$ を考慮しています。(1)
$g * (t) = \sum_{i} \frac{G_{i}}{G_{\infty}} e^{- \frac{t}{τ_{i}}}$

サブスクリプト $i$ は、プロニー級数の $i$ 番目の成分を示します。

$G_{\infty}$

長期偏差弾性係数の入力。

$G_{i}$

偏差弾性係数比。

$τ_{i}$

緩和時間。

この関数のフーリエ変換は次のとおりです：(2)
$\tilde{g} * (ω) = F (g * (t)) = R + i I$

ここで、 $F$ はフーリエ変換です。式 2 は、周波数依存の偏差緩和です。同様に、周波数依存のバルク緩和の方程式を同じような形で記述できます。 (3)
$\tilde{k} * (ω) = F (k * (t)) = R_{k} + i I_{k}$

フーリエ変換の後で、実数部と虚数部 $R$ 、 $I$ および $R_{k}$ 、 $I_{k}$ は周波数依存となり、4とおりの方法（FORMULA、TABLE、PRONY、またはPRELOAD）で入力できます。
TYPE=FORMULAの場合： (4)
$R = R 1 * f^{- a}$
(5)
$I = I 1 * f^{- a}$
(6)
$R_{k} = R 1 k * f^{- b}$
(7)
$I_{k} = I 1 k * f^{- b}$

ここで、

$f$

周波数。

$R$ 、 $I$

偏差緩和関数のフーリエ変換の実数部と虚数部。

$R_{k}$ 、 $I_{k}$

バルク緩和関数のフーリエ変換の実数部と虚数部。
TYPE=TABLEの場合：

TID1

$ω R$ を周波数 $f$ の関数として指定します。

TID2

$ω I$ を周波数 $f$ の関数として指定します。

TID3

$ω R_{k}$ を周波数 $f$ の関数として指定します。

TID4

$ω I_{k}$ を周波数 $f$ の関数として指定します。

ここで、

$ω$

角周波数。

$R$ 、 $I$

偏差緩和関数のフーリエ変換の実数部と虚数部。

$R_{k}$ 、 $I_{k}$

バルク緩和関数のフーリエ変換の実数部と虚数部。
TYPE=PRONYの場合、データ入力はMATVEカードと同じ意味を持ちます。貯蔵弾性係数と損失弾性係数は、時間領域プロニー級数のフーリエ変換の後で、以下のように算定されます。 (8)
$G_{s} = (G_{\infty} + \sum_{i} \frac{G_{i} τ_{i}^{2} ω^{2}}{1 + τ_{i}^{2} ω^{2}})$
(9)
$K_{s} = (K_{\infty} + \sum_{i} \frac{K_{i} τ_{i}^{2} ω^{2}}{1 + τ_{i}^{2} ω^{2}})$
(10)
$G_{l} = (\sum_{i} \frac{G_{i} τ_{i} ω}{1 + τ_{i}^{2} ω^{2}})$
(11)
$K_{l} = (\sum_{i} \frac{K_{i} τ_{i} ω}{1 + τ_{i}^{2} ω^{2}})$

ここで、

$G_{s}$ 、 $G_{l}$

偏差貯蔵弾性係数および偏差損失弾性係数。

$K_{s}$ 、 $K_{l}$

バルク貯蔵弾性係数およびバルク損失弾性係数。

$G_{\infty}$ 、 $K_{\infty}$

長期材料プロパティ入力。

$G_{i}$ 、 $K_{i}$

弾性係数比。

$τ_{i}$

緩和時間。

$ω$

角周波数。

(12)
$ω I = 1 - \frac{G_{s}}{G_{\infty}}$
(13)
$ω R = \frac{G_{l}}{G_{\infty}}$
(14)
$ω I_{k} = 1 - \frac{K_{s}}{K_{\infty}}$
(15)
$ω R_{k} = \frac{K_{l}}{K_{\infty}}$

ここで、

$ω$

角周波数。

$R$ 、 $I$

偏差緩和関数のフーリエ変換の実数部と虚数部。

$R_{k}$ 、 $I_{k}$

バルク緩和関数のフーリエ変換の実数部と虚数部。
TYPE= PRELOADの場合：

EUNIAXI

単軸損失弾性係数、単軸貯蔵弾性係数、周波数 $f$ 、単軸公称ひずみを指定するための試験データ。

EVOLUME

バルク損失弾性係数、バルク貯蔵弾性係数、周波数 $f$ 、体積比（現行体積 / 初期体積）を指定するための試験データ。
すべてのタイプについて、偏差成分とバルク成分の指定を同時に空白にすることはできません。
直接法による周波数応答解析のみがサポートされています。
直接法による周波数応答解析の前に、直接法による周波数解析のベース状態を確立する非線形LGDISP解析を行う必要があります。直接法による周波数解析のサブケース情報セクション内のSTATSUB(PRELOAD)は、この非線形LGDISP解析を指している必要があります。荷重なしのダミーの非線形LGDISP解析は、ベース状態が初期コンフィギュレーションであることを意味します。
時間領域と周波数領域の粘弾性挙動を共に指定することはできません。すなわち、MATVEとMATFVEを共に使用することはできません。

MATFVE

フォーマット1（TYPE=FORMULA）

フォーマット2（TYPE=TABLE）

フォーマット3（TYPE=PRONY）

フォーマット4（TYPE=PRELOAD）

例

定義

コメント

フォーマット1（`TYPE`=FORMULA）

フォーマット2（`TYPE`=TABLE）

フォーマット3（`TYPE`=PRONY）

フォーマット4（`TYPE`=PRELOAD）