JPH08232751A - Suction air amount estimating device for internal combustion engine - Google Patents

Suction air amount estimating device for internal combustion engine

Info

Publication number
JPH08232751A
JPH08232751A JP3885395A JP3885395A JPH08232751A JP H08232751 A JPH08232751 A JP H08232751A JP 3885395 A JP3885395 A JP 3885395A JP 3885395 A JP3885395 A JP 3885395A JP H08232751 A JPH08232751 A JP H08232751A
Authority
JP
Japan
Prior art keywords
air amount
throttle
intake air
chamber
internal combustion
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP3885395A
Other languages
Japanese (ja)
Other versions
JP3337338B2 (en
Inventor
Masuhiro Yoshizaki
益博 吉崎
Isao Komoriya
勲 小森谷
Shinichi Kitajima
真一 北島
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Honda Motor Co Ltd
Original Assignee
Honda Motor Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Honda Motor Co Ltd filed Critical Honda Motor Co Ltd
Priority to JP03885395A priority Critical patent/JP3337338B2/en
Priority to US08/607,866 priority patent/US5597951A/en
Publication of JPH08232751A publication Critical patent/JPH08232751A/en
Application granted granted Critical
Publication of JP3337338B2 publication Critical patent/JP3337338B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Combined Controls Of Internal Combustion Engines (AREA)

Abstract

PURPOSE: To perform accurate estimation of a cylinder suction air amount by a method wherein the fluctuation of a filling air amount are detected by a detecting means and the primary delay correction factor of a throttle opening is individually set according to a case wherein a filling air amount is increased and a case wherein a filling air amount is decreased. CONSTITUTION: When a suction air amount is estimated by a control unit 50, an amount Gc of suction air to a combustion chamber when an engine is in a steady operation state is determined based on the number of revolutions of an engine and a pressure in a chamber. Further, based on a throttle opening and the pressure in the chamber, the first effective opening area A of a throttle is determined. The primary delay value of an effective opening area A is then determined and the value forms the second effective opening area B of the throttle and an actual suction amount Gc is determined by a formula of Gc= Gc'×(A/B). In this case, the fluctuation of a filling air amount is detected and according to a detecting result, the primary delay correction factor of the throttle opening is set. Based on the primary delay correction factor, the primary delay value of the first effective opening area A is determined.

Description

【発明の詳細な説明】Detailed Description of the Invention

【0001】[0001]

【産業上の利用分野】本発明は、内燃機関に吸入される
空気量を推定する内燃機関の吸入空気量推定装置に関す
る。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to an intake air amount estimating apparatus for an internal combustion engine, which estimates the amount of air taken into the internal combustion engine.

【0002】[0002]

【従来の技術】近年、内燃機関の吸気系に流体力学モデ
ルを適用し、モデル式によって正しい気筒吸入空気量を
推定する手法が提案されている。本出願人も、先に、特
開平6−74076号において、スロットル弁をオリフ
ィスと見なし、スロットル弁の前後の差圧から絞り式流
量形の原理式を用いてスロットル通過空気量を求め、こ
の値をもとに気筒吸入空気量を算出する手法を提案し
た。
2. Description of the Related Art In recent years, a method has been proposed in which a fluid dynamics model is applied to the intake system of an internal combustion engine, and a correct cylinder intake air amount is estimated by a model formula. The applicant of the present invention has previously considered that the throttle valve is an orifice in Japanese Patent Laid-Open No. 6-74076, and calculates the throttle passing air amount from the differential pressure before and after the throttle valve using a throttle type flow rate principle formula. We proposed a method to calculate the cylinder intake air amount based on.

【0003】[0003]

【発明が解決しようとする課題】しかし、このような流
体力学モデルは、あくまでも理想状態を前提としてお
り、モデル化の際に使用する各種定数(比熱比など)を
正確に特定することも困難であり、さらに、流体力学の
式はべき乗、平方根などの計算を必要であるため、実用
上は近似値を用いている。これらの理由により、モデル
化誤差を拭いきることはできないのが実情である。
However, such a fluid dynamics model is premised on an ideal state, and it is difficult to accurately specify various constants (specific heat ratio etc.) used in modeling. In addition, since the equation of fluid dynamics requires calculation of exponentiation, square root, etc., an approximate value is practically used. For these reasons, the modeling error cannot be wiped out.

【0004】そこで、このモデルから求められる気筒吸
入空気量と実際の気筒吸入空気量との誤差を如何に少な
くするかが極めて重要となる。この点を改善すべく、本
出願人は、流体力学モデルを前提としながらも複雑な演
算を必要とせずにモデル式の誤差を吸収し、かつ機関運
転の過渡状態や劣化、バラツキ、経年変化などを解消
し、気筒吸入空気量をより正確に推定する手法を提案し
ている(特願平5−208835号)。
Therefore, it is extremely important how to reduce the error between the cylinder intake air amount obtained from this model and the actual cylinder intake air amount. In order to improve this point, the present applicant absorbs the error of the model formula without requiring complicated calculation while assuming the fluid dynamics model, and transient state and deterioration of engine operation, variation, secular change, etc. And a method of more accurately estimating the cylinder intake air amount is proposed (Japanese Patent Application No. 5-208835).

【0005】本発明の目的は、すでに提案されている手
法にさらに改良を加え、流体力学モデルから得られる気
筒吸入空気量と実際の気筒吸入空気量との誤差をより低
減し、気筒吸入空気量を一層正確に推定し得る内燃機関
の吸入空気量推定装置を提供することにある。
The object of the present invention is to further improve the method already proposed, further reduce the error between the cylinder intake air amount obtained from the fluid dynamics model and the actual cylinder intake air amount, and improve the cylinder intake air amount. It is an object of the present invention to provide an intake air amount estimation device for an internal combustion engine, which can estimate more accurately.

【0006】[0006]

【課題を解決するための手段】そこで、請求項1のかか
る内燃機関の吸入空気量推定装置は、機関回転数Ne、
スロットル開度θTH及びチャンバ内圧力Pb を含む、内
燃機関の運転状態を検出する第1の手段と、少なくとも
第1の手段で検出された機関回転数Ne とチャンバ内圧
力Pb とを基に、定常運転状態時における燃焼室への吸
入空気量Gc ’を求める第2の手段と、少なくとも第1
の手段で検出されたスロットル開度θTHとチャンバ内圧
力Pb とを基に、スロットルの第1有効開口面積Aを求
める第3の手段と、このスロットルの第1有効開口面積
Aの1次遅れ値を求め、その値をスロットルの第2有効
開口面積ADELAY とする第4の手段と、第3の手段で求
めた第1有効開口面積Aと第4の手段で求めた第2有効
開口面積ADELAY とを基に、吸入空気量Gc ’を補正す
ることで、実吸入空気量Gc を、Gc =Gc ’×(A/
ADELAY )として求める第5の手段とを備えており、ス
ロットルから燃焼室に至る吸気系をモデル化して1つの
チャンバとみなし、この吸気系モデルに基づき、この燃
焼室へ吸入される実吸入空気量Gc を推定する内燃機関
の吸入空気量推定装置において、この第4の手段は、チ
ャンバ内に充填される充填空気量Gbの増加、減少を検
知する検知手段と、検知手段の検知結果を受けて充填空
気量Gb が増加した場合と減少した場合とで、スロッ
トル開度の1次遅れ補正係数Bを個々に設定する設定手
段と、この設定手段で設定された1次遅れ補正係数Bを
基に、スロットルの第1有効開口面積Aの1次遅れ値を
求める演算手段とを備えて構成する。
Therefore, the intake air amount estimating apparatus for an internal combustion engine according to the first aspect of the present invention is configured so that the engine speed Ne,
Based on the first means for detecting the operating state of the internal combustion engine, which includes the throttle opening θTH and the chamber pressure Pb, and at least the engine speed Ne and the chamber pressure Pb detected by the first means, the steady state is established. Second means for determining the intake air amount Gc 'into the combustion chamber during operation, and at least the first means
Means for obtaining the first effective opening area A of the throttle on the basis of the throttle opening θTH and the chamber pressure Pb detected by the means, and the first-order lag value of the first effective opening area A of the throttle. And the fourth effective means for determining the value as the second effective opening area ADELAY of the throttle, the first effective opening area A determined by the third means and the second effective opening area ADELAY determined by the fourth means. By correcting the intake air amount Gc 'on the basis of the following equation, the actual intake air amount Gc is calculated as Gc = Gc' × (A /
ADELAY) and a fifth means for calculating, and the intake system from the throttle to the combustion chamber is modeled and regarded as one chamber, and based on this intake system model, the actual intake air amount sucked into this combustion chamber In the intake air amount estimating apparatus for an internal combustion engine for estimating Gc, the fourth means receives a detection means for detecting an increase or a decrease in the filling air amount Gb filled in the chamber and a detection result of the detection means. Based on the setting means for individually setting the first-order lag correction coefficient B of the throttle opening depending on whether the filled air amount Gb increases or decreases, and the first-order lag correction coefficient B set by this setting means. And a calculating means for obtaining a first-order lag value of the first effective opening area A of the throttle.

【0007】また、請求項2にかかる内燃機関の吸入空
気量推定装置では、前述した検知手段を、スロットル開
度θTHの増加及び減少に基づいて、充填空気量Gb の増
加、減少を検知する手段として構成する。
Further, in the intake air amount estimating apparatus for an internal combustion engine according to a second aspect of the present invention, the above-mentioned detecting means is means for detecting an increase or decrease in the filled air amount Gb based on an increase or decrease in the throttle opening θTH. Configure as.

【0008】さらに、請求項3にかかる内燃機関の吸入
空気量推定装置では、請求項1及び2における検知手段
を、インテークマニホルドの容積変化、内燃機関のバル
ブタイミングの変化、大気圧の変化及び目標空燃比の変
化のうち、少なくともいずれか1つの変化に基づいて、
充填空気量Gb の増加、減少を検知する手段として構成
する。
Further, in the intake air amount estimating apparatus for an internal combustion engine according to a third aspect, the detecting means according to the first and second aspects is provided with the intake manifold volume change, the valve timing change of the internal combustion engine, the atmospheric pressure change and the target. Based on at least one of the changes in the air-fuel ratio,
It is configured as a means for detecting the increase or decrease of the filling air amount Gb.

【0009】なお、本発明における定常運転状態時と
は、充填空気量Gb の増減変化が実質的に無くほぼ一定
である期間をいい、過渡運転状態時とは、充填空気量G
b が変動している期間をいう。
In the present invention, the steady operation state means a period during which the increase / decrease of the filling air amount Gb is substantially constant and substantially constant, and the transient operation state means the filling air amount Gb.
The period during which b is fluctuating.

【0010】[0010]

【作用】請求項1における内燃機関の吸入空気量推定装
置では、検知手段によって充填空気量Gb の増加、減少
を検知し、かつ、この検知結果をもとに、設定手段によ
って、充填空気量Gb が増加した場合と減少した場合と
に応じて、スロットル開度の1次遅れ補正係数Bを個々
に設定する。このような設定を行うことにより、充填空
気量Gb の増加時と減少時とでは吸入空気の挙動が相異
するが、充填空気量Gb の増加時と減少時において、そ
れぞれ好適に追従するように、スロットル開度の1次遅
れ補正係数Bが設定される。
In the intake air amount estimating apparatus for the internal combustion engine according to the first aspect, the detecting means detects the increase or decrease of the filling air amount Gb, and based on the detection result, the setting means detects the filling air amount Gb. The first-order lag correction coefficient B for the throttle opening is individually set depending on whether the value increases or decreases. By performing such a setting, the behavior of the intake air differs when the filling air amount Gb increases and decreases, but the intake air behaves appropriately when the filling air amount Gb increases and decreases. The first-order delay correction coefficient B of the throttle opening is set.

【0011】また、請求項2にかかる内燃機関の吸入空
気量推定装置では、充填空気量Gbが変動する場合の主
要因となるスロットル弁の動きに着目した。すなわち、
スロットル弁が開き方向に変移した場合には運転状態が
加速状態にあり、この場合には充填空気量Gb が増加傾
向にある。また、スロットル弁が閉じ方向に変移した場
合には運転状態が減速状態にあり、この場合には充填空
気量Gb が減少傾向にある。このように、スロットル弁
の動きを検知手段で検知することで、充填空気量Gb の
増減を検知できる。
Further, in the intake air amount estimating apparatus for the internal combustion engine according to the second aspect of the present invention, attention is paid to the movement of the throttle valve which is a main factor when the filling air amount Gb changes. That is,
When the throttle valve shifts in the opening direction, the operating state is in the accelerating state, and in this case, the filling air amount Gb tends to increase. Further, when the throttle valve shifts in the closing direction, the operating state is in a decelerating state, and in this case, the filling air amount Gb tends to decrease. In this way, by detecting the movement of the throttle valve by the detecting means, it is possible to detect the increase or decrease of the filling air amount Gb.

【0012】また、請求項3にかかる内燃機関の吸入空
気量推定装置では、スロットル弁の動き以外に、インテ
ークマニホルドの容積変化等によっても充填空気量Gb
が増減するため、これらを検出手段で検知することで充
填空気量Gb の増減を検知できる。
Further, in the intake air amount estimating apparatus for the internal combustion engine according to the third aspect, the filling air amount Gb is caused not only by the movement of the throttle valve but also by the volume change of the intake manifold.
Is increased or decreased, the increase or decrease of the filling air amount Gb can be detected by detecting them by the detection means.

【0013】[0013]

【実施例】以下、本発明の実施例を添付図面を参照して
説明する。
Embodiments of the present invention will be described below with reference to the accompanying drawings.

【0014】図1において、符号10は4気筒の内燃機
関を示しており、吸気路12の先端に配置されたエアク
リーナ14から導入された吸気は、スロットル弁16で
その流量を調節されつつサージタンク18とインテーク
マニホルド20を経て第1〜第4気筒に流入される。各
気筒の吸気弁(図示せず)の付近には、インジェクタ2
2が設けれており、ここから燃料を噴射する。噴射され
た燃料と吸気とが一体となった混合気は、各気筒内で点
火プラグで点火され燃焼され、内燃機関のピストンを駆
動する。燃焼後の排気ガスは排気弁を介してエキゾース
トマニホルド24に放出され、エキゾーストパイプ26
を経て三元触媒コンバータ28で浄化され機関外に排出
される。
In FIG. 1, reference numeral 10 indicates a four-cylinder internal combustion engine, in which the intake air introduced from an air cleaner 14 arranged at the tip of an intake passage 12 is adjusted in its flow rate by a throttle valve 16 and a surge tank. It is introduced into the first to fourth cylinders via the intake manifold 18 and the intake manifold 20. An injector 2 is provided near the intake valve (not shown) of each cylinder.
2 is provided and fuel is injected from here. The air-fuel mixture in which the injected fuel and intake air are integrated is ignited by a spark plug in each cylinder and burned to drive a piston of the internal combustion engine. The exhaust gas after combustion is discharged to the exhaust manifold 24 through the exhaust valve, and the exhaust pipe 26
After that, it is purified by the three-way catalytic converter 28 and discharged to the outside of the engine.

【0015】また、内燃機関のディストリビュータ(図
示せず)内には、ピストンのクランク角度位置を検出す
るクランク角センサ34が設けられており、この他、ス
ロットル弁16の開度θTHを検出するスロットル開度セ
ンサ36、スロットル弁16下流の吸気圧力(チャンバ
内圧力)Pb を絶対圧力で検出する吸気圧センサ38も
設けられている。また、スロットル弁16の上流側に
は、大気圧Paを検出する大気圧センサ40、吸入空気
の温度TAを検出する吸気温センサ42、吸入空気の湿
度を検出する湿度センサ44が設けられている。さら
に、吸気路12には、2次空気量の調整用として、スロ
ットル弁16の前後の吸気路をバイパスするバイパス路
32が設けられており、電磁弁90を駆動することによ
り、このバイパス路32の開閉を行う。また、排気系に
おいてエキゾーストマニホルド24の下流側で三元触媒
コンバータ28の上流側には、酸素濃度検出素子からな
る広域空燃比46が設けられ、排気ガスの空燃比を検出
する。これらセンサ34などの出力は、制御ユニット5
0に送られる。
A crank angle sensor 34 for detecting the crank angle position of the piston is provided in the distributor (not shown) of the internal combustion engine. In addition to this, a throttle for detecting the opening degree θTH of the throttle valve 16 is provided. An opening sensor 36 and an intake pressure sensor 38 for detecting the intake pressure (chamber pressure) Pb downstream of the throttle valve 16 as an absolute pressure are also provided. An atmospheric pressure sensor 40 that detects the atmospheric pressure Pa, an intake air temperature sensor 42 that detects the temperature TA of the intake air, and a humidity sensor 44 that detects the humidity of the intake air are provided upstream of the throttle valve 16. . Further, the intake passage 12 is provided with a bypass passage 32 for bypassing the intake passages before and after the throttle valve 16 for adjusting the secondary air amount. By driving the solenoid valve 90, the bypass passage 32 is provided. Open and close. In the exhaust system, a wide-range air-fuel ratio 46 including an oxygen concentration detecting element is provided downstream of the exhaust manifold 24 and upstream of the three-way catalytic converter 28 to detect the air-fuel ratio of the exhaust gas. The outputs of these sensors 34 and the like are output to the control unit 5
Sent to 0.

【0016】図2に制御ユニットの詳細を示す。広域空
燃比センサ46の出力は検出回路52に入力され空燃比
A/Fが検出される。この検出回路52の出力は、A/
D変換回路54を介し、CPU56、ROM58及びR
AM60からなるマイクロ・コンピュータ内に取り込ま
れ、RAM60に格納される。同様に、スロットル開度
センサ36などのアナログ出力は、レベル変換回路6
2、マルチプレクサ64及び第2のA/D変換回路66
を介して、マイクロ・コンピュータ内に入力される。ま
た、クランク角センサ34の出力は、波形整形回路68
で波形整形された後、カウンタ70で出力値がカウント
され、カウント値はマイクロ・コンピュータ内に入力さ
れる。マイクロ・コンピュータにおいて、CPU56は
ROM58に格納された命令に従って、後述するような
制御値を演算し、駆動回路72を介して各気筒のインジ
ェクタ22を駆動する。
FIG. 2 shows the details of the control unit. The output of the wide range air-fuel ratio sensor 46 is input to the detection circuit 52 and the air-fuel ratio A / F is detected. The output of the detection circuit 52 is A /
CPU 56, ROM 58 and R via D conversion circuit 54
It is taken into the microcomputer composed of AM60 and stored in RAM60. Similarly, the analog output of the throttle opening sensor 36 and the like is supplied to the level conversion circuit 6
2, multiplexer 64 and second A / D conversion circuit 66
Is input to the microcomputer via the. In addition, the output of the crank angle sensor 34 is the waveform shaping circuit 68.
After the waveform is shaped by, the output value is counted by the counter 70, and the count value is input into the microcomputer. In the microcomputer, the CPU 56 calculates a control value, which will be described later, according to the instruction stored in the ROM 58, and drives the injector 22 of each cylinder via the drive circuit 72.

【0017】ここで、本実施例において採用する、流体
力学モデルを用いた気筒吸入空気量の推定手法の基本原
理を概略的に説明する(参照:特願平6−197,23
8号)。
Here, the basic principle of the method for estimating the cylinder intake air amount using the fluid dynamics model, which is adopted in this embodiment, will be schematically described (see Japanese Patent Application No. 6-197, 23).
No. 8).

