/MAT/LAW35 (FOAM

この材料則は、シェルおよびソリッド要素のみに適用され、開放セルのフォーム、ポリマーおよびエラストマー、シートクッション、ダミーの詰め物に使用できます。

フォーマット

(1)	(3)	(5)	(7)	(8)	(9)
/MAT/LAW35/`mat_ID`/`unit_ID`または/MAT/FOAM_VISC/`mat_ID`/`unit_ID`
`mat_title`
$ρ_{i}$
`E`	$ν$	`E`₁	`E`₂		`n`
`C`₁	`C`₂	`C`₃		`I`_Flag	`P`_min
`fct_ID`_f	`Fscale`_prs			`F`_smooth	`F`_cut
`E`_t	$ν_{t}$	$η_{0}$	$λ$
`P`₀	$Φ$	$γ_{0}$

定義

フィールド	内容	SI単位の例
`mat_ID`	材料識別子（整数、最大10桁）
`unit_ID`	単位識別子（整数、最大10桁）
`mat_title`	材料のタイトル（文字、最大100文字）
$ρ_{i}$	初期密度（実数）	$[\frac{kg}{m^{3}}]$
`E`	ヤング率（実数）	$[Pa]$
$ν$	ポアソン比（実数）
`E`₁	ヤング率更新の係数 $E = E_{1} \dot{ε} + E_{2}$ （実数）	$[Pa]$
`E`₂	ヤング率更新の係数（実数）	$[Pa]$
`n`	相対体積の指数（実数）
`C`₁	圧力計算の係数（実数）
`C`₂	圧力計算の係数（実数）
`C`₃	圧力計算の係数（実数）
`I`_Flag	開放セルフォームフラグ = 0 入力は圧力対圧縮曲線/MAT/LAW21（DPRAG）の場合と同じ。 = 1 入力は、システム対圧力から削除される“相当気圧”を定義する関数。これは開放セルフォーム定式化と対応します。（整数）
`P`_min	最小圧力（実数）	$[Pa]$
`fct_ID`_f	体積ひずみの関数としての圧力の曲線識別子体積ひずみは $γ = \frac{ρ_{0}}{ρ} - 1$ （整数）
`Fscale`_prs	圧力関数スケールファクターデフォルト = 1.0（実数）
`F`_smooth	ひずみ速度スムージングオプションフラグ。 = 0（デフォルト）ひずみ速度を平滑化しません。 = 1 ひずみ速度スムージングはアクティブ。（整数）
`F`_cut	ひずみ速度フィルタリングのカットオフ周波数。デフォルト = 10000 Hz （実数）	$[Hz]$
`E`_t	接線率（実数）	$[Pa]$
$ν_{t}$	接線ポアソン比（実数）
$η_{0}$	純せん断の粘性係数（Navierの定数）（実数）	$[Pa \cdot s]$
$λ$	Navierの定数（実数）	$[Pa \cdot s]$
`P`₀	初期気圧（実数）	$[Pa]$
$Φ$	発泡とポリマーの密度の比率。（実数）
$γ_{0}$	初期体積ひずみ。（実数）

例（フォーム材）

#RADIOSS STARTER
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/UNIT/1
unit for mat
                  Mg                  mm                   s
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#-  2. MATERIALS:
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/MAT/FOAM_VISC/1/1
Foam
#              RHO_I
             4.3E-11                   
#                  E                  Nu                  E1                  E2                   n
                  .2                   0                   0                   0                   0
#                 C1                  C2                  C3               Iflag                Pmin
                   1                   1                   1                   0                   0
# func_IDf                    Fscale_prs                                 Fsmooth                Fcut
      1076                             0                                       1             15000.0
#                 Et                Nu_t               eta_0               Lamda
                 .25                   0               10000                   0
#                 P0                 Phi             gamma_0
                   0                   0                   0
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#-  3. FUNCTIONS:
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/FUNCT/1076
Pressure function
#                  X                   Y
#  Volumetric strain            Pressure
             -100000               -1000                                                            
                 -10               -1000                                                            
                   0                   0                                                            
                3000               7.633                                                            
              209000               7.633                                                            
              210000                18.5                                                            
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#ENDDATA
/END
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

いずれの状況でも、せん断および体積弾性率の計算ではヤング率の下記の値が使用されます：(1)
$E = \max (E, E_{1} \dot{ε} + E_{2}) . {(\frac{V}{V_{0}})}^{n}$
fct_ID_f = 0の場合、(2)
$\frac{d P}{d t} = C_{1} K {\dot{ε}}_{k k} - C_{2} [\frac{K + K_{t}}{3 λ + 2 η_{0}} σ_{k k}] + C_{3} [\frac{K K_{_{t}}}{3 λ + 2 η_{0}} ε_{k k}]$

ここで、

$K = \frac{E}{3 (1 - 2 v)}$

$K_{t} = \frac{E_{t}}{3 (1 - 2 v_{t})}$

$P = - \frac{1}{3} σ_{k k}$

$ε_{k k} = \ln (\frac{V}{V_{0}})$
fct_ID_f ≠ 0の場合、圧力は曲線から読み取られます。
密閉セルポリウレタンフォームでは、下記の分だけ骨格球面応力が増加する場合があります：(3)
$P_{a i r} = - \frac{P_{0} \cdot γ}{1 + γ - Φ}$
パラメータF_smoothとF_cutは、ひずみ速度のフィルタリングを可能にします。3つのケースが設定できます：
- F_smooth = 0 および F_cut = 0.0の場合、ひずみ速度フィルタリングはオフ。
- F_smooth = 1 および F_cut = 0.0の場合、ひずみ速度フィルタリングは、デフォルトのカットオフ周波数10 kHzを使用。
- F_cut ≠ 0, F_smoothの場合、自動的に1に設定され、ひずみ速度フィルタリングは、指定されたカットオフ周波数を使用。

/MAT/LAW35 (FOAM_VISC)

フォーマット

定義

例（フォーム材）

コメント