(57)【要約】【課題】  ＣＶＴの実際のベルト摩擦特性を考慮して、早期にベルト滑りを予測する。【解決手段】  μ勾配算出回路３２は、エンジントルク、プライマリ回転角速度ω1、セカンダリ回転角速度ω2、ベルト掛かり径Ｒ1，Ｒ2、プライマリ回転角加速度、セカンダリ回転角加速度、プライマリベルト挟圧力Ｆ1、セカンダリベルト挟圧力Ｆ2に基づいて、ベルトμ勾配ｋ1、ｋ2を求める。ベルト滑り予測回路３３は、ベルトμ勾配ｋ1の絶対値、ベルトμ勾配ｋ2の絶対値のそれぞれと、しきい値生成回路２９で生成されたしきい値とを比較して、いずれか一方のベルトμ勾配の絶対値がしきい値以下となったときに、ベルト滑りを予測する。【選択図】  図３
2006002795.tif 000002
【特許請求の範囲】
【請求項１】
入力側プーリ、出力側プーリ、前記入力側プーリ及び前記出力側プーリに掛けられたベルトを有する無段変速機のベルト滑りを予測するベルト滑り予測装置であって、
前記入力側プーリの入力回転速度を検出する入力回転速度検出手段と、
前記出力側プーリの出力回転速度を検出する出力回転速度検出手段と、
前記無段変速機への入力トルクを推定する入力トルク推定手段と、
前記入力側プーリ及び前記出力側プーリのいずれかである第１のプーリのベルト挟圧力を検出する第１のベルト挟圧力検出手段と、
前記入力回転速度検出手段で検出された入力回転速度と、前記出力回転速度検出手段で検出された前記出力回転速度と、前記入力トルク推定手段で推定された入力トルクと、前記第１のベルト挟圧力検出手段で検出されたベルト挟圧力と、に基づいて、前記第１のプーリにおけるベルトμ勾配を推定するベルトμ勾配推定手段と、
前記ベルトμ勾配推定手段で推定されたベルトμ勾配に基づいて、前記ベルト滑りを予測する予測手段と、
を備えたベルト滑り予測装置。
【請求項２】
ベルトμ勾配推定手段は、ベルト摩擦特性における動作点について少なくともベルトμ勾配を用いて表したときのプーリの回転慣性に関する関係式を用いて、前記ベルトμ勾配を推定する
請求項１に記載のベルト滑り予測装置。
【請求項３】
前記予測手段は、前記ベルトμ勾配推定手段で推定されたベルトμ勾配の絶対値がしきい値を超えたときに、ベルト滑りを予測する
請求項１または請求項２に記載のベルト滑り予測装置。
【請求項４】
入力側プーリ及び前記出力側プーリのうち、前記第１のプーリと異なる第２のプーリのベルト挟圧力を検出する第２のベルト挟圧力検出手段をさらに備え、
前記ベルトμ勾配推定手段は、さらに、前記第２のベルト挟圧力検出手段で検出されたベルト挟圧力に基づいて、前記第２のプーリにおけるベルトμ勾配を推定する
請求項１または請求項２に記載のベルト滑り予測装置。
【請求項５】
前記予測手段は、ベルトμ勾配推定手段で推定された第１及び第２のベルトμ勾配の絶対値の少なくとも一方がしきい値を超えたときに、ベルト滑りを予測する
請求項４に記載のベルト滑り予測装置。
【請求項６】
前記ベルト摩擦特性は、線形領域及び非線形領域の両方を含むことを特徴とする
請求項２から請求項５のいずれか１項に記載のベルト滑り予測装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ベルト滑り予測装置に係り、特に無段変速機（ＣＶＴ）のベルト滑りを予測するベルト滑り予測装置に関する。
【背景技術】
【０００２】
従来、ベルト−プーリ間のベルト滑り速度に対する摩擦特性μが飽和する点に着目し、滑り速度に対するベルトμの接線勾配（ベルトμ勾配）を検出し、そのμ勾配が０近傍となった時にベルト滑りを検出する技術が提案されている（特許文献１を参照。）。
【０００３】
特許文献１には、プライマリ、セカンダリにおけるベルト−プーリ間の摩擦トルクの釣り合いに関する運動方程式と、ベルトの摩擦力に関する運動方程式とに基づき、プライマリ、セカンダリの回転速度信号を用いてベルト滑り速度に対するμ勾配を検出し、μ勾配の低下からベルト滑りを検出することが記載されている。
【特許文献１】特開２００３−２１５７８号公報
【発明の開示】
【発明が解決しようとする課題】
【０００４】
特許文献１では、サンプリング期間中にトルク変化、掛かり径変化、挟圧力変化がないという仮定の下でμ勾配算出式を導出していた。しかし、実際には、これらの変化を無視することができない問題があった。
【０００５】
特に、トルクはＣＶＴを駆動する入力値であり、トルク変化時にμ勾配の検出精度が低下してしまう問題があった。また、特許文献１ではベルト運動を接線運動として記述されているが、実際には回転体運動である。したがって、運動方程式に予め誤差が含まれているので、μ勾配の検出精度が高くない問題もあった。
【０００６】
本発明は、上述した課題を解決するために提案されたものであり、ＣＶＴの実際のベルト摩擦特性を考慮してμ勾配の推定精度を向上させるとともに、早期にベルト滑りを予測するベルト滑り予測装置を提供することを目的とする。
【課題を解決するための手段】
【０００７】
本発明に係るベルト滑り予測装置は、入力側プーリ、出力側プーリ、前記入力側プーリ及び前記出力側プーリに掛けられたベルトを有する無段変速機のベルト滑りを予測するベルト滑り予測装置であって、前記入力側プーリの入力回転速度を検出する入力回転速度検出手段と、前記出力側プーリの出力回転速度を検出する出力回転速度検出手段と、前記無段変速機への入力トルクを推定する入力トルク推定手段と、前記入力側プーリ及び前記出力側プーリのいずれかである第１のプーリのベルト挟圧力を検出する第１のベルト挟圧力検出手段と、前記入力回転速度検出手段で検出された入力回転速度と、前記出力回転速度検出手段で検出された前記出力回転速度と、前記入力トルク推定手段で推定された入力トルクと、前記第１のベルト挟圧力検出手段で検出されたベルト挟圧力と、に基づいて、前記第１のプーリにおけるベルトμ勾配を推定するベルトμ勾配推定手段と、前記ベルトμ勾配推定手段で推定されたベルトμ勾配に基づいて、前記ベルト滑りを予測する予測手段と、を備えている。
【０００８】
入力回転速度検出手段は、入力側プーリの入力回転速度に限らず、これと等価な物理量、例えばエンジン回転数、入力回転数を検出してもよい。同様に、出力回転速度検出手段は、出力側プーリの出力回転速度に限らず、これと等価な物理量、例えば出力回転数を検出してもよい。
【０００９】
入力トルク推定手段は、無段変速機の入力側プーリに入力される入力トルクを推定する。入力トルク推定手段は、入力トルクとしてエンジントルクを推定してもよい。つまり、入力トルク推定手段の代わりにエンジントルク推定手段を用いてもよい。エンジントルクは例えばエンジンへの流入空気量や車速に基づき推定可能であるが、推定手法は特に限定されない。
【００１０】
第１のベルト挟圧力検出手段は、前記入力側プーリ、前記出力側プーリのいずれか一方である第１のプーリのベルト挟圧力を検出する。なお、第１のベルト挟圧力検出手段は、ベルト挟圧力に限らず、これに等価な物理量、例えばプーリを構成する１組のシーブのうち、一方のシーブを他方のシーブに押圧するためのシーブ圧を検出してもよい。
【００１１】
ベルトμ勾配推定手段は、入力回転速度、出力回転速度、入力トルク、前記第１のベルト挟圧力検出手段で検出されたベルト挟圧力に基づいて、第１のプーリにおけるベルトμ勾配を推定する。ベルトμ勾配とは、ベルト滑り速度（又はベルト滑り率）に対するベルト摩擦係数を表したベルト摩擦特性上における動作点の接線勾配をいい、ベルトの滑りやすさを表している。例えば、ベルトμ勾配が大きい場合ベルトは滑りにくい状態を表し、ベルトμ勾配が小さい場合ベルトは滑りやすい状態を表す。
【００１２】
そこで、予測手段は、ベルトμ勾配推定手段で推定されたベルトμ勾配に基づいて、ベルト滑りを予測する。
【００１３】
したがって、本発明に係るベルト滑り予測装置は、入力回転速度、出力回転速度、入力トルク、ベルト挟圧力に基づいて、第１のプーリにおけるベルトμ勾配を推定し、ベルトμ勾配に基づいてベルト滑りを予測することにより、早期にベルト滑りを予測することができる。
【００１４】
ベルトμ勾配推定手段は、ベルト摩擦特性における動作点について少なくともベルトμ勾配を用いて表したときのプーリの回転慣性に関する関係式を用いて、前記ベルトμ勾配を推定してもよい。
【００１５】
ベルト摩擦特性の動作点は、ベルトμ勾配のほかにベルト摩擦切片を用いることによって、例えば、ベルト滑り速度（又はベルト滑り率）の１次関数で表される。ベルト滑り速度（又はベルト滑り率）が小さいときは、ベルト摩擦特性は線形性を有し、このときベルトμ勾配は大きな値となる。ベルト滑り速度（又はベルト滑り率）が大きくなると、ベルト摩擦特性は非線形性を有し、このときベルトμ勾配は小さな値となる。ベルトμ勾配がほぼ０なると、ベルト滑りの前兆であるマクロスリップの限界となる。
【００１６】
前記予測手段は、前記ベルトμ勾配推定手段で推定されたベルトμ勾配の絶対値がしきい値を超えたときに、ベルト滑りを予測してもよい。
【００１７】
また、前記ベルト滑り予測装置は、入力側プーリ及び前記出力側プーリのうち、前記第１のプーリと異なる第２のプーリのベルト挟圧力を検出する第２のベルト挟圧力検出手段をさらに備え、前記ベルトμ勾配推定手段は、さらに、前記第２のベルト挟圧力検出手段で検出されたベルト挟圧力に基づいて、前記第２のプーリにおけるベルトμ勾配を推定してもよい。
【００１８】
これにより、前記ベルト滑り予測装置は、入力側及び出力側プーリにおける第１及び第２のベルトμ勾配を推定することができるので、より高精度にベルト滑りを予測することができる。
【００１９】
このとき、予測手段は、ベルトμ勾配推定手段で推定された第１及び第２のベルトμ勾配の絶対値の少なくとも一方がしきい値を超えたときに、ベルト滑りを予測ればよい。
【００２０】
これにより、前記ベルト滑り予測装置は、入力側プーリとベルトの間、出力側プーリとベルトの間のいずれかで発生する可能性のあるベルト滑りを早期に予測することができる。
【００２１】
さらに、前記ベルト摩擦特性は、線形領域及び非線形領域の両方を含むものであってもよい。
【発明の効果】
【００２２】
本発明に係るベルト滑り予測装置は、入力回転速度、出力回転速度、入力トルク、ベルト挟圧力に基づいて、第１のプーリにおけるベルトμ勾配をより精度良く推定し、ベルトμ勾配に基づいてベルト滑りを予測することにより、早期にベルト滑りを予測することができる。
【発明を実施するための最良の形態】
【００２３】
以下、発明を実施するための最良の形態について図面を参照しながら詳細に説明する。
【００２４】
本発明に係るベルト滑り予測装置は、入力軸側のプライマリプーリ（以下「プライマリ」と省略する。）、出力軸側のセカンダリプーリ（以下「セカンダリ」と省略する。）、及びプライマリとセカンダリの間に掛けられたベルトを備えた無段変速機（ＣＶＴ）において、このベルトの滑りを予測するものである。なお、各プーリは、固定シーブと、この固定シーブに向けてシーブ圧に応じて付勢されている可動シーブで構成されている。
【００２５】
最初に発明の原理について説明し、次に第１及び第２の実施形態について説明する。
【００２６】
［発明の原理］
１．μ勾配を用いたベルトμ特性の表現
図１は、無段変速機のプライマリ又はセカンダリに掛けられたベルトのベルトμ特性を示す図である。ベルトμ特性は飽和特性をもつ。ベルト滑り速度に対するベルトμの接線勾配（μ勾配）ｋは、マクロスリップ限界に近づくに従って小さくなる（０近傍の値となる。）。したがって、μ勾配ｋを検出することができれば、μ勾配ｋの０接近時にマクロスリップ限界が検出されて、ベルト滑りを予測できる。
【００２７】
図１の動作点におけるベルトμを（１）式で表す。
【００２８】
【数１】
2006002795.tif 000003
【００２９】
ここで、（１）式のパラメータは、
ｋ  ：動作点におけるベルト滑り速度に対する
ベルトμの接線勾配（μ勾配）［ｓ／ｍ］
Δｖ：動作点における滑り速度［ｍ／ｓ］
μ０：動作点におけるμ切片
Ｖ  ：ベルト速度
ω  ：プーリ回転速度
Ｒ  ：掛かり径
である。
【００３０】
２．μ勾配を用いたプーリ運動方程式
μ勾配は、プライマリプーリ、セカンダリプーリの各プーリ回転運動方程式と検出信号とに基づく同定により検出される。
【００３１】
プライマリ、セカンダリのプーリ回転運動方程式は、（２）式及び（３）式となる。
【００３２】
【数２】
2006002795.tif 000004
【００３３】
ただし、時間微分は“・”で表した。また、添え字の“１”はプライマリ、添え字の“２”はセカンダリであることを表す。（２）式及び（３）式のパラメータは次の通りである。
【００３４】
Ｊ1，Ｊ2：回転慣性［ｋｇ・ｍ2］
ω1，ω2：回転角速度［ｒａｄ／ｓ］
μ1，μ2：動作点におけるベルトμ
Ｆ1，Ｆ2：ベルト挟圧力［Ｎ］
Ｒ1，Ｒ2：ベルト掛かり径［ｍ］
Ｔｉｎ：入力トルク［Ｎｍ］
Ｔｏｕｔ：出力トルク［Ｎｍ］
θ：プーリ傾角
（２）、（３）式中、計測可能な信号は回転角速度ω1，ω2、ベルト挟圧力Ｆ1，Ｆ2である。角加速度ｄω1／ｄｔ、ｄω2／ｄｔは、回転角速度ω1，ω2の近似微分することで算出される。ベルト掛かり径Ｒ1，Ｒ2は、 減速比γ（＝ω1／ω2）に基づいて算出される。入力トルクＴｉｎは、流入空気量及び車速等から算出されるエンジントルクを使用すれば、算出される。出力トルクＴｏｕｔは、入力トルクＴｉｎ、減速比γ、伝達効率ηを用いて、（４）式で算出される。
【００３５】
【数３】
2006002795.tif 000005
【００３６】
伝達効率η＝１．０とみなし、（１）及び（４）式を用いて（２）及び（３）式を書き換えると、（５）及び（６）式となる。
【００３７】
【数４】
2006002795.tif 000006
【００３８】
（５）及び（６）式中、“~”は未知パラメータであることを示す。未知パラメータである滑り速度ΔＶ1、ΔＶ2は、（７）及び（８）式で表される。
【００３９】
【数５】
2006002795.tif 000007
【００４０】
（７）、（８）式において、ベルト速度Ｖを未知パラメータとして取り扱い、同定することにより、直線と回転を含むベルトの運動の記述が不要となる。その結果、モデル精度の低下が抑えられ、μ勾配推定精度向上につながる。（７）及び（８）式を用いて（５）及び（６）式を書き換えると、（９）及び（１０）式となる。
【００４１】
【数６】
2006002795.tif 000008
【００４２】
３．μ勾配の同定方法
（９）及び（１０）式に基づきμ勾配を同定する。（９）及び（１０）式において、既知の入力信号、出力信号及び未知パラメータをまとめた式は、（１１）式のように表現される。
【００４３】
【数７】
2006002795.tif 000009
【００４４】
ただし、
ｙ［ｉ］：既知の出力信号
ξ［ｉ］：既知の入力信号
τ［ｉ］：未知パラメータ
である。（１１）式のｉは更新周期の刻みである。（９）式におけるプライマリ側のｙ、τ、ξは、各々（１２）、（１３）及び（１４）式となる。
【００４５】
【数８】
2006002795.tif 000010
【００４６】
同様にして、セカンダリ側のｙ、τ、ξは、各々（１５）、（１６）及び（１７）式となる。
【００４７】
【数９】
2006002795.tif 000011
【００４８】
そして、（１２）から（１４）式を（１１）式に代入し、例えば重み付け最小自乗法を用いてオンライン同定することでτ、すなわちプライマリ側のベルトμ勾配ｋ1を求めることができる。同様に、（１５）から（１７）式を（１１）式に代入し、例えば重み付け最小自乗法を用いてオンライン同定することでτ、すなわちセカンダリ側のベルトμ勾配ｋ2を同定できる。なお、両方のベルトμ勾配ｋ1、ｋ2を同定してもよいし、いずれか一方を同定してもよい。
【００４９】
上記演算は、ベルト滑りにより時々刻々変化するある一点での動作点の各状態量を用い、２つのサンプル点間での状態量を用いないため、従来法のようなトルクと掛かり径一定といった仮定を設ける必要がないため、μ勾配の推定精度を向上できる。
【００５０】
（１４）及び（１７）式のτに含まれる未知パラメータは、各々２個である。すなわち、エンジントルク（入力トルクＴｉｎ）を使用することで、τに含まれる未知パラメータを２つにすることができ、この結果、同定演算の遅れを大きく低減することができる。その理由は、次の通りである。
【００５１】
セカンダリ側のベルトμ勾配ｋ2を求めるための式、つまり（１１）式、（１５）から（１７）式は２行の行列式となる。１サンプリング時間が例えば１０［ｍｓ］の場合、この行列式に必要なパラメータを取得するための時間は２０［ｍｓ］だけである。なお、プライマリ側のベルトμ勾配ｋ1を求めるための式、つまり（１１）から（１４）式も同様に２行の行列式となる。したがって、本発明は、演算時間をかけることなくベルトμ勾配ｋ1，ｋ2を求めることが可能になる。
【００５２】
そして、本発明は、同定したμ勾配ｋ1あるいはｋ2の絶対値が、所定のしきい値以下となったときに、ベルト滑りを予測する。これにより、ベルト滑りを早期に予測することができる。
【００５３】
［実施形態］
図２は、本発明の実施の形態に係るベルト滑り予測装置の構成を示すブロック図である。
【００５４】
ベルト滑り予測装置は、プライマリの回転速度に応じて回転パルス信号を生成するプライマリ回転速度センサ１１と、セカンダリの回転速度に応じて回転パルス信号を生成するセカンダリ回転速度センサ１２と、エンジンに流入する流入空気量に応じたセンサ信号を生成する流入空気量センサ１３と、車輪の回転角速度に応じたセンサ信号を生成する車輪速センサ１４と、プライマリのシーブ圧を検出するプライマリシーブ圧センサ１５と、セカンダリのシーブ圧を検出するセカンダリシーブ圧センサ１６と、各センサ信号に基づいてＣＶＴのベルト滑りを予測する電子制御ユニット（ＥＣＵ）２０と、を備えている。
【００５５】
図３は、ＥＣＵ２０の構成を示すブロック図である。
【００５６】
ＥＣＵ２０は、プライマリ回転速度検出回路２１と、セカンダリ回転速度検出回路２２と、エンジントルク推定回路２３と、プライマリ挟圧力検出回路２４と、セカンダリ挟圧力検出回路２５と、を備えている。
【００５７】
さらに、ＥＣＵ２０は、プライマリ回転角加速度算出回路２６と、減速比算出回路２７と、セカンダリ回転角加速度算出回路２８と、しきい値生成回路２９と、掛かり径算出回路３０と、伝達効率マップ参照回路３１と、μ勾配算出回路３２と、ベルト滑り予測回路３３と、を備えている。
【００５８】
プライマリ回転速度検出回路２１は、プライマリ回転速度センサ１１からの回転パルス信号に基づいてプライマリ回転角速度ω1［ｒａｄ／ｓ］を検出し、このプライマリ回転角速度ω1をプライマリ回転角加速度算出回路２６、減速比算出回路２７、μ勾配算出回路３２に供給する。
【００５９】
セカンダリ回転速度検出回路２２は、セカンダリ回転速度センサ１２からの回転パルス信号に基づいてセカンダリ回転角速度ω2［ｒａｄ／ｓ］を検出し、このセカンダリ回転角速度ω2を減速比算出回路２７、セカンダリ回転角加速度算出回路２８、μ勾配算出回路３２に供給する。
【００６０】
エンジントルク推定回路２３は、エンジンへの流入空気量及び車速に基づいてエンジントルクを推定し、エンジントルクの推定値をμ勾配算出回路３２に供給する。このエンジントルクの推定値は、ＣＶＴへの入力トルクの値となる。
【００６１】
プライマリ挟圧力検出回路２４は、プライマリシーブ圧センサ１５のセンサ信号に基づいてプライマリのシーブ圧を検出し、これを用いてプライマリのベルト挟圧力Ｆ1を検出する。そして、検出したベルト挟圧力Ｆ1をμ勾配算出回路３２に供給する。
【００６２】
セカンダリ挟圧力検出回路２５は、セカンダリシーブ圧センサ１６のセンサ信号に基づいてセカンダリのシーブ圧を検出し、これを用いてセカンダリのベルト挟圧力Ｆ2を検出する。そして、検出したベルト挟圧力Ｆ2をμ勾配算出回路３２に供給する。
【００６３】
プライマリ回転角加速度算出回路２６は、プライマリ回転角速度ω1を時間微分することでプライマリ回転角加速度を算出し、このプライマリ回転角加速度をμ勾配算出回路３２に供給する。
【００６４】
減速比算出回路２７は、プライマリ回転角速度ω1及びセカンダリ回転角速度ω2に基づいて減速比γ（＝ω1／ω2）を演算し、減速比γをしきい値生成回路２９、掛かり径算出回路３０及び伝達効率マップ参照回路３１に供給する。
【００６５】
セカンダリ回転角加速度算出回路２８は、セカンダリ回転角速度ω2を時間微分することでセカンダリ回転角加速度を算出し、このセカンダリ回転角加速度をμ勾配算出回路３２に供給する。
【００６６】
しきい値生成回路２９は、減速比γの様々な値に対応するしきい値を表したマップを記憶している。しきい値は、例えば、マクロスリップ発生直前におけるベルトμ勾配の値を示している。しきい値生成回路２９は、このマップを参照して、減速比算出回路２７で算出された減速比γに対応するしきい値を求め、求めたしきい値をベルト滑り予測回路３３に供給する。なお、しきい値生成回路２９は、ベルト滑り予測回路３３がベルト滑りを予測したときは、チャタリングを防止するために、しきい値を上げるようになっている。
【００６７】
掛かり径算出回路３０は、減速比算出回路２７で算出された減速比γを用いて、（１８）式に従ってプライマリのベルト掛かり径Ｒ1、（１９）式に従ってセカンダリのベルト掛かり径Ｒ2を算出し、ベルト掛かり径Ｒ1，Ｒ2をμ勾配算出回路３２に供給する。
【００６８】
【数１０】
2006002795.tif 000012
【００６９】
ただし、ａ，ｂ，ｃはそれぞれ定数である。
【００７０】
伝達効率マップ参照回路３１は、減速比γの様々な値に対応する伝達効率を表した伝達効率マップを記憶している。伝達効率とは、プライマリからセカンダリへのトルクの伝達効率を表す係数である。伝達効率マップ参照回路３１は、この伝達効率マップを参照して、減速比算出回路２７で算出された減速比γに対応する伝達効率ηを求め、求めた伝達効率ηをベルト滑り予測回路３３に供給する。なお、本実施形態では、伝達効率マップ参照回路３１で求められる伝達効率ηは１とするが、伝達効率ηはこれに限定されるものではない。
【００７１】
μ勾配算出回路３２は、プライマリ回転角速度ω1、セカンダリ回転角速度ω2、ベルト掛かり径Ｒ1，Ｒ2、プライマリ回転角加速度、セカンダリ回転角加速度、伝達効率η、エンジントルク（入力トルク）Ｔｉｎ、プライマリベルト挟圧力Ｆ1、セカンダリベルト挟圧力Ｆ2に基づいて、ベルトμ勾配を算出する。
【００７２】
具体的には、μ勾配算出回路３２は、各回路から供給されたパラメータを（１２）式のｙ、（１３）式のξにそれぞれ代入する。そして、（１１）式に対して最小自乗法用いて、時々刻々（１４）式のτを同定する。また、μ勾配算出回路３２は、各回路から供給されたパラメータを（１５）式のｙ、（１６）式のξにそれぞれ代入する。そして、（１１）式に対して最小自乗法用いて、時々刻々（１７）式のτを同定する。
【００７３】
このように、μ勾配算出回路３２は、τ（ｉ）を同定することによって、プライマリのベルトμ勾配ｋ1、セカンダリのベルトμ勾配ｋ2を求め、そしてベルトμ勾配ｋ1，ｋ2をベルト滑り予測回路３３に供給する。
【００７４】
ベルト滑り予測回路３３は、ベルトμ勾配ｋ1の絶対値、ベルトμ勾配ｋ2の絶対値のそれぞれと、しきい値生成回路２９で生成されたしきい値とを比較して、いずれか一方のベルトμ勾配の絶対値がしきい値以下となったときに、ベルト滑りを予測する。これにより、ベルト滑り予測回路３３は、プライマリ側、セカンダリ側のいずれであっても、ベルト滑りを予測できる。
【００７５】
ベルト滑り時はｋ1，ｋ2ともに低下するため、μ勾配算出回路３２はベルトμ勾配ｋ1，ｋ2のいずれか一方だけを算出してもよい。そして、ベルト滑り予測回路３３は、μ勾配算出回路３２で算出された１つのベルトμ勾配ｋとしきい値とを比較して、ベルト滑りを予測すればよい。
【００７６】
図４（Ａ）はプライマリへの入力トルクＴｉｎ［Ｎｍ］、（Ｂ）はプライマリ及びセカンダリの回転数［ｒｐｍ］、（Ｃ）は減速比γ、（Ｄ）はセカンダリのベルトμ勾配ｋ2の経時変化を示す図である。また、ベルト滑りの予兆となるマクロスリップ発生点も示した。
【００７７】
図４では、定常運転中、挟圧力一定の状態で入力トルクをあげていき、マクロスリップさせたときのベルトμ勾配ｋ2の同定結果を表している。本法では、マクロスリップ付近でのベルトμ勾配ｋ2の立ち下がりが早く、ベルト滑りを精度良く早期に検出する性能が向上している。
【００７８】
以上のように、本発明の実施の形態に係るベルト滑り予測装置は、ベルトμ勾配を用いたプーリ回転運動方程式、プライマリ回転速度、セカンダリ回転速度、エンジントルク、プライマリ又はセカンダリのベルト挟圧力に基づいて、ある動作点上でのベルトμ勾配ｋ1又はｋ2を推定する。そして、ベルトμ勾配ｋ1，ｋ2としきい値と比較することで、ベルト滑りを精度良く早期に予測することができる。
【００７９】
特に、ベルト滑り時にはかかり径変化、ベルトμ勾配特性変化による入力トルク変化を生じるが、従来法は掛かり径と入力トルクが一定の条件を設けるため、ベルトμ勾配推定精度の低下が避けられなかった。これに対して本方法は上記の変化を盛り込んだ形でモデル表現するため、従来法と比べてベルトμ勾配推定精度が向上する利点がある。
【００８０】
ベルト滑り（マクロスリップ発生）直前は、ＣＶＴ伝達効率の最大点となる。よって、ベルト滑り予測装置は、ベルト滑りを予測すると同時に、ＣＶＴ伝達効率の最大点も検出することができる。
【００８１】
また、上記ベルト滑り予測装置は、エンジントルクを用いることで、プーリ回転運動方程式における未知パラメータの数を少なくできる。このため、ベルトμ勾配を非常に短時間で算出して、ベルト滑りを早期に予測することができる。
【００８２】
なお、本発明は、上述した実施の形態に限定されるものではなく、特許請求の範囲に記載された範囲内で設計上の変更をされたものにも適用可能であるのは勿論である。
【００８３】
例えば、上記実施形態に係るベルト滑り予測装置は、ベルトμ勾配ｋ1，ｋ2の両方を算出したが、いずれか一方のみを算出してもよい。ベルト滑り予測装置は、ベルトμ勾配ｋ1のみを算出する場合は、セカンダリシーブ圧センサ１６及びセカンダリ挟圧力検出回路２５を設けてなくてもよい。一方、ベルト滑り予測装置は、ベルトμ勾配ｋ2のみを算出する場合は、プライマリシーブ圧センサ１５及びプライマリ挟圧力検出回路２４を設けてなくてもよい。
【００８４】
また例えば、上記実施の形態では、「ベルト滑り速度」に対するベルトμの接線勾配（ベルトμ勾配ｋ）を求めたが、「ベルト滑り率」に対するベルトμの接線勾配を求めてもよい。
【００８５】
さらに、ベルト滑り予測装置は、例えばＣＶＴ油圧制御バルブのようなベルト挟圧力を制御するベルト挟圧力制御手段を設けてもよい。そして、ＥＣＵ２０がベルトμ勾配をしきい値付近に維持するようにＣＶＴ油圧制御バルブを制御すればよい。これにより、ＣＶＴ伝達効率をほぼ最大に維持できるので、燃費を向上させることができる。
【図面の簡単な説明】
【００８６】
【図１】無段変速機のベルトμ特性を示す図である。
【図２】本発明の実施の形態に係るベルト滑り予測装置の構成を示すブロック図である。
【図３】ＥＣＵの構成を示すブロック図である。
【図４】（Ａ）はプライマリへの入力トルクＴｉｎ［Ｎｍ］、（Ｂ）はプライマリ及びセカンダリの回転数［ｒｐｍ］、（Ｃ）は減速比γ、（Ｄ）はセカンダリのベルトμ勾配ｋ２の経時変化を示す図である。
【符号の説明】
【００８７】
２０  ＥＣＵ
２１  プライマリ回転速度検出回路
２２  セカンダリ回転速度検出回路
２３  エンジントルク推定回路
２４  プライマリ挟圧力検出回路
２５  セカンダリ挟圧力検出回路
２６  プライマリ回転角加速度算出回路
２７  減速比算出回路
２８  セカンダリ回転角加速度算出回路
２９  しきい値生成回路
３０  掛かり径算出回路
３１  伝達効率マップ参照回路
３２  μ勾配算出回路
３３  ベルト滑り予測回路
【図１】
2006002795.tif 000013
【図２】
2006002795.tif 000014
【図３】
2006002795.tif 000015
【図４】
2006002795.tif 000016