(57)【要約】【課題】  簡易な計算により、内燃機関の筒内に吸入されている筒内ガスの着火前の温度（又は、圧力）を精度良く推定し得る内燃機関の筒内ガス温度（圧力）推定方法の提供。【解決手段】  この方法は、吸気酸素濃度ROXCを含む４つの因子の値が任意のクランク角度CAにおいて筒内ガス温度とそれぞれ線形関係にあることを利用する。即ち、４つの因子を対応する基準値に維持した状態での任意のCAにおける筒内ガス温度（基準筒内ガス温度Tabase(CA)）と、因子毎の任意のCAにおける因子の値の変化に対する筒内ガス温度の変化の割合である各傾き（gradTa_ROXC(CA)等）を予め実験等で求めておき、所望のCA、各因子の実際値と対応する前記基準値との偏差、及び前記傾きから筒内ガス補正量（DTa_ROXC(CA)等）を因子毎に求め、前記Tabase(CA)に各因子についての同筒内ガス補正量を加えた値を同所望のCAにおける筒内ガス温度Ta(CA)として求める。【選択図】  図３
2005090368.tif 000002
【特許請求の範囲】
【請求項１】
筒内に吸入されている筒内ガスの温度に影響を与える因子であってその値が任意のクランク角度において同筒内ガス温度と線形関係にある因子の値を取得する取得手段と、
前記因子の値が所定の基準値であるときの前記筒内ガス温度である基準筒内ガス温度を取得するための基準データをクランク角度に対応させながら予め記憶する基準データ記憶手段と、
前記因子の値の変化に対する前記筒内ガス温度の変化の程度を表す傾きを取得するための傾きデータをクランク角度に対応させながら予め記憶する傾きデータ記憶手段と、
を備えた内燃機関に適用される筒内ガス温度推定方法であって、
前記取得手段により前記因子の値を取得し、
前記取得した因子の値と、前記基準データと、前記傾きデータとに基づいて、任意のクランク角度における前記筒内ガス温度を推定する内燃機関の筒内ガス温度推定方法。
【請求項２】
請求項１に記載の筒内ガス温度推定方法において、
前記筒内ガス温度に影響を与える因子は、少なくとも、吸気管内のガス温度、前記筒内ガスの組成、同筒内ガスの圧縮開始時における筒内容積、及び燃焼室壁面温度を含み、
前記取得した因子の値と同因子に対応する前記所定の基準値との偏差に、同因子についての前記傾きデータから取得される任意のクランク角度における同因子についての前記傾きを乗じることで得られる同因子についての筒内ガス温度補正量を因子毎に求め、
前記基準データから取得される前記任意のクランク角度における前記基準筒内ガス温度に、前記各因子についての筒内ガス温度補正量をそれぞれ加算することで得られる値を同任意のクランク角度における前記筒内ガス温度として推定する内燃機関の筒内ガス温度推定方法。
【請求項３】
筒内に吸入されている筒内ガスの圧力に影響を与える因子であってその値が任意のクランク角度において同筒内ガス圧力と線形関係にある因子の値を取得する取得手段と、
前記因子の値が所定の基準値であるときの前記筒内ガス圧力である基準筒内ガス圧力を取得するための基準データをクランク角度に対応させながら予め記憶する基準データ記憶手段と、
前記因子の値の変化に対する前記筒内ガス圧力の変化の程度を表す傾きを取得するための傾きデータをクランク角度に対応させながら予め記憶する傾きデータ記憶手段と、
を備えた内燃機関に適用される筒内ガス圧力推定方法であって、
前記取得手段により前記因子の値を取得し、
前記取得した因子の値と、前記基準データと、前記傾きデータとに基づいて、任意のクランク角度における前記筒内ガス圧力を推定する内燃機関の筒内ガス圧力推定方法。
【請求項４】
請求項３に記載の筒内ガス圧力推定方法において、
前記筒内ガス圧力に影響を与える因子は、少なくとも、吸気管内のガス温度、吸気管内のガス圧力、前記筒内ガスの組成、同筒内ガスの圧縮開始時における筒内容積、及び燃焼室壁面温度を含み、
前記取得した因子の値と同因子に対応する前記所定の基準値との偏差に、同因子についての前記傾きデータから取得される任意のクランク角度における同因子についての前記傾きを乗じることで得られる同因子についての筒内ガス圧力補正量を因子毎に求め、
前記基準データから取得される前記任意のクランク角度における前記基準筒内ガス圧力に、前記各因子についての筒内ガス圧力補正量をそれぞれ加算することで得られる値を同任意のクランク角度における前記筒内ガス圧力として推定する内燃機関の筒内ガス圧力推定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、内燃機関の筒内に吸入されている筒内ガスの着火前の温度を推定する内燃機関の筒内ガス温度推定方法、及び同筒内ガスの着火前の圧力を推定する内燃機関の筒内ガス圧力推定方法に関する。
【背景技術】
【０００２】
火花点火式内燃機関、ディーゼル機関等の内燃機関から排出されるＮＯｘ等のエミッションの量は、着火後の火炎温度（燃焼温度）に強い相関を有している。従って、ＮＯｘ等のエミッションの排出量を低減するためには火炎温度を所定温度に制御することが有効である。一般に、この火炎温度は検出され得ないから、火炎温度を所定温度に制御するためには同火炎温度を推定する必要がある。一方、火炎温度は、着火前の筒内ガス（即ち、新気、ＥＧＲガス等）の温度、又は圧力（本願においては、単に「筒内ガス温度」、又は「筒内ガス圧力」と云うこともある。）により変化する。よって、火炎温度を推定するためには筒内ガス温度、又は筒内ガス圧力を推定することが効果的である。
【０００３】
また、特に、圧縮による自己着火により混合気が燃焼を開始するディーゼル機関においては、機関の運転状態に応じて着火時期を適切に制御する必要がある。この着火時期は着火前の筒内ガス温度、又は筒内ガス圧力に大きく依存する。従って、着火時期を適切に制御するためにも筒内ガス温度、又は筒内ガス圧力を推定する必要がある。
【０００４】
このような観点に基づき、下記特許文献１に記載のディーゼル機関の制御装置は、圧縮行程（及び膨張行程）における着火前の筒内ガスの状態（即ち、温度、圧力）が、吸気弁閉弁により同筒内ガスが密閉された後に断熱変化するとの仮定のもと、ポリトロープ指数を用いた断熱変化を表す一般的な式に基づいて任意のクランク角度における筒内ガス温度、及び筒内ガス圧力を推定している。
【特許文献１】特開２００１−２０７６５号公報
【０００５】
上記装置においては、筒内ガス温度、及び筒内ガス圧力を推定する際に使用されるポリトロープ指数の値を一定としている。しかしながら、実際には、ポリトロープ指数の値は筒内ガスの組成、或いは、時々刻々と変化する筒内ガス温度等に応じて変化する。従って、筒内ガス温度、及び筒内ガス圧力が精度良く推定され得ないという問題があった。
【０００６】
一方、上記断熱変化を表す一般的な式に基づいて筒内ガス温度、及び筒内ガス圧力を精度良く推定するため、ポリトロープ指数の値を筒内ガスの組成、或いは、時々刻々と変化する筒内ガス温度等に応じて可変とすることも考えられる。しかしながら、この場合、計算が煩雑となってマイクロコンピュータ（ＣＰＵ）の演算負荷が大きくなってしまう。
【発明の開示】
【０００７】
本発明は、かかる課題に対処するためになされたものであって、その目的は、簡易な計算により、内燃機関の筒内に吸入されている筒内ガスの着火前の温度、又は圧力を精度良く推定し得る内燃機関の筒内ガス温度推定方法、又は筒内ガス圧力推定方法を提供することにある。
【０００８】
筒内ガスの温度、又は筒内ガスの圧力に影響を与える因子には様々なものがある。本願出願人は、種々の実験等を重ねた結果、これらの因子の中で、任意のクランク角度において（少なくとも燃料噴射時期、及び／又は着火時期を含んだ圧縮行程、或いは膨張行程における所定の範囲内のいずれのクランク角度においても）その値が着火前の筒内ガス温度、又は筒内ガス圧力と線形関係にあるものが存在することを見出した。本発明は、この線形関係を利用した簡易な計算を行うことで、任意の（所望の）クランク角度における着火前の筒内ガス温度、又は筒内ガス圧力を推定する方法を提供する。
【０００９】
本発明による内燃機関の筒内ガス温度推定方法は、取得手段と、基準データ記憶手段と、傾きデータ記憶手段を備えた内燃機関に適用される。以下、これら各手段について順に説明する。
【００１０】
取得手段は、上記のようにその値が任意のクランク角度において筒内ガス温度と線形関係にある因子の値を取得する手段である。取得手段は、前記因子の値を所定の計算により推定する手段であっても、センサ等により物理的に検出する手段であってもよい。この取得手段は、取得すべき因子が複数種類存在する場合、各因子の値をそれぞれ取得する。
【００１１】
基準データ記憶手段は、前記因子の値が所定の基準値であるときの前記筒内ガス温度である基準筒内ガス温度を取得するための基準データをクランク角度に対応させながら予め記憶する手段である。より具体的には、基準データ記憶手段は、前記因子の実際の値を所定の基準値に維持した状態（前記因子が複数種類存在する場合、各因子の実際の値を対応する所定の各基準値にそれぞれ維持した状態）で内燃機関を運転させたときに得られる実際の筒内ガス温度（即ち、前記基準筒内ガス温度）に関するデータ（即ち、前記基準データ）を、実験等を通して、クランク角度に対応させながら（例えば、所定のクランク角度間隔毎に）予め取得しておき、これらのデータをマイクロコンピュータのＲＯＭ等に記憶した手段である。この基準データから任意の（所望の）クランク角度における基準筒内ガス温度が取得され得る。
【００１２】
傾きデータ記憶手段は、前記因子の値の変化に対する前記筒内ガス温度の変化の程度を表す傾きを取得するための傾きデータをクランク角度に対応させながら予め記憶する手段である。上述したように、上記因子の値は、任意のクランク角度において（少なくとも燃料噴射時期、及び／又は着火時期を含んだ圧縮行程、或いは膨張行程における所定の範囲内のいずれのクランク角度においても）筒内ガス温度と線形関係にあるから、同任意のクランク角度において、同因子の値の変化に対する筒内ガス温度の変化の程度を表す傾きが存在する（因子が複数種類存在する場合、因子毎に各因子についての傾きが存在する）。
【００１３】
即ち、傾きデータ記憶手段は、前記因子の実際の値を変化させたときの実際の筒内ガス温度の変化から得られる前記傾きに関するデータ（即ち、傾きデータ。前記因子が複数種類存在する場合、残りの因子をそれぞれ各所定値（例えば、前記対応する所定の各基準値）に維持した状態で各因子の実際の値を一つずつ変化させていったときの実際の筒内ガス温度の変化から因子毎に得られる同各因子についての前記傾きに関するデータ）を、実験等を通して、クランク角度に対応させながら（例えば、所定のクランク角度間隔毎に）予め取得しておき、これらのデータをマイクロコンピュータのＲＯＭ等に記憶した手段である。この傾きデータから任意の（所望の）クランク角度における傾き（前記因子が複数種類存在する場合、各因子についての因子毎の傾き）が取得され得る。
【００１４】
このような各手段を備えた内燃機関に適用される本発明による筒内ガス温度推定方法は、前記取得手段により前記因子の値を取得し、前記取得した因子の値（因子が複数種類存在する場合、それぞれ取得した各因子の値）と、前記基準データと、前記傾きデータ（因子が複数種類存在する場合、各因子についてのそれぞれの傾きデータ）とに基づいて、任意の（所望の）クランク角度における前記筒内ガス温度を推定する方法である。
【００１５】
これによれば、上記線形関係を利用することにより、前記基準データと前記傾きデータを予め実験等を通して準備しておくだけで、簡易な計算により、任意の（所望の）クランク角度における着火前の筒内ガス温度を精度良く推定することができる。
【００１６】
より具体的には、上記本発明による筒内ガス温度推定方法において、前記筒内ガス温度に影響を与える因子は、少なくとも、吸気管内のガス温度、前記筒内ガスの組成、同筒内ガスの圧縮開始時における筒内容積、及び燃焼室壁面温度を含み、各因子の値をそれぞれ取得するとともに、前記取得した因子の値と同因子に対応する前記所定の基準値との偏差に、同因子についての前記傾きデータから取得される任意の（所望の）クランク角度における同因子についての前記傾きを乗じることで得られる同因子についての筒内ガス温度補正量を因子毎に求め、前記基準データから取得される前記任意の（所望の）クランク角度における前記基準筒内ガス温度に、前記各因子についての筒内ガス温度補正量をそれぞれ加算することで得られる値を同任意の（所望の）クランク角度における前記筒内ガス温度として推定することが好適である。
【００１７】
本発明による内燃機関の筒内ガス圧力推定方法は、上記本発明による内燃機関の筒内ガス温度推定方法と同様であるので、その説明を省略する。この場合、前記取得手段により取得すべき因子（即ち、その値が任意のクランク角度において筒内ガス圧力と線形関係にある因子）としては、上記本発明による内燃機関の筒内ガス温度推定方法に使用される上記列挙した、吸気管内のガス温度、前記筒内ガスの組成、同筒内ガスの圧縮開始時における筒内容積、及び燃焼室壁面温度に加え、吸気管内のガス圧力が挙げられる。
【発明を実施するための最良の形態】
【００１８】
以下、本発明の実施形態に係る内燃機関の筒内ガス温度推定方法を実施する内燃機関（ディーゼル機関）の制御装置について図面を参照しつつ説明する。
【００１９】
図１は、係る内燃機関の制御装置を４気筒内燃機関（ディーゼル機関）１０に適用したシステム全体の概略構成を示している。このシステムは、燃料供給系統を含むエンジン本体２０、エンジン本体２０の各気筒の燃焼室（筒内）にガスを導入するための吸気系統３０、エンジン本体２０からの排ガスを放出するための排気系統４０、排気還流を行うためのＥＧＲ装置５０、及び電気制御装置６０を含んでいる。
【００２０】
エンジン本体２０の各気筒の上部には燃料噴射弁（噴射弁、インジェクタ）２１が配設されている。各燃料噴射弁２１は、図示しない燃料タンクと接続された燃料噴射用ポンプ２２に燃料配管２３を介して接続されている。燃料噴射用ポンプ２２は、電気制御装置６０と電気的に接続されていて、同電気制御装置６０からの駆動信号（後述する最終燃料噴射圧力Pcrfinに応じた指令信号）により燃料の実際の噴射圧力（吐出圧力）が同最終燃料噴射圧力Pcrfinになるように同燃料を昇圧するようになっている。
【００２１】
これにより、燃料噴射弁２１には、燃料噴射用ポンプ２２から前記最終燃料噴射圧力Pcrfinまで昇圧された燃料が供給されるようになっている。また、燃料噴射弁２１は、電気制御装置６０と電気的に接続されていて、同電気制御装置６０からの駆動信号（後述する燃料噴射量qfinに応じた指令信号）により所定時間だけ開弁し、これにより各気筒の燃焼室内に前記最終燃料噴射圧力Pcrfinにまで昇圧された燃料を前記燃料噴射量qfinだけ直接噴射するようになっている。
【００２２】
吸気系統３０は、エンジン本体２０の各気筒の燃焼室にそれぞれ接続された吸気マニホールド３１、吸気マニホールド３１の上流側集合部に接続され同吸気マニホールド３１とともに吸気通路を構成する吸気管３２、吸気管３２内に回動可能に保持されたスロットル弁３３、電気制御装置６０からの駆動信号に応答してスロットル弁３３を回転駆動するスロットル弁アクチュエータ３３ａ、スロットル弁３３の上流において吸気管３２に順に介装されたインタクーラー３４と過給機３５のコンプレッサ３５ａ、及び吸気管３２の先端部に配設されたエアクリーナ３６とを含んでいる。
【００２３】
排気系統４０は、エンジン本体２０の各気筒にそれぞれ接続された排気マニホールド４１、排気マニホールド４１の下流側集合部に接続された排気管４２、排気管４２に配設された過給機３５のタービン３５ｂ、及び排気管４２に介装されたディーゼルパティキュレートフィルタ（以下、「ＤＰＮＲ」と称呼する。）４３を含んでいる。排気マニホールド４１及び排気管４２は排気通路を構成している。
【００２４】
ＤＰＮＲ４３は、コージライト等の多孔質材料から形成されたフィルタ４３ａを備え、通過する排気ガス中のパティキュレートを細孔表面にて捕集するフィルタである。ＤＰＮＲ４３は、担体としてのアルミナに、カリウムＫ，ナトリウムＮａ，リチウムＬｉ，セシウムＣｓのようなアルカリ金属、バリウムＢａ，カルシウムＣａのようなアルカリ土類金属、及びランタンＬａ、イットリウムＹのような希土類金属から選ばれた少なくとも一つを白金とともに担持し、ＮＯｘを吸収した後に同吸収したＮＯｘを放出して還元する吸蔵還元型ＮＯｘ触媒としても機能するようになっている。
【００２５】
ＥＧＲ装置５０は、排気ガスを還流させる通路（ＥＧＲ通路）を構成する排気還流管５１と、排気還流管５１に介装されたＥＧＲ制御弁５２と、ＥＧＲクーラー５３とを備えている。排気還流管５１はタービン３５ｂの上流側排気通路（排気マニホールド４１）とスロットル弁３３の下流側吸気通路（吸気マニホールド３１）を連通している。ＥＧＲ制御弁５２は電気制御装置６０からの駆動信号に応答し、再循環される排気ガス量（排気還流量、ＥＧＲガス流量）を変更し得るようになっている。
【００２６】
電気制御装置６０は、互いにバスで接続されたＣＰＵ６１、ＣＰＵ６１が実行するプログラム、テーブル（ルックアップテーブル、マップ）、及び定数等を予め記憶したＲＯＭ６２、ＣＰＵ６１が必要に応じてデータを一時的に格納するＲＡＭ６３、電源が投入された状態でデータを格納するとともに同格納したデータを電源が遮断されている間も保持するバックアップＲＡＭ６４、並びにＡＤコンバータを含むインターフェース６５等からなるマイクロコンピュータである。
【００２７】
インターフェース６５は、空気流量（新気流量）計測手段であって吸気管３２に配置された熱線式エアフローメータ７１、スロットル弁３３の下流であって排気還流管５１が接続された部位よりも下流の吸気通路に設けられた吸気温センサ７２、スロットル弁３３の下流であって排気還流管５１が接続された部位よりも下流の吸気通路に配設された吸気管圧力センサ７３、クランクポジションセンサ７４、アクセル開度センサ７５、スロットル弁３３の下流であって排気還流管５１が接続された部位よりも下流の吸気通路に設けられた吸気酸素濃度センサ７６、及び、エンジン本体２０のシリンダボアライナー部に配設された壁面温度センサ７７と接続されていて、これらのセンサからの信号をＣＰＵ６１に供給するようになっている。また、インターフェース６５は、燃料噴射弁２１、燃料噴射用ポンプ２２、スロットル弁アクチュエータ３３ａ、及びＥＧＲ制御弁５２と接続されていて、ＣＰＵ６１の指示に応じてこれらに駆動信号を送出するようになっている。
【００２８】
熱線式エアフローメータ７１は、吸気通路内を通過する吸入空気の質量流量（単位時間当りの吸入空気量、単位時間あたりの新気量）を計測し、同質量流量Ga（空気流量Ga）を表す信号を発生するようになっている。吸気温センサ７２は、エンジン１０のシリンダ（即ち、燃焼室、筒内）に吸入されるガスの温度（即ち、吸気温度）を検出し、同吸気温度Tbを表す信号を発生するようになっている。吸気管圧力センサ７３は、エンジン１０のシリンダに吸入されるガスの圧力（即ち、吸気管圧力）を検出し、同吸気管圧力Pbを表す信号を発生するようになっている。
【００２９】
クランクポジションセンサ７４は、各気筒の絶対クランク角度を検出し、クランク角度CAを表すとともにエンジン１０の回転速度であるエンジン回転速度NEをも表す信号を発生するようになっている。アクセル開度センサ７５は、アクセルペダルＡＰの操作量を検出し、アクセル操作量Accpを表す信号を発生するようになっている。吸気酸素濃度センサ７６は、エンジン１０の気筒（即ち、筒内、燃焼室内）に吸入されるガス（従って、筒内ガス）中の酸素濃度（即ち、吸気酸素濃度）を検出し、同吸気酸素濃度ROXCを表す信号を発生するようになっている。壁面温度センサ７７は、燃焼室壁面温度を検出し、同燃焼室壁面温度Twallを表す信号を発生するようになっている。
【００３０】
（筒内ガス温度の推定方法の概要）
次に、上記のように構成された内燃機関の制御装置（以下、「本装置」と云うこともある。）による本発明の実施形態に係る筒内ガス温度の推定方法について説明する。図２は、或る一つの気筒のシリンダ内(筒内)に吸気マニホールド３１からガスが吸入され、筒内に吸入されたガス（即ち、筒内ガス）が排気マニホールド４１へ排出される様子を模式的に示した図である。
【００３１】
図２に示したように、筒内に吸入されるガス（従って、筒内ガス）には、吸気管３２の先端部からスロットル弁３３を介して吸入された新気と、排気還流管５１からＥＧＲ制御弁５２を介して吸入されたＥＧＲガスが含まれる。吸入される新気量（質量）と吸入されるＥＧＲガス量（質量）の和に対するＥＧＲガス量の割合（即ち、ＥＧＲ率）は、運転状態に応じて電気制御装置６０（ＣＰＵ６１）により適宜制御されるスロットル弁３３の開度、及びＥＧＲ制御弁５２の開度に応じて変化する。
【００３２】
かかる新気、及びＥＧＲガスは、吸気行程において開弁している吸気弁Vinを介してピストンの下降に伴って筒内に吸入されて筒内ガスとなる。筒内ガスは、機関１０の運転状態に応じて決定される吸気弁閉弁時にて吸気弁Vinが閉弁することにより筒内に密閉され、その後の圧縮行程においてピストンの上昇に伴って圧縮される。そして、ピストンが上死点近傍に達っすると（具体的には、後述する最終燃料噴射時期finjfinが到来すると）、本装置は、前記燃料噴射量qfinに応じた所定時間だけ燃料噴射弁２１を開弁することで燃料を筒内に直接噴射する。この結果、噴射された燃料は、時間の経過に伴って同筒内ガスと混ざり合いながら混合気となって筒内において円錐状に拡散していき、所定のタイミングにて自己着火が発生することに起因して燃焼する。そして、燃焼後の筒内ガスは、排気行程において開弁している排気弁Voutを介してピストンの上昇に伴って排気マニホールド４１へ排出されて排ガスとなる。
【００３３】
上記のように吸気弁閉弁時以降筒内に密閉される筒内ガスの着火前の温度に影響を与える因子には様々なものがある。本願出願人は、種々の実験等を重ねた結果、これらの因子の中で、吸気弁閉弁時以降、着火前までの間のいずれのクランク角度CAにおいても、その値が筒内ガス温度と線形関係（或いは、略線形関係）にある因子が存在することを見出した。
【００３４】
係る因子としては、筒内ガスの組成としての筒内ガスの酸素濃度（即ち、吸気酸素濃度ROXC）、吸気管（吸気マニホールド３１）内のガス温度（即ち、吸気温度Tb）、同筒内ガスの圧縮開始時（以下、単に「圧縮開始時」と云うこともある。）における筒内容積Vastart、及び燃焼室壁面温度Twallの４つの因子が挙げられる。そして、本願出願人は、これらの各因子と筒内ガス温度との各線形関係をそれぞれ利用して、吸気弁閉弁時以降着火前の任意の（所望の）クランク角度CAにおける筒内ガス温度を推定する方法を見出した。
【００３５】
この方法を使用するには、上記吸気酸素濃度ROXC，吸気温度Tb，筒内容積Vastart，及び燃焼室壁面温度Twallがそれぞれ、対応する所定の基準値である、基準吸気酸素濃度ROXCbase，基準吸気温度Tbbase，基準筒内容積Vastartbase，及び基準燃焼室壁面温度Twallbaseに維持されている状態でのクランク角度CAにおける筒内ガス温度である基準筒内ガス温度Tabase(CA)と、因子毎のクランク角度CAにおける上記因子の値の変化に対する筒内ガス温度の変化の割合である傾き、即ち、後述する吸気酸素濃度についての傾きgradTa_ROXC(CA)，吸気温度についての傾きgradTa_Tb(CA)，圧縮開始時における筒内容積についての傾きgradTa_Vastart(CA)，及び燃焼室壁面温度についての傾きgradTa_Twall(CA)を予めそれぞれ求めておく必要がある。これらの求め方については後述する。
【００３６】
そして、この方法は、概略的に述べると、上記各因子の実際の値を対応するセンサ等からそれぞれ検出し、所望のクランク角度CAと、同検出値と対応する基準値との偏差と、上記傾きとから因子毎の筒内ガス温度補正量、即ち、後述する吸気酸素濃度についての筒内ガス温度補正量DTa_ROXC(CA)，吸気温度についての筒内ガス温度補正量DTa_Tb(CA)，圧縮開始時における筒内容積についての筒内ガス温度補正量DTa_Vastart(CA)，及び燃焼室壁面温度についての筒内ガス温度補正量DTa_Twall(CA)をそれぞれ求め、上記基準筒内ガス温度Tabase(CA)に上記各因子についての筒内ガス温度補正量をそれぞれ加算した値を上記所望のクランク角度CAにおける筒内ガス温度Ta(CA)として求める（推定する）方法である。
【００３７】
以下、先ず、係る方法に従って筒内ガス温度Ta(CA)を推定する際に必要となる基準筒内ガス温度Tabase(CA)と、吸気酸素濃度についての筒内ガス温度補正量DTa_ROXC(CA)について、任意のクランク角度CAにおける吸気酸素濃度ROXCと筒内ガス温度Taとの線形関係を示した図３を参照しながら説明する。
【００３８】
図３に示した曲線Ａは、上記４つの因子が総て対応する上記基準値に維持されている状態（例えば、図３に示したように、吸気酸素濃度ROXCは上記基準吸気酸素濃度ROXCbaseに維持されている状態）におけるクランク角度CAと基準筒内ガス温度Tabase(CA)との関係を示している。従って、例えば、クランク角度CAが圧縮行程においてピストンが上死点に達した状態に対応する角度である場合（以下、「ATDC0°」と称呼する。）について考えると、この場合における基準筒内ガス温度Tabase(ATDC0°)は点Ｂに示した値となる。
【００３９】
上述したように、任意のクランク角度CAにおいて吸気酸素濃度ROXCと筒内ガス温度Taとは線形関係にあるから、例えば、ATDC0°において、直線Ｃに示すように、吸気酸素濃度ROXCと筒内ガス温度Taとは或る傾きを持った線形関係にある。この直線Ｃの傾き（即ち、ATDC0°における吸気酸素濃度についての傾き）をgradTa_ROXC(ATDC0°)と表す。そうすると、吸気酸素濃度ROXCが吸気酸素濃度センサ７６に基づく検出値（検出吸気酸素濃度）ROXC0（図３を参照）であるとき（且つ、残りの３つの因子の値が総て対応する基準値に維持されているとき）、ATDC0°における筒内ガス温度Taは、点Ｄに示すように、ATDC0°における基準筒内ガス温度Tabase(ATDC0°)に、検出吸気酸素濃度ROXC0と基準吸気酸素濃度ROXCbaseとの偏差に直線Ｃの傾きを乗じた値「gradTa_ROXC(ATDC0°)・(ROXC0-ROXCbase)」を加えた値となる。この値「gradTa_ROXC(ATDC0°)・(ROXC0-ROXCbase)」を、ATDC0°における吸気酸素濃度についての筒内ガス温度補正量DTa_ROXC(ATDC0°)と表す。
【００４０】
同様に、例えば、ATDC-αについて考えると、この場合における基準筒内ガス温度Tabase(ATDC-α)は点Ｅに示した値となる。そして、ATDC-αにおいても、直線Ｆに示すように、吸気酸素濃度ROXCと筒内ガス温度Taとは或る傾きを持った線形関係にある。この直線Ｆの傾き（即ち、ATDC-αにおける吸気酸素濃度についての傾き）をgradTa_ROXC(ATDC-α)と表す。そうすると、吸気酸素濃度ROXCが前記検出吸気酸素濃度ROXC0であるとき（且つ、残りの３つの因子の値が総て対応する基準値に維持されているとき）、ATDC-αにおける筒内ガス温度Taは、点Ｇに示すように、ATDC-αにおける基準筒内ガス温度Tabase(ATDC-α)に、検出吸気酸素濃度ROXC0と基準吸気酸素濃度ROXCbaseとの偏差に直線Ｆの傾きを乗じた値「gradTa_ROXC(ATDC-α)・(ROXC0-ROXCbase)」を加えた値となる。この値「gradTa_ROXC(ATDC-α)・(ROXC0-ROXCbase)」を、ATDC-αにおける吸気酸素濃度についての筒内ガス温度補正量DTa_ROXC(ATDC-α)と表す。
【００４１】
即ち、吸気酸素濃度ROXCが前記検出吸気酸素濃度ROXC0であるとき（且つ、残りの３つの因子の値が総て対応する基準値に維持されているとき）、任意のクランク角度CAにおける吸気酸素濃度についての筒内ガス温度補正量DTa_ROXC(CA)は、図４から容易に理解できるように、下記(1)式に従って求められる。下記(1)式において、gradTa_ROXC(CA)は任意のクランク角度CAにおける吸気酸素濃度についての傾きである。
【００４２】
DTa_ROXC(CA)＝gradTa_ROXC(CA)・(ROXC0-ROXCbase)  ・・・(1)
【００４３】
そして、この任意のクランク角度CAにおける吸気酸素濃度についての筒内ガス温度補正量DTa_ROXC(CA)を用いれば、吸気酸素濃度ROXCが前記検出吸気酸素濃度ROXC0であるとき（且つ、残りの３つの因子の値が総て対応する基準値に維持されているとき）、任意のクランク角度CAにおける筒内ガス温度Ta(CA)は、図３の曲線Ｈに示すように、「Tabase(CA)＋DTa_ROXC(CA)」と表すことができる。以上、任意のクランク角度CAにおける基準筒内ガス温度Tabase(CA)と、任意のクランク角度CAにおける吸気酸素濃度についての筒内ガス温度補正量DTa_ROXC(CA)について説明した。
【００４４】
次に、その他の因子についての筒内ガス温度補正量について説明する。これらも吸気酸素濃度についての筒内ガス温度補正量DTa_ROXC(CA)と同様の考え方に基づいて求めることができる。即ち、先ず、吸気温度についての筒内ガス温度補正量DTa_Tb(CA)について説明すると、吸気温度Tbが吸気温センサ７２に基づく検出値（検出吸気温度）Tb0であるとき（且つ、残りの３つの因子の値が総て対応する基準値に維持されているとき）、任意のクランク角度CAにおける吸気温度についての筒内ガス温度補正量DTa_Tb(CA)は、図５から容易に理解できるように、下記(2)式に従って求められる。下記(2)式において、gradTa_Tb(CA)は任意のクランク角度CAにおける吸気温度についての傾きである。
【００４５】
DTa_Tb(CA)＝gradTa_Tb(CA)・(Tb0-Tbbase)  ・・・(2)
【００４６】
同様に、圧縮開始時における筒内容積についての筒内ガス温度補正量DTa_Vastart(CA)について説明すると、圧縮開始時における筒内容積Vastartが検出値（検出圧縮開始時筒内容積）Vastart0であるとき（且つ、残りの３つの因子の値が総て対応する基準値に維持されているとき）、任意のクランク角度CAにおける圧縮開始時における筒内容積についての筒内ガス温度補正量DTa_Vastart(CA)は、図６から容易に理解できるように、下記(3)式に従って求められる。下記(3)式において、gradTa_Vastart(CA)は任意のクランク角度CAにおける圧縮開始時における筒内容積についての傾きである。なお、筒内容積Vaは機関１０の設計諸元に基づいてクランク角度CAの関数Va(CA)として取得することができるから、運転状態に応じて決定される吸気弁閉弁時におけるクランク角度CAに基づいて圧縮開始時における筒内容積Vastartを検出することができる。
【００４７】
DTa_Vastart(CA)＝gradTa_Vastart(CA)・(Vastart0-Vastartbase)  ・・・(3)
【００４８】
同様に、燃焼室壁面温度についての筒内ガス温度補正量DTa_Twall(CA)について説明すると、燃焼室壁面温度Twallが壁面温度センサ７７に基づく検出値（検出燃焼室壁面温度）Twall0であるとき（且つ、残りの３つの因子の値が総て対応する基準値に維持されているとき）、任意のクランク角度CAにおける燃焼室壁面温度についての筒内ガス温度補正量DTa_Twall(CA)は、図７から容易に理解できるように、下記(4)式に従って求められる。下記(4)式において、gradTa_Twall(CA)は任意のクランク角度CAにおける燃焼室壁面温度についての傾きである。
【００４９】
DTa_Twall(CA)＝gradTa_Twall(CA)・(Twall0-Twallbase)  ・・・(4)
【００５０】
以上のようにして、各因子についての筒内ガス温度補正量を求めることができる。そして、この方法では、吸気酸素濃度ROXC，吸気温度Tb，筒内容積Vastart，及び燃焼室壁面温度Twallがそれぞれ、対応する前記検出値である、検出吸気酸素濃度ROXC0，検出吸気温度Tb0，検出圧縮開始時筒内容積Vastart0，及び検出燃焼室壁面温度Twall0となっているときにおける任意のクランク角度CAにおける筒内ガス温度Ta(CA)を、下記(5)式に基づいて求める。以上が、本発明の実施形態に係る筒内ガス温度の推定方法の概要である。
【００５１】
Ta(CA)＝Tabase(CA)+DTa_ROXC(CA)+DTa_Tb(CA)+DTa_Vastart(CA)+DTa_Twall(CA)・・・(5)
【００５２】
（基準筒内ガス温度Tabase(CA)の求め方）
上述したクランク角度CAにおける基準筒内ガス温度Tabase(CA)は、例えば、４つの気筒のうちのいずれか一つの気筒に筒内ガス温度を直接検出できる筒内温度センサを取り付け、所定の方法で上記４つの因子の実際の値を上記各基準値にそれぞれ維持させた状態で内燃機関１０を所定時間に渡って運転させ、同所定時間内に得られた多数の筒内ガス温度検出値（サンプル値）の平均値をクランク角度に対応させながら（例えば、クランク角度１°毎に）取得することで求めることができる。
【００５３】
（各因子についての傾きの求め方）
例えば、吸気酸素濃度についての傾きgradTa_ROXC(CA)は、例えば、上記と同様に筒内温度センサを取り付け、所定の方法で上記４つの因子の実際の値を上記各基準値にそれぞれ維持させた状態で内燃機関１０を運転させ、この状態から、吸気酸素濃度ROXCのみを基準吸気酸素濃度ROXCbaseから変化させたときの筒内ガス温度検出値の変化から得られる傾きをクランク角度に対応させながら（例えば、クランク角度１°毎に）取得することで求めることができる。残りの３つの因子についての傾きについても同様にして求めることができる。
【００５４】
（燃料噴射制御の概要）
上記筒内ガス温度推定方法を実施する本装置は、上記のようにして求めた基準筒内ガス温度Tabase(CA)、並びに、吸気酸素濃度についての傾きgradTa_ROXC(CA)、吸気温度についての傾きgradTa_Tb(CA)、圧縮開始時における筒内容積についての傾きgradTa_Vastart(CA)、及び燃焼室壁面温度についての傾きgradTa_Twall(CA)を取得するためのデータを、それぞれマップとしてＲＯＭ６２に記憶している。
【００５５】
本装置は、目標着火時期を所定時期（ATDCθref、例えば、ATDC10°）に設定し、燃料が噴射される対象となる気筒（以下、「燃料噴射気筒」と称呼する。）の吸気弁閉弁時が到来する毎に、同時点における吸気酸素濃度ROXC，吸気温度Tb，筒内容積Vastart，及び燃焼室壁面温度Twallを、それぞれ、検出吸気酸素濃度ROXC0，検出吸気温度Tb0，検出圧縮開始時筒内容積Vastart0，及び検出燃焼室壁面温度Twall0として検出するとともに、ATDCθrefにおける筒内ガス温度Ta(ATDCθref)を上記(1)〜(5)式に従って推定する。そして、本装置は、混合気の着火時期を目標着火時期ATDCθrefに一致させるため、ATDCθrefにおける筒内ガス温度Ta(ATDCθref)が所定の目標筒内ガス温度Tarefになるように、燃料噴射開始時期、及び燃料噴射圧力をフィードバック制御する。
【００５６】
具体的には、今回の燃料噴射気筒について取得された上記筒内ガス温度Ta(ATDCθref)の値が前記目標筒内ガス温度Tarefよりも高いとき、今回の燃料噴射気筒についての燃料噴射開始時期を基本燃料噴射時期よりも所定量だけ遅らせ、且つ、燃料噴射圧力を基本燃料噴射圧力よりも所定量だけ低くする。これにより、今回の燃料噴射気筒についての目標着火時期ATDCθrefにおける実際の筒内ガス温度が低くなるように制御され、この結果、今回の燃料噴射気筒の実際の着火時期が前記目標着火時期ATDCθrefに一致せしめられる。
【００５７】
一方、今回の燃料噴射気筒について取得された上記筒内ガス温度Ta(ATDCθref)の値が前記目標筒内ガス温度Tarefよりも高いとき、今回の燃料噴射気筒についての燃料噴射開始時期を基本燃料噴射時期よりも所定量だけ早め、且つ、燃料噴射圧力を基本燃料噴射圧力よりも所定量だけ高くする。これにより、今回の燃料噴射気筒についての目標着火時期ATDCθrefにおける実際の筒内ガス温度が高くなるように制御され、この結果、今回の燃料噴射気筒の実際の着火時期が前記目標着火時期ATDCθrefに一致せしめられる。以上が、燃料噴射制御の概要である。
【００５８】
（実際の作動）
次に、上記のように構成された内燃機関の制御装置の実際の作動について説明する。
<筒内ガス温度の計算>
ＣＰＵ６１は、図８にフローチャートにより示した上記目標着火時期ATDCθrefにおける筒内ガス温度の計算を行うためのルーチンを所定時間の経過毎に繰り返し実行するようになっている。従って、所定のタイミングになると、ＣＰＵ６１はステップ８００から処理を開始し、ステップ８０５に進んで燃料噴射気筒についての吸気弁閉弁時が到来したか否かを判定し、「Ｎｏ」と判定する場合、ステップ８９５に直ちに進んで本ルーチンを一旦終了する。
【００５９】
いま、燃料噴射気筒についての吸気弁閉弁時が到来したものとして説明を続けると、ＣＰＵ６１はステップ８０５にて「Ｙｅｓ」と判定してステップ８１０に進み、同吸気弁閉弁時における吸気酸素濃度ROXC，吸気温度Tb，筒内容積Vastart，及び燃焼室壁面温度Twallを、それぞれ、検出吸気酸素濃度ROXC0，検出吸気温度Tb0，検出圧縮開始時筒内容積Vastart0，及び検出燃焼室壁面温度Twall0として検出する。
【００６０】
次に、ＣＰＵ６１はステップ８１５に進んで、ＲＯＭ６２にそれぞれ記憶されている、吸気酸素濃度についての傾きgradTa_ROXC(CA)を取得するためのマップ、吸気温度についての傾きgradTa_Tb(CA)を取得するためのマップ、圧縮開始時における筒内容積についての傾きgradTa_Vastart(CA)を取得するためのマップ、及び燃焼室壁面温度についての傾きgradTa_Twall(CA)を取得するためのマップに基づいて、上記目標着火時期ATDCθrefにおける吸気酸素濃度についての傾きgradTa_ROXC、同目標着火時期ATDCθrefにおける吸気温度についての傾きgradTa_Tb、同目標着火時期ATDCθrefにおける圧縮開始時における筒内容積についての傾きgradTa_Vastart、及び同目標着火時期ATDCθrefにおける燃焼室壁面温度についての傾きgradTa_Twallをそれぞれ求める。
【００６１】
次いで、ＣＰＵ６１はステップ８２０に進んで、ＲＯＭ６２にそれぞれ記憶されている基準筒内ガス温度Tabase(CA)を取得するためのマップに基づいて、上記目標着火時期ATDCθrefにおける基準筒内ガス温度Tabaseを求める。
【００６２】
続いて、ＣＰＵ６１はステップ８２５に進んで、ステップ８１０にて検出した各因子についての検出値、及びステップ８１５にて求めた各因子についての傾きに基づき、上記(1)〜(4)式に従って、上記目標着火時期ATDCθrefにおける吸気酸素濃度についての筒内ガス温度補正量DTa_ROXC、同目標着火時期ATDCθrefにおける吸気温度についての筒内ガス温度補正量DTa_Tb、同目標着火時期ATDCθrefにおける圧縮開始時における筒内容積についての筒内ガス温度補正量DTa_Vastart、及び同目標着火時期ATDCθrefにおける燃焼室壁面温度についての筒内ガス温度補正量DTa_Twallをそれぞれ求める。
【００６３】
そして、ＣＰＵ６１ステップ８３０に進んで、ステップ８２０にて求めた上記目標着火時期ATDCθrefにおける基準筒内ガス温度Tabaseと、ステップ８２５にて求めた同目標着火時期ATDCθrefにおける各因子についての筒内ガス補正量とに基づいて、上記(5)式に従って、同目標着火時期ATDCθrefにおける筒内ガス温度Taを求めた後、ステップ８９５に進んで本ルーチンを終了する。このようにして、燃料噴射気筒の吸気弁閉弁時が到来する毎に、目標着火時期ATDCθrefにおける筒内ガス温度Taが推定される。この値は後述するルーチンにて使用され、この結果、燃料噴射開始時期、及び燃料噴射圧力がフィードバック制御される。
<燃料噴射量等の制御>
ＣＰＵ６１は、図９にフローチャートにより示した燃料噴射量、燃料噴射時期、及び燃料噴射圧力を制御するためのルーチンを所定時間の経過毎に繰り返し実行するようになっている。従って、所定のタイミングになると、ＣＰＵ６１はステップ９００から処理を開始し、ステップ９０５に進んでアクセル開度Accp、エンジン回転速度NE、及び図１０に示したテーブル（マップ）Mapqfinから燃料噴射量qfinを求める。テーブルMapqfinは、アクセル開度Accp及びエンジン回転速度NEと指令燃料噴射量qfinとの関係を規定するテーブルであり、ＲＯＭ６２内に格納されている。
【００６４】
次いで、ＣＰＵ６１はステップ９１０に進み、燃料噴射量qfin、エンジン回転速度NE、及び図１１に示したテーブルMapfinjbaseから基本燃料噴射時期finjbaseを決定する。テーブルMapfinjbaseは、燃料噴射量qfin及びエンジン回転速度NEと基本燃料噴射時期finjbaseとの関係を規定するテーブルであり、ＲＯＭ６２内に格納されている。
【００６５】
その後、ＣＰＵ６１はステップ９１５に進んで、燃料噴射量qfin、エンジン回転速度NE、及び図１２に示したテーブルMapPcrbaseから基本燃料噴射圧力Pcrbaseを決定する。テーブルMapPcrbaseは、燃料噴射量qfin及びエンジン回転速度NEと基本燃料噴射圧力Pcrbaseとの関係を規定するテーブルであり、ＲＯＭ６２内に格納されている。
【００６６】
次に、ＣＰＵ６１はステップ９２０に進み、前記目標筒内ガス温度Tarefから先の図８のルーチンのステップ８３０により求められている最新の（従って、今回の燃料噴射気筒についての）目標着火時期ATDCθrefにおける筒内ガス温度Taを減じた値を筒内ガス温度偏差ΔTaとして格納する。
【００６７】
続いて、ＣＰＵ６１はステップ９２５に進んで、噴射時期補正値Δθを、前記筒内ガス温度偏差ΔTaと図１３に示したテーブルMapΔθとから決定する。テーブルMapΔθは、筒内ガス温度偏差ΔTaと噴射時期補正値Δθとの関係を規定するテーブルであり、ＲＯＭ６２内に格納されている。
【００６８】
その後、ＣＰＵ６１はステップ９３０に進んで、噴射圧力補正値ΔPcrを、前記筒内ガス温度偏差ΔTaと図１４に示したテーブルMapΔPcrとから決定する。テーブルMapΔPcrは、混合気温度偏差ΔTmixと噴射圧力補正値ΔPcrとの関係を規定するテーブルであり、ＲＯＭ６２内に格納されている。
【００６９】
次いで、ＣＰＵ６１はステップ９３５に進み、基本噴射時期finjbaseを噴射時期補正値Δθで補正して最終燃料噴射時期finjfinを決定する。これにより、筒内ガス温度偏差ΔTaに応じて噴射時期が補正されることになる。この場合、図１３から明らかなように、筒内ガス温度偏差ΔTaが正の大きい値になるほど噴射時期補正値Δθが正の大きな値となって最終燃料噴射時期finjfinが進角側となり、同筒内ガス温度偏差ΔTaが負の大きい値（絶対値が大きい値）になるほど噴射時期補正値Δθは負の大きな値となって最終燃料噴射時期finjfinが遅角側に移行される。
【００７０】
続いて、ＣＰＵ６１はステップ９４０に進み、基本燃料噴射圧力Pcrbaseを噴射圧力補正値ΔPcrで補正して最終燃料噴射圧力Pcrfinを決定する。これにより、筒内ガス温度偏差ΔTaに応じて噴射圧力が補正されることになる。この場合、図１４から明らかなように、筒内ガス温度偏差ΔTaが正の大きい値になるほど噴射圧力補正値ΔPcrが正の大きな値となって最終燃料噴射圧力Pcrfinが高圧側となり、同筒内ガス温度偏差ΔTaが負の大きい値（絶対値が大きい値）になるほど噴射圧力補正値ΔPcrは負の大きな値となって最終燃料噴射圧力Pcrfinが低圧側に移行される。この結果、燃料噴射用ポンプ２２の吐出圧力が制御されることで、燃料噴射弁２１には前記決定された最終燃料噴射圧力Pcrfinまで昇圧された燃料が供給される。
【００７１】
そして、ＣＰＵ６１は、続くステップ９４５にて燃料噴射気筒についての噴射開始時期（即ち、上記決定された最終燃料噴射時期finjfin）が到来したか否かを判定し、「Ｙｅｓ」と判定する場合はステップ９５０に進んで、上記決定された燃料噴射量qfinの燃料を燃料噴射気筒についての燃料噴射弁２１から上記決定された最終燃料噴射圧力Pcrfinをもって噴射した後、ステップ９９５に進んで本ルーチンを一旦終了する。また、ステップ９４５にて「Ｎｏ」と判定される場合、直接ステップ９９５に進んで本ルーチンを一旦終了する。以上により、燃料噴射量、燃料噴射時期、及び燃料噴射圧力の制御が達成される。
【００７２】
以上、説明したように、本発明の実施形態に係る内燃機関の筒内ガス温度の推定方法によれば、上記基準筒内ガス温度Tabase(CA)、並びに、吸気酸素濃度についての傾きgradTa_ROXC(CA)、吸気温度についての傾きgradTa_Tb(CA)、圧縮開始時における筒内容積についての傾きgradTa_Vastart(CA)、及び燃焼室壁面温度についての傾きgradTa_Twall(CA)を取得するための各データを予め実験等を通して取得し、それぞれマップとしてＲＯＭ６２に記憶しておくことで、簡易な計算により（従って、ＣＰＵ６１の計算負荷を増大させることなく）、任意の（所望の）クランク角度CAにおける着火前の筒内ガス温度を精度良く推定することができた。
【００７３】
（筒内ガス圧力の推定方法の概要）
次に、上記のように構成された内燃機関の制御装置による本発明の実施形態に係る筒内ガス圧力の推定方法について説明する。本願出願人は、種々の実験等を重ねた結果、上述した筒内ガス温度の場合と同様、吸気弁閉弁時以降、着火前までの間のいずれのクランク角度CAにおいても、その値が筒内ガス圧力と線形関係（或いは、略線形関係）にある因子が存在することを見出した。
【００７４】
係る因子としては、上述した筒内ガス温度の場合において使用される吸気酸素濃度ROXC、吸気温度Tb、圧縮開始時における筒内容積Vastart、及び燃焼室壁面温度Twallの４つに加え、吸気管圧力Pbが挙げられる。そして、本願出願人は、これらの各因子と筒内ガス圧力との各線形関係をそれぞれ利用して、吸気弁閉弁時以降着火前の任意の（所望の）クランク角度CAにおける筒内ガス圧力を推定する方法を見出した。
【００７５】
この方法を使用するには、先に説明した筒内ガス温度の推定方法の場合と同様、上記吸気酸素濃度ROXC，吸気管圧力Pb、吸気温度Tb，筒内容積Vastart，及び燃焼室壁面温度Twallがそれぞれ、対応する所定の基準値である、基準吸気酸素濃度ROXCbase，基準吸気管圧力Pbbase，基準吸気温度Tbbase，基準筒内容積Vastartbase，及び基準燃焼室壁面温度Twallbaseに維持されている状態でのクランク角度CAにおける筒内ガス圧力である基準筒内ガス圧力Pabase(CA)と、因子毎のクランク角度CAにおける上記因子の値の変化に対する筒内ガス圧力の変化の割合である傾き、即ち、吸気酸素濃度についての傾きgradPa_ROXC(CA)，吸気管圧力についての傾きgradPa_Pb(CA)，吸気温度についての傾きgradPa_Tb(CA)，圧縮開始時における筒内容積についての傾きgradPa_Vastart(CA)，及び燃焼室壁面温度についての傾きgradPa_Twall(CA)を予めそれぞれ求めておく必要がある。これらは、先に紹介した基準筒内ガス温度Tabase(CA)の求め方、及び各因子についての傾きの求め方において、４つの気筒のうちのいずれか一つの気筒に筒内温度センサを取り付ける代わりに、同いずれか一つの気筒に筒内ガス圧力を直接検出できる筒内圧力センサを取り付けることで求めることができる。
【００７６】
そして、この方法は、先に説明した筒内ガス温度の推定方法の場合と同様、上記５つの因子の実際の値を対応するセンサ等からそれぞれ検出し、所望のクランク角度CAと、同検出値と対応する基準値との偏差と、上記傾きとから因子毎の筒内ガス圧力補正量、即ち、吸気酸素濃度についての筒内ガス圧力補正量DPa_ROXC(CA)，吸気温度についての筒内ガス圧力補正量DPa_Tb(CA)，吸気管圧力についての筒内ガス圧力補正量DPa_Pb(CA)，圧縮開始時における筒内容積についての筒内ガス圧力補正量DPa_Vastart(CA)，及び燃焼室壁面温度についての筒内ガス圧力補正量DPa_Twall(CA)をそれぞれ求め、上記基準筒内ガス圧力Pabase(CA)に上記各因子についての筒内ガス圧力補正量をそれぞれ加算した値を上記所望のクランク角度CAにおける筒内ガス圧力Pa(CA)として求める（推定する）方法である。係る筒内ガス圧力Pa(CA)の推定方法の詳細については、先に説明した筒内ガス温度Ta(CA)の推定方法と同様であるので、その説明を省略する。
【００７７】
この本発明の実施形態に係る筒内ガス圧力の推定方法を実施する内燃機関の制御装置は、先に説明した図８に示したルーチンに代えて図１５に示したルーチンを所定時間の経過毎に繰り返し実行することで、燃料噴射気筒についての吸気弁閉弁時が到来する毎に、上記目標着火時期ATDCθrefにおける筒内ガス圧力Paを推定する。図１５のステップ１５００〜１５３０は、それぞれ、図８のステップ８００〜８３０に対応するステップであるので、ここでは、これらの詳細についての説明を省略する。
【００７８】
そして、この制御装置は、先に説明した図９に示したルーチンにおいて、ステップ９２０〜９３０にて使用される筒内ガス温度偏差ΔTaに代わりに筒内ガス圧力偏差ΔPa（即ち、上記目標着火時期ATDCθrefにおける目標筒内ガス圧力Parefからステップ１５３０にて算出される上記筒内ガス圧力Paを減じた値）を使用するルーチンを所定時間の経過毎に繰り返し実行することで、ATDCθrefにおける筒内ガス圧力Paが同目標筒内ガス圧力Parefと一致するように燃料噴射開始時期、及び燃料噴射圧力をフィードバック制御する。この結果、混合気の着火時期が目標着火時期ATDCθrefに一致せしめられる。
【００７９】
以上、説明したように、本発明の実施形態に係る内燃機関の筒内ガス圧力の推定方法によれば、上記基準筒内ガス圧力Pabase(CA)、並びに、吸気酸素濃度についての傾きgradPa_ROXC(CA)、吸気管圧力についての傾きgradPa_Pb(CA)、吸気温度についての傾きgradPa_Tb(CA)、圧縮開始時における筒内容積についての傾きgradPa_Vastart(CA)、及び燃焼室壁面温度についての傾きgradPa_Twall(CA)を取得するための各データを予め実験等を通して取得し、それぞれマップとしてＲＯＭ６２に記憶しておくことで、簡易な計算により任意の（所望の）クランク角度CAにおける着火前の筒内ガス圧力を精度良く推定することができた。
【００８０】
本発明は上記実施形態に限定されることはなく、本発明の範囲内において種々の変形例を採用することができる。例えば、上記内燃機関の筒内ガス温度の推定方法の実施形態では、４つの因子の各々と筒内ガス温度との各線形関係を利用して同筒内ガス温度を推定しているが、同４つの因子の一部の各々と筒内ガス温度との各線形関係を利用して同筒内ガス温度を推定するように構成してもよい。また、４つの因子に他の因子を加えた５つ以上の因子の各々と筒内ガス温度との各線形関係を利用して同筒内ガス温度を推定するように構成してもよい。
【００８１】
同様に、上記内燃機関の筒内ガス圧力の推定方法の実施形態では、５つの因子の各々と筒内ガス圧力との各線形関係を利用して同筒内ガス圧力を推定しているが、同５つの因子の一部の各々と筒内ガス圧力との各線形関係を利用して同筒内ガス圧力を推定するように構成してもよい。また、５つの因子に他の因子を加えた６つ以上の因子の各々と筒内ガス圧力との各線形関係を利用して同筒内ガス圧力を推定するように構成してもよい。
【００８２】
また、上記内燃機関の筒内ガス温度の推定方法の実施形態、及び上記内燃機関の筒内ガス圧力の推定方法の実施形態においては、筒内ガスの組成として筒内ガスの酸素濃度を採用しているが、例えば、吸気濃度ＮＯｘセンサにより検出される吸気ＮＯｘ濃度を採用してもよい。
【００８３】
また、上記内燃機関の筒内ガス温度の推定方法の実施形態、及び上記内燃機関の筒内ガス圧力の推定方法の実施形態においては、各因子の値をセンサ等により物理的に検出しているが、所定の計算により推定するように構成してもよい。
【図面の簡単な説明】
【００８４】
【図１】本発明の実施形態に係る内燃機関の筒内ガス温度推定方法を実施する内燃機関の制御装置を４気筒内燃機関（ディーゼル機関）に適用したシステム全体の概略構成図である。
【図２】或る一つの気筒のシリンダ内(筒内)に吸気マニホールドからガスが吸入され、筒内に吸入された筒内ガスが排気マニホールドへ排出される様子を模式的に示した図である。
【図３】任意のクランク角度における吸気酸素濃度と筒内ガス温度との線形関係を示した図である。
【図４】任意のクランク角度における吸気酸素濃度についての筒内ガス温度補正量を求めるための式の理解を容易にするための図である。
【図５】任意のクランク角度における吸気温度についての筒内ガス温度補正量を求めるための式の理解を容易にするための図である。
【図６】任意のクランク角度における圧縮開始時における筒内容積についての筒内ガス温度補正量を求めるための式の理解を容易にするための図である。
【図７】任意のクランク角度における燃焼室壁面温度についての筒内ガス温度補正量を求めるための式の理解を容易にするための図である。
【図８】図１に示したＣＰＵが実行する筒内ガス温度の計算を行うためのルーチンを示したフローチャートである。
【図９】図１に示したＣＰＵが実行する燃料噴射量等を制御するためのルーチンを示したフローチャートである。
【図１０】図１に示したＣＰＵが図９に示したルーチンを実行する際に参照する燃料噴射量を決定するためのテーブルである。
【図１１】図１に示したＣＰＵが図９に示したルーチンを実行する際に参照する基本燃料噴射時期を決定するためのテーブルである。
【図１２】図１に示したＣＰＵが図９に示したルーチンを実行する際に参照する基本燃料噴射圧力を決定するためのテーブルである。
【図１３】図１に示したＣＰＵが図９に示したルーチンを実行する際に参照する噴射時期補正値を決定するためのテーブルである。
【図１４】図１に示したＣＰＵが図９に示したルーチンを実行する際に参照する噴射圧力補正値を決定するためのテーブルである。
【図１５】本発明の実施形態に係る内燃機関の筒内ガス圧力推定方法を実施する内燃機関の制御装置のＣＰＵが実行する筒内ガス圧力の計算を行うためのルーチンを示したフローチャートである。
【符号の説明】
【００８５】
２１…燃料噴射弁、２２…燃料噴射用ポンプ、６０…電気制御装置、６１…ＣＰＵ、７２…吸気温センサ、７３…吸気管圧力センサ、７４…クランクポジションセンサ、７６…吸気酸素濃度センサ、７７…壁面温度センサ
【図１】
2005090368.tif 000003
【図２】
2005090368.tif 000004
【図３】
2005090368.tif 000005
【図４】
2005090368.tif 000006
【図５】
2005090368.tif 000007
【図６】
2005090368.tif 000008
【図７】
2005090368.tif 000009
【図８】
2005090368.tif 000010
【図９】
2005090368.tif 000011
【図１０】
2005090368.tif 000012
【図１１】
2005090368.tif 000013
【図１２】
2005090368.tif 000014
【図１３】
2005090368.tif 000015
【図１４】
2005090368.tif 000016
【図１５】
2005090368.tif 000017