(57)【要約】      （修正有）【課題】火花点火機関の燃焼を制御する方法を提供する。【解決手段】機関制御システムは、燃焼室内のガスと燃料の混合物の燃焼に関連付けられた物理パラメータの値が、燃焼を最適化するように定められたそれらの設定値に等しくなるようにアクチュエータを制御する。設定値は、燃料混合物の点火クランク角について定められる。この設定値は次に、物理パラメータがその設定値に達する前に補正される。したがって、この点火角度設定値に適用される補正値は、クランク角ＣＡyがその設定値に等しくなるように算出される。最後に、機関制御システムは、最適な燃焼を維持するために、クランク角が補正設定値に等しいときに燃焼室内の混合物の点火を制御する。【選択図】図１
2010116917.tif 000002
【特許請求の範囲】
【請求項１】
火花点火機関の燃焼を制御する方法であって、
−燃焼室内のガスと燃料の混合物の燃焼に関連付けられた物理パラメータの設定値と、前記混合物の点火クランク角の設定値が求められ、前記設定値としては、燃焼を最適化するような値が求められ、
−前記物理パラメータの値が前記設定値に等しくなるように、機関制御システムがアクチュエータを制御する方法において、
−消費される燃料質量の経時的な変化を微分方程式によってモデル化する燃焼モデルによって、燃料のｙ％が燃焼中に消費されるクランク角ＣＡｙが最適化された燃焼についてのこの角度の設定値に等しくなるように、点火クランク角の前記設定値に適用すべき補正値ｄθallを算出することによって、前記物理パラメータがそれらの設定値に達する前に点火クランク角の設定値を補正することと、
−前記機関制御システムが、最適な燃焼を維持するために、前記クランク角が前記補正点火クランク角設定値に等しいときに前記燃焼室内の前記混合物の点火を制御することとを含む、ことを特徴とする方法。
【請求項２】
前記補正値ｄθallは、前記物理パラメータの実値ｐと前記物理パラメータの前記設定値ｐrefとの差ｄｐを考慮することによって求められる、請求項１に記載の方法。
【請求項３】
前記補正値ｄθallは、前記燃焼モデルを設定値ｐrefを中心とするｐに線形化し、次に、補正値ｄθallが前記差ｄｐに比例するように導くべき前記補正値の一次解を算出することによって求められる、請求項２に記載の方法。
【請求項４】
前記補正値ｄθallが、
−前記物理パラメータの前記実値を求める段階と、
−前記実値と前記設定値の前記差ｄｐを算出する段階と、
−前記モデルの各パラメータにその設定値を割り当てることにより前記燃焼モデルの数値積分法によってクランク角ＣＡｙの前記設定値を求める段階と、
−燃焼モデルを設定値ｐrefを中心とするｐに線形化することによって燃焼モデルの線形化行列Λを算出する段階と、
−以下の公式、すなわち、
【数１】
2010116917.tif 000016
によって補正値ｄθallを算出する段階とを適用することによって求められ、
（θall）refは前記混合物の前記点火クランク角の前記設定値であり、（ＣＡy）refはクランク角ＣＡｙの前記設定値である、請求項３に記載の方法。
【請求項５】
クランク角ＣＡｙは、燃焼中に前記燃料の５０％が消費されるクランク角である、請求項１から４のいずれか１項に記載の方法。
【請求項６】
前記物理パラメータは、弁が閉じられたときには少なくとも以下のパラメータ、すなわち、前記燃焼室内の圧力（ＰIVC）、前記燃焼室内の温度（ＴIVC）、前記燃焼室内の既燃ガス質量と全ガス質量との比（ＸIVC）、前記気筒内の空気質量（ＭIVC）、および吸気弁の閉角度（θivc）から選択される、請求項１から５のいずれか１項に記載の方法。
【請求項７】
前記物理パラメータがそれらの設定値に達する前に前記燃焼室内に噴射される燃料の質量も、前記リッチ燃焼を制御することによって適合させられる、請求項１から６のいずれか１項に記載の方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、機関制御の分野に関し、特に火花点火機関の燃焼制御に関する。
【背景技術】
【０００２】
（ガソリン式）火花点火機関の動作は、空気と既燃ガスと燃料との混合物の燃焼に基づく動作である。機関周期は４つの相に分類することができる（図１）。
−吸気相（ＡＤＭ）：吸気弁は、空気と既燃ガスの混合物をチャンバＣＨＢ内に入れるのを可能にする。空気は機関の外部環境から取り込まれる。既燃ガスは排気マニフォルドＥＣＨから取り込まれ、吸気マニフォルドに送り返され（排気ガス再循環ＥＧＲ）および／または排気弁によって吸い戻される（内部排気ガス再循環ｉＥＧＲ）。燃料は、吸気相中に噴射される。可変弁タイミング（ＶＶＴ）装置が、吸気弁リフトプロファイル（ＶＶＴａ）および排気弁リフトプロファイル（ＶＶＴｅ）に時間遅れを適用するのを可能にする。このことは、ガスの組成および燃焼室内の乱流に直接影響を与える。
−圧縮相：吸気弁が閉じた後（ＩＶＣ：吸気弁閉動作）、ピストンＰＩＳがガスを圧縮する。
−燃焼相：点火プラグＢＯＵＧが、空気と、既燃ガスと、燃料において得られる化学エネルギーを放出しながら発火し、したがって、ピストンを後方に押す過圧を発生させる燃料との混合物の燃焼を開始する火花を生じさせる。
−膨張相：ピストンが再び下降し始めるとすぐに、排気弁が開き、次にガス混合物が排気マニフォルドを通して排出される。
【０００３】
機関制御の目的は、騒音および汚染物質の放出を最小限に抑えつつ運転者に必要なトルクを供給することである。したがって、様々なガスおよび燃料の量の制御をできるだけ細かく調整する必要がある。
【０００４】
火花点火機関の燃焼制御を実施するにあたって、機関に取り付けられた検出器によって燃焼媒体を判定するのを可能にする公知の方法がある。以下によって行われる測定のような燃焼室内の直接的な測定に基づく方法がある。
−気筒圧力検出器：非特許文献１参照
−燃焼室内のイオン化電流検出器：非特許文献２参照
しかし、標準的な車両でこのような検出器を使用するのは、かなりのコストがかかるため困難である。さらに、これらの検出器は一般に、比較的高速のドリフトの影響を受けやすい。
【０００５】
各量およびタイミングが、各静的作用点（速度およびトルク）に対して各点で理想的な方法が実施されるように最適化される方法もある。したがって、以下の主要な３つのデータセットの最適値が得られるように機関試験のベンチ較正が行われる。
空気ループ：
−燃焼室で必要な空気Ｍairおよび既燃ガスＭgbの質量
−燃焼室内でのこれらのガスの圧力および温度
−可変弁リフト装置（ＶＶＴ）の位置、特にθivcで示される吸気弁の閉角度
これらの熱力学的変数および物理変数Ｘair＝（Ｍair，Ｍbg，Ｐ，Ｔ，θIVC）はＸairで代表される。
燃料ループ：
燃焼室内に噴射される（チャンバ内への直接的な噴射または吸気管内への間接的な噴射）燃料Ｍｆの質量、Ｘfuel＝（Ｍf）。
点火ループ：
Ｘall＝（θall）で表される、（プラグを介して）火花が生じるクランク角θall。
【０００６】
しかし、このような方法は遷移相では不十分である。すなわち、作用点から別の作用点までの遷移相（車両速度または路面形状の変化）の間、機関制御は、機関に存在する様々なアクチュエータを監視し、騒音、汚染物質の排出、および消費量を最小限に抑えつつ所望のトルクを確保する。したがって、これは、初期点の各パラメータの値から最終点の各パラメータの値への変化に変換される。
【０００７】
【数１】
2010116917.tif 000003
【０００８】
次に、機関には２つのタイムスケールがある。より高速のタイムスケール（５０Ｈｚ）は、燃焼現象全体（１機関周期）に相当する。このスケールでは、噴射（Ｘfuel）および点火（Ｘall）方法を変更して燃焼を制御することができる。より低速のタイムスケール（１Ｈｚ）は、機関マニフォルド（吸気、排気、既燃ガス再循環）内の気体力学およびアクチュエータ（ターボ圧縮機ＴＣ）の慣性に相当する。この空気ループ（Ｘair）の方法はより高速に変更することはできない。
【先行技術文献】
【非特許文献】
【０００９】
【非特許文献１】Paljoo Yoon et al., "Closed-loop Control of Spark Advance and Air-fuel Ratio in Sl Engines Using Cylinder Pressure", Society of Automotive Engineering World Congress, 2000-01-0933
【非特許文献２】Lars Eriksen et al., "Closed loop Ignition Control by Ionization Current Interpretation", SAE 1997 Transactions, Journal of Engines, Vol. 106, Section 3, pp. 1216-1223, 1997
【発明の概要】
【発明が解決しようとする課題】
【００１０】
したがって、現在の方法では、この力学の相違のために、各被制御変数（Ｘair，Ｘfuel，Ｘall）がその設定値に同時に達することはない。したがって、静止相ではトルクの生成、消費量、汚染物質、および騒音に関する目標が満たされ（動的ループがその基準値で安定化される）、もし遷移相では予防策が講じられないなら、パラメータの一部がほぼ瞬間的にその最終設定値に達する。一方で他の部分は依然として初期設定値である。したがって、機関は、より多くの汚染物質を放出するかまたはより大きい騒音を発し、場合によっては停止する。
【００１１】
さらに、気筒圧力検出器がない場合、公知の方法では、遷移相中に燃焼タイミングを制御することはできない。次に、図２および３に示されているように、これは、非定常条件の下で機関を適切に動作させるのに十分ではない。
【００１２】
本発明の目的は、従来技術の問題を解消しつつ、特に非定常条件の下での火花点火機関の燃焼を制御する方法に関する。この方法は、一方では３つの動的ループを別個に制御し、他方では、角度ＣＡ50の制御を介して点火角度の基準値を補正することによってこの目的を実現する。
【課題を解決するための手段】
【００１３】
本発明は、火花点火機関の燃焼を制御する方法であって、燃焼室内のガスと燃料の混合物の燃焼に関連付けられた物理パラメータの設定値、および混合物の点火クランク角の設定値（θall）refであって、燃焼を最適化するように定められた各設定値を求めることと、機関制御システムが、物理パラメータの値が設定値に等しくなるようにアクチュエータを制御することとを含む方法に関する。
【００１４】
この方法は、以下の段階、すなわち、
燃焼中に燃料のｙ％が消費されるクランク角ＣＡｙが、最適化された燃焼に関してこの角度の設定値に等しくなるように設定値（θall）refに適用すべき補正値ｄθallを算出することによって物理パラメータがその設定値に達する前に設定値（（θall）ref）を補正する段階と、
機関制御システムが、最適な燃焼を維持するために、クランク角が補正設定値（θall）refに等しいときに燃焼室内の混合物の点火を制御する段階とを含む。
【００１５】
本発明によれば、補正値ｄθallは、物理パラメータの実値ｐと物理パラメータの設定値ｐrefとの差ｄｐを考慮することによって求めることができる。したがって、消費された燃料の質量の経時的な変化をモデル化するのを可能にする微分方程式によって定義された燃焼モデルを使用し、この燃焼モデルを、設定値ｐrefを中心とするｐに線形化し、次に、補正値ｄθallが差ｄｐに比例するように導くべき補正値の一次解を算出することにより可能である。
【００１６】
したがって、一実施態様によれば、補正値ｄθallは、以下の段階、すなわち、
−物理パラメータの実値を求める段階と、
−実値と設定値の差ｄｐを算出する段階と、
−モデルの各パラメータにその設定値を割り当てることにより燃焼モデルの数値積分法によってクランク角ＣＡｙの設定値（ＣＡy）refを求める段階と、
−燃焼モデルを、設定値ｐrefを中心とするｐに線形化することによって燃焼モデルの線形化行列Λを算出する段階と、
−以下の公式、すなわち、
【００１７】
【数２】
2010116917.tif 000004
【００１８】
によって補正値ｄθallを算出する段階とを適用することによって求めることができる。
【００１９】
好ましい実施態様によれば、クランク角ＣＡｙは、燃焼中に燃料の５０％が消費されるクランク角である。
【００２０】
本発明によれば、物理パラメータは、弁が閉じられたときには少なくとも以下のパラメータ、すなわち、燃焼室内の圧力（ＰIVC）、燃焼室内の温度（ＴIVC）、燃焼室内の既燃ガス質量と全ガス質量との比（ＸIVC）、気筒内の空気質量（ＭIVC）、および吸気弁の閉角度（θivc）から選択することができる。
【００２１】
最後に、リッチ燃焼を制御することによって、物理パラメータがそれらの設定値に達する前に、燃焼室内に噴射される燃料の質量を適合させることも可能である。
【図面の簡単な説明】
【００２２】
【図１】花火点火機関の燃焼周期の様々な相を示す図である。
【図２】３つの燃焼制御状況、すなわち、（安定化された相で実施される）最適制御、ＣＡ50制御を伴わない遷移相における現在の制御、ＣＡ50制御を伴う遷移相における所望の制御によるクランク角の関数としての燃焼の経時変化を示す図である。
【図３】図２で説明した３つの状況についてのクランク角θの関数としての３本のエネルギー放出曲線Ｑを示す図である。
【図４】点火角度の補正値ｄθallの算出方式を示す図である。
【発明を実施するための形態】
【００２３】
本発明による方法は、静止相と遷移相において花火点火機関の燃焼進行を制御するのを可能にする。この方法は、空気ループと整合するように高速ループ力学を適合させることによって空気ループ（低速ループ）と燃料ループおよび点火ループ（高速ループ）とを別個にかつ独立に制御することを含む。したがって、この方法は、（運転者のトルク要求によって）要求された燃焼の特性を維持するようにＸfuelおよびＸallを適合させるのを可能にする。したがって、運転者に要求されたトルクを確保しつつ、排出物および騒音に対する影響が制限される。
【００２４】
この方法によれば、花火点火機関の燃焼の制御は５つの段階で実施される。
【００２５】
１−様々な物理パラメータの設定値を求める
ある作用点から別の作用点までの遷移相（車両速度または路面形状の変化）の間、機関制御は、機関に存在する様々なアクチュエータを監視し、騒音、汚染物質の放出、および消費量を最小限に抑えつつ所望のトルクを確保する。したがって、これは、初期点のパラメータＸair、Ｘfuel、およびＸallの値から最終点の各パラメータの値への変化に変換される。
【００２６】
【数３】
2010116917.tif 000005
【００２７】
騒音を最小限に抑えつつ排出および消費量を最小限に抑えるために、最終値は、燃焼を最適化し、すなわち、最大量の燃料を燃焼させるように定められる。燃焼を最適化するこれらの最終値は設定値と呼ばれる。機関制御は、これらの設定値を向上させるものである。
【００２８】
空気ループによって規制される重要な物理パラメータは、チャンバ内のガスの圧力、温度、化学組成と、吸気弁閉角度である。理想的には、これらのパラメータはその設定値に瞬間的に達する。実際には、空気ループが低速であるため、遷移相全体にわたってこれらのパラメータＸairに関する設定値と実値との誤差が生じる。したがって、気筒に吸入される供給ガスの熱力学的パラメータ（質量、圧力、温度、および既燃ガス流量）はその設定値とは異なる。燃料および点火ループ制御は、以下のパラメータに対する誤差に適合させられる。
【００２９】
−Ｐ：燃焼室内の圧力。この圧力はクランク角θによって決まる。
【００３０】
−Ｔ：燃焼室内の温度。この温度はクランク角θによって決まる。
【００３１】
−Ｘ：燃焼室内の既燃ガス質量と全ガス質量との比（０と１の間のパラメータ）。この比はクランク角θによって決まる。
【００３２】
−Ｍair：気筒内に閉じ込められた空気の質量
弁が閉じられたとき（ＩＶＣ）のこれらのパラメータの値は区別される。
【００３３】
−ＰIVC：弁が閉じられたときの燃焼室内の圧力
−ＴIVC：弁が閉じられたときの燃焼室内の温度
−ＸIVC：弁が閉じられたときの燃焼室内の既燃ガス質量と全ガス質量との比
−ＭIVC：弁が閉じられたときの気筒内の空気の質量
−θIVC：吸気弁の閉角度；閉角度は、燃焼室内の乱流に直接影響を及ぼす。
【００３４】
この５つのパラメータの値は連続的に求められる。したがって、弁が閉じられたときの気筒内の組成（ＸIVC）および圧力（ＰIVC）は、（検出器または推定器によって）測定値が得られる吸気マニフォルドの組成および圧力と同じであると仮定する。理想気体法則ＴIVC＝ＰIVCＶIVC／ＲＭIVCによってＴＩＶＣを推定する。この式で、Ｒは理想気体定数（Ｒ＝２８７）であり、ＭIVCは、流量計によって測定される、気筒によって吸入される質量である。
【００３５】
機関の燃焼室内の気体酸化性物質の吸入に関連付けられたこれらの５つのパラメータについて、設定値の値はそれぞれＰref、Ｔref、Ｘref、Ｍref、および（θivc）refによって示される。
【００３６】
これらの設定値は、機関の試験ベンチ上に確立される設定値マップから得られ、これらのパラメータの設定値は、試験ベンチ上にマップされる最適値（これらのパラメータが達しなければならない値）によって与えられる。これらの設定値は、燃焼を最適化するために求められる。これらのパラメータは、理想気体関係（ＰＶ＝ＭＲＴ）によって関係付けられるが、説明を簡単にするため、この関係については直接明示しない。このことが、提供される方法に何らかの影響を及ぼすことはない。
【００３７】
本発明によれば、制御されるパラメータは混合物点火角度：θallである。このパラメータの基準値（試験ベンチ上にマップされる最適値によって与えられる）は（θall）refによって表される。一定に維持されるパラメータはクランク角ＣＡｙ、すなわち、燃焼中に燃料のｙ％が消費される角度である。この角度は、最適な燃焼が行われるようにこの角度の設定値（ＣＡy）refに維持される。好ましい実施形態によれば、２分の１の燃焼角度ＣＡ50が使用される。この燃焼角度は、最適化された燃焼（設定値によって得られる燃焼）中に燃料の５０％が消費されたクランク角である。
【００３８】
２−空気ループ制御（低速ループ）
設定値Ｐref、Ｔref、Ｘref、Ｍref、および（θivc）refが求められると、機関制御システムは、物理パラメータＰIVC、ＴIVC、ＸIVC、およびθivcの値がその設定値Ｐref、Ｔref、Ｘref、Ｍref、および（θivc）refに等しくなるようにアクチュエータを制御する。
【００３９】
３−燃料ループ適合（高速ループ）
噴射される燃料質量の制御を空気ループの力学に適合させることは従来、リッチ燃焼を制御することによって実現されており、すなわち、ガソリン機関からの排気ガスの汚染除去は３元触媒によって行うことができる。３元触媒は、排気ガスが大域的に酸化性でも還元性でもないならば、燃焼によって生成されるＣＯ、ＨＣ、およびＮＯｘを効率的に処理するのを可能にする。リッチ燃焼（λ）は、化学量論的燃焼の場合（この化学量論比はＰＣＯによって示される）と同じ比にされた、燃料質量Ｍfに関して過剰空気質量Ｍairと定義される。したがって、λ＝Ｍf／Ｍair・１／ＰＣＯが得られる。したがって、汚染除去は、リッチネスが１に近いときに効率的になる。したがって、噴射される燃料質量の制御手順は、空気ループパラメータからの、気筒に吸入される空気質量の概算に要約される。この場合、空気質量Ｍairの概算は、コマンドＭf＝１／ＰＣＯ・Ｍairを適用することを可能とする。
【００４０】
４−点火ループ適合
ａ．点火角度設定値（θall）refの補正値の算出
点火角度θallについての従来の制御手順では、機関速度および（マッピングを介した）気筒内に吸入される空気質量の概算に応じて事前に位置決めが行われる。燃料質量制御とは異なり、この手順は最適ではない。すなわち、噴射される燃料質量がトルクポテンシャルを与える場合、このポテンシャルの有効利用を保証するのは点火タイミングである。
−点火が早すぎる場合、燃焼が起こるのも早くなり、圧縮相中に圧力が上昇し、トルクの発生を妨げる。
−点火が遅すぎる場合、燃焼が起こるのも遅くなり、膨張相全体がより低い圧力で生じる。この場合、発生するトルクが低くなる。
【００４１】
ＣＡ50（燃料の５０％が燃焼するクランク角）は、この燃焼タイミングを考慮するのを可能にするクランク角である。従来、各機関は、機関の技術的データに応じて一定のクランク角（ＣＡ50）refを有することができる。そして、ＣＡ50がその基準値（ＣＡ50）refに規制される場合に点火手順は最適なものになる。
【００４２】
したがって、点火角度を最適に制御するには、燃焼に影響を及ぼすすべての熱力学パラメータおよび物理パラメータを考慮に入れる必要がある。したがって、機関速度および吸入空気質量に応じて従来の点火角度の事前の位置決めは、空気ループの低速パラメータに対する点火のループの適応の最初であるが、完全なものではない。すなわち、空気ループの他のパラメータ、すなわち、
−燃焼室内の圧力
−燃焼室内の温度
−燃焼室内の既燃ガス流量
−燃焼室内の乱流（θivcを介する）
も燃焼に影響を与える。
【００４３】
空気ループの制御が完全である場合、５つのパラメータＰIVC、ＴIVC、ＸIVC、ＭIVC、およびθivcは瞬間的にそれら基準値Ｐref、Ｔref、Ｘref、Ｍref、および（θivc）refに達する。実際には、遷移相では、パラメータＰIVC、ＴIVC、ＸIVC、ＭIVC、およびθivcはそれらの基準値とは異なる。したがって、弁が閉じられたときの気筒の内容物は、点火手順がマッピングされた基準内容物とは異なる。
【００４４】
基準燃焼（最適化された燃焼）にできるだけ近い燃焼を維持するように点火角度を修正するために、弁が閉じられたときのこれらのパラメータの誤差を考慮に入れる必要がある。以下のパラメータが使用される。
−ｐベクトル（ＰIVC；ＴIVC；ＭIVC；ＸIVC；θivc）；
−ｐrefベクトル（Ｐref；Ｔref；Ｍref；Ｘref；（θivc）ref）；
−ｄｐベクトルｐ−ｐref＝（ｄＰ；ｄＴ；ｄＭ；ｄＸ；ｄθivc）；
この場合、ｄＰ＝ＰIVC−Ｐref、ｄＴ＝ＴIVC−Ｔref、ｄＭ＝ＭIVC−Ｍref、ｄＸ＝ＸIVC−Ｘref、およびｄθ＝θivc−（θivc）refが成立する。
【００４５】
したがって、角度ＣＡ50がその基準値（ｄＣＡ50＝ＣＡ50−（ＣＡ50）ref＝０）になるように新しい補正点火角度（θall）ｒｅｆ＋ｄθallが求められる。したがって、（３つの状況については図２および３を参照されたい）：
−誤差がない場合、すなわち、すべてのパラメータがそれらの基準値に達した場合（ｄｐ＝０）、ちょうど基準作用点の状況であり、したがって、ｄθall＝０が成立する（状況１）。
−パラメータがそれらの基準値に達していない場合（ｄｐ≠０）、燃焼速度は基準燃焼の燃焼速度と同一ではない。したがって、燃焼の位相ずれが生じ、角度ＣＡ50はそれらの基準値に達していない（状況２）。
−したがって、誤差ｄｐ≠０を相殺するために、同じ位相のＣＡ50が得られるように点火角度に角度補正値ｄθ≠０が導入される（状況３）。
【００４６】
図２は、３つの状況による燃焼の経時変化を示している。状況ごとに、水平軸はクランク角θを表す。これらの軸は、点火角度の設定値（θall）ref、点火角度θall、および補正項ｄθallを有している。図３は、上述の３つの状況（図２）についてのクランク角θの関数としての３本のエネルギー放出曲線Ｑを示している。
【００４７】
補正値ｄθallを求めるため、燃焼系のモデル化が行われる。特定の実施形態によれば、燃焼によって消費される燃料の質量の経時的な変化をモデル化するのを可能にする微分方程式によって定義された燃焼モデルを使用することが可能である。したがって、このような燃焼モデルは、以下のようにコンパクトな形式で記すことができる（付録１はこのような燃焼モデルを示している）。
【００４８】
【数４】
2010116917.tif 000006
【００４９】
上式で、
ｘ：既燃燃料の質量分率（ｘは次元１×１を有する）
ｙ：燃焼モデルに必要である力学の他の変数（圧力、温度、、、）。ｙは次元１×ｎを有し、この場合、ｎ∈である。
【００５０】
ｐ：遷移相中に補償すべき空気ループのパラメータ
ｐは次元１×ｎを有し、この場合、ｎ∈である。
【００５１】
θall：混合物点火角度
ｆ、ｇ、およびｈは全体的に公知の関数である（たとえば、付録１を参照されたい）。
【００５２】
本発明による点火角度制御方法が微分方程式形態で任意の燃焼モデルに適用できることに留意されたい。
点火角度設定値（θall）refの補正値の推定
補正値の計算は、差分ｄｐを導入することによって燃焼モデルを基準値ｐrefを中心とするｐに線形化することによって行われる。
【００５３】
燃焼モデルの複雑な形態を考えると、補正値ｄθallの分析式を求めることは困難である。したがって、導くべき補正値の一次解が求められる。このため、補正値は空気ループ誤差ｄｐに比例する。
【００５４】
補正値は以下のように得られる。
−２つの変数ｐおよびθallに対する微分方程式（１）の感度分析が行われる。したがって、それらの基準値を中心とするｐおよびθallのあらゆる変分に連続する微分方程式（ｄｘ，ｄｙ）の状態の一次変分が得られる。
−誤差ｄｐおよび補正値ｄθallが存在する場合、ＣＡ50が影響を受けるはずはないという事実（ｄＣＡ50＝０）が、状態の感度に対する条件ｘ：ｄｘ（θall）＝０に変換される。したがって、ｄｐおよびｄθallが現れる数式が得られる。この数式を反転することによって、ｄｐおよびモデルに存在するすべての変数の関数としてｄθallの数式が得られる。
【００５５】
得られる補正値は次式のように書かれる。
【００５６】
【数５】
2010116917.tif 000007
【００５７】
この場合、次式が成立する。
【００５８】
【数６】
2010116917.tif 000008
【００５９】
ｂ．点火ループ適合（高速ループ）
機関制御システムは、最適な燃焼を維持するために、クランク角が補正設定値（θall）ref＋ｄθallに等しいときに燃焼室内で燃料点火システムを駆動する。
【００６０】
この方法の１つの利点は、空気ループの誤差が、行列Λを介して点火コマンドに適用される補正値に直接関係付けられることである。補正値は、完全に計算可能であり、燃焼モデル、基準値Ｐref、Ｔref、Ｘref、およびＭref、ならびに一定の数の公知の定数のみに依存する。
【００６１】
上記の補正値を点火角度に適用することによって、角度ＣＡ50がその基準値になるようにすることができる（一次）。空気ループによって誤差ｄＰ、ｄＴ、ｄＭ、およびｄＸが徐々に零に近くなっていき、したがって、安定化された静止相では補正値が消える。この制御方法は、図４に概略的に示されている。この図は、点火角度の補正値ｄθallの計算方式を示している。パラメータｐの実値を推定または測定し（ＥＳＴ−ＡＣＴ）、これらのパラメータの設定値ｐrefおよび（θall）refを求めた（ＤＥＴ−ＣＯＮＳ）後、燃焼モデルの線形化行列が算出される（ＣＡＬ−ＭＡＴ）。次に、以下の係数（（ＣＡ50）ref−（θall）ref）Λが算出される。
【００６２】
したがって、（ＣＡ50）refの値（基準燃焼の燃焼角度の２分の１の値）が必要である。したがって、初期条件として基準値ｐrefおよび（θall）refを有する微分系（１）（燃焼モデル）が使用される。以下の系が得られる。
【００６３】
【数７】
2010116917.tif 000009
【００６４】
この系の数値積分法によって、ｘref＝０．５のときの角度θの値が求められる。この角度は（ＣＡ50）refの値に相当する。
【００６５】
最後に、補正値ｄθallが算出される（ＣＡＬ−ＣＯＲ）。
【００６６】
【数８】
2010116917.tif 000010
【００６７】
したがって、本発明による方法は、３つの動的ループを別個に制御し、点火角度の基準値を補正することによって、花火点火機関の燃焼制御を確実に行うことを可能にする。この補正値としては、角度ＣＡ50がその基準値になるように求められ、基準燃焼（最適化された燃焼）と同じエネルギー放出が結果として得られる。
燃焼モデル
微分方程式によって定義される燃焼モデルの一例が以下の文献に提示されている。
−F.-A. Lafossas, et al., "Application of a new 1d combustion model to gasoline transient engine operation," in Proc. SAE World Congress, no. 2005-01-2107, 2005
このモデルは、（無限に薄い層としてモデル化された）火炎前面によって分離された２つのゾーン（既燃ゾーンおよび未燃ゾーン）内の気筒の体積を表している。燃焼全体にわたって、火炎は既燃ゾーンから未燃ゾーンまで伝播する。このモデルの数式は以下のとおりである。
【００６８】
【数９】
2010116917.tif 000011
【００６９】
上式で、各記号は以下の数量を表す。
【００７０】
【表１】
2010116917.tif 000012
【００７１】
補償すべき空気ループパラメータはすべて、以下のように３つのパラメータでグループ化される。
【００７２】
【数１０】
2010116917.tif 000013
【００７３】
最後に、モデルのパラメータはＣ1、Ｃ2、Ｃ3、Ｃ4、αである。これらのパラメータの数値例（国際単位系ＳＩ）を以下の表に示す。
【００７４】
【表２】
2010116917.tif 000014
【００７５】
したがって、モデルは最終的に、以下のように簡略化された形態で書かれる。
【００７６】
【数１１】
2010116917.tif 000015
【００７７】
上式で、ｐ＝（ｐ1，ｐ2，ｐ3）TかつＣ＝（１，０，０，）Tであり、ｆおよびｇは、付録の最初の部分の燃焼モデルの数式によって定義される。
【００７８】
これらの数式で、項ＣＡ50は、任意の角度ＣＡyで容易に置き換えることができる。
【図１】
2010116917.tif 000017
【図２】
2010116917.tif 000018
【図３】
2010116917.tif 000019
【図４】
2010116917.tif 000020
【外国語明細書】
2010116917.tif 000021
2010116917.tif 000022
2010116917.tif 000023
2010116917.tif 000024
2010116917.tif 000025
2010116917.tif 000026
2010116917.tif 000027
2010116917.tif 000028
2010116917.tif 000029
2010116917.tif 000030
2010116917.tif 000031
2010116917.tif 000032
2010116917.tif 000033
2010116917.tif 000034
2010116917.tif 000035
2010116917.tif 000036
2010116917.tif 000037
2010116917.tif 000038
2010116917.tif 000039
2010116917.tif 000040
2010116917.tif 000041
2010116917.tif 000042
2010116917.tif 000043
2010116917.tif 000044
2010116917.tif 000045
2010116917.tif 000046
2010116917.tif 000047
2010116917.tif 000048
2010116917.tif 000049
2010116917.tif 000050
2010116917.tif 000051
2010116917.tif 000052