/MAT/LAW15 (CHANG)

ブロックフォーマットのキーワード この材料則は、LAW25と同様に、複合シェル要素のモデル化に使用されます。Tsai-Wu基準(Tsai-Wuの説明については/MAT/LAW25 (COMPSH)をご参照ください)に基づく塑性挙動、およびChang-Chang破壊基準に基づく破壊が使用されます。

ただし、LAW15ではなく、LAW25を個別のChang-Chang破壊基準(/MAT/LAW25でキーワード/FAIL/CHANGを指定)と組み合わせて使用することが推奨されます。

フォーマット

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
/MAT/LAW15/mat_ID/unit_IDまたは/MAT/CHANG/mat_ID/unit_ID
mat_title
ρ i                
E11 E22 ν12        
G12 G23 G31        
b n fmax        
Wpmax Wpref Ioff          
σ 1 y t σ 2 y t σ 1 y c σ 2 y c α
σ 12 y c σ 12 y t c ε ˙ 0 ICC  
β τ max MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqiXdq3aaS baaSqaaiGac2gacaGGHbGaaiiEaaqabaaaaa@3ABC@ S1 S2 S12
Fsmooth Fcut C1 C2      

定義

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

(整数、最大10桁)

 
unit_ID 単位識別子

(整数、最大10桁)

 
mat_title 材料のタイトル

(文字、最大100文字)

 
ρ i 初期密度

(実数)

[ kg m 3 ]
E11 方向1のヤング率

(実数)

[ Pa ]
E22 方向2のヤング率

(実数)

[ Pa ]
ν 12 ポアソン比

(実数)

 
G12 せん断係数

(実数)

[ Pa ]
G23 せん断係数

(実数)

[ Pa ]
G311 せん断係数

(実数)

[ Pa ]
b 硬化パラメータ

(実数)

 
n 硬化指数

デフォルト = 1.0(実数)

 
fmax 降伏関数の最大値 2

デフォルト = 1030(実数)

[ Pa ]
W p max 体積単位あたりの最大塑性エネルギー

デフォルト = 1030(実数)

[ J m 3 ]
W p ref 体積単位あたりの基準塑性エネルギー

デフォルト = 1.0(実数)

[ J m 3 ]
Ioff 全要素破壊規準 4
= 0
1つの層で W p * > W p max * の場合にシェルを削除
= 1
すべての層で W p * > W p max * の場合にシェルを削除
= 2
各層で、 W p * > W p max * または方向1の引張破壊の場合
= 3
各層で、 W p * > W p max * または方向2の引張破壊の場合
= 4
各層で、 W p * > W p max * または方向1および2の引張破壊の場合
= 5
すべてのレイヤーの場合: W p * > W p max * または、方向1の引張破壊
または、すべてのレイヤーの場合: W p * > W p max * または、方向2の引張損傷
= 6
各層で、 W p * > W p max * あるいは方向1または方向2の引張破壊の場合

(整数)

 
σ 1 y t 方向1の張力における複合材降伏応力 2

(実数)

[ Pa ]
σ 2 y t 方向2の張力における複合材降伏応力

(実数)

[ Pa ]
σ 1 y c 方向1の圧縮における複合材降伏応力

(実数)

[ Pa ]
σ 2 y c 方向2の圧縮における複合材降伏応力

(実数)

[ Pa ]
α F12減少係数 2

デフォルトは1.0に設定されます(実数)

 
σ 12 y c 方向12のせん断およびひずみ速度圧縮における降伏応力

(実数)

[ Pa ]
σ 12 y t 方向12のせん断およびひずみ速度張力における降伏応力

(実数)

[ Pa ]
c せん断およびひずみ速度の降伏応力係数 2
= 0
ひずみ速度依存なし。

(実数)

 
ε ˙ 0 せん断およびひずみ速度参照の降伏応力

(実数)

[ 1 s ]
ICC ひずみ速度計算フラグ 2
= 1(デフォルト)
fmax にはひずみ速度効果あり、 W p max にはなし
= 2
右記に対するひずみ速度効果なし; fmaxおよび W p max
= 3
右記にに対するひずみ速度効果あり; fmaxおよび W p max
= 4
fmaxにはひずみ速度効果なし、 W p max にはあり

(整数)

 
β せん断スケーリング係数 1

(実数)

 
τ max MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqiXdq3aaS baaSqaaiGac2gacaGGHbGaaiiEaaqabaaaaa@3ABC@ 時間緩和 3

デフォルト = 1030(実数)

[ s ]
S1 縦方向引張り強度 1

デフォルト = 1030(実数)

[ Pa ]
S2 横方向引張り強度

デフォルト = 1030(実数)

[ Pa ]
S12 せん断強度

デフォルト = 1030(実数)

[ Pa ]
Fsmooth ひずみ速度スムージングオプションフラグ
= 0(デフォルト)
ひずみ速度を平滑化しません
= 1
ひずみ速度スムージングはアクティブ

(整数)

 
Fcut ひずみ速度フィルタリングのカットオフ周波数。

デフォルト = 1030(実数)

[Hz]
C1 縦方向圧縮強度 1

デフォルト = 1030(実数)

[ Pa ]
C2 横方向圧縮強度

デフォルト = 1030(実数)

[ Pa ]

例(Carbon)

#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/LAW15/1/1
Carbon
#              RHO_I
              1.8E-6                   0
#                E11                 E22                nu12
                  41                 3.3                  .3
#                G12                 G23                 G31
                 5.2                 1.3                 1.3
#                  b                   n                fmax
                8E-6                   1              100000
#              Wpmax              Wpref       Ioff
              100000                   0         0
