JPH1191409A - 内燃機関の動力伝達制御装置 - Google Patents

内燃機関の動力伝達制御装置

Info

Publication number
JPH1191409A
JPH1191409A JP9257029A JP25702997A JPH1191409A JP H1191409 A JPH1191409 A JP H1191409A JP 9257029 A JP9257029 A JP 9257029A JP 25702997 A JP25702997 A JP 25702997A JP H1191409 A JPH1191409 A JP H1191409A
Authority
JP
Japan
Prior art keywords
cylinder
value
elapsed time
calculated
correction coefficient
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
Application number
JP9257029A
Other languages
English (en)
Inventor
Nobuyuki Shibagaki
信之 柴垣
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Toyota Motor Corp
Original Assignee
Toyota Motor Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Toyota Motor Corp filed Critical Toyota Motor Corp
Priority to JP9257029A priority Critical patent/JPH1191409A/ja
Publication of JPH1191409A publication Critical patent/JPH1191409A/ja
Pending legal-status Critical Current

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)
  • Control Of Fluid Gearings (AREA)

Abstract

(57)【要約】 【課題】 車両の良好な運転性を維持しつつ燃焼消費率
およびNOX の発生量を十分に低減する。 【解決手段】 クランクシャフトの角速度に基づいて各
気筒のトルク変動量を検出し、検出されたトルク変動量
が目標トルク変動量となるように基準空燃比をリーンリ
ミットフィードバック補正係数FLLFBにより補正す
る。圧縮行程末期から爆発行程初期までのクランク角度
領域内に設定されたクランク角度範囲内におけるクラン
クシャフトの角速度を検出し、この角速度の変動に基づ
いて車両振動の代表値AMPを算出する。ロックアップ
機構がオンのときにFLLFBが第1の設定値F1より
も大きくかつAMPが第2の設定値A1よりも大きいと
きにはカウント値CSINPをインクリメントし、CS
INPが設定値CSINP1よりも大きくなったらロッ
クアップ機構がオンになるのを禁止する。

Description

