/INIGRAV

ブロックフォーマットキーワード 重力荷重による静水圧から初期状態を定義します。

フォーマット

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
/INIGRAV/inigrav_ID/unit_ID
inigrav_title
grpart_ID surf_ID grav_ID              
Pref Bx By Bz    

定義

フィールド 内容 SI 単位の例
inigrav_ID 初期重力荷重識別子

(整数、最大10桁)

 
unit_ID 単位の識別子(オプション)

(整数、最大10桁)

 
inigrav_title 初期重力荷重タイトル

(文字、最大100文字)

 
grpart_ID 特定されたパートのグループ定義されない場合(0)、初期状態はすべてのソリッドPARTに適用されます。

デフォルト = 0(整数)

 
surf_ID 自由サーフェスの識別子(オプション)

定義されていない場合、基礎点が読まれます。

(整数)

 
grav_ID 使用される重力カード識別子(必須)

(整数)

 
Pref 参照サーフェス上の参照圧力

(実数)

[ Pa ]
Bx 基礎点のX座標(surf_ID=0の場合にのみ使用)

(実数)

[ m ]
By 基礎点のY座標(surf_ID=0の場合にのみ使用)

(実数)

[ m ]
Bz 基礎点のZ座標(surf_ID=0の場合にのみ使用)

(実数)

[ m ]

例(水の静水圧)

#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/UNIT/1
unit for load
#              MUNIT               LUNIT               TUNIT
                   g                  mm                  ms
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/INIGRAV/1/1
hydrostatique pressure - “ sea level : {Z=0} where Pref=0.101325 MPa
#GRPART_ID   SURF_ID   GRAV_ID    
         0         0         1  
#               Pref                  Bx                  By                  Bz          
            0.101325                   0                   0                   0
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/GRAV/1/1
Terrestrial gravity field
#funct_IDT       DIR   skew_ID sensor_ID  grnod_ID                      Ascale_x            Fscale_Y
         2         Z         0         0         5                             0            -0.00981
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/FUNCT/2
unity
#                  X                   Y
                   0                 1.0                                                            
             1000000                 1.0                                                            
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#enddata

コメント

  1. 重力値と方向は/GRAVカードから読み出され、grav_ID値は必須となります。
  2. サーフェスは、サーフェス識別子surf_IDが与えられない限り、基礎点B=(Bx, By, Bz)と垂直ベクトル gn MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGNbGaey yXICTaaCOBaaaa@3B95@ により定義されます。2018.0リリースで、surf_IDは現時点において/SURF/PLANEオプションと適合性があります。
  3. Pref は参照サーフェス上で測られる参照圧力で、通常これは海抜ゼロにおける圧力です。サーフェス圧力より上であれば減少、下であれば増加しています。圧力変化は、参照サーフェスからの符号付距離に依存します:
    (1)
    Δ P 0 = ρ 0 | g | d i s t MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqGHuoarca WGqbWaaSbaaSqaaiaaicdaaeqaaOGaeyypa0JaeqyWdi3aaSbaaSqa aiaaicdaaeqaaOWaaqWaaeaacaWGNbaacaGLhWUaayjcSdGaamizai aadMgacaWGZbGaamiDaaaa@4620@
    ここで、
    d i s t = n d i s t MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGKbGaam yAaiaadohacaWG0bGaeyypa0JaaCOBaiabgwSixlaahsgacaWHPbGa aC4Caiaahshaaaa@434F@
    符号付距離
    ρ 0
    要素の参照密度
    符号付距離は、要素の重心から参照サーフェスまでで計算されます。指定の要素での初期圧力は:(2)
    P = P r e f + ρ 0 ( d i s t ) | g | d i s t MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGqbGaey ypa0JaamiuamaaBaaaleaacaWGYbGaamyzaiaadAgaaeqaaOGaey4k aSIaeqyWdi3aaSbaaSqaaiaaicdaaeqaaOWaaeWaaeaacaWGKbGaam yAaiaadohacaWG0baacaGLOaGaayzkaaWaaqWaaeaacaWGNbaacaGL hWUaayjcSdGaamizaiaadMgacaWGZbGaamiDaaaa@4DD3@
    ここで、 ρ 0 ( d i s t ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaHbpGCda WgaaWcbaGaaGimaaqabaGcdaqadaqaaiaadsgacaWGPbGaam4Caiaa dshaaiaawIcacaGLPaaaaaa@3F69@ は指定の要素の参照密度。


    図 1.
  4. 同じ圧力値が、要素のすべての積分点に適用されます。
  5. /INIVOLオプションがLAW51で体積比率の初期化に使用される場合、参照密度は初期体積比率から自動的に計算されます。
  6. 現時点では、流体力学則3、4、6、37、49、51(ALE、EULER)とのみ適合性を有します。
  7. grav_IDカードは、重力/GRAVをt=0からの一定重力関数fct_IDTと共に参照する必要があります。または、関数は空白のままにして、一定の重力値を/GRAV, FscaleYで定義することができます。


    図 2.