/LOAD/PBLAST

エアブラストの入射圧力は、実験データからフィッティングされ、爆風圧力はサーフェスの向きから爆発位置へ導かれます。ユーザーは、爆発位置、爆発時間および相当TNT質量を与える必要があります。

これは、到着時間と入射圧力が障害物について調整されないため、簡素化された載荷手法です。閉じ込め効果や地面効果は考慮されません。

フォーマット

定義

コメント

1. モデリング状況はExp_dataフラグを使用して設定されます。ユーザーは、爆発データ（X_det、Y_det、Z_det）、爆発質量（W_TNT）、ターゲットサーフェス（surf_ID）、爆発時間（T_det）を指定します。他のすべてのパラメータとフラグにはデフォルト値があります。
   Exp_data=3の場合は、爆発高さを定義する必要があります。

   
   
   図 1.

   ユーザーサーフェス上の指定されたポイントで、対応する半径 $R$ と爆発質量W_TNTを使用して、衝風波動の特性値が決定されます（到着時間 $t_a$ 、最大圧力P_max、正の期間 $\text{Δ}{t}_{+}$ 、衝撃 $I_{+}$ など）。。入射波と反射波はどちらもFriedlanderの式に従います：

   • I_model = 1（Friedlanderモデル）の場合(1)
     $${\mathrm{P}}_{Friedlander}\left(t\right)=P_{}\cdot e^{\frac{-\left(t-t_a\right)}{\text{Δ}{t}_{+}}}\left(1-\frac{t-t_a}{\text{Δ}{t}_{+}}\right)$$
   • If I_model = 2（修正Friedlanderモデル）の場合(2)
     $${\mathrm{P}}_{Friedlander}\left(t\right)=P_{\max }\cdot e^{\frac{-b(t-t_a)}{\text{Δ}{t}_{+}}}\left(1-\frac{t-t_a}{\text{Δ}{t}_{+}}\right)$$

   ここで、 $P_{\max },\text{Δ}{t}_{+},t_a$ は、スケーリングされた指定の距離 $\frac{R}{W^{\raisebox{1ex}{$1$}\!\left/ \!\raisebox{-1ex}{$3$}\right.}}$ において実験的にわかっています。 ³

   修正Friedlanderモデル（I_model=2）では、‘b’は正のインパルスをフィッテングするために導入された減衰パラメータです。

   'b’は次のように解かれます：(3)

   $${\displaystyle \underset{t_a}{\overset{t_a+\text{Δ}{t}_{+}}{\int }}{\mathrm{P}}_{Friedlander}\left(t\right)dt=I_{+}}$$

   
   
   図 2. Friedlanderの式からの爆風プロファイル

2. このフィッティングされた時刻歴関数 ${\mathrm{P}}_{incident}\left(t\right)$ および ${\mathrm{P}}_{reflected}\left(t\right)$ は、与えられたフェイスの重心Z’（図 3）における爆風の載荷 ${\mathrm{P}}_{BLAST}\left(t\right)$ の計算にも用いられます。 ²(4)
   $${\mathrm{P}}_{BLAST}\left(t\right)=\left\{\begin{array}{c}{\cos }^2\theta \cdot {\mathrm{P}}_{reflected}\left(t\right)+\left(1+{\cos }^{}\theta -2{\cos }^2\theta \right)\cdot {\mathrm{P}}_{incident}\left(t\right)\text{ if }\cos \theta >0\\ {\mathrm{P}}_{incident}\left(t\right)\text{ if }\cos \theta \le 0\end{array}\right.$$

   
   
   図 3. フェイスの重心Z’に付与される爆風圧力. フェイスの向きに依存

   ここで、 $\theta$ はサーフェスのセグメント（重心Z’）と爆発位置への方向との間の角度です。

   これは、セグメントが直接爆発位置に面している場合、爆風圧力は反射圧力に等しく、セグメントが直接爆発位置に面していない場合は、入射圧力に等しいことを示しています。このモデル化は、到着時間と入射圧力が、関連する構造のシャドーイングで調整されていないため、簡易なものです。閉じ込め効果やトンネル効果も考慮されていません。

   また、サーフェスが外を向いた法線ベクトルを有していることが必要とされます。

3. I_z =1の場合、Rは一定で、Starter実行時にtime=0.00において計算されます。I_z =2の場合、 $R=R(t)$ は、Engine計算時のサイクルごとに更新されます。
4. W_TNTが設定されていない場合、質量はゼロで、関連するサーフェスに圧力は付与されません。
5. モデル化された爆発物がTNTではない場合、等価のTNT質量が与えられなければなりません。
6. 実験データは単位系{cm, g, $\mu s$ }を使用します。/BEGINで定義された単位は、実験データの単位をモデル単位に変換するために使用されます。したがって、/BEGINで定義された単位はモデルで使用される単位と精確に一致する必要があります。
7. $T=0$ から $t^{*}=\inf \left(T_{arrival}\right)$ までの計算時間をスキップすることが可能です。シフト値はStarter実行中に自動的に計算されます。 $t^{*}$ までの計算を無効にするには、I_tshiftの値を2にする必要があります。

   
   
   図 4. I_tshiftは最初の波の到着時間までの計算時間をスキップさせることが可能

8. $N_{dt}$ パラメータは、構造の時間ステップが十分大きくない場合に最小時間ステップを課すことが可能です。 $\text{Δ}{t}_{blast}=\frac{\inf \left(T_0\right)}{N_{dt}}$ を課すと、正のフェイズ中、すなわち爆風波の指数関数的減少の間に十分な時間ステップが確保されます。デフォルトでは、 $N_{dt}=100$ です。

   
   
   図 5.

9. Friedlander爆風モデルの正の部分を維持するために、パラメータ $P_{\min }$ が導入されました。(5)
   $${\mathrm{P}}_{}\left(t\right)=\max (P_{BLAST}(t),P_{\min })$$

   図 6. 過剰圧力（載荷の正の部分）のみを維持するために使用されるパラメータ $P_{\min }$