/MONVOL/COMMU1

ブロックフォーマットキーワード 注入気体の混合入力を使用する多室型エアバッグを記述します。このキーワードは、/MONVOL/COMMU(廃止)に似ていますが、より柔軟な入力が可能です。

  • 個別の/MAT/GASカードで指定されることが可能な気体材料
  • インジェクターは、インジェクター用の別個の/PROP/INJECT1で指定できます。
  • 時間または相対圧力に応じて、エアバッグ室間の連結領域をスケーリングできます。

フォーマット

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
/MONVOL/COMMU1/monvol_ID/unit_ID
monvol_title
surf_IDex   Hconv            
Ascalet AscaleP AscaleS AscaleA AscaleD
mat_ID   μμ Pext T0 Iequil Ittf
インジェクタの数
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
Njet                  
各インジェクタについて定義
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
inject_ID sens_ID Ijet node_ID1 node_ID2 node_ID3        
ジェット関数のデータ(Ijet = 1の場合にのみ読み出す)
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
fct_IDPt fct_IDPθfct_IDPθ fct_IDPδfct_IDPδ   FscalePt FscalePθFscalePθ FscalePδFscalePδ
NventベントホールとNporsurf多孔繊維サーフェスを定義
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
Nvent Nporsurf                
各ベントホールについて定義
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
surf_IDv Iform Avent Bvent     vent_title
Tstart Tstop ΔPdefΔPdef ΔtPdefΔtPdef   IdtPdef
fct_IDt fct_IDP fct_IDA   Fscalet FscaleP FscaleA
fct_IDt' fct_IDP' fct_IDA'   Fscalet' FscaleP' FscaleA'
各多孔サーフェスに対して挿入
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
surf_IDps Iformps Iblockage           surface_title
Tstart Tstop ΔPdefΔPdef ΔtPdefΔtPdef   IdtPdef
Iformps =0の場合にのみ挿入
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
Cps Areaps fct_IDcps fct_IDaps Fscalecps Fscaleaps
Chemkinモデルのデータ(Iform = 2またはIformps =2の場合のみ挿入)
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
fct_IDv   Fscalev            
連結エアバッグの数
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
Nbag                  
Nbag個の連結エアバッグを定義(連結エアバッグあたり2行)
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
bag_ID surf_IDc ΔPCdefΔPCdef Acom Tcom ΔtPCdefΔtPCdef
fct_IDCt fct_IDCP FscaleCt FscaleCP    

定義

フィールド 内容 SI単位の例
monvol_ID モニター体積識別子

(整数、最大10桁)

 
unit_ID 単位識別子

(整数、最大10桁)

 
monvol_title モニター体積のタイトル

(文字、最大100文字)

 
surf_IDex 外部サーフェス識別子

(整数)

 
Hconv 熱伝達係数 23

(実数)

[Wm2K][Wm2K]
Ascalet 時間ベースの関数に対する横軸のスケールファクター

デフォルト = 1.0(実数)

[s][s]
AscaleP 圧力ベースの関数に対する横軸のスケールファクター

デフォルト = 1.0(実数)

[Pa][Pa]
AscaleS 面積ベースの関数に対する横軸のスケールファクター

デフォルト = 1.0(実数)

[m2][m2]
AscaleA 角度ベースの関数に対する横軸のスケールファクター

デフォルト = 1.0(実数)

[rad][rad]
AscaleD 距離ベースの関数に対する横軸のスケールファクター

デフォルト = 1.0(実数)

