/MAT/LAW60 (PLAS_T3)

ブロックフォーマットキーワード この材料則は、多様なひずみ速度の応力-ひずみ曲線(塑性ひずみと応力の関係の曲線)の加工硬化部について、ユーザー定義関数を使用して等方性弾塑性材料をモデル化します。

降伏応力が関数からの非線形補間であることを除き、LAW36と同様です。

フォーマット

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
/MAT/LAW60/mat_ID/unit_IDまたは/MAT/PLAS_T3/mat_ID/unit_ID
mat_title
ρ i                
E ν ε p m a x εt εm
Nfunct Fsmooth Chard Fcut      
fct_IDp Fscale fct_IDE Einf CE    
fct_ID1 fct_ID2 fct_ID3 fct_ID4 fct_ID5          
6 ≤ Nfunct ≤ 10の場合のみ読み込み
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
fct_ID6 fct_ID7 fct_ID8 fct_ID9 fct_ID10          
常に読み込み
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
Fscale1 Fscale2 Fscale3 Fscale4 Fscale5
6 < Nfunct < 10の場合のみ読み込み
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
Fscale6 Fscale7 Fscale8 Fscale9 Fscale10
常に読み込み
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
ε ˙ 1 ε ˙ 2 ε ˙ 3 ε ˙ 4 ε ˙ 5
6 ≤ Nfunct ≤ 10の場合のみ読み込み
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
ε ˙ 6 ε ˙ 7 ε ˙ 8 ε ˙ 9 ε ˙ 10

定義

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

(整数、最大10桁)

 
unit_ID 単位識別子

(整数、最大10桁)

 
mat_title 材料のタイトル

(文字、最大100文字)

 
ρ i 初期密度

(実数)

[ kg m 3 ]
E ヤング率

(実数)

[ Pa ]
ν ポアソン比。

(実数)

 
ε p m a x 破壊塑性ひずみ。

デフォルト = 1.0 × 1030(実数)

 
ε t 応力が減少し始める引張破壊ひずみ。

デフォルト = 1.0 x 1030(実数)

 
ε m 要素が削除される最大引張破壊ひずみ

デフォルト = 2.0 x 1030(実数)

 
Nfunct 関数の数4 <Nfunct < 10であることが必要

デフォルト ≤ 10(整数)

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

(整数)

 
Chard 硬化係数。
= 0
硬化は完全等方性モデルです。
= 1
硬化はPrager-Ziegler運動モデル
= 01の値
硬化は2つのモデル間で補間されます

(実数)

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

デフォルト = 1.0 × 1030(実数)

[Hz]
fct_IDp 圧力対降伏係数の関数 9

デフォルト = 0(整数)

 
Fscale fct_IDpの降伏係数のスケールファクター

デフォルト = 1.0(実数)

[ Pa ]
fct_IDE ヤング率のスケールファクターの関数ID(ヤング率が塑性ひずみの関数である場合) 6

デフォルト = 0: この場合はヤング率の進展はEinfおよびCEに応じて決まります。

(整数)

 
Einf 無限の塑性ひずみに対する飽和ヤング率

(実数)

[ Pa ]
CE ヤング率進展のパラメータ

(実数)

 
fct_ID1 ひずみ速度 ε ˙ 1 に対応する降伏応力関数1の識別子

(整数)

 
fct_ID2 ひずみ速度 ε ˙ 2 に対応する降伏応力関数2の識別子

(整数)

 
fct_ID3 ひずみ速度 ε ˙ 3 に対応する降伏応力関数3の識別子

(整数)

 
fct_ID4 ひずみ速度 ε ˙ 4 に対応する降伏応力関数4の識別子

(整数)

 
fct_ID5 ひずみ速度 ε ˙ 5 に対応する降伏応力関数5の識別子

(整数)

 
fct_ID6 ひずみ速度 ε ˙ 6 に対応する降伏応力関数6の識別子

(整数)

 
fct_ID7 ひずみ速度 ε ˙ 7 に対応する降伏応力関数7の識別子

