(57)【特許請求の範囲】
【請求項１】
エンジンのクランク軸の回転変動波形を回転センサで検出して測定データを得る回転変動検出ステップと、回転変動波形データにおける該検出された測定データから基準データを算出する基準作成ステップと、該エンジンの各気筒の上死点近傍において前記測定データと前記基準データとを比較して各気筒毎の燃焼状態を推定する推定ステップと、該推定した各気筒の燃焼状態を基に該当する気筒の燃焼状態を制御する制御パラメータを補正する補正ステップと、該補正された制御パラメータにより各気筒の燃焼状態を制御する制御ステップとを有してなるエンジンの燃焼制御方法において、
前記基準作成ステップが、各気筒におけるクランク角の上死点前と上死点近傍の測定データを基に、前記基準データを算出する数式モデルを作成する数式モデル作成ステップを含むと共に、
前記推定ステップにおいて、前記エンジンの燃焼状態として燃料着火時期を推定するものであって、前記測定データと前記基準データとの差である偏差データが所定の閾値を超える時期を検出して回転変動開始時期とし、この回転変動開始時期から予め作成された燃料着火時期マップに基づいて燃料着火時期を推定すると共に、
前記補正ステップにおいて、前記推定された燃料着火時期と燃料噴射発令時期とを比較して該燃料噴射発令時期を補正し、
更に、前記推定ステップにおいて、燃料着火時期が上死点前であるか、上死点以後であるかを推定し、前記偏差データが所定の閾値を超える時期である回転変動開始時期を検出する際に、推定した前記燃料着火時期が上死点前で有る場合には前記偏差データから算出する変動値を前記所定の閾値として採用し、推定した前記燃料着火時期が上死点以後で有る場合には予め設定した固定値を前記所定の閾値として採用することを特徴とするエンジンの燃焼制御方法。
【請求項２】
エンジンのクランク軸の回転変動波形を回転センサで検出して測定データを得る回転変動検出手段と、回転変動波形データにおける該検出された測定データから基準データを算出する基準データ作成手段と、該エンジンの各気筒の上死点近傍において前記測定データと前記基準データとを比較して各気筒毎の燃焼状態を推定する燃焼状態推定手段と、該推定した各気筒の燃焼状態を基に該当する気筒の燃焼状態を制御する制御パラメータを補正する制御パラメータ補正手段と、該補正された制御パラメータにより各気筒の燃焼状態を制御する燃焼制御手段とを有してなるエンジンの燃焼制御装置であって、
前記基準データ作成手段が、各気筒におけるクランク角の上死点前と上死点近傍の測定データを基に、前記基準データを算出する数式モデルを作成する数式モデル作成手段を含むと共に、
前記燃焼状態推定手段は、前記エンジンの燃焼状態として燃料着火時期を推定するものであって、前記測定データと前記基準データとの差である偏差データが所定の閾値を超える時期を検出して回転変動開始時期とし、この回転変動開始時期から予め作成された燃料着火時期マップに基づいて燃料着火時期を推定すると共に、
前記制御パラメータ補正手段において、前記推定された燃料着火時期と燃料噴射発令時期とを比較して該燃料噴射発令時期を補正し、
更に、前記燃焼状態推定手段は、燃料着火時期が上死点前であるか、上死点以後であるかを推定し、前記偏差データが所定の閾値を超える時期である回転変動開始時期を検出する際に、推定した前記燃料着火時期が上死点前で有る場合には前記偏差データから算出する変動値を前記所定の閾値として採用し、推定した前記燃料着火時期が上死点以後で有る場合には予め設定した固定値を前記所定の閾値として採用することを特徴とするエンジンの燃焼制御装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、エンジン出力軸の回転変動から、各気筒内の燃料の着火時期等の燃焼状態を検出して、各気筒毎のインジェクタの燃料の噴射の制御パラメータを補正して、各気筒間の燃焼のバラツキを減少するエンジンの燃焼制御方法とその装置に関するものである。
【０００２】
【従来の技術】
自動車等の多気筒のディーゼルエンジンの燃料噴射においては、高圧燃料ポンプから供給された高圧の燃料をコモンレールに貯留し、このコモンレールから各気筒のインジェクタへ燃料を供給して燃料噴射するコモンレール式燃料噴射システムが知られている。
【０００３】
このコモンレールから各インジェクタの噴孔に至る燃料経路には高圧の燃料の通過及び遮断を制御するための電磁弁が備えられており、エンジンのコントローラにより、エンジンの運転状態に応じて、最適な条件で燃料噴射ができるように、コモンレールの圧力と各インジェクタ用の電磁弁の制御が行われている。
【０００４】
しかしながら、燃料噴射の圧力と噴射の時期（タイミング）や噴射時間を同じ条件に設定しても、インジェクタの燃料噴射特性にバラツキがあるために、また、経年変化により燃料噴射特性が変化するために、インジェクタ毎に燃料噴射量が異なってしまい、各気筒間における燃焼状態にバラツキを生じるので、エンジンの出力トルクが変動し、振動や騒音が発生するという問題がある。このエンジンの振動と騒音は、特にエンジンのアイドル運転時や低負荷運転時に大きくなり、この問題の解決が重要課題となっている。
【０００５】
この問題の対処法の一つとして、各インジェクタを同一の目標燃料噴射量に基づいて噴射量を制御する場合に、燃料の実噴射量の多少に応じて各気筒におけるクランク軸の回転速度が増減することに着目して、この回転速度の変動から各気筒の燃焼状態を検出して燃焼のバラツキをなくす制御を行う燃焼制御方法が提案されている。
【０００６】
つまり、ピストンが上死点に至るまで吸気弁及び排気弁が共に閉じた状態にあるとすると、ピストンが上死点に達する前の圧縮行程中に燃焼が開始する場合には、燃焼室内の燃料の着火・燃焼による発生熱によってピストンの上昇を抑える方向のトルクが発生するので、クランク軸の回転速度は減少する。反対に、ピストンが上死点に達した後の膨張行程中に燃焼が開始する場合には、燃料が着火・燃焼して燃焼室内で発生した熱を受けた燃焼ガスによってピストンの下降を後押しするトルクが発生するので、クランク軸の回転速度は増加する。
【０００７】
従って、燃焼室内で燃料が燃焼すると、エンジンの出力軸にトルクが発生し、このトルクがクランク軸の回転速度を増減するように作用するので、回転速度はクランク軸の回転角度に応じて時々刻々変動し、平均回転速度の上下に変動を繰り返している。
【０００８】
また、モータリング等の燃料噴射しない運転状態におけるエンジン出力軸の回転速度変動と、実際に燃料を燃焼室に噴射した燃焼時の回転速度変動との偏差は、燃焼室内に発生する熱に起因すると考えられるので、この偏差で判断すると、燃焼時の回転速度変動の性状を捕らえることが容易となり、この偏差に着目して、燃料噴射時期及び燃料噴射量の補正を行うことができる。
【０００９】
この知見を得て、本発明の出願人は、特願平１１−９６９６０号公報（特開２０００−２９１４６６号公報）において、熱発生率に相当するエンジン出力軸の回転速度を検出することにより、各気筒の燃料噴射量のバラツキによる燃焼のバラツキを把握し、また、燃料の着火時期のバラツキに起因する燃焼バラツキも把握して、燃焼のタイミングも考慮して燃焼バラツキをより一層正確に捉えて解消するエンジンの燃料噴射制御装置を提案している。
【００１０】
更に、本発明の発明者の１人は、特願平１１−２５７０８９号において、回転センサで検出した、エンジンのクランク軸の上死点前の所定の回転角度の範囲における燃料着火前の実回転変動データを、上死点を境に反転させて上死点後基準回転変動波形データとし、この基準のデータと上死点後の計測データとの偏差に基づいてエンジンの燃焼状態を算出するエンジンの燃焼状態検出方法を提案している。
【００１１】
【発明が解決しようとする課題】
しかしながら、上記のエンジンの燃焼状態検出方法を使用したエンジンの燃焼制御装置においては、実際に、車載した状態ではモータリング運転状態におけるクランク軸の回転速度のデータを得ることは困難であるため、予め、モータリング状態におけるクランク軸の回転速度に関する基準データを用意してコントローラに記憶させておく必要がある。
【００１２】
つまり、この燃焼状態検出方法では、ピストンが上死点通過後（以下「上死点後」という）に着火する場合には、ピストンが上死点通過前（以下「上死点前」という）の測定データを反転して上死点後の基準データにするので、この部分の基準データに関しては予め記憶しておく必要がないが、上死点前に着火する場合には、測定データがモータリング運転状態における基準データと異なってくるので、この反転の方法で基準データを作ることができなくなる。
【００１３】
そのため、上死点前に燃料が着火する場合には、実際の測定データと同じ運転状態の上死点前基準回転変動波形データを計測する等して、予め用意しなければならないという問題が生じる。
【００１４】
しかも、このクランク軸の回転変動の波形は、エンジンの回転速度や負荷状態等に応じて変化するため、この基準データの記憶には多量のメモリが必要になるという問題が生じる。
【００１５】
更に、上死点前後において燃料の着火が行われる場合には、燃料の着火が上死点前と上死点後とのどちらかで行われているかによって、燃焼状態検出方法が異なってくるが、上記の燃焼状態判別方法では、測定データと予め用意した基準データとに基づいて、燃料の着火が上死点前と上死点後とのどちらかで行われているかを自動判別することまで考慮されて折らず、そのため、上死点前に着火が発生する場合には、上記の燃焼状態判別方法を適用できないという問題がある。
【００１６】
本発明は、上述の知見を得て、従来技術の問題を解決するためになされたものであり、その目的は、クランク軸の回転変動波形を回転センサで検出し、この検出した測定データから数式モデルを作成し、この数式モデルから基準データを作成することによって、この測定データと基準データとから、各気筒における燃料の着火時期がピストンの上死点通過の前後のいずれかにあるかを判別して、燃料の着火時期を推定し、更に、予め作成された燃料着火時期マップに基づいて算出される基準噴射発令時期との差を求めて、各気筒における噴射発令時期を補正して、各気筒内の燃焼のバラツキを減少させるエンジンの燃焼制御方法とその装置を提供することにある。
【００１７】
【課題を解決するための手段】
以上のような目的を達成するためのエンジンの燃焼制御方法は、以下のように構成される。
【００１８】
１）エンジンのクランク軸の回転変動波形を回転センサで検出して測定データを得る回転変動検出ステップと、回転変動波形データにおける該検出された測定データから基準データを算出する基準作成ステップと、該エンジンの各気筒の上死点近傍において前記測定データと前記基準データとを比較して各気筒毎の燃焼状態を推定する推定ステップと、該推定した各気筒の燃焼状態を基に該当する気筒の燃焼状態を制御する制御パラメータを補正する補正ステップと、該補正された制御パラメータにより各気筒の燃焼状態を制御する制御ステップとを有してなるエンジンの燃焼制御方法において、前記基準作成ステップが、各気筒におけるクランク角の上死点前と上死点近傍の測定データを基に、前記基準データを算出する数式モデルを作成する数式モデル作成ステップを含むと共に、前記推定ステップにおいて、前記エンジンの燃焼状態として燃料着火時期を推定するものであって、前記測定データと前記基準データとの差である偏差データが所定の閾値を超える時期を検出して回転変動開始時期とし、この回転変動開始時期から予め作成された燃料着火時期マップに基づいて燃料着火時期を推定すると共に、前記補正ステップにおいて、前記推定された燃料着火時期と燃料噴射発令時期とを比較して該燃料噴射発令時期を補正し、更に、前記推定ステップにおいて、燃料着火時期が上死点前であるか、上死点以後であるかを推定し、前記偏差データが所定の閾値を超える時期である回転変動開始時期を検出する際に、推定した前記燃料着火時期が上死点前で有る場合には前記偏差データから算出する変動値を前記所定の閾値として採用し、推定した前記燃料着火時期が上死点以後で有る場合には予め設定した固定値を前記所定の閾値として採用するように構成される。
【００１９】
この「上死点近傍」とはピストンが上死点を通過する前後の時期をいい、以下で使用する「上死点前」とはピストンが上死点を通過する前の時期をいい、「上死点後」とはピストンが上死点を通過する後の時期をいう。
【００２０】
この構成によれば、制御しようとしている運転状態で、燃焼状態の検出に使用する基準データを、数式モデル（数学理論モデル）を用いることにより、その運転状態で計測されるエンジンのクランク軸の回転変動の測定データから作成できる。この基準データは、その運転状態で、仮に燃料を着火しないあるいは燃料を噴射しない場合（モータリング運転状態）に発生すると考えられる回転変動波形データに相当するデータである。
【００２１】
数式モデルを介在させることにより、各気筒における燃料の着火が上死点前後のどの間に発生するかには拘らず、この基準データを推測できるので、モータリング運転状態における回転変動波形の基準データを予め用意しなくても良いので、この基準データの収集の手間が省けると共に、この基準データのための記憶メモリも不要になる。
この構成によれば、測定データとこの測定データから数式モデルを介して得られた基準データとの比較により、燃料着火時期を精度良く推定でき、また、この推定した燃料着火時期から各気筒の燃料噴射発令時期を補正できるので、各気筒における燃焼状態を、エンジンの運転状態から得られる目標とする最適な燃焼状態に近づけることができる。
また、この構成によれば、燃料着火時期が上死点前であるか、上死点以後であるかを推定し、回転変動開始時期を検出する際の所定の閾値をそれぞれの状態に適した値にすることができるので、この回転変動開始時期から得られる燃料着火時期を精度良く推定できるので、各気筒における燃焼状態のバラツキをより効率的に減少できる。
【００２２】
２）上記のエンジンの燃焼制御方法の前記数式モデル作成ステップにおいて、前記数式モデルを作成する際に使用するデータを、通常の運転において着火が発生しない上死点前の所定のクランク角度範囲内における前記測定データからサンプリングしたデータと、上死点近傍における代表値とで構成する。
【００２３】
この通常の運転とは、故障などが生じない正常の運転状態における運転のことをいい、この着火が発生しない上死点前の所定のクランク角度範囲内は、正常な運転状態における燃料噴射発令時期の設定時期と各インジェクタ特性のバラツキ範囲とから予め予測できるので、この予測に基づいて設定される。
【００２４】
この構成によれば、燃料の着火が予想される上死点近傍以外では、測定データは、モータリング運転状態の回転変動波形データと略同じであるので、この測定データからのサンプリング値を採用し、燃料の着火が予想される上死点近傍では、測定データは、モータリング運転状態の回転変動波形データから大きく外れる恐れがあるので、測定データそのものを使用せずに、測定データから計算される代表値を使用するので、精度良く、基準データを求めることができる。
【００２５】
３）上記のエンジンの燃焼制御方法の前記数式モデル作成ステップにおいて、前記代表値を上死点近傍における前記測定データの最低値と上死点における測定データとの平均値とする。
【００２６】
この構成によれば、代表値を簡単に用意することができ、数式モデルの決定が簡便にできる。
【００２７】
４）上記のエンジンの燃焼制御方法の前記数式モデル作成ステップにおいて、前記数式モデルが、θをクランク角度、Ｙを回転速度、ｍ＝１，２，３，・・・Ｍとして、Ｙ＝Ｆ（θ）＝Ａ＋Ｂ×ｅｘｐ（Σm ａm ×ｃｏｓ（ｍθ））で表される。
【００２８】
この構成によれば、エンジンの出力軸であるクランク軸の回転速度の変動波形データを少ないパラメータで効率よく近似できる。そのため、演算時間やメモリーの節約をできる。
【００３３】
そして、上記のエンジンの燃焼制御方法を実現するためのエンジンの燃焼制御装置は、以下のように構成される。
【００３４】
１）エンジンのクランク軸の回転変動波形を回転センサで検出して測定データを得る回転変動検出手段と、回転変動波形データにおける該検出された測定データから基準データを算出する基準データ作成手段と、該エンジンの各気筒の上死点近傍において前記測定データと前記基準データとを比較して各気筒毎の燃焼状態を推定する燃焼状態推定手段と、該推定した各気筒の燃焼状態を基に該当する気筒の燃焼状態を制御する制御パラメータを補正する制御パラメータ補正手段と、該補正された制御パラメータにより各気筒の燃焼状態を制御する燃焼制御手段とを有してなるエンジンの燃焼制御装置であって、前記基準データ作成手段が、各気筒におけるクランク角の上死点前と上死点近傍の測定データを基に、前記基準データを算出する数式モデルを作成する数式モデル作成手段を含むと共に、前記燃焼状態推定手段は、前記エンジンの燃焼状態として燃料着火時期を推定するものであって、前記測定データと前記基準データとの差である偏差データが所定の閾値を超える時期を検出して回転変動開始時期とし、この回転変動開始時期から予め作成された燃料着火時期マップに基づいて燃料着火時期を推定すると共に、前記制御パラメータ補正手段において、前記推定された燃料着火時期と燃料噴射発令時期とを比較して該燃料噴射発令時期を補正し、更に、前記燃焼状態推定手段は、燃料着火時期が上死点前であるか、上死点以後であるかを推定し、前記偏差データが所定の閾値を超える時期である回転変動開始時期を検出する際に、推定した前記燃料着火時期が上死点前で有る場合には前記偏差データから算出する変動値を前記所定の閾値として採用し、推定した前記燃料着火時期が上死点以後で有る場合には予め設定した固定値を前記所定の閾値として採用するように構成される。
【００３５】
２）上記のエンジンの燃焼制御装置において、前記数式モデル作成手段が、前記数式モデルを作成する際に使用するデータを、通常の運転において着火が発生しない上死点前の所定のクランク角度範囲内における前記測定データからサンプリングしたデータと、上死点近傍における代表値とで構成する。
【００３６】
３）上記のエンジンの燃焼制御装置の前記数式モデル作成手段において、前記代表値を上死点近傍における前記測定データの最低値と上死点における測定データとの平均値とする。
【００３７】
４）上記のエンジンの燃焼制御装置の前記数式モデル作成手段において、前記数式モデルが、θをクランク角度、Ｙを回転速度、ｍ＝１，２，３，・・・Ｍとして、Ｙ＝Ｆ（θ）＝Ａ＋Ｂ×ｅｘｐ（Σm ａm ×ｃｏｓ（ｍθ））で表されるように構成する。
【００４０】
なお、上記の数式モデルとしては、Ｙ＝Ｆ（θ）とは別の、Ｙ＝Ｆ’（θ）＝Ａ＋Σm ａm ×ｃｏｓ（ｍθ）のフーリエ関数の数式モデル等でも実施可能である。但し、Ｙ＝Ｆ（θ）の場合には、ｍ＝１で十分に正確な基準データを作成できるが、Ｙ＝Ｆ’（θ）の場合にはｍ≧２が必要となる。
【００４１】
【発明の実施の形態】
以下、図面を用いて、本発明に係る実施の形態のエンジンの燃焼制御方法及びその装置について、推定する燃焼状態を燃料着火時期とし、補正される制御パラメータと燃料噴射発令時期として説明する。
【００４２】
〔装置の構成〕
このエンジンの燃焼制御装置１は、図１に示すように、回転変動検出手段１０と、基準データ作成手段２０と、燃焼状態推定手段３０と、制御パラメータ補正手段４０と、燃焼制御手段５０とを有して構成される。
【００４３】
この回転変動検出手段１０は、エンジンのクランク軸の回転状態を回転センサで、クランク軸が一定角度回転する毎に発生するクランク角センサの出力パルスの時間間隔をカウンタによって求め、時々刻々のクランク軸の回転速度を検出し、この回転速度の変動を示す回転変動波形の時系列の測定データを得るものである。
【００４４】
この回転センサは、エンジンの出力軸であるクランク軸の回転変動波形データを検出するが、高分解能を有する光学式や磁気式のロータリエンコーダを使用することができる。この光学式のロータリエンコーダは、クランク軸の回転と共に回転する回転円板に多数のスリットを設け、このスリットを通過する発光素子からの光を受光素子で受けて回転を検出する。
【００４５】
また、基準データ作成手段２０は、検出された測定データから基準データを算出するものであり、エンジンの各気筒におけるクランク角の上死点前と上死点近傍の測定データを基にサンプリング値や代表値を選出する手段２１や、基準データを算出する数式モデルを作成する数式モデル作成手段２２や、数式モデルから基準データを算出する手段２３等を含んで構成される。
【００４６】
そして、燃焼状態推定手段３０は、エンジンの各気筒の上死点近傍において測定データと基準データとを比較して各気筒毎の燃焼状態を推定するものであり、燃料着火時期が上死点前後のいずれかにあるかを判断する手段３１や、回転変動開始時期を検出する手段３２や、燃料着火時期を推定する手段３３等を含んで構成される。
【００４７】
また、制御パラメータ補正手段４０は、推定した各気筒の燃焼状態を基に該当する気筒の燃焼状態を制御する制御パラメータを補正する手段であり、ここでは推定された燃料着火時期と、制御の指令を出している燃料噴射発令時期との比較に基づいて補正量を計算する手段４１や、制御パラメータを補正する手段４２等を含んでおり、個々では、制御パラメータの一つである燃料噴射発令時期を補正している。
【００４８】
そして、燃焼制御手段５０は、燃焼状態を制御するための燃料噴射弁制御手段５１やコモンレール圧を調整する手段や吸気量をする手段やＥＧＲ量を調整する手段等の種々の制御手段を有して構成される。
【００４９】
これらの各手段は、エンジンを制御するエンジンコントロールユニット（ＥＣＵ）と呼ばれる制御装置の制御プログラムの一部として構成される。また、ここで使用される回転センサからの測定データもエンジンコントロールで使用している回転センサの測定データが利用される。
【００５０】
〔エンジンの燃焼制御方法〕
上記の構成のエンジンの燃料制御装置を使用して、エンジンの燃焼制御を行う。以下に、このエンジンの燃焼制御の概略について説明するが、この燃焼制御は、図２に示すような、フローに従って行われる。
【００５１】
先ず、ステップＳ１で、燃料の噴射制御を行う運転状態において、クランク軸の回転変動波形データの測定データＲｍｉを得て、次のステップＳ２で、この測定データＲｍｉから数式モデルを介在させて、基準データＲｓｉを作成する。
【００５２】
そして、ステップＳ３で、この基準データＲｓｉと測定データＲｍｉとから、燃料着火時期を推定し、ステップＳ４で制御パラメータを補正し、ステップＳ５でこの補正した制御パラメータにより燃焼を制御する。
【００５３】
以下、各ステップに関して説明する。
【００５４】
〔回転変動検出ステップＳ１〕
回転変動検出ステップＳ１では、燃料の噴射制御を行う運転状態において、エンジンの出力軸であるクランク軸の回転を回転センサで検出し、この検出データからクランク軸の回転変動波形データの測定データＲｍｎを得る。
【００５５】
〔基準作成ステップＳ２〕
次に、基準作成ステップＳ２では、図３に示すサブフローに従って、燃料の噴射制御を行う運転状態において、クランク角度でＴＤＣ（上死点）からマイナス約９０度くらい前の範囲（Ｚａ）内のエンジンの出力軸の回転変動波形データの測定データＲｍｉとピストンが上死点（ＴＤＣ）を通過する時期Ｔtdc （上死点位置情報）とから、ピストンが上死点を通過する前後のクランク角度表示で約±９０度くらいの範囲内の基準データＲｓｎを、数学理論モデルに基づいて求める。
【００５６】
なお、この基準データＲｓｎは、同じ運転状態において、モータリング運転のような、燃料を着火しないあるいは燃料を噴射しない状態、即ち、燃料の燃焼による発熱が無い状態の回転変動波形データに対応する比較用のデータである。
【００５７】
そして、この基準作成ステップＳ２が実行されると、図３に示すように、ステップＳ２１で、測定データＲｍｉから基準データＲｓｎを算出するために使用するサンプリング値ＲｄｋをＫ個（図５では５点）サンプリングする。
【００５８】
このサンプリング値Ｒｄｋは、クランク角度でＴＤＣ（上死点）からマイナス約９０度くらい前の範囲（Ｚａ）内で、燃料の着火が発生する可能性のあるＴＤＣ近傍（Ｚｂ）を除いて、Ｋ点採用される。このサンプリングはクランク角度θに関して等間隔の方がデータ処理が楽になるので好ましいが必ずしも等間隔でなくてもよい。
【００５９】
次に、ステップＳ２２で、燃料の着火が発生する可能性のあるＴＤＣ近傍（Ｚｂ）における代表値Ｒｄtdc を算定する。この算定は、図５に示すように、測定データＲｍｉから最小値Ｒｍmin とＴＤＣにおける実測値Ｒｍtdc を検出して、両者の平均からＴＤＣにおける代表値Ｒｄtdc （＝（Ｒｍmin ＋Ｒｍtdc ）／２）を算出する。ちなみに、燃料の着火がＴＤＣの後で行われる場合には、図６や図７に示すように、代表値Ｒｄtdc （≒Ｒｍmin ≒Ｒｍtdc ≒０）は略ゼロとなる。
【００６０】
そして、次に、ステップＳ２３でサンプリング値Ｒｄｋと代表値Ｒｄtdc から、基準データＲｓｎの関数Ｆ（θｎ）、即ち、サンプリング値Ｒｄｋと代表値Ｒｄtdc を近似表現する非線型数学モデルで作った回帰曲線或いは代数曲線である（１）式のパラメータＡ，Ｂ，ａm を求める。
【００６１】
つまり、数式モデルＦ（θ）を、θを上死点をゼロ度とするクランク角度、Ｙを回転速度として、（１）式で表す。そして、この関数のパラメータＡ，Ｂ，ａm を未知数として、サンプリング値Ｒｄｋ、Ｔｄｋと、代表値データＲｄtdc 、Ｔtdc を代入して、（２）式及び（３）式のＫ＋１個の連立方程式を作成し、これを最小二乗法等により解いてパラメータＡ，Ｂ，ａm の各数値を求める。
【００６２】
Ｙ＝Ｆ（θ）＝Ａ＋Ｂ×ｅｘｐ（Σm ａm ×ｃｏｓ（ｍθ））
ｍ＝１，２，３，・・・Ｍ                              （１）
Ｒｄｋ＝Ａ＋Ｂ×ｅｘｐ（Σm ａm ×ｃｏｓ（ｍ×Ｔｄｋ））
ｋ＝１，２，３，・・・Ｋ        （Ｋ＋１≧Ｍ＋２）    （２）
Ｒｄtdc ＝Ａ＋Ｂ×ｅｘｐ（Σm ａm ×ｃｏｓ（ｍ×Ｔtdc ））  （３）
ステップＳ２４では、このパラメータ、Ａ，Ｂ，ａm の数値を使って、（４）式により、上死点以前の基準データＲｓｉを算出する。
【００６３】
Ｒｓｉ＝Ａ＋Ｂ×ｅｘｐ（Σm ａm ×ｃｏｓ（ｍ×Ｔｓｉ））
ｉ＝１，２，３，・・・Ｉ                              （４）
そして、ＴＤＣより後の範囲（Ｚｃ）の基準データＲｓｊも算出する。この算出は（１）式の代数式より算出することができるが、ＴＤＣより前（θのマイナス側）の範囲（Ｚａ）の基準データＲｓｉをＴＤＣ軸に対して反転してＴＤＣより後（θのプラス側）の範囲（Ｚｃ）の基準データＲｓｊを算出することもできる。
【００６４】
このＴＤＣより前の範囲（Ｚａ）の基準データＲｓｉとＴＤＣより後の範囲（Ｚｃ）の基準データＲｓｊとを合わせて、基準データＲｓｎとする。
【００６５】
この図３のフローによって算定された基準データＲｓｎと測定データＲｍｎ等を図５〜図７に示す。
【００６６】
〔推定ステップＳ３〕
次に、図２の推定ステップＳ３について説明するが、このステップＳ３は，燃料の着火時期がＴtdc の前後にあるかの判定するステップＳ３Ａと、上死点前に燃料が着火する場合の着火時期を推定するステップＳ３Ｂと、上死点後に燃料が着火する場合の着火時期を推定するステップＳ３Ｃとからなる。
〔燃料の着火時期判定ステップＳ３Ａ〕
この燃料の着火時期ＴｆがＴtdc （ピストンがＴＤＣを通過する時期）の前後のどちらにあるかを判定するステップＳ３Ａは、図４に示すようなフロー及びアルゴリズムに従って行われ、測定データＲｍｎとその偏差データΔＲｄｎに基づいて、上死点Ｔtdc 前後のどの時期に実際の着火時期があるかを自動的に判別する。
【００６７】
先ず、最初にステップＳ３１において、測定データＲｍｎの最小値Ｒｍmin とその発生時期Ｔmin を求める。
【００６８】
そして、ステップＳ３２において、発生時期Ｔmin ＋Δｔ’がＴtdc より小さいか（前か）を判定し、小さければ、Ｔtdc より前に着火していると判定し、ステップＳ５１に行く。このΔｔ’は噴射の発令から着火までの時間遅れに関する調整項であり、例えばクランク角で３度である。
【００６９】
また、ステップＳ３２において、発生時期Ｔmin ＋Δｔ’がＴtdc より大きければ、更に詳細に調べるために、ステップＳ３３に行く。
【００７０】
このステップＳ３３では、図５〜図７に示すような測定データＲｍｎと基準データＲｓｎから、図８〜図１２に示すような偏差データΔＲｄｎ（＝Ｒｍｎ−Ｒｓｎ）を作成する。この偏差データΔＲｄｎの図８、図１０、図１２は、それぞれ図５、図６、図７に対応する。
【００７１】
次に、ステップＳ３４では、図８〜図１２に示すような偏差データΔＲｄｎから、現在の噴射発令時期Ｔｏより所定の時間ΔＴ前の範囲（Ｚｄ）の間における最大偏差値ΔＲｄＢと、現在の噴射発令時期Ｔｏより後（Ｚｅ）の最小値ΔＲｄＣを求める。
【００７２】
そして、ステップＳ３５で、最小値ΔＲｄＣの絶対値｜ΔＲｄＣ｜と最大偏差値ΔＲｄＢとを比較し、絶対値｜ΔＲｄＣ｜が大きければ、最小値ΔＲｄＣの発生時期Ｔｃを主体に判定するとしてステップＳ３６に行き、絶対値｜ΔＲｄＣ｜が小さければ、最大偏差値ΔＲｄＢを主体に判定するとしてステップＳ４１に行く。
【００７３】
ステップＳ３６では、更に、最小値ΔＲｄＣの発生時期Ｔｃの直前において、偏差データΔＲｄｎがΔＲｄC2＝ΔＲｄＣ／２の値となる第１参照時期Ｔc2を求める。
【００７４】
そして、ステップＳ３７で、最小値ΔＲｄＣの発生時期Ｔｃと第１参照時期Ｔc2の両方が、噴射発令時期ＴｏとＴtdc の間にあるか否かを判定する。間にあれば、Ｔtdc より前に着火していると判定し、ステップＳ５１に行く。間に無ければ、ステップＳ４２に行く。
【００７５】
ステップＳ４１では、最大偏差値ΔＲｄＢの発生位置Ｔｂが噴射発令時期ＴｏとＴtdc の間で且つ、後ろから第１番目の第１極小値ΔＲｄＤがあるか否か（第１極小値ΔＲｄＤがシステムの分解能よりも大きくて検出されるか否か）を判定する。無ければ、Ｔtdc より後に着火していると判定し、ステップＳ５２に行く。あれば、ステップＳ４２に行く。
【００７６】
ステップＳ４２では、第１極小値ΔＲｄＤが発生する時期Ｔｅの直前の偏差データΔＲｄｎがΔＲｄD2＝ΔＲｄＤ／２の値となる第２参照時期ＴD2を求める。そして、ステップＳ４３で第２参照時期ＴD2が最小噴射発令時期ＴｏとＴtdc の間にあるか否かを判定する。有れば、Ｔtdc より前に着火していると判定し、ステップＳ５１に行く。無ければ、Ｔtdc より後に着火していると判定し、ステップＳ５２に行く。
【００７７】
以上のフローに従うと、Ｔtdc より前に着火していると判定する場合は、次の３つの場合であり、それ以外は、Ｔtdc より後で着火していると判定する。
【００７８】
先ず第１の場合は、最小値Ｒｍmin の発生時期Ｔmin がＴtdc −Δｔ’より小さい図５に示すような場合である。
【００７９】
また、第２の場合は、現在の噴射発令時期Ｔｏより後（Ｚｅ）の最小値ΔＲｄＣの絶対値が現在の噴射発令時期Ｔｏより所定の時間ΔＴ前の範囲（Ｚｄ）の間における最大偏差値ΔＲｄＢよりも大きく、且つ、最小値ΔＲｄＣの発生時期Ｔｃの直前において、偏差データΔＲｄｎがΔＲｄC2＝ΔＲｄＣ／２の値となる第１参照時期ＴC2と、最小値ΔＲｄＣの発生時期Ｔｃとの両方が、噴射発令時期ＴｏとＴtdc の間にある場合である。
【００８０】
また、第３の場合は、噴射発令時期ＴｏとＴtdc の間で且つ、後ろから第１番目の第１極小値ΔＲｄＤに対して、第１極小値ΔＲｄＤが発生する直前の偏差データΔＲｄｎがΔＲｄD2＝ΔＲｄＤ／２の値となる第２参照時期ＴD2が最小噴射発令時期ＴｏとＴtdc の間にある場合である。
【００８１】
〔燃料の着火時期の推定〕
次に、上死点前後のどの間に燃料の着火があるかが推定できたので、それぞれに適したアルゴリズムによって、偏差データΔＲｄｎが所定の版定置ΔＲＣを超える回転変化開始時期Ｔｖを検出し、回転変化開始時期Ｔｖと燃料着火時期Ｔｆとを対応させた図１３に示すような燃料着火時期マップ（ＴｖとＴｆの関係）に基づいて燃料着火時期Ｔｆを求める。
【００８２】
この図１３の燃料着火時期マップは、予め作成してコントローラに記憶させておくものであり、燃料噴射時期Ｔｏを変化させた場合の気筒の筒内圧力から計算した熱発生率から求めた燃料着火時期Ｔｆ（横軸）と、図８〜図１２に示すような偏差データΔＲｄｎから求められる回転変化開始時期Ｔｖ（縦軸）との関係を示すものである。
【００８３】
この図１３から分かるように、燃料着火時期Ｔｆと回転変化開始時期Ｔｖとは、実質的に一次式の相関関係にあるので、この相関関係を考慮してこれらのデータから一次回帰で得られた回帰直線（Ｔｆ＝ａｉ×Ｔｖ＋ｂｉ）のパラメータａ１，ｂ１，ａ２，ｂ２をマップデータとして記憶しておき、このマップデータに基づいて、偏差データΔＲｄｎから求められる回転変化開始時期のデータＴｖＳをこの一次回帰直線にあてはめて、燃料着火時期ＴｆＳ（＝ａｉ×ＴｖＳ＋ｂｉ）を算出する。この燃料着火時期ＴｆＳが実際に発生していると推定された燃料着火時期Ｔｆである。
【００８４】
〔上死点前に燃料の着火する場合の着火時期の推定〕
上死点前に燃料着火の場合には、図８に示すように、偏差データΔＲｄｎに雑音が多いので、先ず、噴射発令時期Ｔｏの後とＴtdc の前の範囲において、最小値ΔＲｄＣを求め、その半分ΔＲｄＤ＝ΔＲｄＣ／２を第１閾値ΔＲＣ１とする。
【００８５】
そして、噴射発令時期Ｔｏの後における、第１閾値ΔＲＣ１より小さい偏差データΔＲｄｎが発生する時期を回転変化開始時期ＴｖＳとし、この回転変化開始時期ＴｖＳの値から、燃料着火時期マップに基づいて燃料着火時期ＴｆＳを求める。
【００８６】
〔上死点以後に燃料に着火する場合の着火時期の推定〕
上死点以後に燃料の着火する場合の燃料着火時期Ｔｆの推定は、図６、図７、図１０〜図１２に示すように、特定の正の値を所定の判定値ΔＲＣ２とし、これを偏差データΔＲｄｎに対する第２閾値ΔＲＣ２にして、この第２閾値ΔＲＣ２を超える時期を回転変化開始時期ＴｖＳとして、この回転変化開始時期ＴｖＳの値から、燃料着火時期マップ（Ｔｖ−Ｔｆ）に基づいて燃料着火時期ＴｆＳを求め、燃料着火時期Ｔｆとする。
【００８７】
この第２閾値ΔＲＣ２は固定値であり、上死点前に燃料着火がある場合の第１閾値ΔＲＣ１にように変動しない。
【００８８】
〔補正ステップＳ４〕
次に、図２の補正ステップＳ４について説明する。
【００８９】
このステップＳ４は，各気筒毎に推定された実際の燃料着火時期Ｔｆと各、エンジンの運転状態から目標として設定される燃料着火時期Ｔｆｏを比較し、この両者の差を、各気筒のインジェクタ毎に設定される燃料噴射発令時期Ｔｏに反映させて、実際の燃料着火時期Ｔｆを目標として設定される燃料着火時期Ｔｆｏに近付けるためのものである。
【００９０】
そして、この推定された実際の燃料着火時期Ｔｆと目標として設定される燃料着火時期Ｔｆｏとの差ΔＴｆ＝Ｔｆ−Ｔｆｏを使用して、新たな燃料噴射発令時期ＴｏＮをＴｏＮ＝Ｔｏ＋ΔＴｆとして補正する。
【００９１】
〔制御ステップＳ５〕
そして、図２の制御ステップＳ５では、新たな燃料噴射発令時期ＴｏＮを噴射発令時期Ｔｏとして、各気筒のインジェクタの制御弁等を制御し、各気筒の燃焼状態を制御する。
【００９２】
〔効果〕
以上のフローに従ったエンジンの燃焼制御方法によれば、燃焼制御すべき運転状態において、エンジンのクランク軸の実際の回転変動を回転センサによって検出し、燃料の着火時期Ｔｆが上死点Ｔtdc の前後のいずれかの時期にあるかを問わずに、回転変動波形の測定データＲｍｎと燃料の噴射発令時期（燃料の噴射命令を出すタイミング）Ｔｏとピストンが上死点を通過する上死点通過時期Ｔtdc の情報とから、測定データＲｍｎを基に作成される数学理論モデル（Ｆ（θ））に基づいて基準データＲｓｎを作成し、この基準データＲｓｎと測定データＲｍｎとの比較から、燃料着火時期Ｔｆが上死点Ｔtdc の前後のいずれかにあるかを判別できる。
【００９３】
そして、燃料着火時期Ｔｆが上死点Ｔtdc 以前にある場合には、偏差データΔＲｄｎから求まる第１閾値ΔＲＣ１を使用し、また、燃料着火時期Ｔｆが上死点Ｔtdc の後にある場合には、偏差データΔＲｄｎとは関係なく事前に設定した固定値である第２閾値ΔＲＣ２を使用して、偏差データΔＲｄｎがこの所定の判定値である第１閾値ΔＲＣ１又は第２閾値ΔＲＣ２を超える回転変化開始時期ＴｖＳを検出して、この検出した時期ＴｖＳから予め作成された燃料着火時期マップ（Ｔｖ−Ｔｆ）に基づいて、燃料着火時期Ｔｆを推定することができる。
【００９４】
この燃料着火時期Ｔｆと、エンジンの回転数や負荷等から決まるエンジンの運転状態における目標とする燃料着火時期Ｔｆｏとの差ΔＴｆをチェックして、現状の噴射発令時期Ｔｏを補正することができる。
【００９５】
この補正した噴射発令時期ＴｏＮで各インジェクタの燃料噴射を制御することにより、各気筒における燃料着火時期Ｔｆを，目標の燃料着火時期Ｔｆｏにフィードバック制御することができるので、各気筒における燃料着火時期Ｔｆのバラツキを防止できる。
【００９６】
【発明の効果】
以上説明したように、本発明に係るエンジンの燃焼制御方法とその装置によれば、次のような効果を奏することができる。
【００９７】
エンジンのクランク軸の回転センサの出力から回転変動波形を検出し、この回転変動波形の測定データから、仮に燃料を着火しない又は燃料を噴射しないモーターリング運転を行った場合に出現する回転変動波形に相当すると考えられる基準データを算出する際に、上死点以前の測定データからサンプリングしたデータを基に数式で表現した数式モデルを作成し、この数式モデルから基準データを算出しているので、燃料着火時期が上死点前である場合も、また、上死点後である場合も、基準データを容易に作成することができる。
【００９８】
そのため、基準データをモータリング運転から作成してそのデータを予めコントローラに入力しておく必要がなくなり、測定データとこの測定データから得られる基準データとから燃焼状態をきめ細かく推定できる。
【００９９】
また、数式モデル作成において、数式モデルを作成する際に使用するデータを、通常の運転において着火が発生しない上死点前の所定のクランク角度範囲内からサンプリングしたデータと、上死点近傍の測定データから算出した代表値とで構成し、燃料の着火が発生する可能性がある上死点近傍では、測定データそのものの使用を避けて上死点近傍における代表値を使用するので、燃料の着火の影響を受けない基準データを得ることができる。そのため、従来技術では適用できなかった燃料着火が上死点前の場合でも燃焼状態を推定できる。
【０１００】
そして、エンジンの燃焼状態の一つである燃料着火時期を推定する場合に、測定データと基準データとの差である偏差データが所定の閾値を超える回転変化開始時期を検出し、この回転変化開始時期から予め作成された燃料着火時期マップに基づいて燃料着火時期を推定するので、簡便に燃料着火時期を推定できる。
【０１０１】
エンジンの燃料着火時期を推定する場合に、燃料着火時期が上死点以前であるか、上死点以後であるかを判定して、偏差データから回転変動開始時期を検出する際に使用する所定の閾値を、上死点以前では偏差データから算出する変動値とし、上死点以後では予め設定した固定値とすることにより、より精度良く簡便に燃料着火時期を推定することができる。
【０１０２】
また、この精度良く推定した燃料着火時期から各気筒の燃料噴射発令時期を補正して、燃焼状態を制御できるので、各気筒における燃焼状態を、エンジンの運転状態から得られる目標とする最適な燃焼状態に近づけることができ、各気筒における燃焼状態のバラツキをより効率的に減少できる。
【図面の簡単な説明】
【図１】本発明の実施の形態のエンジンの燃焼制御装置のシステム構成図である。
【図２】本発明の実施の形態のエンジンの燃焼制御方法の基本的なフローを示す図である。
【図３】図２の基準作成ステップの詳細なフローを示す図である。
【図４】図２の推定ステップの燃料の着火時期が上死点Ｔtdc の前後のどちらにあるかを判定するステップの詳細なフローを示す図である。
【図５】クランク軸の回転変動波形の例を示す図で、燃料着火が上死点前１０度の場合の例である。
【図６】クランク軸の回転変動波形の例を示す図で、燃料着火が上死点の場合の例である。
【図７】クランク軸の回転変動波形の例を示す図で、燃料着火が上死点後１０度の場合の例である。
【図８】クランク軸の回転変動波形の偏差データの例を示す図で、燃料着火が上死点前１０度の場合の例である。
【図９】クランク軸の回転変動波形の偏差データの例を示す図で、燃料着火が上死点前５度の場合の例である。
【図１０】クランク軸の回転変動波形の偏差データの例を示す図で、燃料着火が上死点の場合の例である。
【図１１】クランク軸の回転変動波形の偏差データの例を示す図で、燃料着火が上死点後５度の場合の例である。
【図１２】クランク軸の回転変動波形の偏差データの例を示す図で、燃料着火が上死点後１０度の場合の例である。
【図１３】回転変化開始時期と燃料着火時期との関係を示す燃料着火時期マップの例を示す図である。
【符号の説明】
１    エンジンの燃焼制御装置
１０    回転変動検出手段
２０    基準データ作成手段
２２    数式モデル作成手段
３０    燃焼状態推定手段
４０    制御パラメータ補正手段
５０    燃焼制御手段
Ｓ１    回転変動検出ステップ
Ｓ２    基準作成ステップ
Ｓ２３  数式モデル作成ステップ
Ｓ３    推定ステップ
Ｓ４    補正ステップ
Ｓ５    制御ステップ
Ｒｄｋ    サンプリングしたデータ
Ｒｄtdc   代表値
Ｒｍｎ    測定データ
Ｒｍmin   最低値
Ｒｍtdc   上死点における測定データ
Ｒｓｎ    基準データ
ΔＲｄｎ  偏差データ
ΔＲＣ１  第１閾値（変動値）
ΔＲＣ２  第２閾値（固定値）
Ｙ＝Ｆ（θ）  数式モデル
Ｔｆ，ＴｆＳ    燃料着火時期
Ｔｖ，ＴｖＳ  回転変動開始時期
Ｔｏ  燃料噴射発令時期
【図１】
0004075312.tif 000002
【図２】
0004075312.tif 000003
【図３】
0004075312.tif 000004
【図４】
0004075312.tif 000005
【図５】
0004075312.tif 000006
【図６】
0004075312.tif 000007
【図７】
0004075312.tif 000008
【図８】
0004075312.tif 000009
【図９】
0004075312.tif 000010
【図１０】
0004075312.tif 000011
【図１１】
0004075312.tif 000012
【図１２】
0004075312.tif 000013
【図１３】
0004075312.tif 000014