(57)【特許請求の範囲】
【請求項１】
ブースト圧と吸気温度とエンジン回転数を入力変数とし体積効率を定数とする基本式により吸入空気量を演算する吸入空気量演算方法において、
あらかじめ基準運転状態にて測定された体積効率を前記基本式に使用し、
エンジンパラメータの変動による体積効率の変動分を補正するために、あらかじめ基準運転状態にて測定された当該吸気温度と現在の当該吸気温度との比の冪乗である吸気温度補正項、あらかじめ基準運転状態にて測定された冷却水水温と現在の当該冷却水水温との比の冪乗である冷却水水温補正項、及びあらかじめ基準運転状態にて測定された吸気対排気圧力比と現在の当該吸気対排気圧力比との比の冪乗である吸気対排気圧力比補正項の何れか２つ以上の補正項により、前記基本式の演算値を補正し、
各補正項の補正乗数（冪数）は、前記基準運転状態から吸気温度と冷却水水温と吸気対排気圧力比の何れか１つを変動させた運転状態にて吸入空気量を測定し、その吸入空気量から前記基本式を逆算して体積効率を求めると共に吸気温度と冷却水水温と吸気対排気圧力比の何れか１つの変動分に対する体積効率の変動の割合から各補正項の補正乗数の仮値を求め、前記基本式と各補正項の積を吸入空気量とする演算式に各補正項の補正乗数の仮値を代入し、その演算式を逆算して体積効率を求めると共に吸気温度と冷却水水温と吸気対排気圧力比の何れか１つの変動分に対する体積効率の変動の割合から各補正項の補正乗数の仮値を１回以上更新して計算する
ことを特徴とする吸入空気量演算方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、エンジンの運転状態によらず吸入空気量が精度よく演算できる吸入空気量演算方法に関する。
【背景技術】
【０００２】
吸入空気量は、エンジン制御に重要なエンジンパラメータの一つである。吸入空気量を検出するために、エンジンの吸気管に空気流量センサ（Mass Airflow sensor；以下、ＭＡＦセンサという）が設置される。これとは別に吸入空気量以外のエンジンパラメータを使用して、演算により吸入空気量を求めるスピードデンシティ方式が知られている。スピードデンシティ方式を導入すると、スピードデンシティ方式による演算値とＭＡＦセンサによる検出値とを比較してＭＡＦセンサの診断や校正をしたり、演算値を検出値に置き換えてエンジン制御に使用したりすることができる。さらには、ＭＡＦセンサを搭載しないで演算値のみをエンジン制御に使用する廉価な車両を提供することも可能となる。
【０００３】
従来技術では、吸入空気量は、演算式（１）で与えられる。
【０００４】
【数１】
0005899685.tif 000002
【０００５】
ただし、
Mair   ＝吸入空気量
Rair   ＝空気のガス定数
TAIm  ＝吸気温度（インマニ温度）
PAIm  ＝吸気圧力（ブースト圧力、インマニ圧力）
PA_A  ＝大気圧
Vcyl   ＝総排気量
Ne    ＝エンジン回転数
η    ＝体積効率
２７３＝摂氏温度を絶対温度に変換する値
２    ＝４気筒における１回転での吸気行程回数
６０  ＝ｒｐｍで表されるエンジン回転数を１秒当たりに変換する値
である。
【０００６】
空気のガス定数Rairは固定値である。吸気温度TAImは吸気マニホールドにおける吸気の温度であり、吸気温度センサにより検出される。吸気温度TAImは摂氏で検出されるので２７３を加えて絶対温度とする。吸気圧力PAImは吸気マニホールドにおける吸気の圧力であり、ブースト圧センサにより検出される。大気圧PA_Aは大気圧センサにより検出される。吸気圧力PAImは、ゲージ圧力として検出されるので、大気圧PA_Aを足すことで絶対圧力となる。総排気量Vcyl はエンジンに固有の値である。エンジン回転数Neはクランク角センサにより検出される。体積効率ηは、シリンダ断面積とピストンストローク長で決まる機械的なシリンダの容量と実際にシリンダに吸入される空気量との比率で表される定数である。
【０００７】
このように、吸入空気量の演算式（１）は、吸気温度TAImと吸気圧力PAImとエンジン回転数Neとを入力変数とし、体積効率ηを定数とする。
【先行技術文献】
【特許文献】
【０００８】
【特許文献１】特開２０１０−１１６８５７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００９】
ところが、本発明者は、前述の入力変数とそれ以外も含む各種のエンジンパラメータを変動させてエンジンを運転し、そのとき得られる入力変数を用いて吸入空気量演算を行い、演算値を評価する実験を行ったところ、特定のエンジンパラメータが変動すると吸入空気量の誤差が変動することを見出し、この原因は、当該エンジンパラメータが表しているエンジンの運転状態に応じて体積効率が変動するからではないかと推測した。体積効率がエンジンの運転状態に応じて変動するので、吸入空気量の演算値に誤差が生じると考えられる。
【００１０】
そこで、本発明の目的は、上記課題を解決し、エンジンの運転状態によらず吸入空気量が精度よく演算できる吸入空気量演算方法を提供することにある。
【課題を解決するための手段】
【００１１】
上記目的を達成するための本発明は、ブースト圧と吸気温度とエンジン回転数を入力変数とし体積効率を定数とする基本式により吸入空気量を演算する吸入空気量演算方法において、あらかじめ基準運転状態にて測定された体積効率を前記基本式に使用し、エンジンパラメータの変動による体積効率の変動分を補正するために、あらかじめ基準運転状態にて測定された当該吸気温度と現在の当該吸気温度との比の冪乗である吸気温度補正項、あらかじめ基準運転状態にて測定された冷却水水温と現在の当該冷却水水温との比の冪乗である冷却水水温補正項、及びあらかじめ基準運転状態にて測定された吸気対排気圧力比と現在の当該吸気対排気圧力比との比の冪乗である吸気対排気圧力比補正項の何れか２つ以上の補正項により、前記基本式の演算値を補正し、各補正項の補正乗数（冪数）は、前記基準運転状態から吸気温度と冷却水水温と吸気対排気圧力比の何れか１つを変動させた運転状態にて吸入空気量を測定し、その吸入空気量から前記基本式を逆算して体積効率を求めると共に吸気温度と冷却水水温と吸気対排気圧力比の何れか１つの変動分に対する体積効率の変動の割合から各補正項の補正乗数の仮値を求め、前記基本式と各補正項の積を吸入空気量とする演算式に各補正項の補正乗数の仮値を代入し、その演算式を逆算して体積効率を求めると共に吸気温度と冷却水水温と吸気対排気圧力比の何れか１つの変動分に対する体積効率の変動の割合から各補正項の補正乗数の仮値を１回以上更新して計算する吸入空気量演算方法である。
【発明の効果】
【００１６】
本発明は次の如き優れた効果を発揮する。
【００１７】
（１）エンジンの運転状態によらず吸入空気量が精度よく演算できる。
【図面の簡単な説明】
【００１８】
【図１】本発明の吸入空気量演算方法が適用されるエンジン及び吸排気系の構成図である。
【図２】従来の吸入空気量演算方法による演算値の誤差を示すグラフである。
【図３】本発明の吸入空気量演算方法により、吸気温度補正項、冷却水水温補正項、吸気対排気圧力比補正項を全て使用した場合の演算値の誤差を示すグラフである。
【図４】本発明の吸入空気量演算方法により、吸気温度補正項、冷却水水温補正項を使用した場合の演算値の誤差を示すグラフである。
【発明を実施するための形態】
【００１９】
以下、本発明の一実施形態を添付図面に基づいて詳述する。
【００２０】
図１に示されるように、エンジン１の吸気管２にＭＡＦセンサ３が設置される。吸気管２には、ＭＡＦセンサ３の下流にターボチャージャ４のコンプレッサ５が接続され、コンプレッサ５の下流にインタークーラ６が接続される。インタークーラ６の下流では、吸気管２に吸気スロットル７が設置され、吸気スロットル７の下流で吸気管２がエンジン１の吸気マニホールド８に接続される。エンジン１の排気マニホールド９に排気管１０が接続され、排気管１０の下流はターボチャージャ４のタービン１１に接続される。タービン１１には、過給量の調整弁１２が設けられる。排気マニホールド９と吸気マニホールド８の間には、ＥＧＲ管１３が設けられる。ＥＧＲ管１３には、ＥＧＲクーラ１４とＥＧＲ弁１５が設けられる。
【００２１】
吸気管２には、吸気温度センサ１６が設けられる。吸気温度センサ１６の設置場所は、図示のように吸気スロットル７の入口側でもよいし、吸気スロットル７の出口側でもよい。吸気温度を検出する吸気マニホールド８には、吸気圧力を検出するブースト圧センサ１７が設けられる。エンジン１のクランク軸（図示せず）には、クランク軸の回転を検出するクランク角センサ１８が設けられる。エンジン１の冷却水循環系の、例えば、ウォータジャケット（図示せず）に、冷却水の温度を検出する冷却水水温センサ１９が設けられる。排気マニホールド９には、排気圧力を検出する排気圧力センサ２０が設けられる。大気圧を検出する大気圧センサ２１とエンジン１の電子制御を行う電子制御装置（Electronical Control Unit；ＥＣＵ）２２は、エンジン１を搭載した車両内の適宜な場所に設置される。ＥＣＵ２２には、吸入空気量の演算プログラムとその演算に使用する定数やマップが格納される。
【００２２】
本発明の吸入空気量演算方法における吸入空気量の演算式（２）は、基本式と３つの補正項とからなる。基本式は、従来技術で述べた演算式（１）と同等である。補正項は、体積効率の変動をもたらすエンジンパラメータごとに独立に設定され、いずれも基本式に対する係数の項（かけ算する項）となる。本実施形態では、体積効率の変動をもたらすエンジンパラメータとして、吸気温度と冷却水水温と吸気対排気圧力比を使用している。したがって、３つの補正項は、吸気温度補正項、冷却水水温補正項、吸気対排気圧力比補正項となる。
【００２３】
【数２】
0005899685.tif 000003
【００２４】
ただし、
Mair   ＝吸入空気量
Rair   ＝空気のガス定数
TAIm  ＝吸気温度（インマニ温度）
PAIm  ＝吸気圧力（ブースト圧力、インマニ圧力）
PA_A  ＝大気圧
Vcyl   ＝総排気量
Ne    ＝エンジン回転数
η0  ＝基準体積効率
２７３＝摂氏温度を絶対温度に変換する値
２    ＝４気筒における１回転での吸気行程回数
６０  ＝ｒｐｍで表されるエンジン回転数を１秒当たりに変換する値
TAIm0  ＝基準吸気温度（基準インマニ温度）
m      ＝吸気温度補正乗数
THW   ＝冷却水水温
THW0  ＝基準冷却水水温
a       ＝冷却水水温補正乗数
PETbI   ＝排気圧力（タービン入口圧力）
PAIm0  ＝基準吸気圧力（基準ブースト圧力、基準インマニ圧力）
PETbI0  ＝基準排気圧力（基準タービン入口圧力）
PA_A0   ＝基準大気圧
n       ＝圧力比補正乗数
である。
【００２５】
基準体積効率η0、基準吸気温度TAIm0、基準冷却水水温THW0、基準吸気圧力PAIm0、基準排気圧力PETbI0、基準大気圧PA_A0には、実験により測定及び演算して求められた値が設定される。吸気温度補正乗数m、冷却水水温補正乗数a、圧力比補正乗数nには、実験により演算して求められた値が設定される。ただし、ＥＣＵ２２に基本演算命令として冪乗演算命令が搭載されていない場合は、冪乗演算については、近似式あるいはマップをＥＣＵ２２に格納することになる。
【００２６】
設定する数値を求めるための実験は、燃料量、エンジン回転数、吸入空気温度、冷却水水温など所望するエンジンパラメータをエンジン外部から任意の値に制御できるエンジンテストベンチ等の実験設備にて行われる。
【００２７】
基準体積効率η0とその他の基準値を求める方法を詳しく説明する。
【００２８】
エンジン１が定常運転状態（過渡的でない運転状態）のときに、精密に校正された測定器で吸入空気量Mairを実測し、この吸入空気量Mairと吸気温度TAImと吸気圧力PAImとエンジン回転数Neとから基本式を逆算して基準体積効率η0を求める。燃料量とエンジン回転数とをそれぞれ適宜な刻みで変更し、定常運転状態にて実測した吸入空気量Mairから基本式を逆算して基準体積効率η0を求め、各基準体積効率η0を燃料量とエンジン回転数を参照軸としてマップに設定する。ただし、燃料量がない状態では、外力でエンジン１を回転させるモータリング運転により強制的に吸排気させ、実測した吸入空気量Mairから基本式を逆算して基準体積効率η0を求める。
【００２９】
基準体積効率η0を求めたときの運転状態を基準運転状態とする。基準運転状態の際の吸入空気温度、冷却水水温、吸気圧力、排気圧力、大気圧の値を基準値とする。これにより、基準吸気温度TAIm0、基準冷却水水温THW0、基準吸気圧力PAIm0、基準排気圧力PETbI0、基準大気圧PA_A0が得られる。これらの基準値は、基準体積効率η0と同様に燃料量とエンジン回転数を参照軸としてマップに設定する。
【００３０】
次に、吸気温度補正乗数m、冷却水水温補正乗数a、圧力比補正乗数nを求める方法を詳しく説明する。
【００３１】
まず、ステップＳ１にて、エンジン１を基準運転状態で運転する。
【００３２】
次いで、ステップＳ２にて、吸気温度、冷却水水温、吸気対排気圧力比のうち所望する１つのエンジンパラメータが変動するよう、運転状態を変動させる。当該エンジンパラメータが変動した運転状態で、精密に校正された測定器で吸入空気量Mairを実測し、吸入空気量Mairから基本式を逆算して体積効率を求め、当該エンジンパラメータの変動分に対する体積効率の変動の割合（例えば、グラフの傾き）から補正乗数の仮値を求める。
【００３３】
例えば、吸気温度が変動するように運転状態を変動させたとき、吸気温度の変動分に対する体積効率の変動の割合から吸気温度補正乗数mの仮値を求めることができる。冷却水水温が変動するように運転状態を変動させたとき、冷却水水温の変動分に対する体積効率の変動の割合から冷却水水温補正乗数aの仮値を求めることができる。吸気対排気圧力比が変動するように運転状態を変動させたとき、吸気対排気圧力比の変動分に対する体積効率の変動の割合から圧力比補正乗数nの仮値を求めることができる。
【００３４】
しかし、ステップＳ２で補正乗数の仮値を求めるにあたり、所望する１つのエンジンパラメータを変動させようと運転状態を変動させた際、純粋に所望したエンジンパラメータのみを変動させることは不可能であり、所望しないエンジンパラメータも随伴して変動してしまう。つまり、ステップＳ２で求めた補正乗数仮値は、所望しないエンジンパラメータの変動の影響を含んだものである。
【００３５】
ステップＳ３では、ステップＳ２と同様に所望する１つのエンジンパラメータが変動するよう、運転状態を変動させ、当該エンジンパラメータが変動した運転状態で、精密に校正された測定器で吸入空気量Mairを実測する。この吸入空気量Mairを用い、基本式に当該補正乗数を求める対象以外の補正項が付いた演算式を逆算して体積効率を求める。
【００３６】
例えば、吸気温度補正乗数mを求める場合は、基本式と冷却水水温補正項と吸気対排気圧力比補正項の積が吸入空気量Mairである演算式を用い、各補正項には既に求めた補正乗数仮値を代入し、演算式を逆算して体積効率を求める。この体積効率からステップＳ２と同様に補正乗数仮値を求めるが、このとき求まる補正乗数仮値はステップＳ２で求めた補正乗数仮値に比べて、随伴して変動してしまった所望しないエンジンパラメータの変動の影響が少ない補正乗数仮値に更新されている。
【００３７】
冷却水水温補正乗数aを求める場合は、基本式と吸気温度補正項と吸気対排気圧力比補正項の積が吸入空気量Mairである演算式を用い、吸気温度補正乗数mを求めた場合と同様にして補正乗数仮値を更新する。圧力比補正乗数nを求める場合は、基本式と吸気温度補正項と冷却水水温補正項の積が吸入空気量Mairである演算式を用い、同様に補正乗数仮値を更新する。
【００３８】
ステップＳ３は、１回に限らず、複数回繰り返して各補正乗数仮値を順次更新してもよい。最終的に更新された補正乗数仮値を補正乗数として終了する。これにより、随伴して変動してしまった所望しないエンジンパラメータの変動の影響を可能な限り取り除いた状態にて、精度よく補正乗数を求めることができる。繰り返しの回数を多くすることで精度が増すが、回数は要求精度に応じて何回にしてもよく、特に回数は限定しない。
【００３９】
以上の手順により、各補正乗数が燃料量とエンジン回転数を参照軸とするマップに設定される。
【００４０】
次に、本発明の吸入空気量演算方法の原理と作用効果を説明する。
【００４１】
本発明者は、堆積効率変動の要因となるエンジンパラメータが変動したときの演算値の誤差を評価する実験を行った。具体的には、吸気温度と冷却水水温と吸気対排気圧力比の３つのエンジンパラメータを個別に変動させてエンジン１を運転し、吸気温度TAImと吸気圧力PAImとエンジン回転数Neの入力変数を用いて従来の演算式（１）による吸入空気量演算を行い、一方で精密に校正された測定器で吸入空気量を測定し、測定値と演算値の誤差を評価する実験を行った。例として、表１に、実験番号＃０〜＃１２におけるエンジンパラメータの設定を示す。このように、冷却水水温については、標準（８０℃）と８５℃と４５℃の３段階を設定し、吸気温度については厳密な数値で設定できないので、標準と標準より高側と標準より低側の３段階を設定し、吸気対排気圧力比については排気抵抗の増加を想定し、無負荷最高回転時の排気圧力を調整し、標準と＋１０ｋＰａの２段階に設定するものとした。
【００４２】
【表１】
0005899685.tif 000004
【００４３】
図２に示されるように、実験番号＃０〜＃１２の各実験において、測定値と演算値には誤差が生じた。各実験における誤差は、−３％から＋２％の広い幅で分布している。これは、体積効率が実際はエンジン１の運転状態に応じて変動するにもかかわらず、演算式（１）では体積効率ηを固定値としているからである。
【００４４】
吸気温度について考察すると、気体は吸気マニホールド８からエンジン１の各シリンダに吸気ポートを経由して吸入される。このような吸気の動作がエンジン回転数に応じ１秒間に数十回、繰り返し行われる。このとき、気体の吸入されやすさは気体の粘性の影響を受ける。気体の粘性は、温度に依存するので、吸気温度によって気体の粘性が変動して、シリンダへの気体の吸入されやすさが変動する。このため、吸気温度の変動によって体積効率が変動することになる。
【００４５】
冷却水水温について考察すると、シリンダや吸気マニホールド８の壁温度は、通常、大気温度より高い。吸気管から吸気マニホールド８を経由してシリンダに入ってきた気体は、壁に暖められて膨張するため、その後から入ろうとする気体を阻害する。壁温度によって、シリンダへ入ろうとする気体への阻害の大きさが異なる。冷却水水温は壁温度を反映しており、したがって、冷却水水温の変動によって体積効率が変動することになる。
【００４６】
吸気対排気圧力比について考察すると、排気管１０には、排気ガスフィルタなどの後処理装置（図示せず）が設けられる。後処理装置の詰まり具合によって、シリンダから出ていく排気量が変動する。すなわち、後処理装置に捕集物が多く溜まっていると、排気抵抗が大きくなるため、シリンダから出ていかない気体が存在することになる。このような現象を内部ＥＧＲという。内部ＥＧＲが起きると、シリンダに気体が入りにくくなり、体積効率が変動する。内部ＥＧＲの程度は、吸気圧力と排気圧力の比で表すことができる。つまり、吸気対排気圧力比の変動によって体積効率が変動することになる。
【００４７】
ここまでの考察によれば、体積効率は、吸気温度と冷却水水温と吸気対排気圧力比のいずれの変動によっても変動し、しかも各々は個別の事象である。そこで、本発明者は、基準運転状態にて測定された体積効率（基準値）を使用する基本式で演算を行い、その演算値をこれらのエンジンパラメータごとに補正することを考えた。補正項は、基本式に対して係数となる項とし、あらかじめ基準運転状態にて測定された当該エンジンパラメータ（基準値）と現在の当該エンジンパラメータ（センサ値）との比に基づくものとした。比で表された項は、分母である基準値より分子であるセンサ値が大きければ１より大きく、基準値よりセンサ値が小さければ１より小さくなる。この比の冪乗をとることで、補正項を体積効率の変動分に合わせ込むことができる。すなわち、冪乗に用いる補正乗数を、実験時の測定結果と演算結果が一致するよう近似させて求める。具体的な補正乗数の求め方は既に説明した通りである。
【００４８】
このようにして、本発明の吸入空気量演算方法では、基本式に補正項を付加した吸入空気量の演算式（２）を使用することになった。
【００４９】
図３の実験結果は、吸気温度補正項、冷却水水温補正項、吸気対排気圧力比補正項を全て使用した場合のものである。前述と同様の実験番号＃０〜＃１２の実験を実施し、演算式（２）による吸入空気量演算を行い、測定値と演算値の誤差を評価したところ、図３に示されるように、全ての実験番号＃０〜＃１２の誤差が−０．５〜＋０．８％以内となった。図２と比較すると、エンジンパラメータの設定によらず、誤差が顕著に小さいことが分かる。これより、本発明の吸入空気量演算方法は、体積効率に影響を与える吸気温度、冷却水水温、吸気対排気圧力比が変動しても吸入空気量が精度よく演算できることが確認できた。
【００５０】
図４の実験結果は、吸気温度補正項と冷却水水温補正項を使用し、吸気対排気圧力比補正項は使用しなかった場合のものである。図３の実験結果と比較すると、エンジン回転数が低い領域では誤差は同程度であるが、高い領域において誤差の範囲が拡がる傾向にある。言い換えると、吸気対排気圧力比補正は、エンジン回転数が低い領域では大きな効果は示さないがエンジン回転数が高い領域で効果が顕著である。このことから、エンジン１の用途が高速回転を必要としない用途である場合、吸気対排気圧力比補正を省略しても吸入空気量の演算は誤差が小さいことが分かる。
【００５１】
この他にも、吸気温度補正のみを行う実験、冷却水水温補正のみを行う実験、吸気対排気圧力比補正のみを行う実験を行ったところ、冷却水水温補正のみの補正効果として、基本式による誤差が−５．１％であったものが０．９％に改善され、吸気温度補正のみの補正効果として、 基本式による誤差が１．５％であったものが０．２％に改善され、吸気対排気圧力比補正のみの補正効果として、基本式による誤差が２．１％であったものが−０．１％に改善された。
【００５２】
以上説明したように、従来の吸入空気量演算方法では、エンジン１の運転状態によらず演算式（１）の体積効率ηが固定値であったため、エンジン１の運転状態に応じて実際の体積効率が変動すると、吸入空気量が精度よく演算できなかったのに対し、本発明の吸入空気量演算方法によれば、基準体積効率η0を使用する基本式に、エンジンパラメータの変動による体積効率の変動分を補正する補正項を付加したので、体積効率の変動に対応でき、吸入空気量が精度よく演算できる。
【００５３】
本発明の吸入空気量演算方法によれば、１つのエンジンパラメータごとに補正項を設けたので、複数のエンジンパラメータの変動に個別に対応でき、しかも、重ね合わせることが可能となる。なお、実施形態では、体積効率の変動に寄与するエンジンパラメータとして、吸気温度と冷却水水温と吸気対排気圧力比とを取り扱ったが、この他にも体積効率の変動に寄与するエンジンパラメータがあれば、同様の補正項を設けることで対応できる。
【００５４】
本発明の吸入空気量演算方法によれば、補正項が基準運転状態にて測定されたエンジンパラメータと現在のエンジンパラメータとの比を補正乗数で冪乗して求められるので、エンジンパラメータの変動による体積効率の変動が補正項に適切に反映される。なお、複数機種のエンジンに対する実験から、補正乗数はエンジン機種に依存せず、同じような値になることが分かった。
【符号の説明】
【００５５】
１  エンジン
２  吸気管
３  ＭＡＦセンサ
８  吸気マニホールド
９  排気マニホールド
１６  吸気温度センサ
１７  ブースト圧センサ
１８  クランク角センサ
１９  冷却水水温センサ
２０  排気圧力センサ
２１  大気圧センサ
【図１】
0005899685.tif 000005
【図２】
0005899685.tif 000006
【図３】
0005899685.tif 000007
【図４】
0005899685.tif 000008