/MAT/LAW6 (HYDROまたはHYD_VISC)

ブロックフォーマットキーワード 流体材料を記述します。圧力は、/EOSオプションの定義によって与えられる状態方程式を用いて計算されます。

フォーマット

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
/MAT/LAW6/mat_ID/unit_IDまたは/MAT/HYDRO/mat_ID/unit_IDまたは/MAT/HYD_VISC/mat_ID/unit_ID
mat_title
ρ i ρ 0            
ν Pmin            

定義

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

(整数、最大10桁)

 
unit_ID 単位識別子

(整数、最大10桁)

 
mat_title 材料のタイトル

(文字、最大100文字)

 
ρ i 初期密度

(実数)

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

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

[ kg m 3 ]
ν 運動粘性

(実数)

[ m 2 s ]
Pmin 圧力のカットオフ

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

[ Pa ]

例(空気)

#RADIOSS STARTER
/UNIT/1
unit for mat
                  kg                   m                   s
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/MAT/HYDRO/4/1
AIR
#              RHO_I               RHO_0
                1.22                   0 
#                Knu                Pmin
              1.5E-5                   0
/EOS/POLYNOMIAL/4/1
AIR
#                 C0                  C1                  C2                  C3
                   0                   0                   0                   0
#                 C4                  C5                  E0                 Psh               RHO_0
                 0.4                 0.4              253300                   0                1.22
/EULER/MAT/4
#     Modif. factor.
                   0
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#enddata
/END
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

コメント

  1. S ij = 2 ρ ν eq e ˙ ij
    ここで、
    ν eq = ν
    乱流なし
    S i j
    偏差応力テンソル
    e ˙ ij
    偏差ひずみテンソル
  2. 流体力学的圧力の状態方程式は、/EOSカードで指定されていなければなりません。
  3. 体積膨張を伴う線形材料の場合:

    C 1 = E 3 ( 1 2 ν ) および C 4 = α ν C 1 ρ C ν T

    C 4 = C 5 = γ 1 および C 0 = C 2 = C 3 = 0

    したがって:(1)
    p = C 1 μ + ( C 4 + C 5 μ ) E = C 1 μ + C 4 ( 1 + μ ) E = C 1 μ + C 4 ( 1 + μ ) ρ 0 e = C 1 μ + C 4 ( 1 + μ ) ρ 0 C ν T
    (2)
    p = C 1 μ + C 4 ρ C ν T = C 1 μ + α ν T

    p = c s t = 0 の場合、 C 1 μ + α ν T = 0 l=、したがって μ = α ν T C 1

    ここで、
    μ
    膨張係数
    μ < 0
    膨張
    この場合パラメータC2C3は考慮されません。
  4. すべての熱データ( ρ 0 C p , T 0 , A , and B )はキーワード/HEAT/MATで定義できます。
  5. 液相(気相を含まない)に対してLAW37とLAW6を組み合わせて使用する場合は、液体EOSの適合性が以下のようになります:
    • Δ P 1 = C 1 μ --右記の場合; /MAT/LAW37 (BIPHAS)
    • p = C 0 + C 1 μ + C 2 μ 2 + C 3 μ 3 + ( C 4 + C 5 μ ) E --LAW6の場合 (上記の例で定義されている多項式EOSを介し)

    ここで、 C 0 = C 2 = C 3 = C 4 = C 5 =E=0 MathType@MTEF@5@5@+= feaahqart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGdbWaaS baaSqaaiaaicdaaeqaaOGaeyypa0Jaam4qamaaBaaaleaacaaIYaaa beaakiabg2da9iaadoeadaWgaaWcbaGaaG4maaqabaGccqGH9aqpca WGdbWaaSbaaSqaaiaaisdaaeqaaOGaeyypa0Jaam4qamaaBaaaleaa caaI1aaabeaakiabg2da9iaadweacqGH9aqpcaaIWaaaaa@46A9@

    したがって: p = C 1 μ

  6. 気相(液相を含まない)に対してLAW37とLAW6を組み合わせて使用する場合は、気体EOSの適合性が以下のようになります:
    • P V γ = c o n s t . --LAW37の場合
    • p = ( γ 1 ) ( μ + 1 ) E --LAW6の場合(/EOS/IDEAL-GAS状態方程式を介し)

    ここで、 E MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGfbaaaa@3832@ は単位体積当たりのエネルギーです。