/MONVOL/FVMBAG2

• 個別の/MAT/GASカードで指定される気体材料
• 個別の/PROP/INJECT1または/PROP/INJECT2カードで指定される注入気体混合物とインジェクタープロパティ

フォーマット

定義

コメント

1. エアバッグの外部サーフェスは、4節点シェル要素と3節点シェル要素のみで構築されている必要があります。/SURF/SEGオプションでエアバッグの外部サーフェスを定義することはできません（/SURF/SEGでサブサーフェスを定義した場合は、/SURF/SURFで定義することもできません）。
2. 外部サーフェスは、法線が外向きの閉じたボリュームで構成される必要があります。
3. モデルの正しい単位を、/BEGINで定義するか、unit_IDで参照しているローカルの/UNITで定義する必要があります。気体定数、注入速度、事前定義の気体材料は、モデルで定義している単位に基づいて設定されます。
4. エアバッグ外側の空気の圧力と温度と、エアバッグ内部の空気の初期圧力と温度は、P_extおよびT₀に設定されます。
5. ベントホールを通じた排気：

   I_form = 1の場合、通気速度は、エアバッグ内の局所圧力を使用してBernoulli式によって計算されます。

   流出速度は次の式で与えられます。(1)

   $$u^2=\frac{2\gamma }{\gamma -1}\frac{P}{\rho }\left(1-{\left(\frac{P_{\mathit{ext}}}{P}\right)}^{\frac{\gamma -1}{\gamma }}\right)$$

   質量の流出速度は次の式によって定義されます。

   I_form = 2の場合、通気速度はChemkin式によって計算されます。(2)

   $$v=Fscal{e}_v\cdot f_v(P-P_{ext})$$

   ここで、 $f_v$ はfct_ID_vによって定義されます。

   I_form = 3の場合、通気速度は、ベントホールサーフェスに垂直な局所流体速度の成分に等しくなります。局所密度と局所エネルギーが、穴を通過する流出質量と流出エネルギーの計算に使用されます。

6. 気体注入をアクティブ化するセンサーが存在しない場合、時間TがT_startより大きくなるか、 $\text{Δ}t{P}_{def}$ で指定された時間よりも長い時間、圧力PがP_def値を超えると、ベントホールと空隙がアクティブになります。
7. 少なくとも1つのインジェクターがセンサーでアクティブになっている場合、排気オプションと空隙オプションのアクティブ化はI_ttfで制御されます。

   Tinjは、センサーによって最初のインジェクターがアクティブになる時間です。

   I_ttf = 0

   	排気、空隙
   アクティブ化	$P>\text{Δ}{P}_{def}$ の場合に $\text{Δ}t{P}_{def}$ よりも長い時間、または $T>T_{start}$
   非アクティブ化	Tstop
   時間依存関数	シフトなし

   I_ttf = 3

   	排気、空隙
   アクティブ化	$T>T_{inj}$ および $P>\text{Δ}{P}_{def}$ の場合、 $\text{Δ}t{P}_{def}$ よりも長い時間、または $T>T_{inj}+T_{start}$
   非アクティブ化	$T_{inj}+T_{stop}$
   時間依存関数	右記によってシフト； $T_{inj}+T_{start}$

   関連する他のすべての曲線は、対応する排気、空隙、または連結のオプションがアクティブな場合にアクティブになります。

   I_ttfの値の多様性は、これまでの経緯に起因しています。I_ttf =1および 2の値は廃止されているので、使用しないでください。通常使用する値は、I_ttf=0（シフトなし）またはI_ttf=3（関連するすべてのオプションがT_injだけシフト）です。

8. surf_ID_v ≠ 0（surf_ID_vが定義されている）の場合、ベントホール面積は次のように計算されます：(3)
   $$vent\_holes\_area=A_{vent}\cdot {\mathrm{f}}_A\left(\frac{A}{A_0}\right)\cdot {\mathrm{f}}_t\left(t\right)\cdot {\mathrm{f}}_P\left(P-P_{ext}\right)$$
   ここで、

   $A$

   サーフェスsurf_IDの面積_v

   $A_0$

   サーフェスsurf_IDの初期面積_v

   ${\mathrm{f}}_t$ 、 ${\mathrm{f}}_P$ 、 ${\mathrm{f}}_A$

   fct_ID_t、fct_ID_Pおよびfct_IDの関数_A

9. ベント閉鎖がアクティブになった場合、べントホールのサーフェスは次のように計算されます：(4)
   $$vent\_holes\_area=A_{vent}\cdot A_{non\_impacted}\cdot {\mathrm{f}}_t\left(t\right)\cdot {\mathrm{f}}_P\left(P-P_{ext}\right)\cdot {\mathrm{f}}_A\left(\frac{A_{non\_impacted}}{A_0}\right)$$
   (5)
   $$+B_{\mathit{vent}}\cdot A_{\mathit{impacted}}\cdot {\text{f}}_{t^{\prime }}\left(t\right)\cdot {\text{f}}_{P^{\prime }}\left(P-P_{\mathit{ext}}\right)\cdot {\text{f}}_{A^{\prime }}\left(\frac{A_{\mathit{impacted}}}{A_0}\right)$$
   ここで、衝撃を受けるサーフェスは：(6)

   $$A_{\mathit{impacted}}=\sum _{e\in S_{\mathit{vent}}}\frac{n_c\left(e\right)}{n\left(e\right)}{A}_e$$
   衝撃を受けないサーフェスは：(7)

   $$A_{\mathit{non}\_\mathit{impacted}}=\sum _{e\in S_{\mathit{vent}}}\left(1-\frac{n_c\left(e\right)}{n\left(e\right)}\right)A_e$$

   
   
   図 1.

   ここで、ベントホールsurf_ID_vの各要素eに対して、 $n_c\left(e\right)$ は、要素を定義している節点 $n\left(e\right)$ のうち衝撃を受ける節点の数を表しています。

   A₀はサーフェスsurf_IDの初期面積_v

   f_t、f_P、およびf_Aは、fct_ID_t、fct_ID_P、およびfct_IDの関数。_A

   f_t'、f_P'、およびf_A'は、fct_ID_t'、fct_ID_P'、およびfct_IDの関数。_A'

