JPH0617681A

JPH0617681A - 内燃機関の燃料噴射量制御装置

Info

Publication number: JPH0617681A
Application number: JP4200331A
Authority: JP
Inventors: Hidetaka Maki; 秀隆牧; Yusuke Hasegawa; 祐介長谷川; Naosuke Akasaki; 修介赤崎; Isao Komoriya; 勲小森谷; Toshiaki Hirota; 俊明廣田
Original assignee: Honda Motor Co Ltd
Current assignee: Honda Motor Co Ltd
Priority date: 1992-07-03
Filing date: 1992-07-03
Publication date: 1994-01-25
Anticipated expiration: 2012-12-03
Also published as: JP2683985B2

Abstract

(57)【要約】【構成】壁面付着プラントの伝達特性の逆伝達特性を
持つ補償器をプラントに直列に接続して１つの仮想プラ
ントとする。補償器のパラメータはマップで補正すると
共に、仮想プラントの伝達特性が１及びその付近から逸
脱するとき、その逆伝達特性を持つ様に適応制御器を動
作させる。【効果】吸気圧力などの速い変化にも即応して、かつ
経時変化などにも適応し、筒内実吸入燃料量を常に目標
値に一定させることができる。

Description

【発明の詳細な説明】

【０００１】

【産業上の利用分野】この発明は内燃機関の燃料噴射量
制御装置に関し、より具体的には噴射した燃料が吸気管
などに付着することで生ずる燃料の輸送遅れを適応的に
補正して筒内実吸入燃料量を目標値に常に一致させる様
にした内燃機関の燃料噴射量制御装置に関する。

【０００２】

【従来の技術】機関の過渡運転時には実際の筒内実吸入
燃料量が目標値に一致せず、リーンスパイクやリッチス
パイクが生じることがある。その原因の一つに、燃料が
吸気管の壁面などに付着することによって起こる燃料の
輸送遅れを挙げることができる。この燃料の輸送遅れ
は、機関の運転状態や、初期バラツキや、更には吸気管
へのデポジット付着などによる経時変化によってその挙
動が変化する。そこで、特開平２−１７３３３４号公報
などに記載の様に、適応制御を導入して付着パラメータ
を直接算出または推定し、内燃機関の燃料噴射量制御を
行う技術が提案されている。

【０００３】

【発明が解決しようとする課題】ところで燃料の輸送遅
れをもたらす運転状態には、水温、吸気温などの比較的
変化の遅いもののみならず、吸気圧力などの様に時間的
な変化が速いものも含まれる。例えば低回転時にアクセ
ルペダルを急激に踏み込んだ様なとき吸気圧力は速やか
に立ち上がり、それに応じて燃料の付着状態も速やかに
変化するが、前記した従来技術の場合には付着プラント
の入出力の応答のみを見ているため、この様な速い変化
に十分追従することができなかった。即ち、実際の燃料
の挙動は付着プラントの出力変化になって表れてくる前
に変化しているが、前記した従来技術においてはプラン
ト出力が変化して始めて付着パラメータを算出（推定）
するため、追従性に欠ける不都合があった。

【０００４】従って、この発明の目的は上記した欠点を
解消し、燃料の挙動をリアルタイムに把握し、燃料の輸
送遅れの変化に迅速に対応して筒内実吸入燃料量を目標
値に一致させる様にした内燃機関の燃料噴射量制御装置
を提供することを目的とする。

【０００５】

【課題を解決するための手段】上記の目的を解決するた
めに本発明に係る内燃機関の燃料噴射量制御装置は、内
燃機関の吸気管に付着する燃料量を状態変数として燃料
の挙動をモデル化し、筒内実吸入燃料量が目標筒内吸入
燃料量に常に一致する様に、該プラントの伝達関数のパ
ラメータを同定・調整する適応制御器を備えてなる内燃
機関の燃料噴射量制御装置において、後に図２に示す様
に、前記プラントの伝達特性の逆伝達特性を持つ補償器
を前記プラントに直列に接続し、該補償器の伝達特性の
パラメータを前記内燃機関の運転状態に応じて予め設定
した特性に従って補正すると共に、該補償器を含んで１
つのプラントと仮想し、その仮想プラントの伝達特性が
１及びその付近から逸脱したとき、前記仮想プラント及
び適応制御器の伝達特性が全体として１及びその付近に
入る様に前記適応制御器を動作させる、如く構成した。

【０００６】

【作用】上記の如く、付着プラントの伝達特性の逆伝達
特性を持つ補償器を付着プラントに直列に接続し、その
補償器の伝達特性のパラメータ自体は運転状態に応じて
予め設定する特性に従って補正することから、例えば吸
気圧力の様に時間的変化の速いものに起因して燃料の付
着状態が変化したときも、その変化に良く追従して制御
値を変化させることができる。また、その補償器を含ん
で１つの仮想的なプラントとみなし、その仮想プラント
の持つ伝達特性、即ち、付着プラントと補償器のそれぞ
れの伝達関数の積が１及びその付近から逸脱したとき、
前記仮想プラント及び適応制御器の伝達特性が全体とし
てその範囲に入る様に適応制御器を動作させる様にした
ことから、予め設定する特性が経年変化などで実際の特
性とずれることがあっても、そのずれを良く補償するこ
とができ、変化に良く追従し筒内実吸入燃料量を目標値
に一致させることができる。

【０００７】

