(57)【要約】【課題】多くの適合工数を必要とすることなく、適切な点火時期制御を行うことが可能な内燃機関の点火時期制御装置を提供する。【解決手段】内燃機関の運転状態を検出する運転状態検出手段５１と、運転状態検出手段で検出した運転状態に基づいて、筒内ガスの着火遅れを検出する着火遅れ検出手段と、着火遅れ検出手段で検出した着火遅れに基づいてノック発生の指標となるノック発生指標を算出するノック発生指標算出手段５３と、ノック発生指標算出手段で算出したノック発生指標に基づいてノックを発生しない限界の点火時期を算出するノック限界点火時期算出手段５４〜５５と、を有することを特徴とする。【選択図】図１
2005233111.tif 000002
【特許請求の範囲】
【請求項１】
内燃機関の運転状態を検出する運転状態検出手段と、
前記運転状態検出手段で検出した運転状態に基づいて、筒内ガスの着火遅れを検出する着火遅れ検出手段と、
前記着火遅れ検出手段で検出した着火遅れに基づいてノック発生の指標となるノック発生指標を算出するノック発生指標算出手段と、
前記ノック発生指標算出手段で算出したノック発生指標に基づいてノックを発生しない限界の点火時期を算出するノック限界点火時期算出手段と、
を有する内燃機関の点火時期制御装置。
【請求項２】
前記着火遅れ検出手段は、前記運転状態検出手段で検出した運転状態に基づいて筒内温度及び筒内圧力を算出する温度圧力算出手段を有し、その筒内温度及び筒内圧力に基づいて、筒内ガスの着火遅れを検出する、
ことを特徴とする請求項１に記載の内燃機関の点火時期制御装置。
【請求項３】
前記着火遅れ検出手段は、前記運転状態検出手段で検出した運転状態と、前記温度圧力算出手段で算出した筒内温度及び筒内圧力と、に基づいて筒内の未燃ガスの温度を算出する未燃ガス温度算出手段を有し、その未燃ガス温度及び筒内圧力に基づいて筒内ガスの着火遅れを検出する、
ことを特徴とする請求項２に記載の内燃機関の点火時期制御装置。
【請求項４】
前記ノック発生指標算出手段は、前記ノックが発生する時期の直前の筒内ガスの質量燃焼割合を算出するノック直前筒内ガス質量燃焼割合算出手段と、前記着火遅れ検出手段で検出した着火遅れに基づいてノックが発生する時期を算出するノック時期算出手段とを備え、そのノック時期算出手段で算出したノック発生時期に基づいてノック発生の指標となるノック発生指標を算出し、
前記ノック限界点火時期算出手段は、前記ノック発生直前の筒内ガスの質量燃焼割合及び前記ノック発生時期に基づいてノックを発生しない限界の点火時期を算出する、
ことを特徴とする請求項１から請求項３までのいずれか１項に記載の内燃機関の点火時期制御装置。
【請求項５】
前記ノック発生指標に対するしきい値を設定するしきい値設定手段を備え、
前記ノック時期算出手段は、前記しきい値と前記ノック発生指標とが略一致するクランク角をノック発生時期として算出する、
ことを特徴とする請求項４に記載の内燃機関の点火時期制御装置。
【請求項６】
運転状態に基づいて基本点火時期を算出する基本点火時期算出手段と、
少なくともノック発生時期又はノック発生直前の筒内ガスの質量燃焼割合に基づいて、前記基本点火時期を補正する点火時期補正手段と、
備え、
前記点火時期補正手段は、ノック発生時期が進角側であるほど、又はノック発生直前の筒内ガスの質量燃焼割合が小さいほど、点火時期を遅角側に補正する、
ことを特徴とする請求項４又は請求項５に記載の内燃機関の点火時期制御装置。
【請求項７】
前記温度圧力算出手段は、前記運転状態に基づいて筒内ガスの熱発生パターンを算出する熱発生パターン算出手段を備え、その熱発生パターンに基づいて筒内温度及び筒内圧力を算出し、
前記ノック直前筒内ガス質量燃焼割合算出手段は、前記熱発生パターンと、前記ノック発生時期と、に基づいて、ノックが発生する直前の筒内ガスの質量燃焼割合を算出する、
ことを特徴とする請求項４から請求項６までのいずれか１項に記載の内燃機関の点火時期制御装置。
【請求項８】
前記ノック発生指標算出手段は、所定の少なくとも一点以上の演算クランク角を設定する演算クランク角設定手段を有し、
前記温度圧力算出手段は、前記演算クランク角の各点における筒内温度及び筒内圧力を算出し、
前記着火遅れ検出手段は、前記各点の瞬時着火遅れを算出し、
前記ノック発生指標算出手段は、前記瞬時着火遅れの逆数を時間積分してノック発生指標を算出する、
ことを特徴とする請求項２から請求項７までのいずれか１項に記載の内燃機関の点火時期制御装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ノックを発生しないように点火時期を制御する内燃機関の点火時期制御装置に関するものである。
【背景技術】
【０００２】
エンジンの点火時期は燃費や出力に影響し、また点火時期がずれるとノッキングや失火を生じるおそれがあるので、その制御は非常に重要である。通常は、１０〜１５degATDCでシリンダ内圧力が最大圧になるように点火時期を制御している。しかし、点火時期設定の基本となる最大トルク最小進角点火時期(Minimmum Spark Advance for Best Torque；以下「ＭＢＴ」と略す)は、エンジン回転速度、負荷、空燃比、ＥＧＲ率等によって時々刻々と変動する。そこで従来は、エンジン回転速度や負荷に応じて設定した基本点火時期マップと、さまざまな運転状態を想定した補正値マップとを備え、運転状態に応じてそれらのマップを参照して制御していた。しかし、そのような方法では、精度を向上させるためにマップの格子数やマップそのものを増やさなければならず、マップを作成するために膨大な予備実験を要する。
【０００３】
また、従来は、ノックをセンサによって検出したら、それをフィードバックして点火時期を遅らせている。しかし、この方法では実際に発生したノックを検出してからフィードバック制御を行うので、制御遅れを回避できない。そこで特許文献１では、クランク角に対する気筒内圧力及び演算開始クランク角の気筒内ガス温度を入力し、それらに基づいて、気筒内ガスが断熱圧縮されるとの仮定のもとで気筒内ガス温度を計算して、その気筒内ガス温度が１２００Ｋ以上に達した時点がノック発生時期であるとして、ノック発生時期を予測している。
【特許文献１】特開平７−３３２１４９号公報（第３−５頁、第１，２，５図等）
【発明の開示】
【発明が解決しようとする課題】
【０００４】
ところで、１／τマップから自着火時期を求めるものが知られている。これによれば圧縮開始から燃焼に至るまでの燃焼室内の未燃ガスの圧力温度をトレースして知ることで、より精度よく自着火時期を求めることができる。しかしながら、特許文献１では圧力センサを用いて筒内圧力を直接検出している。このようにセンサを用いる方式では、運転条件が種々変化した場合には、その都度適合する必要があるので、適合工数が多大なものとなっていた。
【０００５】
本発明は、このような従来の問題点に着目してなされたものであり、多くの適合工数を必要とすることなく、適切な点火時期制御を行うことが可能な内燃機関の点火時期制御装置を提供することを目的としている。
【課題を解決するための手段】
【０００６】
本発明は以下のような解決手段によって前記課題を解決する。なお、理解を容易にするために本発明の実施形態に対応する符号を付するが、これに限定されるものではない。
【０００７】
本発明は、内燃機関の運転状態を検出する運転状態検出手段（５１）と、前記運転状態検出手段で検出した運転状態に基づいて、筒内ガスの着火遅れを検出する着火遅れ検出手段（５３）と、前記着火遅れ検出手段（５３）で検出した着火遅れに基づいてノック発生の指標となるノック発生指標を算出するノック発生指標算出手段（５３）と、前記ノック発生指標算出手段で算出したノック発生指標に基づいてノックを発生しない限界の点火時期を算出するノック限界点火時期算出手段（５４〜５５）と、を有することを特徴とする。
【発明の効果】
【０００８】
可燃性ガスは、その種類及び状態によって定まる所定の着火遅れ時間経過後に爆発的に燃焼することが知られている。また着火遅れ時間とノック発生のしやすさには、強い相関がある。そこで本発明では着火遅れに基づいてノック限界点火時期を算出する構成とすることで、適切にノック限界点火時期を算出できるようになったのである。
【発明を実施するための最良の形態】
【０００９】
以下では図面等を参照して本発明の実施の形態についてさらに詳しく説明する。
【００１０】
図１は本発明のシステムを説明するための概略図である。
【００１１】
空気は吸気コレクタ２に蓄えられた後、吸気マニホールド３を介して各気筒の燃焼室５に導入される。燃料は各気筒の吸気ポート４に配置された燃料インジェクタ２１より噴射供給される。空気中に噴射された燃料は気化しつつ空気と混合してガス（混合気）を作り、燃焼室５に流入する。この混合気は吸気弁１５が閉じることで燃焼室５内に閉じこめられ、ピストン６の上昇によって圧縮される。
【００１２】
この圧縮混合気に対して高圧火花により点火を行うため、パワートランジスタ内蔵の点火コイルを各気筒に配した電子配電システムの点火装置１１を備える。すなわち、点火装置１１は、バッテリからの電気エネルギーを蓄える点火コイル１３と、点火コイル１３の一次側への通電、遮断を行うパワートランジスタと、燃焼室５の天井に設けられ点火コイル１３の一次電流の遮断によって点火コイル１３の二次側に発生する高電圧を受けて、火花放電を行う点火プラグ１４とからなっている。
【００１３】
圧縮上死点より少し手前で点火プラグ１４により火花が飛ばされ圧縮混合気に着火されると、火炎が広がりやがて爆発的に燃焼し、この燃焼によるガス圧がピストン６を押し下げる仕事を行う。この仕事はクランクシャフト７の回転力として取り出される。燃焼後のガス（排気）は排気弁１６が開いたとき排気通路８へと排出される。
【００１４】
排気通路８には三元触媒９を備える。三元触媒９は排気の空燃比が理論空燃比を中心とした狭い範囲（ウインドウ）にあるとき、排気に含まれるＨＣ、ＣＯ、ＮＯｘといった有害三成分を同時に効率よく除去できる。空燃比は吸入空気量と燃料量の比であるので、エンジンの１サイクル（４サイクルエンジンではクランク角で７２０°区間）当たりに燃焼室５に導入される吸入空気量と、燃料インジェクタ２１からの燃料噴射量との比が理論空燃比となるように、エンジンコントローラ５０ではエアフローメータ３２からの吸入空気流量の信号とクランク角センサ（３３、３４）からの信号に基づいて燃料インジェクタ２１からの燃料噴射量を定めると共に、三元触媒９の上流に設けたＯ2センサ３５からの信号に基づいて空燃比をフィードバック制御している。
【００１５】
吸気コレクタ２の上流には絞り弁２３がスロットルモータ２４により駆動される、いわゆる電子制御スロットル２２を備える。運転者が要求するトルクはアクセルペダル４１の踏み込み量（アクセル開度）に現れるので、エンジンコントローラ５０ではアクセルセンサ４２からの信号に基づいて目標トルクを定め、この目標トルクを実現するための目標空気量を定め、この目標空気量が得られるようにスロットルモータ２４を介して絞り弁２３の開度を制御する。
【００１６】
吸気弁用カムシャフト２５、排気弁用カムシャフト２６及びクランクシャフト７の各前部にはそれぞれカムスプロケット、クランクスプロケットが取り付けられ、これらスプロケットにタイミングチェーン（不図示）を掛け回すことで、カムシャフト２５、２６がエンジンのクランクシャフト７により駆動されるのであるが、このカムスプロケットと吸気弁用カムシャフト２５との間に介在して、作動角一定のまま吸気弁用カムの位相を連続的に制御し得る吸気バルブタイミングコントロール機構（以下、「吸気ＶＴＣ機構」という。）２７と、カムスプロケットと排気弁用カムシャフト２６との間に介在して、作動角一定のまま排気弁用カムの位相を連続的に制御し得る排気バルブタイミングコントロール機構（以下、「排気ＶＴＣ機構」という。）２８とを備える。吸気弁１５の開閉時期や排気弁１６の開閉時期を変えると燃焼室５に残留する不活性ガスの量が変化する。燃焼室５内の不活性ガスの量が増えるほどポンピングロスが減って燃費がよくなるので、運転条件によりどのくらいの不活性ガスが燃焼室５内に残留したらよいかを目標吸気弁閉時期や目標排気弁閉時期にして予め定めており、エンジンコントローラ５０ではそのときの運転条件（エンジンの負荷と回転速度）より目標吸気弁閉時期と目標排気弁閉時期を定め、それら目標値が得られるように吸気ＶＴＣ機構２７、排気ＶＴＣ機構２８の各アクチュエータを介して吸気弁閉時期と排気弁閉時期を制御する。
【００１７】
吸気温度センサ４３からの吸気温度の信号、吸気圧力センサ４４からの吸気圧力の信号、排気温度センサ４５からの排気温度の信号、排気圧力センサ４６からの排気圧力の信号が、水温センサ３７からの冷却水温の信号と共に入力されるエンジンコントローラ５０では、パワートランジスタ１３を介して点火プラグ１４の一次側電流の遮断時期である点火時期を制御する。
【００１８】
続いて、本発明の特徴である制御の基本的な考え方について説明する。
【００１９】
図２、図３はノック発生メカニズムについて説明する図である。
【００２０】
混合ガスに点火され、燃焼が開始すると、熱発生率が急激に上昇する（図２（Ａ））。すると、筒内温度（図２（Ｂ））、筒内圧力（図２（Ｃ））も上昇する。そして、点火された火が伝わる前に未燃ガスが自着火して急激な燃焼を起こすと、ノックを発生する（図２（Ｃ））。自着火時間τの逆数（１／τ）が着火遅れを表す。この１／τの時間積分値が１となるクランク角でノックが発生する（図２（Ｄ））。すなわち、∫１／τ・ｄｔがノック発生指標である。
【００２１】
図３を参照してさらに詳しく説明する。なお図３（Ａ）は横軸に筒内温度、縦軸に筒内圧力をとったときの１／τマップである。筒内温度及び筒内圧力がわかれば、このマップから１／τの値が求まる。１／τを筒内圧力、筒内温度の履歴に沿って求めて時間積分して∫１／τ・ｄｔを求める。∫１／τ・ｄｔ＝１となるクランク角でノックが発生すると推定できる。
【００２２】
その推定の実証結果を図３（Ｂ）に示す。図３（Ｂ）の横軸にノック発生クランク角の実測値をとり、一方、縦軸に上記のようにして算出したノック発生クランク角の推定値をとる。すると図３（Ｂ）からも明らかなように推定値が実測値と略一致し、上記推定の正しいことが立証できた。
【００２３】
本件発明者は、ノック発生指標として∫１／τ・ｄｔを導入した。これにより、ノック発生指標を筒内圧力及び筒内温度に帰着させることができたのである。筒内圧力及び筒内温度は筒内ガスの最も基本的な状態量であるので、機関機種によらずに略同一の特性とすることが可能になり、またエンジンの暖機状態や環境条件にかかわらず、精度よくリアルタイムにノック発生時期を推定することが可能になったのである。
【００２４】
図４は本発明による点火時期制御装置の制御の全体を示す制御ブロック図である。
【００２５】
エンジンコントローラ５０は、運転状態５１に基づいて燃焼速度５２１、燃焼期間５２２、基本点火時期５２３を順次算出し、また、ノック発生指標５３を演算し、進角補正／限界（リタード幅）５４を求め、それらよりノック限界点火時期５５を算出する。
【００２６】
続いて、各構成ブロックについてさらに詳しく説明する。なお、基本点火時期算出部５２の具体的な内容については特開２００３−１４８２３６号公報において詳細に説明されているので本件では説明を省略する。
【００２７】
図５はノック発生指標ブロック５３を詳述する図である。
【００２８】
ノック発生指標ブロック５３においては、基本点火時期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を超えるまで、割込演算を要求する。
【００２９】
ノック発生指標演算ブロック５３０については図６を参照しながら説明する。ノック発生指標演算ブロック５３０は、筒内平均温度・圧力演算ブロック５３１と、ノック発生指標演算ブロック５３２とを有する。筒内平均温度・圧力演算ブロック５３１においては、クランク角度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]n+1]    …(1)
θS  ：実燃焼開始クランク角(deg)
θB  ：実燃焼期間(deg)
ａ，ｎ：定数
と記述することができる。上記ａ，ｎは燃焼室形状や点火プラグ位置、および燃焼室内のガス流動特性等、エンジン機種によって定まる。これは実験等によって予め決めておくことができる。この(1)式がWiebe関数である。
【００３８】
上記(1)式をθで微分すると熱発生率を表す関数となり、
ｄＸ／ｄθ
＝ａ(ｎ＋１)/θB・[(θ-θS)/θB]n・exp[-ａ・[(θ-θS)/θB]n+1]    …(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の式について説明する。
【００５０】
冷却損失QLは以下により算出される、
QL = ∫t1t2 h・A・( T - TW )・dt            …(4)
h = 110・d-0.2・P0.8・T-0.53・[C1・cm + C2・(VS・TI)/(PI・VI)・(P - PM)]0.8  …(5)
h  ：熱伝達係数      [ kcal/(m2・h・K) ]
d  ：シリンダボア径      [ m ]
VS  ：シリンダ容積      [ m3 ]
cm ：平均ピストン速度      [ m/s ]
PI  ：吸気弁閉時筒内圧    [ kgf/cm2 ]
VI  ：吸気弁閉時筒内容積    [ m3 ]
TI  ：吸気弁閉時筒内温度    [ K ]
PM  ：モータリング筒内圧    [ kgf/cm2 ]
P  ：筒内圧        [ kgf/cm2 ]
T  ：筒内温度      [ K ]
C1  ：6.18（排気），2.28（圧縮＆膨張）  [−]
C2  ：3.24×10-3      [ m/(s・K) ]
なお、上記でC1、C2は一般に用いられている定数である。
【００５１】
図１６は(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関数定数設定ブロック（熱発生パターン算出手段）
５３１２２  熱発生率算出ブロック
５３１２３  質量燃焼割合算出ブロック（ノック直前筒内ガス質量燃焼割合算出手段）
５３１２４  発熱量／１シリンダ／１サイクル算出ブロック
５３１２５  発熱量／１シリンダ／演算間隔算出ブロック
５３２  ノック発生指標演算ブロック（ノック発生指標算出手段）
５３２２  未燃ガス温度演算ブロック（未燃ガス温度算出手段）
５３２３  瞬時着火遅れ逆数演算ブロック（ノック時期算出手段）
５４  リタード幅算出ブロック（ノック限界点火時期算出手段（点火時期補正量算出手段））
５４１  しきい値設定手段
５５  ノック限界点火時期算出ブロック（ノック限界点火時期算出手段）
【図１】
2005233111.tif 000003
【図２】
2005233111.tif 000004
【図３】
2005233111.tif 000005
【図４】
2005233111.tif 000006
【図５】
2005233111.tif 000007
【図６】
2005233111.tif 000008
【図７】
2005233111.tif 000009
【図８】
2005233111.tif 000010
【図９】
2005233111.tif 000011
【図１０】
2005233111.tif 000012
【図１１】
2005233111.tif 000013
【図１２】
2005233111.tif 000014
【図１３】
2005233111.tif 000015
【図１４】
2005233111.tif 000016
【図１５】
2005233111.tif 000017
【図１６】
2005233111.tif 000018
【図１７】
2005233111.tif 000019
【図１８】
2005233111.tif 000020
【図１９】
2005233111.tif 000021
【図２０】
2005233111.tif 000022
【図２１】
2005233111.tif 000023
【図２２】
2005233111.tif 000024
【図２３】
2005233111.tif 000025
【図２４】
2005233111.tif 000026
【図２５】
2005233111.tif 000027
【図２６】
2005233111.tif 000028