JPH0327758B2 - - Google Patents
Info
- Publication number
- JPH0327758B2 JPH0327758B2 JP57125885A JP12588582A JPH0327758B2 JP H0327758 B2 JPH0327758 B2 JP H0327758B2 JP 57125885 A JP57125885 A JP 57125885A JP 12588582 A JP12588582 A JP 12588582A JP H0327758 B2 JPH0327758 B2 JP H0327758B2
- Authority
- JP
- Japan
- Prior art keywords
- pnax
- fuel ratio
- ignition timing
- air
- 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 - Lifetime
Links
- 239000000446 fuel Substances 0.000 claims description 35
- 238000002485 combustion reaction Methods 0.000 claims description 18
- 238000000034 method Methods 0.000 claims description 17
- 238000012545 processing Methods 0.000 claims description 5
- 238000002347 injection Methods 0.000 claims description 3
- 239000007924 injection Substances 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 6
- 239000011159 matrix material Substances 0.000 description 5
- 238000001514 detection method Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000001052 transient effect Effects 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000012546 transfer Methods 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 230000006399 behavior Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02P—IGNITION, OTHER THAN COMPRESSION IGNITION, FOR INTERNAL-COMBUSTION ENGINES; TESTING OF IGNITION TIMING IN COMPRESSION-IGNITION ENGINES
- F02P5/00—Advancing or retarding ignition; Control therefor
- F02P5/04—Advancing or retarding ignition; Control therefor automatically, as a function of the working conditions of the engine or vehicle or of the atmospheric conditions
- F02P5/145—Advancing or retarding ignition; Control therefor automatically, as a function of the working conditions of the engine or vehicle or of the atmospheric conditions using electrical means
- F02P5/1455—Advancing or retarding ignition; Control therefor automatically, as a function of the working conditions of the engine or vehicle or of the atmospheric conditions using electrical means by using a second control of the closed loop type
-
- 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/40—Engine management systems
Landscapes
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Electrical Control Of Ignition Timing (AREA)
- Combined Controls Of Internal Combustion Engines (AREA)
Description
【発明の詳細な説明】
(技術分野)
この発明は、内燃機関において、気筒内圧力情
報から点火時期を適宜調整して機関の出力又は効
率が最大となる(これを通常MBT制御:
Minimun advance for Best Torqueと称する)
ように制御し、かつ所定空燃比(A/F)に同時
に制御する燃焼制御方法に関する。
報から点火時期を適宜調整して機関の出力又は効
率が最大となる(これを通常MBT制御:
Minimun advance for Best Torqueと称する)
ように制御し、かつ所定空燃比(A/F)に同時
に制御する燃焼制御方法に関する。
(背景技術)
機関の燃焼制御方法としては、燃焼によるシリ
ンダ内圧を検出し、この圧力が最大となるクラン
ク角位置が、予め定められた所定値となるように
点火時期を制御する方法がある。機関によつて多
少の差異はあるが、燃焼室内圧が最大となる位置
(以後θPnaxとする)が、クランク角で上死点後
(以後ATDCとする)10゜〜20゜の位置にくるよう
に点火時期を設定した場合に、その機関の発生ト
ルクが最大となる。従つて、点火時期やフイード
バツク制御して、燃焼室内圧が最大となるクラン
ク角位置がATDC10゜〜20゜となるように調整する
ことで、発生トルクを最大にすることができる。
ンダ内圧を検出し、この圧力が最大となるクラン
ク角位置が、予め定められた所定値となるように
点火時期を制御する方法がある。機関によつて多
少の差異はあるが、燃焼室内圧が最大となる位置
(以後θPnaxとする)が、クランク角で上死点後
(以後ATDCとする)10゜〜20゜の位置にくるよう
に点火時期を設定した場合に、その機関の発生ト
ルクが最大となる。従つて、点火時期やフイード
バツク制御して、燃焼室内圧が最大となるクラン
ク角位置がATDC10゜〜20゜となるように調整する
ことで、発生トルクを最大にすることができる。
第1図は、この方法を用いた点火時期制御装置
の原理を示すブロツク図である。同図において、
1はシリンダ内圧検出器、2はクランク角位置検
出器、3はθPnax検出回路、4はメモリ、5はP
nax演算回路、6は比較器、7は点火時期調整器
である。シリンダ内圧検出器1で検出されるシリ
ンダ内圧Pとクランク角位置検出器2で検出され
るクランク角回転位置θとの関係を第2図に示
す。同図において、Pnaxは1サイクルについてシ
リンダ内圧Pが最大となる値を示し、θPnaxはこ
の時の回転位置を示す。
の原理を示すブロツク図である。同図において、
1はシリンダ内圧検出器、2はクランク角位置検
出器、3はθPnax検出回路、4はメモリ、5はP
nax演算回路、6は比較器、7は点火時期調整器
である。シリンダ内圧検出器1で検出されるシリ
ンダ内圧Pとクランク角位置検出器2で検出され
るクランク角回転位置θとの関係を第2図に示
す。同図において、Pnaxは1サイクルについてシ
リンダ内圧Pが最大となる値を示し、θPnaxはこ
の時の回転位置を示す。
θPnax検出回路3は、シリンダ内圧検出器1及び
クランク角位置検出器2からの信号を入力し、連
続するサイクルについて各サイクルごとに、
(θPnax)1、(θPnax)2、…(θPnax)k、…(θPna
x)oを
求めメモリ4は、各サイクルごとの(θPnax)k、
k=1〜nを記憶する。Pnax演算回路5は、所
定回数n個の(θPnax)kの平均値Pnax Pnax=1/no 〓k=1 (θPnax)k ……(1) を演算する。
クランク角位置検出器2からの信号を入力し、連
続するサイクルについて各サイクルごとに、
(θPnax)1、(θPnax)2、…(θPnax)k、…(θPna
x)oを
求めメモリ4は、各サイクルごとの(θPnax)k、
k=1〜nを記憶する。Pnax演算回路5は、所
定回数n個の(θPnax)kの平均値Pnax Pnax=1/no 〓k=1 (θPnax)k ……(1) を演算する。
前述したように、機関の出力又は効率が最大に
なるように制御するためには、このPnaxが基準
値Cnio〜Cnax度ATDCの範囲内にあるように、点
火時期を調整する必要がある。従つて、比較器6
はこのPnaxをCnio〜Cnaxと比較し、 Cnio PnaxCnax ……(2) であれば、機関はMBT制御状態であると判定す
る。また Pnax<Cnio ……(3) であれば、機関はMBT制御状態ではないと判定
し、点火時期を ΔADV=f1(Pnax−Cnio) ……(4) (但し、f1は関数であることを示す)だけ遅らせ
る。同様に Pnax>Cnax であれば、機関はMBT制御状態ではないと判定
し、点火時期を ΔADV=f2(Pnax−Cnax) ……(5) だけ進ませる。これらの点火時期の調整は、点火
時期調整器7で行なう。
なるように制御するためには、このPnaxが基準
値Cnio〜Cnax度ATDCの範囲内にあるように、点
火時期を調整する必要がある。従つて、比較器6
はこのPnaxをCnio〜Cnaxと比較し、 Cnio PnaxCnax ……(2) であれば、機関はMBT制御状態であると判定す
る。また Pnax<Cnio ……(3) であれば、機関はMBT制御状態ではないと判定
し、点火時期を ΔADV=f1(Pnax−Cnio) ……(4) (但し、f1は関数であることを示す)だけ遅らせ
る。同様に Pnax>Cnax であれば、機関はMBT制御状態ではないと判定
し、点火時期を ΔADV=f2(Pnax−Cnax) ……(5) だけ進ませる。これらの点火時期の調整は、点火
時期調整器7で行なう。
第3図に、これらの動作のフローチヤートを示
す。
す。
一方、空燃比(燃料と空気の混合気)の制御
は、先ず機関の回転速度Nと吸入空気流量Qから
基本燃料供給量TPをTP=KQ/N(但しKは定数)
によつて求める。そして、排気混合気の酸素濃度
に応じて空燃比に応じた信号を出力するO2セン
サの出力値に基づいて、基本燃料供給量TPに対
する補正率αをPI制御(比例積分制御)するこ
とにより、実際の空燃比A/Fを所定空燃比にフ
イードバツク制御している。
は、先ず機関の回転速度Nと吸入空気流量Qから
基本燃料供給量TPをTP=KQ/N(但しKは定数)
によつて求める。そして、排気混合気の酸素濃度
に応じて空燃比に応じた信号を出力するO2セン
サの出力値に基づいて、基本燃料供給量TPに対
する補正率αをPI制御(比例積分制御)するこ
とにより、実際の空燃比A/Fを所定空燃比にフ
イードバツク制御している。
しかしながら、このような燃焼制御方法にあつ
ては、確率的にばらつくθPnaxを平均操作した量
を見て行なうが、平均操作での時間的な遅れ、あ
るいは点火時期とθPnaxの動的挙動(dynamic
behaivor)を考慮せずにフイードバツク制御す
る構成となつている上、空燃比フイードバツク制
御とは全く独立に行なう構成となつているため、
過渡時にはθPnaxを所定値にレギユレートするこ
とが困難で、しかも所定空燃比からのズレが激し
く、燃費及び運転性の悪化を招くという問題点が
あつた。
ては、確率的にばらつくθPnaxを平均操作した量
を見て行なうが、平均操作での時間的な遅れ、あ
るいは点火時期とθPnaxの動的挙動(dynamic
behaivor)を考慮せずにフイードバツク制御す
る構成となつている上、空燃比フイードバツク制
御とは全く独立に行なう構成となつているため、
過渡時にはθPnaxを所定値にレギユレートするこ
とが困難で、しかも所定空燃比からのズレが激し
く、燃費及び運転性の悪化を招くという問題点が
あつた。
(発明の目的)
この発明は、このような従来の問題点に着目し
てなされたもので、過渡時に特にその振舞いが生
じる燃焼のダイナミクス及び燃料と空燃比間のダ
イナミクスを考慮して点火時期燃料噴射制御をす
ることにより、上記問題点を解決することを目的
とする。
てなされたもので、過渡時に特にその振舞いが生
じる燃焼のダイナミクス及び燃料と空燃比間のダ
イナミクスを考慮して点火時期燃料噴射制御をす
ることにより、上記問題点を解決することを目的
とする。
(発明の構成及び作用)
そこでこの発生は、従来のPID制御とは異な
り、以下に詳細に説明するδITとδFuelとを制御
入力とし、δθ^PnaxとδA/Fとを制御出力とする制
御対象である機関のダイナミツクモデルに基づい
て、上記各入出力とで多変数制御することを特徴
とするものである。
り、以下に詳細に説明するδITとδFuelとを制御
入力とし、δθ^PnaxとδA/Fとを制御出力とする制
御対象である機関のダイナミツクモデルに基づい
て、上記各入出力とで多変数制御することを特徴
とするものである。
以下、この発明を図面に基づいて説明する。
第4図はこの発明の一実施例を示す図である。
まず構成を説明する。8は、ある運動条件で
MBTとなる点火時期からの摂動分を表わすδIT
及び所定空燃比となる燃料噴射量からの摂動分
δFuelを入力とし、MBTとなる気筒内圧力が最
大となるクランク角度からの摂動分を統計処理
(例えば移動平均)した量を表わすδθ^Pnax及び所
定空燃比からのズレ分δA/Fを出力とする制御
対象である機関の線形化されたシステムの回転同
期サンプル値系の、実験から得られた伝達関数行
列T(z)、 T(z)=T11(z)、T12(z) T21(z)、T22(z) ……(6) を示す。
まず構成を説明する。8は、ある運動条件で
MBTとなる点火時期からの摂動分を表わすδIT
及び所定空燃比となる燃料噴射量からの摂動分
δFuelを入力とし、MBTとなる気筒内圧力が最
大となるクランク角度からの摂動分を統計処理
(例えば移動平均)した量を表わすδθ^Pnax及び所
定空燃比からのズレ分δA/Fを出力とする制御
対象である機関の線形化されたシステムの回転同
期サンプル値系の、実験から得られた伝達関数行
列T(z)、 T(z)=T11(z)、T12(z) T21(z)、T22(z) ……(6) を示す。
9は制御対象である機関のダイナミツク状態情
報を記憶していて、入力δITとδFuel及び出力
δθ^PnaxとδA/Fから機関のダイナミツクな内部状
態を推定する状態観測器であり、内部状態を代表
する状態変数量δx(例えばδx1、δx2、δx3、δx4を
要素とする状態変数ベクトル)の推定値δを計
算する。
報を記憶していて、入力δITとδFuel及び出力
δθ^PnaxとδA/Fから機関のダイナミツクな内部状
態を推定する状態観測器であり、内部状態を代表
する状態変数量δx(例えばδx1、δx2、δx3、δx4を
要素とする状態変数ベクトル)の推定値δを計
算する。
状態観測器9は、制御対象である機関をシミユ
レーシヨンするもので、ダイナミツクな内部状態
を状態変数δx(n次のベクトルδx1〜δxo)で代表
させる。この状態変数δxは、実際の内部状態を
表わす種々の物理量に対応させる必要はない。ま
た状態変数δxの次数nは、nが大きい程シミユ
レーシヨンが精確になるが、反面計算が複雑にな
る。この発明の2入力2出力の場合には、n=4
程度が妥当である。
レーシヨンするもので、ダイナミツクな内部状態
を状態変数δx(n次のベクトルδx1〜δxo)で代表
させる。この状態変数δxは、実際の内部状態を
表わす種々の物理量に対応させる必要はない。ま
た状態変数δxの次数nは、nが大きい程シミユ
レーシヨンが精確になるが、反面計算が複雑にな
る。この発明の2入力2出力の場合には、n=4
程度が妥当である。
10は、状態観測器9により推定して得られた
状態情報をゲイン−K倍して入力フイードバツク
するゲインブロツクである。
状態情報をゲイン−K倍して入力フイードバツク
するゲインブロツクである。
次に作用を説明する。
一般に機関は、同一運転条件で同じ点火時期、
燃料噴射量であつても、燃焼挙動は確率的にばら
つき、気筒内圧の最大となるクランク角度θPnax
も確率的に分散するこのようなシステムに対し確
定的に制御を行なうには、何らかの統計処理をし
て確定量を抽出しなければならない。
燃料噴射量であつても、燃焼挙動は確率的にばら
つき、気筒内圧の最大となるクランク角度θPnax
も確率的に分散するこのようなシステムに対し確
定的に制御を行なうには、何らかの統計処理をし
て確定量を抽出しなければならない。
この発明は、ある機関に対してはMBTとなる
θPnaxが一定になるという事実に基づき、点火時
期制御により常にθPnaxがMBTとなるようレギユ
レータ制御し、かつ所定(A/F)*にレギユレー
タ制御する方法であり、確率的にばらつくθPnax
に統計処理し、又O2センサのリツチリーン信号
からδA/Fを算出し制御出力として用いる。統
計処理の例としては平均があるが、単純加算平均
では長時間の観察が必要となり、過渡時には不利
であり、ここでは移動平均δθ^Pnaxをとる。移動平
均することや、燃料供給から空燃比を検出するま
での間でも当然遅れ時間を生ずることになるが、
この遅れ時間は機関の伝達関数行列T(z)の中
に含まれる。
θPnaxが一定になるという事実に基づき、点火時
期制御により常にθPnaxがMBTとなるようレギユ
レータ制御し、かつ所定(A/F)*にレギユレー
タ制御する方法であり、確率的にばらつくθPnax
に統計処理し、又O2センサのリツチリーン信号
からδA/Fを算出し制御出力として用いる。統
計処理の例としては平均があるが、単純加算平均
では長時間の観察が必要となり、過渡時には不利
であり、ここでは移動平均δθ^Pnaxをとる。移動平
均することや、燃料供給から空燃比を検出するま
での間でも当然遅れ時間を生ずることになるが、
この遅れ時間は機関の伝達関数行列T(z)の中
に含まれる。
まず、T(z)から機関の動的な振舞いを記述
する状態変数モデル δx(n)=Aδx(n-1)+Bδu(n-1) ……(7) δy(n-1)=Cδx(n-1) ……(8) を導く。各量カツコ内の(n)は現時点を、また
(n−1)は1つ前の時点を表わす。δu(n−1)
は制御入力で、δIT(n−1)、δFuel(n−1)を
要素とする。δy(n−1)は制御出力で、δθPnax
(n−1)、δA/Fを要素とする。行列A、B、
Cは運転条件ごとの伝達関数行列T(z)の係数
から決まる定数行列である。
する状態変数モデル δx(n)=Aδx(n-1)+Bδu(n-1) ……(7) δy(n-1)=Cδx(n-1) ……(8) を導く。各量カツコ内の(n)は現時点を、また
(n−1)は1つ前の時点を表わす。δu(n−1)
は制御入力で、δIT(n−1)、δFuel(n−1)を
要素とする。δy(n−1)は制御出力で、δθPnax
(n−1)、δA/Fを要素とする。行列A、B、
Cは運転条件ごとの伝達関数行列T(z)の係数
から決まる定数行列である。
ここで、次のようなアルゴリズムを持つ状態観
測器を構成する。
測器を構成する。
δ(n)=(A−GC)δ(n−1)
+Bδu(n−1)+Gδy(n−1) ……(9)
ここに、Gは任意に与えられる行列で、δ
(・)((・)は任意のサンプル点を表わす)は推
定された機関の内部状態変数である。(7)(8)(9)式よ
り変形すると 〔 δx(n)−δ(n)〕 =(A−GC)〔δx(n−1) −δ(n−1)〕 ……(10) となり、(A−GC)行列の固有値が単位円内にあ
るようにGを選べば n→大で δ(n)→δx(n) となり、内部状態δx(n)を入出力より推定する
ことができる。
(・)((・)は任意のサンプル点を表わす)は推
定された機関の内部状態変数である。(7)(8)(9)式よ
り変形すると 〔 δx(n)−δ(n)〕 =(A−GC)〔δx(n−1) −δ(n−1)〕 ……(10) となり、(A−GC)行列の固有値が単位円内にあ
るようにGを選べば n→大で δ(n)→δx(n) となり、内部状態δx(n)を入出力より推定する
ことができる。
また、Gを適当に選び、(A−GC)の固有値が
全て零とすることも可能で、この時状態観測器は
有限整定状態観測器となる。
全て零とすることも可能で、この時状態観測器は
有限整定状態観測器となる。
このように推定された状態情報から、θPnax及
び空燃比の一定値制御を行なうため、次のような
状態フイードバツクを施す。
び空燃比の一定値制御を行なうため、次のような
状態フイードバツクを施す。
δu=−Kδ(n) ……(11)
ここで、ゲインKを、次の評価関数Jが最小に
なるように決定する。
なるように決定する。
J=∞
〓k=0
〔δytQδy+δutRu〕 ……(12)
ただし、Q、Rは重みパラメータである。R
は、制御入力の制約項の重みを調整するもので、
δIT、δFuelの制御範囲、過度特性を見て決定す
る。(12)式の評価関数は、所定θ* Pnaxと所定(A/
F)*との偏差δθ^2 nax、δA/F2を常に最小にするも
ので、このようなゲインKは次のリカツテイ
(Riccati)方程式 P=AtPA−AtPB(BtPB+R) -1BtPA+Q ……(13) を解いたPを用いて、 K=(R+BtPB)-1BtPA ……(14) となる。従つて、運転条件ごとのA、B、Cを用
いて解いて求めたゲインKを記憶しておき、その
運転条件の合つたゲインをルツクアツプすればよ
い。
は、制御入力の制約項の重みを調整するもので、
δIT、δFuelの制御範囲、過度特性を見て決定す
る。(12)式の評価関数は、所定θ* Pnaxと所定(A/
F)*との偏差δθ^2 nax、δA/F2を常に最小にするも
ので、このようなゲインKは次のリカツテイ
(Riccati)方程式 P=AtPA−AtPB(BtPB+R) -1BtPA+Q ……(13) を解いたPを用いて、 K=(R+BtPB)-1BtPA ……(14) となる。従つて、運転条件ごとのA、B、Cを用
いて解いて求めたゲインKを記憶しておき、その
運転条件の合つたゲインをルツクアツプすればよ
い。
以上の燃焼制御手順を示したのが第5図であ
る。手順を説明すると、まず気筒内圧検出器及び
クランク角位置からの情報により、筒内圧最大と
なるクランク角度θPnaxを検出し、またO2センサ
からのリツチ・リーンのオン・オフ信号を検出す
る(ステツプ11)。次にθPnaxを移動平均し、
MBTとなる。θ* Pnaxと比較し、δθ^Pnaxを計算する
(ステツプ12)。また、第6図に示すように、O2
センサからのリツチ・リーンのオン・オフ信号
(O2センサは所定空燃比を境にして燃料のリツチ
(濃)側でオン信号を、リーン(淡)側でオフ信
号をそれぞれ出力する)を各制御周期ごと観測す
る。例えば最初の周期(0〜1)でオンの時間と
オフの時間を計測し、オン(リツチ)信号を
(+)、オフ(リーン)信号を(−)として加算
し、 δA/F=−t1+t2−t3 ……(15) により得られたδA/Fの値をもつて、その制御
周期内で空燃比が所定値(A/F)*よりどれだけ
ズレているかを表わす量とすればよい(ステツプ
12)。次に(9)式のアルゴリズムに従い状態推定を
行なう(ステツプ13)。この例では、 A−GC=0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 0 ……(16) となるGを選択し、有限整定状態観測器を構成し
ている。次にステツプ16で推定されたダイナミツ
ク状態情報δ1、δ2、δ3、δ4を用いて
フイードバツクゲイン−Kを乗じ、最適制御入力
の基準値からの増量分を計算する。以上の手順中
の係数B、G、Kは、運転条件ごとに求めてお
き、記憶し、適宜ルツクアツプする。
る。手順を説明すると、まず気筒内圧検出器及び
クランク角位置からの情報により、筒内圧最大と
なるクランク角度θPnaxを検出し、またO2センサ
からのリツチ・リーンのオン・オフ信号を検出す
る(ステツプ11)。次にθPnaxを移動平均し、
MBTとなる。θ* Pnaxと比較し、δθ^Pnaxを計算する
(ステツプ12)。また、第6図に示すように、O2
センサからのリツチ・リーンのオン・オフ信号
(O2センサは所定空燃比を境にして燃料のリツチ
(濃)側でオン信号を、リーン(淡)側でオフ信
号をそれぞれ出力する)を各制御周期ごと観測す
る。例えば最初の周期(0〜1)でオンの時間と
オフの時間を計測し、オン(リツチ)信号を
(+)、オフ(リーン)信号を(−)として加算
し、 δA/F=−t1+t2−t3 ……(15) により得られたδA/Fの値をもつて、その制御
周期内で空燃比が所定値(A/F)*よりどれだけ
ズレているかを表わす量とすればよい(ステツプ
12)。次に(9)式のアルゴリズムに従い状態推定を
行なう(ステツプ13)。この例では、 A−GC=0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 0 ……(16) となるGを選択し、有限整定状態観測器を構成し
ている。次にステツプ16で推定されたダイナミツ
ク状態情報δ1、δ2、δ3、δ4を用いて
フイードバツクゲイン−Kを乗じ、最適制御入力
の基準値からの増量分を計算する。以上の手順中
の係数B、G、Kは、運転条件ごとに求めてお
き、記憶し、適宜ルツクアツプする。
(発明の効果)
以上説明してきたように、この発明によれば、
その構成を点火時期とθPnax、供給燃料と空燃比、
更には確率的にばらつくθPnaxの統計処理で生じ
る遅れ等のダイナミクスを考慮した上で最適制御
を行なう方法としたため、過渡時もθPnaxを所定
値にレギユレートでき、しかも所定空燃比にレギ
ユレートでき、運転性、燃費を考慮した最大トル
ク運転ができるという効果が得られる。第7図の
加速時に従来法(A)と本方法(B)を適用した時の効果
を示した。本方法の制御性の良さが示されてい
る。
その構成を点火時期とθPnax、供給燃料と空燃比、
更には確率的にばらつくθPnaxの統計処理で生じ
る遅れ等のダイナミクスを考慮した上で最適制御
を行なう方法としたため、過渡時もθPnaxを所定
値にレギユレートでき、しかも所定空燃比にレギ
ユレートでき、運転性、燃費を考慮した最大トル
ク運転ができるという効果が得られる。第7図の
加速時に従来法(A)と本方法(B)を適用した時の効果
を示した。本方法の制御性の良さが示されてい
る。
第1図は従来の内燃機関における燃焼制御装置
のブロツク図、第2図は機関軸回転位置とシリン
ダ内圧との関係を示すグラフ、第3図は第1図の
装置の動作を示すフローチヤート、第4図はこの
発明による内燃機関における燃焼制御方法を実現
する制御装置の構成図、第5図はこの発明による
制御方法を説明するフローチヤート、第6図は
O2センサの出力波形図である。第7図は本方法
の効果を示した図である。 1……シリンダ内圧検出器、2……クランク角
位置検出器、3……θPnax検出回路、4……メモ
リ、5……Pnax演算回路、6……比較器、7…
…点火時期調整器、8……内燃機関制御対象、9
……状態観測器、10……ゲインブロツク、
θPnax……気筒内圧力が最大となるクランク角位
置、Cnax,Cnio……θPnaxの所定値。
のブロツク図、第2図は機関軸回転位置とシリン
ダ内圧との関係を示すグラフ、第3図は第1図の
装置の動作を示すフローチヤート、第4図はこの
発明による内燃機関における燃焼制御方法を実現
する制御装置の構成図、第5図はこの発明による
制御方法を説明するフローチヤート、第6図は
O2センサの出力波形図である。第7図は本方法
の効果を示した図である。 1……シリンダ内圧検出器、2……クランク角
位置検出器、3……θPnax検出回路、4……メモ
リ、5……Pnax演算回路、6……比較器、7…
…点火時期調整器、8……内燃機関制御対象、9
……状態観測器、10……ゲインブロツク、
θPnax……気筒内圧力が最大となるクランク角位
置、Cnax,Cnio……θPnaxの所定値。
Claims (1)
- 1 内燃機関の気筒内圧力情報に基づき、点火時
期及び燃料供給量により機関の燃焼を制御する方
法において、点火時期及び燃料供給量と、気筒内
圧力が最大となるクランク角度を統計処理した量
θ^Pnaxと空燃比A/Fとの間の確定的ダイナミツ
クモデルに基づき、ある運転条件でMBT
(Minimun advance for Best Torque)となる
点火時期との摂動分δIT及び所定空燃比(A/
F)*となる燃料噴射量からの摂動分δFuelと、あ
る運転条件でMBTとなる気筒内圧力が最大とな
るクランク角度θ* Pnaxに対する気筒内圧検出器及
びクランク角位置検出器により検出された気筒内
圧力が最大となる各サイクルについて求めた所定
回数のクランク角度θPnaxを統計処理した量θ^Pnax
の摂動分δθ^Pnax及び所定空燃比(A/F)*からの
摂動分δA/Fとから、前記内燃機関のダイナミ
ツクな内部状態を代表する適当な次数の状態変数
量δxi(i=1、2、…n)を前記確定的ダイナミ
ツクモデルが推定し、該推定された状態変数量δ
xiをゲインki(i=1、2、…n)倍して点火
時期及び燃料供給量にフイードバツクすることに
より前記θ^Pnax及び空燃比A/Fが所定量となる
ように制御することを特徴とする燃焼制御方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP57125885A JPS5918269A (ja) | 1982-07-21 | 1982-07-21 | 燃焼制御方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP57125885A JPS5918269A (ja) | 1982-07-21 | 1982-07-21 | 燃焼制御方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPS5918269A JPS5918269A (ja) | 1984-01-30 |
JPH0327758B2 true JPH0327758B2 (ja) | 1991-04-16 |
Family
ID=14921326
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP57125885A Granted JPS5918269A (ja) | 1982-07-21 | 1982-07-21 | 燃焼制御方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JPS5918269A (ja) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS62195463A (ja) * | 1986-02-19 | 1987-08-28 | Honda Motor Co Ltd | 内燃機関の点火時期制御装置 |
US5070842A (en) * | 1989-07-19 | 1991-12-10 | Mitsubishi Denki Kabushiki Kaisha | Apparatus for controlling ignition timing in internal combustion engine |
-
1982
- 1982-07-21 JP JP57125885A patent/JPS5918269A/ja active Granted
Also Published As
Publication number | Publication date |
---|---|
JPS5918269A (ja) | 1984-01-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US5497329A (en) | Prediction method for engine mass air flow per cylinder | |
US6382198B1 (en) | Individual cylinder air/fuel ratio control based on a single exhaust gas sensor | |
US5531208A (en) | Air-fuel ratio feedback control system for internal combustion engine | |
US5080064A (en) | Adaptive learning control for engine intake air flow | |
US4969439A (en) | Adaptive control system for an internal combustion engine | |
GB2329979A (en) | Controlling running smoothness of an internal combustion engine | |
JPH01182552A (ja) | 空燃比適応制御装置 | |
JPH0240044A (ja) | 内燃機関のスロットル開度制御装置 | |
US4469074A (en) | Electronic control for internal combustion engine | |
US5690087A (en) | EGO based adaptive transient fuel compensation for a spark ignited engine | |
US5569847A (en) | Air-fuel ratio estimator for internal combustion engine | |
JPH11159377A (ja) | エンジン制御装置 | |
US5638788A (en) | Automotive actuator interface | |
JPH09287513A (ja) | エンジンのトルク制御装置 | |
US5762054A (en) | Ego based adaptive transient fuel compensation for a spark ignited engine | |
EP0210766B1 (en) | Adaptive control system for an internal combustion engine | |
JPH0327758B2 (ja) | ||
JP2689362B2 (ja) | 内燃機関の空燃比検出方法 | |
JPH0320591B2 (ja) | ||
JPS61232355A (ja) | 内燃機関のアイドリング速度を制御する方法 | |
US5101787A (en) | Ignition timing control apparatus for an engine | |
JP6543509B2 (ja) | 内燃機関の制御装置 | |
JPH0320590B2 (ja) | ||
JPH0368224B2 (ja) | ||
JP3055378B2 (ja) | 内燃機関の制御装置 |