/BEM/DAA

ブロックフォーマットキーワード 流体の質量マトリックスが境界要素法で計算される水中爆発の二重漸近近似。

フォーマット

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
/BEM/DAA/daa_ID/unit_ID
daa_title
surf_ID grav_ID            
ρ C Pmin        
Xs Ys Zs        
Iform Ipri Ipres     Kform Freesurf Afterflow Integr  
Ipres=1の場合挿入
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
Pm θ a a θ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGHbWaaS baaSqaaiabeI7aXbqabaaaaa@3A30@    
Ipres=2の場合挿入
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
fct_IDP   FscaleP        
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
Xc Yc Zc        
grav_ID > 0またはFreesurf=2の場合挿入
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
XA YA ZA        
Dir-X Dir-Y Dir-Z        

定義

フィールド 内容 SI 単位の例
daa_ID DAAブロック識別子

(整数、最大10桁)

 
unit_ID 単位識別子

(整数、最大10桁)

 
daa_title DAAブロックタイトル

(文字、最大100文字)

 
surf_ID 湿潤サーフェス識別子 2 3

(整数)

 
grav_ID /GRAVオプション識別子

(整数)

 
ρ 流体密度。

(実数)

[ kg m 3 ]
C 流体の音速

(実数)

[ m s ] MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaamWaaeaada Wcaaqaaiaab2gaaeaacaqGZbaaaaGaay5waiaaw2faaaaa@39DE@
Pmin 圧力のカットオフ( < 0 )

デフォルト = -1030(実数)

[ Pa ]
Xs スタンドオフポイントのX座標 3

(実数)

[ m ]
Ys スタンドオフポイントのY座標 3

(実数)

[ m ]
Zs スタンドオフポイントのZ座標 3

(実数)

[ m ]
Iform BEMソリューションフラグ
= 1(デフォルト)
ガウス積分
= 2
解析的積分

(整数)

 
Ipri 出力フラグレベル
= 1(デフォルト)
低減出力
= 2
フル出力

(整数)

 
Ipres 圧力荷重フラグ 6
=1
圧力は、 P i ( t ) = P m e t θ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaaciGGqbWaaS baaSqaaiaacMgaaeqaaOWaaeWaaeaacaWG0baacaGLOaGaayzkaaGa eyypa0JaamiuamaaBaaaleaacaWGTbaabeaakiaadwgadaahaaWcbe qaaiabgkHiTmaalaaabaGaamiDaaqaaiabeI7aXbaaaaaaaa@43A9@ として計算されます。
= 2
関数による入力

(整数)

 
Kform 解析フラグ
= 1(デフォルト)
DAA1定式化
= 2
高周波数

(整数)

 
Freesurf 自由サーフェスフラグ 6
= 1(デフォルト)
×
= 2
あり

(整数)

 
Afterflow アフターフローの計算 7
=1
×
= 2(デフォルト)
あり

(整数)

 
Integr 時間整数フラグ
=1
1次
= 2(デフォルト)
2次

(整数)

 
Pm 最大圧力 5

(実数)

[ Pa ]
θ 減衰時間

(実数)

[ s ]
a 最大圧力定数 5

(実数)

 
a θ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGHbWaaS baaSqaaiabeI7aXbqabaaaaa@3A30@ 圧力減衰時間定数 5

(実数)

 
fct_IDP 入射圧力関数識別子

(整数)

 
FscaleP fct_IDPの縦軸(圧力)のスケールファクター

(実数)

[ Pa ]
XC 装薬のX座標

(実数)

[ m ]
YC 装薬のY座標

(実数)

[ m ]
ZC 装薬のZ座標

(実数)

[ m ]
XA 自由サーフェス上のポイントAのX座標

(実数)

[ m ]
YA 自由サーフェス上のポイントAのY座標

(実数)

[ m ]
ZA 自由サーフェス上のポイントAのZ座標

(実数)

[ m ]
Dir-X 自由サーフェス平面への法線のX座標

(実数)

 
Dir-Y 自由サーフェス平面への法線のY座標

(実数)

 
Dir-Z 自由サーフェス平面への法線のZ座標

(実数)

 