【発明の構成】図１はこの発明に係る制御装置を全体的
に示すブロック図であり、同図を参照して説明すると、
この制御装置は、機関回転数Ｎｅ、吸気圧力Ｐｂなどか
ら予め設定され、マップ化された特性を検索して目標筒
内吸入燃料量（Ｔｉ）を決定するＭＡＰブロック、スロ
ットル開度θTH、吸気圧力Ｐb 、大気圧Ｐa などから実
吸入空気量（Ｇａｉｒ）の動的な挙動を推定するＧａｉ
ｒモデルブロック、排気系集合部の空燃比より各気筒の
空燃比を推定するＡ／Ｆオブザーバブロック、及び燃料
噴射量（Ｔｏｕｔ）を決定する燃料噴射制御ブロックか
ら構成される。この構成において、推定された筒内実吸
入空気量Ｇａｉｒと各気筒の空燃比Ａ／Ｆとから、各時
刻（燃焼サイクル）における気筒吸入燃料量Ｇｆｕｅｌ
が推定され、それが目標燃料量Ｔｉに一致する様に、前
記燃料噴射ブロックのパラメータが調整され、噴射燃料
量Ｔｏｕｔが決定される。

【０００８】以下、それぞれについて詳細に説明する。

【０００９】先ず燃料噴射制御について説明する

【００１０】燃料噴射制御について図１を書き直すと図
２の様になる。ここで入力パラメータは、次の通りとす
る。（１）目標筒内吸入燃料量Ｔｉ．．．．．各種センサか
らの入力を用いて算出または推定した筒内実吸入空気量
を目標空燃比で割った値。尚、筒内実吸入空気量の算出
については後で詳述する。（２）筒内実吸入燃料量Ｇｆｕｅｌ．．．同様に上記の
筒内実吸入空気量を広域空燃比センサ測定値から求めた
気筒ごとの実空燃比で割った値。気筒ごとの実空燃比の
算出についても後述する。（３）その他．．．．．．．．．．．．．壁面付着補正
補償器に必要な各種測定値または推定値（例えば水温Ｔ
ｗ、吸気圧力Ｐb 、機関回転数Ｎｅなど）。

【００１１】即ち、上に述べた如く、ある時刻（ｋ−
ｎ）の燃焼サイクルの筒内実吸入空気量Ｇａｉｒ（ｋ−
ｎ）を求めて目標空燃比Ａ／Ｆ（ｋ−ｎ）で除算して目
標筒内吸入燃料量Ｔｉ（ｋ−ｎ）が決定される。また同
時刻の燃焼サイクルの筒内実吸入空気量Ｇａｉｒ（ｋ−
ｎ）を当該気筒の測定空燃比Ａ／Ｆ（ｋ−ｎ）で除算
し、筒内実吸入燃料量Ｇｆｕｅｌ（ｋ−ｎ）が決定され
る。そして適応制御器において、目標筒内吸入燃料量Ｔ
ｉ（ｋ−ｎ）に筒内実吸入燃料量Ｇｆｕｅｌ（ｋ−ｎ）
が常に一致する様に補償器の調整がなされ、制御値（噴
射燃料量）Ｔｏｕｔが決定される。ここで、前記した付
着パラメータの変化に即応するために、壁面付着プラン
トの前に、それと逆の伝達関数を持つ壁面付着補正補償
器を直列に挿入する。この壁面付着補正補償器の付着パ
ラメータは、予め機関運転状態との対応関係に基づいて
決定したマップにより検索する。もし壁面付着補正補償
器の持つ付着パラメータと実機の持つ真の付着パラメー
タとが等しければ、両者は外から見ると伝達関数が１と
なり、即ちプラントと補償器の伝達関数の積が１とな
り、目標筒内吸入燃料量＝筒内実吸入燃料量となるの
で、完全な補正が行われるはずである。ところが、付着
パラメータは一般に機関運転状態によって複雑に変化す
るため、完全に一致させることは難しい。また実機には
初期バラツキがあり、更にはデポジットの付着などによ
って経時変化も生じる。それらの理由から、両者の付着
パラメータが異なった場合、伝達関数は１（またはその
付近）以外の値となり、時間応答を生じるので、目標筒
内吸入燃料量と筒内実吸入燃料量とは等しくならない。
そこで、この付着補正補償器を含んで１つの仮想プラン
トとみなし、その仮想プラントの伝達特性が１（または
その付近の値）以外になったとき、仮想プラント及び適
応制御器の伝達特性が全体として１（またはその付近の
値）になる様に、より具体的には、仮想プラントの逆伝
達特性を持つ様に適応制御器を動作させる様にした。適
応制御器には、目標値として目標筒内吸入燃料量が入力
され、仮想プラントの出力である筒内実吸入燃料量が目
標値と一致する様に変化する適応パラメータが用いられ
る。適応制御器のパラメータは、適応パラメータ調整
（同定）器によって計算される。適応パラメータ調整
（同定）器には、仮想プラントへの過去の値を含む入出
力値が用いられる。尚、上記において適応制御器は、筒
内吸入空気量の持つ誤差を吸収する働きも行う。即ち、
筒内吸入空気量の推定に誤差があっても、結果的に筒内
吸入空気量を測定空燃比と目標空燃比でそれぞれ除算し
て得られる筒内実吸入燃料量と目標筒内吸入燃料量とが
常に一致する様に適応パラメータが調整されるので、空
気量推定誤差は吸収されることになる。

【００１２】以下、詳述する。

【００１３】壁面付着プラントとしては数１に示す様な
１次系のモデルを用いる。ここではパラメータは２個と
する。

【００１４】

【数１】

【００１５】それを離散形の伝達関数で示すと、数２の
様になる。またブロック線図で示すと、図３の様にな
る。

【００１６】

【数２】

【００１７】また数３に壁面付着補正補償器の伝達関数
を示す。前記の如く、壁面付着プラントのそれの逆伝達
関数をとる。

