JP4437626B2 - 内燃機関の空燃比制御装置 - Google Patents
内燃機関の空燃比制御装置 Download PDFInfo
- Publication number
- JP4437626B2 JP4437626B2 JP2001143800A JP2001143800A JP4437626B2 JP 4437626 B2 JP4437626 B2 JP 4437626B2 JP 2001143800 A JP2001143800 A JP 2001143800A JP 2001143800 A JP2001143800 A JP 2001143800A JP 4437626 B2 JP4437626 B2 JP 4437626B2
- Authority
- JP
- Japan
- Prior art keywords
- value
- fuel ratio
- air
- deterioration
- output
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F01—MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
- F01N—GAS-FLOW SILENCERS OR EXHAUST APPARATUS FOR MACHINES OR ENGINES IN GENERAL; GAS-FLOW SILENCERS OR EXHAUST APPARATUS FOR INTERNAL COMBUSTION ENGINES
- F01N11/00—Monitoring or diagnostic devices for exhaust-gas treatment apparatus, e.g. for catalytic activity
- F01N11/007—Monitoring or diagnostic devices for exhaust-gas treatment apparatus, e.g. for catalytic activity the diagnostic devices measuring oxygen or air concentration downstream of the exhaust apparatus
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F01—MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
- F01N—GAS-FLOW SILENCERS OR EXHAUST APPARATUS FOR MACHINES OR ENGINES IN GENERAL; GAS-FLOW SILENCERS OR EXHAUST APPARATUS FOR INTERNAL COMBUSTION ENGINES
- F01N2550/00—Monitoring or diagnosing the deterioration of exhaust systems
- F01N2550/02—Catalytic activity of catalytic converters
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
Landscapes
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Combustion & Propulsion (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Exhaust Gas After Treatment (AREA)
- Electrical Control Of Air Or Fuel Supplied To Internal-Combustion Engine (AREA)
- Combined Controls Of Internal Combustion Engines (AREA)
Description
【発明の属する技術分野】
本発明は、内燃機関の空燃比制御装置に関し、より詳しくは排ガス浄化用の触媒装置の劣化状態を評価することができる空燃比制御装置に関する。
【0002】
【従来の技術】
内燃機関の排気通路に設けられた三元触媒等の排ガス浄化用の触媒装置の適正な浄化性能を確保するために、内燃機関で燃焼させる混合気の空燃比を制御する技術が従来より知られている。例えば本願出願人が、特開平9−324681号公報や特開平11−153051号公報に提案した技術が知られている。この技術では、触媒装置の下流側にO2センサ(酸素濃度センサ)を配置し、そのO2センサの出力をあらかじめ定めた所定の目標値(一定値)に収束させるように、フィードバック制御(具体的にはスライディングモード制御)の処理により混合気の空燃比を操作するものである。
【0003】
一方、前記触媒装置は、内燃機関の累積的な運転に伴い次第に劣化し、排ガスの浄化能力が低下していくので、その劣化がある程度進行した段階で新品の触媒装置に交換する等の処置を採る必要がある。このため、触媒装置の劣化状態を適正に評価し得る技術が望まれ、この種の技術としては、例えば特許公報第2526640号や、特開平7−19033号公報に見られる技術が従来より知られている。これらの従来技術では、内燃機関で燃焼させる混合気の空燃比を積極的にリーン側からリッチ側、あるいはリッチ側からリーン側に変化させ、このとき、触媒装置の下流側に設けたO2センサの出力が反転するまでの時間や、該O2センサの出力の反転周期を計測する。そして、該従来技術は、その計測値に基づいて触媒装置の劣化状態を評価するものである。
【0004】
しかし、これらの従来技術は、触媒装置の劣化状態を評価するために積極的に空燃比を変化させる必要があることから、該劣化状態の評価の際には、触媒装置の適正な浄化性能を確保することが困難である。換言すれば、触媒装置の適正な浄化性能を確保しながら、触媒装置の劣化状態を評価することが困難なものである。
【0005】
そこで、本願出願人は、前述のように触媒装置の下流側のO2センサの出力を所定の目標値に収束させるように空燃比を操作しつつ、換言すれば、触媒装置の良好な浄化性能を確保しつつ、触媒装置の劣化状態を評価することが可能な装置の開発を試みている(例えば本願出願人による特願2000-139860号)。この技術は、上記のような空燃比制御を実行している際における触媒装置の下流側のO2センサの出力のデータから、ある適当なパラメータを求めたとき、該パラメータの値が触媒装置の劣化の進行に伴い特徴的な特性を持つことに着目し、そのパラメータの値に基づき触媒装置の劣化状態を評価するものである。このような技術によれば、触媒装置の良好な浄化性能を確保しつつ、触媒装置の劣化状態を評価することが可能となる。
【0006】
ところで、前述のように触媒装置の下流側のO2センサの出力を所定の目標値に収束させるように空燃比を操作する技術(特許公報第2526640号や特開平7−19033号公報等)では、下流側のO2センサの出力の目標値を基本的には所定の一定値としているが、触媒装置の種類等によっては、該触媒装置の好適な浄化性能を確保する上で、内燃機関の運転状態に応じてO2センサの出力の目標値を可変的に設定することが好ましいと考えられる場合もある。
【0007】
そして、このように触媒装置の下流側のO2センサの出力の目標値を可変的に設定して空燃比制御を行なうものでは、その空燃比制御と並行して触媒装置の劣化状態の評価を行なうようにしたとき、その劣化状態の評価の信頼性に関して次のような不都合を生じる虞れがある。
【0008】
すなわち、酸素濃度に対するO2センサの出力は非線形性を有するため(図2の実線aを参照)、O2センサの出力の目標値が異なると、該O2センサの出力を目標値に収束制御している際の該O2センサの出力の変化の形態が異なるものになると考えられる。この結果、該O2センサの出力のデータから一義的に求めた前記パラメータの値が、触媒装置の劣化状態以外の要因の影響を受けるようになり、該パラメータの値に基づく触媒装置の劣化状態の評価の信頼性が低下する虞れがある。
【0009】
【発明が解決しようとする課題】
本発明はかかる背景に鑑みてなされたものであり、触媒装置の下流側のO2センサの出力の目標値を内燃機関の運転状態に応じて可変的に設定して、該O2センサの出力をその目標値に収束させるように空燃比を操作するものにおいて、該目標値の可変化の影響を補償しつつ、触媒装置の劣化状態の評価を適正に行うことができる内燃機関の空燃比制御装置を提供することを目的とする。
【0010】
【課題を解決するための手段】
本発明の内燃機関の空燃比制御装置はかかる目的を達成するために、二つの態様があり、その第1の態様は、内燃機関の排気通路に設けた触媒装置の下流側に配置された酸素濃度センサと、該酸素濃度センサの出力を内燃機関の運転状態に応じて設定された所定の目標値に収束させるように該内燃機関で燃焼させる混合気の空燃比を操作する空燃比操作手段と、該空燃比操作手段の処理の実行中における前記酸素濃度センサの出力のデータからあらかじめ定めた所定のアルゴリズムにより前記触媒装置の劣化状態を評価するための劣化評価用パラメータを生成するパラメータ生成手段とを備えた内燃機関の空燃比制御装置において、前記劣化評価用パラメータは、前記酸素濃度センサの出力の時系列データを変数成分として表した所定の劣化評価用線形関数の値のばらつき度合いを表すデータであり、前記劣化評価用パラメータを前記所定の目標値又は前記酸素濃度センサの出力の平均値に応じて補正し、その補正後の劣化評価用パラメータの値に基づき前記触媒装置の劣化状態を評価する劣化状態評価手段を備えたことを特徴とするものである。
【0011】
また、第2の態様は、内燃機関の排気通路に設けた触媒装置の下流側に配置された酸素濃度センサと、該酸素濃度センサの出力を内燃機関の運転状態に応じて設定された所定の目標値に収束させるように該内燃機関で燃焼させる混合気の空燃比を操作する空燃比操作手段と、該空燃比操作手段の処理の実行中における前記酸素濃度センサの出力のデータからあらかじめ定めた所定のアルゴリズムにより前記触媒装置の劣化状態を評価するための劣化評価用パラメータを生成するパラメータ生成手段と、該劣化評価用パラメータの値を所定の判定値と比較することにより前記触媒装置の劣化状態を評価する劣化状態評価手段とを備えた内燃機関の空燃比制御装置において、前記劣化評価用パラメータは、前記酸素濃度センサの出力の時系列データを変数成分として表した所定の劣化評価用線形関数の値のばらつき度合いを表すデータであり、前記劣化状態評価手段は、前記判定値を前記所定の目標値又は前記酸素濃度センサの出力の平均値に応じて設定する手段を備え、その設定した判定値と前記劣化評価用パラメータの値を比較することにより前記触媒装置の劣化状態を評価することを特徴とするものである。
【0012】
かかる本発明によれば、前記酸素濃度センサの出力の目標値又は該酸素濃度センサの出力の平均値に応じて、換言すれば、前記空燃比操作手段の処理によって酸素濃度センサの出力がどのようなレベル値に制御されるかに応じて、前記酸素濃度センサの出力データから生成される劣化評価用パラメータを補正し(第1の態様)、あるいは、該劣化評価用パラメータと比較する前記判定値を設定する(第2の態様)ので、酸素濃度センサの出力の目標値を内燃機関の運転状態に応じて可変的に設定しても、その目標値によらずに、触媒装置の劣化状態を適正に評価することが可能となる。
【0013】
すなわち、本発明の第1の態様では、劣化評価用パラメータを前記目標値又は酸素濃度センサの出力の平均値に応じて補正することで、その補正後の劣化評価用パラメータの値を、酸素濃度センサの出力の目標値によらずに触媒装置の劣化状態のみに応じて変化するようなものにすることができ、その結果、該補正後の劣化評価用パラメータの値に基づいて、触媒装置の劣化状態の評価を適正に行なうことができる。また、本発明の第2の態様では、劣化評価用パラメータの値は、酸素濃度センサの出力の目標値の影響を受けるものの、該劣化評価用パラメータの値と比較する判定値を前記目標値又は酸素濃度センサの出力の平均値に応じて設定することにより、前記第1の態様と同等の効果を奏することができる。そして、本発明では、前記空燃比操作手段の処理の実行中、すなわち、前記酸素濃度センサの出力の前記目標値への収束制御を行なっている状態における酸素濃度センサの出力のデータから前記劣化評価用パラメータを生成するので、触媒装置の良好な浄化性能を確保しながら触媒装置の劣化状態を評価することができる。
【0014】
さらに本発明では、前記劣化評価用パラメータは、前記酸素濃度センサの出力の時系列データを変数成分として表した所定の劣化評価用線形関数の値のばらつき度合いを表すデータである。
【0015】
すなわち、酸素濃度センサの出力の前記目標値への収束制御を行っている状態において、該酸素濃度センサの時系列データから、その時系列データを変数成分とする、ある適当な所定の線形関数(酸素濃度センサの出力の時系列データの線形結合の形に表現されるような関数)の値を求めたとき、その線形関数の値は、触媒装置の劣化の進行度合いとの間に特徴的な相関性を呈する傾向がある。例えば、触媒装置がほぼ新品状態であるときには、前記線形関数の値は、ある所定値の近傍に集積するが、触媒装置の劣化が進行していくと、前記線形関数の値は、上記所定値から離間した値を採り易くなる傾向がある。つまり、触媒装置の劣化の進行に伴い該線形関数の値のばらつき度合いが大きくなる。
【0016】
従って、前記線形関数の値のばらつき度合いを表すデータを劣化評価用パラメータとして生成することで、該劣化評価用パラメータの値と触媒装置の劣化状態との相関性を高めることができ、ひいては、触媒装置の劣化状態の評価の信頼性を高めることができる。
【0017】
この場合、前記劣化評価用パラメータは、例えば前記劣化評価用線形関数の値と所定値との偏差の二乗値や該偏差の絶対値等であってもよいが、好ましくは、前記パラメータ生成手段は、前記劣化評価用線形関数の値の時系列データの各データ値と該劣化評価用線形関数の値の中心値としてあらかじめ定めた所定値との偏差の二乗値又は絶対値にローパス特性のフィルタリング処理を施すことにより前記劣化評価用パラメータを求める。
【0018】
このように前記偏差の二乗値又は絶対値にローパス特性のフィルタリング処理を施して劣化評価用パラメータを求めたとき、その劣化評価用パラメータの値は、劣化評価用線形関数の値のばらつき度合いを表すものとして適正なものとなると共に、触媒装置の劣化の進行に伴い単調的に増加していくようになり、該触媒装置の劣化状態との相関性が明白になる。従って、この劣化評価用パラメータの値に基づいて、触媒装置の劣化状態の評価を高い信頼性で確実に行うことが可能となる。
【0019】
また、前述のように劣化評価用パラメータを求めるために前記劣化評価用線形関数を用いる本発明では、前記空燃比操作手段は、前記酸素濃度センサの出力を前記目標値に収束させるようにスライディングモード制御の処理により前記空燃比を操作するための操作量を逐次生成しつつ、該操作量に応じて前記空燃比を操作する手段であり、前記劣化評価用線形関数は、前記スライディングモード制御の処理に用いる切換関数に応じて定められた線形関数であることが好適である。
【0020】
すなわち、前記劣化評価用線形関数の値の前述のような触媒装置の劣化状態に応じた傾向は、前記空燃比操作手段により空燃比を操作するための前記操作量(例えば目標空燃比等)をフィードバック制御の一手法であるスライディングモード制御の処理により生成した場合に、好適に顕著なものとなりやすい。そして、このようにスライディングモード制御の処理を用いて空燃比を操作する場合には、特に、触媒装置の劣化状態との相関性が高い劣化評価用線形関数は該スライディングモード制御の処理に用いる切換関数に密接に関連しており、該切換関数に応じて定まる線形関数を劣化評価用線形関数とすることが好適である。
【0021】
この場合、より具体的には、前記スライディングモード制御の処理では、例えば前記酸素濃度センサの出力と前記目標値との偏差の時系列データを変数成分として表した線形関数を前記切換関数として用いる。そして、このような切換関数をスライディングモード制御の処理に用いる場合には、前記劣化評価用線形関数は、その変数成分に係る係数値を前記切換関数の変数成分に係る係数値と同一とした線形関数とすることが好適である。尚、この線形関数は、スライディングモード制御用の前記切換関数そのものであってもよい。
【0022】
このようにスライディングモード制御用の切換関数に応じて定まる線形関数を前記劣化評価用線形関数とすることで、該劣化評価用線形関数の値と触媒装置の劣化状態との相関性が顕著なものとなり、該劣化評価用線形関数の値に基づく触媒装置の劣化状態の評価を適正に行うことができる。
【0023】
【発明の実施の形態】
本発明の第1実施形態を図1〜図20を参照して説明する。
【0024】
図1は本実施形態の装置の全体的なシステム構成を示すブロック図である。図中、1は、例えば自動車やハイブリッド車にその車両の推進源として搭載された4気筒のエンジン(内燃機関)である。このエンジン1の各気筒毎に燃料及び空気の混合気の燃焼により生成される排ガスは、エンジン1の近傍で共通の排気管2(排気通路)に集合され、該排気管2を介して大気中に放出される。排気管2には、排ガスを浄化するために、例えば三元触媒により構成された二つの触媒装置3,4が上流側から順に介装されている。
【0025】
尚、本実施形態で劣化状態を評価する触媒装置は、上流側の触媒装置3であり、下流側の触媒装置4は省略してもよい。
【0026】
本実施形態の装置では、基本的には、触媒装置3の最適な浄化性能を確保するようにエンジン1で燃焼させる混合気の空燃比を操作する。また、この空燃比制御を行いながら、触媒装置3の劣化状態を評価する。
【0027】
そして、このような処理を行うために、触媒装置3の上流側(より詳しくはエンジン1の各気筒毎の排ガスの集合箇所)で排気管2に設けられた空燃比センサ5と、触媒装置3の下流側(触媒装置4の上流側)で排気管2に設けられたO2センサ(酸素濃度センサ)6と、これらのセンサ5,6の出力等に基づき後述の制御処理や触媒装置3の劣化状態の評価を行う制御ユニット7とを備えている。
【0028】
尚、制御ユニット7には、前記空燃比センサ5やO2センサ6の出力の他に、エンジン1の運転状態を検出するための図示しない回転数センサや吸気圧センサ、冷却水温センサ等、各種のセンサの出力が与えられる。
【0029】
O2センサ6は、触媒装置3を通過した排ガス中の酸素濃度に応じたレベルの出力VO2/OUT(酸素濃度の検出値を示す出力)を生成する通常的なO2センサである。ここで、排ガス中の酸素濃度は、その排ガスを燃焼により生成した混合気の空燃比に応じたものとなる。そして、このO2センサ6の出力VO2/OUTは、図2に実線aで示す如く、排ガス中の酸素濃度に対応する空燃比が理論空燃比近傍の範囲Δに存するような状態で、該排ガス中の酸素濃度に対して高感度な変化を生じるものとなる。また、その範囲Δを逸脱した空燃比に対応する酸素濃度では、O2センサ6の出力VO2/OUTは飽和して、ほぼ一定のレベルとなる。この場合、上記範囲Δにおいて、酸素濃度に対するO2センサ6の出力VO2/OUTの感度(酸素濃度の変化に対する出力VO2/OUTの変化率)は、該出力VO2/OUTの上限レベル及び下限レベルの間の中央レベルVO2/OUT1の近傍で最大となり、この中央レベルVO2/OUTよりも上側あるいは下側のレベル(例えば図2のVO2/OUT2や、VO2/OUT3のレベル)では、該感度は若干低下する。
【0030】
空燃比センサ5は、触媒装置3に進入する排ガスの酸素濃度により把握される空燃比(以下、この空燃比を単にエンジン1の空燃比ということがある)の検出値を表す出力KACTを生成するものである。この空燃比センサ5は、例えば本願出願人が特開平4−369471号公報にて詳細に説明した広域空燃比センサにより構成されたものであり、図2に実線bで示す如く、O2センサ6よりも排ガス中の酸素濃度の広範囲にわたってそれに比例したレベルの出力を生成するものである。換言すれば、該空燃比センサ5(以下、LAFセンサ5という)は、排ガス中の酸素濃度に対応した空燃比の広範囲にわたってそれに比例したレベルの出力KACTを生成するものである。
【0031】
制御ユニット7は、エンジン1の空燃比を規定する操作量としての目標空燃比KCMD(LAFセンサ5で検出される空燃比の目標値)を算出するための処理、並びに触媒装置3の劣化状態を評価するための処理を担う制御ユニット7a(以下、排気側制御ユニット7aという)と、上記目標空燃比KCMDに応じてエンジン1の燃料噴射量(燃料供給量)を調整することでエンジン1で燃焼させる混合気の空燃比を操作する処理を担う制御ユニット7b(以下、機関側制御ユニット7bという)とから構成されている。
【0032】
これらの制御ユニット7a,7bはマイクロコンピュータを用いて構成されたものであり、それぞれの処理を所定の制御サイクルで実行する。ここで、本実施形態では、排気側制御ユニット7aがその処理(目標空燃比KCMDの生成処理や触媒装置3の劣化状態の評価処理)を実行する制御サイクルは、触媒装置3に起因した後述の無駄時間や演算負荷等を考慮して、あらかじめ定めた一定周期(例えば30〜100ms)としている。
【0033】
また、機関側制御ユニット7bが実行する処理(燃料噴射量の調整処理)は、エンジン1の回転数(詳しくはエンジン1の燃焼サイクル)に同期させて行う必要がある。このため、機関側制御ユニット7bがその処理を実行する制御サイクルは、エンジン1のクランク角周期(所謂TDC)に同期した周期としている。
【0034】
尚、排気側制御ユニット7aの制御サイクルの一定周期は、前記クランク角周期(TDC)よりも長いものとされている。
【0035】
前記機関側制御ユニット7bを図1を参照してさらに説明する。機関側制御ユニット7bは、その機能的構成として、エンジン1への基本燃料噴射量Timを求める基本燃料噴射量算出部8と、基本燃料噴射量Timを補正するための第1補正係数KTOTAL及び第2補正係数KCMDMをそれぞれ求める第1補正係数算出部9及び第2補正係数算出部10とを具備する。
【0036】
前記基本燃料噴射量算出部8は、エンジン1の回転数NEと吸気圧PBとから、それらにより規定されるエンジン1の基準の燃料噴射量(燃料供給量)をあらかじめ設定されたマップを用いて求め、その基準の燃料噴射量をエンジン1の図示しないスロットル弁の有効開口面積に応じて補正することで基本燃料噴射量Timを算出するものである。
【0037】
また、第1補正係数算出部9が求める第1補正係数KTOTALは、エンジン1の排気還流率(エンジン1の吸入空気中に含まれる排ガスの割合)や、エンジン1の図示しないキャニスタのパージ時にエンジン1に供給される燃料のパージ量、エンジン1の冷却水温、吸気温等を考慮して前記基本燃料噴射量Timを補正するためのものである。
【0038】
また、第2補正係数算出部10が求める第2補正係数KCMDMは、排気側制御ユニット7aが後述の如く算出する目標空燃比KCMDに対応してエンジン1へ流入する燃料の冷却効果による吸入空気の充填効率を考慮して基本燃料噴射量Timを補正するためのものである。
【0039】
これらの第1補正係数KTOTAL及び第2補正係数KCMDMによる基本燃料噴射量Timの補正は、第1補正係数KTOTAL及び第2補正係数KCMDMを基本燃料噴射量Timに乗算することで行われ、この補正によりエンジン1の要求燃料噴射量Tcylが得られる。
【0040】
尚、前記基本燃料噴射量Timや、第1補正係数KTOTAL、第2補正係数KCMDMのより具体的な算出手法は、特開平5−79374号公報等に本願出願人が開示しているので、ここでは詳細な説明を省略する。
【0041】
機関側制御ユニット7bは、上記の機能的構成の他、さらに、排気側制御ユニット7aが算出する目標空燃比KCMDにLAFセンサ5の出力KACT(エンジン1の空燃比の検出値)を収束させるようにフィードバック制御によりエンジン1の燃料噴射量を調整することでエンジン1で燃焼させる混合気の空燃比を操作するフィードバック制御部14を備えている。
【0042】
このフィードバック制御部14は、本実施形態では、エンジン1の各気筒の全体的な空燃比をフィードバック制御する大局的フィードバック制御部15と、エンジン1の各気筒毎の空燃比をフィードバック制御する局所的フィードバック制御部16とに分別される。
【0043】
前記大局的フィードバック制御部15は、LAFセンサ5の出力KACTが前記目標空燃比KCMDに収束するように、前記要求燃料噴射量Tcylを補正する(要求燃料噴射量Tcylに乗算する)フィードバック補正係数KFBを逐次求めるものである。
【0044】
この大局的フィードバック制御部15は、LAFセンサ5の出力KACTと目標空燃比KCMDとの偏差に応じて周知のPID制御を用いて前記フィードバック補正係数KFBとしてのフィードバック操作量KLAFを生成するPID制御器17と、LAFセンサ5の出力KACTと目標空燃比KCMDとからエンジン1の運転状態の変化や特性変化等を考慮して前記フィードバック補正係数K FBを規定するフィードバック操作量KSTRを適応的に求める適応制御器18(図ではSTRと称している)とをそれぞれ独立的に具備している。
【0045】
ここで、本実施形態では、前記PID制御器17が生成するフィードバック操作量KLAFは、LAFセンサ5の出力KACT(空燃比の検出値)が目標空燃比KCMDに一致している状態で「1」となり、該操作量KLAFをそのまま前記フィードバック補正係数KFBとして使用できるようになっている。一方、適応制御器18が生成するフィードバック操作量KSTRはLAFセンサ5の出力KACTが目標空燃比KCMDに一致する状態で「目標空燃比KCMD」となるものである。このため、該フィードバック操作量KSTRを除算処理部19で目標空燃比KCMDにより除算してなるフィードバック操作量kstr(=KSTR/KCMD)が前記フィードバック補正係数KFBとして使用できるようになっている。
【0046】
大局的フィードバック制御部15は、PID制御器17により生成されるフィードバック操作量KLAFと、適応制御器18が生成するフィードバック操作量KSTRを目標空燃比KCMDにより除算してなるフィードバック操作量kstrとを切換部20で適宜、択一的に選択する。そして、いずれか一方のフィードバック操作量KLAF又はkstrを前記フィードバック補正係数KFBとして使用し、該補正係数KFBを前記要求燃料噴射量Tcylに乗算することにより該要求燃料噴射量Tcylを補正する。尚、かかる大局的フィードバック制御部15(特に適応制御器18)については後にさらに詳細に説明する。
【0047】
前記局所的フィードバック制御部16は、LAFセンサ5の出力KACTから各気筒毎の実空燃比#nA/F(n=1,2,3,4)を推定するオブザーバ21と、このオブザーバ21により推定された各気筒毎の実空燃比#nA/Fから各気筒毎の空燃比のばらつきを解消するよう、PID制御を用いて各気筒毎の燃料噴射量のフィードバック補正係数#nKLAFをそれぞれ求める複数(気筒数個)のPID制御器22とを具備する。
【0048】
ここで、オブザーバ21は、それを簡単に説明すると、各気筒毎の実空燃比#nA/Fの推定を次のように行うものである。すなわち、エンジン1からLAFセンサ5の箇所(各気筒毎の排ガスの集合部)にかけての系を、エンジン1の各気筒毎の実空燃比#nA/FからLAFセンサ5で検出される空燃比を生成する系と考える。そして、この系を、LAFセンサ5の検出応答遅れ(例えば一次遅れ)や、LAFセンサ5で検出される空燃比に対するエンジン1の各気筒毎の空燃比の時間的寄与度を考慮してモデル化する。そして、そのモデルの基で、LAFセンサ5の出力KACTから、逆算的に各気筒毎の実空燃比#nA/Fを推定する。
【0049】
尚、このようなオブザーバ21は、本願出願人が例えば特開平7−83094号公報に詳細に開示しているので、ここでは、さらなる説明を省略する。
【0050】
また、局所的フィードバック制御部16の各PID制御器22は、LAFセンサ5の出力KACTを、前回の制御サイクルで各PID制御器22により求められたフィードバック補正係数#nKLAFの全気筒についての平均値により除算してなる値を各気筒の空燃比の目標値とする。さらに、その目標値とオブザーバ21により求められた各気筒毎の実空燃比#nA/Fの推定値との偏差が解消するように、今回の制御サイクルにおける各気筒毎のフィードバック補正係数#nKLAFを求める。
【0051】
そして、局所的フィードバック制御部16は、前記要求燃料噴射量Tcylに大局的フィードバック制御部15のフィードバック補正係数KFBを乗算してなる値に、各気筒毎のフィードバック補正係数#nKLAFを乗算することで、各気筒の出力燃料噴射量#nTout(n=1,2,3,4)を求める。
【0052】
このようにして求められる各気筒の出力燃料噴射量#nToutは、機関側制御ユニット7bに備えた各気筒毎の付着補正部23により、吸気管での燃料の壁面付着を考慮した補正が各気筒毎になされた後、エンジン1の図示しない燃料噴射装置に与えられる。そして、その付着補正がなされた出力燃料噴射量#nToutに従って、エンジン1の各気筒への燃料噴射が行われるようになっている。
【0053】
尚、上記付着補正については、本願出願人が例えば特開平8−21273号公報に詳細に開示しているので、ここではさらなる説明を省略する。また、図1において、参照符号24を付したセンサ出力選択処理部は、前記オブザーバ21による各気筒毎の実空燃比#nA/Fの推定に適したLAFセンサ5の出力KACTをエンジン1の運転状態に応じて選択するもので、これについては、本願出願人が特開平7−259588号公報にて詳細に開示しているので、ここではさらなる説明を省略する。
【0054】
一方、前記排気側制御ユニット7aは、LAFセンサ5の出力K ACTと所定の基準値FLAF/BASEとの偏差kact(=KACT−FLAF/BASE)を求める減算処理部11と、O2センサ6の出力VO2/OUTとその目標値VO2/TARGETとの偏差VO2(=VO2/OUT−VO2/TARGET)を求める減算処理部12とを備えている。
【0055】
ここで、O2センサ6の出力VO2/OUTに対する上記目標値VO2/TARGETは、触媒装置3の最適な浄化性能が得られるようなO2センサ6の出力値として排気側制御ユニット7aがエンジン1の運転状態に応じて逐次決定するものである。そして、この場合、本実施形態では、排気側制御ユニット7aは、エンジン1の回転数NE及び吸気圧PBからあらかじめ設定されたマップを用いてO2センサ6の目標値VO2/TARGETを制御サイクル毎に決定するようにしている。該目標値VO2/TARGETは、例えば図2の出力値VO2/OUT2,VO2/OUT3の間の値である。
【0056】
尚、LAFセンサ5の出力KACTに対する前記基準値FLAF/BASEは、本実施形態では、「理論空燃比」に設定されている。また、以下の説明において、前記減算処理部11,12がそれぞれ求める偏差kact,V O2をそれぞれLAFセンサ5の偏差出力kact及びO2センサ6の偏差出力VO2と称する。
【0057】
排気側制御ユニット7aはさらに、上記の偏差出力kact,VO2のデータがそれぞれLAFセンサ5の出力及びO2センサ6の出力のデータとして与えられる排気側主演算処理部13を備えている。
【0058】
この排気側主演算処理部13は、その機能として、図3に示すように、前記偏差出力kact,VO2のデータに基づいてエンジン1の目標空燃比KCMDを逐次算出する目標空燃比算出手段13aと、O2センサ6の偏差出力VO2のデータに基づいて触媒装置3の劣化状態を評価する劣化状態評価手段13bとを有する。
【0059】
この場合、より詳しくは目標空燃比算出手段13aは、排気管2のLAFセンサ5の箇所からO2センサ6の箇所にかけての触媒装置3を含む排気系(図1で参照符号Eを付した部分)を制御対象とする。そして、この対象排気系Eが有する無駄時間や、前記エンジン1及び機関側制御ユニット7bが有する無駄時間、対象排気系Eの挙動変化等を考慮しつつ、スライディングモード制御(詳しくは適応スライディングモード制御)を用いてO2センサ6の出力VO2/OUTをその目標値VO2/TARGETに収束させるようにエンジン1の目標空燃比KCMDを逐次算出するものである。
【0060】
また、劣化状態評価手段13bは、O2センサ6の偏差出力VO2の時系列データから定まる後述の劣化評価用線形関数の値に基づいて触媒装置3の劣化状態を評価すると共に、その評価結果に応じて本実施形態の装置に備えた劣化報知器29(図1参照)の動作を制御するものである。該劣化報知器29は、ランプの点灯もしくは点滅、あるいはブザーの鳴動、あるいは文字もしくは図形の表示等により触媒装置3の劣化状態を外部に報知するものである。
【0061】
尚、前記目標空燃比算出手段13aは、前記機関側制御ユニット7bと併せて本発明における空燃比操作手段に相当するものである。また、本実施形態における前記劣化状態評価手段13bは、本発明におけるパラメータ生成手段としての機能を含むものである。
【0062】
前記目標空燃比算出手段13a及び劣化状態評価手段13bをさらに説明する。
【0063】
まず、目標空燃比算出手段13aに関し、本実施形態では、該手段13aの処理を行うために、前記対象排気系Eを、前記LAFセンサ5の出力KACT(エンジン1の空燃比の検出値)から無駄時間要素及び応答遅れ要素を介してO2センサ6の出力VO2/OUT(触媒装置3を通過した排ガス中の酸素濃度の検出値)を生成する系と見なし、その挙動をあらかじめ離散時間系でモデル化している。
【0064】
さらに、前記エンジン1及び機関側制御ユニット7bから成る系を、目標空燃比KCMDから無駄時間要素を介してLAFセンサ5の出力KACTを生成する系(以下、この系を空燃比操作系と称する)と見なし、その挙動をあらかじめ離散時間系でモデル化している。
【0065】
この場合、本実施形態では、対象排気系Eの挙動を離散時間系で表現するモデル(以下、排気系モデルという)は、LAFセンサ5の出力KACT及びO2センサ6の出力VO2/OUTの代わりに、LAFセンサ5の前記偏差出力kact(=KACT−FLAF/BASE)とO2センサ6の前記偏差出力VO2(=VO2/OUT−VO2/TARGET)とを用いて、次式(1)により表現する。
【0066】
【数1】
【0067】
この式(1)は、対象排気系EがLAFセンサ5の偏差出力kactから、無駄時間要素及び応答遅れ要素を介してO2センサ6の偏差出力VO2を生成する系であると見なして、該対象排気系Eの挙動を離散時間系のモデル(より詳しくは対象排気系Eの入力量としての偏差出力kactに無駄時間を有する自己回帰モデル)で表現したものである。
【0068】
ここで、上式(1)において、「k」は前記排気側制御ユニット7aの離散時間的な制御サイクルの番数を示し(以下、同様)、「d1」は対象排気系Eに存する無駄時間(詳しくはLAFセンサ5が検出する各時点の空燃比がO2センサ6の出力VO2/OUTに反映されるようになるまでに要する無駄時間)を制御サイクル数で表したものである。この場合、対象排気系Eの無駄時間は、排気側制御ユニット7aの制御サイクルの周期(これは本実施形態では一定である)を30〜100msとしたとき、一般的には、3〜10制御サイクル分の時間(d1=3〜10)である。そして、本実施形態では、式(1)により表した排気系モデルにおける無駄時間d1の値として、対象排気系Eの実際の無駄時間と等しいか、もしくはそれよりも若干長いものにあらかじめ設定した所定の一定値(本実施形態では例えばd1=7)を設定している。
【0069】
また、式(1)の右辺第1項及び第2項はそれぞれ対象排気系Eの応答遅れ要素に対応するもので、第1項は1次目の自己回帰項、第2項は2次目の自己回帰項である。そして、「a1」、「a2」はそれぞれ1次目の自己回帰項のゲイン係数、2次目の自己回帰項のゲイン係数である。これらのゲイン係数a1,a2は別の言い方をすれば、対象排気系Eの出力量としてのO2センサ6の偏差出力VO2に係る係数である。
【0070】
さらに、式(1)の右辺第3項は対象排気系Eの入力量としてのLAFセンサ5の偏差出力kactに対象排気系Eの無駄時間d1を含めて表現したものであり、「b1」はその入力量(=LAFセンサ5の偏差出力kact)に係るゲイン係数である。これらのゲイン係数a1,a2,b1は排気系モデルの挙動を規定するパラメータであり、本実施形態では後述の同定器によって逐次同定するものである。
【0071】
一方、エンジン1及び機関側制御ユニット7bからなる前記空燃比操作系の離散時間系のモデル(以下、空燃比操作系モデルという)は、本実施形態では、排気系モデルの場合と同様にLAFセンサ5の出力KACTの代わりにLAFセンサ5の前記偏差出力kact(=KACT−FLAF/BASE)を用いると共に、これに対応させて目標空燃比KCMDの代わりに該目標空燃比KCMDの前記基準値FLAF/BASEに対する偏差kcmd(=KCMD−FLAF/BASE。これはLAFセンサ5の偏差出力kactの目標値に相当する。以下、これを目標偏差空燃比kcmdという)とを用い、次式(2)により表す。
【0072】
【数2】
【0073】
この式(2)は空燃比操作系が前記目標偏差空燃比kcmdから無駄時間要素を介してLAFセンサ5の偏差出力kactを生成する系(各制御サイクルにおける偏差出力kactが無駄時間前の目標偏差空燃比kcmdに一致するような系)であると見なして、該空燃比操作系を離散時間系のモデルで表現したものである。
【0074】
ここで、式(2)において、「d2」が空燃比操作系の無駄時間(詳しくは各時点の目標空燃比KCMDがLAFセンサ5の出力KACTに反映されるようになるまでに要する無駄時間)を排気側制御ユニット7aの制御サイクル数で表したものである。この場合、空燃比操作系の無駄時間は、エンジン1の回転数NEによって変化し、エンジン1の回転数が低くなる程、長くなる。そして、本実施形態では、式(2)により表した空燃比操作系モデルにおける無駄時間d2の値としては、上記のような空燃比操作系の無駄時間の特性を考慮し、例えばエンジン1の低速回転域の回転数であるアイドリング回転数において実際の空燃比操作系が有する無駄時間(これは、エンジン1の任意の回転数において空燃比操作系が採り得る最大側の無駄時間である)と等しいか、もしくはそれよりも若干長いものにあらかじめ設定した所定の一定値(本実施形態では例えばd2=3)を用いる。
【0075】
尚、前記空燃比操作系には、実際には、無駄時間要素の他、エンジン1に起因した応答遅れ要素も含まれる。しかるに、目標空燃比KCMDに対するLAFセンサ5の出力KACTの応答遅れは、基本的には前記機関側制御ユニット7bのフィードバック制御部14(特に適応制御器18)によって補償されるため、排気側制御ユニット7aから見た空燃比操作系では、エンジン1に起因する応答遅れ要素を考慮せずとも支障はない。
【0076】
前記排気側主演算処理部13の目標空燃比算出手段13aは、式(1)及び式(2)によりそれぞれ表現した排気系モデル及び空燃比操作系モデルに基づく目標空燃比KCMDの算出処理を排気側制御ユニット7aの制御サイクルで実行するものである。そして、この処理を行うために、図3に示すような機能的構成を具備している。
【0077】
すなわち、目標空燃比算出手段13aは、前記排気系モデル(式(1))の設定すべきパラメータである前記ゲイン係数a1,a2,b1の同定値a1ハット,a2ハット,b1ハット(以下、同定ゲイン係数a1ハット,a2ハット,b1ハットという)を制御サイクル毎に逐次求める同定器25と、対象排気系Eの無駄時間d1及び空燃比操作系の無駄時間d2を合わせた合計無駄時間d(=d1+d2)後のO2センサ6の偏差出力VO2の推定値VO2バー(以下、推定偏差出力VO2バーという)を制御サイクル毎に逐次求める推定器26と、適応スライディングモード制御の処理により前記目標空燃比KCMDを制御サイクル毎に逐次算出するスライディングモード制御器27とを具備する。
【0078】
これらの同定器25、推定器26及びスライディングモード制御器27による演算処理のアルゴリズムは以下のように構築されている。
【0079】
まず、前記同定器25は、前記式(1)により表現した排気系モデルの実際の対象排気系Eに対するモデル化誤差を極力小さくするように前記ゲイン係数a1,a2,b1の値をリアルタイムで逐次同定するものであり、その同定処理を次のように行う。
【0080】
すなわち、同定器25は、排気側制御ユニット7aの制御サイクル毎に、まず、今現在設定されている排気系モデルの同定ゲイン係数a1ハット,a2ハット,b1ハット、すなわち前回の制御サイクルで決定した同定ゲイン係数a1(k-1)ハット,a2(k-1)ハット,b1(k-1)ハットの値と、LAFセンサ5の偏差出力kact及びO2センサ6の偏差出力VO2の過去値のデータkact(k-d1-1),VO2(k-1),VO2(k-2)とを用いて、次式(3)により排気系モデル上でのO2センサ6の偏差出力VO2(排気系モデルの出力)の値VO2(k)ハット(以下、同定偏差出力VO2(k)ハットという)を求める。
【0081】
【数3】
【0082】
この式(3)は、排気系モデルを表す前記式(1)を1制御サイクル分、過去側にシフトし、ゲイン係数a1,a2,b1として同定ゲイン係数a1(k-1)ハット,a2(k-1)ハット,b1(k-1)ハットを用いたものである。また、式(3)の第3項で用いる対象排気系Eの無駄時間d1の値は、前述の如く設定した一定値(本実施形態ではd1=7)を用いる。
【0083】
ここで、次式(4),(5)で定義されるベクトルΘ及びξを導入すると(式(4),(5)中の添え字「T」は転置を意味する。以下同様。)、
【0084】
【数4】
【0085】
【数5】
【0086】
前記式(3)は、次式(6)により表される。
【0087】
【数6】
【0088】
さらに同定器25は、前記式(3)あるいは式(6)により求められるO2センサ6の同定偏差出力VO2(k)ハットと今現在のO2センサ6の偏差出力VO2(k)との偏差id/e(k)を排気系モデルの実際の対象排気系Eに対するモデル化誤差を表すものとして次式(7)により求める(以下、偏差id/eを同定誤差id/eという)。
【0089】
【数7】
【0090】
そして、同定器25は、上記同定誤差id/eを最小にするように新たな同定ゲイン係数a1(k)ハット,a2(k)ハット,b1(k)ハット、換言すれば、これらの同定ゲイン係数を要素とする新たな前記ベクトルΘ(k)(以下、このベクトルを同定ゲイン係数ベクトルΘという)を求めるものであり、その算出を、次式(8)により行う。すなわち、同定器25は、前回の制御サイクルで決定した同定ゲイン係数a1(k-1)ハット,a2(k-1)ハット,b1(k-1)ハットを、同定誤差id/e(k)に比例させた量だけ変化させることで新たな同定ゲイン係数a1(k)ハット,a2(k)ハット,b1(k)ハットを求める。
【0091】
【数8】
【0092】
ここで、式(8)中の「Kθ」は次式(9)により決定される三次のベクトル(各同定ゲイン係数a1ハット,a 2ハット,b1ハットの同定誤差id/eに応じた変化度合いを規定するゲイン係数ベクトル)である。
【0093】
【数9】
【0094】
また、上式(9)中の「P」は次式(10)の漸化式により決定される三次の正方行列である。
【0095】
【数10】
【0096】
尚、式(10)中の「λ1」、「λ2」は0<λ1≦1及び0≦λ2<2の条件を満たすように設定され、また、「P」の初期値P(0)は、その各対角成分を正の数とする対角行列である。
【0097】
この場合、式(10)中の「λ1」、「λ2」の設定の仕方によって、固定ゲイン法、漸減ゲイン法、重み付き最小二乗法、最小二乗法、固定トレース法等、各種の具体的な同定アルゴリズムが構成され、本実施形態では、例えば最小二乗法(この場合、λ1=λ2=1)を採用している。
【0098】
本実施形態における同定器25は基本的には前述のようなアルゴリズム(演算処理)によって、前記同定誤差id/eを最小化するように排気系モデルの前記同定ゲイン係数a1ハット,a2ハット,b1ハットを制御サイクル毎に逐次求めるものである。このような処理によって、実際の対象排気系Eに適合した同定ゲイン係数a1ハット,a2ハット,b1ハットが逐次得られる。
【0099】
以上説明したアルゴリズムが同定器25が実行する基本的なアルゴリズムである。
【0100】
次に、前記推定器26は、後に詳細を説明するスライディングモード制御器27による目標空燃比KCMDの算出処理に際しての対象排気系Eの無駄時間d1及び前記空燃比操作系の無駄時間d2の影響を補償するために、前記合計無駄時間d(=d1+d2)後のO2センサ6の偏差出力VO2の推定値である前記推定偏差出力VO2バーを制御サイクル毎に逐次求めるものである。その推定処理のアルゴリズムは次のように構築されている。
【0101】
まず、排気系モデルを表す前記式(1)に、空燃比操作系モデルを表す式(2)を適用すると、次式(11)が得られる。
【0102】
【数11】
【0103】
この式(11)は、対象排気系E及び空燃比操作系を合わせた系を、目標偏差空燃比kcmdから対象排気系E及び空燃比操作系の両者の無駄時間要素と対象排気系Eの応答遅れ要素とを介してO2センサ6の偏差出力VO2を生成する系であるとして、該系の挙動を離散時間系のモデルで表現したものである。
【0104】
そして、この式(11)を用いることで、各制御サイクルにおける前記合計無駄時間d後のO2センサ6の偏差出力VO2(k+d)の推定値である前記推定偏差出力VO2 (k+d)バーは、O2センサ6の偏差出力VO2の現在値及び過去値の時系列データVO2(k)及びVO2(k-1)と、スライディングモード制御器27が求める目標空燃比KCMD(具体的な求め方は後述する)に相当する目標偏差空燃比kcmd(=KCMD−FLAF/BASE)の過去値の時系列データkcmd(k-j)(j=1,2,…,d)とを用いて次式(12)により表される。
【0105】
【数12】
【0106】
ここで、式(12)において、α1,α2は、それぞれ同式(12)中のただし書きで定義した行列Aの巾乗Ad(d:合計無駄時間)の第1行第1列成分、第1行第2列成分である。また、βj(j=1,2,…,d)は、それぞれ行列Aの巾乗Aj-1(j=1,2,…,d)と同式(12)中のただし書きで定義したベクトルBとの積Aj-1・Bの第1行成分である。
【0107】
さらに、式(12)中の目標偏差空燃比kcmdの過去値の時系列データkcmd(k-j)(j=1,2,…,d)のうち、現在から空燃比操作系の無駄時間d2以前の目標偏差空燃比kcmdの過去値の時系列データkcmd(k-d2),kcmd(k-d2-1),…,kcmd(k-d)は前記式(2)によって、それぞれ、LAFセンサ5の偏差出力kactの現在以前に得られるデータkact(k),kact(k-1),…,kact(k-d+d2)に置き換えることができる。そして、この置き換えを行うことで、次式(13)が得られる。
【0108】
【数13】
【0109】
この式(13)が本実施形態において、推定器26が前記推定偏差出力VO2(k+d)バーを算出するための基本式である。つまり、本実施形態では、推定器26は、制御サイクル毎に、O2センサ6の偏差出力VO2の時系列データVO2(k)及びVO2(k-1)と、スライディングモード制御器27が過去に求めた目標空燃比KCMDを表す目標偏差空燃比kcmdの過去値の時系列データkcmd(k-j)(j=1,…,d2-1)と、LAFセンサ5の偏差出力kactの現在値及び過去値の時系列データkact(k-i)(i=0,…,d1)とを用いて式(13)の演算を行うことによって、O2センサ6の推定偏差出力VO2(k+d)バーを求める。
【0110】
この場合、本実施形態では、式(13)により推定偏差出力VO2(k+d)バーを算出するために必要となる係数α1,α2及びβj(j=1,2,…,d)の値は、基本的には、前記ゲイン係数a1,a2,b1(これらは式(12)のただし書きで定義した行列A及びベクトルBの成分である)の同定値である前記同定ゲイン係数a1ハット,a2ハット,b1ハットを用いて算出する。また、式(13)の演算で必要となる無駄時間d1,d2の値は、前述の如く設定した値を用いる。
【0111】
尚、推定偏差出力VO2(k+d)バーは、LAFセンサ5の偏差出力kactのデータを使用せずに、式(12)の演算により求めるようにしてもよいが、推定偏差出力VO2(k+d)バーの信頼性を高める上では、エンジン1等の実際の挙動が反映されるLAFセンサ5の偏差出力kactのデータを用いた式(13)の演算により推定偏差出力VO2(k+d)バーを求めることが好ましい。また、空燃比操作系の無駄時間d2を「1」に設定できるような場合には、式(12)中の目標偏差空燃比kcmdの過去値の時系列データkcmd(k-j)(j=1,2,…,d)の全てをそれぞれ、LAFセンサ5の偏差出力kactの現在以前に得られる時系列データkact(k),kact(k-1),…,kact(k-d+d2)に置き換えることができる。このため、この場合には、推定偏差出力VO2(k+d)バーは、目標偏差空燃比kcmdのデータを含まない次式(14)により求めることができる。
【0112】
【数14】
【0113】
次に、前記スライディングモード制御器27を説明する。
【0114】
本実施形態のスライディングモード制御器27は、通常的なスライディングモード制御に、外乱等の影響を極力排除するための適応則(適応アルゴリズム)を加味した適応スライディングモード制御によりO2センサ6の出力VO2/OUTをその目標値VO2/TARGETに整定させるように(O2センサ6の偏差出力VO2を「0」に収束させるように)、制御対象である前記対象排気系Eに与えるべき入力量(詳しくは、LAFセンサ5の出力KACT(空燃比の検出値)と前記基準値FLAF/BASEとの偏差の目標値で、これは前記目標偏差空燃比kcmdに等しい。以下、この入力量をSLD操作入力Uslと称する)を決定し、その決定したSLD操作入力Uslから前記目標空燃比KCMDを決定するものである。そして、その処理のためのアルゴリズムは次のように構築されている。
【0115】
まず、スライディングモード制御器27が実行する適応スライディングモード制御のアルゴリズムに必要な切換関数と、この切換関数により定義される超平面(これはすべり面とも言われる)とについて説明する。
【0116】
本実施形態におけるスライディングモード制御の基本的な考え方としては、制御すべき状態量(制御量)として、例えば各制御サイクルで得られたO2センサ6の偏差出力VO2(k)と、その1制御サイクル前に得られた偏差出力VO2(k-1)とを用い、スライディングモード制御用の切換関数σを、次式(15)のように、これらの偏差出力VO2(k),VO2(k-1)を変数成分とする線形関数として定義する。尚、前記偏差出力VO2(k),VO2(k-1)を成分とするベクトルとして式(15)中の但し書きで定義したベクトルXを以下、状態量Xという。
【0117】
【数15】
【0118】
この場合、切換関数σの係数s1,s2は、次式(16)の条件を満たすように設定する。
【0119】
【数16】
【0120】
尚、本実施形態では、簡略化のために係数s1=1とし(この場合、s2/s1=s2である)、−1<s2<1の条件を満たすように係数s2の値を設定している。
【0121】
このような切換関数σに対して、スライディングモード制御用の超平面はσ=0なる式によって定義されるものである。この場合、状態量Xは二次系であるので超平面σ=0は図4に示すように直線となる。該超平面は、位相空間の次数によって、切換線又は切換面とも言われる。
【0122】
尚、本実施形態では、スライディングモード制御用の切換関数の変数成分である状態量として、実際には前記推定器26により求められる前記推定偏差出力VO2バーの時系列データを用いるのであるが、これについては後述する。
【0123】
本実施形態で用いる適応スライディングモード制御は、状態量X=(VO2(k),VO2(k-1))を上記の如く設定した超平面σ=0に収束させるための制御則である到達則と、その超平面σ=0への収束に際して外乱等の影響を補償するための制御則である適応則(適応アルゴリズム)とにより該状態量Xを超平面σ=0に収束させる(図4のモード1)。そして、該状態量Xを所謂、等価制御入力によって超平面σ=0に拘束しつつ、該状態量Xを超平面σ=0上の平衡点であるVO2(k)=VO2(k-1)=0となる点、すなわち、O2センサ6の出力VO2/OUTの時系列データVO2/OUT(k),VO2/OUT(k-1)が目標値VO2/TARGETに一致するような点に収束させる(図4のモード2)。
【0124】
上記のように状態量Xを超平面σ=0の平衡点に収束させるためにスライディングモード制御器27が生成する前記SLD操作入力Usl(=目標偏差空燃比kcmd)は、状態量Xを超平面σ=0上に拘束するための制御則に従って対象排気系Eに与えるべき入力成分である等価制御入力Ueqと、前記到達則に従って対象排気系Eに与えるべき入力成分Urch(以下、到達則入力Urchという)と、前記適応則に従って対象排気系Eに与えるべき入力成分Uadp(以下、適応則入力Uadpという)との総和により表される(次式(17))。
【0125】
【数17】
【0126】
そして、これらの等価制御入力Ueq、到達則入力Urch及び適応則入力Uad pは、本実施形態では、前記式(11)により表される離散時間系のモデル(式(1)中のLAFセンサ5の偏差出力kact(k-d1)を合計無駄時間dを用いた目標偏差空燃比kcmd(k-d)で置き換えたモデル)に基づいて、次のように決定する。
【0127】
まず、状態量Xを超平面σ=0に拘束するために対象排気系Eに与えるべき入力成分である前記等価制御入力Ueqは、σ(k+1)=σ(k)=0なる条件を満たす目標偏差空燃比kcmdである。そして、このような条件を満たす等価制御入力Ueqは、式(11)と式(15)とを用いて次式(18)により与えられる。
【0128】
【数18】
【0129】
この式(18)が本実施形態において、制御サイクル毎に等価制御入力Ueq(k)を求めるための基本式である。
【0130】
次に、前記到達則入力Urchは、本実施形態では、基本的には次式(19)により決定するものとする。
【0131】
【数19】
【0132】
すなわち、到達則入力Urchは、前記合計無駄時間dを考慮し、合計無駄時間d後の切換関数σの値σ(k+d)に比例させるように決定する。
【0133】
この場合、式(19)中の係数F(これは到達則のゲインを規定する)は、次式(20)の条件を満たすように設定する。
【0134】
【数20】
【0135】
尚、切換関数σの値の挙動に関しては、該切換関数σの値が超平面σ=0に対して振動的な変化(所謂チャタリング)を生じる虞れがあり、このチャタリングを抑制するためには、到達則入力Urchに係わる係数Fは、さらに次式(21)の条件を満たすように設定することが好ましい。
【0136】
【数21】
【0137】
次に、前記適応則入力Uadpは、本実施形態では、基本的には次式(22)により決定するものとする(式(22)中のΔTは排気側制御ユニット7aの制御サイクルの周期である)。
【0138】
【数22】
【0139】
すなわち、適応則入力Uadpは、合計無駄時間dを考慮し、該合計無駄時間d後までの切換関数σの値と排気側制御ユニット7aの周期ΔTとの積の制御サイクル毎の積算値(これは切換関数σの値の積分値に相当する)に比例させるように決定する。
【0140】
この場合、式(22)中の係数G(これは適応則のゲインを規定する)は、次式(23)の条件を満たすように設定する。
【0141】
【数23】
【0142】
尚、前記式(16)、(20)、(21)、(23)の設定条件のより具体的な導出の仕方については、本願出願人が既に特開平11−93741号公報等にて詳細に説明しているので、ここでは詳細な説明を省略する。
【0143】
本実施形態におけるスライディングモード制御器27は、基本的には前記式(18)、(19)、(22)により決定される等価制御入力Ueq、到達則入力Urch及び適応則入力Uadpの総和(Ueq+Urch+Uadp)を対象排気系Eに与えるべきSLD操作入力Uslとして決定するのであるが、前記式(18)、(19)、(22)で使用するO2センサ6の偏差出力VO2(k+d),VO2(k+d-1)や、切換関数σの値σ(k+d)等は未来値であるので直接的には得られない。
【0144】
そこで、本実施形態では、スライディングモード制御器27は、実際には、前記式(18)により前記等価制御入力Ueqを決定するためのO2センサ6の偏差出力VO2(k+d),VO2(k+d-1)の代わりに、前記推定器26で求められる推定偏差出力VO2(k+d)バー,VO2(k+d-1)バーを用い、次式(24)により制御サイクル毎の等価制御入力Ueqを算出する。
【0145】
【数24】
【0146】
また、本実施形態では、実際には、推定器26により前述の如く逐次求められた推定偏差出力VO2バーの時系列データを制御すべき状態量とし、前記式(15)により設定された切換関数σに代えて、次式(25)によりスライディングモード制御用の切換関数σバーを定義する(この切換関数σバーは、前記式(15)の偏差出力VO2の時系列データを推定偏差出力VO2バーの時系列データで置き換えたものに相当する)。
【0147】
【数25】
【0148】
そして、スライディングモード制御器27は、前記式(19)により前記到達則入力Urchを決定するための切換関数σの値の代わりに、前記式(25)により表される切換関数σバーの値を用いて次式(26)により制御サイクル毎の到達則入力Urchを算出する。
【0149】
【数26】
【0150】
同様に、スライディングモード制御器27は、前記式(22)により前記適応則入力Uadpを決定するための切換関数σの値の代わりに、前記式(25)により表される切換関数σバーの値を用いて次式(27)により制御サイクル毎の適応則入力Uadpを算出する。
【0151】
【数27】
【0152】
尚、前記式(24),(26),(27)により等価制御入力Ueq、到達則入力Urch及び適応則入力Uadpを算出する際に必要となる前記ゲイン係数a1,a2,b1としては、本実施形態では基本的には前記同定器25により求められた最新の同定ゲイン係数a1(k)ハット,a2(k)ハット,b1(k)ハットを用いる。
【0153】
そして、スライディングモード制御器27は、前記式(24),(26),(27)によりそれぞれ求められる等価制御入力Ueq、到達則入力Urch及び適応則入力Uadpの総和を対象排気系Eに与えるべき前記SLD操作入力Uslとして求める(前記式(17)を参照)。尚、この場合において、前記式(24)、(26)、(27)中で用いる前記係数s1,s2,F,Gの設定条件は前述の通りである。
【0154】
これが、本実施形態において、スライディングモード制御器27により、対象排気系Eに与えるべきSLD操作入力Usl(=目標偏差空燃比kcmd)を制御サイクル毎に決定するための基本的な演算処理(アルゴリズム)である。このようにしてSLD操作入力Uslを決定することで、該SLD操作入力Uslは、O2センサ6の推定偏差出力VO2バーを「0」に収束させるように(結果的にはO2センサ6の出力VO2/OUTを目標値VO2/TARGETに収束させるように)決定される。
【0155】
ところで、本実施形態におけるスライディングモード制御器27は最終的には前記目標空燃比KCMDを制御サイクル毎に逐次求めるものあるが、前述のように求められるSLD操作入力U slは、LAFセンサ5で検出される空燃比と前記基準値FLAF/BASEとの偏差の目標値、すなわち前記目標偏差空燃比kcmdである。このため、スライディングモード制御器27は、最終的には、次式(28)に示すように、制御サイクル毎に、前述の如く求めたSLD操作入力Uslに前記基準値FLAF/BASEを加算することで、目標空燃比KCMDを決定する。
【0156】
【数28】
【0157】
以上が本実施形態でスライディングモード制御器27により目標空燃比KCMDを決定するための基本的アルゴリズムである。
【0158】
尚、本実施形態では、スライディングモード制御器27による適応スライディングモード制御の処理の安定性を判別して、前記SLD操作入力Uslの値を制限したりするのであるが、これについては後述する。
【0159】
次に、前記劣化状態評価手段13bの処理を説明する。
【0160】
前述のようにスライディングモード制御器27により目標空燃比KCMDを逐次求め、この目標空燃比KCMDにLAFセンサ5の出力KACT(エンジン1の空燃比の検出値)を収束させるように前記機関側制御ユニット7bによりエンジン1の燃料噴射量を調整したとき、O2センサ6の出力VO2/OUTの時系列データは、例えば前記超平面σ=0に対して、触媒装置3の劣化状態に応じた特徴的な変化を呈する。
【0161】
この様子を図5〜図7を参照して説明する。これらの図5〜図7は、それぞれ新品の触媒装置3、比較的小さな劣化度合いの触媒装置3、比較的大きな劣化度合いの触媒装置3について上記のようにエンジン1の燃料噴射量を目標空燃比KCMDに応じて調整した場合に、排気側制御ユニット7aの制御サイクル毎に得られるO2センサ6の偏差出力VO2の時系列データVO2(k),VO2(k-1)の組、すなわち前記状態量Xのサンプリングデータを点描で示したものである。また、これらの図5〜図7に併記した直線は、前記超平面σ=0である。尚、これらの図5〜図7では、O2センサ6の出力VO2/OUTの目標値VO2/TARGETは一定値としている。
【0162】
これらの図5〜図7のうち、図5に見られるように、状態量Xは、触媒装置3の新品状態では、超平面σ=0の近傍に集積する傾向がある。そして、触媒装置3の劣化が進行すると、図6や図7に見られるように、状態量Xは、超平面σ=0から離間する側に分散する(超平面σ=0のまわりのばらつきが大きくなる)傾向がある。しかも、その分散化は、触媒装置3の劣化度合いが大きい程(劣化が進行している程)、高まる傾向がある。換言すれば、前記式(15)により定めた切換関数σの値が、触媒装置3の劣化の進行に伴い「0」から離間した値を採りやすくなり、切換関数σの値の「0」に対するばらつきが大きくなる。これは、触媒装置3の劣化の進行に伴い、前記式(1)の排気系モデルの誤差が生じやすくなるため、状態量Xの超平面σ=0への収束性が低下するためであると考えられる。
【0163】
尚、上記のような傾向は、O2センサ6の出力VO2/OUTの目標値VO2/TARGETは一定値とした場合、推定器26により求められる推定偏差出力VO2バーを変数成分として前記式(25)により定めた切換関数σバー(本実施形態でスライディングモード制御用の切換関数として実際に用いる切換関数)の値についても同様に見られる傾向である。但し、切換関数σバーは、O2センサ6の偏差出力VO2の推定値を用いるものであるのに対し、式(15)の切換関数はO2センサ6の実際の偏差出力VO2を用いるものであるため、触媒装置3の実際の劣化状態がよりよく反映されるものと考えられる。
【0164】
このようなことから、本実施形態では、前記式(15)により定めた切換関数σの値を用いて、触媒装置3の劣化状態を評価するための劣化評価用パラメータを後述のように生成し、その劣化評価用パラメータの値に基づいて触媒装置3の劣化状態を評価する。尚、本実施形態では、スライディングモード制御用の実際の切換関数は、前述の如く、推定器26が求める推定偏差出力VO2バーを変数成分として前記式(25)により定義した切換関数σバーである。つまり、厳密に言えば、式(15)により定めた切換関数σは、本実施形態でのスライディングモード制御用の切換関数ではない。そこで、以下の説明では、式(15)により定めた関数σを、劣化評価用線形関数σという。
【0165】
前記劣化状態評価手段13bが上記劣化評価用線形関数σを用いて触媒装置3の劣化状態の評価を行うための基本的アルゴリズムは、本実施形態では次のように構築されている。
【0166】
すなわち、触媒装置3の劣化の進行に伴う前述のような劣化評価用線形関数σの値の変化の傾向を考慮し、劣化状態評価手段13bは、該劣化評価用線形関数σの値の二乗値σ2を制御サイクル毎に逐次求める。
【0167】
そして、この二乗値σ2にローパス特性のフィルタリング処理を施すことにより該二乗値σ2の中心値(以下、これに参照符号LSσ2を付する)を基本的な劣化評価用パラメータとして求める。
【0168】
この場合、劣化評価用パラメータLSσ2を求めるための上記フィルタリング処理は、逐次型の統計処理アルゴリズムにより構成され、次式(29)により与えられる。
【0169】
【数29】
【0170】
すなわち、劣化評価用パラメータLSσ2は、排気側制御ユニット7aの制御サイクル毎に、該劣化評価用パラメータLSσ2の前回値LSσ2(k-1)と、前記二乗値σ2の現在値σ2(k)と、次式(30)の漸化式により制御サイクル毎に更新されるゲインパラメータBPとから逐次更新されつつ求められる。
【0171】
【数30】
【0172】
ここで、式(30)中の「η1」、「η2」は0<η1≦1及び0≦η2<2の条件を満たすように設定され、その設定の仕方によって、固定ゲイン法、漸減ゲイン法、重み付き最小二乗法、最小二乗法、固定トレース法等、各種の具体的なアルゴリズムが構成される。本実施形態では、例えばη1を「1」よりも小さな正の所定値に設定する(0<η1<1)と共に、η2=1とし、重み付最小二乗法のアルゴリズムを採用している。
【0173】
このようにして、劣化評価用線形関数σの二乗値σ2の中心値(本実施形態では最小二乗中心値)としての劣化評価用パラメータLSσ2を求めたとき、該劣化評価用パラメータLSσ2の値は、触媒装置3の劣化状態に対して、図8に示すような傾向を呈する(但し、ここではO2センサ6の出力VO2/OUTの目標値VO2/TARGETは一定であるとする)。図8は前記図5〜図7の触媒装置3の各劣化状態において、前述のように求められる劣化評価用パラメータLSσ2と触媒装置3を流れる排ガスの流量(以下、排ガスボリュームという)との関係を示すものである。
【0174】
図8に見られるように、劣化評価用パラメータLSσ2は、触媒装置3の各劣化状態において前記排ガスボリュームによらずにほぼ一定となると共に、触媒装置3の劣化の進行に伴い、値が大きくなる。従って、劣化評価用パラメータLSσ2の値が触媒装置3の劣化度合いを表すものとなる。
【0175】
また、本実施形態では、例えば触媒装置3の交換を要するかもしくはその交換時期が近い程度に該触媒装置3が劣化した状態(以下、劣化進行状態という)と、それほどには劣化していない状態(以下、未劣化状態という)とに分けて触媒装置3の劣化状態を把握することとし、上記「劣化進行状態」において、その旨を前記劣化報知器29により報知する。
【0176】
そこで、本実施形態では、図8の破線で示すように、劣化評価用パラメータLSσ2に対する閾値CATAGELMTをあらかじめ設定しておく。そして、本実施形態の基本的な考え方では、劣化評価用パラメータLSσ2が上記閾値CATAGELMT以上となったときに、触媒装置3の劣化状態が「劣化進行状態」であると判断し、劣化評価用パラメータLSσ2が上記閾値CATAGELMTに満たない場合は、触媒装置3の劣化状態が「未劣化状態」であると判断する。
【0177】
ところで、上述した劣化状態評価手段13bのアルゴリズムの説明では、O2センサ6の出力VO2/OUTの目標値VO2/TARGETは一定であるとして説明したが、本実施形態では、該目標値VO2/TARGETは、エンジン1の運転状態(具体的には回転数NE及び吸気圧PB)に応じて逐次可変的に決定される。そして、この場合、O2センサ6の出力特性は非線形であるため(図2参照)、O2センサ6の出力VO2/OUTを目標値VO2/TARGETに収束させる制御を行なっている状態における前記劣化評価用線形関数σの値の「0」に対するばらつき度合いは、触媒装置3の劣化状態だけでなく、目標値VO2/TARGETの影響も受ける。例えば、該目標値VO2/TARGETが、図2のVO2/OUT1である状態とVO2/OUT2である状態とを比較すると、前者の状態での目標値VO2/TARGET=VO2/OUT1の近傍におけるO2センサ6の出力の感度が高いため、O2センサ6の出力が目標値VO2/TARGETに対して変動を生じ易い。この結果、VO2/TARGET=VO2/OUT1である状態と、VO2/TARGET=VO2/OUT2である状態とでは、触媒装置3の劣化状態が同じであっても、前者の状態の方が後者の状態よりも、前記劣化評価用線形関数σの値のばらつき度合いが大きくなり易く、ひいては劣化評価用パラメータLSσ2の値が大きくなり易い。
【0178】
そこで、本実施形態では、劣化状態評価手段13bは、前記式(29)により劣化評価用パラメータLSσ2を算出するのと並行して、O2センサ6の出力VO2/OUTの平均値VO2/OUTAVEを逐次更新しつつ算出する。そして、劣化状態評価手段13bは、触媒装置3の劣化状態を判断するに際しては、実際には、前記式(29)により最終的に得られる劣化評価用パラメータLSσ2を、O2センサ6の出力VO2/OUTの平均値VO2/OUTAVE(最新値)に応じて補正し、その補正後の劣化評価用パラメータを前記閾値CATAGELMTと比較することで、触媒装置3の劣化状態が「未劣化状態」であるか「劣化進行状態」であるかを判断する(以下の説明では、補正前の劣化評価用パラメータLSσ2を基本劣化評価用パラメータLSσ2と称し、補正後の劣化評価パラメータを補正後劣化評価用パラメータCLSσ2と称する)。
【0179】
この場合、前記基本劣化評価用パラメータLSσ2の補正は例えば次のように行なわれる。すなわち、劣化状態評価手段13bは、図9に示すようにあらかじめ定められたデータテーブルに基づいて基本劣化評価用パラメータLSσ2を補正する)ための補正係数KRESS(>0)を求め、この求めた補正係数KRESSを基本劣化評価用パラメータLSσ2に乗算することにより補正後劣化評価用パラメータCLSσ2を求める。ここで、図9のデータテーブルでは、補正係数KRESSの値は、O2センサ6の出力VO2/OUTの平均値VO2/OUTAVEがO2センサ6の出力VO2/OUTの感度がほぼ最大となる値、すなわち図2のVO2/OUT1であるとき(これは、目標値VO2/TARGETの平均値がVO2/OUT1の近傍の値であることと同等である)、最小値となり、平均値VO2/OUTAVEのVO2/OUT1からの差|VO2/OUTAVE−VO2/OUT1|が大きくなるに伴って補正係数KRESSの値が大きくなるように定められている。
【0180】
このような補正係数KRESSを基本劣化評価用パラメータLSσ2を乗算してなる補正後劣化評価用パラメータCLSσ2は、基本的には、O2センサ6の出力VO2/OUTの目標値VO2/TARGETによらずに触媒装置3の劣化状態のみに応じて変化するものとなる。
【0181】
以上説明したアルゴリズムが、劣化状態評価手段13bによる触媒装置3の劣化状態の評価の基本的なアルゴリズムである。尚、劣化状態評価手段13bは、触媒装置3の劣化状態の評価に際して排ガスボリュームの変化状態の把握等、付加的な処理も行うのであるが、これについては後述する。
【0182】
次に、前記機関側制御ユニット7bの大局的フィードバック制御部15、特に前記適応制御器18をさらに説明する。
【0183】
前記図1を参照して、大局的フィードバック制御部15は、前述のようにLAFセンサ5の出力KACT(空燃比の検出値)を目標空燃比KCMDに収束させるようにフィードバック制御を行うものであるが、このとき、このようなフィードバック制御を周知のPID制御だけで行うようにすると、エンジン1の運転状態の変化や経年的特性変化等、動的な挙動変化に対して、安定した制御性を確保することが困難である。
【0184】
前記適応制御器18は、上記のようなエンジン1の動的な挙動変化を補償したフィードバック制御を可能とする漸化式形式の制御器であり、I.D.ランダウ等により提唱されているパラメータ調整則を用いて、図10に示すように、複数の適応パラメータを設定するパラメータ調整部30と、設定された適応パラメータを用いて前記フィードバック操作量KSTRを算出する操作量算出部31とにより構成されている。
【0185】
ここで、パラメータ調整部30について説明すると、ランダウ等の調整則では、離散系の制御対象の伝達関数B(Z-1)/A(Z-1)の分母分子の多項式を一般的に下記の式(31),(32)のようにおいたとき、パラメータ調整部30が設定する適応パラメータθハット(j)(jは制御サイクルの番数を示す)は、式(33)のようにベクトル(転置ベクトル)で表される。また、パラメータ調整部30への入力ζ(j)は、式(34)のように表される。この場合、本実施形態では、大局的フィードバック制御部15の制御対象であるエンジン1が一次系で3制御サイクル分の無駄時間dp(エンジン1の燃焼サイクルの3サイクル分の時間)を持つプラントと考え、式(31)〜式(34)でm=n=1,dp=3とし、設定する適応パラメータはs0,r1,r2,r3,b0の5個とした(図10参照)。尚、式(34)の上段式及び中段式におけるus,ysは、それぞれ、制御対象への入力(操作量)及び制御対象の出力(制御量)を一般的に表したものであるが、本実施形態では、上記入力をフィードバック操作量KSTR、制御対象(エンジン1)の出力を前記LAFセンサ5の出力KACT(空燃比の検出値)とし、パラメータ調整部30への入力ζ(j)を、式(34)の下段式により表す(図10参照)。
【0186】
【数31】
【0187】
【数32】
【0188】
【数33】
【0189】
【数34】
【0190】
ここで、前記式(33)に示される適応パラメータθハットは、適応制御器18のゲインを決定するスカラ量要素b0ハット(j)、操作量を用いて表現される制御要素BRハット(Z-1,j)、及び制御量を用いて表現される制御要素Sハット(Z- 1,j)からなり、それぞれ、次式(35)〜(37)により表現される(図10の操作量算出部31のブロック図を参照)。
【0191】
【数35】
【0192】
【数36】
【0193】
【数37】
【0194】
パラメータ調整部30は、これらのスカラ量要素や制御要素の各係数を設定して、それを式(33)に示す適応パラメータθハットとして操作量算出部31に与えるもので、現在から過去に渡るフィードバック操作量KSTRの時系列データとLAFセンサ5の出力KACTとを用いて、該出力KACTが前記目標空燃比KCMDに一致するように、適応パラメータθハットを算出する。
【0195】
この場合、具体的には、適応パラメータθハットは、次式(38)により算出する。
【0196】
【数38】
【0197】
同式(38)において、Γ(j)は、適応パラメータθハットの設定速度を決定するゲイン行列(この行列の次数はm+n+dp)、eアスタリスク(j)は、適応パラメータθハットの推定誤差を示すもので、それぞれ式(39),(40)のような漸化式で表される。
【0198】
【数39】
【0199】
【数40】
【0200】
ここで、式(40)中の「D(Z-1)」は、収束性を調整するための、漸近安定な多項式であり、本実施形態ではD(Z-1)=1としている。
【0201】
尚、式(39)のλ1(j),λ2(j)の選び方により、漸減ゲインアルゴリズム、可変ゲインアルゴリズム、固定トレースアルゴリズム、固定ゲインアルゴリズム等の種々の具体的なアルゴリズムが得られる。エンジン1の燃料噴射あるいは空燃比等の時変プラントでは、漸減ゲインアルゴリズム、可変ゲインアルゴリズム、固定ゲインアルゴリズム、および固定トレースアルゴリズムのいずれもが適している。
【0202】
前述のようにパラメータ調整部30により設定される適応パラメータθハット(s0,r1,r2,r3,b0)と、前記排気側主演算処理部13の目標空燃比算出手段13aにより決定される目標空燃比KCMDとを用いて、操作量算出部31は、次式(41)の漸化式により、フィードバック操作量KSTRを求める。図10の操作量算出部31は、同式(41)の演算をブロック図で表したものである。
【0203】
【数41】
【0204】
尚、式(41)により求められるフィードバック操作量KSTRは、LAFセンサ5の出力KACTが目標空燃比KCMDに一致する状態において、「目標空燃比KCMD」となる。このために、前述の如く、フィードバック操作量KSTRを除算処理部19によって目標空燃比KCMDで除算することで、前記フィードバック補正係数KFBとして使用できるフィードバック操作量kstrを求めるようにしている。
【0205】
このように構築された適応制御器18は、前述したことから明らかなように、制御対象であるエンジン1の動的な挙動変化を考慮した漸化式形式の制御器であり、換言すれば、エンジン1の動的な挙動変化を補償するために、漸化式形式で記述された制御器である。そして、より詳しくは、漸化式形式の適応パラメータ調整機構を備えた制御器と定義することができる。
【0206】
尚、この種の漸化式形式の制御器は、所謂、最適レギュレータを用いて構築する場合もあるが、この場合には、一般にはパラメータ調整機構は備えられておらず、エンジン1の動的な挙動変化を補償する上では、前述のように構成された適応制御器18が好適である。
【0207】
以上が、本実施形態で採用した適応制御器18の詳細である。
【0208】
尚、適応制御器18と共に、大局的フィードバック制御部15に具備したPID制御器17は、一般のPID制御と同様に、LAFセンサ5の出力KACTと、その目標空燃比KCMDとの偏差から、比例項(P項)、積分項(I項)及び微分項(D項)を算出し、それらの各項の総和をフィードバック操作量KLA Fとして算出する。この場合、本実施形態では、積分項(I項)の初期値を“1”とすることで、LAFセンサ5の出力KACTが目標空燃比KCMDに一致する状態において、フィードバック操作量KLAFが“1”になるようにし、該フィードバック操作量KLAFをそのまま燃料噴射量を補正するための前記フィードバック補正係数KFBとして使用することができるようしている。また、比例項、積分項及び微分項のゲインは、エンジン1の回転数NEと吸気圧PBとから、あらかじめ定められたマップを用いて決定される。
【0209】
また、大局的フィードバック制御部15の前記切換部20は、エンジン1の冷却水温の低温時や、高速回転運転時、吸気圧の低圧時等、エンジン1の燃焼が不安定なものとなりやすい場合、あるいは、目標空燃比KCMDの変化が大きい時や、空燃比のフィードバック制御の開始直後等、これに応じたLAFセンサ5の出力KACTが、そのLAFセンサ5の応答遅れ等によって、信頼性に欠ける場合、あるいは、エンジン1のアイドル運転時のようエンジン1の運転状態が極めて安定していて、適応制御器18による高ゲイン制御を必要としない場合には、PID制御器17により求められるフィードバック操作量KLAFを燃料噴射量を補正するためのフィードバック補正量数KFBとして出力する。そして、上記のような場合以外の状態で、適応制御器18により求められるフィードバック操作量KSTRを目標空燃比KCMDで除算してなるフィードバック操作量kstrを燃料噴射量を補正するためのフィードバック補正係数KFBとして出力する。これは、適応制御器18が、高ゲイン制御で、LAFセンサ5の出力KACTを急速に目標空燃比KCMDに収束させるように機能するため、上記のようにエンジン1の燃焼が不安定となったり、LAFセンサ5の出力KACTの信頼性に欠ける等の場合に、適応制御器18のフィードバック操作量KSTRを用いると、かえって空燃比の制御が不安定なものとなる虞れがあるからである。
【0210】
このような切換部20の作動は、例えば特開平8−105345号公報に本願出願人が詳細に開示しているので、ここでは、さらなる説明を省略する。
【0211】
次に本実施形態の装置の作動の詳細を説明する。
【0212】
まず、図11のフローチャートを参照して、前記機関側制御ユニット7bによるエンジン1の空燃比の制御のためのエンジン1の各気筒毎の出力燃料噴射量#nTout(n=1,2,3,4)の算出処理について説明する。機関側制御ユニット7bは、各気筒毎の出力燃料噴射量#nToutの算出処理をエンジン1のクランク角周期(TDC)と同期した制御サイクルで次のように行う。
【0213】
機関側制御ユニット7bは、まず、前記LAFセンサ5及びO2センサ6を含む各種センサの出力を読み込む(STEPa)。この場合、LAFセンサ5の出力KACT及びO2センサ6の出力VO2/OUTはそれぞれ過去に得られたものを含めて時系列的に図示しないメモリに記憶保持される。
【0214】
次いで、基本燃料噴射量算出部8によって、前述の如くエンジン1の回転数NE及び吸気圧PBに対応する燃料噴射量をスロットル弁の有効開口面積に応じて補正してなる基本燃料噴射量Timが求められる(STEPb)。さらに、第1補正係数算出部9によって、エンジン1の冷却水温やキャニスタのパージ量等に応じた第1補正係数KTOTALが算出される(STEPc)。
【0215】
次いで、機関側制御ユニット7bは、エンジン1の運転モードが排気側主演算処理部13の目標空燃比算出手段13aが生成する目標空燃比KCMDを使用して燃料噴射量の調整を行う運転モード(以下、通常運転モード)であるか否かの判別処理を行って、該運転モードが通常運転モードであるか否かをそれぞれ値「1」,「0」で表すフラグf/prism/onの値を設定する(STEPd)。
【0216】
上記の判別処理では、図12に示すように、O2センサ6及びLAFセンサ5が活性化しているか否かの判別が行われる(STEPd−1,d−2)。このとき、いずれかが活性化していない場合には、排気側主演算処理部13の処理に使用するO2センサ6やLAFセンサ5の検出データを精度よく得ることができない。従って、エンジン1の運転モードは通常運転モードではなく、フラグf /prism/onの値を「0」にセットする(STEPd−10)。
【0217】
また、エンジン1のリーン運転中(希薄燃焼運転中)であるか否か(STEPd−3)、エンジン1の始動直後の触媒装置3の早期活性化を図るためにエンジン1の点火時期が遅角側に制御されているか否か(STEPd−4)、エンジン1のスロットル弁が略全開であるか否か(STEPd−5)、及びエンジン1への燃料供給の停止中(フュエルカット中)であるか否か(STEPd−6)の判別が行われる。これらのいずれかの条件が成立している場合には、排気側主演算処理部13が生成する目標空燃比KCMDを使用してエンジン1の燃料供給を制御することは好ましくないか、もしくは制御することができない。従って、エンジン1の運転モードは通常運転モードではなく、フラグf/prism/onの値を「0」にセットする(STEPd−10)。
【0218】
さらに、エンジン1の回転数NE及び吸気圧PBがそれぞれ所定範囲内(正常な範囲内)にあるか否かの判別が行われる(STEPd−7,d−8)。このとき、いずれかが所定範囲内にない場合には、排気側主演算処理部13が生成する目標空燃比KCMDを使用してエンジン1の燃料供給を制御することは好ましくない。従って、エンジン1の運転モードは通常運転モードではなく、フラグf/prism/onの値を「0」にセットする(STEPd−10)。
【0219】
そして、STEPd−1,d−2,d−7,d−8の条件が満たされ、且つ、STEPd−3,d−4,d−5,d−6の条件が成立していない場合(これは、エンジン1の通常的な運転状態である)に、エンジン1の運転モードが通常運転モードであるとして、フラグf/prism/onの値を「1」にセットする(STEPd−9)。
【0220】
図11の説明に戻って、上記のようにフラグf/prism/onの値を設定した後、機関側制御ユニット7bは、フラグf/prism/onの値を判断し(STEPe)、f/prism/on=1である場合には、排気側主演算処理部13で生成された最新の目標空燃比KCMDを読み込む(STEPf)。また、f/prism/on=0である場合には、目標空燃比KCMDを所定値に設定する(STEPg)。この場合、目標空燃比KCMDとして設定する所定値は、例えばエンジン1の回転数NEや吸気圧PBからあらかじめ定めたマップ等を用いて決定する。
【0221】
次いで、機関側制御ユニット7bは、前記局所的フィードバック制御部16において、前述の如くオブザーバ21によりLAFセンサ5の出力KACTから推定した各気筒毎の実空燃比#nA/F(n=1,2,3,4)に基づき、PID制御器22により、各気筒毎のばらつきを解消するようにフィードバック補正係数#nKLAFを算出する(STEPh)。さらに、大局的フィードバック制御部15により、フィードバック補正係数KFBを算出する(STEPi)。
【0222】
この場合、大局的フィードバック制御部15は、前述の如く、PID制御器17により求められるフィードバック操作量KL AFと、適応制御器18により求められるフィードバック操作量KSTRを目標空燃比KCMDで除算してなるフィードバック操作量kstrとから、切換部20によってエンジン1の運転状態等に応じていずれか一方のフィードバック操作量KLAF又はkstrを選択する(通常的には適応制御器18側のフィードバック操作量kstrを選択する)。そして、選択したフィードバック操作量KLAF又はkstrを燃料噴射量を補正するためのフィードバック補正量数KFBとして出力する。
【0223】
尚、フィードバック補正係数KFBを、PID制御器17側のフィードバック操作量KLAFから適応制御器18側のフィードバック操作量kstrに切り換える際には、該補正係数KFBの急変を回避するために、適応制御器18は、その切換えの際の制御サイクルに限り、補正係数KFBを前回の補正係数KFB(=KLAF)に保持するように、フィードバック操作量KSTRを求める。同様に、補正係数KFBを、適応制御器18側のフィードバック操作量kstrからPID制御器17側のフィードバック操作量KLAFに切り換える際には、PID制御器17は、自身が前回の制御サイクルで求めたフィードバック操作量K LAFが、前回の補正係数KFB(=kstr)であったものとして、今回の補正係数KLAFを算出する。
【0224】
上記のようにしてフィードバック補正係数KFBが算出された後、さらに、前記STEPfあるいはSTEPgで決定された目標空燃比KCMDに応じた第2補正係数KCMDMが第2補正係数算出部10により算出される(STEPj)。
【0225】
次いで、機関側制御ユニット7bは、前述のように求められた基本燃料噴射量Timに、第1補正係数KTOTAL、第2補正係数KCMDM、フィードバック補正係数KFB、及び各気筒毎のフィードバック補正係数#nKLAFを乗算することで、各気筒毎の出力燃料噴射量#nToutを求める(STEPk)。そして、この各気筒毎の出力燃料噴射量#nToutが、付着補正部23によって、エンジン1の吸気管における燃料の壁面付着を考慮した補正を施された後(STEPm)、エンジン1の図示しない燃料噴射装置に出力される(STEPn)。
【0226】
そして、エンジン1にあっては、各気筒毎の出力燃料噴射量#nToutに従って、各気筒への燃料噴射が行われる。
【0227】
以上のような各気筒毎の出力燃料噴射量#nToutの算出及びそれに応じたエンジン1への燃料噴射がエンジン1のクランク角周期に同期したサイクルタイムで逐次行われ、これによりLAFセンサ5の出力KACT(空燃比の検出値)が、目標空燃比KCMDに収束するように、エンジン1の空燃比が制御される。この場合、特に、フィードバック補正係数KFBとして、適応制御器18側のフィードバック操作量kstrを使用している状態では、エンジン1の運転状態の変化や特性変化等の挙動変化に対して、高い安定性を有して、LAFセンサ5の出力KACTが迅速に目標空燃比KCMDに収束制御される。また、エンジン1が有する応答遅れの影響も適正に補償される。
【0228】
一方、前述のようなエンジン1の燃料供給の制御と並行して、前記排気側制御ユニット7aは、一定周期の制御サイクルで図13のフローチャートに示すメインルーチン処理を実行する。
【0229】
すなわち、図12を参照して、排気側制御ユニット7aは、まず、前記排気側主演算処理部13における演算処理を実行するか否かの判別処理を行って、その実行の可否を規定するフラグf/prism/calの値を設定する(STEP1)。尚、フラグf/prism/calの値は、それが「0」のとき、排気側主演算処理部13における演算処理を行わないことを意味し、「1」のとき、排気側主演算処理部13における演算処理を行うことを意味する。
【0230】
上記の判別処理は、図14のフローチャートに示すように行われる。
【0231】
すなわち、O2センサ6及びLAFセンサ5が活性化しているか否かの判別が行われる(STEP1−1,1−2)。このとき、いずれかが活性化していない場合には、排気側主演算処理部13の処理に使用するO2センサ6及びLAFセンサ5の検出データを精度よく得ることができないため、フラグf/prism/calの値を「0」にセットする(STEP1−6)。さらにこのとき、同定器25の後述する初期化を行うために、その初期化を行うか否かを規定するフラグf/id/resetの値を「1」にセットする(STEP1−7)。ここで、フラグf/id/resetの値は、それが「1」であるとき、同定器25の初期化を行うことを意味し、「0」であるとき、初期化を行わないことを意味する。
【0232】
また、エンジン1のリーン運転中(希薄燃焼運転中)であるか否か(STEP1−3)、及びエンジン1の始動直後の触媒装置3の早期活性化を図るためにエンジン1の点火時期が遅角側に制御されているか否か(STEP1−4)の判別が行われる。これらのいずれかの条件が成立している場合には、O2センサ6の出力VO2/OUTを目標値VO2/TARGETに整定させるような目標空燃比KCMDを算出しても、それをエンジン1の燃料供給の制御に使用することはないので、フラグf/prism/calの値を「0」にセットする(STEP1−6)。さらにこのとき、同定器25の初期化を行うために、フラグf/id/resetの値を「1」にセットする(STEP1−7)。
【0233】
図13の説明に戻って、上記のような判別処理を行った後、排気側制御ユニット7aは、さらに、同定器25による前記ゲイン係数a1,a2,b1の同定(更新)処理を実行するか否かの判別処理を行って、その実行の可否を規定するフラグf/id/calの値を設定する(STEP2)。尚、フラグf/id/calの値は、それが「0」のとき、同定器25による前記ゲイン係数a1,a2,b1の同定(更新)処理を行わないことを意味し、「1」のとき、同定(更新)処理を行うことを意味する。
【0234】
このSTEP2の判別処理では、エンジン1のスロットル弁が略全開であるか否か、及びエンジン1への燃料供給の停止中(フュエルカット中)であるか否かの判別が行われる。これらのいずれかの条件が成立している場合には、前記ゲイン係数a1,a2,b1を適正に同定することが困難であるため、フラグf/id/calの値を「0」にセットする。そして、上記のいずれの条件も成立していない場合には、同定器25による前記ゲイン係数a1,a2,b1の同定(更新)処理を実行すべくフラグf/id/calの値を「1」にセットする。
【0235】
次いで、排気側制御ユニット7aは、エンジン1の現在の回転数NEと吸気圧PBとからあらかじめ定めたマップによりO2センサ6の出力VO2/OUTの目標値VO2/TARGETの今回値を求めた後(STEP3)、前記減算処理部11,12によりそれぞれ最新の前記偏差出力kact(k)(=KACT(k)−FLAF/BASE)及びVO2(k)(=VO2/OUT(k)−VO2/TARGET)を算出する(STEP4)。この場合、減算処理部11,12は、前記図11のSTEPaにおいて取り込まれて図示しないメモリに記憶されたLAFセンサ5の出力KACT及びO2センサ6の出力VO2/OUTの時系列データの中から、最新のものを選択して前記偏差出力kact(k)及びVO 2(k)を算出する。また、特に、偏差出力VO2(k)の算出においては、今回の制御サイクルにおいてSTEP4で求めた目標値VO2/TARGETが用いられる。そしてこの偏差出力kact(k)及びVO2(k)は、排気側制御ユニット7aにおいて、過去に算出したものを含めて時系列的に図示しないメモリに記憶保持される。
【0236】
次いで、排気側制御ユニット7aは、前記STEP1で設定されたフラグf/prism/calの値を判断する(STEP5)。このとき、f/prism/cal=0である場合、すなわち、排気側主演算処理部13の演算処理を行わない場合には、スライディングモード制御器27で求めるべき前記SLD操作入力Usl(目標偏差空燃比kcmd)を強制的に所定値に設定する(STEP14)。この場合、該所定値は、例えばあらかじめ定めた固定値(例えば「0」)あるいは前回の制御サイクルで決定したSLD操作入力Uslの値である。
【0237】
尚、このようにSLD操作入力Uslを所定値とした場合において、排気側制御ユニット7aは、その所定値のSLD操作入力Uslに前記基準値FLAF/BASE を加算することで、今回の制御サイクルにおける目標空燃比KC MDを決定し(STEP15)、今回の制御サイクルの処理を終了する。
【0238】
一方、STEP4の判断で、f/prism/cal=1である場合、すなわち、排気側主演算処理部13の演算処理を行う場合には、排気側制御ユニット7aは、まず、前記目標空燃比算出手段13aの同定器25による演算処理を行う(STEP6)。
【0239】
この同定器25による演算処理は図15のフローチャートに示すように行われる。
【0240】
すなわち、同定器25は、まず、前記STEP2で設定されたフラグf/id/calの値を判断する(STEP6−1)。このときf/id/cal=0であれば、同定器25によるゲイン係数a1,a2,b1の同定処理を行わないので、直ちに図13のメインルーチンに復帰する。
【0241】
一方、f/id/cal=1であれば、同定器25は、さらに該同定器25の初期化に係わる前記フラグf/id/resetの値(これは、前記STEP1でその値が設定される)を判断し(STEP6−2)、f/id/reset=1である場合には、同定器25の初期化を行う(STEP6−3)。この初期化では、前記同定ゲイン係数a1ハット,a2ハット,b1ハットの各値があらかじめ定めた初期値に設定され(式(4)の同定ゲイン係数ベクトルΘの初期化)、また、前記式(9)で用いる行列P(対角行列)の各成分があらかじめ定めた初期値に設定される。さらに、フラグf/id/resetの値は「0」にリセットされる。
【0242】
次いで、同定器25は、現在の同定ゲイン係数a1(k-1)ハット,a2(k-1)ハット,b1(k-1)ハットの値と、前記STEP3で制御サイクル毎に算出される偏差出力VO2及びkactの過去値のデータVO2(k-1),VO2(k-2),kact(k-d-1)とを用いて、前記式(3)により前記同定偏差出力VO2(k)ハットを算出する(STEP6−4)。
【0243】
さらに同定器25は、新たな同定ゲイン係数a1ハット,a2ハット,b1ハットを決定する際に使用する前記ベクトルKθ(k)を式(9)により算出した後(STEP6−5)、前記同定誤差id/e(k)(前記同定偏差出力VO2ハットと、実際の偏差出力VO2との偏差。式(7)参照)を算出する(STEP6−6)。
【0244】
ここで、前記同定誤差id/e(k)は、基本的には、前記式(7)に従って算出すればよいのであるが、本実施形態では、前記図13のSTEP4で制御サイクル毎に算出する偏差出力VO2と、前記STEP6−4で制御サイクル毎に算出する同定偏差出力VO2ハットとから式(7)の演算により得られた値(=VO2(k)−VO2(k)ハット)に、さらにローパス特性のフィルタリングを施すことで同定誤差id/e(k)を求める。
【0245】
これは、触媒装置3を含む対象排気系Eの挙動は一般にローパス特性を有するため、前記排気系モデルのゲイン係数a1,a2,b1を適正に同定する上では、対象排気系Eの低周波数側の挙動を重視することが好ましいからである。
【0246】
尚、このようなフィルタリングは、結果的に、偏差出力VO2及び同定偏差出力VO2ハットの両者に同じローパス特性のフィルタリングが施されていればよく、例えば偏差出力VO2及び同定偏差出力VO2ハットにそれぞれ各別にフィルタリングを施した後に式(7)の演算を行って同定誤差id/e(k)を求めるようにしてもよい。また、前記のフィルタリングは、例えばディジタルフィルタの一手法である移動平均処理によって行われる。
【0247】
次いで、同定器25は、STEP6−6で求めた同定誤差id/e(k)と、前記STEP6−5で算出したKθ(k)とを用いて前記式(8)により新たな同定ゲイン係数ベクトルΘ(k)、すなわち、新たな同定ゲイン係数a1(k)ハット,a2(k)ハット,b1(k)ハットを算出する(STEP6−7)。
【0248】
このようにして新たな同定ゲイン係数a1(k)ハット,a2(k)ハット,b1(k)ハットを算出した後、同定器25は、該同定ゲイン係数a1ハット,a2ハット,b1ハット(同定ゲイン係数ベクトルΘの要素)の値を、所定の条件を満たすように制限する処理を行う(STEP6−8)。そして、同定器25は次回の制御サイクルの処理のために前記行列P(k)を前記式(10)により更新した後(STEP6−9)、図13のメインルーチンの処理に復帰する。
【0249】
この場合、上記STEP6−8において同定ゲイン係数a1ハット,a2ハット,b1ハットの値を制限する処理は、同定ゲイン係数a1ハット,a2ハットの値の組み合わせを所定の組み合わせに制限する処理(同定ゲイン係数a1ハット,a2ハットを成分とする座標平面上の所定の領域内に点(a1ハット,a2ハット)を制限する処理)と、同定ゲイン係数b1ハットの値を所定の範囲内に制限する処理とからなる。前者の処理では、STEP6−7で算出した同定ゲイン係数a1(k)ハット,a 2(k)ハットにより定まる上記座標平面上の点(a1(k)ハット,a2(k)ハット)が該座標平面上にあらかじめ定めた所定の領域から逸脱している場合に同定ゲイン係数a1(k)ハット,a2(k)ハットの値を強制的に上記所定の領域内の点の値に制限する。また、後者の処理では、前記STEP6−7で算出した同定ゲイン係数b1(k)ハットの値が所定の範囲の上限値あるいは下限値を超えている場合に、該同定ゲイン係数b1(k)ハットの値を強制的にその上限値あるいは下限値に制限する。
【0250】
このような同定ゲイン係数a1ハット,a2ハット,b1ハットの制限処理は、スライディングモード制御器27が算出するSLD操作入力Usl(目標偏差空燃比kcmd)、ひいては目標空燃比KCMDの安定性を確保するためのものである。
【0251】
尚、このような同定ゲイン係数a1ハット,a2ハット,b1ハットの制限処理のより具体的な手法については、本願出願人が例えば特開平11−153051号公報にて詳細に説明しているので、ここでは詳細な説明を省略する。
【0252】
また、図15のSTEP6−7で新たな同定ゲイン係数a1(k)ハット,a2(k)ハット,b1(k)ハットを求めるために使用する同定ゲイン係数の前回値a1(k-1)ハット,a2(k-1)ハット,b1(k-1)ハットは、前回の制御サイクルにおけるSTEP6−8の制限処理を行った後の同定ゲイン係数の値である。
【0253】
以上が図13のSTEP6における同定器25の演算処理の詳細である。
【0254】
図13の説明に戻って、上記のように同定器25の演算処理を行った後、排気側制御ユニット7aは、ゲイン係数a1,a2,b1の値を決定する(STEP7)。この処理では、前記STEP2で設定されたフラグf/id/calの値が「1」である場合、すなわち、同定器25によるゲイン係数a1,a2,b1の同定処理を行った場合には、ゲイン係数a1,a2,b1の値として、それぞれ前記STEP5で前述の通り同定器25により求められた最新の同定ゲイン係数a1(k)ハット,a2(k)ハット,b1(k)ハット(STEP6−8の制限処理を施したもの)を設定する。また、f/id/cal=0である場合、すなわち、同定器25によるゲイン係数a1,a2,b1の同定処理を行わなかった場合には、ゲイン係数a1,a2,b1の値をそれぞれあらかじめ定めた所定値とする。
【0255】
次いで、排気側制御ユニット7aは、前記推定器26による演算処理(推定偏差出力VO2バーの算出処理)を行う(STEP8)。
【0256】
すなわち、推定器26は、まず、前記STEP7で決定されたゲイン係数a1,a2,b1(これらの値は基本的には、前記同定ゲイン係数a1ハット,a2ハット,b1ハットである)を用いて、前記式(13)で使用する係数値α1,α2,βj(j=1,…,d)を、式(12)の但し書きの定義に従って算出する。
【0257】
次いで、推定器26は、前記STEP4で制御サイクル毎に算出されるO2センサ6の偏差出力VO2の現在の制御サイクル以前の時系列データVO2(k),VO2(k-1)、並びにLAFセンサ5の偏差出力kactの現在の制御サイクル以前の時系列データkact(k-j)(j=0,…,d1)と、スライディングモード制御器27から制御サイクル毎に与えられる前記目標偏差空燃比kcmd(=SLD操作入力Usl)の前回の制御サイクル以前の時系列データkcmd(k-j)(=Usl(k-j)。j=1,…,d2−1)と、上記の如く算出した係数値α1,α2,βjとを用いて前記式(13)により、推定偏差出力VO2(k+d)バー(今回の制御サイクルの時点から前記合計無駄時間d後の偏差出力VO2の推定値)を算出する。
【0258】
排気側制御ユニット7aは、次に、スライディングモード制御器27によって、前記SLD操作入力Usl(=目標偏差空燃比kcmd)を算出する(STEP9)。
【0259】
すなわち、スライディングモード制御器27は、まず、前記STEP8で推定器26により求められた推定偏差出力VO2バーの時系列データVO2(k+d)バー,VO2(k+d-1)バーを用いて、前記式(25)により定義された切換関数σバーの今回の制御サイクルから前記合計無駄時間d後の値σ(k+d)バー(これは、式(15)で定義された線形関数σの合計無駄時間d後の推定値に相当する)を算出する。
【0260】
尚、この場合、切換関数σバーの値があらかじめ定めた所定の許容範囲内に収まるようにし、上記の如く求められるσ(k+d)バーがその許容範囲の上限値又は下限値を超えた場合には、それぞれσバーの値σ(k+d)バーを強制的に該上限値又は下限値に制限する。これは、切換関数σバーの値が過大になると、前記到達則入力Urchが過大になると共に、前記適応則入力Uadpの急変を生じるために、O2センサ6の出力VO2/OUTの目標値VO2/TARGETへの収束制御の安定性が損なわれる虞れがあるからである。
【0261】
さらに、スライディングモード制御器27は、上記切換関数σバーの値σ(k+d)バーに、排気側制御ユニット7aの制御サイクルの周期ΔT(一定周期)を乗算したものσ(k+d)バー・ΔTを累積的に加算していく、すなわち、前回の制御サイクルで求められた加算結果に今回の制御サイクルで算出されたσ(k+d)バーと周期ΔTとの積σ(k+d)バー・ΔTを加算することで、前記式(27)のΣ(σバー・ΔT)の項の演算結果であるσバーの積算値(以下、この積算値をΣσバーにより表す)を算出する。
【0262】
尚、この場合、本実施形態では、上記積算値Σσバーがあらかじめ定めた所定の許容範囲内に収まるようにし、該積算値Σσバーがその許容範囲の上限値又は下限値を超えた場合には、それぞれ該積算値Σσバーを強制的に該上限値又は下限値に制限する。これは、積算値Σσバーが過大になると、前記式(27)により求められる適応則入力Uadpが過大となって、O2センサ6の出力VO2/OUTの目標値VO2/TARGETへの収束制御の安定性が損なわれる虞れがあるからである。
【0263】
次いで、スライディングモード制御器27は、前記STEP8で推定器26により求められた推定偏差出力VO2バーの現在値及び過去値の時系列データVO2(k+d)バー,VO2(k+d-1)バーと、上記の如く求めた切換関数σバーの値σ(k+d)バー及びその積算値Σσバーと、STEP6で決定したゲイン係数a1,a2,b1(これらの値は基本的には、最新の同定ゲイン係数a1(k)ハット,a2(k)ハット,b1(k)ハットである)とを用いて、前記式(24),(26),(27)に従って、それぞれ等価制御入力Ueq、到達則入力Urch及び適応則入力Uadpを算出する。
【0264】
そして、スライディングモード制御器27は、この等価制御入力Ueq、到達則入力Urch及び適応則入力Uadpを加算することで、前記SLD操作入力Usl、すなわち、O2センサ6の出力VO2/OUTを目標値VO2/TARGETに収束させるために必要な対象排気系Eへの入力量(=目標偏差空燃比kcmd)を算出する。
【0265】
上記のようにSLD操作入力Uslを算出した後、排気側制御ユニット7aは、スライディングモード制御器27による適応スライディングモード制御の安定性(より詳しくは、適応スライディングモード制御に基づくO2センサ6の出力VO2/OUTの制御状態(以下、SLD制御状態という)の安定性)を判別する処理を行って、該SLD制御状態が安定であるか否かをそれぞれ値「1」、「0」で表すフラグf/sld/stbの値を設定する(STEP10)。
【0266】
この安定性の判別処理は図16のフローチャートに示すように行われる。
【0267】
すなわち、排気側制御ユニット7aは、まず、前記STEP8で算出される切換関数σバーの今回値σ(k+d)バーと前回値σ(k+d-1)バーとの偏差Δσバー(これは切換関数σバーの変化速度に相当する)を算出する(STEP10−1)。
【0268】
次いで、排気側制御ユニット7aは、上記偏差Δσバーと切換関数σバーの今回値σ(k+d)バーとの積Δσバー・σ(k+d)バー(これはσバーに関するリアプノフ関数σバー2/2の時間微分関数に相当する)があらかじめ定めた所定値ε(≧0)以下であるか否かを判断する(STEP10−2)。
【0269】
ここで、上記積Δσバー・σ(k+d)バー(以下、これを安定判別パラメータPstbという)について説明すると、この安定判別パラメータPstbの値がPstb>0となる状態は、基本的には、切換関数σバーの値が「0」から離間しつつある状態である。また、安定判別パラメータPstbの値がPstb≦0となる状態は、基本的には、切換関数σバーの値が「0」に収束しているか、もしくは収束しつつある状態である。そして、一般に、スライディングモード制御ではその制御量を目標値に安定に収束させるためには、切換関数の値が安定に「0」に収束する必要がある。従って、基本的には、前記安定判別パラメータPstbの値が「0」以下であるか否かによって、それぞれ前記SLD制御状態が安定、不安定であると判断することができる。
【0270】
但し、安定判別パラメータPstbの値を「0」と比較することでSLD制御状態の安定性を判断すると、切換関数σバーの値に僅かなノイズが含まれただけで、安定性の判別結果に影響を及ぼしてしまう。このため、本実施形態では、前記STEP10−2で安定判別パラメータPstbと比較する所定値εは、「0」よりも若干大きな正の値としている。
【0271】
そして、STEP10−2の判断で、Pstb>εである場合には、SLD制御状態が不安定であるとし、前記STEP9で算出されるSLD操作入力Uslを用いた目標空燃比KCMDの決定を所定時間、禁止するためにタイマカウンタtm(カウントダウンタイマ)の値を所定の初期値TMにセットする(タイマカウンタtmの起動。STEP10−4)。さらに、前記フラグf/sld/stbの値を「0」に設定した後(STEP10−5)、図13のメインルーチンの処理に復帰する。
【0272】
一方、前記STEP10−2の判断で、Pstb≦εである場合には、排気側制御ユニット7aは、切換関数σバーの今回値σ(k+d)バーがあらかじめ定めた所定範囲内にあるか否かを判断する(STEP10−3)。
【0273】
この場合、切換関数σバーの今回値σ(k+d)バーが、所定範囲内に無い状態は、該今回値σ(k+d)バーが「0」から大きく離間している状態であるので、SLD制御状態が不安定であると考えられる。このため、STEP10−3の判断で、切換関数σバーの今回値σ(k+d)バーが、所定範囲内に無い場合には、SLD制御状態が不安定であるとして、前述の場合と同様に、STEP10−4及び10−5の処理を行って、タイマカウンタtmを起動すると共に、フラグf/sld/stbの値を「0」に設定する。
【0274】
尚、本実施形態では、前述のSTEP9の処理において、切換関数σバーの値を所定の許容範囲内に制限するので、STEP10−3の判断処理は省略してもよい。
【0275】
また、STEP10−3の判断で、切換関数σバーの今回値σ(k+d)バーが、所定範囲内にある場合には、排気側制御ユニット7aは、前記タイマカウンタtmを所定時間Δtm分、カウントダウンする(STEP10−6)。そして、このタイマカウンタtmの値が「0」以下であるか否か、すなわち、タイマカウンタtmを起動してから前記初期値TM分の所定時間が経過したか否かを判断する(STEP10−7)。
【0276】
このとき、tm>0である場合、すなわち、タイマカウンタtmが計時動作中でまだタイムアップしていない場合は、STEP10−2あるいはSTEP10−3の判断でSLD制御状態が不安定であると判断されてから、さほど時間を経過していないので、SLD制御状態が不安定なものとなりやすい。このため、STEP10−7でtm>0である場合には、前記STEP10−5の処理を行って前記フラグf/sld/stbの値を「0」に設定する。
【0277】
そして、STEP10−7の判断でtm≦0である場合、すなわち、タイマカウンタtmがタイムアップしている場合には、SLD制御状態が安定であるとして、フラグf/sld/stbの値を「1」に設定する(STEP10−8)。
【0278】
以上のような処理によって、SLD制御状態の安定性が判断され、不安定であると判断した場合には、フラグf/sld/stbの値が「0」に設定され、安定であると判断した場合には、フラグf/sld/stbの値が「1」に設定される。
【0279】
尚、以上説明したSLD制御状態の安定性の判断の手法は例示的なもので、この他の手法によって安定性の判断を行うようにすることも可能である。例えば制御サイクルよりも長い所定期間毎に、各所定期間内における前記安定判別パラメータPstbの値が前記所定値εよりも大きくなる頻度を計数する。そして、その頻度があらかじめ定めた所定値を超えるような場合にSLD制御状態が不安定であると判断し、逆の場合に、SLD制御状態が安定であると判断するようにしてもよい。
【0280】
図13の説明に戻って、上記のようにSLD制御状態の安定性を示すフラグf/sld/stbの値を設定した後、排気側制御ユニット7aは、フラグf/sld/stbの値を判断する(STEP11)。このとき、f/sld/stb=1である場合、すなわち、SLD制御状態が安定であると判断された場合には、スライディングモード制御器27が前記STEP9で算出したSLD操作入力Uslのリミット処理を行う(STEP12)。このリミット処理では、STEP9で算出されたSLD操作入力Uslの今回値Usl(k)が所定の許容範囲内にあるか否かが判断され、該今回値Usl(k)がその許容範囲の上限値又は下限値を超えている場合には、それぞれ、SLD操作入力Uslの今回値Usl(k)が強制的に該上限値又は下限値に制限される。
【0281】
尚、STEP12のリミット処理を経たSLD操作入力Usl(=目標偏差空燃比kcmd)は、図示しないメモリに時系列的に記憶保持され、それが、推定器26の前述の演算処理のために使用される。
【0282】
次いで、排気側制御ユニット7aは、触媒装置3の劣化状態を評価する処理(詳細は後述する)を排気側主演算処理部13の劣化状態評価手段13bによって実行した後(STEP13)、スライディングモード制御器27によって、STEP12のリミット処理を経たSLD操作入力Uslに前記基準値FLAF/BASEを加算することで、前記目標空燃比KCMDを算出し(STEP15)、今回の制御サイクルの処理を終了する。
【0283】
また、前記STEP11の判断でf/sld/stb=0である場合、すなわち、SLD制御状態が不安定であると判断された場合には、排気側制御ユニット7aは、今回の制御サイクルにおけるSLD操作入力Uslの値を強制的に所定値(固定値あるいはSLD操作入力Uslの前回値)に設定した後(STEP14)、前記式(28)に従って前記目標空燃比KCMDを算出し(STEP15)、今回の制御サイクルの処理終了する。
【0284】
尚、STEP15で最終的に決定される目標空燃比KCMDは、制御サイクル毎に図示しないメモリに時系列的に記憶保持される。そして、前記大局的フィードバック制御器15等が、排気側制御ユニット7aで決定された目標空燃比KCMDを用いるに際しては(図11のSTEPfを参照)、上記のように時系列的に記憶保持された目標空燃比KCMDの中から最新のものが選択される。
【0285】
前記STEP13における触媒装置3の劣化状態の評価処理を次に説明する。
【0286】
この評価処理は、排気側制御ユニット7aの劣化状態評価手段13bによって図17のフローチャートに示すように行われる。
【0287】
すなわち、劣化状態評価手段13bは、まず、前記式(15)により定めた劣化評価用線形関数σの値を、前記図13のSTEP4で算出されるO2センサ6の偏差出力VO2の時系列データVO2(k),VO2(k-1)(偏差出力VO2の現在値と1制御サイクル前の過去値)から算出する(STEP13−1)。
【0288】
尚、この場合に使用する係数s1,s2の値は、前記STEP9でスライディングモード制御器27が切換関数σバーの値を求めるために用いる係数s1,s2の値と同一である。
【0289】
次いで、劣化状態評価手段13bは、フラグF/DONEの値を判断する(STEP13−2)。ここで、このフラグF/DONEは、エンジン1の現在の運転中に、触媒装置3の劣化状態の評価を完了したか否かをそれぞれ値「1」、「0」で表すフラグであり、後述のSTEP13−5の処理においてその値が「1」に設定されるものである。そして、該フラグF/DONEは、エンジン1の始動時にその値が「0」に初期化されるものである。
【0290】
このとき、F/DONE=0である場合(触媒装置3の劣化状態の評価が未完了の場合)には、劣化状態評価手段13bは、排ガスボリューム(排気管2を流れる排ガスの流量)の変動状態を判断する処理を行う(STEP13−3)。この処理は、より詳しくは、排ガスボリュームがほぼ一定に維持されている状態(以下、クルーズ状態という)であるか否かを判断し、該クルーズ状態であるか否かをそれぞれ値「1」、「0」で表すフラグF/CRSの値を設定する処理である。この場合、本実施形態では、この処理は、排気側制御ユニット7aの制御サイクルの周期(30〜100ms)よりも長い周期(例えば1秒。以下、ここでは排ガスボリューム変動判断周期という)で、図18のフローチャートに示すように行われる。
【0291】
すなわち、まず、エンジン1の現在の回転数NE及び吸気圧PBの検出データから、次式(42)によって、現在の排ガスボリュームの推定値A BSV(以下、推定排ガスボリュームという)を算出する(STEP13−3−1)。
【0292】
【数42】
【0293】
ここで、本実施形態では、エンジン1の回転数が1500rpmであるときの排ガスボリュームを基準としているため、上記式(42)では、回転数NE(検出値)を「1500」で除算している。また、式(42)において、SVPRAは、エンジン1の排気量等に応じてあらかじめ定めた定数である。
【0294】
尚、排ガスボリュームは、上記のように推定する他、例えばエンジン1の燃料供給量や吸入空気量から推定したり、あるいはフローセンサを用いて直接的に検出するようにしてもよい。
【0295】
次いで、前記排ガスボリューム変動判断周期毎にSTEP13−3−1で算出される推定排ガスボリュームABSVに所定のフィルタリング処理を施すことによって、排ガスボリュームの変動状態を表す排ガスボリューム変動パラメータSVMAを求める(STEP13−3−2)。
【0296】
この場合、上記フィルタリング処理は、次式(43)により与えられる。
【0297】
【数43】
【0298】
すなわち、排ガスボリューム変動判断周期毎の推定排ガスボリュームABSVの変化量の複数周期分(本実施形態では3周期分)の移動平均を求めることにより、排ガスボリューム変動パラメータSVMAを算出する。尚、式(43)中の「n」は、排ガスボリューム変動判断周期のサイクル番数を表すものである。
【0299】
このようにして排ガスボリューム変動パラメータSVMAを算出したとき、該変動パラメータSVMAは、推定排ガスボリュームABSVの変化速度を表すものとなる。従って、排ガスボリューム変動パラメータSVMAは、その値が「0」に近い程、推定排ガスボリュームABSVの経時的変化が小さい状態(推定排ガスボリュームABSVがほぼ一定である状態)であることを意味する。
【0300】
劣化状態評価手段13bは次に、上記排ガスボリューム変動パラメータSVMAの値を二乗したもの、すなわち、該変動パラメータSVMAの二乗値SVMA2をあらかじめ定めた所定値δと比較する(STEP13−3−3)。ここで、該所定値δは、「0」近傍の正の値である。
【0301】
このとき、SVMA2≧δである場合(現在の排ガスボリュームの変動が比較的大きい場合)には、タイマカウンタ(カウントダウンタイマ)TMCRSJUDの値をあらかじめ定めた初期値X/TMCRSJSTに設定する(STEP13−3−4)。さらに、排ガスボリュームの変動状態が前記クルーズ状態(略一定に維持されている状態)でないとして、前記フラグF/CRSの値を「0」に設定した後(STEP13−3−5)、図17のルーチン処理に復帰する。
【0302】
一方、STEP13−3−3の判断処理において、SVMA2<δである場合(現在の排ガスボリュームの変動が比較的小さい場合)には、この状態が継続する限り前記タイマカウンタTMCRSJUDの値を排ガスボリューム変動判断周期毎に所定値ずつカウントダウンする(STEP13−3−6)。そして、このタイマカウンタTMCRSJUDの値が「0」以下になったか否か、すなわち、該タイマカウンタTMCRSJUDがタイムアップしたか否かを判断する(STEP13−3−7)。
【0303】
このとき、TMCRSJUD≦0で、タイマカウンタTMCRSJUDがタイムアップしている場合には、排ガスボリュームの変動状態が前記クルーズ状態であると判断して、タイマカウンタTMCRSJUDの値を「0」に保持すると共に(STEP13−3−8)、前記フラグF/CRSの値を「1」に設定する(STEP13−3−9)。そして、図17のルーチン処理に復帰する。
【0304】
また、STEP13−3−7の判断処理で、TMCRSJUD>0で、タイマカウンタTMCRSJUDがタイムアップしていない場合には、前記STEP13−3−5でフラグF/CRSの値を「0」に設定した後、図17のルーチン処理に復帰する。
【0305】
以上説明した処理が図17のSTEP13−3の処理である。このような処理によって、前記排ガスボリューム変動パラメータSVMAの二乗値SVMA2がSVMA2<δとなる状態、すなわち、排ガスボリュームの変動が小さい状態が前記タイマカウンタTMCRSJUDの初期値X/TMCRSJSTに相当する時間(例えば10〜15秒)、継続した場合に、前記クルーズ状態であるとして、フラグF/CRSの値が「1」に設定される。そして、これ以外の場合には、排ガスボリュームの変動状態はクルーズ状態ではないとしてフラグF/CRSの値が「0」に設定される。
【0306】
このようなSTEP13−3の処理により、排ガスボリュームがほぼ一定に維持されているような状態を適正に把握することができる。尚、前記排ガスボリューム変動判断周期の1周期内における前記排気側制御ユニット7aの各制御サイクルでは、フラグF/CRSの値は一定に維持される。
【0307】
図17の説明に戻って、前記劣化状態評価手段13bは、次に、前記基本劣化評価用パラメータLSσ2とO2センサ6の出力VO2/OUTの平均値VO2/OUTAVEとを算出する処理を実行する(STEP13−4)。この処理は、図19のフローチャートに示すように行われる。
【0308】
すなわち、劣化状態評価手段13bは、基本劣化評価用パラメータLSσ2を算出するための所定の条件が成立しているか否かを判断する(STEP13−4−1)。ここで判断する条件としては、前記STEP13−3で設定されるフラグF/CRSの値や、前記図11のSTEPdで機関側制御ユニッット7bが設定するフラグf/prism/onの値がある。
【0309】
この場合、F/CRS=1である場合、すなわち排ガスボリュームの変動状態が前記クルーズ状態である場合には、基本劣化評価用パラメータLSσ2を算出するための条件(以下、劣化評価条件という)が成立していないとして、基本劣化評価用パラメータLSσ2や平均値VO2/OUTAVEの算出処理を行うことなく、直ちに図17のルーチン処理に復帰する。
【0310】
このように前記クルーズ状態、すなわち排ガスボリュームがほぼ一定に維持されている状態で、基本劣化評価用パラメータLSσ2の算出を行わないのは、次の理由による。すなわち、前記クルーズ状態では、O2センサ6の出力VO2/OUTは目標値VO2/TARGETに安定して維持されやすいため、前記劣化評価用線形関数σの値は、触媒装置3の劣化が進行した状態であってもその変化が生じにくい。このため、前記クルーズ状態では、劣化評価用線形関数σの値は、前記図5〜図7を参照して説明したような触媒装置3の劣化状態に応じた傾向を生じにくい。そこで、本実施形態では、前記クルーズ状態では、基本劣化評価用パラメータLSσ2の算出は行わないこととし、また、これに伴い、平均値V O2/OUTAVEの算出も行なわない。
【0311】
また、STEP13−4−1でf/prism/on=0である場合、すなわち排気側制御ユニット7aのスライディングモード制御器27が求める目標空燃比KCMDに応じてエンジン1の燃料供給の制御を行う通常運転モード以外の運転モードである場合にも、前記劣化評価条件が成立していないとして基本劣化評価用パラメータLSσ2や平均値VO2/OUTAVEの算出処理を行うことなく、直ちに図17のルーチン処理に復帰する。これは、劣化評価用パラメータLSσ2による触媒装置3の劣化状態の評価を適正に行う上では、前記スライディングモード制御器27が適応スラーディングモード制御により生成する目標空燃比KCMDに応じてエンジン1の空燃比を制御している状態で得られるO2センサ6の出力VO2/OUTのデータを用いて基本劣化評価用パラメータLSσ2を求めることが好ましいからである。
【0312】
尚、STEP13−4−1の条件判断では、フラグF/CRS及びフラグf/prism/onの値の他、例えばエンジン1を搭載している車両の車速が所定の範囲内にあるか否か、エンジン1の始動後、ある程度の時間が経過した状態であるか否か、触媒装置3が活性化しているか否か等の条件判断もなされ、これらの条件が満たされていない場合には前記劣化評価条件が成立していないとして基本劣化評価用パラメータLSσ2や平均値VO2/OUTAVEの算出処理を行うことなく、直ちに図17のルーチン処理に復帰する。
【0313】
一方、STEP13−4−1で劣化評価条件が成立している場合(このときF/CRS=0且つf/prism/on=1である)には、劣化状態評価手段13bは、前記図17のSTEP13−1で排気側制御ユニット7aの制御サイクル毎に求められる劣化評価用線形関数σの値の二乗値σ2を算出する(STEP13−4−2)。
【0314】
そして、この二乗値σ2の現在値σ2(k)と、基本劣化評価用パラメータLSσ2の現在値LSσ2(k-1)と、前記式(30)の漸化式により定まるゲインパラメータBPの現在値BP(k-1)とから、前記式(29)によって、新たな基本劣化評価用パラメータLSσ2(k)を算出する(STEP13−4−3)。
【0315】
さらに、ゲインパラメータBPの値を式(30)により更新した後(STEP13−4−4)、基本劣化評価用パラメータLSσ2及びゲインパラメータBPの更新回数(これは基本劣化評価用パラメータLSσ2を求めるために使用した劣化評価用線形関数σの値の個数に相当する)をカウントするカウンタCB1Pの値を「1」だけ増加させる(STEP13−4−5)。尚、カウンタCB1Pの値は、エンジン1の始動時に「0」に初期化される。
【0316】
さらに、劣化状態評価手段13bは、O2センサ6の出力VO2/OUTの平均値VO2/OUTAVEを算出(更新)した後(STEP13−4−6)、図17のルーチン処理に復帰する。この場合、平均値VO2/OUTAVEを算出は、STEP13−4−6の処理を実行する毎にO2センサ6の出力VO2/OUTの現在値を累積加算し、その累積加算値をその累積加算の総回数で除算することにより算出される。従って、平均値VO2/OUTAVEは、基本劣化評価用パラメータLSσ2の算出に関与したO2センサ6の出力VO2/OUTの平均値である。
【0317】
尚、上述のようにSTEP13−4−3,13−4−4でそれぞれ求める基本劣化評価用パラメータLSσ2、ゲインパラメータBP、並びに、STEP13−4−6で平均値VO2/OUTAVEを算出する際に用いるO2センサ6の出力VO2/OUTの累積加算値及びその累積加算の総回数の値は、エンジン1の運転停止中も失われることがないように、エンジン1の運転終了時に図示しないEEPROM等の不揮発性メモリに記憶される。そして、次回のエンジン1の運転時には、その記憶された基本劣化評価用パラメータLSσ2、ゲインパラメータBP、累積加算値及び累積加算の基本総回数の値がそれらの初期値として用いられる。また、エンジン1の最初の運転時における基本劣化評価用パラメータLSσ2及びゲインパラメータBPの初期値はそれぞれ「0」、「1」であり、O2センサ6の出力VO2/OUTの累積加算値及びその累積加算の総回数の値の初期値はそれぞれ「0」である。さらに前記カウンタCB1Pの値は、エンジン1の始動時に「0」に初期化される。
【0318】
図17の説明に戻って、上記のように基本劣化評価用パラメータLSσ2の値や平均値VO2/OUTAVEを算出(更新)した後、劣化状態評価手段13bは、基本劣化評価用パラメータLSσ2の補正と、その補正により得られる補正後劣化評価用パラメータCLSσ2に基づく触媒装置3の劣化状態の評価とを行う処理を実行する(STEP13−5)。この処理は図20のフローチャートに示すように行われる。
【0319】
すなわち、劣化状態評価手段13bは、前記ゲインパラメータBPの現在値BP(k)と前回値BP(k-1)とが略等しいか否か(ゲインパラメータBPがほぼ収束したか否か)の判断と、前記カウンタCB1Pの値が所定値CB1CAT以上になったか否か(劣化評価用パラメータLSσ2を求めるために使用した劣化評価用線形関数σの値の個数が所定値CB1CATに達したか否か)の判断とを行う(STEP13−5−1,13−5−2)。
【0320】
尚、本実施形態では、エンジン1の始動前に図示しない車両のバッテリが一旦取り外された場合やエンジン1の初回の運転時等のように、エンジン1の始動時に前回の運転時の劣化評価用パラメータLSσ2及びゲインパラメータBPのデータが保持されていない場合(それらの値が「0」に初期化されている場合)には、前記STEP13−5−2で前記カウンタCB1Pの値と比較する所定値は、上記のパラメータLSσ2及びBPが保持されている場合よりも、大きな値に設定される。
【0321】
上記STEP13−5−1,13−5−2の判断において、それらのいずれかの条件が満たされていない場合には、現在の制御サイクルにおいてSTEP13−4で求められた基本劣化評価用パラメータLSσ2は、劣化評価用線形関数σの値の二乗値σ2の中心値に未だ十分に収束していないと考えられるので、触媒装置3の劣化状態の評価を行うことなく、STEP13−5の処理を終了する。
【0322】
一方、STEP13−5−1,13−5−2のいずれの条件も満たされている場合には、STEP13−4で求められた基本劣化評価用パラメータLSσ2の今回値は、劣化評価用線形関数σの値の二乗値σ2の中心値を表すものとなっている。そして、この場合には、劣化状態評価手段13bは、O2センサ6の出力VO2/OUTの平均値VO2/OUTAVEの現在値(STEP13−4−6で算出した最新値)から、前記図9のデータテーブルに基づいて補正係数KRESSを求める(STEP13−5−3)。さらに劣化状態評価手段13bは、この補正係数KRESSを基本劣化評価用パラメータLSσ2の現在値に乗算することで該基本劣化評価用パラメータLSσ2を補正し、補正後劣化評価用パラメータCLSσ2を求める(STEP13−5−4)。
【0323】
次いで、劣化状態評価手段13bは、この補正後劣化評価用パラメータCLSσ2を前記図9に示した所定の閾値CATAGELMTと比較する(STEP13−5−5)。
【0324】
このときCLSσ2≧CATAGELMTである場合には、触媒装置3の劣化状態が前記「劣化進行状態」(触媒装置3の交換を要するかもしくはその交換時期が近い程度に触媒装置3が劣化した状態)であると判断して、その旨を前記劣化報知器29に報知させる(STEP13−5−6)。そして、今回のエンジン1の運転中における触媒装置3の劣化状態の評価が終了したか否かをそれぞれ値「1」、「0」で表す前記フラグF/DONEの値を「1」に設定した後(STEP13−5−7)、STEP13−5の処理を終了する。
【0325】
また、STEP13−5−5で、CLSσ2<CATAGELMTである場合には、触媒装置3の劣化状態は前記「未劣化状態」であるので、劣化報知器29による報知を行うことなく、前記STEP13−5−7でフラグF/DONEの値を「1」に設定し、STEP13−5の処理を終了する。
【0326】
以上説明した処理が、図13のSTEP13で劣化状態評価手段13bが行う処理である。
【0327】
以上説明した本実施形態の装置によれば、排気側主演算処理部13の目標空燃比算出手段13aによって、触媒装置3の下流側のO2センサ6の出力VO2/OUTを、エンジン1の運転状態(回転数NE、吸気圧PBの状態)に応じて設定された目標値VO2/TARGETに収束(整定)させるように、適応スライディングモード制御の処理を用いてエンジン1の目標空燃比(LAFセンサ5が検出する空燃比の目標値)が逐次決定される。さらに、この目標空燃比KCMDにLAFセンサ5の出力KACTを収束させるようにエンジン1の燃料噴射量を調整することによって、エンジン1の空燃比が目標空燃比KCMDにフィードバック制御される。これにより、O2センサ6の出力VO2/OUTが目標値VO2/TARGETに収束制御され、ひいては触媒装置3の最適な排ガス浄化性能を確保することができる。
【0328】
また、このようなエンジン1の空燃比の制御と並行して、排気側主演算処理部13の劣化状態評価手段13bが、O2センサ6の偏差出力VO2の時系列データから劣化評価用線形関数σを逐次求める。さらにその劣化評価用線形関数σの二乗値σ2の中心値(本実施形態では最小二乗中心値)としての基本劣化評価用パラメータLSσ2を、逐次型の統計処理アルゴリズム(本実施形態では重み付最小二乗法のアルゴリズム)を用いて求める。また、この基本劣化評価用パラメータLSσ2の算出に用いられたO2センサ6の出力VO2/OUTの平均値VO2/OUTAVEに応じて補正して、補正後劣化評価用パラメータCLSσ2を求め、この補正後劣化評価用パラメータCLSσ2をあらかじめ定めた閾値CATAGELMTと比較することによって、触媒装置3の劣化状態を評価する。
【0329】
これにより、エンジン1の運転状態に合わせて触媒装置3の最適な浄化性能を確保しながら触媒装置3の劣化状態を適正に評価することができることとなる。また、この評価のために用いる補正後劣化評価用パラメータCLSσ2は、劣化評価用線形関数σの二乗値σ2の中心値としての基本劣化評価用パラメータCLSσ2をO2センサ6の出力VO2/OUTの平均値VO2/OUTAVEに応じて、換言すれば、O2センサ6の出力VO2/OUTが平均的にどのようなレベルの目標値VO2/TARGETを中心にして制御されたかに応じて補正したものであるので、O2センサ6の出力VO2/OUTを制御する目標値VO2/TARGETによらずに、触媒装置3の劣化状態との相関性が高いものとなる。この結果、該補正後劣化評価用パラメータCLSσ2に基づく触媒装置3の劣化状態の評価を適正に行うことができる。
【0330】
また、本実施形態では、排ガスボリュームがほぼ一定に維持されるようなクルーズ状態、すなわち、排ガスボリュームの変動が小さく、劣化評価用線形関数σの値の変化も生じにくいような状況では、基本劣化評価用パラメータLSσ2の算出やO2センサ6の出力VO2/OUTの平均値VO2/OUTAVEの算出を行わない。そして、このような状況以外の状況で、基本劣化評価用パラメータLSσ2や平均値VO2/OUTAVEを算出し、それらの値により定まる補正後劣化評価用パラメータCLSσ2に基づいて触媒装置3の劣化状態を評価する。このため、触媒装置3の劣化状態を表すものとしての補正後劣化評価用パラメータCLSσ2の信頼性が高まり、触媒装置3の劣化状態の評価を正確に行うことができる。
【0331】
従って、本実施形態の装置によれば、触媒装置3の良好な浄化性能を確保しながら、該触媒装置3の劣化状態の評価を高い信頼性で良好に行うことができる。
【0332】
尚、以上説明した実施形態では、基本劣化評価用パラメータLSσ2を、O2センサ6の出力VO2/OUTの平均値VO2/OUTAVEに応じて補正するようにしたが、エンジン1の空燃比は、O2センサ6の出力VO2/OUTを目標値VO2/TARGETに収束させるように制御されるので、O2センサ6の出力VO2/OUTの平均値VO2/OUTAVEは、目標値VO2/TARGETの平均値とほぼ同等のものとなる。従って、O2センサ6の出力VO2/OUTの平均値VO2/OUTAVEの代りに、目標値VO2/TARGETの平均値に応じて基本劣化評価用パラメータLSσ2を補正するようにしてもよい。この場合の実施形態に関しては、前記第1実施形態におけるSTEP13−4−6(図19参照)で、出力VO2/OUTの平均値VO2 /OUTAVEを求める代りに目標値VO2/TARGETの平均値を求めるようにすると共に、STEP13−5−3で、該目標値VO2/TARGETの平均値から前記図9と同様の形態のデータテーブルに基づいて補正係数KRESSを求めるようにすればよい。そして、その他の処理は前記第1実施形態と同一でよい。
【0333】
次に本発明の第2実施形態を図21〜図23を参照して説明する。尚、本実施形態は、触媒装置3の劣化状態の評価に係わる処理の一部のみが前記第1実施形態と相違するものであるので、第1実施形態と同一構成部分及び同一処理部分については第1実施形態と同一の図面及び同一の参照符号を用いて詳細な説明を省略する。
【0334】
本実施形態では、前記排気側主演算処理部13の劣化状態評価手段13bは、前記図13のSTEP13における処理(触媒装置3の劣化状態の評価処理)を次のように行なう。
【0335】
すなわち、図21のフローチャートを参照して、劣化状態評価手段13bは、前記第1実施形態と全く同様にSTEP13−1〜13−4の処理を行なった後、前記触媒装置3の劣化状態を判別するための閾値CATAGELMT2の設定処理と、該基本劣化評価用パラメータLSσ2に基づく触媒装置3の劣化状態の評価を行なう処理とを実行する(STEP13−5’)。この処理は図22のフローチャートに示すように行われる。
【0336】
すなわち、劣化状態評価手段13bは、前記第1実施形態と全く同様にSTEP13−5−1及びSTEP13−5−2の判断処理を行なった後、STEP13−4(図21参照)で求めたO2センサ6の出力VO2/OUTの平均値VO2/OUTAVEの現在値(最新値)から、図23に示すようにあらかじめ定められたデータテーブルにより閾値CATAGELMT2を求める(STEP13−5−3’)。
【0337】
ここで、閾値CATAGELMT2は、触媒装置3の劣化状態が前記「未劣化状態」であるか「劣化進行状態」であるかを判断するために前記基本劣化評価用パラメータLSσ2と比較する閾値である。そして、この閾値CATAGELMT2を設定するための図23のデータテーブルでは、前記基本劣化評価用パラメータLSσ2がO2センサ6の出力VO2/OUTの目標値VO2/TARGETの変化によって受ける影響を考慮し、閾値CATAGELMT2の値は、O2センサ6の出力VO2/OUTの平均値VO2/OUTAVEがO2センサ6の出力VO2/OUTの感度がほぼ最大となる値、すなわち図2のVO2/OUT1であるとき、最大値となり、平均値VO2/OUTAVEのVO2/OUT1からの差|VO2/OUTAVE−VO2/OUT1|が大きくなるに伴って閾値CATAGELMTの値が小さくなるように定められている。つまり、目標値VO2/TARGETが図2のVO2/OUT1に近い程、O2センサ6の出力VO2/OUTの変動、ひいては前記劣化評価用線形関数σの値のばらつきが生じ易く、基本劣化評価用パラメータLSσ2の値が大きくなりやすいという傾向を考慮し、この傾向に合わせて閾値CATAGELMT2の値を増減させるように該閾値CATAGELMT2を設定している。
【0338】
そして、本実施形態では、劣化評価手段13bは、このように設定したCATAGELMT2を前記基本劣化評価用パラメータLSσ2とSTEP13−5−5’で比較し、LSσ2≧CATAGELMT2であるときには、触媒装置3の劣化状態が「劣化進行状態」であるとして、劣化報知器29による報知を行わしめる(STEP13−5−6)。そして、第1実施形態と同様に、STEP13−5−7でフラグF/DONEの値を「1」に設定した後、前記図13のルーチン処理に復帰する。また、LSσ2<CATAGELMT2であるときには、触媒装置3の劣化状態が「未劣化状態」であるとして、劣化報知器29による報知を行うことなく、STEP13−5−7の処理を実行した後、図13のルーチン処理に復帰する。
【0339】
以上説明した以外の構成及び処理は前記第1実施形態と全く同一である。
【0340】
かかる本実施形態の装置によれば、基本劣化評価用パラメータLSσ2は、触媒装置3の劣化状態だけでなく、O2センサ6の出力VO2/OUTの目標値VO2/TARGETの変化の影響も受けるが、該目標値VO2/TARGETの変化の影響に合わせて閾値CATAGELMT2を設定することにより、前記第1実施形態と同様に、目標値VO2/TARGETの変化の影響によらずに適正に触媒装置3の劣化状態の評価を行なうことができる。従って、本実施形態においても、触媒装置3の良好な浄化性能を確保しながら、触媒装置3の劣化状態の評価を適正に精度よく行なうことができる。
【0341】
尚、本実施形態では、閾値CATAGELMT2を設定するために、O2センサ6の出力VO2/OUTの平均値VO2/OUTAVEを用いたが、該平均値VO2/OUTAVEの代りに、O2センサ6の出力VO2/OUTの目標値VO2/TARGETの平均値を用いてもよい。この場合には、図21のSTEP13−4で、目標値VO2/TARGETの平均値を求めておき、図22のSTEP13−5−3’では、該目標値VO2/TARGETの平均値から、図23のデータテーブルと同様の形態のデータテーブルにより、閾値CATAGELMT2を求めるようにすればよい。
【0342】
また、本発明は前述した第1及び第2実施形態に限定されるものではなく、例えば次のような変形態様も可能である。
【0343】
前記第1及び第2実施形態では、劣化評価用線形関数σの二乗値σ2の最小二乗中心値を基本劣化評価用パラメータLSσ2としたが、例えば劣化評価用線形関数σの絶対値の最小二乗中心値を基本劣化評価用パラメータとして求めるようにしてもよい。この場合には、前記図17あるいは図21のSTEP12−4において劣化評価用線形関数σの二乗値σ2の代わりに該線形関数σの絶対値を求め、前記式(29)中の「σ2」をその絶対値で置き換えた演算処理を行うことによって、触媒装置3の劣化状態に対して前記基本劣化評価用パラメータLSσ2と同様の傾向を呈する基本劣化評価用パラメータを得ることができる。そして、この基本劣化評価用パラメータをO2センサ6の出力VO2/OUTあるいは目標値VO2/TARGETの平均値に応じて補正したものをあらかじめ定めた所定の閾値と比較し、もしくは、該基本劣化評価用パラメータを、O2センサ6の出力VO2/OUTあるいは目標値VO2/TARGETの平均値に応じて設定した閾値と比較することで、前記第1又は第2実施形態と同様に、触媒装置3の劣化状態を評価することができる。
【0344】
さらには、劣化評価用線形関数σの二乗値σ2あるいは絶対値の最小二乗中心値に限らず、それらの二乗値σ2もしくは絶対値の平均値等の中心値を基本劣化評価用パラメータとしてもよい。あるいは、例えば劣化評価用線形関数σの値の分散(より正確には「0」に対する分散で、劣化評価用線形関数σの値の二乗値σ2の平均値)や、標準偏差(分散の平方根)を基本劣化評価用パラメータとして求めるようにしてもよい。このようにしても、該基本劣化評価用パラメータは、触媒装置3の劣化状態に対して、前記基本劣化評価用パラメータLSσ2と同様の傾向を呈する。
【0345】
また、前記第1及び第2実施形態では、劣化評価用線形関数σは、O2センサ6の偏差出力VO2の二つの時系列データを変数成分とする前記式(15)により定めたが、さらに多くの偏差出力VO2の時系列データを変数成分とする線形関数により劣化評価用線形関数を定義してもよい。この場合、スライディングモード制御用の切換関数は劣化評価用線形関数に含まれる偏差出力VO2の時系列データを推定偏差出力VO2バーの時系列データで置き換えた線形関数により定義することが好適である。
【0346】
また、劣化評価用線形関数は、例えば式(15)の偏差出力VO2(k),VO2(k-1)をO2センサ6の出力VO2/OUT(k),VO2/OUT(k-1)に置き換えた式によって定めてもよい。この場合には、この劣化評価用線形関数の値の中心値は、基本的には「(s1+s2)・VO2/TARGET」となる。そして、例えばこの中心値(s1+s2)・VO2/TARGETと劣化評価用線形関数の値との偏差の二乗値あるいは絶対値の最小二乗中心値等、該中心値(s1+s2)・VO2/TARGETに対する劣化評価用線形関数の値のばらつき度合いを表すパラメータを基本劣化評価用パラメータとして求めれば、前記第1又第2実施形態と同様に触媒装置3の劣化状態を評価することができる。
【0347】
さらには、例えば前記式(25)の切換関数σバー、すなわち、O2センサ6の推定偏差出力VO2バーの時系列データを変数成分とする線形関数を劣化評価用線形関数として用いてもよい。但し、O2センサ6の偏差出力VO2の前記合計無駄時間d後の推定値である推定偏差出力VO2バーを変数成分とする切換関数σバーよりも、O2センサ6の実際の偏差出力VO2を変数成分とする式(15)の劣化評価用線形関数σを用いた方が、触媒装置3の実際の状態が該線形関数σにより良く反映されるので、評価結果の信頼性を高める上で好ましいと考えられる。
【0348】
また、前記第1実施形態では、触媒装置3の劣化状態を評価するために劣化評価用線形関数σの二乗値σ2を用いたが、例えば該線形関数σの値とその変化速度との積(これは前記SLD制御状態の安定性の判別のために前記STEP10で用いた安定判別パラメータPstbである)を用いて触媒装置3の劣化状態を評価するようにすることも可能である。この場合においても、例えば上記積の分散(より一般的には該積の値のばらつき度合いを表すもの)を劣化評価用パラメータとして求めれば、その劣化評価用パラメータの値に基づいて触媒装置3の劣化状態を評価することが可能である。
【0349】
また、前記第1実施形態では、触媒装置3の劣化状態を「劣化進行状態」と「未劣化状態」との二つに分けて評価するようにしたが、前記補正後劣化評価用パラメータCLSσ2と比較する閾値、あるいは、前記基本劣化評価用パラメータLSσ2と比較する閾値をさらに多くすれば、触媒装置3の劣化状態をさらに多くの劣化度合いに分別して評価するようにすることも可能である。そして、それぞれの劣化度合いに応じて各別の報知をするようにすることも可能である。
【0350】
また、前記第1及び第2実施形態では、目標空燃比KCMDを算出するために、適応スライディングモード制御の処理を用いたが、適応則(適応アルゴリズム)を用いないスライディングモード制御の処理を用いてもよい。この場合には、前記式(28)から適応則入力Uadpの項を除去した式によって、目標空燃比KCMDを求めるようにすればよい。
【0351】
また、前記第1及び第2実施形態では、目標空燃比KCMDの算出に際しては、推定器26により前記合計無駄時間dの影響を補償するようにしたが、前記空燃比操作系の無駄時間が無視できる程小さいような場合には、対象排気系Eの無駄時間d1の影響のみを補償するようにしてもよい。この場合には、推定器26は、前記式(12)の「kcmd」及び「d」をそれぞれ「kact」及び「d1」に置き換えた次式(44)を用いて、前記実施形態と同様にO2センサ6の偏差出力VO2の無駄時間d1後の推定値VO2(k+d1)を制御サイクル毎に逐次求める。
【0352】
【数44】
【0353】
さらに、この場合には、スライディングモード制御器27は、前記式(24)〜(27)で「d」を「d1」に置き換えた式によって、等価制御入力Usl、到達則入力Urch及び適応則入力Uadpを制御サイクル毎に求め、それらを加算することで、目標偏差空燃比kcmdを求める。このようにすることで、対象排気系Eの無駄時間d1の影響を補償した目標空燃比KCMDを求めることができる。
【0354】
尚、この場合における同定器25や劣化状態評価手段13b、機関側制御ユニット7bの処理は、前述の第1又は第2実施形態のものと同一でよい。
【0355】
また、空燃比操作系の無駄時間d1だけでなく対象排気系Eの無駄時間d1も無視できるほど小さいような場合には、前記推定器26を省略してもよい。この場合には、スライディングモード制御器27や同定器25の演算処理は、d=d1=0として行えばよい。
【0356】
また、前記第1及び第2実施形態では、同定器25を備えたが、前記排気系モデルのゲイン係数a1,a2,b1をあらかじめ定めた固定値としたり、エンジン1の回転数や吸気圧等からマップ等を用いて適宜、ゲイン係数a1,a2,b1の値を設定するようにしてもよい。
【0357】
また、前記第1及び第2実施形態では、O2センサ6の出力VO2/OUTを目標値VO2/TARGETに収束させるためのフィードバック制御手法としてスライディングモード制御の処理を用いたが、他のフィードバック制御手法の処理によってO2センサ6の出力VO2/OUTを目標値VO2/TARGETに収束制御しつつ、触媒装置3の劣化状態の評価を行うことも可能である。
【図面の簡単な説明】
【図1】本発明の内燃機関の空燃比制御装置の第1実施形態の全体的システム構成を示すブロック図。
【図2】図1の装置で使用するO2センサ及び空燃比センサの出力特性図。
【図3】図1の装置の要部の基本構成を示すブロック図。
【図4】図1の装置で用いるスライディングモード制御を説明するための説明図。
【図5】図1の装置で用いる触媒装置の劣化状態の評価手法を説明するための図。
【図6】図1の装置で用いる触媒装置の劣化状態の評価手法を説明するための図。
【図7】図1の装置で用いる触媒装置の劣化状態の評価手法を説明するための図。
【図8】図1の装置で用いる触媒装置の劣化状態の評価手法を説明するための図。
【図9】図1の装置で用いる触媒装置の劣化状態の評価手法を説明するための図。
【図10】図1の装置で用いる適応制御器を説明するためのブロック図。
【図11】図1の装置の機関側制御ユニットの処理を示すフローチャート。
【図12】図11のフローチャートのサブルーチン処理を示すフローチャート。
【図13】図1の装置の排気側制御ユニットの処理を示すフローチャート。
【図14】図13のフローチャートのサブルーチン処理を示すフローチャート。
【図15】図13のフローチャートのサブルーチン処理を示すフローチャート。
【図16】図13のフローチャートのサブルーチン処理を示すフローチャート。
【図17】図13のフローチャートのサブルーチン処理を示すフローチャート。
【図18】図17のフローチャートのサブルーチン処理を示すフローチャート。
【図19】図17のフローチャートのサブルーチン処理を示すフローチャート。
【図20】図17のフローチャートのサブルーチン処理を示すフローチャート。
【図21】本発明の内燃機関の空燃比制御装置の第2実施形態の処理を示すフローチャート。
【図22】図21のフローチャートのサブルーチン処理を示すフローチャート。
【図23】図22のフローチャートの要部の処理を説明するための線図。
【符号の説明】
1…エンジン(内燃機関)、2…排気管(排気通路)、3…触媒装置、6…O2センサ、7b…機関側制御ユニット(空燃比操作手段)、13a…目標空燃比算出手段(空燃比操作手段)、13b…劣化状態評価手段(パラメータ生成手段)。
Claims (4)
- 内燃機関の排気通路に設けた触媒装置の下流側に配置された酸素濃度センサと、該酸素濃度センサの出力を内燃機関の運転状態に応じて設定された所定の目標値に収束させるように該内燃機関で燃焼させる混合気の空燃比を操作する空燃比操作手段と、該空燃比操作手段の処理の実行中における前記酸素濃度センサの出力のデータからあらかじめ定めた所定のアルゴリズムにより前記触媒装置の劣化状態を評価するための劣化評価用パラメータを生成するパラメータ生成手段とを備えた内燃機関の空燃比制御装置において、
前記劣化評価用パラメータは、前記酸素濃度センサの出力の時系列データを変数成分として表した所定の劣化評価用線形関数の値のばらつき度合いを表すデータであり、
前記劣化評価用パラメータを前記所定の目標値又は前記酸素濃度センサの出力の平均値に応じて補正し、その補正後の劣化評価用パラメータの値に基づき前記触媒装置の劣化状態を評価する劣化状態評価手段を備えたことを特徴とする内燃機関の空燃比制御装置。 - 内燃機関の排気通路に設けた触媒装置の下流側に配置された酸素濃度センサと、該酸素濃度センサの出力を内燃機関の運転状態に応じて設定された所定の目標値に収束させるように該内燃機関で燃焼させる混合気の空燃比を操作する空燃比操作手段と、該空燃比操作手段の処理の実行中における前記酸素濃度センサの出力のデータからあらかじめ定めた所定のアルゴリズムにより前記触媒装置の劣化状態を評価するための劣化評価用パラメータを生成するパラメータ生成手段と、該劣化評価用パラメータの値を所定の判定値と比較することにより前記触媒装置の劣化状態を評価する劣化状態評価手段とを備えた内燃機関の空燃比制御装置において、
前記劣化評価用パラメータは、前記酸素濃度センサの出力の時系列データを変数成分として表した所定の劣化評価用線形関数の値のばらつき度合いを表すデータであり、
前記劣化状態評価手段は、前記判定値を前記所定の目標値又は前記酸素濃度センサの出力の平均値に応じて設定する手段を備え、その設定した判定値と前記劣化評価用パラメータの値を比較することにより前記触媒装置の劣化状態を評価することを特徴とする内燃機関の空燃比制御装置。 - 前記パラメータ生成手段は、前記劣化評価用線形関数の値の時系列データの各データ値と該劣化評価用線形関数の値の中心値としてあらかじめ定めた所定値との偏差の二乗値又は絶対値にローパス特性のフィルタリング処理を施すことにより前記劣化評価用パラメータを求めることを特徴とする請求項1又は2記載の内燃機関の空燃比制御装置。
- 前記空燃比操作手段は、前記酸素濃度センサの出力を前記目標値に収束させるようにスライディングモード制御の処理により前記空燃比を操作するための操作量を逐次生成しつつ、該操作量に応じて前記空燃比を操作する手段であり、
前記劣化評価用線形関数は、前記スライディングモード制御の処理に用いる切換関数に応じて定められた線形関数であることを特徴とする請求項1〜3のいずれか1項に記載の内燃機関の空燃比制御装置。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2001143800A JP4437626B2 (ja) | 2001-05-14 | 2001-05-14 | 内燃機関の空燃比制御装置 |
US10/142,831 US6711891B2 (en) | 2001-05-14 | 2002-05-13 | Apparatus for controlling air-fuel ratio of internal combustion engine |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2001143800A JP4437626B2 (ja) | 2001-05-14 | 2001-05-14 | 内燃機関の空燃比制御装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2002339783A JP2002339783A (ja) | 2002-11-27 |
JP4437626B2 true JP4437626B2 (ja) | 2010-03-24 |
Family
ID=18989878
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2001143800A Expired - Fee Related JP4437626B2 (ja) | 2001-05-14 | 2001-05-14 | 内燃機関の空燃比制御装置 |
Country Status (2)
Country | Link |
---|---|
US (1) | US6711891B2 (ja) |
JP (1) | JP4437626B2 (ja) |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7826954B2 (en) * | 2004-06-25 | 2010-11-02 | Honda Motor Co., Ltd. | System for monitoring sensor outputs of a gas turbine engine |
DE102005045932A1 (de) * | 2005-09-26 | 2007-03-29 | Robert Bosch Gmbh | Gerät und Verfahren zur Diagnose einer technischen Vorrichtung |
JP4697201B2 (ja) * | 2007-07-19 | 2011-06-08 | トヨタ自動車株式会社 | 内燃機関の異常検出装置 |
US7983542B2 (en) * | 2007-10-29 | 2011-07-19 | Smiths Medical Asd, Inc. | PID coefficient adjustment for respiratory heater closed loop control |
JP2009115012A (ja) * | 2007-11-08 | 2009-05-28 | Denso Corp | 内燃機関の空燃比制御装置 |
US8600647B2 (en) * | 2009-01-30 | 2013-12-03 | Toyota Jidosha Kabushiki Kaisha | Air-fuel ratio control apparatus of a multi-cylinder internal combustion engine |
US8511651B2 (en) | 2011-03-29 | 2013-08-20 | Smiths Medical Asd, Inc. | Heater unit humidification chamber monitor |
JP5616274B2 (ja) * | 2011-03-31 | 2014-10-29 | 本田技研工業株式会社 | 空燃比制御装置 |
US10227940B2 (en) * | 2012-07-17 | 2019-03-12 | Honda Motor Co., Ltd. | Control device for internal combustion engine |
DE102016006328A1 (de) * | 2016-05-24 | 2017-11-30 | GM Global Technology Operations LLC (n. d. Ges. d. Staates Delaware) | Verfahren und Vorrichtung zum Überprüfen eines Sauerstoffsensors |
WO2019074623A1 (en) * | 2017-10-09 | 2019-04-18 | Exxonmobil Upstream Research Company | AUTOMATIC TUNING CONTROL DEVICE AND METHOD |
JP6624319B1 (ja) * | 2019-02-18 | 2019-12-25 | トヨタ自動車株式会社 | 触媒劣化検出装置、触媒劣化検出システム、データ解析装置、内燃機関の制御装置、および中古車の状態情報提供方法 |
FR3119643B1 (fr) * | 2021-02-08 | 2022-12-23 | Psa Automobiles Sa | Procede de regulation du fonctionnement d'un catalyseur de gaz d'echappement produits par le moteur thermique d'un vehicule automobile |
US11624333B2 (en) | 2021-04-20 | 2023-04-11 | Kohler Co. | Exhaust safety system for an engine |
Family Cites Families (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2526640B2 (ja) | 1988-07-20 | 1996-08-21 | トヨタ自動車株式会社 | 内燃機関の触媒劣化判別装置 |
JP3065127B2 (ja) | 1991-06-14 | 2000-07-12 | 本田技研工業株式会社 | 酸素濃度検出装置 |
JP2678985B2 (ja) | 1991-09-18 | 1997-11-19 | 本田技研工業株式会社 | 内燃エンジンの空燃比制御装置 |
JP3075025B2 (ja) | 1993-06-30 | 2000-08-07 | スズキ株式会社 | 内燃機関の触媒劣化判定装置 |
JP3162553B2 (ja) | 1993-09-13 | 2001-05-08 | 本田技研工業株式会社 | 内燃機関の空燃比フィードバック制御装置 |
EP0643213B1 (en) | 1993-09-13 | 1998-03-04 | Honda Giken Kogyo Kabushiki Kaisha | Air-fuel ratio detection system for internal combustion engine |
JP3162585B2 (ja) | 1993-09-13 | 2001-05-08 | 本田技研工業株式会社 | 内燃機関の空燃比検出装置 |
JP3325392B2 (ja) | 1994-07-06 | 2002-09-17 | 本田技研工業株式会社 | 内燃エンジンの燃料噴射制御装置 |
JP3269945B2 (ja) | 1994-08-12 | 2002-04-02 | 本田技研工業株式会社 | 内燃機関の燃料噴射制御装置 |
DE69527478T2 (de) | 1994-08-12 | 2002-11-14 | Honda Giken Kogyo K.K., Tokio/Tokyo | Kraftstoffzuteilungssteuersystem für eine Brennkraftmaschine |
US5852930A (en) | 1996-04-05 | 1998-12-29 | Honda Giken Kogyo Kabushiki Kaisha | Air-fuel ratio control system for internal combustion engines |
JP3261059B2 (ja) | 1996-04-05 | 2002-02-25 | 本田技研工業株式会社 | 内燃機関の空燃比制御装置 |
US5966930A (en) * | 1996-08-22 | 1999-10-19 | Honda Giken Kogyo Kabushiki Kaisha | Catalyst deterioration-determining system for internal combustion engines |
JP3331159B2 (ja) | 1997-09-16 | 2002-10-07 | 本田技研工業株式会社 | プラントの制御装置 |
JP3592519B2 (ja) | 1997-09-16 | 2004-11-24 | 本田技研工業株式会社 | 内燃機関の排気系の空燃比制御装置及びプラントの制御装置 |
CA2337865C (en) * | 1998-07-17 | 2008-12-16 | Honda Giken Kogyo Kabushiki Kaisha | Method of determining the deterioration of a catalytic converter for purifying an exhaust gas |
-
2001
- 2001-05-14 JP JP2001143800A patent/JP4437626B2/ja not_active Expired - Fee Related
-
2002
- 2002-05-13 US US10/142,831 patent/US6711891B2/en not_active Expired - Lifetime
Also Published As
Publication number | Publication date |
---|---|
US20030024514A1 (en) | 2003-02-06 |
US6711891B2 (en) | 2004-03-30 |
JP2002339783A (ja) | 2002-11-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3967524B2 (ja) | 内燃機関の空燃比制御装置 | |
JP3992925B2 (ja) | 排ガス浄化用触媒装置の劣化判別方法 | |
JP3592519B2 (ja) | 内燃機関の排気系の空燃比制御装置及びプラントの制御装置 | |
JP3484088B2 (ja) | プラントの制御装置 | |
JP3621839B2 (ja) | プラントの制御装置 | |
JP3354088B2 (ja) | 内燃機関の排気系の空燃比制御装置 | |
JP3261038B2 (ja) | 内燃機関の空燃比制御装置 | |
JP3773684B2 (ja) | 内燃機関の空燃比制御装置 | |
JP4354068B2 (ja) | 内燃機関の排ガスの空燃比制御装置 | |
JP4437626B2 (ja) | 内燃機関の空燃比制御装置 | |
JP4490000B2 (ja) | 内燃機関の空燃比制御装置 | |
JP4312325B2 (ja) | 排ガス浄化用触媒装置の劣化状態評価方法 | |
JP3655145B2 (ja) | 多気筒内燃機関の空燃比制御装置 | |
JPH1193741A (ja) | 内燃機関の排気系の空燃比制御装置及びプラントの制御装置 | |
JP3913940B2 (ja) | 内燃機関の空燃比制御装置 | |
EP1403491A1 (en) | DEVICE, METHOD, AND PROGRAM RECORDING MEDIUM FOR CONTROL OF AIR−FUEL RATIO OF INTERNAL COMBUSTION ENGINE | |
JP3773859B2 (ja) | プラントの制御装置及び内燃機関の空燃比制御装置 | |
JP2006233973A (ja) | 制御装置 | |
JP3696570B2 (ja) | プラントの制御装置 | |
JP3748418B2 (ja) | 内燃機関の空燃比制御装置 | |
JP4074142B2 (ja) | 内燃機関の空燃比制御装置 | |
JP2003195907A (ja) | 制御装置 | |
JP2003020983A (ja) | 内燃機関の排気系の空燃比制御装置及びプラントの制御装置 | |
JP2003195906A (ja) | 制御装置 | |
JP2006177369A (ja) | 制御装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20071126 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20090630 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20090721 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20090908 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20091222 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20091224 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130115 Year of fee payment: 3 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130115 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140115 Year of fee payment: 4 |
|
LAPS | Cancellation because of no payment of annual fees |