/EOS/LINEAR

ブロックフォーマットキーワード 初期圧力と体積弾性率を擁する線形状態方程式 P ( μ ) = P 0 + B μ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGqbWaae WaaeaacqaH8oqBaiaawIcacaGLPaaacqGH9aqpcaWGqbWaaSbaaSqa aiaaicdaaeqaaOGaey4kaSIaamOqaiabeY7aTbaa@41A6@ を記述します。

フォーマット

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
/EOS/LINEAR/mat_ID/unit_ID
eos_title
P 0 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGqbWaaS baaSqaaiaaicdaaeqaaaaa@3923@ B Psh ρ 0  

定義

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

(整数、最大10桁)

 
unit_ID 単位識別子

(整数、最大10桁)

 
eos_title 状態方程式のタイトル。

(文字、最大100文字)

 
P 0 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGqbWaaS baaSqaaiaaicdaaeqaaaaa@3923@ 初期圧力

(実数)

[ Pa ]
B 体積弾性率。

(実数)

[ Pa ]
Psh 圧力シフト

(実数)

[ Pa ]
ρ 0 参照密度

(実数)

[ kg m 3 ]

#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/UNIT/1
unit for mat
                  kg                   m                   s
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/MAT/LAW06/7/1
law06
#              RHO_I        
             2.33e-3          
#                 NU                PMIN
                0.22               -0.02
/EOS/LINEAR/7/1
linear EOS (Artificial data)
#                 P0                   B                 PSH                RHO0          
                   1                10.0                 0.0             2.33e-3  
/ALE/MAT/7
#               Flrd
                   0
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#enddata

コメント

  1. 線形状態方程式は、以下のフォームを有します:(1)
    P ( μ ) = P 0 + B μ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGqbWaae WaaeaacqaH8oqBaiaawIcacaGLPaaacqGH9aqpcaWGqbWaaSbaaSqa aiaaicdaaeqaaOGaey4kaSIaamOqaiabeY7aTbaa@41A6@
    ここで、(2)
    µ = ρ ρ 0 1 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaaqaaaaaaaaa WdbiaadwlacqGH9aqpdaWcaaWdaeaapeGaeqyWdihapaqaa8qacqaH bpGCpaWaaSbaaSqaa8qacaaIWaaapaqabaaaaOWdbiabgkHiTiaaig daaaa@3F63@
    これは下記の多項式状態方程式から導き出されることが可能:(3)
    P = C 0 + C 1 μ + C 2 μ 2 + C 3 μ 3 + ( C 4 + C 5 μ ) E 0 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiuaiabg2 da9iaadoeadaWgaaWcbaGaaGimaaqabaGccqGHRaWkcaWGdbWaaSba aSqaaiaaigdaaeqaaOGaeqiVd0Maey4kaSIaam4qamaaBaaaleaaca aIYaaabeaakiabeY7aTnaaCaaaleqabaGaaGOmaaaakiabgUcaRiaa doeadaWgaaWcbaGaaG4maaqabaGccqaH8oqBdaahaaWcbeqaaiaaio daaaGccqGHRaWkdaqadaqaaiaadoeadaWgaaWcbaGaaGinaaqabaGc cqGHRaWkcaWGdbWaaSbaaSqaaiaaiwdaaeqaaOGaeqiVd0gacaGLOa GaayzkaaGaamyramaaBaaaleaacaaIWaaabeaaaaa@5292@
    ここで、
    C 0 = P 0 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4qamaaBa aaleaacaaIWaaabeaakiabg2da9iaadcfadaWgaaWcbaGaaGimaaqa baaaaa@3A6F@
    C 1 = B MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4qamaaBa aaleaacaaIXaaabeaakiabg2da9iaadkeaaaa@397C@
    C 2 = C 3 = C 4 = C 5 = 0 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4qamaaBa aaleaacaaIYaaabeaakiabg2da9iaadoeadaWgaaWcbaGaaG4maaqa baGccqGH9aqpcaWGdbWaaSbaaSqaaiaaisdaaeqaaOGaeyypa0Jaam 4qamaaBaaaleaacaaI1aaabeaakiabg2da9iaaicdaaaa@41B6@
  2. 体積弾性率は通常、次のように推定されます:(4)
    B = ρ 0 c 0 2 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOqaiabg2 da9iabeg8aYnaaBaaaleaacaaIWaaabeaakiabgwSixlaadogadaWg aaWcbaGaaGimaaqabaGcdaahaaWcbeqaaiaaikdaaaaaaa@3F7E@

    ここで、 c 0 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4yamaaBa aaleaacaaIWaaabeaaaaa@37C4@ は初期音速。

  3. Psh パラメータは、出力圧力のシフトを可能にします。出力圧力は、相対圧力 Δ P ( μ ) = P P s h MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeuiLdqKaam iuamaabmaabaGaeqiVd0gacaGLOaGaayzkaaGaeyypa0Jaamiuaiab gkHiTiaadcfadaWgaaWcbaGaam4CaiaadIgaaeqaaaaa@411E@ にもなり得ます。
  4. 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)