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