(57)【要約】【課題】この発明は、内燃機関の筒内圧推定方法および筒内圧推定装置に関し、細かいクランク角刻みで筒内圧推定値を精度良く求めることを目的とする。【解決手段】本発明の内燃機関の筒内圧推定方法は、内燃機関の吸気弁閉時期と排気弁開時期との間の複数点の筒内圧を取得する筒内圧取得ステップ１３２と、筒内圧取得ステップ１３２で取得された筒内圧に基づいて、筒内の熱発生パターンと筒内圧との関係を表すモデル式に含まれる複数の未定パラメータを決定するパラメータ決定ステップ１３４，１３６と、パラメータ決定ステップ１３４，１３６で決定されたパラメータが代入されたモデル式に基づいて、所定クランク角での筒内圧推定値を算出する筒内圧推定ステップ１４０と、を備える。【選択図】図３
2008255932.tif 000002
【特許請求の範囲】
【請求項１】
内燃機関の吸気弁閉時期と排気弁開時期との間の複数点の筒内圧を取得する筒内圧取得ステップと、
前記筒内圧取得ステップで取得された筒内圧に基づいて、筒内の熱発生パターンと筒内圧との関係を表すモデル式に含まれる複数の未定パラメータを決定するパラメータ決定ステップと、
前記パラメータ決定ステップで決定されたパラメータが代入された前記モデル式に基づいて、所定クランク角での筒内圧推定値を算出する筒内圧推定ステップと、
を備えることを特徴とする内燃機関の筒内圧推定方法。
【請求項２】
前記モデル式は、Wiebe関数に基づくものであり、
前記複数のパラメータには、形状パラメータm、効率パラメータk、燃焼期間θpおよび熱発生開始クランク角θbが含まれることを特徴とする請求項１記載の内燃機関の筒内圧推定方法。
【請求項３】
前記筒内圧取得ステップは、
クランク角速度を取得するクランク角速度取得ステップと、
前記内燃機関の機械損失により生ずる機械損失トルクを算出する機械損失トルク算出ステップと、
前記内燃機関の運動部材により生ずる慣性質量トルクを算出する慣性質量トルク算出ステップと、
前記クランク角速度と、前記機械損失トルクと、前記慣性質量トルクとに基づいて、筒内圧により生ずる筒内圧トルクを算出する筒内圧トルク算出ステップと、
筒内容積変化率を算出する筒内容積変化率算出ステップと、
前記筒内圧トルクと前記筒内容積変化率とから筒内圧を算出する筒内圧算出ステップと、
を含むことを特徴とする請求項１または２記載の内燃機関の筒内圧推定方法。
【請求項４】
前記内燃機関は、クランク角を検出するクランク角センサを備え、
前記筒内圧取得ステップでは、前記クランク角センサの計測刻み毎の筒内圧を取得することを特徴とする請求項１乃至３の何れか１項記載の内燃機関の筒内圧推定方法。
【請求項５】
前記内燃機関は、複数気筒を有し、
前記筒内圧取得ステップは、
クランク角速度を取得するクランク角速度取得ステップと、
前記内燃機関の機械損失により生ずる機械損失トルクを算出する機械損失トルク算出ステップと、
前記内燃機関の運動部材により生ずる慣性質量トルクを算出する慣性質量トルク算出ステップと、
前記クランク角速度と、前記機械損失トルクと、前記慣性質量トルクとに基づいて、全気筒の筒内圧により生ずる全気筒筒内圧トルクを算出する全気筒筒内圧トルク算出ステップと、
各気筒の筒内容積変化率を算出する筒内容積変化率算出ステップと、
前記全気筒筒内圧トルクと、前記各気筒の筒内容積変化率とから、気筒別筒内圧を算出する気筒別筒内圧算出ステップと、
を含み、
前記パラメータ決定ステップでは、前記気筒別筒内圧に基づいて、前記未定パラメータを気筒別に決定し、
前記筒内圧推定ステップでは、気筒別に決定された前記パラメータが代入された気筒別の前記モデル式に基づいて、所定クランク角での筒内圧推定値を気筒別に算出することを特徴とする請求項１乃至４の何れか１項記載の内燃機関の筒内圧推定方法。
【請求項６】
前記気筒別筒内圧算出ステップでは、気筒別筒内圧が所定時間変化しないとの仮定に基づいて、前記気筒別筒内圧を算出することを特徴とする請求項５記載の内燃機関の筒内圧推定方法。
【請求項７】
前記内燃機関は、クランク角を検出するクランク角センサを備え、
前記クランク角速度取得ステップは、前記クランク角センサの計測刻み毎に取得されるクランク角速度を補間することにより、前記計測刻みより小さい刻み幅の補間後クランク角速度を算出する補間後クランク角速度算出ステップを含み、
前記気筒別筒内圧算出ステップでは、気筒別筒内圧が、前記小さい刻み幅に応じた所定の時間変化しないとの仮定に基づいて、前記気筒別筒内圧を算出することを特徴とする請求項５記載の内燃機関の筒内圧推定方法。
【請求項８】
内燃機関の吸気弁閉時期と排気弁開時期との間の複数点の筒内圧を取得する筒内圧取得手段と、
前記筒内圧取得手段により取得された筒内圧に基づいて、筒内の熱発生パターンと筒内圧との関係を表すモデル式に含まれる複数の未定パラメータを決定するパラメータ決定手段と、
前記パラメータ決定手段で決定されたパラメータが代入された前記モデル式に基づいて、所定クランク角での筒内圧推定値を算出する筒内圧推定手段と、
を備えることを特徴とする内燃機関の筒内圧推定装置。
【請求項９】
前記モデル式は、Wiebe関数に基づくものであり、
前記複数のパラメータには、形状パラメータm、効率パラメータk、燃焼期間θpおよび熱発生開始クランク角θbが含まれることを特徴とする請求項８記載の内燃機関の筒内圧推定装置。
【請求項１０】
前記筒内圧取得手段は、
クランク角速度を取得するクランク角速度取得手段と、
前記内燃機関の機械損失により生ずる機械損失トルクを算出する機械損失トルク算出手段と、
前記内燃機関の運動部材により生ずる慣性質量トルクを算出する慣性質量トルク算出手段と、
前記クランク角速度と、前記機械損失トルクと、前記慣性質量トルクとに基づいて、筒内圧により生ずる筒内圧トルクを算出する筒内圧トルク算出手段と、
筒内容積変化率を算出する筒内容積変化率算出手段と、
前記筒内圧トルクと前記筒内容積変化率とから筒内圧を算出する筒内圧算出手段と、
を含むことを特徴とする請求項８または９記載の内燃機関の筒内圧推定装置。
【請求項１１】
前記内燃機関は、クランク角を検出するクランク角センサを備え、
前記筒内圧取得手段は、前記クランク角センサの計測刻み毎の筒内圧を取得することを特徴とする請求項８乃至１０の何れか１項記載の内燃機関の筒内圧推定装置。
【請求項１２】
前記内燃機関は、複数気筒を有し、
前記筒内圧取得手段は、
クランク角速度を取得するクランク角速度取得手段と、
前記内燃機関の機械損失により生ずる機械損失トルクを算出する機械損失トルク算出手段と、
前記内燃機関の運動部材により生ずる慣性質量トルクを算出する慣性質量トルク算出手段と、
前記クランク角速度と、前記機械損失トルクと、前記慣性質量トルクとに基づいて、全気筒の筒内圧により生ずる全気筒筒内圧トルクを算出する全気筒筒内圧トルク算出手段と、
各気筒の筒内容積変化率を算出する筒内容積変化率算出手段と、
前記全気筒筒内圧トルクと、前記各気筒の筒内容積変化率とから、気筒別筒内圧を算出する気筒別筒内圧算出手段と、
を含み、
前記パラメータ決定手段は、前記気筒別筒内圧に基づいて、前記未定パラメータを気筒別に決定し、
前記筒内圧推定手段は、気筒別に決定された前記パラメータが代入された気筒別の前記モデル式に基づいて、所定クランク角での筒内圧推定値を気筒別に算出することを特徴とする請求項８乃至１１の何れか１項記載の内燃機関の筒内圧推定装置。
【請求項１３】
前記気筒別筒内圧算出手段は、気筒別筒内圧が所定時間変化しないとの仮定に基づいて、前記気筒別筒内圧を算出することを特徴とする請求項１２記載の内燃機関の筒内圧推定装置。
【請求項１４】
前記内燃機関は、クランク角を検出するクランク角センサを備え、
前記クランク角速度取得手段は、前記クランク角センサの計測刻み毎に取得されるクランク角速度を補間することにより、前記計測刻みより小さい刻み幅の補間後クランク角速度を算出する補間後クランク角速度算出手段を含み、
前記気筒別筒内圧算出手段は、気筒別筒内圧が、前記小さい刻み幅に応じた所定の時間変化しないとの仮定に基づいて、前記気筒別筒内圧を算出することを特徴とする請求項１２記載の内燃機関の筒内圧推定装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、内燃機関の筒内圧推定方法および筒内圧推定装置に関する。
【背景技術】
【０００２】
特開２００４−２４５１７３号公報には、火花点火式エンジンの運転条件から、圧縮上死点後の所定クランク位置に設定した基準クランク角での筒内圧力と未燃ガス温度とを求め、これら筒内圧力と未燃ガス温度とを用いてノッキング指標値を算出する技術が開示されている。
【０００３】
【特許文献１】特開２００４−２４５１７３号公報
【特許文献２】特開２００５−２３３１１０号公報
【発明の開示】
【発明が解決しようとする課題】
【０００４】
上記従来技術では、基準クランク角での筒内圧を運転条件から推定している。しかしながら、この方法では、筒内圧を精度良く推定することは困難である。
【０００５】
また、近年では、例えば燃料性状検出、失火検出、筒内状態量（燃焼割合等）をベースとする制御などを実現するために、筒内圧を細かいクランク角の刻み幅で精度良く求めたいという要望がある。
【０００６】
しかしながら、クランク角センサの計測刻みは、通常、１０度あるいは３０度と粗い。また、クランク角センサは、センサ特性の個体差や、ロータ寸法の個体差などの影響により、ノイズも大きい。このため、細かいクランク角刻みで筒内圧を精度良く求めることは困難である。
【０００７】
また、筒内圧センサは高価であるので、なるべくならば筒内圧センサを設置することなしに筒内圧を求めたいという要望もある。
【０００８】
この発明は、上述のような課題を解決するためになされたもので、細かいクランク角刻みで筒内圧推定値を精度良く求めることのできる内燃機関の筒内圧推定方法および筒内圧推定装置を提供することを目的とする。
【課題を解決するための手段】
【０００９】
第１の発明は、上記の目的を達成するため、内燃機関の筒内圧推定方法であって、
内燃機関の吸気弁閉時期と排気弁開時期との間の複数点の筒内圧を取得する筒内圧取得ステップと、
前記筒内圧取得ステップで取得された筒内圧に基づいて、筒内の熱発生パターンと筒内圧との関係を表すモデル式に含まれる複数の未定パラメータを決定するパラメータ決定ステップと、
前記パラメータ決定ステップで決定されたパラメータが代入された前記モデル式に基づいて、所定クランク角での筒内圧推定値を算出する筒内圧推定ステップと、
を備えることを特徴とする。
【００１０】
また、第２の発明は、第１の発明において、
前記モデル式は、Wiebe関数に基づくものであり、
前記複数のパラメータには、形状パラメータm、効率パラメータk、燃焼期間θpおよび熱発生開始クランク角θbが含まれることを特徴とする。
【００１１】
また、第３の発明は、第１または第２の発明において、
前記筒内圧取得ステップは、
クランク角速度を取得するクランク角速度取得ステップと、
前記内燃機関の機械損失により生ずる機械損失トルクを算出する機械損失トルク算出ステップと、
前記内燃機関の運動部材により生ずる慣性質量トルクを算出する慣性質量トルク算出ステップと、
前記クランク角速度と、前記機械損失トルクと、前記慣性質量トルクとに基づいて、筒内圧により生ずる筒内圧トルクを算出する筒内圧トルク算出ステップと、
筒内容積変化率を算出する筒内容積変化率算出ステップと、
前記筒内圧トルクと前記筒内容積変化率とから筒内圧を算出する筒内圧算出ステップと、
を含むことを特徴とする。
【００１２】
また、第４の発明は、第１乃至第３の発明の何れかにおいて、
前記内燃機関は、クランク角を検出するクランク角センサを備え、
前記筒内圧取得ステップでは、前記クランク角センサの計測刻み毎の筒内圧を取得することを特徴とする。
【００１３】
また、第５の発明は、第１乃至第４の発明の何れかにおいて、
前記内燃機関は、複数気筒を有し、
前記筒内圧取得ステップは、
クランク角速度を取得するクランク角速度取得ステップと、
前記内燃機関の機械損失により生ずる機械損失トルクを算出する機械損失トルク算出ステップと、
前記内燃機関の運動部材により生ずる慣性質量トルクを算出する慣性質量トルク算出ステップと、
前記クランク角速度と、前記機械損失トルクと、前記慣性質量トルクとに基づいて、全気筒の筒内圧により生ずる全気筒筒内圧トルクを算出する全気筒筒内圧トルク算出ステップと、
各気筒の筒内容積変化率を算出する筒内容積変化率算出ステップと、
前記全気筒筒内圧トルクと、前記各気筒の筒内容積変化率とから、気筒別筒内圧を算出する気筒別筒内圧算出ステップと、
を含み、
前記パラメータ決定ステップでは、前記気筒別筒内圧に基づいて、前記未定パラメータを気筒別に決定し、
前記筒内圧推定ステップでは、気筒別に決定された前記パラメータが代入された気筒別の前記モデル式に基づいて、所定クランク角での筒内圧推定値を気筒別に算出することを特徴とする。
【００１４】
また、第６の発明は、第５の発明において、
前記気筒別筒内圧算出ステップでは、気筒別筒内圧が所定時間変化しないとの仮定に基づいて、前記気筒別筒内圧を算出することを特徴とする。
【００１５】
また、第７の発明は、第５の発明において、
前記内燃機関は、クランク角を検出するクランク角センサを備え、
前記クランク角速度取得ステップは、前記クランク角センサの計測刻み毎に取得されるクランク角速度を補間することにより、前記計測刻みより小さい刻み幅の補間後クランク角速度を算出する補間後クランク角速度算出ステップを含み、
前記気筒別筒内圧算出ステップでは、気筒別筒内圧が、前記小さい刻み幅に応じた所定の時間変化しないとの仮定に基づいて、前記気筒別筒内圧を算出することを特徴とする。
【００１６】
また、第８の発明は、内燃機関の筒内圧推定装置であって、
内燃機関の吸気弁閉時期と排気弁開時期との間の複数点の筒内圧を取得する筒内圧取得手段と、
前記筒内圧取得手段により取得された筒内圧に基づいて、筒内の熱発生パターンと筒内圧との関係を表すモデル式に含まれる複数の未定パラメータを決定するパラメータ決定手段と、
前記パラメータ決定手段で決定されたパラメータが代入された前記モデル式に基づいて、所定クランク角での筒内圧推定値を算出する筒内圧推定手段と、
を備えることを特徴とする。
【００１７】
また、第９の発明は、第８の発明において、
前記モデル式は、Wiebe関数に基づくものであり、
前記複数のパラメータには、形状パラメータm、効率パラメータk、燃焼期間θpおよび熱発生開始クランク角θbが含まれることを特徴とする。
【００１８】
また、第１０の発明は、第８または第９の発明において、
前記筒内圧取得手段は、
クランク角速度を取得するクランク角速度取得手段と、
前記内燃機関の機械損失により生ずる機械損失トルクを算出する機械損失トルク算出手段と、
前記内燃機関の運動部材により生ずる慣性質量トルクを算出する慣性質量トルク算出手段と、
前記クランク角速度と、前記機械損失トルクと、前記慣性質量トルクとに基づいて、筒内圧により生ずる筒内圧トルクを算出する筒内圧トルク算出手段と、
筒内容積変化率を算出する筒内容積変化率算出手段と、
前記筒内圧トルクと前記筒内容積変化率とから筒内圧を算出する筒内圧算出手段と、
を含むことを特徴とする。
【００１９】
また、第１１の発明は、第８乃至第１０の発明の何れかにおいて、
前記内燃機関は、クランク角を検出するクランク角センサを備え、
前記筒内圧取得手段は、前記クランク角センサの計測刻み毎の筒内圧を取得することを特徴とする。
【００２０】
また、第１２の発明は、第８乃至第１１の発明の何れかにおいて、
前記内燃機関は、複数気筒を有し、
前記筒内圧取得手段は、
クランク角速度を取得するクランク角速度取得手段と、
前記内燃機関の機械損失により生ずる機械損失トルクを算出する機械損失トルク算出手段と、
前記内燃機関の運動部材により生ずる慣性質量トルクを算出する慣性質量トルク算出手段と、
前記クランク角速度と、前記機械損失トルクと、前記慣性質量トルクとに基づいて、全気筒の筒内圧により生ずる全気筒筒内圧トルクを算出する全気筒筒内圧トルク算出手段と、
各気筒の筒内容積変化率を算出する筒内容積変化率算出手段と、
前記全気筒筒内圧トルクと、前記各気筒の筒内容積変化率とから、気筒別筒内圧を算出する気筒別筒内圧算出手段と、
を含み、
前記パラメータ決定手段は、前記気筒別筒内圧に基づいて、前記未定パラメータを気筒別に決定し、
前記筒内圧推定手段は、気筒別に決定された前記パラメータが代入された気筒別の前記モデル式に基づいて、所定クランク角での筒内圧推定値を気筒別に算出することを特徴とする。
【００２１】
また、第１３の発明は、第１２の発明において、
前記気筒別筒内圧算出手段は、気筒別筒内圧が所定時間変化しないとの仮定に基づいて、前記気筒別筒内圧を算出することを特徴とする。
【００２２】
また、第１４の発明は、第１２の発明において、
前記内燃機関は、クランク角を検出するクランク角センサを備え、
前記クランク角速度取得手段は、前記クランク角センサの計測刻み毎に取得されるクランク角速度を補間することにより、前記計測刻みより小さい刻み幅の補間後クランク角速度を算出する補間後クランク角速度算出手段を含み、
前記気筒別筒内圧算出手段は、気筒別筒内圧が、前記小さい刻み幅に応じた所定の時間変化しないとの仮定に基づいて、前記気筒別筒内圧を算出することを特徴とする。
【発明の効果】
【００２３】
第１の発明によれば、内燃機関の吸気弁閉時期と排気弁開時期との間の複数点の筒内圧を取得し、その取得された筒内圧に基づいて、筒内の熱発生パターンと筒内圧との関係を表すモデル式に含まれる複数の未定パラメータを決定することができる。そして、その決定したパラメータを代入したモデル式に基づいて、所定クランク角での筒内圧推定値を算出することができる。このため、第１の発明によれば、クランク角センサ信号のノイズ等の影響を、上記モデル式を用いることで補正することができるので、高精度な筒内圧推定値を得ることができる。また、クランク角センサの計測刻みにかかわらず、より細かいクランク角刻みで筒内圧推定値を求めることができる。
【００２４】
第２の発明によれば、形状パラメータm、効率パラメータk、燃焼期間θpおよび熱発生開始クランク角θbを含むWiebe関数に基づいたモデル式を利用することにより、推定精度を更に向上することができる。
【００２５】
第３の発明によれば、クランク角速度と、機械損失トルクと、慣性質量トルクとに基づいて、筒内圧により生ずる筒内圧トルクを算出し、その筒内圧トルクと筒内容積変化率とから、上記モデル式のパラメータを決定するための筒内圧を精度良く算出することができる。
【００２６】
第４の発明によれば、上記モデル式のパラメータを決定するための筒内圧をクランク角センサの計測刻み毎に取得することができる。このため、上記モデル式のパラメータをより高い精度で決定することができる。
【００２７】
第５の発明によれば、気筒別筒内圧を算出し、その気筒別筒内圧に基づいて、上記モデル式のパラメータを気筒別に決定することができる。このため、気筒別のモデル式に基づいて、筒内圧推定値を気筒別に精度良く算出することができる。
【００２８】
第６の発明によれば、気筒別筒内圧が所定時間変化しないとの仮定に基づいて気筒別筒内圧を算出する。このため、簡便な方法で、気筒別筒内圧を精度良く算出することができる。
【００２９】
第７の発明によれば、クランク角センサの計測刻み毎に取得されるクランク角速度を補間することにより、その計測刻みより小さい刻み幅の補間後クランク角速度を算出することができる。そして、気筒別筒内圧が、その小さい刻み幅に応じた所定の時間変化しないとの仮定に基づいて気筒別筒内圧を算出する。このため、第７の発明によれば、クランク角センサの計測刻みが大きい場合であっても、上記仮定を精度良く成立させることができる。よって、気筒別筒内圧をより高い精度で算出することができる。
【００３０】
第８の発明によれば、内燃機関の吸気弁閉時期と排気弁開時期との間の複数点の筒内圧を取得し、その取得された筒内圧に基づいて、筒内の熱発生パターンと筒内圧との関係を表すモデル式に含まれる複数の未定パラメータを決定することができる。そして、その決定したパラメータを代入したモデル式に基づいて、所定クランク角での筒内圧推定値を算出することができる。このため、第８の発明によれば、クランク角センサ信号のノイズ等の影響を、上記モデル式を用いることで補正することができるので、高精度な筒内圧推定値を得ることができる。また、クランク角センサの計測刻みにかかわらず、より細かいクランク角刻みで筒内圧推定値を求めることができる。
【００３１】
第９の発明によれば、形状パラメータm、効率パラメータk、燃焼期間θpおよび熱発生開始クランク角θbを含むWiebe関数に基づいたモデル式を利用することにより、推定精度を更に向上することができる。
【００３２】
第１０の発明によれば、クランク角速度と、機械損失トルクと、慣性質量トルクとに基づいて、筒内圧により生ずる筒内圧トルクを算出し、その筒内圧トルクと筒内容積変化率とから、上記モデル式のパラメータを決定するための筒内圧を精度良く算出することができる。
【００３３】
第１１の発明によれば、上記モデル式のパラメータを決定するための筒内圧をクランク角センサの計測刻み毎に取得することができる。このため、上記モデル式のパラメータをより高い精度で決定することができる。
【００３４】
第１２の発明によれば、気筒別筒内圧を算出し、その気筒別筒内圧に基づいて、上記モデル式のパラメータを気筒別に決定することができる。このため、気筒別のモデル式に基づいて、筒内圧推定値を気筒別に精度良く算出することができる。
【００３５】
第１３の発明によれば、気筒別筒内圧が所定時間変化しないとの仮定に基づいて気筒別筒内圧を算出する。このため、簡便な方法で、気筒別筒内圧を精度良く算出することができる。
【００３６】
第１４の発明によれば、クランク角センサの計測刻み毎に取得されるクランク角速度を補間することにより、その計測刻みより小さい刻み幅の補間後クランク角速度を算出することができる。そして、気筒別筒内圧が、その小さい刻み幅に応じた所定の時間変化しないとの仮定に基づいて気筒別筒内圧を算出する。このため、第１４の発明によれば、クランク角センサの計測刻みが大きい場合であっても、上記仮定を精度良く成立させることができる。よって、気筒別筒内圧をより高い精度で算出することができる。
【発明を実施するための最良の形態】
【００３７】
実施の形態１．
［システム構成の説明］
図１は、本発明の実施の形態１のシステム構成を説明するための図である。図１に示すように、本実施形態のシステムは、内燃機関１０を備えている。
【００３８】
内燃機関１０の筒内には、吸気通路１２および排気通路１４が連通している。吸気通路１２には、吸入空気量Gaを検出するエアフローメータ１６が配置されている。エアフローメータ１６の下流には、スロットル弁１８が配置されている。スロットル弁１８の開度は、スロットルモータ２０の作動によって調整される。スロットル弁１８の近傍には、スロットル弁１８の開度を検出するためのスロットルポジションセンサ２２が配置されている。また、アクセルペダルの近傍には、アクセルペダルの踏み込み量を検出するアクセルポジションセンサ２４が設けられている。
【００３９】
内燃機関１０の気筒には、吸気ポート１１内に燃料を噴射するための燃料インジェクタ２６が配置されている。内燃機関１０の気筒には、更に、吸気弁２８、点火プラグ３０、排気弁３２およびピストン３４が設けられている。なお、本発明における内燃機関は、図示のようなポート噴射式のものに限らず、燃料を筒内に直接噴射する筒内直接噴射式のものでもよく、更にはポート噴射と筒内噴射とを併用するものでもよい。
【００４０】
また、本実施形態のシステムは、クランク軸３６の回転角度（回転位置）、すなわちクランク角θを検出するクランク角センサ３８と、エンジンオイルの温度を検出する油温センサ４２と、ＥＣＵ(Electronic Control Unit)５０とを備えている。ＥＣＵ５０には、上述した各種のセンサおよびアクチュエータが電気的に接続されている。
【００４１】
［実施の形態１の特徴］
本実施形態では、クランク角センサ３８の出力に加えて、熱発生パターンモデルであるWiebe関数を利用したモデル式を用いることにより、高精度な筒内圧推定値を算出する。以下、その具体的手順について説明する。図２および図３は、筒内圧推定値を算出するために本実施形態においてＥＣＵ５０が実行するルーチンのフローチャートである。
【００４２】
これらのルーチンによれば、まず、クランク角速度ωの履歴が検出される。具体的には、クランク角センサ３８の出力に基づいて、クランク角θの履歴が時間tの関数θ(t)として取得される（図２のステップ１００，１０２）。次いで、θ(t)を時間tで微分することにより、クランク角速度履歴ω(t)が算出され（ステップ１０４）、そのω(t)をクランク角θの関数に変換することにより、クランク角速度履歴ω(θ)が算出される（ステップ１０６）。
【００４３】
これとともに、負荷率KLの履歴が検出される。具体的には、エアフローメータ１６で検出される吸入空気量Gaに基づいて負荷率履歴KL(t)が算出され（ステップ１０８，１１０）、このKL(t)をθの関数に変換することにより、負荷率履歴KL(θ)が算出される（ステップ１１２）。
【００４４】
更に、エンジンオイル温度Toの履歴が検出される。具体的には、油温センサ４２の出力に基づいて、エンジンオイル温度履歴To(t)が算出され（ステップ１１４，１１６）、このTo(t)をθの関数に変換することにより、エンジンオイル温度履歴To(θ)が算出される（ステップ１１８）。
【００４５】
続いて、以上算出されたクランク角速度履歴ω(θ)、負荷率履歴KL(θ)およびエンジンオイル温度履歴To(θ)を用いて、内燃機関１０の機械損失（ピストン３４、クランク軸３６等のフリクション）によって生ずる機械損失トルクTfの履歴が算出される（ステップ１２０）。機械損失トルクTfは、クランク角速度ωが高いほど大きく、また、負荷率KLが高いほど大きく、更に、エンジンオイル温度Toが低いほど大きい。ＥＣＵ５０には、これらの傾向を表したマップが予め記憶されており、ステップ１２０では、そのマップに基づいて、機械損失トルク履歴Tf(θ)が算出される。
【００４６】
また、図２に示すルーチンによれば、更に、内燃機関１０の慣性モーメントJの履歴と、内燃機関１０の運動部材（ピストン３４等）の慣性により生ずる慣性質量トルクTmの履歴とが算出される（ステップ１２２，１２４）。慣性モーメントJおよび慣性質量トルクTmは、設計値より、クランク角θの関数として求めることができる。そこで、ステップ１２２，１２４においては、ＥＣＵ５０に予め記憶されたマップあるいは数式に基づいて、慣性モーメント履歴J(θ)および慣性質量トルク履歴Tm(θ)が算出される。
【００４７】
ここで、内燃機関１０の筒内圧Pcがピストン３４に及ぼす力によって生ずるトルク（以下「筒内圧トルク」という）をTpとすると、クランク軸３６の運動方程式は、下記（１）式で表される。よって、筒内圧トルクTpは、下記（２）式で表される。
【００４８】
【数１】
2008255932.tif 000003
【００４９】
図２に示すルーチンによれば、上述の処理によって算出されたクランク角速度履歴ω(θ)、機械損失トルク履歴Tf(θ)、慣性モーメント履歴J(θ)および慣性質量トルク履歴Tm(θ)を、上記（２）式に代入することにより、筒内圧トルク履歴Tp(θ)が算出される（ステップ１２６）。
【００５０】
続いて、クランク角θに対する筒内容積Vの変化率dV(θ)/dθの履歴が算出される（ステップ１２８）。筒内容積変化率dV(θ)/dθは、クランク角θに応じて幾何学的に決まるので、設計値から求めることができる。そこで、ステップ１２８においては、ＥＣＵ５０に予め記憶されたマップあるいは数式に基づいて、筒内容積変化率履歴dV(θ)/dθが算出される。
【００５１】
続いて、筒内圧履歴Pc(θ)が算出される（ステップ１３０）。筒内圧Pc、筒内圧トルクTpおよび筒内容積変化率dV(θ)/dθの間には、力学上、下記（３）式が成立する。そこで、このステップ１３０では、下記（３）式に、上記ステップ１２６で求められたTp(θ)と、上記ステップ１２８で求められたdV(θ)/dθとを代入することにより、筒内圧履歴Pc(θ)が算出される。
【００５２】
【数２】
2008255932.tif 000004
【００５３】
ここで、本実施形態において用いる、Wiebe関数を利用したモデル式について説明する。Wiebe関数によれば、クランク角θにおける質量燃焼割合Xbは、下記（４）式で表すことができる。
【００５４】
【数３】
2008255932.tif 000005
【００５５】
ただし、θbは熱発生が開始するクランク角、θpは燃焼期間、mは形状パラメータ、aは所定の定数である。
【００５６】
また、燃料の燃焼により発生する熱量をQとすると、熱発生率dQ/dθは、上記（４）式を微分した式を用いることにより、次式（５）で表すことができる。
【００５７】
【数４】
2008255932.tif 000006
【００５８】
ただし、Qfは筒内に供給された燃料が有する熱量、kは燃料の熱量Qfが実際に熱に変換される効率を表す効率パラメータである。
【００５９】
上記（５）式と、下記（６）式で表される熱力学第１法則とから、下記（７）式が得られる。
【００６０】
【数５】
2008255932.tif 000007
【００６１】
ただし、γは比熱比である。上記（７）式が本実施形態で用いるモデル式である。このモデル式において、燃料の熱量Qfは、空燃比あるいは燃料噴射量等に基づいて算出することができる。比熱比γは、既知の値である。また、本実施形態おいて、aは、前述したように所定の定数（例えば６．９）とする。そして、dθ/dt（＝ω）、筒内容積変化率dV(θ)/dθ、筒内容積V(θ)は、前述した処理で算出されたクランク角速度ωの履歴および筒内容積変化率dV(θ)/dθの履歴から求めることができる。よって、上記モデル式中の未定パラメータは、効率パラメータk、形状パラメータm、熱発生開始クランク角θbおよび燃焼期間θpの４つとなる。
【００６２】
上記ステップ１３０で筒内圧履歴Pc(θ)が算出された後は、上記未定パラメータk、m、θbおよびθpの値を決定するべく、最後の（直近の）サイクルの吸気弁閉時期（以下「IVC」という）から排気弁開時期（以下「EVO」という）までの間の複数点の筒内圧Pc(θ)が抽出される（図３のステップ１３２）。
【００６３】
クランク角センサ３８は、所定の計測刻みΔθ（例えば１０度または３０度）毎に、クランク角θを検出する。ステップ１３２では、IVCからEVOまでの計測刻みΔθ毎の筒内圧Pc(θ)を、上記ステップ１３０で算出された筒内圧履歴Pc(θ)から抽出する。よって、ここで抽出されるデータ数nは、IVCからEVOまでのクランク角期間を計測刻みΔθで割ったときの商に等しい。本実施形態では、このデータ数nは、上記モデル式の未定パラメータの数以上、つまり４以上であればよい。
【００６４】
上記ステップ１３２によってn個の筒内圧Pc(θ)のデータが抽出されたら、その各データを上記（７）式に代入することにより、k、m、θbおよびθpの４つの未定パラメータを変数として含むn個（n≧４）の式が生成される（ステップ１３４）。これらn個の式を、連立方程式として、あるいは最小二乗法を用いて解くことにより、k、m、θbおよびθpがそれぞれ算出される（ステップ１３６）。
【００６５】
上記ステップ１３６で算出されたk、m、θbおよびθpは、下記（８）式に代入される（ステップ１３８）。
【００６６】
【数６】
2008255932.tif 000008
【００６７】
上記（８）式は、上記（７）式中の筒内圧Pc(θ)を補正後筒内圧Pw(θ)で置き換えた式である。この（８）式をPw(θ)の微分方程式として解くことにより、補正後筒内圧履歴Pw(θ)が算出される（ステップ１４０）。
【００６８】
このようにして算出された補正後筒内圧履歴Pw(θ)は、高精度な筒内圧推定値として用いることができる。すなわち、上記ステップ１３０で算出された筒内圧履歴Pc(θ)には、クランク角センサ３８のノイズに起因する誤差が含まれているが、Wiebe関数を利用した上記モデル式を用いて算出された補正後筒内圧履歴Pw(θ)においては、クランク角センサ３８のノイズによる誤差が精度良く補正されている。また、補正後筒内圧履歴Pw(θ)によれば、所望の（任意の）クランク角θにおける補正後筒内圧Pw(θ)を精度良く算出することができる。このため、クランク角センサ３８の計測刻みΔθの大きさにかかわらず、より細かいクランク角刻みで補正後筒内圧Pw(θ)を精度良く算出することができる。
【００６９】
よって、本実施形態では、補正後筒内圧Pw(θ)を筒内圧推定値として用いることにより、例えば燃料性状検出、失火検出、筒内状態量（燃焼割合等）をベースとする制御など、筒内圧を利用した制御を精度良く実行することができる。
【００７０】
また、上述した実施の形態１においては、ＥＣＵ５０が、上記ステップ１００〜１３２の処理を実行することにより前記第８の発明における「筒内圧取得手段」が、上記ステップ１３４および１３６の処理を実行することにより前記第８の発明における「パラメータ決定手段」が、上記ステップ１３８および１４０の処理を実行することにより前記第８の発明における「筒内圧推定手段」が、上記ステップ１００〜１０６の処理を実行することにより前記第１０の発明における「クランク角速度取得手段」が、上記ステップ１０８〜１２０の処理を実行することにより前記第１０の発明における「機械損失トルク算出手段」が、上記ステップ１２４の処理を実行することにより前記第１０の発明における「慣性質量トルク算出手段」が、上記ステップ１２６の処理を実行することにより前記第１０の発明における「筒内圧トルク算出手段」が、上記ステップ１２８の処理を実行することにより前記第１０の発明における「筒内容積変化率算出手段」が、上記ステップ１３０の処理を実行することにより前記第１０の発明における「筒内圧算出手段」が、それぞれ実現されている。
【００７１】
実施の形態２．
次に、図４乃至図６を参照して、本発明の実施の形態２について説明するが、上述した実施の形態１との相違点を中心に説明し、同様の事項については、その説明を簡略化または省略する。
【００７２】
［実施の形態２の特徴］
複数の気筒を有する内燃機関１０においては、上記（２）式で算出される筒内圧トルクTpは、それら複数気筒の筒内圧トルクが足し合わせられたものである。このため、上記（３）式では、筒内圧Pcを気筒別に算出することはできない。本実施形態は、気筒別の筒内圧を算出できるように、実施の形態１に改良を施したものである。なお、以下の説明において、内燃機関１０は、＃１〜＃４の４つの気筒を有しているものとする。また、以下の説明において、各記号の添え字のｉは、気筒番号を表す。すなわち、（ｉ＝１，２，３，４）である。
【００７３】
図４乃至図６は、気筒別の筒内圧推定値を算出するために本実施形態においてＥＣＵ５０が実行するルーチンのフローチャートである。
【００７４】
これらのルーチンによれば、まず、筒内圧トルク履歴Tp(θ)が算出される（図４のステップ１４２）。このステップ１４２の処理は、前述した実施の形態１のステップ１００〜１２６の処理と同様である。
【００７５】
続いて、＃１気筒の筒内容積変化率履歴dV1(θ)/dθと、＃２気筒の筒内容積変化率履歴dV2(θ)/dθと、＃３気筒の筒内容積変化率履歴dV3(θ)/dθと、＃４気筒の筒内容積変化率履歴dV4(θ)/dθとがそれぞれ算出される（ステップ１４４，１４６，１４８，１５０）。この場合、各気筒の筒内容積変化率dV(θ)/dθは、ＥＣＵ５０に予め記憶された気筒毎のマップあるいは数式に基づいて算出することができる。
【００７６】
ここで、＃１〜＃４の各気筒の筒内圧をそれぞれPc1，Pc2，Pc3，Pc4とし、現在のクランク角をθ0とすると、現在の筒内圧トルクTp(θ0)は、下記（９）式で表すことができる。
【００７７】
【数７】
2008255932.tif 000009
【００７８】
同様にして、クランク角センサ３８の１計測刻み前の筒内圧トルクTp(θ0−Δθ)、２計測刻み前の筒内圧トルクTp(θ0−２Δθ)、３計測刻み前の筒内圧トルクTp(θ0−３Δθ)は、それぞれ、下記（１０）〜（１２）式で表すことができる。
【００７９】
【数８】
2008255932.tif 000010
【００８０】
なお、上記（９）〜（１２）式中のdVi(θ0)/dθ，dVi(θ0−Δθ)/dθ，dVi(θ0−２Δθ)/dθ，dVi(θ0−３Δθ)/dθ（ｉ＝１，２，３，４）は、それぞれ、上記ステップ１４４〜１５０で算出された筒内容積変化率dVi(θ)/dθから抽出することができる。
【００８１】
ここで、本実施形態では、一定時間、具体的には(θ0−３Δθ)≦θ≦θ0の間、各気筒の筒内圧Pciが変化しないと仮定する。すなわち、下記（１３）式を仮定する。
【００８２】
【数９】
2008255932.tif 000011
【００８３】
本実施形態では、上記（１３）式の仮定の下で、各気筒の筒内圧Pci(θ0)を算出する。すなわち、まず、上記（１３）式を上記（９）〜（１２）式に代入することにより、Pci(θ0−Δθ)，Pci(θ0−２Δθ)，Pci(θ0−３Δθ)（ｉ＝１，２，３，４）を消去する（図５のステップ１５２）。
【００８４】
これにより、上記（９）〜（１２）式は、４つの未知数Pci(θ0)（ｉ＝１，２，３，４）を含む４つの式となるので、これらを連立方程式として解くことにより、Pci(θ0)（ｉ＝１，２，３，４）、すなわち、＃１〜＃４気筒の現在の気筒別筒内圧を算出することができる（ステップ１５４）。
【００８５】
続いて、現在クランク角θ0が下記（１４）式により更新され（ステップ１５６）、上記ステップ１５４の処理が再度実行される。すなわち、現在クランク角θ0が、クランク角センサ３８の１計測刻みΔθ分だけ進められた後、気筒別筒内圧Pci(θ0)（ｉ＝１，２，３，４）が再度算出される。
θ0＝θ0＋Δθ      ・・・（１４）
【００８６】
そして、上記ステップ１５２〜１５６のループが、内燃機関１０の１サイクル分、すなわちクランク角７２０度分の気筒別筒内圧Pci(θ0)（ｉ＝１，２，３，４）を取得するまで、繰り返し実行される。
【００８７】
続いて、正確な気筒別筒内圧推定値を求めたい気筒（ここでは、最後にEVOを迎えた気筒とする）を＃ｉ気筒としたとき、この＃ｉ気筒の最後のサイクルのIVCからEVOまでの間のni点の筒内圧Pci(θ)が、上記１５２〜１５６の処理により算出されたデータの中から抽出される（ステップ１５８）。このデータ数niは、実施の形態１と同様に、＃ｉ気筒の最後IVCからEVOまでのクランク角期間を計測刻みΔθで割ったときの商であり、モデル式の未定パラメータの数以上、つまり４以上であればよい。
【００８８】
続いて、上記ステップ１５８によって抽出されたni個のデータが下記（１５）式に代入される（図６のステップ１６０）。
【００８９】
【数１０】
2008255932.tif 000012
【００９０】
上記（１５）式は、本実施形態で用いるモデル式であり、実施の形態１における（７）式と同様である。上記ステップ１６０の処理により、ki、mi、θbiおよびθpiの４つの未定パラメータを変数として含むni個（ni≧４）の式が生成される。これらni個の式を、連立方程式として、あるいは最小二乗法を用いて解くことにより、＃ｉ気筒のモデル式におけるki、mi、θbiおよびθpiがそれぞれ算出される（ステップ１６２）。
【００９１】
上記ステップ１６２で算出されたki、mi、θbiおよびθpiは、下記（１６）式に代入される（ステップ１６４）。
【００９２】
【数１１】
2008255932.tif 000013
【００９３】
上記（１６）式は、上記（１５）式中の筒内圧Pci(θ)を補正後筒内圧Pwi(θ)で置き換えた式である。この（１６）式をPwi(θ)の微分方程式として解くことにより、＃ｉ気筒の補正後筒内圧履歴Pwi(θ)が算出される（ステップ１６６）。
【００９４】
以上説明したように、本実施形態によれば、複数気筒を有する内燃機関１０の気筒別の筒内圧推定値（補正後筒内圧履歴Pwi(θ)）を高精度に算出することができる。
【００９５】
また、上述した実施の形態２においては、ＥＣＵ５０が、上記ステップ１４２の処理を実行することにより前記第１２の発明における「全気筒筒内圧トルク算出手段」が、上記ステップ１４４〜１５０の処理を実行することにより前記第１２の発明における「筒内容積変化率算出手段」が、上記ステップ１５２〜１５６の処理を実行することにより前記第１２の発明における「気筒別筒内圧算出手段」が、それぞれ実現されている。
【００９６】
実施の形態３．
次に、図７乃至図１０を参照して、本発明の実施の形態３について説明するが、上述した実施の形態１との相違点を中心に説明し、同様の事項については、その説明を簡略化または省略する。
【００９７】
［実施の形態３の特徴］
前述した実施の形態２における上記（１３）式の仮定は、クランク角センサ３８の計測刻みΔθが大きい場合には、成立しにくくなる。このため、気筒別筒内圧推定値（Pwi）の推定精度を十分に確保できない場合もある。本実施形態では、そのような場合であっても十分な推定精度を得るべく、クランク角速度ω(θ)の波形を適切な関数を用いて補間することにより、計算上のクランク角刻みを小さくすることとした。
【００９８】
図７乃至図１０は、気筒別の筒内圧推定値を算出するために本実施形態においてＥＣＵ５０が実行するルーチンのフローチャートである。
【００９９】
これらのルーチンによれば、まず、クランク角速度履歴ω(θ)が算出される（図７のステップ１６８）。この処理は、前述した実施の形態１のステップ１００〜１０６の処理と同様である。
【０１００】
本実施形態では、クランク角速度ω(θ)の波形を、下記（１７）式の３次関数により補間することとした。
【０１０１】
【数１２】
2008255932.tif 000014
【０１０２】
上記ステップ１６８の処理に続いて、上記（１７）式中の４つの未定係数α3，α2，α1，α0を決定するため、現在のクランク角θ0、クランク角センサ３８の１計測刻み前のクランク角(θ0−Δθ)、２計測刻み前のクランク角(θ0−２Δθ)、３計測刻み前のクランク角(θ0−３Δθ)、の４つの時点におけるクランク角速度ωが、上記クランク角速履歴から抽出される（ステップ１７０）。すなわち、ω(θ0−３Δθ)，ω(θ0−２Δθ)，ω(θ0−Δθ)，ω(θ0)が抽出される。
【０１０３】
続いて、上記ステップ１７０で抽出されたω(θ0−３Δθ)，ω(θ0−２Δθ)，ω(θ0−Δθ)，ω(θ0)の値が上記（１７）式に代入される（ステップ１７２）。これにより、未知数４、式数４の連立方程式が得られるので、この連立方程式を解くことにより、係数α3，α2，α1，α0の値を算出することができる（ステップ１７４）。
【０１０４】
次いで、算出された係数α3，α2，α1，α0の値が、下記（１８）式に代入される（ステップ１７６）。
【０１０５】
【数１３】
2008255932.tif 000015
【０１０６】
上記（１８）式中のωcは、補間後クランク角速度である。この（１８）式によれば、クランク角センサ３８の計測刻みΔθよりも大幅に小さい計算刻みΔθc（例えば、Δθc＝０．１度）で、補間後クランク角速度ωcを精度良く算出することができる。
【０１０７】
本実施形態では、更に精度を向上するため、(θ0−３Δθ)、(θ0−２Δθ)、(θ0−Δθ)、θ0の４点を用いて求めた上記（１８）式の補間後クランク角速度ωcのうち、(θ0−２Δθ)から(θ0−Δθ)までを使用することとした。そこで、本実施形態では、上記ステップ１７６の処理に続いて、上記（１８）式により、(θ0−２Δθ)から(θ0−Δθ)まで、計算刻みΔθc毎に補間後クランク角速度ωcが算出される（ステップ１７８）。すなわち、θ＝(θ0−２Δθ，θ0−２Δθ＋Δθc，θ0−２Δθ＋２Δθc，・・・，θ0−Δθ)のときの補間後クランク角速度ωcが算出される。
【０１０８】
続いて、現在クランク角θ0が下記（１９）式により更新され（ステップ１８０）、上記ステップ１７０〜１８０の処理が再度実行される。すなわち、現在クランク角θ0が、クランク角センサ３８の１計測刻みΔθ分だけ進められた後、(θ0−２Δθ)から(θ0−Δθ)までの補間後クランク角速度ωcが計算刻みΔθc毎に算出される。
θ0＝θ0＋Δθ      ・・・（１９）
【０１０９】
そして、上記ステップ１７０〜１８０のループが、内燃機関１０の１サイクル分、すなわちクランク角７２０度分、繰り返し実行される。これにより、計算刻みΔθc毎の７２０度分の補間後クランク角速度ωcの履歴を取得することができる。
【０１１０】
上記の処理とともに、負荷率履歴KL(θ)が算出される（図８のステップ１８２）。この算出処理は、前述した実施の形態１のステップ１０８〜１１２の処理と同様である。また、エンジンオイル履歴To(θ)が算出される（ステップ１８４）。この算出処理は、前述した実施の形態１のステップ１４８〜１１８の処理と同様である。
【０１１１】
続いて、以上の処理によって算出された補間後クランク角速度ωcの履歴と、負荷率履歴KL(θ)と、エンジンオイル温度履歴To(θ)とに基づいて、機械損失トルクTfの履歴が計算刻みΔθc毎に算出される（ステップ１８６）。この算出処理は、前述した実施の形態１のステップ１２０の処理と同様である。
【０１１２】
続いて、下記（２０）式に基づいて、計算刻みΔθc毎の筒内圧トルクTp(θ)の履歴が算出される（ステップ１８８）。なお、下記（２０）式は、実施の形態１で説明した（２）式と同様の式である。
【０１１３】
【数１４】
2008255932.tif 000016
【０１１４】
＃１〜＃４の各気筒の筒内圧をそれぞれPc1，Pc2，Pc3，Pc4とすると、クランク角θ00における筒内圧トルクTp(θ00)、その１計算刻みΔθc前の筒内圧トルクTp(θ00−Δθc)、２計算刻みΔθc前の筒内圧トルクTp(θ00−２Δθc)、３計算刻みΔθc前の筒内圧トルクTp(θ00−３Δθc)は、それぞれ、下記（２１）〜（２４）式で表すことができる。  
【０１１５】
【数１５】
2008255932.tif 000017
【０１１６】
ここで、本実施形態では、一定時間、具体的には(θ00−３Δθc)≦θ≦θ00の間、各気筒の筒内圧Pciが変化しないと仮定する。すなわち、下記（２５）式を仮定する。
【０１１７】
【数１６】
2008255932.tif 000018
【０１１８】
本実施形態では、上述したように、計算刻みΔθcを極めて小さくすることができる。このため、上記（２５）式の仮定は、クランク角センサ３８の計測刻みΔθの大きさにかかわらず、精度良く成立する。
【０１１９】
上記（２５）式の仮定の下で各気筒の筒内圧Pci(θ0)を算出するには、上記（２５）式を上記（２１）〜（２４）式に代入することにより、Pci(θ00−Δθc)，Pci(θ00−２Δθc)，Pci(θ00−３Δθc)（ｉ＝１，２，３，４）を消去する（図９のステップ１９０）。
【０１２０】
これにより、上記（２１）〜（２４）式は、４つの未知数Pci(θ00)（ｉ＝１，２，３，４）を含む４つの式となるので、これらを連立方程式として解くことにより、クランク角θ00における＃１〜＃４気筒の気筒別筒内圧Pci(θ00)（ｉ＝１，２，３，４）を算出することができる（ステップ１９２）。
【０１２１】
続いて、クランク角θ00が下記（２６）式により更新され（ステップ１９４）、上記ステップ１９０の処理が再度実行される。すなわち、クランク角θ00が１計算刻みΔθc分だけ進められた後、各気筒の気筒別筒内圧Pci(θ00)（ｉ＝１，２，３，４）が再度算出される。
θ00＝θ00＋Δθc      ・・・（２６）
【０１２２】
そして、上記ステップ１９０〜１９４のループが、内燃機関１０の１サイクル、すなわちクランク角７２０度分の気筒別筒内圧Pci(θ00)（ｉ＝１，２，３，４）を取得するまで、繰り返し実行される。
【０１２３】
続いて、正確な気筒別筒内圧推定値を求めたい気筒（ここでは、最後にEVOを迎えた気筒とする）を＃ｉ気筒としたとき、この＃ｉ気筒の最後のサイクルのIVCからEVOまでの間のni点の筒内圧Pci(θ)が、上記１９０〜１９４の処理により算出されたデータの中から抽出される（ステップ１９６）。このデータ数niは、＃ｉ気筒の最後IVCからEVOまでのクランク角期間を計算刻みΔθcで割ったときの商である。
【０１２４】
続いて、上記ステップ１９６によって抽出されたni個のデータが、前述した実施の形態２の（１５）式と同様のモデル式に代入される（図１０のステップ１９８）。
【０１２５】
上記ステップ１９８の処理により、ki、mi、θbiおよびθpiの４つの未定パラメータを変数として含むni個（ni≧４）の式が生成される。これらni個の式を、連立方程式として、あるいは最小二乗法を用いて解くことにより、＃ｉ気筒のモデル式におけるki、mi、θbiおよびθpiがそれぞれ算出される（ステップ２００）。
【０１２６】
上記ステップ２００で算出されたki、mi、θbiおよびθpiは、前述した実施の形態２の（１６）式と同様の式に代入される（ステップ２０２）。この式をPwi(θ)の微分方程式として解くことにより、＃ｉ気筒の補正後筒内圧履歴Pwi(θ)が算出される（ステップ２０４）。
【０１２７】
以上説明したように、本実施形態では、計算刻みΔθcを極めて小さくすることができるので、上記（２５）式の仮定が精度良く成立する。このため、複数気筒を有する内燃機関１０の気筒別の筒内圧推定値（補正後筒内圧履歴Pwi(θ)）をより高い精度で算出することができる。
【０１２８】
また、上述した実施の形態３においては、ＥＣＵ５０が、上記ステップ１６８の〜１８０処理を実行することにより前記第１４の発明における「クランク角速度取得手段」が、上記ステップ１７０〜１８０の処理を実行することにより前記第１４の発明における「補間後クランク角速度算出手段」が、上記ステップ１８２の〜１８８処理を実行することにより前記第１４の発明における「全気筒筒内圧トルク算出手段」が、上記ステップ１９０〜１９４の処理を実行することにより前記第１４の発明における「気筒別筒内圧算出手段」が、それぞれ実現されている。
【０１２９】
なお、上述した実施の形態３においては、クランク角速度ω(θ)の波形を３次関数を用いて補間しているが、この補間に用いる関数はこれに限定されるものではなく、例えばスプライン関数などを用いてもよい。
【０１３０】
また、上述した各実施の形態では、上記（１）式のようなクランク軸３６の運動方程式を利用して求めた筒内圧トルクTp(θ)から筒内圧履歴Pc(θ)を求めることとしているが、本発明では、筒内圧センサを設け、その筒内圧センサの出力から筒内圧履歴Pc(θ)を求めるようにしてもよい。
【図面の簡単な説明】
【０１３１】
【図１】本発明の実施の形態１のシステム構成を説明するための図である。
【図２】本発明の実施の形態１において実行されるルーチンのフローチャートである。
【図３】本発明の実施の形態１において実行されるルーチンのフローチャートである。
【図４】本発明の実施の形態２において実行されるルーチンのフローチャートである。
【図５】本発明の実施の形態２において実行されるルーチンのフローチャートである。
【図６】本発明の実施の形態２において実行されるルーチンのフローチャートである。
【図７】本発明の実施の形態３において実行されるルーチンのフローチャートである。
【図８】本発明の実施の形態３において実行されるルーチンのフローチャートである。
【図９】本発明の実施の形態３において実行されるルーチンのフローチャートである。
【図１０】本発明の実施の形態３において実行されるルーチンのフローチャートである。
【符号の説明】
【０１３２】
１０  内燃機関
１１  吸気ポート
１２  吸気通路
１４  排気通路
１６  エアフローメータ
１８  スロットル弁
２６  燃料インジェクタ
３０  点火プラグ
３６  クランク軸
３８  クランク角センサ
４２  水温センサ
５０  ＥＣＵ
【図１】
2008255932.tif 000019
【図２】
2008255932.tif 000020
【図３】
2008255932.tif 000021
【図４】
2008255932.tif 000022
【図５】
2008255932.tif 000023
【図６】
2008255932.tif 000024
【図７】
2008255932.tif 000025
【図８】
2008255932.tif 000026
【図９】
2008255932.tif 000027
【図１０】
2008255932.tif 000028