JP3358449B2 - Exhaust system pressure estimation device for internal combustion engine - Google Patents

Exhaust system pressure estimation device for internal combustion engine

Info

Publication number
JP3358449B2
JP3358449B2 JP17891996A JP17891996A JP3358449B2 JP 3358449 B2 JP3358449 B2 JP 3358449B2 JP 17891996 A JP17891996 A JP 17891996A JP 17891996 A JP17891996 A JP 17891996A JP 3358449 B2 JP3358449 B2 JP 3358449B2
Authority
JP
Japan
Prior art keywords
exhaust
temperature
engine
intake air
detecting
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Lifetime
Application number
JP17891996A
Other languages
Japanese (ja)
Other versions
JPH1026049A (en
Inventor
浩之 糸山
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nissan Motor Co Ltd
Original Assignee
Nissan Motor Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Priority to JP17891996A priority Critical patent/JP3358449B2/en
Application filed by Nissan Motor Co Ltd filed Critical Nissan Motor Co Ltd
Priority to US08/867,980 priority patent/US6012431A/en
Priority to GB9825391A priority patent/GB2329040B/en
Priority to DE19758641A priority patent/DE19758641B4/en
Priority to GB9711483A priority patent/GB2313927B/en
Priority to DE19723210A priority patent/DE19723210B4/en
Publication of JPH1026049A publication Critical patent/JPH1026049A/en
Priority to US09/427,004 priority patent/US6167342A/en
Priority to US09/699,342 priority patent/US6298299B1/en
Application granted granted Critical
Publication of JP3358449B2 publication Critical patent/JP3358449B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Classifications

    • Y02T10/121
    • Y02T10/47

Description

【発明の詳細な説明】DETAILED DESCRIPTION OF THE INVENTION

【0001】[0001]

【発明の属する技術分野】本発明は、内燃機関のEGR
制御等に必要な排気系の圧力を推定する技術に関する。
The present invention relates to an EGR for an internal combustion engine.
The present invention relates to a technique for estimating an exhaust system pressure required for control and the like.

【0002】[0002]

【従来の技術】従来の内燃機関の吸気系や排気系の圧力
検出方式としては、これら圧力を直接検出するセンサを
設けることが一般的である。また、特願平7−9263
1号に開示されるように、エアフローメータの出力等か
ら推定する技術もある。
2. Description of the Related Art As a conventional pressure detection method for an intake system and an exhaust system of an internal combustion engine, a sensor for directly detecting these pressures is generally provided. In addition, Japanese Patent Application No. 7-9263
As disclosed in No. 1, there is also a technique of estimating from an output of an air flow meter or the like.

【0003】[0003]

【発明が解決しようとする課題】しかしながら、上記の
ようにセンサを設ける場合は、コストが高くつくことは
勿論であるが、排気圧検出用のセンサの場合、EGR等
排気を循環させたり、高温雰囲気にさらしたりするとい
うことにより、耐久性の面でも非常に厳しい。また、耐
久性を確保できるようにすると、センサの時定数が非常
に大きくなり、過渡時には正確な値を表示できなくな
り、制御への反映は困難であることが判明した。
However, in the case where the sensor is provided as described above, the cost is of course high. However, in the case of the sensor for detecting the exhaust pressure, the exhaust gas such as EGR is circulated, Due to exposure to the atmosphere, durability is extremely severe. Further, it has been found that when the durability can be ensured, the time constant of the sensor becomes extremely large, and it is difficult to display an accurate value in a transient state, and it is difficult to reflect the value on the control.

【0004】また、排気圧に影響を与える排気温度は、
噴射時期等の着火時期を代表するパラメータや、シリン
ダ内の空気流動によって変化するため、これに対して前
記特願平7−92631号の技術では、対応しきれない
という問題があった。
[0004] The exhaust temperature which affects the exhaust pressure is:
Since it changes depending on parameters representing the ignition timing such as the injection timing and the air flow in the cylinder, the technique of Japanese Patent Application No. 7-92631 cannot cope with this.

【0005】[0005]

【課題を解決するための手段】このため、請求項1に係
る発明は図1に示すように、機関に噴射される燃料の噴
射量を検出する燃料噴射量検出手段と、機関に吸入され
る空気の温度を検出する吸入空気温度検出手段と、前記
燃料の噴射時期を検出する燃料噴射時期検出手段と、前
記吸入空気の流量を検出する吸入空気流量検出手段と、
機関の回転速度を検出する機関回転速度検出手段と、
料噴射量に基づいて基本排気温度を演算し、該基本排気
温度を吸入空気温度と前回推定された排気圧力と燃料噴
射時期とで補正してシリンダからの排気温度を演算し、
吸入空気流量と機関回転速度とから求めたシリンダから
の排気流量と前記シリンダからの排気温度とに基づいて
排気圧を推定する排気系圧力推定手段と、を含んで構成
したことを特徴とする。
Therefore, according to the present invention, as shown in FIG. 1, a fuel injection amount detecting means for detecting an injection amount of fuel injected into an engine, and a fuel injected into the engine. Intake air temperature detecting means for detecting the temperature of the air, fuel injection timing detecting means for detecting the fuel injection timing, intake air flow rate detecting means for detecting the flow rate of the intake air,
And the engine rotational speed detecting means for detecting a rotational speed of the engine, fuel
The basic exhaust temperature is calculated based on the fuel injection amount, and the basic exhaust temperature is calculated.
The intake air temperature, the previously estimated exhaust pressure, and the fuel injection
Calculate the exhaust temperature from the cylinder by correcting with the firing timing,
From the cylinder obtained from the intake air flow rate and the engine speed
Based on the exhaust gas flow rate and the exhaust gas temperature from the cylinder.
Exhaust pressure estimating means for estimating exhaust pressure.

【0006】(作用・効果)燃料噴射量検出手段によって
検出された燃料噴射量によって燃焼時の発熱量が決まる
ので、基本的な排気温度が定まる。また、吸入空気温度
が変わると燃焼温度が変化するため、排気温度も変化す
る。
(Operation / Effect) Since the amount of heat generated during combustion is determined by the fuel injection amount detected by the fuel injection amount detecting means, the basic exhaust temperature is determined. Further, when the intake air temperature changes, the combustion temperature changes, so that the exhaust gas temperature also changes.

【0007】更に、燃料の噴射時期についても噴射時期
を早める (進角量を大きくする) と、燃焼状態が改善さ
れ、シリンダ排出時の温度が下げられるので、排気温度
が低下し、逆に噴射時期を遅くすると前記と逆の作用で
排気温度が上昇するというように排気温度に影響を与え
る。そこで、燃料噴射量と吸入空気温度と燃料噴射時期
とに基づいて排気温度が精度良く推定される。
Further, when the injection timing of the fuel is advanced (the advance amount is increased), the combustion state is improved, and the temperature at the time of discharging the cylinder is lowered. If the timing is delayed, the exhaust gas temperature is affected such that the exhaust gas temperature increases due to the opposite effect. Therefore, the exhaust gas temperature is accurately estimated based on the fuel injection amount, the intake air temperature, and the fuel injection timing.

【0008】一方、吸入空気流量と機関回転速度とに基
づいてシリンダ排出時の排気流量を求めることができ
る。なお、定常状態であれば、 (単位時間当りの) 吸入
空気流量がシリンダから排出される排気流量と略一致す
るが、過渡状態では吸気系の容積等で異なってくる。そ
こで、吸入空気流量と機関回転速度とに基づいてシリン
ダから排出される1回当りの排気量を求めてから、再度
機関回転速度を用いてシリンダから排出される (単位時
間当りの) 排気流量に換算することにより、精度良く排
気流量を演算することができる。
On the other hand, the exhaust flow rate at the time of cylinder discharge can be obtained based on the intake air flow rate and the engine speed. In the steady state, the flow rate of the intake air (per unit time) substantially matches the flow rate of the exhaust gas discharged from the cylinder, but in the transient state, it differs depending on the volume of the intake system. Therefore, the amount of exhaust per cylinder discharged from the cylinder is calculated based on the intake air flow rate and the engine speed, and then the exhaust flow rate (per unit time) discharged from the cylinder using the engine speed is calculated again. By performing the conversion, the exhaust flow rate can be accurately calculated.

【0009】そして、前記推定した排気温度と、演算し
た排気流量と、に基づいてベルヌーイの式から排気圧力
を求めることができる。また、排気温度の推定に際し、
燃料噴射量に基づいて求めた基本排気温度を、前記吸入
空気流量と燃料噴射時期の他、前回推定された排気圧力
でも補正して排気温度を推定する。すなわち、断熱変化
の関係から排気圧力が上昇すると排気温度も所定の割合
で上昇するので、前回の排気圧力の推定値による補正も
同時に行い、より精度良く排気温度を推定する。
Then, the exhaust pressure can be obtained from Bernoulli's equation based on the estimated exhaust temperature and the calculated exhaust flow rate. Also, when estimating the exhaust gas temperature,
The basic exhaust temperature obtained based on the fuel injection amount is
In addition to the air flow and fuel injection timing, the previously estimated exhaust pressure
However, the exhaust gas temperature is estimated after correction. That is, adiabatic change
As the exhaust pressure rises, the exhaust temperature also increases
, So the correction based on the previous estimated exhaust pressure
Simultaneously, the exhaust gas temperature is more accurately estimated.

【0010】[0010]

【0011】また、請求項2に係る発明は、図2に示す
ように、機関に噴射される燃料の噴射量を検出する燃料
噴射量検出手段と、機関に吸入される空気の温度を検出
する吸入空気温度検出手段と、機関の吸気系に備えられ
たスワール制御弁によるスワール制御状態を検出するス
ワール制御状態検出手段と、前記吸入空気の流量を検出
する吸入空気流量検出手段と、機関の回転速度を検出す
る機関回転速度検出手段と、燃料噴射量に基づいて基本
排気温度を演算し、該基本排気温度を吸入空気温度と前
回推定された排気圧力とスワール制御状態とで補正して
シリンダからの排気温度を演算し、吸入空気流量と機関
回転速度とから求めたシリンダからの排気流量と前記シ
リンダからの排気温度とに基づいて排気圧を推定する
気系圧力推定手段と、を含んで構成したことを特徴とす
る。
Further, as shown in FIG. 2, the invention according to claim 2 detects fuel injection amount detecting means for detecting the amount of fuel injected into the engine, and detects the temperature of air taken into the engine. Intake air temperature detecting means, swirl control state detecting means for detecting a swirl control state by a swirl control valve provided in an intake system of the engine, intake air flow rate detecting means for detecting a flow rate of the intake air, and rotation of the engine. Engine speed detection means for detecting the speed and basic based on the fuel injection amount
Calculate the exhaust gas temperature, and calculate the basic exhaust gas temperature as the intake air temperature.
Correction with the estimated exhaust pressure and swirl control state
Calculates the exhaust temperature from the cylinder, and calculates the intake air flow rate and engine
Exhaust flow rate from the cylinder obtained from the rotational speed
Exhaust pressure estimating means for estimating the exhaust pressure based on the exhaust temperature from the cylinder .

【0012】(作用・効果)吸気系にスワール制御弁を備
えたものでは、スワール制御弁によってスワール強さが
制御され、それによって燃焼性が変化するため、排気温
度も変わってくる。そこで、該スワール制御弁によるス
ワールの制御状態を検出し、前記燃料噴射量と、吸入空
気温度に加えて該スワール制御状態も考慮して排気温度
を推定する。なお、噴射時期も含めて排気温度を推定し
てもよいことは、勿論であり、より推定精度が向上す
る。
(Operation / Effect) In a system having a swirl control valve in the intake system, the swirl control valve controls the swirl intensity and thereby changes the flammability, so that the exhaust gas temperature also changes. Therefore, the control state of the swirl by the swirl control valve is detected, and the exhaust temperature is estimated in consideration of the swirl control state in addition to the fuel injection amount and the intake air temperature. In addition, it goes without saying that the exhaust gas temperature may be estimated including the injection timing, and the estimation accuracy is further improved.

【0013】このようにして精度良く推定された排気温
度と、吸入空気流量と機関回転速度とから求めた排気流
量とに基づいて排気圧を精度良く推定することができる
ことがわかる。
[0013] it is possible to accurately estimate the exhaust pressure on the basis of the exhaust temperature estimated accurately in this manner, the exhaust flow rate determined from the intake air flow rate and the engine rotational speed
You can see that.

【0014】そこで、燃料噴射量によって燃料燃焼時に
発熱量に応じた基本排気温度を求め、これを、吸入空気
温度と前回推定された排気圧力とスワール制御状態とで
補正することにより、精度良く排気温度を推定演算する
ことができ、更にシリンダからの排気流量と排気温度と
に基づいて排気温度を精度良く推定することができる。
Therefore, the basic exhaust temperature corresponding to the calorific value at the time of fuel combustion is obtained from the fuel injection amount, and the basic exhaust temperature is corrected by the intake air temperature, the exhaust pressure estimated last time, and the swirl control state, so that the exhaust gas is accurately discharged. The temperature can be estimated and calculated, and the exhaust temperature can be accurately estimated based on the exhaust flow rate from the cylinder and the exhaust temperature.

【0015】また、請求項3に係る発明は、前記基本排
気温度のスワール制御状態による補正は、スワール制御
弁の開度と機関回転速度とに基づいて設定される補正係
数により行うことを特徴とする。(作用・効果)スワール
制御弁の開度が大きいときは、機関回転速度が変化して
もスワールの影響は全域にわたって小さいが、スワール
制御弁の開度が小さいときは、低速域では適度なスワー
ルの発生により燃焼性が改善されて排気温度が下がるの
に対し高速域ではスワールが強すぎて燃焼性が低下し、
また、絞り作用によりシリンダ吸入空気量が減少する等
によって、排気温度は増大する。
Further, the invention according to claim 3 is characterized in that the correction of the basic exhaust gas temperature in the swirl control state is performed by a correction coefficient set based on the opening degree of the swirl control valve and the engine speed. I do. (Operation / Effect) When the swirl control valve opening is large, the effect of the swirl is small even if the engine speed changes, but when the swirl control valve opening is small, an appropriate swirl is effective at low speeds. In the high-speed range, the swirl is too strong and the flammability decreases.
Further, the exhaust gas temperature increases due to a decrease in the cylinder intake air amount due to the throttling action.

【0016】そこで、スワール弁開度と機関回転速度と
に基づいて排気温度の補正係数を設定し、該補正係数を
用いて排気温度を補正することにより、排気温度の推定
精度をより高めることができる。また、請求項4に係る
発明は、前記基本排気温度のスワール制御状態による補
正は、吸入空気流量と機関回転速度とスワール制御弁の
開度とから演算したスワール流速相当値に基づいて設定
される補正係数により行うことを特徴とする。
Therefore, a correction coefficient for the exhaust gas temperature is set based on the swirl valve opening and the engine speed, and the exhaust gas temperature is corrected using the correction coefficient, thereby further improving the estimation accuracy of the exhaust gas temperature. it can. Further, in the invention according to claim 4 , the correction of the basic exhaust gas temperature by the swirl control state is set based on a swirl flow rate equivalent value calculated from the intake air flow rate, the engine rotation speed, and the opening of the swirl control valve. The correction is performed by using a correction coefficient.

【0017】(作用・効果)吸入空気流量と機関回転速度
とでシリンダに吸入されるときの吸入空気流量を演算で
き、このシリンダ吸入空気流量をスワール制御弁の開度
で除算することにより、スワール流速に相当する値が求
められる。そして、スワール流速がある値のときに燃焼
性が最も良く改善されて排気温度が低下し、スワール流
速がそれより大きいときも、小さいときも燃焼性が低下
して排気温度は上昇する傾向がある。
(Operation / Effect) The intake air flow rate when the intake air is taken into the cylinder can be calculated based on the intake air flow rate and the engine rotation speed. By dividing this cylinder intake air flow rate by the opening of the swirl control valve, the swirl can be calculated. A value corresponding to the flow velocity is determined. When the swirl flow velocity is at a certain value, the flammability is improved most and the exhaust temperature is reduced. When the swirl flow velocity is higher or lower, the flammability lowers and the exhaust temperature tends to increase. .

【0018】そこで、前記吸入空気流量と機関回転速度
とスワール制御弁の開度とから演算したスワール流速相
当値に基づいて排気温度の補正係数を設定し、該補正係
数を用いて排気温度を補正することにより、排気温度の
推定精度をより高めることができる。
Therefore, a correction coefficient for the exhaust gas temperature is set based on the swirl flow velocity equivalent value calculated from the intake air flow rate, the engine rotation speed, and the opening of the swirl control valve, and the exhaust gas temperature is corrected using the correction coefficient. By doing so, the estimation accuracy of the exhaust gas temperature can be further improved .

【0019】[0019]

【0020】[0020]

【0021】[0021]

【0022】[0022]

【発明の実施の形態】以下に本発明の実施形態を図に基
づいて説明する。一実施形態の全体構成を示す図4にお
いて、過給機1は、エアフィルタ2でダストを除去され
て吸気通路3に吸入された空気を吸気コンプレッサ1A
により圧縮過給して下流側の吸気マニホールド4へ送り
込む。
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS An embodiment of the present invention will be described below with reference to the drawings. In FIG. 4 showing the overall configuration of one embodiment, a supercharger 1 is configured to remove air that has been dust-removed by an air filter 2 and that is taken into an intake passage 3 by an intake compressor 1A.
, And is fed into the downstream intake manifold 4.

【0023】一方、機関 (ディーゼル機関) 5の燃焼室
に装着された燃料噴射ノズル6には、噴射ポンプ7から
各気筒に分配して燃料が圧送供給され、該燃料噴射ノズ
ル6から燃焼室に向けて燃料が噴射され、該噴射された
燃料は圧縮行程末期に着火して燃焼される。また、排気
マニホールド8と吸気マニホールド4とを結んでEGR
制御弁9を介装したEGR通路10が接続されると共
に、前記吸気通路3の吸気コンプレッサ1Aの上流側に
EGR制御時に吸気を絞って排気圧と吸気圧との差圧を
拡大してEGRしやすくするためのスロットル弁31が
介装され、主としてアイドル時や低負荷時に排気改善,
騒音対策のために前記スロットル弁31を絞ると同時に
EGR制御弁9の開度を制御してEGR制御を行う。具
体的には、バキュームポンプ11からの負圧を電磁弁3
2を介してダイアフラム装置33に導いて前記絞り弁3
1を絞ると同時に、前記負圧をデューティ制御される電
磁弁12で大気との希釈割合を制御することによって前
記EGR制御弁9の圧力室に導かれる圧力を制御し、も
って開度を制御することによりEGR率を制御してい
る。これらEGR率や燃料噴射制御は、コントロールユ
ニット13により行われる。
On the other hand, fuel is distributed to each cylinder from an injection pump 7 and supplied to the fuel injection nozzle 6 mounted in the combustion chamber of the engine (diesel engine) 5 under pressure, and the fuel is supplied from the fuel injection nozzle 6 to the combustion chamber. Fuel is injected toward the fuel, and the injected fuel is ignited and burned at the end of the compression stroke. The EGR is connected to the exhaust manifold 8 and the intake manifold 4.
An EGR passage 10 provided with a control valve 9 is connected. At the time of EGR control, the intake air is throttled to the upstream side of the intake compressor 1A in the intake passage 3 to expand the differential pressure between the exhaust pressure and the intake pressure to perform EGR. A throttle valve 31 is provided for ease of exhaustion, and exhaust is improved mainly at idle or low load.
The EGR control is performed by controlling the opening of the EGR control valve 9 at the same time as the throttle valve 31 is throttled to reduce noise. Specifically, the negative pressure from the vacuum pump 11 is
2 to the diaphragm device 33 through the throttle valve 3
At the same time as reducing the value of 1, the solenoid valve 12 that duty-controls the negative pressure controls the dilution ratio with the atmosphere, thereby controlling the pressure guided to the pressure chamber of the EGR control valve 9, thereby controlling the opening. Thus, the EGR rate is controlled. The EGR rate and the fuel injection control are performed by the control unit 13.

【0024】前記EGR制御弁9には、弁体のリフト量
を検出するリフトセンサ34が設置されている。燃焼後
の排気は、排気マニホールド8より前記過給機1の排気
タービン1Bを回転駆動させた後、排気中に含まれるパ
ーティキュレート (排気微粒子) 等がフィルタ14で捕
集され、マフラー15で消音された後に大気中に放出さ
れる。
The EGR control valve 9 is provided with a lift sensor 34 for detecting a lift amount of the valve body. After the combustion, the exhaust turbine 1B of the supercharger 1 is driven to rotate by the exhaust manifold 8, and the particulates (exhaust particulates) and the like contained in the exhaust are collected by the filter 14 and silenced by the muffler 15. After being released into the atmosphere.

【0025】前記過給機1の吸気コンプレッサ1A上流
の吸気通路3には、吸入空気流量を検出するエアフロー
メータ16が設けられ、また、機関回転速度Neを検出
する回転速度センサ17、前記燃料噴射ポンプ7のコン
トロールレバー開度 (アクセル開度) を検出するレバー
開度センサ18、水温を検出する水温センサ19等が設
けられ、これらの検出値に基づいて吸気系圧力,排気系
圧力を検出しつつシリンダ吸入空気量に見合った燃料の
許容最大噴射量が後述するようにして設定される。
In the intake passage 3 upstream of the intake compressor 1A of the supercharger 1, an air flow meter 16 for detecting an intake air flow rate is provided, a rotation speed sensor 17 for detecting an engine rotation speed Ne, and the fuel injection. A lever opening sensor 18 for detecting the control lever opening (accelerator opening) of the pump 7, a water temperature sensor 19 for detecting the water temperature, and the like are provided, and the intake system pressure and the exhaust system pressure are detected based on these detected values. Meanwhile, the allowable maximum injection amount of fuel corresponding to the cylinder intake air amount is set as described later.

【0026】以下、コントロールユニット13による各
種演算について説明する。図5は、吸入空気流量Qas
0の演算のフローに示す。ステップ (図ではSと記す。
以下同様) 1では、前記エアフローメータ16の出力電
圧Usを読み込み、ステップ2で図6に示すような特性
テーブル (電圧−吸入空気流量変換テーブル) を用いて
吸入空気流量Qas0 dに変換し、ステップ3で平均
化処理を行い、Qas0として処理を終了する。
Hereinafter, various calculations by the control unit 13 will be described. FIG. 5 shows the intake air flow rate Qas
This is shown in the flow of the operation of 0. Step (denoted by S in the figure)
In step 1, the output voltage Us of the air flow meter 16 is read, and in step 2, the intake air flow rate Qas0 is calculated using a characteristic table (voltage-intake air flow rate conversion table) as shown in FIG. Then, the averaging process is performed in step 3, and the process ends as Qas0.

【0027】図7は、シリンダへ吸入される1吸気行程
当りの吸入空気量の演算フローであり、機関回転に同期
したタイミングで演算する。ステップ11では、機関回転
速度Neを読み込み、ステップ12では、図示する式によ
り前記図5で処理された吸入空気流量Qas0から1吸
気行程当りの吸入空気量Qac0を演算する。
FIG. 7 is a calculation flow of the amount of intake air per one intake stroke drawn into the cylinder, and is calculated at a timing synchronized with the engine rotation. In step 11, the engine rotational speed Ne is read, and in step 12, the intake air amount Qac0 per intake stroke is calculated from the intake air flow rate Qas0 processed in FIG.

【0028】ステップ13では、所定回数前に演算された
吸入空気量Qac0をQacnとして処理を終了する。
この所定回前のQac0をQacnとするのは、エアフ
ローメータ16の測定位置からコレクタ入口までの輸送
遅れが存在するためである。ステップ14では、図示する
ような一次遅れ処理を行って、シリンダ吸入空気量Qa
cを演算して処理を終了する。
In step 13, the process ends with the intake air amount Qac0 calculated a predetermined number of times before as Qacn.
The reason why Qac0 before this predetermined time is Qacn is that there is a transport delay from the measurement position of the air flow meter 16 to the collector inlet. In step 14, a first-order lag process is performed as shown in FIG.
c is calculated and the process is terminated.

【0029】図8は、シリンダへ吸入、供給されてから
排出されるまでサイクル遅れ (デッドタイム処理) を行
うフローであり、シリンダ吸入空気量Qac,実噴射時
期ltist,燃料噴射量Qsol,吸気温度Tnに対
して各々所定サイクル分のデッドタイム処理を行う。こ
れは、求めようとする排気圧の元となるシリンダ吸入前
のデータを得るため行われる。
FIG. 8 is a flow chart for performing a cycle delay (dead time processing) from the intake and supply to the cylinder to the discharge from the cylinder. The cylinder intake air amount Qac, the actual injection timing ltist, the fuel injection amount Qsol, and the intake air temperature are shown. A dead time process for a predetermined cycle is performed for each Tn. This is performed in order to obtain data before cylinder intake, which is the basis of the exhaust pressure to be obtained.

【0030】図9は、以上のようにして求められた各値
を用いて、排気圧を演算するフローである。ステップ31
では、前記図8で処理された燃料噴射量サイクル処理値
Qsoldから、図10に示すテーブルから基本排気温度
Tehiを検索する。なお、燃料噴射量の増大に伴い燃
焼熱量が増大するから基本排気温度Texhiも比例的
に増大する特性としている。
FIG. 9 is a flow chart for calculating the exhaust pressure using the values obtained as described above. Step 31
Then, based on the fuel injection amount cycle processing value Qsold processed in FIG. 8, the basic exhaust temperature Tehi is searched from the table shown in FIG. Since the amount of combustion heat increases with an increase in the fuel injection amount, the basic exhaust temperature Texhi has a characteristic of proportionally increasing.

【0031】ステップ32では、吸気温度サイクル処理値
Tneにより、図11に示すテーブルから吸気温度による
排温補正係数Ktmpnを検索する。なお、吸気温度が
高いときは、排気温度も高くなるから、吸気温度の増大
に伴い排温補正係数Ktmpnも増大する特性としてい
る。ステップ33では、前回演算の排気圧Pexhn−1
により、図12に示すテーブルから排気圧力による排温補
正係数Kkmppを演算する。なお、排気圧が高いとき
ほど排気温度も増大するので、排温補正係数Ktmpp
もそのような特性としてある。
In step 32, the exhaust temperature correction coefficient Ktmpn based on the intake air temperature is searched from the table shown in FIG. 11 based on the intake air temperature cycle processing value Tne. When the intake air temperature is high, the exhaust gas temperature is also high. Therefore, the exhaust gas temperature correction coefficient Ktmpn increases with the intake air temperature. In step 33, the exhaust pressure Pexhn-1 of the previous calculation is calculated.
Thus, the exhaust temperature correction coefficient Kkmpp based on the exhaust pressure is calculated from the table shown in FIG. Since the exhaust gas temperature increases as the exhaust pressure increases, the exhaust temperature correction coefficient Ktmppp
Is also such a property.

【0032】ステップ34では、実噴射時期サイクル処理
値ltistdにより、図13に示すテーブルから噴射時
期による排温補正係数Ktmpiを演算する。なお、噴
射時期が遅くなるほど、燃焼終了が遅れる等の理由で排
気温度は高くなるから、排温補正係数Ktmpiもその
特性に合わせて設定してある。ステップ35では、前記ス
テップ31で求めた基本排気温度Texhiを、ステップ
32,33,34で求めた各排温補正係数Ktmpn,Ktm
pp,Ktmpiを用いて、次式で補正することによ
り、シリンダ排出温度Texhcを演算する。
In step 34, the exhaust temperature correction coefficient Ktmpi based on the injection timing is calculated from the table shown in FIG. 13 based on the actual injection timing cycle processing value ltistd. It should be noted that since the exhaust temperature increases as the injection timing becomes later, for example, the end of combustion is delayed, the exhaust temperature correction coefficient Ktmpi is also set in accordance with the characteristic. In step 35, the basic exhaust gas temperature Texhi obtained in step 31 is
Each exhaust temperature correction coefficient Ktmpn, Ktm obtained in 32, 33, 34
By using pp and Ktmpi, the cylinder discharge temperature Texhc is calculated by correcting the following equation.

【0033】 Texhc=Texhi×Ktmpn×Ktmpp×Ktmpi ステップ36では、前記シリンダ排出温度Texhiに、
次式に示す一次遅れ処理を行って排気温度Texhとす
る。 Texh=Texhn-1 × (1−KO) +Texhc×KO KOは定数 ステップ37では、シリンダ当りの排出空気量Qexh,
機関回転速度Ne,排気温度Texhと定数KPEX
H,OPEXH (標準状態における大気圧に相当する
値) を用いて、次式により排気圧基本値Pehbを演算
する。 Pehb= (Qexh×Ne/KC) 2 ×Texh×KPEXH+OPEXH ステップ38では、前記排気圧基本値Pexhbに、次式
に示す一次遅れ処理を施して排気圧Pexhとする。
Texhc = Texhi × Ktmpn × Ktmpp × Ktmpi In step 36, the cylinder discharge temperature Texhi is
The exhaust temperature Texh is obtained by performing a first-order lag process represented by the following equation. Texh = Texh n-1 × (1-KO) + Texhc × KO KO is a constant In step 37, the amount of exhaust air per cylinder Qexh,
Engine speed Ne, exhaust temperature Texh and constant KPEX
Using H, OPEXH (a value corresponding to the atmospheric pressure in the standard state), a basic exhaust pressure value Pehb is calculated by the following equation. Pehb = (Qexh × Ne / KC) 2 × Texh × KEXH + OPEXH In step 38, the exhaust pressure basic value Pexhb is subjected to a first-order lag process represented by the following equation to obtain the exhaust pressure Pexh.

【0034】 Pexh=Pexhn-1 × (1−KP) +Pexhb×KP KPは定数 次に、第2の実施形態について説明する。このもので
は、図14に示すように、吸気マニホールド4の各気筒の
吸気ポート部分に、低速時に吸気スワールを強化して燃
焼性を高めるためのスワール制御弁21が介装されてお
り、該スワール制御弁21の開度がコントロールユニット
13により機関の運転状態に応じて制御されるようになっ
ている。そして、本実施形態では、該スワール制御弁21
の開度に応じて排気圧を補正して推定する。
Pexh = Pexh n-1 × (1-KP) + Pexhb × KP KP is a constant Next, a second embodiment will be described. In this embodiment, as shown in FIG. 14, a swirl control valve 21 for strengthening the intake swirl at low speed to enhance the combustion performance is interposed at the intake port portion of each cylinder of the intake manifold 4, and the swirl control valve 21 is provided. The opening of control valve 21 is the control unit
The control is performed according to the operating state of the engine by means of 13. In the present embodiment, the swirl control valve 21
The exhaust pressure is corrected and estimated in accordance with the opening degree.

【0035】図3〜図12については、第1の実施形態と
共通に用いる。図15は、図8と同様に行われるサイクル
遅れ (デッドタイム処理) を行うフローであり、シリン
ダ吸入空気量Qac,スワール制御弁開度Riscv,
燃料噴射量Qsold,吸気温度Tneに対して各々所
定サイクル分のデッドタイム処理を行う。
3 to 12 are used in common with the first embodiment. FIG. 15 is a flow for performing a cycle delay (dead time processing) performed in the same manner as in FIG. 8, and includes a cylinder intake air amount Qac, a swirl control valve opening Riscv,
Dead time processing for a predetermined cycle is performed for each of the fuel injection amount Qsold and the intake air temperature Tne.

【0036】図16は、排気圧を演算するフローである。
ステップ41〜ステップ43は、前記ステップ31〜ステップ
33までと同様であるので説明を省略する。ステップ44で
は、前記図15でのスワール制御弁開度サイクル処理値R
iscvdと機関回転速度Neとから、図17に示すテー
ブルからスワール制御弁開度による排温補正係数Ktm
pscを検索する。なお、スワール制御弁21が全開のと
き(実質的にスワール制御弁を備えない場合に相当)
は、排温補正係数Ktmpsc=1で固定 (実質的に補
正を行わない) とし、スワール制御弁21の開度が小さい
ときは、機関が低速のときは適切なスワール強さとなっ
て燃焼性が改善される結果、排気温度は低下するが、機
関の高速時はスワールが過大となって燃焼性が悪化し、
また、絞り作用によりシリンダ吸入空気量が減少する等
の理由によって、排気温度は上昇するので、その特性に
合わせて排温補正係数Ktmpscを設定してある。
FIG. 16 is a flowchart for calculating the exhaust pressure.
Steps 41 to 43 correspond to steps 31 to
The description is omitted because it is the same as 33. In step 44, the swirl control valve opening cycle processing value R shown in FIG.
Based on the iscvd and the engine speed Ne, the exhaust temperature correction coefficient Ktm based on the swirl control valve opening is obtained from the table shown in FIG.
Search for psc. When the swirl control valve 21 is fully opened (corresponding to a case where the swirl control valve is not substantially provided)
Is fixed at the exhaust temperature correction coefficient Ktmpsc = 1 (substantially no correction is performed). When the opening degree of the swirl control valve 21 is small, the swirl strength becomes appropriate when the engine is at low speed, and the combustibility is reduced. As a result, the exhaust temperature decreases, but at high engine speeds, the swirl becomes excessive and the combustibility deteriorates.
Further, the exhaust gas temperature rises due to, for example, a decrease in the cylinder intake air amount due to the throttling action. Therefore, the exhaust temperature correction coefficient Ktmpsc is set in accordance with the characteristic.

【0037】ステップ45では、前記基本排気温度Tex
hiを、ステップ42,43,44で求めた各排温補正係数K
tmpn,Ktmpp,Ktmpscを用いて、次式で
補正することにより、シリンダ排出温度Texhcを演
算する。 Texhc=Texhi×Ktmpn×Ktmpp×K
tmpsc 以下、ステップ46〜ステップ48については、前記図9の
ステップ36〜ステップ38と同様であり、ステップ46で、
前記シリンダ排出温度Texhiに一次遅れ処理を行っ
て排気温度Texhとし、ステップ47で、図示の式によ
り排気圧基本値Pehbを演算し、ステップ48では、前
記排気圧基本値Pexhbに一次遅れ処理を施して排気
圧Pexhとする。
In step 45, the basic exhaust gas temperature Tex
hi is the exhaust temperature correction coefficient K obtained in steps 42, 43 and 44.
Using the tmpn, Ktmpp, and Ktmpsc, the cylinder discharge temperature Texhc is calculated by correcting the following equation. Texhc = Texhi × Ktmpn × Ktmpp × K
tmpsc Hereinafter, steps 46 to 48 are the same as steps 36 to 38 in FIG.
A first-order lag process is performed on the cylinder discharge temperature Texhi to obtain an exhaust gas temperature Texh. In step 47, a basic exhaust pressure value Pexhb is calculated by an equation shown in the drawing, and in step 48, a first-order lag process is performed on the basic exhaust pressure value Pexhb. To the exhaust pressure Pexh.

【0038】前記実施形態では、スワール制御弁の開度
と機関回転速度とから、3次元テーブルからの検索によ
って排温補正係数Ktmpscを設定したが、シリンダ
吸入空気量Qac,機関回転速度Ne,スワール弁開度
Rascvからスワール弁開度に対する吸入空気流量の
割合、つまり、スワールの流速に相当する値Vscを算
出し、該流速相当値Vscに応じて排温補正係数Ktm
pscを設定するようにしてもよい。
In the above-described embodiment, the exhaust temperature correction coefficient Ktmpsc is set by searching the three-dimensional table from the opening degree of the swirl control valve and the engine speed, but the cylinder intake air amount Qac, the engine speed Ne, the swirl The ratio of the intake air flow rate to the swirl valve opening degree, that is, the value Vsc corresponding to the swirl flow velocity is calculated from the valve opening degree Rascv, and the exhaust temperature correction coefficient Ktm is calculated according to the flow velocity equivalent value Vsc.
The psc may be set.

【0039】以下に、上記の実施形態を示す。図3〜図
12については、第1の実施形態と共通に用いる。図18
は、前記各実施形態と同様に行われるサイクル遅れ (デ
ッドタイム処理)を行うフローであり、シリンダ吸入空
気量Qac,スワール流速Vsc,燃料噴射量Qsol
d,吸気温度Tneに対して各々所定サイクル分のデッ
ドタイム処理を行う。
The above embodiment will be described below. FIG. 3 to FIG.
12 is used in common with the first embodiment. Fig. 18
Is a flow for performing a cycle delay (dead time processing) performed in the same manner as in the above-described embodiments, and includes a cylinder intake air amount Qac, a swirl flow velocity Vsc, and a fuel injection amount Qsol.
d, dead time processing for a predetermined cycle is performed for each intake air temperature Tne.

【0040】ここで、スワール流速Vscは、図示の式
のようにシリンダ部分の吸入空気流量Qac×Neをス
ワール制御弁の開度Rascvで除算して求められる。
図19は、排気圧を演算するフローである。ステップ51〜
ステップ53は、前記ステップ31〜ステップ33までと同様
であるので説明を省略する。
Here, the swirl flow rate Vsc is obtained by dividing the intake air flow rate Qac × Ne in the cylinder portion by the swirl control valve opening Rascv, as shown in the equation in the figure.
FIG. 19 is a flowchart for calculating the exhaust pressure. Step 51 ~
Step 53 is the same as steps 31 to 33, and a description thereof will not be repeated.

【0041】ステップ54では、前記図18でのスワール流
速相当値のサイクル処理値Vscから、図20に示すテー
ブルからスワール制御弁開度による排温補正係数Ktm
pscを検索する。これは、スワール流速相当値Vsc
がある値のときに燃焼状態が最もよくなって排気温度が
下がり、そのVscからどちらにずれても排気温度は上
昇するという特性に合わせて排温補正係数Ktmpsc
を設定してある。
In step 54, the exhaust gas temperature correction coefficient Ktm based on the swirl control valve opening is obtained from the cycle processing value Vsc corresponding to the swirl flow velocity in FIG.
Search for psc. This is the swirl flow velocity equivalent value Vsc
At a certain value, the combustion state becomes the best and the exhaust gas temperature decreases, and the exhaust gas temperature rises regardless of the deviation from Vsc.
Is set.

【0042】ステップ55では、前記基本排気温度Tex
hiを、ステップ52,53,54で求めた各排温補正係数K
tmpn,Ktmpp,Ktmpscを用いて図示の式
で補正することにより、シリンダ排出温度Texhiを
演算する。以下、ステップ56〜ステップ58については、
前記図9のステップ36〜ステップ38と同様であり、ステ
ップ56で、前記シリンダ排出温度Texhiに一次遅れ
処理を行って排気温度Texhとし、ステップ57で、図
示の式により排気圧基本値Pehbを演算し、ステップ
58では、前記排気圧基本値Pexhbに一次遅れ処理を
施して排気圧Pexhとする。
In step 55, the basic exhaust gas temperature Tex
hi is set to each exhaust temperature correction coefficient K obtained in steps 52, 53 and 54.
The cylinder discharge temperature Texhi is calculated by performing correction using the equations shown in the figure using tmpn, Ktmpp, and Ktmpsc. Hereinafter, regarding steps 56 to 58,
9 is the same as steps 36 to 38 in FIG. 9. In step 56, a first-order lag process is performed on the cylinder discharge temperature Texhi to obtain an exhaust temperature Texh. In step 57, a basic exhaust pressure value Pehb is calculated by the illustrated equation. Then step
In step 58, the exhaust pressure basic value Pexhb is subjected to a first-order lag process to obtain the exhaust pressure Pexh.

【0043】次に、吸気圧推定について説明する。図3
は、吸気圧推定の機能ブロック図を示し、図21は、吸気
圧を演算するフローであり、以下、図21のフローにした
がって説明する。ステップ61では、シリンダ吸入空気量
Qacと機関回転速度Neより、図22に示すテーブルか
ら体積効率相当値基本値Kinbを検索する。ステップ
62では、体積効率相当値Kinを、吸気温度Tintを
用いて図示の式 (Kin=Kinb×TA/Tint)
で演算する。
Next, the intake pressure estimation will be described. FIG.
Shows a functional block diagram of intake pressure estimation, and FIG.
This is a flow for calculating the pressure.
This will be explained. In step 61, a volume efficiency equivalent value basic value Kinb is retrieved from the table shown in FIG. 22 based on the cylinder intake air amount Qac and the engine rotation speed Ne. Steps
In 62, the volume efficiency equivalent value Kin is calculated using the intake air temperature Tint using the equation (Kin = Kinb × TA / Tint) shown in FIG.
Calculate with.

【0044】ステップ63では、シリンダ吸入空気量Qa
c,体積効率相当値Kinを変数として、図示の式 (P
m=Qac/Kin×TA×Ra/VCYL) で吸気圧
Pmを演算して処理を終了する。なお、このフローで出
てくる定数TA,RA,VCYLは各々標準状態温度,
空気ガス定数,シリンダ容積である。図23は、基本燃料
噴射量Qsol1を演算するフローである。
In step 63, the cylinder intake air amount Qa
c, using the volume efficiency equivalent value Kin as a variable, the equation (P
m = Qac / Kin × TA × Ra / VCYL), and the process ends. The constants TA, RA, and VCYL appearing in this flow are the standard state temperature,
Air gas constant and cylinder volume. FIG. 23 is a flowchart for calculating the basic fuel injection amount Qsol1.

【0045】ステップ71では、機関回転速度Ne、ステ
ップ72でアクセル開度 (コントロールレバー開度) CL
を読み込む。ステップ73では、機関回転速度Neとアク
セル開度CLとから、図24に示すテーブルから基本燃料
噴射量Mqdrvを検索する。ステップ74では、前記基
本燃料噴射量Mqdrvを水温等の各種補正係数によっ
て補正して基本燃料噴射量Qsol1を求める。
In step 71, the engine speed Ne is set, and in step 72, the accelerator opening (control lever opening) CL
Read. In step 73, the basic fuel injection amount Mqdrv is retrieved from the table shown in FIG. 24 based on the engine speed Ne and the accelerator opening CL. In step 74, the basic fuel injection amount Mqdrv is corrected by various correction coefficients such as a water temperature to obtain a basic fuel injection amount Qsol1.

【0046】図25は、最大燃料噴射量Qfulを演算す
るフローで、回転に同期したタイミングで処理される。
ステップ81では、機関回転速度Neを読み込み、ステッ
プ82では、機関回転速度Neから図26に示すスモークリ
ミット相当の値が設定されたテーブルから、限界空気過
剰率相当値Klambを検索する。
FIG. 25 is a flowchart for calculating the maximum fuel injection amount Qful, which is processed at a timing synchronized with the rotation.
In step 81, the engine rotational speed Ne is read, and in step 82, a limit excess air ratio equivalent value Klamb is retrieved from the table in which values corresponding to the smoke limit shown in FIG. 26 are set from the engine rotational speed Ne.

【0047】ステップ83では、シリンダ吸入空気量Qa
cを読み込み、ステップ84で図示の式 (Qful=Qa
c/Klamb/14.7) で最大燃料噴射量Qfulを演
算して、処理を終了する。図26は、燃料噴射量を最終的
に設定するフローであり、機関回転同期で実行される。
In step 83, the cylinder intake air amount Qa
c, and at step 84 the equation (Qfull = Qa)
c / Klamb / 14.7) is used to calculate the maximum fuel injection amount Qful, and the process ends. FIG. 26 is a flowchart for finally setting the fuel injection amount, which is executed in synchronization with the engine rotation.