【0018】この手法は、内燃機関の吸気系における空
気の挙動を物理式によりモデル化し(図3)、スロット
ルを通過する空気量Gthとチャンバへ充填される空気量
Gbとに基づいて、気筒内に吸入される吸入空気量(実
吸入空気量)Gc を推定する手法である。なお、ここで
「チャンバ」とは、いわゆるサージタンク相当部位のみ
ならず、スロットルから吸気ポートに至る間において空
気が流通する全ての部位を含み、実際にチャンバとして
働く実効容積を意味するものとする。
In this method, the behavior of air in the intake system of an internal combustion engine is modeled by a physical formula (FIG. 3), and the amount of air passing through the throttle Gth and the amount of air filled in the chamber Gb are used to calculate the cylinder interior. This is a method of estimating the intake air amount (actual intake air amount) Gc that is taken into the air. It should be noted that the term "chamber" as used herein means not only a so-called surge tank-corresponding portion but also all portions through which air flows from the throttle to the intake port, and means an effective volume that actually acts as a chamber. .

【0019】まず、図4に示すように、スロットル開度
θTHから予め設定した特性に従ってスロットルの投影面
積(吸気管の長手方向へのスロットルの投影面積)Sを
求める。他方、図5に示すように、スロットル開度θTH
と吸気圧力(チャンバ内圧力)Pbについて設定した別
の特性に従って係数C(流量係数αと気体の膨張補正係
数εの積)を求め、両者を乗じてスロットルの有効開口
面積Aを求める。なお、いわゆるスロットル全開領域で
は、スロットルが絞りとして機能しなくなるため、機関
回転数ごとにそれぞれスロットル全開領域を臨界値とし
て求めておき、検出されたスロットル開度がこの臨界値
を超えたときは、この臨界値をスロットル開度とする。
First, as shown in FIG. 4, a throttle projection area (throttle projection area in the longitudinal direction of the intake pipe) S is obtained from the throttle opening θTH according to a preset characteristic. On the other hand, as shown in FIG. 5, the throttle opening θTH
And a coefficient C (product of flow coefficient α and expansion correction coefficient ε of gas) according to another characteristic set for the intake pressure (chamber pressure) Pb, and multiplied by both to obtain the effective opening area A of the throttle. Since the throttle does not function as a throttle in the so-called throttle fully open region, the throttle fully open region is obtained as a critical value for each engine speed, and when the detected throttle opening exceeds this critical value, This critical value is the throttle opening.