【００１８】

【数３】

【００１９】先に述べた如く、壁面付着補正補償器の直
接率Ａハットと持ち去り率Ｂハットは、機関運転状態、
例えば水温Ｔｗ、吸気圧力Ｐb 、機関回転数Ｎｅなどの
関数とし、予め特性をマップに設定しておき、それから
検索する。尚、この明細書で「ハット」は、推定値を意
味する。

【００２０】次いで、適応制御器について説明する。壁
面付着補正に求められる条件としては、常に輸送遅れを
小さくする方向で働くこと、Ａ，Ｂ項の変化に追従でき
ること、が挙げられるが、その様に時変（時間的に変
化）であるプラントに追従して適応制御する手法として
は、ＭＲＡＣＳ（モデル規範形適応制御）が良く知られ
ている。そこで壁面付着補正にＭＲＡＣＳを適用する
と、図４に示す様になる。このとき規範モデルは時変プ
ラントの中央値付近に取っても良いし、壁面付着補正補
償器が制御し易い様に取っても良い。尚、ＭＲＡＣＳは
無駄時間を持つプラントにのみ有効なため、付着プラン
トへの入力を１サイクル遅らせることによって見掛け上
無駄時間を挿入し、仮想プラントとしている（挿入ブロ
ックに「仮想」なる語を付す）。

【００２１】ここで気付くことは、仮想付着補正補償器
と仮想規範モデルが直列に並んでいることである。これ
らは互いに逆伝達関数の関係にあるので、キャンセルす
ることができる。その結果、仮想規範モデルの直後のｚ
^d＝ｚ（ｄ＝１）のブロックと、Ｄ（ｚ^-1）が残るが、
ｚは未来値を出力する伝達関数であり、そのままでは存
在できない。よってＤ（ｚ^-1）をＤ（ｚ^-1）＝ｚ^-1とお
くことにより、この両者もキャンセルする。通常、Ｄ
（ｚ^-1）は、Ｄ（ｚ^-1）＝１＋ｄ₁ｚ^-1＋・・・＋ｄ_n
ｚ^-n、の様におかれるが、Ｄ（ｚ^-1）＝ｚ^-1とおいても
安定性には問題がない。よって、図４を整理すると、図
５の様になる（これによって適応制御器はレギュレータ
問題を扱うこととなってＳＴＲ（セルフチューニングレ
ギュレータ）に変形する）。ここで適応制御器はパラメ
ータ同定機構によって同定された係数ベクトルを受け取
ってフィードバック補償器を形成する。但し、この動作
自体は公知なものであって、例えば「コンピュートロー
ル」Ｎｏ．２７、『ディジタル適応制御』、２８頁から
４１頁に詳述されているので、説明は省略する。

【００２２】図示した構成についてのシミュレーション
による応答結果を図６に示す。この図より、ＭＲＡＣＳ
のパラメータ同定機構は、前述の構成において正常に動
作していることが分かる。しかし空燃比の挙動には暴れ
が残ってしまっている。これを微視的にみるために、図
７の（ａ）の様な目標吸入燃料量を入力してみると、プ
ラント出力と空燃比は同図（ｂ）（ｃ）の様になる。こ
れから、プラント出力が１サイクル遅れているのが分か
る。仮想プラントとしてプラントに無駄時間が挿入され
ているため、この遅れが原因で、過渡運転時に目標燃料
量とプラントの時間差が１サイクル分生じるため、空燃
比にスパイクが生じることが分かる。

【００２３】ここで、仮想プラントを外から見た場合、
無駄時間ｚ^-1をプラントの前に挿入しても後に挿入して
も等価なので、先ずこれをプラントの後に付けることに
する。そして、プラントの出力ｙ′（ｋ）はプラントの
直後から取り出し、その後に挿入した無駄時間ｚ^-1の後
からパラメータ同定機構が必要とする仮想プラント出力
ｙ（ｋ）を取り出すことにする。こうすれば入力ｒ
（ｋ）からプラント出力ｙ′（ｋ）の経路に無駄時間が
存在せず、かつパラメータ同定機構は無駄時間の入った
仮想プラント出力ｙ（ｋ）を用いることができる。その
構成を図８に示す。また図９に図８の構成のシミュレー
ション結果を示す。収束後は図９（ｃ）の様に、実吸入
燃料量はほぼ目標燃料量となり、そのときの空燃比も１
４．７付近でフラットに推移している。また先と同じく
同定終了後の微視的応答を同じスケールでみると、図１
０の様になる（破線は比較のための無駄時間対策前の応
答）。これをみても、無駄時間対策後の応答は、同定が
終了すれば非常にフラットな空燃比になっているのが分
かる。

【００２４】尚、ここで無駄時間の挿入は上に述べたも
のに限られるものではなく、図８に想像線で示す様に、
プラントの出力の次数に対応して入力及び／又は出力に
適宜挿入することとする。

【００２５】ここで、パラメータ同定則について述べる
と、図４などに示したパラメータ同定機構においてラン
ダウなどの手法を用いるとき、ゲイン行列は数４の様に
表される。

【００２６】

【数４】

