JP4376563B2

JP4376563B2 - 内燃機関の制御装置

Info

Publication number: JP4376563B2
Application number: JP2003201213A
Authority: JP
Inventors: 晴文武藤
Original assignee: Toyota Motor Corp
Current assignee: Toyota Motor Corp
Priority date: 2003-07-24
Filing date: 2003-07-24
Publication date: 2009-12-02
Anticipated expiration: 2023-07-24
Also published as: JP2005042573A

Description

【０００１】
【発明の属する技術分野】
本発明は、内燃機関の制御装置に関する。
【０００２】
【従来の技術】
正確な空燃比制御を実現するためには吸入空気量を正確に把握することが必要である。吸入空気量を検出するために機関吸気系にエアフローメータを配置することが一般的であるが、エアフローメータは応答遅れを有するために、機関定常時において比較的正確な吸入空気量を検出することはできても、機関過渡時においては不正確な吸入空気量しか検出することができない。
【０００３】
機関過渡時においても比較的正確な吸入空気量を把握するために、機関吸気系をモデル化して吸入空気量を逐次計算することが提案されている（例えば、特許文献１参照）。
【０００４】
吸入空気量を算出するためには、スロットル弁の開口面積により変化するスロットル弁通過空気量を算出しなければならない。スロットル弁通過空気量は、絞りの式を使用して、スロットル弁の下流側圧力と上流側圧力との比を変数とする関数とスロットル弁の開口面積を変数とする関数との積の形の計算式によって表される。しかしながら、スロットル弁の開口面積を変数とする関数において、スロットル弁の開口形状等が考慮されていないために、この計算式を使用して実機に対して正確なスロットル弁通過空気量を算出することは困難である。
【０００５】
ところで、吸入空気量は、スロットル弁の下流側圧力に応じて変化するものであり、実機に適合させた下流側圧力の一次式により近似することができる。機関定常時において吸入空気量とスロットル弁通過空気量とは等しくなるために、現在のスロットル弁の開口面積（又は開度）と、これに対する定常時の下流側圧力（現在のスロットル弁開度に対して収束した場合の下流側圧力）とに基づき前述の計算式により算出される定常時のスロットル弁通過空気量は、この定常下流側圧力に基づき前述の一次式により算出される定常吸入空気量と等しいとすることができる。この関係を使用して、前述の計算式におけるスロットル弁の開口面積を変数とする関数は、定常下流側圧力に基づき算出された定常吸入空気量と、下流側圧力と上流側圧力との比を変数とする前述の関数において下流側圧力を定常下流側圧力とした関数との比によって置換することができ、ここで、定常下流側圧力を現在のスロットル弁開度に対する実機の適合値を使用すれば、定常吸入空気量を算出するのに使用される一次式も実機に適合しているために、スロットル弁通過空気量を正確に算出することができる。
【０００６】
【特許文献１】
特開２００１−４１０９５号公報
【特許文献２】
特開２００１−１３００３９号公報
【０００７】
【発明が解決しようとする課題】
前述の従来技術において、スロットル弁通過空気量を算出するのに使用される定常下流側圧力と定常吸入空気量を算出する一次式を特定するための係数及び定数とは、バルブタイミング等を含めた現在の機関状態毎に変化させており、そのためには、これらを実機に対して機関状態毎の適合値として設定しておかなければならず、この設定作業のために多大な工数が必要となっている。
【０００８】
従って、本発明の目的は、吸入空気量を推定するために、スロットル弁通過空気量が、スロットル弁の下流側圧力と上流側圧力との比を変数とする第一関数と、現在のスロットル弁の開口面積又は開度に対するスロットル弁の定常下流側圧力に基づき一次式により算出される定常吸入空気量と第一関数において下流側圧力を定常下流側圧力とした第二関数との比との積を使用して算出される内燃機関の制御装置において、スロットル弁通過空気量の算出に際して必要な実機に対する機関状態毎の適合値を減少することにより、適合値を設定するための作業工数を低減することである。
【０００９】
【課題を解決するための手段】
本発明による請求項１に記載の内燃機関の制御装置は、スロットル弁モデルにより現在のスロットル弁通過空気量を逐次的に算出すると共に、吸気弁モデル及び吸気管モデルにより気筒内へ供給される現在の吸入空気量を逐次的に算出する内燃機関の制御装置であって、
前記吸気弁モデルは、気筒内へ供給される吸入空気量がスロットル弁下流側圧力を変数とする一次式で表され、前記一次式の係数及び定数は、前記スロットル弁下流側圧力に基づき前記吸入空気量を算出する際に前記吸入空気量に影響する機関状態毎に設定され、
前記スロットル弁モデルは、現在のスロットル弁通過空気量を、スロットル弁下流側圧力に対するスロットル弁上流側圧力の比を変数とする関数と、スロットル弁の開口面積又は開度を変数とするスロットル弁の流量係数を含むもう一つの関数との積により表される式であって、
いずれの前記機関状態においても、任意のスロットル弁の開口面積又は開度に対する定常時の定常スロットル弁下流側圧力に基づき前記スロットル弁モデルにより算出される前記定常時の定常スロットル弁通過空気量が、前記定常時の定常スロットル弁下流側圧力に基づき前記吸気弁モデルにより算出される定常時の定常吸入空気量に等しいことに基づき、前記もう一つの関数を、前記関数において前記定常スロットル弁下流側圧力に対するスロットル弁上流側圧力の比を変数として代入した値と、前記吸気弁モデルから算出される前記定常吸入空気量とを使用して置換した式により算出する内燃機関の制御装置において、
現在のスロットル弁通過空気量は、前記スロットル弁モデルを使用して、前記関数において前記関数の前記比として現在のスロットル弁上流側圧力に対する現在のスロットル弁の開口面積又は開度における現在のスロットル弁下流側圧力の比を代入して算出される値と、
前記関数において前記関数の前記比として現在のスロットル弁上流側圧力に対する現在のスロットル弁の開口面積又は開度における定常時の定常スロットル弁下流側圧力の比を代入して算出される値、に対する前記一次式において前記スロットル弁下流側圧力として現在のスロットル弁の開口面積又は開度における定常時の定常スロットル弁下流側圧力を代入して算出される定常時の定常吸入空気量の比、との積により算出され、
前記定常スロットル弁下流側圧力は、予め定められた一つの機関状態に対してスロットル弁の開口面積毎又は開度毎に予め求められ、現在のスロットル弁の開口面積又は開度に対して予め求められた前記定常スロットル弁下流側圧力は、現在の機関状態が前記予め定められた一つの機関状態以外の機関状態の時においても、現在のスロットル弁通過空気量を算出するために必要な前記定常スロットル弁下流側圧力としてそのまま使用されると共に、現在のスロットル弁通過空気量を算出するために必要な前記定常吸入空気量を前記予め定められた一つの機関状態に対して係数及び定数が設定された前記一次式により算出するために前記スロットル弁下流側圧力として代入されることを特徴とする。
【００１０】
【発明の実施の形態】
図１は、本発明による制御装置が取り付けられる内燃機関を示す概略図である。同図において、１は機関本体であり、２は各気筒共通のサージタンクである。３はサージタンク２と各気筒とを連通する吸気枝管であり、４はサージタンク２の上流側の吸気通路である。各吸気枝管３には燃料噴射弁５が配置され、吸気通路４におけるサージタンク２の直上流側にはスロットル弁６が配置されている。ここで、スロットル弁６は、アクセルペダルに連動するものではなく、ステップモータ等の駆動装置によって自由に開度設定可能なものとしたが、アクセルペダルに連動するものでも良い。スロットル弁６には、スロットル弁開度を測定するためのスロットル弁開度センサが取り付けられている。機関本体１において、８は吸気弁であり、９は排気弁であり、１０はピストンである。
【００１１】
内燃機関１における燃焼空燃比を、例えば、理論空燃比等の所望空燃比にするためには、機関過渡時を含めて気筒内へ流入する吸入空気量を正確に推定することが必要とされる。本実施形態は、機関吸気系をモデル化して計算により吸入空気量を推定するようにしている。
【００１２】
先ず、スロットル弁６をモデル化することにより、吸気がスロットル弁６を通過する際のエネルギ保存則、運動量保存則、及び、状態方程式を使用して、今回のスロットル弁通過空気量ｍｔ_(i)（ｇ／ｓｅｃ）を、次式（１）によって表すことができる。以下の式を含めて、スロットル弁通過空気量等の変数の添え字（ｉ）は今回を示し、（ｉ−１）は前回を示している。
【数１】

