JP3621731B2 - エンジンの空燃比制御方法 - Google Patents
エンジンの空燃比制御方法 Download PDFInfo
- Publication number
- JP3621731B2 JP3621731B2 JP28097194A JP28097194A JP3621731B2 JP 3621731 B2 JP3621731 B2 JP 3621731B2 JP 28097194 A JP28097194 A JP 28097194A JP 28097194 A JP28097194 A JP 28097194A JP 3621731 B2 JP3621731 B2 JP 3621731B2
- Authority
- JP
- Japan
- Prior art keywords
- fuel
- intake
- cylinder
- mass
- intake air
- 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.)
- Expired - Fee Related
Links
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/12—Improving ICE efficiencies
Landscapes
- Output Control And Ontrol Of Special Type Engine (AREA)
- Electrical Control Of Air Or Fuel Supplied To Internal-Combustion Engine (AREA)
- Combined Controls Of Internal Combustion Engines (AREA)
Description
【０００１】
【産業上の利用分野】
本発明は、過渡時の燃料付着量を補填する燃料追加質量を、筒内に実際に吸入される燃料質量に加算して、燃料噴射の際に必要な燃料噴射質量を求めるエンジンの空燃比制御方法に関する。
【０００２】
【従来の技術】
一般に、インジェクタから噴射された燃料は、その全てが筒内に吸込まれることはなく、一部は吸気ポート内壁等に付着して燃料液膜となる。定常的には燃料付着率が一定であれば、インジェクタからの燃料噴射量と筒内に吸込まれる燃料量とはほぼ一致するが、過渡的にはインジェクタからの噴射量が急増しても、燃料付着量が相対的に増加する分、筒内へ吸入される燃料量は要求燃料量よりも少なくなる。そのため、筒内混合気の空燃比が一時的にリーン化し、上記吸気ポート内壁等に付着した燃料は、ある遅れを有して筒内に吸込まれるため、筒内への吸入燃料量はゆっくりと増加することになる。
【０００３】
従って、エンジンの空燃比制御では、このような過渡的な空燃比のリーン化を防止し、安定した出力を得るために上記燃料付着量分を補填するフィードフォワード制御が必要となる。
【０００４】
例えば、特開平５−１８２８９号公報では、理論的な燃料付着モデルに基づき、まず、１回転当りの吸入空気量Ｑ／Ｎとエンジン回転数Ｎｅとによりマップを参照して暖機完了後の飽和状態の壁面付着量ＱＭＷを補間計算し、前回求めた壁面付着量ＱＭＷとの差分から、変化量ΔＱＭＷを算出し、この変化量ΔＱＭＷを、即時に変化する即時補正と、徐々に変化するテーリング補正とで補正して、壁面付着補正量ＦＭＷを求め、この壁面付着補正量ＦＭＷを、吸入空気量Ｑとエンジン回転数Ｎｅとに基づいて求めた基本燃料噴射量ＴＡＵＰに加算し、必要な補正項で補正して、インジェクタからの噴射量ＴＡＵを算出している。
【０００５】
上記変化量ΔＱＭＷを、即時補正とテーリング補正とに分けて補正しているのは、過渡的な上記壁面付着量ＱＭＷは、計算上は瞬時に変化するが、実際には徐々に変化するためである。
【０００６】
【発明が解決しようとする課題】
上述したように、一般的には、理論的な燃料付着モデルから壁面付着補正量ＦＭＷを求める場合が多いが、吸入空気量Ｑを計測する吸入空気量センサ自体の応答遅れや、ノイズの影響等で壁面付着量ＱＭＷの検出性能にばらつきが生じ易く、充分な信頼性が得られない。
【０００７】
また、古くはスロットル弁の開度変化から過渡的な付着補正を試みる、いわゆる三角増量法がある。この方式は、スロットル弁の開度変化に基づくため応答性は良いが、燃料付着による遅れを三角形で近似的に求めているため、空燃比が部分的にリーン化したり、リッチ化したりし易く、空燃比制御精度に問題がある。
【０００８】
本発明は、上記事情に鑑みてなされたもので、応答遅れや、ノイズの影響を受けること無く、高精度の付着補正量を求めることができて、信頼性の高い過渡補正を行うことのできるエンジンの空燃比制御方法を提供することを目的としている。
【０００９】
【課題を解決するための手段】
上記目的を達成するため本発明は、筒内への１吸気行程当りの吸入空気質量を基本として目標空燃比に対応する筒内への１吸気行程当りの吸入燃料質量を設定するエンジンの空燃比制御方法において、スロットル開度とエンジン回転数とに基づき静的な筒内吸入空気質量に相当する静的指標値を設定し、上記静的指標値に吸気チャンバによる遅れを想定した一次遅れ処理を施して動的な筒内吸入空気質量に相当する動的指標値を設定し、更に上記動的指標値を加重平均して燃料付着による遅れ指標値を算出し、上記動的指標値から上記遅れ指標値を減算して、この差分に基づいて燃料追加質量を設定し、上記燃料追加質量を上記吸入燃料質量に加算して１サイクル当りの燃料噴射量を設定することを特徴とする。
【００１０】
【作用】
本発明は、スロットル開度とエンジン回転数とに基づき静的な筒内吸入空気質量に相当する静的指標値を設定し、静的指標値に吸気チャンバによる遅れを想定した一次遅れ処理を施して動的な筒内吸入空気質量に相当する動的指標値を設定する。そして、更に動的指標値を加重平均して燃料付着による遅れ指標値を算出し、上記動的指標値から遅れ指標値を減算して、この差分に基づき燃料追加質量を設定する。そして、基本とする吸入燃料質量に、動的指標値と遅れ指標値との差分に基づいて設定した燃料追加質量を加算して１サイクル当りの燃料噴射量を設定する。
【００１４】
【実施例】
以下、図面に基づき本発明の実施例を説明する。図２５にはエンジンの全体概略図が示されており、図中の符号１はエンジン（図においては水平対向４気筒型エンジンを示す）を示し、シリンダヘッド２の吸気ポート２ａにインテークマニホルド３が連通され、このインテークマニホルド３の上流にエアーチャンバ４を介してスロットル通路５が連通されている。このスロットル通路５の上流側には、吸気管６を介してエアクリーナ７が取付けられ、このエアクリーナ７が吸入空気の取り入れ口であるエアインテークチャンバ８に連通されている。また、上記排気ポート２ｂにエキゾーストマニホルド９を介して排気管１０が連通され、この排気管１０に触媒コンバータ１１が介装されてマフラ１２に連通されている。
【００１５】
一方、上記スロットル通路５にスロットル弁５ａが設けられ、このスロットル通路５の直上流の上記吸気管６にインタークーラ１３が介装され、さらに、上記吸気管６の上記エアクリーナ７の下流側にレゾネータチャンバ１４が介装されている。
【００１６】
また、上記レゾネータチャンバ１４と上記インテークマニホルド３とを連通して上記スロットル弁５ａの上流側と下流側とをバイパスするバイパス通路１５に、アイドル空気量を調整するアイドル制御（ＩＳＣ）弁１６が介装されている。さらに、上記ＩＳＣ弁１６の直下流側に、吸気圧が負圧のとき開弁し、またターボチャージャ１８によって過給されて吸気圧が正圧になったとき閉弁するチェックバルブ１７が介装されている。
【００１７】
上記ターボチャージャ１８は、上記吸気管６の上記レゾネータチャンバ１４の下流側にコンプレッサが介装され、タービンが上記排気管１０に介装されている。さらに、上記ターボチャージャ１８のタービンハウジング流入口には、ウエストゲート弁１９が介装され、このウエストゲート弁１９には、ウエストゲート弁作動用アクチュエータ２０が連設されている。
【００１８】
上記ウエストゲート弁作動用アクチュエータ２０は、ダイヤフラムにより２室に仕切られ、一方がウエストゲート弁制御用デューティソレノイド弁２１に連通される圧力室を形成し、他方が上記ウエストゲート弁１９を閉方向に付勢するスプリングを収納したスプリング室を形成している。
【００１９】
上記ウエストゲート弁制御用デューティソレノイド弁２１は、上記レゾネータチャンバ１４と上記吸気管６の上記ターボチャージャ１８のコンプレッサ下流とを連通する通路に介装されており、後述する電子制御装置５０（ＥＣＵ；図２８参照）から出力される制御信号のデューティ比に応じて、上記レゾネータチャンバ１４側の圧力と上記コンプレッサ下流側の圧力とを調圧し、上記ウエストゲート弁作動用アクチュエータ２０の圧力室に供給する。
【００２０】
すなわち、上記ＥＣＵ５０によって上記ウエストゲート弁制御用デューティソレノイド弁２１を制御し、上記ウエストゲート弁作動用アクチュエータ２０を作動させて上記ウエストゲート弁１９による排気ガスリリーフを調整することにより、上記ターボチャージャ１８による過給圧を制御するようになっている。
【００２１】
また、上記インテークマニホルド３に、吸気管圧力Ｐを絶対圧で検出する吸気管圧力センサ２２が通路２３を介して連通され、さらに、上記インテークマニホルド３の各気筒の各吸気ポート２ａの直上流側にインジェクタ２５が臨まされている。また、上記シリンダヘッド２の各気筒毎に、その先端を燃焼室に露呈する点火プラグ２６ａが取付けられ、この点火プラグ２６ａに各気筒毎に配設された点火コイル２６ｂを介してイグナイタ２７が接続されている。
【００２２】
上記インジェクタ２５には、燃料タンク２８内に設けたインタンク式の燃料ポンプ２９から燃料フィルタ３０を経て燃料が圧送され、プレッシャレギュレータ３１にてインジェクタ２５への燃料圧力が調圧される。
【００２３】
また、上記吸気管６の上記エアクリーナ７の直下流に、ホットワイヤ或はホットフィルム等を用いた熱式の吸入空気量センサ３２が介装され、上記スロットル弁５ａに、スロットル開度センサ３３ａとアイドルスイッチ３３ｂとを内蔵したスロットルセンサ３３が連設されている。さらに、上記エアーチャンバ４に吸気温センサ４６が臨まされている。
【００２４】
また、上記エンジン１のシリンダブロック１ａにノックセンサ３４が取付けられるとともに、このシリンダブロック１ａの左右両バンクを連通する冷却水通路３５に水温センサ３６が臨まされ、上記排気管１０の上記エキゾーストマニホルド９の集合部にＯ２ センサ３７が臨まされている。
【００２５】
また、上記シリンダブロック１ａに支承されたクランクシャフト１ｂにクランクロータ３８が軸着され、このクランクロータ３８の外周に、電磁ピックアップなどからなるクランク角センサ３９が対設されている。さらに、上記エンジン１のカムシャフト１ｃに連設するカムロータ４０に、電磁ピックアップなどからなる気筒判別用のカム角センサ４１が対設されている。尚、上記クランク角センサ３９及び上記カム角センサ４１は、電磁ピックアップなどの磁気センサに限らず、光センサなどでも良い。
【００２６】
上記クランクロータ３８は、図２６に示すように、その外周に突起３８ａ，３８ｂ，３８ｃが形成され、これらの各突起３８ａ，３８ｂ，３８ｃが、各気筒（＃１，＃２と＃３，＃４）の圧縮上死点前（ＢＴＤＣ）θ１，θ２，θ３ の位置に形成されており、本実施例においては、θ１ ＝９７゜ＣＡ、θ２ ＝６５゜ＣＡ、θ３ ＝１０゜ＣＡである。
【００２７】
上記クランクロータ３８の各突起は、上記クランク角センサ３９によって検出され、ＢＴＤＣ９７゜，６５゜，１０゜のクランクパルスがエンジン１／２回転毎（１８０゜ＣＡ毎）に出力される。そして、各信号の入力間隔時間がタイマによって計時され、エンジン回転数Ｎｅが算出される。
【００２８】
また、図２７に示すように、上記カムロータ４０の外周には、気筒判別用の突起４０ａ，４０ｂ，４０ｃが形成され、突起４０ａが＃３，＃４気筒の圧縮上死点後（ＡＴＤＣ）θ４ の位置に形成され、突起４０ｂが３個の突起で構成されて最初の突起が＃１気筒のＡＴＤＣθ５ の位置に形成されている。さらに、突起４０ｃが２個の突起で形成され、最初の突起が＃２気筒のＡＴＤＣθ６ の位置に形成されている。本実施例においては、θ４ ＝２０゜ＣＡ、θ５ ＝５゜ＣＡ、θ６ ＝２０゜ＣＡである。
【００２９】
そして、上記カムロータ４０の各突起が上記カム角センサ４１によって検出され、各気筒の燃焼行程順を＃１→＃３→＃２→＃４とした場合、この燃焼行程順と、上記カム角センサ４１からのカムパルスをカウンタによって計数した値とのパターンに基づいて、気筒判別がなされる。
【００３０】
一方、図２８において、符号５０は、エンジン系を制御する電子制御装置（ＥＣＵ）５０であり、このＥＣＵ５０は、燃料噴射制御、点火時期制御などを行なうメインコンピュータ５１と、ノック検出処理を行なう専用のサブコンピュータ５２との２つのコンピュータを中心として構成され、各部に所定の安定化電源を供給する定電圧回路５３や各種の周辺回路が組込まれている。
【００３１】
上記定電圧回路５３は、ＥＣＵリレー５４のリレー接点を介してバッテリ５５に接続されており、このバッテリ５５に、上記ＥＣＵリレー５４のリレーコイルがイグニッションスイッチ５６を介して接続され、上記イグニッションスイッチ５６がＯＮされて上記ＥＣＵリレー５４のリレー接点が閉となったときには制御用電源を上記各コンピュータ５１，５２に供給すると共に、上記定電圧回路５３は上記バッテリ５５に直接接続されており、イグニッションスイッチ５６のＯＮ，ＯＦＦに拘らず、バックアップＲＡＭ６１にバックアップ用の電源を供給する。
【００３２】
尚、上記バッテリ５５には、燃料ポンプリレー５７のリレー接点を介して燃料ポンプ２９が接続されている。
【００３３】
上記メインコンピュータ５１は、ＣＰＵ５８、ＲＯＭ５９、ＲＡＭ６０、バックアップＲＡＭ６１、カウンタ・タイマ群６２、シリアル通信インターフェースであるＳＣＩ６３、及び、Ｉ／Ｏインターフェース６４がバスライン６５を介して接続されたマイクロコンピュータであり、上記バックアップＲＡＭ６１には、上記イグニッションスイッチ５６のＯＮ／ＯＦＦに拘らず、バッテリ５５に直接接続する上記定電圧回路５３からバックアップ電源が常時供給されてデータが保持される。
【００３４】
尚、上記カウンタ・タイマ群６２は、フリーランカウンタ、カム角センサ信号の入力計数用カウンタなどの各種カウンタ、燃料噴射タイマ、点火タイマ、定期割込みを発生させるための定期割込みタイマ、クランク角センサ信号の入力間隔計時用タイマ、及び、システム異常監視用のウオッチドッグタイマなどの各種タイマを便宜上総称するものであり、上記メインコンピュータ５１においては、その他、各種のソフトウエアカウンタ・タイマが用いられる。
【００３５】
また、上記サブコンピュータ５２も、上記メインコンピュータ５１と同様、ＣＰＵ７１、ＲＯＭ７２、ＲＡＭ７３、カウンタ・タイマ群７４、ＳＣＩ７５、及び、Ｉ／Ｏインターフェース７６がバスライン７７を介して接続されたマイクロコンピュータであり、上記メインコンピュータ５１とサブコンピュータ５２とは、上記ＳＣＩ６３，７５を介してシリアル通信ラインにより互いに接続されている。
【００３６】
上記メインコンピュータ５１のＩ／Ｏインターフェース６４には、入力ポートに、吸入空気量センサ３２、スロットル開度センサ３３ａ、水温センサ３６、Ｏ２センサ３７、吸気管圧力センサ２２、大気圧センサ４４、車速センサ４２、吸気温センサ４６、及び、バッテリ５５が、Ａ／Ｄ変換器６６を介して接続されるとともに、アイドルスイッチ３３ｂ、スタータスイッチ４３、及びクランク角センサ３９、カム角センサ４１などが接続され、さらに、図示しない各種センサ及びスイッチ類が接続されている。
【００３７】
また、上記Ｉ／Ｏインターフェース６４の出力ポートには、イグナイタ２７が接続されるとともに、ＩＳＣ弁１６、インジェクタ２５、燃料ポンプリレー５７のリレーコイル、ウエストゲート弁制御用デューティソレノイド弁２１が駆動回路６７を介して接続されており、さらに、図示しない各種のアクチュエータ類が接続されている。
【００３８】
一方、上記サブコンピュータ５２のＩ／Ｏインターフェース７６は、入力ポートに、クランク角センサ３９、カム角センサ４１が接続されるとともに、Ａ／Ｄ変換器７８、周波数フィルタ７９、アンプ８０を介してノックセンサ３４が接続されており、上記ノックセンサ３４からのノック検出信号が上記アンプ８０で所定のレベルに増幅された後に上記周波数フィルタ７９により必要な周波数成分が抽出され、上記Ａ／Ｄ変換器７８にてデジタル信号に変換されて入力されるようになっている。
【００３９】
上記メインコンピュータ５１では、各センサ類からの検出信号を処理し、燃料噴射量制御、点火時期制御、アイドル制御などを行い、一方、上記サブコンピュータ５２では、エンジン回転数とエンジン負荷とに基づいてノックセンサ３４からの信号のサンプル区間を設定し、このサンプル区間でノックセンサ３４からの信号を高速にＡ／Ｄ変換して振動波形を忠実にデジタルデータに変換し、このデータに基づきノック発生の有無を判定する。
【００４０】
上記サブコンピュータ５２のＩ／Ｏインターフェース７６の出力ポートは、上記メインコンピュータ５１のＩ／Ｏインターフェース６４の入力ポートに接続されており、上記サブコンピュータ５２でのノック判定結果がＩ／Ｏインターフェース７６に出力される。そして、上記メインコンピュータ５１では、上記サブコンピュータ５２からノック発生有りの判定結果が出力されると、ＳＣＩ６３を介してシリアル通信ラインよりサブコンピュータ５２からノックデータを読込み、このノックデータに基づいて直ちに該当気筒の点火時期を遅らせ、ノックを回避する。
【００４１】
このようなエンジン制御において、上記メインコンピュータ５１では、センサ・スイッチ類からの信号入力処理、燃料噴射制御、点火時期制御、アイドル制御に係わる各種ジョブが、一つのオペレーティングシステム（ＯＳ）の管理下で効率的に実行される。このＯＳは、車輌制御のための各種マネジメント機能、及び、このマネジメント機能に密着した内部ストラテジーを有し、各種ジョブを体系的に結合し、等時間間隔処理により各種ジョブを効率的に実行するようになっている。
【００４２】
以下、上記メインコンピュータ５１による燃料噴射制御について、図１〜図９のルーチンに従って説明する。尚、サブコンピュータ５２はノック検出処理専用のコンピュータであるため、その動作説明を省略する。
【００４３】
本実施例では、ＯＳ側で処理したＡ／Ｄ変換結果、クランク位置情報、エンジン回転数等に基づき、ユーザ側ジョブで、燃料噴射量、点火時期等を設定すると共に、これらの指示値がＯＳによって噴射タイマ、点火タイマにセットされる。そして、１０ｍｓ毎ジョブ，５０ｍｓ毎ジョブ等の定期割込み要求を出力し、またクランク角センサ３９からの信号入力毎（ＢＴＤＣ９７゜，６５゜，１０゜ＣＡ毎のエンジン１回転に６回）に起動されるクランクパルス入力毎ジョブの割込みを許可する。
【００４４】
１０ｍｓ毎ジョブでは、図１〜図３に示す燃料噴射有効パルス幅及び燃料噴射無駄パルス幅設定ルーチンと、図５に示す吸入空気量設定ルーチンが実行される。また、５０ｍｓ毎ジョブでは、図６に示す係数設定ルーチン、及び図７〜図９に示す要求当量比設定ルーチンが実行される。
【００４５】
尚、以下の説明では、図２４に示すように、吸気系のスロットル弁５ａの下流から吸気弁上流までを吸気チャンバ６Ａと総称する。従って、吸気チャンバ６Ａは、スロットル通路５、エアーチャンバ４、インテークマニホルド３、吸気ポート２ａの総称となる。
【００４６】
ここで、以下の説明で使用する用語について簡単に説明する。Ｑは質量流量［ｍｇ／ｓｅｃ］、Ｇは１サイクル当りの質量［ｍｇ／ｃｙｃｌｅ］を意味する。ただし、後述する筒内吸入空気質量Ｇａｉｒ、筒内吸入燃料質量Ｇｆｕｅｌ等は、吸気行程時にのみ吸入されるので、実質的には、１吸気行程当りの質量と同じ意味になる。また、図２４の左から、Ｑａは吸入空気量センサ３２を通過する空気の質量流量（センサ通過空気質量流量）、Ｍは吸気チャンバ６Ａの空気質量、Ｐは吸気チャンバ６Ａの吸気管内圧力、Ｇｉｎｊは燃料噴射質量、Ｑｃは筒内吸入空気質量流量、Ｇａｉｒは筒内吸入空気質量、Ｇｆｕｅｌは筒内吸入燃料質量を示す。
【００４７】
まず、図１〜図３に示す燃料噴射有効パルス幅及び燃料噴射無駄パルス幅設定ルーチンを説明する前に、このルーチンで取入れる各パラメータを設定するルーチンについて説明する。
【００４８】
図５の１０ｍｓ毎ジョブで実行されるルーチンでは、Ｄジェトロニック方式により１吸気行程当りの筒内吸入空気質量Ｇａｉｒを算出する。まず、ステップＳ４１で、吸気管圧力センサ２２の検出値に基づいて算出したスロットル弁５ａの下流側の吸気管圧力Ｐと、吸気温センサ４６の検出値に基づいて算出した吸気温度Ｔとを読込み、ステップＳ４２で、現在の吸気温度Ｔが５０℃以下、すなわち通常の吸気温度領域にあるか判断し、Ｔ＜５０℃のときには、ステップＳ４３へ進み、現在の吸気温度Ｔで、吸気温度設定値Ｔ’を設定してステップＳ４５へ進む。また、Ｔ≧５０℃のときには、ステップＳ４４で、上記吸気温度設定値Ｔ’を固定値である５０℃に設定した後、ステップＳ４５へ進む。
【００４９】
ところで、筒内吸入空気質量Ｇａｉｒは、上記吸気管圧力Ｐが一定であっても吸気温度Ｔに影響されて変動するため、空気密度σから筒内吸入空気質量Ｇａｉｒを算出する必要がある。筒内の空気密度σが吸気行程時における吸気チャンバ６Ａ内の空気密度と同じであると考えれば、上記筒内の空気密度αは、空気の状態方程式により以下の式から求めることができる。
【００５０】
σ＝Ｐ×（１−Ｒｅｇｒ）／（Ｒ×Ｔ） …（１）
Ｒｅｇｒ；ＥＧＲ率
Ｒ；ガス定数
図１１に筒内吸入空気質量Ｇａｉｒと吸気温度Ｔとの関係を示す。実線は実測により得た吸入空気質量（実測吸気量）の特性、二点鎖線は上記（１）式により算出した空気密度σに従い求めた吸入空気質量（密度比例吸気量）、一点鎖線は吸気温度Ｔをある温度に固定して算出した空気密度σに従い求めた吸入空気質量（圧力比例吸気量）を示す。図に示すように、吸気温度Ｔが５０℃以下の通常の吸気温度領域では、密度比例吸気量が実測吸気量にほぼ沿う特性を示す。しかし、ターボチャージャ付エンジンでの実測結果によると、吸気温度Ｔが５０℃を越えた領域では、上記筒内吸入空気質量Ｇａｉｒは空気密度σに比例して吸込まれておらず、密度比例吸気量は、実測吸気量に比し過小になり、結果として、空燃比がリーン化してしまう。ところで、本来なら、吸気温度Ｔが５０℃以上の領域では、温度毎の補正係数を用いて補正すべきであるが、制御が複雑化するのと、圧力比例吸気量がほぼ実測吸気量に沿う特性を示しているので、本ルーチンでは吸気温度設定値Ｔ’を５０℃に固定して、筒内吸入空気質量Ｇａｉｒを求めるようにした。尚、上記５０℃は実験により求めた臨界温度であり、採用するエンジン形式により相違するため、個別に設定する必要がある。また、本実施例のようにＥＧＲ装置が搭載されていないエンジンでは、上記ＥＧＲ率はゼロである。
【００５１】
そして、ステップＳ４５へ進むと、上記（１）式の吸気温度Ｔを上記吸気温度設定値Ｔ’に代えて空気密度σを算出し、ステップＳ４６で、上記空気密度σと行程容積Ｄとに基づき次式から１吸気行程当りの理論吸気質量Ｇｔｈを算出する。
【００５２】
Ｇｔｈ←Ｄ×σ …（２）
従って、この理論吸気質量Ｇｔｈは、吸気温度設定値Ｔ’を考慮しなければ、上記吸気管圧力Ｐに比例する。
【００５３】
次いで、ステップＳ４７で、上記理論吸気質量Ｇｔｈを基本として筒内吸入空気質量Ｇａｉｒを次式から算出して、ルーチンを抜ける。
【００５４】
Ｇａｉｒ←（Ｇｔｈ−ηｂ）×ηｖ×Ｋｔｒｍ …（３）
ηｂ：吸気損失質量
ηｖ：体積効率
Ｋｔｒｍ：吸入空気量誤差補正係数
上記理論吸気質量Ｇｔｈと筒内吸入空気質量Ｇａｉｒとの関係は、図１２に示す通りであり、理論吸気質量Ｇｔｈがゼロ、すなわち吸気管圧力Ｐが完全な真空になる前に、上記筒内吸入空気質量Ｇａｉｒがゼロになるような、横軸接点ηｂ、傾きηｖの一次式で表される。但し、実際には、筒内吸入空気質量Ｇａｉｒがゼロになるような運転条件は存在しない。また、１吸気行程当りの吸気損失質量ηｂと体積効率ηｖとは、後述するように、エンジン回転数Ｎｅに基づき一次元マップを補間計算付で参照して設定する。さらに、上記吸入空気量誤差補正係数Ｋｔｒｍは、上記理論吸気質量Ｇｔｈと筒内吸入空気質量Ｇａｉｒとの関係が完全には一次式に乗らないこともあるため、運転条件による空気量計測の定常的な誤差分を補正するためのもので、例えば、スロットル開度αとエンジン回転数Ｎｅとに基づき、二次元マップを補間計算付で参照して設定する。図１３に、この吸入空気量誤差補正係数Ｋｔｒｍを設定する二次元マップの特性を例示する。この図に示すように、この吸入空気量誤差補正係数Ｋｔｒｍの値は基本的には、１．０である。
【００５５】
次に、図６に示す５０ｍｓ毎ジョブで実行される係数設定ルーチンについて説明する。このルーチンでは、上記吸気損失質量ηｂ、及び上記体積効率ηｖを設定する。
【００５６】
まず、ステップＳ５１でエンジン回転数Ｎｅを読込み、ステップＳ５２、Ｓ５３で、このエンジン回転数Ｎｅに基づき一次元マップを補間計算付で参照して、体積効率ηｖと吸気損失質量ηｂとを各々設定してルーチンを抜ける。
【００５７】
上記体積効率ηｖ、上記吸気損失質量ηｂの値は、理論的には、１．０であるが、エンジン回転数毎にカムの同調などの影響によって変化する値であり、また上記体積効率ηｖには、高地補正、排気圧補正等の因子を加味させることができる。
【００５８】
このように、本ルーチンでは、吸気管圧力Ｐを基本に筒内吸入空気質量Ｇａｉｒを簡単な計算式に基づいて算出し、また、上記体積効率ηｖ、上記吸気損失質量ηｂは、一次元マップから設定することができるため、コンピュータに掛かる負荷が軽くなり、しかも、上記吸気管圧力Ｐさえ正確に計測されれば、始動時を含む全ての領域で適用することができる。
【００５９】
ところで、上記筒内吸入空気質量Ｇａｉｒは、Ｌジェトロニック方式により算出することもできる。このＬジェトロニック方式を利用して筒内吸入空気質量Ｇａｉｒを算出する手順を、図１０の４ｍｓ毎ジョブで実行されるルーチンに従って説明する。尚、上記吸入空気量センサ３２の検出時刻は、スロットル弁５ａまでの吸気管路長、及び流速によっても異なるが、スロットル弁５ａでの検出時刻とほぼ同じと考えて良い。
【００６０】
まず、ステップＳ１０１で、吸入空気量センサ３２からの出力電圧ＶＡＦＭを読込み、次のステップＳ１０２，Ｓ１０３で上記出力電圧ＶＡＦＭの過渡補正をする。すなわち、図１４に示すように、実験によれば、過渡時の吸入空気量の急激な変化に対して、吸入空気量センサ３２の出力電圧ＶＡＦＭの変化は、理想の変化の２／３の電圧までは充分に速い応答を示すが、残りの１／３は、約２００ｍｓの時定数を持つ遅れで応答することが明らかにされている。このときのセンサ通過空気質量流量Ｑａと出力電圧ＶＡＦＭとの関係を、時定数τ（ｓ）の伝達関数で表せば、
ＶＡＦＭ＝（２／３）×Ｑａ＋｛１／（３τ（ｓ）＋３）｝×Ｑａ …（４）
となり、応答遅れは、この（４）式の逆関数で補正すれば良く、この補正式は、
Ｑａ＝（１／２）×［３−１／｛（２／３）τ（ｓ）＋１｝］×ＶＡＦＭ …（５）
となる。
【００６１】
上記時定数τ（ｓ）は２００ｍｓであるため、（２／３）τ（ｓ）＋１≒１３０とすれば、上記（５）式は、
と単純化される。すなわち、吸入空気量センサ３２の出力電圧ＶＡＦＭの３倍から、この出力電圧ＶＡＦＭを１３０ｍｓの時定数で一次遅れさせた出力電圧ＶＡＦＭを減算し、その値を１／２にすれば、この吸入空気量センサ３２の過渡時の応答遅れを、簡単に補正することができる。
【００６２】
このような考えに基づき、上記ステップＳ１０２では、前回求めた一次遅れ出力電圧ＶＡＦＭａｖと今回の出力電圧ＶＡＦＭとから、次式に示す所定時定数（本ルーチンでは、１３０ｍｓの近似値として３２×４ｍｓ）の加重平均に基づいて、今回の一次遅れ出力電圧ＶＡＦＭａｖを算出する。
【００６３】
ＶＡＦＭａｖ←（３１・ＶＡＦＭａｖ＋ＶＡＦＭ）／３２ …（７）
次いで、ステップＳ１０３で、今回の出力電圧ＶＡＦＭ及び一次遅れ出力電圧ＶＡＦＭａｖから、過渡補正後の出力電圧（過渡補正出力電圧）Ｖ’ＡＦＭを次式から算出する。
【００６４】
Ｖ’ＡＦＭ←（３・ＶＡＦＭ−ＶＡＦＭａｖ）／２ …（８）
このように、上記吸入空気量センサ３２の過渡応答遅れを簡単な一次式で導き出すことができるため、コンピュータの負荷を軽減させることができる。
【００６５】
そして、ステップＳ１０４で、過渡補正された上記過渡補正出力電圧Ｖ’ＡＦＭを、等間隔３２格子のワード補間テーブルを用いてセンサ通過空気質量流量Ｑａに変換し、ステップＳ１０５で、このセンサ通過空気質量流量Ｑａに吸入空気量誤差補正係数Ｋｔｒｍを乗算して、今回のセンサ通過空気質量流量Ｑ’ａを算出する。
【００６６】
この吸入空気量誤差補正係数Ｋｔｒｍは、運転条件による空気量計測の定常的な誤差分を補正する係数で、前述のＤジェトロニック方式と同様に、例えば、スロットル開度αとエンジン回転数Ｎｅとに基づき、二次元マップを補間計算付で参照して設定する（図１３参照）。
【００６７】
次いで、ステップＳ１０６で、エンジンが始動後かをエンジン回転数Ｎｅに基づいて判断し、Ｎｅ＜３００ｒｐｍのときは、始動時と判断してステップＳ１０７へ進み、Ｎｅ≧３００ｒｐｍのときは、始動後と判断してステップＳ１０８へ進む。
【００６８】
ステップＳ１０７へ進むと、始動時の筒内吸入空気質量Ｇａｉｒを初期値設定してルーチンを抜ける。
【００６９】
始動時に、筒内吸入空気質量Ｇａｉｒを初期値設定するのは、以下の理由による。
【００７０】
１）筒内吸入空気質量Ｇａｉｒは、定常的には、空気流量／エンジン回転数で算出されるが、エンジン回転数及び空気流量が共に小さい場合を考えると、計算結果は不定になり易い。
【００７１】
２）後述のステップＳ１０８で説明するように、筒内吸入空気質量Ｇａｉｒは、始動後は履歴に依存して逐次計算されるため、始動後最初のルーチン実行時に初期値が設定されている必要がある。
【００７２】
３）始動時に、ある始動時噴射幅を設定するのではなく、筒内吸入空気質量Ｇａｉｒを初期値設定することで、始動時、始動後に拘らず全ての運転領域で筒内吸入空気質量Ｇａｉｒを利用した演算が可能になる。
【００７３】
そこで、始動時には、筒内吸入空気質量Ｇａｉｒを推定した値で、初期値設定する。この初期値は、例えば、吸気チャンバ６Ａ内の吸気管圧力Ｐを大気圧とみなし、また吸気温度を冷却水温度により代用して次式から算出する。
【００７４】
初期値＝固定値×大気圧／冷却水温（゜Ｋ） …（９）
一方、始動後と判断されてステップＳ１０８へ進むと、上記ステップＳ１０５で算出したセンサ通過空気質量流量Ｑ’ａに基づき、筒内吸入空気質量Ｇａｉｒを次式に示す吸気チャンバモデルから算出して、ルーチンを抜ける。
Ｇａｉｒ←Ｑ’ａ×Ｍｔｃｈ＋（ＭＴＣＳ−Ｍｔｃｈ）×Ｇａｉｒ／ＭＴＣＳ …（１０）
ＭＴＣＳ：半回転時間
この（１０）式は、以下の考えに基づく。すなわち、吸気チャンバモデルの計算は、次の２式から算出することができる。
【００７５】
１）Ｇａｉｒ＝σ×Ｄ＝Ｍ×Ｄ／Ｖ
σ：吸気チャンバ６Ａ内の空気密度
Ｄ：行程容積
Ｖ：吸気チャンバ６Ａ内の容積（チャンバ容積）
２）Ｍ（ｎｅｗ）＝Ｇａｉｒ（ｏｌｄ）＋Ｑａ×Δｔ−Ｑｃ×Δｔ
（ｎｅｗ）：現在の値
（ｏｌｄ）：Δｔ前の値
Δｔ：演算周期
上記２式からＭとＱｃとを消去すると、次式に示すように、上記センサ通過空気質量流量Ｑａから筒内吸入空気質量Ｇａｉｒを算出することができる。
Ｇａｉｒ（ｎｅｗ）＝Ｇａｉｒ（ｏｌｄ）＋Ｑａ×Δｔ×Ｄ／Ｖ−Ｑｃ×Δｔ×Ｄ／Ｖ …（１０−１）
上記筒内吸入空気質量流量Ｑｃは、単位時間当りの吸気行程数が４サイクルエンジンの場合、気筒数×Ｎｅ／２であるため、
Ｑｃ＝（気筒数／２）×Ｎｅ×Ｇａｉｒ
となり、上記（１０−１）式は、
…（１０−２）
となる。ここで、Δｔ×Ｄ／Ｖ＝固定値＝Ｍｔｃｈとおき、また、採用するエンジンが４気筒であれば、上記（１０−２）式は、
Ｇａｉｒ（ｎｅｗ）＝Ｑａ×Ｍｔｃｈ＋（１−２Ｎｅ×Ｍｔｃｈ）×Ｇａｉｒ（ｏｌｄ）…（１０−３）
となる。尚、本ルーチンでは、計算の便宜上、エンジン回転数Ｎｅ［ｒｐｓ］ではなく半回転数時間（Ｎｅ／２）を用いており、この半回転数時間をＭＴＣＳとおけば、上記（１０−３）式は、
Ｇａｉｒ（ｎｅｗ）＝Ｑａ×Ｍｔｃｈ＋Ｇａｉｒ（ｏｌｄ）×（ＭＴＣＳ−Ｍｔｃｈ）／ＭＴＣＳ …（１０’）
となり、上記ステップＳ１０８で示した（１０）式が導き出される。
【００７６】
そして、この（１０）式を演算周期Δｔ毎に行えば、筒内吸入空気質量Ｇａｉｒを履歴に依存して逐次的に算出することができる。
【００７７】
このように、本実施例では、Ｄジェトロニック方式とＬジェトロニック方式との違いは、単に吸気管圧力センサ２２の検出値に基づいて筒内吸入空気質量Ｇａｉｒを算出するか、吸入空気量センサ３２の検出値に基づいて筒内吸入空気質量Ｇａｉｒを算出するか、すなわち、この筒内吸入空気質量Ｇａｉｒの算出方式の違いに過ぎず、以下のルーチンでは、筒内吸入空気質量Ｇａｉｒが算出されていれば、Ｄジェトロニック方式とＬジェトロニック方式とのいずれを採用するエンジンであっても共通となる。
【００７８】
次に、運転条件による筒内混合気の要求当量比設定ルーチンについて、図７〜図９のルーチンに従って説明する。このルーチンは５０ｍｓ毎のジョブで実行され、定常的な要求増量である可燃限界当量比φｔｗと、最大出力当量比φｆｕｌと、排気ガス温度限界当量比φｔｅｘと、触媒浄化率最良当量比φｇａｓと、燃費率最良当量比φｅｃｏとを設定する。これら各要求当量比φｔｗ，φｆｕｌ，φｔｅｘ，φｇａｓ，φｅｃｏは定常的な要求増量であるため、過渡的な燃料遅れ等によって生じる空燃比変動をフィードフォワード的に補正する係数（エアコン増量補正係数、Ｄレンジ増量補正係数等）、或はセンサやインジェクタの誤差を補正する係数（混合比補正係数等）は含まない。また、このルーチンで設定する各要求当量比は、後述する図１のルーチン実行時に比較され、その最大値が目標当量比（理論燃空比Ｆ／Ａに対する増量係数）ＣＯＥＦとして採用される。尚、この当量比φは、空気過剰率（（Ａ／Ｆ）実際／（Ａ／Ｆ）理論）λに対して、φ＝１／λの関係にある。
【００７９】
まず、ステップＳ６１で冷却水温Ｔｗとエンジン回転数Ｎｅとを読込む。そして、ステップＳ６２〜Ｓ６５或はＳ６６で可燃限界当量比φｔｗを設定する。ステップＳ６２では、上記冷却水温Ｔｗとエンジン回転数Ｎｅとに基づき二次元マップを補間計算付で参照して、可燃限界当量比φｔｗを設定する際の基準となる基準可燃限界当量比Ｍｔｗを設定する。この基準可燃限界当量比Ｍｔｗは、冷却水温Ｔｗとの関係によるリーン側の可燃限界を示すもので、図１５に、この基準可燃限界当量比Ｍｔｗを設定する二次元マップの特性を例示する。この図１５に示すように、冷却水温Ｔｗ側の軸では、従来の水温増量に近い設定がなされており、また、ある程度の高水温では、上記基準可燃限界当量比Ｍｔｗはリーン当量比に設定される。また、エンジン回転数Ｎｅ側の軸では、クランキング回転数を含む全ての回転数領域を検索できるように設定されている。筒内吸入空気質量Ｇａｉｒは始動時においても設定されることから（図５或は図１０のフローチャート参照）、このときの可燃限界当量比φｔｗが始動時の噴射パルス幅を決定する要因になる。また、始動時において、上記基準可燃限界当量比Ｍｔｗは大きな値に設定されているが、これが始動後にエンジン回転数の増加に従って小さな値に変化して行く際に、後述するように、ある一次遅れ時定数で上記可燃限界当量比φｔｗの変化を遅らせることで始動後の増量が確保される。
【００８０】
次いで、ステップＳ６３へ進むと、エンジン回転数等からエンスト状態にあるかを判断し、エンスト状態にあると判断した場合には、ステップＳ６６へ分岐し、上記基準可燃限界当量比Ｍｔｗで今回の可燃限界当量比φｔｗを設定してステップＳ６７へジャンプする。一方、エンスト状態にないと判断した場合には、ステップＳ６４へ進み、上記基準可燃限界当量比Ｍｔｗと前回の可燃限界当量比φｔｗとを比較し、Ｍｔｗ≧φｔｗのときは、要求当量比が増量されているので、ステップＳ６６へ分岐し、この基準可燃限界当量比Ｍｔｗで今回の可燃限界当量比φｔｗを設定して、増量要求に対応させた後、ステップＳ６７へジャンプする。一方、Ｍｔｗ＜φｔｗのときは、要求当量比が減少されているため、ステップＳ６５へ進み、次式に示す所定時定数（本ルーチンでは、１６×５０ｍｓ）の加重平均により一次遅れ処理を行い、今回の可燃限界当量比φｔｗを設定する。その結果、始動後増量が実質的に確保される。
【００８１】
φｔｗ←（１５φｔｗ＋Ｍｔｗ）／１６ …（１１）
次いで、ステップＳ６７〜Ｓ７０或はＳ７１で上記最大出力当量比φｆｕｌが設定される。まず、ステップＳ６７では、スロットル開度αを読込み、ステップＳ６８で、エンジン回転数Ｎｅに基づき、一次元マップを補間計算付で参照してフル増量基準値αｔｈを設定する。このフル増量基準値αｔｈは、エンジンが最大出力を要求しているかを判断する基準値で、図１６に示すように、エンジン回転数Ｎｅとスロットル開度αとの関係から特性を求めてマップ化したものである。
【００８２】
そして、ステップＳ６９で上記スロットル開度αと上記フル増量基準値αｔｈとを比較し、α≦αｔｈのときは、フル増量条件不成立と判断し、ステップＳ７０へ進み、最大出力当量比φｆｕｌを０として、ステップＳ７２へ進む。一方、α＞αｔｈのときは、フル増量条件成立と判断し、ステップＳ７１へ分岐し、上記最大出力当量比φｆｕｌを１．２、すなわち、当該エンジンの最大出力となる当量比（フル増量）に設定する。
【００８３】
尚、目標当量比ＣＯＥＦは、本ルーチンにより設定する各要求当量比φｔｗ，φｆｕｌ，φｔｅｘ，φｇａｓ，φｅｃｏの中から最大値を選択的に設定するものであるため、上記ステップＳ７０で最大出力当量比φｆｕｌが０に設定されると、この最大出力当量比φｆｕｌが目標当量比ＣＯＥＦの選択対象から外れることになる。また、排気ガス温度の上昇を抑える目的の増量は、排気ガス温度限界当量比φｔｅｘとして別に設定されているため、上記最大出力当量比φｆｕｌはエンジン回転数Ｎｅによらず、理論空燃比Ａ／Ｆ（＝１４．６）に対して、おおよそ１５〜２０％増量となるはずである。
【００８４】
次いで、ステップＳ７２〜Ｓ７９或はＳ８０で排気ガス温度限界当量比φｔｅｘを設定する。まず、ステップＳ７２で、負荷の一例である筒内吸入空気質量Ｇａｉｒとエンジン回転数Ｎｅと基づき二次元マップを補間計算付で参照して基本排気ガス温度限界当量比Ｍｔｅｘを設定する。この基本排気ガス温度限界当量比Ｍｔｅｘは、排気ガス温度の上昇を設計限界以下に抑え、エンジン及び排気系を保護するために燃料冷却を行う増量である。図１７に、この基本排気ガス温度限界当量比Ｍｔｅｘを設定する二次元マップの特性を例示する。図に示すように、高負荷（高Ｇａｉｒ）、高回転で大きな増量を行うように設定されている。尚、エンジン負荷を検出するパラメータとして上記筒内吸入空気質量Ｇａｉｒに代えて吸気管圧力Ｐを採用し、二次元マップは吸気管圧力Ｐとエンジン回転数Ｎｅとの関係による特性に基づいて設定するようにしても良い。
【００８５】
そして、ステップＳ７３へ進むと、上記基本排気ガス温度限界当量比Ｍｔｅｘと、上記最大出力当量比φｆｕｌのフル増量値である１．２とを比較し、Ｍｔｅｘ≧１．２のときはステップＳ７４へ進み、またＭｔｅｘ＜１．２のときはステップＳ７６へ進む。
【００８６】
そして、上記ステップＳ７６へ進むと、今回の平均排気ガス温度限界当量比φｔｅｘＡＶを次式に示す加重平均により、所定時定数（本ルーチンでは、１６×５０ｍｓ）の一次遅れ処理を行い、ステップＳ７８へ進む。
【００８７】
φｔｅｘＡＶ←（１５・φｔｅｘＡＶ＋Ｍｔｅｘ）／１６ …（１２）
また、ステップＳ７４へ進むと、上記基本排気ガス温度限界当量比Ｍｔｅｘと前回までの平均排気ガス温度限界当量比φｔｅｘＡＶとを比較し、Ｍｔｅｘ＜φｔｅｘＡＶのとき、すなわち、１．２≦Ｍｔｅｘ＜φｔｅｘＡＶのときは、ステップＳ７５へ進み、またＭｔｅｘ≧φｔｅｘＡＶのときは、ステップＳ７７へ進む。
【００８８】
上記ステップＳ７５では、今回の平均排気ガス温度限界当量比φｔｅｘＡＶを次式に示す加重平均により所定時定数（本ルーチンでは、３２×５０ｍｓ）の一次遅れ処理を行い、ステップＳ７８へ進む。
【００８９】
φｔｅｘＡＶ←（３１・φｔｅｘＡＶ＋Ｍｔｅｘ）／３２ …（１３）
一方、ステップＳ７７へ進むと、今回の平均排気ガス温度限界当量比φｔｅｘＡＶを次式に示す加重平均により所定時定数（本ルーチンでは１２８×５０ｍｓ）の一次遅れ処理を行って、ステップＳ７８へ進む。
【００９０】
φｔｅｘＡＶ←（１２７・φｔｅｘＡＶ＋Ｍｔｅｘ）／１２８ …（１４）
従って、Ｍｔｅｘ＜１．２のいわゆる要求当量比が薄い場合には、比較的早い時定数で平均排気ガス温度限界当量比φｔｅｘＡＶが設定され、また、Ｍｔｅｘ≧１．２、且つＭｔｅｘ≧φｔｅｘＡＶの要求当量比が濃い場合には、ゆっくりとした時定数で上記平均排気ガス温度限界当量比φｔｅｘＡＶが設定される。
【００９１】
そして、ステップＳ７８へ進むと、上記ステップＳ７５，Ｓ７６或はＳ７７で算出した平均排気ガス温度限界当量比φｔｅｘＡＶと上記基本排気ガス温度限界当量比Ｍｔｅｘとを比較し、Ｍｔｅｘ＜φｔｅｘＡＶのときは、ステップＳ７９で、上記基本排気ガス温度限界当量比Ｍｔｅｘを今回の排気ガス温度限界当量比φｔｅｘとして設定し、また、Ｍｔｅｘ≧φｔｅｘＡＶのときは、ステップＳ８０で、今回の排気ガス温度限界当量比φｔｅｘを平均排気ガス温度限界当量比φｔｅｘＡＶで設定する。
【００９２】
従って、排気ガス温度限界当量比φｔｅｘは平均排気ガス温度限界当量比φｔｅｘＡＶで上限リミットされ、また、基本排気ガス温度限界当量比Ｍｔｅｘで下限リミットされる。その結果、上記排気ガス温度限界当量比φｔｅｘの上限リミット値は、エンジン回転数Ｎｅと筒内吸入空気質量Ｇａｉｒとに応じて、すなわち、上記ステップＳ７５，Ｓ７６或はＳ７７で設定される平均排気ガス温度限界当量比φｔｅｘＡＶに従って、ゆっくり変化したり、あるいは比較的早く変化したりするようになる。その結果、長時間の全開運転において排気ガス温度が実際に設計限界値に近づくまではφｔｅｘ＜φｆｕｌとなり、従来のように、最大出力当量比φｆｕｌの中に、燃料冷却成分を含ませていた場合に比し、出力の増加が期待できる。尚、この排気ガス温度限界当量比φｔｅｘを排気ガス温度を計測或は推測した値に基づいて設定するようにしても良い。
【００９３】
次いで、上記ステップＳ７９或はＳ８０からステップＳ８１へ進み、このステップＳ８１〜Ｓ８５或はＳ８６で触媒浄化率最良当量比φｇａｓが設定される。この触媒浄化率最良当量比φｇａｓは、触媒コンバータ１１として３元触媒が採用されている場合に対象となり、排気ガスを浄化しなければならない条件が成立したときに、φｇａｓ←１．０とし、それ以外では、φｇａｓ←０とする。
【００９４】
ステップＳ８１〜Ｓ８４では、排気ガスを浄化しなければならない運転条件にあるかを判断する。この判断条件は、スロットル開度αとエンジン回転数Ｎｅとで行い、スロットル開度αが、１０゜≦α＜７０゜（ステップＳ８１，Ｓ８２）で、且つエンジン回転数Ｎｅが、８００ｒｐｍ≦Ｎｅ＜３０００ｒｐｍ（ステップＳ８３，Ｓ８４）のとき、触媒浄化率最良当量比φｇａｓを１．０に設定し（ステップＳ８５）、それ以外のときは、φｇａｓを０に設定する（ステップＳ８６）。尚、上記判断条件は一例であり、排気ガス浄化領域はエンジンの特性等に応じて任意に設定でき、また例えばスロットル開度αに代えて、吸気管圧力Ｐや筒内吸入空気質量Ｇａｉｒを採用しても良い。また、上記触媒コンバータ１１としてリーンＮｏｘ触媒が採用されている場合には、触媒浄化率最良当量比φｇａｓは全運転領域でφｇａｓ←０に設定する。
【００９５】
その後、ステップＳ８７で、負荷の一例としての吸気管圧力Ｐとエンジン回転数Ｎｅとに基づいて二次元マップを補間計算付で参照して、燃費率最良当量比φｅｃｏを設定して、ルーチンを抜ける。
【００９６】
図１８に、この燃費率最良当量比φｅｃｏを設定する二次元マップの特性を例示する。図に示すように、この二次元マップの各領域には、当該エンジン運転条件下での最良の燃費率を得ることのできる当量比が予め実験などから求めて格納されている。
【００９７】
ここで、この燃費率最良当量比φｅｃｏと前述の可燃限界当量比φｔｗとの関係について簡単に説明する。
【００９８】
空燃比を可燃限界当量比φｔｗまでリーン化させると一般には最良の燃費を得ることができなくなる。しかし、暖機途中の低水温時などでは空燃比を燃費率最良当量比φｅｃｏまでリーン化するよりも、上記可燃限界当量比φｔｗにより、リーン化を抑制した方が良好な燃焼を得ることができる場合もあるので、上記可燃限界当量比φｔｗとは別に燃費率最良当量比φｅｃｏが設定されている。尚、上記触媒浄化率最良当量比φｇａｓが全運転領域で、φｇａｓ←１．０の場合には、上記燃費率最良当量比φｅｃｏは隠れてしまうため、その意味を失うが、例えば、排気浄化システムがリーンＮｏｘ触媒を採用している場合には、前述のように、触媒浄化率最良当量比φｇａｓが全運転領域でφｇａｓ←０に設定されるため、この燃費率最良当量比φｅｃｏは触媒浄化率最良当量比としての要素も合わせ持つことになる。また、ＥＧＲを用いて理論空燃比制御により燃費向上を行う際には、Ｄジェトロニック方式では、上記（１），（２）式から明らかなように、理論吸気質量Ｇｔｈには、ＥＧＲ量が除かれた新気量のみであり、またＬジェトロニック方式では、吸入空気量センサ３２の出力電圧ＶＡＦＭに基づいて算出した筒内吸入空気質量Ｇａｉｒは、結果としてＥＧＲ量が除かれた新気量のみとなるため、全運転領域で、φｅｃｏ←１．０としたほうが、空燃比制御性が良くなる。また、この燃費率最良当量比φｅｃｏを、吸気管圧力Ｐに代えて筒内吸入空気質量Ｇａｉｒと、エンジン回転数Ｎｅとの関係により作成した二次元マップを参照して設定するようにしても良い。
【００９９】
そして、上記各ルーチンで求めた設定値を用いて、図１〜図３の１０ｍｓ毎ジョブで実行される燃料噴射有効パルス幅及び燃料噴射無駄パルス幅設定ルーチンのステップＳ１〜Ｓ１０で１サイクル当りの燃料噴射質量を算出し、ステップＳ１１〜Ｓ２４で燃料噴射有効パルス幅Ｔｅと燃料噴射無駄パルス幅Ｔｓとを算出する。
【０１００】
まず、ステップＳ１では、上記筒内混合気の要求当量比設定ルーチンで設定した各要求当量比φｔｗ，φｆｕｌ，φｔｅｘ、φｇａｓ，φｅｃｏを比較し、これらの要求値の最大値を目標当量比ＣＯＥＦとして決定する。
【０１０１】
次いで、ステップＳ２で筒内への１吸気行程当りの筒内吸入燃料質量Ｇｆｕｅｌを次式から算出する。
【０１０２】
Ｇｆｕｅｌ←Ｇａｉｒ×Ｆ／Ａ×ＣＯＥＦ×Ｋｆｂ …（１５）
Ｆ／Ａ：理論燃空比
Ｋｆｂ：空燃比フィードバック補正係数
ここで、理論燃空比Ｆ／Ａは理論空燃比Ａ／Ｆの逆数で、理論空燃比が１４．６であれば、理論燃空比は１／１４．６となる。また、この理論空燃比Ａ／Ｆは、燃料を完全燃焼するに必要な最小の空気量と燃料量との比であり、多種燃料に対応したエンジン等では変数となる。
【０１０３】
一方、上記空燃比フィードバック補正係数Ｋｆｂの初期値は１．０であり、空燃比フィードバック制御、及び空燃比学習制御を行う外部のストラテジーによって書換えられる。また、理論燃空比Ｆ／Ａに目標当量比ＣＯＥＦを掛けることで、運転条件に応じた筒内混合気の目標空燃比が設定される。従って、リニアＡ／Ｆセンサでフィードバック制御する際に設定される上記空燃比フィードバック補正係数Ｋｆｂの目標値は上記理論燃空比Ｆ／Ａを１とした場合の上記目標当量比ＣＯＥＦである。
【０１０４】
ところで、上記筒内吸入燃料質量Ｇｆｕｅｌは、筒内へ１行程で吸込んで欲しい燃料量であり、定常的にはインジェクタ２５からの噴射量と一致するが、過渡時には応答遅れがあるため一致しない。すなわち、過渡的にはインジェクタ２５からの燃料噴射量が急増しても、燃料の一部が吸気ポート内壁に付着してから筒内に吸込まれるので、その分遅れが生じ、筒内に吸込まれる燃料量はゆっくりと増加する。従って、上記筒内吸入燃料質量Ｇｆｕｅｌが急増したときには、インジェクタ２５からの燃料量を、この筒内吸入燃料質量Ｇｆｕｅｌよりも多めにフィードフォワード噴射させて、筒内への吸入燃料量を筒内吸入燃料質量Ｇｆｕｅｌと一致させる必要がある。この追加増量分に相当する過渡時の１吸気行程当りの燃料追加質量Ｇａｃｃを、次のステップＳ３〜Ｓ９で算出する。
【０１０５】
まず、ステップＳ３では、スロットル開度αとエンジン回転数Ｎｅとに基づき二次元マップ（図１９参照）を補間計算付で参照して、静的な筒内吸入空気質量Ｇａｉｒに相当する指標値（静的指標値）Ｍａｃｃを設定する。この静的指標Ｍａｃｃは、負荷変動の少ない定常走行で、スロットル開度αとエンジン回転数Ｎｅを一定に保持すれば、最終的には筒内吸入空気質量Ｇａｉｒが一定になるということを前提にしている。従って、走行中の負荷変動が無視できるとすれば、図４のステップＳ３１に示すように、上記静的指標Ｍａｃｃをスロットル開度αに基づき一次元マップから補間計算付で設定することも可能である。
【０１０６】
また、過渡的には、上記筒内吸入空気質量Ｇａｉｒにも遅れがあり、この遅れに対応した燃料追加質量Ｇａｃｃを求める必要がある。
【０１０７】
そこで、ステップＳ４では、まず、吸気チャンバ６Ａ（図２４参照）内の吸気遅れを想定した所定時定数（運転条件によって異なるが本ルーチンでは４×１０ｍｓに固定している）の一次遅れ処理を、次式に示す加重平均に基づいて行い、動的な筒内吸入空気質量Ｇａｉｒに相当する指標値（動的指標値）Ｓａｃｃを算出する。
【０１０８】
Ｓａｃｃ←（３・Ｓａｃｃ＋Ｍａｃｃ）／４ …（１６）
次いで、ステップＳ５で、上記動的指標値Ｓａｃｃと、前回のルーチン実行時に算出した燃料付着による遅れ指標値Ｔａｃｃとを比較する。尚、この遅れ指標値Ｔａｃｃは上記筒内吸入燃料質量Ｇｆｕｅｌに相当する。
【０１０９】
そして、Ｓａｃｃ＜Ｔａｃｃ、すなわち、前回の筒内吸入燃料質量Ｇｆｕｅｌに相当する値（Ｔａｃｃ）に対して今回の筒内吸入空気質量Ｇａｉｒに相当する値（Ｓａｃｃ）が少ないときには、ステップＳ６へ進み、上記遅れ指標値Ｔａｃｃを上記動的指標値Ｓａｃｃとして、ステップＳ８へ進む。その結果、Ｓａｃｃ＝Ｔａｃｃとなり、後述のステップＳ９では、上記燃料追加質量Ｇａｃｃ＝０となり、燃料が追加増量されない。
【０１１０】
一方、Ｓａｃｃ≧Ｔａｃｃ、すなわち、前回の筒内吸入燃料質量Ｇｆｕｅｌに相当する値（Ｔａｃｃ）に対して今回の筒内吸入空気質量Ｇａｉｒに相当する値（Ｓａｃｃ）が増加しているか、或は変化していないときには、ステップＳ７へ分岐して、この動的指標Ｓａｃｃを、次式に示す加重平均により、燃料付着による時定数（本ルーチンでは、８×１０ｍｓ）の一次遅れ処理を行って、今回の燃料付着による遅れ指標値Ｔａｃｃを算出して、ステップＳ８へ進む。
【０１１１】
Ｔａｃｃ←（７・Ｔａｃｃ＋Ｓａｃｃ）／８ …（１７）
そして、ステップＳ８へ進むと、エンジン回転数Ｎｅに基づき一次元マップを補間計算付で参照して係数Ｒａｃｃを設定する。
【０１１２】
その後、ステップＳ９へ進み、過渡時の１吸気行程当りの燃料追加質量Ｇａｃｃを次式から算出する。
【０１１３】
Ｇａｃｃ←Ｒａｃｃ×（Ｓａｃｃ−Ｔａｃｃ） …（１８）
この（１８）式に示すように、今回の筒内吸入空気質量Ｇａｉｒに相当する値（Ｓａｃｃ）と今回の筒内吸入燃料質量Ｇｆｕｅｌに相当する値（Ｔａｃｃ）との差分が、燃料付着による不足分を補填する増量パターンということになる。
【０１１４】
このように、この燃料追加質量Ｇａｃｃは、スロットル開度αとエンジン回転数Ｎｅのみから、簡易的に筒内吸入空気質量Ｇａｉｒの変化量だけを遅れなく検出し、この値を利用して算出している。
【０１１５】
図２０にスロットル弁５ａを走行中に僅かに開弁させたときの上記各指標値Ｍａｃｃ，Ｓａｃｃ，Ｔａｃｃの特性を示す。静的指標値Ｍａｃｃはスロットル開度αに追従して増加し、一方、動的指標値Ｓａｃｃは上記静的指標値Ｍａｃｃに対して４×１０ｍｓの一次遅れで増加する。さらに、上記遅れ指標値Ｔａｃｃは、この動的指標値Ｓａｃｃに対して８×１０ｍｓの時定数で増加する。この遅れ指標値Ｔａｃｃが筒内吸入燃料質量Ｇｆｕｅｌに対応していると考えれば、付着分の遅れは、上記動的指標値Ｓａｃｃと上記遅れ指標値Ｔａｃｃとの差分（図のハッチングで示す領域）に相当することになり、上記ステップＳ９では、この差分指標値（Ｓａｃｃ−Ｔａｃｃ）に上記係数Ｒａｃｃを掛けて、燃料追加質量Ｇａｃｃを求めている。ところで、図２１に示すように、この差分指標値（Ｓａｃｃ−Ｔａｃｃ）が、従来のスロットル弁５ａの開弁直後に生じるリーンスパイクのパターンに符合することが、実験から明らかになった。従って、この差分指標値（Ｓａｃｃ−Ｔａｃｃ）相応分を加速増量すれば、排気空燃比は過渡時においてもリーン化することなく一定となる。
【０１１６】
上記燃料追加質量Ｇａｃｃは加速増量・加速時追加パルスに近似するが、この燃料追加質量Ｇａｃｃはあくまでも筒内への吸入燃料量を筒内吸入燃料質量Ｇｆｕｅｌにするためのフィードフォワード分であり、加速時に空燃比Ａ／Ｆを一時的にリッチにするためのものではない。従って、もし、加速時に空燃比Ａ／Ｆをリッチにする必要があれば、前記目標当量比ＣＯＥＦの選択枝として、前記図７〜図９に示す要求当量比設定ルーチン中に、加速当量比φａｃｃを第６の増量要因として加えることも、当然考えられる。
【０１１７】
また、先に述べた、従来の三角増量法では、筒内へ吸入される燃料量の一次遅れを三角形で近似させながら増量するので、上記図２１に示すリーンスパイクに対して正確に符合せず、過渡時の空燃比が部分的にリッチになったりリーンになったりしてしまう。
【０１１８】
そして、ステップＳ１０へ進むと、上記筒内吸入燃料質量Ｇｆｕｅｌに上記燃料追加質量Ｇａｃｃを加算して、１サイクル当りの燃料噴射質量Ｇｉｎｊを算出する。ところで、噴射時期が早期に設定されていて、噴射が終了した後でまだ吸気行程になる前のタイミングで上記燃料噴射質量Ｇｉｎｊが急増した場合には、自動的に追加噴射することも可能で、この燃料噴射質量Ｇｉｎｊは、ある吸気行程から１回転半前、すなわち吸気行程終了直後から次の吸気行程開始までの期間にインジェクタ２５から噴射される燃料の総量が算出される。
【０１１９】
次に、ステップＳ１１〜Ｓ２４で、気筒別燃料噴射有効パルス幅Ｔｅ１〜Ｔｅ４、及び気筒別燃料噴射無駄パルス幅Ｔｓ１〜Ｔｓ４が算出される（ここで、１，２，３，４は気筒番号を示す）。
【０１２０】
まず、ステップＳ１１で、最大出力当量比φｆｕｌの値を参照し、ステップＳ１２で排気ガス温度限界当量比φｔｅｘの値を参照する。そして、φｆｕｌ＝０のエンジンが最大出力を要求しておらず、しかも、φｔｅｘ＝０の排気ガス温度限界による増量も要求していない場合には、ステップＳ１３へ進み、空燃比の気筒別のばらつきを補正する気筒別噴射量補正係数Ｋｔｎ１〜Ｋｔｎ４を１．０に設定して、ステップＳ１８へ進む。一方、上記ステップＳ１１で、φｆｕｌ＝１．２（フル増量）、あるいは、ステップＳ１２で、φｔｅｘ≠０と判断されている場合には、ステップＳ１４へ分岐し、ステップＳ１４〜Ｓ１７でエンジン回転数Ｎｅに基づき、気筒別の一次元マップを補間計算付で参照して、上記気筒別噴射量補正係数Ｋｔｎ１〜Ｋｔｎ４を各々設定して、ステップＳ１８へ進む。噴射量補正係数Ｋｔｎは燃料噴射量を気筒別に増減量して空燃比を均一に設定するもので、上記各一次元マップは、各気筒の吸気特性に対応する値が格納されている。
【０１２１】
ところで、運転条件によって燃料噴射量を気筒毎に変えなければならない原因は、本来、筒内吸入空気質量Ｇａｉｒが気筒毎に相違するからであり、従って、正確には筒内吸入空気質量Ｇａｉｒを吸気特性を考慮して気筒毎に算出し、この算出結果から筒内吸入燃料質量Ｇｆｕｅｌ及び燃料噴射質量Ｇｉｎｊを気筒毎に算出すべきであるが、ＣＰＵの計算負荷が重くなるため、上記気筒別噴射量補正係数Ｋｔｎ１〜Ｋｔｎ４を用いて、見かけ上、気筒別燃料噴射有効パルス幅Ｔｅ１〜Ｔｅ４［ｍｓ］を補正するようにした。
【０１２２】
また、空燃比の気筒毎のばらつきが特に問題となるのは、高負荷運転時のノッキングである。中低負荷運転での空燃比のばらつきは排気ガスに余り影響せず、エンジン振動の観点からすると全気筒の燃料量が一致している方が却って各気筒の出力が均等になるので望ましい。従って、フル増量が要求されている場合には、全気筒が最大出力が得られるようにする必要があり、また、排気ガス温度限界による増量が要求されている場合には、全気筒の排気系等の高温化を設計限界以内に抑制してエンジンを保護する必要があるため、このような場合、すなわち、気筒間の空燃比が均一であることを前提とする要求増量が設定されている場合のみ、エンジン回転数毎に各気筒の燃料噴射量を適正に増減させる補正係数（気筒別噴射量補正係数）を設定し、各気筒の筒内混合気の空燃比が目標空燃比となるように個別に制御する。
【０１２３】
そして、ステップＳ１８へ進むと、全体としての燃料噴射有効パルス幅Ｔｅ＿ａｌｌを次式から算出する。
【０１２４】
Ｔｅ＿ａｌｌ←Ｋｍｒ×Ｇｉｎｊ×Ｋｃｏｎ …（１９）
Ｋｍｒ：パルス幅誤差補正係数
Ｋｃｏｎ：インジェクタ容量係数［ｓｅｃ／ｇ］
ここで、パルス幅誤差補正係数Ｋｍｒは、インジェクタ２５の動的な流量特性の非線型性の補正であり、本ルーチンでは、エンジン回転数Ｎｅと燃料噴射質量Ｇｉｎｊとに基づき、図２２に示す二次元マップを補間計算付で参照して設定する。また、上記インジェクタ容量係数Ｋｃｏｎは、インジェクタ２５の静的質量流量特性の逆数であり、１グラムの燃料を０．１秒間に噴射するインジェクタでは、Ｋｃｏｎ＝０．１［ｓｅｃ／ｇ］である。なお、このインジェクタ容量係数Ｋｃｏｎは、本ルーチンでは固定値としているが、多種燃料に対応可能なＦＦＶエンジン等で、燃料の比重や粘性が変化する場合には変数として用いる。
【０１２５】
次いで、ステップＳ１９〜Ｓ２２で、上記燃料噴射有効パルス幅Ｔｅ＿ａｌｌを上記気筒別噴射量補正係数Ｋｔｎ１〜Ｋｔｎ４で補正して、気筒別燃料噴射有効パルス幅Ｔｅ１〜Ｔｅ４をそれぞれ設定する。
【０１２６】
その後、ステップＳ２３へ進み、バッテリ電圧ＶＢに基づき、一次元マップを補間計算付で参照し燃料噴射無駄パルス幅Ｔｓ＿ａｌｌ［ｍｓ］を設定し、ステップＳ２４で、気筒別燃料噴射無駄パルス幅Ｔｓ１〜Ｔｓ４を、上記燃料噴射無駄パルス幅Ｔｓ＿ａｌｌで設定してルーチンを抜ける。
【０１２７】
ＯＳ側では、上記気筒別燃料噴射有効パルス幅Ｔｅ１〜Ｔｅ４と上記気筒別燃料噴射無駄パルス幅Ｔｓ１〜Ｔｓ４とをそれぞれ加算して、燃料噴射パルス幅を気筒別に算出する。
【０１２８】
インジェクタ２５に与える電圧パルスの幅（燃料噴射パルス幅）と、燃料噴射質量Ｇｉｎｊとの関係は、図２３に示す通りであり、燃料噴射質量Ｇｉｎｊが決定されれば、基本的に、この燃料噴射質量Ｇｉｎｊにインジェクタ容量係数Ｋｃｏｎを掛けることで上記有効燃料噴射パルス幅Ｔｅが求められ、この有効燃料噴射パルス幅Ｔｅに無効噴射パルス幅である燃料噴射無駄パルス幅Ｔｓを加算することで燃料噴射パルス幅が算出される。
【０１２９】
本実施例では、燃料噴射質量Ｇｉｎｊ分の燃料を１回で噴射するか、２回に分割して噴射するかはＯＳ内部で決定されるため、上記気筒別燃料噴射有効パルス幅Ｔｅ１〜Ｔｅ４と上記気筒別燃料噴射無駄パルス幅Ｔｓ１〜Ｔｓ４とを、予め加算せずに別々にＯＳに渡し、このＯＳ内部で噴射設定の直前に加算して、各インジェクタに与える燃料噴射パルス幅を算出している。こうすることで、ＯＳ内部で、燃料噴射質量Ｇｉｎｊ分の燃料を２回に分割して噴射すると決定した場合でも、１／２の燃料噴射有効パルス幅Ｔｅに燃料噴射無駄パルス幅Ｔｓを加算することで、１回当りの燃料噴射パルス幅を簡単に与えることができる。
【０１３０】
このように、本実施例によれば、例えば新しい燃料付着モデルを開発した際は、筒内吸入空気質量Ｇａｉｒを算出する前記（３）式、或は（１０）式中に演算式を追加すれば良く、またインジェクタの流量が非線型であるならば燃料噴射有効パルス幅Ｔｅ，燃料噴射無駄パルス幅Ｔｓを補正するだけで良いことになる。その結果、ある運転領域での空燃比制御性に不具合が生じた場合に、どの部分を変更すべきかが明確となるばかりか、その変更が他の部分に影響を及ぼし難いため、ベース制御ソフトとしての活用が期待できる。
【０１３１】
【発明の効果】
以上、説明したように本発明によれば、スロットル開度とエンジン回転数とに基づき静的な筒内吸入空気質量に相当する静的指標値を設定し、静的指標値に吸気チャンバによる遅れを想定した一次遅れ処理を施して動的な筒内吸入空気質量に相当する動的指標値を設定する。そして、更に動的指標値を加重平均して燃料付着による遅れ指標値を算出し、上記動的指標値から遅れ指標値を減算して、この差分に基づき燃料追加質量を設定する。そして、基本とする吸入燃料質量に、動的指標値と遅れ指標値との差分に基づいて設定した燃料追加質量を加算して１サイクル当りの燃料噴射量を設定するので、従来の燃料付着モデルを用いた場合に比し、応答遅れや、ノイズの影響を受けることが無く、過渡時の燃料追加質量を高精度に求めることができて、信頼性の高い空燃比制御性能を得ることができる。
特に、燃料付着補正については、スロットル開度とエンジン回転数とに基づいて設定した静的な筒内吸入空気質量に相当する静的指標値に対し、一次遅れ処理を施して、更に加重平均処理を施し、一次遅れ処理により設定した動的な筒内吸入空気質量に相当する動的指標値から、該動的指標値を加重平均処理して設定した遅れ指標値を減算し、この差分に基づき燃料追加質量を設定するので、マップやテーブル、或いは、マップやテーブルに相当する補正回路等を必要とせず、スロットル開度とエンジン回転数のみから、非常に簡素な処理で、応答遅れを解消して、信頼性の高い高精度な燃料付着補正を実現できる。
【図面の簡単な説明】
【図１】１０ｍｓ毎ジョブで実行される燃料噴射有効パルス幅及び燃料噴射無駄パルス幅設定ルーチンを示すフローチャート
【図２】１０ｍｓ毎ジョブで実行される燃料噴射有効パルス幅及び燃料噴射無駄パルス幅設定ルーチンを示すフローチャート（続き）
【図３】１０ｍｓ毎ジョブで実行される燃料噴射有効パルス幅及び燃料噴射無駄パルス幅設定ルーチンを示すフローチャート（続き）
【図４】他の態様による燃料噴射有効パルス幅及び燃料噴射無駄パルス幅設定ルーチンの要部を示すフローチャート
【図５】１０ｍｓ毎ジョブで実行される筒内吸入空気質量設定ルーチンを示すフローチャート
【図６】５０ｍｓ毎ジョブで実行される係数設定ルーチンを示すフローチャート
【図７】５０ｍｓ毎ジョブで実行される要求当量比設定ルーチンを示すフローチャート
【図８】５０ｍｓ毎ジョブで実行される要求当量比設定ルーチンを示すフローチャート（続き）
【図９】５０ｍｓ毎ジョブで実行される要求当量比設定ルーチンを示すフローチャート（続き）
【図１０】他の態様による４ｍｓ毎ジョブ実行される筒内吸入空気質量設定ルーチンを示すフローチャート
【図１１】筒内吸入空気質量と吸気温度との関係を示す特性図
【図１２】理論吸気質量と筒内吸入空気質量との関係を示す特性図
【図１３】吸入空気量誤差補正係数を設定する二次元マップの概念図
【図１４】過渡時における吸入空気量センサの出力電圧を示す特性図
【図１５】基本可燃限界当量比を設定する二次元マップの概念図
【図１６】最大出力当量比の増量領域を示す特性図
【図１７】基本排気ガス温度限界当量比を設定する二次元マップの特性図
【図１８】燃費率最良当量比を設定する二次元マップの特性図
【図１９】静的な筒内吸入空気質量に相当する指標値を設定する二次元マップの概念図
【図２０】静的指標値と動的指標値とこの動的指標値に対する遅れ指標値との関係を示す特性図
【図２１】過渡運転時のリーンスパイクと増量パターンとの関係を示す特性図
【図２２】パルス幅誤差補正係数を設定する二次元マップの特性図
【図２３】燃料噴射パルス幅と燃料噴射質量との関係を示す特性図
【図２４】エンジンの吸気系モデルを示す説明図
【図２５】エンジンの全体構成図
【図２６】クランクロータとクランク角センサの正面図
【図２７】カムロータとカム角センサの正面図
【図２８】電子制御装置の回路構成図
【符号の説明】
６Ａ…吸気チャンバ
Ｇａｉｒ…吸入空気質量
Ｇａｃｃ…燃料追加質量
Ｇｆｕｅｌ…吸入燃料質量
Ｇｉｎｊ…燃料噴射質量
Ｍａｃｃ…静的指標値
Ｓａｃｃ…動的指標値
Ｔａｃｃ…遅れ指標値
α…スロットル開度
【産業上の利用分野】
本発明は、過渡時の燃料付着量を補填する燃料追加質量を、筒内に実際に吸入される燃料質量に加算して、燃料噴射の際に必要な燃料噴射質量を求めるエンジンの空燃比制御方法に関する。
【０００２】
【従来の技術】
一般に、インジェクタから噴射された燃料は、その全てが筒内に吸込まれることはなく、一部は吸気ポート内壁等に付着して燃料液膜となる。定常的には燃料付着率が一定であれば、インジェクタからの燃料噴射量と筒内に吸込まれる燃料量とはほぼ一致するが、過渡的にはインジェクタからの噴射量が急増しても、燃料付着量が相対的に増加する分、筒内へ吸入される燃料量は要求燃料量よりも少なくなる。そのため、筒内混合気の空燃比が一時的にリーン化し、上記吸気ポート内壁等に付着した燃料は、ある遅れを有して筒内に吸込まれるため、筒内への吸入燃料量はゆっくりと増加することになる。
【０００３】
従って、エンジンの空燃比制御では、このような過渡的な空燃比のリーン化を防止し、安定した出力を得るために上記燃料付着量分を補填するフィードフォワード制御が必要となる。
【０００４】
例えば、特開平５−１８２８９号公報では、理論的な燃料付着モデルに基づき、まず、１回転当りの吸入空気量Ｑ／Ｎとエンジン回転数Ｎｅとによりマップを参照して暖機完了後の飽和状態の壁面付着量ＱＭＷを補間計算し、前回求めた壁面付着量ＱＭＷとの差分から、変化量ΔＱＭＷを算出し、この変化量ΔＱＭＷを、即時に変化する即時補正と、徐々に変化するテーリング補正とで補正して、壁面付着補正量ＦＭＷを求め、この壁面付着補正量ＦＭＷを、吸入空気量Ｑとエンジン回転数Ｎｅとに基づいて求めた基本燃料噴射量ＴＡＵＰに加算し、必要な補正項で補正して、インジェクタからの噴射量ＴＡＵを算出している。
【０００５】
上記変化量ΔＱＭＷを、即時補正とテーリング補正とに分けて補正しているのは、過渡的な上記壁面付着量ＱＭＷは、計算上は瞬時に変化するが、実際には徐々に変化するためである。
【０００６】
【発明が解決しようとする課題】
上述したように、一般的には、理論的な燃料付着モデルから壁面付着補正量ＦＭＷを求める場合が多いが、吸入空気量Ｑを計測する吸入空気量センサ自体の応答遅れや、ノイズの影響等で壁面付着量ＱＭＷの検出性能にばらつきが生じ易く、充分な信頼性が得られない。
【０００７】
また、古くはスロットル弁の開度変化から過渡的な付着補正を試みる、いわゆる三角増量法がある。この方式は、スロットル弁の開度変化に基づくため応答性は良いが、燃料付着による遅れを三角形で近似的に求めているため、空燃比が部分的にリーン化したり、リッチ化したりし易く、空燃比制御精度に問題がある。
【０００８】
本発明は、上記事情に鑑みてなされたもので、応答遅れや、ノイズの影響を受けること無く、高精度の付着補正量を求めることができて、信頼性の高い過渡補正を行うことのできるエンジンの空燃比制御方法を提供することを目的としている。
【０００９】
【課題を解決するための手段】
上記目的を達成するため本発明は、筒内への１吸気行程当りの吸入空気質量を基本として目標空燃比に対応する筒内への１吸気行程当りの吸入燃料質量を設定するエンジンの空燃比制御方法において、スロットル開度とエンジン回転数とに基づき静的な筒内吸入空気質量に相当する静的指標値を設定し、上記静的指標値に吸気チャンバによる遅れを想定した一次遅れ処理を施して動的な筒内吸入空気質量に相当する動的指標値を設定し、更に上記動的指標値を加重平均して燃料付着による遅れ指標値を算出し、上記動的指標値から上記遅れ指標値を減算して、この差分に基づいて燃料追加質量を設定し、上記燃料追加質量を上記吸入燃料質量に加算して１サイクル当りの燃料噴射量を設定することを特徴とする。
【００１０】
【作用】
本発明は、スロットル開度とエンジン回転数とに基づき静的な筒内吸入空気質量に相当する静的指標値を設定し、静的指標値に吸気チャンバによる遅れを想定した一次遅れ処理を施して動的な筒内吸入空気質量に相当する動的指標値を設定する。そして、更に動的指標値を加重平均して燃料付着による遅れ指標値を算出し、上記動的指標値から遅れ指標値を減算して、この差分に基づき燃料追加質量を設定する。そして、基本とする吸入燃料質量に、動的指標値と遅れ指標値との差分に基づいて設定した燃料追加質量を加算して１サイクル当りの燃料噴射量を設定する。
【００１４】
【実施例】
以下、図面に基づき本発明の実施例を説明する。図２５にはエンジンの全体概略図が示されており、図中の符号１はエンジン（図においては水平対向４気筒型エンジンを示す）を示し、シリンダヘッド２の吸気ポート２ａにインテークマニホルド３が連通され、このインテークマニホルド３の上流にエアーチャンバ４を介してスロットル通路５が連通されている。このスロットル通路５の上流側には、吸気管６を介してエアクリーナ７が取付けられ、このエアクリーナ７が吸入空気の取り入れ口であるエアインテークチャンバ８に連通されている。また、上記排気ポート２ｂにエキゾーストマニホルド９を介して排気管１０が連通され、この排気管１０に触媒コンバータ１１が介装されてマフラ１２に連通されている。
【００１５】
一方、上記スロットル通路５にスロットル弁５ａが設けられ、このスロットル通路５の直上流の上記吸気管６にインタークーラ１３が介装され、さらに、上記吸気管６の上記エアクリーナ７の下流側にレゾネータチャンバ１４が介装されている。
【００１６】
また、上記レゾネータチャンバ１４と上記インテークマニホルド３とを連通して上記スロットル弁５ａの上流側と下流側とをバイパスするバイパス通路１５に、アイドル空気量を調整するアイドル制御（ＩＳＣ）弁１６が介装されている。さらに、上記ＩＳＣ弁１６の直下流側に、吸気圧が負圧のとき開弁し、またターボチャージャ１８によって過給されて吸気圧が正圧になったとき閉弁するチェックバルブ１７が介装されている。
【００１７】
上記ターボチャージャ１８は、上記吸気管６の上記レゾネータチャンバ１４の下流側にコンプレッサが介装され、タービンが上記排気管１０に介装されている。さらに、上記ターボチャージャ１８のタービンハウジング流入口には、ウエストゲート弁１９が介装され、このウエストゲート弁１９には、ウエストゲート弁作動用アクチュエータ２０が連設されている。
【００１８】
上記ウエストゲート弁作動用アクチュエータ２０は、ダイヤフラムにより２室に仕切られ、一方がウエストゲート弁制御用デューティソレノイド弁２１に連通される圧力室を形成し、他方が上記ウエストゲート弁１９を閉方向に付勢するスプリングを収納したスプリング室を形成している。
【００１９】
上記ウエストゲート弁制御用デューティソレノイド弁２１は、上記レゾネータチャンバ１４と上記吸気管６の上記ターボチャージャ１８のコンプレッサ下流とを連通する通路に介装されており、後述する電子制御装置５０（ＥＣＵ；図２８参照）から出力される制御信号のデューティ比に応じて、上記レゾネータチャンバ１４側の圧力と上記コンプレッサ下流側の圧力とを調圧し、上記ウエストゲート弁作動用アクチュエータ２０の圧力室に供給する。
【００２０】
すなわち、上記ＥＣＵ５０によって上記ウエストゲート弁制御用デューティソレノイド弁２１を制御し、上記ウエストゲート弁作動用アクチュエータ２０を作動させて上記ウエストゲート弁１９による排気ガスリリーフを調整することにより、上記ターボチャージャ１８による過給圧を制御するようになっている。
【００２１】
また、上記インテークマニホルド３に、吸気管圧力Ｐを絶対圧で検出する吸気管圧力センサ２２が通路２３を介して連通され、さらに、上記インテークマニホルド３の各気筒の各吸気ポート２ａの直上流側にインジェクタ２５が臨まされている。また、上記シリンダヘッド２の各気筒毎に、その先端を燃焼室に露呈する点火プラグ２６ａが取付けられ、この点火プラグ２６ａに各気筒毎に配設された点火コイル２６ｂを介してイグナイタ２７が接続されている。
【００２２】
上記インジェクタ２５には、燃料タンク２８内に設けたインタンク式の燃料ポンプ２９から燃料フィルタ３０を経て燃料が圧送され、プレッシャレギュレータ３１にてインジェクタ２５への燃料圧力が調圧される。
【００２３】
また、上記吸気管６の上記エアクリーナ７の直下流に、ホットワイヤ或はホットフィルム等を用いた熱式の吸入空気量センサ３２が介装され、上記スロットル弁５ａに、スロットル開度センサ３３ａとアイドルスイッチ３３ｂとを内蔵したスロットルセンサ３３が連設されている。さらに、上記エアーチャンバ４に吸気温センサ４６が臨まされている。
【００２４】
また、上記エンジン１のシリンダブロック１ａにノックセンサ３４が取付けられるとともに、このシリンダブロック１ａの左右両バンクを連通する冷却水通路３５に水温センサ３６が臨まされ、上記排気管１０の上記エキゾーストマニホルド９の集合部にＯ２ センサ３７が臨まされている。
【００２５】
また、上記シリンダブロック１ａに支承されたクランクシャフト１ｂにクランクロータ３８が軸着され、このクランクロータ３８の外周に、電磁ピックアップなどからなるクランク角センサ３９が対設されている。さらに、上記エンジン１のカムシャフト１ｃに連設するカムロータ４０に、電磁ピックアップなどからなる気筒判別用のカム角センサ４１が対設されている。尚、上記クランク角センサ３９及び上記カム角センサ４１は、電磁ピックアップなどの磁気センサに限らず、光センサなどでも良い。
【００２６】
上記クランクロータ３８は、図２６に示すように、その外周に突起３８ａ，３８ｂ，３８ｃが形成され、これらの各突起３８ａ，３８ｂ，３８ｃが、各気筒（＃１，＃２と＃３，＃４）の圧縮上死点前（ＢＴＤＣ）θ１，θ２，θ３ の位置に形成されており、本実施例においては、θ１ ＝９７゜ＣＡ、θ２ ＝６５゜ＣＡ、θ３ ＝１０゜ＣＡである。
【００２７】
上記クランクロータ３８の各突起は、上記クランク角センサ３９によって検出され、ＢＴＤＣ９７゜，６５゜，１０゜のクランクパルスがエンジン１／２回転毎（１８０゜ＣＡ毎）に出力される。そして、各信号の入力間隔時間がタイマによって計時され、エンジン回転数Ｎｅが算出される。
【００２８】
また、図２７に示すように、上記カムロータ４０の外周には、気筒判別用の突起４０ａ，４０ｂ，４０ｃが形成され、突起４０ａが＃３，＃４気筒の圧縮上死点後（ＡＴＤＣ）θ４ の位置に形成され、突起４０ｂが３個の突起で構成されて最初の突起が＃１気筒のＡＴＤＣθ５ の位置に形成されている。さらに、突起４０ｃが２個の突起で形成され、最初の突起が＃２気筒のＡＴＤＣθ６ の位置に形成されている。本実施例においては、θ４ ＝２０゜ＣＡ、θ５ ＝５゜ＣＡ、θ６ ＝２０゜ＣＡである。
【００２９】
そして、上記カムロータ４０の各突起が上記カム角センサ４１によって検出され、各気筒の燃焼行程順を＃１→＃３→＃２→＃４とした場合、この燃焼行程順と、上記カム角センサ４１からのカムパルスをカウンタによって計数した値とのパターンに基づいて、気筒判別がなされる。
【００３０】
一方、図２８において、符号５０は、エンジン系を制御する電子制御装置（ＥＣＵ）５０であり、このＥＣＵ５０は、燃料噴射制御、点火時期制御などを行なうメインコンピュータ５１と、ノック検出処理を行なう専用のサブコンピュータ５２との２つのコンピュータを中心として構成され、各部に所定の安定化電源を供給する定電圧回路５３や各種の周辺回路が組込まれている。
【００３１】
上記定電圧回路５３は、ＥＣＵリレー５４のリレー接点を介してバッテリ５５に接続されており、このバッテリ５５に、上記ＥＣＵリレー５４のリレーコイルがイグニッションスイッチ５６を介して接続され、上記イグニッションスイッチ５６がＯＮされて上記ＥＣＵリレー５４のリレー接点が閉となったときには制御用電源を上記各コンピュータ５１，５２に供給すると共に、上記定電圧回路５３は上記バッテリ５５に直接接続されており、イグニッションスイッチ５６のＯＮ，ＯＦＦに拘らず、バックアップＲＡＭ６１にバックアップ用の電源を供給する。
【００３２】
尚、上記バッテリ５５には、燃料ポンプリレー５７のリレー接点を介して燃料ポンプ２９が接続されている。
【００３３】
上記メインコンピュータ５１は、ＣＰＵ５８、ＲＯＭ５９、ＲＡＭ６０、バックアップＲＡＭ６１、カウンタ・タイマ群６２、シリアル通信インターフェースであるＳＣＩ６３、及び、Ｉ／Ｏインターフェース６４がバスライン６５を介して接続されたマイクロコンピュータであり、上記バックアップＲＡＭ６１には、上記イグニッションスイッチ５６のＯＮ／ＯＦＦに拘らず、バッテリ５５に直接接続する上記定電圧回路５３からバックアップ電源が常時供給されてデータが保持される。
【００３４】
尚、上記カウンタ・タイマ群６２は、フリーランカウンタ、カム角センサ信号の入力計数用カウンタなどの各種カウンタ、燃料噴射タイマ、点火タイマ、定期割込みを発生させるための定期割込みタイマ、クランク角センサ信号の入力間隔計時用タイマ、及び、システム異常監視用のウオッチドッグタイマなどの各種タイマを便宜上総称するものであり、上記メインコンピュータ５１においては、その他、各種のソフトウエアカウンタ・タイマが用いられる。
【００３５】
また、上記サブコンピュータ５２も、上記メインコンピュータ５１と同様、ＣＰＵ７１、ＲＯＭ７２、ＲＡＭ７３、カウンタ・タイマ群７４、ＳＣＩ７５、及び、Ｉ／Ｏインターフェース７６がバスライン７７を介して接続されたマイクロコンピュータであり、上記メインコンピュータ５１とサブコンピュータ５２とは、上記ＳＣＩ６３，７５を介してシリアル通信ラインにより互いに接続されている。
【００３６】
上記メインコンピュータ５１のＩ／Ｏインターフェース６４には、入力ポートに、吸入空気量センサ３２、スロットル開度センサ３３ａ、水温センサ３６、Ｏ２センサ３７、吸気管圧力センサ２２、大気圧センサ４４、車速センサ４２、吸気温センサ４６、及び、バッテリ５５が、Ａ／Ｄ変換器６６を介して接続されるとともに、アイドルスイッチ３３ｂ、スタータスイッチ４３、及びクランク角センサ３９、カム角センサ４１などが接続され、さらに、図示しない各種センサ及びスイッチ類が接続されている。
【００３７】
また、上記Ｉ／Ｏインターフェース６４の出力ポートには、イグナイタ２７が接続されるとともに、ＩＳＣ弁１６、インジェクタ２５、燃料ポンプリレー５７のリレーコイル、ウエストゲート弁制御用デューティソレノイド弁２１が駆動回路６７を介して接続されており、さらに、図示しない各種のアクチュエータ類が接続されている。
【００３８】
一方、上記サブコンピュータ５２のＩ／Ｏインターフェース７６は、入力ポートに、クランク角センサ３９、カム角センサ４１が接続されるとともに、Ａ／Ｄ変換器７８、周波数フィルタ７９、アンプ８０を介してノックセンサ３４が接続されており、上記ノックセンサ３４からのノック検出信号が上記アンプ８０で所定のレベルに増幅された後に上記周波数フィルタ７９により必要な周波数成分が抽出され、上記Ａ／Ｄ変換器７８にてデジタル信号に変換されて入力されるようになっている。
【００３９】
上記メインコンピュータ５１では、各センサ類からの検出信号を処理し、燃料噴射量制御、点火時期制御、アイドル制御などを行い、一方、上記サブコンピュータ５２では、エンジン回転数とエンジン負荷とに基づいてノックセンサ３４からの信号のサンプル区間を設定し、このサンプル区間でノックセンサ３４からの信号を高速にＡ／Ｄ変換して振動波形を忠実にデジタルデータに変換し、このデータに基づきノック発生の有無を判定する。
【００４０】
上記サブコンピュータ５２のＩ／Ｏインターフェース７６の出力ポートは、上記メインコンピュータ５１のＩ／Ｏインターフェース６４の入力ポートに接続されており、上記サブコンピュータ５２でのノック判定結果がＩ／Ｏインターフェース７６に出力される。そして、上記メインコンピュータ５１では、上記サブコンピュータ５２からノック発生有りの判定結果が出力されると、ＳＣＩ６３を介してシリアル通信ラインよりサブコンピュータ５２からノックデータを読込み、このノックデータに基づいて直ちに該当気筒の点火時期を遅らせ、ノックを回避する。
【００４１】
このようなエンジン制御において、上記メインコンピュータ５１では、センサ・スイッチ類からの信号入力処理、燃料噴射制御、点火時期制御、アイドル制御に係わる各種ジョブが、一つのオペレーティングシステム（ＯＳ）の管理下で効率的に実行される。このＯＳは、車輌制御のための各種マネジメント機能、及び、このマネジメント機能に密着した内部ストラテジーを有し、各種ジョブを体系的に結合し、等時間間隔処理により各種ジョブを効率的に実行するようになっている。
【００４２】
以下、上記メインコンピュータ５１による燃料噴射制御について、図１〜図９のルーチンに従って説明する。尚、サブコンピュータ５２はノック検出処理専用のコンピュータであるため、その動作説明を省略する。
【００４３】
本実施例では、ＯＳ側で処理したＡ／Ｄ変換結果、クランク位置情報、エンジン回転数等に基づき、ユーザ側ジョブで、燃料噴射量、点火時期等を設定すると共に、これらの指示値がＯＳによって噴射タイマ、点火タイマにセットされる。そして、１０ｍｓ毎ジョブ，５０ｍｓ毎ジョブ等の定期割込み要求を出力し、またクランク角センサ３９からの信号入力毎（ＢＴＤＣ９７゜，６５゜，１０゜ＣＡ毎のエンジン１回転に６回）に起動されるクランクパルス入力毎ジョブの割込みを許可する。
【００４４】
１０ｍｓ毎ジョブでは、図１〜図３に示す燃料噴射有効パルス幅及び燃料噴射無駄パルス幅設定ルーチンと、図５に示す吸入空気量設定ルーチンが実行される。また、５０ｍｓ毎ジョブでは、図６に示す係数設定ルーチン、及び図７〜図９に示す要求当量比設定ルーチンが実行される。
【００４５】
尚、以下の説明では、図２４に示すように、吸気系のスロットル弁５ａの下流から吸気弁上流までを吸気チャンバ６Ａと総称する。従って、吸気チャンバ６Ａは、スロットル通路５、エアーチャンバ４、インテークマニホルド３、吸気ポート２ａの総称となる。
【００４６】
ここで、以下の説明で使用する用語について簡単に説明する。Ｑは質量流量［ｍｇ／ｓｅｃ］、Ｇは１サイクル当りの質量［ｍｇ／ｃｙｃｌｅ］を意味する。ただし、後述する筒内吸入空気質量Ｇａｉｒ、筒内吸入燃料質量Ｇｆｕｅｌ等は、吸気行程時にのみ吸入されるので、実質的には、１吸気行程当りの質量と同じ意味になる。また、図２４の左から、Ｑａは吸入空気量センサ３２を通過する空気の質量流量（センサ通過空気質量流量）、Ｍは吸気チャンバ６Ａの空気質量、Ｐは吸気チャンバ６Ａの吸気管内圧力、Ｇｉｎｊは燃料噴射質量、Ｑｃは筒内吸入空気質量流量、Ｇａｉｒは筒内吸入空気質量、Ｇｆｕｅｌは筒内吸入燃料質量を示す。
【００４７】
まず、図１〜図３に示す燃料噴射有効パルス幅及び燃料噴射無駄パルス幅設定ルーチンを説明する前に、このルーチンで取入れる各パラメータを設定するルーチンについて説明する。
【００４８】
図５の１０ｍｓ毎ジョブで実行されるルーチンでは、Ｄジェトロニック方式により１吸気行程当りの筒内吸入空気質量Ｇａｉｒを算出する。まず、ステップＳ４１で、吸気管圧力センサ２２の検出値に基づいて算出したスロットル弁５ａの下流側の吸気管圧力Ｐと、吸気温センサ４６の検出値に基づいて算出した吸気温度Ｔとを読込み、ステップＳ４２で、現在の吸気温度Ｔが５０℃以下、すなわち通常の吸気温度領域にあるか判断し、Ｔ＜５０℃のときには、ステップＳ４３へ進み、現在の吸気温度Ｔで、吸気温度設定値Ｔ’を設定してステップＳ４５へ進む。また、Ｔ≧５０℃のときには、ステップＳ４４で、上記吸気温度設定値Ｔ’を固定値である５０℃に設定した後、ステップＳ４５へ進む。
【００４９】
ところで、筒内吸入空気質量Ｇａｉｒは、上記吸気管圧力Ｐが一定であっても吸気温度Ｔに影響されて変動するため、空気密度σから筒内吸入空気質量Ｇａｉｒを算出する必要がある。筒内の空気密度σが吸気行程時における吸気チャンバ６Ａ内の空気密度と同じであると考えれば、上記筒内の空気密度αは、空気の状態方程式により以下の式から求めることができる。
【００５０】
σ＝Ｐ×（１−Ｒｅｇｒ）／（Ｒ×Ｔ） …（１）
Ｒｅｇｒ；ＥＧＲ率
Ｒ；ガス定数
図１１に筒内吸入空気質量Ｇａｉｒと吸気温度Ｔとの関係を示す。実線は実測により得た吸入空気質量（実測吸気量）の特性、二点鎖線は上記（１）式により算出した空気密度σに従い求めた吸入空気質量（密度比例吸気量）、一点鎖線は吸気温度Ｔをある温度に固定して算出した空気密度σに従い求めた吸入空気質量（圧力比例吸気量）を示す。図に示すように、吸気温度Ｔが５０℃以下の通常の吸気温度領域では、密度比例吸気量が実測吸気量にほぼ沿う特性を示す。しかし、ターボチャージャ付エンジンでの実測結果によると、吸気温度Ｔが５０℃を越えた領域では、上記筒内吸入空気質量Ｇａｉｒは空気密度σに比例して吸込まれておらず、密度比例吸気量は、実測吸気量に比し過小になり、結果として、空燃比がリーン化してしまう。ところで、本来なら、吸気温度Ｔが５０℃以上の領域では、温度毎の補正係数を用いて補正すべきであるが、制御が複雑化するのと、圧力比例吸気量がほぼ実測吸気量に沿う特性を示しているので、本ルーチンでは吸気温度設定値Ｔ’を５０℃に固定して、筒内吸入空気質量Ｇａｉｒを求めるようにした。尚、上記５０℃は実験により求めた臨界温度であり、採用するエンジン形式により相違するため、個別に設定する必要がある。また、本実施例のようにＥＧＲ装置が搭載されていないエンジンでは、上記ＥＧＲ率はゼロである。
【００５１】
そして、ステップＳ４５へ進むと、上記（１）式の吸気温度Ｔを上記吸気温度設定値Ｔ’に代えて空気密度σを算出し、ステップＳ４６で、上記空気密度σと行程容積Ｄとに基づき次式から１吸気行程当りの理論吸気質量Ｇｔｈを算出する。
【００５２】
Ｇｔｈ←Ｄ×σ …（２）
従って、この理論吸気質量Ｇｔｈは、吸気温度設定値Ｔ’を考慮しなければ、上記吸気管圧力Ｐに比例する。
【００５３】
次いで、ステップＳ４７で、上記理論吸気質量Ｇｔｈを基本として筒内吸入空気質量Ｇａｉｒを次式から算出して、ルーチンを抜ける。
【００５４】
Ｇａｉｒ←（Ｇｔｈ−ηｂ）×ηｖ×Ｋｔｒｍ …（３）
ηｂ：吸気損失質量
ηｖ：体積効率
Ｋｔｒｍ：吸入空気量誤差補正係数
上記理論吸気質量Ｇｔｈと筒内吸入空気質量Ｇａｉｒとの関係は、図１２に示す通りであり、理論吸気質量Ｇｔｈがゼロ、すなわち吸気管圧力Ｐが完全な真空になる前に、上記筒内吸入空気質量Ｇａｉｒがゼロになるような、横軸接点ηｂ、傾きηｖの一次式で表される。但し、実際には、筒内吸入空気質量Ｇａｉｒがゼロになるような運転条件は存在しない。また、１吸気行程当りの吸気損失質量ηｂと体積効率ηｖとは、後述するように、エンジン回転数Ｎｅに基づき一次元マップを補間計算付で参照して設定する。さらに、上記吸入空気量誤差補正係数Ｋｔｒｍは、上記理論吸気質量Ｇｔｈと筒内吸入空気質量Ｇａｉｒとの関係が完全には一次式に乗らないこともあるため、運転条件による空気量計測の定常的な誤差分を補正するためのもので、例えば、スロットル開度αとエンジン回転数Ｎｅとに基づき、二次元マップを補間計算付で参照して設定する。図１３に、この吸入空気量誤差補正係数Ｋｔｒｍを設定する二次元マップの特性を例示する。この図に示すように、この吸入空気量誤差補正係数Ｋｔｒｍの値は基本的には、１．０である。
【００５５】
次に、図６に示す５０ｍｓ毎ジョブで実行される係数設定ルーチンについて説明する。このルーチンでは、上記吸気損失質量ηｂ、及び上記体積効率ηｖを設定する。
【００５６】
まず、ステップＳ５１でエンジン回転数Ｎｅを読込み、ステップＳ５２、Ｓ５３で、このエンジン回転数Ｎｅに基づき一次元マップを補間計算付で参照して、体積効率ηｖと吸気損失質量ηｂとを各々設定してルーチンを抜ける。
【００５７】
上記体積効率ηｖ、上記吸気損失質量ηｂの値は、理論的には、１．０であるが、エンジン回転数毎にカムの同調などの影響によって変化する値であり、また上記体積効率ηｖには、高地補正、排気圧補正等の因子を加味させることができる。
【００５８】
このように、本ルーチンでは、吸気管圧力Ｐを基本に筒内吸入空気質量Ｇａｉｒを簡単な計算式に基づいて算出し、また、上記体積効率ηｖ、上記吸気損失質量ηｂは、一次元マップから設定することができるため、コンピュータに掛かる負荷が軽くなり、しかも、上記吸気管圧力Ｐさえ正確に計測されれば、始動時を含む全ての領域で適用することができる。
【００５９】
ところで、上記筒内吸入空気質量Ｇａｉｒは、Ｌジェトロニック方式により算出することもできる。このＬジェトロニック方式を利用して筒内吸入空気質量Ｇａｉｒを算出する手順を、図１０の４ｍｓ毎ジョブで実行されるルーチンに従って説明する。尚、上記吸入空気量センサ３２の検出時刻は、スロットル弁５ａまでの吸気管路長、及び流速によっても異なるが、スロットル弁５ａでの検出時刻とほぼ同じと考えて良い。
【００６０】
まず、ステップＳ１０１で、吸入空気量センサ３２からの出力電圧ＶＡＦＭを読込み、次のステップＳ１０２，Ｓ１０３で上記出力電圧ＶＡＦＭの過渡補正をする。すなわち、図１４に示すように、実験によれば、過渡時の吸入空気量の急激な変化に対して、吸入空気量センサ３２の出力電圧ＶＡＦＭの変化は、理想の変化の２／３の電圧までは充分に速い応答を示すが、残りの１／３は、約２００ｍｓの時定数を持つ遅れで応答することが明らかにされている。このときのセンサ通過空気質量流量Ｑａと出力電圧ＶＡＦＭとの関係を、時定数τ（ｓ）の伝達関数で表せば、
ＶＡＦＭ＝（２／３）×Ｑａ＋｛１／（３τ（ｓ）＋３）｝×Ｑａ …（４）
となり、応答遅れは、この（４）式の逆関数で補正すれば良く、この補正式は、
Ｑａ＝（１／２）×［３−１／｛（２／３）τ（ｓ）＋１｝］×ＶＡＦＭ …（５）
となる。
【００６１】
上記時定数τ（ｓ）は２００ｍｓであるため、（２／３）τ（ｓ）＋１≒１３０とすれば、上記（５）式は、
【００６２】
このような考えに基づき、上記ステップＳ１０２では、前回求めた一次遅れ出力電圧ＶＡＦＭａｖと今回の出力電圧ＶＡＦＭとから、次式に示す所定時定数（本ルーチンでは、１３０ｍｓの近似値として３２×４ｍｓ）の加重平均に基づいて、今回の一次遅れ出力電圧ＶＡＦＭａｖを算出する。
【００６３】
ＶＡＦＭａｖ←（３１・ＶＡＦＭａｖ＋ＶＡＦＭ）／３２ …（７）
次いで、ステップＳ１０３で、今回の出力電圧ＶＡＦＭ及び一次遅れ出力電圧ＶＡＦＭａｖから、過渡補正後の出力電圧（過渡補正出力電圧）Ｖ’ＡＦＭを次式から算出する。
【００６４】
Ｖ’ＡＦＭ←（３・ＶＡＦＭ−ＶＡＦＭａｖ）／２ …（８）
このように、上記吸入空気量センサ３２の過渡応答遅れを簡単な一次式で導き出すことができるため、コンピュータの負荷を軽減させることができる。
【００６５】
そして、ステップＳ１０４で、過渡補正された上記過渡補正出力電圧Ｖ’ＡＦＭを、等間隔３２格子のワード補間テーブルを用いてセンサ通過空気質量流量Ｑａに変換し、ステップＳ１０５で、このセンサ通過空気質量流量Ｑａに吸入空気量誤差補正係数Ｋｔｒｍを乗算して、今回のセンサ通過空気質量流量Ｑ’ａを算出する。
【００６６】
この吸入空気量誤差補正係数Ｋｔｒｍは、運転条件による空気量計測の定常的な誤差分を補正する係数で、前述のＤジェトロニック方式と同様に、例えば、スロットル開度αとエンジン回転数Ｎｅとに基づき、二次元マップを補間計算付で参照して設定する（図１３参照）。
【００６７】
次いで、ステップＳ１０６で、エンジンが始動後かをエンジン回転数Ｎｅに基づいて判断し、Ｎｅ＜３００ｒｐｍのときは、始動時と判断してステップＳ１０７へ進み、Ｎｅ≧３００ｒｐｍのときは、始動後と判断してステップＳ１０８へ進む。
【００６８】
ステップＳ１０７へ進むと、始動時の筒内吸入空気質量Ｇａｉｒを初期値設定してルーチンを抜ける。
【００６９】
始動時に、筒内吸入空気質量Ｇａｉｒを初期値設定するのは、以下の理由による。
【００７０】
１）筒内吸入空気質量Ｇａｉｒは、定常的には、空気流量／エンジン回転数で算出されるが、エンジン回転数及び空気流量が共に小さい場合を考えると、計算結果は不定になり易い。
【００７１】
２）後述のステップＳ１０８で説明するように、筒内吸入空気質量Ｇａｉｒは、始動後は履歴に依存して逐次計算されるため、始動後最初のルーチン実行時に初期値が設定されている必要がある。
【００７２】
３）始動時に、ある始動時噴射幅を設定するのではなく、筒内吸入空気質量Ｇａｉｒを初期値設定することで、始動時、始動後に拘らず全ての運転領域で筒内吸入空気質量Ｇａｉｒを利用した演算が可能になる。
【００７３】
そこで、始動時には、筒内吸入空気質量Ｇａｉｒを推定した値で、初期値設定する。この初期値は、例えば、吸気チャンバ６Ａ内の吸気管圧力Ｐを大気圧とみなし、また吸気温度を冷却水温度により代用して次式から算出する。
【００７４】
初期値＝固定値×大気圧／冷却水温（゜Ｋ） …（９）
一方、始動後と判断されてステップＳ１０８へ進むと、上記ステップＳ１０５で算出したセンサ通過空気質量流量Ｑ’ａに基づき、筒内吸入空気質量Ｇａｉｒを次式に示す吸気チャンバモデルから算出して、ルーチンを抜ける。
Ｇａｉｒ←Ｑ’ａ×Ｍｔｃｈ＋（ＭＴＣＳ−Ｍｔｃｈ）×Ｇａｉｒ／ＭＴＣＳ …（１０）
ＭＴＣＳ：半回転時間
この（１０）式は、以下の考えに基づく。すなわち、吸気チャンバモデルの計算は、次の２式から算出することができる。
【００７５】
１）Ｇａｉｒ＝σ×Ｄ＝Ｍ×Ｄ／Ｖ
σ：吸気チャンバ６Ａ内の空気密度
Ｄ：行程容積
Ｖ：吸気チャンバ６Ａ内の容積（チャンバ容積）
２）Ｍ（ｎｅｗ）＝Ｇａｉｒ（ｏｌｄ）＋Ｑａ×Δｔ−Ｑｃ×Δｔ
（ｎｅｗ）：現在の値
（ｏｌｄ）：Δｔ前の値
Δｔ：演算周期
上記２式からＭとＱｃとを消去すると、次式に示すように、上記センサ通過空気質量流量Ｑａから筒内吸入空気質量Ｇａｉｒを算出することができる。
Ｇａｉｒ（ｎｅｗ）＝Ｇａｉｒ（ｏｌｄ）＋Ｑａ×Δｔ×Ｄ／Ｖ−Ｑｃ×Δｔ×Ｄ／Ｖ …（１０−１）
上記筒内吸入空気質量流量Ｑｃは、単位時間当りの吸気行程数が４サイクルエンジンの場合、気筒数×Ｎｅ／２であるため、
Ｑｃ＝（気筒数／２）×Ｎｅ×Ｇａｉｒ
となり、上記（１０−１）式は、
となる。ここで、Δｔ×Ｄ／Ｖ＝固定値＝Ｍｔｃｈとおき、また、採用するエンジンが４気筒であれば、上記（１０−２）式は、
Ｇａｉｒ（ｎｅｗ）＝Ｑａ×Ｍｔｃｈ＋（１−２Ｎｅ×Ｍｔｃｈ）×Ｇａｉｒ（ｏｌｄ）…（１０−３）
となる。尚、本ルーチンでは、計算の便宜上、エンジン回転数Ｎｅ［ｒｐｓ］ではなく半回転数時間（Ｎｅ／２）を用いており、この半回転数時間をＭＴＣＳとおけば、上記（１０−３）式は、
Ｇａｉｒ（ｎｅｗ）＝Ｑａ×Ｍｔｃｈ＋Ｇａｉｒ（ｏｌｄ）×（ＭＴＣＳ−Ｍｔｃｈ）／ＭＴＣＳ …（１０’）
となり、上記ステップＳ１０８で示した（１０）式が導き出される。
【００７６】
そして、この（１０）式を演算周期Δｔ毎に行えば、筒内吸入空気質量Ｇａｉｒを履歴に依存して逐次的に算出することができる。
【００７７】
このように、本実施例では、Ｄジェトロニック方式とＬジェトロニック方式との違いは、単に吸気管圧力センサ２２の検出値に基づいて筒内吸入空気質量Ｇａｉｒを算出するか、吸入空気量センサ３２の検出値に基づいて筒内吸入空気質量Ｇａｉｒを算出するか、すなわち、この筒内吸入空気質量Ｇａｉｒの算出方式の違いに過ぎず、以下のルーチンでは、筒内吸入空気質量Ｇａｉｒが算出されていれば、Ｄジェトロニック方式とＬジェトロニック方式とのいずれを採用するエンジンであっても共通となる。
【００７８】
次に、運転条件による筒内混合気の要求当量比設定ルーチンについて、図７〜図９のルーチンに従って説明する。このルーチンは５０ｍｓ毎のジョブで実行され、定常的な要求増量である可燃限界当量比φｔｗと、最大出力当量比φｆｕｌと、排気ガス温度限界当量比φｔｅｘと、触媒浄化率最良当量比φｇａｓと、燃費率最良当量比φｅｃｏとを設定する。これら各要求当量比φｔｗ，φｆｕｌ，φｔｅｘ，φｇａｓ，φｅｃｏは定常的な要求増量であるため、過渡的な燃料遅れ等によって生じる空燃比変動をフィードフォワード的に補正する係数（エアコン増量補正係数、Ｄレンジ増量補正係数等）、或はセンサやインジェクタの誤差を補正する係数（混合比補正係数等）は含まない。また、このルーチンで設定する各要求当量比は、後述する図１のルーチン実行時に比較され、その最大値が目標当量比（理論燃空比Ｆ／Ａに対する増量係数）ＣＯＥＦとして採用される。尚、この当量比φは、空気過剰率（（Ａ／Ｆ）実際／（Ａ／Ｆ）理論）λに対して、φ＝１／λの関係にある。
【００７９】
まず、ステップＳ６１で冷却水温Ｔｗとエンジン回転数Ｎｅとを読込む。そして、ステップＳ６２〜Ｓ６５或はＳ６６で可燃限界当量比φｔｗを設定する。ステップＳ６２では、上記冷却水温Ｔｗとエンジン回転数Ｎｅとに基づき二次元マップを補間計算付で参照して、可燃限界当量比φｔｗを設定する際の基準となる基準可燃限界当量比Ｍｔｗを設定する。この基準可燃限界当量比Ｍｔｗは、冷却水温Ｔｗとの関係によるリーン側の可燃限界を示すもので、図１５に、この基準可燃限界当量比Ｍｔｗを設定する二次元マップの特性を例示する。この図１５に示すように、冷却水温Ｔｗ側の軸では、従来の水温増量に近い設定がなされており、また、ある程度の高水温では、上記基準可燃限界当量比Ｍｔｗはリーン当量比に設定される。また、エンジン回転数Ｎｅ側の軸では、クランキング回転数を含む全ての回転数領域を検索できるように設定されている。筒内吸入空気質量Ｇａｉｒは始動時においても設定されることから（図５或は図１０のフローチャート参照）、このときの可燃限界当量比φｔｗが始動時の噴射パルス幅を決定する要因になる。また、始動時において、上記基準可燃限界当量比Ｍｔｗは大きな値に設定されているが、これが始動後にエンジン回転数の増加に従って小さな値に変化して行く際に、後述するように、ある一次遅れ時定数で上記可燃限界当量比φｔｗの変化を遅らせることで始動後の増量が確保される。
【００８０】
次いで、ステップＳ６３へ進むと、エンジン回転数等からエンスト状態にあるかを判断し、エンスト状態にあると判断した場合には、ステップＳ６６へ分岐し、上記基準可燃限界当量比Ｍｔｗで今回の可燃限界当量比φｔｗを設定してステップＳ６７へジャンプする。一方、エンスト状態にないと判断した場合には、ステップＳ６４へ進み、上記基準可燃限界当量比Ｍｔｗと前回の可燃限界当量比φｔｗとを比較し、Ｍｔｗ≧φｔｗのときは、要求当量比が増量されているので、ステップＳ６６へ分岐し、この基準可燃限界当量比Ｍｔｗで今回の可燃限界当量比φｔｗを設定して、増量要求に対応させた後、ステップＳ６７へジャンプする。一方、Ｍｔｗ＜φｔｗのときは、要求当量比が減少されているため、ステップＳ６５へ進み、次式に示す所定時定数（本ルーチンでは、１６×５０ｍｓ）の加重平均により一次遅れ処理を行い、今回の可燃限界当量比φｔｗを設定する。その結果、始動後増量が実質的に確保される。
【００８１】
φｔｗ←（１５φｔｗ＋Ｍｔｗ）／１６ …（１１）
次いで、ステップＳ６７〜Ｓ７０或はＳ７１で上記最大出力当量比φｆｕｌが設定される。まず、ステップＳ６７では、スロットル開度αを読込み、ステップＳ６８で、エンジン回転数Ｎｅに基づき、一次元マップを補間計算付で参照してフル増量基準値αｔｈを設定する。このフル増量基準値αｔｈは、エンジンが最大出力を要求しているかを判断する基準値で、図１６に示すように、エンジン回転数Ｎｅとスロットル開度αとの関係から特性を求めてマップ化したものである。
【００８２】
そして、ステップＳ６９で上記スロットル開度αと上記フル増量基準値αｔｈとを比較し、α≦αｔｈのときは、フル増量条件不成立と判断し、ステップＳ７０へ進み、最大出力当量比φｆｕｌを０として、ステップＳ７２へ進む。一方、α＞αｔｈのときは、フル増量条件成立と判断し、ステップＳ７１へ分岐し、上記最大出力当量比φｆｕｌを１．２、すなわち、当該エンジンの最大出力となる当量比（フル増量）に設定する。
【００８３】
尚、目標当量比ＣＯＥＦは、本ルーチンにより設定する各要求当量比φｔｗ，φｆｕｌ，φｔｅｘ，φｇａｓ，φｅｃｏの中から最大値を選択的に設定するものであるため、上記ステップＳ７０で最大出力当量比φｆｕｌが０に設定されると、この最大出力当量比φｆｕｌが目標当量比ＣＯＥＦの選択対象から外れることになる。また、排気ガス温度の上昇を抑える目的の増量は、排気ガス温度限界当量比φｔｅｘとして別に設定されているため、上記最大出力当量比φｆｕｌはエンジン回転数Ｎｅによらず、理論空燃比Ａ／Ｆ（＝１４．６）に対して、おおよそ１５〜２０％増量となるはずである。
【００８４】
次いで、ステップＳ７２〜Ｓ７９或はＳ８０で排気ガス温度限界当量比φｔｅｘを設定する。まず、ステップＳ７２で、負荷の一例である筒内吸入空気質量Ｇａｉｒとエンジン回転数Ｎｅと基づき二次元マップを補間計算付で参照して基本排気ガス温度限界当量比Ｍｔｅｘを設定する。この基本排気ガス温度限界当量比Ｍｔｅｘは、排気ガス温度の上昇を設計限界以下に抑え、エンジン及び排気系を保護するために燃料冷却を行う増量である。図１７に、この基本排気ガス温度限界当量比Ｍｔｅｘを設定する二次元マップの特性を例示する。図に示すように、高負荷（高Ｇａｉｒ）、高回転で大きな増量を行うように設定されている。尚、エンジン負荷を検出するパラメータとして上記筒内吸入空気質量Ｇａｉｒに代えて吸気管圧力Ｐを採用し、二次元マップは吸気管圧力Ｐとエンジン回転数Ｎｅとの関係による特性に基づいて設定するようにしても良い。
【００８５】
そして、ステップＳ７３へ進むと、上記基本排気ガス温度限界当量比Ｍｔｅｘと、上記最大出力当量比φｆｕｌのフル増量値である１．２とを比較し、Ｍｔｅｘ≧１．２のときはステップＳ７４へ進み、またＭｔｅｘ＜１．２のときはステップＳ７６へ進む。
【００８６】
そして、上記ステップＳ７６へ進むと、今回の平均排気ガス温度限界当量比φｔｅｘＡＶを次式に示す加重平均により、所定時定数（本ルーチンでは、１６×５０ｍｓ）の一次遅れ処理を行い、ステップＳ７８へ進む。
【００８７】
φｔｅｘＡＶ←（１５・φｔｅｘＡＶ＋Ｍｔｅｘ）／１６ …（１２）
また、ステップＳ７４へ進むと、上記基本排気ガス温度限界当量比Ｍｔｅｘと前回までの平均排気ガス温度限界当量比φｔｅｘＡＶとを比較し、Ｍｔｅｘ＜φｔｅｘＡＶのとき、すなわち、１．２≦Ｍｔｅｘ＜φｔｅｘＡＶのときは、ステップＳ７５へ進み、またＭｔｅｘ≧φｔｅｘＡＶのときは、ステップＳ７７へ進む。
【００８８】
上記ステップＳ７５では、今回の平均排気ガス温度限界当量比φｔｅｘＡＶを次式に示す加重平均により所定時定数（本ルーチンでは、３２×５０ｍｓ）の一次遅れ処理を行い、ステップＳ７８へ進む。
【００８９】
φｔｅｘＡＶ←（３１・φｔｅｘＡＶ＋Ｍｔｅｘ）／３２ …（１３）
一方、ステップＳ７７へ進むと、今回の平均排気ガス温度限界当量比φｔｅｘＡＶを次式に示す加重平均により所定時定数（本ルーチンでは１２８×５０ｍｓ）の一次遅れ処理を行って、ステップＳ７８へ進む。
【００９０】
φｔｅｘＡＶ←（１２７・φｔｅｘＡＶ＋Ｍｔｅｘ）／１２８ …（１４）
従って、Ｍｔｅｘ＜１．２のいわゆる要求当量比が薄い場合には、比較的早い時定数で平均排気ガス温度限界当量比φｔｅｘＡＶが設定され、また、Ｍｔｅｘ≧１．２、且つＭｔｅｘ≧φｔｅｘＡＶの要求当量比が濃い場合には、ゆっくりとした時定数で上記平均排気ガス温度限界当量比φｔｅｘＡＶが設定される。
【００９１】
そして、ステップＳ７８へ進むと、上記ステップＳ７５，Ｓ７６或はＳ７７で算出した平均排気ガス温度限界当量比φｔｅｘＡＶと上記基本排気ガス温度限界当量比Ｍｔｅｘとを比較し、Ｍｔｅｘ＜φｔｅｘＡＶのときは、ステップＳ７９で、上記基本排気ガス温度限界当量比Ｍｔｅｘを今回の排気ガス温度限界当量比φｔｅｘとして設定し、また、Ｍｔｅｘ≧φｔｅｘＡＶのときは、ステップＳ８０で、今回の排気ガス温度限界当量比φｔｅｘを平均排気ガス温度限界当量比φｔｅｘＡＶで設定する。
【００９２】
従って、排気ガス温度限界当量比φｔｅｘは平均排気ガス温度限界当量比φｔｅｘＡＶで上限リミットされ、また、基本排気ガス温度限界当量比Ｍｔｅｘで下限リミットされる。その結果、上記排気ガス温度限界当量比φｔｅｘの上限リミット値は、エンジン回転数Ｎｅと筒内吸入空気質量Ｇａｉｒとに応じて、すなわち、上記ステップＳ７５，Ｓ７６或はＳ７７で設定される平均排気ガス温度限界当量比φｔｅｘＡＶに従って、ゆっくり変化したり、あるいは比較的早く変化したりするようになる。その結果、長時間の全開運転において排気ガス温度が実際に設計限界値に近づくまではφｔｅｘ＜φｆｕｌとなり、従来のように、最大出力当量比φｆｕｌの中に、燃料冷却成分を含ませていた場合に比し、出力の増加が期待できる。尚、この排気ガス温度限界当量比φｔｅｘを排気ガス温度を計測或は推測した値に基づいて設定するようにしても良い。
【００９３】
次いで、上記ステップＳ７９或はＳ８０からステップＳ８１へ進み、このステップＳ８１〜Ｓ８５或はＳ８６で触媒浄化率最良当量比φｇａｓが設定される。この触媒浄化率最良当量比φｇａｓは、触媒コンバータ１１として３元触媒が採用されている場合に対象となり、排気ガスを浄化しなければならない条件が成立したときに、φｇａｓ←１．０とし、それ以外では、φｇａｓ←０とする。
【００９４】
ステップＳ８１〜Ｓ８４では、排気ガスを浄化しなければならない運転条件にあるかを判断する。この判断条件は、スロットル開度αとエンジン回転数Ｎｅとで行い、スロットル開度αが、１０゜≦α＜７０゜（ステップＳ８１，Ｓ８２）で、且つエンジン回転数Ｎｅが、８００ｒｐｍ≦Ｎｅ＜３０００ｒｐｍ（ステップＳ８３，Ｓ８４）のとき、触媒浄化率最良当量比φｇａｓを１．０に設定し（ステップＳ８５）、それ以外のときは、φｇａｓを０に設定する（ステップＳ８６）。尚、上記判断条件は一例であり、排気ガス浄化領域はエンジンの特性等に応じて任意に設定でき、また例えばスロットル開度αに代えて、吸気管圧力Ｐや筒内吸入空気質量Ｇａｉｒを採用しても良い。また、上記触媒コンバータ１１としてリーンＮｏｘ触媒が採用されている場合には、触媒浄化率最良当量比φｇａｓは全運転領域でφｇａｓ←０に設定する。
【００９５】
その後、ステップＳ８７で、負荷の一例としての吸気管圧力Ｐとエンジン回転数Ｎｅとに基づいて二次元マップを補間計算付で参照して、燃費率最良当量比φｅｃｏを設定して、ルーチンを抜ける。
【００９６】
図１８に、この燃費率最良当量比φｅｃｏを設定する二次元マップの特性を例示する。図に示すように、この二次元マップの各領域には、当該エンジン運転条件下での最良の燃費率を得ることのできる当量比が予め実験などから求めて格納されている。
【００９７】
ここで、この燃費率最良当量比φｅｃｏと前述の可燃限界当量比φｔｗとの関係について簡単に説明する。
【００９８】
空燃比を可燃限界当量比φｔｗまでリーン化させると一般には最良の燃費を得ることができなくなる。しかし、暖機途中の低水温時などでは空燃比を燃費率最良当量比φｅｃｏまでリーン化するよりも、上記可燃限界当量比φｔｗにより、リーン化を抑制した方が良好な燃焼を得ることができる場合もあるので、上記可燃限界当量比φｔｗとは別に燃費率最良当量比φｅｃｏが設定されている。尚、上記触媒浄化率最良当量比φｇａｓが全運転領域で、φｇａｓ←１．０の場合には、上記燃費率最良当量比φｅｃｏは隠れてしまうため、その意味を失うが、例えば、排気浄化システムがリーンＮｏｘ触媒を採用している場合には、前述のように、触媒浄化率最良当量比φｇａｓが全運転領域でφｇａｓ←０に設定されるため、この燃費率最良当量比φｅｃｏは触媒浄化率最良当量比としての要素も合わせ持つことになる。また、ＥＧＲを用いて理論空燃比制御により燃費向上を行う際には、Ｄジェトロニック方式では、上記（１），（２）式から明らかなように、理論吸気質量Ｇｔｈには、ＥＧＲ量が除かれた新気量のみであり、またＬジェトロニック方式では、吸入空気量センサ３２の出力電圧ＶＡＦＭに基づいて算出した筒内吸入空気質量Ｇａｉｒは、結果としてＥＧＲ量が除かれた新気量のみとなるため、全運転領域で、φｅｃｏ←１．０としたほうが、空燃比制御性が良くなる。また、この燃費率最良当量比φｅｃｏを、吸気管圧力Ｐに代えて筒内吸入空気質量Ｇａｉｒと、エンジン回転数Ｎｅとの関係により作成した二次元マップを参照して設定するようにしても良い。
【００９９】
そして、上記各ルーチンで求めた設定値を用いて、図１〜図３の１０ｍｓ毎ジョブで実行される燃料噴射有効パルス幅及び燃料噴射無駄パルス幅設定ルーチンのステップＳ１〜Ｓ１０で１サイクル当りの燃料噴射質量を算出し、ステップＳ１１〜Ｓ２４で燃料噴射有効パルス幅Ｔｅと燃料噴射無駄パルス幅Ｔｓとを算出する。
【０１００】
まず、ステップＳ１では、上記筒内混合気の要求当量比設定ルーチンで設定した各要求当量比φｔｗ，φｆｕｌ，φｔｅｘ、φｇａｓ，φｅｃｏを比較し、これらの要求値の最大値を目標当量比ＣＯＥＦとして決定する。
【０１０１】
次いで、ステップＳ２で筒内への１吸気行程当りの筒内吸入燃料質量Ｇｆｕｅｌを次式から算出する。
【０１０２】
Ｇｆｕｅｌ←Ｇａｉｒ×Ｆ／Ａ×ＣＯＥＦ×Ｋｆｂ …（１５）
Ｆ／Ａ：理論燃空比
Ｋｆｂ：空燃比フィードバック補正係数
ここで、理論燃空比Ｆ／Ａは理論空燃比Ａ／Ｆの逆数で、理論空燃比が１４．６であれば、理論燃空比は１／１４．６となる。また、この理論空燃比Ａ／Ｆは、燃料を完全燃焼するに必要な最小の空気量と燃料量との比であり、多種燃料に対応したエンジン等では変数となる。
【０１０３】
一方、上記空燃比フィードバック補正係数Ｋｆｂの初期値は１．０であり、空燃比フィードバック制御、及び空燃比学習制御を行う外部のストラテジーによって書換えられる。また、理論燃空比Ｆ／Ａに目標当量比ＣＯＥＦを掛けることで、運転条件に応じた筒内混合気の目標空燃比が設定される。従って、リニアＡ／Ｆセンサでフィードバック制御する際に設定される上記空燃比フィードバック補正係数Ｋｆｂの目標値は上記理論燃空比Ｆ／Ａを１とした場合の上記目標当量比ＣＯＥＦである。
【０１０４】
ところで、上記筒内吸入燃料質量Ｇｆｕｅｌは、筒内へ１行程で吸込んで欲しい燃料量であり、定常的にはインジェクタ２５からの噴射量と一致するが、過渡時には応答遅れがあるため一致しない。すなわち、過渡的にはインジェクタ２５からの燃料噴射量が急増しても、燃料の一部が吸気ポート内壁に付着してから筒内に吸込まれるので、その分遅れが生じ、筒内に吸込まれる燃料量はゆっくりと増加する。従って、上記筒内吸入燃料質量Ｇｆｕｅｌが急増したときには、インジェクタ２５からの燃料量を、この筒内吸入燃料質量Ｇｆｕｅｌよりも多めにフィードフォワード噴射させて、筒内への吸入燃料量を筒内吸入燃料質量Ｇｆｕｅｌと一致させる必要がある。この追加増量分に相当する過渡時の１吸気行程当りの燃料追加質量Ｇａｃｃを、次のステップＳ３〜Ｓ９で算出する。
【０１０５】
まず、ステップＳ３では、スロットル開度αとエンジン回転数Ｎｅとに基づき二次元マップ（図１９参照）を補間計算付で参照して、静的な筒内吸入空気質量Ｇａｉｒに相当する指標値（静的指標値）Ｍａｃｃを設定する。この静的指標Ｍａｃｃは、負荷変動の少ない定常走行で、スロットル開度αとエンジン回転数Ｎｅを一定に保持すれば、最終的には筒内吸入空気質量Ｇａｉｒが一定になるということを前提にしている。従って、走行中の負荷変動が無視できるとすれば、図４のステップＳ３１に示すように、上記静的指標Ｍａｃｃをスロットル開度αに基づき一次元マップから補間計算付で設定することも可能である。
【０１０６】
また、過渡的には、上記筒内吸入空気質量Ｇａｉｒにも遅れがあり、この遅れに対応した燃料追加質量Ｇａｃｃを求める必要がある。
【０１０７】
そこで、ステップＳ４では、まず、吸気チャンバ６Ａ（図２４参照）内の吸気遅れを想定した所定時定数（運転条件によって異なるが本ルーチンでは４×１０ｍｓに固定している）の一次遅れ処理を、次式に示す加重平均に基づいて行い、動的な筒内吸入空気質量Ｇａｉｒに相当する指標値（動的指標値）Ｓａｃｃを算出する。
【０１０８】
Ｓａｃｃ←（３・Ｓａｃｃ＋Ｍａｃｃ）／４ …（１６）
次いで、ステップＳ５で、上記動的指標値Ｓａｃｃと、前回のルーチン実行時に算出した燃料付着による遅れ指標値Ｔａｃｃとを比較する。尚、この遅れ指標値Ｔａｃｃは上記筒内吸入燃料質量Ｇｆｕｅｌに相当する。
【０１０９】
そして、Ｓａｃｃ＜Ｔａｃｃ、すなわち、前回の筒内吸入燃料質量Ｇｆｕｅｌに相当する値（Ｔａｃｃ）に対して今回の筒内吸入空気質量Ｇａｉｒに相当する値（Ｓａｃｃ）が少ないときには、ステップＳ６へ進み、上記遅れ指標値Ｔａｃｃを上記動的指標値Ｓａｃｃとして、ステップＳ８へ進む。その結果、Ｓａｃｃ＝Ｔａｃｃとなり、後述のステップＳ９では、上記燃料追加質量Ｇａｃｃ＝０となり、燃料が追加増量されない。
【０１１０】
一方、Ｓａｃｃ≧Ｔａｃｃ、すなわち、前回の筒内吸入燃料質量Ｇｆｕｅｌに相当する値（Ｔａｃｃ）に対して今回の筒内吸入空気質量Ｇａｉｒに相当する値（Ｓａｃｃ）が増加しているか、或は変化していないときには、ステップＳ７へ分岐して、この動的指標Ｓａｃｃを、次式に示す加重平均により、燃料付着による時定数（本ルーチンでは、８×１０ｍｓ）の一次遅れ処理を行って、今回の燃料付着による遅れ指標値Ｔａｃｃを算出して、ステップＳ８へ進む。
【０１１１】
Ｔａｃｃ←（７・Ｔａｃｃ＋Ｓａｃｃ）／８ …（１７）
そして、ステップＳ８へ進むと、エンジン回転数Ｎｅに基づき一次元マップを補間計算付で参照して係数Ｒａｃｃを設定する。
【０１１２】
その後、ステップＳ９へ進み、過渡時の１吸気行程当りの燃料追加質量Ｇａｃｃを次式から算出する。
【０１１３】
Ｇａｃｃ←Ｒａｃｃ×（Ｓａｃｃ−Ｔａｃｃ） …（１８）
この（１８）式に示すように、今回の筒内吸入空気質量Ｇａｉｒに相当する値（Ｓａｃｃ）と今回の筒内吸入燃料質量Ｇｆｕｅｌに相当する値（Ｔａｃｃ）との差分が、燃料付着による不足分を補填する増量パターンということになる。
【０１１４】
このように、この燃料追加質量Ｇａｃｃは、スロットル開度αとエンジン回転数Ｎｅのみから、簡易的に筒内吸入空気質量Ｇａｉｒの変化量だけを遅れなく検出し、この値を利用して算出している。
【０１１５】
図２０にスロットル弁５ａを走行中に僅かに開弁させたときの上記各指標値Ｍａｃｃ，Ｓａｃｃ，Ｔａｃｃの特性を示す。静的指標値Ｍａｃｃはスロットル開度αに追従して増加し、一方、動的指標値Ｓａｃｃは上記静的指標値Ｍａｃｃに対して４×１０ｍｓの一次遅れで増加する。さらに、上記遅れ指標値Ｔａｃｃは、この動的指標値Ｓａｃｃに対して８×１０ｍｓの時定数で増加する。この遅れ指標値Ｔａｃｃが筒内吸入燃料質量Ｇｆｕｅｌに対応していると考えれば、付着分の遅れは、上記動的指標値Ｓａｃｃと上記遅れ指標値Ｔａｃｃとの差分（図のハッチングで示す領域）に相当することになり、上記ステップＳ９では、この差分指標値（Ｓａｃｃ−Ｔａｃｃ）に上記係数Ｒａｃｃを掛けて、燃料追加質量Ｇａｃｃを求めている。ところで、図２１に示すように、この差分指標値（Ｓａｃｃ−Ｔａｃｃ）が、従来のスロットル弁５ａの開弁直後に生じるリーンスパイクのパターンに符合することが、実験から明らかになった。従って、この差分指標値（Ｓａｃｃ−Ｔａｃｃ）相応分を加速増量すれば、排気空燃比は過渡時においてもリーン化することなく一定となる。
【０１１６】
上記燃料追加質量Ｇａｃｃは加速増量・加速時追加パルスに近似するが、この燃料追加質量Ｇａｃｃはあくまでも筒内への吸入燃料量を筒内吸入燃料質量Ｇｆｕｅｌにするためのフィードフォワード分であり、加速時に空燃比Ａ／Ｆを一時的にリッチにするためのものではない。従って、もし、加速時に空燃比Ａ／Ｆをリッチにする必要があれば、前記目標当量比ＣＯＥＦの選択枝として、前記図７〜図９に示す要求当量比設定ルーチン中に、加速当量比φａｃｃを第６の増量要因として加えることも、当然考えられる。
【０１１７】
また、先に述べた、従来の三角増量法では、筒内へ吸入される燃料量の一次遅れを三角形で近似させながら増量するので、上記図２１に示すリーンスパイクに対して正確に符合せず、過渡時の空燃比が部分的にリッチになったりリーンになったりしてしまう。
【０１１８】
そして、ステップＳ１０へ進むと、上記筒内吸入燃料質量Ｇｆｕｅｌに上記燃料追加質量Ｇａｃｃを加算して、１サイクル当りの燃料噴射質量Ｇｉｎｊを算出する。ところで、噴射時期が早期に設定されていて、噴射が終了した後でまだ吸気行程になる前のタイミングで上記燃料噴射質量Ｇｉｎｊが急増した場合には、自動的に追加噴射することも可能で、この燃料噴射質量Ｇｉｎｊは、ある吸気行程から１回転半前、すなわち吸気行程終了直後から次の吸気行程開始までの期間にインジェクタ２５から噴射される燃料の総量が算出される。
【０１１９】
次に、ステップＳ１１〜Ｓ２４で、気筒別燃料噴射有効パルス幅Ｔｅ１〜Ｔｅ４、及び気筒別燃料噴射無駄パルス幅Ｔｓ１〜Ｔｓ４が算出される（ここで、１，２，３，４は気筒番号を示す）。
【０１２０】
まず、ステップＳ１１で、最大出力当量比φｆｕｌの値を参照し、ステップＳ１２で排気ガス温度限界当量比φｔｅｘの値を参照する。そして、φｆｕｌ＝０のエンジンが最大出力を要求しておらず、しかも、φｔｅｘ＝０の排気ガス温度限界による増量も要求していない場合には、ステップＳ１３へ進み、空燃比の気筒別のばらつきを補正する気筒別噴射量補正係数Ｋｔｎ１〜Ｋｔｎ４を１．０に設定して、ステップＳ１８へ進む。一方、上記ステップＳ１１で、φｆｕｌ＝１．２（フル増量）、あるいは、ステップＳ１２で、φｔｅｘ≠０と判断されている場合には、ステップＳ１４へ分岐し、ステップＳ１４〜Ｓ１７でエンジン回転数Ｎｅに基づき、気筒別の一次元マップを補間計算付で参照して、上記気筒別噴射量補正係数Ｋｔｎ１〜Ｋｔｎ４を各々設定して、ステップＳ１８へ進む。噴射量補正係数Ｋｔｎは燃料噴射量を気筒別に増減量して空燃比を均一に設定するもので、上記各一次元マップは、各気筒の吸気特性に対応する値が格納されている。
【０１２１】
ところで、運転条件によって燃料噴射量を気筒毎に変えなければならない原因は、本来、筒内吸入空気質量Ｇａｉｒが気筒毎に相違するからであり、従って、正確には筒内吸入空気質量Ｇａｉｒを吸気特性を考慮して気筒毎に算出し、この算出結果から筒内吸入燃料質量Ｇｆｕｅｌ及び燃料噴射質量Ｇｉｎｊを気筒毎に算出すべきであるが、ＣＰＵの計算負荷が重くなるため、上記気筒別噴射量補正係数Ｋｔｎ１〜Ｋｔｎ４を用いて、見かけ上、気筒別燃料噴射有効パルス幅Ｔｅ１〜Ｔｅ４［ｍｓ］を補正するようにした。
【０１２２】
また、空燃比の気筒毎のばらつきが特に問題となるのは、高負荷運転時のノッキングである。中低負荷運転での空燃比のばらつきは排気ガスに余り影響せず、エンジン振動の観点からすると全気筒の燃料量が一致している方が却って各気筒の出力が均等になるので望ましい。従って、フル増量が要求されている場合には、全気筒が最大出力が得られるようにする必要があり、また、排気ガス温度限界による増量が要求されている場合には、全気筒の排気系等の高温化を設計限界以内に抑制してエンジンを保護する必要があるため、このような場合、すなわち、気筒間の空燃比が均一であることを前提とする要求増量が設定されている場合のみ、エンジン回転数毎に各気筒の燃料噴射量を適正に増減させる補正係数（気筒別噴射量補正係数）を設定し、各気筒の筒内混合気の空燃比が目標空燃比となるように個別に制御する。
【０１２３】
そして、ステップＳ１８へ進むと、全体としての燃料噴射有効パルス幅Ｔｅ＿ａｌｌを次式から算出する。
【０１２４】
Ｔｅ＿ａｌｌ←Ｋｍｒ×Ｇｉｎｊ×Ｋｃｏｎ …（１９）
Ｋｍｒ：パルス幅誤差補正係数
Ｋｃｏｎ：インジェクタ容量係数［ｓｅｃ／ｇ］
ここで、パルス幅誤差補正係数Ｋｍｒは、インジェクタ２５の動的な流量特性の非線型性の補正であり、本ルーチンでは、エンジン回転数Ｎｅと燃料噴射質量Ｇｉｎｊとに基づき、図２２に示す二次元マップを補間計算付で参照して設定する。また、上記インジェクタ容量係数Ｋｃｏｎは、インジェクタ２５の静的質量流量特性の逆数であり、１グラムの燃料を０．１秒間に噴射するインジェクタでは、Ｋｃｏｎ＝０．１［ｓｅｃ／ｇ］である。なお、このインジェクタ容量係数Ｋｃｏｎは、本ルーチンでは固定値としているが、多種燃料に対応可能なＦＦＶエンジン等で、燃料の比重や粘性が変化する場合には変数として用いる。
【０１２５】
次いで、ステップＳ１９〜Ｓ２２で、上記燃料噴射有効パルス幅Ｔｅ＿ａｌｌを上記気筒別噴射量補正係数Ｋｔｎ１〜Ｋｔｎ４で補正して、気筒別燃料噴射有効パルス幅Ｔｅ１〜Ｔｅ４をそれぞれ設定する。
【０１２６】
その後、ステップＳ２３へ進み、バッテリ電圧ＶＢに基づき、一次元マップを補間計算付で参照し燃料噴射無駄パルス幅Ｔｓ＿ａｌｌ［ｍｓ］を設定し、ステップＳ２４で、気筒別燃料噴射無駄パルス幅Ｔｓ１〜Ｔｓ４を、上記燃料噴射無駄パルス幅Ｔｓ＿ａｌｌで設定してルーチンを抜ける。
【０１２７】
ＯＳ側では、上記気筒別燃料噴射有効パルス幅Ｔｅ１〜Ｔｅ４と上記気筒別燃料噴射無駄パルス幅Ｔｓ１〜Ｔｓ４とをそれぞれ加算して、燃料噴射パルス幅を気筒別に算出する。
【０１２８】
インジェクタ２５に与える電圧パルスの幅（燃料噴射パルス幅）と、燃料噴射質量Ｇｉｎｊとの関係は、図２３に示す通りであり、燃料噴射質量Ｇｉｎｊが決定されれば、基本的に、この燃料噴射質量Ｇｉｎｊにインジェクタ容量係数Ｋｃｏｎを掛けることで上記有効燃料噴射パルス幅Ｔｅが求められ、この有効燃料噴射パルス幅Ｔｅに無効噴射パルス幅である燃料噴射無駄パルス幅Ｔｓを加算することで燃料噴射パルス幅が算出される。
【０１２９】
本実施例では、燃料噴射質量Ｇｉｎｊ分の燃料を１回で噴射するか、２回に分割して噴射するかはＯＳ内部で決定されるため、上記気筒別燃料噴射有効パルス幅Ｔｅ１〜Ｔｅ４と上記気筒別燃料噴射無駄パルス幅Ｔｓ１〜Ｔｓ４とを、予め加算せずに別々にＯＳに渡し、このＯＳ内部で噴射設定の直前に加算して、各インジェクタに与える燃料噴射パルス幅を算出している。こうすることで、ＯＳ内部で、燃料噴射質量Ｇｉｎｊ分の燃料を２回に分割して噴射すると決定した場合でも、１／２の燃料噴射有効パルス幅Ｔｅに燃料噴射無駄パルス幅Ｔｓを加算することで、１回当りの燃料噴射パルス幅を簡単に与えることができる。
【０１３０】
このように、本実施例によれば、例えば新しい燃料付着モデルを開発した際は、筒内吸入空気質量Ｇａｉｒを算出する前記（３）式、或は（１０）式中に演算式を追加すれば良く、またインジェクタの流量が非線型であるならば燃料噴射有効パルス幅Ｔｅ，燃料噴射無駄パルス幅Ｔｓを補正するだけで良いことになる。その結果、ある運転領域での空燃比制御性に不具合が生じた場合に、どの部分を変更すべきかが明確となるばかりか、その変更が他の部分に影響を及ぼし難いため、ベース制御ソフトとしての活用が期待できる。
【０１３１】
【発明の効果】
以上、説明したように本発明によれば、スロットル開度とエンジン回転数とに基づき静的な筒内吸入空気質量に相当する静的指標値を設定し、静的指標値に吸気チャンバによる遅れを想定した一次遅れ処理を施して動的な筒内吸入空気質量に相当する動的指標値を設定する。そして、更に動的指標値を加重平均して燃料付着による遅れ指標値を算出し、上記動的指標値から遅れ指標値を減算して、この差分に基づき燃料追加質量を設定する。そして、基本とする吸入燃料質量に、動的指標値と遅れ指標値との差分に基づいて設定した燃料追加質量を加算して１サイクル当りの燃料噴射量を設定するので、従来の燃料付着モデルを用いた場合に比し、応答遅れや、ノイズの影響を受けることが無く、過渡時の燃料追加質量を高精度に求めることができて、信頼性の高い空燃比制御性能を得ることができる。
特に、燃料付着補正については、スロットル開度とエンジン回転数とに基づいて設定した静的な筒内吸入空気質量に相当する静的指標値に対し、一次遅れ処理を施して、更に加重平均処理を施し、一次遅れ処理により設定した動的な筒内吸入空気質量に相当する動的指標値から、該動的指標値を加重平均処理して設定した遅れ指標値を減算し、この差分に基づき燃料追加質量を設定するので、マップやテーブル、或いは、マップやテーブルに相当する補正回路等を必要とせず、スロットル開度とエンジン回転数のみから、非常に簡素な処理で、応答遅れを解消して、信頼性の高い高精度な燃料付着補正を実現できる。
【図面の簡単な説明】
【図１】１０ｍｓ毎ジョブで実行される燃料噴射有効パルス幅及び燃料噴射無駄パルス幅設定ルーチンを示すフローチャート
【図２】１０ｍｓ毎ジョブで実行される燃料噴射有効パルス幅及び燃料噴射無駄パルス幅設定ルーチンを示すフローチャート（続き）
【図３】１０ｍｓ毎ジョブで実行される燃料噴射有効パルス幅及び燃料噴射無駄パルス幅設定ルーチンを示すフローチャート（続き）
【図４】他の態様による燃料噴射有効パルス幅及び燃料噴射無駄パルス幅設定ルーチンの要部を示すフローチャート
【図５】１０ｍｓ毎ジョブで実行される筒内吸入空気質量設定ルーチンを示すフローチャート
【図６】５０ｍｓ毎ジョブで実行される係数設定ルーチンを示すフローチャート
【図７】５０ｍｓ毎ジョブで実行される要求当量比設定ルーチンを示すフローチャート
【図８】５０ｍｓ毎ジョブで実行される要求当量比設定ルーチンを示すフローチャート（続き）
【図９】５０ｍｓ毎ジョブで実行される要求当量比設定ルーチンを示すフローチャート（続き）
【図１０】他の態様による４ｍｓ毎ジョブ実行される筒内吸入空気質量設定ルーチンを示すフローチャート
【図１１】筒内吸入空気質量と吸気温度との関係を示す特性図
【図１２】理論吸気質量と筒内吸入空気質量との関係を示す特性図
【図１３】吸入空気量誤差補正係数を設定する二次元マップの概念図
【図１４】過渡時における吸入空気量センサの出力電圧を示す特性図
【図１５】基本可燃限界当量比を設定する二次元マップの概念図
【図１６】最大出力当量比の増量領域を示す特性図
【図１７】基本排気ガス温度限界当量比を設定する二次元マップの特性図
【図１８】燃費率最良当量比を設定する二次元マップの特性図
【図１９】静的な筒内吸入空気質量に相当する指標値を設定する二次元マップの概念図
【図２０】静的指標値と動的指標値とこの動的指標値に対する遅れ指標値との関係を示す特性図
【図２１】過渡運転時のリーンスパイクと増量パターンとの関係を示す特性図
【図２２】パルス幅誤差補正係数を設定する二次元マップの特性図
【図２３】燃料噴射パルス幅と燃料噴射質量との関係を示す特性図
【図２４】エンジンの吸気系モデルを示す説明図
【図２５】エンジンの全体構成図
【図２６】クランクロータとクランク角センサの正面図
【図２７】カムロータとカム角センサの正面図
【図２８】電子制御装置の回路構成図
【符号の説明】
６Ａ…吸気チャンバ
Ｇａｉｒ…吸入空気質量
Ｇａｃｃ…燃料追加質量
Ｇｆｕｅｌ…吸入燃料質量
Ｇｉｎｊ…燃料噴射質量
Ｍａｃｃ…静的指標値
Ｓａｃｃ…動的指標値
Ｔａｃｃ…遅れ指標値
α…スロットル開度
Claims (1)
- 筒内への１吸気行程当りの吸入空気質量を基本として目標空燃比に対応する筒内への１吸気行程当りの吸入燃料質量を設定するエンジンの空燃比制御方法において、
スロットル開度とエンジン回転数とに基づき静的な筒内吸入空気質量に相当する静的指標値を設定し、
上記静的指標値に吸気チャンバによる遅れを想定した一次遅れ処理を施して動的な筒内吸入空気質量に相当する動的指標値を設定し、
更に上記動的指標値を加重平均して燃料付着による遅れ指標値を算出し、
上記動的指標値から上記遅れ指標値を減算して、この差分に基づいて燃料追加質量を設定し、
上記燃料追加質量を上記吸入燃料質量に加算して１サイクル当りの燃料噴射量を設定することを特徴とするエンジンの空燃比制御方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP28097194A JP3621731B2 (ja) | 1994-11-15 | 1994-11-15 | エンジンの空燃比制御方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP28097194A JP3621731B2 (ja) | 1994-11-15 | 1994-11-15 | エンジンの空燃比制御方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH08144805A JPH08144805A (ja) | 1996-06-04 |
| JP3621731B2 true JP3621731B2 (ja) | 2005-02-16 |
Family
ID=17632454
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP28097194A Expired - Fee Related JP3621731B2 (ja) | 1994-11-15 | 1994-11-15 | エンジンの空燃比制御方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3621731B2 (ja) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2965584B1 (fr) * | 2010-10-05 | 2013-06-28 | Renault Sas | Procede de determination d'un taux de gaz d'echappement recircules a l'entree d'un cylindre d'un moteur a combustion interne et moteur mettant en oeuvre un tel procede |
-
1994
- 1994-11-15 JP JP28097194A patent/JP3621731B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JPH08144805A (ja) | 1996-06-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP2884472B2 (ja) | 内燃機関の燃料性状検出装置 | |
| JP3314294B2 (ja) | 内燃機関の制御装置 | |
| JPH0953487A (ja) | 内燃機関の燃料噴射制御装置 | |
| JP2927074B2 (ja) | 内燃機関の空燃比制御装置 | |
| JP3323542B2 (ja) | 吸気温度予測方法 | |
| JP3621731B2 (ja) | エンジンの空燃比制御方法 | |
| JP3691092B2 (ja) | エンジンの空燃比制御方法 | |
| JP3686437B2 (ja) | エンジンの空燃比制御方法 | |
| JP3593162B2 (ja) | エンジンの吸入空気量検出方法 | |
| JP3963993B2 (ja) | エンジンの空燃比制御装置 | |
| JP2916804B2 (ja) | 内燃機関の空燃比制御装置 | |
| JPH08144834A (ja) | エンジンの吸入空気量検出方法 | |
| JP4263258B2 (ja) | リーンバーンエンジンの燃料噴射制御装置 | |
| JPH08144835A (ja) | エンジンの吸入空気量検出方法 | |
| JPS61157741A (ja) | 吸入空気量検出装置 | |
| JP2543763B2 (ja) | 内燃機関の燃料供給制御装置 | |
| JP3817830B2 (ja) | 内燃機関の点火時期制御装置 | |
| JPS63268951A (ja) | 内燃機関の燃料供給制御装置 | |
| JPH0710048Y2 (ja) | 内燃機関の燃料供給制御装置 | |
| JP2500946Y2 (ja) | 内燃機関の電子制御燃料供給装置 | |
| JPH0656112B2 (ja) | 内燃機関の燃料噴射制御装置 | |
| JPH068298Y2 (ja) | 内燃機関の点火時期制御装置 | |
| JPH1193782A (ja) | リーンバーンエンジンの蒸発燃料パージ制御装置 | |
| JPS61201842A (ja) | 内燃機関のアイドル時稀薄空燃比制御装置 | |
| JPH06101564A (ja) | 内燃機関の出力変動検出装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20040812 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20040824 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20041004 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20041116 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20041119 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| LAPS | Cancellation because of no payment of annual fees |