JP3772584B2 - 燃料性状検出装置 - Google Patents
燃料性状検出装置 Download PDFInfo
- Publication number
- JP3772584B2 JP3772584B2 JP09870999A JP9870999A JP3772584B2 JP 3772584 B2 JP3772584 B2 JP 3772584B2 JP 09870999 A JP09870999 A JP 09870999A JP 9870999 A JP9870999 A JP 9870999A JP 3772584 B2 JP3772584 B2 JP 3772584B2
- Authority
- JP
- Japan
- Prior art keywords
- fuel
- model
- fuel property
- frequency
- supply amount
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Landscapes
- Combined Controls Of Internal Combustion Engines (AREA)
- Testing Of Engines (AREA)
Description
【発明の属する技術分野】
この発明は燃料性状を検出する装置、特にガソリンエンジン用燃料の燃料性状を検出するものに関する。
【0002】
【従来の技術】
エンジンが不安定となる冷間始動時に、始動後増量補正係数により燃料増量を行い、空燃比を理論空燃比よりもリッチ側の値とすることによってエンジンを安定させるようにしたものがある(特開平6−105129号公報参照。)
【0003】
【発明が解決しようとする課題】
ところで、上記の噴射量補正が必要となる理由は次の通りである。吸気ポートに臨んで設けられている燃料噴射弁で具体的に述べると、噴射弁からの噴射燃料は、その全てが噴霧のまま気流に乗ってシリンダに導入されるのではなく、吸気弁傘裏部や吸気ポート壁に付着して液状となる。この吸気ポート壁に付着した燃料は、液状のままポート壁をゆっくりと伝いつつシリンダに到達する。つまり、この分の燃料(壁流燃料)だけは燃料供給遅れが生じてしまうのである。
【0004】
この場合、冷却水温が低くなるほど大きくなるこの壁流燃料量は、さらに燃料の性状(特に揮発性)にも依存し、揮発性が悪い燃料ほど壁流燃料量が多くなる。こうした燃料の揮発性の違いによる燃料壁流量差を考慮し、従来の各種噴射量補正では、市販されている燃料のうち最も揮発性の悪い燃料(最重質ガソリン)が使用される場合でも、冷間時のエンジン回転が不安定とならないようにマッチングされている。上記の始動後増量補正係数であれば、始動後増量補正係数を演算する際に用いるテーブル値を最重質ガソリンに対してマッチングしておくわけである。
【0005】
しかしながら、最重質ガソリンよりも揮発性の良い燃料が使用されるときは、各種噴射量補正量が大きくなり過ぎ、これによって最重質ガソリン使用時よりも空燃比がリッチ側に傾くため、排気エミッション(特にCO、HC)が悪くなる。
【0006】
そこで本発明は、過渡時に燃料噴射量に対する排気空燃比の応答波形をサンプリングし、これら過渡時データに基づいて、予めECM(エレクトロニックコントロールモジュール)上に構築したプラントモデルのパラメータを、基準燃料に対するプラントモデルである規範モデルとの予測誤差が最小となるように調整することにより、使用燃料に対するプラントモデルを同定し、この同定したプラントモデルのカットオフ周波数と規範モデルに対するカットオフ周波数とを比較することにより、使用燃料の燃料性状の推定を可能とすることを目的とする。
【0007】
【課題を解決するための手段】
第1の発明は、図28に示すように、エンジンの運転条件に応じた燃料供給量を演算する手段51と、この燃料供給量をエンジンに供給する手段52と、エンジンの排気空燃比を検出する手段53と、過渡時に前記燃料供給量を入力、前記排気空燃比を出力として現在の使用燃料での前記燃料供給量に対する排気空燃比の応答波形のデータをサンプリングする手段54と、これら入出力データに基づいて、予め構築したプラントモデルのパラメータを規範モデルとの予測誤差が小さくなるように調整することにより、前記使用燃料に対するプラントモデルを同定する手段55と、この同定されたプラントモデルのカットオフ周波数fcRealを演算する手段56と、このプラントモデルのカットオフ周波数fcRealと規範モデルのカットオフ周波数fcRefとを比較して前記使用燃料の燃料性状を推定する手段57とを設けた。
【0008】
第2の発明は、図29に示すように、エンジンの運転条件に応じた燃料供給量を演算する手段51と、この燃料供給量をエンジンに供給する手段52と、エンジンの排気空燃比を検出する手段53と、過渡時に前記燃料供給量を入力、前記排気空燃比を出力として現在の使用燃料での前記燃料供給量に対する排気空燃比の応答波形のデータをサンプリングする手段54と、これら入出力データに基づいて、予め構築したプラントモデルのパラメータを規範モデルとの予測誤差が小さくなるように調整することにより、前記使用燃料に対するプラントモデルを同定する手段55と、この同定されたプラントモデルのカットオフ周波数fcRealを演算する手段56と、このプラントモデルのカットオフ周波数fcRealと規範モデルのカットオフ周波数fcRefの差と許容範囲とを比較して前記使用燃料の燃料性状を推定する手段61とを設けた。
【0009】
第3の発明は、図30に示すように、エンジンの運転条件に応じた燃料供給量を演算する手段51と、この燃料供給量をエンジンに供給する手段52と、エンジンの排気空燃比を検出する手段53と、過渡時に前記燃料供給量を入力、前記排気空燃比を出力として現在の使用燃料での前記燃料供給量に対する排気空燃比の応答波形のデータをサンプリングする手段54と、これら入出力データに基づいて、予め構築したプラントモデルのパラメータを規範モデルとの予測誤差が小さくなるように調整することにより、前記使用燃料に対するプラントモデルを同定する手段55と、この同定されたプラントモデルのカットオフ周波数fcRealを演算する手段56と、カットオフ周波数に対する燃料性状推定値の特性を予め設定する手段71と、前記演算されたカットオフ周波数fcRealからこの特性を検索することにより燃料性状推定値を演算する手段72とを設けた。
【0010】
第4の発明では、第1または第2の発明において前記燃料性状の推定結果を不揮発性メモリ(たとえばEEPROM)に記憶させておく。
【0011】
第5の発明では、第3の発明において前記燃料性状推定値を不揮発性メモリ(たとえばEEPROM)に記憶させておく。
【0012】
第6の発明では、第1の発明において前記燃料性状の推定が、前記規範モデルを基準燃料に対してマッチングした場合に、前記同定したプラントモデルのカットオフ周波数が規範モデルのカットオフ周波数より高いとき、前記基準燃料よりも軽質であると推定することである。
【0013】
第7の発明では、第2の発明において前記燃料性状の推定が、前記規範モデルを基準燃料に対してマッチングした場合に、前記同定したプラントモデルのカットオフ周波数と規範モデルのカットオフ周波数の差が許容範囲外であり、かつプラントモデルのカットオフ周波数が規範モデルのカットオフ周波数よりも大きいとき、前記基準燃料よりも軽質であると推定することである。
【0014】
第8の発明では、第1から第7までのいずれか一つの発明において前記予測誤差が小さくなるように調整することが、予測誤差が最小となるように調整することである。
【0015】
第9の発明では、第1から第8までのいずれか一つの発明において前記プラントモデルを、燃料挙動モデル、排気モデルおよび無駄時間の3つより構成する。
【0016】
第10発明では、第9の発明において前記燃料挙動モデルを分母2次、分子2次の2次遅れ系モデルとする。
【0017】
第11の発明では、第9の発明において前記排気モデルを分母1次の1次遅れ系モデルとする。
【0018】
第12の発明では、第9の発明において前記無駄時間を適合項と演算項と判定項の和とする。
【0019】
第13の発明では、第1から第12までのいずれか一つの発明において前記同定の手法がARXモデルを用いた一括処理最小2乗法である。
【0020】
第14の発明では、第1から第13までのいずれか一つの発明において前記燃料性状の推定が運転の開始から停止までの間に1回だけである。
【0021】
【発明の効果】
プラントモデルのパラメータを規範モデルとの予測誤差が小さくなるように入出力データに基づいて調整することで、プラントモデルが同定される。このときのパラメータよりプラントモデルの伝達関数がわかり、これよりプラントモデルのカットオフ周波数が定まる。
【0022】
ここで、規範モデルを基準燃料に対してマッチングしてある場合に、使用燃料が基準燃料よりも軽質であるときは、プラントモデルのカットオフ周波数のほうが基準燃料に対するよりも高くなる(この逆に、使用燃料が基準燃料よりも重質であるときは、プラントモデルのカットオフ周波数のほうが基準燃料に対するよりも低くなる)。これは、基準燃料よりも軽質の燃料のほうが基準燃料よりも燃料輸送遅れが小さくなるため、その燃料応答性が高くなり、基準燃料と比べ、高周波域まで応答ゲインを維持できるためである。
【0023】
したがって、第1、第8、第9、第12、第13の発明によれば、プラントモデルのカットオフ周波数と規範モデルのカットオフ周波数とを比較することで、使用燃料の燃料性状を推定することが可能となる。たとえば、第6の発明のように規範モデルを基準燃料に対してマッチングしてある場合に、同定したプラントモデルのカットオフ周波数が規範モデルのカットオフ周波数より高いとき、また第7の発明のように規範モデルを基準燃料に対してマッチングしてある場合に、同定したプラントモデルのカットオフ周波数と規範モデルのカットオフ周波数の差が許容範囲外であり、かつプラントモデルのカットオフ周波数が規範モデルのカットオフ周波数よりも大きいとき、使用燃料は基準燃料よりも軽質であると推定すればよい。このようにして、燃料性状を推定することが可能となると、これを燃料供給量の補正に活かすことができる。
【0024】
第2の発明によれば、基準燃料に対するカットオフ周波数がバラツクことがあっても、燃料性状の推定を安定して行うことができる。
【0025】
第3の発明によれば、燃料供給量の演算精度を高めることができる。
【0026】
第4、第5の発明によれば、次回の運転時に始動当初より燃料性状の推定結果や燃料性状推定値を利用できる。
【0027】
燃料挙動に伴う無駄時間を燃料挙動の数学モデルに取り込んだのでは、モデル次数が上昇してしまうが、第10の発明によれば、無駄時間を燃料挙動の数学モデルに取り込んでないので、モデル次数の上昇を抑えることができる。
【0028】
排気系は遅れ系として、排気ガス輸送遅れ+ガス混合遅れ+センサ応答遅れが考えられ、それぞれが1次遅れ以上のモデルであるので、これを忠実に表現するとすれば、モデル次数が高次となってしまうのであるが、第11の発明によれば、排気モデルを分母1次の1次遅れ系モデルとしたので、モデル次数の上昇を抑えることができる。また、排気系の無駄時間を排気モデルに取り込んだのでは、モデル次数が上昇してしまうが、第11の発明によれば、無駄時間を排気モデルに取り込んでないので、モデル次数の上昇を抑えることができる。
【0029】
運転中に何度も燃料性状の推定を行わせるのでは、演算負荷が大きくなるが、第14の発明によれば1回だけの推定であるので、演算負荷を小さくすることができる。
【0030】
【発明の実施の形態】
図1において、1はエンジン本体で、吸入空気はエアクリーナから吸気管8を通ってシリンダに供給される。燃料は、運転条件に応じて所定の空燃比となるようにECM2よりの噴射信号に基づき燃料噴射弁7からエンジン1の吸気ポートに向けて噴射される。
【0031】
ECM2にはクランク角センサ4からのREF信号(4気筒では180°ごと、6気筒では120°ごとに発生する気筒を識別するための信号)と1°信号、エアフローメータ6からの吸入空気量信号、三元触媒10の上流側(排気マニフォールド集合部)に設置した広域空燃比センサ(以下単に「A/Fセンサ」という)3からの空燃比信号、水温センサ11からの冷却水温信号、スロットルセンサ12からの絞り弁5開度信号等が入力され、これらに基づいてECM2では、吸入空気量Qaとエンジン回転数Neとから基本噴射パルス幅Tpを演算するとともに、加減速時にはこのTpに過渡補正量Kathosを加算することによって壁流燃料に関する補正を行う。過渡補正量Kathosは、加減速時に限らず、壁流燃料が大きく変化する始動時や燃料リカバー時、さらには後述する目標当量比Tfbyaの切換時にも働く。
【0032】
ECM2ではまた、冷間始動時のエンジン安定性をよくしたり高負荷時の要求出力に応えるため目標当量比Tfbyaを用いて燃料補正を行うほか、トランスミッションのギヤ位置センサ13からのギヤ位置信号、車速センサ(図示しない)からの車速信号等に基づいて運転状態を判断しながら条件に応じてリーン空燃比と理論空燃比との制御を行う。排気管9には三元触媒10が設置され、理論空燃比の運転時に最大の転換効率をもって、排気中のNOxの還元とHC、COの酸化を行う。この三元触媒10はリーン空燃比のときはHC、COは酸化するが、NOxの還元効率は低い。しかしながら、空燃比がリーン側に移行すればするほどNOxの発生量は少なくなり、所定の空燃比以上では三元触媒10で浄化するのと同じ程度にまで下げることができ、同時に、リーン空燃比になるほど燃費が改善される。したがって、負荷のそれほど大きくない所定の運転領域においては目標当量比Tfbyaを1.0より小さな値とすることによってリーン空燃比による運転を行い、それ以外の運転領域ではTfbyaを1.0とすることにより空燃比を理論空燃比に制御するのである。
【0033】
このように目標当量比Tfbyaは運転条件の変化に応じて切換わるのであるが、上記の過渡補正量KathosをTfbya=1.0(つまり理論空燃比)に対する値として計算しているのでは、出力空燃比域(このときTfbyaは1.0より大きい)からの減速時などTfbyaの切換時に過渡補正量Kathosに不足を生じて空燃比が一時的にオーバーリッチやオーバーリーンになり、制御空燃比の追従性が悪くなるので、これに対処するためECM2では、平衡付着量Mfhを目標当量比Tfbyaをもパラメータとして演算している。これを演算式で表すと、
Mfh=Avtp×Mfhtvo×Tfbya×CYLDRN#
ただし、Mfh :全気筒トータルでの平衡付着量、
Avtp :噴射弁部空気量相当パルス幅、
Mfhtvo :付着倍率、
CYLDRN#:シリンダ数、
である(特開平10−18882号公報参照)。そして、この平衡付着量(壁流燃料の平衡値のこと)Mfhと現時点での付着量Mfとの差に分量割合Kmfを乗じる演算により、つまり
Vmf=(Mfh−Mf)×Kmf
の式により付着速度(単位周期当たりの付着量のこと)Vmf(1サイクル毎の値)を計算し、平衡付着量Mfhが増加する場合(たとえば加速時)にはこれを過渡補正量Kathosとし、
ただし、Kathos:過渡補正量(1サイクル毎の値)、
α :空燃比フィードバック補正係数、
Ts :無効噴射パルス幅、
Chosn1 :気筒別壁流補正量(各気筒1サイクル毎の値)、の式により、シーケンシャル噴射(4気筒ではエンジン2回転毎に1回、各気筒の点火順序に合わせて噴射)時に燃料噴射弁に与える実噴射パルス幅CTInを気筒別に計算する。なお、「1サイクル毎の値」とは、1REF信号の入力毎の値、「各気筒1サイクル毎の値」とは4REF信号(4気筒の場合)の入力毎の値である。CTInとChosn1の「n」は気筒番号を表す。
【0034】
ここで、上記の気筒別壁流補正量Chosnについて説明すると、壁流燃料には直接にシリンダに流入される分が少なく比較的応答の遅いもの(低周波成分という)と、直接にシリンダに流入される分が主で応答の速いもの(高周波成分という)とがあり、上記のVmfが低周波成分を対象とする壁流補正量であるのに対して、Chosnは高周波成分を対象とする補正量である。つまり、Vmfだけでは高周波成分に対して対処不可能なため、高周波成分に対する補正量であるChosnを導入する必要があるわけである。具体的には、前回噴射からの噴射弁部空気量相当パルス幅Avtpの変化量であるΔAvtpnを用いて、Avtpが増えているとき(加速時)であれば、
Chosn=ΔAvtpn×Gztwp
ただし、Gztwp:増量ゲイン、
の式により、またAvtpが減少しているとき(減速時)は
Chosn=ΔAvtpn×Gztwm
ただし、Gztwm:減量ゲイン、
の式により計算し、これを気筒別に同期噴射の燃料噴射パルス幅に加算することによって、高周波成分に対する壁流補正を行っている。なお、上記の増量ゲインGztwp、減量ゲインGztwmは水温補正を行うためのものである。また、ΔAvtpnの最後に添付されているnは、CTInの場合と同じに気筒番号を表す。
【0035】
このように低周波成分に対する壁流補正量に加えて高周波成分に対する壁流補正量を導入しているものにおいても、Chosnの演算にTfbyaが考慮されていないのでは、特に出力空燃比域からの減速時などTfbyaの切換時にChosnに不足を生じて一時的にオーバーリッチやオーバーリーンが生じるので、これに対処するためECM2では、高周波成分に対する壁流補正量であるChosnについてもTfbyaに応じた値としている(特開平10−18882号公報参照)。これを演算式で表すと、
ただし、Chosn1 :1サイクル目のChosn、
Kathos-4Ref:各気筒1サイクル前(4REF信号前)のKathos、
Gztwc :増量ゲインGztwpまたは減量ゲインGztwm、
A :低周波成分の1サイクル目の応答ゲイン、
である。
【0036】
なお、燃料カット(燃料カットを気筒別に行う場合と全気筒同時に行う場合とがある)を考慮して燃料リカバー時の壁流補正量(Chosn1 、Vmf)を計算することにより、目標当量比Tfbyaの切換を伴う燃料リカバー時にも、最適な壁流補正量を与えるようにしてもかまわない(特開平10−18882号公報参照)。
【0037】
一方、エンジンが不安定となる冷間始動時には、始動後増量補正係数KASにより燃料増量を行い、空燃比を理論空燃比よりもリッチ側の値とすることによってエンジンを安定させている(特開平6−101529号参照)。たとえば、後述する数36式によりKASを計算しており、このときのKASは、概略、スタータスイッチのON時の値を初期値としてスタータスイッチのOFFタイミングより急な勾配で直線的に減少し、途中からはより緩やかな勾配となって減少し、最後に0になる値である。
【0038】
そして、この始動後増量補正係数KASは、上記の目標当量比Tfbyaの一部であり、たとえば、
Tfbya=Kml+KAS
ただし、Dml:燃空比補正係数、
の式により目標当量比Tfbyaを算出している。
【0039】
ここで、Kmlは、運転条件に応じた目標空燃比を定めるもので、エンジンの回転数と負荷をパラメータとするマップを検索することにより求められる。なお、目標空燃比の切換時にKmlのマップ値を検索した上で所定のダンパ操作を行わせるものもある。
【0040】
また、始動時には特別な燃料噴射パルス幅TISTを設定している(特開平7−63082号公報参照)。急加速時のように特に大きなトルクが要求されるときは、加速時割り込み噴射パルス幅IJSETnが演算され、同期噴射の途中でも割り込み噴射が行われることがある。
【0041】
さて、上記のKASなど各種噴射量補正が必要となる理由は、壁流燃料の分だけの燃料供給遅れが生じてしまうからであるが、この壁流燃料量は、さらに燃料性状(特に揮発性)にも依存し、揮発性が悪い燃料ほど壁流燃料量が多くなる。こうした燃料の揮発性の違いによる壁流燃料量差を考慮し、従来の各種噴射量補正(たとえば従来の始動後増量補正)では、最重質ガソリンが使用される場合でも、冷間時のエンジン回転が不安定とならないようにKASの演算に用いるテーブル値がマッチングされている。
【0042】
しかしながら、最重質ガソリンよりも揮発性の良い燃料が使用されるときは、各種噴射量補正量が大きくなり過ぎ、これによって最重質ガソリンの使用時に期待される空燃比目標値よりも実際の空燃比がリッチ側に傾くため、排気エミッション(特にCO、HC)が悪くなる。
【0043】
そこで、壁流燃料量が多い領域での過渡時に、燃料噴射量に対する排気空燃比の応答波形をサンプリングし、これら過渡時データに基づいて予めECM2上に構築したプラントモデルのパラメータを、規範モデルとの推定誤差が最小となるように調整することによりプラントモデルを同定し、前記調整されたパラメータに基づいてプラントモデルのカットオフ周波数を求め、これと基準燃料(Ref燃料)に対するプラントモデルである規範モデルのカットオフ周波数とを比較することにより、燃料性状を推定し、この推定した燃料性状に応じて各種噴射量補正量を演算する。
【0044】
ECM2で実行されるこの最適化制御を次に説明する。
【0045】
図2は最適化制御の制御システムのブロック図である。
【0046】
本制御システムは大きく分けて、プラント同定部21、燃料性状推定部22、トリガリング機能23、コントローラ24から構成されている。
【0047】
図2を用いて制御の概要を説明し、その後にプラント同定について詳述する。
【0048】
まず、プラント同定部21の主な構成要素は、プラントモデル31、誤差検出手段32、最適化計算手段33、入力(実噴射パルス幅)のバッファリング手段34、出力(排気空燃比)のバッファリング手段35からなっており、エンジンパラメータより判定されたトリガより、実噴射パルス幅CTInおよびA/Fセンサ3出力電圧をサンプリングし、それぞれをプラントモデル31の入出力信号としてその領域においてのプラントモデル31の同定を行う。モデルの形式(次数)は予め物理モデルから設定してあり、実際の入出力信号に対して最適となるようにモデルパラメータを調整している。
【0049】
ここで、プラントモデル31は、燃料挙動特性を分母2次、分子2次の2次遅れ系モデル、排気動特性を分母1次の1次遅れ系モデルとした離散系カスケード結合のモデル(分母3次、分子3次の物理モデル)である。同定手法は最も一般的であるARXモデルを用いた一括処理最小2乗法としている。
【0050】
燃料性状推定部22は、基準燃料(Ref燃料)に対するプラントモデルである規範モデル37、カットオフ周波数比較手段38からなり、プラント同定部21により同定されたプラントモデル31のカットオフ周波数fcRealと規範モデル37のカットオフ周波数fcRefとを比較手段38において比較することにより燃料性状を判定する。
【0051】
ここでは、説明の簡単化のため、使用燃料が2種類しかない場合(揮発性の悪い方の燃料を重質ガソリン、揮発性の良いほうの燃料を軽質ガソリンとする)で説明すると、重質ガソリンに対してマッチングした規範モデルを用いた場合、同定したプラントモデル31のカットオフ周波数が規範モデルのカットオフ周波数より高いとき、軽質ガソリンが使用されていると判定できる。これは、軽質ガソリンのほうが重質ガソリンより燃料輸送遅れが小さくなるため、軽質ガソリンのほうが燃料応答性が高くなり、重質ガソリンと比べ、高周波域まで応答ゲインを維持できるためである。なお、燃料性状の推定前には、燃料性状の判定値を重質ガソリンに初期設定しておく。
【0052】
本制御では1トリップ(エンジンの運転の開始から停止まで)の間にたとえば1回、燃料性状を推定し、その推定結果を不揮発性メモリ(本実施形態においてはEEPROM14、図1参照)に記憶させておき、コントローラ24により次回始動時の燃料制御に反映させる。
【0053】
トリガリング機能23は、プラント同定に必要な入出力信号をサンプリングするとともに、プラント同定を開始するトリガを発生させる条件を判定するものである。
【0054】
ここで、一般的にシステムを同定するためには、広い帯域の周波数を含んだ入力が必要である。エンジンにおいては、M系列等の入力を生成するのは現実的でないので、入力がステップ的に変化する点をサンプリングトリガとする。また、EGR(排気還流)やスワールコントロール等、燃料挙動に大きく影響を及ぼす条件も排除する必要がある。上記のような条件が同定するときのモード中に存在するかまたは同定用入力信号を生成する必要がある。
【0055】
コントローラ24は、具体的には始動後増量補正係数調整手段で、1トリップ中に推定された燃料性状に応じて、次回の始動直後の始動後増量補正量KASを適切な値とする。調整する項目はKASの初期値と減衰割合である。
【0056】
ここでは、前述のように制御の簡略化のため、重質ガソリン用か軽質ガソリン用かの2値切換で説明するが(したがって、始動後増量補正係数調整手段は始動後増量補正係数切換手段40となる)、燃料性状の分離性能とエンジンの要求から燃料性状の切換の段数を決定すればよい。
【0057】
燃料噴射量演算手段15は、前述したVmf、CTIn、Chosn1を演算するものである。
【0058】
次に、プラントモデル31の同定について、項を分けて詳述する。
【0059】
1.プラントモデル
燃料性状を推定するためには、エンジンのダイナミクスから燃料挙動のダイナミクスだけを抽出する必要がある。4サイクルエンジン(プラント)は、図3に示すような要素から構成されると考えられる。観測できる入出力は吸入空気量Qaに基づいて基本噴射パルス幅Tpを演算し、これに各種補正を行って求めた実噴射パルス幅CTIn(ただし、始動時は始動時噴射パルス幅TIST、加速時割り込み噴射時はIJSETn)と、A/D変換された排気マニフォールド集合部のA/Fセンサ出力値である。この入出力より得られるプラントモデルは、
1.1:燃料挙動モデル(無駄時間+遅れ系)、
1.2:排気モデル(無駄時間+遅れ系)および
1.3:各種演算および燃焼サイクルに依存する無駄時間
の3つより構成される。
【0060】
1.1燃料挙動モデル
燃料噴射弁から噴射される燃料の挙動は、図4のようにモデル化でき、その数学モデルは次のように表される。
【0061】
【数1】
Ffc=(1−kWW)・Ffi+Ffe
Ffe=e-t/TWW・kWW・Ffi/TWW
GWW(s)=(1−kWW)+kWW/(sTWW+1)
ただし、GWW:燃料挙動の伝達関数、
Ffi:燃料噴射分、
Ffe:燃料蒸発分、
Ffc:シリンダ吸入燃料、
kWW:付着率、
TWW:蒸発の時定数、
この数学モデルは、一つの時定数(TWW)および一つのゲイン(kWW)で表されているが、燃料の振る舞いには、一般的に燃料の付着、蒸発による時定数と、シリンダ吸入遅れによる時定数とがあり、前述の表現でいえば、応答が遅いほうが低周波成分、応答が速いほうが高周波成分である。
【0062】
そこで、応答の異なる2種燃料の挙動に対しては、数1式の数学モデルを並列結合すればよい。このときの数学モデルは次のように表せる。
【0063】
【数2】
ただし、Tsample:サンプル周期(空燃比をサンプリングする周期)、
T1 :低周波成分の時定数、
T2 :高周波成分の時定数、
k1 :低周波成分のゲイン、
k2 :高周波成分のゲイン、
A1 :e-Tsample/T1、
A2 :e-Tsample/T2、
B1 :k1、
B2 :k2、
数2式において、1番目の式(連続値系)を離散値系に変換したものが2番目の式、この2番目の式にz=esTsampleを代入して整理したものが3番目の式である。また、3番目の式をブロック図で示したのが図5である。
【0064】
なお、燃料挙動に伴う無駄時間については燃料挙動の数学モデルに取り込まず、出力信号を時系列的にオフセットさせることで、モデル次数の上昇を抑える(詳細は1.3で後述する)。
【0065】
1.2排気モデル
排気モデルは、図6のように各気筒毎の排気ダイナミクス、排気マニフォールド集合部におけるガス混合ダイナミクス、センサ特性の3要素から構成されると考える。これらを総合すると、「無駄時間+遅れ系」の物理モデルで表すことができる。遅れ系としては、排気ガス輸送遅れ+ガス混合遅れ+センサ応答遅れが考えられ、それぞれが1次遅れ以上のモデルである。ただし、今回はオンボードで(ECM2上で)プラントモデルを同定することを考えると、なるべく高次となることを避ける必要があり、排気モデルとして1つの時定数で代表し、次のように1次遅れ系の数学モデルで記述する。
【0066】
【数3】
Gex(s)=1/(sTex+1)
Gex(z)=(1−e-Tsample/Tex)/(z−e-Tsample/Tex)
Gex(z)=(1−A3)/(z−A3)
ただし、Gex :排気動特性の伝達関数、
Tex :蒸発時定数、
Tsample:サンプル周期、
A3 :e-Tsample/Tex、
数3式においても、1番目の式を離散値系に変換したものが2番目の式、2番目の式にz=esTsampleを代入して整理したものが3番目の式である。
【0067】
また、無駄時間については、燃料挙動モデルと同様に時系列的にオフセットさせることとし、排気モデルに取り込まない(1.3で後述する)。
【0068】
1.3無駄時間モデル
本プラントモデル31は、入力に実噴射パルス幅CTIn、出力にA/Fセンサ3出力電圧読み込み値とおいているため、図3に示したように、実噴射パルス幅CTInを演算してからA/Fセンサ出力値をECM2が読み込むまでにはいくつかの無駄時間(図では「Delay」で表示)が存在する。そこで、改めて入出力間の無駄時間を図7に詳細に示す。
【0069】
ここで、各無駄時間を説明する。
【0070】
1)Delay1:噴射量演算から実噴射タイミングまでのディレイ、
実噴射パルス幅CTInは10ms毎に演算しており、実際に燃料を噴射するタイミングまでは必ずしも毎サイクル同一ではない。そこで演算された噴射タイミング(燃料噴き始め)角度とそのときの回転数より、パルス幅演算タイミングからの時間を算出することにより無駄時間Delay1を求める。
【0071】
2)Delay2:実噴射タイミングからIVO(吸気弁開)までのディレイ、これは、燃料を噴射し、吸気弁が開いて燃料がシリンダに吸入されるまでの時間である。本ディレイは燃料挙動特性により決まり、各運転条件および燃料性状によって設定する。たとえば、市販されている燃料のうち揮発性がほぼ中間の燃料を用いて、エンジン回転数と負荷毎に燃料をステップ的に変化させて実際に応答時間を計測しながら適切な値を設定する。
【0072】
3)Delay3:IVOからEVO(排気弁開)まで(吸入→圧縮→燃焼→排気)のディレイ、
これは、燃料が吸気弁よりシリンダ内に吸入され、燃焼ガスが排気弁から排出されるまでの時間である。このディレイは回転数およびカムプロフィールから求めることができ、設計仕様より求められる。
【0073】
4)Delay4:燃焼ガスが排気弁を出てからA/Fセンサに到達するまでのディレイ、
これは、シリンダ内の燃焼ガスが排気弁より排出されてから排気管を通りA/Fセンサ部に到達するまでの時間である。排気の流速(回転数、負荷等に依存)と排気長さ、A/Fセンサ取り付け位置等によって設定する。なお、計算が複雑になるが、運転条件およびハードスペックから求めてもよい。
【0074】
5)Delay5:センサ応答ディレイ、
これは、A/Fセンサ部にガスが到達してからA/Fセンサが電圧を出力するまでの時間である。A/Fセンサ自体は酸素量が変動すると、数msで反応するが、センサカバーによるガス混合遅れが支配的であり、その時間はセンサカバーの形状により大きく異なる。このため、本ディレイはDelay4と同様にして適切な値を設定する。
【0075】
6)Delay6:センサ出力値をA/D変換しECM2に取り込むまでのディレイ、
これは、A/Fセンサが出力する電圧をECM2に取り込むためにA/D変換を行うことによるディレイである。現在のハード構成ではA/FセンサのA/D変換は2ms毎に行っており、最大無駄時間は2msである。
【0076】
7)Delay7:センサ出力値をメモリにバッファリングするまでのディレイ、
A/Fセンサ出力値を10msでサンプリングする場合に、サンプリングタイミングによっては最大10msの無駄時間が生じるおそれがある。このため2msでA/D変換した値より算出した応答開始からサンプリングタイミングまでの無駄時間をDelay7として算出する。
【0077】
8)Delay8:バラツキ分、
これは、Delay1〜7より求めた無駄時間以外に生じるバラツキ分である。機種間バラツキや適合のバラツキ等が考えられ、運転毎にも不確定であるため、Delay1〜7経過後の空燃比信号の立ち上がりを判定し決定する。
【0078】
これで、各無駄時間の説明を終える。
【0079】
上記のDelay1〜8を分析すると、無駄時間は
▲1▼運転条件により決定する項、
▲2▼演算タイミングより決定する事項、
▲3▼燃料性状等により変動する項
に分類することができる(図8参照)。よって、実際の無駄時間は、次の式で表すことができる。
【0080】
無駄時間
=適合項(運転条件により決定、テーブルやマップの検索)
+演算項(演算タイミングより決定、タイミング信号を等を用いて演算)
+判定項(燃料性状等のバラツキ等により変動するため、随時判定する)
2.プラントモデルの同定
2.1同定するモデルの作成
実際のエンジンは強い非線形性を有するが、本制御では、ある動作点近傍では線形でありかつ時不変である、いわゆる線形時不変システム(LTI:Linear Time-Invariant System)であると仮定する。
【0081】
また、離散時間系LTIシステムを、Z領域ではなく時間領域で入出力を記述するため、シフトオペレータq-1を以下のように定義する。
【0082】
【数4】
q-1x(k)=x(k−1)
ただし、離散時間=kT(T:サンプリング周期、k=0,1,2,・・・)である。
【0083】
これを用いて、離散値系の入力u(t)、出力y(t)のシステム伝達関数を記述すると、
【0084】
【数5】
y(k)=G(q,θ)・u(k)
となる。θはモデルを記述するパラメータにより構成される。しかし、これは理想的な入出力であり、外部からの雑音を考慮すると、
【0085】
【数6】
y(k)=G(q,θ)・u(k)+H(q,θ)・w(k)
と記述できる。ここで、H(q,θ)は雑音モデルであり、一般的な離散時間系LTIシステムは、数6式で表すことができる。同システムのブロック図は図9である。
【0086】
ここで、同システムの伝達関数G(q)は、数2式の3番目の式と数3式の3番目の式の積であり、さらにZ-1で記述されたものをシフトオペレータq-1で記述したものとなって、
【0087】
【数7】
と表せる。このシステムの伝達関数G(q)を、
【0088】
【数8】
G(q,θ)=B(q,θ)/A(q,θ)
の式で定義すると、システムの出力値y(k)は、
【0089】
【数9】
y(k,θ)={B(q,θ)/A(q,θ)}・u(k)+H(q,θ)・w(k)
と表すことができる。このように、同定するモデルとしては、プラントモデルであるG(q)と、雑音モデルであるH(q)を適切な形としたものとを組み合わせたものを採用する。
【0090】
2.2同定手法
数6式で定義した離散時間系LTIシステムにおいて、時刻(k−1)までに測定された入出力データに基づいた出力y(k)の一段先予測値y(k|θ)は、
【0091】
【数10】
の式で表される。これにより、時刻kにおける出力を(k−1)までに取得したデータで記述することができる。
【0092】
予測誤差ε(k|θ)は
【0093】
【数11】
ε(k|θ)=y(k)−y(k|θ)
の式で表すことができる。
【0094】
さて、パラメータ推定のための評価規範JN(θ)として、
【0095】
【数12】
を設定する。ここで、関数l(k,θ,ε(k,θ))は予測誤差ε(k,θ)の大きさを測る任意のスカラ値関数であり、どのようなノルムを選択するかは、同定結果の利用目的に依存する(2乗ノルムや対数尤度など)。このような評価規範を定義することによって、未知パラメータθの推定値(θ(N)とする)が決定される。つまり、
【0096】
【数13】
となるθを求めることである。
【0097】
一般的に同定手法には様々な手法が提案されているが、エンジンのようなものは間欠的なイベント(燃焼サイクル等)であり、非常に非線形性が強い制御対象である。しかしながら、本制御ではアルゴリズムの簡略化のため、動作点周りでは線形時不変(LTI)システムであると仮定している。
【0098】
今回は演算量の少なさ、同定精度、対外乱性能を考慮し、線形モデルの同定手法の代表的なものである「パラメトリックモデル同定であるARXモデルを用いた一括同定手法」を採用する。
【0099】
2.3ARXモデルの同定手法
ARXモデルは式誤差モデルと呼ばれ、次のように差分方程式の右辺に外乱項e(k)(ARXモデルでは白色雑音として仮定しており、w(k)とする)が入っている。
【0100】
【数14】
モデルを記述するパラメータベクトルθは、
【0101】
【数15】
θ=[a1,・・・,ana,b1,・・・,bnb]T
となる。データベクトル(回帰ベクトル)ψ(k)を、
【0102】
【数16】
と定義すると、出力y(k)は次式のように表現できる。
【0103】
【数17】
y(k)=θTψ(k)+w(k)
ARXモデルの一段先予測値y(k|θ)は、数10式より求めると、θに関して線形であり、
【0104】
【数18】
y(k|θ)=θTψ(k)
と表される。このときの予測誤差ε(k,θ)は、
【0105】
【数19】
ε(k,θ)=y(k)−θTψ(k)
と表わすことができる。この線形回帰モデルに対して最小2乗法を適用すると、スカラ値関数l(k,θ,ε(k,θ))は、
【0106】
【数20】
l(k,θ,ε(k,θ))=ε2(k,θ)
となり、パラメータ推定の評価規範JN(θ)は、
【0107】
【数21】
となる。数21式をさらに計算すると、
【0108】
【数22】
【0109】
【数23】
JN(θ)=c(N)−2θTf(N)+θTR(N)θ
とおくことができる。
【0110】
ただし、数23式のc(N)、f(N)、R(N)は次の通りである。
【0111】
【数24】
【0112】
【数25】
【0113】
【数26】
評価規範JN(θ)が最小となるのは、JN(θ)がθに関する二次関数であるため、最高次の係数が正であれば、JN(θ)の微分値がゼロとなるところである。数23式の微分値=0とすると、次の正規方程式(θに関する連立一次方程式)が得られる。
【0114】
【数27】
これより、R(N)が正定値行列であれば、JN(θ)は微分値がゼロのとき最小となり(JN(θ)は下に凸の関数、図10参照)、
【0115】
【数28】
の式によりパラメータθ(N)を推定することができる。以上の同定手順を図11に示す。
【0116】
なお、上記の正定値行列の条件には次の3つがある。
【0117】
1)同定対象がn次の場合は、入力信号u(k)はn個以上の正弦波を含んでいなければならない(ステップ入力信号に十分な周波数成分を含ませる)。
【0118】
2)同定対象は安定である(エンジンは定常では安定系と考えて差し支えない)。
【0119】
3)同定対象は可観測である。すなわち、A(q,θ)とB(q,θ)は共通因子を持たない(本モデルは離散系であるためB(q,θ)のほうが次数が高いが問題なし)。
【0120】
2.4実際のARXモデルの同定
本モデルは、数7式より分母3次、分子3次の離散系モデルであり、
【0121】
【数29】
A(q)=1+a1・q-1+a2・q-2+a3・q-3
【0122】
【数30】
B(q)=b1・q-1+b2・q-2+b3・q-3
と表すことができる。よってパラメータベクトルθおよびデータベクトルψ(k)は、以下のように表すことができる。
【0123】
【数31】
θ=[a1,a2,a3,b1,b2,b3]T
【0124】
【数32】
エンジン回転数が1200rpm時のサンプリング総数NをN=128(1280ms)とすると、数24式〜数26式は、
【0125】
【数33】
【0126】
【数34】
【0127】
【数35】
と表すことができる。
【0128】
2.5プラントモデル同定に必要な入力信号
システムの同定を行うためには、入力信号が、対象のもつ全てのモードを励起している必要がある。つまり、入力信号が多数の周波数成分を含んでいる必要がある。システムの同定においては、理想的には白色性入力が望ましいが、実際には疑似白色2値信号(M系列)が用いられる。しかし、エンジンの壁流応答のようなものでは、有効な周波数帯域は非常に低いところであり(応答が遅い)、M系列のような入力を加えてもほとんど応答波形を得ることができない。そこで、ステップ入力を与えることにより得られる波形(図12参照)をもとに、システムを同定する。なお、ステップ入力のラプラス変換は1/sであるので、周波数ゲインは周波数に対して反比例で減少するため、パワースペクトルより有効な周波数域を決めておく必要がある。
【0129】
2.6実験結果
このようにして求めたパラメータθを用いれば、システムの伝達関数G(q,θ)が定まるので、同定結果と実データを重ねたボード線図を図13に示す(燃料性状が軽質であるほどカットオフ周波数が高くなる傾向がある)。実験結果によれば、吸気ポートに設けたスワールコントロールバルブが開状態、エンジン回転数が1200rpm近傍、冷却水温が40℃近傍かつ低負荷域において±3σで燃料性状の異なる2つのガソリンを分離することができた。
【0130】
これで、項分け説明を終える。
【0131】
次に、ECM2で実行される制御内容を、フローチャートにしたがって説明する。
【0132】
図14は燃料性状を推定するためのもので、一定時間毎(10ms毎)に実行する。ここでは、図14をメインルーチン、図15、図16を図14のサブルーチンとして構成しており、したがって、以下ではメインルーチンの説明途中でサブルーチンのあるステップになると、サブルーチンを説明する。
【0133】
図14においてステップ1では燃料性状判定済みフラグをみる。まだ燃料性状を判定していないときは、ステップ2以降に進む。
【0134】
ステップ2〜6は排気の空燃比(出力データ)をサンプリングする部分である。ステップ2では、A/Fセンサ3で検出される空燃比を読み込み、ステップ3でこの空燃比の読み込み数(サンプリング数)S1とサンプリング総数N(たとえば128)を比較する。S1≦Nであるときは、ステップ4に進んで、空燃比をバッファリングして今回の処理を終了する。S1>Nとなる前はステップ2、4の処理を繰り返す。
【0135】
S1>Nとなったタイミングでステップ3よりステップ5に進む。このとき、バッファにはN個の出力データが格納されている。たとえば、今回値をy(1)に、前回値をy(2)に、2回前の値をy(3)に、・・・、N−1回前の値をy(N)にというようにして、合計でN個の出力データが格納されている。
【0136】
ステップ5、6ではN回前の出力データをバッファから捨て、今回読み込んだ空燃比をバッファリングする。つまり、今回値をy(1)に格納する。なお、過去の出力データは1回ずつ古い側にシフトして格納されることはいうまでもない。
【0137】
図示しないが、ステップ応答時の実噴射パルス幅CTIn(入力データ)も、今回値をu(1)に、前回値をu(2)に、2回前の値をu(3)に、・・・、N−1回前の値をu(N)にというようにして合計でN個の入力データが格納されており、これらの入力データと上記の出力データとは前述した遅れ時間を考慮して対応付けられる。
【0138】
ステップ7、8では、これらN個ずつの入出力データを用いて実噴射パルス幅がステップ変化したときの排気空燃比の応答波形を解析し、その解析結果から使用燃料の燃料性状を推定する。
【0139】
ここで、「排気空燃比の応答波形の解析」とは、前述の表現によれば、入出力データに基づいてARXモデル(プラントモデル)のパラメータθを規範モデルとの予測誤差が最小となるように調整してARXモデルを同定することである。そこでARXモデルの同定を図15のサブルーチンにより、また燃料性状の推定について図16のサブルーチンにより説明する。
【0140】
まず図15において、ステップ21では、バッファにある入出力データ(入力についてu(1)〜u(128)、出力について−y(1)〜−y(128))より数32式を用いてデータベクトルψ(k)を作成する。
【0141】
ステップ22ではこのデータベクトルψ(k)から上記の数33式を用いてR(N)を、またステップ23では出力データy(k)とこのデータベクトルψ(k)から上記の数34式を用いてf(N)を演算し、これらR(N)、f(N)からステップ24において上記の数28式を用いてモデルパラメータθを演算する。
【0142】
次に、図16に移り、このようにして求めたモデルパラメータθからステップ31において離散時間系LTIシステムの伝達関数G(q,θ)を演算する(θから上記の数29式、数30式を用いてA(q)、B(q)を作成し、この2つよりG(q,θ)(=B(q)/A(q))を算出する)。
【0143】
このシステム伝達関数G(q,θ)からステップ32においてARXモデルのカットオフ周波数fcRealを演算する。ステップ33ではこのカットオフ周波数fcRealと規範モデルのカットオフ周波数fcRefを比較する。
【0144】
ここで、基準燃料に重質ガソリンを用いているので、軽質ガソリンが使用されていればfcReal>fcRefとなり、重質ガソリンが使用されているときはfcR eal≦fcRefとなる。したがって、fcReal>fcRefのとき(軽質ガソリンの使用時)はステップ34に進んで燃料性状切換フラグ=1とし、これに対して、fcReal≦fcRefのとき(重質ガソリンの使用時)は、ステップ33よりステップ35に進んで燃料性状切換フラグ=0とする。
【0145】
このようにしていずれの燃料が使用されているのかの判定が終了したら、図14のステップ9に進み、燃料性状切換フラグの値(燃料性状の判定結果)をEEPROMに格納したあと、ステップ10において燃料性状判定済みフラグ=1とする。この燃料性状判定済みフラグ=1の処理により、次回以降は、図14のステップ2以降に進むことができない(燃料性状の判定回数が1回だけとなる)。
【0146】
このようにして燃料性状の判定が可能になると、燃料噴射量の各種補正量や始動時燃料噴射量を燃料性状の違いに応じて与えることができる。これを具体的にKASの場合で説明する。
【0147】
図17は、始動後増量補正係数KASを演算するためのもので、一定時間毎(たとえば10ms毎)に実行する。
【0148】
ステップ41でEEPROMに格納されている燃料性状切換フラグを読み込む。ステップ42ではスタータスイッチをみてこれがONのときは、ステップ43に進み、冷却水温TWと回転数Neを読み込む。このうち冷却水温TWと燃料性状切換フラグの値から、ステップ44、45においてそれぞれ図18、図19を内容とするテーブルを検索することにより始動後増量水温補正値TKAS、第2始動後増量補正係数KASSを、また回転数Neと燃料性状切換フラグの値からステップ46において図20を内容とするテーブルを検索することにより始動後増量回転補正値TNKASを演算し、これらの値を用いステップ47において
【0149】
【数36】
KAS=TKAS×TNKAS+KASS
の式により始動後増量補正係数KASを算出する。
【0150】
ステップ48では、スタータスイッチがOFFになってからの処理に備えるため、TKASの値をTKASn-1に、KASSの値をKASSn-1に移して今回の処理を終了する。TKASn-1、KASSn-1は前回値を保持するためのメモリである。
【0151】
やがてスタータスイッチがOFFになると(始動完爆)、ステップ42よりステップ49以降の減衰操作に進む。
【0152】
ステップ49では、燃料性状切換フラグの値と始動後時間tとから図21を内容とするテーブルを検索することにより始動後増量減少時間割合TMKASを演算し、ステップ50でこの値だけ前回値を減少させた値を今回のTKAS(=TKASn-1−TMKAS)として算出する。スタータスイッチがOFFになった直後はTKAS>0であるので、そのままステップ53に進む。
【0153】
ステップ53〜56はステップ49〜52と同様である。ステップ53で燃料性状判定フラグの値と始動後時間tから図22を内容とするテーブルを検索することにより第2始動後増量減少時間割合TMKASSを演算し、ステップ54でこの値だけ前回値を減少させた値を今回のKASS(=KASSn-1−TMKASS)として算出する。このときも、KASS>0であるので、そのままステップ46に進んでステップ46以降の処理を実行する。
【0154】
次回以降はステップ49、50、53、54を繰り返すことになるので、やがてTKASやKASSが0以下となり、このときはステップ52やステップ56に進んでTKASやKASSを0に制限する。
【0155】
この結果、TKAS、KASSとも、スタータスイッチのOFF時の値を初期値として、スタータスイッチのOFF後に一定の割合で減衰して0になる(ただし、TNKASが一定のとき)。ただし、TKASの初期値のほうがKASSの初期値より大きく、かつTKASの減少時間割合のほうがKASSの減少時間割合より大きい。したがって、TKASとKASSを加算した値であるKASは、スタータスイッチOFF時のTKASの値とKASSの値の合計を初期値として、スタータスイッチOFF後にまず急激な勾配で小さくなり、TKASが0になったタイミングからは緩やかな勾配に切換わって減少していく。
【0156】
この場合に、軽質ガソリンの使用時のほうが重質ガソリンの使用時よりもTKAS、KASSの各初期値(図18、図19のテーブル値)を小さく、かつTKAS、KASSの各減少時間割合(図21、図22のテーブル値)を大きくしているので、軽質ガソリン使用時のKASは、重質ガソリン使用時のKASより小さくなる(図23参照)。つまり、燃料性状の判定を行っていないものでは、軽質ガソリンの使用時にも、重質ガソリンに対してマッチングしたテーブル値を用いることによる空燃比のリッチ化を招くのであるが、このように、軽質ガソリンの使用であることを判定したときは、次回の始動時のKASの演算に際して、軽質ガソリン用のKASを演算することで、軽質ガソリンの使用時にも空燃比がリッチ側に偏ることがなくなるのである。
【0157】
図24のフローチャートは第2実施形態で、第1実施形態の図16に対応する。図16と同一部分には同一のステップ番号を付している。
【0158】
基準燃料に重質ガソリンを用いた第1実施形態に対して、第2実施形態は、市販されている燃料のうち揮発性が悪くもなく良くもないほぼ中間の燃料を基準燃料として、またこの基準燃料よりも揮発性の良い燃料を軽質ガソリン、この逆に基準燃料よりも揮発性の悪い燃料を重質ガソリンとして設定しておき、プラントモデルのカットオフ周波数と規範モデルのカットオフ周波数(つまり基準燃料に対するカットオフ周波数)の差を演算し、この周波数差と許容範囲とを比較して使用されている燃料の燃料性状を判定するようにしたもので、これによって、基準燃料に対するカットオフ周波数がバラツクことがあっても、燃料性状の推定を安定して行うことができる。
【0159】
図24において、図16と相違する部分を主に述べると、ステップ61で
【0160】
【数37】
Δfc=fcReal−fcRef
の式により基準燃料とのカットオフ周波数差Δfcを計算し、この周波数差Δfcの絶対値と許容範囲を定める所定値a(>0)、あるいはΔfcとaをステップ62、63において比較する。|Δfc|≦a(つまり基準燃料が使用されている)であれば、ステップ62よりステップ64に進んで燃料性状切換2フラグ=0とし、|Δfc|>aかつΔfc>a(つまり軽質ガソリンが使用されている)であるときはステップ63よりステップ65に進んで燃料性状切換2フラグ=1とし、それ以外(つまり重質ガソリンが使用されている)のときはステップ63よりステップ66に進んで燃料性状切換2フラグ=2とする。
【0161】
このように第2実施形態では、使用されている燃料が、基準燃料、軽質ガソリン、重質ガソリンのいずれであるかが判定された。
【0162】
ただし、第2実施形態のように、燃料性状判定値が3つの値になると、始動後増量補正係数KASを演算するに際し、図18〜図22に対応して3種類のテーブル値を用意する必要がある。
【0163】
次に、図25、図26のフローチャートは第3実施形態で、それぞれ第1実施形態の図14、図16に対応する。図25において図14と同一部分に、また図26において図16と同一部分に同一のステップ番号を付している。
【0164】
前述の2つの実施形態が燃料性状判定値(つまり燃料性状切換フラグや燃料性状切換2フラグの値)が2値あるいは3値であったのに対して、第3実施形態は、燃料性状推定値を演算し、これをEEPROMに格納するようにしたものである(図25のステップ71、72)。この燃料性状推定値を用いることで、各種噴射量補正量の演算精度を高めることができる。
【0165】
燃料性状推定値の演算について具体的に図26により説明すると、ステッ31、32で第1実施形態と同じにプラントモデルのカットオフ周波数を演算し、その演算したカットオフ周波数からステップ81において図27を内容とするテーブルを検索することにより燃料性状推定値を演算する。プラントモデルのカットオフ周波数と燃料性状推定値との関係は図27のようになるので、同特性を予めマッチングにより定めておけば、プラントモデルのカットオフ周波数から使用燃料の燃料性状を推定できるのである。
【0166】
実施形態では、過渡時に燃料噴射量に対する排気空燃比の応答波形をサンプリングし、これら過渡時データに基づいて、予めECM上に構築したプラントモデルのパラメータを、基準燃料に対するプラントモデルである規範モデルとの予測誤差が最小となるように調整することにより、使用燃料に対するプラントモデルを同定する場合で説明したが、これに限られるものでなく、燃料噴射量に代えて燃料供給量を用いることもできる。また、予測誤差が最小となるように調整するほか、予測誤差が小さくなるように調整することでもかまわない。
【図面の簡単な説明】
【図1】エンジン制御の制御システム図。
【図2】燃料性状の推定に関係する制御システム図。
【図3】エンジンプラントモデルのブロック図。
【図4】燃料挙動のモデル図。
【図5】燃料挙動のパラレルパスブロック図。
【図6】排気モデル図。
【図7】入出力間の無駄時間を表す波形図。
【図8】無駄時間を分類した表図。
【図9】LTIシステムの一般的なブロック図。
【図10】評価関数(評価規範)の特性図。
【図11】ARXモデルの同定手法を示すフローチャート。
【図12】モデル同定に必要な入力信号とその応答を示す波形図。
【図13】同定結果と実データを重ねて示すボード線図。
【図14】燃料性状の推定を説明するためのフローチャート。
【図15】ARXモデルの同定を説明するためのフローチャート。
【図16】燃料性状の切換判定を説明するためのフローチャート。
【図17】始動後増量補正係数KASの演算を説明するためのフローチャート。
【図18】始動後増量水温補正値(初期値)の特性図。
【図19】第2始動後増量補正係数(初期値)の特性図。
【図20】始動後増量回転補正値の特性図。
【図21】始動後増量減少時間割合の特性図。
【図22】第2始動後増量減少時間割合の特性図。
【図23】始動後増量補正係数KASの時系列イメージを示す波形図。
【図24】第2実施形態の燃料性状の切換判定を説明するためのフローチャート。
【図25】第3実施形態の燃料性状の推定を説明するためのフローチャート。
【図26】第3実施形態の燃料性状推定値の演算を説明するためのフローチャート。
【図27】カットオフ周波数に対する燃料性状推定値の特性図。
【図28】第1の発明のクレーム対応図。
【図29】第2の発明のクレーム対応図。
【図30】第3の発明のクレーム対応図。
【符号の説明】
2 ECM
3 A/Fセンサ
7 燃料噴射弁
14 EEPROM
21 プラント同定部
22 燃料性状推定部
23 トリガリング機能
24 コントローラ
31 プラントモデル
37 規範モデル
38 比較手段
Claims (14)
- エンジンの運転条件に応じた燃料供給量を演算する手段と、
この燃料供給量をエンジンに供給する手段と、
エンジンの排気空燃比を検出する手段と、
過渡時に前記燃料供給量を入力、前記排気空燃比を出力として現在の使用燃料での前記燃料供給量に対する排気空燃比の応答波形のデータをサンプリングする手段と、
これら入出力データに基づいて、予め構築したプラントモデルのパラメータを規範モデルとの予測誤差が小さくなるように調整することにより、前記使用燃料に対するプラントモデルを同定する手段と、
この同定されたプラントモデルのカットオフ周波数を演算する手段と、
このプラントモデルのカットオフ周波数と規範モデルのカットオフ周波数とを比較して前記使用燃料の燃料性状を推定する手段と
を設けたことを特徴とする燃料性状検出装置。 - エンジンの運転条件に応じた燃料供給量を演算する手段と、
この燃料供給量をエンジンに供給する手段と、
エンジンの排気空燃比を検出する手段と、
過渡時に前記燃料供給量を入力、前記排気空燃比を出力として現在の使用燃料での前記燃料供給量に対する排気空燃比の応答波形のデータをサンプリングする手段と、
これら入出力データに基づいて、予め構築したプラントモデルのパラメータを規範モデルとの予測誤差が小さくなるように調整することにより、前記使用燃料に対するプラントモデルを同定する手段と、
この同定されたプラントモデルのカットオフ周波数を演算する手段と、
このプラントモデルのカットオフ周波数と規範モデルのカットオフ周波数の差と許容範囲とを比較して前記使用燃料の燃料性状を推定する手段と
を設けたことを特徴とする燃料性状検出装置。 - エンジンの運転条件に応じた燃料供給量を演算する手段と、
この燃料供給量をエンジンに供給する手段と、
エンジンの排気空燃比を検出する手段と、
過渡時に前記燃料供給量を入力、前記排気空燃比を出力として現在の使用燃料での前記燃料供給量に対する排気空燃比の応答波形のデータをサンプリングする手段と、
これら入出力データに基づいて、予め構築したプラントモデルのパラメータを規範モデルとの予測誤差が小さくなるように調整することにより、前記使用燃料に対するプラントモデルを同定する手段と、
この同定されたプラントモデルのカットオフ周波数を演算する手段と、
カットオフ周波数に対する燃料性状推定値の特性を予め設定する手段と、
前記演算されたカットオフ周波数からこの特性を検索することにより燃料性状推定値を演算する手段と
を設けたことを特徴とする燃料性状検出装置。 - 前記燃料性状の推定結果を不揮発性メモリに記憶させておくことを特徴とする請求項1または2に記載の燃料性状検出装置。
- 前記燃料性状推定値を不揮発性メモリに記憶させておくことを特徴とする請求項3に記載の燃料性状検出装置。
- 前記燃料性状の推定は、前記規範モデルを基準燃料に対してマッチングした場合に、前記同定したプラントモデルのカットオフ周波数が規範モデルのカットオフ周波数より高いとき、前記基準燃料よりも軽質であると推定することであることを特徴とする請求項1に記載の燃料性状検出装置。
- 前記燃料性状の推定は、前記規範モデルを基準燃料に対してマッチングした場合に、前記同定したプラントモデルのカットオフ周波数と規範モデルのカットオフ周波数の差が許容範囲外であり、かつプラントモデルのカットオフ周波数が規範モデルのカットオフ周波数よりも大きいとき、前記基準燃料よりも軽質であると推定することであることを特徴とする請求項2に記載の燃料性状検出装置。
- 前記予測誤差が小さくなるように調整することは、予測誤差が最小となるように調整することであることを特徴とする請求項1から7までのいずれか一つに記載の燃料性状検出装置。
- 前記プラントモデルを、燃料挙動モデル、排気モデルおよび無駄時間の3つより構成することを特徴とする請求項1から8までのいずれか一つに記載の燃料性状検出装置。
- 前記燃料挙動モデルを分母2次、分子2次の2次遅れ系モデルとすることを特徴とする請求項9に記載の燃料性状検出装置。
- 前記排気モデルを分母1次の1次遅れ系モデルとすることを特徴とする請求項9に記載の燃料性状検出装置。
- 前記無駄時間を適合項と演算項と判定項の和とすることを特徴とする請求項9に記載の燃料性状検出装置。
- 前記同定の手法はARXモデルを用いた一括処理最小2乗法であることを特徴とする請求項1から12までのいずれか一つに記載の燃料性状検出装置。
- 前記燃料性状の推定は運転の開始から停止までの間に1回だけであることを特徴とする請求項1から13までのいずれか一つに記載の燃料性状検出装置。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP09870999A JP3772584B2 (ja) | 1999-04-06 | 1999-04-06 | 燃料性状検出装置 |
US09/545,166 US6363313B1 (en) | 1999-04-06 | 2000-04-06 | Fuel property detecting system |
DE10017100A DE10017100B4 (de) | 1999-04-06 | 2000-04-06 | System zur Erfassung einer Kraftstoffeigenschaft |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP09870999A JP3772584B2 (ja) | 1999-04-06 | 1999-04-06 | 燃料性状検出装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2000291484A JP2000291484A (ja) | 2000-10-17 |
JP3772584B2 true JP3772584B2 (ja) | 2006-05-10 |
Family
ID=14227051
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP09870999A Expired - Fee Related JP3772584B2 (ja) | 1999-04-06 | 1999-04-06 | 燃料性状検出装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP3772584B2 (ja) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6619277B2 (en) | 2001-07-12 | 2003-09-16 | Nissan Motor Co., Ltd. | Engine air-fuel ratio control |
JP3965947B2 (ja) | 2001-07-25 | 2007-08-29 | 日産自動車株式会社 | エンジンの空燃比制御装置 |
CN102472176B (zh) | 2010-04-05 | 2013-06-19 | 丰田自动车株式会社 | 内燃机的控制装置 |
-
1999
- 1999-04-06 JP JP09870999A patent/JP3772584B2/ja not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JP2000291484A (ja) | 2000-10-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6327850B1 (en) | Air-fuel ratio control apparatus for multicylinder internal combustion engine | |
JPH01216047A (ja) | エンジンの空燃比制御方法および装置 | |
JP5273183B2 (ja) | 内燃機関の制御装置 | |
US20130220284A1 (en) | Air amount estimating apparatus for internal combustion engine with supercharger | |
JP3768780B2 (ja) | 内燃機関の空燃比制御装置 | |
JPS63314339A (ja) | 空燃比制御装置 | |
JP3039162B2 (ja) | 内燃機関の空燃比制御装置 | |
US10240516B2 (en) | System for wastegate control | |
CN101187341A (zh) | 内燃机的空气量运算装置及燃料控制装置 | |
US9644549B2 (en) | Control apparatus for internal combustion engine | |
JP3832137B2 (ja) | エンジンの燃料供給制御装置 | |
JP3772584B2 (ja) | 燃料性状検出装置 | |
JP3832136B2 (ja) | 燃料性状検出装置 | |
KR920002456B1 (ko) | 엔진의 연료제어장치 | |
JP3697934B2 (ja) | 燃料性状検出装置およびエンジンの燃料供給制御装置 | |
JP2001132521A (ja) | 内燃機関の圧力検出器異常時制御方法 | |
JPH06108901A (ja) | 内燃機関の空燃比制御装置 | |
JP4385542B2 (ja) | 内燃機関の空燃比制御装置 | |
Deur et al. | Influence of intake manifold heat transfer effects on accuracy of si engine air-charge prediction | |
JP2002089338A (ja) | 内燃機関の燃料噴射制御装置 | |
JPS62255554A (ja) | エンジン制御装置 | |
JP2528279B2 (ja) | 内燃機関の電子制御燃料噴射装置 | |
JPH0392567A (ja) | スロットル開度算出方法 | |
JPH08128348A (ja) | エンジンの制御装置 | |
JPH0763101A (ja) | 内燃機関の燃料噴射量制御装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20040819 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20050215 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20050411 |
|
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: 20060124 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20060206 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090224 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100224 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100224 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110224 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120224 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120224 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130224 Year of fee payment: 7 |
|
LAPS | Cancellation because of no payment of annual fees |