【0048】ステップ91では、前記基本燃料噴射量Qs
ol1と最大燃料噴射量Qfulとを比較し、前者が大
のときはステップ92へ進み、燃料噴射量Qsolに最大
燃料噴射量Qfulを用い、前者が小のときはステップ
93へ進んで燃料噴射量Qsolに基本燃料噴射量Qso
l1を設定して処理を終了する。以下では、本発明によ
り推定した排気圧と吸気圧とを用いて行われるEGR制
御について説明する。
In step 91, the basic fuel injection amount Qs
ol1 is compared with the maximum fuel injection amount Qful, and when the former is large, the routine proceeds to step 92, where the maximum fuel injection amount Qful is used as the fuel injection amount Qsol, and when the former is small, the step is performed.
Proceed to 93 to add the basic fuel injection amount Qso to the fuel injection amount Qsol.
l1 is set and the process ends. Hereinafter, the EGR control performed using the exhaust pressure and the intake pressure estimated according to the present invention will be described.

【0049】図28は、EGR弁の開口面積を演算するフ
ローである。ステップ101 では、目標のEGR量を演算
する。この方法については、後述する。ステップ102 で
は、前記推定した吸気圧Pmを読み込み、ステップ103
では同じく推定した排気圧Pexhを読み込む。
FIG. 28 is a flowchart for calculating the opening area of the EGR valve. In step 101, a target EGR amount is calculated. This method will be described later. At step 102, the estimated intake pressure Pm is read, and at step 103
Then, the estimated exhaust pressure Pexh is read.