【００１３】
ここで、μ_(i)は流量係数であり、Ａ_(i)はスロットル弁６の開口面積（ｍ³）である。もちろん、機関吸気系にアイドルスピードコントロールバルブ（ＩＳＣ弁）が設けられている時には、Ａ_(i)には、ＩＳＣ弁の開口面積が加えられる。流量係数及びスロットル弁の開口面積は、それぞれがスロットル弁開度ＴＡ_(i)（度）の関数となっており、図２及び３には、それぞれのスロットル弁開度ＴＡに対するマップが図示されている。Ｒは気体定数であり、Ｔａはスロットル弁上流側の吸気温度（Ｋ）であり、Ｐａはスロットル弁上流側の吸気通路圧力（ｋＰａ）（以下、上流側圧力と称する）であり、Ｐｍ_(i)はスロットル弁下流側の吸気管圧力（ｋＰａ）（以下、下流側圧力と称する）である。また、関数Φ（Ｐｍ_(i)／Ｐａ）は、比熱比κを使用して次式（２）によって表されるものであり、図４にはＰｍ／Ｐａに対するマップが図示されている。
【数２】

【００１４】
次いで、吸気弁をモデル化する。気筒内へ供給される吸入空気量ｍｃ_(i)（ｇ／ｓｅｃ）は、下流側圧力Ｐｍ_(i)に基づきほぼ線形に変化するものであるために、次式（３）に示す一次式によって表すことができる。
【数３】

