(57)【要約】【課題】ＥＧＲガスの導入がある場合においてもスロットル弁通過空気流量を正確に算出することのできるスロットル弁通過空気流量算出装置を提供する。【解決手段】燃焼室から排気通路に排出された排気ガスの少なくとも一部を吸気通路１３に再循環させる排気再循環通路２１と、該排気再循環通路を通る排気ガスの流量を調整する制御弁２２とを具備する内燃機関１におけるスロットル弁通過空気流量ｍｔを算出する装置であって、吸気管内圧力Ｐｍと、排気ガス再循環中の定常運転時における吸気管内圧力Ｐｍｅｇｒであって、スロットル弁の開度及び上記制御弁の開度に基づき算出される吸気管内圧力Ｐｍｅｇｒと、スロットル弁上流の圧力もしくは大気圧Ｐａとに基づきスロットル弁通過空気流量ｍｔを算出する。【選択図】図１
2008069790.tif 000002
【特許請求の範囲】
【請求項１】
燃焼室から排気通路に排出された排気ガスの少なくとも一部を吸気通路に再循環させる排気再循環通路と、該排気再循環通路を通る排気ガスの流量を調整する制御弁とを具備する内燃機関におけるスロットル弁通過空気流量を算出する装置であって、吸気管内圧力Ｐｍと、排気ガス再循環中の定常運転時における吸気管内圧力Ｐｍｅｇｒであって、スロットル弁の開度及び上記制御弁の開度に基づき算出される吸気管内圧力Ｐｍｅｇｒと、スロットル弁上流の圧力もしくは大気圧Ｐａとに基づきスロットル弁通過空気流量ｍｔを算出するスロットル弁通過空気流量算出装置。
【請求項２】
上記制御弁の開度をスロットル弁の開度に換算する手段を有していて、上記排気ガス再循環中の定常運転時における吸気管内圧力Ｐｍｅｇｒが、機関回転数と、上記制御弁の開度分を含めたスロットル弁の開度とに基づいて定められる、請求項１に記載のスロットル弁通過空気流量算出装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、内燃機関のスロットル弁通過空気流量算出装置に関する。
【背景技術】
【０００２】
内燃機関の燃焼室において燃焼される混合気の空燃比を最適な値とするためには、吸気弁が閉じたときに燃焼室内に充填されている吸気ガスの量（以下、「筒内充填ガス量」と称す）を正確に推定する必要がある。通常、筒内充填ガス量は、流量センサ（エアフローメータ）等の多数のセンサ、およびこれらセンサからの出力値を引数とした多数のマップから推定される。ところが、このようにマップを用いて筒内充填ガス量を推定すると、必要なマップの数およびその引数の数が多くなり、これにより筒内充填ガス量算出時のマップ検索操作が増大し、制御負荷が増大してしまう問題がある。
【０００３】
さらに、上記のようにマップに基づいて吸入ガス量を算出するためには、各内燃機関の形式毎に実際の測定に基づいて各マップを作成する、いわゆる適合作業を行わなければならないが、この適合作業における測定点はマップの引数の数に応じて増大するため、引数の数が多くなると適合を行う際の工数が多大になり、また引数の数が一つ増えるだけで、その工数は何十倍にも増加してしまう問題がある。
【０００４】
そこで、近年において、流体力学等に基づく式で表される吸入ガス量モデルを用いることにより、マップの数および引数を減らして筒内充填ガス量を算出することが検討されている。このような吸入ガス量モデルにおいては、例えば、スロットルモデル、吸気管モデル、吸気弁モデル等を構築し、これら各モデルを用いることによりスロットル弁開度、大気圧、および大気温度から、筒内充填ガス量を求めるようにしている。
【０００５】
より詳細には、スロットルモデルは、エネルギ保存則、運動量保存則、質量保存則、および状態方程式に基づいて、スロットル弁開度、大気圧、および大気温度から、スロットル弁を通過する空気流量（以下、「スロットル弁通過空気流量」と称す）を算出するようになっている。また、吸気管モデルは、エネルギ保存則、質量保存則および状態方程式に基づいて、スロットルモデルにおいて算出されたスロットル弁通過空気流量と後述する吸気弁モデルにおいて算出された筒内吸入ガス流量とから、吸気管内圧力および吸気管内温度を算出するようになっている。さらに、吸気弁モデルは、筒内充填ガス量が吸気管内圧力にほぼ比例するという経験則から導かれた式（実験式）に基づいて、吸気管モデルにおいて算出された吸気管内圧力、吸気管内温度、および大気温度から、燃焼室内に流入する吸気ガスの流量（以下、「筒内吸入ガス流量」と称す）および筒内充填ガス量を算出するようになっている。
【０００６】
このように、吸入ガス量モデルを用いて筒内充填ガス量を算出することにより、ＥＣＵに記憶させるべきマップの数が減少し、適合を行う際の工数を削減すことができる。また、各モデルに用いられる計算式は比較的単純なものであるため、マップを使用する場合に比べて制御負荷を小さくすることができる。
【０００７】
ところで、上述したように、このような吸入ガス量モデルを用いて筒内充填ガス量を算出するためには、スロットル弁通過空気流量を求める必要があり、特許文献１には、このスロットル弁通過空気流量をスロットル弁開度等の機関運転状態に基づいて求める方法が開示されている。すなわち、スロットル弁通過空気流量ｍｔは、以下の数２によって求められるとしている。
【０００８】
【数１】
2008069790.tif 000003
【０００９】
この式は、スロットル弁の開口部分を一種のオリフィスと仮定すれば、エネルギ保存則や運動量保存則等に基づいて得られるオリフィスを通過する流体の量を表す式がスロットル弁通過空気流量を表す式として利用できること、及び上述した筒内充填ガス量と吸気管内圧力との関係を示す実験式に基づいて筒内吸入ガス流量を求める場合において、定常運転時では筒内吸入ガス流量とスロットル弁通過空気流量とが等しくなること等から得られるものである。式中、ｍｔはスロットル弁通過空気流量、Ｔｈａはスロットル弁上流側温度（大気温度）、κは空気の比熱比（κ＝Ｃｐ（等圧比熱）／Ｃｖ（等容比熱））、Ｐｍは吸気管内圧力（スロットル弁より下流で吸気弁より上流の部分の圧力）、Ｐａは大気圧、ＰｍＴＡは定常運転時における吸気管内圧力をそれぞれ表している。
【００１０】
ここで、ＰｍＴＡはスロットル弁開度θｔ及び機関回転数ＮＥ等に基づいて決定することができる。ＰｍＴＡと、スロットル弁開度θｔ及び機関回転数ＮＥ等との関係は予め実験によって求められ、関数マップとして記憶しておく。また、ａ、ｂは何れも定数（適合パラメータ）であり、機関回転数ＮＥ等に基づいて設定され得る。これら定数ａあるいはｂと機関回転数ＮＥ等との関係についても予め実験により定められ、関数マップとして記憶しておく。更に、Ｔ０は、こうした実験により各定数ａ、ｂを設定した時の吸入空気の基準温度を表している。
【００１１】
この数２から明らかなように、大気温度Ｔｈａ、大気圧Ｐａ、吸気管内圧力Ｐｍを求め、スロットル弁開度θｔ及び機関回転数ＮＥ等から関数マップに基づいて定常運転時における吸気管内圧力ＰｍＴＡを求め、機関回転数ＮＥ等から関数マップに基づいて定数ａ、ｂを設定すれば、スロットル弁通過空気流量ｍｔが求まることになる。
【００１２】
ところで、従来より、燃焼温度を下げてＮＯｘ生成を低減する等の目的で、内燃機関の排気側から吸気側へ排気ガスを再循環させ、この再循環ガス（ＥＧＲガス）を新気と共に筒内へ吸入させるようにした排気還流機構（ＥＧＲ機構）付きの内燃機関が知られている。そしてこのような内燃機関では、通常、スロットル弁よりも下流であって吸気弁よりも上流である吸気管部分にＥＧＲガスが導入されるのであるが、この場合、ＥＧＲガス導入の影響で吸気管内の温度及び圧力が上昇することになる。
【００１３】
一方、上記数２は、このようなＥＧＲ機構付きの内燃機関に対応したものではないので、上記数２をＥＧＲ機構付きの内燃機関に適用した場合には、求められるスロットル弁通過空気流量ｍｔは、ＥＧＲガス導入に伴う吸気管内温度及び圧力の上昇により大きな誤差を含むことになる。
【００１４】
【特許文献１】特開２００１−４１０９５号公報
【特許文献２】特開２００２−１４７２７９号公報
【発明の開示】
【発明が解決しようとする課題】
【００１５】
本発明は、上記問題に鑑みてなされたものであり、その目的は、ＥＧＲガスの導入がある場合においてもスロットル弁通過空気流量を正確に算出することのできるスロットル弁通過空気流量算出装置を提供することである。
【課題を解決するための手段】
【００１６】
上記課題を解決するために、１番目の発明によれば、燃焼室から排気通路に排出された排気ガスの少なくとも一部を吸気通路に再循環させる排気再循環通路と、該排気再循環通路を通る排気ガスの流量を調整する制御弁とを具備する内燃機関におけるスロットル弁通過空気流量を算出する装置であって、吸気管内圧力Ｐｍと、排気ガス再循環中の定常運転時における吸気管内圧力Ｐｍｅｇｒであって、スロットル弁の開度及び上記制御弁の開度に基づき算出される吸気管内圧力Ｐｍｅｇｒと、スロットル弁上流の圧力もしくは大気圧Ｐａとに基づきスロットル弁通過空気流量ｍｔを算出するスロットル弁通過空気流量算出装置が提供される。
【００１７】
また２番目の発明によれば１番目の発明において、上記制御弁の開度をスロットル弁の開度に換算する手段を有していて、上記排気ガス再循環中の定常運転時における吸気管内圧力Ｐｍｅｇｒが、機関回転数と、上記制御弁の開度分を含めたスロットル弁の開度とに基づいて定められる。
【発明の効果】
【００１８】
ＥＧＲガスの導入がある場合においてもスロットル弁通過空気流量を正確に算出することができる。
【発明を実施するための最良の形態】
【００１９】
以下、図面を参照して本発明の実施形態について説明する。図１に概略的に示した機関本体１は筒内噴射型火花点火式内燃機関を示す。しかしながら、本発明を別の火花点火式内燃機関や圧縮自着火式内燃機関に適用してもよい。
【００２０】
図１に示したように、本発明の第１の実施形態では機関本体１はシリンダブロック２と、シリンダブロック２内で往復動するピストン３と、シリンダブロック２上に固定されたシリンダヘッド４とを具備する。ピストン３とシリンダヘッド４との間には燃焼室５が形成される。シリンダヘッド４には各気筒毎に吸気弁６と、吸気ポート７と、排気弁８と、排気ポート９とが配置される。さらに、図１に示したようにシリンダヘッド４の内壁面の中央部には点火プラグ１０が配置され、シリンダヘッド４内壁面周辺部には燃料噴射弁１１が配置される。またピストン３の頂面には燃料噴射弁１１の下方から点火プラグ１０の下方まで延びるキャビティ１２が形成されている。
【００２１】
各気筒の吸気ポート７は下流側の吸気管１３を介してサージタンク１４に連結され、サージタンク１４は上流側の吸気管１５を介してエアクリーナ１６に連結される。吸気管１５内にはステップモータ１７によって駆動されるスロットル弁１８が配置される。一方、各気筒の排気ポート９は排気管１９に連結され、この排気管１９は排気浄化装置２０に連結される。排気通路（排気ポート、排気管等）と吸気通路（吸気ポート、吸気管）とは排気再循環通路（以下、「ＥＧＲ通路」と称す）２１を介して互いに連結され、この排気再循環通路２１内には排気再循環通路を通る排気ガスの流量を調整するための制御弁（以下、「ＥＧＲ弁」と称す）２２が配置される。
【００２２】
電子制御ユニット（ＥＣＵ）３１はディジタルコンピュータからなり、双方向性バス３２を介して相互に接続されたＲＡＭ（ランダムアクセスメモリ）３３、ＲＯＭ（リードオンリメモリ）３４、ＣＰＵ（マイクロプロセッサ）３５、入力ポート３６および出力ポート３７を具備する。吸気管１３には、吸気管内の吸気ガス（新気およびＥＧＲガス）の圧力を検出するための吸気管内圧力センサ４０が設けられており、吸気管内圧力センサ４０は吸気管内圧力に比例した出力電圧を発生し、この出力電圧が対応するＡＤ変換器３８を介して入力ポート３６に入力される。
【００２３】
また、スロットル弁１８の開度を検出するためのスロットル弁開度センサ４３と、内燃機関の周囲の大気の圧力、または吸気管１５に吸入される空気の圧力（吸気圧）を検出するための大気圧センサ４４と、内燃機関の周囲の大気の温度、または吸気管１５に吸入される空気の温度（吸気温）を検出するための大気温センサ４５が設けられ、これらセンサの出力電圧は対応するＡＤ変換器３８を介して入力ポート３６に入力される。また、アクセルペダル４６にはアクセルペダル４６の踏込み量に比例した出力電圧を発生する負荷センサ４７が接続され、負荷センサ４７の出力電圧は対応するＡＤ変換器３８を介して入力ポート３６に入力される。クランク角センサ４８は例えばクランクシャフトが３０度回転する毎に出力パルスを発生し、この出力パルスが入力ポート３６に入力される。ＣＰＵ３５ではこのクランク角センサ４８の出力パルスから機関回転数が計算される。一方、出力ポート３７は対応する駆動回路３９を介して点火プラグ１０、燃料噴射弁１１、ステップモータ１７およびＥＧＲ弁２２に接続される。
【００２４】
なお、ＥＧＲ弁２２の開度については、ＥＧＲ弁２２へ送信されたステップ信号の累積値（ステップ数）に基づいてＣＰＵ３５において計算により求められる。
【００２５】
ところで、内燃機関の制御装置として、燃焼室において燃焼される混合気の空燃比を目標空燃比にするために、吸気弁が閉じたときに燃焼室内に充填されている吸気ガスの量（筒内充填ガス量）Ｍｃを推定し、推定された筒内充填ガス量Ｍｃに基づいて混合気の空燃比が目標空燃比となるように燃料噴射弁によって噴射する燃料の量（以下、「燃料噴射量」と称す）を定めるものが公知である。そしてこのような場合、内燃機関の燃焼室において燃焼される混合気の空燃比を正確に目標空燃比とするためには、筒内充填ガス量Ｍｃを正確に推定する必要がある。
【００２６】
通常、このような筒内充填ガス量Ｍｃは、流量センサ（エアフローメータ）等の多数のセンサ、およびこれらセンサからの出力値を引数とした多数のマップから推定されるが、このようにマップを用いて筒内充填ガス量Ｍｃを推定する場合、推定される筒内充填ガス量Ｍｃの値をより正確なものにするためには、必要なマップの数およびその引数の数が多くなる。このようにマップの数が多くなると、マップを保存するためのＥＣＵのＲＯＭを記憶容量の大きいものにしなければならず、内燃機関の制御装置の製造コストが高くなってしまう。さらに、各マップを作成するにはマップが搭載される車種毎に適合を行わなければならないため、マップの数およびその引数の数が多くなると適合を行う際の労力が多大となってしまう。
【００２７】
そこで、マップを用いずに様々なモデルを用いて、数値計算により筒内充填ガス量Ｍｃを算出する内燃機関の制御装置が検討されている。このような制御装置では、数値計算を多用することにより必要なマップの数を極力減らすようにしており、これにより適合作業を行う際の工数を大幅に削減しながらも、筒内充填ガス量Ｍｃを正確に算出することができる。このような制御装置のうち本願出願人により提案されたものの一つに、図２に示した吸入ガス量モデルＭ２０を搭載した制御装置がある（特願２００１−３１６３５０号）。図示した吸入ガス量モデルＭ２０はＥＧＲ機構を装備していない内燃機関に適用される最も単純なモデルであるが、モデルを用いた内燃機関制御（より詳細には筒内充填ガス量Ｍｃの算出）の一例として以下でこの吸入ガス量モデルＭ２０について説明する。
【００２８】
吸入ガス量モデルＭ２０は、図２に示したようにスロットルモデルＭ２１、吸気管モデルＭ２２、吸気弁モデルＭ２３を備える。スロットルモデルＭ２１には、スロットル弁開度センサによって検出されたスロットル弁の開度（スロットル弁開度）θｔと、大気圧センサによって検出された内燃機関周囲の大気圧（または、吸気管に吸入されるスロットル弁上流側の空気の圧力）Ｐａと、大気温センサによって検出された内燃機関周囲の大気温度（または、吸気管に吸入される空気の温度）Ｔａと、後述する吸気管モデルＭ２２において算出された吸気管内の圧力（吸気管内圧力）Ｐｍとが入力され、これら入力された各パラメータの値を後述するスロットルモデルＭ２１のモデル式に代入することで、単位時間当たりにスロットル弁を通過する空気の流量（以下、「スロットル弁通過空気流量ｍｔ」と称す）が算出される。スロットルモデルＭ２１において算出されたスロットル弁通過空気流量ｍｔは、吸気管モデルＭ２２へ入力される。
【００２９】
吸気管モデルＭ２２には、スロットルモデルＭ２１において算出されたスロットル弁通過空気流量ｍｔと、以下で詳述する単位時間当たりに燃焼室内に流入する吸気ガスの流量（以下、「筒内吸入ガス流量ｍｃ」と称す。なお、筒内吸入ガス流量ｍｃの定義については、吸気弁モデルＭ２３において詳述する）とが入力され、これら入力された各パラメータの値を後述する吸気管モデルＭ２２のモデル式に代入することで、吸気管内に存在する吸気ガスの圧力、すなわち吸気管内圧力Ｐｍと吸気管内に存在する吸気ガスの温度、すなわち吸気管内温度Ｔｍとが算出される。吸気管モデルＭ２２において算出された吸気管内圧力Ｐｍと吸気管内温度Ｔｍは共に吸気弁モデルＭ２３へ入力され、さらに吸気管内圧力ＰｍはスロットルモデルＭ２１にも入力される。
【００３０】
吸気弁モデルＭ２３には、吸気管モデルＭ２２において算出された吸気管内圧力Ｐｍおよび吸気管内温度Ｔｍの他に大気温度Ｔａが入力され、これら入力された各パラメータの値を後述する吸気弁モデルＭ２３のモデル式に代入することで、筒内吸入ガス流量ｍｃが算出される。算出された筒内吸入ガス流量ｍｃは、筒内充填ガス量Ｍｃに変換され、この筒内充填ガス量Ｍｃに基づいて燃料噴射弁からの燃料噴射量が決定される。また、吸気弁モデルＭ２３において算出された筒内吸入ガス流量ｍｃは吸気管モデルＭ２２に入力される。
【００３１】
図２から分かるように、吸入ガス量モデルＭ２０ではあるモデルにおいて算出されたパラメータの値が別のモデルへの入力値として利用されるので、吸入ガス量モデルＭ２０全体では、実際に入力される値はスロットル弁開度θｔ、大気圧Ｐａ、および大気温度Ｔａの三つのパラメータのみであり、これら三つのパラメータから筒内充填ガス量Ｍｃが算出される。
【００３２】
次に、吸入ガス量モデルＭ２０の各モデルＭ２１〜Ｍ２３について説明する。
【００３３】
スロットルモデルＭ２１では、大気圧Ｐａ、大気温度Ｔａ、吸気管内圧力Ｐｍ、スロットル弁開度θｔから、下記数４に基づいてスロットル弁通過空気流量ｍｔが算出される。ここで、数４におけるμはスロットル弁における流量係数で、スロットル弁開度θｔの関数であり、図３に示したようなマップから定まる。また、Ａｔはスロットル弁の開口断面積（以下、「スロットル開口面積」と称す）を示し、スロットル弁開度θｔの関数である。なお、これら流量係数μおよび開口断面積Ａｔをまとめたμ・Ａｔをスロットル弁開度θｔから一つのマップで求めるようにしてもよい。また、Ｒａは気体定数に関連する定数であり、気体定数Ｒを１mol当たりの空気の質量Ｍａで除算した値である（Ｒａ＝Ｒ／Ｍａ）
【００３４】
【数２】
2008069790.tif 000004
【００３５】
また、Φ（Ｐｍ／Ｐａ）は下記数５に示した関数であり、この数５におけるκは比熱比（κ＝Ｃｐ（等圧比熱）／Ｃｖ（等容比熱）であり、一定値とする）である。この関数Φ（Ｐｍ／Ｐａ）は図４に示したようなグラフに表すことができるので、このようなグラフをマップとしてＥＣＵのＲＯＭに保存し、実際には数５を用いて計算するのではなくマップからΦ（Ｐｍ／Ｐａ）の値を求めるようにしてもよい。
【００３６】
【数３】
2008069790.tif 000005
【００３７】
これらスロットルモデルＭ２１の数４および数５は、スロットル弁１８上流の気体の圧力を大気圧Ｐａ、スロットル弁１８上流の気体の温度を大気温度Ｔａ、スロットル弁１８を通過する気体の圧力を吸気管内圧力Ｐｍとして、図５に示したようなスロットル弁１８のモデルに対して、質量保存則、エネルギ保存則および運動量保存則を適用し、さらに気体の状態方程式、比熱比の定義式、およびマイヤーの関係式を利用することによって得られる。
【００３８】
吸気管モデルＭ２２では、スロットル弁通過空気流量ｍｔ、筒内吸入ガス流量ｍｃ、および大気温度Ｔａから、下記数６および数７に基づいて吸気管内圧力Ｐｍおよび吸気管内温度Ｔｍが算出される。なお、数６および数７におけるＶｍはスロットル弁から吸気弁までの吸気管等の部分（以下、「吸気管部分」と称す）１３´の容積に等しい定数である。
【００３９】
【数４】
2008069790.tif 000006
【００４０】
【数５】
2008069790.tif 000007
【００４１】
ここで、吸気管モデルＭ２２について図６を参照して説明する。吸気管部分１３´の総気体量（総吸気ガス量）をＭとすると、総気体量Ｍの時間的変化は、吸気管部分１３´に流入する気体の流量、すなわちスロットル弁通過空気流量ｍｔと、吸気管部分１３´から流出する気体の流量、すなわち筒内吸入ガス流量ｍｃとの差に等しいため、質量保存則により下記数８が得られ、この数８および気体の状態方程式（Ｐｍ・Ｖｍ＝Ｍ・Ｒａ・Ｔｍ）より、数６が得られる。
【００４２】
【数６】
2008069790.tif 000008
【００４３】
また、吸気管部分１３´の気体のエネルギＭ・Ｃｖ・Ｔｍの時間的変化量は、吸気管部分１３´に流入する気体のエネルギと吸気管部分１３´から流出する気体のエネルギとの差に等しい。このため、吸気管部分１３´に流入する気体の温度を大気温度Ｔａ、吸気管部分１３´から流出する気体の温度を吸気管内温度Ｔｍとすると、エネルギ保存則により下記数９が得られ、この数９および上記気体の状態方程式より、数７が得られる。
【００４４】
【数７】
2008069790.tif 000009
【００４５】
吸気弁モデルＭ２３では、吸気管内圧力Ｐｍ、吸気管内温度Ｔｍ、および大気温度Ｔａから、下記数１０に基づいて、筒内吸入ガス流量ｍｃが算出される。なお、数１０におけるａ、ｂは、機関回転数ＮＥから、さらに吸気弁の位相角（バルブタイミング）および作用角を変更できる可変動弁機構を備えた内燃機関の場合には吸気弁６の位相角、作用角から定まる適合パラメータである。
【００４６】
【数８】
2008069790.tif 000010
【００４７】
上述した吸気弁モデルＭ２３について図７を参照して説明する。一般に、吸気弁６が閉じたときに燃焼室５内に充填されている吸気ガスの量である筒内充填ガス量Ｍｃは、吸気弁６が閉弁するとき（吸気弁閉弁時）に確定し、吸気弁閉弁時の燃焼室５内の圧力に比例する。また、吸気弁閉弁時の燃焼室５内の圧力は吸気弁上流の気体の圧力、すなわち吸気管内圧力Ｐｍと等しいとみなすことができる。したがって、筒内充填ガス量Ｍｃは、吸気管内圧力Ｐｍに比例すると近似することができる。
【００４８】
ここで、単位時間当たりに吸気管部分１３´から流出する全吸気ガスの量を平均化したもの、または単位時間当たりに吸気管部分１３´から全ての燃焼室５に吸入される吸気ガスの量を一つの気筒の吸気行程に亘って（後述するように本実施形態ではクランク角１８０°分）平均化したものを筒内吸入ガス流量ｍｃ（以下で詳述する）とすると、筒内充填ガス量Ｍｃが吸気管内圧力Ｐｍに比例することから、筒内吸入ガス流量ｍｃも吸気管内圧力Ｐｍに比例すると考えられる。このことから、理論および経験則に基づいて、上記数１０が得られる。なお、数１０における適合パラメータａは比例係数であり、適合パラメータｂは排気弁閉弁時において燃焼室５内に残存している既燃ガス量に関連する値（以下で説明する）である。また、実際の運転では過渡時に吸気管内温度Ｔｍが大きく変化する場合があるため、これに対する補正として理論および経験則に基づいて導かれたＴａ／Ｔｍが乗算されている。
【００４９】
ここで、筒内吸入ガス流量ｍｃについて、図８を参照して内燃機関が４気筒である場合について説明する。なお、図８は横軸がクランクシャフトの回転角度、縦軸が単位時間当たりに吸気管部分１３´から燃焼室５に実際に流入する吸気ガスの量である。図８に示したように、４気筒の内燃機関では、吸気弁６が例えば１番気筒、３番気筒、４番気筒、２番気筒の順に開弁し、各気筒に対応する吸気弁６の開弁量に応じて吸気管部分１３´から各気筒の燃焼室５内へ吸気ガスが流入する。吸気管部分１３´から各気筒の燃焼室５内に流入する吸気ガスの流量の変位は図８に破線で示した通りであり、これを総合して吸気管部分１３´から全気筒の燃焼室５に流入する吸気ガスの流量は図８に実線で示した通りである。また、例えば１番気筒への筒内充填ガス量Ｍｃは図８に斜線で示した部分に相当する。
【００５０】
これに対して、実線で示した吸気管部分１３´から全ての気筒の燃焼室５に流入する吸気ガスの量を平均化したものが筒内吸入ガス流量ｍｃであり、図中に一点鎖線で示されている。そして、この一点鎖線で示した筒内吸入ガス流量ｍｃに、４気筒の場合にはクランクシャフトが１８０°（すなわち、４ストローク式内燃機関において１サイクル中にクランクシャフトが回転する角度７２０°を気筒数で割った角度）回転するのにかかる時間ΔＴ180°を乗算したものが筒内充填ガス量Ｍｃとなる。したがって、吸気弁モデルＭ２３で算出された筒内吸入ガス流量ｍｃにΔＴ180°を乗算することで、筒内充填ガス量Ｍｃを算出することができる（Ｍｃ＝ｍｃ・ΔＴ180°）。より詳細には、筒内充填ガス量Ｍｃは吸気弁閉弁時の圧力に比例することから、吸気弁閉弁時の筒内吸入ガス流量ｍｃにΔＴ180°を乗算したものが筒内充填ガス量Ｍｃとされる。なお、以上の説明からも明らかなように、数１０における値ｂにΔＴ180°を乗算し、更にＴａ／Ｔｍを乗算すると、排気弁８閉弁時において燃焼室５内に残存している既燃ガス量が得られると考えられる。
【００５１】
次に、上記吸入ガス量モデルＭ２０を内燃機関の制御装置に実装して、実際に筒内充填ガス量Ｍｃを算出する場合について説明する。筒内充填ガス量Ｍｃは吸入ガス量モデルＭ２０を用いて、上記数４、数６、数７、および数１０を解くことにより表される。この場合、ＥＣＵで処理するために、これらの式を離散化する必要がある。時刻ｔ、計算間隔Δｔを用いて数４、数６、数７、および数１０を離散化すると、それぞれ下記数１１、数１２、数１３、および数１４が得られる。なお、吸気管内温度Ｔｍ（ｔ＋Δｔ）は、数１２および数１３によってそれぞれ算出されたＰｍ／Ｔｍ（ｔ＋Δｔ）およびＰｍ（ｔ＋Δｔ）から、数１５によって算出される。
【００５２】
【数９】
2008069790.tif 000011
【００５３】
【数１０】
2008069790.tif 000012
【００５４】
【数１１】
2008069790.tif 000013
【００５５】
【数１２】
2008069790.tif 000014
【００５６】
【数１３】
2008069790.tif 000015
【００５７】
このようにして実装された吸入ガス量モデルＭ２０では、スロットルモデルＭ２１の数１１で算出された時刻ｔにおけるスロットル弁通過空気流量ｍｔ（ｔ）と、吸気弁モデルＭ２３の数１４で算出された時刻ｔにおける筒内吸入ガス流量ｍｃ（ｔ）とが、吸気管モデルＭ２２の数１２および数１３に代入され、これにより時刻ｔ＋Δｔにおける吸気管内圧力Ｐｍ（ｔ＋Δｔ）および吸気管内温度Ｔｍ（ｔ＋Δｔ）が算出される。次いで、算出されたＰｍ（ｔ＋Δｔ）およびＴｍ（ｔ＋Δｔ）は、スロットルモデルＭ２１および吸気弁モデルＭ２３の数１１および数１４に代入され、これにより時刻ｔ＋Δｔにおけるスロットル弁通過空気流量ｍｔ（ｔ＋Δｔ）および筒内吸入ガス流量ｍｃ（ｔ＋Δｔ）が算出される。そして、このような計算を繰り返すことによって、スロットル弁開度θｔ、大気圧Ｐａ、および大気温度Ｔａから、任意の時刻ｔにおける筒内吸入ガス流量ｍｃが算出され、算出された筒内吸入ガス流量ｍｃに上記時間ΔＴ180°を乗算することで、任意の時刻ｔにおける筒内充填ガス量Ｍｃが算出される。
【００５８】
なお、内燃機関の始動時には、すなわち時刻ｔ＝０においては、吸気管内圧力Ｐｍは大気圧と等しい（Ｐｍ（０）＝Ｐａ）とされ、吸気管内温度Ｔｍは大気温度と等しい（Ｔｍ（０）＝Ｔａ）とされて、各モデルＭ２１〜Ｍ２３における計算が開始される。
【００５９】
なお、上記吸入ガス量モデルＭ２０では、大気温度Ｔａおよび大気圧Ｐａが一定であるとしているが、時刻によって変化する値としてもよく、例えば、大気温度を検出するための大気温センサによって時刻ｔにおいて検出された値を大気温度Ｔａ（ｔ）、大気圧を検出するための大気圧センサによって時刻ｔにおいて検出された値を大気圧Ｐａ（ｔ）として上記数１１、数１３、および数１４に代入するようにしてもよい。
【００６０】
ところで、以上で説明した吸入ガス量モデルＭ２０は、上述したようにＥＧＲ機構を備えていない内燃機関に適用される最も単純なモデルであり、これをそのまま図１に例示したようなＥＧＲ機構を備えた内燃機関に適用すると十分な精度が得られず制御上の不都合が生じてしまう。その原因の一つは、ＥＧＲガスの導入の影響で吸気管内温度Ｔｍ及び圧力Ｐｍが上昇し、吸気管部分１３´内への新気流入量であるスロットル弁通過空気流量ｍｔがＥＧＲガスの導入がない場合とは大きく異なるためであると考えられる。
【００６１】
そこで本発明のスロットル弁通過空気流量算出装置では、ＥＧＲガスの導入がある場合においてもスロットル弁通過空気流量が正確に算出できるように、以下で説明するような方法でスロットル弁通過空気流量を算出する。これによって、モデルを用いた内燃機関制御における制御精度の向上を図ることが可能となる。
【００６２】
すなわち、本発明のスロットル弁通過空気流量算出装置においては、下記の数１６に基づいてスロットル弁通過空気流量ｍｔが算出される。
【００６３】
【数１４】
2008069790.tif 000016
【００６４】
ここで、Ｐｍｅｇｒは、その時点での機関回転数ＮＥ、スロットル弁開度θｔ及びＥＧＲ弁開度ＳＴＰを保持した場合に収束する吸気管内圧力、すなわち排気ガス再循環中の定常運転時における吸気管内圧力を示しており、機関回転数ＮＥ、スロットル弁開度θｔ及びＥＧＲ弁開度ＳＴＰに基づいて定められる（Ｐｍｅｇｒ＝ｆ１（ＮＥ，θｔ，ＳＴＰ））。ｅ、ｄ、ｒは適合パラメータを示しており、これらについては後述する。また、Φ（Ｐｍ／Ｐａ）は数５で示された関数と同様の関数であり、Φ（Ｐｍｅｇｒ／Ｐａ）は数５で示された関数のＰｍにＰｍｅｇｒを代入したものである。
【００６５】
次に、数１６の導出過程について説明する。一般にスロットル弁通過空気流量ｍｔは上述した数４によって求めることができるが、数４に基づいて定常運転時におけるスロットル弁通過空気流量ｍｔＴＡを求めると、定常運転時における吸気管内圧力はＰｍｅｇｒであるので以下の数１７のようになる。
【００６６】
【数１５】
2008069790.tif 000017
【００６７】
そして更に、数４の各辺をそれぞれ数１７の対応する各辺で除して変形すると以下の数１８が得られる。
【００６８】
【数１６】
2008069790.tif 000018
【００６９】
一方、ＥＧＲガスが導入される場合の定常運転時における筒内吸入空気流量ｍｃａは、以下の数１９で表すことができる（数１９の導出過程については後述する）。
【００７０】
【数１７】
2008069790.tif 000019
【００７１】
ここで、ｅは機関回転数ＮＥとＥＧＲ弁開度ＳＴＰに基づいて定まる適合パラメータ（ｅ＝ｆ２（ＮＥ，ＳＴＰ））を示し、ｄは機関回転数に基づいて定まる適合パラメータ（ｄ＝ｆ３（ＮＥ））を示し、ｒは機関回転数とＥＧＲ弁開度ＳＴＰに基づいて定まる適合パラメータ（ｒ＝ｆ４（ＮＥ，ＳＴＰ））を示す。そして適合パラメータｅについては、より詳細には、以下の数２０で示されるように、吸気管内圧力Ｐｍと適合パラメータｄとの大小関係によって異なる値をとる。
【００７２】
【数１８】
2008069790.tif 000020
【００７３】
すなわち、適合パラメータｅは機関回転数ＮＥとＥＧＲ弁開度ＳＴＰによって上記ｅ１またはｅ２の何れかの値をとることが定められ、吸気管内圧力Ｐｍと適合パラメータｄとの大小関係によって最終的に上記ｅ１またはｅ２のうちの何れの値をとるかが定められる。
【００７４】
ここで、数１９は定常運転時に関して示していることを考慮すると、定常運転時においては、筒内吸入空気量ｍｃａはスロットル弁通過空気流量ｍｔＴＡと等しくなり、吸気管内圧力Ｐｍは定常運転時における吸気管内圧力Ｐｍｅｇｒとなるので、数１９より以下の数２１が得られる。
【００７５】
【数１９】
2008069790.tif 000021
【００７６】
この時、上記数２０に対応して適合パラメータｅは以下の数２２のように示される。
【００７７】
【数２０】
2008069790.tif 000022
【００７８】
そしてこの数２１を数１８に代入することによって、上述の数１６が得られる。
【００７９】
本実施形態のスロットル弁通過空気流量算出装置においては、定常運転時における吸気管内圧力Ｐｍｅｇｒは、機関回転数ＮＥ、スロットル弁開度θｔ及びＥＧＲ弁開度ＳＴＰの関数（Ｐｍｅｇｒ＝ｆ１（ＮＥ，θｔ，ＳＴＰ））として予め実験等によって求められ、マップの形でＲＯＭ３４に記憶されている。同様に適合パラメータｅは、機関回転数ＮＥとＥＧＲ弁開度ＳＴＰの関数（ｅ＝ｆ２（ＮＥ，ＳＴＰ））として予め実験等によって求められると共に、Ｐｍｅｇｒと適合パラメータｄとの大小関係をも加味して一つの値が定められるようにマップの形でＲＯＭ３４に記憶されている。適合パラメータｄは、機関回転数ＮＥの関数（ｄ＝ｆ３（ＮＥ））として予め実験等によって求められ、マップの形でＲＯＭ３４に記憶されている。適合パラメータｒは、機関回転数ＮＥとＥＧＲ弁開度ＳＴＰの関数（ｒ＝ｆ４（ＮＥ，ＳＴＰ））として予め実験等によって求められ、マップの形でＲＯＭ３４に記憶されている。
【００８０】
本実施形態のスロットル弁通過空気流量算出装置では、上述の各マップから機関回転数ＮＥ、スロットル弁開度θｔ及びＥＧＲ弁開度ＳＴＰ等に基づいてＰｍｅｇｒ及び適合パラメータｅ、ｄ、ｒが定められ、吸気管内圧力センサ４０の検出値に基づいて吸気管内圧力Ｐｍが求められて、数１６に基づいてスロットル弁通過空気流量ｍｔが算出される。なお、吸気管内圧力Ｐｍは他の手段によって推定もしくは算出されてもよく、例えば、スロットル弁開度θｔに基づいて、あるいはエアフローメータを有している場合にはその検出値に基づいて算出されてもよい。
【００８１】
以上の説明から明らかなように、数１６に基づいて算出されるスロットル弁通過空気流量ｍｔは、ＥＧＲガスの導入の影響、例えばＥＧＲガスの導入に伴う吸気管内温度Ｔｍ及び圧力Ｐｍの上昇を加味したものとなる。また、数１６は過渡状態におけるスロットル弁通過空気流量ｍｔを求めるものであるので、数１６に基づいてスロットル弁通過空気流量ｍｔを求める本実施形態のスロットル弁通過空気流量算出装置によれば、ＥＧＲガスの導入がある場合における過渡状態のスロットル弁通過空気流量が正確に算出できる。
【００８２】
なお、スロットル弁１８の形状が複雑であることや脈動の影響があるために、数１６において関数Φから得られる値の精度が低下し、数１６から得られるスロットル弁通過空気流量ｍｔの精度が低下する場合が考えられるが、このような場合であっても、定常運転時におけるスロットル弁通過空気流量は正確に算出することができる。これは、定常運転時においてはＰｍ＝Ｐｍｅｇｒとなるために、数１６は以下の数２３のようになり、関数Φが式中から無くなるためである。
【００８３】
【数２１】
2008069790.tif 000023
【００８４】
次に本発明の別の実施形態のスロットル弁通過空気流量算出装置について説明する。このスロットル弁通過空気流量算出装置は、上述した実施形態のスロットル弁通過空気流量算出装置と比較して、数１６において用いられる定常運転時における吸気管内圧力Ｐｍｅｇｒの求め方のみが異なっているので、以下では主にこの点について説明する。
【００８５】
上述した実施形態のスロットル弁通過空気流量算出装置においては、数１６において用いられるＰｍｅｇｒは、予め作成されてＲＯＭ３４に記憶されたマップから求められる。そして、このマップは機関回転数ＮＥと、スロットル弁開度θｔと、ＥＧＲ弁開度ＳＴＰとに基づいてＰｍｅｇｒが定められるようになっている。
【００８６】
これに対し、本実施形態のスロットル弁通過空気流量算出装置においては、Ｐｍｅｇｒがマップに基づいて求められる点では同様であるが、Ｐｍｅｇｒを定めるためのマップの引数は三つから二つに低減されている。
【００８７】
そして、このようなマップに関しては、その引数の数によって、マップ作成時の適合工数や制御負荷等が大きく異なり、例えば引数の数が一つ増えるだけで、適合工数が何十倍にも増加してしまう場合がある。したがって、マップの引数は一般に少ない方が好ましい。しかしながら、引数の数を減らすとマップから得られるパラメータの値の正確性が低下する場合もあるため、引数の数を減らす場合にはこの点を考慮する必要がある。
【００８８】
ところで、上述した実施形態のスロットル弁通過空気流量算出装置においては、Ｐｍｅｇｒは機関回転数ＮＥと、スロットル弁開度θｔと、ＥＧＲ弁開度ＳＴＰとに基づいて定められていたが、これはＰｍｅｇｒが主にスロットル弁１８を介して吸気管部分１３´に流入する空気量と、ＥＧＲ弁２２を介して吸気官部分１３´に流入するＥＧＲガス量と、吸気弁６を介して吸気管部分１３´から燃焼室５内へ流出していくガス（空気とＥＧＲガスの混合ガス）の量とのバランスによって定まるためである。概略的には、機関回転数ＮＥが吸気弁６を介した流出ガス量に対応し、スロットル弁開度θｔがスロットル弁１８を介した流入空気量に対応し、ＥＧＲ弁開度ＳＴＰがＥＧＲ弁２２を介した流入ＥＧＲガス量に対応すると考えられる。
【００８９】
一般に、弁を通過する気体の量は、弁の開度（より詳細には、開口面積）と、弁の上流側及び下流側の圧力により決定される。上述した吸気管部分１３´への流入空気量と流入ＥＧＲガス量について考えてみると、流入空気量についてはスロットル弁開度θｔと大気圧Ｐａ及び吸気管内圧力Ｐｍによって、流入ＥＧＲガス量についてはＥＧＲ弁開度ＳＴＰと排気圧Ｐｅ及び吸気管内圧力Ｐｍによって決定される。
【００９０】
ここで、ＥＧＲガスの導入が行われるのは、通常、低負荷の運転状態であり、このような場合、排気圧Ｐｅは大気圧Ｐａとほぼ同じである。したがって、スロットル弁１８とＥＧＲ弁２２の上下流の圧力は共に同じであると考えられ、各弁における開口面積の増減は吸気管部分１３´に流入するガス（空気またはＥＧＲガス）の量の増減に対してほぼ同様の効果があると考えられる。このためＥＧＲ弁開度ＳＴＰを相当するスロットル弁開度に換算し、その開度をもとのスロットル弁開度θｔに加えることで、二つの弁を一つの弁とみなし、モデルを単純化することが可能であると考えられる。
【００９１】
本実施形態のスロットル弁通過空気流量算出装置は、このような考えにより、Ｐｍｅｇｒを求めるマップの引数を機関回転数ＮＥとＥＧＲ弁開度分を含めたスロットル弁開度θｔの二つのみとするものであり、この引数の低減に起因して、得られるパラメータの値、すなわちＰｍｅｇｒの正確性が低下することは殆どないと考えられる。
【００９２】
また、本実施形態のスロットル弁通過空気流量算出装置においては、Ｐｍｅｇｒは機関回転数ＮＥ及びスロットル弁開度θｔの関数（Ｐｍｅｇｒ＝ｆ５（ＮＥ，θｔ））として予め実験等によって求められ、マップの形でＲＯＭ３４に記憶されているが、このマップはＥＧＲ機構を備えていない内燃機関の場合の機関回転数ＮＥ及びスロットル弁開度θｔから定常運転時における吸気管内圧力ＰｍＴＡを求めるマップと同様のものとなる。
【００９３】
以下、ＥＧＲ弁開度ＳＴＰをスロットル弁開度に換算する具体的な方法について説明する。
【００９４】
第１の方法では、まず、図９に示すように、ある機関回転数ＮＥに対してＥＧＲ弁開度ＳＴＰが種々の開度である場合におけるスロットル弁開度θｔとＰｍｅｇｒとの関係が測定される。図９の例では、全閉（０）の場合を含めステップ数で表された７種類のＥＧＲ弁開度ＳＴＰに対するスロットル弁開度θｔとＰｍｅｇｒとの関係が示されている。そして、所定の基準スロットル弁開度（好ましくは、アイドリング時に相当するスロットル弁開度とする。）θｔｒにおける各ＥＧＲ弁開度ＳＴＰの場合のＰｍｅｇｒ（以下、Ｐｍｅｇｒｂと称する）が求められ、次いで、ＥＧＲ弁開度ＳＴＰを全閉とした場合に、その各Ｐｍｅｇｒｂとなる相当スロットル弁開度θｔｂがそれぞれ求められる。
【００９５】
そして、各ＥＧＲ弁開度ＳＴＰの場合について上記相当スロットル弁開度θｔｂと上記基準スロットル弁開度θｔｒとの差が求められ、図１０で示すようなＥＧＲ弁開度ＳＴＰをスロットル弁開度（すなわち、換算スロットル弁開度θｔｃ１）に換算するマップが作成される。図１０のマップは横軸がＥＧＲ弁開度ＳＴＰを示し、縦軸がＥＧＲ弁開度ＳＴＰに相当する換算スロットル弁開度θｔｃ１を示している。例えば、図１０の例においてＥＧＲ弁開度ＳＴＰが１０ステップである場合にはスロットル弁開度の２．５°分に相当し、換算スロットル弁開度θｔｃ１は２．５°であるので、もとのスロットル弁開度に２．５°を加えれば、ＥＧＲ弁開度１０ステップ分を含めたスロットル弁開度が求められる。
【００９６】
このようにして、図１０に示したような換算マップを各機関回転数ＮＥについて予め作成しておけば、スロットル弁開度θｔと、ＥＧＲ弁開度ＳＴＰと、機関回転数ＮＥが与えられた場合に、ＥＧＲ弁開度ＳＴＰをスロットル弁開度に換算し、それに基づいてＥＧＲ弁開度分を含めたスロットル弁開度が求められる。そしてそのスロットル弁開度と機関回転数ＮＥに基づいて、上述したようなＰｍｅｇｒを機関回転数ＮＥ及びスロットル弁開度θｔの関数として示したマップ、すなわちＥＧＲ機構を備えていない内燃機関の場合の定常運転時における吸気管内圧力ＰｍＴＡを求めるマップからＰｍｅｇｒを求めることが可能になる。
【００９７】
そして、このようにして求められたＰｍｅｇｒを用いて、第１の実施形態のスロットル弁通過空気流量算出装置と同様に、数１６からスロットル弁通過空気流量ｍｔが求められる。
【００９８】
次にＥＧＲ弁開度ＳＴＰをスロットル弁開度に換算する第２の方法について説明する。この方法は、基本的に上述した第１の方法と同様であるが、スロットル弁開度θｔとスロットル弁開口面積Ａｔとの間の非線形な関係を考慮している点で異なっており、これによりＥＧＲ弁開度ＳＴＰをスロットル弁開度に換算する精度が向上される。
【００９９】
すなわち、通常、スロットル弁では、スロットル弁開度θｔとスロットル弁開口面積Ａｔとの間に非線形な関係があり、例えば、図１１に示すように開度の小さいところにおいては開度の中程度以降のところに比べ、同じ開度増加に対する開口面積Ａｔの増加が少なくなっている。このような場合にＥＧＲ弁開度分を、その時のスロットル弁開度にかかわらず同様にスロットル弁開度に換算すると、その換算スロットル弁開度では意図するスロットル弁での開口面積Ａｔの増加が得られない場合がある。
【０１００】
より具体的には、例えば、上述した図１０のマップを使用する場合、このマップの作成において基準とした上記基準スロットル弁開度θｔｒよりもその時のスロットル弁開度が小さい時には図１０のマップで得られる換算スロットル弁開度θｔｃ１をその時のスロットル弁開度に加えたスロットル弁開度としても意図する開口面積には足りず、逆に上記基準スロットル弁開度θｔｒよりもその時のスロットル弁開度が大きい時には図１０のマップで得られる換算スロットル弁開度θｔｃ１をその時のスロットル弁開度に加えたスロットル弁開度とすると意図する開口面積よりも大きくなってしまう。
【０１０１】
この方法では、以下で説明するような方法で図１０のマップで得られる換算スロットル弁開度θｔｃ１を補正することにより、スロットル弁開度θｔとスロットル弁開口面積Ａｔとの間の非線形な関係を考慮し、ＥＧＲ弁開度ＳＴＰを精度良くスロットル弁開度に換算する。
【０１０２】
すなわち、図１２は図１１のスロットル弁開度θｔの小さい範囲について拡大したものであるが、まずここで図１０のマップの作成において基準とした上記基準スロットル弁開度θｔｒにおける傾きＳｂを求める。そして、この傾きＳｂを図１１のスロットル弁開度θｔに対するスロットル開口面積Ａｔを示す曲線の傾きＳで除して各スロットル弁開度θｔに対するＳｂ／Ｓを求めマップとする。こうして求めたマップが図１３の実線で示されるものであり、この値が原則として各スロットル弁開度θｔに対する補正係数Ｃｃとなる。
【０１０３】
ここで、図１１に示されているようにスロットル弁開度θｔの極めて大きな部分では、スロットル弁開度θｔが大きくなってもスロットル開口面積Ａｔはそれ以上大きくならない部分があり、この部分のＳｂ／Ｓを求めると、Ｓが零になるために無限大になってしまう。一方、このようなスロットル弁開度θｔが極めて大きい部分においては、スロットル弁開度θｔが変わってもＰｍｅｇｒは殆ど変化せず、また、ＥＧＲ弁開度ＳＴＰが変化してもＰｍｅｇｒは殆ど変化しないことがわかっているので、実質的にはＥＧＲ弁開度分をスロットル弁開度に換算して加える必要のない部分である。このようなことから、本実施形態においては、スロットル弁開度θｔの大きな部分における補正係数Ｃｃを図１３に点線で示したように単純なＳｂ／Ｓとは異なる値とし、制御上の不都合を回避すると共に制御の連続性を維持し、且つ、実際の現象に合致するようにしている。つまり、図１３において補正係数Ｃｃは、スロットル弁開度θｔの大きい部分以外では実線で、スロットル弁開度θｔの大きい部分では点線で示された値となる。
【０１０４】
このような図１３のマップから得られる補正係数Ｃｃを図１０から得られる換算スロットル弁開度θｔｃ１に乗ずることにより、スロットル弁開度θｔとスロットル弁開口面積Ａｔとの間の非線形な関係を考慮した換算スロットル弁開度θｔｃ２が得られる。
【０１０５】
そして、上述した実施形態の場合と同様に、この換算スロットル弁開度θｔｃ２に基づいてＥＧＲ弁開度分を含めたスロットル弁開度が求められ、そのスロットル弁開度と機関回転数ＮＥに基づいて、上述したようなＰｍｅｇｒを機関回転数ＮＥ及びスロットル弁開度θｔの関数として示したマップからＰｍｅｇｒを求めることが可能になる。そして、このようにして求められたＰｍｅｇｒを用いて、第１の実施形態のスロットル弁通過空気流量算出装置と同様に、数１６からスロットル弁通過空気流量ｍｔが求められる。特に本実施形態においては、換算スロットル弁開度θｔｃ２が精度良く求められるので、最終的に算出されるスロットル弁通過空気流量ｍｔの精度も向上する。
【０１０６】
最後に、上述した数１９の導出過程について説明する。
【０１０７】
ＥＧＲ弁２２が開弁されしたがってＥＧＲガスが供給されているときには、各気筒の筒内に新気とＥＧＲガスとの混合ガスが吸入される。したがって、吸気弁７が開弁し次いで閉弁したときに各気筒の筒内に充填されているＥＧＲガスの量を筒内充填ＥＧＲガス量Ｍｃｅｇｒと称すると、筒内充填ガス量Ｍｃは筒内充填新気量Ｍｃａと筒内充填ＥＧＲガス量Ｍｃｅｇｒとの和で表されることになる（Ｍｃ＝Ｍｃａ＋Ｍｃｅｇｒ）。また、同様に、筒内吸入ガス流量ｍｃは、そのうちの新気分である筒内吸入空気流量ｍｃａとＥＧＲガス分である筒内吸入ＥＧＲガス流量ｍｃｅｇｒとの和で表される（ｍｃ＝ｍｃａ＋ｍｃｅｇｒ）。
【０１０８】
ところで、上述したように、筒内充填ガス量Ｍｃは吸気弁７が閉弁したときの吸気管内圧力Ｐｍに比例し、したがってその一次式で表すことができる。また、同様に筒内吸入ガス流量ｍｃも吸気管内圧力Ｐｍに比例し、したがってその一次式で表すことができる。そしてＥＧＲガスが供給されていないときには筒内に新気のみが充填されるので、筒内充填新気量Ｍｃａ及び筒内吸入空気流量ｍｃａを吸気管内圧力Ｐｍの一次式で表すことができる。
【０１０９】
ところが、ＥＧＲガスが供給されているときには状況が全く異なり、筒内には新気だけでなくＥＧＲガスも充填される。このため、従来では、筒内充填新気量Ｍｃａや筒内吸入空気流量ｍｃａを吸気管内圧力Ｐｍの一次式で表すことは到底できないと考えられていたのである。
【０１１０】
ここで、筒内吸入ＥＧＲガス流量ｍｃｅｇｒを吸気管内圧力Ｐｍの一次式で表すことができるならば、筒内吸入ガス流量ｍｃを吸気管内圧力Ｐｍの一次式で表すことができること、筒内吸入ガス流量ｍｃが筒内吸入空気流量ｍｃａと筒内吸入ＥＧＲガス流量ｍｃｅｇｒとの和であることを考えれば、筒内吸入空気流量ｍｃａを吸気管内圧力Ｐｍの一次式で表すことができる。
【０１１１】
しかしながら、従来では、筒内吸入ＥＧＲガス流量ｍｃｅｇｒも吸気管内圧力Ｐｍの一次式で表すことができないと考えられていたのである。このことを図１４を参照しながら説明する。
【０１１２】
まず、図１４（Ａ）に示されるように、ＥＧＲ弁２２上流のＥＧＲガス圧力が排気管１９内の排気圧Ｐｅであり、ＥＧＲ制御弁上流のＥＧＲガス温度が排気管１９内の排気温Ｔｅであり、ＥＧＲ弁２２を通過するＥＧＲガスの圧力が吸気管内圧力Ｐｍであると考えると、ＥＧＲ弁２２を通過するＥＧＲガスの流量であるＥＧＲ弁通過ガス流量ｍｅｇｒは、上述した数４の場合と同様にして以下の数２４により表すことができる。
【０１１３】
【数２２】
2008069790.tif 000024
【０１１４】
ここで、μはＥＧＲ弁２２における流量係数を、ＡｅはＥＧＲ弁２２における開口断面積を、Ｒｅは気体定数Ｒに関する定数を、Φ（Ｐｍ／Ｐｅ）はＰｍ／Ｐｅの関数を、それぞれ表している。なお、流量係数μおよび開口断面積ＡｅはＥＧＲ弁２２の開度ＳＴＰによって定まる値であり、定数Ｒｅは気体定数Ｒを１ｍｏｌ当たりの排気ガスないしＥＧＲガスの質量Ｍｅで除算した値である（Ｒｅ＝Ｒ／Ｍｅ）。
【０１１５】
また、関数Φ（Ｐｍ／Ｐｅ）は、上述した数５と同様に、比熱比κを用いて以下の数２５により表される。
【０１１６】
【数２３】
2008069790.tif 000025
【０１１７】
ここで、計算を簡単にするために排気圧Ｐｅが大気圧Ｐａであるとすると、数２４により表されるＥＧＲ弁通過ガス流量ｍｅｇｒは図１４（Ｂ）のようになる。すなわち、ＥＧＲ弁通過ガス流量ｍｅｇｒは吸気管内圧力Ｐｍが小さいときにはほぼ一定に維持され、吸気管内圧力Ｐｍが高くなると図１４（Ｂ）においてＮＲで示されるように吸気管内圧力Ｐｍに対し非線形性を示しながら大気圧Ｐａに向けて減少する。なお、この非線形性部分ＮＲは数２４のうちＰｅ／√Ｔｅの部分及び関数Φ（Ｐｍ／Ｐｅ）によるものである。
【０１１８】
したがって、ＥＧＲ弁通過ガス流量ｍｅｇｒ、とりわけ非線形性部分ＮＲを吸気管内圧力Ｐｍの一次式により表すことはできないものと考えられていたのである。
【０１１９】
ところが、本願発明者らによれば、ＥＧＲ弁通過ガス流量ｍｅｇｒを吸気管内圧力Ｐｍの二つの一次式で表すことができ、したがって筒内吸入空気量ｍｃａを吸気管内圧力Ｐｍの二つの一次式で表すことができることが判明したのである。
【０１２０】
すなわち、まず、図１５に示されるように、排気温Ｔｅは吸気管内圧力Ｐｍの増大に対し、排気圧Ｐｅが増大するよりも大幅に増大し、その結果Ｐｅ／√Ｔｅを吸気管内圧力Ｐｍの一次式で表すことができるのである。
【０１２１】
また、関数Φ（Ｐｍ／Ｐｅ）も吸気管内圧力Ｐｍの一次式で表すことができる。これについて図１６を参照して説明する。排気圧Ｐｅが一定の大気圧Ｐａに維持されるのではなく、吸気管内圧力Ｐｍに応じて変動することを考慮すると、図１６（Ａ）に示されるように、吸気管内圧力ＰｍがＰｍ１のときの関数Φ（Ｐｍ／Ｐｅ）は大気圧Ｐａに収束する曲線Ｃａ上にあるのではなく、排気圧Ｐｅ１に収束する曲線Ｃ１上にあり、これがプロット（○）で表されている。同様に、Ｐｍ＝Ｐｍ２（＞Ｐｍ１）のときのΦ（Ｐｍ／Ｐｅ）は排気圧Ｐｅ２（＞Ｐｅ１）に収束する曲線Ｃ２上にあり、Ｐｍ＝Ｐｍ３（＞Ｐｍ２）のときのΦ（Ｐｍ／Ｐｅ）は排気圧Ｐｅ３（＞Ｐｅ２）に収束する曲線Ｃ３上にある。
【０１２２】
このようにして得られるプロットは図１６（Ｂ）に示されるように、直線Ｌ２で結ぶことができる。したがって、関数Φ（Ｐｍ／Ｐｅ）は吸気管内圧力Ｐｍが小さいときには直線Ｌ１に相当する吸気管内圧力Ｐｍの一次式により、吸気管内圧力Ｐｍが大きいときには直線Ｌ２に相当する吸気管内圧力Ｐｍの一次式により表すことができ、斯くして吸気管内圧力Ｐｍの二つの一次式で表すことができることになる。すなわち、ＥＧＲ弁通過ガス流量ｍｅｇｒを吸気管内圧力Ｐｍの二つの一次式で表すことができるのである。
【０１２３】
そうすると、定常運転時にはＥＧＲ弁通過ガス流量ｍｅｇｒが筒内吸入ＥＧＲガス流量ｍｃｅｇｒに等しいことを考慮すれば、定常運転時の筒内吸入ＥＧＲガス流量ｍｃｅｇｒを吸気管内圧力Ｐｍの一次式で表すことができるということになる。
【０１２４】
したがって、定常運転時の筒内吸入空気流量ｍｃａを吸気管内圧力Ｐｍの二つの一次式で表すことができるということになり、これが基本的な考え方となる。
【０１２５】
図１７には、機関回転数ＮＥ及びＥＧＲ弁開度ＳＴＰがそれぞれ一定であるときの、定常運転時の筒内吸入空気流量ｍｃａを表す吸気管内圧力Ｐｍの二つの一次式の一例が示されている。図１７に示されるように、筒内吸入空気流量ｍｃａは、勾配が互いに異なりかつ接続点ＣＰにおいて連続している、吸気管内圧力Ｐｍの二つの一次式により表される。すなわち、吸気管内圧力Ｐｍが小さいときには勾配ｅ１の一次式により、吸気管内圧力Ｐｍが高いときには勾配ｅ２の一次式により、筒内吸入空気流量ｍｃａが表される。
【０１２６】
ここで、一次式の勾配をｅとし、接続点ＣＰにおける吸気管内圧力及び筒内吸入空気流量をそれぞれｄ、ｒとすると、これら二つの一次式は以下の数２６により表すことができる。
【０１２７】
【数２４】
2008069790.tif 000026
【０１２８】
そして、これらをひとまとめにして表すと上述の数１９及び数２０が得られる。
【図面の簡単な説明】
【０１２９】
【図１】本発明の一実施形態のスロットル弁通過空気流量算出装置を備えた内燃機関全体を示す説明図である。
【図２】ＥＧＲ機構を備えていない内燃機関に適用可能な吸入ガス量モデルを示す図である。
【図３】スロットル弁開度と流量係数との関係を示す図である。
【図４】関数Φ（Ｐｍ／Ｐａ）を示す図である。
【図５】スロットルモデルの基本概念を示す図である。
【図６】ＥＧＲ機構を備えていない内燃機関に適用可能な吸気管モデルの基本概念を示す図である。
【図７】吸気弁モデルの基本概念を示す図である。
【図８】筒内充填ガス量及び筒内吸入ガス流量の定義に関する図である。
【図９】スロットル弁開度と定常運転時における吸気管内圧力との関係をＥＧＲ弁開度をパラメータとして示した図である。
【図１０】ＥＧＲ弁開度と換算スロットル弁開度との関係を示す図である。
【図１１】スロットル弁開度とスロットル開口面積との関係を示す図である。
【図１２】図１１のスロットル弁開度の小さい範囲について拡大した図である。
【図１３】スロットル弁開度と図１０から得られる換算スロットル弁開度を補正する補正係数との関係を示す図である。
【図１４】ＥＧＲ弁通過ガス流量ｍｅｇｒを説明するための図である。
【図１５】排気圧Ｐｅ、排気温Ｔｅ、及びＰｅ／√Ｔｅを示す線図である。
【図１６】関数Φ（Ｐｍ／Ｐｅ）を示す線図である。
【図１７】筒内吸入空気流量と吸気管内圧力との関係の一例を示す線図である。
【符号の説明】
【０１３０】
１    機関本体
５    燃焼室
６    吸気弁
７    吸気ポート
８    排気弁
１１    燃料噴射弁
１３    吸気管
１８    スロットル弁
２２    ＥＧＲ弁
【図１】
2008069790.tif 000027
【図２】
2008069790.tif 000028
【図３】
2008069790.tif 000029
【図４】
2008069790.tif 000030
【図５】
2008069790.tif 000031
【図６】
2008069790.tif 000032
【図７】
2008069790.tif 000033
【図８】
2008069790.tif 000034
【図９】
2008069790.tif 000035
【図１０】
2008069790.tif 000036
【図１１】
2008069790.tif 000037
【図１２】
2008069790.tif 000038
【図１３】
2008069790.tif 000039
【図１４】
2008069790.tif 000040
【図１５】
2008069790.tif 000041
【図１６】
2008069790.tif 000042
【図１７】
2008069790.tif 000043
【手続補正書】
【提出日】平成19年12月27日(2007.12.27)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
燃焼室から排気通路に排出された排気ガスの少なくとも一部を吸気通路に再循環させる排気再循環通路と、該排気再循環通路を通る排気ガスの流量を調整する制御弁とを具備する内燃機関におけるスロットル弁通過空気流量を算出する装置であって、吸気管内圧力と、排気ガス再循環中の定常運転時における吸気管内圧力であって、スロットル弁の開度及び上記制御弁の開度に基づき算出される吸気管内圧力と、スロットル弁上流の圧力もしくは大気圧とに基づいて、過渡状態におけるスロットル弁通過空気流量を算出するスロットル弁通過空気流量算出装置。
【請求項２】
上記制御弁の開度をスロットル弁の開度に換算する手段を有していて、上記排気ガス再循環中の定常運転時における吸気管内圧力が、機関回転数と、上記制御弁の開度分を含めたスロットル弁の開度とに基づいて定められる、請求項１に記載のスロットル弁通過空気流量算出装置。
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】００１６
【補正方法】変更
【補正の内容】
【００１６】
上記課題を解決するために、１番目の発明によれば、燃焼室から排気通路に排出された排気ガスの少なくとも一部を吸気通路に再循環させる排気再循環通路と、該排気再循環通路を通る排気ガスの流量を調整する制御弁とを具備する内燃機関におけるスロットル弁通過空気流量を算出する装置であって、吸気管内圧力と、排気ガス再循環中の定常運転時における吸気管内圧力であって、スロットル弁の開度及び上記制御弁の開度に基づき算出される吸気管内圧力と、スロットル弁上流の圧力もしくは大気圧とに基づいて、過渡状態におけるスロットル弁通過空気流量を算出するスロットル弁通過空気流量算出装置が提供される。
【手続補正３】
【補正対象書類名】明細書
【補正対象項目名】００１７
【補正方法】変更
【補正の内容】
【００１７】
また２番目の発明によれば１番目の発明において、上記制御弁の開度をスロットル弁の開度に換算する手段を有していて、上記排気ガス再循環中の定常運転時における吸気管内圧力が、機関回転数と、上記制御弁の開度分を含めたスロットル弁の開度とに基づいて定められる。