特許 6472371 エンジン特性試験方法およびエンジン特性試験システム

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】6472371

(24)【登録日】2019年2月1日

(45)【発行日】2019年2月20日

(54)【発明の名称】エンジン特性試験方法およびエンジン特性試験システム

(51)【国際特許分類】

   G01M 15/10 20060101AFI20190207BHJP

   F02D 45/00 20060101ALI20190207BHJP

【ＦＩ】

   G01M15/10

   F02D45/00 372Z

【請求項の数】4

【全頁数】12

(21)【出願番号】特願2015-239142(P2015-239142)

(22)【出願日】2015年12月8日

(65)【公開番号】特開2017-106763(P2017-106763A)

(43)【公開日】2017年6月15日

【審査請求日】2018年7月3日

(73)【特許権者】

【識別番号】000145806

【氏名又は名称】株式会社小野測器

(74)【代理人】

【識別番号】100094330

【弁理士】

【氏名又は名称】山田 正紀

(74)【代理人】

【識別番号】100109689

【弁理士】

【氏名又は名称】三上 結

(72)【発明者】

【氏名】村瀬 道夫

(72)【発明者】

【氏名】翁 志強

(72)【発明者】

【氏名】三瓶 祐一郎

【審査官】 小野 郁磨

(56)【参考文献】

【文献】 特開２００３−０１３７９４（ＪＰ，Ａ）

【文献】 特開２０１４−２２４７８６（ＪＰ，Ａ）

【文献】 特開２００９−９９０５０（ＪＰ，Ａ）

【文献】 特開２００９−２１０４２６（ＪＰ，Ａ）

【文献】 特開２００９−３６０５１（ＪＰ，Ａ）

【文献】 特開２００９−２５８０６８（ＪＰ，Ａ）

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ｇ０１Ｍ １５／００−１５／１０

Ｆ０２Ｄ ４５／００

(57)【特許請求の範囲】

【請求項1】

エンジンの回転速度とトルクとを変数とする二次元平面上の離散的な複数の巡回点を、同一の巡回点の複数回の通過を許容して順次辿る巡回経路を生成するにあたって、
試験対象のエンジンの回転速度とトルクとを説明変数とし、該エンジンの回転速度とトルクの変化の履歴に応じて定まる排気温度を目的変数とするエンジンモデルを参照するとともに、
同一の巡回点を複数回通過したときの該同一の巡回点における排気温度のうちの最高温度と最低温度との差分が大きく、かつ、巡回点から巡回点への1回の移動を単位巡回経路と呼んだとき、単位巡回経路の始点である移動前巡回点と単位巡回経路の終点である移動後巡回点とがともに一致する単位巡回経路を複数回通過した場合であっても通過した単位巡回経路の数をその単位巡回経路について１つとして計数したときの通過した単位巡回経路の数が多いことをもって高い評価を与える評価関数を用いて、
巡回経路の生成と、前記評価関数を用いた評価による淘汰とを繰り返すことにより、該評価関数を用いた評価の高い高評価巡回経路を生成する第１ステップと、
試験対象のエンジンの回転速度とトルクが前記高評価巡回経路を辿るようにエンジン制御パラメータを調整することによって該エンジンの運転を制御して、該エンジンの排気成分と成分量とを含むエンジン特性を測定する第２ステップとを有することを特徴とするエンジン特性試験方法。

【請求項2】

前記巡回経路が、同一の巡回点を複数回通過した場合であっても１回の通過ごとに、巡回点の通過数を累積加算したときの巡回点の数が予め定められている巡回経路であることを特徴とする請求項１に記載のエンジン特性試験方法。

【請求項3】

前記第１ステップが、前記巡回経路上の巡回点を遺伝子として捉えた遺伝アルゴリズムを利用して該巡回経路の生成と淘汰を繰り返すステップであることを特徴とする請求項１または２に記載のエンジン特性試験方法。

【請求項4】