【００１５】
ここで、Ｔｍ_(i)はスロットル弁下流側の吸気温度（Ｋ）であり、式（３）を特定するための係数及び定数ａ、ｂ、及びｃは、内燃機関毎に適合させて設定される適合値である。ここで、ａｂ−ｃは気筒内の残留既燃ガス量に相当する値であり、バルブオーバーラップ大きいほど吸気管へ既燃ガスが逆流するために増加する。このように、正確な吸入空気量ｍｃが算出するためには、これらの係数及び定数は、機関状態毎に設定されることが好ましい。すなわち、機関回転数毎だけでなく、内燃機関が可変バルブタンミング機構を有する場合にはバルブオーバーラップ量毎及び吸気弁の最大リフト量毎、また、その他の吸入空気量に影響する制御装置が設けられている場合にはその制御量毎に、ａ、ｂ及びｃの値をマップ化することが好ましい。
【００１６】
次いで、吸気管をモデル化する。吸気管内に存在する吸気の質量保存則、エネルギ保存則、及び、状態方程式を使用して、下流側圧力Ｐｍとスロットル弁下流側の吸気温度Ｔｍとの比における時間変化率は次式（４）によって表され、また、下流側圧力Ｐｍの時間変化率は次式（５）によって表される。ここで、Ｖは、機関吸気系のスロットル弁６下流側における吸気管の容積（ｍ³）であり、具体的には、吸気通路４のスロットル弁６下流側とサージタンク２と吸気枝管３との合計容積である。
【数４】

