(57)【特許請求の範囲】
【請求項１】
気筒内に燃料を直接、噴射するとともに、前記気筒から排気通路に排出された排ガスの一部を外部ＥＧＲガスとして吸気通路に還流させるＥＧＲ装置を備える内燃機関の制御装置であって、
前記気筒内に充填される混合気のＥＧＲ率の目標となる目標ＥＧＲ率を設定する目標ＥＧＲ率設定手段と、
前記気筒内の圧力を筒内圧として検出する筒内圧センサと、
前記内燃機関の運転状態を検出する運転状態検出手段と、
当該検出された内燃機関の運転状態に応じて、混合気の燃焼が開始される直前のクランク角を、基準クランク角として設定する基準クランク角設定手段と、
混合気のＥＧＲ率が前記設定された目標ＥＧＲ率に等しく且つ混合気の空燃比が理論空燃比であるという条件で、混合気の比熱比の温度特性に基づき、前記設定された基準クランク角において発生する気筒内の圧力を、基準筒内圧として算出する基準筒内圧算出手段と、
前記基準クランク角において前記筒内圧センサにより検出された実筒内圧と前記算出された基準筒内圧との圧力差に基づき、混合気のＥＧＲ率を推定するＥＧＲ率推定手段と、
当該推定されたＥＧＲ率に応じて前記内燃機関を制御する制御手段と、
を備えることを特徴とする内燃機関の制御装置。
【請求項２】
前記ＥＧＲ装置は、前記排気通路における過給機のタービンよりも下流側から前記吸気通路における前記過給機のコンプレッサよりも上流側に、外部ＥＧＲガスを還流させるように構成されていることを特徴とする、請求項１に記載の内燃機関の制御装置。
【請求項３】
前記運転状態検出手段は、前記内燃機関の運転状態として、点火時期、前記気筒に吸入される吸気の圧力、及び前記内燃機関の回転数を検出し、
前記基準クランク角設定手段は、前記検出された点火時期、吸気圧力及び内燃機関の回転数に応じて、前記基準クランク角を設定することを特徴とする、請求項１又は２に記載の内燃機関の制御装置。
【請求項４】
圧縮行程において混合気の圧縮が開始される圧縮開始時のクランク角を、初期クランク角として取得する初期クランク角取得手段と、
前記圧縮開始時における前記気筒内の温度を、初期筒内温度として取得する初期筒内温度取得手段と、
前記圧縮開始時における前記気筒内の圧力を、初期筒内圧として取得する初期筒内圧取得手段と、をさらに備え、
前記基準筒内圧算出手段は、前記基準クランク角と、前記取得された初期クランク角、初期筒内温度及び初期筒内圧と、前記設定された目標ＥＧＲ率に応じ、混合気の比熱比の温度特性に基づいて、前記基準筒内圧を算出することを特徴とする、請求項１ないし３のいずれかに記載の内燃機関の制御装置。
【請求項５】
前記内燃機関の回転数を検出する回転数検出手段と、
前記内燃機関を冷却する冷却水の温度を検出する冷却水温度検出手段と、をさらに備え、
前記基準筒内圧算出手段は、前記検出された内燃機関の回転数及び冷却水の温度に応じて、前記基準筒内圧を補正することを特徴とする、請求項４に記載の内燃機関の制御装置。
【請求項６】
前記ＥＧＲ率推定手段は、前記基準クランク角と、前記初期クランク角、前記初期筒内温度及び前記初期筒内圧に応じ、混合気の比熱比の温度特性に基づいて、前記圧力差に対する前記ＥＧＲ率の傾きを表すＥＧＲ係数を算出するとともに、当該算出されたＥＧＲ係数を前記圧力差に乗算した値に、前記目標ＥＧＲ率を加算することによって、前記ＥＧＲ率を算出することを特徴とする、請求項４又は５に記載の内燃機関の制御装置。
【請求項７】
前記筒内圧センサは、前記筒内圧を検出するための圧力検出素子と、当該圧力検出素子から出力される信号を増幅し、出力する増幅回路とを有し、当該圧力検出素子及び増幅回路が、前記気筒内に燃料を噴射する燃料噴射弁に一体に設けられていることを特徴とする、請求項１ないし６のいずれかに記載の内燃機関の制御装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、内燃機関から排出された排ガスの一部を吸気通路に還流させるＥＧＲ装置を備える内燃機関の制御装置に関し、特にＥＧＲ率を推定し、ＥＧＲ率に応じて内燃機関を制御する制御装置に関する。
【背景技術】
【０００２】
従来のＥＧＲ率の推定方法として、気筒に吸入される新気量をエアフローセンサで検出するとともに、吸気圧センサで検出された吸気圧に基づいて気筒に吸入される総ガス量を算出し、これらの新気量と総ガス量からＥＧＲ率を推定することが知られている。しかし、この推定方法を低圧ＥＧＲ装置（ターボチャージャのタービンよりも下流側から排ガスを取り出し、吸気通路のコンプレッサよりも上流側に還流させるＥＧＲ装置）に用いた場合には、気筒までの外部ＥＧＲガスの流路が比較的長く、外部ＥＧＲガスが気筒内に到達するまでの遅れが大きいため、ＥＧＲ率を精度良く推定することは困難である。
【０００３】
また、ＥＧＲ率を推定する従来の他の方法として、例えば特許文献１に開示されたものが知られている。この推定方法は、内燃機関の圧縮行程における混合気の状態変化がポリトロープ変化であることや、混合気の組成に応じて比熱比が変化することに着目したものである。具体的には、圧縮行程中の所定の２つのクランク角ＣＡ１、ＣＡ２において筒内圧Ｐ１、Ｐ２を筒内圧センサで検出し、これらの筒内圧Ｐ１、Ｐ２とクランク角ＣＡ１、ＣＡ２に対応する気筒容積Ｖ１、Ｖ２から、混合気の比熱比κを次式によって算出する。
κ  ＝  ｌｏｇ（Ｐ１／Ｐ２）／ｌｏｇ（Ｖ２／Ｖ１）
そして、算出された比熱比κに基づき、比熱比とＥＧＲ率（ＥＧＲガス濃度）との関係を規定した所定のテーブルを参照することによって、ＥＧＲ率が算出される。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００８−２３１９９５号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
上述した従来のＥＧＲ率の推定方法では、圧縮行程における混合気の状態変化がポリトロープ変化であることを前提とし、圧縮行程中に検出された筒内圧Ｐ１、Ｐ２の変化量に基づいて、混合気の比熱比が算出される。一方、実際の内燃機関の制御では、内燃機関の負荷や回転数などに応じて、点火時期が変更され、それに伴って混合気の燃焼の開始タイミングが変化するのが通常である。
【０００６】
これに対し、従来の推定方法では、筒内圧Ｐ１、Ｐ２の検出タイミングとして、所定のクランク角ＣＡ１、ＣＡ２が一律に設定されているため、後側のクランク角ＣＡ２が混合気の実際の燃焼開始タイミングよりも遅くなる場合がある。その場合には、混合気の燃焼による圧力上昇によって、混合気の状態変化がポリトロープ変化に沿わなくなるため、比熱比の算出精度が低下してしまう。また、このような不具合を回避するために、例えば、後側のクランク角ＣＡ２をより進角側に設定した場合には、筒内圧Ｐ１、Ｐ２の圧力差を十分に確保できないおそれがあり、やはり比熱比を精度良く算出することができない。
【０００７】
また、混合気の比熱比は、温度に応じて変化するという温度特性を有するため、混合気の組成やＥＧＲ率が同じでも、圧縮行程での圧縮による混合気の温度上昇に伴って変化する。これに対し、従来の推定方法では、比熱比を算出する際にも、比熱比からＥＧＲ率を算出する際にも、比熱比の温度特性は考慮されておらず、この点からも比熱比及びＥＧＲ率の算出を精度良く行えない。
【０００８】
本発明は、以上のような課題を解決するためになされたものであり、圧縮行程における混合気の状態変化に基づき、比熱比の温度特性を良好に反映させながら、ＥＧＲ率を精度良く推定できるとともに、推定したＥＧＲ率を用いて内燃機関を適切に制御することができる内燃機関の制御装置を提供することを目的とする。
【課題を解決するための手段】
【０００９】
上記の目的を達成するために、請求項１に係る発明は、気筒３ａ内に燃料を直接、噴射するとともに、気筒３ａから排気通路７に排出された排ガスの一部を外部ＥＧＲガスとして吸気通路６に還流させるＥＧＲ装置１４を備える内燃機関の制御装置であって、気筒３ａ内に充填される混合気のＥＧＲ率Ｒ＿ＥＧＲの目標となる目標ＥＧＲ率ＴＧＴ＿ＥＧＲを設定する目標ＥＧＲ率設定手段（ＥＣＵ２、図２０）と、気筒３ａ内の圧力を筒内圧ＰＣＹＬとして検出する筒内圧センサ５１と、内燃機関３の運転状態（点火時期ＩＧＬＯＧ、吸気圧ＰＢＡ、エンジン回転数ＮＥ）を検出する運転状態検出手段（吸気圧センサ５６、クランク角センサ５２、ＥＣＵ２）と、検出された内燃機関３の運転状態に応じて、気筒３ａ内に充填された混合気の燃焼が開始される直前のクランク角を、基準クランク角ＣＡ＿ＲＥＦとして設定する基準クランク角設定手段（ＥＣＵ２、図３のステップ５、図４）と、混合気のＥＧＲ率Ｒ＿ＥＧＲが設定された目標ＥＧＲ率ＴＧＴ＿ＥＧＲに等しく且つ混合気の空燃比が理論空燃比であるという条件で、混合気の比熱比の温度特性に基づき、設定された基準クランク角ＣＡ＿ＲＥＦにおいて発生する気筒３ａ内の圧力を、基準筒内圧Ｐ＿ＲＥＦとして算出する基準筒内圧算出手段（ＥＣＵ２、図３のステップ６、図６）と、基準クランク角ＣＡ＿ＲＥＦにおいて筒内圧センサ５１により検出された実筒内圧Ｐ＿ＣＰＳと算出された基準筒内圧Ｐ＿ＲＥＦとの圧力差ΔＰに基づき、混合気のＥＧＲ率Ｒ＿ＥＧＲを推定するＥＧＲ率推定手段（ＥＣＵ２、図３のステップ９〜１１）と、推定されたＥＧＲ率Ｒ＿ＥＧＲに応じて内燃機関３を制御する制御手段（ＥＣＵ２、図１８）と、を備えることを特徴とする。
【００１０】
この内燃機関では、気筒内に燃料が直接、噴射されるとともに、気筒から排気通路に排出された排ガスの一部が、外部ＥＧＲガスとして吸気通路に還流する。本発明の内燃機関の制御装置では、混合気のＥＧＲ率の目標となる目標ＥＧＲ率が設定されるとともに、筒内圧（気筒内の圧力）が筒内圧センサによって検出される。また、検出された内燃機関の運転状態に応じて、混合気の燃焼が開始される直前のクランク角が、基準クランク角として設定され、さらに、この基準クランク角において発生する気筒内の圧力が、基準筒内圧として算出される。この基準筒内圧の算出は、ＥＧＲ率が目標ＥＧＲ率に等しく且つ空燃比が理論空燃比であるという混合気の組成の条件で、混合気の比熱比の温度特性に基づいて行われる。
【００１１】
前述したように、混合気の比熱比は、基本的に混合気の組成に応じて定まるとともに、混合気の温度に応じて変化するという温度特性を有する。したがって、混合気のＥＧＲ率が目標ＥＧＲ率に等しく且つ混合気の空燃比が理論空燃比であるという条件で、基準筒内圧の算出を、比熱比の温度特性に基づいて行うことによって、混合気の比熱比の温度特性を反映させながら、基準筒内圧を一義的かつ適切に求めることができる。
【００１２】
また、本発明によれば、基準クランク角において検出された筒内圧を実筒内圧として求め、実筒内圧と基準筒内圧との圧力差に基づいて、ＥＧＲ率が推定される。この実筒内圧には、外部ＥＧＲガスを含む混合気の実際の組成とそれに応じた温度や比熱比が反映される。したがって、実筒内圧と基準筒内圧との圧力差に基づいて、ＥＧＲ率を推定できる。
【００１３】
また、基準クランク角は、混合気の燃焼の開始直前のクランク角であり、検出された内燃機関の運転状態に応じて設定される。このような基準クランク角の設定により、燃焼がまだ行われず、混合気の状態変化がポリトロープ変化に保たれた状態で、実筒内圧を取得するとともに、実筒内圧と基準筒内圧との大きな圧力差を確保できる。したがって、この圧力差に基づき、比熱比の温度特性を良好に反映させながら、ＥＧＲ率を精度良く推定することができる。また、精度良く推定されたＥＧＲ率に応じて内燃機関を制御するので、内燃機関の制御を適切に行うことができる。
【００１４】
請求項２に係る発明は、請求項１に記載の内燃機関の制御装置において、ＥＧＲ装置１４は、排気通路７における過給機（ターボチャージャ１３）のタービン２３よりも下流側から吸気通路６における過給機のコンプレッサ２１よりも上流側に、外部ＥＧＲガスを還流させるように構成されていることを特徴とする。
【００１５】
上記のように構成されるＥＧＲ装置は、いわゆる低圧ＥＧＲ装置であり、外部ＥＧＲガスの流路が比較的長いため、外部ＥＧＲガスが気筒に到達するまでの遅れが大きい傾向がある。上述したように、本発明では、ＥＧＲ率の推定を、気筒内の実際の圧力である実筒内圧とその基準筒内圧をパラメータとして行うので、低圧ＥＧＲ装置の場合にも、外部ＥＧＲガスの遅れの影響を受けることなくＥＧＲ率を精度良く推定でき、したがって、本発明による効果を特に有効に得ることができる。
【００１６】
請求項３に係る発明は、請求項１又は２に記載の内燃機関の制御装置において、運転状態検出手段は、内燃機関３の運転状態として、点火時期ＩＧＬＯＧ、気筒３ａに吸入される吸気の圧力（吸気圧ＰＢＡ）、及び内燃機関３の回転数ＮＥを検出し、基準クランク角設定手段は、検出された点火時期ＩＧＬＯＧ、吸気圧力及び内燃機関３の回転数ＮＥに応じて、基準クランク角ＣＡ＿ＲＥＦを設定すること（図４）を特徴とする。
【００１７】
前述したように、基準クランク角は、混合気の燃焼の開始直前のクランク角として設定される。また、混合気の燃焼の開始タイミングは、点火時期の影響をダイレクトに受けるとともに、吸気圧力に応じて変化し、また、クランク角で表される場合には、内燃機関の回転数に応じて変化する。この構成によれば、検出されたこれらの３つのパラメータに応じて、基準クランク角を設定するので、基準クランク角の設定を、内燃機関の実際の運転状態に応じて適切に行うことができ、基準クランク角における基準筒内圧及び実筒内圧を適切に得ることができる。なお、本明細書における各種のパラメータの「検出」には、パラメータを、センサなどで直接、検出することの他、演算によって推定することなども含まれる。
【００１８】
請求項４に係る発明は、請求項１ないし３のいずれかに記載の内燃機関の制御装置において、圧縮行程において混合気の圧縮が開始される圧縮開始時のクランク角を、初期クランク角（吸気閉弁タイミングＩＶＣ）として取得する初期クランク角取得手段（吸気位相センサ５３、ＥＣＵ２）と、圧縮開始時における気筒３ａ内の温度を、初期筒内温度Ｔ＿ＳＴＲＴとして取得する初期筒内温度取得手段（吸気温センサ５７、吸気位相センサ５３、排気位相センサ５４、ＥＣＵ２、図６のステップ３２）と、圧縮開始時における気筒３ａ内の圧力を、初期筒内圧（吸気圧ＰＢＡ）として取得する初期筒内圧取得手段（吸気圧センサ５６）と、をさらに備え、基準筒内圧算出手段は、基準クランク角ＣＡ＿ＲＥＦと、取得された初期クランク角、初期筒内温度Ｔ＿ＳＴＲＴ及び初期筒内圧と、目標ＥＧＲ率ＴＧＴ＿ＥＧＲに応じ、混合気の比熱比の温度特性に基づいて、基準筒内圧Ｐ＿ＲＥＦを算出すること（図６のステップ３３）を特徴とする。
【００１９】
基準筒内圧は、混合気の燃焼の開始直前に相当する基準クランク角において発生する筒内圧であるため、基準クランク角に応じて変化するとともに、混合気の圧縮の開始タイミングや、圧縮開始時における混合気の温度及び圧力に応じて変化する。また、ＥＧＲ率が変更されると、混合気の組成が変化し、それに応じて比熱比が変化するため、基準筒内圧はＥＧＲ率に応じて変化する。以上の関係から、この構成によれば、基準筒内圧を、基準クランク角と、圧縮開始時の初期クランク角、初期筒内温度及び初期筒内圧と、目標ＥＧＲ率に応じて算出するので、基準筒内圧の算出を適切に行うことができる。
【００２０】
請求項５に係る発明は、請求項４に記載の内燃機関の制御装置において、内燃機関の回転数ＮＥを検出する回転数検出手段（クランク角センサ５２）と、内燃機関３を冷却する冷却水の温度ＴＷを検出する冷却水温度検出手段（水温センサ５９）と、をさらに備え、基準筒内圧算出手段は、検出された内燃機関３の回転数ＮＥ及び冷却水温度ＴＷに応じて基準筒内圧Ｐ＿ＲＥＦを補正すること（図６のステップ３４、３５）を特徴とする。
【００２１】
この構成によれば、検出された内燃機関の回転数及び冷却水温度に応じて基準筒内圧を補正することによって、気筒内と外部との間で授受される熱の影響を補償することができる。
【００２２】
請求項６に係る発明は、請求項４又は５に記載の内燃機関の制御装置において、ＥＧＲ率推定手段は、基準クランク角ＣＡ＿ＲＥＦと、初期クランク角、初期筒内温度Ｔ＿ＳＴＲＴ及び初期筒内圧に応じ、混合気の比熱比の温度特性に基づいて、圧力差ΔＰに対するＥＧＲ率Ｒ＿ＥＧＲの傾きを表すＥＧＲ係数Ｃ＿ＥＧＲを算出するとともに、算出されたＥＧＲ係数Ｃ＿ＥＧＲを圧力差ΔＰに乗算した値に、目標ＥＧＲ率ＴＧＴ＿ＥＧＲを加算することによって、ＥＧＲ率Ｒ＿ＥＧＲを算出すること（図３のステップ７、１１、図１３）を特徴とする。
【００２３】
基準筒内圧は、ＥＧＲ率が目標ＥＧＲ率のときに発生する筒内圧である。このため、実際のＥＧＲ率が目標ＥＧＲ率と等しくなったときには、実筒内圧が基準筒内圧に一致し、両者の圧力差は０になり、実ＥＧＲ率と目標ＥＧＲ率との差であるＥＧＲ率差が大きいほど、圧力差は大きくなる。また、図５に示すように、この圧力差はＥＧＲ率と線形関係（ＥＧＲ率差と比例関係）にあるとともに、その傾き（比例定数）が吸気条件及び圧縮条件に応じて変化するという特性が認められている。
【００２４】
以上の関係から、この構成によれば、ＥＧＲ率を推定する際、まず、圧力差に対するＥＧＲ率の傾きを表すＥＧＲ係数を、基準クランク角、初期クランク角、初期筒内温度及び初期筒内圧に応じて算出する。これにより、基準筒内圧の算出に用いられるパラメータのうち、目標ＥＧＲ率以外のパラメータを用い、混合気の吸気・圧縮条件を反映させながら、ＥＧＲ係数を適切に算出できる。また、算出したＥＧＲ係数を圧力差に乗算した値に目標ＥＧＲ率を加算することによって、ＥＧＲ率を算出するので、ＥＧＲ率を精度良く推定することができる。
【００２５】
請求項７に係る発明は、請求項１ないし６のいずれかに記載の内燃機関の制御装置において、筒内圧センサ５１は、筒内圧を検出するための圧力検出素子と、圧力検出素子から出力される信号を増幅し、出力する増幅回路とを有し、圧力検出素子及び増幅回路が、気筒３ａ内に燃料を噴射する燃料噴射弁４に一体に設けられていることを特徴とする。
【００２６】
上記のように構成される筒内圧センサは、その圧力検出素子及び増幅回路が燃料噴射弁に一体に設けられているため、点火動作によるノイズや他の気筒の燃料噴射弁の噴射動作によるノイズの影響を受けにくい。このため、筒内圧センサによる実筒内圧の検出精度が高められることで、ＥＧＲ率の推定精度をさらに向上させることができる。
【図面の簡単な説明】
【００２７】
【図１】本発明を適用した内燃機関の構成を概略的に示す図である。
【図２】制御装置の概略構成を示すブロック図である。
【図３】ＥＧＲ率の推定処理のメインフローである。
【図４】基準クランク角の設定処理を示すサブルーチンである。
【図５】実筒内圧−基準筒内圧の圧力差とＥＧＲ率及び目標ＥＧＲ率との関係を示す図である。
【図６】基準筒内圧の算出処理を示すサブルーチンである。
【図７】基準筒内圧マップの入出力関係を示す図である。
【図８】混合気の各成分の比熱比の温度特性を示す図である。
【図９】圧縮行程におけるクランク角に対する混合気の比熱比の関係の例を、ＥＧＲガスが存在しない場合と存在する場合について示す図である。
【図１０】基準筒内圧マップにおける、基準クランク角及び吸気閉弁タイミングに対する基準筒内圧の設定例を示す図である。
【図１１】基準筒内圧マップにおける、初期筒内温度に対する基準筒内圧の設定例を示す図である。
【図１２】基準筒内圧マップにおける、吸気圧に対する基準筒内圧の設定例を示す図である。
【図１３】ＥＧＲ係数の算出処理を示すサブルーチンである。
【図１４】ＥＧＲ係数マップの入出力関係を示す図である。
【図１５】ＥＧＲ係数マップにおける、基準クランク角及び吸気閉弁タイミングに対するＥＧＲ係数の設定例を示す図である。
【図１６】ＥＧＲ係数マップにおける、初期筒内温度に対するＥＧＲ係数の設定例を示す図である。
【図１７】ＥＧＲ係数マップにおける、吸気圧に対するＥＧＲ係数の設定例を示す図である。
【図１８】ＥＧＲ率を用いた点火時期の制御処理を示すフローチャートである。
【図１９】変形例による基準筒内圧の算出処理を示すサブルーチンである。
【図２０】目標ＥＧＲ率の設定処理を示すフローチャートである。
【図２１】基準筒内圧マップにおける、目標ＥＧＲ率に対する基準筒内圧の設定例を示す図である。
【発明を実施するための形態】
【００２８】
以下、図面を参照しながら、本発明の好ましい実施形態を詳細に説明する。図１に示すように、本発明を適用した内燃機関（以下「エンジン」という）３は、例えば４つの気筒３ａを有するガソリンエンジンであり、車両（図示せず）に動力源として搭載されている。
【００２９】
エンジン３の各気筒３ａには、燃料噴射弁（以下「インジェクタ」という）４及び点火プラグ５が、気筒３ａの燃焼室（図示せず）に臨むように設けられている。インジェクタ４は、燃焼室内に燃料を直接、噴射するタイプのものである。点火プラグ５からの火花の放電によって、燃料と空気との混合気が点火され、燃焼が行われる。インジェクタ４からの燃料噴射量及び燃料噴射時期と点火プラグ５の点火時期ＩＧＬＯＧは、電子制御ユニット（以下「ＥＣＵ」という）２からの制御信号によって制御される（図２参照）。
【００３０】
なお、本実施形態において「混合気」は、気筒３ａに充填され、燃焼に供される筒内ガスであり、後述するＥＧＲ装置１４による排気還流（外部ＥＧＲ）が行われる場合には、外部ＥＧＲガスを含むものである。
【００３１】
エンジン３の各気筒３ａには、その内部の圧力（筒内圧）を検出する筒内圧センサ５１が設けられている。本実施形態では、筒内圧センサ５１は、インジェクタ一体型のものであり、図示しないが、燃焼室に臨み、筒内圧をピックアップする圧力検出素子や、圧力検出素子からの信号を増幅し、出力する増幅回路などが、インジェクタ４に一体に組み付けられている。筒内圧センサ５１で検出された筒内圧ＰＣＹＬを表す検出信号は、ＥＣＵ２に入力される。
【００３２】
また、エンジン３は、可変吸気位相機構１１、可変排気位相機構１２、ターボチャージャ１３及びＥＧＲ装置１４などを備えている。
【００３３】
可変吸気位相機構１１は、エンジン３のクランクシャフトに対する吸気弁（いずれも図示せず）の相対的な位相（以下「吸気位相」という）ＣＡＩＮを無段階に変更するものであり、吸気位相制御モータ１１ａ（図２参照）などを備えている。吸気位相制御モータ１１ａは、ＥＣＵ２からの制御信号に応じて、クランクシャフトに対して吸気カムシャフト（図示せず）を相対的に回転させ、両者の相対角度を変化させることによって、吸気位相ＣＡＩＮを無段階に変更する。
【００３４】
同様に、可変排気位相機構１２は、クランクシャフトに対する排気弁（図示せず）の相対的な位相（以下「排気位相」という）ＣＡＥＸを無段階に変更するものであり、排気位相制御モータ１２ａ（図２参照）などを備えている。排気位相制御モータ１２ａは、ＥＣＵ２からの制御信号に応じて、クランクシャフトに対して排気カムシャフト（図示せず）を相対的に回転させ、両者の相対角度を変化させることによって、排気位相ＣＡＥＸを無段階に変更する。
【００３５】
これらの可変吸気位相機構１１及び可変排気位相機構１２は、吸気位相ＣＡＩＮと排気位相ＣＡＥＸの変更によって、吸気弁及び排気弁の開閉弁タイミングをそれぞれ制御するとともに、吸気弁と排気弁とのバルブオーバーラップによる内部ＥＧＲを制御するのに用いられる。
【００３６】
ターボチャージャ１３は、吸気通路６に設けられたコンプレッサ２１と、排気通路７に設けられ、シャフト２２を介してコンプレッサ２１に一体に連結されたタービン２３を備えている。排気通路７を流れる排ガスによってタービン２３が駆動され、それと一体にコンプレッサ２１が回転することによって、吸気が過給される。また、ウェイストゲートバルブ（図示せず）などをＥＣＵ２からの制御信号で制御することで、過給圧が調整される。
【００３７】
吸気通路６には、上流側から順に、吸気絞り弁２５、ターボチャージャ１３のコンプレッサ２１、過給によって昇温した吸気を冷却するためのインタークーラ２６、及びスロットル弁２７が設けられている。吸気絞り弁２５は、その下流側に外部ＥＧＲガスを導入するための負圧を発生させるものであり、その開度は、ＥＣＵ２からの制御信号に応じ、ＬＰアクチュエータ２５ａを介して制御される。
【００３８】
スロットル弁２７は、吸気通路６の吸気マニホルド６ａよりも上流側に配置されている。スロットル弁２７の開度は、ＥＣＵ２からの制御信号に応じ、ＴＨアクチュエータ２７ａを介して制御され、それにより、気筒３ａに吸入される筒内ガス量が制御される。
【００３９】
排気通路７のタービン２３よりも下流側には、三元触媒２８が設けられている。三元触媒２８は、活性状態において、排ガス中のＨＣやＣＯを酸化するとともに、ＮＯｘを還元することによって、排ガスを浄化する。
【００４０】
ＥＧＲ装置１４は、気筒３ａから排気通路７に排出された排ガスの一部を、ＥＧＲ通路４１を介し、外部ＥＧＲガスとして吸気通路６に還流させるものである。図１に示すように、ＥＧＲ通路４１は、排気通路７のタービン２３及び三元触媒２８よりも下流側と、吸気通路６のコンプレッサ２１と吸気絞り弁２５との間に接続されている。この構成により、外部ＥＧＲガスは、排ガスがタービン２３に対して仕事を行った後の状態で取り出されるため、比較的低圧になる。すなわち、ＥＧＲ装置１４は、いわゆる低圧ＥＧＲ装置として構成されている。
【００４１】
ＥＧＲ通路４１の途中には、ＥＧＲ弁４２と、外部ＥＧＲガスを冷却するためのＥＧＲクーラ４３が設けられている。ＥＧＲ弁４２の開度は、ＥＣＵ２からの制御信号に応じ、ＥＧＲアクチュエータ４２ａを介して制御され、それにより、外部ＥＧＲガス量が制御される。
【００４２】
また、エンジン３には、その運転状態を検出するために、前述した筒内圧センサ５１に加えて、以下のような各種のセンサが設けられている（図２参照）。
【００４３】
クランク角センサ５２は、クランクシャフトの回転に伴い、所定のクランク角度ごとに、パルス信号であるＣＲＫ信号及びＴＤＣ信号をＥＣＵ２に出力する。ＣＲＫ信号は、所定のクランク角度（例えば０．５度）ごとに出力される。ＥＣＵ２は、このＣＲＫ信号に基づき、エンジン３の回転数（以下「エンジン回転数」という）ＮＥを算出する。
【００４４】
また、ＴＤＣ信号は、いずれかの気筒３ａにおいて、エンジン３のピストン（図示せず）が吸気ＴＤＣ（上死点）付近の所定のクランク角度位置にあることを表す信号であり、本実施形態のようにエンジン３が４気筒の場合には、クランク角度１８０度ごとに出力される。ＥＣＵ２は、ＴＤＣ信号およびＣＲＫ信号に応じて、ＴＤＣ信号の出力タイミングを基準とするクランク角ＣＡを、気筒３ａごとに算出する。また、ＥＣＵ２は、ＴＤＣ信号及びＣＲＫ信号に応じて、所定のクランク角度（例えば３０度）ごとに、クランク角ステージＦＩＳＴＧ（＝０〜２３）を算出し、割り当てる。
【００４５】
また、可変吸気位相機構１１を取り付けた吸気カムシャフト、及び可変排気位相機構１２を取り付けた排気カムシャフトには、吸気位相センサ５３及び排気位相センサ５４がそれぞれ設けられている。吸気位相センサ５３は、吸気カムシャフトの回転に伴い、所定のカム角度（例えば０．５度）ごとに、パルス信号であるＣＡＭＩＮ信号をＥＣＵ２に出力する。ＥＣＵ２は、このＣＡＭＩＮ信号とＣＲＫ信号に基づき、吸気位相ＣＡＩＮを算出する。同様に、排気位相センサ５４は、排気カムシャフトの回転に伴い、所定のカム角度（例えば０．５度）ごとに、ＣＡＭＥＸ信号をＥＣＵ２に出力する。ＥＣＵ２は、このＣＡＭＥＸ信号とＣＲＫ信号に基づき、排気位相ＣＡＥＸを算出する。
【００４６】
また、吸気通路６には、吸気絞り弁２５の上流側にエアフローセンサ５５が設けられ、スロットル弁２７の下流側の吸気チャンバ６ｂに、吸気圧センサ５６及び吸気温センサ５７が設けられている。エアフローセンサ５５は、気筒３ａに吸入される空気（新気）の量（吸入空気量）ＧＡＩＲを検出し、吸気圧センサ５６は、吸気圧ＰＢＡを絶対圧として検出し、吸気温センサ５７は、気筒３ａに吸入される、外部ＥＧＲガスを含む吸気の温度（吸気温）ＴＡを検出する。これらの検出信号はＥＣＵ２に入力される。
【００４７】
排気通路７のタービン２３と三元触媒２８との間には、ＬＡＦセンサ５８が設けられている。ＬＡＦセンサ５８は、理論空燃比を含む広い空燃比領域において、三元触媒２８に流入する排ガス中の酸素濃度を連続的に検出し、その検出信号をＥＣＵ２に出力する。ＥＣＵ２は、この検出信号に基づき、排ガスの当量比ＫＡＣＴを算出する。
【００４８】
さらに、ＥＣＵ２には、水温センサ５９からエンジン３を冷却する冷却水の温度（以下「エンジン水温」という）ＴＷを表す検出信号が、アクセル開度センサ６０から、車両のアクセルペダル（図示せず）の踏込み量（以下「アクセル開度」という）ＡＰを表す検出信号が、それぞれ入力される。
【００４９】
ＥＣＵ２は、ＣＰＵ、ＲＡＭ、ＲＯＭ及びＩ／Ｏインターフェース（いずれも図示せず）などから成るマイクロコンピュータで構成されている。ＥＣＵ２は、前述した各種のセンサの検出信号などに応じて、エンジン３の運転状態を判別し、インジェクタ４の燃料噴射量や点火プラグ５の点火時期ＩＧＬＯＧの制御などを含むエンジン制御を実行する。
【００５０】
また、本実施形態では特に、ＥＣＵ２は、気筒３ａ内に充填される混合気のＥＧＲ率Ｒ＿ＥＧＲを推定するとともに、推定したＥＧＲ率Ｒ＿ＥＧＲに応じて点火時期ＩＧＬＯＧを制御する。なお、混合気のＥＧＲ率Ｒ＿ＥＧＲは、混合気（筒内ガス）の総量に対するＥＧＲガス量の比率として定義される。
【００５１】
本実施形態では、ＥＣＵ２が、目標ＥＧＲ率設定手段、基準クランク角設定手段、基準筒内圧算出手段、ＥＧＲ率推定手段、制御手段、初期クランク角取得手段、及び初期筒内温度取得手段に相当する。
【００５２】
図２０は、ＥＣＵ２で実行される、目標ＥＧＲ率ＴＧＴ＿ＥＧＲの設定処理を示す。本処理では、まずステップ６１において、アクセル開度ＡＰ及びエンジン回転数ＮＥに応じ、所定のマップ（図示せず）を検索することによって、要求トルクＴＲＱＣＭＤを算出する。次に、要求トルクＴＲＱＣＭＤ及びエンジン回転数ＮＥに応じ、所定のマップ（図示せず）を検索することによって、目標ＥＧＲ率ＴＧＴ＿ＥＧＲを算出し（ステップ６２）、本処理を終了する。目標ＥＧＲ率ＴＧＴ＿ＥＧＲが設定されると、それに応じた制御信号がＥＧＲアクチュエータ４２ａに出力され、ＥＧＲ弁４２の開度が制御されることで、ＥＧＲ率Ｒ＿ＥＧＲが目標ＥＧＲ率ＴＧＴ＿ＥＧＲに制御される。
【００５３】
なお、上述した目標ＥＧＲ率ＴＧＴ＿ＥＧＲの設定処理は一例であり、他の適当な設定方法を採用することが可能である。例えば、上記の要求トルクＴＲＱＣＭＤに代えて、エンジン３の負荷を表す他のパラメータ、例えば吸入空気量ＧＡＩＲや燃料噴射量などを用いて、目標ＥＧＲ率ＴＧＴ＿ＥＧＲを設定してもよく、さらにノッキングの発生状況などに応じた補正を加えてもよい。
【００５４】
次に、ＥＣＵ２で実行されるＥＧＲ率Ｒ＿ＥＧＲの推定処理について説明する。図３は、そのメインフローを示す。本処理は、気筒３ａごとに、前述したクランク角ステージＦＩＳＴＧの切替周期と同じ周期（例えばクランク角度３０度ごと）で、繰り返し実行される。なお、筒内圧センサ５１で検出された筒内圧ＰＣＹＬに直接、関連する処理は、本処理とは別個に、ＣＲＫ信号の発生周期と同じ周期（例えばクランク角度０．５度ごと）で実行され、例えば、検出された筒内圧ＰＣＹＬがクランク角ＣＡに対応して記憶される。
【００５５】
図３の推定処理では、まずステップ１（「Ｓ１」と図示。以下同じ）において、クランク角ステージＦＩＳＴＧが、吸気ＴＤＣ（上死点）に相当する第１所定値ＳＴＧ１に等しいか否かを判別する。この判別結果がＹＥＳで、当該気筒３ａが吸気行程に移行した直後の段階にあるときには、吸気関連パラメータを取得する（ステップ２）。具体的には、吸気関連パラメータとして、吸気温ＴＡ、エンジン水温ＴＷ及び排気位相ＣＡＥＸを読み出すとともに、ＥＣＵ２のＲＡＭの所定領域に記憶する。その後、本処理を終了する。
【００５６】
前記ステップ１の判別結果がＮＯのときには、クランク角ステージＦＩＳＴＧが、圧縮ＢＤＣ（下死点）に相当する第２所定値ＳＴＧ２に等しいか否かを判別する（ステップ３）。この判別結果がＹＥＳで、当該気筒３ａが圧縮行程に移行した直後の段階にあるときには、圧縮関連パラメータを取得する（ステップ４）。具体的には、圧縮関連パラメータとして、検出された吸気圧ＰＢＡ、エンジン回転数ＮＥ及び吸気位相ＣＡＩＮと、その時点で設定されている点火時期ＩＧＬＯＧを読み出すとともに、ＥＣＵ２のＲＡＭの所定領域に記憶する。
【００５７】
次に、基準クランク角ＣＡ＿ＲＥＦの設定処理を実行する（ステップ５）。この設定処理は、混合気の燃焼が開始される直前のタイミングを予測し、基準クランク角ＣＡ＿ＲＥＦとして設定するものである。図４はそのサブルーチンを示す。
【００５８】
本処理では、まずステップ２１において、前記ステップ４で取得した吸気圧ＰＢＡ及びエンジン回転数ＮＥに応じ、所定のマップ（図示せず）を検索することによって、遅角補正量ΔＣ＿ＣＡを算出する。この遅角補正量ΔＣ＿ＣＡは、点火時期ＩＧＬＯＧで点火プラグ５による点火動作が行われた後、混合気が着火し、燃焼が開始されるまでの着火遅れ時間に相当し、クランク角度で表される。吸気圧ＰＢＡが低いほど、混合気が着火しにくくなり、また、エンジン回転数ＮＥが高いほど、同じ着火遅れ時間に対応するクランク角度は大きくなる。このため、上記のマップでは、遅角補正量ΔＣ＿ＣＡは、吸気圧ＰＢＡが低いほど、また、エンジン回転数ＮＥが高いほど、より大きくなるように設定されている。
【００５９】
次に、前記ステップ４で取得した点火時期ＩＧＬＯＧから遅角補正量ΔＣ＿ＣＡを減算した値を、基準クランク角ＣＡ＿ＲＥＦとして設定し（ステップ２２）、本処理を終了する。なお、この基準クランク角ＣＡ＿ＲＥＦは、各気筒３ａの圧縮ＴＤＣを原点（０度）とし、進角側を正として表される（図１０参照）。
【００６０】
図３に戻り、上記ステップ５に続くステップ６では、基準筒内圧Ｐ＿ＲＥＦの算出処理を実行する。この基準筒内圧Ｐ＿ＲＥＦは、混合気のＥＧＲ率Ｒ＿ＥＧＲが図２０の処理で設定された目標ＥＧＲ率ＴＧＴ＿ＥＧＲに等しく、且つ混合気の空燃比が理論空燃比であるという条件で、上記の基準クランク角において発生する筒内圧である。その算出処理の詳細については、後述する。
【００６１】
次に、ＥＧＲ係数Ｃ＿ＥＧＲの算出処理を実行し（ステップ７）、本処理を終了する。図５に示すように、このＥＧＲ係数Ｃ＿ＥＧＲは、圧力差ΔＰ（後述する実筒内圧Ｐ＿ＣＰＳと基準筒内圧Ｐ＿ＲＥＦとの差）とＥＧＲ率差ΔＥＧＲ（実ＥＧＲ率Ｒ＿ＥＧＲと目標ＥＧＲ率ＴＧＴ＿ＥＧＲとの差）の間に、比例関係が認められることから、その傾き（＝ΔＥＧＲ／ΔＰ）をＥＧＲ係数Ｃ＿ＥＧＲと定義したものである。その算出処理の詳細については、後述する。
【００６２】
前記ステップ３の判別結果がＮＯのときには、クランク角ステージＦＩＳＴＧが、圧縮ＴＤＣ（上死点）に相当する第３所定値ＳＴＧ３に等しいか否かを判別する（ステップ８）。この判別結果がＮＯのときには、そのまま本処理を終了する。一方、ステップ８の判別結果がＹＥＳで、当該気筒３ａが圧縮行程が終了した直後の段階にあるときには、ステップ５で設定した基準クランク角ＣＡ＿ＲＥＦにおいて検出された筒内圧ＰＣＹＬを、ＲＡＭから読み出し、実筒内圧Ｐ＿ＣＰＳとして取得する（ステップ９）。
【００６３】
次に、取得した実筒内圧Ｐ＿ＣＰＳとステップ６で算出された基準筒内圧Ｐ＿ＲＥＦとの差（＝Ｐ＿ＣＰＳ−Ｐ＿ＲＥＦ）を、圧力差ΔＰとして算出する（ステップ１０）。次に、これまでに算出した圧力差ΔＰ及びＥＧＲ係数Ｃ＿ＥＧＲと目標ＥＧＲ率ＴＧＴ＿ＥＧＲを用い、次式（Ａ）によって、混合気のＥＧＲ率Ｒ＿ＥＧＲを算出し（ステップ１１）、本処理を終了する。
Ｒ＿ＥＧＲ  ＝  ΔＰ×Ｃ＿ＥＧＲ＋ＴＧＴ＿ＥＧＲ      ・・・（Ａ）
なお、この式（Ａ）は、上述したＥＧＲ係数Ｃ＿ＥＧＲの定義を表す次式（Ａ）’から直接、導き出される（図５参照）。
Ｃ＿ＥＧＲ  ＝  ΔＥＧＲ／ΔＰ  
＝（Ｒ＿ＥＧＲ−ＴＧＴ＿ＥＧＲ）／ΔＰ    ・・・（Ａ）’
【００６４】
次に、図６を参照しながら、図３のステップ６で実行される基準筒内圧Ｐ＿ＲＥＦの算出処理について説明する。本処理では、まずステップ３１において、前記ステップ２で取得した吸気位相ＣＡＩＮから、吸気弁の閉弁タイミング（以下「吸気閉弁タイミング」という）ＩＶＣを算出する。この吸気閉弁タイミングＩＶＣは、前述した基準クランク角ＣＡ＿ＲＥＦと同様、圧縮ＴＤＣを原点（０度）とし、進角側を正とするクランク角で表される。
【００６５】
この吸気閉弁タイミングＩＶＣが圧縮行程中に設定される場合、混合気の圧縮は実質的に吸気弁の閉弁時から開始されるので、吸気閉弁タイミングＩＶＣは、圧縮開始時のクランク角（初期クランク角）に相当する。また、吸気圧ＰＢＡは、圧縮開始時における筒内圧（初期筒内圧）に相当する。
【００６６】
次に、吸気温ＴＡ、吸気位相ＣＡＩＮ及び排気位相ＣＡＥＸに応じ、所定のマップ（図示せず）を検索することによって、圧縮開始時における気筒３ａ内の温度である初期筒内温度Ｔ＿ＳＴＲＴを算出する（ステップ３２）。上記のパラメータのうち、吸気位相ＣＡＩＮ及び排気位相ＣＡＥＸは、吸気弁と排気弁とのバルブオーバーラップによる内部ＥＧＲが実行される場合に、内部ＥＧＲ量に応じた筒内温度の上昇を反映させるためのものである。このため、上記のマップでは、初期筒内温度Ｔ＿ＳＴＲＴは、吸気温ＴＡが高いほど、また、吸気位相ＣＡＩＮ及び排気位相ＣＡＥＸに対しては、バルブオーバーラップが大きい側に位置するほど、より高い値に設定されている。
【００６７】
次のステップ３３では、基準クランク角ＣＡ＿ＲＥＦ、吸気閉弁タイミングＩＶＣ、初期筒内温度Ｔ＿ＳＴＲＴ及び吸気圧ＰＢＡと、目標ＥＧＲ率ＴＧＴ＿ＥＧＲに応じ、図７に示す基準筒内圧マップを検索することによって、基準筒内圧Ｐ＿ＲＥＦを算出する。以下、この基準筒内圧マップについて説明する。
【００６８】
まず、気筒３ａ内に充填された混合気（筒内ガス）の比熱比と圧縮行程における状態変化について説明する。混合気の比熱比κは、定圧比熱Ｃｐと気体定数Ｒを用いて次式（１）で表され、定圧比熱Ｃｐは次式（２）で表される。
【数１】
0006456256.tif 000002
【数２】
0006456256.tif 000003
【００６９】
式（２）に示されるように、混合気の比熱比κは、その組成（成分と各成分のモル数）に応じて変化する。また、図８に示すように、混合気の各成分の比熱比は、温度が上昇するにつれて低下するという温度特性を有し、これらの成分で構成される混合気の比熱比κもまた、同様の温度特性を有する。さらに、図９に示すように、混合気中にＥＧＲガスが含まれる場合には、混合気の組成が変化し、ＥＧＲガスのＣＯ２成分が加わることで、混合気の比熱比κが増大するという特性を有する。
【００７０】
また、圧縮行程における混合気の状態変化が、断熱圧縮変化であり、ポリトロープ変化とみなされることから、クランク角ＣＡ＝ａのときの筒内温度Ｔａは、次式（３）で表される。
【数３】
0006456256.tif 000004
【００７１】
式（３）に示されるように、筒内温度Ｔは比熱比κの関数であり、また、上述したように、混合気の比熱比κは筒内温度Ｔの関数である。このため、比熱比κ及び筒内温度Ｔを正確に求めるために、式（１）（２）と式（３）の演算結果を逐次、相互に適用した逐次計算が行われる。その結果、クランク角ＣＡ＝最終クランク角θのときの筒内温度（最終筒内温度）Ｔθは、次式（４）で表される。
【数４】
0006456256.tif 000005
【００７２】
また、クランク角ＣＡ＝ａのときの筒内圧Ｐａは、次式（５）で表され、さらに、この式（５）から、クランク角ＣＡ＝θのときの筒内圧（最終筒内圧）Ｐθは、次式（６）で表される。
【数５】
0006456256.tif 000006
【数６】
0006456256.tif 000007
【００７３】
式（６）に示されるように、最終筒内圧Ｐθは、初期筒内圧Ｐ０、初期気筒容積Ｖ０、最終気筒容積Ｖθ及び逐次計算される比熱比κの関数である。また、比熱比κは、逐次計算される筒内温度Ｔの関数であり、筒内温度Ｔは、初期筒内温度Ｔ０と比熱比κの関数である。また、気筒容積Ｖはクランク角ＣＡから一義的に求められるため、初期気筒容積Ｖ０と最終気筒容積Ｖθは、初期クランク角ＣＡ０と最終クランク角ＣＡθにそれぞれ置き換えられる。
【００７４】
以上から、最終筒内圧Ｐθは、式（２）における混合気の組成が与えられた条件で、初期筒内圧Ｐ０、初期筒内温度Ｔ０、初期クランク角ＣＡ０及び最終クランク角ＣＡθの関数として、求められる。
【００７５】
前述した基準筒内圧マップは、以上の関係に基づいており、図７に示すように、初期筒内圧Ｐ０、初期筒内温度Ｔ０及び初期クランク角ＣＡ０にそれぞれ相当する吸気圧ＰＢＡ、初期筒内温度Ｔ＿ＳＴＲＴ及び吸気閉弁タイミングＩＶＣと、最終クランク角ＣＡθに相当する基準クランク角ＣＡ＿ＲＥＦを入力パラメータとし、最終筒内圧Ｐθに相当する基準筒内圧Ｐ＿ＲＥＦを出力として得るものである。
【００７６】
また、混合気の組成の条件として、その時点で設定されている目標ＥＧＲ率ＴＧＴ＿ＥＧＲが入力パラメータに加えられ、この目標ＥＧＲ率ＴＧＴ＿ＥＧＲと、空燃比が理論空燃比であるという条件から、式（２）の各成分のモル数ｎｘが定まる。具体的には、上記の５つの入力パラメータの様々な条件に対し、式（１）〜（６）に基づいて基準筒内圧Ｐ＿ＲＥＦをあらかじめ算出し、その算出結果を入力パラメータに対してマップ化することによって、基準筒内圧マップが作成されている。
【００７７】
図１０〜図１２及び図２１は、基準筒内圧マップにおける、各入力パラメータに対する基準筒内圧Ｐ＿ＲＥＦの設定例を示す。まず、図２１に示すように、基準筒内圧Ｐ＿ＲＥＦは、目標ＥＧＲ率Ｃ＿ＥＧＲが高いほど、より大きな値に設定されている。これは、ＥＧＲ率が高いほど、混合気の比熱比κが高くなり（図９参照）、それに伴って最終的な筒内圧が大きくなるためである。
【００７８】
また、図１０に示すように、基準筒内圧Ｐ＿ＲＥＦは、基準クランク角ＣＡ＿ＲＥＦの値が０に近いほど、すなわち基準クランク角ＣＡ＿ＲＥＦが圧縮ＴＤＣに近いほど、より大きな値に設定され、また、吸気閉弁タイミングＩＶＣの値が大きいほど、すなわち圧縮行程における吸気弁の閉弁タイミングが早いほど、より大きな値に設定されている。これは、基準クランク角ＣＡ＿ＲＥＦが圧縮ＴＤＣに近いほど、また、吸気弁の閉弁タイミングが早いほど、混合気の実質的な圧縮期間が長くなることで、最終的な筒内圧が大きくなるためである。
【００７９】
また、図１１に示すように、基準筒内圧Ｐ＿ＲＥＦは、初期筒内温度Ｔ＿ＳＴＲＴが高いほど、より小さな値に設定されている。これは、初期筒内温度Ｔ＿ＳＴＲＴが高いほど、筒内温度がより高くなるのに応じて混合気の比熱比κが低下する結果、筒内圧の上昇度合が低下するためである。
【００８０】
さらに、図１２に示すように、基準筒内圧Ｐ＿ＲＥＦは、吸気圧ＰＢＡに比例するように設定されている。これは、基準筒内圧Ｐ＿ＲＥＦ及び吸気圧ＰＢＡは、最終筒内圧Ｐθ及び初期筒内圧Ｐ０にそれぞれ相当し、両者が比例関係にあるためである（式（６）参照）。
【００８１】
前述したように、図６のステップ３３では、上記の５つのパラメータに応じ、上述した基準筒内圧マップを検索することによって、基準筒内圧Ｐ＿ＲＥＦが算出される。次のステップ３４では、エンジン回転数ＮＥ及びエンジン水温ＴＷに応じ、所定のマップを検索することによって、伝熱補正係数Ｋ＿ＨＴを算出する。この伝熱補正係数Ｋ＿ＨＴは、気筒３ａ内と外部との間で授受される熱の影響を補償するためのものである。
【００８２】
次に、ステップ３３で算出された基準筒内圧Ｐ＿ＲＥＦに伝熱補正係数Ｋ＿ＨＴを乗算することによって、最終的な基準筒内圧Ｐ＿ＲＥＦを算出し（ステップ３５）、本処理を終了する。
【００８３】
次に、図１３を参照しながら、図３のステップ７で実行されるＥＧＲ係数Ｃ＿ＥＧＲの算出処理について説明する。前述したように、ＥＧＲ係数Ｃ＿ＥＧＲは、圧力差ΔＰ（実筒内圧Ｐ＿ＣＰＳと基準筒内圧Ｐ＿ＲＥＦとの差）に対するＥＧＲ率差ΔＥＧＲ（実ＥＧＲ率Ｒ＿ＥＧＲと目標ＥＧＲ率ＴＧＴ＿ＥＧＲとの差）の傾きとして定義され（図５参照）、ＥＧＲ率Ｒ＿ＥＧＲの算出に用いられる。また、上記の傾きが吸気条件及び圧縮条件に応じて変化するという特性が認められたため、本処理において、ＥＧＲ係数Ｃ＿ＥＧＲを算出するものである。
【００８４】
本処理では、まずステップ４１において、基準クランク角ＣＡ＿ＲＥＦ、吸気閉弁タイミングＩＶＣ、初期筒内温度Ｔ＿ＳＴＲＴ及び吸気圧ＰＢＡを取得する。これらのパラメータは、上記の吸気条件及び圧縮条件を表すものであり、前述した基準筒内圧マップの４つの入力パラメータと共通である。このため、ステップ４１におけるパラメータの取得は、図６の基準筒内圧Ｐ＿ＲＥＦの算出処理で得られたデータを読み出すことによって、行われる。
【００８５】
次に、取得した４つのパラメータに応じ、図１４に示すＥＧＲ係数マップを検索することによって、ＥＧＲ係数Ｃ＿ＥＧＲを算出し（ステップ４２）、本処理を終了する。このＥＧＲ係数マップは、上記の４つの入力パラメータの様々な条件に対し、式（１）〜（６）に基づいてＥＧＲ係数Ｃ＿ＥＧＲをあらかじめ算出し、その結果を入力パラメータに対してマップ化したものである。なお、基準筒内圧Ｐ＿ＲＥＦと異なり、ＥＧＲ係数Ｃ＿ＥＧＲは目標ＥＧＲ率ＴＧＴ＿ＥＧＲの影響を受けないので、目標ＥＧＲ率ＴＧＴ＿ＥＧＲは、ＥＧＲ係数マップの入力パラメータには含まれない。
【００８６】
図１５〜図１７は、ＥＧＲ係数マップにおける、各入力パラメータに対するＥＧＲ係数Ｃ＿ＥＧＲの設定例を示す。図１５に示すように、ＥＧＲ係数Ｃ＿ＥＧＲは、基準クランク角ＣＡ＿ＲＥＦが圧縮ＴＤＣに近いほど、また、圧縮行程における吸気弁の閉弁タイミングが早いほど、より小さな値に設定されている。これは、基準クランク角ＣＡ＿ＲＥＦが圧縮ＴＤＣに近いほど、また、吸気弁の閉弁タイミングが早いほど、混合気の実質的な圧縮期間が長くなることで、圧力差ΔＰが大きくなり、それに応じてＥＧＲ係数Ｃ＿ＥＧＲが小さくなるためである。
【００８７】
また、図１６に示すように、ＥＧＲ係数Ｃ＿ＥＧＲは、初期筒内温度Ｔ＿ＳＴＲＴが高いほど、より小さな値に設定されている。これは、以下の理由による。すなわち、混合気の成分の中で、燃料は、他の成分と比較して、定圧比熱Ｃｐの温度変化が大きく、混合気の比熱比κの温度特性への寄与度が大きい。一方、ＥＧＲ率Ｒ＿ＥＧＲが増加すると、それに伴って燃料の割合が低下し、その寄与度が低くなることで、温度に応じた比熱比κの変化度合は小さくなる。このため、初期筒内温度Ｔ＿ＳＴＲＴが高いほど、圧縮中に比熱比κがより高いレベルで変化することで、圧力差ΔＰが大きくなり、ＥＧＲ係数Ｃ＿ＥＧＲが小さくなるためである。
【００８８】
さらに、図１７に示すように、ＥＧＲ係数Ｃ＿ＥＧＲは、吸気圧ＰＢＡが高いほど、より小さな値に設定されている。これは、初期筒内圧である吸気圧ＰＢＡが高いほど、それに比例して実筒内圧Ｐ＿ＣＰＳ及び圧力差ΔＰが増大し、それに応じてＥＧＲ係数Ｃ＿ＥＧＲが小さくなるためである。
【００８９】
ここで、図５を再び参照しながら、上述したＥＧＲ率Ｒ＿ＥＧＲの推定処理によって得られる動作について説明する。例えば、ＥＧＲが停止された状態（目標ＥＧＲ率ＴＧＴ＿ＥＧＲ＝０）から、図５に示す目標ＥＧＲ率ＴＧＴ＿ＥＧＲが設定され、外部ＥＧＲが開始された場合、その開始時から低圧ＥＧＲ装置に顕著な遅れ時間が経過するまでは、実際のＥＧＲ率Ｒ＿ＥＧＲが０であることで、ＥＧＲ率差ΔＥＧＲは負値であり、また、実筒内圧Ｐ＿ＣＰＳが基準筒内圧ＰＲＥＦよりも小さいことで、圧力差ΔＰも負値である。
【００９０】
上記の遅れ時間が経過し、気筒３ａに外部ＥＧＲガスが導入されるようになると、ＥＧＲ率Ｒ＿ＥＧＲは０から増加し、最終的に目標ＥＧＲ率ＴＧＴ＿ＥＧＲに達する。これに伴って実筒内圧Ｐ＿ＣＰＳが増加することで、圧力差ΔＰは、同図に矢印Ａで示すように、ＥＧＲ係数Ｃ＿ＥＧＲに対応する傾きで線形に増加し、ＥＧＲ率Ｒ＿ＥＧＲが目標ＥＧＲ率ＴＧＴ＿ＥＧＲに達したときに、０になる。したがって、この間の任意の点において圧力差ΔＰ＝ΔＰ１が得られた場合、ＥＧＲ係数Ｃ＿ＥＧＲと目標ＥＧＲ率ＴＧＴ＿ＥＧＲを用い、前記式（Ａ）によって、そのときの実ＥＧＲ率Ｒ＿ＥＧＲ１を精度良く算出することができる。
【００９１】
上記とは逆に、目標ＥＧＲ率ＴＧＴ＿ＥＧＲがより大きな値から小さな値に変更された場合にも、上記と同様にＥＧＲ率Ｒ＿ＥＧＲを算出できる。すなわち、この場合には、目標ＥＧＲ率ＴＧＴ＿ＥＧＲの変更開始時から遅れ時間が経過するまでは、実際のＥＧＲ率Ｒ＿ＥＧＲが変更前の目標ＥＧＲ率のままであることで、ＥＧＲ率差ΔＥＧＲは正値であり、それに応じて圧力差ΔＰもまた正値である。
【００９２】
上記の遅れ時間が経過し、気筒３ａに導入される外部ＥＧＲガス量が減少するようになると、ＥＧＲ率Ｒ＿ＥＧＲが減少し、最終的に目標ＥＧＲ率ＴＧＴ＿ＥＧＲに達する。これに伴い、実筒内圧Ｐ＿ＣＰＳが減少することで、圧力差ΔＰは、同図に矢印Ｂで示すように、ＥＧＲ係数Ｃ＿ＥＧＲに対応する傾きで線形に減少し、ＥＧＲ率Ｒ＿ＥＧＲ＝目標ＥＧＲ率ＴＧＴ＿ＥＧＲになったときに、０になる。したがって、この間の任意の点において圧力差ΔＰ＝ΔＰ２が得られた場合、上記のΔＰ１の場合と同様、式（Ａ）によって、そのときの実ＥＧＲ率Ｒ＿ＥＧＲ２を精度良く算出することができる。
【００９３】
次に、図１８を参照しながら、ＥＧＲ率Ｒ＿ＥＧＲを用いた点火時期の制御処理について説明する。本処理は、ＴＤＣ信号の発生に同期して、気筒３ａごとに実行される。本処理では、まずステップ５１において、エンジン回転数ＮＥ及び要求トルクＴＲＱＣＭＤに応じ、所定のマップ（図示せず）を検索することによって、基本点火時期ＩＧ＿ＢＡＳＥを算出する。
【００９４】
次に、推定されたＥＧＲ率Ｒ＿ＥＧＲに応じ、所定のマップ（図示せず）を検索することによって、ＥＧＲ補正量ΔＩＧＥＧＲを算出する（ステップ５２）。
【００９５】
次に、エンジン水温ＴＷやエンジン回転数ＮＥなどに応じて、ＥＧＲ率Ｒ＿ＥＧＲ以外の要因による補正量ΔＩＧＴＴＬを算出する（ステップ５３）。
【００９６】
最後に、基本点火時期ＩＧ＿ＢＡＳＥに、ＥＧＲ補正量ΔＩＧＥＧＲ及び補正量ΔＩＧＴＴＬを加算することによって、点火時期ＩＧＬＯＧを算出し（ステップ５４）、本処理を終了する。
【００９７】
以上のように、本実施形態によれば、ＥＧＲ率Ｒ＿ＥＧＲが目標ＥＧＲ率ＴＧＴ＿ＥＧＲに等しく且つ空燃比が理論空燃比であるという混合気の組成の条件で、混合気の比熱比κの温度特性に基づき、基準クランク角ＣＡ＿ＲＥＦにおいて発生する基準筒内圧Ｐ＿ＲＥＦを算出する。そして、基準クランク角ＣＡ＿ＲＥＦにおいて検出された実筒内圧Ｐ＿ＣＰＳと基準筒内圧Ｐ＿ＲＥＦとの圧力差ΔＰに基づいて、ＥＧＲ率Ｒ＿ＥＧＲを算出するので、混合気の比熱比κの温度特性を反映させながら、ＥＧＲ率Ｒ＿ＥＧＲを推定できる。
【００９８】
また、基準クランク角ＣＡ＿ＲＥＦは、混合気の燃焼の開始直前のクランク角であるので、燃焼がまだ行われず、混合気の状態変化がポリトロープ変化に保たれた状態で、実筒内圧Ｐ＿ＣＰＳを取得するとともに、実筒内圧Ｐ＿ＣＰＳと基準筒内圧Ｐ＿ＲＥＦとの大きな圧力差ΔＰを確保できる。したがって、この圧力差ΔＰに基づき、混合気の比熱比κの温度特性を良好に反映させながら、ＥＧＲ率ＥＧＲ率Ｒ＿ＥＧＲを精度良く推定することができる。また、精度良く推定されたＥＧＲ率Ｒ＿ＥＧＲを用いて、点火時期ＩＧＬＯＧの制御を適切に行うことができる。
【００９９】
また、ＥＧＲ率Ｒ＿ＥＧＲの推定を、気筒３ａ内の実際の圧力である実筒内圧Ｐ＿ＣＰＳとその基準筒内圧Ｐ＿ＲＥＦをパラメータとして行うので、ＥＧＲ装置１４が低圧ＥＧＲ装置の場合においても、外部ＥＧＲガスの遅れの影響を受けることなく、ＥＧＲ率Ｒ＿ＥＧＲを精度良く推定することができる。
【０１００】
さらに、点火時期ＩＧＬＯＧ、吸気圧ＰＢＡ及びエンジン回転数ＮＥを用いて、基準クランク角ＣＡ＿ＲＥＦを設定するので、その設定を、エンジン３の実際の運転状態に応じて適切に行うことができ、基準クランク角ＣＡ＿ＲＥＦにおける基準筒内圧Ｐ＿ＲＥＦ及び実筒内圧Ｐ＿ＣＰＳを適切に得ることができる。
【０１０１】
また、基準筒内圧Ｐ＿ＲＥＦを、基準クランク角ＣＡ＿ＲＥＦ、圧縮開始時の初期クランク角に相当する吸気閉弁タイミングＩＶＣ、初期筒内温度Ｔ＿ＳＴＲＴ、及び初期筒内圧に相当する吸気圧ＰＢＡと、目標ＥＧＲ率ＴＧＴ＿ＥＧＲに応じて、適切に算出することができる。さらに、算出された基準筒内圧Ｐ＿ＲＥＦをエンジン回転数ＮＥ及びエンジン水温ＴＷに応じて補正することにより、気筒３ａ内と外部との間で授受される熱の影響を補償することができる。
【０１０２】
また、基準筒内圧Ｐ＿ＲＥＦの算出に用いたパラメータのうち、目標ＥＧＲ率ＴＧＴ＿ＥＧＲ以外の４つのパラメータ（基準クランク角ＣＡ＿ＲＥＦ、吸気閉弁タイミングＩＶＣ、初期筒内温度Ｔ＿ＳＴＲＴ及び吸気圧ＰＢＡ）に応じ、混合気の吸気・圧縮条件を反映させながら、ＥＧＲ係数Ｃ＿ＥＧＲを適切に算出できる。そして、算出したＥＧＲ係数Ｃ＿ＥＧＲを圧力差ΔＰに乗算し、さらに目標ＥＧＲ率ＴＧＴ＿ＥＧＲを加算することによって、ＥＧＲ率Ｒ＿ＥＧＲを算出するので、ＥＧＲ率Ｒ＿ＥＧＲの推定を精度良く行うことができる。
【０１０３】
さらに、筒内圧センサ５１は、その圧力検出素子及び増幅回路がインジェクタ４に一体に設けられているため、点火動作によるノイズや他の気筒３ａのインジェクタ４の噴射動作によるノイズの影響を受けにくい。このため、筒内圧センサ５１による実筒内圧Ｐ＿ＣＰＳの検出精度が高められることで、ＥＧＲ率Ｒ＿ＥＧＲの推定精度をさらに向上させることができる。
【０１０４】
次に、図１９を参照しながら、基準筒内圧Ｐ＿ＲＥＦの算出処理の変形例について説明する。この変形例は、前述したように、基準筒内圧Ｐ＿ＲＥＦが吸気圧ＰＢＡに比例するという関係（図１２）から、基準筒内圧マップの入力パラメータから吸気圧ＰＢＡを除外し、基準筒内圧マップで得られたマップ値を吸気圧ＰＢＡで補正するようにしたものである。本処理は、図６の処理に代えて実行される。また、図１９において、図６と同じ処理内容のステップには、同じステップ番号が付されている。
【０１０５】
本処理では、図６の処理と同じステップ３１及びステップ３２を実行し、吸気閉弁タイミングＩＶＣ及び初期筒内温度Ｔ＿ＳＴＲＴを算出する。次に、基準クランク角ＣＡ＿ＲＥＦ、吸気閉弁タイミングＩＶＣ、初期筒内温度Ｔ＿ＳＴＲＴ及び目標ＥＧＲ率ＴＧＴ＿ＥＧＲに応じ、基準筒内圧マップ（図示せず）を検索することによって、基準筒内圧Ｐ＿ＲＥＦを算出する（ステップ３０１）。なお、この基準筒内圧マップでは、圧縮開始時の初期筒内圧は、定数として扱われ、基準大気圧ＰＡＴＭ（７６０ｍｍＨｇ）が用いられている。
【０１０６】
次に、吸気圧ＰＢＡを基準大気圧ＰＡＴＭで除した値を、吸気圧補正係数Ｋ＿ＰＢとして設定する（ステップ３０２）とともに、この吸気圧補正係数Ｋ＿ＰＢをステップ３０１で算出された基準筒内圧Ｐ＿ＲＥＦに乗算することによって、補正された基準筒内圧Ｐ＿ＲＥＦを算出する（ステップ３０３）。
【０１０７】
その後の処理内容は、図６と同じであり、ステップ３０３で算出された基準筒内圧Ｐ＿ＲＥＦに、エンジン回転数ＮＥ及びエンジン水温ＴＷに応じて算出した伝熱補正係数Ｋ＿ＨＴを乗算することによって、最終的な基準筒内圧Ｐ＿ＲＥＦを算出し（ステップ３４、３５）、本処理を終了する。
【０１０８】
以上の変形例によれば、図６の算出処理の場合と同等の基準筒内圧Ｐ＿ＲＥＦを算出できるとともに、入力パラメータが削減されることで、基準筒内圧マップを容易に作成でき、その負荷を軽減することができる。
【０１０９】
なお、本発明は、説明した実施形態に限定されることなく、種々の態様で実施することができる。例えば、実施形態では、基準クランク角ＣＡ＿ＲＥＦを算出する際のパラメータとして、点火時期ＩＧＬＯＧ、吸気圧ＰＢＡ及びエンジン回転数ＮＥを用いているが、他の適当なパラメータを併せて用いてもよい。
【０１１０】
また、実施形態では、基準筒内圧Ｐ＿ＲＥＦ及びＥＧＲ係数Ｃ＿ＥＧＲの算出に用いる初期筒内温度Ｔ＿ＳＴＲＴを、吸気温ＴＡ、吸気位相ＣＡＩＮ及び排気位相ＣＡＥＸに応じて算出しているが、吸排気弁のバルブオーバーラップによる内部ＥＧＲが実行されない場合には、吸気温ＴＡをそのまま初期筒内温度としてもよい。さらに、初期筒内圧として、吸気圧ＰＢＡを用いたが、圧縮開始時に筒内圧センサ５１で検出された筒内圧ＰＣＹＬを用いることも可能である。
【０１１１】
また、基準筒内圧Ｐ＿ＲＥＦを、エンジン回転数ＮＥ及びエンジン水温ＴＷに応じて補正しているが、気筒３ａの内外間の熱の授受に影響を及ぼす他の適当なパラメータをさらに用いて、補正を行ってもよい。
【０１１２】
さらに、実施形態では、推定したＥＧＲ率Ｒ＿ＥＧＲに応じて、点火時期制御を実行しているが、これに代えて又はこれとともに、他のエンジン制御、例えばＥＧＲ弁４２を介したＥＧＲ制御、スロットル弁２７を介した吸入空気量制御や、インジェクタ４を介した燃料噴射制御などを実行してもよい。
【０１１３】
また、実施形態では、ＥＧＲ装置１４は、低圧ＥＧＲ装置で構成されているが、これに代えて又はこれとともに、高圧ＥＧＲ装置を用いてもよく、その場合にも前述した効果を同様に得ることができる。さらに、筒内圧センサ５１は、インジェクタ４と一体型のものであるが、インジェクタ４と分離して配置される別体型のものでもよいことは、もちろんである。
【０１１４】
さらに、実施形態では、エンジン３は車両用のエンジンであるが、本発明は、他の用途のエンジン、例えばクランクシャフトを鉛直方向に配置した船外機用のエンジンなどにも適用可能である。その他、本発明の趣旨の範囲内で、細部の構成を適宜、変更することができる。
【符号の説明】
【０１１５】
２  ＥＣＵ（目標ＥＧＲ率設定手段、基準クランク角設定手段、基準筒内圧算出手段、      ＥＧＲ率推定手段、制御手段、初期クランク角取得手段、初期筒内温度取得手段）
３  内燃機関
３ａ  気筒
４  燃料噴射弁
６  吸気通路
７  排気通路
１３  ターボチャージャ（過給機）
１４  ＥＧＲ装置
２１  コンプレッサ
２３  タービン
５１  筒内圧センサ
５２  クランク角センサ（運転状態検出手段、回転数検出手段）
５３  吸気位相センサ（初期クランク角取得手段、初期筒内温度取得手段）
５４  排気位相センサ（初期筒内温度取得手段）
５６  吸気圧センサ（運転状態検出手段、初期筒内圧取得手段）
５７  吸気温センサ（初期筒内温度取得手段）
５９  水温センサ（冷却水温度検出手段）
κ    混合気の比熱比
ＴＧＴ＿ＥＧＲ    目標ＥＧＲ率
ＰＣＹＬ    筒内圧
ＣＡ    クランク角
ＣＡ＿ＲＥＦ    基準クランク角
Ｐ＿ＲＥＦ    基準筒内圧
Ｐ＿ＣＰＳ    実筒内圧
ΔＰ    実筒内圧と基準筒内圧との圧力差
Ｒ＿ＥＧＲ    ＥＧＲ率
ＩＧＬＯＧ    点火時期
ＰＢＡ    吸気圧（吸気の圧力、初期筒内圧）
ＮＥ    エンジン回転数（内燃機関の回転数）
ＩＶＣ    吸気閉弁タイミング（初期クランク角）
Ｔ＿ＳＴＲＴ    初期筒内温度
ＴＷ    エンジン水温（冷却水の温度）
Ｃ＿ＥＧＲ    ＥＧＲ係数
【図１】
0006456256.tif 000008
【図２】
0006456256.tif 000009
【図３】
0006456256.tif 000010
【図４】
0006456256.tif 000011
【図５】
0006456256.tif 000012
【図６】
0006456256.tif 000013
【図７】
0006456256.tif 000014
【図８】
0006456256.tif 000015
【図９】
0006456256.tif 000016
【図１０】
0006456256.tif 000017
【図１１】
0006456256.tif 000018
【図１２】
0006456256.tif 000019
【図１３】
0006456256.tif 000020
【図１４】
0006456256.tif 000021
【図１５】
0006456256.tif 000022
【図１６】
0006456256.tif 000023
【図１７】
0006456256.tif 000024
【図１８】
0006456256.tif 000025
【図１９】
0006456256.tif 000026
【図２０】
0006456256.tif 000027
【図２１】
0006456256.tif 000028