/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 | ||||||
cpai | cpbi | cpci |
(1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | (10) |
---|---|---|---|---|---|---|---|---|---|
Njet |
(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 |
(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 |
(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(実数) |
|
AscaleP | 圧力ベースの関数に対する横軸のスケールファクター デフォルト = 1.0(実数) |
|
AscaleS | 面積ベースの関数に対する横軸のスケールファクター デフォルト = 1.0(実数) |
|
AscaleA | 角度ベースの関数に対する横軸のスケールファクター デフォルト = 1.0(実数) |
|
AscaleD | 距離ベースの関数に対する横軸のスケールファクター デフォルト = 1.0(実数) |
|
Pext | 外圧 (実数) |
|
T0 | 初期温度。 デフォルト = 295(実数) |
[K] |
Iequi | 初期熱力学的つり合いフラグ
(整数) |
|
Ittf | 排気時間シフトフラグ。注入センサーが指定されている場合のみアクティブ。
|
|
初期温度における比熱の比 5
(実数) |
||
cpai | cpa 係数(関係式cpi(T)の) (実数) |
|
cpbi | cpb 係数(関係式cpi(T)の) (実数) |
|
cpci | cpc 係数(関係式cpi(T)の) (実数) |
|
Njet | インジェクタの数 (整数) |
|
比熱の比
(実数) |
||
cpa | cpa 係数(関係式cp(T)の) (実数) |
|
cpb | cpb 係数(関係式cp(T)の) (実数) |
|
cpc | cpc 係数(関係式cp(T)の) (実数) |
|
fct_IDmas | 注入気体の質量対時間関数の識別子 (整数) |
|
Iflow | 質量対時間関数の入力タイプフラグ
(整数) |
|
Fscalemas | 質量関数に対するスケールファクター デフォルト = 1.0(実数) |
または |
fct_IDT | 注入気体の温度対時間識別子 (整数) |
|
FscaleT | 温度スケールファクター デフォルト = 1.0(実数) |
|
sens_ID | 注入を開始するセンサー識別子 (整数) |
|
Isjet | インジェクタサーフェス識別子(インジェクタごとに異なる必要があります) (整数) |
|
fct_IDvel | 注入気体速度識別子 (整数) |
|
Fscalevel | 注入気体スケールファクター デフォルト = 1.0(実数) |
|
Nvent | ベントホールの数 (整数) |
|
surf_IDv | ベントホール膜サーフェスまたは多孔サーフェス識別子 (整数) |
|
Avent | surf_IDv ≠ 0の場合: サーフェスに対するスケールファクター デフォルト = 1.0 surf_IDv = 0の場合: ベントホールのサーフェス デフォルト = 0.0(実数) |
、surf_IDV = 0の場合 |
Bvent | surf_IDv ≠ 0の場合: 衝撃を受けるサーフェスに対するスケールファクター デフォルト = 1.0 surf_IDv = 0の場合: Bvent は0にリセットされます。 デフォルト = 0.0(実数) |
、surf_IDV = 0の場合 |
Itvent | 排気の定式化 7
(整数) |
|
Tvent | 排気の開始時間 デフォルト = 0.0(実数) |
|
ベントホール膜を開口する圧力差(
= Pdef - Pext) (実数) |
||
圧力がPdefを超えてベントホール膜が開口するまでの最短時間 (実数) |
||
fct_IDV | 流出速度関数識別子 (整数) |
|
FscaleV | 右記のスケールファクター; fct_IDV デフォルト = 1.0(実数) |
|
IdtPdef |
に到達した場合の時間遅延フラグ:
|
|
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 (実数) |
|
L2 | 長さL2 (実数) |
|
L3 | 長さL3 (実数) |
|
Nb1 | 方向1における有限体積の数 デフォルト = 1(整数) |
|
Nb2 | 方向2における有限体積の数 デフォルト = 1(整数) |
|
Nb3 | 方向3における有限体積の数 デフォルト = 1(整数) |
|
grbric_ID | ユーザー定義のソリッドグループの識別子 (整数) |
|
surf_IDin | 内部サーフェス識別子 25 (整数) |
|
Iref | 参照ジオメトリ上に自動FVMメッシュを適用するためのフラグ 24
(整数) |
|
Igmerg | 全体結合定式化のフラグ 19 デフォルト = 1(整数) |
|
Cgmerg | 全体結合の係数 19 デフォルト = 0.02(実数) |
|
Cnmerg | 近傍結合の係数 19 (実数) |
|
Ptole | 有限体積識別用の許容値 デフォルト = 10-5(実数) |
|
qa | 2次体積粘性 デフォルト = 0.0(実数) |
|
qb | 線形体積粘性 デフォルト = 0.0(実数) |
|
Hmin | 浸透性を有する3角形の最小の高さ 21 (実数) |
|
Ilvout | 出力レベル: 0または1 デフォルト = 0(整数) |
|
Nlayer | V3方向に折りたたまれているエアバッグの層数の推定値(コメント23) 22 デフォルト = 10(整数) |
|
Nfacmax | 最初の自動メッシングステップにおける有限体積の影響を受けるエアバッグセグメントの推定最大数 デフォルト = 20(整数) |
|
Nppmax | 多角形の頂点の推定最大数 デフォルト = 20(整数) |
|
Ifvani | 有限体積のRadioss StarterアニメーションA000ファイルへの書き込みフラグ
(整数) |
コメント
- surf_IDex は、3節点または4節点シェル要素(ボイド要素の場合もあります)に関連付けられたセグメントを使用して定義する必要があります。
- 体積は、閉じていて法線が外向きである必要があります。
- 横軸のスケールファクターは、エアバッグ関数の横軸の単位を変換するために使用されます。たとえば:
(1) ここで、tは時間です。(2) ここで、pは圧力です。
- 初期圧力はPextに設定されます。
- = 0の場合、エアバッグに初期充填されている気体の特性は、第1インジェクターによって供給される気体の特性に設定されます。
- FVMBAG内の気体流は有限体積を使用して求められます。
これらの有限体積の一部は、エアバッグ内部に配置され、特定のパートまたは内部体積全体を満たしているソリッドのグループを通してユーザーが入力できます。ユーザー定義ソリッドによって離散化されていない内部体積部分がまだ存在する場合、残りの体積は自動メッシュ生成プロシージャによって生成されます。この機能を使用して、例えばキャニスターをモデル化できます。
有限体積は、三角端面のセット内に存在します。これらの頂点は、必ずしもエアバッグの節点と一致する必要はありません。エアバッグの外皮をモデル化する場合、4節点膜または3節点膜を使用できますが、3節点膜の使用を推奨します。 - 流出速度は次の式で与えられます。
(3) 質量の流出速度は次の式によって定義されます。(4) エネルギーの流出速度は次の式によって定義されます:(5) 通気速度は、次の式で計算されます:(6) - T > Tventの場合、または よりも長い時間、圧力がPdefを超えた場合は、ベントホール膜は収縮します。
- surf_IDv ≠ 0(surf_IDvが定義されている)場合。
(7) ここで、A歯サーフェスsurf_IDの面積、A0はサーフェスsurf_IDvの初期面積です。
- surf_IDv = 0の場合(surf_IDvは未定義)、ベントホールは無視されます。
(8) - 関数fct_IDt'とfct_IDP’を指定しなかった場合(識別子が0)、これらは1に等しいとみなされます。
- 関数fct_IDAを指定しなかった場合、fct_IDA(A/A0) = 1であるとみなされます。
- ベントホールのサーフェスは、次のように計算されます。
(9) ここで、衝撃を受けるサーフェスは:(10) 衝撃を受けないサーフェスは:(11) ここで、ベントホールsurf_IDvの各要素eに対して、nc(e)は、要素を定義している節点n(e)のうち衝撃を受ける節点の数を表しています。
- 関数fct_IDt'とfct_IDP’を指定しなかった場合(識別子が0)、これらは1に等しいとみなされます。
- 接触中に空隙率を使用するためには、関係するインターフェース(インターフェースタイプ5および7の行3)において、フラグIBAGを1に設定する必要があります。そうしなかった場合は、衝撃でインターフェースに押し込まれる節点は、AimpactedとAnon_impactedに関する前の式の被衝撃節点と見なされません。
- 自動有限体積メッシングパラメータ
- 有限体積は2つのステップで生成されます。
- 第1ステップでは、エアバッグの外皮上に排他的に配置された頂点が生成されます。これにより、外皮の変形に合わせて有限体積を更新したり、以下の手順に対応したりできます(分かりやすくするために2次元で表示):
この手順では、方向V3(指定された切断方向)と方向V1の入力が必要です。切断方向に垂直な平面内の2つ目の方向V2が計算されます。V1とV2の両方向で有限体積を配置して、切断幅を決定するためには、原点Oだけでなく、原点から正と負の両方向でカウントされた長さLiとステップ数Niも指定する必要があります。カッティング幅は次のように与えられます:Wi = 2Li / Ni
原点Oおよび長さLi(Oからの正方向と負方向の両方を含む)によって水平面(V3に対して垂直)内に描かれるボックスには、水平面に投影されたメッシュにその体積の外皮の境界ボックスが含まれている必要があります。これは、この体積全体が有限体積に分割されることを保証するために必要です。
- 第2ステップでは、有限体積の水平切断が実行されますが、堅く折りたたまれたエアバッグの場合はほとんど役に立たない可能性があります。この作業は、特に、エアバッグを広げる前に注入気体によって満たされたキャニスター内で注入を行う場合に必要です。この第2ステップでは、エアバッグ内部に配置された頂点が生成される場合があります。エアバッグの膨張に合わせて頂点を移動するために、エアバッグの外皮上に配置された2つの頂点間の垂直セグメント(方向V3と並行)にそれぞれの頂点が取り付けられます(図 6)。その参照セグメント内部の頂点の局所座標は、膨張プロセスを通して一定に保たれます。
水平切断幅は次の式で与えられます:W3 =2L3 /N3原点Oおよび長さL3(正方向と負方向の両方を含む)によってV3方向に与えられるセグメントには、V3方向のメッシュ投影に体積の外皮の境界ボックスが含まれる必要はありません。これは、2番目のステップでは既存の有限体積のみがカットされるためです。
- 第1ステップでは、エアバッグの外皮上に排他的に配置された頂点が生成されます。これにより、外皮の変形に合わせて有限体積を更新したり、以下の手順に対応したりできます(分かりやすくするために2次元で表示):
- 自動メッシングに実際に使用されるベクトルV1は、ベクトルV3に対して入力ベクトルが直交化されたあとで取得されます。
- エアバッグの膨張プロセス中に有限体積がエラーになる(体積が負の内部質量になるか、エネルギーが負になる)場合は、計算を継続するために、その近傍のいずれかと結合されます。次の2つの結合アプローチが使用されます。
- 全体結合: 有限体積は、その体積がすべての有限体積の平均体積に特定の係数をかけた値より小さくなった場合に結合されます。Igmergフラグは、平均体積として現在の平均体積を使用するか(Igmerg =1)、初期平均体積を使用するか(Igmerg =2)を指定します。平均体積から最小体積を求めるための係数は、Cgmergです。
- 近傍結合: 有限体積は、その体積が近傍の平均体積に特定の係数をかけた値より小さくなった場合に結合されます。平均体積から最小体積を求めるための係数は、Cnmergです。
- CgmergとCnmerg の両方が0でない場合は、両方の結合アプローチが同時に使用されることを意味します。強い衝撃を与える場合は、qa = 1.1およびqb = 0.05を設定することを推奨します。
- 2層の布が物理的に接触している場合は、有限体積間に流れが発生しないはずです。インターフェースギャップが原因で数値的にはこのようになりません。Hminは、浸透性を有する三角面の最小の高さを表しています。この値よりも小さい三角面は、不浸透性を有します。その値は、エアバッグのセルフ接触インターフェースのギャップに近づける必要があります。
- Nlayer、Nfacmax、Nppmaxは、有限体積作成プロセスを支援するメモリパラメータです。これらの値を変更しても、計算が停止することはありません。値を増やすと、自動メッシングのためのメモリ量が増え、計算時間が短縮されます。
- 有限体積作成プロセスでは、最初に、平面多角形が作成され、その後で、閉じた多面体に組み立てられ、三角端面に分解されます。 Nppmax は、これらの多角形に含まれる頂点の最大数です。
- Iref = 1の設定はモニター体積タイプFVMBAG または FVMBAG1に対して/REFSTAに基づいた参照ジオメトリにのみ働きます(/XREFに基づいた参照ジオメトリはまだサポートされていません)。
- 内部サーフェスまたはバッフルをモニター体積内の気体流れの障害物として考慮に用いるには、有限体積法のみ可能です。内部サーフェスは、モニター体積の外部および内部サーフェス節点と一致する節点があるTETRA4(HEXAおよびPENTAの場合もあります)のように、モニター体積がソリッド要素で満たされている場合のみ、FVMで考慮されます。(これらのソリッドはgrbrick_IDで定義する必要があります)。材料モデルがLAW19の内部サーフェス繊維にのみ、空隙率を、0(空隙率なし)から1(空隙率完全)(ベント)の範囲で適用できます。気体流の方向が内部サーフェスの法線方向と反対の場合は、内部サーフェス上でインジェクターサーフェスも定義できます。
- ベントホールのサーフェスの要素(surf_IDv)がインジェクター(surf_IDinj)に属している場合、この要素はベントホールから無視されます。t = 0の時点で計算される定数cの補正係数f がベントホールのサーフェス全体に適用されます
(12)