【００２７】このλ1(k), λ2(k)の選び方により、具体
的なパラメータ同定則が決定する。ＭＲＡＣＳの代表的
な同定則は、固定ゲイン法、漸減ゲイン法（最小二乗法
を含む）、可変ゲイン法（重みつき最小二乗法を含
む）、固定トレース法の４種に大別される。それらにつ
いて図４に示した構成に基づき以下の条件でシミュレー
ションを行った。即ち、実機への適用においては時変プ
ラントが対象となる可能性が高いため、時変プラントを
使用した。図１１〜図１４にそのシミュレーション結果
を示す。シミュレーション結果から分かる様に、プラン
トが時変の場合に、固定ゲイン法（図１１）ではプラン
ト出力値が目標値を中心に激しくハンチングを起こして
いる。これは、目標値が変化しているとき（過渡運転
時）に特に顕著である。過渡運転時には規範モデル出力
と、その目標値であるプラント出力値との差が大きくな
るので、ＭＲＡＣＳのパラメータ同定機構は一度に大幅
にパラメータを変えようとする。このため、プラントの
変化が速すぎるときなどの場合にはオーバシュートを起
こし、ハンチングしてしまう。それに対し、漸減ゲイン
法（図１２）と可変ゲイン法（図１３）と固定トレース
法（図１４）とでは、目標値である規範モデルに対して
プラント出力はきちんと追従している。部分的に振動が
見られるが、目標値に収束していることが分かる。この
程度の振動は、収束のためのパラメータを調整すれば、
例えばゲイン行列の数値やＤ（ｚ^-1）を変えることによ
り、収束スピードを犠牲にすることなく、抑え込むこと
が可能である。よって、この３つの同定則は収束スピー
ドが固定ゲイン法に比べて速く、プラントが時変であっ
ても追従可能であることを示している。

【００２８】続いて、筒内実吸入空気量Ｇａｉｒの推定
について説明する。

【００２９】先に述べた様に、筒内実吸入燃料量Ｇｆｕ
ｅｌを正確に求めるためには、吸入空気量を精度良く求
める必要がある。従来より、吸入空気量を直接的に計測
するマスフロー方式、吸気チャンバ内圧力によって間接
的に推定するスピードデンシティ方式などが提案されて
いるが、これらの従来手法は基本的には気筒吸入空気量
と相関の高いパラメータを用いてマッピング（マップ
化）しておき、それを検索して求めるものであるため、
マッピング時に考慮されていなかったパラメータの変化
に対しては全く無力であり、劣化、バラツキ、経年変化
などに対してタフネスがなかった。またマッピングは基
本的には定常状態でしか行うことができず、過渡運転状
態を表現していないため、過渡時の気筒吸入空気量はセ
ッティングに頼る他はなかった。従って、この発明にお
いては、吸気系に各種条件下において吸入空気量の変化
を反映することができる流体力学モデルを適用し、計測
自体は従来通り間接的であるにしても、マッピング、セ
ッティングを廃して精度良く求める様にした。即ち、ス
ロットルをオリフィスとみなし、スロットルまわりの流
体力学モデルを構築してスロットル通過空気量を推定
し、チャンバ充填遅れを考慮して動的に実吸入空気量を
推定する様にした。以下、説明する。

【００３０】先ず、図１５の吸気系モデルに示す様に、
スロットルをオリフィスとみなすと、数５に示すベルヌ
ーイの式、数６に示す連続の式、数７に示す断熱変化の
関係式より、数８に示す絞り式流量計などで使用される
圧縮性流体についての流量の計算式を導出でき、単位時
間当たりのスロットル通過空気量Ｇｔｈを求めることが
できる。

【００３１】

【数５】

【００３２】

【数６】

【００３３】

【数７】

【００３４】

【数８】

【００３５】次いで、気体の状態式に基づく数９に示す
式からチャンバ内空気量Ｇｂを求める。尚、ここで「チ
ャンバ」はいわゆるサージタンク相当部位のみならず、
スロットル弁下流から吸気ポートに至る間の全ての部位
を意味する。

【００３６】

【数９】

【００３７】従って、今回充填された空気量の変化分Δ
Ｇｂは、圧力変化分より数１０の式から求めることがで
きる。

【００３８】

【数１０】

【００３９】即ち、定常運転状態であればＧｔｈ＝Ｇａ
ｉｒとなる。他方、過渡状態で例えばスロットル弁が急
に開かれたとき吸入圧力が高くなるのは、空気がチャン
バを充填しているからである。逆に言えば、チャンバ内
を充填した空気量とスロットル弁を通過した空気量とが
分かれば気筒内に流入した空気量も分かることになる。
即ち、チャンバに充填された空気量分は当然ながら気筒
へ吸入されないものとすれば、単位時間ΔＴ当たりの筒
内実吸入空気量Ｇａｉｒは数１１で表すことができ、こ
れによって実吸入空気量の動的な挙動が推定可能とな
る。図１６にこの算出手法によるシミュレーション結果
を示す。

【００４０】

【数１１】

【００４１】上記について実験結果を示す。図１７にテ
スト装置の概略を示す。

【００４２】実験においては、スロットル開度一定と
し、空気量を変化させ、そのときのスロットル上、下流
の圧力を測定し、また同時に空気流量も測定した。テス
トはスロットル上流側に関して１０種のスロットル開度
について行った。そのうちスロットル開度３１．６度に
ついての結果を図１８に示す。図示したものも含めた実
験結果から以下が判明した。（１）スロットル下流の圧力は、スロットルから１Ｄ〜
２Ｄ（Ｄ：スロットル弁の径）で一度落ち込み、３Ｄ〜
４Ｄで回復し、再び下流へ向かうほど緩やかに圧力が下
がっていく（下流の落ち込みは、スロットル弁により流
れの縮流、渦流、剥離が発生するため）。（２）落ち込んだ圧力を測定すると、スロットルの前後
差圧が実際よりも大きくなってしまうので、回復した圧
力値を用いてスロットル通過空気量を算出する必要があ
る。

【００４３】またスロットル上流については、スロット
ル手前で圧力が低下することが判明した。

