(57)【特許請求の範囲】
【請求項１】
クランクシャフトに連結された複数の気筒を備えた内燃機関を制御する内燃機関の制御装置において、
前記内燃機関におけるクランク角度と前記クランク角度での時刻とを検出するクランク角度検出手段と、
前記クランク角度検出手段から前記内燃機関の燃焼を含む期間に対応する前記クランク角度と前記時刻を読み出し、クランク角速度及び角加速度を推定する角速度及び角加速度推定手段と、
前記内燃機関のピストン、コンロッド及び前記クランクシャフトを含む系の運動方程式から、燃焼により発生するトルクが零として上死点近傍の複数の標本点での演算を行い、上死点上の外部負荷トルクを推定する外部負荷トルク推定手段と、
前記外部負荷トルク推定手段で推定した前記外部負荷トルクとクランク角速度及び角加速度を前記運動方程式に適用して推定した燃焼によるガス圧トルクから、前記内燃機関の筒内圧力を推定する筒内圧力推定手段と、
を備えたことを特徴とする内燃機関の制御装置。
【請求項２】
前記外部負荷トルク推定手段は、上死点上の外部負荷トルクの推定を線形近似または多項式近似を利用して行うことを特徴とする請求項１に記載の内燃機関の制御装置。
【請求項３】
前記外部負荷トルク推定手段は、上死点近傍の外部負荷トルクの分布において、上死点が進角側の２点間に存在する隣接した３点の標本点の演算値から２次近似により算出した外部負荷トルクと、上死点が遅角側の２点間に存在する隣接した３点の標本点から２次近似により算出した外部負荷トルクの加重平均を求めて最終的な外部負荷トルクを推定することを特徴とする請求項１または２に記載の内燃機関の制御装置。
【請求項４】
前記クランク角度検出手段は、あらかじめ前記クランク角度検出手段により割り付けられたクランク角位置と、実際のクランク位置との誤差量を把握して補正を掛け、適正な上死点上での外部負荷トルク推定を行うクランク角度補正手段を備えたことを特徴とする請求項１から３の何れか一項に記載の内燃機関の制御装置。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、内燃機関の制御装置に係り、特には内燃機関のトルクや筒内圧などの燃焼状態を推定する内燃機関の制御装置に関するものである。
【背景技術】
【０００２】
内燃機関の燃費性能やエミッション性能を向上させるために、内燃機関の燃焼状態を計測し、その計測結果をフィードバックさせて制御する方法が有効である。そのためには、内燃機関の燃焼状態を正確に計測することが重要である。内燃機関の燃焼状態は筒内圧を計測することにより、正確に計測できることが広く知られている。筒内圧の計測方法では筒内圧センサ信号から直接測定する方法の他に、クランク角度信号などのエンジンにおける各機構の情報からガス圧トルクを推定する方法がある。
【０００３】
従来では、例えば特開２００９−２７５６１８号公報（特許文献１）に、クランク角度センサの出力信号から燃焼状態を推定する燃焼状態推定装置が開示されている。特許文献１に開示された燃焼状態推定装置では、クランク角度センサの信号からクランク角速度及び角加速度を推定し、上下死点における外力トルクと推定した角加速度から求めた軸トルクとの関係により外力トルクを推定し、それを基にガス圧トルクを推定する。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００９−２７５６１８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、クランク角度センサの信号から得られるクランク角度とその時刻データは離散的であり、連続的な角加速度演算を行うことが出来ない。また、クランク軸とピストンの中心位置にはオフセットが存在することから上死点はクランク軸の頂点から微小にずれるため、標本点と上死点が重ならない場合が多い。そのため上死点が標本点上に存在しない場合、外部負荷トルクを推定する際に上死点上での演算を行うことが出来ない。このため、筒内圧力を推定し、筒内圧から燃焼状態を推定する場合において筒内圧波形に上死点近傍の誤差が大きくなりスパイクや段差が発生する。
【０００６】
この発明は、従来の内燃機関の制御装置に於ける前述のような課題を解決するためになされたものであり、複数気筒を有する内燃機関に於いて、気筒ごとの燃焼状態を精度よく推定することのできる内燃機関の制御装置を提供することを目的とするものである。
【課題を解決するための手段】
【０００７】
この発明に係る内燃機関の制御装置は、クランクシャフトに連結された複数の気筒を備えた内燃機関を制御する内燃機関の制御装置において、
前記内燃機関におけるクランク角度と前記クランク角度での時刻とを検出するクランク角度検出手段と、前記クランク角度検出手段から前記内燃機関の燃焼を含む期間に対応する前記クランク角度と前記時刻を読み出し、クランク角速度及び角加速度を推定する角速度及び角加速度推定手段と、前記内燃機関のピストン、コンロッド及び前記クランクシャフトを含む系の運動方程式から、燃焼により発生するトルクが零として上死点近傍の複数の標本点での演算を行い、上死点上の外部負荷トルクを推定する外部負荷トルク推定手段と、前記外部負荷トルク推定手段で推定した前記外部負荷トルクとクランク角速度及び角加速度を前記運動方程式に適用して推定した燃焼によるガス圧トルクから、前記内燃機関
の筒内圧力を推定する筒内圧力推定手段とを備えたものである。
【発明の効果】
【０００８】
この発明に係る内燃機関の制御装置によれば、標本点上に上死点が存在しない場合でも、上死点で燃焼により発生するトルクが零となるように精度良く演算を行うことが可能となる。このため、精度の良いガス圧トルクと共に外部負荷トルクを推定する外部負荷推定処理、燃焼パラメータ推定処理を含む内燃機関の制御装置を提供することができる。
【図面の簡単な説明】
【０００９】
【図１】この発明の実施の形態１に係る内燃機関の制御装置を含む内燃機関の概略構成図である。
【図２】この発明の実施の形態１に係る内燃機関の制御装置における燃焼状態推定処理部の動作を示すフローチャートである。
【図３】この発明の実施の形態１に係る内燃機関のピストンとクランクシャフトの力学的な関係を示す断面図である。
【図４】この発明の実施の形態１に係る内燃機関の制御装置での燃焼状態推定処理部の外部負荷トルク推定手段における近似方法を説明する図である。
【図５】この発明の実施の形態１に係る内燃機関の制御装置の燃焼状態推定処理部で推定した筒内圧力と他の内燃機関の制御装置の燃焼状態推定処理部により推定した筒内圧力を比較した図である。
【図６】この発明の実施の形態１に係る内燃機関の制御装置の燃焼状態推定処理部で推定した熱発生率と他の内燃機関の燃焼状態推定処理部により推定した熱発生率を比較した図である。
【発明を実施するための形態】
【００１０】
以下、この発明に係る内燃機関の制御装置の好適な実施の形態について図面を参照して詳細に説明する。
【００１１】
実施の形態１．
図１は、この発明の実施の形態１に係る内燃機関の制御装置を含む内燃機関の概略構成図である。図１において、内燃機関１の第１気筒＃１、第２気筒＃２及び第３気筒＃３のそれぞれの気筒で発生した燃焼ガス圧は、ピストン２の上面で受け止められ、コンロッド３を介して一つのクランクシャフト４を回転させる動力としてクランクシャフト４に伝達される。そして、クランクシャフト４の回転によって動力として出力される出力トルクに基づいて車輛が駆動される。
【００１２】
クランク角度センサ５は、クランクシャフト４の軸方向の前端部に取り付けられたクランク角度センサ用のシグナルロータ６の側面に対して垂直方向に設置されている。シグナルロータ６には複数の歯が刻まれており、歯の凹凸によるクランク角度センサ５とシグナルロータ６の側面との距離の変化を読み取り、矩形波の信号としてＥＣＵ（Engine Control Unit）７に入力される。この矩形波の周期時間や立下り位置などによってクランク軸情報を得ることが出来る。
【００１３】
カムシャフト８の角度を検知するため、カム角度センサ用のシグナルロータ９の歯の位置を示すカム角度センサ１０の信号をＥＣＵ７に入力し、この２種類の矩形波信号の関係からクランクシャフト４の角度を検出する。カムシャフト８はクランクシャフト４とチェーン１１で連結されており、クランクシャフト４が２回転する間にカムシャフト８は１回転するように構成されており、カム角度を検出することによりエンジンの行程を見極めることが出来る。
【００１４】
また、フライホイール１２におけるリングギアについても、第２クランク角度センサ１３を用いてクランク角度センサ５と同様に歯の信号を計測するようにしてもよい。これはフライホイール１２の質量がシグナルロータ６が位置するエンジン前端部と比較して大きいため高周波振動が抑えられる他、シグナルロータ６よりも歯数が多く、欠け歯もないため色々な特定のクランク角度領域における標本点数が多く、精度の高い演算を行うことが出来るためである。
【００１５】
インマニ圧センサ１４ではインテークマニホールドの圧力であるインマニ圧を測定し、大気圧センサ１５では大気圧を計測する。大気圧センサ１５はＥＣＵ７に内蔵してもよい。
【００１６】
ＥＣＵ７には各種センサからの信号が入力され、得られた情報から運転状態を把握し、ＥＣＵ７がそれに応じて各種アクチュエータの制御を行うことにより燃費性能とエミッション性能の向上を実現する。
【００１７】
図２は、実施の形態１に係る内燃機関の制御装置における燃焼状態推定処理部の動作を示すフローチャートである。図２において、ＥＣＵ７に設けられた燃焼状態推定処理部は大まかに分けて、クランク角度検出手段１６と、クランク角度補正手段１７と、角速度及び角加速度推定手段１８と、外部負荷トルク推定手段１９と、筒内圧力及び燃焼パラメータ推定手段２０を備えて構成されている。
【００１８】
クランク角度検出手段１６は、クランク角度センサ５における矩形波信号に現れる欠け歯の位置と、その時のカム角度センサ１０で測定した矩形波信号の関係から基準となるクランク角度を検出する。そして第２クランク角度センサ１３で測定した矩形波信号における基準角度からの立下り、或いは立上がりの回数とその時刻からクランク角度とその時刻を把握し、所定期間だけ記憶する。
【００１９】
前記所定期間は、燃焼状態を推定している対象の気筒が、燃焼により熱量を発生する期間が含まれていれば問題ないが、ＦＦＴ（高速フーリエ変換）による周波数操作や角速度角加速度を推定する場合は２の乗数の標本点数で行うため、その点数を超える角度領域の時刻を記憶する必要がある。
【００２０】
クランク角度補正手段１７は、あらかじめクランク角度検出手段１６により割り付けられたクランク角度と、実際のクランク角度との誤差量を把握して補正を掛け、適正な上死点上での外部負荷トルク推定を行う。即ち、クランク角度補正手段１７は、燃焼させずにクランクシャフト４を回転させ、シリンダ内部におけるクランク角度に対する筒内圧力を計測し、その圧力波形の頂点におけるクランク角度と幾何学的に推定した圧縮過程における上死点でのクランク角度との差をあらかじめ把握しておき、適正なクランク角度になるように補正を掛ける。
【００２１】
角速度及び角加速度推定手段１８は、まず、クランク角度検出手段１６により計測したクランク角度を引数とする時刻データから、燃焼状態推定の演算対象となる燃焼を含む期間の時刻データを読み出す。その時刻データを基に生成した周期関数をＦＦＴにより周波数領域の関数に変換し、ノイズとみなす高周波成分を除去、もしくは成分量を補正する。ＦＩＲ（Finite Impulse Response）フィルタや、ＩＩＲ（Infinite Impulse Response）フィルタなどのデジタルフィルタによる周波数操作を用いても良い。
【００２２】
そして、周波数空間における代数的微分を利用して、前記クランク角度を引数とするクランク角速度及び角加速度を推定する。またはＦＩＲフィルタやＩＩＲフィルタなどのデジタルフィルタを用いた場合は、差分法による微分を行い、前記クランク角度を引数とするクランク角速度及び角加速度を推定する。
【００２３】
外部負荷トルク推定手段１９は、演算対象の気筒における上死点近傍では燃焼により発生するトルクは零であるという仮定の元、上死点近傍における複数の標本点についてクランク軸周りの運動方程式から各標本点上の外部負荷トルクを推定し、線形近似または多項式近似を利用した演算を行い、上死点上での外部負荷トルクを推定する。
【００２４】
筒内圧力及び燃焼パラメータ推定手段２０は、クランク軸周りの運動方程式に外部負荷トルク推定手段１９により推定した外部負荷トルクを適用して燃焼ガス圧トルクを推定し、クランク角度とピストン２の位置関係から燃焼によってピストン２を押し下げる力を推定し、ピストン面積で割ることで筒内圧を推定する。
【００２５】
更に、筒内圧力と筒内容積を基にクランク角度位置に対する発生熱量を推定し内燃機関１の燃焼状態を示す燃焼パラメータである質量燃焼割合（ＭＦＢ）を推定する。
【００２６】
次に、図１に示した本実施の形態に係る内燃機関の制御装置が燃焼状態推定処理の過程で、クランク角速度および角加速度の推定、エンジンから出力されるトルクの推定、エンジンの燃焼状態の推定を実施する方法について、図２のフローチャートを用いて具体的に説明する。
【００２７】
図２に於いて、ステップＳ２０１では、クランク角度センサ５から読み取った信号の欠け歯位置、及びその際のカム角度センサ１０の矩形波信号からクランクシャフト４の角度位置を検出する。同時に第２クランク角度センサ１３によりフライホイール１２のリングギアの信号を検出する。クランク角度センサ５から読み取っているシグナルロータ６の歯は例えば３６歯システムでは１０［ｄｅｇ］間隔で並んでいるが、リングギア１２は例えば４［ｄｅｇ］間隔の周期で歯が９０歯並んでおり、より詳細な角度位置を検出することが出来る。
【００２８】
次に、ステップＳ２０２では演算対象の気筒について、燃焼を含む期間に対応するクランク角度と検出時刻を読み出し、クランク角度θｊに対する時刻のデータ系列t（θｊ）を構成する。このデータ系列は燃焼状態を推定している対象の気筒が、燃焼により熱量を発生する期間が含まれていることが条件であり、読み出す期間を広げ、サンプル点数を拡大して燃焼状態を推定することも可能である。
【００２９】
次に，ステップＳ２０３では、あらかじめピストン２の内部に筒内圧センサを設置し、モータリング計測などの方法から未燃時の筒内圧のピーク位置に対するクランク角度を検出する。未燃時はシリンダ容積が最小となる上死点において筒内圧が最大となることを利用して、その検出したクランク角度と上死点位置との差異δθを調べ、適正なクランク角位置になるように式（１）のように補正を掛ける。クランク角度θｊは補正後のクランク角、θｏｕｔ_ｊは補正前である検出直後のクランク角を示す。もしくはクランク角度と上死点位置との差異δθを把握するために燃焼解析装置により上死点位置を把握してもよい。
【００３０】
【数１】
0006029726.tif 000002
【００３１】
次に、ステップＳ２０４では、クランク角度θｊとそれに対する時刻のデータ系列ｔ（θｊ）についてインデックスをｊ＝１，２，・・・，ｋ，・・・，Ｎとする。Ｎは標本点数であり、ＦＦＴによる周波数処理、クランク角速度及び角加速度の推定を行う場合は２の乗数点数で行う必要がある。本実施の形態においてはＦＦＴによる周波数処理、クランク角速度及び角加速度の推定を行う。ＦＦＴにより周波数操作を行うため、データ系列ｔ（θｊ）を式（２）のように周期関数と見なせる関数ｘ（θｊ）に変換する。
【００３２】
【数２】
0006029726.tif 000003
ここで、αは次式（３）で示される。
【００３３】
【数３】
0006029726.tif 000004
【００３４】
式（２）（３）から求めた周期関数ｘ（θｊ）についてＦＦＴを行い、周波数領域の関数Ｘ(ω)（|ω|≦π／Δθ）に変換する。その際にノイズが存在する周波数帯域の成分量を除去又は補正することによってノイズを低減させる。周波数成分調整を行ったＸ(ω)を＜hat＞Ｘ(ω)（Ｘの上に“ ^ ” を乗せた文字を文章中では＜hat＞Ｘと表記する）とおく。
【００３５】
次にＳ２０５では、ステップＳ２０４で周波数成分を調整して求めた＜hat＞Ｘ(ω)から代数的微分を行うことによりクランク角速度、角加速度を求める。
クランク角度θｊの関数としての角速度＜dot＞θ（θｊ）（θの上に“ ・ ” を乗せた文字を文章中では＜dot＞θと表記する）は式（４）から求めることが出来る。
【００３６】
【数４】
0006029726.tif 000005
ｘ（θｊ）のθに関する微分は、ＦＦＴにより算出した＜hat＞Ｘ(ω)に（ｉ・ω）を掛けた後に逆ＦＦＴすることにより求められる（ｉを虚数単位とおく）。
【００３７】
クランク角度θjの関数としての角加速度＜２dot＞θ（θｊ）（θの上に“ ・・ ” を乗せた文字を文章中では＜２dot＞θと表記する）は式（５）から求めることが出来る。
【００３８】
【数５】
0006029726.tif 000006
【００３９】
ｘ（θｊ）のθに関する２階微分は、ＦＦＴにより算出した＜hat＞Ｘ(ω)に（ｉ・ω）２を掛けた後に逆ＦＦＴすることにより求められる。
【００４０】
次に、ステップＳ２０６、ステップＳ２０７では、ステップＳ２０４で推定したクランク角速度と角加速度をピストン２、コンロッド３、及びクランクシャフト４を含む系の運動方程式（６）に適用し、上死点における条件として燃焼により発生するトルクは零となることから外部負荷トルクを推定する。図３はクランクシャフト４に対して垂直なエンジンの断面図である。図３に示すように各変数を設定する。
【００４１】
【数６】
0006029726.tif 000007
【００４２】
ここでＩは慣性モ−メント、Ｓはピストン２の頂面の投影面積、ｍｐはピストン２の質量（往復質量）、ｒはクランクシャフト４の半径とする。また、ｉは気筒番号を示し、ｎは気筒数（本実施の形態はｎ＝３）、Ｐｃｙｌ_ｉはｉ番気筒の筒内圧力、ａｉはｉ番気筒のピストン２の往復方向の加速度、θｉはｉ番気筒のクランク角度、φｉはｉ番気筒のコンロッド３の角度を表す。また、外部負荷トルクをＴｌｏａｄとする。式（６）における右辺の第一項は各気筒において筒内圧変化によりピストン２を押し下げる力、ピストン２の往復慣性により発生する力が変換されたトルクを推定しており、第二項はフリクションや補機負荷、走行抵抗などエンジンの外部からクランクシャフト４にかかる負荷を外部負荷トルクとして演算を行う。なお、Ｒｉはピストン２にかかる力に対するトルク変換係数を示している。
【００４３】
前記式（６）のＰｃｙｌ_ｉにおいて、燃焼が行われていない気筒の筒内圧力に関してはインマニ圧センサ１４で計測した吸気行程のインマニ圧、大気圧センサ１５で計測した排気行程の排気圧（大気圧）、初期筒内圧から推定することが出来、ピストン２の往復方向における加速度はクランクシャフト４の角速度から幾何学的に推定することが出来る。
ｂ番気筒において燃焼が行われている場合、式（７）のように表すことが出来、式（７）は式（７Ａ）のように変換出来る。この式（７Ａ）の右辺第一項は「気筒ｂの燃焼によるガス圧トルク」を示し、第二項は「未燃焼の気筒におけるガス圧トルクの総和」を示している。また、第三項は、「ピストンの往復慣性により発生するトルク」を示し、第四項が「外部負荷トルク成分」となる。筒内のガス圧にピストン面積を掛けてピストン面の垂直方向に働く力が求まり、トルク変換係数Ｒを乗算することによってクランクシャフトにおける回転方向のトルク成分に変換される。従って、式（７Ａ）の右辺第一項のＰｃｙｌｂ・Ｓ・Ｒｂがガス圧トルクに相当するパラメータとなる。
【００４４】
【数７】
0006029726.tif 000008
【００４５】
ここで、Ｐｕｂ_ｉはｉ番気筒の未燃時の筒内圧力を表す。式（７）はｂ番気筒における上死点上では式（８）として表すことが出来る。
【００４６】
【数８】
0006029726.tif 000009
【００４７】
上記の通り未燃時の筒内圧、ピストン２の加速度は各入力情報から推定できるため、式（８）により上死点における外部負荷トルクを推定することが出来る。外部負荷トルクは１サイクルの燃焼の期間において大きな変化がないことから、本実施の形態では外部負荷トルクは１回の燃焼における演算毎に上死点で算出した値を一定値として扱う。
【００４８】
ただし、ピストン２がクランクシャフト４に対して微小にオフセットされており、上死点はクランク角度が零［ｄｅｇ］より遅角側となる。上死点が標本点上に位置しないため、式（８）から直接的にＴｌｏａｄを算出できないので、上死点近傍においても式（８）が成り立つと仮定し、上死点近傍における複数の標本点での演算結果の分布から近似を行い、上死点上での適正なＴｌｏａｄを推定する。
【００４９】
近似方法は上死点を跨ぐ隣接した２点間から線形近似による推定方法である式（９）で最も簡単にＴｌｏａｄを求めることが出来る。ここでのｉは上死点直前の標本点におけるクランク角度、Δθは検出するクランク角度の周期間隔、θＴＤＣは上死点のクランク角度を表す。
【００５０】
【数９】
0006029726.tif 000010
【００５１】
しかし、式（８）のクランク角度に対する振る舞いが直線的でないため、精度が粗く、推定した筒内圧力波形に上死点近傍で段差が生じる。また、隣接した３点における２次近似からＴｌｏａｄの推定を行った場合でも３点間で大きな曲率変化がある場合には筒内圧
波形に段差が生じる。
【００５２】
出来る限り少ない演算量で実測した筒内圧波形との誤差を解消するために、本実施の形態では、図４のように２通りの上死点を領域に含む隣接した３点から２次近似を行い、推定したＴｌｏａｄの重み平均を取ることで非常に精度良くＴｌｏａｄを推定する。
【００５３】
次に、Ｔｌｏａｄ推定の詳細を説明する。まず遅角側の２点間が上死点を跨ぐ場合の隣接した３点において２次近似を行う。各標本点において、式（８）から算出したＴｌｏａｄ(θ)を用いて式(１０)のような二次方程式における各項の係数を式（１１）から式（１３）により算出する。算出した係数を適用した式（１０）に上死点であるクランク角度を代入し、Ｔｌｏａｄ（θＴＤＣ）を推定し、これをＴｌｏａｄ_１とする。
【００５４】
【数１０】
0006029726.tif 000011
【００５５】
次に参照する標本点を遅角側にずらし、進角側の２点間が上死点を跨ぐ場合の隣接した３点において同様に２次近似を行い、Ｔｌｏａｄ（θＴＤＣ）を推定する。これをＴｌｏａｄ_２とする。そして式（１４）にＴｌｏａｄ_１、Ｔｌｏａｄ_２を適用し、２つの演算結果における加重平均を取り、それをＴｌｏａｄとして筒内圧を推定する。
【００５６】
【数１１】
0006029726.tif 000012
【００５７】
ここで、ｇは０≦ｇ≦１の係数であり、基本的に０．５として扱うが適宜調整が可能とする。
【００５８】
次に、ステップＳ２０８では、ステップＳ２０７により推定した外部負荷トルクを、式（７）に適用することで燃焼によるガス圧トルクを推定し、ガス圧によりピストン２を押し下げる力のベクトルと、クランク角度に対してトルクがかかるベクトル（クランクシャフト４の接線方向）の関係から筒内圧力を推定する。式（７）を推定対象であるＰｃｙｌ
_ｂについて変形すると式（１５）になる。
【００５９】
【数１２】
0006029726.tif 000013
【００６０】
次に、ステップＳ２０９ではステップＳ２０８で得られた燃焼気筒の筒内圧力Ｐｃｙｌ_ｂを用いて燃焼パラメータを推定する。クランク角度毎に得られた筒内圧力Ｐｃｙｌ_ｂ(θ)より、式（１６）を用いてクランク角度に対するシリンダの内部の気体に与えられた熱量を示す熱発生率ｄＱ（θ）を推定する。
【００６１】
【数１３】
0006029726.tif 000014
【００６２】
ここでκは燃料と空気の混合気の比熱比であり、具体的には１．３から１．４の値である。
【００６３】
次に、燃焼パラメータである質量燃焼割合ＭＦＢ(θ)を求める。ＭＦＢ(θ)は、熱発生率の積算値の割合であり、ＭＦＢ(θ)を推定するための積算領域において、熱発生率における瞬時積算値Ｑ(θ)を全積算値Ｑ０で割った値である。よって、ＭＦＢ(θ)は、次式（１７）で表される。
【００６４】
【数１４】
0006029726.tif 000015
【００６５】
図５は外部負荷トルク推定における近似処理別の推定筒内圧力波形図であり、図６は燃焼パラメータ推定における近似処理別の熱発生率波形図である。
破線で示す近似処理ＡはＳ２０７において線形近似により外部負荷トルクを推定した場合であり、実線で示す近似処理Ｂは本実施の形態における推定処理である。上死点近傍において近似処理Ａでは筒内圧波形に段差が生じている。このように上死点において推定したガス圧トルクが零とならない場合、筒内圧にスパイクや段差が生じ、燃焼パラメータ推定にも影響を与える。それに対し本実施の形態における外部負荷トルク推定によれば、筒内圧力波形に不自然な段差やスパイクが解消される。
【００６６】
前記において説明した実施の形態１では、３気筒エンジンにおける推定処理を説明したが、４気筒エンジンについてもほぼ同様の処理で外部負荷トルクを推定することが可能である。４気筒では各気筒の上下死点における位置が重なるため、式（８）は式（１８）と
なる。
【００６７】
【数１５】
0006029726.tif 000016
【００６８】
しかしながらそれぞれの気筒で燃焼により熱量が発生する期間が重なる場合は、筒内圧及び燃焼パラメータを推定する際に燃焼が重なっていない領域で推定した筒内圧の結果から圧力の変化を予測する、もしくは燃焼が重なる領域を除外して筒内圧力推定を行う必要がある。
【００６９】
燃焼パラメータの推定結果からＭＦＢ(θ)が５０％となるクランク角度を推定し、その位置が最適になるように点火時期を制御したり、ＭＦＢ(θ)が１０％から８０％となる期間を推定し、その値を基にＥＧＲ量を調整することで燃費性能やエミッション性能を向上させることができる。
【００７０】
また、これにより気筒別での筒内圧力推定から上記の制御を行うことで、各気筒が最良の燃焼を維持することが出来る。
【００７１】
以上のように、実施の形態１による内燃機関の制御装置によれば、精度良く外部負荷トルクを推定することが出来、筒内圧力推定において上死点近傍でも滑らかな波形となるため、燃焼パラメータ推定における精度の向上とばらつきの低減につながる。
【００７２】
なお、前記においてはこの発明の実施の形態１に係る内燃機関の制御装置について説明したが、この発明は、その発明の範囲において、実施の形態を適宜、変形、省略することが可能である。
【符号の説明】
【００７３】
１  内燃機関、２  ピストン、３  コンロッド、４  クランクシャフト、５  クランク角度センサ、６、９  シグナルロータ、７  ＥＣＵ、８  カムシャフト、１０  カム角度センサ、１１  チェーン、１２  フライホイール、１３  第２クランク角度センサ、１４  インマニ圧センサ、１５  大気圧センサ、１６  クランク角度検出手段、１７  クランク角度補正手段、１８  角速度及び角加速度推定手段、１９  外部負荷トルク推定手段、２０  筒内圧力及び燃焼パラメータ推定手段
【要約】
【課題】複数気筒を有する内燃機関に於いて、気筒ごとの燃焼状態を精度よく推定することのできる内燃機関の制御装置を提供する。
【解決手段】上死点においてピストンにかかる力によるトルクが発生しないことを利用して外部負荷トルクを推定する。上死点が標本点上に存在しない場合でも、上死点近傍における複数の標本点での演算を行い、線形近似や多項式近似を利用して上死点上の値を求めることにより外部負荷トルクの誤差を低減させ、高精度な筒内圧及び燃焼状態を示すパラメータを推定する。
【選択図】図２【図１】
0006029726.tif 000017
【図２】
0006029726.tif 000018
【図３】
0006029726.tif 000019
【図４】
0006029726.tif 000020
【図５】
0006029726.tif 000021
【図６】
0006029726.tif 000022