/MONVOL/FVMBAG (廃止)

ブロックフォーマットキーワード FVMBAGタイプのエアバッグを記述します。入力はAIRBAGタイプと同様です。

フォーマット

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
/MONVOL/FVMBAG/monvol_ID/unit_ID
monvol_title
surf_IDex                  
Ascalet AscaleP AscaleS AscaleA AscaleD
        Pext T0 Iequi Ittf
γ i cpai cpbi cpci    
インジェクタの数
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
Njet                  
Njet個のインジェクターを定義(インジェクターあたり3行)
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
γ cpa cpb cpc    
fct_IDmas Iflow Fscalemas fct_IDT FscaleT sens_ID    
Isjet                  
fct_IDvel   Fscalevel            
ベントホールの数
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
Nvent                  
Nvent個のベントホールを定義(ベントホールあたり4行)
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
surf_IDv Avent Bvent     Itvent    
Tvent Δ Pdef Δ tPdef fct_IDV FscaleV 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)
Vx3 Vy3 Vz3        
Vx1 Vy1 Vz1        
X0 Y0 Z0        
L1 L2 L3        
Nb1 Nb2 Nb3 grbrc_ID surf_IDin Iref        
その他のFVMBAGパラメータ
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
Igmerg   Cgmerg Cnmerg Ptole    
qa qb Hmin        
Ilvout Nlayer Nfacmax Nppmax Ifvani          

定義

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

(整数、最大10桁)

 
unit_ID 単位識別子

(整数、最大10桁)

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

(文字、最大100文字)

 
surf_IDex 外部サーフェス識別子 1

(整数)

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

デフォルト = 1.0(実数)

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

デフォルト = 1.0(実数)

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

デフォルト = 1.0(実数)

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

デフォルト = 1.0(実数)

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

デフォルト = 1.0(実数)

[ m ]
Pext 外圧

(実数)

[ Pa ]
T0 初期温度。

デフォルト = 295(実数)

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

(整数)

 
Ittf 排気時間シフトフラグ。注入センサーが指定されている場合のみアクティブ。
= 0または1
時間依存多孔性曲線は、注入センサーのアクティブ化時間ではシフトされません。TventTstopは無視されます。
= 2
時間依存多孔性曲線はTinjTinjは、第1インジェクターがセンサーでアクティブになる時間として定義されます)でシフトされます。
TventTstopは無視されます。
= 3
時間依存多孔性曲線はTinj +Tventでシフトされます。Tstopが指定されている場合、排気はTinj + Tstopで停止します。
 
γ i 初期温度における比熱の比 5

γ i = Cp i / Cv i

(実数)

 
cpai cpa 係数(関係式cpi(T)の)

(実数)

[ J kgK ] MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeaaciGaaiaabeqaamaabaabaaGcbaWaamWaaeaada WcaaqaaiaabQeaaeaacaqGRbGaae4zaiabgwSixlaabUeaaaaacaGL BbGaayzxaaaaaa@3DB3@
cpbi cpb 係数(関係式cpi(T)の)

(実数)

[ J kg K 2 ]
cpci cpc 係数(関係式cpi(T)の)

(実数)

[ J kg K 3 ]
Njet インジェクタの数

(整数)

 
γ 比熱の比

γ = C p / C v

(実数)

 
cpa cpa 係数(関係式cp(T)の)

(実数)

[ J kgK ] MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeaaciGaaiaabeqaamaabaabaaGcbaWaamWaaeaada WcaaqaaiaabQeaaeaacaqGRbGaae4zaiabgwSixlaabUeaaaaacaGL BbGaayzxaaaaaa@3DB3@
cpb cpb 係数(関係式cp(T)の)

(実数)

[ J kg K 2 ]
cpc cpc 係数(関係式cp(T)の)

(実数)

[ J kg K 3 ]
fct_IDmas 注入気体の質量対時間関数の識別子

(整数)

 
Iflow 質量対時間関数の入力タイプフラグ
= 0
質量が入力されます
= 1
質量流量が入力されます。

(整数)

 
Fscalemas 質量関数に対するスケールファクター

デフォルト = 1.0(実数)

[ kg ] または [ kg s ]
fct_IDT 注入気体の温度対時間識別子

(整数)

 
FscaleT 温度スケールファクター

デフォルト = 1.0(実数)

[ K ]
sens_ID 注入を開始するセンサー識別子

(整数)

 
Isjet インジェクタサーフェス識別子(インジェクタごとに異なる必要があります)

(整数)

 
fct_IDvel 注入気体速度識別子

(整数)

 
Fscalevel 注入気体スケールファクター

デフォルト = 1.0(実数)

[ m s ] MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaamWaaeaada Wcaaqaaiaab2gaaeaacaqGZbaaaaGaay5waiaaw2faaaaa@39DE@
Nvent ベントホールの数