エンジンの回転速度とトルクとを変数とする二次元平面上の離散的な複数の巡回点を、同一の巡回点の複数回の通過を許容して順次辿る巡回経路を生成するにあたって、
試験対象のエンジンの回転速度とトルクとを説明変数とし、該エンジンの回転速度とトルクの変化の履歴に応じて定まる排気温度を目的変数とするエンジンモデルを参照するとともに、
同一の巡回点を複数回通過したときの該同一の巡回点における排気温度のうちの最高温度と最低温度との差分が大きく、かつ、巡回点から巡回点への1回の移動を単位巡回経路と呼んだとき、単位巡回経路の始点である移動前巡回点と単位巡回経路の終点である移動後巡回点とがともに一致する単位巡回経路を複数回通過した場合であっても通過した単位巡回経路の数をその単位巡回経路について１つとして計数したときの通過した単位巡回経路の数が多いことをもって高い評価を与える評価関数を用いて、
巡回経路の生成と、前記評価関数を用いた評価による淘汰とを繰り返すことにより、該評価関数を用いた評価の高い高評価巡回経路を生成する演算装置と、
試験対象のエンジンに連結されるダイナモを備え、該エンジンの回転速度とトルクが前記高評価巡回経路を辿るようにエンジン制御パラメータを調整することによって該エンジンの運転を制御する制御装置と、
前記エンジンの排気成分と成分量とを含むエンジン特性を測定する測定装置とを備えたことを特徴とするエンジン特性試験システム。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、エンジンの様々な運転条件と排気ガス中のＮＯｘや燃費などのエンジン特性との関係を表わすエンジンモデルを作成するためのエンジン特性試験方法およびエンジンモデル試験システムに関する。

【背景技術】

【0002】

近年の自動車の普及と環境科学の進歩に伴い、燃費や排ガスの規制が激しく強化されてきている。このため、精密なエンジンモデルを作成し、エンジンの様々な運転条件下での排気ガス中のＮＯｘや燃費などを高精度にシミュレーションすることが求められている。

【0003】

ここで、エンジンモデルの作成にあたっては、エンジンにダイナモを連結した実験システム上で実験を行ない、その実験結果に基づいてエンジンモデルを作成することが従来より行なわれている（特許文献１参照）。

【0004】

エンジンモデルを作成して、自動車に搭載されたエンジンの制御パラメータをその作成されたエンジンモデルに基づいて調整することで、エンジンの回転速度やトルクを所要の値に保ちながら排気ガス中のＮＯｘの低減や燃費の向上を図る努力がなされている。

【0005】

従来は、エンジンの回転速度とトルクとからなる２次元平面上の各点での定常運転によるエンジンモデルが作成され、活用されている。しかしながら、この定常運転によるエンジンモデルの場合、エンジンの回転速度やトルクが変化しつつある過渡運転時、例えば運転開始時などについては、十分に最適化されたエンジンモデルにはなっていない。すなわち、エンジンモデルから予測されるＮＯｘの値や燃費等のエンジン特性値が、エンジンを実際に自動車に搭載して運転して実測したときのエンジン特性値に高精度に近似しているとは言い難く、それらの間の誤差が十分に小さいとは言い難い。

【0006】

この誤差が生じる大きな原因の１つは、過渡運転時の排気温度がエンジンモデルに正しく組み込まれていないことにある。すなわち、定常運転によるエンジンモデルの場合、排気温度は、エンジンの現在の回転速度やトルクによって定まるのに対し、過渡運転の場合、排気温度は、エンジンの現在の回転速度やトルクだけでは定まらず、現在に至るまでの回転速度やトルクの変化の影響が時間的に遅れて現在の排気温度に反映される。

【0007】

例えば、排気温度が低いときは、エンジンの壁温も低いことが多く、再燃焼が妨げられることによる煙の増加や、燃焼温が下がることによるＮＯｘの減少が起こる。また、例えば排気温度が低いときは、媒体の効きが悪くなり、未燃炭化水素（ＨＣ）等の除去率が低下する。したがって、正確なシミュレーションを行なうためには、排気温度が大きく変化するような様々な過渡運転を行なって、ＮＯｘや燃費等のエンジン特性を測定し、その測定値に基づいて高精度なエンジンモデルを作成する必要がある。

