(57)【要約】【課題】内燃機関の運転条件、運転状況を反映した火花点火式内燃機関の最適点火時期設定を行うことである。【解決手段】最適点火時期を求める方法は、点火時期を初期値ＳＡ１から少しずつ進角して（Ｓ１０、Ｓ１２）、燃焼シミュレーションを実行して平均指圧特性曲線とその進角におけるトルクτを求め（Ｓ１４）、ノッキングが発生するか否かについて点火時期予測手法を用いて判断し（Ｓ１６，Ｓ１８）、ノッキングが発生する限界の点火時期をもって最適点火時期とする（Ｓ２２）ものである。ここで、点火時期予測手法（Ｓ１６）は、燃焼シミュレーションによって算出される平均指圧特性曲線から出発し、サイクル変動モデル、自着火予測モデル、ノック判定値変動モデルを用い、その点火時期についてノッキング発生の有無を予測する。【選択図】図１
2007332934.tif 000002
【特許請求の範囲】
【請求項１】
火花点火式内燃機関について、任意の運転条件の下における点火時期を与えて、燃焼シミュレーションを実行し、その点火時期の下で、クランク角度に対する平均指圧の関係を示す平均指圧特性曲線を算出し、算出された平均指圧特性曲線に対し、ノックする可能性のあるノック可能性指圧特性曲線を求める工程と、
ノック可能性指圧特性曲線を満足する条件の下で燃焼シミュレーションを実行し、クランク角度に対する質量燃焼割合特性曲線及び筒内温度特性曲線を算出する手段と、
ノック可能性指圧特性曲線及び筒内温度特性曲線に基づき、温度と圧力の関数である着火遅れ時間を用いる所定の自着火予測モデルに従って自着火時期を求める工程と、
質量燃焼割合特性曲線から自着火時期に対応する質量燃焼割合である自着火質量燃焼割合を求め、求められた自着火質量燃焼割合に基づきノック発生の有無を判断するノック有無判断工程と、
を含み、点火時期を変化させてノック発生の限界の点火時期を最適点火時期とすることを特徴とする火花点火式内燃機関の最適点火時期設定方法。
【請求項２】
請求項１に記載の火花点火式内燃機関の最適点火時期設定方法において、
ノック可能性指圧特性曲線を求める工程は、
平均指圧特性曲線に対し、予め記憶されているサイクル変動下での指圧特性曲線群を検索し、その中で最大の指圧特性を示す最大指圧特性曲線を特定し、これをノック可能性指圧特性曲線とすることを特徴とする火花点火式内燃機関の最適点火時期設定方法。
【請求項３】
請求項１に記載の火花点火式内燃機関の最適点火時期設定方法において、
自着火時期を求める工程は、
中間温度領域においてその前後の温度領域に比べ温度に関する着火遅れ時間の変化が少なくなる自着火予測モデルを用いて自着火時期を求めることを特徴とする火花点火式内燃機関の最適点火時期設定方法。
【請求項４】
請求項１に記載の火花点火式内燃機関の最適点火時期設定方法において、
着火遅れ時間は、指圧の大きさ、又は筒内温度、又は（理論空燃比／実際の空燃比）の値、又はオクタン価、又は残留ガス量の少なくとも１つの関数であるとして自着火時期を求めることを特徴とする火花点火式内燃機関の最適点火時期設定方法。
【請求項５】
請求項１に記載の火花点火式内燃機関の最適点火時期設定方法において、
ノック有無判断工程は、
求められた自着火質量燃焼割合と、予め定めたノック発生質量燃焼割合とを比較して、ノック発生の有無を判断することを特徴とする火花点火式内燃機関の最適点火時期設定方法。
【請求項６】
火花点火式内燃機関について、任意の運転条件の下において点火時期を与えて燃焼シミュレーションを実行し、点火時期とトルクの初期値を与える工程と、
初期値の点火時期を進角し、燃焼シミュレーションを実行して、その進角の下でのトルクと、クランク角度に対する平均指圧特性曲線を算出する進角工程と、
算出された平均指圧特性曲線に対し、予め記憶されているサイクル変動下での指圧特性曲線群を検索し、その中で最大の指圧特性を示す最大指圧特性曲線をノック可能性指圧特性曲線として特定する工程と、
ノック可能性指圧特性曲線を満足する条件の下で燃焼シミュレーションを実行し、クランク角度に対する質量燃焼割合特性曲線と、筒内温度特性曲線とを算出する工程と、
ノック可能性指圧特性曲線及び筒内温度特性曲線に基づき、温度と圧力の関数である着火遅れ時間を用いる所定の自着火予測モデルに従って自着火時期を求める工程と、
質量燃焼割合特性曲線から自着火時期に対応する質量燃焼割合を求め、予め定めたノック発生質量燃焼割合と比較し、ノック発生の有無を判断するノック有無判断工程と、
ノック有無判断工程において、ノックが発生すると判断する場合、及び進角工程で求められたトルクが進角前のトルク以下の場合には、進角工程で進角を行う前の点火時期を最適点火時期とする最適点火時期設定工程と、
を含み、ノック有無判断工程において、ノックが発生すると判断されない場合、進角工程に戻り、さらに進角して、以後の各工程を繰り返すことを特徴とする火花点火式内燃機関の最適点火時期設定方法。
【請求項７】
火花点火式内燃機関の最適点火時期をコンピュータによって求める最適点火時期設定装置であって、
火花点火式内燃機関について、任意の運転条件の下において点火時期を与えて燃焼シミュレーションを実行し、点火時期とトルクの初期値を与える処理手段と、
初期値の点火時期を進角し、燃焼シミュレーションを実行して、その進角の下でのトルクと、クランク角度に対する平均指圧特性曲線を算出する進角処理手段と、
算出された平均指圧特性曲線に対し、予め記憶されているサイクル変動下での指圧特性曲線群を検索し、その中で最大の指圧特性を示す最大指圧特性曲線をノック可能性指圧特性曲線として特定する処理手段と、
ノック可能性指圧特性曲線を満足する条件の下で燃焼シミュレーションを実行し、クランク角度に対する質量燃焼割合特性曲線と、筒内温度特性曲線とを算出する処理手段と、
ノック可能性指圧特性曲線及び筒内温度特性曲線に基づき、温度と圧力の関数である着火遅れ時間を用いる所定の自着火予測モデルに従って自着火時期を求める処理手段と、
質量燃焼割合特性曲線から自着火時期に対応する質量燃焼割合を求め、予め定めたノック発生質量燃焼割合と比較し、ノック発生の有無を判断するノック有無判断処理手段と、
ノック有無判断工程において、ノックが発生すると判断する場合、及び進角工程で求められたトルクが進角前のトルク以下の場合には、進角工程で進角を行う前の点火時期を最適点火時期とする最適点火時期設定処理手段と、
を含み、ノック有無判断工程において、ノックが発生すると判断されない場合、進角工程に戻り、さらに進角して、以後の各処理を繰り返すことを特徴とする火花点火式内燃機関の最適点火時期設定装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、火花点火式内燃機関の最適点火時期設定方法及び火花点火式内燃機関の最適点火時期設定装置に係り、特に計算モデルとシミュレーション等を用いて最適点火時期を求める火花点火式内燃機関の最適点火時期設定方法及び火花点火式内燃機関の最適点火時期設定装置に関する。
【背景技術】
【０００２】
火花点火式内燃機関においては、圧縮比を大きくすることが熱効率向上のために有効であるが、ノッキングの発生はその阻害要因の１つに数えられる。そこで機関にノックセンサを設けてノッキングの発生を検知し、点火時期にフィードバック制御を行う試みがなされており、すでに多くの機関にこの仕組みが取り入れられている。このようにノッキングの発生を検知してから回避制御を行うだけでなく、より確実にノッキングの回避を行うために、ノッキングの発生を計算モデルに従って予測することが行われている。なお、ノッキングは単にノックと呼ばれることも多く、以下では、ノッキングあるいはノックとして述べる。
【０００３】
例えば、特許文献１には、火花点火式内燃機関におけるノッキングの発生を計算モデルによって予測制御する燃焼制御装置が開示されている。ここでは、燃焼室について火炎面を境に未燃部と既燃部とに分割した２領域モデルを基本とし、機関の運転状態の情報から推定された各状態値を初期値とし、燃焼に関するエネルギー保存式等を連立させ、既燃部の質量割合をＷｉｅｂｅ関数形でクランク角度の関数として、クランク角度に応じて方程式を解いて、未燃部・既燃部の状態変化を予測することが述べられている。そして、自着火モデルにこのシミュレーションを適用して、未燃ガスの自着火発生時期及び自着火による発熱量を予測してノッキング強度を推定し、このノッキング強度が所定値以下となるように燃焼制御が行われる。
【０００４】
ここで自着火モデルとしては３つ述べられている。１つは、未燃ガスが所定温度・圧力状態にあるときの着火遅れ時間を与える式に基づき、ある時点での未燃ガスの温度・圧力から着火遅れ時間を求め、その遅れ時間経過後に自着火が発生するというものである。ここで用いられる着火遅れ時間の式は、温度の逆数（１／Ｔ）について指数関数であり、圧力について（−ｎ）乗の関数である。第２の自着火モデルは、時間経過と共に着火遅れ時間が変化することを考慮したＬｉｖｅｎｇｏｏｄ−Ｗｕ積分を用い、ある瞬間における温度・圧力等からその状態下での着火遅れ時間を求め、この逆数を時間推移に沿って積分し、積分値が１となったときに自着火が発生するというものである。第３の自着火モデルは、燃料室内のガスを化学種ｉが質量分率ｙｉの割合で混合したガスとし、これら化学種ｉの質量分率ｙｉが反応の進行に伴って変化するとする化学素反応モデルに基づくもので、これにより未燃ガスの温度・圧力の時間推移に伴う変化を予測し、未燃ガス部の温度が所定値以上になったときに自着火が発生するというものである。
【０００５】
【特許文献１】特開２００４−３３２５８４号公報
【発明の開示】
【発明が解決しようとする課題】
【０００６】
特許文献１では、いくつかの自着火モデルを用い、未燃ガスの自着火発生時期及び自着火による発熱量を予測してノッキング強度を推定し、このノッキング強度が所定値以下となるように燃焼制御が行われることが述べられている。ここでは機関の運転状態の情報から推定された各状態値を初期値としてシミュレーションが行われている。機関の運転状態でノッキングの発生に密接なものは、クランク角度に対する気筒内の圧力の関係で、これは一般的に指圧履歴あるいは単に指圧と呼ばれている。
【０００７】
特許文献１をはじめ、一般的な内燃機関の燃焼シミュレーションでは、計算の基礎となる指圧履歴として、平均指圧履歴が用いられている。平均指圧履歴は、内燃機関の運転条件を与えたときの典型的な指圧履歴であるが、実際の内燃機関の運転状況はかなりの変動があり、同じ運転条件の下での指圧履歴についてもかなりのサイクル変動がある。ノッキングは、指圧履歴が高い圧力を示す方が発生しやすいが、従来技術においては平均指圧履歴を用いているので、ノッキングの発生を低めに評価していることになる。すなわち、圧縮比の向上の面から見れば、安全側とは逆側にとっていることになる。
【０００８】
また、特許文献１において第１及び第２の自着火予測モデルは、未燃ガスが所定温度・圧力状態にあるときの着火遅れ時間を与える式に基づいて算出される着火遅れ時間を用いている。ここで用いられる着火遅れ時間の式は、温度の逆数（１／Ｔ）について指数関数であり、圧力について（−ｎ）乗の関数である。したがって、温度、すなわち筒内温度と、圧力すなわち指圧については上記の関数形を仮定して考慮されているが、実際の運転状況においては、運転条件によって自着火遅れ時間が変化し、また、自着火遅れ時間も上記のような関数形に従わない場合もある。
【０００９】
このように、従来技術では、内燃機関の実際の運転状況で生じるサイクル変動、自着火遅れ時間の状況を反映していないことがあり、自着火時期の推定が不正確となり、ノッキング抑制のための最適点火時期の設定が不正確となる。
【００１０】
本発明の目的は、内燃機関の運転条件、運転状況を反映した火花点火式内燃機関の最適点火時期設定方法及び火花点火式内燃機関の最適点火時期設定装置を提供することである。
【課題を解決するための手段】
【００１１】
本発明に係る火花点火式内燃機関の最適点火時期設定方法は、火花点火式内燃機関について、任意の運転条件の下における点火時期を与えて、燃焼シミュレーションを実行し、その点火時期の下で、クランク角度に対する平均指圧の関係を示す平均指圧特性曲線を算出し、算出された平均指圧特性曲線に対し、ノックする可能性のあるノック可能性指圧特性曲線を求める工程と、ノック可能性指圧特性曲線を満足する条件の下で燃焼シミュレーションを実行し、クランク角度に対する質量燃焼割合特性曲線及び筒内温度特性曲線を算出する手段と、ノック可能性指圧特性曲線及び筒内温度特性曲線に基づき、温度と圧力の関数である着火遅れ時間を用いる所定の自着火予測モデルに従って自着火時期を求める工程と、質量燃焼割合特性曲線から自着火時期に対応する質量燃焼割合である自着火質量燃焼割合を求め、求められた自着火質量燃焼割合に基づきノック発生の有無を判断するノック有無判断工程と、を含み、点火時期を変化させてノック発生の限界の点火時期を最適点火時期とすることを特徴とする。
【００１２】
また、本発明に係る火花点火式内燃機関の最適点火時期設定方法において、ノック可能性指圧特性曲線を求める工程は、平均指圧特性曲線に対し、予め記憶されているサイクル変動下での指圧特性曲線群を検索し、その中で最大の指圧特性を示す最大指圧特性曲線を特定し、これをノック可能性指圧特性曲線とすることが好ましい。
【００１３】
また、本発明に係る火花点火式内燃機関の最適点火時期設定方法において、自着火時期を求める工程は、中間温度領域においてその前後の温度領域に比べ温度に関する着火遅れ時間の変化が少なくなる自着火予測モデルを用いて自着火時期を求めることが好ましい。
【００１４】
また、本発明に係る火花点火式内燃機関の最適点火時期設定方法において、着火遅れ時間は、指圧の大きさ、又は筒内温度、又は（理論空燃比／実際の空燃比）の値、又はオクタン価、又は残留ガス量の少なくとも１つの関数であるとして自着火時期を求めることが好ましい。
【００１５】
また、本発明に係る火花点火式内燃機関の最適点火時期設定方法において、ノック有無判断工程は、求められた自着火質量燃焼割合と、予め定めたノック発生質量燃焼割合とを比較して、ノック発生の有無を判断することが好ましい。
【００１６】
また、本発明に係る火花点火式内燃機関の最適点火時期設定方法は、火花点火式内燃機関について、任意の運転条件の下において点火時期を与えて燃焼シミュレーションを実行し、点火時期とトルクの初期値を与える工程と、初期値の点火時期を進角し、燃焼シミュレーションを実行して、その進角の下でのトルクと、クランク角度に対する平均指圧特性曲線を算出する進角工程と、算出された平均指圧特性曲線に対し、予め記憶されているサイクル変動下での指圧特性曲線群を検索し、その中で最大の指圧特性を示す最大指圧特性曲線をノック可能性指圧特性曲線として特定する工程と、ノック可能性指圧特性曲線を満足する条件の下で燃焼シミュレーションを実行し、クランク角度に対する質量燃焼割合特性曲線と、筒内温度特性曲線とを算出する工程と、ノック可能性指圧特性曲線及び筒内温度特性曲線に基づき、温度と圧力の関数である着火遅れ時間を用いる所定の自着火予測モデルに従って自着火時期を求める工程と、質量燃焼割合特性曲線から自着火時期に対応する質量燃焼割合を求め、予め定めたノック発生質量燃焼割合と比較し、ノック発生の有無を判断するノック有無判断工程と、ノック有無判断工程において、ノックが発生すると判断する場合、及び進角工程で求められたトルクが進角前のトルク以下の場合には、進角工程で進角を行う前の点火時期を最適点火時期とする最適点火時期設定工程と、を含み、ノック有無判断工程において、ノックが発生すると判断されない場合、進角工程に戻り、さらに進角して、以後の各工程を繰り返すことを特徴とする。
【００１７】
また、本発明に係る火花点火式内燃機関の最適点火時期設定装置は、火花点火式内燃機関の最適点火時期をコンピュータによって求める最適点火時期設定装置であって、火花点火式内燃機関について、任意の運転条件の下において点火時期を与えて燃焼シミュレーションを実行し、点火時期とトルクの初期値を与える処理手段と、初期値の点火時期を進角し、燃焼シミュレーションを実行して、その進角の下でのトルクと、クランク角度に対する平均指圧特性曲線を算出する進角処理手段と、算出された平均指圧特性曲線に対し、予め記憶されているサイクル変動下での指圧特性曲線群を検索し、その中で最大の指圧特性を示す最大指圧特性曲線をノック可能性指圧特性曲線として特定する処理手段と、ノック可能性指圧特性曲線を満足する条件の下で燃焼シミュレーションを実行し、クランク角度に対する質量燃焼割合特性曲線と、筒内温度特性曲線とを算出する処理手段と、ノック可能性指圧特性曲線及び筒内温度特性曲線に基づき、温度と圧力の関数である着火遅れ時間を用いる所定の自着火予測モデルに従って自着火時期を求める処理手段と、質量燃焼割合特性曲線から自着火時期に対応する質量燃焼割合を求め、予め定めたノック発生質量燃焼割合と比較し、ノック発生の有無を判断するノック有無判断処理手段と、ノック有無判断工程において、ノックが発生すると判断する場合、及び進角工程で求められたトルクが進角前のトルク以下の場合には、進角工程で進角を行う前の点火時期を最適点火時期とする最適点火時期設定処理手段と、を含み、ノック有無判断工程において、ノックが発生すると判断されない場合、進角工程に戻り、さらに進角して、以後の各処理を繰り返すことを特徴とする。
【発明の効果】
【００１８】
上記構成の少なくとも１つにより、与えられた運転条件の下で平均指圧特性曲線を算出し、算出された平均指圧特性曲線に対し、ノックする可能性のあるノック可能性指圧特性曲線を求め、以後、このノック可能性指圧特性曲線に基づいて燃焼シミュレーションの実行等を行って自着火時期を求め、その時期における質量燃焼割合に基づいてノック有無を判断する。そして、点火時期を変化させてノック発生の限界の点火時期を最適点火時期とする。
【００１９】
このように、平均指圧特性曲線、すなわち平均指圧履歴を用いるのでなく、平均指圧履歴に対応するノック可能性指圧特性曲線を求めて、これについて自着火時期等を求める。したがって、内燃機関の実際の運転条件、運転状況を反映して最適点火時期を設定することが可能となる。
【００２０】
また、ノック可能性指圧特性曲線を求めるのは、平均指圧特性曲線に対し、予め記憶されているサイクル変動下での指圧特性曲線群を検索し、その中で最大の指圧特性を示す最大指圧特性曲線を特定して行われる。したがって、内燃機関の実際の運転条件、運転状況における指圧履歴のサイクル変動を反映して最適点火時期を設定することが可能となる。
【００２１】
また、中間温度領域においてその前後の温度領域に比べ温度に関する着火遅れ時間の変化が少なくなる自着火予測モデルを用いて自着火時期を求めるので、自着火遅れ時間を、単に、温度の逆数についての指数関数とするのでなく、内燃機関の実際の運転状況における自着火遅れ時間の変化を反映して最適点火時期を設定することが可能となる。
【００２２】
また、着火遅れ時間は、指圧の大きさ、又は筒内温度、又は（理論空燃比／実際の空燃比）の値、又はオクタン価、又は残留ガス量の少なくとも１つの関数であるとして自着火時期を求めるので、内燃機関の実際の運転状況を反映して最適点火時期を設定することが可能となる。
【００２３】
また、求められた自着火質量燃焼割合と、予め定めたノック発生質量燃焼割合とを比較して、ノック発生の有無を判断する。自着火時の質量燃焼割合が低いほどノッキングが発生しやすいことが知られており、運転条件及びその変動条件を定めれば、ノッキングが発生するときの質量燃焼割合を予め予測できる。この予測されたノック発生質量燃焼割合よりも求められた自着火質量燃焼割合が低ければノッキングが発生する。したがって、両者の比較で、ノッキング発生の限界の自着火時期を知ることができ、これを最適点火時期とすることで、内燃機関の実際の運転条件、運転状況を反映して最適点火時期を設定することが可能となる。
【００２４】
また、任意の運転条件の下において点火時期を与えて燃焼シミュレーションを実行し、点火時期とトルクの初期値を与え、その初期値から点火時期を進角させて、平均指圧特性曲線算出、最大指圧特性曲線特定、質量燃焼割合特性曲線算出、筒内温度特性曲線算出、自着火時期算出、ノック発生の有無判断等を行い、ノックが発生すると判断する場合、及び進角工程で求められたトルクが進角前のトルク以下の場合には、進角工程で進角を行う前の点火時期を最適点火時期とする。そして、ノックが発生すると判断されない場合、さらに進角して、以後の各工程を繰り返す。このように、点火時期をずらしていって、最大指圧特性曲線に基づいて最適点火時期を求めることができる。したがって、内燃機関の実際の運転条件、運転状況を反映して最適点火時期を設定することが可能となる。
【発明を実施するための最良の形態】
【００２５】
以下に図面を用いて本発明の実施の形態につき、詳細に説明する。以下において図示される特性曲線等は、説明の便宜のために用いるものであって、実際の内燃機関の運転条件、運転状況を反映するものであれば、これ以外の特性を有する特性曲線等であっても構わない。
【００２６】
図１は、最適点火時期を求める方法の原理的流れを説明するフローチャートであり、図２は、原理的流れの中で用いられるモデルの役割について説明するフローチャートである。なお、具体的にコンピュータ上で処理される手順のフローチャートについては後に詳述する。
【００２７】
図１に示されるように、最適点火時期を求める方法は、点火時期を少しずつ進角して、ノッキングが発生するか否かについて点火時期予測手法を用いて判断し、ノッキングが発生する限界の点火時期をもって最適点火時期とするものである。すなわち、対象となる内燃機関の与えられた運転条件の下で、点火時期の初期値ＳＡ１を与え、その点火時期において燃焼シミュレーションを実行してトルクの初期値τ１を算出する。図１では、与えられた点火時期の初期値がＳＡ０、算出されたトルクの初期値がτ０として示されている（Ｓ６，Ｓ８，Ｓ１０）。点火時期の初期値ＳＡ１は、ノッキングが発生しないことがわかっている点火時期にする。一般的に、点火時期とトルクの特性は、ある点火時期において最大のトルクを示すことが知られており、また、ノッキングは、点火時期をその最大トルクを示す点火時期に向かって進角するととともに発生しやすくなることが知られている。ここで、点火時期を進角させる、とは、クランクの回転において、先行する角度、すなわち、現在より前の時刻のクランク角度に対応する点火時期に設定することである。これに対し現在より後の時刻のクランク角度に対応する点火時期に設定することは、点火時期を遅角させると呼ばれる。したがって、初期値の点火時期ＳＡ１＝ＳＡ０は、点火時期を十分に遅角させることがよい。
【００２８】
そして、点火時期をΔＳＡだけ進角させ、ＳＡとする（Ｓ１２）。この進角させた状態のもとで燃焼シミュレーションを実行し、平均指圧特性、すなわち、各クランク角度に対する筒内圧力である各指圧の典型的値を求め、また、求められた平均指圧特性曲線から、点火時期ＳＡにおけるトルクτを求める（Ｓ１４）。このように、指圧特性とは、横軸にクランク角度、縦軸に指圧を取った特性曲線であり、クランクが回転するに伴い内燃機関の筒内圧力の変化履歴を示すものであるので、指圧履歴と呼ばれる。平均指圧特性曲線、すなわち、平均指圧履歴を求める燃焼シミュレーションプログラムは、すでに様々なものが提案されているので、それらの中で適当なものを用いることができる。
【００２９】
平均指圧特性曲線が求められると、これに基づいて点火時期予測手法を用い（Ｓ１６）、その進角ＳＡの下でノッキングが発生するか否かが判断され、これにより、点火時期をさらに進角するか否かが判定される（Ｓ１８）。そしてＳ１６の予測の結果、ノッキングが発生すると判断されると、進角前の点火時期ＳＡ１を最適点火時期ＳＡ＊とすることになる（Ｓ２２）。
【００３０】
上記のように、初期値ＳＡ１＝ＳＡ０の設定は、トルク最大を示す点火時期より十分に遅角されているので、進角の大きさΔＳＡが適当であれば、通常はノック無しと判断され、またそのときのトルクτは、進角前のときのトルクτ１より大きくなることになる。そこで、τがτ１より大きいか否かの判定（Ｓ２０）においてτがτ１より大きいと判断されない場合を除き、ノッキング無しと判定される場合、進角されたＳＡの状態を初期値に置き換え（Ｓ２４）、再びＳ１２に戻ることになる。すなわち、点火時期がさらにΔＳＡだけ進角される。そして以後の処理が繰り返される。点火時期が順次進角されると、次第にノッキングが発生しやすくなるので、ノッキングが初めて発生した進角状態の前の進角ＳＡ１の状態をもって、最適点火時期ＳＡ＊とされる。Ｓ２０においてτがτ１より大きいと判断されない場合も、進角状態の前の進角ＳＡ１の状態をもって、最適点火時期ＳＡ＊とされる。このようにして、点火時期を少しずつ進角して、ノッキングが発生するか否かについて点火時期予測手法を用いて判断し、ノッキングが発生する限界の点火時期をもって最適点火時期とすることができる。
【００３１】
図２には、点火時期予測手法（Ｓ１６）の詳細な内部フローチャートが示されている。点火時期予測手法は、燃焼シミュレーションによって算出される平均指圧特性曲線から出発し、サイクル変動モデル、自着火予測モデル、ノック判定値変動モデルを用い、その点火時期についてノッキング発生の有無を予測する手法である。
【００３２】
サイクル変動モデルとは、与えられた運転条件の下で、例えば４００サイクルの指圧履歴の変動をモデル化し、その運転条件の下での平均指圧特性曲線と対応付けたものである。このサイクル変動モデルは、運転条件ごとに作成されて適当な記憶装置に記憶され、平均指圧特性曲線を検索キーとして、検索されることができる。
【００３３】
自着火予測モデルは、指圧特性曲線に基づき燃焼シミュレーションを行って得られる質量燃焼割合と筒内温度の履歴を用い、自着火時期と、そのときにおける質量燃焼割合を求める際に用いられるもので、いつ自着火が生じるかのモデルである。ここでは、特許文献１の第２の自着火モデルを改良して用いる。特許文献１の第２の自着火モデルは、時間経過と共に着火遅れ時間が変化することを考慮したＬｉｖｅｎｇｏｏｄ−Ｗｕ積分を用い、ある瞬間における温度・圧力等からその状態下での着火遅れ時間を求め、この逆数を時間推移に沿って積分し、積分値が１となったときに自着火が発生するというものである。特許文献１で、着火遅れ時間ｆを規定する式は、温度の逆数（１／Ｔ）について指数関数であり、圧力について（−ｎ）乗の関数である。図２における自着火予測モデルは、後に詳細に説明するように、温度についての関数形を実際のものに適合させ、また、温度・圧力以外の運転条件パラメータの効果を考慮したものを用いる。これらのパラメータの効果を考慮した自着火予測モデルは、適当な記憶装置に記憶され、１又は複数のパラメータを検索キーとして、着火遅れ時間ｆを読み出すことができる。
【００３４】
ノック判定値変動モデルは、ノックが発生するときの質量燃焼割合をモデル化したものである。ここで、質量燃焼割合とは、（現在までの発熱量）／（投入された燃料量全体の低位発熱量）で定義される。ここで、低位発熱量とは、水分が気体（水蒸気）のままでいる場合の発熱量で、実際に利用できる熱量、いわゆる真発熱量に当る。その単位は、Ｊ／ｋｇである。投入された燃料量によって総発熱量が計算される。単位から理解されるように、発熱量は燃焼質量と関係付けられ、その観点から、質量燃焼割合と呼ばれるものである。質量燃焼割合は、ＢＭＦ（Ｂｕｒｎｅｄ  Ｍａｓｓ  Ｆｒａｃｔｉｏｎ）とも呼ばれる。
【００３５】
一般的に質量燃焼割合が低いほどノッキングが発生しやすいことが知られ、運転条件によって、ノッキングが発生する質量燃焼割合を求めておくことができる。そこで、ノッキングが発生する質量燃焼割合と、自着火時期における質量燃焼割合とを比較することで、ノック判定をすることができる。ノック判定変動モデルは、このような観点のモデルである。
【００３６】
図２に戻り、点火時期予測手法（Ｓ１６）においては、まず平均指圧特性曲線が与えられる（Ｓ３０）。次に、この平均指圧特性曲線に対応するサイクル変動モデルが取得される（Ｓ３２）。そして、取得されたサイクル変動モデルに含まれる複数の指圧特性曲線の中から、ノックする可能性のある指圧特性曲線を特定する。一般的には指圧が高いほどノッキングしやすいので、サイクル変動モデルに含まれる複数の指圧特性曲線の中の最大指圧特性曲線がノックする可能性のある指圧特性曲線として特定される（Ｓ３４）。
【００３７】
この最大指圧特性曲線の条件の下で、質量燃焼割合の履歴曲線と、筒内温度の履歴曲線、すなわちクランク角度に対する質量燃焼割合の変化を示す特性曲線と、同じくクランク角度に対する筒内温度の変化を示す特性曲線とが求められる（Ｓ３６）。
【００３８】
そして、最大指圧特性曲線と、筒内温度特性曲線とに基づいて、時間経過、すなわちクランク角度に対応する着火遅れ時間ｆが算出される。着火遅れ時間ｆは、上記のように、運転条件の各パラメータについてモデル化されているので、各パラメータを検索キーとして、着火遅れ時間ｆを時間経過と共に求めることができる。そして、Ｌｉｖｅｎｇｏｏｄ−Ｗｕ積分を実行し、積分値が１となる時期が自着火時期ｔｋである。そして、質量燃焼割合特性曲線上でｔｋに対応する質量燃焼割合ｐＢＭＦを求める（Ｓ３８）。
【００３９】
そして、ノック判定値変動モデルにしたがって、ノックするｋＢＭＦと、Ｓ３８で求められたｐＢＭＦとを比較し（Ｓ４２）、ノックの有無を判定する。このようにして、燃焼シミュレーションによって算出される平均指圧特性曲線から出発し、サイクル変動モデル、自着火予測モデル、ノック判定値変動モデルを用い、その点火時期についてノッキング発生の有無を予測することができる。
【００４０】
図３は、火花点火式内燃機関の最適点火時期設定装置１０の構成を示すブロック図である。最適点火時期設定装置１０は、計算モデルとシミュレーション等を用いて最適点火時期を求める機能を有する装置で、ＣＰＵ１２と、初期条件等を入力するキーボード等の入力部１４と、計算結果等を出力するプリンタ、ディスプレイ等の出力部１６と、計算モデル及びプログラム等を記憶する記憶部１８とを含んで構成される。これらの要素は相互に内部バスで接続される。かかる最適点火時期設定装置１０は、シミュレーションに適したパーソナルコンピュータあるいは科学計算用コンピュータ等で構成することができる。
【００４１】
記憶部１８には、図２で説明したサイクル変動モデル３４、自着火予測モデル３６、ノック判定値変動モデル３８、及び燃焼シミュレーションプログラム４０等が格納される。上記のように、サイクル変動モデル３４は、平均指圧特性曲線を検索キーとして、それに対応するサイクル変動による複数の指圧特性曲線を検索でき、特に最大指圧特性曲線を検索することができる。また自着火予測モデル３６は、筒内温度、指圧、運転状況のパラメータ等を検索キーとして、着火遅れ時間ｆを検索することができる。ノック判定値変動モデル３８は、運転条件等を検索キーとしてｋＢＭＦを検索することができる。燃焼シミュレーションプログラム４０としては、運転条件を与えることで平均指圧特性曲線を算出でき、また、指圧特性曲線に基づいて、筒内温度特性曲線、質量燃焼割合特性曲線等を算出する機能を有するものが格納される。
【００４２】
ＣＰＵ１２は、最適点火時期を求めるために、点火時期を少しずつ進角させて演算を進めさせる進角処理部２０、入力部１４から入力される運転条件を取得し、燃焼シミュレーションプログラムを実行させて平均指圧特性曲線を算出する平均指圧算出部２２、サイクル変動モデルを用い、平均指圧特性曲線に対応する最大指圧特性曲線を求める最大指圧算出部２４、最大指圧特性曲線を満足する条件の下で燃焼シミュレーションプログラムを実行させ、質量燃焼割合特性曲線と、筒内温度特性曲線とを求める質量燃焼割合・筒内温度算出部２６、自着火予測モデルを用いて自着火時期ｔｋを求める自着火時期算出部２８、ノック判定値変動モデルを用いてノック有無を判断するノック有無判断部３０、進角を繰り返し、ノッキングが生じる限界の点火時期を最適点火時期として設定する最適点火時期設定部３２とを含んで構成される。これらの機能は、ソフトウェアを実行することで実現でき、具体的には、対応する最適点火時期設定プログラムを実行することで実現できる。また、これらの機能の一部をハードウェアで実現するものとしてもよい。
【００４３】
上記構成の最適点火時期設定装置１０の作用、特に、ＣＰＵ１２の各機能に付き、図４のフローチャート、及び図５から図１５に示される各特性曲線に従って、以下に詳細に説明する。図４は、最適点火時期設定のための手順を示すフローチャートで、各手順は、対応する最適点火時期設定プログラムの各処理手順に該当する。
【００４４】
最適点火時期設定装置１０を立ち上げ、最適点火時期設定プログラムを立ち上げた後、入力部１４から、対象とする内燃機関の運転条件を入力する。運転条件としては、その内燃機関の諸元等のほか、燃料のオクタン価等が含まれる。そして、最適点火時期を求めるための出発点である点火時期の初期値ＳＡ１＝ＳＡ０が与えられる。そしてその条件の下で燃焼シミュレーションが実行され、トルクの初期値τ１＝τ０が算出される（Ｓ４６，Ｓ４８，Ｓ５０）。
【００４５】
その様子を図５に示す。図５は、横軸に点火時期ＳＡをとり、縦軸にトルクτを取って、点火時期ＳＡを変化させたときのトルクτの変化を示すもので、ある点火時期でトルクが最大になることが示されている。この点火時期−トルク特性曲線５０において、点火時期をこの最大トルクを示す時期より遅角側に設定すればするほどノッキングが起こりにくく、反対に、最大トルクを示す時期に向かって進角側に設定するほどノッキングが発生しやすくなることが知られている。そこで、ノッキングが発生する限界で最適点火時期を設定する計算を行う初期値としては、この最大トルクを示す時期より十分に遅角側に点火時期の初期値を設定することが好ましい。図５では、点火時期の初期値としてＳＡ１、トルクの初期値として、ＳＡ１におけるトルクτ１が示されている。なお、図５におけるＴＤＣは、クランクの上死点に相当する点火時期を表している。
【００４６】
次に、ＣＰＵ１２の進角処理部２０の機能により、点火時期の初期値ＳＡ１をΔＳＡだけ進角する（Ｓ５２）。その様子は、上記図５において矢印方向にて示されている。その進角した条件の下で、燃焼シミュレーションを実行し、平均指圧特性曲線と、その進角した点火時期におけるトルクτを求める（Ｓ５４）。この機能は、ＣＰＵ１２の平均指圧算出部２２によって実行される。具体的には、Ｓ５０で与えられた運転条件の下で点火時期をＳＡ＝ＳＡ１＋ΔＳＡとし、記憶部１８の燃焼シミュレーションプログラムを用いて燃焼シミュレーションを実行する。
【００４７】
燃焼シミュレーションの結果得られた平均指圧特性曲線５２の例を図６に示す。平均指圧特性曲線５２は、与えられた運転条件の下での指圧履歴の典型的なデータを示すものである。図６に示されるように、平均指圧特性曲線５２は、横軸にクランク角度を、縦軸に指圧を取った特性曲線である。クランク角度は、上死点（ＴＤＣ）を０°として、そこからの角度（Ａｆｔｅｒ  ＴＤＣ；ＡＴＤＣ）で定義される。なお、内燃機関のクランクの回転速度が与えられれば、クランク角度はそのまま時間に対応する。後述するように、着火遅れ時間ｆは時間の単位を有し、１／ｆは時間について積分するので、以後では機械的角度を示すクランク角度は、これに対応する時間としての性格を有するものとして扱う。すなわち、クランク角度に対応する時間とは、そのような意味で用いている。したがって、点火時期といえば時間を示すが、この時間はクランク角度に容易に変換できるので、点火時期をその時間に対応するクランク角度で示すこともできるのである。
【００４８】
平均指圧特性曲線５２が求まると、その点火時期ＳＡ＝ＳＡ１＋ΔＳＡにおけるトルクτを求め、点火時期ＳＡ１の下でのトルクτ１と比較する（Ｓ５６）。図５で説明したように、Ｓ５０における点火時期ＳＡ１は、点火時期−トルク特性におけるトルク最大となる点火時期より十分遅角側に設定されるので、ΔＳＡが適切であれば、進角後の点火時期におけるトルクτは、進角前のトルクτ１より大きくなる。したがって、１回目の進角後においては、通常であれば、Ｓ５６はＹＥＳと判断され、Ｓ５８に進む。
【００４９】
なお、後述のように、複数回進角を繰り返すと、進角後の点火時期におけるトルクτが次第に最大トルクに近づき、ついには最大トルクを超え、進角後の点火時期におけるトルクτが、進角前のトルクτ１より小さくなる。このときＳ５６はＮＯと判断され、Ｓ７４に進み、その最後の進角の１つ前の進角における点火時期が最適点火時期とされる。
【００５０】
Ｓ５６においてＹＥＳと判断されるとＳ５８に進み、平均指圧特性曲線５２の条件に対応するサイクル変動モデルが取得される。そして、取得されたサイクル変動モデルの複数の指圧特性曲線の中で最大指圧のものが、ノックする可能性の最も高い指圧特性として特定される。すなわち、平均指圧特性曲線５２に対応する最大指圧特性曲線が特定される（Ｓ６０）。この機能は、最大指圧算出部２４によって実行される。具体的には、記憶部１８のサイクル変動モデル３４において、Ｓ５４によって算出された平均指圧特性曲線５２を検索キーとして、それに対応するサイクル変動モデルを検索し、それに含まれる複数の指圧特性曲線データの中で、最大の指圧を有する指圧特性曲線がノックする可能性の高い最大指圧特性曲線として特定される。その様子を図７に示す。図７は、平均指圧特性曲線５２に対応するサイクル変動モデルに属する複数の指圧特性曲線を示すもので、その中で最大トルクを有するものが最大指圧特性曲線５４として示される。
【００５１】
以後の処理は、この最大指圧特性曲線５４に基づいて実行される。このように、平均指圧特性曲線５２を用いずに、最大指圧特性曲線５４を用いることで、実際の内燃機関の運転状況の変動を反映し、最もノッキングの発生しやすい条件の下で、ノッキング発生の限界の点火時期を最適点火時期とできる。すなわち、ノッキング発生の限界の評価をぎりぎりまで詰めることができて、その分、内燃機関の圧縮比を大きくし、熱効率向上を図りながら、ノッキングを抑制できる点火時期を設定することが可能となる。
【００５２】
最大指圧特性曲線５４が算出されると、その条件を満たすように燃焼シミュレーションが実行され、質量燃焼割合特性曲線と筒内温度特性曲線とが算出される（Ｓ６２）。質量燃焼割合特性曲線は、質量燃焼割合の時間経過に伴う変化を示すもので、具体的には、クランク角度に対する質量燃焼割合の変化として示される。なお、後述の図１４に質量燃焼割合特性曲線７０が示されている。また、筒内温度特性曲線は、内燃機関のシリンダ内温度の時間経過に伴う変化を示すもので、具体的には、クランク角度に対する筒内温度の変化として示される。なお、後述の図１２に筒内温度特性曲線６６が示されている。
【００５３】
筒内温度特性曲線が得られると、最大指圧特性曲線と合わせて、各クランク角度に対する圧力と温度とが得られるので、各クランク角度について着火遅れ時間ｆを求めることができる。具体的には、次のような手順が実行される。すなわち、まず、ｆ−１／Ｔ特性を求める（Ｓ６４）。ここでｆは上記のように着火遅れ時間であり、（１／Ｔ）は、筒内温度Ｔの逆数である。ｆ−１／Ｔ特性は、実際の運転条件、運転状況に関するパラメータごとに、記憶部１８の自着火予測モデル３６に記憶されているので、パラメータを検索キーとして、必要なｆ−１／Ｔ特性を読み出すことで取得することができる。
【００５４】
図８から図１１は、自着火予測モデル３６に記憶されている着火遅れ時間ｆに関する特性の例を示すものである。図８は、横軸に１／Ｔの１０００倍の値をとり、縦軸に対数目盛でｆをとって、ｆ−１／Ｔ特性５６を示したものである。ｆ−１／Ｔ特性５６は、図８に示されるように、中間温度領域においてその前後の温度領域に比べ、温度Ｔに対する着火遅れ時間ｆの変化が少なくなっている。これに対し、上記特許文献１においては、着火遅れ時間ｆを単純に（１／Ｔ）の指数関数として扱っており、このような中間温度領域における特徴的な特性を用いていない。なお、図８には、着火遅れ時間ｆに対する指圧、つまり圧力による効果を、特性曲線５８で示してある。ここでは、指圧が大きくなるほど着火遅れ時間ｆが短くなることが示されている。
【００５５】
図９は、着火遅れ時間ｆに対する当量比の効果を特性曲線６０で示すものである。横軸、縦軸のとり方は図８と同様である。ここで当量比は、（理論空燃比／実際の空燃比）で定義される値である。上記特許文献１において、着火遅れ時間ｆに対する当量比の効果は示されていないが、実際の運転状況の下では、図９に示されるように、当量比によって着火遅れ時間ｆが変化する。ここでは、当量比が大きくなるほど着火遅れ時間ｆが短くなることが示されている。
【００５６】
図１０は、着火遅れ時間ｆに対する燃料オクタン価の効果を特性曲線６２で示すものである。横軸、縦軸のとり方は図８と同様である。上記特許文献１において、着火遅れ時間ｆに対するオクタン価の効果は示されていないが、実際の運転状況の下では、図９に示されるように、オクタン価によって着火遅れ時間ｆが変化する。ここでは、オクタン価が大きくなるほど着火遅れ時間ｆが長くなることが示されている。
【００５７】
図１１は、着火遅れ時間ｆに対する筒内残留ガス量の効果を特性曲線６４で示すものである。ここでは、横軸に筒内温度Ｔをとり、縦軸に、着火遅れ時間ｆと、残留ガス割合が０であるときの着火遅れ時間ｆ０との比であるｆ／ｆ０をとってある。図１１に示されるように、残留ガス量の影響度ｆ／ｆ０は、ある筒内温度において最大となる分布特性を有し、また、残留ガス量が大きいほど、着火遅れ時間に対する影響度ｆ／ｆ０が大きくなる。
【００５８】
これらの着火遅れ時間ｆに関するデータを用いて、１／ｆをバルブ閉に対応するクランク角度に対応する時間から上死点側の方向に任意のクランク角度に対応する時間に渡って時間積分し（Ｓ６６）、その時間積分値が１となるときの時間又はこれに対応するクランク角度を自着火時期ｔｋとする（Ｓ６８）。これは、上記のように、Ｌｉｖｅｎｇｏｏｄ−Ｗｕ積分として知られる手法である。この機能は、自着火時期算出部２８の機能により実行される。
【００５９】
その様子は図１２、図１３に示される。すなわち、図１２は、Ｓ６２によって求められた筒内温度特性曲線６６を示すもので、横軸にクランク角度、縦軸に筒内温度Ｔの１／１０００の値がとられている。また、Ｓ６０において算出された最大指圧特性曲線５４も、同じ横軸に対応させて示されている。この筒内温度特性曲線６６の各クランク角度における温度Ｔと、最大指圧特性曲線５４の各クランク角度における指圧と、実際の運転条件、運転状況下における当量比、オクタン価、残留ガス量とを用い、図８から図１１で説明した着火遅れ時間ｆの特性を参照して、実際の運転状況下における各クランク角度における着火遅れ時間ｆを求める。そして各クランク角度を対応する時間に変換し、着火遅れ時間ｆの逆数を時間積分する。積分範囲は、上記のようにバルブ閉に対応するクランク角度に対応する時間が始点で、上死点側の方向に任意のクランク角度に対応する時間が終点である。図１３は、１／ｆの時間積分曲線６８の様子を示すもので、積分範囲の始点から上死点の方向に向かって時間積分するに従い積分値が大きくなる。そこで、積分値が１となる時期を、自着火時期ｔｋとする（Ｓ６８）。
【００６０】
次に、Ｓ６２で求めた質量燃焼割合特性曲線において、自着火時期ｔｋに対応する値であるｐＢＭＦを求める（Ｓ７０）。そして、ノック判定値変動モデルに従い、予め定めてあるノック発生質量燃焼割合ｋＢＭＦと比較し、ノックの有無を判定する（Ｓ７２）。この機能は、ノック有無判断部３０により実行される。具体的には、記憶部１８からｋＢＭＦを読み出し、Ｓ７０で求められたｐＢＭＦとの大小関係を判断し、ｐＢＭＦがｋＢＭＦ以下であるか否かを判断し、判断がＹＥＳであればノック有りとし、ＮＯであればノック無しとする。
【００６１】
その様子を図１４に示す。図１４は、Ｓ６２で求められた質量燃焼割合特性曲線７０を示す図で、横軸にクランク角度、縦軸に質量燃焼割合ＢＭＦをとって示してある。自着火時期ｔｋに対応するクランク角度における質量燃焼割合がｐＢＭＦであり、図１４にｋＢＭＦとして示されているものがノック発生質量燃焼割合である。図１４の例では、ｐＢＭＦがｋＢＭＦより大きいので、Ｓ７２の判断はＮＯとなる。
【００６２】
Ｓ７２で判断がＮＯということは、この点火時期においては、まだノックが発生していないので、進角処理部２０の機能により、再びＳ５２に戻り、点火時期をΔＳＡだけ進角し、さらにＳ５４以下の手順を実行する。このように、点火時期が順次進角されると、次第にノッキングが発生しやすくなるので、ノッキングが初めて発生した進角状態の前の進角ＳＡ１の状態をもって、最適点火時期ＳＡ＊とされる。また、Ｓ５６に関連して説明したように、Ｓ５６においてτがτ１より大きいと判断されない場合も、進角状態の前の進角ＳＡ１の状態をもって、最適点火時期ＳＡ＊とされる（Ｓ７４）。この機能は、最適点火時期設定部３２によって実行される。このようにして、点火時期を少しずつ進角して、ノッキングが発生するか否かについて点火時期予測手法を用いて判断し、ノッキングが発生する限界の点火時期をもって最適点火時期とすることができる。
【００６３】
なお、図１５は、図１２、図１３、図１４の内容を、横軸を共通として、１枚の図にまとめたもので、総合的な点火時期最適化の様子を示す図である。
【図面の簡単な説明】
【００６４】
【図１】本発明に係る実施の形態において、最適点火時期を求める方法の原理的流れを説明するフローチャートである。
【図２】本発明に係る実施の形態において、点火時期予測手法の詳細な内部フローチャートである。
【図３】本発明に係る実施の形態において、火花点火式内燃機関の最適点火時期設定装置の構成を示すブロック図である。
【図４】本発明に係る実施の形態において、最適点火時期設定のための手順を示すフローチャートである。
【図５】本発明に係る実施の形態において、点火時期−トルク特性曲線を示す図である。
【図６】本発明に係る実施の形態において、平均指圧特性曲線の例を示す図である。
【図７】本発明に係る実施の形態において、平均指圧特性曲線に対応する最大指圧特性曲線を示す図である。
【図８】本発明に係る実施の形態において、着火遅れ時間ｆと筒内温度Ｔの関係を示す図である。
【図９】本発明に係る実施の形態において、着火遅れ時間ｆに対する当量比の効果を示す図である。
【図１０】本発明に係る実施の形態において、着火遅れ時間ｆに対する燃料オクタン価の効果を示す図である。
【図１１】本発明に係る実施の形態において、着火遅れ時間ｆに対する筒内残留ガス量の効果を示す図である。
【図１２】本発明に係る実施の形態において、筒内温度特性曲線の例を示す図である。
【図１３】本発明に係る実施の形態において、着火遅れ時間ｆの逆数の時間積分と、積分値が１となる時期を自着火時期ｔｋとする様子を示す図である。
【図１４】本発明に係る実施の形態において、自着火時期ｔｋにおける質量燃焼割合ｐＢＭＦとノック発生質量燃焼割合ｋＢＭＦとの関係を説明する図である。
【図１５】本発明に係る実施の形態において、総合的な点火時期最適化の様子を示す図である。
【符号の説明】
【００６５】
１０  最適点火時期設定装置、１２  ＣＰＵ、１４  入力部、１６  出力部、１８  記憶部、２０  進角処理部、２２  平均指圧算出部、２４  最大指圧算出部、２６  質量燃焼割合・筒内温度算出部、２８  自着火時期算出部、３０  ノック有無判断部、３２  最適点火時期設定部、３４  サイクル変動モデル、３６  自着火予測モデル、３８  ノック判定値変動モデル、４０  燃焼シミュレーションプログラム、５０  点火時期−トルク特性曲線、５２  平均指圧特性曲線、５４  最大指圧特性曲線、５６  ｆ−１／Ｔ特性曲線、５８，６０，６２，６４  着火遅れ時間ｆに対する効果を示す特性曲線、６６  筒内温度特性曲線、６８  １／ｆの時間積分曲線、７０  質量燃焼割合特性曲線。
【図１】
2007332934.tif 000003
【図２】
2007332934.tif 000004
【図３】
2007332934.tif 000005
【図４】
2007332934.tif 000006
【図５】
2007332934.tif 000007
【図６】
2007332934.tif 000008
【図７】
2007332934.tif 000009
【図８】
2007332934.tif 000010
【図９】
2007332934.tif 000011
【図１０】
2007332934.tif 000012
【図１１】
2007332934.tif 000013
【図１２】
2007332934.tif 000014
【図１３】
2007332934.tif 000015
【図１４】
2007332934.tif 000016
【図１５】
2007332934.tif 000017