10. Chemkin排気の定式化（I_form=2）でsurf_ID_v = 0（surf_ID_vが定義されていない）の場合、Starterエラーが出力されてRadiossが終了します。
11. 関数fct_ID_tとfct_ID_Pを指定しなかった場合（識別子が0）、これらは1に等しいとみなされます。
12. 関数fct_ID_Aは、指定されなかった場合、1であるとみなされます。
13. ベントホールと多孔サーフェス領域の接触ブロックについては、対応するインターフェース（インターフェース/INTER/TYPE7または/INTER/TYPE23の３行目）でフラグI_BAGを1に設定する必要があります。そうしなかった場合は、衝撃でインターフェースに押し込まれる節点は、A_impactedとA_{non_impacted}に関する前の式の被衝撃節点と見なされません。
14. 空隙率定式化による漏れについて、質量の流出速度は次のように計算されます：
    • Iform_ps = 1 ${\dot{m}}_{\mathit{out}}=A_{\mathit{eff}}\sqrt{2P\rho }{Q}^{\frac{1}{\gamma }}\sqrt{\frac{\gamma }{\gamma -1}\left[1-Q^{\frac{\gamma -1}{\gamma }}\right]}$ (Isentropic - Wang Nefske)
    • Iform_ps = 2 ${\dot{m}}_{\mathit{out}}=A_{\mathit{eff}}\rho v(P-P_{\mathit{ext}})$

      ここで、vは気体の流出速度（Chemkin）

    • Iformps= 3 ${\dot{m}}_{\mathit{out}}=A_{\mathit{eff}}\sqrt{2\rho (P-P_{\mathit{ext}})}$ (Graefe)

    有効な排気面積A_effは、TYPE19またはTYPE58の繊維材料に関する/LEAK/MATの入力に従って計算されます。

15. 漏れのブロックをアクティブにすると（Iblockage=1）、有効な排気面積は次のように変更されます：(8)
    $$A_{\mathit{eff}}=A_{\mathit{non}\_\mathit{impacted}}$$

    $A_{\mathit{non}\_\mathit{impacted}}$ は衝撃を受けないサーフェスです。

    ブロックがアクティブになるのは、関係する接触インターフェース（インターフェースTYPE7、TYPE19、およびTYPE23の行3）でフラグIBAGを1に設定している場合のみです。

16. エアバッグの膨張プロセス中に有限体積がエラーになる（体積が負の内部質量になるか、エネルギーが負になる）場合は、計算を継続するために、その近傍のいずれかと結合されます。次の2つの結合アプローチが使用されます。
    • 全体結合： 有限体積は、その体積がすべての有限体積の平均体積に特定の係数をかけた値より小さくなった場合に結合されます。フラグI_gmergは、平均体積として現在の平均体積を使用するか（I_gmerg =1）、初期平均体積を使用するか（I_gmerg =2）を指定します。平均体積から最小体積を求めるための係数は、C_gmergです。
    • 時間ステップで決まる結合： 有限体積に対する時間ステップが、/DT/FVMBAGに定義した値未満の場合、その有限体積は隣接する有限体積と結合されます。
17. 失われた熱流量は次の式によって定義されます：(9)
    $$\dot{\text{Q}}\left(x,t\right)=H_{\mathit{conv}}\cdot \text{Area}\left(x,t\right)\cdot \left(\text{T}\left(x,t\right)-T_0\right)$$
18. ベントホールのサーフェスの要素（surf_ID_v）がインジェクター（surf_ID_inj）に属している場合、この要素はベントホールから無視されます。t = 0の時点で計算される定数の補正係数fがベントホールのサーフェス全体に適用されます：(10)
    $$f=\frac{S_{\mathit{vent}}}{S_{\mathit{vent}}-{\text{S}}_{\mathit{injector}}}$$
19. 多孔サーフェスの要素もインジェクタ（surf_ID_inj）に属している場合、この要素は多孔サーフェスから無視されます。
20. T_switchを均等圧力に切り替える時間は、始動時間に比例します。
21. オプションI_switch=2を使用すると、有限体積が1つしか残っていない場合でも、エアバッグは常に有限体積法で計算されます。気体パラメータは、単一の有限体積に切り替える前と後で同じです。一定圧力法（I_switch=1）に切り替えた場合、圧力または気体に多少の変動が見られる場合があります。
22. Pswitchは、エアバッグの平均圧力に対する有限体積圧力の平均偏差の比率です。(11)
    $$P_{switch}=\frac{\text{SD(FV pressure)}}{\text{Average pressure}}$$

    この比率は、/TH/MONVOL変数UPCRITを用いて出力することができます。P_switch は各有限体積内の圧力がエアバッグ内の平均圧力に近づくにつれてゼロに近づきます。