JP2004251123A - Exhaust emission control device for internal combustion engine - Google Patents
Exhaust emission control device for internal combustion engine Download PDFInfo
- Publication number
- JP2004251123A JP2004251123A JP2003039170A JP2003039170A JP2004251123A JP 2004251123 A JP2004251123 A JP 2004251123A JP 2003039170 A JP2003039170 A JP 2003039170A JP 2003039170 A JP2003039170 A JP 2003039170A JP 2004251123 A JP2004251123 A JP 2004251123A
- Authority
- JP
- Japan
- Prior art keywords
- value
- amount
- air
- fuel ratio
- catalyst
- 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.)
- Pending
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/40—Engine management systems
Abstract
Description
【0001】
【発明の属する技術分野】
本発明は、内燃機関の排気通路に配設された触媒の上流側に空燃比センサを備え、同空燃比センサの出力値に基づいて同触媒に流入するガスの空燃比をフィードバック制御する内燃機関の排気浄化装置に関する。
【0002】
【従来の技術】
従来より、内燃機関の排ガスを浄化するための三元触媒(本明細書においては、単に「触媒」と云うこともある。)が、同機関の排気通路に配設されている。この三元触媒は、酸素を吸蔵するO2ストレージ機能(以下、この機能を「酸素吸蔵機能」と称呼し、触媒内に吸蔵されている酸素量を「酸素吸蔵量」と称呼する。)を有していて、流入する排ガスの空燃比がリッチである場合には吸蔵している酸素によりHC,CO等の未燃成分を酸化するとともに、流入するガスの空燃比がリーンである場合には窒素酸化物(NOx)を還元して同NOxから奪った酸素を内部に吸蔵する。換言すれば、三元触媒に流入する排ガスの空燃比がリッチである場合には同三元触媒の酸素吸蔵量が減少するとともに三元触媒に流入する排ガスの空燃比がリーンである場合には同三元触媒の酸素吸蔵量が増加する。これにより、三元触媒は、機関の空燃比が理論空燃比から偏移した場合でも、未燃成分や窒素酸化物等の有害成分を浄化することができる。
【0003】
ところで、内燃機関の運転状態は時々刻々変化するから、機関の空燃比はリッチとなったり、リーンとなったりする。一方、触媒の酸素吸蔵量は、「0」から吸蔵し得る最大の酸素量(以下、「最大酸素吸蔵量」と称呼する。)の間で変化する。従って、仮に、触媒内の酸素吸蔵量が「0」近傍であるときに機関の空燃比がリッチとなると、触媒内で未燃成分を十分に酸化することができず同触媒から未燃成分が流出してしまう。反対に、触媒内の酸素吸蔵量が最大酸素吸蔵量近傍であるときに機関の空燃比がリーンとなると、触媒内で窒素酸化物を十分に還元することができず同触媒から窒素酸化物が流出してしまう。以上のことから、未燃成分や窒素酸化物を効率よく浄化するため触媒内の酸素吸蔵量が所定の量(例えば、最大酸素吸蔵量の略半分の量)になるように機関の空燃比(従って、触媒に流入する排ガスの空燃比)が制御されることが好ましい。
【0004】
このため、例えば、特許文献1に開示された内燃機関の排気浄化装置(空燃比制御装置)は、内燃機関の排気通路に配設された触媒よりも上流側に空燃比センサを介装していて、空燃比センサの出力値と理論空燃比に相当する所定の目標値との偏差(実際には、筒内吸入空気量を空燃比センサ出力値に相当する検出空燃比で除した値である実際の筒内燃料供給量と同筒内吸入空気量を所定の目標値に相当する目標空燃比(従って、理論空燃比)で除した値である目標筒内燃料供給量との偏差)、及び同偏差の時間積分値に基づいて(同偏差を比例・積分処理(PI処理)して)フィードバック制御量を算出し、同フィードバック制御量に基づいて燃料噴射量を補正することで、前記偏差がゼロになるように機関の空燃比をフィードバック制御するようになっている。
【0005】
これによれば、前記偏差の時間積分値(前記積分処理(I処理))に基づくフィードバック制御量の値は、触媒に流入する排ガスの空燃比がリッチとなる前記偏差が発生している間は同偏差を時間積分していくことで同空燃比がリーン側に補正される方向に更新されていき、同触媒に流入する排ガスの空燃比がリーンとなる前記偏差が発生している間は同偏差を時間積分していくことで同空燃比がリッチ側に補正される方向に更新されていく。従って、内燃機関が前記偏差がゼロとなり得る定常運転状態から燃料の過不足(即ち、空気の過不足)が発生し易い(従って、前記偏差が発生し易い)過渡運転状態に一旦移行した後、再び定常運転状態に復帰するまでの期間における触媒に流入する排ガスの空燃比の平均値は目標空燃比(従って、理論空燃比)になることが保障され得る。換言すれば、前記偏差の時間積分値に基づく空燃比フィードバック制御により、触媒に流入する排ガスの空燃比は前記偏差の収支がゼロになるように制御される。従って、上記期間の前後において触媒の酸素吸蔵量は変化することなく、この結果、触媒の酸素吸蔵量が前記所定の量近傍に維持され得る。
【0006】
【特許文献1】
特開平6−280648号公報
【0007】
【発明が解決しようとする課題】
上記開示された装置においては、例えば、急加速運転時等、内燃機関が過渡運転状態にあって触媒に流入する空燃比が或る程度の期間に渡って理論空燃比から相当程度偏移する過度のリーン空燃比になった場合、同空燃比が理論空燃比に復帰した時点における前記偏差の時間積分値に基づくフィードバック制御量の値は、同空燃比をリッチ側に大きく補正する値に更新されている。従って、この場合、触媒に流入する排ガスの空燃比は前記過度のリーン空燃比から理論空燃比に復帰した後、過度なリッチ空燃比に不必要に補正されることになり、この結果、触媒の酸素吸蔵量が「0」となって触媒からの未燃HC,CO等の有害成分(エミッション)の排出量が増大し、また、燃料噴射量の過度の変動に基づく機関出力の変動によりドライバビリティが悪化するという問題がある。
【0008】
従って、本発明の目的は、内燃機関の排気通路に配設された触媒の上流側の空燃比センサの出力値の目標値からの偏差に基づく値の時間積分値に少なくとも基づいて触媒に流入するガスの空燃比をフィードバック制御する内燃機関の排気浄化装置において、同空燃比を不必要に過補正することがないものを提供することにある。
【0009】
【本発明の概要】
本発明の特徴は、内燃機関の排気通路に配設された触媒と、前記触媒よりも上流の前記排気通路に配設された空燃比センサと、少なくとも前記空燃比センサの出力値と所定の目標値との偏差に基づく値の時間積分値に基づいてフィードバック制御量を算出するフィードバック制御量算出手段と、前記フィードバック制御量に基づいて(前記偏差に基づく値がゼロになるように)前記触媒に流入するガスの空燃比をフィードバック制御する空燃比制御手段とを備えた内燃機関の排気浄化装置において、前記フィードバック制御量算出手段が、前記偏差に基づく値の時間積分値が所定の範囲を超える場合、同偏差に基づく値の時間積分値の代わりに同所定の範囲内の値に基づいて前記フィードバック制御量を算出するように構成されたことにある。
【0010】
ここにおいて、前記所定の目標値は理論空燃比に相当する値に設定されることが好適である。また、「空燃比センサの出力値と所定の目標値との偏差に基づく値」は、例えば、空燃比センサの出力値と所定の目標値との偏差、空燃比センサの出力値に相当する検出空燃比(実空燃比)と所定の目標値に相当する目標空燃比との偏差、筒内吸入空気量を空燃比センサの出力値に相当する検出空燃比で除した値である実際の筒内燃料供給量と同筒内吸入空気量を所定の目標値に相当する目標空燃比で除した値である目標筒内燃料供給量との偏差であって、これらに限定されない。
【0011】
また、前記フィードバック制御量に基づいて(前記偏差に基づく値がゼロになるように)前記触媒に流入するガスの空燃比をフィードバック制御する空燃比制御手段は、例えば、機関に供給される混合気の空燃比を制御する手段であってもよく、或いは、同機関に供給される混合気の空燃比の制御に加え、又は同制御とは独立して、触媒上流の排気通路に備えられたノズル等から酸化剤(空気等)又は還元剤(未燃HC等)を所定量供給することで触媒に流入するガスの空燃比を制御する手段であってもよい。なお、機関に供給される混合気の空燃比を制御すれば、触媒に流入するガスの空燃比を制御することができる。
【0012】
これによれば、前記偏差に基づく値の時間積分値に基づく空燃比フィードバック制御により、先に述べたように、触媒に流入するガスの空燃比が前記偏差に基づく値の収支がゼロになるように制御される。従って、内燃機関が前記偏差に基づく値がゼロとなり得る定常運転状態になる毎に触媒の酸素吸蔵量が所定の量(例えば、最大酸素吸蔵量の半分の量)近傍に維持され得るから、未燃成分や窒素酸化物を効率よく浄化することができる。
【0013】
更には、上記排気浄化装置によれば、前記偏差に基づく値の時間積分値が所定の範囲を超える場合、同偏差に基づく値の時間積分値の代わりに同所定の範囲内の値に基づいてフィードバック制御量が算出される。換言すれば、前記偏差に基づく値の時間積分値において前記所定の範囲を超える分の値は前記フィードバック制御量として反映されないように構成されている。
【0014】
従って、例えば、触媒に流入するガスの空燃比が或る程度の期間に渡って理論空燃比から相当程度偏移する過度のリーン空燃比又はリッチ空燃比のいずれか一方の空燃比となることで前記偏差に基づく値の時間積分値が、同空燃比をいずれか他方の空燃比側に補正する方向の値であって前記所定の範囲を超える値となる場合、同偏差に基づく値の時間積分値の代わりに同所定の範囲内の値(例えば、同所定の範囲の上限値又は下限値のいずれか一方の値)に基づいてフィードバック制御量が決定されていく。よって、前記所定の範囲を適切な範囲に設定することにより、触媒に流入する空燃比は前記過度のいずれか一方の空燃比から理論空燃比に復帰した後に過度ないずれか他方の空燃比に不必要に補正されることがなく、この結果、触媒からのエミッションの排出量を少なくし、また、燃料噴射量の変動に基づくドライバビリティの悪化の程度を少なくすることができる。
【0015】
この場合、前記排気浄化装置は、少なくとも前記偏差に基づく値の時間積分値に基づいて前記触媒が吸蔵している酸素量である酸素吸蔵量を推定する酸素吸蔵量推定手段を備え、前記フィードバック制御量算出手段は、前記偏差に基づく値の時間積分値に基づく値としての前記推定された酸素吸蔵量に基づいて前記フィードバック制御量を算出するように構成されることが好適である。
【0016】
前記空燃比センサの出力値と所定の目標値との偏差に基づく値は筒内に供給される燃料の過不足量(従って、触媒に流入するガス中の空気(酸素)の過不足量)に相当する値である。従って、前記偏差に基づく値の時間積分値は触媒に流入するガス中の酸素の過不足量の時間積分値に相当する値である。一方、触媒の酸素吸蔵量の変化量は同触媒に流入するガス中の酸素の過不足量に相当する値となるから触媒の酸素吸蔵量は同触媒に流入するガス中の酸素の過不足量の時間積分値に相当する値となる。以上のことから、触媒の酸素吸蔵量は前記偏差に基づく値の時間積分値に基づいて計算され得る値となる。
【0017】
従って、上記のように、前記偏差に基づく値の時間積分値に基づく値としての前記酸素吸蔵量推定手段により推定された酸素吸蔵量に基づいて空燃比フィードバック制御を実行しても、前記偏差に基づく値の時間積分値に基づく空燃比フィードバック制御の場合と同様、触媒に流入する空燃比が前記偏差に基づく値の収支がゼロになるように制御され得る。従って、内燃機関が前記偏差がゼロとなり得る定常運転状態になる毎に触媒の酸素吸蔵量が前記所定の量(例えば、最大酸素吸蔵量の半分の量)近傍に維持され得ることになる。
【0018】
更には、触媒の酸素吸蔵量は「0」から最大酸素吸蔵量の間で変化する。従って、前記酸素吸蔵量推定手段は「0」から最大酸素吸蔵量の間で変化する値となるように触媒の酸素吸蔵量を推定するから、このように推定される酸素吸蔵量に基づくフィードバック制御量は、「0」から最大酸素吸蔵量の間の幅を有する所定の範囲内の値に基づいて算出され得ることになる。この結果、前記所定の範囲を「0」から最大酸素吸蔵量の間の幅を有する所定の範囲に設定することと同等の効果が得られることになり、触媒に流入するガスの空燃比が前記所定の目標値に相当する目標空燃比から大きく偏移することで前記偏差に基づく値の時間積分値が触媒の酸素吸蔵量が「0」から最大酸素吸蔵量の範囲外の量になることに相当する値になる場合であっても、上記と同様、その後において同空燃比が過度に不必要に補正されることがなくなる。
【0019】
なお、前記酸素吸蔵量推定手段は、「0」以上の所定の下限値から最大酸素吸蔵量以下の所定の上限値の間で変化する値となるように触媒の酸素吸蔵量を推定してもよい。この場合、前記所定の範囲を前記所定の下限値から前記所定の上限値の間の幅を有する所定の範囲に設定することと同等の効果が得られることになる。
【0020】
前記フィードバック制御量算出手段が前記推定された酸素吸蔵量に基づいて前記フィードバック制御量を算出するように構成されている場合、同フィードバック制御量算出手段は、所定の目標酸素吸蔵量と前記推定された酸素吸蔵量との偏差に基づいて(同偏差がゼロになるように)前記フィードバック制御量を算出するように構成されることが更に好適である。この場合、前記所定の目標酸素吸蔵量は最大酸素吸蔵量の略半分の量であることが好ましい。
【0021】
これによれば、内燃機関が定常運転状態になる毎に触媒の酸素吸蔵量が確実に目標酸素吸蔵量に維持され得るようになる。また、所定の目標酸素吸蔵量を最大酸素吸蔵量の半分の量に設定した場合、所定の目標酸素吸蔵量と前記推定された酸素吸蔵量との偏差は、−「最大酸素吸蔵量の半分の値」から「最大酸素吸蔵量の半分の値」の間で変化する値となる。従って、かかる偏差に基づいて算出されるフィードバック制御量における触媒に流入するガスの空燃比をリッチ側に補正する方向の値の最大値と同空燃比をリーン側に補正する方向の値の最大値とが等しくなる。この結果、触媒に流入するガスの空燃比がリッチ方向に過度に補正されることが防止される程度と同空燃比がリーン方向に過度に補正されることが防止される程度とを等しくすることができる。
【0022】
また、前記フィードバック制御量算出手段が前記推定された酸素吸蔵量に基づいて前記フィードバック制御量を算出するように構成されている場合、前記酸素吸蔵量推定手段は、前記触媒内における酸素吸蔵反応の速度を考慮して構築された触媒モデルに基づいて前記酸素吸蔵量を推定するように構成されることが好適である。
【0023】
触媒に流入するガス中の酸素の過不足量が同一であっても、触媒の酸素吸蔵量の変化量は触媒の酸素吸蔵(放出)反応の速度に応じて異なる。従って、触媒の酸素吸蔵量を正確に推定するためには触媒内における酸素吸蔵反応の速度を考慮する必要がある。従って、上記のように構成すれば、前記触媒内における酸素吸蔵反応の速度を考慮して構築された触媒モデルに基づいて前記酸素吸蔵量が推定されるから、より一層正確に触媒の酸素吸蔵量を推定することができ、この結果、より正確に触媒に流入するガスの空燃比をフィードバック制御することができる。
【0024】
また、前記フィードバック制御量算出手段が前記所定の目標酸素吸蔵量と前記推定された酸素吸蔵量との偏差に基づいて前記フィードバック制御量を算出するように構成されている場合であって、前記排気浄化装置が前記内燃機関の運転状態に応じて目標空燃比に相当する前記所定の目標値を設定(変更)する目標値設定手段を備える場合、前記フィードバック制御量算出手段は、更に前記空燃比センサの出力値と前記所定の目標値との偏差に基づく値に基づいて(同偏差がゼロになるように)前記フィードバック制御量を算出するように構成されているとともに、前記所定の目標値が理論空燃比と異なる空燃比に相当する値に設定されている場合、前記所定の目標酸素吸蔵量と前記推定された酸素吸蔵量との偏差の代わりに前記空燃比センサの出力値と同所定の目標値との偏差に基づく値の時間積分値に基づいて前記フィードバック制御量を算出するように構成されることが好適である。
【0025】
即ち、この場合、前記フィードバック制御量算出手段は、前記所定の目標値が理論空燃比に相当する値に設定されている場合、前記空燃比センサの出力値と前記所定の目標値との偏差に基づく値、及び前記所定の目標酸素吸蔵量と前記推定された酸素吸蔵量との偏差に基づいて前記フィードバック制御量を算出するとともに、前記所定の目標値が理論空燃比と異なる空燃比に相当する値に設定されている場合、前記空燃比センサの出力値と前記所定の目標値との偏差に基づく値、及び同偏差に基づく値の時間積分値に基づいて(前記開示した従来の装置と同様に、同偏差に基づく値を比例・積分処理(PI処理)して)フィードバック制御量を算出するように構成される。
【0026】
所定の目標値に相当する目標空燃比が理論空燃比と異なる空燃比となっている場合、内燃機関が定常運転状態にあって触媒に流入するガスの空燃比が目標空燃比に収束していても、触媒に流入するガス中には酸素の過不足量が発生しているから触媒の酸素吸蔵量は変化し、その結果、触媒の酸素吸蔵量は目標酸素吸蔵量に収束し得ない。換言すれば、触媒に流入するガスの空燃比を目標空燃比に収束させる制御と触媒の酸素吸蔵量を目標酸素吸蔵量に収束させる制御とを同時に行うと互いの制御目標が相違し、両制御が互いに干渉することになるから良好な空燃比フィードバック制御を行うことができない。これに対し、上記のように構成すれば、所定の目標値に相当する目標空燃比が理論空燃比と異なる空燃比となっている場合、前記偏差に基づく値を比例・積分処理(PI処理)してフィードバック制御量が算出されるように構成されるから、制御目標が前記偏差に基づく値をゼロに収束させること(従って、触媒に流入するガスの空燃比を目標空燃比に収束させること)のみとなって前述した制御の干渉の発生を回避することができる。
【0027】
また、前記何れかの排気浄化装置が適用される内燃機関が指示に応じて燃料噴射を行う燃料噴射手段を備えている場合であって、前記空燃比制御手段が、前記内燃機関に供給される混合気の空燃比を前記所定の目標値に相当する目標空燃比とするために必要な基本燃料噴射量を同内燃機関の運転状態に応じて決定するとともに同基本燃料噴射量を前記フィードバック制御量に基づいて補正することにより得られる燃料噴射量の燃料の噴射指示を前記燃料噴射手段に対して行うことで前記触媒に流入するガスの空燃比をフィードバック制御するように構成されている場合、前記何れかの排気浄化装置は、更に、前記フィードバック制御量算出手段が前記空燃比センサの出力値と前記所定の目標値との偏差に基づく値の時間積分値が前記所定の範囲を超えるか否かに拘わらず同偏差に基づく値の時間積分値に基づいて前記フィードバック制御量を算出した場合であって、且つ前記内燃機関が定常運転状態にある場合における同偏差に基づく値の時間積分値に基づいて前記基本燃料噴射量を補正する基本燃料噴射量補正手段を備えることが好適である。
【0028】
前記排気浄化装置において、燃料噴射手段(例えば、インジェクタ)の誤差(指示燃料噴射量と実際の燃料噴射量の差)、(指示)燃料噴射量を計算するために使用される吸入空気流量センサ(例えば、エアフローメータ)の誤差(吸入空気流量計測値と実際の吸入空気流量の差)等の誤差が発生している場合、内燃機関が定常運転状態にあって前記偏差に基づく値がゼロになっているとき(従って、機関の空燃比(従って、触媒に流入するガスの空燃比)が所定の目標値に相当する目標空燃比に収束しているとき)における前記フィードバック制御量(従って、同偏差に基づく値の時間積分値)は前記誤差の大きさに応じた値となる。換言すれば、かかる誤差は、前記偏差に基づく値の時間積分値が同誤差の大きさに応じた値に収束することにより吸収され得る。
【0029】
しかしながら、前述のごとく、前記偏差に基づく値の時間積分値が所定の範囲を超える場合に同偏差に基づく値の時間積分値の代わりに同所定の範囲内の値に基づいて前記フィードバック制御量を算出するように構成すると、定常運転状態において前記偏差に基づく値の時間積分値が前記所定の範囲外の値になることに相当する大きな前記誤差が発生したとき、かかる誤差の総てを吸収することができず、前記偏差に基づく値がゼロに収束し得なくなる事態が発生する。
【0030】
これに対し、上記のような基本燃料噴射量補正手段を設けることにより、前記誤差の大きさに拘わらず、同基本燃料噴射量補正手段が基本燃料噴射量を補正する際に使用する定常運転状態における前記偏差に基づく値の時間積分値は正確に同誤差の大きさに応じた値となっている。従って、前記偏差に基づく値の時間積分値に基づいて前記誤差の大きさに応じた分だけ正確に基本燃料噴射量を補正することができる。換言すれば、前記偏差に基づく値の時間積分値が前記所定の範囲外の値になることに相当する大きな前記誤差が発生している場合であっても、前記誤差が基本燃料噴射量の補正により確実に補償され得る。
【0031】
【発明の実施の形態】
以下、本発明による内燃機関の排気浄化装置を含む空燃比制御装置の実施形態について図面を参照しつつ説明する。
【0032】
図1は、本発明の実施形態による空燃比制御装置を火花点火式多気筒(4気筒)内燃機関10に適用したシステムの概略構成を示している。この内燃機関10は、シリンダブロック、シリンダブロックロワーケース、及びオイルパン等を含むシリンダブロック部20と、シリンダブロック部20の上に固定されるシリンダヘッド部30と、シリンダブロック部20にガソリン混合気を供給するための吸気系統40と、シリンダブロック部20からの排ガスを外部に放出するための排気系統50とを含んでいる。
【0033】
シリンダブロック部20は、シリンダ21、ピストン22、コンロッド23、及びクランク軸24を含んでいる。ピストン22はシリンダ21内を往復動し、ピストン22の往復動がコンロッド23を介してクランク軸24に伝達され、これにより同クランク軸24が回転するようになっている。シリンダ21とピストン22のヘッドは、シリンダヘッド部30とともに燃焼室25を形成している。
【0034】
シリンダヘッド部30は、燃焼室25に連通した吸気ポート31、吸気ポート31を開閉する吸気弁32、吸気弁32を駆動するインテークカムシャフトを含むとともに同インテークカムシャフトの位相角を連続的に変更する可変吸気タイミング装置33、可変吸気タイミング装置33のアクチュエータ33a、燃焼室25に連通した排気ポート34、排気ポート34を開閉する排気弁35、排気弁35を駆動するエキゾーストカムシャフト36、点火プラグ37、点火プラグ37に与える高電圧を発生するイグニッションコイルを含むイグナイタ38、及び燃料を吸気ポート31内に噴射するインジェクタ(燃料噴射手段)39を備えている。
【0035】
吸気系統40は、吸気ポート31に連通し同吸気ポート31とともに吸気通路を形成するインテークマニホールドを含む吸気管41、吸気管41の端部に設けられたエアフィルタ42、吸気管41内にあって吸気通路の開口断面積を可変とするスロットル弁43、スロットル弁駆動手段を構成するDCモータからなるスロットル弁アクチュエータ43a、スワールコントロールバルブ(以下、「SCV」と称呼する。)44、及びDCモータからなるSCVアクチュエータ44aを備えている。
【0036】
排気系統50は、排気ポート34に連通したエキゾーストマニホールド51、エキゾーストマニホールド51(実際には、各排気ポート34に連通した各々のエキゾーストマニホールド51が集合した集合部)に接続されたエキゾーストパイプ(排気管)52、エキゾーストパイプ52に配設(介装)された上流側の三元触媒53(上流側触媒コンバータ、又はスタート・キャタリティック・コンバータとも云うが、以下「第1触媒53」と称呼する。)、及びこの第1触媒53の下流のエキゾーストパイプ52に配設(介装)された下流側の三元触媒54(車両のフロア下方に配設されるため、アンダ・フロア・キャタリティック・コンバータとも云うが、以下「第2触媒54」と称呼する。)を備えている。排気ポート34、エキゾーストマニホールド51、及びエキゾーストパイプ52は、排気通路を構成している。
【0037】
一方、このシステムは、熱線式エアフローメータ61、スロットルポジションセンサ62、カムポジションセンサ63、クランクポジションセンサ64、水温センサ65、第1触媒53の上流の排気通路(本例では、上記各々のエキゾーストマニホールド51が集合した集合部)に配設された空燃比センサ66(以下、「上流側空燃比センサ66」と称呼する。)、第1触媒53の下流であって第2触媒54の上流の排気通路に配設された空燃比センサ67(以下、「下流側空燃比センサ67」と称呼する。)、及びアクセル開度センサ68を備えている。
【0038】
熱線式エアフローメータ61は、吸気管41内を流れる吸入空気の単位時間あたりの質量流量に応じた電圧Vgを出力するようになっている。かかるエアフローメータ61の出力Vgと、計測された吸入空気量(流量)Gaとの関係は、図2に示したとおりである。スロットルポジションセンサ62は、スロットル弁43の開度を検出し、スロットル弁開度TAを表す信号を出力するようになっている。カムポジションセンサ63は、インテークカムシャフトが90°回転する毎に(即ち、クランク軸24が180°回転する毎に)一つのパルスを有する信号(G2信号)を発生するようになっている。クランクポジションセンサ64は、クランク軸24が10°回転する毎に幅狭のパルスを有するとともに同クランク軸24が360°回転する毎に幅広のパルスを有する信号を出力するようになっている。この信号は、エンジン回転速度NEを表す。水温センサ65は、内燃機関10の冷却水の温度を検出し、冷却水温THWを表す信号を出力するようになっている。
【0039】
上流側空燃比センサ66は、限界電流式の酸素濃度センサであり、図3に示したように、空燃比A/Fに応じた電流を出力し、この電流に応じた電圧である出力値vabyfsを出力するようになっていて、特に、空燃比が理論空燃比であるときには出力値vabyfsは上流側目標値vstoichになる。図3から明らかなように、上流側空燃比センサ66によれば、広範囲にわたる空燃比A/Fを精度良く検出することができる。
【0040】
下流側空燃比センサ67は、起電力式(濃淡電池式)の酸素濃度センサであり、図4に示したように、理論空燃比近傍において急変する電圧である出力値Voxsを出力するようになっている。より具体的に述べると、下流側空燃比センサ67は、空燃比が理論空燃比よりもリーンのときは略0.1(V)、空燃比が理論空燃比よりもリッチのときは略0.9(V)、及び空燃比が理論空燃比のときは0.5(V)の電圧を出力するようになっている。アクセル開度センサ68は、運転者によって操作されるアクセルペダル81の操作量を検出し、同アクセルペダル81の操作量Accpを表す信号を出力するようになっている。
【0041】
電気制御装置70は、互いにバスで接続されたCPU71、CPU71が実行するルーチン(プログラム)、テーブル(ルックアップテーブル、マップ)、及び定数等を予め記憶したROM72、CPU71が必要に応じてデータを一時的に格納するRAM73、電源が投入された状態でデータを格納するとともに同格納したデータを電源が遮断されている間も保持するバックアップRAM74、並びにADコンバータを含むインターフェース75等からなるマイクロコンピュータである。インターフェース75は、前記センサ61〜68と接続され、CPU71にセンサ61〜68からの信号を供給するとともに、同CPU71の指示に応じて可変吸気タイミング装置33のアクチュエータ33a、イグナイタ38、インジェクタ39、スロットル弁アクチュエータ43a、及びSCVアクチュエータ44aに駆動信号を送出するようになっている。
【0042】
(空燃比フィードバック制御の概要)
次に、上記のように構成された排気浄化装置を含んだ空燃比制御装置が行う機関の空燃比のフィードバック制御の概要について説明する。
【0043】
第1触媒53(第2触媒54も同様である。)は、同第1触媒53に流入するガスの空燃比が理論空燃比であるときに、HC,COを酸化するとともにNOxを還元し、これらの有害成分を高い効率で浄化する。また、第1触媒53は、酸素を吸蔵・放出する機能(酸素吸蔵機能、酸素吸蔵・放出機能)を有し、この酸素吸蔵・放出機能により、空燃比が理論空燃比からある程度まで偏移したとしても、HC,CO、及びNOxを浄化することができる。即ち、機関の空燃比がリーンとなって第1触媒53に流入するガスにNOxが多量に含まれると、第1触媒53はNOxから酸素分子を奪って同酸素分子を吸蔵するとともに同NOxを還元し、これによりNOxを浄化する。また、機関の空燃比がリッチになって第1触媒53に流入するガスにHC,COが多量に含まれると、第1触媒53はこれらに吸蔵している酸素分子を与えて(放出して)酸化し、これによりHC,COを浄化する。換言すれば、第1触媒53に流入する排ガスの空燃比がリッチである場合には同第1触媒53の酸素吸蔵量が減少するとともに第1触媒53に流入する排ガスの空燃比がリーンである場合には同第1触媒53の酸素吸蔵量が増加する。
【0044】
従って、第1触媒53が連続的に流入するリッチ空燃比の排ガス中にある多量のHC,COを効率的に浄化するためには、同第1触媒53が酸素を多量に貯蔵していなければならず、逆に連続的に流入するリーン空燃比の排ガス中にある多量のNOxを効率的に浄化するためには、同第1触媒53が酸素を十分に貯蔵し得る状態になければならないことになる。以上のことから、HC,COやNOxを効率よく浄化するため第1触媒53内の酸素吸蔵量が最大酸素吸蔵量の半分の量(目標酸素吸蔵量)になるように機関の空燃比(従って、第1触媒53に流入する排ガスの空燃比)が制御されることが好ましい。
【0045】
第1触媒53内の酸素吸蔵量が最大酸素吸蔵量の半分の量に維持され得るように制御するためには、少なくとも、第1触媒53に流入する排ガスの空燃比の平均値を理論空燃比と等しくするように制御すること(即ち、第1触媒53に流入する排ガスの空燃比の理論空燃比からの偏差の収支をゼロにするように制御すること)が要求される。このため、上流側空燃比センサ66の出力値vabyfsと同出力値の目標値である上流側目標値(所定の目標値)vstoichとの偏差の時間積分値に基づいてフィードバック制御量を求め、同フィードバック制御量に基づいて同偏差が「0」になるように第1触媒53に流入する排ガスの空燃比(従って、機関10に供給される混合気の空燃比)をフィードバック制御する必要がある。
【0046】
一方、上流側空燃比センサ66の出力値vabyfsと上流側目標値vstoichとの偏差は、シリンダ内(筒内)に供給される燃料の過不足量(従って、第1触媒53に流入するガス中の空気(酸素)の過不足量)に相当する値である。一方、第1触媒53の酸素吸蔵量の変化量も同第1触媒53に流入する排ガス中の酸素の過不足量に相当する値となる。以上のことから、第1触媒53の酸素吸蔵量の変化量は前記偏差に相当する値となるから、第1触媒53の酸素吸蔵量は前記偏差の時間積分値に基づく値となる。
【0047】
そこで、本実施形態の空燃比制御装置(排気浄化装置、以下、単に「本装置」と云うこともある。)は、後述する触媒モデルを使用して上流側空燃比センサ66の出力値vabyfsと上流側目標値vstoichとの偏差に基づく値である筒内燃料供給量偏差DFcの時間積分値に基づいて第1触媒53の酸素吸蔵量OSAallを逐次求める。また、本装置は、原則的に、下記数1に従って、前記筒内燃料供給量偏差DFcに比例ゲインGpを乗じた比例項(P項)DFipと、酸素吸蔵量OSAallから最大酸素吸蔵量Cmaxallの半分の量を減じた値に積分ゲインGiを乗じた積分項(I項)DFiiとの和であるフィードバック制御量DFiを求める。ここで、筒内燃料供給量偏差DFcは、後述する筒内吸入空気量Mcを上流側空燃比センサ66の出力値vabyfsに相当する検出空燃比abyfsで除した値である実際の筒内燃料供給量Fcと同筒内吸入空気量Mcを目標空燃比abyfrで除した値である目標筒内燃料供給量Fcrとの偏差である。また、下記数1において、Cmaxallは第1触媒53の最大酸素吸蔵量であって、この最大酸素吸蔵量Cmaxallの値は、別途、後述するように求められる。
【0048】
【数1】
DFi=Gp・DFc+Gi・(OSAall−(Cmaxall/2))
【0049】
そして、本装置は、内燃機関10の運転状態に応じて決定される後述する(補正後)基本燃料噴射量Fbaseに前記フィードバック制御量DFiを加えた量の燃料を最終燃料噴射量Fiとして設定し、同最終燃料噴射量Fiの燃料を吸気行程を迎える気筒のインジェクタ39から噴射することにより機関10に供給される混合気の空燃比をフィードバック制御する。この結果、第1触媒53に流入する排ガスの空燃比の理論空燃比からの偏差の収支がゼロになるように制御されるとともに、筒内燃料供給量偏差DFcが「0」になるように(従って、上流側空燃比センサ66の出力値vabyfsが上流側目標値vstoichになるように)、且つ、第1触媒53の酸素吸蔵量OSAallが最大酸素吸蔵量Cmaxallの半分の量になるようにそれぞれ制御される。
【0050】
(空燃比フィードバック制御の具体的内容)
前述した空燃比フィードバック制御についてより具体的に述べると、この空燃比制御装置は、機能ブロック図である図5に示したように、A1〜A16の各手段等を含んで構成されている。以下、図5を参照しながら各手段等について説明していく。
【0051】
<(補正前)基本燃料噴射量の算出>
先ず、筒内吸入空気量算出手段A1は、エアフローメータ61が計測している吸入空気流量Gaと、クランクポジションセンサ64の出力に基づいて得られるエンジン回転速度NEと、ROM72が記憶しているテーブルMAPMcとに基づき今回の吸気行程を迎える気筒の吸入空気量である筒内吸入空気量Mc(k)を求める。ここで、添え字の(k)は、今回の吸気行程に対する値であることを示している(以下、他の物理量についても同様。)。筒内吸入空気量Mcは、各気筒の吸気行程に対応されながらRAM73に記憶されていく。
【0052】
目標値設定手段としての上流側目標空燃比設定手段A2は、内燃機関10の運転状態であるエンジン回転速度NE、及びスロットル弁開度TA等に基づいて上流側目標値に相当する上流側目標空燃比abyfr(k)を決定する。この上流側目標空燃比abyfr(k)は、例えば、内燃機関10の暖機運転中、スロットル弁開度TAが所定値以上である急加速運転中等の特殊な場合を除き理論空燃比に設定されている。また、上流側目標空燃比abyfrは、各気筒の吸気行程に対応されながらRAM73に記憶されていく。
【0053】
補正前基本燃料噴射量算出手段A3は、筒内吸入空気量算出手段A1により求められた筒内吸入空気量Mc(k)を上流側目標空燃比設定手段A2により設定された上流側目標空燃比abyfr(k)で除することにより、機関の空燃比を同上流側目標空燃比abyfr(k)とするための今回の吸気行程に対する補正前基本燃料噴射量Fbasebを求める。
【0054】
このように、本装置は、筒内吸入空気量算出手段A1、上流側目標空燃比設定手段A2、及び補正前基本燃料噴射量算出手段A3を利用して、補正前基本燃料噴射量Fbasebを求める。
【0055】
<最終燃料噴射量の算出>
先ず、補正後基本燃料噴射量算出手段A4は、補正前基本燃料噴射量算出手段A3により求められた補正前基本燃料噴射量Fbasebに後述する基本燃料噴射量補正量取得手段A16による最新の基本燃料噴射量補正量算出テーブルΔFbase(Ga)に従って求められる基本燃料噴射量補正量ΔFbaseを加算することで補正後基本燃料噴射量Fbase(即ち、目標筒内燃料供給量Fcr(k))を求める。目標筒内燃料供給量Fcrは、各気筒の吸気行程に対応されながらRAM73に記憶されていく。
【0056】
最終燃料噴射量算出手段A5は、前記補正後基本燃料噴射量Fbaseに前述したフィードバック制御量DFiを加えることで、下記数2に基いて最終燃料噴射量Fiを求める。本装置は、このようにして、補正後基本燃料噴射量算出手段A4、及び最終燃料噴射量算出手段A5により補正後基本燃料噴射量Fbaseをフィードバック制御量DFiに基づいて補正することにより最終燃料噴射量Fiを求め、同最終燃料噴射量Fiの燃料を今回の吸気行程を迎える気筒に対してインジェクタ39により噴射する。このようにして、機関10に供給される混合気の空燃比(従って、第1触媒53に流入する排ガスの空燃比)をフィードバック制御量DFiに基づいてフィードバック制御する手段が空燃比制御手段に相当する。
【0057】
【数2】
Fi=Fbase+DFi
【0058】
<フィードバック制御量の算出>
先ず、テーブル変換手段A6は、上流側空燃比センサ66の出力値vabyfsと、先に説明した図3に示した上流側空燃比センサ出力値vabyfsと空燃比A/Fとの関係を規定したテーブルとに基づいて、上流側空燃比センサ66が検出する現時点における検出空燃比abyfsを求める。
【0059】
筒内吸入空気量遅延手段A7は、筒内吸入空気量算出手段A1により吸気行程毎に求められRAM73に記憶されている筒内吸入空気量Mcのうち、現時点からNストローク(N回の吸気行程)前に吸気行程を迎えた気筒の筒内吸入空気量McをRAM73から読み出し、これを筒内吸入空気量Mc(k−N)として設定する。
【0060】
筒内燃料供給量算出手段A8は、筒内吸入空気量遅延手段A7により求められた現時点からNストローク前の筒内吸入空気量Mc(k−N)をテーブル変換手段A6により求められた現時点における検出空燃比abyfsで除することで、現時点からNストローク前の実際の筒内燃料供給量Fc(k−N)を求める。ここで、前記値Nは、内燃機関10の排気量、及び燃料室25から上流側空燃比センサ66までの距離等により異なる値である。
【0061】
このように、現時点からNストローク前の実際の筒内燃料供給量Fc(k−N)を求めるために、現時点からNストローク前の筒内吸入吸気量Mc(k−N)を現時点における検出空燃比abyfsで除するのは、燃焼室25内で燃料された混合気が上流側空燃比センサ66に到達するまでには、Nストロークに相当する時間Lを要しているからである。
【0062】
目標筒内燃料供給量遅延手段A9は、補正後基本燃料噴射量算出手段A4により吸気行程毎に求められRAM73に記憶されている目標筒内燃料供給量Fcrのうち、現時点からNストローク前の目標筒内燃料供給量FcrをRAM73から読み出し、これを目標筒内燃料供給量Fcr(k−N)として設定する。
【0063】
筒内燃料供給量偏差算出手段A10は、下記数3に基づいて、目標筒内燃料供給量遅延手段A9により設定された現時点からNストローク前の目標筒内燃料供給量Fcr(k−N)から筒内燃料供給量算出手段A8により求められた現時点からNストローク前の実際の筒内燃料供給量Fc(k−N)を減じることにより、前述した筒内燃料供給量偏差DFcを求める。この筒内燃料供給量偏差DFcは、Nストローク前の時点で筒内に供給された燃料の過不足分を表す量であって、上流側空燃比センサ66の出力値vabyfsと上流側目標値(上流側目標空燃比abyfrが理論空燃比のときは図3に示すvstoich)との偏差に基づく値である。
【0064】
【数3】
DFc=Fcr(k−N)−Fc(k−N)
【0065】
比例項算出手段A11は、前記筒内燃料供給量偏差DFcを比例処理(P処理)することで、下記数4に基づいて比例項(P項)に基づくフィードバック制御量DFipを求める。下記数4において、Gpは予め設定された比例ゲイン(比例定数(正の値))である。
【0066】
【数4】
DFip=Gp・DFc
【0067】
触媒モデルA12は、前記筒内燃料供給量偏差DFcの時間積分値、後述する最大酸素吸蔵量取得手段A15による求められている第1触媒53の最大酸素吸蔵量Cmaxallの値等に基づいて第1触媒53の酸素吸蔵量OSAallを求めるモデルであるが、これについては後述する。
【0068】
積分項算出手段A13は、触媒モデルA12により求められた第1触媒53の酸素吸蔵量OSAallの値と前記第1触媒53の最大酸素吸蔵量Cmaxallの値とから下記数5に基づいて積分項(I項)に基づくフィードバック制御量DFiiを求める。下記数5において、Giは予め設定された積分ゲイン(積分定数(正の値))である。
【0069】
【数5】
DFii=Gi・(OSAall−(Cmaxall/2))
【0070】
フィードバック制御量算出手段A14は、前記比例項に基づくフィードバック制御量DFipに前記積分項に基づくフィードバック制御量DFiiを加えることで、前述のフィードバック制御量DFiを求める。このフィードバック制御量DFiは、先に述べたように最終燃料噴射量算出手段A5により最終燃料噴射量Fiを求める際に使用される。
【0071】
このようにして、本装置は、上流側空燃比センサ66の出力値vabyfsと上流側目標空燃比abyfrに対応する上流側目標値との偏差に基づく値である筒内燃料供給量偏差DFc、及び、同筒内燃料供給量偏差DFcの時間積分値に基づく値としての第1触媒53の酸素吸蔵量OSAallと目標酸素吸蔵量(第1触媒53の最大酸素吸蔵量Cmaxallの半分の量)との偏差に基づいてフィードバック制御量DFiを求める。そして、補正後基本燃料噴射量Fbaseにフィードバック制御量DFiを加えることで、同補正後基本燃料噴射量Fbaseを補正する。
【0072】
例えば、機関の空燃比が急変してリーンとなると、テーブル変換手段A6にて求められる検出空燃比abyfsは上流側目標空燃比設定手段A2により設定されている上流側目標空燃比abyfrよりもリーンな値(より大きな値)として求められる。このため、筒内燃料供給量算出手段A8にて求められる実際の筒内燃料供給量Fc(k−N)は目標筒内燃料供給量遅延手段A9にて求められる目標筒内燃料供給量Fcr(k−N)よりも小さい値となり、筒内燃料供給量偏差DFcは大きい正の値として求められる。また、筒内燃料供給量偏差DFcの値が正の値であるとき(第1触媒53に流入する排ガスの空燃比がリーンであるとき)第1触媒53内の酸素吸蔵量OSAallは増加するように算出されていく。従って、比例項に基づくフィードバック制御量DFipも積分項に基づくフィードバック制御量DFiiも正の値となってフィードバック制御量DFiが大きい正の値となる。これにより、最終燃料噴射量算出手段A5にて求められる燃料噴射量Fiは、補正後基本燃料噴射量Fbaseよりも大きくなって、機関の空燃比がリッチとなるように制御される。
【0073】
反対に、機関の空燃比が急変してリッチとなると、検出空燃比abyfsは上流側目標空燃比abyfrよりもリッチな値(より小さな値)として求められる。このため、実際の筒内燃料供給量Fc(k−N)は目標筒内燃料供給量Fcr(k−N)よりも大きい値となり、筒内燃料供給量偏差DFcは負の値として求められる。また、筒内燃料供給量偏差DFcの値が負の値であるとき(第1触媒53に流入する排ガスの空燃比がリッチであるとき)第1触媒53内の酸素吸蔵量OSAallは減少するように算出されていく。従って、比例項に基づくフィードバック制御量DFipも積分項に基づくフィードバック制御量DFiiも負の値となってフィードバック制御量DFiが負の値となる。これにより、燃料噴射量Fiは、補正後基本燃料噴射量Fbaseよりも小さくなって、機関の空燃比がリーンとなるように制御される。以上、テーブル変換手段A6、筒内吸入空気量遅延手段A7、筒内燃料供給量算出手段A8、目標筒内燃料供給量遅延手段A9、筒内燃料供給量偏差算出手段A10、比例項算出手段A11、触媒モデルA12、積分項算出手段A13、及びフィードバック制御量算出手段A14はフィードバック制御量算出手段に相当する。
【0074】
図6は、本装置を適用した機関10の空燃比(従って、第1触媒53に流入する排ガスの空燃比)が一時的に理論空燃比から大きく偏移する場合における第1触媒53の酸素吸蔵量OSAall、積分項に基づくフィードバック制御量DFii(以下、「積分項DFii」と云うこともある。)の変化を実線にて示したタイムチャートである。また、図6は、比較例として、図7に機能ブロック図にて示した前記開示した従来の排気浄化装置(空燃比制御装置)を内燃機関10に適用した場合の変化も破線にて示している。なお、図7において図5と同一又は類似の手段等については図5の符号と同一の符号を付している。
【0075】
従来の装置は、フィードバック制御量DFiの算出のみが本装置と相違していて、従来の装置ではフィードバック制御量DFiはPIコントローラA11により下記数6に従って算出される。下記数6において、Gpは上記数1におけるものと同一の比例ゲインであり、Gi’は積分ゲインである。即ち、従来の装置において、比例項に基づくフィードバック制御量DFipは上記数1におけるものと同一のGp・DFcであり、積分項に基づくフィードバック制御量(積分項)DFiiはGi’・ΣDFcである。
【0076】
【数6】
DFi=Gp・DFc+Gi’・ΣDFc
【0077】
図6に示すように、時刻t0から時刻t1の間、(a)に示すように運転者によるアクセルペダル81の操作量Accpが所定の一定値に維持されていることで内燃機関10が定常運転状態にあり、(b)に示すように機関の空燃比abyfsが理論空燃比近傍に維持され、(c)に示すように第1触媒53の酸素吸蔵量OSAallが目標酸素吸蔵量(Cmaxall/2)の近傍に維持されていて、(d)に示すように積分項(I項)が略「0」に維持されているものとする。
【0078】
この状態から、時刻t1においてアクセルペダル81の操作量Accpが急に最大値まで増加したものとすると、エアフローメータ61の遅れによる吸入空気流量Gaの計測誤差等により実際の吸入空気量に対する燃料噴射量Fiが不足して機関の空燃比abyfsが時刻t1以降、一時的に相当のリーン空燃比になる。この時刻t1以降の変化について、先ず、従来の装置を適用した場合について説明する。なお、説明の便宜上、比例項に基づくフィードバック制御量DFipの作用についての説明は省略する。
【0079】
第1触媒53の酸素吸蔵量OSAallは、第1触媒53に流入する排ガスの空燃比(従って、機関の空燃比abyfs)がリーンのとき増加し、機関の空燃比abyfsがリッチのとき減少する。また、積分項(Gi’・ΣDFc)(積分項(Gi・(OSAall−(Cmaxall/2)))についても同様である。)も、機関の空燃比abyfsがリーンのとき筒内燃料供給量偏差DFcが正の値となることから増加し、機関の空燃比abyfsがリッチのとき筒内燃料供給量偏差DFcが負の値となることから減少する。
【0080】
従って、時刻t1以降、機関の空燃比abyfsがリーンとなっているから酸素吸蔵量OSAallは増加していき、時刻t2になると最大酸素吸蔵量Cmaxallに達するとともに時刻t2以降Cmaxallに維持される。一方、積分項(Gi’・ΣDFc)の値も時刻t1以降増加していき、更に時刻t2以降も機関の空燃比abyfsがリーンとなっているから増加し続ける。これにより、積分項の値は正の値となって機関の空燃比abyfsをリッチ側に補正する値になっていくから、機関の空燃比abyfsは或る時点から減少し始め、時刻t3になるとリッチに変化する。この時刻t3まで積分項の値は増加し続けるから同時刻t3の時点において積分項の値は相当に大きい正の値(最大値)となって機関の空燃比abyfsをリッチ側に大きく補正する値になっている。
【0081】
従って、時刻t3以降、機関の空燃比abyfsは不必要に大きくリッチ側に補正される。この結果、時刻t3以降、機関の空燃比abyfsが相当のリッチとなることで最大酸素吸蔵量Cmaxallから減少を開始した酸素吸蔵量OSAallは時刻t5になると「0」に達するとともに時刻t5以降「0」に維持される。一方、積分項の値も時刻t3以降、機関の空燃比abyfsが再びリーンに変化する時刻t7までの間、空燃比abyfsがリッチとなっているから減少し続ける。これにより、積分項の値は時刻t7の時点でも相当に絶対値が大きい負の値となって機関の空燃比abyfsをリーン側に過度に大きく補正する値になっている。
【0082】
よって、機関の空燃比abyfsは時刻t7以降もリーン側に或る程度大きく補正され、再びリッチ側に変化する時刻t9の時点以降にて理論空燃比近傍に収束する。このように、従来の装置を適用した場合、積分項(Gi’・ΣDFc)の値に制限がないことから、機関の空燃比abyfs(即ち、最終燃料噴射量Fi)は同積分項の作用により不必要に過度に補正されることで振動的となって、第1触媒53から多量のエミッション(時刻t2〜時刻t3の間はNOx、時刻t5〜時刻t7の間はHC,CO)が排出されるとともに、ドライバビリティが悪化する。
【0083】
これに対し、本装置を適用した場合、積分項(Gi・(OSAall−(Cmaxall/2)))の値は酸素吸蔵量OSAallと同様に変化し、酸素吸蔵量OSAallが最大酸素吸蔵量Cmaxallに達する時刻t2以降、積分項の値は最大値(Cmaxall/2)に維持されて同最大値より大きい値になることはない。従って、機関の空燃比abyfsがリーンからリッチに変化する時刻t4の時点において積分項の値(Cmaxall/2)は機関の空燃比abyfsをリッチ側に補正する値になっているもののリッチ側に過度に大きく補正する値とはなっていない。
【0084】
従って、機関の空燃比abyfsは時刻t4以降、リーン側に過度に補正されることはなく、再びリーン側に変化する時刻t6の時点(又は、再びリッチ側に変化する時刻t8の時点)以降にて理論空燃比近傍に収束する。また、酸素吸蔵量OSAallは、最大酸素吸蔵量Cmaxallとなっている時刻t4以降「0」にまで達することはなく、時刻t8の時点でほぼ目標酸素吸蔵量(Cmax/2)に収束する。更に、積分項の値も、時刻t4以降大きく変動することはなく、時刻t8の時点でほぼ「0」に収束する。
【0085】
このように、本装置を適用した場合、積分項(Gi・(OSAall−(Cmaxall/2)))の値に制限(−Gi・(Cmaxall/2)≦積分項≦Gi・(Cmaxall/2))があるから、機関の空燃比abyfs(即ち、最終燃料噴射量Fi)は同積分項の作用により不必要に過度に補正されることがなく(即ち、図6(d)に斜線にて示した面積に相当する分だけ積分項によるリッチ側への補正量が少ないから)、比較的速やかに収束する。従って、第1触媒53からのエミッションの排出量が少なく(時刻t2〜時刻t4の間にのみNOxが多量に排出される。)、また、ドライバビリティが悪化の程度が少なくなる。
【0086】
以上、本装置のように積分項(Gi・(OSAall−(Cmaxall/2)))の値に基づいてフィードバック制御量DFiを算出することは、上流側空燃比センサ66の出力値vabyfsと上流側目標空燃比abyfrに対応する上流側目標値との偏差に基づく値である筒内燃料供給量偏差DFcの時間積分値(に基づく値)である積分項(Gi’・ΣDFc)の値に基づいてフィードバック制御量DFiを算出するとともに同積分項(Gi’・ΣDFc)の値が所定の範囲(−Gi・(Cmaxall/2)≦積分項≦Gi・(Cmaxall/2))を超える場合、同積分項(Gi’・ΣDFc)の代わりに同所定の範囲の上限値Gi・(Cmaxall/2)又は下限値−Gi・(Cmaxall/2)のいずれか一方に基づいてフィードバック制御量DFiを算出することと実質的に等価である。
【0087】
(目標空燃比が理論空燃比以外のときの処置)
上記数1に示したように、本装置においては、比例項DFip(=Gp・DFc)に基づいて筒内燃料供給量偏差DFcが「0」になるように(従って、上流側空燃比センサ66の検出空燃比abyfsが上流側目標空燃比abyfrになるように)、且つ、積分項DFii(=Gi・(OSAall−(Cmaxall/2)))に基づいて第1触媒53の酸素吸蔵量OSAallが目標酸素吸蔵量(Cmaxall/2)になるようにそれぞれ制御される。
【0088】
しかしながら、先に述べたように、目標空燃比abyfrが理論空燃比と異なる空燃比となる場合があり、この場合、内燃機関10が定常運転状態にあって第1触媒53に流入するガスの前記検出空燃比abyfsが目標空燃比abyfrに収束していても、第1触媒53に流入するガス中には酸素の過不足量が発生しているから第1触媒53の酸素吸蔵量OSAallは変化し、同酸素吸蔵量OSAallは目標酸素吸蔵量(Cmaxall/2)に収束し得ない。換言すれば、比例項DFipに基づく制御と積分項DFiiに基づく制御とを同時に行うと互いの制御目標が相違し、両制御が互いに干渉する。
【0089】
そこで、本装置は、目標空燃比abyfrが理論空燃比と異なる空燃比となる場合、上記数1の代わりに上記数6に従って、フィードバック制御量DFiを算出する。即ち、目標空燃比abyfrが理論空燃比と異なる空燃比となる場合、本装置は、図5に示した機能ブロック図の代わりに図7に示した前記開示した従来の排気浄化装置の機能ブロック図に基づいてフィードバック制御量DFiを算出する。これにより、制御目標が筒内燃料供給量偏差DFcをゼロに収束させること(従って、第1触媒53に流入するガスの検出空燃比abyfsを目標空燃比abyfrに収束させること)のみとなって前述した制御の干渉の発生を回避することができる。
【0090】
(基本燃料噴射量の補正)
前述した空燃比制御装置において、例えば、インジェクタ39の誤差((指示)燃料噴射量Fiと実際の燃料噴射量の差)、(指示)燃料噴射量Fiを計算するために使用されるエアフローメータ61の誤差(吸入空気流量計測値Gaと実際の吸入空気流量の差)等の誤差が発生している場合((指示)燃料噴射量Fiの値が目標空燃比abyfrを得るための適切な値と異なる値となる場合)がある。この場合、内燃機関10が定常運転状態にあって筒内燃料供給量偏差DFcが「0」になっているとき(従って、機関10の検出空燃比abyfsが目標空燃比abyfrに収束しているとき)における上記数1に基づくフィードバック制御量DFiの値(即ち、比例項DFipの値が「0」になっていることから積分項DFiiの値そのものの値)は、「0」とはならず、前記誤差の大きさに応じた値に収束する。従って、前記誤差が発生している場合、第1触媒53の酸素吸蔵量OSAallは目標酸素吸蔵量(Cmaxall/2)以外の値に収束することになる。
【0091】
更には、前述のごとく、積分項DFiiの値は前記所定の範囲(−Gi・(Cmaxall/2)≦積分項≦Gi・(Cmaxall/2))を超えることはないから、定常運転状態において積分項DFiiの値が同所定の範囲外の値に収束することに相当する大きな前記誤差が発生したとき、かかる誤差の総てを吸収することができず、筒内燃料供給量偏差DFcがゼロに収束し得なくなる事態(従って、機関10の検出空燃比abyfsが目標空燃比abyfrに収束し得なくなる事態)が発生する。
【0092】
そこで、本装置は、所定の条件が成立する毎に、前述した目標空燃比abyfrが理論空燃比と異なる空燃比となる場合と同様、上記数1の代わりに上記数6に従ってフィードバック制御量DFiを算出することとし、前記誤差の大きさを表す値である内燃機関10が定常運転状態にあるときにおけるフィードバック制御量DFiの値(従って、積分項Gi’・ΣDFcの値が前記所定の範囲を超えるか否かに拘わらず同積分項Gi’・ΣDFcの値)に基づいて前述した基本燃料噴射量補正量取得手段A16により基本燃料噴射量補正量ΔFbaseを求めるための基本燃料噴射量補正量算出テーブルΔFbase(Ga)を求める(更新する)。この基本燃料噴射量補正量算出テーブルΔFbase(Ga)は、更新される毎にその最新のものがバックアップRAM74に記憶される。
【0093】
そして、本装置は、基本燃料噴射量補正量算出テーブルΔFbase(Ga)を更新した後、再び上記数1に従ってフィードバック制御量DFiを算出することを開始するとともに、補正前基本燃料噴射量Fbasebに基本燃料噴射量補正量取得手段A16による最新の基本燃料噴射量補正量算出テーブルΔFbase(Ga)から求めた基本燃料噴射量補正量ΔFbaseを加えることで補正後基本燃料噴射量Fbaseを求める。これにより、前述の大きな誤差が発生している場合であっても前記誤差が基本燃料噴射量の補正により確実に補償されるとともに、第1触媒53の酸素吸蔵量OSAallが確実に目標酸素吸蔵量(Cmaxall/2)に収束するようになる。以上、このように、基本燃料噴射量を補正する手段が基本燃料噴射量補正手段に相当する。
【0094】
(触媒モデル)
次に、第1触媒53の酸素吸蔵量OSAallを求めるための触媒モデルA12について説明する。一般に、触媒にリーンな空燃比のガスが流入したときには、同触媒の上流側においてより多くの酸素が吸蔵され、同触媒にリッチな空燃比のガスが流入したときには、同触媒の上流側から吸蔵されている酸素が消費されていく。従って、触媒内に吸蔵されている酸素は同触媒の排ガスの流れ方向において均一に分布しているわけではない。よって、触媒内の酸素吸蔵量を正確に求めるためには、かかる吸蔵酸素の分布を考慮した計算を行う必要がある。
【0095】
また、触媒の酸素吸蔵量は同触媒内において発生する酸素吸蔵・放出反応の程度に応じて変化する。この酸素吸蔵・放出反応の程度は、触媒に流入する排ガスに含まれる上記酸素吸蔵・放出反応に関係する特定成分の量に依存する。従って、触媒の酸素吸蔵量を正確に求めるためには、前記特定成分の量を考慮に入れた計算を行う必要がある。そこで、本装置は以下に説明する触媒モデルを第1触媒53に適用することにより、同第1触媒53内の酸素吸蔵量OSAallを算出する。
【0096】
この触媒モデルにおいては、図8に模式的に示したように、軸線に直交する断面形状が一定である柱状の触媒を同軸線に直交する面によりN個の(複数の)領域(「ブロック」とも称呼する。)に分割する。即ち、触媒モデルが対象とする触媒は排ガスの流れ方向に沿ってN個のブロックに分割されている。分割された各ブロックの軸線方向の長さはLである。なお、説明の便宜上、各ブロックには、排ガスの流れ方向に沿って上流側から順に図8に示すように番号が付されている。また、任意のi番目のブロックに関連する変数・記号等には、それらの末尾に(i)が付されている。
【0097】
この触媒モデルにおいては、図9に示したように、分割されたブロックのうちのi番目のブロック(i)(特定領域)に注目し、同ブロック(i)における酸素吸蔵・放出反応に関係する特定成分のCPU71の演算周期あたりの収支を考える。このとき、触媒での酸化・還元反応である三元反応は瞬時かつ完全に終了するものと仮定し、その結果としての酸素の過不足に基く酸素の吸蔵・放出反応のみに着目するものとする。この仮定(触媒モデル)は、現実的でありかつ計算精度のよいものである。なお、図9に示した排ガス相は排ガスが通過する空間であり、コート層は触媒機能を発生せしめる白金(Pt)等の貴金属からなる活性成分及び酸素吸蔵機能を発生せしめるセリア(CeO2)等の成分が担持された層である。
【0098】
特定成分は、例えば、酸素(分子)O2、窒素酸化物NOx、一酸化炭素CO、炭化水素HCから選択された成分であってもよいが、この触媒モデルでは、上記三元反応が瞬時かつ完全に終了するものと仮定した状態における排ガスに含まれる酸素(酸素分子及び窒素酸化物の酸素。本明細書では、酸素分子及び窒素酸化物の酸素を総称して「酸素」と称呼する。)(の過不足)を特定成分として選択している。この酸素の量である酸素量CgO2は、同酸素が過剰であるとき(即ち、排ガス中にO2及びNOxが存在する場合)に正の値となり、同酸素が不足しているとき(即ち、排ガス中に未燃HC,COが存在する場合)に負の値となるように計算される。
【0099】
また、注目するブロック(i)において、CPU71の演算周期あたり同ブロック(i)の排ガス相に流入する酸素量CgO2を流入酸素量CginO2(i)、同演算周期あたり同ブロック(i)の排ガス相から流出する酸素量CgO2を流出酸素量CgoutO2(i)と称呼し、同演算周期あたり同ブロック(i)のコート層に吸蔵され又は同コート層から放出される酸素量CgO2を酸素吸蔵量変化量δOSA(i)と称呼する。この酸素吸蔵量変化量δOSA(i)は、酸素がコート層に吸蔵されるときに正の値となり、酸素がコート層から放出されるときに負の値となるように計算される。また、現時点におけるブロック(i)のコート層における酸素吸蔵量を酸素吸蔵量OSA(i)と称呼し、現時点におけるブロック(i)のコート層における最大酸素吸蔵量を最大酸素吸蔵量Cmax(i)と称呼する。
【0100】
いま、図9に示すブロック(i)における酸素量CgO2の上記演算周期あたりの収支を考えると、同ブロック(i)の排ガス相に流入した流入酸素量CginO2(i)のうち酸素吸蔵量変化量δOSA(i)だけがコート層に吸蔵され、同流入酸素量CginO2(i)のうちコート層に吸蔵されなかった残りの酸素量CgO2が流出酸素量CgoutO2(i)となるから、流入酸素量CginO2(i),流出酸素量CgoutO2(i)及び酸素吸蔵量変化量δOSA(i)の間には下記数7に示した関係が成立する。この下記数7に示した関係が本触媒モデルの基本式である。
【0101】
【数7】
CgoutO2(i)=CginO2(i)−δOSA(i)
【0102】
次に、酸素吸蔵量変化量δOSA(i)について考える。流入酸素量CginO2(i)が正の値のときはブロック(i)の排ガス相に流入する排ガス中の酸素が過剰であることを意味し、同排ガス中の酸素の一部はブロック(i)のコート層に吸蔵されるから酸素吸蔵量変化量δOSA(i)は正の値となる。このときの酸素吸蔵反応の量、即ち酸素吸蔵量変化量δOSA(i)は、流入酸素量CginO2(i)の値に比例するとともにブロック(i)の現時点での最大酸素吸蔵量Cmax(i)と現時点での酸素吸蔵量OSA(i)との差の値に比例すると考えられる。従って、流入酸素量CginO2(i)が正の値のとき、酸素吸蔵量変化量δOSA(i)は下記数8,及び下記数9に基づいて算出することができる。
【0103】
【数8】
δOSA(i)=H(i)・CginO2(i)
【0104】
【数9】
H(i)=h・((Cmax(i)−OSA(i))/Cmax(i)) (0 ≦ H(i) < 1)
【0105】
上記数8及び上記数9において、H(i)はブロック(i)における流入酸素量CginO2(i)に対する吸蔵される酸素量(δOSA(i))の割合を示す反応率である。hは反応速度定数であり、本モデルでは正の一定値としているが触媒の温度に応じて変化する正の値(例えば、触媒の温度の増加に応じて単調増加する正の値)としてもよい。また、上記数9における現時点での最大酸素吸蔵量Cmax(i)と現時点での酸素吸蔵量OSA(i)との差の値(Cmax(i)−OSA(i))は、ブロック(i)における現時点での酸素吸蔵余裕量を示している。このように、本触媒モデルでは、少なくとも触媒内の酸素吸蔵量に基いて同触媒が内部に流入する排ガスから吸蔵する酸素量を算出する。
【0106】
一方、流入酸素量CginO2(i)が負の値のときはブロック(i)の排ガス相に流入する排ガスの酸素が不足していることを意味し、同排ガスにはブロック(i)のコート層から放出された酸素が与えられるから酸素吸蔵量変化量δOSA(i)は負の値となる。このときの酸素放出反応の量、即ち酸素吸蔵量変化量δOSA(i)(の絶対値)は、流入酸素量CginO2(i)の値に比例するとともにブロック(i)の現時点での酸素吸蔵量OSA(i)の値に比例すると考えられる。従って、流入酸素量CginO2(i)が負の値のとき、酸素吸蔵量変化量δOSA(i)は上記数8と同一の関係を示す下記数10,及び下記数11に基いて算出することができる。
【0107】
【数10】
δOSA(i)=H(i)・CginO2(i)
【0108】
【数11】
H(i)=h・(OSA(i)/Cmax(i)) (0 ≦ H(i) < 1)
【0109】
上記数10及び上記数11において、H(i)はブロック(i)における流入酸素量CginO2(i)(負の値)に対する放出される酸素量(δOSA(i),負の値)の割合を示す反応率である。hは反応速度定数であり上記数9にて使用されているものと同様である。また、上記数11における現時点での酸素吸蔵量OSA(i)の値は、ブロック(i)における現時点での酸素放出余裕量を示している。このように、本触媒モデルでは、少なくとも触媒内の酸素吸蔵量に基いて同触媒が内部に吸蔵している酸素から放出する酸素量を算出する。
【0110】
なお、上記数9及び上記数11にて使用するブロック(i)における最大酸素吸蔵量Cmax(i)は後述する手法により予め求められている。また、上記数9及び上記数11にて使用するブロック(i)における現時点での酸素吸蔵量OSA(i)は、初期値が付与された時点から現時点までの酸素吸蔵量変化量δOSA(i)の積算値であるから下記数12に基いて算出することができる。
【0111】
【数12】
OSA(i)=ΣδOSA(i) (0 ≦ OSA(i) ≦ Cmax(i))
【0112】
次に、各ブロック間での境界条件について考えると、図8に示したように、互いに隣接する2つのブロックのうちの上流側のブロックの排ガス相の流出面と下流側のブロックの排ガス相の流入面は互いに連続しているから、図9に示したように、ブロック(i)に流入する流入酸素量CginO2(i)は、ブロック(i)に隣接する上流側のブロック(i−1)から流出する流出酸素量CgoutO2(i−1)と等しく、また、ブロック(i)から流出する流出酸素量CgoutO2(i)は、ブロック(i)に隣接する下流側のブロック(i+1)に流入する流入酸素量CginO2(i+1)と等しい。従って、下記数13に示した関係が成立する。換言すると、任意のi番目のブロック(i)の流出酸素量CgoutO2(i)が求まればブロック(i)に隣接する下流側のブロック(i+1)の流入酸素量CginO2(i+1)が求まる。
【0113】
【数13】
CginO2(i+1)=CgoutO2(i)
【0114】
以上のことから、最上流のブロック(1)における流入酸素量CginO2(1)が境界条件として付与されれば、上記数8又は上記数10によりブロック(1)における酸素吸蔵量変化量δOSA(1)が求まり、その結果、上記数12によりブロック(1)における酸素吸蔵量OSA(1)を更新できるとともに上記数7によりブロック(1)における流出酸素量CgoutO2(1)が求まる。ブロック(1)における流出酸素量CgoutO2(1)が求まれば、上記数13によりブロック(2)における流入酸素量CginO2(2)が求まり、その結果、上記数8又は上記数10によりブロック(2)における酸素吸蔵量変化量δOSA(2)が求まる。これにより、上記数12によりブロック(2)における酸素吸蔵量OSA(2)を更新できるとともに、上記数7によりブロック(2)における流出酸素量CgoutO2(2)が求まる。
【0115】
CPU71は、このような処理を所定の演算周期毎に繰り返し実行する。従って、CPU71の演算周期が経過する毎に最上流のブロック(1)における流入酸素量CginO2(1)が境界条件として付与されれば、上記数7〜上記数13より、最上流のブロック(1)から、順次、各ブロック(i) (i=1,2,・・・,N) における酸素吸蔵量OSA(i),流入酸素量CginO2(i),及び流出酸素量CgoutO2(i)を全て算出することができる。これにより、触媒内部の酸素吸蔵量の分布が精度よく計算される。また、各ブロックの酸素吸蔵量OSA(i) (i=1,2,・・・,N) を触媒全体について積算すれば、同触媒全体の酸素吸蔵量OSAallについても精度よく計算することができる。
【0116】
次に、上記数9及び数11において反応率H(i)を求める際に必要となるブロック(i)における最大酸素吸蔵量Cmax(i)の求め方について説明する。図10は、本触媒モデルにおいて、最大酸素吸蔵量Cmax(i)を求める考え方を示した最大酸素吸蔵量分布マップであり、斜線で示された部分の面積は触媒全体の最大酸素吸蔵量Cmaxallの値に対応している。
【0117】
このように、触媒の各ブロック毎の各最大酸素吸蔵量Cmax(n) (n=1,・・・,N)は、同各最大酸素吸蔵量Cmax(n)の総和が同触媒全体の最大酸素吸蔵量Cmaxallの値となるように設定されるとともに、上流側のブロックから下流側のブロックに推移するに従い、所定の勾配をもって線形的に増加するとして取り扱われる。これは、触媒の上流側部分の方が下流側部分に比して、内部に流入する排ガス中に含まれる鉛や硫黄等により被毒し易いので、同上流側部分の最大酸素吸蔵量が同下流側部分に比して低下し易くなるからである。
【0118】
具体的には、本装置は、触媒モデルが対象とする第1触媒53の各ブロック(i)毎の各最大酸素吸蔵量Cmax(i) (i=1,・・・,N)を、図10に示した最大酸素吸蔵量分布マップに基づいた下記数14に基づいて算出する。
【0119】
【数14】
Cmax(i)=A・(i−(N/2))+(Cmaxall/N) (i=1,・・・,N)
【0120】
上記数14において、Aは正の定数であって、各ブロック毎の最大酸素吸蔵量の分布の上記勾配を決定する値である。また、第1触媒53の最大酸素吸蔵量Cmaxallは後述するルーチンにより別途求められている。なお、第1触媒53の各ブロック毎の各最大酸素吸蔵量は、上流側のブロックから下流側のブロックに推移するに従い増加するように設定されていればよく、例えば、非線形的に増加するように設定されていてもよい。このようにして、第1触媒53の各ブロック毎の各最大酸素吸蔵量が算出される。
【0121】
(触媒モデルの適用)
次に、以上説明した触媒モデルを、図11に示したように、第1触媒53に適用し、各種値を求める例について説明する。
【0122】
以下、第1触媒53のi番目のブロックであるブロック(i)における流入酸素量を流入酸素量CginO2(i),流出酸素量を流出酸素量CgoutO2(i),酸素吸蔵量をOSA(i),最大酸素吸蔵量をCmax(i)とそれぞれ称呼する。また、各ブロックの酸素吸蔵量OSA(i) (i=1,2,・・・,N) を積算することにより得られる第1触媒53全体の酸素吸蔵量を酸素吸蔵量OSAallと称呼し、各ブロックの最大酸素吸蔵量Cmax(i) (i=1,2,・・・,N) を積算した値である第1触媒53全体の最大酸素吸蔵量を最大酸素吸蔵量Cmaxallと称呼する。
【0123】
この触媒モデルにおいては、図11に示したように、第1触媒53の各ブロックにおける酸素吸蔵量OSA(i) (i=1,2,・・・,N)の初期値を初期条件として付与するとともに、CPU71の演算周期が経過する毎に、第1触媒53の最上流のブロック(1)における流入酸素量CginO2(1)を境界条件として付与すれば、第1触媒53の各ブロック(i) (i=1,2,・・・,N) における酸素吸蔵量OSA(i),流入酸素量CginO2(i),及び流出酸素量CgoutO2(i)を全て算出することができる。これにより、第1触媒53全体の酸素吸蔵量OSAallも取得・算出することができる。
【0124】
そこで、先ず、第1触媒53の各ブロックにおける酸素吸蔵量の初期値を付与する手法について説明すると、本装置は、下流側空燃比センサ67の出力Voxsが0.7(V)より大きい値を示したとき、即ち、第1触媒53の下流側の空燃比が明白なリッチ空燃比となったときは、同第1触媒53内に酸素が全く存在せず未燃HC,COが浄化されない状態となったことを意味するから、第1触媒53の各ブロックにおける酸素吸蔵量OSA(i) (i=1,2,・・・,N) 及び第1触媒53全体の酸素吸蔵量OSAallを全て「0」に設定する。このようにして、第1触媒53の各ブロックにおける酸素吸蔵量の初期値「0」が初期条件として付与される。
【0125】
次に、第1触媒53の最上流のブロック(1)における流入酸素量CginO2(1)を付与する手法について説明すると、本装置は、下記数15に基づいてCPU71の演算周期毎の流入酸素量CginO2(1)を算出する。
【0126】
【数15】
CginO2(1)=0.23・DFc・stoich・t(NE,Δt)
【0127】
上記数15において、値「0.23」は大気中に含まれる酸素の重量割合であり、DFcは前述した筒内燃料供給量偏差であり、stoichは理論空燃比(例えば、14.7)である。また、t(NE,Δt)は、エンジン回転速度NEと上記演算周期Δtから一吸気行程あたりの時間を一演算周期Δtに変換する関数である。
【0128】
この数15に示したように、筒内燃料供給量偏差DFcに理論空燃比stoichを乗じることで一吸気行程あたりの空気の過剰量が求められ、この一吸気行程あたりの空気の過剰量にt(NE,Δt)を乗じることで、一演算周期Δtにおける空気の過剰量が求められ、この一演算周期Δtにおける空気の過剰量に酸素の重量割合を乗じることで一演算周期Δtにおける酸素の過剰量、即ち現時点での流入酸素量CginO2(1)が求められる。
【0129】
このようにして算出される流入酸素量CginO2(1)は、上記数15から明らかなように、酸素が過剰であるとき(即ち、空燃比がリーンであってDFcが正の値となるとき)に正の値となり、酸素が不足しているとき(即ち、空燃比がリッチであってDFcが負の値となるとき)に負の値となるように計算される。このようにして、CPU71の演算周期毎に第1触媒53の最上流のブロック(1)における流入酸素量CginO2(1)が境界条件として付与される。
【0130】
なお、上記数15に示すように、流入酸素量CginO2(1)(従って、流入酸素量CginO2(i))は、上流側空燃比センサ66の出力値vabyfsと上流側目標値との偏差に基づく値である筒内燃料供給量偏差DFcに基づいて計算されている。また、上記数8及び上記数12に示すように、酸素吸蔵量OSA(1)(従って、酸素吸蔵量OSA(i))は流入酸素量CginO2(1)の時間積分値に基づいて計算されている。従って、第1触媒53の酸素吸蔵量OSAallは、実質的に、上流側空燃比センサ66の出力値vabyfsと上流側目標値との偏差に基づく値の時間積分値に基づいて計算されていると云える。
【0131】
以上のようにして、本装置は、CPU71の演算周期毎に、第1触媒53の各ブロック(i) (i=1,2,・・・,N) における酸素吸蔵量OSA(i),流入酸素量CginO2(i),及び流出酸素量CgoutO2(i)を全て算出するとともに、第1触媒53全体の酸素吸蔵量OSAallの値を算出(推定)する。このように、触媒内における酸素吸蔵反応の速度(前記反応速度定数h)を考慮して構築された触媒モデルA12に基づいて第1触媒53の酸素吸蔵量OSAallを推定する手段が酸素吸蔵量推定手段に相当する。
【0132】
(実際の作動)
次に、上記空燃比制御装置(排気浄化装置)の実際の作動について説明する。(基本燃料噴射量補正量の更新)
CPU71は、図12にフローチャートにより示した基本燃料噴射量補正量の更新を行うルーチンを所定時間の経過毎に繰り返し実行している。従って、所定のタイミングになると、CPU71はステップ1200から処理を開始し、ステップ1205に進んでの値が「0」であって、且つ、基本燃料噴射量補正量ΔFbaseの更新条件が成立しているか否かを判定する。基本燃料噴射量補正量更新処理実行中フラグLEARNは、その値が「1」であるとき基本燃料噴射量補正量更新処理が実行されていることを示し、その値が「0」であるとき基本燃料噴射量補正量更新処理が実行されていないことを示す。なお、基本燃料噴射量補正量更新処理実行中フラグLEARNは、図示しないイグニッションスイッチが「OFF」から「ON」に変更されたときに「0」に設定される。
【0133】
基本燃料噴射量補正量ΔFbaseの更新条件は、内燃機関の暖機運転が終了した場合(具体的には、水温センサ65が検出する冷却水温THWが所定温度未満から同所定温度以上に変化した場合)、燃料の屈折率等を検出して燃料の性状を検出する図示しない周知の燃料性状センサの出力に基づいて燃料の性状に変化があったと判定された場合等に成立する。燃料性状の変化は、例えば、燃料を補給したとき等に発生する。燃料性状の変化があることを更新条件に加えているのは、燃料性状に変化があると、理論空燃比そのものが変化する可能性があるからである。
【0134】
いま、内燃機関10の暖気運転が終了した直後であるものとして説明を続けると、CPU71はステップ1205にて「Yes」と判定してステップ1210に進み、基本燃料噴射量補正量更新処理実行中フラグLEARNの値を「1」に」設定するとともに、続くステップ1215に進んで後述するルーチンにて使用する筒内燃料供給量偏差の積分値SDFc(=ΣDFc)の値を「0」に初期化した後、ステップ1295に進んで本ルーチンを一旦終了する。以降、基本燃料噴射量補正量更新処理実行中フラグLEARNの値は「1」になっているから、CPU71はステップ1205にて「No」と判定してステップ1295に直ちに進んで本ルーチンを一旦終了するようになる。
【0135】
(目標空燃比の設定)
次に、上流側目標空燃比abyfrを設定する際の作動について説明すると、CPU71は図13にフローチャートにより示したルーチンを所定時間の経過毎に繰り返し実行している。従って、所定のタイミングになると、CPU71はステップ1300から処理を開始してステップ1305に進み、現時点での内燃機関の運転状態に応じて上流側目標空燃比abyfr(k)を設定する。
【0136】
次に、CPU71はステップ1310に進み、前回の本ルーチン実行時にて設定された上流側目標空燃比abyfr(k)が理論空燃比(例えば、14.7)であって、今回の本ルーチン実行時にて設定された上流側目標空燃比abyfr(k)が理論空燃比以外の空燃比であるか否か(上流側目標空燃比abyfr(k)が理論空燃比から理論空燃比以外に変化したか否か)を判定する。
【0137】
現時点は、内燃機関10の暖気運転が終了した直後であるから、上流側目標空燃比abyfr(k)は理論空燃比以外から理論空燃比に変化している。従って、CPU71はステップ1310にて「No」と判定してステップ1315に進み、上流側目標空燃比abyfr(k)が理論空燃比以外から理論空燃比に変化したか否かを判定し、同ステップ1315にて「Yes」と判定してステップ1320に進む。
【0138】
CPU71はステップ1320に進むと、空燃比フラグXABYFの値を「0」に設定するとともに、ステップ1395に進んで本ルーチンを一旦終了する。ここで、空燃比フラグXABYFは、その値が「1」のとき上流側目標空燃比abyfr(k)が理論空燃比以外の空燃比に設定されていることを示し、その値が「0」のとき上流側目標空燃比abyfr(k)が理論空燃比に設定されていることを示す。以降、上流側目標空燃比abyfr(k)が理論空燃比に維持されている限りにおいて、CPU71はステップ1310、1315にて共に「0」と判定してステップ1395に進むから空燃比フラグXABYFの値は「0」に維持される。
【0139】
(燃料噴射量の計算)
また、CPU71は、図14にフローチャートにより示した燃料噴射量Fiの計算、及び燃料噴射の指示を行うルーチンを、各気筒のクランク角が各吸気上死点前の所定クランク角度(例えば、BTDC90°CA)となる毎に、繰り返し実行するようになっている。従って、任意の気筒のクランク角度が前記所定クランク角度になると、CPU71はステップ1400から処理を開始してステップ1405に進み、エアフローメータ61により計測された吸入空気流量Ga、エンジン回転速度NE、及び前記上流側目標空燃比abyfr(k)に基づいて、機関の空燃比を同上流側目標空燃比abyfr(k)とするための補正前基本燃料噴射量Fbasebを求める。
【0140】
次いで、CPU71はステップ1410に進み、前記補正前基本燃料噴射量Fbasebに、後述するルーチンにて更新されている最新の基本燃料噴射量補正量算出テーブルΔFbase(Ga)と前記吸入空気流量Gaとから求めた基本燃料噴射量補正量ΔFbaseを加えることで補正後基本燃料噴射量Fbaseを求め、続くステップ1415にて同補正後基本燃料噴射量Fbaseの値を目標筒内燃料供給量Fcr(k)として格納する。
【0141】
次に、CPU71はステップ1420にて、補正後基本燃料噴射量Fbaseに係数Kを乗じた値に後述するフィードバック制御量DFiを加えた値を最終燃料噴射量Fiとして設定する。この係数Kの値は、通常は「1.00」であり、後述するように、最大酸素吸蔵量Cmaxallを求めるために強制的に空燃比を変更するとき、「1.00」以外の所定値に設定される。
【0142】
次に、CPU71はステップ1425に進み、同ステップ1425にて燃料噴射量Fiの燃料を噴射するための指示をインジェクタ39に対して行う。その後、CPU71はステップ1430に進み、その時点の燃料噴射量合計量mfrに最終燃料噴射量Fiを加えた値を新たな燃料噴射量合計量mfrに設定する。この燃料噴射量合計量mfrは、後述する酸素吸蔵量を算出する際に用いられる。その後、CPU71はステップ1495に進み、本ルーチンを一旦終了する。以上により、フィードバック制御量DFiによりフィードバック補正された後の燃料噴射量Fiの燃料が吸気行程を迎える気筒に対して噴射される。
【0143】
(従来の装置を適用した場合のフィードバック制御量の計算)
次に、上記従来の装置を適用した場合のフィードバック制御量DFiを算出する際の作動について説明すると、CPU71は図15にフローチャートにより示したルーチンを所定時間の経過毎に繰り返し実行している。従って、所定のタイミングになると、CPU71はステップ1500から処理を開始し、ステップ1505に進んでフィードバック制御条件が成立しているか否かを判定する。このフィードバック制御条件は、例えば、機関の冷却水温THWが第1所定温度以上であって、機関の一回転当りの吸入空気量(負荷、筒内吸入空気量Mc)が所定値以下であり、上流側空燃比センサ66が正常(活性状態であることを含む。)であり、且つ、後述する最大酸素吸蔵量取得制御実行中フラグXHANの値が「0」のときに成立する。なお、最大酸素吸蔵量取得制御実行中フラグXHANは、その値が「1」のとき最大酸素吸蔵量Cmaxallを取得するために強制的に空燃比を変更する空燃比制御(アクティブ制御)を実行していることを示し、その値が「0」のとき最大酸素吸蔵量Cmaxallを取得するための空燃比制御を実行していないことを示す。
【0144】
いま、フィードバック制御条件が成立しているものとして説明を続けると、CPU71はステップ1505にて「Yes」と判定してステップ1510に進み、基本燃料噴射量補正量更新処理実行中フラグLEARNと空燃比フラグXABYFの少なくとも一つの値が「1」になっているか否かを判定する。
【0145】
現時点では、基本燃料噴射量補正量更新処理実行中フラグLEARNの値が「1」に、空燃比フラグXABYFの値が「0」に設定されているから、CPU71はステップ1510にて「Yes」と判定してステップ1515に進み、現時点の上流側空燃比センサ66の出力値vabyfsを図3に示したテーブルに基いて変換することにより、現時点における検出空燃比abyfsを求める。
【0146】
次に、CPU71はステップ1520に進み、現時点からNストローク(N回の吸気行程)前に吸気行程を迎えた気筒の吸入空気量である筒内吸入空気量Mc(k−N)を前記求めた検出空燃比abyfsで除することにより、現時点からNストローク前の実際の筒内燃料供給量Fc(k−N)を求める。
【0147】
そして、CPU71はステップ1525に進んで、図14のステップ1415にて既に求めてある現時点からNストローク前の目標筒内燃料供給量Fcr(k−N)から前記筒内燃料供給量Fc(k−N)を減じた値を筒内燃料供給量偏差DFcとして設定する。つまり、筒内燃料供給量偏差DFcは、Nストローク前の時点で筒内に供給された燃料の過不足分を表す量となる。次に、CPU71はステップ1530に進み、上記数6に基づくステップ1530内に示した式に従ってフィードバック制御量DFiを求め、続くステップ1535にてその時点における筒内燃料供給量偏差DFcの積分値SDFcに前記ステップ1525にて求めた筒内燃料供給量偏差DFcを加えて、新たな筒内燃料供給量偏差の積分値SDFcを求めた後、続くステップ1595に進んで本ルーチンを一旦終了する。
【0148】
このように、基本燃料噴射量補正量更新処理が実行されている間(基本燃料噴射量補正量更新処理実行中フラグLEARNの値が「1」になっている間)、CPU71がステップ1510にて「Yes」と判定し続けるから、前記従来の装置を適用した場合のフィードバック制御量DFiがステップ1530にて求められ、このフィードバック制御量DFiが前述した図14のステップ1420により燃料噴射量Fiに反映されることで上述した従来の装置を適用した場合の機関の空燃比制御が実行される。
【0149】
一方、ステップ1505の判定時において、フィードバック制御条件が不成立であると、CPU71は同ステップ1505にて「No」と判定してステップ1540に進んでフィードバック制御量DFiの値を「0」に設定し、その後、ステップ1595に進んで本ルーチンを一旦終了する。このように、フィードバック制御条件が不成立であるときは、フィードバック制御量DFiを「0」として前記従来の装置に基づくフィードバック制御による機関の空燃比の補正を行わない。
【0150】
(基本燃料噴射量補正量更新処理終了判定)
また、CPU71は、基本燃料噴射量補正量更新処理終了判定を行うための図16にフローチャートにより示したルーチンを所定時間の経過毎に繰り返し実行している。従って、所定のタイミングになると、CPU71はステップ1600から処理を開始してステップ1605に進んで、基本燃料噴射量補正量更新処理実行中フラグLEARNの値が「1」であり、且つ、内燃機関10が所定時間以上定常運転状態にあるか否かを判定し、同ステップ1605にて「No」と判定する限りにおいてステップ1695に直ちに進んで本ルーチンを一旦終了する。なお、内燃機関10が定常運転状態にあると判定される条件は、冷却水温THWが所定温度以上であり、図示しない車速センサにより得られた車速が所定の高車速以上であり、スロットル弁開度TAの単位時間あたりの変化量が所定量以下である等の場合に成立する。
【0151】
現時点では、基本燃料噴射量補正量更新処理実行中フラグLEARNの値が「1」になっている。従って、いま、内燃機関10が所定時間以上定常運転状態にあるものとすると、CPU71はステップ1605にて「Yes」と判定してステップ1610に進み、現時点にて図15のステップ1530にて算出されているフィードバック制御量DFiの値を更新用フィードバック制御量DFilearnに格納する。なお、現時点では、筒内燃料供給量偏差DFcの値が「0」になっていることから更新用フィードバック制御量DFilearnの値は積分項(Gi’・SDFc)の値と等しくなる。この更新用フィードバック制御量DFilearnの値は、現時点での前述したインジェクタ39等の誤差の大きさを表す。
【0152】
次に、CPU71はステップ1615に進んで前記更新用フィードバック制御量DFilearnの値に基づいて、基本燃料噴射量補正量算出テーブルΔFbase(Ga)を更新し、更新後の最新の基本燃料噴射量補正量算出テーブルΔFbase(Ga)をバックアップRAM74に記憶した後、ステップ1620に進んで基本燃料噴射量補正量更新処理実行中フラグLEARNの値を「0」に設定し、ステップ1695に進んで本ルーチンを一旦終了する。以上の処理により、基本燃料噴射量補正量算出テーブルΔFbase(Ga)が更新される。以降、基本燃料噴射量補正量更新処理実行中フラグLEARNの値が「0」であるから、CPU71はステップ1605にて「No」と判定してステップ1695に直ちに進むようになる。
【0153】
以降、基本燃料噴射量補正量更新処理実行中フラグLEARNの値が「0」であるから、CPU71はステップ1605にて「No」と判定してステップ1695に直ちに進むようになる。また、この結果、現時点では、基本燃料噴射量補正量更新処理実行中フラグLEARNの値も空燃比フラグXABYFの値も「0」になっているから、CPU71は図15のステップ1510にて「No」と判定してステップ1595に直ちに進み、前記従来の装置を適用した場合のフィードバック制御量DFiを算出する処理を実行しないようになる。
【0154】
(本装置を適用した場合のフィードバック制御量の計算)
一方、CPU71は、本装置を適用した場合のフィードバック制御量DFiを算出するための図17にフローチャートにより示したルーチンを所定時間の経過毎に繰り返し実行している。従って、所定のタイミングになると、CPU71はステップ1700から処理を開始し、ステップ1705に進んで前述のステップ1505と同一のフィードバック制御条件が成立しているか否かを判定する。
【0155】
いま、フィードバック制御条件が成立しているものとして説明を続けると、CPU71はステップ1705にて「Yes」と判定してステップ1710に進み、基本燃料噴射量補正量更新処理実行中フラグLEARNの値と空燃比フラグXABYFの値が共に「0」になっているか否かを判定する。
【0156】
先に述べたように、現時点では基本燃料噴射量補正量更新処理実行中フラグLEARNの値も空燃比フラグXABYFの値も「0」に設定されているから、CPU71はステップ1710にて「Yes」と判定してステップ1715〜ステップ1725に順次進み、前述のステップ1515〜ステップ1525と同様にして筒内燃料供給量偏差DFcを求める。
【0157】
次に、CPU71はステップ1730に進み、上記数1に基づくステップ1730内に示した式に従ってフィードバック制御量DFiを求めた後、ステップ1795に進んで本ルーチンを一旦終了する。なお、酸素吸蔵量OSAall、及び最大酸素吸蔵量Cmaxallは、それぞれ、後述するルーチンにより別途求められている。
【0158】
このように、以降、基本燃料噴射量補正量更新処理が実行されておらず(基本燃料噴射量補正量更新処理実行中フラグLEARNの値が「0」になっていて)、且つ、目標空燃比abyfr(k)が理論空燃比に維持されている(空燃比フラグXABYFの値が「0」になっている)限りにおいて、CPU71がステップ1710にて「Yes」と判定し続けるから、本装置を適用した場合のフィードバック制御量DFiがステップ1730にて求められ、このフィードバック制御量DFiが前述した図14のステップ1420により燃料噴射量Fiに反映されることで本装置を適用した場合の機関の空燃比制御が実行される。
【0159】
一方、ステップ1705の判定時において、フィードバック制御条件が不成立であると、CPU71は同ステップ1705にて「No」と判定してステップ1735に進んでフィードバック制御量DFiの値を「0」に設定し、その後、ステップ1595に進んで本ルーチンを一旦終了する。このように、フィードバック制御条件が不成立であるときは、フィードバック制御量DFiを「0」として本装置に基づくフィードバック制御による機関の空燃比の補正を行わない。
【0160】
次に、この状態(基本燃料噴射量補正量更新処理実行中フラグLEARNの値も空燃比フラグXABYFの値も「0」に設定されている状態)から、内燃機関10が急加速運転状態等になって上流側目標空燃比abyfr(k)が理論空燃比以外の空燃比に設定された場合について説明する。この場合、CPU71は図13のステップ1310に進んだとき「Yes」と判定してステップ1325に進み、空燃比フラグXABYFの値を「1」に設定するとともに、続くステップ1330にて図15のステップ1535にて算出される筒内燃料供給量偏差DFcの積分値SDFcの値を「0」に初期化した後、ステップ1395に進んで本ルーチンを一旦終了する。以降、上流側目標空燃比abyfr(k)が理論空燃比以外の空燃比に維持されている限りにおいて、CPU71はステップ1310、1315にて共に「0」と判定してステップ1395に進むから空燃比フラグXABYFの値は「1」に維持される。
【0161】
この結果、現時点では、基本燃料噴射量補正量更新処理実行中フラグLEARNの値は「0」に、空燃比フラグXABYFの値は「1」になっているから、CPU71は図17のステップ1710に進んだとき「No」と判定して本装置を適用した場合のステップ1730にて算出されるフィードバック制御量DFiを算出する処理を実行しないようになる一方で、図15のステップ1510に進んだとき「Yes」と判定して前記従来の装置を適用した場合のステップ1530にて実行されるフィードバック制御量DFiを算出する処理を再び開始するようになる。
【0162】
このようにして、基本燃料噴射量補正量算出テーブルΔFbase(Ga)の更新処理が実行されている間(基本燃料噴射量補正量更新処理実行中フラグLEARNの値が「1」に設定されている間)、及び、上流側目標空燃比abyfr(k)が理論空燃比以外の空燃比に設定されている間(空燃比フラグXABYFの値が「1」に設定されている間)にのみ、本装置を適用した場合のフィードバック制御量DFiの代わりに、前記従来の装置を適用した場合のフィードバック制御量DFiが算出される。
【0163】
(最大酸素吸蔵量の算出)
次に、最大酸素吸蔵量算出のために強制的に空燃比を変更する最大酸素吸蔵量取得制御について説明する。CPU71は図18〜図22のフローチャートにより示された各ルーチンを所定時間の経過毎に実行するようになっている。
【0164】
従って、所定のタイミングになると、CPU71は図18のステップ1800から処理を開始し、ステップ1805に進んで最大酸素吸蔵量取得制御実行中フラグXHANの値が「0」であるか否かを判定する。いま、最大酸素吸蔵量算出のための最大酸素吸蔵量取得制御を行っておらず、且つ、最大酸素吸蔵量取得制御開始条件が成立していないとして説明を続けると、最大酸素吸蔵量取得制御実行中フラグXHANの値は「0」となっている。従って、CPU71はステップ1805にて「Yes」と判定してステップ1810に進み、先に説明した図14のステップ1420にて使用される係数Kの値を1.00に設定する。
【0165】
次いで、CPU71はステップ1815にて最大酸素吸蔵量取得制御の開始判定条件が成立しているか否かを判定する。この最大酸素吸蔵量取得制御の開始判定条件は、冷却水温THWが所定温度以上であり、図示しない車速センサにより得られた車速が所定の高車速以上であって、スロットル弁開度TAの単位時間あたりの変化量が所定量以下であること等の機関が定常運転されている条件が成立し、且つ、下流側空燃比センサ67の出力Voxsが理論空燃比よりもリッチな空燃比に相当する値となっていて、且つ、前回の最大酸素吸蔵量取得制御の終了時点から所定時間以上が経過している場合等に成立する。現段階では、上述したように、最大酸素吸蔵量取得制御の開始判定条件は成立していないから、CPU71はステップ1815にて「No」と判定してステップ1895に進み、本ルーチンを一旦終了する。
【0166】
次に、現時点では最大酸素吸蔵量取得制御を行っていないが、最大酸素吸蔵量取得制御の開始判定条件が成立したものとして説明を続けると、この場合、CPU71はステップ1805にて「Yes」と判定してステップ1810に進み、同ステップ1810にて係数Kの値を1.00に設定する。次いで、CPU71は、前記開始判定条件が成立しているので、ステップ1815にて「Yes」と判定してステップ1820に進み、最大酸素吸蔵量取得制御実行中フラグHANの値を「1」に設定する。
【0167】
そして、CPU71はステップ1825に進み、第1モードに移行するためにModeの値を「1」に設定するとともに、続くステップ1830にて係数Kの値を0.98に設定し、ステップ1895に進んで本ルーチンを一旦終了する。これにより、前述のフィードバック制御条件が成立しなくなるから、CPU71は図15のステップ1505又は図17のステップ1705にて「No」と判定するようになり、フィードバック制御量DFiの値は「0」に設定される。この結果、図14のステップ1420の実行により、補正後基本燃料噴射量Fbaseが0.98倍された値が最終燃料噴射量Fiとして算出され、この最終燃料噴射量Fiの燃料が噴射されるから、機関の空燃比(従って、第1触媒53に流入する排ガスの空燃比)は理論空燃比よりもリーンな所定のリーン空燃比に制御される。
【0168】
以降、CPU71は図18のルーチンの処理をステップ1800から繰り返し実行するが、最大酸素吸蔵量取得制御実行中フラグXHANの値が「1」となっていることから、ステップ1805にて「No」と判定して直ちにステップ1895に進み、本ルーチンを一旦終了するようになる。
【0169】
一方、CPU71は図19に示した第1モード制御ルーチンを所定時間の経過毎に繰り返し実行している。従って、所定のタイミングとなると、CPU71はステップ1900から処理を開始してステップ1905に進み、Modeの値が「1」であるか否かを判定する。このとき、Modeの値が「1」でなければ、CPU71は直ちにステップ1995に進んで本ルーチンを一旦終了する。
【0170】
以下、図18のステップ1825の処理によりModeの値が「1」に変更された直後であるものとして説明を続けると、この場合、CPU71はステップ1905にて「Yes」と判定してステップ1910に進み、下流側空燃比センサ67の出力Voxsが理論空燃比よりもリッチな空燃比を示す値から同理論空燃比よりもリーンな空燃比を示す値に変化したか否かを判定する。
【0171】
現時点では、機関の空燃比を所定のリーン空燃比に変更した直後であるから、下流側空燃比センサ出力Voxsは理論空燃比よりもリッチな空燃比を示している。従って、CPU71はステップ21910にて「No」と判定し、ステップ1995にて本ルーチンを一旦終了する。
【0172】
以降、CPU71は図19のステップ1900〜1910を繰り返し実行する。また、空燃比は所定のリーン空燃比に維持されているから、所定の時間が経過すると第1触媒53の酸素吸蔵量が最大酸素吸蔵量に到達する。従って、これに応じて下流側空燃比センサ出力Voxsはリッチを示す値からリーンを示す値に変化する。これにより、CPU71はステップ1910に進んだとき、「Yes」と判定してステップ1915に進み、Modeの値を「2」に設定するとともに、続くステップ1920にて係数Kの値を1.02に設定し、その後ステップ1995にて本ルーチンを一旦終了する。この結果、機関の空燃比は理論空燃比よりもリッチな上記所定のリッチ空燃比に制御される。
【0173】
また、CPU71は図20に示した第2モード制御ルーチンを所定時間の経過毎に繰り返し実行している。従って、所定のタイミングとなると、CPU71はステップ2000から処理を開始してステップ2005に進み、Modeの値が「2」であるか否かを判定する。いま、先の図19のステップ1915の処理によりModeの値が「2」に変更された直後であるものとして説明を続けると、CPU71はステップ2005にて「Yes」と判定してステップ2010に進み、下流側空燃比センサ67の出力Voxsが理論空燃比よりリーンな空燃比を示す値から同理論空燃比よりリッチな空燃比を示す値に変化したか否かを判定する。この時点では、機関の空燃比を所定のリッチ空燃比に変更した直後であるから第1触媒53から下流側空燃比センサ出力Voxsは理論空燃比よりもリーンな空燃比を示している。従って、CPU71はステップ2010にて「No」と判定し、ステップ2095に進んで本ルーチンを一旦終了する処理を行う。
【0174】
その後、機関の空燃比は引き続いて前記所定のリッチ空燃比に維持されているから、第1触媒53に貯蔵されている酸素が消費されていき、所定の時間が経過すると第1触媒53の酸素吸蔵量が「0」に到達する。この結果、下流側空燃比センサ67の出力Voxsが理論空燃比よりリーンな空燃比を示す値から同理論空燃比よりリッチな空燃比を示す値に変化する。これにより、CPU71はステップ2010からステップ2015に進んでModeの値を「0」に再設定し、続くステップ2020にて最大酸素吸蔵量取得制御実行中フラグXHANの値を「0」に設定した後、ステップ2095に進んで本ルーチンを一旦終了する。
【0175】
これにより、CPU71は図18のルーチンを実行する際、ステップ1805にて「Yes」と判定してステップ1810に進むので、係数Kの値が1.00に戻される。また、前述のフィードバック制御条件が成立していれば、CPU71は図15のステップ1505又は図17のステップ1705にて「Yes」と判定するから、空燃比フィードバック制御が再開される。
【0176】
以上、説明したように、最大酸素吸蔵量所得制御の開始判定条件が成立すると、機関の空燃比(従って、第1触媒53に流入する排ガスの空燃比)が所定のリーン空燃比、及び所定のリッチ空燃比の順に強制的に制御される。
【0177】
次に、最大酸素吸蔵量取得のための酸素吸蔵量の算出における作動について説明する。CPU71は図21のフローチャートにより示されたルーチンを所定時間の経過毎に実行するようになっている。従って、所定のタイミングになると、CPU71は図21のステップ2100から処理を開始し、ステップ2105に進んでModeの値が「2」であるか否かを判定し、Modeの値が「2」であれば同ステップ2105にて「Yes」と判定してステップ2110に進み、下記数16により酸素吸蔵量変化量ΔO2を求める。
【0178】
【数16】
ΔO2=0.23・mfr・(stoich − abyfsave)
【0179】
上記数16において、値「0.23」は大気中に含まれる酸素の重量割合である。mfrは所定時間tsample内の燃料噴射量Fiの合計量であり、stoichは理論空燃比(例えば、14.7)である。abyfsaveは所定時間tsampleにおいて上流側空燃比センサ66により検出された空燃比A/Fの平均値である。この数16に示したように、所定時間tsample内の噴射量の合計量mfrに、検出された空燃比A/Fの平均値の理論空燃比からの偏移(stoich − abyfsave)を乗じることで、同所定時間tsampleにおける空気の消費量(不足量)が求められる。この空気の消費量に酸素の重量割合を乗じることで前記所定時間tsampleにおける酸素の消費量(酸素吸蔵量変化量ΔO2)が求められる。
【0180】
そして、CPU71は、ステップ2115にてその時点の酸素吸蔵量OSA2に上記酸素吸蔵量変化量ΔO2の値を加えた値を新たな酸素吸蔵量OSA2として設定し、その後ステップ2120に進む。
【0181】
このような処置(ステップ2100〜2115)は、Modeの値が「2」である限り繰り返し実行される。この結果、第1触媒53の上流の空燃比が所定のリッチ空燃比とされる第2モード(Mode=2)において、第1触媒53の酸素吸蔵量OSA2が算出されて行く。第2モードのおいては、第1触媒53に貯蔵されている酸素が消費されていくからである。なお、ステップ2105の判定において「No」と判定される場合、CPU71は同ステップ2105からステップ2120に直接進む。そして、CPU71は、ステップ2120に進むと、燃料噴射量Fiの合計量mfrを「0」に設定し、その後ステップ2195に進んで本ルーチンを一旦終了する。
【0182】
次に、最大酸素吸蔵量算出における作動について説明する。CPU71は図22のフローチャートにより示されたルーチンを所定時間の経過毎に実行するようになっている。従って、所定のタイミングになると、CPU71は図22のステップ2200から処理を開始し、ステップ2205に進んで最大酸素吸蔵量取得制御実行中フラグXHANの値が「1」から「0」に変化したか否かをモニタする。このとき、第2モードが終了して、先に説明した図20のステップ2020にて最大酸素吸蔵量取得制御実行中フラグXHANの値が「0」に変更されると、CPU71はステップ2205にて「Yes」と判定してステップ2210に進む。ここで、最大酸素吸蔵量取得制御実行中フラグXHANの値が変化していなければ、CPU71はステップ2205からステップ2295に直接進んで本ルーチンを一旦終了する。
【0183】
いま、第2モードが終了した直後であるとすると、最大酸素吸蔵量取得制御実行中フラグXHANの値が「1」から「0」に変更された直後であるから、CPU71はステップ2205からステップ2210に進み、その時点の酸素吸蔵量OSA2を、第1触媒53の最大酸素吸蔵量Cmaxallとして格納する。
【0184】
次いで、CPU71はステップ2215に進み、カウンタ値nの値を「0」に設定した後、ステップ2220に進んで第1触媒53の各ブロック毎の最大酸素吸蔵量を算出する処理を開始する。まず、CPU71はステップ2220においてカウンタ値nの値を「1」だけ増大して「1」に設定した後、ステップ2225に進んで前記ステップ2210にて取得した最大酸素吸蔵量Cmaxallの値と、カウンタ値nの値と、上記数14に基づくステップ2225内に記載した式とに基づいて第1触媒53のブロック(n)における最大酸素吸蔵量Cmax(n)を算出する。この時点ではカウンタ値nの値は「1」であるから、ブロック(1)における最大酸素吸蔵量Cmax(1)が算出される。
【0185】
そして、CPU71はステップ2230に進んでカウンタ値nの値が第1触媒53のブロック数Nと等しいか否かを判定する。現時点ではカウンタ値nの値は「1」であるから、CPU71はステップ2230にて「No」と判定し、再びステップ2220に戻ってカウンタ値nの値を「1」だけ増大した後ステップ2225及びステップ2230の処理を実行する。即ち、ステップ2225及びステップ2230の処理は、カウンタ値nの値が第1触媒53のブロック数Nと等しくなるまで繰り返し実行される。これにより、第1触媒53の最上流のブロック(1)から最下流のブロック(N)までの各ブロック(n)の最大酸素吸蔵量Cmax(n)の値が順次算出されていく。
【0186】
前述のステップ2220の処理が繰り返されることによりカウンタ値nの値が第1触媒53のブロック数Nと等しくなると、CPU71はステップ2230にて「Yes」と判定してステップ2235に進み、酸素吸蔵量OSA2の値を「0」に設定した後、ステップ2295に進んで本ルーチンを一旦終了する。
【0187】
(各ブロック毎の流出酸素量及び酸素吸蔵量、並びに触媒全体の酸素吸蔵量の計算)
次に、第1触媒53の各ブロック毎の流出酸素量,酸素吸蔵量等の算出における作動について説明する。CPU71は図23及び図24のフローチャートにより示された各ルーチンを所定時間の経過毎に実行するようになっている。
【0188】
従って、所定のタイミングになると、CPU71は、第1触媒53の各ブロック毎の流出酸素量,酸素吸蔵量等を算出するための図23に示したルーチンのステップ2300から処理を開始し、ステップ2305に進んで、図15のステップ1525又は図17のステップ1725にて逐次更新されている筒内燃料供給量偏差DFcの値と、エンジン回転速度NEと、関数t(NE,Δt)と、上記数15に基づくステップ2305内に記載した式とに基づいて、先に説明したように境界条件である第1触媒53のブロック(1)における流入酸素量CginO2(1)を算出する。
【0189】
次いで、CPU71はステップ2310に進んでカウンタ値nの値、及び第1触媒53の酸素吸蔵量OSAallの値をそれぞれ「0」に設定した後、ステップ2315に進んで第1触媒53の各ブロック毎の流出酸素量,酸素吸蔵量等を算出する処理を開始する。まず、CPU71はステップ2315においてカウンタ値nの値を「1」だけ増大して「1」に設定する。カウンタ値nは第1触媒53のブロックの番号を示している。この時点ではカウンタ値nの値は「1」であり、続くステップ2320からステップ2375までの今回の処理においてカウンタ値nの値は「1」に維持されるので、今回の同ステップ2320〜ステップ2375までの処理においては最上流のブロック(1)における計算が実行される。
【0190】
まず、CPU71はステップ2320に進んで、流入酸素量CginO2(1)の値が「0」以上であるか否かを判定し、流入酸素量CginO2(1)の値が「0」以上であれば同ステップ2320において「Yes」と判定するとともにステップ2325に進んで、図22のステップ2225にて既に算出されているブロック(1)の最大酸素吸蔵量Cmax(1)の値と、後述するステップ2360にて前回本ルーチンが実行されたときに計算(更新)されたブロック(1)の酸素吸蔵量OSA(1)の値と、上記数9(の右辺)に基づくステップ2325内に記載した式とに基づいてブロック(1)における反応率Hを算出する。
【0191】
また、ステップ2320における判定において、流入酸素量CginO2(1)の値が「0」以上でなければCPU71は同ステップ2320において「No」と判定するとともにステップ2330に進んで、上記最大酸素吸蔵量Cmax(1)の値と、上記酸素吸蔵量OSA(1)の値と、上記数11(の右辺)に基づくステップ2330内に記載した式とに基づいてブロック(1)における反応率Hを算出する。
【0192】
次いで、CPU71はステップ2335に進み、ステップ2325又はステップ2330にて算出した反応率Hの値と、ステップ2305にて算出したブロック(1)における流入酸素量CginO2(1)の値と、上記数8(の右辺)又は上記数10(の右辺)に基づくステップ2335内に記載した式とに基づいてブロック(1)における酸素吸蔵量変化量δOSA(1)を算出する。
【0193】
次に、CPU71はステップ2340に進んで、後述するステップ2360にて前回本ルーチンが実行されたときに計算されたブロック(1)の酸素吸蔵量OSA(1)の値とステップ2335にて今回算出したブロック(1)の酸素吸蔵量変化量δOSA(1)の値とを積算した値がブロック(1)における最大酸素吸蔵量Cmax(1)の値以下であるか否かを判定する。
【0194】
ここで、前記積算した値が最大酸素吸蔵量Cmax(1)の値以下であればCPU71は同ステップ2340にて「Yes」と判定してステップ2345に進んで同積算した値が「0」以上であるか否かを判定するとともに、同積算した値が「0」以上であれば同ステップ2345にて「Yes」と判定してステップ2360に進んで、同積算した値を新たな酸素吸蔵量OSA(1)として設定する。このように、前記積算した値が「0」以上最大酸素吸蔵量Cmax(1)以下であればステップ2335にて算出された酸素吸蔵量変化量δOSA(1)の値がそのままブロック(1)における酸素吸蔵量変化量として使用される。
【0195】
一方、ステップ2340の判定において、前記積算した値が最大酸素吸蔵量Cmax(1)の値を超えていればCPU71は同ステップ2340にて「No」と判定してステップ2350に進み、最大酸素吸蔵量Cmax(1)の値から前回算出した酸素吸蔵量OSA(1)の値を減算した値を酸素吸蔵量変化量δOSA(1)に格納した後、ステップ2360に進む。このように、前記積算した値がブロック(1)における最大酸素吸蔵量Cmax(1)を超えていれば今回ステップ2360にて算出されるブロック(1)における酸素吸蔵量OSA(1)の値が最大酸素吸蔵量Cmax(1)を超えてしまうことを意味するので、今回ステップ2360にて算出される酸素吸蔵量OSA(1)の値が最大酸素吸蔵量Cmax(1)の値と等しくなるように酸素吸蔵量変化量δOSA(1)が調整される。
【0196】
同様に、ステップ2345の判定において、前記積算した値が「0」未満(負の値)となっていればCPU71は同ステップ2345にて「No」と判定してステップ2355に進み、前回算出した酸素吸蔵量OSA(1)の値に対して符号を反転させた値を酸素吸蔵量変化量δOSA(1)に格納した後、ステップ2360に進む。このように、前記積算した値が「0」未満であれば今回ステップ2360にて算出されるブロック(1)における酸素吸蔵量OSA(1)の値が「0」未満(負の値)になってしまうことを意味するので、今回ステップ2360にて算出される酸素吸蔵量OSA(1)の値が「0」になるように酸素吸蔵量変化量δOSA(1)が調整される。
【0197】
ステップ2360にてブロック(1)における今回の酸素吸蔵量OSA(1)を算出した後、CPU71はステップ2365に進んでステップ2305にて算出したブロック(1)における流入酸素量CginO2(1)の値と、上記調整後のブロック(1)における酸素吸蔵量変化量δOSA(1)の値と、上記数7(の右辺)に基づくステップ2365内に記載した式とに基づいてブロック(1)における流出酸素量CgoutO2(1)を算出する。
【0198】
次いで、CPU71はステップ2370に進んで、(現時点ではステップ2310の実行により「0」となっている)現時点における第1触媒53の酸素吸蔵量OSAallの値にステップ2360にて算出されたブロック(1)における今回の酸素吸蔵量OSA(1)の値を加えた値を新たな酸素吸蔵量OSAallとして格納した後、ステップ2375に進んで、ステップ2365にて算出されたブロック(1)における流出酸素量CgoutO2(1)の値と、上記数13に基づいてブロック(1)に隣接する下流側のブロック(2)における流入酸素量CginO2(2)を算出する。
【0199】
そして、CPU71はステップ2380に進んでカウンタ値nの値が第1触媒53のブロック数Nと等しいか否かを判定する。現時点ではカウンタ値nの値は「1」であるから、CPU71はステップ2380にて「No」と判定し、再びステップ2315に戻ってカウンタ値nの値を「1」だけ増大して「2」に設定した後、続くステップ2320〜ステップ2375までの処理を実行することで次のブロックであるブロック(2)における計算を実行する。このとき、ステップ2365における流入酸素量CginO2(2)の値としては前回ステップ2375にて算出した流入酸素量CginO2(2)の値が使用される。
【0200】
このようにして、ステップ2320〜ステップ2375までの処理は、カウンタ値nの値が第1触媒53のブロック数Nと等しくなるまで繰り返し実行される。これにより、第1触媒53の最上流のブロック(1)から最下流のブロック(N)までの各ブロック(n)の流入酸素量CginO2(n)、流出酸素量CgoutO2(n)、酸素吸蔵量変化量δOSA(n)、及び酸素吸蔵量OSA(n)の値が順次算出されていく。また、ステップ2370の処理が繰り返し実行されることにより、第1触媒53の酸素吸蔵量OSAallも算出される。
【0201】
ステップ2315の処理が繰り返されることによりカウンタ値nの値が第1触媒53のブロック数Nと等しくなると、CPU71はステップ2380にて「Yes」と判定した後、ステップ2395に進んで本ルーチンを一旦終了する。
【0202】
また、CPU71は、所定のタイミングになると、第1触媒53の各ブロック毎の酸素吸蔵量の値、並びに第1触媒53の酸素吸蔵量の値を初期化(クリア)するための図24にフローチャートにより示したルーチンのステップ2400から処理を開始し、ステップ2405に進んで、下流側空燃比センサ67の出力Voxsの値が0.7(V)より大きいか否かをモニタする。このとき、下流側空燃比センサ67の出力Voxsの値が0.7(V)より大きければ、即ち、第1触媒53の下流空燃比がリッチ空燃比であれば、同第1触媒53内全体に吸蔵されている酸素量が「0」であることを意味するので、CPU71はステップ2410に進んで第1触媒53の各ブロック毎の酸素吸蔵量の値及び第1触媒53の酸素吸蔵量の値を総て「0」に設定する処理を開始する。一方、ステップ2405の判定において、下流側空燃比センサ67の出力Voxsの値が0.7(V)以下であれば、CPU71はステップ2405からステップ2495に直接進んで本ルーチンを一旦終了する。
【0203】
いま、ステップ2405の判定において下流側空燃比センサ67の出力Voxsの値が0.7(V)より大きくなっているとすると、CPU71はステップ2410に進み、カウンタ値nの値を「0」に設定した後、ステップ2415に進んでカウンタ値nの値を「1」だけ増大して「1」に設定する。次いで、CPU71はステップ2420に進んで、第1触媒53のブロック(n)における酸素吸蔵量OSA(n)の値を「0」に設定する。この時点ではカウンタ値nの値は「1」であるから、最上流のブロック(1)における酸素吸蔵量OSA(1)の値が「0」に設定される。
【0204】
そして、CPU71はステップ2425に進んでカウンタ値nの値が第1触媒53のブロック数Nと等しいか否かを判定する。現時点ではカウンタ値nの値は「1」であるから、CPU71はステップ2425にて「No」と判定し、再びステップ2415に戻ってカウンタ値nの値を「1」だけ増大した後ステップ2420及びステップ2425の処理を実行する。即ち、ステップ2420及びステップ2425の処理は、カウンタ値nの値が第1触媒53のブロック数Nと等しくなるまで繰り返し実行される。これにより、第1触媒53の最上流のブロック(1)から最下流のブロック(N)までの各ブロック(n)における酸素吸蔵量OSA(n)の値が総て「0」にクリアされる。
【0205】
前述のステップ2415の処理が繰り返されることによりカウンタ値nの値が触媒53のブロック数Nと等しくなると、CPU71はステップ2425にて「Yes」と判定してステップ2430に進み、第1触媒53の酸素吸蔵量OSAallの値を「0」に設定した後、ステップ2495に進んで本ルーチンを一旦終了する。
【0206】
以上、説明したように、本発明による内燃機関の排気浄化装置の実施形態によれば、補正後基本燃料噴射量Fbaseを補正して空燃比フィードバック制御を行うためのフィードバック制御量DFiが、原則的に、DFi=Gp・DFc+Gi・(OSAall−(Cmaxall/2))なる式に基づいて算出される。この結果、比例項であるGp・DFcに基づいて筒内燃料供給量偏差DFcが「0」になるように(従って、上流側空燃比センサ66の検出空燃比abyfsが上流側目標空燃比abyfrになるように)制御されるとともに、第1触媒53の酸素吸蔵量OSAallの値が筒内燃料供給量偏差DFcの時間積分値に基づく値であることから積分項であるGi・(OSAall−(Cmaxall/2))に基づいて、第1触媒53に流入する排ガスの空燃比の理論空燃比からの偏差の収支がゼロになるように制御されるとともに第1触媒53の酸素吸蔵量OSAallが目標酸素吸蔵量(最大酸素吸蔵量Cmaxallの半分の量)になるように制御される。
【0207】
また、前記積分項(Gi・(OSAall−(Cmaxall/2)))の値に制限(−Gi・(Cmaxall/2)≦積分項≦Gi・(Cmaxall/2))があるから、機関の空燃比abyfs(即ち、最終燃料噴射量Fi)は同積分項の作用により不必要に過度に補正されることがなく、空燃比フィードバック制御において大きな外乱が発生しても比較的速やかに収束する。従って、第1触媒53からのエミッションの排出量が少なくなり、また、ドライバビリティが悪化の程度が少なくなった。
【0208】
本発明は上記実施形態に限定されることはなく、本発明の範囲内において種々の変形例を採用することができる。例えば、上記各実施形態においては、筒内燃料供給量偏差DFcに基づいてフィードバック制御量DFiを計算しているが、上流側目標空燃比abyfrと上流側空燃比センサ66による検出空燃比abyfsとの偏差に基づいてフィードバック制御量DFiを計算してもよい。また、上流側空燃比センサ66の出力の目標値である上流側目標値を設定することにより、同上流側目標値と同上流側空燃比センサ66の出力値vabyfsとの偏差に基づいてフィードバック制御量DFiを計算してもよい。
【0209】
ここで、例えば、上流側目標空燃比abyfrと上流側空燃比センサ66による検出空燃比abyfsとの偏差に基づいてフィードバック制御量DFiを計算する場合、触媒モデルが境界条件として使用する第1触媒53の最上流のブロック(1)における流入酸素量CginO2(1)は、上記数15の代わりに、上記数16と類似の下記数17に基づいて計算できる。
【0210】
【数17】
ΔO2=0.23・mfr1・(abyfsave−stoich)
【0211】
上記数17において、mfr1は上記mfrと同様、演算周期Δt内の燃料噴射量Fiの合計量であり、abyfsaveは演算周期Δtにおいて上流側空燃比センサ66により検出された空燃比A/Fの平均値である。この数17に示したように、演算周期Δt内の噴射量の合計量mfr1に、検出された空燃比A/Fの平均値の理論空燃比からの偏移(abyfsave−stoich)を乗じることで、同演算周期Δtにおける空気の過剰量が求められる。この空気の過剰量に酸素の重量割合「0.23」を乗じることで前記演算周期Δtにおける酸素の過剰量、即ち、流入酸素量CginO2(1)が求められる。
【0212】
また、上記各実施形態においては、正負の値を採りえるフィードバック制御量DFiを補正後基本燃料噴射量Fbaseに加算することによりフィードバック制御を実行しているが、同フィードバック制御量DFiに相当するフィードバック制御係数(>0)を設定し、同補正後基本燃料噴射量Fbaseに同フィードバック制御係数を乗算することによりフィードバック制御を実行してもよい。
【0213】
また、上記各実施形態においては、触媒モデルを第1触媒53に適用するにあたり、同第1触媒53をN個のブロックに分割しているが、同第1触媒53を分割することなく触媒モデルを一つのブロックとしての第1触媒53に適用してもよい。
【0214】
また、上記各実施形態においては、フィードバック制御量DFiが、比例項Gp・DFcと積分項Gi・(OSAall−(Cmaxall/2))との和として計算されているが、フィードバック制御量DFiを積分項Gi・(OSAall−(Cmaxall/2))のみで計算してもよい。
【図面の簡単な説明】
【図1】本発明の実施形態に係る空燃比制御装置(排気浄化装置)を適用した内燃機関の概略図である。
【図2】図1に示したエアフローメータの出力電圧と計測された吸入空気流量との関係を示したグラフである。
【図3】図1に示した上流側空燃比センサの出力電圧と空燃比との関係を示したグラフである。
【図4】図1に示した下流側空燃比センサの出力電圧と空燃比との関係を示したグラフである。
【図5】図1に示した空燃比制御装置が空燃比フィードバック制御を実行する際の機能ブロック図である。
【図6】アクセル操作量が急変した場合の機関の空燃比、第1触媒の酸素吸蔵量、及びフィードバック制御量を算出するための積分項の値の変化を、図1に示した空燃比制御装置と従来の装置とで比較しながら示したタイムチャートである。
【図7】従来の空燃比制御装置が空燃比フィードバック制御を実行する際の機能ブロック図である。
【図8】本発明の排気浄化装置が採用する触媒モデルを模式的に示した図である。
【図9】本発明の排気浄化装置が採用する触媒モデルの特定領域に注目したときの同特定領域における酸素吸蔵・放出反応に関係する特定成分の収支を示した図である。
【図10】本発明の排気浄化装置が採用する触媒モデルが対象とする触媒全体の最大酸素吸蔵量から同触媒の各ブロック毎の各最大酸素吸蔵量の分布を求めるためのマップである。
【図11】本発明の排気浄化装置が採用する触媒モデルを第1触媒に適用した場合の模式図である。
【図12】図1に示したCPUが実行する基本燃料噴射量補正量の更新処理の開始判定を行うためのルーチンを示したフローチャートである。
【図13】図1に示したCPUが実行する目標空燃比の設定を行うためのルーチンを示したフローチャートである。
【図14】図1に示したCPUが実行する燃料噴射量の計算等を行うためのルーチンを示したフローチャートである。
【図15】図1に示したCPUが実行する従来の空燃比制御装置を適用した場合のフィードバック制御量を計算するためのルーチンを示したフローチャートである。
【図16】図1に示したCPUが実行する基本燃料噴射量補正量の更新処理の終了判定を行うためのルーチンを示したフローチャートである。
【図17】図1に示したCPUが実行する本発明の排気浄化装置(空燃比制御装置)を適用した場合のフィードバック制御量を計算するためのルーチンを示したフローチャートである。
【図18】図1示したCPUが実行する最大酸素吸蔵量取得制御を開始するか否かを決定するためのルーチンを示したフローチャートである。
【図19】図1に示したCPUが実行する第1モードのルーチンを示したフローチャートである。
【図20】図1に示したCPUが実行する第2モードのルーチンを示したフローチャートである。
【図21】図1に示したCPUが実行する酸素吸蔵量を算出するためのルーチンを示したフローチャートである。
【図22】図1に示したCPUが実行する最大酸素吸蔵量を算出するためのルーチンを示したフローチャートである。
【図23】図1に示したCPUが実行する第1触媒の各ブロック毎の流出酸素量、酸素吸蔵量等を算出するためのルーチンを示したフローチャートである。
【図24】図1に示したCPUが実行する第1触媒の各ブロック毎の酸素吸蔵量をクリアするためのルーチンを示したフローチャートである。
【符号の説明】
10…内燃機関、25…燃焼室、39…インジェクタ、52…エキゾーストパイプ(排気管)、53…三元触媒(第1触媒)、54…三元触媒(第2触媒)、66…上流側空燃比センサ、67…下流側空燃比センサ、70…電気制御装置、71…CPU[0001]
TECHNICAL FIELD OF THE INVENTION
The present invention provides an internal combustion engine including an air-fuel ratio sensor upstream of a catalyst disposed in an exhaust passage of the internal combustion engine, and performing feedback control of an air-fuel ratio of gas flowing into the catalyst based on an output value of the air-fuel ratio sensor. An exhaust purification device.
[0002]
[Prior art]
2. Description of the Related Art Conventionally, a three-way catalyst for purifying exhaust gas of an internal combustion engine (hereinafter, may be simply referred to as a “catalyst”) is disposed in an exhaust passage of the engine. This three-way catalyst has a function of storing oxygen. 2 It has a storage function (hereinafter, this function is referred to as an “oxygen storage function”, and the amount of oxygen stored in the catalyst is referred to as an “oxygen storage amount”). When rich, the stored oxygen oxidizes unburned components such as HC and CO, and when the air-fuel ratio of the inflowing gas is lean, it reduces nitrogen oxides (NOx) to reduce NOx. Oxygen taken from is stored inside. In other words, when the air-fuel ratio of the exhaust gas flowing into the three-way catalyst is rich, the oxygen storage amount of the three-way catalyst decreases, and when the air-fuel ratio of the exhaust gas flowing into the three-way catalyst is lean, The oxygen storage amount of the three-way catalyst increases. Thereby, the three-way catalyst can purify harmful components such as unburned components and nitrogen oxides even when the air-fuel ratio of the engine deviates from the stoichiometric air-fuel ratio.
[0003]
By the way, since the operating state of the internal combustion engine changes every moment, the air-fuel ratio of the engine becomes rich or lean. On the other hand, the oxygen storage amount of the catalyst varies between “0” and the maximum amount of oxygen that can be stored (hereinafter, referred to as “maximum oxygen storage amount”). Therefore, if the air-fuel ratio of the engine becomes rich when the oxygen storage amount in the catalyst is near “0”, the unburned components cannot be sufficiently oxidized in the catalyst, and the unburned components are removed from the catalyst. Will leak out. Conversely, if the air-fuel ratio of the engine becomes lean when the oxygen storage amount in the catalyst is near the maximum oxygen storage amount, nitrogen oxides cannot be sufficiently reduced in the catalyst, and nitrogen oxides are generated from the catalyst. Will leak out. From the above, the air-fuel ratio of the engine (for example, approximately half the maximum oxygen storage amount) is set so that the amount of oxygen stored in the catalyst becomes a predetermined amount (for example, approximately half of the maximum oxygen storage amount) in order to efficiently purify unburned components and nitrogen oxides. Therefore, it is preferable to control the air-fuel ratio of the exhaust gas flowing into the catalyst.
[0004]
For this reason, for example, the exhaust gas purification device (air-fuel ratio control device) for an internal combustion engine disclosed in
[0005]
According to this, the value of the feedback control amount based on the time integral value of the deviation (the integration process (I process)) is determined while the air-fuel ratio of the exhaust gas flowing into the catalyst is rich. By integrating the deviation over time, the air-fuel ratio is updated in a direction in which the air-fuel ratio is corrected to the lean side, and the air-fuel ratio of the exhaust gas flowing into the catalyst becomes the same while the deviation occurs in which the air-fuel ratio becomes lean. By integrating the deviation with time, the air-fuel ratio is updated in a direction in which the air-fuel ratio is corrected to the rich side. Therefore, after the internal combustion engine temporarily transitions from a steady operation state in which the deviation can be zero to a transient operation state in which an excess or deficiency of fuel (that is, an excess or deficiency of air) is likely to occur (therefore, the deviation is likely to occur), The average value of the air-fuel ratio of the exhaust gas flowing into the catalyst during the period before returning to the steady operation state can be guaranteed to be the target air-fuel ratio (therefore, the stoichiometric air-fuel ratio). In other words, the air-fuel ratio of the exhaust gas flowing into the catalyst is controlled by the air-fuel ratio feedback control based on the time integral value of the deviation so that the balance of the deviation becomes zero. Therefore, the oxygen storage amount of the catalyst does not change before and after the above period, and as a result, the oxygen storage amount of the catalyst can be maintained near the predetermined amount.
[0006]
[Patent Document 1]
JP-A-6-280648
[0007]
[Problems to be solved by the invention]
In the disclosed device, for example, when the internal combustion engine is in a transient operation state, such as during a rapid acceleration operation, the air-fuel ratio flowing into the catalyst deviates considerably from the stoichiometric air-fuel ratio for a certain period of time. When the air-fuel ratio becomes lean, the value of the feedback control amount based on the time integrated value of the deviation at the time when the air-fuel ratio returns to the stoichiometric air-fuel ratio is updated to a value that largely corrects the air-fuel ratio to the rich side. ing. Therefore, in this case, the air-fuel ratio of the exhaust gas flowing into the catalyst is returned to the stoichiometric air-fuel ratio from the excessive lean air-fuel ratio, and then is unnecessarily corrected to the excessive rich air-fuel ratio. The oxygen storage amount becomes “0”, the emission of harmful components (emissions) such as unburned HC and CO from the catalyst increases, and the drivability increases due to the fluctuation of the engine output based on the excessive fluctuation of the fuel injection amount. Is worse.
[0008]
Therefore, an object of the present invention is to flow into a catalyst based at least on a time integral of a value based on a deviation from a target value of an output value of an air-fuel ratio sensor upstream of a catalyst disposed in an exhaust passage of an internal combustion engine. An object of the present invention is to provide an exhaust gas purifying apparatus for an internal combustion engine that performs feedback control of an air-fuel ratio of a gas and that does not unnecessarily overcorrect the air-fuel ratio.
[0009]
[Overview of the present invention]
A feature of the present invention is that a catalyst disposed in an exhaust passage of an internal combustion engine, an air-fuel ratio sensor disposed in the exhaust passage upstream of the catalyst, and at least an output value of the air-fuel ratio sensor and a predetermined target Feedback control amount calculating means for calculating a feedback control amount based on a time integration value of a value based on a deviation from the value; and a feedback control amount calculating unit based on the feedback control amount (so that a value based on the deviation becomes zero). In the exhaust gas purifying apparatus for an internal combustion engine, comprising: an air-fuel ratio control unit that feedback-controls an air-fuel ratio of an inflowing gas, when the feedback control amount calculating unit determines that a time integration value of a value based on the deviation exceeds a predetermined range. The feedback control amount is calculated based on a value within the same predetermined range instead of a time integration value of a value based on the deviation.
[0010]
Here, it is preferable that the predetermined target value is set to a value corresponding to a stoichiometric air-fuel ratio. The “value based on the deviation between the output value of the air-fuel ratio sensor and the predetermined target value” is, for example, a deviation between the output value of the air-fuel ratio sensor and the predetermined target value, and a detection value corresponding to the output value of the air-fuel ratio sensor. The difference between the air-fuel ratio (actual air-fuel ratio) and a target air-fuel ratio corresponding to a predetermined target value, the actual in-cylinder value obtained by dividing the in-cylinder intake air amount by the detected air-fuel ratio corresponding to the output value of the air-fuel ratio sensor. It is a deviation from the target in-cylinder fuel supply amount, which is a value obtained by dividing the fuel supply amount and the in-cylinder intake air amount by a target air-fuel ratio corresponding to a predetermined target value, and is not limited thereto.
[0011]
The air-fuel ratio control means for feedback-controlling the air-fuel ratio of the gas flowing into the catalyst based on the feedback control amount (so that the value based on the deviation becomes zero) includes, for example, an air-fuel mixture supplied to the engine. Means for controlling the air-fuel ratio of the air-fuel ratio of the air-fuel mixture supplied to the engine, or in addition to or independently of the control, a nozzle provided in the exhaust passage upstream of the catalyst. A means for controlling the air-fuel ratio of the gas flowing into the catalyst by supplying a predetermined amount of an oxidizing agent (air or the like) or a reducing agent (unburned HC or the like) from the above. If the air-fuel ratio of the air-fuel mixture supplied to the engine is controlled, the air-fuel ratio of the gas flowing into the catalyst can be controlled.
[0012]
According to this, the air-fuel ratio feedback control based on the time integration value of the value based on the deviation causes the air-fuel ratio of the gas flowing into the catalyst to have a balance of zero based on the deviation as described above. Is controlled. Accordingly, the oxygen storage amount of the catalyst can be maintained near a predetermined amount (for example, half of the maximum oxygen storage amount) every time the internal combustion engine enters a steady operation state where the value based on the deviation can be zero. Combustion components and nitrogen oxides can be efficiently purified.
[0013]
Furthermore, according to the exhaust gas purification device, when the time integral of the value based on the deviation exceeds a predetermined range, the time integral based on the value within the predetermined range is used instead of the time integral of the value based on the deviation. A feedback control amount is calculated. In other words, it is configured such that a value exceeding the predetermined range in the time integration value of the value based on the deviation is not reflected as the feedback control amount.
[0014]
Therefore, for example, the air-fuel ratio of the gas flowing into the catalyst becomes one of an excessive lean air-fuel ratio or a rich air-fuel ratio that deviates considerably from the stoichiometric air-fuel ratio for a certain period of time. If the time integration value of the value based on the deviation is a value in a direction in which the air-fuel ratio is corrected to one of the other air-fuel ratios and exceeds the predetermined range, the time integration of the value based on the deviation is performed. Instead of the value, the feedback control amount is determined based on a value within the predetermined range (for example, one of an upper limit value and a lower limit value of the predetermined range). Therefore, by setting the predetermined range to an appropriate range, the air-fuel ratio flowing into the catalyst returns from one of the excessive air-fuel ratios to the stoichiometric air-fuel ratio, and then does not match the excessive one of the other air-fuel ratios. As a result, the amount of emission from the catalyst is reduced, and the degree of deterioration in drivability due to a change in the fuel injection amount can be reduced.
[0015]
In this case, the exhaust gas purification device includes oxygen storage amount estimation means for estimating an oxygen storage amount, which is an oxygen amount stored by the catalyst, based on a time integration value of at least a value based on the deviation, and the feedback control Preferably, the amount calculation means is configured to calculate the feedback control amount based on the estimated oxygen storage amount as a value based on a time integration value of the value based on the deviation.
[0016]
The value based on the deviation between the output value of the air-fuel ratio sensor and a predetermined target value is determined by the excess / deficiency of the fuel supplied into the cylinder (accordingly, the excess / deficiency of air (oxygen) in the gas flowing into the catalyst). The corresponding value. Therefore, the time integral of the value based on the deviation is a value corresponding to the time integral of the excess or deficiency of oxygen in the gas flowing into the catalyst. On the other hand, the amount of change in the amount of oxygen stored in the catalyst is a value corresponding to the amount of oxygen in the gas flowing into the catalyst. Therefore, the amount of oxygen stored in the catalyst is the amount of oxygen in the gas flowing into the catalyst. Becomes a value corresponding to the time integral value of. From the above, the oxygen storage amount of the catalyst is a value that can be calculated based on the time integrated value of the value based on the deviation.
[0017]
Therefore, as described above, even if the air-fuel ratio feedback control is performed based on the oxygen storage amount estimated by the oxygen storage amount estimation means as a value based on the time integration value of the value based on the deviation, As in the case of the air-fuel ratio feedback control based on the time integral value based on the value, the air-fuel ratio flowing into the catalyst can be controlled such that the balance of the value based on the deviation becomes zero. Therefore, each time the internal combustion engine enters a steady operation state in which the deviation can be zero, the oxygen storage amount of the catalyst can be maintained near the predetermined amount (for example, half of the maximum oxygen storage amount).
[0018]
Further, the oxygen storage amount of the catalyst varies between "0" and the maximum oxygen storage amount. Therefore, the oxygen storage amount estimating means estimates the oxygen storage amount of the catalyst so as to be a value that varies between “0” and the maximum oxygen storage amount. Therefore, the feedback control based on the estimated oxygen storage amount is performed. The amount could be calculated based on a value within a predetermined range having a range between “0” and the maximum oxygen storage amount. As a result, an effect equivalent to setting the predetermined range to a predetermined range having a width between “0” and the maximum oxygen storage amount can be obtained, and the air-fuel ratio of the gas flowing into the catalyst can be reduced. By greatly deviating from the target air-fuel ratio corresponding to the predetermined target value, the time integration value of the value based on the deviation becomes a value outside the range of the oxygen storage amount of the catalyst from “0” to the maximum oxygen storage amount. Even in the case of a corresponding value, the air-fuel ratio is not excessively and unnecessarily corrected thereafter, as in the above case.
[0019]
The oxygen storage amount estimating means may estimate the oxygen storage amount of the catalyst such that the oxygen storage amount changes from a predetermined lower limit value equal to or greater than “0” to a predetermined upper limit value equal to or less than the maximum oxygen storage amount. Good. In this case, an effect equivalent to setting the predetermined range to a predetermined range having a width between the predetermined lower limit value and the predetermined upper limit value can be obtained.
[0020]
When the feedback control amount calculating means is configured to calculate the feedback control amount based on the estimated oxygen storage amount, the feedback control amount calculating means determines that the predetermined target oxygen storage amount is the predetermined target oxygen storage amount. It is further preferable that the feedback control amount is calculated based on a deviation from the oxygen storage amount (so that the deviation becomes zero). In this case, the predetermined target oxygen storage amount is preferably substantially half the maximum oxygen storage amount.
[0021]
According to this, the oxygen storage amount of the catalyst can be reliably maintained at the target oxygen storage amount every time the internal combustion engine enters the steady operation state. Further, when the predetermined target oxygen storage amount is set to half the maximum oxygen storage amount, the deviation between the predetermined target oxygen storage amount and the estimated oxygen storage amount is-`` half of the maximum oxygen storage amount. The value varies between “value” and “half the maximum oxygen storage amount”. Therefore, the maximum value of the value in the direction of correcting the air-fuel ratio of the gas flowing into the catalyst to the rich side and the maximum value of the value in the direction of correcting the air-fuel ratio to the lean side in the feedback control amount calculated based on the deviation are calculated. Is equal to As a result, the extent to which the air-fuel ratio of the gas flowing into the catalyst is prevented from being excessively corrected in the rich direction is made equal to the extent to which the air-fuel ratio is prevented from being excessively corrected in the lean direction. Can be.
[0022]
Further, when the feedback control amount calculating means is configured to calculate the feedback control amount based on the estimated oxygen storage amount, the oxygen storage amount estimating means is configured to calculate an oxygen storage reaction in the catalyst. It is preferable to be configured to estimate the oxygen storage amount based on a catalyst model constructed in consideration of speed.
[0023]
Even if the amount of oxygen in the gas flowing into the catalyst is the same, the amount of change in the oxygen storage amount of the catalyst differs depending on the speed of the oxygen storage (release) reaction of the catalyst. Therefore, in order to accurately estimate the oxygen storage amount of the catalyst, it is necessary to consider the speed of the oxygen storage reaction in the catalyst. Therefore, with the above configuration, the oxygen storage amount is estimated based on the catalyst model constructed in consideration of the speed of the oxygen storage reaction in the catalyst. Can be estimated, and as a result, the air-fuel ratio of the gas flowing into the catalyst can be feedback controlled more accurately.
[0024]
The feedback control amount calculating means may be configured to calculate the feedback control amount based on a deviation between the predetermined target oxygen storage amount and the estimated oxygen storage amount, and When the purification device includes target value setting means for setting (changing) the predetermined target value corresponding to a target air-fuel ratio in accordance with an operation state of the internal combustion engine, the feedback control amount calculation means further includes the air-fuel ratio sensor Is configured to calculate the feedback control amount based on a value based on a deviation between the output value and the predetermined target value (so that the deviation becomes zero). When a value corresponding to an air-fuel ratio different from the air-fuel ratio is set, the air-fuel ratio sensor is used instead of the deviation between the predetermined target oxygen storage amount and the estimated oxygen storage amount. It is configured to calculate the feedback control amount based on the time integral value of the value based on the deviation between the output value and the predetermined target value is suitable.
[0025]
That is, in this case, when the predetermined target value is set to a value corresponding to the stoichiometric air-fuel ratio, the feedback control amount calculation means calculates a deviation between the output value of the air-fuel ratio sensor and the predetermined target value. The feedback control amount is calculated based on a value based on the predetermined target oxygen storage amount and the estimated oxygen storage amount, and the predetermined target value corresponds to an air-fuel ratio different from the stoichiometric air-fuel ratio. When the value is set to a value, a value based on a deviation between the output value of the air-fuel ratio sensor and the predetermined target value, and a time integration value of the value based on the deviation (similar to the conventional device disclosed above). Then, a feedback control amount is calculated by performing a proportional / integral process (PI process) on a value based on the deviation.
[0026]
If the target air-fuel ratio corresponding to the predetermined target value is an air-fuel ratio different from the stoichiometric air-fuel ratio, the air-fuel ratio of the gas flowing into the catalyst converges to the target air-fuel ratio when the internal combustion engine is in a steady operation state. However, since an excess or deficiency of oxygen is generated in the gas flowing into the catalyst, the oxygen storage amount of the catalyst changes, and as a result, the oxygen storage amount of the catalyst cannot converge on the target oxygen storage amount. In other words, when the control for converging the air-fuel ratio of the gas flowing into the catalyst to the target air-fuel ratio and the control for converging the oxygen storage amount of the catalyst to the target oxygen storage amount are performed simultaneously, the control targets differ from each other, and both control targets are different. Interfere with each other, it is not possible to perform good air-fuel ratio feedback control. On the other hand, according to the above configuration, when the target air-fuel ratio corresponding to the predetermined target value is different from the stoichiometric air-fuel ratio, the value based on the deviation is proportionally / integrally processed (PI processing). The control target converges the value based on the deviation to zero (therefore, the air-fuel ratio of the gas flowing into the catalyst converges to the target air-fuel ratio). Thus, the above-described control interference can be avoided.
[0027]
Further, the internal combustion engine to which any one of the exhaust purification devices is applied includes a fuel injection unit that injects fuel according to an instruction, and the air-fuel ratio control unit is supplied to the internal combustion engine. A basic fuel injection amount necessary for setting an air-fuel ratio of an air-fuel mixture to a target air-fuel ratio corresponding to the predetermined target value is determined according to an operation state of the internal combustion engine, and the basic fuel injection amount is determined by the feedback control amount. If the configuration is such that the air-fuel ratio of the gas flowing into the catalyst is feedback-controlled by instructing the fuel injection means to inject the fuel of the fuel injection amount obtained by correcting based on In any one of the exhaust emission control devices, the feedback control amount calculation means may further include a time integration value of a value based on a deviation between an output value of the air-fuel ratio sensor and the predetermined target value. Where the feedback control amount is calculated based on the time integral value of the value based on the deviation regardless of whether or not exceeds the value, and the value based on the deviation when the internal combustion engine is in a steady operation state. It is preferable to include a basic fuel injection amount correcting means for correcting the basic fuel injection amount based on a time integration value.
[0028]
In the exhaust gas purifying apparatus, an error (a difference between an indicated fuel injection amount and an actual fuel injection amount) of a fuel injection unit (for example, an injector) and an (instruction) intake air flow rate sensor ( For example, when an error such as an error of the air flow meter (a difference between the measured value of the intake air flow rate and the actual intake air flow rate) occurs, the value based on the deviation becomes zero when the internal combustion engine is in a steady operation state. (Therefore, when the air-fuel ratio of the engine (therefore, the air-fuel ratio of the gas flowing into the catalyst) converges to the target air-fuel ratio corresponding to the predetermined target value), Is a value corresponding to the magnitude of the error. In other words, such an error can be absorbed by the time integral of the value based on the deviation converging to a value corresponding to the magnitude of the error.
[0029]
However, as described above, when the time integral of the value based on the deviation exceeds a predetermined range, the feedback control amount is determined based on the value within the predetermined range instead of the time integral of the value based on the deviation. When configured to calculate, when a large error corresponding to the time integral value of the value based on the deviation becoming a value outside the predetermined range occurs in the steady operation state, all such errors are absorbed. And the value based on the deviation cannot converge to zero.
[0030]
On the other hand, by providing the basic fuel injection amount correcting means as described above, regardless of the magnitude of the error, the basic fuel injection amount correcting means uses the steady operation state used when correcting the basic fuel injection amount. The time integration value of the value based on the deviation is exactly a value corresponding to the magnitude of the error. Therefore, the basic fuel injection amount can be accurately corrected by an amount corresponding to the magnitude of the error based on the time integration value of the value based on the deviation. In other words, even when a large error corresponding to the time integrated value of the value based on the deviation being outside the predetermined range occurs, the error is corrected by the correction of the basic fuel injection amount. Can be surely compensated.
[0031]
BEST MODE FOR CARRYING OUT THE INVENTION
Hereinafter, an embodiment of an air-fuel ratio control device including an exhaust gas purification device for an internal combustion engine according to the present invention will be described with reference to the drawings.
[0032]
FIG. 1 shows a schematic configuration of a system in which an air-fuel ratio control device according to an embodiment of the present invention is applied to a spark ignition type multi-cylinder (four cylinder)
[0033]
The
[0034]
The
[0035]
The
[0036]
The
[0037]
On the other hand, this system includes a hot-wire
[0038]
The hot-wire type
[0039]
The upstream air-
[0040]
The downstream air-fuel ratio sensor 67 is an electromotive force type (concentration cell type) oxygen concentration sensor, and outputs an output value Voxs, which is a voltage that rapidly changes near the stoichiometric air-fuel ratio, as shown in FIG. ing. More specifically, the downstream air-fuel ratio sensor 67 outputs approximately 0.1 (V) when the air-fuel ratio is leaner than the stoichiometric air-fuel ratio and approximately 0.1 V when the air-fuel ratio is richer than the stoichiometric air-fuel ratio. 9 (V), and when the air-fuel ratio is the stoichiometric air-fuel ratio, a voltage of 0.5 (V) is output. The
[0041]
The
[0042]
(Overview of air-fuel ratio feedback control)
Next, an outline of feedback control of the air-fuel ratio of the engine performed by the air-fuel ratio control device including the exhaust gas purification device configured as described above will be described.
[0043]
When the air-fuel ratio of the gas flowing into the
[0044]
Therefore, in order for the
[0045]
In order to control the amount of oxygen stored in the
[0046]
On the other hand, the difference between the output value vabyfs of the upstream air-
[0047]
Therefore, the air-fuel ratio control device (exhaust gas purifying device, hereinafter sometimes simply referred to as “the present device”) of the present embodiment uses the catalyst model described later to calculate the output value vabyfs of the upstream-side air-
[0048]
(Equation 1)
DFi = Gp · DFc + Gi · (OSAall− (Cmaxall / 2))
[0049]
Then, the present apparatus sets the final fuel injection amount Fi to an amount obtained by adding the feedback control amount DFi to a later-described (corrected) basic fuel injection amount Fbase determined according to the operation state of the
[0050]
(Specific contents of air-fuel ratio feedback control)
The air-fuel ratio feedback control described above will be described more specifically. As shown in FIG. 5 which is a functional block diagram, the air-fuel ratio control device is configured to include each unit of A1 to A16. Hereinafter, each means and the like will be described with reference to FIG.
[0051]
<Calculation of basic fuel injection amount (before correction)>
First, the in-cylinder intake air amount calculating means A1 calculates the intake air flow rate Ga measured by the
[0052]
The upstream target air-fuel ratio setting means A2 as the target value setting means is configured to output the upstream target air-fuel ratio corresponding to the upstream target value based on the engine rotational speed NE, the throttle valve opening TA, and the like, which are the operating states of the
[0053]
The pre-correction basic fuel injection amount calculating means A3 converts the in-cylinder intake air amount Mc (k) obtained by the in-cylinder intake air amount calculating means A1 into the upstream target air-fuel ratio set by the upstream target air-fuel ratio setting means A2. By dividing by the abyfr (k), the pre-correction basic fuel injection amount Fbaseb for the current intake stroke for obtaining the air-fuel ratio of the engine at the same upstream-side target air-fuel ratio abyfr (k) is obtained.
[0054]
As described above, the present apparatus obtains the pre-correction basic fuel injection amount Fbaseb using the in-cylinder intake air amount calculation means A1, the upstream target air-fuel ratio setting means A2, and the pre-correction basic fuel injection amount calculation means A3. .
[0055]
<Calculation of final fuel injection amount>
First, the post-correction basic fuel injection amount calculation means A4 adds the latest basic fuel to the pre-correction basic fuel injection amount Fbaseb obtained by the pre-correction basic fuel injection amount calculation means A3 by the basic fuel injection amount correction amount acquisition means A16 described later. The corrected basic fuel injection amount Fbase (that is, the target in-cylinder fuel supply amount Fcr (k)) is obtained by adding the basic fuel injection amount correction amount ΔFbase obtained according to the injection amount correction amount calculation table ΔFbase (Ga). The target in-cylinder fuel supply amount Fcr is stored in the
[0056]
The final fuel injection amount calculating means A5 obtains the final fuel injection amount Fi based on the
[0057]
(Equation 2)
Fi = Fbase + DFi
[0058]
<Calculation of feedback control amount>
First, the table conversion means A6 is a table defining the relationship between the output value vabyfs of the upstream air-
[0059]
The in-cylinder intake air amount delay means A7 includes, among the in-cylinder intake air amount Mc obtained by the in-cylinder intake air amount calculation means A1 for each intake stroke and stored in the
[0060]
The in-cylinder fuel supply amount calculation means A8 calculates the in-cylinder intake air amount Mc (kN) N strokes before the current time obtained by the in-cylinder intake air amount delay means A7 at the current time obtained by the table conversion means A6. By dividing by the detected air-fuel ratio abyfs, the actual in-cylinder fuel supply amount Fc (k−N) N strokes before the current time is obtained. Here, the value N varies depending on the displacement of the
[0061]
As described above, in order to obtain the actual in-cylinder fuel supply amount Fc (k−N) N strokes before the current time, the in-cylinder intake air amount Mc (k−N) N strokes before the current time is detected at the present time. The reason why the fuel ratio is divided by the fuel ratio abyfs is that a time L corresponding to N strokes is required until the fuel-air mixture in the
[0062]
The target in-cylinder fuel supply amount delay means A9 calculates the target in-cylinder fuel supply amount Fcr obtained by the corrected basic fuel injection amount calculation means A4 for each intake stroke and stored in the
[0063]
The in-cylinder fuel supply amount deviation calculating means A10 calculates the target in-cylinder fuel supply amount Fcr (k-N) N strokes before the current time set by the target in-cylinder fuel supply amount delay means A9 based on the following equation (3). By subtracting the actual in-cylinder fuel supply amount Fc (k−N) N strokes before the current time obtained by the in-cylinder fuel supply amount calculation means A8, the above-described in-cylinder fuel supply amount deviation DFc is obtained. The in-cylinder fuel supply amount deviation DFc is an amount representing an excess or deficiency of the fuel supplied into the cylinder at the time before the N stroke, and is an output value vabyfs of the upstream air-
[0064]
[Equation 3]
DFc = Fcr (k−N) −Fc (k−N)
[0065]
The proportional term calculating means A11 performs a proportional process (P process) on the in-cylinder fuel supply amount deviation DFc to obtain a feedback control amount DFip based on the proportional term (P term) based on the following equation (4). In Equation 4 below, Gp is a preset proportional gain (proportional constant (positive value)).
[0066]
(Equation 4)
DFip = Gp · DFc
[0067]
The catalyst model A12 is based on the time integral value of the in-cylinder fuel supply amount deviation DFc, the value of the maximum oxygen storage amount Cmaxall of the
[0068]
The integral term calculating means A13 calculates the integral term based on the following equation 5 based on the value of the oxygen storage amount OSAall of the
[0069]
(Equation 5)
DFii = Gi · (OSAall− (Cmaxall / 2))
[0070]
The feedback control amount calculating means A14 obtains the feedback control amount DFi by adding the feedback control amount DFii based on the integral term to the feedback control amount DFip based on the proportional term. The feedback control amount DFi is used when the final fuel injection amount Fi is obtained by the final fuel injection amount calculation means A5 as described above.
[0071]
In this manner, the present apparatus provides the in-cylinder fuel supply deviation DFc, which is a value based on the deviation between the output value vabyfs of the upstream air-
[0072]
For example, when the air-fuel ratio of the engine suddenly changes and becomes lean, the detected air-fuel ratio abyfs obtained by the table conversion means A6 is leaner than the upstream target air-fuel ratio abyfr set by the upstream target air-fuel ratio setting means A2. It is determined as a value (larger value). For this reason, the actual in-cylinder fuel supply amount Fc (kN) obtained by the in-cylinder fuel supply amount calculation means A8 is equal to the target in-cylinder fuel supply amount Fcr (kr) obtained by the target in-cylinder fuel supply amount delay means A9. kN), and the in-cylinder fuel supply amount deviation DFc is obtained as a large positive value. Further, when the value of the in-cylinder fuel supply amount deviation DFc is a positive value (when the air-fuel ratio of the exhaust gas flowing into the
[0073]
Conversely, when the air-fuel ratio of the engine suddenly changes and becomes rich, the detected air-fuel ratio abyfs is obtained as a value richer (smaller value) than the upstream target air-fuel ratio abyfr. Therefore, the actual in-cylinder fuel supply amount Fc (k-N) becomes a value larger than the target in-cylinder fuel supply amount Fcr (k-N), and the in-cylinder fuel supply amount deviation DFc is obtained as a negative value. Further, when the value of the in-cylinder fuel supply amount deviation DFc is a negative value (when the air-fuel ratio of the exhaust gas flowing into the
[0074]
FIG. 6 shows the oxygen storage of the
[0075]
The conventional device differs from the present device only in the calculation of the feedback control amount DFi. In the conventional device, the feedback control amount DFi is calculated by the PI controller A11 according to the following equation (6). In the following equation 6, Gp is the same proportional gain as in the
[0076]
(Equation 6)
DFi = Gp · DFc + Gi ′ · ΣDFc
[0077]
As shown in FIG. 6, between the time t0 and the time t1, the operation amount Accp of the
[0078]
From this state, assuming that the operation amount Accp of the
[0079]
The oxygen storage amount OSAall of the
[0080]
Therefore, the oxygen storage amount OSAall increases since the air-fuel ratio abyfs of the engine is lean after time t1, reaches the maximum oxygen storage amount Cmaxall at time t2, and is maintained at Cmaxall after time t2. On the other hand, the value of the integral term (Gi′ΣΔDFc) also increases after time t1, and further increases after time t2 because the air-fuel ratio abyfs of the engine is lean. Thereby, the value of the integral term becomes a positive value and becomes a value for correcting the air-fuel ratio abyfs of the engine to the rich side. Therefore, the air-fuel ratio abyfs of the engine starts to decrease from a certain point in time, and at time t3. Change rich. Since the value of the integral term continues to increase until time t3, the value of the integral term becomes a considerably large positive value (maximum value) at the time t3, and the air-fuel ratio abyfs of the engine is largely corrected to the rich side. It has become.
[0081]
Therefore, after time t3, the air-fuel ratio abyfs of the engine is corrected to an unnecessarily large rich side. As a result, after time t3, the oxygen storage amount OSAall, which has started to decrease from the maximum oxygen storage amount Cmaxall because the air-fuel ratio abyfs of the engine becomes considerably rich, reaches “0” at time t5 and becomes “0” after time t5. Is maintained. On the other hand, the value of the integral term also continues to decrease from time t3 to time t7 when the air-fuel ratio abyfs of the engine changes to lean again because the air-fuel ratio abyfs is rich. As a result, the value of the integral term becomes a negative value whose absolute value is considerably large even at the time t7, and is a value that excessively corrects the air-fuel ratio abyfs of the engine to the lean side.
[0082]
Accordingly, the air-fuel ratio abyfs of the engine is corrected to a certain degree to the lean side even after time t7, and converges to the vicinity of the stoichiometric air-fuel ratio after time t9 when it changes to the rich side again. As described above, when the conventional apparatus is applied, the value of the integral term (Gi ′ · ΣDFc) is not limited. Therefore, the air-fuel ratio abyfs (ie, the final fuel injection amount Fi) of the engine is determined by the action of the integral term. Unnecessarily excessive correction results in oscillation, and a large amount of emission (NOx between time t2 and time t3 and HC and CO between time t5 and time t7) is emitted from the
[0083]
On the other hand, when this apparatus is applied, the value of the integral term (Gi · (OSAall− (Cmaxall / 2))) changes in the same manner as the oxygen storage amount OSAall, and the oxygen storage amount OSAall becomes the maximum oxygen storage amount Cmaxall. After the reaching time t2, the value of the integral term is maintained at the maximum value (Cmaxall / 2) and does not become a value larger than the maximum value. Therefore, at time t4 when the air-fuel ratio abyfs of the engine changes from lean to rich, the value of the integral term (Cmaxall / 2) is a value that corrects the air-fuel ratio abyfs of the engine to the rich side, but the value of the integral term is excessively rich. It is not a value that greatly corrects the above.
[0084]
Therefore, the air-fuel ratio abyfs of the engine is not excessively corrected to the lean side after the time t4, and after the time t6 when the engine again changes to the lean side (or the time t8 when the engine again changes to the rich side). Converges near the stoichiometric air-fuel ratio. Further, the oxygen storage amount OSAall does not reach “0” after time t4 when the maximum oxygen storage amount Cmaxall is reached, and converges substantially to the target oxygen storage amount (Cmax / 2) at time t8. Further, the value of the integral term does not fluctuate greatly after time t4, and converges to almost “0” at time t8.
[0085]
As described above, when the present apparatus is applied, the value of the integral term (Gi · (OSAall− (Cmaxall / 2))) is limited to (−Gi · (Cmaxall / 2) ≦ integral term ≦ Gi · (Cmaxall / 2)). ), The air-fuel ratio abyfs of the engine (that is, the final fuel injection amount Fi) is not unnecessarily and excessively corrected by the action of the same integral term (that is, shown by the hatched portion in FIG. 6D). (The amount of correction to the rich side by the integral term is small by an amount corresponding to the area), and converges relatively quickly. Therefore, the emission amount of the emission from the
[0086]
As described above, calculating the feedback control amount DFi based on the value of the integral term (Gi · (OSAall− (Cmaxall / 2))) as in the present device is based on the output value vabyfs of the upstream air-
[0087]
(Treatment when target air-fuel ratio is other than stoichiometric air-fuel ratio)
As shown in the above equation (1), in the present apparatus, the in-cylinder fuel supply amount deviation DFc becomes “0” based on the proportional term DFip (= Gp · DFc) (therefore, the upstream air-fuel ratio sensor 66). So that the detected air-fuel ratio abyfs becomes the upstream target air-fuel ratio abyfr), and the oxygen storage amount OSAall of the
[0088]
However, as described above, the target air-fuel ratio abyfr may have an air-fuel ratio different from the stoichiometric air-fuel ratio. In this case, when the
[0089]
Therefore, when the target air-fuel ratio abyfr has an air-fuel ratio different from the stoichiometric air-fuel ratio, the present device calculates the feedback control amount DFi according to the above-described equation 6 instead of the above-described
[0090]
(Correction of basic fuel injection amount)
In the air-fuel ratio control device described above, for example, the
[0091]
Further, as described above, since the value of the integral term DFii does not exceed the predetermined range (−Gi · (Cmaxall / 2) ≦ integral term ≦ Gi (Cmaxall / 2)), the integral When a large error corresponding to the value of the term DFii converging to a value outside the predetermined range occurs, all such errors cannot be absorbed, and the in-cylinder fuel supply amount deviation DFc becomes zero. A situation in which convergence cannot be achieved (accordingly, a situation in which the detected air-fuel ratio abyfs of the
[0092]
Therefore, the present apparatus sets the feedback control amount DFi in accordance with Equation 6 instead of
[0093]
Then, after updating the basic fuel injection amount correction amount calculation table ΔFbase (Ga), the present apparatus starts calculating the feedback control amount DFi again according to the
[0094]
(Catalyst model)
Next, the catalyst model A12 for obtaining the oxygen storage amount OSAall of the
[0095]
Further, the oxygen storage amount of the catalyst changes according to the degree of the oxygen storage / release reaction generated in the catalyst. The degree of the oxygen storage / release reaction depends on the amount of the specific component related to the oxygen storage / release reaction contained in the exhaust gas flowing into the catalyst. Therefore, in order to accurately determine the oxygen storage amount of the catalyst, it is necessary to perform a calculation in consideration of the amount of the specific component. Therefore, the present apparatus calculates the oxygen storage amount OSAall in the
[0096]
In this catalyst model, as schematically shown in FIG. 8, a columnar catalyst having a constant cross-sectional shape orthogonal to the axis is divided into N (plural) regions ("blocks") by a surface orthogonal to the coaxial line. ). That is, the catalyst targeted by the catalyst model is divided into N blocks along the flow direction of the exhaust gas. The length in the axial direction of each divided block is L. For convenience of explanation, each block is numbered sequentially from the upstream side along the flow direction of the exhaust gas as shown in FIG. In addition, variables, symbols, and the like related to an arbitrary i-th block have (i) at the end thereof.
[0097]
In this catalyst model, as shown in FIG. 9, the i-th block (i) (specific region) among the divided blocks is focused on and related to the oxygen storage / release reaction in the block (i). Consider the balance of the specific component per calculation cycle of the
[0098]
The specific component is, for example, oxygen (molecule) O 2 May be a component selected from nitrogen oxides NOx, carbon monoxide CO, and hydrocarbon HC. However, in this catalyst model, the exhaust gas in a state where the three-way reaction is assumed to end instantaneously and completely is described. Oxygen (oxygen of oxygen molecules and nitrogen oxides. In this specification, oxygen of oxygen molecules and nitrogen oxides are collectively referred to as “oxygen”) (excess or insufficient) is selected as a specific component. I have. The oxygen amount CgO2, which is the amount of oxygen, is determined when the oxygen is excessive (that is, Og in the exhaust gas). 2 And when NOx is present), and becomes a negative value when the oxygen is insufficient (ie, when unburned HC and CO are present in the exhaust gas).
[0099]
In the block (i) of interest, the amount of oxygen CgO2 flowing into the exhaust gas phase of the block (i) per calculation cycle of the
[0100]
Considering the balance of the oxygen amount CgO2 in the block (i) shown in FIG. 9 per the above calculation cycle, the oxygen storage amount change amount in the inflowing oxygen amount CginO2 (i) flowing into the exhaust gas phase of the block (i) is considered. Only δOSA (i) is occluded in the coat layer, and the remaining oxygen amount CgO2 not occluded in the coat layer of the inflow oxygen amount CginO2 (i) becomes the outflow oxygen amount CgoutO2 (i). (I) The relationship shown in the following equation 7 is established between the outflow oxygen amount CgoutO2 (i) and the oxygen storage amount change amount δOSA (i). The relationship shown in Equation 7 below is the basic equation of the present catalyst model.
[0101]
(Equation 7)
CgoutO2 (i) = CginO2 (i) -δOSA (i)
[0102]
Next, the oxygen storage amount change amount δOSA (i) will be considered. When the inflowing oxygen amount CginO2 (i) is a positive value, it means that the oxygen in the exhaust gas flowing into the exhaust gas phase of the block (i) is excessive, and a part of the oxygen in the exhaust gas is in the block (i). Therefore, the oxygen storage amount change amount δOSA (i) becomes a positive value because it is occluded by the coat layer. The amount of the oxygen storage reaction at this time, that is, the oxygen storage amount change amount δOSA (i) is proportional to the value of the inflowing oxygen amount CginO2 (i), and the current maximum oxygen storage amount Cmax (i) of the block (i). And the oxygen storage amount OSA (i) at this time. Therefore, when the inflowing oxygen amount CginO2 (i) is a positive value, the oxygen storage amount change amount δOSA (i) can be calculated based on the following
[0103]
(Equation 8)
δOSA (i) = H (i) · CginO2 (i)
[0104]
(Equation 9)
H (i) = h · ((Cmax (i) −OSA (i)) / Cmax (i)) (0 ≦ H (i) <1)
[0105]
In the
[0106]
On the other hand, when the inflowing oxygen amount CginO2 (i) is a negative value, it means that the exhaust gas flowing into the exhaust gas phase of the block (i) is short of oxygen, and the exhaust gas in the exhaust gas phase of the block (i) has a coating layer of the block (i). , The amount of change in the oxygen storage amount δOSA (i) becomes a negative value. The amount of oxygen release reaction at this time, that is, the oxygen storage amount change amount δOSA (i) (absolute value) is proportional to the value of the inflowing oxygen amount CginO2 (i) and the oxygen storage amount at the present time in the block (i). It is considered to be proportional to the value of OSA (i). Therefore, when the inflowing oxygen amount CginO2 (i) is a negative value, the oxygen storage amount change amount δOSA (i) can be calculated based on the following
[0107]
(Equation 10)
δOSA (i) = H (i) · CginO2 (i)
[0108]
[Equation 11]
H (i) = h · (OSA (i) / Cmax (i)) (0 ≦ H (i) <1)
[0109]
In the above equations (10) and (11), H (i) represents the ratio of the released oxygen amount (δOSA (i), negative value) to the inflowing oxygen amount CginO2 (i) (negative value) in block (i). It is a reaction rate shown. h is a reaction rate constant, which is the same as that used in Equation 9 above. Further, the value of the current oxygen storage amount OSA (i) in the above equation 11 indicates the current oxygen release allowance in the block (i). As described above, in the present catalyst model, the amount of oxygen released from the oxygen stored inside the catalyst is calculated based on at least the amount of oxygen stored in the catalyst.
[0110]
Note that the maximum oxygen storage amount Cmax (i) in the block (i) used in Expressions 9 and 11 is obtained in advance by a method described later. Further, the oxygen storage amount OSA (i) at the present time in the block (i) used in Expressions 9 and 11 is the oxygen storage amount change amount δOSA (i) from the time when the initial value is given to the present time. Can be calculated based on the following equation (12).
[0111]
(Equation 12)
OSA (i) = ΣδOSA (i) (0 ≦ OSA (i) ≦ Cmax (i))
[0112]
Next, considering the boundary conditions between the blocks, as shown in FIG. 8, the outflow surface of the exhaust gas phase of the upstream block and the exhaust gas phase of the downstream block of the two blocks adjacent to each other, as shown in FIG. Since the inflow surfaces are continuous with each other, as shown in FIG. 9, the inflowing oxygen amount CginO2 (i) flowing into the block (i) is equal to the upstream block (i-1) adjacent to the block (i). The amount of oxygen CgoutO2 (i) flowing out of the block (i) is equal to the amount of oxygen CgoutO2 (i-1) flowing out of the block (i), and flows into the downstream block (i + 1) adjacent to the block (i). It is equal to the inflowing oxygen amount CginO2 (i + 1). Therefore, the relationship shown in the following Expression 13 is established. In other words, if the outflow oxygen amount CgoutO2 (i) of an arbitrary i-th block (i) is obtained, the inflow oxygen amount CginO2 (i + 1) of the downstream block (i + 1) adjacent to the block (i) is obtained.
[0113]
(Equation 13)
CginO2 (i + 1) = CgoutO2 (i)
[0114]
From the above, if the inflowing oxygen amount CginO2 (1) in the most upstream block (1) is given as a boundary condition, the oxygen storage amount change amount δOSA (1) in the block (1) is calculated by the above equation (8) or (10). ) Is obtained. As a result, the oxygen storage amount OSA (1) in the block (1) can be updated by the above equation (12), and the outflow oxygen amount CgoutO2 (1) in the block (1) can be obtained by the above equation (7). If the outflow oxygen amount CgoutO2 (1) in the block (1) is obtained, the inflow oxygen amount CginO2 (2) in the block (2) is obtained from the above equation (13), and as a result, the block (2) is obtained from the above equation (8) or (10). ) Is obtained. As a result, the oxygen storage amount OSA (2) in the block (2) can be updated by the equation (12), and the outflow oxygen amount CgoutO2 (2) in the block (2) can be obtained by the equation (7).
[0115]
The
[0116]
Next, a method of obtaining the maximum oxygen storage amount Cmax (i) in the block (i) necessary for obtaining the reaction rate H (i) in the above equations 9 and 11 will be described. FIG. 10 is a maximum oxygen storage amount distribution map showing the concept of obtaining the maximum oxygen storage amount Cmax (i) in the present catalyst model. The area of the hatched portion indicates the maximum oxygen storage amount Cmaxall of the entire catalyst. Corresponds to the value.
[0117]
As described above, the maximum oxygen storage amount Cmax (n) (n = 1,..., N) for each block of the catalyst is determined by the sum of the maximum oxygen storage amounts Cmax (n) being the maximum of the entire catalyst. The oxygen storage amount Cmaxall is set to be equal to the value, and is treated as increasing linearly with a predetermined gradient as the block moves from the upstream block to the downstream block. This is because the upstream portion of the catalyst is more susceptible to poisoning by lead or sulfur contained in the exhaust gas flowing into the catalyst than the downstream portion, and therefore the maximum oxygen storage amount of the upstream portion is the same. This is because it tends to decrease as compared with the downstream portion.
[0118]
Specifically, the present apparatus calculates the maximum oxygen storage amount Cmax (i) (i = 1,..., N) for each block (i) of the
[0119]
[Equation 14]
Cmax (i) = A · (i− (N / 2)) + (Cmaxall / N) (i = 1,..., N)
[0120]
In the above equation (14), A is a positive constant and is a value for determining the gradient of the distribution of the maximum oxygen storage amount for each block. Further, the maximum oxygen storage amount Cmaxall of the
[0121]
(Application of catalyst model)
Next, an example in which the above-described catalyst model is applied to the
[0122]
Hereinafter, in the block (i) that is the i-th block of the
[0123]
In this catalyst model, as shown in FIG. 11, an initial value of the oxygen storage amount OSA (i) (i = 1, 2,..., N) in each block of the
[0124]
Therefore, first, a method of giving an initial value of the oxygen storage amount in each block of the
[0125]
Next, a method of giving the inflowing oxygen amount CginO2 (1) in the most upstream block (1) of the
[0126]
(Equation 15)
CginO2 (1) = 0.23 · DFc · stoich · t (NE, Δt)
[0127]
In the above Expression 15, the value “0.23” is the weight ratio of oxygen contained in the atmosphere, DFc is the above-described cylinder fuel supply amount deviation, and stoich is the stoichiometric air-fuel ratio (for example, 14.7). is there. Further, t (NE, Δt) is a function for converting the time per intake stroke from the engine rotation speed NE and the above-mentioned calculation cycle Δt to one calculation cycle Δt.
[0128]
As shown in Equation 15, the excess amount of air per intake stroke is obtained by multiplying the in-cylinder fuel supply amount deviation DFc by the stoichiometric air-fuel ratio stoich, and the excess amount of air per intake stroke is t By multiplying by (NE, Δt), the excess amount of air in one operation period Δt is obtained. By multiplying the excess amount of air in one operation period Δt by the weight ratio of oxygen, the excess amount of oxygen in one operation period Δt is obtained. The amount, that is, the current inflowing oxygen amount CginO2 (1) is obtained.
[0129]
As is apparent from Equation 15, the inflowing oxygen amount CginO2 (1) calculated as described above is obtained when oxygen is excessive (that is, when the air-fuel ratio is lean and DFc is a positive value). Is calculated to be a positive value, and to be a negative value when oxygen is insufficient (that is, when the air-fuel ratio is rich and DFc is a negative value). In this way, the inflowing oxygen amount CginO2 (1) in the most upstream block (1) of the
[0130]
As shown in the above equation 15, the inflowing oxygen amount CginO2 (1) (accordingly, the inflowing oxygen amount CginO2 (i)) is based on the deviation between the output value vabyfs of the upstream air-
[0131]
As described above, in this apparatus, the oxygen storage amount OSA (i) in each block (i) (i = 1, 2,..., N) of the
[0132]
(Actual operation)
Next, the actual operation of the air-fuel ratio control device (exhaust gas purification device) will be described. (Update of basic fuel injection amount correction amount)
The
[0133]
The basic fuel injection amount correction amount ΔFbase is updated when the warm-up operation of the internal combustion engine is completed (specifically, when the cooling water temperature THW detected by the
[0134]
Now, assuming that it is immediately after the warm-up operation of the
[0135]
(Setting of target air-fuel ratio)
Next, the operation when setting the upstream target air-fuel ratio abyfr will be described. The
[0136]
Next, the
[0137]
At this point, since the warm-up operation of the
[0138]
When the
[0139]
(Calculation of fuel injection amount)
In addition, the
[0140]
Next, the
[0141]
Next, in
[0142]
Next, the
[0143]
(Calculation of feedback control amount when conventional device is applied)
Next, the operation of calculating the feedback control amount DFi when the above-described conventional device is applied will be described. The
[0144]
Now, assuming that the feedback control condition is satisfied, the
[0145]
At this time, since the value of the basic fuel injection amount correction amount update processing execution flag LEARN is set to “1” and the value of the air-fuel ratio flag XABYF is set to “0”, the
[0146]
Next, the
[0147]
Then, the
[0148]
As described above, while the basic fuel injection amount correction amount updating process is being executed (while the value of the basic fuel injection amount correction amount updating process execution flag LEARN is “1”), the
[0149]
On the other hand, if the feedback control condition is not satisfied at the time of determination in
[0150]
(Basic fuel injection amount correction amount update process end determination)
In addition, the
[0151]
At this time, the value of the flag LEARN during execution of the basic fuel injection amount correction amount update process is “1”. Accordingly, assuming that the
[0152]
Next, the
[0153]
Thereafter, since the value of the basic fuel injection amount correction amount updating process execution flag LEARN is “0”, the
[0154]
(Calculation of feedback control amount when this device is applied)
On the other hand, the
[0155]
Now, assuming that the feedback control condition is satisfied, the
[0156]
As described above, since both the value of the basic fuel injection amount correction amount updating process executing flag LEARN and the value of the air-fuel ratio flag XABYF are currently set to “0”, the
[0157]
Next, the
[0158]
As described above, thereafter, the basic fuel injection amount correction amount update processing has not been executed (the value of the basic fuel injection amount correction amount update processing execution flag LEARN has been set to “0”), and the target air-fuel ratio has been changed. As long as abyfr (k) is maintained at the stoichiometric air-fuel ratio (the value of the air-fuel ratio flag XABYF is “0”), the
[0159]
On the other hand, if the feedback control condition is not satisfied at the time of the determination in
[0160]
Next, from this state (a state in which the value of the basic fuel injection amount correction amount updating process execution flag LEARN and the value of the air-fuel ratio flag XABYF are set to “0”), the
[0161]
As a result, at this time, the value of the basic fuel injection amount correction amount updating process execution flag LEARN is “0” and the value of the air-fuel ratio flag XABYF is “1”, so that the
[0162]
In this way, while the updating process of the basic fuel injection amount correction amount calculation table ΔFbase (Ga) is being executed (the value of the basic fuel injection amount correction amount updating process execution flag LEARN is set to “1”). Only when the upstream target air-fuel ratio abyfr (k) is set to an air-fuel ratio other than the stoichiometric air-fuel ratio (while the value of the air-fuel ratio flag XABYF is set to “1”). Instead of the feedback control amount DFi when the device is applied, the feedback control amount DFi when the conventional device is applied is calculated.
[0163]
(Calculation of maximum oxygen storage amount)
Next, the maximum oxygen storage amount acquisition control for forcibly changing the air-fuel ratio for calculating the maximum oxygen storage amount will be described. The
[0164]
Therefore, at a predetermined timing, the
[0165]
Next, at
[0166]
Next, although the maximum oxygen storage amount acquisition control is not performed at the present time, the description is continued assuming that the start determination condition of the maximum oxygen storage amount acquisition control is satisfied. In this case, the
[0167]
Then, the
[0168]
Thereafter, the
[0169]
On the other hand, the
[0170]
Hereinafter, the description will be continued assuming that it is immediately after the value of Mode is changed to “1” by the processing of
[0171]
At present, since the air-fuel ratio of the engine has just been changed to the predetermined lean air-fuel ratio, the downstream-side air-fuel ratio sensor output Voxs indicates an air-fuel ratio richer than the stoichiometric air-fuel ratio. Therefore, the
[0172]
Thereafter, the
[0173]
The
[0174]
Thereafter, since the air-fuel ratio of the engine is continuously maintained at the predetermined rich air-fuel ratio, the oxygen stored in the
[0175]
Accordingly, when executing the routine in FIG. 18, the
[0176]
As described above, when the start determination condition of the maximum oxygen storage amount income control is satisfied, the air-fuel ratio of the engine (therefore, the air-fuel ratio of the exhaust gas flowing into the first catalyst 53) becomes the predetermined lean air-fuel ratio and the predetermined air-fuel ratio. Forced control is performed in the order of the rich air-fuel ratio.
[0177]
Next, the operation in calculating the oxygen storage amount for obtaining the maximum oxygen storage amount will be described. The
[0178]
(Equation 16)
ΔO2 = 0.23 · mfr · (stoich-abyfsave)
[0179]
In Equation 16, the value “0.23” is the weight ratio of oxygen contained in the atmosphere. mfr is the total amount of the fuel injection amount Fi within the predetermined time tsample, and stoich is the stoichiometric air-fuel ratio (for example, 14.7). abyfsave is an average value of the air-fuel ratio A / F detected by the upstream air-
[0180]
Then, in
[0181]
Such a process (steps 2100 to 2115) is repeatedly executed as long as the value of Mode is “2”. As a result, in the second mode (Mode = 2) in which the air-fuel ratio upstream of the
[0182]
Next, the operation in calculating the maximum oxygen storage amount will be described. The
[0183]
Now, assuming that immediately after the second mode is ended, immediately after the value of the maximum oxygen storage amount acquisition control execution flag XHAN has been changed from “1” to “0”, the
[0184]
Next, the
[0185]
Then, the
[0186]
When the value of the counter value n becomes equal to the number N of blocks of the
[0187]
(Calculation of outflow oxygen amount and oxygen storage amount for each block, and oxygen storage amount of the entire catalyst)
Next, the operation of the
[0188]
Therefore, at a predetermined timing, the
[0189]
Next, the
[0190]
First, the
[0191]
If the value of the inflowing oxygen amount CginO2 (1) is not equal to or greater than “0” in the determination in
[0192]
Next, the
[0193]
Next, the
[0194]
Here, if the integrated value is equal to or less than the value of the maximum oxygen storage amount Cmax (1), the
[0195]
On the other hand, if it is determined in
[0196]
Similarly, in the determination in
[0197]
After calculating the current oxygen storage amount OSA (1) in block (1) in
[0198]
Next, the
[0199]
Then, the
[0200]
In this way, the processing from
[0201]
When the value of the counter value n becomes equal to the number N of blocks of the
[0202]
At a predetermined timing, the
[0203]
If it is determined in
[0204]
Then, the
[0205]
When the value of the counter value n becomes equal to the number of blocks N of the
[0206]
As described above, according to the embodiment of the exhaust gas purifying apparatus for an internal combustion engine according to the present invention, the feedback control amount DFi for performing the air-fuel ratio feedback control by correcting the corrected basic fuel injection amount Fbase is basically the same. DFi = Gp な る DFc + Gi ・ (OSAall- (Cmaxall / 2)). As a result, the in-cylinder fuel supply amount deviation DFc becomes “0” based on the proportional term Gp · DFc (therefore, the air-fuel ratio abyfs detected by the upstream air-
[0207]
In addition, since there is a limit (−Gi · (Cmaxall / 2) ≦ integration term ≦ Gi · (Cmaxall / 2)) in the value of the integral term (Gi · (OSAall− (Cmaxall / 2))), the engine is empty. The fuel ratio abyfs (that is, the final fuel injection amount Fi) is not unnecessarily and excessively corrected by the action of the same integral term, and converges relatively quickly even if a large disturbance occurs in the air-fuel ratio feedback control. Therefore, the emission amount of the emission from the
[0208]
The present invention is not limited to the above embodiments, and various modifications can be adopted within the scope of the present invention. For example, in each of the above embodiments, the feedback control amount DFi is calculated based on the in-cylinder fuel supply amount deviation DFc, but the upstream target air-fuel ratio abyfr and the air-fuel ratio abyfs detected by the upstream air-
[0209]
Here, for example, when calculating the feedback control amount DFi based on the deviation between the upstream target air-fuel ratio abyfr and the air-fuel ratio abyfs detected by the upstream air-
[0210]
[Equation 17]
ΔO2 = 0.23 · mfr1 · (abyfsave-stoich)
[0211]
In the above formula 17, mfr1 is the total amount of the fuel injection amount Fi within the calculation cycle Δt, similarly to the above mfr, and abyfsave is the average of the air-fuel ratio A / F detected by the upstream air-
[0212]
Further, in each of the above embodiments, the feedback control is performed by adding the feedback control amount DFi that can take a positive or negative value to the corrected basic fuel injection amount Fbase, but the feedback control amount corresponding to the feedback control amount DFi is added. The feedback control may be performed by setting a control coefficient (> 0) and multiplying the corrected basic fuel injection amount Fbase by the feedback control coefficient.
[0213]
Further, in each of the above embodiments, when the catalyst model is applied to the
[0214]
In each of the above embodiments, the feedback control amount DFi is calculated as the sum of the proportional term Gp · DFc and the integral term Gi · (OSAall− (Cmaxall / 2)). The calculation may be performed only with the term Gi · (OSAall− (Cmaxall / 2)).
[Brief description of the drawings]
FIG. 1 is a schematic diagram of an internal combustion engine to which an air-fuel ratio control device (exhaust gas purification device) according to an embodiment of the present invention is applied.
FIG. 2 is a graph showing a relationship between an output voltage of the air flow meter shown in FIG. 1 and a measured intake air flow rate.
FIG. 3 is a graph showing a relationship between an output voltage of an upstream air-fuel ratio sensor shown in FIG. 1 and an air-fuel ratio.
FIG. 4 is a graph showing a relationship between an output voltage of a downstream air-fuel ratio sensor shown in FIG. 1 and an air-fuel ratio.
FIG. 5 is a functional block diagram when the air-fuel ratio control device shown in FIG. 1 executes the air-fuel ratio feedback control.
FIG. 6 shows changes in the value of the integral term for calculating the air-fuel ratio of the engine, the oxygen storage amount of the first catalyst, and the feedback control amount when the accelerator operation amount changes suddenly. 6 is a time chart showing a comparison between a device and a conventional device.
FIG. 7 is a functional block diagram when a conventional air-fuel ratio control device executes air-fuel ratio feedback control.
FIG. 8 is a diagram schematically showing a catalyst model adopted by the exhaust gas purification apparatus of the present invention.
FIG. 9 is a diagram showing a balance of a specific component related to an oxygen storage / release reaction in a specific region of the catalyst model employed in the exhaust gas purification apparatus of the present invention when focusing on the specific region.
FIG. 10 is a map for obtaining a distribution of each maximum oxygen storage amount for each block of the catalyst from a maximum oxygen storage amount of the entire catalyst which is a target of a catalyst model employed in the exhaust purification device of the present invention.
FIG. 11 is a schematic diagram when a catalyst model adopted by the exhaust gas purification apparatus of the present invention is applied to a first catalyst.
FIG. 12 is a flowchart showing a routine for performing a start determination of an update process of a basic fuel injection amount correction amount executed by the CPU shown in FIG. 1;
FIG. 13 is a flowchart illustrating a routine for setting a target air-fuel ratio executed by the CPU illustrated in FIG. 1;
FIG. 14 is a flowchart showing a routine for calculating a fuel injection amount and the like executed by a CPU shown in FIG. 1;
15 is a flowchart showing a routine for calculating a feedback control amount when the conventional air-fuel ratio control device executed by the CPU shown in FIG. 1 is applied.
FIG. 16 is a flowchart showing a routine for determining whether to end a process of updating a basic fuel injection amount correction amount executed by the CPU shown in FIG. 1;
17 is a flowchart showing a routine executed by the CPU shown in FIG. 1 for calculating a feedback control amount when the exhaust gas purification device (air-fuel ratio control device) of the present invention is applied.
FIG. 18 is a flowchart showing a routine for determining whether to start a maximum oxygen storage amount acquisition control executed by the CPU shown in FIG. 1;
FIG. 19 is a flowchart showing a first mode routine executed by the CPU shown in FIG. 1;
20 is a flowchart showing a second mode routine executed by the CPU shown in FIG. 1;
FIG. 21 is a flowchart showing a routine for calculating an oxygen storage amount executed by the CPU shown in FIG. 1;
FIG. 22 is a flowchart showing a routine for calculating a maximum oxygen storage amount executed by the CPU shown in FIG. 1;
23 is a flowchart showing a routine executed by the CPU shown in FIG. 1 for calculating an outflow oxygen amount, an oxygen storage amount, and the like for each block of the first catalyst.
24 is a flowchart showing a routine executed by the CPU shown in FIG. 1 for clearing the oxygen storage amount of each block of the first catalyst.
[Explanation of symbols]
10 internal combustion engine, 25 combustion chamber, 39 injector, 52 exhaust pipe (exhaust pipe), 53 three-way catalyst (first catalyst), 54 three-way catalyst (second catalyst), 66 air upstream Fuel ratio sensor 67 67 Downstream air-
Claims (6)
前記触媒よりも上流の前記排気通路に配設された空燃比センサと、
少なくとも前記空燃比センサの出力値と所定の目標値との偏差に基づく値の時間積分値に基づいてフィードバック制御量を算出するフィードバック制御量算出手段と、
前記フィードバック制御量に基づいて前記触媒に流入するガスの空燃比をフィードバック制御する空燃比制御手段と、を備えた内燃機関の排気浄化装置において、
前記フィードバック制御量算出手段は、前記偏差に基づく値の時間積分値が所定の範囲を超える場合、同偏差に基づく値の時間積分値の代わりに同所定の範囲内の値に基づいて前記フィードバック制御量を算出するように構成された内燃機関の排気浄化装置。A catalyst disposed in an exhaust passage of the internal combustion engine;
An air-fuel ratio sensor disposed in the exhaust passage upstream of the catalyst,
Feedback control amount calculating means for calculating a feedback control amount based on a time integration value of at least a value based on a deviation between an output value of the air-fuel ratio sensor and a predetermined target value,
Air-fuel ratio control means for feedback-controlling the air-fuel ratio of the gas flowing into the catalyst based on the feedback control amount,
When the time integrated value of the value based on the deviation exceeds a predetermined range, the feedback control amount calculating means performs the feedback control based on the value within the predetermined range instead of the time integrated value of the value based on the deviation. An exhaust gas purification device for an internal combustion engine configured to calculate an amount.
少なくとも前記偏差に基づく値の時間積分値に基づいて前記触媒が吸蔵している酸素量である酸素吸蔵量を推定する酸素吸蔵量推定手段を備え、
前記フィードバック制御量算出手段は、前記偏差に基づく値の時間積分値に基づく値としての前記推定された酸素吸蔵量に基づいて前記フィードバック制御量を算出するように構成された内燃機関の排気浄化装置。An exhaust gas purification device for an internal combustion engine according to claim 1,
An oxygen storage amount estimating unit that estimates an oxygen storage amount that is an amount of oxygen stored by the catalyst based on a time integration value of at least a value based on the deviation,
The feedback control amount calculating means is configured to calculate the feedback control amount based on the estimated oxygen storage amount as a value based on a time integration value of the value based on the deviation. .
前記フィードバック制御量算出手段は、所定の目標酸素吸蔵量と前記推定された酸素吸蔵量との偏差に基づいて前記フィードバック制御量を算出するように構成された内燃機関の排気浄化装置。The exhaust gas purification device for an internal combustion engine according to claim 2,
The exhaust gas purification device for an internal combustion engine, wherein the feedback control amount calculating means is configured to calculate the feedback control amount based on a deviation between a predetermined target oxygen storage amount and the estimated oxygen storage amount.
前記酸素吸蔵量推定手段は、前記触媒内における酸素吸蔵反応の速度を考慮して構築された触媒モデルに基づいて前記酸素吸蔵量を推定するように構成された内燃機関の排気浄化装置。The exhaust gas purification apparatus for an internal combustion engine according to claim 2 or 3,
The exhaust gas purifying apparatus for an internal combustion engine, wherein the oxygen storage amount estimating means is configured to estimate the oxygen storage amount based on a catalyst model constructed in consideration of a speed of an oxygen storage reaction in the catalyst.
前記内燃機関の運転状態に応じて目標空燃比に相当する前記所定の目標値を設定する目標値設定手段を備え、
前記フィードバック制御量算出手段は、
更に前記空燃比センサの出力値と前記所定の目標値との偏差に基づく値に基づいて前記フィードバック制御量を算出するように構成されているとともに、
前記所定の目標値が理論空燃比と異なる空燃比に相当する値に設定されている場合、前記所定の目標酸素吸蔵量と前記推定された酸素吸蔵量との偏差の代わりに前記空燃比センサの出力値と同所定の目標値との偏差に基づく値の時間積分値に基づいて前記フィードバック制御量を算出するように構成された内燃機関の排気浄化装置。An exhaust gas purification device for an internal combustion engine according to claim 3 or claim 4,
Target value setting means for setting the predetermined target value corresponding to a target air-fuel ratio according to the operating state of the internal combustion engine,
The feedback control amount calculation means,
Furthermore, it is configured to calculate the feedback control amount based on a value based on a deviation between the output value of the air-fuel ratio sensor and the predetermined target value,
When the predetermined target value is set to a value corresponding to an air-fuel ratio different from the stoichiometric air-fuel ratio, the air-fuel ratio sensor replaces the deviation between the predetermined target oxygen storage amount and the estimated oxygen storage amount. An exhaust emission control device for an internal combustion engine configured to calculate the feedback control amount based on a time integration value of a value based on a deviation between an output value and the predetermined target value.
前記排気浄化装置が適用される内燃機関は、指示に応じて燃料噴射を行う燃料噴射手段を備えていて、
前記空燃比制御手段は、前記内燃機関に供給される混合気の空燃比を前記所定の目標値に相当する目標空燃比とするために必要な基本燃料噴射量を同内燃機関の運転状態に応じて決定するとともに同基本燃料噴射量を前記フィードバック制御量に基づいて補正することにより得られる燃料噴射量の燃料の噴射指示を前記燃料噴射手段に対して行うことで前記触媒に流入するガスの空燃比をフィードバック制御するように構成されていて、
前記フィードバック制御量算出手段が前記空燃比センサの出力値と前記所定の目標値との偏差に基づく値の時間積分値が前記所定の範囲を超えるか否かに拘わらず同偏差に基づく値の時間積分値に基づいて前記フィードバック制御量を算出した場合であって、且つ前記内燃機関が定常運転状態にある場合における同偏差に基づく値の時間積分値に基づいて前記基本燃料噴射量を補正する基本燃料噴射量補正手段を備えた内燃機関の排気浄化装置。An exhaust gas purification apparatus for an internal combustion engine according to any one of claims 1 to 5, wherein
The internal combustion engine to which the exhaust purification device is applied includes a fuel injection unit that performs fuel injection according to an instruction,
The air-fuel ratio control means determines a basic fuel injection amount necessary for setting an air-fuel ratio of an air-fuel mixture supplied to the internal combustion engine to a target air-fuel ratio corresponding to the predetermined target value according to an operating state of the internal combustion engine. The fuel injection amount is obtained by correcting the basic fuel injection amount based on the feedback control amount and instructing the fuel injection means to inject the fuel into the fuel injection means. It is configured to feedback control the fuel ratio,
The feedback control amount calculating means determines whether the time integration value of the value based on the deviation between the output value of the air-fuel ratio sensor and the predetermined target value exceeds the predetermined range regardless of whether the time integration value is based on the deviation. Basically correcting the basic fuel injection amount based on a time integration value of a value based on the deviation when the feedback control amount is calculated based on an integrated value and the internal combustion engine is in a steady operation state. An exhaust gas purification device for an internal combustion engine, comprising a fuel injection amount correction means.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2003039170A JP2004251123A (en) | 2003-02-18 | 2003-02-18 | Exhaust emission control device for internal combustion engine |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2003039170A JP2004251123A (en) | 2003-02-18 | 2003-02-18 | Exhaust emission control device for internal combustion engine |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2004251123A true JP2004251123A (en) | 2004-09-09 |
Family
ID=33023413
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2003039170A Pending JP2004251123A (en) | 2003-02-18 | 2003-02-18 | Exhaust emission control device for internal combustion engine |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2004251123A (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015105012A1 (en) | 2014-01-10 | 2015-07-16 | Toyota Jidosha Kabushiki Kaisha | Control System of Internal Combustion Engine |
WO2015177996A1 (en) | 2014-05-23 | 2015-11-26 | Toyota Jidosha Kabushiki Kaisha | Control system of internal combustion engine |
-
2003
- 2003-02-18 JP JP2003039170A patent/JP2004251123A/en active Pending
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015105012A1 (en) | 2014-01-10 | 2015-07-16 | Toyota Jidosha Kabushiki Kaisha | Control System of Internal Combustion Engine |
US10221789B2 (en) | 2014-01-10 | 2019-03-05 | Toyota Jidosha Kabushiki Kaisha | Control system of internal combustion engine |
WO2015177996A1 (en) | 2014-05-23 | 2015-11-26 | Toyota Jidosha Kabushiki Kaisha | Control system of internal combustion engine |
US10132257B2 (en) | 2014-05-23 | 2018-11-20 | Toyota Jidosha Kabushiki Kaisha | Control system of internal combustion engine |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3846375B2 (en) | Catalyst degradation judgment method | |
JP3963130B2 (en) | Catalyst deterioration judgment device | |
JP4973807B2 (en) | Air-fuel ratio control device for internal combustion engine | |
JP3972748B2 (en) | Exhaust gas purification device for internal combustion engine | |
JP4380745B2 (en) | Device for acquiring degree of deterioration of catalyst of internal combustion engine | |
JP2009162139A (en) | Air-fuel ratio control device for internal combustion engine | |
JP3922091B2 (en) | Air-fuel ratio control device for internal combustion engine | |
JP2007239700A (en) | Fuel injection control device for internal combustion engine | |
JP2004360591A (en) | Exhaust emission control device for internal combustion engine | |
JP2007113437A (en) | Air-fuel ratio control device of internal combustion engine | |
US8949000B2 (en) | Control device and control method for internal combustion engine | |
JP2004003405A (en) | Catalyst degradation determining device | |
JP3815386B2 (en) | Catalyst degradation judgment method | |
JP2004251123A (en) | Exhaust emission control device for internal combustion engine | |
JP2004076681A (en) | Catalyst deterioration determination device | |
JP2003314334A (en) | Air-fuel ratio control device for internal combustion engine | |
JP2004263591A (en) | Exhaust emission control device for internal combustion engine | |
JP2005009401A (en) | Exhaust emission control device for internal combustion engine | |
JP4239521B2 (en) | Exhaust gas purification device for internal combustion engine | |
JP2004300981A (en) | Device for judging catalyst degradation | |
JP4103522B2 (en) | Exhaust gas purification device for internal combustion engine | |
JP3879596B2 (en) | Air-fuel ratio sensor state determination device | |
JP2004092472A (en) | Exhaust emission control device for internal combustion engine | |
JP2009197683A (en) | Air-fuel ratio controller of internal combustion engine | |
JP2013002374A (en) | Apparatus for estimating maximum oxygen occlusion amount of catalyst |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20060125 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20071203 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20071204 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20080121 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20080610 |