JP4182897B2

JP4182897B2 - 内燃機関の点火時期制御装置

Info

Publication number: JP4182897B2
Application number: JP2004044495A
Authority: JP
Inventors: 太容吉野
Original assignee: Nissan Motor Co Ltd
Current assignee: Nissan Motor Co Ltd
Priority date: 2004-02-20
Filing date: 2004-02-20
Publication date: 2008-11-19
Anticipated expiration: 2024-02-20
Also published as: JP2005233114A

Description

本発明は、ノックを発生しないように点火時期を制御する内燃機関の点火時期制御装置に関するものである。

エンジンの点火時期は燃費や出力に影響し、また点火時期がずれるとノッキングや失火を生じるおそれがあるので、その制御は非常に重要である。通常は、１０〜１５degATDCでシリンダ内圧力が最大圧になるように点火時期を制御している。しかし、点火時期設定の基本となる最大トルク最小進角点火時期(Minimmum Spark Advance for Best Torque；以下「ＭＢＴ」と略す)は、エンジン回転速度、負荷、空燃比、ＥＧＲ率等によって時々刻々と変動する。そこで従来は、エンジン回転速度や負荷に応じて設定した基本点火時期マップと、さまざまな運転状態を想定した補正値マップとを備え、運転状態に応じてそれらのマップを参照して制御していた。しかし、そのような方法では、精度を向上させるためにマップの格子数やマップそのものを増やさなければならず、マップを作成するために膨大な予備実験を要する。

また、従来は、ノックをセンサによって検出したら、それをフィードバックして点火時期を遅らせている。しかし、この方法では実際に発生したノックを検出してからフィードバック制御を行うので、制御遅れを回避できない。そこで特許文献１では、クランク角に対する気筒内圧力及び演算開始クランク角の気筒内ガス温度を入力し、それらに基づいて、気筒内ガスが断熱圧縮されるとの仮定のもとで気筒内ガス温度を計算して、その気筒内ガス温度が１２００Ｋ以上に達した時点がノック発生時期であるとして、ノック発生時期を予測している。
特開平７−３３２１４９号公報（第３−５頁、第１，２，５図等）

しかし、特許文献１には、基本としてノックしない点火時期に設定しておき、他の制御量について「ノックの発生しない最適状態」を適切に設定するものである。基本となる「ノックしない点火時期」は運転条件毎にあらかじめ求めておかなければならず、多大な適合工数を要していた。

本発明は、このような従来の問題点に着目してなされたものであり、「ノックしない点火時期」を演算によって導くことで少ない適合工数で適切な点火時期制御を行うことができるとともに、「ノックしない点火時期」の演算にあたって少ない演算負荷と十分な演算精度とを両立する内燃機関の点火時期制御装置を提供することを目的としている。

本発明は以下のような解決手段によって前記課題を解決する。なお、理解を容易にするために本発明の実施形態に対応する符号を付するが、これに限定されるものではない。

本発明は、エンジンの運転状態に基づいて筒内温度及び筒内圧力を推定する筒内温度圧力演算手段（５３１）と、点火された火が伝わる前に筒内温度及び筒内圧力によって未燃ガスが自着火する時間τの逆数１／τを算出する自着火時間算出手段（５３２３）と、自着火時間τの逆数１／τを時間積分してノック発生指標を算出するノック発生指標算出手段（５３２４）と、前記ノック発生指標に基づいて点火時期を制御する点火時期制御手段（５４，５５）と、を備える内燃機関の点火時期制御装置であって、前記ノック発生指標算出手段は、運転状態に基づいて積分区間を変化させる、ことを特徴とする。

本発明によれば、ノック発生指標∫１／τ・ｄｔを導入することで少ない適合工数で「ノックしない点火時期」を判定できるようにした。そしてそのノック発生指標の算出にあたっては、運転状態に応じて積分区間を適切に設定するようにすることで、少ない演算負荷と十分な演算精度との両立を図ることができるようになったのである。

以下では図面等を参照して本発明の実施の形態についてさらに詳しく説明する。

図１は本発明のシステムを説明するための概略図である。

空気は吸気コレクタ２に蓄えられた後、吸気マニホールド３を介して各気筒の燃焼室５に導入される。燃料は各気筒の吸気ポート４に配置された燃料インジェクタ２１より噴射供給される。空気中に噴射された燃料は気化しつつ空気と混合してガス（混合気）を作り、燃焼室５に流入する。この混合気は吸気弁１５が閉じることで燃焼室５内に閉じこめられ、ピストン６の上昇によって圧縮される。