【0008】

なお、後述の実施形態の説明のために非特許文献１を挙げておく。

【先行技術文献】

【特許文献】

【0009】

【特許文献1】特開２００９−２１０４２６号公報

【非特許文献】

【0010】

【非特許文献1】自動車技術会論文集 Ｖｏｌ．４５，Ｎｏ．６ ２０１４４７８４ 「ディーゼルエンジンの排気マニフォールド内ガス温度を推定するバーチャルセンサの開発」

【発明の概要】

【発明が解決しようとする課題】

【0011】

本発明は、上記事情に鑑み、高精度なエンジンモデル作成のために好適なエンジン特性試験方法およびエンジン特性試験システムを提供することを目的とする。

【課題を解決するための手段】

【0012】

上記目的を達成する本発明のエンジン特性試験方法は、
エンジンの回転速度とトルクとを変数とする二次元平面上の離散的な複数の巡回点を、同一の巡回点の複数回の通過を許容して順次辿る巡回経路を生成するにあたって、
試験対象のエンジンの回転速度とトルクとを説明変数とし、エンジンの回転速度とトルクの変化の履歴に応じて定まる排気温度を目的変数とするエンジンモデルを参照するとともに、
同一の巡回点を複数回通過したときの同一の巡回点における排気温度のうちの最高温度と最低温度との差分が大きく、かつ、巡回点から巡回点への1回の移動を単位巡回経路と呼んだとき、単位巡回経路の始点である移動前巡回点と単位巡回経路の終点である移動後巡回点とがともに一致する単位巡回経路を複数回通過した場合であっても通過した単位巡回経路の数をその単位巡回経路について１つとして計数したときの通過した単位巡回経路の数が多いことをもって高い評価を与える評価関数を用いて、
巡回経路の生成と、上記評価関数を用いた評価による淘汰とを繰り返すことにより、その評価関数を用いた評価の高い高評価巡回経路を生成する第１ステップと、
試験対象のエンジンの回転速度とトルクが上記高評価巡回経路を辿るようにエンジン制御パラメータを調整することによってエンジンの運転を制御して、エンジンの排気成分と成分量とを含むエンジン特性を測定する第２ステップとを有することを特徴とする。

【0013】

本発明のエンジン特性試験方法によれば、上記の第１ステップにおいて、同じ回転速度および同じトルクであっても、それ以前の運転条件の履歴に応じた様々な排気温度を得る高評価循回経路が定められる。そして、第２ステップにおいて、その高評価巡回経路を辿るようにエンジン制御パラメータを調整することによってエンジンの運転を制御し、幅広い様々な排気温度下でのエンジン特性を測定する。こうすることで、過渡運転を含む幅広い運転条件下での実運転を高精度にシミュレートしたエンジンモデルを作成することができる。

【0014】

ここで、本発明のエンジン特性試験方法において、上記巡回経路が、巡回点を１つ通過するごとに、同一の巡回点を複数回通過した場合であっても１回の通過ごとに、巡回点の通過数を累積加算したときの巡回点の数が予め定められている巡回経路であることが好ましい。

【0015】

巡回経路の長さ、すなわち１つの巡回経路を構成する巡回点の数を予め定めておくことにより、高評価巡回経路を算出するまでの無駄な演算が抑えられる。ここで、巡回経路の長さは、唯一に定められていてもよく、長さＡ〜長さＢの間の長さ、というようにある範囲を持って定められていてもよい。

【0016】

また、本発明のエンジン特性試験装置において、上記第１ステップが、前記巡回経路上の巡回点を遺伝子として捉えた遺伝アルゴリズムを利用して巡回経路の生成と淘汰を繰り返すステップであることが好ましい。

【0017】

遺伝アルゴリズムを適用すると、高評価巡回経路を容易に算出することができる。

【0018】

