(57)【要約】本発明は、自動車に装備された内燃エンジンのシリンダに吸入される空気質量の計算方法及び該計算方法を利用する噴射計算機に関する。本発明によれば、吸気マニホールドにおける圧力の推定値は、内燃エンジンの各シリンダについて、単数または複数の吸気弁の閉鎖時期（ｔ3）に対して、エンジンの動作を描写するパラメータの測定値の基づいて予測される。吸気マニホールドにおける圧力の予測値（Ｐres pred）は、本発明の対象である、独創的な吸気マニホールドのモデルを実現することから誘導される。
2004528504.tif 000002
【特許請求の範囲】
【請求項１】
シリンダへ噴射するべき燃料の量を決めるための、内燃エンジンのシリンダに吸入される空気質量の計算方法であって、上記エンジンは、複数の上記シリンダの中へ空気を導入する吸気マニホールドの中に設けられた圧力センサから供給される値に基づいて燃料噴射器の動作を制御する噴射計算機を含んでなり、吸入される上記空気質量は、上記吸気マニホールドにおける圧力Ｐcollに基づいて計算される、内燃エンジンのシリンダに吸入される空気質量の計算方法において、
上記シリンダに取り付けられた吸気弁の閉鎖の瞬間（ｔ3）に上記シリンダの中へ吸入される上記空気質量の予測計算を可能にするために、
−計算の瞬間における上記エンジンの実際の動作を表すパラメータ（alpha_pap、Ｎ、Ｔcoll、Ｐamont、ＰresＭes）を計測または推定し、次いで
−計算の瞬間における、スロットルバルブにおける空気流量（Ｄpap）と、上記エンジンの吸気行程における空気流量（Ｄmot）とを知るために、上記吸気マニホールドの動作モデルを算定し、次いで
−そこから上記吸気弁の閉鎖の瞬間の、上記吸気マニホールドにおける圧力の予測値（Ｐres pred）を得る、
ことの繰り返しを含んでなることを特徴とする、内燃エンジンのシリンダに吸入される空気質量の計算方法。
【請求項２】
サイクルiのときの、上記エンジンへの上記空気流量（Ｄmot）と上記スロットルバルブにおける上記空気流量（Ｄpap）の値の予測を、次式の形の上記吸気マニホールドの動作モデル：
【数１】
2004528504.tif 000013
ここに、τｉは、エンジン回転数Ｎと、上記吸気マニホールド及び上記シリンダの幾何学的形状と、上記エンジンの容積効率rempと、噴射計算の繰り返しとの関数；
Ｋiは、上記エンジン回転数Ｎと、上記エンジンの上記容積効率remp及び幾何学的形状と、上記吸気マニホールドにおける温度と、上記スロットルバルブの断面との関数；
ｆｂｓｖは、上記スロットルバルブの流量係数を定義するための関数発生器によって予め定義された関数；
から得られた、上記吸気マニホールドにおける圧力（Ｐcoll）と上記吸気マニホールドの上流における圧力（Ｐamont）との比を表す状態変数（Ｘｉ）に基づいて実行することを特徴とする、請求項１に記載の、内燃エンジンのシリンダに吸入される空気質量の計算方法。
【請求項３】
上記スロットルバルブの上記流量係数を定義するために、Ｘiの小さい値についての水平な第１の区間と、Ｘｉのより大きい値についてのほぼ垂直な第３の区間と、Ｘｉの中間の値についての単調に減少する第２の区間とによって、上記状態変数Ｘｉの関数としてグラフで表示された上記関数ｆｂｓｖを使用することを特徴とする、請求項２に記載の、内燃エンジンのシリンダに吸入される空気質量の計算方法。
【請求項４】
次式：
【数２】
2004528504.tif 000014
を用いて上記状態変数（Ｘｉ）の値を決めることが可能になるように、実際の上記関数ｆｂｓｖを表す曲線を近似する線分の勾配（pente）と原点の縦座標（Ｙ0）を決めて、上記スロットルバルブの上記流量係数を近似的に定義する、直線化された関数ｆｒｅｄを生成するために、先行するサイクルで決められた状態変数（Ｘｉ−１）の値を使用することを特徴とする、請求項３に記載の、内燃エンジンのシリンダに吸入される空気質量の計算方法。
【請求項５】
実際の上記関数ｆｂｓｖの代わりに、点（Ｘｉ，ｆｒｅｄ（Ｘｉ））を点（Ｘｉ，ｖｒａｉ，ｆｂｓｖ（Ｘｉ，ｖｒａｉ））に関連づける関係に従って、上記直線化された関数ｆｒｅｄによって計算された値Ｘｉに対応する、状態変数の真の値（Ｘｉ，ｖｒａｉ）を求めることを特徴とする、請求項４に記載の、内燃エンジンのシリンダに吸入される空気質量の計算方法。
【請求項６】
上記吸気マニホールドのモデルは、上記エンジンの動作を描写する：
−上記スロットルバルブの角度センサ（１）によって測定された上記スロットルバルブの開口角度alpha_pap；
−エンジン回転速度すなわちエンジン回転数のセンサ（２）によって測定された、上記エンジン回転数Ｎすなわち上記エンジンの回転速度；
−上記吸気マニホールドに設けられた温度センサ（３）によって測定された、上記吸気マニホールドにおける空気温度Ｔcoll；
−圧力センサ（５）によって測定、または推定された、上記スロットルバルブの上流における上記圧力Ｐamont；
−圧力センサ（６）によって測定された、上記吸気マニホールドにおける測定圧力（Ｐres Ｍes）；
の各パラメータを入力されることを特徴とする、請求項５に記載の、内燃エンジンのシリンダに吸入される空気質量の計算方法。
【請求項７】
上記吸気マニホールドにおける圧力の予測値から推定された上記空気質量を予測する段階をさらに含むことを特徴とする、請求項１〜６のいずれか１つに記載の、内燃エンジンのシリンダに吸入される空気質量の計算方法。
【請求項８】
少なくとも１つの上記シリンダと上記燃料噴射器を含んでなる上記内燃エンジンの動作を制御するために、請求項１〜７のいずれか１つに記載の内燃エンジンのシリンダに吸入される空気質量の計算方法を実行する噴射計算機において、
−モデル化された上記吸気マニホールドにおける上記圧力Ｐcolの、計算の瞬間における上流における上記圧力に対する比を表す上記状態変数（Ｘｉ）を生じる、考慮対象の上記シリンダの上記吸気マニホールドのモデルを実現するモジュール（７）；
−上記直線化された関数（ｆｒｅｄ）と、上記スロットルバルブの上記空気流量係数の実際の関数（ｆｂｓｖ）との間の関係から、上記状態変数（Ｘｉ）を修正するモジュール（８）；
−上記スロットルバルブにおける上記空気流量（Ｄpap）と考慮対象の上記シリンダへの上記吸入マニホールドにおける上記空気流量（Ｄmot）とを計算するモジュール（９）；
−考慮対象の上記シリンダの上記吸気弁の閉鎖の瞬間の、上記吸気マニホールドにおける上記圧力の値（Ｐres pred）を予測するモジュール（１０）；
を含んでなることを特徴とする噴射計算機。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、自動車に装備された内燃エンジンのシリンダに吸入される空気質量の計算方法及び該計算方法を利用する噴射計算機に関する。
【背景技術】
【０００２】
自動車を駆動する熱エンジンのための噴射計算機に主として適用される計算方法が、既に現在の技術において開示されている。特に、本出願人の名において出願された、フランス国特許公報ＦＲ−Ａ−２．７０９．１５１を参照する。
この現在の技術には、シリンダの中における最良の燃焼に必要な空気質量を如何にして予想するかが示されている。そのために、空気質量を、時間ｔにおける吸気マニホールドにおける圧力Ｐcolの測定に基づいて、予想継続時間Δｔの間に、次式によって予想するべきであることが明確にされた。
【０００３】
【数３】
2004528504.tif 000003
【特許文献１】
フランス国特許公報ＦＲ−Ａ−２．７０９．１５１
【発明の開示】
【発明が解決しようとする課題】
【０００４】
特許公報ＦＲ−Ａ−２．７０９．１５１に明示された方法においては、ステップiにおけるｄＰcol／ｄｔを計算するために、先行する計算段階（i−１）におけるＰcolのモデル化された値が用いられている。
その結果、負荷が強いと（Ｐcolは、上流における圧力の値に近い）、計算段階i−１においてモデル化された吸気マニホールドにおける圧力は、計算段階iのためのモデルを計算する圧力と異なり、それゆえモデルは振動し更に発散するようになるという問題がある。
換言すれば、特許公報ＦＲ−Ａ−２．７０９．１５１に示された方法は、熱エンジンが低負荷の回転数で作動するときには正確で、熱エンジンが最大負荷の近辺で作動するときには、ますますデリケートな、実時間での修正措置を必要とするという問題がある。
【課題を解決するための手段】
【０００５】
本発明は、従来技術における上記の問題を解決するためのものであって、
シリンダへ噴射するべき燃料の量を決めるための、内燃エンジンのシリンダに吸入される空気質量の計算方法であって、上記エンジンは、複数の上記シリンダの中へ空気を導入する吸気マニホールドの中に設けられた圧力センサから供給される値に基づいて燃料噴射器の動作を制御する噴射計算機を含んでなり、吸入される上記空気質量は、上記吸気マニホールドにおける測定された圧力Ｐcollに基づいて計算される、内燃エンジンのシリンダに吸入される空気質量の計算方法において、
上記シリンダに取り付けられた吸気弁の閉鎖の瞬間に上記シリンダの中へ吸入される上記空気質量の予測計算を可能にするために、
−計算の瞬間における上記エンジンの実際の動作を表すパラメータ（alpha_pap、Ｎ、Ｔcoll、Ｐamont、ＰresＭes）であって、それらのいくつかのものは、測定された大きさに対する測定の遅れを含むパラメータを計測または推定し、次いで、
−考慮対象の計算の瞬間における、スロットルバルブにおける空気流量と、エンジンの吸気行程における空気流量とを知るために、上記吸気マニホールドの動作モデルを算定し、次いで
−そこから上記吸気弁の閉鎖の瞬間の、上記吸気マニホールドにおける圧力の予測値を得る、
ことの繰り返しを含んでなることを特徴とする、内燃エンジンのシリンダに吸入される空気質量の計算方法に関する。
【０００６】
その他の特徴によれば、本発明の内燃エンジンのシリンダに吸入される空気質量の計算方法は、
サイクルiのときの、上記エンジンへの上記空気流量と上記スロットルバルブにおける上記空気流量の値の予測を、次式の形の吸気マニホールドの動作モデル：
【０００７】
【数４】
2004528504.tif 000004
ここに、τｉは、エンジン回転数Ｎと、上記吸気マニホールド及び上記シリンダの幾何学的形状と、上記エンジンの容積効率rempと、噴射計算の繰り返しとの関数；
Ｋiは、上記エンジン回転数Ｎと、上記エンジンの上記容積効率remp及び幾何学的形状と、上記吸気マニホールドにおける温度と、上記スロットルバルブの断面との関数；
ｆｂｓｖは、上記スロットルバルブの流量係数を定義するための関数発生器によって予め定義された関数；
から得られる、上記吸気マニホールドにおける圧力と上記吸気マニホールドの上流における圧力との比を表す状態変数に基づいて実行する。
【０００８】
その他の特徴によれば、本発明の内燃エンジンのシリンダに吸入される空気質量の計算方法は、上記スロットルバルブの上記流量係数を定義するために、Ｘiの小さい値についての水平な第１の区間と、Ｘｉのより大きい値についてのほぼ垂直な第３の区間と、Ｘｉの中間の値についての単調に減少する第２の区間とによって、上記状態変数Ｘｉの関数としてグラフで表示された上記関数ｆｂｓｖを使用する。
【０００９】
その他の特徴によれば、本発明の内燃エンジンのシリンダに吸入される空気質量の計算方法は、
次式：
【００１０】
【数５】
2004528504.tif 000005
を用いて上記状態変数Ｘｉの値を決めることが可能になるように、実際の上記関数ｆｂｓｖを表す曲線を近似する線分の勾配（pente）と原点の縦座標（Ｙ0）を決めて、上記スロットルバルブの上記流量係数を近似的に定義する、直線化された関数ｆｒｅｄを生成するために、先行するサイクルで決められた状態変数の値を使用する。
【００１１】
その他の特徴によれば、本発明の内燃エンジンのシリンダに吸入される空気質量の計算方法は、実際の上記関数ｆｂｓｖの代わりに、点（Ｘｉ，ｆｒｅｄ（Ｘｉ））を点（Ｘｉ，ｖｒａｉ，ｆｂｓｖ（Ｘｉ，ｖｒａｉ））に関連づける関係に従って、上記直線化された関数ｆｒｅｄによって計算された値Ｘｉに対応する、状態変数の真の値（Ｘｉ，ｖｒａｉ）を求める。
【００１２】
その他の特徴によれば、本発明の内燃エンジンのシリンダに吸入される空気質量の計算方法によって、上記吸気マニホールドは、上記エンジンの動作を描写する：
−上記スロットルバルブの角度センサによって測定された上記スロットルバルブの開口角度alpha_pap；
−エンジン回転速度すなわちエンジン回転数のセンサによって測定された、上記エンジン回転数Ｎすなわちエンジンの回転速度；
−上記吸気マニホールドに設けられた温度センサによって測定された、上記吸気マニホールドにおける空気温度Ｔcoll；
−圧力センサによって測定、または推定された、上記スロットルバルブの上流における上記圧力Ｐamont；
−圧力センサによって測定された、上記吸気マニホールドにおける測定圧力（Ｐres Ｍes）；
の各パラメータによってモデル化される。
【００１３】
その他の特徴によれば、本発明の内燃エンジンのシリンダに吸入される空気質量の計算方法は、上記吸気マニホールドにおける圧力の予測から推定された上記空気質量を予測する段階をさらに含む。
【００１４】
本発明は、本発明による内燃エンジンのシリンダに吸入される空気質量の計算方法を実行することを特徴とする噴射計算機にも関する。
【００１５】
その他の特徴によれば、本発明の噴射計算機は、
−モデル化された上記吸気マニホールドにおける上記圧力の、計算の瞬間における上流における上記圧力に対する比を表す上記状態変数（Ｘｉ）を生じる、上記吸気マニホールドのモデルを実現するモジュール；
−上記直線化された上記関数（ｆｒｅｄ）と、上記スロットルバルブの上記空気流量係数の実際の関数（ｆｂｓｖ）との間の関係から、上記状態変数（Ｘｉ）を修正するモジュール；
−上記スロットルバルブにおける上記空気流量（Ｄpap）と上記シリンダへの上記空気流量（Ｄmot）とを計算するモジュール；
−考慮対象の上記シリンダの上記吸気弁の閉鎖の瞬間の、上記吸気マニホールドにおける上記圧力の値（Ｐres pred）を予測するモジュール；
を含んでなる。
【発明を実施するための最良の形態】
【００１６】
本発明のその他の特徴及び利点は、理解のために添付図面を参照する、以下の詳細な記述を読むことによって明らかとなるであろう。これらの図において：
−図１は、本発明が根拠とする原理を説明するグラフ；
−図２は、本発明の方法を実行する装置のブロック図；
−図３は、本発明の方法によってもたらされた改良を説明するグラフ；
である。
【００１７】
図１に、本発明が根拠とする原理を説明するグラフを示す。この技術についての補足的な説明として、本出願人の名において交付された、フランス国特許公報ＦＲ  ２．７０９．１５１を参照する。
燃エンジンを制御するために用いられる燃料噴射技術においては、シリンダとその弁の動作を考察する。制御下にあるシリンダの吸気弁の閉鎖の瞬間における空気圧力の推定は、シリンダへ持ち込まれる空気質量の値を知ることと、それから最適の燃焼のために噴射するガソリンの量を導きだすこととを可能にする。しかしながら、間接噴射においては、気化という物理現象のために、このガソリンの量は、吸気弁の開口より前に導入されねばならず、また、あらゆる噴射タイプにおいても、測定には大きな遅延が伴い、測定に満足できるような吸気マニホールドにおける圧力の最新の値を知ることはできない。
【００１８】
図１に、時間ｔ1と時間ｔ3との間における吸気マニホールドにおける圧力Ｐcollの変化を示した。時間ｔ1において、吸気マニホールドに設けられた圧力センサのような適当な手段によって、吸気マニホールドにおける圧力Ｐcoll（ｔ1）の測定を実施した。
噴射計算の時間である時間ｔ2において、吸気マニホールドの物理モデルを用いて、吸気マニホールドにおける圧力の瞬間的な変化を推定する。この物理モデルは、吸気弁閉鎖の時間ｔ3における吸気マニホールドにおける圧力の予測を、時間間隔Δｔ＝ｔ3−ｔ1として、Ｐcollの事前の測定からのＰcollの予測に該当する次式によって可能にする。
【００１９】
【数６】
2004528504.tif 000006
図１のグラフを参照すると、時間ｔ1において測定された吸気マニホールドにおける圧力Ｐcollの値と、増分または減分ΔＰcollとの合計に相当する圧力Ｐcoll predは、吸気マニホールドにおける圧力の瞬間的な変化の実際の曲線の点に一致することができないことが認められる。
【００２０】
図２に、本発明の方法を実行する噴射計算機の種々の構成要素を示すブロック図を示す。本発明の方法においては、吸気弁の閉鎖時期における吸気マニホールドにおける圧力の予測値Ｐcoll predは、計算時に利用可能な測定データに基づいて、次ぎの噴射を準備するために実行される。下記のパラメータ（alpha_pap、Ｎ、Ｔcoll、Ｐamont、ＰresＭes）は、計算の瞬間におけるエンジンの実際の動作を表す。これらのパラメータは：
−スロットルバルブの角度センサ１によって測定されたスロットルバルブの開口角度alpha_pap；
−エンジン回転速度すなわちエンジン回転数のセンサ２によって測定された、エンジン回転数Ｎすなわちエンジンの回転速度；
−吸気マニホールドに設けられた温度センサ３によって測定された、吸気マニホールドにおける空気温度Ｔcoll；
−圧力センサ５によって測定、または上流における空気圧力の他のモデルに基づいて推定された、スロットルバルブの上流における空気の圧力Ｐamont；
−圧力センサ６によって測定された、ここではＰresＭesと名付けられた、吸気マニホールドにおける測定圧力。この測定は、無視できない測定時定数を有し、ある時間以前のエンジンの動作を表す；
である。
【００２１】
最初の４つのデータは、コンピュータからなるモジュール７へ供給される。このコンピュータの中では、本発明によって定義された吸気マニホールドの物理モデルがプログラムされ、コンピュータは、制御されるエンジンに適応化される。本発明において使用される吸気マニホールドのモデルは、スロットルバルブにおける空気流量Ｄpapを有して上流から流入する空気流量によって満たされ、シリンダへの空気流量Ｄmotを有して下流から流出する空気流量によって空になる、吸気マニホールドに典型的な体積の原理に基づいている。
【００２２】
１実施の形態においては、モジュール７は、次式によって定義されるスロットルバルブにおける空気流量Ｄpapを表す関数を生成するための手段を含んでなる。
【００２３】
【数７】
2004528504.tif 000007
この関数は、考慮されるモデルに関して、Ｐcollに本質的に依存し、この関数の中で、ｇは、スロットルバルブの空気力学的な作用を表す関数であり、ｆｂｓｖは、後に説明する図３に示される曲線に応じて定義される関数発生器によって、モジュール７の中で定義される関数である。関数ｇは、スロットルバルブの断面を表すＳpapと、それぞれスロットルバルブの上流における圧力と温度を表すＰamontとＴamontとのパラメータに依存する。これらのデータは、モジュール７の中に記録されているか、入力パラメータの検出モジュール１〜５の１つから受けられる。
【００２４】
一般に、図３にも示すように、スロットルバルブからの流量の流量係数を定義するための関数ｆｂｓｖは、吸気マニホールドにおける圧力が上流における圧力に関係付けされる比Ｐcoll／Ｐamontを表す状態変数Ｘｉに依存する。関数ｆｂｓｖのグラフ化された表示は、Ｘｉの小さい値についての水平な第１の区間と、Ｘｉのより大きい値についてのほぼ垂直な第３の区間と、Ｘｉの中間の値についての単調に減少する第２の区間とを含む。
【００２５】
１実施の形態においては、モジュール７は、次式によって定義されるエンジンのシリンダへの空気流量Ｄmotを表す関数を生成するための手段を含んでなる。
【００２６】
【数８】
2004528504.tif 000008
この関数は、考慮されるモデルに関して、Ｐcollに本質的に依存し、この関数の中で、ｈは、関数ｈの発生器（図示しない）によって作られ、入力パラメータの検出モジュール２、３からパラメータＮ、Ｔcollを入力として受ける、吸気中のシリンダの空気力学的な作用を表す関数であり、Ｐ0は、エンジンの考慮対象のシリンダへの空気流量を保証する吸気マニホールドの最低圧力である。関数ｈは、噴射計算機の作用対象である熱エンジンの、容積効率に特有の係数rempにも依存する。
【００２７】
吸気マニホールドのモデルは、モデルから導きだされた微分方程式と上述の関数ｆ、ｇ、ｈによって、モジュール７の中へインプリメーションされ、考慮対象のシリンダに対して、予想進行中のサイクル番号iについて、次式によって定義される。
【００２８】
【数９】
2004528504.tif 000009
ここに、
Ｘｉは、サイクル番号ｉにおける、モデル化された吸気マニホールドにおける圧力の、計測または推定された上流における圧力に対する比；
τｉは、エンジン回転数Ｎと、吸気マニホールド及びシリンダの幾何学的形状と、エンジンの容積効率（透過性）rempと、噴射計算の繰り返しとの関数；
Ｋｉは、エンジン回転数Ｎと、エンジンの容積効率（透過性）remp及び幾何学的形状と、吸気マニホールドにおける温度と、スロットルバルブの断面との関数；
である。
【００２９】
関数ｆｂｓｖに特有の形のために、この関数の逆関数を得ることはできないので、上の関係式によって、実時間で解析的に、ＸｉをＸｉ−１の関数として導きだすことはできない。
本発明に従ってこの問題を解くために、関数ｆｂｓｖを、予め定められた一連の線分を介して関数ｆｂｓｖに近づく関数ｆｒｅｄに置き換える。このため、計算のこの段階において、関数ｆｂｓｖの断片による線形化を実行する。図３に示された関数ｆｒｅｄの各線分は、次式を有する。
【００３０】
【数１０】
2004528504.tif 000010
ここに、「pente」は、前の予測の際に得られ、横座標の点（Ｘｉ−１）を用いて特定された、関数ｆｒｅｄの線分の勾配であり、Ｙ0は、この線分を有する直線の原点の縦座標である。これらの値は、直線化された関数ｆｒｅｄの発生手段（図示しない）の中の表に記載することができる。従って、モデルを次式の形のＸｉによって解くことができる。
【００３１】
【数１１】
2004528504.tif 000011
ここに、penteとＹ0は、本発明の近似法において、関数ｆｂｓｖを置き換える直線化された関数ｆｒｅｄに特有のものである。
【００３２】
従って、１実施の形態においては、吸気マニホールドのモデル化のモジュール７は、図１の予測の時間ｔ2における状態変数Ｘｉを発生し、状態変数Ｘｉは、時間ｔ3における吸気マニホールドの予想される圧力の計算手段であるモジュール８へ伝達される。
【００３３】
モジュール８は、横座標の点Ｘｉが最大値Ｘ＝１．０（図３を見よ）の極めて近傍にある時に、状態変数Ｘｉの修正を可能にする、修正手段を含んでなる。実際には、直線化された関数ｆｒｅｄは、ほぼ垂直な単数または複数の線分上でモデルがループ化することを避けるように、理論上の最大の横座標（Ｘ＝１）の近傍における現実の関数ｆｂｓｖからかなり遠ざかる。関数ｆｒｅｄの最後の線分は、従って、無限大でない勾配を有する。この最後の線分上におけるモデルの最終段階は、従って修正を要する不正確なＸｉ及びｆｒｅｄ（Ｘｉ）の値に導く。
【００３４】
本発明の方法においては、段階（Ｅ1）のときに、変数Ｘｉ−１の値から、直線化された関数ｆｒｅｄに関する線分（従って、直線化された関数ｆｒｅｄの「pente」及びＹ0の値）を選ぶ。
【００３５】
次に、段階（Ｅ2）のときに、上記の吸気マニホールドのモデルを実行するモジュール７は、先に取得された値Ｘｉ−１から値Ｘｉを発生する。値Ｘｉは、直線化された関数ｆｒｅｄを表す直線上の点Ｍｉに対応する。次ぎに、吸気マニホールドのモデルのループ（モジュール７）の外で、先行する直線に近い実際の関数ｆｂｓｖを表す曲線に一致する点Ｍｉ'を得る。そのため、修正用のモジュール８（図２）は、段階（Ｅ3）（図３）のときに、曲線上の点［Ｘｉ，ｆｒｅｄ（Ｘｉ）］を、真の曲線ｆｂｓｖ上の対応する点［Ｘｉ，ｖｒａｉ，ｆｂｓｖ（Ｘｉ，ｖｒａｉ）］に関連づける表に示された関数を用いて、状態変数Ｘｉの真の値と、実際の関数ｆｂｓｖ（Ｘｉ）の値とを発生することを可能にする。
【００３６】
そこで、真のＸｉ，ｖｒａｉと真のｆｂｓｖ（Ｘｉ，ｖｒａｉ）の値は、エンジンへの空気流量Ｄmotとスロットルバルブにおける空気流量Ｄpapの計算モジュール９へ伝達される。計算モジュール９は、それぞれ予め記録された上述の関数ｇおよびｈを実行する、Ｄmot関数発生器とＤpap関数発生器とを含んでなる。これらの２つの関数発生器は、計算された真のＸｉと真のｆｂｓｖ（Ｘｉ，ｖｒａｉ）の値、並びに、センサ５によって計測または推定された、スロットルバルブの上流における圧力Ｐamontの値を入力として受ける。これらの値は、吸気マニホールドにおける圧力を予測する時間ｔ2において利用可能な値である。
【００３７】
計算モジュール９によって発生された、エンジンへの空気流量Ｄmotとスロットルバルブにおける空気流量Ｄpapの予測値は、そこで最後の、吸気マニホールドにおける圧力の予測値の計算モジュール１０の入力へ伝達される。計算モジュール１０は、関数発生器によって、関数ＧＰ（Ｄmot，Ｄpap，ＰresＭes）の計算を実行する。ＰresＭesの値は、図１のグラフのｔ2の瞬間に測定され、ｔ1の瞬間の値に相当する、吸気マニホールドにおける圧力のために利用可能な、最新の測定値に該当する。
【００３８】
１実施の形態においては、関数ＧＰは、関数ｇｐ（）を、モジュール９によって発生された、エンジンへの空気流量Ｄmotとスロットルバルブにおける空気流量Ｄpaphaの予測値の、予め定められた関数として、ＧＰ＝ＰresＭes＋ｇｐ（Ｄmot，Ｄpap）の形で書き表すことができる。１実施の形態においては、関数ｇｐ（）は、次式によって表すことができる。
【００３９】
【数１２】
2004528504.tif 000012
ここに、ｒは乗算係数、ＴcolとＶcolは吸気マニホールドにおける温度と体積の状態、Δｔは、積分時間すなわちモデル化時間である。
【００４０】
計算モジュール１０の出力の値は、時間ｔ3(図1)における吸気マニホールドにおける圧力の予測値を表し、従って、噴射計算機の残りの部分（図示しない）へ伝達されたこの値は、エンジンの各シリンダに対する吸入空気の質量の改良された計算を可能にする。
【００４１】
本発明は、特に、直線化された関数ｆｒｅｄの直線性によってモデルの安定性が獲得され、安定性に関して妥協することがないので、従来技術のモデルに対してより良好な動力学を確保することを可能にする。
【００４２】
本発明は、負荷（吸気マニホールドにおける圧力）にかかわらないモデルの安定性によって容易にされる、噴射計算機の調整の実行を可能にする。
【００４３】
本発明は、従来技術の解決策に対して、測定された物理的な現実からの吸気マニホールドのモデルの遅れを減らすことを可能にし、このことが時間ｔ3における吸気マニホールドにおける圧力の予測の良好な妥当性を保証する。
【００４４】
最後に、圧力の予測、従って噴射の正確さは、強い負荷において改良される。従来技術の噴射計算機においては、最適化手段から生じた特有で非物理的なパラメータは、従来の解決法の安定性の拘束を考慮に入れたのに対して、本発明は、熱エンジンの物理的な大きさに対応する、測定された、または噴射計算器によって実行される他の計算で計算可能な、既知の調整パラメータの利用を可能にする。
【図面の簡単な説明】
【００４５】
【図１】本発明が根拠とする原理を説明するグラフである。
【図２】本発明の方法を実行する装置のブロック図である。
【図３】本発明の方法によってもたらされた改良を説明するグラフである。
【図２】
2004528504.tif 000015
【国際公開パンフレット】
2004528504.tif 000016
2004528504.tif 000017
2004528504.tif 000018
2004528504.tif 000019
2004528504.tif 000020
2004528504.tif 000021
2004528504.tif 000022
2004528504.tif 000023
2004528504.tif 000024
2004528504.tif 000025
2004528504.tif 000026
2004528504.tif 000027
2004528504.tif 000028
2004528504.tif 000029
2004528504.tif 000030
2004528504.tif 000031
2004528504.tif 000032
2004528504.tif 000033
2004528504.tif 000034
【国際調査報告】
2004528504.tif 000035
2004528504.tif 000036
2004528504.tif 000037
2004528504.tif 000038
2004528504.tif 000039
2004528504.tif 000040