この圧縮混合気に対して高圧火花により点火を行うため、パワートランジスタ内蔵の点火コイルを各気筒に配した電子配電システムの点火装置１１を備える。すなわち、点火装置１１は、バッテリからの電気エネルギーを蓄える点火コイル１３と、点火コイル１３の一次側への通電、遮断を行うパワートランジスタと、燃焼室５の天井に設けられ点火コイル１３の一次電流の遮断によって点火コイル１３の二次側に発生する高電圧を受けて、火花放電を行う点火プラグ１４とからなっている。

圧縮上死点より少し手前で点火プラグ１４により火花が飛ばされ圧縮混合気に着火されると、火炎が広がりやがて爆発的に燃焼し、この燃焼によるガス圧がピストン６を押し下げる仕事を行う。この仕事はクランクシャフト７の回転力として取り出される。燃焼後のガス（排気）は排気弁１６が開いたとき排気通路８へと排出される。

排気通路８には三元触媒９を備える。三元触媒９は排気の空燃比が理論空燃比を中心とした狭い範囲（ウインドウ）にあるとき、排気に含まれるＨＣ、ＣＯ、ＮＯｘといった有害三成分を同時に効率よく除去できる。空燃比は吸入空気量と燃料量の比であるので、エンジンの１サイクル（４サイクルエンジンではクランク角で７２０°区間）当たりに燃焼室５に導入される吸入空気量と、燃料インジェクタ２１からの燃料噴射量との比が理論空燃比となるように、エンジンコントローラ５０ではエアフローメータ３２からの吸入空気流量の信号とクランク角センサ（３３、３４）からの信号に基づいて燃料インジェクタ２１からの燃料噴射量を定めると共に、三元触媒９の上流に設けたＯ₂センサ３５からの信号に基づいて空燃比をフィードバック制御している。

吸気コレクタ２の上流には絞り弁２３がスロットルモータ２４により駆動される、いわゆる電子制御スロットル２２を備える。運転者が要求するトルクはアクセルペダル４１の踏み込み量（アクセル開度）に現れるので、エンジンコントローラ５０ではアクセルセンサ４２からの信号に基づいて目標トルクを定め、この目標トルクを実現するための目標空気量を定め、この目標空気量が得られるようにスロットルモータ２４を介して絞り弁２３の開度を制御する。

吸気弁用カムシャフト２５、排気弁用カムシャフト２６及びクランクシャフト７の各前部にはそれぞれカムスプロケット、クランクスプロケットが取り付けられ、これらスプロケットにタイミングチェーン（不図示）を掛け回すことで、カムシャフト２５、２６がエンジンのクランクシャフト７により駆動されるのであるが、このカムスプロケットと吸気弁用カムシャフト２５との間に介在して、作動角一定のまま吸気弁用カムの位相を連続的に制御し得る吸気バルブタイミングコントロール機構（以下、「吸気ＶＴＣ機構」という。）２７と、カムスプロケットと排気弁用カムシャフト２６との間に介在して、作動角一定のまま排気弁用カムの位相を連続的に制御し得る排気バルブタイミングコントロール機構（以下、「排気ＶＴＣ機構」という。）２８とを備える。吸気弁１５の開閉時期や排気弁１６の開閉時期を変えると燃焼室５に残留する不活性ガスの量が変化する。燃焼室５内の不活性ガスの量が増えるほどポンピングロスが減って燃費がよくなるので、運転条件によりどのくらいの不活性ガスが燃焼室５内に残留したらよいかを目標吸気弁閉時期や目標排気弁閉時期にして予め定めており、エンジンコントローラ５０ではそのときの運転条件（エンジンの負荷と回転速度）より目標吸気弁閉時期と目標排気弁閉時期を定め、それら目標値が得られるように吸気ＶＴＣ機構２７、排気ＶＴＣ機構２８の各アクチュエータを介して吸気弁閉時期と排気弁閉時期を制御する。

吸気温度センサ４３からの吸気温度の信号、吸気圧力センサ４４からの吸気圧力の信号、排気温度センサ４５からの排気温度の信号、排気圧力センサ４６からの排気圧力の信号が、水温センサ３７からの冷却水温の信号と共に入力されるエンジンコントローラ５０では、パワートランジスタ１３を介して点火プラグ１４の一次側電流の遮断時期である点火時期を制御する。

続いて、本発明の特徴である制御の基本的な考え方について説明する。

図２、図３はノック発生メカニズムについて説明する図である。