(整数)

 
fct_ID8 ひずみ速度 ε ˙ 8 に対応する降伏応力関数8の識別子

(整数)

 
fct_ID9 ひずみ速度 ε ˙ 9 に対応する降伏応力関数9の識別子

(整数)

 
fct_ID10 ひずみ速度 ε ˙ 10 に対応する降伏応力関数10の識別子

(整数)

 
Fscale1 fct_ID1の縦軸(応力)のスケールファクター

デフォルト = 1.0(実数)

[ Pa ]
Fscale2 fct_ID2の縦軸(応力)のスケールファクター

デフォルト = 1.0(実数)

[ Pa ]
Fscale3 fct_ID3の縦軸(応力)のスケールファクター

デフォルト = 1.0(実数)

[ Pa ]
Fscale4 fct_ID4の縦軸(応力)のスケールファクター

デフォルト = 1.0(実数)

[ Pa ]
Fscale5 fct_ID5の縦軸(応力)のスケールファクター

デフォルト = 1.0(実数)

[ Pa ]
Fscale6 fct_ID6の縦軸(応力)のスケールファクター

デフォルト = 1.0(実数)

[ Pa ]
Fscale7 fct_ID7の縦軸(応力)のスケールファクター

デフォルト = 1.0(実数)

[ Pa ]
Fscale8 fct_ID8の縦軸(応力)のスケールファクター

デフォルト = 1.0(実数)

[ Pa ]
Fscale9 fct_ID9の縦軸(応力)のスケールファクター

デフォルト = 1.0(実数)

[ Pa ]
Fscale10 fct_ID10の縦軸(応力)のスケールファクター

デフォルト = 1.0(実数)

[ Pa ]
ε ˙ 1 ひずみ速度1

(実数)

[ 1 s ]
ε ˙ 2 ひずみ速度2

(実数)

[ 1 s ]
ε ˙ 3 ひずみ速度3

(実数)

[ 1 s ]
ε ˙ 4 ひずみ速度4

(実数)

[ 1 s ]
ε ˙ 5 ひずみ速度5

(実数)

[ 1 s ]
ε ˙ 6 ひずみ速度6

(実数)

[ 1 s ]
ε ˙ 7 ひずみ速度7

(実数)

[ 1 s ]
ε ˙ 8 ひずみ速度8

(実数)

[ 1 s ]
ε ˙ 9 ひずみ速度9

(実数)

[ 1 s ]
ε ˙ 10 ひずみ速度10

(実数)

[ 1 s ]

例(アルミニウム)

#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/LAW60/1/1
Aluminium_example
#              RHO_I
               .0027 
#                  E                  Nu           Eps_p_max               Eps_t               Eps_m
               60400                 .33                   0                   0                   0
#  N_funct  F_smooth              C_hard               F_cut                
         4         0                   0                   0                    
#  fct_IDp              Fscale   Fct_IDE                EInf                  CE
         0                   0         0                   0                   0
# Funtions
         1         2         3         4
# Scale factors          Fscale_5
                   1                 1.2                 1.4                 1.6
# Strain rates          Eps_dot_5
                   0                  20                  30                  40
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#-  3. FUNCTIONS:
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/FUNCT/1
function_36
#                  X                   Y
                   0                  90                                                            
              2.5E-4                 100                                                           
                .001                 104                                                           
                .009                 121                                                            
                .017                 136                                                            
                .021                 143                                                           
                .036                 156                                                           
                .045                 162                                                         
                .055                 165                                                           
                .072                 170                                                           
                .075                 170                                                             
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/FUNCT/2
function_36
#                  X                   Y
                   0                  90                                                            
              2.5E-4                 100                                                           
                .001                 104                                                           
                .009                 121                                                            
                .017                 136                                                            
                .021                 143                                                           
                .036                 156                                                           
                .045                 162                                                         
                .055                 165                                                           
                .072                 170                                                           
                .075                 170                                                             
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/FUNCT/3
function_36
#                  X                   Y
                   0                  90                                                            
              2.5E-4                 100                                                           
                .001                 104                                                           
                .009                 121                                                            
                .017                 136                                                            
                .021                 143                                                           
                .036                 156                                                           
                .045                 162                                                         
                .055                 165                                                           
                .072                 170                                                           
                .075                 170                                                             
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/FUNCT/4
function_36
#                  X                   Y
                   0                  90                                                            
              2.5E-4                 100                                                           
                .001                 104                                                           
                .009                 121                                                            
                .017                 136                                                            
                .021                 143                                                           
                .036                 156                                                           
                .045                 162                                                         
                .055                 165                                                           
                .072                 170                                                           
                .075                 170                                                             
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#ENDDATA
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