【0020】次いで、気体の状態方程式に基づく数1式
から、チャンバ内の空気量Gb を求め、チャンバ内の圧
力変化ΔPから数2式に従って、今回チャンバ内に充填
された空気量ΔGb を求める。今回チャンバ内に充填さ
れた空気量は、当然ながら気筒燃焼室に吸入されないも
のとすれば、単位時間ΔT当たりの吸入空気量Gc は、
数3式のように表すことができる。
Next, the air amount Gb in the chamber is obtained from the equation 1 based on the equation of state of the gas, and the air amount ΔGb filled in the chamber this time is obtained from the pressure change ΔP in the chamber according to the equation 2. If the amount of air filled in the chamber this time is naturally not taken into the cylinder combustion chamber, the amount of intake air Gc per unit time ΔT is
It can be expressed as in Equation 3.

【0021】[0021]

【数1】 [Equation 1]

【0022】[0022]

【数2】 [Equation 2]

【0023】[0023]

【数3】 (Equation 3)

【0024】他方、前述したROM58には、図6に示
すように、定常運転状態時の燃料噴射量Timapを、いわ
ゆるスピードデンシティ方式に基づいて機関回転数Ne
と吸気圧力Pb とから検索できるように予め設定してマ
ップ化して格納している。また、ここで検索された燃料
噴射量Timapは機関回転数Ne と吸気圧力Pb に応じて
決定される目標空燃比A/Fに応じて設定されることか
ら、図7に示すように、目標空燃比A/Fの基本値KBS
も機関回転数Ne と吸気圧力Pb とから自在に検索でき
るように、予めマップ化して格納しておく。なお、燃料
噴射量Timapは、前述した流体力学モデルで、定常運転
状態時を満足するように設定する。直接的にはインジェ
クタ22の開弁時間を単位として設定する。
On the other hand, as shown in FIG. 6, the above-mentioned ROM 58 shows the fuel injection amount Timap in the steady operation state as the engine speed Ne based on the so-called speed density method.
It is set in advance so that it can be retrieved from the intake pressure Pb and the intake pressure Pb, and is stored as a map. Further, since the fuel injection amount Timap retrieved here is set according to the target air-fuel ratio A / F which is determined according to the engine speed Ne and the intake pressure Pb, as shown in FIG. Basic value of fuel ratio A / F KBS
Is stored in advance as a map so that it can be freely searched from the engine speed Ne and the intake pressure Pb. The fuel injection amount Timap is set so as to satisfy the steady operation state in the fluid dynamics model described above. The valve opening time of the injector 22 is directly set as a unit.

【0025】ここで、定常運転状態時のある条件下(機
関回転数Ne1と吸気圧力Pb1によって規定される条件
下)において、マップ検索によって決定した燃料噴射量
Timap1は数4式に示す通りとなる。
Here, under a certain condition in a steady operation state (condition defined by the engine speed Ne1 and the intake pressure Pb1), the fuel injection amount Timap1 determined by the map search becomes as shown in the equation (4). .

【0026】[0026]

【数4】 [Equation 4]

【0027】前述したモデル式を満足するようにこのマ
ップ値を作成しておくことで、流体力学モデルに基づい
て決定される、定常運転状態時の燃料噴射量Timap1’
は、マップ検索によって決定した燃料噴射量Timap1と
当然に一致する。
By preparing this map value so as to satisfy the above-mentioned model formula, the fuel injection amount Timap1 'in the steady operation state determined based on the fluid dynamics model.
Naturally agrees with the fuel injection amount Timap1 determined by the map search.

【0028】一方、流体力学モデルに基づいて決定され
る、定常運転状態時の燃料噴射量Timap1’及び過渡運
転状態時の燃料噴射量Timap2’は、目標空燃比を理論
空燃比(14.7:1)とすると、数5式及び数6式で
表される。
On the other hand, the fuel injection amount Timap1 'in the steady operation state and the fuel injection amount Timap2' in the transient operation state, which are determined based on the fluid dynamics model, are obtained by changing the target air-fuel ratio to the theoretical air-fuel ratio (14.7: 1), it can be expressed by Equation 5 and Equation 6.

【0029】[0029]

【数5】 (Equation 5)

【0030】[0030]

【数6】 (Equation 6)

【0031】この両式において、定常運転状態時のスロ
ットル通過空気量Gth1 と過渡運転状態時のスロットル
通過空気量Gth2 とを比較すると、スロットルの有効開
口面積A1とAのみが異なっている。従って、過渡運転
状態時のスロットル通過空気量Gth2 は、数7式で表す
ことができる。
In both of these equations, comparing the throttle passing air amount Gth1 in the steady operation state and the throttle passing air amount Gth2 in the transient operation state, only the effective opening areas A1 and A of the throttle differ. Therefore, the throttle passing air amount Gth2 in the transient operation state can be expressed by the equation (7).

【0032】[0032]

【数7】 (Equation 7)

【0033】数3式及び数7式に基づき、定常時のスロ
ットル弁の有効開口面積と過渡時のスロットル弁の有効
開口面積との比を用いることによって、定常時のスロッ
トル通過空気量Gth1 を基に、過渡運転状態時のスロッ
トル通過空気量Gth2 を表現することができる。
Based on the equations (3) and (7), the ratio of the effective opening area of the throttle valve in the steady state to the effective opening area of the throttle valve in the transient state is used to determine the throttle passing air amount Gth1 in the steady state. In addition, the throttle passing air amount Gth2 in the transient operation state can be expressed.

【0034】さらに、現在のスロットル弁の有効開口面
積をAとし、定常運転状態時のスロットル弁の有効開口
面積をA1とすると、定常運転状態時のスロットル弁の
有効開口面積A1は、現在のスロットル弁の有効開口面
積Aの一次遅れとして把握できる(図8)。即ち、Aの
一次遅れを「ADELAY 」と呼ぶと、A1とADELAY と
は、ほぼ同様の値となっていることが分る。従って、流
体力学モデルの考え方を用いてモデルを近似する際に
は、A/(Aの一次遅れ)を用いればよい。
Further, assuming that the current effective opening area of the throttle valve is A and the effective opening area of the throttle valve in the steady operation state is A1, the effective opening area A1 of the throttle valve in the steady operation state is the current throttle value. It can be grasped as the primary delay of the effective opening area A of the valve (Fig. 8). That is, when the first-order lag of A is called "ADELAY", it can be seen that A1 and ADELAY have almost the same value. Therefore, when approximating the model using the concept of the fluid dynamics model, A / (first-order lag of A) may be used.