混合ガスに点火され、燃焼が開始すると、熱発生率が急激に上昇する（図２（Ａ））。すると、筒内温度（図２（Ｂ））、筒内圧力（図２（Ｃ））も上昇する。そして、点火された火が伝わる前に未燃ガスが自着火して急激な燃焼を起こすと、ノックを発生する（図２（Ｃ））。自着火時間τの逆数（１／τ）が着火遅れを表す。この１／τの時間積分値が１となるクランク角でノックが発生する（図２（Ｄ））。すなわち、∫１／τ・ｄｔがノック発生指標である。なお本実施形態では積分開始位置を点火位置にしてある。積分計算を厳密に行うには積分開始位置は吸気弁閉位置とすべきであるが、点火が開始されるまでは１／τの値が微少であるため演算精度上無視できるので演算負担を軽減するために点火が開始した後から積分するように、本実施形態では積分開始位置を点火位置としたのである。また図２では積分終了位置を∫１／τ・ｄｔ＝１を超えたところにしている。図２（Ｄ）の黒丸は一定の演算クランク角間隔毎の演算ポイントを示している。上述の通り∫１／τ・ｄｔ＝１が成立するとノックが発生するので、それ以降について演算を継続するのは意味がない。そこで黒丸の演算ポイントにおいて∫１／τ・ｄｔ＝１を超えたら積分を終了することで余分な演算を省略するようにしたのである。なお積分終了位置は燃焼割合が燃焼終了判定値を超えたときにしてもよい。すなわち燃焼割合が１００％になったら燃焼が終了しているのでノックが発生しないのであるから、燃焼割合に基づいて積分終了位置を決めてもよいのである。ただし質量燃焼割合が１００％となる位置を見極めるのは非常に困難である。そこで燃焼終了判定値として例えば９５％として設定しておくとよい。なお本実施形態では、後述のように図５の５３−２において燃焼割合X_burn(i)が燃焼終了判定値X_burn_endを超えたか否かを判定するとともに、図５の５３−３においてノック発生指標idx_kocr(i)が１を超えたか否かを判定し、いずれかが成立したら割込演算要求を終了することで積分演算を終了している。

図３を参照してさらに詳しく説明する。なお図３（Ａ）は横軸に筒内温度、縦軸に筒内圧力をとったときの１／τマップである。筒内温度及び筒内圧力がわかれば、このマップから１／τの値が求まる。１／τを筒内圧力、筒内温度の履歴に沿って求めて時間積分して∫１／τ・ｄｔを求める。∫１／τ・ｄｔ＝１となるクランク角でノックが発生すると推定できる。

その推定の実証結果を図３（Ｂ）に示す。図３（Ｂ）の横軸にノック発生クランク角の実測値をとり、一方、縦軸に上記のようにして算出したノック発生クランク角の推定値をとる。すると図３（Ｂ）からも明らかなように推定値が実測値と略一致し、上記推定の正しいことが立証できた。

本件発明者は、ノック発生指標として∫１／τ・ｄｔを導入した。これにより、ノック発生指標を筒内圧力及び筒内温度に帰着させることができたのである。筒内圧力及び筒内温度は筒内ガスの最も基本的な状態量であるので、機関機種によらずに略同一の特性とすることが可能になり、またエンジンの暖機状態や環境条件にかかわらず、精度よくリアルタイムにノック発生時期を推定することが可能になったのである。

図４は本発明による点火時期制御装置の制御の全体を示す制御ブロック図である。

エンジンコントローラ５０は、運転状態５１に基づいて燃焼速度５２１、燃焼期間５２２、基本点火時期５２３を順次算出し、また、ノック発生指標５３を演算し、進角補正／限界（リタード幅）５４を求め、それらよりノック限界点火時期５５を算出する。

続いて、各構成ブロックについてさらに詳しく説明する。なお、基本点火時期算出部５２の具体的な内容については特開２００３−１４８２３６号公報において詳細に説明されているので本件では説明を省略する。

図５はノック発生指標ブロック５３を詳述する図である。

ノック発生指標ブロック５３においては、基本点火時期MBTCAL及び演算クランク角間隔delta_thetaを入力し、ノック発生指標演算開始要求JOB_REQに基づき、演算カウンタ(i)ごとに演算クランク角CA_calc(i)を順次計算する。そしてその演算クランク角CA_calc(i)に基づいてノック発生指標５３０を演算し、演算カウンタ(i)ごとにノック発生指標idx_kocr(i)及び質量燃焼割合X_burn(i)を順次計算する。なお質量燃焼割合X_burn(i)が燃焼終了時燃焼割合X_burn_endを超えるまでは、またはノック発生指標idx_kocr(i)が１を超えるまでは割込演算を要求する。

