(57)【特許請求の範囲】
【請求項１】
【数１】
0004577211.tif 000017
（ただし、
Ｑ[J]：内燃機関の気筒内の発熱量
Ｑtotal[J]：気筒内に供給された燃料の熱量
θ[degCA]：熱発生開始後の経過クランク角度
θｐ[degCA]：燃焼期間
ｍ：形状パラメータ
ｋ：効率）
上記（Ｉ）式またはこれと等価な式で表されるWiebe関数のパラメータを決定する方法であって、
内燃機関の気筒内におけるクランク角度毎の実熱発生率を求めるステップと、
実熱発生率が最大となったクランク角度での実燃焼割合αを求めるステップと、
前記実燃焼割合αに基づいて、下記（II）式またはこれと等価な式により、形状パラメータｍの値を決定するステップと、
を備えることを特徴とするWiebe関数パラメータの決定方法。
【数２】
0004577211.tif 000018
【請求項２】
実熱発生率の積分値を総発熱量ｋＱtotalとして求めるステップと、
実熱発生率の最大値をβとしたとき、この実熱発生率最大値βと、前記総発熱量ｋＱtotalと、前記形状パラメータｍの値とに基づいて、下記（III）式またはこれと等価な式により、ａ／θｐｍ＋１の値を決定するステップと、
を更に備えることを特徴とする請求項１記載のWiebe関数パラメータの決定方法。
【数３】
0004577211.tif 000019
【請求項３】
実熱発生率が最大となったときのクランク角度をθATDC＊[degATDC]とし、内燃機関の実点火時期をSA[degBTDC]としたとき、前記クランク角度θATDC＊と、前記実点火時期SAと、前記形状パラメータｍの値と、前記ａ／θｐｍ＋１の値とに基づいて、下記（IV）および（Ｖ）式、またはこれらと等価な式により、実点火時期とWiebe関数の熱発生開始点とのズレ量を表す熱発生開始点ズレ量θｂ[degCA]の値を決定するステップを更に備えることを特徴とする請求項２記載のWiebe関数パラメータの決定方法。
【数４】
0004577211.tif 000020
【請求項４】
実熱発生率の積分値を総発熱量ｋＱtotalとして求めるステップと、
気筒内に供給された燃料の熱量Ｑtotalの値を求めるステップと、
前記総発熱量ｋＱtotalを前記Ｑtotalの値で除することにより、効率ｋの値を決定するステップと、
を更に備えることを特徴とする請求項１乃至３の何れか１項記載のWiebe関数パラメータの決定方法。
【請求項５】
【数５】
0004577211.tif 000021
（ただし、
Ｑ[J]：内燃機関の気筒内の発熱量
Ｑtotal[J]：気筒内に供給された燃料の熱量
θ[degCA]：熱発生開始後の経過クランク角度
θｐ[degCA]：燃焼期間
ｍ：形状パラメータ
ｋ：効率）
上記（VI）式またはこれと等価な式で表されるWiebe関数のパラメータを決定する装置であって、
内燃機関の気筒内におけるクランク角度毎の実熱発生率を求める実熱発生率取得手段と、
実熱発生率が最大となったクランク角度での実燃焼割合αを求める実燃焼割合取得手段と、
前記実燃焼割合αに基づいて、下記（VII）式またはこれと等価な式により、形状パラメータｍを決定するパラメータ決定手段と、
を備えることを特徴とするWiebe関数パラメータの決定装置。
【数６】
0004577211.tif 000022
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、Wiebe関数パラメータの決定方法および決定装置に係り、内燃機関の気筒内の熱発生を表すWiebe関数のパラメータを決定するのに好適なWiebe関数パラメータの決定方法および決定装置に関する。
【背景技術】
【０００２】
近年、内燃機関の気筒内の熱発生をWiebe関数によって推定する試みがなされている。Wiebe関数によれば、クランク角度毎の熱発生率や燃焼割合等を算出することができ、更に、それらの値から筒内圧や図示トルク等を算出することができる。ゆえに、気筒内の熱発生をWiebe関数によって正確に推定することができれば、Wiebe関数による推定値を内燃機関の制御等に利用することができる。
【０００３】
内燃機関の気筒内の熱発生をWiebe関数によって推定するためには、Wiebe関数に含まれるいくつかのパラメータの値を精度良く決定する必要がある。それらのパラメータの値は、空燃比、機関回転数、負荷率、点火時期などの運転条件によって異なり、また、エンジンの機種によっても異なるものと考えられる。
【０００４】
特開２００４−３３２６５９号公報には、Wiebe関数を利用して点火時期をＭＢＴ（Minimum advance for the Best Torque）付近に制御する技術が開示されている。この公報では、同一エンジンのＭＢＴにおいては、燃焼期間の長さにかかわらず熱発生率の形状が同じになり、かつ熱発生率が最大となるときのクランク角度および燃焼割合も一定になることを仮定している。そして、Wiebe関数パラメータのうち、熱発生率形状に関与するａおよびｍを定数とし、それらの値を実験で求めることとしている（上記公報の段落番号００３０参照）。なお、上記公報でのパラメータ「ｎ」が本明細書における「ｍ」に相当する。
【０００５】
【特許文献１】特開２００４−３３２６５９号公報
【発明の開示】
【発明が解決しようとする課題】
【０００６】
しかしながら、上記公報には、Wiebe関数パラメータのａおよびｍを実験データから求める具体的な方法については何ら開示されていない。また、本発明者の知見によれば、同一エンジンのＭＢＴにおいても、ａやｍの値は運転条件によって変化し、一定とはならない。また、種々の制御が行われる現代のエンジンにおいては、必ずしも点火時期をＭＢＴに合わせればよいとは限らない。このため、ＭＢＴ以外の点火時期の下でのトルク等を推定する要請も存在し、そのためには点火時期の変化によるａやｍの値の変化傾向を把握する必要がある。
【０００７】
このように、Wiebe関数を内燃機関の制御等に実際に利用するには、様々な運転条件の下でWiebe関数のパラメータを精度良く決定して、それらの値が運転条件に応じてどのように変化するかを正確に把握することが必要である。しかしながら、従来、Wiebe関数のパラメータを精度良く決定する手法は確立されていないのが実状である。
【０００８】
この発明は、上述のような課題を解決するためになされたもので、内燃機関の気筒内の熱発生パターンを表すためのWiebe関数のパラメータを精度良く決定することのできるWiebe関数パラメータの決定方法および決定装置を提供することを目的とする。
【課題を解決するための手段】
【０００９】
第１の発明は、上記の目的を達成するため、
【数１】
0004577211.tif 000002
（ただし、
Ｑ[J]：内燃機関の気筒内の発熱量
Ｑtotal[J]：気筒内に供給された燃料の熱量
θ[degCA]：熱発生開始後の経過クランク角度
θｐ[degCA]：燃焼期間
ｍ：形状パラメータ
ｋ：効率）
上記（Ｉ）式またはこれと等価な式で表されるWiebe関数のパラメータを決定する方法であって、
内燃機関の気筒内におけるクランク角度毎の実熱発生率を求めるステップと、
実熱発生率が最大となったクランク角度での実燃焼割合αを求めるステップと、
前記実燃焼割合αに基づいて、下記（II）式またはこれと等価な式により、形状パラメータｍの値を決定するステップと、
を備えることを特徴とする。
【数２】
0004577211.tif 000003
【００１０】
また、第２の発明は、第１の発明において、
実熱発生率の積分値を総発熱量ｋＱtotalとして求めるステップと、
実熱発生率の最大値をβとしたとき、この実熱発生率最大値βと、前記総発熱量ｋＱtotalと、前記形状パラメータｍの値とに基づいて、下記（III）式またはこれと等価な式により、ａ／θｐｍ＋１の値を決定するステップと、
を更に備えることを特徴とする。
【数３】
0004577211.tif 000004
【００１１】
また、第３の発明は、第２の発明において、
実熱発生率が最大となったときのクランク角度をθATDC＊[degATDC]とし、内燃機関の実点火時期をSA[degBTDC]としたとき、前記クランク角度θATDC＊と、前記実点火時期SAと、前記形状パラメータｍの値と、前記ａ／θｐｍ＋１の値とに基づいて、下記（IV）および（Ｖ）式、またはこれらと等価な式により、実点火時期とWiebe関数の熱発生開始点とのズレ量を表す熱発生開始点ズレ量θｂ[degCA]の値を決定するステップを更に備えることを特徴とする。
【数４】
0004577211.tif 000005
【００１２】
また、第４の発明は、第１乃至第３の発明の何れかにおいて、
実熱発生率の積分値を総発熱量ｋＱtotalとして求めるステップと、
気筒内に供給された燃料の熱量Ｑtotalの値を求めるステップと、
前記総発熱量ｋＱtotalを前記Ｑtotalの値で除することにより、効率ｋの値を決定するステップと、
を更に備えることを特徴とする。
【００１３】
また、第５の発明は、
【数５】
0004577211.tif 000006
（ただし、
Ｑ[J]：内燃機関の気筒内の発熱量
Ｑtotal[J]：気筒内に供給された燃料の熱量
θ[degCA]：熱発生開始後の経過クランク角度
θｐ[degCA]：燃焼期間
ｍ：形状パラメータ
ｋ：効率）
上記（VI）式またはこれと等価な式で表されるWiebe関数のパラメータを決定する装置であって、
内燃機関の気筒内におけるクランク角度毎の実熱発生率を求める実熱発生率取得手段と、
実熱発生率が最大となったクランク角度での実燃焼割合αを求める実燃焼割合取得手段と、
前記実燃焼割合αに基づいて、下記（VII）式またはこれと等価な式により、形状パラメータｍを決定するパラメータ決定手段と、
を備えることを特徴とする。
【数６】
0004577211.tif 000007
【発明の効果】
【００１４】
第１の発明によれば、内燃機関の気筒内での実際の熱発生パターンに、Wiebe関数が高い精度で適合するように、形状パラメータｍの値を精度良く決定することができる。すなわち、最小二乗法や近似計算などの処理を行う必要がないので、誤差が生じにくく、高精度に形状パラメータｍを算出することができる。また、演算負荷が小さく、容易に算出処理を行うことができる。
【００１５】
第２の発明によれば、内燃機関の気筒内での実際の熱発生パターンに、Wiebe関数が更に高い精度で適合するように、ａ／θｐｍ＋１の値を精度良く決定することができる。
【００１６】
第３の発明によれば、内燃機関の気筒内での実際の熱発生パターンに、Wiebe関数が更に高い精度で適合するように、熱発生開始点ズレ量θｂの値を精度良く決定することができる。
【００１７】
第４の発明によれば、内燃機関の気筒内での実際の熱発生パターンに、Wiebe関数が更に高い精度で適合するように、効率ｋの値を精度良く決定することができる。
【００１８】
第５の発明によれば、内燃機関の気筒内での実際の熱発生パターンに、Wiebe関数が高い精度で適合するように、形状パラメータｍの値を精度良く決定することができる。すなわち、最小二乗法や近似計算などの処理を行う必要がないので、誤差が生じにくく、高精度に形状パラメータｍを算出することができる。また、演算負荷が小さく、容易に算出処理を行うことができる。
【発明を実施するための最良の形態】
【００１９】
実施の形態１．
［システム構成の説明］
図１は、本発明の実施の形態１で用いられるシステム構成を説明するための図である。図１に示すように、本実施形態のシステムは、火花点火式の内燃機関１０を備えている。内燃機関１０には、クランク角を検出するクランク角センサ１２が組み込まれている。クランク角センサ１２は、クランク軸が所定回転角だけ回転する毎に、Hi出力とLo出力を反転させるセンサである。クランク角センサ１２の出力によれば、クランク角度（クランク軸の回転位置）や、機関回転数NEなどを検知することができる。また、内燃機関１０には、筒内圧センサ１８が組み込まれている。筒内圧センサ１８によれば、筒内（燃焼室内）に生ずる圧力を検出することができる。
【００２０】
内燃機関１０の吸気通路１９の途中には、サージタンク２０が設けられている。サージタンク２０には、その内部の圧力、すなわち吸気管圧力を検出する吸気圧センサ２１が設置されている。吸気圧センサ２１の出力によれば、内燃機関１０の負荷率KL［％］などを取得することができる。
【００２１】
また、吸気通路１９には、その内部を流通する吸入空気量GAを検出するエアフローメータ２２が配置されている。エアフローメータ２２の下流には、スロットル弁２４が配置されている。スロットル弁２４の近傍には、スロットル開度TAを検出するスロットル開度センサ２６が組み付けられている。
【００２２】
内燃機関１０の吸気ポートには、ガソリン等の燃料を噴射するための燃料噴射弁２８が配設されている。なお、本発明における内燃機関は、このようなポート噴射式の内燃機関に限定されるものではなく、筒内直接噴射式の内燃機関でも良く、あるいはポート噴射と筒内噴射とを併用する内燃機関でも良い。
【００２３】
また、内燃機関１０には、燃焼室内の混合気に点火するための点火プラグ３０が設置されている。更に、内燃機関１０の排気通路３２には、排気ガスの空燃比を検出する空燃比センサ３３が設置されている。また、排気通路３２には、排気ガスを浄化するための触媒３４が組み込まれている。
【００２４】
本実施形態のシステムは、制御装置およびパラメータ決定装置としてのECU(Electronic Control Unit)５０を備えている。ECU５０には、上述した各種のセンサからセンサ信号が供給されている。ECU５０は、それらのセンサ信号に基づいて、燃料噴射弁２８や、点火プラグ３０などの各種アクチュエータを制御することができる。更に、ECU５０は、クランク角センサ１２、筒内圧センサ１８等のセンサ信号に基づいて、後述するWiebe関数パラメータの決定処理を行うことができる。
【００２５】
［Wiebe関数］
Wiebe関数によれば、内燃機関の気筒内の熱発生率dＱ／dθおよび燃焼割合Ｑ／ｋＱtotalは、それぞれ、下記（１）および（２）式で表される。
【数７】
0004577211.tif 000008
上記（１）および（２）式中の各記号の意味は、下記の通りである。
Ｑ[J]：内燃機関の気筒内の発熱量
Ｑtotal[J]：気筒内に供給された燃料の熱量
θ[degCA]：熱発生開始後の経過クランク角度
θｐ[degCA]：燃焼期間
ｍ：形状パラメータ
ｋ：効率
【００２６】
図２は、クランク角度毎の熱発生率を示す図である。図２中、実線で示すグラフは、上記（１）式のWiebe関数による熱発生率dＱ／dθであり、破線で示すグラフは、実験により求められた実熱発生率である。なお、実熱発生率は、後述するように、筒内圧の実測データに基づいて算出することができる。
【００２７】
本発明によれば、後述する方法によってWiebe関数のパラメータを精度良く決定することができ、その結果、図２に示すように、実熱発生率をWiebe関数によって精度良く近似することができる。以下、Wiebe関数の各パラメータについて説明する。
【００２８】
（燃焼期間θｐ）
図２に示すように、燃焼期間θｐは、燃焼による熱の発生が継続する期間をクランク角度で表したものとしての物理的意味を有している。
【００２９】
（形状パラメータｍ）
形状パラメータｍは、図２のようなグラフにおける熱発生率dＱ／dθの形状に大きく関与するパラメータであり、特に、熱発生率が最大値をとるクランク角度に大きく関与するパラメータである。
【００３０】
（熱発生開始点ズレ量θｂ）
図２に示すように、上記Wiebe関数においては、クランク角度θ＝０のときに熱発生率dＱ／dθ＝０であり、θが０を超えたところから熱発生が開始する。つまり、本実施形態では、Wiebe関数上での熱発生開始点をθ＝０としており、θは熱発生開始後の経過クランク角度を表す値である。
【００３１】
本発明者の知見によれば、Wiebe関数の熱発生開始点を実際の点火時期SAに一致させることとすると、Wiebe関数を実熱発生率に必ずしも精度良く適合させることができない。そこで、本実施形態では、Wiebe関数の熱発生開始点（θ＝０）を実際の点火時期SAと異なる点とすることとした。そして、その両者のズレ量を熱発生開始点ズレ量θｂとして表し、これをWiebe関数のパラメータとして導入することとした。すなわち、熱発生開始点ズレ量θｂ[degCA]は、Wiebe関数の熱発生開始点と、実点火時期SAとのズレ量をクランク角度で表すパラメータである（図２参照）。なお、本実施形態のWiebe関数は、０≦θ≦θｐの範囲で定義されるものであるので、この熱発生開始点ズレ量θｂは、上記（１）および（２）式中には表れない。
【００３２】
（効率ｋ）
上記（１）および（２）式中、気筒内に供給された燃料の熱量Ｑtotalは、筒内に供給された燃料量に、その燃料の低位発熱量を乗ずることにより算出することができる。なお、低位発熱量は、真発熱量とも呼ばれる物性値である。低位発熱量とは、単位量の燃料が完全燃焼したときに発生する熱量から、燃料中に含まれる水分および燃焼によって生じる水分を蒸発させるのに必要な熱量（潜熱）を差し引いた残りの熱量を意味する。図１に示すシステムにおいては、Ｑtotalの値は、燃料の低位発熱量が既知であるものとして、燃料噴射弁２８からの燃料噴射量に基づいて算出することができる。あるいは、空燃比A/Fおよび筒内空気量（負荷率KL）からＱtotalを算出することもできる。
【００３３】
内燃機関における燃焼では、冷却損失や燃料の燃え残りなどに起因して、何らかの熱損失を伴うのが普通である。このため、気筒内に供給された燃料の熱量Ｑtotalが全部そのまま気筒内の発熱量Ｑに変換されることは実際上はあり得ない。本実施形態では、このことをWiebe関数に反映させるパラメータとして、効率ｋを導入することとした。すなわち、効率ｋは、気筒内に供給された燃料の熱量Ｑtotalが発熱量Ｑに変換される効率としての物理的意味を有しており、０＜ｋ＜１なる範囲の数である。
【００３４】
効率ｋを用いると、気筒内での総発熱量は、Ｑtotalにｋを乗じたｋＱtotalとして表すことができる。総発熱量ｋＱtotalは、図２において、実熱発生率を示す破線と、熱発生率dＱ／dθ＝０の直線とがグラフ上で囲む面積に相当する。つまり、総発熱量ｋＱtotalの値は、実熱発生率を積分することによって算出することが可能である。そのようにして算出された総発熱量ｋＱtotalを、上述の方法で算出したＱtotalで除することにより、効率ｋの値を求めることができる。
【００３５】
上記のように総発熱量がｋＱtotalとして表されるので、あるクランク角度θにおける燃焼割合は、そのクランク角度θまでの発熱量Ｑと、総発熱量ｋＱtotalとの比Ｑ／ｋＱtotalとして表すことができる。この燃焼割合Ｑ／ｋＱtotalは、Wiebe関数によれば、上記（２）式で表される。上記（２）式をθで微分すると、上記（１）式が得られる。この意味で、上記（２）式は、上記（１）式と等価である。
【００３６】
（パラメータａ）
燃焼割合の本来の定義からすると、燃焼割合の値は、燃焼終了時に１００％、すなわち１となるべきである。これに対し、Wiebe関数では、その性質上、燃焼終了時に燃焼割合の値が１とならない。すなわち、Wiebe関数では、上記（２）式にθ＝θｐを代入して得られる通り、燃焼終了時の燃焼割合の値は、下記式で表される。この下記式の右辺の値は１とすることができない。
Ｑ／ｋＱtotal（θ＝θｐ）＝１−exp（−ａ）
【００３７】
このように、パラメータａは、Wiebe関数上での燃焼終了時の燃焼割合の値に関与するパラメータであり、図２のようなグラフにおける熱発生率dＱ／dθの形状にも関与するパラメータである。
【００３８】
従来、次のようにして、このパラメータａを定数と仮定する方法がある。すなわち、Wiebe関数上での燃焼終了時の燃焼割合の値｛１−exp（−ａ）｝を１に近い所定の値（例えば０．９９９）と等しいと仮定すれば、ａは定数（例えば６．９）として求まる。しかしながら、本発明者の知見によれば、このようにしてａを定数とする方法では、Wiebe関数による熱発生率dＱ／dθのグラフ上での形状を、実熱発生率の形状に必ずしも精度良く適合させることはできない。このため、ａの値は、定数とすべきでなく、機種や運転条件に応じて変化するものであると考えられる。
【００３９】
本発明では、後述するように、パラメータａを含むａ／θｐｍ＋１の値を決定することができる。このａ／θｐｍ＋１の値を決定することで、機種や運転条件の違いによるａの値の変動をWiebe関数に精度良く反映させることができる。
【００４０】
（形状パラメータｍの決定方法）
本実施形態では、熱発生率が最大となったときの燃焼割合に着目することにより、形状パラメータｍを決定することができる。熱発生率が最大となったときの燃焼割合は、図２に示すように、熱発生率のグラフが囲む全部の面積に対する、熱発生率が最大となるクランク角度までの面積（図２中の実線ハッチング部分）の割合で表される。
【００４１】
まず、上記（１）式をθで微分し、d２Ｑ／dθ２＝０とおくことにより、Wiebe関数による熱発生率dＱ／dθが最大となるクランク角度θ＊を求める。下記の式展開を行うことにより、クランク角度θ＊を表す式が得られる。
【数８】
0004577211.tif 000009
【００４２】
上記（３）式で表されるクランク角度θ＊を上記（２）式に代入することにより、Wiebe関数による熱発生率dＱ／dθが最大となるクランク角度θ＊での燃焼割合は、下記（４）式で表すことができる。
【数９】
0004577211.tif 000010
【００４３】
上記（４）式に示されるように、熱発生率が最大となるクランク角度での燃焼割合は、Wiebe関数上では形状パラメータｍのみの関数で表される。そこで、これに対応する実測値を求め、その実測値を上記（４）式の右辺と等しいとおけば、ｍを決定することができる。
【００４４】
すなわち、実熱発生率が最大となったクランク角度での実燃焼割合をαとすると、このαを上記（４）式の右辺と等しいとおくことにより、形状パラメータｍを下記（５）式で決定することができる。
【数１０】
0004577211.tif 000011
【００４５】
（ａ／θｐｍ＋１の決定方法）
上記（３）式で表されるクランク角度θ＊を上記（１）式に代入することにより、Wiebe関数による熱発生率の最大値（dＱ／dθ）maxは、下記（６）式で表すことができる。
【数１１】
0004577211.tif 000012
【００４６】
よって、実熱発生率の最大値をβとすると、この最大値βが上記（６）式の下段と等しいとおくことにより、下記式が得られる。
【数１２】
0004577211.tif 000013
【００４７】
上記（７）式の右辺において、ａ／θｐｍ＋１以外の部分は、総発熱量ｋＱtotalおよび形状パラメータｍのみで表されている。総発熱量ｋＱtotalの値は、前述したように、実熱発生率を積分することによって算出することが可能である。そして、ｍの値は、上記（５）式により決定することができる。よって、それらの値と、実熱発生率最大値βとを上記（７）式に代入することにより、ａ／θｐｍ＋１の値を決定することができる。
【００４８】
（熱発生開始点ズレ量θｂの決定方法）
図３は、熱発生開始点ズレ量θｂを決定する方法を説明するための図である。本実施形態では、Wiebe関数による熱発生率が最大となる位置と、実熱発生率が最大となる位置とを一致させることにより、熱発生開始点ズレ量θｂを決定する。すなわち、実熱発生率が最大となったときのクランク角度をθATDC＊[degATDC]としたとき、Wiebe関数上で熱発生率が最大となるクランク角度θ＊の位置がθATDC＊に一致するようにする。
【００４９】
図３から判るように、実点火時期をSA[degBTDC]とし、θ＊とθATDC＊との位置が一致するものとすると、SA、θATDC＊、θｂ、およびθ＊の間に、下記（９）式が成り立つ。また、上記（３）式は、下記（８）式のように変形することができる。
【数１３】
0004577211.tif 000014
【００５０】
上記（８）式によれば、先に決定されたｍとａ／θｐｍ＋１との値に基づいてクランク角度θ＊を算出することができる。そして、そのθ＊と、実点火時期SAと、実測データから得られるθATDC＊とを上記（９）式に代入することにより、熱発生開始点ズレ量θｂを決定することができる。
【００５１】
以上の方法により、ｍ、ｋ、ａ／θｐｍ＋１、およびθｂの値をそれぞれ決定することができる。ここで、上記（１）式のWiebe関数は、下記式のように変形することができる。
【数１４】
0004577211.tif 000015
【００５２】
上記（１０）式によれば、ｍ、ｋ、ａ／θｐｍ＋１、およびθｂの値が決定されれば、Wiebe関数を具体的に記述することができることが判る。すなわち、本実施形態では、Wiebe関数パラメータとして、ｍ、ｋ、ａ／θｐｍ＋１、およびθｂの４つの値を決定すればよく、燃焼期間θｐおよびパラメータａの値は個別には決定しなくてもよいことが判る。
【００５３】
［実施の形態１における具体的処理］
図４は、上述した原理に基づいてWiebe関数パラメータを決定するために本実施形態においてECU５０が実行するルーチンのフローチャートである。図４に示すルーチンによれば、まず、内燃機関１０の運転中に、クランク角センサ１２および筒内圧センサ１８の出力に基づいて、クランク角度毎の筒内圧（燃焼圧）Ｐ[Pa]が計測される（ステップ１００）。このステップ１００において、筒内圧Ｐを計測するクランク角度間隔は、特に限定されないが、例えば１degCA間隔とすることができる。また、筒内圧Ｐの計測は、全クランク角度範囲（０〜７２０°）に渡って行わなくてもよく、実際の燃焼期間が含まれる範囲で行えばよい。
【００５４】
次に、上記ステップ１００で実測されたクランク角度毎の筒内圧Ｐに基づいて、下記の方法により、クランク角度毎の実熱発生率が算出される（ステップ１０２）。熱力学上のエネルギー保存則によれば、実熱発生率dＱ／dθは、筒内圧Ｐと筒内容積Ｖ[m3]とを用いて、次式で表される。なお、次式中のクランク角度θは、前記Wiebe関数におけるθではなく、１サイクルが０〜７２０°で表される通常のクランク角度を示すものとする。
【数１５】
0004577211.tif 000016
【００５５】
上記（１１）式中、κは比熱比であり、燃焼ガスの組成などに基づいて定まる既知の値である。また、筒内容積Ｖおよびその変化率dＶ／dθは、クランク角度θに応じて幾何学的に決まる既知の値である。よって、それらの既知の値と、上記ステップ１００で実測されたクランク角度毎の筒内圧Ｐの値とを上記（１１）式に代入することにより、クランク角度毎の実熱発生率を算出することができる。
【００５６】
次いで、総発熱量ｋＱtotalが算出される（ステップ１０４）。具体的には、上記ステップ１０２で算出されたクランク角度毎の実熱発生率を積分（積算）することにより、総発熱量ｋＱtotalが算出される。
【００５７】
次いで、効率ｋの値が決定される（ステップ１０６）。具体的には、まず、燃料噴射弁２８からの燃料噴射量等に基づいて、気筒内に供給された燃料の熱量Ｑtotalが算出される。そして、上記ステップ１０４で算出された総発熱量ｋＱtotalをこのＱtotalで除することにより、効率ｋが算出される。
【００５８】
次いで、上記ステップ１０２で算出されたクランク角度毎の実熱発生率のうちの最大値およびそのときのクランク角度を見つけ出すことにより、実熱発生率最大値βと、そのときのクランク角度θATDC＊とが取得される（ステップ１０８）。
【００５９】
次いで、実熱発生率が最大となったクランク角度θATDC＊での実燃焼割合αが算出される（ステップ１１０）。具体的には、上記ステップ１０２で算出されたクランク角度毎の実熱発生率をクランク角度θATDC＊まで積分（積算）した値を、上記ステップ１０４で算出された総発熱量ｋＱtotalで除することにより、実燃焼割合αが算出される。
【００６０】
次いで、形状パラメータｍの値が決定される（ステップ１１２）。具体的には、上記ステップ１１０で算出された実燃焼割合αを上記（５）式に代入することにより、ｍが算出される。
【００６１】
次いで、ａ／θｐｍ＋１の値が決定される（ステップ１１４）。具体的には、上記ステップ１０４で算出された総発熱量ｋＱtotalと、上記ステップ１０８で取得された実熱発生率最大値βと、上記ステップ１１２で算出された形状パラメータｍとを上記（７）式に代入することにより、ａ／θｐｍ＋１が算出される。
【００６２】
次いで、熱発生開始点ズレ量θｂの値が決定される（ステップ１１６）。具体的には、ECU５０が設定している実点火時期SAと、上記ステップ１０８で取得された実熱発生率が最大となるクランク角度θATDC＊と、上記ステップ１１２で算出された形状パラメータｍと、上記ステップ１１４で算出されたａ／θｐｍ＋１とを上記（８）および（９）式に代入することにより、熱発生開始点ズレ量θｂが算出される。
【００６３】
以上説明したルーチンにより、内燃機関１０での実測データに基づいて、Wiebe関数パラメータであるｍ、ｋ、ａ／θｐｍ＋１、およびθｂの値をそれぞれ決定することができる。本発明によれば、内燃機関１０での実測データから上記の数式による計算処理を行うことによって直接的に各Wiebe関数パラメータを算出することができる。このため、誤差が生じ易い最小二乗法や近似計算などを行う必要がない。それゆえ、各Wiebe関数パラメータを極めて精度良く決定することができ、実機での熱発生パターンに高精度に適合したWiebe関数を得ることができる。また、Wiebe関数パラメータの算出に要する演算負荷も少なく、容易に各Wiebe関数パラメータを算出することができる。
【００６４】
各Wiebe関数パラメータの値は、内燃機関１０の運転条件（点火時期SA、機関回転数NE、空燃比A/F、負荷率KL等）によって異なるものと考えられる。本発明によれば、例えば内燃機関１０の制御ロジックの開発段階において、図４に示すルーチンに従って各Wiebe関数パラメータを決定する操作を、複数の異なる運転条件に対して行うことにより、各Wiebe関数パラメータと、点火時期SA、機関回転数NE、空燃比A/F、負荷率KL等の運転条件パラメータとの関係を正確に把握することができる。そして、そのようにして把握した関係を利用すれば、あらゆる運転条件に対してWiebe関数を精度良く決定することができ、そのWiebe関数を利用して内燃機関１０の筒内圧Ｐやトルクなどを精度良く予測（推定）することができる。よって、内燃機関１０の制御ロジックの検討などを行う際、従来は多数の運転条件について筒内圧Ｐやトルクなどを実測する作業が必要であるが、本発明によればそのような実測作業をWiebe関数に基づく推定処理で置き換えることができる。このため、実機データの計測点数を大幅に削減することができ、開発工数、開発コストを大幅に削減することができる。
【００６５】
また、上記のようにして把握された各Wiebe関数パラメータと運転条件との関係は、マップや演算式の形にして、車載用ECUに記憶させておくこともできる。そのようにすれば、内燃機関１０を車両に搭載した状態で、種々の運転条件下における内燃機関１０の筒内圧Ｐやトルクをシミュレーションにより精度良く予測することができる。そして、その予測結果を利用して、例えばトルクデマンド制御やドライバビリティ改善制御などを精度良く行うことが可能となる。
【００６６】
また、図４に示すようなルーチンは、車載用ECUに実行させることも可能である。例えば、図４に示すようなルーチンを車載用ECUによって実行することによって算出されたｍ、ｋ、ａ／θｐｍ＋１、θｂの値を学習することにより、ECUに記憶されたそれらのマップを、経時変化や個体差に応じて補正する制御を行うことができる。
【００６７】
また、本実施形態では、前述したように、（ａ）実熱発生率最大値β、（ｂ）実熱発生率が最大値βをとったときのクランク角度θATDC＊、（ｃ）実熱発生率が最大値βをとったときの実燃焼割合α、に基づいて各Wiebe関数パラメータを決定している。本発明者の知見によれば、これら（ａ）〜（ｃ）の値は、点火時期SA、機関回転数NE、空燃比A/F、負荷率KL等の運転条件に応じて連続的に変化することが実験的に確かめられている。この点からも、本実施形態のWiebe関数パラメータの決定方法が極めて妥当であることが認められる。
【００６８】
なお、本実施形態では、前記のような各種の数式を利用したが、本発明では、それらの数式は前記のものと完全に同じである必要はなく、それらと等価な数式を利用してもよい。また、本発明では、ｍ、ｋ、ａ／θｐｍ＋１、θｂの全部の値を前述した方法によって決定しなくても良く、ｋ、ａ／θｐｍ＋１、θｂの値については他の方法によって決定してもよい。
【００６９】
また、上述した実施の形態１においては、ECU５０が、上記ステップ１００および１０２の処理を実行することにより前記第５の発明における「実熱発生率取得手段」が、上記ステップ１０４，１０８および１１０の処理を実行することにより前記第５の発明における「実燃焼割合取得手段」が、上記ステップ１１２の処理を実行することにより前記第５の発明における「パラメータ決定手段」が、それぞれ実現されている。
【図面の簡単な説明】
【００７０】
【図１】本発明の実施の形態１で用いられるシステム構成を説明するための図である。
【図２】クランク角度毎の熱発生率を示す図である。
【図３】熱発生開始点ズレ量θｂを決定する方法を説明するための図である。
【図４】本発明の実施の形態１において実行されるルーチンのフローチャートである。
【符号の説明】
【００７１】
１０  内燃機関
１２  クランク角センサ
１８  筒内圧センサ
２１  吸気圧センサ
３０  点火プラグ
３３  空燃比センサ
３４  触媒
５０  ECU(Electronic Control Unit)
【図１】
0004577211.tif 000023
【図２】
0004577211.tif 000024
【図３】
0004577211.tif 000025
【図４】
0004577211.tif 000026