【0050】ステップ104 では、図示の式[Cqe=
{K (Pexh−Pm) }1/2 ;Kは定数]で、前記E
GR流速相当値Cqeを演算し、ステップ105 で、要求
EGR量Tqekと前記EGR流速相当値Cqeを用
い、図示の式 (Aev=Tqek/Cqe) でEGR弁
の開口面積Aevを演算する。ステップ106 では、前記
EGR流速相当値Cqeにより、図29に示すテーブルか
ら、後述するEGR弁開口面積を加重平均処理する際の
重み定数Nlkを検索する。なお、EGR流速相当値C
qeが小さいときに重み定数Nlkを大きく、EGR流
速相当値Cqeが大きいときに重み定数Nlkを小さく
設定している。これは、流速が小さいときは微小な流速
変化でも要求の開口面積は大きく変化させる必要があ
り、安定しづらくなるため加重平均を多く (重く) して
いる。また、流速が大きいときは逆の現象となり、一般
的に過渡時は流速が大きくなる (吸気圧と排気圧の差圧
が大きくなる) ため、過渡の追従性から加重平均はなる
べく行わない方が望ましいため、加重平均定数は小さく
なるように設定している。なお、反比例のような特性と
したのは、ステップ104 で図示した式に示されるよう
に、差圧に対して流速は平方根の特性を持ち、加重平均
定数の要求は、その逆数になるためである。
In step 104, the equation [Cqe =
{K (Pexh-Pm)} 1/2 ; K is a constant]
The GR flow rate equivalent value Cqe is calculated, and in step 105, the required EGR amount Tqek and the EGR flow rate equivalent value Cqe are used to calculate the opening area Aev of the EGR valve by the illustrated equation (Aev = Tqek / Cqe). In step 106, a weight constant Nlk for performing a weighted averaging process on an EGR valve opening area, which will be described later, is searched from the table shown in FIG. 29 based on the EGR flow velocity equivalent value Cqe. The EGR flow rate equivalent value C
When qe is small, the weight constant Nlk is set large, and when the EGR flow velocity equivalent value Cqe is large, the weight constant Nlk is set small. This is because when the flow velocity is small, the required opening area must be largely changed even with a small change in the flow velocity, and the weighted average is increased (heavy) because it becomes difficult to stabilize. In addition, when the flow velocity is high, the opposite phenomenon occurs.In general, the flow velocity increases during a transient (the differential pressure between the intake pressure and the exhaust pressure increases). For this reason, the weighted average constant is set to be small. The reason why the characteristic is inversely proportional is that, as shown in the equation shown in step 104, the flow velocity has a square root characteristic with respect to the differential pressure, and the requirement of the weighted average constant is the reciprocal thereof. is there.

