(57)【要約】【課題】演算負荷の増加を抑制しつつ、点火時期が操作されても、燃焼室の出口ガス温度を精度よく推定できる内燃機関の制御装置及び制御方法を提供する。【解決手段】燃焼室ガス量Ｑａｌｌ、燃焼熱量ｈｌ＿ｂ、冷却損失率ηｃ、及び図示熱効率ηｉに基づいて、燃焼室の出口における排気ガスの温度である出口ガス温度Ｔｏｕｔを算出し、点火時期ＳＡの操作による図示熱効率の変化量Δηｉを算出し、図示熱効率の変化量Δηｉに基づいて冷却損失率の基本値ηｃ０を補正して、最終的な冷却損失率ηｃを算出する内燃機関の制御装置。【選択図】図５
2019132152.tif 000002
【特許請求の範囲】
【請求項１】
燃焼室に流入したガス量である燃焼室ガス量を算出する燃焼室ガス量算出部と、
燃料の燃焼により前記燃焼室で発生した燃焼熱量を算出する燃焼熱量算出部と、
前記燃焼熱量のうち、前記燃焼室の壁面に放熱される熱量の割合である冷却損失率を算出する冷却損失率算出部と、
前記燃焼熱量のうち、前記燃焼室内の筒内圧による仕事として取り出せる熱量の割合である図示熱効率を算出する図示熱効率算出部と、
前記燃焼室ガス量、前記燃焼熱量、前記冷却損失率、及び前記図示熱効率に基づいて、前記燃焼室の出口における排気ガスの温度である出口ガス温度を算出する熱勘定温度算出部と、
点火時期を操作する点火制御部と、を備え、
前記図示熱効率算出部は、前記点火時期の操作による前記図示熱効率の変化量を算出し、
前記冷却損失率算出部は、前記冷却損失率の基本値を算出し、前記図示熱効率の変化量に基づいて前記冷却損失率の基本値を補正して、最終的な前記冷却損失率を算出する内燃機関の制御装置。
【請求項２】
前記冷却損失率算出部は、前記図示熱効率の変化量に基づいて、前記点火時期の操作による前記冷却損失率の変化量を算出し、前記冷却損失率の基本値を前記冷却損失率の変化量により変化させて、最終的な前記冷却損失率を算出する請求項１に記載の内燃機関の制御装置。
【請求項３】
前記図示熱効率算出部は、予め設定された基本点火時期からの前記点火時期の遅角による前記図示熱効率の変化量を算出し、
前記冷却損失率算出部は、前記冷却損失率の基本値として、前記点火時期が前記基本点火時期に設定された場合の前記冷却損失率を算出し、前記図示熱効率の変化量に基づいて、前記基本点火時期からの前記点火時期の遅角による前記冷却損失率の変化量を算出する請求項２に記載の内燃機関の制御装置。
【請求項４】
前記基本点火時期は、トルクが最大になる前記点火時期に設定されている請求項３に記載の内燃機関の制御装置。
【請求項５】
前記冷却損失率算出部は、前記燃焼室から排気管に排出される排気ガス流量に基づいて、前記冷却損失率の基本値を算出する請求項１から４のいずれか一項に記載の内燃機関の制御装置。
【請求項６】
空燃比が理論空燃比よりリッチである場合に、エンリッチ量に応じて、前記出口ガス温度の低下補正を行う空燃比温度補正部を更に備えた請求項１から５のいずれか一項に記載の内燃機関の制御装置。
【請求項７】
前記燃焼室の出口から推定位置までの排気管の放熱による前記排気ガスの温度低下量を算出する放熱量算出部と、
前記出口ガス温度から前記温度低下量を減算して前記推定位置の排気ガス温度を推定する排気温度推定部と、を更に備えた請求項１から６のいずれか一項に記載の内燃機関の制御装置。
【請求項８】
前記放熱量算出部は、前記排気管内の排気ガス流量に基づいて、前記排気管内の排気ガスを加熱流体とし、前記排気管外の空気を受熱流体とした単流熱交換器としての前記排気管の温度効率を算出し、
前記温度効率に基づいて前記温度低下量を算出する請求項７に記載の内燃機関の制御装置。
【請求項９】
前記放熱量算出部は、前記排気管の温度効率をηとし、前記排気管内の排気ガス流量をＱｅｘとし、演算定数をＫηとして、
η＝１−ｅｘｐ（−Ｋη／Ｑｅｘ）
の算出式により、前記温度効率を算出し、
前記温度効率に基づいて前記温度低下量を算出する請求項７又は８に記載の内燃機関の制御装置。
【請求項１０】
前記放熱量算出部は、前記排気管内の排気ガス流量と前記排気管の温度効率との関係が予め設定された温度効率特性データを参照し、現在の前記排気ガス流量に対応する前記温度効率を算出し、
前記温度効率に基づいて前記温度低下量を算出する請求項７又は８に記載の内燃機関の制御装置。
【請求項１１】
前記放熱量算出部は、前記出口ガス温度から外気温を減算した値に、前記温度効率を乗算した値を、前記温度低下量として算出する請求項８から１０のいずれか一項に記載の内燃機関の制御装置。
【請求項１２】
前記排気ガス温度の推定値を用いて前記排気ガス温度の制御を行う排気温度制御、前記排気ガス温度の推定値を用いて排気ガスが流れるバルブの流量特性を算出するバルブ流量特性演算、前記排気ガス温度の推定値を用いて前記排気管に設けられた排気ガス温度センサの異常診断を行う排気温度センサ異常診断、及び前記排気ガス温度の推定値を用いて前記排気管に設けられた過給機のタービン出力を算出するタービン出力演算のいずれか１つ以上を実行する排気温度利用制御部を更に備えた請求項７から１１のいずれか一項に記載の内燃機関の制御装置。
【請求項１３】
燃焼室に流入したガス量である燃焼室ガス量を算出する流入ガス量算出ステップと、
燃料の燃焼により前記燃焼室で発生した燃焼熱量を算出する燃焼熱量算出ステップと、
前記燃焼熱量のうち、前記燃焼室の壁面に放熱される熱量の割合である冷却損失率を算出する冷却損失率算出ステップと、
前記燃焼熱量のうち、前記燃焼室内の筒内圧による仕事として取り出せる熱量の割合である図示熱効率を算出する図示熱効率算出ステップと、
前記燃焼室ガス量、前記燃焼熱量、前記冷却損失率、及び前記図示熱効率に基づいて、前記燃焼室の出口における排気ガスの温度である出口ガス温度を算出する出口ガス温度算出ステップと、
点火時期を操作する点火制御ステップと、を備え、
前記図示熱効率算出ステップでは、前記点火時期の操作による前記図示熱効率の変化量を算出し、
前記冷却損失率算出ステップでは、前記冷却損失率の基本値を算出し、前記図示熱効率の変化量に基づいて前記冷却損失率の基本値を補正して、最終的な前記冷却損失率を算出する内燃機関の制御方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本開示は、内燃機関の排気ガス温度を推定する内燃機関の制御装置及び制御方法に関する。
【背景技術】
【０００２】
内燃機関を好適に運転するためには、内燃機関の運転状態に応じて変化する各部の圧力や温度等の状態量を正確に把握することが重要である。近年では、排気系の状態量も用いて内燃機関の制御が行われるようになっている。排気系の状態量を用いる制御としては、例えば、高回転高負荷運転時に排気ガス温度が上昇し、排気系に設けられた触媒や空燃比センサ等を損傷するおそれがある場合は、空燃比をリッチ化して排気ガス温度を低減させるエンリッチ制御が行われている。この制御を行うために用いる排気ガス温度は、排気ガス温度センサを用いて検出してもよいが、センサコストを低減するため、推定する方法もある。例えば、回転速度及び充填効率と排気ガス温度との関係が、予め実験データ等に基づいて設定されたマップデータを参照して、排気ガス温度を推定する方法が用いられる。
【０００３】
排気ガス温度を用いる他の例としては、例えば特許文献１において、体積効率相当値を用いて燃焼室の吸入空気量や内部ＥＧＲ率を算出する方法が開示されているが、この体積効率相当値の算出に排気ガス温度が用いられている。この文献では、前述同様のマップデータを用いた排気ガス温度の推定方法が用いられている。これら以外にも、ＥＧＲバルブの開度に基づいて外部ＥＧＲの流量を算出するために排気ガス温度が用いられる場合もある。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特許第５４０９８３２号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
前述のように、排気ガス温度の推定方法としてマップデータを用いた推定方法があるが、排気ガス温度は、点火時期やＥＧＲ量等の影響を受けやすく、例えば回転速度と充填効率をマップ軸としたマップ設定だけでは、点火時期やＥＧＲ量等の運転条件が変化すると精度良く推定することが難しいという問題がある。点火時期やＥＧＲ量が変化した場合の排気ガス温度をすべてマップとして記憶するとしても、マップの数が膨大になるし、適合に必要な工数も膨大となるため、マップによる排気ガス温度の推定では推定精度や適合工数の面で問題がある。
【０００６】
排気ガス温度の推定精度が低いと、内部ＥＧＲ量や外部ＥＧＲ量の推定精度が低くなるため、ＥＧＲ量に基づく点火時期制御が適正に行われず、燃費効果を十分に得られないという問題がある。排気ガス温度の推定誤差により、過給圧を制御するウェイストゲート開度の制御が適正に行われない場合にも、ドライバの要求する加速感が得られないためにドライバビリティが悪化するという問題がある。これらの問題は、排気ガス温度の推定精度を向上することができれば解決され得る。
【０００７】
排気管内の排気ガス温度の推定精度を向上するためには、燃焼室の出口における排気ガス温度である出口ガス温度の推定精度を向上する必要がある。そこで、本願の発明者は、燃料の燃焼により発生する燃焼熱量の分配を考慮する熱勘定により、燃焼室の出口ガス温度の推定精度を向上する処理を検討している。具体的には、本願の発明者は、燃焼熱量のうち、燃焼室内の筒内圧による仕事として取り出せる熱量の割合を図示熱効率とし、膨張行程中に燃焼室の壁面へ放熱される熱量の割合を冷却損失率とし、残りが排気損失率として排気ガスの温度上昇に使われる熱量になる割合とした熱勘定を検討している。しかし、ノック制御等の点火時期の操作により、図示熱効率だけでなく、冷却損失率も変化するため、点火時期の操作による冷却損失率の変化を考慮しないと、出口ガス温度の推定精度が悪化する問題があった。具体的には、点火時期が遅角されると、図示熱効率が低下し、図示熱効率の低下分の熱量が出口ガス温度の上昇に用いられるが、図示熱効率の低下分の熱量の一部が冷却損失率の上昇にも用いられる。しかし、簡単な演算で、点火時期の操作から、直接的に、冷却損失率の変化量を精度よく算出することが容易でなかった。
【０００８】
そこで、演算負荷の増加、定数適合工数の増加を抑制しつつ、点火時期が操作されても、燃焼室の出口ガス温度を精度よく推定できる内燃機関の制御装置及び制御方法が望まれる。
【課題を解決するための手段】
【０００９】
本開示に係る内燃機関の制御装置は、燃焼室に流入したガス量である燃焼室ガス量を算出する燃焼室ガス量算出部と、
燃料の燃焼により前記燃焼室で発生した燃焼熱量を算出する燃焼熱量算出部と、
前記燃焼熱量のうち、前記燃焼室の壁面に放熱される熱量の割合である冷却損失率を算出する冷却損失率算出部と、
前記燃焼熱量のうち、前記燃焼室内の筒内圧による仕事として取り出せる熱量の割合である図示熱効率を算出する図示熱効率算出部と、
前記燃焼室ガス量、前記燃焼熱量、前記冷却損失率、及び前記図示熱効率に基づいて、前記燃焼室の出口における排気ガスの温度である出口ガス温度を算出する熱勘定温度算出部と、
点火時期を操作する点火制御部と、を備え、
前記図示熱効率算出部は、前記点火時期の操作による前記図示熱効率の変化量を算出し、
前記冷却損失率算出部は、前記冷却損失率の基本値を算出し、前記図示熱効率の変化量に基づいて前記冷却損失率の基本値を補正して、最終的な前記冷却損失率を算出するものである。
【００１０】
本開示に係る内燃機関の制御方法は、燃焼室に流入したガス量である燃焼室ガス量を算出する流入ガス量算出ステップと、
燃料の燃焼により前記燃焼室で発生した燃焼熱量を算出する燃焼熱量算出ステップと、
前記燃焼熱量のうち、前記燃焼室の壁面に放熱される熱量の割合である冷却損失率を算出する冷却損失率算出ステップと、
前記燃焼熱量のうち、前記燃焼室内の筒内圧による仕事として取り出せる熱量の割合である図示熱効率を算出する図示熱効率算出ステップと、
前記燃焼室ガス量、前記燃焼熱量、前記冷却損失率、及び前記図示熱効率に基づいて、前記燃焼室の出口における排気ガスの温度である出口ガス温度を算出する出口ガス温度算出ステップと、
点火時期を操作する点火制御ステップと、を備え、
前記図示熱効率算出ステップでは、前記点火時期の操作による前記図示熱効率の変化量を算出し、
前記冷却損失率算出ステップでは、前記冷却損失率の基本値を算出し、前記図示熱効率の変化量に基づいて前記冷却損失率の基本値を補正して、最終的な前記冷却損失率を算出するものである。
【発明の効果】
【００１１】
本開示に係る内燃機関の制御装置及び制御方法によれば、点火時期の操作による冷却損失率の変化と強い相関性のある図示熱効率の変化量に基づいて、冷却損失率を変化させるので、点火時期の操作により変化する冷却損失率を精度よく算出することができる。また、図示熱効率算出部による図示熱効率の算出を利用できるので、点火時期の操作から直接的に冷却損失率の変化を算出する場合よりも、演算負荷、定数適合工数を低減することができる。
【図面の簡単な説明】
【００１２】
【図１】実施の形態１に係る内燃機関及び制御装置の概略構成図である。
【図２】実施の形態１に係る制御装置のブロック図である。
【図３】実施の形態１に係る制御装置のハードウェア構成図である。
【図４】実施の形態１に係る熱勘定の考え方を示す図である。
【図５】実施の形態１に係る出口ガス温度算出部の詳細ブロック図である。
【図６】実施の形態１に係る変化量特性データの設定例を説明する図である。
【図７】実施の形態１に係る損失率特性データの設定例を説明する図である。
【図８】実施の形態１に係る排気管の単流熱交換器モデルを概略的に示す図である。
【図９】実施の形態１に係る制御装置の処理を示すフローチャート図である。
【発明を実施するための形態】
【００１３】
実施の形態１．
実施の形態１に係る内燃機関の制御装置５０（以下、単に制御装置５０と称す）について図面を参照して説明する。図１は、本実施の形態に係る内燃機関１の概略構成図であり、図２は、本実施の形態に係る制御装置５０のブロック図である。内燃機関１及び制御装置５０は、車両に搭載され、内燃機関１は、車両（車輪）の駆動力源となる。
【００１４】
１．内燃機関１の構成
図１に示すように、内燃機関１は、空気と燃料の混合気を燃焼する燃焼室２５を備えている。内燃機関１は、燃焼室２５に空気を供給する吸気管２３と、燃焼室２５で燃焼した排気ガスを排出する排気管１７とを備えている。燃焼室２５は、シリンダ（気筒）とピストンにより構成されている。以下では、燃焼室２５を気筒とも称す。内燃機関１は、ガソリンエンジンとされている。内燃機関１は、吸気管２３を開閉するスロットルバルブ６を備えている。スロットルバルブ６は、制御装置５０により制御される電気モータにより開閉駆動される電子制御式スロットルバルブとされている。スロットルバルブ６には、スロットルバルブ６の開度に応じた電気信号を出力するスロットル開度センサ７が設けられている。
【００１５】
スロットルバルブ６の上流側の吸気管２３には、吸気管２３に吸入される吸入空気流量に応じた電気信号を出力するエアフローセンサ３と、吸入空気の温度に応じた電気信号を出力する吸入空気温度センサ４と、が設けられている。吸入空気温度センサ４に検出された吸入空気の温度は、外気温Ｔａに等しいとみなすことができる。
【００１６】
内燃機関１は、排気管１７から吸気マニホールド１２に排気ガスを還流するＥＧＲ流路２１と、ＥＧＲ流路２１を開閉するＥＧＲバルブ２２と、を有している。吸気マニホールド１２は、スロットルバルブ６の下流側の吸気管２３の部分である。ＥＧＲバルブ２２は、制御装置５０により制御される電気モータにより開閉駆動される電子制御式ＥＧＲバルブとされている。ＥＧＲバルブ２２には、ＥＧＲバルブ２２の開度に応じた電気信号を出力するＥＧＲ開度センサ２７が設けられている。なお、ＥＧＲは、排気ガス再循環、すなわち、Exhaust Gas Recirculationの頭文字である。ＥＧＲバルブ２２を介して排気ガスが再循環するＥＧＲを、外部ＥＧＲといい、吸排気バルブのバルブオーバーラップにより燃焼室内に排気ガスが残留するＥＧＲを、内部ＥＧＲという。以下、外部ＥＧＲを単にＥＧＲと称す。
【００１７】
吸気マニホールド１２には、吸気マニホールド１２内の気体の圧力であるマニホールド圧に応じた電気信号を出力するマニホールド圧センサ８と、吸気マニホールド１２内の気体の温度であるマニホールド温度Ｔｉｎに応じた電気信号を出力するマニホールド温度センサ９と、が設けられている。
【００１８】
内燃機関１には、燃焼室２５に燃料を供給するインジェクタ１３が設けられている。インジェクタ１３は、燃焼室２５内に直接燃料を噴射するように設けられている。インジェクタ１３は、吸気マニホールド１２の下流側の部分に燃料を噴射するように設けられてもよい。内燃機関１には、大気圧に応じた電気信号を出力する大気圧センサ２が設けられている。
【００１９】
燃焼室２５の頂部には、空気と燃料の混合気に点火する点火プラグと、点火プラグに点火エネルギーを供給する点火コイル１６と、が設けられている。また、燃焼室２５の頂部には、吸気管２３から燃焼室２５内に吸入される吸入空気量を調節する吸気バルブ１４と、燃焼室２５から排気管１７に排出される排気ガス量を調節する排気バルブ１５と、が設けられている。吸気バルブ１４には、そのバルブ開閉タイミングを可変にする吸気可変バルブタイミング機構が設けられている。排気バルブ１５には、そのバルブ開閉タイミングを可変にする排気可変バルブタイミング機構が設けられている。可変バルブタイミング機構１４、１５は、電動アクチュエータを有している。内燃機関１のクランク軸には、その回転角に応じた電気信号を出力するクランク角センサ２０が設けられている。
【００２０】
排気管１７には、排気ガス中の空気と燃料との比率である空燃比ＡＦ（Ａｉｒ／Ｆｕｅｌ）に応じた電気信号を出力する空燃比センサ１８が設けられている。また、排気管１７には、排気ガスを浄化する触媒１９が設けられている。
【００２１】
シリンダブロックにはノックセンサ２８が固定されている。ノックセンサ２８は、内燃機関１の振動に応じた信号（振動波形信号）を出力する。ノックセンサ２８は、圧電素子等により構成される。
【００２２】
２．制御装置５０の構成
次に、制御装置５０について説明する。制御装置５０は、内燃機関１を制御対象とする制御装置である。図２のブロック図に示すように、制御装置５０は、運転状態検出部５１、出口ガス温度算出部５２、放熱量算出部５３、排気温度推定部５４、排気温度利用制御部５５、及び点火制御部５６等の制御部を備えている。制御装置５０の各制御部５１〜５６等は、制御装置５０が備えた処理回路により実現される。具体的には、制御装置５０は、図３に示すように、処理回路として、ＣＰＵ（Central Processing Unit）等の演算処理装置９０（コンピュータ）、演算処理装置９０とデータのやり取りをする記憶装置９１、演算処理装置９０に外部の信号を入力する入力回路９２、及び演算処理装置９０から外部に信号を出力する出力回路９３等を備えている。
【００２３】
演算処理装置９０として、ＡＳＩＣ（Application Specific Integrated Circuit）、ＩＣ（Integrated Circuit）、ＤＳＰ（Digital Signal Processor）、ＦＰＧＡ（Field Programmable Gate Array）、各種の論理回路、及び各種の信号処理回路等が備えられてもよい。また、演算処理装置９０として、同じ種類のもの又は異なる種類のものが複数備えられ、各処理が分担して実行されてもよい。記憶装置９１として、演算処理装置９０からデータを読み出し及び書き込みが可能に構成されたＲＡＭ（Random Access Memory）や、演算処理装置９０からデータを読み出し可能に構成されたＲＯＭ（Read Only Memory）等が備えられている。入力回路９２は、各種のセンサやスイッチが接続され、これらセンサやスイッチの出力信号を演算処理装置９０に入力するＡ／Ｄ変換器等を備えている。出力回路９３は、電気負荷が接続され、これら電気負荷に演算処理装置９０から制御信号を出力する駆動回路等を備えている。
【００２４】
そして、制御装置５０が備える各制御部５１〜５６等の各機能は、演算処理装置９０が、ＲＯＭ等の記憶装置９１に記憶されたソフトウェア（プログラム）を実行し、記憶装置９１、入力回路９２、及び出力回路９３等の制御装置５０の他のハードウェアと協働することにより実現される。なお、各制御部５１〜５６等が用いる特性データ、定数等の設定データは、ソフトウェア（プログラム）の一部として、ＲＯＭ等の記憶装置９１に記憶されている。
【００２５】
本実施の形態では、入力回路９２には、大気圧センサ２、エアフローセンサ３、吸入空気温度センサ４、スロットル開度センサ７、マニホールド圧センサ８、マニホールド温度センサ９、空燃比センサ１８、クランク角センサ２０、アクセルポジションセンサ２６、ＥＧＲ開度センサ２７、及びノックセンサ２８等が接続されている。出力回路９３には、スロットルバルブ６（電気モータ）、インジェクタ１３、吸気可変バルブタイミング機構１４、排気可変バルブタイミング機構１５、点火コイル１６、及びＥＧＲバルブ２２（電動アクチュエータ）等が接続されている。なお、制御装置５０には、図示していない各種のセンサ、スイッチ、及びアクチュエータ等が接続されている。
【００２６】
運転状態検出部５１は、内燃機関１の運転状態を検出する。運転状態検出部５１は、各種のセンサの出力信号等に基づいて各種の運転状態を検出する。具体的には、運転状態検出部５１は、大気圧センサ２の出力信号に基づいて大気圧を検出し、エアフローセンサ３の出力信号に基づいて吸入空気流量を検出し、吸入空気温度センサ４の出力信号に基づいて外気温Ｔａを検出し、スロットル開度センサ７の出力信号に基づいてスロットル開度を検出し、マニホールド圧センサ８の出力信号に基づいてマニホールド圧を検出し、マニホールド温度センサ９の出力信号等に基づいて吸気マニホールド１２内の気体の温度であるマニホールド温度Ｔｉｎを検出し、空燃比センサ１８の出力信号に基づいて、排気ガスの空燃比ＡＦを検出し、クランク角センサ２０の出力信号に基づいてクランク角度及び回転速度Ｎｅを検出し、アクセルポジションセンサ２６の出力信号に基づいてアクセル開度を検出し、ＥＧＲ開度センサ２７の出力信号に基づいてＥＧＲ開度を検出し、ノックセンサ２８の出力信号に基づいてノック強度ＫＮＫを検出する。
【００２７】
運転状態検出部５１は、吸入空気流量、回転速度Ｎｅ等に基づいて燃焼室２５内に流入した空気量である吸入空気量Ｑｃ［ｇ／ｓｔｒｏｋｅ］、充填効率Ｅｃ［％］を算出する。例えば、運転状態検出部５１は、吸入空気流量［ｇ／ｓ］に、回転速度Ｎｅに応じた行程周期を乗算した値に、吸気マニホールドの遅れを模擬したフィルタ処理を行った値を、吸入空気量Ｑｃ［ｇ／ｓｔｒｏｋｅ］として算出する。行程周期は、３気筒エンジンであれば、２４０ｄｅｇＣＡ間の周期となり、４気筒エンジンであれば、１８０ｄｅｇＣＡ間の周期となる。或いは、運転状態検出部５１は、マニホールド圧、回転速度Ｎｅ等に基づいて吸入空気量Ｑｃ［ｇ／ｓｔｒｏｋｅ］、充填効率Ｅｃ［％］を算出してもよい。
【００２８】
運転状態検出部５１は、ＥＧＲ開度等に基づいて、燃焼室２５内に流入した排気ガス再循環量であるＥＧＲ量Ｑｃｅ［ｇ／ｓｔｒｏｋｅ］を算出する。例えば、運転状態検出部５１は、ＥＧＲ開度及びマニホールド圧等に基づいて、ＥＧＲバルブ２２を通過するＥＧＲ流量［ｇ／ｓ］を算出し、ＥＧＲ流量に行程周期を乗算した値に、フィルタ処理を行った値を、ＥＧＲ量Ｑｃｅ［ｇ／ｓｔｒｏｋｅ］として算出する。運転状態検出部５１は、吸入空気量Ｑｃに対するＥＧＲ量Ｑｃｅの比率であるＥＧＲ率Ｒｅｇｒ［％］を算出する。
【００２９】
制御装置５０は、算出された回転速度Ｎｅ、吸入空気量Ｑｃ、充填効率Ｅｃ、ＥＧＲ率Ｒｅｇｒ等に基づいて、目標空燃比、燃料噴射量Ｑｆ、点火時期ＳＡ等を算出し、インジェクタ１３及び点火コイル１６等を駆動制御する。例えば、制御装置５０は、回転速度Ｎｅ及び充填効率Ｅｃに基づいて目標空燃比を算出し、吸入空気量Ｑｃを目標空燃比で除算して燃料噴射量の基本値を算出する。そして、制御装置５０は、空燃比センサ１８を用いた空燃比フィードバック制御を行う場合は、燃料噴射量の基本値を補正して、最終的な燃料噴射量Ｑｆを算出し、空燃比フィードバック制御を行わない場合は、燃料噴射量の基本値をそのまま最終的な燃料噴射量Ｑｆに設定する。
【００３０】
制御装置５０は、アクセル開度等に基づいて、目標スロットル開度を算出し、スロットル開度が目標スロットル開度に近づくように、スロットルバルブ６の電気モータを駆動制御する。制御装置５０は、回転速度Ｎｅ及び充填効率Ｅｃ等に基づいてＥＧＲバルブ２２の目標ＥＧＲ開度を算出し、ＥＧＲ開度が目標ＥＧＲ開度に近づくように、ＥＧＲバルブ２２の電動アクチュエータを駆動制御する。制御装置５０は、回転速度Ｎｅ及び充填効率Ｅｃ等に基づいて、吸気バルブ１４及び排気バルブ１５のそれぞれの目標開閉タイミング（位相）を算出し、吸気バルブ１４及び排気バルブ１５の開閉タイミングが、それぞれの目標開閉タイミングに近づくように、吸気及び排気可変バルブタイミング機構の電動アクチュエータを駆動制御する。
【００３１】
トルクベース制御を行う場合は、制御装置５０は、アクセル開度等に基づいて算出された内燃機関１の要求出力トルク、又は変速機制御装置等の外部の制御装置から要求された内燃機関１の要求出力トルクを実現するように、スロットル開度、点火時期ＳＡ、ＥＧＲ開度、吸気バルブ１４及び排気バルブ１５の開閉タイミングを制御する。具体的には、制御装置５０は、適合により予め設定された充填効率Ｅｃと図示熱効率ηｉとの関係に基づいて、要求出力トルクに対応する図示熱効率を実現するために必要な目標充填効率を算出し、目標充填効率で運転した場合に、燃費、排出ガスが最適となる目標ＥＧＲ率を算出し、目標充填効率及び目標ＥＧＲ率を達成するような目標スロットル開度、目標ＥＧＲ開度、吸気バルブ１４の目標開閉タイミング、及び排気バルブ１５の目標開閉タイミングを算出し、これらの目標値に基づいて、スロットルバルブ６の電気モータ、ＥＧＲバルブ２２の電動アクチュエータ、吸気及び排気可変バルブタイミング機構の電動アクチュエータを駆動制御する。
【００３２】
２−１．点火制御部５６
点火時期ＳＡの制御について、より詳細に説明する。点火制御部５６は、点火時期ＳＡを操作する。本実施の形態では、点火制御部５６は、回転速度Ｎｅ、充填効率Ｅｃ、及びＥＧＲ率Ｒｅｇｒ等の運転状態と、基本点火時期ＳＡ０との関係が予め設定された点火時期特性データを参照し、現在の運転状態に対応する基本点火時期ＳＡ０を算出する。基本点火時期ＳＡ０は、基本的に、トルクが最大になる点火時期（ＭＢＴ：Minimum advance for the Best Torque）に設定されている。
【００３３】
点火制御部５６は、ノックセンサ２８により検出したノック強度ＫＮＫに応じてノック遅角量を算出し、基本点火時間ＳＡ０からノック遅角量だけ遅角した点火時期を最終的な点火時期ＳＡとして算出する。また、点火制御部５６は、内燃機関１の暖機運転のため、触媒１９の昇温のため、及び内燃機関１のトルク低下のために、点火時期ＳＡを、基本点火時間ＳＡ０よりも遅角させる遅角制御を行う。点火制御部５６は、点火時期ＳＡで点火プラグに火花が生じるように、点火時期ＳＡ及びクランク角度に基づいて、点火コイル１６への通電制御を行う。
【００３４】
２−２．出口ガス温度算出部５２
＜出口ガス温度Ｔｏｕｔの算出方法＞
出口ガス温度Ｔｏｕｔの算出方法について説明する。燃焼室内に供給された燃料が燃焼することにより発生する燃焼熱量ｈｌ＿ｂがどのように分配されるかを示したものを熱勘定といい、熱勘定の考え方を図４に示す。燃焼熱量ｈｌ＿ｂのうち、燃焼室内の筒内圧による仕事として取り出せる熱量の割合を図示熱効率ηｉ［％］とし、膨張行程中に燃焼室の壁面へ放熱される熱量の割合を冷却損失率ηｃ［％］とし、残りが排気損失率ηｅｘ［％］として排気ガスの温度上昇に使われる熱量になる割合とする。なお、ポンピングロスや機械損失は、一旦仕事として取り出したものの、軸出力以外に使われる仕事と考えて、図示熱効率ηｉに含むものとする。ここで、図示熱効率ηｉ、冷却損失率ηｃが分かるものとすると、排気損失率ηｅｘは、次式で表せる。
【数１】
2019132152.tif 000003
【００３５】
この時、吸気マニホールド１２に導入された気体の温度であるマニホールド気温をＴｉｎ［Ｋ］とし、今回の燃焼により発生する燃焼熱量をｈｌ＿ｂ［Ｊ／ｓｔｒｏｋｅ］とし、今回の燃焼における燃焼室内ガス量をＱａｌｌ［ｇ／ｓｔｒｏｋｅ］とし、燃焼室内ガス（排気ガス）の比熱をＣｅｘ［Ｊ／（ｇ・Ｋ）］とすると、出口ガス温度の基本値Ｔｏｕｔ０［Ｋ］は、次式で算出できる。
【数２】
2019132152.tif 000004
【００３６】
ここで、（２）式で算出する出口ガス温度を基本値としたのは、空燃比が理論空燃比又はリーンである場合は、（２）式で算出した出口ガス温度をそのまま用いることができるが、空燃比がリッチである場合は、エンリッチにより、出口ガス温度が低下するため、その分の補正を行う必要があるためである。エンリッチによる温度低下の理由は、燃料の気化による蒸発熱や、未燃燃料が分子量の小さい炭化水素に分解される際に消費されるエネルギーによるものと考えられる。結局、エンリッチによる温度低下量をエンリッチ温度低下量ΔＴｒｉｃｈとすると、出口ガス温度Ｔｏｕｔは次式により算出することができる。以上により出口ガス温度Ｔｏｕｔが算出できるので、排気ガス温度Ｔｅｘも算出できることになる。
【数３】
2019132152.tif 000005
【００３７】
＜出口ガス温度算出部５２＞
以上で導出した算出方法に基づいて設計された制御装置５０について説明する。出口ガス温度算出部５２は、燃焼室の出口における排気ガスの温度である出口ガス温度Ｔｏｕｔを算出する。図５に、出口ガス温度算出部５２の詳細ブロック図を示すように、出口ガス温度算出部５２は、燃焼室ガス量算出部５２ａ、燃焼熱量算出部５２ｂ、冷却損失率算出部５２ｃ、図示熱効率算出部５２ｄ、熱勘定温度算出部５２ｅ、及び空燃比温度補正部５２ｆを備えている。
【００３８】
＜燃焼室ガス量算出部５２ａ＞
燃焼室ガス量算出部５２ａは、燃焼室に流入したガス量である燃焼室ガス量Ｑａｌｌを算出する。本実施の形態では、燃焼室ガス量算出部５２ａは、次式に示すように、燃焼室に流入した空気量である吸入空気量Ｑｃ［ｇ／ｓｔｒｏｋｅ］、燃焼室に流入した排気ガス再循環量であるＥＧＲ量Ｑｃｅ［ｇ／ｓｔｒｏｋｅ］、及び燃焼室に供給した燃料噴射量Ｑｆ［ｇ／ｓｔｒｏｋｅ］を合計した値を、燃焼室ガス量Ｑａｌｌ［ｇ／ｓｔｒｏｋｅ］として算出する。なお、吸入空気量Ｑｃ及びＥＧＲ量Ｑｃｅは、上述したように、運転状態検出部５１により算出される。
【数４】
2019132152.tif 000006
【００３９】
また、燃焼室ガス量算出部５２ａは、次式に示すように、燃焼室ガス量Ｑａｌｌを行程周期ΔＴｓｇｔで除算した値を、排気行程で排出される排気ガス流量Ｑｅｘ［ｇ／ｓ］として算出する。行程周期ΔＴｓｇｔは、３気筒エンジンであれば、２４０ｄｅｇＣＡ間の周期となり、４気筒エンジンであれば、１８０ｄｅｇＣＡ間の周期となる。ここで、（ｎ）は、現在の行程を表し、（ｎ−３）は、現在の行程から３行程前の行程を表している。Ｑｅｘ（ｎ）は、現在排気行程である気筒から排出されている排気ガス流量を表し、Ｑａｌｌ（ｎ−３）は、現在排気行程である気筒が、３行程前に吸気行程であった時に流入した燃焼室ガス量を表している。
【数５】
2019132152.tif 000007
【００４０】
＜燃焼熱量算出部５２ｂ＞
燃焼熱量算出部５２ｂは、燃料の燃焼により燃焼室で発生した燃焼熱量ｈｌ＿ｂを算出する。燃焼熱量算出部５２ｂは、燃料噴射量Ｑｆ［ｇ／ｓｔｒｏｋｅ］及び空燃比ＡＦに基づいて、燃料噴射量Ｑｆのうち、実際に燃焼した燃焼燃料量Ｑｆｂ［ｇ／ｓｔｒｏｋｅ］を算出する。空燃比ＡＦが理論空燃比又はリーンである場合は、燃料噴射量Ｑｆの全てが燃焼すると仮定できるが、空燃比ＡＦがリッチである場合は、理論空燃比分の燃料量は燃焼するが、理論空燃比を超えたリッチ分の燃料量は燃焼しない。そこで、燃焼熱量算出部５２ｂは、次式に示すように、空燃比ＡＦが理論空燃比又はリーンである場合は、燃料噴射量Ｑｆをそのまま燃焼燃料量Ｑｆｂに設定し、空燃比ＡＦがリッチである場合は、燃料噴射量Ｑｆに、空燃比ＡＦを理論空燃比ＡＦ０で除算した値を、燃焼燃料量Ｑｆｂとして算出する。空燃比ＡＦには、燃料噴射量Ｑｆの演算に用いられる目標空燃比が用いられてもよいし、空燃比センサ１８を用いて検出された排気ガスの空燃比が用いられてもよい。
【数６】
2019132152.tif 000008
【００４１】
そして、燃焼熱量算出部５２ｂは、燃焼燃料量Ｑｆｂ［ｇ／ｓｔｒｏｋｅ］に単位発熱量を乗算した値を、燃焼熱量ｈｌ＿ｂ［Ｊ／ｓｔｒｏｋｅ］として算出する。ガソリンの単位発熱量は、４４０００［Ｊ／ｇ］程度に設定される。なお、燃料の気化熱量を、単位発熱量から減算してもよいが、ガソリンの気化熱量は２７２［Ｊ／ｇ］程度であるので無視してもよい。
【００４２】
＜図示熱効率算出部５２ｄ＞
図示熱効率算出部５２ｄは、燃焼熱量ｈｌ＿ｂのうち、燃焼室内の筒内圧による仕事として取り出せる熱量の割合である図示熱効率ηｉを算出する。図示熱効率算出部５２ｄは、回転速度Ｎｅ、充填効率Ｅｃ、ＥＧＲ率Ｒｅｇｒ、及び点火時期ＳＡ等の運転状態と、図示熱効率ηｉとの関係が予め設定された図示熱効率特性データを参照し、現在の運転状態に対応する現在の図示熱効率ηｉを算出する。ここで、現在の充填効率Ｅｃ、ＥＧＲ率Ｒｅｇｒには、現在排気行程である気筒が、３行程前に吸気行程であった時の値が用いられる。
【００４３】
図示熱効率特性データは、トルクベース制御において出力トルクの演算に用いるものが用いられ、実験データに基づいて予め設定されている。図示熱効率ηｉの計測は、全運転ポイントを計測すると膨大となるため、ＭＢＣ（Model Based Calibration）、ＤｏＥ（Design of Experiments）といった手法を用いて計測ポイントを抑制し、計測したポイントでのデータを近似式で繋ぐようにして全運転ポイントに相当する近似データを作成する。図示熱効率特性データには、更に簡易化した近似式が用いられる。なお、トルクベース制御が行われない場合は、図示熱効率特性データは、回転速度Ｎｅ及び充填効率Ｅｃと、図示熱効率ηｉとの関係が設定されたものとされてもよい。なお、図示熱効率ηｉには、燃焼により発生する熱量を分母にし、図示平均有効圧に相当する熱量を分子にして算出したものが用いられる。
【００４４】
図示熱効率算出部５２ｄは、点火時期ＳＡの操作による図示熱効率の変化量Δηｉを算出する。本実施の形態では、図示熱効率算出部５２ｄは、予め設定された基本点火時期ＳＡ０からの点火時期の遅角による図示熱効率の変化量Δηｉを算出するように構成されている。本実施の形態では、基本点火時期ＳＡ０は、トルクが最大になる点火時期（ＭＢＴ）に設定されている。図示熱効率算出部５２ｄは、点火制御部５６が算出した基本点火時期ＳＡ０を用いてもよいし、点火制御部５６と同様に、運転状態に基づいて、基本点火時期ＳＡ０を算出してもよい。
【００４５】
図示熱効率算出部５２ｄは、基本点火時期ＳＡ０に対応する図示熱効率ηｉ０（以下、基本図示熱効率ηｉ０と称す）を算出する。図示熱効率算出部５２ｄは、上述した図示熱効率特性データを参照し、点火時期以外の現在の運転状態（回転速度Ｎｅ、充填効率Ｅｃ、及びＥＧＲ率Ｒｅｇｒ等）、及び基本点火時期ＳＡ０に対応する図示熱効率を基本図示熱効率ηｉ０として算出する。
【００４６】
図示熱効率算出部５２ｄは、次式に示すように、現在の点火時期ＳＡに対応する図示熱効率ηｉから、基本点火時期ＳＡ０に対応する基本図示熱効率ηｉ０を減算して、図示熱効率の変化量Δηｉを算出する。点火時期ＳＡを基本点火時期ＳＡ０から遅角すると、図示熱効率ηｉが減少するので、図示熱効率の変化量Δηｉは負の値となる。
【数７】
2019132152.tif 000009
【００４７】
＜冷却損失率算出部５２ｃ＞
冷却損失率算出部５２ｃは、燃焼により発生した燃焼熱量ｈｌ＿ｂのうち、燃焼室の壁面に放熱される熱量の割合である冷却損失率ηｃを算出する。
【００４８】
ところで、ノック制御等の点火時期ＳＡの操作により、図示熱効率だけでなく、冷却損失率も変化するため、点火時期ＳＡの操作による冷却損失率の変化を考慮しないと、出口ガス温度の推定精度が悪化する問題があった。具体的には、点火時期ＳＡが遅角されると、図示熱効率が低下し、図示熱効率の低下分の熱量が出口ガス温度の上昇に用いられるが、図示熱効率の低下分の熱量の一部が冷却損失率の上昇にも用いられる。
【００４９】
しかし、簡単な演算で、点火時期ＳＡの操作から、直接的に、冷却損失率の変化量を精度よく算出することが容易でなかった。発明者が実施した実験結果によれば、冷却損失率の変化は、図示熱効率の変化と強い相関性があるため、冷却損失率の変化を精度よく算出するためには、図示熱効率の変化を考慮する必要があることが分かった。これは、点火時期ＳＡの操作による図示熱効率の変化量に応じて、燃焼室内のガス温度が上昇し、燃焼室内のガス温度の上昇量に応じて、燃焼室の壁面へ放熱される冷却損失率も上昇するためだと考えられる。すなわち、点火時期ＳＡの操作による図示熱効率の変化量に応じて、冷却損失率の変量が定まることとなる。なお、図示熱効率は、回転速度Ｎｅ、充填効率Ｅｃ、ＥＧＲ率Ｒｅｇｒ、及び点火時期ＳＡ等の複数のパラメータによって変化するため、直接的に冷却損失率の変化を精度よく算出するためには、複雑な演算及び定数適合が必要であることがわかる。
【００５０】
そこで、本実施の形態では、冷却損失率算出部５２ｃは、冷却損失率の基本値ηｃ０を算出し、点火時期ＳＡの操作による図示熱効率の変化量Δηｉに基づいて冷却損失率の基本値ηｃ０を補正して、最終的な冷却損失率ηｃを算出するように構成されている。
【００５１】
この構成によれば、点火時期の変化による冷却損失率の変化と強い相関性のある、図示熱効率の変化量Δηｉに基づいて冷却損失率を変化させるので、点火時期の操作により変化する冷却損失率を精度よく算出することができる。また、図示熱効率算出部５２ｄによる図示熱効率の算出を利用できるので、直接的に冷却損失率の変化を算出する場合よりも、演算負荷、定数適合工数を低減することができる。
【００５２】
本実施の形態では、冷却損失率算出部５２ｃは、冷却損失率の基本値ηｃ０として、点火時期ＳＡが基本点火時期ＳＡ０に設定された場合の冷却損失率を算出する。そして、冷却損失率算出部５２ｃは、基本点火時期ＳＡ０からの点火時期の遅角による図示熱効率の変化量Δηｉに基づいて、基本点火時期ＳＡ０からの点火時期の遅角による冷却損失率の変化量Δηｃを算出する。そして、冷却損失率算出部５２ｃは、次式に示すように、冷却損失率の基本値ηｃ０を冷却損失率の変化量Δηｃにより変化させて、最終的な冷却損失率ηｃを算出する。
【数８】
2019132152.tif 000010
【００５３】
この構成によれば、基本点火時期ＳＡ０を共通の基準にして、冷却損失率の基本値ηｃ０、図示熱効率の変化量Δηｉ、及び冷却損失率の変化量Δηｃを算出するので、演算処理及びデータ設定を簡単化することができる。
【００５４】
発明者が実施した実験結果によって、基本点火時期ＳＡ０からの点火時期の遅角による図示熱効率の変化量Δηｉと、基本点火時期ＳＡ０からの点火時期の遅角による冷却損失率の変化量Δηｃとの間には、強い相関性があることが見出された。そこで、冷却損失率算出部５２ｃは、図示熱効率の変化量Δηｉと冷却損失率の変化量Δηｃとの関係が予め設定された変化量特性データを参照し、図示熱効率算出部５２ｄによって算出された図示熱効率の変化量Δηｉに対応する冷却損失率の変化量Δηｃを算出する。図６に示すように、変化量特性データには、図示熱効率の変化量Δηｉが０から減少するに従って、冷却損失率の変化量Δηｃが０から増加する特性が、実験データに基づいて予め設定されている。なお、各特性データには、データマップ、データテーブル、多項式、数式等が用いられ、その設定データは、記憶装置９１に記憶されている。上述したように、基本点火時期ＳＡ０は、トルクが最大になる点火時期（ＭＢＴ）に設定されている。
【００５５】
発明者が実施した実験結果から、排気ガス流量Ｑｅｘと、逆算により算出した冷却損失率の基本値ηｃ０との間には、運転条件に依らず、強い相関関係があることが見いだされた。そこで、冷却損失率算出部５２ｃは、排気管内の排気ガス流量Ｑｅｘに基づいて冷却損失率の基本値ηｃ０を算出する。具体的には、冷却損失率算出部５２ｃは、排気ガス流量Ｑｅｘと冷却損失率の基本値ηｃ０との関係が予め設定された損失率特性データを参照し、現在の排気ガス流量Ｑｅｘに対応する冷却損失率の基本値ηｃ０を算出する。図７に示すように、損失率特性データには、排気ガス流量Ｑｅｘが増加するに従って、冷却損失率の基本値ηｃ０が低下する特性が、実験データに基づいて予め設定されている。
【００５６】
＜熱勘定温度算出部５２ｅ＞
熱勘定温度算出部５２ｅは、燃焼室ガス量Ｑａｌｌ［ｇ／ｓｔｒｏｋｅ］、燃焼熱量ｈｌ＿ｂ［Ｊ／ｓｔｒｏｋｅ］、冷却損失率ηｃ、及び図示熱効率ηｉに基づいて、燃焼室の出口における排気ガスの温度である出口ガス温度Ｔｏｕｔを算出する。具体的には、熱勘定の考え方に基づく（１）式、（２）式を用いて上述したように、熱勘定温度算出部５２ｅは、次式に示すように、１００［％］から図示熱効率ηｉ［％］及び冷却損失率ηｃ［％］を減算した値を、燃焼熱量のうち排気ガスの温度上昇に使われる熱量の割合である排気損失率ηｅｘ［％］として算出する。そして、熱勘定温度算出部５２ｅは、燃焼熱量ｈｌ＿ｂ［Ｊ／ｓｔｒｏｋｅ］に、排気損失率ηｅｘ［％］を乗算して、排気ガスの温度上昇に用いられる温度上昇熱量［Ｊ／ｓｔｒｏｋｅ］を算出し、温度上昇熱量［Ｊ／ｓｔｒｏｋｅ］を、燃焼室ガス量Ｑａｌｌ［ｇ／ｓｔｒｏｋｅ］と排気ガスの比熱Ｃｅｘ［Ｊ／（ｇ・Ｋ）］を乗算した水当量［Ｊ／（ｓｔｒｏｋｅ・ｇ）］で除算して温度上昇量［Ｋ］を算出し、温度上昇量［Ｋ］をマニホールド温度Ｔｉｎ［Ｋ］に加算して、出口ガス温度の基本値Ｔｏｕｔ０［Ｋ］を算出する。
【数９】
2019132152.tif 000011
【００５７】
排気ガスの比熱Ｃｅｘには、１．１［Ｊ／（ｇ・Ｋ）］程度の値が設定される。排気ガスの比熱Ｃｅｘは、厳密には空燃比に応じて変化するため、空気及び燃料の比熱を用いて補正するなど、空燃比ＡＦに応じて変化させてもよいが、簡単のために、同じ固定値が用いられてもよい。マニホールド温度温Ｔｉｎとして、マニホールド温度センサ９により検出された吸気マニホールド１２内のガス温度がそのまま用いられてもよいが、燃焼室内に流入したガス温度に近づくように、吸気マニホールド１２内のガス温度に所定値が加算された温度がマニホールド温Ｔｉｎとして用いられてもよい。或いは、吸排気バルブの開閉タイミングに基づいて算出された内部ＥＧＲ率に応じて吸気マニホールド１２内のガス温度が補正された温度が、マニホールド温Ｔｉｎとして用いられてもよい。マニホールド温度温Ｔｉｎとして、外気温Ｔａに基づいて推定されたガス温度が用いられてもよい。
【００５８】
＜空燃比温度補正部５２ｆ＞
空燃比温度補正部５２ｆは、空燃比ＡＦが理論空燃比よりリッチである場合に、エンリッチ量ΔＡＦｒに応じて、出口ガス温度Ｔｏｕｔの低下補正を行う。本願の発明者が実施した実験結果から、運転条件に依らず、エンリッチ量ΔＡＦｒが１増加すると、出口ガス温度Ｔｏｕｔが３５〜４０℃低下する特性となった。空燃比温度補正部５２ｆは、エンリッチ量ΔＡＦｒ（＝ＡＦ０−ＡＦ）とエンリッチ温度低下量ΔＴｒｉｃｈとの関係が予め設定された温度低下特性データを参照し、現在のエンリッチ量ΔＡＦｒに対応するエンリッチ温度低下量ΔＴｒｉｃｈを算出する。温度低下特性データには、エンリッチ量ΔＡＦｒが増加するに従って、正の値のエンリッチ温度低下量ΔＴｒｉｃｈが増加する特性が、実験データに基づいて予め設定されている。
【００５９】
空燃比温度補正部５２ｆは、次式に示すように、空燃比ＡＦが理論空燃比又はリーンである場合は、出口ガス温度の基本値Ｔｏｕｔ０をそのまま出口ガス温度Ｔｏｕｔに設定し、空燃比ＡＦがリッチである場合は、出口ガス温度の基本値Ｔｏｕｔ０から、正の値のエンリッチ温度低下量ΔＴｒｉｃｈを減算した値を、出口ガス温度Ｔｏｕｔとして算出する。
【数１０】
2019132152.tif 000012
【００６０】
２−３．排気ガス温度の推定演算
制御装置５０は、排気ガス温度Ｔｅｘの推定を行うように構成されている。推定された排気ガス温度Ｔｅｘは、後述する排気温度利用制御部５５において、ＥＧＲ率の算出等に用いられる。
【００６１】
＜単流熱交換器モデルを用いた排気ガス温度の算出方法＞
まず、単流熱交換器モデルを用いた、排気ガス温度の算出方法について説明する。なお、単流熱交換器は、詳しくは、「大学講義  伝熱工学」（丸善株式会社、Ｐ２２４〜２２６、１９８３年）に記載されている。
【００６２】
図８に、単流熱交換器モデルを概略的に示す。図中の推定位置に、適合用の排気温度センサを取り付けて、実験データを収集する。適合用の排気温度センサは、試験用の内燃機関１にのみ装着されて出口ガス温度算出部５２、放熱量算出部５３、排気温度推定部５４等に含まれる各種定数の適合に用いられる。生産された内燃機関１では、出口ガス温度算出部５２、放熱量算出部５３、排気温度推定部５４により推定された排気ガス温度Ｔｅｘを用いて各種制御が行われる。
【００６３】
また、燃焼室出口及び排気管入口は、燃焼室と排気ポートとの境界部分に相当し、排気バルブ１５が配置された位置に相当する。燃焼室出口から推定位置（適合用の排気温度センサの取付位置）までを円筒と仮定する。燃焼室出口から推定位置までの排気管の内部表面積が、排気ガスから排気管に伝熱される全伝熱面積Ａ０［ｍ２］になる。なお、複数気筒の場合は、全気筒の燃焼室出口から推定位置までの全内部表面積Ａ０を、気筒数で割ったものを用いることができる。排気管を流れる排気ガス流量をＱｅｘ［ｇ／ｓ］とする。排気管の外部は外気（排気管近傍の空気）により冷却されており、その外気温Ｔａ［Ｋ］は、一定であると仮定する。
【００６４】
次に、燃焼室出口の排気ガスの温度である出口ガス温度をＴｏｕｔ［Ｋ］とし、この出口ガス温度Ｔｏｕｔと外気温Ｔａとの差をθ１とする。推定位置における排気ガス温度をＴｅｘ［Ｋ］とし、排気ガス温度Ｔｅｘと外気温Ｔａとの差をθ２とする。なお、排気ガス温度Ｔｅｘは、排気温度センサの応答遅れの影響がない温度であり、排気ガス温度の瞬時値となる。また、全伝熱面積Ａ０のうち燃焼室出口から任意の位置までの面積をＡ［ｍ２］とし、その任意の位置における瞬時排気温度をＴ［Ｋ］とし、任意の位置の微小伝熱面積ｄＡにおける排気ガス温度の変化をｄＴとする。瞬時排気温度Ｔと外気温Ｔａとの差をθとすると、ｄθ＝ｄＴである。さて、この微小伝熱面積ｄＡで単位時間に交換される熱量ｄＱは、熱貫流率（熱通過率）Ｋｈｔ［Ｗ／（ｍ２・Ｋ）］を用いて表すと次式となる。
【数１１】
2019132152.tif 000013
【００６５】
排気ガス流量Ｑｅｘ［ｇ／ｓ］と排気ガスの比熱Ｃｅｘ［Ｊ／（ｇ・Ｋ）］の積を水当量という。この水当量の排気ガスがｄＱ［Ｊ］の熱を失って、ｄＴだけ温度が下がるので、次式で表せる。
【数１２】
2019132152.tif 000014
【００６６】
（１１）式、（１２）式から、ｄＱを削除して整理し、積分すると次式となる。ここで、Ｃｏｎｓｔは、積分定数である。
【数１３】
2019132152.tif 000015
【００６７】
燃焼室出口では、Ａ＝０、θ＝θ１であり、推定位置では、Ａ＝Ａ０、θ＝θ２であるので、式（１３）に適用して、式変形すると次式を得る。
【数１４】
2019132152.tif 000016
【００６８】
この排気管において外気に奪われる熱量Ｑ［Ｊ］は次式で表される。
【数１５】
2019132152.tif 000017
【００６９】
また、この排気管において外気に奪われる熱量の最大値Ｑｍａｘは、θ２＝０の場合（排気ガス温度Ｔｅｘが外気温Ｔａまで冷却された場合）なので、次式となる。
【数１６】
2019132152.tif 000018
【００７０】
（１４）式、（１５）式、（１６）式から、排気管の温度効率ηは次式で表される。
【数１７】
2019132152.tif 000019
【００７１】
（１７）式において、熱貫流率Ｋｈｔ、燃焼室出口から推定位置までの全伝熱面積Ａ０、及び排気ガスの比熱Ｃｅｘが、一定値であると仮定とすると、排気管の温度効率ηは排気ガス流量Ｑｅｘの関数であることがわかる。結局、排気ガス流量Ｑｅｘに対応する排気管の温度効率η［％］、外気温Ｔａ、出口ガス温度Ｔｏｕｔがわかれば、排気ガス温度Ｔｅｘは、次式により推定することができる。
【数１８】
2019132152.tif 000020
【００７２】
（１８）式の第２式の右辺第２項が、排気系に単流熱交換器モデルを適用した場合の温度低下量ΔＴｄを表すこととなる。以上により、出口ガス温度Ｔｏｕｔ、外気温Ｔａ、排気管の温度効率ηから排気ガス温度Ｔｅｘを算出する方法が示された。
【００７３】
＜放熱量算出部５３＞
放熱量算出部５３は、燃焼室の出口から推定位置までの排気管の放熱による排気ガスの温度低下量ΔＴｄを算出する。上述したように、排気管を単流熱交換器にモデル化して導出した（１７）式から、排気管の温度効率ηは、排気ガス流量Ｑｅｘの関数であることがわかり、（１８）式から、排気管の温度効率ηに基づいて、温度低下量ΔＴｄを算出できることがわかる。そこで、放熱量算出部５３は、排気管内の排気ガス流量Ｑｅｘに基づいて、排気管内の排気ガスを加熱流体とし、排気管外の空気を受熱流体とした単流熱交換器としての排気管の温度効率ηを算出し、排気管の温度効率ηに基づいて温度低下量ΔＴｄを算出する。
【００７４】
放熱量算出部５３は、（１７）式と同様の次式を用いて、排気ガス流量Ｑｅｘに基づいて、排気管の温度効率ηを算出する。
【数１９】
2019132152.tif 000021
【００７５】
Ｋηは、演算定数である。演算定数Ｋηは、排気管の熱貫流率Ｋｈｔに、燃焼室出口から推定位置までの全伝熱面積Ａ０を乗算し、排気ガスの比熱Ｃｅｘで除算することにより設定できる。排気管の熱貫流率Ｋｈｔは、適合値であるが、例えば、１０〜１５［ｗ／（ｍ２・Ｋ）］程度の値になる。全伝熱面積Ａ０は、排気管の構造から算出できる。排気ガスの比熱Ｃｅｘには、上述した値が用いられる。演算定数Ｋηは、固定値とされてもよいが、空燃比ＡＦに応じて変化する排気ガスの比熱Ｃｅｘに応じて変化されてもよい。また、演算定数Ｋηは、実験による適合値とされてもよい。
【００７６】
或いは、放熱量算出部５３は、排気ガス流量Ｑｅｘと排気管の温度効率ηとの関係が予め設定された温度効率特性データを参照し、現在の排気ガス流量Ｑｅｘに対応する排気管の温度効率ηを算出するように構成されてもよい。温度効率特性データは、（１９）式の特性を表すものであるが、実験により適合されてもよい。例えば、複数の排気ガス流量Ｑｅｘの動作点において測定した、適合用の排気温度センサを用いた出口ガス温度Ｔｏｕｔの測定値、適合用の排気温度センサを用いた推定位置の排気ガス温度Ｔｅｘ０の測定値、外気温Ｔａの測定値に基づいて、（１８）式の第１式を用い、排気管の温度効率ηを算出する。そして、複数の動作点の排気ガス流量Ｑｅｘ及び排気管の温度効率ηの実験データを近似して、温度効率特性データを設定する。同様の方法により、演算定数Ｋηが設定されてもよい。
【００７７】
放熱量算出部５３は、（１８）式の第２式と同様の次式に示すように、出口ガス温度算出部５２により算出された出口ガス温度Ｔｏｕｔから外気温Ｔａを減算した値に、排気管の温度効率ηを乗算した値を、排気ガスの温度低下量ΔＴｄとして算出する。
【数２０】
2019132152.tif 000022
【００７８】
＜排気温度推定部５４＞
排気温度推定部５４は、次式に示すように、出口ガス温度Ｔｏｕｔから温度低下量ΔＴｄを減算して推定位置の排気ガス温度Ｔｅｘを推定する。
【数２１】
2019132152.tif 000023
【００７９】
推定位置は、後述する排気温度利用制御部５５において必要になる排気ガス温度の位置に設定される。例えば、推定位置は、触媒１９の上流側の位置、排気管１７とＥＧＲ流路２１との接続位置、排気管に過給機が設けられる場合はタービンの上流側の位置、排気管に排気ガスセンサが設けられる場合は排気ガスセンサの位置等に設定される。推定位置に合わせて、温度効率ηの算出に必要な演算定数Ｋη又は温度効率特性データが設定される。
【００８０】
排気温度推定部５４は、複数の推定位置の排気ガス温度を推定するように構成されてもよい。この場合は、放熱量算出部５３は、推定位置に応じて演算定数Ｋη又は温度効率特性データの設定値を切り替えて、各推定位置の温度効率ηを算出し、各推定位置の温度低下量ΔＴｄを算出する。そして、排気温度推定部５４は、各推定位置の温度低下量ΔＴｄを用いて、各推定位置の排気ガス温度Ｔｅｘを推定する。このように、推定位置に応じて演算定数Ｋη又は温度効率特性データを切り替えるだけで、容易に複数の推定位置の排気ガス温度を推定することができる。
【００８１】
後述するように、排気管に排気ガス温度センサが設けられ、排気ガス温度センサの異常診断を行う場合は、排気ガス温度センサにより検出した排気ガス温度と、推定した排気ガス温度Ｔｅｘとを比較する。しかし、排気ガス温度センサにより検出した排気ガス温度には、センサの熱容量等による応答遅れが生じる。或いは、排気ガス温度には、排気管の熱容量による応答遅れが生じる。そこで、排気温度推定部５４は、排気ガス温度Ｔｅｘに対して、応答遅れ処理を行うことにより、応答遅れ処理後の排気ガス温度Ｔｅｘｆｔを算出する。例えば、排気温度推定部５４は、次式に示す一次遅れフィルタ処理を行って、応答遅れ処理後の排気ガス温度Ｔｅｘｆｔを算出する。ここで、フィルタ定数Ｋｆは、センサの時定数τと演算周期Δｔから設定される。（ｎ）は、今回の演算周期の値であることを表し、（ｎ−１）は、前回の演算周期の値であることを表す。
【数２２】
2019132152.tif 000024
【００８２】
＜排気温度利用制御部５５＞
排気温度利用制御部５５は、排気ガス温度の推定値Ｔｅｘを用いた、排気温度制御、バルブ流量特性演算、排気温度センサ異常診断、及びタービン出力演算のいずれか１つ以上を実行する。
【００８３】
排気温度制御は、排気ガス温度の推定値Ｔｅｘを用いて排気ガス温度の制御を行う処理である。排気温度推定部５４は、排気ガス温度の制御を行う位置を推定位置に設定して、排気ガス温度の推定を行う。例えば、排気温度利用制御部５５は、排気ガスの温度を低下させるためのエンリッチ制御を行っている場合に、排気ガス温度の推定値Ｔｅｘが目標温度に近づくように、燃料噴射のエンリッチ量を変化させる。また、排気温度利用制御部５５は、排気ガス温度の推定値Ｔｅｘが目標温度に近づくように、点火時期ＳＡを変化させたり、排気行程の燃料噴射量を変化させたりする。
【００８４】
バルブ流量特性演算は、排気ガス温度の推定値Ｔｅｘを用いて排気ガスが流れるバルブの流量特性を算出する処理である。排気温度推定部５４は、バルブの上流側の位置を推定位置に設定して、排気ガス温度の推定を行う。排気温度利用制御部５５は、排気ガス温度の推定値Ｔｅｘに基づいて、バルブの流量特性としてバルブの上流側の排気ガスの音速及び密度を算出する。排気ガスが流れるバルブは、ＥＧＲバルブ２２、過給機のタービンを迂回するウェイストゲートバルブ等とされる。排気温度利用制御部５５は、ＥＧＲバルブ２２の流量特性を用いて、ＥＧＲ流量を算出し、ＥＧＲ量Ｑｃｅ、ＥＧＲ率Ｒｅｇｒを算出する。排気温度利用制御部５５は、ウェイストゲートバルブの流量特性を用いて、ウェイストゲートバルブを通過するバイパス流量を算出し、排気ガス流量からバイパス流量を減算して、タービン通過流量を算出し、タービン通過流量はタービン出力に比例するため、タービン通過流量を用いてタービの出力を算出する。タービン出力は、過給圧の制御に用いられる。
【００８５】
排気温度センサ異常診断は、排気ガス温度の推定値Ｔｅｘを用いて排気管に設けられた排気ガス温度センサの異常診断を行う処理である。排気ガス温度センサは、触媒、微粒子捕集フィルタ等の排気ガスの浄化装置に流入する排気ガス温度を管理するために設けられる。排気温度推定部５４は、排気ガス温度センサの位置を推定位置に設定して、排気ガス温度の推定を行う。排気温度利用制御部５５は、応答遅れ処理後の排気ガス温度の推定値Ｔｅｘｆｔと、排気ガス温度センサによる排気ガス温度の検出値とを比較し、両者の差が大きい場合に、排気ガス温度センサに異常が発生したと判定する。
【００８６】
タービン出力演算は、排気ガス温度の推定値Ｔｅｘを用いて排気管に設けられた過給機のタービン出力を算出する演算である。排気温度推定部５４は、タービンの上流側の位置を推定位置に設定して、排気ガス温度の推定を行う。排気温度利用制御部５５は、タービンに流入する排気ガス温度はタービン出力に比例するため、排気ガス温度の推定値Ｔｅｘを用いてタービン出力を算出する。
【００８７】
２−４．フローチャート
本実施の形態に係る制御装置５０の概略的な処理の手順（内燃機関１の制御方法）について、図９に示すフローチャートに基づいて説明する。図９のフローチャートの処理は、演算処理装置９０が記憶装置９１に記憶されたソフトウェア（プログラム）を実行することにより、所定の演算周期毎に繰り返し実行される。
【００８８】
ステップＳ０１で、運転状態検出部５１は、上記のように、内燃機関１の各種の運転状態を検出する運転状態検出処理（運転状態ステップ）を実行する。
【００８９】
ステップＳ０２で、燃焼室ガス量算出部５２ａは、上記のように、燃焼室に流入したガス量である燃焼室ガス量Ｑａｌｌを算出する燃焼室ガス量算出処理（燃焼室ガス量算出ステップ）を実行する。ステップＳ０３で、燃焼熱量算出部５２ｂは、上記のように、燃料の燃焼により燃焼室で発生した燃焼熱量ｈｌ＿ｂを算出する燃焼熱量算出処理（燃焼熱量算出ステップ）を実行する。
【００９０】
ステップＳ０４で、図示熱効率算出部５２ｄは、上記のように、燃焼熱量ｈｌ＿ｂのうち、燃焼室内の筒内圧による仕事として取り出せる熱量の割合である図示熱効率ηｉを算出する図示熱効率算出処理（図示熱効率算出ステップ）を実行する。また、図示熱効率算出部５２ｄは、上記のように、点火時期ＳＡの操作による図示熱効率の変化量Δηｉを算出する。本実施の形態では、図示熱効率算出部５２ｄは、上記にように、予め設定された基本点火時期ＳＡ０からの点火時期の遅角による図示熱効率の変化量Δηｉを算出する。本実施の形態では、基本点火時期ＳＡ０は、トルクが最大になる点火時期（ＭＢＴ）に設定されている。
【００９１】
ステップＳ０５で、冷却損失率算出部５２ｃは、上記のように、燃焼により発生した燃焼熱量ｈｌ＿ｂのうち、燃焼室の壁面に放熱される熱量の割合である冷却損失率ηｃを算出する冷却損失率算出処理（冷却損失率算出ステップ）を実行する。本実施の形態では、冷却損失率算出部５２ｃは、冷却損失率の基本値ηｃ０を算出し、点火時期ＳＡの操作による図示熱効率の変化量Δηｉに基づいて冷却損失率の基本値ηｃ０を補正して、最終的な冷却損失率ηｃを算出する。より具体的には、冷却損失率算出部５２ｃは、冷却損失率の基本値ηｃ０として、点火時期ＳＡが基本点火時期ＳＡ０に設定された場合の冷却損失率を算出する。冷却損失率算出部５２ｃは、基本点火時期ＳＡ０からの点火時期の遅角による図示熱効率の変化量Δηｉに基づいて、基本点火時期ＳＡ０からの点火時期の遅角による冷却損失率の変化量Δηｃを算出する。そして、冷却損失率算出部５２ｃは、冷却損失率の基本値ηｃ０を冷却損失率の変化量Δηｃにより変化させて、最終的な冷却損失率ηｃを算出する。
【００９２】
ステップＳ０６で、熱勘定温度算出部５２ｅは、上記のように、燃焼室ガス量Ｑａｌｌ、燃焼熱量ｈｌ＿ｂ、冷却損失率ηｃ、及び図示熱効率ηｉに基づいて、燃焼室の出口における排気ガスの温度である出口ガス温度Ｔｏｕｔを算出する熱勘定温度算出処理（熱勘定温度算出ステップ）を実行する。ステップＳ０７で、空燃比温度補正部５２ｆは、上記のように、空燃比ＡＦが理論空燃比よりリッチである場合に、エンリッチ量ΔＡＦｒに応じて、出口ガス温度Ｔｏｕｔの低下補正を行う空燃比温度補正処理（空燃比温度補正ステップ）を実行する。
【００９３】
ステップＳ０８で、放熱量算出部５３は、上記のように、燃焼室の出口から推定位置までの排気管の放熱による排気ガスの温度低下量ΔＴｄを算出する放熱量算出処理（放熱量算出ステップ）を実行する。本実施の形態では、放熱量算出部５３は、排気管内の排気ガス流量Ｑｅｘに基づいて、単流熱交換器としての排気管の温度効率ηを算出し、排気管の温度効率ηに基づいて温度低下量ΔＴｄを算出する。放熱量算出部５３は、（１９）式を用い、排気ガス流量Ｑｅｘに基づいて、排気管の温度効率ηを算出する。
【００９４】
ステップＳ０９で、排気温度推定部５４は、上記のように、出口ガス温度Ｔｏｕｔから温度低下量ΔＴｄを減算して推定位置の排気ガス温度Ｔｅｘを推定する排気温度推定処理（排気温度推定ステップ）を実行する。
【００９５】
ステップＳ１０で、排気温度利用制御部５５は、上記のように、排気ガス温度の推定値Ｔｅｘを用いた、排気温度制御、バルブ流量特性演算、排気温度センサ異常診断、及びタービン出力演算のいずれか１つ以上を行う排気温度利用制御処理（排気温度利用制御ステップ）を実行する。
【００９６】
〔その他の実施の形態〕
最後に、本開示のその他の実施の形態について説明する。なお、以下に説明する各実施の形態の構成は、それぞれ単独で適用されるものに限られず、矛盾が生じない限り、他の実施の形態の構成と組み合わせて適用することも可能である。
【００９７】
（１）上記の実施の形態１においては、内燃機関１は、ガソリンエンジンとされている場合を例として説明した。しかし、本開示の実施の形態はこれに限定されない。すなわち、内燃機関１は、ディーゼルエンジン、ＨＣＣＩ燃焼（Homogeneous-Charge Compression Ignition Combustion）を行うエンジン等の各種の内燃機関とされてもよい。この場合でも、制御装置５０は、筒内圧センサ又は振動センサ等により検出した、点火時期（燃焼開始時期）の操作（変化）に応じて、点火時期の操作（変化）による図示熱効率の変化量Δηｉを算出し、冷却損失率ηｃを補正するように構成されてもよい。
【００９８】
（２）上記の実施の形態１において、例示的に説明したように、内燃機関１は、過給機を備えてもよい。過給機は、排気管に設けられたタービンと、吸気管におけるスロットルバルブの上流側に設けられ、タービンと一体的に回転する圧縮機と、タービンを迂回するタービンバイパス通路に設けられたウェイストゲートバルブと、を備えている。また、内燃機関１は、触媒に加えて、微粒子捕集フィルタを備えてもよく、触媒、微粒子捕集フィルタの上流側に排気ガス温度センサを備えてもよい。
【００９９】
（３）上記の実施の形態１において、基本点火時期ＳＡ０は、トルクが最大になる点火時期（ＭＢＴ）に設定されている場合を例として説明した。しかし、本開示の実施の形態はこれに限定されない。すなわち、基本点火時期ＳＡ０は、ＭＢＴ以外の点火時期に設定されてもよい。例えば、ＭＢＴにおいて、ノックが発生する運転状態では、基本点火時期ＳＡ０は、ＭＢＴから遅角された点火時期設定される。
【０１００】
なお、本開示は、その開示の範囲内において、実施の形態を適宜、変形、省略したりすることが可能である。
【符号の説明】
【０１０１】
１  内燃機関、１７  排気管、２５  燃焼室、５０  内燃機関の制御装置、５１  運転状態検出部、５２  出口ガス温度算出部、５３  放熱量算出部、５４  排気温度推定部、５５  排気温度利用制御部、５２ａ  燃焼室ガス量算出部、５２ｂ  燃焼熱量算出部、５２ｃ  冷却損失率算出部、５２ｄ  図示熱効率算出部、５２ｅ  熱勘定温度算出部、５２ｆ  空燃比温度補正部、Ｋη  演算定数、Ｑａｌｌ  燃焼室ガス量、Ｑｅｘ  排気ガス流量、Ｒｅｇｒ  ＥＧＲ率、ＳＡ  点火時期、ＳＡ０  基本点火時期、Ｔａ  外気温、Ｔｅｘ  排気ガス温度、Ｔｏｕｔ  出口ガス温度、ｈｌ＿ｂ  燃焼熱量、ΔＴｄ  温度低下量、η  温度効率、ηｃ  冷却損失率、ηｃ０  冷却損失率の基本値、Δηｃ  冷却損失率の変化量、ηｅｘ  排気損失率、ηｉ  図示熱効率、ηｉ０  基本図示熱効率、Δηｉ  図示熱効率の変化量
【図１】
2019132152.tif 000025
【図２】
2019132152.tif 000026
【図３】
2019132152.tif 000027
【図４】
2019132152.tif 000028
【図５】
2019132152.tif 000029
【図６】
2019132152.tif 000030
【図７】
2019132152.tif 000031
【図８】
2019132152.tif 000032
【図９】
2019132152.tif 000033