JP4422716B2 - Air-fuel ratio control device for internal combustion engine - Google Patents
Air-fuel ratio control device for internal combustion engine Download PDFInfo
- Publication number
- JP4422716B2 JP4422716B2 JP2006346578A JP2006346578A JP4422716B2 JP 4422716 B2 JP4422716 B2 JP 4422716B2 JP 2006346578 A JP2006346578 A JP 2006346578A JP 2006346578 A JP2006346578 A JP 2006346578A JP 4422716 B2 JP4422716 B2 JP 4422716B2
- Authority
- JP
- Japan
- Prior art keywords
- air
- fuel ratio
- catalyst
- amount
- specific
- 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
-
- Y02T10/47—
Landscapes
- 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
本発明は、内燃機関の空燃比制御装置に関する。 The present invention relates to an air-fuel ratio control apparatus for an internal combustion engine.
従来より、内燃機関の排ガスを浄化するため、排気浄化触媒(三元触媒)を排気通路上に配設し、排気通路に設けた空燃比センサにより検出される空燃比に基づいて機関に供給される混合気の空燃比が理論空燃比となるようにフィードバック制御を行っている。これにより、窒素酸化物NOx、一酸化炭素CO、及び炭化水素HCが同時に低減される。内燃機関から排出される排ガスの浄化率を向上させるには、上述したフィードバック制御を精度良く行うことが有効である。 Conventionally, in order to purify exhaust gas from an internal combustion engine, an exhaust purification catalyst (three-way catalyst) is disposed on the exhaust passage, and is supplied to the engine based on an air-fuel ratio detected by an air-fuel ratio sensor provided in the exhaust passage. Feedback control is performed so that the air-fuel ratio of the air-fuel mixture becomes the stoichiometric air-fuel ratio. Thereby, nitrogen oxide NOx, carbon monoxide CO, and hydrocarbon HC are simultaneously reduced. In order to improve the purification rate of the exhaust gas discharged from the internal combustion engine, it is effective to perform the above-described feedback control with high accuracy.
また、更なる排気浄化性能の向上のため、排気浄化触媒の酸素吸蔵作用を効果的に利用する方法が検討されて来ている。この方法は、排ガス中の浄化すべき物質を酸化還元する際に、排気浄化触媒に吸蔵した酸素や排気浄化触媒の吸蔵機能を利用する。 In order to further improve the exhaust purification performance, methods for effectively utilizing the oxygen storage function of the exhaust purification catalyst have been studied. This method uses oxygen stored in the exhaust purification catalyst or the storage function of the exhaust purification catalyst when oxidizing and reducing the substance to be purified in the exhaust gas.
従来の内燃機関の空燃比制御装置は、このような酸素吸蔵作用に着目してなされた装置であって、排気浄化触媒全体に吸蔵される酸素量(酸素吸蔵量、酸素の実際充填度)を、吸入空気流量、空気中の酸素の含有量、及びラムダ値1(理論空燃比)からの偏差の積に基づいて推定し、この酸素吸蔵量をある固定された目標値となるように空燃比を制御するようになっている(特許文献1を参照。)。
しかしながら、上述した制御は酸素吸蔵量を目標値近傍に維持しようとする制御であり、排ガス中の着目した一つ又は複数の特定の成分の排出量を所望の値以下とするような空燃比制御、或いは、触媒全体の酸素吸蔵量ではなく同触媒の上流部分までにおける酸素吸蔵量を所望の値近傍に維持しようとする空燃比制御を行うことができない等の問題がある。 However, the control described above is control for maintaining the oxygen storage amount close to the target value, and the air-fuel ratio control is performed so that the emission amount of one or more specific components of interest in the exhaust gas is not more than a desired value. Alternatively, there is a problem that air-fuel ratio control cannot be performed in which the oxygen storage amount up to the upstream portion of the catalyst, not the oxygen storage amount of the entire catalyst, is maintained near a desired value.
そこで、発明者は、排気浄化触媒から流出する(排気浄化触媒により浄化される)排ガス中に含まれる特定成分の排出量(あるいは、排ガスの状態を表す代表値)を推定する手法を開発するとともに、この推定値を所定の目標状態とするように空燃比制御を行うことにより排気浄化性能を向上しえることを見い出した。 Accordingly, the inventor has developed a method for estimating the emission amount (or representative value representing the state of exhaust gas) of a specific component contained in the exhaust gas flowing out of the exhaust purification catalyst (purified by the exhaust purification catalyst). The present inventors have found that exhaust purification performance can be improved by performing air-fuel ratio control so that this estimated value is set to a predetermined target state.
本発明の目的は、排気浄化触媒から流出する排ガス中に含まれる特定成分の排出量(あるいは、排ガスの状態を表す代表値)を推定し、その推定値に基づいて空燃比の制御を行うことにより、排気浄化触媒による排ガスの浄化効率を向上し得る内燃機関の空燃比制御装置を提供することにある。 An object of the present invention is to estimate an emission amount (or a representative value representing the state of exhaust gas) of a specific component contained in exhaust gas flowing out from an exhaust purification catalyst, and to control an air-fuel ratio based on the estimated value. Thus, an object of the present invention is to provide an air-fuel ratio control device for an internal combustion engine that can improve the purification efficiency of exhaust gas by an exhaust purification catalyst.
本発明は、流入するガスが通過する空間と、同空間に露呈するとともに触媒機能をもたらす物質及び酸素吸蔵放出機能をもたらす物質を担持したコート層と、を含んでなる排気浄化触媒を排気通路に配設した内燃機関の空燃比制御装置であって、前記排気浄化触媒の一部又は全部を通過した前記機関の排ガスに含まれる特定成分の排出量に応じた値を推定値として推定する推定手段と、前記推定手段によって推定される推定値の少なくとも一つが所定の目標状態となるように、前記排気浄化触媒に流入する排ガスの空燃比を制御する空燃比制御手段と、を備えた内燃機関の空燃比制御装置を提供する。 The present invention, exhaust and space gas inflows passes, the exhaust gas purifying catalyst comprising a carrier and a coating layer material which results in material and oxygen storage function results in a catalyst function while exposed to the space passage An air-fuel ratio control apparatus for an internal combustion engine disposed in the engine, wherein estimation is performed by estimating a value according to an emission amount of a specific component contained in exhaust gas of the engine that has passed a part or all of the exhaust purification catalyst. And an air-fuel ratio control means for controlling the air-fuel ratio of the exhaust gas flowing into the exhaust purification catalyst so that at least one of the estimated values estimated by the estimating means is in a predetermined target state An air-fuel ratio control apparatus is provided.
これによれば、推定された排ガス中の特定成分の排出量(濃度等)に応じた空燃比制御がなされるので、排気浄化性能が向上する。 According to this, since the air-fuel ratio control is performed according to the estimated emission amount (concentration, etc.) of the specific component in the exhaust gas, the exhaust purification performance is improved.
この場合、前記特定成分は、前記排気浄化触媒に流入する前記機関の排ガスに含まれるとともに還元機能を有するリダクション成分及び同排ガスに含まれるとともに同リダクション成分に酸素を供給し得るストレージ成分のうちの少なくとも一つの成分であり、前記推定手段は、前記特定成分の物質収支を考慮して形成された推定モデルに基づいて前記推定値を推定するように構成され得る。 In this case, the specific component is included in the exhaust gas of the engine that flows into the exhaust purification catalyst and has a reduction function and a storage component that is included in the exhaust gas and can supply oxygen to the reduction component. It is at least one component, and the estimation means may be configured to estimate the estimated value based on an estimation model formed in consideration of a material balance of the specific component.
さらに、前記推定手段のモデルは、前記排気浄化触媒を排ガスの流れ方向に複数の領域に分割したうちの特定領域に着目し、同特定領域の前記空間に流入する前記特定成分の量と同特定領域の空間から流出する同特定成分の量、及び同特定領域の空間から同特定領域の前記コート層に拡散する同特定成分の量に基づいて形成され得る。 Further, the model of the estimating means focuses on a specific area of the exhaust purification catalyst divided into a plurality of areas in the exhaust gas flow direction, and has the same specific quantity as the amount of the specific component flowing into the space of the specific area. It may be formed based on the amount of the specific component flowing out from the space of the region and the amount of the specific component diffusing from the space of the specific region to the coat layer in the specific region.
さらに、前記推定手段の推定モデルは、前記特定領域の空間から前記特定領域のコート層に拡散する前記特定成分の量と同コート層で消費される同特定成分の量とに基づいて形成され得る。 Furthermore, the estimation model of the estimation means may be formed based on the amount of the specific component diffusing from the space of the specific region to the coat layer of the specific region and the amount of the specific component consumed in the coat layer. .
このようなモデルにより特定成分の排出量を推定する場合、前記特定成分は、酸素、又は酸素と一酸化炭素とすることが計算上簡便である。 When estimating the emission amount of a specific component using such a model, it is simple in calculation that the specific component is oxygen or oxygen and carbon monoxide.
以下、本発明による空燃比制御装置の実施形態について図面を参照しつつ説明する。図1は、かかる空燃比制御装置を有する内燃機関の構成を示している。この空燃比制御装置は内燃機関であるエンジン1を制御するものである。
Hereinafter, embodiments of an air-fuel ratio control apparatus according to the present invention will be described with reference to the drawings. FIG. 1 shows the configuration of an internal combustion engine having such an air-fuel ratio control device. This air-fuel ratio control apparatus controls the
エンジン1は多気筒エンジンであり、図1においては、その一気筒のみが断面図で示されている。このエンジン1は、点火プラグ2によって各シリンダ3内の混合気に対して点火を行うことにより駆動力を発生するようになっている。エンジン1は、吸気通路4を介して外部から吸入された空気とインジェクタ5から噴射された燃料とを混合することで形成される混合気をシリンダ3内に吸入するようになっている。吸気バルブ6は、シリンダ3の内部と吸気通路4との間を開閉するようになっている。排気通路7は、シリンダ3の内部で燃焼された混合気を排ガスとして排気する通路である。排気バルブ8は、シリンダ3の内部と排気通路7との間を開閉するようになっている。
The
スロットルバルブ9は、シリンダ3内に吸入される吸入空気量を調節するためのものであって、吸気通路4に配設されている。このスロットルバルブ9には、その開度を検出するスロットルポジションセンサ10が接続されている。また、スロットルバルブ9は、スロットルモータ11と連結されており、スロットルモータ11の駆動力により開閉されるようになっている。アクセルポジションセンサ12は、スロットルバルブ9の近傍に配設され、アクセルペダルの操作量(アクセル開度)を検出するようになっている。これらの構成に基き、スロットルバルブ9の開度は電子的に制御されるようになっている。即ち、本エンジン1には、電子制御スロットル方式が採用されている。
The throttle valve 9 is for adjusting the amount of intake air sucked into the cylinder 3, and is disposed in the intake passage 4. A
また、エンジン1は、吸入空気量(吸入空気流量)を検出するエアフローメータ13、クランクシャフトの位置を検出することでシリンダ3内のピストン15の位置やエンジン回転速度NEを求める信号を発生するクランクポジションセンサ14、同エンジン1のノッキングを検出するノックセンサ16、及び冷却水温を検出する水温センサ17を備えている。
The
排気浄化触媒(触媒コンバータ、以下、単に「触媒」と称呼することもある。)19は、排気通路7に配設されている。このような触媒は、排気通路上に複数設けられる場合もある。例えば、排気の流れに対して直列に複数設けられる場合や、排気通路が分岐されている構成を備えるエンジンに対しては、各分岐した排気通路に備えられる場合(より具体的には、4気筒エンジンにおいて、2気筒の排気管が一つにまとめられた箇所に触媒が一つ配置され、残りの2気筒の排気管が一つにまとめられた箇所にもう一つの触媒が配置される場合)である。本実施形態においては、各シリンダ3の排気管が一つにまとめられた箇所よりも下流側に一つの触媒19が配設されている。
An exhaust purification catalyst (catalytic converter, hereinafter may be simply referred to as “catalyst”) 19 is disposed in the exhaust passage 7. A plurality of such catalysts may be provided on the exhaust passage. For example, in the case where a plurality of exhaust passages are provided in series, or for an engine having a configuration in which the exhaust passage is branched, more specifically, in a case where each branch exhaust passage is provided (more specifically, four cylinder (In the engine, one catalyst is arranged at a place where the exhaust pipes of two cylinders are combined into one, and another catalyst is arranged at a place where the exhaust pipes of the remaining two cylinders are combined into one) It is. In the present embodiment, one
エンジン1は、更に、触媒19の温度を測定する触媒温度センサ21、チャコールキャニスタ23、チャコールキャニスタ23に捕集された燃料タンクからの蒸発燃料を吸気通路4にパージさせるパージコントロールバルブ24、触媒19の上流側に取り付けられた上流側空燃比センサ25、及び触媒19の下流側に取り付けられた下流側空燃比センサ26を備えている。
The
前記空燃比センサ25,26は、各取付位置における排ガス中の酸素濃度から同排ガスの空燃比(排ガスの排気空燃比、排気空燃比)をそれぞれ検出するようになっている。空燃比センサ25は、空燃比をリニアに検出するリニア空燃比センサであり、空燃比センサ26は、空燃比がリッチであるかリーンであるかを識別する濃淡電池型センサである。
The air-
上述した点火プラグ2、インジェクタ5、スロットルポジションセンサ10、スロットルモータ11、アクセルポジションセンサ12、エアフローメータ13、クランクポジションセンサ14、ノックセンサ16、水温センサ17、触媒温度センサ21、パージコントロールバルブ24、上流側空燃比センサ25、及び下流側空燃比センサ26は、エンジン1を制御する電子制御ユニット(ECU)18と接続されていて、ECU18からの信号に基づいて制御され、或いは、検出結果をECU18に対して送出するようになっている。
Spark plug 2, injector 5,
ECU18は、演算を行うCPU、演算結果などの各種情報を記録するRAM、バッテリによって記憶内容が保持されるバックアップRAM、及び制御プログラム等を格納したROMを備えていて、各種演算を行い、インジェクタ5から噴射する燃料噴射量の制御、点火時期の制御、酸素吸蔵量の算出、後述するモデル補正、及び上記センサの診断等を行うようになっている。
The
(触媒の酸素吸蔵作用)
次に、触媒19の構成及び酸素吸蔵作用について説明する。
(Catalyst oxygen storage)
Next, the configuration of the
触媒19は、図2に外観を示したように、断面が楕円形(断面積がdA一定)の柱状のモノリス触媒コンバータと称呼される三元触媒であり、軸に直交する平面で同触媒19を切断した拡大断面図である図3に示したように、セラミックの一種であるコージェライトからなる担体19aにより、その内部が軸方向に延在する軸線方向空間に細分されている。各軸線方向空間は、軸線に垂直な平面で切断すると略正方形状を有していて、セルとも称呼される。担体19aは、アルミナのコート層19bによりコーティングされていて、同コート層19bは白金(Pt)等の貴金属からなる活性成分(触媒成分)及びセリア(CeO2)等の成分を担持している。
As shown in FIG. 2, the
かかる触媒19は、同触媒19に流入するガスの空燃比が理論空燃比のときに未燃成分(HC,CO)を酸化し、同時に窒素酸化物(NOx)を還元する機能(これを「触媒機能」又は「酸化還元機能」と云う。)を有する。また、触媒19は、上記セリア等の成分を担持することにより、同触媒19に流入する排ガス中の酸素分子を吸蔵(貯蔵,吸着)及び放出する性質(以下、「酸素吸蔵機能」と称呼する。)を有していて、この酸素吸蔵機能により、空燃比が理論空燃比からある程度まで偏移したとしても、HC,CO、及びNOxを浄化することができる。
The
即ち、触媒19は、機関の吸入する混合気の空燃比(以下、単に「機関の空燃比」とも称呼する。この機関の空燃比は触媒に流入するガスの空燃比である「排気空燃比」と等しい。)がリーンとなって触媒19に流入するガスに過剰の酸素及び窒素酸化物NOxが多量に含まれると、同過剰の酸素を吸蔵するとともに、同窒素酸化物NOxから酸素を奪って(NOxを還元して)同酸素を吸蔵し、これによりNOxを浄化する。また、触媒19は、機関の空燃比がリッチになって同触媒19に流入するガスに炭化水素HCや一酸化炭素CO等の未燃成分が多量に含まれると、内部に吸蔵している酸素分子をこれらの未燃成分に与えて同未燃成分を酸化し、これによりHC,COを浄化する。
That is, the
従って、触媒19が酸素を吸蔵し得る限界まで酸素を吸蔵していると(換言すれば、酸素吸蔵量OSAが最大酸素吸蔵量OSAmax=Cmaxに到達していると)、同触媒19に流入するガスの排気空燃比がリーンとなったときに酸素を吸蔵することができないので、排ガス中のNOxを充分に浄化できなくなる。一方、触媒19が酸素を放出しきっていて酸素を全く吸蔵していなければ(換言すれば、酸素吸蔵量OSAが「0」になっていると)、同触媒19に流入するガスの排気空燃比がリッチとなったときに酸素を放出することができないので、排ガス中のHCやCOを充分に浄化できなくなる。このため、触媒19に流入するガスの空燃比が過渡的に相当のリーン空燃比又は相当のリッチ空燃比となった場合であっても、同ガス中の有害成分を充分に浄化できるように、触媒19の酸素吸蔵量OSAを精度良く推定するとともに、同酸素吸蔵量OSAを所定の値に維持するように同触媒19に流入するガスの空燃比を制御しておくことが望ましい。
Therefore, when the
また、触媒19にリーンな空燃比のガスが流入したときには、同触媒19の上流側においてより多くの酸素が吸蔵され、同触媒19にリッチな空燃比のガスが流入したときには、同触媒19の上流側から吸蔵されている酸素が消費されて行く。従って、同触媒19の最上流位置から任意の位置までの酸素吸蔵量OSAの合計量を推定し、同推定値に基く空燃比制御を行えば、触媒19全体の酸素吸蔵量が「0」または最大酸素吸蔵量Cmaxに到達することが回避し易くなるので、同空燃比の制御に不可避的な制御遅れが存在していても、エミッションを効果的に低減できる。さらに、触媒19から流出するガス(又は触媒19の一部又は全部により浄化されたガス)に含まれる特定のガス成分の濃度を推定できれば、同ガス成分に着目した空燃比制御を行うことで、同ガス成分の排出を高精度に抑制することができる。
Further, when a lean air-fuel ratio gas flows into the
以上の要求から、本空燃比制御装置は、推定モデル(触媒モデル)を使用することにより、特定ガス成分の濃度(排出量)と酸素吸蔵量を推定する。これらの推定値が、排気浄化触媒の全部又はその内部の所定領域(一部)から流出する排ガス中に含まれる少なくとも一つの特定成分の排出量又は流出する排ガスの状態を示す少なくとも一つの代表値である。 From the above requirements, the air-fuel ratio control apparatus estimates the concentration (discharge amount) and oxygen storage amount of a specific gas component by using an estimation model (catalyst model). These estimated values are at least one representative value indicating the emission amount of at least one specific component contained in the exhaust gas flowing out of the entire exhaust purification catalyst or a predetermined region (part) inside the exhaust purification catalyst or the state of the exhaust gas flowing out. It is.
(触媒モデル)
以下、触媒モデルについて説明すると、先ず、図4に模式的に示したように、触媒19をガスの入口(流入側)Frから出口(流出側)Rrに向う軸線に直交する面により複数の領域(「ブロック」とも称呼する。)に分割する。即ち、触媒19を排ガスの流れ方向に沿って、複数の領域に分割する。分割された各領域の軸線方向の長さはL(微小の長さでありdxとも書く。)である。上述したように、触媒19の断面積はdA一定である。
(Catalyst model)
Hereinafter, the catalyst model will be described. First, as schematically shown in FIG. 4, the
なお、この触媒モデルは、触媒を複数の領域に分割して構築されたモデルであるが、触媒全体を以下に述べる一つの領域として考えることにより、同触媒19を複数の領域に分割することなく触媒モデルを構築することもできる。
This catalyst model is a model constructed by dividing the catalyst into a plurality of regions. However, the
次に、分割された領域のうちの任意の特定領域に注目し、同特定領域を通過する特定の化学種(特定成分)の物質の収支を考える。化学種は、排ガス中に含まれる成分であり、例えば、酸素O2、一酸化炭素CO、炭化水素HC、及び窒素酸化物NOxである。なお、化学種は、触媒に流入する排ガスの空燃比がリッチのときに同排ガス中に含まれる成分を総合したもの(リッチ成分)、或いは、触媒に流入する排ガスの空燃比がリーンのときに同排ガス中に含まれる成分を総合したもの(リーン成分)とすることもできる。 Next, paying attention to an arbitrary specific area among the divided areas, a balance of a substance of a specific chemical species (specific component) passing through the specific area is considered. The chemical species are components contained in the exhaust gas, and are, for example, oxygen O 2 , carbon monoxide CO, hydrocarbon HC, and nitrogen oxide NOx. The chemical species is a combination of the components contained in the exhaust gas when the air-fuel ratio of the exhaust gas flowing into the catalyst is rich (rich component), or when the air-fuel ratio of the exhaust gas flowing into the catalyst is lean It can also be a combination of components contained in the exhaust gas (lean component).
ここで、本触媒モデルで使用する種々の値を次の表1のように定義する。 Here, various values used in the present catalyst model are defined as shown in Table 1 below.
いま、時刻t〜t+Δtの所与の期間における特定領域での化学種の収支を考えると、図5に示したように、特定領域の排ガス相(単に、「ガス相」とも称呼する。)における化学種の変化量ΔMは、下記の(1)式に示したとおり、同特定領域に流入した同化学種の量Minから、同特定領域から流出した同化学種の量Moutとコート層に奪われた同化学種の量Mcoatとを減算した量と等しい。このように、触媒モデルは特定領域における特定成分の物質収支に基づいて構築される。
以下、(1)式の各項について個別に検討する。先ず、(1)式の左辺にある化学種の変化量ΔMは、下記(2)式により求めることができる。(2)式は、上記所与の期間における化学種の濃度変化量(化学種の濃度Cgの時間変化量を所与の期間に渡り積分した量)に微小体積σ・dA・dxを乗じた値を特定領域に渡って(軸方向に)積分したものである。
(1)式の右辺第1項のMinは、単位時間あたりに特定領域に流入するガスの体積に相当する値である「特定領域に流入するガスの流速vginと同特定領域の断面積dAの積vgin・dA(実際には、断面積dAで開口率σの触媒内に流速Vginのガスが流れ込むので、触媒内部でのガスの流速はVgin/σとなり、この実際の流速Vgin/σと触媒の実質的な断面積σ・dAの積)」に同流入するガス中の化学種の濃度Cginを乗じた値Cgin・Vgin・dAを所与の期間に渡り積分した値である。また、(1)式の右辺第2項のMoutは、特定領域から流出するガスの流速vgoutと同特定領域の断面積dAの積vgout・dA(実際には、ガスの流速Vgout/σと実質的な断面積σ・dAの積)に同流出するガス中の同化学種の濃度Cgoutを乗じた値Cgout・vgout・dAを所与の期間に渡り積分した値である。即ち、上記(1)式の右辺第1項及び第2項は下記(3)式のように記述することができる。
ところで、特定領域に流入するガスの流速vginと同特定領域から流出するガスの流速vgoutとの間に大きな差異はないので、vg=vgin=vgoutとおくと、(3)式は、下記(4)式のように変形される。
次に、(1)式の右辺第3項のコート層に伝達される(移動する)化学種の量Mcoatについて検討する。幾何学的表面積Sgeoは触媒の単位体積あたりの化学種の反応に寄与する表面積であるから、特定領域において化学種の反応に寄与する表面積はSgeo・dA・dxであり、同特定領域の単位長あたりに同反応に寄与する面積はSgeo・dAとなる。また、コート層に伝達される化学種の量は、フィックの法則から、ガス相の化学種の濃度Cgとコート層の化学種の濃度Cwとの差に比例すると考えることができる。これらから、下記の(5)式が得られる。なお、hDは比例定数であるが、上記の表に示したように、物質伝達率と称呼される値である。
従って、上記(1),(2),(4)式、及び(5)式から、以下の(6)式が得られる。
この(6)式に準定常(quasi state)近似を適用すると、(6)式の左辺は「0」(∂Cg/∂t=0)であると考えることができるから(即ち、濃度Cgは瞬間的に定常値に至ると考えられるから)、下記の(7)式が得られる。
ここで、見かけの拡散速度(実質的な拡散速度)RDを(8)式のようにおけば、(7)式は(9)式に書き直される。
次に、特定領域のコート層における化学種の収支(特定成分の物質収支)を上記と同様に考えると、下記(10)式に示したように、コート層内における化学種の時間的変化量(単位時間あたりの変化量)ΔMcは、単位時間あたりに排ガス相からコート層へ伝達される同化学種の量Mdから、同単位時間あたりにコート層にて反応により消費される同化学種の量Mrを減じた量である。
(10)式の左辺(コート層内における化学種の時間的変化量)ΔMcは、下記(11)式に示したように、化学種の濃度変化(∂Cw/∂t)に体積((1−σ)・dA・dx)を乗じることにより求められ、右辺第1項(単位時間あたりに排ガス相からコート層へ伝達される化学種の量Md)は(5)式で説明した理由と同じ理由により、即ち、フィックの法則から考えると、下記(12)式のように記述することができる。
また、(10)式の右辺第2項(単位時間あたりにコート層にて反応により消費される化学種の量Mr)は、コート層での化学種の消費速度Rを用いた下記(13)式により求められる。
従って、(10)〜(13)式から、下記の(14)式が得られる。
この(14)式に準定常(quasi state)近似を適用すると(∂Cw/∂t=0)、下記の(15)式が得られる。
ここで、(15)式に(8)式を適用すれば、下記の(16)式が得られる。
以上を要約すると、(9)式及び(16)式が触媒モデルの基本式である。(9)式は、ある化学種の(特定領域への流入量)と(排ガス相からコート層への拡散量+特定領域からの流出量)とが釣り合っていることを示し、(16)式は、同化学種の(排ガス相からコート層への拡散量)と(コート層での消費量)とが釣り合っていることを示している。 In summary, Equations (9) and (16) are basic equations of the catalyst model. Equation (9) shows that (inflow amount into a specific area) of a certain chemical species is balanced with (amount of diffusion from the exhaust gas phase to the coating layer + outflow amount from the specific area). Indicates that the (amount of diffusion from the exhaust gas phase to the coat layer) and (consumption in the coat layer) of the same chemical species are balanced.
次に、かかる触媒モデルを使用して特定領域から流出する特定の化学種iの濃度Cgoutを実際に算出するための方法について説明する。先ず、(9)式を離散化すると、下記(17)式が得られる。なお、以下においては上記dxをLとして表す。
ここで、図6に概念的に示したように、特定領域Iから流出する化学種の濃度Cgoutは同特定領域Iの化学種の濃度Cg(I)の影響を強く受けると考えられるので、下記の(18)式のように置くことができる。かかる考え方は「風上法」と称呼される。換言すると、風上法とは、「特定領域Iに隣接する上流側の領域(I−1)における濃度Cg(I-1)の化学種が、特定領域Iに流入する。」という考え方であり、下記(19)式のように記述することもできる。
ところで、反応速度論に基けば、ある化学種の消費速度Rは、その化学種のコート層の平均濃度Cwの関数fcw(例えば、Cwのn乗)となるので、この関数fcwを最も簡便となるようにfcw(x)=xとおけば、同消費速度Rは(20)式にて示したように置くことができる。なお、以下において、(20)式中のR*を便宜上「消費速度定数」と称呼する。
この(20)式を上記(16)式(R=RD・(Cg-Cw)…(16))に適用すると下記(21)式が得られ、同(21)式を変形することにより下記(22)式が得られる。
また、上述した風上法によれば、Cg=Cgoutであるから、(22)式は下記(23)式に書き換えられる。
そして、Cg=Cgoutなる関係を上記(17)式に適用してCgを消去するとともに、同(17)式と上記(23)式とからCwを消去すると、下記(24)式が得られる。
そこで、値SPを下記(25)式のようにおけば、(24)式は(26)式のように書き直すことができる。値SPは、見かけの拡散速度RDと消費速度定数R*のうちの小さい方の値に強い影響を受ける値であるから、Cgoutの変化が物質の伝達(RD)又は化学的反応(R*)の何れにより律速されているかを示す値であり、従って、「反応律速因子」と呼ぶこともできる。
以上のことから、消費速度定数R*と見かけの拡散速度RDとを決定できれば、特定領域に流入する化学種濃度Cginを与えることにより、(25)式と(26)式とに基づいて同特定領域から流出する化学種の濃度Cgoutを求めることができる。以上が、化学種の濃度Cgoutを算出する基本的考え方である。 From the above, if the consumption rate constant R * and the apparent diffusion rate RD can be determined, the concentration of the chemical species flowing into the specific region Cgin is given, based on the equations (25) and (26). The concentration Cgout of the chemical species flowing out from the specific region can be obtained. The above is the basic concept for calculating the concentration Cgout of the chemical species.
次に、上記消費速度定数R*と見かけの拡散速度RDを決定するとともに、特定領域から流出する化学種濃度Cgoutを求める際のより具体的な方法の一例について説明する。この例(触媒モデル)では、触媒での酸化・還元反応である三元反応は瞬時に且つ完全に終了するものと仮定し、その結果としての酸素の過不足に基く酸素の吸蔵・放出反応に着目することとする。なお、この仮定(触媒モデル)は、現実的であり且つ精度の良いものである。 Next, an example of a more specific method when determining the consumption rate constant R * and the apparent diffusion rate R D and determining the chemical species concentration Cgout flowing out from the specific region will be described. In this example (catalyst model), it is assumed that the three-way reaction, which is an oxidation / reduction reaction at the catalyst, ends instantaneously and completely, and the resulting oxygen storage / release reaction is based on the excess or deficiency of oxygen. Pay attention. This assumption (catalyst model) is realistic and accurate.
この場合、着目する化学種iは、例えば、酸素O2や窒素酸化物の一つである一酸化窒素NOのように酸素を生成する(酸素をもたらす)化学種(ストレージ・エージェント)、及び、一酸化炭素COや炭化水素HCのように酸素を消費する化学種(リダクション・エージェント)から選ばれた化学種である。 In this case, the chemical species i of interest is, for example, a chemical species (storage agent) that generates oxygen (provides oxygen) such as oxygen O 2 or nitrogen monoxide NO, which is one of nitrogen oxides, and It is a chemical species selected from chemical species (reduction agents) that consume oxygen, such as carbon monoxide CO and hydrocarbon HC.
また、以下において、ストレージ・エージェントの化学種iのCgoutをCgout,stor,i、同化学種iのCwをCw,stor,i、同化学種iのCginをCgin,stor,i、同化学種iの見かけの拡散速度RDをRD,i、同化学種iの消費速度をRstor,i、同化学種iの消費速度定数をR*stor,i、及び同化学種iの反応律速因子SPstor,i(この場合、化学種iはO2又はNO)と表す。 Also, in the following, Cgout of storage agent chemical species i is Cgout, stor, i, Cw of chemical species i is Cw, stor, i, Cgin of chemical species i is Cgin, stor, i, chemical species i of the diffusion rate R D of the apparent R D, i, Rstor consumption rate of the chemical species i, i, the reaction rate-limiting factor for the consumption rate constant of the chemical species i R * stor, i, and the species i SPstor, i (in this case, chemical species i is O 2 or NO).
同様に、リダクション・エージェントの化学種iのCgoutをCgout,reduc,i、同化学種iのCwをCw,reduc,i、同化学種iのCginをCgin,reduc,i、同化学種iの見かけの拡散速度RDをRD,i、同化学種iの消費速度をRreduc,i、同化学種iの消費速度定数をR*reduc,i、及び同化学種iの反応律速因子SPreduc,i(この場合、化学種iはCO又はHC等)と表す。このように各値を表すと、上記(20),(23),(25),(26)式から以下の(27)〜(34)式が得られる。
これらの式に基づいて、Cgout,sotr,i(具体的には、特定領域から流出する酸素の濃度Cgout,O2、特定領域から流出する一酸化窒素の濃度Cgout,NO)及びCgout,reduc,i(具体的には、特定領域から流出する一酸化炭素の濃度Cgout,CO、特定領域から流出する炭化水素の濃度Cgout,HC)を求めるため、先ず、消費速度定数R*stor,i及び消費速度定数R*reduc,iを求める。 Based on these equations, Cgout, sotr, i (specifically, the concentration Cgout, O2 of oxygen flowing out from the specific region, the concentration Cgout, NO of nitrogen monoxide flowing out from the specific region) and Cgout, reduc, i (Specifically, in order to obtain the concentration Cgout, CO of carbon monoxide flowing out from the specific region and the concentration Cgout, HC of hydrocarbon flowing out from the specific region), first, the consumption rate constant R * stor, i and the consumption rate Find the constant R * reduc, i.
ところで、反応速度論によれば、特定領域のコート層で酸素が吸蔵される速度(酸素の吸蔵速度)Rstor,iは、同コート層のストレージ・エージェント(O2、NOx等)の濃度Cw,stor,i(例えば、Cw,O2、Cw,NO)の関数f1(Cw,stor,i)の値に比例するとともに、特定領域のコート層の最大酸素吸蔵密度と実際の(その時点における)酸素吸蔵密度との差(Ostmax-Ost)の関数f2(Ostmax-Ost)の値とに比例すると考えられる。この最大酸素吸蔵密度と酸素吸蔵密度との差(Ostmax-Ost)は、着目している特定領域における酸素吸蔵余裕量を表す。 By the way, according to the reaction kinetics, the rate at which oxygen is occluded in the coat layer in a specific region (oxygen occlusion rate) Rstor, i is the concentration Cw, storage agent (O 2 , NOx, etc.) of the coat layer. It is proportional to the value of the function f1 (Cw, stor, i) of stor, i (for example, Cw, O2, Cw, NO), and the maximum oxygen storage density of the coating layer in a specific region and the actual (at that time) oxygen It is considered to be proportional to the value of the function f2 (Ostmax-Ost) of the difference (Ostmax-Ost) with the occlusion density. The difference (Ostmax−Ost) between the maximum oxygen storage density and the oxygen storage density represents the oxygen storage margin in the specific region of interest.
そこで、簡単のために関数f1(x)=f2(x)=xとすると、下記の(35)式が得られる。下記(35)式のkstor,iは酸素吸蔵速度係数(吸蔵側反応速度係数,ストレージ・エージェントの消費速度係数)であって、よく知られたアレニウスの式で表される温度に依存して変化する係数であり、別途検出又は推定される触媒温度Tempと所定の関数(酸素吸蔵速度係数kstor,iと触媒温度Tempとの間の関係を規定したマップでも良い。)とに基づいて求めることができる。なお、酸素吸蔵速度係数kstor,iは、触媒劣化程度に応じても変化するので、同触媒劣化程度に応じて求めてもよい。
従って、(27)式と(35)式とから、消費速度定数R*stor,iは下記(36)式により求めることができる。
また、酸素の吸蔵(吸着)と放出のみに着目しているこのモデルにおいては、還元剤であるリダクション・エージェントはコート層に吸蔵されている酸素の放出のみに使用されるから、同リダクション・エージェントの消費速度Rredcu,iはコート層に吸蔵されている酸素が放出される速度(酸素の放出速度)Rrel,iと等しい。 In this model, which focuses only on oxygen storage (adsorption) and release, the reduction agent, which is a reducing agent, is used only to release oxygen stored in the coat layer. Consumption rate Rredcu, i is equal to the rate at which oxygen stored in the coat layer is released (oxygen release rate) Rrel, i.
そこで、酸素の放出速度Rrel,iについて検討すると、同放出速度Rrel,iは、酸素の吸蔵速度Rstor,iと同様に反応速度論に基き、同コート層において酸素を消費する化学種(例えば、CO,HC)の濃度Cw,reduc,i(例えば、Cw,CO、Cw,HC)の関数g1(Cw,reduc,i)の値に比例するとともに、酸素吸蔵密度Ostの関数g2(Ost)の値とに比例すると考えられる。 Therefore, when the oxygen release rate Rrel, i is examined, the release rate Rrel, i is based on the reaction kinetics in the same manner as the oxygen storage rate Rstor, i, and is a chemical species that consumes oxygen in the coat layer (for example, (CO, HC) concentration Cw, reduc, i (for example, Cw, CO, Cw, HC) is proportional to the value of function g1 (Cw, reduc, i), and oxygen storage density Ost of function g2 (Ost) It is considered to be proportional to the value.
そこで、簡単のために関数g1(x)=g2(x)=xとすると、下記の(37)式が得られる。下記(37)式のkrel,iは酸素放出速度係数(吸脱側反応速度係数)であって、酸素吸蔵速度係数kstor,iと同様にアレニウスの式で表される温度に依存して変化する係数であり、別途検出又は推定される触媒温度Tempに基づいて所定の関数(酸素放出速度係数krel,iと触媒温度Tempとの間の関係を規定したマップでも良い。)に基づいて求めることができる。なお、酸素放出速度係数krel,iは、触媒劣化程度に応じても変化するので、同触媒劣化程度に応じて求めてもよい。
この結果、上述したようにリダクション・エージェントの消費速度Rredcu,iはコート層の酸素の放出速度Rrel,iと等しいから、消費速度定数R*reduc,iは(31)式と(37)式とを比較することにより得られる下記(38)式に基づいて求めることができる。
以上のことから、酸素吸蔵密度Ostが求められれば(酸素吸蔵密度Ostの求め方については、後述する。)、(36)式から消費速度定数R*stor,i(例えば、R*O2)を求めることができる。一方、見かけの拡散速度RD,i(例えば、RD,O2)は(8)式のようにSgeo・hD,iであるから、温度と流量の関数(触媒の温度と同触媒を通過するガスの流量の関数)として実験的に求めておくことができる。この結果、(29)式からSPstor,i(例えば、SPstor,O2)が決定されるので、境界条件としてCgin,stor,i(例えば、Cgin,O2)が与えられるとき、(30)式からCgout,stor,i(例えば、Cout,O2)が求められる。そして、新たなCw,stor,i(例えば、Cw,O2)が(28)式により求められる。 From the above, if the oxygen storage density Ost is obtained (the method for obtaining the oxygen storage density Ost will be described later), the consumption rate constant R * stor, i (for example, R * O2) is calculated from the equation (36). Can be sought. On the other hand, the apparent diffusion rate R D , i (for example, R D , O 2) is Sgeo · h D , i as shown in equation (8), so it is a function of temperature and flow rate (the temperature of the catalyst passes through the same catalyst) As a function of the gas flow rate). As a result, SPstor, i (for example, SPstor, O2) is determined from the expression (29). Therefore, when Cgin, stor, i (for example, Cgin, O2) is given as the boundary condition, Cgout from the expression (30) , stor, i (eg, Cout, O2). Then, new Cw, stor, i (for example, Cw, O2) is obtained by the equation (28).
同様に、酸素吸蔵密度Ostが求められれば、(38)式から消費速度定数R*reduc,i(例えば、R*reduc,CO)を求めることができる。一方、見かけの拡散速度RD,i(例えば、RD,CO)は(8)式のようにSgeo・hD,iであるから、温度と流量の関数(触媒の温度と同触媒を通過するガスの流量の関数)として実験的に求めておくことができる。この結果、(33)式からSPreduc,i(例えば、SPreduc,CO)が決定されるので、境界条件としてCgin,reduc,i(例えば、Cgin,CO)が与えられるとき、(34)式からCgout,reduc,i(例えば、Cgout,CO)が求められる。そして、新たなCw,reudc,i(例えば、Cw,CO)が(32)式により求められる。 Similarly, if the oxygen storage density Ost is obtained, the consumption rate constant R * reduc, i (for example, R * reduc, CO) can be obtained from the equation (38). On the other hand, the apparent diffusion rate R D , i (for example, R D , CO) is Sgeo · h D , i as shown in Equation (8), so it is a function of temperature and flow rate (the temperature of the catalyst passes through the catalyst) As a function of the gas flow rate). As a result, SPreduc, i (for example, SPreduc, CO) is determined from the expression (33). Therefore, when Cgin, reduc, i (for example, Cgin, CO) is given as the boundary condition, Cgout from the expression (34) , reduce, i (eg, Cgout, CO). Then, a new Cw, reudc, i (for example, Cw, CO) is obtained by the equation (32).
次に、Cgout,stor,i、Cgout,reduc,iを求めるために必要となる酸素吸蔵密度Ostの求め方について説明する。 Next, how to obtain the oxygen storage density Ost required for obtaining Cgout, stor, i and Cgout, reduc, i will be described.
先ず、コート層での化学種としての酸素の収支について着目すると、同収支はコート層での酸素の吸蔵分と酸素の放出分の差であるから、下記(39)式により記述される。(39)式でdA・Lは特定領域の体積dVである。
この(39)式を変形すると、下記(40)式が得られる。
この(40)式を、(35)式と(37)式とを用いながら離散化すると、下記の(41)式が得られる。
この(41)式を変形すると、下記(42)式〜(44)式が得られ、これらから酸素吸蔵密度Ostを求めること(更新して行くこと)ができる。
このように、式(42)〜(44)式から酸素吸蔵密度Ostが求められるので、上述したようにCgout,stor,i、Cgout,reduc,iを求めることができる。また、酸素吸蔵密度Ostが求められるから、下記(45)式に基づいて特定領域の酸素吸蔵量OSAを求めることができる。
従って、触媒に流入する化学種濃度Cgin,iが境界条件として与えられたとき、触媒上流のブロック(特定領域)から、順次、(45)式を用いて各ブロックの酸素吸蔵量OSAを求めることができ、これにより、触媒内部の酸素吸蔵量の分布が精度良く推定される。また、各ブロックの酸素吸蔵量OSAを触媒全体について積算すれば、同触媒全体の酸素吸蔵量についても精度良く推定することができる。 Therefore, when the chemical species concentration Cgin, i flowing into the catalyst is given as a boundary condition, the oxygen storage amount OSA of each block is obtained sequentially from the block upstream of the catalyst (specific region) using equation (45). As a result, the distribution of the oxygen storage amount inside the catalyst is accurately estimated. Further, if the oxygen storage amount OSA of each block is integrated for the entire catalyst, the oxygen storage amount of the entire catalyst can also be accurately estimated.
次に、実際の触媒における酸素O2、一酸化炭素CO、炭化水素HC、及び窒素酸化物(ここでは一酸化窒素)NOの各化学種濃度Cgout,O2、Cgout,CO、Cgout,HC、及びCgout,NOを求めるための具体的方法の一例について、フローチャートを用いて説明する。この例においても、触媒での酸化・還元反応である三元反応は瞬時に且つ完全に終了するものと仮定し、その結果としての酸素の過不足に基く酸素の吸蔵・放出反応にのみ着目することとする。 Next, the chemical species concentrations Cgout, O2, Cgout, CO, Cgout, HC, and Cgout of oxygen O2, carbon monoxide CO, hydrocarbon HC, and nitrogen oxide (here, nitric oxide) NO in the actual catalyst , An example of a specific method for obtaining NO will be described using a flowchart. In this example as well, it is assumed that the three-way reaction, which is an oxidation / reduction reaction at the catalyst, is completed instantaneously and focuses only on the oxygen storage / release reaction based on the excess or deficiency of oxygen as a result. I will do it.
ECU18のCPUは、図7から図11の一連のフローチャートによりに示されたプログラムを所定時間の経過毎に図7のプログラムから順に実行するようになっている。また、これらのプログラムは、触媒19を多数のブロック(特定領域)に分割し、そのうちの一つのブロックIについてのプログラムであって、全く同様なプログラムが他のブロックについても同時並行的にCPUにより実行されている。
The CPU of the
CPUは所定のタイミングになると図7のステップ700から処理を開始し、ステップ705に進んで前回の本プログラムの演算時において後述するステップ745にて算出されたコート層の酸素濃度Cw,O2(k+1)を今回のコート層の酸素濃度Cw,O2の値であるCw,O2(k)に設定し、続くステップ710にて前回の本プログラムの演算時において後述する図11のステップ1125にて算出された酸素吸蔵密度Ost(k+1)を今回の酸素吸蔵密度Ostの値Ost(k)に設定する。
The CPU starts the process from
次に、CPUは、ステップ715にて触媒19の温度Tempと同触媒19の劣化程度を表す劣化指標値REKKAと図12に示したようなマップ(ルックアップテーブル)MapkstorO2とから酸素吸蔵速度係数kstor,O2(k)を決定する。
Next, in
触媒温度Tempは、温度センサ21によって検出してもよいし、エンジン1の運転状態(例えば、吸入空気量Gaとエンジン回転速度NE)に応じて推定してもよい。
The catalyst temperature Temp may be detected by the
劣化指標値REKKAは、触媒19の最大酸素吸蔵量Cmaxに応じて求められる値(例えば、最大酸素吸蔵量Cmaxの単調増加関数)である。最大酸素吸蔵量Cmaxは、次のようにして求められる。即ち、エンジン1が所定の定常運転状態にある場合に、下流側センサ26が理論空燃比よりリーンな空燃比を検出しているとき、触媒19に流入するガスの空燃比を所定のリッチ空燃比に維持し、これにより触媒19内に貯蔵されている酸素を総べて消費する。
The deterioration index value REKKA is a value obtained according to the maximum oxygen storage amount Cmax of the catalyst 19 (for example, a monotonically increasing function of the maximum oxygen storage amount Cmax). The maximum oxygen storage amount Cmax is obtained as follows. That is, when the
その結果、下流側センサ26が理論空燃比よりリーンな空燃比に対応する値から理論空燃比よりリッチな空燃比に対応する値を出力すると、その時点t1から触媒19に流入するガスの空燃比を所定のリーン空燃比に設定し、次に下流側センサ26が理論空燃比よりリッチな空燃比に対応する値から理論空燃比よりリーンな空燃比に対応する値を出力する時点t2までに、同触媒19に流入するガス中に含まれる酸素の量を下記(46)式、及び(47)式に基づいて求める。この(47)式により求められる積算値O2storage(=Cmax1)を最大酸素吸蔵量Cmaxとして採用する。
なお、時点t2の後、触媒19に流入するガスの空燃比を所定のリッチ空燃比に維持し、次に下流側センサ26が理論空燃比よりリーンな空燃比に対応する値から理論空燃比よりリッチな空燃比に対応する値を出力する時点t3までに、同触媒19に流入するガス中に含まれる単位時間あたりの酸素の不足量を上記(46)式と同様な式により求め、上記(47)式と同様な式により同単位時間あたりの酸素の不足量を時点t2〜t3に渡って積算し、これにより求められる積算値Cmax2を最大酸素吸蔵量Cmaxとして採用してもよく、前記最大酸素吸蔵量Cmax1とこの最大酸素吸蔵量Cmax2との平均値を最大酸素吸蔵量Cmaxとして採用してもよい。
Note that, after time t2, the air-fuel ratio of the gas flowing into the
次いで、CPUは、ステップ720にて同ステップ720内に記載した式(上記(36)式を参照。)に従って酸素の消費速度定数R*stor,O2(k)を求める。なお、ステップ720にて用いる最大酸素吸蔵密度Ostmaxは、一定値としてもよいが、前記触媒劣化指標値REKKA(又は、最大酸素吸蔵量Cmax)に応じて決定されることが望ましい(以下、同じ。)。その後、CPUはステップ725にて見かけの拡散速度RD,O2(k)を触媒温度TempとマップMapRDO2とから決定する。
Next, in
続いて、CPUはステップ730にて酸素の反応律速因子SPstor,O2を同ステップ730内に記載した式(上記(29)式を参照。)により求め、ステップ735にて、このプログラムが対象としているブロックIよりも前の(即ち、上流の)ブロックI−1から流出する酸素濃度Cgout,O2(k)を、同対象としているブロックIに流入する酸素濃度Cgin,O2(k)として取り込む。
Subsequently, the CPU obtains the oxygen reaction rate-determining factor SPstor, O2 in
このとき、対象としているブロックIが触媒19の最も上流のブロックである場合、前のブロックI−1は存在しないから、ステップ735における前のブロックのCgout,O2(k)は、同触媒19に流入するガスの酸素濃度Cgin,O2である。この触媒19に流入するガスの酸素濃度Cgin,O2は、同触媒19に流入するガスの空燃比A/Fと同ガスの流量とに基く関数fO2により求められる。下記の(48)式の右辺は、この関数fO2の具体例である。
上記(48)式の導出過程を簡単に述べると、触媒19に流入する排ガスの空燃比AFはGa/Gfであり、Gfに対して理論空燃比を得るために必要な空気質量をGastoichとすると、理論空燃比AFstoichはGastoich/Gfとなる。一方、供給燃料質量がGaであるときに空燃比がAFとなったとき、理論空燃比AFstoichを得るために必要な空気質量に対する過剰な空気質量はGa−Gastoichであるから、酸素の質量をMassO2とおくと、下記(49)式が得られ、この(49)式から上記(48)式が得られる。
次に、CPUはステップ740に進み、同ステップ740に記述した式(上記(30)式を参照。)に従ってCgout,O2(k+1)を求める。vgの値はエアフローメータ13が検出した吸入空気流量とする。このように、ステップ740では、対象としているブロックIから流出する酸素濃度Cgout,O2を新たに算出する。次いで、CPUはステップ745に進み、同ステップ745に記述した式(上記(28)式を参照。)に従ってCw,O2(k+1)を求める。即ち、CPUは、ステップ745にて対象としているブロックIのコート層の酸素濃度Cw,O2を新たに算出し、ステップ795を経由して図8に示したステップ800に進む。このように、図7により示したプログラムは、特定領域Iにおける排ガス相の酸素濃度推定手段、及びコート層の酸素濃度推定手段を構成している。
Next, the CPU proceeds to step 740 and obtains Cgout, O2 (k + 1) according to the equation described in step 740 (see the above equation (30)). The value of vg is the intake air flow rate detected by the
次に、CPUはステップ800からステップ805に進んで、前回の本プログラムの演算時において後述するステップ840にて算出されたコート層の一酸化炭素Cw,CO(k+1)を今回のコート層の一酸化炭素濃度Cw,COの値であるCw,CO(k)に設定する。
Next, the CPU proceeds from
次に、CPUは、ステップ810にて触媒19の温度Tempと同触媒19の劣化指標値REKKAと図12に示したようなマップMapkrelCOとから係数krel,CO(k)を決定し、続く、ステップ815にて同ステップ815内に記載した式(上記(38)式を参照。)に従って消費速度定数R*reduc,CO(k)を求め、その後、ステップ820にて見かけの拡散速度RD,CO(k)を触媒温度TempとマップMapRDCOとから決定する。
Next, in
続いて、CPUはステップ825にて一酸化炭素の反応律速因子SPreduc,COを同ステップ825内に記載した式(上記(33)式を参照。)により求め、ステップ830にて、このプログラムが対象としているブロックIよりも前の(即ち、上流の)ブロックI−1から流出する一酸化炭素濃度Cgout,CO(k)を、同対象としているブロックに流入する一酸化炭素濃度Cgin,CO(k)として取り込む。
Subsequently, the CPU obtains the carbon monoxide reaction rate-determining factor SPreduc, CO in
このとき、対象としているブロックIが触媒19の最も上流のブロックである場合、前のブロックI−1は存在しないから、ステップ830における前のブロックのCgout,CO(k)は、同触媒19に流入するガスの一酸化炭素濃度Cgin,COであり、同触媒19に流入するガスの空燃比A/Fと一酸化炭素濃度Cgin,COとの関係を規定した図13に示したマップに基づいて求められる。
At this time, if the target block I is the most upstream block of the
ここでは、触媒に流入するガスの空燃比A/Fと一酸化炭素濃度Cgin,COとの関係を予め実験により求めておき、この関係と同触媒に流入する実際の空燃比A/Fとに基づいて上記一酸化炭素濃度Cgin,COを求めているが、空燃比A/Fの他、点火時期、エンジン1の冷却水温、及び触媒に流入するガス量(例えば、エアフローメータ13により検出される吸入空気流量と等しい量)等のパラメータの一つ以上と一酸化炭素濃度Cgin,COとの関係を予め実験により求めておき、この関係と実際のパラメータとにも基づいて一酸化炭素濃度Cgin,COを更に精度良く求めてもよい。
Here, the relationship between the air-fuel ratio A / F of the gas flowing into the catalyst and the carbon monoxide concentration Cgin, CO is obtained in advance by experiment, and this relationship is compared with the actual air-fuel ratio A / F flowing into the catalyst. The carbon monoxide concentration Cgin, CO is calculated based on the above, but in addition to the air-fuel ratio A / F, the ignition timing, the cooling water temperature of the
次に、CPUはステップ835に進み、同ステップ835に記述した式(上記(34)式を参照。)に従ってCgout,CO(k+1)を求める。即ち、対象としているブロックIから流出する一酸化炭素濃度Cgout,COを新たに算出する。次いで、CPUはステップ840に進み、同ステップ840に記述した式(上記(32)式を参照。)に従ってCw,CO(k+1)を求める。即ち、CPUは、ステップ840にて対象としているブロックIのコート層の一酸化炭素濃度Cw,COを新たに算出し、ステップ895を経由して図9に示したステップ900に進む。このように、図8により示したプログラムは、特定領域Iにおける排ガス相の一酸化炭素濃度推定手段、及びコート層の一酸化炭素濃度推定手段を構成している。 Next, the CPU proceeds to step 835 and obtains Cgout, CO (k + 1) according to the equation described in step 835 (see the above equation (34)). That is, the carbon monoxide concentration Cgout, CO flowing out from the target block I is newly calculated. Next, the CPU proceeds to step 840 and obtains Cw, CO (k + 1) according to the equation described in step 840 (see the above equation (32)). That is, the CPU newly calculates the carbon monoxide concentration Cw, CO of the coating layer of the target block I in step 840, and proceeds to step 900 shown in FIG. As described above, the program shown in FIG. 8 constitutes the carbon monoxide concentration estimating means for the exhaust gas phase in the specific region I and the carbon monoxide concentration estimating means for the coat layer.
図9に示したプログラムは、炭化水素HCについての演算を行うプログラムであり、一酸化炭素COについての演算を行うための先に説明した図8のプログラムと同様なプログラムである。 The program shown in FIG. 9 is a program that performs an operation on hydrocarbon HC, and is a program similar to the program of FIG. 8 described above for performing an operation on carbon monoxide CO.
簡単に説明すると、CPUはステップ900からステップ905に進んで、前回の本プログラムの演算時において後述するステップ940にて算出されたコート層の一酸化炭素Cw,HC(k+1)を今回のコート層の一酸化炭素濃度Cw,HCの値であるCw,HC(k)に設定する。
Briefly, the CPU proceeds from
次に、CPUは、ステップ910にて触媒19の温度Tempと同触媒19の劣化指標値REKKAと図12に示したようなマップMapkrelHCとから係数krel,HC(k)を決定し、続く、ステップ915にて同ステップ915内に記載した式(上記(38)式を参照。)に従って消費速度定数R*reduc,HC(k)を求め、その後、ステップ920にて見かけの拡散速度RD,HC(k)を触媒温度TempとマップMapRDHCとから決定する。
Next, in
続いて、CPUはステップ925にて炭化水素の反応律速因子SPreduc,HCを同ステップ925内に記載した式(上記(33)式を参照。)により求め、ステップ930にて、このプログラムが対象としているブロックIよりも前の(即ち、上流の)ブロックI−1から流出する炭化水素濃度Cgout,HC(k)を、同対象としているブロックに流入する炭化水素濃度Cgin,HC(k)として取り込む。
Subsequently, the CPU obtains the hydrocarbon reaction rate-determining factor SPreduc, HC in
このとき、このプログラムが対象としているブロックIが触媒19の最も上流のブロックである場合、前のブロックI−1は存在せず、ステップ930におけるCgout,HC(k)は、同触媒19に流入するガスの炭化水素濃度Cgin,HCである。この場合、触媒19に流入するガスの空燃比A/Fと炭化水素濃度Cgin,HCとの関係は図14のグラフに示したようであり、エンジン1の排気温度と炭化水素濃度Cgin,HCとの関係は図15のグラフに示したようであるから、触媒19に流入するガスの空燃比A/F及びエンジン1の排気温度と炭化水素濃度Cgin,HCとの関係を予め実験により求めてマップとして記憶しておき、実際のガスの空燃比A/F及び実際のエンジン1の排気温度と同マップとから炭化水素濃度Cgin,HCを求める。
At this time, if the block I targeted by this program is the most upstream block of the
なお、上記炭化水素濃度Cgin,HCを求めるにあたり、空燃比A/F及びエンジン1の排気温度の他、点火時期、エンジン1の冷却水温、及び触媒に流入するガス量(例えば、エアフローメータ13により検出される吸入空気流量と等しい量)等のパラメータの一つ以上と炭化水素濃度Cgin,HCとの関係を予め実験により求めておき、この関係と実際のパラメータとにも基づいて上記炭化水素濃度Cgin,HCを更に精度良く求めてもよい。
In determining the hydrocarbon concentration Cgin, HC, in addition to the air-fuel ratio A / F and the exhaust temperature of the
次に、CPUはステップ935に進み、同ステップ935に記述した式(上記(34)式を参照。)に従ってCgout,HC(k+1)を求める。即ち、対象としているブロックIから流出する一酸化炭素濃度Cgout,HCを新たに算出する。次いで、CPUはステップ940に進み、同ステップ940に記述した式(上記(32)式を参照。)に従ってCw,HC(k+1)を求める。即ち、CPUは、ステップ940にて対象としているブロックIのコート層の一酸化炭素濃度Cw,HCを新たに算出し、ステップ995を経由して図10に示したステップ1000に進む。このように、図9により示したプログラムは、特定領域Iにおける排ガス相の炭化水素濃度推定手段、及びコート層の炭化水素濃度推定手段を構成している。 Next, the CPU proceeds to step 935 and obtains Cgout, HC (k + 1) according to the equation described in step 935 (see the above equation (34)). That is, the carbon monoxide concentration Cgout, HC flowing out from the target block I is newly calculated. Next, the CPU proceeds to step 940 and obtains Cw, HC (k + 1) according to the equation described in step 940 (see the above equation (32)). That is, the CPU newly calculates the carbon monoxide concentration Cw, HC of the coating layer of the target block I in step 940, and proceeds to step 1000 shown in FIG. As described above, the program shown in FIG. 9 constitutes the hydrocarbon concentration estimation means for the exhaust gas phase in the specific region I and the hydrocarbon concentration estimation means for the coat layer.
図10に示したプログラムは、窒素酸化物(ここでは、窒素酸化物の代表として一酸化窒素NOを選択した。)について、上記図7〜図9に示したプログラムに基く演算と同様の演算を行うためのプログラムである。 The program shown in FIG. 10 performs an operation similar to the operation based on the program shown in FIGS. 7 to 9 above for nitrogen oxide (here, nitric oxide NO is selected as a representative of nitrogen oxide). It is a program to do.
簡単に説明すると、CPUはステップ1000からステップ1005に進んで、前回の本プログラムの演算時において後述するステップ1040にて算出されたコート層の一酸化窒素Cw,NO(k+1)を今回のコート層の一酸化窒素濃度Cw,NOの値であるCw,NO(k+1)に設定する。
Briefly, the CPU proceeds from
次に、CPUは、ステップ1010にて触媒19の温度Tempと同触媒19の劣化指標値REKKAと図12に示したようなマップMapkstorNOとから酸素吸蔵速度係数kstor,NO(k)を決定し、続く、ステップ1015にて同ステップ1015内に記載した式(上記(36)式を参照。)に従って消費速度定数R*stor,NO(k)を求め、その後、ステップ1020にて見かけの拡散速度RD,NO(k)を触媒温度TempとマップMapRDNOとから決定する。
Next, in
続いて、CPUはステップ1025にて一酸化窒素の反応律速因子SPstor,NOを同ステップ1025内に記載した式(上記(29)式を参照。)により求め、ステップ1030にて、このプログラムが対象としているブロックIよりも前の(即ち、上流の)ブロックI−1から流出する一酸化窒素濃度Cgout,NO(k)を、同対象としているブロックに流入する一酸化窒素濃度Cgin,NO(k)として取り込む。
Subsequently, the CPU obtains the reaction rate limiting factor SPstor, NO of nitric oxide in
このとき、このプログラムが対象としているブロックIが触媒19の最も上流のブロックである場合、前のブロックI−1は存在せず、ステップ1030における前のブロックのCgout,NO(k)は、同触媒19に流入するガスの一酸化窒素濃度Cgin,NOである。この場合、触媒19に流入するガスの空燃比A/Fと一酸化窒素濃度Cgin,NOとの関係は図16のグラフに示したようであり、エンジン1の一気筒の一吸気行程あたりの吸入空気流量Mc(筒内吸入空気量)と一酸化窒素濃度Cgin,NOとの関係は図17のグラフに示したようであるから、触媒19に流入するガスの空燃比A/F及び筒内吸入空気量と一酸化窒素濃度Cgin,NOとの関係を予め実験により求めてマップとして記憶しておき、実際のガスの空燃比A/F及び実際の筒内吸入空気量と同マップとから一酸化窒素濃度Cgin,NOを求める。
At this time, if the block I targeted by this program is the most upstream block of the
なお、上記一酸化窒素濃度Cgin,NOを求めるにあたり、触媒に流入するガスの空燃比A/F及び筒内吸入空気量の他、点火時期、エンジン1の冷却水温、及び触媒に流入するガス量(例えば、エアフローメータ13により検出される吸入空気流量と等しい量)等のパラメータの一つ以上と一酸化窒素濃度Cgin,NOとの関係を予め実験により求めておき、この関係と実際のパラメータとにも基づいて上記一酸化窒素濃度Cgin,NOを更に精度良く求めてもよい。
In determining the nitrogen monoxide concentration Cgin, NO, in addition to the air-fuel ratio A / F of the gas flowing into the catalyst and the in-cylinder intake air amount, the ignition timing, the cooling water temperature of the
次に、CPUはステップ1035に進み、同ステップ1035に記述した式(上記(30)式を参照。)に従ってCgout,NO(k+1)を求める。即ち、対象としているブロックIから流出する一酸化窒素濃度Cgout,NOを新たに算出する。次いで、CPUはステップ1040に進み、同ステップ1040に記述した式(上記(28)式を参照。)に従ってCw,NO(k+1)を求める。即ち、CPUは、ステップ1040にて対象としているブロックIのコート層の一酸化窒素濃度Cw,NOを新たに算出し、ステップ1095を経由して図11に示したステップ1100に進む。このように、図10により示したプログラムは、特定領域Iにおける排ガス相の一酸化窒素濃度推定手段、及びコート層の一酸化窒素濃度推定手段を構成している。 Next, the CPU proceeds to step 1035 and obtains Cgout, NO (k + 1) according to the equation described in step 1035 (see equation (30) above). That is, the nitric oxide concentration Cgout, NO flowing out from the target block I is newly calculated. Next, the CPU proceeds to step 1040 to obtain Cw, NO (k + 1) according to the equation described in step 1040 (see the above equation (28)). That is, the CPU newly calculates the nitric oxide concentration Cw, NO of the coating layer of the target block I in step 1040, and proceeds to step 1100 shown in FIG. As described above, the program shown in FIG. 10 constitutes the nitric oxide concentration estimating means for the exhaust gas phase in the specific region I and the nitric oxide concentration estimating means for the coat layer.
図11に示したプログラムは、酸素吸蔵密度(酸素吸蔵濃度)Ostを算出するためのプログラムであり、具体的に述べると、CPUはステップ1105にて上記(43)式に基く同ステップ1105内に記述した式により係数Pを求めるとともに、続くステップ1120にて上記(44)式に基く同ステップ1120に記述した式により係数Qを求める。次いで、CPUはステップ1125にて上記(42)式に基く同ステップ1125に記述した式により酸素吸蔵密度Ost(k+1)を求め、ステップ1195にて図7〜図11に示した本プログラムを一旦終了する。なお、ステップ1125の次に、上記(45)式に基づいてこのブロックの酸素吸蔵量OSA,Iを求めてから、ステップ1195に進むように構成してもよい。このように、図11に示したプログラムは、酸素吸蔵密度算出手段、及び酸素吸蔵量算出手段を構成している。
The program shown in FIG. 11 is a program for calculating the oxygen storage density (oxygen storage concentration) Ost. Specifically, the CPU executes
以上、説明したように、対象としているブロックIの化学種iの濃度Cgout,iが求められ、上述した風上法を用いることで、順次隣接するブロックIの化学種iの濃度Cgout,iが求められて行く。同様に、各ブロックI毎に酸素吸蔵密度Ostが求められ、上述した(45)式を使用することで、同ブロックIの酸素吸蔵量OSAが求められる。また、このブロックIの酸素吸蔵量OSAを触媒の入口から任意のブロックKまで積算すれば、同ブロックKまでの触媒内の積算酸素吸蔵量OSA,Kが求められ、ブロックKが触媒の出口のブロックであれば、同触媒の酸素吸蔵量OSAallが求められる。 As described above, the concentration Cgout, i of the chemical species i of the target block I is obtained, and by using the above-described upwind method, the concentration Cgout, i of the chemical species i of the adjacent block I is successively determined. Go for it. Similarly, the oxygen storage density Ost is obtained for each block I, and the oxygen storage amount OSA of the block I is obtained by using the above-described equation (45). Further, if the oxygen storage amount OSA of this block I is integrated from the inlet of the catalyst to an arbitrary block K, the integrated oxygen storage amount OSA, K in the catalyst up to the block K is obtained, and the block K is the outlet of the catalyst. If it is a block, the oxygen storage amount OSAall of the catalyst is obtained.
次に、上述した空燃比制御装置を用いた制御の実施形態について、順次説明する。 Next, embodiments of control using the above-described air-fuel ratio control apparatus will be sequentially described.
まず、第一実施形態について、その制御のフローチャートを示した図18を参照しながら説明する。本実施形態では、排ガス中の酸素量(実際には、触媒19から排出される酸素の濃度)と一酸化炭素量(実際には触媒19から排出される一酸化炭素の濃度)を上述した推定値として求める。また、酸素量は空燃比制御において酸素の過不足量として扱う。つまり、酸素が過剰に存在する場合に酸素量は正の値、酸素が不足している場合に酸素量は負の値とする。更に、本実施形態では、触媒19全体を一つの領域(特定領域)として考えている。
First, the first embodiment will be described with reference to FIG. 18 showing a flowchart of the control. In the present embodiment, the amount of oxygen in exhaust gas (actually the concentration of oxygen discharged from the catalyst 19) and the amount of carbon monoxide (actually the concentration of carbon monoxide discharged from the catalyst 19) are estimated as described above. Calculate as a value. Further, the oxygen amount is treated as an excess / deficiency amount of oxygen in the air-fuel ratio control. That is, the oxygen amount is a positive value when oxygen is excessive, and the oxygen amount is a negative value when oxygen is insufficient. Further, in the present embodiment, the
触媒19での排ガスの浄化反応においては、排ガス中の浄化すべき成分を酸化させるために(即ち、未燃燃料である炭化水素HCや一酸化炭素COを酸化させるために)酸素が消費される。このため、触媒19から流出する排ガス中にこれらの酸化すべき成分が含まれているときは、酸素量が不足しているとみなす。即ち、推定値は負の値をとる。なお、ここでは、未燃成分を一酸化炭素COで代表させる。
In the exhaust gas purification reaction at the
これに対し、排ガス中の浄化すべき成分(窒素酸化物NOx)が還元されることにより生じる酸素、及び流入する排ガス中に含まれる酸素が、触媒19の酸素吸蔵作用によって同触媒19に吸蔵しきれない場合、酸素成分が同触媒19から流出する。このときは、酸素量が過剰であるとみなす。即ち、推定値は正の値をとる。なお、ここでは、酸素吸蔵作用に係る成分を酸素O2で代表させる。
On the other hand, oxygen generated by the reduction of the component to be purified (nitrogen oxide NOx) in the exhaust gas and oxygen contained in the inflowing exhaust gas are stored in the
本空燃比制御装置は、この酸素に関する推定値を用いて空燃比制御を行う。具体的に述べると、同装置のCPUは、図18のフローチャートに示したように、ステップ10にて上流側空燃比センサ25が検出する触媒19に流入する排ガスの空燃比(排気空燃比)を取得し、次のステップ30にて前回の(本プログラムの前回の実行時における後述するステップ55にて求めた)酸素吸蔵密度Ostを用いてRosc(即ちR*stor,O2、R*reudc,CO)及びRD(RD,O2、RD,CO)を算出する。
The air-fuel ratio control apparatus performs air-fuel ratio control using the estimated value related to oxygen. Specifically, as shown in the flowchart of FIG. 18, the CPU of the apparatus determines the air-fuel ratio (exhaust air-fuel ratio) of the exhaust gas flowing into the
このとき、触媒19に流入する排ガスの空燃比が理論空燃比よりリーンの場合、RDの算出には既に説明した(8)式に基く以下の(50)式を用い、Roscの算出には、既に説明した(36)式に基く以下の(51)式を用いる。
また、触媒19に流入する排ガスの空燃比が理論空燃比よりリッチの場合、RDの算出には既に説明した(8)式に基く以下の(52)式を用い、Roscの算出には、既に説明した(38)式に基く以下の(53)式を用いる。
Roscについては、例えば温度についての関数となる。ここでは、例えば、hDが温度の関数であり、ksotr,O2、krel,COは図12に示したように、温度と触媒劣化度の関数である。温度は温度センサ21によって検出され、検出された温度に基づいてECU18内のROMに格納された関数式やマップからRoscが算出される。
For Rosc, for example, it is a function of temperature. Here, for example, h D is a function of temperature, and ksotr, O2, krel, and CO are functions of temperature and the degree of catalyst deterioration, as shown in FIG. The temperature is detected by the
その後、CPUはステップ40にて、触媒19に流入する排ガスの空燃比が理論空燃比よりリーンの場合、上述した(30)式を用いてCgout,O2(=CgoutO2)、触媒19に流入する排ガスの空燃比が理論空燃比よりリッチの場合、上述した(34)式を用いてCgout,COを算出する。そして、Cgout,COが求められたときには、Cgout,O2=−2・Cgout,COとして、最終的なCgout,O2を求める。
After that, in
また、ステップ40でCgout,O2(Cgout,COで置換する前のCgout,O2)及びCgout,COを算出するには、Cgin,O2及びCgin,COがそれぞれ必要となる。ステップ40では、Cgin,O2を上述した(48)式から求めている。このとき、排気空燃比が理論空燃比よりもリッチであり、(48)式に従って計算したCgin,O2が正の値であればその値をそのままCin,O2として採用し、(48)式に従って計算したCgin,O2が負の値であれば、その値の絶対値の1/2をCgin,COをとして採用する。
Further, in order to calculate Cgout, O2 (Cgout, O2 before substitution with Cgout, CO) and Cgout, CO in
次いで、CPUはステップ50にてコート層の各化学種濃度Cw(Cw,O2、Cw,CO)を(28)式及び(32)式に基づいて求め、ステップ55にて(42)式〜(44)式によりOst(今回値)を算出するとともに、ステップ60にて、空燃比制御用のフィードバック補正量(F/B補正量)を、(Cgout,O2−O2ref)・Gなどから求める。ここで、O2refは制御目標(目標状態)、Gは制御ゲインである。図19は、このCgout,O2とF/B補正量との関係を示したグラフである。本実施形態では、制御目標は所定の範囲として定められている。決定されたF/B補正量は、空燃比を決定する燃料噴射量の補正係数の一つとして空燃比制御に反映される。
Next, the CPU obtains the chemical species concentration Cw (Cw, O2, Cw, CO) of the coat layer based on the equations (28) and (32) in
次に、第二実施形態について、その制御のフローチャートを示した図20を参照しながら説明する。本実施形態では、上述した推定値として、排ガス中のリッチ成分とリーン成分を用いている。また、本実施形態では、触媒19全体を一つの領域(特定領域)として考えている。
Next, a second embodiment will be described with reference to FIG. 20 showing a flowchart of the control. In the present embodiment, rich components and lean components in the exhaust gas are used as the estimated values described above. In the present embodiment, the
リッチ成分とは、排気空燃比がリッチのときに排ガス中の含有量が増加する成分を総合したもので、触媒19から流出する排ガスの状態を示す代表値の一つである。具体的には、排ガスに含まれるのCOやHCの量を一つの指標としてまとめた代表値である。一方、リーン成分とは、排気空燃比がリーンのときに排ガス中の含有量が増加する成分を総合したもので、触媒19から流出する排ガスの状態を示す代表値の一つである。具体的には、排ガスに含まれるNOxやO2の量を一つの指標としてまとめた代表値である。本実施形態では、上述した二つの推定値を用いて空燃比制御を行うので、Cgoutに関しては、リッチ成分に関するものをCgoutRとおく。つまり、
まず、図20のフローチャートに示されるように、CPUはステップ70にてエンジン1の燃焼後の排ガス(触媒19に流入する排ガス)中のリッチ成分CginRとリーン成分CginLとを推定する。かかるリッチ成分CginRとリーン成分CginLの推定においては、触媒19での酸化・還元反応(三元反応)及び酸素吸蔵放出反応を考慮する。三元反応に関連する化学式を(54)式に示す。また、酸素吸蔵放出反応に関連する化学式を(55)式に示す。なお、この酸素吸蔵放出反応における各反応の反応速度係数を、KO2(=Kstor,O2)、KNO(=Kstor,NO),KCO(=Krel,CO),KHC(=Krel,HC)とする。(55)式での酸素の吸蔵・放出においては、酸素吸蔵・放出機能を生じせしめる担持金属をセリウム(Ce)として記載してある。
次いで、CPUは、ステップ100に進んで上述した第一実施形態と同様に、Rosc(即ちR*stor,O2、R*stor,NO、R*reduc,CO、R*reudc,HC)及びRD(即ちRD,O2、RD,NO、RD,CO、RD,HC)を算出し、続くステップ110にて、上述した(34)式を用いてCgoutR(=Cgout,reduc,CO+Cgout,reduc,HC)を算出し、(30)式を用いてCgoutL(=Cgout,stor,O2+Cgout,stor,NO)を算出する。
Next, the CPU proceeds to step 100 and, similarly to the first embodiment described above, Rosc (ie, R * stor, O2, R * stor, NO, R * reduc, CO, R * reudc, HC) and R D (That is, R D , O 2, R D , NO, R D , CO, R D , HC) is calculated, and in the
さらに、CPUはステップ120にてコート層の各化学種濃度Cw(即ちCw,O2、Cw,NO、Cw,CO、Cw,HC)を(28)式及び(32)式を用いて求め、ステップ125にて(42)式〜(44)式によりOst(今回値)を算出するとともに、ステップ130にて空燃比制御用のフィードバック補正量(F/B補正量)を、(CgoutR−Ref)×GR+(CgoutL−Ref)×GLから求める。ここで、(CgoutR−Ref)×GRがリッチ成分に対応する部分で、(CgoutL−Ref)×GLがリーン成分に対応する部分である。Refは制御目標(目標状態)、GR,GLは各成分毎に設定される制御ゲインである。
Further, the CPU obtains the chemical species concentration Cw (ie, Cw, O2, Cw, NO, Cw, CO, Cw, HC) of the coat layer in
ここで、CgoutRとCgoutLとが等しくなるように空燃比制御が行われる。即ち、CgoutR及びCgoutLが何れも目標値Refとなるように、各成分に対応するF/B補正量(CgoutR−Ref)×GR,(CgoutL−Ref)×GLが決定される。また、本実施形態では、制御目標は所定の範囲として定められている。決定された補正量は、空燃比を決定する燃料噴射量の補正係数の一つとして空燃比制御に反映される。 Here, air-fuel ratio control is performed so that CgoutR and CgoutL are equal. That is, F / B correction amounts (CgoutR−Ref) × GR and (CgoutL−Ref) × GL corresponding to each component are determined so that both CgoutR and CgoutL become the target value Ref. In the present embodiment, the control target is defined as a predetermined range. The determined correction amount is reflected in the air-fuel ratio control as one of the fuel injection amount correction coefficients for determining the air-fuel ratio.
なお、より具体的には、下記(56)式によりフィードバック補正量を計算することになる。
次に、第三実施形態について、その制御のフローチャートを示した図21を参照しながら説明する。本実施形態では、上述した推定値として、排ガス中の酸素排出量、一酸化炭素排出量及び一酸化窒素排出量(実際には、それらの各濃度)を用いている。また、本実施形態では、触媒19全体を一つの領域(特定領域)として考えている。本実施形態では、上述した三つの推定値を用いて空燃比制御を行うので、Cgoutに関しては、酸素排出量に関するものをCgoutO2(Cgout,O2)とし、一酸化炭素排出量に関するものをCgoutCO(Cgout,CO)とし、一酸化窒素排出量に関するものをCgoutNO(Cgout,NO)として説明する。
Next, a third embodiment will be described with reference to FIG. 21 showing a flowchart of the control. In this embodiment, the oxygen emission amount, the carbon monoxide emission amount, and the nitrogen monoxide emission amount (actually their respective concentrations) in the exhaust gas are used as the estimated values described above. In the present embodiment, the
まず、CPUは、図21のフローチャートに示されるように、ステップ140にて、エンジン1の燃焼後の排ガス(触媒19に流入する排ガス)中のリッチ成分CginR(Cgin,CO)とリーン成分CginL(Cgin,O2, Cgin,NO)とを前述したように推定し、ステップ170にて、Rosc(即ちR*stor,O2、R*stor,NO、R*reduc,CO)及びRD(即ちRD,O2、RD,NO、RD,CO)を算出する。その後、CPUはステップ180にて上述した(34)式を用いてCgoutR(Cgout,reduc,CO)を算出し、(30)式を用いてCgoutL(Cgout,stor,O2=CgoutO2とCgout,stor,NO=CgoutNO)を算出する。
First, as shown in the flowchart of FIG. 21, in
さらに、CPUはステップ190にてコート層における各化学種濃度Cw(即ちCw,O2、Cw,NO、Cw,CO)を(28)式及び(32)式を用いて求め、ステップ195にて(42)式〜(44)式によりOst(今回値)を算出する。そして、CPUは、ステップ200にて空燃比制御にフィードバックするためのフィードバック補正量(F/B補正量)を、(CgoutO2−O2ref)×G02+(CgoutNO−N0ref)×GNO+(CgoutCO−COref)×GCOから求める。ここで、(CgoutO2−O2ref)×G02が酸素成分に対応する部分で、(CgoutNO−NOref)×GNOが一酸化窒素成分に対応する部分で、(CgoutCO−COref)×GCOが一酸化炭素成分に対応する部分である。O2ref,NOref,COrefは各成分毎の制御目標、G02,GNO,GCOは各成分毎に設定される制御ゲインである。決定された補正量は、空燃比を決定する燃料噴射量の補正係数の一つとして空燃比制御に反映される。なお、本実施形態においても、制御目標を所定の範囲として定めていてもよい。
Further, the CPU obtains the chemical species concentration Cw (that is, Cw, O2, Cw, NO, Cw, CO) in the coat layer in
次に、第四実施形態について説明する。本実施形態では、第一実施形態と同様に、上述した推定値として排ガス中の酸素量(触媒19から排出される酸素の濃度)を用いるが、その推定値は所定時間後の推定値(現時点から所定時間先の時点の推定値)である。すなわち、本実施形態においては、いわゆる「先読み」による制御を行う。これに伴って、推定値の推定手法も、第一実施形態とは異なる。ここでの推定値(酸素量)も過不足量として把握され、触媒19全体が一つの領域(特定領域)として考えられている。
Next, a fourth embodiment will be described. In the present embodiment, as in the first embodiment, the amount of oxygen in the exhaust gas (the concentration of oxygen discharged from the catalyst 19) is used as the estimated value described above. The estimated value is an estimated value after a predetermined time (current time). (Estimated value at a predetermined time ahead). That is, in the present embodiment, control based on so-called “look ahead” is performed. Accordingly, the estimation value estimation method is also different from the first embodiment. The estimated value (oxygen amount) here is also grasped as an excess / deficiency amount, and the
推定値(排出酸素過不足量)の推定には上述した第一実施形態と同様の推定モデル(推定方法)が用いられる。本実施形態ではさらに、吸入空気量や燃焼に寄与する燃料量などもモデルを用いて推定する。吸入空気量に関するモデルを、ここでは空気モデルと言うこととし、燃料量に関するモデルを、ここでは燃料挙動モデルと言うこととする。 For estimation of the estimated value (exhaust oxygen excess / deficiency), the same estimation model (estimation method) as in the first embodiment described above is used. In the present embodiment, the intake air amount and the fuel amount contributing to combustion are also estimated using a model. Here, the model relating to the intake air amount is referred to as an air model, and the model relating to the fuel amount is referred to herein as a fuel behavior model.
第四実施形態における制御のフローチャートを図22に示す。まず、図22のフローチャートに示されるように、CPUは上述した空気モデルを用いて、所定時間後の吸入空気量を推定する(ステップ210)。この空気モデルについての詳しい説明はここでは省略するが、一例としては、吸入空気量を検出する吸気管内圧センサやエアフローメータ13の出力に加えて、スロットルバルブ9の開度の時間変化やその微積分値、エンジン回転数の時間的変化や微積分値などを用いて、所定時間後の筒内吸入空気量を推定する手法などがある。
FIG. 22 shows a flowchart of control in the fourth embodiment. First, as shown in the flowchart of FIG. 22, the CPU estimates the intake air amount after a predetermined time using the air model described above (step 210). Although a detailed description of this air model is omitted here, as an example, in addition to the output of the intake pipe internal pressure sensor for detecting the intake air amount and the output of the
以下、スロットルバルブ9の開度を用いた吸入空気流量の推定方法の一例について簡単に説明する。
(1)CPUは、アクセルポジションセンサ12から現時点のアクセル開度を検出する。
(2)CPUは、検出した現時点のアクセル開度に基づいて暫定的な目標スロットルバルブ開度を決定する。
(3)CPUは、決定された暫定的な目標スロットルバルブ開度を所定時間だけ遅延して実際の目標スロットルバルブ開度とする。換言すると、現時点から所定時間前の暫定的な目標スロットルバルブ開度が実際の目標スロットルバルブ開度として設定される。暫定的な目標スロットルバルブ開度は、それが実際の目標スロットルバルブ開度として使用されるまで時系列的にRAM内に格納される。
(4)CPUは、設定された実際の目標スロットルバルブ開度とスロットルバルブ9の開度が一致するようにスロットルモータ11を制御する。
(5)CPUは、現時点から次に吸気行程を迎える気筒の吸気弁が閉じる時点(吸気弁閉弁時)に至るまでの時間を推定し、その推定した時間と前記RAMに格納されている暫定的な目標スロットルバルブ開度とに基づいて、同気筒の吸気弁閉弁時における目標スロットルバルブ開度を算出し、同目標スロットルバルブ開度により同吸気弁閉弁時の推定スロットルバルブ開度を算出する。。
(6)CPUは、次に吸気行程を迎える気筒の吸気弁閉弁時の吸入空気流量(筒内吸入空気量)Mcを、前記推定スロットルバルブ開度とエンジン回転速度とマップとに基づいて推定する。
Hereinafter, an example of a method for estimating the intake air flow rate using the opening of the throttle valve 9 will be briefly described.
(1) The CPU detects the current accelerator opening from the
(2) The CPU determines a temporary target throttle valve opening based on the detected current accelerator opening.
(3) The CPU delays the determined provisional target throttle valve opening by a predetermined time to obtain the actual target throttle valve opening. In other words, the provisional target throttle valve opening a predetermined time before the current time is set as the actual target throttle valve opening. The provisional target throttle valve opening is stored in the RAM in time series until it is used as the actual target throttle valve opening.
(4) The CPU controls the
(5) The CPU estimates the time from the current time until the time when the intake valve of the cylinder that will reach the next intake stroke closes (when the intake valve is closed), and the estimated time and the provisional stored in the RAM The target throttle valve opening when the intake valve of the cylinder is closed is calculated based on the target throttle valve opening, and the estimated throttle valve opening when the intake valve is closed is calculated based on the target throttle valve opening. calculate. .
(6) The CPU estimates the intake air flow rate (in-cylinder intake air amount) Mc when the intake valve of the cylinder that reaches the next intake stroke is closed based on the estimated throttle valve opening, the engine speed, and the map. To do.
次いで、CPUは、ステップ210の推定結果を踏まえて、基本燃料噴射量(例えば、エンジン1に吸入される混合気の空燃比を理論空燃比とするために必要な噴射量)を算出する(ステップ220)。この基本燃料噴射量が各種補正係数により補正され、最終的な燃料噴射量が決定される。
Next, the CPU calculates a basic fuel injection amount (for example, an injection amount necessary for setting the air-fuel ratio of the air-fuel mixture sucked into the
さらに、この基本燃料噴射量(あるいは基本燃料噴射量に対して補正が加えられた燃料噴射量)と燃料挙動モデルとに基づいて排ガスの空燃比を推定する(ステップ230)。この燃料挙動モデルについての詳しい説明もここでは省略するが、一例としては、基本燃料噴射量自体に加えて、吸気ポートやシリンダ内壁等の部材への燃料の付着分や同内壁に付着している燃料の剥離分などをも考慮して、所定時間後の排気空燃比を推定する手法などがある。 Further, the air-fuel ratio of the exhaust gas is estimated based on the basic fuel injection amount (or the fuel injection amount corrected for the basic fuel injection amount) and the fuel behavior model (step 230). A detailed description of this fuel behavior model is also omitted here, but as an example, in addition to the basic fuel injection amount itself, the fuel adheres to the members such as the intake port and the cylinder inner wall, and adheres to the inner wall. There is a method for estimating the exhaust air-fuel ratio after a predetermined time in consideration of the amount of fuel separation.
例えば、この燃料モデルによれば、fi(k)だけの燃料量をインジェクタ5から噴射した後の燃料付着量fw(k+1)は下記(57)式により求められる。
fw(k+1)=R・fi(k)+P・fw(k) …(57)
For example, according to this fuel model, the fuel adhesion amount fw (k + 1) after injecting the fuel amount of fi (k) from the injector 5 is obtained by the following equation (57).
fw (k + 1) = R ・ fi (k) + P ・ fw (k)… (57)
上記(57)式において、fw(k)はfi(k)だけの燃料量を噴射する前の吸気ポート等への燃料付着量、Pは吸気ポート等に既に付着していた燃料のうち一吸気行程を経た後に同吸気ポート等に付着したまま残留している燃料の割合(残留率)、Rは噴射された燃料のうち吸気ポート等へ直接付着する燃料の割合(付着率)である。 In the above equation (57), fw (k) is the amount of fuel adhering to the intake port etc. before injecting the fuel amount of fi (k), and P is one intake of the fuel already adhering to the intake port etc. The ratio (residual ratio) of the fuel remaining after adhering to the intake port after passing the stroke, and R is the ratio (adhesion ratio) of the fuel directly adhering to the intake port etc. of the injected fuel.
一方、今回の燃料噴射量fi(k)のうち燃焼室(気筒)内に吸入される燃料量は(1−R)・fi(k)となり、既に付着している燃料の量(燃料付着量)fw(k)のうち燃焼室(気筒)内に吸入される燃料量は(1−P)・fw(k)となる。即ち、今回の燃料噴射により、気筒内には下記(58)式に示した燃料量fc(k)が吸入されることになる。
fc(k)=(1−R)・fi(k)+(1−P)・fw(k) …(58)
On the other hand, the amount of fuel sucked into the combustion chamber (cylinder) in this fuel injection amount fi (k) is (1-R) · fi (k), and the amount of fuel already attached (fuel attachment amount) ) Of fw (k), the amount of fuel sucked into the combustion chamber (cylinder) is (1-P) · fw (k). That is, with this fuel injection, the fuel amount fc (k) shown in the following equation (58) is sucked into the cylinder.
fc (k) = (1-R) ・ fi (k) + (1-P) ・ fw (k) (58)
従って、実際には、基本噴射量の燃料をインジェクタ5から噴射し、上記(58)式のfi(k)を基本噴射量と等しいとおいて実際に気筒に吸入される燃料量fc(k)を求め、前記推定した筒内吸入空気量Mcをこの燃料量fc(k)で除することにより、排気空燃比を予測する。また、次回の計算に備え、(57)式のfi(k)を基本噴射量と等しいとおいて燃料付着量fw(k+1)を求めておく。 Therefore, in actuality, the fuel of the basic injection amount is injected from the injector 5, and the fuel amount fc (k) actually sucked into the cylinder is set by assuming that fi (k) in the above equation (58) is equal to the basic injection amount. The exhaust air / fuel ratio is predicted by dividing the estimated in-cylinder intake air amount Mc by the fuel amount fc (k). Further, in preparation for the next calculation, the fuel adhesion amount fw (k + 1) is obtained assuming that fi (k) in the equation (57) is equal to the basic injection amount.
次に、CPUは、ステップ235に進み、前記ステップ230にて推定された排気空燃比を上記(48)式のAFとおき、同(48)式からCginO2(=Cgin,O2)を推定するとともに、前記ステップ230にて推定された排気空燃比、及び触媒に流入するガスの空燃比A/Fと一酸化炭素濃度CginCOとの関係を規定した図13に示したマップに基づいて一酸化炭素濃度CginCOを推定する。その後、CPUは、第一実施形態で述べた手法と同様にしてステップ240にてCgoutO2を推定する。このCgoutO2は、前記排気空燃比が現時点から所定時間後の空燃比を予測(先読み)したものであり、従って、CginO2及びCginCOもそれぞれ現時点から所定時間後のCginO2及びCginCOとなるから、現時点から所定時間後のCgoutO2を予測した値となっている。
Next, the CPU proceeds to step 235, sets the exhaust air-fuel ratio estimated in
次いで、CPUはステップ250にて前記CgoutO2が0より大きいか否かを判定し、同ステップにて「Yes」と判定したときはリッチ制御(ステップ260)、「No」と判定したときはリーン制御(ステップ270)を実行する。即ち、本実施形態におけるCgoutO2の目標状態(目標値)は0(ゼロ)であり、同CgoutO2が0となるように空燃比制御が実行される。
Next, the CPU determines in
リッチ制御は空燃比(機関の空燃比)をリッチ側に補正する空燃比制御(噴射量を増量する制御)で、リーン制御は空燃比をリーン側に補正する制御(噴射量を減量する制御)である。本実施形態では、その補正量はリッチ制御・リーン制御の各制御毎に一定である。なお、推定したCgoutO2とその目標状態との乖離度に応じて補正量を変更してもよい。また、ステップ260では空燃比をリッチ側に補正する補正係数を算出し、ステップ270では空燃比をリーン側に補正する補正係数を算出し、これらの補正係数を基本噴射量に乗ずるなどの処理を行って燃料噴射制御に反映させるようにしてもよい。
Rich control is air-fuel ratio control (control that increases the injection amount) that corrects the air-fuel ratio (engine air-fuel ratio) to the rich side, and lean control is control that corrects the air-fuel ratio to the lean side (control that reduces the injection amount). It is. In the present embodiment, the correction amount is constant for each control of rich control and lean control. Note that the correction amount may be changed according to the degree of deviation between the estimated CgoutO2 and the target state. In
次に、第五実施形態について説明する。本実施形態でも、上述した第四実施形態と同様に、いわゆる「先読み」による制御を行う。また、本実施形態においては、排ガス中の複数の排出成分についての推定値が算出される。ここでは、これら複数の推定値(具体的にはCgout,O2やCgout,COなど)をCgout**として記載することとする。同様に、CginについてはCgin**として記載する。また、ここでも触媒19全体は、一つの領域(特定領域)として扱われている。
Next, a fifth embodiment will be described. Also in the present embodiment, control by so-called “prefetching” is performed as in the fourth embodiment described above. Further, in the present embodiment, estimated values for a plurality of emission components in the exhaust gas are calculated. Here, these multiple estimated values (specifically, Cgout, O2, Cgout, CO, etc.) are described as Cgout **. Similarly, Cgin is described as Cgin **. Also here, the
第五実施形態における制御のフローチャートを図23に示す。まず、図23のフローチャートに示されるように、CPUは、上述した空気モデルを用いて、所定時間後の吸入空気量を推定する(ステップ280)。次いで、CPUは、上述した燃料挙動モデルを用いてシリンダ3内に実際に供給されるであろう燃料量を推定し(ステップ290)、ステップ280及びステップ290の結果から空燃比を算出する(ステップ300)。
A flowchart of control in the fifth embodiment is shown in FIG. First, as shown in the flowchart of FIG. 23, the CPU estimates the intake air amount after a predetermined time using the air model described above (step 280). Next, the CPU estimates the amount of fuel that will actually be supplied into the cylinder 3 using the fuel behavior model described above (step 290), and calculates the air-fuel ratio from the results of
そして、CPUは、算出した空燃比に基づいて触媒19に流入する各成分毎のCgin**を算出する(ステップ310)。さらに、触媒19内の反応を考慮して、各成分毎のCgout**と触媒19の酸素吸蔵量OSAとを推定する(ステップ320)。Cgin**とCgout**の推定については、上述した実施形態の手法による。また、OSAは、上述した(42)〜(45)式に基づいて推定する。
Then, the CPU calculates Cgin ** for each component flowing into the
次に、CPUは、排気浄化に不利な状況となりそうかどうかを、スロットル駆動要求やその時点でのOSAなどの運転条件から予測する。スロットル駆動要求の有無は、運転者によるアクセルベグルの操作量(アクセル開度)の時間変化(アクセル開度の時間微分値)に基づいて判断され得る。また、排気浄化に不利な状況であるか否かは、たとえば、スロットル開度が急激に変化する過渡状態時であるか否かに基づいて判定され得る。過渡状態時には、触媒19の酸素吸蔵状態が不安定になりやすく、空燃比制御による補正が間に合わずに排気浄化率が悪化する可能性が高い。このようなとき、CPUは、スロットルバルブ9の開閉動作を電子制御によって遅延させ(開閉動作速度を小さくし)、これにより浄化すべき成分の触媒19下流側への流出を抑止しつつ、スロットル動作の遅延による要求トルク未達分を点火時期制御によって補う。
Next, the CPU predicts whether or not the situation is likely to be disadvantageous for exhaust purification based on the throttle drive request and the operating conditions such as OSA at that time. The presence / absence of a throttle drive request can be determined based on the time change (time differential value of the accelerator opening) of the operation amount (accelerator opening) of the accelerator beggle by the driver. In addition, whether or not the situation is disadvantageous for exhaust purification can be determined based on, for example, whether or not the engine is in a transient state in which the throttle opening changes rapidly. During the transient state, the oxygen storage state of the
具体的に述べると、CPUは、スロットル駆動要求があるか否かを判定し(ステップ330)、このようなスロットル駆動要求がなければ、排気浄化に不利な状況にはなりにくいと判定して、図23のフローチャートに示される制御を終了する。一方、スロットル駆動供給がある場合、CPUは、その時点でのOSA、算出したCgout**、スロットル開度予測値、及びアクセル開度の時間変化量などの運転状態に基づいて、排気浄化率が悪化することが予想されるか否かを判定する(ステップ340)。そして、CPUは排気浄化率が悪化することが予想されないと判定すると、駆動要求通りにスロットルバルブ9を駆動する(ステップ360)。 Specifically, the CPU determines whether or not there is a throttle drive request (step 330). If there is no such throttle drive request, the CPU determines that it is unlikely to be a disadvantageous condition for exhaust purification, The control shown in the flowchart of FIG. On the other hand, when there is a throttle drive supply, the CPU determines the exhaust gas purification rate based on operating conditions such as the OSA at that time, the calculated Cgout **, the predicted throttle opening, and the amount of time change in the accelerator opening. It is determined whether or not deterioration is expected (step 340). If the CPU determines that the exhaust gas purification rate is not expected to deteriorate, it drives the throttle valve 9 as requested (step 360).
―方、CPUは、排気浄化率の悪化が予測されるとき、スロットル駆動動作の遅延量および点火時期によるトルク補正量とを決定するとともに、上述した推定値Cgout**を目標状態に一致させるための空燃比制御を行い(ステップ350)、スロットル駆動を行う(ステップ360)。スロットル駆動動作の遅延量は、例えば、スロットルバルブ開度の急激な変化を抑制するため、上述した実際の目標スロットルバルブ開度に対して一次遅れ処理を施した開度を実際のスロットルバルブ開度として設定し直すために必要な量である。点火時期の補正量は、この実際のスロットルバルブ開度を設定し直したことに伴うエンジントルクの低下分を補うための補正量(進角量)である。また、推定値Cgout**を目標状態とするための空燃比制御は、上述した他の実施形態と同様である。 -On the other hand, when the exhaust gas purification rate is predicted to deteriorate, the CPU determines the delay amount of the throttle drive operation and the torque correction amount based on the ignition timing, and matches the above estimated value Cgout ** with the target state. The air-fuel ratio is controlled (step 350), and the throttle is driven (step 360). The amount of delay in the throttle drive operation is, for example, the actual throttle valve opening that is obtained by performing a first-order delay process on the actual target throttle valve opening described above in order to suppress a rapid change in the throttle valve opening. This is the amount necessary to reset as. The ignition timing correction amount is a correction amount (advance amount) for compensating for a decrease in engine torque caused by resetting the actual throttle valve opening. The air-fuel ratio control for setting the estimated value Cgout ** to the target state is the same as in the other embodiments described above.
ここまでに説明した各実施形態は、いずれも触媒19全体を一つの領域として考えて各推定値を算出するものであつた。これに対し、以下の実施形態は、触媒19を排ガスの流れ方向に複数の領域(ブロック)に分割し、その複数の領域のうちの一つの領域(特定領域J)についての各推定値に対して目標状態を設定して、同各推定値が目標状態となるように空燃比制御を行う。
In each of the embodiments described so far, each estimated value is calculated considering the
まず、触媒19を複数領域に分割する実施形態のうち、第六実施形態について、その制御のフローチャートを示した図24を参照しながら説明する。この実施形態は、触媒19をn段に分割して取り扱う。ここで、触媒19の分割された複数の領域について、上流側からI=1〜nの番号を順に付するとともに、各領域I毎の推定値等にI=1〜nの添え字を付する。
First, among the embodiments in which the
また、本実施形態では、排ガス中の酸素量(触媒19の各領域から排出される酸素の濃度)と一酸化炭素量(触媒19の各領域から排出される一酸化炭素の濃度)とを上述した推定値として求める。更に、酸素量は空燃比制御において第一実施形態と同様に過不足量として扱う。なお、既に述べたように、上記(27)式〜(34)式、(36)式、及び(38)式等による各数値は、これらの式に基づいて触媒19の上流側から順次計算することによって、すべての領域について算出できる。
In the present embodiment, the amount of oxygen in exhaust gas (the concentration of oxygen discharged from each region of the catalyst 19) and the amount of carbon monoxide (the concentration of carbon monoxide discharged from each region of the catalyst 19) are described above. As an estimated value. Furthermore, the oxygen amount is handled as an excess / deficiency amount in the air-fuel ratio control as in the first embodiment. As already described, the numerical values according to the equations (27) to (34), (36), and (38) are calculated sequentially from the upstream side of the
以下、この実施形態の制御方法について図24に示したフローチャートを参照しながら具体的に説明すると、まず、CPUは上流側空燃比センサ25が検出する触媒19に流入する排ガスの空燃比を取り込み(ステップ370)、上述した特定領域を表す値Iを「0」にリセットする(ステップ380)。次いで、CPUは、I←I+1として、Iの値を「1」だけインクリメントする(ステップ390)。この段階では、Iの値は1となる。
Hereinafter, the control method of this embodiment will be described in detail with reference to the flowchart shown in FIG. 24. First, the CPU takes in the air-fuel ratio of the exhaust gas flowing into the
次いで、CPUは、領域Iに対するRoscI(即ちR*stor,O2,I、R*reudc,CO,I)及びRDI(RD,O2,I、RD,CO,I)を算出する(ステップ410)。 Then, CPU is, RoscI for the region I (i.e. R * stor, O2, I, R * reudc, CO, I) and R D I (R D, O2 , I, R D, CO, I) is calculated ( Step 410).
このとき、触媒19に流入する排ガスの空燃比が理論空燃比よりリーンの場合、RD,Iの算出には既に説明した(50)式を用い、Rosc,Iの算出には既に説明した(51)式を用いる。また、触媒19に流入する排ガスの空燃比が理論空燃比よりリッチの場合、RD,Iの算出には既に説明した(52)式を用い、Rosc,Iの算出には既に説明した(53)式を用いる。
At this time, when the air-fuel ratio of the exhaust gas flowing into the
その後、CPUは上述した第一実施形態の図18のステップ40と同様に、(30)式と(34)式とを用いてCgoutI(Cgout,O2,IとCgout,CO,I)を算出する(ステップ420)。このとき、CPUは、ステップ420の処理を行う際に必要とするCin,O2,I等をステップ40と同様に上記(48)式に基づいて算出する。なお、Cgout,CO,Iが求められたときには、Cgout,O2,Iに−2・Cgout,CO,Iを設定する。
After that, the CPU calculates CgoutI (Cgout, O2, I and Cgout, CO, I) using the equations (30) and (34) as in
次いで、CPUはコート層の化学種濃度CwI(Cw,O2,I、Cw,CO,I)を(28)式及び(32)式に基づいて求め(ステップ430)、(42)式〜(44)式によりOstI(今回値)を算出するとともに(ステップ435)、現在のIの値がn以上であるか否かを判定する(ステップ440)。このとき、Iの値がn未満であれば、CPUはステップ390に戻り、Iの値を「1」だけ大きくして次の下流側の領域について同様の算出を行う(ステップ410〜430,ステップ435)。一方、Iの値がn以上となると、すべての領域についての各種値の算出が終了したことになるので、この場合、CPUはステップ440からステップ450に進む。
Next, the CPU obtains the chemical species concentration CwI (Cw, O2, I, Cw, CO, I) of the coat layer based on the equations (28) and (32) (step 430), and the equations (42) to (44). ) Calculates OstI (current value) (step 435) and determines whether the current value of I is n or more (step 440). At this time, if the value of I is less than n, the CPU returns to step 390, increases the value of I by “1”, and performs the same calculation for the next downstream region (
そして、CPUはステップ450にて空燃比のフィードバック補正量(F/B補正量)を、特定領域J(Jはnより小さいか等しい)についての推定値CgoutO2,Jに基づいて、(CgoutO2,J−O2ref)×Gから求める。ここで、O2refは制御目標、Gは制御ゲインである。決定されたフィードバック補正量は、空燃比を決定する燃料噴射量の補正係数の一つとして空燃比制御に反映される。
In
このように、本実施形態においては、触媒19を複数の領域に分割し、特定領域Jに関する上述した推定値を求めているので、同特定領域Jの推定値に基づいて、触媒19の状況をより正確に把握できる。また、特定領域Jは、最も下流側の領域よりも上流側の領域に設定することができ、その結果、排ガスの浄化性能をより一層向上させることが可能となる。換言すると、特定領域J(J<n)の状態を理想の状態に近づけておけば、予期しない多量の未燃成分や窒素酸化物が触媒19に流入したとしても、特定領域J+1〜nの領域にて、これらを浄化する可能性が高まる。即ち、空燃比制御遅れによる弊害(エミッションの悪化)を小さくすることができる。
Thus, in the present embodiment, the
次に、第七実施形態について、その制御のフローチャートを示した図25を参照しながら説明する。本実施形態は、各領域毎の各種数値の算出に関する限り、第六実施形態と何ら変わるところがない。このため、同一の制御ステップには同一の符号を付して、その説明を省略する。 Next, a seventh embodiment will be described with reference to FIG. 25 showing a flowchart of the control. This embodiment is not different from the sixth embodiment as far as various numerical values are calculated for each region. For this reason, the same code | symbol is attached | subjected to the same control step, and the description is abbreviate | omitted.
本実施形態において、CPUがステップ440にて「Yes」と判定すると(即ち、すべての領域についての各数値の算出が終了すると)、同CPUは特定領域(制御するべき領域、以下、特定制御領域Jとも言う。)を決定するための運転条件を取得し(ステップ460)、取得した運転条件に基づいて特定制御領域Jを決定する(ステップ470)。上述した運転条件としては、様々な条件を設定することができる。また、一つの条件だけでなく、複数の条件に基づいて特定制御領域Jを決定することもできる。ここでは、一つの条件に基づいて特定制御領域Jを決定する具体例を四つ挙げる。図26〜図29は、各例において使用するマップを示している。 In the present embodiment, when the CPU determines “Yes” in step 440 (that is, when the calculation of each numerical value for all areas is completed), the CPU determines that the specific area (the area to be controlled, hereinafter, the specific control area). (Also referred to as J) is acquired (step 460), and the specific control region J is determined based on the acquired operating condition (step 470). Various conditions can be set as the operating conditions described above. In addition, the specific control region J can be determined based not only on one condition but also on a plurality of conditions. Here, four specific examples of determining the specific control region J based on one condition are given. 26 to 29 show maps used in each example.
図26は、上述した運転条件として吸入空気量を用いた場合に使用する「吸入空気量と特定領域(特定制御領域)Jとの関係」を示すマップである。吸入空気量は、エアフロメータ13によって検出される。このマップによれば、吸入空気量が多いほど特定領域(特定制御領域)Jが上流側(上流段)に設定される。
FIG. 26 is a map showing the “relationship between the intake air amount and the specific region (specific control region) J” used when the intake air amount is used as the operating condition described above. The intake air amount is detected by the
図27は、上述した運転条件としてアクセル開度の時間変化量(Δアクセル開度)を用いた場合に使用する「Δアクセル開度と特定領域(特定制御領域)との関係」を示すマップである。Δアクセル開度は、アクセルポジシヨンセンサ12の検出結果を加工することによって得られる。このマップによれば、Δアクセル開度が大きいほど特定領域(特定制御領域)Jが上流側に設定される。
FIG. 27 is a map showing the “relationship between the Δ accelerator opening and the specific area (specific control area)” used when the time change amount of the accelerator opening (Δ accelerator opening) is used as the operation condition described above. is there. The Δ accelerator opening is obtained by processing the detection result of the
吸入空気量やΔアクセル開度に関しては、その値が大きくなるほど特定領域(特定制御領域)Jの位置を上流側に設定し、その値が小さくなるほど特定領域(特定制御領域)Jの位置を下流側に設定する。これは、いわゆる「吹き抜け現象」が発生しやすい場合には空燃比制御に用いられる特定領域(特定制御領域)Jを上流側に設定し、「吹き抜け現象」が発生しにくい場合には特定領域(特定制御領域)Jを下流側に設定するためである。「吹き抜け現象」とは、触媒19自体は浄化する能力を有しているにもかかわらず、流速が速かったり、流量が多かったりという理由によって、排ガス中の浄化されるべき成分が十分に浄化される前に下流側に流出してしまう現象である。
As for the intake air amount and the Δ accelerator opening, the position of the specific area (specific control area) J is set upstream as the value increases, and the position of the specific area (specific control area) J decreases as the value decreases. Set to the side. This is because the specific region (specific control region) J used for air-fuel ratio control is set upstream when the so-called “blow-through phenomenon” is likely to occur, and the specific region ( This is because the specific control area J is set on the downstream side. “Blow-through phenomenon” means that although the
図28は、上述した運転条件として触媒19の活性状態を用いた場合に使用する「触媒の活性状態と特定領域(特定制御領域)Jとの関係」を示すマップである。触媒19の活性状態は、例えば、触媒温度センサ21によって検出した触媒19の温度から決定することができる。このマップによれば、触媒19の活性状態が低い(まだ十分に活性化していない)ほど特定領域(特定制御領域)Jが上流側に設定される。これにより、特定領域(特定制御領域)Jよりも下流側の領域をバッファ(排気浄化のための予備的な領域)として使用することができるので、排ガス中の浄化すべき成分が十分に浄化され得る。
FIG. 28 is a map showing the “relationship between the active state of the catalyst and the specific region (specific control region) J” used when the active state of the
図29は、上述した運転条件として空燃比(触媒19に流入するガスの空燃比、即ち、エンジン1に供給される混合気の空燃比)の時間的変化量(Δ空燃比)を用いた場合に使用する「Δ空燃比と特定領域(特定制御領域)Jとの関係」を示すマップである。Δ空燃比は、空燃比制御を司っているECU18によって算出することができる。このマップによれば、Δ空燃比が小さい(空燃比が余り変化していない場合)ほど特定領域(特定制御領域)Jが下流側に設定される。
FIG. 29 shows the case where the temporal change amount (Δair-fuel ratio) of the air-fuel ratio (the air-fuel ratio of the gas flowing into the
このように、CPUは図25に示したステップ470にて特定領域(特定制御領域)Jを決定すると、ステップ480に進んで同決定された特定領域(特定制御領域)Jの推定値Cgout,O2,Jに基づいて、空燃比フィードバック補正量(F/B補正量)を、(Cgout,O2,J−O2ref)×Gに基づいて求める。ここで、O2refは制御目標、Gは制御ゲインである。決定された補正量は、空燃比を決定する燃料噴射量の補正係数の一つとして空燃比制御に反映される。このように、触媒19を複数の領域に分割し、運転状態に応じて特定領域(特定制御領域)Jを変更することによって、排ガスの浄化性能をより一層向上させることが可能となる。
As described above, when the CPU determines the specific area (specific control area) J in
上述した第六実施形態および第七実施形態においては、特定領域(特定制御領域)は一つだけであった。しかし、特定領域(特定制御領域)は一つでなくても良く、複数設定されてもよい。このようにすることによって、排ガスの浄化をより一層効果的に行うことができるようになる場合がある。次に説明する第八実施形態は、特定領域(特定制御領域)を複数設定している。 In the sixth embodiment and the seventh embodiment described above, there is only one specific area (specific control area). However, the number of specific areas (specific control areas) is not limited to one, and a plurality of specific areas may be set. By doing so, the exhaust gas may be more effectively purified in some cases. In an eighth embodiment to be described next, a plurality of specific areas (specific control areas) are set.
第八実施形態における制御のフローチャートを図30に示す。本実施形態は、各領域毎の各種数値の算出に関して、上述した第六実施形態及び第七実施形態と何ら変わるところがない。このため、同一の制御ステップには同一の符号を付して、その説明を省略する。 A flowchart of control in the eighth embodiment is shown in FIG. The present embodiment is not different from the above-described sixth and seventh embodiments with respect to the calculation of various numerical values for each area. For this reason, the same code | symbol is attached | subjected to the same control step, and the description is abbreviate | omitted.
本実施形態においては、CPUがステップ440にて「Yes」と判定すると(即ち、すべての領域についての各数値の算出が終了すると)、CPUはステップ490に進み、あらかじめ決定されている二つの特定領域(特定制御領域)(領域I1及び領域I2)の推定値Cgout,O2,I1及びCgout,O2,I2に基づいて空燃比フィードバック補正量(F/B補正量)を[(Cgout,O2,I1−O2ref)×G1+(Cgout,O2,I2−O2ref)×G2]から求める。ここで、O2refは制御目標であり、本実施形態においては、特定領域(特定制御領域)I1及び特定領域(特定制御領域)I2の双方について同一の制御目標が用いられている。このステップ490において決定された補正量は、空燃比を決定する燃料噴射量の補正係数の一つとして空燃比制御に反映される。なお、特定領域(特定制御領域)I1と特定領域(特定制御領域)I2とで異なる制御目標を設定することもできる。
In the present embodiment, when the CPU determines “Yes” in step 440 (that is, when the calculation of each numerical value for all the regions is completed), the CPU proceeds to step 490 and performs two specific determinations that are determined in advance. region (specific control region) estimates CGout, O2, I 1 and CGout, O2, air-fuel ratio feedback correction amount based on I 2 of (region I 1 and region I 2) to (F / B correction amount) [(CGout , obtained from O2, I 1 -O2ref) × G 1 + (Cgout, O2, I 2 -O2ref) × G 2]. Here, O2ref is the control target, in the present embodiment, the same control target is used for both the specific region (specific control region) I 1 and the specific region (specific control region) I 2. The correction amount determined in
また、G1,G2はそれぞれ特定領域(特定制御領域)I1,特定領域(特定制御領域)I2に関しての制御ゲインである。本実施形態では、制御ゲインG1,G2を異ならせて、各特定領域(特定制御領域)の推定値の空燃比制御への影響度についての重み付けをしている。このようにすることによって、特定領域(特定制御領域)を複数設定する場合に、それらの特定領域(特定制御領域)についての状況が空燃比制御に与える影響度を個別に設定でき、排気浄化性能をより向上させることができる場合がある。このように、触媒19を複数の領域に分割し、複数の特定領域(特定制御領域)を設定することによって、排ガスの浄化性能をより一層向上させることが可能となる。
G 1 and G 2 are control gains relating to the specific area (specific control area) I 1 and the specific area (specific control area) I 2 , respectively. In this embodiment, the control gains G 1 and G 2 are made different to weight the degree of influence of the estimated value of each specific region (specific control region) on the air-fuel ratio control. In this way, when multiple specific areas (specific control areas) are set, it is possible to individually set the degree of influence of the conditions of these specific areas (specific control areas) on the air-fuel ratio control, and the exhaust purification performance. In some cases. Thus, the exhaust gas purification performance can be further improved by dividing the
上述した第八実施形態においては、制御ゲインG1,G2はあらかじめ決定された固定値とされていた。しかし、特定領域(特定制御領域)を複数設定した場合の各特定領域(特定制御領域)に関する制御ゲインはあらかじめ固定値として決定されていなくても良く、運転状況などに応じて変動する変数値として設定されてもよい。次に説明する第九実施形態は、各特定領域(特定制御領域)に対応する制御ゲインを可変としている。 In the above-described eighth embodiment, the control gains G 1 and G 2 are fixed values determined in advance. However, when a plurality of specific areas (specific control areas) are set, the control gain for each specific area (specific control area) does not have to be determined in advance as a fixed value. It may be set. In the ninth embodiment described below, the control gain corresponding to each specific area (specific control area) is variable.
第九実施形態における制御のフローチャートを図31に示す。本実施形態は、各領域毎の各種数値の算出に関して、上述した第六実施形態〜第八実施形態と何ら変わるところがない。このため、同一の制御ステップには同一の符号を付して、その説明を省略する。 FIG. 31 shows a flowchart of control in the ninth embodiment. The present embodiment is not different from the above-described sixth to eighth embodiments with respect to calculation of various numerical values for each region. For this reason, the same code | symbol is attached | subjected to the same control step, and the description is abbreviate | omitted.
本実施形態においては、CPUがステップ440にて「Yes」と判定すると(即ち、すべての領域についての各数値の算出が終了すると)、CPUはステップ500に進んで各特定領域に対応する制御ゲインを決定するための運転状態である吸入空気流量Gaをエアフロメータ13から取得する。そして、CPUは、ステップ510にて、前記取得した吸入空気流量Gaと図32に示したマップと基づいて、各特定領域毎(特定領域I1及び特定領域I2)に対する制御ゲイン(G1、G2)を決定する。
In the present embodiment, when the CPU determines “Yes” in step 440 (that is, when the calculation of each numerical value for all the regions is completed), the CPU proceeds to step 500 and performs control gain corresponding to each specific region. The intake air flow rate Ga which is the operation state for determining the air flow is acquired from the
図32に示したマップにおいては、横軸が特定領域の位置を示し、縦軸がG1,G2として設定される制御ゲインを示している。また、このマップでは、特定領域I1は特定領域I2よりも下流側に位置し、吸入空気量Gaの大小に応じた二本の曲線が設定されている。各曲線は、上流側の特定領域I2の制御ゲインG2よりも下流側の特定領域I1の制御ゲインG1の方が大きくなるように設定されている。すなわち、このマップにより制御ゲインを決定すると、下流側の特定領域に関する推定値の空燃比制御に対する影響度が大きくなる。なお、これらの曲線は、適宜定めることができ、例えば、これらの曲線を、上流側に向うに連れて制御ゲインを大きくするように設定すれば、上流側の特定領域の推定値が重視されることになる。 In the map shown in FIG. 32, the horizontal axis indicates the position of the specific region, and the vertical axis indicates the control gain set as G 1 and G 2 . Further, in this map, the specific area I 1 is located on the downstream side of the specific area I 2 , and two curves are set according to the magnitude of the intake air amount Ga. Each curve is set to than the control gain G 2 of the specific region I 2 of the upstream side toward the control gain G 1 of the specific region I 1 downstream increases. That is, when the control gain is determined using this map, the degree of influence on the air-fuel ratio control of the estimated value related to the downstream specific region increases. These curves can be determined as appropriate. For example, if these curves are set so that the control gain increases toward the upstream side, the estimated value of the specific area on the upstream side is emphasized. It will be.
さらに、図32に示すマップでは、吸入空気量Gaが大きい場合に適用する曲線が、吸入空気量Gaが小さい場合に適用する曲線よりも上に位置されている。このため、制御ゲインG1,G2は、吸入空気量Gaが多いほど大きくなるように設定される。 Further, in the map shown in FIG. 32, the curve applied when the intake air amount Ga is large is positioned above the curve applied when the intake air amount Ga is small. Therefore, the control gains G 1 and G 2 are set so as to increase as the intake air amount Ga increases.
CPUは、このようにして制御ゲインを決定すると、空燃比フィードバック補正量(F/B補正量)を[(Cgout,O2,I1−O2ref)×G1+(Cgout,O2,I2−O2ref)×G2]なる式に基づいて求める(ステップ520)。O2refは制御目標である。 When the CPU determines the control gain in this manner, the air-fuel ratio feedback correction amount (F / B correction amount) is set to [(Cgout, O2, I 1 −O2ref) × G 1 + (Cgout, O2, I 2 −O2ref ) × G 2 ] (step 520). O2ref is the control target.
このように、触媒を複数の領域に分割して複数の特定領域を設定し、各特定領域毎の空燃比制御への影響度をエンジン1の運転状態に応じて変更することによって、排ガスの浄化性能をより一層向上させることが可能となる。
In this way, the catalyst is divided into a plurality of regions, a plurality of specific regions are set, and the degree of influence on the air-fuel ratio control for each specific region is changed according to the operating state of the
なお、制御ゲインを決定するための上述した運転状態としては、様々な運転状態を採用することができる。また、一つの運転状態だけでなく、複数の運転状態に基づいて各特定領域に対する制御ゲインを決定することもできる。また、本実施形態においては、特定領域I1及び特定領域I2の双方について同一の制御目標が用いられていたが、特定領域I1と特定領域I2とで異なる制御目標を設定してもよい。 In addition, various driving | running states are employable as the driving | running state mentioned above for determining a control gain. In addition, the control gain for each specific region can be determined based on not only one operating state but also a plurality of operating states. In the present embodiment, although the same control target has been used for both the specific region I 1 and the specific region I 2, be set different control targets the specific region I 1 and the specific region I 2 Good.
次に、上述した推定値に基づく空燃比制御を実行するとともに、下流側空燃比センサ26の出力を用いて上述した推定モデルを補正する第十実施形態について説明する。本実施形態における制御のフローチャートを図33に示す。図33のフローチャートは、空燃比を制御するためのプログラムを示したものではなく、推定モデルを補正する制御を行うためのプログラムを示したものである。なお、本実施形態においては、第一実施形態と同様に、触媒19全体を一つの領域として考えて推定値を算出し、排ガス中の酸素量(触媒19から流出する酸素濃度)を同推定値の一つとして算出する。
Next, a description will be given of a tenth embodiment that executes air-fuel ratio control based on the estimated value described above and corrects the estimated model described above using the output of the downstream air-
具体的に説明すると、CPUは、ステップ530にて、上流側空燃比センサ25が検出する触媒19に流入する排ガスの空燃比(排気空燃比)を取得し、ステップ540に進んで同取得した排気空燃比に基づいてCgout,O2(=CgoutO2)等を算出する。これらのステップ530,540に基く処理は、実際には、図18に示したステップ10〜ステップ50、及びステップ55に基く処理と同様な処理を実行することで達成される。
Specifically, in
次いで、CPUはステップ550に進み、ステップ540で算出されたCgoutO2に基づいて下流側空燃比センサ26が出力するであろう出力を推定(予測)し、ステップ560にて(下流側空燃比センサ26の実際の出力)−(ステップ550において推定されたセンサ出力推定値)をErrorO2mdlとして算出する。
Next, the CPU proceeds to step 550 to estimate (predict) an output that the downstream air-
そして、CPUは、続くステップ570にて、このErrorO2mdlの絶対値が、所定のモデル許容誤差Emdlよも大きいか否かを判定する。このとき、下流側空燃比センサ26の実際の出力とステップ550において推定された出力推定値との差が許容範囲内であれば、CPUはステップ570に「No」と判定し、この場合、上述した推定モデルに対して補正を行う必要はないので、図33に示したフローチャートによる制御を一旦終了する。
Then, in
これに対し、下流側空燃比センサ26の実際の出力とステップ550において推定された出力推定値との差ErrorO2mdlが許容範囲(モデル許容誤差Emdl)を越えていると、CPUはステップ570にて「Yes」と判定し、ステップ580に進んで、推定モデルで使用する(上述した(36)式及び(38)式等で使用する)kstor(Kstor,i)及びkrel(krel,i)を修正する。
On the other hand, if the difference ErrorO2mdl between the actual output of the downstream side air-
このように、上述した推定値を用いた空燃比制御を実行するとともに、この推定値と下流側空燃比センサ26の出力とを用いて推定モデルを補正することによって、同推定モデルの推定精度をより一層向上させることができる。
As described above, by executing the air-fuel ratio control using the estimated value described above and correcting the estimated model using the estimated value and the output of the downstream air-
上述した第十実施形態では、下流側空燃比センサ26の出力に基づいて推定モデルを補正したが、次に説明する第十一実施形態においては、推定モデルによる推定結果に基づいて下流側空燃比センサ26の異常判定(診断)を行う。本実施形態における制御のフローチャートを図34に示す。図34のフローチャートは、空燃比を制御するためのプログラムを示したものではなく、センサの診断(異常判定)を行うためのプログラムを示したものである。なお、本制御においては、第一実施形態と同様に、触媒19全体を一つの領域として考えて推定値を算出し、排ガス中の酸素量(触媒19からの排出される酸素の濃度)を同推定値の一つとして算出する。
In the tenth embodiment described above, the estimation model is corrected based on the output of the downstream side air-
具体的に説明すると、CPUは、ステップ530にて、空燃比センサ25が検出する触媒19に流入する排ガスの空燃比を検出取得し、ステップ540に進んで同取得した排気空燃比に基づいてCgoutO2を算出するとともに、ステップ550にて同算出されたCgoutO2に基づいて下流側空燃比センサ26が出力するであろう出力を推定(予測)する。これらのステップ530〜550は、上述した第十実施形態におけるステップ530〜550と同一である。
More specifically, the CPU detects and acquires the air-fuel ratio of the exhaust gas flowing into the
次いで、CPUはステップ590に進み、(下流側空燃比センサ26の実際の出力)−(ステップ550において推定されたセンサ出力推定値)をErrorO2snsとして算出する。 Next, the CPU proceeds to step 590 and calculates (actual output of the downstream air-fuel ratio sensor 26) − (sensor output estimated value estimated in step 550) as ErrorO2sns.
そして、CPUは、続くステップ600にて、このErrorO2snsの絶対値が、所定のセンサ許容誤差Esnsよりも大きいか否かを判定する。このとき、下流側空燃比センサ26の実際の出力とステップ550において推定された出力推定値との差ErrorO2snsが許容範囲(センサ許容誤差Esns)内であれば、CPUはステップ600にて「No」と判定し、この場合、下流側空燃比センサ26は正常であると判断し、図34に示したフローチャートによる制御を一旦終了する。
Then, in
これに対し、下流側空燃比センサ26の実際の出力とステップ550において推定されたセンサ出力推定値との差ErrorO2snsが許容範囲(センサ許容誤差Esns)を越えていると、CPUはステップ600にて「Yes」と判定し、ステップ610に進んで、下流側空燃比センサ26が異常であると判定する。
On the other hand, if the difference ErrorO2sns between the actual output of the downstream side air-
この実施形態によれば、上述した推定値を用いて空燃比制御が実行されるとともに、この推定値と下流側空燃比センサ26の出力とに基づいて同空燃比センサ26の診断が行われ得る。
According to this embodiment, air-fuel ratio control is executed using the estimated value described above, and diagnosis of the air-
第十実施形態では、下流側空燃比センサ26の出力に基づいて推定値を求めるための推定モデルを補正した。これに対し、第十一実施形態では、推定モデルにより算出される推定値に基づいて下流側空燃比センサ26の診断を行った。これらは、互いに矛盾する考え方ではあるが、推定値と下流側空燃比センサ26の出力とのうち何れが信頼性が高いかに基づいて選択すればよい。また、エンジン1の運転状態に応じて、信頼できる方が変わるような場合は、エンジン1の運転状態に応じて推定モデル補正制御と下流側空燃比センサ26の診断制御とが排他的に行われるようにしてもよい。また、第十一実施形態のように、下流側空燃比センサ26の正常異常を判定するだけでなく、ErrorO2snsに応じて下流側空燃比センサ26の出力を補正するように構成してもよい。
In the tenth embodiment, the estimation model for obtaining the estimated value is corrected based on the output of the downstream side air-
これまで説明してきた実施形態は、図1に示されるような、排気通路7上に一つの触媒19のみを有するエンジン1に適用されていた。これに対し、以下に説明する第十二実施形態では、全体構成である図35に示したように、エンジン1がその排気通路7上に複数(以下の実施形態では二つ)の排気浄化触媒19a,19bを有している。なお、図35において、図1に示した構成と同一の構成には同一の符号を付し、その説明を省略するとともに、以下においては、図1に示した内燃機関と相違している部分についてのみ説明する。
The embodiment described so far has been applied to an
図35に示されるこの実施形態にかかる内燃機関に関しては、排気通路7上の二ヶ所に排気浄化触媒19a,19bが配設されている。上流側の触媒を上流側排気浄化触媒19a(以下、上流側触媒19aとも言う。)、下流側の触媒を下流側排気浄化触媒19b(以下、下流側触媒19bとも言う。)と称呼して区別する。
In the internal combustion engine according to this embodiment shown in FIG. 35,
上流側触媒19aは、例えば、四気筒のエンジンにおいて、ある二気筒の排気管が一つにまとめられた箇所に一つ、残りの二気筒の排気管が一つにまとめられた箇所にもう一つというように、並列に計二つ設けられるような場合もある。なお、この場合、下流側触媒19bは、すべての排気管が一つにまとめられた部分よりも下流の排気管に設けられることが多い。
For example, in a four-cylinder engine, the
上流側触媒19aは、始動時触媒とも呼ばれ、容量の小さなものがシリンダ3にできるだけ近い位置に配設され、排ガスの熱によって冷間始動後に早期に活性化温度にまで昇温され、早期に排気浄化性能を発揮することを目的として配設されることが多い。
The
これに対して、下流側触媒19bは、アンダーフロア触媒とも呼ばれ、十分な容量を有するものが車両の床下に配設され、排ガス中の浄化すべき成分を確実に浄化することを目的として配設される場合が多い。そして、上流側触媒19aの上流側には、上流側触媒19aに流入する排ガスの空燃比を検出する上流側空燃比センサ25が配設されている。また、下流側触媒19bの下流側には、下流側触媒19bから流出する排ガスの空燃比を検出する下流側空燃比センサ26が配設されている。さらに、上流側触媒19aと下流側触媒19bとの間には、上流側触媒19aから流出し、下流側触媒19bに流入する排ガスの空燃比を検出する中間部空燃比センサ27が配設されている。
On the other hand, the
これらの空燃比センサ25,26,27は、ECU18に接続されており、その検出結果をECU18に送出すると共に、ECU18から供給される電力によって早期に昇温されるべくヒータが内蔵されている。また、上流側触媒19a及び下流側触媒19bには、それらの温度を検出する温度センサ21a,21bがそれぞれ取り付けられている。
These air-
次に、本実施形態の作動について説明すると、この実施形態では、上流側触媒19aと下流側触媒19bとの双方で、上述した推定モデルによる推定値の算出が行われ、これらに基づく空燃比制御が行われる。これらの推定値と、中間部空燃比センサ27の検出結果とに基づいて、上流側触媒19aから流出して下流側触媒19bに流入する排ガスの空燃比も制御される。なお、上流側触媒19aに適用する推定モデルの基本構成と、下流側触媒19bに適用する推定モデルの基本構成は、それぞれのモデルで使用するパラメータの数値が異なっている点を除き、同一である。
Next, the operation of the present embodiment will be described. In this embodiment, the estimated value is calculated by the above estimation model in both the
本実施形態における制御のフローチャートを図36に示す。図36のフローチャートは、推定値を算出するためのプログラムのみを示している。また、本実施形態においては、上流側触媒19a全体を一つの領域として考えて推定値を算出し、下流側触媒19b全体を一つの領域として考えて別の推定値を算出する。この実施形態の推定値の一つは、排ガス中の酸素量(触媒19a,19bから流出する酸素の濃度)である。
FIG. 36 shows a flowchart of control in this embodiment. The flowchart in FIG. 36 shows only the program for calculating the estimated value. In the present embodiment, the estimated value is calculated considering the entire
また、上述したCgin,Cgoutに関しては、各触媒毎に以下のように表すこととする。即ち、上流側触媒19aに流入する化学種(ここでは酸素)濃度はCginO2SC、同触媒19aから流出する化学種濃度はCgoutO2SCと記述する。同様に、下流側触媒19bに流入する化学種(ここでは酸素)濃度はCginO2UF、同触媒19bから流出する化学種濃度はCgoutO2UFと記述する。上流側触媒19aから流出した化学種は下流側触媒19bに流入するので、CgoutO2SC=CginO2UFである。
In addition, Cgin and Cgout described above are expressed as follows for each catalyst. That is, the chemical species (in this case, oxygen) concentration flowing into the
具体的に、図36のフローチャートにより示されたプログラムに従って説明すると、CPUは、ステップ650にて、上流側空燃比センサ25が検出する上流側触媒19aに流入する排ガスの空燃比を取得し、ステップ652に進んで、同取得した排気空燃比に基づいて上流側触媒19aに関するCgoutO2SCを算出する。これらのステップ650,652に基く処理は、実際には、図18に示したステップ10〜ステップ50、及びステップ55に基く処理と同様な処理を実行することで達成される。次いで、CPUは、ステップ654にて、前記ステップ652で算出されたCgoutO2SCから中間部空燃比センサ27が出力するであろう出力を推定(予測)する。
Specifically, according to the program shown in the flowchart of FIG. 36, the CPU acquires the air-fuel ratio of the exhaust gas flowing into the
そして、CPUは、ステップ656に進み、(中間部空燃比センサ27の実際の出力)−(ステップ654において推定された出力推定値)をErrorO2mdlとして算出し、続くステップ658にて、このErrorO2mdlの絶対値が、所定のモデル許容誤差Emdlよりも大きいか否かを判定する。このとき、中間部空燃比センサ27の実際の出力とステップ654において推定された出力推定値との差が許容範囲を越えていると、CPUはステップ658にて「Yes」と判定してステップ660に進み、同ステップ660にてCgoutO2SCを補正し、ステップ662に進む。
Then, the CPU proceeds to step 656 to calculate (actual output of intermediate air / fuel ratio sensor 27) − (output estimated value estimated in step 654) as ErrorO2mdl, and in
このステップ660におけるCgoutO2SCの補正は、まず推定モデルを補正し、次いで補正した推定モデルに基づいて再度CgoutO2SCを算出することにより達成される。推定モデルの補正は、上述した図33のフローチャートにおけるステップ580での補正と同様に行われる。
The correction of CgoutO2SC in
一方、中間部空燃比センサ27の実際の出力とステップ654において推定された出力推定値との差が許容範囲内であると、CPUはステップ658にて「No」と判定し、上述した推定モデルに対して補正を行う必要がないので、ステップ662に直接進んでCgoutO2SCをCginO2UFに設定する。なお、ステップ660を経由してステップ662に到達した場合、ステップ660にて補正されたCgoutO2SCがCginO2UFに設定される。
On the other hand, if the difference between the actual output of the intermediate air-
次いで、CPUは、ステップ664にて、上記算出した下流側触媒に関するCginO2UFに基づいてCgoutO2UFを、ステップ652と同様に算出する。そして、CPUは図示しないステップにて、このように算出されたCgoutO2SC及びCgoutO2UFに基づいて空燃比フィードバック補正量を算出し、空燃比制御に反映する。このように、排気通路7上の複数の触媒19a,19bの双方に関して推定値をそれぞれ算出し、これらの推定値を用いた空燃比制御を行うことによって、排ガスの浄化性能をより一層向上させることが可能となる。
Next, in
次に、第十三実施形態について説明する。本実施形態は、上述した第十二実施形態と類似しているが、CgoutO2SCの補正は行わなわれない。なお、第十三実施形態の制御を第十二実施形態の制御と組み合わせて行うようにしても構わない。第十三実施形態は、下流側触媒19bの浄化率を高く維持するため、算出した推定値を用いて、下流側触媒19bに流入する排ガスの空燃比を制御する。具体的には、下流側触媒19bの浄化率を高く維持するために、同下流側触媒19bに流入する排ガス中の酸素過不足量の積算値がゼロとなるように、空燃比制御を行う。
Next, a thirteenth embodiment will be described. This embodiment is similar to the twelfth embodiment described above, but no correction of CgoutO2SC is performed. Note that the control of the thirteenth embodiment may be performed in combination with the control of the twelfth embodiment. In the thirteenth embodiment, in order to maintain the purification rate of the
上記下流側触媒19bに流入する酸素過不足量の積算値とは、下流側触媒19bの浄化率を高く維持させるために下流側触媒19bに流入する酸素量の収支(総和・平均値)をゼロにする際の余剰・不足状態を示すものである。本実施形態における制御のフローチャートを図37に示す。
The integrated value of the oxygen excess / deficiency flowing into the
このフローチャートに沿って、本装置の作動を説明すると、CPUはステップ650にて上流側空燃比センサ25が検出する上流側触媒19aに流入する排ガスの空燃比を取得し、ステップ652にて前記取得した空燃比に基づいて、上流側触媒19aに関するCgoutO2SCを算出する。これらのステップ650,652に基く処理は、実際には、図18に示したステップ10〜ステップ50、及びステップ55に基く処理と同様な処理を実行することで達成される。次いで、CPUはステップ662に進み、ステップ652で算出したCgoutO2SCをCginO2UFに設定する。
The operation of the present apparatus will be described along this flowchart. The CPU acquires the air-fuel ratio of the exhaust gas flowing into the
次いで、CPUはステップ670に進み、ステップ662において算出されたCginO2UFの積算値(ΣCginO2UFと表現する)を求める。即ち、CPUは、これまでの積算値ΣCginO2UF(k−1)をストックしており、これに対して新たにステップ662で算出されたCginO2UF(k)を加えることによって、新たな積算値ΣCginO2UF(k)を得る。
Next, the CPU proceeds to step 670 to obtain the integrated value of CginO2UF calculated in step 662 (expressed as ΣCginO2UF). That is, the CPU stocks the accumulated value ΣCginO2UF (k−1) so far, and newly adds CginO2UF (k) calculated in
そして、CPUは、ステップ672にて、前記算出された積算値ΣCginO2UF(k)がゼロより大きいか否かを判定し、同ステップ672にて「Yes」と判定するとき、ステップ674に進んで下流側触媒19bに流入する排ガスの空燃比がリッチとなるような制御(リッチ制御)を行う。一方、ステップ672の実行時に、前記算出された積算値ΣCginO2UF(k)がゼロ以下であれば、CPUはステップ672にて「No」と判定してステップ676に進み、下流側触媒19bに流入する排ガスの空燃比をリーンとするように制御(リーン制御)を行う。
Then, in
このように、上流側触媒19aと下流側触媒19bとの双方に関して推定値を算出し、これらの推定値に基づいて、下流側触媒19bを浄化率の高い状態に維持する空燃比制御を行うことによって、排ガスの浄化性能をより一層向上させることが可能となる。
In this way, the estimated values are calculated for both the
以上、説明したように、本発明の各実施形態によれば、排気浄化触媒に流入する排ガスの空燃比に基づいて、排気浄化触媒から排出される(同触媒の全部又は同触媒の一部である特定領域の一つから排出される)特定成分の排出量(或いは、排出される排ガスの状態を示す代表値)を推定し、推定した推定値に対する目標状態を設定し、この推定値が目標状態となるように空燃比制御を行う。これにより、排気浄化触媒或いは同触媒を分割した場合の特定領域(特定ブロック)から排出される特定成分の排出量(触媒又は特定領域から排出される排ガスの状態、触媒全体又は同触媒の最上流位置から特定領域までの酸素吸蔵量等で代表される触媒の状態を含む。)を所望の状態に制御でき、排気浄化性能を向上させることができる。 As described above, according to each embodiment of the present invention, the exhaust gas is exhausted from the exhaust gas purification catalyst based on the air-fuel ratio of the exhaust gas flowing into the exhaust gas purification catalyst (all of the catalyst or a part of the catalyst). Estimate the emission amount of a specific component (or a representative value indicating the state of exhaust gas emitted) from one specific area, and set a target state for the estimated value. Air-fuel ratio control is performed so that the state is reached. As a result, the exhaust amount of the specific component discharged from the exhaust purification catalyst or the specific region (specific block) when the catalyst is divided (the state of the exhaust gas discharged from the catalyst or the specific region, the entire catalyst or the uppermost stream of the catalyst) The state of the catalyst represented by the oxygen storage amount from the position to the specific region, etc.) can be controlled to a desired state, and the exhaust purification performance can be improved.
1…エンジン、7…排気通路、19…排気浄化触媒、25,26…空燃比センサ、18…電子制御ユニット。
DESCRIPTION OF
Claims (4)
前記排気浄化触媒の一部又は全部を通過した前記機関の排ガスに含まれる特定成分の排出量に応じた値を推定値として推定する推定手段と、
前記推定手段によって推定される推定値の少なくとも一つが所定の目標状態となるように、前記排気浄化触媒に流入する排ガスの空燃比を制御する空燃比制御手段と、
を備えた内燃機関の空燃比制御装置において、
前記特定成分は、前記排気浄化触媒に流入する前記機関の排ガスに含まれるとともに還元機能を有するリダクション成分及び同排ガスに含まれるとともに同リダクション成分に酸素を供給し得るストレージ成分のうちの少なくとも一つの成分であり、
前記推定手段は、前記特定成分の物質収支を考慮して形成された推定モデルに基づいて前記推定値を推定するように構成され、
前記推定手段の推定モデルは、前記排気浄化触媒を排ガスの流れ方向に複数の領域に分割したうちの特定領域に着目し、同特定領域の前記空間に流入する前記特定成分の量と同特定領域の空間から流出する同特定成分の量、及び同特定領域の空間から同特定領域の前記コート層に拡散する同特定成分の量に基づいて形成されたことを特徴とする空燃比制御装置。 An internal combustion engine in which an exhaust gas purification catalyst comprising a space through which an inflowing gas passes, and a coating layer that is exposed to the space and has a catalytic function and a substance that provides an oxygen storage / release function is disposed in the exhaust passage An air-fuel ratio control device for an engine,
Estimation means for estimating a value according to the emission amount of the specific component contained in the exhaust gas of the engine that has passed a part or all of the exhaust purification catalyst, as an estimated value;
Air-fuel ratio control means for controlling the air-fuel ratio of the exhaust gas flowing into the exhaust purification catalyst so that at least one of the estimated values estimated by the estimation means is in a predetermined target state;
In the air-fuel ratio control equipment of an internal combustion engine having a
The specific component is included in the exhaust gas of the engine flowing into the exhaust purification catalyst and has a reduction component having a reduction function and a storage component that is included in the exhaust gas and can supply oxygen to the reduction component. Ingredients,
The estimating means is configured to estimate the estimated value based on an estimation model formed in consideration of a mass balance of the specific component;
The estimation model of the estimation means pays attention to a specific region among the exhaust purification catalyst divided into a plurality of regions in the exhaust gas flow direction, and the same specific region as the amount of the specific component flowing into the space of the specific region. An air-fuel ratio control device formed on the basis of the amount of the specific component flowing out from the space and the amount of the specific component diffusing from the space in the specific region to the coat layer in the specific region.
前記推定手段のモデルは、前記特定領域の空間から前記特定領域のコート層に拡散する前記特定成分の量と同コート層で消費される同特定成分の量とに基づいて形成されたことを特徴とする空燃比制御装置。 The air-fuel ratio control apparatus for an internal combustion engine according to claim 1 ,
The model of the estimating means is formed based on the amount of the specific component diffusing from the space of the specific region to the coat layer of the specific region and the amount of the specific component consumed in the coat layer. An air-fuel ratio control device.
前記特定成分は酸素である空燃比制御装置。 The air-fuel ratio control apparatus for an internal combustion engine according to claim 1 or 2 ,
An air-fuel ratio control apparatus in which the specific component is oxygen.
前記特定成分は酸素と一酸化炭素である空燃比制御装置。 The air-fuel ratio control apparatus for an internal combustion engine according to claim 1 or 2 ,
The air-fuel ratio control apparatus in which the specific components are oxygen and carbon monoxide.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2006346578A JP4422716B2 (en) | 2001-06-18 | 2006-12-22 | Air-fuel ratio control device for internal combustion engine |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2001183493 | 2001-06-18 | ||
JP2006346578A JP4422716B2 (en) | 2001-06-18 | 2006-12-22 | Air-fuel ratio control device for internal combustion engine |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2003505470A Division JP3968075B2 (en) | 2001-06-18 | 2002-06-18 | Air-fuel ratio control device for internal combustion engine |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2009250530A Division JP4821905B2 (en) | 2001-06-18 | 2009-10-30 | Air-fuel ratio control device for internal combustion engine |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2007077997A JP2007077997A (en) | 2007-03-29 |
JP4422716B2 true JP4422716B2 (en) | 2010-02-24 |
Family
ID=37938543
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2006346578A Expired - Fee Related JP4422716B2 (en) | 2001-06-18 | 2006-12-22 | Air-fuel ratio control device for internal combustion engine |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP4422716B2 (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008240528A (en) | 2007-03-24 | 2008-10-09 | Hitachi Ltd | Variable displacement vane pump |
JP5124541B2 (en) * | 2009-08-05 | 2013-01-23 | 本田技研工業株式会社 | Engine oil replacement judgment device |
US9297288B2 (en) | 2011-07-15 | 2016-03-29 | Toyota Jidosha Kabushiki Kaisha | Exhaust purification system for internal combustion engine |
-
2006
- 2006-12-22 JP JP2006346578A patent/JP4422716B2/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JP2007077997A (en) | 2007-03-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4821905B2 (en) | Air-fuel ratio control device for internal combustion engine | |
JP4088412B2 (en) | Air-fuel ratio control device for internal combustion engine | |
EP2061958B1 (en) | Catalyst deterioration monitoring system and catalyst deterioration monitoring method | |
KR20070089240A (en) | Method and device for determining an oxygen storage capacity of the exhaust catalyst of an internal combustion engine and method and device for determining a duration for exhaust probes of an internal combustion engine | |
JP3239699B2 (en) | Catalyst deterioration determination device for internal combustion engine | |
JP3912054B2 (en) | Exhaust gas purification device for internal combustion engine | |
JP4422716B2 (en) | Air-fuel ratio control device for internal combustion engine | |
JP4973370B2 (en) | Exhaust gas purification device for internal combustion engine | |
JP5464391B2 (en) | Air-fuel ratio control device for internal combustion engine | |
JP2004169607A (en) | Control device for internal combustion engine | |
JP5545631B2 (en) | Air-fuel ratio control device | |
JP3787913B2 (en) | Exhaust gas purification device for internal combustion engine | |
JP2009091921A (en) | Catalyst deterioration diagnosis device for internal combustion engine | |
JP2000337130A (en) | Exhaust emission control system for internal combustion engine | |
JP4069367B2 (en) | Exhaust gas purification device for internal combustion engine | |
JP5679839B2 (en) | Air-fuel ratio control device | |
JPH05272329A (en) | Degradation detecting method and device for engine exhaust gas purifying catalyst | |
WO2023238361A1 (en) | Control method and control device for internal combustion engine | |
JP2004232576A (en) | Exhaust emission control device for internal combustion engine | |
JP4743637B2 (en) | Catalyst deterioration detection device for internal combustion engine | |
JP5088632B2 (en) | Catalyst deterioration diagnosis device | |
JP2007327475A (en) | Catalyst typical temperature acquiring device of internal combustion engine | |
JP2004116295A (en) | Exhaust emission control device of internal combustion engine | |
JP2009041515A (en) | Catalyst deterioration diagnosis device for internal combustion engine | |
JP2004270581A (en) | Emission control device for internal combustion engine |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20061225 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20070117 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20090901 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20091030 |
|
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: 20091125 |
|
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: 20091204 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121211 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121211 Year of fee payment: 3 |
|
S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313532 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121211 Year of fee payment: 3 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131211 Year of fee payment: 4 |
|
LAPS | Cancellation because of no payment of annual fees |