JPH1182086A - Method for controlling air-fuel ratio of internal combustion engine - Google Patents
Method for controlling air-fuel ratio of internal combustion engineInfo
- Publication number
- JPH1182086A JPH1182086A JP24884397A JP24884397A JPH1182086A JP H1182086 A JPH1182086 A JP H1182086A JP 24884397 A JP24884397 A JP 24884397A JP 24884397 A JP24884397 A JP 24884397A JP H1182086 A JPH1182086 A JP H1182086A
- Authority
- JP
- Japan
- Prior art keywords
- fuel ratio
- cylinder
- air
- correction coefficient
- torque fluctuation
- 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.)
- Pending
Links
Landscapes
- Control Of Driving Devices And Active Controlling Of Vehicle (AREA)
- Control Of Vehicle Engines Or Engines For Specific Uses (AREA)
- Electrical Control Of Air Or Fuel Supplied To Internal-Combustion Engine (AREA)
- Combined Controls Of Internal Combustion Engines (AREA)
Abstract
Description
【0001】[0001]
【発明の属する技術分野】本発明は内燃機関の空燃比制
御方法に関する。The present invention relates to a method for controlling an air-fuel ratio of an internal combustion engine.
【0002】[0002]
【従来の技術】従来より、内燃機関のトルク変動量を検
出し、このトルク変動量が目標トルク変動量となるよう
に空燃比を制御する内燃機関が知られている。即ち、空
燃比がリーンになるほど機関で発生するNOX 量が減少
し、しかも燃料消費率が小さくなるので空燃比をできる
だけリーンにするのが好ましい。ところが空燃比が或る
程度以上リーンになると燃焼が不安定となり、その結果
トルク変動量が大きくなって車両振動が増大し、斯くし
て車両の運転性が悪化する。そこでこの内燃機関では、
排出NOX 量および燃料消費率をできるだけ低減し、同
時に車両の運転性の悪化をできるだけ低減しうるトルク
変動量領域内に目標トルク変動量を定め、トルク変動量
が目標トルク変動量となるように空燃比を制御してい
る。2. Description of the Related Art There has been known an internal combustion engine which detects a torque fluctuation amount of an internal combustion engine and controls an air-fuel ratio so that the torque fluctuation amount becomes a target torque fluctuation amount. That will reduce the amount of NO X produced by the engine as the air fuel ratio becomes lean, yet preferably as much as possible a lean air-fuel ratio because the fuel consumption rate is reduced. However, when the air-fuel ratio becomes lean to a certain degree or more, combustion becomes unstable, and as a result, the amount of torque fluctuation increases, vehicle vibration increases, and the drivability of the vehicle deteriorates. So in this internal combustion engine,
A target torque variation is set within a torque variation range in which the emission NO X amount and the fuel consumption rate are reduced as much as possible, and at the same time, the deterioration of the drivability of the vehicle is reduced as much as possible. The air-fuel ratio is controlled.
【0003】ところが、この内燃機関にロックアップ機
構を備えたトルクコンバータを適用すると次のような問
題を生ずる。即ち、ロックアップがオンのときのときに
はクランクシャフトとトルクコンバータの入力軸とが直
結状態とされるのでこのときトルク変動が大きくなると
車両の運転性が悪化しやすく、これに対しロックアップ
がオフのときにはクランクシャフトの回転力がポンプカ
バーやタービンランナなどを介してトルクコンバータの
入力軸に伝達されるのでこのときトルク変動が大きくな
っても車両の運転性はそれほど悪化しない。従って、ロ
ックアップがオフのときに最適となるように上述の目標
トルク変動量を定めるとロックアップがオンのときに車
両の運転性が悪化し、ロックアップがオンのときに最適
となるように上述の目標トルク変動量を定めるとロック
アップがオフのときに排出NOX量および燃料消費量を
十分に低減することができない。そこで、ロックアップ
がオンのときの目標トルク変動量をロックアップがオフ
のときの目標トルク変動量よりも小さく定めた空燃比制
御方法が公知である(特開平4−224245号公報参
照)。However, when a torque converter having a lock-up mechanism is applied to this internal combustion engine, the following problem occurs. That is, when lock-up is on, the crankshaft is directly connected to the input shaft of the torque converter. If the torque fluctuations are large at this time, the drivability of the vehicle is likely to deteriorate, whereas the lock-up is off. Sometimes, the rotational force of the crankshaft is transmitted to the input shaft of the torque converter via a pump cover, a turbine runner, or the like, so that the operability of the vehicle does not deteriorate so much even if the torque fluctuation increases. Therefore, if the target torque fluctuation amount is determined so as to be optimal when the lock-up is off, the drivability of the vehicle is deteriorated when the lock-up is on and is optimized when the lock-up is on. lockup when determining the target torque variation amount described above can not be sufficiently reduced emissions amount of nO X and the fuel consumption in the off. Therefore, there is known an air-fuel ratio control method in which the target torque fluctuation amount when lock-up is on is smaller than the target torque fluctuation amount when lock-up is off (see Japanese Patent Application Laid-Open No. 4-224245).
【0004】[0004]
【発明が解決しようとする課題】しかしながら、ロック
アップクラッチのスリップ制御時におけるトルク変動に
基づく車両振動はロックアップがオンのときよりも小さ
く、ロックアップがオフのときよりも大きい。従って、
例えばロックアップクラッチのスリップ制御時の目標ト
ルク変動量をロックアップがオフのときのための目標ト
ルク変動量とすると車両の運転性が悪化し、ロックアッ
プクラッチのスリップ制御時の目標トルク変動量をロッ
クアップがオンのときのための目標トルク変動量とする
と排出NO X 量および燃料消費量を十分に低減すること
ができないという問題点がある。云い換えるとロックア
ップクラッチのスリップ量、即ちクランクシャフトとト
ルクコンバータの入力軸との回転数差に応じて目標トル
ク変動量を定める必要がある。上述の空燃比制御方法は
このことについて何ら示唆していない。However, the lock
For torque fluctuation during slip control of up clutch
Based vehicle vibration is smaller than when lockup is on
Larger than when lockup is off. Therefore,
For example, the target torque during slip control of the lock-up clutch
Set the target torque for when the lockup is off.
If the amount of torque fluctuation is considered, the drivability of the vehicle is
Locks the target torque fluctuation during clutch slip control.
Target torque fluctuation for when the backup is on
And emission NO XReduce fuel consumption and fuel consumption
There is a problem that can not be. In other words, Rockua
Slip amount, that is, the crankshaft and
Target torque according to the rotational speed difference from the input shaft of the Luc converter.
It is necessary to determine the amount of fluctuation. The air-fuel ratio control method described above
Nothing is suggested about this.
【0005】[0005]
【課題を解決するための手段】上記問題点を解決するた
めに1番目の発明によれば、内燃機関の出力変動を検出
し、出力変動が目標出力変動となるように空燃比を制御
する内燃機関の空燃比制御方法において、ロックアップ
クラッチのスリップ量を検出し、このスリップ量に基づ
いて目標出力変動を補正するようにしている。即ち1番
目の発明によれば、ロックアップクラッチのスリップ量
に基づいて目標出力変動が補正されるのでスリップ量に
関わらず、即ちスリップ制御中にもスリップ制御中でな
いときにも車両の運転性の悪化が阻止されつつ燃料消費
率およびNOX の発生量が低減される。According to a first aspect of the present invention, there is provided an internal combustion engine for detecting an output fluctuation of an internal combustion engine and controlling an air-fuel ratio so that the output fluctuation becomes a target output fluctuation. In the air-fuel ratio control method of the engine, the slip amount of the lock-up clutch is detected, and the target output fluctuation is corrected based on the slip amount. That is, according to the first aspect, the target output fluctuation is corrected based on the slip amount of the lock-up clutch. while deterioration is prevented generation of the fuel consumption rate and NO X is reduced.
【0006】また、上記問題点を解決するために2番目
の発明によれば、空燃比を検出し、空燃比が目標リーン
空燃比となるように空燃比を制御する内燃機関の空燃比
制御方法において、ロックアップクラッチのスリップ量
を検出し、このスリップ量に基づいて目標リーン空燃比
を補正するようにしている。即ち2番目の発明では、ロ
ックアップクラッチのスリップ量に基づいて目標リーン
空燃比が補正されるのでスリップ量に関わらず、即ちス
リップ制御中にもスリップ制御中でないときにも車両の
運転性の悪化が阻止されつつ燃料消費率およびNOX の
発生量が低減される。According to a second aspect of the present invention, there is provided an air-fuel ratio control method for an internal combustion engine which detects an air-fuel ratio and controls the air-fuel ratio so that the air-fuel ratio becomes a target lean air-fuel ratio. , The slip amount of the lock-up clutch is detected, and the target lean air-fuel ratio is corrected based on the detected slip amount. That is, in the second aspect, the target lean air-fuel ratio is corrected based on the slip amount of the lock-up clutch. Is reduced, and the fuel consumption rate and the amount of generated NO X are reduced.
【0007】[0007]
【発明の実施の形態】図1を参照すると、1は1番気筒
#1、2番気筒#2、3番気筒#3、4番気筒#4から
なる4つの気筒を具備した機関本体を示す。各気筒#
1,#2,#3,#4は夫々対応する吸気枝管2を介し
てサージタンク3に連結され、各吸気枝管2内には夫々
対応する吸気ポート内に向って燃料を噴射する燃料噴射
弁4が取付けられる。サージタンク3は吸気ダクト5お
よびエアフローメータ6を介してエアクリーナ7に連結
され、吸気ダクト5内にはスロットル弁8が配置され
る。一方、各気筒#1,#2,#3,#4は排気マニホ
ルド9および排気管10を介してNOx 吸収剤11を内
蔵したケーシング12に連結される。このNOx 吸収剤
11は空燃比がリーンのときに排気ガス中に含まれるN
Ox を吸収し、空燃比が理論空燃比又はリッチになると
吸収したNOx を放出しかつ還元する機能を有する。DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS Referring to FIG. 1, reference numeral 1 denotes an engine main body having four cylinders including a first cylinder # 1, a second cylinder # 2, a third cylinder # 3, and a fourth cylinder # 4. . Each cylinder #
The fuel tanks # 1, # 2, # 3, and # 4 are connected to the surge tank 3 through the corresponding intake branch pipes 2 and inject fuel into the corresponding intake ports in the respective intake branch pipes 2. The injection valve 4 is attached. The surge tank 3 is connected to an air cleaner 7 via an intake duct 5 and an air flow meter 6, and a throttle valve 8 is arranged in the intake duct 5. On the other hand, the cylinders # 1, # 2, # 3, # 4 is connected to the casing 12 with a built-in the NO x absorbent 11 via an exhaust manifold 9 and an exhaust pipe 10. This NO x absorbent 11 contains N contained in exhaust gas when the air-fuel ratio is lean.
It has a function of absorbing O x and releasing and reducing the absorbed NO x when the air-fuel ratio becomes stoichiometric or rich.
【0008】電子制御ユニット20はディジタルコンピ
ュータからなり、双方向性バス21によって相互に接続
されたROM(リードオンリメモリ)22、RAM(ラ
ンダムアクセスメモリ)23、CPU(マイクロプロセ
ッサ)24、常時電源に接続されたバックアップRAM
25、入力ポート26および出力ポート27を具備す
る。機関の出力軸またはクランクシャフト13には外歯
付ロータ14が取付けられ、ロータ14の外歯に対面し
て電磁ピックアップからなるクランク角センサ15が配
置される。第1実施例ではロータ14の外周上に30°
クランク角度毎に外歯が形成されており、例えば1番気
筒の圧縮上死点を検出するために一部の外歯が削除され
ている。従ってこの外歯が削除された部分、即ち欠歯部
分を除いてクランク角センサ15は出力軸13が30°
クランク角度回転する毎に出力パルスを発生し、この出
力パルスが入力ポート26に入力される。The electronic control unit 20 is composed of a digital computer, and is connected to a ROM (read only memory) 22, a RAM (random access memory) 23, a CPU (microprocessor) 24, and a power supply connected to each other by a bidirectional bus 21. Backup RAM connected
25, an input port 26 and an output port 27. A rotor 14 with external teeth is attached to an output shaft or a crankshaft 13 of the engine. In the first embodiment, 30 °
External teeth are formed for each crank angle. For example, some external teeth are deleted to detect the compression top dead center of the first cylinder. Accordingly, the crank angle sensor 15 has an output shaft 13 of 30 ° except for the portion where the external teeth are deleted, that is, the tooth missing portion.
An output pulse is generated each time the crank angle rotates, and the output pulse is input to the input port 26.
【0009】エアフローメータ6は吸入空気量に比例し
た出力電圧を発生し、この出力電圧が対応するAD変換
器28を介して入力ポート26に入力される。また、ス
ロットル弁8にはスロットル弁8がアイドリング開度に
あることを検出するためのアイドルスイッチ16が取付
けられ、このアイドルスイッチ16の出力信号が入力ポ
ート26に入力される。また、排気マニホルド9内には
空燃比を検出するための空燃比センサ(O2 センサ)1
7が配置されており、この空燃比センサ17の出力信号
が対応するAD変換器28を介して入力ポート26に入
力される。一方、出力ポート27は対応する駆動回路2
9を介して各燃料噴射弁4に接続される。The air flow meter 6 generates an output voltage proportional to the amount of intake air, and this output voltage is input to an input port 26 via a corresponding AD converter 28. An idle switch 16 for detecting that the throttle valve 8 is at the idling opening is attached to the throttle valve 8, and an output signal of the idle switch 16 is input to an input port 26. An air-fuel ratio sensor (O 2 sensor) 1 for detecting the air-fuel ratio is provided in the exhaust manifold 9.
The output signal of the air-fuel ratio sensor 17 is input to the input port 26 via the corresponding AD converter 28. On the other hand, the output port 27 is connected to the corresponding drive circuit 2
9 is connected to each fuel injection valve 4.
【0010】図1に示されるように機関の出力軸13は
自動変速機30に連結され、自動変速機30の出力軸3
1が駆動輪に連結される。自動変速機30はトルクコン
バータ32を備えており、このトルクコンバータ32内
にロックアップ機構33が設けられている。即ち、トル
クコンバータ32は出力軸13に連結されて出力軸13
と共に回転するポンプカバー34と、ポンプカバー34
により支承されたポンプインペラ35と、自動変速機3
0の入力軸36に取付けられたタービンランナ37と、
ステータ37aとを具備し、出力軸13の回転運動がポ
ンプカバー34、ポンプインペラ35およびタービンラ
ンナ37を介して入力軸36に伝達される。As shown in FIG. 1, an output shaft 13 of the engine is connected to an automatic transmission 30, and an output shaft 3 of the automatic transmission 30 is provided.
1 is connected to the drive wheels. The automatic transmission 30 includes a torque converter 32, and a lock-up mechanism 33 is provided in the torque converter 32. That is, the torque converter 32 is connected to the output shaft 13 and
A pump cover 34 that rotates with the pump cover 34
Impeller 35 supported by the automatic transmission 3
A turbine runner 37 attached to the input shaft 36 of
The rotary motion of the output shaft 13 is transmitted to the input shaft 36 via the pump cover 34, the pump impeller 35, and the turbine runner 37.
【0011】一方、ロックアップ機構33は入力軸36
に対してその軸線方向に移動可能に取付けられかつ入力
軸36と共に回転するロックアップクラッチ板38を具
備する。通常は、即ちロックアップオフ時には入力軸3
6内のオイル通路を介してロックアップクラッチ板38
とポンプカバー34間の部屋39内に加圧オイルが供給
され、次いでこの部屋39から流出した加圧オイルはポ
ンプインペラ35およびタービンランナ37の周りの部
屋40内に送り込まれた後、入力軸36内のオイル通路
を介して排出される。このときロックアップクラッチ板
38両側の部屋39,40間の圧力差はほとんど生じな
いためにロックアップクラッチ板38はポンプカバー3
4の内壁面から離れており、従ってこのときにはクラン
クシャフト13の回転力はポンプカバー34、ポンプイ
ンペラ35およびタービンランナ37を介して入力軸3
6に伝達される。なお、このときのロックアップクラッ
チのスリップ量NSLP、即ちクランクシャフトの回転
数Nと自動変換器30の入力軸36即ちタービンランナ
37の回転数NTとの差(N−NT)はかなり大きくな
っている。On the other hand, the lock-up mechanism 33 has an input shaft 36.
And a lock-up clutch plate 38 movably mounted in the axial direction with respect to the input shaft 36 and rotating together with the input shaft 36. Normally, that is, when the lock-up is off, the input shaft 3
6 through the oil passage in the lock-up clutch plate 38.
The pressurized oil is supplied into a room 39 between the pump cover 34 and the pump cover 34, and then the pressurized oil flowing out of the room 39 is fed into a room 40 around the pump impeller 35 and the turbine runner 37, and then the input shaft 36. It is discharged through an oil passage inside. At this time, since the pressure difference between the chambers 39 and 40 on both sides of the lock-up clutch plate 38 hardly occurs, the lock-up clutch plate 38 is
4 and, at this time, the rotational force of the crankshaft 13 is transmitted through the pump cover 34, the pump impeller 35 and the turbine runner 37 to the input shaft 3.
6 is transmitted. At this time, the difference (N-NT) between the slip amount NSLP of the lock-up clutch, that is, the rotational speed N of the crankshaft, and the rotational speed NT of the input shaft 36 of the automatic converter 30, that is, the turbine runner 37, becomes considerably large. I have.
【0012】一方、ロックアップクラッチのスリップ制
御を行うべきときには入力軸36内のオイル通路を介し
て部屋40内に加圧オイルが供給され、部屋39内のオ
イルは入力軸36内のオイル通路を介して排出される。
このとき部屋40内の圧力が部屋39内の圧力よりも高
くなり、斯くしてロックアップクラッチ板38がポンプ
カバー34の内周面上に圧接されてクランクシャフトと
入力軸36とがスリップしつつ回転する状態となる。こ
のときポンプカバー34に対するロックアップクラッチ
板38の圧接力が強くなればなるほどスリップ量NSL
Pは小さくなり、従ってスリップ量NSLPは部屋3
9,40内に供給されるオイルを制御することによって
制御できる。なお、ロックアップクラッチのスリップ制
御時にはスリップ量NSLPは約50〜150r.p.m に
制御される。On the other hand, when the slip control of the lock-up clutch is to be performed, pressurized oil is supplied into the room 40 via the oil passage in the input shaft 36, and the oil in the room 39 passes through the oil passage in the input shaft 36. Is discharged through.
At this time, the pressure in the room 40 becomes higher than the pressure in the room 39, and thus the lock-up clutch plate 38 is pressed against the inner peripheral surface of the pump cover 34, and the crankshaft and the input shaft 36 slip. It will be in a rotating state. At this time, as the pressure contact force of the lock-up clutch plate 38 against the pump cover 34 increases, the slip amount NSL increases.
P becomes smaller, so that the slip amount NSLP becomes room 3
It can be controlled by controlling the oil supplied into 9,40. During the slip control of the lock-up clutch, the slip amount NSLP is controlled to about 50 to 150 rpm.
【0013】ロックアップクラッチのスリップ制御時に
比べて部屋40内に加圧オイルがさらに供給され、部屋
39からオイルがさらに排出されるとロックアップクラ
ッチ板38がポンプカバー34の内周面上に強く圧接さ
れてクランクシャフトと入力軸36とが同速度で回転す
る直結状態となり、斯くしてロックアップがオンとな
る。このときのスリップ量NSLPは零である。部屋3
9,40内へのオイルの供給制御、即ちロックアップ機
構33のオン・オフ制御およびロックアップクラッチの
スリップ制御は自動変速機30内に設けられた制御弁に
よって制御され、この制御弁は電子制御ユニット20の
出力信号に基づいて制御される。また、自動変速機20
内には変速作用を行うための多数のクラッチが設けられ
ており、これらのクラッチも電子制御ユニット20の出
力信号に基づいて制御される。When the pressurized oil is further supplied into the room 40 and the oil is further discharged from the room 39 as compared with the slip control of the lock-up clutch, the lock-up clutch plate 38 is strongly placed on the inner peripheral surface of the pump cover 34. When pressed, the crankshaft and the input shaft 36 rotate directly at the same speed, so that the lockup is turned on. The slip amount NSLP at this time is zero. Room 3
The control of oil supply to the insides 9, 40, that is, the on / off control of the lock-up mechanism 33 and the slip control of the lock-up clutch are controlled by a control valve provided in the automatic transmission 30, and this control valve is controlled electronically. It is controlled based on the output signal of the unit 20. Also, the automatic transmission 20
A number of clutches for performing a shifting operation are provided therein, and these clutches are also controlled based on an output signal of the electronic control unit 20.
【0014】また、自動変速機30内には入力軸36、
即ちタービンランナ37の回転速度を表わす出力パルス
を発生する回転速度センサ41および出力軸31の回転
速度を表わす出力パルスを発生する回転速度センサ42
が配置されており、これら回転速度センサ41,42の
出力パルスは入力ポート26に入力される。第1実施例
では燃料噴射時間TAUが次式に基づいて算出される。The automatic transmission 30 has an input shaft 36,
That is, a rotation speed sensor 41 that generates an output pulse indicating the rotation speed of the turbine runner 37 and a rotation speed sensor 42 that generates an output pulse indicating the rotation speed of the output shaft 31
The output pulses of the rotation speed sensors 41 and 42 are input to the input port 26. In the first embodiment, the fuel injection time TAU is calculated based on the following equation.
【0015】TAU=TP・FLEAN・FLLFB・
FAF+TAUV ここでTPは基本燃料噴射時間を、FLEANはリーン
補正係数を、FLLFBはリーンリミットフィードバッ
ク補正係数を、FAFは理論空燃比フィードバック補正
係数を、TAUVは無効噴射時間を夫々示している。基
本燃料噴射時間TPは空燃比を理論空燃比とするのに必
要な噴射時間を示している。この基本燃料噴射時間TP
は実験により求められ、この基本燃料噴射時間TPは機
関負荷Q/N(吸入空気量Q/機関回転数N)および機
関回転数Nの関数として図2に示すマップの形で予めR
OM22内に記憶されている。TAU = TP / FLEAN / FLLFFB /
FAF + TAUV Here, TP indicates a basic fuel injection time, FLEAN indicates a lean correction coefficient, FLLFB indicates a lean limit feedback correction coefficient, FAF indicates a stoichiometric air-fuel ratio feedback correction coefficient, and TAUV indicates an invalid injection time. The basic fuel injection time TP indicates an injection time required to make the air-fuel ratio a stoichiometric air-fuel ratio. This basic fuel injection time TP
Is determined experimentally. The basic fuel injection time TP is calculated as a function of the engine load Q / N (intake air amount Q / engine speed N) and the engine speed N in advance in the form of a map shown in FIG.
It is stored in the OM 22.
【0016】リーン補正係数FLEANは空燃比をリー
ン空燃比とするための補正係数であり、このリーン補正
係数FLEANは機関負荷Q/Nおよび機関回転数Nの
関数として図4に示すマップの形で予めROM22内に
記憶されている。リーンリミットフィードバック補正係
数FLLFBは空燃比をリーン限界に維持するための補
正係数である。第1実施例では吸入空気量Qと機関回転
数Nに対してリーン空燃比フィードバック制御に対する
学習領域が図5に示されるように例えば9つの領域で分
けられており、各学習領域に対して夫々リーンリミット
フィードバック補正係数FLLFB11〜FLLFB33が
設定されている。The lean correction coefficient FLEAN is a correction coefficient for making the air-fuel ratio a lean air-fuel ratio. The lean correction coefficient FLEAN is a function of the engine load Q / N and the engine speed N in the form of a map shown in FIG. It is stored in the ROM 22 in advance. The lean limit feedback correction coefficient FLLFB is a correction coefficient for maintaining the air-fuel ratio at the lean limit. In the first embodiment, the learning region for the lean air-fuel ratio feedback control is divided into, for example, nine regions with respect to the intake air amount Q and the engine speed N as shown in FIG. The lean limit feedback correction coefficients FLLFB 11 to FLLFB 33 are set.
【0017】理論空燃比フィードバック補正係数FAF
は空燃比を理論空燃比に維持するための係数である。理
論空燃比フィードバック補正係数FAFは空燃比を理論
空燃比に維持すべきときに空燃比センサ17の出力信号
に基づいて制御され、このとき理論空燃比フィードバッ
ク補正係数FAFはほぼ1.0を中心として上下動す
る。The stoichiometric air-fuel ratio feedback correction coefficient FAF
Is a coefficient for maintaining the air-fuel ratio at the stoichiometric air-fuel ratio. The stoichiometric air-fuel ratio feedback correction coefficient FAF is controlled based on the output signal of the air-fuel ratio sensor 17 when the air-fuel ratio is to be maintained at the stoichiometric air-fuel ratio. Move up and down.
【0018】図4に示されるように破線により囲まれた
運転領域内については機関の運転状態に応じてリーン補
正係数FLEANが定められており、この運転領域内で
は空燃比がリーン空燃比に維持される。これに対して図
4の破線で囲まれた領域外の運転領域では空燃比が理論
空燃比に維持される。空燃比を理論空燃比に維持すべき
ときにはリーン補正係数FLEANおよびリーンリミッ
トフィードバック補正係数FLLFBは1.0に固定さ
れ、理論空燃比フィードバック補正係数FAFが空燃比
センサ17の出力信号に基づいて制御される。As shown in FIG. 4, a lean correction coefficient FLEAN is determined in accordance with the operating state of the engine in an operating region surrounded by a broken line, and the air-fuel ratio is maintained at a lean air-fuel ratio in this operating region. Is done. On the other hand, the air-fuel ratio is maintained at the stoichiometric air-fuel ratio in the operation region outside the region surrounded by the broken line in FIG. When the air-fuel ratio should be maintained at the stoichiometric air-fuel ratio, the lean correction coefficient FLEAN and the lean limit feedback correction coefficient FLLFB are fixed to 1.0, and the stoichiometric air-fuel ratio feedback correction coefficient FAF is controlled based on the output signal of the air-fuel ratio sensor 17. You.
【0019】一方、空燃比をリーン空燃比に維持すべき
ときには理論空燃比フィードバック補正係数FAFが
1.0に固定され、即ち空燃比センサ17の出力信号に
基づくフィードバック制御が停止され、リーン補正係数
FLEANとリーンリミットフィードバック補正係数F
LLFBとにより空燃比がリーン空燃比に制御される。
次に図3を参照しつつリーンリミットフィードバック制
御について説明する。図3は機関出力トルク変動量およ
びNOx 発生量と空燃比との関係を示している。空燃比
がリーンになるほど燃料消費率は小さくなり、また空燃
比がリーンになるほどNOx の発生量が少なくなる。従
ってこれらの点からみると空燃比はできるだけリーンに
することが好ましいことになる。ところが空燃比が或る
程度以上リーンになると燃焼が不安定となり、その結果
図3に示されるようにトルク変動量が大きくなる。そこ
で第1実施例では図3に示されるようにトルク変動が増
大し始める空燃比制御領域内に空燃比を維持するように
している。On the other hand, when the air-fuel ratio is to be maintained at the lean air-fuel ratio, the stoichiometric air-fuel ratio feedback correction coefficient FAF is fixed at 1.0, that is, the feedback control based on the output signal of the air-fuel ratio sensor 17 is stopped, and the lean correction coefficient FLEAN and lean limit feedback correction coefficient F
The air-fuel ratio is controlled to a lean air-fuel ratio by the LLFB.
Next, the lean limit feedback control will be described with reference to FIG. Figure 3 shows the relationship between the engine output torque variation amount and NO x generation amount and the air-fuel ratio. Fuel consumption rate as the air-fuel ratio is lean becomes small, and the generation amount of the more NO x air-fuel ratio becomes lean is reduced. Therefore, from these points, it is preferable that the air-fuel ratio be as lean as possible. However, when the air-fuel ratio becomes lean to a certain degree or more, combustion becomes unstable, and as a result, the amount of torque fluctuation increases as shown in FIG. Therefore, in the first embodiment, as shown in FIG. 3, the air-fuel ratio is maintained within the air-fuel ratio control region where the torque fluctuation starts to increase.
【0020】即ち具体的に云うとリーン補正係数FLE
ANはリーンリミットフィードバック補正係数FLLF
BをFLLFB=1.0としたときに空燃比が図3に示
される空燃比制御領域の中央部となるように定められて
いる。一方、リーンリミットフィードバック補正係数F
LLFBはトルク変動量に応じて図3に示されるトルク
変動制御領域内において制御され、トルク変動量が大き
くなればリーンリミットフィードバック補正係数FLL
FBが増大せしめられ、即ち空燃比が小さくされ、トル
ク変動量が小さくなればリーンリミットフィードバック
補正係数FLLFBが減少せしめられ、即ち空燃比が大
きくされる。このようにして空燃比が図3に示される空
燃比制御領域内に制御される。That is, specifically, the lean correction coefficient FLE
AN is the lean limit feedback correction coefficient FLLF
When B is set to FLLFFB = 1.0, the air-fuel ratio is determined to be at the center of the air-fuel ratio control region shown in FIG. On the other hand, the lean limit feedback correction coefficient F
LLFB is controlled in the torque fluctuation control region shown in FIG. 3 in accordance with the torque fluctuation amount, and when the torque fluctuation amount becomes large, the lean limit feedback correction coefficient FLL
When FB is increased, that is, the air-fuel ratio is reduced, and the torque fluctuation amount is reduced, the lean limit feedback correction coefficient FLLFB is reduced, that is, the air-fuel ratio is increased. Thus, the air-fuel ratio is controlled within the air-fuel ratio control region shown in FIG.
【0021】なお、リーンリミットフィードバック補正
係数FLLFBはリーン補正係数FLEANが定められ
ている機関運転領域をカバーするように設定されてい
る。また、第1実施例では図3のトルク変動制御領域内
に目標トルク変動量を定め、トルク変動量が目標トルク
変動量となるように空燃比を制御するようにしている。Note that the lean limit feedback correction coefficient FLLFB is set so as to cover the engine operating region in which the lean correction coefficient FLEAN is determined. Further, in the first embodiment, the target torque fluctuation amount is determined in the torque fluctuation control region of FIG. 3, and the air-fuel ratio is controlled so that the torque fluctuation amount becomes the target torque fluctuation amount.
【0022】ところが冒頭で述べたようにトルク変動が
車両振動に与える影響はロックアップクラッチのスリッ
プ量NSLPに応じて変動し、従ってロックアップクラ
ッチのスリップ量NSLPに応じて図3のトルク変動制
御領域を定める必要があり、云い換えるとスリップ量N
SLPに応じて目標トルク変動量を定める必要がある。
即ち、同一のトルク変動に対しスリップ量NSLPが小
さいときほど車両振動が大きくなる。そこで、第1実施
例ではスリップ量NSLPが小さいときほど目標トルク
変動量が小さくなるようにしている。However, as described at the beginning, the influence of the torque fluctuation on the vehicle vibration fluctuates according to the slip amount NSLP of the lock-up clutch. In other words, the slip amount N
It is necessary to determine the target torque fluctuation amount according to the SLP.
That is, the vehicle vibration increases as the slip amount NSLP decreases with respect to the same torque fluctuation. Therefore, in the first embodiment, the smaller the slip amount NSLP is, the smaller the target torque fluctuation amount is.
【0023】第1実施例では目標トルク変動量LVLL
FBが次式に基づいて算出される。 LVLLFB=LVLLFBB−DLTF ここでLVLLFBBは基本目標トルク変動量を、DL
TFはスリップ量補正係数を夫々示している。基本目標
トルク変動量LVLLFBBはロックアップがオフのと
きに最適な予め実験により求められた目標トルク変動量
であり、等しい変動値を実線で示した図6(A)に示さ
れるように機関負荷Q/Nが高くなるほど大きくなり、
機関回転数Nが高くなるほど大きくなる。この基本目標
トルク変動値LVLLFBBは図6(B)に示されるよ
うに機関負荷Q/Nおよび機関回転数Nの関数としてマ
ップの形で予めROM22内に記憶されている。In the first embodiment, the target torque fluctuation amount LVLL
FB is calculated based on the following equation. LVLLFB = LVLLFB-DLTF Here, LVLLFB is the basic target torque variation, DL
TF indicates a slip amount correction coefficient. The basic target torque variation LVLLFBB is a target torque variation that is optimally obtained by an experiment in advance when the lockup is off, and the engine load Q as shown in FIG. The higher the / N, the larger the
It increases as the engine speed N increases. The basic target torque fluctuation value LVLLFBB is stored in the ROM 22 in advance in the form of a map as a function of the engine load Q / N and the engine speed N as shown in FIG.
【0024】スリップ量補正係数DLTFはスリップ量
NSLPに応じて基本目標トルク変動量LVLLFBB
を補正するものである。このスリップ量補正係数DLT
Fはロックアップがオフのときには零に維持され、ロッ
クアップがオンまたはロックアップクラッチのスリップ
制御時には図7に示されるマップから算出される。即
ち、スリップ量補正係数DLTFは図7に示されるよう
にスリップ量NSLPが小さくなるほど大きくなる。な
お、図7においてNSLP=0はロックアップがオンの
場合を示している。スリップ量補正係数DLTFが大き
くなるにつれて目標トルク変動量LVLLFBが小さく
なるのでスリップ量NSLPが小さくなるにつれて目標
トルク変動量LVLLFBが小さくされることになる。
このようにするとスリップ量NSLPに関わらず、即ち
スリップ制御中であるか否かに関わらず目標トルク変動
量LVLLFBを最適にすることができる。なお、スリ
ップ量補正係数DLTFはスリップ量NSLPの関数と
して図7に示すマップの形で予めROM22内に記憶さ
れている。The slip amount correction coefficient DLTF is calculated based on the slip amount NSLP based on the basic target torque fluctuation amount LVLLFBB.
Is to be corrected. This slip amount correction coefficient DLT
F is maintained at zero when the lockup is off, and is calculated from the map shown in FIG. 7 when the lockup is on or the slip control of the lockup clutch is performed. That is, the slip amount correction coefficient DLTF increases as the slip amount NSLP decreases as shown in FIG. In FIG. 7, NSLP = 0 indicates a case where lock-up is on. As the slip amount correction coefficient DLTF increases, the target torque fluctuation amount LVLLFB decreases. Therefore, as the slip amount NSLP decreases, the target torque fluctuation amount LVLLFB decreases.
This makes it possible to optimize the target torque variation LVLLFB regardless of the slip amount NSLP, that is, regardless of whether the slip control is being performed. The slip amount correction coefficient DLTF is stored in the ROM 22 in advance in the form of a map shown in FIG. 7 as a function of the slip amount NSLP.
【0025】ところで、ロックアップがオフのときには
前述したようにスリップ量NSLPがかなり大きくなる
のでスリップ量NSLPがしきい値例えば150r.
p.m.よりも高くなったときにスリップ量補正係数D
LTFが零となるように図7に示すマップを定め、ロッ
クアップがオフのときにも図7のマップからスリップ量
補正係数DLTFを算出するようにしてもよい。しかし
ながら、過渡運転時や変速時にはロックアップがオフで
あってもスリップ量NSLPがしきい値よりも低くなる
場合があり、或いはスリップ制御中であってもスリップ
量NSLPがしきい値よりも高くなる場合がある。そこ
で第1実施例ではロックアップがオフのときにはスリッ
プ量NSLPに関わらずスリップ量補正係数DLTFを
零に維持し、ロックアップがオンのときまたはスリップ
制御時に図7のマップを用いてスリップ量NSLPに応
じスリップ量補正係数DLTFを算出するようにしてい
る。When the lock-up is off, the slip amount NSLP becomes considerably large as described above.
p. m. The slip amount correction coefficient D
The map shown in FIG. 7 may be determined so that the LTF becomes zero, and the slip amount correction coefficient DLTF may be calculated from the map of FIG. 7 even when lockup is off. However, the slip amount NSLP may be lower than the threshold value during the transient operation or the shift even if the lockup is off, or the slip amount NSLP may be higher than the threshold value even during the slip control. There are cases. Therefore, in the first embodiment, when the lockup is off, the slip amount correction coefficient DLTF is maintained at zero regardless of the slip amount NSLP. The slip amount correction coefficient DLTF is calculated accordingly.
【0026】トルク変動量が目標トルク変動量になるよ
うに制御されると、即ち図3に示されるトルク変動制御
領域内に制御されると良好な車両の運転性を確保しつつ
燃料消費率およびNOx の発生量を大巾に低減すること
ができる。ただし、このようにトルク変動量をトルク変
動制御領域内に制御するためにはトルク変動量を検出し
なければならない。When the torque fluctuation amount is controlled to be equal to the target torque fluctuation amount, that is, when the torque fluctuation amount is controlled within the torque fluctuation control region shown in FIG. the generation amount of the NO x can be reduced by a large margin. However, in order to control the amount of torque fluctuation within the torque fluctuation control region, the amount of torque fluctuation must be detected.
【0027】ところでトルク変動量を算出する方法は従
来より種々の方法が提案されている。代表的な例を挙げ
ると燃焼室内に燃焼圧センサを取付けてこの燃焼圧セン
サの出力信号に基づきトルク変動量を算出する方法や、
或いは例えば特公平7−33809号公報に開示されて
いるようにクランクシャフトの第1の角速度ωaの2乗
とクランクシャフトの第2の角速度ωbの2乗との差に
基づいてトルク変動量を算出する方法が挙げられる。Various methods have conventionally been proposed for calculating the amount of torque fluctuation. A typical example is a method of mounting a combustion pressure sensor in a combustion chamber and calculating a torque fluctuation amount based on an output signal of the combustion pressure sensor,
Alternatively, as disclosed in, for example, Japanese Patent Publication No. 7-33809, the torque fluctuation amount is calculated based on the difference between the square of the first angular velocity ωa of the crankshaft and the square of the second angular velocity ωb of the crankshaft. Method.
【0028】即ち、各気筒において燃焼が行われると燃
焼圧によってクランクシャフトの角速度は第1の角速度
ωaから第2の角速度ωbへ上昇せしめられる。このと
き、機関の回転慣性モーメントをIとすると燃焼圧によ
って運動エネルギが(1/2)・Iωa2 から(1/
2)・Iωb2 へ上昇せしめられる。概略的に云うとこ
の運動エネルギの上昇量(1/2)・I・(ωb2 −ω
a2 )によってトルクが発生するので発生トルクは(ω
b2 −ωa2 )に比例することになる。従って発生トル
クは第1の角速度ωaの2乗と第2の角速度ωbの2乗
との差から求まることになり、従ってこのようにして求
めた発生トルクからトルク変動量を算出することができ
る。That is, when combustion is performed in each cylinder, the angular velocity of the crankshaft is increased from the first angular velocity ωa to the second angular velocity ωb by the combustion pressure. At this time, assuming that the rotational inertia moment of the engine is I, the kinetic energy is changed from (1/2) · Iωa 2 to (1 /
2) It is raised to Iωb 2 . Roughly speaking, the amount of increase in kinetic energy (1/2) · I · (ωb 2 −ω
a 2 ), the generated torque is (ω
b 2 −ωa 2 ). Accordingly, the generated torque is obtained from the difference between the square of the first angular velocity ωa and the square of the second angular velocity ωb. Therefore, the amount of torque fluctuation can be calculated from the generated torque thus obtained.
【0029】燃焼圧センサを用いると燃焼圧センサを取
付けた気筒が発生するトルクを確実に検出することがで
きるという利点がある反面、燃焼圧センサが必要である
という欠点を有している。これに対して角速度ωa,ω
bは従来より内燃機関が備えているクランク角センサの
出力信号から算出することができるので角速度ωa,ω
bに基づき出力トルクを算出するようにした場合には新
たなセンサを設ける必要がないという利点がある。ただ
し、この場合機関駆動系が捩り振動を生ずるとトルク変
動量を正確に検出できなくなるという問題を有してい
る。しかしながらこの問題を解決しさえすれば新たなセ
ンサを必要としない角速度に基づくトルク算出方法の方
が好ましいことは明らかである。そこで第1実施例では
発生トルクを角速度に基づき算出するようにし、その際
機関駆動系が捩り振動を生じたとしてもトルク変動量を
正確に検出しうるようにしている。The use of a combustion pressure sensor has the advantage that the torque generated by the cylinder to which the combustion pressure sensor is attached can be reliably detected, but has the disadvantage of requiring a combustion pressure sensor. On the other hand, the angular velocities ωa, ω
b can be calculated from the output signal of the crank angle sensor provided in the internal combustion engine, so that the angular velocities ωa, ω
When the output torque is calculated based on b, there is an advantage that it is not necessary to provide a new sensor. However, in this case, if the engine drive system generates torsional vibration, there is a problem that the torque fluctuation amount cannot be accurately detected. However, it is clear that a torque calculation method based on angular velocity that does not require a new sensor is preferable as long as this problem is solved. Therefore, in the first embodiment, the generated torque is calculated based on the angular velocity, and in this case, even if the engine drive system generates torsional vibration, the amount of torque fluctuation can be accurately detected.
【0030】次に機関の出力変動およびトルク変動を算
出する方法について説明する。まず初めに、機関駆動系
が捩り振動を生じていない定常運転時を示す図8
(A),(B)を参照しつつ各気筒が発生する駆動力お
よび各気筒が発生するトルクを算出する方法について説
明する。前述したようにクランク角センサ15はクラン
クシャフトが30°クランク角度回転する毎に出力パル
スを発生し、更にクランク角センサ15は各気筒#1,
#2,#3,#4の圧縮上死点TDCにおいて出力パル
スを発生するように配置されている。従ってクランク角
センサ15は各気筒#1,#2,#3,#4の圧縮上死
点TDCから30°クランク角毎に出力パルスを発生す
ることになる。なお、第1実施例において用いられてい
る内燃機関の点火順序は1−3−4−2である。Next, a method for calculating the output fluctuation and the torque fluctuation of the engine will be described. First, FIG. 8 shows a time of steady operation in which the engine drive system does not generate torsional vibration.
A method of calculating the driving force generated by each cylinder and the torque generated by each cylinder will be described with reference to (A) and (B). As described above, the crank angle sensor 15 generates an output pulse every time the crankshaft rotates by 30 ° crank angle.
The output pulses are generated at the compression top dead center TDC of # 2, # 3, and # 4. Therefore, the crank angle sensor 15 generates an output pulse every 30 ° crank angle from the compression top dead center TDC of each of the cylinders # 1, # 2, # 3, and # 4. The ignition sequence of the internal combustion engine used in the first embodiment is 1-3-4-2.
【0031】図8(A),(B)において縦軸T30は
クランク角センサ15が出力パルスを発生してから次の
出力パルスを発生するまでの30°クランク角度の経過
時間を表わしている。また、Ta(i)はi番気筒の圧
縮上死点(以下TDCと称す)から圧縮上死点後(以下
ATDCと称す)30°までの経過時間を示しており、
Tb(i)はi番気筒のATDC60°からATDC9
0°までの経過時間を示している。従って例えばTa
(1)は1番気筒のTDCからATDC30°までの経
過時間を示しており、Tb(1)は1番気筒のATDC
60°からATDC90°までの経過時間を示している
ことになる。一方、30°クランク角度を経過時間T3
0で除算するとこの除算結果は角速度ωを表わしてい
る。第1実施例では30°クランク角度/Ta(i)を
i番気筒における第1の角速度ωaと称し、30°クラ
ンク角度/Tb(i)をi番気筒における第2の角速度
ωbと称する。従って30°クランク角度/Ta(1)
は1番気筒の第1の角速度ωaを表わし、30°クラン
ク角度/Tb(1)は1番気筒の第2の角速度ωbを表
わすことになる。In FIGS. 8A and 8B, the vertical axis T30 represents the elapsed time of the 30 ° crank angle from the time when the crank angle sensor 15 generates an output pulse to the time when the next output pulse is generated. Further, Ta (i) indicates the elapsed time from the compression top dead center (hereinafter, referred to as TDC) of the i-th cylinder to 30 ° after the compression top dead center (hereinafter, referred to as ATDC),
Tb (i) is from ATDC60 ° of the i-th cylinder to ATDC9
The elapsed time up to 0 ° is shown. Therefore, for example, Ta
(1) indicates the elapsed time from TDC of the first cylinder to 30 ° ATDC, and Tb (1) indicates the ATDC of the first cylinder.
This indicates the elapsed time from 60 ° to ATDC 90 °. On the other hand, the 30 ° crank angle is changed to the elapsed time T3.
When divided by 0, the result of the division represents the angular velocity ω. In the first embodiment, 30 ° crank angle / Ta (i) is referred to as a first angular velocity ωa in the i-th cylinder, and 30 ° crank angle / Tb (i) is referred to as a second angular velocity ωb in the i-th cylinder. Therefore, 30 ° crank angle / Ta (1)
Represents the first angular velocity ωa of the first cylinder, and 30 ° crank angle / Tb (1) represents the second angular velocity ωb of the first cylinder.
【0032】図8(A),(B)の1番気筒に注目して
みると、燃焼が開始されて燃焼圧が高まると経過時間が
Ta(1)からTb(1)まで低下し、次いでTb
(1)から再び上昇する。云い換えるとクランクシャフ
トの角速度ωが第1の角速度ωaから第2の角速度ωb
まで上昇し、次いで第2の角速度ωbから再び下降す
る。即ち、燃焼圧によってクランクシャフトの角速度ω
が第1の角速度ωaから第2の角速度ωbへと増大せし
められたことになる。図8(A)は燃焼圧が比較的高い
場合を示しており、図8(B)は燃焼圧が比較的低い場
合を示している。図8(A),(B)から燃焼圧が高い
場合には燃焼圧が低い場合に比べて経過時間の減少量
(Ta(i)−Tb(i))が大きくなり、従って角速
度ωの増大量(ωb−ωa)が大きくなる。燃焼圧が高
くなればその気筒の発生する駆動力が大きくなり、従っ
て角速度ωの増大量(ωb−ωa)が大きくなれば気筒
の発生する駆動力が大きくなることになる。従って第1
の角速度ωaと第2の角速度ωbとの差(ωb−ωa)
から気筒の発生する駆動力を算出することができる。Looking at the first cylinder in FIGS. 8A and 8B, when the combustion starts and the combustion pressure increases, the elapsed time decreases from Ta (1) to Tb (1), and then. Tb
It rises again from (1). In other words, the angular velocity ω of the crankshaft is changed from the first angular velocity ωa to the second angular velocity ωb
And then fall again from the second angular velocity ωb. That is, the angular velocity ω of the crankshaft is determined by the combustion pressure.
Has been increased from the first angular velocity ωa to the second angular velocity ωb. FIG. 8A shows a case where the combustion pressure is relatively high, and FIG. 8B shows a case where the combustion pressure is relatively low. 8A and 8B, when the combustion pressure is high, the amount of decrease in the elapsed time (Ta (i) -Tb (i)) is larger than when the combustion pressure is low, and therefore the angular velocity ω is increased. The large amount (ωb−ωa) increases. When the combustion pressure increases, the driving force generated by the cylinder increases. Therefore, when the increase amount (ωb−ωa) of the angular velocity ω increases, the driving force generated by the cylinder increases. Therefore the first
(Ωb−ωa) between the second angular velocity ωb and the second angular velocity ωa
, The driving force generated by the cylinder can be calculated.
【0033】一方、機関の回転慣性モーメントをIとす
ると燃焼圧によって運動エネルギが(1/2)Iωa2
から(1/2)Iωb2 に増大せしめられる。この運動
エネルギの増大量(1/2)・I・(ωb2 −ωa2 )
はその気筒が発生するトルクを表わしており、従って第
1の角速度ωaの2乗と第2の角速度ωbの2乗との差
(ωb2 −ωa2 )から気筒の発生するトルクを算出で
きることになる。On the other hand, assuming that the rotational inertia moment of the engine is I, the kinetic energy is (1/2) Iωa 2 due to the combustion pressure.
From () Iωb 2 . Increased amount of kinetic energy (1/2) · I · (ωb 2 -ωa 2)
It is capable of calculating the square and the torque generated by the cylinder from the difference (ωb 2 -ωa 2) of the square of the second angular velocity [omega] b for represents the torque that cylinder will occur, thus the first angular velocity ωa Become.
【0034】このように第1の角速度ωaと第2の角速
度ωbを検出すればこれらの検出値から対応する気筒の
発生する駆動力および対応する気筒の発生するトルクを
算出できることになる。なお、図8(A),(B)に示
される経過時間T30の変化は機関によって若干異な
り、従って第1の角速度ωaを検出すべきクランク角度
範囲および第2の角速度ωbを検出すべきクランク角度
範囲は機関に応じて(ωb−ωa)が機関の発生する駆
動力を最もよく表わすように、或いは(ωb2 −ωa
2 )が機関の発生するトルクを最もよく表わすように定
められる。従って機関によっては第1の角速度ωaを検
出すべきクランク角度範囲が圧縮上死点前BTDC30
°からTDCであり、第2の角速度ωbを検出すべきク
ランク角度範囲がATDC90°からATDC120°
となることもあり得る。Thus, by detecting the first angular velocity ωa and the second angular velocity ωb, the driving force generated by the corresponding cylinder and the torque generated by the corresponding cylinder can be calculated from these detected values. The changes in the elapsed time T30 shown in FIGS. 8A and 8B slightly differ depending on the engine, and therefore, the crank angle range for detecting the first angular velocity ωa and the crank angle for detecting the second angular velocity ωb The range is set so that (ωb−ωa) best represents the driving force generated by the engine, or (ωb 2 −ωa) depending on the engine.
2 ) is determined to best represent the torque generated by the engine. Therefore, depending on the engine, the crank angle range in which the first angular velocity ωa should be detected is the BTDC 30 before the compression top dead center.
° to TDC, and the crank angle range for detecting the second angular velocity ωb is from ATDC 90 ° to ATDC 120 °.
It is possible that
【0035】従って各角速度ωa,ωbの検出のしかた
について一般的に表現すると、圧縮行程末期から爆発行
程初期までのクランク角度領域内に第1のクランク角度
範囲を設定し、第1のクランク角度範囲から一定のクラ
ンク角を隔てた爆発行程中期のクランク角度領域内に第
2のクランク角度範囲を設定し、第1のクランク角度範
囲内におけるクランクシャフトの第1の角速度ωaを検
出し、第2のクランク角度範囲内におけるクランクシャ
フトの第2の角速度ωbを検出するということになる。Therefore, in general terms, how to detect each of the angular velocities ωa and ωb is as follows. The first crank angle range is set within the crank angle range from the end of the compression stroke to the beginning of the explosion stroke. A second crank angle range is set in a crank angle region in a middle stage of an explosion stroke separated by a certain crank angle from the first crank angle range, a first angular velocity ωa of the crankshaft in the first crank angle range is detected, and a second This means that the second angular velocity ωb of the crankshaft within the crank angle range is detected.
【0036】上述したように角速度ωa,ωbを検出す
れば検出値に基づいて対応する気筒の発生する駆動力お
よびトルクを算出することができる。しかしながら機関
駆動系には各気筒において順次行われる爆発作用により
駆動系の固有振動数でもって振動する捩り振動が発生し
ており、このように機関駆動系に捩り振動が発生してい
ると角速度ωa,ωbに基づいて気筒の発生する駆動力
およびトルクを正確に算出することができなくなる。次
にこのことについて図9および図10を参照しつつ説明
する。As described above, if the angular velocities ωa and ωb are detected, the driving force and torque generated by the corresponding cylinder can be calculated based on the detected values. However, in the engine drive system, torsional vibration oscillating at the natural frequency of the drive system is generated by the explosion effect sequentially performed in each cylinder, and when the torsional vibration is generated in the engine drive system, the angular velocity ωa , Ωb, the driving force and torque generated by the cylinder cannot be accurately calculated. Next, this will be described with reference to FIGS.
【0037】図9は機関駆動系に捩り振動が発生してい
るときに各気筒に対し順次算出される経過時間Ta
(i)の変化を示している。機関駆動系に捩り振動が発
生するとこの捩り振動によってクランクシャフトの角速
度が周期的に増大減少せしめられるので経過時間Ta
(i)は図9に示されるように周期的に増大減少するこ
とになる。FIG. 9 shows an elapsed time Ta sequentially calculated for each cylinder when torsional vibration occurs in the engine drive system.
The change of (i) is shown. When torsional vibration occurs in the engine drive system, the torsional vibration causes the angular velocity of the crankshaft to periodically increase and decrease, so that the elapsed time Ta
(I) periodically increases and decreases as shown in FIG.
【0038】一方、図10は図9において経過時間Ta
(i)が減少している部分を拡大して示している。図1
0に示されるように経過時間Ta(i)はTa(1)と
Ta(3)との間でho時間だけ減少しており、このh
o時間の減少は捩り振動による捩れ量の増大によるもの
と考えられる。この場合、Ta(1)とTa(3)との
間では捩り振動による経過時間の減少量は時間の経過と
共にほぼ直線的に増大するものと考えられ、従ってこの
捩り振動による経過時間の減少量はTa(1)およびT
a(3)を結ぶ破線とTa(1)を通る水平線との差で
表わされることになる。従ってTa(1)とTb(1)
との間では捩り振動によって経過時間がhだけ減少して
いることになる。FIG. 10 shows the elapsed time Ta in FIG.
The portion where (i) is decreasing is shown in an enlarged manner. FIG.
As shown by 0, the elapsed time Ta (i) decreases by the ho time between Ta (1) and Ta (3), and this h
It is considered that the decrease in the o time is due to an increase in the amount of torsion due to torsional vibration. In this case, between Ta (1) and Ta (3), the amount of decrease in the elapsed time due to the torsional vibration is considered to increase almost linearly with the passage of time, and therefore, the amount of decrease in the elapsed time due to this torsional vibration. Are Ta (1) and T
It is represented by the difference between the dashed line connecting a (3) and the horizontal line passing through Ta (1). Therefore, Ta (1) and Tb (1)
The elapsed time is reduced by h due to the torsional vibration.
【0039】このように機関駆動系に捩り振動が発生す
るとTb(1)はTa(1)に対して経過時間が減少
し、この減少した経過時間は燃焼圧による経過時間の減
少量fと捩り振動による経過時間の減少量hとを含んで
いることになる。従って燃焼圧により減少した経過時間
Tb(1)だけを求めるためにはTb(1)にhを加算
しなければならないことになる。即ち、検出された経過
時間Ta(i)およびTb(i)に基づいて各気筒が発
生する駆動力或いはトルクを求めても真の駆動力或いは
トルクを求めることができず、斯くして真の機関の出力
変動或いはトルク変動を求めることができない。As described above, when torsional vibration occurs in the engine drive system, the elapsed time of Tb (1) decreases with respect to Ta (1). This includes the decrease amount h of the elapsed time due to the vibration. Therefore, in order to obtain only the elapsed time Tb (1) reduced by the combustion pressure, h must be added to Tb (1). That is, even if the driving force or torque generated by each cylinder is obtained based on the detected elapsed times Ta (i) and Tb (i), the true driving force or torque cannot be obtained, and thus the true driving force or torque cannot be obtained. The engine output fluctuation or torque fluctuation cannot be determined.
【0040】更に、多気筒内燃機関ではこのような機関
駆動系の捩り振動に加えてクランクシャフト自体の捩り
振動を発生し、このようなクランクシャフト自体の捩り
振動が発生した場合にも真の機関の出力変動或いはトル
ク変動を求めることができなくなる。次にこのことにつ
いて図11を参照しつつ説明する。多気筒内燃機関、例
えば図1に示されるような4気筒内燃機関では1番気筒
および2番気筒において大きなクランクシャフト自体の
捩り振動が発生する。即ち、クランクシャフト自体の捩
り振動がほとんど発生しない気筒、例えば3番気筒#3
では図11に示されるようにTa(3)からTb(3)
に向けて経過時間は徐々に減少するが1番気筒#1にお
いてはTa(1)からTb(1)に向けて経過時間は徐
々に減少せず、クランクシャフト自体の捩り振動によっ
てTb(1)の経過時間が長くなってしまう。その結
果、1番気筒#1については検出された経過時間Ta
(1)およびTb(1)に基づいて1番気筒#1が発生
する駆動力或いはトルクを求めても真の駆動力或いはト
ルクを求めることができず、斯くして真の機関の出力変
動或いはトルク変動を求めることができない。Further, in a multi-cylinder internal combustion engine, in addition to such torsional vibration of the engine drive system, a torsional vibration of the crankshaft itself is generated. Output fluctuation or torque fluctuation cannot be obtained. Next, this will be described with reference to FIG. In a multi-cylinder internal combustion engine, for example, a four-cylinder internal combustion engine as shown in FIG. 1, a large torsional vibration of the crankshaft itself occurs in the first and second cylinders. That is, a cylinder in which torsional vibration of the crankshaft itself hardly occurs, for example, the third cylinder # 3
Now, from Ta (3) to Tb (3) as shown in FIG.
, The elapsed time decreases gradually from Ta (1) to Tb (1) in the first cylinder # 1, and Tb (1) due to torsional vibration of the crankshaft itself. The elapsed time becomes longer. As a result, for the first cylinder # 1, the detected elapsed time Ta
Even if the driving force or torque generated by the first cylinder # 1 is determined based on (1) and Tb (1), the true driving force or torque cannot be determined, and thus the true engine output fluctuation or The torque fluctuation cannot be determined.
【0041】そこで第1実施例では機関駆動系の捩り振
動が発生しても、またクランクシャフト自体の捩り振動
が発生しても真の機関の出力変動或いはトルク変動を算
出することができる方法を採用している。次にこの算出
方法について図12を参照しつつ説明する。図12にお
いてTa(1)j-1 およびTb(1)j-1 は1番気筒#
1の先の燃焼時における経過時間を表しており、Ta
(1)j およびTb(1)j は1番気筒#1の次の燃焼
時における経過時間を表している。一方、Ta(3)
j-1 は1番気筒#1の先の燃焼にひき続いて行われる3
番気筒#3の先の燃焼時における経過時間を表してお
り、Ta(3)j は3番気筒#3の次の燃焼時における
経過時間を表わしている。Therefore, in the first embodiment, a method capable of calculating a true engine output fluctuation or torque fluctuation even if torsional vibration of the engine drive system or torsion vibration of the crankshaft itself occurs. Has adopted. Next, this calculation method will be described with reference to FIG. In FIG. 12, Ta (1) j-1 and Tb (1) j-1 are cylinders # 1
1 represents the elapsed time during the preceding combustion, and Ta
(1) j and Tb (1) j represents the elapsed time during the next combustion of the first cylinder # 1. On the other hand, Ta (3)
j-1 is performed following the combustion of the first cylinder # 1 3
The elapsed time during the previous combustion of the # 3 cylinder # 3 is shown, and Ta (3) j represents the elapsed time during the next combustion of the # 3 cylinder # 3.
【0042】まず初めにTa(1)j とTa(1)j-1
との差DTa(1)(=Ta(1) j −Ta(1)
j-1 )、およびTa(3)j とTb(3)j-1 との差D
Ta(3)(=Ta(3)j −Ta(3)j-1 )とを求
める。次いで1番気筒#1の先の燃焼時の出力トルクと
次の燃焼時の出力トルクとが同一であったと仮定したと
きの1番気筒#1の次の燃焼時におけるATDC60°
からATDC90°までの仮想の経過時間をTb′
(1)j とし、このTb′(1)j とTb(1)j-1 と
の差Kb(1)(=Tb′(1)j −Tb(1)j-1 )
を求める。これらの差DTa(1),Kb(1)および
DTa(3)を一直線上における高さとして書き直すと
図13に示されるようになる。First, Ta (1)jAnd Ta (1)j-1
DTa (1) (= Ta (1) j-Ta (1)
j-1 ), And Ta (3)jAnd Tb (3)j-1 And the difference D
Ta (3) (= Ta (3)j-Ta (3)j-1 ) And
Confuse. Next, the output torque of the first cylinder # 1 during the preceding combustion and
Assuming that the output torque during the next combustion is the same
ATDC at the time of the next combustion of the first cylinder # 1
Is the virtual elapsed time from Tb 'to ATDC 90 °
(1)jAnd Tb '(1)jAnd Tb (1)j-1 When
Difference Kb (1) (= Tb '(1)j-Tb (1)j-1 )
Ask for. These differences DTa (1), Kb (1) and
When DTa (3) is rewritten as a height on a straight line,
As shown in FIG.
【0043】図13に示されるように経過時間差DTa
(1)とDTa(3)との間では経過時間差がI0 時間
だけ減少している。クランクシャフト自体に捩り振動が
発生したとしてもこの捩り振動による影響は経過時間差
には表れず、従って経過時間差のI0 時間の減少は機関
駆動系の捩り振動によるものである。この場合、DTa
(1)とDTa(3)との間では機関駆動系の捩り振動
による経過時間差の減少量は時間の経過と共にほぼ直線
的に増大するものと考えられる。従って1番気筒#1の
先の燃焼時の出力トルクと次の燃焼時の出力トルクとが
同一であると仮定すると経過時間差DTa(1)とKb
(1)との間では機関駆動系の捩り振動によって経過時
間差がIだけ減少するものと考えられる。従って図13
からわかるように1番気筒#1の先の燃焼時の出力トル
クと次の燃焼時の出力トルクとが同一であると仮定した
場合に経過時間差Kb(1)は次式で表される。As shown in FIG. 13, the elapsed time difference DTa
The elapsed time difference between (1) and DTa (3) is reduced by I0 time. Even if torsional vibration occurs in the crankshaft itself, the effect of the torsional vibration does not appear in the elapsed time difference, and therefore, the decrease in the I0 time of the elapsed time difference is due to the torsional vibration of the engine drive system. In this case, DTa
It is considered that the amount of decrease in the elapsed time difference between (1) and DTa (3) due to torsional vibration of the engine drive system increases almost linearly with time. Therefore, assuming that the output torque of the first cylinder # 1 during the previous combustion is the same as the output torque during the next combustion, the elapsed time difference DTa (1) and Kb
It is considered that the elapsed time difference decreases by I due to the torsional vibration of the engine drive system between (1). Therefore, FIG.
As can be seen from the above, the elapsed time difference Kb (1) is expressed by the following equation, assuming that the output torque of the first cylinder # 1 during the previous combustion and the output torque during the next combustion are the same.
【0044】 Kb(1)=(2DTa(1)+DTa(3))/3 ところで1番気筒#1の先の燃焼時の出力トルクと次の
燃焼時の出力トルクとが同一であると仮定したときの1
番気筒#1の仮想の経過時間Tb′(1)j は次式で表
される。 Tb′(1)j =Tb(1)j-1 +Kb(1) 従って1番気筒#1の先の燃焼時の出力トルクと次の燃
焼時における出力トルクが同一であると仮定したときの
1番気筒#1の次の燃焼時における仮想の出力トルクD
NS(1)は次式で表される。 DNS(1)={30°/Tb′(1)j }2 −{30°/Ta(1)j }2 ={30°/(Tb(1)j-1 +Kb(1))}2 −{30°/Ta(1)j }2 一方、1番気筒#1の次の燃焼時におけるATDC60
°からATDC90°までの実際の経過時間Tb(1)
j を用いて実際の出力トルクを算出するとこの出力トル
クDN(1)は次式で表される。Kb (1) = (2DTa (1) + DTa (3)) / 3 By the way, it is assumed that the output torque of the first cylinder # 1 during the previous combustion is the same as the output torque during the next combustion. Time 1
The virtual elapsed time Tb ′ (1) j of the cylinder # 1 is expressed by the following equation. Tb '(1) j = Tb (1) j-1 + Kb (1) Accordingly, when the output torque of the first cylinder # 1 during the previous combustion is assumed to be the same as the output torque during the next combustion, 1 Virtual output torque D during the next combustion of cylinder # 1
NS (1) is represented by the following equation. DNS (1) = {30 ° / Tb ′ (1) j } 2 − {30 ° / Ta (1) j } 2 = {30 ° / (Tb (1) j−1 + Kb (1))} 2 − {30 ° / Ta (1) j } 2 ATDC60 at the time of the next combustion of the first cylinder # 1
Actual elapsed time Tb (°) to ATDC 90 ° (1)
When the actual output torque is calculated using j , this output torque DN (1) is expressed by the following equation.
【0045】DN(1)={30°/Tb(1)j }2
−{30°/Ta(1)j }2 この場合、1番気筒#1の仮想の出力トルクDNS
(1)と実際の出力トルクDN(1)との差は1番気筒
#1のトルク変動量を表しており、従って1番気筒#1
のトルク変動量DLN(1)は次式で表されることにな
る。 DLN(1)=DNS(1)−DN(1)={30°/
Tb′(1)j }2−{30°/Tb(1)j }2 一般的に表すとi番気筒のトルク変動量DLN(i)は
次式で表すことができる。DN (1) = {30 ° / Tb (1) j } 2
− {30 ° / Ta (1) j } 2 In this case, the virtual output torque DNS of the first cylinder # 1
The difference between (1) and the actual output torque DN (1) indicates the amount of torque fluctuation of the first cylinder # 1, and therefore, the first cylinder # 1
The torque fluctuation amount DLN (1) is expressed by the following equation. DLN (1) = DNS (1) −DN (1) = {30 ° /
Tb ′ (1) j } 2 − {30 ° / Tb (1) j } 2 Generally, the torque fluctuation amount DLN (i) of the i-th cylinder can be expressed by the following equation.
【0046】DLN(i)=DNS(i)−DN(i)
={30°/Tb′(i)}2−{30°/Tb
(i)}2 即ち、i番気筒のトルク変動量はi番気筒の仮想の第2
の角速度の2乗と実際の第2の角速度の2乗との偏差か
ら求まることになる。このような方法でもってトルク変
動量を算出することによって機関駆動系に捩り振動が発
生したとしても、更にクランクシャフト自体に捩り振動
が発生したとしても各気筒のトルク変動量を正確に検出
できることになる。一方、このような方法でもってトル
ク変動量を算出する場合、ロータ14(図1)の外周に
沿って形成されている外歯の間隔にばらつきがあると各
気筒のトルク変動量を正確に検出することができない。
そこで第1実施例ではロータ14の外歯の間隔にばらつ
きがあったとしても各気筒のトルク変動量を正確に検出
しうるように検出された各経過時間Ta(i),Tb
(i)を補正するようにしている。次に経過時間Ta
(i)の補正方法について図14を参照しつつ説明す
る。DLN (i) = DNS (i) -DN (i)
= {30 ° / Tb '(i)} 2- {30 ° / Tb
(I)} 2 That is, the torque fluctuation amount of the i-th cylinder is the virtual second
And the square of the actual second angular velocity. By calculating the amount of torque fluctuation by such a method, it is possible to accurately detect the amount of torque fluctuation of each cylinder even if torsional vibration occurs in the engine drive system or even if torsional vibration occurs in the crankshaft itself. Become. On the other hand, when calculating the torque fluctuation amount by such a method, if there is a variation in the interval between the external teeth formed along the outer periphery of the rotor 14 (FIG. 1), the torque fluctuation amount of each cylinder is accurately detected. Can not do it.
Therefore, in the first embodiment, each elapsed time Ta (i), Tb is detected so that the amount of torque fluctuation of each cylinder can be accurately detected even if the interval between the external teeth of the rotor 14 varies.
(I) is corrected. Next, the elapsed time Ta
The correction method (i) will be described with reference to FIG.
【0047】第1実施例では経過時間Ta(i)が次式
に基づいて算出される。 Ta(i)=(i番気筒のTDCからATDC30°ま
での所要時間)・(1+KTa(i)) ここでKTa(i)はi番気筒に対する補正係数を表し
ており、これらの補正係数KTa(i)は機関駆動系の
捩り振動が発生せず、クランクシャフト自体の捩り振動
も発生しない減速運転時の燃料供給停止時に算出され
る。In the first embodiment, the elapsed time Ta (i) is calculated based on the following equation. Ta (i) = (time required from TDC of cylinder i to 30 ° ATDC). (1 + KTa (i)) where KTa (i) represents a correction coefficient for the i-th cylinder, and these correction coefficients KTa ( i) is calculated when fuel supply is stopped during deceleration operation in which no torsional vibration of the engine drive system occurs and no torsional vibration of the crankshaft itself occurs.
【0048】図14は減速運転時における経過時間Ta
(i)の変化を示している。減速運転時にはクランクシ
ャフトの回転速度が次第に低下するので経過時間Ta
(i)はTa(1)j-1 ,Ta(3)j-1 ,Ta(4)
j-1 ,Ta(2)j-1 ,Ta(1)j ,Ta(3)j ,
Ta(4)j ,Ta(2)j で表されるように次第に増
大する。このときロータ14の外歯の間隔が同一であっ
たとすると各経過時間Ta(i)は図14の実線に沿っ
て変化するものと考えられ、検出された経過時間Ta
(i)が図14の実線からずれている場合には図14の
実線からずれている経過時間Ta(i)を検出するため
の外歯の間隔が正規の間隔からずれているものと考えら
れる。そこで第1実施例では経過時間Ta(i)が図1
4に示す直線からずれている場合にはその経過時間Ta
(i)が図14の実線上に位置するように経過時間Ta
(i)を補正係数KTa(i)によって補正するように
している。FIG. 14 shows the elapsed time Ta during the deceleration operation.
The change of (i) is shown. Since the rotation speed of the crankshaft gradually decreases during the deceleration operation, the elapsed time Ta
(I) is Ta (1) j−1 , Ta (3) j−1 , Ta (4)
j-1 , Ta (2) j-1 , Ta (1) j , Ta (3) j ,
It gradually increases as represented by Ta (4) j and Ta (2) j . At this time, if the intervals between the external teeth of the rotor 14 are the same, each elapsed time Ta (i) is considered to change along the solid line in FIG. 14, and the detected elapsed time Ta
When (i) deviates from the solid line in FIG. 14, it is considered that the interval between the external teeth for detecting the elapsed time Ta (i) deviating from the solid line in FIG. 14 deviates from the regular interval. . Therefore, in the first embodiment, the elapsed time Ta (i) is different from that in FIG.
In the case of deviation from the straight line shown in FIG.
The elapsed time Ta is set so that (i) is located on the solid line in FIG.
(I) is corrected by the correction coefficient KTa (i).
【0049】具体的に云うと、第1実施例では720°
クランク角範囲(j−1),(j)…における4気筒の
経過時間Ta(i)の平均値TaAVj-1 {=(Ta
(1) j-1 +Ta(3)j-1 +Ta(4)j-1 +Ta
(2)j-1 )/4},TaAVj{=(Ta(1)j +
Ta(3)j +Ta(4)j +Ta(2)j )/4},
…を順次算出する。次いで各経過時間Ta(4)j-1 ,
Ta(2)j-1 ,Ta(1)j ,Ta(3)j が平均値
TaAVj-1 から平均値TaAVj まで平均的な一定の
増大率でもって増大したと考え、平均的な一定の増大率
でもって増大したと考えたときの各気筒についての経過
時間AAV(4),AAV(2),AAV(1),AA
V(3)を次式より求める。More specifically, in the first embodiment, 720 ° is used.
The four cylinders in the crank angle range (j-1), (j) ...
Average value TaAV of elapsed time Ta (i)j-1 {= (Ta
(1) j-1 + Ta (3)j-1 + Ta (4)j-1 + Ta
(2)j-1 ) / 4}, TaAVj{= (Ta (1)j+
Ta (3)j+ Ta (4)j+ Ta (2)j) / 4},
Are sequentially calculated. Next, each elapsed time Ta (4)j-1 ,
Ta (2)j-1 , Ta (1)j, Ta (3)jIs the average
TaAVj-1 From the average value TaAVjAverage constant up to
Considering that it increased with the growth rate, the average constant growth rate
The progress of each cylinder when it is considered that it has increased
Time AAV (4), AAV (2), AAV (1), AA
V (3) is obtained from the following equation.
【0050】AAV(4)=(TaAVj −TaAV
j-1 )・(1/8)+TaAVj-1 AAV(2)=(TaAVj −TaAVj-1 )・(3/
8)+TaAVj-1 AAV(1)=(TaAVj −TaAVj-1 )・(5/
8)+TaAVj-1 AAV(3)=(TaAVj −TaAVj-1 )・(7/
8)+TaAVj-1 次いで次式に基づき各気筒についてこれら経過時間AA
V(4),AAV(2),AAV(1),AAV(3)
に対する実際に検出された経過時間Ta(4) j-1 ,T
a(2)j-1 ,Ta(1)j ,Ta(3)j のずれ割合
KAFC(4),KAFC(2),KAFC(1),K
AFC(3)を求める。AAV (4) = (TaAVj-TaAV
j-1 ) ・ (1/8) + TaAVj-1 AAV (2) = (TaAVj-TaAVj-1 ) ・ (3 /
8) + TaAVj-1 AAV (1) = (TaAVj-TaAVj-1 ) ・ (5 /
8) + TaAVj-1 AAV (3) = (TaAVj-TaAVj-1 ) ・ (7 /
8) + TaAVj-1 Next, the elapsed time AA for each cylinder is calculated based on the following equation.
V (4), AAV (2), AAV (1), AAV (3)
Elapsed time Ta (4) actually detected with respect to j-1 , T
a (2)j-1 , Ta (1)j, Ta (3)jDeviation ratio
KAFC (4), KAFC (2), KAFC (1), K
Find AFC (3).
【0051】KAFC(4)={Ta(4)j-1 −AA
V(4)}/AAV(4) KAFC(2)={Ta(2)j-1 −AAV(2)}/
AAV(2) KAFC(1)={Ta(1)j −AAV(1)}/A
AV(1) KAFC(3)={Ta(3)j −AAV(3)}/A
AV(3) 次に次式を用いて基本補正係数KTaB(i)を各気筒
毎に求める。KAFC (4) = {Ta (4) j−1 −AA
V (4)} / AAV (4) KAFC (2) = {Ta (2) j−1 −AAV (2)} /
AAV (2) KAFC (1) = {Ta (1) j −AAV (1)} / A
AV (1) KAFC (3) = {Ta (3) j −AAV (3)} / A
AV (3) Next, a basic correction coefficient KTaB (i) is obtained for each cylinder using the following equation.
【0052】KTaB(i)=KTa(i)j-1 −KA
FC(i)/4 即ち、各気筒に対するこれまでの補正係数KTa(i)
j-1 からずれ割合KAFC(i)の1/4を減算するこ
とによって基本補正係数KTaB(i)が算出される。
次いで基本補正係数KTaB(i)の平均値tKTaM
が次式に基づき算出される。KTaB (i) = KTa (i) j-1 -KA
FC (i) / 4 That is, the previous correction coefficient KTa (i) for each cylinder
The basic correction coefficient KTaB (i) is calculated by subtracting 1/4 of the deviation ratio KAFC (i) from j-1 .
Next, the average value tKTaM of the basic correction coefficient KTaB (i)
Is calculated based on the following equation.
【0053】tKTaM={KTaB(1)+KTaB
(2)+KTaB(3)+KTaB(4)}/4 次いで次式に示されるように各基本補正係数KTaB
(i)から平均値tKTaMを減算することによって各
気筒に対する補正係数KTa(i)が算出される。TKTaM = {KTaB (1) + KTaB
(2) + KTaB (3) + KTaB (4)} / 4 Then, as shown in the following equation, each basic correction coefficient KTaB
The correction coefficient KTa (i) for each cylinder is calculated by subtracting the average value tKTaM from (i).
【0054】 KTa(i)=KTaB(i)−tKTaM 上述したように補正係数KTa(i)j-1 をずれ割合K
AFC(i)に基づいて直接補正せず、一旦基本補正係
数KTaB(i)を求めるようにしているのは次の理由
による。例えば1番気筒#1についてのみずれ割合KA
FC(1)が存在し、このずれ割合KAFC(1)に基
づいて1番気筒#1の補正係数KTa(1)のみを補正
したとする。この場合、1番気筒#1のみの経過時間T
a(1)が増大又は減少せしめられる。しかしながら1
番気筒#1のみの経過時間Ta(1)が増大又は減少せ
しめられると今度は残りの気筒#2,#3,#4の補正
係数KTa(2),KTa(3),KTa(4)がずれ
を生じることになる。KTa (i) = KTaB (i) −tKTaM As described above, the correction coefficient KTa (i) j−1 is set to the deviation ratio K.
The reason why the basic correction coefficient KTaB (i) is once obtained without directly correcting based on AFC (i) is as follows. For example, the deviation ratio KA only for the first cylinder # 1
It is assumed that FC (1) exists and only the correction coefficient KTa (1) of the first cylinder # 1 is corrected based on the deviation ratio KAFC (1). In this case, the elapsed time T of only the first cylinder # 1
a (1) is increased or decreased. However 1
When the elapsed time Ta (1) of only the cylinder # 1 is increased or decreased, the correction coefficients KTa (2), KTa (3), and KTa (4) of the remaining cylinders # 2, # 3, and # 4 are increased. A shift will occur.
【0055】このような問題が生じないようにするため
に基本補正係数KTaB(i)を一旦求め、この基本補
正係数KTaB(i)から基本補正係数の平均値tKT
aMを減算することによって最終的な補正係数KTa
(i)を求めるようにしている。即ち、例えば1番気筒
#1の基本補正係数KTaB(1)が増大せしめられた
場合には全ての気筒#1,#2,#3,#4の基本補正
係数KTaB(i)から基本補正係数の平均値tKTa
Mが減算される。このように全ての気筒の基本補正係数
KTaB(i)からtKTaMが減算されると1番気筒
#1の補正係数KTa(1)は増大するが残りの気筒#
2,#3,#4の補正係数KTa(2),KTa
(3),KTa(4)は減少せしめられ、KTa
(1),KTa(2),KTa(3)およびKTa
(4)の和は常に零に維持される。このようにKTa
(1),KTa(2),KTa(3)およびKTa
(4)の和が常に零になるように各補正係数KTa
(i)を補正するといずれか一つの補正係数KTa
(i)が補正されても他の補正係数KTa(i)はずれ
を生じなくなる。In order to prevent such a problem from occurring, a basic correction coefficient KTaB (i) is once obtained, and an average value tKT of the basic correction coefficient is calculated from the basic correction coefficient KTaB (i).
The final correction coefficient KTa is obtained by subtracting aM.
(I) is obtained. That is, for example, when the basic correction coefficient KTaB (1) of the first cylinder # 1 is increased, the basic correction coefficient KTaB (i) of all the cylinders # 1, # 2, # 3, and # 4 is calculated. Average value of tKTa
M is subtracted. As described above, when tKTaM is subtracted from the basic correction coefficient KTaB (i) of all cylinders, the correction coefficient KTa (1) of the first cylinder # 1 increases, but the remaining cylinders ##
2, # 3, # 4 correction coefficients KTa (2), KTa
(3), KTa (4) is reduced and KTa
(1), KTa (2), KTa (3) and KTa
The sum of (4) is always maintained at zero. Thus, KTa
(1), KTa (2), KTa (3) and KTa
Each correction coefficient KTa is such that the sum of (4) is always zero.
When (i) is corrected, one of the correction coefficients KTa
Even if (i) is corrected, other correction coefficients KTa (i) do not shift.
【0056】図14に示される実施例では例えば1番気
筒#1の経過時間Ta(1)j が平均的な一定の増大率
でもって増大したと考えたときの経過時間AAV(1)
よりも大きいときには1番気筒#1に対するずれ割合K
AFC(1)が正の値となる。その結果補正係数KTa
(1)は、ずれ割合KAFC(1)が零となるまで減少
せしめられ、ずれ割合KAFC(1)が零になると補正
係数KTa(1)の値は一定値に落ち着くことになる。
このとき経過時間Ta(1)は経過時間AVV(1)に
一致する。全ての気筒についての補正係数KTa(i)
が一定値に落ち着くとロータ14の回転速度が一定のと
きには補正係数KTa(i)を用いて補正された各気筒
の経過時間Ta(i)は全て等しくなる。従ってロータ
14の外歯の間隔にばらつきがあったとしても各気筒の
トルク変動を正確に検出することができる。In the embodiment shown in FIG. 14, for example, the elapsed time AAV (1) when the elapsed time Ta (1) j of the first cylinder # 1 is considered to have increased at an average constant increasing rate.
Is larger than the displacement ratio K with respect to the first cylinder # 1.
AFC (1) has a positive value. As a result, the correction coefficient KTa
(1) is decreased until the deviation ratio KAFC (1) becomes zero, and when the deviation ratio KAFC (1) becomes zero, the value of the correction coefficient KTa (1) becomes a constant value.
At this time, the elapsed time Ta (1) matches the elapsed time AVV (1). Correction coefficient KTa (i) for all cylinders
Reaches a constant value, when the rotation speed of the rotor 14 is constant, the elapsed time Ta (i) of each cylinder corrected using the correction coefficient KTa (i) becomes equal. Therefore, even if the interval between the external teeth of the rotor 14 varies, the torque variation of each cylinder can be accurately detected.
【0057】再びトルク変動の算出に話しを戻すと機関
駆動系に捩り振動が発生したときに経過時間差DTa
(i)は図12に示されるように変動する。しかしなが
らこの経過時間差DTa(i)は車両が凸凹道を走行し
たときにも変動し、しかもこのときにはDTa(i)の
変動巾が極めて大きくなる場合がある。図15は車両が
凸凹道を走行したときのDTa(i)の変動を示してお
り、図15のAMPは最小のDTa(i)と最大のDT
a(i)との差、即ち振幅を示している。この振幅AM
Pが小さいときにはこれまで述べた方法によって各気筒
の出力変動およびトルク変動を正確に検出することがで
きる。Returning again to the calculation of the torque fluctuation, when the torsional vibration occurs in the engine drive system, the elapsed time difference DTa
(I) varies as shown in FIG. However, the elapsed time difference DTa (i) also fluctuates when the vehicle travels on an uneven road, and at this time, the fluctuation width of DTa (i) may become extremely large. FIG. 15 shows the variation of DTa (i) when the vehicle travels on an uneven road, and the AMP of FIG.
a (i), that is, the amplitude. This amplitude AM
When P is small, the output fluctuation and torque fluctuation of each cylinder can be accurately detected by the method described above.
【0058】しかしながら振幅AMPが大きくなるとD
Ta(i)が最大又は最小となる気筒の出力変動又はト
ルク変動を正確に検出できなくなる。即ち、図15にお
いて例えば最初にDTa(i)が最大になる気筒が1番
気筒#1であったとすると1番気筒#1の仮想の経過時
間差Kb(1)の捩り振動による減少量Iは図13のD
Ta(1)とDTa(3)とを結ぶ鎖線の傾きから求め
られる。しかしながらDTa(1)が最大となる付近で
は捩り振動による経過時間の増大量又は減少量はDTa
(2),DTa(1),DTa(3)を通る滑らかな曲
線で変化しており、従って1番気筒#1のKb(1)を
DTa(1)とDTa(3)から求めるとKb(1)の
値は実際の値よりもかなり小さく計算される。その結
果、Kb(1)が正規の値を示さなくなり、斯くして出
力変動量およびトルク変動量を正確に検出できなくな
る。振幅AMPが大きくなるとDTa(i)が最小とな
る気筒においても同じことが生ずる。However, when the amplitude AMP increases, D
The output fluctuation or torque fluctuation of the cylinder in which Ta (i) becomes maximum or minimum cannot be detected accurately. That is, in FIG. 15, for example, if the first cylinder having the maximum DTa (i) is the first cylinder # 1, the decrease I due to the torsional vibration of the virtual elapsed time difference Kb (1) of the first cylinder # 1 is shown in FIG. 13 D
It is determined from the inclination of the chain line connecting Ta (1) and DTa (3). However, near the time when DTa (1) becomes maximum, the amount of increase or decrease of the elapsed time due to torsional vibration is DTa (1).
(2), DTa (1), and a smooth curve passing through DTa (3). Therefore, when Kb (1) of the first cylinder # 1 is obtained from DTa (1) and DTa (3), Kb (1) is obtained. The value of 1) is calculated to be much smaller than the actual value. As a result, Kb (1) does not show a regular value, and thus the output fluctuation amount and the torque fluctuation amount cannot be accurately detected. When the amplitude AMP increases, the same occurs in a cylinder in which DTa (i) is minimized.
【0059】また、一つ前に燃焼が行われた気筒のDT
a(i)に対してDTa(i)が急変した気筒において
もKb(i)の値が実際の値からずれ、斯くして出力変
動およびトルク変動を正確に検出できなくなる。そこで
第1実施例では振幅AMPが大きいときにはDTa
(i)が最大又は最小となる気筒については出力変動量
又はトルク変動量を求めずに、更に一つ前に燃焼が行わ
れた気筒のDTa(i)に対してDTa(i)が急変し
た気筒についても出力変動量又はトルク変動量を求めな
いようにしている。The DT of the cylinder in which the combustion was performed immediately before
Even in a cylinder in which DTa (i) suddenly changes with respect to a (i), the value of Kb (i) deviates from the actual value, so that output fluctuation and torque fluctuation cannot be detected accurately. Therefore, in the first embodiment, when the amplitude AMP is large, DTa
For the cylinder in which (i) becomes the maximum or minimum, the output fluctuation amount or the torque fluctuation amount is not obtained, and DTa (i) suddenly changes with respect to DTa (i) of the cylinder in which the combustion was performed immediately before. The output fluctuation amount or the torque fluctuation amount is not determined for the cylinder.
【0060】次に図16から図27を参照しつつ各気筒
のトルク変動量を求めるためのルーチンについて説明す
る。なお、図27は各ルーチンにおいて行われる各値の
計算タイミングを示している。図16は30°クランク
角度毎に行われる割込みルーチンを示している。図16
を参照するとまず初めに経過時間差DTa(i)および
経過時間Tb(i)を算出するためのルーチン(ステッ
プ100)に進む。このルーチンは図17から図20に
示されている。次いでトルク変動の算出を許可するか否
かをチェックするためのルーチン(ステップ200)に
進む。このルーチンは図21から図23に示されてい
る。次いでトルク変動を算出するためのルーチン(ステ
ップ300)に進む。このルーチンは図25に示されて
いる。次いでトルク変動値の算出に用いるカウンタCD
LNIXの処理ルーチンに進む。このルーチンは図26
に示されている。Next, a routine for obtaining the torque fluctuation amount of each cylinder will be described with reference to FIGS. FIG. 27 shows the calculation timing of each value performed in each routine. FIG. 16 shows an interrupt routine performed every 30 ° crank angle. FIG.
First, the process proceeds to a routine (step 100) for calculating the elapsed time difference DTa (i) and the elapsed time Tb (i). This routine is shown in FIGS. Next, the routine proceeds to a routine (step 200) for checking whether calculation of torque fluctuation is permitted. This routine is shown in FIGS. Next, the routine proceeds to a routine (step 300) for calculating the torque fluctuation. This routine is shown in FIG. Next, a counter CD used for calculating the torque fluctuation value
The process proceeds to the LNIX processing routine. This routine is shown in FIG.
Is shown in
【0061】経過時間差DTa(i)および経過時間T
b(i)の算出ルーチンを示す図17から図20を参照
すると、まず初めにステップ101において時刻TIM
EがTIMEOとされる。電子制御ユニット20は時刻
を表わすフリーランカウンタを備えており、このフリー
ランカウンタのカウント値から時刻TIMEが算出され
る。次いでステップ102では現在の時刻TIMEが取
込まれる。従ってステップ101のTIMEOは30°
クランク角度前の時刻を表わしていることになる。Elapsed time difference DTa (i) and elapsed time T
Referring to FIGS. 17 to 20 showing the routine for calculating b (i), first, at step 101, the time TIM
E is set to TIMEO. The electronic control unit 20 includes a free-run counter that indicates the time, and the time TIME is calculated from the count value of the free-run counter. Next, at step 102, the current time TIME is taken. Therefore, the TIMEO of step 101 is 30 °
This indicates the time before the crank angle.
【0062】次いでステップ103では現在i番気筒の
ATDC30°であるか否かが判別される。現在i番気
筒のATDC30°でない場合にはステップ111にジ
ャンプして現在i番気筒のATDC90°であるか否か
が判別される。現在i番気筒のATDC90°でない場
合には経過時間差DTa(i)および経過時間Tb
(i)の算出ルーチンを完了する。Next, at step 103, it is determined whether the i-th cylinder is at ATDC 30 ° or not. If the ATDC of the i-th cylinder is not 30 °, the process jumps to step 111 to determine whether or not the ATDC of the i-th cylinder is 90 °. If the current i-th cylinder is not at 90 ° ATDC, the elapsed time difference DTa (i) and the elapsed time Tb
The calculation routine of (i) is completed.
【0063】これに対してステップ103において現在
i番気筒のATDC30°であると判別されたときには
ステップ104に進んで720°クランク角度前に算出
された経過時間Ta(i)がTaO(i)とされる。次
いでステップ105では次式に基づいてi番気筒のTD
CからATDC30°までの最終的な経過時間Ta
(i)が算出される。On the other hand, when it is determined in step 103 that the ATDC of the i-th cylinder is 30 °, the routine proceeds to step 104, where the elapsed time Ta (i) calculated before the 720 ° crank angle is TaO (i). Is done. Next, at step 105, the TD of the i-th cylinder is calculated based on the following equation.
Final elapsed time Ta from C to ATDC 30 °
(I) is calculated.
【0064】Ta(i)=(TIME−TIMEO)・
(1+KTa(i)) 即ち、例えば現在1番気筒#1のATDC30°である
とすると1番気筒#1のTDCからATDC30°まで
の最終的な経過時間Ta(1)が(TIME−TIME
O)・(1+KTa(1))から算出される。ここで
(TIME−TIMEO)はクランク角センサ15によ
り実測された経過時間Ta(1)を表わしており、KT
a(1)はロータ13の外歯間隔による誤差を補正する
ための補正係数であり、従って(TIME−TIME
O)に(1+KTa(1))を乗算することによって得
られた最終的な経過時間Ta(1)はクランクシャフト
が30°クランク角度回転する間の経過時間を正確に表
わしていることになる。Ta (i) = (TIME-TIMEO).
(1 + KTa (i)) That is, for example, if it is the ATDC 30 ° of the first cylinder # 1 now, the final elapsed time Ta (1) from TDC of the first cylinder # 1 to ATDC 30 ° is (TIME−TIME)
O) · (1 + KTa (1)). Here, (TIME-TIMEO) represents the elapsed time Ta (1) actually measured by the crank angle sensor 15, and KT
a (1) is a correction coefficient for correcting an error due to the external tooth interval of the rotor 13, and accordingly, (TIME-TIME)
The final elapsed time Ta (1) obtained by multiplying (O) by (1 + KTa (1)) accurately represents the elapsed time during which the crankshaft rotates by 30 ° crank angle.
【0065】次いでステップ106では今回算出された
経過時間Ta(i)から720°クランク角度前に算出
された経過時間TaO(i)を減算することによって経
過時間差DTa(i)(=Ta(i)−TaO(i))
が算出される。次いでステップ107では現在2番気筒
#2のATDC30°であるか否かが判別される。現在
2番気筒#2のATDC30°でないときにはステップ
110にジャンプし、一つ前に燃焼が行われた(i−
1)番気筒のトルク変動量を算出すべきことを示すフラ
グXCAL(i−1)がセット(XCAL(i−1)←
“1”)される。次いでステップ111に進む。第1実
施例では前述したように点火順序が1−3−4−2であ
るので現在1番気筒#1のATDC30°であるとする
と一つ前に燃焼が行われた2番気筒#2のトルク変動量
を算出すべきことを示すフラグXCAL(2)がセット
される。同様に図27に示される如く最終的な経過時間
Ta(3)が算出されるとフラグXCAL(1)がセッ
トされ、最終的な経過時間Ta(4)が算出されるとフ
ラグXCAL(3)がセットされ、最終的な経過時間T
a(2)が算出されるとフラグXCAL(4)がセット
される。Next, at step 106, the elapsed time difference DTa (i) (= Ta (i) is obtained by subtracting the elapsed time TaO (i) calculated 720 ° crank angle before from the currently calculated elapsed time Ta (i). -TaO (i))
Is calculated. Next, at step 107, it is determined whether or not the ATDC of the second cylinder # 2 is 30 °. When the ATDC of the second cylinder # 2 is not 30 °, the process jumps to step 110, and the combustion is performed immediately before (i-
1) A flag XCAL (i-1) indicating that the torque variation of the cylinder No. should be calculated is set (XCAL (i-1) ←
"1"). Next, the routine proceeds to step 111. In the first embodiment, as described above, since the ignition order is 1-3-4-2, assuming that the ATDC of the first cylinder # 1 is 30 ° at present, the ignition of the second cylinder # 2 in which the combustion was performed immediately before is performed. A flag XCAL (2) indicating that the torque variation should be calculated is set. Similarly, as shown in FIG. 27, the flag XCAL (1) is set when the final elapsed time Ta (3) is calculated, and the flag XCAL (3) is calculated when the final elapsed time Ta (4) is calculated. Is set, and the final elapsed time T
When a (2) is calculated, the flag XCAL (4) is set.
【0066】一方、ステップ111において現在i番気
筒のATDC90°であると判別されたときにはステッ
プ112に進んで720°クランク角度前に算出された
経過時間Tb(i)がTbO(i)とされる。次いでス
テップ113では次式に基づいてi番気筒のATDC6
0°からATDC90°までの最終的な経過時間Tb
(i)が算出される。On the other hand, when it is determined in step 111 that the ATDC of the i-th cylinder is 90 °, the routine proceeds to step 112, where the elapsed time Tb (i) calculated before the 720 ° crank angle is set to TbO (i). . Next, at step 113, the ATDC6 of the i-th cylinder is determined based on the following equation.
Final elapsed time Tb from 0 ° to 90 ° ATDC
(I) is calculated.
【0067】Tb(i)=(TIME−TIMEO)・
(1+KTb(i)) 即ち、例えば現在1番気筒#1のATDC90°である
とすると1番気筒#1のATDC60°からATDC9
0°までの最終的な経過時間Tb(1)が(TIME−
TIMEO)・(1+KTb(1))から算出される。
この場合にもロータ13の外歯間隔による誤差を補正す
るための値(1+KTb(1))が(TIME−TIM
EO)に乗算されているので最終的な経過時間Tb
(1)はクランクシャフトが30°クランク角度回転す
る間の経過時間を正確に表わしていることになる。次い
でステップ114では現在2番気筒#2のATDC90
°であるか否かが判別される。現在2番気筒#2のAT
DC90°でないときには経過時間差DTa(i)およ
び経過時間Tb(i)の算出ルーチンを完了する。Tb (i) = (TIME-TIMEO).
(1 + KTb (i)) That is, for example, if the ATDC of the first cylinder # 1 is currently 90 °, the ATDC of the first cylinder # 1 is 60 ° to ATDC9.
The final elapsed time Tb (1) up to 0 ° is (TIME-
TIMEO) · (1 + KTb (1)).
Also in this case, the value (1 + KTb (1)) for correcting the error due to the external tooth interval of the rotor 13 is (TIME-TIM).
EO), the final elapsed time Tb
(1) accurately represents the elapsed time during which the crankshaft rotates by 30 ° crank angle. Next, at step 114, the ATDC 90 of the current second cylinder # 2
° is determined. AT of current cylinder # 2
If DC is not 90 °, the calculation routine of the elapsed time difference DTa (i) and the elapsed time Tb (i) is completed.
【0068】一方、ステップ107において現在2番気
筒#2のATDC30°であると判別されたときにはス
テップ108に進んで次式に基づき図14の720°ク
ランク角度範囲(j)における経過時間Ta(i)の平
均値TaAVj が算出される。 TaAVj =(Ta(1)+Ta(3)+Ta(4)+
Ta(2))/4 ここでTa(1),Ta(3),Ta(4),Ta
(2)は図14におけるTa(1)j ,Ta(3)j ,
Ta(4)j ,Ta(2)j に夫々相当する。次いでス
テップ109の補正係数KTa(i)の算出ルーチンに
進む。このルーチンは図19に示されている。On the other hand, when it is determined in step 107 that the ATDC of the second cylinder # 2 is 30 ° at present, the routine proceeds to step 108, where the elapsed time Ta (i) in the 720 ° crank angle range (j) of FIG. average TaAV j of) is calculated. TaAV j = (Ta (1) + Ta (3) + Ta (4) +
Ta (2)) / 4 where Ta (1), Ta (3), Ta (4), Ta
(2) is Ta (1) j , Ta (3) j ,
Ta (4) j and Ta (2) j respectively. Next, the routine proceeds to a routine for calculating the correction coefficient KTa (i) in step 109. This routine is shown in FIG.
【0069】図19を参照するとまず初めにステップ1
20において減速運転中の燃料供給停止時であるか否か
が判別される。減速運転中の燃料供給停止時でないとき
には処理サイクルを完了し、減速運転中の燃料供給停止
時にはステップ121に進む。ステップ121ではステ
ップ108において算出された経過時間Ta(i)の平
均値TaAVj と、既に算出されている図14の720
°クランク角度範囲(j−1)における経過時間Ta
(i)の平均値TaAVj-1 (=(Ta(1)+Ta
(3)+Ta(4)+Ta(2))/4)(ここでTa
(1),Ta(3),Ta(4),Ta(2)は夫々図
14に示されるTa(1)j-1 ,Ta(3) j-1 ,Ta
(4)j-1 ,Ta(2)j-1 に相当する)から、平均的
な一定の増大率でもって増大したと考えたときの図14
に示す経過時間AAV(4),AAV(2),AAV
(1),AAV(3)が次のようにして算出される。Referring to FIG. 19, first, in step 1
In 20 whether or not the fuel supply is stopped during deceleration operation
Is determined. When fuel supply is not stopped during deceleration operation
Completes the processing cycle and stops fuel supply during deceleration operation
Sometimes the process proceeds to step 121. Step 121
Of the elapsed time Ta (i) calculated in step 108
Average value TaAVjAnd 720 of FIG.
° Elapsed time Ta in crank angle range (j-1)
Average value TaAV of (i)j-1 (= (Ta (1) + Ta
(3) + Ta (4) + Ta (2)) / 4) (where Ta
(1), Ta (3), Ta (4), Ta (2)
Ta (1) shown in 14j-1 , Ta (3) j-1 , Ta
(4)j-1 , Ta (2)j-1 Is equivalent to)
FIG. 14 when it is considered that it has increased at a certain constant increase rate
Elapsed time AAV (4), AAV (2), AAV
(1), AAV (3) is calculated as follows.
【0070】AAV(4)=(TaAVj −TaAV
j-1 )・(1/8)+TaAVj-1 AAV(2)=(TaAVj −TaAVj-1 )・(3/
8)+TaAVj-1 AAV(1)=(TaAVj −TaAVj-1 )・(5/
8)+TaAVj-1 AAV(3)=(TaAVj −TaAVj-1 )・(7/
8)+TaAVj-1 次いでステップ122では次式に基づいてこれら経過時
間AAV(4),AAV(2),AAV(1),AAV
(3)に対する実際に検出された経過時間TaO
(4),TaO(2),Ta(1),Ta(3)のずれ
割合KAFC(4),KAFC(2),KAFC
(1),KAFC(3)が算出される。ここでTaO
(4),TaO(2),Ta(1),Ta(3)は夫々
図14に示されるTa(4)j-1 ,Ta(2)j-1 ,T
a(1)j ,Ta(3)j に相当する。AAV (4) = (TaAV j −TaAV
j-1 ). (1/8) + TaAV j-1 AAV (2) = (TaAV j -TaAV j-1 ). (3 /
8) + TaAV j−1 AAV (1) = (TaAV j −TaAV j−1 ) · (5 /
8) + TaAV j−1 AAV (3) = (TaAV j −TaAV j−1 ) · (7 /
8) + TaAV j-1 Next, at step 122, these elapsed times AAV (4), AAV (2), AAV (1), AAV
Elapsed time TaO actually detected for (3)
(4), TaO (2), Ta (1), Ta (3) deviation ratio KAFC (4), KAFC (2), KAFC
(1), KAFC (3) is calculated. Here TaO
(4), TaO (2), Ta (1), Ta (3) are Ta (4) j−1 , Ta (2) j−1 , T
a (1) j and Ta (3) j .
【0071】KAFC(4)={TaO(4)−AAV
(4)}/AAV(4) KAFC(2)={TaO(2)−AAV(2)}/A
AV(2) KAFC(1)={Ta(1)−AAV(1)}/AA
V(1) KAFC(3)={Ta(3)−AAV(3)}/AA
V(3) 次いでステップ123では補正係数KTa(i)の学習
条件が成立しているか否かが判別される。例えば燃料供
給停止から一定時間経過しており、機関回転数が予め定
められた範囲内にあり、かつ車両が悪路を走行していな
いとき、例えば図15に示される経過時間差DTa
(i)の振幅AMPが設定値を越えていないときに学習
条件が成立していると判断される。学習条件が成立して
いないときにはステップ128に進んでずれ量KAFC
(i)の積算値KAFCI(i)が零とされ、次いでス
テップ129において積算カウント値CKAFCがクリ
アされる。KAFC (4) = {TaO (4) -AAV
(4)} / AAV (4) KAFC (2) = {TaO (2) −AAV (2)} / A
AV (2) KAFC (1) = {Ta (1) -AAV (1)} / AA
V (1) KAFC (3) = {Ta (3) -AAV (3)} / AA
V (3) Next, in step 123, it is determined whether or not the learning condition of the correction coefficient KTa (i) is satisfied. For example, when a certain time has elapsed since the fuel supply was stopped, the engine speed is within a predetermined range, and the vehicle is not traveling on a rough road, for example, the elapsed time difference DTa shown in FIG.
When the amplitude AMP in (i) does not exceed the set value, it is determined that the learning condition is satisfied. If the learning condition is not satisfied, the routine proceeds to step 128, where the deviation amount KAFC
The integrated value KAFCI (i) of (i) is set to zero, and then at step 129, the integrated count value CKAFC is cleared.
【0072】これに対して学習条件が成立している場合
にはステップ124に進んで各気筒に対するずれ量KA
FC(i)が対応する積算値KAFCI(i)に加算さ
れ、次いでステップ125において積算カウント値CK
AFCが1だけインクリメントされる。次いでステップ
126では積算カウント値CKAFCが8になったか否
かが判別される。積算カウント値CKAFCが8でない
ときには処理サイクルを完了し、積算カウント値CKA
FCが8になるとステップ127に進んで補正係数KT
a(i)が算出される。即ち、各気筒についてずれ量K
AFC(i)が8回積算されるとステップ127に進ん
で補正係数KTa(i)が算出される。On the other hand, if the learning condition is satisfied, the routine proceeds to step 124, where the deviation KA for each cylinder is determined.
FC (i) is added to the corresponding integrated value KAFCI (i), and then at step 125, the integrated count value CK
AFC is incremented by one. Next, at step 126, it is determined whether or not the integrated count value CKAFC has reached 8. When the accumulated count value CKAFC is not 8, the processing cycle is completed, and the accumulated count value CKA
When FC becomes 8, the routine proceeds to step 127, where the correction coefficient KT
a (i) is calculated. That is, the deviation amount K for each cylinder
When AFC (i) is integrated eight times, the routine proceeds to step 127, where the correction coefficient KTa (i) is calculated.
【0073】ステップ127では次のようにして補正係
数KTa(i)が算出される。即ち、まず初めに各積算
値KAFCI(i)がKAFCE(i)に置き換えられ
る。次いで次式に基づき基本補正係数KTaB(i)が
算出される。 KTaB(i)=KTa(i)−(KAFCI(i)/
8)/4 即ち、実際に検出された経過時間Ta(i)が平均的な
一定の増大率でもって増大したと考えたときの経過時間
AAV(i)に対してずれを生じている場合には現在の
補正係数KTa(i)を積算値KAFCI(i)の平均
値の1/4だけ修正した値が基本補正係数KTaB
(i)とされる。次いで次式に基づき全気筒に対する基
本補正係数KTaB(i)の平均値tKTaMが算出さ
れる。In step 127, the correction coefficient KTa (i) is calculated as follows. That is, first, each integrated value KAFCI (i) is replaced with KAFCE (i). Next, a basic correction coefficient KTaB (i) is calculated based on the following equation. KTaB (i) = KTa (i)-(KAFCI (i) /
8) / 4 That is, when the elapsed time Ta (i) actually detected is shifted from the elapsed time AAV (i) when it is considered that the elapsed time Ta (i) increases at an average constant increase rate. Is the basic correction coefficient KTaB obtained by correcting the current correction coefficient KTa (i) by 4 of the average value of the integrated value KAFCI (i).
(I). Next, the average value tKTaM of the basic correction coefficient KTaB (i) for all cylinders is calculated based on the following equation.
【0074】tKTaM=(KTaB(1)+KTaB
(2)+KTaB(3)+KTaB(4))/4 次いで次式に示されるように基本補正係数KTaB
(i)をその平均値tKTaMでもって修正することに
より補正係数KTa(i)が更新される。 KTa(i)=KTaB(i)−tKTaM このようにして減速運転中の燃料噴射停止時に各気筒に
対する補正係数KTa(i)の更新が行われる。TKTaM = (KTaB (1) + KTaB
(2) + KTaB (3) + KTaB (4)) / 4 Then, as shown in the following equation, the basic correction coefficient KTaB
By correcting (i) with its average value tKTaM, the correction coefficient KTa (i) is updated. KTa (i) = KTaB (i) -tKTaM In this way, the correction coefficient KTa (i) for each cylinder is updated when fuel injection is stopped during deceleration operation.
【0075】一方、図18のステップ114において現
在2番気筒#2のATDC90°であると判別されたと
きにはステップ115に進んでTaAVj を求めたとき
と同じ方法により次式に基づいて720°クランク角度
範囲における経過時間Tb(i)の平均値TbAVj が
算出される。 TbAVj =(Tb(1)+Tb(3)+Tb(4)+
Tb(2))/4 次いでステップ116の補正係数KTb(i)の算出ル
ーチンに進む。このルーチンは図20に示されている。On the other hand, if it is determined in step 114 of FIG. 18 that the ATDC of the second cylinder # 2 is 90 ° at present, the routine proceeds to step 115, and the 720 ° crank angle is obtained by the same method as used when TaAV j is obtained, based on the following equation. average TbAV j of the elapsed time in the angular range Tb (i) is calculated. TbAV j = (Tb (1) + Tb (3) + Tb (4) +
Tb (2)) / 4 Next, the routine proceeds to a routine for calculating the correction coefficient KTb (i) in step 116. This routine is shown in FIG.
【0076】図20を参照するとまず初めにステップ1
30において減速運転中の燃料供給停止時であるか否か
が判別される。減速運転中の燃料供給停止時でないとき
には処理サイクルを完了し、減速運転中の燃料供給停止
時にはステップ131に進む。ステップ131ではステ
ップ115において算出された経過時間Tb(i)の平
均値TbAVj と、一つ前の720°クランク角度範囲
における経過時間Tb(i)の平均値TbAVj-1 (=
(Tb(1)+Tb(3)+Tb(4)+Tb(2))
/4)から、平均的な一定の増大率でもって増大したと
考えたときの経過時間BAV(4),BAV(2),B
AV(1),BAV(3)が算出される。Referring to FIG. 20, first, in step 1
At 30, it is determined whether or not the fuel supply is stopped during the deceleration operation. When the fuel supply is not stopped during the deceleration operation, the processing cycle is completed. When the fuel supply is stopped during the deceleration operation, the process proceeds to step 131. The average value TbAV j in step 131 the elapsed time Tb which is calculated in step 115 (i), the average value of the time elapsed in 720 ° crank angle range immediately preceding Tb (i) TbAV j-1 (=
(Tb (1) + Tb (3) + Tb (4) + Tb (2))
/ 4), the elapsed times BAV (4), BAV (2), B
AV (1) and BAV (3) are calculated.
【0077】BAV(4)=(TbAVj −TbAV
j-1 )・(1/8)+TbAVj-1 BAV(2)=(TbAVj −TbAVj-1 )・(3/
8)+TbAVj-1 BAV(1)=(TbAVj −TbAVj-1 )・(5/
8)+TbAVj-1 BAV(3)=(TbAVj −TbAVj-1 )・(7/
8)+TbAVj-1 次いでステップ132では次式に基づいてこれら経過時
間BAV(4),BAV(2),BAV(1),BAV
(3)に対する実際に検出された経過時間TbO
(4),TbO(2),Tb(1),Tb(3)のずれ
割合KBFC(4),KBFC(2),KBFC
(1),KBFC(3)が算出される。BAV (4) = (TbAV j −TbAV
j-1 ). (1/8) + TbAV j-1 BAV (2) = (TbAV j -TbAV j-1 ). (3 /
8) + TbAV j−1 BAV (1) = (TbAV j −TbAV j−1 ) · (5 /
8) + TbAV j−1 BAV (3) = (TbAV j −TbAV j−1 ) · (7 /
8) + TbAV j-1 Next, at step 132, these elapsed times BAV (4), BAV (2), BAV (1), BAV
Elapsed time TbO actually detected for (3)
(4), deviation ratio of TbO (2), Tb (1), Tb (3) KBFC (4), KBFC (2), KBFC
(1), KBFC (3) is calculated.
【0078】KBFC(4)={TbO(4)−BAV
(4)}/BAV(4) KBFC(2)={TbO(2)−BAV(2)}/B
AV(2) KBFC(1)={Tb(1)−BAV(1)}/BA
V(1) KBFC(3)={Tb(3)−BAV(3)}/BA
V(3) 次いでステップ133では補正係数KTb(i)の学習
条件が成立しているか否かが判別される。前述したよう
に、例えば燃料供給停止から一定時間経過しており、機
関回転数が予め定められた範囲内にあり、かつ車両が悪
路を走行していないとき、例えば図15に示される経過
時間差DTa(i)の振幅AMPが設定値を越えていな
いときに学習条件が成立していると判断される。学習条
件が成立していないときにはステップ138に進んでず
れ量KBFC(i)の積算値KBFCI(i)が零とさ
れ、次いでステップ139において積算カウント値CK
BFCがクリアされる。KBFC (4) = {TbO (4) -BAV
(4)} / BAV (4) KBFC (2) = {TbO (2) -BAV (2)} / B
AV (2) KBFC (1) = {Tb (1) -BAV (1)} / BA
V (1) KBFC (3) = {Tb (3) -BAV (3)} / BA
V (3) Next, at step 133, it is determined whether or not the learning condition of the correction coefficient KTb (i) is satisfied. As described above, for example, when a certain time has elapsed since the fuel supply was stopped, the engine speed is within a predetermined range, and the vehicle is not traveling on a bad road, for example, the elapsed time difference shown in FIG. When the amplitude AMP of DTa (i) does not exceed the set value, it is determined that the learning condition is satisfied. If the learning condition is not satisfied, the routine proceeds to step 138, where the integrated value KBFCI (i) of the deviation amount KBFC (i) is made zero, and then at step 139, the integrated count value CK
BFC is cleared.
【0079】これに対して学習条件が成立している場合
にはステップ134に進んで各気筒に対するずれ量KB
FC(i)が対応する積算値KBFCI(i)に加算さ
れ、次いでステップ135において積算カウント値CK
BFCが1だけインクリメントされる。次いでステップ
136では積算カウント値CKBFCが8になったか否
かが判別される。積算カウント値CKBFCが8でない
ときには処理サイクルを完了し、積算カウント値CKB
FCが8になるとステップ137に進んで補正係数KT
b(i)が算出される。即ち、各気筒についてずれ量K
BFC(i)が8回積算されるとステップ137に進ん
で補正係数KTb(i)が算出される。On the other hand, if the learning condition is satisfied, the routine proceeds to step 134, where the deviation amount KB for each cylinder is determined.
FC (i) is added to the corresponding integrated value KBFCI (i), and then at step 135, the integrated count value CK
BFC is incremented by one. Next, at step 136, it is determined whether or not the integrated count value CKBFC has reached 8. When the accumulated count value CKBFC is not 8, the processing cycle is completed, and the accumulated count value CKB
When FC becomes 8, the routine proceeds to step 137, where the correction coefficient KT
b (i) is calculated. That is, the deviation amount K for each cylinder
When BFC (i) is integrated eight times, the routine proceeds to step 137, where the correction coefficient KTb (i) is calculated.
【0080】ステップ137では次のようにして補正係
数KTb(i)が算出される。即ち、まず初めに各積算
値KBFCI(i)がKBFCE(i)に置き換えられ
る。次いで次式に基づき基本補正係数KTbB(i)が
算出される。 KTbB(i)=KTb(i)−(KBFCI(i)/
8)/4 即ち、実際に検出された経過時間Tb(i)が平均的な
一定の増大率でもって増大したと考えたときの経過時間
BAV(i)に対してずれを生じている場合には現在の
補正係数KTb(i)を積算値KBFCI(i)の平均
値の1/4だけ修正した値が基本補正係数KTbB
(i)とされる。次いで次式に基づき全気筒に対する基
本補正係数KTbB(i)の平均値tKTbMが算出さ
れる。In step 137, the correction coefficient KTb (i) is calculated as follows. That is, first, each integrated value KBFCI (i) is replaced with KBFCE (i). Next, the basic correction coefficient KTbB (i) is calculated based on the following equation. KTbB (i) = KTb (i)-(KBFCI (i) /
8) / 4 That is, when the actually detected elapsed time Tb (i) is shifted from the elapsed time BAV (i) when it is considered that the elapsed time Tb (i) increases at an average constant increase rate. Is the basic correction coefficient KTbB obtained by correcting the current correction coefficient KTb (i) by 4 of the average value of the integrated value KBFCI (i).
(I). Next, an average value tKTbM of the basic correction coefficients KTbB (i) for all cylinders is calculated based on the following equation.
【0081】tKTbM=(KTbB(1)+KTbB
(2)+KTbB(3)+KTbB(4))/4 次いで次式に示されるように基本補正係数KTbB
(i)をその平均値tKTbMでもって修正することに
より補正係数KTb(i)が更新される。 KTb(i)=KTbB(i)−tKTbM このようにして減速運転中の燃料噴射停止時に各気筒に
対する補正係数KTb(i)の更新が行われる。TKTbM = (KTbB (1) + KTbB
(2) + KTbB (3) + KTbB (4)) / 4 Then, as shown in the following equation, the basic correction coefficient KTbB
By correcting (i) with its average value tKTbM, the correction coefficient KTb (i) is updated. KTb (i) = KTbB (i) -tKTbM In this manner, when the fuel injection is stopped during the deceleration operation, the correction coefficient KTb (i) for each cylinder is updated.
【0082】次に図21から図23に示されるトルク変
動算出許可チェックルーチンについて図24を参照しつ
つ説明する。このルーチンは車両が凸凹道を走行するこ
とにより経過時間差DTa(i)の変動の振幅AMP
(図15)が大きくなったときには特定の気筒について
のトルク変動量の算出を禁止するために設けられてい
る。Next, the torque variation calculation permission check routine shown in FIGS. 21 to 23 will be described with reference to FIG. In this routine, the amplitude AMP of the fluctuation of the elapsed time difference DTa (i) due to the vehicle traveling on the uneven road
It is provided to prohibit the calculation of the torque fluctuation amount for a specific cylinder when (FIG. 15) becomes large.
【0083】即ち、図21から図23を参照すると、ま
ず初めにステップ201において現在いずれかの気筒の
ATDC30°であるか否かが判別される。現在いずれ
かの気筒のATDC30°でないときには処理サイクル
を完了し、現在いずれかの気筒のATDC30°である
ときにはステップ202に進む。ステップ202からス
テップ204では経過時間差DTa(i)が増大し次い
で減少する際の最大経過時間差DT30maxが算出さ
れる。即ち、ステップ202では図17のステップ10
6において算出されたDTa(i)が最大経過時間差D
T30maxよりも大きいか否かが判別される。DT3
0max>DTa(i)のときにはステップ205にジ
ャンプし、これに対してDT30max≦DTa(i)
のときにはステップ203に進んでDTa(i)がDT
30maxとされる。次いでステップ204ではDTa
(i)が増大していることを示す増大フラグXMXRE
Cがセット(XMXREC←“1”)され、次いでステ
ップ205に進む。That is, referring to FIGS. 21 to 23, first, at step 201, it is determined whether or not the ATDC of any of the cylinders is currently 30 °. If the ATDC of any of the cylinders is not 30 °, the processing cycle is completed. If the ATDC of any of the cylinders is currently 30 °, the process proceeds to step 202. In steps 202 to 204, the maximum elapsed time difference DT30max when the elapsed time difference DTa (i) increases and then decreases is calculated. That is, in step 202, step 10 in FIG.
6, the DTa (i) calculated is the maximum elapsed time difference D
It is determined whether it is larger than T30max. DT3
If 0max> DTa (i), the routine jumps to step 205, whereas DT30max ≦ DTa (i)
In the case of, the routine proceeds to step 203, where DTa (i) is DT
30 max. Next, at step 204, DTa
An increase flag XMXRE indicating that (i) is increasing
C is set (XMXREC ← “1”), and then the routine proceeds to step 205.
【0084】ステップ205からステップ207では経
過時間差DTa(i)が減少し次いで増大する際の最小
経過時間差DT30minが算出される。即ち、ステッ
プ205では図17のステップ106において算出され
たDTa(i)が最小経過時間差DT30minよりも
小さいか否かが判別される。DT30min<DTa
(i)のときにはステップ208にジャンプし、これに
対してDT30min≧DTa(i)のときにはステッ
プ206に進んでDTa(i)がDT30minとされ
る。次いでステップ207ではDTa(i)が減少して
いることを示す減少フラグXMNRECがセット(XM
NREC←“1”)され、次いでステップ208に進
む。In steps 205 to 207, the minimum elapsed time difference DT30min when the elapsed time difference DTa (i) decreases and then increases is calculated. That is, in step 205, it is determined whether or not DTa (i) calculated in step 106 in FIG. 17 is smaller than the minimum elapsed time difference DT30min. DT30min <DTa
In the case of (i), the process jumps to step 208. On the other hand, when DT30min ≧ DTa (i), the process proceeds to step 206, where DTa (i) is set to DT30min. Next, at step 207, a decrease flag XMNREC indicating that DTa (i) is decreasing is set (XM
NREC ← “1”), and then proceed to step 208.
【0085】ステップ208からステップ214ではD
Ta(i)の変動の振幅AMP(図15)が設定値A0
を越えたときにはDTa(i)が最大となった気筒につ
いてのトルク変動量の算出を禁止する禁止フラグがセッ
トされる。即ち、ステップ208ではDT30max>
DTa(i)でかつXMXREC=“1”であるか否か
が判別される。DT30max≦DTa(i)である
か、又は増大フラグXMXRECがリセット(XMXR
EC=“0”)されているときにはステップ215にジ
ャンプし、これに対してDT30max>DTa(i)
でかつXMXREC=“1”のときにはステップ209
に進む。In steps 208 to 214, D
The amplitude AMP of the fluctuation of Ta (i) (FIG. 15) is equal to the set value A 0.
Is exceeded, a prohibition flag is set to prohibit the calculation of the amount of torque fluctuation for the cylinder in which DTa (i) is maximized. That is, in step 208, DT30max>
It is determined whether or not DTa (i) and XMXREC = "1". DT30max ≦ DTa (i) or the increase flag XMXREC is reset (XMXR
If EC = “0”), the process jumps to step 215, where DT30max> DTa (i)
And if XMXREC = "1", step 209
Proceed to.
【0086】即ち、図24に示されるように時刻t1 に
おいて1番気筒#1の経過時間差DTa(1)が最大に
なったとする。この場合、時刻t1 において行われる割
込みルーチンではステップ202からステップ203に
進んでDTa(1)がDT30maxとされ、次いでス
テップ204において最大フラグXMXRECがセット
される。一方、図24の時刻t2 において行われる割込
みルーチンではステップ202からステップ205にジ
ャンプする。このときステップ208ではDT30ma
x>DTa(3)であり、かつXMXREC=“1”で
あると判断されるのでステップ209に進む。即ち、ス
テップ209に進むのは経過時間差DTa(i)が減少
しはじめる時刻t2 である。[0086] That is, the No. 1 cylinder # 1 of the difference of elapsed times DTa at time t 1 as shown in FIG. 24 (1) is to become the maximum. In this case, in the interruption routine performed at the time t 1 proceeds from step 202 to step 203 DTa (1) is a DT30max, then the maximum flag XMXREC is set in step 204. On the other hand, in the interruption routine performed at the time t 2 in FIG. 24 jumps from step 202 to step 205. At this time, in step 208, DT30ma
Since it is determined that x> DTa (3) and XMXREC = "1", the flow proceeds to step 209. That is, the time t 2 the flow proceeds to step 209 where the elapsed time DTa (i) starts to decrease.
【0087】ステップ209では最大経過時間差DT3
0maxがTMXRECとされる。次いでステップ21
0では最大経過時間差TMXRECから最小経過時間差
TMNREC(後述するステップ216で求められる)
を減算することによってDTa(i)の変動の振幅AM
Pが算出される。次いでステップ211では最小経過時
間差DT30minの初期値がDTa(i)とされる。
次いでステップ212では増大フラグXMXRECがリ
セット(XMXREC←“0”)される。次いでステッ
プ213では振幅AMPが設定値A0 よりも大きいか否
かが判別される。AMP<A0 のときにはステップ21
5にジャンプする。これに対してAMP≧A0 のときに
はステップ214に進んでトルク変動算出禁止フラグX
NOCALがセット(XNOCAL←“1”)される。
即ち、図24の時刻t2 において行われる割込みルーチ
ンでは前述したように1番気筒#1のトルク変動量が算
出される。従ってこの割込みルーチンにおいてAMP≧
A0 となり、トルク変動算出禁止フラグXNOCALが
セットされると1番気筒#1のトルク変動量の算出、即
ち、DTa(i)が最大となる気筒のトルク変動量の算
出が禁止される。In step 209, the maximum elapsed time difference DT3
0max is set to TMXREC. Then step 21
At 0, the minimum elapsed time difference TMNREC from the maximum elapsed time difference TMXREC (determined in step 216 described later)
Is subtracted to obtain the amplitude AM of the variation of DTa (i).
P is calculated. Next, at step 211, the initial value of the minimum elapsed time difference DT30min is set to DTa (i).
Next, at step 212, the increase flag XMXREC is reset (XMXREC ← “0”). Next, at step 213 whether the amplitude AMP is greater than the set value A 0 is judged. Step 21 when AMP <A 0
Jump to 5. On the other hand, when AMP ≧ A 0, the routine proceeds to step 214, where the torque fluctuation calculation prohibition flag X
NOCAL is set (XNOCAL ← “1”).
That is, 1st amount of torque fluctuation of the cylinders # 1, as described above in the interruption routine performed at the time t 2 in FIG. 24 is calculated. Therefore, in this interrupt routine, AMP ≧
A 0 becomes, when the torque fluctuation calculation prohibition flag XNOCAL is set determining of the first cylinder # 1 of the amount of torque fluctuation, i.e., DTa (i) to calculate the amount of torque fluctuation of the cylinder giving the maximum is prohibited.
【0088】ステップ215からステップ221ではD
Ta(i)の変動の振幅AMPが設定値A0 を越えたと
きにはDTa(i)が最小となった気筒についてのトル
ク変動量の算出を禁止する禁止フラグがセットされる。
即ち、ステップ215ではDT30min<DTa
(i)でかつXMNREC=“1”であるか否かが判別
される。DT30min≧DTa(i)であるか、又は
減少フラグXMNRECがリセット(XMNREC=
“0”)されているときにはステップ222にジャンプ
し、これに対してDT30min<DTa(i)でかつ
XMNREC=“1”のときにはステップ216に進
む。In steps 215 to 221, D
Amplitude AMP of the variation of ta (i) is set prohibition flag for prohibiting calculation of the amount of torque fluctuation for cylinder DTa (i) is minimized when it exceeds the set value A 0.
That is, in step 215, DT30min <DTa
It is determined whether or not (i) and XMNREC = "1". DT30min ≧ DTa (i) or the decrease flag XMNREC is reset (XMREC =
If "0"), the process jumps to step 222, whereas if DT30min <DTa (i) and XMNREC = "1", the process proceeds to step 216.
【0089】即ち、図24に示されるように時刻t3 に
おいて1番気筒#1の経過時間差DTa(1)が最小に
なったとする。この場合、時刻t3 において行われる割
込みルーチンではステップ205からステップ206に
進んでDTa(1)がDT30minとされ、次いでス
テップ207において減少フラグXMNRECがセット
される。一方、図24の時刻t4 において行われる割込
みルーチンではステップ205からステップ208にジ
ャンプする。このときステップ215ではDT30mi
n<DTa(3)であり、かつXMNREC=“1”で
あると判断されるのでステップ216に進む。即ち、ス
テップ216に進むのは経過時間差DTa(i)が増大
しはじめる時刻t4 である。[0089] That is, the first cylinder # 1 of the difference of elapsed times DTa at time t 3 as shown in FIG. 24 (1) is minimized. In this case, in the interruption routine performed at the time t 3 proceeds from step 205 to step 206 DTa (1) is a DT30min, then decrease flag XMNREC is set in step 207. On the other hand, in the interruption routine performed at the time t 4 in FIG. 24 jumps from step 205 to step 208. At this time, in step 215, DT30mi
Since it is determined that n <DTa (3) and XMNREC = "1", the process proceeds to step 216. That is, the time t 4 when the elapsed time difference DTa proceed to step 216 (i) starts to increase.
【0090】ステップ216では最小経過時間差DT3
0minがTMNRECとされる。次いでステップ21
7では最大経過時間差TMXRECから最小経過時間差
TMNRECを減算することによってDTa(i)の変
動の振幅AMPが算出される。次いでステップ218で
は最大経過時間差DT30maxの初期値がDTa
(i)とされる。次いでステップ219では減少フラグ
XMNRECがリセット(XMNREC←“0”)され
る。次いでステップ220では振幅AMPが設定値A0
よりも大きいか否かが判別される。AMP<A0 のとき
にはステップ222にジャンプする。これに対してAM
P≧A0 のときにはステップ221に進んでトルク変動
算出禁止フラグXNOCALがセット(XNOCAL←
“1”)される。即ち、図24の時刻t4 において行わ
れる割込みルーチンでは1番気筒#1のトルク変動量が
算出される。従ってこの割込みルーチンにおいてAMP
≧A0となり、トルク変動算出禁止フラグXNOCAL
がセットされると1番気筒#1のトルク変動量の算出、
即ち、DTa(i)が最小となる気筒のトルク変動量の
算出が禁止される。In step 216, the minimum elapsed time difference DT3
0 min is set as TMNREC. Then step 21
In step 7, the amplitude AMP of the fluctuation of DTa (i) is calculated by subtracting the minimum elapsed time difference TMNREC from the maximum elapsed time difference TMXREC. Next, at step 218, the initial value of the maximum elapsed time difference DT30max is set to DTa.
(I). Next, at step 219, the decrease flag XMNREC is reset (XMNREC ← “0”). Next, at step 220, the amplitude AMP is set to the set value A 0.
It is determined whether or not it is greater than At the time of the AMP <A 0 jumps to step 222. AM
When P ≧ A 0, the routine proceeds to step 221, and the torque fluctuation calculation prohibition flag XNOCAL is set (XNOCAL ←
"1"). That is, the torque variation of the first cylinder # 1 is in the interruption routine performed at the time t 4 in FIG. 24 is calculated. Therefore, in this interrupt routine, AMP
≧ A 0 and the torque fluctuation calculation prohibition flag XNOCAL
Is set, the torque fluctuation amount of the first cylinder # 1 is calculated,
That is, calculation of the torque fluctuation amount of the cylinder in which DTa (i) is minimized is prohibited.
【0091】ステップ222およびステップ223では
経過時間差DTa(i)が急変した気筒についてのトル
ク変動量の算出が禁止される。即ち、ステップ222で
は|DTa(i−2)−DTa(i−1)|がKo ・|
DTa(i−1)−DTa(i)|よりも大きいか否か
が判別される。ここで定数Ko は3.0から4.0程度
の値である。ステップ222において|DTa(i−
2)−DTa(i−1)|<Ko ・|DTa(i−1)
−DTa(i)|であると判別されたときには処理ルー
チンを完了し、|DTa(i−2)−DTa(i−1)
|≧Ko ・|DTa(i−1)−DTa(i)|である
と判別されたときにはステップ223に進んでトルク変
動算出禁止フラグXNOCALがセットされる。In steps 222 and 223, the calculation of the torque fluctuation amount for the cylinder in which the elapsed time difference DTa (i) has suddenly changed is prohibited. That is, at step 222 | DTa (i-2) -DTa (i-1) | is K o · |
DTa (i-1) -DTa (i) | is determined. Here, the constant Ko is a value of about 3.0 to 4.0. In step 222, | DTa (i−
2) -DTa (i-1) | <K o · | DTa (i-1)
-DTa (i) |, the processing routine is completed, and | DTa (i-2) -DTa (i-1).
| ≧ K o · | DTa ( i-1) -DTa (i) | a is judged to be the torque fluctuation calculation prohibition flag XNOCAL proceeds to step 223 when the are set.
【0092】即ち、今図24の時刻t3 における割込み
ルーチンであるとするとこのときには|DTa(4)−
DTa(2)|≧Ko ・|DTa(2)−DTa(1)
|であるか否かが判別される。図24に示されるように
DTa(4)に対してDTa(2)が急変すると|DT
a(4)−DTa(2)|はKo ・|DTa(2)−D
Ta(1)|よりも大きくなる。このときトルク変動算
出禁止フラグがセットされ、経過時間差DTa(i)が
急変した2番気筒#2のトルク変動量の算出が禁止され
る。[0092] That is, when an interrupt routine at time t 3 now Figure 24 at this time | DTa (4) -
DTa (2) | ≧ K o · | DTa (2) −DTa (1)
| Is determined. As shown in FIG. 24, when DTa (2) changes rapidly with respect to DTa (4), | DT
a (4) -DTa (2) | is K o · | DTa (2) -D
Ta (1) |. At this time, the torque fluctuation calculation prohibition flag is set, and the calculation of the torque fluctuation amount of the second cylinder # 2 in which the elapsed time difference DTa (i) has suddenly changed is prohibited.
【0093】次に図25に示すトルク変動算出ルーチン
について説明する。図25を参照すると、まず初めにス
テップ301において一つ前に燃焼が行われた(i−
1)番気筒のトルク変動量を算出すべきことを示すフラ
グXCAL(i−1)がセットされているか否かが判別
される。フラグXCAL(i−1)=“0”のとき、即
ちフラグXCAL(i−1)がセットされていないとき
には処理サイクルを完了する。これに対してフラグXC
AL(i−1)=“1”のとき、即ちフラグXCAL
(i−1)がセットされているときにはステップ302
に進んでフラグXCAL(i−1)がリセットされ、次
いでステップ303に進む。Next, the torque fluctuation calculation routine shown in FIG. 25 will be described. Referring to FIG. 25, first, in step 301, combustion was performed immediately before (i-
1) It is determined whether or not a flag XCAL (i-1) indicating that the torque variation of the cylinder No. should be calculated is set. When the flag XCAL (i-1) = "0", that is, when the flag XCAL (i-1) is not set, the processing cycle is completed. On the other hand, the flag XC
When AL (i-1) = "1", that is, when the flag XCAL
When (i-1) is set, step 302
The flag XCAL (i-1) is reset, and then the routine proceeds to step 303.
【0094】ステップ303では一つ前に燃焼が行われ
た気筒についてのトルク変動量の算出を禁止する禁止フ
ラグXNOCALがリセット(XNOCAL=“0”)
されているか否かが判別される。この禁止フラグがセッ
ト(XNOCAL=“1”)されているときにはステッ
プ311に進んで禁止フラグXNOCALがリセットさ
れる。これに対して禁止フラグがリセットされていると
きにはステップ304に進む。即ち、フラグXCALが
セットされており、かつ禁止フラグXNOCALがリセ
ットされているときのみステップ304に進む。ステッ
プ304では一つ前に燃焼が行われた(i−1)番気筒
の仮想の経過時間差Kb(i−1)(図12および図1
3参照)が次式に基づいて算出される。In step 303, the prohibition flag XNOCAL for prohibiting the calculation of the amount of torque fluctuation for the cylinder in which combustion was performed immediately before is reset (XNOCAL = "0").
It is determined whether or not it has been performed. When this prohibition flag is set (XNOCAL = "1"), the routine proceeds to step 311 where the prohibition flag XNOCAL is reset. On the other hand, when the prohibition flag has been reset, the routine proceeds to step 304. That is, the process proceeds to step 304 only when the flag XCAL is set and the prohibition flag XNOCAL is reset. In step 304, the virtual elapsed time difference Kb (i-1) of the (i-1) th cylinder in which combustion was performed immediately before (FIGS. 12 and 1)
3) is calculated based on the following equation.
【0095】Kb(i−1)={2・DTa(i−1)
+DTa(i)}/3 次いでステップ305では次式に基づいて一つ前に燃焼
が行われた(i−1)番気筒の実際の出力トルクDN
(i−1)が算出される。 DN(i−1)={30°/Tb(i−1)}2−{3
0°/Ta(i−1)}2 次いでステップ306では次式に基づいて一つ前に燃焼
が行われた(i−1)番気筒の仮想の出力トルクDNS
(i−1)が算出される。Kb (i−1) = {2 · DTa (i−1)
+ DTa (i)} / 3 Next, at step 305, based on the following equation, the actual output torque DN of the (i-1) th cylinder in which combustion was performed immediately before:
(I-1) is calculated. DN (i−1) = {30 ° / Tb (i−1)} 2 − {3
0 ° / Ta (i−1)} 2 Then, in step 306, the virtual output torque DNS of the (i−1) th cylinder in which combustion was performed immediately before based on the following equation:
(I-1) is calculated.
【0096】DNS(i−1)={30°/(TbO
(i−1)+Kb(i−1))}2−{30°/Ta
(i−1)}2 次いでステップ307では次式に示すように仮想の出力
トルクDNS(i−1)から実際の出力トルクDN(i
−1)を減算することによって一つ前に燃焼が行われた
(i−1)番気筒のトルク変動量DLN(i−1)が算
出される。DNS (i-1) = {30 ° / (TbO)
(I-1) + Kb (i-1))} 2- {30 ° / Ta
(I-1)} 2 Next, in step 307, the virtual output torque DNS (i-1) is converted to the actual output torque DN (i
By subtracting -1), the torque fluctuation amount DLN (i-1) of the (i-1) th cylinder in which combustion was performed immediately before is calculated.
【0097】 DLN(i−1)=DNS(i−1)−DN(i−1) 即ち、例えば今3番気筒#3のATDC30°であって
フラグXCAL(1)がセットされているとするとステ
ップ304において1番気筒#1に対する仮想の経過時
間差Kb(1)が算出され、次いでステップ305にお
いて1番気筒#1の実際の出力トルクDN(1)が算出
され、ステップ306において1番気筒#1の仮想の出
力トルクDNS(1)が算出され、次いでステップ30
7において1番気筒#1のトルク変動量DLN(1)が
算出される。DLN (i-1) = DNS (i-1) -DN (i-1) That is, for example, suppose that the ATDC30 ° of the third cylinder # 3 is set and the flag XCAL (1) is set. In step 304, a virtual elapsed time difference Kb (1) with respect to the first cylinder # 1 is calculated. Next, in step 305, the actual output torque DN (1) of the first cylinder # 1 is calculated. 1 virtual output torque DNS (1) is calculated, then step 30
At 7, the torque fluctuation amount DLN (1) of the first cylinder # 1 is calculated.
【0098】なお、各気筒の出力変動量GLN(i−
1)を求める場合にはステップ305からステップ30
7において以下のような計算を行えばよい。即ち、ステ
ップ305では実際のトルク変動量DN(i−1)を求
める代りに次式に基づいて実際の出力変動量GN(i−
1)を算出する。 GN(i−1)={30°/Tb(i−1)}−{30
°/Ta(i−1)} 次いでステップ306では仮想のトルク変動量DNS
(i−1)を求める代りに次式に基づいて仮想の出力変
動量GNS(i−1)を算出する。The output fluctuation amount GLN (i-
When determining 1), steps 305 to 30
7, the following calculation may be performed. That is, in step 305, instead of obtaining the actual torque fluctuation DN (i-1), the actual output fluctuation GN (i-
1) is calculated. GN (i-1) = {30 ° / Tb (i-1)}-{30
° / Ta (i-1)} Next, at step 306, the virtual torque fluctuation amount DNS
Instead of obtaining (i-1), the virtual output fluctuation amount GNS (i-1) is calculated based on the following equation.
【0099】GNS(i−1)={30°/(TbO
(i−1)+Kb(i−1))}−{30°/Ta(i
−1)} 次いでステップ307ではトルク変動量DLN(i−
1)を求める代りに次式に基づいて出力変動量GLN
(i−1)を算出する。 GLN(i−1)=GNS(i−1)−GN(i−1) さて、ステップ307においてトルク変動量DLN(i
−1)が算出されるとステップ308に進んでトルク変
動量DLN(i−1)が正であるか否かが判別される。
DLN(i−1)≧0であればステップ310にジャン
プして一つ前に燃焼が行われた気筒のトルク変動量DL
N(i−1)を積算すべきことを示す積算要求フラグX
DLNI(i−1)がセット(XDLNI(i−1)←
“1”)される。これに対してDLN(i−1)<0で
あればステップ309に進んでDLN(i−1)が零と
され、次いでステップ310に進む。なお、各気筒のト
ルクは上昇と低下を繰返し、この場合トルク変動量を求
めるにはトルクの上昇分かトルクの減少分のいずれかを
積算すればよい。図25に示すルーチンではトルクの減
少分のみを積算するようにしており、従って上述したよ
うにDLN(i−1)<0のときにはDLN(i−1)
を零にしている。GNS (i-1) = {30 ° / (TbO)
(I-1) + Kb (i-1))}-{30 ° / Ta (i
−1)} Next, at step 307, the torque fluctuation amount DLN (i−
Instead of obtaining 1), the output fluctuation amount GLN is calculated based on the following equation.
(I-1) is calculated. GLN (i-1) = GNS (i-1) -GN (i-1) In step 307, the torque fluctuation amount DLN (i
When -1) is calculated, the routine proceeds to step 308, where it is determined whether or not the torque fluctuation amount DLN (i-1) is positive.
If DLN (i−1) ≧ 0, the routine jumps to step 310 and the torque fluctuation amount DL of the cylinder in which combustion was performed immediately before is performed.
Integration request flag X indicating that N (i-1) should be integrated
DLNI (i-1) is set (XDLNI (i-1) ←
"1"). On the other hand, if DLN (i-1) <0, the routine proceeds to step 309, where DLN (i-1) is set to zero, and then proceeds to step 310. Note that the torque of each cylinder repeatedly increases and decreases. In this case, the amount of torque fluctuation may be obtained by integrating either the increase in torque or the decrease in torque. In the routine shown in FIG. 25, only the amount of decrease in torque is integrated. Therefore, when DLN (i-1) <0, as described above, DLN (i-1)
Is set to zero.
【0100】次に図26を参照しつつカウンタCDLN
IXの処理について説明する。このカウンタCDLNI
Xのカウント値は後に説明する機関のトルク変動値を算
出する際に使用される。図26を参照すると、まず初め
にステップ401において現在3番気筒#3のATDC
30°であるか否かが判別される。現在3番気筒#3の
ATDC30°でないときには処理サイクルを完了し、
現在3番気筒#3のATDC30°であるときにはステ
ップ402に進む。ステップ402では機関のトルク変
動値を算出するためのトルク変動値算出条件が成立して
いるか否かが判別される。例えば空燃比をリーンとする
条件が成立していないか、或いは機関負荷Q/Nの単位
時間当りの変化量ΔQ/Nが設定値以上であるか、或い
は機関回転数の単位時間当りの変化量ΔNが設定値以上
であるときにはトルク変動値算出条件が成立していない
と判断され、それ以外のときにはトルク変動値算出条件
が成立していると判断される。Next, the counter CDLN will be described with reference to FIG.
IX processing will be described. This counter CDLNI
The count value of X is used when calculating an engine torque fluctuation value described later. Referring to FIG. 26, first, at step 401, the ATDC of the current third cylinder # 3
It is determined whether it is 30 °. If the current ATDC of the third cylinder # 3 is not 30 °, the processing cycle is completed.
If it is currently ATDC 30 ° of the third cylinder # 3, the routine proceeds to step 402. In step 402, it is determined whether a torque fluctuation value calculation condition for calculating a torque fluctuation value of the engine is satisfied. For example, the condition for making the air-fuel ratio lean is not satisfied, or the variation ΔQ / N of the engine load Q / N per unit time is equal to or greater than a set value, or the variation of the engine speed per unit time. When ΔN is equal to or larger than the set value, it is determined that the torque fluctuation value calculation condition is not satisfied. Otherwise, it is determined that the torque fluctuation value calculation condition is satisfied.
【0101】ステップ402においてトルク変動値算出
条件が成立していると判断されたときにはステップ40
8に進んでカウント値CDLNIXが1だけインクリメ
ントされる。このカウント値CDLNIXのインクリメ
ント作用は3番気筒#3がATDC30°となる毎に、
即ち720°クランク角度毎に行われる。次いでステッ
プ409ではカウント値CDLNIXのインクリメント
作用が開始されてからカウント値CDLNIXがクリア
されるまでの間の機関回転数の平均値NAVE および吸入
空気量Qの平均値QAVE が算出される。If it is determined in step 402 that the torque fluctuation value calculation condition is satisfied, the process proceeds to step 40.
Proceeding to 8, the count value CDLNIX is incremented by one. The increment operation of the count value CDLNIX is performed every time the third cylinder # 3 reaches ATDC 30 °.
That is, it is performed every 720 ° crank angle. Next, at step 409, the average value N AVE of the engine speed and the average value Q AVE of the intake air amount Q from the start of the increment operation of the count value CDLNIX until the count value CDLNIX is cleared are calculated.
【0102】一方、ステップ402においてトルク変動
値算出条件が成立していないと判断されたときにはステ
ップ403に進んでカウント値CDLNIXがクリアさ
れる。次いでステップ404では各気筒に対するトルク
変動量DLN(i)の積算値DLNI(i)(この積算
値は後に説明するルーチンにおいて算出される)がクリ
アされ、次いでステップ405では各気筒に対する積算
カウント値CDLNI(i)(この積算カウント値は後
に説明するルーチンにおいて算出される)がクリアされ
る。On the other hand, when it is determined in step 402 that the torque fluctuation value calculation condition is not satisfied, the routine proceeds to step 403, where the count value CDLNIX is cleared. Next, at step 404, the integrated value DLNI (i) of the torque fluctuation amount DLN (i) for each cylinder (this integrated value is calculated in a routine described later) is cleared, and then at step 405, the integrated count value CDLNI for each cylinder. (I) (this integrated count value is calculated in a routine described later) is cleared.
【0103】次いでステップ406の目標トルク変動値
LVLLFBの算出ルーチンに進む。このルーチンは図
28に示されている。図28を参照すると、まず初めに
ステップ431ではロックアップをオフとすべき条件が
成立しているか否かが判断される。例えば自動変速機3
0の3速または4速において定常運転または緩加速運転
が行われているか或いは減速運転時にロックアップをオ
フとすべき条件が成立していないと判断され、それ以外
はロックアップをオフとすべき条件が成立していると判
断される。ロックアップをオフとすべき条件が成立して
いると判断されたときには次いでステップ432に進
み、スリップ量補正係数DLTFが零とされる。次いで
ステップ435に進む。これに対しロックアップをオフ
とすべき条件が成立していないと判断されたとき、即ち
ロックアップをオンとすべきか或いはスリップ制御を行
うべきときにはステップ433に進み、機関回転数Nと
タービンランナ37の回転数NTとの差(N−NT)と
してスリップ量NSLPが算出される。次いでステップ
434では図7に示すマップからスリップ量補正係数D
LTFが算出される。次いでステップ435に進む。Next, the routine proceeds to a routine for calculating the target torque fluctuation value LVLLFB in step 406. This routine is shown in FIG. Referring to FIG. 28, first, in step 431, it is determined whether a condition for turning off lock-up is satisfied. For example, automatic transmission 3
It is determined that steady-state operation or slow acceleration operation is being performed in the third or fourth speed of 0, or the condition for turning off lock-up is not satisfied during deceleration operation, and lock-up should be turned off otherwise. It is determined that the condition is satisfied. When it is determined that the condition for turning off the lockup is satisfied, the process proceeds to step 432, where the slip amount correction coefficient DLTF is set to zero. Next, the routine proceeds to step 435. On the other hand, when it is determined that the condition for turning off the lockup is not satisfied, that is, when the lockup should be turned on or the slip control should be performed, the process proceeds to step 433, where the engine speed N and the turbine runner 37 are determined. The slip amount NSLP is calculated as the difference (N-NT) from the rotational speed NT of the vehicle. Next, at step 434, the slip amount correction coefficient D is obtained from the map shown in FIG.
LTF is calculated. Next, the routine proceeds to step 435.
【0104】ステップ435では図6(B)に示すマッ
プから基本目標トルク変動量LVLLFBBが算出され
る。次いでステップ436では基本目標トルク変動量L
VLLFBBからスリップ量補正係数DLTFを減算す
ることにより目標トルク変動量LVLLFBが算出され
る。第1実施例では後に詳細に説明するように算出され
たトルク変動値がこの目標トルク変動値LVLLFBと
なるように空燃比がフィードバック制御される。In step 435, the basic target torque variation LVLLFBB is calculated from the map shown in FIG. Next, at step 436, the basic target torque variation L
The target torque variation LVLLFB is calculated by subtracting the slip amount correction coefficient DLTF from VLLFBB. In the first embodiment, the air-fuel ratio is feedback-controlled such that the torque fluctuation value calculated as described in detail later becomes the target torque fluctuation value LVLLFB.
【0105】次いでステップ407では各気筒のトルク
変動値DLNISM(i)(このトルク変動値は後に説
明するルーチンにおいて算出される)が図28のルーチ
ンから算出された目標トルク変動値LVLLFBとされ
る。図29は繰返し実行されるメインルーチンを示して
いる。このメインルーチンではまず初めにトルク変動値
の算出ルーチン(ステップ600)が実行される。この
ルーチンが図30および図31に示されている。次いで
リーンリミットフィードバック補正係数FLLFBの算
出ルーチン(ステップ700)が実行される。このルー
チンが図32に示されている。次いで予め定められたク
ランク角になったときに噴射時間算出ルーチン(ステッ
プ800)が実行される。このルーチンが図34に示さ
れている。次いでその他のルーチン(ステップ900)
が実行される。Next, at step 407, the torque fluctuation value DLNISM (i) of each cylinder (this torque fluctuation value is calculated in a routine to be described later) is set as the target torque fluctuation value LVLLFB calculated from the routine of FIG. FIG. 29 shows a main routine that is repeatedly executed. In this main routine, first, a torque fluctuation value calculation routine (step 600) is executed. This routine is shown in FIGS. Next, a calculation routine (step 700) of the lean limit feedback correction coefficient FLLFB is executed. This routine is shown in FIG. Next, when a predetermined crank angle is reached, an injection time calculation routine (step 800) is executed. This routine is shown in FIG. Then, other routines (step 900)
Is executed.
【0106】次に図30および図31に示されるトルク
変動値の算出ルーチンについて説明する。図30および
図31を参照すると、まず初めにステップ601におい
てトルク変動量DLN(i)を積算すべきことを示す積
算要求フラグXDLNI(i)がセット(XDLNI
(i)=“1”)されているか否かが判別される。積算
要求フラグXDLNI(i)がセットされていないとき
にはステップ609にジャンプし、積算要求フラグXD
LNI(i)がセットされているときにはステップ60
2に進む。ステップ602では積算要求フラグXDLN
I(i)がリセットされる。次いでステップ603では
トルク変動量DLN(i)がトルク変動量積算値DLN
I(i)に加算される。次いでステップ604では積算
カウント値CDLNI(i)が1だけインクリメントさ
れる。即ち、例えばステップ601において1番気筒に
ついての積算要求フラグXDLNI(1)がセットされ
たとするとステップ602においてこのフラグXDLN
I(1)がリセットされ、ステップ603においてトル
ク変動量積算値DLNI(1)が算出され、ステップ6
04において積算カウント値CDLNI(1)が1だけ
インクリメントされる。Next, a description will be given of a torque variation value calculation routine shown in FIGS. Referring to FIGS. 30 and 31, first, at step 601, an integration request flag XDLNI (i) indicating that the torque variation DLN (i) should be integrated is set (XDLNI).
(I) = “1”) is determined. If the integration request flag XDLNI (i) has not been set, the process jumps to step 609, where the integration request flag XD
Step 60 when LNI (i) is set
Proceed to 2. In step 602, the integration request flag XDLN
I (i) is reset. Next, at step 603, the torque variation DLN (i) is changed to the torque variation integrated value DLN.
It is added to I (i). Next, at step 604, the integrated count value CDLNI (i) is incremented by one. That is, for example, assuming that the integration request flag XDLNI (1) for the first cylinder is set in step 601, this flag XDLN is set in step 602.
I (1) is reset, and in step 603, the torque variation integrated value DLNI (1) is calculated.
At 04, the integrated count value CDLNI (1) is incremented by one.
【0107】次いでステップ605では積算カウント値
CDLNI(i)が“8”になったか否かが判別され
る。CDLNI(i)が“8”でないときにはステップ
609にジャンプし、CDLNI(i)が“8”になる
とステップ606に進んで次式から各気筒のトルク変動
値DLNISM(i)が算出される。 DLNISM(i)=DLNISM(i)+{DLNI
(i)−DLNISM(i)}/4 次いでステップ607では各気筒に対するトルク変動量
積算値DLNI(i)がクリアされ、次いでステップ6
08では積算カウント値CDLNI(i)がリセットさ
れる。Next, at step 605, it is determined whether or not the integrated count value CDLNI (i) has become "8". If CDLNI (i) is not "8", the routine jumps to step 609, and if CDLNI (i) becomes "8", the routine proceeds to step 606, where the torque fluctuation value DLNISM (i) of each cylinder is calculated from the following equation. DLNISM (i) = DLNISM (i) + {DLNI
(I) -DLNISM (i)} / 4 Next, in step 607, the torque variation integrated value DLNI (i) for each cylinder is cleared, and then step 6
At 08, the integrated count value CDLNI (i) is reset.
【0108】即ち、算出されたトルク変動量積算値DL
NI(i)とこれまで用いられてきたトルク変動値DL
NISM(i)との間に差があるときにはこの差{DL
NI(i)−DLNISM(i)}に1/4を乗算した
値がトルク変動値DLNISM(i)に加算される。従
って例えば1番気筒#1についての積算カウント値CD
LNI(1)が“8”になるとステップ606において
トルク変動値DLNISM(1)が算出されることにな
る。That is, the calculated torque variation integrated value DL
NI (i) and torque fluctuation value DL used so far
When there is a difference from NISM (i), this difference {DL
NI (i) -DLNISM (i)} multiplied by 1/4 is added to torque fluctuation value DLNISM (i). Therefore, for example, the integrated count value CD for the first cylinder # 1
When LNI (1) becomes “8”, the torque fluctuation value DLNISM (1) is calculated in step 606.
【0109】次いでステップ609では図26に示すル
ーチンにおいて算出されたカウント値CDLNIXが
“8”になったか否かが判別される。CDLNIXが
“8”でないときには処理サイクルを完了し、CDLN
IXが“8”になるとステップ610に進んで各気筒の
トルク変動値DLNISM(i)の平均値である平均ト
ルク変動値DLNISM(={DLNISM(1)+D
LNISM(2)+DLNISM(3)+DLNISM
(4)}/4)が算出される。次いでステップ611で
はカウント値CDLNIXがクリアされる。このように
して機関のトルク変動量を代表する値DLNISMが算
出される。Next, at step 609, it is determined whether or not the count value CDLNIX calculated in the routine shown in FIG. 26 has become "8". When CDLNIX is not "8", the processing cycle is completed and CDLNIX is not executed.
When IX becomes "8", the routine proceeds to step 610, where the average torque fluctuation value DLNISM (= {DLNISM (1) + D), which is the average of the torque fluctuation values DLNISM (i) of each cylinder.
LNISM (2) + DLNISM (3) + DLNISM
(4)} / 4) is calculated. Next, at step 611, the count value CDLNIX is cleared. In this way, the value DLNISM representing the torque fluctuation amount of the engine is calculated.
【0110】なお、前述したようにカウント値CDLN
IXは720°クランク角度毎に1だけインクリメント
され、いずれの気筒についてもトルクの算出が禁止され
たことがなければカウント値CDLNIXが“8”にな
ったときには全ての気筒に対する積算カウント値CDL
NI(1),CDLNI(2),CDLNI(3),C
DLNI(4)は既に“8”となっている。従ってこの
場合には全ての気筒についてトルク変動値DLNISM
(i)が算出される。一方、例えば1番気筒#1につい
てトルク変動量の算出が禁止されたとするとカウント値
CDLNIXが“8”になったときに1番気筒#1の積
算カウント値CDLNI(1)だけは“8”になってお
らず、斯くして1番気筒#1については新たなトルク変
動量積算値DLNI(1)は算出されていない。従って
この場合、ステップ610において平均トルク変動値D
LNISMを求める際には1番気筒#1だけについては
以前に算出されたトルク変動値DLNISM(1)が使
用される。As described above, the count value CDLN
IX is incremented by 1 every 720 ° crank angle, and unless the calculation of torque is prohibited for any cylinder, when the count value CDLNIX becomes “8”, the integrated count value CDL for all cylinders
NI (1), CDLNI (2), CDLNI (3), C
DLNI (4) is already "8". Therefore, in this case, the torque fluctuation value DLNISM for all cylinders
(I) is calculated. On the other hand, for example, if the calculation of the torque fluctuation amount is prohibited for the first cylinder # 1, when the count value CDLNIX becomes "8", only the integrated count value CDLNI (1) of the first cylinder # 1 becomes "8". As a result, a new integrated torque variation amount DLNI (1) is not calculated for the first cylinder # 1. Therefore, in this case, in step 610, the average torque fluctuation value D
When obtaining the LNISM, the torque fluctuation value DLNISM (1) previously calculated is used for only the first cylinder # 1.
【0111】次に図32を参照しつつFLLFB算出ル
ーチンについて説明する。図32を参照すると、まず初
めにステップ701においてリーンリミットフィードバ
ック補正係数FLLFBの更新条件が成立しているか否
かが判別される。例えば暖機運転時であるとき、或いは
機関の運転状態が図5において破線で囲まれた学習領域
にないときには更新条件が成立していないと判断され、
その他のときには更新条件が成立していると判断され
る。更新条件が成立していないときには処理サイクルを
完了し、更新条件が成立しているときにはステップ70
2に進む。Next, the FLLFB calculation routine will be described with reference to FIG. Referring to FIG. 32, first, at step 701, it is determined whether or not the condition for updating the lean limit feedback correction coefficient FLLFB is satisfied. For example, when the engine is warming up, or when the operating state of the engine is not in the learning area surrounded by the broken line in FIG. 5, it is determined that the update condition is not satisfied,
At other times, it is determined that the update condition is satisfied. When the update condition is not satisfied, the processing cycle is completed. When the update condition is satisfied, step 70 is executed.
Proceed to 2.
【0112】ステップ702では図28に示すルーチン
から目標トルク変動値LVLLFBが算出される。次い
でステップ703およびステップ704では目標トルク
変動値LVLLFBに応じた変動量判別値DH(n),
DL(n)に基づいて次式に示されるトルク変動レベル
LVLH(n),LVLL(n)が算出される。 LVLH(n)=LVLLFB+DH(n) LVLL(n)=LVLLFB+DL(n) ここで、変動量判別値DH(n)およびDL(n)は図
33(A)に示されるように予め定められている。即
ち、図33(A)からわかるようにDH(n)について
は3つの正の値が定められており、DH(3)>DH
(2)>DH(1)の関係を有する。更に、これらDH
(1),DH(2),DH(3)は目標トルク変動値L
VLLFBが大きくなるにつれて次第に増大する。一
方、DL(n)については3つの負の値が定められてお
り、DL(1)>DL(2)>DL(3)の関係を有す
る。更に、これらDL(1),DL(2),DL(3)
の絶対値は目標トルク変動値LVLLFBが大きくなる
につれて次第に増大する。In step 702, the target torque fluctuation value LVLLFB is calculated from the routine shown in FIG. Next, at step 703 and step 704, a fluctuation amount discrimination value DH (n) corresponding to the target torque fluctuation value LVLLFB,
Based on DL (n), torque fluctuation levels LVLH (n) and LVLL (n) represented by the following equations are calculated. LVLH (n) = LVLLFB + DH (n) LVLL (n) = LVLLFB + DL (n) Here, the fluctuation amount determination values DH (n) and DL (n) are predetermined as shown in FIG. . That is, as can be seen from FIG. 33A, three positive values are defined for DH (n), and DH (3)> DH
(2)> DH (1) Furthermore, these DH
(1), DH (2), DH (3) are the target torque fluctuation values L
It gradually increases as VLLFB increases. On the other hand, three negative values are defined for DL (n), and have a relationship of DL (1)> DL (2)> DL (3). Further, DL (1), DL (2), DL (3)
Of the target torque fluctuation value LVLLFB gradually increases.
【0113】ところで今、ステップ702において算出
された目標トルク変動値LVLLFBが破線で示される
値だったとする。この場合、ステップ703では破線上
のDH(1),DH(2),DH(3)を目標トルク変
動値LVLLFBに加算した値が夫々トルク変動レベル
LVLH(1),LVLH(2),LVLH(3)とさ
れ、ステップ704では破線上のDL(1),DL
(2),DL(3)を目標トルク変動値LVLLFBに
加算した値が夫々トルク変動レベルLVLL(1),L
VLL(2),LVLL(3)とされる。Now, suppose that the target torque fluctuation value LVLLFB calculated in step 702 is a value indicated by a broken line. In this case, in step 703, the values obtained by adding the DH (1), DH (2), and DH (3) on the broken line to the target torque fluctuation value LVLFB are the torque fluctuation levels LVLH (1), LVLH (2), LVLH ( 3), and in step 704, DL (1), DL on the broken line
The values obtained by adding (2) and DL (3) to the target torque fluctuation value LVLLFB are the torque fluctuation levels LVLL (1) and LLL, respectively.
VLL (2) and LVLL (3).
【0114】一方、図33(B)に示されるように各ト
ルク変動レベルLVLH(n),LVLL(n)間の領
域に対してフィードバック補正値+a1 ,+a2 ,+a
3 ,+a4 ,−b1 ,−b2 ,−b3 ,−b4 が予め定
められており、例えばトルク変動レベルがLVLH
(1)とLVLH(2)の間の領域に対してはフィード
バック補正値は+a2 となる。これらフィードバック補
正値は+a4 >+a3 >+a2 >+a1 でありかつ−b
1 >−b2 >−b3 >−b4 である。図33(B)に示
す各フィードバック補正値+a1 ,+a2 ,+a3 ,+
a4 ,−b1 ,−b 2 ,−b3 ,−b4 が図33(A)
の対応する領域に示されている。On the other hand, as shown in FIG.
Area between the torque fluctuation levels LVLH (n) and LVLL (n)
Feedback correction value + a1 , + ATwo , + A
Three , + AFour , -B1 , -BTwo , -BThree , -BFour Is predetermined
For example, when the torque fluctuation level is LVLH
Feed for the area between (1) and LVLH (2)
The back correction value is + aTwo Becomes These feedback supplements
Positive value is + aFour > + AThree > + ATwo > + A1 And -b
1 > -BTwo > -BThree > -BFour It is. As shown in FIG.
Each feedback correction value + a1 , + ATwo , + AThree , +
aFour , -B1 , -B Two , -BThree , -BFour Is shown in FIG.
In the corresponding area.
【0115】ステップ703およびステップ704にお
いて夫々トルク変動レベルLVLH(n),LVLL
(n)が算出されるとステップ705に進んで図30お
よび図31に示すトルク変動値の算出ルーチンにより算
出された平均トルク変動値DLNISMが図33(B)
に示されるどのトルク変動レベルLVLH(n),LV
LL(n)の間にあるかが判別される。次いでステップ
706では対応するフィードバック補正値DLFBが算
出される。例えば今、目標トルク変動値LVLLFBが
図33(A)において破線で示される値であり、算出さ
れた平均トルク変動値DLNISMが図33(B)のL
VLH(1)とLVLH(2)との間である場合、即ち
目標トルク変動値LVLLFBに対する平均トルク変動
値DLNISMの偏差が図33(A)の破線上において
DH(1)とDH(2)の間にある場合にはフィードバ
ック補正値DLFBは+a2 とされる。In steps 703 and 704, the torque fluctuation levels LVLH (n) and LVLL are respectively
When (n) is calculated, the routine proceeds to step 705, where the average torque fluctuation value DLNISM calculated by the torque fluctuation value calculation routine shown in FIG. 30 and FIG.
Which torque fluctuation level LVLH (n), LV
It is determined whether it is between LL (n). Next, at step 706, the corresponding feedback correction value DLFB is calculated. For example, now, the target torque fluctuation value LVLLFB is a value indicated by a broken line in FIG. 33 (A), and the calculated average torque fluctuation value DLNISM is represented by L in FIG. 33 (B).
When the difference between the average torque fluctuation value DLNISM and the target torque fluctuation value LVLLFB is between VLH (1) and LVLH (2), the deviation of DH (1) and DH (2) on the broken line in FIG. If it is between them, the feedback correction value DLFB is set to + a 2 .
【0116】次いでステップ707では図26に示すC
DLNIXの処理ルーチンのステップ409において求
められた機関回転数の平均値NAVE および吸入空気量Q
の平均値QAVE に基づいて更新すべきリーンリミットフ
ィードバック補正係数FLLFBijが図5に示されるど
の学習領域のリーンリミットフィードバック補正係数で
あるかが決定される。次いでステップ708ではステッ
プ707において決定されたリーンリミットフィードバ
ック補正係数FLLFBijにフィードバック補正値DL
FBが加算される。Next, at step 707, C shown in FIG.
Average engine speed N AVE and intake air amount Q obtained in step 409 of the processing routine of DLNIX
The lean limit feedback correction coefficient FLLFFB ij to be updated is determined based on the average value Q AVE of the learning region shown in FIG. Next, at step 708, a feedback correction value DL is added to the lean limit feedback correction coefficient FLLFB ij determined at step 707.
FB is added.
【0117】即ち、上述したように例えばDLNISM
>LVLLFBであって、LVLH(1)<DLNIS
M<LVLH(2)である場合にはリーンリミットフィ
ードバック補正係数FLLFBijに+a2 が加算され
る。その結果、空燃比が小さくなるので各気筒のトルク
変動量が減少せしめられる。一方、DLNISM<LV
LLFBであってLVLL(1)>DLNISM>LV
LL(2)である場合にはリーンリミットフィードバッ
ク補正係数FLLFBijに−b2 が加算される。その結
果、空燃比が大きくなるので各気筒のトルク変動量が増
大せしめられる。このようにして全気筒の平均トルク変
動値DLNISMが目標トルク変動値LVLLFBとな
るようにリーン運転時の空燃比が制御される。That is, as described above, for example, DLNISM
> LVLLFB and LVLH (1) <DLNIS
If M <LVLH (2), + a 2 is added to the lean limit feedback correction coefficient FLLFB ij . As a result, the air-fuel ratio becomes smaller, so that the torque fluctuation amount of each cylinder is reduced. On the other hand, DLNISM <LV
LLFB and LVLL (1)>DLNISM> LV
-B 2 is added to the lean limit feedback correction coefficient FLLFB ij if a LL (2). As a result, the air-fuel ratio increases, so that the amount of torque fluctuation in each cylinder is increased. In this way, the air-fuel ratio during the lean operation is controlled such that the average torque fluctuation value DLNISM of all cylinders becomes the target torque fluctuation value LVLLFB.
【0118】なお、図26に示すルーチンに示されるよ
うにトルク変動値の算出条件が成立しないときにはステ
ップ407において全ての気筒に対するDLNISM
(i)がLVLLFBとされ、斯くして平均トルク変動
値DLNISMも目標トルク変動値LVLLFBとされ
る。従ってこのときにはリーンリミットフィードバック
補正係数FLLFBijの更新は行われない。When the condition for calculating the torque fluctuation value is not satisfied as shown in the routine shown in FIG.
(I) is LVLLFB, and thus the average torque fluctuation value DLNISM is also the target torque fluctuation value LVLLFB. Therefore, at this time, the lean limit feedback correction coefficient FLLFB ij is not updated.
【0119】次に図34を参照しつつ燃料噴射時間の算
出ルーチンについて説明する。図34を参照すると、ま
ず初めにステップ801において図2に示すマップから
基本燃料噴射時間TPが算出される。次いでステップ8
02ではリーン運転を行うべき運転状態か否かが判別さ
れる。リーン運転を行うべき運転状態のときにはステッ
プ803に進んで理論空燃比フィードバック補正係数F
AFの値が1.0に固定される。次いでステップ804
では図4に示すマップからリーン補正係数FLEANが
算出され、次いでステップ805では図5に示すマップ
からリーンリミットフィードバック補正係数FLLFB
が読込まれる。次いでステップ809では次式に基づい
て燃料噴射時間TAUが算出される。Next, a routine for calculating the fuel injection time will be described with reference to FIG. Referring to FIG. 34, first, in step 801, the basic fuel injection time TP is calculated from the map shown in FIG. Then step 8
In 02, it is determined whether or not the operating state is such that the lean operation should be performed. When the vehicle is in the operating state in which the lean operation is to be performed, the routine proceeds to step 803, where the stoichiometric air-fuel ratio feedback correction coefficient F
The value of AF is fixed at 1.0. Then step 804
Then, the lean correction coefficient FLEAN is calculated from the map shown in FIG. 4, and then in step 805, the lean limit feedback correction coefficient FLLFFB is calculated from the map shown in FIG.
Is read. Next, at step 809, the fuel injection time TAU is calculated based on the following equation.
【0120】TAU=TP・FLEAN・FLLFB・
FAF+TAUV これに対し、ステップ802においてリーン運転を行う
べき運転状態でないと判別されたとき、即ち空燃比を理
論空燃比にすべきときにはステップ806に進んでリー
ン補正係数FLEANが1.0に固定され、次いでステ
ップ807においてリーンリミットフィードバック補正
係数FLLFBが1.0に固定される。次いでステップ
808では空燃比センサ17の出力信号に基づいて空燃
比が理論空燃比となるように理論空燃比フィードバック
補正係数FAFが制御される。次いでステップ809に
進み、燃料噴射時間TAUが算出される。TAU = TP • FLEAN • FLLFB •
FAF + TAUV On the other hand, when it is determined in step 802 that the operating state is not such that lean operation should be performed, that is, when the air-fuel ratio is to be set to the stoichiometric air-fuel ratio, the routine proceeds to step 806, where the lean correction coefficient FLEAN is fixed at 1.0. Next, at step 807, the lean limit feedback correction coefficient FLLFB is fixed to 1.0. Next, at step 808, the stoichiometric air-fuel ratio feedback correction coefficient FAF is controlled based on the output signal of the air-fuel ratio sensor 17 so that the air-fuel ratio becomes the stoichiometric air-fuel ratio. Next, the routine proceeds to step 809, where the fuel injection time TAU is calculated.
【0121】これまで述べてきた実施例ではクランクシ
ャフトの角速度に基づいて機関の出力変動またはトルク
変動を求めるようにしているが、機関の出力変動または
トルク変動をどのように求めてもよく、例えば燃焼室内
に配置した燃焼圧センサの出力信号に基づいて機関の出
力変動またはトルク変動を求めるようにすることもでき
る。また、機関駆動系またはクランクシャフト自体の捩
り振動が生じたとき、またはロータ14の外歯の間隔に
ばらつきがある場合に機関の出力変動またはトルク変動
をどのように補正してもよい。In the above-described embodiments, the output fluctuation or torque fluctuation of the engine is obtained based on the angular velocity of the crankshaft. However, the output fluctuation or torque fluctuation of the engine may be obtained in any manner. The output fluctuation or torque fluctuation of the engine may be determined based on the output signal of the combustion pressure sensor arranged in the combustion chamber. In addition, when the torsional vibration of the engine drive system or the crankshaft itself occurs, or when the interval between the external teeth of the rotor 14 varies, the output fluctuation or torque fluctuation of the engine may be corrected in any manner.
【0122】次に第2実施例を説明する。図3を参照し
て上述したように空燃比を空燃比制御領域内に制御でき
ると車両の良好な運転性を確保しつつ燃料消費率および
NOX の発生量を大幅に低減することができる。そこ
で、第2実施例では図3に示される空燃比制御領域内に
目標リーン空燃比を定め、空燃比がこの目標リーン空燃
比となるようにしている。Next, a second embodiment will be described. As described above with reference to FIG. 3, when the air-fuel ratio can be controlled within the air-fuel ratio control region, the fuel consumption rate and the generation amount of NO X can be significantly reduced while ensuring good driving performance of the vehicle. Therefore, in the second embodiment, a target lean air-fuel ratio is determined within the air-fuel ratio control region shown in FIG. 3 so that the air-fuel ratio becomes the target lean air-fuel ratio.
【0123】ところが、上述したようにトルク変動制御
領域はロックアップクラッチのスリップ量NSLPに応
じて変動するので空燃比制御領域もスリップ量NSLP
に応じて変動する。即ち、スリップ量NSLPが小さい
ときほどトルク変動制御領域を小さく定める必要がある
のでスリップ量NSLPが小さいときほど空燃比制御領
域をリッチ側に定める必要がある。そこで第2実施例で
はスリップ量NSLPが小さいときほど目標リーン空燃
比が小さくなるようにしている。However, as described above, the torque fluctuation control region fluctuates according to the slip amount NSLP of the lock-up clutch, so that the air-fuel ratio control region also has the slip amount NSLP.
It fluctuates according to. That is, the smaller the slip amount NSLP, the smaller the torque fluctuation control region needs to be set. Therefore, the smaller the slip amount NSLP, the more the air-fuel ratio control region needs to be set to the rich side. Therefore, in the second embodiment, the smaller the slip amount NSLP is, the smaller the target lean air-fuel ratio is.
【0124】第2実施例では燃料噴射時間TAUが次式
に基づいて算出される。 TAU=TP・FLNTG・FAFLN・FAF+TA
UV ここでFLNTGはリーン補正係数を、FAFLNはリ
ーン空燃比フィードバック補正係数を夫々示している。
基本燃料噴射時間TPは上述したように空燃比を理論空
燃比とするのに必要な噴射時間を示しており、図2に示
すマップの形で予めROM22内に記憶されている。In the second embodiment, the fuel injection time TAU is calculated based on the following equation. TAU = TP ・ FLNTG ・ FAFLN ・ FAF + TA
UV Here, FLNTG indicates a lean correction coefficient, and FAFLN indicates a lean air-fuel ratio feedback correction coefficient.
As described above, the basic fuel injection time TP indicates the injection time required for setting the air-fuel ratio to the stoichiometric air-fuel ratio, and is stored in the ROM 22 in advance in the form of a map shown in FIG.
【0125】リーン補正係数FLNTGは空燃比を目標
リーン空燃比とするための補正係数である。このリーン
補正係数FLNTGは目標リーン空燃比が図3の空燃比
制御領域の中央部となるように定められており、目標リ
ーン空燃比が大きくなるにつれて小さくなる。このリー
ン補正係数FLNTGは次式に基づいて算出される。 FLNTG=FLNTGB+DLAF ここでFLNTGBは基本リーン補正係数を、DLAF
はスリップ量補正係数を夫々示している。基本リーン補
正係数FLNTGBはロックアップがオフのときに最適
な予め実験により求められたリーン補正係数であり、図
35に示されるように機関負荷Q/Nおよび機関回転数
Nの関数としてマップの形で予めROM22内に記憶さ
れている。The lean correction coefficient FLNTG is a correction coefficient for setting the air-fuel ratio to the target lean air-fuel ratio. The lean correction coefficient FLNTG is determined such that the target lean air-fuel ratio is located at the center of the air-fuel ratio control region in FIG. 3, and decreases as the target lean air-fuel ratio increases. This lean correction coefficient FLNTG is calculated based on the following equation. FLNTG = FLNTGB + DLAF where FLNTGB represents the basic lean correction coefficient and DLAF
Indicates slip amount correction coefficients. The basic lean correction coefficient FLNTGB is an optimal lean correction coefficient obtained by an experiment in advance when the lock-up is OFF, and as shown in FIG. Is stored in the ROM 22 in advance.
【0126】スリップ量補正係数DLAFはスリップ量
NSLPに応じて基本リーン補正係数FLNTGBを補
正するものである。このスリップ量補正係数DLAFは
ロックアップがオフのときには零に維持され、ロックア
ップがオンまたはロックアップクラッチのスリップ制御
時には図36に示されるマップから算出される。即ち、
スリップ量補正係数DLAFは図36に示されるように
スリップ量NSLPが小さくなるほど大きくなる。な
お、図36においてNSLP=0はロックアップがオン
の場合を示している。スリップ量補正係数DLAFが大
きくなるにつれてリーン補正係数FLNTGが大きくな
るのでスリップ量NSLPが小さくなるにつれてリーン
補正係数FLNTGが大きくされることになり、従って
スリップ量NSLPが小さくなるにつれて目標リーン空
燃比が小さくされることになる。このようにするとスリ
ップ量NSLPに関わらず、即ちスリップ制御中である
か否かに関わらず目標リーン空燃比を最適にすることが
できる。なお、スリップ量補正係数DLAFはスリップ
量NSLPの関数として図36に示すマップの形で予め
ROM22内に記憶されている。The slip amount correction coefficient DLAF is for correcting the basic lean correction coefficient FLNTGB according to the slip amount NSLP. This slip amount correction coefficient DLAF is maintained at zero when lockup is off, and is calculated from the map shown in FIG. 36 when lockup is on or when the lockup clutch is slip controlled. That is,
As shown in FIG. 36, the slip amount correction coefficient DLAF increases as the slip amount NSLP decreases. In FIG. 36, NSLP = 0 indicates a case where lockup is on. Since the lean correction coefficient FLNTG increases as the slip amount correction coefficient DLAF increases, the lean correction coefficient FLNTG increases as the slip amount NSLP decreases. Therefore, the target lean air-fuel ratio decreases as the slip amount NSLP decreases. Will be done. In this manner, the target lean air-fuel ratio can be optimized regardless of the slip amount NSLP, that is, regardless of whether the slip control is being performed. The slip amount correction coefficient DLAF is stored in the ROM 22 in advance in the form of a map shown in FIG. 36 as a function of the slip amount NSLP.
【0127】リーン空燃比フィードバック補正係数FA
FLNは空燃比を目標リーン空燃比に維持するための係
数である。リーン空燃比フィードバック補正係数FAF
LNは空燃比を目標リーン空燃比に維持すべきときに空
燃比センサ17の出力信号に基づいて制御され、このと
きリーン空燃比フィードバック補正係数FAFLNはほ
ぼ1.0を中心として上下動する。Lean air-fuel ratio feedback correction coefficient FA
FLN is a coefficient for maintaining the air-fuel ratio at the target lean air-fuel ratio. Lean air-fuel ratio feedback correction coefficient FAF
LN is controlled based on the output signal of the air-fuel ratio sensor 17 when the air-fuel ratio is to be maintained at the target lean air-fuel ratio. At this time, the lean air-fuel ratio feedback correction coefficient FAFLN moves up and down about 1.0.
【0128】理論空燃比フィードバック補正係数FAF
は空燃比を理論空燃比に維持するための係数である。理
論空燃比フィードバック補正係数FAFは空燃比を理論
空燃比に維持すべきときに空燃比センサ17の出力信号
に基づいて制御され、このとき理論空燃比フィードバッ
ク補正係数FAFはほぼ1.0を中心として上下動す
る。Theoretical air-fuel ratio feedback correction coefficient FAF
Is a coefficient for maintaining the air-fuel ratio at the stoichiometric air-fuel ratio. The stoichiometric air-fuel ratio feedback correction coefficient FAF is controlled based on the output signal of the air-fuel ratio sensor 17 when the air-fuel ratio is to be maintained at the stoichiometric air-fuel ratio. Move up and down.
【0129】図35に示されるように破線により囲まれ
た運転領域内については機関の運転状態に応じてリーン
補正係数FLNTGが定められており、この運転領域内
では空燃比が目標リーン空燃比に維持される。これに対
して図35の破線で囲まれた領域外の運転領域では空燃
比が理論空燃比に維持される。空燃比を理論空燃比に維
持すべきときにはリーン補正係数FLNTGおよびリー
ン空燃比フィードバック補正係数FAFLNは1.0に
固定され、理論空燃比フィードバック補正係数FAFが
空燃比センサ17の出力信号に基づいて制御される。As shown in FIG. 35, a lean correction coefficient FLNTG is determined according to the operating state of the engine in an operating region surrounded by a broken line, and in this operating region, the air-fuel ratio becomes equal to the target lean air-fuel ratio. Will be maintained. On the other hand, the air-fuel ratio is maintained at the stoichiometric air-fuel ratio in the operation region outside the region surrounded by the broken line in FIG. When the air-fuel ratio should be maintained at the stoichiometric air-fuel ratio, the lean correction coefficient FLNTG and the lean air-fuel ratio feedback correction coefficient FAFLN are fixed to 1.0, and the stoichiometric air-fuel ratio feedback correction coefficient FAF is controlled based on the output signal of the air-fuel ratio sensor 17. Is done.
【0130】一方、空燃比を目標リーン空燃比に維持す
べきときには理論空燃比フィードバック補正係数FAF
が1.0に固定され、リーン空燃比フィードバック補正
係数FAFLNが空燃比センサ17の出力信号に基づい
て制御される。即ち、第2実施例において空燃比センサ
17は空燃比に対応した出力電圧を発生し、従って空燃
比センサ17の出力電圧から空燃比を検出することがで
きるようになっている。検出された空燃比が目標リーン
空燃比に対しリーンになるとリーン空燃比フィードバッ
ク補正係数FAFLNが増大せしめられ、検出された空
燃比が目標リーン空燃比に対しリッチになるとリーン空
燃比フィードバック補正係数FAFLNが減少せしめら
れ、斯くして空燃比が目標リーン空燃比に維持される。On the other hand, when the air-fuel ratio should be maintained at the target lean air-fuel ratio, the stoichiometric air-fuel ratio feedback correction coefficient FAF
Is fixed at 1.0, and the lean air-fuel ratio feedback correction coefficient FAFLN is controlled based on the output signal of the air-fuel ratio sensor 17. That is, in the second embodiment, the air-fuel ratio sensor 17 generates an output voltage corresponding to the air-fuel ratio, so that the air-fuel ratio can be detected from the output voltage of the air-fuel ratio sensor 17. When the detected air-fuel ratio becomes lean with respect to the target lean air-fuel ratio, the lean air-fuel ratio feedback correction coefficient FAFLN is increased. The air-fuel ratio is maintained at the target lean air-fuel ratio.
【0131】図37は第2実施例における燃料噴射時間
の算出ルーチンを示している。このルーチンはメインル
ーチン内で実行される。図37を参照すると、まず初め
にステップ1001において図2に示すマップから基本
燃料噴射時間TPが算出される。次いでステップ100
2ではリーン運転を行うべき運転状態か否かが判別され
る。リーン運転を行うべき運転状態のときにはステップ
1003に進んで理論空燃比フィードバック補正係数F
AFの値が1.0に固定される。次いでステップ100
4のリーン補正係数FLNTGの算出ルーチンに進む。
このルーチンが図38に示されている。FIG. 37 shows a routine for calculating the fuel injection time in the second embodiment. This routine is executed in the main routine. Referring to FIG. 37, first, in step 1001, the basic fuel injection time TP is calculated from the map shown in FIG. Then step 100
In 2, it is determined whether or not the vehicle is in an operating state in which the lean operation is to be performed. When the engine is in the operating state in which the lean operation is to be performed, the routine proceeds to step 1003, where the stoichiometric air-fuel ratio feedback correction coefficient F
The value of AF is fixed at 1.0. Then step 100
The routine proceeds to a routine for calculating a lean correction coefficient FLNTG of 4.
This routine is shown in FIG.
【0132】図38を参照すると、まず初めにステップ
1051ではロックアップをオフとすべき条件が成立し
ているか否かが判断される。上述したように例えば自動
変速機30の3速または4速において定常運転または緩
加速運転が行われているか或いは減速運転時にロックア
ップをオフとすべき条件が成立していないと判断され、
それ以外はロックアップをオフとすべき条件が成立して
いると判断される。ロックアップをオフとすべき条件が
成立していると判断されたときには次いでステップ10
52に進み、スリップ量補正係数DLAFが零とされ
る。次いでステップ1055に進む。これに対しロック
アップをオフとすべき条件が成立していないと判断され
たとき、即ちロックアップをオンとすべきか或いはスリ
ップ制御を行うべきときにはステップ1053に進み、
機関回転数Nとタービンランナ37の回転数NTとの差
(N−NT)としてスリップ量NSLPが算出される。
次いでステップ1054では図36に示すマップからス
リップ量補正係数DLAFが算出される。次いでステッ
プ1055に進む。Referring to FIG. 38, first, at step 1051, it is determined whether or not a condition for turning off lock-up is satisfied. As described above, for example, it is determined that the steady operation or the gradual acceleration operation is being performed at the third or fourth speed of the automatic transmission 30 or the condition for turning off the lockup during the deceleration operation is not satisfied,
Otherwise, it is determined that the condition for turning off the lockup is satisfied. If it is determined that the condition for turning off the lock-up is satisfied, then step 10 is executed.
Proceeding to 52, the slip amount correction coefficient DLAF is made zero. Next, the routine proceeds to step 1055. On the other hand, when it is determined that the condition for turning off the lockup is not satisfied, that is, when the lockup should be turned on or the slip control should be performed, the process proceeds to step 1053,
The slip amount NSLP is calculated as the difference (N-NT) between the engine speed N and the speed NT of the turbine runner 37.
Next, at step 1054, a slip amount correction coefficient DLAF is calculated from the map shown in FIG. Next, the routine proceeds to step 1055.
【0133】ステップ1055では図35に示すマップ
から基本リーン補正係数FLNTGBが算出される。次
いでステップ1056では基本リーン補正係数FLNT
GBにスリップ量補正係数DLAFを加算することによ
り目標リーン空燃比を表すリーン補正係数FLNTGが
算出される。次いでステップ1005に進み、空燃比セ
ンサ17の出力信号に基づいて空燃比が目標リーン空燃
比となるようにリーン空燃比フィードバック補正係数F
AFLNが制御される。次いでステップ1009では次
式に基づいて燃料噴射時間TAUが算出される。In step 1055, the basic lean correction coefficient FLNTGB is calculated from the map shown in FIG. Next, at step 1056, the basic lean correction coefficient FLNT
A lean correction coefficient FLNTG representing a target lean air-fuel ratio is calculated by adding the slip amount correction coefficient DLAF to GB. Next, the routine proceeds to step 1005, where the lean air-fuel ratio feedback correction coefficient F is set so that the air-fuel ratio becomes the target lean air-fuel ratio based on the output signal of the air-fuel ratio sensor 17.
AFLN is controlled. Next, at step 1009, the fuel injection time TAU is calculated based on the following equation.
【0134】TAU=TP・FLNTG・FAFLN・
FAF+TAUV これに対し、ステップ1002においてリーン運転を行
うべき運転状態でないと判別されたとき、即ち空燃比を
理論空燃比にすべきときにはステップ1006に進んで
リーン補正係数FLNTGが1.0に固定され、次いで
ステップ1007においてリーン空燃比フィードバック
補正係数FAFLNが1.0に固定される。次いでステ
ップ1008では空燃比センサ17の出力信号に基づい
て空燃比が理論空燃比となるように理論空燃比フィード
バック補正係数FAFが制御される。次いでステップ1
009に進み、燃料噴射時間TAUが算出される。TAU = TP / FLNTG / FAFLN /
FAF + TAUV On the other hand, when it is determined in step 1002 that the operating state is not such that lean operation should be performed, that is, when the air-fuel ratio should be set to the stoichiometric air-fuel ratio, the routine proceeds to step 1006, where the lean correction coefficient FLNTG is fixed at 1.0. Next, at step 1007, the lean air-fuel ratio feedback correction coefficient FAFLN is fixed at 1.0. Next, at step 1008, the stoichiometric air-fuel ratio feedback correction coefficient FAF is controlled based on the output signal of the air-fuel ratio sensor 17 so that the air-fuel ratio becomes the stoichiometric air-fuel ratio. Then step 1
Proceeding to 009, the fuel injection time TAU is calculated.
【0135】[0135]
【発明の効果】スリップ量に関わらず車両の運転性の悪
化を阻止しつつ燃料消費率およびNO X の発生量を低減
することができる。The driving performance of the vehicle is poor regardless of the slip amount.
Fuel consumption rate and NO while preventing XReduces the amount of
can do.
【図1】内燃機関の全体図である。FIG. 1 is an overall view of an internal combustion engine.
【図2】基本燃料噴射時間のマップを示す図である。FIG. 2 is a diagram showing a map of a basic fuel injection time.
【図3】NOx の発生量とトルク変動を示す図である。3 is a diagram showing a generation amount and the torque fluctuation of the NO x.
【図4】リーン補正係数FLEANのマップを示す図で
ある。FIG. 4 is a diagram showing a map of a lean correction coefficient FLEAN.
【図5】リーンリミットフィードバック補正係数のマッ
プを示す図である。FIG. 5 is a diagram showing a map of a lean limit feedback correction coefficient.
【図6】目標トルク変動値を示す図である。FIG. 6 is a diagram showing a target torque fluctuation value.
【図7】スリップ量補正係数DLTFを示す図である。FIG. 7 is a diagram showing a slip amount correction coefficient DLTF.
【図8】30°クランク角度の経過時間Ta(i),T
b(i)の変化を示すタイムチャートである。FIG. 8 shows elapsed times Ta (i) and T of a 30 ° crank angle.
It is a time chart which shows the change of b (i).
【図9】30°クランク角度の経過時間Ta(i)の変
化を示すタイムチャートである。FIG. 9 is a time chart showing a change in elapsed time Ta (i) at a 30 ° crank angle.
【図10】30°クランク角度の経過時間Ta(i),
Tb(i)の変化を示すタイムチャートである。FIG. 10 shows the elapsed time Ta (i) of a 30 ° crank angle,
It is a time chart which shows change of Tb (i).
【図11】30°クランク角度の経過時間Ta(i),
Tb(i)の変化を示すタイムチャートである。FIG. 11 shows the elapsed time Ta (i) of a 30 ° crank angle,
It is a time chart which shows change of Tb (i).
【図12】30°クランク角度の経過時間Ta(i)の
変化を示すタイムチャートである。FIG. 12 is a time chart showing a change in elapsed time Ta (i) at a 30 ° crank angle.
【図13】経過時間差DTa(i)とKb(i)との関
係を示す図である。FIG. 13 is a diagram showing a relationship between elapsed time differences DTa (i) and Kb (i).
【図14】減速運転時における経過時間Ta(i)の変
化を示すタイムチャートである。FIG. 14 is a time chart showing changes in elapsed time Ta (i) during deceleration operation.
【図15】経過時間差DTa(i)の変化を示すタイム
チャートである。FIG. 15 is a time chart showing changes in the elapsed time difference DTa (i).
【図16】割込みルーチンを示すフローチャートであ
る。FIG. 16 is a flowchart showing an interrupt routine.
【図17】経過時間差DTa(i)および経過時間Tb
(i)を算出するためのフローチャートである。FIG. 17 shows an elapsed time difference DTa (i) and an elapsed time Tb.
It is a flowchart for calculating (i).
【図18】経過時間差DTa(i)および経過時間Tb
(i)を算出するためのフローチャートである。FIG. 18 shows an elapsed time difference DTa (i) and an elapsed time Tb.
It is a flowchart for calculating (i).
【図19】KTa(i)を算出するためのフローチャー
トである。FIG. 19 is a flowchart for calculating KTa (i).
【図20】KTb(i)を算出するためのフローチャー
トである。FIG. 20 is a flowchart for calculating KTb (i).
【図21】トルク変動量算出の許可をチェックするため
のフローチャートである。FIG. 21 is a flowchart for checking permission of calculation of a torque fluctuation amount.
【図22】トルク変動量算出の許可をチェックするため
のフローチャートである。FIG. 22 is a flowchart for checking permission of calculation of a torque fluctuation amount.
【図23】トルク変動量算出の許可をチェックするため
のフローチャートである。FIG. 23 is a flowchart for checking permission for calculating a torque fluctuation amount.
【図24】経過時間差DTa(i)の変化とフラグXM
XREC,XMNRECの変化を示すタイムチャートで
ある。FIG. 24 shows changes in elapsed time difference DTa (i) and flag XM.
It is a time chart which shows a change of XREC and XMREC.
【図25】トルク変動量を算出するためのフローチャー
トである。FIG. 25 is a flowchart for calculating a torque fluctuation amount.
【図26】カウンタCDLNIXを処理するためのフロ
ーチャートである。FIG. 26 is a flowchart for processing a counter CDLNIX.
【図27】種々の値の計算タイミングを示す図である。FIG. 27 is a diagram showing calculation timings of various values.
【図28】目標トルク変動量を算出するためのフローチ
ャートである。FIG. 28 is a flowchart for calculating a target torque fluctuation amount.
【図29】メインルーチンを示すフローチャートであ
る。FIG. 29 is a flowchart showing a main routine.
【図30】トルク変動値を算出するためのフローチャー
トである。FIG. 30 is a flowchart for calculating a torque fluctuation value.
【図31】トルク変動値を算出するためのフローチャー
トである。FIG. 31 is a flowchart for calculating a torque fluctuation value.
【図32】リーンリミットフィードバック補正係数を算
出するためのフローチャートである。FIG. 32 is a flowchart for calculating a lean limit feedback correction coefficient.
【図33】変動量判別値DH(n),DL(n)および
トルク変動レベルLVLH(n),LVLL(n)を示
す図である。FIG. 33 is a diagram showing fluctuation amount determination values DH (n) and DL (n) and torque fluctuation levels LVLH (n) and LVLL (n).
【図34】燃料噴射時間を算出するためのフローチャー
トである。FIG. 34 is a flowchart for calculating a fuel injection time.
【図35】基本リーン補正係数FLNTGBのマップを
示す図である。FIG. 35 is a diagram showing a map of a basic lean correction coefficient FLNTGB.
【図36】スリップ量補正係数DLAFを示す図であ
る。FIG. 36 is a diagram showing a slip amount correction coefficient DLAF.
【図37】第2実施例において燃料噴射時間を算出する
ためのフローチャートである。FIG. 37 is a flowchart for calculating a fuel injection time in a second embodiment.
【図38】リーン補正係数FLNTGを算出するための
フローチャートである。FIG. 38 is a flowchart for calculating a lean correction coefficient FLNTG.
3…サージタンク 4…燃料噴射弁 8…スロットル弁 14…ロータ 15…クランク角センサ 32…トルクコンバータ 33…ロックアップ機構 3 ... Surge tank 4 ... Fuel injection valve 8 ... Throttle valve 14 ... Rotor 15 ... Crank angle sensor 32 ... Torque converter 33 ... Lock-up mechanism
Claims (2)
動が目標出力変動となるように空燃比を制御する内燃機
関の空燃比制御方法において、ロックアップクラッチの
スリップ量を検出し、該スリップ量に基づいて目標出力
変動を補正するようにした内燃機関の空燃比制御方法。1. An air-fuel ratio control method for an internal combustion engine that detects an output fluctuation of an internal combustion engine and controls an air-fuel ratio so that the output fluctuation becomes a target output fluctuation. An air-fuel ratio control method for an internal combustion engine which corrects a target output fluctuation based on a slip amount.
空燃比となるように空燃比を制御する内燃機関の空燃比
制御方法において、ロックアップクラッチのスリップ量
を検出し、該スリップ量に基づいて目標リーン空燃比を
補正するようにした内燃機関の空燃比制御方法。2. An air-fuel ratio control method for an internal combustion engine that detects an air-fuel ratio and controls the air-fuel ratio so that the air-fuel ratio becomes a target lean air-fuel ratio. An air-fuel ratio control method for an internal combustion engine, wherein a target lean air-fuel ratio is corrected on the basis of a target lean air-fuel ratio.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP24884397A JPH1182086A (en) | 1997-09-12 | 1997-09-12 | Method for controlling air-fuel ratio of internal combustion engine |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP24884397A JPH1182086A (en) | 1997-09-12 | 1997-09-12 | Method for controlling air-fuel ratio of internal combustion engine |
Publications (1)
Publication Number | Publication Date |
---|---|
JPH1182086A true JPH1182086A (en) | 1999-03-26 |
Family
ID=17184253
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP24884397A Pending JPH1182086A (en) | 1997-09-12 | 1997-09-12 | Method for controlling air-fuel ratio of internal combustion engine |
Country Status (1)
Country | Link |
---|---|
JP (1) | JPH1182086A (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100853048B1 (en) * | 2000-05-18 | 2008-08-19 | 로베르트 보쉬 게엠베하 | Device and method for controlling a drive unit of vehiche |
JP2009156293A (en) * | 2007-12-25 | 2009-07-16 | Toyota Motor Corp | Vehicular control device |
JP2010144686A (en) * | 2008-12-22 | 2010-07-01 | Honda Motor Co Ltd | Torque control device for internal combustion engine |
-
1997
- 1997-09-12 JP JP24884397A patent/JPH1182086A/en active Pending
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100853048B1 (en) * | 2000-05-18 | 2008-08-19 | 로베르트 보쉬 게엠베하 | Device and method for controlling a drive unit of vehiche |
JP2009156293A (en) * | 2007-12-25 | 2009-07-16 | Toyota Motor Corp | Vehicular control device |
JP2010144686A (en) * | 2008-12-22 | 2010-07-01 | Honda Motor Co Ltd | Torque control device for internal combustion engine |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3303739B2 (en) | Air-fuel ratio control method for internal combustion engine | |
US6199426B1 (en) | Method of detection of output fluctuation in internal combustion engine | |
JP3791032B2 (en) | Fuel injection control device for internal combustion engine | |
JP3303732B2 (en) | Control device for internal combustion engine | |
JPH1182086A (en) | Method for controlling air-fuel ratio of internal combustion engine | |
JP3216577B2 (en) | Air-fuel ratio control method for internal combustion engine | |
EP0811758B1 (en) | Method of controlling an air-fuel ratio of an engine | |
JP3262003B2 (en) | Output fluctuation detection method for multi-cylinder internal combustion engine | |
JPH09296743A (en) | Internal combustion engine control device for vehicle | |
JP3085220B2 (en) | Air-fuel ratio control method for internal combustion engine | |
JPH1182119A (en) | Output control device for internal combustion engine | |
EP0799983B1 (en) | Method of detection of angular velocity and torque in an internal combustion engine | |
JP3246325B2 (en) | Detection method in internal combustion engine | |
JP3303669B2 (en) | Air-fuel ratio control method for internal combustion engine | |
JP3246328B2 (en) | Detection method in internal combustion engine | |
JP3279179B2 (en) | Detection method in a multi-cylinder internal combustion engine | |
JPH1191409A (en) | Motive power transmission control device for internal combustion engine | |
JP3156588B2 (en) | Air-fuel ratio control method for internal combustion engine | |
JP3218970B2 (en) | Detection method in internal combustion engine | |
JPH10331685A (en) | Air-fuel ratio controller for internal combustion engine | |
JPH1130148A (en) | Method for controlling air-fuel ratio of internal combustion engine | |
JPH09281006A (en) | Change detection method in internal-combustion engine | |
JPH0719090A (en) | Stability controller of engine | |
JP2002180873A (en) | Driving force controller for vehicle | |
JP2005207366A (en) | Control device for internal combustion engine |