Loading [MathJax]/jax/output/CommonHTML/jax.js

/MAT/LAW95 (BERGSTROM_BOYCE)

ブロックフォーマットキーワード この材料則は、エラストマー状の材料の非線形時間依存を予測するための構成モデルです。

非線形粘弾性時間依存材料の応答を表すために、超弾性材料応答およびBergstrom-Boyce材料モデル1に多項式材料モデルを使用します。この材料則はソリッド要素とのみ適合性があります。

フォーマット

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
/MAT/LAW95/mat_ID/unit_IDまたは/MAT/BERGSTROM_BOYCE/mat_ID/unit_ID
mat_title
ρi                
C10 C01 C20 C11 C02
C30 C21 C12 C03 sb
D1 D2 D3    
A C M ξ Tau_ref

定義

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

(整数、最大10桁)

 
unit_ID 単位識別子

(整数、最大10桁)

 
mat_title 材料のタイトル

(文字、最大100文字)

 
ρi 初期密度

(実数)

[kgm3]
C10 超弾性モデル用の材料パラメータ

デフォルト = 0.0(実数)

[Pa]
C01 超弾性モデル用の材料パラメータ

デフォルト = 0.0(実数)

[Pa]
C20 超弾性モデル用の材料パラメータ

デフォルト = 0.0(実数)

[Pa]
C11 超弾性モデル用の材料パラメータ

デフォルト = 0.0(実数)

[Pa]
C02 超弾性モデル用の材料パラメータ

デフォルト = 0.0(実数)

[Pa]
C30 超弾性モデル用の材料パラメータ

デフォルト = 0.0(実数)

[Pa]
C21 超弾性モデル用の材料パラメータ

デフォルト = 0.0(実数)

[Pa]
C12 超弾性モデル用の材料パラメータ

デフォルト = 0.0(実数)

[Pa]
C03 超弾性モデル用の材料パラメータ

デフォルト = 0.0(実数)

[Pa]
Sb ネットワークB用の応力スケールファクター

デフォルト = 0.0(実数)

 
D1 体積弾性率の計算に使用する体積材料パラメータ1

K=2D1

デフォルト = 0.0(実数)

[1Pa]
D2 体積材料パラメータ2

デフォルト = 0.0(実数)

[1Pa]
D3 体積材料パラメータ3

デフォルト = 0.0(実数)

[1Pa]
A 有効クリープひずみ速度

デフォルト = 0.0(正の実数)

[1s]
C ネットワークBでの有効クリープひずみ速度のクリープひずみ依存性を特性化する指数(-1 < C < 0)

デフォルト = -0.7(実数)

 
M ネットワークBでの有効クリープひずみ速度の有効応力依存性を特性化する正の指数( M1.0

デフォルト = 1.0(実数)

 
ξ 非変形状態近傍のクリープひずみ速度の正則化の定数

デフォルト = 0.01(実数)

 
Tau_ref 2番目のネットワークでの有効クリープひずみ速度の基準応力

デフォルト = 1.0(実数)

[Pa]

コメント

  1. 材料の応答は、2つの並列ネットワークAとBを用いて表すことができます。ネットワークAは、非線形超弾性コンポーネントを伴う均衡ネットワークです。ネットワークBでは、非線形超弾性コンポーネントは非線形粘弾性流れ要素と直列であり、したがって、時間依存のネットワークです。


    図 1. 時間依存のネットワーク
  2. 両方のネットワークで、超弾性コンポーネントに同じ多項式ひずみエネルギーポテンシャルが用いられます。ネットワークBでは、このポテンシャルは係数Sbによってスケーリングされます。ひずみエネルギー密度はその後、ネットワークの超弾性コンポーネント用に書き出されます:(1)
    WA=3i+j=1Cij(ˉI13)i(ˉI23)j+3i=11Di(J1)2i
    および(2)
    WB=SbWA
    ここで、
    • ˉI1=ˉλ21+ˉλ22+ˉλ23
    • ˉI2=ˉλ21+ˉλ22+ˉλ23
    • ˉλi=J13λi
  3. Cij の特別な値について、多項式モデルは下記の材料モデルに縮小することが可能です:
    • Yeoh: j=0

      ここで、C10,C20C30は0以外

    • Mooney-Rivlin: i+j =1

      ここで、C10C01は0以外で、D2 =D3=0

    • Neo-Hookean:

      C10D1のみが0以外

  4. 初期せん断弾性係数と体積弾性係数は次のように計算されます:(3)
    μ=2(Sb+1)(C10+C01)
    および(4)
    K=2D1(1+Sb)
  5. D1 = 0の場合、非圧縮性材料が考慮されます。
  6. A =0の場合、超弾性多項式材料モデルのみが粘弾性時間依存材料の応答なしで使用されます。
  7. ネットワーク B での有効クリープひずみ速度は、次の式で与えられます:(5)
    ˙εvB=A(˜λ1+ξ)C(ˉσBτref)M
    ここで、
    ˜λ=ˉI13 および ˉσB
    ネットワークBでの有効応力
    ξ M C および τref
    入力材料パラメータ
1
Bergström, J. S., and M. C. Boyce.Constitutive modeling of the large strain time-dependent behavior of elastomers."Journal of the Mechanics and Physics of Solids" 46, No. 5 (1998): 931-954