【0051】ステップ107 では、前記ステップ106 で検
索した重み定数Nlkを用いて、次式によりステップ10
5 で求められた開口面積Aevに加重平均を行い、その
結果を目標のEGR弁開口面積Aevfとして、処理を
終了する。 Aevf=Aev/2Nlk + (1−1/2Nlk ) ×Ae
vfn-1 なお、ステップ104,ステップ105 の式は理論状態であ
る。実際のEGR弁駆動装置には、目標開口面積Aev
fを、図30に示すアクチュエータ特性で変換して指令値
となる。
In step 107, using the weight constant Nlk retrieved in step 106, step 10
The weighted average is performed on the opening area Aev obtained in step 5, and the result is set as the target EGR valve opening area Aevf, and the process is terminated. Aevf = Aev / 2 Nlk + (1-1 / 2 Nlk ) × Ae
vf n-1 The equations in steps 104 and 105 are theoretical states. The actual EGR valve driving device has a target opening area Aev
f is converted into a command value by the actuator characteristics shown in FIG.

【0052】図31は、目標のEGR量を演算するフロー
であり、機関回転又はそれ相当に同期したタイミングで
演算する。ステップ111 では、シリンダ吸入空気量Qa
cを読み込む。ステップ112 では、目標のEGR率Me
grを演算する。ステップ113 では、図示の式 (Mqe
c=Qac×Megr) で1吸気当りの目標のEGR量
Mqecを演算する。
FIG. 31 is a flowchart for calculating the target EGR amount, which is calculated at a timing synchronized with the engine rotation or equivalent. In step 111, the cylinder intake air amount Qa
Read c. In step 112, the target EGR rate Me
Calculate gr. In step 113, the equation (Mqe
The target EGR amount Mqec per intake air is calculated by c = Qac × Megr).

