(57)【要約】【課題】タイヤにスリップ率が与えられて接地面に滑り域を形成するときのタイヤの摩耗特性を予測する際、トレッド部材におけるゴム物性の評価のみならず、トレッド部材の路面に対する滑り量を用いてタイヤの摩耗特性を予測する。【解決手段】タイヤにスリップ率を与えたときのタイヤ回転軸に作用するタイヤ軸力及びセルフアライニングトルクの少なくとも一方の、スリップ率依存性の特性曲線を取得する。次に、この特性曲線から、複数のタイヤ力学要素パラメータを用いて構成されるタイヤ力学モデルに基づいてタイヤ力学要素パラメータの値を導出する。さらに、導出したタイヤ力学要素パラメータの値をタイヤ力学モデルに与えて、タイヤ接地面の凝着域と滑り域を求め、この滑り域からタイヤ滑り量を求める。このタイヤ滑り量を、トレッドゴム部材の摩耗特性データと共に用いて、スリップ率が与えられたときのタイヤのトレッド部材の摩耗特性を評価する。【選択図】図１０
2006232011.tif 000002
【特許請求の範囲】
【請求項１】
タイヤにスリップ率が与えられて接地面に滑り域を形成するときのタイヤ滑り量を求めることにより、タイヤの摩耗特性を予測するタイヤの摩耗予測方法であって、
タイヤにスリップ率を与えたときのタイヤ回転軸に作用するタイヤ軸力の、前記スリップ率に対して変化する特性曲線を取得するステップと、
前記特性曲線から、複数のタイヤ力学要素パラメータを用いて構成されるタイヤ力学モデルに基づいて、前記特性曲線を定めるタイヤ力学要素パラメータの値を導出するステップと、
前記タイヤ力学要素パラメータの値を前記タイヤ力学モデルに与えて、タイヤの接地面に形成される凝着域と滑り域を求め、この滑り域からタイヤ滑り量を求めるステップと、
このタイヤ滑り量を、トレッドゴム部材の摩耗特性データと共に用いて、前記スリップ率が与えられたときのタイヤのトレッド部材の摩耗特性を予測するステップと、を有することを特徴とするタイヤの摩耗予測方法。
【請求項２】
前記スリップ率は、タイヤにスリップ角が与えられたときのスリップ率及び制駆動方向のスリップ率の少なくとも一方を含む請求項１に記載のタイヤの摩耗予測方法。
【請求項３】
前記タイヤ軸力は、前記タイヤにスリップ角が与えられて、タイヤ回転軸に対して平行な方向に作用する横力であり、
前記特性曲線を取得するステップでは、前記横力の特性曲線の他に、前記横力により生じるセルフアライニングトルクのスリップ角依存性を表す特性曲線を取得する請求項２に記載のタイヤの摩耗予測方法。
【請求項４】
前記タイヤ力学モデルは、タイヤにスリップ角が与えられたときの横力を算出するとともに、セルフアライニングトルクをタイヤの接地面に作用する横力によって生じる横力トルク成分と、タイヤの接地面に作用する前後力によって生じる前後力トルク成分とに分けてセルフアライニングトルクを算出するモデルである請求項３に記載のタイヤの摩耗予測方法。
【請求項５】
前記タイヤ力学要素パラメータの値を導出する際、前記横力の特性曲線と前記タイヤ力学モデルで算出される横力の対応する曲線との二乗残差和と、前記セルフアライニングトルクの特性曲線と前記タイヤ力学モデルで算出されるセルフアライニングトルクの対応する曲線との二乗残差和とを、重み付け係数を用いて重み付け加算した値であって、前記重み付け係数として、前記横力および前記セルフアライニングトルクのそれぞれの特性曲線の、スリップ角に依存して変化する値のばらつきの情報から求められる係数を用いた複合二乗残差和の値が、所定値以下となるように、前記タイヤ力学要素パラメータの値を導出する請求項３又は４に記載のタイヤの摩耗予測方法。
【請求項６】
前記タイヤに対してスリップ角及び制駆動方向のスリップ率が与えられ、前記タイヤ軸力は、前記タイヤ回転軸に対して平行な方向に作用する横力及び前記タイヤ回転軸に対して直交する方向に作用する前後力であり、
前記特性曲線を取得するステップでは、前記横力のスリップ角依存性の特性曲線のほかに、前記横力により生じるセルフアライニングトルクのスリップ角依存性を表す特性曲線と、前記前後力のスリップ率依存性を表す特性曲線を取得する請求項２に記載のタイヤの摩耗予測方法。
【請求項７】
導出された前記タイヤ力学要素パラメータの値を用いて、所定のスリップ角及び制駆動方向のスリップ率におけるタイヤ滑り量を、前記タイヤ力学モデルに基づいて算出し、
算出したタイヤ滑り量を用いて、前記所定のスリップ角及び制駆動方向のスリップ率におけるタイヤのトレッド摩耗特性を予測するステップ、をさらに有する請求項６に記載のタイヤの摩耗予測方法。
【請求項８】
前記タイヤ力学要素パラメータの値を導出する際、前記前後力の特性曲線と前記タイヤ力学モデルで算出される前後力の対応する曲線との二乗残差和と、前記横力の特性曲線と前記タイヤ力学モデルで算出される横力の対応する曲線との二乗残差和と、前記セルフアライニングトルクの特性曲線と前記タイヤ力学モデルで算出されるセルフアライニングトルクの対応する曲線との二乗残差和とを、重み付け係数を用いて重み付け加算した値であって、前記重み付け係数として、前記横力および前記セルフアライニングトルクのそれぞれの特性曲線の、スリップ角に依存して変化する値のばらつきの情報から求められる係数を用いた複合二乗残差和の値が、所定値以下となるように、前記タイヤ力学要素パラメータの値を導出する請求項６又は７に記載のタイヤの摩耗予測方法。
【請求項９】
前記特性曲線から、前記タイヤ力学モデルに基づいて、前記タイヤ力学要素パラメータの値を導出する際、セルフアライニングトルクにより発生するタイヤの捩じり変形によって、付与されるスリップ角が修正された実効スリップ角を用いてタイヤ力学要素パラメータの値を導出する請求項３〜８のいずれか１項に記載のタイヤの摩耗予測方法。
【請求項１０】
導出される前記タイヤ力学要素パラメータの値は、タイヤのトレッド部材と路面との間の凝着摩擦係数およびすべり摩擦係数と接地圧分布の形状を規定する形状規定係数を含む請求項１〜９のいずれか１項に記載のタイヤの摩耗予測方法。
【請求項１１】
前記凝着摩擦係数、前記すべり摩擦係数および前記形状規定係数は、予め求められた、タイヤの剪断変形に対する剛性パラメータ、タイヤの横曲げ変形に対する剛性パラメータおよびタイヤの捩じり変形に対する剛性パラメータの少なくとも1つを用いて導出する請求項１０に記載のタイヤの摩耗予測方法。
【請求項１２】
前記タイヤに対して制駆動方向のスリップ率が与えられ、前記タイヤ軸力は、前記タイヤ回転軸に対して直交する方向に作用する前後力であり、
前記特性曲線を取得するステップでは、前記前後力のスリップ率依存性を表す特性曲線を取得する請求項２に記載のタイヤの摩耗予測方法。
【請求項１３】
請求項１〜１２のいずれか１項に記載のタイヤの摩耗予測方法を用いてトレッド部材の摩耗特性を予測するステップと、
トレッド部材の摩耗特性の予測結果が目標を達成しない場合、タイヤ力学要素パラメータの値を修正するステップと、
トレッド部材の摩耗特性の予測結果が目標を達成する場合、設定されたタイヤ力学要素パラメータの値をタイヤ設計仕様特性として決定するステップと、を有することを特徴とするタイヤの設計方法。
【請求項１４】
修正する前記タイヤ力学要素パラメータは、前記タイヤ力学モデルにおける剛性を表す剛性パラメータを含み、前記タイヤ設計仕様特性として決定したタイヤ力学要素パラメータの値に基づいて、タイヤの構造部材の設計を行うステップを、さらに有する請求項１３に記載のタイヤの設計方法。
【請求項１５】
修正する前記タイヤ力学要素パラメータは、前記タイヤ力学モデルにおける凝着摩擦係数又は滑り摩擦係数を含み、前記タイヤ設計仕様特性として決定したタイヤ力学要素パラメータの値に基づいて、タイヤのゴム部材の材料設計を行うステップを、さらに有する請求項１３又は１４に記載のタイヤの設計方法。
【請求項１６】
請求項１３〜１５のいずれか１項に記載のタイヤの設計方法により定められたタイヤの構成部材又はタイヤのゴム部材を用いて、タイヤを製造することを特徴とするタイヤの製造方法。
【請求項１７】
タイヤにスリップ率が与えられて接地面に滑り域を形成するときのタイヤ滑り量を求めることにより、タイヤの摩耗特性を予測するタイヤの摩耗評価予測システムであって、
タイヤにスリップ率を与えたときのタイヤ回転軸に作用するタイヤ軸力の、前記スリップ率に対して変化する特性曲線を取得する手段と、
前記タイヤ軸力の特性曲線から、複数のタイヤ力学要素パラメータを用いて構成されるタイヤ力学モデルに基づいて、前記特性曲線を定めるタイヤ力学要素パラメータの値を導出する手段と、
前記タイヤ力学要素パラメータの値を前記タイヤ力学モデルに与えて、タイヤの接地面に形成される凝着域と滑り域を求め、この滑り域からタイヤ滑り量を求める手段と、
このタイヤ滑り量を、トレッドゴム部材の摩耗特性データと共に用いて、前記スリップ率が与えられたときのタイヤのトレッド摩耗特性を予測する手段と、を有することを特徴とするタイヤの摩耗予測システム。
【請求項１８】
タイヤにスリップ率が与えられて接地面に滑り域を形成するときのタイヤ滑り量を求めることによりタイヤの摩耗特性を予測させるための、コンピュータが実行可能なプログラムであって、
タイヤにスリップ率を与えたときのタイヤ回転軸に作用するタイヤ軸力の、前記スリップ率に対して変化する特性曲線を取得し、コンピュータのメモリに記憶させる手順と、
前記メモリに記憶された前記タイヤ軸力の特性曲線から、複数のタイヤ力学要素パラメータを用いて構成されるタイヤ力学モデルに基づいて、前記特性曲線を定めるタイヤ力学要素パラメータの値を、前記コンピュータの演算手段に導出させる手順と、
導出された前記タイヤ力学要素パラメータの値を前記タイヤ力学モデルに与えて、タイヤの接地面に形成される凝着域と滑り域を前記演算手段に求めさせ、この滑り域からタイヤ滑り量を前記演算手段に求めさせる手順と、
このタイヤ滑り量を、トレッドゴム部材の摩耗特性データと共に用いて、前記スリップ率が与えられたときのタイヤのトレッド摩耗特性を、前記演算手段に予測させる手順と、を有することを特徴とするプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、タイヤにスリップ率が与えられて接地面に滑り域を形成するときのタイヤの摩耗特性を予測するタイヤの摩耗予測方法、この方法を用いてタイヤを設計するタイヤの設計方法、この設計方法で設計されたタイヤを製造するタイヤの製造方法、前記タイヤの摩耗予測方法を行うタイヤの摩耗予測システム及びこれを実行させるプログラムに関する。
【背景技術】
【０００２】
車両に装着されるタイヤは、路面に接触して転動するため、走行とともにトレッド部材が摩耗する消耗品である。このため、タイヤのトレッド部材は、他性能（操縦性能や振動乗心地性能）を犠牲にしない範囲で摩耗特性が良好なタイヤが望まれている。
一般に、タイヤのトレッド部材の摩耗特性、例えばトレッド部材の摩耗量は路面に大して部分的に滑るトレッド部材の滑り量と、トレッド部材におけるゴムの耐摩耗物性とによって決定される。
一般に、トレッド部材の使用条件に対する摩耗割合は、大きい順に並べるとコーナリング時、制動時、駆動時となっていることが知られている。
【０００３】
ところで、摩耗割合が大きいコーナリング時には、トレッド部材の路面に対する滑り量がスリップ角によって一意的に定まるものではなく、トレッド部材の凝着摩擦係数、滑り摩擦係数のみならず、ベルト部分の構造やサイド部分の構造の影響を受けて定まる。  
又、摩耗割合が比較的大きい制動時においても、トレッド部材の路面に対する滑り量は、制動時のスリップ率によって一意的に定まるものではなく、トレッド部材の凝着摩擦係数、滑り摩擦係数のみならず、ベルト部分の構造やサイド部分の構造の影響を受けて定まる。
【０００４】
一方、タイヤの摩耗特性を予測する場合、トレッド部材におけるゴムの破断物性やゴムの室内ラボ試験（摩耗量測定試験）結果で評価される。室内ラボ試験として、例えば、新ＪＩＳ  Ｋ６２６４に記されている各種試験法が用いられる。
しかし、上述したようにタイヤにおけるトレッド部材の摩耗特性は、トレッド部材のゴム物性のみによって左右されるものではない。摩耗特性は、トレッド部材の凝着摩擦係数、滑り摩擦係数の他に、ベルト部分の構造やサイド部分の構造の影響を受けるものであり、ゴム物性の評価だけではタイヤの摩耗特性を精度よく予測することはできない。又、実車による摩耗試験を行うことで、摩耗特性を精度良く予測することはできるが、試験に時間がかかり過ぎるといった問題があった。
【発明の開示】
【発明が解決しようとする課題】
【０００５】
そこで、本発明は、タイヤにスリップ率が与えられて接地面に滑り域を形成するときのタイヤの摩耗特性を予測する際、トレッド部材におけるゴム物性の評価のみならず、トレッド部材の路面に対する滑り量を用いてタイヤの摩耗特性を予測するタイヤの摩耗予測方法、この予測を用いて行うタイヤの設計方法、このタイヤの設計方法を用いて行うタイヤの製造方法、及びタイヤの摩耗特性を予測するタイヤの摩耗予測システム及びこのシステムを実行させるプログラムを提供することを目的とする。
【課題を解決するための手段】
【０００６】
本発明は、タイヤにスリップ率が与えられて接地面に滑り域を形成するときのタイヤ滑り量を求めることにより、タイヤの摩耗特性を予測するタイヤの摩耗予測方法であって、タイヤにスリップ率を与えたときのタイヤ回転軸に作用するタイヤ軸力の、前記スリップ率に対して変化する特性曲線を取得するステップと、前記特性曲線から、複数のタイヤ力学要素パラメータを用いて構成されるタイヤ力学モデルに基づいて、前記特性曲線を定めるタイヤ力学要素パラメータの値を導出するステップと、前記タイヤ力学要素パラメータの値を前記タイヤ力学モデルに与えて、タイヤの接地面に形成される凝着域と滑り域を求め、この滑り域からタイヤ滑り量を求めるステップと、このタイヤ滑り量を、トレッドゴム部材の摩耗特性データと共に用いて、前記スリップ率が与えられたときのタイヤのトレッド部材の摩耗特性を予測するステップと、を有することを特徴とするタイヤの摩耗予測方法を提供する。
【０００７】
ここで、前記スリップ率は、タイヤにスリップ角が与えられたときのスリップ率及び制駆動方向のスリップ率の少なくとも一方を含むことが好ましい。
前記タイヤ軸力は、前記タイヤにスリップ角が与えられて、タイヤ回転軸に対して平行な方向に作用する横力であり、この場合、前記特性曲線を取得するステップでは、前記横力の特性曲線の他に、前記横力により生じるセルフアライニングトルクのスリップ角依存性を表す特性曲線を取得することが好ましい。
その際、前記タイヤ力学モデルは、タイヤにスリップ角が与えられたときの横力を算出するとともに、セルフアライニングトルクをタイヤの接地面に作用する横力によって生じる横力トルク成分と、タイヤの接地面に作用する前後力によって生じる前後力トルク成分とに分けてセルフアライニングトルクを算出するモデルであることが好ましい。
又、前記タイヤ力学要素パラメータの値を導出する際、前記横力の特性曲線と前記タイヤ力学モデルで算出される横力の対応する曲線との二乗残差和と、前記セルフアライニングトルクの特性曲線と前記タイヤ力学モデルで算出されるセルフアライニングトルクの対応する曲線との二乗残差和とを、重み付け係数を用いて重み付け加算した値であって、前記重み付け係数として、前記横力および前記セルフアライニングトルクのそれぞれの特性曲線の、スリップ角に依存して変化する値のばらつきの情報から求められる係数を用いた複合二乗残差和の値が、所定値以下となるように、前記タイヤ力学要素パラメータの値を導出することが好ましい。
【０００８】
又、前記タイヤに対してスリップ角及び制駆動方向のスリップ率が与えられ、前記タイヤ軸力は、前記タイヤ回転軸に対して平行な方向に作用する横力及び前記タイヤ回転軸に対して直交する方向に作用する前後力であり、前記特性曲線を取得するステップでは、前記横力のスリップ角依存性の特性曲線のほかに、前記横力により生じるセルフアライニングトルクのスリップ角依存性を表す特性曲線と、前記前後力のスリップ率依存性を表す特性曲線を取得することも好ましい。
その際、導出された前記タイヤ力学要素パラメータの値を用いて、所定のスリップ角及び制駆動方向のスリップ率におけるタイヤ滑り量を、前記タイヤ力学モデルに基づいて算出し、算出したタイヤ滑り量を用いて、前記所定のスリップ角及び制駆動方向のスリップ率におけるタイヤのトレッド摩耗特性を予測するステップ、をさらに有することが好ましい。前記タイヤ力学要素パラメータの値を導出する際、前記前後力の特性曲線と前記タイヤ力学モデルで算出される前後力の対応する曲線との二乗残差和と、前記横力の特性曲線と前記タイヤ力学モデルで算出される横力の対応する曲線との二乗残差和と、前記セルフアライニングトルクの特性曲線と前記タイヤ力学モデルで算出されるセルフアライニングトルクの対応する曲線との二乗残差和とを、重み付け係数を用いて重み付け加算した値であって、前記重み付け係数として、前記横力および前記セルフアライニングトルクのそれぞれの特性曲線の、スリップ角に依存して変化する値のばらつきの情報から求められる係数を用いた複合二乗残差和の値が、所定値以下となるように、前記タイヤ力学要素パラメータの値を導出することが好ましい。
【０００９】
さらに、前記特性曲線から、前記タイヤ力学モデルに基づいて、前記タイヤ力学要素パラメータの値を導出する際、セルフアライニングトルクにより発生するタイヤの捩じり変形によって、付与されるスリップ角が修正された実効スリップ角を用いてタイヤ力学要素パラメータの値を導出することが好ましい。
【００１０】
導出される前記タイヤ力学要素パラメータの値は、タイヤのトレッド部材と路面との間の凝着摩擦係数およびすべり摩擦係数と接地圧分布の形状を規定する形状規定係数を含むことが好ましい。前記凝着摩擦係数、前記すべり摩擦係数および前記形状規定係数は、予め求められた、タイヤの剪断変形に対する剛性パラメータ、タイヤの横曲げ変形に対する剛性パラメータおよびタイヤの捩じり変形に対する剛性パラメータの少なくとも1つを用いて導出することが好ましい。
【００１１】
又、前記タイヤに対して制駆動方向のスリップ率が与えられ、前記タイヤ軸力は、前記タイヤ回転軸に対して直交する方向に作用する前後力であり、前記特性曲線を取得するステップでは、前記前後力のスリップ率依存性を表す特性曲線を取得することが好ましい。
【００１２】
さらに、本発明は、前記タイヤの摩耗予測方法を用いてトレッド部材の摩耗特性を予測するステップと、トレッド部材の摩耗特性の予測結果が目標を達成しない場合、タイヤ力学要素パラメータの値を修正するステップと、トレッド部材の摩耗特性の予測結果が目標を達成する場合、設定されたタイヤ力学要素パラメータの値をタイヤ設計仕様特性として決定するステップと、を有することを特徴とするタイヤの設計方法を提供する。
【００１３】
この場合、修正する前記タイヤ力学要素パラメータは、前記タイヤ力学モデルにおける剛性を表す剛性パラメータを含み、前記タイヤ設計仕様特性として決定したタイヤ力学要素パラメータの値に基づいて、タイヤの構造部材の設計を行うステップを、さらに有することが好ましい。又、修正する前記タイヤ力学要素パラメータは、前記タイヤ力学モデルにおける凝着摩擦係数又は滑り摩擦係数を含み、前記タイヤ設計仕様特性として決定したタイヤ力学要素パラメータの値に基づいて、タイヤのゴム部材の材料設計を行うステップを、さらに有することが好ましい。
【００１４】
さらに、本発明は、前記タイヤの設計方法により定められたタイヤの構成部材又はタイヤのゴム部材を用いて、タイヤを製造することを特徴とするタイヤの製造方法を提供する。
【００１５】
さらに、本発明は、タイヤにスリップ率が与えられて接地面に滑り域を形成するときのタイヤ滑り量を求めることにより、タイヤの摩耗特性を予測するタイヤの摩耗評価予測システムであって、タイヤにスリップ率を与えたときのタイヤ回転軸に作用するタイヤ軸力の、前記スリップ率に対して変化する特性曲線を取得する手段と、前記タイヤ軸力の特性曲線から、複数のタイヤ力学要素パラメータを用いて構成されるタイヤ力学モデルに基づいて、前記特性曲線を定めるタイヤ力学要素パラメータの値を導出する手段と、前記タイヤ力学要素パラメータの値を前記タイヤ力学モデルに与えて、タイヤの接地面に形成される凝着域と滑り域を求め、この滑り域からタイヤ滑り量を求める手段と、このタイヤ滑り量を、トレッドゴム部材の摩耗特性データと共に用いて、前記スリップ率が与えられたときのタイヤのトレッド摩耗特性を予測する手段と、を有することを特徴とするタイヤの摩耗予測システムを提供する。
【００１６】
さらに、本発明は、タイヤにスリップ率が与えられて接地面に滑り域を形成するときのタイヤ滑り量を求めることによりタイヤの摩耗特性を予測させるための、コンピュータが実行可能なプログラムであって、タイヤにスリップ率を与えたときのタイヤ回転軸に作用するタイヤ軸力の、前記スリップ率に対して変化する特性曲線を取得し、コンピュータのメモリに記憶させる手順と、前記メモリに記憶された前記タイヤ軸力の特性曲線から、複数のタイヤ力学要素パラメータを用いて構成されるタイヤ力学モデルに基づいて、前記特性曲線を定めるタイヤ力学要素パラメータの値を、前記コンピュータの演算手段に導出させる手順と、導出された前記タイヤ力学要素パラメータの値を前記タイヤ力学モデルに与えて、タイヤの接地面に形成される凝着域と滑り域を前記演算手段に求めさせ、この滑り域からタイヤ滑り量を前記演算手段に求めさせる手順と、このタイヤ滑り量を、トレッドゴム部材の摩耗特性データと共に用いて、前記スリップ率が与えられたときのタイヤのトレッド摩耗特性を、前記演算手段に予測させる手順と、を有することを特徴とするプログラムを提供する。
【発明の効果】
【００１７】
本発明では、タイヤにスリップ率が与えられて接地面に滑り域を形成するときのタイヤ滑り量を、複数のタイヤ力学要素パラメータを用いて構成されるタイヤ力学モデルに基づいて算出する。導出されるタイヤ力学要素パラメータの値は、特性曲線を精度高く再現するので、タイヤ滑り量も精度高く予測することができる。このため、予測されるタイヤ滑り量を用いて精度高くタイヤ摩耗特性を評価、予測することができる。
【発明を実施するための最良の形態】
【００１８】
以下、添付の図面に示す実施形態に基づいて、本発明のタイヤの摩耗予測方法、タイヤの設計方法、タイヤの製造方法、タイヤの摩耗予測システム及びプログラムを詳細に説明する。
【００１９】
本発明のタイヤの摩耗予測方法の下記に示す実施形態では、まず、タイヤにスリップ角を与えたときのタイヤ回転軸に作用するタイヤ軸力のスリップ角依存性を表す特性曲線と、タイヤ回転軸に作用するセルフアライニングトルクのスリップ角依存性を表す特性曲線とを、種々のスリップ率の下に取得する。これらの特性曲線から、複数のタイヤ力学要素パラメータを用いて構成されるタイヤ力学モデルに基づいて、前記特性曲線を定めるタイヤ力学要素パラメータの値を導出する。次に、このタイヤ力学要素パラメータの値をタイヤ力学モデルに与えて、タイヤの接地面に形成される凝着域と滑り域を求め、この滑り域からタイヤ滑り量を求める。この滑り量により、タイヤの摩耗特性を予測する。すなわち、タイヤにスリップ角と制駆動方向のスリップ率を自在に与えることにより、そのときのトレッド部材の滑り量を求め、この滑り量と、ゴム物性の情報とを用いてタイヤの摩耗特性を予測する。
なお、スリップ角をαとするとｔａｎαはスリップ率であるので、本発明における「スリップ率」は制駆動方向のスリップ率の他に、スリップ角αとしたときのｔａｎαも含む。
以降では、スリップ角と制駆動方向のスリップ率を区別して説明する。
このようなタイヤ滑り量を求めるタイヤ力学モデルは下記に詳述する。
【００２０】
図１は、本発明のタイヤの摩耗予測方法を実施する装置１の構成図である。装置１は、各種プログラムを実行させることで車両の設計方法を行うコンピュータで構成される。
装置１は、前後力Ｆx、横力Ｆｙ、セルフアライニングトルク（以降、単にトルクという）Ｍｚ等のタイヤ軸力・トルクのデータの入力を受けることによって後述するタイヤ力学モデルに基づいて複数のタイヤ力学要素パラメータ（以降、パラメータという）の値を導出する、あるいは、タイヤ力学モデルにおけるタイヤ力学要素パラメータの値の入力を受けることによって、タイヤ力学モデルを用いて前後力、横力およびトルク等のタイヤ軸力・トルクのデータを算出する装置である。
【００２１】
装置１は、コンピュータの各部位及び各プログラムの実行を管理、制御するＣＰＵ２と、バス３を介して各種条件や演算結果を記憶するメモリ４と、各種条件や各種情報を指示入力するマウスやキーボード等の入力操作系５と、入力操作系５をバス３に接続するインターフェース６と、各種条件や情報の入力画面や各種プログラムの処理結果を表示し、プリント出力する出力装置７と、後述する各種プログラムを有し、本装置１の機能を実行するプログラム群８と、を有して構成される。
【００２２】
ここでプログラム群８は、統合・管理プログラム１０、タイヤ力学モデル演算プログラム１２、タイヤ力学要素パラメータ導出プログラム１４、タイヤ軸力・トルクデータ算出プログラム１６、滑り量算出プログラム１８及び摩耗特性予測プログラム２０を有する。
タイヤ力学モデル演算プログラム１２は、タイヤ力学モデルを解析式で表し、タイヤ力学要素パラメータを用いて前後力、横力およびトルク等のタイヤ軸力・トルクを算出するプログラムである。
タイヤ力学モデル演算プログラム１４は、与えられたタイヤ力学要素パラメータの値を用いて、所定のスリップ角及び制駆動方向のスリップ率の条件下、前後力、横力およびトルクの対応算出データ（前後力Ｆx’、横力Ｆｙ'およびトルクＭｚ’）を処理結果として、タイヤ力学要素パラメータ導出プログラム１４、タイヤ軸力・トルクデータ算出プログラム１６に返す演算部である。
【００２３】
図２、図３、図４（ａ）〜（ｃ）、図５（ａ）〜（ｄ）および図６（ａ）〜（ｃ）はタイヤ力学モデルを説明する図である。
タイヤ力学モデルは、図２に示すように、剛体の円筒部材にサイドウォールのばね特性を表す複数のばね要素からなるサイドウォールモデルと、これらのばね要素に接続された弾性リング体からなるベルトモデルと、この弾性リング体の表面に接続されたトレッドモデルを表す弾性要素からなるトレッドモデルとを有して構成されるものである。
【００２４】
タイヤ力学モデルに基づいて導出される力学要素パラメータは、以下のものが例示される。
（ａ）タイヤの横方向、前後方向の剪断剛性によって定められる前後剛性・横剛性Ｋt、
（ｂ）路面とタイヤ間の滑り速度０における滑り摩擦係数μｄ0、
（ｃ）路面とタイヤ間の凝着摩擦係数μｓ、
（ｄ）ベルト部材の横方向曲げ係数ε、
（ｅ）タイヤのタイヤ中心軸周りのねじり剛性の逆数であるねじりコンプライアンス（１／Ｇｍｚ）、
（ｆ）横力発生中の接地面の接地圧力分布を規定する係数ｎ、
（ｇ）接地圧力分布の偏向の程度を表す係数Ｃｑ、
（ｈ）接地面におけるタイヤ中心位置の前後方向への移動の程度を示す移動係数Ｃｘｃ、
（ｉ）横力発生時の実効接地長ｌｅ、
（ｊ）滑り摩擦係数μｄの転動速度依存性係数ｂv、等である。
【００２５】
ここで、前後剛性・横剛性Ｋt、横方向曲げ係数ε、ねじりコンプライアンス（１／Ｇｍｚ）のＧｍｚは、それぞれタイヤの剪断変形に対する剛性パラメータ、横曲げ変形に対する剛性パラメータおよびタイヤの捩じり変形に対する剛性パラメータである。また、横力の発生する方向である横方向とは、タイヤ回転軸の軸線方向を意味する。したがって、タイヤが直進状態で転動する場合の横方向は転動方向に対する左右方向となって方向が一致するが、スリップ角が付いた場合の横方向はタイヤの転動方向に対してスリップ角分ずれる。前後方向とは、タイヤの接地する路面に平行であり、かつタイヤ回転軸の軸線方向に対して直交する方向をいう。また、タイヤ中心軸（図５（ａ），（ｂ）中の軸ＣＬをいう）は、タイヤの転動する回転軸に直交し、かつタイヤの幅方向の中心面を通る、路面に垂直な軸である。
【００２６】
図３に示すように、線形パラメータおよびベルト部材の横方向曲げ係数εや係数Ｃｑ等の非線形パラメータからなる力学要素パラメータが設定され、スリップ角α、制駆動方向のスリップ率Ｓ、転動速度Ｖ、及び前後力Ｆx、横力Ｆｙ、トルクＭｚを入力することで、図３中の式（１）〜（８）によって処理された前後力、横力およびトルクの値（以降、前後力Ｆx’、横力Ｆｙ’、トルクＭｚ’とする）が算出されるように構成されている。勿論、入力された前後力Ｆx、横力ＦｙおよびトルクＭｚの値と、処理された前後力Ｆx’、横力Ｆｙ’およびトルクＭｚ’の値との誤差が所定値以下、すなわち略一致した（収束した、タイヤ力学モデルで力が釣り合い状態となった）場合にのみ、前後力Ｆx’、横力Ｆｙ’およびトルクＭｚ’の値が力の釣り合い状態を実現するタイヤの横力およびトルクの値として決定される。
なお、線形パラメータとは、式（６）〜（８）において線形の形式で表されている力学要素パラメータをいい、非線形パラメータとは、式（６）〜（８）において非線形の形式で表されている力学要素パラメータをいう。
【００２７】
タイヤ力学要素演算プログラム１２は、式（１）に示すように、入力されたトルクＭｚとねじりコンプライアンス（１／Ｇｍｚ）とによって求められるねじり戻し角を算出し、このねじり戻し角を、付与されたスリップ角αから差し引くことにより、実効スリップ角αｅを算出する。このように実効スリップ角αｅを算出するのは、トルクＭｚが０より大きい場合、トルクは付与されたスリップ角を低減するようにタイヤ自身に作用し捩じり戻す作用を有するからである。したがって、トルクＭｚが０より大きい場合、図４（ａ）に示すように、実際に付与されたスリップ角αに対して実効スリップ角αｅは小さくなる。
【００２８】
さらに、式（２）により、前後力Ｆxから接地圧分布の形状を規定する偏向係数ｑを算出する。偏向係数ｑとは、スリップ角α＝０の直進状態の接地圧分布（図５（ａ）参照）が、図５（ｂ）に示すように横力Ｆyが発生して接地圧分布が進行方向前方（接地面における踏込み端）に向かって偏向した接地圧分布の形状を表すパラメータである。この接地圧分布をｐ（ｔ）（ｔは、図５（ａ），（ｂ）中の進行方向の後方向に向かってｔ軸をとった場合の接地長さで規格化した座標位置）とすると、接地圧分布ｐ（ｔ）の形状は、図５（ｂ）中の式（９）で表される関数Ｄｇｓｐ(ｔ；ｎ，ｑ）で規定される。
ここで、関数Ｄｇｓｐ(ｔ；ｎ，ｑ）中の係数ｎは横力発生中の接地面の接地圧分布を規定するもので、図５（ｃ）に示すように接地圧分布の踏込み端および蹴りだし端付近で角張る（曲率が大きくなる）ように接地圧分布を規定する係数である。また、図５（ｄ）に示すように係数ｑが０から１になるにしたがって接地圧分布のピーク位置は踏込み端側に移動するように設定されている。このように係数qおよび係数ｎは、接地圧分布の形状を規定する形状規定係数である。
【００２９】
さらに、式（３）により、横力Ｆｙの発生時のタイヤ中心位置が踏込み端側へ移動する程度を表す値（xc／ｌ）をトルクＭｚと関連づけて算出する。ここでｌは、接地長である。このように式（３）においてタイヤ中心位置Ｏの移動を定めるのは、図５（ｂ）に示すように、トルクＭｚの回転中心となるタイヤ中心位置Ｏが横力Ｆｙの発生により接地面の踏込み側に移動するためである。
【００３０】
さらに、式（４）により、スリップ角αが大きいときに起こる接地面内での滑り摩擦と凝着摩擦との境界位置（ｌh／ｌ)を算出する。境界位置（ｌh／ｌ)は、以下のように定義される。
図６（ａ）〜（ｃ）に示される最大摩擦曲線は、凝着摩擦係数μｓに接地圧分布ｐ（ｔ）を乗算したものである。踏込み端で路面と接地したタイヤのトレッド部材は、蹴りだし端に移動するにつれてスリップ角αによって徐々に路面から横方向の剪断を受け、トレッド部材に横方向剪断力（凝着摩擦力）が発生する。さらに、路面の移動速度とタイヤの転動速度と差によって生じる制駆動方向のスリップ率Ｓによって、トレッド部材は徐々に路面から前後方向に剪断を受け、トレッド部材に前後方向剪断力（凝着摩擦力）が発生する。タイヤと路面との間に発生する剪断力は、横方向剪断力と前後方向剪断力との合力により表される。
この剪断力の合力は、徐々に大きくなって最大摩擦曲線に達すると、路面に凝着していたタイヤトレッド部材は滑り出し、滑り摩擦係数μｄに接地圧分布ｐ（ｔ）を乗算したすべり摩擦曲線に従って滑り摩擦力が発生する。図６（ａ）では、境界位置（ｌh／ｌ）より踏込み端側の領域がタイヤトレッド部材が路面に凝着した凝着域となり、蹴りだし側の領域がタイヤトレッド部材が路面に対して滑る滑り域となる。境界位置（ｌh／ｌ）は、式（４）により定まる。
図６（ｂ）は、スリップ角αが図６（ａ）に示すスリップ角αよりも大きくなった状態を示している。境界位置（ｌh／ｌ)は図６（ａ）に比べて踏込み端側に移動している。さらに、スリップ角αが大きくなると、図６（ｃ）に示すように接地面の踏込み端の位置からすべり摩擦が発生する状態となる。
【００３１】
図６（ａ）〜（ｃ）からわかるように、スリップ角αによって凝着域とすべり域の割合が大きく変化する。このような凝着域およびすべり域の横方向摩擦力、すなわち横力成分をタイヤ幅方向に沿って積分することによって横力Ｆｙ’を算出することができ、さらにタイヤ中心Ｏ周りのモーメントを算出することによってトルクＭｚ’を算出することができる。
同様に、前後方向剪断力についても、凝着域および滑り域の前後方向摩擦力、すなわち前後力成分をタイヤ幅方向に沿って積分することによって前後力Ｆx’を算出することができる。
式（７）〜（９）では、上述の凝着域および滑り域に分けて、実効スリップ角αｅを用いて前後力Ｆx’、横力Ｆｙ’およびトルクＭｚ’を算出する。
なお、滑り摩擦係数μdは、式（５）に示すように、滑り速度依存性を有するように規定されている。
【００３２】
式（６）は、２つの項（２つの前後力成分）の和によって前後力Ｆx’を算出する。第１項は積分範囲が０〜（ｌh／ｌ）の積分であって、凝着域に発生する凝着前後力成分を表す。第２項は積分範囲が（ｌh／ｌ）〜１の積分であって滑り域に発生する滑り前後力成分を表す。
式（７）では２つの項（２つの横力成分）の和によって横力Ｆｙ’を算出する。第１項は積分範囲が０〜（ｌh／ｌ）の積分であって、凝着域に発生する凝着横力成分を表す。第２項は積分範囲が（ｌh／ｌ）〜１の積分であってすべり域に発生する滑り横力成分を表す。
【００３３】
また、式（８）中、第１項は積分範囲が０〜（ｌh／ｌ）の積分であって、凝着域に発生する凝着横力成分によって生じるトルク成分を表し、第２項は積分範囲が（ｌh／ｌ）〜１の積分であってすべり域に発生するすべり横力成分によって生じるトルク成分を表す。なお、式（８）中では、上記２つのトルク成分の他に別のトルク成分、すなわち、第３、４項が設けられている。第３、４項は（ｌh／ｌ）・tanαｅに比例した項であり、後述するようにタイヤの接地面がスリップ角αによって横方向に移動し、この時の移動量とタイヤの前後力とによって生じるタイヤ中心Ｏ周りのトルク成分を表す。すなわち、トルクＭｚ’は、凝着横力によって生じるトルク成分、すべり横力によって生じるトルク成分および前後力によって生じるトルク成分の３つの成分の合計によって算出される。
【００３４】
式（７）中の第１項の凝着横力成分は凝着域における横力であり、式（７）では、実効スリップ角αｅによって生じるトレッド部材の横方向変位がベルトの横曲げ変形によって緩和された状態を表すことによって凝着横力成分を算出する。
第２項のすべり横力成分はすべり域における横力であり、式（７）では、実効スリップ角αｅによって生じる接地圧分布ｐ（ｔ）の形状を関数Ｄｇｓｐ(ｔ；ｎ，ｑ)で表してすべり横力成分を算出する。
なお、本実施形態では、横剛性及び前後剛性を同じ値Ｋtとしたタイヤ力学モデルであるが、本発明では、トレッド部材に設けられるトレッドパターンによって横剛性及び前後剛性は異なるものとして、異なるパラメータとして設定してもよい。この場合、横剛性をＫtとし、前後剛性をＫxとした場合、式（６）中の第１項に含まれるＫtをＫxとし、式（８）中の第３項に含まれるＫtをＫxとして、別々の非線形パラメータとして定めるとよい。
又、滑り域における滑り摩擦係数μdについても、トレッドパターンによって横方向と前後方向で異なるパラメータとして扱ってもよい。
【００３５】
図４（ａ）〜（ｃ）は、実効スリップ角αｅ、ベルトの変形によって生じる緩和された凝着横力成分および前後力成分とトルク成分との関係を、接地面を模式的に表して図示したものである。
図４（ａ）は、スリップ角αが付与された際、スリップ角αによって生じるトルクによってスリップ角αを減ずるようにタイヤ自身に作用し、実効スリップ角αｅとなっている状態を示している。図４（ｂ）は、この実効スリップ角αｅによって生じる横方向変位とベルトの横曲げ変形によって生じる横方向変位の関係を示している。図４（ｃ）はタイヤの接地面が横力によって横方向に移動することによって生じる、前後力分布がトルクＭｚ'に寄与するメカニズムを示している。図４(ｃ)中、Ｍｚ１およびＭｚ２は凝着横力成分によるトルク成分および滑り横力成分によるトルク成分を、Ｍｚ３は接地面に作用する前後力によるトルク成分を示している。
【００３６】
図７は、制駆動方向のスリップ率Ｓ、スリップ角αが付与されタイヤ力学モデルに基づいて前後力Ｆx'、横力Ｆｙ'およびトルクＭｚ'が算出されるまでの処理ブロック図である。図７からわかるように、本発明におけるタイヤ力学モデルは、前後力Ｆx'、横力Ｆｙ'およびトルクＭｚ'の算出の際、ベルトの横曲げ変形、接地圧分布の形状変化およびタイヤの捩じり変形がフィードバックされて式（６）〜（８）において算出される。なお、前後力Ｆx'、横力Ｆｙ'およびトルクＭｚ'を算出する際に用いるベルトの横曲げ変形、接地圧分布の形状変化およびタイヤの捩じり変形には、付与される横力ＦｙおよびトルクＭｚが用いられる。
【００３７】
なお、タイヤ力学モデル演算部１２において算出される前後力Ｆx'、横力Ｆｙ'およびトルクＭｚ'は、付与された前後力Ｆx、横力ＦｙおよびトルクＭｚと必ずしも一致しない。しかし、後述するタイヤ力学要素パラメータ導出プログラム１４、タイヤ軸力・トルクデータ算出プログラム１６において行なわれるシーケンス処理により、タイヤ力学モデルにおいて、付与される前後力Ｆx、横力ＦｙおよびトルクＭｚと算出される前後力Ｆx'、横力Ｆｙ'およびトルクＭｚ'とが略一致する（力の釣り合い状態となる）ように、前後力Ｆx、横力ＦｙおよびトルクＭｚが探索され、タイヤ力学モデルにおいて釣り合い状態にある前後力、横力およびトルクが算出される。
【００３８】
タイヤ力学要素パラメータ導出プログラム１４は、制駆動方向のスリップ率Ｓを例えば０〜１まで変化させ、さらにスリップ角αを例えば０〜２０度まで変化させ、接地面に凝着域と滑り域が存在する時の前後力Ｆx’、横力Ｆｙ’およびトルクＭｚ’が、与えられた特性曲線の前後力Ｆx、横力ＦｙおよびトルクＭｚに一致するように上述の線形パラメータおよび非線形パラメータを導出する部分である。
【００３９】
図８は、タイヤ力学要素パラメータ導出プログラム１４において行なわれる処理の流れを示している。
具体的には、図８に示すように、一定の負荷荷重において制駆動方向のスリップ率を種々変えた条件の下、スリップ角を０〜２０度変化させて前後力Ｆx、横力ＦｙおよびトルクＭｚの特性曲線を取得する（ステップＳ２０２）。
次に、非線形パラメータである横方向曲げ係数ε、ねじりコンプライアンス（１／Ｇｍｚ）、係数ｎ、前後剛性・横剛性と凝着摩擦係数との比Ｋt／μｓ、の比係数Ｃｑ、移動係数Ｃｘｃ等を所定の値に初期設定する（ステップＳ２０４）。
【００４０】
次に、前後力Ｆx、横力ＦｙおよびトルクＭｚの特性曲線と初期設定された非線形パラメータを用いて線形最小二乗回帰を行なう（ステップＳ２０６）。具体的には、滑り速度０における滑り摩擦係数μd0、係数ｂv、前後剛性・横剛性Ｋtの線形パラメータに関する正規方程式を作成し、この正規方程式を解くことによって線形パラメータの値を算出する。すなわち、線形最小二乗回帰を行なう。ここで、正規方程式とは残差二乗和を上記線形パラメータのそれぞれで偏微分して偏微分値を０とした線形パラメータの個数分作成される線形パラメータに関する方程式である。
こうして初期設定された非線形パラメータの値および正規方程式を用いて算出された線形パラメータの値および前後力Ｆx、横力ＦｙおよびトルクＭｚの特性曲線のデータをタイヤ力学モデル演算プログラム１２に付与する。この付与によって図７のブロック図の流れに従って制駆動方向のスリップ率Ｓにおけるスリップ角αにおける前後力Ｆx’、横力Ｆｙ’およびトルクＭｚ’が算出される。
【００４１】
次に、これらの前後力Ｆx’、横力Ｆｙ’およびトルクＭｚ’の算出データと、付与されたスリップ角αにおける前後力Ｆx、横力ＦｙおよびトルクＭｚの特性曲線のデータとを用いて、下記式（１０）で表される複合二乗残差和Ｑｃを計算する（ステップＳ２０８）。この場合、式（１０）中のＮは付与される制駆動方向のスリップ率Ｓ、スリップ角αの条件設定数である。また、このときの重み付け係数ｇx,ｇy，ｇｍは、Ｎ個のスリップ角、制駆動方向のスリップ率の条件における前後力Ｆx、横力Ｆｙ，トルクＭｚの分散から求められたものである。
【００４２】
【数１】
2006232011.tif 000003
【００４３】
また、ｇx,ｇyおよびｇｍは、前後力Ｆx，横力ＦｙおよびトルクＭｚのデータについてＮ個のスリップ角、スリップ率の条件における前後力Ｆx,横力ＦｙおよびトルクＭｚの分散をσx２,σy２およびσｍ２としたとき、下記式で表される係数であり、複合二乗残差和Ｑｃを求める際に用いる重み付け係数である。
ｇx＝１／σx２
ｇy＝１／σy２ 
ｇｍ＝１／σｍ２
すなわち、複合二乗残差和Ｑｃは特性曲線の値のばらつきの情報である分散の逆数を重み付け係数とし、前後力、横力およびトルクのそれぞれの二乗残差和を重み付け加算したたものである。
【００４４】
このように、前後力Ｆxの特性曲線の値とタイヤ力学モデル演算プログラム１２にて算出された前後力Ｆx’の値との残差二乗和に重み付け係数ｇxを乗算したものと、横力Ｆｙの特性曲線の値とタイヤ力学モデル演算プログラム１２にて算出された横力Ｆｙ’の値との残差二乗和に重み付け係数ｇyを乗算したものと、トルクＭｚの特性曲線の値とタイヤ力学モデル演算プログラム１２にて算出されたトルクＭｚ’の値との残差二乗和に重み付け係数ｇｍを乗算したものとを加算することによって複合二乗残差和を計算する。ここで、複合二乗残差和を用いるのは、上述の場合と同様に、非線形パラメータの算出において複数のスリップ角、スリップ率の条件における前後力Ｆx’、横力Ｆｙ’およびトルクＭｚ’を前後力Ｆx、横力ＦｙおよびトルクＭｚのそれぞれに同時に一致させるためである。
複合二乗残差和は、所定値以下となって収束しているか否かを判別する（ステップＳ２１０）。
収束していないと判別すると、ステップＳ２０４で初期設定された非線形パラメータの調整を行う（ステップＳ２１２）。この非線形パラメータの調整は、例えばNewton-Raphson法に従って行なわれる。
【００４５】
この非線形パラメータはステップＳ２１０において収束すると判別されるまで調整されるが、この調整が行なわれる度に、線形パラメータに関する線形最小二乗回帰（ステップＳ２０６）および複合二乗残差和の計算（ステップＳ２０８）を行なって、上記式（１０）による複合残差二乗和を求める。そして、複合二乗残差和が所定値以下になるまで、非線形パラメータの調整が行なわれる。複合二乗残差和が所定値以下になると、線形最小二乗回帰で算出された各非線形パラメータを決定し（ステップＳ２１４）、これらの力学要素パラメータをメモリ４に記憶する。
以上が、タイヤ力学要素パラメータ導出プログラム１４の行なう、タイヤ力学モデルを用いた制駆動方向のスリップ率Ｓ，スリップ角αにおける線形パラメータおよび非線形パラメータの導出の流れである。
【００４６】
タイヤ軸力・トルクデータ算出プログラム１６は、所定の負荷荷重における、上記タイヤ力学モモデルにおけるタイヤ力学要素パラメータである線形パラメータおよび非線形パラメータを用いて、所定の制駆動方向のスリップ率及び所定のスリップ角におけるタイヤ軸力・トルクデータを算出する部分である。例えば、制駆動方向のスリップ率Ｓを一定として、スリップ角αを０〜２０度まで変化させたときの特性曲線を求めたり、スリップ角αを一定として、制駆動方向のスリップ率Ｓを０〜１まで変化させたときの特性曲線を求めたり、また、スリップ角αを一定にして、縦軸を横力、横軸を前後力としたときの摩擦楕円を求めたりする。
【００４７】
図９は、タイヤ軸力・トルクデータ算出プログラム１６において行なわれる処理の流れの一例を示している。この例は、前後力、横力およびトルクのスリップ角αの依存性を表す特性曲線を算出する例である。
タイヤ軸力・トルクデータ算出プログラム１６は、まず、導出された線形パラメータおよび非線形パラメータをメモリ４から読み出して設定する（ステップＳ４００）。
さらに、負荷荷重Ｆzにおける前後力Ｆx、横力ＦｙおよびトルクＭｚを初期設定する（ステップＳ４０２）。
この後、スリップ角依存性を表す特性曲線を算出する場合、設定されたスリップ角α＝Δαとともに線形パラメータおよび非線形パラメータおよび初期設定された前後力Ｆx、横力ＦｙおよびトルクＭｚをタイヤ力学モデル演算プログラム１２に付与する。タイヤ力学モデル演算プログラム１２では、付与された線形パラメータおよび非線形パラメータと、初期設定された前後力Ｆx、横力ＦｙおよびトルクＭｚが用いられて式（６）〜（８）に従って前後力Ｆx’、横力Ｆｙ'、トルクＭｚ’が算出される（ステップＳ４０４）。
【００４８】
こうして算出された前後力Ｆx’、横力Ｆｙ'、トルクＭｚ’はタイヤ軸力・トルクデータ算出プログラム１６に返される。タイヤ軸力・トルクデータ算出プログラム１６は、タイヤ力学モデル演算プログラム１２に付与した前後力Ｆx、横力ＦｙおよびトルクＭｚの設定値と算出された前後力Ｆx’、横力Ｆｙ'、トルクＭｚ’の計算値との複合二乗残差和を式（１０）に従って算出する（ステップＳ４０６）。
次に、算出された複合二乗残差和が所定値以下となって収束しているか否かを判別する（ステップＳ４０８）。
収束していないと判別すると、先に設定された前後力Ｆx、横力ＦｙおよびトルクＭｚの設定値を調整する（ステップＳ４１０）。この調整された前後力Ｆx、横力ＦｙおよびトルクＭｚと線形パラメータおよび非線形パラメータとが再度タイヤ力学モデル演算プログラム１２に付与される。
こうして、複合二乗残差和が所定値以下となって収束するまで、前後力Ｆx、横力ＦｙおよびトルクＭｚの設定値を調整する。この設定値の調整は、例えば上述したNewton-Raphson法に従って行なわれる。こうして、前後力Ｆx’、横力Ｆｙ'、トルクＭｚ'を決定する（ステップＳ４１２）。
【００４９】
次に、スリップ角αが所定のスリップ角以下であるか否かを判別する（ステップＳ４１６）。
スリップ角αが所定のスリップ角以下であると判別した場合、スリップ角αの条件が変更される（α→α＋Δα）（ステップＳ４１４）。そして、変更されたスリップ角αにおける前後力Ｆx、横力Ｆｙ、トルクＭｚの初期値が設定され（ステップＳ４０２）、前後力Ｆx’、横力Ｆｙ'およびトルクＭｚ’が算出され（ステップＳ４０４）、複合二乗残差和が算出され（ステップＳ４０６）、この複合二乗残差和の収束が判別される（ステップＳ４０８）。
こうして、スリップ角αが所定スリップ角となるまで繰り返し変更される（ステップＳ４１６）。このスリップ角の変更の度に前後力Ｆx’、横力Ｆｙ’およびトルクＭｚ’を算出し、収束する前後力Ｆx’、横力Ｆｙ’およびトルクＭｚ’を決定する。決定された前後力Ｆx’、横力Ｆｙ’およびトルクＭｚ’はメモリ４に記憶される。
このようにして、スリップ角αに依存する前後力、横力およびトルクの特性曲線を求める。
【００５０】
滑り量算出プログラム１８は、タイヤ軸力・トルクデータ算出プログラム１６において、収束した前後力、横力およびトルクが算出されるたびに、このときの式（４）中のｌ・[Ｓ2＋｛tan(αe)  −(ε・l) ・Ｆy・（１−ｌh／ｌ）｝2](1/2)を算出し、この算出結果をタイヤの滑り量とする部分である。
図１０に示すように、凝着域と滑り域との境界位置（ｌh／ｌ）における変形量であるｌ・[Ｓ2＋｛tan(αe)  −(ε・l) ・Ｆy・（１−ｌh／ｌ）｝2](1/2)がタイヤ滑り量となる。算出されたタイヤ滑り量はメモリ４に記憶される。
【００５１】
摩耗特性予測プログラム２０は、滑り量算出プログラム１８で算出され、メモリ４に記憶されたタイヤの滑り量を、トレッドゴム部材の摩耗特性データと共に用いて、所定の横力が発生時のタイヤのトレッド部材の摩耗特性を予測する部分である。
トレッドゴム部材の摩耗特性は、予め新ＪＩＳ  Ｋ６２６４で規定される試験法にて求められるゴム部材の物性データをデータベースとして保有しており、この物性データと所定の横力発生時のタイヤの滑り量とを乗算した結果をタイヤの摩耗特性として予測する部分である。
ゴム部材の物性データとして、例えばリング状試験片をＤＩＮ摩耗試験法によって摩耗量を測定して得られる単位走行距離当たりの摩耗量を求める。このゴム物性としての摩耗量は、予め基準ゴム部材を定め、この基準ゴム部材の摩耗量を基準として指数として表すこともできる。
ゴム部材のゴム物性データとして、この他に破断伸びＥbや破断強度Ｔbや破断エネルギーＥb×Ｔbを用いることもできる。この場合、基準ゴム部材の値を基準値として指数を用いることが好ましい。
このようにして算出された摩耗特性の予測結果としての数値或いは指数が、出力装置７に画面表示、或いはプリント出力される。
装置１は、以上のように構成される。
【００５２】
装置１は、タイヤに制駆動方向のスリップ率が与えられて接地面に滑り域を形成するときのタイヤ滑り量を求めることによりタイヤの摩耗特性を予測するタイヤの摩耗予測を行わせるための装置として機能するが、このような装置は、コンピュータに所定の機能を動作させる下記プログラムを実行することによって構成される。
このようなプログラムは、タイヤに制駆動方向のスリップ率を与えたときのタイヤ回転軸に作用するタイヤ軸力及びセルフアライニングトルクの、制駆動方向のスリップ率に対して変化する特性曲線を取得し、コンピュータのメモリ４に記憶させる手順と、
メモリ４に記憶された特性曲線から、複数のタイヤ力学要素パラメータを用いて構成されるタイヤ力学モデルに基づいて、この特性曲線を定めるタイヤ力学要素パラメータの値を、コンピュータの演算ユニットに導出させる手順と、
導出されたタイヤ力学要素パラメータの値をタイヤ力学モデルに与えて、タイヤと接地面の間に形成される凝着域と滑り域を前記演算手段に求めさせ、この滑り域からタイヤ滑り量をコンピュータの演算ユニットに求めさせる手順と、
このタイヤ滑り量を、トレッドゴム部材の摩耗特性データと共に用いて、制駆動方向のスリップ率が与えられたときのタイヤのトレッド摩耗  特性を、コンピュータの演算ユニットに評価させる手順と、を有するとよい。
【００５３】
図１１は、このような装置１を用いて行われる本発明のタイヤの摩耗特性予測方法を示すフローチャートである。
まず、タイヤ力学要素パラメータ導出プログラム８では、メモリ４に記憶されている特性曲線の計測データが呼び出されて、制駆動方向のスリップ率を種種変えた条件の下、スリップ角を０〜２０度に変化させたときの特性曲線が取得される（ステップＳ３００）。特性曲線は、例えばFLAT TRAC-I,FLAT TRAC-II,FLAT TRAC-III（ＭＴＳ社製商品名）等の室内試験機を用いて計測データが取得される。
【００５４】
次に、タイヤ力学モデル演算プログラム１２が実行され、取得した特性曲線の値を用いて、タイヤ力学モデルで定義されている線形パラメータ及び非線形パラメータを含むタイヤ力学要素パラメータの値が導出される（ステップＳ３０２）。
次に、導出されたタイヤ力学要素パラメータを用いて、タイヤ軸力・トルクデータ算出プログラム１６において、所定の制駆動方向のスリップ率Ｓにおける特性曲線が算出されるとともに、所定の横力発生時のタイヤ滑り量が算出される（ステップＳ３０４）。タイヤ滑り量は、境界位置（ｌh／ｌ）  における変形量であるｌ・[Ｓ2＋｛tan(αe)  −(ε・l) ・Ｆy・（１−ｌh／ｌ）｝2](1/2)である。
【００５５】
次に、算出されたタイヤ滑り量と、タイヤの接地幅と、メモリ４に記憶されているトレッド部材のゴム物性データをデータベースから呼び出して、ゴム物性データとを乗算して、摩耗特性が予測される（ステップＳ３０６）。
ゴム物性データは、例えばゴム部材の室内摩耗試験結果のデータであり、又破断エネルギー（Ｅb×Ｔb）が所定の基準ゴム部材を基準とした指数で表したものである。
こうして、タイヤ滑り量、ゴム物性データおよび接地幅を乗算した結果が出力装置７に出力される。
【００５６】
このようなタイヤの摩耗予測は、タイヤ設計に用いることができる。
図１２は、このようなタイヤ設計の流れを示すフローチャートである。
まず、タイヤ力学要素パラメータの値が設定される（ステップＳ３２０）。
この設定されたタイヤ力学要素パラメータの値を用いてタイヤ軸力・トルクデータ算出プログラム１６において、所定の制駆動方向のスリップ率Ｓ、スリップ角における、タイヤ軸力、トルクの値が算出されるとともに、所定のタイヤ軸力（横力）発生時のタイヤ滑り量が算出される（ステップＳ３２２）。このタイヤ滑り量を用いて、上述したステップＳ３０６のように摩耗特性の予測が行われる（ステップＳ３２４）。予測された摩耗特性が予め定められた目標を達成するか、否かが判別される（ステップＳ３２６）。判別の結果、目標を達成しない場合、タイヤ力学要素パラメータが修正される（ステップＳ３２８）。修正されたタイヤ力学要素パラメータは、ステップＳ３２２に戻され、再度タイヤ滑り量が算出される。こうして、ステップＳ３２６において、目標を達成するまで繰り返しタイヤ力学要素パラメータが修正される。修正の方法は特に限定されないが、例えば予め定められた変化幅を用いて段階的に修正する。
最後に、目標を達成するタイヤ力学要素パラメータの値は、タイヤ設計仕様特性として決定され（ステップＳ３３０）、この特性に基づいて、タイヤ形状を含めたタイヤの構造設計及びゴム部材の配合設計を含めた材料設計が行われる。
こうして設計されたタイヤは、タイヤ設計仕様特性を満足するように製造される。
【００５７】
本発明のタイヤ摩耗予測方法について、タイヤに制駆動方向のスリップ率、スリップ角を与えたときの特性曲線から、所定の制駆動方向のスリップ率、スリップ角におけるタイヤの摩耗特性を予測する方法を説明したが、本発明においては、制駆動方向のスリップ率を０として、スリップ角依存性の特性曲線からコーナリング時のタイヤの摩耗特性を予測することもできる。
図１３は、このときのタイヤ力学モデルを説明する図である。
図１３に示すタイヤ力学モデルを表す式（１１）〜（１６）は、図３に示す力学モデルのうち、式（１）〜（８）に対して制駆動方向のスリップ率Ｓ＝０とすることにより得られる式と対比して、式（２），（３）中のＦxをＭzとしたことと、式（８）中の第３、第４項が異なること以外は同じである。
ここで、図１３中の式（１６）におけるＡxは、接地面内の前後剛性Ａｘであり、接地面が横方向に移動することによりトルクＭｚ'に影響を与える線形パラメータである。
また、図１４は、スリップ角αが付与されタイヤ力学モデルに基づいて横力Ｆｙ'およびトルクＭｚ'が算出されるまでの処理ブロック図である。図１４からわかるように、タイヤ力学モデルは、横力Ｆｙ'およびトルクＭｚ'の算出の際、ベルトの横曲げ変形、接地圧分布の形状変化およびタイヤの捩じり変形がフィードバックされて式（１５），（１６）において算出される。なお、横力Ｆｙ'およびトルクＭｚ'を算出する際に用いるベルトの横曲げ変形、接地圧分布の形状変化およびタイヤの捩じり変形には、付与される横力ＦｙおよびトルクＭｚが用いられる。
【００５８】
このようなタイヤ力学モデルを用いてタイヤ力学要素パラメータを導出するタイヤ力学要素パラメータ導出プログラム１４では、横力及びトルクの特性曲線が与えられると、この２つの特性曲線に合うようにタイヤ力学要素パラメータの値が導出される。
この場合、複合二乗残差和Ｑｃは、下記式（１７）に従って算出される。この場合、式（１７）中のＮは付与されるスリップ角αの条件設定数である。また、このときの重み付け係数ｇy，ｇｍは、Ｎ個のスリップ角の条件における横力Ｆｙ，トルクＭｚの分散から求められたものである。
【００５９】
【数２】
2006232011.tif 000004
【００６０】
さらに、本発明のタイヤ摩耗予測方法においては、スリップ角を０として、制駆動方向のスリップ率依存性の特性曲線から、タイヤ制動時、又は駆動時におけるタイヤの摩耗特性を予測することもできる。
図１５は、このときのタイヤ力学モデルを説明する図である。
図１５に示すタイヤ力学モデルを表す式（２１）〜（２４）は、図３に示す力学モデルのうち、式（１）〜（８）に対してスリップ角α＝０とすることにより得られる式と同じである。
図１６（ａ）は、制駆動方向のスリップ率Ｓ、転動速度Ｖが付与されタイヤ力学モデルに基づいて前後力Ｆx'が算出されるまでの処理ブロック図である。図１６（ａ）からわかるように、タイヤ力学モデルは、前後力Ｆx'の算出の際、接地圧分布の形状変化がフィードバックされて式（２４）において算出される。
図１６（ｂ）は、このときの接地面内における凝着域及び滑り域を説明する図である。
図１６（ｂ）では、踏み込み端から５ｃｍの位置がｌhに対応し、この位置におけるトレッド部材の前後方向の、路面に対する変位量が滑り量とされる。
【００６１】
このようなタイヤ力学モデルを用いてタイヤ力学要素パラメータを導出するタイヤ力学要素パラメータ導出プログラム１４では、前後力の制駆動方向のスリップ率依存性の特性曲線が与えられると、この特性曲線に合うようにタイヤ力学要素パラメータの値が導出される。
この場合、収束しているか否かの判別に用いる複合二乗残差和Ｑｃは、特性曲線における前後力Ｆxの値と、タイヤ力学モデルを用いて算出された前後力Ｆx’の値との二乗残差和である。
【００６２】
図１７（ａ）〜（ｃ）は、図１３及び図１４に示すタイヤ力学モデルを用いて、横力のスリップ角依存性の特性曲線からタイヤ力学要素パラメータの値を導出したときのタイヤ滑り量を表す図である。
図１７（ａ）は、与えられた横力の特性曲線を表した図であり、図１７（ｂ）は、この特性曲線から導出されたタイヤ力学要素パラメータの値を用いて算出された横力の算出データを示した図である。図１７（ｃ）は、導出されたタイヤ力学要素パラメータの値を用いて算出されるタイヤ滑り量を示している。
図１７（ｂ）に示すように、横力の算出データは、図１７（ａ）に示す特性曲線を忠実に再現しており、タイヤ力学要素パラメータの値が精度高く導出されたことを意味し、また、タイヤ力学モデルが特性曲線を忠実に再現できるモデルであることがわかる。
【００６３】
図１７（ｃ）は、このようにして求められたタイヤ力学要素パラメータの値を用いて算出されるタイヤ滑り量の横力依存性を表した図である。タイヤＡ〜Ｅの５種類である。これからわかるように、タイヤ滑り量は、タイヤによって種々異なっていることがわかる。
このようなタイヤ滑り量のうち、所定の横力におけるタイヤ滑り量を代表値として、上述したように、トレッド部材のゴム物性データとタイヤの接地幅とを乗算して、タイヤの摩耗特性を予測する。具体的には、ゴム物性データとして所定の摩耗試験法によって摩耗量を測定して得られる単位走行距離当たりの摩耗量の場合、摩耗特性の値が大きいほど、摩耗特性は劣っているものと判断される。
【００６４】
図１８（ａ）、（ｂ）は、図１５及び図１６に示すタイヤ力学モデルを用いて、前後力の制駆動方向のスリップ率依存性の特性曲線に従ってタイヤ力学要素パラメータの値を導出したときの特性曲線と算出データを示す図である。
図１８（ａ）は、与えられた前後力の特性曲線を表した図であり、図１８（ｂ）は、この特性曲線から導出されたタイヤ力学要素パラメータの値を用いて算出された前後力の算出データを示した図である。
図１８（ｂ）に示すように、前後力の算出データは、図１８（ａ）に示す特性曲線を忠実に再現しており、タイヤ力学要素パラメータの値が精度高く導出されたことを意味し、また、タイヤ力学モデルが特性曲線を忠実に再現できるモデルであることがわかる。
【００６５】
さらに、タイヤ力学モデルを構成するタイヤ力学要素パラメータを一旦導出すると、このパラメータを用いて、異なる制駆動方向のスリップ率及びスリップ角における横力及び前後力を算出することができる。図１９はこのようにして算出される前後力と横力の関係の一例（コーナリング制動時の摩擦円）を示している。
タイヤ滑り量は、このようなタイヤ力学モデルと、このモデルを構成するタイヤ力学要素パラメータを用いて算出される。
【００６６】
図２０（ａ）は、タイヤＡ〜Ｃにおけるトレッドゴム部材のゴム物性データとして、破断伸びＥbと破断強度Ｔbのデータを示している。指数が低いほうが破断伸び及び破談強度が小さく、破断エネルギーＥb×Ｔbが大きいほどゴム部材としての耐摩耗性がよいという点からすると、ゴム物性から判断されるタイヤの摩耗特性はタイヤＣ（優）＞タイヤＢ＞タイヤＡ（劣）と予測される。しかし、実際のタイヤでの摩耗寿命は、タイヤＡ（６００００ｋｍ以上）＞タイヤＢ（４０３００ｋｍ）＞タイヤＣ（２８３００ｋｍ）となっており、ゴム物性データでは説明ができない。
【００６７】
しかし、本発明に従って、図２０（ｄ）に示すようにタイヤ滑り量を算出することにより、タイヤの摩耗特性を予測することができる。
図２０（ｂ），（ｃ）は、室内試験で計測される横力及びトルクのスリップ角依存性の特性曲線を示す図であり、導出されたタイヤ力学要素パラメータの値を用いて特性曲線の対応する横力及びトルクを算出して●、▲、■でプロットしている。横力、トルクの特性曲線のいずれにおいても、●、▲、■でプロットは特性曲線上に乗っており、精度よくタイヤ力学要素パラメータの値が導出されていることがわかる。
この導出されたパラメータの値を用いて算出されるタイヤ滑り量は、図２０（ｄ）に示すように、タイヤＡはタイヤＢ，Ｃに比べて極めて小さく、滑り量の小さい順番は、タイヤＡ、タイヤＢ、タイヤＣの順になっている。したがって、このタイヤ滑り量と上述した破断エネルギーＥb×Ｔbの指数とを乗算した結果は、タイヤＡ（摩耗特性：優）＞タイヤＢ＞タイヤＣ（摩耗特性：劣）となる。これがタイヤ摩耗特性の予測である。この予測結果は、ゴム部材のゴム物性データでは説明のできない上述の摩耗寿命の順を説明することができる。
【００６８】
このように、タイヤ力学モデルを構成するタイヤ力学要素パラメータの値を特性曲線から導出し、これを用いて路面に対するタイヤ滑り量を算出することで、タイヤの摩耗特性を精度よく評価することができる。またタイヤ力学要素パラメータの値を修正することにより、修正したパラメータの値を有するタイヤの摩耗特性を容易に予測することができる。
【００６９】
以上、本発明のタイヤの摩耗予測方法、タイヤの設計方法、タイヤの製造方法、タイヤの摩耗予測システム及びプログラムについて詳細に説明したが、本発明は上記実施形態に限定されず、本発明の主旨を逸脱しない範囲において、種々の改良や変更をしてもよいのはもちろんである。
【図面の簡単な説明】
【００７０】
【図１】本発明のタイヤの摩耗予測方法を実施する、本発明のタイヤの摩耗予測システムの一例の装置のブロック図である。
【図２】本発明のタイヤの摩耗予測方法において用いられるタイヤ力学モデルを説明する図である。
【図３】本発明のタイヤの摩耗予測方法において用いられるタイヤ力学モデルを説明する他の説明図である。
【図４】（ａ）〜（ｃ）は、本発明のタイヤの摩耗予測方法において用いられるタイヤ力学モデルを説明する他の説明図である。
【図５】（ａ）〜（ｄ）は、本発明のタイヤの摩耗予測方法において用いられるタイヤ力学モデルを説明する他の説明図である。
【図６】（ａ）〜（ｃ）は、本発明のタイヤの摩耗予測方法において用いられるタイヤ力学モデルを説明する他の説明図である。
【図７】本発明のタイヤの摩耗予測方法において用いられるタイヤ力学モデルにおいて前後力、横力およびトルクを算出するまでの処理ブロック図である。
【図８】本発明のタイヤの摩耗予測方法において実施される１つの処理の流れを示すフローチャートである。
【図９】本発明のタイヤの摩耗予測方法において実施される１つの処理の流れを示すフローチャートである。
【図１０】本発明のタイヤの摩耗予測方法において用いられるタイヤ力学モデルを説明する他の説明図である。
【図１１】本発明のタイヤの摩耗予測方法の処理の流れを示すフローチャートである。
【図１２】本発明のタイヤの設計方法の処理の流れを示すフローチャートである。
【図１３】本発明のタイヤの摩耗予測方法において用いられる他の実施形態として用いられるタイヤ力学モデルを説明する図である。
【図１４】図１３に示すタイヤ力学モデルにおいて横力およびトルクを算出するまでの処理ブロック図である。
【図１５】本発明のタイヤの摩耗予測方法において用いられる他の実施形態として用いられるタイヤ力学モデルを説明する図である。
【図１６】（ａ）は、図１５に示すタイヤ力学モデルにおいて前後力を算出するまでの処理ブロック図であり、（ｂ）は、図１５に示すタイヤ力学モデルにおける凝着域及び滑り域を説明する図である。
【図１７】（ａ）は、本発明のタイヤの摩耗予測方法において用いられる特性曲線の例を示す図であり、（ｂ）は本発明のタイヤの摩耗予測方法において算出される特性曲線に対応する算出データの例を示す図であり、（ｃ）は、本発明のタイヤの摩耗予測方法において算出されるタイヤの滑り量の例を示す図である。
【図１８】（ａ）は、本発明のタイヤの摩耗予測方法において用いられる特性曲線の他の例を示す図であり、（ｂ）は本発明のタイヤの摩耗予測方法において算出される特性曲線に対応する算出データの他の例を示す図である。
【図１９】本発明のタイヤの摩耗予測方法において算出される摩擦楕円の例を示す図である。
【図２０】（ａ）〜（ｄ）は、本発明のタイヤの摩耗予測方法において求められるタイヤの摩耗特性の例を説明する図である。
【符号の説明】
【００７１】
１  装置
２  ＣＰＵ
３  バス
４  メモリ
５  入力操作系
６  インターフェース
７  出力装置
８  プログラム群
１０  統合・管理プログラム
１２  タイヤ力学モデル演算プログラム
１４  タイヤ力学要素パラメータ導出プログラム
１６  タイヤ軸力・トルク算出プログラム
１８  滑り量算出プログラム
２０  摩耗特性予測プログラム
【図１】
2006232011.tif 000005
【図２】
2006232011.tif 000006
【図３】
2006232011.tif 000007
【図４】
2006232011.tif 000008
【図５】
2006232011.tif 000009
【図６】
2006232011.tif 000010
【図７】
2006232011.tif 000011
【図８】
2006232011.tif 000012
【図９】
2006232011.tif 000013
【図１０】
2006232011.tif 000014
【図１１】
2006232011.tif 000015
【図１２】
2006232011.tif 000016
【図１３】
2006232011.tif 000017
【図１４】
2006232011.tif 000018
【図１５】
2006232011.tif 000019
【図１６】
2006232011.tif 000020
【図１７】
2006232011.tif 000021
【図１８】
2006232011.tif 000022
【図１９】
2006232011.tif 000023
【図２０】
2006232011.tif 000024