ノック発生指標演算ブロック５３０については図６を参照しながら説明する。ノック発生指標演算ブロック５３０は、筒内平均温度・圧力演算ブロック５３１と、ノック発生指標演算ブロック５３２とを有する。筒内平均温度・圧力演算ブロック５３１においては、クランク角度CA、エンジン回転数NE、基本点火時期MBTCAL、噴射パルス幅TP、シリンダ容積初期値V_cyl_ini、筒内温度初期値T_cyl_ini、筒内圧初期値P_cyl_ini、点火むだ時間IGNDEAD、燃焼期間BURN、基準燃焼割合X_ref、残ガス率MRESFRに基づいて、筒内温度T_cyl、筒内圧P_cyl、質量燃焼割合X_burnを算出する。ノック発生指標演算ブロック５３２においては、クランク角CA、エンジン回転数NE、基本点火時期MBTCAL、筒内温度T_cyl、筒内圧P_cylに基づいて、ノック発生指標idx_kocrを算出する。筒内平均温度・圧力演算ブロック５３１については図７を参照しながら詳述し、ノック発生指標演算ブロック５３２については図２１を参照しながら詳述する。

まず、図７を参照して、筒内平均温度・圧力演算ブロック５３１について説明する。

筒内平均温度・圧力演算ブロック５３１は、シリンダ容積演算ブロック５３１１と、発熱量演算ブロック５３１２と、冷却損失演算ブロック５３１３と、筒内温度・圧力演算ブロック５３１４とを有する。

シリンダ容積演算ブロック５３１１においては、クランク角度CAに基づいて、シリンダ容積V_cyl、ピストン変位x_pisを算出する。発熱量演算ブロック５３１２においては、クランク角度CA、噴射パルス幅TP、基本点火時期MBTCAL、点火むだ時間IGNDEAD、燃焼期間BURN、基準燃焼割合X_ref、残ガス率MRESFRに基づいて、発熱率Q_burn、物性演算用質量燃焼割合X_burn_r、質量燃焼割合X_burnを算出する。冷却損失演算ブロック５３１３においては、シリンダ容積V_cyl、ピストン変位x_pis、残ガス率MRESFR、筒内温度初期値T_cyl_ini、シリンダ容積初期値V_cyl_ini、筒内圧初期値P_cyl_ini、エンジン回転数NE、筒内圧前回値P_cyl_z、筒内温度前回値T_cyl_zに基づいて、冷却損失Q_lossを算出する。筒内温度・圧力演算ブロック５３１４においては、噴射パルス幅TP、シリンダ容積V_cyl、発熱率Q_burn、物性演算用質量燃焼割合X_burn_r、冷却損失Q_loss、筒内圧前回値P_cyl_z、筒内温度前回値T_cyl_zに基づいて、筒内圧P_cyl、筒内温度T_cylを算出する。

続いて、図７に示した筒内平均温度・圧力演算ブロック５３１を構成する各ブロックについてさらに詳しく説明する。

はじめに、シリンダ容積演算ブロック５３１１の詳細について、図８を参照して説明する。シリンダ容積演算ブロック５３１１は、ピストン変位演算ブロック５３１１１と、シリンダ容積演算ブロック５３１１２とを有する。ピストン変位演算ブロック５３１１１において、クランク角度CAからピストン変位x_pisを算出する。シリンダ容積演算ブロック５３１１２において、シリンダ容積V_cylを算出する。

次に、発熱量演算ブロック５３１２の詳細について、図９を参照して説明する。

発熱量演算ブロック５３１２は、Wiebe関数定数設定ブロック５３１２１と、熱発生率算出ブロック５３１２２と、質量燃焼割合算出ブロック５３１２３と、発熱量／１シリンダ／１サイクル算出ブロック５３１２４と、発熱量／１シリンダ／演算間隔算出ブロック５３１２５とを有する。

ここでWiebe関数について説明する。熱発生率の特性を表すものとして、質量燃焼割合X(%)をクランク角θ(deg)に対する関数として一般に
Ｘ＝１−exp[-ａ・[(θ-θ_S)/θ_B]ⁿ⁺¹] …(1)
θ_S ：実燃焼開始クランク角(deg)
θ_B ：実燃焼期間(deg)
ａ，ｎ：定数
と記述することができる。上記ａ，ｎは燃焼室形状や点火プラグ位置、および燃焼室内のガス流動特性等、エンジン機種によって定まる。これは実験等によって予め決めておくことができる。この(1)式がWiebe関数である。

上記(1)式をθで微分すると熱発生率を表す関数となり、
ｄＸ／ｄθ
＝ａ(ｎ＋１)/θ_B・[(θ-θ_S)/θ_B]ⁿ・exp[-ａ・[(θ-θ_S)/θ_B]ⁿ⁺¹] …(2)
(1)式を実燃焼期間θ_Bについて解けば、
θ_B＝θ/[-log(１−Ｘ)/ａ]^1/(n+1) …(3)
となる。