【００１７】
式（４）及び式（５）は離散化され、それぞれ、次式（６）及び（７）が得られ、式（７）によって今回の下流側圧力Ｐｍ_(i)が得られれば、式（６）によって今回の吸気管内の吸気温度Ｔｍ_(i)を得ることができ、また、式（１）によって現在のスロットル弁の開口面積Ａ_(i)及び現在の流量係数μ_(i)に基づき今回のスロットル弁通過空気量ｍｔ_(i)を得ることができる。さらに、式（３）によって、今回の吸気温度Ｔｍ_(i)及び今回の下流側圧力Ｐｍ_(i)に基づき今回の吸入空気量ｍｃ_(i)を得ることができる。次の計算時期においては、下流側圧力、吸気温度、スロットル弁通過空気量、吸入空気量をそれぞれ前回値として、式（７）により今回の下流側圧力Ｐｍ_(i)を算出する。これを繰り返して吸入空気量ｍｃを逐次計算する。式（６）及び（７）において、離散時間Δｔは、吸入空気量ｍｃを算出する時間間隔とされ、例えば８ｍｓである。また、スロットル弁通過空気量ｍｔ_(i)の算出に使用する今回のスロットル弁の開口面積Ａ_(i)及び流量係数μ_(i)は、現在のアクセルペダルの踏み込み量に対してスロットル弁の駆動装置（ステップモータ）の応答遅れ等が考慮されて推定される。
【数５】

【００１８】
このようにして吸入空気量ｍｃを算出することとなるが、実機に対して正確な吸入空気量を算出するためには、前述の各式が実機に対して各値を正確に算出しなければならない。そのためには、各式を式（３）のように実機に対して適合させることが好ましい。特に、スロットル弁通過空気量ｍｔを算出するための式（１）において、実機のスロットル弁は、開口面積が変化すると同時に開口形状も変化するものであり、この吸気流れの特性を、スロットル弁の開口面積Ａと流量係数μとの積によって正確に表すことは困難である。それにより、本実施形態では、式（１）を実機に対して適合させるように変形してスロットル弁通過空気量を算出するようにしている。
【００１９】
式（１）において、現在のスロットル弁の開口面積Ａ_(i)が暫く持続する定常時を考えれば、下流側吸気圧及びスロットル弁通過空気量は、それぞれ、収束して、定常下流側圧力Ｐｍｔａ及び定常スロットル弁通過空気量ｍｔｔａとなる。それにより、式（１）に定常下流側圧力Ｐｍｔａ及び定常スロットル弁通過空気量ｍｔｔａを代入して次式（８）を得ることができ、式（８）を次式（９）のように変形することができる。また、この定常スロットル弁通過空気量ｍｔｔａは、この定常状態における定常吸入空気量と等しくなるために、下流側圧力に基づく吸入空気量を算出する式（３）を使用すれば、式（９）は、定常下流側圧力Ｐｍｔａに基づき次式（１０）のように変形することができる。
【数６】

