(57)【要約】  自動車両内燃機関のシリンダの燃焼室内に封入される質量を推定する方法であって、自動車両内燃機関が、シリンダ内の圧力を測定する手段、参照位置に対してのクランク軸の角度を測定する手段、および燃焼室冷却液の温度を測定する手段を備える。方法は、機関の回転速度がクランク軸角度の測定値から決定されるステップと、機関負荷が決定されるステップとを含む。方法はさらに、以下のステップ、すなわち、シリンダの弁が完全に閉じられるとき、および、燃料が最初に噴射される前に、シリンダの体積が、クランク軸の角度の関数として決定されるステップと、壁での熱損失に関連付けられる較正パラメータが、機関回転速度に依存的なマップから決定されるステップと、シリンダ内でモデリングされる温度が、燃焼室内の圧力および温度の非線形無次元モデルを使用して、シリンダ内の測定される圧力の、シリンダ内の気体に対する比気体定数の、シリンダ内の気体の一定体積での比熱容量の、燃焼室の体積の、および、壁での熱損失に関連付けられる較正パラメータの関数として決定されるステップと、シリンダ内に封入される質量が決定されるステップとを含む。【選択図】なし【特許請求の範囲】
【請求項１】
自動車両内燃機関のシリンダの燃焼室内に封入される質量を推定するための方法であって、前記自動車両内燃機関は、前記シリンダ内の圧力を測定するための手段、参照位置に対してのクランク軸角度を測定するための手段、および燃焼室冷却材の温度を測定するための手段を備え、
前記機関の回転速度が前記クランク軸角度の測定値から決定されるステップと、
機関負荷が決定されるステップと
を含む方法において、
以下のステップ、すなわち
前記シリンダの弁が全閉されるとき、および、燃料の最初の噴射の前に、
前記シリンダの体積が、前記クランク軸角度の関数として決定されるステップと、
壁での熱損失に関連付けられる較正パラメータが、前記機関の前記回転速度の関数であるマップの関数として決定されるステップと、
前記シリンダ内のモデリングされる温度が、前記燃焼室内の圧力および温度の非線形無次元モデルを使用して、前記シリンダ圧力の測定値、前記シリンダ内の気体に対する比気体定数、前記シリンダ内の前記気体の一定体積での比熱、前記燃焼室の体積、および、前記壁での前記熱損失に関連付けられる較正パラメータの関数として決定されるステップと、
前記のことから、前記シリンダ内に封入される前記質量に対する推定値が導出されるステップと
を含むことを特徴とする方法。
【請求項２】
前記非線形無次元モデルを求解するために、スライディングオブザーバが使用される、請求項１に記載の方法。
【請求項３】
前記壁での前記熱損失に関連付けられる前記較正パラメータのマップは、Ｗｏｓｃｈｎｉの相関を適用することにより、前記壁での前記熱損失の測定値の、シリンダ温度の測定値の、および、シリンダ冷却材温度の関数として、テストベッド上であらかじめ決定される、請求項１または２に記載の方法。
【請求項４】
前記シリンダ内に封入される前記質量は、前記シリンダ内の前記モデリングされる温度、前記シリンダ圧力の前記測定値、前記シリンダの前記体積、および、前記シリンダに収容される前記気体に対する前記比気体定数の関数として決定される、請求項１から３のいずれか一項に記載の方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の技術分野は、内燃機関管理である。
【背景技術】
【０００２】
汚染制御基準によって、自動車両製造者は、特には、汚染物質排出物、特に窒素酸化物ＮＯｘを推定および制御するために、軽油機関およびガソリン機関での燃焼の質をより良好に制御するための新しい戦略を開発することを余儀なくされている。
【０００３】
ＮＯｘを推定するための１つの解決策は、所定の車両についてふさわしい標準品である、シリンダ圧力センサを使用することである。このセンサを使用すると、シリンダ内の圧力を、あらゆる時点で知ることが可能である。残念ながら、窒素酸化物ＮＯｘの形成は、シリンダ内の温度による、次数１の関数として変動する。気体混合物の温度（Ｔ）は、シリンダ圧力（Ｐ）、および、シリンダ内に封入される質量（ｍ）が知られるならば、完全気体の状態方程式（Ｐ＊Ｖ＝ｍ＊ｒ＊Ｔ）を適用することにより、第１近似まで決定され得る。体積（Ｖ）は、クランク軸角度の幾何学的関数であり、Ｒは、完全気体定数である。
【０００４】
今は、燃焼室内に封入される質量に対するセンサは存在しない。
【０００５】
この封入される質量は従来では、充填マップを使用して推定されるが、これらのマップは、汚染物質排出物を推定することの目的に対して十分には正確でない。
【０００６】
さらには、排気ラインに適合させられ得る、汚染物質排出物センサ、特に、窒素酸化物ＮＯｘのセンサが存在する。しかしながら、そのようなセンサのコストは、この解決策が経済的に現実的ではないということを意味するものである。
【０００７】
さらには、文書ＷＯ２００９０４７４１２およびＦ２９４５３２０で説明されているように、シリンダ圧力Ｐｃｙｌセンサを伴う窒素酸化物ＮＯｘ推定器を使用することにより、ただし、封入される質量を推定することなく、窒素酸化物排出物を推定することが可能である。
【０００８】
あるいは、さらには、シリンダ圧力センサを伴うバージョンより高価でない、シリンダ圧力Ｐｃｙｌセンサを伴わないＮＯｘ推定器を使用することが可能である。そのような推定器は、例えば、文書ＷＯ２０１５／０５９０３４で説明されている。
【０００９】
しかしながら、両方の実例では、排出される窒素酸化物ＮＯｘの推定値は、今日の汚染制御基準を満たすのに十分には、確固たるもの、および正確ではない。
【発明の概要】
【００１０】
それゆえに、温度、および、窒素酸化物ＮＯｘの生み出される量を決定するように、シリンダ圧力を測定することにより、シリンダ内に封入される質量を決定するための方法に対する必要性が存在する。
【００１１】
本発明の主題は、自動車両内燃機関のシリンダの燃焼室内に封入される質量を推定するための方法であって、自動車両内燃機関は、シリンダ内の圧力を測定するための手段、参照位置に対してのクランク軸角度を測定するための手段、および燃焼室冷却材の温度を測定するための手段を備える。方法は、機関の回転速度がクランク軸角度の測定値から決定されるステップと、機関負荷が決定されるステップとを含む。
【００１２】
方法はさらには、以下のステップを含む。シリンダの弁が全閉されるとき、および、燃料の最初の噴射の前に、
シリンダの体積が、クランク軸角度の関数として決定されるステップと、
壁での熱損失に関連付けられる較正パラメータが、機関の回転速度の関数であるマップの関数として決定されるステップと、
シリンダ内のモデリングされる温度が、燃焼室内の圧力および温度の非線形無次元モデルを使用して、シリンダ圧力の測定値、シリンダ内の気体に対する比気体定数、シリンダ内の気体の一定体積での比熱、燃焼室の体積、および、壁での熱損失に関連付けられる較正パラメータの関数として決定されるステップと、
シリンダ内に封入される質量が決定されるステップ。
【００１３】
非線形無次元モデルを求解するために、スライディングオブザーバが使用されうる。
【００１４】
壁での熱損失に関連付けられる較正パラメータのマップは、Ｗｏｓｃｈｎｉの相関を適用することにより、壁での熱損失の測定値の、シリンダ温度の測定値の、および、シリンダ冷却材温度の関数として、テストベッド上であらかじめ決定されうる。
【００１５】
シリンダ内に封入される質量は、シリンダ内のモデリングされる温度、シリンダ圧力の測定値、シリンダの体積、および、シリンダに収容される気体に対する比気体定数の関数として決定されうる。
【００１６】
さらに、本発明の趣意、特徴、および利点は、単に非制限的な例ということで与えられる、後に続く説明を読むことから明らかになるであろう。
【発明を実施するための形態】
【００１７】
燃焼室内に封入される質量を推定するための方法は、燃焼室内の圧力および温度の動態を説明する０Ｄモデルに基づく。式の、後に続く系が、このモデルを例解する。
2019501324.tif 000002
2019501324.tif 000003
ここで、
ｐｃｙｌ：燃焼室内の圧力、
Ｔｃｙｌ：燃焼室内の温度、
ｒ：シリンダ内の気体に対する比気体定数、
ｃｖ：一定体積での比熱容量、
Ｖ：燃焼室の体積（クランク軸角度：シータの関数である幾何学的パラメータ）、
2019501324.tif 000004
：吸気に関連付けられるエンタルピーの変化のレート、
2019501324.tif 000005
：排気に関連付けられるエンタルピーの変化のレート、
2019501324.tif 000006
：燃焼に関連付けられるエンタルピーの変化のレート、
ｍｃｙｌ：燃焼室内の質量、および、
δＱｐａｒ：壁損失である。
【００１８】
封入される質量は、シリンダに、後に続く式により説明される完全気体法則を適用することにより、式（１）および（２）から推定され得る。
2019501324.tif 000007
【００１９】
それゆえに、封入される質量を推定することの問題は、式２での温度Ｔｃｙｌを、ｐｃｙｌの測定値から推定することの問題であると表現され得る。確認されうるように、式１および式２の系は、観測レベルで管理するのが困難な、きわめて複雑な系である。この理由で、機関の所定の段階のみが、温度を、推定することをより容易にするために考慮される。
【００２０】
それゆえに、後に続く想定が適用される。
・弁が全閉される間の時間の区間のみが、考慮に入れられる。
・質量は、弁の全閉と、燃料の最初の噴射との間に推定される。
【００２１】
これらの想定によって、エンタルピーの変化のレートに関連付けられる項を無視することが可能になり、そのことにより、温度を推定することの問題が、著しく単純化されることが可能となる。
【００２２】
それゆえに、式１および式２は、後に続くように再度定式化され得る。
2019501324.tif 000008
2019501324.tif 000009
【００２３】
並行して、壁損失δＱｐａｒは、次式により与えられるＷｏｓｃｈｎｉの近似を使用してモデリングされる。
δＱｐａｒ（ｔ）＝Ａｗｗｐｃｙｌ（ｔ）（ｋ１Ｔｃｙｌ（ｔ）−ｋ０Ｔｗ）  （式６）
ここで、
Ａｗ：伝熱面積、
ｗ：ｒａｄ／ｓでの機関の回転速度、
Ｔｗ：壁温度、
ｋ０およびｋ１：較正パラメータである。
【００２４】
壁温度Ｔｗは、冷却材温度と実質的に等しいと考えられるということに注目されたい。
【００２５】
Ｗｏｓｃｈｎｉの相関は、文献から知られている。式６により与えられる壁損失δＱｐａｒを、式４および式５に代入することにより、後に続く系を得る。
2019501324.tif 000010
2019501324.tif 000011
【００２６】
確認され得るように、式７および式８の系は、それでもなお、非線形性およびパラメトリック変動が原因で、求解するにはかなり複雑である。
【００２７】
それゆえに、提案は、燃焼室内の温度を推定するために、スライディングモード観測戦略（非線形戦略）を使用することであり、その戦略によって、封入される質量を推定することが可能になる。
【００２８】
そのことを行うために、燃焼室内の温度を推定するために後に続くようなスライディングモードオブザーバ（文献から知られている技法）が使用される。
2019501324.tif 000012
2019501324.tif 000013
上式で、
2019501324.tif 000014
であり、
2019501324.tif 000015
：燃焼室内の圧力に対する推定される値、
2019501324.tif 000016
：燃焼室内の温度の推定される値、ならびに、
λ１およびλ２：オブザーバパラメータであり、
モデリングされる温度
2019501324.tif 000017
は、後に続くように規定される。
2019501324.tif 000018
【００２９】
後に続く条件が設定される。
2019501324.tif 000019
【００３０】
式１３の条件が満たされるならば、温度の推定値についての誤差が、ゼロに向かって漸近的に収束し、このことによって、封入される質量が推定されることが可能となるということが論証され得る。
【００３１】
パラメータλ１およびλ２は、式１３の条件が充足される度に、オブザーバの収束のレートを調節するように較正され得る。
【００３２】
燃焼室内に封入される質量を推定するための方法は、後に続くステップを含む。
【００３３】
第１のステップ１の間、シリンダ圧力ｐｃｙｌ測定値、機関角度シータ、および、冷却材温度Ｔｗが決定される。
【００３４】
第２のステップの間、壁での熱損失に対するパラメータｋ０およびｋ１が、機関回転速度の関数としての、壁での熱損失のマップの関数として決定される。マップは、式６から、パラメータｋ０およびｋ１の較正により決定される。
【００３５】
第４のステップ４の間、シリンダ内のモデリングされる温度
2019501324.tif 000020
は、シリンダ内の測定される圧力ｐｃｙｌの、パラメータｋ０およびｋ１の、壁での熱損失の、ならびに、シリンダに関係のある構造パラメータの関数である、シリンダ内の圧力および温度の０Ｄモデルを使用して決定される。そのことを行うために、式７および式８が適用される。
【００３６】
第５のステップ５の間、シリンダ内に封入される質量は、式３を適用することにより、シリンダ内のモデリングされる温度
2019501324.tif 000021
の、シリンダ内の圧力ｐｃｙｌの測定値の、機関角度シータの関数としてのシリンダの体積Ｖ（ｔ）の、および、シリンダ内の気体に対する比気体定数ｒの関数として決定される。
【００３７】
この方策では、シリンダ内に封入される質量は、開ループで、すなわち、燃焼室内の圧力ｐｃｙｌの測定値についての、燃焼室内の温度Ｔｃｙｌのモデルの何らのループバックすることもなく推定される。そのことは、燃焼室内の温度の、およびそれゆえに、封入される質量の発散のおそれを生むので、代替案として、オブザーバパラメータλ１およびλ２に依存的なスライディングオブザーバを使用することにより、質量ｍｃｙｌの推定値をより確固たるものにすることが可能である。式９および式８が次いで、式７および式８の代わりに適用される。パラメータλ１およびλ２は、オブザーバの収束のレートを調節するような方策で決定される。収束のより高いレートはオブザーバの安定性を妨げるので、パラメータの選択は、収束の速度と、安定性との間の折り合いの結果である。
【国際調査報告】
2019501324.tif 000022
2019501324.tif 000023
2019501324.tif 000024
2019501324.tif 000025
2019501324.tif 000026
2019501324.tif 000027