JP4151605B2 - エンジンのノック制御装置 - Google Patents
エンジンのノック制御装置 Download PDFInfo
- Publication number
- JP4151605B2 JP4151605B2 JP2004131299A JP2004131299A JP4151605B2 JP 4151605 B2 JP4151605 B2 JP 4151605B2 JP 2004131299 A JP2004131299 A JP 2004131299A JP 2004131299 A JP2004131299 A JP 2004131299A JP 4151605 B2 JP4151605 B2 JP 4151605B2
- Authority
- JP
- Japan
- Prior art keywords
- knock
- combustion
- ignition timing
- updated
- calculated
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
Landscapes
- Electrical Control Of Ignition Timing (AREA)
- Electrical Control Of Air Or Fuel Supplied To Internal-Combustion Engine (AREA)
- Combined Controls Of Internal Combustion Engines (AREA)
Description
例えば(1)式より目標空燃比が理論空燃比のときTFBYA=1.0となり、目標空燃比が例えば22.0といったリーン側の値であるとき、TFBYAは1.0未満の正の値である。
Vc=(π/4)D2・Hx/(ε−1) …(3)
Hivc={(CND+ST2/2)−(CRoff−PISoff)2}1/2
−{(ST/2)・cos(θivc+θoff)}
+(CND2−X2)1/2 …(4)
X =(ST/2)・sin(θivc+θoff)−CRoff+PISoff …(5)
θoff=arcsin{(CRoff−PISoff)/(CND・(ST/2))}…(6)
ただし、Vc :隙間容積[m3]、
ε :圧縮比、
D :シリンダボア径[m]、
ST :ピストンの全ストローク[m]、
Hivc :吸気弁閉時期におけるピストンピン76の
TDCからの距離[m]、
Hx :ピストンピン76のTDCからの距離の最大値と最小値の
差[m]、
CND :コネクティングロッド74の長さ[m]、
CRoff :結節点75のシリンダ中心軸73からのオフセット距離
[m]、
PISoff:クランクシャフト回転中心72のシリンダ中心軸73から
のオフセット距離[m]、
θivc :吸気弁閉時期のクランク角[degATDC]、
θoff :ピストンピン76とクランクシャフト回転中心72とを結
ぶ線がTDCにおいて垂直線となす角度[deg]、
X :結節点75とピストンピン76との水平距離[m]、
吸気弁閉時期のクランク角θivcは前述のように、エンジンコントローラ31から吸気VTC機構27への指令信号によって決まるので、既知である。式(2)〜(6)にこのときのクランク角θivc(=IVC)を代入すれば、燃焼室5の吸気弁閉時期における容積VIVCを算出することができる。したがって、実用上は燃焼室5の吸気弁閉時期における容積VIVCは吸気弁閉時期IVCをパラメータとするテーブルで設定したものを用いる。吸気VTC機構27を備えないときには定数で与えることができる。
ステップ14では燃焼室5の吸気弁閉時期IVCにおける圧力(つまり圧縮開始時期圧力)PINI[Pa]を算出する。すなわち、吸気弁閉時期IVCになったタイミングでのコレクタ内圧力PCOLを吸気弁閉時期IVCにおける圧力PINIとして取り込む。
具体的に説明すると、MRESFR、TWK、TFBYAの3つのパラメータの組み合わせによって得られる反応確率の最大値を100%とし、これらのパラメータと反応確率RPROBAの関係を実験的に求め、求めた反応確率RPROBAをパラメータに応じたテーブルとしてエンジンコントローラ31のメモリに予め格納しておく。ステップ14ではパラメータに応じてこのテーブルを検索することにより反応確率RPROBAを求める。
具体的にはエンジン回転速度NRPMから、エンジンコントローラ31のメモリに予め格納された図9に示す特性のテーブルを検索することにより基準クランク角θPMAXを求める。算出を容易にするために、基準クランク角θPMAXを一定とみなすことも可能である。
図10、図12はどちらを先に実行してもかまわない。
具体的には前回燃焼開始時期MBTCYCLにおけるピストン6のストローク位置と、燃焼室5のボア径から、燃焼室5のMBTCYCLにおける容積V0を算出する。図5のステップ12では、燃焼室5の吸気弁閉時期IVCにおける容積VIVCを、吸気弁閉時期をパラメータとする吸気弁閉時期容積のテーブルを検索することにより求めたが、ここではMBTCYCLをパラメータとする前回燃焼開始時期容積のテーブルを検索することにより、燃焼室5の前回燃焼開始時期MBTCYCLにおける容積V0を求めればよい。
=V0/VIVC …(12)
ステップ164では吸気弁閉時期IVCから燃焼開始時期に至る間の燃焼室5内の温度上昇率TCOMPを次式に示すように有効圧縮比Ecに基づいて算出する。
ただし、κ:比熱比、
(13)式は断熱圧縮されるガスの温度上昇率の式である。なお、(13)式右辺の「^」は累乗計算を表している。この記号は後述する式でも使用する。
T0=TINI×TCOMP …(14)
の式により算出する。
ただし、κ:比熱比、
(41)式も(13)式と同じに断熱圧縮されるガスの圧力上昇率の式である。(41)式右辺の「^」も(13)式と同じに累乗計算を表している。
P0=PINI×PCOMP …(42)
の式により算出する。
=SLstd×(T0/Tstd)2.18×(P0/Pstd)-0.16 …(15)
ただし、Tstd :基準温度[K]、
Pstd :基準圧力[Pa]、
SLstd:基準温度Tstdと基準圧力Pstdにおける基準層流燃焼
速度[m/sec]、
T0 :燃焼室5の燃焼開始時期における温度[K]、
P0 :燃焼室5の燃焼開始時期における圧力[Pa]、
層流燃焼速度(層流火炎速度)は気体の流れがない状態での火炎の伝播速度のことであり、燃焼室5内の圧縮速度、燃焼室5内の吸気流速に因らず、燃焼室5の温度及び圧力の関数となることが知られていることから、初期燃焼期間における層流燃焼速度を燃焼開始時温度T0と燃焼開始時圧力P0の関数として、また後述するように主燃焼期における層流燃焼速度を圧縮上死点温度TTDCと圧縮上死点圧力PTDCの関数としている。これは、層流燃焼速度は一般的に、エンジン負荷、燃焼室5内の不活性ガス率、吸気弁閉時期、比熱比、吸気温度により変化するのであるが、これらは燃焼室5内の温度Tと圧力Pに影響する因子であるので、層流燃焼速度は最終的に燃焼室5内の温度Tと圧力Pにより規定できるとするものである。
=1.0×0.7×(T0/550)2.18 …(16)
ステップ169では、初期燃焼期間におけるガス流動の乱れ強さST1を算出する。このガス流動の乱れ強さST1は無次元の値であり、燃焼室5に流入する新気の流速と燃料インジェクタ21の噴射燃料のペネトレーションとに依存する。
ただし、C1:定数、
乱れ強さST1を回転速度NRPMをパラメータとするテーブルから求めることも可能である。
燃焼室5内にガス乱れがあるとガスの燃焼速度が変化する。(18)式はこのガス乱れに伴う燃焼速度への寄与(影響)を考慮したものである。
/(RPROBA×AF1×FLAME1) …(19)
ただし、AF1 :火炎核の反応面積(固定値)[m2]、
この(19)式および後述する(22)式は、燃焼ガス質量を燃焼速度で割ると燃焼期間が得られるとする次の基本式より導いたものであるが、(19)、(22)式右辺の分子、分母ががただちに燃焼ガス質量、燃焼速度を表すものではない。
×火炎表面積[m2]×火炎速度[m/sec])
…(補1)
(補1)式右辺分母の未燃ガス密度は、未燃ガス質量[g]を未燃ガス体積[m3]で割った値であるので、従来装置のように質量に相当する充填効率ITACのみの関数では未燃ガス密度を正確に計算できているとはいえない。そこで、(補1)式に対して実験結果とを照らし合わせつつ所定の近似を導入して初めて得られたのが上記(19)式及び後述する(22)式に示す実験式である。
…(43)
(43)式において燃焼室5の圧縮上死点における容積VTDCは運転条件によらず一定であり、予めエンジンコントローラ31のメモリに格納しておけばよい。
=Ec 2^(κ−1)…(44)
ただし、κ:比熱比、
図11と同様の特性のテーブルを予めエンジンコントローラ31のメモリに格納しておき、有効圧縮比Ec 2から当該テーブルを検索することにより温度上昇率TCOMP 2を求めることも可能である。
…(45)
(45)式右辺の括弧内の「1」は新気分、「TFBYA/14.7」は燃料分である。
/(Cv×MGAS)…(46)
ただし、Q :燃料の定発熱量、
BRk:シリンダ内燃料の燃焼質量割合、
Cv :定積比熱、
(46)式右辺の分子はシリンダ内燃料による発生総熱量[J]、分母は単位発生熱量当たりの温度上昇率[J/K]を意味している。すなわち、(46)式は熱力学の公式に当てはめた近似式である。
…(47)
ステップ187では、この燃焼室5の圧縮上死点における温度TTDCと容積VTDC及び燃焼室5の吸気弁閉時期における圧力PINI、容積VIVC及び温度TINIから次式により燃焼室5の圧縮上死点における圧力PTDC[K]を算出する。
…(48)
(48)式は状態方程式を用いて得たものである。すなわち、吸気弁閉時期における圧力、容積及び温度(PINI、VIVC、TINI)を用いて次の状態方程式が成立する。
ただし、n:モル数、
R:ガス定数、
圧縮上死点近傍では容積はほぼ等しいので、圧縮上死点での圧力、容積及び温度(PTDC、VTDC、TTDC)を用いて次の状態方程式が成立する。
この(補3)式と上記(補2)との両式からn・Rを消去しPTDCについて解くと、上記(48)式が得られる。
=SLstd×(TTDC/Tstd)2.18×(PTDC/Pstd)-0.16
…(49)
ただし、Tstd :基準温度[K]、
Pstd :基準圧力[Pa]、
SLstd:基準温度Tstdと基準圧力Pstdにおける基準層流燃焼速度
[m/sec]、
TTDC:燃焼室5の圧縮上死点における温度[K]、
PTDC:燃焼室5の圧縮上死点における圧力[Pa]、
(49)式の解説は上記(16)式と同様ある。すなわち、(49)式の基準温度Tstdと基準圧力Pstdと基準層流燃焼速度SLstdは実験により予め定められる値である。燃焼室5の通常の圧力である2bar以上の圧力下では、(49)式の圧力項(PTDC/Pstd)-0.16は小さな値となる。従って、圧力項(PTDC/Pstd)-0.16を一定値として、基準層流燃焼速度SLstdを基準温度Tstdのみで規定することも可能である。よって、基準温度Tstdが550[K]で、基準層流燃焼速度SLstdが1.0[m/sec]で、圧力項が0.7である場合の圧縮上死点における温度TTDCと層流燃焼速度SL2との関係は近似的に次式で定義することができる。
=1.0×0.7×(TTDC/550)2.18
…(50)
ステップ189では主燃焼期間におけるガス流動の乱れ強さST2を算出する。このガス流動の乱れ強さST2も初期燃焼期間におけるガス流動の乱れ強さST1と同様に、エンジン回転速度NRPMの関数として次式で表すことができる。
ただし、C2:定数、
乱れ強さST2を回転速度をパラメータとするテーブルから求めることも可能である。
ただし、SL2:層流燃焼速度[m/sec]、
(21)式は(18)式と同様、ガス乱れに伴う燃焼速度への寄与を考慮したものである。
/(RPROBA×AF2×FLAME2) …(22)
ただし、AF2:火炎核の反応面積[m2]、
ここで、(22)式右辺のBR2は主燃焼期間の開始時期より終了時期までの燃焼質量割合の変化代である。初期燃焼期間の終了時期に燃焼質量割合BRが2%になり、その後、主燃焼期間が開始し、燃焼質量割合BRが60%に達して主燃焼期間が終了すると考えているので、BR2=60%−2%=58%を設定している。AF2は火炎核の成長行程における平均の反応面積であり、(19)式のAF1と同様に、予め実験的に定めた固定値とする。
P0=PINI×Ec^κ …(10−2)
ただし、κ:比熱比、
ステップ416では燃焼室5の燃焼開始時期における圧力P0の関数として点火無駄時間基本値DEADTIME0[μsec]を次式により算出する。
例えば燃焼室5の燃焼開始時期における圧力P0から図44を内容とするテーブルを検索することにより点火無駄時間基本値DEADTIME0を求めればよい。点火無駄時間基本値DEADTIME0は理論空燃比のときかつ燃焼室5の燃焼開始時期における温度が基準温度である場合の値で、図44のように圧力P0が高いほど大きくなる。高圧で点火無駄時間が長くなるのは、高圧下では混合気が気化しにくくなるためである。
例えば燃焼室5の燃焼開始時期における温度T0から図45を内容とするテーブルを検索することにより温度補正係数Kdtmtを求めればよい。温度補正係数Kdtmtは図45のように基準温度以上のとき1.0であり、基準温度より低くなるほど大きくなる値である。すなわち、温度T0が低いほど点火無駄時間は長くなる。
ステップ419では、この目標空気過剰率TLAMBDAの関数として点火無駄時間の空気過剰率補正係数Kdtmlmbを次式により算出する。
例えば目標空気過剰率TLAMBDAから図46を内容とするテーブルを検索することにより空気過剰率補正係数Kdtmlmbを求めればよい。空気過剰率補正係数Kdtmlmbは図46のように目標空気過剰率TLAMBDAが1のとき(つまり理論空燃比のとき)最小の1.0となり、これより大きくても小さくても大きくなる値である。すなわち、空燃比が理論空燃比よりリッチ側にずれてもリーン側にずれても点火無駄時間は長くなる。
…(10−7)
ここでは、圧力P0によるものを基本値(DEADTIME0)として温度T0と空燃比(TLAMBDA)により補正する形式としたが、これに限られるものでなく、温度T0によるものを基本値として圧力P0と空燃比(TLAMBDA)により補正したり、空燃比(TLAMBDA)によるものを基本値として温度T0と圧力P0により補正するようにしてもかまわない。また、点火無駄時間に影響する因子はこれに限られるものでなく、充填効率や不活性ガス率(内部不活性ガス率MRESや外部不活性ガス率)があるので、これらの因子をも考慮して点火無駄時間を算出することができる。ただし、5つの因子全てを考慮する必要は必ずしもなく、少なくとも1つの因子をパラメータとして点火無駄時間を算出すればよい。
(10−8)式は、エンジン回転速度NRPMから点火無駄時間DEADTIMEに相当するクランク角である点火無駄時間相当クランク角IGNDEADを算出するためのものである。
ステップ44では、この基本点火時期MBTCALから点火無駄時間相当クランク角IGNDEADを差し引いた値を前回燃焼開始時期MBTCYCL[degBTDC]として算出する。
このようにして燃焼室5の排気弁閉時期EVCにおける不活性ガス量MRESCYLの算出を終了したら図15に戻り、ステップ62で吸排気弁15、16のオーバーラップ(図では「O/L」と略記する)中に排気側から吸気側へ吹き返す不活性ガス量であるオーバーラップ中吹き返し不活性ガス量MRESOLを算出する。
例えば、吸気VTC機構27用アクチュエータへの非通電時に吸気弁開時期IVOが吸気上死点位置にあり、吸気VTC機構27用アクチュエータへの通電時に吸気弁開時期が吸気上死点より進角する特性であり、かつ排気VTC機構28用アクチュエータへの非通電時に排気弁閉時期EVCが排気上死点にあり、排気弁VTC機構28用アクチュエータへの通電時に排気弁閉時期EVCが排気上死点より進角する特性である場合には、IVOとEVCの合計が吸排気弁のオーバーラップ量VTCOLとなる。
この吸気排気圧力比PINBYEXは無名数であり、これと1をステップ103で比較する。吸気排気圧力比PINBYEXが1以下の場合には過給無しと判断し、ステップ104に進んで過給判定フラグTBCRG(ゼロに初期設定)=0とする。
^{SHEATR/(SHEATR−1)} …(27a)
SLCHOKEH={−2/(SHEATR+1)}
^{−SHEATR/(SHEATR−1)}…(27b)
これらのチョーク判定しきい値SLCHOKEL、SLCHOKEHは、チョークする限界値を算出している。
〈2〉過給判定フラグTBCRG=0かつチョーク判定フラグCHOKE=0のとき
〈3〉過給判定フラグTBCRG=0かつチョーク判定フラグCHOKE=1のとき
〈4〉過給判定フラグTBCRG=1かつチョーク判定フラグCHOKE=0のとき
そして、上記〈1〉のときにはステップ89に進んで、過給無しかつチョーク無し時のオーバーラップ中の平均吹き返し不活性ガス流量MRESOLtmp1を、上記〈2〉のときにはステップ90に進んで過給無しかつチョーク有り時のオーバーラップ中の吹き返し不活性ガス流量MRESOLtmp2を、上記〈3〉のときにはステップ91に進んで過給有りかつチョーク無し時のオーバーラップ中の平均吹き返し不活性ガス流量MRESOLtmp3を、上記〈4〉のときにはステップ92に進んで過給有りかつチョーク有り時の吹き返し不活性ガス流量MRESOLtmp4をそれぞれ算出し、算出結果をオーバーラップ中の吹き返し不活性ガス流量MRESOLtmpに移す。
図19(図17ステップ89のサブルーチン)においてステップ121では、図16のステップ73、75で算出されている不活性ガスのガス定数REX、燃焼室5の排気弁閉時期における圧力PEVCを読み込む。
ここで、(28)式右辺の「SQRT」はすぐ右のカッコ内の値の平方根を計算させる関数である。
×{PTNBYEX^(2/SHEATR)
−PTNBYEX^((SHEATR+1)/SHEATR)}]…(29)
ステップ124では、これら密度項MRSOLD、圧力差項MRSOLPと、燃焼室5の排気弁閉時期における圧力PEVCとから、過給無しかつチョーク無し時のオーバーラップ中の吹き返し不活性ガス流量MRESOLtmp1を次式(ガス流量の算出式)により算出し、その算出値をステップ125でオーバーラップ中の吹き返し不活性ガス流量MRESOLtmpに移す。
次に、過給無しかつチョーク有り時の吹き返し不活性ガス流量の算出について図20のフローにより説明する
図20(図17ステップ90のサブルーチン)においてステップ131、132では、図19のステップ121、122と同様にして、不活性ガスのガス定数REX、燃焼室5の排気弁閉時期における圧力PEVCを読み込み、これらから前述の(28)式により密度項MRSOLDを算出する。
なお、(31)式の累乗計算と平方根計算とが困難な場合には、(31)式の算出結果を、チョーク時圧力差項MRSOLPCのテーブルとしてエンジンコントローラ31のメモリに予めに記憶しておき、不活性ガスの比熱比SHEATRからそのテーブルを検索することにより求めてもよい。
次に、過給有りかつチョーク無し時の吹き返しガス流量の算出について図21のフローにより説明する
図21(図17ステップ91のサブルーチン)においてステップ141では、吸気圧力センサ44により検出される吸気圧力PINを読み込む。
×{PINBYEX^(−2/SHEATR)
−PINBYEX^(−(SHEATR+1)/SHEATR)}]…(33)
なお、(33)式の累乗計算と平方根計算とが困難な場合は、(33)式の算出結果を、過給時圧力差項MRSOLPTのマップとしてエンジンコントローラ31のメモリに予め記憶しておき、不活性ガスの比熱比SHEATRと吸気排気圧力比PINBYEXとからそのマップを検索することにより求めてもよい。
ここで、(34)式の吹き返し不活性ガス流量MRESOLtmp3は負の値とすることで、オーバーラップ中に吸気系から排気系へ吹き抜ける混合気のガス流量を表すことができる。
図22(図17ステップ92のサブルーチン)においてステップ151、152では、図21のステップ141と同じく吸気圧力センサ44により検出される吸気圧力PINを読み込むと共に、図20のステップ132と同じくチョーク時圧力差項MRSOLPCを前述の(31)式により算出する。
ここで、(35)式の吹き返し不活性ガス流量MRESOLtmp4も、MRESOLtmp3と同様、負の値とすることで、オーバーラップ中に吸気側から排気側へ吹き抜ける混合気のガス流量を表すことができる。
/(NRPM×360) …(36)
このようにしてオーバーラップ中の吹き返し不活性ガス量MRESOLの算出を終了したら図15に戻り、ステップ63において燃焼室5内の排気弁閉時期EVCにおける不活性ガス量MRESCYLと、このオーバーラップ中吹き返しガス量MRESOLとを加算して、つまり次式により内部不活性ガス量MRESを算出する。
前述のように、過給有り時にはオーバーラップ中吹き返し不活性ガス流量(MRESOLtmp3、MRESOLtmp4)が負となるため、上記(36)式のオーバーラップ中の吹き返し不活性ガス量MRESOLも負となり、このとき(37)式によれば、オーバーラップ中の吹き返し不活性ガス量MRESOLの分だけ内部不活性ガス量が減じられる。
/{MRES+MACYL×(1+TFBYA/14.7)}…(38)
これで内部不活性ガス率MRESFRの算出を総て終了する。
ただし、CF#:燃料の低位発熱量、
一方、この発熱量Qにより燃焼室5内のガス温度が上昇するので、この温度上昇分をΔTとすれば次式が成立する。
ただし、M:燃焼室5内の全てのガスの質量、
Cv:既燃ガスの定容比熱、
(補4)、(補5)の両式は等しいと置いて温度上昇分ΔTについて解くと次式が得られる。
気体の状態方程式PV=nRTの両辺を微分する(ただし、定容変化であるためVは一定)。
ここで、ノックが発生するような状態ではモル数nの変化は小さいため、(補7)式右辺のdn=0として次式を得る。
(補8)、(補6)の両式より温度上昇分dT(=ΔT)を消去し圧力上昇分dPについて整理すると最終的に次式を得る。
すなわち、(補9)式は、未燃燃料量MUB、燃焼室5の自着火時期における容積V、既燃ガスの定容比熱Cv、燃焼室5内の全てのガスの質量M、燃焼室5内の全てのガスの総モル数nが分かれば、圧力上昇分dPを計算式により求めることができることを示している。
ただし、図31では計算を簡単にするため、燃焼開始遅れ期間(IGNDEADref1)、初期燃焼期間(BURN1ref1)、主燃焼期間(BURN2ref1)の3つ分け、各期間の特性を直線で近似している。
∴E=Cp×T …(補12)
(補12)式より定圧比熱Cpは次式で与えられる。
一方、理想気体で等圧変化のときにはCp−Cv=Rが成立するので、この式と(補13)式より定圧比熱Cpを消去し定容比熱Cvについて整理すると最終的に次式が得られる。
ただし、E:エンタルピ、
T:燃焼室5内の自着火時時期における平均温度、
上記(補9)式の燃焼室5内の全てのガスの質量Mは次式により計算できる。
ただし、MRES:内部不活性ガス量、
MACYL:シリンダ新気量、
QINJ:燃料量、
このように、未燃燃料量MUB、既燃ガスの定容比熱Cv、燃焼室5内の全てのガスの質量Mも(補10)、(補14)、(補15)の計算式によりそれぞれ求めることができることがわかる。残る未知数は、上記(補9)式の燃焼室5内の全てのガスの総モル数n、(補14)式のエンタルピE及び燃焼室5内の自着火時期における平均温度T(=TE)である。
/(WIDRY+QINJ+MASSZ) …(60)
ここでは、不活性ガスと新気の比熱を等しくして式を簡略化している。
×(1/τ100−1/τ80)/(100−80)
…(補17)
ただし、1/τEST:オクタン価推定値OCTESTの燃料での1/τ、
1/τ100:オクタン価100の燃料での1/τ、
1/τ80 :オクタン価80の燃料での1/τ、
ここで、オクタン価推定値OCTESTの算出については図37により後述する。
BR=0 …(61)
第1初期燃焼期間;
BR=SS1×(Θ+MBTCAL1−IGNDEADref1)
…(62)
第1主燃焼期間 ;
BR=0.02+SS2×(Θ+MBTCAL1−IGNDEADref1
−BURN1ref1)
…(63)
ただし、SS1:0.02/BURN1ref1、
SS2:0.58/BURN2ref1、
ここで述べた燃焼質量割合BRの設定(初期設定)は、後述する図39、図40のステップ308の内容であるため、第1初期燃焼期間は「BURN1ref1」、第1主燃焼期間は「BURN2ref1」、第1点火無駄時間相当クランク角は「IGNDEADref1」、第1基本点火時期は「MBTCAL1」としてあり、図10のステップ171で算出されているBURN1、図12のステップ191で算出されているBURN2、図51のステップ421で算出されている点火無駄時間相当値IGNDEAD、図13のステップ191で算出されている基本点火時期MBTCALと区別する。ただし、物理的な意味としては変わりない。
+CF#×QINJ×BR/(MASSZ+WIDRY+QINJ)
…(64)
PC=PC0×εθ^1.35×TC/TC0/εθ^0.35…(65)
ただし、εθ :瞬間圧縮比、
CF#:燃料の低位発熱量、
(64)、(65)式は燃焼室5内でガスが断熱圧縮されると共に定容変化で燃焼すると仮定したときの式である。すなわち、(64)式右辺第1項が断熱圧縮後の温度を、(65)式右辺のPC0×εθ^1.35が断熱圧縮後の圧力を、これに対して(64)式右辺第2項が定容変化で燃焼により温度上昇した分を、(65)式右辺のTC/TC0/εθ^0.35が定容変化での燃焼による圧力上昇率を表している。
^(0.35/1.35) …(66)
(66)式は燃焼室5内でガスが断熱圧縮されると共に、(64)式とは相違して、可逆断熱変化で燃焼すると仮定したときの式である。すなわち、(66)式右辺のTC0×εθ^0.35が断熱圧縮後の温度を、(66)式右辺の(PC/PC0/εθ^1.35)^(0.35/1.35)が可逆断熱変化での燃焼による温度上昇率を表している。
(67)式は上記(補10)式そのものである。
ステップ243では燃焼室5内の燃料が全て燃焼した(つまりBR=1)ときの各ガス成分のモル数を算出する。ただし、ガス成分としては、燃料のほかは、O2、N2、CO2、CO、H2Oに限定する。また、ガソリンの燃料組成をC7H14で近似している。
/(B#×AC#+A#×AH)×(A#/4) …(69.1)
酸素 ;XEO2 ={MIDRY#×WlDRY×0.21−QINJ
/(B#×AC#+A#×AH)×(B#+A#/4)}
/WEDRY …(69.2)
二酸化炭素;XECO2={QINJ/(B#×AC#+A#×AH#)×B#}
/WEDRY …(69.3)
一酸化炭素;XECO =0 …(69.4)
窒素 ;XEN2 =1−XEO2−XECO2−XECO …(69.5)
水 ;XEH2O={MIDRY#×WIDRY×15/745
+QINJ/(B#×AC#+A#×AH#)
×A#/2}/WEDRY …(69.6)
ただし、MIDRY#:1g当たりの新気ガスのモル数
AH#:水素のモル質量、
AC#:炭素のモル質量、
A#、B#:定数、
ここで、ガソリンの組成をC7H14で近似しているため、定数A#は14、定数B#は7である。
…(70.1)
不活性ガス;WEGR=MIDRY#×WIDRY×RTOEGR…(70.2)
酸素 ;WO2 =MIDRY#×WIDRY×0.21+WEGR×XEO2
…(70.3)
窒素 ;WN2 =MIDRY#×WIDRY×0.89+WEGR×XEN2
…(70.4)
二酸化炭素;WCO2=WEGR×XECO2 …(70.5)
一酸化炭素;WCO =WEGR×XECO …(70.6)
水 ;WH2O=MIDRY#×WIDRY×15/745
+WEGR×XEH2O …(70.7)
次に、全て燃焼した(つまりBR=1)ときの各ガス成分のモル数MLGAS[mol]、MLO2[mol]、MLN2[mol]、MLCO2[mol]、MLCO[mol]、MLH2O[mol]を次のように算出する。
…(71.1)
酸素 ;MLO2 =WO2−(B#+A#/4)×QlNJ
/(B#×AC#+A#×AH#) …(71.2)
窒素 ;MLN2 =WN2 …(71.3)
二酸化炭素;MLCO2=WCO2+B#×QINJ
/(B#×AC#+A#×AH#) …(71.4)
一酸化炭素;MLCO =WCO …(71.5)
水 ;MLH2O=WH2O+A#/2×QINJ
/(B#×AC#+A#×AH#) …(71.6)
これで、燃焼室5内の燃料が全て燃焼した(つまりBR=1)ときの各ガス成分のモル数の算出を終了するので、ステップ244に進み各ガス成分のモル数の総和を燃焼室5内の燃料が全て燃焼したときの総ガスモル数MLALLとして、つまり次式により総ガスモル数MLALLを算出する。
…(71.7)
このようにして総ガスモル数MLALLの算出を終わったら図33のステップ224に戻り、ガスエンタルピ(自着火する燃料ガスのエンタルピ)E[cal/mol]を算出する。このガスエンタルピの算出については図35のフローにより説明する。図35(図33ステップ224のサブルーチン)においてステップ251では、図33のステップ220で算出されている燃焼室5の自着火時平均温度TE、図34のステップ243、244で算出されている各ガス成分のモル数MLGAS、MLO2、MLN2、MLCO2、MLCO、MLH2O、総ガスモル数MLALLを読み込む。
(1)TE<1200Kの場合
E=A0#+1000×(A1#×(TE/1000)
+A2#/2×(TE/1000)^2
+A3#/3×(TE/1000)^3
+A4#/4×(TE/1000)^4
+A5#/5×(TE/1000)^5)+HDL#…(72.1)
(2)TE>1200Kの場合
E=B0#+1000×(B1#×(TE/1000)
+B2#×LN(TE/1000)
−B3#/(TE/1000)
−B4#/2/(TE/1000)^2
−B5#/3/(TE/1000)^3)+HDL#…(72.2)
ただし、A0#〜A5#、B0〜B5#、HDL#は実験により求める適合値、
ステップ253では燃料のエンタルピEGを次式により算出する。
+(B#/AC#+A#/AH#/4)×EO2…(72.3)
ステップ254では各ガス成分の平均エンタルピEを次式により計算して図35の処理を終了し、図33のステップ225に戻る。
+MLCO2×ECO2+MLCO×ECO+MLH2O×EH2O)
/MLALL …(72.4)
図33のステップ225では既燃ガスの定容比熱Cv[J/K・g]をガスエンタルピE、燃焼室5の自着火時平均温度TEを用いて次式により計算する。
ただし、R#:普遍気体定数、
(73)式は上記(補14)式においてT→TE、R→R#の置き換えにより得られる式である。
/{Cv×Vknk×(MASSZ+QINJ+WIDRY)}
…(74)
ただし、CF#:燃料の低位発熱量、
(74)式の圧力上昇分DPは図29に示したようにノック発生により燃焼室5内の圧力がステップ的に大きくなるので、この圧力上昇分を計算式により求めるようにしたものである。
ここで、(75)式右辺の相関係数1はノック強度との相関を表すための係数である。この場合、ノックによる圧力上昇分DPが大きいほどノック強度推定値基本値KIC0が大きくなるようにしている。
ここで、回転速度補正係数KNは、エンジン回転速度NRPMが低いときのほうが回転速度が高いときよりノックによる圧力振動をドライバーが強く感じるので、この違いを第1ノック強度推定値に反映させるためのものである。すなわち、図36のようにKNの値は、基準回転速度NRPM0のときを1.0としてこれより低い回転速度域では1.0を超える値、この逆に基準回転速度NRPM0より高い回転速度域で1.0未満の値である。実際の値は実験により適合する。
ただし、OCTEST(new):更新後のオクタン価推定値、
OCTEST(old):更新前のオクタン価推定値、
const03 :小さくする側への更新量、
ノックが検出されない場合にはステップ261よりステップ263に進み図2のステップ3において算出されている点火時期最小値PADV[degBTDC]と図2のステップ1において算出されている基本点火時期MBTCAL[degBTDC]を比較する。点火時期最小値PADVが基本点火時期MBTCALと一致するときにはオクタン価推定値OCTESTと実際のオクタン価が一致しており、従ってオクタン価推定値を変更する必要がないので、そのまま今回の処理を終了する。
ただし、OCTEST(new):更新後のオクタン価推定値、
OCTEST(old):更新前のオクタン価推定値、
const05 :大きくする側への更新量、
オクタン価推定値OCTESTの更新は、カウンタ値countが所定値const04となる毎であるので、ステップ267ではカウンタ値countをゼロにリセットする。
…(81)
Vmb31p=f25(θign−IGNDEAD−BURN1
−BURN2/2)
…(82)
(81)式右辺の変数としてのθign−IGNDEAD−BURN1×(1/2)は燃焼質量割合が1%のときのクランク角[degBTDC]、(82)式右辺の変数としてのθign−IGNDEAD−BURN1−BURN2×(1/2)は燃焼質量割合が31%のときのクランク角[degBTDC]で、(81)、(82)式はこれら各クランク角の関数として燃焼室5のその各クランク角における容積Vmb1p、Vmb31pを求めるものである。すなわち、図31に示したように、燃焼質量割合とクランク角Θとの間には関係があり、また、図5のステップ12、図10のステップ162で求めたように燃焼室5の容積はそのときのクランク角により定まる(クランク角の関数である)ので、燃焼質量割合が1%、31%のときの各クランク角Θ[degATDC]を図31に示す燃焼質量割合の特性と同様の燃焼質量割合の特性(予め設定してある)を用いて求め、その求めたクランク角Θをクランク角θ[degBTDC]に変換し、その変換したクランク角θから、クランク角θをパラメータとするテーブルを検索することにより燃焼室5の燃焼質量割合1%時における容積Vmb1pと、燃焼室5の燃焼質量割合31%時における容積Vmb31pとを求めることができる。
ステップ343〜347によりTCとPCを算出するための操作は図32のステップ202、204、205、215、216によりTCとPCを算出するための操作と同様である。異なるのは図32のステップ215、216においてはそのときのクランク角θでの圧縮比であるεθと、そのときのクランク角θでの燃焼質量割合であるBRとを用いたのに対して、図41のステップ346、347では燃焼質量割合1%時の瞬間圧縮比であるEmb1pと、燃焼質量割合1%時の燃焼質量割合であるBR1pとを用いる点である。すなわち、図41のステップ346、347では燃焼室5内の燃料が燃焼したときの平均温度TC[K]と平均圧力PC[Pa]を次式により算出する。
+CF#×QINJ×BR1p/(MASSZ+WIDRY+QINJ)
…(84)
PC=PC0×Emb1p^1.35×TC/TC0/Emb1p^0.35
…(85)
ただし、CF#:燃料の低位発熱量、
ステップ348〜351により第1初期燃焼期間を算出するための操作は図10のステップ168〜171により初期燃焼期間を算出するための操作と同様である。異なるのは図10のステップ168において燃焼室5の燃焼開始時期における温度T0と圧力P0とを用いて層流燃焼速度SL1を算出したのに対して、図41のステップ348では上記(84)式、(85)式で求めた温度TCと圧力PCを用いて層流燃焼速度SL1を算出する点、また図10のステップ171において初期燃焼期間BURN1の終了時期における燃焼質量割合であるBR1(=2%)と、燃焼室5の燃焼開始時期における容積であるV0とを用いて初期燃焼期間BURN1を算出したのに対して、図41のステップ351では第1初期燃焼期間の中間時期における燃焼質量割合である1%(=BR1p)と、燃焼室5の燃焼質量割合1%時における容積であるVmb1pとを用いて第1初期燃焼期間BURN1ref1を算出する点である。すなわち、ステップ348では第1初期燃焼期間における層流燃焼速度SL1[m/sec]を、またステップ351では第1初期燃焼期間BURN1ref1[deg]をそれぞれ次式により算出する。
ただし、Tstd :基準温度[K]、
Pstd :基準圧力[Pa]、
SLstd:基準温度Tstdと基準圧力Pstdにおける基準層流燃焼
速度[m/sec]、
BURN1ref1={(NRPM×6)×BR1p×Vmb1p}
/(RPROBA×AF1×FLAME1) …(87)
ただし、AF1 :火炎核の反応面積(固定値)[m2]、
なお、図41のステップ351に示す「BURN1refi」は、第1初期燃焼期間である「BURN1ref1」と第2初期燃焼期間である「BURN1ref2」とをまとめて示している(つまりiは1または2)。これは、図41のフローが図39のステップ305と図40のステップ321の2箇所で用いられ、このうち図39のステップ305では第1初期燃焼期間BURN1ref1を、図40のステップ321では第2初期燃焼期間BURN1ref2を算出するようにしているためである。
ステップ363〜367によりTCとPCを算出するための操作は図32のステップ202、204、205、215、216によりTCとPCを算出するための操作と同様である。異なるのは図32のステップ215、216においてはそのときのクランク角θでの圧縮比であるεθと、そのときのクランク角θでの燃焼質量割合であるBRとを用いたのに対して、図42のステップ366、367では燃焼質量割合31%時の圧縮比であるEmb31pと、燃焼質量割合31%時の燃焼質量割合であるBR31pとを用いる点である。すなわち、図42のステップ366、367では燃焼室5内の燃料が燃焼したときの平均温度TC[K]と平均圧力PC[Pa]を次式により算出する。
+CF#×QINJ×BR31p/(MASSZ+WIDRY+QINJ)
…(89)
PC=PC0×Emb31p^1.35×TC/TC0/Emb31p^0.35
…(90)
ただし、CF#:燃料の低位発熱量、
ステップ368〜371により第1主燃焼期間を算出するための操作は図12のステップ188〜191により主燃焼期間を算出するための操作と同様である。異なるのは図12のステップ188において燃焼室5の圧縮上死点における温度TTDCと圧力PTDCとを用いたのに対して、図42のステップ368では上記(89)式、(90)式で求めた温度TCと圧力PCとを用いる点、また図12のステップ191において主燃焼期間BURN2の終了時期における燃焼質量割合であるBR2(=58%)と、燃焼室5の圧縮上死点における容積であるVTDCとを用いたのに対して、図42のステップ371では第1主燃焼期間の中間時期の近傍における燃焼質量割合である31%(=BR31p)と、燃焼室5の燃焼質量割合31%時における容積であるVmb31pとを用いる点である。すなわち、ステップ368では第1主燃焼期間における層流燃焼速度SL2[m/sec]を、またステップ371では第1主燃焼期間BURN2ref1[deg]をそれぞれ次式により算出する。
ただし、Tstd :基準温度[K]、
Pstd :基準圧力[Pa]、
SLstd:基準温度Tstdと基準圧力Pstdにおける基準層流燃焼
速度[m/sec]、
BURN2ref1={(NRPM×6)×BR31p×Vmb31p}
/(RPROBA×AF2×FLAME1) …(92)
ただし、AF2 :火炎核の反応面積(固定値)[m2]、
なお、図42のステップ371に示す「BURN2refi」も、第1主燃焼期間である「BURN2ref1」と第2主燃焼期間である「BURN2ref2」とをまとめて示している(つまりiは1または2)。これは、図42のフローが図39のステップ305と図40のステップ321の2箇所で用いられ、このうち図39のステップ305では第1主燃焼期間BURN2ref1を、図40のステップ321では第2主燃焼期間BURN2ref2を算出するようにしているためである。
+IGNDEADref1
…(93)
ここでは図13のフローにより算出されている基本点火時期MBTCALと区別するため、添字としての「1」をつけている。従ってこの(93)式は上記(23)式と基本的に同じである。
図40のステップ314〜327は、第2基本点火時期MBTCAL2より遅角した点火時期で点火して燃焼すると仮定した場合のノック強度推定推定値を第2ノック強度推定値KIC2として算出する部分である。図40においてまずステップ314では点火時期θign[degBTDC](=MBTCAL1)の値を前回値を表すθignz[degBTDC]に移し、この点火時期前回値θignzからノックリタード量KNRTを差し引いた値(つまり第1基本点火時期MBTCAL1よりノックリタード量KNRTだけ遅角側の値)を点火時期θign[degBTDC]として算出する。これは初回だけであり、2回目からは点火時期θignの値を前回値を表すθignzに移し、この点火時期前回値θignzからノックリタード量KNRTを差し引いた値(つまり図40のステップ323で算出される第2基本点火時期MBTCAL2よりノックリタード量KNRTだけ遅角側の値)を点火時期θignとして算出する。
−BURN1ref2/2)
…(95)
ステップ317〜320は燃焼室5の燃焼質量割合31%時における容積Vmb31pを算出する部分である。まずステップ317では燃焼質量割合が31%のときのクランク角の前回値θmb31pz[degBTDC]を次式により算出する。
−BURNN2ref2/2
…(96)
ステップ318では燃焼質量割合が2%のときのクランク角の前回値と、燃焼質量割合が2%のときのクランク角の今回値との差θ2pdif[deg]を次式により算出する。
−(θign−IGNDEADref2−BURN1ref2) …(97)
ステップ319では燃焼質量割合が31%のときのクランク角の今回値θmb31p[degBTDC]を次式により算出する。
ステップ320ではこのクランク角θmb31pから燃焼室5の燃焼質量割合31%時における容積Vmb31p[m3]を次式により算出する。この操作は図39のステップ304の操作と同様である。
ステップ321、322では上記(95)式、(99)式で得られた燃焼室容積Vmb1p、Vmb31pを用いて第2初期燃焼期間BURN1ref2、第2主燃焼期間BURN2ref2を、また初回は第1基本点火時期MBTCAL1より遅角された点火時期(MBTCAL1−KNRT)を用いて燃焼室5の燃焼開始時期における温度T0及び圧力P0を算出し第2点火無駄時間相当クランク角IGNDEADref2を算出する。これら第2燃焼期間BURN1ref2、BURN2ref2及び第2点火無駄時間相当クランク角IGNDEADref2の算出には図41、図42、図43のフローを用いればよい。
+IGNDEADref2
…(100)
これで、燃焼質量波形BRの再設定に必要となる第2基本点火時期MBTCAL2、第2初期燃焼期間BURN1ref2、第2主燃焼期間BURN2ref2、第2点火無駄時間相当クランク角IGNDEADref2が全て求まったので、ステップ324ではこれらの値を用いて燃焼質量波形を再設定する。すなわち、上記(61)式〜(63)式と同様の次式により第2燃焼遅れ期間、第2初期燃焼期間、第2主燃焼期間における燃焼質量割合BRの一次式を算出し、算出した式は記憶しておく。
BR=0 …(101)
第2初期燃焼期間;
BR=SS1×(Θ+MBTCAL2−IGNDEADref2)
…(102)
第2主燃焼期間 ;
BR=0.02+SS2×(Θ+MBTCAL2−IGNDEADref2
−BURN1ref2)
…(103)
ただし、SS1:0.02/BURN1ref2、
SS2:0.58/BURN2ref2、
図48は初期設定の燃焼質量割合の波形(実線参照)と、再設定される燃焼質量割合の波形(一点鎖線参照)との違いをモデルで示したものである。第1基本点火時期MBTCAL1よりノックリタード量KNRTだけ遅角させた点火時期(MBTCAL1−KNRT)で点火して燃焼する場合には第1基本点火時期MBTCAL1で点火して燃焼する場合よりも燃焼速度が遅くなるため、そのぶん燃焼質量割合の波形が変化する。
そこで、改めて、第1基本点火時期MBTCAL1よりノックリタード量KNRTだけ遅角させた点火時期(MBTCAL1−KNRT)で点火して燃焼する場合の値である第2燃焼期間BURN1ref2、BURN2ref2、第2点火無駄時間相当クランク角IGNDEADref2、第2基本点火時期MBTCAL2を算出し、これらより燃焼質量割合を再設定したのが一点鎖線の特性である。この一点鎖線の特性によれば、基本点火時期からの遅角に伴う燃焼速度の遅れの分だけ変化した実際の燃焼質量割合の変化をよく近似できるものとなっている。
ただし、KNRT(new):更新後のノックリタード量、
KNRT(old):更新前のノックリタード量、
const1 :大きくする側への更新量、
このようにしてノックリタード量を更新した後はステップ314に戻ってステップ314〜326の操作を再び実行する。すなわち、ステップ314では更新後のノックリタード量KNRTを用いて点火時期θignを算出し、この新たな点火時期θignを用いてステップ315〜323の操作を実行して、燃焼質量波形の設定に必要となるMBTCAL2、BURN1ref2、BURN2ref2、IGNDEADref2を算出しこれらを用いステップ324において燃焼質量波形を再設定し、その再設定した燃焼質量割合の波形を用いてその新たな点火時期θignでの第2ノック強度推定値KIC2をステップ325で算出し、その算出した第2ノック強度推定値KIC2とトレースノック強度をステップ326で比較する。比較の結果、第2ノック強度推定値KIC2がトレースノック強度と一致すればステップ326よりステップ327に進んでそのときの点火時期であるMBTCAL2−KNRT(=θign)の値をノック限界点火時期KNOCKcalに移して今回の処理を終了する。
ただし、KNRT(new):更新後のノックリタード量、
KNRT(old):更新前のノックリタード量、
const2 :小さくする側への更新量、
このようにしてノックリタード量を更新した後はステップ314に戻ってステップ314〜326の操作を再び実行する。すなわち、ステップ314では更新後のノックリタード量KNRTを用いて点火時期θignを算出し、この新たな点火時期θignを用いてステップ315〜323の操作を実行して、燃焼質量波形の設定に必要となるMBTCAL2、BURN1ref2、BURN2ref2、IGNDEADref2を算出しこれらを用いステップ324において燃焼質量波形を再設定し、その再設定した燃焼質量割合の波形を用いてその新たな点火時期θignでの第2ノック強度推定値KIC2をステップ325で算出し、その算出した第2ノック強度推定値KIC2とトレースノック強度をステップ326で比較する。比較の結果、第2ノック強度推定値KIC2がトレースノック強度と一致すればステップ326よりステップ327に進んでそのときの点火時期であるMBTCAL2−KNRT(=θign)の値をノック限界点火時期KNOCKcalに移して今回の処理を終了する。
5 燃焼室
11 点火装置(火花点火手段)
21 燃料インジェクタ
31 エンジンコントローラ
Claims (15)
- 所定の燃焼期間に基づいて得られる基本点火時期を第1基本点火時期として算出する第1基本点火時期算出手段と、
ノック限界点火時期を設定するノック限界点火時期設定手段と、
これら第1基本点火時期、ノック限界点火時期のいずれかのうち遅角側の値で火花点火を行う火花点火手段と
を備えるエンジンのノック制御装置において、
前記ノック限界点火時期設定手段が、
前記第1基本点火時期で点火した場合に燃焼室内の燃料が自着火する時期を推定する自着火時期推定手段と、
この算出した自着火時期に基づいて第1ノック強度推定値を算出する第1ノック強度推定値算出手段と、
この第1ノック強度推定値とスライスレベルを比較する比較手段と、
この比較結果より第1ノック強度推定値がスライスレベル未満である場合に前記第1基本点火時期をそのままノック限界点火時期として設定するノック限界点火時期設定手段と、
前記比較結果より第1ノック強度推定値がスライスレベル以上である場合に所定のノックリタード量を算出するノックリタード量算出手段と、
前記第1基本点火時期よりこのノックリタード量だけ遅角させた点火時期で点火して燃焼する場合の燃焼開始から所定クランク角までの燃焼期間を算出する燃焼期間算出手段と、
この燃焼期間に基づいて得られる基本点火時期を第2基本点火時期として算出する第2基本点火時期算出手段と、
この第2基本点火時期より前記ノックリタード量だけ遅角させた点火時期で点火して燃焼する場合のノック強度を第2ノック強度推定値として算出する第2ノック強度推定値算出手段と、
この第2ノック強度推定値が前記スライスレベルと一致するとき前記第2基本点火時期より前記ノックリタード量だけ遅角させた点火時期をノック限界点火時期として設定するノック限界点火時期設定手段と、
前記第2ノック強度推定値が前記スライスレベルと一致せず第2ノック強度推定値がスライスレベルより大きい場合には、
前記算出したノックリタード量を所定値だけ大きくなる側に更新し、
前記第1基本点火時期よりこの更新後のノックリタード量だけ遅角させた点火時期で点火して燃焼する場合の燃焼開始から所定クランク角までの燃焼期間を更新後の燃焼期間として算出し、
この更新後の燃焼期間に基づいて得られる基本点火時期を更新後の第2基本点火時期として算出し、
この更新後の第2基本点火時期より前記更新後のノックリタード量だけ遅角させた点火時期で点火して燃焼する場合のノック強度を更新後の第2ノック強度推定値として算出し、
この更新後の第2ノック強度推定値と前記スライスレベルとを比較する
ことを繰り返し、この繰り返しにより更新後の第2ノック強度推定値がスライスレベルと一致するとき、更新後の第2基本点火時期より更新後のノックリタード量だけ遅角させた点火時期をノック限界点火時期として、
また前記第2ノック強度推定値が前記スライスレベルと一致せず第2ノック強度推定値がスライスレベルより小さい場合には、
前記算出したノックリタード量を所定値だけ小さくなる側に更新し、
前記第1基本点火時期よりこの更新後のノックリタード量だけ遅角させた点火時期で点火して燃焼する場合の燃焼開始から所定クランク角までの燃焼期間を更新後の燃焼期間として算出し、
この更新後の燃焼期間に基づいて得られる基本点火時期を更新後の第2基本点火時期として算出し、
この更新後の第2基本点火時期より前記更新後のノックリタード量だけ遅角させた点火時期で点火して燃焼する場合のノック強度を更新後の第2ノック強度推定値として算出し、
この更新後の第2ノック強度推定値と前記スライスレベルとを比較する
ことを繰り返し、この繰り返しにより更新後の第2ノック強度推定値がスライスレベルと一致するとき、更新後の第2基本点火時期より更新後のノックリタード量だけ遅角させた点火時期をノック限界点火時期としてそれぞれ設定するノック限界点火時期設定手段と
を備えることを特徴とするエンジンのノック制御装置。 - 前記第1ノック強度推定値算出手段は、前記第1基本点火時期で点火して燃焼する場合の燃焼質量割合の特性に基づいて前記自着火時期における燃焼質量割合を算出する自着火時燃焼質量割合算出手段と、この算出した自着火時期における燃焼質量割合と燃料量とに基づいて第1ノック強度推定値を算出するノック強度推定値算出手段とからなり、
前記第2ノック強度推定値算出手段は、前記燃焼期間に基づいて燃焼質量割合の特性を、前記第2基本点火時期より前記ノックリタード量だけ遅角させた点火時期で点火して燃焼する場合の燃焼質量割合の特性へと再設定する再設定手段と、この再設定後の燃焼質量割合の特性を用いて前記第2基本点火時期より前記ノックリタード量だけ遅角させた点火時期で点火して燃焼する場合のノック強度を第2ノック強度推定値として算出するノック強度推定値算出手段とからなることを特徴とする請求項1に記載のエンジンのノック制御装置。 - 前記自着火時期推定手段は、燃焼室内の燃料が自着火する時期を燃焼室内の燃料が自着火にいたるまでの時間の逆数の分布を表す特性に基づいて推定すると共に、
前記第1ノック強度推定値算出手段は、前記第1基本点火時期で点火して燃焼する場合の燃焼質量割合の特性に基づいて前記自着火時期における燃焼質量割合を算出する自着火時燃焼質量割合算出手段と、この算出した自着火時期における燃焼質量割合と燃料量とに基づいて未燃燃料量または未燃燃料割合を算出する未燃燃料量算出手段と、この未燃燃料量または未燃燃料割合に基づいて燃焼室内のノックによる圧力上昇量を算出する圧力上昇量算出手段と、この燃焼室内のノックによる圧力上昇量に基づいて第1ノック強度推定値を算出するノック強度推定値算出手段とからなり、
前記第2ノック強度推定値算出手段は、前記燃焼期間に基づいて燃焼質量割合の特性を、前記第2基本点火時期より前記ノックリタード量だけ遅角させた点火時期で点火して燃焼する場合の燃焼質量割合の特性へと再設定する再設定手段と、この再設定後の燃焼質量割合の特性を用いて前記第2基本点火時期より前記ノックリタード量だけ遅角させた点火時期で点火して燃焼する場合のノック強度を第2ノック強度推定値として算出するノック強度推定値算出手段とからなることを特徴とする請求項1に記載のエンジンのノック制御装置。 - 燃焼ガスの層流状態での燃焼速度である層流燃焼速度を算出する層流燃焼速度算出手段と、
燃焼室内の燃焼ガス体積に相当する容積を算出する燃焼ガス体積相当容積算出手段と、
所定のクランク角までに前記燃焼室内で燃焼するガスの燃焼質量割合を算出する燃焼質量割合算出手段と、
所定運転条件での燃焼ガスの燃焼のしやすさを示す反応確率を算出する反応確率算出手段と、
これら層流燃焼速度、燃焼ガス体積相当容積、燃焼質量割合及び反応確率に基づいて燃焼開始から所定クランク角までの燃焼期間を第1燃焼期間として算出する第1燃焼期間算出手段と、
この第1燃焼期間に基づいて得られる基本点火時期を第1基本点火時期として算出する第1基本点火時期算出手段と、
ノック限界点火時期を設定するノック限界点火時期設定手段と、
これら第1基本点火時期、ノック限界点火時期のいずれかのうち遅角側の値で火花点火を行う火花点火手段と
を備えるエンジンのノック制御装置において、
前記ノック限界点火時期設定手段が、
燃焼室内の燃料が自着火にいたるまでの時間の逆数の分布を表す特性に基づいて燃焼室内の燃料が自着火する時期を推定する自着火時期推定手段と、
前記第1基本点火時期で点火して燃焼するときの燃焼質量割合の特性に基づいて前記自着火時期における燃焼質量割合を算出する自着火時燃焼質量割合算出手段と、
この算出した自着火時期における燃焼質量割合と燃料量とに基づいて未燃燃料量または未燃燃料割合を算出する未燃燃料量算出手段と、
この未燃燃料量または未燃燃料割合に基づいて燃焼室内のノックによる圧力上昇量を算出する圧力上昇量算出手段と、
この燃焼室内のノックによる圧力上昇量に基づいて第1ノック強度推定値を算出する第1ノック強度推定値算出手段と、
この第1ノック強度推定値とスライスレベルを比較する比較手段と、
この比較結果より第1ノック強度推定値がスライスレベル未満である場合に前記第1基本点火時期をそのままノック限界点火時期として設定するノック限界点火時期設定手段と、
前記比較結果より第1ノック強度推定値がスライスレベル以上である場合に所定のノックリタード量を算出するノックリタード量算出手段と、
前記第1基本点火時期よりこのノックリタード量だけ遅角させた点火時期で点火して燃焼すると仮定した場合の燃焼開始から所定クランク角までの燃焼期間を第2燃焼期間として算出する第2燃焼期間算出手段と、
この第2燃焼期間に基づいて得られる基本点火時期を第2基本点火時期として算出する第2基本点火時期算出手段と、
前記第2燃焼期間に基づいて燃焼質量割合の特性を、前記第2基本点火時期より前記ノックリタード量だけ遅角させた点火時期で点火して燃焼する場合の燃焼質量割合の特性へと再設定する再設定手段と、
この再設定後の燃焼質量割合の特性を用いて前記第2基本点火時期より前記ノックリタード量だけ遅角させた点火時期で点火して燃焼すると仮定した場合のノック強度を第2ノック強度推定値として算出する第2ノック強度推定値算出手段と、
この第2ノック強度推定値が前記スライスレベルと一致するとき前記第2基本点火時期より前記ノックリタード量だけ遅角させた点火時期をノック限界点火時期として設定するノック限界点火時期設定手段と、
前記第2ノック強度推定値が前記スライスレベルと一致せず第2ノック強度推定値がスライスレベルより大きい場合には、
前記算出したノックリタード量を所定値だけ大きくなる側に更新し、
前記第1基本点火時期よりこの更新後のノックリタード量だけ遅角させた点火時期で点火して燃焼すると仮定した場合の燃焼開始から所定クランク角までの燃焼期間を更新後の第2燃焼期間として算出し、
この更新後の第2燃焼期間に基づいて得られる基本点火時期を更新後の第2基本点火時期として算出し、
前記更新後の第2燃焼期間に基づいて燃焼質量割合の特性を、前記更新後の第2基本点火時期より前記更新後のノックリタード量だけ遅角させた点火時期で点火して燃焼する場合の燃焼質量割合の特性へと再設定し、
この再設定後の燃焼質量割合の特性を用いて前記更新後の第2基本点火時期より前記更新後のノックリタード量だけ遅角させた点火時期で点火して燃焼すると仮定した場合のノック強度を更新後の第2ノック強度推定値として算出し、
この更新後の第2ノック強度推定値と前記スライスレベルとを比較する
ことを繰り返し、この繰り返しにより更新後の第2ノック強度推定値がスライスレベルと一致するとき、更新後の第2基本点火時期より更新後のノックリタード量だけ遅角させた点火時期をノック限界点火時期として、
また前記第2ノック強度推定値が前記スライスレベルと一致せず第2ノック強度推定値がスライスレベルより小さい場合には、
前記算出したノックリタード量を所定値だけ小さくなる側に更新し、
前記第1基本点火時期よりこの更新後のノックリタード量だけ遅角させた点火時期で点火して燃焼すると仮定した場合の燃焼開始から所定クランク角までの燃焼期間を更新後の第2燃焼期間として算出し、
この更新後の第2燃焼期間に基づいて得られる基本点火時期を更新後の第2基本点火時期として算出し、
前記更新後の第2燃焼期間に基づいて燃焼質量割合の特性を、前記更新後の第2基本点火時期より前記更新後のノックリタード量だけ遅角させた点火時期で点火して燃焼する場合の燃焼質量割合の特性へと再設定し、
この再設定後の燃焼質量割合の特性を用いて前記更新後の第2基本点火時期より前記更新後のノックリタード量だけ遅角させた点火時期で点火して燃焼すると仮定した場合のノック強度を更新後の第2ノック強度推定値として算出し、
この更新後の第2ノック強度推定値と前記スライスレベルとを比較する
ことを繰り返し、この繰り返しにより更新後の第2ノック強度推定値がスライスレベルと一致するとき、更新後の第2基本点火時期より更新後のノックリタード量だけ遅角させた点火時期をノック限界点火時期としてそれぞれ設定するノック限界点火時期設定手段と
を備えることを特徴とするエンジンのノック制御装置。 - 前記第2燃焼期間算出手段が、前記第2燃焼期間を、前記第1基本点火時期で点火して燃焼する場合の層流燃焼速度と異なる層流燃焼速度を用いて算出することを特徴とする請求項4に記載のエンジンのノック制御装置。
- 前記第2燃焼期間算出手段が、前記第2燃焼期間を、前記第1基本点火時期で点火して燃焼する場合の層流燃焼速度と同じ層流燃焼速度を用いて算出することを特徴とする請求項4に記載のエンジンのノック制御装置。
- 前記第1、第2の燃焼期間算出手段が前記層流燃焼速度を乱れ強さで補正すると共に、
前記第2燃焼期間算出手段が用いる乱れ強さは前記第1燃焼期間算出手段が用いる乱れ強さと同じであることを特徴とする請求項5または6に記載のエンジンのノック制御装置。 - 前記第1及び第2の燃焼期間を前記燃焼質量割合の変化に応じて複数に分割することを特徴とする請求項4に記載のエンジンのノック制御装置。
- 前記第2燃焼期間算出手段が、前記第2燃焼期間を算出する際に、燃焼室のその燃焼期間の中間時期における容積及び燃焼室のその燃焼期間の中間時期における燃焼質量割合を用いることを特徴とする請求項4に記載のエンジンのノック制御装置。
- 点火信号を前記火花点火手段に出力してからこの火花点火手段により燃焼室内で点火が行われるまでのクランク角区間である点火無駄時間相当クランク角を燃焼に関するパラメータに応じて算出する点火無駄時間相当クランク角算出手段を備え、前記再設定手段がこの点火無駄時間相当クランク角に基づいても前記燃焼質量割合の特性を再設定することを特徴とする請求項4に記載のエンジンのノック制御装置。
- 前記燃焼に関するパラメータは、空燃比、燃焼室の燃焼開始時期における圧力、燃焼室の燃焼開始時期における温度、充填効率、不活性ガス率の少なくとも一つであることを特徴とする請求項10に記載のエンジンのノック制御装置。
- 前記所定のノックリタード量は前記第1ノック強度推定値と前記スライスレベルの差に基づく値であることを特徴とする請求項1から4までのいずれか一つに記載のエンジンのノック制御装置。
- 前記自着火時期に基づいて燃焼室の自着火時期における容積を算出する自着火時容積算出手段を備え、
この燃焼室の自着火時期における容積に基づいても前記ノックによる圧力上昇量を算出することを特徴とする請求項4に記載のエンジンのノック制御装置。 - 燃焼室の前記自着火する燃料ガスの比熱を推定する燃料ガス比熱推定手段を備え、この燃焼室の前記自着火する燃料ガスの比熱に基づいても前記ノックによる圧力上昇量を算出することを特徴とする請求項4に記載のエンジンのノック制御装置。
- 前記燃料ガス比熱算出手段は、前記自着火する燃料ガスのエンタルピを計算式により算出する燃料ガスエンタルピ算出手段と、燃焼室の前記自着火時期における平均温度を計算式により算出する自着火時平均温度算出手段と、これら燃料ガスのエンタルピと自着火時期における平均温度とに基づいて前記自着火する燃料ガスの比熱を算出する燃料ガス比熱算出手段とからなることを特徴とする請求項14に記載のエンジンのノック制御装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2004131299A JP4151605B2 (ja) | 2004-04-27 | 2004-04-27 | エンジンのノック制御装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2004131299A JP4151605B2 (ja) | 2004-04-27 | 2004-04-27 | エンジンのノック制御装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2005315099A JP2005315099A (ja) | 2005-11-10 |
| JP4151605B2 true JP4151605B2 (ja) | 2008-09-17 |
Family
ID=35442800
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2004131299A Expired - Fee Related JP4151605B2 (ja) | 2004-04-27 | 2004-04-27 | エンジンのノック制御装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP4151605B2 (ja) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2014010022A1 (ja) * | 2012-07-09 | 2014-01-16 | トヨタ自動車株式会社 | 内燃機関の点火時期制御装置 |
| WO2019049878A1 (ja) | 2017-09-06 | 2019-03-14 | 株式会社Ihi | エンジン制御システム |
-
2004
- 2004-04-27 JP JP2004131299A patent/JP4151605B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JP2005315099A (ja) | 2005-11-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US7255089B2 (en) | Engine control device and control method | |
| JP3975936B2 (ja) | ノッキング指標値算出装置 | |
| JP4241511B2 (ja) | エンジンのノック制御装置 | |
| JP4120620B2 (ja) | エンジンの制御装置 | |
| JP4251069B2 (ja) | エンジンのノック検出装置及びノック制御装置 | |
| JP2005233110A (ja) | 内燃機関の点火時期制御装置 | |
| JP4774893B2 (ja) | エンジンの点火時期制御方法及びエンジンの点火時期制御装置 | |
| JP4151605B2 (ja) | エンジンのノック制御装置 | |
| JP4254605B2 (ja) | エンジンの点火時期制御装置 | |
| JP4241510B2 (ja) | エンジンのノック制御装置及びエンジンの点火時期制御装置。 | |
| JP4604903B2 (ja) | エンジンの点火時期制御装置 | |
| JP2005226481A (ja) | エンジンのデポジット量検出装置及びノック制御装置 | |
| JP4075862B2 (ja) | エンジンの点火時期制御装置 | |
| JP2008196506A (ja) | エンジンのノック制御装置及びエンジンの点火時期制御装置 | |
| JP4055648B2 (ja) | 内燃機関の点火時期制御装置 | |
| JP4055647B2 (ja) | エンジンの点火時期制御装置 | |
| JP4135655B2 (ja) | 内燃機関の点火時期制御装置 | |
| JP4239607B2 (ja) | エンジンの点火制御装置 | |
| JP4244654B2 (ja) | エンジンの点火時期制御装置 | |
| JP4182914B2 (ja) | 内燃機関の点火時期制御装置 | |
| JP4055632B2 (ja) | 内燃機関の点火時期制御装置 | |
| JP4055646B2 (ja) | エンジンの点火時期制御装置 | |
| JP4200777B2 (ja) | エンジンの制御装置 | |
| JP4075819B2 (ja) | 内燃機関の点火時期制御装置 | |
| JP4747763B2 (ja) | エンジンの点火時期制御方法及びエンジンの点火時期制御装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20060224 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20071206 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20080205 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20080401 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20080610 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20080623 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110711 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120711 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120711 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130711 Year of fee payment: 5 |
|
| LAPS | Cancellation because of no payment of annual fees |