また、上記目的を達成する本発明のエンジン特性試験システムは、
エンジンの回転速度とトルクとを変数とする二次元平面上の離散的な複数の巡回点を、同一の巡回点の複数回の通過を許容して順次辿る巡回経路を生成するにあたって、
試験対象エンジンの回転速度とトルクとを説明変数とし、そのエンジンの回転速度とトルクの変化の履歴に応じて定まる排気温度を目的変数とするエンジンモデルを参照するとともに、
同一の巡回点を複数回通過したときの同一の巡回点における排気温度のうちの最高温度と最低温度との差分が大きく、かつ、巡回点から巡回点への1回の移動を単位巡回経路と呼んだとき、単位巡回経路の始点である移動前巡回点と単位巡回経路の終点である移動後巡回点とがともに一致する単位巡回経路を複数回通過した場合であっても通過した単位巡回経路の数をその単位巡回経路について１つとして計数したときの通過した単位巡回経路の数が多いことをもって高い評価を与える評価関数を用いて、
巡回経路の生成と、上記評価関数を用いた評価による淘汰とを繰り返すことにより、その評価関数を用いた評価の高い高評価巡回経路を生成する演算装置と、
試験対象エンジンに連結されるダイナモを備え、そのエンジンの回転速度とトルクが上記高評価巡回経路を辿るようにエンジン制御パラメータを調整することによってエンジンの運転を制御する制御装置と、
エンジンの排気成分と成分量とを含むエンジン特性を測定する測定装置とを備えたことを特徴とする。

【発明の効果】

【0019】

以上の本発明によれば、高精度のエンジンモデル作成のためのエンジン特性試験方法およびエンジン特性試験システムが実現する。

【図面の簡単な説明】

【0020】

【図1】本発明の一実施形態としてのエンジン特性試験システムの構成図である。

【図2】ＪＥ０５モード走行パターンを示した図である。

【図3】運転パターンの概念図である。

【図4】エンジン回転速度とトルクと排気温度の時間変化の一例を示した図である。

【図5】回転速度とトルクとからなる２次元平面を示した模式図である。

【図6】ＪＥ０５モード走行時の、ＮＯｘ排出予測精度を示した図である。

【発明を実施するための形態】

【0021】

以下、本発明の実施の形態を説明する。

【0022】

図１は、本発明の一実施形態としてのエンジン特性試験システムの構成図である。

【0023】

ここに示すエンジン特性試験システム１は、エンジン特性試験対象のエンジン１０に連結されたダイナモ２０を備えている。また、このエンジン特性試験システム１は、制御・計測用パソコン３０、スロットルアクチュエータ４０、ＥＣＵ（エンジンコントロールユニット）５０、排ガス分析計６０、およびダイナモ制御盤７０を備えている。また、ここには、エンジン１０の排気管１１内を流れる排ガスの排気温度を測定する温度センサ８１や、排ガスの流量を測定する流量センサ８２や、エンジン１０の回転軸１２のトルクを測定するトルクセンサ８３や、回転軸の回転速度を測定する回転速度センサ８４などの様々なセンサが備えられている。

【0024】

制御・計測用パソコン３０は、スロットルアクチュエータ４０を介してのエンジン１０のスロットル（不図示）の開度の制御によるトルク制御、ダイナモ制御盤７０を介してのダイナモ２０の制御によるエンジン１０の回転軸１２の回転速度制御、ＥＣＵ５０を介しての噴射進角の制御等を行なう。また、この制御・計測用パソコン３０では、エンジン１０の回連速度、トルク、排気温度、排気流量等のデータの収集が行なわれる。排ガス分析計６０では、排気管１１内を流れる排ガスの成分分析が行なわれ、制御・計測用パソコン３０には、その分析結果も収集される。さらに、この制御・計測用パソコン３０には、ＥＣＵ５０によるエンジン１０の制御パラメータの値も収集される。さらに、この制御・計測用パソコン３０は、燃料の供給、排ガスの処理等を行なう様々な環境設備９０の制御も行なう。

【0025】