【0053】ステップ114 では、中間変数Rqecを次
式により演算する。 Rqec=Mqec×KIN×KVOl+Rqec
n-1(1−KIN×KVOl) 但し、KINは体積効率相当値, KVOl=VE/NC
/VMであり、VEは排気量,NCは気筒数,VMは吸
気系容積である。ステップ115 では、次式により進み補
正処理を行い、その結果をTqecとする。この式は、
通常の進み処理を簡易化したものである。
In step 114, the intermediate variable Rqec is calculated by the following equation. Rqec = Mqec × KIN × KVOL + Rqec
n-1 (1−KIN × KVOL) where KIN is a value corresponding to volumetric efficiency , KVOL = VE / NC
/ VM, VE is the displacement, NC is the number of cylinders, and VM is the intake system volume. In step 115, advance correction processing is performed by the following equation, and the result is set as Tqec. This expression is
This is a simplified version of the normal advance processing.

【0054】Tqec=GKQEC×Mqec− (GK
QEC−1) ×Rqecn-1 ステップ116 では、前記進み処理後の目標EGR量Tq
ecを、次式により単位時間当りの目標EGR量Tqe
kに変換して処理を終了する。 Tqek=Tqec×Ne/KCON 図32は、目標EGR率Megrを演算するフローで、機
関回転に同期したタイミングで実行される。
Tqec = GKQEC × Mqec− (GK
QEC-1) × Rqec n-1 In step 116, the target EGR amount Tq after the advance processing is performed.
ec is calculated by the following equation as the target EGR amount Tqe per unit time.
After that, the process is terminated. Tqek = Tqec × Ne / KCON FIG. 32 is a flow for calculating the target EGR rate Megr, which is executed at a timing synchronized with the engine rotation.