以上を踏まえて、Wiebe関数定数設定ブロック５３１２１について、図１０を参照しながら説明する。なお図１０は上述の(3)式を制御ブロック図として記載したものである。すなわち、(3)式において、θ＝BURN、X＝X_refを代入すれば実燃焼期間θ_B(BURN_r)を算出することができる。Wiebe関数定数設定ブロック５３１２１は、基本点火時期MBTCAL、点火むだ時間IGNDEADに基づいて、熱発生開始時期BURN_iniを算出するとともに、燃焼期間BURN、基準燃焼割合X_refに基づいて、実燃焼期間BURN_rを算出する。

また熱発生率算出ブロック５３１２２については、図１１を参照しながら説明する。なお図１１は上述の(2)式を制御ブロック図として記載したものである。すなわち、(2)式において、θ＝CA、θ_S＝BURN_ini、θ_B＝BURN_rを代入すれば熱発生率ｄＸ／ｄθを算出することができる。熱発生率算出ブロック５３１２２は、クランク角度CA、熱発生開始時期BURN_ini、実燃焼期間BURN_rに基づいて、熱発生率を算出する。

質量燃焼割合算出ブロック５３１２３については、図１２を参照しながら説明する。なお図１２は上述の(1)式を制御ブロック図として記載したものである。すなわち、(1)式において、θ＝CA、θ_S＝BURN_ini、θ_B＝BURN_rを代入すれば質量燃焼割合Ｘ(X_burn)を算出することができる。質量燃焼割合算出ブロック５３１２３は、クランク角度CA、熱発生開始時期BURN_ini、実燃焼期間BURN_r、残ガス率MRESFRに基づいて、質量燃焼割合X_burn、物性演算用質量燃焼割合X_burn_rを算出する。

発熱量／１シリンダ／１サイクル算出ブロック５３１２４については、図１３を参照しながら説明する。発熱量／１シリンダ／１サイクル算出ブロック５３１２４は、噴射パルス幅TP、定数KCONST、目標空燃比TABYF、低発熱量Heat_Lowerに基づいて、１サイクル当たりの１シリンダ毎の発熱量を算出する。

発熱量／１シリンダ／演算間隔算出ブロック５３１２５については、図１４を参照しながら説明する。発熱量／１シリンダ／演算間隔算出ブロック５３１２５は、熱発生率及び１サイクル当たりの１シリンダ毎の発熱量に基づいて、発熱量Q_burnを算出する。

以上が発熱量演算ブロック５３１２の説明である。続いて、図７の冷却損失演算ブロック５３１３の詳細については、図１５を参照して説明する。

冷却損失演算ブロック５３１３は、燃焼室表面積算出ブロック５３１３１と、平均ピストン速度算出ブロック５３１３２と、モータリング−ガス物性値算出ブロック５３１３３と、モータリング−筒内圧算出ブロック５３１３４と、モータリング−筒内温度算出ブロック５３１３５と、熱伝達係数算出ブロック５３１３６と、冷却損失算出ブロック５３１３７とを有する。

燃焼室表面積算出ブロック５３１３１は、ピストン変位x_pisに基づいて燃焼室表面積を算出する。平均ピストン速度算出ブロック５３１３２は、エンジン回転数に基づいて平均ピストン速度を算出する。

モータリング−ガス物性値算出ブロック５３１３３は、残ガス率MRESFRから求めた質量燃焼割合と、筒内温度とに基づいて、比熱比κを算出する。具体的な算出方法は、ガス物性値算出ブロック５３１４１（図１９）と同じであるので、ここでは説明を省略する。

モータリング−筒内圧算出ブロック５３１３４は、筒内圧前回値、シリンダ容積前回値、シリンダ容積今回値、比熱比κに基づいて、筒内圧P_cylを算出する。なお本演算では断熱変化であるとして計算している。モータリング−筒内温度算出ブロック５３１３５は、筒内温度前回値、シリンダ容積前回値、シリンダ容積今回値、比熱比κに基づいて、筒内温度T_cylを算出する。なお本演算では断熱変化であるとして計算している。

次に、熱伝達係数算出ブロック５３１３６について説明する前に、冷却損失Woschniの式について説明する。

冷却損失Q_Lは以下により算出される、
Q_L = ∫_t1 ^t2h・A・( T - T_W )・dt …(4)
h = 110・d^-0.2・P^0.8・T^-0.53・[C₁・c_m + C₂・(V_S・T_I)/(P_I・V_I)・(P - P_M)]^0.8 …(5)
h ：熱伝達係数 [ kcal/(m²・h・K) ]
d ：シリンダボア径 [ m ]
V_S：シリンダ容積 [ m³ ]
c_m ：平均ピストン速度 [ m/s ]
P_I ：吸気弁閉時筒内圧 [ kgf/cm² ]
V_I ：吸気弁閉時筒内容積 [ m³ ]
T_I ：吸気弁閉時筒内温度 [ K ]
P_M ：モータリング筒内圧 [ kgf/cm² ]
P ：筒内圧 [ kgf/cm² ]
T ：筒内温度 [ K ]
C₁ ：6.18（排気），2.28（圧縮＆膨張） [−]
C₂ ：3.24×10^-3 [ m/(s・K) ]
なお、上記でC₁、C₂は一般に用いられている定数である。

