(57)【要約】  内燃機関の制御に必要な機関パラメータ、特に、筒内圧センサの検出結果に基づく燃焼パラメータをリアルタイムに精度良く演算でき、それにより内燃機関の制御性を向上させることができる内燃機関の制御装置を提供する。本発明の内燃機関の制御装置は、気筒（３ａ）内の圧力（筒内圧（ＰＣＹＬ））を検出する筒内圧センサ（２１）と、電子制御ユニット（２）に設けられ、検出された筒内圧（ＰＣＹＬ）を用いて熱発生率（ｄＱｄθ）を演算する燃焼モデルを有し、熱発生率（ｄＱｄθ）を含む、内燃機関（３）の状態を表す機関パラメータ（熱発生率（ｄＱｄθ）、インマニ圧（Ｐｉｎ）、ＥＧＲ温度（Ｔｅｇｒ）、ＥＧＲ圧（Ｐｅｇｒ））を演算するプラントモデル（４２）と、電子制御ユニット（２）に設けられ、プラントモデル（４２）で演算された機関パラメータを用いて、内燃機関（３）を制御するエンジン制御部（４３）と、を備える。
2017094349.tif 000002
【特許請求の範囲】
【請求項１】
気筒内の圧力を検出する筒内圧センサと、  電子制御ユニットに設けられ、前記筒内圧センサの検出結果を用いて前記気筒内の燃焼状態を表す燃焼パラメータを演算する燃焼モデルを有し、前記燃焼パラメータを含む、前記内燃機関の状態を表す機関パラメータを演算するプラントモデルと、  前記電子制御ユニットに設けられ、前記プラントモデルで演算された機関パラメータを用いて、前記内燃機関を制御する制御器と、を備えることを特徴とする内燃機関の制御装置。
【請求項２】
前記燃焼パラメータは熱発生率であり、前記燃焼モデルは、熱発生率の近似関数であるＷiebe関数を複数の１次関数によって近似した１次関数モデル式を用いて、熱発生率を演算することを特徴とする、請求項１に記載の内燃機関の制御装置。
【請求項３】
前記１次関数モデル式は複数のモデルパラメータを有し、  前記燃焼モデルは、前記筒内圧センサの検出結果に基づいて、前記複数のモデルパラメータをリアルタイムに同定する同定手段を有することを特徴とする、請求項２に記載の内燃機関の制御装置。
【請求項４】
前記電子制御ユニットは複数のプロセッサコアを有し、前記筒内圧センサの検出結果を用いて燃焼演算を実行する燃焼演算器、前記プラントモデル及び制御器が、前記複数のプロセッサコアにそれぞれ別個に搭載されていることを特徴とする、請求項１ないし３のいずれかに記載の内燃機関の制御装置。
【請求項５】
前記内燃機関は、前記気筒内に燃料を直接、噴射する燃料噴射弁を有し、  前記筒内圧センサは、前記燃料噴射弁と一体に設けられていることを特徴とする、請求項１ないし４のいずれかに記載の内燃機関の制御装置。
【発明の詳細な説明】
【技術分野】
【０００１】  本発明は、プラントモデルを用いて内燃機関の状態を表す機関パラメータを演算し、その演算結果に応じて内燃機関を制御する内燃機関の制御装置に関する。
【背景技術】
【０００２】  近年、内燃機関の排ガスや燃費に関する規制が厳しくなるのに応じて、この規制をクリアするために設けられるセンサやデバイスの数や種類が増加し、それに伴ってＥＣＵ（電子制御ユニット）への検出信号の入力数が増加する傾向にあり、コスト上昇の原因になっている。このため、低コスト化と制御性の向上を図るために、内燃機関の状態をセンサで検出するのに代えて、演算によって推定するためのプラントモデル（バーチャルセンサ）が開発されている。【０００３】  例えば特許文献１に開示された内燃機関の制御装置では、ＥＧＲ制御のために、ＥＧＲ率の複数の目標値候補が設定され、複数の目標値候補のそれぞれについて、ＥＧＲ弁開度の将来値がプラントモデルを用いて予測されるとともに、プラントモデルの複数のモデルパラメータが個々に学習される。また、制御装置は、多数のプロセッサコアを搭載したマルチコアプロセッサを備えており、各目標値候補に対する将来値の予測タスクと各モデルパラメータに対する学習タスクが、別々のコアに割り当てられる。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１３−２２８８５９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】  上述したように、従来の制御装置では、将来値の予測タスクが複数の目標値候補ごとに、学習タスクが複数のモデルパラメータごとに、プロセッサコアにそれぞれ割り当てられ、将来値の予測とモデルパラメータの学習に多数のプロセッサコアが使用される。このようなプロセッサコアの使用は、車両に通常、搭載されている、限られたマイコン性能に対して現実的でないとともに、内燃機関の制御に必要な、内燃機関の状態を表す他のパラメータの演算に支障を及ぼすおそれがある。【０００６】  例えば、内燃機関の気筒内の圧力を筒内圧センサで検出し、その検出結果から、燃焼によって発生する気筒内の圧力、熱やエネルギなどの燃焼状態を表す燃焼パラメータを取得することが知られている。このような燃焼パラメータは、気筒内での実際の燃焼状態を良好に反映するため、内燃機関の制御に非常に有効である。一方、有効な燃焼パラメータを得るためには、筒内圧センサの検出結果を、燃焼サイクルごとに逐次、解析することが望ましく、その場合の演算負荷は非常に大きい。上述した構成の従来の制御装置では、演算性能が制限されるため、燃焼パラメータの演算を良好に行えないおそれがある。【０００７】  本発明は、以上のような課題を解決するためになされたものであり、内燃機関の制御に必要な機関パラメータ、特に、筒内圧センサの検出結果に基づく燃焼パラメータをリアルタイムに精度良く演算でき、それにより内燃機関の制御性を向上させることができる内燃機関の制御装置を提供することを目的とする。
【課題を解決するための手段】
【０００８】  上記の目的を達成するため、請求項１に係る内燃機関の制御装置は、気筒３ａ内の圧力を検出する筒内圧センサ２１と、電子制御ユニット２に設けられ、筒内圧センサ２１の検出結果（実施形態における（以下、本項において同じ）筒内圧ＰＣＹＬ）を用いて気筒３ａ内の燃焼状態を表す燃焼パラメータ（熱発生率ｄＱｄθ）を演算する燃焼モデルを有し、燃焼パラメータを含む、内燃機関３の状態を表す機関パラメータ（熱発生率ｄＱｄθ、インマニ圧Ｐｉｎ、ＥＧＲ温度Ｔｅｇｒ、ＥＧＲ圧Ｐｅｇｒ）を演算するプラントモデル（モデル演算部４２）と、電子制御ユニット２に設けられ、プラントモデルで演算された機関パラメータを用いて、内燃機関３を制御する制御器（エンジン制御部４３）と、を備えることを特徴とする。【０００９】  この構成によれば、プラントモデルに含まれる燃焼モデルが、筒内圧センサの検出結果を用いて燃焼パラメータを演算するので、気筒内に発生する実際の圧力を反映させながら、燃焼パラメータを精度良く演算できる。また、燃焼モデルと、燃焼パラメータを用いて内燃機関を制御する制御器が、１つの電子制御ユニット内に設けられているので、制御器は、燃焼モデルで演算された燃焼パラメータを、通信ディレイがない状態でリアルタイムに用いることができる。以上から、燃焼パラメータを用いた内燃機関の制御性を向上させることができる。さらに、燃焼モデル以外のプラントモデルが燃焼パラメータ以外の機関パラメータを演算するので、その検出のために設けられていたセンサを省略でき、低コスト化を図ることができる。【００１０】  請求項２に係る発明は、請求項１に記載の内燃機関の制御装置において、燃焼パラメータは熱発生率ｄＱｄθであり、燃焼モデルは、熱発生率ｄＱｄθの近似関数であるＷiebe関数を複数の１次関数によって近似した１次関数モデル式（図５、図６）を用いて、熱発生率ｄＱｄθを演算すること（図９のステップ１２、１３）を特徴とする。【００１１】  この構成によれば、燃焼パラメータとしての熱発生率が、Ｗiebe関数を複数の１次関数で近似した１次関数モデル式を用いて演算される。Ｗiebe関数は、熱発生率の近似関数として知られており、全体的な形状が比較的単純で、直線に近い部分を多く有する。このため、Ｗiebe関数を複数の１次関数で精度良く近似できる。また、複数の１次関数で構成される１次関数モデル式の演算負荷は、Ｗiebe関数と比べて非常に小さい。したがって、熱発生率の演算を、その精度を維持しながら、短時間で応答良く行うことができ、熱発生率を用いた内燃機関の制御性をさらに向上させることができる。【００１２】  請求項３に係る発明は、請求項２に記載の内燃機関の制御装置において、１次関数モデル式は複数のモデルパラメータ（第１〜第４モデル基準点ＰＭ１〜ＰＭ４）を有し、燃焼モデルは、筒内圧センサ２１の検出結果に基づいて、複数のモデルパラメータをリアルタイムに同定する同定手段（モデル演算部４２、図１５）を有することを特徴とする。【００１３】  この構成によれば、１次関数モデル式の複数のモデルパラメータが、筒内圧センサの検出結果にに基づいてリアルタイムに同定されるので、燃焼状態のばらつきや経年劣化などによる１次関数モデル式のモデル誤差を随時、適切に補償でき、熱発生率の演算精度を良好に維持することができる。【００１４】  請求項４に係る発明は、請求項１ないし３のいずれかに記載の内燃機関の制御装置において、電子制御ユニット２は複数のプロセッサコア４１〜４３を有し、筒内圧センサ２１の検出結果を用いて燃焼演算を実行する燃焼演算器（ＣＰＳ演算部４１）、プラントモデル（モデル演算部４２）及び制御器（エンジン制御部４３）が、複数のプロセッサコア４１〜４３にそれぞれ別個に搭載されていることを特徴とする。【００１５】  この構成によれば、筒内圧センサの検出結果を用いて燃焼演算を実行する燃焼演算器と、プラントモデル及び制御器が、電子制御ユニットの複数のプロセッサコアにそれぞれ別個に搭載されている。これにより、燃焼演算器による燃焼演算、プラントモデルによる機関パラメータの演算と、制御器による内燃機関の制御を、それぞれ高い演算速度又は制御速度で行うとともに、互いの間のデータの授受も応答良く行うことができるので、内燃機関の制御性をさらに向上させることができる。【００１６】  請求項５に係る発明は、請求項１ないし４のいずれかに記載の内燃機関の制御装置において、内燃機関３は、気筒３ａ内に燃料を直接、噴射する燃料噴射弁４を有し、筒内圧センサ２１は、燃料噴射弁４と一体に設けられていることを特徴とする。【００１７】  この構成によれば、筒内圧センサが燃料噴射弁と一体に設けられているので、燃料噴射弁や点火プラグなどの機器とシリンダヘッドとの間に座金状の検出部を配置したタイプの筒内圧センサと比較して、シリンダヘッドの振動の影響を抑制しながら、筒内圧をより精度良く検出できる。これにより、筒内圧センサの検出結果を用いた燃焼パラメータの演算精度がさらに高められることで、内燃機関の制御性をさらに向上させることができる。
【図面の簡単な説明】
【００１８】
【図１】本発明を適用した内燃機関を概略的に示す図である。
【図２】内燃機関の制御装置を示すブロック図である。
【図３】図２の制御装置を詳細に示す図である。
【図４】エア系モデルの概念を示す図である。
【図５】熱発生率を算出するための燃焼モデルを示す図である。
【図６】燃焼モデルの設定方法を説明するための図である。
【図７】モデル演算処理を示すフローチャートである。
【図８】インマニモデルを気体のパラメータの入出力関係と併せて示す図である。
【図９】筒内温度の算出処理を示すフローチャートである。
【図１０】エキマニモデルを気体のパラメータの入出力関係と併せて示す図である。
【図１１】ＥＧＲ弁の上流側のＥＧＲ通路モデルを気体のパラメータの入出力関係と併せて示す図である。
【図１２】ＥＧＲ制御処理を示すフローチャートである。
【図１３】補正基準点を算出するための燃焼演算処理を示すフローチャートである。
【図１４】補正基準点の算出方法を説明するための図である。
【図１５】モデル基準点の同定処理を示すフローチャートである。
【図１６】筒内圧センサの故障判定処理を示すフローチャートである。
【図１７】燃料噴射弁及びそれと一体の筒内圧センサの外観を示す図である。
【発明を実施するための形態】
【００１９】  以下、図面を参照しながら、本発明の好ましい実施形態を詳細に説明する。図１は、本発明を適用した内燃機関（以下「エンジン」という）３を示す。このエンジン３は、車両（図示せず）に搭載された、例えば４気筒のガソリンエンジンである。各気筒３ａ（１つのみ図示）のピストン３ｂとシリンダヘッド３ｃとの間には、燃焼室３ｄが形成されている。
【００２０】  各気筒３ａには、吸気コレクタ部６ａを有する吸気マニホルド６ｂを介して、吸気通路６が接続され、吸気弁８が設けられるとともに、排気コレクタ部７ａを有する排気マニホルド７ｂを介して、排気通路７が接続され、排気弁９が設けられている。また、各気筒３ａには、燃焼室３ｄに臨むように、燃料噴射弁４及び点火プラグ５（図２参照）が設けられている。燃料噴射弁４による燃料の噴射量・噴射時期、及び点火プラグ５の点火時期は、後述する電子制御ユニット（以下「ＥＣＵ」という）２からの制御信号によって制御される。
【００２１】  さらに、各気筒３ａには、その内部の圧力である筒内圧ＰＣＹＬを検出するための筒内圧センサ２１が設けられている（図２参照）。図１７に示すように、筒内圧センサ２１は、燃料噴射弁４と一体型のものであり、燃料噴射弁４の先端部に取り付けられたリング状の圧力検出素子２１ａと、増幅回路ユニット（図示せず）などを備えている。圧力検出素子２１ａは、筒内圧ＰＣＹＬの変化率を検出し、増幅回路ユニットは、圧力検出素子２１ａの検出信号をフィルタリング及び増幅するとともに、筒内圧ＰＣＹＬに変換した後、その検出信号をＥＣＵ２に出力する。このように、筒内圧センサ２１が燃料噴射弁４の先端部に一体に設けられているので、一般的な座金タイプの筒内圧センサと比較して、シリンダヘッド３ｃの振動の影響を抑制しながら、筒内圧ＰＣＹＬをより精度良く検出することができる。
【００２２】  吸気通路６の吸気コレクタ部６ａよりも上流側には、スロットル弁機構１０が設けられている。このスロットル弁機構１０は、吸気通路６内に配置されたバタフライ式のスロットル弁１０ａと、スロットル弁１０ａを駆動するＴＨアクチュエータ１０ｂを有する。スロットル弁１０ａの開度（以下「スロットル弁開度」という）θＴＨは、ＴＨアクチュエータ１０ｂに供給される電流をＥＣＵ２で制御することによって制御され、それにより、燃焼室３ｄに供給される新気量が調整される。
【００２３】  また、エンジン３には、燃焼室３ｄから排気通路７に排出された排ガスの一部を、ＥＧＲガスとして、吸気通路６に還流させるためのＥＧＲ装置１１が設けられている。ＥＧＲ装置１１は、ＥＧＲ通路１２と、ＥＧＲ通路１２の途中に設けられたＥＧＲ弁機構１３及びＥＧＲクーラ１４などで構成されている。ＥＧＲ通路１２は、排気通路７の排気コレクタ部７ａと吸気通路６の吸気コレクタ部６ａに接続されている。
【００２４】  ＥＧＲ弁機構１３は、ＥＧＲ通路１２内に配置されたポペット式のＥＧＲ弁１３ａと、ＥＧＲ弁１３ａを駆動するＥＧＲアクチュエータ１３ｂを有する。ＥＧＲ弁１３ａのリフト量（以下「ＥＧＲ弁開度」という）ＬＥＧＲは、ＥＧＲアクチュエータ１３ｂに供給される電流をＥＣＵ２で制御することによって制御され、それにより、吸気通路６に還流するＥＧＲ量が調整される。
【００２５】  エンジン３のクランクシャフトには、クランク角センサ２２が設けられている（図２参照）。クランク角センサ２２は、クランクシャフトの回転に伴い、パルス信号であるＣＲＫ信号及びＴＤＣ信号をＥＣＵ２に出力する。ＣＲＫ信号は、所定のクランク角度（例えば１度）ごとに出力される。ＥＣＵ２は、このＣＲＫ信号に基づき、エンジン３の回転数（以下「エンジン回転数」という）ＮＥを算出する。
【００２６】  ＴＤＣ信号は、いずれかの気筒３ａにおいてエンジン３のピストン３ｂが吸気行程の開始時のＴＤＣ位置にあることを表す信号であり、本実施形態のようにエンジン３が４気筒の場合には、クランク角度１８０度ごとに出力される。ＥＣＵ２は、ＴＤＣ信号およびＣＲＫ信号に応じて、ＴＤＣ信号の出力タイミングを基準とするクランク角θを、気筒３ａごとに算出する。
【００２７】  また、吸気通路６のスロットル１０ａよりも上流側には、大気圧センサ２３及び外気温センサ２４が設けられている。大気圧センサ２３は大気圧ＰＡを検出し、外気温センサ２４は吸気通路６に導入される外気（新気）の温度ＴＡを検出し、それらの検出信号をＥＣＵ２に出力する。
【００２８】  さらに、ＥＣＵ２には、スロットル弁開度センサ２５からスロットル弁開度θＴＨを表す検出信号が入力され、ＥＧＲ弁開度センサ２６からＥＧＲ弁開度ＬＥＧＲを表す検出信号が入力される。
【００２９】  図３に示すように、ＥＣＵ２は、入出力部３１と、マルチコア型の演算処理ユニット（以下「ＭＣＵ」という）３２を備えている。入出力部３１は、上述した各種のセンサ２１〜２６から検出信号が入力されるとともに、燃料噴射弁４、点火プラグ５及びＥＧＲアクチュエータ１３ｂなどへ駆動信号を出力するものである。
【００３０】  ＭＣＵ３２は、第１〜第３のプロセッサコア４１〜４３と、プロセッサコア４１〜４３にそれぞれ対応して設けられたキャッシュメモリ４４〜４６と、プロセッサコア４１〜４３によって共通に使用される共有メモリ４７を有する。キャッシュメモリ４４〜４６、共有メモリ４７及び入出力部３１は、バス５０を介して互いに接続されている。入出力部３１に入力されたデータは、まず共有メモリ４７に格納される。プロセッサコア４１〜４３は、演算処理に必要なデータを共有メモリ４７から読み出して、キャッシュメモリ４４〜４６に一時的に格納し、演算処理を実行する。
【００３１】  より具体的には、第１のプロセッサコア４１（以下「ＣＰＳ演算部４１」という）は、筒内圧センサ２１で検出された筒内圧ＰＣＹＬとクランク角θに基づき、気筒３ａ内の燃焼状態を表す、熱発生率ｄＱｄθなどの燃焼パラメータを演算する燃焼演算処理を実行する。
【００３２】  第２のプロセッサコア４２（以下「モデル演算部４２」という）は、後述するプラントモデルに基づき、エンジン３の状態を表すエンジンパラメータを演算するモデル演算処理を実行する。このエンジンパラメータには、吸気通路６、排気通路７及びＥＧＲ通路１２を吸気、排ガス及びＥＧＲガスのそれぞれの質量流量、温度や圧力などが含まれる。
【００３３】  また、第３のプロセッサコア４３（以下「エンジン制御部４３」という）は、モデル演算部４２で演算されたエンジンパラメータを用い、エンジン３の燃料噴射弁４、点火プラグ５、スロットル弁１０ａやＥＧＲ弁１３ａなどのデバイスを制御するための制御パラメータを演算するエンジン制御処理を実行する。演算された制御パラメータは、入出力部３１に送られ、入出力部３１で駆動信号に変換された後、デバイスに出力される。
【００３４】  なお、本実施形態では、ＣＰＳ演算部４１が燃焼演算器に、モデル演算部４２がプラントモデル及び同定手段に、エンジン制御部４３が制御器に、それぞれ相当する。
【００３５】  上述したモデル演算処理の基礎となるプラントモデルは、エア系モデルと燃焼モデルに分類される。図４に示すように、エア系モデルは、吸気、排ガスやＥＧＲが流れるエンジン３の流路（吸気通路６、排気通路７及びＥＧＲ通路１２など）の構成を、スロットル弁１０ａやＥＧＲ弁１３ａなどが存在する「オリフィス」の部分とそれ以外の「レシーバ」の部分との組み合わせとして、モデル化したものである。また、レシーバに連続の式（質量保存則、エネルギ保存則）や気体の状態方程式などを適用し、オリフィスにオリフィスの式を適用することによって、エンジン３の流路の各部位における流体の質量流量、温度や圧力などが算出される。
【００３６】  一方、燃焼モデルは、熱発生率の近似関数として一般に知られているＷiebe関数を、演算負荷の軽減のために簡略化し、モデル化したものである。より具体的には、図５及び図６に示すように、燃焼モデルは、Ｗiebe関数（点線）を熱発生率の発生パターンに応じて４つの期間（第１気化期間ｅｈ１、第２気化期間ｅｈ２、第１燃焼期間ｂｈ１、第２燃焼期間ｂｈ２）に区分し、これらの４つの期間をそれぞれ第１〜第４の１次関数Ｉ〜IVで近似したものである。また、１次関数Ｉ〜IVを設定するために、次の４つのモデル基準点ＰＭ１〜ＰＭ４が用いられる。
【００３７】  第１モデル基準点ＰＭ１は、燃焼の開始直前に熱発生率ｄＱｄθが最小になる点に相当し、そのときの最小熱発生率ｄＱｄθｍｉｎとそれに対応するクランク角θｍｉｎから、ＰＭ１＝（θｍｉｎ，ｄＱｄθｍｉｎ）で表される。第２モデル基準点ＰＭ２は、熱発生率ｄＱｄθをクランク角θで微分した熱発生率微分値ｄＱｄ２θが最大になる点に相当し、そのときの熱発生率である最大微分値対応熱発生率ｄＱｄθｍａｘ２とそれに対応するクランク角θｍａｘ２から、ＰＭ２＝（θｍａｘ２，ｄＱｄθｍａｘ２）で表される。
【００３８】  第３モデル基準点ＰＭ３は、熱発生率ｄＱｄθが最大になる点に相当し、そのときの最大熱発生率ｄＱｄθｍａｘとそれに対応するクランク角θｍａｘから、ＰＭ３＝（θｍａｘ，ｄＱｄθｍａｘ）で表される。また、第４モデル基準点ＰＭ４は、熱発生率微分値ｄＱｄ２θが最小になる点に相当し、そのときの熱発生率である最小微分値対応熱発生率ｄＱｄθｍｉｎ２とそれに対応するクランク角θｍｉｎ２から、ＰＭ４＝（θｍｉｎ２，ｄＱｄθｍｉｎ２）で表される。
【００３９】  以上の４つのモデル基準点ＰＭ１〜ＰＭ４が定まると、それに基づき、１次関数Ｉ〜IVが以下のようにして設定される。まず、図６（ａ）に示すように、第３の１次関数III は、第２モデル基準点ＰＭ２と第３モデル基準点ＰＭ３を通る直線（１次式）として、一義的に設定される。
【００４０】  具体的には、この１次式をｄＱｄθ＝Ａ・θ＋Ｂ（Ａ：傾き、Ｂ：切片）とすると、傾きＡは、Ａ＝（ｄＱｄθｍａｘ−ｄＱｄθｍａｘ２）／（θｍａｘ−θｍａｘ２）で算出され、切片Ｂは、Ｂ＝ｄＱｄθ−Ａ・θに、ｄＱｄθｍａｘ、θｍａｘと、算出したＡを代入することによって算出される。このような算出手法は、後述する他の１次関数を設定する場合にも同様に適用される。また、第３の１次関数III において熱発生率ｄＱｄθ＝０になるクランク角θが、燃焼開始角θｂｓとして算出され、燃焼開始点Ｐｂｓ（θｂｓ，０）が設定される。
【００４１】  また、図６（ｂ）に示すように、第４の１次関数IVは、第３モデル基準点ＰＭ３と第４モデル基準点ＰＭ４を通る直線として設定される。また、設定された第４の１次関数IVにおいて熱発生率ｄＱｄθ＝０になるクランク角θが、燃焼終了角θｂｅとして算出され、燃焼終了点Ｐｂｅ（θｂｅ，０）が設定される。さらに、図６（ｃ）に示すように、第３モデル基準点ＰＭ３と燃焼開始点Ｐｂｓとの間のクランク角θの差（＝θｍａｘ−θｂｓ）が、第１燃焼期間ｂｈ１として算出され、燃焼終了点Ｐｂｅと第３モデル基準点ＰＭ３との間のクランク角θの差（＝θｂｅ−θｍａｘ）が、第２燃焼期間ｂｈ２として算出される。
【００４２】  図６（ｄ）に示すように、第１の１次関数Ｉは、気化開始点Ｐｅｓ（θｅｓ，０）と第１モデル基準点ＰＭ１を通る直線として設定される。この気化開始点Ｐｅｓは、燃焼前に混合気が気化し始める点であり、気化開始角θｅｓは所定の一定値に設定されている。また、第２の１次関数IIは、第１モデル基準点ＰＭ１と燃焼開始点Ｐｂｓを通る直線として設定される。さらに、第１モデル基準点ＰＭ１と気化開始点Ｐｅｓとの間のクランク角θの差（＝θｍｉｎ−θｅｓ）が、第１気化期間ｅｈ１として算出され、燃焼開始点Ｐｂｓと第１モデル基準点ＰＭ１との間のクランク角θの差（＝θｂｓ−θｍｉｎ）が、第２気化期間ｅｈ２として算出される。
【００４３】  以上のように、燃焼モデルが第１〜第４の１次関数Ｉ〜IVによって簡略化して設定される結果、この燃焼モデルを用いて熱発生率ｄＱｄθを算出する際の演算負荷は、Ｗiebe関数を用いる場合と比較して大幅に軽減される。
【００４４】  次に、図７を参照しながら、モデル演算部４２で実行されるモデル演算処理について説明する。この処理は、前述したプラントモデルに基づき、ＥＧＲ制御に必要なインマニ圧Ｐｉｎと、ＥＧＲ弁１３のすぐ上流側における排ガスの温度及び圧力を、それぞれＥＧＲ温度Ｔｅｇｒ及びＥＧＲ圧Ｐｅｇｒとして推定するものである。本処理は、ＣＲＫ信号の発生に同期して気筒３ａごとに実行される。
【００４５】  本処理では、まずステップ１（「Ｓ１」と図示。以下同じ」）において、吸気通路６のスロットル弁１０ａよりも下流側の吸気圧であるインマニ圧Ｐｉｎを算出する。このインマニ圧Ｐｉｎの算出は、図８に示すように、図１の吸気通路６のスロットル弁１３ａよりも下流側から吸気チャンバ６ａまでの部分とＥＧＲ通路１２の接続部を、インマニモデル（レシーバ）として設定するとともに、インマニモデルにおいて成立する、以下に述べるようなパラメータ間の関係に基づいて行われる。
【００４６】  具体的には、図８に示すように、レシーバに、ポートＰＯ１を介して流入する新気の質量流量、温度、定圧比熱、定積比熱及びエネルギを、それぞれｍｄｏｔ１、Ｔ１、Ｃｐ１、Ｃｖ１及びＥ１とし、ポートＰＯ３を介して流入するＥＧＲガスの質量流量、温度、定圧比熱、定積比熱及びエネルギを、それぞれｍｄｏｔ３、Ｔ３、Ｃｐ３、Ｃｖ３及びＥ３とする。また、レシーバ内の気体の質量（新気分質量、ＥＧＲガス分質量）、温度、定圧比熱、定積比熱、圧力及びＥＧＲ率を、それぞれＭ（Ｍ１，Ｍ３）、Ｔ、Ｃｐ、Ｃｖ、Ｐ及びｒＰｏｒｔ３とし、レシーバからポートＰＯ２を介して流出する気体の質量流量及びエネルギをそれぞれｍｄｏｔ２及びＥ２とすると、これらのパラメータの間には、連続の式や気体の状態方程式などから、以下の式（１）〜（１６）の関係が成立する。なお、これらのパラメータは時間の関数であるが、式中では便宜上、その表示は省略されている。
【００４７】  まず、レシーバ内の気体に対する質量保存則から、次式（１）が成立し、ポートＰＯ１及びＰＯ３を介して流入する新気及びＥＧＲガスに対する質量保存則から、式（２）及び（３）がそれぞれ成立する。【数１】
2017094349.tif 000003
【数２】
2017094349.tif 000004
【数３】
2017094349.tif 000005
【００４８】  また、レシーバに流入する新気及びＥＧＲガスの定圧熱容量Ｍ１・Ｃｐ１、Ｍ３・Ｃｐ３がそれぞれ保存される関係から、式（４）及び（５）が成立し、新気及びＥＧＲガスの定積熱容量Ｍ１・Ｃｖ１、Ｍ３・Ｃｖ３がそれぞれ保存される関係から、式（６）及び（７）が成立する。【数４】
2017094349.tif 000006
【数５】
2017094349.tif 000007
【数６】
2017094349.tif 000008
【数７】
2017094349.tif 000009
【００４９】  また、レシーバに流入する新気のエネルギ（エンタルピー）Ｅ１、レシーバから流出する気体のエネルギＥ２、及びレシーバに流入するＥＧＲガスのエネルギＥ３は、それぞれ式（８）〜（１０）で表される。【数８】
2017094349.tif 000010
【数９】
2017094349.tif 000011
【数１０】
2017094349.tif 000012
レシーバ内の気体の定圧熱容量Ｃｐ・Ｍ及び定積熱容量Ｃｖ・Ｍはそれぞれ、新気分及びＥＧＲガス分の和であるという関係から、式（１１）及び（１２）が成立する。
【数１１】
2017094349.tif 000013
【数１２】
2017094349.tif 000014
【００５０】  レシーバから外部への放熱量Ｑｗａｌｌは、次式（１３）で表される。【数１３】
2017094349.tif 000015
ここで、Ｔｗａｌｌはレシーバの壁温度、Ｓｗａｌｌは壁面積（定数）、Ｋは熱伝達係数（定数）である。  また、レシーバ内の気体に対するエネルギ保存則から、式（１４）が成立する。
【数１４】
2017094349.tif 000016
さらに、レシーバ内に気体の状態方程式を適用することによって、式（１５）が成立する。
【数１５】
2017094349.tif 000017
また、レシーバ内の気体のＥＧＲ率ｒＰｏｒｔ３は、式（１６）で表される。
【数１６】
2017094349.tif 000018
【００５１】  以上の式（１）〜（１６）に例えば連立方程式を適用することにより、新気及びＥＧＲガスに関するパラメータを既知として、レシーバ内の気体に関するパラメータや、レシーバから流出する気体の質量流量ｍｄｏｔ２及びエネルギＥ２を算出することができ、前記ステップ１では、レシーバ内の圧力Ｐがインマニ圧Ｐｉｎとして算出される。
【００５２】  図７に戻り、前記ステップ１に続くステッ２では、筒内温度Ｔｃｙｌを算出する。この算出処理は、前述した燃焼モデルを設定するとともに、設定した燃焼モデルに基づいて筒内温度Ｔｃｙｌを算出するものであり、図９に示すサブルーチンに従って実行される。
【００５３】  図９の処理では、まずステップ１１において、燃焼モデルのモデル基準点ＰＭ１〜ＰＭ４を算出する。その算出は、エンジン３の運転状態、例えばエンジン回転数ＮＥ、混合気の空燃比、点火時期及びＥＧＲ率に応じ、所定のマップ（図示せず）から、モデル基準点ＰＭ１〜ＰＭ４のマップ値をそれぞれ検索するとともに、これらのマップ値を後述する補正項で補正することによって、行われる。また、上記のＥＧＲ率としては、例えば、インマニモデルにおいて前記式（１６）によって算出されたＥＧＲ率ｒＰｏｒｔ３が用いられる。
【００５４】  次に、算出したモデル基準点ＰＭ１〜ＰＭ４を用い、前述した手法によって、４つの１次関数Ｉ〜IVから成る燃焼モデルを設定する（ステップ１２）とともに、設定した燃焼モデルを用いて、熱発生率ｄＱｄθを算出する（ステップ１３）。次に、算出した熱発生率ｄＱｄθを用い、次式（１７）によって、燃焼モデルによる推定筒内圧Ｐｍを算出する（ステップ１４）。【数１７】
2017094349.tif 000019
ここで、筒内容積Ｖの変化量ｄＶは、クランク角θに応じて一義的に定まり、筒内圧変化量ｄＰｍは、２つの演算タイミング間の差分として求められる。比熱比κは定数である。
【００５５】  次に、気筒３ａ内に気体の状態方程式を適用し、推定筒内圧Ｐｍを用い、次式（１８）によって、筒内温度Ｔｃｙｌを算出し（ステップ１５）、本処理を終了する。【数１８】
2017094349.tif 000020
【００５６】  図７に戻り、上記ステップ２に続くステップ３では、排気マニホルド７ｂ内の温度であるエキマニ温度Ｔｅｘを算出する。このエキマニ温度Ｔｅｘの算出は、図１０に示すように、排気通路７の排気マニホルド７ｂから排気チャンバ７ａまでの部分とＥＧＲ通路１２の分岐部を、エキマニモデル（レシーバ）として設定するとともに、エキマニモデルにおいて成立する、以下に述べるようなパラメータ間の関係に基づいて行われる。
【００５７】  図８との比較から明らかなように、インマニモデルでは、入力ポートが２つ、出力ポートが１つ（２入力／１出力）であるのに対し、エキマニモデルでは、入力ポートが１つで出力ポートが２つ（１入力／２出力）であるため、インマニモデルの場合と部分的に異なる、以下のようなパラメータの関係が成立する。
【００５８】  まず、レシーバ内の気体に対する質量保存則に基づく関係式は、インマニモデルにおける前記式（１）に代えて、次式（１）’で表される。【数１９】
2017094349.tif 000021
また、レシーバに流入する燃焼ガスの定圧熱容量Ｍ１・Ｃｐ１及び定積熱容量Ｍ１・Ｃｖ１の保存に関する関係式は、次式（４）’及び（６）’によって表される。
【数２０】
2017094349.tif 000022
【数２１】
2017094349.tif 000023
【００５９】  レシーバに流入する燃焼ガスのエネルギＥ１、レシーバから排気通路７の下流側に流出する排ガスのエネルギＥ２、及びレシーバからＥＧＲ通路１２に流出するＥＧＲガスのエネルギＥ３は、前記式（８）〜（１０）で同様に表される。また、レシーバ内の気体の定圧熱容量Ｃｐ・Ｍ及び定積熱容量Ｃｖ・Ｍはそれぞれ、次式（１１）’及び（１２）’によって表され、レシーバからの放熱量Ｑｗａｌｌは、前記式（１３）で同様に表される。【数２２】
2017094349.tif 000024
【数２３】
2017094349.tif 000025
【００６０】  また、レシーバ内の気体に対するエネルギ保存則に基づく関係式は、次式（１４）’によって表され、レシーバ内の気体の状態方程式は、前記式（１５）で同様に表される。【数２４】
2017094349.tif 000026
【００６１】  以上の式（１）’（４）’（６）’（８）〜（１０）（１１）’（１２）’（１３）（１４）’及び（１５）に、連立方程式を適用することにより、レシーバに流入する燃焼ガスに関するパラメータを既知として、レシーバ内の気体に関するパラメータや、レシーバから排気通路７側に流出する排ガスの質量流量ｍｄｏｔ２及びエネルギＥ２と、ＥＧＲ通路１２側に流出するＥＧＲガスの質量流量ｍｄｏｔ３及びエネルギＥ３を算出することができ、前記ステップ３では、レシーバ内の温度Ｔがエキマニ温度Ｔｅｘとして算出される。
【００６２】  図７に戻り、前記ステップ３に続くステップ４及び５では、ＥＧＲ弁１３ａのすぐ上流側におけるＥＧＲガスの温度及び圧力を、ＥＧＲ温度Ｔｅｇｒ及びＥＧＲ圧Ｐｅｇｒとしてそれぞれ算出し、本処理を終了する。このＥＧＲ温度Ｔｅｇｒ及びＥＧＲ圧Ｐｅｇｒの算出は、図１１に示すように、ＥＧＲ通路１２の排気通路７からの分岐部からＥＧＲ弁１３ａすぐ上流側までの部分を、ＥＧＲ通路モデル（レシーバ）として設定するとともに、ＥＧＲ通路モデルにおいて成立する、以下に述べるようなパラメータの関係に基づいて行われる。
【００６３】  図１１に示すように、ＥＧＲ通路モデルでは、入力ポート及び出力ポートがいずれも１つ（１入力／１出力）であるため、前述したインマニモデル及びエキマニモデルの場合と部分的に異なる、以下のようなパラメータの関係が成立する。
【００６４】  レシーバ内の気体に対する質量保存則に基づく関係式は、次式（１）”で表される。【数２５】
2017094349.tif 000027
また、レシーバに流入するＥＧＲガスの定圧熱容量Ｍ１・Ｃｐ１及び定積熱容量Ｍ１・Ｃｖ１の保存に関する関係式は、次式（４）”及び（６）”によって表される。
【数２６】
2017094349.tif 000028
【数２７】
2017094349.tif 000029
【００６５】  レシーバに流入するＥＧＲガスのエネルギＥ１と、レシーバから流出するＥＧＲガスのエネルギＥ２は、前記式（８）及び（９）によって同様に表される。また、レシーバ内の気体の定圧熱容量Ｃｐ・Ｍ及び定積熱容量Ｃｖ・Ｍと、レシーバからの放熱量Ｑｗａｌｌは、前記式（１１）’（１２）’及び（１３）で同様に表される。
【００６６】  また、レシーバ内の気体に対するエネルギ保存則に基づく関係式は、次式（１４）”によって表され、レシーバ内の気体の状態方程式は、前記式（１５）で同様に表される。【数２８】
2017094349.tif 000030
【００６７】  以上の式（１）”（４）”（６）”（８）（９）（１１）’（１２）’（１３）（１４）”及び（１５）に、連立方程式を適用することにより、レシーバに流入するＥＧＲガスに関するパラメータを既知として、レシーバ内の気体に関するパラメータを算出することができ、前記ステップ４では、レシーバ内の温度ＴがＥＧＲ温度Ｔｅｇｒとして算出され、前記ステップ５では、レシーバ内の圧力ＰがＥＧＲ圧Ｐｅｇｒとして算出される。
【００６８】  次に、図１２を参照しながら、エンジン制御部４３で実行されるＥＧＲ制御処理について説明する。本処理は、ＴＤＣ信号の発生に同期して気筒３ａごとに実行される。本処理では、まずステップ２１において、目標ＥＧＲ量ＧＥＧＲＣＭＤを設定する。この設定は、例えば目標トルク及びエンジン回転数ＮＥに応じ、所定のマップ（図示せず）を検索することによって、行われる。
【００６９】  次に、図７のステップ１で算出したインマニ圧Ｐｉｎとステップ５で算出したＥＧＲ圧Ｐｅｇｒを用い、次式（１９）によって、圧力関数Ψを算出する（ステップ２２）。【数２９】
2017094349.tif 000031
【００７０】  次に、ＥＧＲ圧Ｐｅｇｒ及び圧力関数Ψと、図７のステップ４で算出したＥＧＲ温度Ｔｅｇｒを用い、次式（２０）によって、ＥＧＲ弁１３ａを通過するＥＧＲガスの質量流量（以下「実ＥＧＲ量」という）ＧＥＧＲＡＣＴを算出する（ステップ２３）。【数３０】
2017094349.tif 000032
この式（２０）は、ＥＧＲ弁１３ａにオリフィスの式を適用したものであり、Ｒは気体定数、Ｃｄは流量係数で、いずれも定数である。また、Ａは、ＥＧＲ弁１３ａの開口面積であり、ＥＧＲ弁開度ＬＥＧＲに基づいて算出される。
【００７１】  次に、次式（２１）によって、ＥＧＲ弁１３ａの開口面積Ａの目標値である目標開口面積ＡＣＭＤを設定する（ステップ２４）。【数３１】
2017094349.tif 000033
この式（２１）は、上記のオリフィスの式（２０）を開口面積Ａについて表すとともに、実ＥＧＲ量ＧＥＧＲＡＣＴを目標ＥＧＲ量ＧＥＧＲＣＭＤに、開口面積Ａを目標開口面積ＡＣＭＤに、それぞれ置き換えたものである。
【００７２】  次に、目標ＥＧＲ量ＧＥＧＲＣＭＤと実ＥＧＲ量ＧＥＧＲＡＣＴとの差を、ＥＧＲ量偏差ΔＧＥＧＲとして算出する（ステップ２５）とともに、このＥＧＲ量偏差ΔＧＥＧＲに応じて、フィードバック補正項ΔＡＦＢを算出する（ステップ２６）。そして、このフィードバック補正項ΔＡＦＢを目標開口面積ＡＣＭＤに加算することによって、これを補正する（ステップ２７）。
【００７３】  次に、補正された目標開口面積ＡＣＭＤに応じて、ＥＧＲ弁１３ａを駆動するＥＧＲアクチュエータ１３ｂの目標電流値ＩＣＭＤを設定する（ステップ２８）。また、目標開口面積ＡＣＭＤとＥＧＲ弁開度ＬＥＧＲに基づいて算出される実際の開口面積Ａとの差を、開口面積偏差ΔＡとして算出する（ステップ２９）とともに、この開口面積偏差ΔＡに応じて、フィードバック補正項ΔＩＦＢを算出する（ステップ３０）。そして、このフィードバック補正項ΔＩＦＢを目標電流値ＩＣＭＤに加算することによって、これを補正し（ステップ３１）、本処理を終了する。
【００７４】  次に、燃焼モデルの同定処理について説明する。この同定処理は、燃焼モデルのモデル基準点ＰＭ１〜ＰＭ４を、筒内圧センサ２１で検出された実際の筒内圧ＰＣＹＬに基づいてリアルタイムに同定（補正）するものであり、ＣＰＳ演算部４１で実行される燃焼演算処理と、その結果を用いてモデル演算部４２で実行される同定演算処理で構成される。
【００７５】  図１３に示す燃焼演算処理は、筒内圧ＰＣＹＬに基づき、燃焼モデルのモデル基準点ＰＭ１〜ＰＭ４を同定する際の基準となる補正基準点ＰＣ１〜ＰＣ４を算出するものであり、ＣＲＫ信号の発生に同期して気筒３ａごとに実行される。
【００７６】  本処理では、まずステップ４１において、筒内圧ＰＣＹＬ及びクランク角θに基づき、次式（２２）によって熱発生率ｄＱｄθを算出する。【数３２】
2017094349.tif 000034
これにより、例えば図１４（ａ）の筒内圧ＰＣＹＬを表す曲線から、同図（ｂ）に示す熱発生率ｄＱｄθを表す曲線が得られる。
【００７７】  次に、ステップ４２において、熱発生率ｄＱｄθをクランク角θで微分することによって、熱発生率微分値ｄＱｄ２θを算出する。これにより、図１４（ｃ）に示す熱発生率微分値ｄＱｄ２θを表す曲線が得られる。
【００７８】  次に、ステップ４３〜４６において、図１４（ｂ）〜（ｄ）に示すように、熱発生率ｄＱｄθ及び熱発生率微分値ｄＱｄ２θに基づき、モデル基準点ＰＭ１〜ＰＭ４に対応する補正基準点ＰＣ１〜ＰＣ４をそれぞれ算出し、本処理を終了する。
【００７９】  具体的には、ステップ４３では、ステップ４１で算出された熱発生率ｄＱｄθのうちの燃焼の開始直前に生じる最小値を、最小熱発生率ｄＱｄθｍｉｎａとして抽出するとともに、それに対応するクランク角θｍｉｎａとの組み合わせから成る点（θｍｉｎａ，ｄＱｄθｍｉｎａ）を、第１補正基準点ＰＣ１として設定する。
【００８０】  ステップ４４では、ステップ４２で算出された熱発生率微分値ｄＱｄ２θの最大値が得られるときの熱発生率ｄＱｄθを、最大微分値対応熱発生率ｄＱｄθｍａｘ２ａとして抽出するとともに、それに対応するクランク角θｍａｘ２ａとの組み合わせから成る点（θｍａｘ２ａ，ｄＱｄθｍａｘ２ａ）を、第２補正基準点ＰＣ２として設定する。
【００８１】  ステップ４５では、熱発生率ｄＱｄθの最大値を、最大熱発生率ｄＱｄθｍａｘａとして抽出するとともに、それに対応するクランク角θｍａｘａとの組み合わせから成る点（θｍａｘａ，ｄＱｄθｍａｘａ）を、第３補正基準点ＰＣ３として設定する。
【００８２】  また、ステップ４６では、熱発生率微分値ｄＱｄ２θの最小値が得られるときの熱発生率ｄＱｄθを、最小微分値対応熱発生率ｄＱｄθｍｉｎ２ａとして抽出するとともに、それに対応するクランク角θｍｉｎ２ａとの組み合わせから成る点（θｍｉｎ２ａ，ｄＱｄθｍｉｎ２ａ）を、第４補正基準点ＰＣ４として設定する。
【００８３】  次に、図１５を参照しながら、モデル演算部４２で実行される同定演算処理について説明する。本処理は、燃焼モデルのモデル基準点ＰＭ１〜ＰＭ４を、同じ燃焼サイクルで得られた補正基準点ＰＣ１〜ＰＣ４にそれぞれ近似するように同定（補正）するものである。本処理は、ＴＤＣ信号の発生に同期して気筒３ａごとに実行される。
【００８４】  本処理では、まずステップ５１において、第１補正基準点ＰＣ１のクランク角要素であるθｍｉｎａと第１モデル基準点ＰＭ１のクランク角要素であるθｍｉｎとの差（＝θｍｉｎａ−θｍｉｎ）に、クランク角補正用の所定の補正係数Ｋθを乗算することによって、第１モデル基準点ＰＭ１のクランク角補正項ΔθＣ１を算出する。
【００８５】  また、ステップ５２では、第１補正基準点ＰＣ１の熱発生率要素であるｄＱｄθｍｉｎａと第１モデル基準点ＰＭ１の熱発生率要素であるｄＱｄθｍｉｎとの差（＝ｄＱｄθｍｉｎａ−ｄＱｄθｍｉｎ）に、熱発生率補正用の所定の補正係数ＫｄＱを乗算することによって、第１モデル基準点ＰＭ１の熱発生率補正項ΔｄＱＣ１を算出する。
【００８６】  以下同様に、第２モデル基準点ＰＭ２については、ステップ５３において、第２補正基準点ＰＣ２のクランク角θｍａｘ２ａと第２モデル基準点ＰＭ２のクランク角θｍａｘ２との差（＝θｍａｘ２ａ−θｍａｘ２）に、補正係数Ｋθを乗算することによって、クランク角補正項ΔθＣ２を算出し、ステップ５４において、第２補正基準点ＰＣ２の熱発生率ｄＱｄθｍａｘ２ａと第２モデル基準点ＰＭ２の熱発生率ｄＱｄθｍａｘ２との差（＝ｄＱｄθｍａｘ２ａ−ｄＱｄθｍａｘ２）に、補正係数ＫｄＱを乗算することによって、熱発生率補正項ΔｄＱＣ２を算出する。
【００８７】  第３モデル基準点ＰＭ３については、ステップ５５において、第３補正基準点ＰＣ３のクランク角θｍａｘａと第３モデル基準点ＰＭ３のクランク角θｍａｘとの差（＝θｍａｘａ−θｍａｘ）に、補正係数Ｋθを乗算することによって、クランク角補正項ΔθＣ３を算出し、ステップ５６において、第３補正基準点ＰＣ３の熱発生率ｄＱｄθｍａｘａと第３モデル基準点ＰＭ３の熱発生率ｄＱｄθｍａｘとの差（＝ｄＱｄθｍａｘａ−ｄＱｄθｍａｘ）に、補正係数ＫｄＱを乗算することによって、熱発生率補正項ΔｄＱＣ３を算出する。
【００８８】  また、第４モデル基準点ＰＭ４については、ステップ５７において、第４補正基準点ＰＣ４のクランク角θｍｉｎ２ａと第４モデル基準点ＰＭ４のクランク角θｍｉｎ２との差（＝θｍｉｎ２ａ−θｍｉｎ２）に、補正係数Ｋθを乗算することによって、クランク角補正項ΔθＣ４を算出し、ステップ５８において、第４補正基準点ＰＣ４の熱発生率ｄＱｄθｍｉｎ２ａと第４モデル基準点ＰＭ４の熱発生率ｄＱｄθｍａｘとの差（＝ｄＱｄθｍｉｎ２ａ−ｄＱｄθｍｉｎ２）に、補正係数ＫｄＱを乗算することによって、熱発生率補正項ΔｄＱＣ４を算出し、本処理を終了する。
【００８９】  以上のように算出されたクランク角補正項ΔθＣ１〜ΔθＣ４及び熱発生率補正項ΔｄＱＣ１〜ΔｄＱＣ４は、次回の燃焼サイクルにおいて、エンジン３の運転状態に応じてマップ検索された第１〜第４モデル基準点ＰＭ１〜ＰＭ４の対応するクランク角要素及び熱発生率要素に加算され、それにより、第１〜第４モデル基準点ＰＭ１〜ＰＭ４がリアルタイムで同定（補正）される。
【００９０】  次に、図１６を参照しながら、モデル演算部４２で実行される故障判定処理について説明する。本処理は、第１〜第４モデル基準点ＰＭ１〜ＰＭ４と第１補正基準点ＰＣ１〜ＰＣ４との比較結果に基づいて、筒内圧センサ２１の故障の有無を判定するものである。本処理は、ＴＤＣ信号の発生に同期して気筒３ａごとに実行される。
【００９１】  本処理では、まずステップ６１において、第１〜第４補正基準点ＰＣ１〜ＰＣ４の各クランク角要素と、それに対応する第１〜第４モデル基準点ＰＭ１〜ＰＭ４のクランク角要素との差の絶対値を、それぞれクランク角偏差Δθ１〜Δθ４として算出する。次に、算出したクランク角偏差Δθ１〜Δθ４がいずれも、クランク角判定用の所定のしきい値θＲＥＦ以下であるか否かを判別する（ステップ６２）。この答えがＮＯで、クランク角偏差Δθ１〜Δθ４の少なくとも１つがしきい値θＲＥＦを超えているときには、筒内圧センサ２１に故障が発生していると判定し、故障フラグＦ＿ＣＹＬＮＧを「１」にセットし（ステップ６３）、本処理を終了する。
【００９２】  ステップ６２の答えがＹＥＳのときには、ステップ６４において、第１〜第４補正基準点ＰＣ１〜ＰＣ４の各熱発生率要素と、それに対応する第１〜第４モデル基準点ＰＭ１〜ＰＭ４の熱発生率要素との差の絶対値を、それぞれ熱発生率偏差ΔｄＱ１〜ΔｄＱ４として算出する。次に、算出した熱発生率偏差ΔｄＱ１〜ΔｄＱ４がいずれも、熱発生率判定用の所定のしきい値ｄＱＲＥＦ以下であるか否かを判別する（ステップ６５）。この答えがＮＯで、熱発生率偏差ΔｄＱ１〜ΔｄＱ４の少なくとも１つがしきい値ｄＱＲＥＦを超えているときには、筒内圧センサ２１に故障が発生していると判定し、前記ステップ６３に進み、故障フラグＦ＿ＣＹＬＮＧを「１」にセットし、本処理を終了する。
【００９３】  一方、前記ステップ６５の答えがＹＥＳのときには、筒内圧センサ２１に故障が発生していないと判定し、故障フラグＦ＿ＣＹＬＮＧを「０」にセットし（ステップ６６）、本処理を終了する。上記のように、筒内圧センサ２１が故障していると判定され、故障フラグＦ＿ＣＹＬＮＧが「１」にセットされた場合には、筒内圧センサ２１の検出結果に基づく図１３の燃焼演算と図１４の同定演算が禁止される。
【００９４】  以上のように、本実施形態によれば、筒内圧センサ２１の検出結果を用いて設定される、プラントモデルの燃焼モデルに基づき、モデル演算部４２によって熱発生率ｄＱｄθを演算するので、気筒３ａ内に発生する実際の圧力を反映させながら、熱発生率ｄＱｄθを精度良く演算することができる。
【００９５】  また、モデル演算部４２と、エンジン３を制御するエンジン制御部４３が、いずれもプロセッサコアで構成され、１つのＥＣＵ２内に設けられているので、エンジン制御部４３は、モデル演算部４２で演算された熱発生率ｄＱｄθを、通信ディレイがない状態でリアルタイムに用いることができる。以上から、熱発生率ｄＱｄθを用いたＥＧＲ制御の制御性を向上させることができる。
【００９６】  さらに、プラントモデルのエア系モデルに基づき、ＥＧＲ制御に必要なインマニ圧Ｐｉｎ、ＥＧＲ温度ＴｅｇｒやＥＧＲ圧Ｐｅｇｒが演算によって求められるので、その検出のために設けられていたセンサを省略でき、低コスト化を図ることができる。
【００９７】  また、Ｗiebe関数を複数の１次関数で近似した１次関数モデル式によって燃焼モデルを設定し、この燃焼モデルを用いて熱発生率ｄＱｄθを演算するので、この熱発生率ｄＱｄθの演算を、その精度を維持しながら、短時間で応答良く行うことができ、熱発生率ｄＱｄθを用いたＥＧＲ制御の制御性をさらに向上させることができる。
【００９８】  さらに、この燃焼モデルのモデルパラメータであるモデル基準点ＰＭ１〜ＰＭ４が、筒内圧センサ２１の検出結果に基づいて算出された補正基準点ＰＣ１〜ＰＣ４によってリアルタイムに同定されるので、燃焼状態のばらつきや経年劣化などによる燃焼モデルのモデル誤差を随時、適切に補償でき、熱発生率ｄＱｄθの演算精度を良好に維持することができる。
【００９９】  また、モデル基準点ＰＭ１〜ＰＭ４と補正基準点ＰＣ１〜ＰＣ４との比較結果に基づいて、筒内圧センサ２１の故障を判定するので、筒内圧センサ２１の故障の判定を、燃焼モデルの設定と同定に用いられるパラメータを利用しながら、効率良く適切に行うことができる。
【０１００】  また、筒内圧センサ２１の検出結果を用いて補正基準点ＰＣ１〜ＰＣ４を演算するＣＰＳ演算部４１、モデル演算部４２及びエンジン制御部４３が、ＥＣＵ２のプロセッサコアにそれぞれ別個に搭載されているので、ＣＰＳ演算部４１による補正基準点ＰＣ１〜ＰＣ４の演算、モデル演算部４２による熱発生率ｄＱｄθ及び他のエンジンパラメータの演算と、エンジン制御部４３によるエンジン３の制御を、それぞれ高い演算速度又は制御速度で行うとともに、互いの間のデータの授受も応答良く行うことができるので、エンジン３の制御性をさらに向上させることができる。
【０１０１】  さらに、筒内圧センサ２１が燃料噴射弁４の先端部に一体に設けられているので、一般的な座金タイプの筒内圧センサと比較して、シリンダヘッド３ｃの振動の影響を抑制しながら、筒内圧ＰＣＹＬをより精度良く検出でき、したがって、筒内圧ＰＣＹＬを用いた熱発生率ｄＱｄθの演算精度をさらに向上させることができる。
【０１０２】  なお、本発明は、説明した実施形態に限定されることなく、種々の態様で実施することができる。例えば、実施形態は、燃焼パラメータとして熱発生率ｄＱｄθを算出し、熱発生率ｄＱｄθを用いて、エンジン制御としてＥＧＲ制御を実行する例である。本発明は、これに限らず、例えば、燃焼パラメータとして、図示平均有効圧力や燃焼トルク、あるいは筒内圧が最大になる最大筒内圧角、所定の燃焼質量割合が得られるクランク角（例えばＭＦＢ５０）や実着火時期などを算出してもよい。また、それらの算出結果に応じ、エンジン制御として、燃料噴射量や点火時期などを制御してもよい。
【０１０３】  また、実施形態では、ＣＰＳ演算部４１、モデル演算部４２及びエンジン制御部４３が、ＥＣＵ２内の複数のプロセッサコアにそれぞれ別個に搭載されているが、それらの全部又は一部を統合し、単一のＥＣＵ２内に設けてもよい。
【０１０４】  さらに、実施形態のエンジン３は４気筒のガソリンエンジンであるが、エンジン３のタイプや気筒３ａの数は任意である。また、実施形態では、筒内圧センサ２１は、すべての気筒３ａに設けられているが、一部の気筒３ａに設けられていてもよい。また、実施形態のエンジン３は、車両用のものであるが、本発明は、これに限らず、他の用途のエンジン、例えばクランクシャフトを鉛直方向に配置した船外機のような船舶推進機用エンジンなどにも適用できる。その他、本発明の趣旨の範囲内で、細部の構成を適宜、変更することが可能である。
【符号の説明】
【０１０５】  ２  ＥＣＵ（電子制御ユニット）  ３  エンジン（内燃機関）３ａ  気筒  ４  燃料噴射弁２１  筒内圧センサ４１  ＣＰＳ演算部（プロセッサコア、燃焼演算部）４２  モデル演算部（プロセッサコア、プラントモデル、同定手段）４３  エンジン制御部（プロセッサコア、制御器）      ＰＣＹＬ  筒内圧（筒内圧センサの検出結果）      ｄＱｄθ  熱発生率（燃焼パラメータ）        Ｐｉｎ  インマニ圧（機関パラメータ）      Ｔｅｇｒ  ＥＧＲ温度（機関パラメータ）      Ｐｅｇｒ  ＥＧＲ圧（機関パラメータ）ＰＭ１〜ＰＭ４  モデル基準点（モデルパラメータ）ＰＣ１〜ＰＣ４  補正基準点（筒内圧センサの検出結果）            θ  クランク角
【図１】
2017094349.tif 000035
【図２】
2017094349.tif 000036
【図３】
2017094349.tif 000037
【図４】
2017094349.tif 000038
【図５】
2017094349.tif 000039
【図６】
2017094349.tif 000040
【図７】
2017094349.tif 000041
【図８】
2017094349.tif 000042
【図９】
2017094349.tif 000043
【図１０】
2017094349.tif 000044
【図１１】
2017094349.tif 000045
【図１２】
2017094349.tif 000046
【図１３】
2017094349.tif 000047
【図１４】
2017094349.tif 000048
【図１５】
2017094349.tif 000049
【図１６】
2017094349.tif 000050
【図１７】
2017094349.tif 000051
【手続補正書】
【提出日】平成30年6月5日(2018.6.5)
【手続補正１】
【補正対象書類名】明細書
【補正対象項目名】００１３
【補正方法】変更
【補正の内容】
【００１３】
この構成によれば、１次関数モデル式の複数のモデルパラメータが、筒内圧センサの検出結果に基づいてリアルタイムに同定されるので、燃焼状態のばらつきや経年劣化などによる１次関数モデル式のモデル誤差を随時、適切に補償でき、熱発生率の演算精度を良好に維持することができる。
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】００２７
【補正方法】変更
【補正の内容】
【００２７】
また、吸気通路６のスロットル弁１０ａよりも上流側には、大気圧センサ２３及び外気温センサ２４が設けられている。大気圧センサ２３は大気圧ＰＡを検出し、外気温センサ２４は吸気通路６に導入される外気（新気）の温度ＴＡを検出し、それらの検出信号をＥＣＵ２に出力する。
【手続補正３】
【補正対象書類名】明細書
【補正対象項目名】００３２
【補正方法】変更
【補正の内容】
【００３２】
第２のプロセッサコア４２（以下「モデル演算部４２」という）は、後述するプラントモデルに基づき、エンジン３の状態を表すエンジンパラメータを演算するモデル演算処理を実行する。このエンジンパラメータには、吸気通路６、排気通路７及びＥＧＲ通路１２を流れる吸気、排ガス及びＥＧＲガスのそれぞれの質量流量、温度や圧力などが含まれる。
【手続補正４】
【補正対象書類名】明細書
【補正対象項目名】００３５
【補正方法】変更
【補正の内容】
【００３５】
上述したモデル演算処理の基礎となるプラントモデルは、エア系モデルと燃焼モデルに分類される。図４に示すように、エア系モデルは、吸気、排ガスやＥＧＲガスが流れるエンジン３の流路（吸気通路６、排気通路７及びＥＧＲ通路１２など）の構成を、スロットル弁１０ａやＥＧＲ弁１３ａなどが存在する「オリフィス」の部分とそれ以外の「レシーバ」の部分との組み合わせとして、モデル化したものである。また、レシーバに連続の式（質量保存則、エネルギ保存則）や気体の状態方程式などを適用し、オリフィスにオリフィスの式を適用することによって、エンジン３の流路の各部位における流体の質量流量、温度や圧力などが算出される。
【手続補正５】
【補正対象書類名】明細書
【補正対象項目名】００４５
【補正方法】変更
【補正の内容】
【００４５】
本処理では、まずステップ１（「Ｓ１」と図示。以下同じ」）において、吸気通路６のスロットル弁１０ａよりも下流側の吸気圧であるインマニ圧Ｐｉｎを算出する。このインマニ圧Ｐｉｎの算出は、図８に示すように、図１の吸気通路６のスロットル弁１０ａよりも下流側から吸気コレクタ部６ａまでの部分とＥＧＲ通路１２の接続部を、インマニモデル（レシーバ）として設定するとともに、インマニモデルにおいて成立する、以下に述べるようなパラメータ間の関係に基づいて行われる。
【手続補正６】
【補正対象書類名】明細書
【補正対象項目名】００５６
【補正方法】変更
【補正の内容】
【００５６】
図７に戻り、上記ステップ２に続くステップ３では、排気マニホルド７ｂ内の温度であるエキマニ温度Ｔｅｘを算出する。このエキマニ温度Ｔｅｘの算出は、図１０に示すように、排気通路７の排気マニホルド７ｂから排気コレクタ部７ａまでの部分とＥＧＲ通路１２の分岐部を、エキマニモデル（レシーバ）として設定するとともに、エキマニモデルにおいて成立する、以下に述べるようなパラメータ間の関係に基づいて行われる。
【手続補正７】
【補正対象書類名】明細書
【補正対象項目名】００８８
【補正方法】変更
【補正の内容】
【００８８】
また、第４モデル基準点ＰＭ４については、ステップ５７において、第４補正基準点ＰＣ４のクランク角θｍｉｎ２ａと第４モデル基準点ＰＭ４のクランク角θｍｉｎ２との差（＝θｍｉｎ２ａ−θｍｉｎ２）に、補正係数Ｋθを乗算することによって、クランク角補正項ΔθＣ４を算出し、ステップ５８において、第４補正基準点ＰＣ４の熱発生率ｄＱｄθｍｉｎ２ａと第４モデル基準点ＰＭ４の熱発生率ｄＱｄθｍｉｎ２との差（＝ｄＱｄθｍｉｎ２ａ−ｄＱｄθｍｉｎ２）に、補正係数ＫｄＱを乗算することによって、熱発生率補正項ΔｄＱＣ４を算出し、本処理を終了する。
【手続補正８】
【補正対象書類名】図面
【補正対象項目名】図１６
【補正方法】変更
【補正の内容】
【図１６】
2017094349.tif 000052
【国際調査報告】
2017094349.tif 000053
2017094349.tif 000054
2017094349.tif 000055
2017094349.tif 000056