【0055】ステップ121 では、機関回転速度Ne, 燃
料噴射量Qsol, 機関冷却水温度Twを読み込む。ス
テップ122 では、機関回転速度Neと燃料噴射量Qso
lとにより、図33に示すテーブルから、基本目標EGR
率Megrを演算する。ステップ123 では、機関冷却水
温度Twにより、図34に示すテーブルから目標EGR率
Megrを補正するための補正係数Kegr twとす
る。
In step 121, the engine speed Ne, the fuel injection amount Qsol, and the engine coolant temperature Tw are read. In step 122, the engine speed Ne and the fuel injection amount Qso
and the basic target EGR from the table shown in FIG.
Calculate the rate Megr. In step 123, a correction coefficient Kegr for correcting the target EGR rate Megr from the table shown in FIG. 34 based on the engine coolant temperature Tw. tw.

【0056】ステップ124 では、図示の式 (Megr=
Megrb×Kegr tw) により、目標EGR率M
egrを演算する。ステップ125 では、機関の状態が完
爆状態か否かを判定する。この方法は、後述する。ステ
ップ126 で完爆と判定されたときは、そのまま処理を終
了し、完爆と判定されない場合は、ステップ127 へ進
み、目標EGR率Megrを0として処理を終了する。
この場合、同時にスロットル弁31は全開に制御される。
In step 124, the equation (Megr =
Megrb × Kegr tw), the target EGR rate M
Calculate egr. In step 125, it is determined whether or not the state of the engine is a complete explosion state. This method will be described later. If it is determined in step 126 that the explosion is complete, the process is terminated as it is. If it is not determined that the explosion is complete, the process proceeds to step 127, where the target EGR rate Megr is set to 0 and the process is terminated.
In this case, the throttle valve 31 is simultaneously controlled to be fully opened.

【0057】図35は、機関の完爆を判定するフローで、
10msec等の時間に同期したタイミングで実行され
る。ステップ131 では、機関回転速度Neを読み込み、
ステップ132 では、完爆判定スライスレベルNRPMK
と比較し、機関回転速度Neの方が大きいときには、ス
テップ133 へ進む。
FIG. 35 is a flow chart for judging a complete explosion of the engine.
It is executed at a timing synchronized with a time such as 10 msec. In step 131, the engine speed Ne is read, and
In step 132, the complete explosion determination slice level NRPMK
If the engine rotation speed Ne is higher than the above, the routine proceeds to step 133.

【0058】ステップ133 では、機関回転速度Neによ
る完爆判定後のカウンタの値Tmrkbと所定時間TM
RKBPとを比較し、Tmrkbの方が大きいときはス
テップ134 へ進み、完爆状態であるとして処理を終了す
る。ステップ132 で機関回転速度Neの方が小さいと判
定されたときは、ステップ136 へ進んでTmrkbをク
リアし、ステップ137 で完爆状態でないとして処理を終
了する。
In step 133, the counter value Tmrkb after the complete explosion determination based on the engine speed Ne and the predetermined time TM
RKBP is compared, and if Tmrkb is larger, the routine proceeds to step 134, where it is determined that the explosion is complete and the processing is terminated. If it is determined in step 132 that the engine rotation speed Ne is lower, the routine proceeds to step 136, where Tmrkb is cleared, and in step 137, the process is terminated because it is not in the complete explosion state.

【0059】ステップ133 でTmrkbの方が小さいと
判定されたときは、ステップ135 へ進んでTmrkbを
インクリメントした後ステップ137 へ進み、完爆状態で
ないとして処理を終了する。当処理では、機関回転速度
Neが所定以上 (例えば400 回転以上) となり、所定時
間経過したときに完爆状態であると判定するという処理
を行っている。
If it is determined in step 133 that Tmrkb is smaller, the routine proceeds to step 135, where Tmrkb is incremented. Then, the routine proceeds to step 137, where it is determined that the state is not the complete explosion state, and the processing is terminated. In this process, when the engine speed Ne becomes equal to or higher than a predetermined value (for example, 400 rotations or more), and when a predetermined time has elapsed, it is determined that the combustion state is the complete explosion state.

【0060】以上説明してきたように、本発明によれ
ば、エアフローメータを始めとする信号から排気圧を
算することにより、コストアップを招くことなく耐久
性, 過渡応答性を補償する圧力推定が可能となり、更に
噴射時期, 着火時期やスワール制御弁開度, 開口面積に
対する吸気の流速の比で排気圧を補正することにより、
推定精度が向上し、EGR制御等の制御精度が向上す
る。
As described above, according to the present invention, by calculating the exhaust pressure from a signal such as an air flow meter, the durability and the transient response can be improved without increasing the cost. By compensating the pressure estimation, it is possible to correct the exhaust pressure based on the injection timing, ignition timing, swirl control valve opening, and the ratio of the intake flow velocity to the opening area.
The estimation accuracy is improved, and the control accuracy of EGR control and the like is improved.