図１６は(5)式を制御ブロックとして図示したものである。この図１６に示すように、熱伝達係数算出ブロック５３１３６は、筒内圧、筒内温度、平均ピストン速度、モータリング筒内圧、筒内温度初期値、シリンダ容積初期値、筒内圧初期値に基づいて、熱伝達係数を算出する。

図１７は(4)式を制御ブロックとして図示したものである。この図１７に示すように、冷却損失算出ブロック５３１３７は、筒内温度、燃焼室壁温度、燃焼室表面積、熱伝達係数、演算間隔、エンジン回転数に基づいて、冷却損失Q_lossを算出する。

以上が冷却損失演算ブロック５３１３の説明である。次に図７の筒内温度・圧力演算ブロック５３１４について、図１８を参照して説明する。

筒内温度・圧力演算ブロック５３１４は、ガス物性値算出ブロック５３１４１と、全ガス量算出ブロック５３１４２と、断熱変化筒内圧算出ブロック５３１４３と、断熱変化筒内温度算出ブロック５３１４４と、筒内温度算出ブロック５３１４５と、筒内圧算出ブロック５３１４６とを有する。

ガス物性値算出ブロック５３１４１は、質量燃焼割合、筒内温度に基づいて、ガス物性値（比熱比、定容比熱）を算出する。ここで、既燃ガスの割合（「物性演算用質量燃焼割合」を適用）に基づいてガス物性を算出する場合、
(ガス物性値)=(既燃ガス物性値)*(既燃ガス割合)+(未燃ガス物性値)*(1-既燃ガス割合) …(6)
であるので、それを図１９のブロック図に示した。図１９に示すように、モータリング−ガス物性値算出ブロック５３１３３は、筒内温度T_cyl(初期値はT_cyl_ini)に基づいて、既燃分の定容比熱及び比熱比を求めるとともに、同じく筒内温度T_cyl(初期値はT_cyl_ini)に基づいて、未燃分の定容比熱及び比熱比を求め、既燃分及び未燃分の割合を掛け合わせて筒内ガスの物性値を算出する。

全ガス量算出ブロック５３１４２は、図２０に示すように、噴射パルス幅TP、KCONST、空燃比ABYFに基づいて全ガス量を算出する。

断熱変化筒内圧算出ブロック５３１４３は、筒内圧前回値、シリンダ容積前回値、シリンダ容積今回値、比熱比κに基づいて、断熱変化筒内圧を算出する。断熱変化筒内温度算出ブロック５３１４４は、筒内温度前回値、シリンダ容積前回値、シリンダ容積今回値、比熱比κに基づいて、断熱変化筒内温度を算出する。筒内温度算出ブロック５３１４５は、断熱変化筒内温度、発熱量、冷却損失、全ガス量、定容比熱に基づいて筒内温度T_cylを算出する。筒内圧算出ブロック５３１４６は、断熱変化筒内圧、筒内温度、断熱変化筒内温度に基づいて、筒内圧P_cylを算出する。

以上が筒内温度・圧力演算ブロック５３１４の説明であり、これによって図６の筒内平均温度・圧力演算ブロック５３１についてすべて説明した。次に、図６のノック発生指標演算ブロック５３２について、図２１を参照しながら詳述する。

ノック発生指標演算ブロック５３２は、resetフラグ演算ブロック５３２１と、未燃ガス温度演算ブロック５３２２と、瞬時着火遅れ逆数演算ブロック５３２３と、ノック発生指標演算ブロック５３２４とを有する。

resetフラグ演算ブロック５３２１の詳細を図２２に示す。resetフラグ演算ブロック５３２１は、クランク角度CAが基本点火時期MBTCALを超えたらresetフラグを立てる。未燃ガス温度演算ブロック５３２２は、図２３に示すように未燃ガス物性値算出ブロック５３２２１と、未燃ガス温度算出ブロック５３２２２とを有する。未燃ガス物性値算出ブロック５３２２１は、ガス物性値算出ブロック５３１４１（図１９）と同じため、説明を省略する。質量燃焼割合としてゼロを入力して未燃ガスだけの物性値を算出している。また未燃ガス温度算出ブロック５３２２２は、未燃ガス温度前回値、筒内圧今回値、筒内圧前回値、比熱比κに基づいて未燃ガス温度T_ubを算出する。

