JP3593162B2 - エンジンの吸入空気量検出方法 - Google Patents
エンジンの吸入空気量検出方法 Download PDFInfo
- Publication number
- JP3593162B2 JP3593162B2 JP28097394A JP28097394A JP3593162B2 JP 3593162 B2 JP3593162 B2 JP 3593162B2 JP 28097394 A JP28097394 A JP 28097394A JP 28097394 A JP28097394 A JP 28097394A JP 3593162 B2 JP3593162 B2 JP 3593162B2
- Authority
- JP
- Japan
- Prior art keywords
- intake
- intake air
- cylinder
- fuel
- amount
- 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
- Combined Controls Of Internal Combustion Engines (AREA)
- Output Control And Ontrol Of Special Type Engine (AREA)
- Electrical Control Of Air Or Fuel Supplied To Internal-Combustion Engine (AREA)
Description
【０００１】
【産業上の利用分野】
本発明は、Ｄジェトロニック方式により吸入空気量を算出するエンジンに関し、詳しくはスロットル弁下流の空気密度から筒内吸入空気量を設定するエンジンの吸入空気量検出方法に関する。
【０００２】
【従来の技術】
一般に、この種の燃料噴射制御では、吸入空気量を吸入空気量センサの計測値に基づいて算出するＬジェトロニック方式と、スロットル弁下流の吸気管圧力に基づいて算出するＤジェトロニック方式とがあり、その中でも、最近のＤジェトロニック方式を採用するエンジンでは、筒内へ実際に吸入される空気量を過渡時であっても正確に推定することのできる技術が種々提案されている。
【０００３】
例えば、特開平６−１８５３９１号公報では、筒内への吸入空気量を空気の状態方程式を用いて推定している。すなわち、スロットル弁下流のチャンバ容積に対する吸入空気の入出力関係に着目し、まず、スロットル弁を単位時間（Δｔ）当りに通過する吸気量（スロットル通過空気量）Ｇｔｒ・Δｔを設定し、次いで、スロットル弁下流のチャンバ内空気量Ｇｂの圧力変化分に基づいて今回チャンバに充填された空気量の変化分ΔＧｂを設定する。そして、今回チャンバに充填された空気量分は、筒内へ吸入されなかったものと見なして、単位時間当りの筒内吸入空気量Ｇｃを算出している。これを式に表せば以下の通である。
【０００４】
Ｇｃ＝Ｇｔｒ×Δｔ−ΔＧｂ
尚、上記スロットル通過空気量Ｇｔｒは、スロットル弁の上下流の圧力差、すなわち大気圧と上記吸気管圧力との偏差、スロットル有効開口面積、及びスロットル弁上流側の空気密度等に基づき近似的に算出している。
【０００５】
【発明が解決しようとする課題】
ところで、上記先行技術によるＤジェトロニック方式を用いた筒内吸入空気量の算出過程では、以下の理由からコンピュータに掛かる負荷が重くなる。
【０００６】
（１）スロットル下流側チャンバへの流入吸気量と、このチャンバに充填された吸入空気量とを個々に求めるなど、複雑な演算を必要としている。
【０００７】
（２）例えばスロットル通過空気量Ｇｔｒを設定する関数として必要なスロットル有効開口面積Ａは、係数Ｃとスロットル投影面積Ｓとの積で設定しているが、この係数Ｃを設定する際に、スロットル開度θＴＨと吸気管圧力Ｐとに基づいて二次元マップを参照するようにしている。
【０００８】
一方、上記先行技術では、チャンバ内空気量Ｇｂ及びこのチャンバ内空気量の変化分ΔＧｂは、吸気管圧力Ｐの変化分ΔＰが一定であれば、ほぼ空気温度Ｔに反比例した値に設定されるが、実測結果によれば、ターボエンジンなどでは、チャンバ内の吸気温度がある温度を越えると、吸入空気量は吸気温度に依存して減少せず、従って、全運転領域の吸入空気量を吸気温度に反比例した特性で設定すると、高回転、高負荷運転領域などでは、計算により得られた吸入空気量が実際の吸入空気量に対して過小になり、この過小な吸入空気量に基づいて燃料噴射量が設定されるため、空燃比がリーン化してしまう課題がある。
【０００９】
本発明は、上記事情に鑑みてなされたもので、一次式及び一次元マップにより筒内への吸入空気量を簡単に算出することができて、コンピュータに対する負担を軽減することができ、しかも、スロットル弁下流の吸気温度が高温になっても筒内への吸入空気量を実際の吸入空気量とほぼ等しい値に設定することができて、空燃比のリーン化を防止するとともに、過渡応答性の良いエンジンの吸入空気量検出方法を提供することを目的としている。
【００１０】
【課題を解決するための手段】
上記目的を達成するため第１発明によるエンジンの吸入空気量検出方法は、少なくともスロットル弁下流の吸気管圧力と吸気温度とに基づき空気の状態方程式から空気密度を設定し、この空気密度と行程容積とから１吸気行程当りの筒内への理論吸気量を設定し、またエンジン回転数に基づき１吸気行程当りの筒内への吸気損失量と体積効率とを各々に対応する一次元マップを参照して設定し、上記理論吸気量と上記吸気損失量との差分と上記体積効率とから１吸気行程当りの筒内への吸入空気量を設定することを特徴とする。
【００１１】
上記目的を達成するため第２発明によるエンジンの吸入空気量検出方法は、前記吸気温度が、前記筒内への吸入空気量が吸気温度に依存して変動しない温度として予め設定した臨界温度以上のときには、上記吸気温度を固定値で設定して上記吸気密度を設定することを特徴とする。
【００１２】
【作 用】
第１発明では、筒内への吸入空気量を算出するに際して、まず、吸気行程時の筒内の空気密度を、少なくともスロットル弁下流の吸気管圧力と吸気温度とに基づき空気の状態方程式から算出する。次いでこの空気密度と行程容積との積から、１吸気行程当りの筒内への理論吸気量を算出する。一方、エンジン回転数に基づいて１吸気行程当りの吸気損失量と体積効率を各々に対応する一次元マップを参照して設定し、上記理論吸気量と上記吸気損失量との差と、上記体積効率との積から、１吸気行程当りの筒内吸入空気量を算出する。
【００１３】
第２発明では、筒内吸入空気量が吸気温度に依存して変動しなくなる温度を臨界温度として予め設定し、上記吸気温度が臨界温度を越えたときには、実際の吸気温度に拘らず、上記固定値と上記吸気管圧力とに基づいて吸気密度を設定する。
【００１４】
【実施例】
以下、図面に基づき本発明の実施例を説明する。図２３にはエンジンの全体概略図が示されており、図中の符号１はエンジン（図においては水平対向４気筒型エンジンを示す）を示し、シリンダヘッド２の吸気ポート２ａにインテークマニホルド３が連通され、このインテークマニホルド３の上流にエアーチャンバ４を介してスロットル通路５が連通されている。このスロットル通路５の上流側には、吸気管６を介してエアクリーナ７が取付けられ、このエアクリーナ７が吸入空気の取り入れ口であるエアインテークチャンバ８に連通されている。また、上記排気ポート２ｂにエキゾーストマニホルド９を介して排気管１０が連通され、この排気管１０に触媒コンバータ１１が介装されてマフラ１２に連通されている。
【００１５】
一方、上記スロットル通路５にスロットル弁５ａが設けられ、このスロットル通路５の直上流の上記吸気管６にインタークーラ１３が介装され、さらに、上記吸気管６の上記エアクリーナ７の下流側にレゾネータチャンバ１４が介装されている。
【００１６】
また、上記レゾネータチャンバ１４と上記インテークマニホルド３とを連通して上記スロットル弁５ａの上流側と下流側とをバイパスするバイパス通路１５に、アイドル空気量を調整するアイドル制御（ＩＳＣ）弁１６が介装されている。さらに、上記ＩＳＣ弁１６の直下流側に、吸気圧が負圧のとき開弁し、またターボチャージャ１８によって過給されて吸気圧が正圧になったとき閉弁するチェックバルブ１７が介装されている。
【００１７】
上記ターボチャージャ１８は、上記吸気管６の上記レゾネータチャンバ１４の下流側にコンプレッサが介装され、タービンが上記排気管１０に介装されている。さらに、上記ターボチャージャ１８のタービンハウジング流入口には、ウエストゲート弁１９が介装され、このウエストゲート弁１９には、ウエストゲート弁作動用アクチュエータ２０が連設されている。
【００１８】
上記ウエストゲート弁作動用アクチュエータ２０は、ダイヤフラムにより２室に仕切られ、一方がウエストゲート弁制御用デューティソレノイド弁２１に連通される圧力室を形成し、他方が上記ウエストゲート弁１９を閉方向に付勢するスプリングを収納したスプリング室を形成している。
【００１９】
上記ウエストゲート弁制御用デューティソレノイド弁２１は、上記レゾネータチャンバ１４と上記吸気管６の上記ターボチャージャ１８のコンプレッサ下流とを連通する通路に介装されており、後述する電子制御装置５０（ＥＣＵ；図２６参照）から出力される制御信号のデューティ比に応じて、上記レゾネータチャンバ１４側の圧力と上記コンプレッサ下流側の圧力とを調圧し、上記ウエストゲート弁作動用アクチュエータ２０の圧力室に供給する。
【００２０】
すなわち、上記ＥＣＵ５０によって上記ウエストゲート弁制御用デューティソレノイド弁２１を制御し、上記ウエストゲート弁作動用アクチュエータ２０を作動させて上記ウエストゲート弁１９による排気ガスリリーフを調整することにより、上記ターボチャージャ１８による過給圧を制御するようになっている。
【００２１】
また、上記インテークマニホルド３に、吸気管圧力Ｐを絶対圧で検出する吸気管圧力センサ２２が通路２３を介して連通され、さらに、上記インテークマニホルド３の各気筒の各吸気ポート２ａの直上流側にインジェクタ２５が臨まされている。また、上記シリンダヘッド２の各気筒毎に、その先端を燃焼室に露呈する点火プラグ２６ａが取付けられ、この点火プラグ２６ａに各気筒毎に配設された点火コイル２６ｂを介してイグナイタ２７が接続されている。
【００２２】
上記インジェクタ２５には、燃料タンク２８内に設けたインタンク式の燃料ポンプ２９から燃料フィルタ３０を経て燃料が圧送され、プレッシャレギュレータ３１にてインジェクタ２５への燃料圧力が調圧される。
【００２３】
また、上記スロットル弁５ａに、スロットル開度センサ３３ａとアイドルスイッチ３３ｂとを内蔵したスロットルセンサ３３が連設され、上記エアーチャンバ４に吸気温センサ４６が臨まされている。
【００２４】
また、上記エンジン１のシリンダブロック１ａにノックセンサ３４が取付けられるとともに、このシリンダブロック１ａの左右両バンクを連通する冷却水通路３５に水温センサ３６が臨まされ、上記排気管１０の上記エキゾーストマニホルド９の集合部にＯ２ センサ３７が臨まされている。
【００２５】
また、上記シリンダブロック１ａに支承されたクランクシャフト１ｂにクランクロータ３８が軸着され、このクランクロータ３８の外周に、電磁ピックアップなどからなるクランク角センサ３９が対設されている。さらに、上記エンジン１のカムシャフト１ｃに連設するカムロータ４０に、電磁ピックアップなどからなる気筒判別用のカム角センサ４１が対設されている。尚、上記クランク角センサ３９及び上記カム角センサ４１は、電磁ピックアップなどの磁気センサに限らず、光センサなどでも良い。
【００２６】
上記クランクロータ３８は、図２４に示すように、その外周に突起３８ａ，３８ｂ，３８ｃが形成され、これらの各突起３８ａ，３８ｂ，３８ｃが、各気筒（＃１，＃２と＃３，＃４）の圧縮上死点前（ＢＴＤＣ）θ１，θ２，θ３ の位置に形成されており、本実施例においては、θ１ ＝９７゜ＣＡ、θ２ ＝６５゜ＣＡ、θ３ ＝１０゜ＣＡである。
【００２７】
上記クランクロータ３８の各突起は、上記クランク角センサ３９によって検出され、ＢＴＤＣ９７゜，６５゜，１０゜のクランクパルスがエンジン１／２回転毎（１８０゜ＣＡ毎）に出力される。そして、各信号の入力間隔時間がタイマによって計時され、エンジン回転数Ｎｅが算出される。
【００２８】
また、図２５に示すように、上記カムロータ４０の外周には、気筒判別用の突起４０ａ，４０ｂ，４０ｃが形成され、突起４０ａが＃３，＃４気筒の圧縮上死点後（ＡＴＤＣ）θ４ の位置に形成され、突起４０ｂが３個の突起で構成されて最初の突起が＃１気筒のＡＴＤＣθ５ の位置に形成されている。さらに、突起４０ｃが２個の突起で形成され、最初の突起が＃２気筒のＡＴＤＣθ６ の位置に形成されている。本実施例においては、θ４ ＝２０゜ＣＡ、θ５ ＝５゜ＣＡ、θ６ ＝２０゜ＣＡである。
【００２９】
そして、上記カムロータ４０の各突起が上記カム角センサ４１によって検出され、各気筒の燃焼行程順を＃１→＃３→＃２→＃４とした場合、この燃焼行程順と、上記カム角センサ４１からのカムパルスをカウンタによって計数した値とのパターンに基づいて、気筒判別がなされる。
【００３０】
一方、図２６において、符号５０は、エンジン系を制御する電子制御装置（ＥＣＵ）５０であり、このＥＣＵ５０は、燃料噴射制御、点火時期制御などを行なうメインコンピュータ５１と、ノック検出処理を行なう専用のサブコンピュータ５２との２つのコンピュータを中心として構成され、各部に所定の安定化電源を供給する定電圧回路５３や各種の周辺回路が組込まれている。
【００３１】
上記定電圧回路５３は、ＥＣＵリレー５４のリレー接点を介してバッテリ５５に接続されており、このバッテリ５５に、上記ＥＣＵリレー５４のリレーコイルがイグニッションスイッチ５６を介して接続され、上記イグニッションスイッチ５６がＯＮされて上記ＥＣＵリレー５４のリレー接点が閉となったときには制御用電源を上記各コンピュータ５１，５２に供給すると共に、上記定電圧回路５３は上記バッテリ５５に直接接続されており、イグニッションスイッチ５６のＯＮ，ＯＦＦに拘らずバックアップＲＡＭ６１にバックアップ用の電源を供給する。
【００３２】
尚、上記バッテリ５５には、燃料ポンプリレー５７のリレー接点を介して燃料ポンプ２９が接続されている。
【００３３】
上記メインコンピュータ５１は、ＣＰＵ５８、ＲＯＭ５９、ＲＡＭ６０、バックアップＲＡＭ６１、カウンタ・タイマ群６２、シリアル通信インターフェースであるＳＣＩ６３、及び、Ｉ／Ｏインターフェース６４がバスライン６５を介して接続されたマイクロコンピュータであり、上記バックアップＲＡＭ６１には、上記イグニッションスイッチ５６のＯＮ／ＯＦＦに拘らず、バッテリ５５に直接接続する上記定電圧回路５３からバックアップ電源が常時供給されてデータが保持される。
【００３４】
尚、上記カウンタ・タイマ群６２は、フリーランカウンタ、カム角センサ信号の入力計数用カウンタなどの各種カウンタ、燃料噴射タイマ、点火タイマ、定期割込みを発生させるための定期割込みタイマ、クランク角センサ信号の入力間隔計時用タイマ、及び、システム異常監視用のウオッチドッグタイマなどの各種タイマを便宜上総称するものであり、上記メインコンピュータ５１においては、その他、各種のソフトウエアカウンタ・タイマが用いられる。
【００３５】
また、上記サブコンピュータ５２も、上記メインコンピュータ５１と同様、ＣＰＵ７１、ＲＯＭ７２、ＲＡＭ７３、カウンタ・タイマ群７４、ＳＣＩ７５、及び、Ｉ／Ｏインターフェース７６がバスライン７７を介して接続されたマイクロコンピュータであり、上記メインコンピュータ５１とサブコンピュータ５２とは、上記ＳＣＩ６３，７５を介してシリアル通信ラインにより互いに接続されている。
【００３６】
上記メインコンピュータ５１のＩ／Ｏインターフェース６４には、入力ポートに、スロットル開度センサ３３ａ、水温センサ３６、Ｏ２センサ３７、吸気管圧力センサ２２、大気圧センサ４４、車速センサ４２、吸気温センサ４６、及び、バッテリ５５が、Ａ／Ｄ変換器６６を介して接続されるとともに、アイドルスイッチ３３ｂ、スタータスイッチ４３、及びクランク角センサ３９、カム角センサ４１などが接続され、さらに、図示しない各種センサ及びスイッチ類が接続されている。
【００３７】
また、上記Ｉ／Ｏインターフェース６４の出力ポートには、イグナイタ２７が接続されるとともに、ＩＳＣ弁１６、インジェクタ２５、燃料ポンプリレー５７のリレーコイル、ウエストゲート弁制御用デューティソレノイド弁２１が駆動回路６７を介して接続されており、さらに、図示しない各種のアクチュエータ類が接続されている。
【００３８】
一方、上記サブコンピュータ５２のＩ／Ｏインターフェース７６は、入力ポートに、クランク角センサ３９、カム角センサ４１が接続されるとともに、Ａ／Ｄ変換器７８、周波数フィルタ７９、アンプ８０を介してノックセンサ３４が接続されており、上記ノックセンサ３４からのノック検出信号が上記アンプ８０で所定のレベルに増幅された後に上記周波数フィルタ７９により必要な周波数成分が抽出され、上記Ａ／Ｄ変換器７８にてデジタル信号に変換されて入力されるようになっている。
【００３９】
上記メインコンピュータ５１では、各センサ類からの検出信号を処理し、燃料噴射量制御、点火時期制御、アイドル制御などを行い、一方、上記サブコンピュータ５２では、エンジン回転数とエンジン負荷とに基づいてノックセンサ３４からの信号のサンプル区間を設定し、このサンプル区間でノックセンサ３４からの信号を高速にＡ／Ｄ変換して振動波形を忠実にデジタルデータに変換し、このデータに基づきノック発生の有無を判定する。
【００４０】
上記サブコンピュータ５２のＩ／Ｏインターフェース７６の出力ポートは、上記メインコンピュータ５１のＩ／Ｏインターフェース６４の入力ポートに接続されており、上記サブコンピュータ５２でのノック判定結果がＩ／Ｏインターフェース７６に出力される。そして、上記メインコンピュータ５１では、上記サブコンピュータ５２からノック発生有りの判定結果が出力されると、ＳＣＩ６３を介してシリアル通信ラインよりサブコンピュータ５２からノックデータを読込み、このノックデータに基づいて直ちに該当気筒の点火時期を遅らせ、ノックを回避する。
【００４１】
このようなエンジン制御において、上記メインコンピュータ５１では、センサ・スイッチ類からの信号入力処理、燃料噴射制御、点火時期制御、アイドル制御に係わる各種ジョブが、一つのオペレーティングシステム（ＯＳ）の管理下で効率的に実行される。このＯＳは、車輌制御のための各種マネジメント機能、及び、このマネジメント機能に密着した内部ストラテジーを有し、各種ジョブを体系的に結合し、等時間間隔処理により各種ジョブを効率的に実行するようになっている。
【００４２】
以下、上記メインコンピュータ５１による燃料噴射制御について、図１〜図９のルーチンに従って説明する。尚、サブコンピュータ５２はノック検出処理専用のコンピュータであるため、その動作説明を省略する。
【００４３】
本実施例では、ＯＳ側で処理したＡ／Ｄ変換結果、クランク位置情報、エンジン回転数等に基づき、ユーザ側ジョブで、燃料噴射量、点火時期等を設定すると共に、これらの指示値がＯＳによって噴射タイマ、点火タイマにセットされる。そして、１０ｍｓ毎ジョブ，５０ｍｓ毎ジョブ等の定期割込み要求を出力し、またクランク角センサ３９からの信号入力毎（ＢＴＤＣ９７゜，６５゜，１０゜ＣＡ毎のエンジン１回転に６回）に起動されるクランクパルス入力毎ジョブの割込みを許可する。
【００４４】
１０ｍｓ毎ジョブでは、図１〜図３に示す燃料噴射有効パルス幅及び燃料噴射無駄パルス幅設定ルーチンと、図５に示す吸入空気量設定ルーチンが実行される。また、５０ｍｓ毎ジョブでは、図６に示す係数設定ルーチン、及び図７〜図９に示す要求当量比設定ルーチンが実行される。
【００４５】
尚、以下の説明では、図２２に示すように、吸気系のスロットル弁５ａの下流から吸気弁上流までを吸気チャンバ６Ａと総称する。従って、吸気チャンバ６Ａは、スロットル通路５、エアーチャンバ４、インテークマニホルド３、吸気ポート２ａの総称となる。
【００４６】
ここで、以下の説明で使用する用語について簡単に説明する。Ｑは質量流量［ｍｇ／ｓｅｃ］、Ｇは１サイクル当りの質量［ｍｇ／ｃｙｃｌｅ］を意味する。ただし、後述する筒内吸入空気質量Ｇａｉｒ、筒内吸入燃料質量Ｇｆｕｅｌ等は、吸気行程時にのみ吸入されるので、実質的には、１吸気行程当りの質量と同じ意味になる。また、Ｍは吸気チャンバ６Ａの空気質量、Ｐは吸気チャンバ６Ａの吸気管内圧力、Ｇｉｎｊは燃料噴射質量、Ｑｃは筒内吸入空気質量流量、Ｇａｉｒは筒内吸入空気質量、Ｇｆｕｅｌは筒内吸入燃料質量を示す。
【００４７】
まず、図１〜図３に示す燃料噴射有効パルス幅及び燃料噴射無駄パルス幅設定ルーチンを説明する前に、このルーチンで取入れる各パラメータを設定するルーチンについて説明する。
【００４８】
図５の１０ｍｓ毎ジョブで実行されるルーチンでは、１吸気行程当りの筒内吸入空気質量Ｇａｉｒを算出する。まず、ステップＳ４１で、吸気管圧力センサ２２の検出値に基づいて算出したスロットル弁５ａの下流側の吸気管圧力Ｐと、吸気温センサ４６の検出値に基づいて算出した吸気温度Ｔとを読込み、ステップＳ４２で、現在の吸気温度Ｔが５０℃以下、すなわち通常の吸気温度領域にあるか判断し、Ｔ＜５０℃のときには、ステップＳ４３へ進み、現在の吸気温度Ｔで、吸気温度設定値Ｔ’を設定してステップＳ４５へ進む。また、Ｔ≧５０℃のときには、ステップＳ４４で、上記吸気温度設定値Ｔ’を固定値である５０℃に設定した後、ステップＳ４５へ進む。
【００４９】
ところで、筒内吸入空気質量Ｇａｉｒは、上記吸気管圧力Ｐが一定であっても吸気温度Ｔに影響されて変動するため、空気密度σから筒内吸入空気質量Ｇａｉｒを算出する必要がある。筒内の空気密度σが吸気行程時における吸気チャンバ６Ａ内の空気密度と同じであると考えれば、上記筒内の空気密度αは、空気の状態方程式により以下の式から求めることができる。
【００５０】
σ＝Ｐ×（１−Ｒｅｇｒ）／（Ｒ×Ｔ） …（１）
Ｒｅｇｒ；ＥＧＲ率
Ｒ；ガス定数
図１０に筒内吸入空気質量Ｇａｉｒと吸気温度Ｔとの関係を示す。実線は実測により得た吸入空気質量（実測吸気量）の特性、二点鎖線は上記（１）式により算出した空気密度σに従い求めた吸入空気質量（密度比例吸気量）、一点鎖線は吸気温度Ｔをある温度に固定して算出した空気密度σに従い求めた吸入空気質量（圧力比例吸気量）を示す。図に示すように、吸気温度Ｔが５０℃以下の通常の吸気温度領域では、密度比例吸気量が実測吸気量にほぼ沿う特性を示す。しかし、ターボチャージャ付エンジンでの実測結果によると、吸気温度Ｔが５０℃を越えた領域では、上記筒内吸入空気質量Ｇａｉｒは空気密度σに比例して吸込まれておらず、密度比例吸気量は、実測吸気量に比し過小になり、結果として、空燃比がリーン化してしまう。ところで、本来なら、吸気温度Ｔが５０℃以上の領域では、温度毎の補正係数を用いて補正すべきであるが、制御が複雑化するのと、圧力比例吸気量がほぼ実測吸気量に沿う特性を示しているので、本ルーチンでは吸気温度設定値Ｔ’を５０℃に固定して、筒内吸入空気質量Ｇａｉｒを求めるようにした。尚、上記５０℃は実験により求めた臨界温度であり、採用するエンジン形式により相違するため、個別に設定する必要がある。また、本実施例のようにＥＧＲ装置が搭載されていないエンジンでは、上記ＥＧＲ率はゼロである。
【００５１】
そして、ステップＳ４５へ進むと、上記（１）式の吸気温度Ｔを上記吸気温度設定値Ｔ’に代えて空気密度σを算出し、ステップＳ４６で、上記空気密度σと行程容積Ｄとに基づき次式から１吸気行程当りの理論吸気質量Ｇｔｈを算出する。
【００５２】
Ｇｔｈ←Ｄ×σ …（２）
従って、この理論吸気質量Ｇｔｈは、吸気温度設定値Ｔ’を考慮しなければ、上記吸気管圧力Ｐに比例する。
【００５３】
次いで、ステップＳ４７で、上記理論吸気質量Ｇｔｈを基本として筒内吸入空気質量Ｇａｉｒを次式から算出して、ルーチンを抜ける。
【００５４】
Ｇａｉｒ←（Ｇｔｈ−ηｂ）×ηｖ×Ｋｔｒｍ …（３）
ηｂ：吸気損失質量
ηｖ：体積効率
Ｋｔｒｍ：吸入空気量誤差補正係数
上記理論吸気質量Ｇｔｈと筒内吸入空気質量Ｇａｉｒとの関係は、図１１に示す通りであり、理論吸気質量Ｇｔｈがゼロ、すなわち吸気管圧力Ｐが完全な真空になる前に、上記筒内吸入空気質量Ｇａｉｒがゼロになるような、横軸接点ηｂ、傾きηｖの一次式で表される。但し、実際には、筒内吸入空気質量Ｇａｉｒがゼロになるような運転条件は存在しない。また、１吸気行程当りの吸気損失質量ηｂと体積効率ηｖとは、後述するように、エンジン回転数Ｎｅに基づき一次元マップを補間計算付で参照して設定する。さらに、上記吸入空気量誤差補正係数Ｋｔｒｍは、上記理論吸気質量Ｇｔｈと筒内吸入空気質量Ｇａｉｒとの関係が完全には一次式に乗らないこともあるため、運転条件による空気量計測の定常的な誤差分を補正するためのもので、例えば、スロットル開度αとエンジン回転数Ｎｅとに基づき、二次元マップを補間計算付で参照して設定する。図１２に、この吸入空気量誤差補正係数Ｋｔｒｍを設定する二次元マップの特性を例示する。この図に示すように、この吸入空気量誤差補正係数Ｋｔｒｍの値は基本的には、１．０である。
【００５５】
次に、図６に示す５０ｍｓ毎ジョブで実行される係数設定ルーチンについて説明する。このルーチンでは、上記吸気損失質量ηｂ、及び上記体積効率ηｖを設定する。
【００５６】
まず、ステップＳ５１でエンジン回転数Ｎｅを読込み、ステップＳ５２、Ｓ５３で、このエンジン回転数Ｎｅに基づき一次元マップを補間計算付で参照して、体積効率ηｖと吸気損失質量ηｂとを各々設定してルーチンを抜ける。
【００５７】
上記体積効率ηｖ、上記吸気損失質量ηｂの値は、理論的には、１．０であるが、エンジン回転数毎にカムの同調などの影響によって変化する値であり、また上記体積効率ηｖには、高地補正、排気圧補正等の因子を加味させることができる。
【００５８】
このように、本発明では、吸気管圧力Ｐを基本に筒内吸入空気質量Ｇａｉｒを簡単な計算式に基づいて算出し、また、上記体積効率ηｖ、上記吸気損失質量ηｂは、一次元マップから設定することができるため、コンピュータに掛かる負荷が軽くなり、しかも、上記吸気管圧力Ｐさえ正確に計測されれば、始動時を含む全ての領域で適用することができる。
【００５９】
次に、運転条件による筒内混合気の要求当量比設定ルーチンについて、図７〜図９のルーチンに従って説明する。このルーチンは５０ｍｓ毎のジョブで実行され、定常的な要求増量である可燃限界当量比φｔｗと、最大出力当量比φｆｕｌと、排気ガス温度限界当量比φｔｅｘと、触媒浄化率最良当量比φｇａｓと、燃費率最良当量比φｅｃｏとを設定する。これら各要求当量比φｔｗ，φｆｕｌ，φｔｅｘ，φｇａｓ，φｅｃｏは定常的な要求増量であるため、過渡的な燃料遅れ等によって生じる空燃比変動をフィードフォワード的に補正する係数（エアコン増量補正係数、Ｄレンジ増量補正係数等）、或はセンサやインジェクタの誤差を補正する係数（混合比補正係数等）は含まない。また、このルーチンで設定する各要求当量比は、後述する図１のルーチン実行時に比較され、その最大値が目標当量比（理論燃空比Ｆ／Ａに対する増量係数）ＣＯＥＦとして採用される。尚、この当量比φは、空気過剰率（（Ａ／Ｆ）実際／（Ａ／Ｆ）理論）λに対して、φ＝１／λの関係にある。
【００６０】
まず、ステップＳ６１で冷却水温Ｔｗとエンジン回転数Ｎｅとを読込む。そして、ステップＳ６２〜Ｓ６５或はＳ６６で可燃限界当量比φｔｗを設定する。ステップＳ６２では、上記冷却水温Ｔｗとエンジン回転数Ｎｅとに基づき二次元マップを補間計算付で参照して、可燃限界当量比φｔｗを設定する際の基準となる基準可燃限界当量比Ｍｔｗを設定する。この基準可燃限界当量比Ｍｔｗは、冷却水温Ｔｗとの関係によるリーン側の可燃限界を示すもので、図１３に、この基準可燃限界当量比Ｍｔｗを設定する二次元マップの特性を例示する。この図１３に示すように、冷却水温Ｔｗ側の軸では、従来の水温増量に近い設定がなされており、また、ある程度の高水温では、上記基準可燃限界当量比Ｍｔｗはリーン当量比に設定される。また、エンジン回転数Ｎｅ側の軸では、クランキング回転数を含む全ての回転数領域を検索できるように設定されている。筒内吸入空気質量Ｇａｉｒは始動時においても設定されることから（図５のフローチャート参照）、このときの可燃限界当量比φｔｗが始動時の噴射パルス幅を決定する要因になる。また、始動時において、上記基準可燃限界当量比Ｍｔｗは大きな値に設定されているが、これが始動後にエンジン回転数の増加に従って小さな値に変化して行く際に、後述するように、ある一次遅れ時定数で上記可燃限界当量比φｔｗの変化を遅らせることで始動後の増量が確保される。
【００６１】
次いで、ステップＳ６３へ進むと、エンジン回転数等からエンスト状態にあるかを判断し、エンスト状態にあると判断した場合には、ステップＳ６６へ分岐し、上記基準可燃限界当量比Ｍｔｗで今回の可燃限界当量比φｔｗを設定してステップＳ６７へジャンプする。一方、エンスト状態にないと判断した場合には、ステップＳ６４へ進み、上記基準可燃限界当量比Ｍｔｗと前回の可燃限界当量比φｔｗとを比較し、Ｍｔｗ≧φｔｗのときは、要求当量比が増量されているので、ステップＳ６６へ分岐し、この基準可燃限界当量比Ｍｔｗで今回の可燃限界当量比φｔｗを設定して、増量要求に対応させた後、ステップＳ６７へジャンプする。一方、Ｍｔｗ＜φｔｗのときは、要求当量比が減少されているため、ステップＳ６５へ進み、次式に示す所定時定数（本ルーチンでは、１６×５０ｍｓ）の加重平均により一次遅れ処理を行い、今回の可燃限界当量比φｔｗを設定する。その結果、始動後増量が実質的に確保される。
【００６２】
φｔｗ←（１５φｔｗ＋Ｍｔｗ）／１６ …（４）
次いで、ステップＳ６７〜Ｓ７０或はＳ７１で上記最大出力当量比φｆｕｌが設定される。まず、ステップＳ６７では、スロットル開度αを読込み、ステップＳ６８で、エンジン回転数Ｎｅに基づき、一次元マップを補間計算付で参照してフル増量基準値αｔｈを設定する。このフル増量基準値αｔｈは、エンジンが最大出力を要求しているかを判断する基準値で、図１４に示すように、エンジン回転数Ｎｅとスロットル開度αとの関係から特性を求めてマップ化したものである。
【００６３】
そして、ステップＳ６９で上記スロットル開度αと上記フル増量基準値αｔｈとを比較し、α≦αｔｈのときは、フル増量条件不成立と判断し、ステップＳ７０へ進み、最大出力当量比φｆｕｌを０として、ステップＳ７２へ進む。一方、α＞αｔｈのときは、フル増量条件成立と判断し、ステップＳ７１へ分岐し、上記最大出力当量比φｆｕｌを１．２、すなわち、当該エンジンの最大出力となる当量比（フル増量）に設定する。
【００６４】
尚、目標当量比ＣＯＥＦは、本ルーチンにより設定する各要求当量比φｔｗ，φｆｕｌ，φｔｅｘ，φｇａｓ，φｅｃｏの中から最大値を選択的に設定するものであるため、上記ステップＳ７０で最大出力当量比φｆｕｌが０に設定されると、この最大出力当量比φｆｕｌが目標当量比ＣＯＥＦの選択対象から外れることになる。また、排気ガス温度の上昇を抑える目的の増量は、排気ガス温度限界当量比φｔｅｘとして別に設定されているため、上記最大出力当量比φｆｕｌはエンジン回転数Ｎｅによらず、理論空燃比Ａ／Ｆ（＝１４．６）に対して、おおよそ１５〜２０％増量となるはずである。
【００６５】
次いで、ステップＳ７２〜Ｓ７９或はＳ８０で排気ガス温度限界当量比φｔｅｘを設定する。まず、ステップＳ７２で、負荷の一例である筒内吸入空気質量Ｇａｉｒとエンジン回転数Ｎｅと基づき二次元マップを補間計算付で参照して基本排気ガス温度限界当量比Ｍｔｅｘを設定する。この基本排気ガス温度限界当量比Ｍｔｅｘは、排気ガス温度の上昇を設計限界以下に抑え、エンジン及び排気系を保護するために燃料冷却を行う増量である。図１５に、この基本排気ガス温度限界当量比Ｍｔｅｘを設定する二次元マップの特性を例示する。図に示すように、高負荷（高Ｇａｉｒ）、高回転で大きな増量を行うように設定されている。尚、エンジン負荷を検出するパラメータとして上記筒内吸入空気質量Ｇａｉｒに代えて吸気管圧力Ｐを採用し、二次元マップは吸気管圧力Ｐとエンジン回転数Ｎｅとの関係による特性に基づいて設定するようにしても良い。
【００６６】
そして、ステップＳ７３へ進むと、上記基本排気ガス温度限界当量比Ｍｔｅｘと、上記最大出力当量比φｆｕｌのフル増量値である１．２とを比較し、Ｍｔｅｘ≧１．２のときはステップＳ７４へ進み、またＭｔｅｘ＜１．２のときはステップＳ７６へ進む。
【００６７】
そして、上記ステップＳ７６へ進むと、今回の平均排気ガス温度限界当量比φｔｅｘＡＶを次式に示す加重平均により、所定時定数（本ルーチンでは、１６×５０ｍｓ）の一次遅れ処理を行い、ステップＳ７８へ進む。
【００６８】
φｔｅｘＡＶ←（１５・φｔｅｘＡＶ＋Ｍｔｅｘ）／１６ …（５）
また、ステップＳ７４へ進むと、上記基本排気ガス温度限界当量比Ｍｔｅｘと前回までの平均排気ガス温度限界当量比φｔｅｘＡＶとを比較し、Ｍｔｅｘ＜φｔｅｘＡＶのとき、すなわち、１．２≦Ｍｔｅｘ＜φｔｅｘＡＶのときは、ステップＳ７５へ進み、またＭｔｅｘ≧φｔｅｘＡＶのときは、ステップＳ７７へ進む。
【００６９】
上記ステップＳ７５では、今回の平均排気ガス温度限界当量比φｔｅｘＡＶを次式に示す加重平均により所定時定数（本ルーチンでは、３２×５０ｍｓ）の一次遅れ処理を行い、ステップＳ７８へ進む。
【００７０】
φｔｅｘＡＶ←（３１・φｔｅｘＡＶ＋Ｍｔｅｘ）／３２ …（６）
一方、ステップＳ７７へ進むと、今回の平均排気ガス温度限界当量比φｔｅｘＡＶを次式に示す加重平均により所定時定数（本ルーチンでは１２８×５０ｍｓ）の一次遅れ処理を行って、ステップＳ７８へ進む。
【００７１】
φｔｅｘＡＶ←（１２７・φｔｅｘＡＶ＋Ｍｔｅｘ）／１２８ …（７）
従って、Ｍｔｅｘ＜１．２のいわゆる要求当量比が薄い場合には、比較的早い時定数で平均排気ガス温度限界当量比φｔｅｘＡＶが設定され、また、Ｍｔｅｘ≧１．２、且つＭｔｅｘ≧φｔｅｘＡＶの要求当量比が濃い場合には、ゆっくりとした時定数で上記平均排気ガス温度限界当量比φｔｅｘＡＶが設定される。
【００７２】
そして、ステップＳ７８へ進むと、上記ステップＳ７５，Ｓ７６或はＳ７７で算出した平均排気ガス温度限界当量比φｔｅｘＡＶと上記基本排気ガス温度限界当量比Ｍｔｅｘとを比較し、Ｍｔｅｘ＜φｔｅｘＡＶのときは、ステップＳ７９で、上記基本排気ガス温度限界当量比Ｍｔｅｘを今回の排気ガス温度限界当量比φｔｅｘとして設定し、また、Ｍｔｅｘ≧φｔｅｘＡＶのときは、ステップＳ８０で、今回の排気ガス温度限界当量比φｔｅｘを平均排気ガス温度限界当量比φｔｅｘＡＶで設定する。
【００７３】
従って、排気ガス温度限界当量比φｔｅｘは平均排気ガス温度限界当量比φｔｅｘＡＶで上限リミットされ、また、基本排気ガス温度限界当量比Ｍｔｅｘで下限リミットされる。その結果、上記排気ガス温度限界当量比φｔｅｘの上限リミット値は、エンジン回転数Ｎｅと筒内吸入空気質量Ｇａｉｒとに応じて、すなわち、上記ステップＳ７５，Ｓ７６或はＳ７７で設定される平均排気ガス温度限界当量比φｔｅｘＡＶに従って、ゆっくり変化したり、あるいは比較的早く変化したりするようになる。その結果、長時間の全開運転において排気ガス温度が実際に設計限界値に近づくまではφｔｅｘ＜φｆｕｌとなり、従来のように、最大出力当量比φｆｕｌの中に、燃料冷却成分を含ませていた場合に比し、出力の増加が期待できる。尚、この排気ガス温度限界当量比φｔｅｘを排気ガス温度を計測或は推測した値に基づいて設定するようにしても良い。
【００７４】
次いで、上記ステップＳ７９或はＳ８０からステップＳ８１へ進み、このステップＳ８１〜Ｓ８５或はＳ８６で触媒浄化率最良当量比φｇａｓが設定される。この触媒浄化率最良当量比φｇａｓは、触媒コンバータ１１として３元触媒が採用されている場合に対象となり、排気ガスを浄化しなければならない条件が成立したときに、φｇａｓ←１．０とし、それ以外では、φｇａｓ←０とする。
【００７５】
ステップＳ８１〜Ｓ８４では、排気ガスを浄化しなければならない運転条件にあるかを判断する。この判断条件は、スロットル開度αとエンジン回転数Ｎｅとで行い、スロットル開度αが、１０゜≦α＜７０゜（ステップＳ８１，Ｓ８２）で、且つエンジン回転数Ｎｅが、８００ｒｐｍ≦Ｎｅ＜３０００ｒｐｍ（ステップＳ８３，Ｓ８４）のとき、触媒浄化率最良当量比φｇａｓを１．０に設定し（ステップＳ８５）、それ以外のときは、φｇａｓを０に設定する（ステップＳ８６）。尚、上記判断条件は一例であり、排気ガス浄化領域はエンジンの特性等に応じて任意に設定でき、また例えばスロットル開度αに代えて、吸気管圧力Ｐや筒内吸入空気質量Ｇａｉｒを採用しても良い。また、上記触媒コンバータ１１としてリーンＮｏｘ触媒が採用されている場合には、触媒浄化率最良当量比φｇａｓは全運転領域でφｇａｓ←０に設定する。
【００７６】
その後、ステップＳ８７で、負荷の一例としての吸気管圧力Ｐとエンジン回転数Ｎｅとに基づいて二次元マップを補間計算付で参照して、燃費率最良当量比φｅｃｏを設定して、ルーチンを抜ける。
【００７７】
図１６に、この燃費率最良当量比φｅｃｏを設定する二次元マップの特性を例示する。図に示すように、この二次元マップの各領域には、当該エンジン運転条件下での最良の燃費率を得ることのできる当量比が予め実験などから求めて格納されている。
【００７８】
ここで、この燃費率最良当量比φｅｃｏと前述の可燃限界当量比φｔｗとの関係について簡単に説明する。
【００７９】
空燃比を可燃限界当量比φｔｗまでリーン化させると一般には最良の燃費を得ることができなくなる。しかし、暖機途中の低水温時などでは空燃比を燃費率最良当量比φｅｃｏまでリーン化するよりも、上記可燃限界当量比φｔｗにより、リーン化を抑制した方が良好な燃焼を得ることができる場合もあるので、上記可燃限界当量比φｔｗとは別に燃費率最良当量比φｅｃｏが設定されている。尚、上記触媒浄化率最良当量比φｇａｓが全運転領域で、φｇａｓ←１．０の場合には、上記燃費率最良当量比φｅｃｏは隠れてしまうため、その意味を失うが、例えば、排気浄化システムがリーンＮｏｘ触媒を採用している場合には、前述のように、触媒浄化率最良当量比φｇａｓが全運転領域でφｇａｓ←０に設定されるため、この燃費率最良当量比φｅｃｏは触媒浄化率最良当量比としての要素も合わせ持つことになる。また、ＥＧＲを用いて理論空燃比制御により燃費向上を行う際には、上記（１），（２）式から明らかなように、理論吸気質量Ｇｔｈには、ＥＧＲ量が除かれた新気量のみであるため、全運転領域で、φｅｃｏ←１．０としたほうが、空燃比制御性が良くなる。また、この燃費率最良当量比φｅｃｏを、吸気管圧力Ｐに代えて筒内吸入空気質量Ｇａｉｒと、エンジン回転数Ｎｅとの関係により作成した二次元マップを参照して設定するようにしても良い。
【００８０】
そして、上記各ルーチンで求めた設定値を用いて、図１〜図３の１０ｍｓ毎ジョブで実行される燃料噴射有効パルス幅及び燃料噴射無駄パルス幅設定ルーチンのステップＳ１〜Ｓ１０で１サイクル当りの燃料噴射質量を算出し、ステップＳ１１〜Ｓ２４で燃料噴射有効パルス幅Ｔｅと燃料噴射無駄パルス幅Ｔｓとを算出する。
【００８１】
まず、ステップＳ１では、上記筒内混合気の要求当量比設定ルーチンで設定した各要求当量比φｔｗ，φｆｕｌ，φｔｅｘ、φｇａｓ，φｅｃｏを比較し、これらの要求値の最大値を目標当量比ＣＯＥＦとして決定する。
【００８２】
次いで、ステップＳ２で筒内への１吸気行程当りの筒内吸入燃料質量Ｇｆｕｅｌを次式から算出する。
【００８３】
Ｇｆｕｅｌ←Ｇａｉｒ×Ｆ／Ａ×ＣＯＥＦ×Ｋｆｂ …（８）
Ｆ／Ａ：理論燃空比
Ｋｆｂ：空燃比フィードバック補正係数
ここで、理論燃空比Ｆ／Ａは理論空燃比Ａ／Ｆの逆数で、理論空燃比が１４．６であれば、理論燃空比は１／１４．６となる。また、この理論空燃比Ａ／Ｆは、燃料を完全燃焼するに必要な最小の空気量と燃料量との比であり、多種燃料に対応したエンジン等では変数となる。
【００８４】
一方、上記空燃比フィードバック補正係数Ｋｆｂの初期値は１．０であり、空燃比フィードバック制御、及び空燃比学習制御を行う外部のストラテジーによって書換えられる。また、理論燃空比Ｆ／Ａに目標当量比ＣＯＥＦを掛けることで、運転条件に応じた筒内混合気の目標空燃比が設定される。従って、リニアＡ／Ｆセンサでフィードバック制御する際に設定される上記空燃比フィードバック補正係数Ｋｆｂの目標値は上記理論燃空比Ｆ／Ａを１とした場合の上記目標当量比ＣＯＥＦである。
【００８５】
ところで、上記筒内吸入燃料質量Ｇｆｕｅｌは、筒内へ１行程で吸込んで欲しい燃料量であり、定常的にはインジェクタ２５からの噴射量と一致するが、過渡時には応答遅れがあるため一致しない。すなわち、過渡的にはインジェクタ２５からの燃料噴射量が急増しても、燃料の一部が吸気ポート内壁に付着してから筒内に吸込まれるので、その分遅れが生じ、筒内に吸込まれる燃料量はゆっくりと増加する。従って、上記筒内吸入燃料質量Ｇｆｕｅｌが急増したときには、インジェクタ２５からの燃料量を、この筒内吸入燃料質量Ｇｆｕｅｌよりも多めにフィードフォワード噴射させて、筒内への吸入燃料量を筒内吸入燃料質量Ｇｆｕｅｌと一致させる必要がある。この追加増量分に相当する過渡時の１吸気行程当りの燃料追加質量Ｇａｃｃを、次のステップＳ３〜Ｓ９で算出する。
【００８６】
まず、ステップＳ３では、スロットル開度αとエンジン回転数Ｎｅとに基づき二次元マップ（図１７参照）を補間計算付で参照して、静的な筒内吸入空気質量Ｇａｉｒに相当する指標値（静的指標値）Ｍａｃｃを設定する。この静的指標Ｍａｃｃは、負荷変動の少ない定常走行で、スロットル開度αとエンジン回転数Ｎｅを一定に保持すれば、最終的には筒内吸入空気質量Ｇａｉｒが一定になるということを前提にしている。従って、走行中の負荷変動が無視できるとすれば、図４のステップＳ３１に示すように、上記静的指標Ｍａｃｃをスロットル開度αに基づき一次元マップから補間計算付で設定することも可能である。
【００８７】
また、過渡的には、上記筒内吸入空気質量Ｇａｉｒにも遅れがあり、この遅れに対応した燃料追加質量Ｇａｃｃを求める必要がある。
【００８８】
そこで、ステップＳ４では、まず、吸気チャンバ６Ａ（図２２参照）内の吸気遅れを想定した所定時定数（運転条件によって異なるが本ルーチンでは４×１０ｍｓに固定している）の一次遅れ処理を、次式に示す加重平均に基づいて行い、動的な筒内吸入空気質量Ｇａｉｒに相当する指標値（動的指標値）Ｓａｃｃを算出する。
【００８９】
Ｓａｃｃ←（３・Ｓａｃｃ＋Ｍａｃｃ）／４ …（９）
次いで、ステップＳ５で、上記動的指標値Ｓａｃｃと、前回のルーチン実行時に算出した燃料付着による遅れ指標値Ｔａｃｃとを比較する。尚、この遅れ指標値Ｔａｃｃは上記筒内吸入燃料質量Ｇｆｕｅｌに相当する。
【００９０】
そして、Ｓａｃｃ＜Ｔａｃｃ、すなわち、前回の筒内吸入燃料質量Ｇｆｕｅｌに相当する値（Ｔａｃｃ）に対して今回の筒内吸入空気質量Ｇａｉｒに相当する値（Ｓａｃｃ）が少ないときには、ステップＳ６へ進み、上記遅れ指標値Ｔａｃｃを上記動的指標値Ｓａｃｃとして、ステップＳ８へ進む。その結果、Ｓａｃｃ＝Ｔａｃｃとなり、後述のステップＳ９では、上記燃料追加質量Ｇａｃｃ＝０となり、燃料が追加増量されない。
【００９１】
一方、Ｓａｃｃ≧Ｔａｃｃ、すなわち、前回の筒内吸入燃料質量Ｇｆｕｅｌに相当する値（Ｔａｃｃ）に対して今回の筒内吸入空気質量Ｇａｉｒに相当する値（Ｓａｃｃ）が増加しているか、或は変化していないときには、ステップＳ７へ分岐して、この動的指標Ｓａｃｃを、次式に示す加重平均により、燃料付着による時定数（本ルーチンでは、８×１０ｍｓ）の一次遅れ処理を行って、今回の燃料付着による遅れ指標値Ｔａｃｃを算出して、ステップＳ８へ進む。
【００９２】
Ｔａｃｃ←（７・Ｔａｃｃ＋Ｓａｃｃ）／８ …（１０）
そして、ステップＳ８へ進むと、エンジン回転数Ｎｅに基づき一次元マップを補間計算付で参照して係数Ｒａｃｃを設定する。
【００９３】
その後、ステップＳ９へ進み、過渡時の１吸気行程当りの燃料追加質量Ｇａｃｃを次式から算出する。
【００９４】
Ｇａｃｃ←Ｒａｃｃ×（Ｓａｃｃ−Ｔａｃｃ） …（１１）
この（１１）式に示すように、今回の筒内吸入空気質量Ｇａｉｒに相当する値（Ｓａｃｃ）と今回の筒内吸入燃料質量Ｇｆｕｅｌに相当する値（Ｔａｃｃ）との差分が、燃料付着による不足分を補填する増量パターンということになる。
【００９５】
このように、この燃料追加質量Ｇａｃｃは、スロットル開度αとエンジン回転数Ｎｅのみから、簡易的に筒内吸入空気質量Ｇａｉｒの変化量だけを遅れなく検出し、この値を利用して算出している。
【００９６】
図１８にスロットル弁５ａを走行中に僅かに開弁させたときの上記各指標値Ｍａｃｃ，Ｓａｃｃ，Ｔａｃｃの特性を示す。静的指標値Ｍａｃｃはスロットル開度αに追従して増加し、一方、動的指標値Ｓａｃｃは上記静的指標値Ｍａｃｃに対して４×１０ｍｓの一次遅れで増加する。さらに、上記遅れ指標値Ｔａｃｃは、この動的指標値Ｓａｃｃに対して８×１０ｍｓの時定数で増加する。この遅れ指標値Ｔａｃｃが筒内吸入燃料質量Ｇｆｕｅｌに対応していると考えれば、付着分の遅れは、上記動的指標値Ｓａｃｃと上記遅れ指標値Ｔａｃｃとの差分（図のハッチングで示す領域）に相当することになり、上記ステップＳ９では、この差分指標値（Ｓａｃｃ−Ｔａｃｃ）に上記係数Ｒａｃｃを掛けて、燃料追加質量Ｇａｃｃを求めている。ところで、図１９に示すように、この差分指標値（Ｓａｃｃ−Ｔａｃｃ）が、従来のスロットル弁５ａの開弁直後に生じるリーンスパイクのパターンに符合することが、実験から明らかになった。従って、この差分指標値（Ｓａｃｃ−Ｔａｃｃ）相応分を加速増量すれば、排気空燃比は過渡時においてもリーン化することなく一定となる。
【００９７】
上記燃料追加質量Ｇａｃｃは加速増量・加速時追加パルスに近似するが、この燃料追加質量Ｇａｃｃはあくまでも筒内への吸入燃料量を筒内吸入燃料質量Ｇｆｕｅｌにするためのフィードフォワード分であり、加速時に空燃比Ａ／Ｆを一時的にリッチにするためのものではない。従って、もし、加速時に空燃比Ａ／Ｆをリッチにする必要があれば、前記目標当量比ＣＯＥＦの選択枝として、前記図７〜図９に示す要求当量比設定ルーチン中に、加速当量比φａｃｃを第６の増量要因として加えることも、当然考えられる。
【００９８】
また、従来、上記燃料追加質量Ｇａｃｃを算出する方法として三角増量法がある。この三角増量法はスロットル弁の動きを基本に増量値が設定されるため応答性はよいが、筒内へ吸入される燃料量の一次遅れを三角形で近似させながら増量するので、上記図１９に示すリーンスパイクに対して正確に符合せず、過渡時の空燃比が部分的にリッチになったりリーンになったりしてしまう。また、理論的な燃料付着モデルは吸入空気量を計測するセンサ自体の応答遅れやノイズがあるため、性能にばらつきが生じ易く、充分な信頼性が得られない。
【００９９】
そして、ステップＳ１０へ進むと、上記筒内吸入燃料質量Ｇｆｕｅｌに上記燃料追加質量Ｇａｃｃを加算して、１サイクル当りの燃料噴射質量Ｇｉｎｊを算出する。ところで、噴射時期が早期に設定されていて、噴射が終了した後でまだ吸気行程になる前のタイミングで上記燃料噴射質量Ｇｉｎｊが急増した場合には、自動的に追加噴射することも可能で、この燃料噴射質量Ｇｉｎｊは、ある吸気行程から１回転半前、すなわち吸気行程終了直後から次の吸気行程開始までの期間にインジェクタ２５から噴射される燃料の総量が算出される。
【０１００】
次に、ステップＳ１１〜Ｓ２４で、気筒別燃料噴射有効パルス幅Ｔｅ１〜Ｔｅ４、及び気筒別燃料噴射無駄パルス幅Ｔｓ１〜Ｔｓ４が算出される（ここで、１，２，３，４は気筒番号を示す）。
【０１０１】
まず、ステップＳ１１で、最大出力当量比φｆｕｌの値を参照し、ステップＳ１２で排気ガス温度限界当量比φｔｅｘの値を参照する。そして、φｆｕｌ＝０のエンジンが最大出力を要求しておらず、しかも、φｔｅｘ＝０の排気ガス温度限界による増量も要求していない場合には、ステップＳ１３へ進み、空燃比の気筒別のばらつきを補正する気筒別噴射量補正係数Ｋｔｎ１〜Ｋｔｎ４を１．０に設定して、ステップＳ１８へ進む。一方、上記ステップＳ１１で、φｆｕｌ＝１．２（フル増量）、あるいは、ステップＳ１２で、φｔｅｘ≠０と判断されている場合には、ステップＳ１４へ分岐し、ステップＳ１４〜Ｓ１７でエンジン回転数Ｎｅに基づき、気筒別の一次元マップを補間計算付で参照して、上記気筒別噴射量補正係数Ｋｔｎ１〜Ｋｔｎ４を各々設定して、ステップＳ１８へ進む。噴射量補正係数Ｋｔｎは燃料噴射量を気筒別に増減量して空燃比を均一に設定するもので、上記各一次元マップは、各気筒の吸気特性に対応する値が格納されている。
【０１０２】
ところで、運転条件によって燃料噴射量を気筒毎に変えなければならない原因は、本来、筒内吸入空気質量Ｇａｉｒが気筒毎に相違するからであり、従って、正確には筒内吸入空気質量Ｇａｉｒを吸気特性を考慮して気筒毎に算出し、この算出結果から筒内吸入燃料質量Ｇｆｕｅｌ及び燃料噴射質量Ｇｉｎｊを気筒毎に算出すべきであるが、ＣＰＵの計算負荷が重くなるため、上記気筒別噴射量補正係数Ｋｔｎ１〜Ｋｔｎ４を用いて、見かけ上、気筒別燃料噴射有効パルス幅Ｔｅ１〜Ｔｅ４［ｍｓ］を補正するようにした。
【０１０３】
また、空燃比の気筒毎のばらつきが特に問題となるのは、高負荷運転時のノッキングである。中低負荷運転での空燃比のばらつきは排気ガスに余り影響せず、エンジン振動の観点からすると全気筒の燃料量が一致している方が却って各気筒の出力が均等になるので望ましい。従って、フル増量が要求されている場合には、全気筒が最大出力が得られるようにする必要があり、また、排気ガス温度限界による増量が要求されている場合には、全気筒の排気系等の高温化を設計限界以内に抑制してエンジンを保護する必要があるため、このような場合、すなわち、気筒間の空燃比が均一であることを前提とする要求増量が設定されている場合のみ、エンジン回転数毎に各気筒の燃料噴射量を適正に増減させる補正係数（気筒別噴射量補正係数）を設定し、各気筒の筒内混合気の空燃比が目標空燃比となるように個別に制御する。
【０１０４】
そして、ステップＳ１８へ進むと、全体としての燃料噴射有効パルス幅Ｔｅ＿ａｌｌを次式から算出する。
【０１０５】
Ｔｅ＿ａｌｌ←Ｋｍｒ×Ｇｉｎｊ×Ｋｃｏｎ …（１２）
Ｋｍｒ：パルス幅誤差補正係数
Ｋｃｏｎ：インジェクタ容量係数［ｓｅｃ／ｇ］ここで、パルス幅誤差補正係数Ｋｍｒは、インジェクタ２５の動的な流量特性の非線型性の補正であり、本ルーチンでは、エンジン回転数Ｎｅと燃料噴射質量Ｇｉｎｊとに基づき、図２０に示す二次元マップを補間計算付で参照して設定する。また、上記インジェクタ容量係数Ｋｃｏｎは、インジェクタ２５の静的質量流量特性の逆数であり、１グラムの燃料を０．１秒間に噴射するインジェクタでは、Ｋｃｏｎ＝０．１［ｓｅｃ／ｇ］である。なお、このインジェクタ容量係数Ｋｃｏｎは、本ルーチンでは固定値としているが、多種燃料に対応可能なＦＦＶエンジン等で、燃料の比重や粘性が変化する場合には変数として用いる。
【０１０６】
次いで、ステップＳ１９〜Ｓ２２で、上記燃料噴射有効パルス幅Ｔｅ＿ａｌｌを上記気筒別噴射量補正係数Ｋｔｎ１〜Ｋｔｎ４で補正して、気筒別燃料噴射有効パルス幅Ｔｅ１〜Ｔｅ４をそれぞれ設定する。
【０１０７】
その後、ステップＳ２３へ進み、バッテリ電圧ＶＢに基づき、一次元マップを補間計算付で参照し燃料噴射無駄パルス幅Ｔｓ＿ａｌｌ［ｍｓ］を設定し、ステップＳ２４で、気筒別燃料噴射無駄パルス幅Ｔｓ１〜Ｔｓ４を、上記燃料噴射無駄パルス幅Ｔｓ＿ａｌｌで設定してルーチンを抜ける。
【０１０８】
ＯＳ側では、上記気筒別燃料噴射有効パルス幅Ｔｅ１〜Ｔｅ４と上記気筒別燃料噴射無駄パルス幅Ｔｓ１〜Ｔｓ４とをそれぞれ加算して、燃料噴射パルス幅を気筒別に算出する。
【０１０９】
インジェクタ２５に与える電圧パルスの幅（燃料噴射パルス幅）と、燃料噴射質量Ｇｉｎｊとの関係は、図２１に示す通りであり、燃料噴射質量Ｇｉｎｊが決定されれば、基本的に、この燃料噴射質量Ｇｉｎｊにインジェクタ容量係数Ｋｃｏｎを掛けることで上記有効燃料噴射パルス幅Ｔｅが求められ、この有効燃料噴射パルス幅Ｔｅに無効噴射パルス幅である燃料噴射無駄パルス幅Ｔｓを加算することで燃料噴射パルス幅が算出される。
【０１１０】
本実施例では、燃料噴射質量Ｇｉｎｊ分の燃料を１回で噴射するか、２回に分割して噴射するかはＯＳ内部で決定されるため、上記気筒別燃料噴射有効パルス幅Ｔｅ１〜Ｔｅ４と上記気筒別燃料噴射無駄パルス幅Ｔｓ１〜Ｔｓ４とを、予め加算せずに別々にＯＳに渡し、このＯＳ内部で噴射設定の直前に加算して、各インジェクタに与える燃料噴射パルス幅を算出している。こうすることで、ＯＳ内部で、燃料噴射質量Ｇｉｎｊ分の燃料を２回に分割して噴射すると決定した場合でも、１／２の燃料噴射有効パルス幅Ｔｅに燃料噴射無駄パルス幅Ｔｓを加算することで、１回当りの燃料噴射パルス幅を簡単に与えることができる。
【０１１１】
このように、本実施例によれば、例えば新しい燃料付着モデルを開発した際は、筒内吸入空気質量Ｇａｉｒを算出する前記（３）式中に演算式を追加すれば良く、またインジェクタの流量が非線型であるならば燃料噴射有効パルス幅Ｔｅ，燃料噴射無駄パルス幅Ｔｓを補正するだけで良いことになる。その結果、ある運転領域での空燃比制御性に不具合が生じた場合に、どの部分を変更すべきかが明確となるばかりか、その変更が他の部分に影響を及ぼし難いため、ベース制御ソフトとしての活用が期待できる。
【０１１２】
【発明の効果】
以上、説明したように第１発明によれば、吸気管圧力さえ正確に計測されれば、少なくとも、この吸気管圧力と吸気温度、及び吸気損失量と体積効率とに基づいて一次式から筒内吸入空気量を簡単に導くことができ、しかも、上記吸気損失量と体積効率とがエンジン回転数に対する特性を求めて一次元マップ化されているので、二次元マップを必要とする従来の吸入空気量の演算過程に比し、コンピュータの負荷を軽減し、しかもメモリ容量の有効利用を図ることができる。
【０１１３】
また、第２発明によれば、上記吸気温度が、上記筒内吸入空気量が吸気温度に依存して変動しない温度として予め設定した臨界温度以上のときには、この吸気温度を固定値で設定して吸気密度を設定するようにしたので、ターボチャージャ付エンジン等で、上記吸気温度が高温になっても筒内吸入空気量を実際の吸入空気量にほぼ等しい値に設定することができて空燃比のリーン化を防止することができるばかりでなく、良好な過渡応答性を得ることができる。
【図面の簡単な説明】
【図１】１０ｍｓ毎ジョブで実行される燃料噴射有効パルス幅及び燃料噴射無駄パルス幅設定ルーチンを示すフローチャート
【図２】１０ｍｓ毎ジョブで実行される燃料噴射有効パルス幅及び燃料噴射無駄パルス幅設定ルーチンを示すフローチャート（続き）
【図３】１０ｍｓ毎ジョブで実行される燃料噴射有効パルス幅及び燃料噴射無駄パルス幅設定ルーチンを示すフローチャート（続き）
【図４】他の態様による燃料噴射有効パルス幅及び燃料噴射無駄パルス幅設定ルーチンの要部を示すフローチャート
【図５】１０ｍｓ毎ジョブで実行される筒内吸入空気質量設定ルーチンを示すフローチャート
【図６】５０ｍｓ毎ジョブで実行される係数設定ルーチンを示すフローチャート
【図７】５０ｍｓ毎ジョブで実行される要求当量比設定ルーチンを示すフローチャート
【図８】５０ｍｓ毎ジョブで実行される要求当量比設定ルーチンを示すフローチャート（続き）
【図９】５０ｍｓ毎ジョブで実行される要求当量比設定ルーチンを示すフローチャート（続き）
【図１０】筒内吸入空気質量と吸気温度との関係を示す特性図
【図１１】理論吸気質量と筒内吸入空気質量との関係を示す特性図
【図１２】吸入空気量誤差補正係数を設定する二次元マップの概念図
【図１３】基準可燃限界当量比を設定する二次元マップの概念図
【図１４】最大出力当量比の増量領域を示す特性図
【図１５】基本排気ガス温度限界当量比を設定する二次元マップの特性図
【図１６】燃費率最良当量比を設定する二次元マップの特性図
【図１７】静的な筒内吸入空気質量に相当する指標値を設定する二次元マップの概念図
【図１８】静的指標値と動的指標値とこの動的指標値に対する遅れ指標値との関係を示す特性図
【図１９】過渡運転時のリーンスパイクと増量パターンとの関係を示す特性図
【図２０】パルス幅誤差補正係数を設定する二次元マップの特性図
【図２１】燃料噴射パルス幅と燃料噴射質量との関係を示す特性図
【図２２】エンジンの吸気系モデルを示す説明図
【図２３】エンジンの全体構成図
【図２４】クランクロータとクランク角センサの正面図
【図２５】カムロータとカム角センサの正面図
【図２６】電子制御装置の回路構成図
【符号の説明】
５ａ…スロットル弁
Ｄ…行程容積
Ｇａｉｒ…筒内吸入空気質量（吸入空気量）
Ｇｆｕｅｌ…筒内吸入燃料質量
Ｇｔｈ…理論吸気質量（理論吸気量）
Ｎｅ…エンジン回転数
Ｐ…吸気管圧力
Ｔ…吸気温度
σ…空気密度
ηｂ…吸気損失質量（吸気損失量）
ηｖ…体積効率
５０℃…臨界温度、固定値
【産業上の利用分野】
本発明は、Ｄジェトロニック方式により吸入空気量を算出するエンジンに関し、詳しくはスロットル弁下流の空気密度から筒内吸入空気量を設定するエンジンの吸入空気量検出方法に関する。
【０００２】
【従来の技術】
一般に、この種の燃料噴射制御では、吸入空気量を吸入空気量センサの計測値に基づいて算出するＬジェトロニック方式と、スロットル弁下流の吸気管圧力に基づいて算出するＤジェトロニック方式とがあり、その中でも、最近のＤジェトロニック方式を採用するエンジンでは、筒内へ実際に吸入される空気量を過渡時であっても正確に推定することのできる技術が種々提案されている。
【０００３】
例えば、特開平６−１８５３９１号公報では、筒内への吸入空気量を空気の状態方程式を用いて推定している。すなわち、スロットル弁下流のチャンバ容積に対する吸入空気の入出力関係に着目し、まず、スロットル弁を単位時間（Δｔ）当りに通過する吸気量（スロットル通過空気量）Ｇｔｒ・Δｔを設定し、次いで、スロットル弁下流のチャンバ内空気量Ｇｂの圧力変化分に基づいて今回チャンバに充填された空気量の変化分ΔＧｂを設定する。そして、今回チャンバに充填された空気量分は、筒内へ吸入されなかったものと見なして、単位時間当りの筒内吸入空気量Ｇｃを算出している。これを式に表せば以下の通である。
【０００４】
Ｇｃ＝Ｇｔｒ×Δｔ−ΔＧｂ
尚、上記スロットル通過空気量Ｇｔｒは、スロットル弁の上下流の圧力差、すなわち大気圧と上記吸気管圧力との偏差、スロットル有効開口面積、及びスロットル弁上流側の空気密度等に基づき近似的に算出している。
【０００５】
【発明が解決しようとする課題】
ところで、上記先行技術によるＤジェトロニック方式を用いた筒内吸入空気量の算出過程では、以下の理由からコンピュータに掛かる負荷が重くなる。
【０００６】
（１）スロットル下流側チャンバへの流入吸気量と、このチャンバに充填された吸入空気量とを個々に求めるなど、複雑な演算を必要としている。
【０００７】
（２）例えばスロットル通過空気量Ｇｔｒを設定する関数として必要なスロットル有効開口面積Ａは、係数Ｃとスロットル投影面積Ｓとの積で設定しているが、この係数Ｃを設定する際に、スロットル開度θＴＨと吸気管圧力Ｐとに基づいて二次元マップを参照するようにしている。
【０００８】
一方、上記先行技術では、チャンバ内空気量Ｇｂ及びこのチャンバ内空気量の変化分ΔＧｂは、吸気管圧力Ｐの変化分ΔＰが一定であれば、ほぼ空気温度Ｔに反比例した値に設定されるが、実測結果によれば、ターボエンジンなどでは、チャンバ内の吸気温度がある温度を越えると、吸入空気量は吸気温度に依存して減少せず、従って、全運転領域の吸入空気量を吸気温度に反比例した特性で設定すると、高回転、高負荷運転領域などでは、計算により得られた吸入空気量が実際の吸入空気量に対して過小になり、この過小な吸入空気量に基づいて燃料噴射量が設定されるため、空燃比がリーン化してしまう課題がある。
【０００９】
本発明は、上記事情に鑑みてなされたもので、一次式及び一次元マップにより筒内への吸入空気量を簡単に算出することができて、コンピュータに対する負担を軽減することができ、しかも、スロットル弁下流の吸気温度が高温になっても筒内への吸入空気量を実際の吸入空気量とほぼ等しい値に設定することができて、空燃比のリーン化を防止するとともに、過渡応答性の良いエンジンの吸入空気量検出方法を提供することを目的としている。
【００１０】
【課題を解決するための手段】
上記目的を達成するため第１発明によるエンジンの吸入空気量検出方法は、少なくともスロットル弁下流の吸気管圧力と吸気温度とに基づき空気の状態方程式から空気密度を設定し、この空気密度と行程容積とから１吸気行程当りの筒内への理論吸気量を設定し、またエンジン回転数に基づき１吸気行程当りの筒内への吸気損失量と体積効率とを各々に対応する一次元マップを参照して設定し、上記理論吸気量と上記吸気損失量との差分と上記体積効率とから１吸気行程当りの筒内への吸入空気量を設定することを特徴とする。
【００１１】
上記目的を達成するため第２発明によるエンジンの吸入空気量検出方法は、前記吸気温度が、前記筒内への吸入空気量が吸気温度に依存して変動しない温度として予め設定した臨界温度以上のときには、上記吸気温度を固定値で設定して上記吸気密度を設定することを特徴とする。
【００１２】
【作 用】
第１発明では、筒内への吸入空気量を算出するに際して、まず、吸気行程時の筒内の空気密度を、少なくともスロットル弁下流の吸気管圧力と吸気温度とに基づき空気の状態方程式から算出する。次いでこの空気密度と行程容積との積から、１吸気行程当りの筒内への理論吸気量を算出する。一方、エンジン回転数に基づいて１吸気行程当りの吸気損失量と体積効率を各々に対応する一次元マップを参照して設定し、上記理論吸気量と上記吸気損失量との差と、上記体積効率との積から、１吸気行程当りの筒内吸入空気量を算出する。
【００１３】
第２発明では、筒内吸入空気量が吸気温度に依存して変動しなくなる温度を臨界温度として予め設定し、上記吸気温度が臨界温度を越えたときには、実際の吸気温度に拘らず、上記固定値と上記吸気管圧力とに基づいて吸気密度を設定する。
【００１４】
【実施例】
以下、図面に基づき本発明の実施例を説明する。図２３にはエンジンの全体概略図が示されており、図中の符号１はエンジン（図においては水平対向４気筒型エンジンを示す）を示し、シリンダヘッド２の吸気ポート２ａにインテークマニホルド３が連通され、このインテークマニホルド３の上流にエアーチャンバ４を介してスロットル通路５が連通されている。このスロットル通路５の上流側には、吸気管６を介してエアクリーナ７が取付けられ、このエアクリーナ７が吸入空気の取り入れ口であるエアインテークチャンバ８に連通されている。また、上記排気ポート２ｂにエキゾーストマニホルド９を介して排気管１０が連通され、この排気管１０に触媒コンバータ１１が介装されてマフラ１２に連通されている。
【００１５】
一方、上記スロットル通路５にスロットル弁５ａが設けられ、このスロットル通路５の直上流の上記吸気管６にインタークーラ１３が介装され、さらに、上記吸気管６の上記エアクリーナ７の下流側にレゾネータチャンバ１４が介装されている。
【００１６】
また、上記レゾネータチャンバ１４と上記インテークマニホルド３とを連通して上記スロットル弁５ａの上流側と下流側とをバイパスするバイパス通路１５に、アイドル空気量を調整するアイドル制御（ＩＳＣ）弁１６が介装されている。さらに、上記ＩＳＣ弁１６の直下流側に、吸気圧が負圧のとき開弁し、またターボチャージャ１８によって過給されて吸気圧が正圧になったとき閉弁するチェックバルブ１７が介装されている。
【００１７】
上記ターボチャージャ１８は、上記吸気管６の上記レゾネータチャンバ１４の下流側にコンプレッサが介装され、タービンが上記排気管１０に介装されている。さらに、上記ターボチャージャ１８のタービンハウジング流入口には、ウエストゲート弁１９が介装され、このウエストゲート弁１９には、ウエストゲート弁作動用アクチュエータ２０が連設されている。
【００１８】
上記ウエストゲート弁作動用アクチュエータ２０は、ダイヤフラムにより２室に仕切られ、一方がウエストゲート弁制御用デューティソレノイド弁２１に連通される圧力室を形成し、他方が上記ウエストゲート弁１９を閉方向に付勢するスプリングを収納したスプリング室を形成している。
【００１９】
上記ウエストゲート弁制御用デューティソレノイド弁２１は、上記レゾネータチャンバ１４と上記吸気管６の上記ターボチャージャ１８のコンプレッサ下流とを連通する通路に介装されており、後述する電子制御装置５０（ＥＣＵ；図２６参照）から出力される制御信号のデューティ比に応じて、上記レゾネータチャンバ１４側の圧力と上記コンプレッサ下流側の圧力とを調圧し、上記ウエストゲート弁作動用アクチュエータ２０の圧力室に供給する。
【００２０】
すなわち、上記ＥＣＵ５０によって上記ウエストゲート弁制御用デューティソレノイド弁２１を制御し、上記ウエストゲート弁作動用アクチュエータ２０を作動させて上記ウエストゲート弁１９による排気ガスリリーフを調整することにより、上記ターボチャージャ１８による過給圧を制御するようになっている。
【００２１】
また、上記インテークマニホルド３に、吸気管圧力Ｐを絶対圧で検出する吸気管圧力センサ２２が通路２３を介して連通され、さらに、上記インテークマニホルド３の各気筒の各吸気ポート２ａの直上流側にインジェクタ２５が臨まされている。また、上記シリンダヘッド２の各気筒毎に、その先端を燃焼室に露呈する点火プラグ２６ａが取付けられ、この点火プラグ２６ａに各気筒毎に配設された点火コイル２６ｂを介してイグナイタ２７が接続されている。
【００２２】
上記インジェクタ２５には、燃料タンク２８内に設けたインタンク式の燃料ポンプ２９から燃料フィルタ３０を経て燃料が圧送され、プレッシャレギュレータ３１にてインジェクタ２５への燃料圧力が調圧される。
【００２３】
また、上記スロットル弁５ａに、スロットル開度センサ３３ａとアイドルスイッチ３３ｂとを内蔵したスロットルセンサ３３が連設され、上記エアーチャンバ４に吸気温センサ４６が臨まされている。
【００２４】
また、上記エンジン１のシリンダブロック１ａにノックセンサ３４が取付けられるとともに、このシリンダブロック１ａの左右両バンクを連通する冷却水通路３５に水温センサ３６が臨まされ、上記排気管１０の上記エキゾーストマニホルド９の集合部にＯ２ センサ３７が臨まされている。
【００２５】
また、上記シリンダブロック１ａに支承されたクランクシャフト１ｂにクランクロータ３８が軸着され、このクランクロータ３８の外周に、電磁ピックアップなどからなるクランク角センサ３９が対設されている。さらに、上記エンジン１のカムシャフト１ｃに連設するカムロータ４０に、電磁ピックアップなどからなる気筒判別用のカム角センサ４１が対設されている。尚、上記クランク角センサ３９及び上記カム角センサ４１は、電磁ピックアップなどの磁気センサに限らず、光センサなどでも良い。
【００２６】
上記クランクロータ３８は、図２４に示すように、その外周に突起３８ａ，３８ｂ，３８ｃが形成され、これらの各突起３８ａ，３８ｂ，３８ｃが、各気筒（＃１，＃２と＃３，＃４）の圧縮上死点前（ＢＴＤＣ）θ１，θ２，θ３ の位置に形成されており、本実施例においては、θ１ ＝９７゜ＣＡ、θ２ ＝６５゜ＣＡ、θ３ ＝１０゜ＣＡである。
【００２７】
上記クランクロータ３８の各突起は、上記クランク角センサ３９によって検出され、ＢＴＤＣ９７゜，６５゜，１０゜のクランクパルスがエンジン１／２回転毎（１８０゜ＣＡ毎）に出力される。そして、各信号の入力間隔時間がタイマによって計時され、エンジン回転数Ｎｅが算出される。
【００２８】
また、図２５に示すように、上記カムロータ４０の外周には、気筒判別用の突起４０ａ，４０ｂ，４０ｃが形成され、突起４０ａが＃３，＃４気筒の圧縮上死点後（ＡＴＤＣ）θ４ の位置に形成され、突起４０ｂが３個の突起で構成されて最初の突起が＃１気筒のＡＴＤＣθ５ の位置に形成されている。さらに、突起４０ｃが２個の突起で形成され、最初の突起が＃２気筒のＡＴＤＣθ６ の位置に形成されている。本実施例においては、θ４ ＝２０゜ＣＡ、θ５ ＝５゜ＣＡ、θ６ ＝２０゜ＣＡである。
【００２９】
そして、上記カムロータ４０の各突起が上記カム角センサ４１によって検出され、各気筒の燃焼行程順を＃１→＃３→＃２→＃４とした場合、この燃焼行程順と、上記カム角センサ４１からのカムパルスをカウンタによって計数した値とのパターンに基づいて、気筒判別がなされる。
【００３０】
一方、図２６において、符号５０は、エンジン系を制御する電子制御装置（ＥＣＵ）５０であり、このＥＣＵ５０は、燃料噴射制御、点火時期制御などを行なうメインコンピュータ５１と、ノック検出処理を行なう専用のサブコンピュータ５２との２つのコンピュータを中心として構成され、各部に所定の安定化電源を供給する定電圧回路５３や各種の周辺回路が組込まれている。
【００３１】
上記定電圧回路５３は、ＥＣＵリレー５４のリレー接点を介してバッテリ５５に接続されており、このバッテリ５５に、上記ＥＣＵリレー５４のリレーコイルがイグニッションスイッチ５６を介して接続され、上記イグニッションスイッチ５６がＯＮされて上記ＥＣＵリレー５４のリレー接点が閉となったときには制御用電源を上記各コンピュータ５１，５２に供給すると共に、上記定電圧回路５３は上記バッテリ５５に直接接続されており、イグニッションスイッチ５６のＯＮ，ＯＦＦに拘らずバックアップＲＡＭ６１にバックアップ用の電源を供給する。
【００３２】
尚、上記バッテリ５５には、燃料ポンプリレー５７のリレー接点を介して燃料ポンプ２９が接続されている。
【００３３】
上記メインコンピュータ５１は、ＣＰＵ５８、ＲＯＭ５９、ＲＡＭ６０、バックアップＲＡＭ６１、カウンタ・タイマ群６２、シリアル通信インターフェースであるＳＣＩ６３、及び、Ｉ／Ｏインターフェース６４がバスライン６５を介して接続されたマイクロコンピュータであり、上記バックアップＲＡＭ６１には、上記イグニッションスイッチ５６のＯＮ／ＯＦＦに拘らず、バッテリ５５に直接接続する上記定電圧回路５３からバックアップ電源が常時供給されてデータが保持される。
【００３４】
尚、上記カウンタ・タイマ群６２は、フリーランカウンタ、カム角センサ信号の入力計数用カウンタなどの各種カウンタ、燃料噴射タイマ、点火タイマ、定期割込みを発生させるための定期割込みタイマ、クランク角センサ信号の入力間隔計時用タイマ、及び、システム異常監視用のウオッチドッグタイマなどの各種タイマを便宜上総称するものであり、上記メインコンピュータ５１においては、その他、各種のソフトウエアカウンタ・タイマが用いられる。
【００３５】
また、上記サブコンピュータ５２も、上記メインコンピュータ５１と同様、ＣＰＵ７１、ＲＯＭ７２、ＲＡＭ７３、カウンタ・タイマ群７４、ＳＣＩ７５、及び、Ｉ／Ｏインターフェース７６がバスライン７７を介して接続されたマイクロコンピュータであり、上記メインコンピュータ５１とサブコンピュータ５２とは、上記ＳＣＩ６３，７５を介してシリアル通信ラインにより互いに接続されている。
【００３６】
上記メインコンピュータ５１のＩ／Ｏインターフェース６４には、入力ポートに、スロットル開度センサ３３ａ、水温センサ３６、Ｏ２センサ３７、吸気管圧力センサ２２、大気圧センサ４４、車速センサ４２、吸気温センサ４６、及び、バッテリ５５が、Ａ／Ｄ変換器６６を介して接続されるとともに、アイドルスイッチ３３ｂ、スタータスイッチ４３、及びクランク角センサ３９、カム角センサ４１などが接続され、さらに、図示しない各種センサ及びスイッチ類が接続されている。
【００３７】
また、上記Ｉ／Ｏインターフェース６４の出力ポートには、イグナイタ２７が接続されるとともに、ＩＳＣ弁１６、インジェクタ２５、燃料ポンプリレー５７のリレーコイル、ウエストゲート弁制御用デューティソレノイド弁２１が駆動回路６７を介して接続されており、さらに、図示しない各種のアクチュエータ類が接続されている。
【００３８】
一方、上記サブコンピュータ５２のＩ／Ｏインターフェース７６は、入力ポートに、クランク角センサ３９、カム角センサ４１が接続されるとともに、Ａ／Ｄ変換器７８、周波数フィルタ７９、アンプ８０を介してノックセンサ３４が接続されており、上記ノックセンサ３４からのノック検出信号が上記アンプ８０で所定のレベルに増幅された後に上記周波数フィルタ７９により必要な周波数成分が抽出され、上記Ａ／Ｄ変換器７８にてデジタル信号に変換されて入力されるようになっている。
【００３９】
上記メインコンピュータ５１では、各センサ類からの検出信号を処理し、燃料噴射量制御、点火時期制御、アイドル制御などを行い、一方、上記サブコンピュータ５２では、エンジン回転数とエンジン負荷とに基づいてノックセンサ３４からの信号のサンプル区間を設定し、このサンプル区間でノックセンサ３４からの信号を高速にＡ／Ｄ変換して振動波形を忠実にデジタルデータに変換し、このデータに基づきノック発生の有無を判定する。
【００４０】
上記サブコンピュータ５２のＩ／Ｏインターフェース７６の出力ポートは、上記メインコンピュータ５１のＩ／Ｏインターフェース６４の入力ポートに接続されており、上記サブコンピュータ５２でのノック判定結果がＩ／Ｏインターフェース７６に出力される。そして、上記メインコンピュータ５１では、上記サブコンピュータ５２からノック発生有りの判定結果が出力されると、ＳＣＩ６３を介してシリアル通信ラインよりサブコンピュータ５２からノックデータを読込み、このノックデータに基づいて直ちに該当気筒の点火時期を遅らせ、ノックを回避する。
【００４１】
このようなエンジン制御において、上記メインコンピュータ５１では、センサ・スイッチ類からの信号入力処理、燃料噴射制御、点火時期制御、アイドル制御に係わる各種ジョブが、一つのオペレーティングシステム（ＯＳ）の管理下で効率的に実行される。このＯＳは、車輌制御のための各種マネジメント機能、及び、このマネジメント機能に密着した内部ストラテジーを有し、各種ジョブを体系的に結合し、等時間間隔処理により各種ジョブを効率的に実行するようになっている。
【００４２】
以下、上記メインコンピュータ５１による燃料噴射制御について、図１〜図９のルーチンに従って説明する。尚、サブコンピュータ５２はノック検出処理専用のコンピュータであるため、その動作説明を省略する。
【００４３】
本実施例では、ＯＳ側で処理したＡ／Ｄ変換結果、クランク位置情報、エンジン回転数等に基づき、ユーザ側ジョブで、燃料噴射量、点火時期等を設定すると共に、これらの指示値がＯＳによって噴射タイマ、点火タイマにセットされる。そして、１０ｍｓ毎ジョブ，５０ｍｓ毎ジョブ等の定期割込み要求を出力し、またクランク角センサ３９からの信号入力毎（ＢＴＤＣ９７゜，６５゜，１０゜ＣＡ毎のエンジン１回転に６回）に起動されるクランクパルス入力毎ジョブの割込みを許可する。
【００４４】
１０ｍｓ毎ジョブでは、図１〜図３に示す燃料噴射有効パルス幅及び燃料噴射無駄パルス幅設定ルーチンと、図５に示す吸入空気量設定ルーチンが実行される。また、５０ｍｓ毎ジョブでは、図６に示す係数設定ルーチン、及び図７〜図９に示す要求当量比設定ルーチンが実行される。
【００４５】
尚、以下の説明では、図２２に示すように、吸気系のスロットル弁５ａの下流から吸気弁上流までを吸気チャンバ６Ａと総称する。従って、吸気チャンバ６Ａは、スロットル通路５、エアーチャンバ４、インテークマニホルド３、吸気ポート２ａの総称となる。
【００４６】
ここで、以下の説明で使用する用語について簡単に説明する。Ｑは質量流量［ｍｇ／ｓｅｃ］、Ｇは１サイクル当りの質量［ｍｇ／ｃｙｃｌｅ］を意味する。ただし、後述する筒内吸入空気質量Ｇａｉｒ、筒内吸入燃料質量Ｇｆｕｅｌ等は、吸気行程時にのみ吸入されるので、実質的には、１吸気行程当りの質量と同じ意味になる。また、Ｍは吸気チャンバ６Ａの空気質量、Ｐは吸気チャンバ６Ａの吸気管内圧力、Ｇｉｎｊは燃料噴射質量、Ｑｃは筒内吸入空気質量流量、Ｇａｉｒは筒内吸入空気質量、Ｇｆｕｅｌは筒内吸入燃料質量を示す。
【００４７】
まず、図１〜図３に示す燃料噴射有効パルス幅及び燃料噴射無駄パルス幅設定ルーチンを説明する前に、このルーチンで取入れる各パラメータを設定するルーチンについて説明する。
【００４８】
図５の１０ｍｓ毎ジョブで実行されるルーチンでは、１吸気行程当りの筒内吸入空気質量Ｇａｉｒを算出する。まず、ステップＳ４１で、吸気管圧力センサ２２の検出値に基づいて算出したスロットル弁５ａの下流側の吸気管圧力Ｐと、吸気温センサ４６の検出値に基づいて算出した吸気温度Ｔとを読込み、ステップＳ４２で、現在の吸気温度Ｔが５０℃以下、すなわち通常の吸気温度領域にあるか判断し、Ｔ＜５０℃のときには、ステップＳ４３へ進み、現在の吸気温度Ｔで、吸気温度設定値Ｔ’を設定してステップＳ４５へ進む。また、Ｔ≧５０℃のときには、ステップＳ４４で、上記吸気温度設定値Ｔ’を固定値である５０℃に設定した後、ステップＳ４５へ進む。
【００４９】
ところで、筒内吸入空気質量Ｇａｉｒは、上記吸気管圧力Ｐが一定であっても吸気温度Ｔに影響されて変動するため、空気密度σから筒内吸入空気質量Ｇａｉｒを算出する必要がある。筒内の空気密度σが吸気行程時における吸気チャンバ６Ａ内の空気密度と同じであると考えれば、上記筒内の空気密度αは、空気の状態方程式により以下の式から求めることができる。
【００５０】
σ＝Ｐ×（１−Ｒｅｇｒ）／（Ｒ×Ｔ） …（１）
Ｒｅｇｒ；ＥＧＲ率
Ｒ；ガス定数
図１０に筒内吸入空気質量Ｇａｉｒと吸気温度Ｔとの関係を示す。実線は実測により得た吸入空気質量（実測吸気量）の特性、二点鎖線は上記（１）式により算出した空気密度σに従い求めた吸入空気質量（密度比例吸気量）、一点鎖線は吸気温度Ｔをある温度に固定して算出した空気密度σに従い求めた吸入空気質量（圧力比例吸気量）を示す。図に示すように、吸気温度Ｔが５０℃以下の通常の吸気温度領域では、密度比例吸気量が実測吸気量にほぼ沿う特性を示す。しかし、ターボチャージャ付エンジンでの実測結果によると、吸気温度Ｔが５０℃を越えた領域では、上記筒内吸入空気質量Ｇａｉｒは空気密度σに比例して吸込まれておらず、密度比例吸気量は、実測吸気量に比し過小になり、結果として、空燃比がリーン化してしまう。ところで、本来なら、吸気温度Ｔが５０℃以上の領域では、温度毎の補正係数を用いて補正すべきであるが、制御が複雑化するのと、圧力比例吸気量がほぼ実測吸気量に沿う特性を示しているので、本ルーチンでは吸気温度設定値Ｔ’を５０℃に固定して、筒内吸入空気質量Ｇａｉｒを求めるようにした。尚、上記５０℃は実験により求めた臨界温度であり、採用するエンジン形式により相違するため、個別に設定する必要がある。また、本実施例のようにＥＧＲ装置が搭載されていないエンジンでは、上記ＥＧＲ率はゼロである。
【００５１】
そして、ステップＳ４５へ進むと、上記（１）式の吸気温度Ｔを上記吸気温度設定値Ｔ’に代えて空気密度σを算出し、ステップＳ４６で、上記空気密度σと行程容積Ｄとに基づき次式から１吸気行程当りの理論吸気質量Ｇｔｈを算出する。
【００５２】
Ｇｔｈ←Ｄ×σ …（２）
従って、この理論吸気質量Ｇｔｈは、吸気温度設定値Ｔ’を考慮しなければ、上記吸気管圧力Ｐに比例する。
【００５３】
次いで、ステップＳ４７で、上記理論吸気質量Ｇｔｈを基本として筒内吸入空気質量Ｇａｉｒを次式から算出して、ルーチンを抜ける。
【００５４】
Ｇａｉｒ←（Ｇｔｈ−ηｂ）×ηｖ×Ｋｔｒｍ …（３）
ηｂ：吸気損失質量
ηｖ：体積効率
Ｋｔｒｍ：吸入空気量誤差補正係数
上記理論吸気質量Ｇｔｈと筒内吸入空気質量Ｇａｉｒとの関係は、図１１に示す通りであり、理論吸気質量Ｇｔｈがゼロ、すなわち吸気管圧力Ｐが完全な真空になる前に、上記筒内吸入空気質量Ｇａｉｒがゼロになるような、横軸接点ηｂ、傾きηｖの一次式で表される。但し、実際には、筒内吸入空気質量Ｇａｉｒがゼロになるような運転条件は存在しない。また、１吸気行程当りの吸気損失質量ηｂと体積効率ηｖとは、後述するように、エンジン回転数Ｎｅに基づき一次元マップを補間計算付で参照して設定する。さらに、上記吸入空気量誤差補正係数Ｋｔｒｍは、上記理論吸気質量Ｇｔｈと筒内吸入空気質量Ｇａｉｒとの関係が完全には一次式に乗らないこともあるため、運転条件による空気量計測の定常的な誤差分を補正するためのもので、例えば、スロットル開度αとエンジン回転数Ｎｅとに基づき、二次元マップを補間計算付で参照して設定する。図１２に、この吸入空気量誤差補正係数Ｋｔｒｍを設定する二次元マップの特性を例示する。この図に示すように、この吸入空気量誤差補正係数Ｋｔｒｍの値は基本的には、１．０である。
【００５５】
次に、図６に示す５０ｍｓ毎ジョブで実行される係数設定ルーチンについて説明する。このルーチンでは、上記吸気損失質量ηｂ、及び上記体積効率ηｖを設定する。
【００５６】
まず、ステップＳ５１でエンジン回転数Ｎｅを読込み、ステップＳ５２、Ｓ５３で、このエンジン回転数Ｎｅに基づき一次元マップを補間計算付で参照して、体積効率ηｖと吸気損失質量ηｂとを各々設定してルーチンを抜ける。
【００５７】
上記体積効率ηｖ、上記吸気損失質量ηｂの値は、理論的には、１．０であるが、エンジン回転数毎にカムの同調などの影響によって変化する値であり、また上記体積効率ηｖには、高地補正、排気圧補正等の因子を加味させることができる。
【００５８】
このように、本発明では、吸気管圧力Ｐを基本に筒内吸入空気質量Ｇａｉｒを簡単な計算式に基づいて算出し、また、上記体積効率ηｖ、上記吸気損失質量ηｂは、一次元マップから設定することができるため、コンピュータに掛かる負荷が軽くなり、しかも、上記吸気管圧力Ｐさえ正確に計測されれば、始動時を含む全ての領域で適用することができる。
【００５９】
次に、運転条件による筒内混合気の要求当量比設定ルーチンについて、図７〜図９のルーチンに従って説明する。このルーチンは５０ｍｓ毎のジョブで実行され、定常的な要求増量である可燃限界当量比φｔｗと、最大出力当量比φｆｕｌと、排気ガス温度限界当量比φｔｅｘと、触媒浄化率最良当量比φｇａｓと、燃費率最良当量比φｅｃｏとを設定する。これら各要求当量比φｔｗ，φｆｕｌ，φｔｅｘ，φｇａｓ，φｅｃｏは定常的な要求増量であるため、過渡的な燃料遅れ等によって生じる空燃比変動をフィードフォワード的に補正する係数（エアコン増量補正係数、Ｄレンジ増量補正係数等）、或はセンサやインジェクタの誤差を補正する係数（混合比補正係数等）は含まない。また、このルーチンで設定する各要求当量比は、後述する図１のルーチン実行時に比較され、その最大値が目標当量比（理論燃空比Ｆ／Ａに対する増量係数）ＣＯＥＦとして採用される。尚、この当量比φは、空気過剰率（（Ａ／Ｆ）実際／（Ａ／Ｆ）理論）λに対して、φ＝１／λの関係にある。
【００６０】
まず、ステップＳ６１で冷却水温Ｔｗとエンジン回転数Ｎｅとを読込む。そして、ステップＳ６２〜Ｓ６５或はＳ６６で可燃限界当量比φｔｗを設定する。ステップＳ６２では、上記冷却水温Ｔｗとエンジン回転数Ｎｅとに基づき二次元マップを補間計算付で参照して、可燃限界当量比φｔｗを設定する際の基準となる基準可燃限界当量比Ｍｔｗを設定する。この基準可燃限界当量比Ｍｔｗは、冷却水温Ｔｗとの関係によるリーン側の可燃限界を示すもので、図１３に、この基準可燃限界当量比Ｍｔｗを設定する二次元マップの特性を例示する。この図１３に示すように、冷却水温Ｔｗ側の軸では、従来の水温増量に近い設定がなされており、また、ある程度の高水温では、上記基準可燃限界当量比Ｍｔｗはリーン当量比に設定される。また、エンジン回転数Ｎｅ側の軸では、クランキング回転数を含む全ての回転数領域を検索できるように設定されている。筒内吸入空気質量Ｇａｉｒは始動時においても設定されることから（図５のフローチャート参照）、このときの可燃限界当量比φｔｗが始動時の噴射パルス幅を決定する要因になる。また、始動時において、上記基準可燃限界当量比Ｍｔｗは大きな値に設定されているが、これが始動後にエンジン回転数の増加に従って小さな値に変化して行く際に、後述するように、ある一次遅れ時定数で上記可燃限界当量比φｔｗの変化を遅らせることで始動後の増量が確保される。
【００６１】
次いで、ステップＳ６３へ進むと、エンジン回転数等からエンスト状態にあるかを判断し、エンスト状態にあると判断した場合には、ステップＳ６６へ分岐し、上記基準可燃限界当量比Ｍｔｗで今回の可燃限界当量比φｔｗを設定してステップＳ６７へジャンプする。一方、エンスト状態にないと判断した場合には、ステップＳ６４へ進み、上記基準可燃限界当量比Ｍｔｗと前回の可燃限界当量比φｔｗとを比較し、Ｍｔｗ≧φｔｗのときは、要求当量比が増量されているので、ステップＳ６６へ分岐し、この基準可燃限界当量比Ｍｔｗで今回の可燃限界当量比φｔｗを設定して、増量要求に対応させた後、ステップＳ６７へジャンプする。一方、Ｍｔｗ＜φｔｗのときは、要求当量比が減少されているため、ステップＳ６５へ進み、次式に示す所定時定数（本ルーチンでは、１６×５０ｍｓ）の加重平均により一次遅れ処理を行い、今回の可燃限界当量比φｔｗを設定する。その結果、始動後増量が実質的に確保される。
【００６２】
φｔｗ←（１５φｔｗ＋Ｍｔｗ）／１６ …（４）
次いで、ステップＳ６７〜Ｓ７０或はＳ７１で上記最大出力当量比φｆｕｌが設定される。まず、ステップＳ６７では、スロットル開度αを読込み、ステップＳ６８で、エンジン回転数Ｎｅに基づき、一次元マップを補間計算付で参照してフル増量基準値αｔｈを設定する。このフル増量基準値αｔｈは、エンジンが最大出力を要求しているかを判断する基準値で、図１４に示すように、エンジン回転数Ｎｅとスロットル開度αとの関係から特性を求めてマップ化したものである。
【００６３】
そして、ステップＳ６９で上記スロットル開度αと上記フル増量基準値αｔｈとを比較し、α≦αｔｈのときは、フル増量条件不成立と判断し、ステップＳ７０へ進み、最大出力当量比φｆｕｌを０として、ステップＳ７２へ進む。一方、α＞αｔｈのときは、フル増量条件成立と判断し、ステップＳ７１へ分岐し、上記最大出力当量比φｆｕｌを１．２、すなわち、当該エンジンの最大出力となる当量比（フル増量）に設定する。
【００６４】
尚、目標当量比ＣＯＥＦは、本ルーチンにより設定する各要求当量比φｔｗ，φｆｕｌ，φｔｅｘ，φｇａｓ，φｅｃｏの中から最大値を選択的に設定するものであるため、上記ステップＳ７０で最大出力当量比φｆｕｌが０に設定されると、この最大出力当量比φｆｕｌが目標当量比ＣＯＥＦの選択対象から外れることになる。また、排気ガス温度の上昇を抑える目的の増量は、排気ガス温度限界当量比φｔｅｘとして別に設定されているため、上記最大出力当量比φｆｕｌはエンジン回転数Ｎｅによらず、理論空燃比Ａ／Ｆ（＝１４．６）に対して、おおよそ１５〜２０％増量となるはずである。
【００６５】
次いで、ステップＳ７２〜Ｓ７９或はＳ８０で排気ガス温度限界当量比φｔｅｘを設定する。まず、ステップＳ７２で、負荷の一例である筒内吸入空気質量Ｇａｉｒとエンジン回転数Ｎｅと基づき二次元マップを補間計算付で参照して基本排気ガス温度限界当量比Ｍｔｅｘを設定する。この基本排気ガス温度限界当量比Ｍｔｅｘは、排気ガス温度の上昇を設計限界以下に抑え、エンジン及び排気系を保護するために燃料冷却を行う増量である。図１５に、この基本排気ガス温度限界当量比Ｍｔｅｘを設定する二次元マップの特性を例示する。図に示すように、高負荷（高Ｇａｉｒ）、高回転で大きな増量を行うように設定されている。尚、エンジン負荷を検出するパラメータとして上記筒内吸入空気質量Ｇａｉｒに代えて吸気管圧力Ｐを採用し、二次元マップは吸気管圧力Ｐとエンジン回転数Ｎｅとの関係による特性に基づいて設定するようにしても良い。
【００６６】
そして、ステップＳ７３へ進むと、上記基本排気ガス温度限界当量比Ｍｔｅｘと、上記最大出力当量比φｆｕｌのフル増量値である１．２とを比較し、Ｍｔｅｘ≧１．２のときはステップＳ７４へ進み、またＭｔｅｘ＜１．２のときはステップＳ７６へ進む。
【００６７】
そして、上記ステップＳ７６へ進むと、今回の平均排気ガス温度限界当量比φｔｅｘＡＶを次式に示す加重平均により、所定時定数（本ルーチンでは、１６×５０ｍｓ）の一次遅れ処理を行い、ステップＳ７８へ進む。
【００６８】
φｔｅｘＡＶ←（１５・φｔｅｘＡＶ＋Ｍｔｅｘ）／１６ …（５）
また、ステップＳ７４へ進むと、上記基本排気ガス温度限界当量比Ｍｔｅｘと前回までの平均排気ガス温度限界当量比φｔｅｘＡＶとを比較し、Ｍｔｅｘ＜φｔｅｘＡＶのとき、すなわち、１．２≦Ｍｔｅｘ＜φｔｅｘＡＶのときは、ステップＳ７５へ進み、またＭｔｅｘ≧φｔｅｘＡＶのときは、ステップＳ７７へ進む。
【００６９】
上記ステップＳ７５では、今回の平均排気ガス温度限界当量比φｔｅｘＡＶを次式に示す加重平均により所定時定数（本ルーチンでは、３２×５０ｍｓ）の一次遅れ処理を行い、ステップＳ７８へ進む。
【００７０】
φｔｅｘＡＶ←（３１・φｔｅｘＡＶ＋Ｍｔｅｘ）／３２ …（６）
一方、ステップＳ７７へ進むと、今回の平均排気ガス温度限界当量比φｔｅｘＡＶを次式に示す加重平均により所定時定数（本ルーチンでは１２８×５０ｍｓ）の一次遅れ処理を行って、ステップＳ７８へ進む。
【００７１】
φｔｅｘＡＶ←（１２７・φｔｅｘＡＶ＋Ｍｔｅｘ）／１２８ …（７）
従って、Ｍｔｅｘ＜１．２のいわゆる要求当量比が薄い場合には、比較的早い時定数で平均排気ガス温度限界当量比φｔｅｘＡＶが設定され、また、Ｍｔｅｘ≧１．２、且つＭｔｅｘ≧φｔｅｘＡＶの要求当量比が濃い場合には、ゆっくりとした時定数で上記平均排気ガス温度限界当量比φｔｅｘＡＶが設定される。
【００７２】
そして、ステップＳ７８へ進むと、上記ステップＳ７５，Ｓ７６或はＳ７７で算出した平均排気ガス温度限界当量比φｔｅｘＡＶと上記基本排気ガス温度限界当量比Ｍｔｅｘとを比較し、Ｍｔｅｘ＜φｔｅｘＡＶのときは、ステップＳ７９で、上記基本排気ガス温度限界当量比Ｍｔｅｘを今回の排気ガス温度限界当量比φｔｅｘとして設定し、また、Ｍｔｅｘ≧φｔｅｘＡＶのときは、ステップＳ８０で、今回の排気ガス温度限界当量比φｔｅｘを平均排気ガス温度限界当量比φｔｅｘＡＶで設定する。
【００７３】
従って、排気ガス温度限界当量比φｔｅｘは平均排気ガス温度限界当量比φｔｅｘＡＶで上限リミットされ、また、基本排気ガス温度限界当量比Ｍｔｅｘで下限リミットされる。その結果、上記排気ガス温度限界当量比φｔｅｘの上限リミット値は、エンジン回転数Ｎｅと筒内吸入空気質量Ｇａｉｒとに応じて、すなわち、上記ステップＳ７５，Ｓ７６或はＳ７７で設定される平均排気ガス温度限界当量比φｔｅｘＡＶに従って、ゆっくり変化したり、あるいは比較的早く変化したりするようになる。その結果、長時間の全開運転において排気ガス温度が実際に設計限界値に近づくまではφｔｅｘ＜φｆｕｌとなり、従来のように、最大出力当量比φｆｕｌの中に、燃料冷却成分を含ませていた場合に比し、出力の増加が期待できる。尚、この排気ガス温度限界当量比φｔｅｘを排気ガス温度を計測或は推測した値に基づいて設定するようにしても良い。
【００７４】
次いで、上記ステップＳ７９或はＳ８０からステップＳ８１へ進み、このステップＳ８１〜Ｓ８５或はＳ８６で触媒浄化率最良当量比φｇａｓが設定される。この触媒浄化率最良当量比φｇａｓは、触媒コンバータ１１として３元触媒が採用されている場合に対象となり、排気ガスを浄化しなければならない条件が成立したときに、φｇａｓ←１．０とし、それ以外では、φｇａｓ←０とする。
【００７５】
ステップＳ８１〜Ｓ８４では、排気ガスを浄化しなければならない運転条件にあるかを判断する。この判断条件は、スロットル開度αとエンジン回転数Ｎｅとで行い、スロットル開度αが、１０゜≦α＜７０゜（ステップＳ８１，Ｓ８２）で、且つエンジン回転数Ｎｅが、８００ｒｐｍ≦Ｎｅ＜３０００ｒｐｍ（ステップＳ８３，Ｓ８４）のとき、触媒浄化率最良当量比φｇａｓを１．０に設定し（ステップＳ８５）、それ以外のときは、φｇａｓを０に設定する（ステップＳ８６）。尚、上記判断条件は一例であり、排気ガス浄化領域はエンジンの特性等に応じて任意に設定でき、また例えばスロットル開度αに代えて、吸気管圧力Ｐや筒内吸入空気質量Ｇａｉｒを採用しても良い。また、上記触媒コンバータ１１としてリーンＮｏｘ触媒が採用されている場合には、触媒浄化率最良当量比φｇａｓは全運転領域でφｇａｓ←０に設定する。
【００７６】
その後、ステップＳ８７で、負荷の一例としての吸気管圧力Ｐとエンジン回転数Ｎｅとに基づいて二次元マップを補間計算付で参照して、燃費率最良当量比φｅｃｏを設定して、ルーチンを抜ける。
【００７７】
図１６に、この燃費率最良当量比φｅｃｏを設定する二次元マップの特性を例示する。図に示すように、この二次元マップの各領域には、当該エンジン運転条件下での最良の燃費率を得ることのできる当量比が予め実験などから求めて格納されている。
【００７８】
ここで、この燃費率最良当量比φｅｃｏと前述の可燃限界当量比φｔｗとの関係について簡単に説明する。
【００７９】
空燃比を可燃限界当量比φｔｗまでリーン化させると一般には最良の燃費を得ることができなくなる。しかし、暖機途中の低水温時などでは空燃比を燃費率最良当量比φｅｃｏまでリーン化するよりも、上記可燃限界当量比φｔｗにより、リーン化を抑制した方が良好な燃焼を得ることができる場合もあるので、上記可燃限界当量比φｔｗとは別に燃費率最良当量比φｅｃｏが設定されている。尚、上記触媒浄化率最良当量比φｇａｓが全運転領域で、φｇａｓ←１．０の場合には、上記燃費率最良当量比φｅｃｏは隠れてしまうため、その意味を失うが、例えば、排気浄化システムがリーンＮｏｘ触媒を採用している場合には、前述のように、触媒浄化率最良当量比φｇａｓが全運転領域でφｇａｓ←０に設定されるため、この燃費率最良当量比φｅｃｏは触媒浄化率最良当量比としての要素も合わせ持つことになる。また、ＥＧＲを用いて理論空燃比制御により燃費向上を行う際には、上記（１），（２）式から明らかなように、理論吸気質量Ｇｔｈには、ＥＧＲ量が除かれた新気量のみであるため、全運転領域で、φｅｃｏ←１．０としたほうが、空燃比制御性が良くなる。また、この燃費率最良当量比φｅｃｏを、吸気管圧力Ｐに代えて筒内吸入空気質量Ｇａｉｒと、エンジン回転数Ｎｅとの関係により作成した二次元マップを参照して設定するようにしても良い。
【００８０】
そして、上記各ルーチンで求めた設定値を用いて、図１〜図３の１０ｍｓ毎ジョブで実行される燃料噴射有効パルス幅及び燃料噴射無駄パルス幅設定ルーチンのステップＳ１〜Ｓ１０で１サイクル当りの燃料噴射質量を算出し、ステップＳ１１〜Ｓ２４で燃料噴射有効パルス幅Ｔｅと燃料噴射無駄パルス幅Ｔｓとを算出する。
【００８１】
まず、ステップＳ１では、上記筒内混合気の要求当量比設定ルーチンで設定した各要求当量比φｔｗ，φｆｕｌ，φｔｅｘ、φｇａｓ，φｅｃｏを比較し、これらの要求値の最大値を目標当量比ＣＯＥＦとして決定する。
【００８２】
次いで、ステップＳ２で筒内への１吸気行程当りの筒内吸入燃料質量Ｇｆｕｅｌを次式から算出する。
【００８３】
Ｇｆｕｅｌ←Ｇａｉｒ×Ｆ／Ａ×ＣＯＥＦ×Ｋｆｂ …（８）
Ｆ／Ａ：理論燃空比
Ｋｆｂ：空燃比フィードバック補正係数
ここで、理論燃空比Ｆ／Ａは理論空燃比Ａ／Ｆの逆数で、理論空燃比が１４．６であれば、理論燃空比は１／１４．６となる。また、この理論空燃比Ａ／Ｆは、燃料を完全燃焼するに必要な最小の空気量と燃料量との比であり、多種燃料に対応したエンジン等では変数となる。
【００８４】
一方、上記空燃比フィードバック補正係数Ｋｆｂの初期値は１．０であり、空燃比フィードバック制御、及び空燃比学習制御を行う外部のストラテジーによって書換えられる。また、理論燃空比Ｆ／Ａに目標当量比ＣＯＥＦを掛けることで、運転条件に応じた筒内混合気の目標空燃比が設定される。従って、リニアＡ／Ｆセンサでフィードバック制御する際に設定される上記空燃比フィードバック補正係数Ｋｆｂの目標値は上記理論燃空比Ｆ／Ａを１とした場合の上記目標当量比ＣＯＥＦである。
【００８５】
ところで、上記筒内吸入燃料質量Ｇｆｕｅｌは、筒内へ１行程で吸込んで欲しい燃料量であり、定常的にはインジェクタ２５からの噴射量と一致するが、過渡時には応答遅れがあるため一致しない。すなわち、過渡的にはインジェクタ２５からの燃料噴射量が急増しても、燃料の一部が吸気ポート内壁に付着してから筒内に吸込まれるので、その分遅れが生じ、筒内に吸込まれる燃料量はゆっくりと増加する。従って、上記筒内吸入燃料質量Ｇｆｕｅｌが急増したときには、インジェクタ２５からの燃料量を、この筒内吸入燃料質量Ｇｆｕｅｌよりも多めにフィードフォワード噴射させて、筒内への吸入燃料量を筒内吸入燃料質量Ｇｆｕｅｌと一致させる必要がある。この追加増量分に相当する過渡時の１吸気行程当りの燃料追加質量Ｇａｃｃを、次のステップＳ３〜Ｓ９で算出する。
【００８６】
まず、ステップＳ３では、スロットル開度αとエンジン回転数Ｎｅとに基づき二次元マップ（図１７参照）を補間計算付で参照して、静的な筒内吸入空気質量Ｇａｉｒに相当する指標値（静的指標値）Ｍａｃｃを設定する。この静的指標Ｍａｃｃは、負荷変動の少ない定常走行で、スロットル開度αとエンジン回転数Ｎｅを一定に保持すれば、最終的には筒内吸入空気質量Ｇａｉｒが一定になるということを前提にしている。従って、走行中の負荷変動が無視できるとすれば、図４のステップＳ３１に示すように、上記静的指標Ｍａｃｃをスロットル開度αに基づき一次元マップから補間計算付で設定することも可能である。
【００８７】
また、過渡的には、上記筒内吸入空気質量Ｇａｉｒにも遅れがあり、この遅れに対応した燃料追加質量Ｇａｃｃを求める必要がある。
【００８８】
そこで、ステップＳ４では、まず、吸気チャンバ６Ａ（図２２参照）内の吸気遅れを想定した所定時定数（運転条件によって異なるが本ルーチンでは４×１０ｍｓに固定している）の一次遅れ処理を、次式に示す加重平均に基づいて行い、動的な筒内吸入空気質量Ｇａｉｒに相当する指標値（動的指標値）Ｓａｃｃを算出する。
【００８９】
Ｓａｃｃ←（３・Ｓａｃｃ＋Ｍａｃｃ）／４ …（９）
次いで、ステップＳ５で、上記動的指標値Ｓａｃｃと、前回のルーチン実行時に算出した燃料付着による遅れ指標値Ｔａｃｃとを比較する。尚、この遅れ指標値Ｔａｃｃは上記筒内吸入燃料質量Ｇｆｕｅｌに相当する。
【００９０】
そして、Ｓａｃｃ＜Ｔａｃｃ、すなわち、前回の筒内吸入燃料質量Ｇｆｕｅｌに相当する値（Ｔａｃｃ）に対して今回の筒内吸入空気質量Ｇａｉｒに相当する値（Ｓａｃｃ）が少ないときには、ステップＳ６へ進み、上記遅れ指標値Ｔａｃｃを上記動的指標値Ｓａｃｃとして、ステップＳ８へ進む。その結果、Ｓａｃｃ＝Ｔａｃｃとなり、後述のステップＳ９では、上記燃料追加質量Ｇａｃｃ＝０となり、燃料が追加増量されない。
【００９１】
一方、Ｓａｃｃ≧Ｔａｃｃ、すなわち、前回の筒内吸入燃料質量Ｇｆｕｅｌに相当する値（Ｔａｃｃ）に対して今回の筒内吸入空気質量Ｇａｉｒに相当する値（Ｓａｃｃ）が増加しているか、或は変化していないときには、ステップＳ７へ分岐して、この動的指標Ｓａｃｃを、次式に示す加重平均により、燃料付着による時定数（本ルーチンでは、８×１０ｍｓ）の一次遅れ処理を行って、今回の燃料付着による遅れ指標値Ｔａｃｃを算出して、ステップＳ８へ進む。
【００９２】
Ｔａｃｃ←（７・Ｔａｃｃ＋Ｓａｃｃ）／８ …（１０）
そして、ステップＳ８へ進むと、エンジン回転数Ｎｅに基づき一次元マップを補間計算付で参照して係数Ｒａｃｃを設定する。
【００９３】
その後、ステップＳ９へ進み、過渡時の１吸気行程当りの燃料追加質量Ｇａｃｃを次式から算出する。
【００９４】
Ｇａｃｃ←Ｒａｃｃ×（Ｓａｃｃ−Ｔａｃｃ） …（１１）
この（１１）式に示すように、今回の筒内吸入空気質量Ｇａｉｒに相当する値（Ｓａｃｃ）と今回の筒内吸入燃料質量Ｇｆｕｅｌに相当する値（Ｔａｃｃ）との差分が、燃料付着による不足分を補填する増量パターンということになる。
【００９５】
このように、この燃料追加質量Ｇａｃｃは、スロットル開度αとエンジン回転数Ｎｅのみから、簡易的に筒内吸入空気質量Ｇａｉｒの変化量だけを遅れなく検出し、この値を利用して算出している。
【００９６】
図１８にスロットル弁５ａを走行中に僅かに開弁させたときの上記各指標値Ｍａｃｃ，Ｓａｃｃ，Ｔａｃｃの特性を示す。静的指標値Ｍａｃｃはスロットル開度αに追従して増加し、一方、動的指標値Ｓａｃｃは上記静的指標値Ｍａｃｃに対して４×１０ｍｓの一次遅れで増加する。さらに、上記遅れ指標値Ｔａｃｃは、この動的指標値Ｓａｃｃに対して８×１０ｍｓの時定数で増加する。この遅れ指標値Ｔａｃｃが筒内吸入燃料質量Ｇｆｕｅｌに対応していると考えれば、付着分の遅れは、上記動的指標値Ｓａｃｃと上記遅れ指標値Ｔａｃｃとの差分（図のハッチングで示す領域）に相当することになり、上記ステップＳ９では、この差分指標値（Ｓａｃｃ−Ｔａｃｃ）に上記係数Ｒａｃｃを掛けて、燃料追加質量Ｇａｃｃを求めている。ところで、図１９に示すように、この差分指標値（Ｓａｃｃ−Ｔａｃｃ）が、従来のスロットル弁５ａの開弁直後に生じるリーンスパイクのパターンに符合することが、実験から明らかになった。従って、この差分指標値（Ｓａｃｃ−Ｔａｃｃ）相応分を加速増量すれば、排気空燃比は過渡時においてもリーン化することなく一定となる。
【００９７】
上記燃料追加質量Ｇａｃｃは加速増量・加速時追加パルスに近似するが、この燃料追加質量Ｇａｃｃはあくまでも筒内への吸入燃料量を筒内吸入燃料質量Ｇｆｕｅｌにするためのフィードフォワード分であり、加速時に空燃比Ａ／Ｆを一時的にリッチにするためのものではない。従って、もし、加速時に空燃比Ａ／Ｆをリッチにする必要があれば、前記目標当量比ＣＯＥＦの選択枝として、前記図７〜図９に示す要求当量比設定ルーチン中に、加速当量比φａｃｃを第６の増量要因として加えることも、当然考えられる。
【００９８】
また、従来、上記燃料追加質量Ｇａｃｃを算出する方法として三角増量法がある。この三角増量法はスロットル弁の動きを基本に増量値が設定されるため応答性はよいが、筒内へ吸入される燃料量の一次遅れを三角形で近似させながら増量するので、上記図１９に示すリーンスパイクに対して正確に符合せず、過渡時の空燃比が部分的にリッチになったりリーンになったりしてしまう。また、理論的な燃料付着モデルは吸入空気量を計測するセンサ自体の応答遅れやノイズがあるため、性能にばらつきが生じ易く、充分な信頼性が得られない。
【００９９】
そして、ステップＳ１０へ進むと、上記筒内吸入燃料質量Ｇｆｕｅｌに上記燃料追加質量Ｇａｃｃを加算して、１サイクル当りの燃料噴射質量Ｇｉｎｊを算出する。ところで、噴射時期が早期に設定されていて、噴射が終了した後でまだ吸気行程になる前のタイミングで上記燃料噴射質量Ｇｉｎｊが急増した場合には、自動的に追加噴射することも可能で、この燃料噴射質量Ｇｉｎｊは、ある吸気行程から１回転半前、すなわち吸気行程終了直後から次の吸気行程開始までの期間にインジェクタ２５から噴射される燃料の総量が算出される。
【０１００】
次に、ステップＳ１１〜Ｓ２４で、気筒別燃料噴射有効パルス幅Ｔｅ１〜Ｔｅ４、及び気筒別燃料噴射無駄パルス幅Ｔｓ１〜Ｔｓ４が算出される（ここで、１，２，３，４は気筒番号を示す）。
【０１０１】
まず、ステップＳ１１で、最大出力当量比φｆｕｌの値を参照し、ステップＳ１２で排気ガス温度限界当量比φｔｅｘの値を参照する。そして、φｆｕｌ＝０のエンジンが最大出力を要求しておらず、しかも、φｔｅｘ＝０の排気ガス温度限界による増量も要求していない場合には、ステップＳ１３へ進み、空燃比の気筒別のばらつきを補正する気筒別噴射量補正係数Ｋｔｎ１〜Ｋｔｎ４を１．０に設定して、ステップＳ１８へ進む。一方、上記ステップＳ１１で、φｆｕｌ＝１．２（フル増量）、あるいは、ステップＳ１２で、φｔｅｘ≠０と判断されている場合には、ステップＳ１４へ分岐し、ステップＳ１４〜Ｓ１７でエンジン回転数Ｎｅに基づき、気筒別の一次元マップを補間計算付で参照して、上記気筒別噴射量補正係数Ｋｔｎ１〜Ｋｔｎ４を各々設定して、ステップＳ１８へ進む。噴射量補正係数Ｋｔｎは燃料噴射量を気筒別に増減量して空燃比を均一に設定するもので、上記各一次元マップは、各気筒の吸気特性に対応する値が格納されている。
【０１０２】
ところで、運転条件によって燃料噴射量を気筒毎に変えなければならない原因は、本来、筒内吸入空気質量Ｇａｉｒが気筒毎に相違するからであり、従って、正確には筒内吸入空気質量Ｇａｉｒを吸気特性を考慮して気筒毎に算出し、この算出結果から筒内吸入燃料質量Ｇｆｕｅｌ及び燃料噴射質量Ｇｉｎｊを気筒毎に算出すべきであるが、ＣＰＵの計算負荷が重くなるため、上記気筒別噴射量補正係数Ｋｔｎ１〜Ｋｔｎ４を用いて、見かけ上、気筒別燃料噴射有効パルス幅Ｔｅ１〜Ｔｅ４［ｍｓ］を補正するようにした。
【０１０３】
また、空燃比の気筒毎のばらつきが特に問題となるのは、高負荷運転時のノッキングである。中低負荷運転での空燃比のばらつきは排気ガスに余り影響せず、エンジン振動の観点からすると全気筒の燃料量が一致している方が却って各気筒の出力が均等になるので望ましい。従って、フル増量が要求されている場合には、全気筒が最大出力が得られるようにする必要があり、また、排気ガス温度限界による増量が要求されている場合には、全気筒の排気系等の高温化を設計限界以内に抑制してエンジンを保護する必要があるため、このような場合、すなわち、気筒間の空燃比が均一であることを前提とする要求増量が設定されている場合のみ、エンジン回転数毎に各気筒の燃料噴射量を適正に増減させる補正係数（気筒別噴射量補正係数）を設定し、各気筒の筒内混合気の空燃比が目標空燃比となるように個別に制御する。
【０１０４】
そして、ステップＳ１８へ進むと、全体としての燃料噴射有効パルス幅Ｔｅ＿ａｌｌを次式から算出する。
【０１０５】
Ｔｅ＿ａｌｌ←Ｋｍｒ×Ｇｉｎｊ×Ｋｃｏｎ …（１２）
Ｋｍｒ：パルス幅誤差補正係数
Ｋｃｏｎ：インジェクタ容量係数［ｓｅｃ／ｇ］ここで、パルス幅誤差補正係数Ｋｍｒは、インジェクタ２５の動的な流量特性の非線型性の補正であり、本ルーチンでは、エンジン回転数Ｎｅと燃料噴射質量Ｇｉｎｊとに基づき、図２０に示す二次元マップを補間計算付で参照して設定する。また、上記インジェクタ容量係数Ｋｃｏｎは、インジェクタ２５の静的質量流量特性の逆数であり、１グラムの燃料を０．１秒間に噴射するインジェクタでは、Ｋｃｏｎ＝０．１［ｓｅｃ／ｇ］である。なお、このインジェクタ容量係数Ｋｃｏｎは、本ルーチンでは固定値としているが、多種燃料に対応可能なＦＦＶエンジン等で、燃料の比重や粘性が変化する場合には変数として用いる。
【０１０６】
次いで、ステップＳ１９〜Ｓ２２で、上記燃料噴射有効パルス幅Ｔｅ＿ａｌｌを上記気筒別噴射量補正係数Ｋｔｎ１〜Ｋｔｎ４で補正して、気筒別燃料噴射有効パルス幅Ｔｅ１〜Ｔｅ４をそれぞれ設定する。
【０１０７】
その後、ステップＳ２３へ進み、バッテリ電圧ＶＢに基づき、一次元マップを補間計算付で参照し燃料噴射無駄パルス幅Ｔｓ＿ａｌｌ［ｍｓ］を設定し、ステップＳ２４で、気筒別燃料噴射無駄パルス幅Ｔｓ１〜Ｔｓ４を、上記燃料噴射無駄パルス幅Ｔｓ＿ａｌｌで設定してルーチンを抜ける。
【０１０８】
ＯＳ側では、上記気筒別燃料噴射有効パルス幅Ｔｅ１〜Ｔｅ４と上記気筒別燃料噴射無駄パルス幅Ｔｓ１〜Ｔｓ４とをそれぞれ加算して、燃料噴射パルス幅を気筒別に算出する。
【０１０９】
インジェクタ２５に与える電圧パルスの幅（燃料噴射パルス幅）と、燃料噴射質量Ｇｉｎｊとの関係は、図２１に示す通りであり、燃料噴射質量Ｇｉｎｊが決定されれば、基本的に、この燃料噴射質量Ｇｉｎｊにインジェクタ容量係数Ｋｃｏｎを掛けることで上記有効燃料噴射パルス幅Ｔｅが求められ、この有効燃料噴射パルス幅Ｔｅに無効噴射パルス幅である燃料噴射無駄パルス幅Ｔｓを加算することで燃料噴射パルス幅が算出される。
【０１１０】
本実施例では、燃料噴射質量Ｇｉｎｊ分の燃料を１回で噴射するか、２回に分割して噴射するかはＯＳ内部で決定されるため、上記気筒別燃料噴射有効パルス幅Ｔｅ１〜Ｔｅ４と上記気筒別燃料噴射無駄パルス幅Ｔｓ１〜Ｔｓ４とを、予め加算せずに別々にＯＳに渡し、このＯＳ内部で噴射設定の直前に加算して、各インジェクタに与える燃料噴射パルス幅を算出している。こうすることで、ＯＳ内部で、燃料噴射質量Ｇｉｎｊ分の燃料を２回に分割して噴射すると決定した場合でも、１／２の燃料噴射有効パルス幅Ｔｅに燃料噴射無駄パルス幅Ｔｓを加算することで、１回当りの燃料噴射パルス幅を簡単に与えることができる。
【０１１１】
このように、本実施例によれば、例えば新しい燃料付着モデルを開発した際は、筒内吸入空気質量Ｇａｉｒを算出する前記（３）式中に演算式を追加すれば良く、またインジェクタの流量が非線型であるならば燃料噴射有効パルス幅Ｔｅ，燃料噴射無駄パルス幅Ｔｓを補正するだけで良いことになる。その結果、ある運転領域での空燃比制御性に不具合が生じた場合に、どの部分を変更すべきかが明確となるばかりか、その変更が他の部分に影響を及ぼし難いため、ベース制御ソフトとしての活用が期待できる。
【０１１２】
【発明の効果】
以上、説明したように第１発明によれば、吸気管圧力さえ正確に計測されれば、少なくとも、この吸気管圧力と吸気温度、及び吸気損失量と体積効率とに基づいて一次式から筒内吸入空気量を簡単に導くことができ、しかも、上記吸気損失量と体積効率とがエンジン回転数に対する特性を求めて一次元マップ化されているので、二次元マップを必要とする従来の吸入空気量の演算過程に比し、コンピュータの負荷を軽減し、しかもメモリ容量の有効利用を図ることができる。
【０１１３】
また、第２発明によれば、上記吸気温度が、上記筒内吸入空気量が吸気温度に依存して変動しない温度として予め設定した臨界温度以上のときには、この吸気温度を固定値で設定して吸気密度を設定するようにしたので、ターボチャージャ付エンジン等で、上記吸気温度が高温になっても筒内吸入空気量を実際の吸入空気量にほぼ等しい値に設定することができて空燃比のリーン化を防止することができるばかりでなく、良好な過渡応答性を得ることができる。
【図面の簡単な説明】
【図１】１０ｍｓ毎ジョブで実行される燃料噴射有効パルス幅及び燃料噴射無駄パルス幅設定ルーチンを示すフローチャート
【図２】１０ｍｓ毎ジョブで実行される燃料噴射有効パルス幅及び燃料噴射無駄パルス幅設定ルーチンを示すフローチャート（続き）
【図３】１０ｍｓ毎ジョブで実行される燃料噴射有効パルス幅及び燃料噴射無駄パルス幅設定ルーチンを示すフローチャート（続き）
【図４】他の態様による燃料噴射有効パルス幅及び燃料噴射無駄パルス幅設定ルーチンの要部を示すフローチャート
【図５】１０ｍｓ毎ジョブで実行される筒内吸入空気質量設定ルーチンを示すフローチャート
【図６】５０ｍｓ毎ジョブで実行される係数設定ルーチンを示すフローチャート
【図７】５０ｍｓ毎ジョブで実行される要求当量比設定ルーチンを示すフローチャート
【図８】５０ｍｓ毎ジョブで実行される要求当量比設定ルーチンを示すフローチャート（続き）
【図９】５０ｍｓ毎ジョブで実行される要求当量比設定ルーチンを示すフローチャート（続き）
【図１０】筒内吸入空気質量と吸気温度との関係を示す特性図
【図１１】理論吸気質量と筒内吸入空気質量との関係を示す特性図
【図１２】吸入空気量誤差補正係数を設定する二次元マップの概念図
【図１３】基準可燃限界当量比を設定する二次元マップの概念図
【図１４】最大出力当量比の増量領域を示す特性図
【図１５】基本排気ガス温度限界当量比を設定する二次元マップの特性図
【図１６】燃費率最良当量比を設定する二次元マップの特性図
【図１７】静的な筒内吸入空気質量に相当する指標値を設定する二次元マップの概念図
【図１８】静的指標値と動的指標値とこの動的指標値に対する遅れ指標値との関係を示す特性図
【図１９】過渡運転時のリーンスパイクと増量パターンとの関係を示す特性図
【図２０】パルス幅誤差補正係数を設定する二次元マップの特性図
【図２１】燃料噴射パルス幅と燃料噴射質量との関係を示す特性図
【図２２】エンジンの吸気系モデルを示す説明図
【図２３】エンジンの全体構成図
【図２４】クランクロータとクランク角センサの正面図
【図２５】カムロータとカム角センサの正面図
【図２６】電子制御装置の回路構成図
【符号の説明】
５ａ…スロットル弁
Ｄ…行程容積
Ｇａｉｒ…筒内吸入空気質量（吸入空気量）
Ｇｆｕｅｌ…筒内吸入燃料質量
Ｇｔｈ…理論吸気質量（理論吸気量）
Ｎｅ…エンジン回転数
Ｐ…吸気管圧力
Ｔ…吸気温度
σ…空気密度
ηｂ…吸気損失質量（吸気損失量）
ηｖ…体積効率
５０℃…臨界温度、固定値
Claims (2)
- 少なくともスロットル弁下流の吸気管圧力と吸気温度とに基づき空気の状態方程式から空気密度を設定し、
この空気密度と行程容積とから１吸気行程当りの筒内への理論吸気量を設定し、
またエンジン回転数に基づき１吸気行程当りの筒内への吸気損失量と体積効率とを各々に対応する一次元マップを参照して設定し、
上記理論吸気量と上記吸気損失量との差分と上記体積効率とから１吸気行程当りの筒内への吸入空気量を設定することを特徴とするエンジンの吸入空気量検出方法。 - 前記吸気温度が、前記筒内への吸入空気量が吸気温度に依存して変動しない温度として予め設定した臨界温度以上のときには、上記吸気温度を固定値で設定して上記吸気密度を設定することを特徴とする請求項１記載のエンジンの吸入空気量検出方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP28097394A JP3593162B2 (ja) | 1994-11-15 | 1994-11-15 | エンジンの吸入空気量検出方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP28097394A JP3593162B2 (ja) | 1994-11-15 | 1994-11-15 | エンジンの吸入空気量検出方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH08144833A JPH08144833A (ja) | 1996-06-04 |
| JP3593162B2 true JP3593162B2 (ja) | 2004-11-24 |
Family
ID=17632483
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP28097394A Expired - Fee Related JP3593162B2 (ja) | 1994-11-15 | 1994-11-15 | エンジンの吸入空気量検出方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3593162B2 (ja) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP6512078B2 (ja) * | 2015-11-19 | 2019-05-15 | 株式会社デンソー | 噴射制御装置および噴射制御システム |
-
1994
- 1994-11-15 JP JP28097394A patent/JP3593162B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JPH08144833A (ja) | 1996-06-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP2855952B2 (ja) | 内燃エンジンのアイドル回転数制御方法 | |
| JP3314294B2 (ja) | 内燃機関の制御装置 | |
| JPH08114142A (ja) | エンジンのアイドル制御方法 | |
| JPH0953487A (ja) | 内燃機関の燃料噴射制御装置 | |
| CN111720224B (zh) | 一种充气效率的修正方法、系统 | |
| JP3323542B2 (ja) | 吸気温度予測方法 | |
| JP3593162B2 (ja) | エンジンの吸入空気量検出方法 | |
| JP2927074B2 (ja) | 内燃機関の空燃比制御装置 | |
| JP2003065111A (ja) | エンジンの燃料噴射制御装置 | |
| JP3691092B2 (ja) | エンジンの空燃比制御方法 | |
| JP3621731B2 (ja) | エンジンの空燃比制御方法 | |
| JP3686437B2 (ja) | エンジンの空燃比制御方法 | |
| JP3963993B2 (ja) | エンジンの空燃比制御装置 | |
| JP4263258B2 (ja) | リーンバーンエンジンの燃料噴射制御装置 | |
| JPH08144834A (ja) | エンジンの吸入空気量検出方法 | |
| JP2696444B2 (ja) | 内燃機関の燃料供給制御装置 | |
| JPH08144835A (ja) | エンジンの吸入空気量検出方法 | |
| JP2543763B2 (ja) | 内燃機関の燃料供給制御装置 | |
| JPH0988786A (ja) | エンジンの点火時期制御装置 | |
| JPH1193782A (ja) | リーンバーンエンジンの蒸発燃料パージ制御装置 | |
| JPH0710048Y2 (ja) | 内燃機関の燃料供給制御装置 | |
| JP2592327B2 (ja) | 内燃機関の燃料供給制御装置 | |
| JP2636298B2 (ja) | 多気筒内燃機関の空燃比制御装置 | |
| JPS61201842A (ja) | 内燃機関のアイドル時稀薄空燃比制御装置 | |
| JPS63268951A (ja) | 内燃機関の燃料供給制御装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| 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: 20040810 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20040827 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313113 |
|
| R360 | Written notification for declining of transfer of rights |
Free format text: JAPANESE INTERMEDIATE CODE: R360 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20070903 Year of fee payment: 3 |
|
| R370 | Written measure of declining of transfer procedure |
Free format text: JAPANESE INTERMEDIATE CODE: R370 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20070903 Year of fee payment: 3 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313113 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20070903 Year of fee payment: 3 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100903 Year of fee payment: 6 |
|
| LAPS | Cancellation because of no payment of annual fees |