【００４４】以上から、スロットル下流の圧力Ｐthdown
（数５等の式においてＰ₂）は圧力の回復した位置、即
ち、スロットル弁から約３Ｄ（理想的には３Ｄ〜４Ｄ）
離れた位置で、スロットル上流の圧力Ｐthup（数５等の
式においてＰ₁）はスロットル弁の影響が及ばない位置
で、かつスロットル弁にできるだけ近い位置、即ち、ス
ロットル弁から約１Ｄ（ないしはそれ以上）離れた位置
で測定するのが望ましいことが分かった。尚、スロット
ル下流の圧力はその意味でチャンバ（サージタンク）内
の圧力と等価とみなすことができる。よって、後で述べ
る様にサージタンク内に圧力センサを設け、その検出値
であるサージタンク内の圧力をスロットル下流圧Ｐthdo
wnとしても良い。

【００４５】数８に示した式において、流量係数αのみ
を未知数とし、前記したテストによって流量係数を同定
した。尚、ここで同定は、測定したスロットル前後の圧
力を用い、数８の式からスロットル通過空気量Ｇｔｈを
算出し（初期値は適宜設定）、続いて算出値と実測値と
を比較し、両者が一致する様に流量係数の値を変化さ
せ、以上を繰り返して行って誤差が最小になる値を流量
係数とする、手法を採った。図１９にその手法によって
同定したスロットル開度に対する流量係数の値を示す。
また図２０に同定した流量係数を用いて推定した値と実
測値とを比較して示す（スロットル開度３１．６度につ
いてのみ示す）。

【００４６】図２１に、以上の手法で求めた流量係数を
使用し、またスロットル弁から下流に４Ｄ、上流に１Ｄ
離れた位置で測定した値を用いてシミュレーションで算
出した値と、実際に測定した値とを対比して示す。尚、
同図に示すのは、スロットル開度を７〜２０度に変化さ
せたときのデータである。またＰb は吸気圧センサ実測
値を、Ｇｔｈはエアフロメータによる実測値を示す。同
図から、この算出手法が実効的であることが見てとれよ
う。

【００４７】尚、数８において測定値はＰthup
（Ｐ₁）, Ｐthdown（Ｐ₂）, θTHであるので、それら
をマップ化（マッピング）しておいて特に平方根の演算
時間を短縮することとする。それを数１２に示す。尚、
同式でＭはマップデータを示す。

【００４８】

【数１２】

【００４９】次いで、センサの分解能との関係を述べ
る。図２２は縦軸に一定の計測誤差に対する制御誤差
を、横軸にスロットル開度をとった測定データである。
同図から、一定の計測誤差に対しては、低開度となるほ
ど、制御誤差が大きくなることが分かる。従って、セン
サは低開度側ほど計測誤差を少なくするもの、別言すれ
ば低開度側ほど高分解能を有するものであることが望ま
しい。また図２３は縦軸に同様に制御誤差をとると共
に、横軸にスロットル弁前後の圧力比をとった測定デー
タである。これから同様に、吸気圧力センサについても
高負荷側（大気圧側、図に１で示す）に高分解能を備え
るものを使用することが望ましいことが分かる。従っ
て、実機での応用に際してはスロットル開度センサと吸
気圧力センサとは共に、あるいは少なくとも一方は、こ
の様な分解能を備えたものを用いることが望ましい。

【００５０】尚、前記した空気量の測定について幾つか
付言すると、スロットル前後の圧力比が所定値以下の場
合には、流速が音速となっていることから、所定値（例
えば０．５２８）に固定する。また、吸気温センサは算
出精度を向上させる意味で、スロットル弁の上流の近傍
に設ける様にする。また湿度センサを備えて数８の式の
空気の比重量を補正するのが望ましい。

【００５１】次いで、気筒別の空燃比の検出について説
明する。多気筒内燃機関においてはコストまたは耐久性
上の問題から、一般的には排気系の集合部に１個の空燃
比センサのみを配置している。従って、集合部の空燃比
から各気筒の空燃比を特定する必要がある。そこで、集
合部の空燃比の挙動をモデル化することにより、逆に集
合部の空燃比から数値演算によって各気筒の空燃比を推
定する様にした。

【００５２】先ず、広域空燃比センサの応答遅れを１次
遅れと擬似的にモデル化し、その状態方程式を求め、周
期ΔＴで離散化すると、数１３の様になる。ここで、Ｌ
ＡＦ：広域空燃比センサ出力、Ａ／Ｆ：入力空燃比であ
る。

【００５３】

【数１３】

【００５４】数１３をＺ変換を用いて伝達関数で示せば
数１４の様になる。即ち、図２４に示す如く、数１４の
逆伝達関数を今回（時刻ｋ）のセンサ出力ＬＡＦに乗じ
ることによって前回（時刻ｋ−１）の空燃比を求めるこ
とができる。

【００５５】

【数１４】

【００５６】次いで、上記の如く遅れ補正して求めた空
燃比から各気筒の空燃比を分離抽出する手法について説
明すると、先ず図２５に示す様に内燃機関の排気系をモ
デル化する（これを図１において「ＥＸＭＮＰＬＡＮ
Ｔ」と示した）。尚、このモデル（プラント）ではＦ
（燃料量）を制御量とするため、燃空比Ｆ／Ａを用いて
いる。

【００５７】ここで、集合部の空燃比は発明者達の知見
によれば各気筒の空燃比の時間的な寄与度を考慮した加
重平均として表すことができ、式で示すと、数１５の様
になる。

【００５８】

【数１５】

【００５９】また各気筒の空燃比を漸化式の形式で表す
と、数１６の様になる。

【００６０】

【数１６】

