(57)【特許請求の範囲】
【請求項１】
内燃機関を制御するための制御装置であって、
前記内燃機関の燃焼室に設けられた圧力検出手段と、
前記圧力検出手段の出力をサンプリングするサンプリング回路と、
前記サンプリング回路が生成するサンプル値を移動平均した値に対して該移動平均による遅れを補正して平滑化したサンプル値を算出する平滑化手段と、
前記平滑化したサンプル値に基づいて検出される前記筒内圧力が最大となる時点のクランク角度θPmaxを検出するθPmax検出手段と、
燃料カット運転中に前記θPmax検出手段によって検出される前記θPmaxとモータリング圧力における上死点のクランク角度との差分を検出する差分検出手段と、
を備える、制御装置。
【請求項２】
前記サンプリング回路は、サンプリング開始時点を、前記差分検出手段により検出される差分に基づいて補正する、請求項１記載の制御装置。
【請求項３】
前記サンプリング回路は、前記圧力検出手段のサンプリング開始時点を、エンジン回転数に応じた前記圧力検出手段の出力遅延に応じて遅延させる、請求項１又は２に記載の制御装置。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、クランク角センサにおける圧縮トップ位置のずれを筒内圧センサを用いて検出する技術に関する。
【背景技術】
【０００２】
内燃機関の制御において、クランク角センサの出力は、制御タイミングの基準となるものである。例えば正確な失火検知や着火遅れ検出などを行うためには、クランク角センサの出力が正確であることが望まれる。しかし、このクランク角センサの出力は経年劣化等により狂いが生じることがあり、狂いが生じた場合は上記のような検出を精度よく行うことができない。
【０００３】
クランク角センサの出力を補正する技術が特許文献１に提案されている。特許文献１には、クランク角位置検出手段と、機関燃焼室内の圧力を検出する圧力センサと、圧力センサの出力から機関モータリング時の最高圧力を判定してそのときのクランク角位置を求める手段と、クランク角位置検出手段によるクランク角位置と圧力センサ出力からクランク角位置を求める手段により求められたクランク角位置との相対的なずれ量を演算する手段とを備えた内燃機関のクランク角位置検出補正装置が記載されている。
【特許文献１】特開昭60-114734号
【発明の開示】
【発明が解決しようとする課題】
【０００４】
上記の従来技術では、圧力センサの出力から燃料カット運転中（モータリング時）の最高圧力を判定して、そのときのクランク角位置から相対的なずれ量を演算している。しかしながら、燃料カット運転中の筒内圧は比較的小さくて正確な値を得ることが困難であり、上記従来技術では、圧力センサへ混入されるノイズ等の影響により、モータリング時の最高圧力時におけるクランク角位置の正確な検出が不可能である。
【０００５】
そこで本発明の課題は、燃料カット中にクランク角センサにおける圧縮トップ位置のずれを精度よく検出することができる技術を提供することにある。
【課題を解決するための手段】
【０００６】
この発明の制御装置は、一形態（請求項１）では、内燃機関の燃焼室に設けられた圧力検出手段と、圧力検出手段の出力をサンプリングするサンプリング回路と、サンプリング回路が生成するサンプル値を平滑化する平滑化手段と、平滑化したサンプル値に基づいて検出される筒内圧力が最大となる時点のクランク角度θPmaxを検出するθPmax検出手段と、燃料カット運転中にθPmax検出手段によって検出されるθPmaxとモータリング圧力における上死点のクランク角度との差分を検出する差分検出手段と、を備えるものである。
【０００７】
この発明によると、燃料カット中にクランク角センサにおける圧縮トップ位置ずれを精度よく検出することができる。
【０００８】
この発明の制御装置は、一形態（請求項２）では、サンプリング回路は、サンプリング開始時点を、差分検出手段により検出される差分に基づいて補正するよう構成されている。この発明によると、クランク角センサの狂いやその他の各種の要因によって発生するクランク角位置のずれを補正することができる。
【０００９】
この発明の制御装置は、一形態（請求項３）では、平滑化手段は、サンプル値を移動平均することによる遅れを補正するよう構成されている。この発明によると、移動平均による遅れに起因したクランク角位置のずれが発生しない。
【００１０】
この発明の制御装置は、一形態（請求項４）では、サンプリング回路は、前記圧力検出手段のサンプリング開始時点を、エンジン回転数に応じた前記圧力検出手段の出力遅延に応じて遅延させるよう構成されている。この発明によると、ローパスフィルタ処理による遅延の影響のない圧力サンプル値を得ることができる。
【発明を実施するための最良の形態】
【００１１】
次に図面を参照して、この発明の実施例を説明する。図１は、この発明の制御装置の全体的構成を示すブロック図である。電子制御ユニット１０は、中央演算装置（CPU）を備えたコンピュータである。電子制御ユニットは、コンピュータ・プログラムを格納する読取専用メモリ（ROM）およびプロセサに作業領域を提供し、データおよびプログラムを一時記憶するランダムアクセス・メモリ（RAM）を備えている。入出力インタフェイス１１は、エンジンの各部から検出信号を受け取って、A/D（アナログ・ディジタル）変換を行って次の段階に渡す。また、入出力インタフェイス１１は、CPUの演算結果に基づく制御信号をエンジンの各部に送る。図１では、電子制御ユニットをこの発明に関連する機能を示す機能ブロックで示している。
【００１２】
図１において、筒内圧力センサ１２は、圧電素子であり、エンジンの各気筒（シリンダ）の点火プラグ付近に設けられている。圧力センサ１２は、気筒内の圧力に応じた電荷信号を出力する。この信号をチャージアンプ３１により電圧信号に変換して出力し、ローパスフィルタ３３を介して入出力インタフェイス１１に出力する。入出力インタフェイス１１は、圧力センサ１２からの信号をサンプリング部１３に送る。また、入出力インタフェイス１１は、クランク角センサ１８の出力もサンプリング部１３に送る。
【００１３】
まず図面２を参照して、筒内圧力センサ１２の特性変化を相殺するための補正について説明する。この補正処理は、特に本発明の要部に関連するものではない。ここでは、補正量を通常運転時に決定するものとして説明する。
【００１４】
図２は、通常運転時におけるクランク角度-180度から180度の領域における気筒の燃焼室の圧力を示しており、およそクランク角度-180度から０度の範囲が圧縮行程であり、０度から180度までが膨張（燃焼）行程である。曲線１は、エンジンの１つの気筒のモータリング圧力（燃焼がない時の圧力）の推移を示し、曲線３は、同じ気筒において正常な燃焼が行われたときの筒内圧力の推移を示す。クランク角０度が上死点であり、モータリング圧力は上死点でピークとなり、燃焼時の筒内圧力（曲線３）は、上死点手前で点火後、着火を経て上死点を少し過ぎた時点でピークとなる。
【００１５】
まず、圧縮行程において上死点に達する前の期間、たとえば図２に”a”で示す期間において、筒内圧力検出手段（図１の筒内圧力センサ１２）の検出出力を補正する補正式のパラメータを同定する。黒色のドット５は、筒内圧力センサ１２による検出出力を示す。筒内圧力センサ１２は、エンジンの燃焼室という過酷な環境に置かれており、温度の影響、経年変化などによって特性が変化する。そこで、筒内圧力センサ１２の検出出力がほぼモータリング圧力の曲線１上にくるよう、検出出力を補正する。こうして補正された検出出力を白色のドット７で示す。
【００１６】
検出出力の補正は、筒内圧力センサの検出出力PD（θ）に、補正式  PS = PD（θ）k1 + C1 を適用することによって行われる。k1  は補正係数であり、C1 は定数である。この補正式の２つのパラメータk1およびC1は、圧縮行程の、たとえば図２に”a”で示す期間において、モータリング圧力の推定値PMと、筒内圧力センサの検出出力を上述の補正式によって補正した値PSとの差(PM - PS)の二乗が最小になるよう、最小二乗法により演算して算出される。
【００１７】
センサ出力補正部１７は、上述の補正式  PS = PD（θ）k1 + C1 に従って、センサ出力PD（θ）を補正する。センサ出力補正部１７は、補正されたセンサ出力値PSをθＰｍａｘ算出部４１に渡す。
【００１８】
一方において、燃焼室容積計算部１９は、クランク角θに応じた気筒の燃焼室の容積Vcを次の数式により計算する。
【数１】
0004684944.tif 000002
【００１９】
上の式で、mは、図３の関係から計算される、ピストン8の上死点からの変位を示す。rをクランク半径、l（エル）をコンロッド長とすると、λ=l/r である。Vdeadは、ピストンが上死点にあるときの燃焼室の容積、Apstnは、ピストンの断面積である。
【００２０】
一般に燃焼室の状態方程式は、次の（３）式で表されることが知られている。
【数２】
0004684944.tif 000003
【００２１】
(3)式で、Gは、例えばエアフローメータ、またはエンジン回転数および吸気圧に基づいて得られる吸入空気量、Rは気体定数、Tは、例えば吸気温度センサ、またはエンジン水温などの運転状態に基づいて得られる吸気温度である。ｋは補正係数、Cは定数である。
【００２２】
ここでは、燃焼室の気体状態方程式に基づくモータリング圧力推定値として、たとえばセンサ取り付け部の温度変化等の影響を受けない水晶圧電式の圧力センサを用いて燃焼室の圧力を実測し、この実測値を（３）式と対応させることにより算出されたｋの値ｋ0およびＣの値Ｃ0を求め、これを（３）式に代入して得られる次の（４）式を用いてモータリング圧力を推定する。
【数３】
0004684944.tif 000004
【００２３】
モータリング圧力推定部２０は、基本モータリング圧力計算部２１およびモータリング圧力補正部２２から構成される。基本モータリング圧力計算部２１が（３）式の中の基本項目である基本モータリング圧力GRT/Vを計算する。モータリング圧力補正部２２は、上述のようにして予め求められているパラメータk0およびC0を用いて、基本モータリング圧力を補正する。このパラメータk0およびC0は、吸気管圧力またはエンジン回転数などエンジンの負荷状態を表すパラメータに従って参照することができるテーブルとして用意されている。
【００２４】
なお、モータリング圧力推定部２０は、代替的に、基本モータリング圧力計算部２１のみで構成する形式でもよい。この場合、モータリング圧力ＰＭは、基本モータリング圧力計算部２１が算出した基本モータリング圧力ＧＲＴ／Ｖとなる。
【００２５】
パラメータ同定部２３は、圧縮行程においてモータリング圧力推定部２０が算出するモータリング圧力推定値PMとセンサ出力補正部１７が出力する筒内圧力センサ１２に基づく筒内圧力PSとの誤差(PM-PS)が最小になるよう、最小二乗法によりセンサ出力を補正する補正式のパラメータk1およびC1を同定する。センサ出力検出部１５は、たとえば10kHz分の１の周期で圧力センサの出力をサンプリングし、クランク角度に同期したタイミングでサンプル値の平均値をセンサ出力値PD（θ）として、パラメータ同定部２３に渡す。パラメータ同定部２３は、気筒の圧縮行程において補正式のパラメータを同定する演算を実行する。モータリング圧力補正部から得られるクランク角度に応じたモータリング圧力推定値PM（θ）と、同じクランク角度におけるセンサ出力値PD（θ）に補正式PS = PD（θ）k1 + C1 を適用した値PSとの差の二乗、すなわち（PM（θ） - PD（θ）k1 - C1）2  が最小になる k1 および C1を既知の最小二乗法により求める。
【００２６】
PMの離散値をy(i)で表し、筒内圧力センサから得られる筒内圧力PDのサンプル値（離散値）をx(i)で表すと、P’T=[p’(0),p’(1), …,p’(n)] 、PT=[p(0), p(1), …,p(n)]、X(i)T=[x(0), x(1), …,x(n)]と表される。誤差（P’-P）の離散値の二乗の和は、次の式（５）で表される。サンプル値は、10kHz分の１の周期でとられ、iの値は、たとえば100までの値であり、たとえば10としてもよい。
【数４】
0004684944.tif 000005
【００２７】
このFの値を最小にするkおよびCを求めるには、F(k,C)のkおよびCに関する偏微分が０となるkおよびCを求めればよい。これを数式で表すと、次のようになる。
【数５】
0004684944.tif 000006
【００２８】
式(6)および(7)の右辺を整理すると、次のようになる。
【数６】
0004684944.tif 000007
【００２９】
これを行列で表現すると、次のようになる。
【数７】
0004684944.tif 000008
この式を逆行列を使って変形すると、次のようになる。
【数８】
0004684944.tif 000009
【００３０】
ここで、右辺の逆行列は、次の式で表される。
【数９】
0004684944.tif 000010
【００３１】
補正部１７は、こうして同定されたパラメータを用いてセンサ出力PD（θ）を補正する。補正された所定クランク角ごとのセンサ出力PS（θ）がθＰｍａｘ算出部４１に渡される。このように補正することによって、θＰｍａｘ算出部４１に入力される際には、筒内圧力センサ１２の特性変化が相殺される。以上が筒内圧力センサ１２の特性変化を相殺するための補正についての説明である。
【００３２】
次に、本発明の概要について図１を参照しながら説明する。図１において、サンプリング部１３は、入出力インタフェイス１１から渡される信号を所定の周期、たとえば１０kHz分の１の周期でサンプリングする。サンプリング部１３は、ローパスフィルタ３３のフィルタ処理による筒内圧力センサ１２のセンシング遅れ分だけサンプリング開始タイミングを遅らせてセンシング遅れを補正する（第１の補正処理）。
【００３３】
また、サンプリング部１３は、燃料カット運転時以外の時には、上記第１の補正処理に加え、後述する燃料カット運転時に検出した圧縮トップ位置ずれの分をさらに補正している（第２の補正処理）。この燃料カット運転時の圧縮トップ位置ずれを、燃料カット運転時に検出できるように、サンプリング部１３は、第１の補正処理を加えた圧力センサのサンプル値をセンサ出力検出部１５とセンサ出力補正部１７とを介してθＰｍａｘ算出部４１に生成したサンプル値を送っている。
【００３４】
θＰｍａｘ算出部４１は、燃料カット運転時に筒内圧力センサ１２に混入されるノイズに対してより精度よく圧力を検出するために、筒内圧力センサ１２の出力、すなわちサンプリング部１３から出力されるサンプル値の平滑化（移動平均処理）を行ってから最高圧力クランク角位置の検出を行なう。このとき、併せてクランク角位置に対して平滑化処理による遅れ角度の補正（第３の補正処理）行い、検出した最高圧力クランク角位置をトップ位置ずれ算出部４３に送る。
【００３５】
トップ位置ずれ算出部４３は、燃料カット運転時に、平滑化処理後の圧力サンプル値から求められた最高圧力クランク角位置と、所定の上死点位置（圧縮トップ位置のクランク角位置、例えば３６０度）との「ずれ」を算出する。また、算出した圧縮トップ位置のずれを、サンプリング部１３に送る。
【００３６】
次に、本発明の処理の詳細について図５、図６のフローチャートに沿って説明する。図５は、クランク角位置の補正処理を示すフローチャートであり、図６は、圧縮トップ位置を検出する処理を示すフローチャートである。
【００３７】
まず、サンプリング部１３は、ローパスフィルタ処理によるセンシング遅れについての補正（第１の補正処理）を行うために、エンジン回転数を受け取り、内蔵する図４に示すテーブルからエンジン回転数ｋに対応する筒内圧力センサ１２のサンプリング取得開始ポイントの吸気トップからのオフセット角度Δｋ１を取得する（図７参照）。図４に示すテーブルには、エンジン回転数ｋと、オフセット角度ｋ１との対応関係が予め記憶されており、オフセット角度Δｋ１には、エンジン回転数に対応したローパスフィルタ処理によるセンシング遅れが含まれている。
【００３８】
図７は燃料カット運転の時におけるローパスフィルタ処理による筒内圧力センサ１２のセンシング遅れ角度がない場合の理想的なモデル波形（モータリング圧力の波形）を示す曲線Ｍと、燃料カット運転時における前述のセンシング遅れがある場合の筒内圧力センサ１２の実際の検出出力を示す曲線Ｃを示している。
【００３９】
サンプリング部１３は、エンジン回転数に対応するローパスフィルタ処理によるセンシング遅れが含まれるオフセット角度Δｋ１に対応する時点Ｔ２をサンプリング開始時点として特定し（Ｓ１０１）、ローパスフィルタ３３を経由した筒内圧力センサ１２の出力のサンプリングを開始する。オフセット角度Δｋ１には、エンジン回転数に対応するローパスフィルタ処理によるセンシング遅れについての補正（第１の補正処理）が行われることになる。
【００４０】
さらにサンプリング部１３は、モータリング圧力推定部２０からモータリング圧力（モデル波形）の算出開始ポイントＴ１の吸気トップからのオフセット角度Δｋ２を受け取る。このモデル波形についてのオフセット角度Δｋ２と、先のオフセット角度Δｋ１との差分をモデル波形の算出開始ポイントＴ１からのオフセット量Ａとして算出している（Ｓ１０３）。
【００４１】
また、クランク角センサ１８の狂いやその他の各種要因によって、モデル波形の曲線Ｍと実際の検出出力の曲線Ｃとの間にずれが生じており、そのためさらに以下の処理を行なっている。
【００４２】
サンプリング部１３は、燃料カット運転中であるか否かを判断し（Ｓ１０５）、燃料カット運転中であるならば（Ｓ１０５：Ｙｅｓ）、前述した第１の補正処理を行ったサンプリング出力をセンサ出力検出部１５とセンサ出力補正部１７とを介してθＰｍａｘ算出部４１に送る。
【００４３】
θＰｍａｘ算出部４１は、受け取った第１の補正処理後のサンプリング出力を平滑化処理した後、筒内圧力最大時のクランク角位置を図６に示すサブフローに従って算出する（Ｓ１０７）。まず、θＰｍａｘ算出部４１は、移動平均を算出する（Ｓ１２１）。具体的には、受け取った検出出力のサンプル値Ｐ１を平滑化をするために、受け取った検出出力のサンプル値Ｐ１に次式（１１）を適用して移動平均後のサンプル値Ｐ２を算出する。
【数１０】
0004684944.tif 000011
【００４４】
上記のように移動平均を算出することによって、図８に示すように移動平均前のサンプリング出力を示す曲線Ｆは平滑化され、曲線Ａのように平滑化される。さらに移動平均後のサンプリング出力に基づいて筒内圧力最大のときのクランク角位置である圧縮トップ位置Ａθを検出する。この圧縮トップ位置Ａθの検出を複数サイクルにわたって実行する（Ｓ１２３）。
【００４５】
さらに、この複数サイクルについての圧縮トップ位置を平均して得られるクランク角位置を移動平均による圧縮トップ位置Ｈθとして特定する（Ｓ１２５）。
【００４６】
また、図８の下図は、移動平均前のサンプリング出力に基づく圧縮トップ位置の検出頻度分布と、移動平均前のサンプリング出力に基づく圧縮トップ位置の検出頻度分布を示すものである。移動平均を取ることにより、圧縮トップ位置の検出値のばらつき低減効果を得ることができるが、その反面所定の検出遅れ角度（ずれ量）Ｄが発生する。
【００４７】
そこで、θＰｍａｘ算出部４１は、特定した圧縮トップ位置Ｈθと移動平均によるずれ量Ｄとを次式（１２）に適用することにより遅れ分を補正（第３の補正処理）して、正しい圧縮トップ位置θＰｍａｘを算出する（Ｓ１２９）。算出した圧縮トップ位置θＰｍａｘはトップ位置ずれ算出部４３に送られる。
【数１１】
0004684944.tif 000012
【００４８】
以上のように圧縮トップ位置θＰｍａｘが算出されると、図５に戻って、トップ位置ずれ算出部４３は、受け取った圧縮トップ位置θＰｍａｘとモータリング圧力における上死点のときの所定のクランク角位置θＰＭｍａｘ（たとえば３６０度）とを次式（１３）に適用することによって気筒毎の圧縮トップ位置のずれ量δθを算出する（Ｓ１１１）。
【数１２】
0004684944.tif 000013
【００４９】
圧縮トップ位置のずれ量δθの算出は、全気筒について算出される。燃料カット運転時は燃焼行程がないので、正確な検出ができていれば圧縮トップ位置（上死点）のクランク角位置はモータリング圧力のものと一致するはずである。したがって、モータリング圧力が上死点となるときのクランク角位置との差がトップ位置のずれ量δθとして算出される。
【００５０】
さらに、トップ位置ずれ算出部４３は、算出した各気筒の圧縮トップ位置のずれ量δθを用いて全気筒についての平均の圧縮トップ位置のずれ量Δθを算出する（Ｓ１１３）。全気筒についての平均の圧縮トップ位置のずれ量Δθは、例えば、算出した各気筒ごとの圧縮トップ位置のずれ量について最大値と最小値を除いてから平均をとることによって算出することができる。算出した圧縮トップ位置のずれ量Δθは、サンプリング部１３に送られる。
【００５１】
一方、図５のＳ１０５に戻って、サンプリング部１３において、燃料カット運転中でないと判断されると（Ｓ１０５：ＮＯ）、上記の燃料カット時に算出したオフセット量の補正量Δθ（＝Ｂ）とオフセット量Ａ（＝Δｋ１−Δｋ２）とを次式（１４）に適用することにより最終的なオフセット量Ｃ（Ｃ＝Ｔ３―Ｔ１）を算出する。
【数１３】
0004684944.tif 000014
【００５２】
したがって燃料カット運転中でない場合は、サンプリング部１３は、吸気トップからΔｋ２＋Ｃの時点Ｔ３、あるいはモデル波形のサンプリング開始点Ｔ１からオフセット量Ｃだけオフセットしたクランク角位置Ｔ３をサンプリング開始時点として、ローパスフィルタを経由した筒内圧力センサ１２の出力のサンプリングを開始する。このようにサンプリング開始点をずらすことによって、それ以降のサンプリング時点もずらすことになる。
【００５３】
以上にこの発明を具体的な実施例について説明したが、この発明はこのような実施例に限定されるものでなく、また、ガソリンエンジン、ディーゼルエンジンのいずれに対しても使用することができる。
【図面の簡単な説明】
【００５４】
【図１】この発明の実施形態の機能ブロック図。
【図２】モータリング圧力曲線および点火を生じたときの圧力曲線を表す図。
【図３】ピストン位置を計算するための概念図。
【図４】回転数とサンプリング開始点のクランク角との設定値を示す図。
【図５】正確なオフセット角を設定する処理を示すフローチャート。
【図６】トップ位置クランク角の算出方法を示すフローチャート。
【図７】燃料カット運転時におけるトップ位置ずれについて説明する図。
【図８】移動平均後のトップ位置検出を説明する図。
【符号の説明】
【００５５】
１０      電子制御ユニット（ECU）
１２      筒内圧力センサ
１５      センサ出力検出部
１９      燃焼室（シリンダ）容積計算部
２０      モータリング圧力推定部
４１      θPmax算出部
４３      トップ位置ずれ算出部
【図１】
0004684944.tif 000015
【図２】
0004684944.tif 000016
【図３】
0004684944.tif 000017
【図４】
0004684944.tif 000018
【図５】
0004684944.tif 000019
【図６】
0004684944.tif 000020
【図７】
0004684944.tif 000021
【図８】
0004684944.tif 000022