(57)【特許請求の範囲】
【請求項１】
理論燃空比を目標当量比で補正して定常時の目標空燃比を設定するエンジンの空燃比制御方法において、
前記目標当量比を、エンジン運転状態に基づいて設定する複数の異なる要求当量比の中の最大値で設定することを特徴とするエンジンの空燃比制御方法。
【請求項２】
前記要求当量比は、少なくともエンジン回転数とエンジン温度とに基づき空燃比リーン側の可燃限界増量係数を設定する可燃限界当量比と、エンジン回転数とエンジン負荷とに基づき燃料冷却用増量係数を設定する排気ガス温度限界当量比とであることを特徴とする請求項１記載のエンジンの空燃比制御方法。
【請求項３】
前記可燃限界当量比は、前記エンジン回転数と前記エンジン温度とに基づいて設定した基準可燃限界当量比と前回の演算時に設定した可燃限界当量比とを比較し、該基準可燃限界当量比が前回の演算時に設定した上記可燃限界当量比以上のときは該基準可燃限界当量で設定され、又上記基準可燃限界当量比が前回の演算時に設定した上記可燃限界当量比未満のときは前回の演算時に設定した該可燃限界当量比に一次遅れ処理を施して設定されることを特徴とする請求項２記載のエンジンの空燃比制御方法。
【請求項４】
前記排気ガス温度限界当量比は、前記エンジン回転数と前記エンジン負荷とに基づいて設定した基準可燃限界当量比と前回の演算時までの平均排気ガス温度限界当量比とを比較し、該基準基準可燃限界当量比が該平均排気ガス温度限界当量比以上のときは該平均排気ガス温度限界当量比で設定され、又上記基準可燃限界当量比が上記平均排気ガス温度限界当量比未満のときは基準可燃限界当量比で設定されることを特徴とする請求項２或いは３記載のエンジンの空燃比制御方法。
【発明の詳細な説明】
【０００１】
【産業上の利用分野】
本発明は、理論燃空比を基準としたときの目標空燃比を運転条件に応じて可変設定するエンジンの空燃比制御方法に関する。
【０００２】
【従来の技術】
一般に、この種の燃料噴射制御では、Ｌジェトロニック方式と、Ｄジェトロニック方式とがあるが、この両者は、単に、吸入空気量センサの計測値に基づいて吸入空気量を算出するか、或はスロットル弁下流の吸気管圧力に基づいて吸入空気量を算出するかの違いだけであり、Ｌジェトロニック方式とＤジェトロニック方式のいずれを採用するエンジンであっても吸入空気量算出後は、共通の演算過程で、目標空燃比に対応する燃料噴射量を求めることは可能である。
【０００３】
Ｌジェトロニック方式を採用するエンジンとしては、例えば、特開平２−５７４５号公報に開示されているものがある。このエンジンでは、吸入空気量センサで検出した吸入空気質量流量[mg/sec2]がスロットル弁を通過する時刻と同時刻に計測されたとみなし、このスロットル弁下流のチャンバ内への吸入空気の入出力関係から各気筒に吸込まれた吸入空気量を算出している。
【０００４】
一方、Ｄジェトロニック方式を採用するエンジンとしては、特開平６−１８５３９１号公報に開示されているものがある。このエンジンでは、筒内への吸入空気量を空気の状態方程式を用いて推定している。すなわち、スロットル弁の上下流の圧力差、すなわち大気圧と上記吸気管圧力との偏差、スロットル有効開口面積、及びスロットル弁上流側の空気密度等に基づきスロットル弁を単位時間（Δｔ）当りに通過する吸気量（スロットル通過空気量）を近似的に算出し、次いで、スロットル弁下流のチャンバ内空気量Ｇｂの圧力変化分に基づき今回チャンバに充填された空気量の変化分を設定する。そして、今回チャンバに充填された空気量分は、筒内へ吸入されなかったものと見なして、単位時間当りの筒内吸入空気量を算出している。
【０００５】
そして、上記筒内吸入空気量とエンジン回転数とにより決定される基本噴射量に、そのときのエンジン状態に応じて適切な空燃比にするための補正係数をかけ合わせて各インジェクタに対する燃料噴射量を設定する。
【０００６】
【発明が解決しようとする課題】
このように、従来は、筒内吸入空気量が求められると、この筒内吸入空気量を基礎として燃料噴射量を算出し、空燃比は基本噴射量を種々の補正係数により補正した値のトータルとして、結果的に制御しているに過ぎない。
【０００７】
そのため、ある運転状態での空燃比制御性に不具合が生じた場合、どの補正項を他の部分には影響を及ぼすことなく変更すれば良好な空燃比制御性が得られるかの基準が不明確である。すなわち、各気筒への吸入空気量が正確に推定されても、この筒内吸入空気量を基礎とした空燃比制御体系が従来の制御体系を踏襲するものであるため、エンジン制御ソフトを部分的に変更することの可能なベース制御ソフトとしての利便性に欠け、トータル的な変更が余儀なくされてしまう。
【０００８】
本発明は、上記事情に鑑みてなされたもので、ある運転領域で空燃比制御性に不具合が生じた場合でも、どの補正項を変更すべきかの基準が明確になり、且つその変更が他の部分に影響を及ぼし難く、ベース制御ソフトとしての利便性を充分に備えたエンジンの空燃比制御方法を提供することを目的としている。
【０００９】
【課題を解決するための手段】
上記目的を達成するため本発明によるエンジンの空燃比制御方法は、理論燃空比を目標当量比で補正して定常時の目標空燃比を設定するエンジンの空燃比制御方法において、前記目標当量比を、エンジン運転状態に基づいて設定する複数の異なる要求当量比の中の最大値で設定することを特徴とする。また、好ましくは、前記要求当量比は、少なくともエンジン回転数とエンジン温度とに基づき空燃比リーン側の可燃限界増量係数を設定する可燃限界当量比と、エンジン回転数とエンジン負荷とに基づき燃料冷却用増量係数を設定する排気ガス温度限界当量比とであることを特徴とする。
【００１０】
【作  用】
本発明では、エンジン運転状態に応じて複数の異なる要求当量比を設定し、次いで、この各要求当量比の中の最大値を目標当量比として設定する。そして、この目標当量比で、予め設定され、或は多種燃料に対応して設定された理論燃空比を補正して定常時の目標空燃比を設定する。従って、この理論燃空比を１とすれば、上記目標空燃比は上記目標当量比で決定される。
【００１１】
また、上記要求当量比として、少なくともエンジン回転数とエンジン温度とに基づき空燃比リーン側の可燃限界増量係数を設定する可燃限界当量比と、エンジン回転数とエンジン負荷とに基づき燃料冷却用増量係数を設定する排気ガス温度限界当量比とを採用することで、可燃限界当量比では、低温始動時等にリッチ当量比に設定し、高温、高回転へ移行するに従って、次第にリーン当量比に設定し、一方、排気ガス温度限界当量比では、低回転、低負荷運転時にはリーン当量比或は０に設定し、高回転、高負荷運転へ移行するに従って、次第にリッチ当量比に設定することで、低回転域では、上記可燃限界当量比が目標当量比として決定され、一方、高回転、高負荷域では排気ガス温度限界当量比が目標当量比として決定される等、エンジン運転状態に応じて最適な空燃比が設定される。
【００１２】
【実施例】
以下、図面に基づき本発明の実施例を説明する。図２５にはエンジンの全体概略図が示されており、図中の符号１はエンジン（図においては水平対向４気筒型エンジンを示す）を示し、シリンダヘッド２の吸気ポート２ａにインテークマニホルド３が連通され、このインテークマニホルド３の上流にエアーチャンバ４を介してスロットル通路５が連通されている。このスロットル通路５の上流側には、吸気管６を介してエアクリーナ７が取付けられ、このエアクリーナ７が吸入空気の取り入れ口であるエアインテークチャンバ８に連通されている。また、上記排気ポート２ｂにエキゾーストマニホルド９を介して排気管１０が連通され、この排気管１０に触媒コンバータ１１が介装されてマフラ１２に連通されている。
【００１３】
一方、上記スロットル通路５にスロットル弁５ａが設けられ、このスロットル通路５の直上流の上記吸気管６にインタークーラ１３が介装され、さらに、上記吸気管６の上記エアクリーナ７の下流側にレゾネータチャンバ１４が介装されている。
【００１４】
また、上記レゾネータチャンバ１４と上記インテークマニホルド３とを連通して上記スロットル弁５ａの上流側と下流側とをバイパスするバイパス通路１５に、アイドル空気量を調整するアイドル制御（ＩＳＣ）弁１６が介装されている。さらに、上記ＩＳＣ弁１６の直下流側に、吸気圧が負圧のとき開弁し、またターボチャージャ１８によって過給されて吸気圧が正圧になったとき閉弁するチェックバルブ１７が介装されている。
【００１５】
上記ターボチャージャ１８は、上記吸気管６の上記レゾネータチャンバ１４の下流側にコンプレッサが介装され、タービンが上記排気管１０に介装されている。さらに、上記ターボチャージャ１８のタービンハウジング流入口には、ウエストゲート弁１９が介装され、このウエストゲート弁１９には、ウエストゲート弁作動用アクチュエータ２０が連設されている。
【００１６】
上記ウエストゲート弁作動用アクチュエータ２０は、ダイヤフラムにより２室に仕切られ、一方がウエストゲート弁制御用デューティソレノイド弁２１に連通される圧力室を形成し、他方が上記ウエストゲート弁１９を閉方向に付勢するスプリングを収納したスプリング室を形成している。
【００１７】
上記ウエストゲート弁制御用デューティソレノイド弁２１は、上記レゾネータチャンバ１４と上記吸気管６の上記ターボチャージャ１８のコンプレッサ下流とを連通する通路に介装されており、後述する電子制御装置５０（ＥＣＵ；図２８参照）から出力される制御信号のデューティ比に応じて、上記レゾネータチャンバ１４側の圧力と上記コンプレッサ下流側の圧力とを調圧し、上記ウエストゲート弁作動用アクチュエータ２０の圧力室に供給する。
【００１８】
すなわち、上記ＥＣＵ５０によって上記ウエストゲート弁制御用デューティソレノイド弁２１を制御し、上記ウエストゲート弁作動用アクチュエータ２０を作動させて上記ウエストゲート弁１９による排気ガスリリーフを調整することにより、上記ターボチャージャ１８による過給圧を制御するようになっている。
【００１９】
また、上記インテークマニホルド３に、吸気管圧力Ｐを絶対圧で検出する吸気管圧力センサ２２が通路２３を介して連通され、さらに、上記インテークマニホルド３の各気筒の各吸気ポート２ａの直上流側にインジェクタ２５が臨まされている。また、上記シリンダヘッド２の各気筒毎に、その先端を燃焼室に露呈する点火プラグ２６ａが取付けられ、この点火プラグ２６ａに各気筒毎に配設された点火コイル２６ｂを介してイグナイタ２７が接続されている。
【００２０】
上記インジェクタ２５には、燃料タンク２８内に設けたインタンク式の燃料ポンプ２９から燃料フィルタ３０を経て燃料が圧送され、プレッシャレギュレータ３１にてインジェクタ２５への燃料圧力が調圧される。
【００２１】
また、上記吸気管６の上記エアクリーナ７の直下流に、ホットワイヤ或はホットフィルム等を用いた熱式の吸入空気量センサ３２が介装され、上記スロットル弁５ａに、スロットル開度センサ３３ａとアイドルスイッチ３３ｂとを内蔵したスロットルセンサ３３が連設されている。さらに、上記エアーチャンバ４に吸気温センサ４６が臨まされている。
【００２２】
また、上記エンジン１のシリンダブロック１ａにノックセンサ３４が取付けられるとともに、このシリンダブロック１ａの左右両バンクを連通する冷却水通路３５に水温センサ３６が臨まされ、上記排気管１０の上記エキゾーストマニホルド９の集合部にＯ２  センサ３７が臨まされている。
【００２３】
また、上記シリンダブロック１ａに支承されたクランクシャフト１ｂにクランクロータ３８が軸着され、このクランクロータ３８の外周に、電磁ピックアップなどからなるクランク角センサ３９が対設されている。さらに、上記エンジン１のカムシャフト１ｃに連設するカムロータ４０に、電磁ピックアップなどからなる気筒判別用のカム角センサ４１が対設されている。尚、上記クランク角センサ３９及び上記カム角センサ４１は、電磁ピックアップなどの磁気センサに限らず、光センサなどでも良い。
【００２４】
上記クランクロータ３８は、図２６に示すように、その外周に突起３８ａ，３８ｂ，３８ｃが形成され、これらの各突起３８ａ，３８ｂ，３８ｃが、各気筒（＃１，＃２と＃３，＃４）の圧縮上死点前（ＢＴＤＣ）θ１，θ２，θ３  の位置に形成されており、本実施例においては、θ１  ＝９７゜ＣＡ、θ２  ＝６５゜ＣＡ、θ３  ＝１０゜ＣＡである。
【００２５】
上記クランクロータ３８の各突起は、上記クランク角センサ３９によって検出され、ＢＴＤＣ９７゜，６５゜，１０゜のクランクパルスがエンジン１／２回転毎（１８０゜ＣＡ毎）に出力される。そして、各信号の入力間隔時間がタイマによって計時され、エンジン回転数Ｎｅが算出される。
【００２６】
また、図２７に示すように、上記カムロータ４０の外周には、気筒判別用の突起４０ａ，４０ｂ，４０ｃが形成され、突起４０ａが＃３，＃４気筒の圧縮上死点後（ＡＴＤＣ）θ４  の位置に形成され、突起４０ｂが３個の突起で構成されて最初の突起が＃１気筒のＡＴＤＣθ５  の位置に形成されている。さらに、突起４０ｃが２個の突起で形成され、最初の突起が＃２気筒のＡＴＤＣθ６  の位置に形成されている。本実施例においては、θ４  ＝２０゜ＣＡ、θ５  ＝５゜ＣＡ、θ６  ＝２０゜ＣＡである。
【００２７】
そして、上記カムロータ４０の各突起が上記カム角センサ４１によって検出され、各気筒の燃焼行程順を＃１→＃３→＃２→＃４とした場合、この燃焼行程順と、上記カム角センサ４１からのカムパルスをカウンタによって計数した値とのパターンに基づいて、気筒判別がなされる。
【００２８】
一方、図２８において、符号５０は、エンジン系を制御する電子制御装置（ＥＣＵ）５０であり、このＥＣＵ５０は、燃料噴射制御、点火時期制御などを行なうメインコンピュータ５１と、ノック検出処理を行なう専用のサブコンピュータ５２との２つのコンピュータを中心として構成され、各部に所定の安定化電源を供給する定電圧回路５３や各種の周辺回路が組込まれている。
【００２９】
上記定電圧回路５３は、ＥＣＵリレー５４のリレー接点を介してバッテリ５５に接続されており、このバッテリ５５に、上記ＥＣＵリレー５４のリレーコイルがイグニッションスイッチ５６を介して接続され、上記イグニッションスイッチ５６がＯＮされて上記ＥＣＵリレー５４のリレー接点が閉となったときには制御用電源を上記各コンピュータ５１，５２に供給すると共に、上記定電圧回路５３は上記バッテリ５５に直接接続されており、イグニッションスイッチ５６のＯＮ，ＯＦＦに拘らず、バックアップＲＡＭ６１にバックアップ用の電源を供給する。
【００３０】
尚、上記バッテリ５５には、燃料ポンプリレー５７のリレー接点を介して燃料ポンプ２９が接続されている。
【００３１】
上記メインコンピュータ５１は、ＣＰＵ５８、ＲＯＭ５９、ＲＡＭ６０、バックアップＲＡＭ６１、カウンタ・タイマ群６２、シリアル通信インターフェースであるＳＣＩ６３、及び、Ｉ／Ｏインターフェース６４がバスライン６５を介して接続されたマイクロコンピュータであり、上記バックアップＲＡＭ６１には、上記イグニッションスイッチ５６のＯＮ／ＯＦＦに拘らず、バッテリ５５に直接接続する上記定電圧回路５３からバックアップ電源が常時供給されてデータが保持される。
【００３２】
尚、上記カウンタ・タイマ群６２は、フリーランカウンタ、カム角センサ信号の入力計数用カウンタなどの各種カウンタ、燃料噴射タイマ、点火タイマ、定期割込みを発生させるための定期割込みタイマ、クランク角センサ信号の入力間隔計時用タイマ、及び、システム異常監視用のウオッチドッグタイマなどの各種タイマを便宜上総称するものであり、上記メインコンピュータ５１においては、その他、各種のソフトウエアカウンタ・タイマが用いられる。
【００３３】
また、上記サブコンピュータ５２も、上記メインコンピュータ５１と同様、ＣＰＵ７１、ＲＯＭ７２、ＲＡＭ７３、カウンタ・タイマ群７４、ＳＣＩ７５、及び、Ｉ／Ｏインターフェース７６がバスライン７７を介して接続されたマイクロコンピュータであり、上記メインコンピュータ５１とサブコンピュータ５２とは、上記ＳＣＩ６３，７５を介してシリアル通信ラインにより互いに接続されている。
【００３４】
上記メインコンピュータ５１のＩ／Ｏインターフェース６４には、入力ポートに、吸入空気量センサ３２、スロットル開度センサ３３ａ、水温センサ３６、Ｏ２センサ３７、吸気管圧力センサ２２、大気圧センサ４４、車速センサ４２、吸気温センサ４６、及び、バッテリ５５が、Ａ／Ｄ変換器６６を介して接続されるとともに、アイドルスイッチ３３ｂ、スタータスイッチ４３、及びクランク角センサ３９、カム角センサ４１などが接続され、さらに、図示しない各種センサ及びスイッチ類が接続されている。
【００３５】
また、上記Ｉ／Ｏインターフェース６４の出力ポートには、イグナイタ２７が接続されるとともに、ＩＳＣ弁１６、インジェクタ２５、燃料ポンプリレー５７のリレーコイル、ウエストゲート弁制御用デューティソレノイド弁２１が駆動回路６７を介して接続されており、さらに、図示しない各種のアクチュエータ類が接続されている。
【００３６】
一方、上記サブコンピュータ５２のＩ／Ｏインターフェース７６は、入力ポートに、クランク角センサ３９、カム角センサ４１が接続されるとともに、Ａ／Ｄ変換器７８、周波数フィルタ７９、アンプ８０を介してノックセンサ３４が接続されており、上記ノックセンサ３４からのノック検出信号が上記アンプ８０で所定のレベルに増幅された後に上記周波数フィルタ７９により必要な周波数成分が抽出され、上記Ａ／Ｄ変換器７８にてデジタル信号に変換されて入力されるようになっている。
【００３７】
上記メインコンピュータ５１では、各センサ類からの検出信号を処理し、燃料噴射量制御、点火時期制御、アイドル制御などを行い、一方、上記サブコンピュータ５２では、エンジン回転数とエンジン負荷とに基づいてノックセンサ３４からの信号のサンプル区間を設定し、このサンプル区間でノックセンサ３４からの信号を高速にＡ／Ｄ変換して振動波形を忠実にデジタルデータに変換し、このデータに基づきノック発生の有無を判定する。
【００３８】
上記サブコンピュータ５２のＩ／Ｏインターフェース７６の出力ポートは、上記メインコンピュータ５１のＩ／Ｏインターフェース６４の入力ポートに接続されており、上記サブコンピュータ５２でのノック判定結果がＩ／Ｏインターフェース７６に出力される。そして、上記メインコンピュータ５１では、上記サブコンピュータ５２からノック発生有りの判定結果が出力されると、ＳＣＩ６３を介してシリアル通信ラインよりサブコンピュータ５２からノックデータを読込み、このノックデータに基づいて直ちに該当気筒の点火時期を遅らせ、ノックを回避する。
【００３９】
このようなエンジン制御において、上記メインコンピュータ５１では、センサ・スイッチ類からの信号入力処理、燃料噴射制御、点火時期制御、アイドル制御に係わる各種ジョブが、一つのオペレーティングシステム（ＯＳ）の管理下で効率的に実行される。このＯＳは、車輌制御のための各種マネジメント機能、及び、このマネジメント機能に密着した内部ストラテジーを有し、各種ジョブを体系的に結合し、等時間間隔処理により各種ジョブを効率的に実行するようになっている。
【００４０】
以下、上記メインコンピュータ５１による燃料噴射制御について、図１〜図９のルーチンに従って説明する。尚、サブコンピュータ５２はノック検出処理専用のコンピュータであるため、その動作説明を省略する。
【００４１】
本実施例では、ＯＳ側で処理したＡ／Ｄ変換結果、クランク位置情報、エンジン回転数等に基づき、ユーザ側ジョブで、燃料噴射量、点火時期等を設定すると共に、これらの指示値がＯＳによって噴射タイマ、点火タイマにセットされる。そして、１０ｍｓ毎ジョブ，５０ｍｓ毎ジョブ等の定期割込み要求を出力し、またクランク角センサ３９からの信号入力毎（ＢＴＤＣ９７゜，６５゜，１０゜ＣＡ毎のエンジン１回転に６回）に起動されるクランクパルス入力毎ジョブの割込みを許可する。
【００４２】
１０ｍｓ毎ジョブでは、図１〜図３に示す燃料噴射有効パルス幅及び燃料噴射無駄パルス幅設定ルーチンと、図５に示す吸入空気量設定ルーチンが実行される。また、５０ｍｓ毎ジョブでは、図６に示す係数設定ルーチン、及び図７〜図９に示す要求当量比設定ルーチンが実行される。
【００４３】
尚、以下の説明では、図２４に示すように、吸気系のスロットル弁５ａの下流から吸気弁上流までを吸気チャンバ６Ａと総称する。従って、吸気チャンバ６Ａは、スロットル通路５、エアーチャンバ４、インテークマニホルド３、吸気ポート２ａの総称となる。
【００４４】
ここで、以下の説明で使用する用語について簡単に説明する。Ｑは質量流量[mg/sec]、Ｇは１サイクル当りの質量[mg/cycle]を意味する。ただし、後述する筒内吸入空気質量Ｇａｉｒ、筒内吸入燃料質量Ｇｆｕｅｌ等は、吸気行程時にのみ吸入されるので、実質的には、１吸気行程当りの質量と同じ意味になる。また、図２４の左から、Ｑａは吸入空気量センサ３２を通過する空気の質量流量（センサ通過空気質量流量）、Ｍは吸気チャンバ６Ａの空気質量、Ｐは吸気チャンバ６Ａの吸気管内圧力、Ｇｉｎｊは燃料噴射質量、Ｑｃは筒内吸入空気質量流量、Ｇａｉｒは筒内吸入空気質量、Ｇｆｕｅｌは筒内吸入燃料質量を示す。
【００４５】
まず、図１〜図３に示す燃料噴射有効パルス幅及び燃料噴射無駄パルス幅設定ルーチンを説明する前に、このルーチンで取入れる各パラメータを設定するルーチンについて説明する。
【００４６】
図５の１０ｍｓ毎ジョブで実行されるルーチンでは、Ｄジェトロニック方式により１吸気行程当りの筒内吸入空気質量Ｇａｉｒを算出する。まず、ステップＳ４１で、吸気管圧力センサ２２の検出値に基づいて算出したスロットル弁５ａの下流側の吸気管圧力Ｐと、吸気温センサ４６の検出値に基づいて算出した吸気温度Ｔとを読込み、ステップＳ４２で、現在の吸気温度Ｔが５０℃以下、すなわち通常の吸気温度領域にあるか判断し、Ｔ＜５０℃のときには、ステップＳ４３へ進み、現在の吸気温度Ｔで、吸気温度設定値Ｔ’を設定してステップＳ４５へ進む。また、Ｔ≧５０℃のときには、ステップＳ４４で、上記吸気温度設定値Ｔ’を固定値である５０℃に設定した後、ステップＳ４５へ進む。
【００４７】
ところで、筒内吸入空気質量Ｇａｉｒは、上記吸気管圧力Ｐが一定であっても吸気温度Ｔに影響されて変動するため、空気密度σから筒内吸入空気質量Ｇａｉｒを算出する必要がある。筒内の空気密度σが吸気行程時における吸気チャンバ６Ａ内の空気密度と同じであると考えれば、上記筒内の空気密度αは、空気の状態方程式により以下の式から求めることができる。
【００４８】
σ＝Ｐ×（１−Ｒｅｇｒ）／（Ｒ×Ｔ）            …（１）
Ｒｅｇｒ；ＥＧＲ率
Ｒ；ガス定数
図１１に筒内吸入空気質量Ｇａｉｒと吸気温度Ｔとの関係を示す。実線は実測により得た吸入空気質量（実測吸気量）の特性、二点鎖線は上記（１）式により算出した空気密度σに従い求めた吸入空気質量（密度比例吸気量）、一点鎖線は吸気温度Ｔをある温度に固定して算出した空気密度σに従い求めた吸入空気質量（圧力比例吸気量）を示す。図に示すように、吸気温度Ｔが５０℃以下の通常の吸気温度領域では、密度比例吸気量が実測吸気量にほぼ沿う特性を示す。しかし、ターボチャージャ付エンジンでの実測結果によると、吸気温度Ｔが５０℃を越えた領域では、上記筒内吸入空気質量Ｇａｉｒは空気密度σに比例して吸込まれておらず、密度比例吸気量は、実測吸気量に比し過小になり、結果として、空燃比がリーン化してしまう。ところで、本来なら、吸気温度Ｔが５０℃以上の領域では、温度毎の補正係数を用いて補正すべきであるが、制御が複雑化するのと、圧力比例吸気量がほぼ実測吸気量に沿う特性を示しているので、本ルーチンでは吸気温度設定値Ｔ’を５０℃に固定して、筒内吸入空気質量Ｇａｉｒを求めるようにした。尚、上記５０℃は実験により求めた臨界温度であり、採用するエンジン形式により相違するため、個別に設定する必要がある。また、本実施例のようにＥＧＲ装置が搭載されていないエンジンでは、上記ＥＧＲ率はゼロである。
【００４９】
そして、ステップＳ４５へ進むと、上記（１）式の吸気温度Ｔを上記吸気温度設定値Ｔ’に代えて空気密度σを算出し、ステップＳ４６で、上記空気密度σと行程容積Ｄとに基づき次式から１吸気行程当りの理論吸気質量Ｇｔｈを算出する。
【００５０】
Ｇｔｈ←Ｄ×σ        …（２）
従って、この理論吸気質量Ｇｔｈは、吸気温度設定値Ｔ’を考慮しなければ、上記吸気管圧力Ｐに比例する。
【００５１】
次いで、ステップＳ４７で、上記理論吸気質量Ｇｔｈを基本として筒内吸入空気質量Ｇａｉｒを次式から算出して、ルーチンを抜ける。
【００５２】
Ｇａｉｒ←（Ｇｔｈ−ηｂ）×ηｖ×Ｋｔｒｍ        …（３）
ηｂ：吸気損失質量
ηｖ：体積効率
Ｋｔｒｍ：吸入空気量誤差補正係数
上記理論吸気質量Ｇｔｈと筒内吸入空気質量Ｇａｉｒとの関係は、図１２に示す通りであり、理論吸気質量Ｇｔｈがゼロ、すなわち吸気管圧力Ｐが完全な真空になる前に、上記筒内吸入空気質量Ｇａｉｒがゼロになるような、横軸接点ηｂ、傾きηｖの一次式で表される。但し、実際には、筒内吸入空気質量Ｇａｉｒがゼロになるような運転条件は存在しない。また、１吸気行程当りの吸気損失質量ηｂと体積効率ηｖとは、後述するように、エンジン回転数Ｎｅに基づき一次元マップを補間計算付で参照して設定する。さらに、上記吸入空気量誤差補正係数Ｋｔｒｍは、上記理論吸気質量Ｇｔｈと筒内吸入空気質量Ｇａｉｒとの関係が完全には一次式に乗らないこともあるため、運転条件による空気量計測の定常的な誤差分を補正するためのもので、例えば、スロットル開度αとエンジン回転数Ｎｅとに基づき、二次元マップを補間計算付で参照して設定する。図１３に、この吸入空気量誤差補正係数Ｋｔｒｍを設定する二次元マップの特性を例示する。この図に示すように、この吸入空気量誤差補正係数Ｋｔｒｍの値は基本的には、１．０である。
【００５３】
次に、図６に示す５０ｍｓ毎ジョブで実行される係数設定ルーチンについて説明する。このルーチンでは、上記吸気損失質量ηｂ、及び上記体積効率ηｖを設定する。
【００５４】
まず、ステップＳ５１でエンジン回転数Ｎｅを読込み、ステップＳ５２、Ｓ５３で、このエンジン回転数Ｎｅに基づき一次元マップを補間計算付で参照して、体積効率ηｖと吸気損失質量ηｂとを各々設定してルーチンを抜ける。
【００５５】
上記体積効率ηｖ、上記吸気損失質量ηｂの値は、理論的には、１．０であるが、エンジン回転数毎にカムの同調などの影響によって変化する値であり、また上記体積効率ηｖには、高地補正、排気圧補正等の因子を加味させることができる。
【００５６】
このように、本ルーチンでは、吸気管圧力Ｐを基本に筒内吸入空気質量Ｇａｉｒを簡単な計算式に基づいて算出し、また、上記体積効率ηｖ、上記吸気損失質量ηｂは、一次元マップから設定することができるため、コンピュータに掛かる負荷が軽くなり、しかも、上記吸気管圧力Ｐさえ正確に計測されれば、始動時を含む全ての領域で適用することができる。
【００５７】
ところで、上記筒内吸入空気質量Ｇａｉｒは、Ｌジェトロニック方式により算出することもできる。このＬジェトロニック方式を利用して筒内吸入空気質量Ｇａｉｒを算出する手順を、図１０の４ｍｓ毎ジョブで実行されるルーチンに従って説明する。尚、上記吸入空気量センサ３２の検出時刻は、スロットル弁５ａまでの吸気管路長、及び流速によっても異なるが、スロットル弁５ａでの検出時刻とほぼ同じと考えて良い。
【００５８】
まず、ステップＳ１０１で、吸入空気量センサ３２からの出力電圧ＶAFMを読込み、次のステップＳ１０２，Ｓ１０３で上記出力電圧ＶAFMの過渡補正をする。すなわち、図１４に示すように、実験によれば、過渡時の吸入空気量の急激な変化に対して、吸入空気量センサ３２の出力電圧ＶAFMの変化は、理想の変化の２／３の電圧までは充分に速い応答を示すが、残りの１／３は、約２００ｍｓの時定数を持つ遅れで応答することが明らかにされている。このときのセンサ通過空気質量流量Ｑａと出力電圧ＶAFMとの関係を、時定数τ(s)の伝達関数で表せば、
ＶAFM＝（２／３）×Ｑａ＋｛１／（３τ(s)＋３）}×Ｑａ        …（４）
となり、応答遅れは、この（４）式の逆関数で補正すれば良く、この補正式は、
Ｑａ＝(１／２)×［３−１／｛（２／３）τ(s)＋１｝］×ＶAFM   …（５）
となる。
【００５９】
上記時定数τ(s)は２００ｍｓであるため、（２／３）τ(s)＋１≒１３０とすれば、上記（５）式は、
0003691092.tif 000002
と単純化される。すなわち、吸入空気量センサ３２の出力電圧ＶAFMの３倍から、この出力電圧ＶAFMを１３０ｍｓの時定数で一次遅れさせた出力電圧ＶAFMを減算し、その値を１／２にすれば、この吸入空気量センサ３２の過渡時の応答遅れを、簡単に補正することができる。
【００６０】
このような考えに基づき、上記ステップＳ１０２では、前回求めた一次遅れ出力電圧ＶAFMavと今回の出力電圧ＶAFMとから、次式に示す所定時定数（本ルーチンでは、１３０ｍｓの近似値として３２×４ｍｓ）の加重平均に基づいて、今回の一次遅れ出力電圧ＶAFMavを算出する。
【００６１】
ＶAFMav←（３１・ＶAFMav＋ＶAFM）／３２               …（７）
次いで、ステップＳ１０３で、今回の出力電圧ＶAFM及び一次遅れ出力電圧ＶAFMavから、過渡補正後の出力電圧（過渡補正出力電圧）Ｖ'AFMを次式から算出する。
【００６２】
Ｖ'AFM←（３・ＶAFM−ＶAFMav）／２                         …（８）
このように、上記吸入空気量センサ３２の過渡応答遅れを簡単な一次式で導き出すことができるため、コンピュータの負荷を軽減させることができる。
【００６３】
そして、ステップＳ１０４で、過渡補正された上記過渡補正出力電圧Ｖ'AFMを、等間隔３２格子のワード補間テーブルを用いてセンサ通過空気質量流量Ｑａに変換し、ステップＳ１０５で、このセンサ通過空気質量流量Ｑａに吸入空気量誤差補正係数Ｋｔｒｍを乗算して、今回のセンサ通過空気質量流量Ｑ'ａを算出する。
【００６４】
この吸入空気量誤差補正係数Ｋｔｒｍは、運転条件による空気量計測の定常的な誤差分を補正する係数で、前述のＤジェトロニック方式と同様に、例えば、スロットル開度αとエンジン回転数Ｎｅとに基づき、二次元マップを補間計算付で参照して設定する（図１３参照）。
【００６５】
次いで、ステップＳ１０６で、エンジンが始動後かをエンジン回転数Ｎｅに基づいて判断し、Ｎｅ＜３００ｒｐｍのときは、始動時と判断してステップＳ１０７へ進み、Ｎｅ≧３００ｒｐｍのときは、始動後と判断してステップＳ１０８へ進む。
【００６６】
ステップＳ１０７へ進むと、始動時の筒内吸入空気質量Ｇａｉｒを初期値設定してルーチンを抜ける。
【００６７】
始動時に、筒内吸入空気質量Ｇａｉｒを初期値設定するのは、以下の理由による。
【００６８】
1)筒内吸入空気質量Ｇａｉｒは、定常的には、空気流量／エンジン回転数で算出されるが、エンジン回転数及び空気流量が共に小さい場合を考えると、計算結果は不定になり易い。
【００６９】
2)後述のステップＳ１０８で説明するように、筒内吸入空気質量Ｇａｉｒは、始動後は履歴に依存して逐次計算されるため、始動後最初のルーチン実行時に初期値が設定されている必要がある。
【００７０】
3)始動時に、ある始動時噴射幅を設定するのではなく、筒内吸入空気質量Ｇａｉｒを初期値設定することで、始動時、始動後に拘らず全ての運転領域で筒内吸入空気質量Ｇａｉｒを利用した演算が可能になる。
【００７１】
そこで、始動時には、筒内吸入空気質量Ｇａｉｒを推定した値で、初期値設定する。この初期値は、例えば、吸気チャンバ６Ａ内の吸気管圧力Ｐを大気圧とみなし、また吸気温度を冷却水温度により代用して次式から算出する。
【００７２】
初期値＝固定値×大気圧／冷却水温（゜Ｋ）            …（９）
一方、始動後と判断されてステップＳ１０８へ進むと、上記ステップＳ１０５で算出したセンサ通過空気質量流量Ｑ'ａに基づき、筒内吸入空気質量Ｇａｉｒを次式に示す吸気チャンバモデルから算出して、ルーチンを抜ける。
0003691092.tif 000003
ＭＴＣＳ：半回転時間
この（１０）式は、以下の考えに基づく。すなわち、吸気チャンバモデルの計算は、次の２式から算出することができる。
【００７３】
0003691092.tif 000004
上記２式からＭとＱｃとを消去すると、次式に示すように、上記センサ通過空気質量流量Ｑａから筒内吸入空気質量Ｇａｉｒを算出することができる。
0003691092.tif 000005
上記筒内吸入空気質量流量Ｑｃは、単位時間当りの吸気行程数が４サイクルエンジンの場合、気筒数×Ｎｅ／２であるため、
Ｑｃ＝（気筒数／２）×Ｎｅ×Ｇａｉｒ
となり、上記（１０−１）式は、
0003691092.tif 000006
となる。ここで、Δｔ×Ｄ／Ｖ＝固定値＝Ｍｔｃｈとおき、また、採用するエンジンが４気筒であれば、上記（１０−２）式は、
Ｇａｉｒ(new)＝Ｑａ×Ｍｔｃｈ＋（１−２Ｎｅ×Ｍｔｃｈ）×Ｇａｉｒ(old)…（１０−３）
となる。尚、本ルーチンでは、計算の便宜上、エンジン回転数Ｎｅ[rps]ではなく半回転数時間（Ｎｅ／２）を用いており、この半回転数時間をＭＴＣＳとおけば、上記（１０−３）式は、
0003691092.tif 000007
となり、上記ステップＳ１０８で示した（１０）式が導き出される。
【００７４】
そして、この（１０）式を演算周期Δｔ毎に行えば、筒内吸入空気質量Ｇａｉｒを履歴に依存して逐次的に算出することができる。
【００７５】
このように、本実施例では、Ｄジェトロニック方式とＬジェトロニック方式との違いは、単に吸気管圧力センサ２２の検出値に基づいて筒内吸入空気質量Ｇａｉｒを算出するか、吸入空気量センサ３２の検出値に基づいて筒内吸入空気質量Ｇａｉｒを算出するか、すなわち、この筒内吸入空気質量Ｇａｉｒの算出方式の違いに過ぎず、以下のルーチンでは、筒内吸入空気質量Ｇａｉｒが算出されていれば、Ｄジェトロニック方式とＬジェトロニック方式とのいずれを採用するエンジンであっても共通となる。
【００７６】
次に、運転条件による筒内混合気の要求当量比設定ルーチンについて、図７〜図９のルーチンに従って説明する。このルーチンは５０ｍｓ毎のジョブで実行され、定常的な要求増量である可燃限界当量比φｔｗと、最大出力当量比φｆｕｌと、排気ガス温度限界当量比φｔｅｘと、触媒浄化率最良当量比φｇａｓと、燃費率最良当量比φｅｃｏとを設定する。これら各要求当量比φｔｗ，φｆｕｌ，φｔｅｘ，φｇａｓ，φｅｃｏは定常的な要求増量であるため、過渡的な燃料遅れ等によって生じる空燃比変動をフィードフォワード的に補正する係数（エアコン増量補正係数、Ｄレンジ増量補正係数等）、或はセンサやインジェクタの誤差を補正する係数（混合比補正係数等）は含まない。また、このルーチンで設定する各要求当量比は、後述する図１のルーチン実行時に比較され、その最大値が目標当量比（理論燃空比Ｆ／Ａに対する増量係数）ＣOEFとして採用される。尚、この当量比φは、空気過剰率（(A/F)実際／(A/F)理論）λに対して、φ＝１／λの関係にある。
【００７７】
まず、ステップＳ６１で冷却水温Ｔｗとエンジン回転数Ｎｅとを読込む。そして、ステップＳ６２〜Ｓ６５或はＳ６６で空燃比リーン側の可燃限界増量係数を設定する可燃限界当量比φｔｗを設定する。ステップＳ６２では、上記冷却水温Ｔｗとエンジン回転数Ｎｅとに基づき二次元マップを補間計算付で参照して、可燃限界当量比φｔｗを設定する際の基準となる基準可燃限界当量比Ｍｔｗを設定する。この基準可燃限界当量比Ｍｔｗは、冷却水温Ｔｗとの関係によるリーン側の可燃限界を示すもので、図１５に、この基準可燃限界当量比Ｍｔｗを設定する二次元マップの特性を例示する。この図１５に示すように、冷却水温Ｔｗ側の軸では、従来の水温増量に近い設定がなされており、また、ある程度の高水温では、上記基準可燃限界当量比Ｍｔｗはリーン当量比に設定される。また、エンジン回転数Ｎｅ側の軸では、クランキング回転数を含む全ての回転数領域を検索できるように設定されている。筒内吸入空気質量Ｇａｉｒは始動時においても設定されることから（図５或は図１０のフローチャート参照）、このときの可燃限界当量比φｔｗが始動時の噴射パルス幅を決定する要因になる。また、始動時において、上記基準可燃限界当量比Ｍｔｗは大きな値に設定されているが、これが始動後にエンジン回転数の増加に従って小さな値に変化して行く際に、後述するように、ある一次遅れ時定数で上記可燃限界当量比φｔｗの変化を遅らせることで始動後の増量が確保される。
【００７８】
次いで、ステップＳ６３へ進むと、エンジン回転数等からエンスト状態にあるかを判断し、エンスト状態にあると判断した場合には、ステップＳ６６へ分岐し、上記基準可燃限界当量比Ｍｔｗで今回の可燃限界当量比φｔｗを設定してステップＳ６７へジャンプする。一方、エンスト状態にないと判断した場合には、ステップＳ６４へ進み、上記基準可燃限界当量比Ｍｔｗと前回の可燃限界当量比φｔｗとを比較し、Ｍｔｗ≧φｔｗのときは、要求当量比が増量されているので、ステップＳ６６へ分岐し、この基準可燃限界当量比Ｍｔｗで今回の可燃限界当量比φｔｗを設定して、増量要求に対応させた後、ステップＳ６７へジャンプする。一方、Ｍｔｗ＜φｔｗのときは、要求当量比が減少されているため、ステップＳ６５へ進み、次式に示す所定時定数（本ルーチンでは、１６×５０ｍｓ）の加重平均により一次遅れ処理を行い、今回の可燃限界当量比φｔｗを設定する。その結果、始動後増量が実質的に確保される。
【００７９】
φｔｗ←（１５φｔｗ＋Ｍｔｗ）／１６    …（１１）
次いで、ステップＳ６７〜Ｓ７０或はＳ７１で上記最大出力当量比φｆｕｌが設定される。まず、ステップＳ６７では、スロットル開度αを読込み、ステップＳ６８で、エンジン回転数Ｎｅに基づき、一次元マップを補間計算付で参照してフル増量基準値αｔｈを設定する。このフル増量基準値αｔｈは、エンジンが最大出力を要求しているかを判断する基準値で、図１６に示すように、エンジン回転数Ｎｅとスロットル開度αとの関係から特性を求めてマップ化したものである。
【００８０】
そして、ステップＳ６９で上記スロットル開度αと上記フル増量基準値αｔｈとを比較し、α≦αｔｈのときは、フル増量条件不成立と判断し、ステップＳ７０へ進み、最大出力当量比φｆｕｌを０として、ステップＳ７２へ進む。一方、α＞αｔｈのときは、フル増量条件成立と判断し、ステップＳ７１へ分岐し、上記最大出力当量比φｆｕｌを１．２、すなわち、当該エンジンの最大出力となる当量比（フル増量）に設定する。
【００８１】
尚、目標当量比ＣOEFは、本ルーチンにより設定する各要求当量比φｔｗ，φｆｕｌ，φｔｅｘ，φｇａｓ，φｅｃｏの中から最大値を選択的に設定するものであるため、上記ステップＳ７０で最大出力当量比φｆｕｌが０に設定されると、この最大出力当量比φｆｕｌが目標当量比ＣOEFの選択対象から外れることになる。また、排気ガス温度の上昇を抑える目的の増量は、排気ガス温度限界当量比φｔｅｘとして別に設定されているため、上記最大出力当量比φｆｕｌはエンジン回転数Ｎｅによらず、理論空燃比Ａ／Ｆ（＝１４．６）に対して、おおよそ１５〜２０％増量となるはずである。
【００８２】
次いで、ステップＳ７２〜Ｓ７９或はＳ８０で燃料冷却用増量係数を設定する排気ガス温度限界当量比φｔｅｘを設定する。まず、ステップＳ７２で、負荷の一例である筒内吸入空気質量Ｇａｉｒとエンジン回転数Ｎeと基づき二次元マップを補間計算付で参照して基本排気ガス温度限界当量比Ｍｔｅｘを設定する。この基本排気ガス温度限界当量比Ｍｔｅｘは、排気ガス温度の上昇を設計限界以下に抑え、エンジン及び排気系を保護するために燃料冷却を行う増量である。図１７に、この基本排気ガス温度限界当量比Ｍｔｅｘを設定する二次元マップの特性を例示する。図に示すように、高負荷（高Ｇａｉｒ）、高回転で大きな増量を行うように設定されている。尚、エンジン負荷を検出するパラメータとして上記筒内吸入空気質量Ｇａｉｒに代えて吸気管圧力Ｐを採用し、二次元マップは吸気管圧力Ｐとエンジン回転数Ｎｅとの関係による特性に基づいて設定するようにしても良い。
【００８３】
そして、ステップＳ７３へ進むと、上記基本排気ガス温度限界当量比Ｍｔｅｘと、上記最大出力当量比φｆｕｌのフル増量値である１．２とを比較し、Ｍｔｅｘ≧１．２のときはステップＳ７４へ進み、またＭｔｅｘ＜１．２のときはステップＳ７６へ進む。
【００８４】
そして、上記ステップＳ７６へ進むと、今回の平均排気ガス温度限界当量比φｔｅｘAVを次式に示す加重平均により、所定時定数（本ルーチンでは、１６×５０ｍｓ）の一次遅れ処理を行い、ステップＳ７８へ進む。
【００８５】
φｔｅｘAV←（１５・φｔｅｘAV＋Ｍｔｅｘ）／１６      …（１２）
また、ステップＳ７４へ進むと、上記基本排気ガス温度限界当量比Ｍｔｅｘと前回までの平均排気ガス温度限界当量比φｔｅｘAVとを比較し、Ｍｔｅｘ＜φｔｅｘAVのとき、すなわち、１．２≦Ｍｔｅｘ＜φｔｅｘAVのときは、ステップＳ７５へ進み、またＭｔｅｘ≧φｔｅｘAVのときは、ステップＳ７７へ進む。
【００８６】
上記ステップＳ７５では、今回の平均排気ガス温度限界当量比φｔｅｘAVを次式に示す加重平均により所定時定数（本ルーチンでは、３２×５０ｍｓ）の一次遅れ処理を行い、ステップＳ７８へ進む。
【００８７】
φｔｅｘAV←（３１・φｔｅｘAV＋Ｍｔｅｘ）／３２      …（１３）
一方、ステップＳ７７へ進むと、今回の平均排気ガス温度限界当量比φｔｅｘAVを次式に示す加重平均により所定時定数（本ルーチンでは１２８×５０ｍｓ）の一次遅れ処理を行って、ステップＳ７８へ進む。
【００８８】
φｔｅｘAV←（１２７・φｔｅｘAV＋Ｍｔｅｘ）／１２８  …（１４）
従って、Ｍｔｅｘ＜１．２のいわゆる要求当量比が薄い場合には、比較的早い時定数で平均排気ガス温度限界当量比φｔｅｘAVが設定され、また、Ｍｔｅｘ≧１．２、且つＭｔｅｘ≧φｔｅｘAVの要求当量比が濃い場合には、ゆっくりとした時定数で上記平均排気ガス温度限界当量比φｔｅｘAVが設定される。
【００８９】
そして、ステップＳ７８へ進むと、上記ステップＳ７５，Ｓ７６或はＳ７７で算出した平均排気ガス温度限界当量比φｔｅｘAVと上記基本排気ガス温度限界当量比Ｍｔｅｘとを比較し、Ｍｔｅｘ＜φｔｅｘAVのときは、ステップＳ７９で、上記基本排気ガス温度限界当量比Ｍｔｅｘを今回の排気ガス温度限界当量比φｔｅｘとして設定し、また、Ｍｔｅｘ≧φｔｅｘAVのときは、ステップＳ８０で、今回の排気ガス温度限界当量比φｔｅｘを平均排気ガス温度限界当量比φｔｅｘAVで設定する。
【００９０】
従って、排気ガス温度限界当量比φｔｅｘは平均排気ガス温度限界当量比φｔｅｘAVで上限リミットされ、また、基本排気ガス温度限界当量比Ｍｔｅｘで下限リミットされる。その結果、上記排気ガス温度限界当量比φｔｅｘの上限リミット値は、エンジン回転数Ｎｅと筒内吸入空気質量Ｇａｉｒとに応じて、すなわち、上記ステップＳ７５，Ｓ７６或はＳ７７で設定される平均排気ガス温度限界当量比φｔｅｘAVに従って、ゆっくり変化したり、あるいは比較的早く変化したりするようになる。その結果、長時間の全開運転において排気ガス温度が実際に設計限界値に近づくまではφｔｅｘ＜φｆｕｌとなり、従来のように、最大出力当量比φｆｕｌの中に、燃料冷却成分を含ませていた場合に比し、出力の増加が期待できる。尚、この排気ガス温度限界当量比φｔｅｘを排気ガス温度を計測或は推測した値に基づいて設定するようにしても良い。
【００９１】
次いで、上記ステップＳ７９或はＳ８０からステップＳ８１へ進み、このステップＳ８１〜Ｓ８５或はＳ８６で触媒浄化率最良当量比φｇａｓが設定される。この触媒浄化率最良当量比φｇａｓは、触媒コンバータ１１として３元触媒が採用されている場合に対象となり、排気ガスを浄化しなければならない条件が成立したときに、φｇａｓ←１．０とし、それ以外では、φｇａｓ←０とする。
【００９２】
ステップＳ８１〜Ｓ８４では、排気ガスを浄化しなければならない運転条件にあるかを判断する。この判断条件は、スロットル開度αとエンジン回転数Ｎｅとで行い、スロットル開度αが、１０゜≦α＜７０゜（ステップＳ８１，Ｓ８２）で、且つエンジン回転数Ｎｅが、８００ｒｐｍ≦Ｎｅ＜３０００ｒｐｍ（ステップＳ８３，Ｓ８４）のとき、触媒浄化率最良当量比φｇａｓを１．０に設定し（ステップＳ８５）、それ以外のときは、φｇａｓを０に設定する（ステップＳ８６）。尚、上記判断条件は一例であり、排気ガス浄化領域はエンジンの特性等に応じて任意に設定でき、また例えばスロットル開度αに代えて、吸気管圧力Ｐや筒内吸入空気質量Ｇａｉｒを採用しても良い。また、上記触媒コンバータ１１としてリーンＮｏｘ触媒が採用されている場合には、触媒浄化率最良当量比φｇａｓは全運転領域でφｇａｓ←０に設定する。
【００９３】
その後、ステップＳ８７で、負荷の一例としての吸気管圧力Ｐとエンジン回転数Ｎeとに基づいて二次元マップを補間計算付で参照して、燃費率最良当量比φｅｃｏを設定して、ルーチンを抜ける。
【００９４】
図１８に、この燃費率最良当量比φｅｃｏを設定する二次元マップの特性を例示する。図に示すように、この二次元マップの各領域には、当該エンジン運転条件下での最良の燃費率を得ることのできる当量比が予め実験などから求めて格納されている。
【００９５】
ここで、この燃費率最良当量比φｅｃｏと前述の可燃限界当量比φｔｗとの関係について簡単に説明する。
【００９６】
空燃比を可燃限界当量比φｔｗまでリーン化させると一般には最良の燃費を得ることができなくなる。しかし、暖機途中の低水温時などでは空燃比を燃費率最良当量比φｅｃｏまでリーン化するよりも、上記可燃限界当量比φｔｗにより、リーン化を抑制した方が良好な燃焼を得ることができる場合もあるので、上記可燃限界当量比φｔｗとは別に燃費率最良当量比φｅｃｏが設定されている。尚、上記触媒浄化率最良当量比φｇａｓが全運転領域で、φｇａｓ←１．０の場合には、上記燃費率最良当量比φｅｃｏは隠れてしまうため、その意味を失うが、例えば、排気浄化システムがリーンＮｏｘ触媒を採用している場合には、前述のように、触媒浄化率最良当量比φｇａｓが全運転領域でφｇａｓ←０に設定されるため、この燃費率最良当量比φｅｃｏは触媒浄化率最良当量比としての要素も合わせ持つことになる。また、ＥＧＲを用いて理論空燃比制御により燃費向上を行う際には、Ｄジェトロニック方式では、上記（１），（２）式から明らかなように、理論吸気質量Ｇｔｈには、ＥＧＲ量が除かれた新気量のみであり、またＬジェトロニック方式では、吸入空気量センサ３２の出力電圧ＶAFMに基づいて算出した筒内吸入空気質量Ｇａｉｒは、結果としてＥＧＲ量が除かれた新気量のみとなるため、全運転領域で、φｅｃｏ←１．０としたほうが、空燃比制御性が良くなる。また、この燃費率最良当量比φｅｃｏを、吸気管圧力Ｐに代えて筒内吸入空気質量Ｇａｉｒと、エンジン回転数Ｎｅとの関係により作成した二次元マップを参照して設定するようにしても良い。
【００９７】
そして、上記各ルーチンで求めた設定値を用いて、図１〜図３の１０ｍｓ毎ジョブで実行される燃料噴射有効パルス幅及び燃料噴射無駄パルス幅設定ルーチンのステップＳ１〜Ｓ１０で１サイクル当りの燃料噴射質量を算出し、ステップＳ１１〜Ｓ２４で燃料噴射有効パルス幅Ｔｅと燃料噴射無駄パルス幅Ｔｓとを算出する。
【００９８】
まず、ステップＳ１では、上記筒内混合気の要求当量比設定ルーチンで設定した各要求当量比φｔｗ，φｆｕｌ，φｔｅｘ、φｇａｓ，φｅｃｏを比較し、これらの要求値の最大値を目標当量比ＣOEFとして決定する。
【００９９】
次いで、ステップＳ２で筒内への１吸気行程当りの筒内吸入燃料質量Ｇｆｕｅｌを次式から算出する。
【０１００】
Ｇｆｕｅｌ←Ｇａｉｒ×Ｆ／Ａ×ＣOEF×Ｋｆｂ  …（１５）
Ｆ／Ａ：理論燃空比
Ｋｆｂ：空燃比フィードバック補正係数
ここで、理論燃空比Ｆ／Ａは理論空燃比Ａ／Ｆの逆数で、理論空燃比が１４．６であれば、理論燃空比は１／１４．６となる。また、この理論空燃比Ａ／Ｆは、燃料を完全燃焼するに必要な最小の空気量と燃料量との比であり、多種燃料に対応したエンジン等では変数となる。
【０１０１】
一方、上記空燃比フィードバック補正係数Ｋｆｂの初期値は１．０であり、空燃比フィードバック制御、及び空燃比学習制御を行う外部のストラテジーによって書換えられる。また、理論燃空比Ｆ／Ａに目標当量比ＣOEFを掛けることで、運転条件に応じた筒内混合気の目標空燃比が設定される。従って、リニアＡ／Ｆセンサでフィードバック制御する際に設定される上記空燃比フィードバック補正係数Ｋｆｂの目標値は上記理論燃空比Ｆ／Ａを１とした場合の上記目標当量比ＣOEFである。
【０１０２】
ところで、上記筒内吸入燃料質量Ｇｆｕｅｌは、筒内へ１行程で吸込んで欲しい燃料量であり、定常的にはインジェクタ２５からの噴射量と一致するが、過渡時には応答遅れがあるため一致しない。すなわち、過渡的にはインジェクタ２５からの燃料噴射量が急増しても、燃料の一部が吸気ポート内壁に付着してから筒内に吸込まれるので、その分遅れが生じ、筒内に吸込まれる燃料量はゆっくりと増加する。従って、上記筒内吸入燃料質量Ｇｆｕｅｌが急増したときには、インジェクタ２５からの燃料量を、この筒内吸入燃料質量Ｇｆｕｅｌよりも多めにフィードフォワード噴射させて、筒内への吸入燃料量を筒内吸入燃料質量Ｇｆｕｅｌと一致させる必要がある。この追加増量分に相当する過渡時の１吸気行程当りの燃料追加質量Ｇａｃｃを、次のステップＳ３〜Ｓ９で算出する。
【０１０３】
まず、ステップＳ３では、スロットル開度αとエンジン回転数Ｎｅとに基づき二次元マップ（図１９参照）を補間計算付で参照して、静的な筒内吸入空気質量Ｇａｉｒに相当する指標値（静的指標値）Ｍａｃｃを設定する。この静的指標Ｍａｃｃは、負荷変動の少ない定常走行で、スロットル開度αとエンジン回転数Ｎｅを一定に保持すれば、最終的には筒内吸入空気質量Ｇａｉｒが一定になるということを前提にしている。従って、走行中の負荷変動が無視できるとすれば、図４のステップＳ３１に示すように、上記静的指標Ｍａｃｃをスロットル開度αに基づき一次元マップから補間計算付で設定することも可能である。
【０１０４】
また、過渡的には、上記筒内吸入空気質量Ｇａｉｒにも遅れがあり、この遅れに対応した燃料追加質量Ｇａｃｃを求める必要がある。
【０１０５】
そこで、ステップＳ４では、まず、吸気チャンバ６Ａ（図２４参照）内の吸気遅れを想定した所定時定数（運転条件によって異なるが本ルーチンでは４×１０ｍｓに固定している）の一次遅れ処理を、次式に示す加重平均に基づいて行い、動的な筒内吸入空気質量Ｇａｉｒに相当する指標値（動的指標値）Ｓａｃｃを算出する。
【０１０６】
Ｓａｃｃ←（３・Ｓａｃｃ＋Ｍａｃｃ）／４    …（１６）
次いで、ステップＳ５で、上記動的指標値Ｓａｃｃと、前回のルーチン実行時に算出した燃料付着による遅れ指標値Ｔａｃｃとを比較する。尚、この遅れ指標値Ｔａｃｃは上記筒内吸入燃料質量Ｇｆｕｅｌに相当する。
【０１０７】
そして、Ｓａｃｃ＜Ｔａｃｃ、すなわち、前回の筒内吸入燃料質量Ｇｆｕｅｌに相当する値（Ｔａｃｃ）に対して今回の筒内吸入空気質量Ｇａｉｒに相当する値（Ｓａｃｃ）が少ないときには、ステップＳ６へ進み、上記遅れ指標値Ｔａｃｃを上記動的指標値Ｓａｃｃとして、ステップＳ８へ進む。その結果、Ｓａｃｃ＝Ｔａｃｃとなり、後述のステップＳ９では、上記燃料追加質量Ｇａｃｃ＝０となり、燃料が追加増量されない。
【０１０８】
一方、Ｓａｃｃ≧Ｔａｃｃ、すなわち、前回の筒内吸入燃料質量Ｇｆｕｅｌに相当する値（Ｔａｃｃ）に対して今回の筒内吸入空気質量Ｇａｉｒに相当する値（Ｓａｃｃ）が増加しているか、或は変化していないときには、ステップＳ７へ分岐して、この動的指標Ｓａｃｃを、次式に示す加重平均により、燃料付着による時定数（本ルーチンでは、８×１０ｍｓ）の一次遅れ処理を行って、今回の燃料付着による遅れ指標値Ｔａｃｃを算出して、ステップＳ８へ進む。
【０１０９】
Ｔａｃｃ←（７・Ｔａｃｃ＋Ｓａｃｃ）／８    …（１７）
そして、ステップＳ８へ進むと、エンジン回転数Ｎｅに基づき一次元マップを補間計算付で参照して係数Ｒａｃｃを設定する。
【０１１０】
その後、ステップＳ９へ進み、過渡時の１吸気行程当りの燃料追加質量Ｇａｃｃを次式から算出する。
【０１１１】
Ｇａｃｃ←Ｒａｃｃ×（Ｓａｃｃ−Ｔａｃｃ）  …（１８）
この（１８）式に示すように、今回の筒内吸入空気質量Ｇａｉｒに相当する値（Ｓａｃｃ）と今回の筒内吸入燃料質量Ｇｆｕｅｌに相当する値（Ｔａｃｃ）との差分が、燃料付着による不足分を補填する増量パターンということになる。
【０１１２】
このように、この燃料追加質量Ｇａｃｃは、スロットル開度αとエンジン回転数Ｎｅのみから、簡易的に筒内吸入空気質量Ｇａｉｒの変化量だけを遅れなく検出し、この値を利用して算出している。
【０１１３】
図２０にスロットル弁５ａを走行中に僅かに開弁させたときの上記各指標値Ｍａｃｃ，Ｓａｃｃ，Ｔａｃｃの特性を示す。静的指標値Ｍａｃｃはスロットル開度αに追従して増加し、一方、動的指標値Ｓａｃｃは上記静的指標値Ｍａｃｃに対して４×１０ｍｓの一次遅れで増加する。さらに、上記遅れ指標値Ｔａｃｃは、この動的指標値Ｓａｃｃに対して８×１０ｍｓの時定数で増加する。この遅れ指標値Ｔａｃｃが筒内吸入燃料質量Ｇｆｕｅｌに対応していると考えれば、付着分の遅れは、上記動的指標値Ｓａｃｃと上記遅れ指標値Ｔａｃｃとの差分（図のハッチングで示す領域）に相当することになり、上記ステップＳ９では、この差分指標値（Ｓａｃｃ−Ｔａｃｃ）に上記係数Ｒａｃｃを掛けて、燃料追加質量Ｇａｃｃを求めている。ところで、図２１に示すように、この差分指標値（Ｓａｃｃ−Ｔａｃｃ）が、従来のスロットル弁５ａの開弁直後に生じるリーンスパイクのパターンに符合することが、実験から明らかになった。従って、この差分指標値（Ｓａｃｃ−Ｔａｃｃ）相応分を加速増量すれば、排気空燃比は過渡時においてもリーン化することなく一定となる。
【０１１４】
上記燃料追加質量Ｇａｃｃは加速増量・加速時追加パルスに近似するが、この燃料追加質量Ｇａｃｃはあくまでも筒内への吸入燃料量を筒内吸入燃料質量Ｇｆｕｅｌにするためのフィードフォワード分であり、加速時に空燃比Ａ／Ｆを一時的にリッチにするためのものではない。従って、もし、加速時に空燃比Ａ／Ｆをリッチにする必要があれば、前記目標当量比ＣOEFの選択枝として、前記図７〜図９に示す要求当量比設定ルーチン中に、加速当量比φａｃｃを第６の増量要因として加えることも、当然考えられる。
【０１１５】
また、従来、上記燃料追加質量Ｇａｃｃを算出する方法として三角増量法がある。この三角増量法はスロットル弁の動きを基本に増量値が設定されるため応答性はよいが、筒内へ吸入される燃料量の一次遅れを三角形で近似させながら増量するので、上記図２１に示すリーンスパイクに対して正確に符合せず、過渡時の空燃比が部分的にリッチになったりリーンになったりしてしまう。また、理論的な燃料付着モデルは吸入空気量を計測するセンサ自体の応答遅れやノイズがあるため、性能にばらつきが生じ易く、充分な信頼性が得られない。
【０１１６】
そして、ステップＳ１０へ進むと、上記筒内吸入燃料質量Ｇｆｕｅｌに上記燃料追加質量Ｇａｃｃを加算して、１サイクル当りの燃料噴射質量Ｇｉｎｊを算出する。ところで、噴射時期が早期に設定されていて、噴射が終了した後でまだ吸気行程になる前のタイミングで上記燃料噴射質量Ｇｉｎｊが急増した場合には、自動的に追加噴射することも可能で、この燃料噴射質量Ｇｉｎｊは、ある吸気行程から１回転半前、すなわち吸気行程終了直後から次の吸気行程開始までの期間にインジェクタ２５から噴射される燃料の総量が算出される。
【０１１７】
次に、ステップＳ１１〜Ｓ２４で、気筒別燃料噴射有効パルス幅Ｔｅ１〜Ｔｅ４、及び気筒別燃料噴射無駄パルス幅Ｔｓ１〜Ｔｓ４が算出される（ここで、１，２，３，４は気筒番号を示す）。
【０１１８】
まず、ステップＳ１１で、最大出力当量比φｆｕｌの値を参照し、ステップＳ１２で排気ガス温度限界当量比φｔｅｘの値を参照する。そして、φｆｕｌ＝０のエンジンが最大出力を要求しておらず、しかも、φｔｅｘ＝０の排気ガス温度限界による増量も要求していない場合には、ステップＳ１３へ進み、空燃比の気筒別のばらつきを補正する気筒別噴射量補正係数Ｋｔｎ１〜Ｋｔｎ４を１．０に設定して、ステップＳ１８へ進む。一方、上記ステップＳ１１で、φｆｕｌ＝１．２（フル増量）、あるいは、ステップＳ１２で、φｔｅｘ≠０と判断されている場合には、ステップＳ１４へ分岐し、ステップＳ１４〜Ｓ１７でエンジン回転数Ｎｅに基づき、気筒別の一次元マップを補間計算付で参照して、上記気筒別噴射量補正係数Ｋｔｎ１〜Ｋｔｎ４を各々設定して、ステップＳ１８へ進む。噴射量補正係数Ｋｔｎは燃料噴射量を気筒別に増減量して空燃比を均一に設定するもので、上記各一次元マップは、各気筒の吸気特性に対応する値が格納されている。
【０１１９】
ところで、運転条件によって燃料噴射量を気筒毎に変えなければならない原因は、本来、筒内吸入空気質量Ｇａｉｒが気筒毎に相違するからであり、従って、正確には筒内吸入空気質量Ｇａｉｒを吸気特性を考慮して気筒毎に算出し、この算出結果から筒内吸入燃料質量Ｇｆｕｅｌ及び燃料噴射質量Ｇｉｎｊを気筒毎に算出すべきであるが、ＣＰＵの計算負荷が重くなるため、上記気筒別噴射量補正係数Ｋｔｎ１〜Ｋｔｎ４を用いて、見かけ上、気筒別燃料噴射有効パルス幅Ｔｅ１〜Ｔｅ４[ms]を補正するようにした。
【０１２０】
また、空燃比の気筒毎のばらつきが特に問題となるのは、高負荷運転時のノッキングである。中低負荷運転での空燃比のばらつきは排気ガスに余り影響せず、エンジン振動の観点からすると全気筒の燃料量が一致している方が却って各気筒の出力が均等になるので望ましい。従って、フル増量が要求されている場合には、全気筒が最大出力が得られるようにする必要があり、また、排気ガス温度限界による増量が要求されている場合には、全気筒の排気系等の高温化を設計限界以内に抑制してエンジンを保護する必要があるため、このような場合、すなわち、気筒間の空燃比が均一であることを前提とする要求増量が設定されている場合のみ、エンジン回転数毎に各気筒の燃料噴射量を適正に増減させる補正係数（気筒別噴射量補正係数）を設定し、各気筒の筒内混合気の空燃比が目標空燃比となるように個別に制御する。
【０１２１】
そして、ステップＳ１８へ進むと、全体としての燃料噴射有効パルス幅Ｔｅ＿ａｌｌを次式から算出する。
【０１２２】
Ｔｅ＿ａｌｌ←Ｋｍｒ×Ｇｉｎｊ×Ｋｃｏｎ    …（１９）
Ｋｍｒ：パルス幅誤差補正係数
Ｋｃｏｎ：インジェクタ容量係数[sec/g]
ここで、パルス幅誤差補正係数Ｋｍｒは、インジェクタ２５の動的な流量特性の非線型性の補正であり、本ルーチンでは、エンジン回転数Ｎｅと燃料噴射質量Ｇｉｎｊとに基づき、図２２に示す二次元マップを補間計算付で参照して設定する。また、上記インジェクタ容量係数Ｋｃｏｎは、インジェクタ２５の静的質量流量特性の逆数であり、１グラムの燃料を０．１秒間に噴射するインジェクタでは、Ｋｃｏｎ＝０．１[sec/g]である。なお、このインジェクタ容量係数Ｋｃｏｎは、本ルーチンでは固定値としているが、多種燃料に対応可能なＦＦＶエンジン等で、燃料の比重や粘性が変化する場合には変数として用いる。
【０１２３】
次いで、ステップＳ１９〜Ｓ２２で、上記燃料噴射有効パルス幅Ｔｅ＿ａｌｌを上記気筒別噴射量補正係数Ｋｔｎ１〜Ｋｔｎ４で補正して、気筒別燃料噴射有効パルス幅Ｔｅ１〜Ｔｅ４をそれぞれ設定する。
【０１２４】
その後、ステップＳ２３へ進み、バッテリ電圧ＶＢに基づき、一次元マップを補間計算付で参照し燃料噴射無駄パルス幅Ｔｓ＿ａｌｌ[ms]を設定し、ステップＳ２４で、気筒別燃料噴射無駄パルス幅Ｔｓ１〜Ｔｓ４を、上記燃料噴射無駄パルス幅Ｔｓ＿ａｌｌで設定してルーチンを抜ける。
【０１２５】
ＯＳ側では、上記気筒別燃料噴射有効パルス幅Ｔｅ１〜Ｔｅ４と上記気筒別燃料噴射無駄パルス幅Ｔｓ１〜Ｔｓ４とをそれぞれ加算して、燃料噴射パルス幅を気筒別に算出する。
【０１２６】
インジェクタ２５に与える電圧パルスの幅（燃料噴射パルス幅）と、燃料噴射質量Ｇｉｎｊとの関係は、図２３に示す通りであり、燃料噴射質量Ｇｉｎｊが決定されれば、基本的に、この燃料噴射質量Ｇｉｎｊにインジェクタ容量係数Ｋｃｏｎを掛けることで上記有効燃料噴射パルス幅Ｔｅが求められ、この有効燃料噴射パルス幅Ｔｅに無効噴射パルス幅である燃料噴射無駄パルス幅Ｔｓを加算することで燃料噴射パルス幅が算出される。
【０１２７】
本実施例では、燃料噴射質量Ｇｉｎｊ分の燃料を１回で噴射するか、２回に分割して噴射するかはＯＳ内部で決定されるため、上記気筒別燃料噴射有効パルス幅Ｔｅ１〜Ｔｅ４と上記気筒別燃料噴射無駄パルス幅Ｔｓ１〜Ｔｓ４とを、予め加算せずに別々にＯＳに渡し、このＯＳ内部で噴射設定の直前に加算して、各インジェクタに与える燃料噴射パルス幅を算出している。こうすることで、ＯＳ内部で、燃料噴射質量Ｇｉｎｊ分の燃料を２回に分割して噴射すると決定した場合でも、１／２の燃料噴射有効パルス幅Ｔｅに燃料噴射無駄パルス幅Ｔｓを加算することで、１回当りの燃料噴射パルス幅を簡単に与えることができる。
【０１２８】
このように、本実施例によれば、例えば新しい燃料付着モデルを開発した際は、筒内吸入空気質量Ｇａｉｒを算出する前記（３）式、或は（１０）式中に演算式を追加すれば良く、またインジェクタの流量が非線型であるならば燃料噴射有効パルス幅Ｔｅ，燃料噴射無駄パルス幅Ｔｓを補正するだけで良いことになる。その結果、ある運転領域での空燃比制御性に不具合が生じた場合に、どの部分を変更すべきかが明確となるばかりか、その変更が他の部分に影響を及ぼし難いため、ベース制御ソフトとしての活用が期待できる。
【０１２９】
【発明の効果】
以上、説明したように本発明によれば、目標当量比を、エンジン運転状態に基づいて設定する複数の異なる要求当量比の中の最大値で設定するようにしたので、例えば、ある運転領域で空燃比制御性に不具合が生じた場合に、その運転特性に応じて、どの補正項を変更すべきかの基準が明確になるばかりでなく、その変更が他の補正項に対して影響を及ぼし難く、ベース制御ソフトとしての利便性が備えられる。
【０１３０】
また、上記要求当量比として、少なくとも可燃限界当量比と、排気ガス温度限界当量比とを備えることで、可燃限界当量比では、低温始動時等にリッチ当量比に設定し、高温、高回転へ移行するに従って、次第にリーン当量比に設定し、一方、排気ガス温度限界当量比では、低回転、低負荷運転時にはリーン当量比或は０に設定し、高回転、高負荷運転へ移行するに従って、次第にリッチ当量比に設定することで、低回転域では、上記可燃限界当量比が目標当量比として設定され、一方、高回転、高負荷域では排気ガス温度限界当量比が目標当量比として設定される等、エンジン運転状態に応じて最適な空燃比を設定することができる。
【図面の簡単な説明】
【図１】１０ｍｓ毎ジョブで実行される燃料噴射有効パルス幅及び燃料噴射無駄パルス幅設定ルーチンを示すフローチャート
【図２】１０ｍｓ毎ジョブで実行される燃料噴射有効パルス幅及び燃料噴射無駄パルス幅設定ルーチンを示すフローチャート（続き）
【図３】１０ｍｓ毎ジョブで実行される燃料噴射有効パルス幅及び燃料噴射無駄パルス幅設定ルーチンを示すフローチャート（続き）
【図４】他の態様による燃料噴射有効パルス幅及び燃料噴射無駄パルス幅設定ルーチンの要部を示すフローチャート
【図５】１０ｍｓ毎ジョブで実行される筒内吸入空気質量設定ルーチンを示すフローチャート
【図６】５０ｍｓ毎ジョブで実行される係数設定ルーチンを示すフローチャート
【図７】５０ｍｓ毎ジョブで実行される要求当量比設定ルーチンを示すフローチャート
【図８】５０ｍｓ毎ジョブで実行される要求当量比設定ルーチンを示すフローチャート（続き）
【図９】５０ｍｓ毎ジョブで実行される要求当量比設定ルーチンを示すフローチャート（続き）
【図１０】他の態様による４ｍｓ毎ジョブ実行される筒内吸入空気質量設定ルーチンを示すフローチャート
【図１１】筒内吸入空気質量と吸気温度との関係を示す特性図
【図１２】理論吸気質量と筒内吸入空気質量との関係を示す特性図
【図１３】吸入空気量誤差補正係数を設定する二次元マップの概念図
【図１４】過渡時における吸入空気量センサの出力電圧を示す特性図
【図１５】基本可燃限界当量比を設定する二次元マップの概念図
【図１６】最大出力当量比の増量領域を示す特性図
【図１７】基本排気ガス温度限界当量比を設定する二次元マップの特性図
【図１８】燃費率最良当量比を設定する二次元マップの特性図
【図１９】静的な筒内吸入空気質量に相当する指標値を設定する二次元マップの概念図
【図２０】静的指標値と動的指標値とこの動的指標値に対する遅れ指標値との関係を示す特性図
【図２１】過渡運転時のリーンスパイクと増量パターンとの関係を示す特性図
【図２２】パルス幅誤差補正係数を設定する二次元マップの特性図
【図２３】燃料噴射パルス幅と燃料噴射質量との関係を示す特性図
【図２４】エンジンの吸気系モデルを示す説明図
【図２５】エンジンの全体構成図
【図２６】クランクロータとクランク角センサの正面図
【図２７】カムロータとカム角センサの正面図
【図２８】電子制御装置の回路構成図
【符号の説明】
Ｆ／Ａ…理論燃空比
ＣOEF…目標当量比
Ｇａｉｒ…筒内吸入空気質量（エンジン負荷）
Ｎｅ…エンジン回転数
Ｐ…吸気管圧力（エンジン負荷）
Ｔｗ…冷却水温（エンジン温度）
φｅｃｏ…燃費率最良当量比（要求当量比）
φｆｕｌ…最大出力当量比（要求当量比）
φｇａｓ…触媒浄化率最良当量比（要求当量比）
φｔｅｘ…排気ガス温度限界当量比（要求当量比）
φｔｗ…可燃限界当量比（要求当量比）
【図１】
0003691092.tif 000008
【図２】
0003691092.tif 000009
【図３】
0003691092.tif 000010
【図４】
0003691092.tif 000011
【図５】
0003691092.tif 000012
【図６】
0003691092.tif 000013
【図７】
0003691092.tif 000014
【図８】
0003691092.tif 000015
【図９】
0003691092.tif 000016
【図１０】
0003691092.tif 000017
【図１１】
0003691092.tif 000018
【図１２】
0003691092.tif 000019
【図１３】
0003691092.tif 000020
【図１４】
0003691092.tif 000021
【図１５】
0003691092.tif 000022
【図１６】
0003691092.tif 000023
【図１７】
0003691092.tif 000024
【図１８】
0003691092.tif 000025
【図１９】
0003691092.tif 000026
【図２０】
0003691092.tif 000027
【図２１】
0003691092.tif 000028
【図２２】
0003691092.tif 000029
【図２３】
0003691092.tif 000030
【図２４】
0003691092.tif 000031
【図２５】
0003691092.tif 000032
【図２６】
0003691092.tif 000033
【図２７】
0003691092.tif 000034
【図２８】
0003691092.tif 000035