(57)【要約】【課題】少ない定数による燃焼モデルにより精度良く筒内圧の挙動を推定することができる内燃機関の筒内圧推定装置を提供する。【解決手段】内燃機関の運転状態を検出する運転状態検出手段２０と、前記運転状態検出手段２０により検出された運転状態に基づいて、点火後に燃焼による熱発生が開始するまでの期間である着火遅れを算出する着火遅れ算出手段３０２と、前記運転状態に基づいて、燃焼速度を算出する燃焼速度算出手段３０３と、を備え、点火後の着火遅れ期間の経過後に、点火プラグのギャップ部で発生した火炎面が、前記点火プラグのギャップ部を中心とした楕円体状に前記燃焼速度にて筒内壁面まで拡大することを模擬し、前記火炎面の内側を既燃部、外側を未燃部として前記既燃部と前記未燃部の筒内圧をそれぞれ算出し、前記既燃部と前記未燃部の筒内圧に基づいて筒内全体の筒内圧を推定する。【選択図】図６
2015096700.tif 000002
【特許請求の範囲】
【請求項１】
内燃機関の運転状態を検出する運転状態検出手段と、
前記運転状態検出手段により検出された運転状態に基づいて、点火後に燃焼による熱発生が開始するまでの期間である着火遅れを算出する着火遅れ算出手段と、
前記運転状態に基づいて、燃焼速度を算出する燃焼速度算出手段と、を備え、
点火後の着火遅れ期間の経過後に、点火プラグのギャップ部で発生した火炎が、前記点火プラグのギャップ部を中心とした楕円体状に前記燃焼速度にて筒内壁面まで拡大することを模擬し、
前記火炎の内側を既燃部、外側を未燃部として前記既燃部と前記未燃部の筒内圧をそれぞれ算出し、
前記既燃部と前記未燃部の筒内圧に基づいて筒内全体の筒内圧を推定することを特徴とする内燃機関の筒内圧推定装置。
【請求項２】
前記火炎の体積の変化量に基づいて熱発生率を算出することを特徴とする請求項１に記載の内燃機関の筒内圧推定装置。
【請求項３】
前記火炎の体積の変化量に基づいて熱発生率を算出し、
前記熱発生率に基づいて算出された熱量に基づいて前記既燃部の筒内温度を算出し、
前記既燃部の筒内温度に基づいて前記既燃部の筒内圧を算出し、
前記未燃部の筒内圧と前記既燃部の筒内圧と燃焼割合に基づいて筒内全体の筒内圧を推定することを特徴とする請求項１に記載の内燃機関の筒内圧推定装置。
【請求項４】
前記燃焼速度は、前記熱発生率がピークとなった後は漸減することを特徴とする請求項２又は３に記載の内燃機関の筒内圧推定装置。
【請求項５】
前記火炎の体積と筒内容積に基づいて燃焼割合を算出することを特徴とする請求項１から４の何れか一項に記載の内燃機関の筒内圧推定装置。
【請求項６】
自在に変更できる仮燃焼速度を用いて火炎が楕円体状に拡大することを模擬して仮燃焼割合を算出する仮燃焼割合算出手段と、
予め計測しておいた筒内圧から算出された実燃焼割合と前記仮燃焼割合がほぼ一致するまで仮燃焼速度を変更し、前記実燃焼割合と前記仮燃焼割合がほぼ一致した時の前記仮燃焼速度を実燃焼速度として導出する実燃焼速度導出手段と、
点火後の前記実燃焼速度がゼロである期間を実着火遅れとして算出する実着火遅れ期間算出手段と、
前記実燃焼速度及び前記実着火遅れを前記運転状態に基づいて近似できる燃焼速度近似手段及び着火遅れ近似手段を燃焼速度と着火遅れの適合方法として有することを特徴とする請求項１乃至５の何れか一項に記載の内燃機関の筒内圧推定装置。
【請求項７】
前記適合方法で生じる実燃焼速度の初期サージ部について、燃焼速度が一定の場合に燃焼量が同じとなるように着火遅れを補正することを特徴とする請求項６に記載の内燃機関の筒内圧推定装置。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、内燃機関の筒内圧推定装置に係り、詳しくは、リアルタイムに現在の筒内圧をシミュレーションできる内燃機関の筒内圧推定装置に関するものである。
【背景技術】
【０００２】
近年の環境問題、エネルギー問題を背景に、内燃機関の制御装置に対して更なる燃費の向上が要求されている。この要求に対し、従来から点火時期はＭＢＴ（Minimum Advance for Best Torque）に近づくように最適化され、ＶＶＴ（Variable Valve Timing）機構によりポンピングロスや内部ＥＧＲ量を最適化するとともに、このＶＶＴ機構を用いて吸気バルブの遅閉じによるアトキンソンサイクル化も行われており、更にＥＧＲ（Exhaust Gas Recirculation）機構によってもポンピングロスや外部ＥＧＲ量の最適化が行われている。
【０００３】
このように多くの制御パラメータを最適に制御するために、筒内圧を計測して最適な運転状態を調査し、燃費の向上が図られて来たが、最適化を行うためには本来必要となるデータ量も膨大であるため、ＤｏＥ（Design of Experiments）等を用いて計測ポイントを削減するとともに、制御に必要なマップ数を削減するように、ＭＢＣ（Model-Based Calibration）による統計モデリングや数値最適化等も使用されている。
【０００４】
前記方法は統計モデルを用いた制御パラメータの最適化であるが、統計モデルではなく物理モデルを用いた点火時期等の最適化も考えられている。つまり、ＶＶＴやＥＧＲの制御量及び各種センサ出力値から筒内の吸気量やＥＧＲ量を物理モデルにより算出し、更に算出した吸気量、ＥＧＲ量、点火時期に基づいて物理モデルに基づいた燃焼シミュレーションを実施して筒内圧を推定することで、筒内圧を計測していない運転条件下においても、燃費が最良となるＶＶＴ、ＥＧＲ、点火時期の制御量を算出して各制御パラメータの最適化を行うものである。
【０００５】
更にはこの筒内圧の推定や各種制御量の最適化を、内燃機関の制御装置に用いるマイコンの性能の向上に伴い、内燃機関の制御を実行しながらリアルタイムに燃焼シミュレーションを行って、実施することも考えられている。
【０００６】
物理モデルの中でも特に燃焼シミュレーションに用いられる燃焼モデルの例として、例えば特開２００４−２９３５４１号公報（特許文献１）に開示された方法は、乱流燃焼による火炎の成長と層流燃焼による火炎の成長を用いて火炎伝播を予測することで燃焼の物理モデルを構築している。また、例えば特開２００４−３３２６５８号公報（特許文献２）に開示された方法は、層流燃焼速度に基づいて主燃焼期間を算出し、更に主燃焼期間に基づいてＭＢＴ点火時期の算出を行っている。このような、燃焼の物理モデルに忠実な燃焼シミュレーションとは異なるアプローチとして、例えば特開２００８−２１５２０４号公報（特許文献３）に記載の方法は、燃料割合を近似するＷｉｅｂｅ関数を用いて熱発生率のシミュレーションを実施している。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開２００４−２９３５４１号公報
【特許文献２】特開２００４−３３２６５８号公報
【特許文献３】特開２００８−２１５２０４号公報
【発明の概要】
【発明が解決しようとする課題】
【０００８】
しかしながら、前記特許文献１及び２に開示された方法は、現在知られている燃焼の物理モデルにある程度忠実な燃焼シミュレーションを行っているものの、実際の筒内での燃焼現象は非常に複雑であり、内燃機関の制御装置に実装してリアルタイムに燃焼をシミュレーションするレベルの演算量では高精度な推定は困難であると考えられる。
【０００９】
また、特許文献３に開示された方法では、燃焼の進行状況を近似したものであるので、演算量は少ないものの筒内の形状や運転条件に応じてＷｉｅｂｅ関数の定数も変更する必要があり、筒内形状に起因して熱発生率の形状が複雑（例えば、熱発生率のピークを中心に左右が非対称である場合等）であれば、複数のＷｉｅｂｅ関数を用いてフィッティングする必要がある等、高精度にシミュレーションするには多くの定数とその適合が必要になってしまう課題がある。
【００１０】
この発明は、前記のような問題点を解決するためになされたもので、少ない定数による燃焼モデルにより精度良く筒内圧の挙動を推定することができる内燃機関の筒内圧推定装置を提供することを目的としている。
【課題を解決するための手段】
【００１１】
この発明による内燃機関の筒内圧推定装置は、内燃機関の運転状態を検出する運転状態検出手段と、前記運転状態検出手段により検出された運転状態に基づいて、点火後に燃焼による熱発生が開始するまでの期間である着火遅れを算出する着火遅れ算出手段と、前記運転状態に基づいて、燃焼速度を算出する燃焼速度算出手段と、を備え、
点火後の着火遅れ期間の経過後に、点火プラグのギャップ部で発生した火炎面が、前記点火プラグのギャップ部を中心とした楕円体状に前記燃焼速度にて筒内壁面まで拡大することを模擬し、前記火炎面の内側を既燃部、外側を未燃部として前記既燃部と前記未燃部の筒内圧をそれぞれ算出し、前記既燃部と前記未燃部の筒内圧に基づいて筒内全体の筒内圧を推定するものである。
【発明の効果】
【００１２】
この発明による内燃機関の筒内圧推定装置によれば、少ない制御定数により精度良く筒内圧の挙動を推定することができる。
【図面の簡単な説明】
【００１３】
【図１】この発明の実施の形態１による筒内圧推定装置を搭載した内燃機関を概略的に示す構成図である。
【図２】この発明の実施の形態１による筒内圧推定装置を搭載した内燃機関の制御装置を概略的に示す構成図である。
【図３】燃焼パラメータ近似手段を示すブロック図である。
【図４】初期筒内圧を示す筒内圧−筒内容積線図である。
【図５】着火遅れと燃焼速度を示す図である。
【図６】燃焼シミュレーションの演算内容を示すフローチャートである。
【図７】筒内における火炎の成長を模式的に示す図である。
【図８】火炎の表面積を演算する方法を模式的に示す図である。
【図９】火炎の体積を演算する方法を模式的に示す図である。
【図１０】既燃部の筒内温度を算出する方法を模式的に示す図である。
【図１１】筒内圧から燃焼速度を算出する方法を示すブロック図である。
【発明を実施するための形態】
【００１４】
以下、この発明による内燃機関の筒内圧推定装置の好適な実施の形態について図面を用いて説明する。
【００１５】
実施の形態１．
図１は、この発明の実施の形態１による筒内圧推定装置を搭載した内燃機関（以下、エンジンとも言う。）を概略的に示す構成図であり、図２は、その制御装置を概略的に示す構成図である。
【００１６】
図１において、エンジン１の吸気系の上流に吸入空気量を検出するエアフロセンサ（以下、ＡＦＳ：Air Flow Sensorと言う。）２が設けられている。ＡＦＳ２の下流のエンジン１側には、吸入空気量を調整するために電気的に制御することができる電子制御スロットル３が設けられている。また、電子制御スロットル３の開度を測定するために、スロットル開度センサ４が設けられている。なお、吸入空気量の検出にはＡＦＳ２を用いる代わりに、後述するインマニ圧に基づいて吸入空気量を推定するなど、別の吸入空気量を測定する手段を用いるようにしても良い。
【００１７】
更に、電子制御スロットル３の下流に設けられたサージタンク５及びインテークマニホールド６内を含む空間（以下、インマニと言う。）の圧力（以下、インマニ圧と言う。）を測定するインマニ圧センサ７と、インマニ内の温度（以下、インマニ温と言う。）を測定する吸気温センサ８が設けられている。なお、インマニ温を計測する吸気温センサ８を設ける代わりに、厳密には異なる温度となるが、近似的に外気を計測する温度センサ、例えば、ＡＦＳ２に内蔵されている温度センサを用いて、外気温よりインマニ温を推定することもできる。
【００１８】
更に、サージタンク５には、排気系からの外部ＥＧＲを導入するためのＥＧＲバルブ９からの経路が接続されている。インテークマニホールド６及び筒内を含む吸気バルブの近傍には燃料を噴射するためのインジェクタ１０が設けられ、吸気バルブ及び排気バルブには、バルブタイミングを可変して内部ＥＧＲを操作するための吸気ＶＶＴ１１と排気ＶＶＴ１２がそれぞれ設けられており、シリンダヘッドにはシリンダ内で火花を発生させる点火プラグを駆動するための点火コイル１３が設けられている。エキゾーストマニホールド１４には図示しないＯ２センサや触媒、ＥＧＲバルブ９への経路が設けられている。更に、エンジン１のクランク軸に取り付けられた回転プレート１５のエッジからエンジン回転速度やクランク角度を検出するクランク角度センサ１６が設けられている。
【００１９】
また、電子制御ユニット（以下、ＥＣＵ：Electric Control Unit）１７が設けられ、ＥＣＵ１７には、図２に示すように、ＡＦＳ２で測定された吸入空気量と、スロットル開度センサ４で測定された電子制御スロットル３の開度と、インマニ圧センサ７で測定されたインマニ圧と、吸気温センサ８で測定されたインマニ温と、クランク角度センサ１６で検出された回転プレート１５のエッジと、大気圧センサ１８で測定された大気圧が入力される。なお、大気圧を測定する大気圧センサ１８の代わりに、大気圧を推定する手段を用いても良いし、ＥＣＵ１７に内蔵された大気圧センサを用いても良い。また前記以外の各種センサ１９（図示しないアクセル開度センサ、Ｏ２センサを含む）からもＥＣＵ１７に測定値が入力される。
【００２０】
ＥＣＵ１７内の運転状態検出手段２０では、回転プレート１５のエッジ周期からエンジン回転速度やクランク角度を算出し、更に吸入空気量やインマニ圧等からシリンダ吸入空気量やＥＧＲ率が算出される。ここで算出されたシリンダ吸入空気量やＥＧＲ率に基づいて、インジェクタ１０、点火コイル１３等が駆動される。また、アクセル開度等の入力された各種データより目標トルクが算出され、算出された目標トルクを達成する目標シリンダ吸入空気量が算出され、目標シリンダ吸入空気量及び目標ＥＧＲ率を達成するように目標スロットル開度、目標吸気ＶＶＴ位相角、目標排気ＶＶＴ位相角、目標ＥＧＲ開度を算出し、これらを目標値として電子制御スロットル３の開度、吸気ＶＶＴ１１及び排気ＶＶＴ１２の位相角、ＥＧＲバルブ９の開度が制御される。また、その他の各種アクチュエータ２１も必要に応じて制御される。
【００２１】
ＥＣＵ１７内では以上のような通常のエンジン制御の他に、燃焼シミュレーションも同時に実行される。これは、前述の運転状態検出手段２０で求めた運転状態に基づいて、燃焼パラメータ近似手段２２にて燃焼シミュレーションに必要な燃焼パラメータを算出し、運転条件と燃焼パラメータに基づいて燃焼シミュレーション手段２３にて燃焼シミュレーションを行っている。
【００２２】
次に図３を参照しながら、燃焼パラメータ近似手段２２について説明する。図３は燃焼パラメータ近似手段２２を示すブロック図である。運転状態検出手段２０では、前述の通り、通常のエンジン制御で用いられる運転状態である、エンジン回転速度（Ｎｅ）、充填効率（Ｅｃ）、点火時期（ＩＧ）、空燃比（ＡＦＲ）、ＥＧＲ率（ＩｎＥＧＲ：内部ＥＧＲ率、ＥｘＥＧＲ：外部ＥＧＲ率、ＥＧＲ：総ＥＧＲ率（＝ＩｎＥＧＲ＋ＥｘＥＧＲ））、ＶＶＴ位相角（ＩｎＶＶＴ：吸気側、ＥｘＶＶＴ：排気側）、インマニ圧（Ｐｂ）等が算出される。ところで、燃焼シミュレーションを実施するにはこれらのパラメータだけでは足りず、筒内圧に関するパラメータや燃焼に関するパラメータが必要になる。そこで、燃焼パラメータ近似手段２２では、ブロック３０１にて初期筒内圧（Ｐｉｎ）を、ブロック３０２にて着火遅れ（ＩｇＤｌｙ）を、ブロック３０３にて燃焼速度（Ｖｃｏｍｂ）を算出する。なお、ブロック３０１は初期筒内圧算出手段であり、ブロック３０２は着火遅れ算出手段であり、ブロック３０３は燃焼速度算出手段となる。
【００２３】
まずブロック３０１での初期筒内圧（Ｐｉｎ）の算出方法について説明する。初期筒内圧（Ｐｉｎ）とは、圧縮行程の開始（下死点）時点における筒内圧である。一般にエンジンの圧縮行程においては、吸気バルブ閉から点火時期までの間の空気の出入りのない状態変化中において、筒内容積（Ｖ）と筒内圧（Ｐ）は次式（１）で示されるポリトロープ変化として記述できることが知られている。この式において、ｎはポリトロープ指数を示している。
【００２４】
【数１】
2015096700.tif 000003
【００２５】
この式を用いて図４に示すように、図中の×印の筒内容積と筒内圧から下死点（ＢＴＤＣ１８０ｄｅｇＣＡ）における筒内圧を算出したものを初期筒内圧と定義している。この初期筒内圧は、例えばＤｏＥ等により予め決めた運転条件で計測した筒内圧を用いて、各運転条件での初期筒内圧を算出しておき、算出した初期筒内圧を運転条件（エンジン回転速度、インマニ圧、ＶＶＴ位相角等）を用いて近似しておくことで、図３のブロック３０１において算出することができる。初期筒内圧（Ｐｉｎ）を算出するための近似式の例として、エンジン回転速度（Ｎｅ）、インマニ圧（Ｐｂ）、吸気ＶＶＴ位相角（ＩｎＶＶＴ）、排気ＶＶＴ位相角（ＥｘＶＶＴ）を用いた１次近似式は次式（２）のように書くことができる。この式において、Ｋ＊＿Ｐｉｎは比例係数を示し、Ｋｃｎｔ＿Ｐｉｎは定数項を示している。
【００２６】
【数２】
2015096700.tif 000004
【００２７】
なお、初期筒内圧（Ｐｉｎ）を算出する近似式としては上式に限るものではなく、より多くの運転状態を示すパラメータを用いるようにしても良いし、２次式、３次式のように次数を上げた近似式を用いるようにしても良い。
【００２８】
次に、ブロック３０２に示す着火遅れ（ＩｇＤｌｙ）とブロック３０３に示す燃焼速度（Ｖｃｏｍｂ）の算出方法について説明する。
まず、図５（ａ）、（ｂ）（ｃ）（ｄ）に示す後述の方法により算出された燃焼速度の例に基づいて着火遅れと燃焼速度の定義を説明する。ここで、図５（ａ）はクランク角と燃焼割合（ＭＦＢ）を、図５（ｂ）はクランク角と熱発生率（ｄＱ）を、図５（ｃ）はクランク角と燃焼速度（Ｖｃｏｍｂ(i)）を、図５（ｄ）はクランク角と補正後燃焼速度（Ｖｃｏｍｂ）の関係を示している。着火遅れは、図５（ｄ）に示されるように点火時期から燃焼速度が立ち上がるまでの遅れ期間として定義し、この遅れ期間を得る操作が実着火遅れ期間算出手段となる。なお、前記遅れ期間の単位としては、時間［ｍｓ］を用いても良いし、クランク角度［ｄｅｇＣＡ］を用いても良い。
なお、図５（ｃ）に示すように、燃焼開始時に燃焼速度のサージが生じる場合があり、後述するように、このサージを考慮して着火遅れを補正するようにしても良い。また、燃焼速度は、燃焼開始時の初期サージ部及び燃焼後半の漸減部を除いた部分の平均値を燃焼速度と定義する。これらも初期筒内圧の算出と同様に、ＤｏＥ等により予め決めた運転条件で計測した筒内圧を用いて、着火遅れと燃焼速度を算出しておき、算出した着火遅れと燃焼速度を運転条件（エンジン回転速度、充填効率、点火時期、空燃比、ＥＧＲ率等）を用いて近似することで算出することができる。着火遅れ（ＩｇＤｌｙ）と燃焼速度（Ｖｃｏｍｂ）を算出する近似式の例として、エンジン回転速度（Ｎｅ）、充填効率（Ｅｃ）、点火時期（ＩＧ）、空燃比（ＡＦＲ）、ＥＧＲ率（ＥＧＲ）を用いた１次近似は次式（３）（４）のように書くことができる。この式において、Ｋ＊＿＊は比例係数、Ｋｃｎｔ＿＊は定数項を示している。
【００２９】
【数３】
2015096700.tif 000005
【００３０】
【数４】
2015096700.tif 000006
【００３１】
なお、前記（３）式が着火遅れ算出手段の実現形態である着火遅れ近似手段に、また、前記（４）式が燃焼速度算出手段の実現形態である燃焼速度近似手段にそれぞれ相当するが、着火遅れ（ＩｇＤｌｙ）と燃焼速度（Ｖｃｏｍｂ）を算出する近似式としては上式に限るものではなく、より多くの運転状態を示すパラメータを用いるようにしても良いし、２次式、３次式のように次数を上げた近似式を用いるようにしても良い。また、燃焼速度は、図５（ｃ）に示すように、熱発生率のピーク位置近傍から漸減する傾向があるため、後述する燃焼シミュレーション時において、熱発生率がピークとなったことを検出した後は燃焼速度を漸減補正するようにしても良い。
【００３２】
燃焼シミュレーション手段２３においては、運転状態検出手段２０により算出される運転状態と、燃焼パラメータ近似手段２２により算出される初期筒内圧、着火遅れ、燃焼速度を用いて燃焼シミュレーションが実施される。
【００３３】
次に、燃焼シミュレーション手段２３にて実行される演算について、図６を参照しながら説明する。図６は、燃焼シミュレーションの演算内容を示すフローチャートである。この演算は、圧縮行程の開始（下死点）から膨張行程の終了（下死点）までの３６０［ｄｅｇＣＡ］のクランク角度内に実行される。より具体的には、所定のクランク角度における割り込み処理内で演算しても良いし、複数のクランク角度における割り込み処理内で分散して演算しても良い。また、所定時間毎における割り込み処理内で分散して演算しても良い。なお、理想的には全運転域でクランク角度３６０［ｄｅｇＣＡ］の間に演算が終了できれば良いが、高回転域で演算処理が終了できない場合には、例えばクランク角度７２０［ｄｅｇＣＡ］の間、つまり、エンジン１が２回転する間に１燃焼を演算するようにしても良い。
【００３４】
ステップ６０１ではまず燃焼シミュレーションの初期化処理が実行される。ここでは、筒内圧（Ｐｃｙｌ）の初期値として初期筒内圧（Ｐｉｎ）が代入され、クランク角度（Ａｎｇ）の初期値として圧縮行程の開始点（下死点）に相当する−１８０［ｄｅｇＣＡ］が代入され、次式（５）に基づいて着火時期（ＩＡ）が点火時期（ＩＧ）と着火遅れ（ＩｇＤｌｙ）から算出される。この式において、（０）は初期値を表している。
【００３５】
【数５】
2015096700.tif 000007
【００３６】
更に、充填効率（Ｅｃ）と総ＥＧＲ率（ＥＧＲ）と筒内容積（Ｖ（０））から筒内ガスの密度（ρ（０））を算出し、次式（６）の状態方程式を用いて筒内温度の初期値（Ｔｃｙｌ（０））が算出される。
【００３７】
【数６】
2015096700.tif 000008
【００３８】
ここで、Ｒはガス定数である。また筒内容積（Ｖ）は、予めクランク角度毎の筒内容積をクランク−ピストン機構に基づいて求めておいて配列に記憶しておき、使用時にクランク角度（Ａｎｇ）に応じた筒内容積（Ｖ）を参照するようにしても良いし、使用する度にクランク−ピストン機構に基づいて求めるようにしても良い。
【００３９】
続くステップ６０２では、クランク角度が更新される。クランク角度の更新量は、例えば、１［ｄｅｇＣＡ］とすれば良い。但し、高回転域で演算量を削減する必要がある場合には、エンジン回転速度に応じてクランク角度の更新量を２［ｄｅｇＣＡ］とするようにしても良いし、逆に低回転域や演算処理能力に余裕があるＥＣＵ１７では０．５［ｄｅｇＣＡ］とするようにしても良い。
【００４０】
ステップ６０３では、クランク角度更新前の筒内圧と筒内温度を用いてポリトロープ変化演算を行って、現在のクランク角度における筒内圧と筒内温度を算出する。具体的には、現在のクランク角度を（ｉ）、更新前のクランク角度を（ｉ−１）すると、ポリトロープ変化は次式（７）で演算できる。
【００４１】
【数７】
2015096700.tif 000009
【００４２】
なお、後述する通りこの演算は筒内の既燃部と未燃部のそれぞれにおいて演算する。
【００４３】
ステップ６０４では、現在のクランク角度（Ａｎｇ（ｉ））が着火時期（ＩＡ）より大きいかどうか、つまり着火時期（ＩＡ）を過ぎたかどうかを判定する。ここで、現在のクランク角度（Ａｎｇ（ｉ））が着火時期以前である場合には、以降の処理は実行せずにステップ６０９へ進むが、着火時期以降である場合には、続くステップ６０５へ進む。ステップ６０５では、現在の既燃部体積（Ｖｂ（ｉ））が現在の筒内容積（Ｖ（ｉ））より小さいかどうか、つまり、燃焼が終了したかどうかを判定する。ここで既燃部体積（Ｖｂ（ｉ））が筒内容積（Ｖ（ｉ））と等しければ燃焼は終了しているため、以降の処理は実行せずにステップ６０９へ進むが、燃焼が終了していない場合には続くステップ６０６へ進む。
【００４４】
ステップ６０６では火炎表面積（Ｓｂ）を算出し、続くステップ６０７では既燃部体積（Ｖｂ）を算出する。ここで、実施の形態１における火炎形状の考え方について図７を用いて説明する。図７（ａ）、（ｂ）、（ｃ）は、筒内における火炎の成長を模式的に示す図である。図７（ａ）、（ｂ）、（ｃ）の円柱は筒内の形状を模式的に表しており、内部の球は火炎を模式的に表している。
【００４５】
図７（ａ）は点火直後の火炎を、図７（ｂ）は燃焼割合（ＭＦＢ）が５０％程度進行した場合の火炎を、図７（ｃ）は燃焼が終了する直前の火炎を示すイメージ図であり、点火プラグのギャップ部を中心として火炎が燃焼速度（Ｖｃｏｍｂ）で拡大していく様子を表している。このように火炎の成長を考えることで筒内の燃焼が模擬できる。なお、火炎は球面状である必要はなく、火炎の赤道方向（図７（ａ）、（ｂ）、（ｃ）では左右方向に相当）と極方向（図７（ａ）、（ｂ）、（ｃ）では上下方向に相当）の半径を別の値とした回転楕円体として考えても良い。この場合は、例えば極方向へ火炎が成長する燃焼速度を燃焼速度（Ｖｃｏｍｂ）に所定値を掛けて算出するようにすれば良い。また、火炎が筒内壁面に接触した場合には、その部分の燃焼は終了したものとみなし、その部分は火炎表面積や既燃部体積の演算からは除外するようにする。
【００４６】
以上の考え方を踏まえ、図８を参照しながら、ステップ６０６における火炎表面積（Ｓｂ）の算出方法について説明する。図８は火炎の表面積を演算する方法を模式的に示す図である。点火プラグのギャップ部を中心とした回転楕円体を考え、図８のように回転楕円体の赤道方向をｘ、ｙ軸としてその半径をａ、極方向をｚ軸としてその半径をｂとすると、ｚ軸に垂直な面で輪切りにした断面の半径Ｒは、次式（８）となるので、これを用いて回転楕円体の表面の微小部の長さｄＬを算出すると、次式（９）となる。
【００４７】
【数８】
2015096700.tif 000010
【００４８】
【数９】
2015096700.tif 000011
【００４９】
更にこれをｚ軸中心に１回転させた部分の面積ｄＳは、次式（１０）となる。
【００５０】
【数１０】
2015096700.tif 000012
【００５１】
ＥＣＵ１７において現在のクランク角度における火炎表面積（Ｓｂ（ｉ））を算出するには、（１０）式を火炎のｚ軸方向の高さまで数値積分すれば良い。なお、前述の通り火炎が筒内の壁面に接触した部分は燃焼が終了しているので、その部分は数値積分から除外する必要がある。
【００５２】
続いて、図９を参照しながら、ステップ６０７における既燃部体積（Ｖｂ）の算出方法について説明する。図９は既燃部の体積を演算する方法を模式的に示す図である。表面積の算出と同様に、ある軸に垂直な面で輪切りにした断面について、その断面積Ｓｖ（ｚ）は、次式（１１）のように書くことができる。
【００５３】
【数１１】
2015096700.tif 000013
【００５４】
ＥＣＵ１７において現在のクランク角度における既燃部体積（Ｖｂ（ｉ））を算出するには、（１１）式に微小高さｄｚを掛けたものを既燃部のｚ軸方向の高さまで数値積分すれば良い。なお、既燃部は筒内の壁面までであるので、これを超える部分は数値積分から除外する必要がある。
【００５５】
以上のようにして、ステップ６０６及びステップ６０７において火炎表面積（Ｓｂ）と既燃部体積（Ｖｂ）が算出される。なお、ここでは火炎の形状を回転楕円体と考えた場合の例を示したが、火炎の形状を楕円体と考えても良いし、球と考えても良い。楕円体とは、ｘ軸、ｙ軸、ｚ軸方向の３つの径の長さがすべて異なるものであり、そのいずれか２つの径の長さの等しいものが回転楕円体、すべての径の長さの等しいものが球となる。ここでは、これらを総称して楕円体状と言う。
【００５６】
続くステップ６０８では筒内圧及び筒内温度の算出を行う。まず、図１０を参照しながら筒内温度の算出方法について説明する。図１０は、既燃部、未燃部、新たに燃焼した部分の筒内温度を算出する方法を模式的に示す図である。図１０に示すように、実施の形態１による方法では、既燃部と未燃部と新たに燃焼した部分を分けて考えており、未燃部と既燃部については、（７）式に示すポリトロープ変化によりクランク角度更新前の未燃部と既燃部の筒内温度（Ｔｃ１（ｉ−１）、Ｔｃ２（ｉ−１））から、次式（１２）に示すように、現在のクランク角度における未燃部と既燃部の筒内温度（Ｔｃ１（ｉ）、Ｔｃ２（ｉ））を算出する。また、新たに燃焼した部分に関しては、未燃部の温度（Ｔｃ１（ｉ））に加え、燃焼により発生した熱量で上昇する温度（ΔＴｃ（ｉ））を考慮して新たに燃焼した部分の筒内温度を算出する。そして、既燃部の筒内温度（Ｔｃ２（ｉ））と新たに燃焼した部分の温度（Ｔｃ１（ｉ）＋ΔＴｃ（ｉ））と燃焼割合（ＭＦＢ）に基づいて、現在のクランク角度における既燃部の筒内温度の最終値（Ｔｃ２’（ｉ））を算出している。
【００５７】
【数１２】
2015096700.tif 000014
【００５８】
ここで、現在のクランク角度における燃焼割合（ＭＦＢ（ｉ））は、既燃部体積（Ｖｂ（ｉ））と筒内容積（Ｖ（ｉ））を用いて次式（１３）で算出できる。
【００５９】
【数１３】
2015096700.tif 000015
【００６０】
また、燃焼により発生した熱量（ΔＱ（ｉ））で上昇する温度（ΔＴｃ（ｉ））は、定容比熱（Ｃｖ）を用いて次式（１４）で算出できる。なお、燃焼により発生した熱量（ΔＱ（ｉ）：熱発生率に相当）は更新前と現在のクランク角度における既燃部体積の変化量（ΔＶｂ（ｉ）＝Ｖｂ（ｉ）―Ｖｂ（ｉ−１））と、筒内燃料密度（ρｆ（ｉ）＝ρ（ｉ）／（１＋ＡＦＲ））と燃料の発熱量（ＨＶＦ）より算出することができる。
【００６１】
【数１４】
2015096700.tif 000016
【００６２】
以上により算出した未燃部と既燃部の筒内温度（Ｔｃ１（ｉ）、Ｔｃ２’（ｉ））と燃焼割合（ＭＦＢ（ｉ））から、筒内全体の筒内温度（Ｔｃｙｌ（ｉ））は次式（１５）により算出される。
【００６３】
【数１５】
2015096700.tif 000017
【００６４】
なお、ここで算出した熱発生率（ΔＱ（ｉ））のピークホールド値（ΔＱｐｈ）を算出しておき、このピークホールド値（ΔＱｐｈ）と熱発生率（ΔＱ（ｉ））の差が所定値以上になった場合に、熱発生率のピーク位置と検出して燃焼速度を漸減補正するようにすれば、図５（ｃ）において説明した燃焼速度の挙動をより正確に模擬することができる。
【００６５】
次に筒内圧の算出方法について説明する。未燃部については、（７）式に示すポリトロープ変化によりクランク角度更新前の未燃部の筒内圧（Ｐｃ１（ｉ−１））から、現在のクランク角度における未燃部の筒内圧（Ｐｃ１（ｉ））を次式（１６）により算出する。
【００６６】
【数１６】
2015096700.tif 000018
【００６７】
既燃部の筒内圧の最終値（Ｐｃ１’（ｉ））は、現在のクランク角度における既燃部の筒内温度の最終値（Ｔｃ２’（ｉ））に基づいて、次式（１７）に示す状態方程式を用いて算出する。
【００６８】
【数１７】
2015096700.tif 000019
【００６９】
以上により算出した未燃部と既燃部の筒内圧（Ｐｃ１（ｉ）、Ｐｃ２’（ｉ））と燃焼割合（ＭＦＢ（ｉ））から、筒内全体の筒内圧（Ｐｃｙｌ（ｉ））は次式（１８）により算出される。
【００７０】
【数１８】
2015096700.tif 000020
【００７１】
このようにして、ステップ６０８では筒内圧及び筒内温度の算出を行う。続くステップ６０９では、現在のクランク角度（Ａｎｇ（ｉ））が膨張行程の終了点（下死点）である終了角度（＋１８０［ｄｅｇＣＡ］）と比較され、終了していれば処理は終了し、終了角度以前であれば、ステップ６０２に戻ってクランク角度の更新が行われる。
【００７２】
以上のようにして燃焼シミュレーション手段２３での燃焼シミュレーションを実行することで、筒内の火炎が楕円体状に拡大することを模擬しているので、筒内形状に応じた筒内圧変化の挙動を容易に再現することが可能となり、熱発生率や燃焼割合も容易に算出することができ、点火プラグの位置に応じた筒内圧変化の挙動も容易に再現することが可能となる。
【００７３】
ところで、これまでの説明は火炎の燃焼速度が既知であることが前提である。以下ではこれまで説明した燃焼シミュレーションを用いて、ＤｏＥ等により予め決めた運転条件で計測した筒内圧から燃焼速度を算出する方法について説明する。
図１１は、筒内圧から燃焼速度を算出する方法を示すブロック図である。ブロック１１０１乃至１１０３では、筒内圧データから、熱発生率（ｄＱ）と燃焼割合（ＭＦＢ）を算出する。筒内圧データから現在のクランク角度の熱発生率（ｄＱ（ｉ））を算出するには、一般に次式（１９）が用いられる。
【００７４】
【数１９】
2015096700.tif 000021
【００７５】
ここで、体積変化量（ｄＶ（ｉ））と筒内圧変化量（ｄＰｃｙｌ（ｉ））は、前後のクランク角データから差分により算出することができる。
【００７６】
現在のクランク角度の燃焼割合（ＭＦＢ（ｉ））は、次式（２０）に示す現在のクランク角度までの熱発生率の積算値（Ｑ（ｉ））と全熱発生率の積算値（Ｑａｌｌ）の比として算出できる。
【００７７】
【数２０】
2015096700.tif 000022
【００７８】
このようにして、予め圧縮行程の開始（下死点）から膨張行程の終了（下死点）までの３６０［ｄｅｇＣＡ］のクランク角度の筒内圧データから、熱発生率（ｄＱ）と燃焼割合（ＭＦＢ）を算出しておく。
【００７９】
次に、図６に示す燃焼シミュレーションを用いて燃焼速度を算出するが、より具体的な燃焼速度の算出方法を図１１のブロック１１０４乃至１１０８に示している。
図１１のブロック１１０４乃至１１０８は、圧縮行程の開始（下死点）から膨張行程の終了（下死点）までの３６０［ｄｅｇＣＡ］のクランク角度毎に実行される。ブロック１１０４では、仮燃焼速度の設定・更新を行う。仮燃焼速度とは、真の燃焼速度が未知の状態で仮値として燃焼速度を設定し、その仮燃焼速度を用いて図６の燃焼シミュレーションを実施するためのものである。ここではまず仮燃焼速度を０として初期設定を行う。続くブロック１１０５では仮燃焼速度を用いて仮既燃部体積の算出を、ブロック１１０６では仮燃焼割合（ＭＦＢ’（ｉ））の算出を行うが、その算出にはこれまで説明してきた図６に示す燃焼シミュレーションの特にステップ６０７の方法が用いられる。そして、ブロック１１０７では、そのクランク角度における筒内圧データから算出した燃焼割合（ＭＦＢ（ｉ））と仮燃焼割合（ＭＦＢ’（ｉ））が比較される。なお、ブロック１１０６は仮燃焼割合算出手段であり、ブロック１１０８は実燃焼速度導出手段である。
【００８０】
ここで燃焼が始まる前（着火時期以前）であれば、筒内圧データから算出した燃焼割合（ＭＦＢ（ｉ））も仮燃焼割合（ＭＦＢ’（ｉ））も０であり、不等式が成立するためブロック１１０８に進んで仮燃焼速度を、現在のクランク角度における燃焼速度（Ｖｃｏｍｂ（ｉ））として格納する。燃焼が始まった後であれば、仮燃焼速度が初期設定のままでは筒内圧データから算出した燃焼割合（ＭＦＢ）の方が仮燃焼割合（ＭＦＢ’（ｉ））よりも大きくなるのでブロック１１０４に戻り、仮燃焼速度に所定値を加算して仮燃焼速度を更新し、ブロック１１０５及びブロック１１０６の演算を行う。これをブロック１１０７の条件が成立するまで繰り返すことで筒内圧データから算出した燃焼割合（ＭＦＢ（ｉ））と仮燃焼割合（ＭＦＢ’（ｉ））が等しくなる燃焼速度を求めることができるようになる。当然のことながら、このクランク角度毎の燃焼速度（Ｖｃｏｍｂ（ｉ））を用いて燃焼シミュレーションを実施すれば、計測した筒内圧データにほぼ等しいシミュレーション結果を得ることができる。
【００８１】
ところで、このような方法で燃焼速度を算出すると、図５（ｃ）に示すような燃焼速度の立ち上がり時に初期サージのようなものが発生する場合がある。そこで、この初期サージ部及び燃焼後半の漸減部を除いた部分の平均値を燃焼速度と定義し、クランク角と燃焼速度で囲まれる部分の面積が等しくなるように着火遅れも補正することで、計測した筒内圧データに対する着火遅れと燃焼速度を算出することができる。
【００８２】
以上のようにすることで、実施の形態１の燃焼シミュレーションに適した着火遅れと燃焼速度を算出できるようになり、筒内圧データを用いた適合を容易に行うことが可能となり、着火遅れをより精度良く適合することが可能となる。
【００８３】
以上詳説したように、実施の形態１による内燃機関の筒内圧推定装置は、内燃機関の運転状態を検出する運転状態検出手段と、前記運転状態検出手段により検出された運転状態に基づいて、点火後に燃焼による熱発生が開始するまでの期間である着火遅れを算出する着火遅れ算出手段と、前記運転状態に基づいて、燃焼速度を算出する燃焼速度算出手段と、を備え、
点火後の着火遅れ期間の経過後に、点火プラグのギャップ部で発生した火炎が、前記点火プラグのギャップ部を中心とした楕円体状に前記燃焼速度にて筒内壁面まで拡大することを模擬し、前記火炎の内側を既燃部、外側を未燃部として前記既燃部と前記未燃部の筒内圧をそれぞれ算出し、前記既燃部と前記未燃部の筒内圧に基づいて筒内全体の筒内圧を推定することを特徴とする。
この構成によれば、筒内の火炎が楕円体状に拡大することを模擬しているので、筒内形状に応じた筒内圧変化の挙動を容易に再現することが可能となる。
【００８４】
また、実施の形態１による内燃機関の筒内圧推定装置は、火炎の体積の変化量に基づいて熱発生率を算出することを特徴とする。
この構成によれば、熱発生率のみを算出したい場合には、筒内圧を算出することなく熱発生率を算出することが可能となる。
【００８５】
更に、実施の形態１による内燃機関の筒内圧推定装置は、火炎の体積の変化量に基づいて熱発生率を算出し、前記熱発生率に基づいて算出された熱量に基づいて既燃部の筒内温度を算出し、前記既燃部の筒内温度に基づいて前記既燃部の筒内圧を算出し、未燃部の筒内圧と前記既燃部の筒内圧と燃焼割合に基づいて筒内全体の筒内圧を推定することを特徴とする。
この構成によれば、算出された熱発生率と燃焼割合に基づいて筒内圧を推定することが可能となる。
【００８６】
また、実施の形態１による内燃機関の筒内圧推定装置は、燃焼速度が熱発生率のピーク時点後は漸減することを特徴とする。
この構成によれば、熱発生率がピークに達したのちに未燃ガスの燃焼が終了するまでの筒内圧変化の挙動を容易に再現することが可能となる。
【００８７】
また、実施の形態１による内燃機関の筒内圧推定装置は、火炎形状の体積と筒内容積に基づいて燃焼割合を算出することを特徴とする。
この構成によれば、燃焼割合のみを算出したい場合には、筒内圧や熱発生率を算出することなく燃焼割合を算出することが可能となる。
【００８８】
また、実施の形態１による内燃機関の筒内圧推定装置は、自在に変更できる仮燃焼速度を用いて火炎が楕円体状に拡大することを模擬して仮燃焼割合を算出する仮燃焼割合算出手段と、予め計測しておいた筒内圧から算出された実燃焼割合と前記仮燃焼割合がほぼ一致するまで仮燃焼速度を変更し、前記実燃焼割合と前記仮燃焼割合がほぼ一致した時の前記仮燃焼速度を実燃焼速度として導出する実燃焼速度導出手段と、点火後の前記実燃焼速度がゼロである期間を実着火遅れとして算出する実着火遅れ期間算出手段と、前記実燃焼速度及び前記実着火遅れを運転状態に基づいて近似できる燃焼速度近似手段及び着火遅れ近似手段を燃焼速度と着火遅れの適合方法として有することを特徴とする。
この構成によれば、筒内圧推定装置の適合を容易に行うことが可能となる。
【００８９】
また、実施の形態１による内燃機関の筒内圧推定装置は、前記適合方法で生じる実燃焼速度の初期サージ部について、燃焼速度が一定の場合に燃焼量が同じとなるように着火遅れを補正することを特徴とする。
この構成によれば、着火遅れをより精度良く適合することが可能となる。
【００９０】
なお、この発明の実施の形態１による内燃機関の筒内圧推定装置について説明したが、この発明は、その発明の範囲内において、実施の形態を適宜、変形、省略することが可能である。例えば、ここで示した燃焼シミュレーション方法の要点は火炎を楕円体状と考えることであり、燃焼により発生する熱量の演算時に例えば化学反応式を考慮して熱解離（燃焼時等の高温下で水や二酸化炭素が部分的に逆反応を起こす現象）や酸化窒素（ＮＯｘ）の生成を考慮し、これらの吸熱反応による温度低下を考慮した燃焼モデルに拡張することも可能である。更には、筒内壁面へ熱が逃げる冷却損失を考慮した燃焼モデルに拡張することも可能であるので、この実施の形態として示した方法に限らず適用可能である。
【符号の説明】
【００９１】
１  エンジン、２  ＡＦＳ、３  電子制御スロットル、４  スロットル開度センサ、５  サージタンク、６  インテークマニホールド、７  インマニ圧センサ、８  吸気温センサ、９  ＥＧＲバルブ、１０  インジェクタ、１１  吸気ＶＶＴ、１２  排気ＶＶＴ、１３  点火コイル、１４  エキゾーストマニホールド、１５  回転プレート、１６  クランク角度センサ、１７  ＥＣＵ、１８  大気圧センサ、１９  各種センサ、２０  運転状態検出手段、２１  各種アクチュエータ、２２  燃焼パラメータ近似手段、２３  燃焼シミュレーション手段。
【図１】
2015096700.tif 000023
【図２】
2015096700.tif 000024
【図３】
2015096700.tif 000025
【図４】
2015096700.tif 000026
【図５】
2015096700.tif 000027
【図６】
2015096700.tif 000028
【図１０】
2015096700.tif 000029
【図１１】
2015096700.tif 000030
【図７】
2015096700.tif 000031
【図８】
2015096700.tif 000032
【図９】
2015096700.tif 000033