【図面の簡単な説明】[Brief description of the drawings]

【図1】請求項1に係る発明の構成・機能を示すブロッ
ク図。
FIG. 1 is a block diagram showing the configuration and functions of the invention according to claim 1 ;

【図2】請求項2〜請求項4に係る発明の構成・機能を
示すブロック図。
FIG. 2 is a block diagram showing the configuration and functions of the invention according to claims 2 to 4 ;

【図3】本発明の一実施形態における吸気圧推定の機能
ブロック図。
FIG. 3 is a function of an intake pressure estimation according to an embodiment of the present invention;
Block Diagram.

【図4】本発明の一実施形態のシステム構成を示す図。FIG. 4 is a diagram showing a system configuration according to an embodiment of the present invention.

【図5】吸入空気流量を演算するルーチンのフローチャ
ート。
FIG. 5 is a flowchart of a routine for calculating an intake air flow rate.

【図6】エアフローメータ出力電圧と吸入空気流量との
関係を示す線図。
FIG. 6 is a diagram showing a relationship between an output voltage of an air flow meter and an intake air flow rate .

【図7】同じくシリンダ吸入空気量を演算するルーチン
のフローチャート。
FIG. 7 is a flowchart of a routine for calculating a cylinder intake air amount.

【図8】同じくシリンダサイクル遅れ処理を示すフロー
チャート。
FIG. 8 is a flowchart showing a cylinder cycle delay process.

【図9】同じく排気系圧力を演算するルーチンのフロー
チャート。
FIG. 9 is a flowchart of a routine for calculating the exhaust system pressure.

【図10】同じく燃料噴射量による排温補正係数を設定し
たテーブル。
FIG. 10 is a table in which an exhaust temperature correction coefficient according to a fuel injection amount is set.

【図11】同じく吸入空気温度による排温補正係数を設定
したテーブル。
FIG. 11 is a table in which exhaust temperature correction coefficients according to the intake air temperature are set.

【図12】同じく前回の排気圧推定値による排温補正係数
を設定したテーブル。
FIG. 12 is a table in which an exhaust gas temperature correction coefficient based on a previous exhaust pressure estimation value is set.

【図13】同じく燃料噴射時期による排温補正係数を設定
したテーブル。
FIG. 13 is a table in which exhaust temperature correction coefficients according to fuel injection timing are set.

【図14】本発明の第2の実施形態のシステム構成を示す
図。
FIG. 14 is a diagram illustrating a system configuration according to a second embodiment of the present invention.

【図15】同じくシリンダサイクル遅れ処理を示すフロー
チャート。
FIG. 15 is a flowchart showing a cylinder cycle delay process.

【図16】同じく排気系圧力を演算するルーチンのフロー
チャート。
FIG. 16 is a flowchart of a routine for calculating the exhaust system pressure.

【図17】同じくスワール制御状態による排温補正係数を
設定したテーブル。
FIG. 17 is a table in which exhaust temperature correction coefficients according to the swirl control state are set.

【図18】同じくシリンダサイクル遅れ処理を示すフロー
チャート。
FIG. 18 is a flowchart showing a cylinder cycle delay process.

【図19】別の実施形態における排気系圧力を演算するル
ーチンのフローチャート。
FIG. 19 is a flowchart of a routine for calculating an exhaust system pressure in another embodiment.

【図20】同じくスワール制御状態による排温補正係数を
設定したテーブル。
FIG. 20 is a table in which exhaust temperature correction coefficients according to the swirl control state are set.

【図21】本発明の別実施形態に係る吸気圧力を演算する
ルーチンのフローチャート。
FIG. 21 is a flowchart of a routine for calculating intake pressure according to another embodiment of the present invention.

【図22】同じく体積効率相当値基本値を設定したテーブ
ル。
FIG. 22 is a table in which basic values for volume efficiency equivalent values are set.

【図23】同じく基本燃料噴射量を演算するルーチンのフ
ローチャート。
FIG. 23 is a flowchart of a routine for calculating a basic fuel injection amount.

【図24】同じく基本燃料噴射量を設定したテーブル。FIG. 24 is a table in which basic fuel injection amounts are set.

【図25】同じく最大燃料噴射量を演算するルーチンのフ
ローチャート。
FIG. 25 is a flowchart of a routine for calculating a maximum fuel injection amount.

【図26】同じく限界空気過剰率を設定したテーブル。FIG. 26 is also a table in which a limit excess air ratio is set.

【図27】同じく燃料噴射量を設定するルーチンのフロー
チャート。
FIG. 27 is a flowchart of a routine for similarly setting a fuel injection amount.

【図28】同じくEGR弁開口面積を演算するルーチンの
フローチャート。
FIG. 28 is a flowchart of a routine for calculating an EGR valve opening area in the same manner.

【図29】同じくEGR弁開口面積を加重平均演算するた
めの重み定数を設定したテーブル。
FIG. 29 is a table in which weight constants for calculating a weighted average of the EGR valve opening area are set.

【図30】同じく目標開口面積とステップモータのステッ
プ数の関係を示すテーブル。
FIG. 30 is a table showing the relationship between the target opening area and the number of steps of the step motor.

【図31】同じく目標EGR量を設定するルーチンのフロ
ーチャート。
FIG. 31 is a flowchart of a routine for similarly setting a target EGR amount.

【図32】同じく目標EGR率を設定するルーチンのフロ
ーチャート。
FIG. 32 is a flowchart of a routine for similarly setting a target EGR rate.

【図33】同じく目標EGR率を設定したテーブル。FIG. 33 is a table in which target EGR rates are also set.

【図34】同じく目標EGR率の水温補正係数を設定した
テーブル。
FIG. 34 is a table in which a water temperature correction coefficient of the target EGR rate is set.

【図35】同じく完爆判定ルーチンを示すフローチャー
ト。
FIG. 35 is a flowchart showing a complete explosion determination routine.

【符号の説明】[Explanation of symbols]

5 ディーゼル機関 6 燃料噴射ノズル 7 燃料噴射ポンプ 11 バキュームポンプ 13 コントロールユニット 16 エアフローメータ 17 回転速度センサ 18 レバー開度センサ 19 水温センサ 21 スワール制御弁 31 スロットル弁 34 リフトセンサ Reference Signs List 5 diesel engine 6 fuel injection nozzle 7 fuel injection pump 11 vacuum pump 13 control unit 16 air flow meter 17 rotation speed sensor 18 lever opening sensor 19 water temperature sensor 21 swirl control valve 31 throttle valve 34 lift sensor

───────────────────────────────────────────────────── フロントページの続き (56)参考文献 特開 平6−229324(JP,A) 特開 平7−158479(JP,A) 特開 平5−180057(JP,A) 特開 平7−332139(JP,A) 特開 平2−157438(JP,A) 特開 平8−100695(JP,A) 特開 平8−121233(JP,A) 特開 平1−315635(JP,A) 特開 平4−153547(JP,A) 特開 平1−271643(JP,A) 特開 平8−284735(JP,A) 特開 昭64−53032(JP,A) 特開 平9−170480(JP,A) 特公 平7−13454(JP,B2) (58)調査した分野(Int.Cl.7,DB名) F02D 41/00 - 45/00 F02D 29/00 - 29/06 F02D 9/00 - 28/00 F01N 9/00 - 11/00 F01N 3/00 - 3/38 F02B 31/00 - 31/02 ──────────────────────────────────────────────────続 き Continuation of front page (56) References JP-A-6-229324 (JP, A) JP-A-7-158479 (JP, A) JP-A-5-180057 (JP, A) JP-A-7- 332139 (JP, A) JP-A-2-157438 (JP, A) JP-A-8-10069 (JP, A) JP-A-8-121233 (JP, A) JP-A-1-315635 (JP, A) Japanese Unexamined Patent Publication No. Hei 4-153547 (JP, A) Japanese Unexamined Patent Publication No. Hei 1-27-16473 (JP, A) Japanese Unexamined Patent Publication No. Hei 8-284735 (JP, A) Japanese Unexamined Patent Publication No. Sho 64-53032 (JP, A) Japanese Unexamined Patent Publication No. Hei 9-170480 (JP, A) JP 7-13454 (JP, B2) (58) Fields investigated (Int. Cl. 7 , DB name) F02D 41/00-45/00 F02D 29/00-29/06 F02D 9 / 00-28/00 F01N 9/00-11/00 F01N 3/00-3/38 F02B 31/00-31/02

Claims (4)