この図１に示すエンジン特性試験システム１において、ダイナモ２０、スロットルアクチュエータ４０、ＥＣＵ５０、およびダイナモ制御盤７０や、制御・計測用パソコン３０のそれらを制御する機能の組合せが、本発明にいう制御装置の一例に相当する。

【0026】

また、この図１に示すエンジン特性試験システム１のうちの、エンジン１０の回転速度やトルクなどを測定する各種センサと、排ガス分析計６０と、制御・計測用パソコン３０のデータ収集機能との組合せが、本発明にいう測定装置の一例に相当する。

【0027】

また、この図１に示すエンジン特性試験システム１の制御・計測用パソコン３０は、
エンジン制御装置にエンジンをどのように制御されるかという運転条件（後で説明する高評価巡回経路）を算出する機能を有している。制御・計測用パソコン３０の、この算出機能は、本発明にいう演算装置の一例に相当する。

【0028】

図２は、ＪＥ０５モード走行パターンを示した図である。

【0029】

この図２の横軸は時間［ｓ］、縦軸は車速（ｋｍ／ｈ）である。

【0030】

例えばこのＪＥ０５モード走行によりＮＯｘや燃費等が計測される。

【0031】

ここで、ポイントＡとポイントＢではいずれも停車した状態から走行を開始している。しかしながら、ポイントＡでは停車直後に走行を開始し、ポイントＢでは暫く停車した後走行を開始している。このため、ポイントＡではエンジンが暖まった状態で走行が開始され、ポイントＢではエンジンが冷えた状態で走行が開始される。この温度条件の違いが排気ガス中のＮＯｘの量や燃費等に大きく影響を与える。したがって幅広い温度条件を再現させたエンジン特性試験を行なってデータを収集し、高精度なエンジンモデルを作成する必要がある。

【0032】

図３は、運転パターンの概念図である。

【0033】

ここには、回転速度［ｒ／ｍｉｎ］とトルク［Ｎ・ｍ］との２軸からなる平面に垂直に排気温度［℃］の軸が示されている。すなわち、ここには、回転速度とトルクと排気温度とからなる３次元空間が示されている。

【0034】

この３次元空間に示されているグラフＡは、エンジンが辿る、回転速度とトルクと排気温度からなる巡回経路を示している。

【0035】

図１に示すエンジン特性試験システム１を使ってエンジン１０の運転を制御するにあたり、直接的には、エンジンの回転速度とトルクが制御対象となる。これら回転速度とトルクとが辿る経路に応じて排気温度が上昇、下降し、例えばこの図３に示すような巡回経路を辿ることになる。

【0036】

図４は、エンジン回転速度とトルクと排気温度の時間変化の一例を示した図である。図４（Ａ），（Ｂ），（Ｃ）の縦軸は、それぞれ、回転速度、トルク、および排気温度である。横軸は、いずれも時間［ｓ］である。

【0037】

回転速度およびトルクを図４（Ａ），（Ｂ）に示すように、時間的に変化させると、それに伴って図４（Ｃ）に示すように排気温度が変化する。

【0038】

そこで、ここでは、以下のようにして、排気温度ができるだけ幅広く変化する巡回経路を算出する。

【0039】

図５は、回転速度とトルクとからなる２次元平面を示した模式図である。この図５の横軸は、回転速度［ｒ／ｍｉｎ］、縦軸は、トルク［Ｎ・ｍ］である。

【0040】

ここでは、この２次元平面上であって、試験対象のエンジン１０での運転が可能な領域内に複数の巡回点を設定する。この図５では、巡回点Ａ〜巡回点Ｕが設定されている。そして、以下のようにして、これら複数の巡回点を、同一の巡回点の複数回の通過を許容して順次辿る巡回経路を生成する。

【0041】

ここでは、エンジンモデルが参照される。このエンジンモデルは、試験対象のエンジン１０（図１参照）の回転速度とトルクとを説明変数とし、そのエンジンの回転速度とトルクの変化の履歴に応じて定まる排気温度を目的変数とするエンジンモデルである。

【0042】