コメント

  1. 降伏応力関数(塑性ひずみと応力の関係)の最初のポイントでは、塑性ひずみの値はゼロである必要があります。最初の(静的)関数の最後のポイントの応力が0に等しい場合、 ε p m a x のデフォルト値は ε p の対応値に設定されます。
  2. ε p (塑性ひずみ)が1つの積分点で ε p m a x に到達すると要素が削除されます。
  3. 最大主ひずみが ε 1 > ε t の場合、応力は次の関係式に従って減少します:(1)
    σ = σ ( ε m ε 1 ε m ε t )
  4. ε 1 > ε m の場合、要素は削除されます。
  5. 全体定式化では移動硬化モデルを使用できません(硬化は完全に等方性です)。
  6. 移動硬化およびひずみ速度依存の場合、降伏応力はひずみ速度に依存します。
  7. ひずみ速度フィルタリング入力(Fcut)は、シェル要素およびソリッド要素のみで使用できます。
  8. ひずみ速度のフィルタリングは、ひずみ速度のスムージングに使用されます。
  9. fct_IDpは、特定の材料における引張と圧縮の挙動の区別(すなわち圧力依存降伏)に使用されます。これは、ソリッド要素でのみ使用できます。したがって、有効降伏応力は公称降伏応力に実際の圧力に対応する降伏係数を乗じることによって得られます。


    図 1.
  10. ε ˙ n ε ˙ ε ˙ n + 1 の場合、降伏応力は、fn-1fnfn+1fn+2の各関数間での3次補間になります。
  11. ε ˙ ε ˙ 1 の場合、降伏応力は関数 f 1 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaciOzamaaBa aaleaacaaIXaaabeaaaaa@37C9@ f 2 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaciOzamaaBa aaleaacaaIXaaabeaaaaa@37C9@ および f 3 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaciOzamaaBa aaleaacaaIXaaabeaaaaa@37C9@ の間で補間されます。
  12. ε ˙ Nfunc 1 ε ˙ ε ˙ Nfunc の場合、降伏応力は、fNfunc-3fNfunc-2fNfunc-1およびfNfuncの各関数間で外挿されます。
  13. ε ˙ > ε ˙ Nfunc の場合、降伏応力は、関数fNfunc-2fNfunc-1およびfNfuncの間で外挿されます。

    mat_law60_yield
    図 2.
  14. ひずみ依存を記述する関数は、ひずみ速度値別に定義する必要があります。
  15. ひずみ速度の値は、厳密に昇順で指定する必要があります。
  16. ヤング率の進展:
    • fct_IDE > 0であれば、この曲線は等価な塑性ひずみによるヤング率の進展に対するスケールファクターを定義します。これは、関数 f ( ε ¯ p ) によって以下のようにヤング率がスケーリングされることを意味します:
      • E ( t ) = E f ( ε ¯ p )

        このスケールファクターの初期値は1で、この値から減少していきます。

      • fct_IDE = 0の場合、ヤング率は次のように計算されます:(2)
        E ( t ) = E ( E E inf ) [ 1 exp ( C E ε ¯ p ) ]

    ここで、

    EEinfはそれぞれ初期と漸近するヤング率の値で、 ε ¯ p は累積の相当塑性ひずみです。
    注:

    fct_IDE = 0でCE = 0の場合、ヤング率Eが一定になります。