(57)【要約】【課題】この発明は、出力軸の回転速度が変動する場合でも、目標の回転角が到来するタイミングを高い精度で予測することを目的とする。【解決手段】内燃機関の燃焼を停止した状態において、機械構造に起因するクランク軸３２の回転変動である機械性回転変動ΔＮc(i)を予め計測しておく。ＥＣＵ４０は、内燃機関の運転中に計測した実回転変動ΔＮr(i)から機械性回転変動ΔＮc(i)を減算することにより、燃焼に起因する燃焼性回転変動ΔＮb(i)を算出し、その算出値を用いて将来の燃焼性回転変動を予測する。そして、クランク軸３２の回転速度Ｎと、機械性回転変動ΔＮc(i)と、燃焼性回転変動ΔＮb(i)の予測値とを用いて、現在のクランク角から所望の目標角度までの回転に必要な目標時間Ｔを予測する。これにより、各種の制御を所望のタイミングで正確に実行することができる。【選択図】図７
2009174483.tif 000002
【特許請求の範囲】
【請求項１】
内燃機関の出力軸が回転するときに、当該出力軸の回転に応じた検出信号を出力する回転検出手段と、
前記回転検出手段の検出信号を用いて前記出力軸の回転速度を算出する回転速度算出手段と、
内燃機関の運転中に前記出力軸の回転速度が単位時間当りに変化する変化量を、実回転変動として算出する実変動算出手段と、
内燃機関の燃焼を停止している状態で前記出力軸の回転速度が単位時間当りに変化する変化量を、内燃機関の機械構造に起因する機械性回転変動として取得する機械性変動取得手段と、
前記実回転変動から前記機械性回転変動を除いた残りの回転変動を、内燃機関の燃焼に起因する燃焼性回転変動として算出する燃焼性変動算出手段と、
前記燃焼性回転変動が算出される毎に当該算出値を記憶する記憶手段と、
前記記憶手段に記憶された過去の燃焼性回転変動に基づいて、将来の燃焼性回転変動を予測する燃焼性変動予測手段と、
前記回転速度と、前記機械性回転変動と、前記燃焼性回転変動の予測値とを用いて、前記出力軸の回転角が現在の角度から目標角度となるまでに必要な目標時間を予測する時間予測手段と、
を備えることを特徴とする内燃機関の制御装置。
【請求項２】
前記記憶手段は、前記出力軸の回転角と前記回転速度とに応じて変化する前記機械性回転変動の特性が予め記憶された機械性変動データを備え、
前記機械性変動取得手段は、前記機械性変動データを参照することにより、任意の回転角及び回転速度における前記機械性回転変動を取得する構成としてなる請求項１に記載の内燃機関の制御装置。
【請求項３】
前記燃焼性変動予測手段は、現在または過去の時点で算出された燃焼性回転変動を基準として、この時点から１燃焼サイクルが経過した後の燃焼性回転変動を予測する構成としてなる請求項１または２に記載の内燃機関の制御装置。
【請求項４】
前記燃焼性変動予測手段は、過去の２つの時点における燃焼性回転変動と、前記２つの時点のうち一方の時点から１燃焼サイクルが経過した後の燃焼性回転変動とを用いて、他方の時点から１燃焼サイクルが経過した後の燃焼性回転変動を予測する構成としてなる請求項１または２に記載の内燃機関の制御装置。
【請求項５】
前記記憶手段は、前記実回転変動が零となる回転角である無変動角度と前記回転速度との関係が予め記憶された無変動角データを備え、
前記回転速度算出手段は、過去の回転速度に応じて前記無変動角データを参照することにより、当該回転速度における無変動角度を取得し、前記出力軸の回転角が前記無変動角度となったときに最新の回転速度を算出する構成としてなる請求項１乃至４のうち何れか１項に記載の内燃機関の制御装置。
【請求項６】
前記燃焼性変動算出手段は、将来の燃焼性回転変動を一定の回転角毎に予測する構成とし、前記時間予測手段は、前記回転速度、前記機械性回転変動及び前記燃焼性回転変動の予測値を合計することにより将来の予測回転速度を前記一定の回転角毎に算出し、当該予測回転速度の逆数を現在の角度から前記目標角度まで積算することにより前記目標時間を算出する構成としてなる請求項１乃至５のうち何れか１項に記載の内燃機関の制御装置。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、例えばクランク角センサの検出信号に基づいて各種制御の実行タイミングを定めるのに好適に用いられる内燃機関の制御装置に関する。
【背景技術】
【０００２】
従来、例えば特許文献１（特開２０００−８９４１号公報）に開示されているように、クランク角センサの検出信号に基づいて、各気筒での燃料噴射制御、点火時期制御等を行う構成とした内燃機関の制御装置が知られている。この種の従来技術による内燃機関の制御装置は、クランク角センサの検出信号に基づいてクランク軸の回転角（クランク角）や回転速度を検出する。そして、これらの検出結果や内燃機関の運転状態等に応じて設定される目標のクランク角が到来したときに、燃料噴射や点火などを実行するものである。
【０００３】
この場合、クランク角の検出・演算処理には一定の時間が必要となるので、目標のクランク角を検出した時点で制御を開始したのでは、特に高回転時などに制御の遅れが生じ易い。このため、従来技術では、所定の基準クランク角（例えば、圧縮上死点など）が到来したときに回転速度を検出し、このときの回転速度等に基づいて目標のクランク角が到来するタイミングを事前に予測する構成としている。
【０００４】
【特許文献１】特開２０００−８９４１号公報
【発明の開示】
【発明が解決しようとする課題】
【０００５】
ところで、上述した従来技術では、例えば基準クランク角で検出した回転速度に基づいて、目標のクランク角が到来するタイミングを予測する構成としている。しかし、クランク軸には、各気筒の爆発行程で発生するトルクが間欠的に付加されるだけでなく、クランク軸自体やこれと連動する部品に加わる慣性力、重力等が作用するので、クランク軸の回転速度は１回転の間に複雑に変動する。
【０００６】
このため、従来技術では、基準クランク角においてタイミングを予測したとしても、このタイミングが到来するまでの間に生じる回転変動によって実際に制御が行われるタイミングが目標のクランク角からずれてしまい、正確なタイミングで制御を実行することができないという問題がある。
【０００７】
この発明は、上述のような課題を解決するためになされたもので、出力軸の回転速度が変動する場合でも、目標の回転角が到来するタイミングをより正確に予測することができ、各種の制御をより正確なタイミングで実行することが可能な内燃機関の制御装置を提供することを目的としている。
【課題を解決するための手段】
【０００８】
第１の発明は、内燃機関の出力軸が回転するときに、当該出力軸の回転に応じた検出信号を出力する回転検出手段と、
前記回転検出手段の検出信号を用いて前記出力軸の回転速度を算出する回転速度算出手段と、
内燃機関の運転中に前記出力軸の回転速度が単位時間当りに変化する変化量を、実回転変動として算出する実変動算出手段と、
内燃機関の燃焼を停止している状態で前記出力軸の回転速度が単位時間当りに変化する変化量を、内燃機関の機械構造に起因する機械性回転変動として取得する機械性変動取得手段と、
前記実回転変動から前記機械性回転変動を除いた残りの回転変動を、内燃機関の燃焼に起因する燃焼性回転変動として算出する燃焼性変動算出手段と、
前記燃焼性回転変動が算出される毎に当該算出値を記憶する記憶手段と、
前記記憶手段に記憶された過去の燃焼性回転変動に基づいて、将来の燃焼性回転変動を予測する燃焼性変動予測手段と、
前記回転速度と、前記機械性回転変動と、前記燃焼性回転変動の予測値とを用いて、前記出力軸の回転角が現在の角度から目標角度となるまでに必要な目標時間を予測する時間予測手段と、
を備えることを特徴とする。
【０００９】
第２の発明によると、前記記憶手段は、前記出力軸の回転角と前記回転速度とに応じて変化する前記機械性回転変動の特性が予め記憶された機械性変動データを備え、
前記機械性変動取得手段は、前記機械性変動データを参照することにより、任意の回転角及び回転速度における前記機械性回転変動を取得する構成としている。
【００１０】
第３の発明によると、前記燃焼性変動予測手段は、現在または過去の時点で算出された燃焼性回転変動を基準として、この時点から１燃焼サイクルが経過した後の燃焼性回転変動を予測する構成としている。
【００１１】
第４の発明によると、前記燃焼性変動予測手段は、過去の２つの時点における燃焼性回転変動と、前記２つの時点のうち一方の時点から１燃焼サイクルが経過した後の燃焼性回転変動とを用いて、他方の時点から１燃焼サイクルが経過した後の燃焼性回転変動を予測する構成としている。
【００１２】
第５の発明によると、前記記憶手段は、前記実回転変動が零となる回転角である無変動角度と前記回転速度との関係が予め記憶された無変動角データを備え、
前記回転速度算出手段は、過去の回転速度に応じて前記無変動角データを参照することにより、当該回転速度における無変動角度を取得し、前記出力軸の回転角が前記無変動角度となったときに最新の回転速度を算出する構成としている。
【００１３】
第６の発明によると、前記燃焼性変動算出手段は、将来の燃焼性回転変動を一定の回転角毎に予測する構成とし、前記時間予測手段は、前記回転速度、前記機械性回転変動及び前記燃焼性回転変動の予測値を合計することにより将来の予測回転速度を前記一定の回転角毎に算出し、当該予測回転速度の逆数を現在の角度から前記目標角度まで積算することにより前記目標時間を算出する構成としている。
【発明の効果】
【００１４】
第１の発明によれば、内燃機関の運転中に生じる実際の回転変動を、機械性回転変動と燃焼性変動とに分けて扱うことができる。これにより、機械性回転変動については、内燃機関の燃焼を停止した状態で予め計測しておき、その計測データを制御装置等に記憶させることができる。この結果、内燃機関の運転中には、実際の計測を行わなくても、運転状態に応じた機械性回転変動を計測データから取得することができ、運転中の計測が難しい機械性回転変動を正確かつ容易に求めることができる。
【００１５】
このように取得した機械性回転変動を、内燃機関の運転中に計測した実回転変動から減算することにより、運転中の計測が難しい燃焼性回転変動をリアルタイムで正確かつ容易に算出することができる。そして、燃焼性回転変動の波形には、内燃機関の燃焼サイクルに同期した周期性が存在するので、この周期性を利用することにより、過去に算出された燃焼性回転変動を用いて将来の燃焼性回転変動を予測することができる。
【００１６】
これにより、時間予測手段は、出力軸が所望の角度まで回転するのに必要な目標時間を予測するときに、出力軸の回転速度と、将来の回転変動とを予測結果に反映させることができる。このため、例えば高回転での運転中や、回転速度が大きく変化するような運転状態でも、目標時間の予測を正確に行うことができる。従って、燃料噴射制御や点火時期制御等からなる各種の制御を所望のタイミングで正確に実行することができる。
【００１７】
第２の発明によれば、機械性変動取得手段は、出力軸の回転角と回転速度とに応じて機械性変動データを参照することにより、任意の回転角及び回転速度における機械性回転変動を容易に取得することができる。この機械性回転変動を実回転変動から減算することにより、任意の回転角及び回転速度における燃焼性変動を容易に算出することができる。
【００１８】
また、機械性変動データを用いることにより、まだ到来していない将来の回転角における機械性回転変動を自由に得ることができる。これにより、目標時間の予測処理には、将来の機械性回転変動と、将来の燃焼性回転変動とを予測結果にそれぞれ反映させることができ、高い予測精度を得ることができる。
【００１９】
第３の発明によれば、燃焼性回転変動の波形には、内燃機関の燃焼サイクルに同期した周期性が存在する。このため、現在または過去の時点で算出された燃焼性回転変動を基準として、この時点から１燃焼サイクルが経過した後の燃焼性回転変動を正確に予測することができる。
【００２０】
第４の発明によれば、燃焼性回転変動の波形には、燃焼サイクルに同期した周期性が存在するので、例えば過去の２つの時点における燃焼性回転変動の比率、差分等は、１燃焼サイクル後でもほぼ一定とみなすことができる。これにより、例えば現在または過去に算出された３つの燃焼性回転変動を用いて、将来の燃焼性回転変動を正確に予測することができる。
【００２１】
第５の発明によれば、無変動角度での計測を行うことにより、回転変動の影響を殆ど受けない平均的な回転速度を正確に計測することができる。これにより、回転変動が原因となって機械性回転変動の取得値に誤差が生じるのを防止でき、この取得値の精度を高めることができる。また、無変動角データは、回転速度に応じてマップデータ化されているので、例えば過去の回転速度に応じて無変動角データを参照することにより、出力軸の回転速度に関係なく、常に適切な無変動角度を得ることができる。
【００２２】
第６の発明によれば、時間予測手段は、出力軸の回転速度と、機械性回転変動と、燃焼性回転変動の予測値を合計することにより、将来の回転角における予測回転速度を算出することができ、この算出処理を一定の回転角毎に行うことができる。この場合、予測回転速度の逆数は、個々の回転角において出力軸が回転するときに必要な時間となる。このため、予測回転速度の逆数を現在の角度から目標角度まで積算することにより、目標時間を算出することができる。
【発明を実施するための最良の形態】
【００２３】
実施の形態１．
［実施の形態１の構成］
以下、図１乃至図６を参照しつつ、本発明の実施の形態１について説明する。まず、図１は、実施の形態１のシステム構成を説明するための全体構成図を示している。本実施の形態のシステムは、例えば多気筒型の内燃機関１０を備えている（１気筒のみ図示）。内燃機関１０は、吸入空気が各気筒内に向けて流入する吸気通路１２と、各気筒内で生じた排気ガスが流出する排気通路１４とを備えている。
【００２４】
吸気通路１２には、内燃機関の吸入空気量を検出するエアフローメータ１６と、エアフローメータ１６の下流側に配置された電子制御式のスロットル弁１８とが設けられている。スロットル弁１８は、アクセル開度等に応じてスロットルモータ２０により開，閉駆動され、その開度に応じて吸入空気量を増減させる。一方、排気通路１４には、排気ガスを浄化する触媒２２が設けられている。また、内燃機関１０の各気筒は、燃料を噴射する燃料噴射弁２４と、気筒内の混合気に着火する点火プラグ２６と、吸気弁２８、排気弁３０等を備えている。
【００２５】
また、内燃機関１０には、クランク軸（出力軸）３２の回転に応じた検出信号を出力する回転検出手段としてのクランク角センサ３４が設けられている。クランク角センサ３４は、例えば一般的に公知な非接触式のセンサ等からなり、クランク軸３２と共に回転するギヤ部品の歯を検出するように構成されている。
【００２６】
そして、クランク角センサ３４は、クランク軸３２が所定角度だけ回転する毎に、後述のＥＣＵ４０に検出信号を出力する。本実施の形態では、この所定角度が１°である場合を例に挙げて説明する。ＥＣＵ４０は、クランク角センサ３４の検出信号をカウントすることにより、クランク軸３２の回転角（クランク角）を検出することができる。
【００２７】
さらに、本実施形態のシステムは、ＥＣＵ(Electronic Control Unit)４０を備えている。ＥＣＵ４０は、ＲＯＭ、ＲＡＭ等の記憶回路と、タイマ機能とを備えたマイクロコンピュータ等からなり、本実施の形態の記憶手段を構成している。ＥＣＵ４０の入力側には、エアフローメータ１６、クランク角センサ３４等を含むセンサ系統が接続されている。
【００２８】
また、ＥＣＵ４０の出力側には、スロットルモータ２０、燃料噴射弁２４、点火プラグ２６等を含めて各種のアクチュエータが接続されている。そして、ＥＣＵ４０は、センサ系統の出力に応じて必要なアクチュエータを作動させることにより、燃料噴射制御、点火時期制御等を含む各種の制御を実行することができる。
【００２９】
この場合、燃料噴射制御では、例えば内燃機関の吸入空気量、運転者のアクセル操作等に応じて燃料の目標噴射量を算出する。また、内燃機関の機関回転数等に応じて燃料を噴射すべき適切なクランク角（目標角度）を算出し、この目標角度が到来するタイミングを、後述の時間予測制御によって予測する。そして、目標のタイミングが到来したときに、燃料噴射弁２４から燃料を噴射させる。
【００３０】
また、点火時期制御では、例えば内燃機関の負荷状態、機関回転数等に応じて点火を行うべき適切なクランク角（目標角度）を算出する。そして、燃料噴射制御の場合とほぼ同様に、目標角度に対応するタイミングが到来したときに、点火プラグ２６によって気筒内の混合気に着火する。
【００３１】
［時間予測制御］
時間予測制御は、クランク軸３２が基準の角度（例えば、現在のクランク角）から所望の目標角度まで回転するのに必要な時間を、目標時間Ｔとして予測するものである。この予測制御では、クランク軸３２の回転速度だけでなく、回転速度が単位時間当りに変化する変化量（以下、回転変動と称す）も考慮する。そこで、まず最初に、図２を参照しつつ、クランク軸３２の回転変動について説明する。
【００３２】
図２は、クランク軸の回転変動とクランク角との関係を示す特性線図である。図２中の実線，点線及び一点鎖線は、実回転変動ΔＮr，機械性回転変動ΔＮc及び燃焼性回転変動ΔＮbをそれぞれ示している。ここで、実回転変動ΔＮrとは、内燃機関の運転中に実際に計測される回転変動である。
【００３３】
クランク軸３２の回転変動は、内燃機関の燃焼行程（爆発行程）で発生するトルクや、構造的に偏心した部分に作用する慣性力、重力等が原因となって生じる。本実施の形態では、実回転変動ΔＮrを、発生原因が異なる機械性回転変動ΔＮcと燃焼性回転変動ΔＮbに分けて扱うことにより、目標時間Ｔを正確に予測することを特徴としている。
【００３４】
機械性回転変動ΔＮcとは、内燃機関の機械構造に起因する回転変動であり、内燃機関の燃焼を停止した状態でも発生する回転変動として定義される。クランク軸３２には、回転中心から偏心した部位が存在する上に、バランサ（カウンタウェイト）等の偏心した部品が取付けられている。これらの偏心構造に作用する慣性力、重力等がクランク軸３２の回転運動に影響を与えることにより、機械性回転変動ΔＮcが発生する。
【００３５】
一方、燃焼性回転変動ΔＮbとは、内燃機関の燃焼に起因する回転変動であり、燃焼行程で発生したトルクがクランク軸３２に付加されることにより発生するものである。この場合、燃焼性回転変動ΔＮbは、実回転変動ΔＮrから機械性回転変動ΔＮcを除いた残りの回転変動として定義することができる。
【００３６】
一般に、燃焼性回転変動ΔＮbは、内燃機関の運転状態（出力トルク）に応じて変化するので、機関運転中に実際の燃焼性回転変動ΔＮbを計測するのは難しい。これに対し、機械性回転変動ΔＮcは、内燃機関の燃焼を停止した状態で、モータ等を用いてクランク軸３２を回転させることにより、予め計測しておくことが可能である。
【００３７】
このため、時間予測制御では、内燃機関の運転中に計測した実回転変動ΔＮrから、予めデータ化しておいた機械性回転変動ΔＮcを減算することにより、燃焼性回転変動ΔＮbをリアルタイムに算出する。そして、過去の燃焼性回転変動ΔＮbに基づいて将来の燃焼性回転変動ΔＮbを予測することにより、目標時間Ｔの予測を正確に行うものである。以下、これらの処理について説明する。
【００３８】
（基準回転速度の算出）
この処理では、まず、目標時間Ｔを予測するときにベースの回転速度として用いられる基準回転速度Ｎを算出する。より詳しく述べると、内燃機関の運転中に特定のクランク角（以下、無変動角度と称す）が到来したときに、クランク軸３２の回転速度を計測し、その計測値を基準回転速度Ｎとして記憶する。
【００３９】
ここで、無変動角度とは、内燃機関の１燃焼サイクルに対応するクランク角の範囲（４サイクルの内燃機関では、１〜７２０°）のうち、実回転変動ΔＮrが零となる特定のクランク角である。無変動角度においては、例えば回転変動の要因が互いに打消し合う現象等が生じることにより、実回転変動ΔＮrが殆ど零となる。
【００４０】
図３は、ある回転速度における実回転変動とクランク角との関係を示す特性線図である。この図に示すように、実回転変動ΔＮrは、クランク角に応じて正値と負値との間で変動しつつ、複数の無変動角度において零となる。また、無変動角度は、クランク軸３２が一定の速度で回転しているときにほぼ一定の角度となるが、回転速度が変化することにより異なる角度となる。
【００４１】
このため、ＥＣＵ４０には、無変動角度と回転速度（基準回転速度）との関係をマップデータ化した無変動角データが予め記憶されている。図４は、無変動角データの一例を示すものである。なお、図４では、１つの基準回転速度について１つの無変動角度が得られるデータを例示したが、実回転変動ΔＮrが零となるクランク角は複数存在するので、１つの基準回転速度について複数の無変動角度が得られるようにしてもよい。
【００４２】
ＥＣＵ４０は、例えば前回の算出処理で得られた基準回転速度Ｎを用いて、図４の無変動角データを参照することにより、当該回転速度において回転変動が殆ど生じない無変動角度を取得することができる。そして、内燃機関の運転中にクランク角が無変動角度となったときには、クランク軸３２の回転速度を基準回転速度Ｎとして算出する。この場合、回転速度の算出は、クランク角センサ３４の検出信号を用いて、例えばクランク軸３２が１°回転する間の時間を計測することにより可能となる。
【００４３】
上記構成によれば、無変動角度での計測を行うことにより、回転変動の影響を殆ど受けない平均的な回転速度を、基準回転速度Ｎとして正確に得ることができる。この場合、無変動角データは、回転速度に応じてマップデータ化されているので、例えば前回の基準回転速度Ｎに応じて無変動角データを参照することにより、クランク軸３２の回転速度に関係なく、常に適切な無変動角度を得ることができる。
【００４４】
そして、基準回転速度Ｎを算出することにより、時間予測を行うときにベースとなる基準回転速度Ｎと、ベースの回転速度に対する回転変動ΔＮc(i)，ΔＮb(i)とを分けて扱うことができる。従って、これらのパラメータをそれぞれ適切な方法で目標時間Ｔに反映させることができ、目標時間Ｔを正確に予測することができる。
【００４５】
（実回転変動の算出）
この処理では、内燃機関の運転中にクランク軸３２に生じる実際の回転変動が、実回転変動ΔＮrとして算出される。具体的に述べると、ＥＣＵ４０は、クランク角センサ３４の検出信号を用いて、例えばクランク軸３２が１°回転する毎に回転速度を検出する。そして、時間的に隣接した２つの回転速度の差分を演算することにより、実回転変動ΔＮrを算出することができる。
【００４６】
また、上述の算出処理は、クランク軸３２が１°回転する毎に繰返される。そして、１°毎に得られる実回転変動ΔＮrの算出値は、例えば最新の１燃焼サイクル分の時系列データΔＮr(i)（i＝1,2,3,…,720）としてＥＣＵ４０に記憶される。ここで、ｉはクランク角を示すものである。
【００４７】
（機械性回転変動の取得）
機械性回転変動ΔＮcは、内燃機関の燃焼を停止した状態で予め計測しておくことが可能である。また、本願発明者の実験等により、機械性回転変動ΔＮcは、クランク角と回転速度に応じて変化することが判明した。このため、本実施の形態では、クランク角と回転速度とに応じて変化する機械性回転変動ΔＮcの特性をマップデータ化し、このマップデータを機械性変動データとしてＥＣＵ４０に予め記憶させる構成としている。
【００４８】
図５は、ある回転速度における機械性回転変動とクランク角との関係を示すマップデータである。このマップデータは、機械性変動データのうち当該回転速度に対応する一部のデータを構成している。即ち、機械性変動データは、図５のマップデータを回転速度に応じて複数種類用意したものである。この機械性変動データによれば、ある基準回転速度Ｎにおける機械性回転変動ΔＮcは、例えばクランク角１°毎に設定されたマップデータΔＮc(i)（i＝1,2,3,…,720）として与えられる。
【００４９】
（燃焼性回転変動の算出）
この処理では、下記（１）式に示すように、最新の１燃焼サイクルにおいて計測された実回転変動ΔＮr(i)から、個々のクランク角における機械性回転変動ΔＮc(i)を減算することにより、最新の１燃焼サイクルにおける燃焼性回転変動ΔＮb(i)（i＝1,2,3,…,720）を算出する。この算出処理を行うことにより、ＥＣＵ４０には、最新の１燃焼サイクル分の燃焼性回転変動ΔＮb(i)が常に記憶されるようになる。
【００５０】
ΔＮb(i)＝ΔＮr(i)−ΔＮc(i)  ・・・（１）
【００５１】
（将来の燃焼性回転変動の予測）
図６は、前記（１）式により得られた燃焼性回転変動ΔＮb(i)とクランク角との関係を示している。図６中に示すＡn，Ｂn（n＝1,2,3,4）は、＃ｎ気筒での爆発に起因する燃焼性回転変動の振幅のピーク値である。なお、図６では、４気筒型の内燃機関における燃焼性回転変動の波形を例示したが、本発明は４気筒に限定されるものではない。
【００５２】
この図に示すように、燃焼性回転変動ΔＮb(i)は、各気筒の爆発行程で発生したトルクがクランク軸３２に付加されることにより増大し、トルクが小さくなった時点で減少に転じる。このため、燃焼性回転変動ΔＮb(i)の波形は、各気筒の点火タイミングを起点としてプラス側とマイナス側に振動する。また、個々のピーク値Ａn，Ｂnは、各気筒での発生トルクに対応した大きさとなる。
【００５３】
個々の気筒で発生するトルクは、例えば１燃焼サイクル程度の短時間であれば、急変することがないと考えられる。従って、燃焼性回転変動も、短時間であれば、トルクに応じて比較的緩やかに変化するとみなすことができる。このように、燃焼状態の変化が緩やかな状況において、燃焼性回転変動の波形には、内燃機関の燃焼サイクルに同期した周期性が存在する。このため、本実施の形態では、現在または過去の時点で算出した燃焼性回転変動を基準として、この時点から１燃焼サイクルが経過した後の燃焼性回転変動を予測する構成としている。
【００５４】
具体的に述べると、爆発間隔が十分に短い２つの気筒（例えば、爆発行程が隣接した＃１，＃２気筒）に着目した場合に、ある時点での当該気筒における燃焼性回転変動のピーク値の比率（Ｂ1／Ｂ2）と、この時点から１燃焼サイクルが経過した後のピーク値の比率（Ｂ1′／Ｂ2′）とは、燃焼状態の急変がない状態であれば、ほぼ等しいとみなすことができる。即ち、ピーク値Ｂ1，Ｂ2，Ｂ1′，Ｂ2′の間には、下記（２）式が成立する。
【００５５】
Ｂ1／Ｂ2＝Ｂ1′／Ｂ2′  ・・・（２）
【００５６】
従って、上述した４つのピーク値Ｂ1，Ｂ2，Ｂ1′，Ｂ2′のうち、例えば３つのピーク値Ｂ1，Ｂ2，Ｂ1′を現在または過去のクランク角で算出された燃焼性回転変動とし、ピーク値Ｂ2′を将来のクランク角における燃焼性回転変動とすれば、将来の燃焼性回転変動を予測することができる。即ち、前記（２）式を下記（３）式のように変形することにより、現在または過去のピーク値Ｂ1，Ｂ2，Ｂ1′を用いて将来のピーク値Ｂ2′を算出することができる。
【００５７】
Ｂ2′＝Ｂ1′×Ｂ2／Ｂ1  ・・・（３）
【００５８】
このように、本実施の形態では、過去の２つの時点における燃焼性回転変動（ピーク値Ｂ1，Ｂ2）と、前記２つの時点のうち一方の時点から１燃焼サイクルが経過した後の燃焼性回転変動（ピーク値Ｂ1′）とを用いて、他方の時点から１燃焼サイクルが経過した後の燃焼性回転変動（ピーク値Ｂ2′）を予測するものである。
【００５９】
換言すれば、ＥＣＵ４０は、爆発間隔が近い２つの気筒の燃焼により生じた燃焼性回転変動の比率（Ｂ1／Ｂ2）が、１燃焼サイクル後にも一定であるとみなすことにより、前記２つの気筒の燃焼性回転変動の比率（Ｂ1／Ｂ2）と、前記２つの気筒のうち一方の気筒における１燃焼サイクル後の燃焼性回転変動（ピーク値Ｂ1′）とを用いて、他方の気筒における１燃焼サイクル後の燃焼性回転変動（ピーク値Ｂ2′）を予測する。
【００６０】
そして、ＥＣＵ４０は、上記の予測方法を用いることにより、例えば過去の燃焼性回転変動のデータΔＮb(i-719)〜ΔＮb(i-540)を用いて、これらのデータから１燃焼サイクル（７２０°）が経過した将来の燃焼性回転変動ΔＮb(i+1)〜ΔＮb(i+180)を予測し、その予測値を記憶しておく。
【００６１】
（目標時間の予測）
これまでの処理により、最新の基準回転速度Ｎと、任意のクランク角における機械性回転変動ΔＮc(i)と、将来のクランク角における燃焼性回転変動ΔＮb(i+1)〜ΔＮb(i+180)とが得られた。そこで、これらの値を下記（４）式で合計することにより、将来のクランク角における予測回転速度Ｎ(i+1)〜Ｎ(i+180)を算出することができる。
【００６２】
Ｎ(i+j)＝Ｎ＋ΔＮc(i+j)＋ΔＮb(i+j)  ・・・（４）
【００６３】
ここで、(ｉ+ｊ)とは、現在のクランク角ｉに対して角度ｊだけ進んだ将来のクランク角を示すものとする（本実施の形態では、j＝1,2,3...,180として例示）。即ち、予測回転速度Ｎ(i+j)とは、クランク軸３２が現在の位置から角度ｊだけ回転したときに生じているはずの回転速度を予測したものである。これにより、将来の任意のクランク角におけるクランク軸３２の回転速度を予測することができる。
【００６４】
ここで、予測回転速度の逆数１／Ｎ(i+j)は、個々の回転角において出力軸が１°回転するときに必要な時間となる。このため、クランク軸３２が現在の位置から目標角度αまで回転するのに必要な時間（目標時間Ｔ）を求める場合には、下記（５）式に示すように、予測回転速度の逆数１／Ｎ(i+j)を、現在のクランク角から目標角度αまでの範囲（ｊ＝１〜α）で積算する。具体的な例を挙げると、例えば目標角度が６０°である場合には、１／Ｎ(i+j)を、ｊ＝１〜６０の範囲で積算することにより、目標時間Ｔを算出することができる。
【００６５】
Ｔ＝１／Ｎ(i+1)＋１／Ｎ(i+2)＋１／Ｎ(i+3)＋ ... ＋１／Ｎ(i+α)  ・・・（５）
【００６６】
以上、詳述したように、本実施の形態では、実回転変動ΔＮr(i)を、機械性回転変動ΔＮc(i)と燃焼性回転変動ΔＮb(i)とに分けて扱うことができる。これにより、機械性回転変動ΔＮc(i)については、内燃機関の燃焼を停止した状態で予め計測しておくことができ、その計測値を機械性変動データとしてＥＣＵ４０に記憶させることができる。
【００６７】
この結果、内燃機関の運転中には、実際の計測を行わなくても、基準回転速度Ｎに応じて機械性変動データから機械性回転変動ΔＮc(i)を取得することができ、運転中の計測が難しい機械性回転変動ΔＮc(i)を正確かつ容易に求めることができる。この場合、基準回転速度Ｎは、回転変動の影響を受けない無変動角度において計測しているので、回転変動が原因となって機械性回転変動ΔＮc(i)の取得値に誤差が生じるのを防止でき、この取得値の精度を高めることができる。
【００６８】
また、機械性変動データには、基準回転速度Ｎとクランク角とに応じて変化する機械性回転変動ΔＮcの特性を記憶させる構成としている。このため、ＥＣＵ４０は、機械性変動データを用いて任意のクランク角及び基準回転速度における機械性回転変動ΔＮc(i)を容易に取得することができる。
【００６９】
また、機械性変動データを用いることにより、まだ到来していない将来のクランク角における機械性回転変動ΔＮc(i+j)を自由に得ることができる。これにより、予測回転速度Ｎ(i+j)には、将来の機械性回転変動と、将来の燃焼性回転変動Ｎb(i)とをそれぞれ反映させることができ、高い予測精度を得ることができる。
【００７０】
また、運転中に計測した実回転変動ΔＮr(i)から機械性回転変動ΔＮc(i)を減算することにより、任意のクランク角における燃焼性回転変動ΔＮb(i)をリアルタイムに算出することができる。これにより、運転中の計測が難しい燃焼性回転変動ΔＮb(i)を、例えばクランク角で１°毎の緻密なデータとして正確かつ容易に得ることができる。
【００７１】
そして、燃焼性回転変動ΔＮb(i)の波形には、内燃機関の燃焼サイクルに同期した周期性が存在するので、現在または過去の時点で算出された燃焼性回転変動（例えば、ΔＮb(i-719)〜ΔＮb(i-540)）を基準として、この時点から１燃焼サイクルが経過した後の燃焼性回転変動ΔＮb(i+1)〜ΔＮb(i+180)を予測することができる。即ち、爆発間隔が短い２つの気筒における燃焼性回転変動の比率（Ｂ1／Ｂ2）が、１燃焼サイクル後でもほぼ一定とみなせることを利用して、将来のクランク角における燃焼性回転変動を正確に予測することができる。
【００７２】
この結果、本実施の形態では、基準回転速度Ｎと、任意のクランク角における機械性回転変動ΔＮc(i)と、将来の燃焼性回転変動ΔＮb(i+1)〜ΔＮb(i+180)とを用いて、将来のクランク角における予測回転速度Ｎ(i+1)〜Ｎ(i+180)を算出することができる。この予測回転速度には、クランク軸３２の回転速度及び将来の回転変動をそれぞれ反映させることができるので、将来の回転速度を高い精度で予測することができる。
【００７３】
従って、本実施の形態によれば、回転速度及び将来の回転変動が反映された予測回転速度Ｎ(i+1)〜Ｎ(i+180)を用いて時間予測を行うことができるので、現在のクランク角から目標角度αまでに必要な目標時間Ｔを正確に予測することができる。これにより、例えば高回転での運転中や、回転速度が大きく変化するような運転状態でも、燃料噴射制御や点火時期制御等からなる各種の制御を所望のタイミングで正確に実行することができる。
【００７４】
［実施の形態１を実現するための具体的な処理］
図７は、本実施の形態のシステム動作を実現するために、ＥＣＵ４０が実行するルーチンのフローチャートである。なお、図７に示すルーチンは、内燃機関の運転中に繰返し実行されるものである。
【００７５】
まず、ステップ１００では、クランク角センサ３４から出力される検出信号（クランク角信号）をカウントすることにより、現在のクランク角を検出する。また、ステップ１０２では、前述したように、クランク角信号を用いて回転速度を検出し、時間的に隣接した２つの回転速度の差分を演算することにより、実回転変動ΔＮr(i)を算出する。
【００７６】
次に、ステップ１０４では、前回の基準回転速度Ｎを用いて図４の無変動角データを参照することにより、当該回転速度における無変動角度を取得する。そして、ステップ１０６では、現在のクランク角が無変動角度と等しいか否かを判定する。ここで「ＹＥＳ」と判定したときには、クランク軸３２の回転速度に変動が生じていないタイミングなので、ステップ１０８では、このタイミングにおいて最新の基準回転速度Ｎを計測する。一方、ステップ１０６で「ＮＯ」と判定したときには、現在のクランク角が無変動角度ではないので、基準回転速度Ｎの更新を行わず、ステップ１１０に移る。
【００７７】
ステップ１１０では、基準回転速度Ｎを用いて図５の機械性変動データを参照することにより、当該回転速度における機械性回転変動ΔＮc(i)を取得する。そして、ステップ１１２では、最新の１燃焼サイクル分のクランク角について、実回転変動ΔＮr(i)と機械性回転変動ΔＮc(i)との差分を求めることにより、燃焼性回転変動ΔＮb(i)を算出し、その算出値を記憶する。
【００７８】
次に、ステップ１１４では、前述したように、過去の燃焼性回転変動ΔＮb(i-719)〜ΔＮb(i-540)を用いて、１燃焼サイクルが経過した将来の燃焼性回転変動ΔＮb(i+1)〜ΔＮb(i+180)を予測し、その予測値を記憶する。
【００７９】
ステップ１１８では、例えば燃料噴射制御、点火時期制御等を行うための時間予測を実施すべきタイミングであるか否かを判定する。ここで「ＹＥＳ」と判定したときには、まず、ステップ１２０において、燃料噴射や点火を行うべき目標クランク角（目標角度）αを読込む。また、ステップ１２２では、前記（４）式を用いることにより、各クランク角における予測回転速度Ｎ(i+j)を算出する。
【００８０】
次に、ステップ１２４では、前記（５）式に示すように、１／Ｎ(i+j)を、現在のクランク角から目標クランク角αまで積算することにより、目標時間Ｔを算出する。一方、ステップ１１８で「ＮＯ」と判定したときには、時間予測が必要ないタイミングなので、今回のルーチン処理を終了する。これにより、適切なタイミングで次回のルーチン処理が起動される。
【００８１】
また、ＥＣＵ４０は、目標時間Ｔを予測した時点で、時間計測用のタイマをスタートさせ、このタイマが目標時間Ｔに一致したときに、燃料噴射や点火等の動作を実施する。これにより、実際のクランク角が目標クランク角αと一致する適切なタイミングにおいて、所定の制御を正確に実施することができる。
【００８２】
なお、前記実施の形態では、図７中のステップ１０８が回転速度算出手段の具体例を示している。また、ステップ１０２は実変動算出手段の具体例を示し、ステップ１１０は機械性変動取得手段の具体例を示している。さらに、ステップ１１２は燃焼性変動算出手段の具体例を示し、ステップ１１４は燃焼性変動予測手段の具体例を示し、ステップ１２２，１２４は時間予測手段の具体例を示している。
【００８３】
実施の形態では、回転変動が生じない無変動角度において、基準回転速度Ｎを計測する構成とした。しかし、本発明では、ベースとなる回転速度を必ずしも無変動角度で計測する必要はない。即ち、本発明では、例えば所定のクランク角で計測した回転速度を、基準回転速度Ｎとして用いる構成としてもよく、または、この回転速度を複数の燃焼サイクル間で平均したものを、基準回転速度Ｎとして用いる構成としてもよい。さらには、実回転変動ΔＮr(i)を算出しつつ、これをモニタすることにより、実回転変動ΔＮr(i)が零となったときの回転速度を、基準回転速度Ｎとして採用する構成としてもよい。
【００８４】
また、実施の形態では、爆発間隔が近い２つの気筒における燃焼性回転変動の比率（Ｂ1／Ｂ2）が、１燃焼サイクル後でもほぼ一定とみなせることを利用して、将来の燃焼性回転変動を容易に予測する構成とした。しかし、本発明はこれ限らず、例えば爆発間隔が近い２つの気筒における燃焼性回転変動の「差分」が、１燃焼サイクル後でもほぼ一定とみなせることを利用してもよい。即ち、図６において、＃１気筒の１燃焼サイクル間における燃焼性回転変動の差分（Ｂ1−Ｂ1′）は、＃２気筒の１燃焼サイクル間における差分（Ｂ2−Ｂ2′）とほぼ等しいとみなすことができる。これにより、将来の燃焼性回転変動Ｂ2′は、Ｂ2′＝Ｂ2−Ｂ1＋Ｂ1′として求めることができる。
【００８５】
また、実施の形態では、燃焼性回転変動ΔＮb(i)の予測演算を行うときに、波形のピーク値Ａn，Ｂnのうちマイナス側のピーク値Ｂnを用いる場合を例示した。しかし、本発明はこれに限らず、プラス側のピーク値Ａnを用いて燃焼性回転変動ΔＮb(i)の予測演算を行う構成としてもよい。
【００８６】
一方、実施の形態では、クランク角センサ３４によりクランク角の１°毎に検出信号を出力する構成とし、回転変動ΔＮr(i)，ΔＮc(i)，ΔＮb(i)と予測回転速度Ｎ(i+j)とを１°毎に算出する構成とした。しかし、本発明はこれに限らず、必要であれば、１°よりも大きなクランク角（例えば３０°等）毎に検出信号を出力するクランク角センサを用いてもよい。同様に、回転変動ΔＮr(i)，ΔＮc(i)，ΔＮb(i)と予測回転速度Ｎ(i+1)についても、１°よりも大きなクランク角毎に算出する構成としてもよい。
【００８７】
実施の形態では、クランク角で１８０°分の範囲において、将来の燃焼性回転変動ΔＮb(i+1)〜ΔＮb(i+180)を予測する構成とした。しかし、本発明は１８０°分の範囲に限らず、必要であればこれよりも広い範囲または狭い範囲において燃焼性回転変動の予測を行う構成としてもよい。
【００８８】
実施の形態では、時間予測制御の適用対象として、燃料噴射制御と点火時期制御を例に挙げて説明した。しかし、本発明の時間予測制御は、これらの制御に限らず、実行タイミングを可変とする任意の制御に広く適用することができる。
【図面の簡単な説明】
【００８９】
【図１】本発明の実施の形態１のシステム構成を説明するための全体図である。
【図２】クランク軸の実回転変動，機械性回転変動及び燃焼性回転変動と、クランク角との関係を示す特性線図である。
【図３】クランク軸の実回転変動とクランク角との関係を示す特性線図である。
【図４】無変動角度と回転速度との関係をマップデータ化した無変動角データを示す説明図である。
【図５】ある回転速度における機械性回転変動とクランク角との関係を示す機械性変動データである。
【図６】燃焼性回転変動とクランク角との関係を示す特性線図である。
【図７】本発明の実施の形態１において実行されるルーチンのフローチャートである。
【符号の説明】
【００９０】
１０  内燃機関
１２  吸気通路
１４  排気通路
１６  エアフローメータ
１８  スロットル弁
２０  スロットルモータ
２２  触媒
２４  燃料噴射弁
２６  点火プラグ
２８  吸気弁
３０  排気弁
３２  クランク軸（出力軸）
３４  クランク角センサ（回転検出手段）
４０  ＥＣＵ（記憶手段）
Ｔ  目標時間
Ｎ  基準回転速度
ΔＮr(i)  実回転変動
ΔＮc(i)  機械性回転変動
ΔＮb(i)  燃焼性回転変動
Ａn，Ｂn  ピーク値
【図１】
2009174483.tif 000003
【図２】
2009174483.tif 000004
【図３】
2009174483.tif 000005
【図４】
2009174483.tif 000006
【図５】
2009174483.tif 000007
【図６】
2009174483.tif 000008
【図７】
2009174483.tif 000009