【0035】図9に示すように、過渡運転状態では、ス
ロットルが開かれた瞬間は、スロットル前後の差圧が大
きいため、スロットル通過空気量が一気に流れ、次第に
定常状態に落ちつくが、その過渡運転状態におけるスロ
ットル通過空気量Gthは、この比A/ADELAY で表現で
き、定常運転状態では、図9の下側のグラフで示される
ように、この値が”1”となる。以下、この比を「RATI
O-A 」と呼ぶ。
As shown in FIG. 9, in the transient operation state, at the moment when the throttle is opened, the differential pressure before and after the throttle is large, so the amount of air passing through the throttle flows all at once, and gradually falls to a steady state. The throttle passing air amount Gth in the state can be expressed by this ratio A / ADELAY, and in the steady operation state, this value becomes "1" as shown in the lower graph of FIG. Below, this ratio is
OA ".

【0036】また、スロットルの有効開口面積は、スロ
ットル開度に大きく依存しており、スロットル開度の変
化にほぼ追随して変化する状態となる(図10)。従っ
て、前述した有効開口面積の一次遅れ値を、スロットル
開度の1次遅れ値と等価的に取り扱うこととした。さら
に、チャンバ充填空気量ΔGb が吸入空気量Gc へ反映
される遅れを解消すべく、このスロットル開度の1次遅
れに加え、値ΔGb の1次遅れも考慮することとした。
Further, the effective opening area of the throttle greatly depends on the throttle opening, and changes in a state of almost following the change of the throttle opening (FIG. 10). Therefore, the first-order lag value of the effective opening area described above is treated equivalently to the first-order lag value of the throttle opening. Furthermore, in order to eliminate the delay in which the chamber filling air amount ΔGb is reflected in the intake air amount Gc, in addition to this first-order delay of the throttle opening, the first-order delay of the value ΔGb is considered.

【0037】このようして、チャンバ充填空気量Gb を
スロットル通過空気量Gthから算出することで、スロッ
トル通過空気量Gthを基に吸入空気量Gc を求めること
ができる。これによって、構成が簡易になると共に、演
算量も削減できる。具体的には、単位時間ΔT当たりの
吸入空気量Gc は、数8式のように表すことができる。
また、数9式及び数10式を伝達関数の形式で表すと、
数11式が導かれる。この数11式に示すように、吸入
空気量Gc はスロットル通過空気量Gthの1次遅れ値か
ら求めることができる。
In this way, the intake air amount Gc can be obtained based on the throttle passing air amount Gth by calculating the chamber filling air amount Gb from the throttle passing air amount Gth. This simplifies the configuration and reduces the amount of calculation. Specifically, the intake air amount Gc per unit time ΔT can be expressed by the equation (8).
Further, when the equations 9 and 10 are expressed in the form of transfer function,
Equation 11 is derived. As shown in the equation (11), the intake air amount Gc can be obtained from the primary delay value of the throttle passing air amount Gth.

【0038】[0038]

【数8】 (Equation 8)

【0039】[0039]

【数9】 [Equation 9]

【0040】[0040]

【数10】 [Equation 10]

【0041】[0041]

【数11】 [Equation 11]

【0042】そこで、このような一連の演算処理を、図
11にブロック図として示す。なお、吸入空気量Gc は
燃料噴射量と同様に取り扱うことできるため、以下に示
す各ブロック図等では、便宜上、吸入空気量を燃料噴射
量として取り扱うこととする。また、図中、「(1−
B)/(z−B)」は離散系の伝達関数で1次遅れを意
味する。
Therefore, such a series of arithmetic processing is shown as a block diagram in FIG. Since the intake air amount Gc can be treated in the same manner as the fuel injection amount, the intake air amount is treated as the fuel injection amount in the following block diagrams and the like for convenience. Further, in the figure, "(1-
B) / (z-B) "is a transfer function of a discrete system and means a first-order delay.

【0043】よって、定常運転状態時のある条件下(機
関回転数Ne と吸気圧力Pb によって規定される条件)
でマップ検索により決定される燃料噴射量をTimapと記
すと、実際に出力すべき燃料噴射量Tout は次式より求
めることとした。 Tout =Timap×RATIO-A さらに、シミュレーションを通じて吸入空気の挙動を検
証したところ、チャンバ充填空気量Gb が増加する際と
減少する際では、その挙動が異なっていた。この様子を
図12に示す。図12では、チャンバ充填空気量Gb が
変動する場合の主要因となるスロットル弁の動きに着目
し、このスロットル弁が開き方向に変移する場合及び閉
じ方向に変移する場合における、吸入空気量(スロット
ル通過空気量、気筒吸入空気量、チャンバ充填空気量)
の変化挙動を示している。また、図13に、スロットル
弁の変動に対するチャンバ内の圧力変動を示す。この場
合もスロットル弁が開く方向に変移する場合と閉じる方
向に変移する場合では、チャンバ内の圧力変化の挙動が
相異している。このようにスロットル弁の開き方向と閉
じ方向で吸入空気の挙動が変化するのは、吸入される空
気が圧縮性流体であることがその要因の1つとなってい
ると思われる。
Therefore, under certain conditions during steady operation (conditions defined by the engine speed Ne and the intake pressure Pb)
When the fuel injection amount determined by the map search is described as Timap, the fuel injection amount Tout to be actually output is determined by the following equation. Tout = Timap × RATIO-A Furthermore, when the behavior of the intake air was verified through simulation, the behavior was different when the chamber filling air amount Gb increased and decreased. This state is shown in FIG. In FIG. 12, focusing on the movement of the throttle valve, which is the main factor in the case where the chamber filling air amount Gb fluctuates, the intake air amount (throttle) when the throttle valve moves in the opening direction and in the closing direction is changed. (Passing air amount, cylinder intake air amount, chamber filling air amount)
The change behavior of is shown. Further, FIG. 13 shows the pressure fluctuation in the chamber with respect to the fluctuation of the throttle valve. In this case as well, the behavior of the pressure change in the chamber is different when the throttle valve changes in the opening direction and when it changes in the closing direction. It is considered that one of the factors that the behavior of the intake air changes in the opening direction and the closing direction of the throttle valve is that the intake air is a compressible fluid.

【0044】そこで、本実施例にかかる内燃機関の吸入
空気量推定装置では、より一層正確に気筒吸入空気量を
推定すべく、チャンバ充填空気量Gb が増加する場合と
減少する場合で、スロットルの1次遅れの補正係数Bを
個々に設定することとした。
Therefore, in the intake air amount estimating apparatus for the internal combustion engine according to the present embodiment, in order to more accurately estimate the cylinder intake air amount, the amount of throttle change depending on whether the chamber filling air amount Gb increases or decreases. The correction coefficient B for the primary delay is set individually.

【0045】以下、図14のブロック図及び図15のメ
イン・フローチャートを参照して、この吸入空気量推定
装置の動作について説明する。なお、このフローチャー
トは各TDC位置で起動される。
The operation of the intake air amount estimating apparatus will be described below with reference to the block diagram of FIG. 14 and the main flowchart of FIG. It should be noted that this flowchart is activated at each TDC position.

【0046】まず、各センサで検出した機関回転数N
e、吸気圧力Pb、スロットル開度θTH、気圧Pa、機
関冷却水温Twなどを読み込む(S10)。また、スロ
ットル開度θTHは、アイドル運転状態のスロットル全閉
開度を学習し、その値を用いる。
First, the engine speed N detected by each sensor
e, intake pressure Pb, throttle opening θTH, atmospheric pressure Pa, engine cooling water temperature Tw, etc. are read (S10). Further, as the throttle opening θTH, the throttle fully closed opening in the idle operation state is learned and the value is used.

【0047】次いで、機関がクランキング(始動)中か
否かを判定し(S20)、クランキング中ではないと判
断された場合には(S20で「NO」)、続いてフュー
エル・カットか否かが判断される(S30)。S30で
「NO」と判断された場合には、機関回転数Neと吸気
圧力PbとをもとしてROM58に格納したマップ(図
6)を検索し、定常運転状態時の燃料噴射量Timap(定
常運転状態時の吸入空気量Gc’)を求める(S4
0)。検索した燃料噴射量Timapは、この後、気圧補正
などを必要に応じて適宜加える(図示せず)。
Then, it is judged whether the engine is cranking (starting) or not (S20). If it is judged that the engine is not cranking ("NO" in S20), then it is judged whether or not the fuel is cut. It is determined whether or not (S30). If "NO" is determined in S30, the map (FIG. 6) stored in the ROM 58 including the engine speed Ne and the intake pressure Pb is searched, and the fuel injection amount Timap in the steady operation state (steady operation The intake air amount Gc 'in the state is calculated (S4)
0). Thereafter, the retrieved fuel injection amount Timap is appropriately added with atmospheric pressure correction and the like (not shown).

【0048】次いで、S10で読み込まれたスロットル
開度θTHの1次遅れ値θTH-Dを演算する(S50)。図
16に、ここで行われる演算処理のフローチャートを示
す。なお、このフローチャートも各TDC位置で起動さ
れる。
Next, the first-order delay value θTH-D of the throttle opening θTH read in S10 is calculated (S50). FIG. 16 shows a flowchart of the arithmetic processing performed here. Note that this flowchart is also activated at each TDC position.

