*MAT_156 (MUSCLE)

LS-DYNA入力インターフェースキーワード安全性用途向けの筋肉モデルを定義します。トラス要素でのみ使用されます。

フォーマット

フィールド	内容	SI単位の例
`mat_ID`	材料識別子（整数）
$ρ_{i}$	初期密度（実数）	$[\frac{kg}{m^{3}}]$
`SR_MAX`	最大ひずみ速度（実数）	$[\frac{1}{s}]$
`STS_MAX`	最大応力	$[Pa]$
`STR`	最大応力時のひずみ。（実数）
`CER`	`FUNCT_4` = 0での解析的指数式のひずみスケールファクター。（実数）
`DAMP`	減衰係数。（実数）	$[Pa \cdot s]$
`FUNCT_1`	時間に対するアクティブ化応力係数。 < 0 時間に対するアクティブ化レベルの曲線識別子。 > 0 アクティブ化レベルとして使用される定数値。（実数）
`FUNCT_2`	ひずみの関数としてのアクティブ応力係数。 < 0 曲線の識別子 > 0 定数値1.0が使用されます。（実数）
`FUNCT_3`	ひずみ速度の関数としてのアクティブ応力係数。 < 0 曲線の識別子 > 0 定数値1.0が使用されます。（整数）
`FUNCT_4`	並列弾性要素のひずみの関数としての応力係数。 < 0 曲線の識別子 = 0 解析方程式が使用されます。 4 > 0 定数値1.0が使用されます。（整数）

このキーワードは/PROP/TYPE46 (SPR_MUSCLE)にマップされます。
このキーワードは、トラス要素でのみ使用されます。トラス面積の入力は必須です。
要素内の全応力は次のように計算されます：(1)
$σ = σ_{1} + σ_{2} + σ_{3}$

ここで、

$σ_{1} = S T S_M A X \cdot F U N C T_1 (t) \cdot F U N C T_2 (Δ l) \cdot F U N C T_3 (\dot{\bar{ε}})$

収縮応力

$σ_{2} = S T S_M A X \cdot F U N C T_4 (Δ l)$

受動応力

$σ_{3} = D M P \cdot Δ l \cdot \dot{ε}$

減衰応力

ここで、 $Δ l = \frac{l_{c u r r e n t}}{l_{o r i g i n}}$ 、 $ε = Δ l - 1$ 、 $\dot{ε} = \frac{Δ ε}{Δ t}$ 、 $\dot{\bar{ε}} = \frac{Δ l \cdot \dot{ε}}{S R_M A X}$
FUNCT_4 = 0の場合、受動応力係数は次のように定義されます：

次の場合； $Δ l < 1$

$F U N C T_4 (Δ l) = 0$

右記の場合； $Δ l \geq 1$ で $C E R \neq 0$

$F U N C T_4 (Δ l) = \frac{1}{\exp (C E R) - 1} [\exp (\frac{C E R \cdot ε}{S T R}) - 1]$

右記の場合； $Δ l \geq 1$ で $C E R = 0$

$F U N C T_4 (Δ l) = \frac{ε}{S T R}$
このキーワードの最後にオプション“_TITLE”を追加することができます。“_TITLE”が含まれている場合、キーワード入力行の後に余分に80文字の長さの行が追加され、エンティティタイトルを定義できるようになります。