【００６１】入力Ｕ（ｋ）は未知なので、４気筒として
４ＴＤＣごとに空燃比が再現するものとして漸化式を構
成すると、数１７の様になり、よって数１８の如き通常
の状態方程式の問題に帰着する。

【００６２】

【数１７】

【００６３】

【数１８】

【００６４】従って、時間的な寄与度Ｃが分かれば、カ
ルマンフィルタを設計して図２６に示す様なオブザーバ
を構成することにより、各時刻のＸ（ｋ）をＹ（ｋ）よ
り推定することが可能となる。即ち、上記の様な状態方
程式に対し、適当なゲイン行列をとって数１９で示す様
な式のＸハット（ｋ）を考える。

【００６５】

【数１９】

【００６６】ここで、（Ａ−ＫＣ）が安定行列であれ
ば、Ｘハット（ｋ）がＸ（ｋ）となり、Ｘ（ｋ）（各気
筒の空燃比）をＹ（ｋ）（集合部の空燃比）から推定す
ることが可能となる。尚、その詳細は本出願人が先に提
案した技術（特願平３−３５９３４０号、平成３年１２
月２７日出願）に述べてあるので、これ以上の説明は省
略する。

【００６７】以上について、実機に適用した場合につい
て実施例を説明する。

【００６８】

【実施例】図２７はそれを全体的に示す概略図である。
図において、符号１０は４気筒の内燃機関を示してお
り、吸気路１２の先端に配置されたエアクリーナ１４か
ら導入された吸気は、スロットル弁１６でその流量を調
節されつつサージタンク（チャンバ）１８とインテーク
マニホルド２０を経て第１〜第４気筒に流入される。各
気筒の吸気弁（図示せず）の付近にはインジェクタ２２
が設けられて燃料を噴射する。噴射され吸気と一体とな
った混合気は各気筒内で図示しない点火プラグで点火さ
れて燃焼してピストン（図示せず）を駆動する。燃焼後
の排気ガスは排気弁（図示せず）を介してエキゾースト
マニホルド２４に排出され、エキゾーストパイプ２６を
経て三元触媒コンバータ２８で浄化されて機関外に排出
される。

【００６９】また内燃機関１０のディストリビュータ
（図示せず）内にはピストン（図示せず）のクランク角
度位置を検出するクランク角センサ３４が設けられると
共に、スロットル弁１６の開度θTHを検出するスロット
ル開度センサ３６、スロットル弁１６下流の吸気圧力Ｐ
b を絶対圧力で検出する吸気圧センサ３８も設けられ
る。またスロットル弁１６の上流側には、大気圧Ｐa を
検出する大気圧センサ４０、吸入空気の温度を検出する
吸気温センサ４２、吸入空気の湿度を検出する湿度セン
サ４４が設けられる。更に、排気系においてエキゾース
トマニホルド２４の下流側で三元触媒コンバータ２８の
上流側には酸素濃度検出素子からなる広域空燃比センサ
４６が設けられ、排気ガスの空燃比を検出する。これら
センサ３４などの出力は、制御ユニット５０に送られ
る。尚、上記において、スロットル上流側の圧力を検出
する大気圧センサ４０は、スロットル弁１６の配置位置
から１Ｄ（Ｄ：吸気路１２の径）以上離れた位置に配置
すると共に、スロットル下流側の圧力を検出する吸気圧
センサ３８は、スロットル弁１６の配置位置から３Ｄ以
上離れたサージタンク１８内に配置する。また吸気温セ
ンサ４２と湿度センサ４４とはスロットル弁１６に可能
な限り接近して配置する。尚、スロットル開度センサ３
６と吸気圧センサ３８の分解能は、それぞれ０．０１
度、０．１ｍｍＨｇ以上とする。

【００７０】図２８は制御ユニット５０の詳細を示すブ
ロック図である。広域空燃比センサ４６の出力は検出回
路５２に入力され、そこで適宜な線型化処理が行われて
リーンからリッチにわたる広い範囲において排気ガス中
の酸素濃度に比例したリニアな特性からなる空燃比（Ａ
／Ｆ）が検出される。尚、その詳細は先に本出願人が提
案した出願（特願平３−１６９４５６号、平成３年６月
１４日出願）に述べてあるので、これ以上の説明は省略
する。検出回路５２の出力はＡ／Ｄ変換回路５４を介し
てＣＰＵ５６，ＲＯＭ５８，ＲＡＭ６０からなるマイク
ロ・コンピュータ内に取り込まれ、ＲＡＭ６０に格納さ
れる。同様にスロットル開度センサ３６などのアナログ
出力はレベル変換回路６２、マルチプレクサ６４及び第
２のＡ／Ｄ変換回路６６を介して、またクランク角セン
サ３４の出力は波形整形回路６８で波形整形された後、
カウンタ７０で出力値がカウントされ、カウント値はマ
イクロ・コンピュータ内に入力される。マイクロ・コン
ピュータにおいてＣＰＵ５６はＲＯＭ５８に格納された
命令に従って前記の如く適応制御手法に基づいて制御値
を演算し、駆動回路７２を介して各気筒のインジェクタ
２２を駆動する。

【００７１】続いて、図２８の制御装置の動作を図２９
フロー・チャートを参照して説明する。

【００７２】先ずＳ１０においてクランク角センサ３４
が検出した機関回転数Ｎｅを読み込み、Ｓ１２で大気圧
センサ４０などが検出した大気圧Ｐa （前記したスロッ
トル上流圧Ｐthupに同じ）、吸気圧力Ｐb （前記したス
ロットル下流圧Ｐthdownに同じ）、スロットル開度θT
H、空燃比Ａ／Ｆなどを読み込む。

