(57)【要約】【課題】  吸入空気量算出時点以後の可変動弁機構の動作状態に基づいて、吸気弁の閉弁までの時間をより精度よく推定すること。【解決手段】  本発明の内燃機関の制御装置は、変位角変動状態取得手段（変位角取得手段、変位角変化率取得手段）と、閉弁時間推定手段と、を備えている。変位角変動状態取得手段は、算出時点から閉弁時期までの、可変動弁機構における変位角の変動状態を取得（検出あるいは推定）する。閉弁時間推定手段は、推定された変位角の変動状態に基づいて、算出時点から閉弁時期までの時間である閉弁時間を推定するようになっている。【選択図】  図３
2011007052.tif 000002
【特許請求の範囲】
【請求項１】
吸気ポートと気筒との連通状態を制御する吸気弁の作動タイミングを調整する可変動弁機構を備えた内燃機関の制御装置であって、
吸入空気量の算出時点にて、前記可変動弁機構における変位角を取得する、変位角取得手段と、
前記算出時点における前記変位角の変化率を取得する、変位角変化率取得手段と、
前記変位角取得手段及び前記変位角変化率取得手段によって取得された前記変位角及び前記変化率に基づいて、前記算出時点から前記吸気弁の閉弁時期までの時間である閉弁時間を推定する、閉弁時間推定手段と、
を備えたことを特徴とする、内燃機関の制御装置。
【請求項２】
吸気ポートと気筒との連通状態を制御する吸気弁の作動タイミングを調整する可変動弁機構を備えた内燃機関の制御装置であって、
吸入空気量の算出時点から前記吸気弁の閉弁時期までの、前記可変動弁機構における変位角の変動状態を取得する、変位角変動状態取得手段と、
前記変位角変動状態取得手段によって取得された前記変位角の前記変動状態に基づいて、前記算出時点から前記閉弁時期までの時間である閉弁時間を推定する、閉弁時間推定手段と、
を備えたことを特徴とする、内燃機関の制御装置。
【請求項３】
請求項２に記載の、内燃機関の制御装置であって、
前記変位角変動状態取得手段は、
前記変位角を取得する、変位角取得手段と、
前記変位角の変化率を取得する、変位角変化率取得手段と、
を備え、
前記閉弁時間推定手段は、前記変位角取得手段及び前記変位角変化率取得手段によって取得された前記変位角及び前記変化率に基づいて、前記閉弁時間を推定することを特徴とする、内燃機関の制御装置。
【請求項４】
請求項３に記載の、内燃機関の制御装置であって、
前記変位角取得手段は、前記算出時点における前記変位角を取得し、
前記変位角変化率取得手段は、前記算出時点における前記変位角の変化率を取得することを特徴とする、内燃機関の制御装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、吸気ポートと気筒との連通状態を制御する吸気弁の作動タイミングを調整する可変動弁機構を備えた内燃機関の制御装置に関する。
【背景技術】
【０００２】
この種の装置として、特開平８−２８４７３７号公報、特開２００１−９８９９８号公報、等に開示されたものが知られている。これらの装置は、吸入空気量を可能な限り正確に推定することで、燃料噴射量や点火時期を精度よく制御するようになっている。
【発明の概要】
【０００３】
この種の装置において、前記吸気弁が実際に閉弁する時期（閉弁時期）は、前記可変動弁機構の動作状態に応じて変化する。すなわち、吸入空気量の算出時点以後の前記可変動弁機構の動作により、前記閉弁時期が変動し得る。このため、前記算出時点から前記閉弁時期までの時間（閉弁時間）も、可変動弁機構の動作状態に応じて変化し得る。
【０００４】
上述のように、前記算出時点以後前記閉弁時期までの前記可変動弁機構の動作によって前記閉弁時期が変動すると、吸入空気量等の推定精度が悪化し、以て燃料噴射量や点火時期の制御の精度が悪化するおそれがある。したがって、本発明の目的は、前記算出時点以後の前記可変動弁機構の動作状態に基づいて、前記閉弁時間をより精度よく推定することにある。
【０００５】
本発明の内燃機関の制御装置は、変位角変動状態取得手段（変位角取得手段、変位角変化率取得手段）と、閉弁時間推定手段と、を備えている。
【０００６】
前記変位角変動状態取得手段は、前記算出時点から前記閉弁時期までの、前記可変動弁機構における変位角の変動状態を取得（検出あるいは推定）するようになっている。
【０００７】
前記変位角取得手段は、前記可変動弁機構における変位角を取得するようになっている。また、前記変位角変化率取得手段は、前記変位角の変化率を取得するようになっている。
【０００８】
なお、前記変位角取得手段は、前記算出時点における前記変位角を取得するようになっていてもよい。また、前記変位角変化率取得手段は、前記算出時点における前記変位角の変化率を取得するようになっていてもよい。
【０００９】
前記閉弁時間推定手段は、推定された前記変位角の前記変動状態に基づいて、前記算出時点から前記閉弁時期までの時間である閉弁時間を推定するようになっている。具体的には、例えば、前記閉弁時間推定手段は、前記変位角取得手段及び前記変位角変化率取得手段によって取得された前記変位角及び前記変化率に基づいて、前記閉弁時間を推定するようになっている。
【００１０】
かかる構成を有する本発明の内燃機関の制御装置においては、前記変位角の変動状態（前記変位角及びその変化率）に基づいて、前記閉弁時間が推定される。このため、前記可変動弁機構が前記算出時点以後に過渡的に（急激に）動作した場合であっても、かかる動作状態が、前記閉弁時間の推定に対して、より良好に反映される。したがって、本発明によれば、前記閉弁時間が、より精度よく推定され得る。
【図面の簡単な説明】
【００１１】
【図１】本発明の一実施形態が適用される火花点火式多気筒（直列６気筒）の内燃機関の全体構成を示す概略図である。
【図２】バルブタイミングＶＴの変動による閉弁時間Ｔｖの変化の様子を示す概念図である。
【図３】図１に示されているＣＰＵによって実行される、閉弁時間推定処理の一具体例を示すフローチャートである。
【図４】バルブタイミングＶＴの変動による閉弁時間Ｔｖの変化の様子を示す概念図である。
【図５】図１に示されているＣＰＵによって実行される、閉弁時間推定処理の他の具体例を示すフローチャートである。
【発明を実施するための形態】
【００１２】
以下、本発明の実施形態について、図面を参照しつつ説明する。なお、以下の実施形態に関する記載は、法令で要求されている明細書の記載要件（記述要件・実施可能要件）を満たすために、本発明の具体化の単なる一例を、可能な範囲で具体的に記述しているものにすぎない。
【００１３】
よって、後述するように、本発明が、以下に説明する実施形態の具体的構成に何ら限定されるものではないことは、全く当然である。本実施形態に対して施され得る各種の変更（modification）は、当該実施形態の説明中に挿入されると、一貫した実施形態の説明の理解が妨げられるので、末尾にまとめて記載されている。
【００１４】
＜構成＞
図１は、本発明の一実施形態が適用される火花点火式多気筒（直列６気筒）の内燃機関１の全体構成を示す概略図である。内燃機関１は、シリンダブロック部２と、シリンダヘッド部３と、吸気系統４と、排気系統５と、を備えている。
【００１５】
シリンダブロック部２は、シリンダブロック、シリンダブロックロワーケース、及びオイルパン等を含んでいる。シリンダヘッド部３は、シリンダブロック部２の上に固定されている。吸気系統４は、シリンダブロック部２及びシリンダヘッド部３の内部に形成された燃焼室ＣＣに燃料混合気を供給するように設けられている。排気系統５は、燃焼室ＣＣからの排気ガスを浄化しつつ外部に排出するように設けられている。
【００１６】
シリンダブロック部２は、シリンダ２１、ピストン２２、コンロッド２３、及びクランクシャフト２４を備えている。ピストン２２は、シリンダ２１内にて往復動可能に収容されている。クランクシャフト２４は、コンロッド２３を介してピストン２２と連結されていて、ピストン２２の往復動によって回転駆動されるようになっている。燃焼室ＣＣは、シリンダ２１の内側の空間と、ピストン２２の頂面と、シリンダヘッド部３の下端に設けられた凹部と、によって形成されている。
【００１７】
シリンダヘッド部３は、吸気ポート３１、吸気弁３２、可変吸気タイミング装置３３（アクチュエータ３３ａ及び吸気カムシャフト３３ｂを含む）、排気ポート３４、排気弁３５、排気カムシャフト３６、点火プラグ３７、イグナイタ３８、及びポート噴射弁３９を備えている。
【００１８】
吸気ポート３１は、燃焼室ＣＣに連通するように設けられている。吸気弁３２は、吸気カムシャフト３３ｂの回転に応じて、吸気ポート３１とシリンダ２１との連通状態を制御するように設けられている。本発明の可変動弁機構としての可変吸気タイミング装置３３は、吸気弁３２を駆動するための上述の吸気カムシャフト３３ｂの変位角（位相角）をアクチュエータ３３ａによって連続的に変更可能に構成されている。
【００１９】
排気ポート３４は、燃焼室ＣＣに連通するように設けられている。排気弁３５は、排気カムシャフト３６の回転に応じて、排気ポート３４とシリンダ２１との連通状態を制御するように設けられている。イグナイタ３８は、点火プラグ３７に与える高電圧を発生するイグニッションコイルを含んでいる。燃料噴射弁であるポート噴射弁３９は、燃料を吸気ポート３１内に噴射するように設けられている。
【００２０】
吸気系統４は、吸気管４１、エアフィルタ４２、及びスロットル弁４３を備えている。吸気管４１は、吸気ポート３１とともに吸気通路を形成する吸気マニホールドを含んでいて、この吸気マニホールドは各シリンダ２１に対応する吸気ポート３１の各々と連通するように設けられている。エアフィルタ４２は、吸気管４１の端部に設けられている。スロットル弁４３は、ＤＣモータからなるスロットル弁アクチュエータ４３ａによって駆動されることで、吸気管４１内における開口断面積を可変とするように設けられている。
【００２１】
排気系統５は、排気マニホールド５１と、排気管５２と、三元触媒５３と、を備えている。排気マニホールド５１は、排気管５２及び排気ポート３４とともに排気通路を構成するものであって、各シリンダ２１に対応する排気ポート３４の各々と連通するように設けられている。排気管５２は、排気マニホールド５１の集合部に接続されている。この排気管５２には、排気ガスを浄化するための三元触媒５３が介装されている。
【００２２】
また、内燃機関１には、スロットルポジションセンサ６１、カムポジションセンサ６２、クランクポジションセンサ６３、水温センサ６４、熱線式エアフローメータ６５、空燃比センサ６７、及び圧力センサ６６を含む、各種センサ類が設けられている。
【００２３】
スロットルポジションセンサ６１は、スロットル弁４３に対応する位置に設けられている。このスロットルポジションセンサ６１は、スロットル弁４３の開度に対応する信号（この信号からスロットル弁開度TAが取得される）を出力するようになっている。
【００２４】
カムポジションセンサ６２は、シリンダヘッド部３に装着されている。このカムポジションセンサ６２は、吸気カムシャフト３３ｂが９０°回転する毎に（すなわちクランクシャフト２４が１８０°回転する毎に）一つのパルスを有する信号（Ｇ２信号：この信号に基づいて吸気カムシャフト３３ｂの変位角すなわち吸気弁３２の開閉タイミングが取得される）を発生するようになっている。
【００２５】
クランクポジションセンサ６３は、シリンダブロック部２に装着されている。このクランクポジションセンサ６３は、クランクシャフト２４が１０°回転する毎に幅狭のパルスを有するとともに同クランクシャフト２４が３６０°回転する毎に幅広のパルスを有する信号（この信号からエンジン回転速度NEが取得される）を出力するようになっている。
【００２６】
水温センサ６４は、シリンダブロック部２に装着されている。この水温センサ６４は、内燃機関１の内部を通流する冷却水の温度に対応する信号（この信号から冷却水温THWが取得される）を出力するようになっている。
【００２７】
熱線式エアフローメータ６５は、吸気管４１に配設されている。この熱線式エアフローメータ６５は、吸気管４１を流れる吸入空気の単位時間あたりの質量流量に応じた出力電圧（この出力から吸入空気量（流量）Gaが取得される）を出力するようになっている。
【００２８】
圧力センサ６６は、吸気通路に配設されている。この圧力センサ６６は、スロットル弁４３よりも下流であって吸気弁３２よりも上流の吸気管４１内のガスの圧力に対応する信号（この信号から吸気管圧力Pmが取得される）を出力するようになっている。
【００２９】
空燃比センサ６７は、三元触媒５３の上流の排気通路（本例では、排気マニホールド５１の集合部）に配設されている。この空燃比センサ６７は、排気ガス中の酸素ガス濃度に応じた信号（この信号から排気ガスの空燃比が取得される）を出力するようになっている。
【００３０】
＜＜制御装置＞＞
本発明の一実施形態である制御装置７０は、ＣＰＵ７１、ＲＯＭ７２、ＲＡＭ７３、バックアップＲＡＭ７４、及びインターフェース７５を含むマイクロコンピュータである。ＣＰＵ７１、ＲＯＭ７２、ＲＡＭ７３、バックアップＲＡＭ７４、及びインターフェース７５は、双方向バス７６によって互いに接続されている。
【００３１】
本発明の変位角変動状態取得手段（変位角取得手段、変位角変化率取得手段）及び閉弁時間推定手段を構成するＣＰＵ７１は、内燃機関１の各部の動作を制御するためのルーチン（プログラム）、及びこのルーチンを実行するためのテーブル（ルックアップテーブル、マップ）やパラメータを、ＲＯＭ７２、ＲＡＭ７３、及びバックアップＲＡＭ７４から読み出して、ルーチンを実行するようになっている。
【００３２】
ＲＯＭ７２には、ＣＰＵ７１が実行するルーチン、テーブル、定数（パラメータの初期値等）、等が予め格納されている。ＲＡＭ７３は、ＣＰＵ７１がルーチンを実行する際に、必要に応じてデータを一時的に格納し得るように構成されている。バックアップＲＡＭ７４は、電源が投入された状態でＣＰＵ７１がルーチンを実行する際にデータを格納するとともに、この格納したデータを電源遮断後も保持し得るように構成されている。
【００３３】
インターフェース７５は、上述の熱線式エアフローメータ６５、カムポジションセンサ６２、クランクポジションセンサ６３、水温センサ６４、空燃比センサ６７、及び圧力センサ６６、等の各種センサ類と電気的に接続されていて、これらのセンサからの検出信号をＣＰＵ７１に供給するように構成されている。また、インターフェース７５は、ＣＰＵ７１によるルーチン実行によって発生する指令に応じて、可変吸気タイミング装置３３におけるアクチュエータ３３ａ、イグナイタ３８、ポート噴射弁３９、及びスロットル弁アクチュエータ４３ａへ駆動信号を送出するようになっている。
【００３４】
＜実施形態の構成による動作＞
図２は、バルブタイミングＶＴの変動による閉弁時間Ｔｖの変化の様子を示す概念図である。ここで、バルブタイミングＶＴは、吸気カムシャフト３３ｂの変位角、すなわち、吸気弁３２の開閉タイミングに相当する。また、閉弁時間Ｔｖは、筒内吸入空気量Ｍｃの算出時点Ｔ０から、吸気弁３２の閉弁時期までの時間である。
【００３５】
また、図中、二点鎖線で示された直線は、筒内吸入空気量Ｍｃの算出時点Ｔ０のバルブタイミング（変位角）ＶＴと、閉弁時間Ｔｖと、の関係を示すものである（以下、これを「定常特性線」と称する。）。この定常特性線は、バルブタイミングＶＴが一定の場合における閉弁時間Ｔｖを算出するための、以下の式（１）に対応するものである。
Tv＝T120・（310−VT）／120＋TAUa＋tTB＋８    ……（１）
【００３６】
式（１）中、T120は、クランク角を時間換算するために用いられる値であり、具体的にはクランクシャフト２４が１２０°ＣＡ回転するのに要する時間である。ＴＡＵａは、ポート噴射弁３９の通電時間である。この通電時間ＴＡＵａは、燃料噴射量に対応するものであり、筒内吸入空気量Ｍｃとともに時点Ｔ０にて算出される。ｔＴＢは燃料遅れ時間であり、より詳しくはポート噴射弁３９の無効噴射時間、開弁遅れ時間、及び噴霧飛行時間の合計である。なお、式（１）の導出過程については、本出願人に係る特開平８−２８４７３７号公報に記載されているので、本明細書においてはその説明は省略されている。
【００３７】
従来（例えば特開平８−２８４７３７号公報に記載の技術）においては、筒内吸入空気量Ｍｃの算出時点Ｔ０以降のバルブタイミングＶＴが一定である（図中一点鎖線参照）ものとして、閉弁時間Ｔｖが推定されていた（図中一点鎖線（バルブタイミングＶＴ一定を示す線）と定常特性線との交点に対応するＴｖ’参照）。
【００３８】
しかしながら、内燃機関１の運転状態が過渡状態であって、バルブタイミングＶＴが時点Ｔ０以降（時点Ｔ０の前後）において変動している場合、バルブタイミングＶＴが一定とした場合の閉弁時間の推定値Ｔｖ’と、正確な閉弁時間Ｔｖと、の間に、大きな誤差が生じることがある。
【００３９】
そこで、本実施形態の制御装置７０（ＣＰＵ７１）は、時点Ｔ０以降のバルブタイミングＶＴの変動状態を取得（推定）し、この結果に基づいて、閉弁時間Ｔｖを推定する。
【００４０】
＜＜第一の具体例＞＞
具体的には、時点Ｔ０以降のバルブタイミングＶＴの変化率が一定であると推定され得る場合（図２において実線で示されている矢印参照）、以下の式（２）ないし（４）によって、閉弁時間Ｔｖが推定される。
DLVT＝（VT−VT0）／ΔT                              ……（２）
Tv’＝T120・(310−VT)／120                            ……（３）
Tv＝tTFWD・120／(120＋T120・Tv’)＋TAUa＋tTB＋８  ……（４）
【００４１】
ここで、式（２）におけるDLVTは、時点Ｔ０におけるバルブタイミングＶＴの変化率（変位角変化率）である。ΔＴは、演算間隔（８ｍｓｅｃ）である。VT0は、時点Ｔ０よりもΔＴだけ前の（すなわち前回演算時の）バルブタイミングＶＴである。
【００４２】
式（３）におけるTv’は、時点Ｔ０以降のバルブタイミングＶＴが一定である場合を想定した「仮の」閉弁時間であって、上記の式（１）における右辺第１項に対応する。式（４）は、以下のようにして導出される。
【００４３】
図２における長さＹは、図中実線で示されている矢印の傾きDLVTとＴｖとの積によって表される。一方、長さＹは、ＴｖとTv’との偏差と、図中二点鎖線で示された定常特性線の傾きと、の積によっても表される。したがって、下記の式（５）が成立する。
Tv・DLVT＝（Tv’−Tv）・120／T120                  ……（５）
【００４４】
この式（５）を整理すると式（４）の第１項が求められ、これにTAUa＋tTB＋８を加算することで、式（４）が求められる。
【００４５】
図３は、図１に示されているＣＰＵ７１によって実行される、閉弁時間推定処理の一具体例を示すフローチャートである。ＣＰＵ７１は、図３に示されている閉弁時間推定ルーチン３００を、ΔＴ＝８ｍｓｅｃ毎に、繰り返し実行する。
【００４６】
このルーチンの処理が開始されると、まず、ステップ３１０にて、時刻Ｔ０におけるバルブタイミングＶＴが取得される。次に、ステップ３２０にて、時刻Ｔ０におけるバルブタイミング変化率DLVTが、上記の式（２）に基づいて算出される。続いて、ステップ３３０にて、閉弁時間Ｔｖの推定値が、上記の式（３）及び（４）に基づいて算出され、本ルーチンが一旦終了する。
【００４７】
＜＜第二の具体例＞＞
図４は、バルブタイミングＶＴの変動による閉弁時間Ｔｖの変化の様子を示す概念図である。実際の過渡運転時におけるバルブタイミングＶＴは、図４における（ｉ）にて実線の矢印で示されているように、一定（直線状）の変化ではなく、複雑に変化し得る。
【００４８】
そこで、本具体例においては、図４における（ii）に示されているように、筒内吸入空気量Ｍｃ及びポート噴射弁３９の通電時間ＴＡＵａの算出時点の到来毎に、将来のバルブタイミングＶＴが推定される（図中実線の折れ線参照）。そして、このバルブタイミングＶＴの推定結果に基づいて、閉弁時間Ｔｖがより精度よく推定される。
【００４９】
図５は、図１に示されているＣＰＵ７１によって実行される、閉弁時間推定処理の一具体例を示すフローチャートである。ＣＰＵ７１は、図５に示されている閉弁時間推定ルーチン３００を、ΔＴ＝８ｍｓｅｃ毎に、繰り返し実行する。
【００５０】
このルーチンの処理が開始されると、まず、ステップ５１０にて、カウンタｉが１にセットされる。次に、ステップ５２０にて、バルブタイミングVTiが推定される。このバルブタイミングVTiの推定は、例えば、現在のバルブタイミング変化率の推定値とΔＴとの積とから求められる補正値で、現在のバルブタイミングを補正することで行われ得る。あるいは、例えば、バルブタイミングVTiの推定は、スロットル弁開度の予測・推定（特開２００３−１８４６１３号公報：特に図８及びその説明参照）と同様の方法で行われ得る。
【００５１】
続いて、ステップ５３０にて、暫定閉弁時間Tviが、下記の式（６）に基づいて算出される。この式（６）は、上記の式（１）と等価である。
Tvi＝T120・（310−VTi）／120＋TAUa＋tTB＋８    ……（６）
【００５２】
続いて、ステップ５４０にて、暫定閉弁時間Tviがカウンタｉの値とΔＴとの積以下であるか否かが判定される。
【００５３】
ルーチン起動からしばらくの間は、暫定閉弁時間Tviがカウンタｉの値とΔＴとの積よりも大きい（ステップ５４０＝Ｎｏ）。よって、この場合、処理がステップ５５０に移行してカウンタｉがインクリメントされた後、処理がステップ５２０に戻る。
【００５４】
バルブタイミングＶＴの推定結果の軌跡と定常特性線とが交わると、暫定閉弁時間Tviがカウンタｉの値とΔＴとの積以下となる（ステップ５４０＝Ｙｅｓ）。
【００５５】
すなわち、図４における（ii）に示されているように、ｉ＝ｎ−１においては、バルブタイミング推定結果VTn-1から求められる暫定閉弁時間Tvn-1は、ΔＴ・（ｎ−１）よりも大きい。そして、バルブタイミングＶＴの推定結果の軌跡と定常特性線との交点が生じた後のｉ＝ｎにおいては、バルブタイミング推定結果VTnから求められる暫定閉弁時間Tvnは、ΔＴ・ｎ以下となる。
【００５６】
このとき、実際の閉弁時間と同視できる仮想閉弁時間Tvv（バルブタイミングＶＴの推定結果の軌跡と定常特性線との交点に対応する）は、Tvn-1とTvnとの間となる。そこで、ステップ５４０の判定結果がＹｅｓとなった場合、処理がステップ５６０に進行し、Tvn-1及び／又はTvnに基づいて、閉弁時間Ｔｖの推定値が決定され、本ルーチンが一旦終了する。
【００５７】
具体的には、ステップ５６０においては、例えば、閉弁時間Ｔｖの推定値がTvn-1とTvnとのうちのいずれか一方（例えば常にTvn）に決定される。あるいは、例えば、Tvn-1とTvnとを統計処理（平均等）することで、閉弁時間Ｔｖの推定値が決定される。
【００５８】
＜変形例の例示列挙＞
なお、上述の実施形態は、上述した通り、出願人が取り敢えず本願の出願時点において最良であると考えた本発明の代表的な実施形態を単に例示したものにすぎない。よって、本発明はもとより上述の実施形態に何ら限定されるものではない。
【００５９】
したがって、本発明の本質的部分を変更しない範囲内において、上述の実施形態に対して種々の変形が施され得ることは、当然である。
【００６０】
以下、代表的な変形例について、幾つか例示する。もっとも、言うまでもなく、変形例とて、以下に列挙されたものに限定されるものではない。また、複数の変形例が、技術的に矛盾しない範囲内において、適宜、複合的に適用され得る。
【００６１】
本発明（特に、本発明の課題を解決するための手段を構成する各構成要素における、作用的・機能的に表現されているもの）は、上述の実施形態や、下記変形例の記載に基づいて限定解釈されてはならない。このような限定解釈は、（先願主義の下で出願を急ぐ）出願人の利益を不当に害する反面、模倣者を不当に利するものであって、許されない。
【００６２】
（Ａ）本発明は、上述の実施形態にて開示された具体的な装置構成に限定されない。例えば、本発明は、ガソリンエンジン、ディーゼルエンジン、メタノールエンジン、バイオエタノールエンジン、その他任意のタイプの内燃機関に適用可能である。気筒数、気筒配列方式（直列、Ｖ型、水平対向）、燃料供給方式も、特に限定はない。
【００６３】
ポート噴射弁３９とともに、あるいはこれに代えて、燃焼室ＣＣ内に燃料を直接噴射するための筒内噴射弁が設けられていてもよい（例えば特開２００７−２７８１３７号公報等参照）。かかる構成に対しても、本発明は好適に適用される。
【００６４】
（Ｂ）本発明は、上記の実施形態にて開示された具体的な処理に限定されない。例えば、閉弁時間の推定ルーチンは、各シリンダ２１におけるクランク角が吸気上死点前の所定クランク角度（例えば、ＢＴＤＣ９０°ＣＡ）となる毎に、繰り返し実行され得る。
【００６５】
図５におけるステップ５６０にて、以下の式（７）ないし（10）により、閉弁時間Ｔｖの推定値が決定されてもよい。なお、下記の式において、ｎはステップ５４０にて判定がＹｅｓとなった時点のカウンタｉの値である。
α＝T120・（310−VTn-1）／120          ……（７）
DLVT＝（VTn−VTn-1）／ΔT          ……（８）
β＝TAUa＋tTB＋8−(n−1)・ΔT        ……（９）
Tv＝α・120／(120＋T120・DLVT)＋β  ……（10）
【００６６】
これにより、バルブタイミングＶＴが急激に変化した場合でも、閉弁時間Ｔｖの推定が、より精度よく行われ得る。
【００６７】
（Ｃ）その他、特段に言及されていない変形例についても、本発明の本質的部分を変更しない範囲内において、本発明の範囲内に含まれることは当然である。
【００６８】
また、本発明の課題を解決するための手段を構成する各要素における、作用・機能的に表現されている要素は、上述の実施形態や変形例にて開示されている具体的構造の他、当該作用・機能を実現可能ないかなる構造をも含む。
【００６９】
さらに、本明細書にて引用した各公報の内容（明細書及び図面を含む）は、本明細書の一部を構成するものとして援用され得る。
【符号の説明】
【００７０】
１…エンジン                    ＣＣ…燃焼室
２…シリンダブロック部          ２１…シリンダ            ２２…ピストン
３…シリンダヘッド部            ３１…吸気ポート          ３２…吸気弁
３３…可変吸気タイミング装置    ３４…排気ポート          ３５…排気弁
３７…点火プラグ                ３８…イグナイタ          ３９…ポート噴射弁
４…吸気系統                    ４１…吸気管              ４３…スロットル弁
５…排気系統                    ７０…制御装置            ７１…ＣＰＵ
【先行技術文献】
【特許文献】
【００７１】
【特許文献１】特開平８−２８４７３７号公報
【特許文献２】特開２００１−９８９９８号公報
【図１】
2011007052.tif 000003
【図２】
2011007052.tif 000004
【図３】
2011007052.tif 000005
【図４】
2011007052.tif 000006
【図５】
2011007052.tif 000007