【0049】まず、スロットル開度センサ36から得ら
れるスロットル開度θTHに基づき、運転状態が加速状態
か減速状態かを判断する(S51)。この加減速の判断
は、今回のタイミングで得られるスロットル開度をθTH
(k)とし、前回のタイミングで得られるスロットル開
度をθTH(k−1)とすると、θTH(k)とθTH(k−
1)とを比較することにより行っている。すなわち、θ
TH(k)≧θTH(k−1)であれば加速状態と判断し、
θTH(k)<θTH(k−1)であれば減速状態と判断す
る。
First, it is determined whether the operating state is the accelerating state or the decelerating state based on the throttle opening degree θTH obtained from the throttle opening degree sensor 36 (S51). This acceleration / deceleration judgment is based on the throttle opening obtained at this timing by θTH
(K) and the throttle opening obtained at the previous timing is θTH (k-1), θTH (k) and θTH (k-
This is done by comparing with 1). That is, θ
If TH (k) ≧ θTH (k−1), it is judged as an acceleration state,
If θTH (k) <θTH (k-1), it is determined that the vehicle is in the deceleration state.

【0050】このとき、加速状態と判断された場合には
(S51で「YES」)、加速時の補正係数B1 (例え
ば、B1 =0.4)を、スロットル1次遅れ補正係数B
として設定する(S52)。一方、減速状態と判断され
た場合には(S51で「NO」)、減速時の補正係数B
2 (例えば、B2 =0.8)を、スロットル1次遅れ補
正係数Bとして設定する(S53)。
At this time, if it is determined that the vehicle is in the accelerated state ("YES" in S51), the correction coefficient B 1 (for example, B 1 = 0.4) at the time of acceleration is set to the throttle primary delay correction coefficient B.
(S52). On the other hand, when it is determined that the vehicle is in the deceleration state (“NO” in S51), the correction coefficient B during deceleration is set.
2 (for example, B 2 = 0.8) is set as the throttle primary delay correction coefficient B (S53).

【0051】次いで、このように設定されたスロットル
1次遅れ補正係数Bをもちいて、スロットル開度θTHの
1次遅れ伝達関数「(1−B)/(z−B)」が設定さ
れる(S54)。そして、設定された1次遅れ伝達関数
「(1−B)/(z−B)」とスロットル開度θTHとに
より、スロットル開度の1次遅れ値θTH-Dを演算する
(S55)。
Next, using the throttle first-order lag correction coefficient B thus set, the first-order lag transfer function "(1-B) / (z-B)" of the throttle opening θTH is set ( S54). Then, a first-order delay value θTH-D of the throttle opening is calculated from the set first-order delay transfer function "(1-B) / (z-B)" and the throttle opening θTH (S55).

【0052】続いて、図15のメイン・フローチャート
に戻り、スロットル開度θTHと吸気圧力Pbをもとに、
前述した手法によって現在のスロットルの有効開口面積
Aを算出する(S60)。次いで、スロットル開度の1
次遅れ値θTH-Dと吸気圧力Pbによりスロットルの有効
開口面積の1次遅れ値ADELAY を算出する(S70)。
Subsequently, returning to the main flow chart of FIG. 15, based on the throttle opening θTH and the intake pressure Pb,
The current effective opening area A of the throttle is calculated by the method described above (S60). Next, throttle opening 1
The primary delay value ADELAY of the effective opening area of the throttle is calculated from the secondary delay value θTH-D and the intake pressure Pb (S70).

【0053】次いで、「RATIO-A 」を RATIO-A =(A+ABYPASS)/(A+ABYPASS)DELAY なる式から算出する(S80)。なお、ABYPASSは、ス
ロットル弁を通過せずに、バイパス路32などを経由し
て各気筒燃焼室に吸入される空気量(図14に「リフト
量」として示す)を意味し、正確に吸入空気量を決定す
るためにはこの空気量も勘案する必要があるために、こ
の空気量を考慮して演算を行うものである。具体的に
は、この空気量に対応する値を所定の特性に従ってスロ
ットル開度ABYPASSに換算して求めておいて有効開口面
積Aに加算すると共に、その和(A+ABYPASS)とその
1次遅れの近似値「(A+ABYPASS)DELAY 」の比を求
め、それをRATIO-A とする。
Then, "RATIO-A" is calculated from the equation RATIO-A = (A + ABYPASS) / (A + ABYPASS) DELAY (S80). It should be noted that ABYPASS means the amount of air taken into each cylinder combustion chamber (shown as “lift amount” in FIG. 14) via the bypass 32 and the like without passing through the throttle valve, and the intake air is accurately expressed. Since it is necessary to consider this air amount in order to determine the amount, the calculation is performed in consideration of this air amount. Specifically, a value corresponding to this air amount is converted into a throttle opening ABYPASS in accordance with a predetermined characteristic and added to the effective opening area A, and the sum (A + ABYPASS) and its first-order lag are approximated. Find the ratio of the value "(A + ABYPASS) DELAY" and call it RATIO-A.

