(57)【要約】  ピストンエンジンにおける放熱率を自動的に算出する方法であって、クランク角に応じてシリンダーの圧力を測定するステップと（図を参照）、前記測定に基づき、圧縮ストロークについて第一のポリトロープ指数を算出するステップと（Ｍにおいて）、前記測定に基づき、膨張ストロークについて第二のポリトロープ指数を算出するステップと（Ｍ’において）、前記圧縮ストロークと前記膨張ストローク間のクランク角の間隔について、前記第一および第二ポリトロープ指数の補間を実行するステップと、および、補間された前記ポリトロープ指数に基づき正味の放熱の算出を実行するステップと、を備えたことを特徴とする方法。【選択図】  図１
2010518384.tif 000002
【特許請求の範囲】
【請求項１】
ピストンエンジンにおける放熱率を自動的に算出する方法であって、
a．クランク角に応じてシリンダーの圧力を測定するステップと、
b．前記測定に基づき、圧縮ストロークについて第一のポリトロープ指数(polytrope exponent)を算出するステップと、
c．前記測定に基づき、膨張ストロークについて第二のポリトロープ指数を算出するステップと、
d．前記圧縮ストロークと前記膨張ストローク間のクランク角の間隔について、前記第一および第二ポリトロープ指数の補間を実行するステップと、および
e．補間された前記ポリトロープ指数に基づき正味の放熱の算出を実行するステップと、を備えたこと、
を特徴とする方法。
【請求項２】
請求項１に基づく方法であって、さらに、圧力センサーのオフセットを算出するステップを備えること、
を特徴とする方法。
【請求項３】
請求項２の方法において、前記圧力センサーのオフセットは、前記ポリトロープ指数を、圧力と体積との間のポリトロープ関係に用いることにより推定すること、
を特徴とする方法。
【請求項４】
請求項１の方法において、前記補間は、線形補間であること、
を特徴とする方法。
【請求項５】
請求項１の方法において、前記補間は、三次スプライン補間(cubic spline interpolation)であること、
を特徴とする方法。
【請求項６】
請求項１の方法において、前記ステップ（b）および（c）は、
i.   前記ポリトロープ指数の初期値を設定するステップと、
ii. 前記ポリトロープ指数の初期値に基づいてモデル圧力の軌跡を算出するステップと、
iii.   各クランク角に関するモデルのミスマッチを算出するステップと、
iv.  二乗されたミスマッチの合計のポリトロープ指数に関するデリバティブ(derivative)を算出するステップと、
v.  前記ステップの結果に基づき、所定の量だけ前記ポリトロープ指数を増加または減少させるステップと、
vi. ポリトロープ指数の値を増加または減少させながらステップiiからivを繰り返すステップと、
vii. 二乗されたミスマッチの合計のポリトロープ指数に関して第二のデリバティブを算出するステップと、
viii.  負の符号を有する前記第二のデリバティブで除された前記デリバティブに対応する量により前記ポリトロープ指数の値を変更するステップと、および
ix.  負の符号を有する前記第二のデリバティブで除された前記デリバティブの絶対値に対応する量が、所定の値を下回るまでステップii、 iii、iv、viiおよびviiiを繰り返すステップと、を備えること、
を特徴とする方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ピストンエンジンにおける放熱率の算出方法に関する。
【背景技術】
【０００２】
ピストンエンジンの研究において、放熱率は、長い間、非常に重要な診断ツールであった。放熱率を測定するための標準的な手順としては、エンジンを駆動させ、クランク角の位置に応じてシリンダにおける圧力を測定し、いわゆる“駆動時(motored)”の圧力の軌跡(pressure trace)を得、その後、エンジンに点火し、“点火時(fired)”の圧力の軌跡を記録する。駆動時の圧力軌跡と点火時の圧力軌跡の間の差は、燃焼からの発散、すなわち、放熱のために、どの圧力が効率的であるかを判断するために用いることができる。
前述の方法は、非常に正確な結果を提供するが、残念なことに、良好な結果を得るためには、シリンダー壁への熱伝達、ガスの漏れ、隙間の量等のいくつかの異なる損失効果について考慮しなければならず、非常に多くの労働力を必要とする。
【０００３】
近年、搭載型の診断システムが、乗用車ならびにトラックメーカーにとって、エンジンのパラメーターの監視を可能とするため必要不可欠なツールとなっており、近い将来、放熱率を監視することにより、エンジンシリンダーの放熱率を監視することが求められるであろうし、可能な限り効率的で環境に優しい燃焼を提供する、点火タイミング、燃料噴射タイミングおよび燃料噴射量等のエンジン動作パラメーターを最適化することも可能となるであろう。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】なし
【非特許文献】
【０００５】
【非特許文献１】なし
【発明の概要】
【発明が解決しようとする課題】
【０００６】
上述のように、放熱の分析を実行するには非常に手間がかかっており、熱伝達の値、漏れおよび隙間の量は、エンジンが古くなったり摩耗したりすると変動するので、放熱分析結果を得るための既知の方法は、搭載型の診断システムにあまり適さないものになってしまう。
【０００７】
本発明の目的は、放熱の算出を実行する方法を提供することである。かかる方法は、“自己調節”すなわち、シリンダー壁への熱伝達の差、漏れおよび隙間の量を補償することにより、放熱を迅速に算出することを可能にするとともに、絶対圧を提供している圧力変換器および圧力変化を提供する圧力変換器の両方に用いることができるものである。
【課題を解決するための手段】
【０００８】
本発明は、以下のステップを備えた方法であって、
１．クランク角に応じてシリンダーの圧力を測定するステップと、
２．前記測定に基づき、圧縮ストロークについて第一のポリトロープ指数(polytrope exponent)を算出するステップと、
３．前記測定に基づき、膨張ストロークについて第二のポリトロープ指数を算出するステップと、
４．前記圧縮ストロークと前記膨張ストローク間のクランク角の間隔について、前記第一および第二ポリトロープ指数の補間を実行するステップと、および
５．補間された前記ポリトロープ指数に基づき正味の放熱の算出を実行するステップと
を特徴とする方法、を提供することにより上述の課題を解決する。
【０００９】
シリンダーの圧力を監視するのに用いられる圧力センサーが、絶対圧スケールを用いないセンサーでなければならない場合、圧力センサーのオフセットを算出するステップをさらに含ませることも可能である。かかる圧力オフセットは、前記ポリトロープ指数を、圧力と体積との間のポリトロープ関係に用いることにより推定することも可能である。
【００１０】
できるだけ簡単で効率的な補間を得るため、かかる補間は、線形補間であってもよい。他の実施形態において、かかる補間は、３次スプライン補間法(cubic spline interpolation)であってもよい。
【００１１】
前記ステップ（２）および（３）が、以下のステップ：
i.   前記ポリトロープ指数の初期値を設定するステップと、
ii. 前記ポリトロープ指数の初期値に基づいてモデル圧力の軌跡を算出するステップと、
iii.   各クランク角に関するモデルのミスマッチを算出するステップと、
iv.  二乗されたミスマッチの合計のポリトロープ指数に関するデリバティブ(derivative)を算出するステップと、
v.  前記ステップの結果に基づき、所定の量だけ前記ポリトロープ指数を増加または減少させるステップと、
vi. ポリトロープ指数の値を増加または減少させながらステップiiからivを繰り返すステップと、
vii. 二乗されたミスマッチの合計のポリトロープ指数に関して第二のデリバティブを算出するステップと、
viii.  負の符号を有する前記第二のデリバティブで除された前記デリバティブに対応する量により前記ポリトロープ指数の値を変更するステップと、および
ix.  前記第二のデリバティブで除された前記デリバティブの絶対値に対応する量が、所定の値を下回るまでステップii、 iii、iv、viiおよびviiiを繰り返すステップと、を備えれば有益である。
【図面の簡単な説明】
【００１２】
以下において、本発明は、唯一の図面を参照しつつ説明される。  
【図１】図１は、ＨＣＣＩエンジンの圧縮／膨張ストロークの圧力軌跡を表すグラフである。
【発明を実施するための形態】
【００１３】
図１は、ピストンエンジンにおけるクランク角に応じたシリンダーの圧力を示している、圧力の軌跡ＰＴを表すグラフを示している。かかるグラフは、１回のエンジン周期、すなわち、−１８０度のクランク角度であるＣＡＤ（下死点すなわちＢＤＣと等しい）から、０度のクランク角度（上死点すなわちＴＤＣと等しい）まで、さらに、ＢＤＣと等しい１８０度のクランク角度に達するまでのピストン位置におけるシリンダーの圧力を示している。
既述のように、圧力の軌跡は、１周期におけるシリンダーの圧力を示しており、かかる周期は、シリンダー内の充てん体の燃焼を伴うものである（４サイクルエンジンの他の周期は、気体の交換に費やされる）。かかる軌跡においていくつかの重要な事象を見ることができる。
【００１４】
クランク角度が−１８０度から０度に達するまで充てん体が圧縮され、シリンダー内の体積が減少することにより圧力が上昇する。理想的なケースにおいては、かかる圧力上昇は、理想気体の法則に沿うが、シリンダー壁に対する、あるいは、そこからの熱伝達、ならびに、シリンダー容積からの気体の漏れは、前記軌跡を相殺してしまう。  
【００１５】
かかる圧力の軌跡ＰＴは、ＨＣＣＩエンジンと呼ばれる、すなわち、充てん体が自動的に点火するまで（ディーゼルエンジンのように）圧縮される、燃料と空気の混合気を含むエンジン（スパーク点火エンジンのような）に記録される。ＨＣＣＩエンジンの背景にある法則は、多数の技術論文（例えば、大西らによる、“アクティブ  サーモーアトモスフェアコンビュージョン（ＡＴＡＣ）−−内燃エンジンの新しい燃焼プロセス”、ＳＡＥ論文７９０５０１、NajtおよびFosterによる、“圧縮−点火された充てん体の均一な燃焼”、ＳＡＥ論文８３０２６４、Thringによる“（ＨＣＣＩ）エンジンにおける均一な−充てん体の圧縮−点火”、ＳＡＥ論文８９２０６８、および、Olsson等による、“ターボチャージ式二元燃料ＨＣＣＩエンジン”ＳＡＥ論文２００１−０１−１８９６参照）ならびに、米国特許番号６、２８６、４８２等の特許に説明されている。
【００１６】
図１において見られるように、圧力の軌跡ＰＴは、ＴＤＣ、つまりクランク角度が０度の後に急上昇する。かかる圧力上昇は、燃焼が原因であり、燃焼がなされない場合、ＴＤＣ後に圧力はその代わりに下降する。圧力がピークに達した後、圧力は低下する、この低下は、そこに気体が含まれる、体積が増加する方向へのピストンの移動から生じる体積の増加が原因である。
【００１７】
前の文で述べたとおり、放熱分析を実行するための最も一般的（かつ正確）な方法は、“点火時”と“駆動時”の圧力軌道、すなわち、全く燃焼していないエンジンおよび燃焼により駆動しているエンジンにおいて得られた圧力軌跡、を比較する方法である。他の全ての基準が等しいと仮定すると、点火時の圧力軌跡と駆動時の圧力軌跡間の圧力軌跡における全ての相違は、燃焼が原因であることが推定できる。確立した物理法則を適用することで、燃焼のレートから控除されるもの、すなわち、放熱を推定することが可能となる。
【００１８】
本発明は、比較のために、エンジンを駆動させ、“駆動時”の圧力軌跡を得る必要なしに放熱の算出を可能にする方法に関する。以下の基本的な推定がなされる：
燃焼によって放出される熱エネルギーは、熱として扱われ、環境に関係する唯一の作業は、それがクランク角度０からクランク角度１８０まで（すなわち、ＴＤＣからＢＤＣまで）移動した場合にピストンを通じて実行される作業である。このことは、以下の式により表される
【００１９】
【数１ａ】
2010518384.tif 000003
ここで、
ｄＵは、内部エネルギーの変化であり、
ｄＱは、熱伝達であり、ならびに、
ｄＷは、燃焼室内の気体によって実行される力学的仕事である。
【００２０】
内部エネルギーの変化、すなわち、ｄＵを理想気体の式に代入することにより、以下のように決定される
【００２１】
【数１ｂ】
2010518384.tif 000004
ここで、
nは、シリンダー内部に捕らえられたモル数であり、
Cvは、一定容積における熱容量であり、
ｄＴは、温度変化であり、
ｐは、圧力であり、ならびに
ｄｖは、体積の変化である。
【００２２】
上記の式にｄＴを代入することにより、圧力並びに体積について、ｄＴ、すなわち、温度変化を表すため理想気体の法則を用いることができる。
【００２３】
【数１】
2010518384.tif 000005
ここで、Ｒは、一般気体定数である。
【００２４】
理解することができるように、燃焼がなされず、あるいは、壁とシリンダー充てん体間に熱の伝達がない場合、ｄＱはゼロと推定可能であり、従来の発熱計算においては、熱伝達の値を推定しなければならならず、さもなければ、その発熱計算が不正確となっていた。
【００２５】
しかしながら、本発明によれば、熱伝達を推定することなく、それに代え、Ｃｖを補正することにより熱伝達が補正される。
【００２６】
本発明に基づく方法は、燃焼の前後にＣｖの“補正”を伴う（充てん体の温度が高いので、無論、燃焼後のほうが熱損失は大きい）。燃焼中、Ｃｖの値は補間される。以下の式（２から３）は、ｄＱ＝０の仮定の下でのみ有効である。簡単にするため、式（１）は、以下のように書き換えられる：
【００２７】
【数２】
2010518384.tif 000006
ここで、γ＝Ｃｐ／Ｃｖ(Ｃｐは、一定圧力における熱容量である)。
【００２８】
分母を削除し(Dropping the denominators）、式（２）における２つの差を組み合わせることにより、ポリトロープ関係が提供される
【００２９】
【数３】
2010518384.tif 000007
したがって、ｄＱが０の場合（燃焼がなく熱伝達がないのと同じ）、ｐＶγは一定である。
【００３０】
いくつかの場合、シリンダー内圧力の測定を実行するため未知の圧力レベルを有する圧力変換器が用いられる；このことは、圧力レベルが特定されなければならないことを意味している。この目的のため、測定された圧力（Ｐｍ）を、以下のように表すことが可能である
【００３１】
【数４】
2010518384.tif 000008
ここで、ｐは実際の圧力であり、Δｐは相殺される圧力である。
【００３２】
式（３）の統合ならびに式（４）の適用は、相殺圧力を特定するために用いられる式を生み出す。（既知の圧力レベルを有する圧力変換器が用いられた場合、このステップは不要であることにご注意ください）。以下の式においては、γの代わりにκが用いられることにご注意ください；γは、ＣｐおよびＣｖ間の実際の比率の認められたシンボルであり、以下の式は、ＣｐおよびＣｖ間の実際の比率を提供することはない。むしろ、以下の式は、実際に測定されたシリンダー圧と一致する曲線のポリトロープ指数を提供する。次に、相殺圧力を特定するために用いられる式は、
【００３３】
【数５】
2010518384.tif 000009
となる。
【００３４】
ここで、
κは、ポリトロープ指数であり、
Ｐｍ０は、基準クランク位置における測定圧力であり、および
Ｖ０ は、前記基準クランク位置の体積である。
【００３５】
圧力測定ノイズの影響を減少させるには、前記基準クランク位置をできるだけＴＤＣに近づくよう選択すれば有益であろう；しかし、前記基準クランク位置が燃焼が行われないクランク位置の上で選択されることが重要である。
【００３６】
物理的な寸法を取り除くため、式（５）は、以下のように書き換えられる
【００３７】
【数６】
2010518384.tif 000010
式（６）は、以下のように表される
【００３８】
【数７】
2010518384.tif 000011
ここで、
【００３９】
【数８】
2010518384.tif 000012
である。
【００４０】
式（７）は、ニュートンの非線型最小二乗法のアルゴリズム法（２００１年にルンド大学工学部（Lund University, Faculty of engineering）でPer  Tunestal氏の博士論文として発表された“圧力センサーを用いることによる内燃エンジンのシリンダー内の空気／燃料の割合の推定”）を可能にする形式のものであり、相違点は、前記κが出力、すなわちｙに依存することだけである。
【００４１】
この依存に適合するため、ここでは、かかるアルゴリズムがどのように補正されるかについて示しており、クランク角のベクトルにより決定された出力ならびに以下で与えられるリグレッサー(regressors)を伴う、同じ表記法が採用される：
【００４２】
【数９】
2010518384.tif 000013
ここで、ｙ１・・・ｙｎは、ｙ（クランク角位置の等間隔の配列については（式（８）を参照）の算出された値であり、
φ１・・・φｎは、φ（クランク角位置の同じ等間隔の配列については（式（８）を参照）の算出された値である。
【００４３】
これは、以下のベクトル式を提供する
【００４４】
【数１０】
2010518384.tif 000014
残りのベクトルは、以下により定義される
【００４５】
【数１０ａ】
2010518384.tif 000015
また、損失関数は、以下のように定義される
【００４６】
【数１０ｂ】
2010518384.tif 000016
ここで、ＤＴは、Ｄの転置行列であり、したがって、Ｊはスカラーである。  ２００１年のルンド大学工学部のPer  Tunestal氏の博士論文“圧力センサーを用いることによる内燃エンジンのシリンダー内の空気／燃料の割合の推定”においては、以下のことを示している
【００４７】
【数１１】
2010518384.tif 000017
ここで、
Ｄ＝ＰＹであり、
Ｐは、以下により説明される左右対称の射影行列（symmetric projection matrix）である
【００４８】
【数１１ａ】
2010518384.tif 000018
当業者によく知られているように、有限差分法は、κについてＪから派生するもの(delivertive)に関する分析的な表現を必要とする。
【００４９】
これは、式（１０）からκを識別し、ＹとＰの両方がκの関数であることに留意しつつ、Ｐが左右対称の行列であるという事実を用いることにより得ることができる。
【００５０】
【数１１ｂ】
2010518384.tif 000019
第一項は、以下のように２００１年のルンド大学工学部のPer  Tunestal氏の博士論文“圧力センサーを用いることによる内燃エンジンのシリンダー内の空気／燃料の割合の推定”と同じである
【００５１】
【数１１ｃ】
2010518384.tif 000020
式（１１）中の残りに関する式を用いて第二項を書き換えることができる。
【００５２】
二乗されたミスマッチＪの合計のκに関するデリバティブは、以下のように書き換えることができる
【００５３】
【数１２】
2010518384.tif 000021
ここで、ｄＹ／ｄκおよびｄθ／ｄκは、以下のように表される：
【００５４】
【数１２ａ】
2010518384.tif 000022
ｄＪ／ｄκの分析的な表現を用いることにより、損失関数を最小化するために有限差分ニュートン法を適用することができる。かかる方法は、超一次収束(superlinear convergence)を特徴とし；実際に、このことは、収束中にいくつかの繰り返しが存することを意味する。繰り返しκは、以下のように定義される
【００５５】
【数１２ｂ】
2010518384.tif 000023
ここで、Ｈκは、繰り返しκにおけるｄ２Ｊ／ｄκ２の有限差分近似（finite difference approximation）である。
【００５６】
理解可能であるように、クランク位置間における前記κの値は、絶対に燃焼が起こらないよう選ばれるように算出され；ＳＩエンジン（すなわち、スパーク点火エンジン）において、燃焼の開始は点火タイミングより早くなることはなく、ＣＩエンジン（圧縮点火またはディーゼルエンジン）について、燃焼の開始は燃料の最初の噴射後になる。
【００５７】
図１において、燃焼前ならびに後にκの値を算出するため、圧力の軌跡ＭおよびＭ’の厚い部分に対応するクランク位置がそれぞれ用いられ；見て判るように、これらのクランク位置は、燃焼タイミングから“安全”な距離上に位置する。しかし、κの値が、燃焼中近傍のκの値と大きく異なる（理解できるように、対象となる(investigated)クランクの位置が燃焼に近く付けば近づくほど、圧縮および膨張中の熱伝達の影響は、燃焼中の熱伝達の影響と近似する）ので、燃焼事象から離れすぎた位置に存するクランク間隔を選択しないことが重要である。
【００５８】
上述の基準クランク位置は、クランク位置ＭならびにＭ’内に位置すべきであることに注意することができる。位置ＭならびにＭ’が、圧力軌跡ＰＴの“曲がった”箇所に位置することが好ましいことにも注意すべきであり；当業者によく知られているように、前記軌跡が曲線を示す場合、指数フィッテング動作(exponent fitting operation)が最も正確である。
【００５９】
上記の式から理解可能なように、本発明に基づく方法は、κについていくつかの入力が必要であり、入力がない場合、かかる方法は機能しない。本発明の好ましい実施形態において、実験は、κの値が約１．２から１．４の範囲であることを示している（きれいな空気、γ、すなわち、熱力学的に正しい等エントロピー指数が１．４に近い）ので、κの初期値は１．３に設定される。合理的にうまく選択されたκの初期値を用いると、かかる方法は機能するが、値を下手に選択すると、正しい結果を得るため必要となる繰り返しの数が大きくなる。
【００６０】
第一および第二の繰り返しは、連続的な繰り返しと若干異なる。上述のように、エラーの第一のデリバティブ（すなわち、式１０のＪ）は、分析的に得ることができる（式１２）が、第二のデリバティブについては、２の連続的な周期の結果を比較しなければならない。上記最後の式、すなわち、κκ＋１＝κκ＋δκ、で定義されるように、以下の繰り返しに関するκの新しい値は、κに対するＪの第二のデリバティブにより除されたκについてのエラーのマイナスのデリバティブである；しかし、第一の繰り返しの前には、第二のデリバティブについての値が存在しない。この問題を解決するため、第一の繰り返しの間（エラーならびにかかるエラーのデリバティブが算出される場合）、κに対するＪの第一のデリバティブの符号により決定された方向において、所定の量だけκの値が変更される。
【００６１】
このことは、以下の繰り返しの基準を提供する。
【００６２】
第一の繰り返しについて、κについての所定値が設定される。本発明の好ましい実施形態において、この値は、１．３に設定されるが、例えば１から２の範囲の、いずれの値に設定するようにしてもよい。第一の繰り返しの間、エラーの第一のデリバティブ（すなわちＪ）が、算出可能であり、第一のデリバティブの符号に基づき、κの値は、所定の値、例えば、０．０１だけ増減される。
【００６３】
第二のならびに以下の繰り返しについては、第二のデリバティブにより除された第一のデリバティブのマイナスによって前記κの値を増減少させることが可能である。
【００６４】
最後に、ＭおよびＭ’厚い部分に対応する２つのクランク位置間のκの値の補間について、述べておかなければならないことがある。この補間は、もちろん、線形補間であってもよいが、いずれのタイプの補間であってもよく；よりよい結果をもたらす補間法は、おそらく“三次スプライン補間（cubic spline interpolation）”である。
【００６５】
放熱は、κにより、周知の放熱式に等エントロピー指数(isentropic exponent)γを代入することにより算出される。
【００６６】
これにより、以下の修正された放熱式が提供される：
【００６７】
【数１２ｃ】
2010518384.tif 000024
上では、本発明のいくつかの実施形態が説明されている。しかし、この説明は、本発明の範囲を限定することを意図したものではなく、それに代え、本発明の範囲は、添付された特許請求の範囲において限定されることを意図するものである。
【図１】
2010518384.tif 000025
【国際調査報告】
2010518384.tif 000026
2010518384.tif 000027
2010518384.tif 000028