/EOS/MURNAGHAN

ブロックフォーマットのキーワード Murnaghan状態方程式を記述します。

フォーマット

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
/EOS/MURNAGHAN/mat_ID/unit_ID
eos_title
K0 K1 P0 Psh ρ 0

定義

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

(整数、最大10桁)

 
unit_ID 単位識別子

(整数、最大10桁)

 
eos_title 状態方程式のタイトル

(文字、最大100文字)

 
K0 材料パラメータ

(実数)

[ Pa ]
K1 材料パラメータ

(実数)

 
P0 初期圧力

(実数)

[ Pa ]
Psh 圧力シフト

(実数)

[ Pa ]
ρ 0 参照密度

デフォルト = 材料密度(実数)

[ kg m 3 ]

#---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----|
/MAT/HYDPLA/7
Articficial Linear Material Law
#              RHO_I               RHO_0
             1.22e-3             1.22e-3
#                  E                  nu
                   0                   0
#                  a                   b                   n             eps_max           sigma_max
                1E30                   0                   0                   0                   0
#               Pmin 
                   0
/EOS/MURNAGHAN/7
EoS for NaCl  at atmospheric pressure
#                 K0                  K1                  P0                 PSH                RHO0      
               24000               5.390                  .1                   0            2.165e-3
/ALE/MAT/7

#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#enddata

コメント

  1. この状態方程式は、Tait状態方程式としても知られています:(1)
    P ( V ) = K 0 K 1 [ ( V V 0 ) K 1 1 ]

    ここで、 K 0 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4samaaBa aaleaacaaIWaaabeaaaaa@37AD@ K 1 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4samaaBa aaleaacaaIXaaabeaaaaa@37AE@ は材料パラメータです。

  2. この式は、次のような形でも求めることができます:(2)
    Δ v V 0 = 1 [ 1 + K 1 K 0 p ] 1 K 1 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaSaaaeaacq GHuoarcaWG2baabaGaamOvamaaBaaaleaacaaIWaaabeaaaaGccqGH 9aqpcaaIXaGaeyOeI0YaamWaaeaacaaIXaGaey4kaSYaaSaaaeaaca WGlbWaaSbaaSqaaiaaigdaaeqaaaGcbaGaam4samaaBaaaleaacaaI WaaabeaaaaGccaWGWbaacaGLBbGaayzxaaWaaWbaaSqabeaadaWcca qaaiabgkHiTiaaigdaaeaacaWGlbWaaSbaaWqaaiaaigdaaeqaaaaa aaaaaa@4895@

    ここで、 Δ v = V 0 V MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeyiLdqKaam ODaiabg2da9iaadAfadaWgaaWcbaGaaGimaaqabaGccqGHsislcaWG wbaaaa@3CF2@ および p=P P 0 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiCaiabg2 da9iaadcfacqGHsislcaWGqbWaaSbaaSqaaiaaicdaaeqaaaaa@3B6F@

  3. 一部の出版物では、材料パラメータ K 0 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4samaaBa aaleaacaaIWaaabeaaaaa@37AD@ および K 1 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4samaaBa aaleaacaaIXaaabeaaaaa@37AE@ c MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4yaaaa@36DF@ および k MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4yaaaa@36DF@ によって K 0 = c MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4samaaBa aaleaacaaIWaaabeaakiabg2da9iaadogaaaa@39A5@ および K 1 = c × k MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4samaaBa aaleaacaaIXaaabeaakiabg2da9iaadogacqGHxdaTcaWGRbaaaa@3CAD@ に置き換えられています。
  4. この式を表現する別の方法として、圧縮率 μ を使用する方法があります。(3)
    P ( μ ) = P 0 + K 0 K 1 [ ( 1 + μ ) K 1 1 ] MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaaciGGqbWaae WaaeaacqaH8oqBaiaawIcacaGLPaaacqGH9aqpcaWGqbWaaSbaaSqa aiaaicdaaeqaaOGaey4kaSYaaSaaaeaacaWGlbWaaSbaaSqaaiaaic daaeqaaaGcbaGaam4samaaBaaaleaacaaIXaaabeaaaaGcdaWadaqa amaabmaabaGaaGymaiabgUcaRiabeY7aTbGaayjkaiaawMcaamaaCa aaleqabaGaam4samaaBaaameaacaaIXaaabeaaaaGccqGHsislcaaI XaaacaGLBbGaayzxaaaaaa@4C17@

    ここで、 μ= ρ ρ 0 1= V 0 V 1 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqiVd0Maey ypa0ZaaSaaaeaacqaHbpGCaeaacqaHbpGCdaWgaaWcbaGaaGimaaqa baaaaOGaeyOeI0IaaGymaiabg2da9maalaaabaGaamOvamaaBaaale aacaaIWaaabeaaaOqaaiaadAfaaaGaeyOeI0IaaGymaaaa@443F@

  5. Murnaghan EOSはエネルギーに依存しません。
  6. Radiossにより流体力学的圧力の計算に用いられ、右記の材料則と適合性のある状態方程式。
    • /MAT/LAW3 (HYDPLA)
    • /MAT/LAW4 (HYD_JCOOK)
    • /MAT/LAW6 (HYDROまたはHYD_VISC)
    • /MAT/LAW10 (DPRAG1)
    • /MAT/LAW12 (3D_COMP)
    • /MAT/LAW49 (STEINB)
    • /MAT/LAW102 (DPRAG2)
    • /MAT/LAW103 (HENSEL-SPITTEL)
1 Murnaghan, F. D. "The compressibility of media under extreme pressures." Proceedings of the National Academy of Sciences 30, no. 9 (1944): 244-247