/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 | |||||||||
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文字) |
|
初期密度 (実数) |
||
C10 | 超弾性モデル用の材料パラメータ デフォルト = 0.0(実数) |
|
C01 | 超弾性モデル用の材料パラメータ デフォルト = 0.0(実数) |
|
C20 | 超弾性モデル用の材料パラメータ デフォルト = 0.0(実数) |
|
C11 | 超弾性モデル用の材料パラメータ デフォルト = 0.0(実数) |
|
C02 | 超弾性モデル用の材料パラメータ デフォルト = 0.0(実数) |
|
C30 | 超弾性モデル用の材料パラメータ デフォルト = 0.0(実数) |
|
C21 | 超弾性モデル用の材料パラメータ デフォルト = 0.0(実数) |
|
C12 | 超弾性モデル用の材料パラメータ デフォルト = 0.0(実数) |
|
C03 | 超弾性モデル用の材料パラメータ デフォルト = 0.0(実数) |
|
Sb | ネットワークB用の応力スケールファクター デフォルト = 0.0(実数) |
|
D1 | 体積弾性率の計算に使用する体積材料パラメータ1
デフォルト = 0.0(実数) |
|
D2 | 体積材料パラメータ2 デフォルト = 0.0(実数) |
|
D3 | 体積材料パラメータ3 デフォルト = 0.0(実数) |
|
A | 有効クリープひずみ速度 デフォルト = 0.0(正の実数) |
|
C | ネットワークBでの有効クリープひずみ速度のクリープひずみ依存性を特性化する指数(-1 < C < 0) デフォルト = -0.7(実数) |
|
M | ネットワークBでの有効クリープひずみ速度の有効応力依存性を特性化する正の指数(
) デフォルト = 1.0(実数) |
|
非変形状態近傍のクリープひずみ速度の正則化の定数 デフォルト = 0.01(実数) |
||
Tau_ref | セカンダリネットワークでの有効クリープひずみ速度の参照応力 デフォルト = 1.0(実数) |
例
#RADIOSS STARTER
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/UNIT/1
unit for mat
kg mm ms
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#- 2. MATERIALS:
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/MAT/LAW95/1/1
BERGSTROM
# RHO_I
1.42E-6
# C10 C01 C20 C11 C22
0.2019 0. 4.43E-5
# C30 C21 C12 C03 Sb
1.295E-4 0. 0. 0. 2.0
# D1 D2 D3
2.1839E-3 8.68E-5 -1.794E-5
# A EXPC EXPM KSI Tau_ref
1.0E-1 -0.7 5 0.01
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#ENDDATA
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
コメント
- 材料の応答は、2つの並列ネットワークAとBを用いて表すことができます。ネットワークAは、非線形超弾性コンポーネントを伴う均衡ネットワークです。ネットワークBでは、非線形超弾性コンポーネントは非線形粘弾性流れ要素と直列であり、したがって、時間依存のネットワークです。
- 両方のネットワークで、超弾性コンポーネントに同じ多項式ひずみエネルギーポテンシャルが用いられます。ネットワークBでは、このポテンシャルは係数Sbによってスケーリングされます。ひずみエネルギー密度はその後、ネットワークの超弾性コンポーネント用に書き出されます:
(1) および(2) ここで、 -
の特別な値について、多項式モデルは下記の材料モデルに縮小することが可能です:
- Yeoh: j=0
ここで、C10,C20、C30は0以外
- Mooney-Rivlin: i+j =1
ここで、C10とC01は0以外で、D2 =D3=0
- Neo-Hookean:
C10とD1のみが0以外
- Yeoh: j=0
- 初期せん断弾性係数と体積弾性係数は次のように計算されます:
(3) および(4) - D1 = 0の場合、非圧縮性材料が考慮されます。
- =0の場合、超弾性多項式材料モデルのみが粘弾性時間依存材料の応答なしで使用されます。
- ネットワーク
での有効クリープひずみ速度は、次式で与えられます:
(5) ここで、- および
- ネットワークBでの有効応力
- 、 、 および
- 入力材料パラメータ