【００２０】
こうして、式（１）において計算誤差を発生し易い式（９）の左辺の関数は、式（１０）のように置換することができる。式（１０）において、定常下流側圧力Ｐｍｔａは、実機において現在のスロットル弁の開口面積（又は開度）に対して現在の機関状態に基づき適合させた適合値とし、また、式（１０）の分子である定常吸入空気量（ａ（Ｐｍｔａ−ｂ）＋ｃ）は、この定常下流側圧力Ｐｍｔａと、実機において現在の機関状態に基づき適合させた係数及び定数（ａ、ｂ、及びｃ）とを使用して算出すれば良く、それにより、実機に適合した正確なスロットル弁通過空気量ｍｔ_(i)を算出することができる。
【００２１】
しかしながら、このようにしてスロットル弁通過空気量ｍｔ_(i)を算出するためには、現在の機関状態に対して吸入空気量を算出するための一次式を特定するための係数及び定数（ａ、ｂ、及びｃ）だけでなく、スロットル弁の開口面積に対する定常下流側圧力Ｐｍｔａも、機関状態毎にマップ化しておかなければならない。前述したように、機関状態とは、機関回転数だけでなく、バルブオーバーラップ量、最大バルブリフト量、その他の吸入空気量に影響する制御量によっても変化するものであり、実機に適合させてスロットル弁の開口面積に対する定常下流側圧力Ｐｍｔａを機関状態毎に設定するのには莫大な工数が必要となってしまう。
【００２２】
本実施形態では、この定常下流側圧力Ｐｍｔａの設定工数を低減することを意図している。式（９）を参照すると、式（９）の左辺は、スロットル弁の開口面積Ａ_(i)及び流量係数μ_(i)だけを変数としており、すなわち、機関回転数やバルブオーバーラップ量等の機関状態とは無関係である。言い換えれば、式（９）の左辺は機関状態によっては変化しない。すなわち、式（９）の右辺において、機関状態が変化すれば、同じスロットル弁の開口面積に対する定常スロットル弁通過空気量ｍｔｔａは変化し、また、同じスロットル弁の開口面積に対する定常下流側圧力Ｐｍｔａも変化してこの定常下流側圧力Ｐｍｔａと上流側圧力Ｐａとの比を変数とする関数値Φも変化するが、定常スロットル弁通過空気量ｍｔｔａと関数値Φとの比は、機関状態が変化しても変化しないと言える。
【００２３】
こうして、式（１０）により式（１）の一部を置換してスロットル弁通過空気量ｍｔ_(i)を算出する際には、現在の機関状態とは関係なく、任意の機関状態における現在のスロットル弁の開口面積（又は開度）に対する定常下流側圧力Ｐｍｔａと、同じ機関状態における吸入空気量を算出するための一次式を使用すれば良い。すなわち、式（１）を次式（１１）のように変形してスロットル弁通過空気量ｍｔ_(i)を算出することができる。式（１１）において、現在のスロットル弁の開口面積（又は開度）に対する定常下流側圧力Ｐｍｔａは、特定の機関状態におけるものであり、この定常下流側圧力Ｐｍｔａに基づき定常吸入空気量を算出するための一次式を特定する係数及び定数（ａ１、ｂ１、及びｃ１）は、同じ特定の機関状態におけるものである。
【数７】