瞬時着火遅れ逆数演算ブロック５３２３においては、未燃ガス温度T_ub、筒内圧P_cylに基づいて瞬時着火遅れ逆数tauinvを求める。なお、このマップは予め実験によって求めらておく。ノック発生指標演算ブロック５３２４は、図２４に示すようにエンジン回転数NEから時間を算出し、その時間によって瞬時着火遅れ逆数tauinvを積分してノック発生指標idx_kocrを算出する。なおresetフラグが立っていたら、初期値ゼロをノック発生指標idx_kocrとする。

以上が図６のノック発生指標演算ブロック５３２の説明であり、これによって図４のノック発生指標５３についてすべて説明した。続いて、図４のリタード幅算出ブロック５４について図２５を参照しながら詳述する。なお図２５はリタード幅算出ブロック５４の詳細な制御ブロックを示す図である。

リタード幅算出ブロック５４は、クランク角算出ブロック５４１と、質量燃焼割合算出ブロック５４２と、ノック強度算出ブロック５４３と、リタード幅算出ブロック５４４とを備える。

クランク角算出ブロック５４１は、予め実験によって求められている線図に基づいて、ノック発生指標idx_kocr(i)、演算クランク角CA_calc(i)からノック発生クランク角CA_knkを求める。

質量燃焼割合算出ブロック５４２は、予め実験によって求められている線図に基づいて、演算クランク角CA_calc(i)、質量燃焼割合X_burn(i)、ノック発生クランク角CA_knkからノック発生質量燃焼割合X_bknkを求める。

ノック強度算出ブロック５４３は、予め実験によって求められている線図に基づいて、ノック発生クランク角CA_knk、ノック発生時質量燃焼割合X_bknkからノック強度I_knkを求める。なお図中の線図において、矢印方向（下方）が遅角側である。

リタード幅算出ブロック５４４は、予め実験によって求められている線図に基づいて、ノック強度I_knkからリタード幅（進角補正限界）CA_rtdを求める。

次に、図４のノック限界点火時期算出ブロック５５について図２６を参照しながら詳述する。なお図２６はノック限界点火時期算出ブロック５５の詳細な制御ブロックを示す図である。基本点火時期MBTCALからリタード幅CA_rtdを遅角させて点火時期ADVを求める。

以上で、図４に示すエンジンコントローラ５０の制御についてすべて説明した。

以上、詳細に説明したように、請求項１に記載した発明によれば、積分区間を運転状態に応じて適切に設定するようにしたので、演算精度と演算負荷との両立を図ることができるのである。

請求項２に記載した発明によれば、点火位置から積分計算を開始するようにした（５３−１）。積分計算を厳密に行うには積分開始位置は吸気弁閉位置とすべきであるが、実際には着火遅れの積分値は燃焼が始まるまでは、それ以降と比較してほとんど増大しない。このため点火位置を積分開始位置としても、演算精度は必要十分に確保でき、かつ無用な区間を積分しないので、演算負荷を軽減できるのである。

請求項３に記載した発明によれば、ノックが発生しない正常燃焼時の燃焼速度に基づいて積分計算を終了するようにした。ノックは正常に燃焼する場合の燃焼終了時期よりも早い時点で未燃ガスが爆発的に燃焼することによって発生する。したがって、正常燃焼と仮定した場合の燃焼期間分だけ積分すれば十分である。燃焼期間は燃焼速度の関数であるので、燃焼速度に基づいて積分区間を決めることができるのである。

請求項４に記載した発明によれば、正常燃焼時の燃焼速度に基づいて算出した燃焼割合が燃焼終了判定値以上となったら積分計算を終了するようにした（５３−２）。燃焼期間（または燃焼終了位置）は燃焼割合に基づいて判断することで、より適切に決めることができる。実際には燃焼割合が１００％となる位置を見極めるのは困難であるので、燃焼割合が燃焼終了判定値以上となる時点を燃焼終了としても精度上問題ない。このようにすることで積分区間を適切に設定できるのである。

請求項５に記載した発明によれば、ノック発生指標がノック指標判定値以上となったら積分計算を終了するようにした（５３−３）。ノック発生指標がノック指標判定値以上となった時点でノックが発生するので、それ以降積分計算を継続するのは意味がない。そこでノック発生指標がノック指標判定値以上となった時点で積分計算を終了することで余分な演算を省略できるのである。

請求項６に記載した発明によれば、燃焼割合が燃焼終了判定値以上となった場合（５３−２）、又はノック発生指標がノック指標判定値以上となった場合（５３−３）、のいずれかの条件が成立したら積分計算を終了するようにした。このようにすることで、ノック発生の有無にかかわらず、運転状態に応じて積分区間を適切かつ必要十分に設定することができるのである。

以上説明した実施形態に限定されることなく、その技術的思想の範囲内において種々の変形や変更が可能であり、それらも本発明と均等であることは明白である。

