(57)【要約】【課題】ＥＧＲ制御弁を介して吸気通路に流入するＥＧＲガスの流量を、ＥＧＲガスとＥＧＲ通路間の管摩擦分を考慮することにより精度良く推定することができる内燃機関のＥＧＲガス流量推定装置を提供すること【解決手段】本装置は、機関の排気通路と吸気通路とに接続された排気還流管と、同排気還流管に介装された絞り部を有するＥＧＲ制御弁とを備えた内燃機関のＥＧＲガス流量推定装置である。本装置は、排気圧力をＰｅｘ、吸気圧をＰｂ、絞り部の有効開口面積をＡｅｇｒ、ＥＧＲガスの密度をρａ、ＥＧＲガスの比熱比をκとするとき、Ｇｅｇｒ０＝Ａｅｇｒ・（２・Ｐｅｘ・ρａ）１／２・Φ、Φ＝（（（κ／（κ−１））・（（Ｐｂ／Ｐｅｘ）２／κ−（Ｐｂ／Ｐｅｘ）（１＋１／κ）））１／２なる一般式により暫定ＥＧＲガス流量Ｇｅｇｒ０を推定する。そして、差圧（Ｐｅｘ−Ｐｂ）に応じる補正値ｄＰｇａｉｎを暫定ＥＧＲガス流量Ｇｅｇｒ０に乗じて吸気通路に流入するＥＧＲガス流量Ｇｅｇｒを推定する。【選択図】  図４
2004156458.tif 000002
【特許請求の範囲】
【請求項１】
内燃機関の排気通路と吸気通路とに接続された排気還流管と、同排気還流管に介装され同排気還流管内を流れるＥＧＲガスの流量を制御するための絞り部を有するＥＧＲ制御弁とを備えた内燃機関のＥＧＲガス流量推定装置であって、
前記ＥＧＲ制御弁の上流側におけるＥＧＲガスの圧力を上流側ガス圧力として取得する上流側ガス圧力取得手段と、
前記ＥＧＲ制御弁の下流側におけるガスの圧力を下流側ガス圧力として取得する下流側ガス圧力取得手段と、
圧縮性流体が絞り部を通過する場合の同圧縮性流体の絞り部通過流量を同絞り部直前の上流側における同圧縮性流体の圧力及び同絞り部直後の下流側における同圧縮性流体の圧力とに基いて表す一般式を用いて、前記ＥＧＲ制御弁を通過するＥＧＲガスの流量を暫定ＥＧＲガス流量として推定する手段であり、同絞り部直前の上流側における同圧縮性流体の圧力として前記取得された上流側ガス圧力を採用するとともに、同絞り部直後の下流側における同圧縮性流体の圧力として前記取得された下流側ガス圧力を採用することで同暫定ＥＧＲガス流量を推定する暫定ＥＧＲガス流量推定手段と、
前記推定された暫定ＥＧＲガス流量に含まれる前記排気還流管と前記ＥＧＲガスとの間の管摩擦に起因する誤差を補正することにより前記排気還流管から前記吸気通路に流入するＥＧＲガスの流量を推定するＥＧＲガス流量推定手段と、
を備えた内燃機関のＥＧＲガス流量推定装置。
【請求項２】
請求項１に記載の内燃機関のＥＧＲガス流量推定装置において、
前記上流側ガス圧力取得手段は、前記排気還流管が接続された前記排気通路におけるガスの圧力を前記上流側ガス圧力として取得するように構成され、
前記下流側ガス圧力取得手段は、前記排気還流管が接続された吸気通路におけるガスの圧力を前記下流側ガス圧力として取得するように構成された内燃機関のＥＧＲガス流量推定装置。
【請求項３】
請求項２に記載の内燃機関のＥＧＲガス流量推定装置において、
前記ＥＧＲガス流量推定手段は、前記取得された上流側ガス圧力と前記取得された下流側ガス圧力との差圧に基いて補正値を算出し、同補正値により前記暫定ＥＧＲガス流量に含まれる誤差を補正するように構成された内燃機関のＥＧＲガス流量推定装置。
【請求項４】
請求項３に記載の内燃機関のＥＧＲガス流量推定装置において、
前記暫定ガス流量推定手段は、前記取得された上流側ガス圧力をＰｅｘ、前記取得された下流側ガス圧力をＰｂ、前記絞り部の有効開口面積をＡｅｇｒ、前記絞り部直前の上流側におけるＥＧＲガスの密度をρａ、ＥＧＲガスの比熱比をκとするとき、
Ｇｅｇｒ０＝Ａｅｇｒ・（２・Ｐｅｘ・ρａ）１／２・Φ
Φ＝（（（κ／（κ−１））・（（Ｐｂ／Ｐｅｘ）２／κ−（Ｐｂ／Ｐｅｘ）（１＋１／κ）））１／２
なる前記一般式に基いて前記暫定ＥＧＲガス流量Ｇｅｇｒ０を推定するように構成され、
前記ＥＧＲガス流量推定手段は、前記差圧（Ｐｅｘ−Ｐｂ）が大きいほど１に向けて大きくなるように前記補正値ｄＰｇａｉｎを求め、同補正値ｄＰｇａｉｎを前記暫定ＥＧＲガス流量Ｇｅｇｒ０に乗じることにより前記吸気通路に流入するＥＧＲガスの流量Ｇｅｇｒを推定するように構成された内燃機関のＥＧＲガス流量推定装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、内燃機関の排気還流管から吸気通路に流入するＥＧＲガスの流量を推定する内燃機関のＥＧＲガス流量推定装置に関する。
【０００２】
【従来の技術】
従来より、内燃機関から排出される窒素酸化物（ＮＯｘ）の量を低減するため、同内燃機関の排ガスの一部を排気還流管を介して吸気通路に循環させるＥＧＲ装置が広く知られている。ＥＧＲ装置は、火花点火式内燃機関及びディーゼル機関の何れにも適用される。この場合、ＥＧＲガス流量は、排気還流管に介装されたＥＧＲ制御弁の開度（有効開口断面積）を可変とすることにより制御されている。
【０００３】
このようなＥＧＲ装置は、吸気通路に流入するＥＧＲガス流量を推定するとともに同吸気通路に流入する新気流量を測定し、同推定したＥＧＲガス流量と同測定した新気流量とに基いて、機関に吸入される総べてのガス（即ち、新気とＥＧＲガスとの混合ガスであり、以下「吸気」とも称呼する。）の流量に対する同ＥＧＲガス流量の比であるＥＧＲ率を求め、このＥＧＲ率に基いてＥＧＲガス流量を制御するようになっている。従って、ＥＧＲガス流量の制御等の機関制御において、吸気通路に流入するＥＧＲガス流量を精度良く推定することは極めて重要である。そこで、従来の装置は、ＥＧＲ制御弁の前後差圧と同ＥＧＲ制御弁の開度とに基いて、前記ＥＧＲガス流量を推定するようになっている（例えば、特許文献１参照。）。
【０００４】
【特許文献１】
特開２００１−２８０２０２号公報（
【０００５】及び第４頁）
【０００６】
【発明が解決しようとする課題】
上記従来の装置は、吸気通路に配設されたブーストセンサ（吸気圧センサ）により検出された圧力とＥＧＲ制御弁の上流側の適宜位置に配設された圧力センサにより検出された圧力との差圧を前記前後差圧として使用し、前記ＥＧＲガス流量を推定するようになっている。しかしながら、ＥＧＲ制御弁の直前及び直後の差圧とＥＧＲガス量の推定に使用する前記前後差圧とが、ＥＧＲガスと排気還流管との間に生じる管摩擦により大きく相違することがあり、この場合、ＥＧＲガス流量を精度良く推定できないという問題がある。
【０００７】
特に、圧縮性流体が絞り部（ＥＧＲ制御弁）を通過する場合の同圧縮性流体の絞り部通過流量を同絞り部直前の上流側における同圧縮性流体の圧力及び同絞り部直後の下流側における同圧縮性流体の圧力とに基いて表す同圧縮性流体の絞り部通過流量に関する一般式を用いることにより、ＥＧＲ制御弁を通過するＥＧＲガスの流量（従って、吸気通路に流入するＥＧＲガス流量）を推定しようとする場合、上記管摩擦の影響によるＥＧＲ流量の推定誤差が大きくなり、ＥＧＲガス流量の推定精度が悪化するという問題がある。
【０００８】
本発明は、かかる課題に対処するためになされたものであって、その目的は、ＥＧＲ制御弁を介して吸気通路に流入するＥＧＲガスの流量を、上記管摩擦分を考慮することにより精度良く推定することができる内燃機関のＥＧＲガス流量推定装置を提供することにある。
【０００９】
【本発明の概要】
本発明による内燃機関のＥＧＲガス流量推定装置は、内燃機関の排気通路と吸気通路とに接続された排気還流管と、同排気還流管に介装され同排気還流管内を流れるＥＧＲガスの流量を制御するための絞り部を有するＥＧＲ制御弁とを備えた内燃機関のＥＧＲガス流量推定装置であって、前記ＥＧＲ制御弁の上流側におけるＥＧＲガスの圧力を上流側ガス圧力として取得する上流側ガス圧力取得手段と、前記ＥＧＲ制御弁の下流側におけるガスの圧力を下流側ガス圧力として取得する下流側ガス圧力取得手段と、圧縮性流体が絞り部を通過する場合の同圧縮性流体の絞り部通過流量を同絞り部直前の上流側における同圧縮性流体の圧力及び同絞り部直後の下流側における同圧縮性流体の圧力とに基いて表す一般式を用いて、前記ＥＧＲ制御弁を通過するＥＧＲガスの流量を暫定ＥＧＲガス流量として推定する手段であり、同絞り部直前の上流側における同圧縮性流体の圧力として前記取得された上流側ガス圧力を採用するとともに、同絞り部直後の下流側における同圧縮性流体の圧力として前記取得された下流側ガス圧力を採用することで同暫定ＥＧＲガス流量を推定する暫定ＥＧＲガス流量推定手段と、前記推定された暫定ＥＧＲガス流量に含まれる前記排気還流管と前記ＥＧＲガスとの間の管摩擦に起因する誤差を補正することにより前記排気還流管から前記吸気通路に流入するＥＧＲガスの流量を推定するＥＧＲガス流量推定手段と、を備えている。
【００１０】
圧縮性流体が絞り部を通過する場合の同圧縮性流体の絞り部通過流量を表す一般式によれば、同絞り部直前の上流側における同圧縮性流体の圧力及び同絞り部直後の下流側における同圧縮性流体の圧力とを用いて同絞り部通過流量を精度良く求めることができる。
【００１１】
そこで、本装置は、ＥＧＲ制御弁が形成する絞り部の上流側の適宜位置（例えば、排気通路）におけるガスの圧力（例えば、排気圧力）を取得し、この取得したガスの圧力を前記一般式における前記絞り部直前の上流側における同圧縮性流体の圧力として採用する。また、本装置は、前記絞り部の下流側の適宜位置（例えば、吸気通路）におけるガスの圧力を取得し、この取得したガスの圧力を前記一般式における前記絞り部直後の下流側における同圧縮性流体の圧力として採用し、以って、暫定的なＥＧＲガス流量を求める。
【００１２】
そして、本装置は、前記推定された暫定ＥＧＲガス流量に含まれる前記排気還流管と前記ＥＧＲガスとの間の管摩擦に起因する誤差を補正し、前記吸気通路に流入するＥＧＲガスの流量を推定する。
【００１３】
上記一般式は、絞り部直前の上流側における同圧縮性流体の圧力及び同絞り部直後の下流側における同圧縮性流体の圧力とに基けば、精度良く同絞り部を通過するガス量を推定することを可能とする。しかし、本発明では、その一般式で使用する絞り部直前の上流側における同圧縮性流体の圧力として前記取得された上流側ガス圧力を採用するとともに、同一般式で使用する同絞り部の直後の下流側における同圧縮性流体の圧力として前記取得された下流側ガス圧力を採用する。
【００１４】
従って、本発明の装置により求められる暫定ＥＧＲガス流量には、管摩擦に起因する誤差が含まれている。そこで、上記のように、本発明による装置は前記誤差を補正する手段を備えている。この結果、本装置は、取得される絞り部上流側のガスの圧力と取得される絞り部下流側のガスの圧力との差が小さくなって、絞りの程度よりも管摩擦がＥＧＲガス量を支配的に決定する傾向が強くなる場合（即ち、絞りの変化に対するＥＧＲガス流量の変化量が小さくなるような場合）においても、吸気通路に流入するＥＧＲガス量を精度良く推定することが可能となる。
【００１５】
この場合、前記上流側ガス圧力取得手段は、前記排気還流管が接続された前記排気通路におけるガスの圧力を前記上流側ガス圧力として取得するように構成され、前記下流側ガス圧力取得手段は、前記排気還流管が接続された吸気通路におけるガスの圧力を前記下流側ガス圧力として取得するように構成されることが好適である。
【００１６】
これによれば、例えば、排気圧力センサ及び吸気圧センサの出力から、上流側ガス圧力及び下流側圧力をそれぞれ取得することができる。また、排気圧力は、所定の計算を行うことにより精度良く、且つ、（例えば、排気還流管内圧力に比べて）比較的簡単に推定可能な値であるから、排気圧力センサを省略することも可能となる。そして、本発明によれば、このように各圧力を取得する場合にも、ＥＧＲガス流量が精度良く推定される。
【００１７】
この場合、前記ＥＧＲガス流量推定手段は、前記取得された上流側ガス圧力と前記取得された下流側ガス圧力との差圧に基いて補正値を算出し、同補正値により前記暫定ＥＧＲガス流量に含まれる誤差を補正するように構成されることが好適である。
【００１８】
一般に、前記取得された上流側ガス圧力と前記取得された下流側ガス圧力との差圧ｄＰは、ガス流量に対する管摩擦の影響度合いと強い相関がある。即ち、差圧ｄＰが大きい場合、絞りの程度が絞り部を通過するガスの流量を支配的に決定すると考えられる。これに対し、差圧ｄＰが小さい場合、絞りの程度に代わり、管摩擦がガスの流量を支配的に決定すると考えられる。従って、上記構成によれば、適正に補正値が求められるので、ＥＧＲガス流量の推定精度が向上する。
【００１９】
より具体的な本発明の態様は、前記暫定ガス流量推定手段は、前記取得された上流側ガス圧力をＰｅｘ、前記取得された下流側ガス圧力をＰｂ、前記絞り部の有効開口面積をＡｅｇｒ、前記絞り部の直前の上流側におけるＥＧＲガスの密度をρａ、ＥＧＲガスの比熱比をκとするとき、
Ｇｅｇｒ０＝Ａｅｇｒ・（２・Ｐｅｘ・ρａ）１／２・Φ
Φ＝（（（κ／（κ‐１））・（（Ｐｂ／Ｐｅｘ）２／κ−（Ｐｂ／Ｐｅｘ）（１＋１／κ）））１／２
なる前記一般式に基いて前記暫定ＥＧＲガス流量Ｇｅｇｒ０を推定するように構成され、前記ＥＧＲガス流量推定手段は、前記差圧（Ｐｅｘ−Ｐｂ）が大きいほど１に向けて大きくなるように前記補正値ｄＰｇａｉｎを求め、同補正値ｄＰｇａｉｎを前記暫定ＥＧＲガス流量Ｇｅｇｒ０に乗じることにより前記吸気通路に流入するＥＧＲガスの流量Ｇｅｇｒを推定するように構成される。
【００２０】
【発明の実施の形態】
以下、本発明による内燃機関のＥＧＲガス流量推定装置及びＥＧＲ制御装置を含む内燃機関（ディーゼル機関）の制御装置の実施形態の一つについて図面を参照しつつ説明する。
【００２１】
図１は、本発明による内燃機関の制御装置を４気筒内燃機関（ディーゼル機関）１０に適用したシステム全体の概略構成を示している。このシステムは、燃料供給系統を含むエンジン本体２０、エンジン本体２０の各気筒の燃焼室にガスを導入するための吸気系統３０、エンジン本体２０からの排ガスを放出するための排気系統４０、排気還流を行うためのＥＧＲ装置５０、及び電気制御装置６０を含んでいる。
【００２２】
エンジン本体２０の各気筒の上部には燃料噴射弁２１が配設されている。各燃料噴射弁２１は、電気制御装置６０と電気的に接続されていて、同電気制御装置６０からの駆動信号（指令燃料噴射量ｑｆｉｎに応じた指令信号）により所定時間だけ開弁し、これにより燃料タンクと接続された燃料噴射用ポンプ（図示省略）から供給される高圧の燃料を噴射するようになっている。
【００２３】
吸気系統３０は、エンジン本体２０の各気筒の燃焼室にそれぞれ接続されたインテークマニホールド３１、インテークマニホールド３１の上流側集合部に接続され同インテークマニホールド３１とともに吸気通路を構成する吸気管３２（インテークマニホールド３１及び吸気管３２を「吸気管」と総称することもある。）、吸気管３２内においてスロットル弁アクチュエータ３３ａにより同吸気管３２に回動可能に保持されたスロットル弁３３、スロットル弁３３の上流において吸気管３２に介装されたインタークーラー３４、同インタークーラー３４の上流において吸気管３２に介装されたターボチャージャ３５のコンプレッサ３５ａ、及び吸気管３２の先端部に配設されたエアクリーナ３６を含んでいる。
【００２４】
排気系統４０は、エンジン本体２０の各シリンダに接続されたエキゾーストマニホールド４１、エキゾーストマニホールド４１の下流側集合部に接続された排気管４２、排気管４２に配設されたターボチャージャ３５のタービン３５ｂ、ターボチャージャ絞り弁３５ｃ、及び排気管４２に介装されたディーゼルパティキュレートフィルタ（以下、「ＤＰＮＲ」と称呼する。）４３を含んでいる。エキゾーストマニホールド４１及び排気管４２は排気通路を構成している。
【００２５】
ターボチャージャ絞り弁３５ｃは、電気制御装置６０と接続されていて、同電気制御装置６０からの駆動信号に応答してターボチャージャ３５の容量を実質的に可変とするようにタービン３５ｂに流入する排ガス通路面積を可変とする弁である。このターボチャージャ絞り弁３５ｃが閉じられてタービン３５ｂに流入する排ガス通路面積を小さくすると過給圧が増大し、逆にターボチャージャ絞り弁３５ｃが開かれてタービン３５ｂに流入する排ガス通路面積を大きくすると過給圧が低下する。
【００２６】
ＤＰＮＲ４３は、コージライト等の多孔質材料から形成されたフィルタを内蔵し、通過する排ガス中のパティキュレートを細孔表面にて捕集するフィルタである。ＤＰＮＲ４３は、担体としてのアルミナに、カリウムＫ，ナトリウムＮａ，リチウムＬｉ，セシウムＣｓのようなアルカリ金属、バリウムＢａ，カルシウムＣａのようなアルカリ土類金属、及びランタンＬａ、イットリウムＹのような希土類金属から選ばれた少なくとも一つを白金とともに担持し、ＮＯｘを吸収した後に同吸収したＮＯｘを放出して還元する吸蔵還元型ＮＯｘ触媒としても機能するようになっている。
【００２７】
ＥＧＲ装置５０は、排ガスを還流させる通路（ＥＧＲ通路）を構成する排気還流管５１と、排気還流管５１に介装されたＥＧＲ制御弁５２と、排気還流管５１に介装されたＥＧＲガス冷却装置（ＥＧＲクーラー）５３とを備えている。
【００２８】
排気還流管５１はタービン３５ｂの上流側排気通路（エキゾーストマニホールド４１）と接続された部分をＥＧＲガス（排ガス）の入口部とし、スロットル弁３３の下流側吸気通路（インテークマニホールド３１）と接続された部分をＥＧＲガスの出口部として、同入口部（排気還流管入口部）と同出口部（排気還流管出口部）とを連通し、同入口部と同出口部の間をＥＧＲガスが通過するガス流通管を構成している。
【００２９】
また、前記インテークマニホールド３１をガス流通管として考えた場合、その入口部は同インテークマニホールド３１と排気還流管５１との接続部ということになり、インテークマニホールド３１の出口部は同インテークマニホールド５１と内燃機関１０の燃焼室（シリンダ、気筒）とが接続された同燃焼室への吸気流入部（吸気弁により開閉される開口部）ということになる。
【００３０】
ＥＧＲ制御弁５２は、図２（Ａ）に示したように、コイル５２ａ、弁体５２ｂ、及び開口５２ｃ等を備えている。弁体５２ｂと開口５２ｃとで、絞りの程度が可変な絞り部が形成される。コイル５２ａは、電気制御装置６０と接続されている。コイル５２ａは、電気制御装置６０からの駆動信号（ＥＧＲ制御弁開度指令値ＳＥＧＲに応じて可変となるデューティ比を有する電圧信号）に応答して弁体５２ｂの軸方向移動量（リフト量）を変更するデューティ制御式の電磁弁である。かかるＥＧＲ制御弁５２は、図２（Ｂ）に示したように、ＥＧＲ制御弁開度指令値ＳＥＧＲに応じて開口５２ｃの有効開口面積Ａｅｇｒを変更し、エキゾーストマニホールド４１からインテークマニホールド３１に再循環される排ガス量（排気還流量、ＥＧＲガス流量）を調整するようになっている。
【００３１】
なお、本装置はデューティ制御式電磁弁をＥＧＲ制御弁５２として採用しているが、例えば、負圧を電気的に制御することで弁体のリフト量を変更する負圧駆動式弁、或いは、ステップモータ駆動式の制御弁等をＥＧＲ制御弁として採用することができる。
【００３２】
ＥＧＲガス冷却装置５３は、同装置の入口部から流入するとともに同装置の出口部から流出するＥＧＲガスの通路を内部に形成している。また、ＥＧＲガス冷却装置５３は、前記ＥＧＲガスの通路に露呈した冷却部を備えている。この冷却部内にはエンジン冷却水が冷媒として循環するようになっている。
【００３３】
電気制御装置６０は、互いにバスで接続されたＣＰＵ６１、ＣＰＵ６１が実行するプログラム、テーブル（ルックアップテーブル、マップ）、及び定数等を予め記憶したＲＯＭ６２、ＣＰＵ６１が必要に応じてデータを一時的に格納するＲＡＭ６３、電源が投入された状態でデータを格納するとともに同格納したデータを電源が遮断されている間も保持するバックアップＲＡＭ６４、並びにＡＤコンバータを含むインターフェース６５等からなるマイクロコンピュータである。
【００３４】
インターフェース６５は、吸気管３２に配置された熱線式エアフローメータ７１、インタークーラー３４とスロットル弁３３との間の吸気通路に設けられた新気温センサ（吸気温センサ）７２、スロットル弁３３の下流であって排気還流管５１が接続された部位よりも上流の吸気通路に配設された吸気圧センサ７３、エンジン回転速度センサ７４、水温センサ７５及びアクセル開度センサ７６と接続されていて、これらのセンサからの信号をＣＰＵ６１に供給するようになっている。また、インターフェース６５は、燃料噴射弁２１、スロットル弁アクチュエータ３３ａ、ターボチャージャ絞り弁３５ｃ、及びＥＧＲ制御弁５２と接続されていて、ＣＰＵ６１の指示に応じてこれらに駆動信号を送出するようになっている。
【００３５】
熱線式エアフローメータ７１は、エアクリーナ３６を介して吸気管３２内に新たに吸入された大気（即ち、新気）の質量流量（単位時間当りの吸入空気量、単位時間あたりの新気量）を計測し、同新気の質量流量に応じた信号（新気流量）Ｇａを発生するようになっている。新気温センサ７２は、エアクリーナ３６を介して吸気管３２内に吸入された新気の温度（即ち、新気温度）を検出し、同新気温度を表す信号Ｔａを発生するようになっている。
【００３６】
吸気圧センサ７３は、吸気通路内の圧力（吸気圧、過給圧）を表す信号Ｐｂを発生するようになっている。即ち、吸気圧センサ７３は、排気還流管５１が接続された吸気通路（インテークマニホールド３１）におけるガスの圧力を吸気通路内ガス圧力として取得する吸気通路内ガス圧力取得手段を構成している。また、吸気圧センサ７３は、ＥＧＲ制御弁５２の下流側におけるガスの圧力を下流側ガス圧力として取得する下流側ガス圧力取得手段でもある。
【００３７】
エンジン回転速度センサ７４は、エンジン１０の回転速度を検出し、エンジン回転速度ＮＥを表す信号を発生するとともに、各気筒の絶対クランク角度を検出し得るようになっている。水温センサ７５は、エンジン１０の冷却水温を検出し、同冷却水温を表す信号ＴＨＷを発生するようになっている。アクセル開度センサ７６は、アクセルペダルＡＰの操作量を検出し、アクセル開度（アクセル操作量）を表す信号Ａｃｃｐを発生するようになっている。
【００３８】
次に、上記のように構成された内燃機関の制御装置の作動について説明する。電気制御装置６０のＣＰＵ６１は、図３に機能ブロック図で示された各値を計算するためのプログラムを所定時間の経過毎に繰り返し実行し、実ＥＧＲ率Ｒａｃｔを計算するようになっている。以下、ブロック毎に説明を加える。なお、以下に述べる各値の幾つかは、図４に概念的に示されている。
【００３９】
＜実ＥＧＲ率Ｒａｃｔの取得＞
実ＥＧＲ率Ｒａｃｔは、エンジン１０のシリンダに吸入された単位時間あたりの実ＥＧＲガス流量Ｇｅｇｒ（実際に排気還流管５１から吸気通路内に流入するＥＧＲガスの質量流量であって、以下、「ＥＧＲガス流量Ｇｅｇｒ」と称呼する。）を同エンジン１０に吸入された単位時間あたりの全ガス量Ｇｃｙｌ（実際の全ガス質量流量であって、以下、「全ガス流量Ｇｃｙｌ」とも称呼する。）で除した値（Ｒａｃｔ＝Ｇｅｇｒ／Ｇｃｙｌ）である。ＥＧＲガス流量Ｇｅｇｒは、全ガス流量Ｇｃｙｌから、エンジン１０に吸入された単位時間あたりの新気量Ｇａａｃｔ（大気の質量、以下、「実新気流量Ｇａａｃｔ」と云う。）を減じた量である。従って、ＣＰＵ６１は、ブロックＢ１に示したように、下記（１）式に基づいて実ＥＧＲ率Ｒａｃｔを計算する。
【００４０】
【数１】
2004156458.tif 000003
【００４１】
＜実新気流量Ｇａａｃｔの取得＞
（１）式で用いられる実新気流量Ｇａａｃｔは、エアフローメータ７１が計測した計測新気流量Ｇａよりも時間的に遅れて変化するから、同計測新気流量Ｇａに一次遅れ処理を施すことにより求められる値と略等しくなる。従って、ＣＰＵ６１は、計測新気流量Ｇａに対して一次遅れ処理を施すブロックＢ２に示した下記（２）式に基づいて実新気流量Ｇａａｃｔを計算する。αは０から１までの値を有する定数である。なお、Ｇａａｃｔ（ｎ）は今回の演算による実新気流量Ｇａａｃｔ、Ｇａａｃｔ（ｎ−１）は所定時間前に実行された前回の演算による実新気流量Ｇａａｃｔ、Ｇａ（ｎ）は今回の演算タイミングにおけるエアフローメータ７１の出力に基く計測新気流量Ｇａである。
【００４２】
【数２】
2004156458.tif 000004
【００４３】
＜全ガス流量Ｇｃｙｌの取得＞
（１）式の計算を行うために更に必要となる全ガス流量Ｇｃｙｌは、気体の状態方程式からも推察されるように、スロットル弁３３の下流の吸気管内の圧力（吸気圧）Ｐｂとエンジン１０のシリンダに吸入されるガスの温度（吸入ガス温度）Ｔｂｏｕｔとに応じた値となる。このエンジン１０のシリンダに吸入されるガスの温度Ｔｂｏｕｔは、以下、「インテークマニホールド出口部ガス温度Ｔｂｏｕｔ」と称呼する。
【００４４】
実際には、全ガス流量Ｇｃｙｌはエンジン１０の気筒内に残留するガス量の影響も受ける。従って、ＣＰＵ６１は、ブロックＢ３に示したように、下記（３）式で示した実験式に基づいて実ＥＧＲ率Ｒａｃｔを計算する。（３）式において、ａ，ｂは実験により定められる適合定数であり、Ｔｂａｓｅはこれらの定数ａ，ｂを定めたときのインテークマニホールド出口部ガス温度（基準温度）である。また、（３）式で用いられる吸気圧（過給圧）Ｐｂは、吸気圧センサ７３から取得される。
【００４５】
【数３】
2004156458.tif 000005
【００４６】
＜実ガス温度Ｔｂｏｕｔの取得＞
（３）式の計算を行うためには、インテークマニホールド出口部ガス温度Ｔｂｏｕｔを取得する必要がある。ＣＰＵ６１は、ブロックＢ４に示したように、このインテークマニホールド出口部ガス温度Ｔｂｏｕｔを下記（４）式に従って計算する。
【００４７】
【数４】
2004156458.tif 000006
【００４８】
（４）式において、
Ｔｂｉｎは、図４に示したように、ＥＧＲ制御弁５２の出口側のインテークマニホールド３１の内部、即ち、ＥＧＲガスと新気とが混合された部位（以下、単に「合流部」又は「インテークマニホールド入口部」と称呼する。）の混合ガス温度であり、以下「インテークマニホールド入口部ガス温度Ｔｂｉｎ」と称呼する。
Ｔｗａｌｌｉｍは、インテークマニホールド入口部から吸気弁までの間のインテークマニホールド３１の壁温であり、以下、「インテークマニホールド壁温Ｔｗａｌｌｉｍ」と称呼する。
ηｉｍは、インテークマニホールド入口部からインテークマニホールド出口部（吸気弁により開閉される部分）までの間のインテークマニホールド３１における熱伝達率（冷却効率）であり、以下、「インテークマニホールド熱伝達率ηｉｍ」と称呼する。
【００４９】
上記（４）式は、インテークマニホールド３１の壁面とシリンダに吸入されるガスとの間での熱の授受、及び同インテークマニホールド３１の壁面と外気（インテークマニホールド３１の外側の空気の温度）Ｔａｉｒとの間での熱の授受を考慮した式であって、これらの熱の授受を右辺第２項（ηｉｍ・（Ｔｂｉｎ−Ｔｗａｌｌｉｍ））で代表させた式である。この値（ηｉｍ・（Ｔｂｉｎ−Ｔｗａｌｌｉｍ））は、インテークマニホールド３１内を吸気（新気＋ＥＧＲガス）が通過する際に、吸気の温度がどれだけ変化するかを示す変化温度対応値である。
【００５０】
ガス（吸気）とガス流通管（インテークマニホールド３１）との間の熱授受は、同ガスの入口部温度と同ガス流通管の壁温の差に強い相関（例えば、比例関係）がある。また、熱伝達率は、ガスとガス流通管壁との間の熱の授受、及びガス流通管壁と外部との間の熱の授受を適切に表すことが可能である。従って、上記構成によれば、前記熱授受を簡単で且つ精度良く推定することができ、前記変化温度対応値も精度良く推定することが可能となる。
【００５１】
ところで、（４）式によりインテークマニホールド出口部ガス温度Ｔｂｏｕｔを取得するためには、同（４）式の右辺の各値（Ｔｂｉｎ，Ｔｗａｌｌｉｍ，ηｉｍ）を取得しなければならない。以下、順に説明する。
【００５２】
＜インテークマニホールド入口部ガス温度Ｔｂｉｎの取得＞
ＣＰＵ６１は、ブロックＢ５に示したように、インテークマニホールド入口部ガス温度Ｔｂｉｎをエネルギー保存則に基く下記（５）式に従って計算する。
【００５３】
【数５】
2004156458.tif 000007
【００５４】
（５）式の右辺における各値について図４を参照しながら説明する。
Ｇａａｃｔは先に説明した実新気流量であり、前記ブロックＢ２により（２）式に従って求められる。
Ｔａは先に説明した新気温度であって新気温センサ７２により検出される。
Ｃａｉｒは新気の比熱（新気ガス比熱）であって予め与えられる定数である。
【００５５】
Ｇｅｇｒは先に説明したＥＧＲガス流量であり、後述する方法により求められる。ＴｅｇｒはＥＧＲガスと新気とが合流部にて混合する直前のＥＧＲガス温度である。即ち、温度Ｔｅｇｒは、排気還流管５１の前記吸気通路との接続部であるＥＧＲガス出口部におけるＥＧＲガスの温度であって、以下「排気還流管出口部ＥＧＲガス温度（ＥＧＲ出口部ＥＧＲガス温度）Ｔｅｇｒ」と称呼する。排気還流管出口部ＥＧＲガス温度Ｔｅｇｒは後述する方法により求められる。
ＣｅｇｒはＥＧＲガスの比熱（ＥＧＲガス比熱）であって予め与えられる定数である。
【００５６】
Ｇａｌｌは、ＥＧＲガスと新気とが混合したガスの総和量、即ち、実新気流量ＧａａｃｔとＥＧＲガス流量Ｇｅｇｒとの和であって、以下、「インテークマニホールド入口部ガス流量Ｇａｌｌ」と称呼する。
Ｃａｖｅは、ＥＧＲガスと新気とが混合したガスの比熱（混合ガス比熱）であり、予め与えられる定数である。
【００５７】
（５）式によりインテークマニホールド入口部ガス温度Ｔｂｉｎを取得するためには、排気還流管出口部ＥＧＲガス温度Ｔｅｇｒ、ＥＧＲガス流量Ｇｅｇｒ、インテークマニホールド入口部ガス流量Ｇａｌｌを取得しなければならない。従って、以下、更に順を追って説明する。
【００５８】
＜排気還流管出口部ＥＧＲガス温度Ｔｅｇｒの取得＞
ＣＰＵ６１は、ブロックＢ６に示したように、排気還流管出口部ＥＧＲガス温度Ｔｅｇｒを下記（６）式に従って計算する。ブロックＢ６は、出口部ＥＧＲガス温度推定手段を構成している。
【００５９】
【数６】
2004156458.tif 000008
【００６０】
（６）式において、
Ｔｅｘは、排気還流管５１とエキゾーストマニホールド４１との接合部近傍の同排気還流管５１入口部におけるＥＧＲガス温度（即ち、エキゾーストマニホールド４１と排気還流管５１との接合部近傍での排ガス温度）であり、以下、「排気還流管入口部ガス温度Ｔｅｘ（ＥＧＲ通路入口部ＥＧＲガス温度Ｔｅｘ）」と称呼する。
ηｅｇｒは、ＥＧＲガス冷却装置５３の冷却効率（熱伝達率）である。
ＴＨＷは、エンジン１０の冷却水温であるが、ＥＧＲ冷却装置５３の冷媒はエンジン冷却水であるから、冷却水温ＴＨＷは冷媒の温度Ｔｒｅｉｂａｉと等しい。
【００６１】
上記（６）式は、ＥＧＲガス冷却装置５３（の冷却部）と同ＥＧＲ冷却装置を通過するＥＧＲガスとの間での熱の授受を考慮した式である。即ち、（６）式の右辺第２項（ηｅｇｒ・（Ｔｅｘ‐ＴＨＷ））は、ＥＧＲガス冷却装置５３内をＥＧＲガスが通過する際に、ＥＧＲガスの温度がどれだけ変化するかを示す変化温度対応値である。
【００６２】
実際には、ＥＧＲガスは、排気還流管入口部から流入して排気還流管出口部に到達するまでの間に、排気還流管５１の壁面との間で熱の授受を行う。しかしながら、ＥＧＲガスと排気還流管５１の壁面との間の熱交換量は、ＥＧＲガスとＥＧＲガス冷却装置５３との間の熱交換量に比べて極めて小さい。従って、（６）式の右辺第２項（ηｅｇｒ・（Ｔｅｘ‐ＴＨＷ））は、ＥＧＲガスが排気還流管入口部から流入して排気還流管出口部に到達するまでの間に、同ＥＧＲガスの温度がどれだけ変化するかを示す値と実質的に等しい。
【００６３】
ところで、（６）式により排気還流管出口部ＥＧＲガス温度Ｔｅｇｒを取得するためには、上記排気還流管入口部ＥＧＲガス温度Ｔｅｘ、及びＥＧＲ冷却装置の冷却効率ηｅｇｒを取得しなければならない。以下、順に説明する。
【００６４】
＜排気還流管入口部のＥＧＲガス温度Ｔｅｘの取得＞
ＣＰＵ６１は、ブロックＢ８及びブロックＢ９示したように、排気還流管入口部のＥＧＲガス温度Ｔｅｘ（排ガス温度Ｔｅｘ）を下記（７）式に従って計算する。ブロックＢ８及びブロックＢ９は、ＥＧＲガス温度取得手段を構成している。
【００６５】
【数７】
2004156458.tif 000009
【００６６】
（７）式において、
単位時間あたりの燃料噴射量Ｇｆは、図５のブロックＢＰ３に示したように、指令燃料噴射量ｑｆｉｎ及びエンジン回転速度ＮＥに基いて求められる。例えば、Ｇｆ＝ｋＧｆ・ｑｆｉｎ・ＮＥ（  ｋＧｆは定数）である。
実新気流量Ｇａａｃｔは、上述したブロックＢ２により（２）式に基いて求められる。
過給圧Ｐｂは、吸気圧Ｐｂであり、吸気圧センサ７３から取得される。
エキゾーストマニホールド内ガス圧力Ｐｅｘ（エキゾーストマニホールド４１が構成する排気通路と排気還流管５１との接続部であるＥＧＲガス入口部におけるＥＧＲガスの圧力）は、後述する手法により求められる。
【００６７】
上記（７）式は、「排気還流管入口部ＥＧＲガス温度Ｔｅｘは、シリンダ内へ投入されるエネルギー（発熱量）と、同シリンダ内で発生した熱のガスへの伝達に大きく依存する。」という知見に基いている。シリンダ内へ投入されるエネルギーは燃料噴射量Ｇｆに強い相関がある。また、シリンダ内で発生した熱のガスへの伝達は、実新気流量Ｇａａｃｔ（実新気流量Ｇａａｃｔは発熱に寄与しない代わりに、排ガス温度を低下させるように機能する。）、又は、ガス比熱に関連した値である等量比Φに強い相関がある。そこで、（７）式では変数ＸＴｅｘに上記したような値を選択している。
【００６８】
なお、変数Ｘｔｅｘのひとつに使用されている値（過給圧Ｐｂ／エキゾーストマニホールド内ガス圧力Ｐｅｘ）は、エキゾーストマニホールド４１内における排ガスの通過のし易さ（留まり易さ）を表している。排ガスがエキゾーストマニホールド４１内に留まるほど、同排ガスとエキゾーストマニホールド４１外部との間でより多くの熱が伝達される。従って、（過給圧Ｐｂ／エキゾーストマニホールド内ガス圧力Ｐｅｘ）をパラメータとして導入することで、排気還流管入口部ＥＧＲガス温度Ｔｅｘの推定精度が高くなる。また、過給圧ＰｂはＥＧＲガス量と相関があり、且つ、ＥＧＲガス量が多ければ燃焼開始時の温度が高くなって排ガス温度Ｔｅｘ（排気還流管入口部ＥＧＲガス温度Ｔｅｘ）も高くなるので、過給圧Ｐｂをパラメータとして考慮することは、この点においても同排ガス温度Ｔｅｘの精度向上に寄与する。
【００６９】
（７）式の関数ｆＴｅｘと定数ａはエンジンの機種毎に決定される。以下は、関数ｆＴｅｘと定数ａを決定する手順の一例である。
（Ｓｔｅｐ１）  関数ｆＴｅｘと定数ａを決定すべきエンジンの運転条件を変更し、各エンジン状態量（Ｇｆ，  Ｇａａｃｔ，  Ｐｂ，  Ｐｅｘ，  Ｔｅｘ）を測定する。
（Ｓｔｅｐ２）  測定結果に基いて、変数ＸＴｅｘとＥＧＲガス温度Ｔｅｘの実測値との相関が高くなるように定数ａを決定する。なお、変数ＸＴｅｘに等量比Φを含む値を用いた場合には、定数ａの値の調整（決定）は省略される。
（Ｓｔｅｐ３）  決定された定数ａで決る変数ＸＴｅｘとＥＧＲガス温度Ｔｅｘの実測値とに基いて、関数ｆＴｅｘを決定する。
図８は、変数ＸＴｅｘとしてＧｆ・Φ・（Ｐｂ／Ｐｅｘ）を選択した場合の同変数ＸＴｅｘと実測された排気還流管入口部ＥＧＲガス温度Ｔｅｘとの関係の一例を示している。この場合、関数ｆＴｅｘは、Ｔｅｘ＝ｆＴｅｘ（ＸＴｅｘ）＝５４５．９・ＸＴｅｘ０．３４８９となった。
【００７０】
＜エキゾーストマニホールド内ガス圧力Ｐｅｘの取得＞
上記（７）式の変数ＸＴｅｘとしてエキゾーストマニホールド内ガス圧力Ｐｅｘを含む変数を使用する場合等において、同エキゾーストマニホールド内ガス圧力Ｐｅｘを求める必要がある。ＣＰＵ６１は、機能ブロック図である図５に示したように、下記（８）式に従ってエキゾーストマニホールド内ガス圧力Ｐｅｘを算出する。なお、エキゾーストマニホールド内ガス圧力（排気還流管入口部ＥＧＲガス圧力）Ｐｅｘは、エキゾーストマニホールド４１が構成する排気通路と排気還流管５１との接続部であるＥＧＲガス入口部におけるＥＧＲガスの圧力Ｐｅｘである。また、図５に示した機能ブロックは、排気通路における前記ＥＧＲガスの圧力を取得する排気還流管入口部ＥＧＲガス圧力取得手段（ＥＧＲガス入口部ＥＧＲガス圧力取得、排気圧力取得手段）、或いは、ＥＧＲ制御弁５２の上流側におけるガスの圧力を上流側ガス圧力として取得する上流側ガス圧力取得手段を構成している。
【００７１】
【数８】
2004156458.tif 000010
【００７２】
（８）式において、
燃料噴射量Ｇｆは、図５のブロックＢＰ３に示したように、指令燃料噴射量ｑｆｉｎ及びエンジン回転速度ＮＥに基いて求められる。
実新気流量Ｇａａｃｔは、上述したブロックＢ２により（２）式に従って求められる。過給圧Ｐｂは、吸気圧Ｐｂであり、吸気圧センサ７３から取得される。
可変容量ターボチャージャ開度Ａｖｎは、図５のブロックＢＰ４内に示したように、指令燃料噴射量ｑｆｉｎ及びエンジン回転速度ＮＥを引数とするテーブルにより求める値である。ＣＰＵ６１は、ターボチャージャ絞り弁３５ｃの開度がこの値ＡＶｎに基いた値となるように、同ターボチャージャ絞り弁３５ｃに対して駆動信号を供給する。また、可変容量ターボチャージャ開度Ａｖｎは、ブロックＢＰ５にて定数ａｖｎが加えられ、上記（８）式の可変容量ターボチャージャ絞り係数Ｋｖｎに変換される。なお、後述するように、目標過給圧を設定して実際の過給圧が同目標過給圧となるように値ＡＶｎを決定してもよい。
【００７３】
上記（８）式は、「エキゾーストマニホールド内ガス圧力Ｐｅｘは、シリンダ内に流入するガス量（Ｇａａｃｔ＋Ｇｆ）、可変容量ターボチャージャ絞り弁３５ｃの開度Ａｖｎ、及びターボチャージャ３５のタービン３５ｂの抵抗を表す過給圧に極めて強い相関を有する。」という知見に基いている。
【００７４】
（８）式の関数ｆＰｅｘと定数ａｖｎはエンジンの機種毎に決定される。以下は、関数ｆＰｅｘと定数ａｖｎを決定する手順の一例である。
（Ｓｔｅｐ１）  関数ｆＰｅｘと定数ａｖｎを決定すべきエンジンの運転条件を変更し、各エンジン状態量（Ｇｆ，  Ｇａａｃｔ，  Ｐｂ，  Ａｖｎ，  Ｐｅｘ）を測定する。
（Ｓｔｅｐ２）  測定結果に基いて、変数ＸＰｅｘとエキゾーストマニホールド内ガス圧力Ｐｅｘとの相関が高くなるように定数ａｖｎを決定する。
（Ｓｔｅｐ３）  決定された定数ａｖｎで決る変数ＸＰｅｘとエキゾーストマニホールド内ガス圧力Ｐｅｘの実測値とに基いて、関数ｆＰｅｘを決定する。
図９は、このようにして関数ｆＰｅｘを決定する際に使用した実測値の一例を示している。この例の場合、関数ｆＰｅｘは下記の（９）式のようになった。このように、本実施形態は、排気圧力センサを用いることなく排気圧力Ｐｅｘを求めることができ、装置のコストを低下させることができる。
【００７５】
【数９】
2004156458.tif 000011
【００７６】
以上により、（７）式の変数ＸＴｅｘを求めるための各値（Ｇｆ，Ｇａａｃｔ，Ｐｂ，Ｐｅｘ）が求められて変数ＸＴｅｘが決定される。従って、ＣＰＵ６１は、（７）式に従って計算を行い、排気還流管入口部ＥＧＲガス温度Ｔｅｘ（排ガス温度Ｔｅｘ）を求める。一方、（６）式に従って排気還流管出口部ＥＧＲガス温度Ｔｅｇｒを求めるためには、更に、ＥＧＲガス冷却装置の冷却効率ηｅｇｒを取得しなければならない。
【００７７】
＜ＥＧＲガス冷却装置の冷却効率ηｅｇｒの取得＞
ＣＰＵ６１は、図３のブロックＢ１０に示したように、ＥＧＲガス冷却装置の冷却効率ηｅｇｒを下記（１０）式に従って計算する。ブロックＢ１０は、冷却装置の冷却効率取得手段（推定手段）を構成している。
【００７８】
【数１０】
2004156458.tif 000012
【００７９】
（１０）式に示したように、ＥＧＲガス冷却装置の冷却効率ηｅｇｒを求めるためには、排気還流管入口部ＥＧＲガス温度ＴｅｘとＥＧＲガス流量Ｇｅｇｒとが必要となる。排気還流管入口部ＥＧＲガス温度Ｔｅｘは、ブロックＢ８及びブロックＢ９により上記（７）式により求められる。ＥＧＲガス流量Ｇｅｇｒは、後述するブロックＢ１２により求められる。
【００８０】
なお、（１０）式でのＥＧＲガス流量Ｇｅｇｒは、ＥＧＲガス流量Ｇｅｇｒに対応した値（ＥＧＲガス流量対応値）であればよく、例えば、排気還流管５１の所定箇所におけるＥＧＲガス流速Ｖｅｇｒで置き換えることができる。ＥＧＲ通路（排気還流管５１とＥＧＲガス冷却装置５３の形成するＥＧＲ通路）の形状は既知であるから、ＥＧＲガス流速Ｖｅｇｒに基いてＥＧＲガス流量Ｇｅｇｒを推定できるからである。ＥＧＲガス流速Ｖｅｇｒは、排気還流管５１内に流速センサを配設し、同流速センサの出力から直接取得してもよい。
【００８１】
（１０）式の関数ｆηｅｇｒはエンジンの機種毎に決定される。以下は、関数ｆηｅｇｒを決定する手順の一例である。
（Ｓｔｅｐ１）  関数ｆηｅｇｒを決定すべきエンジンの運転条件を変更し、各エンジン状態量（Ｇｅｇｒ，  Ｔｅｘ，  ηｅｇｒ）を測定する。
（Ｓｔｅｐ２）  測定結果に基いて、図１０に示したように、ηｅｇｒとＧｅｇｒ／Ｔｅｘとの関係をグラフ化する。
（Ｓｔｅｐ３）  ｓｔｅｐ２のグラフに基いて、関数ｆηｅｇｒを決定する。
【００８２】
図２０に示したように、冷却効率ηｅｇｒとＥＧＲガス流量Ｇｅｇｒとの関係は排気還流管入口部ＥＧＲガス温度Ｔｅｘが異なると変化する。これに対し、冷却効率ηｅｇｒと、ＥＧＲガス流量Ｇｅｇｒを排気還流管入口部ＥＧＲガス温度Ｔｅｘで除した値（Ｇｅｇｒ／Ｔｅｘ）と、の関係は、図１０に示したように、排気還流管入口部ＥＧＲガス温度Ｔｅｘに関わらず一義的に定まる。換言すると、冷却効率ηｅｇｒと排気還流管入口部ＥＧＲガス温度Ｔｅｘとは略反比例の関係にあることが実験的に判明したので、値（Ｇｅｇｒ／Ｔｅｘ）を変数として冷却効率ηｅｇｒを求めることで、より簡便に関数ｆηｅｇｒを求めることができる。
【００８３】
本装置は、上記関数ｆηｅｇｒを関数の形、又は、値（Ｇｅｇｒ／Ｔｅｘ）とηｅｇｒとの組み合わせからなるデータをテーブル（一次元マップ）のテーブル値としてＲＯＭ６２内に格納していて、実際に得られた値（Ｇｅｇｒ／Ｔｅｘ）と、格納されている関数又はテーブルとに基いて実際の冷却効率ηｅｇｒを求める。
【００８４】
なお、電気制御装置６０の計算能力及び／又は記憶容量に余裕がある場合、（Ｇｅｇｒ，Ｔｅｘ，ηｅｇｒ）をエンジンの運転状態を変化させながら実測し、その実測データをテーブルＭａｐηｅｇｒ（二次元マップ）としてＲＯＭ６２内に格納しておき、実際のＥＧＲガス流量Ｇｅｇｒ、実際の排気還流管入口部ＥＧＲガス温度Ｔｅｘ、及び前記格納したテーブルＭａｐηｅｇｒとから実際の冷却効率ηｅｇｒを求めるようにしてもよい。或いは、ηｅｇｒ＝ｇＴｅｘ（Ｇｅｇｒ）なる関数ｇＴｅｘを排気還流管入口部ＥＧＲガス温度Ｔｅｘ毎に求めてＲＯＭ内に格納しておき、実際の排気還流管入口部ＥＧＲガス温度Ｔｅｘに基いて記憶している複数の関数ｇＴｅｘの中から適切な関数ｇＴｅｘを選択し、その選択した関数ｇＴｅｘと実際のＥＧＲガス流量Ｇｅｇｒとから冷却効率ηｅｇｒを求めてもよい。
【００８５】
以上により、ブロックＢ６による（６）式に従う計算に必要な排気還流管入口部ＥＧＲガス温度Ｔｅｘ、ＥＧＲガス冷却装置の冷却効率ηｅｇｒ、及び冷却水温ＴＨＷ（冷媒温度Ｔｒｅｉｂａｉ）が得られるので、同（６）式により排気還流管出口部ＥＧＲガス温度Ｔｅｇｒが求められる。この段階において、（５）式の計算を行うために更に求める必要がある変数は、ＥＧＲガス流量Ｇｅｇｒ及びインテークマニホールド入口部ガス流量Ｇａｌｌである。以下、これらの求め方について説明する。
【００８６】
＜ＥＧＲガス流量Ｇｅｇｒの取得＞
前述したＥＧＲ制御弁５２は弁体５２ｂと開口部５２ｃとにより実質的に可変絞り部を形成し、同絞り部をＥＧＲガスが通過する。このような絞り部をＥＧＲガスが通過するとき、絞り部の下流のＥＧＲ流量Ｇｅｇｒ０は下記の（１１）式により求められる。（１１）式は、絞り部での流速が音速を超えない範囲で成立する圧縮性流体に関する流量についての一般的な式である。ＥＧＲ制御弁５２において、絞り部５２ｃでの流速は音速を超えないので、この（１１）式が成立する。
【００８７】
【数１１】
2004156458.tif 000013
【００８８】
（１１）式に使用される値を、図６を参照しながら説明する。
Ｐｕｐは絞り部直前の上流側圧力である。
Ｐｄｎは絞り部直後の下流側圧力である。
ρａは、絞り部上流のＥＧＲガスの密度である。
κは、ＥＧＲガスの比熱比である。
Ａｅｇｒは、絞り部（ＥＧＲ制御弁５２の弁体５２ｂと開口５２ｃとで形成される絞り部）の有効開口面積である。
【００８９】
実際には、ＣＰＵ６１は、上流側圧力Ｐｕｐとしてエキゾーストマニホールド内ガス圧Ｐｅｘを採用し、下流側圧力Ｐｄｎとして吸気圧（過給圧）Ｐｂを採用し、ＥＧＲ制御弁５２の弁体５２ｂと開口５２ｃとで形成される絞り部直後の下流におけるＥＧＲガス流量Ｇｅｇｒ０を、下記（１２）式に基いて求める。このＥＧＲガス流量Ｇｅｇｒ０を、便宜上、暫定ＥＧＲガス流量Ｇｅｇｒ０と称呼する。
【００９０】
【数１２】
2004156458.tif 000014
【００９１】
（１２）式において、
エキゾーストマニホールド内ガス圧Ｐｅｘは、図５のブロックＢＰ１により上記（８）式（実際には、（９）式）に従って求められる。
過給圧Ｐｂは、上記吸気圧センサ７３から取得される。
ＥＧＲガスの比熱比κは予め与えられている。
有効開口面積Ａｅｇｒは、ＣＰＵ６１からＥＧＲ制御弁５２への指令値であるＥＧＲ制御弁開度指令値ＳＥＧＲと、図２（Ｂ）に示したようなＲＯＭ６２内に格納されているテーブルＭａｐＡｅｇｒとに基いて求められる。この場合、ＥＧＲ制御弁開度指令値ＳＥＧＲに代えて、ＥＧＲ制御弁５２のリフト量を検出するセンサからの信号とテーブルＭａｐＡｅｇｒとに基づいて有効開口面積Ａｅｇｒを用いてもよい。
上流側密度ρａは、エキゾーストマニホールド内ガス圧Ｐｅｘと、排気還流管５１入口部におけるＥＧＲガス温度Ｔｅｘと、気体の状態方程式とに基いて求められる。実際には、ＥＧＲガスはＥＧＲ冷却装置５３により冷却されて温度が低下するが、その際に圧力も低下するので、ＥＧＲ制御弁に流入するＥＧＲガスの密度と上記のようにして求める密度とは実質的に等しい密度となる。
【００９２】
しかしながら、上記のようにして求められるＥＧＲガス流量Ｇｅｇｒ０は、ＥＧＲガスとＥＧＲ通路（排気還流管５１及びＥＧＲガス冷却装置５３が形成する管）との間に生じる管摩擦の影響により、必ずしも精度が良好でない場合がある。
【００９３】
即ち、上記（１２）式は、上記（１１）式において上流側圧力Ｐｕｐとしてエキゾーストマニホールド内ガス圧Ｐｅｘを採用し、下流側圧力Ｐｄｎとして吸気圧（過給圧）Ｐｂを採用した結果、絞り（即ち、有効開口面積Ａｅｇｒ）の変化が絞り部を通過するガスの流量を支配的に決定するような状態においては有効である。しかし、それ以外の状態、即ち、管摩擦がガス流量を支配的に決定するようになる状態では、上記（１２）式は必ずしも有効ではない。
【００９４】
絞り（絞りの程度）又は管摩擦のどちらがガス流量を支配的に決定するかは、本来的には、絞りの有効開口面積Ａｅｇｒとエキゾーストマニホールド内ガス圧Ｐｅｘと吸気圧Ｐｂの差圧ｄＰ（＝Ｐｅｘ−Ｐｂ）との相対関係により定まる筈である。しかし、一般には、差圧ｄＰ（＝Ｐｅｘ−Ｐｂ）が大きい場合、絞り部（絞りの程度）が同絞り部を通過するガスの流量を支配的に決定すると考えられる。これに対し、差圧ｄＰが小さい場合、絞りの変化に代わり、上流側圧力Ｐｕｐ及び下流側圧力Ｐｄｎとしてエキゾーストマニホールド内ガス圧Ｐｅｘ及び吸気圧Ｐｂをそれぞれ採用した結果として、管摩擦（排気還流管５１の長さに応じる。）がガスの流量を支配的に決定すると考えられる。
【００９５】
更に、上記（１１）式において上流側圧力Ｐｕｐとしてエキゾーストマニホールド内ガス圧Ｐｅｘを採用し、下流側圧力Ｐｄｎとして吸気圧（過給圧）Ｐｂを採用した結果、計算に使用される絞り部の上下差圧ｄＰは、実際の絞り部の直前直後における上下差圧より大きくなる。このため、上下差圧ｄＰが小さい領域では絞り部の直前直後における上下差圧は極めて小さく、絞りはガス流量の支配的要因とはならない。このようなことから、（１２）式による暫定ＥＧＲガス量Ｇｅｇｒ０は、実際にＥＧＲ制御弁５２を通過するＥＧＲガス量よりも大きめの値として計算されてしまう傾向にある。
【００９６】
そこで、本装置は、「エキゾーストマニホールド内ガス圧Ｐｅｘと吸気圧（過給圧）Ｐｂの差圧ｄＰ（＝Ｐｅｘ−Ｐｂ）」に応じて定まる補正係数（補正値）ｄＰｇａｉｎを導入し、この補正係数ｄｐｇａｉｎにより上記（１２）により求められるＥＧＲガス流量Ｇｅｇｒ０を補正する。
【００９７】
即ち、差圧ｄＰが減少するほど上記管摩擦のＥＧＲガス流量に及ぼす影響が大きくなるので、差圧ｄＰの減少とともに「０」に近づき、差圧ｄＰの増大とともに「１」に近づく補正係数ｄＰｇａｉｎを、上記（１２）により求められるＥＧＲガス流量Ｇｅｇｒ０に乗じることにより、排気還流管出口部ＥＧＲガス流量Ｇｅｇｒを求める。
【００９８】
具体的には、ＣＰＵ６１は、図７に示したブロックＢＧ１により、差圧ｄＰとＲＯＭ６２内に格納されている同ブロックＢＧ１内に示したテーブルＭａｐｄＰｇａｉｎとに基いて補正係数（補正値）ｄＰｇａｉｎを求める。また、ブロックＢ１２ａにより、ＥＧＲ制御弁開度指令値ＳＥＧＲと図２（Ｂ）に示したＲＯＭ６２内に格納されているテーブルＭａｐＡｅｇｒとに基いて有効開口面積Ａｅｇｒを求める。そして、図３のブロックＢ１２により、上記（１２）式及び下記（１３）式に基いて排気還流管出口部ＥＧＲガス流量（実際にＥＧＲ制御弁５２を通過するＥＧＲガス流量）Ｇｅｇｒを求める。
【００９９】
【数１３】
2004156458.tif 000015
【０１００】
なお、ブロックＢ１２の上記（１２）式を実行する部分は、圧縮性流体が絞り部を通過する場合の同圧縮性流体の絞り部通過流量を同絞り部直前の上流側における同圧縮性流体の圧力及び同絞り部直後の下流側における同圧縮性流体の圧力とに基いて表す一般式（上記（１１）式）を用いて、前記ＥＧＲ制御弁を通過するＥＧＲガスの流量を暫定ＥＧＲガス流量Ｇｅｇｒ０として推定する手段であり、同絞り部直前の上流側における同圧縮性流体の圧力として前記取得された上流側ガス圧力を採用するとともに、同絞り部の直後の下流側における同圧縮性流体の圧力として前記取得された下流側ガス圧力を採用することで同暫定ＥＧＲガス流量Ｇｅｇｒ０を推定する暫定ＥＧＲガス流量推定手段を構成している。
【０１０１】
また、ブロックＢ１２の補正係数（補正値）ｄＰｇａｉｎを求める部分、及び（１３）式を実行して暫定ＥＧＲガス流量に補正係数ｄＰｇａｉｎを乗じて同暫定ＥＧＲガス量Ｇｅｇｒ０を補正する部分は、同推定された暫定ＥＧＲガス流量Ｇｅｇｒ０に含まれる排気還流管５１とＥＧＲガスとの間の管摩擦に起因する誤差を補正することにより排気還流管５１から吸気通路に流入するＥＧＲガスの流量Ｇｅｇｒを推定するＥＧＲガス流量推定手段を構成している。
【０１０２】
＜インテークマニホールド入口部ガス流量Ｇａｌｌ＞
インテークマニホールド入口部ガス流量Ｇａｌｌは、上述したように、実新気流量ＧａａｃｔとＥＧＲガス流量Ｇｅｇｒとの和である。ＣＰＵ６１は、図３のブロックＢ１３に示したように、インテークマニホールド入口部ガス流量Ｇａｌｌを下記（１４）式に従って計算する。
【０１０３】
【数１４】
2004156458.tif 000016
【０１０４】
（１４）式における実新気流量ＧａａｃｔはブロックＢ２により上記（２）式に基いて求められる。ＥＧＲガス流量ＧｅｇｒはブロックＢ１２により上記（１２）及び（１３）式に基いて求められる。
【０１０５】
以上により、上記（５）式に従う計算に必要な各値が求められるので、ＣＰＵ６１はブロックＢ５により同（５）式に基いてインテークマニホールド入口部ガス温度Ｔｂｉｎを求める。一方、この段階において、（４）式（ブロックＢ４）によりインテークマニホールド出口部ガス温度Ｔｂｏｕｔを求めるために必要な変数は、インテークマニホールド壁温Ｔｗａｌｌｉｍ、及びインテークマニホールド熱伝達率ηｉｍである。これらは、次のようにして求められる。
【０１０６】
＜インテークマニホールド壁温Ｔｗａｌｌｉｍの取得＞
インテークマニホールド壁温Ｔｗａｌｌｉｍは、水温センサ７５が検出する冷却水温ＴＨＷと相関が強い。従って、ＣＰＵ６１は、ブロックＢ１４により、冷却水温ＴＨＷの増大とともに増大する関数ｆ１Ｔｗａｌｌｉｍを使用して下記（１５）式に従ってインテークマニホールド壁温Ｔｗａｌｌｉｍを計算する。なお、関数ｆ１Ｔｗａｌｌｉｍは、予め実験により求められ、ＲＯＭ６２内に格納されている。
【０１０７】
【数１５】
2004156458.tif 000017
【０１０８】
＜インテークマニホールド熱伝達率ηｉｍの取得＞
ＣＰＵ６１は、ブロックＢ１５及びブロックＢ１６に示したように、インテークマニホールド熱伝達率ηｉｍを下記（１６）式に従って計算する。
【０１０９】
【数１６】
2004156458.tif 000018
【０１１０】
（１６）式において、Ｖｉｍはインテークマニホールド内ガス流速である。このインテークマニホールド内ガス流速Ｖｉｍは、インテークマニホールド３１の形状が既知であるから、上記（１６）式に示したように、インテークマニホールド入口部ガス流量Ｇａｌｌに基いて求めることができる。インテークマニホールド入口部ガス流量Ｇａｌｌは、ブロックＢ１３により上記（１６）式に従って求められる。
【０１１１】
なお、インテークマニホールド内ガス流速Ｖｉｍは、インテークマニホールド３１に流速センサを配設し、同センサの出力から直接取得してもよい。一方、（１６）式の関数ｆηｉｍはインテークマニホールド内ガス流速Ｖｉｍを変数としているが、これに代えて、インテークマニホールド入口部ガス流量Ｇａｌｌを変数としてもよい。
【０１１２】
上記（１６）式は、「インテークマニホールド熱伝達率ηｉｍは、インテークマニホールド３１内のガス流速Ｖｉｍに大きく影響を受ける。」という知見に基いている。なお、（１６）式においては、冷却水温ＴＨＷをも変数としてインテークマニホールド熱伝達率ηｉｍを求めているが、冷却水温ＴＨＷを省略して単にインテークマニホールド内ガス流速Ｖｉｍの関数（ηｉｍ＝ｆηｉｍ（Ｖｉｍ））、又はインテークマニホールド入口部ガス流量Ｇａｌｌの関数（ηｉｍ＝ｆηｉｍ（Ｇａｌｌ））としてインテークマニホールド熱伝達率ηｉｍを求めて求めてもよい。
【０１１３】
関数ｆηｉｍはエンジンの機種毎に異なるので、その都度実測値との比較に基いて決定する。あるエンジンに対する実測値の例を図１１に示す。図１１の例では、下記（１７）式に示すように関数ｆηｉｍが決定される。
【０１１４】
【数１７】
2004156458.tif 000019
【０１１５】
以上により、上記（４）式に従う計算に必要な各値（Ｔｂｉｎ，ηｉｍ，Ｔｗａｌｌｉｍ）が求められるので、ＣＰＵ６１はブロックＢ４により同（４）式に基いてインテークマニホールド出口部ガス温度Ｔｂｏｕｔを求める。従って、ＣＰＵ６１はブロックＢ３により（３）式に基いてエンジン１０に吸入された単位時間あたりの全ガス流量Ｇｃｙｌを求め、その結果、ブロックＢ１により（１）式に基いて実ＥＧＲ率Ｒａｃｔを求める。
【０１１６】
次に、このように取得される各値を用いたエンジン１０の各種制御について説明する。
【０１１７】
＜燃料噴射量及び燃料噴射時期制御＞
ＣＰＵ６１は、図１２にフローチャートにより示した燃料噴射量及び燃料噴射時期を制御するためのルーチンを所定時間の経過毎に繰り返し実行するようになっている。従って、所定のタイミングになると、ＣＰＵ６１はステップ１２００から処理を開始し、ステップ１２０５に進んでアクセル開度Ａｃｃｐ、エンジン回転速度ＮＥ、及び図１３に示したテーブル（マップ）Ｍａｐｑｆｉｎから指令燃料噴射量ｑｆｉｎを求める。テーブルＭａｐｑｆｉｎは、アクセル開度Ａｃｃｐ及びエンジン回転速度ＮＥと指令燃料噴射量ｑｆｉｎとの関係を規定するテーブルであり、ＲＯＭ６２内に格納されている。
【０１１８】
次いで、ＣＰＵ６１はステップ１２１０に進み、指令燃料噴射量ｑｆｉｎ、エンジン回転速度ＮＥ、及び図１４に示したテーブルＭａｐｆｉｎｊから基本燃料噴射時期ｆｉｎｊを決定する。テーブルＭａｐｆｉｎｊは、指令燃料噴射量ｑｆｉｎ及びエンジン回転速度ＮＥと基本燃料噴射時期ｆｉｎｊとの関係を規定するテーブルであり、ＲＯＭ６２内に格納されている。
【０１１９】
その後、ＣＰＵ６１はステップ１２１５に進んで、インテークマニホールド出口部ガス温度基準値Ｔｂｏｕｔｒｅｆを、指令燃料噴射量ｑｆｉ、エンジン回転速度ＮＥ、及び図１５に示したテーブルＭａｐＴｂｏｕｔｒｅｆから決定する。テーブルＭａｐＴｂｏｕｔｒｅｆは、指令燃料噴射量ｑｆｉｎ及びエンジン回転速度ＮＥとインテークマニホールド出口部ガス温度基準値Ｔｂｏｕｔｒｅｆとの関係を規定するテーブルであり、ＲＯＭ６２内に格納されている。このインテークマニホールド出口部ガス温度基準値Ｔｂｏｕｔｒｅｆは、指令燃料噴射量ｑｆｉｎとエンジン回転速度ＮＥとの組み合わせに対して図１４に示した基本噴射時期ｆｉｎｊを決定したときのインテークマニホールド出口部ガス温度Ｔｂｏｕｔである。
【０１２０】
次いで、ＣＰＵ６１は、ステップ１２２０に進んで、噴射時期補正値Δθを、前記決定したインテークマニホールド出口部ガス温度基準値Ｔｂｏｕｔｒｅｆと図３に示したブロックＢ４により求められている実際のインテークマニホールド出口部ガス温Ｔｂｏｕｔとの差（Ｔｂｏｕｔｒｅｆ−Ｔｂｏｕｔ）と図１６に示したテーブルＭａｐΔθとから決定する。テーブルＭａｐΔθは、差（Ｔｂｏｕｔｒｅｆ−Ｔｂｏｕｔ）と噴射時期補正値Δθとの関係を規定するテーブルであり、ＲＯＭ６２内に格納されている。
【０１２１】
次いで、ＣＰＵ６１はステップ１２２５に進み、基本噴射時期ｆｉｎｊを噴射時期補正値Δθで補正して最終噴射時期ｆｉｎｊｆｉｎａｌを決定する。このように、上記ステップ１２１５〜１２２５により、インテークマニホールド出口部ガス温Ｔｂｏｕｔに応じて噴射時期が補正されることになる。この場合、図１６から明らかなように、インテークマニホールド出口部ガス温度Ｔｂｏｕｔがインテークマニホールド出口部ガス温度基準値Ｔｂｏｕｔｒｅｆより高くなるほど噴射時期補正値Δθが負の大きな値となって最終噴射時期ｆｉｎｊｆｉｎａｌが遅角側となり、インテークマニホールド出口部ガス温度Ｔｂｏｕｔがインテークマニホールド出口部ガス温度基準値Ｔｂｏｕｔｒｅｆより低くなるほど噴射時期補正値Δθは正の大きな値となって最終噴射時期ｆｉｎｊｆｉｎａｌが進角側に移行される。
【０１２２】
このように噴射時期を決定するのは、インテークマニホールド出口部ガス温度Ｔｂｏｕｔが高いときは同温度Ｔｂｏｕｔが低いときよりも燃料の着火性が良好であるから、燃料噴射時期を遅角しても着火性に問題がなく、ＮＯｘの排出量を低減することができるからである。一方、インテークマニホールド出口部ガス温度Ｔｂｏｕｔが低いときは燃料の着火性が悪化するから、燃料噴射時期を進角して着火性を確保するためである。これにより、エンジン１０の出力性能の向上及びＮＯｘ排出量の低減が達成され得る。
【０１２３】
そして、ＣＰＵ６１は、続くステップ１２３０にて現時点が上記決定された最終燃料噴射時期ｆｉｎｊｆｉｎａｌと一致しているか否かを判定し、一致している場合はステップ１２３５に進んで上記決定された指令燃料噴射量ｑｆｉｎの燃料を噴射時期が到来しているシリンダの燃料噴射弁２１から噴射し、ステップ１２９５に進んで本ルーチンを一旦終了する。また、ステップ１２３０にて「Ｎｏ」と判定される場合、直接ステップ１２９５に進んで本ルーチンを一旦終了する。以上により、燃料噴射量と燃料噴射時期の制御が達成される。
【０１２４】
＜ＥＧＲ制御＞
次に、ＥＧＲ率制御について説明する。ＣＰＵ６１は、図１７にフローチャートにより示したＥＧＲ率を制御するためのルーチンを所定時間の経過毎に繰り返し実行するようになっている。従って、所定のタイミングになると、ＣＰＵ６１はステップ１７００から処理を開始し、ステップ１７０５に進んで、その時点の指令燃料噴射量ｑｆｉｎ、その時点のエンジン回転速度ＮＥ、及び同ステップ内に示したテーブルＭａｐＯ２ｔｇｔから目標吸気酸素濃度Ｏ２ｔｇｔを決定する。テーブルＭａｐＯ２ｔｇｔは、指令燃料噴射量ｑｆｉｎ及びエンジン回転速度ＮＥと、目標吸気酸素濃度Ｏ２ｔｇｔとの関係を規定するテーブルであり、ＲＯＭ６２内に格納されている。
【０１２５】
次に、ＣＰＵ６１はステップ１７１０にて指令燃料噴射量ｑｆｉｎ及びエンジン回転速度ＮＥから単位時間当りの供給燃料量Ｑを求め、続くステップ１７１５にて空気過剰率λを同ステップ中に記載した式（λ＝ｋλ・Ｇａａｃｔ／Ｑ）により求める。ｋλは定数である。次いで、ＣＰＵ６１はステップ１７２０にて目標ＥＧＲ率Ｒｔｇｔを、上記ステップ１７０５にて決定した目標吸気酸素濃度Ｏ２ｔｇｔ、上記ステップ１７２０にて求めた空気過剰率λ、及び同ステップ１７２０内に記載した式（Ｒｔｇｔ＝λ・（ｐ・Ｏ２ｔｇｔ＋ｑ）、  ｐ及びｑは定数）に基いて求める。なお、吸気酸素濃度、ＥＧＲ率、及び空気過剰率との関係は、例えば、特開平１０−１４１１４７に詳細に開示されている。
【０１２６】
次いで、ＣＰＵ６１はステップ１７２５にて、図３に示したブロックＢ１にて求めた実ＥＧＲ率Ｒａｃｔが前記ステップ１７２０にて求めた目標ＥＧＲ率Ｒｔｇｔより大きいか否かを判定し、同ステップ１７２０にて「Ｙｅｓ」と判定される場合、ステップ１７３０に進んでＥＧＲ制御弁５２を所定の開度だけ閉じ、ＥＧＲ率を減少させてステップ１７９５に進んで本ルーチンを一旦終了する。他方、ステップ１７２５にて「Ｎｏ」と判定される場合、ＣＰＵ６１はステップ１７３５に進んでＥＧＲ制御弁５２を所定の開度だけ開き、ＥＧＲ率を増大させてステップ１７９５に進む。以上により、実際の酸素吸気濃度が目標酸素吸気濃度Ｏ２ｔｇｔに一致するようにＥＧＲ率が制御され、ＮＯｘの排出量低減及びスモークの低減が達成される。
【０１２７】
なお、上述した図１７に示したルーチンによるＥＧＲ率制御においては、目標酸素吸気濃度Ｏ２ｔｇｔを求め、この目標酸素吸気濃度Ｏ２ｔｇｔを目標ＥＧＲ率Ｒｔｇｔに変換することによりＥＧＲ率を制御していたが、指令燃料噴射量ｑｆｉｎ及びエンジン回転速度ＮＥと目標ＥＧＲ率Ｒｔｇｔとの関係を規定した図１８に示したテーブルＭａｐＲｔｇｔ、実際の指令燃料噴射量ｑｆｉｎ、及び実際のエンジン回転速度ＮＥから直接目標ＥＧＲ率Ｒｔｇｔを求めて、実際のＥＧＲ率Ｒａｃｔが同目標ＥＧＲ率Ｒｔｇｔと一致するようにＥＧＲ制御弁５２の開度を制御してもよい。
【０１２８】
＜過給圧制御＞
次に、過給圧制御について説明する。ＣＰＵ６１は、図示しない過給圧制御のためのルーチンを所定時間の経過毎に繰り返し実行し、同所定時間の経過毎に、その時点の指令燃料噴射量ｑｆｉｎ、その時点のエンジン回転速度ＮＥ、及び図１９示したテーブルＭａｐＰｂｔｇｔから目標過給圧Ｐｂｔｇｔを決定する。テーブルＭａｐＰｂｔｇｔは、指令燃料噴射量ｑｆｉｎ及びエンジン回転速度ＮＥと、目標過給圧Ｐｂｔｇｔとの関係を規定するテーブルであり、ＲＯＭ６２内に格納されている。
【０１２９】
そして、ＣＰＵ６１は、決定された目標過給圧Ｐｂｔｇｔと吸気圧センサ７３から取得される実際の過給圧Ｐｂとを比較し、同目標過給圧Ｐｂｔｇｔと同実際の過給圧Ｐｂとが一致するようにターボチャージャ絞り弁３５ｃの開度を制御する。このようにして、過給圧制御が実行される。
【０１３０】
以上、説明したように、本発明による内燃機関の制御装置の実施形態によれば、ＥＧＲガス冷却装置５３の冷却効率が、排気還流管入口部ＥＧＲガス温度とＥＧＲガス流量対応値とに基いて求められるので、結果として、排気還流管出口部ＥＧＲガス温度の推定精度が向上する。また、新気とＥＧＲガスとの混合ガス（吸気）とインテークマニホールド３１との間の熱の授受が考慮されてインテークマニホールド出口部ガス温度Ｔｂｏｕｔが推定されるので、同温度Ｔｂｏｕｔの推定精度も向上し、その結果、ＥＧＲ率を精度良く推定することができる。
【０１３１】
また、上記実施形態によれば、絞り部を通過する圧縮性流体の流量に関する一般式により求めた暫定ＥＧＲガス流量Ｇｅｇｒ０に含まれている管摩擦に起因する誤差が、補正値ｄＰｇａｉｎにより補正される。この結果、本装置は、取得される絞り部上流側のガスの圧力（Ｐｅｘ）と取得される絞り部下流側のガスの圧力（Ｐｂ）との差が小さくなって、絞りの程度よりも管摩擦がＥＧＲガス量を支配的に決定する傾向が強くなる場合においても、吸気通路に流入する（ＥＧＲ制御弁５２を通過する）ＥＧＲガス量Ｇｅｇｒを精度良く推定することが可能となる。
【０１３２】
なお、本発明は上記実施形態に限定されることはなく、本発明の範囲内において種々の変形例を採用することができる。例えば、エキゾーストマニホールド内ガス圧力Ｐｅｘは、エキゾーストマニホールド４１と排気還流管５１との接続部近傍に排気圧力センサ（排気還流管入口部ＥＧＲガス圧力取得手段）を配設し、この排気圧力センサの出力値に基いて求めるように構成してもよい。
【０１３３】
また、上記補正係数ｄＰｇａｉｎは、有効開口面積Ａｅｇｒと差圧ｄＰとの両者に基いて求めるように構成してもよい。即ち、暫定ＥＧＲ流量Ｇｅｇｒ０を有効開口面積Ａｅｇｒと差圧ｄＰとに応じて補正するように構成してもよい。
【図面の簡単な説明】
【図１】本発明の第１実施形態に係る内燃機関の制御装置を４気筒内燃機関（ディーゼル機関）に適用したシステム全体の概略構成図である。
【図２】図２（Ａ）は図１に示したＥＧＲ制御弁の概略断面図、図２（Ｂ）は、そのＥＧＲ制御弁に付与される駆動信号と有効開口面積の関係を示したグラフ（テーブル）である。
【図３】図１に示したＣＰＵが実行するプログラムの内容を表した機能ブロック図である。
【図４】図１に示したＣＰＵが計算する値を説明するための図である。
【図５】図１に示したＣＰＵが実行するプログラムの内容を表した機能ブロック図である。
【図６】絞り部におけるガス流量の一般式を説明するための図である。
【図７】図１に示したＣＰＵが実行するプログラムの内容を表した機能ブロック図である。
【図８】排気還流管入口部ＥＧＲガス温度Ｔｅｘを求めるための関数ｆＴｅｘを決定する際に使用した実測値の一例を示すグラフである。
【図９】エキゾーストマニホールド内ガス圧力Ｐｅｘを求めるための関数ｆＰｅｘを決定する際に使用した実測値の一例を示すグラフである。
【図１０】ＥＧＲガス冷却装置の冷却効率ηｅｇｒとＥＧＲガス流量Ｇｅｇｒを排気還流管入口部ＥＧＲガス温度Ｔｅｘで除した値（Ｇｅｇｒ／Ｔｅｘ）との関係を示すグラフである。
【図１１】インテークマニホールド熱伝達率ηｉｍを求めるための関数ｆηｉｍを決定する際に使用した実測値の一例を示すグラフである。
【図１２】図１に示したＣＰＵが実行するプログラムを示したフローチャートである。
【図１３】図１に示したＣＰＵが図１２に示したプログラムを実行する際に参照する指令燃料噴射量を決定するためのテーブルである。
【図１４】図１に示したＣＰＵが図１２に示したプログラムを実行する際に参照する基本噴射時期を決定するためのテーブルである。
【図１５】図１に示したＣＰＵが図１２に示したプログラムを実行する際に参照するインテークマニホールド出口部ガス温度基準値を決定するためのテーブルである。
【図１６】図１に示したＣＰＵが図１２に示したプログラムを実行する際に参照する噴射時期補正値を決定するためのテーブルである。
【図１７】図１に示したＣＰＵが実行するプログラムを示したフローチャートである。
【図１８】第１実施形態の変形例に係る内燃機関の制御装置のＣＰＵが目標ＥＧＲ率を決定する際に参照するテーブルである。
【図１９】図１に示したＣＰＵが目標過給圧を決定するために参照するテーブルである。
【図２０】冷却効率ηｅｇｒとＥＧＲガス流量Ｇｅｇｒとの関係を、排気還流管入口部ＥＧＲガス温度Ｔｅｘをパラメータとして示しグラフである。
【符号の説明】
２１…燃料噴射弁、３１…インテークマニホールド、３２…吸気管、３３…スロットル弁、３５…ターボチャージャ、４１…エキゾーストマニホールド、４２…排気管、４３…ＤＰＮＲ、５０…ＥＧＲ装置、５１…排気還流管、５２…ＥＧＲ制御弁、６０…電気制御装置。
【図１】
2004156458.tif 000020
【図２】
2004156458.tif 000021
【図３】
2004156458.tif 000022
【図４】
2004156458.tif 000023
【図５】
2004156458.tif 000024
【図６】
2004156458.tif 000025
【図７】
2004156458.tif 000026
【図８】
2004156458.tif 000027
【図９】
2004156458.tif 000028
【図１０】
2004156458.tif 000029
【図１１】
2004156458.tif 000030
【図１２】
2004156458.tif 000031
【図１３】
2004156458.tif 000032
【図１４】
2004156458.tif 000033
【図１５】
2004156458.tif 000034
【図１６】
2004156458.tif 000035
【図１７】
2004156458.tif 000036
【図１８】
2004156458.tif 000037
【図１９】
2004156458.tif 000038
【図２０】
2004156458.tif 000039