【００２４】
ところで、燃焼空燃比を正確に制御するためには、燃料噴射を開始する以前に気筒内への正確な吸入空気量を推定して、燃料噴射量を決定しなければならない。しかしながら、正確な吸入空気量を推定するためには、厳密には、吸気弁閉弁時における吸入空気流量を算出しなければならない。すなわち、燃料噴射量を決定する時において、現在の吸入空気量ｍｃ_(i)ではなく、吸気弁閉弁時における吸入空気量ｍｃ_(i+n)を算出しなければならない。これは、図１に示すような吸気枝管３に燃料を噴射する内燃機関だけでなく、吸気行程において筒内へ直接燃料を噴射する内燃機関においても同様である。
【００２５】
そのためには、現在において、現在のスロットル弁の開口面積Ａ_(i)（又は開度）だけでなく、吸気弁閉弁時までの時間Δｔ毎のスロットル弁の開口面積Ａ_(i ₊₁₎，Ａ_(i+2)，・・・Ａ_(i+n)に基づき、式（１１）において定常下流側圧力Ｐｍｔａを変化させ、各時間のスロットル弁通過空気量ｍｔを算出することが必要となる。
【００２６】
各時間のスロットル弁の開口面積Ａ（又は開度ＴＡ）は、現在の時間に対するアクセルペダルの踏み込み変化量に基づき、この踏み込み変化量が吸気弁閉弁時まで持続するとして、各時間のアクセルペダルの踏み込み量を推定し、それぞれの推定踏み込み量に対して、スロットル弁アクチュエータの応答遅れを考慮して決定することが考えられる。この方法は、スロットル弁がアクセルペダルと機械的に連結されている場合にも適用することができる。
【００２７】
しかしながら、こうして推定される吸気弁閉弁時におけるスロットル弁の開口面積Ａ_(i+n)は、あくまでも予測であり、実際と一致している保証はない。吸気弁閉弁時におけるスロットル弁の開口面積Ａ_(i+n)を実際と一致させるために、スロットル弁を遅れ制御するようにしても良い。アクセルペダルの踏み込み量が変化した時に、アクチュエータの応答遅れによって、スロットル弁開度は遅れて変化するが、この遅れ制御は、このスロットル弁の応答遅れを意図的に増大させるものである。
【００２８】
例えば、機関過渡時において、燃料噴射量を決定する時における現在のアクセルペダルの踏み込み量に対応するスロットル弁開度が、吸気弁閉弁時に実現されるように、実際の応答遅れ（無駄時間）を考慮してスロットル弁のアクチュエータを制御すれば、現在から吸気弁閉弁時までの時間毎のスロットル弁の開口面積Ａ_(i)，Ａ_(i+1)，・・・Ａ_(i+n)を正確に把握することができる。さらに具体的に言えば、アクセルペダルの踏み込み量が変化する時には、直ぐにアクチュエータへ作動信号を発するのではなく、燃料噴射量を決定する時から吸気弁閉弁時までの時間から無駄時間を差し引いた時間だけ経過した時にアクチュエータへの作動信号を発するようにするのである。もちろん、現在のアクセルペダルの踏み込み量に対応するスロットル弁開度を、吸気弁閉弁時以降に実現するようにスロットル弁の遅れ制御を実施しても良い。
【００２９】
【発明の効果】
本発明による請求項１に記載の内燃機関の制御装置によれば、スロットル弁モデルにより現在のスロットル弁通過空気量を逐次的に算出すると共に、吸気弁モデル及び吸気管モデルにより気筒内へ供給される現在の吸入空気量を逐次的に算出する内燃機関の制御装置であって、
前記吸気弁モデルは、気筒内へ供給される吸入空気量がスロットル弁下流側圧力を変数とする一次式で表され、前記一次式の係数及び定数は、前記スロットル弁下流側圧力に基づき前記吸入空気量を算出する際に前記吸入空気量に影響する機関状態毎に設定され、
前記スロットル弁モデルは、現在のスロットル弁通過空気量を、スロットル弁下流側圧力に対するスロットル弁上流側圧力の比を変数とする関数と、スロットル弁の開口面積又は開度を変数とするスロットル弁の流量係数を含むもう一つの関数との積により表される式であって、
いずれの前記機関状態においても、任意のスロットル弁の開口面積又は開度に対する定常時の定常スロットル弁下流側圧力に基づき前記スロットル弁モデルにより算出される前記定常時の定常スロットル弁通過空気量が、前記定常時の定常スロットル弁下流側圧力に基づき前記吸気弁モデルにより算出される定常時の定常吸入空気量に等しいことに基づき、前記もう一つの関数を、前記関数において前記定常スロットル弁下流側圧力に対するスロットル弁上流側圧力の比を変数として代入した値と、前記吸気弁モデルから算出される前記定常吸入空気量とを使用して置換した式により算出する内燃機関の制御装置において、
現在のスロットル弁通過空気量は、前記スロットル弁モデルを使用して、前記関数において前記関数の前記比として現在のスロットル弁上流側圧力に対する現在のスロットル弁の開口面積又は開度における現在のスロットル弁下流側圧力の比を代入して算出される値と、
前記関数において前記関数の前記比として現在のスロットル弁上流側圧力に対する現在のスロットル弁の開口面積又は開度における定常時の定常スロットル弁下流側圧力の比を代入して算出される値、に対する前記一次式において前記スロットル弁下流側圧力として現在のスロットル弁の開口面積又は開度における定常時の定常スロットル弁下流側圧力を代入して算出される定常時の定常吸入空気量の比、との積により算出され、
前記定常スロットル弁下流側圧力は、予め定められた一つの機関状態に対してスロットル弁の開口面積毎又は開度毎に予め求められ、現在のスロットル弁の開口面積又は開度に対して予め求められた前記定常スロットル弁下流側圧力は、現在の機関状態が前記予め定められた一つの機関状態以外の機関状態の時においても、現在のスロットル弁通過空気量を算出するために必要な前記定常スロットル弁下流側圧力としてそのまま使用されると共に、現在のスロットル弁通過空気量を算出するために必要な前記定常吸入空気量を前記予め定められた一つの機関状態に対して係数及び定数が設定された前記一次式により算出するために前記スロットル弁下流側圧力として代入されるようになっている。
それにより、定常下流側圧力を実機に対して機関状態毎に適合させる必要はなく、この適合値を設定するための設定作業の工数を低減することができる。ここで、予め定められた機関状態は、任意の機関状態として良く、例えば、機関状態を表す一つの因子である機関回転数に関して、現在の機関状態の機関回転数が、予め定められた機関状態の機関回転数より高くなっても低くなっても良い。
【図面の簡単な説明】
【図１】本発明による吸入空気量推定装置が取り付けられる内燃機関の概略図である。
【図２】スロットル弁開度ＴＡと流量係数μとの関係を示すマップである。
【図３】スロットル弁開度ＴＡとスロットル弁の開口面積Ａとの関係を示すマップである。
【図４】吸気管圧力Ｐｍと大気圧Ｐａとの比と、関数Φとの関係を示すマップである。
【符号の説明】
１…機関本体
２…サージタンク
３…吸気枝管
４…吸気通路
６…スロットル弁

Claims

スロットル弁モデルにより現在のスロットル弁通過空気量を逐次的に算出すると共に、吸気弁モデル及び吸気管モデルにより気筒内へ供給される現在の吸入空気量を逐次的に算出する内燃機関の制御装置であって、
前記吸気弁モデルは、気筒内へ供給される吸入空気量がスロットル弁下流側圧力を変数とする一次式で表され、前記一次式の係数及び定数は、前記スロットル弁下流側圧力に基づき前記吸入空気量を算出する際に前記吸入空気量に影響する機関状態毎に設定され、
前記スロットル弁モデルは、現在のスロットル弁通過空気量を、スロットル弁下流側圧力に対するスロットル弁上流側圧力の比を変数とする関数と、スロットル弁の開口面積又は開度を変数とするスロットル弁の流量係数を含むもう一つの関数との積により表される式であって、
いずれの前記機関状態においても、任意のスロットル弁の開口面積又は開度に対する定常時の定常スロットル弁下流側圧力に基づき前記スロットル弁モデルにより算出される前記定常時の定常スロットル弁通過空気量が、前記定常時の定常スロットル弁下流側圧力に基づき前記吸気弁モデルにより算出される定常時の定常吸入空気量に等しいことに基づき、前記もう一つの関数を、前記関数において前記定常スロットル弁下流側圧力に対するスロットル弁上流側圧力の比を変数として代入した値と、前記吸気弁モデルから算出される前記定常吸入空気量とを使用して置換した式により算出する内燃機関の制御装置において、
現在のスロットル弁通過空気量は、前記スロットル弁モデルを使用して、前記関数において前記関数の前記比として現在のスロットル弁上流側圧力に対する現在のスロットル弁の開口面積又は開度における現在のスロットル弁下流側圧力の比を代入して算出される値と、
前記関数において前記関数の前記比として現在のスロットル弁上流側圧力に対する現在のスロットル弁の開口面積又は開度における定常時の定常スロットル弁下流側圧力の比を代入して算出される値、に対する前記一次式において前記スロットル弁下流側圧力として現在のスロットル弁の開口面積又は開度における定常時の定常スロットル弁下流側圧力を代入して算出される定常時の定常吸入空気量の比、との積により算出され、
前記定常スロットル弁下流側圧力は、予め定められた一つの機関状態に対してスロットル弁の開口面積毎又は開度毎に予め求められ、現在のスロットル弁の開口面積又は開度に対して予め求められた前記定常スロットル弁下流側圧力は、現在の機関状態が前記予め定められた一つの機関状態以外の機関状態の時においても、現在のスロットル弁通過空気量を算出するために必要な前記定常スロットル弁下流側圧力としてそのまま使用されると共に、現在のスロットル弁通過空気量を算出するために必要な前記定常吸入空気量を前記予め定められた一つの機関状態に対して係数及び定数が設定された前記一次式により算出するために前記スロットル弁下流側圧力として代入されることを特徴とする内燃機関の制御装置。