JP5691905B2 - インフレータ性能試験装置、方法及びプログラム - Google Patents
インフレータ性能試験装置、方法及びプログラム Download PDFInfo
- Publication number
- JP5691905B2 JP5691905B2 JP2011155166A JP2011155166A JP5691905B2 JP 5691905 B2 JP5691905 B2 JP 5691905B2 JP 2011155166 A JP2011155166 A JP 2011155166A JP 2011155166 A JP2011155166 A JP 2011155166A JP 5691905 B2 JP5691905 B2 JP 5691905B2
- Authority
- JP
- Japan
- Prior art keywords
- inflator
- tank
- internal pressure
- flow rate
- injection
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Air Bags (AREA)
Description
本発明は、インフレータの性能を試験するインフレータ性能試験装置、方法及びプログラムに関する。
代表的な移動体としての車両にはエアバッグ装置が具備されるようになってきている。このエアバッグ装置はインフレータを具備して、インフレータによりガスを発生させてエアバッグを展開させている。このようなエアバッグ装置の性能評価や部品開発のためには、エアバッグを展開させるインフレータの特性を把握することが重要である。ところが、エアバッグの展開時間は数十ミリ秒という短時間のため、インフレータの特性を計測することは困難であった。特に、インフレータから噴射されるガスの流量及び温度の変化等のインフレータ特性は、直接計測することが困難であった。
特許文献１には、タンク試験を実施してタンクの圧力時間特性を求める技術が記載されている。このタンク試験ではタンク内に固定したインフレータに着火してからのタンク内の圧力を逐次計測することで、タンク内の圧力上昇変化を時系列的に計測している。
ところで、インフレータ特性は、インフレータからのガスの噴射流量及び噴射温度の変動を把握することが重要である。しかし、タンク試験によるタンク内圧上昇変化の時系列的な計測結果のみでは、これらを特定することが困難である。例えば、このタンク内圧上昇変化の時系列的な計測結果から噴射流量及び噴射温度のインフレータ特性を得ようとすると、噴射流量及び噴射温度の何れか一方を仮定して他方を同定しなければならない。このため、計算負荷が増大すると共に、高精度の計算結果を得ることができない。
近年、部品開発の効率向上のために、コンピュータによる模擬的な試験が実施されている。例えば、部品の各種試験では、ＣＡＥ（Computer Aided Engineering）解析による性能評価が行われている。このＣＡＥ解析では、正確かつ高速に解析結果を得るためのに、エアバッグ装置に具備される実際のインフレータを高い再現性で模擬できるインフレータ特性を得ることが望まれている。
本発明は上記事実を考慮し、タンク内圧の時間特性からインフレータの噴射流量及び噴射温度のインフレータ特性を同定するインフレータ性能試験装置、方法及びプログラムを得ることを課題とする。
請求項１に記載の発明のインフレータ性能導出装置は、比熱が異なる気体が充填された複数のタンクの各々に同一のインフレータを取り付けて、前記インフレータから気体を噴射させたときの前記複数のタンクの各々の内圧の時間特性を測定した測定結果を取得する取得手段と、前記取得手段により取得した各タンクの内圧の時間特性に基づいて、タンクの内圧とインフレータの噴射流量とタンク内の温度との関係から各タンクの内部エネルギの変化量が一致するようにインフレータの噴射流量を同定すると共に、タンクの内部エネルギの時系列の収支が前記タンクの内圧の時間特性に対応するようにインフレータの噴射温度の時間特性を同定することにより、インフレータ性能としてインフレータの噴射流量及び噴射温度の時間特性を算出する算出手段と、を備える。
本発明のインフレータ性能導出装置では、取得手段によって、比熱が異なる気体が充填された複数のタンクの各々に同一のインフレータを取り付けて、前記インフレータから気体を噴射させたときの前記複数のタンクの各々の内圧の時間特性を測定した測定結果を取得する。算出手段は、取得手段により取得した各タンクの内圧の時間特性に基づいて、インフレータ性能としてインフレータの噴射流量及び噴射温度の時間特性を算出する。具体的には、算出手段は、タンクの内圧とインフレータの噴射流量とタンク内の温度との関係から各タンクの内部エネルギの変化量が一致するようにインフレータの噴射流量を同定する。また、算出手段は、タンクの内部エネルギの時系列の収支が前記タンクの内圧の時間特性に対応するようにインフレータの噴射温度の時間特性を同定する。
従って、本発明によれば、異なる比熱の気体に応じて、相違する各タンクの内圧の時間特性を用いて、各タンクの内部エネルギの変化量が一致するようにタンクの内圧とインフレータの噴射流量とタンク内の温度との関係からインフレータの噴射流量を同定できる。また、タンクの内部エネルギの時系列の収支がタンクの内圧の時間特性に対応するようにインフレータの噴射温度の時間特性を同定できる。このように、タンクの内圧の時間特性からインフレータの噴射流量及び噴射温度の時間特性を得ることができるので、インフレータ特性を得るための演算負荷を軽減できる。すなわち、異なる比熱の気体によるタンクの内圧の時間特性を用いて、インフレータの噴射流量及び噴射温度の時間特性を同定するため、短時間で処理を終了することができる。
請求項２に記載の発明は、請求項１に記載のインフレータ性能導出装置において、前記算出手段は、数１のエネルギ収支式および数２の状態方程式の支配方程式を用いて、インフレータの噴射流量及び噴射温度の時間特性を同定することを特徴とする。
ｊはタンクを示す識別子である。
Ｈinf （ｔ）は時刻０から時刻ｔまでにタンクに流入する気体のエンタルピを示す。
Ｕj （ｔ）−Ｕj （０）は時刻０から時刻ｔまでにタンク内の気体の内部エネルギの変化量を示す。
Ｐj （ｔ）はタンク内圧を示す。
Ｖj はタンク容量を示す。
ｎj はタンク内の初期気体の気体量を示す。
ｎinf （ｔ）はインフレータの気体の噴出量を示す。
Ｒは気体定数を示す。
Ｔj （ｔ）はタンク内の気体温度を示す。
Ｈinf （ｔ）は時刻０から時刻ｔまでにタンクに流入する気体のエンタルピを示す。
Ｕj （ｔ）−Ｕj （０）は時刻０から時刻ｔまでにタンク内の気体の内部エネルギの変化量を示す。
Ｐj （ｔ）はタンク内圧を示す。
Ｖj はタンク容量を示す。
ｎj はタンク内の初期気体の気体量を示す。
ｎinf （ｔ）はインフレータの気体の噴出量を示す。
Ｒは気体定数を示す。
Ｔj （ｔ）はタンク内の気体温度を示す。
算出手段は、エネルギ収支式および状態方程式の支配方程式を用いてインフレータの噴射流量及び噴射温度の時間特性を求めるので、同定するインフレータの噴射流量及び噴射温度の時間特性は実際のインフレータの性能が反映されたものになり、精度が向上する。
請求項３に記載の発明は、請求項１又は請求項２に記載のインフレータ性能導出装置において、前記算出手段は、比熱を多項式で近似することを特徴とする。算出手段では、比熱を多項式で近似することによって、算出負荷を軽減することができる。
請求項４に記載の発明のインフレータ性能導出方法は、比熱が異なる気体が充填された複数のタンクの各々に同一のインフレータを取り付けて、前記インフレータから気体を噴射させたときの前記複数のタンクの各々の内圧の時間特性を測定した測定結果を取得する取得工程と、前記取得工程により取得した各タンクの内圧の時間特性に基づいて、タンクの内圧とインフレータの噴射流量とタンク内の温度との関係から各タンクの内部エネルギの変化量が一致するようにインフレータの噴射流量を同定すると共に、タンクの内部エネルギの時系列の収支が前記タンクの内圧の時間特性に対応するようにインフレータの噴射温度の時間特性を同定することにより、インフレータ性能としてインフレータの噴射流量及び噴射温度の時間特性を算出する算出工程と、を有する。
本発明の方法によれば、タンクの内圧の時間特性からインフレータの噴射流量及び噴射温度の時間特性を得ることができるので、インフレータ特性を得るための演算負荷を軽減できる。
なお、請求項５に記載の発明のように、コンピュータを、請求項１〜３の何れか１項に記載のインフレータ性能導出装置の各手段として機能させるためのプログラムとすることを特徴とする。このプログラムは、記憶媒体に記憶して流通可能なようにしてもよい。
詳細には、この発明のプログラムは、比熱が異なる気体が充填された複数のタンクの各々に同一のインフレータを取り付けて、前記インフレータから気体を噴射させたときの前記複数のタンクの各々の内圧の時間特性を測定した測定結果を取得する取得ステップと、前記取得手段により取得した各タンクの内圧の時間特性に基づいて、タンクの内圧とインフレータの噴射流量とタンク内の温度との関係から各タンクの内部エネルギの変化量が一致するようにインフレータの噴射流量を同定すると共に、タンクの内部エネルギの時系列の収支が前記タンクの内圧の時間特性に対応するようにインフレータの噴射温度の時間特性を同定することにより、インフレータ性能としてインフレータの噴射流量及び噴射温度の時間特性を算出する算出ステップと、の各ステップをコンピュータに実行させることでインフレータ性能を導出させる。
本発明のプログラムによれば、コンピュータにおいて短時間でインフレータ性能の同定を終了させ、演算負荷を軽減させることができる。
本発明によれば、タンクの内圧の時間特性からインフレータの噴射流量及び噴射温度の時間特性を得ることができるので、インフレータ特性を得るための演算負荷を軽減できる。
以下、図面を参照して本発明の実施の形態の一例を詳細に説明する。本実施形態は、エンジン等の内燃機関の動力で走行する車両や、エンジンとモータを備えたハイブリッド車両等に搭載されるエアバッグ装置に具備されるインフレータの性能を同定する場合に、本発明を適用したものである。
＜インフレータ性能導出装置＞
図２に示すように、本実施の形態のインフレータ性能導出装置１０は、ＣＰＵ１４、ＲＡＭ１６，ＲＯＭ１８を含むコンピュータ本体１２を備えている。コンピュータ本体１２は、入出力インタフェース２０を含んでおり、入出力インタフェース２０は、ＣＰＵ１４、ＲＡＭ１６、ＲＯＭ１８とデータやコマンドを授受可能にバス２２を介して接続されている。
図２に示すように、本実施の形態のインフレータ性能導出装置１０は、ＣＰＵ１４、ＲＡＭ１６，ＲＯＭ１８を含むコンピュータ本体１２を備えている。コンピュータ本体１２は、入出力インタフェース２０を含んでおり、入出力インタフェース２０は、ＣＰＵ１４、ＲＡＭ１６、ＲＯＭ１８とデータやコマンドを授受可能にバス２２を介して接続されている。
入出力インタフェース２０には、入出力装置２６として、ポインティングデバイスであるマウスやキーボードが接続されており、また、データやコマンド、その他の各種画像を表示するためのディスプレイ２４が接続されている。また、入出力インタフェース２０には、同定部３０，検証部３２、及び記録装置３４が接続されている。
同定部３０は、後述するインフレータ特性同定処理を担当する処理部であり、インフレータ特性としてインフレータの噴射流量及び噴射温度の時間特性を算出することにより同定する処理部である。
また、同定部３０は、上述のハードウェア資源と後述するソフトウェア資源との複合的な動作により機能する装置である。すなわち、同定部３０は詳細を後述する処理のフローチャートを実行させるためのプログラムを記憶したメモリを具備している。従って、ハードウェア資源とソフトウェア資源との構成及び複合的な処理によって、同定部３０を含むインフレータ性能導出装置１０は、インフレータ特性同定装置として機能することになる。
検証部３２は、ＣＡＥ解析処理を担当する処理部であり、対象部品であるインフレータについて性能解析を含むＣＡＥ解析する処理部である。このＣＡＥ解析部３２は、上述のハードウェア資源と後述するソフトウェア資源との複合的な動作により機能する装置である。すなわち、ＣＡＥ解析部３２は処理実行のためのプログラムを記憶したメモリを具備している。従って、ハードウェア資源とソフトウェア資源との構成及び複合的な動作によって、ＣＡＥ解析部３２を含むインフレータ性能導出装置１０は、解析装置として機能することになる。なお、検証部３２は、後述するように同定されたインフレータ特性を検証するための装置として備えたものであり、本実施形態のインフレータ性能導出装置１０において必須の構成とするものではなく、省略してもよい。検証部３２を省略する場合には、検証部３２として代替する他の装置をインフレータ性能導出装置１０に接続して検証可能に構成することができる。
記録装置３４には、タンク内圧の時間特性（内圧カーブ）の実測値が格納されている。この実測値は、例えばタンクテストによる試験等で計測されるものである（後述）。また、記録装置３４には、インフレータ性能を導出するために用いられる数式が格納されている（後述）。なお、記録装置３４は、解析に用いられる様々な部品について解析時の数値計算用の解析モデルが格納される。これら解析モデルは、様々な部品について、その節点や辺および各種付帯条件などから構成されるモデルの数値計算データが格納される。
（タンクテスト）
次に、タンク内圧の時間特性（内圧カーブ）を得るためのタンクテストの一例を説明する。本実施形態では、２個のタンクを用いてタンク内圧の時間特性（内圧カーブ）を得る場合を説明する。
次に、タンク内圧の時間特性（内圧カーブ）を得るためのタンクテストの一例を説明する。本実施形態では、２個のタンクを用いてタンク内圧の時間特性（内圧カーブ）を得る場合を説明する。
図３に示すように、まず、内容積が既知の複数タンク４４，４８を用意して、各々に性能が同一とされるインフレータ４０，４１を固定する。また、この各タンク内には圧力センサ４２，４６が設置される。各タンク４４，４８は、比熱が相違する気体４３，４７（図３では、初期ガスＡと初期ガスＢ）を充填する。本実施形態では、一例として、窒素ガスとアルゴンガスを比熱が相違する気体として用いている（図４参照）。各々のタンクを密封後に、図示しない着火装置により、インフレータ４０，４１を作動させる。この図示しない着火装置によるインフレータの作動開始からのタンク内圧変動を圧力センサにより検出する。すなわち、インフレータの作動開始から連続時間的または一定時間おきにタンク内圧の計測値を得る。従って、タンク４４，４８毎に、タンク内圧と時間のデータセットが計測値として記録装置３４に格納される。
図５は、タンク内圧と時間のデータセットを視覚化した、比熱が相違する気体４３，４７（初期ガスＡ：窒素ガス、初期ガスＢ：アルゴンガス）を充填したタンク４４，４８の各々のタンク内圧時間特性５０，５２を示した一例である。比熱が相違する気体によるタンクでは、インフレータによるガス噴射時のタンク内圧の変動量に相違が発生する。
なお、入出力インタフェース２０には、記録媒体としての磁気円盤が挿抜可能な記録装置や記録領域を備えた電子回路を具備する補助記憶装置が接続されるようにしてもよい。なお、後述する処理ルーチン等のプログラムは、記録装置や補助記憶装置を用いて磁気円盤や記録領域に対して読み書き可能である。従って、後述する処理ルーチン等のプログラムは、ＲＯＭ１８に記憶することなく、予め磁気円盤や記憶領域に記録しておき、記録装置や補助記憶装置を介して、予め記録された処理プログラムを実行してもよい。また、コンピュータ本体１２には、ハードディスク装置等の大容量記憶装置（図示省略）を接続し、プログラムを大容量記憶装置（図示省略）へ格納（インストール）して実行するようにしてもよい。また、記録媒体としては、ＣＤ−ＲＯＭ，ＭＤ，ＭＯ，ＤＶＤ等のディスク、ＤＡＴ等の磁気テープ、そして記録領域としてはフラッシュメモリ（flash memory）やＳＳＤ(Solid State Drive)などの補助記憶装置があり、これらを用いるときには、ＣＤ−ＲＯＭ装置、ＭＤ装置、ＭＯ装置、ＤＶＤ装置、ＤＡＴ装置、補助記憶装置等を用いればよい。
＜インフレータ特性同定＞
ここで、インフレータ特性同定の概要を説明する。本発明者は、計測により得られるデータがタンク内圧の時間特性（内圧カーブ）のみであるとき、インフレータ特性としてインフレータの噴射温度および噴射流量を、同定することができない、という点から、相違する内圧カーブを得ることができれば、双方について同定可能であるという知見に至った。すなわち、１つの内圧カーブのみのときは、例えば、インフレータによる気体の噴射温度が一定値という仮定をした上で同定を進める必要があり、十分なインフレータ特性を同定することができなかった。そこで、相違する内圧カーブを得ることができれば、一定値による仮定をすることなく、インフレータの噴射温度および噴射流量を、同定できる、という知見に至った。
ここで、インフレータ特性同定の概要を説明する。本発明者は、計測により得られるデータがタンク内圧の時間特性（内圧カーブ）のみであるとき、インフレータ特性としてインフレータの噴射温度および噴射流量を、同定することができない、という点から、相違する内圧カーブを得ることができれば、双方について同定可能であるという知見に至った。すなわち、１つの内圧カーブのみのときは、例えば、インフレータによる気体の噴射温度が一定値という仮定をした上で同定を進める必要があり、十分なインフレータ特性を同定することができなかった。そこで、相違する内圧カーブを得ることができれば、一定値による仮定をすることなく、インフレータの噴射温度および噴射流量を、同定できる、という知見に至った。
まず、比熱が相違する気体をタンク内に充填させた２パターンのタンクテストを行うことで、２パターンの内圧カーブを得ることができる。各々のタンクでは、エネルギ収支についての現象をエネルギ収支式で表現できる。また、各々のタンクでは、気体の状態についての現象を状態方程式で表現できる。従って、これらのエネルギ収支式および状態方程式を支配方程式とすることで、各々のタンクにおける現象を数式で記述でき、各々のタンクに同一性能のインフレータを用いたことを共通点とすることで、そのインフレータ特性を同定することができるという知見に至った。
インフレータ特性同定は、次の４つのプロセスからなる。第１プロセスではタンクテストを実施し、第２プロセスではタンクテスト毎のタンク内圧の時間特性（内圧カーブ）を得て、第３プロセスでは支配方程式を導出し、第４プロセスでは支配方程式からインフレータの気体の噴射流量及び噴射温度を同定する。
まず、第１プロセスでは、タンクテストを実施する。すなわち、比熱が相違する気体が充填された各タンクでインフレータを作動させ、インフレータの作動開始からの時系列的なタンク内圧の計測値を得る（図３参照）。本実施形態では、比熱が相違する気体としての窒素を初期ガスＡとして行うタンクテストを、第１タンクテストとし、アルゴンを初期ガスＢとして行うタンクテストを、第２タンクテストとして説明する。
第２プロセスでは、第１及び第２タンクテストで得られたタンク内圧の計測値から、タンクテスト毎のタンク内圧の時間特性（内圧カーブ）を作成する。この内圧カーブは、ことなる特性となる２種類の内圧カーブになる（図５参照）。
第３プロセスでは、支配方程式を導出する。すなわち、タンクテストの現象を再現記述する数式を定める。
まず、第１タンクテストのエネルギ収支式は、次の（１）式で表すことができ、状態方程式は、（２）式で表すことができる。
但し、
Ｈinf （ｔ）は時刻０から時刻ｔまでにタンクに流入する気体のエンタルピを示す。
Ｕ1 （ｔ）−Ｕ1 （０）は時刻０から時刻ｔまでにタンク４４内の気体の内部エネルギの変化量を示す。
但し、
Ｐ1 （ｔ）はタンク４４の内圧を示す。
Ｖ1 はタンク４４の容量を示す。
ｎ1 はタンク４４内の初期気体（窒素）の気体量を示す。
ｎinf （ｔ）はインフレータ４０の気体の噴出量を示す。
Ｒは気体定数を示す。
Ｔ1 （ｔ）はタンク４４内の気体温度を示す。
次に、第２タンクテストのエネルギ収支式は、次の（３）式で表すことができ、状態方程式は、（４）式で表すことができる。
但し、
Ｈinf （ｔ）は時刻０から時刻ｔまでにタンクに流入する気体のエンタルピを示す。
Ｕ2 （ｔ）−Ｕ2 （０）は時刻０から時刻ｔまでにタンク４８内の気体の内部エネルギの変化量を示す。
但し、
Ｐ2 （ｔ）はタンク４８の内圧を示す。
Ｖ2 はタンク４８の容量を示す（タンク４４と同一）。
ｎ2 はタンク４８内の初期気体（窒素）の気体量を示す（タンク４４と同一）。
ｎinf （ｔ）はインフレータ４１の気体の噴出量を示す（インフレータ４０と同一）。
Ｒは気体定数を示す。
Ｔ2 （ｔ）はタンク４８内の気体温度を示す。
ここで、気体の比熱については、次式で表すことができる。まず、定圧比熱Ｃp、定積比熱Ｃvの各々について、第１タンクテストの初期ガスＡである窒素については、（５）式で表わすことができる。
但し、ａ１，ｂ１，ｃ１，は定数である。
但し、ａ１，ｂ１，ｃ１，は定数である。
また、第２タンクテストの初期ガスＢであるアルゴンについては、（６）式で表すことができる。
但し、ａ２，ｂ２，ｃ２，は定数である。
但し、ａ２，ｂ２，ｃ２，は定数である。
また、インフレータ噴出ガスについては、（７）式で表すことができる。
但し、Ａ，Ｂ，Ｃ，は定数である。
但し、Ａ，Ｂ，Ｃ，は定数である。
上記の（１）式は、（５）式〜（７）式を用いて、書き換えることができる。まず、（１）式の左項は、次の（８）式に書き換えることができる。
また、（１）式の右項は、次の（９）式に書き換えることができる。
また、（３）式の右項は、次の（１０）式に書き換えることができる。
第４プロセスでは、インフレータの気体の噴射流量及び噴射温度を同定する。すなわち、実測値である第２プロセスで得られたタンクテスト毎のタンク内圧の時間特性（内圧カーブ）と、第３プロセスで得られたタンクテストの現象を記述した数式を用いて、インフレータの気体の噴射流量及び噴射温度を同定する。この第４プロセスは、上記の支配方程式を連立し、その結果式から一定の解法によりインフレータの気体の噴射流量及び噴射温度を同定する。第４プロセスは、次の３段階の処理からなる。第１段階処理はインフレータの噴出量を求める処理であり、第２段階処理は流量に変換する処理であり、第３段階処理は噴射温度を求める処理である。
まず、第４プロセスの第１段階処理は、インフレータの気体の噴出量ｎinf （ｔ）を求める処理である。ここでの処理は、「内部エネルギの変化量は同一である」という事項から導出することができる。この事項は、次の（１１）式で示すことができる。この（１１）式は、上記の支配方程式を整理することにより、（１２）式で表すことができる。
上記（１２）式の両辺に次の（１３）式を乗算し、その結果の６次方程式を一定の解法で解くことにより、インフレータの気体の噴出量ｎinf （ｔ）を同定する。ここでは、数値解析的に解きやすくするために、乗算している。この６次方程式についての一定の解法には、二分法やニュートン法などがある。
次に、第４プロセスの第２段階処理は、インフレータの気体の流量を求める処理である。この処理は、構造解析ソフトウェア（例えば、ＬＳ−ＤＹＮＡ）を利用した解析をするとき、時刻毎の流量を入力する必要があるため、気体の噴出量から流量に変換する必要があるため、実行される処理である。このことは、上記で同定したインフレータの気体の噴出量ｎinf （ｔ）を流量に変換する処理に相当する。その変換処理の結果式は、次の（１４）式で示すことができる。この（１４）式は、上記の支配方程式を整理することで表すことができる。なお、（１４）式において、Ｍは質量を示している。
上記（１４）式の中心差分を求めることにより、インフレータの気体の流量を同定する。
次に、第４プロセスの第３段階処理は、インフレータの気体の噴出温度Ｔinf （ｔ）を求める処理である。ここでの処理は、「エネルギ収支式を満たす」という事項から導出することができる。この事項は、上記（１）式、すなわち（８）式と（９）式が成立（または（２）式について適用）することを意味する。このことについて、（８）式を左辺式とし、（９）式を右辺式として次の（１５）式に示す。
この（１５）式の右辺項について、（２）式から求まる次の（１６）式を代入する。これにより、（１５）式の右辺項は既知のデータにより同定できる。
次に、（１５）式の左辺項について、次の（１７）式で示すときを、Ｑｉとし、台形積分を用いて整理すると、（１８）式および（１９）式で表すことができる。
上記（１９）式は、漸化式形式の３次方程式である。この（１９）式では、Ｔｉｎｆ（ｉ−１）及びその２次と３次の変数は既知変数であるが、それ以外のＴｉｎｆの変数は未知である。そこで、漸化式形式の（１９）式について、ｔ＝０から順次一定の解法で解くことにより、インフレータの気体の噴出温度Ｔinf （ｔ）を同定する。この一定の解法には、解析的解法、二分法やニュートン法などの数値解法がある。
以上のプロセスにて、比熱が相違する気体による内圧データを用いて、インフレータ特性を同定できる。
＜インフレータ特性同定処理＞
図１には本実施の形態のインフレータ性能導出装置１０において実行されるインフレータ特性同定処理の流れを示した。なお、本実施形態では、インフレータ性能としてインフレータ特性のうちのインフレータの気体（ガス）の噴射流量及びインフレータの気体（ガス）の噴射温度を同定する処理例を説明する。
図１には本実施の形態のインフレータ性能導出装置１０において実行されるインフレータ特性同定処理の流れを示した。なお、本実施形態では、インフレータ性能としてインフレータ特性のうちのインフレータの気体（ガス）の噴射流量及びインフレータの気体（ガス）の噴射温度を同定する処理例を説明する。
まず、ステップ１００はインフレータ特性同定処理の前処理であり、予めタンクテストが実施される。すなわち、気体４３（窒素：初期ガスＡ），が充填されたタンク４４でインフレータ４０を作動させ、インフレータの作動開始からの時系列的なタンク内圧の計測値を第１タンクテストのデータとして得る。同様に、気体４３とは比熱が相違する気体４７（アルゴン：初期ガスＢ）が充填されたタンク４８でインフレータ４１を作動させ、インフレータの作動開始からの時系列的なタンク内圧の計測値を第２タンクテストのデータとして得る。これらの第１タンクテストのデータ及び第２タンクテストのデータを記録装置３４に格納しておく。
なお、ステップ１００では、予め実施したタンクテストの結果を記録装置３４に記憶させた場合を説明したが、タンクテストの装置をインフレータ性能導出装置１０に接続してタンクテストのデータを得るようにしてもよい。
次のステップ１０２では、記録装置３４から予め実施したタンクテストの結果を取得し、次のステップ１０４において内圧カーブを算出する。すなわち、タンクテストで得られたタンク内圧の計測値から、タンクテスト毎のタンク内圧の時間特性（内圧カーブ）を算出する（図５）。
以上のようにして、比熱が相違する気体によるタンク内圧の時間特性を求めた後には、ステップ１０６へ進み、支配方程式を導出する。すなわち、タンクテストの現象を再現記述する数式を定める。ここでの処理は、記録装置３４に記憶された演算式としての数式を読み取ることにより、タンクテストの現象を再現記述する数式を定める。具体的には、第１タンクテストのエネルギ収支式として（１）式を定め（ステップ１０８）、状態方程式として（２）式を定める（ステップ１１０）。同様に、第２タンクテストのエネルギ収支式として（３）式を定め（ステップ１１２）、状態方程式として（４）式を定める（ステップ１１４）。
このステップ１０６では、上述したように、インフレータの気体の噴出量ｎinf （ｔ）を同定するための（１２）式を定める。また、インフレータの気体の流量を同定するための（１４）式を定める。さらに、インフレータの気体の噴出温度Ｔinf （ｔ）を同定するための（１９）式を定める。
次に、ステップ１１６において、インフレータの気体の噴射流量及び噴射温度を同定する。すなわち、実測値であるステップ１０４で得られたタンクテスト毎のタンク内圧の時間特性（内圧カーブ）と、ステップ１０６で定めたタンクテストの現象を記述した数式を用いて、インフレータの気体の噴射流量及び噴射温度を同定する。
詳細には、まず、ステップ１１８において、（１２）式の６次方程式から数値解法によって、インフレータの気体の噴出量ｎinf （ｔ）を同定する。次に、ステップ１２０において、（１４）式の中心差分を求めることにより、インフレータの気体の流量を同定する。次に、（１９）式による漸化式形式の３次方程式から、解析的手法または数値解法によって、インフレータの気体の噴出温度Ｔinf （ｔ）を同定する。
以上説明したように、比熱が相違する気体を充填したタンクによるインフレータ噴射のタンクテスト結果を用いて、インフレータの気体の流量およびインフレータの気体の噴出温度を同定することができる。これにより、インフレータ特性を把握することができ、性能解析等に容易に用いることができる。また、タンクの内圧の時間特性からインフレータの噴射流量及び噴射温度の時間特性を得ることができるので、インフレータ特性を得るための演算負荷を軽減できる。
＜インフレータ特性の検証＞
次に本発明者は、本実施形態のインフレータ性能導出装置１０により同定したインフレータの気体の噴射流量およびインフレータの気体の噴出温度によるインフレータ特性について検証を行った。ここでは、構造解析ソフトウェア（ＬＳ−ＤＹＮＡ）を利用した数値実験による検証を採用した。まず、検証用として仮想のインフレータ特性を作成した。図６には、インフレータ特性を示し、（Ａ）は噴射流量、（Ｂ）は噴射温度を示す。
次に本発明者は、本実施形態のインフレータ性能導出装置１０により同定したインフレータの気体の噴射流量およびインフレータの気体の噴出温度によるインフレータ特性について検証を行った。ここでは、構造解析ソフトウェア（ＬＳ−ＤＹＮＡ）を利用した数値実験による検証を採用した。まず、検証用として仮想のインフレータ特性を作成した。図６には、インフレータ特性を示し、（Ａ）は噴射流量、（Ｂ）は噴射温度を示す。
なお、初期ガスは、窒素、アルゴンを採用し、それ以外の条件は一致させてタンクテストを行った。
図７には、上記仮想インフレータの特性およびタンクテスト条件で、構造解析ソフトウェア（ＬＳ−ＤＹＮＡ）のＣＶ法（コントロールボリューム法）による数値実験の結果として、タンクテスト毎のタンク内圧の時間特性（内圧カーブ）を示した。
図８には、図７のタンク内圧の時間特性（内圧カーブ）をもとにして、本実施形態のインフレータ性能導出装置１０において実行されるインフレータ特性同定処理によって、同定したインフレータの気体の噴射流量（図８（Ａ））およびインフレータの気体の噴出温度（図８（Ｂ））によるインフレータ特性を示した。図８から理解できるように、検証用に設定したインフレータ特性とほぼ一致するインフレータ特性が得られることを確認できる。
図９には、逆検証として、本実施形態のインフレータ性能導出装置１０のインフレータ特性同定処理により同定したインフレータ特性を用いて、構造解析ソフトウェア（ＬＳ−ＤＹＮＡ）のパーティクル法による数値実験の結果として、タンクテスト解析結果を示した。図から理解されるように、ＣＶ法による数値実験結果と、本実施形態で同定したインフレータ特性による結果とが一致している。
以上説明した検証結果からも理解できるように、本実施形態のインフレータ性能導出装置１０によりインフレータの気体の噴射流量およびインフレータの気体の噴出温度によるインフレータ特性を同定できる。
１０…インフレータ性能導出装置
３０…同定部
３４…記録装置
４０…インフレータ
４２…圧力センサ
４３…気体
４４…タンク
４６…圧力センサ
４７…気体
４８…タンク
５０…タンク内圧時間特性
５２…タンク内圧時間特性
３０…同定部
３４…記録装置
４０…インフレータ
４２…圧力センサ
４３…気体
４４…タンク
４６…圧力センサ
４７…気体
４８…タンク
５０…タンク内圧時間特性
５２…タンク内圧時間特性
Claims (5)
- 比熱が異なる気体が充填された複数のタンクの各々に同一のインフレータを取り付けて、前記インフレータから気体を噴射させたときの前記複数のタンクの各々の内圧の時間特性を測定した測定結果を取得する取得手段と、
前記取得手段により取得した各タンクの内圧の時間特性に基づいて、タンクの内圧とインフレータの噴射流量とタンク内の温度との関係から各タンクの内部エネルギの変化量が一致するようにインフレータの噴射流量を同定すると共に、タンクの内部エネルギの時系列の収支が前記タンクの内圧の時間特性に対応するようにインフレータの噴射温度の時間特性を同定することにより、インフレータ性能としてインフレータの噴射流量及び噴射温度の時間特性を算出する算出手段と、
を備えたインフレータ性能導出装置。 - 前記算出手段は、数１のエネルギ収支式および数２の状態方程式の支配方程式を用いて、インフレータの噴射流量及び噴射温度の時間特性を同定することを特徴とする請求項１に記載のインフレータ性能導出装置。
- 前記算出手段は、比熱を多項式で近似することを特徴とする請求項１または請求項２に記載のインフレータ性能導出装置。
- 比熱が異なる気体が充填された複数のタンクの各々に同一のインフレータを取り付けて、前記インフレータから気体を噴射させたときの前記複数のタンクの各々の内圧の時間特性を測定した測定結果を取得する取得工程と、
前記取得工程により取得した各タンクの内圧の時間特性に基づいて、タンクの内圧とインフレータの噴射流量とタンク内の温度との関係から各タンクの内部エネルギの変化量が一致するようにインフレータの噴射流量を同定すると共に、タンクの内部エネルギの時系列の収支が前記タンクの内圧の時間特性に対応するようにインフレータの噴射温度の時間特性を同定することにより、インフレータ性能としてインフレータの噴射流量及び噴射温度の時間特性を算出する算出工程と、
を有するインフレータ性能導出方法。 - コンピュータを、請求項１〜３の何れか１項に記載のインフレータ性能導出装置の各手段として機能させるための車両用空調制御プログラム。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2011155166A JP5691905B2 (ja) | 2011-07-13 | 2011-07-13 | インフレータ性能試験装置、方法及びプログラム |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2011155166A JP5691905B2 (ja) | 2011-07-13 | 2011-07-13 | インフレータ性能試験装置、方法及びプログラム |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2013019850A JP2013019850A (ja) | 2013-01-31 |
| JP5691905B2 true JP5691905B2 (ja) | 2015-04-01 |
Family
ID=47691418
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2011155166A Expired - Fee Related JP5691905B2 (ja) | 2011-07-13 | 2011-07-13 | インフレータ性能試験装置、方法及びプログラム |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP5691905B2 (ja) |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2264734A1 (en) * | 1996-08-30 | 1998-03-05 | Talley Defense Systems, Inc. | Gas generating compositions |
| JP2960388B2 (ja) * | 1997-03-24 | 1999-10-06 | ダイセル化学工業株式会社 | エアバッグ用ガス発生剤組成物成型体 |
-
2011
- 2011-07-13 JP JP2011155166A patent/JP5691905B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JP2013019850A (ja) | 2013-01-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP3575892B1 (en) | Model parameter value estimation device and estimation method, program, recording medium with program recorded thereto, and model parameter value estimation system | |
| Qureshi et al. | Effect of simulated combustor temperature nonuniformity on hp vane and end wall heat transfer: An experimental and computational investigation | |
| US9483605B2 (en) | Probabilistic high cycle fatigue (HCF) design optimization process | |
| US7505949B2 (en) | Process model error correction method and system | |
| Serrano et al. | Analysis and methodology to characterize heat transfer phenomena in automotive turbochargers | |
| CN105209984A (zh) | 用于确定技术系统输出值的模型的方法 | |
| Sinha et al. | Vibratory parameters of blades from coordinate measurement machine data | |
| JP2004324598A (ja) | エンジン部品の破壊予測方法、破壊予測システム及びその制御プログラム | |
| JP4186691B2 (ja) | エンジン部品の破壊予測システム及びその制御プログラム | |
| JP6181589B2 (ja) | 計測値解析装置及び計測値解析方法 | |
| JPWO2017154214A1 (ja) | Ｗｉｅｂｅ関数パラメータ同定装置、方法、プログラム、内燃機関状態検出装置、及び車載制御システム | |
| CN107169221A (zh) | 动力电池测试方法及装置 | |
| JP6540424B2 (ja) | 推定装置、推定方法、推定プログラム、エンジンおよび移動装置 | |
| JP5691905B2 (ja) | インフレータ性能試験装置、方法及びプログラム | |
| Brown et al. | Surrogate modeling of manufacturing variation effects on unsteady interactions in a transonic turbine | |
| Chwalowski et al. | Fun3d analyses in support of the second aeroelastic prediction workshop | |
| GB2427046A (en) | Predicting telemetry data using coupled fluid mechanics and finite element analysis models | |
| US7925479B2 (en) | Efficient process for evaluating engine cooling airflow performance | |
| Pryor et al. | Development of transient simulation methodologies for underhood hot spot analysis of a truck | |
| CN114729608A (zh) | 用于校准机器控制器的方法和系统 | |
| Gurel et al. | Multi-objective optimization of transient air-fuel ratio limitation of a diesel engine using DoE based Pareto-optimal approach | |
| Saini et al. | The effect of blade count on body force model performance for axial fans | |
| Zieße et al. | Validation of a Synthetic-Eddy Method for Modelling Incoming Wakes in Scale-Resolving Simulations of Turbomachinery Cascades | |
| Burke et al. | Dynamic identification of thermodynamic parameters for turbocharger compressor models | |
| Salehi et al. | Effective component tuning in a diesel engine model using sensitivity analysis |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20140114 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20140825 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20140902 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20141003 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20150106 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20150119 |
|
| R151 | Written notification of patent or utility model registration |
Ref document number: 5691905 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |
|
| LAPS | Cancellation because of no payment of annual fees |