【0054】このように、分子、分母の双方に加算する
結果、スロットル弁を通過しないで各気筒燃焼室に吸入
される空気量の計測に誤りがあっても、決定される燃料
噴射量への影響度が小さくなる。続いて、S40で検索
して求めたTimap(Gc’)にRATIO-A を乗じ、スロッ
トル通過空気量に相当する燃料噴射量Tout (Gc)を
算出する(S90)。
As described above, as a result of adding to both the numerator and the denominator, even if there is an error in the measurement of the amount of air taken into each cylinder combustion chamber without passing through the throttle valve, The impact is small. Then, RATIO-A is multiplied by the Timap (Gc ') obtained by the search in S40 to calculate the fuel injection amount Tout (Gc) corresponding to the throttle passing air amount (S90).

【0055】なお、S20においてクランキング中と判
断された場合には(S20で「YES」)、水温Twか
ら所定のテーブル(図示省略)を検索してクランキング
時の燃料噴射量Ticr を算出し(S21)、次いで始動
モードの式(説明省略)に基づいて燃料噴射量Tout を
決定する(S22)。
When it is determined in S20 that cranking is in progress ("YES" in S20), a predetermined table (not shown) is searched from the water temperature Tw to calculate the fuel injection amount Ticr during cranking. (S21), and then, the fuel injection amount Tout is determined based on the equation (the description is omitted) of the start mode (S22).

【0056】一方、S30においてフューエル・カット
と判断された場合には(S30で「YES」)、燃料噴
射量Tout を零に設定する(S31)。
On the other hand, if the fuel cut is determined in S30 ("YES" in S30), the fuel injection amount Tout is set to zero (S31).

【0057】このように、スロットル1次遅れ補正係数
Bを加速時と減速時とで別持ちにした場合には、スロッ
トル1次遅れ補正係数Bを一定とした場合に較べ、より
正確に吸入空気量が推定できるのかを評価した。この評
価は、機関回転数、吸気圧力等を一定として実施した。
As described above, when the throttle primary delay correction coefficient B is separately provided for acceleration and deceleration, the intake air is more accurate than when the throttle primary delay correction coefficient B is constant. It was evaluated whether the quantity could be estimated. This evaluation was carried out with the engine speed, intake pressure, etc. kept constant.

【0058】図17(a)〜(c)にその評価結果を示
す。図17(a)はスロットル開度θTHの変移状態を示
しており、スロットルを所定の開度まで急開させ(加速
状態)、その後、所定の開度まで急閉した状態(減速状
態)を示している。このようにスロットルを変移させた
状態では、スロットル1次遅れ補正係数Bを一定とした
場合(図17(b))、スロットルの開き側に比べて閉
じ側の空燃比の挙動に大きなバラツキが生じているのが
確認できる。一方、スロットルの開き側と閉じ側とにお
いて、スロットル1次遅れ補正係数Bをそれぞれ適正な
値に設定した場合には(図17(c))、空燃比のバラ
ツキが図17(b)に比べて低減されており、全期間に
渡って、空燃比がほぼフラットに推移していることが分
る。従って、スロットルの開き側と閉じ側とにおいて補
正係数Bを個々に適正な値に設定することで、吸入空気
量Gc (燃料噴射量Tout )をより正確に推定できるこ
とが確認できた。なお、図17(b)、(c)における
縦軸は、当量比、すなわち、Mst/Mを示している
(Mst:理論空燃比、M:空気消費量/燃料消費
量)。
17 (a) to 17 (c) show the evaluation results. FIG. 17A shows a transition state of the throttle opening θTH, showing a state in which the throttle is rapidly opened to a predetermined opening (acceleration state) and then rapidly closed to a predetermined opening (deceleration state). ing. In the state where the throttle is displaced in this way, when the throttle first-order lag correction coefficient B is kept constant (FIG. 17 (b)), there is a large variation in the behavior of the air-fuel ratio on the closed side compared to the open side of the throttle. Can be confirmed. On the other hand, when the throttle first-order lag correction coefficient B is set to an appropriate value on each of the opening side and the closing side of the throttle (Fig. 17 (c)), the variation in the air-fuel ratio is larger than that in Fig. 17 (b). It can be seen that the air-fuel ratio has remained almost flat over the entire period. Therefore, it has been confirmed that the intake air amount Gc (fuel injection amount Tout) can be more accurately estimated by setting the correction coefficient B to an appropriate value individually on the opening side and the closing side of the throttle. The vertical axes in FIGS. 17B and 17C represent the equivalence ratio, that is, Mst / M (Mst: theoretical air-fuel ratio, M: air consumption / fuel consumption).

【0059】ここで、図18を参照して他の実施例を説
明する。前述した実施例では、チャンバ充填空気量Gb
が増減する状況を、その主要因となるスロットル開度の
みに着目して検知したが、この他に以下の要因によって
もチャンバ充填空気量Gb が増減するため、これらを総
合的に勘案してチャンバ充填空気量Gb の増減を検知す
ることとした。すなわち、機関回転数が変化した場合、
インテークマニホルドの形態を機械的に変化させてその
容積を増減させた場合、大気圧が変化した場合、吸気弁
・排気弁の開閉タイミング(バルブタイミング)が変化
した場合、目標となる空燃比が変化した場合などにも、
チャンバ充填空気量Gb が実質的に変化する。
Now, another embodiment will be described with reference to FIG. In the above-described embodiment, the chamber filling air amount Gb
The situation in which the increase / decrease is detected by focusing on only the throttle opening, which is the main factor, is detected. However, the chamber filling air amount Gb also increases / decreases due to the following factors. It was decided to detect an increase / decrease in the filled air amount Gb. That is, when the engine speed changes,
The target air-fuel ratio changes when the form of the intake manifold is mechanically changed to increase or decrease its volume, when atmospheric pressure changes, when the intake / exhaust valve opening / closing timing (valve timing) changes. If you do,
The chamber filling air amount Gb changes substantially.

【0060】そこで、本実施例では、これらの要因を総
合的に勘案してチャンバ充填空気量Gb の増減を検知
し、この結果を基に、スロットル開度の1次遅れ値θTH
-Dを演算する場合を例示する。
Therefore, in the present embodiment, the increase and decrease in the chamber filling air amount Gb is detected by comprehensively considering these factors, and based on this result, the first-order lag value θTH of the throttle opening is detected.
An example of calculating -D is shown.

【0061】図18に、スロットル開度の1次遅れ値θ
TH-Dの演算処理を実行するフローチャートを示す。な
お、このフローチャートは、図15に示すメイン・フロ
ーチャートのうち、S50において行われる演算処置を
示すものであり、各TDC位置で起動される。
FIG. 18 shows the first-order delay value θ of the throttle opening.
The flowchart which performs the calculation process of TH-D is shown. It should be noted that this flowchart shows the arithmetic processing performed in S50 of the main flowchart shown in FIG. 15, and is started at each TDC position.

【0062】まず、バルブタイミングV/Tが「Hig
h」か「Low」かを判断し(S100)、いずれの場
合においも、図16の処理と同様に、スロットル開度θ
THの変化状態から、運転状態が加速状態か減速状態かを
判断する(S101、S102)。
First, the valve timing V / T is "High".
h ”or“ Low ”(S100), and in any case, the throttle opening θ is set in the same manner as the processing of FIG.
From the change state of TH, it is determined whether the operating state is the acceleration state or the deceleration state (S101, S102).

【0063】S101において、加速状態と判断した場
合には(S101で「YES」)、加速時の補正係数B
3 (例えば、B3 =0.4)をスロットル1次遅れ補正
係数Bとして設定し(S103)、減速状態と判断され
た場合には(S101で「NO」)、減速時の補正係数
4 (例えば、B4 =0.8)を、スロットル1次遅れ
補正係数Bとして設定する(S104)。
When it is determined in S101 that the vehicle is in the accelerated state ("YES" in S101), the correction coefficient B during acceleration is set.
3 (for example, B 3 = 0.4) is set as the throttle first-order lag correction coefficient B (S103), and when it is determined that the vehicle is in the deceleration state (“NO” in S101), the correction coefficient B 4 at the time of deceleration is set. (For example, B 4 = 0.8) is set as the throttle primary delay correction coefficient B (S104).

【0064】同様に、S102において、加速状態と判
断された場合には(S102で「YES」)、加速時の
補正係数B5 (例えば、B5 =0.3)をスロットル1
次遅れ補正係数Bとして設定し(S105)、減速状態
と判断された場合には(S102で「NO」)、減速時
の補正係数B6 (例えば、B6 =0.6)を、スロット
ル1次遅れ補正係数Bとして設定する(S106)。こ
のように、バルブタイミングV/Tが「High」のと
き、「Low」の場合に比べ、スロットル1次遅れ補正
係数Bが小さな値となる傾向がある。また、いずれの場
合も、加速時の値が減速時の値に比べて小さな値となる
傾向がある。
[0064] Similarly, in S102, the ( "YES" in S102) if it is determined that the acceleration state, the correction coefficient during acceleration B 5 (e.g., B 5 = 0.3) the throttle 1
Set the next delay correction coefficient B (S105), when it is determined that the deceleration state ( "NO" in S102), the correction coefficient of deceleration B 6 (e.g., B 6 = 0.6), throttle 1 The next delay correction coefficient B is set (S106). As described above, when the valve timing V / T is “High”, the throttle first-order delay correction coefficient B tends to have a smaller value than when the valve timing V / T is “Low”. In any case, the value during acceleration tends to be smaller than the value during deceleration.

【0065】続いて、S103〜S106において設定
したスロットル1次遅れ補正係数Bを、機関回転数、大
気圧などに応じてさらに補正する。この場合、まず、図
19(a)〜(d)に示すように、機関回転数Ne、大
気圧PA、吸気空気温度TA及び目標空燃比AFに応じ
た各補正係数KADNe,KADPA,KADTA及び
KADAFを予め求めておき、現在の機関回転数等に応
じて、対応する各種補正係数KADNe,KADPA,
KADTA及びKADAFを検索して求める(S10
7)。次いで、S103〜S106のいずれかで設定し
たスロットル1次遅れ補正係数Bに対して、求めた各種
補正係数KADNe,KADPA,KADTA及びKA
DAFをそれぞれ乗算し、この値を最終的なスロットル
1次遅れ補正係数Bとして設定する(S108)。
Subsequently, the throttle first-order delay correction coefficient B set in S103 to S106 is further corrected according to the engine speed, atmospheric pressure, and the like. In this case, first, as shown in FIGS. 19A to 19D, the correction coefficients KADNe, KADPA, KADATA, and KADAF according to the engine speed Ne, the atmospheric pressure PA, the intake air temperature TA, and the target air-fuel ratio AF. Is calculated in advance, and various correction coefficients KADNe, KADPA, corresponding to the current engine speed, etc. are obtained.
Search KADA and KADAF to obtain (S10
7). Next, with respect to the throttle first-order lag correction coefficient B set in any of S103 to S106, various calculated correction coefficients KADNe, KADPA, KATTA and KA
Multiply each by DAF and set this value as the final throttle primary delay correction coefficient B (S108).

【0066】次いで、このように設定されたスロットル
1次遅れ補正係数Bをもちいて、スロットル開度θTHの
1次遅れ伝達関数「(1−B)/(z−B)」が設定さ
れ(S109)、設定された1次遅れ伝達関数「(1−
B)/(z−B)」とスロットル開度θTHとにより、ス
ロットル開度の1次遅れ値θTH-Dを演算する(S11
0)。
Next, using the throttle first-order lag correction coefficient B thus set, the first-order lag transfer function "(1-B) / (z-B)" of the throttle opening θTH is set (S109). ), The set first-order lag transfer function "(1-
B) / (z−B) ”and the throttle opening θTH, a first-order delay value θTH-D of the throttle opening is calculated (S11).
0).

【0067】このようにして求めたスロットル開度の1
次遅れ値θTH-Dを用いて、図15のメイン・フローチャ
ートに従って燃料噴射量Tout を算出する。
1 of the throttle opening thus obtained
Using the next delay value θTH-D, the fuel injection amount Tout is calculated according to the main flowchart of FIG.

【0068】なお、図18に示した実施例では、フロー
チャートのS100においてバルブタイミングV/Tが
「High」か「Low」かを判断するとして例示した
が、このS100において、インテークマニホルドの容
積変化を検知してもよい。この場合には、インテークマ
ニホルドの形態を機械的に変化させる制御が伴うため、
この制御信号を検知すればよい。この結果、インテーク
マニホルドの容積が増加する場合には、応答遅れが大き
くなることから、図18のフローチャートのS101に
進み、スロットル1次遅れ補正係数Bは大きな値が設定
される。また、減少する場合には、応答遅れが小さくな
ることから、フローチャートのS102に進み、スロッ
トル1次遅れ補正係数Bは小さな値が設定される。
In the embodiment shown in FIG. 18, it is exemplified that the valve timing V / T is determined to be "High" or "Low" in S100 of the flow chart, but in this S100, the change in the intake manifold volume is determined. It may be detected. In this case, because the control that mechanically changes the form of the intake manifold is involved,
It suffices to detect this control signal. As a result, when the intake manifold volume increases, the response delay increases, so the flow proceeds to S101 in the flowchart of FIG. 18, and the throttle first-order delay correction coefficient B is set to a large value. Further, if it decreases, the response delay becomes small, so the flow proceeds to S102 of the flowchart, and the throttle primary delay correction coefficient B is set to a small value.

【0069】上記した燃料噴射量Tout の各算出手法
は、簡易なアルゴリズムによって定常運転状態から過渡
運転状態までを表現することができ、定常運転状態時の
燃料噴射量をマップ検索によってある程度保証すること
ができると同時に、複雑な演算を必要とせずに燃料噴射
量を最適に決定することができる。しかも、増加時と減
少時で相異する吸入空気の挙動に応じて、スロットル開
度の1次遅れ補正係数を設定するため、より正確に吸入
空気量を推定することができ、制御性や制御制度をより
一層向上させることができる。
Each of the above-described methods for calculating the fuel injection amount Tout can represent from the steady operation state to the transient operation state by a simple algorithm, and the fuel injection amount in the steady operation state can be guaranteed to some extent by map search. At the same time, the fuel injection amount can be optimally determined without requiring complicated calculation. Moreover, since the first-order lag correction coefficient of the throttle opening is set according to the behavior of the intake air that differs between the increase and the decrease, the intake air amount can be estimated more accurately, and controllability and control The system can be further improved.