エンジンの回転速度とトルクの変化の履歴に応じて定まる排気温度は、例えば、定常運転時に測定した排気温度を１次遅れ系に入力して実運転時の排気温度の変化に適合するように時定数を定めることなどにより算出することができる。あるいは、前掲の非特許文献１に記載された演算方法を用いて算出してもよい。

【0043】

また、ここでは、生成された巡回経路を評価する評価関数が用いられる。この評価関数は、
（１）同一の巡回点を複数回通過したときのその同一の巡回点における排気温度のうちの最高温度と最低温度との差分が大きいこと、および
（２）同一の巡回点を複数回通過した場合であっても通過した巡回点の数をその巡回点について１つとして計数したときの通過した巡回点の数が多いこと、
を評価基準とし、（１）と（２）の双方の評価が高いことをもって高い評価を与える評価関数である。すなわち、同じ巡回点、つまり同じ回転速度および同じトルクにおける、排気温度の変化幅が大きいことと、様々な回転速度とトルクにおける試験が行なわれることを評価の対象とする。

【0044】

本実施形態に置いては、巡回経路の長さが予め定められていて、その長さが定められている巡回経路が評価の対象となる。ここで、巡回経路の長さとは、同一の巡回点を複数回通過した場合であっても１回の通過ごとに巡回点の通過数を累積加算した巡回点の数をいう。

【0045】

図１に示す制御・計測用パソコン３０では、上記の、予め長さが定められている巡回経路について、巡回経路の生成と、上記の評価関数を用いた評価による淘汰とを繰り返すことにより、その評価関数による評価の高い高評価巡回経路が生成される。

【0046】

そして、この高評価巡回経路を生成した後、図１に示すエンジン特性試験システム１では、エンジン特性試験対象のエンジン１０の回転速度とトルクが、その生成した高評価巡回経路を辿るようにエンジン制御パラメータを調整することによって該エンジンの運転を制御して、その運転時におけるエンジンの排気成分と成分量とを含むエンジン特性が測定される。

【0047】

そしてさらに、制御・計測用パソコン３０では、その測定データを用いて高精度なエンジンモデルが生成される。

【0048】

ここで、巡回経路の生成と淘汰を繰り返して高評価巡回経路を生成するにあたり、本実施形態では、巡回経路上の巡回点を遺伝子として捉えた、以下の遺伝アルゴリズムを利用している。

【0049】

ここでは、初期値として、
「ＡＣＤＪＰＯＮＭＬＧＣＡ」 ・・・（１）
「ＡＤＥＦＫＪＮＲＭＨＣＡ」 ・・・（２）
なる２つの巡回経路（ＤＮＡ）を考える。

【0050】

なお、実際には、図５に示す約２０の巡回点に対し、巡回経路の長さ、すなわち１つの巡回経路を構成する巡回点の数は例えば１００である。

【0051】

ここでは、上記の（１），（２）のＤＮＡから出発し、ＤＮＡの突然変異や交差等によりＤＮＡを増殖させる。例えば、上記（１）のＤＮＡを構成する遺伝子「・・・ＪＰＯ・・・」のＰをＫに突然変異させて「・・・ＪＫＯ・・・」とし、
新たなＤＮＡ
「ＡＣＤＪＫＯＮＭＬＧＣＡ」 ・・・（３）
を生成する。ただし、この突然変異にあたっては、図５に示す経路が繋がるように突然変異させる。

【0052】

また、上記（１）のＤＮＡの前半部分「ＡＣＤＪＰ・・・」と上記（２）のＤＮＡの後半部分「・・・ＮＲＭＨＣＡ」とを繋げて
「ＡＣＤＪＰＯＮＲＭＨＣＡ」 ・・・（４）
を生成する。

【0053】

このようにしてＤＮＡの数を増殖させ、評価関数にあてはめて評価の低いＤＮＡを取り除く。この評価関数としては、
Ε＝α・Σ（最高温度―最低温度）＋β・通過した単位巡回経路の数 ・・・（５）
が採用される。

【0054】

ここで、α，βは、経験的に定めた各定数である。

【0055】

