(57)【要約】  本発明の課題は、過給エンジンの筒内吸入空気量を予測し、その予測値に基づいて過給エンジンを制御する制御装置において、特に過渡運転時の筒内吸入空気量の予測精度を高めることである。このため、本制御装置は、過給エンジンの物理モデルを用いてスロットル開度の予測値から過給圧の予測値を計算するとともに、その補正量を計算する。補正量の計算では、過給圧センサによる過給圧の計測値が取得されるとともに、過給エンジンの物理モデルを用いてスロットル開度の計測値から過給圧の推定値を計算することが行われる。そして、過給圧の計測値と推定値との差分が計算され、その差分が過給圧の予測値の補正量とされる。本制御装置は、過給圧の補正された予測値とスロットル開度の予測値とに基づいて筒内吸入空気量の予測値を計算する。
2013084318.tif 000002
【特許請求の範囲】
【請求項１】
過給エンジンの筒内吸入空気量を予測し、その予測値に基づいて前記過給エンジンを制御する制御装置において、
現時点よりも所定時間先の時点におけるスロットル開度の予測値を取得する手段と、
前記過給エンジンの物理モデルを用いて、スロットル開度の前記予測値から過給圧の予測値を計算する手段と、
過給圧センサによる過給圧の計測値を取得する手段と、
スロットル開度センサによるスロットル開度の計測値を取得する手段と、
前記物理モデルを用いて、スロットル開度の前記計測値から過給圧の推定値を計算する手段と、
過給圧の前記計測値と前記推定値との差分を計算し、当該差分を補正量として過給圧の前記予測値を補正する手段と、
過給圧の前記補正された予測値とスロットル開度の前記予測値とに基づいて筒内吸入空気量の前記予測値を計算する手段と、
を備えることを特徴とする過給エンジンの制御装置。
【請求項２】
前記の筒内吸入空気量の予測値を計算する手段は、
前記物理モデルの一部を構成する第１のサブモデルを用いて、過給圧の前記補正された予測値とスロットル開度の前記予測値とから吸気マニホールド圧の予測値を計算する手段と、
流量センサによる吸気流量の計測値を取得する手段と、
前記物理モデルの一部を構成する第２のサブモデルを用いて、吸気流量の前記計測値から吸気マニホールド圧の第１の推定値を計算する手段と、
前記流量センサの応答モデルと前記第２のサブモデルとを用いて、過給圧の前記推定値を計算する過程において前記物理モデルで計算される所定の吸気通路内流量の値から吸気マニホールド圧の第２の推定値を計算する手段と、
吸気マニホールド圧の前記第１の推定値と前記第２の推定値との差分を計算し、当該差分を補正量として吸気マニホールド圧の前記予測値を補正する手段と、
を備え、吸気マニホールド圧の前記補正された予測値に基づいて筒内吸入空気量の前記予測値を計算するように構成されていることを特徴とする請求項１に記載の過給エンジンの制御装置。
【請求項３】
前記の筒内吸入空気量の予測値を計算する手段は、
前記第１のサブモデルを用いて、過給圧の前記計測値とスロットル開度の前記計測値とから吸気マニホールド圧の第３の推定値を計算する手段と、
吸気マニホールド圧の前記第３の推定値と、過給圧の前記推定値を計算する過程において前記物理モデルで計算される吸気マニホールド圧の第４の推定値との差分を計算し、当該差分を補正量として吸気マニホールド圧の前記補正された予測値をさらに補正する手段と、
をさらに備え、吸気マニホールド圧の前記さらに補正された予測値に基づいて筒内吸入空気量の前記予測値を計算するように構成されていることを特徴とする請求項２に記載の過給エンジンの制御装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、過給エンジンの制御装置に関し、詳しくは、筒内吸入空気量の予測値に基づいて過給エンジンを制御する制御装置に関する。
【背景技術】
【０００２】
今日の自動車のエンジンには電子制御式のスロットルが用いられている。電子制御式のスロットルを備えたエンジンでは、ドライバのアクセル操作量に基づいて目標スロットル開度が決定され、目標スロットル開度に従ってスロットルが操作される。このとき、決定した目標スロットル開度を直ぐにスロットルに与えるのではなく、ある時間だけ遅らせてスロットルに与えるようにすることもできる。このような演算処理をスロットルディレイ制御という。スロットルディレイ制御によれば、目標スロットル開度の変化に対し実際のスロットル開度はディレイ時間分だけ遅れて変化することになるので、ディレイ時間分だけ将来のスロットル開度を目標スロットル開度から予測することができる。将来のスロットル開度が予測可能であることは、吸気ポートに燃料を噴射するポート噴射型エンジンやポート噴射と筒内直接噴射を併用するエンジンにとって空燃比の制御精度を高める上で有用である。筒内吸入空気量は吸気弁の閉時点において確定するが、ポート噴射を行うエンジンの場合、吸気弁の閉時点の前に燃料噴射の開始時期が到来する。このため、目標空燃比の実現に必要な燃料噴射量を正確に計算するためには、将来確定する筒内吸入空気量を燃料噴射の開始時点において予測することが必要となる。スロットルディレイ制御によれば、将来の吸気弁の閉時点でのスロットル開度を目標スロットル開度から予測することができるので、その予測スロットル開度に基づいて筒内吸入空気量を予測することができる。【０００３】
予測スロットル開度に基づく筒内吸入空気量の計算では、吸気通路内の空気の挙動を物理的にモデル化した物理モデルが用いられる。自然吸気エンジンの場合、筒内吸入空気量の計算のための物理モデルは、スロットルモデル、吸気管モデル、及び吸気弁モデルで構成することができる。スロットルモデルは、スロットルを通過する空気の流量を算出するためのモデルである。具体的には、スロットルの前後の差圧とスロットル開度により決まる流路面積と流量係数とを基本とするオリフィスの流量式がスロットルモデルとして用いられる。吸気管モデルは、吸気管内の空気に関する保存則に基づいて構築されたモデルである。具体的には、エネルギー保存則の式と流量保存則の式とが吸気管モデルとして用いられる。吸気弁モデルは、吸気弁流量と吸気管圧との関係について調べた実験ベースのモデルである。実験で得られた経験則により、吸気弁モデルにおいては吸気弁流量と吸気管圧との関係が１又は複数の直線で近似されている。【０００４】
上記の物理モデルを用いた計算では、モデル化誤差の影響を如何に排除するかが筒内吸入空気量の予測精度を担保する上で重要である。物理モデルによる筒内吸入空気量の予測精度を高める方法の一例としては、特許第３７６０７５７号公報に開示された方法を挙げることができる。この公報に開示された方法によれば、スロットル開度に基づき演算されるスロットル流量に基づいて現在から所定時間経過後の第１吸気管圧力が算出されるとともに、エアフローメータの出力に基づいて第２吸気管圧力が算出される。また、現在のスロットル流量をエアフローメータモデルに入力することによって時間遅れを含むエアフローメータ出力が演算され、そのエアフローメータ出力に基づいて第２吸気管圧力と同応答の吸気管圧力が算出される。そして、第１吸気管圧力と第２吸気管圧力を加算したものから同応答吸気管圧力を減じることにより予測圧力が算出され、その予測圧力に基づいて吸気弁流量の推定値が算出される。この方法により算出された吸気弁流量の推定値は、定常運転時にはエアフローメータの出力に一致する。つまり、上記公報に開示された方法によれば、定常運転時のモデル化誤差が補償されることによって筒内吸入空気量の予測精度が高められる。【０００５】
物理モデルを用いた筒内吸入空気量の予測方法は、自然吸気エンジンだけでなく、ターボ過給機や機械式過給機を備えた過給エンジンにも適用することができる。ただし、過給エンジンと自然吸気エンジンとではスロットル上流圧力に関する前提が大きく異なっている。自然吸気エンジンの場合、スロットル上流圧力は大気圧とほぼ等しいとみなすことができる。ところが、過給エンジンの場合、スロットル上流圧力はコンプレッサの回転数やスロットル開度によって変化する。吸気弁流量の推定にはスロットル流量の計算が必要であり、スロットル流量の計算にはスロットル上流圧力が用いられる。図５は、スロットル上流圧力に対するスロットル下流圧力の圧力比“Pm/Pic”とスロットル流量“mt”との関係を示す図である。この図に示すように、スロットル開度が“TA1”から“TA2”に変化したときのスロットル流量“mt”の変化の感度は圧力比“Pm/Pic”に依存する。このため、計算に用いるスロットル上流圧力の精度が低ければスロットル流量の計算精度が低くなり、ひいては筒内吸入空気量の予測精度も低くなる。よって、物理モデルを用いた筒内吸入空気量の予測手法を過給エンジンに適用する上では、スロットル上流圧力に関する精度の高い情報が必要となる。【０００６】
過給エンジンのスロットル上流圧力、すなわち、過給圧は、吸気通路に取り付けた過給圧センサによって計測することができる。よって、過給圧センサによる過給圧の計測値を取得し、それを物理モデルに入力することが１つの方法として考えられる。しかし、この方法には筒内吸入空気量の予測精度、特に、エンジンの過渡運転時における予測精度の点において問題がある。物理モデルで算出される筒内吸入空気量は将来の予測値であるのに対し、過給圧センサによって得られる過給圧は現時点での計測値である。エンジンの過渡運転時には、スロットルやウエストゲートバルブ等のアクチュエータの動作によって過給圧は大きく変化する。このため、過給圧の計測時点から筒内吸入空気量の予測時点までの間に過給圧が変化する状況は容易に起こりうる。その場合、筒内吸入空気量の予測値と実際値（吸気弁の閉時点での実際値）との間には、過給圧の変化の分だけ誤差が生じることになる。【０００７】
一方、前述の物理モデルを用いれば、予測スロットル開度から将来の過給圧を予測することができる。吸気弁の閉時点での過給圧を予測できれば、その予測値と予測スロットル開度から筒内吸入空気量の予測値を得ることができる。しかし、物理モデルによる計算には、物理モデル作成時の適合誤差、エンジン構成部品の製造ばらつき、環境変化や経年変化、アクチュエータのヒステリシス誤差などの様々な要因によるモデル化誤差が含まれうる。よって、物理モデルによる過給圧の予測値から筒内吸入空気量の予測値を計算する場合、筒内吸入空気量の予測値と実際値（吸気弁の閉時点での実際値）との間にはモデル化誤差に起因する誤差が生じる可能性がある。【０００８】
以上述べたように、物理モデルを用いて過給エンジンの筒内吸入空気量を精度良く予測するためには、エンジンの運転状態によって変化する過給圧に関する情報が必要となる。しかし、単に過給圧センサによる過給圧の計測値を入力したとしても、また、単に物理モデルによる過給圧の予測値を入力したとしても、それのみでは筒内吸入空気量、特に、過渡運転時の筒内吸入空気量を高い精度で予測することは難しい。
【先行技術文献】
【特許文献】
【０００９】
【特許文献１】特許第３７６０７５７号公報
【発明の概要】
【００１０】
本発明の課題は、過給エンジンの筒内吸入空気量を予測し、その予測値に基づいて過給エンジンを制御する制御装置において、特に過渡運転時の筒内吸入空気量の予測精度を高めることである。このような課題を達成するためのアプローチとして、本発明では、筒内吸入空気量の予測値の計算には物理モデルによる過給圧の予測値を使用するが、それに含まれるモデル化誤差を補償するために過給圧センサによる過給圧の計測値を使用する。つまり、過給圧センサによる過給圧の計測値と物理モデルによる過給圧の予測値とを相互補完的に利用する。
【００１１】
具体的には、本発明の１つの形態によれば、本制御装置は、現時点よりも所定時間先の時点におけるスロットル開度の予測値を取得する。そして、過給エンジンの物理モデルを用いて、スロットル開度の予測値から過給圧の予測値を計算する。また、それと並行して、本制御装置は、過給圧センサによる過給圧の計測値を取得するとともに、スロットル開度センサによるスロットル開度の計測値を取得し、前記物理モデルを用いてスロットル開度の計測値から過給圧の推定値を計算する。過給圧センサによる過給圧の計測値と物理モデルによる過給圧の推定値はどちらも同時刻における過給圧情報であるから、理論的には両者は一致するはずである。両者の間に差分が生じるのであれば、その差分はモデル化誤差に起因したものと言える。本制御装置は、過給圧の計測値と推定値との差分を計算し、当該差分を補正量として過給圧の予測値を補正する。計算上、補正された過給圧の予測値は定常運転時には過給圧の計測値に等しくなる。これは、過給圧の予測値に含まれる圧力誤差が前記補正量による補正によって補償されたことを示している。本制御装置は、このように補正された過給圧の予測値とスロットル開度の予測値とに基づいて筒内吸入空気量の予測値を計算する。
【００１２】
筒内吸入空気量の予測値を計算する過程では、流量センサによる吸気流量の計測値をモデル化誤差の補償のための付加的な情報として使用することができる。この場合、本制御装置は、前記物理モデルの一部を構成する第１のサブモデルを用いて、過給圧の補正された予測値とスロットル開度の予測値とから吸気管圧の予測値を計算する。また、流量センサによる吸気流量の計測値を取得し、前記物理モデルの一部を構成する第２のサブモデルを用いて、吸気流量の計測値から吸気管圧の第１の推定値を計算する。さらに、本制御装置は、流量センサの応答モデルと第２のサブモデルとを用いて、前記物理モデルで計算される所定の吸気通路内流量（例えばコンプレッサ流量）の値から吸気管圧の第２の推定値を計算する。流量センサの応答モデルは、流量センサによる吸気流量の計測値の実際値に対する応答遅れをモデル化したものである。よって、物理モデルで計算される吸気通路内流量を流量センサの応答モデルで処理した値は、理論的には流量センサによる吸気流量の計測値と一致するはずである。さらに、それらから計算される吸気管圧の第１の推定値と第２の推定値も、理論的には互いに一致するはずである。吸気管圧の第１の推定値と第２の推定値との間に差分が生じるのであれば、モデル化誤差に起因する流量誤差が原因であると言える。本制御装置は、吸気管圧の第１の推定値と第２の推定値との差分を計算し、当該差分を補正量として吸気管圧の予測値を補正する。本制御装置は、このように補正された吸気管圧の予測値に基づいて筒内吸入空気量の予測値を計算する。
【００１３】
上記の補正によってモデル化誤差に起因する流量誤差は補償されるが、過給圧の計測値と推定値との間に誤差が生じている状況では、吸気管圧の補正された予測値と第１の推定値とが定常運転時において完全には一致しない可能性がある。筒内吸入空気量の予測精度をさらに高めるのであれば、上記の計算処理に加えて次の計算処理を付加的に実施することが好ましい。その付加的な処理によれば、本制御装置は、前記第１のサブモデルを用いて、過給圧の計測値とスロットル開度の計測値とから吸気管圧の第３の推定値を計算する。そして、吸気管圧の第３の推定値と、過給圧の推定値を計算する過程において前記物理モデルで計算される吸気管圧の第４の推定値との差分を計算し、当該差分を補正量として吸気管圧の補正された予測値をさらに補正する。この補正により、定常運転時には吸気管圧の予測値は第３の推定値と相殺され、吸気管圧の第２の推定値は第４の推定値と相殺される。その結果、計算上、補正後の吸気管圧の予測値は流量センサによる計測値に基づく第１の推定値と等しくなる。これは、モデル化誤差に起因する流量誤差と圧力誤差の両方が補正によって補償されたことを示している。本制御装置は、このように補正された吸気管圧の予測値に基づいて筒内吸入空気量の予測値を計算する。
【図面の簡単な説明】
【００１４】
【図１】本発明の実施の形態の過給エンジンの制御装置で用いられている筒内吸入空気量予測モデルを示すブロック図である。
【図２】本発明の実施の形態の過給エンジンの制御装置で用いられている筒内吸入空気量予測モデルを示すブロック図である。
【図３】本発明の実施の形態の過給エンジンの制御装置で用いられている筒内吸入空気量予測モデルを示すブロック図である。
【図４】本発明の実施の形態の過給エンジンの制御装置で用いられている筒内吸入空気量予測モデルを示すブロック図である。
【図５】スロットル開度の変化の直後におけるスロットル流量の変化の感度がスロットル前後圧力比に依存していることを示す図である。
【発明を実施するための形態】
【００１５】
本発明の実施の形態について図を参照して説明する。
【００１６】
本実施の形態の制御装置が適用される過給エンジンは、スロットルによる空気量の調整によってトルクを制御する火花点火式の４サイクルレシプロエンジンである。本実施の形態の過給エンジンが備える過給機は、吸気通路に配置されたコンプレッサを排気通路に配置されたタービンの回転によって駆動するターボ式の過給機である。タービンにはその開度を能動的に制御することのできるウエストゲートバルブが付設されている。コンプレッサとスロットルとの間には、コンプレッサによる圧縮によって温度が上昇した空気を冷却するためのインタークーラが設けられている。また、本実施の形態の過給エンジンは、各気筒の吸気ポートに燃料を噴射するポート噴射型のエンジンでもある。
【００１７】
本実施の形態の制御装置は、過給エンジンを制御するＥＣＵ（Electronic Control Unit）の機能の一部として実現される。ＥＣＵには、流量センサであるエアフローメータ、過給圧センサ、スロットル開度センサ等の各種のセンサから、エンジンの運転状態や運転条件に関する様々な情報や信号が入力される。ＥＣＵは、それら情報や信号に基づいてスロットルやウエストゲートバルブ等の各種のアクチュエータを操作する。
【００１８】
制御装置としてのＥＣＵは、筒内吸入空気量を予測する機能を有している。ポート噴射型エンジンの場合、吸気弁が閉弁して筒内吸入空気量が確定するよりも前に必要な燃料噴射量を計算し、燃料噴射を開始しなければならない。このため、燃料噴射量を計算する時点において、将来確定する筒内吸入空気量を予測することが必要とされる。ＥＣＵによる筒内吸入空気量の予測には、プログラムされている筒内吸入空気量予測モデルが用いられる。筒内吸入空気量予測モデルは、過給エンジンにおける空気の挙動を物理的にモデル化したものであって、その概要は図１、図２、図３及び図４の各ブロック図によって表される。
【００１９】
図１−図４に示すように、本実施の形態で用いられる筒内吸入空気量予測モデルは、符号２，４，６，８，１０，１２，１４，１６，１８で示される８つの計算ブロックから構成されている。以下、図１から順に、各図に示す計算ブロックの構成と機能について説明する。
【００２０】
図１には３つの計算ブロック２，４，６が示されている。計算ブロック２は、スロットル開度の予測値から将来の時点における過給圧の予測値を計算するための計算ブロックである。スロットル開度の予測値は、ＥＣＵにより別途実施されるスロットルディレイ制御において、アクセル操作量に基づき決定された目標スロットル開度から先読みされる。本実施の形態では、現時点よりも所定制御周期分だけ将来の時点におけるスロットル開度が目標スロットル開度から予測されているものとする。
【００２１】
計算ブロック２は１つの物理モデルであり、複数の要素モデル、すなわち、ターボ回転数モデルＭ１、コンプレッサモデルＭ２、インタークーラモデルＭ３、スロットルモデルＭ４、吸気管モデルＭ５、及び吸気弁モデルＭ６によって構成されている。以下、計算ブロック２が備えるこれら要素モデルの内容について説明する。ただし、これらの要素モデルはそれぞれに公知であり、また、それ自体は本発明における特徴点ではないことから、数式やマップ等の各要素モデルの詳細については記載を省略する。
【００２２】
ターボ回転数モデルＭ１は、ターボ過給機の回転挙動のモデルであって、吸気弁流量とウエストゲートバルブ開度とターボ回転数との間に成り立つ関係がモデル化されている。ターボ回転数モデルＭ１は、数式或いは実験データに基づくマップによって構成されている。ターボ回転数モデルＭ１では、ウエストゲートバルブの操作量から推定されたウエストゲートバルブ開度“wgv”と、後述する吸気弁モデルＭ６で算出された吸気弁流量“mc”とが入力され、それらの入力情報からターボ回転数“Ntb”が算出される。
【００２３】
コンプレッサモデルＭ２は、ターボ過給機のコンプレッサのモデルであって、ターボ回転数と過給圧とコンプレッサ流量との間に成り立つ関係がモデル化されている。コンプレッサモデルＭ２は、数式或いは実験データに基づくマップによって構成されている。コンプレッサモデルＭ２では、ターボ回転数モデルＭ１で算出されたターボ回転数“Ntb”と、後述するインタークーラモデルＭ３で算出された過給圧“Pic(epicvlv）”等の情報が入力され、それらの入力情報からコンプレッサ流量“mcp”が算出される。なお、過給圧を示す符号である“Pic”に付された“epicvlv”は、他の計算ブロック内で計算される過給圧と区別するための符号である。図中及び本明細書内では、過給圧以外の他のパラメータに関しても、必要に応じてそのような区別用の符号を付したり、或いは、そのような区別用の符号のみでパラメータの内容を表したりしている。
【００２４】
インタークーラモデルＭ３は、吸気通路におけるインタークーラ内の空気に関する保存則に基づいて構築された物理モデルである。インタークーラモデルＭ３としては、具体的にはエネルギー保存則の式と流量保存則の式とが用いられている。インタークーラモデルＭ３では、コンプレッサモデルＭ２で算出されたコンプレッサ流量“mcp”と後述するスロットルモデルＭ４で算出されたスロットル流量“mt”等の情報が入力され、それらの入力情報からスロットル上流圧力としての過給圧“Pic”が算出される。
【００２５】
スロットルモデルＭ４は、スロットルを通過する空気の流量を算出するためのモデルであって、具体的には、スロットルの前後の差圧、スロットル開度により決まる流路面積、及び流量係数を基本とするオリフィスの流量式が用いられている。スロットルモデルＭ４では、スロットルディレイ制御によって先読みされたスロットル開度予測値“eta0”、インタークーラモデルＭ３で算出された過給圧“Pic”、及び後述する吸気管モデルＭ５で算出された吸気管圧“Pm(epmvlv)”等の情報が入力され、それらの入力情報からスロットル流量“mt”が算出される。
【００２６】
吸気管モデルＭ５は、吸気管内の空気に関する保存則に基づいて構築された物理モデルである。吸気管モデルＭ５としては、具体的にはエネルギー保存則の式と流量保存則の式とが用いられている。吸気管モデルＭ５では、スロットルモデルＭ４で算出されたスロットル流量“mt”、及び後述する吸気弁モデルＭ６で算出された吸気弁流量“mc”等の情報が入力され、それらの入力情報から吸気管圧“Pm”が算出される。
【００２７】
吸気弁モデルＭ６は、吸気弁流量と吸気管圧との関係について調べた実験ベースのモデルである。実験で得られた経験則により、吸気弁モデルＭ６においては吸入空気量と吸気管圧との関係が直線で近似されている。その直線の方程式の係数は定数ではなく、エンジン回転数、ウエストゲートバルブ開度、吸気弁のバルブタイミング、排気弁のバルブタイミング等によって決まる変数である。吸気弁モデルＭ６では、吸気管モデルＭ５で算出された吸気管圧“Pm”の他、エンジン回転数“NE”、ウエストゲートバルブ開度“wgv”、吸気弁のバルブタイミング“InVT”、排気弁のバルブタイミング“ExVT”等の情報が入力され、それらの入力情報から吸気弁流量“mc(eklvlv)”が算出される。
【００２８】
ＥＣＵは、計算ブロック２で計算される各種パラメータのうち、インタークーラモデルＭ３で算出された過給圧“epicvlv”を取り出す。この過給圧“epicvlv”は、スロットル開度の予測値“eta0”から計算されたものであり、所定制御周期分だけ将来の時点における過給圧の予測値である。ＥＣＵは、取り出した過給圧予測値“epicvlv”を後述する計算ブロック６に入力する。
【００２９】
計算ブロック４は、スロットル開度センサによるスロットル開度の計測値から現時点における過給圧の推定値を計算するための計算ブロックである。計算ブロック４は、前述の計算ブロック２と同構成の物理モデルであって、ターボ回転数モデルＭ１、コンプレッサモデルＭ２、インタークーラモデルＭ３、スロットルモデルＭ４、吸気管モデルＭ５、及び吸気弁モデルＭ６によって構成されている。これらの要素モデルＭ１，Ｍ２，Ｍ３，Ｍ４，Ｍ５，Ｍ６間における情報の入出力は、計算ブロック２におけるそれと共通している。ただし、計算ブロック４では、スロットル開度センサによるスロットル開度の計測値“TA”がスロットルモデルＭ４に入力される。
【００３０】
ＥＣＵは、計算ブロック４で計算される各種パラメータのうち、インタークーラモデルＭ３で算出された過給圧“epiccrt”と、吸気管モデルＭ５で算出された吸気管圧“epmcrt”と、コンプレッサモデルＭ２で算出されたコンプレッサ流量“mcpcrt”とを取り出す。過給圧“epiccrt”は、現時点でのスロットル開度の計測値“TA”から計算されたものであり、現時点における過給圧の推定値である。同様に、吸気管圧“epmcrt”は現時点における吸気管圧の推定値であり、コンプレッサ流量“mcpcrt”は現時点におけるコンプレッサ流量の推定値である。ＥＣＵは、取り出した過給圧推定値“epiccrt”は次に説明する計算ブロック６に入力し、吸気管圧推定値“epmcrt”は後述する計算ブロック１４に入力し、コンプレッサ流量推定値“mcpcrt”は後述する計算ブロック１８に入力する。
【００３１】
計算ブロック６には、計算ブロック２で計算された過給圧予測値“epicvlv”と計算ブロック４で計算された過給圧推定値“epiccrt”とに加え、過給圧センサによる過給圧の計測値“gapic”が入力される。計算ブロック６では、過給圧計測値“gapic”と過給圧推定値“epiccrt”との差分が計算され、その差分が過給圧予測値“epicvlv”に加算される。補正後の過給圧予測値“epicvlv2”は、次の式（１）によって表される。
【数１】
2013084318.tif 000003
【００３２】
過給圧センサによる過給圧計測値“gapic”と計算ブロック４で計算された過給圧推定値“epiccrt”とはどちらも同時刻における過給圧を示している。したがって、理論的には両者は一致するはずであるが、計算ブロック４を構成する物理モデルがモデル化誤差を有する場合には、過給圧計測値“gapic”と過給圧推定値“epiccrt”との間には差分が生じる。この差分は、定常・過渡を問わずモデル化誤差に対する補正量として利用することができる。式（１）では、この差分を補正量として過給圧予測値“epicvlv”が補正される。定常運転時にはスロットル開度の計測値“TA”と予測値“eta0”とが等しくなることから、過給圧予測値“epicvlv”は過給圧推定値“epiccrt”に一致する。結果、式（１）から分かるように、定常運転時には補正後の過給圧予測値“epicvlv2”と過給圧計測値“gapic”とを一致させることができる。これにより、過渡運転時においても過給圧予測値が実際の値から大きく乖離することを防ぐことができ、筒内吸入空気量予測のモデル化誤差に対するロバスト性が高められる。
【００３３】
ＥＣＵは、計算ブロック６で算出された補正後の過給圧予測値“epicvlv2”を取り出し、図２に示す計算ブロック８に入力する。
【００３４】
図２には４つの計算ブロック８，１０，１２，１４が示されている。計算ブロック８は、補正後の過給圧予測値とスロットル開度の予測値とから吸気管圧の予測値を計算するための計算ブロックである。計算ブロック８は、計算ブロック２の物理モデルの一部を構成するサブモデル（第１のサブモデル）であって、スロットルモデルＭ４、吸気管モデルＭ５、及び吸気弁モデルＭ６によって構成されている。これらの要素モデルＭ４，Ｍ５，Ｍ６間における情報の入出力は、計算ブロック２におけるそれと共通している。ただし、計算ブロック８では、補正後の過給圧予測値“epicvlv2”とスロットル開度予測値“eta0”とがスロットルモデルＭ４に入力される。
【００３５】
ＥＣＵは、計算ブロック８で計算される各種パラメータのうち、吸気管モデルＭ５で算出された吸気管圧“epmvlv2”を取り出す。この吸気管圧“epmvlv2”は、スロットル開度の予測値“eta0”と過給圧の予測値“epicvlv2”とから計算されたものであり、所定制御周期分だけ将来の時点における吸気管圧の予測値である。ＥＣＵは、取り出した吸気管圧予測値“epmvlv2”を後述する計算ブロック１０に入力する。
【００３６】
計算ブロック１０には、計算ブロック８で計算された吸気管圧予測値“epmvlv2”に加え、後述する計算ブロック１６で計算される吸気管圧推定値“epmafm”と、後述する計算ブロック１８で計算される吸気管圧推定値“epmcrtsm”とが入力される。計算ブロック１０では、第１の推定値である吸気管圧推定値“epmafm”と第２の推定値である吸気管圧推定値“epmcrtsm”との差分が計算され、その差分が吸気管圧予測値“epmvlv2”に加算される。
【００３７】
ここで、計算ブロック１６の構成は図３に示されている。計算ブロック１６は、エアフローメータによる吸気流量の計測値から吸気管圧の推定値を計算するための計算ブロックである。計算ブロック１６は、計算ブロック２の物理モデルの一部を構成するサブモデル（第２のサブモデル）であって、インタークーラモデルＭ３、スロットルモデルＭ４、吸気管モデルＭ５、及び吸気弁モデルＭ６によって構成されている。これらの要素モデルＭ３，Ｍ４，Ｍ５，Ｍ６間における情報の入出力は、計算ブロック２におけるそれと共通している。ただし、計算ブロック１６では、エアフローメータによる吸気流量の計測値“egaafm”がインタークーラモデルＭ３に入力される。そして、吸気流量計測値“egaafm”に基づいて現時点における吸気管圧推定値“epmafm”が算出される。
【００３８】
計算ブロック１８の構成は図４に示されている。計算ブロック１８は、物理モデルにより推定したエアフローメータの出力値から吸気管圧の推定値を計算するための計算ブロックである。計算ブロック１８は、前述の計算ブロック１６が備えるインタークーラモデルＭ３、スロットルモデルＭ４、吸気管モデルＭ５、及び吸気弁モデルＭ６に加えて、エアフローメータモデルＭ７を含んでいる。エアフローメータは固有の応答特性に基づく応答遅れを有している。エアフローメータモデルＭ７は、エアフローメータの応答特性をシミュレートしたモデルであり、入力された吸気通路内流量に対してエアフローメータの応答遅れ分のディレイ処理を施す。計算ブロック１８のエアフローメータモデルＭ７には、計算ブロック４のコンプレッサモデルＭ２で算出されたコンプレッサ流量“mcpcrt”が入力される。エアフローメータモデルＭ７でディレイ処理されたコンプレッサ流量“mcpcrtsm”はインタークーラモデルＭ３に入力される。計算ブロック１８では、ディレイ処理後のコンプレッサ流量“mcpcrtsm”に基づいて現時点における吸気管圧推定値“epmcrtsm”が算出される。
【００３９】
計算ブロック１８で計算されるディレイ処理後のコンプレッサ流量“mcpcrtsm”は、理論的にはエアフローメータによる吸気流量の計測値“egaafm”に一致するはずである。しかし、実際には、ディレイ処理後のコンプレッサ流量“mcpcrtsm”と吸気流量計測値“egaafm”との間には、モデル化誤差に起因する流量誤差が存在する。このため、計算ブロック１６で計算された吸気管圧推定値“epmafm”と、計算ブロック１８で計算された吸気管圧推定値“epmcrtsm”との間には差分が生じる。この差分は、ディレイ処理後のコンプレッサ流量“mcpcrtsm”が吸気流量計測値“egaafm”に対して持っている流量誤差を吸気管圧の誤差に換算した値であるとみなすことができる。図２に示す計算ブロック１０では、上記の差分を補正量として吸気管圧予測値“epmvlv2”に加算する。補正後の吸気管圧予測値“epmfwd”は、次の式（２）によって表される。
【数２】
2013084318.tif 000004
【００４０】
ＥＣＵは、計算ブロック１０で計算された補正後の吸気管圧予測値“epmfwd”を計算ブロック１４に入力する。
【００４１】
計算ブロック１４には、計算ブロック１０で計算された補正後の吸気管圧予測値“epmfwd”に加え、後述する計算ブロック１２で計算される吸気管圧推定値“epmcrt2”と、計算ブロック４の吸気管モデルＭ５で算出された吸気管圧推定値“epmcrt”とが入力される。計算ブロック１４では、第３の推定値である吸気管圧推定値“epmcrt2”と第４の推定値である吸気管圧推定値“epmcrt”との差分が計算され、その差分が補正後の吸気管圧予測値“epmfwd”から減算される。
【００４２】
計算ブロック１２は、過給圧センサによる過給圧の計測値とスロットル開度センサによるスロットル開度の計測値とから吸気管圧の推定値を計算するための計算ブロックである。計算ブロック１２は、前述の計算ブロック８と同構成のサブモデルであって、スロットルモデルＭ４、吸気管モデルＭ５、及び吸気弁モデルＭ６によって構成されている。これらの要素モデルＭ４，Ｍ５，Ｍ６間における情報の入出力は、計算ブロック８におけるそれと共通している。ただし、計算ブロック１２では、過給圧計測値“gapic”とスロットル開度計測値“TA”とがスロットルモデルＭ４に入力される。そして、それらの情報に基づいて現時点における吸気管圧推定値“epmcrt2”が算出される。
【００４３】
過給圧センサによる過給圧計測値“gapic”と計算ブロック４で計算された過給圧推定値“epiccrt”との間に誤差が生じている状況では、計算ブロック１４に入力される吸気管圧推定値“epmcrt2”と吸気管圧推定値“epmcrt”との間にも差分が生じる。この差分は、過給圧推定値“epiccrt”が過給圧計測値“gapic”に対して持っている誤差を吸気管圧の誤差に換算した値であるとみなすことができる。計算ブロック１４では、上記の差分を補正量として吸気管圧予測値“epmfwd”から減算する。再補正後の吸気管圧予測値“epmfwd2”は、次の式（３）によって表される。
【数３】
2013084318.tif 000005
【００４４】
定常運転時にはスロットル開度の計測値“TA”と予測値“eta0”とが等しくなり、また、過給圧予測値“epicvlv2”と過給圧計測値“gapic”とが等しくなることから、同構成のモデルを用いて計算される吸気管圧予測値“epmvlv2”と第３の推定値である吸気管圧推定値 “epmcrt2”とは一致する。また、定常運転時には計算ブロック１８で計算される過給圧推定値“epiccrtsm”と計算ブロック４で計算される過給圧推定値“epiccrt”とが等しくなることから、第２の推定値である吸気管圧推定値 “epmcrtsm”と第４の推定値である吸気管圧推定値 “epmcrt”とは一致する。結果、式（３）から分かるように、定常運転時には、再補正後の吸気管圧予測値“epmfwd2”とエアフローメータによる計測値に基づく吸気管圧推定値“epmafm”とを一致させることができる。これにより、過渡運転時においても吸気管圧予測値が実際の値から大きく乖離することを防ぐことができ、筒内吸入空気量予測のモデル化誤差に対するロバスト性がさらに高められる。
【００４５】
ＥＣＵは、計算ブロック１４で算出された再補正後の吸気管圧予測値“epmfwd2”を吸気弁モデルＭ６に入力する。吸気弁モデルＭ６では再補正後の吸気管圧予測値“epmfwd2”から吸気弁流量“eklfwd2”が算出される。この吸気弁流量“eklfwd2”は、所定制御周期分だけ将来の時点における吸気弁流量の予測値である。ＥＣＵは、燃料噴射量の計算タイミングが到来したとき、吸気弁流量予測値“eklfwd2”に基づいて吸気弁の閉弁時点において予測される筒内吸入空気量を算出する。そして、筒内吸入空気量の予測値と目標空燃比とを用いて必要な燃料噴射量を計算する。
【００４６】
その他．
本発明は上述の実施の形態に限定されるものではなく、本発明の趣旨を逸脱しない範囲で種々変形して実施することができる。例えば、上述の実施の形態で用いている要素モデルは、そのうちのいくつかを統合して１つのモデルとしてもよい。また、本発明はポート噴射型エンジンだけでなく、ポート噴射と筒内直接噴射の両方を併用可能なエンジンにも適用可能である。また、ターボ過給機のみならず機械式過給機を備えた過給エンジンにも本発明を適用可能である。
【符号の説明】
【００４７】
２，４，６，８，１０，１２，１４，１６，１８  計算ブロック
Ｍ１  ターボ回転数モデル
Ｍ２  コンプレッサモデル
Ｍ３  インタークーラモデル
Ｍ４  スロットルモデル
Ｍ５  吸気マニホールドモデル
Ｍ６  吸気弁モデル
Ｍ７  エアフローメータモデル
【図１】
2013084318.tif 000006
【図２】
2013084318.tif 000007
【図３】
2013084318.tif 000008
【図４】
2013084318.tif 000009
【図５】
2013084318.tif 000010
【国際調査報告】
2013084318.tif 000011
2013084318.tif 000012