【0070】また、以上説明した実施例では、燃料噴射
量Timapを予めマップ化する例を示したが、これに代え
て、定常運転状態時における吸入空気量Gc ’をマップ
化し、ROM58に格納しておいてもよい。
In the embodiment described above, the fuel injection amount Timap is mapped in advance. Instead of this, the intake air amount Gc 'in the steady operation state is mapped and stored in the ROM 58. You may keep it.

【0071】さらに、例示した実施例では、推定された
吸入空気量に対応して燃料噴射量を決定することとした
が、この例に限定されるものではなく、この他にも、推
定された吸入空気量に対応し、例えば、点火時期、排気
還流量(EGR量)など、他のエンジン制御パラメータ
を算出することもできる。
Further, in the illustrated embodiment, the fuel injection amount is determined according to the estimated intake air amount, but the present invention is not limited to this example, and other estimations are made. Other engine control parameters, such as ignition timing and exhaust gas recirculation amount (EGR amount), can be calculated corresponding to the intake air amount.

【0072】[0072]

【発明の効果】以上説明したように、請求項1にかかる
内燃機関の吸入空気量推定装置では、検知手段を備えた
ので、充填空気量Gb の増加、減少を検知でき、さら
に、設定手段によって、この検知結果もとにスロットル
開度の1次遅れ補正係数を個々に設定できる。従って、
充填空気量Gb の増加時と減少時とで相異する吸入空気
の挙動に対応した吸入空気量の推定処理を実施すること
ができる。このため、スロットルを通過する空気量に基
づき気筒内に流入する空気量を、従来よりも一層正確に
推定することが可能となり、内燃機関の制御性をさらに
向上させることができる。
As described above, the intake air amount estimating apparatus for an internal combustion engine according to the first aspect of the invention is provided with the detecting means, so that it is possible to detect the increase or decrease of the filling air amount Gb, and further by the setting means. The first-order delay correction coefficient of the throttle opening can be individually set based on the detection result. Therefore,
It is possible to carry out the intake air amount estimation processing corresponding to the behavior of the intake air that differs when the filled air amount Gb increases and when it decreases. Therefore, it is possible to more accurately estimate the amount of air flowing into the cylinder based on the amount of air passing through the throttle, and it is possible to further improve the controllability of the internal combustion engine.

【0073】また、請求項2にかかる内燃機関の吸入空
気量推定装置では、この検知手段を、スロットル開度θ
THの増加、減少に基づいて、充填空気量Gb の増加、減
少を検知する手段として構成したので、充填空気量Gb
が変動する場合の主要因となるスロットル弁の動きに基
づいて吸入空気量の推定処理を実施することができる。
Further, in the intake air amount estimating apparatus for the internal combustion engine according to the second aspect, the detecting means is connected to the throttle opening θ.
Since it is configured as a means for detecting the increase / decrease of the filling air amount Gb based on the increase / decrease of TH, the filling air amount Gb
The intake air amount estimation processing can be performed based on the movement of the throttle valve, which is the main factor in the case of the fluctuation of the.

【0074】また、請求項3にかかる内燃機関の吸入空
気量推定装置では、検知手段により、大気圧の変化など
充填空気量Gb の増減に寄与する他の要因についても検
知することとしたので、吸入空気量をより一層正確に推
定することが可能となる。
Further, in the intake air amount estimating apparatus for the internal combustion engine according to the third aspect, the detecting means also detects other factors contributing to the increase and decrease of the filling air amount Gb such as the change of the atmospheric pressure. It is possible to more accurately estimate the intake air amount.

【図面の簡単な説明】[Brief description of drawings]

【図1】本実施例にかかる内燃機関及び各種センサの配
設位置等を示す概略構成図である。
FIG. 1 is a schematic configuration diagram showing an arrangement position and the like of an internal combustion engine and various sensors according to this embodiment.

【図2】内燃機関の吸入空気量推定装置における制御装
置の構成を示すブロック図である。
FIG. 2 is a block diagram showing a configuration of a control device in an intake air amount estimation device for an internal combustion engine.

【図3】本実施例で採用した流体力学モデルを示す説明
図である。
FIG. 3 is an explanatory diagram showing a fluid dynamics model adopted in this embodiment.

【図4】図3の流体力学モデルを用いて、スロットル弁
の有効開口面積を流量係数などを用いて算出する手法を
示すブロック図である。
FIG. 4 is a block diagram showing a method of calculating an effective opening area of a throttle valve using a flow coefficient and the like using the fluid dynamics model of FIG.

【図5】図4の算出で用いる係数のマップ特性を示す説
明図である。
5 is an explanatory diagram showing map characteristics of coefficients used in the calculation of FIG.

【図6】定常運転状態時の燃料噴射量のマップ特性を示
す説明図である。
FIG. 6 is an explanatory diagram showing map characteristics of a fuel injection amount in a steady operation state.

【図7】目標空燃比のマップを示す説明図である。FIG. 7 is an explanatory diagram showing a map of a target air-fuel ratio.

【図8】スロットルの有効開口面積についてのシミュレ
ーション結果を示すデータ図である。
FIG. 8 is a data diagram showing a simulation result of an effective opening area of a throttle.

【図9】本実施例における定常運転状態と過渡運転状態
とを示す説明図である。
FIG. 9 is an explanatory diagram showing a steady operation state and a transient operation state in the present embodiment.

【図10】スロットル開度とスロットルの有効開口面積
との関係を示す説明図である。
FIG. 10 is an explanatory diagram showing the relationship between the throttle opening and the effective opening area of the throttle.

【図11】燃料噴射量(吸入空気量)の算出処理を示す
ブロック図である。
FIG. 11 is a block diagram showing a process of calculating a fuel injection amount (intake air amount).

【図12】加速・減速時における吸入空気の挙動を示す
データ図である。
FIG. 12 is a data diagram showing the behavior of intake air during acceleration / deceleration.

【図13】スロットル弁の開き側と閉じ側でチャンバ内
圧力の挙動が相異することを示すデータ図である。
FIG. 13 is a data diagram showing that the behavior of the pressure in the chamber is different between the opening side and the closing side of the throttle valve.

【図14】運転状態に応じスロットルの1次遅れの補正
係数を個々に設定する場合における、燃料噴射量(吸入
空気量)の算出処理を示すブロック図である。
FIG. 14 is a block diagram showing a calculation process of a fuel injection amount (intake air amount) in the case of individually setting a correction coefficient of the primary delay of the throttle according to an operating state.

【図15】燃料噴射量(吸入空気量)の算出処理を示す
メイン・フローチャートである。
FIG. 15 is a main flowchart showing a process of calculating a fuel injection amount (intake air amount).

【図16】図15のメイン・フローチャートのうち、θ
THの1次遅れ値θTH-Dの算出処置のみを示すフローチャ
ートである。
16 is a flowchart of the main flowchart of FIG.
It is a flow chart which shows only the calculation procedure of the 1st-order lag value THH-D of TH.

【図17】(a)はスロットル開度の変移状態を示すデ
ータ図である。(b)はスロットル1次遅れ補正係数B
を一定とした場合において、(a)のようにスロットル
が変移した状態に対応する空燃比の推移を示すデータ図
である。(c)はスロットル1次遅れ補正係数Bをスロ
ットルの開き側と閉じ側とで別持ちとした場合におい
て、(a)のようにスロットルが変移した状態に対応す
る空燃比の推移を示すデータ図である。
FIG. 17A is a data diagram showing a change state of the throttle opening. (B) is the throttle first-order lag correction coefficient B
FIG. 6 is a data diagram showing a transition of the air-fuel ratio corresponding to a state where the throttle has changed as shown in FIG. (C) is a data diagram showing the transition of the air-fuel ratio corresponding to the state where the throttle has changed as shown in (a) when the throttle first-order lag correction coefficient B is separately held on the opening side and the closing side of the throttle. Is.

【図18】スロットル開度の1次遅れ値の算出処置に関
する、他の実施例を示すフローチャートである。
FIG. 18 is a flowchart showing another embodiment regarding a procedure for calculating the first-order lag value of the throttle opening.

【図19】(a)は機関回転数と補正係数KADNeと
の関係を示すグラフ、(b)は大気圧と補正係数KAD
PAとの関係を示すグラフ、(c)は吸入空気温度と補
正係数KADTAとの関係を示すグラフ、(d)は目標
空燃比と補正係数KADFAとの関係を示すグラフであ
る。
19 (a) is a graph showing the relationship between the engine speed and the correction coefficient KADNe, and FIG. 19 (b) is the atmospheric pressure and the correction coefficient KAD.
6 is a graph showing the relationship with PA, (c) is a graph showing the relationship between the intake air temperature and the correction coefficient KADTA, and (d) is a graph showing the relationship between the target air-fuel ratio and the correction coefficient KADFA.

【符号の説明】[Explanation of symbols]

10…内燃機関、12…吸気路、16…スロットル弁、
20…インテークマニホルド、36…スロットル開度セ
ンサ。
10 ... Internal combustion engine, 12 ... Intake passage, 16 ... Throttle valve,
20 ... Intake manifold, 36 ... Throttle opening sensor.

