/MAT/LAW35 (FOAM_VISC)

ブロックフォーマットキーワード この材料則は、汎用Maxwell-Kelvin-Voigtモデルを使用して粘弾性フォーム材料を記述します。このモデルの粘性はNavier方程式に従います。

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

フォーマット

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
/MAT/LAW35/mat_ID/unit_IDまたは/MAT/FOAM_VISC/mat_ID/unit_ID
mat_title
ρ i                
E ν E1 E2 n
C1 C2 C3   IFlag Pmin
fct_IDf   Fscaleprs       Fsmooth Fcut
Et νt η0 λ    
P0 Φ γ0        

定義

フィールド 内容 SI単位の例
mat_ID 材料識別子

(整数、最大10桁)

 
unit_ID 単位識別子

(整数、最大10桁)

 
mat_title 材料のタイトル

(文字、最大100文字)

 
ρ i 初期密度

(実数)

[ kg m 3 ]
E ヤング率

(実数)

[ Pa ]
ν ポアソン比

(実数)

 
E1 ヤング率更新の係数 E = E 1 ε ˙ + E 2

(実数)

[ Pa ]
E2 ヤング率更新の係数

(実数)

[ Pa ]
n 相対体積の指数

(実数)

 
C1 圧力計算の係数

(実数)

 
C2 圧力計算の係数

(実数)

 
C3 圧力計算の係数

(実数)

 
IFlag 開放セルフォームフラグ
= 0
入力は圧力対圧縮曲線/MAT/LAW21(DPRAG)の場合と同じ。
= 1
入力は、システム対圧力から削除される“相当気圧”を定義する関数。これは開放セルフォーム定式化と対応します。

(整数)

 
Pmin 最小圧力

(実数)

[ Pa ]
fct_IDf 体積ひずみの関数としての圧力の曲線識別子体積ひずみは γ = ρ 0 ρ 1

(整数)

 
Fscaleprs 圧力関数スケールファクター

デフォルト = 1.0(実数)

 
Fsmooth ひずみ速度スムージングオプションフラグ。
= 0(デフォルト)
ひずみ速度を平滑化しません。
= 1
ひずみ速度スムージングはアクティブ。
(整数)
 
Fcut ひずみ速度フィルタリングのカットオフ周波数。

デフォルト = 10000 Hz (実数)

[Hz]
Et 接線率

(実数)

[ Pa ]
ν t 接線ポアソン比

(実数)

 
η 0 純せん断の粘性係数(Navierの定数)

(実数)

[ Pas ]
λ Navierの定数

(実数)

[ Pas ]
P0 初期気圧

(実数)

[ 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. いずれの状況でも、せん断および体積弾性率の計算ではヤング率の下記の値が使用されます:(1)
    E = max ( E , E 1 ε ˙ + E 2 ) . ( V V 0 ) n
  2. fct_IDf = 0の場合、(2)
    d P d t = C 1 K ε ˙ k k C 2 [ K + K t 3 λ + 2 η 0 σ k k ] + C 3 [ K K t 3 λ + 2 η 0 ε k k ]
    ここで、
    K = E 3 ( 1 2 v )
    K t = E t 3 ( 1 2 v t )
    P = 1 3 σ k k
    ε k k = ln ( V V 0 )
  3. fct_IDf0の場合、圧力は曲線から読み取られます。
  4. 密閉セルポリウレタンフォームでは、下記の分だけ骨格球面応力が増加する場合があります:(3)
    P a i r = P 0 γ 1 + γ Φ
  5. パラメータFsmoothFcutは、ひずみ速度のフィルタリングを可能にします。3つのケースが設定できます:
    • Fsmooth = 0 および Fcut = 0.0の場合、ひずみ速度フィルタリングはオフ。
    • Fsmooth = 1 および Fcut = 0.0の場合、ひずみ速度フィルタリングは、デフォルトのカットオフ周波数10 kHzを使用。
    • Fcut ≠ 0, Fsmoothの場合、自動的に1に設定され、ひずみ速度フィルタリングは、指定されたカットオフ周波数を使用。