コメント

  1. 構造全体がモデル化される必要があります。対称解析はサポートされません。
  2. サーフェスの法線 n MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWHUbaaaa@385F@ は、流体に向いている必要があります。
  3. (Xs, Ys, Zs)で定義されたスタンドオフポイントは、入射圧力波が時間=0において与えられた位置です:


    図 1.
  4. 平面波は球波を使用し、装薬を十分離して置いてシミュレートすることができます。
  5. 時間の関数としてのスタンドオフポイントにおける圧力は:(1)
    P i ( t ) = P m e t θ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaaciGGqbWaaS baaSqaaiaacMgaaeqaaOWaaeWaaeaacaWG0baacaGLOaGaayzkaaGa eyypa0JaamiuamaaBaaaleaacaWGTbaabeaakiaadwgadaahaaWcbe qaaiabgkHiTmaalaaabaGaamiDaaqaaiabeI7aXbaaaaaaaa@43A9@
    ここで、
    P m MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGqbWaaS baaSqaaiaad2gaaeqaaaaa@395B@
    最大圧力
    t MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWG0baaaa@3861@
    時間
    θ
    減衰時間
    最大圧力と減衰時間は、以下を用いて計算できます:(2)
    P m = K [ W 1 3 R ] a MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGqbWaaS baaSqaaiaad2gaaeqaaOGaeyypa0Jaam4samaadmaabaWaaSaaaeaa caWGxbWaaWbaaSqabeaadaWccaqaaiaaigdaaeaacaaIZaaaaaaaaO qaaiaadkfaaaaacaGLBbGaayzxaaWaaWbaaSqabeaacaWGHbaaaaaa @41C4@
    (3)
    θ = K θ W 1 3 [ W 1 3 R ] a θ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaH4oqCcq GH9aqpcaWGlbWaaSbaaSqaaiabeI7aXbqabaGccaWGxbWaaWbaaSqa beaadaWccaqaaiaaigdaaeaacaaIZaaaaaaakmaadmaabaWaaSaaae aacaWGxbWaaWbaaSqabeaadaWccaqaaiaaigdaaeaacaaIZaaaaaaa aOqaaiaadkfaaaaacaGLBbGaayzxaaWaaWbaaSqabeaacaWGHbWaaS baaWqaaiabeI7aXbqabaaaaaaa@47E9@
    W MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGxbaaaa@3844@
    爆発物質量
    R MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGsbaaaa@383F@
    爆発までの距離
    K MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGlbaaaa@3838@ α MathType@MTEF@5@5@+= feaahqart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqySdegaaa@3792@ K θ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGlbWaaS baaSqaaiabeI7aXbqabaaaaa@3A1A@ and a θ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGHbWaaS baaSqaaiabeI7aXbqabaaaaa@3A30@
    爆薬に依存する定数
    W MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGxbaaaa@3844@ (kg)、 R MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGsbaaaa@383F@ (m)、 P m MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGqbWaaS baaSqaaiaad2gaaeqaaaaa@395B@ (MPaおよびms)の場合。
      K MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGlbaaaa@3838@ α MathType@MTEF@5@5@+= feaahqart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqySdegaaa@3792@ K θ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGlbWaaS baaSqaaiabeI7aXbqabaaaaa@3A1A@ a θ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGHbWaaS baaSqaaiabeI7aXbqabaaaaa@3A30@
    TNT 52.12 1.180 0.0895 -0.185
    PETN 56.21 1.194 0.0860 -0.257
    HBX 53.51 1.144 0.0920 -0.247
  6. 自由サーフェスは、1つのポイントとその法線ベクトルによって定義される平面です。
  7. アフターフロー法線速度は次のように計算されます:(4)
    v a f t e r f l o w = cos γ ρ R 0 t P ( t ) d t MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWG2bWaaS baaSqaaiaadggacaWGMbGaamiDaiaadwgacaWGYbGaamOzaiaadYga caWGVbGaam4DaaqabaGccqGH9aqpdaWcaaqaaiGacogacaGGVbGaai 4Caiabeo7aNbqaaiabeg8aYjaadkfaaaWaa8qmaeaacaWGqbWaaeWa aeaacaWG0baacaGLOaGaayzkaaGaamizaiaadshaaSqaaiaaicdaae aacaWG0baaniabgUIiYdaaaa@524C@
    P
    流体ポイント
    C
    装薬ポイント
    S
    スタンドオフポイント


    図 2.
1 Littlewood, J. de Runtz T. 2004."USA Code".Mecalog Workshop, Sophia Antipolis, France, 2004
2 Cole, Robert H. Underwater Explosion.Princeton University Press, 1948