Claims (3)

【特許請求の範囲】[Claims] 【請求項1】機関回転数Ne 、スロットル開度θTH及び
チャンバ内圧力Pb を含む、内燃機関の運転状態を検出
する第1の手段と、 少なくとも、前記第1の手段で検出された前記機関回転
数Ne と前記チャンバ内圧力Pb とを基に、定常運転状
態時における燃焼室への吸入空気量Gc ’を求める第2
の手段と、 少なくとも、前記第1の手段で検出された前記スロット
ル開度θTHと前記チャンバ内圧力Pb とを基に、スロッ
トルの第1有効開口面積Aを求める第3の手段と、 前記スロットルの第1有効開口面積Aの1次遅れ値を求
め、その値をスロットルの第2有効開口面積ADELAY と
する第4の手段と、 前記第3の手段で求めた第1有効開口面積Aと前記第4
の手段で求めた第2有効開口面積ADELAY とを基に、前
記吸入空気量Gc ’を補正することで、実吸入空気量G
c を、 Gc =Gc ’×(A/ADELAY ) として求める第5の手段とを備え、 スロットルから燃焼室に至る吸気系をモデル化して1つ
のチャンバとみなし、この吸気系モデルに基づき、この
燃焼室へ吸入される実吸入空気量Gc を推定する内燃機
関の吸入空気量推定装置において、 前記第4の手段は、 前記チャンバ内に充填される充填空気量Gb の増加、減
少を検知する検知手段と、 前記検知手段の検知結果を受けて、前記充填空気量Gb
が増加した場合と減少した場合とで、前記スロットル開
度の1次遅れ補正係数Bを個々に設定する設定手段と、 前記設定手段で設定された1次遅れ補正係数Bを基に、
前記スロットルの第1有効開口面積Aの1次遅れ値を求
める演算手段と、 を備えることを特徴とする内燃機関の吸入空気量推定装
置。
1. A first means for detecting an operating state of an internal combustion engine including an engine speed Ne, a throttle opening θTH and a chamber pressure Pb, and at least the engine speed detected by the first means. A second method for obtaining the intake air amount Gc 'into the combustion chamber in the steady operation state based on the number Ne and the chamber pressure Pb
Means for determining the first effective opening area A of the throttle based on at least the throttle opening θTH and the chamber pressure Pb detected by the first means; A fourth means for obtaining a first-order lag value of the first effective opening area A and setting the value as a second effective opening area ADELAY of the throttle; a first effective opening area A and a first effective opening area A obtained by the third means; Four
By correcting the intake air amount Gc 'based on the second effective opening area ADELAY obtained by the means described above, the actual intake air amount G
and a fifth means for obtaining c as Gc = Gc '× (A / ADELAY), and the intake system from the throttle to the combustion chamber is modeled as one chamber, and based on this intake system model, In an intake air amount estimation device for an internal combustion engine for estimating an actual intake air amount Gc to be sucked into a chamber, the fourth means is a detection means for detecting an increase or decrease of a filling air amount Gb filled in the chamber. And receiving the detection result of the detection means, the filling air amount Gb
Based on the setting means for individually setting the first-order lag correction coefficient B of the throttle opening depending on whether it increases or decreases, and the first-order lag correction coefficient B set by the setting means.
An intake air amount estimating apparatus for an internal combustion engine, comprising: a calculating unit that obtains a first-order lag value of the first effective opening area A of the throttle.
【請求項2】前記検知手段は、前記充填空気量Gb の増
加、減少を、前記スロットル開度θTHの増加及び減少に
基づいて検知することを特徴とする請求項1記載の内燃
機関の吸入空気量推定装置。
2. The intake air for an internal combustion engine according to claim 1, wherein the detecting means detects an increase or a decrease in the filled air amount Gb based on an increase or a decrease in the throttle opening θTH. Quantity estimation device.
【請求項3】前記検知手段は、インテークマニホルドの
容積変化、内燃機関のバルブタイミングの変化、大気圧
の変化及び目標空燃比の変化のうち、少なくともいずれ
か1つの変化に基づいて、前記充填空気量Gb の増加、
減少を検知することを特徴とする請求項1又は2記載の
内燃機関の吸入空気量推定装置。
3. The charged air is detected based on at least one of a change in intake manifold volume, a change in internal combustion engine valve timing, a change in atmospheric pressure, and a change in target air-fuel ratio. Increase in quantity Gb,
The intake air amount estimating apparatus for an internal combustion engine according to claim 1 or 2, wherein a decrease is detected.
JP03885395A 1995-02-27 1995-02-27 Apparatus for estimating intake air amount of internal combustion engine Expired - Fee Related JP3337338B2 (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP03885395A JP3337338B2 (en) 1995-02-27 1995-02-27 Apparatus for estimating intake air amount of internal combustion engine
US08/607,866 US5597951A (en) 1995-02-27 1996-02-27 Intake air amount-estimating apparatus for internal combustion engines

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP03885395A JP3337338B2 (en) 1995-02-27 1995-02-27 Apparatus for estimating intake air amount of internal combustion engine

Publications (2)

Publication Number Publication Date
JPH08232751A true JPH08232751A (en) 1996-09-10
JP3337338B2 JP3337338B2 (en) 2002-10-21

Family

ID=12536770

Family Applications (1)

Application Number Title Priority Date Filing Date
JP03885395A Expired - Fee Related JP3337338B2 (en) 1995-02-27 1995-02-27 Apparatus for estimating intake air amount of internal combustion engine

Country Status (1)

Country Link
JP (1) JP3337338B2 (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002309993A (en) * 2001-04-13 2002-10-23 Denso Corp Control device for internal combustion engine
JP2003314347A (en) * 2002-04-18 2003-11-06 Denso Corp Device for detecting cylinder filling air amount of internal combustion engine
JP2004211590A (en) * 2002-12-27 2004-07-29 Toyota Motor Corp Device for estimating air intake volume of internal combustion engine
JP2004239173A (en) * 2003-02-06 2004-08-26 Osaka Gas Co Ltd System for controlling combustion of internal combustion engine and method for controlling combustion
JP2005194936A (en) * 2004-01-07 2005-07-21 Toyota Motor Corp State quantity estimating device for internal combustion engine
US7181336B2 (en) 2003-08-26 2007-02-20 Toyota Jidosha Kabushiki Kaisha Control system of internal combustion engine
WO2018116677A1 (en) * 2016-12-19 2018-06-28 日立オートモティブシステムズ株式会社 Control device for internal combustion engine

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002309993A (en) * 2001-04-13 2002-10-23 Denso Corp Control device for internal combustion engine
JP2003314347A (en) * 2002-04-18 2003-11-06 Denso Corp Device for detecting cylinder filling air amount of internal combustion engine
JP2004211590A (en) * 2002-12-27 2004-07-29 Toyota Motor Corp Device for estimating air intake volume of internal combustion engine
JP2004239173A (en) * 2003-02-06 2004-08-26 Osaka Gas Co Ltd System for controlling combustion of internal combustion engine and method for controlling combustion
US7181336B2 (en) 2003-08-26 2007-02-20 Toyota Jidosha Kabushiki Kaisha Control system of internal combustion engine
CN100455787C (en) * 2003-08-26 2009-01-28 丰田自动车株式会社 Control device of internal combustion engine
JP2005194936A (en) * 2004-01-07 2005-07-21 Toyota Motor Corp State quantity estimating device for internal combustion engine
WO2018116677A1 (en) * 2016-12-19 2018-06-28 日立オートモティブシステムズ株式会社 Control device for internal combustion engine

Also Published As

Publication number Publication date
JP3337338B2 (en) 2002-10-21

Similar Documents

Publication Publication Date Title
JP4144272B2 (en) Fuel injection amount control device for internal combustion engine
JP3162553B2 (en) Air-fuel ratio feedback control device for internal combustion engine
US6662640B2 (en) Air amount detector for internal combustion engine
JP3330234B2 (en) Fuel injection control device for internal combustion engine
JP3805840B2 (en) Engine control device
US5597951A (en) Intake air amount-estimating apparatus for internal combustion engines
JP2008138630A (en) Control device of internal combustion engine
JP2005188503A (en) Air-fuel ratio control system by cylinder for internal combustion engine
JP3354304B2 (en) Fuel injection control device for internal combustion engine
JPH05306643A (en) Fuel injection amount calculation device of internal combustion engine
US6591822B2 (en) Air-fuel ratio controller of internal combustion engines
JPH08232751A (en) Suction air amount estimating device for internal combustion engine
JP3559338B2 (en) Apparatus for estimating intake air amount of internal combustion engine
JPH08121211A (en) Fuel control device for internal combustion engine
JP3337339B2 (en) Apparatus for estimating intake air amount of internal combustion engine
CN108999709B (en) Method for calculating the charge of an internal combustion engine
JP3052756B2 (en) Control device for internal combustion engine
JP4025977B2 (en) Engine intake air amount calculation device
JP2005002932A (en) Control system of internal combustion engine
JP4385542B2 (en) Air-fuel ratio control device for internal combustion engine
JP4241560B2 (en) Intake air amount estimation device for internal combustion engine
JP2730760B2 (en) Engine intake air amount detection device
JP2009203868A (en) Intake control device for internal combustion engine
JP2007192091A (en) Control device for internal combustion engine
JP3119465B2 (en) Fuel injection control device for internal combustion engine

Legal Events

Date Code Title Description
LAPS Cancellation because of no payment of annual fees