(57)【要約】    【課題】内部ＥＧＲ量の推定に際し、オーバーラップ期間中の吹抜ガス量を簡易に推定する。    【解決手段】吹抜ガス量は、オーバーラップ期間中の吸気圧力及び排気圧力、並びに吸気弁及び排気弁により形成される有効開口面積ＡＳＵＭＯＬをもとに、算出する。オーバーラップ期間を所定の期間ＤＣＡ毎に分割し、分割した各区間について弁リフト量ＶＬＩＦＴｎに応じた区間開口面積を算出する。算出した区間開口面積は、区間毎に対応させた配列として記憶する。吹抜ガス量の算出では、この配列を参照して、各区間の実質的な開口面積を算出するとともに、算出した区間開口面積を積算して、有効開口面積ＡＳＵＭＯＬを算出する。【選択図】  図７
2005307803.tif 000002
【特許請求の範囲】
【請求項１】
吸気弁開期間と排気弁開期間とのオーバーラップ期間に吸気弁及び排気弁により形成される、ポート部の実質的な開口面積を有効開口面積として算出する装置であって、
オーバーラップ期間を含む所定の期間を複数の区間に分割し、
分割した各区間につき、予め定められた弁作動特性値をもとに、吸気弁及び排気弁により形成される区間開口面積を算出し、
算出した区間開口面積をオーバーラップ期間に亘り積算して、有効開口面積を算出するエンジンの弁開口面積算出装置。
【請求項２】
吸気弁開期間と排気弁開期間とのオーバーラップ期間中に排気側と吸気側との間で吹き抜ける排気の量を吹抜ガス量として算出する吹抜ガス量算出手段と、
算出した吹抜ガス量をもとに、エンジンの内部ＥＧＲ量を算出する内部ＥＧＲ量算出手段と、を含んで構成され、
吹抜ガス量算出手段は、オーバーラップ期間における排気圧力を検出する排気圧力検出手段と、オーバーラップ期間における吸気圧力を検出する吸気圧力検出手段と、吸気弁及び排気弁により形成される、吹抜ガスの通路の実質的な開口面積を有効開口面積として算出する有効開口面積算出手段と、を含んで構成され、少なくとも検出した排気圧力及び吸気圧力、並びに算出した有効開口面積をもとに、吹抜ガス量を算出し、
有効開口面積算出手段は、オーバーラップ期間を含む所定の期間を複数の区間に分割するとともに、各区間につき、予め定められた弁作動特性値をもとに、吸気弁及び排気弁により形成される実質的な区間開口面積を算出し、算出した区間開口面積をオーバーラップ期間に亘り積算して、有効開口面積を算出するエンジンの内部ＥＧＲ量推定装置。
【請求項３】
排気弁閉時期に筒内に残る排気の量を残留ガス量として算出する残留ガス量算出手段を更に含んで構成され、
内部ＥＧＲ量算出手段は、吹抜ガス量算出手段及び残留ガス量算出手段により算出された各ガス量をもとに、エンジンの内部ＥＧＲ量を算出する請求項２に記載のエンジンの内部ＥＧＲ量推定装置。
【請求項４】
残留ガス量算出手段は、排気弁閉時期における筒内圧力を検出する手段と、排気弁閉時期における筒内温度を検出する手段と、排気弁閉時期における実質的なシリンダ容積を検出する手段と、混合気空燃比に応じた排気のガス定数を算出する手段と、を含んで構成され、検出した筒内圧力、筒内温度及びシリンダ容積、並びに算出したガス定数をもとに、残留ガス量を算出する請求項３に記載のエンジンの内部ＥＧＲ量推定装置。
【請求項５】
内部ＥＧＲ量算出手段は、算出された吹抜ガス量及び残留ガス量の和として内部ＥＧＲ量を算出する請求項３又は４に記載のエンジンの内部ＥＧＲ量推定装置。
【請求項６】
有効開口面積算出手段は、区間毎に対応させて、吸気弁により形成される区間開口面積の第１の配列と、排気弁により形成される区間開口面積の第２の配列とを作成するとともに、各区間につき、各配列の区間開口面積のうち小さい方をその区間に関する実質的な区間開口面積として選択し、選択した区間開口面積をオーバーラップ期間に亘り積算して、有効開口面積を算出する請求項２〜５のいずれかに記載のエンジンの内部ＥＧＲ量推定装置。
【請求項７】
吸気弁開期間及び排気弁開期間の少なくとも一方がクランク角に関して進遅可能なエンジンに設けられ、
有効開口面積算出手段は、弁開期間の変更により排気弁閉時期が吸気弁開時期に対して変化したときに、第１の配列に対し、第２の配列をこの変化に応じた量だけずらすとともに、各区間につき、各配列の区間開口面積のうち小さい方をその区間に関する実質的な区間開口面積として選択する請求項６に記載のエンジンの内部ＥＧＲ量推定装置。
【請求項８】
有効開口面積算出手段は、所定の期間をクランク角で１°毎に分割する請求項２〜７のいずれかに記載のエンジンの内部ＥＧＲ量推定装置。
【請求項９】
吹抜ガス量算出手段は、排気温度を検出する手段を更に含んで構成され、検出した排気圧力、吸気圧力及び排気温度、並びに算出した有効開口面積をもとに、吹抜ガス量を算出する請求項２〜８のいずれかに記載のエンジンの内部ＥＧＲ量推定装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、エンジンの弁開口面積算出装置及び内部ＥＧＲ量推定装置に関し、詳細には、内部ＥＧＲ量のうち、吸気弁開期間と排気弁開期間とのオーバーラップ期間中に排気側と吸気側との間で吹き抜ける排気の量（すなわち、吹抜ガス量）を簡易、かつ正確に算出するための技術に関する。
【背景技術】
【０００２】
従来、エンジンでは、燃焼温度の上昇による窒素酸化物の発生を抑制するため、排気の一部を筒内に戻す排気還流（以下「ＥＧＲ」という。）が行われている。
ＥＧＲには、排気管と吸気管との間に接続されたＥＧＲ管を介して行う外部ＥＧＲと、このＥＧＲ管を介さずに行う内部ＥＧＲとがある。これらのうち、内部ＥＧＲによる還流ガス量を推定する装置として、次のものが知られている（特許文献１）。吸気弁開期間と排気弁開期間とがオーバーラップしない作動条件のもと、エンジン回転数及び排気弁閉時期に基づいて内部ＥＧＲ量の基本値を算出するとともに、オーバーラップするときは、オーバーラップ量及び吸気圧力等を考慮して算出したオーバーラップ分の補正値を加算して、内部ＥＧＲ量を算出するものである。
【特許文献１】特開２００１−２２１１０５号公報（段落番号００５９）
【発明の開示】
【発明が解決しようとする課題】
【０００３】
しかしながら、基本値に対し、オーバーラップ分の補正値を加算することによる上記の装置には、次のような問題がある。補正値の算出に際し、吸気圧力等を考慮することとしているが、エンジン回転数や混合気空燃比等の影響をも考慮しようとすれば、与えられたオーバーラップ量に対して加算すべき補正値を一義的に決定することができず、多くの適合工数が必要となることである。
【０００４】
本出願人は、次のような内部ＥＧＲ量推定装置について検討している。内部ＥＧＲガスを、オーバーラップ期間中に排気側から吸気側に吹き抜ける排気（以下「吹抜ガス」という。）と、排気弁閉時期を過ぎても筒内に残る排気（以下「残留ガス」という。）とに分け、各ガス量の和として内部ＥＧＲ量を算出するものである。前者の吹抜ガス量は、オーバーラップ期間中に吸気側又は排気側のポート部に形成される絞りを理想ノズルと見立て、エネルギー保存則を適用することにより算出し、他方、後者の残留ガス量は、排気弁閉時期における筒内温度及び圧力等を推定するとともに、各推定値を気体の状態方程式に代入することにより算出する。吹抜ガス量の算出に際し、ノズルの有効開口面積は、吸気ポートの瞬時開口面積と排気ポートの瞬時開口面積とのうち小さい方を実質的な瞬時開口面積とし、各時点における瞬時開口面積を積算した積算面積として算出する。オーバーラップ量と積算面積との関係をテーブルデータとして予め作成及び記憶しておき、運転時には、検出した実際のオーバーラップ量によりこのテーブルを検索して、積算面積、すなわち、有効開口面積を算出する。
【０００５】
しかしながら、本出願人による上記の装置には、有効開口面積の算出に関して次のような問題がある。オーバーラップ量と積算面積との関係を予め特定し、これをテーブルデータとして記憶させることとしているので、バルブクリアランス等の変化（吸気弁と排気弁とでこのクリアランスが異なる場合もある。）により実質的な瞬時開口面積が変化した場合は、算出される有効開口面積に対してこの変化の影響を正確に反映させるのが困難であり、正確さに欠け、吹抜ガス量を正確に算出することができないことである。
【０００６】
本発明は、有効開口面積を正確に算出し、内部ＥＧＲ量の推定では、オーバーラップ期間中の吹抜ガス量を正確に算出することを目的とする。
【課題を解決するための手段】
【０００７】
本発明は、エンジンの弁開口面積算出装置及び内部ＥＧＲ量推定装置を提供する。本発明に係る弁開口面積算出装置は、吸気弁開期間と排気弁開期間とのオーバーラップ期間に吸気弁及び排気弁により形成される、ポート部の実質的な開口面積を有効開口面積として算出する装置である。オーバーラップ期間を含む所定の期間を複数の区間に分割するとともに、各区間につき、予め定められた弁作動特性値をもとに、吸気弁及び排気弁により形成される区間開口面積を算出し、算出した区間開口面積をオーバーラップ期間に亘り積算して、有効開口面積を算出する。また、本発明に係る内部ＥＧＲ量推定装置は、吸気弁開期間と排気弁開期間とのオーバーラップ期間中に排気側と吸気側との間で吹き抜ける排気の量（以下「吹抜ガス量」という。）を算出し、算出した吹抜ガス量をもとに、エンジンの内部ＥＧＲ量を算出する。吹抜ガス量の算出では、オーバーラップ期間中の排気圧力及び吸気圧力を検出するとともに、吸気弁及び排気弁により形成される、吹抜ガスの通路の実質的な開口面積を有効開口面積として算出し、少なくとも検出した排気圧力及び吸気圧力、並びに算出した有効開口面積をもとに、吹抜ガス量を算出する。有効開口面積の算出では、オーバーラップ期間を含む所定の期間を複数の区間に分割するとともに、各区間につき、予め定められた弁作動特性値をもとに、吸気弁及び排気弁により形成される実質的な区間開口面積を算出し、算出した区間開口面積をオーバーラップ期間に亘り積算して、有効開口面積を算出する。有効開口面積の算出では、区間毎に対応させて、吸気弁により形成される区間開口面積の第１の配列と、排気弁により形成される区間開口面積の第２の配列とを作成するとともに、各区間につき、各配列の区間開口面積のうち小さい方をその区間に関する実質的な区間開口面積として選択し、選択した区間開口面積をオーバーラップ期間に亘り積算して、有効開口面積を算出するとよい。
【発明の効果】
【０００８】
本発明によれば、オーバーラップ期間を含む所定の期間を複数の区間に分割するとともに、各区間について実質的な区間開口面積を算出し、算出した区間開口面積を積算して、有効開口面積を算出することとしたので、バルブクリアランス等が変化したとしても、算出される有効開口面積に対し、その影響を正確に反映させることができる。また、区間毎の開口面積を積算することで、演算負荷を抑えることができる。更に、吸気弁及び排気弁の区間開口面積を配列化することで、バルブタイミングが変更された場合に、これに簡単に対応することができる。
【発明を実施するための最良の形態】
【０００９】
以下に図面を参照して、本発明の実施の形態について説明する。
図１は、本発明の一実施形態に係るエンジン１の構成を示している。
吸気通路１１の導入部には、エアクリーナ１２が取り付けられており、エアクリーナ１２により吸入空気中の粉塵等が除去される。吸気通路１１において、エアクリーナ１２の下流には、電子制御式のスロットル弁１３が設置されている。スロットル弁１３の下流には、サージタンク１４が取り付けられており、サージタンク１４にブランチ１５が取り付けられ、吸気マニホールドが構成されている。サージタンク１４内の吸入空気は、ブランチ１５及びシリンダヘッドに形成された吸気ポート１６を介して筒内に流入する。各気筒の吸気ポート１６には、燃料供給用のインジェクタ１７が設置されている。
【００１０】
エンジン本体において、燃焼室１８は、シリンダヘッド及びピストン１９により挟まれた空間として形成される。燃焼室１８は、気筒中心軸を基準とした一側で吸気ポート１６と接続しており、吸気ポート１６は、吸気弁２０により開放及び遮断される。吸気弁２０は、吸気カム２１により駆動される。また、燃焼室１８は、吸気ポート１６とは反対の一側で排気ポート２２と接続しており、排気ポート２２は、排気弁２３により開放及び遮断される。排気弁２３は、排気カム２４により駆動される。吸気カム２１に対して吸気側可変動弁装置２５が、排気カム２４に対して排気側可変動弁装置２６が設けられており、これらの可変動弁装置２５，２６により吸気カム２１又は排気カム２４の各カムシャフトに対する位相を変化させることで、吸気弁２０又は排気弁２３の作動特性を変化させ得るように構成されている。可変動弁装置２５，２６は、油圧型及びソレノイド型等のいかなる形態のものを採用してもよいが、本実施形態では、吸気弁２０又は排気弁２３の開閉時期（すなわち、バルブタイミング）を変化させることで、吸気弁開期間と排気弁開期間とのオーバーラップ期間（以下、単に「オーバーラップ期間」という。）を変化させ得るものを採用している。シリンダヘッドには、燃焼室１８の上部略中央に臨ませて点火プラグ２７が設置されている。
【００１１】
排気通路２８には、排気マニホールドの直後に第１の触媒コンバータ２９が介装されるとともに、その下流に第２の触媒コンバータ３０が介装されている。排気ポート２２に流出した排気は、これらの触媒コンバータ２９，３０及びマフラー３１を通過して、大気中に放出される。
インジェクタ１７、点火プラグ２７及び各可変動弁装置２５，２６の動作は、エンジンコントロールユニット（以下「ＥＣＵ」という。）４１により制御される。ＥＣＵ４１には、エアフローメータ５１からの吸入空気量検出信号、圧力センサ５２からの吸気圧力検出信号、温度センサ５３からの冷却水温度検出信号、クランク角センサ５４からの単位クランク角及び基準クランク角検出信号（ＥＣＵ４１は、これをもとに、エンジン回転数Ｎｅを算出する。）、圧力センサ５５からの排気圧力検出信号、温度センサ５６からの排気温度検出信号、酸素センサ５７からの空燃比検出信号、アクセルセンサ５８からのアクセル開度検出信号、及びカム角センサ５９，６０からのカム角検出信号（これをもとに、カムとカムシャフトとの実際の位相差を検出可能である。）が入力される。ＥＣＵ４１は、入力した各信号をもとに、上記の各デバイスの制御量を設定する。
【００１２】
本実施形態では、ＥＣＵ４１がエンジン１の内部ＥＧＲ量推定装置としての機能を備えている。
次に、ＥＣＵ４１による内部ＥＧＲ量ＭＲＥＳの推定について説明する。
本実施形態では、内部ＥＧＲ量ＭＲＥＳを、オーバーラップ期間中に排気側から吸気側に吹き抜ける排気の量である吹抜ガス量ＭＲＥＳＯＬと、排気弁閉時期を過ぎても筒内に残る排気の量である残留ガス量ＭＲＥＳＣＹＬとに分け、算出した各ガス量を加算することにより算出する。
【００１３】
ＭＲＥＳ＝ＭＲＥＳＯＬ＋ＭＲＥＳＣＹＬ  ・・・（１）
図２は、内部ＥＧＲ量推定ルーチンのフローチャートである。このルーチンは、一サイクル毎に繰り返され、クランク角センサ５４が出力する基準クランク角信号と同期させて実行される。
Ｓ１０１では、エンジン回転数Ｎｅ、吸気圧力Ｐｉｎ、排気圧力Ｐｅｘ及び排気温度Ｔｅｘ等、各種の運転状態を読み込む。
【００１４】
Ｓ１０２では、次のようにして残留ガス量ＭＲＥＳＣＹＬを算出する。排気弁閉時期にシリンダヘッドとピストン１９とにより画成される空間の容積（以下「シリンダ容積」という。）ＶＥＶＣを算出するとともに、目標燃焼当量比ＴＦＢＹＡに応じた排気のガス定数Ｒｅｘを算出する。算出したＶＥＶＣ及びＲｅｘと、排気弁閉時期における筒内温度及び圧力ＴＥＶＣ，ＰＥＶＣをもとに、次式により残留ガス量ＭＲＥＳＣＹＬを算出する。筒内温度ＴＥＶＣ及び筒内圧力ＰＥＶＣは、圧力センサ５５及び温度センサ５６からの信号に基づいて夫々推定することができ、排気弁閉時期は、排気側カム角センサ６０からの信号に基づいて検出する。
【００１５】
ＭＲＥＳＣＹＬ＝（ＰＥＶＣ×ＶＥＶＣ）／（Ｒｅｘ×ＴＥＶＣ）  ・・・（２）
Ｓ１０３では、有効開口面積ＡＳＵＭＯＬを算出する。有効開口面積ＡＳＵＭＯＬは、オーバーラップ期間中に吸気弁２０及び排気弁２３により形成される実質的な瞬時開口面積を積算したものに相当し、後述する有効開口面積演算ルーチンにより算出される。
Ｓ１０４では、排気の比熱比ＳＨＥＡＴＲを算出する。比熱比ＳＨＥＡＴＲは、目標燃焼当量比ＴＦＢＹＡ及び排気温度Ｔｅｘに応じて各比熱比を割り付けたマップ（図３）を検索して算出する。比熱比ＳＨＥＡＴＲは、排気温度Ｔｅｘを一定としたときに、理論空燃比相当下で最も小さく、目標燃焼当量比ＴＦＢＹＡがこれよりも小さく又は大きくなるほど、大きな値として算出される。また、目標燃焼当量比ＴＦＢＹＡを一定としたときに、リーン側及びリッチ側の領域の双方において、排気温度Ｔｅｘが高くなるほど、小さな値として算出される。
【００１６】
Ｓ１０５では、過給状態にあるか否かを判定する。過給状態の判定は、吸気圧力Ｐｉｎと排気圧力Ｐｅｘとの比ＰＩＮＢＹＥＸ（＝Ｐｉｎ／Ｐｅｘ）を算出するとともに、算出した圧力比ＰＩＮＢＹＥＸが１よりも大きいか否かを判定することにより行う。圧力比ＰＩＮＢＹＥＸが１よりも大きく、過給状態にあるときは、Ｓ１０６へ、これが１以下であり、過給状態にないときは、Ｓ１０７へ進む。
【００１７】
Ｓ１０６では、過給により吸気側から排気側に吹き抜ける混合気の流れにチョークが発生しているか否かを判定する。ここでのチョークの判定は、混合気の比熱比ＭＩＸＳＨＲを算出するとともに、比熱比ＳＨＥＡＴＲを算出した比熱比ＭＩＸＳＨＲに置き換えて、下式により第１及び第２のチョーク判定値ＳＬＣＨＯＫＥＨ，ＳＬＣＨＯＫＥＬを算出し、圧力比ＰＩＮＢＹＥＸがこれらのチョーク判定値ＳＬＣＨＯＫＥＨ，ＳＬＣＨＯＫＥＬを上限及び下限とする所定の範囲内にあるか否かを判定することにより行う。圧力比ＰＩＮＢＹＥＸが所定の範囲内になく、チョークが発生しているときは、Ｓ１０８へ、これが所定の範囲内にあり、チョークが発生していないときは、Ｓ１０９へ進む。なお、混合気の比熱比ＭＩＸＳＨＲは、目標燃焼当量比ＴＦＢＹＡに応じて各比熱比を割り付けたテーブル（図４）を検索して算出する。比熱比ＭＩＸＳＨＲは、目標燃焼当量比ＴＦＢＹＡが大きいときほど、小さな値として算出される。
【００１８】
ＳＬＣＨＯＫＥＨ＝（２／（ＳＨＥＡＴＲ＋１））＾（−ＳＨＥＡＴＲ／（ＳＨＥＡＴＲ−１））  ・・・（３ａ）
ＳＬＣＨＯＫＥＬ＝（２／（ＳＨＥＡＴＲ＋１））＾（ＳＨＥＡＴＲ／（ＳＨＥＡＴＲ−１））  ・・・（３ｂ）
Ｓ１０７では、排気側から吸気側に吹き抜ける排気の流れにチョークが発生しているか否かを判定する。ここでのチョークの判定は、（３ａ）及び（３ｂ）式により第１及び第２のチョーク判定値ＳＬＣＨＯＫＥＨ，ＳＬＣＨＯＫＥＬを算出するとともに、圧力比ＰＩＮＢＹＥＸと算出した各チョーク判定値ＳＬＣＨＯＫＥＨ，ＳＬＣＨＯＫＥＬとを比較することにより行う。圧力比ＰＩＮＢＹＥＸが所定の範囲内になく、チョークが発生しているときは、Ｓ１１０へ、これが所定の範囲内にあり、チョークが発生していないときは、Ｓ１１１へ進む。
【００１９】
Ｓ１０８では、次式により過給チョーク時吹抜ガス流量Ｑａを算出し、これを吹抜ガス流量ＭＲＥＳＯＬｔｍｐとする。なお、吸気温度をＴｉｎとし、吸入空気のガス定数をＲｉｎとする。
ＭＲＳＯＬＰＣ＝√｛（ＭＩＸＳＨＲ×（２／（ＭＩＸＳＨＲ＋１））＾（（ＭＩＸＳＨＲ＋１）／（ＭＩＸＳＨＲ−１）））｝  ・・・（４ａ）
Ｑａ＝ＭＲＥＳＯＬｔｍｐ＝−√｛１／（Ｒｉｎ×Ｔｉｎ）｝×Ｐｉｎ×ＭＲＳＯＬＰＣ  ・・・（４ｂ）
Ｓ１０９では、次式により過給非チョーク時吹抜ガス流量Ｑｂを算出し、これを吹抜ガス流量ＭＲＥＳＯＬｔｍｐとする。
【００２０】
ＭＲＳＯＬＰＴ＝√｛（ＭＩＸＳＨＲ／（ＭＩＸＳＨＲ−１）×（ＰＩＮＢＹＥＸ＾（−２／ＭＩＸＳＨＲ）−ＰＩＮＢＹＥＸ＾（−（ＭＩＸＳＨＲ＋１）／ＭＩＸＳＨＲ）））｝  ・・・（５ａ）
Ｑｂ＝ＭＲＥＳＯＬｔｍｐ＝−１．４×√｛１／（Ｒｉｎ×Ｔｉｎ）｝×Ｐｉｎ×ＭＲＳＯＬＰＴ  ・・・（５ｂ）
Ｓ１１０では、次式により無過給チョーク時吹抜ガス流量Ｑｃを算出し、これを吹抜ガス流量ＭＲＥＳＯＬｔｍｐとする。
【００２１】
ＭＲＳＯＬＤ＝√｛１／（Ｒｅｘ×Ｔｅｘ）｝  ・・・（６ａ）
ＭＲＳＯＬＰＣ＝√｛（ＳＨＥＡＴＲ×（２／（ＳＨＥＡＴＲ＋１））＾（（ＳＨＥＡＴＲ＋１）／（ＳＨＥＡＴＲ−１）））｝  ・・・（６ｂ）
Ｑｃ＝ＭＲＥＳＯＬｔｍｐ＝Ｐｅｘ×ＭＲＳＯＬＤ×ＭＲＳＯＬＰＣ  ・・・（６ｃ）
Ｓ１１１では、次式により通常時吹抜ガス流量Ｑｄを算出し、これを吹抜ガス流量ＭＲＥＳＯＬｔｍｐとする。
【００２２】
ＭＲＳＯＬＤ＝√｛１／（Ｒｅｘ×Ｔｅｘ）｝  ・・・（７ａ）
ＭＲＳＯＬＰ＝√｛（ＳＨＥＡＴＲ／（ＳＨＥＡＴＲ−１）×（ＰＩＮＢＹＥＸ＾（２／ＳＨＥＡＴＲ）−ＰＩＮＢＹＥＸ＾（（ＳＨＥＡＴＲ＋１）／ＳＨＥＡＴＲ）））｝  ・・・（７ｂ）
Ｑｄ＝ＭＲＥＳＯＬｔｍｐ＝１．４×Ｐｅｘ×ＭＲＳＯＬＤ×ＭＲＳＯＬＰ  ・・・（７ｃ）
Ｓ１１２では、算出した吹抜ガス流量ＭＲＥＳＯＬｔｍｐをもとに、次式により吹抜ガス量ＭＲＥＳＯＬを算出する。
【００２３】
ＭＲＥＳＯＬ＝（ＭＲＥＳＯＬｔｍｐ×６０×ＡＳＵＭＯＬ）／（Ｎｅ×３６０）  ・・・（８）
Ｓ１１３では、算出した残留ガス量ＭＲＥＳＣＹＬ及び吹抜ガス量ＭＲＥＳＯＬをもとに、式（１）により内部ＥＧＲ量ＭＲＥＳを算出する。
ＥＣＵ４１は、算出した内部ＥＧＲ量ＭＲＥＳをもとに、次式により内部ＥＧＲ率ＭＲＥＳＦＲを算出するとともに、算出した内部ＥＧＲ率ＭＲＥＳＦＲをインジェクタ１７による燃料噴射制御や、点火プラグ２７による点火制御に反映させる。なお、下式において、エアフローメータ５１により検出される吸入空気量をＭＡＣＹＬとし、目標燃空比をＴＦＢＹＡとする。
【００２４】
ＭＲＥＳＦＲ＝ＭＲＥＳ／（ＭＲＥＳ＋ＭＡＣＹＬ×（１＋ＴＦＢＹＡ／１４．７））  ・・・（９）
次に、有効開口面積ＡＳＵＭＯＬの演算について説明する。
図７は、クランク角ＣＡと、弁作動特性値ＣＡＭＰＦ及び弁リフト量ＶＬＩＦＴとの関係を示している。弁作動特性値ＣＡＭＰＦは、カムプロフィールそのものが与える弁の変位であり、弁リフト量ＶＬＩＦＴは、弁作動特定値ＣＡＭＰＦからバルブクリアランスＶＣＬＲを減じた実際の弁の変位である。弁作動特性値ＣＡＭＰＦ及び弁リフト量ＶＬＩＦＴは、ともに閉弁時を基準（＝０）としている。
【００２５】
本実施形態では、最大オーバーラップ時におけるオーバーラップ期間を所定のクランク角ＤＣＡ（ここでは、1°）毎に分割するとともに、ＥＣＵ４１に対し、分割した各区間の吸気弁２０及び排気弁２３の弁作動特性値ＣＡＭＰＦＩｎ，ＣＡＭＰＦＥｎ（ｎ＝１〜Ｎ）を記憶させておく。これらの弁作動特性値ＣＡＭＰＦＩｎ，ＣＡＭＰＦＥｎは、カムに固有のものであり、カムプロフィールを変更した場合は、その都度適合させる。エンジン１の運転時には、記憶されている各弁作動特性値ＣＡＭＰＦＩｎ，ＣＡＭＰＦＥｎから吸気弁２０又は排気弁２３に関するバルブクリアランスＶＣＬＲＩｎ，ＶＣＬＲＥｎを減算し、吸気弁２０が形成する区間開口面積ＶＡＲＥＡＩ、及び排気弁２３が形成する区間開口面積ＶＡＲＥＡＥを算出する。算出した区間開口面積ＶＡＲＥＡＩ，ＶＡＲＥＡＥは、区間毎に対応させた配列として記憶させる（図９）。前者の吸気側区間開口面積ＶＡＲＥＡＩの配列が「第１の配列」に、後者の排気側区間開口面積ＶＡＲＥＡＥの配列が「第２の配列」に相当する。有効開口面積ＡＳＵＭＯＬの演算では、記憶されている配列を参照して、吸気側及び排気側区間開口面積ＶＡＲＥＡＩ，ＶＡＲＥＡＥのうち小さい方をその区間の実質的な区間開口面積ＶＡＲＥＡｎとして選択し、選択した区間開口面積ＶＡＲＥＡｎをオーバーラップ期間に渡り積算して、有効開口面積ＡＳＵＭＯＬを算出する。
【００２６】
図５は、開口面積配列作成ルーチンのフローチャートである。このルーチンは、所定の時間毎に実行される。
Ｓ２０１では、列番号表示値ｎに１を加算する。列番号表示値ｎは、このルーチンによる配列の作成が終了するたびに、０に設定される。
Ｓ２０２では、列番号表示値ｎにより特定される吸気弁２０及び排気弁２３の弁作動特性値ＣＡＭＰＦＩｎ，ＣＡＭＰＦＥｎを読み込む。
【００２７】
Ｓ２０３では、読み込んだ弁作動特性値ＣＡＭＰＦＩｎ，ＣＡＭＰＦＥｎからそれぞれのバルブクリアランスＶＣＬＲＩｎ，ＶＣＬＲＥｎを減算し、吸気弁２０及び排気弁２３の弁リフト量ＶＬＩＦＴＩｎ、ＶＬＩＦＴＥｎを算出する。なお、バルブクリアランスＶＣＬＲＩｎ，ＶＣＬＲＥｎは、冷却水温度Ｔｗ及び排気温度Ｔｅｘ等をもとに、推定することができる。
【００２８】
ＶＬＩＦＴＩｎ＝ＣＡＭＰＦＩｎ−ＶＣＬＲＩｎ  ・・・（１０ａ）
ＶＬＩＦＴＥｎ＝ＣＡＭＰＦＥｎ−ＶＣＬＲＥｎ  ・・・（１０ｂ）
Ｓ２０４では、算出した弁リフト量ＶＬＩＦＴＩｎ，ＶＬＩＦＴＥｎに対し、流量感度係数Ｃｖに応じた係数ＫＣＶＩ＃，ＫＣＶＥ＃と、吸気弁２０又は排気弁２３の弁体投影面積ＶＡＲＥＡＩ０＃，ＶＡＲＥＡＥ０＃とを乗算し、吸気側及び排気側区間開口面積ＶＡＲＥＡＩ，ＶＡＲＥＡＥを算出する。なお、流量感度係数Ｃｖは、弁リフト量ＶＬＩＦＴに対し、理論的に与えられる流量と、実際に与えられる流量との比で表し、弁開期間がオーバーラップする低リフト領域では、弁リフト量ＶＬＩＦＴにほぼ比例する。係数ＫＣＶＩ＃，ＫＣＶＥ＃は、流量感度係数Ｃｖが描く近似直線の傾きとして算出し（図８）、固定値としてＥＣＵ４１に記憶させる。また、一気筒当たりに設けられる吸気弁２０又は排気弁２３の数をａ，ｂとし、吸気弁２０及び排気弁２３のシート当接部径をＶＤＩ，ＶＤＥとする。
【００２９】
ＶＡＲＥＡＩ＝ＶＬＩＦＴＩｎ×ＫＣＶＩ＃×ＶＡＲＥＡＩ０＃×ａ  ・・・（１１ａ）
ＶＡＲＥＡＥ＝ＶＬＩＦＴＥｎ×ＫＣＶＥ＃×ＶＡＲＥＡＥ０＃×ｂ  ・・・（１１ａ）
ＶＡＲＥＡＩ０＃＝（ＶＤＩ／２）＾２×π  ・・・（１２ａ）
ＶＡＲＥＡＥ０＃＝（ＶＤＥ／２）＾２×π  ・・・（１２ｂ）
Ｓ２０５では、算出した区間開口面積ＶＡＲＥＡＩ，ＶＡＲＥＡＥを列番号表示値ｎと対応させて記憶する。
【００３０】
ＶＡＲＥＡＩｎ＝ＶＡＲＥＡＩ  ・・・（１３ａ）
ＶＡＲＥＡＥｎ＝ＶＡＲＥＡＥ  ・・・（１３ｂ）
Ｓ２０６では、列番号表示値ｎが最終列番号Ｎに達したか否かを判定する。達したときは、Ｓ２０７へ進み、達していないときは、Ｓ２０１へ戻り、次の列について吸気側及び排気側区間開口面積ＶＡＲＥＡＩ，ＶＡＲＥＡＥを算出し、記憶する。
【００３１】
Ｓ２０７では、列番号表示値ｎを０に設定する。
図６は、有効開口面積演算ルーチンのフローチャートである。
Ｓ３０１では、吸気カム捻り角ＡＮＧＩ及び排気カム捻り角ＡＮＧＥを読み込む。
Ｓ３０２では、読み込んだカム捻り角ＡＮＧＩ，ＡＮＧＥをもとに、排気側区間開口面積ＶＡＲＥＡＥの配列を相対的にずらし、現状のオーバーラップ期間に適合させる。すなわち、カム捻り角ＡＮＧＩ，ＡＮＧＥをもとに、最大オーバーラップ時からの吸気弁開時期に対する排気弁閉時期の相対変化量ＳＩＦＴＥＶＣを算出し、吸気側区間開口面積ＶＡＲＥＡＩの配列に対し、排気側区間開口面積ＶＡＲＥＡＥの配列を算出した相対変化量ＳＩＦＴＥＶＣに応じた列数だけ前進させる（図１０）。たとえば、最大オーバーラップ時からの吸気弁開時期及び排気弁閉時期の変化量が夫々クランク角で３０°及び１０°である場合は、相対変化量ＳＩＦＴＥＶＣが４０°であり、配列の作成に当たりオーバーラップ期間を１°毎に分割しているので、排気側区間開口面積ＶＡＲＥＡＥの配列を４０°に応じた列数（＝４０）だけ前進させる。
【００３２】
Ｓ３０３では、列番号表示値ｎに１を加算する。
Ｓ３０４では、吸気側区間開口面積ＶＡＲＥＡＩの配列及び前進させた排気側区間開口面積ＶＡＲＥＡＥの配列から、列番号表示値ｎに対応する列の区間開口面積ＶＡＲＥＡＩｎ，ＶＡＲＥＡＥｎを読み込む。
Ｓ３０５では、吸気側区間開口面積ＶＡＲＥＡＩｎが排気側区間開口面積ＶＡＲＥＡＥｎよりも大きいか否かを判定する。大きいときは、Ｓ３０６へ進み、それ以外のときは、Ｓ３０７へ進む。
【００３３】
Ｓ３０６では、排気側区間開口面積ＶＡＲＥＡＥｎをその区間についての実質的な区間開口面積ＶＡＲＥＡｎとする。
Ｓ３０７では、吸気側区間開口面積ＶＡＲＥＡＩｎをその区間についての実質的な区間開口面積ＶＡＲＥＡｎとする。
Ｓ３０８では、各区間について算出した区間開口面積ＶＡＲＥＡｎをオーバーラップ期間に渡り積算する。
【００３４】
ＳＩＧＭＡ＝ＳＩＧＭＡ＋ＶＡＲＥＡｎ×ＤＣＡ  ・・・（１３）
Ｓ３０９では、列番号表示値ｎが最終列番号Ｎに達したか否かを判定する。達したときは、Ｓ３１０へ進み、達していないときは、Ｓ３０３へ戻り、次の列について区間開口面積ＶＡＲＥＡｎを算出し、積算する。
Ｓ３１０では、算出した積算値ＳＩＧＭＡを有効開口面積ＡＳＵＭＯＬとする。
【００３５】
Ｓ３１１では、列番号表示値ｎ及び積算値ＳＩＧＭＡを０に設定する。
なお、ここでは、区間開口面積ＶＡＲＥＡｎをオーバーラップ期間に渡り積算して有効開口面積ＡＳＵＭＯＬを算出することとしているが、区間開口面積ＶＡＲＥＡｎの選択及び積算は、排気弁リフト量ＶＬＩＦＴＥが０となった時点で中断してもよい。
算出した有効開口面積ＡＳＵＭＯＬは、既述の内部ＥＧＲ量推定ルーチンにおいて、吹抜ガス量ＭＲＥＳＯＬの演算（Ｓ１１２）に用いられる。
【００３６】
本実施形態に関し、図２に示すフローチャートのＳ１０１（吸気圧力Ｐｉｎ及び排気圧力Ｐｅｘの検出），１０３〜１１２の処理が吹抜ガス量算出手段を、同フローチャートのＳ１０２の処理が残留ガス量算出手段を、同フローチャートのＳ１１３の処理が内部ＥＧＲ量算出手段を構成する。また、図５及び６に示すフローチャート全体が有効開口面積算出手段を構成する。
【００３７】
本実施形態によれば、次のような効果を得ることができる。
第１に、最大オーバーラップ時におけるオーバーラップ期間を所定のクランク角ＤＣＡ毎に分割するとともに、分割した各区間について実質的な区間開口面積ＶＡＲＥＡｎを算出し、算出した区間開口面積ＶＡＲＥＡｎを積算して有効開口面積ＡＳＵＭＯＬを算出することとしたので、バルブクリアランス等の変化により実際の弁リフト量ＶＬＩＦＴが変化したとしても、その影響を区間毎に反映させて、正確な有効開口面積ＡＳＵＭＯＬを算出することができる。
【００３８】
第２に、区間毎の開口面積ＶＡＲＥＡｎを積算して有効開口面積ＡＳＵＭＯＬを算出することとしたので、カムプロフィールに対する連続的な演算による場合に対し、有効開口面積ＡＳＵＭＯＬの算出精度を確保しつつ、演算負荷を抑えることができる。
第３に、吸気側及び排気側区間開口面積ＶＡＲＥＡＩ，ＶＡＲＥＡＥを配列化し、バルブタイミングが変更された場合は、これらの配列をバルブタイミングの変化に応じた量だけ相対的にずらすこととしたので、変更後のバルブタイミングのもとで区間開口面積ＶＡＲＥＡＩ，ＶＡＲＥＡＥを改めて演算する必要がなく、あらゆるバルブタイミングに対して簡単に対応することができる。
【００３９】
以上では、カムプロフィール自体は一定とし、バルブタイミングのみを変化させることでオーバーラップ期間を変化させる場合を例に説明したが、本発明は、これに限らず、カムプロフィールの変化を伴ってオーバーラップ期間を変化させるものに適用することもできる。
【図面の簡単な説明】
【００４０】
【図１】本発明の一実施形態に係るエンジンの構成
【図２】内部ＥＧＲ量推定ルーチンのフローチャート
【図３】排気比熱比設定テーブル
【図４】混合気比熱比設定テーブル
【図５】開口面積配列作成ルーチンのフローチャート
【図６】有効開口面積演算ルーチンのフローチャート
【図７】クランク角に対する弁作動特性値及び弁リフト量の変化
【図８】弁リフト量に対する流量感度係数の変化
【図９】最大オーバーラップ時に作成される配列の概念
【図１０】実際のバルブタイミングに対して設定される配列の概念
【符号の説明】
【００４１】
１…エンジン、１１…吸気通路、１２…エアクリーナ、１３…スロットル弁、１４…サージタンク、１６…吸気ポート、１７…インジェクタ、１８…燃焼室、１９…ピストン、２０…吸気弁、２１…吸気カム、２２…排気ポート、２３…排気弁、２４…排気カム、２５…吸気側可変動弁装置、２６…排気側可変動弁装置、２７…点火プラグ、２８…排気通路、４１…エンジンコントロールユニット、５１…エアフローメータ、５２…吸気圧力センサ、５３…冷却水温度センサ、５４…クランク角センサ、５５…排気圧力センサ、５６…排気温度センサ、５７…酸素センサ、５８…アクセルセンサ、５９，６０…カム角センサ。
【図１】
2005307803.tif 000003
【図２】
2005307803.tif 000004
【図３】
2005307803.tif 000005
【図４】
2005307803.tif 000006
【図５】
2005307803.tif 000007
【図６】
2005307803.tif 000008
【図７】
2005307803.tif 000009
【図８】
2005307803.tif 000010
【図９】
2005307803.tif 000011
【図１０】
2005307803.tif 000012