[m][m]
mat_ID 初期気体の材料識別子(/MAT/GAS

(整数)

 
μμ 体積粘性

デフォルト = 0.01(実数)

 
Pext 外部圧力

(実数)

[Pa][Pa]
T0 初期温度

デフォルト = 295K(実数)

[K][K]
Iequil 初期熱力学的つり合いフラグ
= 0
エアバッグに初期充填されている気体の質量は、時間ゼロにおける体積に対して決定されます。
= 1
エアバッグに初期充填されている気体の質量は、注入開始時における体積に対して決定されます。

(整数)

 
Ittf 時間シフトフラグ
少なくとも1つの注入センサーが指定されている場合のみアクティブ。センサーで指定されているTTFの時間に注入が開始されたときの排気、空隙、および連結オプションの時間シフトを指定します。
= 0(デフォルト)
=1および2
廃止
=3
(すべてのオプションをシフト). 17

(整数)

 
Njet インジェクタの数

(整数)

 
inject_ID インジェクタプロパティ識別子

(整数)

 
sens_ID センサーの識別子

(整数)

 
Ijet ジェットフラグ
= 0
ジェット効果なし
= 1
ジェット効果あり

(整数)

 
node_ID1node_ID2node_ID3 ジェット形状を定義するための節点識別子N1N2およびN3

(整数)

 
fct_IDPt ΔP(t)ΔP(t) を定義する関数の識別子番号

(整数)

 
fct_IDPθfct_IDPθ ΔP(θ)ΔP(θ) を定義する関数の識別子番号

(整数)

 
fct_IDPδfct_IDPδ ΔP(δ)ΔP(δ) を定義する関数の識別子番号

(整数)

 
FscalePt fct_IDPtのスケールファクター。

デフォルト = 1.0(実数)

[Pa][Pa]
FscalePθFscalePθ fct_IDPθfct_IDPθ のスケールファクター。

デフォルト = 1.0(実数)

[Pa][Pa]
FscalePδFscalePδ fct_IDPδfct_IDPδ のスケールファクター。

デフォルト = 1.0(実数)

[Pa][Pa]
Nvent ベントホールの数

(整数)

 
Nporsurf 多孔サーフェスの数

(整数)

 
surf_IDv ベントホール領域サーフェス識別子

(整数)

 
Iform 定式化フラグ
= 0(デフォルト)
1に設定されます。
= 1
等エンタルピー(デフォルト)
= 2
Chemkin
= 3
等エンタルピーで気体(mat_ID)流入の可能性あり

(整数)

 
Avent surf_IDv≠ 0の場合: ベントホール面積に対するスケールファクター

デフォルト = 1.0(実数)

 
surf_IDv = 0の場合: ベントホール面積

デフォルト = 0.0(実数)

[m2][m2]
Bvent surf_IDv ≠ 0の場合: 衝撃を受けるベントホール面積に対するスケールファクター

デフォルト = 1.0(実数)

 
surf_IDv = 0の場合: Bvent ベントホール面積に対し0にリセットされます。

デフォルト = 0.0(実数)

[m2][m2]
vent_title ベントホールのタイトル

(文字、最大20文字)

 
Tstart 排気の開始時間

デフォルト = 0(実数)

[s][s]
Tstop 排気の停止時間

デフォルト = 1030(実数)

[s][s]
ΔPdefΔPdef ベントホール膜を開口する圧力差

ΔPdef=PdefPextΔPdef=PdefPext

デフォルト = 0(実数)

[Pa][Pa]
ΔtPdefΔtPdef 圧力がPdefを超えてベントホールが開口するまでの最短時間

デフォルト = 0(実数)

[s][s]
IdtPdef ΔPdefΔPdef に到達した場合の時間遅延フラグ:
= 0
排気をアクティブにするには、 ΔtPdefΔtPdef の累積時間の間、圧力が ΔPdefΔPdef を超える必要があります。
= 1
排気は、 ΔPdefΔPdef に到達した後、 ΔtPdefΔtPdef が経過すると、アクティブになります。

(整数)

 
fct_IDt 空隙率対時間関数識別子

(整数)

 
fct_IDP 空隙率対圧力関数識別子

(整数)

 
fct_IDA 空隙率対面積関数識別子

(整数)

 
Fscalet fct_IDtのスケールファクター。

デフォルト = 1.0(実数)

 
FscaleP fct_IDPのスケールファクター。

デフォルト = 1.0(実数)

 
FscaleA fct_IDAのスケールファクター。

デフォルト = 1.0(実数)

 
fct_IDt' 接触時の空隙率対時間関数識別子

(整数)

 
fct_IDP' 接触時の空隙率対圧力関数識別子

(整数)

 
fct_IDA' 接触時の空隙率対被衝撃サーフェス関数識別子

(整数)

 
Fscalet' fct_IDt'のスケールファクター。

デフォルト = 1.0(実数)

 
FscaleP' fct_IDP'のスケールファクター。

デフォルト = 1.0(実数)

 
FscaleA' fct_IDA'のスケールファクター。

デフォルト = 1.0(実数)

 
surf_IDps 多孔サーフェス識別子(Iformps =0の場合は無視)

(整数)

 
Iformps 空隙率定式化
= 0
Bernouilli(Wang & Nefske)(繊維材料への依存なし)
= 1
Bernouilli(Wang & Nefske)
= 2
Chemkin
= 3
Graefe

(整数)

 
Iblockage 接触時の漏れをブロックするフラグ(Iformps > 0)
= 0
なし
= 1
あり

(整数)

 
surface_title 多孔サーフェスのタイトル

(文字、最大20文字)

 
Cps 漏れ面積に対するスケールファクター(Iformps =0)

(実数)

 
Areaps 漏れ面積(Iformps =0)

(実数)

[m2][m2]
fct_IDcps Cps(t)を定義する関数識別子で、Cpsがゼロでない場合は無視

(整数)

 
fct_IDaps Areaps(P-Pext)を定義する関数識別子で、Areapsがゼロでない場合は無視

(整数)

 
Fscalecps fct_IDcpsのスケールファクター。

デフォルト = 1(実数)

 
Fscaleaps fct_IDのスケールファクターaps

デフォルト = 1(実数)

[m2][m2]
fct_IDv 流出速度関数識別子、Chemkinモデル

(整数)

 
Fscalev fct_IDvに対するスケールファクター

デフォルト = 1.0(実数)

[ms][ms]
Nbag 連結エアバッグの数

(整数)

 
bag_ID エアバッグ識別子 20

(整数)

 
surf_IDC 連結サーフェス識別子

(整数)

 
ΔPCdefΔPCdef 連結サーフェス膜を開口する圧力差

(実数)

[Pa][Pa]
Acom surf_IDC = 0の場合、連結サーフェス [m2][m2]
surf_IDC ≠ 0の場合、サーフェスに対するスケールファクター

デフォルト = 1.0(実数) 23

Tcom 連結の開始時間

(実数)

[s][s]
ΔtPCdefΔtPCdef 圧力差が ΔPCdefΔPCdef を超えて連結サーフェス膜を開口するまでの最短時間

(実数)

[s][s]
fct_IDCt 連結サーフェス対時間関数識別子

(整数)

 
fct_IDCP 連結サーフェス対相対圧力関数識別子

(整数)

 
FscaleCt fct_IDCtのスケールファクター。

デフォルト = 1.0(実数)

 
FscaleCP fct_IDCPのスケールファクター。

デフォルト = 1.0(実数)

 

コメント

  1. エアバッグの外部サーフェスは、4節点シェル要素と3節点シェル要素のみで構築されている必要があります。/SURF/SEGでエアバッグの外部サーフェスを定義することはできません(/SURF/SEGでサブサーフェスを定義した場合は、/SURF/SURFで定義することもできません)。
  2. 体積は、閉じていて法線が外向きである必要があります。
  3. 横軸のスケールファクターは、エアバッグ関数の横軸の単位を変換するために使用されます。たとえば:(1)
    F(t)=ft(tAscalet)
    ここで、
    t
    時間
    ft
    fct_IDの関数t
    (2)
    F(P)=fP(PAscaleP)
    ここで、
    P
    圧力
    fP
    fct_IDの関数P
    これらのオプションは廃止されています。通常は、代わりに曲線のスケーリングパラメータを使用します。
  4. エアバッグ外側の空気の圧力と温度と、エアバッグ内部の空気の初期圧力と温度は、PextおよびT0に設定されます。
  5. エアバッグに初期充填されている気体の特性(温度と圧力)にはデフォルト値がないため、定義する必要があります。また、その特性は各伝達エアバッグで同じでなければなりません。
  6. 時間ゼロ(Iequil =0)または噴射開始時(Iequil =1)の初期熱力学的つり合いを記述します。初期熱力学的つり合いは、時間ゼロの体積または噴射開始時の体積に対して、次の式に基づいて求めます:(3)
    PextV=RM0MiT0
    ここで、 Mi 0はエアバッグに初期充填されている気体の質量、 M0 はエアバッグに初期充填されている気体のモル質量、 R /BEGINカードで指定される単位系に依存する気体定数です。例えば、SI単位系では次のようになります:(4)
    R=8.314JmoleK
  7. ジェットを使用した場合、エアバッグの各要素に追加の ΔPjet 圧力がかかります:(5)
    ΔPjet=ΔP(t)ΔP(θ)ΔP(δ)max(n·m,0)

    関数ΔP(t) はセンサーで指定される時間で自動的にシフトし、注入をアクティブ化します。

  8. ここで、mは要素の中心をセグメント(N1N3)上に投影したポイントから要素の中心までの正規化ベクトル、 θ はベクトルMN2mの間の角度(単位は度)、 δ は要素の中心から、要素の中心をセグメント(N1N3)上に投影したポイントまでの距離です。
    ポイントのセグメント(N1N3)上への投影は、ポイントをMN2の方向に向けてライン(N1N3)上へ投影したものとして定義されます。この場合、投影先がセグメント(N1N3)に含まれる必要があります。投影先がセグメントに含まれない場合、ポイントのセグメント(N1N3)上への投影は、N1またはN3のうち、近い方の節点として定義されます。

    clip0087
    図 1. ジェット機の上反角

    ここで、MN1N3の間

  9. node_ID3 = 0の場合、node_ID3node_ID1に設定され、二面体は円錐形に減少します。
  10. fct_IDv = 0の場合: 等エンタルピー流出が想定されます。0でない場合はChemkinモデルが使用され、流出速度は次の式で求められます:(6)
    v=Fscalevfv(PPext)
    ここで、 fv fct_IDvの関数です。
    • 等エンタルピーモデル

      放出、すなわちモニター体積からの気体の排出は等エンタルピーであると想定されます。

      また、流動は衝撃がなく、大容量の容器から流出し、有効表面積Aの狭い開口部を通ると想定されます。

      エンタルピー保存則から、ベントホールにおける速度uが導出されます。したがって、Bernouilli式は次のように記述されます。

      (モニター体積) γγ1Pρ=γγ1Pextρvent+u22 (ベントホール)

      断熱条件を適用すると、次の式になります。

      (モニター体積) Pργ=Pextρventγ (ベントホール)

      ここで、 P はエアバッグに流入する気体の圧力、 ρ はエアバッグに流入する気体の密度です。

      したがって、流出速度は次の式によって定義されます:(7)
      u2=2γγ1Pρ(1(PextP)γ1γ)

      超音速流れでは、流出速度はSupersonic Outlet Flow理論マニュアルの記述で決められます。

      質量の流出速度は次の式によって定義されます。(8)
      ˙mout=ρvent*vent_holes_surface*u=ρ(PextP)1γ*vent_holes_surface*u
      エネルギーの流出速度は次の式によって定義されます:(9)
      ˙Eout=˙moutEρV=(PextP)1γ*vent_holes_surface*uEV

      ここで、 V はエアバッグの体積、 E はエアバッグに流入する気体の内部エネルギーです。

    • Chemkinモデル(10)
      ˙mout=ρvent_holes_surfacefv(PPext)Fscalev

      ここで、 ρ はエアバッグ内部の気体の密度、 fv fct_IDの関数v

  11. ベントホールの面積は、次のように計算されます:(11)
    vent_holes_area =AventAnon_impactedft(t)fP(PPext)fA(Anon_impactedA0)
    +BventAimpactedft(t)fP(PPext)fA(AimpactedA0)
    ここで、衝撃を受けるサーフェスは:(12)
    Aimpacted=eSventnc(e)n(e)Ae
    衝撃を受けないサーフェスは:(13)
    Anon_impacted=eSvent(1nc(e)n(e))Ae

    Image12
    図 2.

    ここで、エアバッグ材料の各要素eに対して、 nc(e) は、要素を定義している n(e) 節点のうち衝撃を受ける節点の数、 Ae は要素eの面積

    また、

    A0はサーフェスsurf_IDの初期面積v

    ft fP およびfAfct_IDtfct_IDPおよびfct_IDの関数A

    ft'fP'、およびfA'は、fct_IDtfct_IDP、およびfct_IDAの関数。

  12. 関数fct_IDt''fct_IDP'を指定しなかった場合(識別子が0)、これらは1に等しいとみなされます。
  13. 関数fct_IDA'は次のようにみなされます:

    fA(A)=A これが指定されなかった場合

  14. ベントホールと多孔サーフェス領域の接触ブロックについては、対応するインターフェース(インターフェースTYPE7またはTYPE23の3行目)でフラグIBAGを1に設定する必要があります。そうしなかった場合は、衝撃でインターフェースに押し込まれる節点は、AimpactedAnon_impactedに関する前の式の被衝撃節点と見なされません。
  15. ベントホールとチャンバー間コンポーネントは、エアバッグ(チャンバー)の外部サーフェスに含まれている必要があります。
  16. 気体注入をアクティブ化するセンサーが存在しない場合、時間TTstartより大きくなるか、 ΔtPdef で指定された時間よりも長い時間、圧力PPdef値を超えると、ベントホールと空隙がアクティブになります。
  17. 少なくとも1つのインジェクターがセンサーでアクティブになっている場合、排気、空隙、および連結オプションのアクティブ化はIttfで制御されます。

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

    Iinj = 0:
      排気、空隙 連結
    アクティブ化 P>ΔPdef の場合に ΔtPdef よりも長い時間、または T>Tstart P>ΔPCdef の場合に ΔtPCdef よりも長い時間、または T>Tcom
    非アクティブ化 Tstop なし
    時間依存関数 シフトなし シフトなし
    Ittf = 3:
      排気、空隙 連結
    アクティブ化 T>Tinj および P>ΔPdef の場合、 ΔtPdef よりも長い時間、または T>Tinj+Tstart T>Tinj および P>ΔPCdef の場合、 ΔtPCdef よりも長い時間、または T>Tinj+Tcom
    非アクティブ化 Tinj+Tstop なし
    時間依存関数 右記によってシフト; Tinj+Tstart 右記によってシフト; Tinj+Tcom

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

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

  18. 空隙率定式化による漏れについて、質量の流出速度は次のように計算されます:
    • Iformps = 0 ˙mout=Aeff2PρQ1γγγ1[1Qγ1γ] (Isentropic - Wang Nefske)
      ここで、(14)
      Q=PextP

      また、 Aeff=CpsAreaps または Aeff=Cps(t)Areaps(PPext)

      有効な排気面積Aeffがエアバッグの繊維材料によって変化することはありません。

    • Iformps > 0の場合、有効な排気面積Aeffは、TYPE19またはTYPE58の繊維材料に関する/LEAK/MATの入力に従って計算されます。
    • Iformps = 1 ˙mout=Aeff2PρQ1γγγ1[1Qγ1γ] (Isentropic - Wang Nefske)
    • Iformps = 2 ˙mout=Aeffρv(PPext)

      ここで、vは気体の流出速度(Chemkin)

    • Iformps= 3 ˙mout=Aeff2ρ(PPext) (Graefe)
  19. 漏れのブロックをアクティブにすると(Iblockage=1)、有効な排気面積は次のように変更されます:(15)
    Aeff=Anon_impacted

    Anon_impactedは衝撃を受けないサーフェスです。 11

    ブロックは、関係する接触インターフェース(インターフェースTYPE7およびTYPE23の行3)において、フラグIBAGが1に設定された場合にのみアクティブになります。

  20. 1つの多室型エアバッグにおいて、/MONVOL/COMMUカードと/MONVOL/COMMU1カードを組み合わせることはできません。ただし、同じモデル内で、エアバッグごとに異なる多室型エアバッグを使用することはできます(/MONVOL/COMMUまたは/MONVOL/COMMU1に基づく)。
  21. 気体の注入をアクティブ化するセンサーが存在しない場合、T>Tcomとなるか、 ΔtPCdef よりも長い時間、圧力が ΔPCdef を超えると、連結サーフェスが開きます。 17
  22. 連結サーフェスScomは、次のように計算されます:
    • surf_ID=0の場合、(16)
      Scom=AcomfCt(t)fCP(ΔP)
    • surf_ID > 0でAreasurf_IDのサーフェスの場合、(17)
      Scom=AcomAreafCt(t)fCP(ΔP)

    ここで、 ΔP はチャンバーと fCt の圧力差、 fCP fct_IDCtおよびfct_IDの関数CP

  23. 失われた熱流量は次の式によって定義されます:(18)
    ˙Q(x,t)=HconvArea(x,t)(T(x,t)T0)