本発明のシステムを説明するための概略図である。ノック発生メカニズムについて説明する図である。ノック発生メカニズムについて説明する図である。本発明による点火時期制御装置の制御の全体を示す制御ブロック図である。ノック発生指標ブロック５３を詳述する図である。ノック発生指標演算ブロック５３０を示す図である。筒内平均温度・圧力演算ブロック５３１を示す図である。シリンダ容積演算ブロック５３１１を示す図である。発熱量演算ブロック５３１２を示す図である。 Wiebe関数定数設定ブロック５３１２１を示す図である。熱発生率算出ブロック５３１２２を示す図である。質量燃焼割合算出ブロック５３１２３を示す図である。発熱量／１シリンダ／１サイクル算出ブロック５３１２４を示す図である。発熱量／１シリンダ／演算間隔算出ブロック５３１２５を示す図である。冷却損失演算ブロック５３１３の詳細を示す図である。熱伝達係数算出ブロック５３１３６を示す図である。冷却損失算出ブロック５３１３７を示す図である。筒内温度・圧力演算ブロック５３１４を示す図である。モータリング−ガス物性値算出ブロック５３１３３を示す図である。全ガス量算出ブロック５３１４２を示す図である。ノック発生指標演算ブロック５３２を示す図である。 resetフラグ演算ブロック５３２１を示す図である。未燃ガス温度演算ブロック５３２２を示す図である。ノック発生指標演算ブロック５３２４を示す図である。リタード幅算出ブロック５４の詳細な制御ブロックを示す図である。ノック限界点火時期算出ブロック５５を示す図である。

符号の説明

１エンジン
５０エンジンコントローラ
５１運転状態検出手段
５２３基本点火時期算出手段
５３ノック発生指標ブロック
５３０ノック発生指標演算ブロック
５３１筒内平均温度・圧力演算ブロック
５３１２１ Wiebe関数定数設定ブロック（燃焼速度算出手段）
５３１２３質量燃焼割合算出ブロック（燃焼割合算出手段）
５３２ノック発生指標演算ブロック（ノック発生指標算出手段）
５３２３瞬時着火遅れ逆数演算ブロック（自着火時間算出手段）
５３２４ノック発生指標演算ブロック（ノック発生指標算出手段）
５４リタード幅算出ブロック（点火時期制御手段（点火時期補正量算出手段））
５５ノック限界点火時期算出ブロック（点火時期制御手段）

Claims

エンジンの運転状態に基づいて筒内温度及び筒内圧力を推定する筒内温度圧力演算手段と、
点火された火が伝わる前に筒内温度及び筒内圧力によって未燃ガスが自着火する時間τの逆数１／τを算出する自着火時間算出手段と、
自着火時間τの逆数１／τを時間積分してノック発生指標を算出するノック発生指標算出手段と、
前記ノック発生指標に基づいて点火時期を制御する点火時期制御手段と、
を備える内燃機関の点火時期制御装置であって、
前記ノック発生指標算出手段は、運転状態に基づいて積分区間を変化させる、
ことを特徴とする内燃機関の点火時期制御装置。
運転状態に基づいて基本点火時期を算出する基本点火時期算出手段を備え、
前記ノック発生指標算出手段は、前記基本点火時期から積分計算を開始する、
ことを特徴とする請求項１に記載の内燃機関の点火時期制御装置。
ノックが発生しない正常燃焼時の燃焼速度を算出する燃焼速度算出手段を備え、
前記ノック発生指標算出手段は、前記正常燃焼時の燃焼速度に基づいて積分計算を終了する、
ことを特徴とする請求項１又は請求項２に記載の内燃機関の点火時期制御装置。
前記正常燃焼時の燃焼速度に基づいて燃焼割合を算出する燃焼割合算出手段を備え、
前記ノック発生指標算出手段は、前記燃焼割合が燃焼終了判定値以上となったら積分計算を終了する、
ことを特徴とする請求項３に記載の内燃機関の点火時期制御装置。
前記ノック発生指標算出手段は、前記ノック発生指標がノック指標判定値以上となったら積分計算を終了する、
ことを特徴とする請求項１又は請求項２に記載の内燃機関の点火時期制御装置。
ノックが発生しない正常燃焼時の燃焼速度を算出する燃焼速度算出手段と、
前記正常燃焼時の燃焼速度に基づいて燃焼割合を算出する燃焼割合算出手段と、
を備え、
前記ノック発生指標算出手段は、前記燃焼割合が燃焼終了判定値以上となった場合、又は前記ノック発生指標がノック指標判定値以上となった場合、のいずれかの条件が成立したら積分計算を終了する、
ことを特徴とする請求項１又は請求項２に記載の内燃機関の点火時期制御装置。