(57)【要約】【課題】燃焼トルクの瞬時値を精度良く、算出することができるトルク推定装置を提供する。【解決手段】クランク軸の回転速度から、燃焼トルクを推定する内燃機関のトルク推定装置であって、前記クランク軸の回転速度を算出する回転速度算出手段１０１と、前記回転速度の波形をフィルタリングすることにより、前記回転速度の波形のうち所定の周波数成分の波形から燃焼トルクを算出する燃焼トルク算出手段１０２と、前記内燃機関の上死点及び下死点を判定する上下死点判定手段１０３と、前記上死点又は下死点であると判定したタイミングにおける前記燃焼トルクに基づいて、前記燃焼トルクの補正量を算出する燃焼トルク補正量算出手段１０７と、前記燃焼トルク補正量に基づいて、前記燃焼トルクを補正する燃焼トルク補正手段１０９と、を備える。【選択図】図６
2010190174.tif 000002
【特許請求の範囲】
【請求項１】
クランク軸の回転速度から、燃焼トルクを推定する内燃機関のトルク推定装置であって、
前記クランク軸の回転速度を算出する回転速度算出手段と、
前記回転速度の波形をフィルタリングすることにより、前記回転速度の波形のうち所定の周波数成分の波形から燃焼トルクを算出する燃焼トルク算出手段と、
前記内燃機関の上死点及び下死点を判定する上下死点判定手段と、
前記上死点又は下死点であると判定したタイミングにおける前記燃焼トルクに基づいて、前記燃焼トルクの補正量を算出する燃焼トルク補正量算出手段と、
前記燃焼トルク補正量に基づいて、前記燃焼トルクを補正する燃焼トルク補正手段と、を備えることを特徴とする内燃機関のトルク推定装置。
【請求項２】
燃焼トルク補正量算出手段は、前記上下死点判定手段が前記上死点又は下死点であると判定した第一のタイミングにおける燃焼トルクと、前記第一のタイミング後において、前記上下死点判定手段が前記上死点又は下死点であると判定した第二のタイミングにおける燃焼トルクと、に基づいて、第一のタイミングから第二のタイミングまでにおける前記燃焼トルクの補正量を算出することを特徴とする請求項１に記載の内燃機関のトルク推定装置。
【請求項３】
前記燃焼トルク算出手段は、
前記回転速度の前記波形を、前記フィルタリングにより、前記内燃機関の燃焼に同期した周波数の成分の波形と、該周波数の整数倍の周波数の成分の波形と、に分解する手段と、
各周波数成分ごとの波形に対して、該波形の位相と振幅を補正する手段と、
各周波数成分ごとの波形を重ね合わせて、前記燃焼トルクの波形を算出する手段と、
を備えることを特徴とする請求項１又は２に記載の内燃機関のトルク推定装置。
【請求項４】
前記トルク推定装置は、前記クランク軸に取付けられたクランク角センサの出力信号から、前記回転速度を算出する回転速度算出手段を備えることを特徴とする請求項１又は２に記載の内燃機関のトルク推定装置。
【請求項５】
前記回転速度センサは、直流発電機であることを特徴とする請求項４に記載の内燃機関のトルク推定装置。
【請求項６】
請求項１から５のいずれかに記載の燃焼トルク推定装置を備え、前記燃焼トルクから各気筒の燃焼エネルギーを推定する燃焼エネルギー推定装置であって、
該燃焼エネルギー推定装置は、前記クランク軸の角度を算出するクランク角度算出手段と、
前記クランク角度に応じて変化する前記燃焼トルクを積分する積分手段と、
該積分した値に基づいて各気筒の燃焼エネルギーを算出する燃焼エネルギー算出手段と、を備えることを特徴とする燃焼エネルギー推定装置。
【請求項７】
前記積分手段は、前記気筒ごとに前記クランク角度の積分区間を設定し、該積分区間に基づいて、前記燃焼トルクを積分し、前記各気筒の燃焼エネルギーを算出することを特徴とする請求項６に記載の燃焼エネルギー推定装置。
【請求項８】
前記気筒において設定されたクランク角度の積分区間は、前記気筒の筒内圧が他のすべての気筒の筒内圧に比べて高い筒内圧となる区間であることを特徴とする請求項７に記載の燃焼エネルギー推定装置。
【請求項９】
前記クランク角度の積分区間は、圧縮行程の時点から爆発行程の時点までの区間であることを特徴とする請求項７に記載の燃焼エネルギー推定装置。
【発明の詳細な説明】
【技術分野】
【０００１】
内燃機関のクランク軸の回転に基づいて気筒毎の燃焼エネルギーを推定する装置に関する。
【背景技術】
【０００２】
内燃機関では所定の回転速度の周期で燃焼が行われ、その燃焼に伴いクランク軸にトルクが動力として付与される。このトルクの推定は、例えば、内燃機関に筒内圧センサを配置し、筒内圧センサの検出値である筒内圧から、燃焼により発生した熱量を推定することにより、行われる場合がある。このトルクの推定値に基づいて、内燃機関の燃焼に関するフィードバック制御が行われている。
【０００３】
たとえば、内燃機関の燃焼トルクを推定する方法として以下の推定方法が提案されている。この推定方法は、クランク軸の回転角速度を検出し、該角速度の変化からトルク変動を推定するようにした内燃機関のトルク変動推定方法である。
【０００４】
具体的には、内燃機関の運転状態が、クランク軸の固有振動による角速度の変化波形に比べて爆発１次の角速度の変化波形の方がトルクの変動を再現性よく顕著に表している第１の運転状態であるか、或いは爆発１次の角速度の変化波形に比べてクランク軸の固有振動による角速度の変化波形の方がトルクの変動を再現性よく顕著に表している第２の運転状態であるかを判断し、第１の運転状態であると判断されたときには、爆発行程開始時における角速度の最小値と爆発行程中における角速度の最大値との差からトルク変動を推定するものである（特許文献１参照）。
【０００５】
この特許文献１に記載のトルク推定方法は、各速度の最小値と最大値の差からトルク変動を推定するものであり、このトルク変動から、失火検出のような気筒で急激にトルクが落ち込んだことを、前後の気筒と比較することで検出を可能とするものである。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２００７−３２４３３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかしながら、各速度の最小値と最大値を利用してトルクを推定する場合には、上述したように失火のような急激なトルク変動を検出するアプリケーションに対しては有効であるが、例えば、安定した燃焼状態においては、急激なトルク変動はないので、この方法ではトルクの瞬時値を知ることは容易ではなく、精度も高くない。そして、燃焼により発生した熱量を筒内圧から推定して燃焼の制御にフィードバックするような、従来筒内圧センサが使われてきた分野では、トルクの絶対値が瞬時値で知ることが望ましいが、この方法で推定されたトルクでは、燃焼制御を行うには充分であるとは言えない。
【０００８】
そこで、本願の発明は、筒内圧センサなどを用いることなく、内燃機関のクランク軸の回転に基づいて、容易にかつ精度良く燃焼トルクの瞬時値を推定することができる内燃機関のトルク推定装置を提供することにある。
【課題を解決するための手段】
【０００９】
前記課題を解決すべく、本発明に係るトルク推定装置は、クランク軸の回転速度から、燃焼トルクを推定する内燃機関のトルク推定装置であって、前記クランク軸の回転速度を算出する回転速度算出手段と、前記回転速度の波形をフィルタリングすることにより、前記回転速度の波形のうち所定の周波数成分の波形から燃焼トルクを算出する燃焼トルク算出手段と、前記内燃機関の上死点及び下死点を判定する上下死点判定手段と、前記上死点又は下死点であると判定したタイミングにおける前記燃焼トルクに基づいて、前記燃焼トルクの補正量を算出する燃焼トルク補正量算出手段と、前記燃焼トルク補正量に基づいて、前記燃焼トルクを補正する燃焼トルク補正手段と、を備えることを特徴とする。
【００１０】
本発明によれば、クランク軸の回転速度の波形から、燃焼トルクを測定するので、燃焼トルクの瞬時値を推定することができるばかりでなく、現実上、上死点及び下死点では、クランク軸に作用する燃焼トルクは０である点に着眼して、上死点及び下死点で算出される燃焼トルクを誤差とみなすことにより、燃焼トルクを補正することにより、より精度よく燃焼トルクを推定することができる。
【００１１】
本発明に係る内燃機関のトルク推定装置は、燃焼トルク補正量算出手段は、前記上下死点判定手段が前記上死点又は下死点であると判定した第一のタイミングにおける燃焼トルクと、前記第一のタイミング後において、前記上下死点判定手段が前記上死点又は下死点であると判定した第二のタイミングにおける燃焼トルクと、に基づいて、第一のタイミングから第二のタイミングまでにおける前記燃焼トルクの補正量を算出することがより好ましい。
【００１２】
本発明によれば、第一のタイミングと第二のタイミング間における燃焼トルクを用いて、その間の燃焼トルクの補正量を補間して算出することができる。なお、この際には、第一のタイミング及び第二のタイミングは、クランク角度を時間変数として、その間の補正量を、第一のタイミングと第二のタイミング間における燃焼トルクを用いて線形補完してもよい。また、第二のタイミングは、第一のタイミング後の上死点又は下死点タイミングであればよいが、より好ましくは、次に判定される上死点又は下死点のタイミングである。これにより、より精度良く、燃焼トルクの補正を行うことができる。
【００１３】
なお、第一及び第二のタイミングにおける燃焼トルクは、随時、制御装置のメモリに記憶しておき、燃焼トルクの補正量の算出後に、燃焼トルクのデータを更新してもよい。
【００１４】
本発明に係る前記トルク推定装置は、前記クランク軸に取付けられたクランク角センサの出力信号から、前記回転速度を算出する回転速度算出手段を備えることがより好ましい。
【００１５】
本発明によれば、内燃機関の制御に用いるクランク角センサを利用することにより、より低コストで、燃焼トルクの推定を行うことができる。
【００１６】
本発明に係る前記トルク推定装置は、前記波形分解手段が、前記クランク軸に取付けられた回転速度センサで検出された回転速度の波形を、周波数分解することを特徴とすることがより好ましく、回転速度を検出することができるのであれば、回転速度センサは、特に限定されるものではないが、より好ましい回転速度センサは、直流発電機である。
【００１７】
本発明によれば、回転速度センサにより、直接的にクランク軸の回転速度を検出し、この回転速度を利用しているので、より精度の良いトルク推定をすることができる。
【００１８】
本発明に係る燃焼エネルギー推定装置は、前記燃焼トルク推定装置を備え、前記燃焼トルクから各気筒の燃焼エネルギーを推定する燃焼エネルギー推定装置であって、該燃焼エネルギー推定装置は、クランク軸の角度を算出するクランク角度算出手段と、前記クランク角度に応じて変化する前記燃焼トルクを積分する積分手段と、該積分した値に基づいて各気筒の燃焼エネルギーを算出する燃焼エネルギー算出手段と、を備える。
【００１９】
本発明によれば、クランク角（位置）の変化にともない、前記クランク軸の回転から燃焼トルクを推定し、その推定した燃焼トルクから燃焼エネルギーを算出（推定）することができるので、気筒内に筒内圧センサを設置する必要がなく、既存のセンサ（クランク角センサとカム角センサ）から安価に燃焼エネルギーを算出することができる。
【００２０】
本発明に係る燃焼エネルギー推定装置は、前記積分手段が、気筒ごとに設定された前記クランク角度の積分区間に基づいて、前記燃焼トルクを積分し、各気筒の燃焼エネルギーを算出することがより好ましい。
【００２１】
本発明によれば、１サイクルにおける、各気筒の筒内圧力の最高圧力は気筒ごとにそのタイミングが異なるので、この点に基づいてクランク角度により定められる積分区間に基づいて、燃焼エネルギーを算出すれば、各気筒の燃焼エネルギーを容易に推定することができる。
【００２２】
本発明に係る燃焼エネルギー推定装置は、前記気筒において設定されたクランク角度の積分区間は、前記気筒の筒内圧が他のすべての気筒の筒内圧に比べて高い筒内圧となる区間であることがより好ましく、別の態様として本発明に係る燃焼エネルギー推定装置は、前記クランク角度の積分区間は、圧縮行程の時点から爆発行程の時点までの区間に設定されることがより好ましい。
【発明の効果】
【００２３】
本発明によれば、クランク軸の回転速度から、連続して変化する燃焼トルクの瞬時値と、各爆発の燃焼エネルギーと、を精度良く、算出することができ、低いコストで燃焼エネルギーを算出することができる。
【図面の簡単な説明】
【００２４】
【図１】第一実施形態に係るトルク推定装置により燃焼トルクを推定する対象となる内燃機関の模式図。
【図２】図１に示す吸気・圧縮・爆発・排気の４つのサイクルにおける筒内圧Ｐｉ、燃焼トルクτ、クランク軸の回転速度ωの関係を示した図であり、（ａ）は、第１気筒の筒内圧を示す図であり、（ｂ）は、第３気筒、第４気筒、第１気筒、第２気筒の順番で爆発しているときの筒内圧を波形で示した図であり、（ｃ）は、（ｂ）から近似される燃焼トルクの波形を示した図であり、（ｄ）は、時間変化するクランク軸回転速度ωを波形で示した図。
【図３】燃焼トルクの推定を説明するための内燃機関の模式図。
【図４】フィルタ出力の誤差を説明する図であり、（ａ）は、その誤差を説明するための図であり、（ｂ）は、誤差から、燃焼トルクの補正量を算出するための基本的な考え方を示した図。
【図５】上下死点における燃焼トルクを説明する図。
【図６】第一実施形態に係る内燃機関のトルク推定装置の模式概念図。
【図７】第一実施形態に係る内燃機関のトルク推定装置のクランク角センサ信号の立下がりに同期して行われる処理のフロー図。
【図８】第一実施形態に係る内燃機関のトルク推定装置の上下死点を検出したときに行われる処理のフロー図。
【図９】クランク角センサ信号の出力の一例を示した図。
【図１０】（ａ）第一実施形態に係る内燃機関のトルク推定装置の回転速度の波形を示した図であり、（ｂ）は第一実施形態に係る内燃機関のトルク推定装置のフィルタ出力値の波形を示した図。
【図１１】フィルタ出力値から補間結果を得るまでの過程を説明した図。
【図１２】フィルタ出力値から補間結果を引いて燃焼トルクを補正する過程を説明する図。
【図１３】第二実施形態に係る内燃機関のトルク推定装置に係るクランク軸の回転角速度と燃焼トルクの周波数成分の関係を説明する図であり、（ａ）は、演算により算出される回転速度ωと、回転速度ωの波形から算出される燃焼トルクτの波形であり、（ｂ）は、クランク軸の回転速度ωと、燃焼トルクτを燃焼に同期した周波数成分の波形、（ｃ）〜（ｅ）は、燃焼周波数に対して、２、３、４倍の周波数の成分に分解したときの各成分の波形を示した図。
【図１４】本発明に係る第二実施形態に係るトルク推定装置の全体構成図。
【図１５】本発明に係る第三実施形態に係るトルク推定装置の全体構成図。
【図１６】本発明に係る第四実施形態に係るトルク推定装置を含む燃焼エネルギー推定装置の全体構成図。
【発明を実施するための形態】
【００２５】
〔第一実施形態〕
以下に、図面を参照して、本発明に係る内燃機関のトルク推定装置のいくつかの実施形態に基づいて説明する。
【００２６】
図１は、第一実施形態に係るトルク推定装置により燃焼トルクを推定する対象となる内燃機関の模式図である。
【００２７】
本実施形態に係る内燃機関１では、ピストン１２が２往復する間に、吸気、圧縮、爆発、排気の４つのサイクルを実行する。図１に示すように、内燃機関１は、ピストン１２の上死点２３から下死点２２に向かってピストン１２が下降するのに同期して吸気弁１３が開くと、スロットル１４で絞られた吸入空気と、インジェクタ１５から噴射された燃料の混合気が気筒２０の燃焼室内に流入する。ピストン１２が下死点２２に達すると、吸気弁１３は閉じ、ピストン１２は上昇する。気筒２０内に閉じ込められた空気はピストン１２によって圧縮される。
【００２８】
ピストン１２が上死点２３に達すると、点火プラグ１６により気筒２０内の混合気が着火され、爆発が始まる。爆発により発生したエネルギーはピストン１２を押し下げ、ピストン１２への圧力はクランク軸１７に伝達され、クランク軸１７を回転させるトルクとなる。
【００２９】
なお、内燃機関１のクランク軸１７には、金属でできた円盤のリングギアが取り付けられており、リングギアの外周には等間隔Δθで金属製の歯が取り付けられている。クランク軸１７の回転は、クランク角センサによって計測される。具体的には、クランク角センサは、リングギア５１の歯の接近に同期して脈動する信号を出力する。
【００３０】
上述したように、気筒内の圧力は、リンク機構を介してクランク軸１７を回転させる。このことを利用して、クランク軸１７の回転速度の変動から気筒内で発生した燃焼エネルギーを推定することが本実施形態の趣旨である。
【００３１】
以下に、本実施形態において、
（１）筒内圧とクランク軸の回転速度の関係を定式化し、クランク角センサにより得られるクランク軸１７の回転速度から筒内圧を推定する方法を説明し、
（２）筒内圧を推定する方法を実現するソフトウエアー構成を説明し、
（３）最後にその動作について述べる。
【００３２】
図２は、吸気・圧縮・爆発・排気の４つのサイクルにおける筒内圧Ｐｉ、燃焼トルクτ、クランク軸１７の回転速度ωの関係を示した図である。図２の（ａ）は、第１気筒の筒内圧を示す図であり、（ｂ）は、第３気筒、第４気筒、第１気筒、第２気筒の順番で爆発しているときの筒内圧を波形で示した図であり、（ｃ）は、（ｂ）から近似される燃焼トルクの波形を示した図であり、（ｄ）は、時間変化するクランク軸回転速度ωを波形で示した図である。なお、４気筒エンジンでは、第３気筒、第４気筒、第１気筒、第２気筒の順に点火される。
【００３３】
図２（ａ）に示すように、第１気筒の吸気行程では筒内圧は吸気管圧（大気圧力）とほぼ同じか、若干低い。その後、ピストン１２の下死点２２から上死点２３までの区間である圧縮行程では、ピストン１２が上死点２３に近づくにつれて、燃焼室内の圧力は大きくなる。
【００３４】
さらに、爆発行程では、行程の開始時点である上死点２３付近で点火されると、圧力はさらに増大し、この筒内圧力がピストン１２を押し下げ体積が膨張し、圧力は徐々に低下し、行程の終了時点である下死点２２で大気圧に近づく。ピストン１２が下死点２２に達すると排気行程が開始され、排気弁１８が開き、気筒２０内の排気ガスは排出される。このときの筒内圧は、大気圧とほぼ同じか、若干高い程度である。このようにして、各気筒についての筒内圧Ｐ１，Ｐ２，Ｐ３，Ｐ４は、図２（ｂ）に示す如き、圧力波形となる。
【００３５】
一方、各気筒により発生する燃焼トルクτは、図２（ｃ）の如き波形となり、このときのクランク軸１７の回転速度ωは、図２（ｄ）の如き波形となる。
【００３６】
ここで、筒内圧Ｐと燃焼トルクτの関係は、クランク軸の回転角θに依存し、以下の式で表せる。
【００３７】
【数１】
2010190174.tif 000003
【００３８】
一方、燃焼トルクτとクランク軸の回転速度ωの関係は以下の式で表せる。
【００３９】
【数２】
2010190174.tif 000004
【００４０】
ただし、数１及び数２に示すパラメータは、
Ａ:各気筒の断面積
ｒ：クランク軸のアーム長
ｌ：コンロッド長
θ：上死点からのクランク軸の回転角
Ｊ：クランク軸系の慣性モーメント
Ｃ：クランク軸に働く摩擦係数
とする。(図３参照)
【００４１】
従って、クランク軸の回転速度ωから燃焼トルクを求めるには、数２の式を計算すればよい。なお、４気筒のエンジンであれば、爆発は４つの気筒で発生しており、爆発の大きさは気筒内の圧力で評価できるため、爆発毎の燃焼エネルギーを求めるには４つの圧力センサが原理的に必要であるが、ここでは、クランク軸の回転速度という一つの物理量から４つの気筒での燃焼エネルギーを求めようとしたものである。
【００４２】
そして、燃焼トルクτから筒内圧Ｐ１〜Ｐ４を求めるには、数１の式を以下のように変形する。
【００４３】
【数３】
2010190174.tif 000005
【００４４】
さらに、最も大きい筒内圧のみに着目し、他は無視することで
【数４】
2010190174.tif 000006
と考えて筒内圧を算出する。各爆発の燃焼エネルギーを表す指標であるＩＭＥＰを近似的に求めるには、必ずしも４つの筒内圧を求める必要はなく、最も大きい筒内圧のみを求めればよい。
【００４５】
ところで、数２の式に基づいて燃焼トルクを求める際に、クランク軸系の慣性Ｊは運転条件に依存しないため、精度よく同定できるが、摩擦係数Ｃは運転条件やエンジンオイルの粘性によって変化するため、精度よく同定できない。このような結果、以下の図４に示すように、実際の燃焼トルクと、算出（推定）した燃焼トルクとには、誤差が生じることになる。
【００４６】
図４（ａ），（ｂ）は、上述した方法で、実燃焼トルクＴ１と、上記演算により算出したエンジントルクＴ２とを示した図であり、（ａ）は、その誤差を説明するための図であり、（ｂ）は、誤差から、燃焼トルクの補正量を算出するための基本的な考え方を示した図である。
【００４７】
図４（ａ），（ｂ）に示すように、数２の式から、燃焼トルクτを求めた結果を示すが、これら図からも明らかなように、低い周波数の誤差が存在し、このような誤差を解消（補正）し、実燃焼トルクに近づけることが必要である。
【００４８】
そこで、その簡易的かつ精度の高い補正として、以下の補正を行う。まず、上死点(θ＝０°)、下死点(θ＝１８０°)のとき、sinθは０となるので、数１の式にθの値を代入すれば、
【００４９】
【数５】
2010190174.tif 000007
となる。
【００５０】
この物理現象を、図５を用いて説明すれば、各気筒のピストンの上死点、下死点においてはコンロッドとクランクアームのなす角θが０°又は１８０゜になる。これにより、燃焼圧力（筒内圧）がクランク軸に燃焼トルクとして伝わらないため、上死点及び下死点でτは０になる。
【００５１】
したがって、数２の式に基づいて、各気筒のピストンの上死点、下死点における算出した燃焼トルクτの値が、誤差そのものである。ここで、図４（ｂ）に示すように、この燃焼トルクの誤差は、燃焼トルクの周波数が低い誤差であるので、燃焼トルクτの上死点、下死点間の誤差は、燃焼トルクτの上死点、下死点の値を直線補間することができる。これによって、実燃焼トルクＴ１と、算出燃焼トルクＴ２の誤差の推定値ＴＤは図４（ｂ）の破線に示すようになり、これを算出燃焼トルクＴから減算する補正を行うことで、より正確な燃焼トルクが算出される。これが、本実施形態における一連の燃焼トルクの算出である。
【００５２】
図６は、本実施形態に係る内燃機関１のトルク推定装置１０Ａの模式概念図である。トルク推定装置１０Ａは、内燃機関を制御するコントロールユニットに含まれるものであり、ＣＰＵ、ＥＰ−ＲＯＭ、ＲＡＭ、及びＡ／Ｄ変換器を含むＩ／ＯＬＳＩ等で構成されている。そして、トルク推定装置１０Ａは、先に示した各種センサからの信号をＩ／ＯＬＳＩで入力として取り込み、この取り込まれた信号に基づいて以下に示す演算処理を実行し、この演算結果として算定された演算結果を出力するように構成されている。
【００５３】
内燃機関１のトルク推定装置１０Ａは、内燃機関１のクランク軸１７の回転速度に基づいて、内燃機関１の燃焼トルクを推定するものである。具体的には、図６に示すように、クランク角センサ３３が、クランク軸に取付けられたリングギア５１の回転に同期して脈動する出力信号し、この出力信号から以下に示すように燃焼トルクを推定している。
【００５４】
尚、クランク角センサ３３の出力信号の一周期分は、リングギア５１の一つの歯に対応する。クランク角センサ信号が、所定の閾値を上から下へ横切る時間の間隔が、Δθの回転所要時間として認識される。
【００５５】
本実施形態に係るトルク推定装置は、回転速度算出手段１０１、フィルタリング手段（燃焼トルク算出手段）１０２、上下死点判定手段１０３、第一及び第二の上下死点フィルタ出力メモリ（第一および第二上下死点燃焼トルク誤差算出手段）１０４、１０５、角度算出手段１０６、補間手段（補正量算出手段）１０７、シフトメモリ（燃焼トルク波形位相シフト手段）１０８、及び、減算手段（燃焼トルク補正手段）１０９を備えている。
【００５６】
回転速度算出手段１０１は、クランク角センサ３３の信号からクランク軸の回転速度を計算するものであり、アナログ信号が所定の閾値を上から下へ横切る度に、今回の横切った時刻と前回の横切った時刻との時間間隔Δｔを計測し、歯の間隔ΔθをΔｔで割ることでエンジンの回転速度（角速度）ω＝Δθ／Δｔを算出し、時間経過に伴う波形として出力するものである。
【００５７】
フィルタリング手段（燃焼トルク算出手段）１０２は、算出された回転速度の高周波成分の波形をフィルタリングして強調することにより、強調された周波数成分の波形から燃焼トルク（補正前の燃焼トルク）を算出するものである。この燃焼トルクの算出は、強調された周波数成分の波形から、先に示した数式により算出することができる。
【００５８】
上下死点判定手段１０３は、クランク角センサ３３又はカム角センサから内燃機関１の気筒における上死点と下死点（上下死点）を判定するものである。なお、この上死点と下死点を判定は、必ずしも気筒ごとに行わなくてもよく、図２及び図４からも明らかなように、少なくとも１つの気筒における上死点と下死点を判定すればよく、ピストンが、上死点、下死点のいずれかに位置することが判定できれば充分である。
【００５９】
第一の上下死点フィルタ出力メモリ（第一の上下死点燃焼トルク誤差算出手段）１０４は、上下死点判定手段１０３の判定結果、及び、フィルタリング手段（燃焼トルク算出手段）で算出された燃焼トルクに基づいて、上下死点（上死点又は下死点いずれか）と判定したタイミング（第一のタイミング）のフィルタの出力、すなわち、数２の式によって算出された燃焼トルク（誤差）の値を記憶する。
【００６０】
第二の上下死点フィルタ出力メモリ（第二の上下死点燃焼トルク誤差算出手段）１０５は、上下死点判定手段１０３の判定結果、及び、フィルタリング手段（波形分解手段）１０２で算出された燃焼トルクに基づいて、第一の上下死点フィルタ出力メモリで記憶してから時間が経過して、次のタイミング（第二のタイミング）すなわち、上下死点（ここで、第一の上下死点フィルタ出力メモリで記憶した値が下死点である場合には、第二の上下死点フィルタ出力メモリで記憶した値は上死点、また、第一の上下死点フィルタ出力メモリで記憶した値が上死点である場合には、第二の上下死点フィルタ出力メモリで記憶した値は下死点）と判定したタイミングのフィルタの出力、すなわち、数２の式によって算出された燃焼トルク（誤差）の値を記憶する。
【００６１】
角度算出手段１０６は、内燃機関１の吸気弁１５または排気弁１８のカム軸の角度を検出するカム角センサ３５とクランク角センサ３３の出力信号からクランク軸１７の角度を算出し、さらに、上下死点判定手段１０３の判定結果に基づいて、第一及び第二の上下死点フィルタ出力メモリに誤差が記憶されるまでの時間（第一のタイミングから第二のタイミングまでの時間）をクランク軸の回転角度として演算する。すなわち、上死点から下死点まで（または下死点から上死点まで）のクランク角センサ３３が発生したパルス数、を数えることで、クランク軸の回転角度を算出する。
【００６２】
補間手段（補正量算出手段）１０７は、角度算出手段１０６で算出されたクランク軸３３の回転角度に応じた燃焼トルク補正量を算出する。具体的には、第一及び第二の上下死点フィルタ出力メモリで記憶された燃焼トルク（誤差）から、算出した回転角度における誤差を、回転角度の増加に伴い直線的に近似して（図４（ｂ）参照）、第一のタイミングから第二のタイミングまでにおける、クランク軸の回転角度に応じた燃焼トルク補正量を波形として算出するものである。
【００６３】
シフトメモリ（燃焼トルク波形位相シフト手段）１０８は、フィルタリング手段（燃焼トルク算出手段）１０２で算出された燃焼トルクの波形の位相（回転角度）を、燃焼トルク補正量の波形の位相に一致させるべく、燃焼トルクの波形の出力を遅延させるものである。これにより、これまでに算出された第一のタイミングから第二のタイミングまでの燃焼トルクに対して、第一のタイミングから第二のタイミングまでにおける燃焼トルク補正量で、補正を行うことができる。
【００６４】
減算手段（燃焼トルク補正手段）１０９は、シフトメモリ（燃焼トルク波形位相シフト手段）１０８により遅延された基本燃料トルクから、補間手段（補正量算出手段）１０７により算出された燃焼トルク補正量の値を減ずる補正をすることで、燃焼トルクを補正し、燃トルクを算出（推定）する。
【００６５】
このように構成された燃焼トルク推定装置１０Ａの演算を図７、図８に示すフローチャートを用いて説明する。
【００６６】
燃焼トルク推定装置１０Ａの動作は、クランク角センサ信号の立下がりに同期して行われる処理と、上下死点を検出したときに行われる処理の２つに分類できる。図７はクランク角センサ信号の立下がりに同期して行われる処理フロー図であり、図８は、上下死点を検出したときに行われる処理のフロー図である。また、図９は、クランク角センサ信号の出力の一例を示した図である。図９（ａ）は、クランク角センサが、等間隔で歯がリングギアについていることにより、クランク角度を検出したときのクランク角センサ信号であり、図９（ｂ）は、リングギアの所定の位置で歯欠けがある場合のクランク角度センサ信号であり、図９（ｃ）は、波形（ｂ）から判定される上下死点信号である。
【００６７】
まず、クランク角センサ３３の立下がりをエッジ検出した後に、ステップ７０１で、回転速度算出手段１０１は、一定周期で発生されるクロックパルスをカウンタで数えることで回転所要時間ΔＴを計測する。再びクランク角センサ信号の立下りが検出されると、カウンタは零リセットされる。その後、カウンタの値はクロックパルスが発生するたびにインクリメントする。クランク角センサの立下りが検出されると、カウンタの値は回転所要時間ΔＴを示しており、回転角Δθをカウンタ値で割ることで回転速度ωが得られる。回転速度ωを計算した後、カウンタ値は零リセットされる。回転速度計算手段１０１ではこれを繰り返し、クランク角センサ信号の立下りのたびに回転速度ωが算出される（図７のステップ７０１）。
【００６８】
図１０（ａ）は、図９（ｂ）に示すクランク角センサ３３が出力する回転角度θから回転速度計算手段が計算した回転速度の一例を示した図である。この図（グラフ）に示すように、クランク角センサ信号が立ち下がるたびに回転速度ωがプロットされる。横軸はクランク角センサ信号の立ち下がり(イベント)の回数である。
【００６９】
ここで、クランク角センサの信号が立ち下がるタイミングで、フィルタリング手段（燃焼トルク算出手段）１０２は、回転速度の波形の高周波を強調するようにフィルタリングし、フィルタリングにより強調された高周波数成分の波形から、燃焼トルクを算出する（図７のステップ７０２）。このようにして、図１０（ａ）に示す回転速度ωをフィルタリングした結果を図１０（ｂ）に示す。このフィルタリングにより算出された値（フィルタ値）は、燃焼トルクに対応するものであり、この例では過去１７イベントの回転速度ωの信号にフィルタの係数の積和をとることで、フィルタリングを行い、燃料トルク（フィルタ出力）を算出している。そのため、最初の１６点は、フィルタ出力がない。
【００７０】
図９（ａ）では、クランク角センサ３３は、等間隔で歯がリングギアについていることにより、クランク角度を検出するものであるが、多くのエンジンでは図９（ｂ）に示すように、リングギアの所定の位置で歯欠けがあり、これによって、上死点からの回転角度を検知している。そこで、上下死点判定手段１０３は、歯欠けからのクランク角センサ信号のパルス数を数えて、例えば図９（ｃ）に示すような上下死点信号を発生させる。
【００７１】
上下死点信号が立ち下がる瞬間に（すなわち、上下死点判定手段１０３により上死点又は下死点であると判定したタイミングにおける）フィルタ出力値（燃焼トルク（誤差）の値）は、第一の上下死点フィルタ出力メモリ（第一の上下死点燃焼トルク誤差算出手段）１０４により記憶され、次に上下死点信号が立ち下がるまで保持される。このとき（次の立下り時点（第二のタイミング））、第一の上下死点フィルタ出力メモリ１０４に記憶されていた値は、第２の上下死点フィルタ出力メモリ（第二の上下死点燃焼トルク誤差算出手段）１０５にシフトされる（図８のステップ８０１）。
【００７２】
このときのフィルタ出力値、及び第一の上下死点フィルタ出力メモリ１０４、第二の上下死点フィルタ出力メモリ１０５に記憶されている値の関係を、図１１を用いて説明する。
【００７３】
クランク角が１８０°のときに上下死点信号が立ち下がる。このとき、フィルタ出力値（燃焼トルクに相当する誤差）ｆ２は第一の上下死点フィルタ出力メモリ１０４に記憶され、出力値ｆ２の値はクランク角３６０゜まで保持されている。なお、クランク角１８０゜のときの上下死点信号が立ち下がる直前の第１の上下死点フィルタ出力メモリ１０４に記憶されている値はｆ１であるが、この値はｆ２が書き込まれる（クランク角度が１８０°になる）直前に第２の上下死点フィルタ出力メモリ１０５に退避されている。このようにして、直近の上下死点とその１つ前の上下死点におけるフィルタ出力の値は保持される。
【００７４】
次に図７及び図１１に示すように、補間手段（補正量算出手段）１０７は、上死点からのクランク角（回転角度）θＣＡに応じて、第一の上下死点フィルタ出力メモリ１０４と第２の上下死点フィルタ出力メモリ１０５に記憶されているフィルタ出力値（誤差）ｆ１，ｆ２を補間すなわち燃焼トルク補正量を算出する（図７のステップ７０４）。補間（燃焼トルク補正量算出）の式は、例えば線形補間であれば、
【数６】
2010190174.tif 000008
である。
【００７５】
なお、補間は上下死点からの角度θＣＡに基づいて行われるため、補間に先立って、角度算出手段１０６で上下死点からのクランク角センサのパルス数をカウントすることで、上下死点からの角度θＣＡを算出する（ステップ７０３）。
【００７６】
フィルタリングの結果となる燃焼トルクの算出値から補間の結果である燃焼トルク補正量を減算することで燃焼トルクは補正されるが、補間の結果はクランク軸が１８０゜回転するのに要する時間だけ遅れる。すなわち、図１１のフィルタ出力値（燃焼トルクの値）の波形と、一番下の補間結果の値（燃焼トルク補正量の値）の波形を比較すればわかる。
【００７７】
このような点から、シフトメモリ（燃焼トルク波形位相シフト手段）１０８は、クランク軸が１８０°回転する間、フィルタ出力値の波形の位相をクランク角度１８０°遅延させる（図７のステップ７０５）。このように、シフトメモリ１０８で遅延させたフィルタ出力値の波形の結果は、図１２の上から２番目の図のようになる。
【００７８】
減算手段（燃焼トルク補正手段）１０９において、シフトメモリ１０８が出力するフィルタ出力値の波形を１８０゜クランク回転分遅延させた結果から、補間手段の結果を減算することで、燃焼トルクを算出する（図７のステップ７０６）。その結果が図１２の一番下の図で、一番上の図に示すフィルタ出力値は必ずしも上下死点で０でなかったのに対し、一番下の図では上下死点の燃焼トルク推定値が０になっている。
【００７９】
以上のような燃焼トルク推定装置によって、クランク角センサの信号から燃焼トルクの瞬時値をより正確に推定することができる。この推定結果は、上下死点における燃焼トルクの誤差に基づいて補正するものであるため、摩擦係数等の変化にもロバストである。
【００８０】
〔第二実施形態〕
第二実施形態と第一実施形態の相違する点は、フィルタリング手段（燃焼トルク算出手段）である。具体的には、フィルタリング手段は、クランク軸の回転速度と燃焼トルクとの波形を周波数成分に分解し、それぞれの振幅と位相の関係を調べる。
【００８１】
具体的には、クランク軸の回転速度（回転角速度）ωを用い、燃焼トルクτは、正弦波（余弦波）の波形の重ね合わせで近似できることから、この波形に着眼し、数７、数８のようにクランク角速度ωと燃焼トルクτを周波数分解すると、以下に示すような成分に分解することができる。
【００８２】
【数７】
2010190174.tif 000009
【００８３】
【数８】
2010190174.tif 000010
【００８４】
ここで、数７に示すω１ｓｉｎ（ｆｔ＋α１）は、燃焼周波数成分（基本成分）であり、燃焼周波数（燃焼に同期した周波数）とは、単位角度ごとの燃焼頻度を表した角度周波数であり、例えば４気筒の場合には燃焼周期（燃焼角度周期）が１８０°ＣＡであり、内燃機関の気筒数に対応する燃焼角度周期の逆数により決定される周波数である。また、ωｎｓｉｎ（ｆｔ＋αｎ）は、ｎ倍成分の波形で、ｎは自然数であり、燃焼周波数のｎ倍の回転速度の周波数成分の波形である。
【００８５】
図１３（ａ）は、演算により算出される回転速度ωと、回転速度ωの波形から算出される燃焼トルクτの波形であり、（ｂ）は、クランク軸の回転速度ωと、燃焼トルクτを燃焼に同期した周波数成分の波形、（ｃ）〜（ｅ）は、燃焼周波数に対して、２、３、４倍の周波数の成分に分解したときの各成分の波形である。
（τの燃焼周波数成分の振幅／ωの燃焼周波数成分の振幅）をｒとすると，
（τの燃焼周波数の2倍成分の振幅／ωの燃焼周波数の2倍成分の振幅）＝r×2
（τの燃焼周波数の3倍成分の振幅／ωの燃焼周波数の3倍成分の振幅）＝r×3
（τの燃焼周波数の4倍成分の振幅／ωの燃焼周波数の4倍成分の振）＝r×4
となる。
【００８６】
また、
（τの燃焼周波数成分とωの燃焼周波数成分  の位相差）＝45゜CA
（τの燃焼周波数の2倍成分とωの燃焼周波数の2倍成分の位相差）＝45゜CA/2
（τの燃焼周波数の2倍成分とωの燃焼周波数の2倍成分の位相差）＝45゜CA/3
（τの燃焼周波数の2倍成分とωの燃焼周波数の2倍成分の位相差）＝45゜CA/4
となる。
【００８７】
したがって，クランク軸の回転速度ωから燃焼トルクτを推定するには、図１４に示すように、まず燃焼周波数フィルタ１４０１，２倍周波数フィルタ１４０２，３倍周波数フィルタ１４０３，４倍周波数フィルタ１４０４，‥‥からなるフィルタリング手段（周波数分解手段）１４０５によりωを燃焼に同期した成分，燃焼の２倍周期成分，３倍周期成分，４倍周期成分，‥‥に周波数分割し，振幅拡大手段１４０６にてそれぞれの振幅をｒ倍，２ｒ倍，３ｒ倍，４ｒ倍し，それぞれの位相を４５゜ＣＡ，４５゜ＣＡ／２，４５゜ＣＡ／３，４５゜ＣＡ／４進ませたいが、位相は遅らせることしかできないので、位相遅延手段１４０７にて４５゜ＣＡ−４５゜ＣＡ，４５゜ＣＡ−４５゜ＣＡ／２，４５゜ＣＡ−４５゜ＣＡ／３，４５゜ＣＡ−４５゜ＣＡ／４遅らせて、４５゜ＣＡ遅れたトルクを推定する。
【００８８】
そして、これらの燃焼トルクτの波形（図１３（ｂ）〜（ｅ））を重ね合わせることで図１３（ａ）に示すように、クランク軸の回転速度から、燃焼トルクを算出することができる。
【００８９】
〔第三実施形態〕
図１５は、本発明に係る第三実施形態に係るトルク推定装置の全体構成図である。
第一実施形態では、クランク角センサ３３の信号を回転速度算出手段１０１にて処理してΔθの回転所要時間Δｔを算出し、Δθ／Δｔによりクランク軸１７の回転速度を算出していた。このかわりに、第三実施形態では、直接的にクランク軸１７の回転速度を計測している。第三実施形態は、第一実施形態に比べて、回転速度の算出方法が相違し、その他の構成は、第一実施形態と同じであるので、同じ機能を有する構成は、同一の符号を付して、詳細な説明は省略する。
【００９０】
図１５に示すように、第三実施形態に係るトルク推定装置１０Ｂは、図６に示す第一実施形態のトルク推定装置１０Ａに係る内燃機関１に、回転速度センサ３４を加えたものである。回転速度センサ３４として、例えば直流発電機を用いれば、回転速度に比例した電圧が発生する。この場合、従来のクランク角センサ３４からの信号は、フィルタ内に記憶されるデータを更新するタイミングをとるために利用される。つまり、クランク角センサが立ち下がるときに、図７に示すフローは実行される。
【００９１】
このような実施形態をとるのは、第一実施形態のように、Δθ／Δｔを算出すると、Δｔの精度をかなり良くしないとエンジンの回転速度ωは正確に求まらない場合がある。本実施形態では、直接エンジン回転速度を求めるので、回転速度ωに対して、それほどの精度はいらない。そこで、ここでは回転速度センサ３４を用いることで、Δｔへの要求精度を緩和することができ、トルク推定にかかる演算処理を軽減することができる。
【００９２】
〔第四実施形態〕
図１６は、本発明に係る第四実施形態に係るトルク推定装置を含む燃焼エネルギー推定装置の全体構成図である。第一実施形態〜第三実施形態は、燃焼トルクの瞬時値を求めるものであった。気筒ごとの燃焼トルクがわかっていれば、これを１エンジンサイクル分積分することでエンジンの燃焼エネルギーが求められる。本実施形態では、図１６のような構成で、エンジンの燃焼エネルギーを求める。
【００９３】
図１６に示すように、この燃焼エネルギー推定手段１０Ｄは、図６に示す第一実施形態のトルク推定装置１０Ａに、クランク角度算出手段１６１、積分手段１６２、及び、燃焼エネルギー算出手段１６３をさらに備えたものである。
【００９４】
クランク角度算出手段１６１は、内燃機関１の吸気弁１３または排気弁１８のカム軸の角度を検出するカム角センサとクランク角センサ３３の出力信号からクランク軸１７の角度を算出している。具体的には、クランク角度算出手段１６１は、内燃機関１の吸気弁１３又は排気弁１８に配置されたカム軸に対して、その回転角度を検出するカム角センサ３５の信号を検出する。そして、各気筒２０のカム軸のカム角が示す、ピストン１２の上死点２３の位置と、クランク角センサ３３による回転角度の情報から、クランクの角度を算出している。このクランク角度を算出する目的は、各気筒を判別し、推定された燃焼トルクから、各気筒ごとにおける燃焼エネルギーを算出するための時間的なパラメータとして用いるためである。
【００９５】
積分手段１６２は、クランク角度算出手段１６１で算出されたクランク角度に応じて変化する前記燃焼トルクを積分している。具体的には、所与のクランク角からクランク角までの区間、燃焼トルクの瞬時値を積分している。
【００９６】
燃焼エネルギー算出手段１６３は、あらかじめ与えた積分値と燃焼エネルギーの対応関係から燃焼エネルギーを算出している。積分する積分区間の考え方を、図２に戻り説明する。前述したように、図２は、第３気筒、第４気筒、第１気筒、第２気筒の順番で爆発しているときの筒内圧（図２（ｂ））、燃焼トルク（図２（ｃ））、クランク軸回転速度（図２（ｄ））を表している。
【００９７】
例えば、第１気筒の燃焼エネルギーを求めるには、積分手段１６２は、この気筒の筒内圧が、他のすべての気筒（第２気筒から第４気筒）の筒内圧に比べて高い（最も大きい）θ４１からθ１３まで区間を積分するように設定されており、この積分区間における燃焼トルクを積分する。筒内圧を測定できる場合にはその測定値に基づいて、この区間を前記条件により設定してもよい。この区間の燃焼トルクは，第１気筒の寄与が大きく、燃焼トルクの合計≒第１気筒の燃焼トルク、と近似できる。
【００９８】
別の言い方をすれば、積分区間の開始点は圧縮行程の所定の時点に設定され、終了点は爆発行程の所定の時点に設定される。この区間の燃焼トルクは、第１気筒の寄与が大きく、燃焼トルクの合計≒第１気筒の燃焼トルク、と近似できる。
【００９９】
燃焼エネルギーを求めるには１エンジンサイクル分の燃焼トルクを積分する必要があるが、［θ４１、θ１３］のみの積分値Ｗ１をもとめておいて、積分値と燃焼エネルギーの関係Ｗから第１気筒の燃焼エネルギーを求める。多くの場合、Ｗ１とＷは、直線の関係で近似できるので、燃焼エネルギー算出手段には、直線の傾きと切片を記憶しておけばよい。
【０１００】
このように、第一〜第三実施形態のいずれかで算出された燃焼トルクを、筒ごとに設定された前記クランク角度の積分区間に基づいて積分し、この積分された値を、所定の対応関係を示す直線関係等の関係式で変換することで、各気筒の燃焼の燃焼エネルギーが算出できる。
【符号の説明】
【０１０１】
１：内燃機関、１２：ピストン、１３：吸気弁、１４：スロットル、１５：インジェクタ、１６：点火プラグ、１７：クランク軸、１８：排気弁、２０：気筒、２２：下死点、２３：上死点、３３：クランク角センサ、３４：回転角センサ、１０Ａ〜１０Ｃ：トルク推定装置、１０Ｄ：燃焼エネルギー推定装置、１０１：回転速度算出手段、１０２：フィルタリング手段（燃焼トルク算出手段）、１０３：上下死点判定手段、１０４：第一の上下死点フィルタ出力メモリ（第一の上下死点燃焼トルク誤差算出手段）１０５：第二の上下死点フィルタ出力メモリ（第二の上下死点燃焼トルク誤差算出手段）、１０６：角度算出手段、１０７：補間手段（補正量算出手段）、１０８：シフトメモリ（燃焼トルク波形位相シフト手段）、１０９：減算手段（燃焼トルク補正手段）
【図１】
2010190174.tif 000011
【図２】
2010190174.tif 000012
【図３】
2010190174.tif 000013
【図４】
2010190174.tif 000014
【図５】
2010190174.tif 000015
【図６】
2010190174.tif 000016
【図７】
2010190174.tif 000017
【図８】
2010190174.tif 000018
【図９】
2010190174.tif 000019
【図１０】
2010190174.tif 000020
【図１１】
2010190174.tif 000021
【図１２】
2010190174.tif 000022
【図１３】
2010190174.tif 000023
【図１４】
2010190174.tif 000024
【図１５】
2010190174.tif 000025
【図１６】
2010190174.tif 000026