【００７３】次いでＳ１４に進んでクランキングか否か
判断し、否定されるときはＳ１６に進んでフュエルカッ
トか否か判断する。そこでも否定されたときはＳ１８に
進み、機関回転数Ｎｅと吸気圧力Ｐb とからマップを検
索し（図１に示す）、目標筒内吸入燃料量Ｔｉを算出
し、Ｓ２０に進んで基本モードの式によって燃料噴射量
Ｔｏｕｔを算出する。尚、ここで基本モードとは前記し
た適応制御によらない手法であって、従来の公知の手法
によるものを意味する。

【００７４】次いでＳ２２に進み、広域空燃比センサ４
６の活性化が完了したか否か判断し、肯定されるときは
Ｓ２４に進んで先に述べた手法で気筒別空燃比を推定
し、Ｓ２６に進んで実吸入空気量Ｇａｉｒを推定し、Ｓ
２８に進んでそれから筒内実吸入燃料量Ｇｆｕｅｌを推
定し、Ｓ３０に進んで適応制御によって最終Ｔｏｕｔを
決定し、Ｓ３２に進んで駆動回路７２を介して当該気筒
のインジェクタ２２に出力する。尚、Ｓ１４でクランキ
ングと判断されたときはＳ３４，Ｓ３６に進んで始動モ
ードの制御値を算出すると共に、Ｓ１６でフュエルカッ
トと判断されたときはＳ３８に進んでＴｏｕｔを零とす
る。またＳ２２でセンサが活性化していないと判断され
るときは直ちにＳ３２にジャンプして基本モードによる
制御値でインジェクタを駆動する。

【００７５】上記した構成においては気筒別に空燃比を
推定して筒内実吸入燃料量を精度良く求め、筒内実吸入
燃料量が目標値に一致する様に制御器のパラメータを適
応的に制御したので、精度の良い適応制御を実現するこ
とができる。

【００７６】また、付着プラントの伝達関数の逆伝達関
数を持つ補償器を付着プラントに直列に接続し、その補
償器の伝達関数のパラメータ自体は運転状態に応じて予
め設定した特性に従って補正する様にしたことから、例
えば吸気圧力の様に時間的変化の速いものに起因して付
着状態が変化するときも、その変化に良く追従して目標
値に適応制御することができる。また、その補償器を含
んで１つの仮想プラントとみなし、その仮想プラントの
伝達関数が１（またはその付近の値）以外になったと
き、その逆伝達関数を持つ如く適応制御器を動作させる
様にしたので、予め設定する特性が経年変化などで実際
の特性と異なることがあっても良くその変化に追従して
目標値となる様に適応的に制御することができる。

【００７７】尚、上記した構成においてこの発明を図１
に関して説明したが、それに限られるものではなく、こ
の発明は図３０に示す様に、実吸入空気量の動的な挙動
を推定するＧａｉｒモデルブロックを設けず、マップ値
に”１４．７”を乗じて実吸入空気量Ｇａｉｒを推定す
る様にした構成にも妥当し、適応制御を行うことによっ
て吸気系の挙動も吸収する、即ち、前記した如く、実吸
入空気量の推定に誤差があっても、それを吸収すること
ができる。更に、この発明は図３１に示す様に、目標値
Ｔｉをマップ化せず、Ｇａｉｒモデルブロックが推定し
た実吸入空気量Ｇａｉｒに”１／１４．７”を乗じて目
標値Ｔｉを決定する様にした構成にも妥当する。

【００７８】更に、上記した構成において、１個の空燃
比センサを用いて各気筒の空燃比を推定し、目標値に制
御する例を示したが、それに限られるものではなく、気
筒ごとに空燃比センサを設けて各気筒の空燃比を直接検
出しても良い。

【００７９】更には、上記した制御を行うために、吸気
系のモデルを立てて実吸入空気量を算出したが、この手
法は開示した制御に限られるものではなく、一般的な燃
料噴射量制御や点火時期制御などにも利用することがで
きる。

【００８０】

【発明の効果】請求項１項記載の内燃機関の燃料噴射量
制御装置は、付着プラントの伝達特性の逆伝達特性を持
つ補償器を前記プラントに直列に接続し、該補償器の伝
達特性のパラメータを前記内燃機関の運転状態に応じて
予め設定した特性に従って補正すると共に、該補償器を
含んで１つのプラントと仮想し、その仮想プラントの伝
達特性が１及びその付近から逸脱したとき、前記仮想プ
ラント及び適応制御器の伝達特性が全体として１及びそ
の付近に入る様に前記適応制御器を動作させる如く構成
したので、吸気圧力の様な時間的に速い変化に起因して
燃料の付着状態が変化するときも良くそれに追従するこ
とができると共に、経年変化などの他の原因に起因する
付着状態の変化に対しても良く適応して筒内実吸入燃料
量を目標値に一致させることができる。

【００８１】請求項２項記載の装置は、前記仮想プラン
トの出力の次数に対応して入力及び／又は出力に無駄時
間を挿入する様に構成したので、空燃比の計測遅れや演
算遅れを考慮した対応づけを行った適応制御器を構成す
ることができて目標吸入燃料量とプラントの出力（実吸
入燃料量）との間に時間差が生じることがなく、よって
一層精度良く筒内実吸入燃料量を目標値に一致させるこ
とができる。

【００８２】請求項３項記載の装置は、前記パラメータ
同定・調整機構で可変ゲイン法または固定トレース法を
用いる様に構成したので、機関の過渡運転時においても
目標値にプラント出力を良く追従させることができる。

