(57)【特許請求の範囲】
【請求項１】
模擬の対象である空間内の流体物理量を一様なものとして模擬する集中定数系流体モデルと、前記流体物理量を流体の進行方向の位置に応じた量として模擬する１次元CFDモデルとを用いる内燃機関の性能予測解析方法であって、
機関回転数、スロットル開度、及び負荷率の少なくとも１つに基づいて脈動条件の成否を判断するステップと、
前記脈動条件の成立時に、前記１次元CFDモデルを用いて前記流体物理量を模擬するステップと、
前記脈動条件の不成立時に、前記集中定数系流体モデルを用いて前記流体物理量を模擬するステップと、
前記流体物理量の模擬値に基づいて筒内流入空気量を算出するステップとを含み、
前記脈動条件は、前記筒内流入空気量に、前記空間内での脈動の影響が及ぶ条件であることを特徴とする内燃機関の性能予測解析方法。
【請求項２】
前記空間は、内燃機関の吸気通路であり、
前記集中定数系流体モデル、及び前記１次元CFDモデルは、何れも吸気通路を流れる空気の流体物理量を模擬するモデルであることを特徴とする請求項１記載の内燃機関の性能予測解析方法。
【請求項３】
機関回転数が、脈動下限回転数と脈動上限回転数との間にある場合に、前記脈動条件の成立を判定することを特徴とする請求項２記載の内燃機関の性能予測解析方法。
【請求項４】
模擬の対象である空間内の流体物理量を一様なものとして模擬する集中定数系流体モデルと、前記流体物理量を流体の進行方向の位置に応じた量として模擬する１次元CFDモデルとを用いて内燃機関の性能予測解析を行うシステムであって、
機関回転数、スロットル開度、及び負荷率の少なくとも１つに基づいて脈動条件の成否を判断する脈動条件判定手段と、
前記脈動条件の成立時に、前記１次元CFDモデルを用いて前記流体物理量を模擬する１次元模擬手段と、
前記脈動条件の不成立時に、前記集中定数系流体モデルを用いて前記流体物理量を模擬する０次元模擬手段と、
前記流体物理量の模擬値に基づいて筒内流入空気量を算出する空気量算出手段とを含み、
前記脈動条件は、前記筒内流入空気量に、前記空間内での脈動の影響が及ぶ条件であることを特徴とする内燃機関の性能予測解析システム。
【請求項５】
前記空間は、内燃機関の吸気通路であり、
前記集中定数系流体モデル、及び前記１次元CFDモデルは、何れも吸気通路を流れる空気の流体物理量を模擬するモデルであることを特徴とする請求項４記載の内燃機関の性能予測解析システム。
【請求項６】
前記脈動条件判定手段は、機関回転数が、脈動下限回転数と脈動上限回転数との間にある場合に、前記脈動条件の成立を判定することを特徴とする請求項５記載の内燃機関の性能予測解析システム。
【請求項７】
前記集中定数系流体モデル、及び前記１次元CFDモデルは、何れもアクセル開度をパラメータとして含むモデルであり、
前記筒内流入空気量に基づいて内燃機関の発生するトルクを模擬するトルク模擬手段と、
前記トルクの模擬値に基づいて、機関回転数及び車速を模擬する車両模擬手段と、
前記機関回転数の模擬値を前記脈動条件判定手段に提供する機関回転数提供手段と、
規定の走行モードを記憶する走行モード記憶手段と、
前記走行モードによる目標車速と前記車速の模擬値とに基づいて、前記目標車速を実現するためのアクセル開度を算出するアクセル開度算出手段と、
前記アクセル開度の算出値を、前記集中定数系流体モデル及び前記１次元CFDモデルに提供するアクセル開度提供手段と、
を含むことを特徴とする請求項６記載の内燃機関の性能予測解析システム。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、内燃機関の性能予測解析方法及び性能予測解析システムに係り、特に、筒内流入空気量の模擬を前提とする性能予測解析に好適な内燃機関の性能予測解析方法及び性能予測解析システムに関する。
【背景技術】
【０００２】
従来、特開２００５−２４９４２０号公報に開示されているように、CFDを用いた内燃機関の性能予測手法が知られている。CFDは、作動流体の運動を模擬する物理モデルを用いたシミュレーションの手法である。上記従来の性能予測方法によれば、内燃機関の諸元を特定することで、CFDにより、吸気や排気の流れを模擬することができ、その結果から、上記の諸元で特定される内燃機関の性能を予測することができる。
【０００３】
上記従来の方法では、吸気や排気の流れを１次元の流れとして模擬する１次元CFDモデルと、その流れを３次元の流れとして模擬する３次元CFDモデルとが適宜選択的に用いられる。１次元CFDモデルを用いた演算（以後、「１次元CFD演算」と称す）は、３次元CFDモデルを用いた演算（以後、「３次元CFD演算」と称す）に比して負荷が軽く、短時間で行うことができる。他方、作動流体の挙動を高い精度で模擬するうえでは、３次元CFD演算が１次元CFD演算に比して優れている。
【０００４】
上記従来の公報には、内燃機関の性能を予測するにあたり、基本的には１次元CFDを行い、特に精度の要求される部位に限って、３次元CFDを行うことが記載されている。１次元CFDと３次元CFDとを、このように切り換えて用いることによれば、演算負荷が不当に重くなるのを避けつつ、十分な精度で内燃機関の性能を予測することができる。このため、上記従来の手法は、諸元の特定された内燃機関の性能を、短時間で正確に予測する手法として有用である。
【０００５】
【特許文献１】特開２００５−２４９４２０号公報
【特許文献２】特開２００５−２４９４１８号公報
【特許文献３】特開２００５−２４９４１９号公報
【特許文献４】特開２００５−２４９４２２号公報
【発明の開示】
【発明が解決しようとする課題】
【０００６】
内燃機関の吸気や排気の流れには、内燃機関の運転状態に応じて脈動が生ずることがある。このような脈動の影響をも考慮して吸気や排気の流れを模擬するためには、少なくとも、その流れを、吸気や排気の進行方向において、個々の位置毎に模擬することが必要である。そして、このような模擬を実現するためには、１次元以上のCFDモデルを用いることが必要である。上述した従来の方法では、常に１次元以上のCFDモデルが用いられる。このため、この方法によれば、吸気や排気の挙動を、常に脈動を考慮して、精度良く模擬することが可能である。
【０００７】
しかしながら、内燃機関における脈動は、常に出力トルクに影響を与えるものではない。例えば、機関回転数が十分に低い低回転領域や、その回転数が十分に高い高回転領域では、内燃機関のトルクに、吸気の脈動の影響は殆ど及ばない。そして、このような状況下では、吸気や排気の挙動が、進行方向に対して一様であるものとして取り扱っても、内燃機関の性能を正しく評価することができる。つまり、このような状況下では、脈動の影響を無視してモデル演算を行っても、内燃機関の動作を正確に模擬することが可能である。この点、常に１次元以上のCFD演算を行う従来の方法は、演算負荷の軽減を図る意味で、必ずしも最適化されたものではなかった。
【０００８】
この発明は、上述のような課題を解決するためになされたもので、脈動の影響が出力トルクに及ばない条件下では、その影響を無視したモデルで模擬を行うことにより、十分な演算負荷の軽減を可能とする内燃機関の性能予測解析方法を提供することを目的とする。
また、この発明は、上記の性能予測解析方法を自動的に実行する内燃機関の性能予測解析システムを提供することを第２の目的とする。
【課題を解決するための手段】
【０００９】
第１の発明は、上記の目的を達成するため、模擬の対象である空間内の流体物理量を一様なものとして模擬する集中定数系流体モデルと、前記流体物理量を流体の進行方向の位置に応じた量として模擬する１次元CFDモデルとを用いる内燃機関の性能予測解析方法であって、
機関回転数、スロットル開度、及び負荷率の少なくとも１つに基づいて脈動条件の成否を判断するステップと、
前記脈動条件の成立時に、前記１次元CFDモデルを用いて前記流体物理量を模擬するステップと、
前記脈動条件の不成立時に、前記集中定数系流体モデルを用いて前記流体物理量を模擬するステップと、
前記流体物理量の模擬値に基づいて筒内流入空気量を算出するステップとを含み、
前記脈動条件は、前記筒内流入空気量に、前記空間内での脈動の影響が及ぶ条件であることを特徴とする。
【００１０】
また、第２の発明は、第１の発明において、前記空間は、内燃機関の吸気通路であり、
前記集中定数系流体モデル、及び前記１次元CFDモデルは、何れも吸気通路を流れる空気の流体物理量を模擬するモデルであることを特徴とする。
【００１１】
また、第３の発明は、第２の発明において、機関回転数が、脈動下限回転数と脈動上限回転数との間にある場合に、前記脈動条件の成立を判定することを特徴とする。
【００１２】
また、第４の発明は、模擬の対象である空間内の流体物理量を一様なものとして模擬する集中定数系流体モデルと、前記流体物理量を流体の進行方向の位置に応じた量として模擬する１次元CFDモデルとを用いて内燃機関の性能予測解析を行うシステムであって、
機関回転数、スロットル開度、及び負荷率の少なくとも１つに基づいて脈動条件の成否を判断する脈動条件判定手段と、
前記脈動条件の成立時に、前記１次元CFDモデルを用いて前記流体物理量を模擬する１次元模擬手段と、
前記脈動条件の不成立時に、前記集中定数系流体モデルを用いて前記流体物理量を模擬する０次元模擬手段と、
前記流体物理量の模擬値に基づいて筒内流入空気量を算出する空気量算出手段とを含み、
前記脈動条件は、前記筒内流入空気量に、前記空間内での脈動の影響が及ぶ条件であることを特徴とする。
【００１３】
また、第５の発明は、第４の発明において、前記空間は、内燃機関の吸気通路であり、
前記集中定数系流体モデル、及び前記１次元CFDモデルは、何れも吸気通路を流れる空気の流体物理量を模擬するモデルであることを特徴とする。
【００１４】
また、第６の発明は、第５の発明において、前記脈動条件判定手段は、機関回転数が、脈動下限回転数と脈動上限回転数との間にある場合に、前記脈動条件の成立を判定することを特徴とする。
【００１５】
また、第７の発明は、第６の発明において、前記集中定数系流体モデル、及び前記１次元CFDモデルは、何れもアクセル開度をパラメータとして含むモデルであり、
前記筒内流入空気量に基づいて内燃機関の発生するトルクを模擬するトルク模擬手段と、
前記トルクの模擬値に基づいて、機関回転数及び車速を模擬する車両模擬手段と、
前記機関回転数の模擬値を前記脈動条件判定手段に提供する機関回転数提供手段と、
規定の走行モードを記憶する走行モード記憶手段と、
前記走行モードによる目標車速と前記車速の模擬値とに基づいて、前記目標車速を実現するためのアクセル開度を算出するアクセル開度算出手段と、
前記アクセル開度の算出値を、前記集中定数系流体モデル及び前記１次元CFDモデルに提供するアクセル開度提供手段と、
を含むことを特徴とする。
【発明の効果】
【００１６】
第１又は第４の発明によれば、脈動条件の成立時に、１次元CFDモデルを用いて流体物理量を模擬し、脈動条件の不成立時に、集中定数系流体モデルを用いて流体物理量を模擬することができる。集中定数系流体モデルの演算は、模擬の対象である空間内の流体物理量を一様なものとして模擬するため、１次元CFDモデルの演算に比して負荷が軽い。このため、本発明によれば、内燃機関の性能予測解析精度を下げることなく、模擬演算の負荷を軽くすることができる。
【００１７】
第２又は第５の発明によれば、集中定数系流体モデル又は１次元CFDモデルによって、内燃機関の吸気通路を流れる空気の流体物理量を、軽い演算負荷で、かつ正確に模擬し、その結果として、筒内流入空気量を、精度良く模擬することができる。
【００１８】
第３又は第６の発明によれば、機関回転数が、脈動下限回転数と脈動上限回転数との間にあるか否かを見ることで、脈動条件の成否を簡単に判定することができる。
【００１９】
第７の発明によれば、集中定数系流体モデル、及び１次元CFDモデルは、それぞれ、アクセル開度をパラメータとして吸気通路内の空気の流体物理量を模擬することができる。また、本発明によれば、それらの流体物理量に基づいて、内燃機関の発生するトルク、更には、機関回転数及び車速を模擬することができる。そして、脈動条件の成否は、機関回転数の模擬値に基づいて判断することができる。更に、本発明によれば、集中定数系流体モデル及び１次元CFDモデルに、車速の模擬値に基づいて算出されたアクセル開度を提供することができる。その結果、本発明によれば、車両が規定の走行モードで走行する際の内燃機関の状態を自動的に模擬することができる。
【発明を実施するための最良の形態】
【００２０】
実施の形態１．
［実施の形態１の構成］
図１は、本発明の実施の形態１の設計値最適化システムの構成を説明するためのブロック図である。本実施形態のシステムは、通常のコンピュータシステムに、後述する機能を実現するためのソフトウェアをインストールすることにより実現することができる。図１に示す各ブロックは、コンピュータシステムが、上記のソフトウェアに従った処理を行うことで実現される機能ブロックを、模式的に表したものである。
【００２１】
図１に示すように、本実施形態のシステムは、入力インターフェイス（IF)１０、及び出力IF１２を備えている。入力IF１０は、具体的には、キーボードやマウスにより実現される。また、出力IF１２は、ディスプレイやスピーカーにより実現される。
【００２２】
また、このシステムは、空気モデル１４を備えている。空気モデル１４は、内燃機関の設計値や、その運転状態を表す各種のパラメータに基づいて、筒内流入空気量を模擬するためのモデルである。本実施形態において、空気モデル１４には、集中定数系の空気モデル１６（以下、「集中定数系流体モデル」と称す）と、１次元CFDによる空気モデル１８（以下、「1DCFDモデル１８」と称す）とが含まれている。
【００２３】
集中定数系流体モデル１６は、模擬の対象である空間において流体が一様であるものとして、その流体の物理量を模擬するモデルである。集中定数系流体モデル１６によれば、内燃機関の吸気通路の設計値（吸気管長、吸気管径等）や、内燃機関の運転状態を表すパラメータ（アクセル開度など）を特定することで、吸気通路の内部を流れる空気の流体物理量（圧力、流量、流速など）を、一様のものとして模擬することができる。
【００２４】
１DCFDモデル１８は、模擬の対象である空間内の流体の流れを１次元の流れ（流れの方向に対して垂直な平面内では一様であり、流れの方向には位置毎に一様でない流れ）として模擬するモデルである。１DCFDモデル１８によれば、内燃機関の吸気通路の設計値（吸気管長、吸気管径等）や、内燃機関の運転状態を表すパラメータ（アクセル開度など）を特定することで、吸気通路の内部を流れる空気の物理量（圧力、流量、流速など）を、流れ方向の位置毎に模擬することができる。
【００２５】
集中定数系流体モデル１６による模擬の結果と、1DCFDモデル１８による模擬の結果との差は、特に、吸気通路内に吸気の脈動が発生する場合に顕著となる。すなわち、集中定数系流体モデル１６では、空間内の流体が一様なものとして模擬されるため、その結果に、吸気の脈動の影響は表れない。これに対して、1DCFDモデル１８の結果には、脈動の影響が正しく表れる。
【００２６】
空気モデル１４は、上述した２つのモデルの一方を選択して、吸気通路内の空気の流体物理量を模擬することができる。そして、空気モデル１４は、その流体物理量（例えば、空気流量）の模擬値に基づいて、筒内流入空気量を算出することができる。
【００２７】
本実施形態のシステムは、更に、トルクモデル２０、車両モデル２２、及び演算制御ユニット２４を備えている。トルクモデル２０は、内燃機関の挙動を模擬するモデルであり、筒内流入空気量に対して発生するトルクを演算することができる。また、車両モデル２２は、内燃機関のトルクを受けて、車両がどのような挙動を示すかを模擬するモデルであり、トルクの提供を受けることにより、車速や機関回転数Neを演算することができる。そして、演算制御ユニット２４は、システムの全体を制御するためのユニットである。図１に示す各ブロックは、相互にデータが授受できるように構成されている。
【００２８】
［実施の形態１の性能予測解析方法］
上述した通り、本実施形態のシステムは、空気モデル１４の中に、集中定数系流体モデル１６と、1DCFDモデル１８を含んでいる。本実施形態の性能予測解析方法は、これらのモデルを適宜選択的に用いることにより、性能予測解析の精度を下げることなく、模擬のための演算負荷を軽量化する点に特徴を有している。
【００２９】
図２は、トルクモデル２０によって算出されるトルクの結果と、機関回転数Neとの間に成立する関係を示す。より具体的には、図２中に示す実線は、1DCFDモデル１８によって筒内流入空気量を算出した場合に、機関回転数Neとトルクの間に成立する関係を示している。また、図２中に示す破線は、集中定数系流体モデル１６（0Dモデル）によって筒内流入空気量を算出した場合に両者間に成立する関係を示す。
【００３０】
上述した通り、集中定数系流体モデル１６による演算の結果と、1DCFDモデル１８による演算の結果との差は、吸気通路内に脈動が発生する場合に顕著となる。そして、現実の内燃機関において、筒内流入空気量は、機関回転数Neが十分に低い低回転領域、及び機関回転数Neが十分に高い高回転領域において、吸気の脈動の影響を殆ど受けない。換言すると、現実の内燃機関において、筒内流入空気量は、中回転領域においてのみ、吸気の脈動の影響を大きく受ける。以下、筒内流入空気量に吸気の脈動が実質的な影響を与える下限の回転数を「脈動下限回転数Ne1」と称し、また、脈動が実質的な影響を与える上限の回転数を「脈動上限回転数Ne2」と称する。
【００３１】
上記の理由により、機関回転数Neが脈動下限回転数Ne1より低い低回転領域、及び機関回転数Neが脈動上限回転数Ne2より高い高回転領域では、1DCFDモデル１８を用いた場合に得られる筒内流入空気量（つまりトルク）と、集中定数系流体モデル１６を用いた場合に得られる筒内流入空気量とは、ほぼ同様となる。他方、機関回転数Neが脈動下限回転数Ne1と脈動上限回転数Ne2との間にある場合は、集中定数系流体モデル１６によるトルク（図２中の破線）には、1DCFDモデル１８によるトルクに対するずれが重畳する。
【００３２】
集中定数系流体モデル１６の演算は、1DCFDモデル１８の演算に比して軽い負荷で行うことができる。このため、模擬の結果に差がなければ、1DCFDモデル１８ではなく、集中定数系流体モデル１６を用いて演算を行うことが適切である。そこで、本実施形態では、吸気に脈動が生じない低回転領域と高回転領域では、集中定数系流体モデル１６を用いて性能予測解析を行い、吸気の脈動が生ずる中回転領域では、1DCFDモデル１８を用いて性能予測解析を行う。
【００３３】
図３は、本実施形態のシステムの動作を説明するためのブロック図である。本実施形態のシステムは、規定の走行モードを入力することにより、車両がその走行モードで走行した場合に、内燃機関がどのような性能を発揮するかを予測・解析することができる。この機能を実現するため、本実施形態のシステムは、走行モードを記憶し、その走行モードに従って走行するための目標車速を算出し、更に、その目標車速を実現するためのアクセル操作量を算出する機能を有している。
【００３４】
上記の機能によって算出されたアクセル操作量（開度）は、図３に示すように、空気モデル１４に提供される。空気モデル１４は、車両モデル２２によって算出される機関回転数Neに基づいて、筒内流入空気量に脈動の影響が及ぶ条件が成立しているか否かを判断し、その結果に基づいて、集中定数系流体モデル１６及び1DCFDモデル１８の一方を選択する。そして、空気モデル１４は、選択したモデルを用いて、アクセル開度等に基づいて、筒内流入空気量を算出する。尚、空気モデル１４における処理の内容は、後に図４を参照してより詳しく説明する。
【００３５】
空気モデル１４において算出された筒内流入空気量は、トルクモデル２０に提供される。トルクモデル２０は、その筒内流入空気量に基づいて、内燃機関が発生するトルクを演算する。トルクモデル２０によって演算されたトルクは、車両モデル２２に提供される。そして、車両モデル２２は、提供されたトルクに基づいて、機関回転数Neや車速の模擬値を算出する。
【００３６】
車両モデル２２において算出された機関回転数Neの模擬値は、上述した通り空気モデル１４に提供される。また、車両モデル２２において算出された車速の模擬値は、アクセル操作量を算出するための基礎とされる。つまり、空気モデル１４に提供されるアクセル操作量は、車速の模擬値と、目標車速との差に基づいて算出される。
【００３７】
図４は、空気モデル１４が実行する処理の流れを示すフローチャートである。図４に示すように、空気モデル１４は、先ず、アクセル操作量に基づいて決められるスロットル開度TAや、車両モデル２２によって算出される機関回転数Neなど、吸気通路の空気流量を算出するために必要なパラメータを取得する（ステップ１００）。
【００３８】
次に、機関回転数Neが、脈動下限回転数Ne1と、脈動上限回転数Ne2との間にあるかが判別される（ステップ１０２）。脈動下限回転数Ne1及び脈動上限回転数Ne2は、予め定められた固定値としても、或いは、負荷率や筒内流入空気量をパラメータとする変動値としてもよい。
【００３９】
上記の判別の結果、Ne1＜Ne＜Ne2の成立、つまり、脈動条件の成立が認められた場合は、演算の精度を確保するため、1DCFDモデル１６によって筒内流入空気量が算出される（ステップ１０４）。他方、上記の条件が成立しないと判断された場合は、脈動の影響を考慮する必要がないため、集中定数系流体モデル１８（ODモデル）によって筒内流入空気量が算出される（ステップ１０６）。
【００４０】
以上説明した処理によれば、脈動の影響を考慮する必要がない状況下では、負荷の軽い集中定数系流体モデル１６を用いた演算が行われ、負荷の重い1DCFDモデル１８を用いた演算は、脈動の影響を考慮する必要がある場合に限って行われる。このため、本実施形態のシステムによれば、吸気通路における脈動の影響を無視することなく、かつ、演算負荷の十分な軽減を図りつつ、内燃機関の性能を正確に予測解析することができる。
【００４１】
ところで、上述した実施の形態１においては、集中定数系流体モデル１６及び1DCFDモデル１８が、吸気通路内の空気流量を模擬することとしているが、模擬の対象は、これに限定されるものではない。例えば、排気通路における排気流量を模擬の対象としてもよい。また、模擬の対象となる物理量は、流量に限られるものではなく、圧力や流速をその対象としてもよい。
【００４２】
また、上述した実施の形態１においては、脈動条件の成否を、機関回転数Neに基づいて判断することとしているが、本発明はこれに限定されるものではない。例えば、スロットル開度TAや負荷率などに基づいて脈動の有無が判断できる場合には、それらの変数に基づいて脈動条件の成否を判断することとしてもよい。
【００４３】
尚、上述した実施の形態１においては、ステップ１０２が前記第４の発明における「脈動条件判定手段」に、ステップ１０４が前記第４の発明における「１次元模擬手段」及び「空気量算出手段」に、ステップ１０６が前記第４の発明における「０次元模擬手段」及び「空気量算出手段」にそれぞれ相当している。
【００４４】
また、上述した実施の形態１においては、トルクモデル２０が前記第７の発明における「トルク模擬手段」に、車両モデル２２が前記第７の発明における「車両模擬手段」に、車両モデル２２において算出された機関回転数Neの模擬値を空気モデル１４に提供する処理が前記第７の発明における「機関回転数提供手段」に、アクセル操作量の基礎として走行モードを記憶しておくことが前記第７の発明における「走行モード記憶手段」に、アクセル操作量を算出することが前記第７の発明における「アクセル開度算出手段」に、算出したアクセル開度を空気モデル１４に提供する処理が前記第７の発明における「アクセル開度提供手段」に、それぞれ相当している。
【図面の簡単な説明】
【００４５】
【図１】本発明の実施の形態１の設計値最適化システムの構成を説明するためのブロック図である。
【図２】図１に示すトルクモデルによって算出されるトルクの結果と、機関回転数Neとの間に成立する関係を示す。
【図３】本発明の実施の形態１のシステムの動作を説明するためのブロック図である。
【図４】図３に示す空気モデルが実行する処理の流れを示すフローチャートである。
【符号の説明】
【００４６】
１４  空気モデル
１６  集中定数系流体モデル（0Dモデル）
１８  1DCFDモデル
２０  トルクモデル
２２  車両モデル
Ne1  脈動下限回転数
Ne2  脈動上限回転数
【図１】
0004710689.tif 000002
【図２】
0004710689.tif 000003
【図３】
0004710689.tif 000004
【図４】
0004710689.tif 000005