#          sigma_1yt           sigma_2yt           sigma_1yc          sigma_2yc               alpha
                .786               .1566                .786               .1566                   0
#         sigma_12yc          sigma_12yt                   c           Eps_dot_0       ICC
               .0655               .0655                   0                   0         0
#               beta                Tmax                  S1                  S2                 S12
                   1                 .01                   0                   0                   0
#  Fsmooth                Fcut                  C1                 C12
         0                   0                   0                   0
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#ENDDATA
/END
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

コメント

  1. Chang Chang破壊基準
    Chang Chang破壊基準では、2つの異なる破断挙動を表すために6つの材料パラメータが使用されます。
    • 繊維破損の場合、破壊基準は:
      • 引張繊維モード σ 11 > 0 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeq4Wdm3aaS baaSqaaiaaigdacaaIXaaabeaakiabg6da+iaaicdaaaa@3B28@
        (1)
        e f 2 = ( σ 11 S 1 ) 2 + β ( σ 12 S 12 ) 2 1.0
        0 failed < 0 elastic plastic
      • 圧縮繊維モード σ 11 < 0 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeq4Wdm3aaS baaSqaaiaaigdacaaIXaaabeaakiabgYda8iaaicdaaaa@3B24@
        (2)
        e c 2 = ( σ 11 C 1 ) 2 1.0
        0 failed < 0 elastic plastic
    • マトリックス亀裂の破壊基準は:
      • 引張繊維モード σ 22 > 0 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeq4Wdm3aaS baaSqaaiaaikdacaaIYaaabeaakiabg6da+iaaicdaaaa@3B2A@
        (3)
        e m 2 = ( σ 22 S 2 ) 2 + β ( σ 12 S 12 ) 2 1.0
        0 failed < 0 elastic plastic
      • 圧縮マトリックスモード σ 22 < 0 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeq4Wdm3aaS baaSqaaiaaikdacaaIYaaabeaakiabgYda8iaaicdaaaa@3B26@
      (4)
      e d 2 = ( σ 22 2 S 12 ) 2 + [ ( C 2 2 S 12 ) 2 1 ] σ 22 C 2 + ( σ 12 S 12 ) 2 1.0
      0 failed < 0 elastic plastic
  2. 破壊される以前(損傷パラメータ e f 2 , e c 2 , e m 2 , e d 2 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyzamaaBa aaleaacaWGMbaabeaakmaaCaaaleqabaGaaGOmaaaakiaacYcacaWG LbWaaSbaaSqaaiaadogaaeqaaOWaaWbaaSqabeaacaaIYaaaaOGaai ilaiaadwgadaWgaaWcbaGaamyBaaqabaGcdaahaaWcbeqaaiaaikda aaGccaGGSaGaamyzamaaBaaaleaacaWGKbaabeaakmaaCaaaleqaba GaaGOmaaaaaaa@43F6@ が0未満)、材料は弾塑性相にあります。塑性挙動はTsai-Wu基準に基づきます(Tsai-Wu基準の説明については、Tsai-Wu定式化(Iform =0)をご参照ください)。
  3. 破壊された後(損傷パラメータ e f 2 , e c 2 , e m 2 , e d 2 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyzamaaBa aaleaacaWGMbaabeaakmaaCaaaleqabaGaaGOmaaaakiaacYcacaWG LbWaaSbaaSqaaiaadogaaeqaaOWaaWbaaSqabeaacaaIYaaaaOGaai ilaiaadwgadaWgaaWcbaGaamyBaaqabaGcdaahaaWcbeqaaiaaikda aaGccaGGSaGaamyzamaaBaaaleaacaWGKbaabeaakmaaCaaaleqaba GaaGOmaaaaaaa@43F6@ が0以上)、数値の不安定性を回避するため、指数関数を使用して応力を減少します。
    緩和のテクニックが応力を徐々に減少させるために使用されます。(5)
    σ ( t ) = f ( t ) σ d ( t r ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaC4WdiaacI cacaWG0bGaaiykaiabg2da9iGacAgaciGGOaGaamiDaiaacMcacqGH flY1caWHdpWaaSbaaSqaaiaadsgaaeqaaOGaaiikaiaadshadaWgaa WcbaGaamOCaaqabaGccaGGPaaaaa@4615@

    緩和の関数の使用で:

    f ( t ) = exp ( t t r τ max ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaciOzaiGacI cacaWG0bGaaiykaiabg2da9iGacwgacaGG4bGaaiiCamaabmaabaGa eyOeI0YaaSaaaeaacaWG0bGaeyOeI0IaamiDamaaBaaaleaacaWGYb aabeaaaOqaaiabes8a0naaBaaaleaaciGGTbGaaiyyaiaacIhaaeqa aaaaaOGaayjkaiaawMcaaaaa@4878@ および t t r MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWG0bGaey yzImRaamiDamaaBaaaleaacaWGYbaabeaaaaa@3C43@

    ここで、
    t MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiDaaaa@36EF@
    時間
    t r MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiDamaaBa aaleaacaWGYbaabeaaaaa@3812@
    損傷基準が推定される場合における緩和の開始時間
    τ max MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqiXdq3aaS baaSqaaiGac2gacaGGHbGaaiiEaaqabaaaaa@3ABB@
    動的緩和の時間
    σ d ( t r ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaC4WdmaaBa aaleaacaWGKbaabeaakiaacIcacaWG0bWaaSbaaSqaaiaadkhaaeqa aOGaaiykaaaa@3BE4@
    損傷開始時の応力成分
  4. シェルに複数の層が存在し、材料が層別である場合(各層で材料が異なり、Ioffも異なる)、使用されるIoffは、シェル要素定義でシェルに関連付けられたものになります。