(57)【要約】【課題】エンジンの制御パラメータを適合する際に用いるエンジン特性モデルの精度向上（制御パラメータの適合値の精度向上）と工数削減とを両立させる。【解決手段】適合対象となる制御パラメータに対する物理パラメータを選択する（１０１）。適合対象となる制御パラメータがＶＣＴ進角値であれば、それに対する物理パラメータとして、筒内ＥＧＲ率、筒内流速、吸気温度、ポンピングロス、吸気管圧力、実圧縮比の中から選択し、適合対象となる制御パラメータが噴射時期であれば、それに対する物理パラメータとして、噴霧移動距離、霧化時間、蒸発燃料量、噴射時筒内流速の中から選択する。次に、制御パラメータと物理パラメータとの関係を計測データにより算出し（１０２）、制御パラメータの実験計画範囲の境界を定める物理パラメータの判定閾値を生じさせる制御パラメータの値を算出して（１０３）、制御パラメータの実験計画範囲を決定する（１０４）。【選択図】図３
2012013637.tif 000002
【特許請求の範囲】
【請求項１】
計測対象となるエンジンの制御パラメータを変化させて当該制御パラメータの変化に応じて変化する物理量（以下「物理パラメータ」という）を所定数の計測点で計測し、前記所定数の計測点での計測結果に基づいて前記制御パラメータと前記物理パラメータとの関係をモデル化してエンジン特性モデルを作成し、前記エンジン特性モデルを用いて前記制御パラメータの適合値を前記エンジンの要求性能が満たされるように決定するエンジン制御パラメータの適合方法において、
前記制御パラメータを変化させる計測点を配置する実験計画範囲を燃焼に寄与する物理パラメータで設定することを特徴とするエンジン制御パラメータの適合方法。
【請求項２】
請求項１に記載のエンジン制御パラメータの適合方法において、
前記制御パラメータの実験計画範囲の境界を定める物理パラメータの判定閾値を当該物理パラメータが安定燃焼限界を越えないように決定することを特徴とするエンジン制御パラメータの適合方法。
【請求項３】
請求項１に記載のエンジン制御パラメータの適合方法において、
単一の制御パラメータを変化させて当該単一の制御パラメータと物理パラメータとの関係を求め、複数の制御パラメータを同時に変化させる場合の実験計画範囲を前記単一の制御パラメータと物理パラメータとの関係を用いて設定することを特徴とするエンジン制御パラメータの適合方法。
【請求項４】
請求項２に記載のエンジン制御パラメータの適合方法において、
複数の制御パラメータの実験計画範囲の境界を定める物理パラメータの判定閾値を当該物理パラメータが安定燃焼限界を越えないように決定する際に、いずれかの制御パラメータの可動範囲の限界でも前記物理パラメータが前記安定燃焼限界を越えない場合は、当該制御パラメータを可動範囲の限界に固定した状態で他の制御パラメータを前記物理パラメータが前記安定燃焼限界を越えるまで変化させて前記複数の制御パラメータの実験計画範囲を設定することを特徴とするエンジン制御パラメータの適合方法。
【請求項５】
請求項１に記載のエンジン制御パラメータの適合方法において、
前記制御パラメータの適合値を決定する際に前記エンジン特性モデルと前記制御パラメータの実験計画範囲を考慮して当該制御パラメータの適合値を決定することを特徴とするエンジン制御パラメータの適合方法。
【請求項６】
請求項１に記載のエンジン制御パラメータの適合方法において、
代表運転条件での計測結果に基づいて当該代表運転条件での制御パラメータの実験計画範囲を設定すると共に、運転条件の変化による物理パラメータの変化を予測し、前記代表運転条件での制御パラメータの実験計画範囲と前記運転条件の変化による物理パラメータの変化の予測結果に基づいて他の運転条件での制御パラメータの実験計画範囲を設定することを特徴とするエンジン制御パラメータの適合方法。
【請求項７】
計測対象となるエンジンの制御パラメータを変化させて当該制御パラメータの変化に応じて変化する物理量（以下「物理パラメータ」という）を所定数の計測点で計測し、前記所定数の計測点での計測結果に基づいて前記制御パラメータと前記物理パラメータとの関係をモデル化してエンジン特性モデルを作成し、前記エンジン特性モデルを用いて前記制御パラメータの適合値を前記エンジンの要求性能が満たされるように決定するエンジン制御パラメータの適合装置において、
前記制御パラメータを変化させる計測点を配置する実験計画範囲を燃焼に寄与する物理パラメータで設定する実験計画範囲設定手段を備えていることを特徴とするエンジン制御パラメータの適合装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、エンジン特性モデルを用いてエンジンの制御パラメータの適合値をエンジンの要求性能が満たれるように決定するエンジン制御パラメータの適合方法及び適合装置に関する発明である。
【背景技術】
【０００２】
近年の高性能エンジンは、出力向上、排気エミッション低減、燃費節減等を実現するために、可変バルブタイミング機構やＥＧＲシステム等の様々な機能を搭載しているため、適合すべき制御パラメータが噴射時期や点火時期のみではなく、バルブタイミングやＥＧＲ率等も適合する必要があり、適合すべき制御パラメータの数が増加して、制御パラメータの適合作業が非常に面倒なものとなってきている。
【０００３】
そこで、特許文献１（特開２００２−２０６４５６号公報）や、特許文献２（特開２００５−４２６５６号公報）に記載されているように、計測対象となるエンジンの制御パラメータを変化させて所定数の計測点でエンジン特性値（物理パラメータ）を計測して、その計測結果に基づいて各制御パラメータとエンジン特性値との関係をモデル化してエンジン特性モデルを作成し、このエンジン特性モデルを用いて制御パラメータの適合値を算出するようにしたものがある。
【０００４】
しかし、上記特許文献１，２の技術では、エンジンの制御パラメータを変化させる計測点が安定燃焼領域内に配置されるとは限らず、燃焼悪化領域に配置される可能性があり、燃焼悪化領域では、失火や燃焼不安定になったり、異常燃焼によってエンジンが損傷する可能性がある。燃焼悪化領域に計測点が配置されると、当該計測点の計測精度が悪化したり、計測不能となるため、エンジン特性モデルの精度が悪化して、当該エンジン特性モデルを用いて算出する制御パラメータの適合値の精度も悪化する。
【０００５】
そこで、特許文献３（特開２００４−２６３６８０号公報）に記載されているように、エンジン特性の計測データに基づいてエンジンの制御パラメータと着火遅れ・燃焼期間との関係と、着火遅れ・燃焼期間と燃焼安定性との関係をモデル化して失火推定モデルを作成し、この失火推定モデルを用いて失火領域を推定した後、推定した失火領域を除いた燃焼可能領域内に必要数の計測点を配置するようにしたものがある。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２００２−２０６４５６号公報
【特許文献２】特開２００５−４２６５６号公報
【特許文献３】特開２００４−２６３６８０号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかし、上記特許文献３の技術では、失火推定モデルを作成するのに、エンジン特性モデルを作成するのとほぼ同様の工程が必要となるため、制御パラメータの適合工数が増加して、工数削減の要求を満たすことができない。
【０００８】
そこで、本発明が解決しようとする課題は、エンジン特性モデルの精度向上（制御パラメータの適合値の精度向上）と工数削減とを両立できるエンジン制御パラメータの適合方法及び適合装置を提供することである。
【課題を解決するための手段】
【０００９】
上記課題を解決するために、請求項１に係る発明は、エンジンの制御パラメータを変化させて当該制御パラメータの変化に応じて変化する物理量（以下「物理パラメータ」という）を所定数の計測点で計測し、前記所定数の計測点での計測結果に基づいて前記制御パラメータと前記物理パラメータとの関係をモデル化してエンジン特性モデルを作成し、前記エンジン特性モデルを用いて前記制御パラメータの適合値を前記エンジンの要求性能が満たされるように決定するエンジン制御パラメータの適合方法において、前記制御パラメータを変化させる計測点を配置する実験計画範囲を燃焼に寄与する物理パラメータで設定するようにしたものである。
【００１０】
従来技術では、制御パラメータの実験計画範囲を当該制御パラメータで設定するようにしていたが、制御パラメータは、エンジンを動作させる際の指令値であり、実際の燃焼との関係が複雑である。しかも、燃焼の結果によって得られる物理パラメータの計測値は、ばらつきが大きく、制御パラメータとの関係を正確に定義するには多くの計測工数が必要となる。
【００１１】
これに対し、本発明では、制御パラメータの実験計画範囲を燃焼に寄与する物理パラメータで設定するようにしている。燃焼に寄与する物理パラメータは、燃焼状態を表す指標となるため、燃焼に寄与する物理パラメータを用いて実験計画範囲を設定すれば、比較的少ない計測工数で、実験計画範囲を安定燃焼領域内に精度良く設定することが可能となる。また、制御パラメータは、変化する物理パラメータが共通のものがあるため、変化する物理パラメータが共通の複数の制御パラメータを同時に変化させた場合、制御パラメータの実験計画範囲も複合的に変化することがあるが、複数の制御パラメータに対して共通の物理パラメータを用いて実験計画範囲を設定すれば、比較的少ない計測工数で、複合的な実験計画範囲も安定燃焼領域内に精度良く設定することが可能となる。これにより、エンジン特性モデルの精度向上（制御パラメータの適合値の精度向上）と工数削減とを両立させることができる。
【００１２】
この場合、請求項２のように、制御パラメータの実験計画範囲の境界を定める物理パラメータの判定閾値を当該物理パラメータが安定燃焼限界（安定燃焼領域の限界）を越えないように決定するようにすれば良い。これにより、計測点が安定燃焼限界を越えた領域に配置されることを確実に防止でき、全ての計測点で物理パラメータを精度良く計測することができる。
【００１３】
また、請求項３のように、単一の制御パラメータを変化させて当該単一の制御パラメータと物理パラメータとの関係を求め、複数の制御パラメータを同時に変化させる場合の実験計画範囲を前記単一の制御パラメータと物理パラメータとの関係を用いて設定するようにしても良い。このようにすれば、全ての制御パラメータの組み合わせを計測する必要がなくなり、適合工数を更に削減できる。また、物理パラメータで実験計画範囲の境界を設定することで、複数の制御パラメータを同時に変化させる場合でも、実験計画範囲の境界のラインが単調な傾向となり、少ない計測点数でも実験計画範囲の境界を精度良く予測することができる。
【００１４】
また、請求項４のように、複数の制御パラメータの実験計画範囲の境界を定める物理パラメータの判定閾値を当該物理パラメータが安定燃焼限界を越えないように決定する際に、いずれかの制御パラメータの可動範囲の限界でも前記物理パラメータが前記安定燃焼限界を越えない場合は、当該制御パラメータを可動範囲の限界に固定した状態で他の制御パラメータを前記物理パラメータが前記安定燃焼限界を越えるまで変化させて前記複数の制御パラメータの実験計画範囲を設定するようにしても良い。このようにすれば、単一の制御パラメータをその可動範囲の限界まで変化させても物理パラメータが安定燃焼限界を越えない場合に、他の制御パラメータを変化させて物理パラメータが安定燃焼限界を越える計測点を探索することが可能となり、実験計画範囲を安定燃焼限界まで拡大することができる。
【００１５】
また、請求項５のように、制御パラメータの適合値を決定する際にエンジン特性モデルと制御パラメータの実験計画範囲を考慮して当該制御パラメータの適合値を決定するようにしても良い。このようにすれば、例えば、エンジンの機差やばらつきを考慮して実験計画範囲の境界から所定の余裕度を持たせた範囲で制御パラメータの適合値を決定することが可能となり、適合値の精度やロバスト性を高めることができる。
【００１６】
また、請求項６のように、代表運転条件での計測結果に基づいて当該代表運転条件での制御パラメータの実験計画範囲を設定すると共に、運転条件の変化による物理パラメータの変化を予測し、前記代表運転条件での制御パラメータの実験計画範囲と前記運転条件の変化による物理パラメータの変化の予測結果に基づいて他の運転条件での制御パラメータの実験計画範囲を設定するようにしても良い。このようにすれば、代表運転条件でのみ物理パラメータを計測すれば良く、計測工数を削減できる。また、運転条件の変化による物理パラメータの変化を予測するため、危険性のある運転条件での安定燃焼限界を実機で探索する必要がなくなると共に、計測精度にとらわれることなく、他の運転条件での制御パラメータの実験計画範囲を設定することができる。
【００１７】
尚、請求項７は、請求項１に記載の適合方法の発明と実質的に同一の技術思想を適合装置の発明として記載したものである。
【図面の簡単な説明】
【００１８】
【図１】図１は本発明の実施例１における適合システムの構成例を概略的に示す図である。
【図２】図２は適合対象となる制御パラメータの種類を説明する図である。
【図３】図３は実施例１の実験計画範囲設定プログラムの処理の流れを示すフローチャートである。
【図４】図４は実施例１のＥＧＲ開度の実験計画範囲を筒内ＥＧＲ率を用いて設定する方法を説明する図である。
【図５】図５は実施例１のＶＣＴ進角値の実験計画範囲を筒内ＥＧＲ率を用いて設定する方法を説明する図である。
【図６】図６は実施例１の噴射時期遅角値の実験計画範囲を蒸発燃料量を用いて設定する方法を説明する図である。
【図７】図７は実施例１のＥＧＲ開度とＶＣＴ進角値を同時に変化させた場合の実験計画範囲を筒内ＥＧＲ率を用いて設定する方法を説明する図である。
【図８】図８は実施例２の実験計画範囲設定プログラムの処理の流れを示すフローチャートである。
【図９】図９は実施例２の物理パラメータ（例えば筒内ＥＧＲ率）と燃焼安定指標（例えばトルク変動）との関係を示す図である。
【図１０】図１０は実施例２の制御パラメータの実験計画範囲を物理パラメータが判定閾値以下となるように設定する方法を説明する図である。
【図１１】図１１は実施例３の実験計画範囲設定プログラムの処理の流れを示すフローチャートである。
【図１２】図１２は実施例３のＥＧＲ開度と筒内ＥＧＲ率との関係の計測データを示す図である。
【図１３】図１３は実施例３の噴射時期と噴霧移動距離との関係の計測データを示す図である。
【図１４】図１４は実施例３の筒内ＥＧＲ率と噴霧移動距離とをパラメータとする評価関数が同一値となるラインを示す図である。
【図１５】図１５は実施例３のＥＧＲ開度と噴射時期を同時に変化させる場合の実験計画範囲を示す図である。
【図１６】図１６は実施例４の実験計画範囲設定プログラムの処理の流れを示すフローチャートである。
【図１７】図１７は実施例４のＶＣＴ進角値と筒内ＥＧＲ率との関係の計測データを示す図である。
【図１８】図１８は実施例４の筒内ＥＧＲ率とトルク変動と安定燃焼限界との関係を説明する図である。
【図１９】図１９は実施例４のＶＣＴ進角値を可動限界に固定した状態で、ＥＧＲ開度をインクリメント量ずつ変化させて物理パラメータを算出又は計測する処理を説明する図である。
【図２０】図２０は実施例４の実験計画範囲の設定方法で設定した実験計画範囲を示す図である。
【図２１】図２１は実施例５の実験計画範囲と適合範囲との関係を説明する図である。
【図２２】図２２は実施例６の負荷変化と物理パラメータとの関係を説明する図である。
【図２３】図２３は実施例６の回転変化と物理パラメータとの関係を説明する図である。
【図２４】図２４（ａ）は代表運転条件の物理パラメータと評価関数１との関係を示す図、同図（ｂ）は代表運転条件の実験計画範囲を示す図、同図（ｃ）は負荷変化時の実験計画範囲を示す図、同図（ｄ）は回転変化時の実験計画範囲を示す図である。
【発明を実施するための形態】
【００１９】
以下、本発明を実施するための形態を具体化した６つの実施例１〜６を説明する。
【実施例１】
【００２０】
本発明の実施例１を図１乃至図７に基づいて説明する。
まず、適合システムの構成を図１に基づいて説明する。
適合するエンジン１１（内燃機関）をベンチ１２上に取り付け、このエンジン１１のクランク軸を動力計１３に連結する。適合作業中は、エンジン１１に装着された後述する各種のアクチュエータを電子制御ユニット（ＥＣＵ）１４によって制御する。この電子制御ユニット１４は、通信ボックス１６を介して適合制御用コンピュータ１７に接続され、適合作業中は、この適合制御用コンピュータ１７から通信ボックス１６を介して電子制御ユニット１４に適合制御信号を送信することで、電子制御ユニット１４内の各制御パラメータのマップ定数等を変更する。適合作業中のエンジン１１のスロットル開度は、スロットルコントロール装置１５によって調整される。
【００２１】
適合作業中は、動力計制御盤１８によって動力計１３とスロットルコントロール装置１５を制御してエンジン負荷を制御すると共に、動力計１３で計測したエンジントルクを適合制御用コンピュータ１７に送信する。エンジン１１には、クーラント（冷却水）の温度を調整するクーラント温度調整装置１９と、エンジンオイルの温度を調整するオイル温度調整装置２０と、燃料の温度を調整する燃費計測機能付きの燃料温度調整装置２１を接続し、適合作業中は、これら各温度調整装置２１によってクーラント温度、エンジンオイル温度、燃料温度が一定条件に自動調整される。適合作業中にエンジン１１から排出される排出ガスは、排出ガス分析計２２で分析され、排出ガス中のＮＯｘ、ＣＯ、ＨＣ等のエミッションの測定結果が適合制御用コンピュータ１７に送信される。
【００２２】
この適合システムで適合可能なエンジン１１は、吸気ポート噴射エンジン、直噴エンジン、ディーゼルエンジン等のいずれの方式でも良い。適合対象となる制御パラメータは、例えば、燃料噴射弁２３の噴射時期、点火プラグ２４の点火時期、モータ等で駆動されるスロットルバルブ２５の開度（スロットル開度）、スワールコントロールバルブ２６の開度（ＳＣＶ開度）、吸気バルブ及び排気バルブの可変バルブタイミング機構２７の進角値（ＶＣＴ進角値）、排気環流制御バルブ２８の開度（ＥＧＲ開度）等である。
【００２３】
制御パラメータを適合する場合は、後述する方法で決定した制御パラメータの実験計画範囲内に所定数の計測点を配置して、当該制御パラメータを各計測点に変化させて、各計測点で当該制御パラメータの変化に応じて変化する物理量（以下「物理パラメータ」という）を計測し、その計測結果に基づいて前記制御パラメータと前記物理パラメータとの関係をモデル化してエンジン特性モデルを作成し、前記エンジン特性モデルを用いて制御パラメータの適合値をエンジン１１の要求性能が満たされるように決定する。
【００２４】
従来の適合方法では、制御パラメータの実験計画範囲（計測点配置範囲）を当該制御パラメータで設定するようにしていたが、制御パラメータは、エンジン１１を動作させる際の指令値であり、実際の燃焼との関係が複雑である。しかも、燃焼の結果によって得られる物理パラメータの計測値（例えば図示平均有効圧ＩＭＥＰの標準偏差、トルク変動等）は、ばらつきが大きく、制御パラメータとの関係を正確に定義するには多くの計測工数が必要となる。
【００２５】
そこで、本実施例１では、制御パラメータの実験計画範囲を燃焼に寄与する物理パラメータで設定するようにしている。燃焼に寄与する物理パラメータは、燃焼状態を表す指標となるため、燃焼に寄与する物理パラメータを用いて実験計画範囲を設定すれば、比較的少ない計測工数で、実験計画範囲を安定燃焼領域内に精度良く設定することが可能となる。また、制御パラメータは、それにより変化する物理パラメータが共通のものがあるため、変化する物理パラメータが共通の複数の制御パラメータを同時に変化させた場合、制御パラメータの実験計画範囲も複合的に変化することがあるが、実験計画範囲を物理パラメータにて設定すれば、比較的少ない計測工数で、複合的な実験計画範囲も安定燃焼領域内に精度良く設定することが可能となる。これにより、エンジン特性モデルの精度向上（制御パラメータの適合値の精度向上）と工数削減とを両立させることができる。
【００２６】
ここで、制御パラメータがＥＧＲ開度の場合は、それにより変化する物理パラメータは筒内ＥＧＲ率、吸気管圧力等である。
制御パラメータがＶＣＴ進角値の場合は、それにより変化する物理パラメータは、筒内ＥＧＲ率、筒内流速、吸気温度、ポンピングロス、吸気管圧力、実圧縮比等である。
制御パラメータが噴射時期の場合は、それにより変化する物理パラメータは、噴霧移動距離、霧化時間、蒸発燃料量（割合）、噴射時筒内流速等である。
【００２７】
本実施例１では、制御パラメータと物理パラメータとの関係をシミュレーション又は実機で求め、物理パラメータで実験計画範囲を設定する。例えば、筒内ＥＧＲ率が所定％以下又は筒内流速が所定ｍ／ｓ以上等である。
【００２８】
或は、複数の物理パラメータのＡＮＤ（論理積）をとっても良い。例えば、筒内ＥＧＲ率が所定％以下且つ筒内流速が所定ｍ／ｓ以上となる範囲を、ＶＣＴ進角値の実験計画範囲としても良い。
【００２９】
また、ＥＧＲと可変バルブタイミング機構２７（ＶＣＴ）は、排気ガスを筒内に戻して内部ＥＧＲ率を変化させる役割を持つため、この２つの制御パラメータは、相互相関し、ＥＧＲによる外部ＥＧＲ率とＶＣＴによる内部ＥＧＲ率とを合計した筒内ＥＧＲ率により、設定すべき実験計画範囲も変化する。そこで、シミュレーションによりＥＧＲとＶＣＴを変化させた場合に筒内に導入されるＥＧＲ量を算出し、同一の筒内ＥＧＲ率以下となる範囲を実験計画範囲とする。
【００３０】
尚、実験計画法で所定数の計測点を予め広い範囲に配置し、物理パラメータが判定閾値を越える計測点を、計測前に排除すると共に、排除した計測点を物理パラメータが判定閾値以下となる実験計画範囲内に再配置するようにしても良い。
【００３１】
また、設定した実験計画範囲内に配置したいずれかの計測点で計測した物理パラメータが判定閾値を越える場合は、当該実験計画範囲を修正して、物理パラメータが判定閾値を越えた計測点を修正後の実験計画範囲内に再配置するようにしても良い。
実験計画法で物理パラメータを配置する際に、物理パラメータの特性を考慮して物理パラメータを配置するようにしても良い。
【００３２】
本実施例１の実験計画範囲の設定は、図３の実験計画範囲設定プログラム（実験計画範囲設定手段）に従って次の手順で行われる。
まず、ステップ１０１で、適合対象となる制御パラメータに対する物理パラメータを選択する。
【００３３】
例えば、適合対象となる制御パラメータがＥＧＲ開度であれば、それに対する物理パラメータとして、筒内ＥＧＲ率と吸気管圧力のいずれかを選択する（図４の例では筒内ＥＧＲ率を選択している）。
【００３４】
適合対象となる制御パラメータがＶＣＴ進角値であれば、それに対する物理パラメータとして、筒内ＥＧＲ率、筒内流速、吸気温度、ポンピングロス、吸気管圧力、実圧縮比の中から選択する（図５の例では筒内ＥＧＲ率を選択している）。
【００３５】
適合対象となる制御パラメータが噴射時期であれば、それに対する物理パラメータとして、噴霧移動距離、霧化時間、蒸発燃料量（割合）、噴射時筒内流速の中から選択する（図６の例では蒸発燃料量を選択している）。
【００３６】
図７に示すように、ＥＧＲ開度とＶＣＴ進角値を同時に変化させる場合は、両者に共通する物理パラメータとして、筒内ＥＧＲ率と吸気管圧力のいずれかを選択すれば良い（図７の例では筒内ＥＧＲ率を選択している）。
【００３７】
次のステップ１０２で、制御パラメータと物理パラメータとの関係（図４〜図７の実線）をシミュレーション又は実機で求めた計測データにより算出する。この後、ステップ１０３に進み、制御パラメータの実験計画範囲の境界を定める物理パラメータの判定閾値を生じさせる制御パラメータの値を算出して、次のステップ１０４で、制御パラメータの実験計画範囲を決定する。
【００３８】
以上のようにして決定した制御パラメータの実験計画範囲内に所定数の計測点を配置して、当該制御パラメータを各計測点に変化させて、各計測点で物理パラメータを計測し、その計測結果に基づいて制御パラメータと物理パラメータとの関係をモデル化してエンジン特性モデルを作成し、このエンジン特性モデルを用いて制御パラメータの適合値をエンジン１１の要求性能が満たされるように決定する。
【００３９】
以上説明した本実施例１によれば、制御パラメータの実験計画範囲を燃焼状態を表す指標となる物理パラメータで設定するようにしたので、比較的少ない計測工数で、実験計画範囲を安定燃焼領域内に精度良く設定することが可能となる。また、制御パラメータは、変化する物理パラメータが共通のものがあるため、変化する物理パラメータが共通の複数の制御パラメータを同時に変化させた場合、制御パラメータの実験計画範囲も複合的に変化することがあるが、複数の制御パラメータに対して共通の物理パラメータを用いて実験計画範囲を設定すれば、比較的少ない計測工数で、複合的な実験計画範囲も安定燃焼領域内に精度良く設定することが可能となる。これにより、エンジン特性モデルの精度向上（制御パラメータの適合値の精度向上）と工数削減とを両立させることができる。
【実施例２】
【００４０】
次に、図８乃至図１０を用いて本発明の実施例２を説明する。
本実施例２では、制御パラメータの実験計画範囲の境界を定める物理パラメータの判定閾値を当該物理パラメータが安定燃焼限界（安定燃焼領域の限界）を越えないように決定する。ここで、安定燃焼限界は、燃焼安定指標によって判断する。燃焼安定指標としては、例えばＣＯＶ（図示平均有効圧ＩＭＥＰの標準偏差）、トルク変動、点火プラグ近傍当量比等を用いれば良い。これらの燃焼安定指標を用いることで、安定燃焼限界を精度良く判定することができる。
【００４１】
燃焼悪化領域では、失火や燃焼不安定になったり、異常燃焼によってエンジンが損傷する可能性があるため、計測点が安定燃焼限界を越えて燃焼悪化領域に配置されると、当該計測点の計測精度が悪化したり、計測不能となる。
【００４２】
安定燃焼限界を定める物理パラメータの判定閾値は、エンジン機種や製品の仕様により変化することがあるため、適合対象システム毎に物理パラメータと燃焼安定指標との関係を求めて、この関係から物理パラメータの判定閾値を求めるようにしても良い。
【００４３】
適合対象システムの制御パラメータと燃焼安定指標との関係は、実機又はシミュレーションで求めれば良い。この関係を実機で求める場合は、図９に示すように、制御パラメータを変化させたときの物理パラメータ（例えば筒内ＥＧＲ率）と燃焼安定指標（例えばトルク変動）を計測し、計測した２つのデータの関係を求め、燃焼安定指標が安定燃焼限界を越える物理パラメータの値を判定閾値として決定すれば良い。
【００４４】
判定閾値の決定後、図１０に示すように、制御パラメータ（例えばＥＧＲ開度）と物理パラメータ（例えば筒内ＥＧＲ率）との関係において、物理パラメータが図９の判定閾値以下となる範囲を実験計画範囲とすれば良い。これにより、実験計画範囲を安定燃焼領域内に設定することができる。
【００４５】
本実施例２の実験計画範囲の設定は、図８の実験計画範囲設定プログラムに従って次の手順で行われる。まず、ステップ２０１で、前記実施例１と同様の方法で、適合対象となる制御パラメータに対する物理パラメータを選択し、次のステップ２０２で、選択した物理パラメータに対する燃焼安定指標を、ＣＯＶ、トルク変動、点火プラグ近傍当量比等の中から選択する。
【００４６】
次のステップ２０３で、制御パラメータに対する物理パラメータと燃焼安定指標をシミュレーション又は実機により算出又は計測する。この後、ステップ２０４で、物理パラメータと燃焼安定指標との関係から、燃焼安定指標が安定燃焼限界を越える物理パラメータの値を判定閾値として決定する。
【００４７】
この後、ステップ２０５で、制御パラメータと物理パラメータとの関係から、物理パラメータが判定閾値となる制御パラメータの値を実験計画範囲の境界値として算出し、次のステップ２０６で、制御パラメータの実験計画範囲を決定する。
【００４８】
以上説明した本実施例２によれば、制御パラメータの実験計画範囲の境界を定める物理パラメータの判定閾値を当該物理パラメータが安定燃焼限界（安定燃焼領域の限界）を越えないように決定することができるため、計測点が安定燃焼限界を越えた領域に配置されることを確実に防止でき、全ての計測点で物理パラメータを精度良く計測することができる。
【実施例３】
【００４９】
次に、図１１乃至図１５を用いて本発明の実施例３を説明する。
本実施例３では、単一の制御パラメータを変化させて当該単一の制御パラメータと物理パラメータとの関係を実機又はシミュレーションにより求め、複数の制御パラメータを同時に変化させる場合の実験計画範囲を前記単一の制御パラメータと物理パラメータとの関係を用いて設定する。このようにすれば、実機又はシミュレーションにより全ての制御パラメータの組み合わせを計測する必要がなくなり、適合工数を更に削減できる。
【００５０】
また、物理パラメータで実験計画範囲の境界を設定することで、複数の制御パラメータを同時に変化させる場合でも、実験計画範囲の境界のラインが単調な傾向となり、少ない計測点数でも実験計画範囲の境界を精度良く予測することができる。
【００５１】
例えば、ＥＧＲ開度と噴射時期を同時に変化させる場合は、まず、図１２に示すように、ＥＧＲ開度と筒内ＥＧＲ率との関係を実機又はシミュレーションにより求めると共に、図１３に示すように、噴射時期と噴霧移動距離（燃料噴射弁２３の噴射燃料が燃焼するまでに移動する距離）との関係を実機又はシミュレーションにより計測又は算出する。
【００５２】
計測又は算出した２つの物理パラメータに相関関係がある場合は、重み付け、加減乗算等により、１つの評価関数を作成し、相関関係が無い場合には、独立した評価関数を作成する。
【００５３】
図１２、図１３の例では、計測又は算出した２つの物理パラメータは、筒内ＥＧＲ率と噴霧移動距離であり、相関関係があるため、２つの物理パラメータを変数とする１つの評価関数を作成する。
評価関数＝ａ×筒内ＥＧＲ率＋ｂ×噴霧移動距離
上式において、ａ，ｂは係数である。この評価関数は、大きな値になるほど、燃焼安定性が低下することを意味する。
【００５４】
上式のように、評価関数を、２つの物理パラメータ（筒内ＥＧＲ率と噴霧移動距離）を変数とする一次式で定義する場合は、図１４に示すように、評価関数が同一の値になるラインが直線となる。評価関数の値は、燃焼安定指標となるため、評価関数の値が所定の閾値以下になる範囲を安定燃焼領域と判断し、評価関数の値が所定の閾値（安定燃焼限界）となる物理パラメータ（筒内ＥＧＲ率と噴霧移動距離）を用いて、評価関数の値が所定の閾値となる制御パラメータ（ＥＧＲ開度と噴射時期）を、図１２と図１３の関係から算出し、図１５に示すように、評価関数の値が所定の閾値以下となる制御パラメータ（ＥＧＲ開度と噴射時期）の実験計画範囲を設定する。これにより、実験計画範囲を安定燃焼領域内に設定することができる。
【００５５】
本実施例３の実験計画範囲の設定は、図１１の実験計画範囲設定プログラムに従って次の手順で行われる。まず、ステップ３０１で、前記実施例１と同様の方法で、適合対象となる制御パラメータに対する物理パラメータを選択し、次のステップ３０２で、単一の制御パラメータと物理パラメータとの関係をシミュレーション又は実機により算出又は計測する処理を、適合する制御パラメータの数だけ繰り返す。
【００５６】
この後、ステップ３０３で、相関関係がある物理パラメータを用いて評価関数の値を算出する。次のステップ３０４で、評価関数の値が所定の閾値（安定燃焼限界）となる制御パラメータの値を実験計画範囲の境界値として算出し、次のステップ３０５で、制御パラメータの実験計画範囲を決定する。
【実施例４】
【００５７】
次に、図１６乃至図２０を用いて本発明の実施例４を説明する。
本実施例４では、複数の制御パラメータの実験計画範囲の境界を定める物理パラメータの判定閾値を当該物理パラメータが安定燃焼限界を越えないように決定する際に、いずれかの制御パラメータの可動範囲の限界（可動限界）でも前記物理パラメータが前記安定燃焼限界を越えない場合は、当該制御パラメータを可動限界に固定した状態で他の制御パラメータを前記物理パラメータが前記安定燃焼限界を越えるまで変化させて前記複数の制御パラメータの実験計画範囲を設定する。このようにすれば、単一の制御パラメータをその可動限界まで変化させても物理パラメータが安定燃焼限界を越えない場合に、他の制御パラメータを変化させて物理パラメータが安定燃焼限界を越える計測点を探索することができ、実験計画範囲を安定燃焼限界まで拡大することができる。また、単独で可動範囲内を変化させても燃焼安定性に大きな影響を与えない制御パラメータを考慮した実験計画範囲を精度良く設定できる。
【００５８】
例えば、ＶＣＴ進角値とＥＧＲ開度の実験計画範囲を設定する方法を図１７乃至図２０を用いて説明する。図１７及び図１８に示すように、ＶＣＴ進角値を可動限界まで変化させても、物理パラメータである筒内ＥＧＲ率は、安定燃焼限界に相当する判定閾値を越えない。
【００５９】
そこで、図１９に示すように、ＶＣＴ進角値を可動限界に固定した状態で、他の制御パラメータであるＥＧＲ開度を所定のインクリメント量ずつ変化させ（以下この処理を「インクリメント」という）、物理パラメータである筒内ＥＧＲ率をシミュレーション又は実機により算出又は計測する処理を、筒内ＥＧＲ率が安定燃焼限界に相当する判定閾値を越えるまで繰り返す。尚、ＶＣＴ進角値を変化させて求める物理パラメータと、筒内ＥＧＲ率を変化させて求める物理パラメータは、異なっても良い。
【００６０】
ＥＧＲ開度のインクリメント量の算出方法は、安定燃焼限界（判定閾値）に相当する物理パラメータの値を整数で割り算した値をＥＧＲ開度の変化量に換算して得られた値をインクリメント量とすれば良い。
【００６１】
本実施例４の実験計画範囲の設定は、図１６の実験計画範囲設定プログラムに従って次の手順で行われる。まず、ステップ４０１で、前記実施例１と同様の方法で、適合対象となる制御パラメータＡに対する物理パラメータを選択し、次のステップ４０２で、選択した物理パラメータに対する燃焼安定指標を、ＣＯＶ、トルク変動、点火プラグ近傍当量比等の中から選択する。
【００６２】
次のステップ４０３で、制御パラメータＡに対する物理パラメータと燃焼安定指標をシミュレーション又は実機により算出又は計測する。この後、ステップ４０４で、物理パラメータと燃焼安定指標との関係から、燃焼安定指標が安定燃焼限界に達する物理パラメータの値を判定閾値として決定する。
【００６３】
この後、ステップ４０５で、物理パラメータが判定閾値を越えたか否かを判定し、物理パラメータが判定閾値を越えていなければ、ステップ４０９に進み、制御パラメータＡと相関の強い制御パラメータＢを選択する。次のステップ４１０で、制御パラメータＢに対する物理パラメータと燃焼安定指標をシミュレーション又は実機により算出又は計測する。この後、ステップ４１１で、物理パラメータと燃焼安定指標との関係から、燃焼安定指標が安定燃焼限界に達する物理パラメータの値を判定閾値として決定する。
【００６４】
この後、ステップ４１２で、安定燃焼限界（判定閾値）に相当する物理パラメータの値を整数で割り算した値を制御パラメータＢの変化量に換算して得られた値をインクリメント量とする。次のステップ４１３で、制御パラメータＡを可動限界に固定した状態で、制御パラメータＢをインクリメント量だけ変化させて物理パラメータを算出又は計測する。      
【００６５】
この後、ステップ４１４で、物理パラメータが判定閾値を越えたか否かを判定し、物理パラメータが判定閾値を越えていなければ、ステップ４１３に戻る。これにより、物理パラメータが判定閾値を越えるまで、制御パラメータＢをインクリメント量ずつ変化させて物理パラメータを算出又は計測する処理を繰り返す。
【００６６】
その後、上記ステップ４１４又は４０５で、物理パラメータが判定閾値を越えたと判定された時点で、ステップ４０６に進み、制御パラメータと物理パラメータとの関係から、物理パラメータが判定閾値となる制御パラメータの値を実験計画範囲の境界値として算出し、次のステップ４０７で、制御パラメータの実験計画範囲を決定する（図２０参照）。
【実施例５】
【００６７】
次に、図２１を用いて本発明の実施例５を説明する。
本実施例５では、制御パラメータの適合値を決定する際に、エンジン特性モデルと制御パラメータの実験計画範囲を考慮して当該制御パラメータの適合値を決定する。このようにすれば、例えば、製品機差やばらつきを考慮して実験計画範囲の境界から所定の余裕度を持たせた適合範囲で制御パラメータの適合値を決定することが可能となり、適合値の精度やロバスト性を高めることができる。
【実施例６】
【００６８】
次に、図２２乃至図２４を用いて本発明の実施例５を説明する。
本実施例６では、代表運転条件での計測結果に基づいて当該代表運転条件での制御パラメータの実験計画範囲を設定すると共に、運転条件の変化による物理パラメータの変化を予測し、前記代表運転条件での制御パラメータの実験計画範囲と前記運転条件の変化による物理パラメータの変化の予測結果に基づいて他の運転条件での制御パラメータの実験計画範囲を設定する。このようにすれば、代表運転条件でのみ物理パラメータを計測すれば良く、計測工数を削減できる。また、運転条件の変化による物理パラメータの変化を予測するため、危険性のある運転条件での安定燃焼限界を実機で探索する必要がなくなると共に、計測精度にとらわれることなく、他の運転条件での制御パラメータの実験計画範囲を設定することができる。
【００６９】
具体的には、まず、代表運転条件で制御パラメータと物理パラメータとの関係をシミュレーション又は実機により算出又は計測して、例えば、前記実施例３と同様の方法で、代表運転条件での制御パラメータの実験計画範囲［図２４（ｂ）の斜線部分］を設定する。図２４の例では、制御パラメータの実験計画範囲は、評価関数１と評価関数２とで囲まれた領域となるが、代表運転条件では評価関数２は全運転領域を許容するようになっているため、代表運転条件での制御パラメータの実験計画範囲は、評価関数１で囲まれた領域となる。
【００７０】
次に、図２２に示すように、負荷が変化したときの制御パラメータと物理パラメータとの関係を予測して、負荷変化時の制御パラメータの実験計画範囲［図２４（ｃ）の斜線部分］を予測する。負荷変化時には、評価関数２は全運転領域を許容するため、負荷変化時の制御パラメータの実験計画範囲は、評価関数１で囲まれた領域となる。
【００７１】
更に、図２３に示すように、エンジン回転速度が変化したときの制御パラメータと物理パラメータとの関係を予測して、図２４（ｄ）に示す回転変化時の制御パラメータの実験計画範囲［図２４（ｄ）の斜線部分］を予測する。回転変化時の制御パラメータの実験計画範囲は、評価関数１と評価関数２とで囲まれた領域となる。
【００７２】
この後、代表運転条件での制御パラメータの実験計画範囲（評価関数１）と運転条件の変化による物理パラメータの変化の予測結果に基づいて他の運転条件での制御パラメータの実験計画範囲を設定する。
【００７３】
本発明は、上記各実施例１〜６に限定されず、適合対象となる制御パラメータと、当該制御パラメータの変化に応じて変化する物理パラメータを適宜変更しても良い等、要旨を逸脱しない範囲内で種々変更して実施できる。
【符号の説明】
【００７４】
１１…エンジン、１４…電子制御ユニット（ＥＣＵ）、１７…適合制御用コンピュータ（実験計画範囲設定手段）
【図１】
2012013637.tif 000003
【図２】
2012013637.tif 000004
【図３】
2012013637.tif 000005
【図４】
2012013637.tif 000006
【図５】
2012013637.tif 000007
【図６】
2012013637.tif 000008
【図７】
2012013637.tif 000009
【図８】
2012013637.tif 000010
【図９】
2012013637.tif 000011
【図１０】
2012013637.tif 000012
【図１１】
2012013637.tif 000013
【図１２】
2012013637.tif 000014
【図１３】
2012013637.tif 000015
【図１４】
2012013637.tif 000016
【図１５】
2012013637.tif 000017
【図１６】
2012013637.tif 000018
【図１７】
2012013637.tif 000019
【図１８】
2012013637.tif 000020
【図１９】
2012013637.tif 000021
【図２０】
2012013637.tif 000022
【図２１】
2012013637.tif 000023
【図２２】
2012013637.tif 000024
【図２３】
2012013637.tif 000025
【図２４】
2012013637.tif 000026