【発明の詳細な説明】
【0001】
【発明の属する技術分野】本発明は内燃機関の動力伝達
制御装置に関する。
【0002】
【従来の技術】従来より、クランクシャフトの角速度に
基づいて各気筒の出力変動を検出し、検出された出力変
動に基づいて、リーン空燃比である基準空燃比をフィー
ドバック補正係数により補正するようにした内燃機関が
知られている。空燃比がリーンになるほど機関で発生す
るNOX 量が減少し、しかも燃料消費率が小さくなるの
で空燃比をできるだけリーンにするのが好ましい。とこ
ろが空燃比が或る程度以上リーンになると燃焼が不安定
となり、その結果トルク変動量が大きくなって車両振動
が増大し、斯くして車両の運転性が悪化する。そこでこ
の内燃機関では、排出NOX 量および燃料消費率をでき
るだけ低減し、同時に車両の運転性の悪化をできるだけ
低減しうるトルク変動量領域内に目標トルク変動量を定
め、トルク変動量が目標トルク変動量となるように基準
空燃比をフィードバック補正係数により補正するように
している。この場合、検出されたトルク変動量が目標ト
ルク変動量よりも小さいときには基準空燃比はリーン側
に補正され、検出されたトルク変動量が目標トルク変動
量よりも大きいときには基準空燃比はリッチ側に補正さ
れる。
【0003】ところがこのような内燃機関に、トルクコ
ンバータ内にロックアップ機構を設けた自動変速機を適
用すると次のような問題が生ずる。即ち、ロックアップ
機構がオンのときにはクランクシャフトとトルクコンバ
ータの入力軸とが直結状態とされるのでこのとき車両が
悪路を走行すると駆動輪に作用する外力が直接クランク
シャフトに加わることになる。その結果、燃焼が安定し
ているとしてもクランクシャフトの回転速度変動が大き
くなって基準空燃比がリッチ側に補正され、斯くして燃
料消費率およびNOX の発生量を十分に低減することが
できない。ところが、ロックアップ機構がオフとなると
車両が悪路を走行してもクランクシャフトの回転速度変
動はそれほど大きくならない。そこで、ロックアップ機
構がオンのときにフィードバック補正係数による基準空
燃比のリッチ側への補正量が設定時間以上、設定値より
も大きいときには車両が悪路を走行していると判断して
ロックアップ機構をオフにするようにした動力伝達制御
装置が公知である(特開平8−177550号公報参
照)。
【0004】
【発明が解決しようとする課題】ところで、上述したよ
うに燃焼圧が変動したときにもフィードバック補正係数
による基準空燃比の補正量はリッチ側に大きくなりうる
ので誤判定防止のために設定値を小さくすることはでき
ない。ところが、空燃比が急激に変動するのを阻止する
ためにフィードバック補正係数は比較的ゆっくりと変化
するので設定値を大きくすると、車両が実際に悪路を走
行してからフィードバック補正係数による基準空燃比の
リッチ側への補正量が設定時間以上、設定値よりも大き
くなるまでに長時間を要することになる。即ち、車両が
実際に悪路を走行し始めてから、車両が悪路を走行して
いると判断されるまでに長時間を要し、この間はロック
アップ機構がオンに維持されるので燃焼消費率およびN
X の発生量を十分に低減することができない。云い換
えると、フィードバック補正係数による基準空燃比のリ
ッチ側への補正量のみに基づいて車両が悪路を走行して
いるか否かを判断するようにすると速やかにかつ正確に
判断することができないという問題点がある。
【0005】
【課題を解決するための手段】上記問題点を解決するた
めに1番目の発明によれば、自動変速機を備え、自動変
速機のトルクコンバータ内にロックアップ機構を設け、
クランクシャフトの角速度に基づいて各気筒の出力変動
を検出し、出力変動に基づいて基準空燃比をフィードバ
ック補正係数により補正するようにした内燃機関におい
て、クランクシャフトの角速度に基づいて車両振動を代
表する代表値を算出し、ロックアップ機構がオンのとき
に設定時間以上、フィードバック補正係数による基準空
燃比のリッチ側への補正量が第1の設定値よりも大きく
かつ代表値が第2の設定値よりも大きいときにロックア
ップ機構がオンになるのを禁止するようにしている。ク
ランクシャフトの角速度変動に基づいて算出される車両
振動を代表する代表値は燃焼圧が変動したときにはそれ
ほど大きくならないが、車両が悪路を走行すると大きく
なる。そこで1番目の発明では、ロックアップ機構がオ
ンのときに設定時間以上、フィードバック補正係数によ
る基準空燃比のリッチ側への補正量が第1の設定値より
も大きくかつ代表値が第2の設定値よりも大きいときに
車両が悪路を走行していると判断してロックアップ機構
がオンになるのを禁止するようにしている。
【0006】2番目の発明では1番目の発明において、
圧縮行程末期から爆発行程初期までのクランク角度領域
内に設定されたクランク角度範囲内におけるクランクシ
ャフトの角速度を検出し、この角速度の変動に基づいて
前記車両振動の代表値を算出するようにしている。
【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
x を吸収し、空燃比が理論空燃比又はリッチになると
吸収したNOx を放出しかつ還元する機能を有する。
【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に入力される。
【0009】エアフローメータ6は吸入空気量に比例し
た出力電圧を発生し、この出力電圧が対応するAD変換
器28を介して入力ポート26に入力される。また、ス
ロットル弁8にはスロットル弁8がアイドリング開度に
あることを検出するためのアイドルスイッチ16が取付
けられ、このアイドルスイッチ16の出力信号が入力ポ
ート26に入力される。また、排気マニホルド9内には
空燃比を検出するための空燃比センサ(O2 センサ)1
7が配置されており、この空燃比センサ17の出力信号
が対応するAD変換器28を介して入力ポート26に入
力される。一方、出力ポート27は対応する駆動回路2
9を介して各燃料噴射弁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に伝達される。
【0011】一方、ロックアップ機構33は入力軸36
に対してその軸線方向に移動可能に取付けられかつ入力
軸36と共に回転するロックアップクラッチ板38を具
備する。通常は、即ちロックアップ機構33がオフのと
きには入力軸36内のオイル通路を介してロックアップ
クラッチ板38とポンプカバー34間の部屋39内に加
圧オイルが供給され、次いでこの部屋39から流出した
加圧オイルはポンプインペラ35およびタービンランナ
37の周りの部屋40内に送りこまれた後、入力軸36
内のオイル通路を介して排出される。このときロックア
ップクラッチ板38両側の部屋39,40間の圧力差は
ほとんど生じないためにロックアップクラッチ板38は
ポンプカバー34の内壁面から離れており、従ってこの
ときにはクランクシャフト12の回転力はポンプカバー
34、ポンプインペラ35およびタービンランナ37を
介して入力軸36に伝達される。なお、このときのロッ
クアップクラッチのスリップ量NSLP、即ちクランク
シャフトの回転数Nと自動変換器30の入力軸36即ち
タービンランナ37の回転数NTとの差(N−NT)は
かなり大きくなっている。
【0012】一方、ロックアップクラッチのスリップ制
御を行うべきときには入力軸36内のオイル通路を介し
て部屋40内に加圧オイルが供給され、部屋39内のオ
イルは入力軸36内のオイル通路を介して排出される。
このとき部屋40内の圧力が部屋39内の圧力よりも高
くなり、斯くしてロックアップクラッチ板38がポンプ
カバー34の内周面上に圧接されてクランクシャフトと
入力軸36とがスリップしつつ回転する状態となる。こ
のときポンプカバー34に対するロックアップクラッチ
板38の圧接力が強くなればなるほどスリップ量NSL
Pは小さくなり、従ってスリップ量NSLPは部屋3
9,40内に供給されるオイルを制御することによって
制御できる。なお、ロックアップクラッチのスリップ制
御時にはスリップ量NSLPは約50〜150r.p.m.に
制御される。
【0013】ロックアップクラッチのスリップ制御時に
比べて部屋40内に加圧オイルがさらに供給され、部屋
39からオイルがさらに排出されるとロックアップクラ
ッチ板38がポンプカバー34の内周面上に強く圧接さ
れてクランクシャフトと入力軸36とが同速度で回転す
る直結状態となり、斯くしてロックアップ機構33がオ
ンとなる。このときのスリップ量NSLPは零である。
部屋39,40内へのオイルの供給制御、即ちロックア
ップ機構33のオン・オフ制御およびロックアップクラ
ッチのスリップ制御は自動変速機30内に設けられた制
御弁によって制御され、この制御弁は電子制御ユニット
20の出力信号に基づいて制御される。また、自動変速
機20内には変速作用を行うための多数のクラッチが設
けられており、これらのクラッチも電子制御ユニット2
0の出力信号に基づいて制御される。
【0014】また、自動変速機30内には入力軸36、
即ちタービンランナ37の回転速度を表わす出力パルス
を発生する回転速度センサ41および出力軸31の回転
速度を表わす出力パルスを発生する回転速度センサ42
が配置されており、これら回転速度センサ41,42の
出力パルスは入力ポート26に入力される。図1に示す
内燃機関では燃料噴射時間TAUが次式に基づいて算出
される。
【0015】TAU=TP・FLEAN・FLLFB・
FAF+TAUV ここでTPは基本燃料噴射時間を、FLEANはリーン
補正係数を、FLLFBはリーンリミットフィードバッ
ク補正係数を、FAFは理論空燃比フィードバック補正
係数を、TAUVは無効噴射時間を夫々示している。基
本燃料噴射時間TPは空燃比を理論空燃比とするのに必
要な噴射時間を示している。この基本燃料噴射時間TP
は実験により求められ、この基本燃料噴射時間TPは機
関負荷Q/N(吸入空気量Q/機関回転数N)および機
関回転数Nの関数として図2に示すマップの形で予めR
OM22内に記憶されている。
【0016】リーン補正係数FLEANは空燃比をリー
ン空燃比とするための補正係数であり、このリーン補正
係数FLEANは機関負荷Q/Nおよび機関回転数Nの
関数として図4に示すマップの形で予めROM22内に
記憶されている。リーンリミットフィードバック補正係
数FLLFBは空燃比をリーン限界に維持するための補
正係数である。本実施例では吸入空気量Qと機関回転数
Nに対してリーン空燃比フィードバック制御に対する学
習領域が図5に示されるように例えば9つの領域で分け
られており、各学習領域に対して夫々リーンリミットフ
ィードバック補正係数FLLFB11〜FLLFB33が設
定されている。
【0017】理論空燃比フィードバック補正係数FAF
は空燃比を理論空燃比に維持するための係数である。理
論空燃比フィードバック補正係数FAFは空燃比を理論
空燃比に維持すべきときに空燃比センサ17の出力信号
に基づいて制御され、このとき理論空燃比フィードバッ
ク補正係数FAFはほぼ1.0を中心として上下動す
る。
【0018】図4に示されるように破線により囲まれた
運転領域内については機関の運転状態に応じてリーン補
正係数FLEANが定められており、この運転領域内で
は空燃比がリーン空燃比に維持される。これに対して図
4の破線で囲まれた領域外の運転領域では空燃比が理論
空燃比に維持される。空燃比を理論空燃比に維持すべき
ときにはリーン補正係数FLEANおよびリーンリミッ
トフィードバック補正係数FLLFBは1.0に固定さ
れ、理論空燃比フィードバック補正係数FAFが空燃比
センサ17の出力信号に基づいて制御される。
【0019】一方、空燃比をリーン空燃比に維持すべき
ときには理論空燃比フィードバック補正係数FAFが
1.0に固定され、即ち空燃比センサ17の出力信号に
基づくフィードバック制御が停止され、リーン補正係数
FLEANとリーンリミットフィードバック補正係数F
LLFBとにより空燃比がリーン空燃比に制御される。
なお、本実施例ではリーン補正係数により表わされる空
燃比を基準空燃比としている。
【0020】次に図3を参照しつつリーンリミットフィ
ードバック制御について説明する。図3は機関出力トル
ク変動量およびNOx 発生量と空燃比との関係を示して
いる。空燃比がリーンになるほど燃料消費率は小さくな
り、また空燃比がリーンになるほどNOx の発生量が少
なくなる。従ってこれらの点からみると空燃比はできる
だけリーンにすることが好ましいことになる。ところが
空燃比が或る程度以上リーンになると燃焼が不安定とな
り、その結果図3に示されるようにトルク変動量が大き
くなる。そこで本実施例では図3に示されるようにトル
ク変動が増大し始める空燃比制御領域内に空燃比を維持
するようにしている。
【0021】即ち具体的に云うとリーン補正係数FLE
ANはリーンリミットフィードバック補正係数FLLF
BをFLLFB=1.0としたときに空燃比が図3に示
される空燃比制御領域の中央部となるように定められて
いる。一方、リーンリミットフィードバック補正係数F
LLFBはトルク変動量に応じて図3に示されるトルク
変動制御領域内において制御され、トルク変動量が大き
くなればリーンリミットフィードバック補正係数FLL
FBが増大せしめられ、即ち空燃比が小さくされ、トル
ク変動量が小さくなればリーンリミットフィードバック
補正係数FLLFBが減少せしめられ、即ち空燃比が大
きくされる。このようにして空燃比が図3に示される空
燃比制御領域内に制御される。
【0022】なお、リーンリミットフィードバック補正
係数FLLFBはリーン補正係数FLEANが定められ
ている機関運転領域をカバーするように設定されてい
る。トルク変動量が図3に示されるトルク変動制御領域
内に制御されると良好な車両の運転性を確保しつつ燃料
消費率およびNOx の発生量を大巾に低減することがで
きる。ただし、このようにトルク変動量をトルク変動制
御領域内に制御するためにはトルク変動量を検出しなけ
ればならない。
【0023】ところでトルク変動量を算出する方法は従
来より種々の方法が提案されている。代表的な例を挙げ
ると燃焼室内に燃焼圧センサを取付けてこの燃焼圧セン
サの出力信号に基づきトルク変動量を算出する方法や、
或いは第1の角速度ωaの2乗と第2の角速度ωbの2
乗との差に基づいてトルク変動量を算出する方法が挙げ
られる。
【0024】燃焼圧センサを用いると燃焼圧センサを取
付けた気筒が発生するトルクを確実に検出することがで
きるという利点がある反面、燃焼圧センサが必要である
という欠点を有している。これに対して角速度ωa,ω
bは従来より内燃機関が備えているクランク角センサの
出力信号から算出することができるので角速度ωa,ω
bに基づき出力トルクを算出するようにした場合には新
たなセンサを設ける必要がないという利点がある。ただ
し、この場合機関駆動系が捩り振動を生ずるとトルク変
動量を正確に検出できなくなるという問題を有してい
る。しかしながらこの問題を解決しさえすれば新たなセ
ンサを必要としない角速度に基づくトルク算出方法の方
が好ましいことは明らかである。そこで本実施例では発
生トルクを角速度に基づき算出するようにし、その際機
関駆動系が捩り振動を生じたとしてもトルク変動量を正
確に検出しうるようにしている。
【0025】次に機関の出力変動およびトルク変動を算
出するための方法について説明する。まず初めに、機関
駆動系が捩り振動を生じていない定常運転時を示す図6
(A),(B)を参照しつつ各気筒が発生する駆動力お
よび各気筒が発生するトルクを算出する方法について説
明する。前述したようにクランク角センサ15はクラン
クシャフトが30°クランク角度回転する毎に出力パル
スを発生し、更にクランク角センサ15は各気筒#1,
#2,#3,#4の圧縮上死点TDCにおいて出力パル
スを発生するように配置されている。従ってクランク角
センサ15は各気筒#1,#2,#3,#4の圧縮上死
点TDCから30°クランク角毎に出力パルスを発生す
ることになる。なお、本実施例において用いられている
内燃機関の点火順序は1−3−4−2である。
【0026】図6(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で除算するとこの除算結果は角速度ωを表わしてい
る。本発明による実施例では30°クランク角度/Ta
(i)をi番気筒における第1の角速度ωaと称し、3
0°クランク角度/Tb(i)をi番気筒における第2
の角速度ωbと称する。従って30°クランク角度/T
a(1)は1番気筒の第1の角速度ωaを表わし、30
°クランク角度/Tb(1)は1番気筒の第2の角速度
ωbを表わすことになる。
【0027】図6(A),(B)の1番気筒に注目して
みると、燃焼が開始されて燃焼圧が高まると経過時間が
Ta(1)からTb(1)まで低下し、次いでTb
(1)から再び上昇する。云い換えるとクランクシャフ
トの角速度ωが第1の角速度ωaから第2の角速度ωb
まで上昇し、次いで第2の角速度ωbから再び下降す
る。即ち、燃焼圧によってクランクシャフトの角速度ω
が第1の角速度ωaから第2の角速度ωbへと増大せし
められたことになる。図6(A)は燃焼圧が比較的高い
場合を示しており、図6(B)は燃焼圧が比較的低い場
合を示している。図6(A),(B)から燃焼圧が高い
場合には燃焼圧が低い場合に比べて経過時間の減少量
(Ta(i)−Tb(i))が大きくなり、従って角速
度ωの増大量(ωb−ωa)が大きくなる。燃焼圧が高
くなればその気筒の発生する駆動力が大きくなり、従っ
て角速度ωの増大量(ωb−ωa)が大きくなれば気筒
の発生する駆動力が大きくなることになる。従って第1
の角速度ωaと第2の角速度ωbとの差(ωb−ωa)
から気筒の発生する駆動力を算出することができる。
【0028】一方、機関の回転慣性モーメントをIとす
ると燃焼圧によって運動エネルギが(1/2)Iωa2
から(1/2)Iωb2 に増大せしめられる。この運動
エネルギの増大量(1/2)・I・(ωb2 −ωa2
はその気筒が発生するトルクを表わしており、従って第
1の角速度ωaの2乗と第2の角速度ωbの2乗との差
(ωb2 −ωa2 )から気筒の発生するトルクを算出で
きることになる。
【0029】このように第1の角速度ωaと第2の角速
度ωbを検出すればこれらの検出値から対応する気筒の
発生する駆動力および対応する気筒の発生するトルクを
算出できることになる。なお、図6(A),(B)に示
される経過時間T30の変化は機関によって若干異な
り、従って第1の角速度ωaを検出すべきクランク角度
範囲および第2の角速度ωbを検出すべきクランク角度
範囲は機関に応じて(ωb−ωa)が機関の発生する駆
動力を最もよく表わすように、或いは(ωb2 −ωa
2 )が機関の発生するトルクを最もよく表わすように定
められる。従って機関によっては第1の角速度ωaを検
出すべきクランク角度範囲が圧縮上死点前BTDC30
°からTDCであり、第2の角速度ωbを検出すべきク
ランク角度範囲がATDC90°からATDC120°
となることもあり得る。
【0030】従って各角速度ωa,ωbの検出のしかた
について一般的に表現すると、圧縮行程末期から爆発行
程初期までのクランク角度領域内に第1のクランク角度
範囲を設定し、第1のクランク角度範囲から一定のクラ
ンク角を隔てた爆発行程中期のクランク角度領域内に第
2のクランク角度範囲を設定し、第1のクランク角度範
囲内におけるクランクシャフトの第1の角速度ωaを検
出し、第2のクランク角度範囲内におけるクランクシャ
フトの第2の角速度ωbを検出するということになる。
【0031】上述したように角速度ωa,ωbを検出す
れば検出値に基づいて対応する気筒の発生する駆動力お
よびトルクを算出することができる。しかしながら機関
駆動系には各気筒において順次行われる爆発作用により
駆動系の固有振動数でもって振動する捩り振動が発生し
ており、このように機関駆動系に捩り振動が発生してい
ると角速度ωa,ωbに基づいて気筒の発生する駆動力
およびトルクを正確に算出することができなくなる。次
にこのことについて図7および図8を参照しつつ説明す
る。
【0032】図7は機関駆動系に捩り振動が発生してい
るときに各気筒に対し順次算出される経過時間Ta
(i)の変化を示しており、図7のAMQは最小のTa
(i)と最大のTa(i)との差、即ち振幅を示してい
る。機関駆動系に捩り振動が発生するとこの捩り振動に
よってクランクシャフトの角速度が周期的に増大減少せ
しめられるので経過時間Ta(i)は図7に示されるよ
うに周期的に増大減少することになる。
【0033】一方、図8は図7において経過時間Ta
(i)が減少している部分を拡大して示している。図8
に示されるように経過時間Ta(i)はTa(1)とT
a(3)との間でho時間だけ減少しており、このho
時間の減少は捩り振動による捩れ量の増大によるものと
考えられる。この場合、Ta(1)とTa(3)との間
では捩り振動による経過時間の減少量は時間の経過と共
にほぼ直線的に増大するものと考えられ、従ってこの捩
り振動による経過時間の減少量はTa(1)およびTa
(3)を結ぶ破線とTa(1)を通る水平線との差で表
わされることになる。従ってTa(1)とTb(1)と
の間では捩り振動によって経過時間がhだけ減少してい
ることになる。
【0034】このように機関駆動系に捩り振動が発生す
るとTb(1)はTa(1)に対して経過時間が減少
し、この減少した経過時間は燃焼圧による経過時間の減
少量fと捩り振動による経過時間の減少量hとを含んで
いることになる。従って燃焼圧により減少した経過時間
Tb(1)だけを求めるためにはTb(1)にhを加算
しなければならないことになる。即ち、検出された経過
時間Ta(i)およびTb(i)に基づいて各気筒が発
生する駆動力或いはトルクを求めても真の駆動力或いは
トルクを求めることができず、斯くして真の機関の出力
変動或いはトルク変動を求めることができない。
【0035】更に、多気筒内燃機関ではこのような機関
駆動系の捩り振動に加えてクランクシャフト自体の捩り
振動を発生し、このようなクランクシャフト自体の捩り
振動が発生した場合にも真の機関の出力変動或いはトル
ク変動を求めることができなくなる。次にこのことにつ
いて図9を参照しつつ説明する。多気筒内燃機関、例え
ば図1に示されるような4気筒内燃機関では1番気筒お
よび2番気筒において大きなクランクシャフト自体の捩
り振動が発生する。即ち、クランクシャフト自体の捩り
振動がほとんど発生しない気筒、例えば3番気筒#3で
は図9に示されるようにTa(3)からTb(3)に向
けて経過時間は徐々に減少するが1番気筒#1において
はTa(1)からTb(1)に向けて経過時間は徐々に
減少せず、クランクシャフト自体の捩り振動によってT
b(1)の経過時間が長くなってしまう。その結果、1
番気筒#1については検出された経過時間Ta(1)お
よびTb(1)に基づいて1番気筒#1が発生する駆動
力或いはトルクを求めても真の駆動力或いはトルクを求
めることができず、斯くして真の機関の出力変動或いは
トルク変動を求めることができない。
【0036】そこで本実施例では機関駆動系の捩り振動
が発生しても、またクランクシャフト自体の捩り振動が
発生しても真の機関の出力変動或いはトルク変動を算出
することができる算出方法を採用している。次にこの算
出方法について図10を参照しつつ説明する。図10に
おいて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の次の燃焼時における
経過時間を表わしている。
【0037】まず初めに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)を一直線上における高さとして書き直すと
図11に示されるようになる。
【0038】図11に示されるように経過時間差DTa
(1)とDTa(3)との間では経過時間差がI0 時間
だけ減少している。クランクシャフト自体に捩り振動が
発生したとしてもこの捩り振動による影響は経過時間差
には表れず、従って経過時間差のI0 時間の減少は機関
駆動系の捩り振動によるものである。この場合、DTa
(1)とDTa(3)との間では機関駆動系の捩り振動
による経過時間差の減少量は時間の経過と共にほぼ直線
的に増大するものと考えられる。従って1番気筒#1の
先の燃焼時の出力トルクと次の燃焼時の出力トルクとが
同一であると仮定すると経過時間差DTa(1)とKb
(1)との間では機関駆動系の捩り振動によって経過時
間差がIだけ減少するものと考えられる。従って図11
からわかるように1番気筒#1の先の燃焼時の出力トル
クと次の燃焼時の出力トルクとが同一であると仮定した
場合に経過時間差Kb(1)は次式で表される。
【0039】 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)は次式で表される。
【0040】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)は
次式で表すことができる。
【0041】DLN(i)=DNS(i)−DN(i)
={30°/Tb′(i)}2−{30°/Tb
(i)}2 即ち、i番気筒のトルク変動量はi番気筒の仮想の第2
の角速度の2乗と実際の第2の角速度の2乗との偏差か
ら求まることになる。このような方法でもってトルク変
動量を算出することによって機関駆動系に捩り振動が発
生したとしても、更にクランクシャフト自体に捩り振動
が発生したとしても各気筒のトルク変動量を正確に検出
できることになる。
【0042】一方、このような方法でもってトルク変動
量を算出する場合、ロータ14(図1)の外周に沿って
形成されている外歯の間隔にばらつきがあると各気筒の
トルク変動量を正確に検出することができない。そこで
本実施例ではロータ14の外歯の間隔にばらつきがあっ
たとしても各気筒のトルク変動量を正確に検出しうるよ
うに検出された各経過時間Ta(i),Tb(i)を補
正するようにしている。次に経過時間Ta(i)の補正
方法について図12を参照しつつ説明する。
【0043】本実施例では経過時間Ta(i)が次式に
基づいて算出される。 Ta(i)=(i番気筒のTDCからATDC30°ま
での所要時間)・(1+KTa(i)) ここでKTa(i)はi番気筒に対する補正係数を表し
ており、これらの補正係数KTa(i)は機関駆動系の
捩り振動が発生せず、クランクシャフト自体の捩り振動
も発生しない減速運転時の燃料供給停止時に算出され
る。
【0044】図12は減速運転時における経過時間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)は図12の実線に沿っ
て変化するものと考えられ、検出された経過時間Ta
(i)が図12の実線からずれている場合には図12の
実線からずれている経過時間Ta(i)を検出するため
の外歯の間隔が正規の間隔からずれているものと考えら
れる。そこで本実施例では経過時間Ta(i)が図12
に示す直線からずれている場合にはその経過時間Ta
(i)が図12の実線上に位置するように経過時間Ta
(i)を補正係数KTa(i)によって補正するように
している。
【0045】具体的に云うと、本実施例では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)を次式より求める。
【0046】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)を求める。
【0047】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)を各気筒
毎に求める。
【0048】KTaB(i)=KTa(i)j-1 −KA
FC(i)/4 即ち、各気筒に対するこれまでの補正係数KTa(i)
j-1 からずれ割合KAFC(i)の1/4を減算するこ
とによって基本補正係数KTaB(i)が算出される。
次いで基本補正係数KTaB(i)の平均値tKTaM
が次式に基づき算出される。
【0049】tKTaM={KTaB(1)+KTaB
(2)+KTaB(3)+KTaB(4)}/4 次いで次式に示されるように各基本補正係数KTaB
(i)から平均値tKTaMを減算することによって各
気筒に対する補正係数KTa(i)が算出される。
【0050】 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)がずれ
を生じることになる。
【0051】このような問題が生じないようにするため
に基本補正係数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)はずれ
を生じなくなる。
【0052】図12に示される実施例では例えば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の外歯の間隔にばらつきがあったとしても各気筒の
トルク変動を正確に検出することができる。
【0053】再びトルク変動の算出に話しを戻すと機関
駆動系に捩り振動が発生したときに経過時間差DTa
(i)は図10に示されるように変動する。しかしなが
らこの経過時間差DTa(i)は車両が凸凹道を走行し
たときにも変動し、しかもこのときにはDTa(i)の
変動巾が極めて大きくなる場合がある。図13は車両が
凸凹道を走行したときのDTa(i)の変動を示してお
り、図13のAMPは最小のDTa(i)と最大のDT
a(i)との差、即ち振幅を示している。この振幅AM
Pが小さいときにはこれまで述べた方法によって各気筒
の出力変動およびトルク変動を正確に検出することがで
きる。
【0054】しかしながら振幅AMPが大きくなるとD
Ta(i)が最大又は最小となる気筒の出力変動又はト
ルク変動を正確に検出できなくなる。即ち、図13にお
いて例えば最初にDTa(i)が最大になる気筒が1番
気筒#1であったとすると1番気筒#1の仮想の経過時
間差Kb(1)の捩り振動による減少量Iは図11の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)が最小とな
る気筒においても同じことが生ずる。
【0055】また、一つ前に燃焼が行われた気筒のDT
a(i)に対してDTa(i)が急変した気筒において
もKb(i)の値が実際の値からずれ、斯くして出力変
動およびトルク変動を正確に検出できなくなる。そこで
本実施例では振幅AMPが大きいときにはDTa(i)
が最大又は最小となる気筒については出力変動量又はト
ルク変動量を求めずに、更に一つ前に燃焼が行われた気
筒のDTa(i)に対してDTa(i)が急変した気筒
についても出力変動量又はトルク変動量を求めないよう
にしている。
【0056】さらに本実施例では、クランクシャフトの
回転速度変動に基づいて車両振動を代表する代表値を求
め、ロックアップ機構33がオンのときにリーンリミッ
トフィードバック補正係数FLLFBと、この代表値と
に基づいて車両が悪路を走行しているか否かを判別し、
車両が悪路を走行していると判別されたときにはロック
アップ機構33がオンになるのを禁止するようにしてい
る。
【0057】即ち、冒頭で述べたようにロックアップ機
構33がオンのときに車両が悪路を走行すると駆動輪に
作用する外力が直接クランクシャフトに伝達されるので
このとき算出されるトルク変動量が大きくなり、その結
果上述したようにリーンリミットフィードバック補正係
数FLLFBが大きくなる。ところが、リーンリミット
フィードバック補正係数FLLFBが大きくなって空燃
比がリッチになっても車両が悪路を走行し続ける限り、
算出されるトルク変動量は小さくならず、リーンリミッ
トフィードバック補正係数FLLFBは増大し続ける。
即ち、車両が悪路を走行するとリーンリミットフィード
バック補正係数FLLFBが一定時間以上、一定値より
も大きくなる。
【0058】一方、車両が悪路を走行すると車両の振動
が大きくなり、ロックアップ機構33がオンのときには
この車両振動がクランクシャフトに伝達され、その結果
上述した経過時間差DTa(i)の振幅AMPが大きく
なる。ところがこの振幅AMPは燃焼圧が変動したとき
にも増大し、しかしながら燃焼圧が変動したときには振
幅AMPは車両が悪路を走行したときよりも大きくなら
ない。即ち、車両が悪路を走行すると車両振動を代表す
る振幅AMPが一定時間以上、一定値よりも大きくな
る。
【0059】そこで本実施例では、設定時間以上、リー
ンリミットフィードバック補正係数FLLFBが第1の
設定値よりも大きくかつ振幅AMPが第2の設定値より
も大きいときには車両が悪路を走行していると判断する
ようにしている。ロックアップ機構33をオフにすると
駆動輪に作用する外力はポンプインペラを介してクラン
クシャフトに伝達されるので車両が悪路を走行したとし
てもクランクシャフトの回転速度変動は小さくなる。そ
こで、ロックアップ機構33がオンになるのを禁止すべ
きときには常にロックアップ機構33をオフにするよう
にすることもできる。しかしながら、ロックアップ機構
33をオフにすると駆動力の伝達ロスが大きくなるので
好ましくない。一方、ロックアップクラッチのスリップ
制御を行うとクランクシャフトに伝達される外力を小さ
く維持することができ、同時に駆動力の伝達ロスも小さ
く維持することができる。そこで本実施例では、ロック
アップ機構33がオンのときに車両が悪路を走行してい
ると判断されたときにはロックアップクラッチのスリッ
プ制御が可能であるか否かを判断し、ロックアップクラ
ッチのスリップ制御が可能であると判断されたときには
スリップ制御を行い、可能でないと判断されたときには
ロックアップ機構33をオフにするようにしている。
【0060】図14はロックアップ機構33をオンにす
べき運転領域およびロックアップクラッチのスリップ制
御が可能な領域を示している。即ち、図14に示すアク
セルペダルの踏み込み量DEPと車速VVとにより定ま
る運転領域において、領域SOおよび領域SSOはロッ
クアップ機構33をオンにすべき運転領域を示してお
り、それ以外の領域はロックアップ機構33をオフにす
べき運転領域を示している。また、図14の領域SSお
よび領域SSOはロックアップクラッチのスリップ制御
が可能な運転領域を示しており、それ以外の領域はロッ
クアップクラッチのスリップ制御が可能でない運転領域
を示している。
【0061】通常時、即ち車両が平滑路を走行している
と判断されたときにDEPおよびVVにより定まる点が
領域SO内または領域SSO内にあるときにはロックア
ップ機構33がオンにされ、DEPおよびVVにより定
まる点が領域SS内にあるときにはスリップ制御実行条
件が成立しているか否かが判断され、スリップ制御実行
条件が成立していると判断されたときにはスリップ制御
が行われ、スリップ制御実行条件が成立していないと判
断されたときにはロックアップ機構33がオフにされ
る。また、DEPおよびVVにより定まる点がそれ以外
の領域内にあるときにはロックアップがオフにされる。
一方、車両が悪路を走行していると判断されたときにD
EPおよびVVにより定まる点が領域SO内にあるとき
にはロックアップ機構33がオフにされ、DEPおよび
VVにより定まる点が領域SSO内にあるときにはスリ
ップ制御実行条件が成立しているか否かが判断され、ス
リップ制御実行条件が成立していると判断されたときに
はスリップ制御が行われ、スリップ制御実行条件が成立
していないと判断されたときにはロックアップ機構33
がオフにされる。なお、本実施例では例えば機関回転数
Nまたは機関負荷Q/Nがが一定値よりも高いときにス
リップ制御実行条件が成立していないと判断され、それ
以外はスリップ制御実行条件が成立していると判断され
る。
【0062】次に図15から図26を参照しつつ各気筒
のトルク変動量を求めるためのルーチンについて説明す
る。なお、図26は各ルーチンにおいて行われる各値の
計算タイミングを示している。図15は30°クランク
角度毎に行われる割込みルーチンを示している。図15
を参照するとまず初めに経過時間差DTa(i)および
経過時間Tb(i)を算出するためのルーチン(ステッ
プ100)に進む。このルーチンは図16から図19に
示されている。次いでトルク変動の算出を許可するか否
かをチェックするためのルーチン(ステップ200)に
進む。このルーチンは図20から図22に示されてい
る。次いでトルク変動を算出するためのルーチン(ステ
ップ300)に進む。このルーチンは図24に示されて
いる。次いでトルク変動値の算出に用いるカウンタCD
LNIXの処理ルーチンに進む。このルーチンは図25
に示されている。
【0063】経過時間差DTa(i)および経過時間T
b(i)の算出ルーチンを示す図16から図19を参照
すると、まず初めにステップ101において時刻TIM
EがTIMEOとされる。電子制御ユニット20は時刻
を表わすフリーランカウンタを備えており、このフリー
ランカウンタのカウント値から時刻TIMEが算出され
る。次いでステップ102では現在の時刻TIMEが取
込まれる。従ってステップ101のTIMEOは30°
クランク角度前の時刻を表わしていることになる。
【0064】次いでステップ103では現在i番気筒の
ATDC30°であるか否かが判別される。現在i番気
筒のATDC30°でない場合にはステップ111にジ
ャンプして現在i番気筒のATDC90°であるか否か
が判別される。現在i番気筒のATDC90°でない場
合には経過時間差DTa(i)および経過時間Tb
(i)の算出ルーチンを完了する。
【0065】これに対してステップ103において現在
i番気筒のATDC30°であると判別されたときには
ステップ104に進んで720°クランク角度前に算出
された経過時間Ta(i)がTaO(i)とされる。次
いでステップ105では次式に基づいてi番気筒のTD
CからATDC30°までの最終的な経過時間Ta
(i)が算出される。
【0066】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°クランク角度回転する間の経過時間を正確に表
わしていることになる。
【0067】次いでステップ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−3−4−2である
ので現在1番気筒#1のATDC30°であるとすると
一つ前に燃焼が行われた2番気筒#2のトルク変動量を
算出すべきことを示すフラグXCAL(2)がセットさ
れる。同様に図26に示される如く最終的な経過時間T
a(3)が算出されるとフラグXCAL(1)がセット
され、最終的な経過時間Ta(4)が算出されるとフラ
グXCAL(3)がセットされ、最終的な経過時間Ta
(2)が算出されるとフラグXCAL(4)がセットさ
れる。
【0068】一方、ステップ111において現在i番気
筒のATDC90°であると判別されたときにはステッ
プ112に進んで720°クランク角度前に算出された
経過時間Tb(i)がTbO(i)とされる。次いでス
テップ113では次式に基づいてi番気筒のATDC6
0°からATDC90°までの最終的な経過時間Tb
(i)が算出される。
【0069】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)の算出ルーチンを完了する。
【0070】一方、ステップ107において現在2番気
筒#2のATDC30°であると判別されたときにはス
テップ108に進んで次式に基づき図12の720°ク
ランク角度範囲(j)における経過時間Ta(i)の平
均値TaAVj が算出される。 TaAVj =(Ta(1)+Ta(3)+Ta(4)+
Ta(2))/4 ここでTa(1),Ta(3),Ta(4),Ta
(2)は図12におけるTa(1)j ,Ta(3)j
Ta(4)j ,Ta(2)j に夫々相当する。次いでス
テップ109の補正係数KTa(i)の算出ルーチンに
進む。このルーチンは図18に示されている。
【0071】図18を参照するとまず初めにステップ1
20において減速運転中の燃料供給停止時であるか否か
が判別される。減速運転中の燃料供給停止時でないとき
には処理サイクルを完了し、減速運転中の燃料供給停止
時にはステップ121に進む。ステップ121ではステ
ップ108において算出された経過時間Ta(i)の平
均値TaAVj と、既に算出されている図12の720
°クランク角度範囲(j−1)における経過時間Ta
(i)の平均値TaAVj-1 (=(Ta(1)+Ta
(3)+Ta(4)+Ta(2))/4)(ここでTa
(1),Ta(3),Ta(4),Ta(2)は夫々図
12に示されるTa(1)j-1 ,Ta(3) j-1 ,Ta
(4)j-1 ,Ta(2)j-1 に相当する)から、平均的
な一定の増大率でもって増大したと考えたときの図12
に示す経過時間AAV(4),AAV(2),AAV
(1),AAV(3)が次のようにして算出される。
【0072】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)は夫々
図12に示されるTa(4)j-1 ,Ta(2)j-1 ,T
a(1)j ,Ta(3)j に相当する。
【0073】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)の学習
条件が成立しているか否かが判別される。例えば燃料供
給停止から一定時間経過しており、機関回転数が予め定
められた範囲内にあり、かつ車両が悪路を走行していな
いとき、例えば図13に示される経過時間差DTa
(i)の振幅AMPが設定値を越えていないときに学習
条件が成立していると判断される。学習条件が成立して
いないときにはステップ128に進んでずれ量KAFC
(i)の積算値KAFCI(i)が零とされ、次いでス
テップ129において積算カウント値CKAFCがクリ
アされる。
【0074】これに対して学習条件が成立している場合
にはステップ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)が算出される。
【0075】ステップ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が算出さ
れる。
【0076】tKTaM=(KTaB(1)+KTaB
(2)+KTaB(3)+KTaB(4))/4 次いで次式に示されるように基本補正係数KTaB
(i)をその平均値tKTaMでもって修正することに
より補正係数KTa(i)が更新される。 KTa(i)=KTaB(i)−tKTaM このようにして減速運転中の燃料噴射停止時に各気筒に
対する補正係数KTa(i)の更新が行われる。
【0077】一方、図17のステップ114において現
在2番気筒#2のATDC90°であると判別されたと
きにはステップ115に進んでTaAVj を求めたとき
と同じ方法により次式に基づいて720°クランク角度
範囲における経過時間Tb(i)の平均値TbAVj
算出される。 TbAVj =(Tb(1)+Tb(3)+Tb(4)+
Tb(2))/4 次いでステップ116の補正係数KTb(i)の算出ル
ーチンに進む。このルーチンは図19に示されている。
【0078】図19を参照するとまず初めにステップ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)が算出される。
【0079】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)が算出される。
【0080】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)の学習
条件が成立しているか否かが判別される。前述したよう
に、例えば燃料供給停止から一定時間経過しており、機
関回転数が予め定められた範囲内にあり、かつ車両が悪
路を走行していないとき、例えば図13に示される経過
時間差DTa(i)の振幅AMPが設定値を越えていな
いときに学習条件が成立していると判断される。学習条
件が成立していないときにはステップ138に進んでず
れ量KBFC(i)の積算値KBFCI(i)が零とさ
れ、次いでステップ139において積算カウント値CK
BFCがクリアされる。
【0081】これに対して学習条件が成立している場合
にはステップ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)が算出される。
【0082】ステップ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が算出さ
れる。
【0083】tKTbM=(KTbB(1)+KTbB
(2)+KTbB(3)+KTbB(4))/4 次いで次式に示されるように基本補正係数KTbB
(i)をその平均値tKTbMでもって修正することに
より補正係数KTb(i)が更新される。 KTb(i)=KTbB(i)−tKTbM このようにして減速運転中の燃料噴射停止時に各気筒に
対する補正係数KTb(i)の更新が行われる。
【0084】次に図20から図22に示されるトルク変
動算出許可チェックルーチンについて図23を参照しつ
つ説明する。このルーチンは車両が凸凹道を走行するこ
とにより経過時間差DTa(i)の変動の振幅AMP
(図13)が大きくなったときには特定の気筒について
のトルク変動量の算出を禁止するために設けられてい
る。
【0085】即ち、図20から図22を参照すると、ま
ず初めにステップ201において現在いずれかの気筒の
ATDC30°であるか否かが判別される。現在いずれ
かの気筒のATDC30°でないときには処理サイクル
を完了し、現在いずれかの気筒のATDC30°である
ときにはステップ202に進む。ステップ202からス
テップ204では経過時間差DTa(i)が増大し次い
で減少する際の最大経過時間差DT30maxが算出さ
れる。即ち、ステップ202では図16のステップ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に進む。
【0086】ステップ205からステップ207では経
過時間差DTa(i)が減少し次いで増大する際の最小
経過時間差DT30minが算出される。即ち、ステッ
プ205では図16のステップ106において算出され
たDTa(i)が最小経過時間差DT30minよりも
小さいか否かが判別される。DT30min<DTa
(i)のときにはステップ208にジャンプし、これに
対してDT30min≧DTa(i)のときにはステッ
プ206に進んでDTa(i)がDT30minとされ
る。次いでステップ207ではDTa(i)が減少して
いることを示す減少フラグXMNRECがセット(XM
NREC←“1”)され、次いでステップ208に進
む。
【0087】ステップ208からステップ214ではD
Ta(i)の変動の振幅AMP(図13)が設定値A0
を越えたときにはDTa(i)が最大となった気筒につ
いてのトルク変動量の算出を禁止する禁止フラグがセッ
トされる。即ち、ステップ208ではDT30max>
DTa(i)でかつXMXREC=“1”であるか否か
が判別される。DT30max≦DTa(i)である
か、又は増大フラグXMXRECがリセット(XMXR
EC=“0”)されているときにはステップ215にジ
ャンプし、これに対してDT30max>DTa(i)
でかつXMXREC=“1”のときにはステップ209
に進む。
【0088】即ち、図23に示されるように時刻t1
おいて1番気筒#1の経過時間差DTa(1)が最大に
なったとする。この場合、時刻t1 において行われる割
込みルーチンではステップ202からステップ203に
進んでDTa(1)がDT30maxとされ、次いでス
テップ204において最大フラグXMXRECがセット
される。一方、図23の時刻t2 において行われる割込
みルーチンではステップ202からステップ205にジ
ャンプする。このときステップ208ではDT30ma
x>DTa(3)であり、かつXMXREC=“1”で
あると判断されるのでステップ209に進む。即ち、ス
テップ209に進むのは経過時間差DTa(i)が減少
しはじめる時刻t2 である。
【0089】ステップ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”)される。
即ち、図23の時刻t2 において行われる割込みルーチ
ンでは前述したように1番気筒#1のトルク変動量が算
出される。従ってこの割込みルーチンにおいてAMP≧
0 となり、トルク変動算出禁止フラグXNOCALが
セットされると1番気筒#1のトルク変動量の算出、即
ち、DTa(i)が最大となる気筒のトルク変動量の算
出が禁止される。
【0090】ステップ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に進
む。
【0091】即ち、図23に示されるように時刻t3
おいて1番気筒#1の経過時間差DTa(1)が最小に
なったとする。この場合、時刻t3 において行われる割
込みルーチンではステップ205からステップ206に
進んでDTa(1)がDT30minとされ、次いでス
テップ207において減少フラグXMNRECがセット
される。一方、図23の時刻t4 において行われる割込
みルーチンではステップ205からステップ208にジ
ャンプする。このときステップ215ではDT30mi
n<DTa(3)であり、かつXMNREC=“1”で
あると判断されるのでステップ216に進む。即ち、ス
テップ216に進むのは経過時間差DTa(i)が増大
しはじめる時刻t4 である。
【0092】ステップ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”)される。即ち、図23の時刻t4 において行わ
れる割込みルーチンでは1番気筒#1のトルク変動量が
算出される。従ってこの割込みルーチンにおいてAMP
≧A0となり、トルク変動算出禁止フラグXNOCAL
がセットされると1番気筒#1のトルク変動量の算出、
即ち、DTa(i)が最小となる気筒のトルク変動量の
算出が禁止される。
【0093】ステップ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)−DT
a(i−1)|<Ko・|DTa(i−1)−DTa
(i)|であると判別されたときには処理ルーチンを完
了し、|DTa(i−2)−DTa(i−1)|≧Ko
・|DTa(i−1)−DTa(i)|であると判別さ
れたときにはステップ223に進んでトルク変動算出禁
止フラグXNOCALがセットされる。即ち、今図23
の時刻t3 における割込みルーチンであるとするとこの
ときには|DTa(4)−DTa(2)|≧Ko ・|D
Ta(2)−DTa(1)|であるか否かが判別され
る。図23に示されるようにDTa(4)に対してDT
a(2)が急変すると|DTa(4)−DTa(2)|
はKo ・|DTa(2)−DTa(1)|よりも大きく
なる。このときトルク変動算出禁止フラグがセットさ
れ、経過時間差DTa(i)が急変した2番気筒#2の
トルク変動量の算出が禁止される。
【0094】次に図24に示すトルク変動算出ルーチン
について説明する。図24を参照すると、まず初めにス
テップ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に進む。
【0095】ステップ303では一つ前に燃焼が行われ
た気筒についてのトルク変動量の算出を禁止する禁止フ
ラグXNOCALがリセット(XNOCAL=“0”)
されているか否かが判別される。この禁止フラグがセッ
ト(XNOCAL=“1”)されているときにはステッ
プ311に進んで禁止フラグXNOCALがリセットさ
れる。これに対して禁止フラグがリセットされていると
きにはステップ304に進む。即ち、フラグXCALが
セットされており、かつ禁止フラグXNOCALがリセ
ットされているときのみステップ304に進む。
【0096】ステップ304では一つ前に燃焼が行われ
た(i−1)番気筒の仮想の経過時間差Kb(i−1)
(図10および図11参照)が次式に基づいて算出され
る。 Kb(i−1)={2・DTa(i−1)+DTa
(i)}/3 次いでステップ305では次式に基づいて一つ前に燃焼
が行われた(i−1)番気筒の実際の出力トルクDN
(i−1)が算出される。
【0097】DN(i−1)={30°/Tb(i−
1)}2−{30°/Ta(i−1)}2 次いでステップ306では次式に基づいて一つ前に燃焼
が行われた(i−1)番気筒の仮想の出力トルクDNS
(i−1)が算出される。 DNS(i−1)={30°/(TbO(i−1)+K
b(i−1))}2−{30°/Ta(i−1)}2 次いでステップ307では次式に示すように仮想の出力
トルクDNS(i−1)から実際の出力トルクDN(i
−1)を減算することによって一つ前に燃焼が行われた
(i−1)番気筒のトルク変動量DLN(i−1)が算
出される。
【0098】DLN(i−1)=DNS(i−1)−D
N(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)が
算出される。
【0099】なお、各気筒の出力変動量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)を算出する。
【0100】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に進む。なお、各気筒のト
ルクは上昇と低下を繰返し、この場合トルク変動量を求
めるにはトルクの上昇分かトルクの減少分のいずれかを
積算すればよい。図24に示すルーチンではトルクの減
少分のみを積算するようにしており、従って上述したよ
うにDLN(i−1)<0のときにはDLN(i−1)
を零にしている。
【0101】次に図25を参照しつつカウンタCDLN
IXの処理について説明する。このカウンタCDLNI
Xのカウント値は後に説明する機関のトルク変動値を算
出する際に使用される。図25を参照すると、まず初め
にステップ401において現在3番気筒#3のATDC
30°であるか否かが判別される。現在3番気筒#3の
ATDC30°でないときには処理サイクルを完了し、
現在3番気筒#3のATDC30°であるときにはステ
ップ402に進む。ステップ402では機関のトルク変
動値を算出するためのトルク変動値算出条件が成立して
いるか否かが判別される。例えば空燃比をリーンとする
条件が成立していないか、或いは機関負荷Q/Nの単位
時間当りの変化量ΔQ/Nが設定値以上であるか、或い
は機関回転数の単位時間当りの変化量ΔNが設定値以上
であるときにはトルク変動値算出条件が成立していない
と判断され、それ以外のときにはトルク変動値算出条件
が成立していると判断される。
【0102】ステップ402においてトルク変動値算出
条件が成立していると判断されたときにはステップ40
8に進んでカウント値CDLNIXが1だけインクリメ
ントされる。このカウント値CDLNIXのインクリメ
ント作用は3番気筒#3がATDC30°となる毎に、
即ち720°クランク角度毎に行われる。次いでステッ
プ409ではカウント値CDLNIXのインクリメント
作用が開始されてからカウント値CDLNIXがクリア
されるまでの間の機関回転数の平均値NAVE および吸入
空気量Qの平均値QAVE が算出される。
【0103】一方、ステップ402においてトルク変動
値算出条件が成立していないと判断されたときにはステ
ップ403に進んでカウント値CDLNIXがクリアさ
れる。次いでステップ404では各気筒に対するトルク
変動量DLN(i)の積算値DLNI(i)(この積算
値は後に説明するルーチンにおいて算出される)がクリ
アされ、次いでステップ405では各気筒に対する積算
カウント値CDLNI(i)(この積算カウント値は後
に説明するルーチンにおいて算出される)がクリアされ
る。
【0104】次いでステップ406では目標トルク変動
値LVLLFBが算出される。本実施例では後に説明す
るように算出されたトルク変動値がこの目標トルク変動
値LVLLFBとなるように空燃比がフィードバック制
御される。この目標トルク変動値LVLLFBは等しい
変動値を実線で示した図27(A)に示されるように機
関負荷Q/Nが高くなるほど大きくなり、機関回転数N
が高くなるほど大きくなる。この目標トルク変動値LV
LLFBは図27(B)に示されるように機関負荷Q/
Nおよび機関回転数Nの関数としてマップの形で予めR
OM22内に記憶されている。次いでステップ407で
は各気筒のトルク変動値DLNISM(i)(このトル
ク変動値は後に説明するルーチンにおいて算出される)
が図27(B)のマップから算出された目標トルク変動
値LVLLFBとされる。
【0105】図28は繰返し実行されるメインルーチン
を示している。このメインルーチンではまず初めにロッ
クアップ機構33の制御ルーチン(ステップ500)が
実行される。このルーチンが図29に示されている。次
いでトルク変動値の算出ルーチン(ステップ600)が
実行される。このルーチンが図31および図32に示さ
れている。次いでリーンリミットフィードバック補正係
数FLLFBの算出ルーチン(ステップ700)が実行
される。このルーチンが図33に示されている。次いで
予め定められたクランク角になったときに噴射時間算出
ルーチン(ステップ800)が実行される。このルーチ
ンが図35に示されている。次いでその他のルーチン
(ステップ900)が実行される。
【0106】まず図29に示されるロックアップ機構3
3の制御ルーチンについて説明する。図29を参照する
と、まず初めにステップ501ではアクセルペダルの踏
み込み量DEPおよび車速VVにより定まる点が図14
に示される領域SO内または領域SSO内にあるか否か
が判別される。DEPおよびVVにより定まる点が図1
4の領域SO内または領域SSO内にないと判別された
ときには次いでステップ502に進み、DEPおよびV
Vにより定まる点が図14の領域SS内にあるか否かが
判別される。DEPおよびVVにより定まる点が図14
の領域SS内にあると判別されたときには次いでステッ
プ503に進み、スリップ制御実行条件が成立している
か否かが判別される。上述したように本実施例では例え
ば機関回転数Nまたは機関負荷Q/Nがが一定値よりも
高いときにスリップ制御実行条件が成立していないと判
別され、それ以外はスリップ制御実行条件が成立してい
ると判別される。スリップ制御実行条件が成立している
と判別されたときには次いでステップ504に進み、ロ
ックアップクラッチのスリップ制御が行われる。これに
対し、スリップ制御実行条件が成立していないと判別さ
れたときには次いでステップ505に進んでロックアッ
プがオフとされる。一方、ステップ502においてDE
PおよびVVにより定まる点が図14の領域SS内にな
いと判別されたときには次いでステップ505に進んで
ロックアップ機構33がオフとされる。
【0107】ステップ501においてDEPおよびVV
により定まる点が図14の領域SO内または領域SSO
内にあると判別されたときには次いでステップ506に
進み、リーンリミットフィードバック補正係数FLLF
B(後述するルーチンで算出される)が設定値F1より
も大きくかつ振幅AMP(図20から図22のルーチン
で算出される)が設定値A1よりも大きいか否かが判別
される。FLLFB>F1かつAMP>A1のときには
次いでステップ507に進んでカウント値CSINPが
1だけインクリメントされ、FLLFB≦F1またはA
MP≦A1のときには次いでステップ508に進んでカ
ウント値CSINPが1だけデクリメントされる。この
カウント値CSINPはFLLFB>F1かつAMP>
A1となっている時間を表している。次いでステップ5
09ではカウント値CSINPが設定値CSINP1よ
りも大きいか否かが判別される。CSINP≦CSIN
P1のときには次いでステップ510に進んでロックア
ップ機構33がオンとされ、或いはオンに維持される。
これに対し、CSINP>CSINP1のときには次い
でステップ503に進み、スリップ制御が実行されるか
或いはロックアップ機構33がオフとされる。
【0108】即ち、ロックアップ機構33がオンのとき
に図30の時刻t1 において車両が悪路を走行し始めた
とするとこのときからリーンリミットフィードバック補
正係数FLLFBが上昇し始め、振幅AMPも上昇し始
める。次いで図30の時刻t 2 に示されるように振幅A
MPが設定値A1よりも大きくなりかつリーンリミット
フィードバック補正係数FLLFBが設定値F1よりも
大きくなるとカウント値CSINPのインクリメント作
用が開始される。次いで図30の時刻t3 に示されるよ
うにカウント値CSINPが設定値CSINP1よりも
大きくなるとロックアップ機構33がオンになるのが禁
止される。その結果、リーンリミットフィードバック補
正係数FLLFBおよび振幅AMPが低下し始める。こ
のようにリーンリミットフィードバック補正係数FLL
FBが低下すると空燃比がリーン側に移行することにな
る。次いで図30の時刻t4 に示されるようにAMP≦
A1となるとカウント値CSINPのデクリメント作用
が開始される。次いで図30の時刻t5 に示されるよう
にCSINP≦CSINP1となるとロックアップ機構
33がオンになることの禁止作用が解除され、即ちロッ
クアップ機構33が再びオンとされる。
【0109】このように本実施例ではリーンリミットフ
ィードバック補正係数FLLFBおよび振幅AMPに基
づいて車両が悪路を走行しているか否かを判断している
ので車両が悪路を走行しているか否かの判断を正確にか
つ速やかに行うことができる。即ち、例えば燃焼圧変動
に基づきリーンリミットフィードバック補正係数FLL
FBが設定値F1よりも大きくなったとしても振幅AM
Pは設定値A1よりも大きくならず、従って設定値F1
を小さく定めることができる。設定値F1を小さく定め
ることができるとリーンリミットフィードバック補正係
数FLLFBが設定値F1よりも大きくなるまでの時
間、即ち車両が実際に悪路を走行し始めてから、車両が
悪路を走行していると判断されるまでの時間を短くする
ことができ、斯くして燃焼消費率およびNOX の発生量
を十分に低減することができる。
【0110】なお、ステップ506においてFLLFB
≦F1またはAMP≦A1のときにはカウント値CSI
NPをクリアするようにしてもよい。しかしながら、本
実施例のようにFLLFB>F1かつAMP>A1のと
きにカウント値CSINPをインクリメントし、FLL
FB≦F1またはAMP≦A1のときにカウント値CS
INPをデクリメントしてこのカウント値CSINPと
設定値CSINP1とを比較するようにすると振幅AM
Pの平均値を求めてこの平均値と設定値CSINP1と
を比較することと同じになる。従って、振幅AMPが比
較的大きく変動するときにもロックアップ機構33がオ
ンになるのを禁止すべきか否かを正確に判断することが
できる。
【0111】次に図31および図32に示されるトルク
変動値の算出ルーチンについて説明する。図31および
図32を参照すると、まず初めにステップ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だけ
インクリメントされる。
【0112】次いでステップ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)がリセットさ
れる。
【0113】即ち、算出されたトルク変動量積算値DL
NI(i)とこれまで用いられてきたトルク変動値DL
NISM(i)との間に差があるときにはこの差{DL
NI(i)−DLNISM(i)}に1/4を乗算した
値がトルク変動値DLNISM(i)に加算される。従
って例えば1番気筒#1についての積算カウント値CD
LNI(1)が“8”になるとステップ606において
トルク変動値DLNISM(1)が算出されることにな
る。
【0114】次いでステップ609では図25に示すル
ーチンにおいて算出されたカウント値CDLNIXが
“8”になったか否かが判別される。CDLNIXが
“8”でないときには処理サイクルを完了し、CDLN
IXが“8”になるとステップ610に進んで各気筒の
トルク変動値DLNISM(i)の平均値である平均ト
ルク変動値DLNISM(={DLNISM(1)+D
LNISM(2)+DLNISM(3)+DLNISM
(4)}/4)が算出される。次いでステップ611で
はカウント値CDLNIXがクリアされる。このように
して機関のトルク変動量を代表する値DLNISMが算
出される。
【0115】なお、前述したようにカウント値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)が使
用される。
【0116】次に図33を参照しつつFLLFB算出ル
ーチンについて説明する。図33を参照すると、まず初
めにステップ701においてリーンリミットフィードバ
ック補正係数FLLFBの更新条件が成立しているか否
かが判別される。例えば暖機運転時であるとき、或いは
機関の運転状態が図5において破線で囲まれた学習領域
にないときには更新条件が成立していないと判断され、
その他のときには更新条件が成立していると判断され
る。更新条件が成立していないときには処理サイクルを
完了し、更新条件が成立しているときにはステップ70
2に進む。
【0117】ステップ702では機関負荷Q/Nと機関
回転数Nから図27(B)に示すマップに基づいて目標
トルク変動値LVLLFBが算出される。次いでステッ
プ703およびステップ704では目標トルク変動値L
VLLFBに応じた変動量判別値DH(n),DL
(n)に基づいて次式に示されるトルク変動レベルLV
LH(n),LVLL(n)が算出される。
【0118】 LVLH(n)=LVLLFB+DH(n) LVLL(n)=LVLLFB+DL(n) ここで、変動量判別値DH(n)およびDL(n)は図
34(A)に示されるように予め定められている。即
ち、図34(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が大きくなる
につれて次第に増大する。
【0119】ところで今、ステップ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)とされる。
【0120】一方、図34(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 である。図34(B)に示
す各フィードバック補正値+a1 ,+a2 ,+a3 ,+
4 ,−b1 ,−b 2 ,−b3 ,−b4 が図34(A)
の対応する領域に示されている。
【0121】ステップ703およびステップ704にお
いて夫々トルク変動レベルLVLH(n),LVLL
(n)が算出されるとステップ705に進んで図31お
よび図32に示すトルク変動値の算出ルーチンにより算
出された平均トルク変動値DLNISMが図34(B)
に示されるどのトルク変動レベルLVLH(n),LV
LL(n)の間にあるかが判別される。次いでステップ
706では対応するフィードバック補正値DLFBが算
出される。例えば今、目標トルク変動値LVLLFBが
図34(A)において破線で示される値であり、算出さ
れた平均トルク変動値DLNISMが図34(B)のL
VLH(1)とLVLH(2)との間である場合、即ち
目標トルク変動値LVLLFBに対する平均トルク変動
値DLNISMの偏差が図34(A)の破線上において
DH(1)とDH(2)の間にある場合にはフィードバ
ック補正値DLFBは+a2 とされる。
【0122】次いでステップ707では図25に示すC
DLNIXの処理ルーチンのステップ409において求
められた機関回転数の平均値NAVE および吸入空気量Q
の平均値QAVE に基づいて更新すべきリーンリミットフ
ィードバック補正係数FLLFBijが図5に示されるど
の学習領域のリーンリミットフィードバック補正係数で
あるかが決定される。次いでステップ708ではステッ
プ707において決定されたリーンリミットフィードバ
ック補正係数FLLFBijにフィードバック補正値DL
FBが加算される。
【0123】即ち、上述したように例えばDLNISM
>LVLLFBであって、LVLH(1)<DLNIS
M<LVLH(2)である場合にはリーンリミットフィ
ードバック補正係数FLLFBijに+a2 が加算され
る。その結果、空燃比が小さくなるので各気筒のトルク
変動量が減少せしめられる。一方、DLNISM<LV
LLFBであってLVLL(1)>DLNISM>LV
LL(2)である場合にはリーンリミットフィードバッ
ク補正係数FLLFBijに−b2 が加算される。その結
果、空燃比が大きくなるので各気筒のトルク変動量が増
大せしめられる。このようにして全気筒の平均トルク変
動値DLNISMが目標トルク変動値LVLLFBとな
るようにリーン運転時の空燃比が制御される。
【0124】なお、図25に示すルーチンに示されるよ
うにトルク変動値の算出条件が成立しないときにはステ
ップ407において全ての気筒に対するDLNISM
(i)がLVLLFBとされ、斯くして平均トルク変動
値DLNISMも目標トルク変動値LVLLFBとされ
る。従ってこのときにはリーンリミットフィードバック
補正係数FLLFBijの更新は行われない。
【0125】次に図35を参照しつつ燃料噴射時間の算
出ルーチンについて説明する。図35を参照すると、ま
ず初めにステップ801において図2に示すマップから
基本燃料噴射時間TPが算出される。次いでステップ8
02ではリーン運転を行うべき運転状態か否かが判別さ
れる。リーン運転を行うべき運転状態のときにはステッ
プ803に進んで理論空燃比フィードバック補正係数F
AFの値が1.0に固定される。次いでステップ804
では図4に示すマップからリーン補正係数FLEANが
算出され、次いでステップ805では図5に示すマップ
からリーンリミットフィードバック補正係数FLLFB
が読込まれる。次いでステップ809では次式に基づい
て燃料噴射時間TAUが算出される。
【0126】TAU=TP・FLEAN・FLLFB・
FAF+TAUV これに対し、ステップ802においてリーン運転を行う
べき運転状態でないと判別されたとき、即ち空燃比を理
論空燃比にすべきときにはステップ806に進んでリー
ン補正係数FLEANが1.0に固定され、次いでステ
ップ807においてリーンリミットフィードバック補正
係数FLLFBが1.0に固定される。次いでステップ
808では空燃比センサ17の出力信号に基づいて空燃
比が理論空燃比となるように理論空燃比フィードバック
補正係数FAFが制御される。次いでステップ809に
進み、燃料噴射時間TAUが算出される。
【0127】これまで述べてきた実施例では経過時間差
DTa(i)の変動の振幅AMP(図13参照)を車両
振動を代表する代表値としている。しかしながら、経過
時間Ta(i)の変動の振幅AMQ(図7参照)を車両
振動を代表する代表値とし、設定時間以上、リーンリミ
ットフィードバック補正係数FLLFBが第1の設定値
よりも大きくかつ振幅AMQが第2の設定値よりも大き
いときに車両が悪路を走行していると判断してロックア
ップ機構33がオンになるのを禁止するようにすること
もできる。
【0128】
【発明の効果】車両が悪路を走行しているか否かを速や
かにかつ正確に判断することができるので車両の良好な
運転性を維持しつつ燃焼消費率およびNOX の発生量を
十分に低減することができる。
【図面の簡単な説明】
【図1】内燃機関の全体図である。
【図2】基本燃料噴射時間のマップを示す図である。
【図3】NOx の発生量とトルク変動を示す図である。
【図4】リーン補正係数のマップを示す図である。
【図5】リーンリミットフィードバック補正係数のマッ
プを示す図である。
【図6】30°クランク角度の経過時間Ta(i),T
b(i)の変化を示すタイムチャートである。
【図7】30°クランク角度の経過時間Ta(i)の変
化を示すタイムチャートである。
【図8】30°クランク角度の経過時間Ta(i),T
b(i)の変化を示すタイムチャートである。
【図9】30°クランク角度の経過時間Ta(i),T
b(i)の変化を示すタイムチャートである。
【図10】30°クランク角度の経過時間Ta(i)の
変化を示すタイムチャートである。
【図11】経過時間差DTa(i)とKb(i)との関
係を示す図である。
【図12】減速運転時における経過時間Ta(i)の変
化を示すタイムチャートである。
【図13】経過時間差DTa(i)の変化を示すタイム
チャートである。
【図14】ロックアップ機構をオンにすべき運転領域お
よびロックアップクラッチのスリップ制御を行うべき運
転領域を示す線図である。
【図15】割込みルーチンを示すフローチャートであ
る。
【図16】経過時間差DTa(i)および経過時間Tb
(i)を算出するためのフローチャートである。
【図17】経過時間差DTa(i)および経過時間Tb
(i)を算出するためのフローチャートである。
【図18】KTa(i)を算出するためのフローチャー
トである。
【図19】KTb(i)を算出するためのフローチャー
トである。
【図20】トルク変動量算出の許可をチェックするため
のフローチャートである。
【図21】トルク変動量算出の許可をチェックするため
のフローチャートである。
【図22】トルク変動量算出の許可をチェックするため
のフローチャートである。
【図23】経過時間差DTa(i)の変化とフラグXM
XREC,XMNRECの変化を示すタイムチャートで
ある。
【図24】トルク変動量を算出するためのフローチャー
トである。
【図25】カウンタCDLNIXを処理するためのフロ
ーチャートである。
【図26】種々の値の計算タイミングを示す図である。
【図27】目標トルク変動値を示す図である。
【図28】メインルーチンを示すフローチャートであ
る。
【図29】ロックアップ機構を制御するためのフローチ
ャートである。
【図30】リーンリミットフィードバック補正係数FL
LFB、振幅AMP、およびカウント値CSINPの変
化を示すタイムチャートである。
【図31】トルク変動値を算出するためのフローチャー
トである。
【図32】トルク変動値を算出するためのフローチャー
トである。
【図33】リーンリミットフィードバック補正係数を算
出するためのフローチャートである。
【図34】変動量判別値DH(n),DL(n)および
トルク変動レベルLVLH(n),LVLL(n)を示
す図である。
【図35】燃料噴射時間を算出するためのフローチャー
トである。
【符号の説明】
3…サージタンク 4…燃料噴射弁 8…スロットル弁 14…ロータ 15…クランク角センサ 33…ロックアップ機構
───────────────────────────────────────────────────── フロントページの続き (51)Int.Cl.6 識別記号 FI // F16H 59:34

Claims (2)

    【特許請求の範囲】
  1. 【請求項1】 自動変速機を備え、自動変速機のトルク
    コンバータ内にロックアップ機構を設け、クランクシャ
    フトの角速度に基づいて各気筒の出力変動を検出し、該
    出力変動に基づいて基準空燃比をフィードバック補正係
    数により補正するようにした内燃機関において、クラン
    クシャフトの角速度変動に基づいて車両振動を代表する
    代表値を算出し、ロックアップ機構がオンのときに設定
    時間以上、フィードバック補正係数による基準空燃比の
    リッチ側への補正量が第1の設定値よりも大きくかつ代
    表値が第2の設定値よりも大きいときにロックアップ機
    構がオンになるのを禁止するようにした内燃機関の動力
    伝達制御装置。
  2. 【請求項2】 圧縮行程末期から爆発行程初期までのク
    ランク角度領域内に設定されたクランク角度範囲内にお
    けるクランクシャフトの角速度を検出し、該角速度の変
    動に基づいて前記車両振動の代表値を算出するようにし
    た請求項1に記載の内燃機関の動力伝達制御装置。
JP9257029A 1997-09-22 1997-09-22 内燃機関の動力伝達制御装置 Pending JPH1191409A (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP9257029A JPH1191409A (ja) 1997-09-22 1997-09-22 内燃機関の動力伝達制御装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP9257029A JPH1191409A (ja) 1997-09-22 1997-09-22 内燃機関の動力伝達制御装置

Publications (1)

Publication Number Publication Date
JPH1191409A true JPH1191409A (ja) 1999-04-06

Family

ID=17300762

Family Applications (1)

Application Number Title Priority Date Filing Date
JP9257029A Pending JPH1191409A (ja) 1997-09-22 1997-09-22 内燃機関の動力伝達制御装置

Country Status (1)

Country Link
JP (1) JPH1191409A (ja)

Similar Documents

Publication Publication Date Title
JP3303739B2 (ja) 内燃機関の空燃比制御方法
US6199426B1 (en) Method of detection of output fluctuation in internal combustion engine
JP3303732B2 (ja) 内燃機関の制御装置
JP3216577B2 (ja) 内燃機関の空燃比制御方法
JP4135504B2 (ja) 内燃機関の制御装置
JP3262003B2 (ja) 多気筒内燃機関の出力変動検出方法
JPH1182086A (ja) 内燃機関の空燃比制御方法
JPH1191409A (ja) 内燃機関の動力伝達制御装置
EP0811758B1 (en) Method of controlling an air-fuel ratio of an engine
JP2004116366A (ja) 内燃機関のアイドル回転数制御装置
JP3085220B2 (ja) 内燃機関の空燃比制御方法
EP0799983B1 (en) Method of detection of angular velocity and torque in an internal combustion engine
JP3303669B2 (ja) 内燃機関の空燃比制御方法
JPH1182119A (ja) 内燃機関の出力制御装置
JP3279179B2 (ja) 多気筒内燃機関における検出方法
JP3246325B2 (ja) 内燃機関における検出方法
JP3246328B2 (ja) 内燃機関における検出方法
JP3156588B2 (ja) 内燃機関の空燃比制御方法
JPH10331685A (ja) 内燃機関の空燃比制御装置
JPH1130148A (ja) 内燃機関の空燃比制御方法
JP3218970B2 (ja) 内燃機関における検出方法
JPH09281006A (ja) 内燃機関における変動検出方法
JP2005207366A (ja) 内燃機関の制御装置
JP2519973Y2 (ja) エンジンの制御装置
JP2004108289A (ja) 燃料噴射制御装置