(57)【要約】【課題】ノッキング回避制御の精度を向上させ、機関の効率を維持可能な燃焼制御を行う。【解決手段】各種センサからの信号に基づき機関の運転状態を検出し、これに基づいて、既燃部と未燃部との燃焼室２領域サイクルシュミレーションによりエンドガスの自着火発生時期および自着火による発熱量を予測し、予測された自着火発生時期および自着火による発熱量からノッキング強度を推定して、このノッキング強度が所定値以下となるよう、機関の燃焼制御を行う。【選択図】  図２
2004332584.tif 000002
【特許請求の範囲】
【請求項１】
火花点火式内燃機関において、
機関の運転状態を検出する手段と、
検出された機関運転状態に基づいてエンドガスの自着火発生時期および自着火による発熱量を予測する手段と、
機関の燃焼を制御する手段と、を備え、
前記予測手段によるエンドガスの自着火発生時期および自着火による発熱量の関係が、ノッキング強度が所定値以下となるよう、前記燃焼制御手段により機関の燃焼制御を行うことを特徴とする火花点火式内燃機関の燃焼制御装置。
【請求項２】
火花点火式内燃機関において、
機関の運転状態を検出する手段と、
検出された機関運転状態に基づいてエンドガスの自着火発生時期および自着火による発熱量を予測する手段と、
予測された自着火発生時期および自着火による発熱量からノッキング強度を推定する手段と、
機関の燃焼を制御する手段と、を備え、
前記ノッキング強度推定手段によるノッキング強度が所定値以下となるよう、前記燃焼制御手段により機関の燃焼制御を行うことを特徴とする火花点火式内燃機関の燃焼制御装置。
【請求項３】
前記ノッキング強度は、エンドガスの自着火による発熱量が大きい程強く、自着火発生時期が進角側である程強いと推定することを特徴とする請求項１または請求項２に記載の火花点火式内燃機関の燃焼制御装置。
【請求項４】
前記ノッキング強度は、機関回転速度が高回転側である程強いと判定することを特徴とする請求項３記載の火花点火式内燃機関の燃焼制御装置。
【請求項５】
前記ノッキング強度に関する所定値は、トレースノック相当のノッキング強度であることを特徴とする請求項１〜請求項４のいずれか１つに記載の火花点火式内燃機関の燃焼制御装置。
【請求項６】
前記燃焼制御手段は、点火時期を制御する手段であることを特徴とする請求項１〜請求項５のいずれか１つに記載の火花点火式内燃機関の燃焼制御装置。
【請求項７】
前記エンドガスの自着火発生時期および自着火による発熱量を予測する手段は、エンドガスの着火遅れ時間を予測する手段であることを特徴とする請求項１〜請求項６のいずれか１つに記載の火花点火式内燃機関の燃焼制御装置。
【請求項８】
前記エンドガスの自着火発生時期および自着火による発熱量を予測する手段は、着火遅れ時間の逆数を積分することでエンドガスの自着火を予測する手段であることを特徴とする請求項１〜請求項６のいずれか１つに記載の火花点火式内燃機関の燃焼制御装置。
【請求項９】
前記エンドガスの自着火発生時期および自着火による発熱量を予測する手段は、化学素反応モデルに基づく自着火予測手段であることを特徴とする請求項１〜請求項６のいずれか１つに記載の火花点火式内燃機関の燃焼制御装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、火花点火式内燃機関の燃焼制御装置に関する。
【０００２】
【従来の技術】
従来、火花点火式内燃機関において、圧縮比を大きくすることが熱効率向上のために有効であるが、ノッキングの発生はその性能向上の阻害要因の１つであることが知られている。
【０００３】
そこで、機関にノックセンサーを取り付けてノッキングの発生を検知し、点火時期等に対してフィードバック制御を行うことで、機関の効率を出来るだけ高く維持しつつ、ノッキングの発生を回避する試みが多くなされてきており、自動車用の内燃機関においては、既に多くの機関にこの仕組みが装備されている。
【０００４】
しかしながら、ノッキングの発生を検知してから回避制御を行うだけでなく、ノッキングの発生を計算モデルに基づいて事前に予測可能であれば、より確実にノッキングの回避を行うことが可能となるはずである。
【０００５】
また点火時期の制御などでは、機関の回転速度（回転数）や出力に対して点火時期をマップ上に設定し、上記のようなノッキング検知による補正を行う手法が一般的であった。この場合、予め実機において膨大な量の実験を行い、設定すべき点火時期を求めておく必要がある。このため、近年のように、可変動弁機構やＥＧＲ等の可変機構が一般化する中では、点火時期の設定のための実験量もさらに拡大し、内燃機関の開発コストや開発機関の増大を招いているという問題があった。
【０００６】
従って、ノッキングの発生を計算モデルに基づいて予測制御することは、点火時期設定のための実験の大部分もしくは少なくとも一部を省略することに繋がり、開発効率上のメリットがある。
【０００７】
このようなノッキングの発生を予測して、これを回避しつつ高効率を維持するべく燃焼制御を行う内燃機関の燃焼制御装置として、特許文献１には、内燃機関の未燃ガスの自着火、すなわちノッキングの発生時期を予測し、ノッキングが発生しないよう、内燃機関の燃焼を制御することが開示されている。
【０００８】
【特許文献１】
特開平７−３３２１４９号公報
【０００９】
【発明が解決しようとする課題】
ここでノッキングの発生は、その強度が大きい場合、機関の耐久性に問題が生じることや、不快な打音を生じることがあるが、軽度のノッキングではこれらの問題を特に引き起こさず、いわゆるトレースノック相当であれば、これに点火時期を設定して機関を継続的に運転可能である。すなわち、ノッキングの発生を予測してこれを回避しようとする場合、ノッキング強度がトレースノック相当になるまで、点火時期を遅角するなどの制御を行えばよい。
【００１０】
しかしながら、特許文献１においては、未燃ガスの自着火が発生しなくなることがノッキング回避の条件とされており、ノッキング強度に関する考え方が導入されておらず、必要以上の点火時期遅角により、本来不要な機関の効率低下を生じるおそれがある。これは、未燃ガスの自着火の予測と、実際のノッキング現象との相関について知見が不足していたことに一因がある。
【００１１】
本発明は係る問題点に鑑みなされたもので、未燃ガスの自着火を予測してノッキング回避制御を行う場合に、ノッキング回避制御の精度を向上させ、機関の効率を維持可能な燃焼制御を行うことを目的とする。
【００１２】
【課題を解決するための手段】
そのため本発明では、火花点火式内燃機関の運転状態を検出し、これに基づいてエンドガスの自着火発生時期および自着火による発熱量を予測し、これらの自着火発生時期および自着火による発熱量の関係が、ノッキング強度が所定値以下となるよう、機関の燃焼制御を行う。
【００１３】
また本発明では、自着火発生時期および自着火による発熱量からノッキング強度を推定して、このノッキング強度が所定値以下となるよう、機関の燃焼制御を行う。
【００１４】
【発明の効果】
本発明によれば、エンドガスの自着火発生時期および自着火による発熱量の関係に基づいて、ノッキング強度が所定値以下となるように燃焼制御されるため、ノッキング強度を考慮しつつ機関の熱効率を良好に維持することができる。
【００１５】
また本発明によれば、エンドガスの自着火発生時期および自着火による発熱量からノッキング強度を推定することで、より確実にノッキング強度が所定値以下となるように燃焼制御可能となる。
【００１６】
【発明の実施の形態】
以下、図面に基づき、本発明の実施形態について説明する。
図１は、本発明に係る火花点火式内燃機関の第１の実施形態を示す構成図である。
【００１７】
内燃機関は、シリンダヘッド１、シリンダブロック２、及びピストン３から画成される燃焼室４を有し、吸気バルブ５を介して吸気ポート７から新気を導入する。
【００１８】
吸気バルブ５は、バルブタイミング可変機構５ｂによって、その開閉タイミングを変更可能である。
吸気ポート７には、燃料噴射弁９が燃焼室４を臨んで配置され、噴射された燃料は新気と混合しながら燃焼室４内に吸入され、燃焼に適した可燃混合気を形成する。そして、燃焼室４内に形成された混合気は、点火プラグ１０により点火・燃焼され、排気バルブ６を介して排気ポート８から排出される。
【００１９】
内燃機関の運転状態は、エンジンコントロールユニット（以下「ＥＣＵ」と称する）１１によって統合的に制御される。ＥＣＵ１１には、アクセル開度センサ１２、水温センサ１３、及びクランク角センサ１４等から各種信号が入力され、これらの信号に基づいてＥＣＵ１１の内部で必要な処理・演算を行い、燃料噴射弁９、点火プラグ１０、及びバルブタイミング可変機構５ｂ等を制御する。なお、後述する燃焼制御に必要なエンドガスの自着火予測等もＥＣＵ１１内で行われる。
【００２０】
図２は、ＥＣＵ１１内において、機関の運転状態の検出から燃焼制御までを模式的に示した図である。
先ず、各種センサからの信号が運転状態検出手段に読み込まれ、これらの信号に基づいてバルブタイミング及びこれに伴う機関の有効圧縮比などが推定される。これと同時に、吸気バルブ５の閉時期における燃焼室４内の温度、圧力、空燃比、及び残留ガス率などが、エアフローメータ、吸気温センサ、スロットル開度センサ等の信号に基づいて推定される。これらの推定値は、後述する燃焼室２領域サイクルシミュレーションにおける燃焼予測の初期値として使われる。
【００２１】
２領域サイクルシュミレーションでは、エンドガス（未燃ガス）の温度・圧力履歴に自着火モデルを組み込むことにより、エンドガスの自着火発生時期および自着火による発熱量を予測する。
【００２２】
２領域サイクルシュミレーション内で予測された自着火発生時期および自着火による発熱量に基づいて、ノッキング強度を推定する。このノッキング強度に応じてＭＢＴ（Ｍｉｎｉｍｕｍ  ａｄｖａｎｃｅ  ｆｏｒ  ｔｈｅ  Ｂｅｓｔ  Ｔｏｒｑｕｅ）、すなわち出力トルク及び燃料消費率が最良となる点火時期と、ノッキング発生点であるノックトレース点とを演算する。
【００２３】
そして、ＭＢＴ及びトレースノック点に基づいて、燃焼制御手段は、トレースノック相当のノッキング強度以下となるように、例えば点火時期を制御する等の制御により燃焼を制御する。一方、ＭＢＴ及びトレースノック点は、２領域サイクルシュミレーションにおいて熱発生パターンの設定に利用されてフィードバックされる。
【００２４】
ここで、筒内の燃焼状態は、燃焼室４を、火炎面を境として未燃部と既燃部とに分割した、いわゆる２領域モデルを基本としたサイクルシミュレーションによって予測する。
【００２５】
図３は、２領域モデルの概略図であり、燃焼室４を上から見た状態を示す図である。図示の通り、燃焼室４（筒内）を既燃部（Ｂｕｒｎｅｄ  Ｚｏｎｅ）の領域と、未燃部（Ｕｎｂｕｒｎｅｄ  Ｚｏｎｅ）すなわちエンドガスの存在する領域との２つの領域に分割している。なお、図中の太線の円は火炎面を示し、×印は点火プラグ１０の位置を示している。
【００２６】
簡略化した燃焼室４の形状において、点火プラグ１０の位置から球面状に火炎（Ｆｌａｍｅ）が伝播するものと仮定する。火炎の進行に伴って、質量およびエンタルピが未燃部から既燃部へと移動し、既燃部における発熱・膨張によって未燃部が圧縮される。そして、エンドガス（未燃部のガス）の温度・圧力により自己着火を予測するため、未燃部には自着火モデルを組み込み、火炎の進行が完了する前に未燃部において自着火が発生すればノッキングと判定し、自着火が発生した時期と自着火による発熱量とを出力する（図２参照）。
【００２７】
次に、２領域モデルにおいて燃焼予測計算を行う手法について説明する。既燃部及び未燃部のエネルギー保存式は、単位時間で見ると以下のようにそれぞれ記述される。
【００２８】
【数１】
2004332584.tif 000003
【００２９】
【数２】
2004332584.tif 000004
【００３０】
ここで、ｍは質量、ｘは質量割合、ｕは比内部エネルギー（内部エネルギーの微係数）、ｐは圧力、Ｖは容積、ｄｘｂは単位時間当たりに未燃部から既燃部へ移動する質量割合、ｈは比エンタルピ、及びＱは気体の授受する熱量を示している。なお、添え字ｂ，ｕはそれぞれ既燃部、未燃部を表す。第３項（ｍｄｘｂｈｕ）は火炎面を通じて未燃部から既燃部へと移動するエンタルピを示し、第４項（ｄＱｂ，ｗａｌｌ、ｄＱｕ，ｗａｌｌ）は、既燃部・未燃部の燃焼室壁面からの熱伝達を示している。本モデルは２領域モデルで燃焼室４は密封であり、時間毎に容積Ｖが保存されるため、既燃部・未燃部の体積には、以下の関係がある。
【００３１】
【数３】
2004332584.tif 000005
【００３２】
なお燃焼室４の全体の容積Ｖは、クランク角度θの関数として与えられる。また、既燃部・未燃部のそれぞれについて、以下の状態方程式が成り立つ。
【００３３】
【数４】
2004332584.tif 000006
【００３４】
【数５】
2004332584.tif 000007
【００３５】
ここで、Ｒはガス定数、Ｔは温度を示している。既燃部の質量割合ｘｂは、実験的に得られた熱発生パターンをテーブルにて与えてもよいが、本実施形態では以下のＷｉｅｂｅ関数にて与えることとする。
【００３６】
【数６】
2004332584.tif 000008
【００３７】
ここで、θ０は燃焼開始時期のクランク角度、Δθは燃焼開始から終了までのクランク角度（期間）、ａ，ｍはモデル定数を示している。燃焼室壁面からの伝熱については、筒内のガスが断熱状態にあるものと仮定して、（１）、（２）式の第４項（ｄＱｂ，ｗａｌｌ、ｄＱｕ，ｗａｌｌ）はゼロとして扱う。そして、上記（１）〜（６）式を連立し、機関運転状態検出手段による情報から推定された各状態値を初期値として、クランク角θに応じて方程式を解いていくことで、燃焼室４内における既燃部・未燃部の状態変化を予測する。なお、燃焼室壁面からの伝達を、後述する実施形態において説明するように内燃機関を対象とした伝熱モデルを利用して与える方法もある。
【００３８】
ノッキングの発生の有無については、前述の通り、エンドガスすなわち未燃部のガスに対して着火モデルを適用することで予測する（図２参照）。そして、所定の時期におけるエンドガスの着火遅れ時間τを予測し、この着火遅れ時間から自着火発生時期と自着火による発熱量とを予測する。エンドガスの着火遅れ時間τを表すモデルとして、次式に示すものが知られている。
【００３９】
【数７】
2004332584.tif 000009
【００４０】
ここで、Ａ，Ｂ，ｎはモデル定数を示している。（７）式は、エンドガスが所定圧力・温度状態にあるときの着火遅れ時間τを与えるものであり、この式の適用に際しては、エンドガスが既に筒内圧の上昇に伴う圧縮により、ある程度高温・高圧の状態にあることが予測精度上必要である。このため、（１）〜（６）式をクランク角θに応じて解いていく過程で、筒内圧力が最も高くなるθｐｍａｘにおいて、（７）式を適用することが最も妥当である。
【００４１】
そして、クランク角θｐｍａｘにおいて求められた着火遅れ時間τを、時間単位からクランク角単位に変換し、θｐｍａｘに加算すればエンドガス自着火発生時期の予測値θｋｎｏｃｋが得られる。同時に、（６）式にθｋｎｏｃｋを代入することで、その時点での熱発生割合が得られるため、１からこれを減算することで、θｋｎｏｃｋにおける未燃ガス割合、すなわち自着火による発熱量（自着火によって燃焼する発熱割合）の予測値Ｘａｕｔｏ−ｉｇｎが得られる。
【００４２】
ここで、自着火発生時期の予測値θｋｎｏｃｋと自着火による発熱量の予測値Ｘａｕｔｏ−ｉｇｎとの組み合わせに対するノッキング強度について説明する。
従来、点火時期を進角した（上死点に近づけた）場合、エンドガスの自着火による発熱量が大きくなり、ノッキング強度も強くなることが知られており、図４に示す引用図（ＳＡＥ  ８８０１９６（１９８８））にもそのことが示されている。なお図４（イ）の横軸はクランク角θ（ｄｅｇ．ＡＢＤＣ）、縦軸は通常の燃焼質量割合（ＮＯＲＭＡＬＩＺＥＤ  ＢＵＲＮＥＤ  ＭＡＳＳ  ＲＡＴＥ）、（ロ）の横軸はクランク角θ、縦軸は自着火混合部分（ＡＵＴＯ−ＩＧＮＩＴＥＤ  ＭＩＸＴＵＲＥ  ＦＲＡＣＴＩＯＮ）を示している。
【００４３】
しかしながら、図４に示されている知見は、同じ運転条件下において点火時期のみを変更することで得られたものであり、例えば燃料のオクタン価が変わった場合や、可変バルブタイミング機構によって実質的な圧縮比が変わった場合においても、普遍的に成立する関係であるか否かは不明であった。
【００４４】
そこで本発明では、様々な運転条件下においてノッキングに関する実験を実施し、その結果を解析することで、自着火発生時期と自着火による発熱量とに対するノッキング強度の関係について明らかにすることを試みた。
【００４５】
図５は、図４と同様の実験についてサイクル毎の実験解析結果を示す図であり、横軸は上死点（ＴＤＣ）からのクランク角（ＣＡ  ａｔ  Ｋｎｏｃｋ  Ｏｃｃｕｒｒｅｎｃｅ）θｋｎｏｃｋ（ｄｅｇ．ＡＴＤＣ）、すなわち自着火発生時期、縦軸は自着火混合部分（Ａｕｔｏ−Ｉｇｎｉｔｅｄ  Ｍｉｘｔｕｒｅ  Ｆｒａｃｔｉｏｎ）（％）、すなわち自着火による発熱量を示している。なお、図中の○印は燃料のオクタン価（ＲＯＮ）が９０の場合の結果、◇印は燃料のオクタン価が１００の場合の結果を示している。
【００４６】
オクタン価９０の燃料を使用した実験では、図４と定性的・定量的に同様の傾向を示していることが解る。一方、同条件においてオクタン価１００の燃料を使用した実験では、同じノッキング強度においても、自着火の発生時期および自着火による発熱量が共にオクタン価９０の場合とは大きく異なっていることが解る。すなわち、オクタン価が高ければ自着火の発生時期が上死点に近づき、自着火による発熱量が少なくなり、ノッキングの発生が起こり難くなることを示している。
【００４７】
図６は、図５におけるサイクル毎の実験解析結果をサイクル数に対して平均し、さらに燃料や圧縮比の異なる条件を追加して示す図であり、横軸は上死点（ＴＤＣ）からのクランク角（ＣＡ  ａｔ  Ｋｎｏｃｋ  Ｏｃｃｕｒｒｅｎｃｅ）θｋｎｏｃｋ（ｄｅｇ．ＡＴＤＣ）、すなわち自着火発生時期、縦軸は自着火混合部分（Ａｕｔｏ−Ｉｇｎｉｔｅｄ  Ｍｉｘｔｕｒｅ  Ｆｒａｃｔｉｏｎ）（％）、すなわち自着火による発熱量を示している。
【００４８】
同ノッキング強度においてエンドガスの自着火が進角側で発生する程、その時の自着火による発熱量が小さいことが解る。すなわち、ノッキング強度は、エンドガスの自着火による発熱量が大きいほど強く、エンドガスの自着火が進角側で発生する場合、僅かな発熱量であっても比較的強度の強いノッキングが発生するということであり、エンドガスの自着火（発生時期、発熱量）の予測からノッキングの発生を予測する場合には、この関係を考慮することが重要となる。
【００４９】
図７は、１２００ｒｐｍ及び２４００ｒｐｍについて、図６におけるノッキング発生点、すなわちトレースノック点だけを示し、それぞれ指数関数による近似を行った図であり、横軸は上死点（ＴＤＣ）からのクランク角（ＣＡ  ａｔ  Ｋｎｏｃｋ  Ｏｃｃｕｒｒｅｎｃｅ）θｋｎｏｃｋ（ｄｅｇ．ＡＴＤＣ）、すなわち自着火発生時期、縦軸は自着火混合部分（Ａｕｔｏ−Ｉｇｎｉｔｅｄ  Ｍｉｘｔｕｒｅ  Ｆｒａｃｔｉｏｎ）（％）、すなわち自着火による発熱量を示している。なお、◆は１２００ｒｐｍ、◇は２４００ｒｐｍの場合を示している。
【００５０】
自着火の発生時期と自着火による発熱量との相関曲線が、指数関数において良く近似していることが解る。そして、２４００ｒｐｍにおける近似曲線は、図の下側の方に位置しており、高回転側ではより少ない発熱量で、ノッキング強度が同じトレースノックを生じていることが解る。すなわち、回転速度（回転数）が高回転である程、ノッキング強度が強いことを示している。なお、１２００ｒｐｍと２４００ｒｐｍとの中間回転速度における相関曲線は、両者の内挿によって得ることができる。なお、異なるノッキング強度に対する相関曲線も求めることができる。
【００５１】
以上の説明のように、本実施形態では実機における様々な条件下のノッキング現象を解析した結果、エンドガスの自着火発生時期および自着火による発熱量の組み合わせと、ノッキング強度との間に相関があることを見出すことができる。この相関を、エンドガスの自着火発生時期の予測値θｋｎｏｃｋおよび自着火による発熱量の予測値Ｘａｕｔｏ−ｉｇｎに対して適用することで、サイクルシミュレーションの予測値からノッキング強度の予測値を得る。
【００５２】
実際の火花点火式内燃機関において、出力トルク及び燃料消費率が最良となる点火時期ＭＢＴより進角側にノッキング限界が存在する場合は、通常ノッキング回避制御を行う必要は無く、ＭＢＴに点火時期を設定することが熱効率の観点から最も適切である。
【００５３】
ＭＢＴを基準とした設定点火時期は、機関回転速度と負荷とに対するマップ上に用意されることが通常であるが、ＭＢＴに相当する点火時期も、ノッキング限界と同様に機関の運転状況に応じて変化するため、設定点火時期をマップで用意することには限界がある。
【００５４】
そこで本制御装置においては、サイクルシミュレーションによりノッキング強度を予測しており（図２参照）、この時のトルクについて容易に同時演算することが可能であることから、ＭＢＴについてもノッキング強度と同時に予測して求めることが合理的である。
【００５５】
すなわち、Ｗｉｅｂｅ関数にて表現された熱発生パターンに対するノッキング強度及びトルクに関する予測を行い（図２参照）、熱発生パターンの時期をずらしつつ、これを数回〜数十回繰り返すことで、トレースノック相当となる点火時期と、ＭＢＴとなる点火時期とを知ることができる。
【００５６】
図８は、この複数回の予測から、トレースノック相当となる点火時期を得る手法を示した図であり、（イ）の横軸は燃焼開始時期、すなわち点火時期（Ｓｔａｒｔ  ｏｆ  Ｃｏｍｂｕｓｔｉｏｎ）θ０（ｄｅｇ．ＡＴＤＣ）、縦軸は自着火混合部分（Ａｕｔｏ−Ｉｇｎｉｔｅｄ  Ｍｉｘｔｕｒｅ  Ｆｒａｃｔｉｏｎ）（％）、すなわち自着火による発熱量を示し、（ロ）の横軸は上死点（ＴＤＣ）からのクランク角（ＣＡ  ａｔ  Ｋｎｏｃｋ  Ｏｃｃｕｒｒｅｎｃｅ）θｋｎｏｃｋ（ｄｅｇ．ＡＴＤＣ）、すなわち自着火発生時期、縦軸は自着火による発熱量を示している。
【００５７】
図８（ロ）に左上から右下にかけて示される曲線αは、サイクルシミュレーションによる自着火発生時期と自着火による発熱量とをプロットしたものを結んだ曲線である。一方、左下から右上に描かれる曲線βは、トレースノック相当を示す相関曲線（ｙ＝０．８２６ｅｘｐ（０．１０１ｘ）で表す）であり、両曲線α、βの交点がこの運転条件下におけるトレースノック点として求められる。
【００５８】
曲線α、βの交点における自着による発熱量に基づいて、図８（イ）に示すように、横軸を点火時期とした図において交点となる点により求めることで、トレースノック点火時期を求めることができる。ＭＢＴについては、単純に計算範囲内で、トルクが最大となる点火時期を求めればよい。
【００５９】
上記の手順にて求められたトレースノック点火時期及びＭＢＴ点火時期について、遅角側にあるものを設定点火時期として選択し、燃焼制御手段に対して設定する（図２参照）。これにより、機関の熱効率を出来るだけ高く維持しつつ、機関の耐久性に問題を与える可能性のある強いノッキングの発生を回避することが可能となる。
【００６０】
以上の説明のように、サイクルシミュレーションによりエンドガスの自着火発生時期と自着火による発熱量とを予測し、両者の関係からノッキング強度を予測して、これが所定値以下、すなわちトレースノック相当以下となるように点火時期を制御することで、過剰な点火時期遅角制御による本来不必要な熱効率低下を生じることなく、機関の耐久性を維持して運転することが可能となる。
【００６１】
本実施形態によれば、火花点火式内燃機関において、機関の運転状態を検出する手段と、検出された機関運転状態に基づいてエンドガスの自着火発生時期および自着火による発熱量を予測する手段と、機関の燃焼を制御する手段と、を備え、予測手段によるエンドガスの自着火発生時期および自着火による発熱量の関係が、ノッキング強度が所定値以下となるよう、燃焼制御手段により機関の燃焼制御を行う。このため、予測されたエンドガスの自着火発生時期および自着火による発熱量の関係に基づいて、ノッキング強度が所定値以下となるように燃焼制御を行うことができ、ノッキングの強度を考慮しつつ機関の熱効率を良好に維持することができる。
【００６２】
また本実施形態によれば、予測された自着火発生時期および自着火による発熱量からノッキング強度を推定する手段を備え、ノッキング強度推定手段によるノッキング強度が所定値以下となるよう、燃焼制御手段により機関の燃焼制御を行う。このため、エンドガスの自着火発生時期と自着火による発熱量とからノッキング強度を推定することで、より確実にノッキング強度が所定値以下となるように燃焼制御可能となる。
【００６３】
また本実施形態によれば、ノッキング強度は、エンドガスの自着火による発熱量が大きい程強く、自着火発生時期が進角側である程強いと判定する。このため、エンドガスの自着火による発熱量と自着火発生時期とノッキング強度との関係が提供され、燃焼制御の実現がより容易になる。
【００６４】
また本実施形態によれば、ノッキング強度は、機関回転速度（回転数）が高回転側である程強いと判定する。このため、エンドガスの自着火発生時期と自着火による発熱量とノッキング強度との関係が、機関回転速度に対しても提供され、更に燃焼制御の実現がより容易になる。
【００６５】
また本実施形態によれば、ノッキング強度に関する所定値は、トレースノック相当のノッキング強度である。このため、トレースノック相当とするという指針を提供でき、機関の耐久性上の問題を回避できるだけでなく、燃焼制御の実現がより容易になる。
【００６６】
また本実施形態によれば、燃焼制御手段は、点火時期を制御する手段である。このため、ノッキングの発生を予測してから実際の燃焼制御を実施するまでの時間が短くて済み、確実にノッキング回避制御が可能となる。
【００６７】
また本実施形態によれば、エンドガスの自着火発生時期および自着火による発熱量を予測する手段は、エンドガスの着火遅れ時間を予測する手段である。このため、予測手段が簡便であり、簡単な構成で実現できる。
【００６８】
また本実施形態によれば、エンドガスの自着火発生時期および自着火による発熱量を予測する手段は、エンドガスの着火遅れ時間を予測する手段である。このため、予測手段を簡便で簡単な構成で実現できる。
【００６９】
次に、第２の実施形態に係る火花点火式内燃機関の燃焼制御装置ついて説明する。本実施形態においては、図２中に示される自着火モデルについて、第１の実施形態とは異なるモデルを採用するものである。
【００７０】
第１の実施形態においては、（７）式に示される着火遅れ時間τを与える式に基づき、ある時点での未燃ガスの温度・圧力から着火遅れ時間τを求め、その遅れ時間τの経過後に自着火が発生するとして、エンドガスの自着火発生時期と自着火による発熱量とを予測していたが、第２の実施形態においては、着火遅れ時間τの経過中における温度・圧力の変化が大きい場合に、予測精度を十分にするためのものである。
【００７１】
そこで本実施形態においては、自着火予測モデルとして、時間経過と共に着火遅れ時間τが変化することを考慮したＬｉｖｅｎｇｏｏｄ−Ｗｕ積分と呼ばれるものを採用する。このモデルは、ある瞬間における温度・圧力等からその状態下での着火遅れ時間τを求め、この逆数を時間推移に沿って積分するものである。この積分値が１となった時点で、自着火が発生すると判定される。この関係を式に示すと以下のようになる。
【００７２】
【数８】
2004332584.tif 000010
【００７３】
（８）式中の着火遅れ時間τは、（７）式のように温度Ｔ・圧力ｐの関数として表現してもよいし、ＥＣＵ１１内に温度Ｔ・圧力ｐを軸とした着火遅れ時間τのマップを用意し、これを参照してもよい。
【００７４】
時間推移に応じ、サイクルシミュレーションによって未燃ガスの温度Ｔ・圧力ｐ等の変化を予測し、その予測値に基づいて着火遅れ時間τを求め、同時に（８）式に示す着火遅れ時間の逆数の積分を求める。この積分値が１となった時点が、エンドガスの自着火発生時期の予測値θｋｎｏｃｋであり、（６）式にこの予測値θｋｎｏｃｋを代入することで、その時点での熱発生割合を求め、１からこれを減算することで、予測値θｋｎｏｃｋにおける未燃ガス割合、すなわち自着火による発熱量（自着火によって燃焼する発熱割合）の予測値Ｘａｕｔｏ−ｉｇｎが得られる。以降は、第１の実施形態と同様に、トレースノック点火時期とＭＢＴ点火時期とを求め、点火時期制御を行う。
【００７５】
ここで、サイクルシミュレーションにおける伝熱モデルの適用例について説明する。第１の実施形態においては、筒内のガスは断熱状態にあるものと仮定して、（１）、（２）式の第４項（−ｄＱｂ，ｗａｌｌ、−ｄＱｕ，ｗａｌｌ）はゼロであるとしたが、燃焼室壁面からの伝熱については、内燃機関を対象とした伝熱モデルが知られており、これを利用して燃焼室壁面への伝熱をシミュレーション内で考慮することで、予測精度をさらに高精度化できる。
【００７６】
燃焼室壁面への熱伝達モデルとしては、Ｗｏｓｃｈｎｉのモデルを適用する。そして、（１）、（２）式における第４項（−ｄＱｗａｌｌ）は、伝熱面積Ａ、熱伝達率ｈｃ、及び燃焼室壁面の温度Ｔｗａｌｌ用いて、
【００７７】
【数９】
2004332584.tif 000011
【００７８】
と表される。熱伝達率ｈｃは、Ｗｏｓｃｈｎｉのモデルによれば、
【００７９】
【数１０】
2004332584.tif 000012
【００８０】
【数１１】
2004332584.tif 000013
【００８１】
と記述される。ここで、Ｃ、ｍはモデル定数、Ｂはボア径、ｐは筒内圧力、及びＴは筒内ガス平均温度である。（１１）式において、Ｓｐの￣は平均ピストン速度、Ｖｄは排気量、ｐｍはモータリング時の筒内圧力を示している。添え字ｒはリファレンスとなる点を示し、例えば吸気バルブ５の閉時期などに設定される。モデル定数Ｃ１，Ｃ２は圧縮〜燃焼〜膨張というサイクル中のフェーズに応じて変更される。
【００８２】
Ｗｏｓｃｈｎｉのモデルは、燃焼室全体の平均熱伝達を与えるモデルであるので、サイクルシミュレーションのような２領域モデルに適用するためには、何らかの方法で熱伝達量を既燃ガス・未燃ガスに分担する必要がある。
【００８３】
しかしながら、熱発生期間中のエンドガスは、ほぼ断熱状態にあるものと見なせるため、（２）式において−ｄＱｕ，ｗａｌｌ＝０とし、（１）式中の−ｄＱｂ，ｗａｌｌが（９）式に示す平均温度から求めた熱伝達量と等しいとして与えればよい。
【００８４】
以上の説明のように、本実施形態によれば、サイクルシミュレーション内で、未燃ガスの着火遅れ時間τの逆数を積分して自着火を予測する手法により、例えば燃焼速度が大きく、未燃ガスの温度・圧力等の変化が大きい機関においても、エンドガスの自着火発生時期と自着火による発熱量を精度良く予測し、両者の関係からノッキング強度を予測して、これがトレースノック相当以下となるように点火時期を制御することで、過剰な点火時期遅角制御による本来不必要な熱効率低下を生じることなく、機関の耐久性を維持して運転することが可能となる。また、燃焼室壁面への伝熱モデルを採用することで、さらに高精度化が可能となる。
【００８５】
次に、第３の実施形態に係る火花点火式内燃機関の燃焼制御装置ついて説明する。
本実施形態においては、図２中に示される自着火モデルについて、第１の実施形態とは異なるモデルを採用するものである。また、伝熱モデルについては、第２の実施形態と同様に扱うものとし、ここでは説明を省略する。
【００８６】
本実施形態においては、化学素反応モデルに基づく自着火予測手段を採用する。ここで化学素反応モデルとは、通常数十〜数百の素反応の組み合わせによって、全体の化学反応を記述したもので、分子同士の衝突・反応という、より実際の物理現象に基づいたモデルをいう。このモデルにより、未燃ガスの自着火現象に対する、温度・圧力の影響だけでなく、ガス組成や燃料性状の違いについても、高精度な予測を行うことが可能である。そしてこの化学素反応モデルを、サイクルシミュレーションの未燃ガスに対して適用することで、バルブタイミングや圧縮比の変化だけでなく、当量比や残留ガス率、ガソリンオクタン価の変化に対しても、高精度にノッキング現象の予測が可能となる。
【００８７】
以下、化学素反応モデルを適用してエンドガスの自着火を予測する手法について説明する。（１）〜（６）式に示す基礎方程式、及び（９）〜（１１）式に示す伝熱モデルに関しては、既に説明した実施形態と同様である。燃焼室４内のガスは、化学種ｉが質量分率ｙｉの割合で混合したガスであると扱われる。素反応モデルでは、これら化学種ｉの質量分率ｙｉが反応の進行に伴って変化することを、以下の式にて記述する。
【００８８】
【数１２】
2004332584.tif 000014
【００８９】
ここで、ω’ｉ，ｕは化学種ｉの生成速度であり、Ｍｉは化学種ｉの分子量を示している。なお、ω’ｉ，ｕはアレニウス型で反応速度を記述した素反応メカニズムを与えることで、ｐ、Ｔｕ、ｙｉ，ｕの関数として表される。この生成速度は、例えば化学素反応計算用として広く知られているＣＨＥＭＫＩＮライブラリを用いることで、容易に得ることができる。
【００９０】
同様に、反応に伴う内部エネルギの変化ｄｕｕは、化学種の生成速度ω’ｉ，ｕを利用して、以下の式で記述される。
【００９１】
【数１３】
2004332584.tif 000015
【００９２】
また、ガス定数Ｒｕの変化は、Ｒｕｎｉの〜を一般ガス定数、Ｍｕを平均分子量として、
【００９３】
【数１４】
2004332584.tif 000016
【００９４】
【数１５】
2004332584.tif 000017
【００９５】
と表すことができる。
これら素反応モデルに関する式は未燃部（添え字ｕ）についてのみ記述したが、添え字をｂに代えて既燃部にも同様に素反応モデルを適用可能である。これにより未燃部の自着火であるノッキング現象だけでなく、既燃部における熱解離の影響やＮＯの生成などを予測することも可能となる。
【００９６】
（１）〜（６）式及び（９）〜（１１）式に加え、（１２）〜（１５）式の方程式を連立し、クランク角θに応じて解いていく。素反応モデルによる自着火計算は、急激な変化を伴う複雑な系であり、計算には高精度かつ安定なソルバーを使用する。サイクルシミュレーションに素反応モデルを適用したことで、未燃ガスの温度Ｔｕ・圧力ｐｕ等及び化学種ｉの質量分率Ｍｉの変化等が、時間推移に伴って予測される。
【００９７】
素反応モデルにおいては、着火現象も連続した化学反応の一部であるので、着火が発生したか否かについて、何らかの判断基準を設ける必要がある。この基準は、エンドガス部の温度Ｔｕが所定値以上になったこと、とすることが可能である。
【００９８】
上記により、エンドガス部で自着火が発生したことが判断されれば、エンドガスの自着火発生時期の予測値θｋｎｏｃｋが定まり、（６）式にθｋｎｏｃｋを代入することで、その時点での熱発生割合を求め、１からこれを減算することで、θｋｎｏｃｋにおける未燃ガス割合、すなわち自着火による発熱量（自着火によって燃焼する発熱割合）の予測値Ｘａｕｔｏ−ｉｇｎが得られる。以降は、第１の実施形態と同様に、トレースノック点火時期とＭＢＴ点火時期とを求め、点火時期制御を行う。
【００９９】
図９は、本実施形態にて説明した、サイクルシミュレーションによる予測例を示す図であり、横軸は圧縮比、縦軸は点火遅角ΔＡＤＶを示している。そして、機関の圧縮比を変更した場合に、トレースノックとなる点火時期がどの程度変化するかについては、予測値を実線、実験値を破線で示して比較している。なお、◆は１２００ｒｐｍ、◇は２４００ｒｐｍの場合を示している。
【０１００】
異なる機関回転速度において、予測値は定性的・定量的に実験値を良く再現しており、本実施形態による予測の妥当性が示されている。
以上の説明のように、本実施形態によれば、サイクルシミュレーション内で、未燃ガス部に化学素反応モデルを適用して自着火を予測することにより、温度・圧力の変化だけでなく、ガス組成や燃料性状の違いについても、高精度な予測を行うことが可能となり、バルブタイミングや圧縮比の変化だけでなく、当量比や残留ガス率、ガソリンオクタン価の変化に対しても、エンドガスの自着火発生時期と自着火による発熱量とを精度良く予測し、両者の関係からノッキング強度を予測して、これがトレースノック相当以下となるように点火時期を制御することで、過剰な点火時期遅角制御による本来不必要な熱効率低下を生じることなく、機関の耐久性を維持して運転することが可能となる。
【図面の簡単な説明】
【図１】火花点火式内燃機関の燃焼制御装置の構成図
【図２】第１の実施形態における制御フローの概略図
【図３】燃焼室２領域サイクルシミュレーションの概略図
【図４】自着火発生時期と自着火による発熱量との関係についての引用図
【図５】サイクル毎の自着火発生時期と自着火による発熱量とを示す図
【図６】自着火発生時期と自着火による発熱量との関係についてノッキング強度を示す図
【図７】トレースノックにおける自着火発生時期と自着火による発熱量とを示す図
【図８】トレースノック相当となる点火時期を求める概略図
【図９】第３の実施形態における予測手法によるトレースノック感度の予測例を示す図
【符号の説明】
１  シリンダヘッド
２  シリンダブロック
３  ピストン
４  燃焼室
５  吸気バルブ
６  排気バルブ
７  吸気ポート
８  排気ポート
９  燃料噴射弁
１０  点火プラグ
１１  エンジンコントロールユニット
１２  アクセル開度センサ
１３  水温センサ
１４  クランク角センサ
【図１】
2004332584.tif 000018
【図２】
2004332584.tif 000019
【図３】
2004332584.tif 000020
【図４】
2004332584.tif 000021
【図５】
2004332584.tif 000022
【図６】
2004332584.tif 000023
【図７】
2004332584.tif 000024
【図８】
2004332584.tif 000025
【図９】
2004332584.tif 000026