(整数)

 
surf_IDv ベントホール膜サーフェスまたは多孔サーフェス識別子

(整数)

 
Avent surf_IDv0の場合: サーフェスに対するスケールファクター

デフォルト = 1.0

surf_IDv = 0の場合: ベントホールのサーフェス

デフォルト = 0.0(実数)

[ m 2 ] surf_IDV = 0の場合
Bvent surf_IDv0の場合: 衝撃を受けるサーフェスに対するスケールファクター

デフォルト = 1.0

surf_IDv = 0の場合: Bvent は0にリセットされます。

デフォルト = 0.0(実数)

[ m 2 ] surf_IDV = 0の場合
Itvent 排気の定式化 7
= 1
通気速度は、ベントホールサーフェスに垂直な局所流体速度の成分に等しくなります。局所密度と局所エネルギーが、穴を通過する流出質量と流出エネルギーの計算に使用されます。
= 2(デフォルト)
通気速度は、エアバッグ内の局所圧力を使用してBernoulli式によって計算されます。局所密度と局所エネルギーは、流出質量と流出エネルギーを計算するために使用されます。
= 3
通気速度はChemkin式によって計算されます。

(整数)

 
Tvent 排気の開始時間

デフォルト = 0.0(実数)

[ s ]
Δ P d e f ベントホール膜を開口する圧力差( Δ P d e f = Pdef - Pext

(実数)

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

(実数)

[ s ]
fct_IDV 流出速度関数識別子

(整数)

 
FscaleV 右記のスケールファクター; fct_IDV

デフォルト = 1.0(実数)

[ m s ] MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaamWaaeaada Wcaaqaaiaab2gaaeaacaqGZbaaaaGaay5waiaaw2faaaaa@39DE@
IdtPdef Δ P d e f に到達した場合の時間遅延フラグ:
= 0
排気をアクティブにするには、 Δt P d e f の累積時間の間、圧力が Δ P d e f を超える必要があります。
= 1
排気は、 Δ P d e f に到達した後、 Δt P d e f が経過すると、アクティブになります。
 
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(実数)

 
Vx3 X成分(ベクトルV3の(全体フレーム内))

(実数)

 
Vy3 Y 成分(ベクトルV3の(全体フレーム内))

(実数)

 
Vz3 Z 成分(ベクトルV3の(全体フレーム内))

(実数)

 
Vx1 X 成分(ベクトルV1の(全体フレーム内))

(実数)

 
Vy1 Y 成分(ベクトルV1の(全体フレーム内))

(実数)

 
Vz1 Z 成分(ベクトルV1の(全体フレーム内))

(実数)

 
X0 X 座標(ローカル原点Oの(全体フレーム内))

(実数)

 
Y0 Y 座標(ローカル原点Oの(全体フレーム内))

(実数)

 
Z0 Z 座標(ローカル原点Oの(全体フレーム内))

(実数)

 
L1 長さL1

(実数)

[ m ]
L2 長さL2

(実数)

[ m ]
L3 長さL3

(実数)

[ m ]
Nb1 方向1における有限体積の数

デフォルト = 1(整数)

 
Nb2 方向2における有限体積の数

デフォルト = 1(整数)

 
Nb3 方向3における有限体積の数

デフォルト = 1(整数)

 
grbric_ID ユーザー定義のソリッドグループの識別子

(整数)

 
surf_IDin 内部サーフェス識別子 25

(整数)

 
Iref 参照ジオメトリ上に自動FVMメッシュを適用するためのフラグ 24
= 0(デフォルト)
×
= 1

(整数)

 
Igmerg 全体結合定式化のフラグ 19

デフォルト = 1(整数)

 
Cgmerg 全体結合の係数 19

デフォルト = 0.02(実数)

 
Cnmerg 近傍結合の係数 19

(実数)

 
Ptole 有限体積識別用の許容値

デフォルト = 10-5(実数)

 
qa 2次体積粘性

デフォルト = 0.0(実数)

 
qb 線形体積粘性

デフォルト = 0.0(実数)

 
Hmin 浸透性を有する3角形の最小の高さ 21

(実数)

[ m ]
Ilvout 出力レベル: 0または1

デフォルト = 0(整数)

 
Nlayer V3方向に折りたたまれているエアバッグの層数の推定値(コメント23) 22

デフォルト = 10(整数)

 
Nfacmax 最初の自動メッシングステップにおける有限体積の影響を受けるエアバッグセグメントの推定最大数

デフォルト = 20(整数)

 
Nppmax 多角形の頂点の推定最大数

デフォルト = 20(整数)

 
Ifvani 有限体積のRadioss StarterアニメーションA000ファイルへの書き込みフラグ
= 0
×
= 1

(整数)

 

コメント

  1. surf_IDex は、3節点または4節点シェル要素(ボイド要素の場合もあります)に関連付けられたセグメントを使用して定義する必要があります。
  2. 体積は、閉じていて法線が外向きである必要があります。
  3. 横軸のスケールファクターは、エアバッグ関数の横軸の単位を変換するために使用されます。たとえば: (1)
    F ( t ' ) = fct _ ID ( t Ascale t )
    ここで、tは時間です。 (2)
    F ( p ' ) = fct _ ID ( p Ascale p )

    ここで、pは圧力です。

  4. 初期圧力はPextに設定されます。
  5. γ i = 0の場合、エアバッグに初期充填されている気体の特性は、第1インジェクターによって供給される気体の特性に設定されます。
  6. FVMBAG内の気体流は有限体積を使用して求められます。

    これらの有限体積の一部は、エアバッグ内部に配置され、特定のパートまたは内部体積全体を満たしているソリッドのグループを通してユーザーが入力できます。ユーザー定義ソリッドによって離散化されていない内部体積部分がまだ存在する場合、残りの体積は自動メッシュ生成プロシージャによって生成されます。この機能を使用して、例えばキャニスターをモデル化できます。

    有限体積は、三角端面のセット内に存在します。これらの頂点は、必ずしもエアバッグの節点と一致する必要はありません。エアバッグの外皮をモデル化する場合、4節点膜または3節点膜を使用できますが、3節点膜の使用を推奨します。

    monvol_airbag-env
    図 1.

    monvol_airbag2
    図 2.
  7. 流出速度は次の式で与えられます。
    (3)
    u 2 = 2 γ γ 1 P ρ ( 1 ( P ext P ) γ 1 γ )
    質量の流出速度は次の式によって定義されます。(4)
    m ˙ out = ρ v * vent _ holes _ surface * u
    エネルギーの流出速度は次の式によって定義されます:(5)
    E ˙ out = m ˙ out E ρ V = ( P ext P ) 1 γ * vent _ holes _ surface * u E V
    通気速度は、次の式で計算されます:(6)
    m out = ρ * vent _ holes _ surface * fct _ ID V * Fscale V ( P P ext )
  8. T > Tventの場合、または Δt P d e f よりも長い時間、圧力がPdefを超えた場合は、ベントホール膜は収縮します。
  9. surf_IDv0surf_IDvが定義されている)場合。 (7)
    vent _ holes _ surface = A vent * fct _ ID A ( A / A 0 ) * fct _ ID t ( t ) * fct _ ID P ( P P ext )

    ここで、A歯サーフェスsurf_IDの面積、A0はサーフェスsurf_IDvの初期面積です。

  10. surf_IDv = 0の場合(surf_IDvは未定義)、ベントホールは無視されます。(8)
    vent _ holes _ surface = A vent * fct _ ID t ( t ) * fct _ ID P ( P P ext )
  11. 関数fct_IDt'fct_IDPを指定しなかった場合(識別子が0)、これらは1に等しいとみなされます。
  12. 関数fct_IDAを指定しなかった場合、fct_IDA(A/A0) = 1であるとみなされます。
  13. ベントホールのサーフェスは、次のように計算されます。(9)
    vent _ holes _ surface =A vent * A non _ impacted * fct _ ID t ( A non _ impacted / A 0 ) * fct _ ID P ( P P ext ) + B vent * A impacted * fct _ ID t ( A impacted / A 0 ) * fct _ ID P ( P P ext )
    ここで、衝撃を受けるサーフェスは:(10)
    A impacted = e S vent n c ( e ) n ( e ) A e
    衝撃を受けないサーフェスは:(11)
    A non _ impacted = e S vent ( 1 n c ( e ) n ( e ) ) A e

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


    Image12
    図 3. 衝撃を受けるサーフェス / 衝撃を受けないサーフェスと接触する節点
  14. 関数fct_IDt'fct_IDPを指定しなかった場合(識別子が0)、これらは1に等しいとみなされます。
  15. 接触中に空隙率を使用するためには、関係するインターフェース(インターフェースタイプ5および7の行3)において、フラグIBAG1に設定する必要があります。そうしなかった場合は、衝撃でインターフェースに押し込まれる節点は、AimpactedAnon_impactedに関する前の式の被衝撃節点と見なされません。
  16. 自動有限体積メッシングパラメータ

    monvol_finite_vol
    図 4.
  17. 有限体積は2つのステップで生成されます。
    • 第1ステップでは、エアバッグの外皮上に排他的に配置された頂点が生成されます。これにより、外皮の変形に合わせて有限体積を更新したり、以下の手順に対応したりできます(分かりやすくするために2次元で表示):

      monvol_step1
      図 5.

      この手順では、方向V3(指定された切断方向)と方向V1の入力が必要です。切断方向に垂直な平面内の2つ目の方向V2が計算されます。V1V2の両方向で有限体積を配置して、切断幅を決定するためには、原点Oだけでなく、原点から正と負の両方向でカウントされた長さLiとステップ数Niも指定する必要があります。カッティング幅は次のように与えられます:Wi = 2Li / Ni

      原点Oおよび長さLi(Oからの正方向と負方向の両方を含む)によって水平面(V3に対して垂直)内に描かれるボックスには、水平面に投影されたメッシュにその体積の外皮の境界ボックスが含まれている必要があります。これは、この体積全体が有限体積に分割されることを保証するために必要です。

    • 第2ステップでは、有限体積の水平切断が実行されますが、堅く折りたたまれたエアバッグの場合はほとんど役に立たない可能性があります。この作業は、特に、エアバッグを広げる前に注入気体によって満たされたキャニスター内で注入を行う場合に必要です。
      この第2ステップでは、エアバッグ内部に配置された頂点が生成される場合があります。エアバッグの膨張に合わせて頂点を移動するために、エアバッグの外皮上に配置された2つの頂点間の垂直セグメント(方向V3と並行)にそれぞれの頂点が取り付けられます(図 6)。その参照セグメント内部の頂点の局所座標は、膨張プロセスを通して一定に保たれます。

      monvol_fvmbag
      図 6.

      水平切断幅は次の式で与えられます:W3 =2L3 /N3原点Oおよび長さL3(正方向と負方向の両方を含む)によってV3方向に与えられるセグメントには、V3方向のメッシュ投影に体積の外皮の境界ボックスが含まれる必要はありません。これは、2番目のステップでは既存の有限体積のみがカットされるためです。

  18. 自動メッシングに実際に使用されるベクトルV1は、ベクトルV3に対して入力ベクトルが直交化されたあとで取得されます。
  19. エアバッグの膨張プロセス中に有限体積がエラーになる(体積が負の内部質量になるか、エネルギーが負になる)場合は、計算を継続するために、その近傍のいずれかと結合されます。次の2つの結合アプローチが使用されます。
    • 全体結合: 有限体積は、その体積がすべての有限体積の平均体積に特定の係数をかけた値より小さくなった場合に結合されます。Igmergフラグは、平均体積として現在の平均体積を使用するか(Igmerg =1)、初期平均体積を使用するか(Igmerg =2)を指定します。平均体積から最小体積を求めるための係数は、Cgmergです。
    • 近傍結合: 有限体積は、その体積が近傍の平均体積に特定の係数をかけた値より小さくなった場合に結合されます。平均体積から最小体積を求めるための係数は、Cnmergです。
  20. CgmergCnmerg の両方が0でない場合は、両方の結合アプローチが同時に使用されることを意味します。強い衝撃を与える場合は、qa = 1.1およびqb = 0.05を設定することを推奨します。
  21. 2層の布が物理的に接触している場合は、有限体積間に流れが発生しないはずです。インターフェースギャップが原因で数値的にはこのようになりません。Hminは、浸透性を有する三角面の最小の高さを表しています。この値よりも小さい三角面は、不浸透性を有します。その値は、エアバッグのセルフ接触インターフェースのギャップに近づける必要があります。
  22. NlayerNfacmaxNppmaxは、有限体積作成プロセスを支援するメモリパラメータです。これらの値を変更しても、計算が停止することはありません。値を増やすと、自動メッシングのためのメモリ量が増え、計算時間が短縮されます。
  23. 有限体積作成プロセスでは、最初に、平面多角形が作成され、その後で、閉じた多面体に組み立てられ、三角端面に分解されます。 Nppmax は、これらの多角形に含まれる頂点の最大数です。
  24. Iref = 1の設定はモニター体積タイプFVMBAG または FVMBAG1に対して/REFSTAに基づいた参照ジオメトリにのみ働きます(/XREFに基づいた参照ジオメトリはまだサポートされていません)。
  25. 内部サーフェスまたはバッフルをモニター体積内の気体流れの障害物として考慮に用いるには、有限体積法のみ可能です。内部サーフェスは、モニター体積の外部および内部サーフェス節点と一致する節点があるTETRA4HEXAおよびPENTAの場合もあります)のように、モニター体積がソリッド要素で満たされている場合のみ、FVMで考慮されます。(これらのソリッドはgrbrick_IDで定義する必要があります)。材料モデルがLAW19の内部サーフェス繊維にのみ、空隙率を、0(空隙率なし)から1(空隙率完全)(ベント)の範囲で適用できます。気体流の方向が内部サーフェスの法線方向と反対の場合は、内部サーフェス上でインジェクターサーフェスも定義できます。
  26. ベントホールのサーフェスの要素(surf_IDv)がインジェクター(surf_IDinj)に属している場合、この要素はベントホールから無視されます。t = 0の時点で計算される定数cの補正係数f がベントホールのサーフェス全体に適用されます(12)
    f = S vent / ( S vent S injector )