/MAT/LAW52 (GURSON)

ブロックフォーマットキーワード この材料則はGurson構成則に基づいており、粘弾塑性ひずみ速度依存多孔質金属をモデル化する場合に使用します。

フォーマット

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
/MAT/LAW52/mat_ID/unit_IDまたは/MAT/GURSON/mat_ID/unit_ID
mat_title
ρ i                
E ν12 Iflag Fsmooth Fcut Iyield  
A B N c p
q 1 q 2 q 3 SN εN
fI fN fc fF    
Iyield > 0の場合:
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
Tab_ID   XFAC YFAC        

定義

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

(整数、最大10桁)

 
単位識別子 単位識別子

(整数、最大10桁)

 
mat_title 材料のタイトル

(文字、最大100文字)

 
ρ i 初期密度

(実数)

[ kg m 3 ]
E ヤング率

(実数)

[ Pa ]
ν 12 ポアソン比

(実数)

 
Iflag 粘塑性流れフラグ 1
= 0
フォンミーゼス基準
= 1
フォンミーゼス基準
= 2
1 + 圧縮におけるボイド発生をゼロに設定
= 3
0 + 圧縮におけるボイド発生をゼロに設定

(整数)

 
Fsmooth ひずみ速度スムージングの計算
= 0(デフォルト)
ひずみ速度を平滑化しません。
= 1
ひずみ速度スムージングはアクティブ

(整数)

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

デフォルト = 1030(実数)

[Hz]
Iyield 降伏応力を計算するためのフラグ 3
= 0
Cowper-Symonds則を使用
= 1
降伏応力テーブル(降伏応力vs塑性ひずみ)を使用

(整数)

 
A 降伏応力。

(実数)

[ Pa ]
B 硬化パラメータ

(実数)

[ Pa ]
N 硬化指数

(実数)

 
c Cowper-Symond則のひずみ速度係数

(実数)

[ 1 s ]
p Cowper-Symond則のひずみ速度指数

(実数)

 
q 1 q 2 q 3 損傷材料のパラメータ

(実数)

 
SN Gauss標準偏差

(実数)

[ Pa ]
ε N ボイド発生後有効塑性ひずみ

(実数)

 
fI 初期ボイド体積率 2

(実数)

 
fN ボイド発生後体積率

(実数)

 
fc 凝集時臨界ボイド体積率 2

(実数)

 
fF 延性破壊時臨界ボイド体積率 2

(実数)

 
Tab_ID 降伏応力のテーブル識別子(ひずみ速度に対応する応力-ひずみ関数)

(整数)

 
XFAC Tab_IDに使用する関数の最初のエントリ(塑性ひずみ)のスケールファクター

デフォルト = 1.0(実数)

 
YFAC Tab_IDに使用する関数の縦軸(降伏応力)のスケールファクター

デフォルト = 1.0(実数)

 

例(パラメータ入力)

#RADIOSS STARTER
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/UNIT/1
unit for mat
                   g                  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/LAW52/1/1
Steel
#              RHO_I
               .0078                   
#                  E               NU_12     Iflag   Fsmooth                Fcut    Iyield
              200000                  .3         0         0                   0         0
#                  A                   B                   N                   c                   p
                 200                 533                   1                 802               3.585
#                q_1                 q_2                 q_3                  SN                EpsN
                1.25                   1                2.25                  .1                  .2
#                 Fi                  FN                  Fc                  FF
                 .01                 .04                 .12                  .2
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#ENDDATA
/END
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

例(関数入力)

#RADIOSS STARTER
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/UNIT/1
unit for mat
                   g                  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/LAW52/1/1
Steel
#              RHO_I
               .0078                   
#                  E               NU_12     Iflag   Fsmooth                Fcut    Iyield
              200000                  .3         0         0                   0         1
#                  A                   B                   N                   c                   p
                 200                 533                   1                 802               3.585
#                q_1                 q_2                 q_3                  SN                EpsN
                1.25                   1                2.25                  .1                  .2
#                 Fi                  FN                  Fc                  FF
                 .01                 .04                 .12                  .2
#   Tab_ID                XFAC                YFAC
      1000                   0                   0
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/TABLE/1/1000
curve_list with strain rates
         2
# function                    stain rate
     10010                        1.0e-4
     10020                           1.0
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/FUNCT/10010
plastic strain vs yield stress funct dt=1.0e-4
#     plastic strain       yield stress    
              0.0000                200.
              1.0000                733.
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/FUNCT/10020
plastic strain vs yield stress funct dt=1.0
#     plastic strain       yield stress
              0.0000                250.
              1.0000                783.
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#ENDDATA
/END
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

コメント

  1. 粘塑性流れのフォンミーゼス基準:
    Iflag = 0の場合:(1)
    Ω v m = σ e q σ M 1 + q 3 f * 2 2 q 1 f * 2 cosh ( 3 q 2 σ m 2 σ M )
    Iflag = 1の場合:(2)
    Ω v m = σ e q 2 σ M 2 + 2 q 1 f * cosh ( 3 2 q 2 σ m σ M ) ( 1 + q 3 f 2 )

    の、 σ m > 0

    (3)
    Ω v m = σ e q 2 σ M 2 + 2 q 1 f ( 1 + q 3 f 2 )

    の、 σ m 0

    ここで、 σ M は可容応力、 σ m はtrace[ σ ](静水圧応力)、 σ e q はフォンミーゼス応力、 q 1 q 2 および q 3 はGurson材料則の材料パラメータ、

    q 3 = q 1 2

    f * は固有の凝集関数です。

    f * = f の、 f f c

    f * = f c + f u f c f F f c ( f f c ) の、 f > f c

    ここで、

    f u = 1 q 1 は、右記の凝集関数に相当; f u = f * ( f F )

  2. ボイド体積率のパラメータの入力は必須で、 f I < f c < f F MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOzamaaBa aaleaacaWGjbaabeaakiabgYda8iaadAgadaWgaaWcbaGaam4yaaqa baGccqGH8aapcaWGMbWaaSbaaSqaaiaadAeaaeqaaaaa@3DD8@ となる。
  3. いずれかの積分点が f * f F に達した場合、その要素は削除されます。
  4. Iyieldフラグがアクティブでない場合(Iyield=0)、降伏応力はCowper-Symond材料則を使って計算されます。
    (4)
    σ M = ( A + B ε M N ) ( 1 + ( ε ˙ c ) 1 p )

    Iyieldフラグがアクティブである場合(Iyield=1)、降伏応力は、降伏応力曲線(Tab_ID)から直接計算されます。

  5. この材料則は、シェル要素およびソリッド要素に対して使用可能です。
  6. プロットファイル(/TH/SHEL/TH/SH3N、および/TH/BRICK)、またはアニメーションファイル(/ANIM)では下記の変数を使用できます。
    • USR1: 塑性ひずみ ε M
    • USR2: f *
    • USR3: 可容応力 σ M
    • USR4: f
    • USR5: ε