(57)【特許請求の範囲】
【請求項１】
タイヤのコーナリング時の過渡応答の計測データからタイヤ力学モデルで用いる過渡応答パラメータの値を算出することを、コンピュータに実行させるデータ処理方法であって、
スリップ角の時系列データをタイヤ測定条件としてタイヤに与えた時に計測されるコーナリング時のタイヤ回転軸に作用する横力およびセルフアライニングトルクの物理量の過渡応答の計測データをコンピュータに取得させるステップと、
タイヤ回転軸に作用する前記物理量の、スリップ角に対する応答関数であって、コーナリング時のタイヤの変形応答を再現するタイヤ力学モデルを用いて定めた応答関数に、過渡応答パラメータの値を付与し、前記値が付与された前記応答関数と前記タイヤ測定条件として用いた前記スリップ角の時系列データとを用いて積分を行うことによって、前記物理量の時系列データをタイヤの予測データとしてコンピュータに算出させるステップと、
算出した前記予測データが許容範囲内で前記計測データに一致するまで、前記過渡応答パラメータの値の修正と、前記予測データの算出をコンピュータに繰り返し行わせるステップと、
前記予測データが前記計測データに前記許容範囲内で一致したときに付与した前記過渡応答パラメータの値を、前記タイヤの過渡応答を再現する過渡応答パラメータの値としてコンピュータに決定させるステップと、を有し、
前記タイヤ力学モデルを用いて定めた前記応答関数は、前記タイヤのベルト部材の曲げ変形の遅れ応答をむだ時間による遅れ応答として構成され、
前記応答関数は、下記式で表される横力Ｆｙ及びセルフアライニングトルクＭｚのスリップ角αに対する伝達関数Ｆｙ（ｓ）／α（ｓ）及びＭｚ（ｓ）／α（ｓ）を用いて得られる時間の関数を含み、
0005790018.tif 000028
0005790018.tif 000029
前記式中のＫｙ、Ａｓ、Ｋｔ、Ｔ０、Ｔ１、Ｔ２、及びＴ３は、前記過渡応答パラメータであり、Ｔ１は前記タイヤのベルト部材の曲げ変形の遅れ応答のむだ時間であり、ｓは、応答関数を定める複素数である、ことを特徴とするデータ処理方法。
【請求項２】
前記応答関数はステップ応答関数であって、
前記予測データの算出は、前記値が付与された前記応答関数と前記スリップ角の時系列データの時間変化量との畳み込み積分をコンピュータに行わせることによって、行われる、請求項１に記載のデータ処理方法。
【請求項３】
前記タイヤの変形応答は、前記タイヤのトレッド部材の変形応答、前記タイヤのベルト部材の曲げ変形の遅れ応答、及び前記タイヤのサイドウォール部材のねじり変形の遅れ応答を含む、請求項１または２に記載のデータ処理方法。
【請求項４】
さらに、決定された前記過渡応答パラメータの値を用いて、前記タイヤのベルト部材の横曲げ剛性および前記タイヤのサイドウォール部材のねじり剛性をコンピュータに少なくとも算出させるステップを有する、請求項１〜３のいずれか１項に記載のデータ処理方法。
【請求項５】
スリップ角を時系列データとして与えたときのコーナリング時のタイヤの過渡応答データをタイヤ力学モデルを用いて算出することを、コンピュータに実行させるタイヤの過渡応答データの算出方法であって、
スリップ角の時系列データを取得するとともに、コーナリング時のタイヤの変形応答を再現するタイヤ力学モデルを用いて定めた応答関数であって、タイヤ回転軸に作用する横力およびセルフアライニングトルクの物理量の、スリップ角に対する応答関数に用いる過渡応答パラメータの値をコンピュータに取得させるステップと、
前記応答関数に前記過渡応答パラメータの値を付与し、前記値が付与された前記応答関数と取得した前記スリップ角の時系列データとを用いて積分を行うことによって、前記物理量の時系列データをタイヤの過渡応答データとしてコンピュータに算出させるステップと、を有し、
前記タイヤ力学モデルを用いて定めた前記応答関数は、前記タイヤのベルト部材の曲げ変形の遅れ応答をむだ時間による遅れ応答として構成され、
前記応答関数は、下記式で表される横力Ｆｙ及びセルフアライニングトルクＭｚのスリップ角αに対する伝達関数Ｆｙ（ｓ）／α（ｓ）及びＭｚ（ｓ）／α（ｓ）を用いて得られる時間の関数を含み、
0005790018.tif 000030
0005790018.tif 000031
前記式中のＫｙ、Ａｓ、Ｋｔ、Ｔ０、Ｔ１、Ｔ２、及びＴ３は、前記過渡応答パラメータであり、Ｔ１は前記タイヤのベルト部材の曲げ変形の遅れ応答のむだ時間であり、ｓは、応答関数を定める複素数である、ことを特徴とするタイヤの過渡応答データの算出方法。
【請求項６】
前記応答関数はステップ応答関数であって、
前記タイヤの過渡応答データの算出は、前記値が付与された前記応答関数と前記スリップ角の時系列データの時間変化量との畳み込み積分をコンピュータに行わせることによって、行われる、請求項５に記載のタイヤの過渡応答データの算出方法。
【請求項７】
前記タイヤの変形応答は、前記タイヤのトレッド部材の変形応答、前記タイヤのベルト部材の曲げ変形の遅れ応答、及び前記タイヤのサイドウォール部材のねじり変形の遅れ応答を含む、請求項５または６に記載のタイヤの過渡応答データの算出方法。
【請求項８】
タイヤのコーナリング時の過渡応答の計測データからタイヤ力学モデルで用いる過渡応答パラメータの値を算出する処理をコンピュータに実行させるプログラムであって、
スリップ角の時系列データをタイヤ測定条件としてタイヤに与えることにより計測されるコーナリング時のタイヤ回転軸に作用する横力およびセルフアライニングトルクの物理量の過渡応答の計測データを取得し、
タイヤ回転軸に作用する物理量の、スリップ角に対する応答関数であって、コーナリング時のタイヤの変形応答を再現するタイヤ力学モデルを用いて定めた応答関数に、過渡応答パラメータの値を付与し、前記値が付与された前記応答関数と前記タイヤ測定条件として用いた前記スリップ角の時系列データとを用いて積分を行うことによって、前記物理量の時系列データをタイヤの予測データとして算出し、
算出した前記予測データが許容範囲内で前記計測データに一致するまで、前記過渡応答パラメータの値の修正と、前記予測データの算出を繰り返し、
前記予測データが前記計測データに前記許容範囲内で一致したときに付与した前記過渡応答パラメータの値を、前記タイヤの過渡応答を再現する過渡応答パラメータの値として決定する、処理をコンピュータに実行させ、
前記タイヤ力学モデルを用いて定めた前記応答関数は、前記タイヤのベルト部材の曲げ変形の遅れ応答をむだ時間による遅れ応答として構成され、
前記応答関数は、下記式で表される横力Ｆｙ及びセルフアライニングトルクＭｚのスリップ角αに対する伝達関数Ｆｙ（ｓ）／α（ｓ）及びＭｚ（ｓ）／α（ｓ）を用いて得られる時間の関数を含み、
0005790018.tif 000032
0005790018.tif 000033
前記式中のＫｙ、Ａｓ、Ｋｔ、Ｔ０、Ｔ１、Ｔ２、及びＴ３は、前記過渡応答パラメータであり、Ｔ１は前記タイヤのベルト部材の曲げ変形の遅れ応答のむだ時間であり、ｓは、応答関数を定める複素数である、ことを特徴とするプログラム。
【請求項９】
スリップ角を時系列データとして与えたときのコーナリング時のタイヤの過渡応答データをタイヤ力学モデルを用いて算出する処理をコンピュータに実行させるプログラムであって、
スリップ角を時系列データを取得するとともに、コーナリング時のタイヤの変形応答を再現するタイヤ力学モデルを用いて定めた応答関数であって、タイヤ回転軸に作用する横力およびセルフアライニングトルクの物理量の、スリップ角に対する応答関数に用いる過渡応答パラメータの値を取得し、
前記応答関数に、前記過渡応答パラメータの値を付与し、前記値が付与された前記応答関数と取得した前記スリップ角の時系列データとを用いて積分を行うことによって、前記物理量の時系列データをタイヤの過渡応答データとして算出する、処理をコンピュータに実行させ、
前記タイヤ力学モデルを用いて定めた前記応答関数は、前記タイヤのベルト部材の曲げ変形の遅れ応答をむだ時間による遅れ応答として構成され、
前記応答関数は、下記式で表される横力Ｆｙ及びセルフアライニングトルクＭｚのスリップ角αに対する伝達関数Ｆｙ（ｓ）／α（ｓ）及びＭｚ（ｓ）／α（ｓ）を用いて得られる時間の関数を含み、
0005790018.tif 000034
0005790018.tif 000035
前記式中のＫｙ、Ａｓ、Ｋｔ、Ｔ０、Ｔ１、Ｔ２、及びＴ３は、前記過渡応答パラメータであり、Ｔ１は前記タイヤのベルト部材の曲げ変形の遅れ応答のむだ時間であり、ｓは、応答関数を定める複素数である、ことを特徴とするプログラム。
【請求項１０】
タイヤのコーナリング時の過渡応答の計測データからタイヤ力学モデルで用いる過渡応答パラメータの値を算出するデータ処理装置であって、
スリップ角の時系列データをタイヤ測定条件としてタイヤに与えた時に計測されるコーナリング時のタイヤ回転軸に作用する横力およびセルフアライニングトルクの物理量の過渡応答の計測データを取得するデータ取得部と、
タイヤ回転軸に作用する前記物理量の、スリップ角に対する応答関数であって、コーナリング時のタイヤの変形応答を再現するタイヤ力学モデルを用いて定めた応答関数に、過渡応答パラメータの値を付与し、前記値が付与された前記応答関数と前記タイヤ測定条件として用いた前記スリップ角の時系列データとを用いて積分を行うことによって、前記物理量の時系列データをタイヤの予測データとして算出する演算部と、
算出した前記予測データが許容範囲内で前記計測データに一致するまで、前記過渡応答パラメータの値の修正を行うことにより、前記演算部に前記予測データの算出を繰り返し実行させ、前記予測データが前記計測データに前記許容範囲内で一致したときに付与した前記過渡応答パラメータの値を、前記タイヤの過渡応答を再現する過渡応答パラメータの値として決定するパラメータ値決定部と、を有し、
前記タイヤ力学モデルを用いて定めた前記応答関数は、前記タイヤのベルト部材の曲げ変形の遅れ応答をむだ時間による遅れ応答として構成され、
前記応答関数は、下記式で表される横力Ｆｙ及びセルフアライニングトルクＭｚの
スリップ角αに対する伝達関数Ｆｙ（ｓ）／α（ｓ）及びＭｚ（ｓ）／α（ｓ）を用いて得られる時間の関数を含み、
0005790018.tif 000036
0005790018.tif 000037
前記式中のＫｙ、Ａｓ、Ｋｔ、Ｔ０、Ｔ１、Ｔ２、及びＴ３は、前記過渡応答パラメータであり、Ｔ１は前記タイヤのベルト部材の曲げ変形の遅れ応答のむだ時間であり、ｓは、応答関数を定める複素数である、ことを特徴とするデータ処理装置。
【請求項１１】
スリップ角を時系列データとして与えたときのコーナリング時のタイヤの過渡応答データをタイヤ力学モデルを用いて算出するデータ算出装置であって、
スリップ角の時系列データを取得するとともに、コーナリング時のタイヤの変形応答を再現するタイヤ力学モデルを用いて定めた応答関数であって、タイヤ回転軸に作用する横力およびセルフアライニングトルクの物理量の、スリップ角に対する応答関数に用いる過渡応答パラメータの値を取得するデータ取得部と、
前記応答関数に、前記過渡応答パラメータの値を付与し、前記値が付与された前記応答関数と取得した前記スリップ角の時系列データとを用いて積分を行うことによって、前記物理量の時系列データをタイヤの過渡応答データとして算出する演算部と、を有し、
前記タイヤ力学モデルを用いて定めた前記応答関数は、前記タイヤのベルト部材の曲げ変形の遅れ応答をむだ時間による遅れ応答として構成され、
前記応答関数は、下記式で表される横力Ｆｙ及びセルフアライニングトルクＭｚのスリップ角αに対する伝達関数Ｆｙ（ｓ）／α（ｓ）及びＭｚ（ｓ）／α（ｓ）を用いて得られる時間の関数を含み、
0005790018.tif 000038
0005790018.tif 000039
前記式中のＫｙ、Ａｓ、Ｋｔ、Ｔ０、Ｔ１、Ｔ２、及びＴ３は、前記過渡応答パラメータであり、Ｔ１は前記タイヤのベルト部材の曲げ変形の遅れ応答のむだ時間であり、ｓは、応答関数を定める複素数である、ことを特徴とするデータ算出装置。  
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、スリップ角を時系列データとして与えたときのコーナリング時のタイヤの過渡応答データをタイヤ力学モデルを用いて算出するデータ算出装置の行うタイヤの過渡応答データの算出方法、算出装置及びこの算出方法をコンピュータに実行させるプログラムと、タイヤのコーナリング時の過渡応答の計測データからタイヤ力学モデルで用いる過渡応答パラメータの値を算出するデータ処理方法、データ処理装置及びこのデータ処理方法をコンピュータに実行させるプログラムに関する。
【背景技術】
【０００２】
現在、自動車産業において、車両の安全走行、危険回避のための高度な車両制御が求められている。タイヤは、車両と路面との間に介在し、路面からの力を唯一車両に伝達するものであるため、車両の安全走行、高度な車両制御にとってタイヤの役割は重要である。このためタイヤのコーナリング特性の解析が必要である。
【０００３】
タイヤのコーナリング特性を、タイヤ力学モデルを用いて再現することは、タイヤ力学モデルに用いる各種ばね特性等のタイヤ力学要素パラメータを用いてコーナリング特性を議論することができるので、タイヤ設計者にとって有効である。例えば、タイヤ力学モデルを用いてコーナリング中のタイヤの過渡応答データを算出し、また、タイヤ力学モデルに用いる過渡応答を定めるパラメータの値を決定する方法が知られている（特許文献１）。
【０００４】
上記方法では、具体的には、タイヤ力学モデルを構成するタイヤ力学要素パラメータの値を取得してコーナリング中のトレッド部の変形応答を規定する１次遅れ応答の応答関数とタイヤ力学モデルに与えたスリップ角の時系列データの時間変化量との畳み込み積分を行ってタイヤ力学モデルにおいて路面に対するトレッド部のスリップ角の過渡応答の時系列データを算出し、この過渡応答の時系列データに基づいて横力、セルフアライニングトルクの出力データを、コーナリングの過渡応答として算出する。また、算出した出力データが、実測データに略一致するように遅延時定数を修正することで、過渡応答を定めるパラメータの値を決定する。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００７−２９０６９４号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかし、上記コーナリング中の過渡応答データを算出する方法、及びタイヤ力学モデルに用いる過渡応答を定めるパラメータの値を決定する方法では、いずれもタイヤ力学モデルから導出される１次遅れ応答の応答関数を用いるので、算出される過渡応答データは、応答周波数が高いほど計測データとの乖離が大きくなる。一方、過渡応答データを算出するとき、及びタイヤ力学モデルに用いる過渡応答を定めるパラメータの値を決定するとき、上記１次遅れ応答に代えて高次の遅れ応答を含んだ応答関数を用いることもできるが、高次の遅れ応答を含んだ応答関数をタイヤ力学モデルからどのように導出すればよいのか、不明である。このため、タイヤ力学モデルにおける力学要素パラメータの値と、高次の遅れ応答を含んだ応答関数に用いる過渡応答パラメータの値との間の関係は不明である。
【０００７】
そこで、本発明は、上記従来技術の問題を解決するために、算出される予測データと計測データとの乖離が少ない高次の遅れ応答を含み、タイヤ力学要素パラメータと、高次の遅れ応答を含んだ応答関数に用いる過渡応答パラメータとを関係付けたタイヤ力学モデルを用いてタイヤの過渡応答データを算出する方法及び装置と、このタイヤ力学モデルを用いたデータ処理方法及び装置と、プログラムを提供することを目的とする。
【課題を解決するための手段】
【０００８】
本発明の一態様は、タイヤのコーナリング時の過渡応答の計測データからタイヤ力学モデルで用いる過渡応答パラメータの値を算出することを、コンピュータに実行させるデータ処理方法である。当該方法は、
スリップ角の時系列データをタイヤ測定条件としてタイヤに与えた時に計測されるコーナリング時のタイヤ回転軸に作用する横力およびセルフアライニングトルクの物理量の過渡応答の計測データをコンピュータに取得させるステップと、
タイヤ回転軸に作用する前記物理量の、スリップ角に対する応答関数であって、コーナリング時のタイヤの変形応答を再現するタイヤ力学モデルを用いて定めた応答関数に、過渡応答パラメータの値を付与し、前記値が付与された前記応答関数と前記タイヤ測定条件として用いた前記スリップ角の時系列データとを用いて積分を行うことによって、前記物理量の時系列データをタイヤの予測データとしてコンピュータに算出させるステップと、
算出した前記予測データが許容範囲内で前記計測データに一致するまで、前記過渡応答パラメータの値の修正と、前記予測データの算出をコンピュータに繰り返し行わせるステップと、
前記予測データが前記計測データに前記許容範囲内で一致したときに付与した前記過渡応答パラメータの値を、前記タイヤの過渡応答を再現する過渡応答パラメータの値としてコンピュータに決定させるステップと、を有し、
前記タイヤ力学モデルを用いて定めた前記応答関数は、前記タイヤのベルト部材の曲げ変形の遅れ応答をむだ時間による遅れ応答として構成され、
前記応答関数は、下記式で表される横力Ｆｙ及びセルフアライニングトルクＭｚのスリップ角αに対する伝達関数Ｆｙ（ｓ）／α（ｓ）及びＭｚ（ｓ）／α（ｓ）を用いて得られる時間の関数を含み、
0005790018.tif 000002
0005790018.tif 000003
前記式中のＫｙ、Ａｓ、Ｋｔ、Ｔ０、Ｔ１、Ｔ２、及びＴ３は、前記過渡応答パラメータであり、Ｔ１は前記タイヤのベルト部材の曲げ変形の遅れ応答のむだ時間であり、ｓは、応答関数を定める複素数である。
【０００９】
本発明の他の一態様は、スリップ角を時系列データとして与えたときのコーナリング時のタイヤの過渡応答データをタイヤ力学モデルを用いて算出することを、コンピュータに実行させるタイヤの過渡応答データの算出方法であって、
スリップ角の時系列データを取得するとともに、コーナリング時のタイヤの変形応答を再現するタイヤ力学モデルを用いて定めた応答関数であって、タイヤ回転軸に作用する横力およびセルフアライニングトルクの物理量の、スリップ角に対する応答関数に用いる過渡応答パラメータの値をコンピュータに取得させるステップと、
前記応答関数に前記過渡応答パラメータの値を付与し、前記値が付与された前記応答関数と取得した前記スリップ角の時系列データとを用いて積分を行うことによって、前記物理量の時系列データをタイヤの過渡応答データとしてコンピュータに算出させるステップと、を有し、
前記タイヤ力学モデルを用いて定めた前記応答関数は、前記タイヤのベルト部材の曲げ変形の遅れ応答をむだ時間による遅れ応答として構成され、
前記応答関数は、下記式で表される横力Ｆｙ及びセルフアライニングトルクＭｚのスリップ角αに対する伝達関数Ｆｙ（ｓ）／α（ｓ）及びＭｚ（ｓ）／α（ｓ）を用いて得られる時間の関数を含み、
0005790018.tif 000004
0005790018.tif 000005
前記式中のＫｙ、Ａｓ、Ｋｔ、Ｔ０、Ｔ１、Ｔ２、及びＴ３は、前記過渡応答パラメータであり、Ｔ１は前記タイヤのベルト部材の曲げ変形の遅れ応答のむだ時間であり、ｓは、応答関数を定める複素数である。
【００１０】
本発明のさらに他の一態様は、タイヤのコーナリング時の過渡応答の計測データからタイヤ力学モデルで用いる過渡応答パラメータの値を算出する処理をコンピュータに実行させるプログラムである。当該プログラムは、
スリップ角の時系列データをタイヤ測定条件としてタイヤに与えることにより計測されるコーナリング時のタイヤ回転軸に作用する横力およびセルフアライニングトルクの物理量の過渡応答の計測データを取得し、
タイヤ回転軸に作用する物理量の、スリップ角に対する応答関数であって、コーナリング時のタイヤの変形応答を再現するタイヤ力学モデルを用いて定めた応答関数に、過渡応答パラメータの値を付与し、前記値が付与された前記応答関数と前記タイヤ測定条件として用いた前記スリップ角の時系列データとを用いて積分を行うことによって、前記物理量の時系列データをタイヤの予測データとして算出し、
算出した前記予測データが許容範囲内で前記計測データに一致するまで、前記過渡応答パラメータの値の修正と、前記予測データの算出を繰り返し、
前記予測データが前記計測データに前記許容範囲内で一致したときに付与した前記過渡応答パラメータの値を、前記タイヤの過渡応答を再現する過渡応答パラメータの値として決定する、処理をコンピュータに実行させ、
前記タイヤ力学モデルを用いて定めた前記応答関数は、前記タイヤのベルト部材の曲げ変形の遅れ応答をむだ時間による遅れ応答として構成され、
前記応答関数は、下記式で表される横力Ｆｙ及びセルフアライニングトルクＭｚのスリップ角αに対する伝達関数Ｆｙ（ｓ）／α（ｓ）及びＭｚ（ｓ）／α（ｓ）を用いて得られる時間の関数を含み、
0005790018.tif 000006
0005790018.tif 000007
前記式中のＫｙ、Ａｓ、Ｋｔ、Ｔ０、Ｔ１、Ｔ２、及びＴ３
は、前記過渡応答パラメータであり、Ｔ１は前記タイヤのベルト部材の曲げ変形の遅れ応答のむだ時間であり、ｓは、応答関数を定める複素数である。
【００１１】
本発明のさらに他の一態様は、スリップ角を時系列データとして与えたときのコーナリング時のタイヤの過渡応答データをタイヤ力学モデルを用いて算出する処理をコンピュータに実行させるプログラムである。当該プログラムは、
スリップ角を時系列データを取得するとともに、コーナリング時のタイヤの変形応答を再現するタイヤ力学モデルを用いて定めた応答関数であって、タイヤ回転軸に作用する横力およびセルフアライニングトルクの物理量の、スリップ角に対する応答関数に用いる過渡応答パラメータの値を取得し、
前記応答関数に、前記過渡応答パラメータの値を付与し、前記値が付与された前記応答関数と取得した前記スリップ角の時系列データとを用いて積分を行うことによって、前記物理量の時系列データをタイヤの過渡応答データとして算出する、処理をコンピュータに実行させ、
前記タイヤ力学モデルを用いて定めた前記応答関数は、前記タイヤのベルト部材の曲げ変形の遅れ応答をむだ時間による遅れ応答として構成され、
前記応答関数は、下記式で表される横力Ｆｙ及びセルフアライニングトルクＭｚのスリップ角αに対する伝達関数Ｆｙ（ｓ）／α（ｓ）及びＭｚ（ｓ）／α（ｓ）を用いて得られる時間の関数を含み、
0005790018.tif 000008
0005790018.tif 000009
前記式中のＫｙ、Ａｓ、Ｋｔ、Ｔ０、Ｔ１、Ｔ２、及びＴ３は、前記過渡応答パラメータであり、Ｔ１は前記タイヤのベルト部材の曲げ変形の遅れ応答のむだ時間であり、ｓは、応答関数を定める複素数である。
【００１２】
また、本発明の他の一態様は、タイヤのコーナリング時の過渡応答の計測データからタイヤ力学モデルで用いる過渡応答パラメータの値を算出するデータ処理装置である。当該装置は、
スリップ角の時系列データをタイヤ測定条件としてタイヤに与えた時に計測されるコーナリング時のタイヤ回転軸に作用する横力およびセルフアライニングトルクの物理量の過渡応答の計測データを取得するデータ取得部と、
タイヤ回転軸に作用する前記物理量の、スリップ角に対する応答関数であって、コーナリング時のタイヤの変形応答を再現するタイヤ力学モデルを用いて定めた応答関数に、過渡応答パラメータの値を付与し、前記値が付与された前記応答関数と前記タイヤ測定条件として用いた前記スリップ角の時系列データとを用いて積分を行うことによって、前記物理量の時系列データをタイヤの予測データとして算出する演算部と、
算出した前記予測データが許容範囲内で前記計測データに一致するまで、前記過渡応答パラメータの値の修正を行うことにより、前記演算部に前記予測データの算出を繰り返し実行させ、前記予測データが前記計測データに前記許容範囲内で一致したときに付与した前記過渡応答パラメータの値を、前記タイヤの過渡応答を再現する過渡応答パラメータの値として決定するパラメータ値決定部と、を有し、
前記タイヤ力学モデルを用いて定めた前記応答関数は、前記タイヤのベルト部材の曲げ変形の遅れ応答をむだ時間による遅れ応答として構成され、
前記応答関数は、下記式で表される横力Ｆｙ及びセルフアライニングトルクＭｚのスリップ角αに対する伝達関数Ｆｙ（ｓ）／α（ｓ）及びＭｚ（ｓ）／α（ｓ）を用いて得られる時間の関数を含み、
0005790018.tif 000010
0005790018.tif 000011
前記式中のＫｙ、Ａｓ、Ｋｔ、Ｔ０、Ｔ１、Ｔ２、及びＴ３は、前記過渡応答パラメータであり、Ｔ１は前記タイヤのベルト部材の曲げ変形の遅れ応答のむだ時間であり、ｓは、応答関数を定める複素数である。
【００１３】
さらに、本発明の他の一態様は、スリップ角を時系列データとして与えたときのコーナリング時のタイヤの過渡応答データをタイヤ力学モデルを用いて算出するデータ算出装置である。当該装置は、
スリップ角の時系列データを取得するとともに、コーナリング時のタイヤの変形応答を再現するタイヤ力学モデルを用いて定めた応答関数であって、タイヤ回転軸に作用する横力およびセルフアライニングトルクの物理量の、スリップ角に対する応答関数に用いる過渡応答パラメータの値を取得するデータ取得部と、
前記応答関数に、前記過渡応答パラメータの値を付与し、前記値が付与された前記応答関数と取得した前記スリップ角の時系列データとを用いて積分を行うことによって、前記物理量の時系列データをタイヤの過渡応答データとして算出する演算部と、を有し、
前記タイヤ力学モデルを用いて定めた前記応答関数は、前記タイヤのベルト部材の曲げ変形の遅れ応答をむだ時間による遅れ応答として構成され、
前記応答関数は、下記式で表される横力Ｆｙ及びセルフアライニングトルクＭｚの
スリップ角αに対する伝達関数Ｆｙ（ｓ）／α（ｓ）及びＭｚ（ｓ）／α（ｓ）を用いて得られる時間の関数を含み、
0005790018.tif 000012
0005790018.tif 000013
前記式中のＫｙ、Ａｓ、Ｋｔ、Ｔ０、Ｔ１、Ｔ２、及びＴ３は、前記過渡応答パラメータであり、Ｔ１は前記タイヤのベルト部材の曲げ変形の遅れ応答のむだ時間であり、ｓは、応答関数を定める複素数である。
【発明の効果】
【００１４】
上記タイヤの過渡応答データの算出方法及び算出装置と、データ処理方法及びデータ処理装置と、プログラムによれば、算出される予測データと計測データとの乖離が少なく、力学要素パラメータと、高次の遅れ応答を含んだ応答関数に用いる過渡応答パラメータとを容易に関係付けることができる。
【図面の簡単な説明】
【００１５】
【図１】（ａ）は、本実施形態のタイヤの過渡応答データの算出方法を実施するデータ算出装置および本実施形態のデータ処理方法を実施するデータ処理装置を兼ねる装置の構成を示す図であり、（ｂ）は、（ａ）に示す装置の機能ブロック図である。
【図２】図１（ｂ）に示す装置のタイヤ力学モデル演算部を説明する図である。
【図３】（ａ），（ｂ）は、図２に示すタイヤ力学モデル演算部で用いるむだ時間を説明する図である。
【図４】（ａ），（ｂ）は、スリップ角が付いたタイヤのトレッド部材に作用するせん断力の様子を説明する図であり、（ｃ）は、タイヤの変形と横力及びセルフアライニングトルクとの関係を説明する図である。
【図５】図２に示すタイヤ力学モデル演算部で用いる横力及びセルフアライニングトルクの伝達関数の制御ブロック図である。
【図６】本実施形態のデータ処理方法の一例のフローを示すフローチャートである。
【図７】本実施形態のデータ算出方法の一例のフローを示すフローチャートである。
【図８】タイヤの横力の計測データの一例を示す図である。
【図９】本実施形態で得られる予測データと計測データとを示す図である。
【図１０】（ａ），（ｂ）は、本実施形態で得られる予測データから得られる特性の一例を、計測データから得られる特性とともに示す図である。
【図１１】本実施形態で得られる予測データから得られる特性の他の一例を示す図である。
【図１２】（ａ）〜（ｅ）は、複数のタイヤにおける、本実施形態で得られるタイヤ力学要素パラメータの値と特性に関する結果を示す図である。
【発明を実施するための形態】
【００１６】
以下、本発明のタイヤの過渡応答データの算出方法及び算出装置と、データ処理方法及びデータ処理装置と、プログラムについて詳細に説明する。
【００１７】
一般に、転動するタイヤにスリップ角が付いたとき、タイヤのトレッド部材の接地面には、スリップ角に応じたせん断が発生し、このせん断に応じて横力がタイヤ回転軸に発生し、さらにこのとき接地面に発生するせん断力の発生分布に応じたセルフアライニングトルクがタイヤ回転軸に発生する。具体的には、スリップ角が与えられたタイヤのトレッド部材には、ベルト部材に対する路面の相対変位に応じたせん断が発生するが、このとき、トレッド部材の位置に影響を与えるベルト部材、さらには、ベルト部材の位置に影響を与えるサイドウォール部材も変形して変位する。このため、横力及びセルフアライニングトルクは、発生した横力あるいはセルフアライニングトルクの影響を再帰的に受けることになる。すなわち、タイヤにスリップ角が与えられたとき、発生する横力及びセルフアライニングトルクはフィードバック制御を受ける。したがって、時系列でスリップ角が変化して横力及びセルフアライニングトルクが発生するとき、横力及びセルフアライニングトルクの応答には、ベルト部材の横曲げ変形およびサイドウォール部材のねじり変形に起因する遅れ応答が含まれる。また、タイヤのトレッド部材のゴムは、変形速度が速くなるほど硬くなる粘弾性特性を有するので、発生する横力及びセルフアライニングトルクには、タイヤのトレッド部材の上記特性に起因する応答が含まれる。
【００１８】
したがって、本実施形態のタイヤ力学モデルで表される変形応答は、タイヤのトレッド部材の変形応答、タイヤのベルト部材の曲げ変形の遅れ応答、及びタイヤのサイドウォール部材のねじり変形の遅れ応答を含む。このとき、上記変形応答は、タイヤのベルト部材の曲げ変形の遅れ応答をむだ時間による遅れ応答として構成されている。このため、本実施形態で用いる上記３つの応答を含んだステップ応答関数は、比較的単純な形で表され、さらに、ステップ応答関数に用いる過渡応答パラメータ（Ｋy，Ｔ0，Ｔ1，Ｔ2，Ａs，Ｋt，Ｔ3）とタイヤ力学モデルに用いるタイヤ力学要素パラメータ（Ｋdr，Ｋrr，ＫLr）との間の関係は式で表され得る。したがって、計測データを用いて過渡応答パラメータの値を算出したとき、この過渡応答パラメータの値からタイヤ力学要素パラメータ（Ｋdr，Ｋrr，ＫLr）の値に変換することができる。これにより、タイヤ設計者にとって、どのばね定数を修正すれば、よりよい横力及びセルフアライニングトルクの応答を得ることができるか直感的に理解することができる。
以下、本実施形態について詳細に説明する。
【００１９】
（データ算出装置、データ処理装置）
図１（ａ）は、タイヤの過渡応答データの算出方法を実施するデータ算出装置と、データ処理方法を実施するデータ処理装置とを兼ねた装置１０の構成を示す図である。装置１０は、ＣＰＵ（Central Processing Unit）１２と、ＲＯＭ（Read Only Memory）１４と、ＲＡＭ（Random Access Memory）あるいはハードディスク等のメモリ１６と、入出力部１８と、を有するコンピュータである。ＣＰＵ１２と、ＲＯＭ１４と、メモリ１６と、入出力部１８とは、互いにバス１９で接続されている。入出力部１８には、マウスやキーボード等の入力操作系２０と、ディスプレイ２２と、計測装置２４と接続されている。
【００２０】
ＣＰＵ１２は、装置１０の各部分の制御と動作管理を行う。ＲＯＭ１４は、装置１０が立ち上がるとき、ＯＳ（Operating System）として呼び出されるプログラムを記録する。メモリ１６は、本実施形態のタイヤの過渡応答データの算出方法およびデータ処理方法を実行するプログラムを記録し、また、本実施形態のタイヤの過渡応答データの算出方法およびデータ処理方法を実行するときに用いる各種データや情報の他、処理結果のデータを記憶する。
計測装置２４は、コーナリング時のタイヤ回転軸に作用する横力あるいはセルフアライニングトルクの時系列データを計測データとして計測し、この計測データを入出力部１８に出力する。
このような装置１０は、メモリ１６に記録されているプログラムを呼び出して実行することにより、タイヤの過渡応答データの算出方法およびデータ処理方法を実行する装置として機能する。
【００２１】
図１（ｂ）は、メモリ１６に記録されている上記プログラムを呼び出した装置１０の機能ブロック図である。
装置１０は、データ取得部２６と、パラメータ値決定部３０と、Ｆｙ，Ｍｚ算出部３２と、タイヤ力学モデル演算部３４と、を有する。これらの各部分は、ＣＰＵ１２が実質的に各部分の機能を担うソフトウェアモジュールである。
【００２２】
装置１０は、オペレータの指示に応じて切り替える、過渡応答算出処理モードと、データ処理モードと、を備える。
過渡応答算出モードでは、装置１０は、スリップ角を時系列データとして与えたときのコーナリング時のタイヤの過渡応答データを、タイヤ力学モデルを用いて定めた応答関数を用いて算出する。データ処理モードでは、装置１０は、タイヤ力学モデルを用いて、タイヤのコーナリング時の過渡応答の計測データからタイヤ力学モデルで用いる過渡応答パラメータの値を算出する。
【００２３】
データ取得部２６は、過渡応答算出モードでは、スリップ角の時系列データを取得し、さらに、後述するタイヤ力学モデルから導出されるステップ応答関数に用いる過渡応答パラメータの値を取得する。ステップ応答関数は、コーナリング時のタイヤの変形応答を再現するタイヤ力学モデルにおける応答関数であって、タイヤ回転軸に作用する横力およびセルフアライニングトルクの、スリップ角に対する応答関数である。データ取得部２６は、取得したスリップ角の時系列データを、タイヤ力学モデル演算部３４が過渡応答の時系列データを算出できるように、タイヤ力学モデル演算部３４に送る。また、データ取得部２６は、取得した過渡応答パラメータの値をパラメータ値決定部３０に送る。
データ取得部２６は、データ処理モードでは、タイヤのコーナリング時のタイヤ回転軸に作用する横力およびセルフアライニングトルクの過渡応答の計測データ、およびこの計測においてタイヤ測定条件として用いたスリップ角の時系列データを取得する。データ取得部２６は、取得した計測データを、パラメータ値決定部３０に送る。データ取得部２６は、取得したスリップ角の時系列データを、タイヤ力学モデル演算部３４が過渡応答の時系列データを算出できるように、タイヤ力学モデル演算部３４に送る。
【００２４】
パラメータ値決定部３０は、データ処理モードでは、設定した過渡応答パラメータの値をタイヤ力学モデル演算部３４に与えることで算出される横力及びセルフアライニングトルクの予測データが、データ取得部２６で得られた計測データに、許容範囲内で一致するか否かを判定する。判定の結果、予測データが計測データに対して許容範囲内で一致しない場合、パラメータ値決定部３０は、設定された過渡応答パラメータの値を修正し、修正した値をタイヤ力学モデル演算部３４に与える。パラメータ値決定部３０は、予測データが計測データに対して許容範囲内で一致するまで、過渡応答パラメータの値の修正を繰り返し行う。
判定の結果、予め設定された許容範囲内で予測データが計測データに一致する場合、このとき設定された過渡応答パラメータの値を、計測に用いたタイヤの過渡応答を再現する過渡応答パラメータの値として決定してディスプレイ２２や図示されないプリンタに出力する。さらに、パラメータ値決定部３０は、決定した過渡応答パラメータの値からタイヤ力学要素パラメータ（Ｋdr，Ｋrr，ＫLr）の値を算出してディスプレイ２２や図示されないプリンタに出力する。過渡応答パラメータの値からタイヤ力学要素パラメータ（Ｋdr，Ｋrr，ＫLr）の値の算出については後述する。
【００２５】
上記データ処理モードでは、最初に設定される過渡応答パラメータの値、すなわち初期値は、予め設定されたデフォルト値でもよいし、以前決定された過渡応答パラメータの値であってもよい。
また、予測データが計測データと許容範囲内で一致するとは、それぞれの時系列データの差分の二乗和が許容範囲として設定された閾値以下となることである。したがって、パラメータ値決定部３０における判定では、予測データと計測データの差分の二乗和が用いられる。
また、予測データが計測データと許容範囲内で一致するまで、過渡応答パラメータの値を繰り返し修正する方法は特に制限されず、公知の方法が用いられ、例えば、ニュートンラフソン法を用いることができる。
【００２６】
一方、過渡応答算出モードでは、パラメータ値決定部３０は、データ取得部２６から送られた過渡応答パラメータの値を、タイヤ力学モデル演算部３４に送る。
【００２７】
Ｆｙ，Ｍｚ算出部３２は、タイヤ力学モデル演算部３４で算出された横力およびセルフアライニングトルクの時系列データである予測データを、所定の形式に修正してディプレイ２２あるいは図示されないプリンタに出力する。
【００２８】
タイヤ力学モデル演算部３４は、データ取得部２６から送られたスリップ角の時系列データと、パラメータ値決定部３０から送られた過渡応答パラメータの値と、設定されたステップ応答関数を用いて、スリップ角の時系列データに対する横力及びセルフアライニングトルクの時系列データを予測データとして算出する。タイヤ力学モデル演算部３４がパラメータ値決定部３０から過渡応答パラメータの値を受けることは、パラメータ値決定部３０から、過渡応答の時系列データを算出する指示を受けることを意味する。このため、タイヤ力学モデル演算部３４は、パラメータ値決定部３０から過渡応答パラメータの値を受けると、直ちに、以下に示す横力及びセルフアライニングトルクの予測データの算出を開始する。
【００２９】
図２は、タイヤ力学モデル演算部３４を説明する図である。
タイヤ力学モデル演算部３４は、具体的に、図２に示されるように、パラメータ値決定部３０から受けた過渡応答パラメータ（Ｋy，Ｔ0，Ｔ1，Ｔ2，Ａs，Ｋt，Ｔ3）の値をタイヤ力学モデル演算部３４に設定されているステップ応答関数に付与し、値が付与されたステップ応答関数と、データ取得部２６から送られたスリップ角の時系列データα（ｔ）の時間変化量との畳み込み積分を行うことによって、横力及びセルフアライニングトルクの時系列データであるＦy（ｔ）、Ｍz（ｔ）をタイヤの過渡応答データ（予測データ）として算出する。
【００３０】
より具体的には、ステップ応答関数は、ラプラス演算子sを用いて表した図２中の伝達関数Ｆy（ｓ）／α（ｓ）、伝達関数Ｍz（ｓ）／α（ｓ）（α（ｓ）はラプラス演算子で表されるスリップ角）に、スリップ角のステップ入力をラプラス変換した関数を乗算し、この乗算結果に時間領域に変換した関数である。ステップ応答関数は、後述するように、式（１５）及び式（１６）で表されるＫst（ｔ），Ａst（ｔ）で表される。このステップ応答関数に、過渡応答パラメータ（Ｋy，Ｔ0，Ｔ1，Ｔ2，Ａs，Ｋt，Ｔ3）の値が付与される。さらに、タイヤ力学モデル演算部３４は、入力したスリップ角の時系列データの時間変化量（ｄα（ｔ’）／ｄｔ’）とＫst（ｔ−ｔ’），Ａst（ｔ−ｔ’）を乗算して積分を行う（ｔ’は積分変数）。これにより、タイヤ力学モデル演算部３４は、Ｆy（ｔ），Ｍz（ｔ）を算出する。算出されたＦy（ｔ），Ｍz（ｔ）は、データ処理モードではパラメータ値決定部２８に、過渡応答算出モードではＦy，Ｍz算出部３２に送られる。
【００３１】
図２に示すように、ステップ応答関数を定める基準となる伝達関数Ｆy（ｓ）／α（ｓ）では、Ｔ1に関する項がｅｘｐ（−Ｔ1ｓ）で表されている。この項はタイヤのベルト部材の曲げ変形の遅れ応答を、むだ時間であるＴ1による遅れ応答を表している。図３（ａ），（ｂ）はむだ時間を説明する図である。むだ時間であるＴ1による遅れ応答とは、図３（ａ）に示すように、時間０において、スリップ角０度が１度に変化してそれ以降１度を維持するスリップ角αのステップ関数に対して、発生する横力Ｆyは、図３（ｂ）に示すように、横力の立ち上がりがＴ1秒遅れる応答をいう。このような遅れ応答を用いることで、時間Ｔ1だけ遅れて立ち上った以降の応答は、伝達関数Ｆy（ｓ）／α（ｓ）のｅｘｐ（−Ｔ1ｓ）の項を除いた残りの項の応答、すなわち、２次の遅れ応答となる。２次の遅れ応答の応答関数は、後述するように式（１５）、式（１６）により陽に書き表すことができるので、タイヤ力学モデル演算部３４において予め式（１５）、式（１６）を用いてステップ応答関数を用意することができる。したがって、本実施形態の装置１０は、高次の遅れ応答を用いて容易にＦy（ｔ），Ｍz（ｔ）を算出することができる。
【００３２】
（タイヤ力学モデル）
図２に示すタイヤ力学モデル演算部３４で用いるステップ応答関数は、タイヤ力学モデルに基づいて導出されたものである。このタイヤ力学モデルについて説明する。
図４（ａ），（ｂ）は、スリップ角が付いたタイヤのトレッド部材に作用するせん断力の様子を説明する図である。タイヤのトレッド部材が地面と接地する接地面Ｐにおいて、ホイールの中心の方向とタイヤの移動方向との間にスリップ角度αの分ずれがあるため、トレッド部材にはせん断力が発生する。このせん断力は、接地面Ｐの踏み込み端から蹴り出し端に向かって徐々に大きくなる分布を有する。これによって図４（ｂ）に示すように、横力Ｆｙ（せん断力の合計）とセルフアライニングトルクＭy（タイヤ中心回りのモーメント）がタイヤ回転軸に発生する。しかし、このとき、ベルト部材は横力Ｆｙに応じて横力Ｆｙの向きに追従した横曲げの変形（ベルトがタイヤ幅方向に曲がる変形）をする。このため、横力Ｆｙの領域Ｒ2は、図４（ｃ）に示すように、図４（ｂ）に示す領域Ｒ1に対して小さくなる。このときの横力Ｆｙによるベルト部材の横曲げ変形による変位量を（εｌ／３）・Ｆｙとする（εｌ／３は、ベルト部材の横曲げ剛性の逆数を表す定数）。
【００３３】
さらに、図４（ｃ）に示すように、このとき発生するセルフアライニングトルクＭzはサイドウォール部材のねじり変形によりスリップ角αを小さくするように作用する。このときのスリップ角の低減量は、Ｍz／Ｇmzである。
したがって、横力Ｆｙは下記式（５）に示すように表され、セルフアライニングトルクＭzは下記式（６）に示すように表される。ここで、Ｋyo，Ａs0は、定数である。
【００３４】
【数５】
0005790018.tif 000014
【００３５】
【数６】
0005790018.tif 000015
【００３６】
さらに、タイヤにスリップ角が時系列データとして入力されたとき、横力Ｆｙ、セルフアライニングトルクＭzは過渡応答を示す。この過渡応答は、式（５）、（６）を変形して下記式（７），（８）で表される。ここで、Ｋyo，Ａs0をＫy，Ａsと書き直している。式（５）、（６）から式（７）、（８）への変形は、スリップ角が変化したとき、タイヤの変形がタイヤ周上全体に影響を及ぼすのに必要な時間はＬ／ｖ（Ｌはタイヤ周長、ｖはタイヤの回転速度）であるので、この時間に基づいて、ベルト部材及びサイドウォール部材の変形が遅れることにより横力Ｆｙ、セルフアライニングトルクＭzの遅延応答が起こると考えることに基づいている。
【００３７】
【数７】
0005790018.tif 000016
【００３８】
【数８】
0005790018.tif 000017
【００３９】
したがって、上記式（７），（８）をラプラス変換して、ラプラス演算子ｓを用いて伝達関数を表すと、下記式（９），（１０）となる。
【００４０】
【数９】
0005790018.tif 000018
【００４１】
【数１０】
0005790018.tif 000019
【００４２】
さらに、タイヤのトレッド部材のゴムは、変形周波数が高くなるほど硬くなる粘弾性特性を有するので、Ｋy，Ａsを下記式（１１），（１２）に示すように修正することにより、下記式（１３），（１４）を得る。式（１１）、式（１２）におけるＴ0は、トレッド部材の変形応答の時定数である。このとき、式（１４）には、スリップ角が生じた瞬間に値が最大となり、スリップ角の変化が小さくなると値が０になるセルフアライニングトルクの周知のねじり変形要素（ねじりトルク）を第２項として加えている。ここで、ねじり変形要素における時定数をＴ3とする。Ｋtは定数。
【００４３】
【数１１】
0005790018.tif 000020
【００４４】
【数１２】
0005790018.tif 000021
【００４５】
【数１３】
0005790018.tif 000022
【００４６】
【数１４】
0005790018.tif 000023
【００４７】
横力に関する伝達関数の式（１３）の上段の右辺中の分母はｓの４次式となっているので、式（１３）の上段で表される応答は４次の遅れ応答である。しかし、４次の遅れ応答を最下段の式に示すように近似することができる。このとき、ｅｘｐ（−Ｔ1ｓ）の項がむだ時間であるＴ1による遅れ応答として構成されている。これにより、ｅｘｐ（−Ｔ1ｓ）以外の項は２次の遅れ応答となるので、横力のステップ応答関数は、２次の遅れ応答のステップ応答関数とむだ時間Ｔ1とを用いて、下記式（１５）に示すように表される。また、セルフアライニングトルクのステップ応答関数は、下記式（１６）に示すように表される。
【００４８】
【数１５】
0005790018.tif 000024
【００４９】
【数１６】
0005790018.tif 000025
【００５０】
さらに、タイヤ力学モデル演算部３４は、式（１５）及び式（１６）で示されるステップ応答関数Ｋst（ｔ）及びＡst（ｔ）と、データ取得部２６から送られたスリップ角の時系列データα（ｔ）の時間変化量との畳み込み積分を、下記式（１７）及び（１８）に示すように行う。これにより、タイヤ力学モデル演算部３８は横力及びセルフアライニングトルクの時系列データＦｙ(ｔ)，Ｍｚ(ｔ)である予測データを算出する。算出した予測データは、データ処理モードではパラメータ値決定部３０に送られ、過渡応答算出モードではＦy，Ｍz算出部３２に送られる。
【００５１】
【数１７】
0005790018.tif 000026
【００５２】
【数１８】
0005790018.tif 000027
【００５３】
図５は、式（１３），（１４）で表される伝達関数の制御ブロック図を示す。タイヤ力学モデルは、このように、トレッド部材の変形応答、ベルト部材の曲げ変形の遅れ応答、およびサイドウォール部材のねじり変形の遅れ応答を、タイヤの変形応答として少なくとも含む。
【００５４】
ここで、図４（ｃ）に示すベルト部材の変位量（εｌ／３）・Ｆｙを、横力Ｆｙを見かけ上の横曲げ剛性で除算することにより得られる変位量と考えることにより、見かけ上の横曲げ剛性、すなわちベルト部材の変形による横曲げ剛性Ｋdrを下記式のように定めることができる。
Ｋdr  ＝  ３／（（εｌ）・Ｌ）    （Ｌはタイヤの周長）
【００５５】
また、図４（ｃ）に示すサイドウォール部材のねじり変形によるスリップ角の低下量Ｍz／Ｇmzを、セルフアライニングトルクＭzを見かけ上のねじり剛性で除算することにより得られる変位相当量と考えることにより、見かけ上のねじり剛性、すなわちサイドウォール部材のねじり剛性Ｋrrを下記式のように定めることができる。
Ｋrr  ＝  Ｇmz／Ｌ      （Ｌはタイヤの周長）
【００５６】
さらに、横ばね定数ＫLrを下記式によって定めることができる。
ＫLr  ＝  （１／Ｋdr  ＋  （１／Ｋrr  ）・（Ａs／Ｋy）)−１
この横ばね定数ＫLrは、非転動時のタイヤが接地した接地面を横方向に変位させたときに得られるタイヤの周知の横ばね定数に対応する。
なお、Ｋdrは、Ｔ1と式（９）中に示したようなＴ1＝Ｋy・（εｌ／３）・Ｌ／ｖの関係があるので、Ｔ1＝Ｋy／（Ｋdr  ・ｖ）となる。また、Ｋrrは、Ｔ1と式（１０）中に示したようなＴ2＝Ａs・Ｌ／（Ｇmz・ｖ）の関係があるので、Ｔ2＝Ａs／（Ｋrr  ・ｖ）となる。
したがって、パラメータ値決定部３０において計測データと予測データが許容範囲で一致することによって決定された過渡応答パラメータＴ1，Ｔ2，Ｋy，Ａsの値は、タイヤ力学要素パラメータであるベルト部材の変形による横曲げ剛性Ｋdr、サイドウォール部材のねじり剛性Ｋrr、さらには、横ばね定数ＫLrの値に変換され得る。
【００５７】
（データ処理方法）
図６は、本実施形態のデータ処理モードで行われるデータ処理方法のフローを示す図である。
【００５８】
まず、データ取得部２６は、計測装置２４から送られたコーナリング時のタイヤの横力及びセルフアライニングトルクの時系列データを計測データとして取得する（ステップＳ１０）。その際、データ取得部２６は、計測装置２４でタイヤ測定条件として用いられたスリップ角の時系列データα（ｔ）を同時に取得する。計測装置２４が横力及びセルフアライニングトルクの計測データおよびスリップ角の時系列データを装置１０に送ったとき、装置１０は一旦メモリ１６に記憶したのち、データ取得部２６がメモリ１６から呼び出して計測データおよびスリップ角の時系列データを取得してもよい。取得したスリップ角の時系列データはタイヤ力学モデル演算部３４に送られる。一方、取得した計測データは、パラメータ値決定部３０に送られる。
【００５９】
次に、パラメータ値決定部３０は、過渡応答パラメータ（Ｋy，Ｔ0，Ｔ1，Ｔ2，Ａs，Ｋt，Ｔ3）の値を初期設定する（ステップＳ２０）。値の初期設定では、予め定められたデフォルト設定値が用いられてもよいし、以前データ処理モードが行われて決定された過渡応答パラメータの値を用いてもよい。
パラメータ値決定部３０は、設定された過渡応答パラメータ（Ｋy，Ｔ0，Ｔ1，Ｔ2，Ａs，Ｋt，Ｔ3）の値をタイヤ力学モデル演算部３４に送る。
【００６０】
タイヤ力学モデル演算部３４は、パラメータ値決定部３０において設定された過渡応答パラメータ（Ｋy，Ｔ0，Ｔ1，Ｔ2，Ａs，Ｋt，Ｔ3）の値を、式（１５）、式（１６）で示されるステップ応答関数Ｋst（ｔ），Ａst（ｔ）の各過渡応答パラメータに付与する（ステップＳ３０）。これにより、ステップ応答関数は数値計算可能な式となる。
【００６１】
この後、タイヤ力学モデル演算部３４は、数値計算可能になったステップ応答関数Ｋst（ｔ），Ａst（ｔ）と、データ取得部２６から送られたスリップ角の時系列データα（ｔ）とを用いて、横力及びセルフアライニングトルクの時系列データＦy（ｔ），Ｍz（ｔ）である予測データを算出する（ステップＳ４０）。具体的には、式（１７）及び式（１８）に従って、スリップ角の時系列データα（ｔ）の時間変化量ｄα（ｔ’）／ｄｔ’（ｔ’は積分変数）と、横力及びセルフアライニングトルクのステップ応答関数Ｋst（ｔ），Ａst（ｔ）のそれぞれとが乗算されて積分される。すなわち、時間変化量ｄα（ｔ’）／ｄｔ’と、横力及びセルフアライニングトルクのステップ応答関数Ｋst（ｔ），Ａst（ｔ）との畳み込み積分が行われる。算出された横力及びセルフアライニングトルクの時系列データＦy（ｔ），Ｍz（ｔ）は、予測データとしてパラメータ値決定部３０に送られる。
【００６２】
パラメータ値決定部３０は、タイヤ力学モデル演算部３４から送られた横力及びセルフアライニングトルクの予測データが、データ取得部２６から送られた計測データと、許容範囲内で一致するか否かを判定する（ステップＳ５０）。具体的には、予測データと計測データの差分の二乗和を算出し、この二乗和が設定された閾値以下であるか否かを判定する。判定の結果、予測データが計測データに一致しないと判定した場合（ステップＳ５０においてＮＯの場合）、パラメータ値決定部３０は、設定された過渡応答パラメータ（Ｋy，Ｔ0，Ｔ1，Ｔ2，Ａs，Ｋt，Ｔ3）の値を修正し（ステップＳ６０）、修正した値をパラメータ値決定部３０は再度タイヤ力学モデル演算部３４に送り、ステップＳ３０に戻る。こうして、予測データが計測データと許容範囲内で一致するまで、パラメータ値決定部３０における過渡応答パラメータの値の修正と、タイヤ力学モデル演算部３４における予測データの算出とが繰り返し行われる。
【００６３】
ステップＳ５０における判定の結果、予測データが計測データと許容範囲内で一致する場合（ステップＳ５０においてＹＥＳの場合）、パラメータ値決定部３０は、現在設定されている過渡応答パラメータ（Ｋy，Ｔ0，Ｔ1，Ｔ2，Ａs，Ｋt，Ｔ3）の値を、計測データを得たタイヤの過渡応答を再現する過渡応答パラメータの値として決定する（ステップＳ７０）。
さらに、パラメータ値決定部３０は、決定された過渡応答パラメータ（Ｋy，Ｔ0，Ｔ1，Ｔ2，Ａs，Ｋt，Ｔ3）の値を用いて、タイヤ力学要素パラメータの値を算出する（ステップＳ８０）。具体的には、パラメータ値決定部３０は、Ｔ1＝Ｋy／（Ｋdr  ・ｖ）及びＴ2＝Ａs／（Ｋrr  ・ｖ）の関係に従って、ベルト部材の横曲げ剛性Ｋdrの値、サイドウォール部材のねじり剛性Ｋrrの値を算出する。さらに、パラメータ値決定部３０は、横ばね定数ＫLr  ＝  （１／Ｋdr  ＋  （１／Ｋrr  ）・（Ａs／Ｋy）)−１の式に従って横ばね定数ＫLrの値を算出する。ここで、ｖは現在予測しようとするタイヤの回転速度であり、Ｌはタイヤの周長であり、これらの値は、予め入力された値であり既知である。こうして決定された過渡応答パラメータ（Ｋy，Ｔ0，Ｔ1，Ｔ2，Ａs，Ｋt，Ｔ3）の値及び、ベルト部材の横曲げ剛性Ｋdrの値、サイドウォール部材のねじり剛性Ｋrrの値、及び横ばね定数ＫLrの値がディスプレイ２２あるいは図示されないプリンタに出力される。さらに、決定された過渡応答パラメータの値、及びベルト部材の横曲げ剛性Ｋdrの値、サイドウォール部材のねじり剛性Ｋrrの値、及び横ばね定数ＫLrの値がメモリ１６に記憶される。
【００６４】
（タイヤの過渡応答データの算出方法）
図７は、本実施形態の過渡応答算出モードで行われる算出方法のフローを示す図である。過渡応答算出モードでは、スリップ角の時系列データが自在に設定されたときの過渡応答データを算出する。
【００６５】
まず、データ取得部２６は、過渡応答パラメータ（Ｋy，Ｔ0，Ｔ1，Ｔ2，Ａs，Ｋt，Ｔ3）の値を取得する（ステップＳ９０）。取得する過渡応答パラメータ（Ｋy，Ｔ0，Ｔ1，Ｔ2，Ａs，Ｋt，Ｔ3）の値は、入力操作系２０を通してオペレータが入力したものであってもよいし、以前データ処理モードを行ってメモリ１６に記憶された値が用いられてもよい。このときデータ取得部２６は、オペレータが入力操作系２０を通して入力したスリップ角の時系列データ、あるいはメモリ１６に記憶されているスリップ角の時系列データを取得する。取得したスリップ角の時系列データはタイヤ力学モデル演算部３４に送られる。一方、取得した過渡応答パラメータ（Ｋy，Ｔ0，Ｔ1，Ｔ2，Ａs，Ｋt，Ｔ3）の値は、パラメータ値決定部３０に送られる。
パラメータ値決定部３０は、データ取得部２６から送られた過渡応答パラメータ（Ｋy，Ｔ0，Ｔ1，Ｔ2，Ａs，Ｋt，Ｔ3）の値をタイヤ力学モデル演算部３４に送り、予測データを算出するようにタイヤ力学モデル演算部３４に指示をする。
【００６６】
タイヤ力学モデル演算部３４は、パラメータ値決定部３０から送られた過渡応答パラメータ（Ｋy，Ｔ0，Ｔ1，Ｔ2，Ａs，Ｋt，Ｔ3）の値を、式（１５）、式（１６）で示されるステップ応答関数Ｋst（ｔ），Ａst（ｔ）の各過渡応答パラメータに付与する（ステップＳ１００）。これにより、ステップ応答関数Ｋst（ｔ），Ａst（ｔ）は数値計算可能な式となる。
この後、タイヤ力学モデル演算部３４は、数値計算可能になったステップ応答関数Ｋst（ｔ），Ａst（ｔ）と、データ取得部２６から送られたスリップ角の時系列データα（ｔ）とを用いて、横力及びセルフアライニングトルクの時系列データＦy（ｔ），Ｍz（ｔ）である予測データを算出する（ステップＳ１１０）。具体的には、式（１７）及び式（１８）に従って、スリップ角の時系列データα（ｔ）の時間変化量ｄα（ｔ’）／ｄｔ’（ｔ’は被積分変数）と、横力及びセルフアライニングトルクのステップ応答関数とが乗算されて積分される。すなわち、時間変化量ｄα（ｔ’）／ｄｔ’と、Ｋst（ｔ）及びＡst（ｔ）のそれぞれとの畳み込み積分が行われる。横力及びセルフアライニングトルクの時系列データＦy（ｔ），Ｍz（ｔ）は、予測データとしてＦy，Ｍz算出部３２に送られる。
Ｆy，Ｍz算出部３２は、タイヤ力学モデル演算部３４で算出された横力及びセルフアライニングトルクの時系列データＦy（ｔ），Ｍz（ｔ）、すなわち予測データを、所定の形式に修正してディプレイ２２あるいは図示されないプリンタに出力する。
【００６７】
このように本実施形態に用いるステップ応答関数では、式（１５）に示すように、（ｔ−Ｔ1）の時間の項が用いられており、すなわち、タイヤのベルト部材の曲げ変形に対応する遅れ応答を、タイヤ力学モデルにおいてむだ時間Ｔ1による遅れ応答として構成しているので、２次の遅れ応答を用いてステップ応答関数を表すことができる。しかも、ステップ応答関数に用いる過渡応答パラメータ（Ｋy，Ｔ0，Ｔ1，Ｔ2，Ａs，Ｋt，Ｔ3）は、タイヤ力学要素パラメータ、具体的にはベルト部材の横曲げ剛性Ｋdr、サイドウォール部材のねじり剛性Ｋrr、及び横ばね定数ＫLrと式によって関係付けられているので、タイヤ設計者は、タイヤ力学要素パラメータを用いてタイヤの評価を行うとともに、コーナリング特性の改良のための対策方法を容易に見出すことができる。例えば、ベルト部材の横曲げ剛性Ｋdrを増大させる場合、ベルト部材の曲げ剛性を高めるベルト構造を選択する等の対策を講じることができる。
【００６８】
図８は、スリップ角の時系列データα（ｔ）をタイヤ測定条件として計測装置２４に与えて、計測装置２４でタイヤを計測した横力の計測データの一例を示す図である。スリップ角以外のタイヤ測定条件は、タイヤサイズ２０５／５５Ｒ１６であり、内圧２００ｋＰａ、荷重４ｋＮ，タイヤの回転速度ｖ＝４０ｋｍ／時である。
このとき、図９は、パラメータ値決定部３０で決定された過渡応答パラメータＫy，Ａsの値と、ベルト部材の横曲げ剛性Ｋdr、サイドウォール部材のねじり剛性Ｋrr、及び横ばね定数ＫLrの値の一例を示している。また、図９は、横軸にスリップ角を、縦軸に横力をとったときの計測データと予測データの軌跡を示している。図９に示される計測データと予測データの比較から明らかなように、予測データは極めてよく計測データを再現していることがわかる。
【００６９】
図１０（ａ）は、横軸に距離周波数を、縦軸にスリップ角１度における横力のゲイン（Ｆy（１）ゲイン）を示し、図１０（ｂ）は、横軸に距離周波数を、縦軸にスリップ角１度における横力の位相（Ｆy（１）位相）を示している。
距離周波数は、スリップ角αの角周波数ω［ｒａｄ／秒］を、そのときのタイヤの回転速度ｖ［ｍ／秒］で除算した値である。図１０（ａ），（ｂ）に示されるように、予測データは計測データによく一致していることがわかる。従来の一次遅れ応答の応答関数を用いて得られるＦy（１）位相は、図１０（ｂ）に一点鎖線で示すように、距離周波数が高くなると、計測データから乖離するが、本実施形態のステップ応答関数を用いて得られるＦy（１）位相は、距離周波数が高くなっても、計測データとよく一致することがわかる。
【００７０】
図１１は、図９に示す例に用いたタイヤの横ばね定数ＫLrの回転速度ｖの依存性を示す図である。回転速度ｖが大きくなるにつれ、横ばね定数ＫLrの値は減少することが示されている。
【００７１】
図１２（ａ）〜（ｅ）は、異なるタイヤＡ〜Ｅの横ばね定数ＫLr、ベルト部材の横曲げ剛性Ｋdr、及びサイドウォール部材のねじり剛性Ｋrrと、Ｆy（１）ゲイン及びＦy（１）位相を示している。タイヤＡ〜Ｅによって、横ばね定数ＫLr、横曲げ剛性Ｋdr、ねじり剛性Ｋrrが種々異なっており、Ｆy（１）ゲイン及びＦy（１）位相も種々異なっている。例えば、タイヤの過渡応答特性においてＦy（１）ゲインを評価対象とし、Ｆy（１）ゲインを大きくしようとする場合、Ｆy（１）ゲインの大きいタイヤＢ，Ｄ，Ｅに注目する。このタイヤＢ，Ｄ，Ｅは、タイヤＡ，Ｃに対してねじり剛性Ｋrrが大きいので、Ｆy（１）ゲインを大きくするためには、ねじり剛性Ｋrr、すなわち、サイド部材のねじり剛性を高くすればよい、と想定される。したがってタイヤ設計者は、サイド部材のねじり剛性を高めるために、サイド部材の選択あるいはサイド部の構造を再検討する。このように、本実施形態を用いていて算出される横ばね定数ＫLr、ベルト部材の横曲げ剛性Ｋdr、及びサイドウォール部材のねじり剛性Ｋrrの値は、タイヤ設計者にとって有効な情報といえる。
【００７２】
なお、本実施形態では、予測データを式（１７）、式（１８）に示すようにステップ応答関数とスリップ角の時系列データの時間変化量との畳み込み積分を行って算出するが、予測データの算出方法は、この方法に制限されない。例えば、ステップ応答関数の代わりにインパルス応答関数を用いてスリップ角の時系列データとの畳み込み積分を行ってもよい。
【００７３】
本実施形態では、過渡応答算出モードでは、予測データとしてＦy（ｔ），Ｍz（ｔ）を算出するが、Ｆy（ｔ），Ｍz（ｔ）の少なくとも一方を算出してもよい。一方、データ処理モードでは、横力及びセルフアライニングトルクの緩和応答パラメータのうち、少なくとも横力の緩和応答パラメータであるＫy，Ｔ0，Ｔ1，Ｔ2の値を算出すればよい。しかし、横ばね定数ＫLr、ベルト部材の横曲げ剛性Ｋdr、及びサイドウォール部材のねじり剛性Ｋrr等のタイヤ力学要素パラメータの値を算出する点から、横力及びセルフアライニングトルクの緩和応答パラメータ（Ｋy，Ｔ0，Ｔ1，Ｔ2，Ａs，Ｋt，Ｔ3）の値を算出することが好ましい。
【００７４】
以上、本発明のタイヤの過渡応答データの算出方法および算出装置、データ処理方法およびデータ処理装置、及びプログラムについて詳細に説明したが、本発明は上記実施形態に限定されず、本発明の主旨を逸脱しない範囲において、種々の改良や変更をしてもよいのはもちろんである。
【符号の説明】
【００７５】
１０  装置
１２  ＣＰＵ
１４  ＲＯＭ
１６  メモリ
１８  入出力部
２０  入力操作系
２２  ディスプレイ
２４  計測装置
２６  データ取得部
３０  パラメータ値決定部
３２  Ｆy，Ｍz算出部
３４  タイヤ力学モデル演算部
【図１】
0005790018.tif 000040
【図２】
0005790018.tif 000041
【図３】
0005790018.tif 000042
【図４】
0005790018.tif 000043
【図５】
0005790018.tif 000044
【図６】
0005790018.tif 000045
【図７】
0005790018.tif 000046
【図８】
0005790018.tif 000047
【図９】
0005790018.tif 000048
【図１０】
0005790018.tif 000049
【図１１】
0005790018.tif 000050
【図１２】
0005790018.tif 000051