【００８３】請求項４項記載の装置は、前記内燃機関の
運転状態に応じて予め設定した特性が、少なくとも吸気
圧力または機関回転数についての特性を含む様に構成し
たので、吸気圧力の様に時間的に速く変化するものに起
因して燃料の付着状態が変化するときも、良くそれに適
応して筒内実吸入燃料量を目標値に追従制御させること
ができる。

【図面の簡単な説明】

【図１】この発明に係る制御装置を全体的に示すブロッ
ク図である。

【図２】図１中の燃料噴射制御を動作的に示すブロック
図である。

【図３】図２中の壁面付着プラントのブロック線図であ
る。

【図４】図２の壁面付着補正にＭＲＡＣＳ（モデル規範
形適応制御）を適用した状態を示すブロック図である。

【図５】図４に示すブロック図を整理した後の状態を示
すブロック図である。

【図６】図５の構成のシミュレーション結果を示すデー
タ図である。

【図７】図６のデータを微視的に検証したシミュレーシ
ョン結果を示すデータ図である。

【図８】図５の構成に無駄時間対策を施した状態を示す
ブロック図である。

【図９】図８の構成のシミュレーション結果を示すデー
タ図である。

【図１０】図９のデータを微視的に検証したシミュレー
ション結果を示すデータ図である。

【図１１】図４の構成において固定ゲイン法について行
ったシミュレーションを示すデータ図である。

【図１２】図４の構成において漸減ゲイン法について行
ったシミュレーションを示すデータ図である。

【図１３】図４の構成において可変ゲイン法について行
ったシミュレーションを示すデータ図である。

【図１４】図４の構成において固定トレース法について
行ったシミュレーションを示すデータ図である。

【図１５】図１に示したＧａｉｒモデルブロックの筒内
実吸入空気量の算出に用いる吸気系のモデルを示す説明
図である。

【図１６】図１５のモデルの筒内実吸気空気量算出のシ
ミュレーション結果を示すデータ図である。

【図１７】図１５のモデルの筒内実吸気空気量算出手法
のテスト装置を示す説明図である。

【図１８】図１７のテスト装置のテスト結果を示すデー
タ図である。

【図１９】図１７テスト装置を用いて行ったスロットル
開度に対するスロットルの流量係数の同定結果を示すデ
ータ図である。

【図２０】図１９の同定結果を用いて求めた推定値と実
測値とを対比的に示すデータ図である。

【図２１】図１５に示したモデルに基づいてシミュレー
ションで求めた値と実測値とを対比的に示すデータ図で
ある。

【図２２】制御誤差とスロットル開度の関係を示すデー
タ図である。

【図２３】制御誤差とスロットル弁前後の圧力比の関係
を示すデータ図である。

【図２４】空燃比センサの検出遅れをモデル化し、真の
空燃比を推定する推定器のブロック線図である。

【図２５】図１のＥＸＭＮＰＬＡＮＴを示すブロック
線図である。

【図２６】図２５の構成にオブザーバを組み込んだ構成
を示すブロック図である。

【図２７】図１の構成を実機に適用した状態を示す、内
燃機関の燃料噴射量制御装置を全体的に示す概略図であ
る。

【図２８】図２７の制御ユニットの構成を示すブロック
図である。

【図２９】図２８の制御ユニットの動作を示すフロー・
チャートである。

【図３０】この発明の別の構成例を示す図１と同様のブ
ロック図である。

【図３１】この発明の更に別の構成例を示す図１と同様
のブロック図である。

【符号の説明】

１０内燃機関１２吸気路１６スロットル弁１８サージタンク２０インテークマニホルド２２インジェクタ２４エキゾーストマニホルド２６エキゾーストパイプ３６スロットル開度センサ３８吸気圧センサ４０大気圧センサ４２吸気温センサ４４湿度センサ４６広域空燃比センサ５０制御ユニット

───────────────────────────────────────────────────── フロントページの続き (72)発明者小森谷勲埼玉県和光市中央１丁目４番１号株式会社本田技術研究所内 (72)発明者廣田俊明埼玉県和光市中央１丁目４番１号株式会社本田技術研究所内

Claims

【特許請求の範囲】

【請求項１】内燃機関の吸気管に付着する燃料量を状
態変数として燃料の挙動をモデル化し、筒内実吸入燃料
量が目標筒内吸入燃料量に常に一致する様に、そのプラ
ントの伝達関数のパラメータを同定・調整する適応制御
器を備えてなる内燃機関の燃料噴射量制御装置におい
て、前記プラントの伝達特性の逆伝達特性を持つ補償器
を前記プラントに直列に接続し、該補償器の伝達特性の
パラメータを前記内燃機関の運転状態に応じて予め設定
した特性に従って補正すると共に、該補償器を含んで１
つのプラントと仮想し、その仮想プラントの伝達特性が
１及びその付近から逸脱したとき、前記仮想プラント及
び適応制御器の伝達特性が１及びその付近に入る様に前
記適応制御器を動作させる様に構成したことを特徴とす
る内燃機関の燃料噴射量制御装置。
【請求項２】前記仮想プラントの出力の次数に対応し
て入力及び／又は出力に無駄時間を挿入することを特徴
とする請求項１項記載の内燃機関の燃料噴射量制御装
置。
【請求項３】前記パラメータ同定・調整機構で可変ゲ
イン法または固定トレース法を用いることを特徴とする
請求項１項または２項記載の内燃機関の燃料噴射量制御
装置。
【請求項４】前記内燃機関の運転状態に応じて予め設
定した特性が、少なくとも吸気圧力または機関回転数に
ついての特性を含むことを特徴とする請求項１項ないし
３項のいずれかに記載の内燃機関の燃料噴射量制御装
置。