(57)【特許請求の範囲】(57) [Claims] 【請求項1】機関に噴射される燃料の噴射量を検出する
燃料噴射量検出手段と、 機関に吸入される空気の温度を検出する吸入空気温度検
出手段と、 前記燃料の噴射時期を検出する燃料噴射時期検出手段
と、 前記吸入空気の流量を検出する吸入空気流量検出手段
と、 機関の回転速度を検出する機関回転速度検出手段と、燃料噴射量に基づいて基本排気温度を演算し、該基本排
気温度を吸入空気温度と前回推定された排気圧力と燃料
噴射時期とで補正してシリンダからの排気温度を演算
し、吸入空気流量と機関回転速度とから求めたシリンダ
からの排気流量と前記シリンダからの排気温度とに基づ
いて排気圧を推定する 排気系圧力推定手段と、 を含んで構成したことを特徴とする内燃機関の排気系圧
力推定装置。
A fuel injection amount detector for detecting an amount of fuel injected into the engine; an intake air temperature detector for detecting a temperature of air taken into the engine; and detecting an injection timing of the fuel. Fuel injection timing detecting means, intake air flow rate detecting means for detecting the flow rate of the intake air, engine rotational speed detecting means for detecting the rotational speed of the engine , and calculating a basic exhaust temperature based on the fuel injection amount; Basic exhaust
The intake air temperature and the exhaust pressure and fuel estimated previously
Calculate exhaust temperature from cylinder by correcting with injection timing
And the cylinder determined from the intake air flow rate and the engine speed
Based on the exhaust gas flow from the cylinder and the exhaust gas temperature from the cylinder.
An exhaust system pressure estimating device for an internal combustion engine, comprising: exhaust system pressure estimating means for estimating exhaust pressure.
【請求項2】 機関に噴射される燃料の噴射量を検出する
燃料噴射量検出手段と、 機関に吸入される空気の温度を検出する吸入空気温度検
出手段と、 機関の吸気系に備えられたスワール制御弁によるスワー
ル制御状態を検出するスワール制御状態検出手段と、 前記吸入空気の流量を検出する吸入空気流量検出手段
と、 機関の回転速度を検出する機関回転速度検出手段と、燃料噴射量に基づいて基本排気温度を演算し、該基本排
気温度を吸入空気温度と前回推定された排気圧力とスワ
ール制御状態とで補正してシリンダからの排気温度を演
算し、吸入空気流量と機関回転速度とから求めたシリン
ダからの排気流量と前記シリンダからの排気温度とに基
づいて排気圧を推定する 排気系圧力推定手段と、を含ん
で構成したことを特徴とする内燃機関の排気系圧力推定
装置。
And the fuel injection amount detecting means for detecting the amount of fuel injected to 2. A engine, and intake air temperature detecting means for detecting a temperature of air sucked into the engine, provided in an intake system of the engine a swirl control condition detecting means for detecting the swirl control state by the swirl control valve, and the intake air flow rate detecting means for detecting the flow rate of the intake air, and the engine rotational speed detecting means for detecting a rotational speed of the engine, the fuel injection amount The basic exhaust temperature is calculated based on the basic exhaust temperature.
The intake air temperature, the exhaust pressure and the sw
To control the exhaust temperature from the cylinder.
Calculated from the intake air flow rate and the engine speed.
And the exhaust temperature from the cylinder.
An exhaust system pressure estimating means for estimating the exhaust pressure based on the exhaust system pressure.
【請求項3】 前記基本排気温度のスワール制御状態によ
る補正は、スワール制御弁の開度と機関回転速度とに基
づいて設定される補正係数により行うことを特徴とする
請求項2に記載の内燃機関の排気系圧力推定装置。
3. The method according to claim 1, wherein the correction of the basic exhaust gas temperature based on the swirl control state is performed by a correction coefficient set based on the opening degree of the swirl control valve and the engine speed.
An exhaust system pressure estimating device for an internal combustion engine according to claim 2 .
【請求項4】 前記基本排気温度のスワール制御状態によ
る補正は、吸入空気流量と機関回転速度とスワール制御
弁の開度とから演算したスワール流速相当値に基づいて
設定される補正係数により行うことを特徴とする請求項
に記載の内燃機関の排気系圧力推定装置。
Wherein correction by the swirl control state of the basic exhaust gas temperature is carried out by the correction coefficient which is set based on the swirl flow velocity equivalent value calculated from the degree of opening of the intake air flow rate and the engine rotational speed and the swirl control valve Claims characterized by the following :
3. The exhaust system pressure estimation device for an internal combustion engine according to 2 .
JP17891996A 1996-06-03 1996-07-09 Exhaust system pressure estimation device for internal combustion engine Expired - Lifetime JP3358449B2 (en)

Priority Applications (8)

Application Number Priority Date Filing Date Title
JP17891996A JP3358449B2 (en) 1996-07-09 1996-07-09 Exhaust system pressure estimation device for internal combustion engine
GB9825391A GB2329040B (en) 1996-06-03 1997-06-03 Apparatus for estimating pressure in intake system and exhaust system of internal combustion engine
DE19758641A DE19758641B4 (en) 1996-06-03 1997-06-03 Control apparatus for exhaust gas recirculation in internal combustion engine - calculates target EGR rate according to standard rate and required amount of adjustment which is determined according to standard injection volume, target and actual injection timing
GB9711483A GB2313927B (en) 1996-06-03 1997-06-03 EGR control apparatus for internal combustion engine
US08/867,980 US6012431A (en) 1996-06-03 1997-06-03 Control apparatus for internal combustion engine and estimation apparatus for estimating pressure in intake and discharge system of internal combustion engine
DE19723210A DE19723210B4 (en) 1996-06-03 1997-06-03 Control apparatus for exhaust gas recirculation in internal combustion engine - calculates target EGR rate according to standard rate and required amount of adjustment which is determined according to standard injection volume, target and actual injection timing
US09/427,004 US6167342A (en) 1996-06-03 1999-10-26 Control apparatus for internal combustion engine and estimation apparatus for estimating pressure in intake and discharge system of internal combustion engine
US09/699,342 US6298299B1 (en) 1996-06-03 2000-10-31 Control apparatus for internal combustion engine and estimation apparatus for estimating pressure in intake and discharge system of internal combustion engine

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP17891996A JP3358449B2 (en) 1996-07-09 1996-07-09 Exhaust system pressure estimation device for internal combustion engine

Related Child Applications (1)

Application Number Title Priority Date Filing Date
JP2002127146A Division JP3395782B2 (en) 2002-04-26 2002-04-26 Intake system pressure estimation device for internal combustion engine

Publications (2)

Publication Number Publication Date
JPH1026049A JPH1026049A (en) 1998-01-27
JP3358449B2 true JP3358449B2 (en) 2002-12-16

Family

ID=16056949

Family Applications (1)

Application Number Title Priority Date Filing Date
JP17891996A Expired - Lifetime JP3358449B2 (en) 1996-06-03 1996-07-09 Exhaust system pressure estimation device for internal combustion engine

Country Status (1)

Country Link
JP (1) JP3358449B2 (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3856118B2 (en) * 2002-01-31 2006-12-13 日産自動車株式会社 Exhaust purification device
DE102007033678B4 (en) * 2007-07-19 2022-08-11 Robert Bosch Gmbh Method and device for controlling an internal combustion engine
WO2012032627A1 (en) * 2010-09-08 2012-03-15 トヨタ自動車株式会社 Control device for internal combustion engine
JP5929015B2 (en) 2011-06-06 2016-06-01 日産自動車株式会社 Exhaust gas recirculation device for internal combustion engine
JP5373952B1 (en) 2012-09-04 2013-12-18 日野自動車株式会社 Exhaust temperature estimation device for internal combustion engine

Also Published As

Publication number Publication date
JPH1026049A (en) 1998-01-27

Similar Documents

Publication Publication Date Title
US6318083B1 (en) Intake air control device of an engine with a charger and method thereof
US6298299B1 (en) Control apparatus for internal combustion engine and estimation apparatus for estimating pressure in intake and discharge system of internal combustion engine
US9261031B2 (en) Control device for internal combustion engine and method for controlling internal combustion engine
US20050172628A1 (en) Boost pressure estimation apparatus for internal combustion engine with supercharger
US7681442B2 (en) Throttle upstream pressure estimating apparatus and cylinder charged air quantity calculating apparatus for internal combustion engine
JP6381728B1 (en) Control device for internal combustion engine
US6019094A (en) Engine exhaust recirculation device
US6640775B2 (en) Air-fuel ratio control system for internal combustion engine
JP3888024B2 (en) Exhaust gas recirculation device
EP1245818B1 (en) Air-fuel ratio control apparatus and method for internal combustion engine
JP3358449B2 (en) Exhaust system pressure estimation device for internal combustion engine
JP3493981B2 (en) Supercharging pressure control device for internal combustion engine with EGR control device
JP2001073789A (en) Supercharging pressure control system for internal combustion engine
JP3395782B2 (en) Intake system pressure estimation device for internal combustion engine
JP3237480B2 (en) Exhaust gas recirculation control device for internal combustion engine
JP3629752B2 (en) Intake and exhaust system pressure estimation device for internal combustion engine
JPH11182356A (en) Egr controller for internal combustion engine
EP1683953B1 (en) Control apparatus for an internal combustion engine
JPH11200926A (en) Fuel injection control device for internal combustion engine
JP3460394B2 (en) Exhaust gas recirculation control device for internal combustion engine
JP2001295702A (en) Exhaust gas recirculation device
JP3704914B2 (en) Torque control device for diesel engine with EGR control device
JP3622506B2 (en) EGR control device for internal combustion engine
JPH10103162A (en) Egr control device for diesel engine
JP3374619B2 (en) Exhaust gas recirculation control device for internal combustion engine

Legal Events

Date Code Title Description
FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20081011

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20091011

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20101011

Year of fee payment: 8

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20111011

Year of fee payment: 9

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20121011

Year of fee payment: 10

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20121011

Year of fee payment: 10

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20131011

Year of fee payment: 11

EXPY Cancellation because of completion of term