/EOS/TILLOTSON

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

フォーマット

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
/EOS/TILLOTSON/mat_ID/unit_ID
eos_title
C1 C2 a b  
ER ES VS E0 ρ 0
α β      

定義

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

（整数、最大10桁）

  
単位識別子 単位識別子。

（整数、最大10桁）

  
eos_title 状態方程式のタイトル

（文字、最大100文字）

  
C1 C1係数

（実数）

 [ Pa ]
C2 C2係数

（実数）

 [ Pa ]
a A係数

（実数）

  
b b係数

（実数）

  
ER 単位参照体積あたりの内部エネルギー

（実数）

 [ J m 3 ]
ES 単位参照体積あたりの昇華エネルギー

（実数）

 [ J m 3 ]
VS 昇華相対体積

（実数）

 [ m 3 ]
E0 単位参照体積あたりの初期エネルギー

（実数）

 [ J m 3 ]
ρ 0 参照密度

デフォルト = 材料密度（実数）

 [ kg m 3 ]
α α 係数

（実数）

  
β β 係数

（実数）

  

例（アルミニウム）

#RADIOSS STARTER
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/UNIT/1
unit for mat
                   g                  cm                 mus
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#-  2. MATERIALS:
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/MAT/HYD_JCOOK/1/1
Aluminium
#              RHO_I               RHO_0              
                 2.8                   0                   
#                 E                   nu
                .734                 .33
#                  A                   B                   n              epsmax              sigmax
               .0024               .0042                  .8                   0               .0068
#               Pmin
              -.0223
#                  C           EPS_DOT_0                   M               Tmelt               Tmax
                .062                1E-6                   1                1220                   0
#              RHOCP                                                         T_r
             2.59E-5                                                           0
/EOS/TILLOTSON/1/1
Aluminium
#                 C1                  C2                   A                   B
                .752                 .65                  .5                1.63
#                 ER                  ES                  VS                  E0               RHO_0
                .135                .081                 1.1                   0                   0
#              ALPHA                BETA
                   5                   5
/FAIL/JOHNSON/3
#                 D1                  D2                  D3                  D4                  D5
                .112                .123                -1.5                .007                   0
#              EPS_0  Ifail_sh  Ifail_so                                    Dadv               Ixfem
                1E-6         0         1                                       0                   0
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#ENDDATA
/END
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

入力の例

"Metallic Equation of State for Hypervelocity Impact" J.H.Tillotson、General Dynamics社、1962年：
材料 ρ 0 [ g c m 3 ] g/cm3 C1 [ Mbar ] C2 [ Mbar ] a b Er [ Mbar ] α β Es [ Mbar ] Vs
Cu 8.9 1.390 1.10 0.5 1.50 2.892 5 5 0.123 1.18
Fe 7.8 1.279 1.05 0.5 1.50 0.741 5 5 0.190 1.21
Al 2.7 0.752 0.65 0.5 1.63 0.135 5 5 0.081 1.10

コメント

  1. ここで、 μ = ρ ρ 0 − 1

    V = 1 ρ 比容積

    η = 1 + μ

    x = 1 − ρ 0 ρ

    およびEは単位参照体積当たりの初期エネルギー

    圧力は以下で定義されます:

    領域1： μ ≥ 0 (1)
    P = C 1 μ + C 2 μ 2 + ( a + b ω ) η E

    ここで、 ω = 1 + E E r η 2

    領域2： μ < 0 V V 0 < V s および E < E s (2)
    P = C 1 μ + ( a + b ω ) η E
    領域3： μ < 0 、 V V 0 > V s または V V 0 < V s および E ≥ E s (3)
    P = C 1 e β x e − α x 2 μ + ( a + b e − α x 2 ω ) η E
  2. 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)