また、（最高温度―最低温度）は、同じ巡回点を複数回通過したときのその巡回点における最高の排気温度と最低の排気温度との差分である。１回しか通過しなかった巡回点については、その差分はゼロとなる。その差分を、通過した全ての巡回点について加算する。

【0056】

また、「通過した単位巡回経路の数」は、同じ単位巡回経路を複数回通過しても、その単位巡回経路についてはその複数回を持って通過した単位巡回経路の数を１つと数えたときの単位巡回経路の数である。

【0057】

このようにして、ＤＮＡの増殖と評価関数による淘汰を、例えば１０００世代にわたって繰り返し、最終的に上記の評価関数において評価の高い高評価ＤＮＡ（高評価巡回回路）を生成する。

【0058】

この高評価巡回経路を生成した後は、エンジン特性試験対象のエンジン１０の回転速度とトルクが、その生成した高評価巡回経路を辿るようにエンジン制御パラメータを調整することによって該エンジンの運転を制御して、エンジンの排気成分と成分量とを含むエンジン特性が測定される。そして、この測定データに基づいて、さらに高精度なエンジンモデルが生成される。

【0059】

図６は、ＪＥ０５モード走行時の、ＮＯｘ排出予測精度を示した図である。

【0060】

この図６は、（Ａ）〜（Ｄ）のいずれにおいても横軸は時間［ｓ］である。

【0061】

また、図６（Ａ），（Ｂ）の縦軸は、それぞれ、回転速度とトルクである。

【0062】

また、図６（Ｃ），（Ｄ）の縦軸は、いずれもＮＯｘの量を示している。

【0063】

図６（Ｃ），（Ｄ）において、実線は実測値を示していて、（Ｃ）と（Ｄ）とで同一のカーブとなっている。また、図６（Ｃ）の破線は、上記のような高評価巡回経路を算出することなく、適宜行なったエンジン試験で集収したデータから生成したエンジンモデルを使って算出した値である。これに対し、図６（Ｄ）の破線は、上記のようにして算出した高評価巡回経路を辿るエンジン試験で集収したデータから生成したエンジンモデルを使った算出値である。

【0064】

図６（Ｄ）の、一点鎖線で囲った部分は、暖気不十分なタイミングを表わしている。図６（Ｃ）と図６（Ｄ）とを比べると、図６（Ｄ）の方が、この暖気不十分の時のＮＯｘ排出量の過大評価が抑えられていることが分かる。

【0065】

なお、上記の実施形態では、１つの巡回経路の長さを一つに定めているが、長さＡ〜長さＢの範囲内というように長さの幅を持たせておいてもよい。あるいは、長さの上限値だけ決めておいて、長さの下限値については特に定めなくてもよい。

【0066】

また、上記の実施形態では、高評価巡回経路を生成するにあたり、遺伝アルゴリズムを使っているが、高評価巡回経路を生成するためのアルゴリズムはどのようなものであってもよく、特定のアルゴリズムに限定されるものではない。

【0067】

さらに、上記の実施形態では、評価関数として式（５）を採用しているが、式（５）の評価関数に限られるものではなく、同一の巡回点を複数回通過したときのその同一の巡回点における排気温度のうちの最高温度と最低温度との差分が大きく、かつ、巡回点から巡回点への1回の移動を単位巡回経路と呼んだとき、単位巡回経路の始点である移動前巡回点と単位巡回経路の終点である移動後巡回点とがともに一致する単位巡回経路を複数回通過した場合であっても通過した単位巡回経路の数をその単位巡回経路について１つとして計数したときの通過した単位巡回経路の数が多いという条件を満たす評価関数であればよい。

【符号の説明】

【0068】

１ エンジン特性試験システム
１０ エンジン
１１ 排気管
１２ 回転軸
２０ ダイナモ
３０ 制御・計測用パソコン
４０ スロットルアクチュエータ
５０ ＥＣＵ
６０ 排ガス分析計
７０ ダイナモ制御盤
８１ 温度センサ
８２ 流量センサ
８３ トルクセンサ
８４ 回転速度センサ

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】