(57)【要約】【課題】可変バルブタイミング機構に与える指令値は変わらないのに、可変バルブ機構に与える指令値によって吸気バルブの開閉タイミングが進遅角した場合でも、トレース点火時期基本値を精度良く算出し得るエンジンの点火時期制御方法を提供する。【解決手段】吸気バルブの作動角は一定のまま吸気バルブの開閉タイミングを指令値に応じて進遅角する可変バルブタイミング機構（２７）と、吸気バルブのバルブリフト最大量及び吸気バルブの作動角を指令値に応じて可変制御する可変バルブ機構（２８）とを備えるエンジンにおいて、エンジンの回転速度、エンジンの負荷、可変バルブタイミング機構（２７）に与える指令値及び可変バルブ機構（２８）に与える指令値に基づいてトレース点火時期基本値を算出するトレース点火時期基本値算出処理手順と、この算出されたトレース点火時期基本値で火花点火を行う火花点火処理手順とをエンジンコントローラ（３１）が含む。【選択図】図１
2007332926.tif 000002
【特許請求の範囲】
【請求項１】
吸気バルブの作動角は一定のまま吸気バルブの開閉タイミングを指令値に応じて進遅角する可変バルブタイミング機構と、
吸気バルブのバルブリフト最大量及び吸気バルブの作動角を指令値に応じて可変制御する可変バルブ機構と
を備えるエンジンにおいて、
エンジンの回転速度、エンジンの負荷、前記可変バルブタイミング機構に与える指令値及び前記可変バルブ機構に与える指令値に基づいてトレース点火時期基本値を算出するトレース点火時期基本値算出処理手順と、
この算出されたトレース点火時期基本値で火花点火を行う火花点火実行処理手順と
を含むことを特徴とするエンジンの点火時期制御方法。
【請求項２】
前記トレース点火時期基本値を、前記可変バルブタイミング機構に与える指令値、吸気バルブ中心角のいずれか一つまたは前記可変バルブ機構に与える指令値、吸気バルブ作動角、吸気バルブリフト最大量のいずれか一つの一次式で近似することを特徴とする請求項１に記載のエンジンの点火時期制御方法。
【請求項３】
前記トレース点火時期基本値を前記可変バルブタイミング機構に与える指令値、吸気バルブ中心角のいずれか一つの一次式で近似する場合にその一次式の傾きと切片とを前記可変バルブ機構に与える指令値、吸気バルブ作動角、吸気バルブリフト最大量のいずれか一つの一次式で近似し、前記トレース点火時期基本値を前記可変バルブ機構に与える指令値、吸気バルブ作動角、吸気バルブリフト最大量のいずれか一つの一次式で近似する場合にその一次式の傾きと切片とを前記可変バルブタイミング機構に与える指令値、吸気バルブ中心角のいずれか一つの一次式で近似することを特徴とする請求項２に記載のエンジンの点火時期制御方法。
【請求項４】
前記一次式の傾きと切片とを、前記一次式で近似した第１のトレース点火時期基本値が実トレース点火時期と所定の適合点で一致するように設定することを特徴とする請求項３に記載のエンジンの点火時期制御方法。
【請求項５】
前記一次式の傾きと切片とを、前記第１のトレース点火時期基本値がさらに前記適合点以外で実トレース点火時期より進角しないよう前記第１のトレース点火時期基本値を遅角側にオフセットして設定することを特徴とする請求項４に記載のエンジンの点火時期制御方法。
【請求項６】
前記適合点以外で実トレース点火時期より進角しないよう遅角側にオフセットして設定している第２のトレース点火時期基本値とトレース点火時期補正量とからトレース点火時期推定値を算出するトレース点火時期推定値算出処理手順と、
このトレース点火時期推定値を前記第２のトレース点火時期基本値に代えて火花点火を行う火花点火実行処理手順と
を含み、
前記適合点で前記トレース点火時期推定値と実トレース点火時期とが一致するように前記トレース点火時期補正量を設定することを特徴とする請求項５に記載のエンジンの点火時期制御方法。
【請求項７】
前記適合点を中心として所定の補正範囲を設け、この補正範囲において前記トレース点火時期推定値が、前記第２のトレース点火時期基本値の直線へと滑らかにつながるように前記適合点以外でのトレース点時期補正量を算出することを特徴とする請求項６に記載のエンジンの点火時期制御方法。
【請求項８】
吸気バルブの作動角は一定のまま吸気バルブの開閉タイミングを指令値に応じて進遅角する可変バルブタイミング機構と、
吸気バルブのバルブリフト最大量及び吸気バルブの作動角を指令値に応じて可変制御する可変バルブ機構と
を備えるエンジンにおいて、
エンジンの回転速度、エンジンの負荷、前記可変バルブタイミング機構に与える指令値及び前記可変バルブ機構に与える指令値に基づいてトレース点火時期基本値を算出するトレース点火時期基本値算出手段と、
この算出されたトレース点火時期基本値で火花点火を行う火花点火実行手段と
を含むことを特徴とするエンジンの点火時期制御装置。
【請求項９】
前記トレース点火時期基本値を、前記可変バルブタイミング機構に与える指令値、吸気バルブ中心角のいずれか一つまたは前記可変バルブ機構に与える指令値、吸気バルブ作動角、吸気バルブリフト最大量のいずれか一つの一次式で近似することを特徴とする請求項８に記載のエンジンの点火時期制御装置。
【請求項１０】
前記トレース点火時期基本値を前記可変バルブタイミング機構に与える指令値、吸気バルブ中心角のいずれか一つの一次式で近似する場合にその一次式の傾きと切片とを前記可変バルブ機構に与える指令値、吸気バルブ作動角、吸気バルブリフト最大量のいずれか一つの一次式で近似し、前記トレース点火時期基本値を前記可変バルブ機構に与える指令値、吸気バルブ作動角、吸気バルブリフト最大量のいずれか一つの一次式で近似する場合にその一次式の傾きと切片とを前記可変バルブタイミング機構に与える指令値、吸気バルブ中心角のいずれか一つの一次式で近似することを特徴とする請求項９に記載のエンジンの点火時期制御装置。
【請求項１１】
前記一次式の傾きと切片とを、前記一次式で近似した第１のトレース点火時期基本値が実トレース点火時期と所定の適合点で一致するように設定することを特徴とする請求項１０に記載のエンジンの点火時期制御装置。
【請求項１２】
前記一次式の傾きと切片とを、前記第１のトレース点火時期基本値がさらに前記適合点以外で実トレース点火時期より進角しないよう前記第１のトレース点火時期基本値を遅角側にオフセットして設定することを特徴とする請求項１１に記載のエンジンの点火時期制御装置。
【請求項１３】
前記適合点以外で実トレース点火時期より進角しないよう遅角側にオフセットして設定している第２のトレース点火時期基本値とトレース点火時期補正量とからトレース点火時期推定値を算出するトレース点火時期推定値算出手段と、
このトレース点火時期推定値を前記第２のトレース点火時期基本値に代えて火花点火を行う火花点火実行手段と
を含み、
前記適合点で前記トレース点火時期推定値と実トレース点火時期とが一致するように前記トレース点火時期補正量を設定することを特徴とする請求項１２に記載のエンジンの点火時期制御装置。
【請求項１４】
前記適合点を中心として所定の補正範囲を設け、この補正範囲において前記トレース点火時期推定値が、前記第２のトレース点火時期基本値の直線へと滑らかにつながるように前記適合点以外でのトレース点時期補正量を算出することを特徴とする請求項１３に記載のエンジンの点火時期制御装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、エンジン（内燃機関）の点火時期制御方法及びエンジンの点火時期制御装置に関する。
【背景技術】
【０００２】
吸気バルブの作動角は一定のまま吸気バルブの開閉タイミング（開時期と閉時期）を指令値に応じて進遅角する可変バルブタイミング機構を備えるエンジンにおいて、吸気バルブ閉時期ＩＶＣが吸気下死点（吸気行程が終了するときにおけるピストンの下死点のこと）であるのときのトレース点火時期をトレース点火時期基本値とし、これに吸気バルブ閉時期ＩＶＣが吸気下死点からズレたとき、このトレース点火時期基本値に吸気バルブ閉時期ＩＶＣの吸気下死点からのズレ量に応じたトレース点火時期補正量を加えることで、トレース点火時期推定値を算出するものがある（特許文献１参照）。
【特許文献１】特開２００５−２３８０６公報
【発明の開示】
【発明が解決しようとする課題】
【０００３】
ところで、上記の可変バルブタイミング機構に加えて、吸気バルブのバルブリフト最大量及び吸気バルブの作動角を指令値に応じて可変制御する可変バルブ機構を備えるエンジンがある。こうしたエンジンでは、可変バルブタイミング機構に与える指令値は変わらないのに、可変バルブ機構に与える指令値によって吸気バルブの開閉タイミング（開時期や閉時期）が進角したり遅角したりするため、この進遅角後の吸気バルブ開閉タイミングに最適なトレース点火時期基本値が、この進遅角前の吸気バルブ開閉タイミングに最適なトレース点火時期基本値と相違することになる。
【０００４】
しかしながら、上記特許文献１の技術には可変バルブ機構について一切記載がないため、可変バルブタイミング機構と可変バルブ機構とを備えるエンジンに特許文献１の技術をそのまま適用した場合において、可変バルブ機構に与える指令値によって吸気バルブの開閉タイミングが進遅角したときに全く対応できず、トレース点火時期基本値の算出に誤差が生じてしまう。
【０００５】
そこで本発明は、可変バルブタイミング機構と可変バルブ機構とを備えるエンジンにおいて、可変バルブタイミング機構に与える指令値は変わらないのに、可変バルブ機構に与える指令値によって吸気バルブの開閉タイミング（開時期や閉時期）が進角したり遅角したりした場合でも、トレース点火時期基本値を精度良く算出し得るエンジンの点火時期制御方法及びエンジンの点火時期制御装置を提供することを目的とする。
【課題を解決するための手段】
【０００６】
本発明は、吸気バルブの作動角は一定のまま吸気バルブの開閉タイミングを指令値に応じて進遅角する可変バルブタイミング機構と、吸気バルブのバルブリフト最大量及び吸気バルブの作動角を指令値に応じて可変制御する可変バルブ機構とを備えるエンジンにおいて、エンジンの回転速度、エンジンの負荷、前記可変バルブタイミング機構に与える指令値及び前記可変バルブ機構に与える指令値に基づいてトレース点火時期基本値を算出し、この算出されたトレース点火時期基本値で火花点火を行うように構成する。
【発明の効果】
【０００７】
本発明によれば、可変バルブタイミング機構と可変バルブ機構とを備えるエンジンにおいて、エンジンの回転速度、エンジンの負荷、可変バルブタイミング機構に与える指令値及び可変バルブ機構に与える指令値に基づいてトレース点火時期基本値を算出し、この算出されたトレース点火時期基本値で火花点火を行うので、可変バルブタイミング機構に与える指令値は変わらないのに、可変バルブ機構に与える指令値によって吸気バルブの開閉タイミング（開時期や閉時期）が進角したり遅角したりした場合でも、トレース点火時期基本値を精度良く算出することができ、トレース点火時期基本値の算出（推定）精度が向上する。
【発明を実施するための最良の形態】
【０００８】
以下、図面に基づき本発明の実施形態について説明する。
【０００９】
図１は、エンジンの点火時期制御方法の実施に直接使用するエンジンの点火時期制御装置の概略構成を示している。
【００１０】
空気は吸気コレクタ２に蓄えられた後、吸気マニホールド３を介して各気筒の燃焼室５に導入される。燃料は各気筒の吸気ポート４に配置された燃料インジェクタ２１より噴射供給される。空気中に噴射された燃料は気化しつつ空気と混合してガス（混合気）を作り、燃焼室５に流入する。この混合気は吸気弁１５が閉じることで燃焼室５内に閉じこめられ、ピストン６の上昇によって圧縮される。
【００１１】
この圧縮混合気に対して高圧火花により点火を行うため、パワートランジスタ内蔵の点火コイルを各気筒に配した電子配電システムの点火装置１１を備える。すなわち、点火装置１１は、バッテリからの電気エネルギーを蓄える点火コイル１３と、点火コイル１３の一次側への通電、遮断を行うパワートランジスタ（図示しない）と、燃焼室５の天井に設けられ点火コイル１３の一次電流の遮断によって点火コイル１３の二次側に発生する高電圧を受けて、火花放電を行う点火プラグ１４とからなっている。
【００１２】
圧縮上死点より少し手前で点火プラグ１４により火花が飛ばされ圧縮混合気に着火されると、火炎が広がりやがて爆発的に燃焼し、この燃焼によるガス圧がピストン６を押し下げる仕事を行う。この仕事はクランクシャフト７の回転力として取り出される。燃焼後のガス（排気）は排気弁１６が開いたときに排気通路８へと排出される。
【００１３】
排気通路８には三元触媒９、１０を備える。三元触媒９、１０は排気の空燃比が理論空燃比を中心とした狭い範囲（ウインドウ）にあるとき、排気中に含まれるＨＣ、ＣＯ、ＮＯｘといった有害三成分を同時に効率よく除去できる。空燃比は吸入空気量と燃料量の比であるので、エンジンの１サイクル（４サイクルエンジンではクランク角で７２０°区間）当たりに燃焼室５に導入される吸入空気量と、燃料インジェクタ２１からの燃料噴射量との比が理論空燃比となるように、エンジンコントローラ３１ではエアフローセンサ３２からの吸入空気流量の信号とクランク角センサ（３３、３４）からの信号に基づいて燃料インジェクタ２１からの燃料噴射量を定めると共に、三元触媒９の上流に設けたＯ2センサ３５からの信号に基づいて空燃比をフィードバック制御している。
【００１４】
吸気コレクタ２の上流には絞り弁２３がスロットルモータ２４により駆動される、いわゆる電子制御スロットル２２を備える。運転者が要求するトルクはアクセルペダル４１の踏み込み量（アクセル開度）に現れるので、エンジンコントローラ３１ではアクセルセンサ４２からの信号に基づいて目標トルクを定め、この目標トルクを実現するための目標空気量を定め、この目標空気量が得られるようにスロットルモータ２４を介して絞り弁２３の開度を制御する。
【００１５】
吸気バルブ１５、排気バルブ１６は、クランクシャフト７を動力源として、各々吸気側カムシャフト２５及び排気側カムシャフト２６に設けられたカムの動作により開閉駆動される。吸気側には、吸気バルブ１５のバルブリフト最大量及び作動角を連続的に可変制御する多節リンク状の機構で構成される可変バルブ機構（ＶＥＬ機構）２８を備える。このＶＥＬ機構２８には吸気バルブ１５のバルブリフト最大量及び作動角を検出するＶＥＬ角度センサ４３が併設されている。
【００１６】
同じく吸気側には、クランクシャフト７と吸気側カムシャフト２５との回転位相差を連続的に可変制御して、吸気バルブ１５の開閉タイミング（開時期ＩＶＯと閉時期ＩＶＣ）を進遅角する可変バルブタイミング機構（ＶＴＣ機構）２７を備える。また、吸気側カムシャフト２５の他端には吸気側カムシャフト２５の回転位置を検出するためのカム角度センサ３４が併設されている。
【００１７】
これらＶＥＬ機構２８及びＶＴＣ機構２７を備える場合のバルブリフト特性を図２に示す。ＶＴＣ機構２７の非作動時かつＶＥＬ機構２８の非作動時のバルブリフト特性が実線であるとすると、ＶＴＣ機構２７の非作動状態でＶＥＬ機構２８に指令値（ＶＥＬ角度）を与えて作動させたときには吸気バルブ１５の作動角及びバルブリフト最大量がＶＥＬ機構２８の非作動時より大きくなるためバルブリフト特性が破線の特性へと移り、吸気バルブ開時期ＩＶＯはＶＥＬ機構２８の非作動時より進角され、吸気バルブ閉時期ＩＶＣはＶＥＬ機構２８の非作動時より遅角される。またＶＥＬ機構２８の非作動状態でＶＴＣ機構２７に指令値（ＶＴＣ角度）を与えて作動させたときには作動角一定のまま吸気バルブ開時期ＩＶＯ及び吸気バルブ閉時期ＩＶＣが進角されるためバルブリフト特性が一点鎖線の特性へと移る。さらに、ＶＴＣ機構２７に指令値を与えて作動させると共にＶＥＬ機構２８にも指令値を与えて作動させたときは二点鎖線の特性へと移る。このようにして、ＶＥＬ機構２８及びＶＴＣ機構２７を組み合わせて用いることで、任意のクランク角位置における吸気バルブの開閉制御が可能となる。これらＶＥＬ機構２８及びＶＴＣ機構２７（可変動弁装置）の具体的な構成は特開２００３−３８７２号公報により公知であるので、その詳しい説明は省略する。
【００１８】
ＶＥＬ機構２８、ＶＴＣ機構２７の各アクチュエータに指令して、吸気バルブ１５のリフト特性（バルブタイミング（開閉時期）や吸気バルブ１５のバルブリフト最大量）を変えると燃焼室５に残留する不活性ガスの量が変化する。燃焼室５内の不活性ガスの量が増えるほどポンピングロスが減って燃費がよくなるので、運転条件によりどのくらいの不活性ガスが燃焼室５内に残留したらよいかを目標吸気バルブ閉時期や目標バルブリフト最大量にして予め定めており、エンジンコントローラ３１ではそのときの運転条件（エンジンの負荷と回転速度）より目標吸気バルブ閉時期と目標バルブリフト最大量とを定め、それら目標値が得られるようにＶＴＣ機構２７及びＶＥＬ機構２８の各アクチュエータを介して吸気バルブ１５の閉時期とバルブリフト最大量とを制御する。
【００１９】
また、ノックセンサ３７からの信号、エンジンの負荷と回転速度の信号がＶＥＬ角度センサ４３からのＶＥＬ角度の信号、カム角度センサ３４からのＶＴＣ角度の信号と共に入力されるエンジンコントローラ３１では、点火コイル１３を介して点火プラグ１４の一次側電流の遮断時期である点火時期を制御する。ここで、ＶＥＬ角度とは、ＶＥＬ機構２８に与える指令値のことで、このＶＥＬ角度がゼロであるときＶＥＬ機構２８は非作動状態にあり、ＶＥＬ角度が正の値で大きくなるほど吸気バルブ１５の作動角及びバルブリフト最大量が大きくなる。また、ＶＴＣ角度とは、ＶＴＣ機構２７に与える指令値のことで、このＶＴＣ角度がゼロであるときＶＴＣ機構２８は非作動状態の最遅角位置にあり、ＶＴＣ角度が正の値で大きくなるほど吸気バルブ開時期ＩＶＯ及び吸気バルブ閉時期ＩＶＣが進角される。
【００２０】
さて、点火時期は、ＭＢＴ（最大トルクの得られる最小進角値）とトレース点火時期とから定まる。ここで、トレース点火時期とは、エンジンにダメージを与えることのない軽いノックを許容している点火時期のことで、トレース点火時期よりも点火時期を進めることは避けなければならない。
【００２１】
そこで、ＶＴＣ機構２７に加えてＶＥＬ機構２８をも備えるエンジンを対象としてトレース点火時期について検討したところを次に述べる。
【００２２】
図３は、エンジン回転速度ＮＲＰＭと充填効率ＩＴＡＣとを所定値に保った状態でＶＴＣ角度とＶＥＬ角度を様々に相違させたときにトレース点火時期基本値がどうなるかを、横軸にＶＴＣ角度を、縦軸にトレース点火時期を採ったグラフにまとめたものである（実験結果）。
【００２３】
図３より、最小ＶＥＬ角度のときＶＴＣ角度を大きくするほどトレース点火時期基本値が進角側に向かい、ＶＴＣ角度の小さな領域であればＶＥＬ角度を大きくするほどトレース点火時期基本値が進角側に向かい、ＶＴＣ角度の大きな領域ではＶＥＬ角度を大きくするほどトレース点火時期基本値が逆の遅角側に向かっている。このため、本発明では実線で示したようにＶＴＣ角度に対するトレース点火時期基本値の特性を一次式（直線）で近似することとする。つまり、ＶＴＣ角度を変数として次式によりトレース点火時期基本値を算出する。
【００２４】
トレー点火時期基本値＝一次式傾き×ＶＴＣ角度＋一次式切片
…（補１）
ただし、一次式切片：ＶＴＣ角度がゼロのときのトレース点火時期基本値、
このようにトレース点火時期基本値を一次式で近似したとき、その一次式は傾きと切片とで特定できるので、一次式傾き（近似式傾き）、一次式切片（近似式切片）を、横軸をＶＥＬ角度として改めてまとめてみたところ、一次式傾きについて図４の特性が、一次式切片について図５の特性が得られた。図４、図５より、本発明ではＶＥＬ角度に対する一次式傾き、ＶＥＬ角度に対する一次式切片の各特性についても一次式で近似することとする。つまり、ＶＥＬ角度を変数として次式により一次式傾き、一次式切片を算出する。
【００２５】
一次式傾き＝一次式傾きの傾き×ＶＥＬ角度＋一次式傾きの切片
…（補２）
一次式切片＝一次式切片の傾き×ＶＥＬ角度＋一次式切片の切片
…（補３）
ただし、一次式傾きの切片：ＶＥＬ角度がゼロのときの一次式傾き、
一次式切片の切片：ＶＥＬ角度がゼロのときの一次式切片、
すると、一次式傾きを近似する一次式、一次式切片を近似する一次式も、傾きと切片とで特定できることとなる。この場合に、図４、図５に示した一次式の傾きと切片（ｋ、ｌ、ｍ、ｎ）はエンジン回転速度ＮＲＰＭ、充填効率ηｃ毎に異なるはずであるから、図４、図５に示した一次式の傾きと切片（一次式傾きの傾きｋ、一次式傾きの切片ｌ、一次式切片の傾きｍ、一次式切片の切片ｎ）つまり、図４に示した一次式傾きについてその傾きと切片（ｋ、ｌ）、図５に示した一次式切片についてその傾きと切片（ｍ、ｎ）はエンジン回転速度ＮＲＰＭと充填効率ηｃのマップとする。
【００２６】
次に、このトレース点火時期基本値の一次式での近似方法が妥当か否かを検討したところを述べる。
【００２７】
図６（Ａ）は上記一次式傾き、一次式切片の傾きと切片（ｋ、ｌ、ｍ、ｎ）を用い、ＶＴＣ角度が最大値に近くかつＶＥＬ角度がほぼ中央の位置をトレース点火時期基本値の適合点（★印参照）としてトレース点火時期基本値を設定したところを破線で書き入れたものであるが、破線で示した第１のトレース点火時期基本値のままでは、○印で囲んだＶＥＬ角度、ＶＴＣ角度のときつまり適合点以外で、第１のトレース点火時期基本値（破線）が実トレース点火時期より進角するためノックが発生する。
【００２８】
そこで、図６（Ｂ）に示したように実トレース点火時期を辿る位置まで、第１のトレース点火時期基本値を遅角側にオフセットすることで（一点鎖線参照）、ノックの発生を抑制できることとなる。しかしながら、適合点以外で実トレース点火時期より進角しないよう遅角側にオフセットして設定しているトレース点火時期基本値、つまり一点鎖線で示した第２のトレース点火時期基本値によれば、今度は★印で示した適合点のＶＴＣ角度のとき、実トレース点火時期とこの第２のトレース点火時期基本値（一点鎖線）との間に点火時期の誤差が発生してしまう。このため、適合点のＶＴＣ角度のときに、第２のトレース点火時期基本値と実トレース点火時期とが一致する補正機能を追加する必要がある。ただし、適合点のＶＴＣ角度のときのみ補正して第２のトレース点火時期基本値と実トレース点火時期とを一致させると、適合点を外れたＶＴＣ角度のとき第２のトレース点火時期基本値と実トレース点火時期との間に点火時期差が生じトルク段差を誘発する可能性がある。
【００２９】
そこで、図６（Ｃ）に示したように適合点のＶＴＣ角度を中心として所定の補正範囲を設け、この補正範囲において適合点のＶＴＣ角度のときに、第２のトレース点火時期基本値と実トレース点火時期とが一致するように滑らかにつなぐトレース点時期補正量を導入する。すなわち、
トレース点火時期補正量＝適合点の実トレース点火時期
−第２のトレース点火時期基本値の、適合点のＶＴＣ角度のときの値
…（補４）
の式により与えられるトレース点火時期補正量が最大となり、後はＶＴＣ角度が適合点より小さい側に（あるいは大きい側に）離れるほど徐々に小さくなって最後は補正範囲の境界でゼロとなるトレース点火時期補正量を算出させる。
【００３０】
そして、第２のトレース点火時期基本値にこのトレース点時期補正量を加算した値をトレース点火時期推定値として、つまり次式によりトレース点火時期推定値を算出する。
【００３１】
トレース点火時期推定値＝第２のトレース点火時期基本値＋トレース点火時期補正量
…（補５）
（補５）式右辺の第２のトレース点火時期基本値の単位は圧縮上死点より進角側に計測した値［ｄｅｇＢＴＤＣ］であるため、正の値のトレース点火時期補正量［ｄｅｇ］を加算することは、第２のトレース点火時期基本値よりも進角されることを意味する。こうして算出されたトレース点火時期推定値は、図６（Ｃ）に示した破線のように、補正範囲内において、第２のトレース点火時期基本値（一点鎖線）から離れて適合点の実トレース時期に向かう直線（破線）となる。
【００３２】
図３では横軸をＶＴＣ角度に採り、図４、図５では横軸にＶＥＬ角を採った場合で説明したが、実はＶＴＣ角度とＶＥＬ角度を入れ換えても、同様となることを確認している。つまり、図３において横軸にＶＥＬ角度を採っても、図３に示したのと同様の特性が、また、図４、図５において横軸にＶＴＣ角度を採っても、図４、図５に示したのと同様の特性が得られるので、この場合にも図６（Ａ）〜図６（Ｃ）で説明したのと同様にしてトレース点火時期推定値の算出方法を考えることができる。
【００３３】
なお、実施形態では、図６（Ａ）〜図６（Ｃ）に示したように、適合点をＶＴＣ角度が最大値に近くかつＶＥＬ角度がほぼ中央の位置としているが、適合点はこの場合に限られるものでない。
【００３４】
このようにして、トレース点火時期推定値の算出方法に対する考え方が定まったので、次には具体的に考える。
【００３５】
上記のＶＴＣ角度、ＶＥＬ角度はＶＴＣ機構２７、ＶＥＬ機構２８に対する指令値のため、実際の吸気バルブ１５のリフト特性を直接表している値でないので、図１９に示したように、バルブリフト最大量の位置のクランク角より吸気下死点（ＢＤＣ）までのクランク角区間を吸気バルブ１５の「中心角」として定義すると、吸気バルブ１５の作動角及びバルブリフト最大量はＶＥＬ角度により定まり、吸気バルブ１５の中心角はＶＴＣ角度により定まるので、上記の（補１）式〜（補３）式は次のように書き直すことができる。
【００３６】
トレー点火時期基本値＝近似式傾き×中心角＋近似式切片
…（補１’）
近似式傾き＝近似式傾きの傾き×作動角（またはリフト）＋近似式傾きの切片
…（補２’）
近似式切片＝近似式切片の傾き×作動角（またはリフト）＋近似式切片の切片
…（補３’）
（補１’）式、（補２’）式右辺の作動角は、吸気バルブ１５の開時期ＩＶＯより閉時期ＩＶＣまでのクランク角区間のことである。また、吸気バルブ１５のバルブリフト最大量は吸気バルブ１５のバルブリフトが最大となる位置のバルブリフト量であるが、（補１’）式、（補２’）式を含めて、以下ではこれを単に吸気バルブ１５の「リフト」という。
【００３７】
また、（補１’）式〜（補３’）式において、作動角（またはリフト）と中心角とを入れ換えた次式が考えられる。
【００３８】
トレー点火時期基本値＝近似式傾き×作動角（またはリフト）＋近似式切片
…（補１”）
近似式傾き＝近似式傾きの傾き×中心角＋近似式傾きの切片
…（補２”）
近似式切片＝近似式切片の傾き×中心角＋近似式切片の切片
…（補３”）
これら（補１’）式〜（補３’）式や（補１”）式〜（補３”）式では３つの傾きと３つの切片が出てくるので、それぞれを区別するため、（補１’）式〜（補３’）式を次の（補６）式〜（補８）式のように、また（補１”）式〜（補３”）式を次の（補６’）式〜（補８’）式のように書き換える。
【００３９】
トレー点火時期基本値＝作動角感度×作動角＋作動角切片
…（補６）
作動角感度＝作動角感度中心角感度×中心角＋作動角感度中心角切片
…（補７）
作動角切片＝作動角切片中心角感度×中心角＋作動角切片中心角切片
…（補８）
トレー点火時期基本値＝リフト感度×リフト＋リフト切片
…（補６’）
リフト感度＝リフト感度中心角感度×中心角＋リフト感度中心角切片
…（補７’）
リフト切片＝リフト切片中心角感度×中心角＋リフト切片中心角切片
…（補８’）
後述するように、（補６）式〜（補８）式に従う場合が第１実施形態、（補６’）式〜（補８’）式に従う場合が第２実施形態である。
【００４０】
エンジンコントローラ３１で実行されるこの制御について以下のフローチャートに基づいて詳述する。
【００４１】
図７〜図１０は第１実施形態のフローチャートである。このうち、図７はトレース点火時期推定値を算出するためのもので、Ｒｅｆ信号の入力毎に実行する。
【００４２】
ステップ０１では、クランク角センサ（３３、３４）により検出されるエンジン回転数ＮＲＰＭ［ｒｐｍ］、エアフローセンサ３２により検出される吸入空気量ＭＡ［ｋｇ／ｓ］を読み込む。この吸入空気量ＭＡからは次式により充填効率を算出する。
【００４３】
ＩＴＡＣ＝ＭＡ／ＭＡＭＸ                      …（補９）
ただし、ＭＡＭＸ：充填効率最大値の吸入空気量［ｋｇ／ｓ］、実機による適合値              とする
ステップ０２では、基本トレース点火時期推定値ＢＥＳＴＴＲＡＣＥ１［ｄｅｇＢＴＤＣ］を算出する。この値の算出については図８により後述する。
【００４４】
ステップ０３ではノックセンサ３７の出力に基づいてノック補正量ＫＮＫＲＥＴ［ｄｅｇ］を算出する。例えば、ノックを判定したときにはクランク角で１ｄｅｇのリタード量をノック補正量として、またノック判定しないときには一定時間経過ごとにクランク角で１ｄｅｇ進角するような進角量をノック補正量として算出する。この場合、ノックを判定したときのノック補正量ＫＮＫＲＥＴには負の符号を、ノック判定しないときのノック補正量ＫＮＫＲＥＴには正の符号を付ける。
【００４５】
ステップ０４では、基本トレース点火時期推定値ＢＥＳＴＴＲＡＣＥ１にノック補正量ＫＮＫＲＥＴを加算した値をトレース点火時期推定値ＥＳＴＴＲＡＣＥ１として、つまり次式によりトレース点火時期推定値ＥＳＴＴＲＡＣＥ１［ｄｅｇＢＴＤＣ］を算出する。
【００４６】
ＥＳＴＴＲＡＣＥ１＝ＢＥＳＴＴＲＡＣＥ１＋ＫＮＫＲＥＴ  …（２）
このようにして算出したトレース点火時期推定値ＥＳＴＴＲＡＣＥ１は、図示しないフローにおいて、ＭＢＴの得られる基本点火時期ＭＢＴＣＡＬと比較され、小さい方の値が最終的な点火時期ＡＤＶ［ｄｅｇＢＴＤＣ］として選択される。点火時期の単位は、圧縮上死点を基準として進角側に計測したクランク角位置であるので、小さい方を選択するとは、トレース点火時期推定値ＥＳＴＴＲＡＣＥ１と基本点火時期ＭＢＴＣＡＬのうち遅角側の値を選択することである。そして、選択された値が点火時期指令値として点火レジスタに移され、実際のクランク角がこの点火時期指令値と一致したタイミングでエンジンコントローラ３１より一次電流を遮断する点火信号がパワートランジスタに出力される。
【００４７】
基本トレース点火時期推定値ＢＥＳＴＴＲＡＣＥ１が選択されているときにノックが生じたときには、次回の点火時に基本トレース点火時期推定値ＢＥＳＴＴＲＡＣＥ１よりノック補正量ＫＮＫＲＥＴの分だけ遅角された値がトレース点火時期推定値ＥＳＴＴＲＡＣＥ１となる。
【００４８】
図８は基本トレース点火時期推定値を算出するためのもので、図７のステップ０２のサブルーチンである。
【００４９】
ステップ０２１では、カム角度センサ３４により検出される実ＶＴＣ角度（ＲＥＶＴＣ）と、ＶＥＬ角度センサ４３により検出される実ＶＥＬ角度（ＲＥＶＥＬ）を読み込む。
【００５０】
ステップ０２２では、実ＶＥＬ作動角ＲＥＶＥＬから図１１を内容とするテーブルを参照することにより、吸気バルブ１５の作動角ＲＥＥＶＥＮＴを算出する。図１１の作動角のテーブル特性はＶＥＬ機構２８の仕様により決まる。吸気バルブ１５の作動角は図１９を用いて説明したように吸気バルブ１５の開時期ＩＶＯより閉時期ＩＶＣまでのクランク角区間のことである。
【００５１】
ステップ０２２ではまた、実ＶＴＣ角度ＲＥＶＴＣから次式により吸気バルブ１５の中心角ＲＥＴＩＭＩＮＧを算出する。
【００５２】
ＲＥＴＩＭＩＮＧ＝ＲＥＴＩＭＩＮＧ０＋ＲＥＶＴＣ        …（３）
ただし、ＲＥＴＩＭＩＮＧ０：吸気バルブの中心角の初期位置、
（３）式右辺の吸気バルブの中心角の初期位置ＲＥＴＩＭＩＮＧ０は、ＶＴＣ機構２７の非作動時の吸気バルブの中心角のことで、エンジンジオメトリーから決まる。吸気バルブ１５の中心角は、図１９を用いて説明したように、バルブリフト最大量となる位置のクランク角より吸気下死点（ＢＤＣ）までのクランク角区間のことである。
【００５３】
ステップ０２３では、基本トレース点火時期推定基本値ＢＢＥＳＴＴＲＡＣＥ１［ｄｅｇＢＴＤＣ］を、またステップ０２４では、トレース点火時期誤差補正量ＴＲＡＣＥＥＲＲ１［ｄｅｇ］を算出し、ステップ０２５で基本トレース点火時期推定基本値ＢＢＥＳＴＴＲＡＣＥ１にトレース点火時期誤差補正量ＴＲＡＣＥＥＲＲ１を加算した値（トレース点火時期誤差補正量ＴＲＡＣＥＥＲＲ１だけ進角側の値）を基本トレース点火時期推定値ＢＥＳＴＴＲＡＣＥ１［ｄｅｇＢＴＤＣ］として、つまり次式により基本トレース点火時期推定値ＢＥＳＴＴＲＡＣＥ１を算出する。（４）式は上記（補５）式に対応する。
ＢＥＳＴＴＲＡＣＥ１＝ＢＢＥＳＴＴＲＡＣＥ１＋ＴＲＡＣＥＥＲＲ１
…（４）
（４）式は図６（Ｃ）に示したように適合点のＶＴＣ角度を中心とする所定の補正範囲で破線の特性を求める操作に相当する。（４）式右辺の基本トレース点火時期推定基本値ＢＢＥＳＴＴＲＡＣＥ１、トレース点火時期誤差補正量ＴＲＡＣＥＥＲＲ１の算出については後述する。
【００５４】
図９は基本トレース点火時期推定基本値ＢＢＥＳＴＴＲＡＣＥ１を算出するためのもので、図８のステップ０２３のサブルーチンである。
【００５５】
ステップ０２３１ではエンジン回転速度ＮＲＰＭと充填効率ＩＴＡＣから図１２を内容とするマップを参照することにより基本トレース点火時期推定基本値の作動角感度中心角感度ＥＶＳＬＴＩＭＳＬを、ステップ０２３２ではエンジン回転速度ＮＲＰＭと充填効率ＩＴＡＣから図１３を内容とするマップを参照することにより基本トレース点火時期推定基本値の作動角感度中心角切片ＥＶＳＬＴＩＭＩＮをそれぞれ算出し、ステップ０２３３で次式により基本トレース点火時期推定基本値の作動角感度を算出する。（５）式は上記（補７）式と同じである。
【００５６】
ＥＶＳＬ＝ＥＶＳＬＴＩＭＳＬ×ＲＥＴＩＭＩＮＧ＋ＥＶＳＬＴＩＭＩＮ
…（５）
ただし、ＥＶＳＬ          ：作動角感度、
ＥＶＳＬＴＩＭＳＬ：作動角感度中心角感度、
ＲＥＴＩＭＩＮＧ  ：吸気バルブの中心角、
ＥＶＳＬＴＩＭＩＮ：作動角感度中心角切片、
ステップ０２３４ではエンジン回転速度ＮＲＰＭと充填効率ＩＴＡＣから図１４を内容とするマップを参照することにより基本トレース点火時期推定基本値の作動角切片中心角感度ＥＶＩＮＴＩＭＳＬを、ステップ０２３５ではエンジン回転速度ＮＲＰＭと充填効率ＩＴＡＣから図１５を内容とするマップを参照することにより基本トレース点火時期推定基本値の作動角切片中心角切片ＥＶＩＮＴＩＭＩＮをそれぞれ算出し、ステップ０２３６で次式により基本トレース点火時期推定基本値の作動角切片を算出する。（６）式は上記（補８）式と同じである。
【００５７】
ＥＶＩＮ＝ＥＶＩＮＴＩＭＳＬ×ＲＥＴＩＭＩＮＧ＋ＥＶＩＮＴＩＭＩＮ
…（６）
ただし、ＥＶＩＮ          ：作動角切片、
ＥＶＩＮＴＩＭＳＬ：作動角切片中心角感度、
ＲＥＴＩＭＩＮＧ  ：吸気バルブの中心角、
ＥＶＩＮＴＩＭＩＮ：作動角切片中心角切片、
ステップ０２３７では、次式により基本トレース点火時期推定基本値を算出する。（７）式は上記（補６）式と同じである。
【００５８】
ＢＢＥＳＴＴＲＡＣＥ１＝ＥＶＳＬ×ＲＥＥＶＥＮＴ＋ＥＶＩＮ
…（７）
ただし、ＢＢＥＳＴＴＲＡＣＥ１：基本トレース点火時期推定基本値、
ＥＶＳＬ              ：作動角感度、
ＲＥＥＶＥＮＴ        ：吸気バルブの作動角、
ＥＶＩＮ              ：作動角切片、
ここで、図１２に示した作動角感度中心角感度ＥＶＳＬＴＩＭＳＬ、図１３に示した作動角感度中心角切片ＥＶＳＬＴＩＭＩＮ、図１４に示した作動角切片中心角感度ＥＶＩＮＴＩＭＳＬ、図１５に示した作動角切片中心角切片ＥＶＩＮＴＩＭＩＮの各マップ特性はエンジン回転速度ＮＲＰＭ、充填効率ＩＴＡＣ、ＶＥＬ角度、ＶＴＣ角度が変化したとき、基本トレース点火時期推定基本値ＢＢＥＳＴＴＲＡＣＥが実トレースノック点火時期と同等もしくはリタードするよう予め適合しておく。
【００５９】
図１０はトレース点火時期誤差補正量ＴＲＡＣＥＥＲＲ１を算出するためのもので、図８のステップ０２４のサブルーチンである。
【００６０】
ステップ０２４１ではエンジン回転速度ＮＲＰＭと充填効率ＩＴＡＣから図１６を内容とするマップを参照することにより吸気バルブ１５の目標作動角ＴＧＥＶＥＮＴを、ステップ０２４２ではエンジン回転速度ＮＲＰＭと充填効率ＩＴＡＣから図１７を内容とするマップを参照することにより吸気バルブ１５の目標中心角ＴＧＴＩＭＩＮＧをそれぞれ算出する。図１６に示した目標作動角ＴＧＥＶＥＮＴ、図１７に示した目標中心角ＴＧＴＩＭＩＮＧの各マップ特性は、適合したエンジン回転速度ＮＲＰＭ、充填効率ＩＴＡＣの状態でベストな状態となるよう予め適合しておく。ここでベストな状態とは要求により定まるものである。例えば燃費（あるいはエンジン出力）が最適となるように目標作動角ＴＧＥＶＥＮＴ及び目標中心角ＴＧＴＩＭＩＮＧが定められる。
【００６１】
ステップ０２４３では、そのときの吸気バルブ１５の実作動角ＲＥＥＶＥＮＴの目標作動角ＴＧＥＶＥＮＴからの誤差（ＴＧＥＶＥＮＴ−ＲＥＥＶＥＮＴ）と、そのときの吸気バルブ１５の実中心角ＲＥＴＩＭＩＮＧの目標中心角ＴＧＴＩＭＩＮＧからの誤差（ＴＧＴＩＭＩＮＧ−ＲＥＴＩＭＩＮＧ）とを用いて次式により作動角・中心角の誤差ＥＶＴＩＭＥＲＲを算出する。
【００６２】
ＥＶＴＩＭＥＲＲ＝（（ＴＧＥＶＥＮＴ−ＲＥＥＶＥＮＴ）／ａ）＾２
＋（（ＴＧＴＩＭＩＮＧ−ＲＥＴＩＭＩＮＧ）／ｂ）＾２
…（８）
ただし、ａ、ｂ：定数、
定数ａ，ｂはステップ０２５で算出する基本トレース点火時期推定値ＢＥＳＴＴＲＡＣＥ１が実トレースノック点火時期と同等もしくはリタードするよう予め適合しておく。
【００６３】
ステップ０２４４では、この作動角・中心角誤差ＥＶＴＩＭＥＲＲと作動角・中心角誤差の閾値ＴＨ１を比較する。ここで、閾値ＴＨ１は図６（Ｃ）に示したように適合点よりいずれの範囲で、基本トレース点火時期推定基本値ＢＢＥＳＴＴＲＡＣＥ１を補正するかを定める値（一定値）で、予め適合しておく。
【００６４】
従って、作動角・中心角誤差ＥＶＴＩＭＥＲＲが閾値ＴＨ１より大きいときには補正範囲にないので、ステップ０２４５に進んでトレース点火時期誤差補正量ＴＲＡＣＥＥＲＲ１＝０とする。
【００６５】
一方、作動角・中心角誤差ＥＶＴＩＭＥＲＲが閾値ＴＨ１以下であるときには適合点を中心とする所定の補正範囲にあると判断してステップ０２４６に進み、エンジン回転速度ＮＲＰＭと充填効率ＩＴＡＣから図１８を内容とするマップを参照することにより、トレース点火時期推定誤差補正量基本値ＢＴＲＡＣＥＥＲＲ１を算出する。このトレース点火時期推定誤差補正量基本値ＢＴＲＡＣＥＥＲＲ１は、上記（補４）式の値を求めるものである。
【００６６】
ステップ０２４７では、次式によりトレース点火時期誤差補正量ＴＲＡＣＥＥＲＲ１を算出する。
【００６７】
ＴＲＡＣＥＥＲＲ１＝ＢＴＲＡＣＥＥＲＲ１
×（ＴＨ１−ＥＶＴＩＭＥＲＲ）／ＴＨ１
…（９）
（９）式右辺の（ＥＶＴＩＭＥＲＲ／ＴＨ１）は１以下の正の値である。吸気バルブ１５の実作動角ＲＥＥＶＥＮＴが目標作動角ＴＧＥＶＥＮＴと一致し、吸気バルブ１５の実中心角ＲＥＴＩＭＩＮＧも目標中心角ＴＧＴＩＭＩＮＧと一致するとき作動角・中心角誤差ＥＶＴＩＭＥＲＲ＝０となり、このとき、（９）式よりトレース点火時期推定誤差補正量基本値ＢＴＲＡＣＥＥＲＲ１がそのままトレース点火時期誤差補正量ＴＲＡＣＥＥＲＲ１となる。また、吸気バルブ１５の実作動角ＲＥＥＶＥＮＴが目標作動角ＴＧＥＶＥＮＴから離れかつ吸気バルブ１５の実中心角ＲＥＴＩＭＩＮＧが目標中心角ＴＧＴＩＭＩＮＧから離れるほど作動角・中心角誤差ＥＶＴＩＭＥＲＲが大きくなって閾値ＴＨ１に近づいて行き、作動角・中心角誤差ＥＶＴＩＭＥＲＲが閾値ＴＨ１と一致したとき、トレース点火時期誤差補正量ＴＲＡＣＥＥＲＲ１はゼロとなる。すなわち、（９）式は図６（Ｃ）に示したように適合点を中心とする所定の補正範囲で一点鎖線の特性から破線の特性へと進角側に移すための補正量を求める操作である。
【００６８】
上記のトレース点火時期推定誤差補正量基本値ＢＴＲＡＣＥＥＲＲ１の特性は、吸気バルブ１５の実作動角ＲＥＥＶＥＮＴが適合点の吸気バルブ１５の作動角（つまり目標作動角ＴＧＥＶＥＮＴ）に近づいてゆくほどステップ０２５で算出する基本トレース点火時期推定値ＢＥＳＴＴＲＡＣＥ１が、適合点の実トレースノック点火時期に近づいて行くように、つまり図６（Ｃ）に示した破線の特性が得られるように予め適合しておく。
【００６９】
ここで、本実施形態の作用効果を説明する。
【００７０】
本実施形態（請求項１、８に記載の発明）によれば、ＶＴＣ機構２７とＶＥＬ機構２８とを備えるエンジンにおいて、エンジンの回転速度ＮＲＰＭ、エンジンの負荷、ＶＴＣ機構２７に与える指令値及びＶＥＬ機構２８に与える指令値に基づいて基本トレース点火時期推定基本値ＢＢＥＳＴＴＲＡＣＥ１（トレース点火時期基本値）を算出し（図８のステップ０２３、図９参照）、この算出された基本トレース点火時期推定基本値ＢＢＥＳＴＴＲＡＣＥ１（トレース点火時期基本値）で火花点火を行うので、ＶＴＣ角度（ＶＴＣ機構２７に与える指令値）は変わらないのに、ＶＥＬ角度（ＶＥＬ機構２８に与える指令値）によって吸気バルブ１５の開時期や閉時期（開閉タイミング）が進角したり遅角したりした場合でも、基本トレース点火時期推定基本値ＢＢＥＳＴＴＲＡＣＥ１（トレース点火時期基本値）を精度良く算出することができ、基本トレース点火時期推定基本値ＢＢＥＳＴＴＲＡＣＥ１（トレース点火時期基本値）の算出（推定）精度が向上する。
【００７１】
本実施形態（請求項２、９に記載の発明）によれば、基本トレース点火時期推定基本値ＢＢＥＳＴＴＲＡＣＥ１（トレース点火時期基本値）を、吸気バルブ１５の作動角ＲＥＥＶＥＮＴの一次式で近似するので（図９のステップ０２３７参照）、吸気バルブ１５の作動角、ＶＥＬ角度（ＶＥＬ機構２８に与える指令値）が相違する場合でも、基本トレース点火時期推定基本値ＢＢＥＳＴＴＲＡＣＥ１（トレース点火時期基本値）を容易に算出できる。
【００７２】
本実施形態（請求項３、１０に記載の発明）によれば、基本トレース点火時期推定基本値ＢＢＥＳＴＴＲＡＣＥ１（トレース点火時期基本値）を吸気バルブ１５の作動角ＲＥＥＶＥＮＴの一次式で近似する場合に、基本トレース点火時期推定基本値の作動角感度ＥＶＳＬ（一次式の傾き）と基本トレース点火時期推定基本値の作動角切片ＥＶＩＮ（一次式の切片）とを吸気バルブ１５の中心角ＲＥＴＩＭＩＮＧの一次式で近似するので（図９のステップ０２３３、０２３６参照）、吸気バルブ１５の中心角、ＶＴＣ角度（ＶＴＣ機構２７に与える指令値）が相違する場合でも、基本トレース点火時期推定基本値の作動角感度ＥＶＳＬ（一次式の傾き）と基本トレース点火時期推定基本値の作動角切片ＥＶＩＮ（一次式の切片）とを容易に算出できる。
【００７３】
基本トレース点火時期推定基本値の作動角感度ＥＶＳＬ（一次式の傾き）と基本トレース点火時期推定基本値の作動角切片ＥＶＩＮ（一次式の切片）とを、基本トレース点火時期推定基本値ＢＢＥＳＴＴＲＡＣＥ１（一次式で近似した第１のトレース点火時期基本値）が所定の適合点で実トレース点火時期と一致するように設定したとき、基本トレース点火時期推定基本値ＢＢＥＳＴＴＲＡＣＥ１（第１のトレース点火時期基本値）が適合点以外では実トレース点火時期より進角してノックが発生することがあり得るのであるが、本実施形態（請求項５、１２に記載の発明）によれば、基本トレース点火時期推定基本値の作動角感度ＥＶＳＬ（一次式の傾き）と基本トレース点火時期推定基本値の作動角切片ＥＶＩＮ（一次式の切片）とを、基本トレース点火時期推定基本値ＢＢＥＳＴＴＲＡＣＥ１（第１のトレース点火時期基本値）がさらに適合点以外で実トレース点火時期より進角しないよう基本トレース点火時期推定基本値ＢＢＥＳＴＴＲＡＣＥ１（第１のトレース点火時期基本値）を遅角側にオフセットして設定するので、適合点以外でのノックの発生を回避できる。
【００７４】
本実施形態（請求項６、１３に記載の発明）によれば、基本トレース点火時期推定基本値ＢＢＥＳＴＴＲＡＣＥ１（適合点以外で実トレース点火時期より進角しないよう遅角側にオフセットして設定している第２のトレース点火時期基本値）とトレース点火時期誤差補正量ＴＲＡＣＥＥＲＲ１（トレース点火時期補正量）とから基本トレース点火時期推定値ＢＥＳＴＴＲＡＣＥ１（トレース点火時期推定値）を算出し（図８のステップ０２５参照）、この基本トレース点火時期推定値ＢＥＳＴＴＲＡＣＥ１（トレース点火時期推定値）を基本トレース点火時期推定基本値ＢＢＥＳＴＴＲＡＣＥ１（第２のトレース点火時期基本値）に代えて火花点火を行うと共に、適合点で基本トレース点火時期推定値ＢＥＳＴＴＲＡＣＥ１（トレース点火時期推定値）と実トレース点火時期とが一致するようにトレース点火時期推定誤差補正量基本値ＢＴＲＡＣＥＥＲＲ１（トレース点火時期補正量）を設定するので（図１０のステップ０２４６参照）、適合点で基本トレース点火時期推定値ＢＥＳＴＴＲＡＣＥ１（トレース点火時期推定値）と実トレース点火時期との間に点火時期誤差が発生してしまうことを避けることができる。
【００７５】
適合点でのみ基本トレース点火時期推定値ＢＥＳＴＴＲＡＣＥ１（トレース点火時期推定値）と実トレース点火時期とが一致するように基本トレース点火時期推定値ＢＥＳＴＴＲＡＣＥ１（トレース点火時期補正量）を設定したのでは、適合点とそれ以外とでトレース点火時期推定値と実トレース点火時期との間に点火時期差が生じトルク段差を誘発する可能性があるのであるが、本実施形態（請求項７、１４に記載の発明）によれば、適合点を中心として所定の補正範囲を設け、この補正範囲において基本トレース点火時期推定値ＢＥＳＴＴＲＡＣＥ１（トレース点火時期推定値）が、基本トレース点火時期推定基本値ＢＢＥＳＴＴＲＡＣＥ１（第２のトレース点火時期基本値の直線）へと滑らかにつながるように、トレース点火時期誤差補正量ＴＲＡＣＥＥＲＲ１（適合点以外でのトレース点時期補正量）を算出するので（図１０のステップ０２４７参照）、こうしたトルク段差を誘発することがない。
【００７６】
次に、図２０〜図２３は第２実施形態のフローチャートで、図２０、図２１、図２２、図２３はそれぞれ図７、図８、図９、図１０に対応する。
【００７７】
ここでは、第１実施形態との違いを主に説明する。
【００７８】
図２０においてステップ１１、１４で基本トレース点火時期推定値ＢＥＳＴＴＲＡＣＥ２、トレース点火時期推定値ＥＳＴＴＲＡＣＥ２を算出するが、内容的には同様の値を算出している。ただし、第１実施形態で算出する値であることを表すために第１実施形態では記号の最後に「１」を付しているので、第２実施形態で算出する値であることを表すために第２実施形態では記号の最後を「１」に代えて「２」を付している。この「２」は図２１、図２２、図２３においても使用している。
【００７９】
図２１（図２０ステップ１２のサブルーチン）においてステップ１２２では実ＶＥＬ角度ＲＥＶＥＬから図２４を内容とするテーブルを参照することにより、吸気バルブ１５のリフトＲＥＬＩＦＴを算出する。図２４のリフトのテーブル特性はＶＥＬ機構２８の仕様により決まる。
【００８０】
ステップ１２３では、基本トレース点火時期推定基本値ＢＢＥＳＴＴＲＡＣＥ２［ｄｅｇＢＴＤＣ］を、またステップ１２４では、トレース点火時期誤差補正量ＴＲＡＣＥＥＲＲ２［ｄｅｇ］を算出し、ステップ１２５で基本トレース点火時期推定基本値ＢＢＥＳＴＴＲＡＣＥ２にトレース点火時期誤差補正量ＴＲＡＣＥＥＲＲ２を加算した値（トレース点火時期誤差補正量ＴＲＡＣＥＥＲＲ２だけ進角側の値）を基本トレース点火時期推定値ＢＥＳＴＴＲＡＣＥ２［ｄｅｇＢＴＤＣ］として、つまり次式により基本トレース点火時期推定値ＢＥＳＴＴＲＡＣＥ２を算出する。（１０）式は上記（補５）式に対応する。
ＢＥＳＴＴＲＡＣＥ２＝ＢＢＥＳＴＴＲＡＣＥ２＋ＴＲＡＣＥＥＲＲ２
…（１０）
（１０）式は図６（Ｃ）に示したように適合点のＶＴＣ角度を中心とする所定の補正範囲で破線を求める操作に相当する。（１０）式右辺の基本トレース点火時期推定基本値ＢＢＥＳＴＴＲＡＣＥ２、トレース点火時期誤差補正量ＴＲＡＣＥＥＲＲ２の算出については後述する。
【００８１】
図２２（図２１のステップ１２３のサブルーチン）において、ステップ１２３１ではエンジン回転速度ＮＲＰＭと充填効率ＩＴＡＣから図２５を内容とするマップを参照することにより基本トレース点火時期推定基本値のリフト感度中心角感度ＬＩＳＬＴＩＭＳＬを、ステップ１２３２ではエンジン回転速度ＮＲＰＭと充填効率ＩＴＡＣから図２６を内容とするマップを参照することにより基本トレース点火時期推定基本値のリフト感度中心角切片ＬＩＳＬＴＩＭＩＮをそれぞれ算出し、ステップ１２３３で次式により基本トレース点火時期推定基本値のリフト感度を算出する。（１１）式は上記（補７’）式と同じである。
【００８２】
ＬＩＳＬ＝ＬＩＳＬＴＩＭＳＬ×ＲＥＴＩＭＩＮＧ＋ＬＩＳＬＴＩＭＩＮ
…（１１）
ただし、ＬＩＳＬ          ：リフト感度、
ＬＩＳＬＴＩＭＳＬ：リフト感度中心角感度、
ＲＥＴＩＭＩＮＧ  ：吸気バルブの中心角、
ＬＩＳＬＴＩＭＩＮ：リフト感度中心角切片、
ステップ１２３４ではエンジン回転速度ＮＲＰＭと充填効率ＩＴＡＣから図２７を内容とするマップを参照することにより基本トレース点火時期推定基本値のリフト切片中心角感度ＬＩＩＮＴＩＭＳＬを、ステップ１２３５ではエンジン回転速度ＮＲＰＭと充填効率ＩＴＡＣから図２８を内容とするマップを参照することにより基本トレース点火時期推定基本値のリフト切片中心角切片ＬＩＩＮＴＩＭＩＮをそれぞれ算出し、ステップ１２３６で次式により基本トレース点火時期推定基本値のリフト切片を算出する。（１２）式は上記（補８’）式と同じである。
【００８３】
ＬＩＩＮ＝ＬＩＩＮＴＩＭＳＬ×ＲＥＴＩＭＩＮＧ＋ＬＩＩＮＴＩＭＩＮ
…（１２）
ただし、ＬＩＩＮ          ：リフト切片、
ＬＩＩＮＴＩＭＳＬ：リフト切片中心角感度、
ＲＥＴＩＭＩＮＧ  ：吸気バルブの中心角、
ＬＩＩＮＴＩＭＩＮ：リフト切片中心角切片、
ステップ１２３７では、次式により基本トレース点火時期推定基本値を算出する。（１３）式は上記（補６’）式と同じである。
【００８４】
ＢＢＥＳＴＴＲＡＣＥ２＝ＬＩＳＬ×ＲＥＬＩＦＴ＋ＬＩＩＮ
…（１３）
ただし、ＢＢＥＳＴＴＲＡＣＥ２：基本トレース点火時期推定基本値、
ＬＩＳＬ              ：リフト感度、
ＲＥＬＩＦＴ          ：吸気バルブのリフト、
ＬＩＩＮ              ：リフト作動角切片、
ここで、図２５に示したリフト感度中心角感度ＬＩＳＬＴＩＭＳＬ、図２６に示したリフト感度中心角切片ＬＩＳＬＴＩＭＩＮ、図２７に示したリフト切片中心角感度ＬＩＩＮＴＩＭＳＬ、図２８に示したリフト切片中心角切片ＬＩＩＮＴＩＭＩＮの各マップ特性はエンジン回転速度ＮＲＰＭ、充填効率ＩＴＡＣ、ＶＥＬ角度、ＶＴＣ角度が変化したとき、基本トレース点火時期推定基本値ＢＢＥＳＴＴＲＡＣＥが実トレースノック点火時期と同等もしくはリタードするよう予め適合しておく。
【００８５】
図２３（図２１のステップ１２４のサブルーチンの）においてステップ１２４１ではエンジン回転速度ＮＲＰＭと充填効率ＩＴＡＣから図２９を内容とするマップを参照することにより吸気バルブ１５の目標リフトＴＧＬＩＦＴを、ステップ１２４２ではエンジン回転速度ＮＲＰＭと充填効率ＩＴＡＣから図１７を内容とするマップを参照することにより吸気バルブ１５の目標中心角ＴＧＴＩＭＩＮＧをそれぞれ算出する。図２９に示した目標リフトＴＧＬＩＦＴ、図１７に示した目標中心角ＴＧＴＩＭＩＮＧの各マップ特性は、適合したエンジン回転速度ＮＲＰＭ、充填効率ＩＴＡＣの状態でベストな状態となるよう予め適合しておく。ここでベストな状態とは要求により定まるものである。例えば燃費（あるいはエンジン出力）が最適となるように目標リフトＴＧＬＩＦＴ及び目標中心角ＴＧＴＩＭＩＮＧが定められる。
【００８６】
ステップ１２４３では、そのときの吸気バルブ１５の実リフトＲＥＬＩＦＴの目標リフトＴＧＬＩＦＴからの誤差（ＴＧＬＩＦＴ−ＲＥＬＩＦＴ）と、そのときの吸気バルブ１５の実中心角ＲＥＴＩＭＩＮＧの目標中心角ＴＧＴＩＭＩＮＧからの誤差（ＴＧＴＩＭＩＮＧ−ＲＥＴＩＭＩＮＧ）とを用いて次式によりリフト・中心角の誤差ＬＩＴＩＭＥＲＲを算出する。
【００８７】
ＬＩＴＩＭＥＲＲ＝（（ＴＧＬＩＦＴ−ＲＥＬＩＦＴ）／ｃ）＾２
＋（（ＴＧＴＩＭＩＮＧ−ＲＥＴＩＭＩＮＧ）／ｄ）＾２
…（１４）
ただし、ｃ、ｄ：定数、
定数ｃ，ｄはステップ１２５で算出する基本トレース点火時期推定値ＢＥＳＴＴＲＡＣＥ２が実トレースノック点火時期と同等もしくはリタードするよう予め適合しておく。
【００８８】
ステップ１２４４では、このリフト・中心角誤差ＬＩＴＩＭＥＲＲとリフト・中心角誤差の閾値ＴＨ２を比較する。ここで、閾値ＴＨ２は図６（Ｃ）に示したように適合点よりいずれの範囲で、基本トレース点火時期推定基本値ＢＢＥＳＴＴＲＡＣＥ２を補正するかを定める値（一定値）で、予め適合しておく。
【００８９】
従って、リフト・中心角誤差ＬＩＴＩＭＥＲＲが閾値ＴＨ２より大きいときには補正範囲にないので、ステップ１２４５に進んでトレース点火時期誤差補正量ＴＲＡＣＥＥＲＲ２＝０とする。
【００９０】
一方、リフト・中心角誤差ＬＩＴＩＭＥＲＲが閾値ＴＨ２以下であるときには適合点を中心とする所定の補正範囲にあると判断してステップ１２４６に進み、エンジン回転速度ＮＲＰＭと充填効率ＩＴＡＣから図３０を内容とするマップを参照することにより、トレース点火時期推定誤差補正量基本値ＢＴＲＡＣＥＥＲＲ２を算出する。このトレース点火時期推定誤差補正量基本値ＢＴＲＡＣＥＥＲＲ２は、上記（補４）式の値を求めるものである。
【００９１】
ステップ１２４７では、次式によりトレース点火時期誤差補正量ＴＲＡＣＥＥＲＲ２を算出する。
【００９２】
ＴＲＡＣＥＥＲＲ２＝ＢＴＲＡＣＥＥＲＲ２×（ＬＩＴＩＭＥＲＲ／ＴＨ２）
…（１５）
（１５）式右辺の（ＬＩＥＶＴＩＭＥＲＲ／ＴＨ２）は１以下の正の値である。吸気バルブ１５の実リフトＲＥＬＩＦＴが目標リフトＴＧＬＩＦＴと一致し、吸気バルブ１５の実中心角ＲＥＴＩＭＩＮＧも目標中心角ＴＧＴＩＭＩＮＧと一致するときリフト作動角・中心角誤差ＬＩＴＩＭＥＲＲ＝０となり、このとき、（１５）式よりトレース点火時期推定誤差補正量基本値ＢＴＲＡＣＥＥＲＲ２がそのままトレース点火時期誤差補正量ＴＲＡＣＥＥＲＲ２となる。また、吸気バルブ１５の実リフトＲＥＬＩＦＴが目標リフトＴＧＬＩＦＴから離れかつ吸気バルブ１５の実中心角ＲＥＴＩＭＩＮＧが目標中心角ＴＧＴＩＭＩＮＧから離れるほどリフト作動角・中心角誤差ＬＩＴＩＭＥＲＲが大きくなって閾値ＴＨ１に近づいて行き、リフト作動角・中心角誤差ＬＩＴＩＭＥＲＲが閾値ＴＨ１と一致したとき、トレース点火時期誤差補正量ＴＲＡＣＥＥＲＲ２はゼロとなる。すなわち、（１５）式は図６（Ｃ）に示したように適合点を中心とする所定の補正範囲で一点鎖線の特性から破線の特性へと進角側に移すための補正量を求める操作である。
【００９３】
上記のトレース点火時期推定誤差補正量基本値ＢＴＲＡＣＥＥＲＲ２の特性は、吸気バルブ１５の実リフトＲＥＬＩＦＴが適合点の吸気バルブ１５のリフト（つまり目標リフトＴＧＬＩＦＴ）に近づいてゆくほどステップ１２５で算出する基本トレース点火時期推定値ＢＥＳＴＴＲＡＣＥ２が、適合点の実トレースノック点火時期に近づいて行くように、つまり図６（Ｃ）に示した破線の特性が得られるように予め適合しておく。
【００９４】
ここで、第２実施形態の作用効果を説明する。
【００９５】
第２実施形態（請求項１、８に記載の発明）によれば、ＶＴＣ機構２７とＶＥＬ機構２８とを備えるエンジンにおいて、エンジンの回転速度ＮＲＰＭ、エンジンの負荷、ＶＴＣ角度（ＶＴＣ機構２７に与える指令値）及びＶＥＬ角度（ＶＥＬ機構２８に与える指令値）に基づいて基本トレース点火時期推定基本値ＢＢＥＳＴＴＲＡＣＥ２（トレース点火時期基本値）を算出し（図２１のステップ１２３、図２２参照）、この算出された基本トレース点火時期推定基本値ＢＢＥＳＴＴＲＡＣＥ１（トレース点火時期基本値）で火花点火を行うので、ＶＴＣ角度（ＶＴＣ機構２７に与える指令値）は変わらないのに、ＶＥＬ角度（ＶＥＬ機構２８に与える指令値）によって吸気バルブ１５の開時期や閉時期（開閉タイミング）が進角したり遅角したりした場合でも、基本トレース点火時期推定基本値ＢＢＥＳＴＴＲＡＣＥ２（トレース点火時期基本値）を精度良く算出することができ、基本トレース点火時期推定基本値ＢＢＥＳＴＴＲＡＣＥ２（トレース点火時期基本値）の算出（推定）精度が向上する。
【００９６】
第２実施形態（請求項２、９に記載の発明）によれば、基本トレース点火時期推定基本値ＢＢＥＳＴＴＲＡＣＥ２（トレース点火時期基本値）を、吸気バルブ１５のリフトＲＥＬＩＦＴの一次式で近似するので（図２２のステップ１２３７参照）、吸気バルブ１５のリフト、ＶＥＬ角度（ＶＥＬ機構２８に与える指令値）が相違する場合でも、基本トレース点火時期推定基本値ＢＢＥＳＴＴＲＡＣＥ２（トレース点火時期基本値）を容易に算出できる。
【００９７】
第２実施形態（請求項３、１０に記載の発明）によれば、基本トレース点火時期推定基本値ＢＢＥＳＴＴＲＡＣＥ２（トレース点火時期基本値）を吸気バルブ１５のリフトＲＥＬＩＦＴ（バルブリフト最大量）の一次式で近似する場合に基本トレース点火時期推定基本値のリフト感度ＬＩＳＬ（一次式の傾き）と基本トレース点火時期推定基本値のリフト切片ＬＩＩＮ（一次式の切片）とを吸気バルブ１５の中心角ＲＥＴＩＭＩＮＧの一次式で近似するので（図２２のステップ１２３３、１２３６参照）、吸気バルブ１５の中心角ＲＥＴＩＭＩＮＧ、ＶＴＣ角度（ＶＴＣ機構２７に与える指令値）が相違する場合でも、基本トレース点火時期推定基本値のリフト感度ＬＩＳＬ（一次式の傾き）と基本トレース点火時期推定基本値のリフト切片ＬＩＩＮ（一次式の切片）とを容易に算出できる。
【００９８】
基本トレース点火時期推定基本値のリフト感度ＬＩＳＬ（一次式の傾き）と基本トレース点火時期推定基本値のリフト切片ＬＩＩＮ（一次式の切片）とを、基本トレース点火時期推定基本値ＢＢＥＳＴＴＲＡＣＥ２（一次式で近似した第１のトレース点火時期基本値）が所定の適合点で実トレース点火時期と一致するように設定したとき、基本トレース点火時期推定基本値ＢＢＥＳＴＴＲＡＣＥ２（第１のトレース点火時期基本値）が適合点以外では実トレース点火時期より進角してノックが発生することがあり得るのであるが、第２実施形態（請求項５、１２に記載の発明）によれば、基本トレース点火時期推定基本値のリフト感度ＬＩＳＬ（一次式の傾き）と基本トレース点火時期推定基本値のリフト切片ＬＩＩＮ（一次式の切片）とを、基本トレース点火時期推定基本値ＢＢＥＳＴＴＲＡＣＥ２（第１のトレース点火時期基本値）がさらに適合点以外で実トレース点火時期より進角しないよう基本トレース点火時期推定基本値ＢＢＥＳＴＴＲＡＣＥ２（第１のトレース点火時期基本値）を遅角側にオフセットして設定するので、適合点以外でのノックの発生を回避できる。
【００９９】
第２実施形態（請求項６、１３に記載の発明）によれば、基本トレース点火時期推定基本値ＢＢＥＳＴＴＲＡＣＥ２（適合点以外で実トレース点火時期より進角しないよう遅角側にオフセットして設定している第２のトレース点火時期基本値）とトレース点火時期誤差補正量ＴＲＡＣＥＥＲＲ２（トレース点火時期補正量）とから基本トレース点火時期推定値ＢＥＳＴＴＲＡＣＥ２（トレース点火時期推定値）を算出し（図２１のステップ１２５参照）、この基本トレース点火時期推定値ＢＥＳＴＴＲＡＣＥ２（トレース点火時期推定値）を基本トレース点火時期推定基本値ＢＢＥＳＴＴＲＡＣＥ２（第２のトレース点火時期基本値）に代えて火花点火を行うと共に、適合点で基本トレース点火時期推定値ＢＥＳＴＴＲＡＣＥ２（トレース点火時期推定値）と実トレース点火時期とが一致するようにトレース点火時期推定誤差補正量基本値ＢＴＲＡＣＥＥＲＲ２（トレース点火時期補正量）を設定するので（図２３のステップ１２４６参照）、適合点で基本トレース点火時期推定値ＢＥＳＴＴＲＡＣＥ２（トレース点火時期推定値）と実トレース点火時期との間に点火時期誤差が発生してしまうことを避けることができる。
【０１００】
適合点でのみ基本トレース点火時期推定値ＢＥＳＴＴＲＡＣＥ２（トレース点火時期推定値）と実トレース点火時期とが一致するように基本トレース点火時期推定値ＢＥＳＴＴＲＡＣＥ２（トレース点火時期補正量）を設定したのでは、適合点とそれ以外とでトレース点火時期推定値と実トレース点火時期との間に点火時期差が生じトルク段差を誘発する可能性があるのであるが、第２実施形態（請求項７、１４に記載の発明）によれば、適合点を中心として所定の補正範囲を設け、この補正範囲において基本トレース点火時期推定値ＢＥＳＴＴＲＡＣＥ２（トレース点火時期推定値）が、基本トレース点火時期推定基本値ＢＢＥＳＴＴＲＡＣＥ２（第２のトレース点火時期基本値の直線）へと滑らかにつながるように、トレース点火時期誤差補正量ＴＲＡＣＥＥＲＲ２（適合点以外でのトレース点時期補正量）を算出するので（図２３のステップ１２４７参照）、こうしたトルク段差を誘発することがない。
【０１０１】
請求項１に記載の発明において、トレース点火時期基本値算出処理手順は図８のステップ０２３、図９によりまたは図２１のステップ１２３、図２２により果たされている。
【０１０２】
請求項８に記載の発明において、トレース点火時期基本値算出手段の機能は図８のステップ０２３、図９によりまたは図２１のステップ１２３、図２２により果たされている。
【図面の簡単な説明】
【０１０３】
【図１】本発明の第１実施形態のエンジンの点火時期制御装置の概略構成図。
【図２】ＶＥＬ機構及びＶＴＣ機構を備える場合のバルブリフト特性図。
【図３】ＶＴＣ角度、ＶＥＬ角度とトレース点火時期基本値の関係を示す特性図。
【図４】ＶＥＬ角度に対する近似式傾きの特性図。
【図５】ＶＥＬ角度に対する近似式切片の特性図。
【図６（Ａ）】トレース点火時期基本値の直線近似方法が妥当か否かを検討するために用いたＶＴＣ角度、ＶＥＬ角度とトレース点火時期の関係を示す特性図。
【図６（Ｂ）】トレース点火時期基本値の直線近似方法が妥当か否かを検討するために用いたＶＴＣ角度、ＶＥＬ角度とトレース点火時期の関係を示す特性図。
【図６（Ｃ）】トレース点火時期基本値の直線近似方法が妥当か否かを検討するために用いたＶＴＣ角度、ＶＥＬ角度とトレース点火時期の関係を示す特性図。
【図７】第１実施形態のトレース点火時期推定値の算出を説明するためのフローチャート。
【図８】第１実施形態の基本トレース点火時期推定値の算出を説明するためのフローチャート。
【図９】第１実施形態の基本トレース点火時期推定基本値の算出を説明するためのフローチャート。
【図１０】第１実施形態のトレース点火時期誤差補正量の算出を説明するためのフローチャート。
【図１１】第１実施形態の吸気バルブ作動角の特性図。
【図１２】第１実施形態の作動角感度中心角感度の特性図。
【図１３】第１実施形態の作動角感度中心角切片の特性図。
【図１４】第１実施形態の作動角切片中心角感度の特性図。
【図１５】第１実施形態の作動角切片中心角切片の特性図。
【図１６】第１実施形態の目標作動角の特性図。
【図１７】第１、第２の実施形態の目標中心角の特性図。
【図１８】第１実施形態の推定誤差補正量基本値の特性図。
【図１９】吸気バルブのバルブリフト特性図。
【図２０】第２実施形態のトレース点火時期推定値の算出を説明するためのフローチャート。
【図２１】第２実施形態の基本トレース点火時期推定値の算出を説明するためのフローチャート。
【図２２】第２実施形態の基本トレース点火時期推定基本値の算出を説明するためのフローチャート。
【図２３】第２実施形態のトレース点火時期誤差補正量の算出を説明するためのフローチャート。
【図２４】第２実施形態のリフトの特性図。
【図２５】第２実施形態のリフト感度中心角感度の特性図。
【図２６】第２実施形態のリフト感度中心角切片の特性図。
【図２７】第２実施形態のリフト切片中心角感度の特性図。
【図２８】第２実施形態のリフト切片中心角切片の特性図。
【図２９】第２実施形態の目標リフトの特性図。
【図３０】第２実施形態の推定誤差補正量基本値の特性図。
【符号の説明】
【０１０４】
１１  点火装置
２７  ＶＴＣ機構（可変バルブタイミング機構）
２８  ＶＥＬ機構（可変バルブ機構）
３１  エンジンコントローラ
３４  カム角度センサ
３７  ノックセンサ
４３  ＶＥＬ角度センサ
【図１】
2007332926.tif 000003
【図２】
2007332926.tif 000004
【図３】
2007332926.tif 000005
【図４】
2007332926.tif 000006
【図５】
2007332926.tif 000007
【図６（Ａ）】
2007332926.tif 000008
【図６（Ｂ）】
2007332926.tif 000009
【図６（Ｃ）】
2007332926.tif 000010
【図７】
2007332926.tif 000011
【図８】
2007332926.tif 000012
【図９】
2007332926.tif 000013
【図１０】
2007332926.tif 000014
【図１１】
2007332926.tif 000015
【図１２】
2007332926.tif 000016
【図１３】
2007332926.tif 000017
【図１４】
2007332926.tif 000018
【図１５】
2007332926.tif 000019
【図１６】
2007332926.tif 000020
【図１７】
2007332926.tif 000021
【図１８】
2007332926.tif 000022
【図１９】
2007332926.tif 000023
【図２０】
2007332926.tif 000024
【図２１】
2007332926.tif 000025
【図２２】
2007332926.tif 000026
【図２３】
2007332926.tif 000027
【図２４】
2007332926.tif 000028
【図２５】
2007332926.tif 000029
【図２６】
2007332926.tif 000030
【図２７】
2007332926.tif 000031
【図２８】
2007332926.tif 000032
【図２９】
2007332926.tif 000033
【図３０】
2007332926.tif 000034