(57)【特許請求の範囲】
【請求項１】
ディーゼルエンジンの回転速度とトルクを変数とする二次元平面上の離散的な複数の点それぞれにおける定常運転時の排気酸素濃度を取得する第１ステップと、
ディーゼルエンジンの回転速度とトルクを変数とする二次元平面上の複数の点を繋ぐ経路に沿うように運転条件を時間的に変化させた過渡運転を行なうようにエンジンの運転条件を制御して、該過渡運転時の、時間的に変化する第１の排気酸素濃度曲線を含むエンジン特性値を測定する第２ステップと、
ディーゼルエンジンの回転速度を、前記経路に沿って、前記過渡運転時の時間的な変化に合わせて繋いだ回転速度曲線、並びに、エンジン特性値から算出される排気酸素濃度およびＥＧＲ率のそれぞれを前記過渡運転時に測定されたエンジン特性値の時間的な変化に従って時間的に変化させた第２の排気酸素濃度曲線およびＥＧＲ率曲線のうちの、該第２の排気酸素濃度曲線を、回転速度およびＥＧＲ率をパラメータとして含む複数のパラメータで記述される時定数を持つ１次遅れ系に入力し、該時定数のパラメータのうちの回転速度およびＥＧＲ率に前記回転速度曲線および前記ＥＧＲ率曲線をそれぞれ代入するとともに、該１次遅れ系の出力として、前記第１の排気酸素濃度曲線からの誤差が最小となる第３の排気酸素濃度曲線が得られるように、該複数のパラメータのうちの回転速度およびＥＧＲ率を除くパラメータの値を算出する第３ステップと、
エンジン制御パラメータからなる説明変数とエンジン特性値からなる目的変数との関係を、前記第３ステップで算出した、時間的に変化する時定数を持つ１次遅れの排気酸素濃度を含んで記述したエンジンモデルを作成する第４ステップとを有することを特徴とするエンジンモデル作成方法。
【請求項２】
前記第２ステップが、前記二次元平面上の、回転速度一定の複数の点を繋ぐ経路に沿うように運転条件を時間的に変化させた過渡運転を行なうようにエンジンの運転条件を制御して、該過渡運転時の、時間的に変化する第１の排気酸素濃度曲線を含むエンジン特性値を測定するステップであることを特徴とする請求項１に記載のエンジンモデル作成方法。
【請求項３】
前記時定数が、
【数１】
0006579975.tif 000034
で表される、切片と係数と回転速度とＥＧＲ率とをパラメータとする時定数であることを特徴とする請求項１または２に記載のエンジンモデル作成方法。
【請求項４】
前記時定数が、
【数２】
0006579975.tif 000035
で表される、切片と係数と回転速度とＥＧＲ率とをパラメータとする時定数であることを特徴とする請求項１または２に記載のエンジンモデル作成方法。
【請求項５】
前記時定数が、
【数３】
0006579975.tif 000036
で表される、切片と係数と回転速度とＥＧＲ率とをパラメータとする時定数であることを特徴とする請求項１または２に記載のエンジンモデル作成方法。
【請求項６】
前記第２ステップが、エンジンモデル作成対象のエンジンにダイナモを連結した実験システム上で、前記第１の排気酸素濃度曲線の測定を行なうステップであることを特徴とする請求項１から５のうちのいずれか１項に記載のエンジンモデル作成方法。
【請求項７】
エンジン制御パラメータの調整により、エンジンモデル作成対象のディーゼルエンジンの運転を制御する、該ディーゼルエンジンに連結されるダイナモを含むエンジン制御装置、
前記ディーゼルエンジン運転時の、該ディーゼルエンジンの、排気酸素濃度およびＥＧＲ率を含む特性を測定するエンジン特性測定装置、および、
ディーゼルエンジンの回転速度とトルクを変数とする二次元平面上の離散的な複数の点それぞれにおける定常運転時の排気酸素濃度を取得する第１ステップと、
前記エンジン制御装置に、ディーゼルエンジンの回転速度とトルクを変数とする二次元平面上の複数の点を繋ぐ経路に沿うように運転条件を時間的に変化させた過渡運転を行なうようにエンジンの制御パラメータを調整させて、前記エンジン特性測定装置に、該過渡運転時の、時間的に変化する第１の排気酸素濃度曲線を含むエンジン特性値を測定させる第２ステップと、
ディーゼルエンジンの回転速度を、前記経路に沿って、前記過渡運転時の時間的な変化に合わせて繋いだ回転速度曲線、並びに、エンジン特性値から算出される排気酸素濃度およびＥＧＲ率のそれぞれを前記過渡運転時に測定されたエンジン特性値の時間的な変化に従って時間的に変化させた第２の排気酸素濃度曲線およびＥＧＲ率曲線のうちの、該第２の排気酸素濃度曲線を、回転速度およびＥＧＲ率をパラメータとして含む複数のパラメータで記述される時定数を持つ１次遅れ系に入力し、該時定数のパラメータのうちの回転速度およびＥＧＲ率に前記回転速度曲線および前記ＥＧＲ率曲線をそれぞれ代入するとともに、該１次遅れ系の出力として、前記第１の排気酸素濃度曲線からの誤差が最小となる第３の排気酸素濃度曲線が得られるように、該複数のパラメータのうちの回転速度およびＥＧＲ率を除くパラメータの値を算出する第３ステップと、
エンジン制御パラメータからなる説明変数とエンジン特性値からなる目的変数との関係を、前記第３ステップで算出した、時間的に変化する時定数を持つ１次遅れの排気酸素濃度を含んで記述したエンジンモデルを作成する第４ステップとを有するエンジンモデル作成方法を実行するエンジンモデル作成装置とを備えたことを特徴とするエンジンモデル作成システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、エンジンの様々な運転条件と排気ガス中のＮＯｘや燃費などとの関係を表わすエンジンモデルを作成するエンジンモデル作成方法およびエンジンモデル作成システムに関する。
【背景技術】
【０００２】
近年の自動車の普及と環境科学の進歩に伴い、燃費や排ガスの規制が厳しく強化されてきている。このため、精密なエンジンモデルを作成し、エンジンの様々な運転条件下での排気ガス中のＮＯｘや燃費などを高精度にシミュレーションすることが求められている。
【０００３】
ここで、エンジンモデルの作成にあたっては、エンジンにダイナモを連結した実験システム上で実験を行ない、その実験結果に基づいてエンジンモデルを作成することが従来より行なわれている（特許文献１参照）。
【０００４】
エンジンモデルを作成して、自動車に搭載されたエンジンの制御パラメータをその作成されたエンジンパラメータに基づいて調整することで、エンジンの回転速度やトルクを所要の値に保ちながら排気ガス中のＮＯｘの低減や燃費の向上を図る努力がなされている。
【０００５】
従来は、エンジンの回転速度とトルクとからなる２次元平面上の各点での定常運転によるエンジンモデルが作成され、活用されている。しかしながら、この定常運転によるエンジンモデルの場合、エンジンの回転速度やトルクが変化しつつある過渡運転時、例えば運転開始時などについては、十分には最適化されたエンジンモデルになっていない。すなわち、エンジンモデルから予測されるＮＯｘの値や燃費等のエンジン特性値が、エンジンを実際に自動車に搭載して運転して実測したときのエンジン特性値に高精度に近似しているとは言い難く、それらの間の誤差が十分に小さいとは言い難い。
【０００６】
この誤差が生じる大きな原因の１つは、過渡運転時の排気酸素濃度、およびその排気酸素濃度から算出される吸気酸素濃度がエンジンモデルに正しく組み込まれていないことにある。すなわち、定常運転によるエンジンモデルの場合、排気酸素濃度は、エンジンの現在の回転速度やトルクによって定まるのに対し、過渡運転の場合、排気酸素濃度は、エンジンの現在の回転速度やトルクだけでは定まらず、現在に至るまでの回転速度やトルクの変化の影響が時間的に遅れて現在の排気酸素濃度に反映される。
【０００７】
エンジンの吸気酸素濃度や排気酸素濃度はＮＯｘやＰＭ等の排出に強く関連している。例えば吸気酸素濃度が高くなると火炎温度が上がり、ＮＯｘが増大することが知られている。したがって、正確なシミュレーションを行なうためには、過渡運転時を含む様々な運転条件下での排気酸素濃度や吸気酸素濃度をエンジンモデルに高精度に組み入れることが必要である。
【０００８】
ここで、特許文献２には、吸気酸素濃度に応じて燃料噴射時期を切り替えることでＮＯｘやＰＭの抑制を実現する技術が提案されている。
【０００９】
また、特許文献３には、触媒の状態も考慮して噴射形態を吸気酸素濃度の変化に応じて制御することでＮＯｘやＰＭの抑制を実現する技術が提案されている。
【００１０】
しかしながら、これら引用文献２，３に提案された技術は、吸気酸素濃度や排気酸素濃度を高精度に測定あるいは算出した後の制御方法に関する技術であって、吸気酸素濃度や排気酸素濃度の高精度の算出方法を示したものではない。
【００１１】
なお、後述の実施形態における説明のために、非特許文献１を挙げておく。
【先行技術文献】
【特許文献】
【００１２】
【特許文献１】特開２００９−２１０４２６号公報
【特許文献２】特開２００９−１２１３８９号公報
【特許文献３】特開２０１０−２３６４６０号公報
【非特許文献】
【００１３】
【非特許文献１】ｈｔｔｐ：／／ｉｓｍｒｅｐｏ．ｉｓｍ．ａｃ．ｊｐ／ｄｓｐａｃｅ／ｂｉｔｓｔｒｅａｍ／１０７８７／８３６／１／ｏｐｅｎｈｏｕｓｅ２０１１−ａ１５ｓｏｍｅｙａ.ｐｄｆ
【発明の概要】
【発明が解決しようとする課題】
【００１４】
本発明は、上記事情に鑑み、過渡運転時の排気酸素濃度および吸気酸素濃度が高精度に推測されたエンジンモデルを作成するエンジンモデル作成方法およびエンジンモデル作成システムを提供することを目的とする。
【課題を解決するための手段】
【００１５】
上記目的を達成する本発明のエンジンモデル作成方法は、
ディーゼルエンジンの回転速度とトルクとを変数とする二次元平面上の離散的な複数の点それぞれにおける定常運転時の排気酸素濃度を取得する第１ステップと、
ディーゼルエンジンの回転速度とトルクを変数とする二次元平面上の複数の点を繋ぐ経路に沿うように運転条件を時間的に変化させた過渡運転を行なうようにエンジンの運転条件を制御して、その過渡運転時の、時間的に変化する第１の排気酸素濃度曲線を含むエンジン特性値を測定する第２ステップと、
ディーゼルエンジンの回転速度を、上記経路に沿って、過渡運転時の時間的な変化に合わせて繋いだ回転速度曲線、並びに、エンジン特性値から算出される排気酸素濃度およびＥＧＲ率のそれぞれを過渡運転時に測定されたエンジン特性値の時間的な変化に従って時間的に変化させた第２の排気酸素濃度曲線およびＥＧＲ率曲線のうちの、第２の排気酸素濃度曲線を、回転速度およびＥＧＲ率をパラメータとして含む複数のパラメータで記述される時定数を持つ１次遅れ系に入力し、その時定数のパラメータのうちの回転速度およびＥＧＲ率に上記の回転速度曲線およびＥＧＲ率曲線をそれぞれ代入するとともに、その１次遅れ系の出力として、第１の排気酸素濃度曲線からの誤差が最小となる第３の排気酸素濃度曲線が得られるように、それら複数のパラメータのうちの回転速度およびＥＧＲ率を除くパラメータの値を算出する第３ステップと、
エンジン制御パラメータからなる説明変数とエンジン特性値からなる目的変数との関係を、第３ステップで算出した、時間的に変化する時定数を持つ１次遅れの排気酸素濃度を含んで記述したエンジンモデルを作成する第４ステップとを有することを特徴とする。
【００１６】
本発明のエンジンモデル作成方法は、定常運転時の排気酸素濃度から作成した第２の排気酸素濃度曲線が過渡運転時の実測値からなる第１の排気酸素濃度曲線に一致するように、１次遅れ系の時定数を定めたものである。しかもその時定数は、回転速度およびＥＧＲ率に応じて変化する時定数である。
【００１７】
本発明のエンジンモデル作成方法では、このように、回転速度およびＥＧＲ率に応じて変化する時定数を持つ１次遅れ系の出力として、実測値からなる第１の排気酸素濃度曲線からの誤差の小さい第３の排気酸素濃度曲線を算出することで、過渡運転時の排気酸素濃度を高精度に組み入れたエンジンモデルが作成される。
【００１８】
なお、本発明のエンジンモデル作成方法によれば、排気酸素濃度が高精度に求められるため、この排気酸素濃度から算出される吸気酸素濃度についても高精度に求められることになる。
【００１９】
ここで、本発明のエンジンモデル作成方法において、上記第２ステップが、上記二次元平面上の、回転速度一定の複数の点を繋ぐ経路に沿うように運転条件を時間的に変化させた過渡運転を行なうようにエンジンの運転条件を制御して、その過渡運転時の、時間的に変化する第１の排気酸素濃度曲線を含むエンジン特性値を測定するステップであることが好ましい。
【００２０】
回転速度一定の経路を設定すると、第３の排気酸素濃度曲線を算出する際の演算処理が容易となる。
【００２１】
また、本発明のエンジンモデル作成方法において、上記時定数が、
【００２２】
【数１】
0006579975.tif 000002
で表される、切片と係数と回転速度とＥＧＲ率とをパラメータとする時定数であることが好ましい。
【００２３】
あるいは、本発明のエンジンモデル作成方法において、上記時定数が、
【００２４】
【数２】
0006579975.tif 000003
で表される、切片と係数と回転速度とＥＧＲ率とをパラメータとする時定数であることも好ましい態様である。
【００２５】
さらには、本発明のエンジンモデル作成方法において、上記時定数が、
【００２６】
【数３】
0006579975.tif 000004
で表される、切片と係数と回転速度とＥＧＲ率とをパラメータとする時定数であることも好ましい態様である。
【００２７】
後述するように、これら式のいずれによっても、第２の排気酸素濃度曲線に高精度に一致した第３の排気酸素濃度曲線が算出される。
【００２８】
また、本発明のエンジンモデル作成方法において、上記第２ステップが、エンジンモデル作成対象のエンジンにダイナモを連結した実験システム上で、排気酸素濃度曲線の測定を行なうステップであることが好ましい。
【００２９】
自動車に搭載された状態のエンジンについて上記の調整および測定を行なうのは極めて大変な作業であり、実験システム上で行なうことが好ましい。
【００３０】
また、上記目的を達成する本発明のエンジンモデル作成システムは、
エンジン制御パラメータの調整により、エンジンモデル作成対象のディーゼルエンジンの運転を制御する、ディーゼルエンジンに連結されるダイナモを含むエンジン制御装置、
ディーゼルエンジン運転時の、ディーゼルエンジンの、排気酸素濃度およびＥＧＲ率を含む特性を測定するエンジン特性測定装置、および、
ディーゼルエンジンの回転速度とトルクとを変数とする二次元平面上の離散的な複数の点それぞれにおける定常運転時の排気酸素濃度を取得する第１ステップと、
上記エンジン制御装置に、ディーゼルエンジンの回転速度とトルクを変数とする二次元平面上の複数の点を繋ぐ経路に沿うように運転条件を時間的に変化させた過渡運転を行なうようにエンジンの制御パラメータを調整させて、上記エンジン特性測定装置に、過渡運転時の、時間的に変化する第１の排気酸素濃度曲線を含むエンジン特性値を測定させる第２ステップと、
ディーゼルエンジンの回転速度を、上記経路に沿って、過渡運転時の時間的な変化に合わせて繋いだ回転速度曲線、並びに、エンジン特性値から算出される排気酸素濃度およびＥＧＲ率のそれぞれを過渡運転時に測定されたエンジン特性値の時間的な変化に従って時間的に変化させた第２の排気酸素濃度曲線およびＥＧＲ率曲線のうちの、第２の排気酸素濃度曲線を、回転速度およびＥＧＲ率をパラメータとして含む複数のパラメータで記述される時定数を持つ１次遅れ系に入力し、その時定数のパラメータのうちの回転速度およびＥＧＲ率に上記の回転速度曲線およびＥＧＲ率曲線をそれぞれ代入するとともに、その１次遅れ系の出力として、第１の排気酸素濃度曲線からの誤差が最小となる第３の排気酸素濃度曲線が得られるように、それら複数のパラメータのうちの回転速度およびＥＧＲ率を除くパラメータの値を算出する第３ステップと、
エンジン制御パラメータからなる説明変数とエンジン特性値からなる目的変数との関係を、前記第３ステップで算出した、時間的に変化する時定数を持つ１次遅れの排気酸素濃度を含んで記述したエンジンモデルを作成する第４ステップとを有するエンジンモデル作成方法を実行するエンジンモデル作成装置とを備えたことを特徴とする。
【発明の効果】
【００３１】
上記の本発明によれば、過渡運転時の排気酸素濃度や吸気酸素濃度が高精度に推測されたエンジンモデルが作成され、排気ガス中のＮＯｘの量や燃費等のエンジン特性値が高精度に推測される。
【図面の簡単な説明】
【００３２】
【図１】本発明の一実施形態としてのエンジンモデル作成システムの構成図である    。
【図２】ディーゼルエンジンにおける吸気および排気の流れを示した模式図であ    る。
【図３】第１ステップにおける実験のイメージを示した図である。
【図４】第２ステップにおける運転条件の変化のさせ方を示した図である。
【図５】比較例としての排気酸素濃度曲線の一例を示した図である。
【図６】「切片」と「係数」を変数とする二次元空間を示した模式図である。
【図７】排気酸素濃度曲線を示した図である。
【図８】排気に含まれるＮＯｘの変化曲線を示した図である。
【発明を実施するための形態】
【００３３】
以下、本発明の実施の形態を説明する。
【００３４】
図１は、本発明の一実施形態としてのエンジンモデル作成システムの構成図である。
【００３５】
ここに示すエンジンモデル作成システム１は、エンジンモデル作成対象のエンジン１０に連結されたダイナモ２０を備えている。また、このエンジンモデル作成システム１は、制御・計測用パソコン３０、スロットルアクチュエータ４０、ＥＣＵ（エンジンコントロールユニット）５０、排ガス分析計６０、およびダイナモ制御盤７０を備えている。また、ここには、エンジン１０の排気管１１内を流れる排ガスの排気温度を測定する温度センサ８１や、排ガスの流量を測定する流量センサ８２や、エンジン１０の回転軸１２のトルクを測定するトルクセンサ８３や、回転軸の回転速度を測定する回転速度センサ８４などの様々なセンサが備えられている。
【００３６】
制御・計測用パソコン３０は、例えば、スロットルアクチュエータ４０を介してのエンジン１０のスロットル（不図示）の開度の制御によるトルク制御、ダイナモ制御盤７０を介してのダイナモ２０の制御によるエンジン１０の回転軸１２の回転速度制御、ＥＣＵ５０を介しての噴射進角の制御等を行なう。また、この制御・計測用パソコン３０では、エンジン１０の回連速度、トルク、排気温度、排気流量等のデータの収集が行なわれる。排ガス分析計６０では、排気管１１内を流れる排ガスの成分分析が行なわれ、制御・計測用パソコン３０には、その分析結果も収集される。さらに、この制御・計測用パソコン３０には、ＥＣＵ５０による、例えばＥＧＲ率などのエンジン１０の制御パラメータの値も収集される。さらに、この制御・計測用パソコン３０は、燃料の供給、排ガスの処理等を行なう様々な環境設備９０の制御も行なう。
【００３７】
この図１に示すエンジンモデル作成システム１において、ダイナモ２０、スロットルアクチュエータ４０、ＥＣＵ５０、およびダイナモ制御盤７０や、制御・計測用パソコン３０のそれらを制御する機能の組合せが、本発明にいうエンジン制御装置の一例に相当する。
【００３８】
また、この図１に示すエンジンモデル作成システム１のうちの、エンジン１０の回転速度やトルクなどを測定する各種センサと、排ガス分析計６０と、制御・計測用パソコン３０のデータ収集機能との組合せが、本発明にいうエンジン特性測定装置の一例に相当する。
【００３９】
また、この図１に示すエンジンモデル作成システム１の制御・計測用パソコン３０は、収集した各種データに基づいてエンジンモデルが作成する機能を有している。制御・計測用パソコン３０の、エンジンモデル作成機能は、本発明にいうエンジンモデル作成装置の一例に相当する。
【００４０】
図２は、ディーゼルエンジンにおける吸気および排気の流れを示した模式図である。
【００４１】
大気がコンプレッサを経て吸気として取り込まれる。ただし、この吸気は、ディーゼルエンジンに取り込まれる前に、そのディーゼルエンジンからの排気がＥＧＲ率分だけ吸気に混入し、ディーゼルエンジンには、排気が混入した空気が取り込まれる。ディーゼルエンジンに取り込まれた吸気中の酸素は、燃料噴射量に応じて燃焼に使用されて、ディーゼルエンジンから排気される。この排気はＥＧＲ率分だけ再び吸気側に回り込み、残りの排気はコンプレッサ作動用のタービンを動かして大気中に排気される。
【００４２】
したがって、ディーゼルエンジンへの吸気の組成やディーゼルエンジンからの排気の組成は、排気が吸気側に回り込むため、そのディーゼルエンジンのこれまでの運転の履歴が反映されたものとなり、過渡運転と定常運転とでは振る舞いが異なる。
【００４３】
そこで、本実施形態における、図１に示すエンジンモデル作成システム１では、以下に説明する第１ステップから第４ステップにより、過渡運転を考慮した高精度なエンジンモデルを作成する。
【００４４】
第１ステップでは、エンジンの回転速度とトルクとを変数とする二次元平面上の離散的な複数の点それぞれにおける定常運転時の排気酸素濃度を取得する。この第１ステップでは、図１に示したエンジンモデル作成システム１を使ってエンジン１０を実際に動かして、排気酸素濃度を、排ガス分析計６０から取得してもよく、あるいは、定常運転におけるデータを基にして作成したエンジンモデルが既に存在するときは、そのエンジンモデルを用いたシミュレーションにより排気酸素濃度を算出してもよい。ここでは、図１に示すエンジンモデル作成システム１を使って実測する場面について説明する。
【００４５】
図３は、第１ステップにおける実験のイメージを示した図である。
【００４６】
ここには、エンジンの回転速度とトルクとを変数とする二次元平面が示されている。この二次元平面上に示された○印は、その１つ１つが１つの測定点を表わしている。すなわち、ここでは、回転速度とトルクが１つの○印に対応する各値に固定された状態となるようにエンジンを動かし、初期の計測値は無視して安定した後半部の計測値の平均値を求める。これを、○印の各測定点１つ１つについて実行する。これにより、エンジンの回転速度とトルクとからなる二次元平面上の複数の点それぞれの定常運転（エンジンの回転速度やトルクが変化せずに一定のままの運転）における排気酸素濃度が求められる。
【００４７】
第２ステップでは、エンジンの回転速度とトルクを変数とする二次元平面上の複数の点を繋ぐ経路に沿うように運転条件を時間的に変化させた過渡運転を行なうようにエンジンの運転条件を制御して、その過渡運転時の、時間的に変化する排気酸素濃度曲線を含むエンジン特性値を測定する。ここでの測定により得られる排気酸素濃度曲線は、本発明にいう第１の排気酸素濃度曲線の一例に相当する。
【００４８】
図４は、第２ステップにおける運転条件の変化のさせ方を示した図である。
【００４９】
図４（Ａ）は、二次元平面上の経路を示している。ここでは、一例として、エンジン回転速度は各経路ごとに一定とし、各経路内ではトルクが上昇して下降するように経路を設定している。各経路ごとのエンジンの回転速度は、例えば、１０００ｒ／ｍｉｎ、２０００ｒ／ｍｉｎ、３０００ｒ／ｍｉｎである。
【００５０】
また、図４（Ｂ）は、トルクの変化曲線を示している。ここでは、ある低いトルクから出発し、途中で直線的にトルクを上昇させ、一定時間高いトルクを保ち、その後直線的にトルクを下降させ、一定時間低いトルクを保つように、トルクを変化させる。
【００５１】
ここでは、このような複数な経路それぞれに沿うように運転条件を制御することで過渡運転を行ない、その過渡運転時の排気酸素濃度曲線を含むエンジン特性値を測定する。
【００５２】
次に、第３ステップについて説明する。
【００５３】
この第３ステップでは、回転速度曲線と、排気酸素濃度曲線と、ＥＧＲ率曲線が用いられる。これらのうち、回転速度曲線は、ディーゼルエンジンの回転速度を、上記の第２ステップで採用した経路に沿って、第２ステップにおける過渡運転時の時間的な変化に合わせて繋いだ曲線である。また、排気酸素濃度およびＥＧＲ率は、いずれも、エンジン特性値から算出される値であり、この第３ステップでは、それら排気酸素濃度およびＥＧＲ率のそれぞれを、第２ステップにおける過渡運転時に測定されたエンジン特性値の時間的な変化に従って変化させた排気酸素濃度曲線およびＥＧＲ率曲線が用いられる。ここで算出される排気酸素濃度曲線は、第２ステップでの測定により得られた第１の排気酸素濃度曲線と区別するために、ここでは、第２の排気酸素濃度曲線と称する。
【００５４】
ここで、排気酸素濃度を算出するにあたっては、下記の式（１）、（２）に従って排気酸素流量Weo2および排気流量Ｗeを計算する。
【００５５】
【数４】
0006579975.tif 000005
【００５６】
【数５】
0006579975.tif 000006
そして、排気酸素濃度Ｃgwo2は、それらの式（１）、（２）に従って算出された排気酸素流量Weo2および排気流量Ｗeを使って、下記の式（３）に従って算出される。
【００５７】
【数６】
0006579975.tif 000007
ただし、
【００５８】
【数７】
0006579975.tif 000008
である。
【００５９】
また、ＥＧＲ率は、下記の式（４）によって算出することができる。
【００６０】
【数８】
0006579975.tif 000009
あるいは、ＥＧＲ率を、下記の式（５）に従って算出してもよい。
【００６１】
【数９】
0006579975.tif 000010
によって算出してもよい。
【００６２】
ちなみに、吸気酸素濃度は、ＥＧＲ率と排気酸素濃度を用いて、以下の式（６）、（７）に従って算出される。
【００６３】
ここでは先ず、ＥＧＲ流量ＷEGRが、
【００６４】
【数１０】
0006579975.tif 000011
に従って計算され、その計算されたＥＧＲ流量ＷEGRを用いて、吸気酸素濃度Ｃiwo2が、
【００６５】
【数１１】
0006579975.tif 000012
に従って計算される。
【００６６】
この第３ステップでは、さらに、上記の計算により求めた第２の排気酸素濃度曲線を１次遅れ系に入力し、その１次遅れ系の出力としての排気酸素濃度曲線を生成する。この１次遅れ系は、回転速度とＥＧＲ率をパラメータとして含む複数のパラメータで記述される時定数を持つ１次遅れ系である。この１次遅れ系の具体例については、後述する。ここでは、時定数のパラメータである回転速度およびＥＧＲ率に回転速度曲線およびＥＧＲ率曲線を代入するとともに、この１次遅れ系の出力としての排気酸素濃度曲線が第２ステップでの測定により得られた排気酸素濃度曲線（第１の排気酸素濃度曲線）からの誤差が最小となるように回転速度およびＥＧＲ率以外のパラメータの値を定めた時定数が採用される。
【００６７】
すなわち、ここで採用される時定数は、先ずは、上記のようにして求められたＥＧＲ率を繋げることで生成されたＥＧＲ率曲線に応じて時間的に変化する時定数である。また、この時定数を決定づけるパラメータの１つとして回転速度が含まれている。したがって、過渡運転の経路として回転速度が変化するような経路を採用したときは、この時定数は、その経路に沿って回転速度を繋げることで生成された回転速度曲線に応じても時間的に変化することになる。ただし、ここでは、図４（Ａ）に示すように、回転速度は各経路ごとに一定であるため、ここでは回転速度は各経路ごとの定数となる。
【００６８】
また、この時定数を決定づける回転速度およびＥＧＲ率を除く他のパラメータについては、時間的には変化しない固定値が採用される。ただし、その固定値は、この１次遅れ系の出力としての排気酸素濃度曲線が第２ステップでの測定により得られた排気酸素濃度曲線（第１の排気酸素濃度曲線）からの誤差が最小となるように決定した固定値である。このようにして、この１次遅れ系からは、上記の回転速度やＥＧＲ率に応じて変化しつつ、それらを除くパラメータについては、上記の誤差が最小となるように定めた時定数を用いて、定常状態の排気酸素濃度を繋ぐことにより生成された排気酸素濃度曲線（第２の排気酸素濃度曲線）が１次遅れを生じた排気酸素濃度曲線が生成される。上記の誤差が最小となるように定めた時定数を用いたときの、この１次遅れ系の出力としての排気酸素濃度曲線は、本発明にいう第３の排気酸素濃度曲線の一例に相当する。
【００６９】
次の第４ステップでは、エンジン制御パラメータからなる説明変数と、排気ガス成分および成分量や燃費等を含むエンジン特性値からなる目的変数との関係を、上記の第３ステップで算出された、時間的に変化する時定数を持つ１次遅れの排気酸素濃度を含んで記述したエンジンモデルが作成される。この第４ステップで作成されるエンジンモデルは、エンジンの過渡運転時の排気酸素濃度や、この排気酸素濃度から算出される吸気酸素濃度が高精度に再現された、高精度のエンジンモデルであり、実測値に一層近づいたエンジン特性値が算出されるモデルとなる。
【００７０】
次に、第３ステップの説明に戻り、時定数の求め方の具体例について説明する。
【００７１】
図５は、比較例としての排気酸素濃度曲線の一例を示した図である。
【００７２】
ここで、図５（Ａ），（Ｂ）のいずれにおいても、横軸は時間［ｓ］、縦軸は排気酸素濃度である。
【００７３】
また、図５（Ａ）はＥＧＲ全閉（ＥＧＲ率ゼロ％）時の排気酸素濃度を表わしていおり、図５（Ｂ）はＥＧＲ使用時の排気酸素濃度を表わしている。ここで、ＥＧＲ使用時のＥＧＲ率は、制御・計測用コンピュータ３０（図１参照）では制御されずにＥＣＵ５０（図１参照）によって制御されており、制御・計測用コンピュータ３０は、ＥＣＵ５０からＥＧＲ率の情報を受取る仕組みとなっている。このＥＧＲ率は時間的に変化する値である。
【００７４】
また、図５（Ａ），（Ｂ）のいずれにおいても、実線は実測値を表わしている。一方、破線は、図３を参照して説明したような定常運転で得られたデータに基づいて算出した計算値である。実線で示す実測値と破線で示す計算値とではオフセットがかかっているため、この図５における縦軸方向のずれは意味がなく、ここでは横軸（時間軸）方向へのずれに着目する。
【００７５】
ＥＧＲ全閉時の場合、図２に示す「ＥＧＲ」が閉じているため、ディーゼルエンジンからの排気が吸気側に回り込むことはない。すなわち、ＥＧＲ全閉時の場合、大気から取り込んだ空気がそのままディーゼルエンジンへの吸気となり、ディーゼルエンジンからの排気がそのまま大気中に排気される。
【００７６】
このため、ＥＧＲ全閉時は、図５（Ａ）に示すように、実測値に対し時間軸方向に僅かにずれただけのほぼ一致した波形の計算値が得られている。
【００７７】
これに対し、ＥＧＲ使用時は、ディーゼルエンジンからの排気がその時点のＥＧＲ率分だけ吸気側に回り込む。すなわち、一部分が排気と吸気との間を循環する。このため、ＥＧＲ使用時は、図５（Ｂ）に示すように、計測値は、実測値と比べＥＧＲ全閉時（図５（Ａ））よりも大きく、時間的にずれている。また、波形の形もＥＧＲ全閉時と比べ実測値とは異なっている。ＥＧＲ使用時の波形の時間的なずれは、後述する通り、ＥＧＲ率とエンジンの回転速度とに依存している。
【００７８】
このように、定常運転時のデータから作成したエンジンモデルの場合、ＥＧＲ率が時間的に変化するような過渡運転時における排気酸素濃度は高精度には再現されない。本実施形態では、具体的には以下のようにして、実測値からの誤差を低減させた排気酸素濃度曲線（第３の排気酸素濃度曲線）を生成している。
【００７９】
ここで、ＥＧＲ率と排気酸素濃度とを関連づける式を導出する。
【００８０】
新気酸素濃度をａ、
ディーゼルエンジンでの燃焼によって消費される酸素濃度をｂ、
ＥＧＲ率をｒ、
とする。ディーゼルエンジンが停止した状態から運転を開始したと考え、ｎ回目の燃焼前の酸素濃度をｘn、ｎ回目の燃焼後の酸素濃度をｙnと定義し、
【００８１】
【数１２】
0006579975.tif 000013
帰納法により、
【００８２】
【数１３】
0006579975.tif 000014
下記の式（８）を利用して、定常運転における燃焼後の酸素濃度を
ｙ∞
とすると、
【００８３】
【数１４】
0006579975.tif 000015
となり、燃焼後の酸素濃度の漸化式（９）が得られる。
【００８４】
【数１５】
0006579975.tif 000016
入力ｕ，出力ｙ、時定数Ｔの伝達関数は下式で表わされる。
【００８５】
【数１６】
0006579975.tif 000017
ここで、ｓはラプラス変換の微分演算子である。
時間領域の表現に直すと、
【００８６】
【数１７】
0006579975.tif 000018
標本化周期をΔＴとして、後退差分で近似計算する場合、
【００８７】
【数１８】
0006579975.tif 000019
ｙnをステップ入力
【００８８】
【数１９】
0006579975.tif 000020
に対する１次遅れ応答ｙの標本としてとらえ、式（９）と比較すると、
【００８９】
【数２０】
0006579975.tif 000021
同様に前進差分の場合、
【００９０】
【数２１】
0006579975.tif 000022
中心差分の場合、
【００９１】
【数２２】
0006579975.tif 000023
漸化式（９）を拡張し、
【００９２】
【数２３】
0006579975.tif 000024
式（１４）を式（１３）のｙn+1に代入してから、漸化式と比較する。
【００９３】
【数２４】
0006579975.tif 000025
ここで、標本化周期ΔＴはエンジンの回転速度に反比例するので、
【００９４】
【数２５】
0006579975.tif 000026
となる。ただし、αは定数である。
【００９５】
また、ＲＧＲ率ｒを％に変換して
【００９６】
【数２６】
0006579975.tif 000027
とする。
すると、上記の後退差分で導いた式（１１）は、
【００９７】
【数２７】
0006579975.tif 000028
となる。
【００９８】
また、上記の前進差分で導いた式（１２）は、
【００９９】
【数２８】
0006579975.tif 000029
となる。
【０１００】
さらに、上記の中心差分で導いた式（１５）は、
【０１０１】
【数２９】
0006579975.tif 000030
となる。
【０１０２】
さらに、回転速度によらないセンサ等の時定数も考慮して、
【０１０３】
【数３０】
0006579975.tif 000031
とすると、上記の式（１６），（１７），（１８）はそれぞれ、
【０１０４】
【数３１】
0006579975.tif 000032
となる。
【０１０５】
すなわち、これら式（１９），（２０），（２１）のいずれの場合も、「時定数」は、「切片」と「係数」と「回転速度」と「ＥＧＲ率」との４つのパラメータで決定されることになる。
【０１０６】
時定数を算出するにあたり、これら式（１９），（２０），（２１）のいずれの式を用いてもよいが、例えば、一例として式（２０）を用いる。
【０１０７】
また、「回転速度」は、図４に示す各経路ごとに一定値であり、ここでは、その一定値を代入する。ただし、回転速度が時間的に変化する経路を辿る過渡運転を行なったときは、その過渡運転の経路に沿って時間的に変化する回転速度が代入される。
【０１０８】
さらに、「ＥＧＲ率」および「回転速度」以外のパラメータである「切片」と「係数」は、以下のようにして決定される。
【０１０９】
図６は、「切片」と「係数」を変数とする二次元空間を示した模式図である。
【０１１０】
「Ｅ」の文字が記入された○印で示す座標点Ｅ（係数＝９．２５、切片＝０）は、第２ステップでの測定で得られた排気酸素濃度曲線（第１の排気酸素濃度曲線）からの誤差が最小となる排気酸素濃度曲線が得られる「切片」と「係数」からなる座標点Ｅである。ただし、この座標点Ｅの「切片」と「係数」は今から算出するのであって、算出前においては、この座標点Ｅの位置は不明である。
【０１１１】
ここでは、この座標点Ｅの位置は未だ不明ではあるが、この座標点Ｅが含まれていることが確実な、４つの座標点Ａ，Ｂ，Ｃ，Ｄで囲まれた領域を設定する。ここに示した例では、座標点Ａは、「切片」＝０、「係数」＝０、座標点Ｂは、「切片」＝１、「係数」＝０、座標点Ｃは、「切片」＝１、「係数」＝２０、座標点Ｄは、「切片」＝０、「係数」＝２０である。
【０１１２】
図７は、排気酸素濃度曲線を示した図である。図７（Ａ）〜（Ｅ）は、図６に示す各座標点Ａ〜Ｅのそれぞれに対応している。また、図７（Ａ）〜（Ｅ）のそれぞれにおいて、横軸は時間［ｓ］、縦軸は排気酸素濃度である。さらに、図７（Ａ）〜（Ｅ）のそれぞれにおいて、上段側の実線グラフと破線のグラフは、ＥＧＲ全閉時の排気酸素濃度、下段側の実線グラフと破線のグラフはＥＧＲ使用時の排気酸素濃度である。さらに、これら図７（Ａ）〜（Ｅ）における実線のグラフは、上述の第２ステップでの測定で得られた実測値としての排気酸素濃度曲線（第１の排気酸素濃度曲線）である。この実線で示す排気酸素濃度曲線は、図４（Ａ）に示すように、エンジンの回転速度をある一定値に固定し、トルクを図４（Ｂ）に示すように時間的に変化させたときの排気酸素濃度である。この実線で示す排気酸素濃度曲線は、図７（Ａ）〜（Ｅ）のいずれにおいても同一である。一方、これら図７（Ａ）〜（Ｅ）における破線のグラフは、式（２０）で示される時定数を持つ１次遅れ系の出力としての排気酸素濃度曲線である。式（（２０）で示す時定数には、「切片」、「係数」、「回転速度」および「ＥＧＲ率」という４つのパラメータが含まれている。これら４つのパラメータのうちの「ＥＧＲ率」は、時間的には変化するものの、図７（Ａ）〜（Ｅ）のいずれにおいても共通の値を持つパラメータである。また、「回転速度」は経路によって固定値あるいは時間的に変化する値をとるが、この「回転速度」も図７（Ａ）〜（Ｅ）のいずれにおいても共通の値を持つパラメータである。残りの２つのパラメータである「切片」と「係数」は、図６に示す２次元平面上の座標点によって異なっている。
【０１１３】
ここで、図７（Ａ）は、図６の座標点Ａに対応するグラフである。すなわち、図７（Ａ））に示す破線のグラフは、「切片」＝０、かつ「係数」＝０のときの排気酸素濃度曲線である。式（２０）から分かる通り、式（２０）に「切片」＝０、かつ「係数」＝０を代入すると、「時定数」＝０となる。「時定数」＝０の場合、式（１０）から分かる通り、ｙ（ｓ）＝ｕ（ｓ）となる。つまり、図７（Ａ）に示す破線のグラフは、１次遅れ系の入力としての排気酸素濃度曲線と同一である。この１次遅れ系への入力としての排気酸素濃度曲線は、前述の第１ステップで取得した定常状態の排気酸素濃度を、図４に示す運転条件（トルク）の時間的な変化に沿うように繋いだものである。この入力としての排気酸素濃度曲線は、前述の通り、本発明にいう第２の排気酸素濃度曲線の一例に相当する。
【０１１４】
１次遅れ系の出力としての排気酸素濃度曲線は、図７（Ａ）〜（Ｅ）に示すように、図６の二次元平面上の座標点によってその形状が異なる。
【０１１５】
この１次遅れ系の出力としての排気酸素濃度曲線は、入力としての排気酸素濃度曲線から、下記の式（２２）の演算によって算出される。
【０１１６】
【数３２】
0006579975.tif 000033
ここで、ｋは、時間軸上のサンプリング点を表わしており、ｘk，ｙk，Ｔkは、サンプリング点ｋにおける、それぞれ１次遅れ系への入力としての排気酸素濃度、その１次遅れ系からの出力としての排気酸素濃度、およびサンプリング点ｋにおける時定数である。ｘk-1，ｙk-1，Ｔk-1は、サンプリング点ｋよりも１つ前のサンプリングｋ−１における各値である。また、ΔＴは、サンプリング間隔（サンプリング点ｋ−１とｋとの間の時間間隔）を表わしている。図７（Ａ）〜（Ｅ）に示す破線の各グラフは、この式（２２）の演算式に従って算出した排気酸素濃度曲線である。
【０１１７】
ここでは、設定した４隅の座標点Ａ，Ｂ，Ｃ，Ｄで囲まれた領域内における、座標点Ｅ、すなわち、実測値としての排気酸素濃度曲線からの誤差が最小となる「切片」および「係数」が探索される。
【０１１８】
ここでは、実測値としての排気酸素濃度曲線（実線）からの誤差の大きさは、決定係数Ｒ2によって評価している。
【０１１９】
座標点Ｅを求めるための探索は、特定の手法に限定されるものではないが、例えば、座標点Ａ，Ｂ，Ｃ，Ｄで囲まれた領域を複数の小領域に分けて各小領域内の代表点について誤差の大きさを評価して誤差の小さい小領域を見つけ、今度はその小領域をさらに複数の小領域に分けて誤差の小さい小領域を見つけるという探索手法を採用してもよい。あるいは、前掲の非特許文献１に記載されている「実数値遺伝アルゴリズム」を採用して座標点Ｅを探索してもよい。あるいは、最も単純には、座標点Ａ〜Ｄで囲まれた領域内の各点について誤差を評価することで、誤差が最小となる座標点Ｅを捜してもよい。
【０１２０】
このようにして、図７（Ｅ）に示すように、誤差が最小となる「切片」および「係数」を見つけ、それらのパラメータに従って時間的に変化する時定数を持つ１次遅れ系の排気酸素濃度、およびその排気酸素濃度から算出される吸気酸素濃度を含んで記述したエンジンモデルが作成される。
【０１２１】
図８は、排気に含まれるＮＯｘの変化曲線を示した図である。この図８の横軸は時間［ｓ］、縦軸はＮＯｘの量である。ここでは、法定のＪＥ０５の運転パターンに則して運転条件を時間的に変化させている。
【０１２２】
図８（Ａ），（Ｂ），（Ｃ）はそれぞれ、実測値、定常運転のデータから作成したエンジンモデル（定常モデル）を使った算出値、および過渡運転時のデータを使って上記の方法により作成したエンジンモデル（過渡モデル）を使った算出値である。
【０１２３】
また、この図８における、（Ａ），（Ｂ），（Ｃ）に跨る一点鎖線の囲いは、アイドリング開始時、すなわち、ある回転速度でアクセルを離し、回転速度が低下してきてアイドリングに入ったタイミングのデータである。ここでは、アイドリングに入った後、もう１度アクセルを踏んでから離し、再びアイドリングに入るように、２度繰り返している。ここでの実験等で使ったディーゼルエンジンは、アイドリング開始時に図８（Ａ）に示すようにＮＯｘの量が増加する傾向にあったが、（Ｂ）の定常モデルでは、そのＮＯｘの増加がうまく再現できていない。これに対し、図８（Ｃ）の過渡モデルの場合、吸気酸素濃度が高い状態で燃料噴射を開始していることがシミュレーションに組み込まれ、ＮＯｘが増加する現象が再現されていることが分かる。
【０１２４】
なお、本実施形態では、過渡運転にあたり、図４に示すように、回転速度を一定に保ちつつトルクを時間的に変化させているが、それとは逆に、トルクを一定に保ちつつ回転速度を時間的に変化させてもよく、あるいは回転速度とトルクとの双方を変化させる運転条件を採用してもよい。
【符号の説明】
【０１２５】
１  エンジンモデル作成システム
１０  エンジン
１１  排気管
１２  回転軸
２０  ダイナモ
３０  制御・計測用パソコン
４０  スロットルアクチュエータ
５０  ＥＣＵ
６０  排ガス分析計
７０  ダイナモ制御盤
８１  温度センサ
８２  流量センサ
８３  トルクセンサ
８４  回転速度センサ
【図１】
0006579975.tif 000037
【図２】
0006579975.tif 000038
【図３】
0006579975.tif 000039
【図４】
0006579975.tif 000040
【図５】
0006579975.tif 000041
【図６】
0006579975.tif 000042
【図７】
0006579975.tif 000043
【図８】
0006579975.tif 000044