(57)【要約】    【課題】吸気弁の作動特性の変更による吸入空気量の制御において、マップの多用を回避し、演算負荷を軽減する。    【解決手段】ソニック流として吸入した場合に得られるシリンダ吸入空気量を仮想ソニック吸入空気量ＱDとし、かつ吸気の開始から終了までの行程容積を吸気弁上流における吸入空気の密度及び温度で充填した場合に得られるシリンダ吸入空気量を理論最大吸入空気量ＱMAXとして、実際のシリンダ吸入空気量Ｑｃｙｌに関して第１の比ＱD／ＱMAXと第２の比Ｑｃｙｌ／ＱMAXとの間の一義的な関係を設定し、コントロールユニットに格納する。実際の運転時では、目標新気量ｔＱｃｙｌをシリンダ吸入空気量Ｑｃｙｌとして第２の比を算出し、これに基づいて前記関係（Ｂ１０１）により仮想ソニック吸入空気量ｔＱDを算出し、これに基づいて目標吸気作動角ｔθｅｖｅｎｔを設定する。【選択図】  図４
2006105098.tif 000002
【特許請求の範囲】
【請求項１】
エンジンの筒内に吸入される空気の量であるシリンダ吸入空気量を、吸気弁の作動特性を変更して制御する装置であって、
目標とするシリンダ吸入空気量である目標新気量を算出する目標新気量算出手段と、
算出された目標新気量に基づいて吸気弁の作動特性に関する所定の制御パラメータを設定するパラメータ設定手段と、を含んで構成され、
前記パラメータ設定手段は、吸気弁の作動特性に応じた開口面積でソニック流として吸入した場合に得られるシリンダ吸入空気量を仮想ソニック吸入空気量ＱDとし、かつ吸気の開始から終了までの行程容積を吸気弁上流における吸入空気の密度及び温度で充填した場合に得られるシリンダ吸入空気量を理論最大吸入空気量ＱMAXとして、実際のシリンダ吸入空気量Ｑｃｙｌに関して第１の比ＱD／ＱMAXと第２の比Ｑｃｙｌ／ＱMAXとの間の一義的な関係が設定され、エンジンの運転時において、前記目標新気量と前記関係とに基づいて前記仮想ソニック吸入空気量を算出するとともに、算出した仮想ソニック吸入空気量に基づいて前記制御パラメータを演算及び設定するエンジンのシリンダ吸入空気量制御装置。
【請求項２】
前記パラメータ設定手段は、前記目標新気量をシリンダ吸入空気量Ｑｃｙｌとして前記第２の比を算出し、算出した第２の比に基づいて前記関係により前記仮想ソニック吸入空気量を算出する請求項１に記載のエンジンのシリンダ吸入空気量制御装置。
【請求項３】
前記目標新気量は、目標とするエンジントルクに対応するシリンダ吸入空気量である請求項１又は２に記載のエンジンのシリンダ吸入空気量制御装置。
【請求項４】
前記目標新気量算出手段により算出された目標新気量に対し、吸気通路内における気柱振動に起因するシリンダ吸入空気量の変動分を排除する補正を施す目標新気量補正手段を更に含んで構成される請求項１〜３のいずれかに記載のエンジンのシリンダ吸入空気量制御装置。
【請求項５】
前記目標新気量補正手段は、筒内へ向かう吸入空気の流れがチョークする第１の領域とこの第１の領域以外の第２の領域とが定められ、前記第１の領域では、前記変動分に応じた目標新気量の補正量を実質的に０とする一方、前記第２の領域では、前記気柱振動分を加味した実際の吸気圧力の逆数に比例し、かつ前記気柱振動分を加味した実際の吸気温度に比例する特性により、前記補正を行う請求項４に記載のエンジンのシリンダ吸入空気量制御装置。
【請求項６】
前記第１及び第２の比の関係を設定するうえで、前記仮想ソニック吸入空気量の算出期間として、筒内圧力が吸気圧力に低下する実効上死点以降の期間が採用されるとともに、
前記目標新気量算出手段により算出された目標新気量から、吸気弁開期間と排気弁開期間とのオーバーラップ期間に筒内から吸気通路内に吹き返す吹返ガス量を減算する吹返ガス量減算手段を更に含んで構成される請求項１〜５のいずれかに記載のエンジンのシリンダ吸入空気量制御装置。
【請求項７】
前記パラメータ設定手段は、前記理論最大吸入空気量を、筒内圧力が吸気圧力に低下する実効上死点から、筒内で吸入空気の圧縮が実際に開始される吸気弁の実効閉時期までのシリンダ吸入空気量として算出する請求項１〜６のいずれかに記載のエンジンのシリンダ吸入空気量制御装置。
【請求項８】
エンジンの筒内に吸入される空気の量であるシリンダ吸入空気量を、吸気弁の作動特性を変更して制御する方法であって、
吸気弁の作動特性に応じた開口面積でソニック流として吸入した場合に得られるシリンダ吸入空気量を仮想ソニック吸入空気量ＱDとし、かつ吸気の開始から終了までの行程容積を吸気弁上流における吸入空気の密度及び温度で充填した場合に得られるシリンダ吸入空気量を理論最大吸入空気量ＱMAXとして、実際のシリンダ吸入空気量Ｑｃｙｌに関して第１の比ＱD／ＱMAXと第２の比Ｑｃｙｌ／ＱMAXとの間の一義的な関係を予め設定し、実際の運転に際し、目標とするシリンダ吸入空気量をＱｃｙｌとして、この目標新気量Ｑｃｙｌを前記関係により前記仮想ソニック吸入空気量に変換するとともに、変換した仮想ソニック吸入空気量に基づいて吸気弁の作動特性を変更するエンジンのシリンダ吸入空気量制御方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、エンジンのシリンダ吸入空気量制御装置及びシリンダ吸入空気量制御方法に関し、詳細には、吸気弁の作動特性の変更により吸入空気量を制御する場合に、この変更に関する制御パラメータを簡易に設定するための技術に関する。
【背景技術】
【０００２】
エンジンにおける吸入空気量の制御として、吸気弁の作動特性を変更することによるものが知られている。すなわち、スロットル弁の制御によらず、主に吸気弁の作動角（及び弁リフト）及びこの作動角の中心位相（以下「作動中心角」という。）をエンジンの運転条件に応じて変更することで、運転条件毎に要求される吸入空気量（すなわち、シリンダ吸入空気量）を達成しようとするものである（特許文献１）。
【特許文献１】特開２００２−２５６９０５号公報
【発明の開示】
【発明が解決しようとする課題】
【０００３】
しかしながら、吸入空気量の制御を吸気弁の作動特性（以下「弁作動特性」という。）の変更によることとした場合は、次のことが問題となる。すなわち、スロットル弁の制御よる場合は、その動作に対して比較的に線形に近い特性で吸入空気量が変化するのに対し、弁作動特性の変更による場合は、吸入空気量の変化が大きく、運転条件毎の緻密な適合が必要となるため、マップからの検索による場合は、その多用が避けられず、演算負荷が大きいことである。
【０００４】
本発明は、このような実情に鑑みなされたものであり、弁作動特性の変更による吸入空気量の制御において、マップの多用を回避し、演算負荷を軽減しつつ、要求される吸入空気量を達成することを目的とする。
【課題を解決するための手段】
【０００５】
本発明は、エンジンの筒内に吸入される空気の量であるシリンダ吸入空気量を制御する装置及び方法を提供する。本発明では、目標とするシリンダ吸入空気量である目標新気量を算出するとともに、算出した目標新気量に基づいて吸気弁の作動特性に関する所定の制御パラメータを設定する。この制御パラメータの設定に関し、吸気弁の作動特性に応じた開口面積でソニック流として吸入した場合に得られるシリンダ吸入空気量を仮想ソニック吸入空気量ＱDとし、かつ吸気の開始から終了までの行程容積を吸気弁上流における吸入空気の密度及び温度で充填した場合に得られるシリンダ吸入空気量を理論最大吸入空気量ＱMAXとして、実際のシリンダ吸入空気量Ｑｃｙｌに関して第１の比ＱD／ＱMAXと第２の比Ｑｃｙｌ／ＱMAXとの間の一義的な関係を設定する。エンジンの運転時では、目標新気量とこの関係とに基づいて仮想ソニック吸入空気量を算出するとともに、算出した仮想ソニック吸入空気量に基づいて制御パラメータを演算及び設定する。仮想ソニック吸入空気量は、目標新気量をシリンダ吸入空気量Ｑｃｙｌとして第２の比を算出し、算出した第２の比に基づいて前記関係により算出するとよい。
【発明の効果】
【０００６】
本発明によれば、第１及び第２の比ＱD／ＱMAX，Ｑｃｙｌ／ＱMAXを定義するとともに、これらの比の一義的な関係を設定したことで、この関係の採用によりマップの多用を回避するとともに、この関係により算出した仮想ソニック吸入空気量をもとに、要求される吸入空気量（すなわち、目標新気量）を達成することができる。
【発明を実施するための最良の形態】
【０００７】
以下に図面を参照して、本発明の実施の形態について説明する。
【０００８】
図１は、本発明の一実施形態に係る火花点火エンジン（以下、単に「エンジン」という。）１の構成を示している。
【０００９】
このエンジン１の吸気通路１０１には、電子制御式のスロットル弁１０２が設置されている。このスロットル弁１０２により吸気通路１０１に導入される空気の量を制御することも可能であるが、本実施形態では、この吸入空気量の制御を主に、後述する吸気弁１０４の作動特性の変更によることとし、スロットル弁１０２は、この作動特性による制御の前提となる吸気圧力Ｐｍの制御に採用することとしている。また、吸気通路１０１には、燃料供給用のインジェクタ１０３が設置されている。このインジェクタ１０３により、制御された吸入空気量のもとで所定の当量比を達成するのに必要な量の燃料が噴射される。吸気通路のポート部１０１ａには、ポペット型の吸気弁１０４が設置されている。この吸気弁１０４は、その上方に配置された動弁装置（以下「吸気動弁装置」という。）１０５により駆動され、この吸気弁１０４の開期間に、吸入空気及び燃料の混合気が筒内に導入される。本実施形態において、吸気動弁装置１０５は、吸気弁１０４の作動角（以下「吸気作動角」という。）θｅｖｅｎｔ及びリフト量を連続的に変更するとともに、作動中心角を連続的に変更することができる。
【００１０】
エンジン本体において、シリンダヘッドＨには、燃焼室の上部略中央に臨ませて点火プラグ１０６が設置されている。筒内に導入された混合気に対し、この点火プラグ１０６により点火が行われる。
【００１１】
燃焼後、発生した排出ガスは、排気通路１０７に送り出される。この排気通路のポート部１０７ａには、ポペット型の排気弁１０８が設置されている。この排気弁１０８は、その上方に配置された他の動弁装置１０９により駆動され、この排気弁１０８の開期間に、排出ガスの送出が行われる。なお、本実施形態では、吸気弁１０４とは異なり、排気弁１０８の作動角、リフト量及び作動角の中心位相を一定のものとしているが、吸気動弁装置１０５と同様な構成のもの又は他の公知の可変動弁装置を採用して、作動角等を変更可能に構成することもできる。
【００１２】
吸気動弁装置１０５及びスロットル弁１０２等の動作は、電子制御ユニットとして構成されるエンジンコントローラ（以下「ＥＣＵ」という。）２０１により制御される。ＥＣＵ２０１には、アクセルペダルの踏込量（アクセル開度ＡＰＯを示す。）を検出するアクセルセンサ２１１の検出信号、及びクランクシャフトの回転位置を検出するクランク角センサ２１２の検出信号（これに基づいてエンジン回転数Ｎｅを算出する。）が入力されるとともに、吸気通路１０１内（ここでは、サージタンク内）の圧力（以下「吸気圧力」という。）Ｐｍを検出する吸気圧力センサ２１３の検出信号、吸気通路１０１内の温度Ｔｍを検出する吸気温度センサ２１４の検出信号、及び排気通路内の圧力（以下「排気圧力」という。）Ｐｅを検出する排気圧力センサ２１５の検出信号等が入力される。ＥＣＵ２０１は、入力した各種の検出信号をもとに、吸気動弁装置１０５により吸気作動角及び作動中心角を、スロットル弁１０２によりスロットル開度を夫々制御する。なお、本実施形態では、ＥＣＵ２０１がシリンダ吸入空気量制御装置としての機能を備えている。
【００１３】
図２は、本実施形態に係る吸気動弁装置１０５の構成を示している。
【００１４】
この吸気動弁装置１０５は、作動角変更機構Ａと中心角変更機構Ｂとを含んで構成される。
【００１５】
吸気弁１０４の上方に駆動軸１５１が気筒列方向に延在させて設置されており、この駆動軸１５１に揺動カム１５２が相対回転可能に取り付けられている。この揺動カム１５２は、吸気弁のリフタ１４１と当接し、このリフタ１４１を介して吸気弁１０４を上下に駆動する。作動角変更機構Ａは、駆動軸１５１と揺動カム１５２とを繋ぐ後述するリンクの姿勢を変化させて、吸気作動角を変更するものである。他方、中心角変更機構Ｂは、駆動軸１５１の図示しないクランクシャフトに対する位相を変化させることで、作動中心角を変更するものである。
【００１６】
ここで、前者の作動角変更機構Ａの作動原理について、図３を参照して説明する。
【００１７】
作動角変更機構Ａは、駆動軸１５１に固定された円形の偏心駆動カム１５３と、この偏心駆動カム１５３に相対回転可能に外嵌するリング状リンク１５４と、駆動軸１５１と平行に配置された制御軸１５５と、この制御軸１５５に固定された円形の偏心制御カム１５６と、この偏心制御カム１５６に相対回転可能に外嵌し、一端でリング状リンク１５４と連結するロッカーアーム１５７と、このロッカーアーム１５７を揺動カム１５２と連結するロッド状リンク１５８とを含んで構成される。制御軸１５５は、電磁アクチュエータ１６１がギア列を駆動することにより回転する。
【００１８】
この作動角変更機構Ａの動作は、次のようである。クランクシャフトに連動して駆動軸１５１が回転すると、これに伴うリング状リンク１５４の往復動作に併せ、ロッカーアーム１５７が偏心制御カム１５６の軸心周りで揺動し、ロッド状リンク１５８により揺動カム１５２を駆動する。また、電磁アクチュエータ１６１により制御軸１５５を回転させることで、偏心制御カム１５６の軸心位置が変化し、ロッカーアーム１５７の回転中心が変位して、吸気作動角（及び弁リフト）が連続的に変化する。このため、電磁アクチュエータ１６１により制御軸１５５を操作することで、吸気作動角を連続的に変更することができる。
【００１９】
なお、中心角変更機構Ｂには、駆動軸１５１のカムスプロケットに対する位相を変化させ得る、公知のいかなる可変動弁装置が採用されてよい。本実施形態では、カムスプロケットと駆動軸１５１との間に中間ギアを介装して、これらの間にヘリカルギア列を形成し、中間ギアを前後させることにより駆動軸１５１の相対位相を変化させるものを採用している。
【００２０】
以下に、本実施形態に係るＥＣＵ２０１の構成及びこれが行う制御（主に、目標吸気作動角ｔθｅｖｅｎｔの設定）の内容について、ブロック図により説明する。
【００２１】
ＥＣＵ２０１が行う制御は、簡単には次のようである。ＥＣＵ２０１は、アクセル開度ＡＰＯ及びエンジン回転数Ｎｅ等の運転条件に基づいてエンジン１が発生すべき目標トルクｔＴｅを演算及び設定するとともに、この目標トルクｔＴｅに基づいて吸気動弁装置１０５及びスロットル弁１０２を作動させる。すなわち、ＥＣＵ２０１は、目標トルクｔＴｅを達成するのに必要なシリンダ吸入空気量として目標新気量ｔＱｃｙｌを算出するとともに、この目標新気量ｔＱｃｙｌに基づいて吸気弁１０４の目標作動角（以下「目標吸気作動角」という。）ｔθｅｖｅｎｔを設定し、吸気動弁装置１０５を作動させる。また、ＥＣＵ２０１は、実際のシリンダ吸入空気量Ｑｃｙｌを測定し、このシリンダ吸入空気量Ｑｃｙｌの目標新気量ｔＱｃｙｌに対する偏差（＝ｔＱｃｙｌ−Ｑｃｙｌ）に応じ、これを減少させる位置にスロットル弁１０２を作動させ、吸気圧力Ｐｍを調整する。
【００２２】
ＥＣＵ２０１には、目標吸気作動角ｔθｅｖｅｎｔの設定のため、ブロックＢ１０１（図４）に示すテーブルデータ（以下「変換テーブル」という。）が格納されている。この変換テーブルは、筒内圧力Ｐ０のもとで下式（１ａ）及び（１ｂ）により与えられる基準シリンダ吸入空気量ＱDと、実際のシリンダ吸入空気量（以下、単に「シリンダ吸入空気量」という。）Ｑｃｙｌ（条件毎に実験又はシミュレーションにより決定する。）との関係を、理論最大吸入空気量ＱMAXにより無次元化して作成したものである。すなわち、この変換テーブルは、基準シリンダ吸入空気量及び理論最大吸入空気量の比（「第１の比」に相当する。）ＲＱ１＝ＱD／ＱMAXと、シリンダ吸入空気量及び理論最大吸入空気量の比（「第２の比」に相当する。）ＲＱ２＝Ｑｃｙｌ／ＱMAXとの間の一義的な関係を規定する。本実施形態では、筒内圧力Ｐ０として、吸気行程における筒内の状態変化が断熱膨張により行われるものとした場合に、吸気弁１０４の作動中心角ＩＶＣＮＴで得られる筒内圧力Ｐｃｔｒを採用している。この場合の筒内圧力Ｐ０（＝Ｐｃｔｒ）は、熱力学上の理論式により容易に算出することができる。また、（１ａ）及び（１ｂ）式において、吸気弁開期間における吸気ポート１０１ａの総開口面積をΣＡ（単位クランク角毎の開口面積Ａを積算して算出する。）とするとともに、開口面積Ａの積算間隔をΔθとし、吸入空気の比熱比、ガス定数及び温度をκ，Ｒａ，Ｔｍとしている。理論最大吸入空気量ＱMAXは、吸気の開始から終了までの行程容積を吸気弁上流における吸入空気の密度（又は圧力）及び温度で充填した場合に得られるシリンダ吸入空気量であり、上死点ＴＤＣにおけるシリンダ容積をＶＴＤＣとし、吸気弁閉時期ＩＶＣにおけるシリンダ容積をＶＩＶＣとして、下式（２）により与えられる。
【００２３】
ＱD＝ΣＡ×（Δθ／６・Ｎｅ）×（Ｐｍ／√（Ｒａ・Ｔｍ））×Ｘ  ・・・（１ａ）
Ｘ＝√｛（２κ／（κ―１））×（（Ｐ０／Ｐｍ）2/κ−（Ｐ０／Ｐｍ）(κ+1)/κ）｝  ・・・（１ｂ）
ＱMAX＝（Ｐｍ／（Ｒａ・Ｔｍ））×（ＶＩＶＣ−ＶＴＤＣ）  ・・・（２）
本実施形態では、筒内圧力Ｐ０として断熱膨張下での圧力Ｐｃｔｒを採用した関係上、エンジン１の運転領域全体で吸入空気の流れが理論的にチョークし、吸入空気が音速で筒内に流入することとなり、（１ｂ）式の圧力比Ｐ０／Ｐｍは、常に臨界圧力比（＝（２／（κ＋１））κ/(κ―1)＝ｃｏｎｓｔ）を示すこととなる。このため、この基準シリンダ吸入空気量ＱDを、特に「仮想ソニック吸入空気量」と呼ぶこととする。
【００２４】
ここで、本実施形態に係る目標吸気作動角ｔθｅｖｅｎｔの設定原理について、図１２を参照して説明する。
【００２５】
図１２は、吸気弁１０４及び排気弁１０８の作動特性（弁リフトＶＬＩＦＴｉ，ＶＬＩＦＴｅにより示す。）と、この作動特性により得られる筒内圧力Ｐｃｙｌ及び単位クランク角当たりのシリンダ吸入空気量（すなわち、流量）ＤＬＴＱとの関係を示している。
【００２６】
本実施形態では、仮想ソニック吸入空気量ＱD及び理論最大吸入空気量ＱMAXの算出期間ＰＲＤＱＤ，ＰＲＤＱＭＡＸを、夫々次のように定義する。すなわち、仮想ソニック吸入空気量ＱDの算出期間ＰＲＤＱＤは、筒内圧力Ｐｃｙｌが低下して吸気圧力Ｐｍに一致する点（以下「実効上死点」という。）ＴＤＣＲから、設定上の吸気弁閉時期（以下「設定閉時期」という。）ＩＶＣまでの期間とする。このため、（１ａ）式の総開口面積ΣＡとして、単位クランク角当たりの開口面積Ａを実効上死点ＴＤＣＲから設定閉時期ＩＶＣに渡り積算したものを採用している。他方、理論最大吸入空気量ＱMAXの算出期間ＰＲＤＱＭＡＸは、実効上死点ＴＤＣＲから、筒内で吸入空気の圧縮が実質的に開始される点（以下「実効閉時期」という。）ＩＶＣＲまでの期間とする。このため、（２）式のシリンダ容積ＶＴＤＣとして、実効上死点ＴＤＣＲにおけるシリンダ容積ＶＴＤＣＲを採用するとともに、シリンダ容積ＶＩＶＣとして、実効閉時期ＩＶＣＲにおけるシリンダ容積ＶＩＶＣＲを採用している。なお、実効上死点ＴＤＣＲ及び実効閉時期ＩＶＣＲは、制御上、オーバーラップ中心角ＯＶＬＣＮＴからのオフセット量ＴＤＣＯＦＳ、設定閉時期ＩＶＣからのオフセット量ＩＶＣＯＦＳにより表すこととする。
【００２７】
このように定められる仮想ソニック吸入空気量ＱD及び理論最大吸入空気量ＱMAXのもと、変換テーブルで定められるシリンダ吸入空気量Ｑｃｙｌは、吸気圧力Ｐｍ及び筒内圧力Ｐｃｙｌに基づいて与えられる、静的なシリンダ吸入空気量である。これに対し、実際のシリンダ吸入空気量は、この静的なシリンダ吸入空気量に吸気脈動に起因する変動分が重畳したものとなる。本実施形態では、この変動分を次に述べる特性により見積ることとし、目標吸気作動角ｔθｅｖｅｎｔの設定に際し、目標新気量ｔＱｃｙｌからこの変動分を排除する補正を施すこととしている。
【００２８】
シリンダ吸入空気量Ｑｃｙｌは、吸気圧力Ｐｍと筒内圧力Ｐｃｙｌとの比をＲＰ（＝Ｐｃｙｌ／Ｐｍ）として、下式（３）により与えられる。
【００２９】
Ｑｃｙｌ＝ΣＡ×（Ｐｍ／√（Ｒａ・Ｔｍ））×√｛（２κ／（κ−１））×（ＲＰ2/κ−ＲＰ(κ+1)/κ）｝×Δｔ  ・・・（３）
ここで、シリンダ吸入空気量Ｑｃｙｌの特性を、吸入空気の流れがチョークするものとした第１の領域とこれ以外の第２の領域とに分けて定義する。第１の領域において、この特性は、下式（４）により与えられる。これは、（３）式の圧力比ＲＰが臨界圧力比（＝ｃｏｎｓｔ）であることから明らかであり、（４）式は、シリンダ吸入空気量Ｑｃｙｌが吸気圧力Ｐｍに比例し、吸気温度Ｔｍの平方根の逆数に比例することを示している。他方、第２の領域において、この特性は、筒内の状態変化が準静的に進行するものとして、下式（５）により与えられる。これは、吸気弁閉時期ＩＶＣに筒内が吸気通路１０１内の密度及び温度で充填されたものとした場合の気体の状態方程式から明らかであり、（５）式は、シリンダ吸入空気量Ｑｃｙｌが吸気脈動分を加味した実際の吸気圧力Ｐｍ（＝Ｐｍａｖｅ＋ΔＰｍｉｖｃ）に比例し、実際の吸気温度Ｔｍ（＝Ｔｍａｖｅ＋ΔＴｍｉｖｃ）の逆数に比例することを示している。なお、（５）式において、Ｐｍａｖｅ，Ｔｍａｖｅは、平均又は代表の吸気圧力及び吸気温度を示し、ΔＰｍｉｖｃ，ΔＴｍｉｖｃは、吸気弁閉時期ＩＶＣにおける吸気圧力Ｐｍｉｖｃ及び吸気温度Ｔｍｉｖｃの平均吸気圧力Ｐｍａｖｅ及び平均吸気温度Ｔｍａｖｅに対する変化量を示している。また、本実施形態に関し、記号∝は、左辺の値が右辺の値に比例することを示すものとする。
【００３０】
ａ）ＤＬＴＱ≫０：
Ｑｃｙｌ∝Ｐｍ×（Ｔｍ）-1/2  ・・・（４）
ｂ）ＤＬＴＱ≒０：
Ｑｃｙｌ∝Ｐｃｙｌｉｖｃ×（Ｔｃｙｌｉｖｃ）-1
＝Ｐｍｉｖｃ×（Ｔｍｉｖｃ）-1
＝（Ｐｍａｖｅ＋ΔＰｍｉｖｃ）×（Ｔｍａｖｅ＋ΔＴｍｉｖｃ）-1  ・・・（５）
このように得られた２つの特性を内包するものとして、第１及び第２の領域を含む運転領域全体に渡る、最も確からしい１つの特性を近似により設定する。本実施形態では、流れの状態に応じて変化させ得る２つの係数をＫ１，Ｋ２として、この１つの特性を次式（６）により定義する。なお、係数Ｋ１，Ｋ２は、いずれも０以上、かつ１以下の値をとり、シリンダ吸入空気量Ｑｃｙｌと理論最大吸入空気量ＱMAXとの比に応じ、これが大きいときほど大きな値に設定される（図７のブロックＢ２０１ａ，Ｂ２０２ａ）。
【００３１】
Ｑｃｙｌ∝（Ｐｍａｖｅ＋Ｋ１・ΔＰｍｉｖｃ）×（Ｔｍａｖｅ＋Ｋ１・ΔＴｍｉｖｃ）-1/(2-K2)  ・・・（６）
吸気脈動分を加味した実際のシリンダ吸入空気量Ｑｃｙｌは、（６）式の特性をもとに、静的なシリンダ吸入空気量をＱｃｙｌ０とし、かつ吸気脈動分ΔＰｍｉｖｃ，ΔＴｍｉｖｃを加味した場合のものとこの吸気脈動分を０とした場合のものとの比をＲｐｕｌとして、次式（７）及び（８）により与えられる。
【００３２】
Ｑｃｙｌ＝Ｑｃｙｌ０×Ｒｐｕｌ  ・・・（７）
Ｒｐｕｌ＝｛（Ｐｍａｖｅ＋Ｋ１・ΔＰｍｉｖｃ）／Ｐｍａｖｅ｝×｛（Ｔｍａｖｅ＋Ｋ１・ΔＴｍｉｖｃ）／Ｔｍａｖｅ｝-1/(2-K2)  ・・・（８）
この（８）式から圧力補正項及び温度補正項を抽出し、夫々シリンダ吸入空気量の圧力補正係数ＰＲＡＴＥ、温度補正係数ＴＲＡＴＥとして設定する。
【００３３】
ＰＲＡＴＥ＝（Ｐｍａｖｅ＋Ｋ１・ΔＰｍｉｖｃ）／Ｐｍａｖｅ  ・・・（９）
ＴＲＡＴＥ＝｛（Ｔｍａｖｅ＋Ｋ１・ΔＴｍｉｖｃ）／Ｔｍａｖｅ｝-1/(2-K2)  ・・・（１０）
ＥＣＵ２０１は、設定した目標新気量ｔＱｃｙｌを圧力及び温度補正係数ＰＲＡＴＥ，ＴＲＡＴＥで除算して、この目標新気量を静的な成分ｔＱｃｙｌ０に変換するとともに、変換した目標新気量ｔＱｃｙｌ０によりブロックＢ１０１の変換テーブルを検索して、仮想ソニック吸入空気量ＱDを算出する。既述の通り、この仮想ソニック吸入空気量ＱDは、吸入空気の流れがチョークしているものとした場合に得られるシリンダ吸入空気量であるから、吸気ポート１０１ａの総開口面積ΣＡに比例する特性を持つ。このため、この仮想ソニック吸入空気量ＱDをもとに、目標とする総開口面積ｔＡｓｕｍｉｎを算出し、これに基づいて目標吸気作動角ｔθｅｖｅｎｔを設定することができる。
【００３４】
図４は、ＥＣＵ２０１のうち目標吸気作動角ｔθｅｖｅｎｔの設定に関する部分の構成を示している。なお、本実施形態では、目標吸気作動角ｔθｅｖｅｎｔの設定に関する演算において、吸気圧力Ｐｍ及び吸気温度Ｔｍとして、吸気通路１０１内のサイクル毎の平均圧力及び平均温度を採用している。
【００３５】
ＥＣＵ２０１は、エンジン１の運転条件に応じた目標トルク（アクセル開度ＡＰＯ等に基づいて算出する。）をもとに、図５に示す傾向のマップからの検索により目標新気量ｔＱｃｙｌを算出する。
【００３６】
目標吸気作動角設定部（図４）は、設定した目標新気量ｔＱｃｙｌに対して下式（１１）による補正を施して、この目標新気量ｔＱｃｙｌから吸気脈動分を排除し、静的な成分ｔＱｃｙｌ０に変換する。なお、（１１）式において、オーバーラップ期間に筒内から吸気通路１０１内へ吹き返すガスの量（すなわち、吹返ガス量）をＱIFBとしている。この吹返ガス量ＱIFBは、係数をＫ３として、吸気圧力Ｐｍ及び排気圧力Ｐｅに基づいて下式（１２）及び（１３）により算出する。なお、（１２）式において、吸気ポート１０１ａの開口面積ΣＡとして、吸気弁開時期ＩＶＯからオーバーラップ中心角ＯＶＬＣＮＴまでの前半開口面積ΣＡｉｖを採用している。オーバーラップ中心角ＯＶＬＣＮＴは、吸気弁１０４の弁リフトＶＬＩＦＴｉと排気弁１０８の弁リフトＶＬＩＦＴｅとの差が最も小さくなるクランク角をいい、本実施形態では、排気弁１０８の作動特性を一定としているため、吸気作動角及び作動中心角から割り出すことができる。また、係数Ｋ３は、エンジン回転数Ｎｅに比例する、１以上の値として設定する（図６）。排気圧力Ｐｅ及び排気温度Ｔｅとして、排気通路１０７内のサイクル毎の平均圧力及び平均温度を採用している。
【００３７】
ｔＱｃｙｌ０＝（ｔＱｃｙｌ／（ＰＲＡＴＥ×ＴＲＡＲＥ））−ＱIFB  ・・・（１１）
ＱIFB＝ΣＡｉｖ×（Δθ／（６・Ｎｅ））×Ｐｅ／（√（Ｒａ・Ｔｅ））×ＶIFB×Ｋ３  ・・・（１２）
ＶIFB＝√｛（２κ／（κ―１））×（（Ｐｍ／Ｐｅ）2/κ−（Ｐｍ／Ｐｅ）(κ+1)/κ）｝  ・・・（１３）
図７は、補正係数算出部の構成を示している。この部分により圧力補正係数ＰＲＡＴＥ及び温度補正係数ＰＲＡＴＥが算出され、目標吸気作動角設定部（図４）に提供される。補正係数算出部は、ブロックＢ２０１において、後述する脈動変化量算出部により算出された脈動圧力変化量ΔＰｍｉｖｃをもとに、（９）式により圧力補正係数ＰＲＡＴＥを算出する。係数Ｋ１は、０以上、かつ１以下の範囲内で、比Ｑｃｙｌ／ＱMAXが大きくなり、筒内の状態変化が準静的なものに近付くのに従い２次的に増大する（ブロックＢ２０１ａ）。また、ブロックＢ２０２において、後述する脈動変化量算出部により算出された脈動温度変化量ΔＴｍｉｖｃをもとに、（１０）式により温度補正係数ＴＲＡＴＥを算出する。係数Ｋ２も、係数ｋ１と同様に０以上、かつ１以下の範囲内で、比Ｑｃｙｌ／ＱMAXが大きくなるのに従い２次的に増大する（ブロックＢ２０２ａ）。
【００３８】
図８は、脈動変化量算出部の構成を示している。前回の演算実行時に設定した目標吸気作動角ｔθｅｖｅｎｔZ（＝ＩＶＣ−ＩＶＯ）及びエンジン回転数Ｎｅに基づいて脈動圧力変化量の基本値ΔＰｍｉｖｃ０をマップＢ３０１からの検索により算出し、この基本値ΔＰｍｉｖｃ０に対して実際の吸気圧力Ｐｍ（＝Ｐｍａｖｅ）による補正を施して、脈動圧力変化量ΔＰｍｉｖｃを算出する。同様にして、目標吸気作動角ｔθｅｖｅｎｔZ及びエンジン回転数Ｎｅに基づいて脈動温度変化量の基本値ΔＴｍｉｖｃ０をマップＢ３０２からの検索により算出し、この基本値ΔＴｍｉｖｃ０に対して実際の吸気温度Ｔｍ（＝Ｔｍａｖｅ）による補正を施して、脈動温度変化量ΔＴｍｉｖｃを算出する。なお、（１４）及び（１５）式において、大気圧をＰａｔｍとしている。
【００３９】
ΔＰｍｉｖｃ＝ΔＰｍｉｖｃ０×Ｐｍ／Ｐａｔｍ  ・・・（１４）
ΔＴｍｉｖｃ＝ΔＴｍｉｖｃ０×Ｐｍ／Ｐａｔｍ  ・・・（１５）
図９は、理論最大吸入空気量算出部の構成を示している。この部分により理論最大吸入空気量ＱMAXが算出され、目標吸気作動角設定部（図４）に提供される。理論最大吸入空気量算出部は、吸気弁１０４の設定閉時期ＩＶＣ及びオーバーラップ中心角ＯＶＬＣＮＴ等をもとに、（２）式により理論最大吸入空気量ＱMAXを算出する。なお、既述の通り、理論最大吸入空気量ＱMAXの算出期間ＰＲＤＱＭＡＸは、実効上死点ＴＤＣＲから実効閉時期ＩＶＣＲまでの期間に設定している。実効上死点ＴＤＣＲは、オーバーラップ中心角ＯＶＬＣＮＴから上死点のオフセット量ＴＤＣＯＦＳだけ遅角させた時期として算出する。このオフセット量ＴＤＣＯＦＳは、図１０に示す傾向のマップからの検索により推定する。この図１０のマップにおいて、オフセット量ＴＤＣＯＦＳは、エンジン回転数Ｎｅが高く、かつオーバーラップ期間中の開口面積（たとえば、前半開口面積ΣＡｉｖ）が小さいときほど大きな値に設定されている。また、実効閉時期ＩＶＣＲは、設定閉時期ＩＶＣから吸気弁閉時期のオフセット量ＩＶＣＯＦＳだけ進角させた時期として算出され、このオフセット量ＩＶＣＯＦＳは、図１１に示す傾向のマップからの検索により推定する。この図１１のマップにおいて、オフセット量ＩＶＣＯＦＳは、エンジン回転数Ｎｅが高く、かつ吸気弁１０４の弁リフトＶＬＩＦＴｉ（たとえば、最大弁リフト）が小さいときほど大きな値に設定されている。
【００４０】
ＥＣＵ２０１は、静的な成分に変換した目標新気量ｔＱｃｙｌ０を理論最大吸入空気量ＱMAXで除算して、第２の比ＲＱ２（＝ｔＱｃｙｌ０／ＱMAX）を算出するとともに、算出した第２の比ＲＱ２により変換テーブルを検索して、第１の比ＲＱ１（＝ＱD／ＱMAX）を算出する。また、算出した第１の比ＲＱ１に理論最大吸入空気量ＱMAXを乗算して、仮想ソニック吸入空気量ｔＱDを算出し、これに基づいて目標吸気作動角ｔθｅｖｅｎｔを算出する。
【００４１】
すなわち、ＥＣＵ２０１は、臨界圧力比（＝（２／（κ＋１））κ/(κ―1)）のもとで得られる単位開口面積当たりのシリンダ吸入空気量ｑｓｏｎｉｃ＃に対し、吸入空気の状態に応じた変数（＝Ｐｍ／√（Ｒａ＃・Ｔｍ））を乗算するとともに、これで仮想ソニック吸入空気量ｔＱDを除算して、目標総開口面積ｔＡｓｕｍｉｎを算出する。更に、算出した目標総開口面積ｔＡｓｕｍｉｎをクランク角単位に換算するとともに、換算後の目標総開口面積ｔＡｓｕｍｉｎに基づいて目標吸気作動角ｔθｅｖｅｎｔを算出及び設定する。すなわち、ＥＣＵ２０１は、目標総開口面積ｔＡｓｕｍｉｎによりブロックＢ１０３のテーブルを検索して、目標吸気作動角ｔθｅｖｅｎｔを算出する。このテーブルにおいて、吸気作動角θｅｖｅｎｔは、総開口面積Ａｓｕｍｉｎが大きいときほど大きな値に設定されている。
【００４２】
設定した目標吸気作動角ｔθｅｖｅｎｔは、吸気弁１０４の目標開時期ｔＩＶＯ及び目標閉時期ｔＩＶＣの演算に用いられ、吸気動弁装置１０５の制御に反映される。
【００４３】
本実施形態によれば、次のような効果を得ることができる。
【００４４】
すなわち、本実施形態では、第１及び第２の比ＱD／ＱMAX，Ｑｃｙｌ／ＱMAXを定義するとともに、これらの比の一義的な関係を変換テーブル（ブロックＢ１０１）として作成及び設定し、この変換テーブルにより仮想ソニック吸入空気量ＱDを算出し、算出した仮想ソニック吸入空気量ＱDに基づいて目標吸気作動角ｔθｅｖｅｎｔを設定することとした。このため、エンジン１の運転領域全体において、マップを多用することなくこの変換テーブルにより目標吸気作動角ｔθｅｖｅｎｔを設定することが可能となり、演算負荷を軽減して目標新気量ｔＱｃｙｌを達成することができる。
【００４５】
また、変換テーブルの作成に際し、シリンダ吸入空気量Ｑｃｙｌを吸気圧力Ｐｍと筒内圧力Ｐｃｙｌとの比に応じた静的なものとして扱うことで、この変換テーブルの作成を容易、かつ正確なものとすることができる。また、実際の運転に際し、目標新気量ｔＱｃｙｌから吸気脈動分を排除する補正をすることとし、この補正に関し、吸入空気の流れがチョークする第１の領域とこれ以外の第２の領域とを定め、領域毎に補正の方法を異ならせたことで、吸気脈動分を的確、かつ簡易に排除することができる。
【００４６】
更に、目標新気量ｔＱｃｙｌから吹返ガス量ＱIFBを減算することで、変換テーブルにより的確な仮想ソニック吸入空気量ＱDを算出することができる。
【００４７】
なお、吸入空気量の制御は、運転領域全体に渡り弁作動特性の制御によることとしてもよいが、弁作動特性による制御を過渡運転時のみで行い、定常運転時には、スロットル開度の制御によることとしてもよい。
【００４８】
本発明に係るシリンダ吸入空気量制御装置は、直噴ガソリンエンジンに採用することもできる。
【図面の簡単な説明】
【００４９】
【図１】本発明の一実施形態に係るエンジンの構成
【図２】同上エンジンの吸気動弁装置の構成
【図３】同上吸気動弁装置の作動角変更機構の構成
【図４】目標吸気作動角設定部の構成
【図５】目標新気量の算出マップ
【図６】係数Ｋ３の設定マップ
【図７】補正係数算出部の構成
【図８】脈動変化量算出部の構成
【図９】理論最大吸入空気量算出部の構成
【図１０】上死点のオフセット量ＴＤＣＯＦＳの算出マップ
【図１１】吸気弁閉時期のオフセット量ＩＶＣＯＦＳの算出マップ
【図１２】弁作動特性、筒内圧力及び単位クランク角当たりのシリンダ吸入空気量の関係
【符号の説明】
【００５０】
１…エンジン、１０１…吸気通路、１０２…スロットル弁、１０３…インジェクタ、１０４…吸気弁、１０５…吸気動弁装置、１０６…点火プラグ、１０７…排気通路、１０８…排気弁、１５１…駆動軸、１５２…揺動カム、１５３…偏心駆動カム、１５４…リング状リンク、１５５…制御軸、１５６…偏心制御カム、１５７…ロッカーアーム、１５８…ロッド状リンク、１６１…電磁アクチュエータ、１６２…ギア列、２０１…エンジンコントローラ、２１１…アクセルセンサ、２１２…クランク角センサ、２１３…吸気圧力センサ、２１４…吸気温度センサ、２１５…排気圧力センサ、Ａ…吸気動弁装置の作動角変更機構、Ｂ…吸気動弁装置の中心角変更機構。
【図１】
2006105098.tif 000003
【図２】
2006105098.tif 000004
【図３】
2006105098.tif 000005
【図４】
2006105098.tif 000006
【図５】
2006105098.tif 000007
【図６】
2006105098.tif 000008
【図７】
2006105098.tif 000009
【図８】
2006105098.tif 000010
【図９】
2006105098.tif 000011
【図１０】
2006105098.tif 000012
【図１１】
2006105098.tif 000013
【図１２】
2006105098.tif 000014