/MAT/LAW16 (GRAY)

ブロックフォーマットキーワード この材料則は、Gray EOSとJohnson-Cookの降伏基準に基づきます。

フォーマット

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
/MAT/LAW16/mat_IDまたは/MAT/GRAY/mat_ID
mat_title
ρ i ρ 0            
E ν            
a b n ε max σ max
P0 C S γ 0 α MathType@MTEF@5@5@+= feaahqart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqySdegaaa@3792@ e
AW Pmin E0        
c ε ˙ 0 m Tmelt Tmax
γ 0 m α MathType@MTEF@5@5@+= feaahqart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqySdegaaa@3792@ m γ e ge Δ S MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeuiLdqKaam 4uaaaa@3834@
Tm0 Vj Vb        
  Eoh Ay θ    

定義

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

(整数、最大10桁)

 
mat_title 材料のタイトル

(文字、最大100文字)

 
ρ i 初期密度

(実数)

[ kg m 3 ]
ρ 0 E.O.S(状態方程式)で使用される基準密度

デフォルト = ρ 0 = ρ i (実数)

[ kg m 3 ]
E ヤング率

(実数)

[ Pa ]
ν ポアソン比

(実数)

 
a 塑性降伏応力。

(実数)

[ Pa ]
b 塑性硬化パラメータ。

(実数)

[ Pa ]
n 塑性硬化指数。

デフォルト = 1.0001(実数)

 
ε max 破壊塑性ひずみ。

デフォルト = 1030(実数)

 
σ max 塑性最大応力

デフォルト = 1030(実数)

[ Pa ]
P0 初期圧力

(実数)

[ Pa ]
C Hugoniotパラメータ

(実数)

[ m s ] MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaamWaaeaada Wcaaqaaiaab2gaaeaacaqGZbaaaaGaay5waiaaw2faaaaa@39DE@
S U s = C + S U p

(実数)

 
γ 0 格子ガンマ

(実数)

 
σ e γ = γ 0 a x

(実数)

 
AW 原子量

(実数)

[ k g m o l e ] MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeaaciGaaiaabeqaamaabaabaaGcbaWaamWaaeaada WcaaqaaiaadUgacaWGNbaabaGaamyBaiaad+gacaWGSbGaamyzaaaa aiaawUfacaGLDbaaaaa@3D93@
Pmin 圧力のカットオフ

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

[ Pa ]
E0 単位体積あたりの初期エネルギー

(実数)

[ J m 3 ]
c ひずみ速度係数。
= 0(デフォルト)
ひずみ速度効果はなし

(実数)

 
ε ˙ 0 参照ひずみ速度(時間単位)-1

(実数)

[ 1 s ]
m 温度指数。

デフォルト = 1.0(実数)

 
Tmelt 溶融温度。

デフォルト = 1030(実数)

[ K ]
Tmax T > Tmaxの場合:m = 1 が使用されます。

デフォルト = 1030(実数)

[ K ]
γ 0m 溶融ガンマ

デフォルト = γ 0(実数)

 
σ m γ m = γ 0 m a m x

デフォルト = a (実数)

 
γ e 電子ガンマ

デフォルト = 2/3(実数)

 
ge 電子エネルギー係数

デフォルト = 0.0(実数)

[ J m o l e K 2 ] MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeaaciGaaiaabeqaamaabaabaaGcbaWaamWaaeaada WcaaqaaiaadQeaaeaacaWGTbGaam4BaiaadYgacaWGLbGaeyyXICTa am4samaaCaaaleqabaGaaGOmaaaaaaaakiaawUfacaGLDbaaaaa@4093@
Δ S MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeuiLdqKaam 4uaaaa@3834@ 溶融エントロピー

デフォルト = 9.637e+3 [ J kgK ] MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeaaciGaaiaabeqaamaabaabaaGcbaWaamWaaeaada WcaaqaaiaabQeaaeaacaqGRbGaae4zaiabgwSixlaabUeaaaaacaGL BbGaayzxaaaaaa@3DB3@ (実数)

空白または0の場合は、溶融エントロピーのデフォルト値が自動的に作業単位に変換されます。

[ J kgK ] MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeaaciGaaiaabeqaamaabaabaaGcbaWaamWaaeaada WcaaqaaiaabQeaaeaacaqGRbGaae4zaiabgwSixlaabUeaaaaacaGL BbGaayzxaaaaaa@3DB3@
Tm0 溶融温度パラメータ

デフォルト = 1.3 Tmelt(実数)

[ K ]
Vj EOS結合箇所の体積

(実数)

[ m 3 kg ] MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeaaciGaaiaabeqaamaabaabaaGcbaWaamWaaeaada Wcaaqaaiaac2gadaahaaWcbeqaaiaacodaaaaakeaacaGGRbGaai4z aaaaaiaawUfacaGLDbaaaaa@3BAF@
Vb 蒸気相の除外体積

デフォルト = 0.5/r0(実数)

[ m 3 k g ] MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeaaciGaaiaabeqaamaabaabaaGcbaWaamWaaeaada Wcaaqaaiaac2gadaahaaWcbeqaaiaacodaaaaakeaacaGGRbGaai4z aaaaaiaawUfacaGLDbaaaaa@3BAF@
Eoh V=V0T=300K、P=0におけるエネルギー

デフォルト = 0.0(実数)

[ J ]
Ay 引力ポテンシャル係数

(実数)

[ J m 3 mole 2 ] MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaaqaaaaaaaaa WdbmaadmaapaqaamaalaaabaGaciOsaiaab2gadaahaaWcbeqaaiaa iodaaaaakeaacaqGTbGaae4BaiaabYgacaqGLbWaaWbaaSqabeaaca aIYaaaaaaaaOWdbiaawUfacaGLDbaaaaa@4000@
θ 結合パラメータ

デフォルト = 1.0(実数)

 

例(アルミニウム)

#RADIOSS STARTER
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/MAT/LAW16/1
Aluminium - GRAY EOS (unit g_cm_mus)
#              RHO_I               RHO_0
                 2.7                   0
#                  E                  Nu
                0.71                 .33
#                  a                   b                   n             Eps_max              sigmax
               0.002             0.00144                0.62                   0                   0
#                 P0                   C                   S             Gamma_0                 a_e
                   0               0.533               1.338                2.18                 1.7
#                 AW               P_min                  E0
               26.98               -5e-3
#                  c           eps_dot_0                   m              T_melt               T_max
                                                           1                1E30
#           Gamma_0m                 a_m             Gamma_e                 g_e             delta_S
                2.18                 1.7              0.6667              8.7e-9             9.637E5
#               T_m0                 V_j                 V_b
                1220             0.47388             0.19025            
#                                   E_0h                 a_y               Theta
                                       0                  47                   1
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#ENDDATA
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

コメント

  1. 正しいモデル単位を/BEGINカードに定義する必要があります。単位系に基づいて溶融エントロピー Δ S MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeuiLdqKaam 4uaaaa@3834@ のデフォルト値が計算されるからです。
    (1)
    σ = ( a + b ε p n ) ( 1 + c ln ε ˙ ε ˙ 0 ) ( 1 ( T * ) m ) T * = T T 0 T melt T 0
1 E. B. Royce, GRAY, a three-phase equation of state, UCRL-51121, Lawrence Livermore Laboratory, 1971