JP4291561B2

JP4291561B2 - Method and apparatus for predicting physical quantity related to tire wear, and computer program

Info

Publication number: JP4291561B2
Application number: JP2002308985A
Authority: JP
Inventors: 全一郎信田
Original assignee: Yokohama Rubber Co Ltd
Current assignee: Yokohama Rubber Co Ltd
Priority date: 2002-10-23
Filing date: 2002-10-23
Publication date: 2009-07-08
Anticipated expiration: 2022-10-23
Also published as: JP2004142571A

Description

【０００１】
【発明の属する技術分野】
この発明は、タイヤの摩耗に関する物理量の予測方法に関し、さらに詳しくは、タイヤのせん断接触力や摩耗量の分布をより正確に予測できるタイヤの摩耗に関する物理量の予測方法及び予測装置、並びにコンピュータプログラムに関する。
【０００２】
【従来の技術】
タイヤの摩耗を予測する方法としては、例えば、特許文献１、２に開示されている方法が知られている。これらに開示されている予測方法は、次の手順から構成される。すなわち、（１）タイヤ表面のゴムの摩耗指数を測定する手順、（２）種々の運転条件下におけるタイヤの摩擦エネルギーを測定する手順、（３）各運転条件の頻度に応じて各運転条件の重み付けをする手順を経て、タイヤの摩耗を予測する。
【０００３】
【特許文献１】
特開平１１−３２６１４３号公報（第２頁〜２３頁、図１〜１２）
【特許文献２】
特開２００１−１７２３号公報（第２頁〜９頁、図１〜１０）
【０００４】
【発明が解決しようとする課題】
このように、従来のタイヤの摩耗予測方法では、摩擦エネルギーを測定することによってタイヤの摩耗を予測するので、摩擦エネルギーの測定に手間を要していた。また、摩擦エネルギーは点で測定されるので、タイヤ表面の単一ブロック内における偏摩耗を含め、トレッド面すべてを評価するためには、膨大な測定が必要である。さらに、ブロックやリブの端部、あるいは細溝が密集した部分のように測定ポイントが極めて小さい部分においては、摩擦エネルギーを測定することが極めて困難である。その結果、タイヤの摩耗に対する予測精度を十分に高くできず、予測結果を設計に十分反映することができなかった。
【０００５】
そこで、この発明は、上記に鑑みてなされたものであって、摩擦エネルギーの測定を省略することによって、効率よくタイヤの摩耗を予測すること、より精度よくタイヤの摩耗を予測することのうち少なくとも一つを達成できるタイヤの摩耗に関する物理量の予測方法及び予測装置、並びにコンピュータプログラムを提供することを目的とする。
【０００６】
【課題を解決するための手段】
本発明者らは摩擦エネルギーの測定を省略するために、摩擦エネルギーを有限要素法、境界要素法その他の解析手法によって予測する方法を用いた。このような解析手法を用いたシミュレーションによって摩擦エネルギーを求めたところ、シミュレーションと実測値とではタイヤの偏摩耗や摩耗速度の評価結果が異なることが判明した。本発明者らは鋭意研究の結果、この原因を見出した。
【０００７】
有限要素法等の解析手法においては、タイヤのブロック等の予測領域を複数の要素に分割して当該要素内は均一とみなして取り扱う。そして、要素を構成するそれぞれの節点についてせん断接触力とすべり量とを求め、両者から摩擦エネルギーを算出する。しかしながら、一定の大きさを持つ予測領域を複数の要素で分割し、複数の節点が接地する条件では、各節点の分担する垂直接触力、せん断接触力及び摩擦エネルギーは要素の分割数や節点の位置によって異なる。このように、各節点における摩擦エネルギーの大きさが、要素の分割数や節点位置によって影響を受けるので、シミュレーションと実測値とにおける摩耗の評価結果が異なるのである。そこで、要素の分割数や節点位置の影響を排除するため、単位面積当たりの垂直接触力や摩擦エネルギーによって摩耗を評価することとして、上記問題を解決した。
【０００８】
本発明に係るタイヤの摩耗に関する物理量の予測方法は、タイヤの摩耗を有限要素法、境界要素法その他の解析手法によって予測するにあたり、タイヤを複数の分割要素に分割してタイヤモデルを作成する手順と、前記分割要素の表面に存在する節点が接地しているか否かを判定する手順と、を含み、前記節点が接地している場合には、接地していない前記分割要素の表面の面積は除外して、前記節点が分担する前記分割要素の表面の面積を算出する手順と、算出した前記節点が分担する面積によってタイヤの摩耗に関する物理量を単位面積当たりの物理量に換算する手順と、を実行し、前記節点が接地していない場合には、全節点に対する処理が終了したか否かを判定する手順を実行することを特徴とする。
【０００９】
また、本発明に係るタイヤの摩耗に関する物理量の予測プログラムは、タイヤの摩耗に関する物理量を有限要素法、境界要素法その他の解析手法によって予測するにあたり、タイヤを複数の要素に分割してタイヤモデルを作成する手順と、前記分割要素の表面を構成する節点が接地しているか否かを判定する手順と、を含み、前記節点が接地している場合には、接地していない前記分割要素の表面の面積は除外して、前記節点が分担する前記分割要素の表面の面積を算出する手順と、算出した前記節点が分担する面積によってタイヤの摩耗に関する物理量を単位面積当たりの物理量に換算する手順と、を実行し、前記節点が接地していない場合には、全節点に対する処理が終了したか否かを判定する手順を実行することをコンピュータに実行させることを特徴とする。
【００１０】
このタイヤの摩耗に関する物理量の予測方法又は予測プログラムは、タイヤモデルを構成する複数の分割要素の表面に含まれる節点が分担する分割要素の表面の面積によって、タイヤの摩耗に関する物理量を単位面積当たりの物理量に換算してタイヤの摩耗を評価する。このため、せん断接触力や摩擦エネルギー等の物理量は要素分割や節点位置の影響を受けないので、接触面全体における摩耗を精度よく予測することができる。
【００１１】
また、本発明に係るタイヤの摩耗に関する物理量の予測方法の各手順をコンピュータに実行させるタイヤの摩耗に関する物理量の予測プログラムによれば、コンピュータ上でタイヤの摩耗を予測できる。これにより、摩擦エネルギーを測定する必要がないので、効率よくタイヤの摩耗を予測して、開発の効率を向上させることができる。また、ブロックの端部や細溝が密集した部分のように測定ポイントが極めて小さく、摩擦エネルギーを測定することが極めて困難な部分でも容易且つ精度よく摩擦エネルギーを予測することができる。これにより、タイヤの摩耗に対する予測精度を十分に高くできるので、当該予測結果を設計に十分に活用させることができる。
【００１２】
ここで、タイヤの摩耗に関する物理量とは、タイヤ接地面に作用する垂直接触力、垂直接触力に摩擦係数を乗じて算出される最大摩擦力、せん断接触力、摩擦エネルギー、摩耗量、すべり量、摩耗量その他の摩耗に関連する物理量をいう（以下同様）。このうち、本発明においては、タイヤ接地面に作用する垂直接触力、垂直接触力に摩擦係数を乗じて算出される最大摩擦力、せん断接触力、摩擦エネルギーその他の単位面積当たりの物理量に換算できるものについては、単位面積当たりの物理量に換算して取り扱う（以下同様）。
【００１３】
また、タイヤの摩耗を評価する際には、ある瞬間における単位面積当たりの摩擦エネルギーを、タイヤ表面が接地を開始してから止まるまでの全時間について加算したものを用いる必要がある。この場合には、上記発明の有する工程に、求めた単位面積当たりの摩擦エネルギーを、タイヤ表面が接地を開始してからの単位面積当たりの総摩擦エネルギーに順次追加する工程を加えて、タイヤ表面が接地していた全時間における単位面積当たりの摩擦エネルギーを算出する。
【００１４】
また、本発明に係るタイヤの摩耗に関する物理量の予測方法は、上記節点が分担する上記分割要素の表面の面積を算出する手順においては、上記節点を含む上記分割要素の表面の面積を算出する手順と、上記節点を含む上記分割要素の表面の分割形状に基づいて、上記節点が上記分割要素の表面の面積を分担する割合を決定する手順と、上記節点が上記分割要素の表面の面積を分担する面積を、上記節点を含むすべての分割要素の表面について加算する手順と、を含むことを特徴とする。
【００１５】
また、本発明に係るタイヤの摩耗に関する物理量の予測プログラムは、上記タイヤの摩耗に関する物理量の予測プログラムにおいて、上記節点が分担する上記分割要素の表面の面積を算出する手順においては、上記節点を含む上記分割要素の表面の面積を算出する手順と、上記節点を含む上記分割要素の表面の分割形状に基づいて、上記節点が上記分割要素の表面の面積を分担する割合を決定する手順と、上記節点が上記分割要素の表面の面積を分担する面積を、上記節点を含むすべての分割要素の表面について加算する手順と、をコンピュータに実行させることを特徴とする。
【００１６】
この発明では、分割要素の表面の分割形状に基づいて節点が分担する面積を求める。そして、タイヤ接地面に作用する垂直接触力、垂直接触力に摩擦係数を乗じて算出される最大摩擦力、せん断接触力、摩擦エネルギーその他の単位面積当たりの物理量に換算できる物理量を、節点が分担する面積を用いて単位面積当たりの物理量に換算する。これによって、分割要素の表面の形状に関わらず節点が分担する面積を正確に算出できるので、タイヤの摩耗に関する物理量を単位面積当たりの物理量に正しく換算できる。また、本発明の手順をコンピュータに実行させるためのタイヤの摩耗に関する物理量の予測プログラムにおいても、同様の作用・効果を奏する。
【００１７】
また、本発明に係るタイヤの摩耗に関する物理量の予測方法は、上記タイヤの摩耗に関する物理量の予測方法において、上記タイヤモデルの転動時における変形形状を基に、前記分割要素の表面に存在する節点が接地しているか否かを各時間ステップ毎に判定することを特徴とする。
【００１８】
また、本発明に係るタイヤの摩耗に関する物理量の予測プログラムは、上記タイヤの摩耗に関する物理量の予測プログラムにおいて、上記タイヤモデルの転動時における変形形状を基に、前記分割要素の表面に存在する節点が接地しているか否かを各時間ステップ毎に判定することをコンピュータに実行させることを特徴とする。
【００１９】
本発明は、上記発明において分割要素の表面が接地しているか否かを判定する工程を含むものであるので、上記発明と同様の作用・効果を奏する。さらに、分割要素の表面が接地しているか否かを判定する工程を含むので、例えば溝壁等における分割要素の表面のように、当該分割要素の表面に含まれる節点は接地しているが分割要素の表面自体は接地していない分割要素の表面を、物理量の予測に際して排除することができる。これによって、溝に接するブロック端部における摩擦エネルギーを適切に評価できるので、このような部分における摩耗も精度よく予想できる。
【００２０】
また、ブロックの端部は測定ポイントが極めて小さいので、従来摩擦エネルギーを測定することは極めて困難であったが、このような部分でも分割要素の表面の接地状態を判定することによって、精度よく摩擦エネルギーを予測することができる。これにより、タイヤの摩耗に対する予測精度を十分に高くできるので、当該予測結果を設計に十分に活用させることができる。また、本発明の手順をコンピュータに実行させるためのタイヤの摩耗に関する物理量の予測プログラムにおいても、同様の作用・効果を奏する。
【００２１】
また、本発明に係るタイヤの摩耗に関する物理量の予測方法は、上記タイヤの摩耗に関する物理量の予測方法において、上記分割要素の表面に存在する節点が接地しているか否かを判定する手順においては、上記タイヤの初期形状を基準として、上記分割要素の表面が接地しているか否かを判定することを特徴とする。
【００２２】
この発明では、タイヤの初期形状を基準として、節点が分担する面積を算出する前に予め分割要素の表面が接地しているか否かを判定する。このように、接地するか否かは、タイヤの初期形状を基にして１回だけ判定するので、計算数を抑えることができる。特に、コンピュータを用いて本発明に係る予測方法を実行する場合には、ハードウェア資源を節約できるという長所がある。なお、具体的な判定方法としては、例えば、上記分割要素の表面の法線ベクトルと、上記タイヤの中心から上記分割要素の表面の中心へ向かうベクトルとのなす角度が所定の角度よりも小さい場合には、当該分割要素の表面は接地すると判定するものがある。
【００２３】
また、本発明に係るタイヤの摩耗に関する物理量の予測方法は、上記タイヤの摩耗に関する物理量の予測方法において、さらに、上記タイヤが変形した形状に応じて、上記分割要素の表面が接地しているか否かを判定することを特徴とする。この発明では、タイヤの変形形状に応じて、分割要素の表面が接地しているか否かを各瞬間毎に判定するので、よりタイヤの摩耗に関する物理量の予測精度を向上させることができる。なお、この発明をコンピュータで実行させる場合には、分割要素の表面が接地しているか否かは各時間ステップ毎に判定される。
【００２４】
接地の判定手法には、例えば、判定対象である分割要素の表面の法線ベクトルと地面の法線ベクトルとのなす角度が所定の角度以上になったときには、当該分割要素の表面は接地すると判定するものがある。また、他の判定手法としては、判定対象である分割要素の表面の鉛直方向における座標と地面の鉛直方向における座標とを比較して、両者が一致すれば当該分割要素の表面が接地していると判定するものがある。さらに、上述したように、判定対象である分割要素の表面の法線ベクトルと、上記タイヤの中心から上記分割要素の表面の中心へ向かうベクトルとのなす角度が所定の角度よりも小さい場合には、当該分割要素の表面は接地すると判定してもよい。
【００２５】
また、本発明に係るタイヤの摩耗に関する物理量の予測方法は、自由転動、制動、駆動、旋回その他の運転条件下において、前記タイヤの摩耗に関する物理量の予測方法によって、それぞれの前記各運転条件下における単位面積当たりの摩擦エネルギーを予測する工程と、前記各運転条件の発生頻度を考慮して、前記各運転条件に対して重み付けをする工程と、前記各運転条件下における単位面積当たりの摩擦エネルギーと前記各運転条件の重み付けとから、実車走行中における単位面積当たりの摩擦エネルギーを予測する工程と、を含むことを特徴とする。
【００２６】
また、本発明に係るタイヤの摩耗に関する物理量の予測プログラムは、自由転動、制動、駆動、旋回その他の運転条件下において、前記タイヤの摩耗に関する物理量の予測方法によって、それぞれの前記各運転条件下における単位面積当たりの摩擦エネルギーを予測する手順と、前記各運転条件の発生頻度を考慮して、前記各運転条件それぞれの重み付けをする手順と、前記各運転条件下における単位面積当たりの摩擦エネルギーと前記各運転条件の重み付けから、実車走行中における単位面積当たりの摩擦エネルギーを予測する手順と、をコンピュータに実行させることを特徴とする。
【００２７】
この発明では、単位面積当たりの摩擦エネルギーを用いてタイヤの摩耗を予測するので、精度の高い摩耗予測ができる。これによって、タイヤの設計に対する有効な情報を得ることができる。また、摩擦エネルギーを測定する必要がないので、タイヤの摩耗を予測する手間も大幅に軽減でき、タイヤの開発期間を短縮できる。なお、本発明の手順をコンピュータに実行させるためのタイヤの摩耗に関する物理量の予測プログラムにおいても、同様の作用・効果を奏する。
【００２８】
また、本発明に係るタイヤの摩耗に関する物理量の予測方法のように、上記タイヤの摩耗に関する物理量の予測方法において、さらに、タイヤの表面を構成するゴムのゴム摩耗指数を測定する手順を含み、当該ゴム摩耗指数を実車走行中における単位面積当たりの摩擦エネルギーに乗じてタイヤの摩耗を予測してもよい。このように、ゴム摩擦指数を使用すれば、単なる摩耗の分布のみならず、実際の摩耗速度も予測できるので、より幅広くタイヤの摩耗に関する評価ができる。
【００２９】
また、本発明に係るタイヤの摩耗に関する物理量の予測装置は、上記タイヤの摩耗に関する物理量の予測方法における各手順を処理する処理手段と、この処理手段に前記タイヤを構成する材料の物性値、運転条件、境界条件その他のデータを与える入力手段と、前記処理手段による予測結果を表示する表示手段と、を備えたことを特徴とする。
【００３０】
このように、本発明に係るタイヤの摩耗に関する物理量の予測装置では、上記タイヤの摩耗に関する物理量の予測方法を実行させるための処理手段を備えているので、せん断接触力や摩擦エネルギー等の物理量が要素分割や節点位置の影響を受けずに、接触面全体における摩耗を精度よく予測することができる。また、摩擦エネルギーを測定する必要がないので、効率よくタイヤの摩耗を予測して、開発の効率を向上させることができる。さらに、ブロックの端部や細溝が密集した部分のように測定ポイントが極めて小さく、摩擦エネルギーを測定することが極めて困難な部分でも容易且つ精度よく摩擦エネルギーを予測することができる。これにより、タイヤの摩耗に対する予測精度を十分に高くできるので、当該予測結果を設計に十分に反映させることができる。
【００３１】
【発明の実施の形態】
以下、この発明につき図面を参照しつつ詳細に説明する。なお、この実施の形態によりこの発明が限定されるものではない。また、下記実施の形態における構成要素には、当業者が容易に想定できるものあるいは実質的に同一のものが含まれる。
【００３２】
（実施の形態１）
図１は、予測対象であるタイヤを、その回転軸を含む子午面で切った断面を示す一部断面図である。同図を用いて、本実施の形態における予測対象であるタイヤ１０の構造について簡単に説明する。キャップトレッド１１は、タイヤ１０の路面接地部に配置されており、カーカス１５、ベルト１４又はブレーカの外側を覆うゴム層である。キャップトレッド１１は、カット衝撃に対してカーカス１５やベルト１４を保護する役目を持っている。また、雨天走行時の排水性を高めるために、キャップトレッド１１には溝１８が形成されている。この溝１８で仕切られる部分がタイヤ１０のブロックとなる。
【００３３】
アンダトレッド１２は、キャップトレッド１１とベルト１４との間に配置されるゴム層で、発熱性、接着性等を向上させる目的で用いられる。サイドトレッド１３は、サイドウォール部の最も外側に配置されて外からの傷がカーカス１５に達するのを防止するとともに、ラジアルタイヤの場合には、車軸からの駆動力を路面に伝える補助的役割も担っている。
【００３４】
ベルト１４は、キャップトレッド１１とカーカス１５との間に配置されたゴム引きコード層である。なお、バイアスタイヤの場合にはブレーカと呼ぶ。ラジアルタイヤにおいて、ベルト１４はタイヤ１０の形状保持及び強度メンバーとして重要な役割を担っている。カーカス１５はタイヤ１０の骨格をなすゴム引きコード層である。カーカス１５は、タイヤ１０に空気や窒素を充填した際に圧力容器としての役目を果たす強度メンバーであり、その内圧によって荷重を支え、走行中の動的荷重に耐える構造としてある。
【００３５】
ビード１６は、内圧によって発生するカーカス１５のコード張力を支えているスチールワイヤの束を、硬質ゴムで固めたリングである。ビード１６は、タイヤ１０をホイールのリムに固定させる役割を果たす他、カーカス１５、ベルト１４及びトレッドとともに、タイヤ１０の強度部材となる。ビードフィラ１７は、カーカス１５をビードワイヤの周囲に巻き込む際に生ずる空間へ充填するゴムである。ビードフィラ１７は、カーカス１５をビード１６に固定するとともに、その部分の形状を整え、ビード部全体の剛性を高める。
【００３６】
ホイール（図示せず）からの回転力、すなわち、車体を動かす駆動力や制動力、あるいは旋回力は、ビード１６を介してタイヤ１０へ伝えられる。このような車体を動かす駆動力や制動力、あるいは旋回力を路面に伝達する際には、キャップトレッド１１と路面との間にせん断接触力が発生する。同時に、このせん断接触力がタイヤ１０の接地面１１ｐと路面とのすべりを発生させ、タイヤ１０の接地面１１ｐを摩耗させる。次に、本発明に係るタイヤの摩耗に関する物理量の予測方法について説明する。ここでは、タイヤの摩耗に関する物理量として摩擦エネルギーを予測するが、次に説明する予測手順はせん断接触力や垂直接触力等の他の物理量を予測する場合も同様に適用できる。
【００３７】
本発明においては、タイヤの摩耗に関する物理量を予測する解析手法としてＦＥＭ（Finite Element Method：有限要素法）を使用する。なお、本発明に係るタイヤの摩耗に関する物理量の予測方法に適用できる解析手法はＦＥＭに限られず、ＢＥＭ（Boundary Element Method：境界要素法）、ＦＤＭ（Finite Differences Method：有限差分法）等も使用できる。タイヤの種類や境界条件等によって最も適当な解析手法を選択し、又は複数の解析手法を組み合わせて使用することが好ましい。
【００３８】
図２は、予測対象であるタイヤを微小要素に分割した例を示す斜視図である。同図に示すように、予測対象であるタイヤ１０は、有限要素法に基づき、有限個の微小要素１０ａ、１０ｂ等に分割される。有限要素法に基づく微小要素とは、例えば２次元平面においては三角形要素、四角形要素、３次元体としては四面体ソリッド要素、五面体ソリッド要素、六面体ソリッド要素等のソリッド要素や、三角形シェル要素、四角形シェル要素等のシェル要素等、コンピュータで用いうる要素とすることが望ましい。このようにして分割された微小要素は、解析の過程において３次元座標を用いて逐一特定される。
【００３９】
本発明に係るタイヤの摩耗に関する物理量の予測方法では、タイヤ１０の接地面１１ｐにおける摩擦エネルギーの分布を正しく評価するために、単位面積当たりの摩擦エネルギーを用いる。次にこの理由について説明する。図３は、キャップトレッドが地面と接する接地面を微小要素に分割した例を示す説明図である。図３（ａ）に示すように、接地面１１ｐを１分割した場合には、それぞれの節点Ｎ１〜Ｎ４には垂直接触力Ｆが均等に作用するので、各節点Ｎ１〜Ｎ４にはそれぞれＦ／４の垂直接触力が作用する。
【００４０】
一般に、ＦＥＭやＢＥＭ等の解析手法においては、一定とみなす要素の分割数を大きくする程、計算時間は多くなるが計算精度は高くなる。したがって、計算時間とのバランスを考慮した上で要素分割の数をできるだけ多くする。また、ブロック１１Ｂの端部においては摩擦エネルギーが高くなるので、要素分割数が少ないとブロック１１Ｂ内部の評価精度が低下する。このため、接地面１１ｐの摩擦エネルギーを予測する場合においては、図３（ｂ）に示すように、接地面１１ｐの分割数を多くする。図３（ｂ）に示した例では接地面１１ｐを４分割しているので、当該接地面１１ｐにはＮ１〜Ｎ９までの計９節点が含まれる。
【００４１】
ここで、接地面１１ｐの角部における節点Ｎ１、Ｎ３、Ｎ７、Ｎ９には、節点数と節点の場所との関係から、それぞれＦ／１６の垂直接触力が作用する。同様に、接地面１１ｐの辺上における節点Ｎ２、Ｎ４、Ｎ６、Ｎ８にはそれぞれＦ／８の垂直接触力が作用し、接地面１１ｐの中心部における節点Ｎ１にはＦ／４の垂直接触力が作用する。ここで、摩擦係数をμ、すべり量をＬとすると、せん断接触力はμ×垂直接触力で、摩擦エネルギーはμ×垂直接触力×Ｌで求めることができる。したがって、節点Ｎ１、Ｎ３、Ｎ７、Ｎ９における摩擦エネルギーはＬ×μ×Ｆ／１６であり、節点Ｎ２、Ｎ４、Ｎ６、Ｎ８における摩擦エネルギーはＬ×μ×Ｆ／８であり、また、節点Ｎ１における摩擦エネルギーはＬ×μ×Ｆ／４である。
【００４２】
本来、この例の接地面１１ｐでは摩擦エネルギーが均一に分布する。しかし、各節点Ｎ１〜Ｎ９のせん断接触力及びすべり量から摩擦エネルギーを算出すると、要素分割数や節点数あるいは節点位置によって摩擦エネルギーの値が異なることとなり、接地面１１ｐの摩耗を正しく評価することはできない。そこで、本発明においては、単位面積当たりの摩擦エネルギーによって接地面１１ｐの摩耗を評価することとした。同時に、垂直接触力及びせん断接触力についても、単位面積当たりの垂直接触力（すなわち垂直接地圧Ｐｎ）と、単位面積当たりのせん断接触力（Ｐｓ＝μ×Ｐｎ）として評価することとした。このようにすれば、接地面１１ｐ表面における摩擦エネルギーの分布を正しく求めることができるので、接地面１１ｐの摩耗を正しく評価することができる。
【００４３】
単位面積当たりにおける摩擦エネルギー等の物理量によって評価する場合には、各節点が分担する分割要素の表面（以下分割要素表面）の面積を求める必要がある。次にこの手順の一例について説明する。図４は、四角形要素によって９分割した接地面を示す平面図である。接地面１１ｐは９分割されているので、節点はＮ１〜Ｎ１６の１６節点となる。ここでは、節点Ｎ６の面積を算出する場合を説明する。まず、節点Ｎ６を含む分割要素表面の情報を取得する。この例において、Ｎ６を含む分割要素表面とはＳ１、Ｓ２、Ｓ４、Ｓ５である。また、分割要素表面の情報とは、当該分割要素表面を構成する節点、及び当該節点の座標といった情報である。
【００４４】
次に、各分割要素表面Ｓ１、Ｓ２、Ｓ４、Ｓ５それぞれの面積を算出する。これは、取得した各分割要素表面を構成する節点の座標、及び各分割要素表面の形状から算出することができる。各分割要素表面の面積を算出したら、節点Ｎ６が分担する面積を各分割要素表面の形状に応じて算出する。図４（ａ）に示した例では、接地面１１ｐが四角形要素で分割されているので、各分割要素表面Ｓ１、Ｓ２、Ｓ４、Ｓ５の形状は四角形となる。この場合、すなわち分割要素によって分割された分割要素表面の形状が四角形である場合には、分割要素表面Ｓ１の面積の１／４を節点Ｎ６が分担する。
【００４５】
このように、接地面１１ｐにおける分割要素表面の形状が四角形である場合には、当該形状が占める面積の１／４の大きさを対象となる節点が分担することになる。また、図４（ｂ）に示すように、接地面における分割要素表面の形状が三角形である場合には、当該形状が占める面積の１／３の大きさを対象となる節点が分担することになる。次に、節点Ｎ６が分担する各分割要素表面Ｓ１、Ｓ２、Ｓ４、Ｓ５の面積をすべて加算する。この面積が、節点Ｎ６が分担する面積となる。同様の手順を残りの節点Ｎ１、Ｎ２等に対して適用して、接地面１１ｐ上における各節点が分担する面積を算出する。
【００４６】
図５は、タイヤモデルの接地面を含む分割要素の一つを示す斜視図である。上記手順によって各節点が分担する面積を算出すると、図５（ｃ）に示す溝壁１１ｓのような実際には接地していない分割要素表面の面積（図５（ａ）のハッチング部分）も節点Ｎ１が分担する面積に含めてしまうことになる。これは、実際の接地面積よりも過大な面積によって単位面積当たりの摩擦エネルギーを求めることになるので、実際よりも小さい値の単位面積当たりの摩擦エネルギーによって接地面１１ｐの摩耗を評価することになる。この不都合を回避するために、接地しない分割要素表面を含む節点（例えばＮ１）の面積を算出する場合には、接地する分割要素表面２０と接地しない分割要素表面２１とを判定して、接地していない分割要素表面２１の面積は除外して各節点が分担する面積を求める必要がある（図５（ｂ））。次に、この判定方法について説明する。
【００４７】
図６は、分割要素表面の接地を判定する方法の一例を示す説明図である。また、図７は、分割要素表面の接地を判定する方法の他の例を示す説明図である。分割要素表面の接地を判定する方法には、タイヤの初期形状に基づいて、接地する（あるいは接地する可能性のある）分割要素表面と、接地しない（あるいは接地する可能性の低い）分割要素表面とを節点が分担する面積を算出する前に判定する方法がある。
【００４８】
この判定法には、例えばタイヤ１０の中心１０ｃから判定対象である分割要素表面の中心に向かうベクトルを基準ベクトルとし、分割要素表面の法線ベクトルと前記基準ベクトルとがなす角度から接地の有無を判定するものがある。この場合、両ベクトルのなす角度が接地基準角度よりも小さい場合には当該分割要素表面は接地すると判定し、両ベクトルのなす角度が所定値よりも大きい場合には当該分割要素表面は接地しないと判定する。両ベクトルのなす角度は３５度〜４５度を接地基準角度とすれば、実用上十分な精度で接地の有無を判断できるので、接地面１１ｐの摩耗を精度よく評価できる。ここで、接地基準角度とは、当該角度よりも小さい場合には、その分割要素表面は接地しない角度をいうものとする。
【００４９】
例えば、図６に示した例においては、タイヤ１０の中心１０ｃから判定対象である分割要素表面２０の中心に向かうベクトルを基準ベクトルＶ１１とし、分割要素表面２０の法線ベクトルＶ１２と前記基準ベクトルＶ１１とがなす角度θ１とする。接地基準角度を例えば３５度とした場合、分割要素表面２０における角度θ１は接地基準角度よりも小さいので、分割要素表面２０は接地面（１１ｐ）であると判定される。一方、分割要素表面２１は、基準ベクトルがＶ２２であり、法線ベクトルはＶ２１である。両ベクトルのなす角度はθ２であるが、θ２は接地基準角度よりも大きいので、分割要素表面２１は接地しないと判定される。なお、分割要素表面２１は、溝壁１１ｓの一部である。
【００５０】
この判定方法によれば、各節点が分担する面積を求める際に、接地しない分割要素表面を除外できるので、精度良く接地面１１ｐの摩耗を予測することができる。また、タイヤの初期形状に基づいて各分割要素表面の接地状態を判定するので、各時間ステップ毎に各分割要素表面の接地状態を判定する必要はない。これにより、計算量を少なくできるのでコンピュータシミュレーションにおいては計算に要する時間を短くでき、また、ハードウェア資源が限られている場合でも有効な評価ができる。
【００５１】
他の判定方法としては、タイヤ１０の転動時における変形形状を基に、接地する分割要素表面と接地しない分割要素表面とを各時間ステップ毎に判定する方法もある。この判定方法では、タイヤ１０の転動時における変形形状に応じて接地面を判断するので、より高い精度で接地面１１ｐの摩耗を予測することができる。この判定方法においては、基準ベクトルとして上記Ｖ１１やＶ２２を用いてもよいし、図７（ａ）に示すように分割要素表面２０、２１の法線ベクトルＶｎ１、Ｖｎ２と路面３０の法線ベクトルＶｂｒとを用いてもよい。なお、路面３０の法線ベクトルＶｂｒを使用する場合には、接地基準角は１６０度〜１７９度が好ましく、当該接地基準角よりも分割要素表面の法線ベクトルと路面の法線ベクトルとのなす角度が大きいときには接地すると判定される。
【００５２】
例えば、接地基準角を１６０度とした場合には、分割要素表面２０における角度θ１は接地基準角よりも大きいので、当該分割要素表面２０は接地すると判定される。一方、分割要素表面２１における角度θ２は接地基準角よりも小さいので、当該分割要素表面２１は接地しないと判定される。なお、本判定例において、分割要素表面の法線ベクトルと路面の法線ベクトルとのなす角度（図７（ａ）のθ１、θ２）は、最大１８０度である。
【００５３】
また、路面３０の座標Ｒ（Ｘｒ、Ｙｒ、Ｚｒ）と判定対象である分割要素表面２０等の座標Ｔ₁（Ｘｔ１、Ｙｔ１、Ｚｔ１）等との差から、分割要素表面２０の接地状態を判定してもよい。この場合には、鉛直方向の座標（図７中矢印Ｚで示す方向）を用いて分割要素表面２０等の接触状態を判定する。例えば、Ｚｔ１とＺｒとの差δ＝（Ｚｔ１−Ｚｒ）が基準値以下である場合には分割要素表面２０は接地すると判定される。このような判定をするのは、次の理由による。すなわち、コンピュータを用いた数値シミュレーションの場合は数値誤差が避けられないので、δが完全に一致することはありえない。したがって、この数値誤差を考慮して接地するか否かを判定する必要があるからである。ここで、前記基準値は０．０１ｍｍ〜２ｍｍ程度が好ましい。
【００５４】
また、分割要素表面２０の座標には、当該分割要素表面２０の中心における座標Ｎｃを用いることができる。例えば、分割要素表面２０に含まれる節点Ｎ１〜Ｎ４の鉛直方向における座標の平均値（Ｚ１＋Ｚ２＋Ｚ３＋Ｚ４）／４を、当該分割要素表面２０の中心における鉛直方向の座標値Ｚｃとすることができる（図７（ｃ）参照）。そして、路面３０の鉛直方向における座標値Ｚｒと分割要素表面２０の鉛直方向における座標値Ｚｃとを比較して、両者が等しい場合には分割要素表面は接地すると判定される。さらに、分割要素表面２０に含まれる節点Ｎ１〜Ｎ４のうち、１節点の鉛直方向における座標が路面３０の鉛直方向における座標と等しくない場合には、分割要素表面２０は接地しないと判断してもよい。
【００５５】
図８は、この発明の実施の形態１に係るタイヤの摩耗に関する物理量予測装置の一例を示す説明図である。物理量予測装置５０は、処理部５２と記憶部５４とで構成される。また、この物理量予測装置５０には、入出力装置５１が接続されており、ここに備えられた入力手段５３でタイヤ１０の補強コード１やゴムの物性値、あるいは境界条件等を処理部５２や記憶部５４へ入力する。ここで、入力手段５３には、キーボード、マウス等の入力デバイスを使用することができる。
【００５６】
記憶部５４には、本発明に係るタイヤの摩耗に関する物理量予測方法を実現する本発明の予測方法が組み込まれたＦＥＭのプログラムが格納されている。ここで、記憶部５４は、ハードディスク装置や光磁気ディスク装置、又はフラッシュメモリ等の不揮発性のメモリ（ＣＤ−ＲＯＭ等のような読み出しのみが可能な記憶媒体）や、ＲＡＭ（Random Access Memory）のような揮発性のメモリ、あるいはこれらの組み合わせにより構成することができる。
【００５７】
また、上記プログラムは、コンピュータシステムにすでに記録されているプログラムとの組み合わせによって、本発明に係るタイヤの摩耗に関する物理量の予測方法を実現できるものであってもよい。また、図８における処理部５２の機能を実現するための上記プログラムをコンピュータ読み取り可能な記録媒体に記録して、この記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することにより本発明に係る予測方法を実行してもよい。なお、ここでいう「コンピュータシステム」とは、ＯＳや周辺機器などのハードウェアを含むものとする。
【００５８】
処理部５２は、メモリ及びＣＰＵにより構成されている。物理量予測時においては設定されたタイヤモデル及び入力データに基づいて、処理部５２が前記プログラムを当該処理部５２に組み込まれたメモリに読み込んで演算する。その際に処理部５２は、適宜記憶部５４へ演算途中の数値を格納し、また格納した数値を取り出して演算を進める。なお、この処理部５２は、前記プログラムの代わりに専用のハードウェアにより実現されるものであってもよい。予測結果は、入出力装置の表示手段５５に表示される。ここで、表示手段５５には、ＣＲＴ（Cathode Ray Tube）や液晶表示装置等を使用することができる。また、予測結果は、必要に応じて設けられたプリンタ（図示せず）に出力することもできる。また、記憶部５４は、処理部５２に内蔵されるものであっても、他の装置（データベースサーバ）内にあってもよい。このように、上記構造体の物理量予測装置５０は、入出力装置５１を備えた端末装置（図示せず）から通信により処理部５２や記憶部５４にアクセスするものであってもよい。
【００５９】
次に、本発明に係るタイヤの摩耗に関する物理量予測方法の具体的な手順について説明する。図９は、この発明の実施の形態１に係るタイヤの摩耗に関する物理量予測方法の手順例を示すフローチャートである。本発明に係るタイヤの摩耗に関する物理量予測方法においては、まず、摩擦エネルギーの予測対象であるタイヤのモデルを作成する（ステップＳ１０１）。具体的には、補強コード１とゴム２とで構成される複合材層、補強コードの弾性率変化あるいは補強コード１の粘弾性特性等を考慮して、ＦＥＭやＢＥＭ等の解析手法で取り扱い易い形に実際のタイヤを単純化する。タイヤモデルを作成する際には、補強コード１やゴム２等の物性値を同時に入力しておく。
【００６０】
次に、タイヤ１０の境界条件を入力する（ステップＳ１０２）。タイヤ１０の境界条件には、タイヤの空気圧、荷重、スリップ角、速度、あるいはホイールのリム幅といった、タイヤ１０の使用条件を入力する。ステップＳ１０１とステップＳ１０２とによってタイヤ１０の性能予測に必要な条件がすべて入力されるので、作成したタイヤモデル及び入力したデータに基づいて、ＦＥＭやＢＥＭ等の解析手法により、タイヤモデルの走行状態における挙動を解析する（ステップＳ１０３）。
【００６１】
ここで、タイヤ１０全体は、ＦＥＭを用いた場合には、例えば図２で示したようなメッシュによって有限個数の微小要素１０ａ、１０ｂ等に分割されており、計算前には記憶部５４に３次元の座標データとして読み込まれている。この微小要素の座標データと上記入力データとから、各時間における各微小要素のひずみや応力等を算出し、タイヤ１０全体の挙動を求める。ＦＥＭにおいては、各微小要素内における応力やひずみは一定であると仮定し、この仮定に基づいて各微小要素の応力状態を演算する。
【００６２】
次に、タイヤモデルの走行状態における挙動解析結果から、摩擦エネルギーの予測に必要な情報を取得する（ステップＳ１０４）。この情報は、例えばタイヤ１０を構成する節点の座標、変位、接触力、路面の変位、あるいは接地表面に関する節点のデータ等の情報である。これらの情報に基づいて、単位面積当たりの摩擦エネルギーを予測する（ステップＳ１０５）。次に、単位面積当たりの摩擦エネルギーの予測手順について説明する。
【００６３】
図１０は、単位面積当たりの摩擦エネルギーを予測する手順を示すフローチャートである。この予測手順においては、タイヤの初期形状に基づいて、接地する分割要素表面と接地しない分割要素表面とを、各節点が分担する面積を算出する前に判定する。まず、接地する分割要素表面か否かを判定する（ステップＳ２０１）。ここで、この判定手順を説明する。図１１は、接地する分割要素表面か否かを判定する手順を示すフローチャートである。
【００６４】
この接地判定においては、まず、節点座標、及び分割要素表面情報（例えば当該分割要素表面を構成する節点や分割形状等）を取得する（ステップＳ３０１）。次に、当該分割要素表面がタイヤ表面の分割要素表面であるか否かを判定する（ステップＳ３０２）。当該分割要素表面がタイヤ表面である場合には、当該分割要素表面が接地するか否かを判定し、その結果を、例えば物理量予測装置５０の記憶部５４（図８参照）に記憶させる（ステップＳ３０３）。なお、本手順をコンピュータで実行させる場合には、ＲＡＭ上やハードディスクの所定記憶領域に前記判定結果を記憶させる。分割要素表面が接地するか否かの判定には、例えば上述した基準ベクトルと分割要素表面の法線ベクトルとのなす角に基づく手法を用いることができる。タイヤ表面に存在する全分割要素表面について処理が終了したか否かを判断し（ステップＳ３０４）、全分割要素表面について処理が終了したら接地の判定手順が終了する。
【００６５】
次に、図１０に戻って説明する。分割要素表面の接地判定手順が終了したら、分割要素の節点が接地しているかどうかを判断する（ステップＳ２０２）。この判断手法は、例えば当該節点と路面との垂直接触力が０を超える場合には、当該節点は接地すると判断する。判定対象である節点が接地していると判定された場合には、当該節点が分担する面積（以下節点面積）を算出する（ステップＳ２０３）。節点面積の算出手順については後述する。
【００６６】
次に、路面と平行な接触力の成分からせん断接触力Ｆｓを算出し（ステップＳ２０４）、また、当該時間ステップと直前の時間ステップとの間における当該節点の変位Ｄｎと路面変位Ｄｒとの差Ｄｎ−Ｄｒから、接地面１１ｐと路面とのすべり量Ｌとを算出する（ステップＳ２０５）。なお、ステップＳ２０３〜Ｓ２０５の順序は問わない。せん断接触力Ｆｓ、すべり量Ｌ及び節点面積Ａから、式（１）に基づいて単位面積当たりの摩擦エネルギーＥを算出する（ステップＳ２０６）。この値を、当該節点が接地を開始してからにおける単位面積当たりの総摩擦エネルギーに追加する（ステップＳ２０７）。
Ｅ＝Ｆｓ×Ｌ／Ａ・・・（１）
ここでＡは節点面積である。
【００６７】
全節点の処理が終了したか否かを判断し（ステップＳ２０８）、全節点の処理が終了していない場合にはステップＳ２０２〜Ｓ２０７を繰り返す。当該ステップにおいて全節点を処理したら全ステップの処理が終了したか否かを判断し（ステップＳ２０９）、全ステップの処理が終了していない場合にはステップＳ２０２〜Ｓ２０８を繰り返して全ステップを処理して、摩擦エネルギーＥの予測手順が終了する。
【００６８】
次に、節点面積の算出手順について説明する。図１２は、節点面積の算出手順を示すフローチャートである。まず、節点面積を算出するため、面積算出対象の節点を含む分割要素表面情報（当該分割要素表面を構成する節点の座標や分割要素表面の分割形状等）を取得する（ステップＳ４０１）。なお、摩擦エネルギーＥを求める際にはトレッド内部の節点情報は不要なので、表面に存在する節点のみが節点面積算出の対象となる。次に、上述した分割要素表面の接地判定手順で記憶手段に記憶させた、接地する分割要素表面か否かの判定結果を取得する（ステップＳ４０２）。なお、ステップＳ４０１とＳ４０２との順序は問わない。
【００６９】
節点面積を算出する手順においては接地する分割要素表面が対象となるので、当該節点を含む分割要素表面が接地しているか否かを判断する（ステップＳ４０３）。これは、ステップＳ４０２で取得した、接地する分割要素表面か否かの判定結果に基づいて判断することができる。
【００７０】
接地する分割要素表面である場合には、当該分割要素表面の面積と分割形状とから面積算出対象である節点の節点面積を算出して、既に求めた節点面積に加算する（ステップＳ４０４）。そして、接地している面積算出対象の節点を含む各分割要素表面すべてについて面積算出対象の節点面積を算出し、既に求めた節点面積に加算したら、当該節点の面積算出が終了する。そして、全分割要素表面の処理が終了していない場合にはステップＳ４０３及びＳ４０４を繰り返して全分割要素表面を処理して（ステップＳ４０５）、全節点面積の算出が終了する。
【００７１】
実施の形態１に係るタイヤの摩耗に関する物理量予測方法では、タイヤの初期形状に基づいて、節点面積の算出前に分割要素表面の接地を判定し、単位面積当たりの摩擦エネルギーを算出する。これにより、各時間ステップ毎に各分割要素表面の接地状態を判定する必要はないので計算量を少なくできる。その結果、コンピュータシミュレーションにおいては計算に要する時間を短くでき、また、ハードウェア資源が限られている場合でも有効な評価ができる。
【００７２】
（変形例）
次に、実施の形態１の変形例に係るタイヤの摩耗に関する物理量予測方法の手順について説明する。この予測手順は、上記予測手順と略同様であるが、単位面積当たりの摩擦エネルギーを予測する手順において、タイヤ１０の転動時における変形形状を基に接地する分割要素表面と接地しない分割要素表面とを各時間ステップ毎に判定する点が異なる。他の構成は、上記予測手順と同様なのでその説明を省略する。
【００７３】
この変形例に係るタイヤの摩耗に関する物理量予測方法において、タイヤモデルの作成から単位面積当たりの摩擦エネルギーを予測するまでにおける一連の流れは、実施の形態１で説明した通りなのでその説明を省略する（図９参照）。図１３は、実施の形態１の変形例に係る単位面積当たりの摩擦エネルギーを予測する手順を示すフローチャートである。なお、本予測手順のステップＳ５０１〜ステップＳ５０８は、実施の形態１に係る単位面積当たりの摩擦エネルギーを予測する手順のステップＳ２０２〜Ｓ２０９（図１０参照）と同様なので、説明を省略する。
【００７４】
本変形例に係る節点面積の算出手順（ステップＳ５０２）は次の通りである。図１４は、節点面積の算出手順を示すフローチャートである。まず、節点面積を算出するため、面積算出対象の節点を含む分割要素表面情報（当該分割要素表面を構成する節点の座標や分割要素表面の分割形状等）を取得する（ステップＳ６０１）。次に、各時間ステップ毎に当該節点を含む分割要素表面が接地しているか否かを判断する（ステップＳ６０２）。節点面積を算出する手順においては接地する分割要素表面が対象となるからである。この判断手法には上述した手法を用いることができ、例えば、分割要素表面の法線ベクトル及び路面の法線ベクトルのなす角と所定の接地基準角とに基づく手法を用いることができる。
【００７５】
接地する分割要素表面である場合には、当該分割要素表面の面積と分割形状とから面積算出対象である節点の節点面積を算出して、既に求めた節点面積に加算する（ステップＳ６０３）。そして、接地している面積算出対象の節点を含む各分割要素表面すべてについて面積算出対象の節点面積を算出し、既に求めた節点面積に加算したら、当該節点の面積算出が終了する。そして、全分割要素表面の処理が終了していない場合にはステップＳ６０２及びＳ６０３を繰り返して全分割要素表面を処理して（ステップＳ６０４）、全節点面積の算出が終了する。このようにして求めた各時間ステップにおける節点面積に基づいて、単位面積当たりの摩擦エネルギーを算出する。
【００７６】
本変形例に係るタイヤの摩耗に関する物理量予測方法では、各時間ステップ毎に接地する分割要素表面を判定し、単位面積当たりの摩擦エネルギーを算出する。これにより、タイヤの転動時における変形形状に応じて接地面を判断できるので、より高い精度でタイヤ接地面の摩耗を予測することができる。
【００７７】
実施の形態１に係る本発明のタイヤの摩耗に関する物理量の予測方法は、上記構造体の物理量予測装置５０を使用しなくとも、予め用意されたプログラムをパーソナル・コンピュータやワークステーションなどのコンピュータで実行することによって実現することができる。このプログラムは、インターネットなどのネットワークを介して配布することができる。また、このプログラムは、ハードディスク、フレキシブルディスク（ＦＤ）、ＣＤ−ＲＯＭ、ＭＯ、ＤＶＤなどのコンピュータで読み取り可能な記録媒体に記録され、コンピュータによって記録媒体から読み出されることによって実行することもできる。以下の実施の形態においても同様である。
【００７８】
（実施例）
ここでは、実施の形態１で説明した摩耗エネルギーの予測方法を用いて、タイヤの摩耗に関する物理量のうち垂直接触力と摩擦エネルギーとを評価した。タイヤのサイズが２１５／４５Ｒ１７の乗用車用ラジアルタイヤをモデル化し、ＦＥＭによって垂直接触力と摩擦エネルギーとを予測した。図１５は、比較例１の垂直接触力分布を示す説明図、図１６は、実施例１の垂直接触力分布を示す説明図、図１７は、実施例２の垂直接触力分布を示す説明図である。また、図１８は、比較例２の摩擦エネルギー分布を示す説明図、図１９は、実施例３の摩擦エネルギー分布を示す説明図、図２０は、実施例４の摩擦エネルギー分布を示す説明図である。
【００７９】
実施例１、２及び比較例１では、内圧２３０ｋＰａ、荷重３．８ｋＮで負荷した際の垂直接触力を評価した。また、実施例３、４及び比較例２では、内圧２３０ｋＰａ、荷重３．８ｋＮ、速度３６ｋｍ／ｈ、キャンバー角１度、旋回力０．２Ｇにおける摩擦エネルギーを評価した。実施例１、３は単位面積当たりの垂直接触力又は摩擦エネルギーによるものであるが、溝壁１１ｓも含めて計算したものである。実施例２、４は単位面積当たりの垂直接触力又は摩擦エネルギーによるものであり、さらに節点面積の算出においては溝壁１１ｓを除外したものである。比較例１、２は単位面積当たりの垂直接触力又は摩擦エネルギーに換算していないものである。
【００８０】
比較例１及び２においては、単位面積当たりの垂直接触力又は摩擦エネルギーを用いていないので、垂直接触力又は摩擦エネルギーの分布に所々ムラが発生していることがわかる（図１５、１８中Ｂで示す領域）。また、図１５〜１７からわかるように、実施例１及び３においては、単位面積当たりの垂直接触力又は摩擦エネルギーとして算出しているので、比較例１、及び２と同じ場所においては、これらよりも垂直接触力又は摩擦エネルギーが一様に分布していることがわかる。ただし、溝壁１１ｓを含めてこれらを算出しているので、溝１８に接するブロック端部の値が小さくなっている。
【００８１】
これに対して、実施例２及び４においては、単位面積当たりの垂直接触力又は摩擦エネルギーに換算し、節点面積の算出においては溝壁１１ｓを除外しているので、上記例と同じ場所において垂直接触力又は摩擦エネルギーの分布にムラは発生せず、また溝１８に接するブロック端部においても適切な値を示す。このように、本発明に係る摩擦エネルギーの予測方法によれば、地面と接触するタイヤ１０の表面のみならず、溝１８に接するブロック端部においても適切に摩擦エネルギーを評価できるようになった。これにより、タイヤ１０の摩耗や偏摩耗を高い精度で予測できるので、タイヤ１０の実設計に対して有効に使用できる。
【００８２】
（実施の形態２）
実施の形態２においては、実施の形態１で説明した摩擦エネルギーの予測方法を利用して、タイヤの偏摩耗及びタイヤの摩耗を予測する方法について説明する。図２１は、この発明の実施の形態２に係るタイヤの偏摩耗を予測する手順を示すフローチャートである。まず、自由転動、制動、駆動、旋回等の実運転条件下において、それぞれ摩擦エネルギーＥ１、Ｅ２・・・Ｅｎを予測する（ステップＳ７０１）。摩擦エネルギーの予測手順は実施の形態１で説明した通りである。
【００８３】
次に、これらの運転条件がどの程度の頻度、あるいはどの程度の大きさで発生するかに基づいて、それぞれの条件の重み付けをする（ステップＳ７０２）。この重み付けは、各運転条件の発生頻度割合に応じて所定の大きさに各重み付け係数Ｃ１、Ｃ２・・・Ｃｎを決定する。そして、各重み付け係数を前記各運転条件下における摩擦エネルギーに乗じ、これらの和Ｃ１×Ｅ１＋Ｃ２×Ｅ２＋・・・＋Ｃｎ×Ｅｎを求める。これが実走行中における摩擦エネルギーＥａとなる（ステップＳ７０３）。この摩擦エネルギーＥａに基づいて、実車走行中におけるタイヤの偏摩耗を予測する（ステップＳ７０４）。
【００８４】
図２２は、この発明の実施の形態２に係るタイヤの偏摩耗及び摩耗を予測する手順を示すフローチャートである。まず、ゴム摩擦指数Ｇを、例えばＪＩＳＫ６２４６に規定する通常のランボーン摩耗試験により求める（ステップＳ７０１'）。次に、自由転動、制動、駆動等の実運転条件下において、それぞれ摩擦エネルギーＥ１、Ｅ２・・・Ｅｎを予測し（ステップＳ７０２'）、各条件の発生頻度に応じて各条件の重み付けをして（ステップＳ７０３'）、実車走行中における摩擦エネルギーＥａを予測する（ステップＳ７０４'）。これらの手順は上述した通りなので説明を省略する。このようにして求めた摩擦エネルギーＥａに基づいて、実車走行中におけるタイヤの偏摩耗及び摩耗を予測する（ステップＳ７０５'）。例えば、ゴム摩擦指数ＧとＥａとの積Ｇ×Ｅａを求め、基準となるタイヤのＧ×Ｅａ_bと比較することにより摩耗速度が予測できる。
【００８５】
このように、ゴム摩擦指数Ｇを使用すれば、摩擦エネルギーの分布のみならず摩擦速度も予測できるので、幅広いタイヤの摩耗評価に対して適用できる。なお、実施の形態２で説明したタイヤの摩耗に関する物理量の予測方法は、実施の形態１と同様にコンピュータによって実行できる。
【００８６】
【発明の効果】
以上説明したように、この発明に係るタイヤの摩耗に関する物理量の予測方法、又はこの発明に係るタイヤの摩耗に関する物理量の予測プログラムでは、タイヤモデルを構成する複数の分割要素表面に含まれる節点が分担する分割要素表面の面積によって、タイヤの摩耗に関する物理量を単位面積当たりの物理量に換算してタイヤの摩耗を評価するようにした。このため、せん断接触力や摩擦エネルギー等の物理量は要素分割や節点位置の影響を受けないので、接触面全体における摩耗を精度よく予測することができる。また、摩擦エネルギーを測定する必要がないので、効率よくタイヤの摩耗を予測して、開発の効率を向上させることができる。
【００８７】
また、この発明に係るタイヤの摩耗に関する物理量の予測方法、又はこの発明に係るタイヤの摩耗に関する物理量の予測プログラムでは、分割要素表面の分割形状に基づいて各節点が分担する面積を求めるようにした。これによって、分割要素表面の形状に関わらず節点が分担する面積を正確に算出して、タイヤの摩耗に関する物理量を精度よく予測できる。
【００８８】
また、この発明に係るタイヤの摩耗に関する物理量の予測方法、又はこの発明に係るタイヤの摩耗に関する物理量の予測プログラムでは、分割要素表面が接地しているか否かを判定する工程を含むようにした。これによって、例えば溝壁等における分割要素表面のように、接地しない分割要素表面を物理量の予測に際して排除することができる。その結果、溝に接するブロック端部における摩擦エネルギーを適切に評価できるので、このような部分における摩耗も精度よく予想できる。従来摩擦エネルギーを測定することが極めて困難であったブロックの端部等も、分割要素表面の接地状態を判定することによって、精度よく摩擦エネルギーを予測することができる。
【００８９】
また、この発明に係るタイヤの摩耗に関する物理量の予測方法では、タイヤの初期形状を基準として、各節点が分担する面積を算出する前に予め分割要素表面が接地しているか否かを判定するようにした。このため、分割要素表面が接地するか否かは、タイヤの初期形状を基にして１回だけ判定するので、計算数を抑えることができる。特に、コンピュータを用いて本発明に係る予測方法を実行する場合には、ハードウェア資源を節約できるという長所がある。
【００９０】
また、この発明に係るタイヤの摩耗に関する物理量の予測方法では、タイヤの変形形状に応じて、分割要素表面が接地しているか否かを各瞬間毎に判定するようにしたので、タイヤの摩耗に関する物理量の予測精度をさらに向上させることができる。
【００９１】
また、この発明に係るタイヤの摩耗に関する物理量の予測方法、又はこの発明に係るタイヤの摩耗に関する物理量の予測プログラムでは、単位面積当たりの摩擦エネルギーを用いてタイヤの摩耗を予測するようにしたので、精度の高い摩耗予測ができる。これによって、タイヤの設計に対する有効な情報を得ることができる。また、摩擦エネルギーを測定する必要がないので、タイヤの摩耗を予測する手間も大幅に軽減でき、タイヤの開発期間を短縮できる。
【００９２】
また、この発明に係るタイヤの摩耗に関する物理量の予測方法では、ゴム摩擦指数を使用するようにしたので、単なる摩耗の分布だけでなく実際の摩耗速度も予測でき、より幅広くタイヤの摩耗に関する評価ができる。
【００９３】
また、この発明に係るタイヤの摩耗に関する物理量の予測方法では、上記タイヤの摩耗に関する物理量の予測方法を実行させるための処理手段を備えるようにした。これによって、せん断接触力や摩擦エネルギー等の物理量は要素分割や節点位置の影響を受けないで求めることができるので、接触面全体における摩耗を精度よく予測することができる。また、摩擦エネルギーを測定する必要がないので、効率よくタイヤの摩耗を予測して、開発の効率を向上させることができる。
【図面の簡単な説明】
【図１】予測対象であるタイヤを、その回転軸を含む子午面で切った断面を示す一部断面図である。
【図２】予測対象であるタイヤを微小要素に分割した例を示す斜視図である。
【図３】キャップトレッドが地面と接する接地面を微小要素に分割した例を示す説明図である。
【図４】四角形要素によって９分割した接地面を示す平面図である。
【図５】タイヤモデルの接地面を含む分割要素の一つを示す斜視図である。
【図６】分割要素表面の接地を判定する方法の一例を示す説明図である。
【図７】分割要素表面の接地を判定する方法の他の例を示す説明図である。
【図８】この発明の実施の形態１に係るタイヤの摩耗に関する物理量予測装置の一例を示す説明図である。
【図９】この発明の実施の形態１に係るタイヤの摩耗に関する物理量予測方法の手順例を示すフローチャートである。
【図１０】単位面積当たりの摩擦エネルギーを予測する手順を示すフローチャートである。
【図１１】接地する分割要素表面か否かを判定する手順を示すフローチャートである。
【図１２】節点面積の算出手順を示すフローチャートである。
【図１３】実施の形態１の変形例に係る単位面積当たりの摩擦エネルギーを予測する手順を示すフローチャートである。
【図１４】節点面積の算出手順を示すフローチャートである。
【図１５】比較例１の垂直接触力分布を示す説明図である。
【図１６】実施例１の垂直接触力分布を示す説明図である。
【図１７】実施例２の垂直接触力分布を示す説明図である。
【図１８】比較例２の摩擦エネルギー分布を示す説明図である。
【図１９】実施例３の摩擦エネルギー分布を示す説明図である。
【図２０】実施例４の摩擦エネルギー分布を示す説明図である。
【図２１】この発明の実施の形態２に係るタイヤの偏摩耗を予測する手順を示すフローチャートである。
【図２２】この発明の実施の形態２に係るタイヤの偏摩耗及び摩耗を予測する手順を示すフローチャートである。
【符号の説明】
１０タイヤ
１１キャップトレッド
１１Ｂブロック
１１ｐ接地面
１１ｓ溝壁
１８溝
２０、２１分割要素表面
３０路面
５０物理量予測装置
５２処理部
５３入力手段
５４記憶部[0001]
BACKGROUND OF THE INVENTION
The present invention relates to a method for predicting a physical quantity related to tire wear, and more particularly to a method and apparatus for predicting a physical quantity related to tire wear, and a computer program, which can more accurately predict the shear contact force and wear distribution of the tire. .
[0002]
[Prior art]
As a method for predicting tire wear, for example, methods disclosed in Patent Documents 1 and 2 are known. The prediction methods disclosed in these are composed of the following procedures. That is, (1) a procedure for measuring the rubber wear index of the tire surface, (2) a procedure for measuring the friction energy of the tire under various operating conditions, and (3) each operating condition according to the frequency of each operating condition. Tire wear is predicted through a weighting procedure.
[0003]
[Patent Document 1]
Japanese Patent Laid-Open No. 11-326143 (pages 2 to 23, FIGS. 1 to 12)
[Patent Document 2]
JP 2001-1723 A (pages 2-9, FIGS. 1-10)
[0004]
[Problems to be solved by the invention]
As described above, in the conventional tire wear prediction method, since the wear of the tire is predicted by measuring the friction energy, it takes time to measure the friction energy. In addition, since the frictional energy is measured by points, enormous measurement is required to evaluate all tread surfaces, including uneven wear within a single block on the tire surface. Furthermore, it is very difficult to measure the frictional energy at a portion where the measurement point is extremely small, such as an end portion of a block or rib, or a portion where fine grooves are densely packed. As a result, the prediction accuracy for tire wear could not be sufficiently increased, and the prediction result could not be sufficiently reflected in the design.
[0005]
Therefore, the present invention has been made in view of the above, and by omitting the measurement of frictional energy, at least one of predicting tire wear efficiently and predicting tire wear more accurately is provided. An object of the present invention is to provide a method and apparatus for predicting a physical quantity related to tire wear that can achieve one, and a computer program.
[0006]
[Means for Solving the Problems]
In order to omit the measurement of the frictional energy, the present inventors used a method for predicting the frictional energy by a finite element method, a boundary element method, or other analysis methods. When the frictional energy was obtained by a simulation using such an analysis method, it was found that the evaluation results of the partial wear and the wear rate of the tire differ between the simulation and the actual measurement value. The present inventors found this cause as a result of earnest research.
[0007]
In an analysis method such as a finite element method, a prediction region such as a tire block is divided into a plurality of elements and the elements are regarded as uniform. Then, the shear contact force and the slip amount are obtained for each node constituting the element, and the friction energy is calculated from both. However, when the prediction area having a certain size is divided by multiple elements and the multiple nodes are in contact with each other, the vertical contact force, shear contact force, and friction energy shared by each node are the number of element divisions and the number of nodes. It depends on the position. Thus, since the magnitude of the frictional energy at each node is affected by the number of element divisions and the node position, the wear evaluation results in the simulation and the actual measurement value are different. Thus, in order to eliminate the influence of the number of element divisions and node positions, the above problem was solved by evaluating wear based on vertical contact force and frictional energy per unit area.
[0008]
The present invention The method for predicting physical quantities related to tire wear is to create a tire model by dividing the tire into a plurality of divided elements in order to predict tire wear by the finite element method, boundary element method, or other analysis methods. procedure And determining whether or not the node existing on the surface of the dividing element is grounded procedure When, If the node is grounded, exclude the area of the surface of the dividing element that is not grounded, Calculate the area of the surface of the dividing element shared by the nodes procedure And the physical quantity related to tire wear is converted into a physical quantity per unit area according to the area shared by the calculated nodes. procedure When If the nodes are not grounded, a procedure for determining whether or not the processing for all the nodes has been completed is executed. And features.
[0009]
Also, In the present invention The physical quantity prediction program related to tire wear is a procedure for creating a tire model by dividing a tire into a plurality of elements when predicting a physical quantity related to tire wear by a finite element method, a boundary element method, or other analysis methods, A procedure for determining whether or not a node constituting the surface of the dividing element is grounded; If the node is grounded, exclude the area of the surface of the dividing element that is not grounded, A procedure for calculating the area of the surface of the dividing element shared by the nodes, and a procedure for converting a physical quantity related to tire wear into a physical quantity per unit area according to the calculated area shared by the nodes; If the nodes are not grounded, a procedure for determining whether or not the processing for all the nodes has been completed is executed. And making the computer execute.
[0010]
This physical quantity prediction method or prediction program for tire wear is based on the area of the surface of each of the divided elements shared by the nodes included in the surfaces of the plurality of divided elements constituting the tire model. Tire wear is evaluated in terms of physical quantities. For this reason, physical quantities such as shear contact force and frictional energy are not affected by element division or node position, so that wear on the entire contact surface can be accurately predicted.
[0011]
Also, a physical quantity prediction program relating to tire wear that causes a computer to execute each procedure of the physical quantity prediction method relating to tire wear according to the present invention. To According to this, tire wear can be predicted on a computer. Thereby, since it is not necessary to measure friction energy, tire wear can be predicted efficiently and development efficiency can be improved. Further, the friction energy can be predicted easily and accurately even in a portion where the measurement point is extremely small, such as a portion where the end of the block or the narrow groove is dense, and it is extremely difficult to measure the friction energy. Thereby, since the prediction accuracy with respect to the wear of the tire can be sufficiently increased, the prediction result can be sufficiently utilized for the design.
[0012]
Here, the physical quantity related to tire wear is the normal contact force acting on the tire contact surface, the maximum friction force calculated by multiplying the vertical contact force by the friction coefficient, shear contact force, friction energy, wear amount, slip amount, Amount of wear and other physical quantities related to wear (hereinafter the same). Of these, in the present invention, the vertical contact force acting on the tire contact surface, the maximum contact force calculated by multiplying the vertical contact force by the friction coefficient, shear contact force, friction energy, and other physical quantities per unit area can be converted. Items are handled in terms of physical quantities per unit area (the same applies hereinafter).
[0013]
Further, when evaluating the wear of the tire, it is necessary to use the friction energy per unit area at a certain moment added for the entire time from when the tire surface starts to ground until it stops. In this case, in addition to the step of the invention, the step of sequentially adding the obtained friction energy per unit area to the total friction energy per unit area after the tire surface starts to contact Friction energy per unit area is calculated for the entire time when is in contact with the ground.
[0014]
Also, The present invention The method for predicting a physical quantity related to tire wear according to the present invention calculates the area of the surface of the dividing element shared by the nodes. procedure In, calculate the area of the surface of the dividing element including the node procedure And determining the ratio of the nodes sharing the surface area of the dividing element based on the division shape of the surface of the dividing element including the node. procedure And the area where the nodes share the surface area of the dividing element is added to the surfaces of all the dividing elements including the nodes. procedure It is characterized by including these.
[0015]
Also, The present invention In the physical quantity prediction program related to tire wear according to the above-described physical quantity prediction program related to tire wear, in the procedure for calculating the surface area of the dividing element shared by the nodes, the dividing element including the nodes A procedure for calculating the area of the surface, a procedure for determining a ratio of the nodes to share the area of the surface of the dividing element, based on a division shape of the surface of the dividing element including the nodes, and And causing the computer to execute a procedure of adding the area sharing the surface area of the element to the surfaces of all the divided elements including the nodes.
[0016]
In the present invention, the area shared by the nodes is obtained based on the division shape of the surface of the division element. Then, the nodes share the physical quantity that can be converted into the normal contact force acting on the tire contact surface, the maximum contact force calculated by multiplying the contact force by the friction coefficient, the shear contact force, the friction energy, and other physical quantities per unit area. The area to be converted into a physical quantity per unit area. Thus, the area shared by the nodes can be accurately calculated regardless of the shape of the surface of the dividing element, so that the physical quantity related to tire wear can be correctly converted into the physical quantity per unit area. In addition, a physical quantity prediction program related to tire wear for causing a computer to execute the procedure of the present invention. To Even in this case, the same actions and effects can be achieved.
[0017]
Also, The present invention According to the method for predicting a physical quantity related to tire wear, the method for predicting a physical quantity related to tire wear described above, Based on the deformation shape at the time of rolling of the tire model, it is determined at each time step whether or not the node existing on the surface of the dividing element is grounded. Features.
[0018]
Also, The present invention The physical quantity prediction program related to tire wear according to the above-described physical quantity prediction program related to tire wear is Based on the deformation shape at the time of rolling of the tire model, it is determined at each time step whether or not the node existing on the surface of the dividing element is grounded. Is executed by a computer.
[0019]
Since the present invention includes a step of determining whether or not the surface of the dividing element is grounded in the above invention, the same operations and effects as the above invention are achieved. Furthermore, since the step of determining whether or not the surface of the dividing element is grounded is included, for example, the node included in the surface of the dividing element is grounded, such as the surface of the dividing element in a groove wall or the like. In the prediction of the physical quantity, the surface of the divided element that is not grounded by the element surface itself can be excluded. As a result, the frictional energy at the block end contacting the groove can be appropriately evaluated, so that wear at such a portion can be accurately predicted.
[0020]
In addition, since the measurement points at the end of the block are extremely small, it was extremely difficult to measure the friction energy in the past. Energy can be predicted. Thereby, since the prediction accuracy with respect to the wear of the tire can be sufficiently increased, the prediction result can be sufficiently utilized for the design. Also, a physical quantity prediction program related to tire wear for causing a computer to execute the procedure of the present invention. In Even if it is, it has the same operation and effect.
[0021]
Also, The present invention According to the method for predicting a physical quantity related to tire wear, the method for predicting a physical quantity related to tire wear described above, Procedure for judging whether or not the node existing on the surface of the above-mentioned dividing element is grounded In ,Up It is characterized by determining whether the surface of the said division | segmentation element is earth | grounding on the basis of the initial shape of the said tire.
[0022]
In the present invention, it is determined in advance whether or not the surface of the dividing element is in contact with the ground before calculating the area shared by the nodes based on the initial shape of the tire. In this way, whether or not to ground is determined only once based on the initial shape of the tire, so the number of calculations can be reduced. In particular, when the prediction method according to the present invention is executed using a computer, there is an advantage that hardware resources can be saved. In addition, as a specific determination method, for example, when an angle formed by a normal vector of the surface of the dividing element and a vector directed from the center of the tire toward the center of the surface of the dividing element is smaller than a predetermined angle In some cases, it is determined that the surface of the dividing element is grounded.
[0023]
Also, The present invention According to the method for predicting a physical quantity related to tire wear according to the present invention, in the method for predicting a physical quantity related to tire wear, it is further determined whether or not the surface of the dividing element is grounded according to a deformed shape of the tire. It is characterized by that. In the present invention, it is determined at each moment whether or not the surface of the dividing element is in contact with the ground according to the deformed shape of the tire, so that it is possible to further improve the physical quantity prediction accuracy related to tire wear. When the present invention is executed by a computer, whether or not the surface of the dividing element is grounded is determined at each time step.
[0024]
For example, when the angle formed between the normal vector of the surface of the division element to be determined and the normal vector of the ground is equal to or greater than a predetermined angle, the determination method of grounding is determined to contact the surface of the division element. There is something to do. As another determination method, the vertical coordinate of the surface of the division element to be determined is compared with the vertical coordinate of the ground, and if the two match, the surface of the division element is grounded. There is something to judge. Furthermore, as described above, when the angle formed by the normal vector of the surface of the dividing element to be determined and the vector from the center of the tire toward the center of the surface of the dividing element is smaller than a predetermined angle The surface of the dividing element may be determined to be grounded.
[0025]
Also, The present invention The physical quantity prediction method related to the wear of the tire is based on free rolling, braking, driving, turning and other driving conditions. Above A method of predicting frictional energy per unit area under each of the above operating conditions by a method of predicting physical quantities related to tire wear, and weighting each operating condition in consideration of the frequency of occurrence of each operating condition And a step of predicting friction energy per unit area during actual vehicle travel from the friction energy per unit area under each driving condition and the weighting of each driving condition. .
[0026]
Also, The present invention The physical quantity prediction program related to tire wear is based on free rolling, braking, driving, turning and other driving conditions. Above In consideration of the procedure for predicting the friction energy per unit area under each of the above operating conditions and the frequency of occurrence of each of the operating conditions by the method for predicting the physical quantity related to tire wear, the weight of each of the operating conditions is set. And a procedure for predicting the friction energy per unit area during actual vehicle travel from the friction energy per unit area under each driving condition and the weighting of each driving condition. To do.
[0027]
In the present invention, since wear of a tire is predicted using friction energy per unit area, it is possible to predict wear with high accuracy. This makes it possible to obtain effective information for the tire design. Further, since it is not necessary to measure frictional energy, the labor for predicting tire wear can be greatly reduced, and the tire development period can be shortened. Note that a physical quantity prediction program related to tire wear for causing a computer to execute the procedure of the present invention. To Even in this case, the same actions and effects can be achieved.
[0028]
Also, The present invention In the method for predicting a physical quantity related to tire wear, such as the method for predicting a physical quantity related to tire wear according to the present invention, the method further includes the step of measuring the rubber wear index of rubber constituting the surface of the tire, Tire wear may be predicted by multiplying the frictional energy per unit area during actual vehicle travel. As described above, if the rubber friction index is used, not only the distribution of wear but also the actual wear rate can be predicted, so that the evaluation on the wear of the tire can be performed more widely.
[0029]
Also, The present invention The apparatus for predicting a physical quantity related to tire wear according to the present invention includes a processing means for processing each procedure in the method for predicting a physical quantity related to tire wear, and physical property values, operating conditions, boundary conditions of the material constituting the tire in the processing means. Input means for providing other data and display means for displaying a prediction result by the processing means.
[0030]
As described above, the physical quantity prediction device related to tire wear according to the present invention includes processing means for executing the physical quantity prediction method related to tire wear, so that physical quantities such as shear contact force and friction energy are present. Wear on the entire contact surface can be accurately predicted without being affected by element division and node position. Further, since it is not necessary to measure the friction energy, it is possible to efficiently predict tire wear and improve the development efficiency. Furthermore, the friction energy can be predicted easily and accurately even in a portion where the measurement point is extremely small, such as a portion where the end of the block or the narrow groove is dense, and it is extremely difficult to measure the friction energy. Thereby, since the prediction accuracy with respect to tire wear can be sufficiently increased, the prediction result can be sufficiently reflected in the design.
[0031]
DETAILED DESCRIPTION OF THE INVENTION
Hereinafter, the present invention will be described in detail with reference to the drawings. Note that the present invention is not limited to the embodiments. In addition, constituent elements in the following embodiments include those that can be easily assumed by those skilled in the art or that are substantially the same.
[0032]
(Embodiment 1)
FIG. 1 is a partial cross-sectional view showing a cross section of a tire to be predicted cut along a meridian plane including a rotation axis thereof. The structure of the tire 10 that is a prediction target in the present embodiment will be briefly described with reference to FIG. The cap tread 11 is a rubber layer that is disposed on the road surface grounding portion of the tire 10 and covers the outside of the carcass 15, the belt 14, or the breaker. The cap tread 11 serves to protect the carcass 15 and the belt 14 against cut impact. In addition, a groove 18 is formed in the cap tread 11 in order to enhance drainage during rainy weather travel. A portion partitioned by the groove 18 becomes a block of the tire 10.
[0033]
The undertread 12 is a rubber layer disposed between the cap tread 11 and the belt 14 and is used for the purpose of improving heat generation, adhesion, and the like. The side tread 13 is disposed on the outermost side of the sidewall portion to prevent external scratches from reaching the carcass 15 and, in the case of a radial tire, has an auxiliary role of transmitting driving force from the axle to the road surface. I'm in charge.
[0034]
The belt 14 is a rubberized cord layer disposed between the cap tread 11 and the carcass 15. In the case of a bias tire, it is called a breaker. In the radial tire, the belt 14 plays an important role as a shape maintaining and strength member of the tire 10. The carcass 15 is a rubberized cord layer that forms the skeleton of the tire 10. The carcass 15 is a strength member that serves as a pressure vessel when the tire 10 is filled with air or nitrogen. The carcass 15 has a structure that supports a load by its internal pressure and withstands a dynamic load during traveling.
[0035]
The bead 16 is a ring in which a bundle of steel wires supporting the cord tension of the carcass 15 generated by internal pressure is hardened with hard rubber. The bead 16 serves as a strength member of the tire 10 together with the carcass 15, the belt 14, and the tread in addition to the role of fixing the tire 10 to the rim of the wheel. The bead filler 17 is a rubber that fills a space generated when the carcass 15 is wound around the bead wire. The bead filler 17 fixes the carcass 15 to the bead 16 and adjusts the shape of the portion to increase the rigidity of the entire bead portion.
[0036]
A rotational force from a wheel (not shown), that is, a driving force, a braking force, or a turning force for moving the vehicle body is transmitted to the tire 10 via the bead 16. When such driving force, braking force, or turning force for moving the vehicle body is transmitted to the road surface, a shear contact force is generated between the cap tread 11 and the road surface. At the same time, the shear contact force causes the ground surface 11p of the tire 10 and the road surface to slip, and the ground surface 11p of the tire 10 is worn. Next, the physical quantity prediction method related to tire wear according to the present invention will be described. Here, the friction energy is predicted as a physical quantity related to tire wear. However, the prediction procedure described below can be similarly applied to a case where other physical quantities such as shear contact force and vertical contact force are predicted.
[0037]
In the present invention, FEM (Finite Element Method) is used as an analysis method for predicting physical quantities related to tire wear. Note that the analysis method applicable to the physical quantity prediction method related to tire wear according to the present invention is not limited to FEM, and BEM (Boundary Element Method), FDM (Finite Differences Method), and the like can also be used. . It is preferable to select the most appropriate analysis method according to the type of tire, boundary conditions, or the like, or to use a combination of a plurality of analysis methods.
[0038]
FIG. 2 is a perspective view showing an example in which a prediction target tire is divided into minute elements. As shown in the figure, the prediction target tire 10 is divided into a finite number of microelements 10a, 10b, etc. based on the finite element method. Microelements based on the finite element method are, for example, a triangular element in a two-dimensional plane, a quadrilateral element, a three-dimensional body such as a tetrahedral solid element, a pentahedral solid element, a hexahedral solid element, a triangular shell element, It is desirable to use an element usable in a computer, such as a shell element such as a rectangular shell element. The microelements divided in this way are identified one by one using the three-dimensional coordinates in the process of analysis.
[0039]
In the physical quantity prediction method relating to tire wear according to the present invention, friction energy per unit area is used in order to correctly evaluate the distribution of friction energy on the contact surface 11p of the tire 10. Next, the reason will be described. FIG. 3 is an explanatory diagram showing an example in which the ground surface on which the cap tread is in contact with the ground is divided into minute elements. As shown in FIG. 3 (a), when the ground contact surface 11p is divided into one, the vertical contact force F acts equally on the nodes N1 to N4, so that F / A vertical contact force of 4 acts.
[0040]
In general, in an analysis technique such as FEM or BEM, the calculation time increases but the calculation accuracy increases as the number of divisions of elements regarded as constant increases. Therefore, the number of element divisions is increased as much as possible in consideration of a balance with calculation time. In addition, since the friction energy becomes high at the end of the block 11B, the evaluation accuracy inside the block 11B decreases if the number of element divisions is small. For this reason, when predicting the frictional energy of the ground plane 11p, as shown in FIG. 3B, the number of divisions of the ground plane 11p is increased. In the example shown in FIG. 3B, since the ground plane 11p is divided into four, the ground plane 11p includes a total of nine nodes from N1 to N9.
[0041]
Here, a vertical contact force of F / 16 acts on the nodes N1, N3, N7, and N9 at the corners of the ground contact surface 11p due to the relationship between the number of nodes and the location of the nodes. Similarly, a vertical contact force of F / 8 acts on each of the nodes N2, N4, N6, and N8 on the side of the ground surface 11p, and a vertical contact force of F / 4 is applied to the node N1 in the center of the ground surface 11p. Works. Here, when the friction coefficient is μ and the slip amount is L, the shear contact force can be obtained by μ × vertical contact force, and the friction energy can be obtained by μ × vertical contact force × L. Therefore, the friction energy at the nodes N1, N3, N7, and N9 is L × μ × F / 16, the friction energy at the nodes N2, N4, N6, and N8 is L × μ × F / 8, and the node N1 The frictional energy at is L × μ × F / 4.
[0042]
Originally, friction energy is uniformly distributed on the ground contact surface 11p in this example. However, if the friction energy is calculated from the shear contact force and the slip amount of each of the nodes N1 to N9, the value of the friction energy varies depending on the number of element divisions, the number of nodes, or the position of the nodes, and the wear of the ground contact surface 11p is correctly evaluated. I can't. Therefore, in the present invention, the wear of the ground contact surface 11p is evaluated based on the friction energy per unit area. At the same time, the vertical contact force and the shear contact force were evaluated as a vertical contact force per unit area (that is, a vertical ground pressure Pn) and a shear contact force per unit area (Ps = μ × Pn). In this way, since the distribution of friction energy on the surface of the ground contact surface 11p can be obtained correctly, the wear of the ground contact surface 11p can be correctly evaluated.
[0043]
When evaluating by a physical quantity such as friction energy per unit area, it is necessary to obtain the area of the surface of the dividing element (hereinafter referred to as the dividing element surface) shared by each node. Next, an example of this procedure will be described. FIG. 4 is a plan view showing a ground plane divided into nine by rectangular elements. Since the ground plane 11p is divided into nine, the nodes are 16 nodes N1 to N16. Here, a case where the area of the node N6 is calculated will be described. First, information on the surface of the dividing element including the node N6 is acquired. In this example, the division element surfaces including N6 are S1, S2, S4, and S5. Further, the information on the surface of the dividing element is information such as the nodes constituting the surface of the dividing element and the coordinates of the nodes.
[0044]
Next, the respective areas of the respective divided element surfaces S1, S2, S4, and S5 are calculated. This can be calculated from the acquired coordinates of the nodes constituting each surface of the divided elements and the shape of each surface of the divided elements. When the area of each division element surface is calculated, the area shared by the node N6 is calculated according to the shape of each division element surface. In the example shown in FIG. 4A, since the ground contact surface 11p is divided by quadrangular elements, the shapes of the respective divided element surfaces S1, S2, S4, and S5 are quadrangular. In this case, that is, when the shape of the surface of the divided element divided by the divided elements is a quadrangle, the node N6 shares ¼ of the area of the divided element surface S1.
[0045]
Thus, when the shape of the surface of the dividing element on the ground contact surface 11p is a quadrangle, the target nodes share a size of 1/4 of the area occupied by the shape. In addition, as shown in FIG. 4B, when the shape of the surface of the dividing element on the ground contact surface is a triangle, the target node shares the size of 1/3 of the area occupied by the shape. Become. Next, all areas of the respective divided element surfaces S1, S2, S4, and S5 shared by the node N6 are added. This area is an area shared by the node N6. The same procedure is applied to the remaining nodes N1, N2, etc., and the area shared by each node on the ground plane 11p is calculated.
[0046]
FIG. 5 is a perspective view showing one of the dividing elements including the contact surface of the tire model. When the area shared by each node is calculated by the above procedure, the area of the surface of the dividing element that is not actually grounded (the hatched portion in FIG. 5A) such as the groove wall 11s shown in FIG. It will be included in the area shared by N1. In this case, the friction energy per unit area is obtained by an area that is larger than the actual ground contact area. Therefore, the wear of the ground contact surface 11p is evaluated by the friction energy per unit area that is smaller than the actual area. . In order to avoid this inconvenience, when calculating the area of the node (for example, N1) including the surface of the split element that is not grounded, the split element surface 20 that is grounded and the split element surface 21 that is not grounded are determined and grounded. It is necessary to determine the area shared by each node, excluding the area of the split element surface 21 that is not (FIG. 5B). Next, this determination method will be described.
[0047]
FIG. 6 is an explanatory diagram illustrating an example of a method for determining the contact of the surface of the split element. FIG. 7 is an explanatory diagram showing another example of a method for determining the contact of the surface of the dividing element. Based on the initial shape of the tire, the method of determining the ground contact of the split element surface is based on the split element surface that is grounded (or may be grounded) and the split element surface that is not grounded (or is unlikely to be grounded). Is calculated before calculating the area shared by the nodes.
[0048]
In this determination method, for example, a vector heading from the center 10c of the tire 10 to the center of the surface of the division element to be determined is used as a reference vector, and the presence / absence of ground contact is determined from the angle formed by the normal vector of the division element surface and the reference vector There is something to judge. In this case, when the angle formed by both vectors is smaller than the ground contact angle, it is determined that the surface of the divided element is grounded. When the angle formed by both vectors is larger than a predetermined value, the surface of the divided element is not grounded. judge. If the angle formed by both vectors is set to a ground contact angle of 35 to 45 degrees, it is possible to determine the presence or absence of ground contact with sufficient accuracy for practical use, so that the wear of the ground contact surface 11p can be accurately evaluated. Here, the ground contact reference angle means an angle at which the surface of the dividing element does not touch when the angle is smaller than the angle.
[0049]
For example, in the example shown in FIG. 6, a vector from the center 10c of the tire 10 toward the center of the division element surface 20 to be determined is a reference vector V11, and the normal vector V12 of the division element surface 20 and the reference vector V11 are used. The angle θ1 is defined as When the ground contact reference angle is set to 35 degrees, for example, the angle θ1 at the split element surface 20 is smaller than the ground reference angle, so the split element surface 20 is determined to be the ground contact surface (11p). On the other hand, the dividing element surface 21 has a reference vector V22 and a normal vector V21. The angle formed by both vectors is θ2, but since θ2 is larger than the ground contact reference angle, it is determined that the split element surface 21 is not grounded. The dividing element surface 21 is a part of the groove wall 11s.
[0050]
According to this determination method, when the area shared by each node is obtained, it is possible to exclude the surface of the divided element that does not come into contact with the ground, so that it is possible to accurately predict the wear of the ground contact surface 11p. Further, since the ground contact state of each split element surface is determined based on the initial shape of the tire, it is not necessary to determine the ground contact state of each split element surface for each time step. As a result, the amount of calculation can be reduced, so that the time required for calculation can be shortened in computer simulation, and effective evaluation can be performed even when hardware resources are limited.
[0051]
As another determination method, there is also a method of determining, on each time step, the surface of the divided element that contacts the ground and the surface of the divided element that does not contact based on the deformed shape at the time of rolling of the tire 10. In this determination method, since the ground contact surface is determined according to the deformation shape at the time of rolling of the tire 10, the wear of the ground contact surface 11p can be predicted with higher accuracy. In this determination method, V11 or V22 may be used as the reference vector, or the normal vectors Vn1 and Vn2 of the division element surfaces 20 and 21 and the normal vector Vbr of the road surface 30 as shown in FIG. And may be used. When the normal vector Vbr of the road surface 30 is used, the contact reference angle is preferably 160 degrees to 179 degrees, and the normal vector of the surface of the dividing element and the normal vector of the road surface are formed with respect to the contact reference angle. When the angle is large, it is determined to be grounded.
[0052]
For example, when the ground contact reference angle is set to 160 degrees, the angle θ1 on the split element surface 20 is larger than the ground reference angle, and therefore it is determined that the split element surface 20 is grounded. On the other hand, since the angle θ2 at the dividing element surface 21 is smaller than the ground contact reference angle, it is determined that the dividing element surface 21 is not grounded. In this determination example, the angle (θ1, θ2 in FIG. 7A) formed by the normal vector of the surface of the dividing element and the normal vector of the road surface is 180 degrees at the maximum.
[0053]
Further, the coordinates R (Xr, Yr, Zr) of the road surface 30 and the coordinates T of the division element surface 20 or the like to be determined. ₁ You may determine the grounding state of the division | segmentation element surface 20 from a difference with (Xt1, Yt1, Zt1) etc. In this case, the contact state of the dividing element surface 20 or the like is determined using the vertical coordinate (the direction indicated by the arrow Z in FIG. 7). For example, when the difference δ = (Zt1−Zr) between Zt1 and Zr is equal to or less than a reference value, it is determined that the dividing element surface 20 is grounded. Such a determination is made for the following reason. That is, in the case of numerical simulation using a computer, numerical errors are unavoidable, and δ cannot be perfectly matched. Therefore, it is necessary to determine whether to ground in consideration of this numerical error. Here, the reference value is preferably about 0.01 mm to 2 mm.
[0054]
Further, the coordinates Nc at the center of the division element surface 20 can be used as the coordinates of the division element surface 20. For example, the average value (Z1 + Z2 + Z3 + Z4) / 4 of the coordinates in the vertical direction of the nodes N1 to N4 included in the dividing element surface 20 can be set as the vertical coordinate value Zc at the center of the dividing element surface 20 (FIG. 7). (See (c)). Then, the coordinate value Zr in the vertical direction of the road surface 30 and the coordinate value Zc in the vertical direction of the dividing element surface 20 are compared, and if both are equal, it is determined that the dividing element surface is grounded. Further, among the nodes N1 to N4 included in the dividing element surface 20, if the coordinates of one node in the vertical direction are not equal to the coordinates in the vertical direction of the road surface 30, it is determined that the dividing element surface 20 is not grounded. Good.
[0055]
FIG. 8 is an explanatory diagram showing an example of a physical quantity prediction device relating to tire wear according to Embodiment 1 of the present invention. The physical quantity prediction device 50 includes a processing unit 52 and a storage unit 54. Further, an input / output device 51 is connected to the physical quantity prediction device 50, and the reinforcing cord 1 of the tire 10, physical property values of rubber, boundary conditions, and the like are input to the processing unit 52 and the input unit 53 provided therein. Input to the storage unit 54. Here, an input device such as a keyboard and a mouse can be used for the input means 53.
[0056]
The storage unit 54 stores an FEM program in which the prediction method of the present invention for realizing the physical quantity prediction method related to tire wear according to the present invention is incorporated. Here, the storage unit 54 is a hard disk device, a magneto-optical disk device, a non-volatile memory such as a flash memory (a storage medium that can be read only such as a CD-ROM), or a RAM (Random Access Memory). Such a volatile memory or a combination thereof can be used.
[0057]
Moreover, the said program may be what can implement | achieve the prediction method of the physical quantity regarding the abrasion of the tire which concerns on this invention by the combination with the program already recorded on the computer system. Further, the present invention is realized by recording the above-described program for realizing the function of the processing unit 52 in FIG. 8 on a computer-readable recording medium, causing the computer system to read and execute the program recorded on the recording medium. You may perform the prediction method concerning. The “computer system” here includes an OS and hardware such as peripheral devices.
[0058]
The processing unit 52 includes a memory and a CPU. At the time of physical quantity prediction, based on the set tire model and input data, the processing unit 52 reads the program into a memory incorporated in the processing unit 52 and performs calculation. At that time, the processing unit 52 appropriately stores the numerical value in the middle of the calculation in the storage unit 54, and advances the calculation by taking out the stored numerical value. The processing unit 52 may be realized by dedicated hardware instead of the program. The prediction result is displayed on the display means 55 of the input / output device. Here, a CRT (Cathode Ray Tube), a liquid crystal display device or the like can be used for the display means 55. The prediction result can also be output to a printer (not shown) provided as necessary. The storage unit 54 may be built in the processing unit 52 or may be in another device (database server). As described above, the physical quantity prediction device 50 for the structure may access the processing unit 52 and the storage unit 54 by communication from a terminal device (not shown) including the input / output device 51.
[0059]
Next, a specific procedure of the physical quantity prediction method related to tire wear according to the present invention will be described. FIG. 9 is a flowchart showing a procedure example of the physical quantity prediction method relating to tire wear according to Embodiment 1 of the present invention. In the physical quantity prediction method related to tire wear according to the present invention, first, a tire model that is a friction energy prediction target is created (step S101). Specifically, taking account of the composite material layer composed of the reinforcing cord 1 and the rubber 2, the elastic modulus change of the reinforcing cord, or the viscoelastic characteristics of the reinforcing cord 1, it is easy to handle with an analysis method such as FEM or BEM. Simplify the actual tire into shape. When creating a tire model, physical property values such as the reinforcement cord 1 and the rubber 2 are input simultaneously.
[0060]
Next, the boundary condition of the tire 10 is input (step S102). In the boundary condition of the tire 10, a use condition of the tire 10 such as tire air pressure, load, slip angle, speed, or wheel rim width is input. Since all the conditions necessary for the performance prediction of the tire 10 are input in step S101 and step S102, based on the created tire model and the input data, an analysis method such as FEM or BEM can be used in the running state of the tire model. The behavior is analyzed (step S103).
[0061]
Here, when the FEM is used, the entire tire 10 is divided into a finite number of microelements 10a, 10b, etc., for example, by a mesh as shown in FIG. It is read as dimension coordinate data. From the coordinate data of the minute elements and the input data, the strain and stress of each minute element at each time are calculated, and the behavior of the entire tire 10 is obtained. In FEM, it is assumed that the stress and strain in each microelement are constant, and the stress state of each microelement is calculated based on this assumption.
[0062]
Next, information required for prediction of friction energy is acquired from the behavior analysis result in the running state of the tire model (step S104). This information is information such as the coordinates of the nodes constituting the tire 10, displacement, contact force, displacement of the road surface, or data of nodes related to the ground contact surface. Based on these pieces of information, the friction energy per unit area is predicted (step S105). Next, a procedure for predicting frictional energy per unit area will be described.
[0063]
FIG. 10 is a flowchart showing a procedure for predicting the friction energy per unit area. In this prediction procedure, based on the initial shape of the tire, the surface of the divided element that contacts the ground and the surface of the divided element that does not contact are determined before calculating the area shared by each node. First, it is determined whether or not the surface of the split element to be grounded (step S201). Here, this determination procedure will be described. FIG. 11 is a flowchart showing a procedure for determining whether or not the surface of the split element to be grounded.
[0064]
In this contact determination, first, node coordinates and division element surface information (for example, nodes and division shapes constituting the division element surface) are acquired (step S301). Next, it is determined whether or not the division element surface is a division element surface of the tire surface (step S302). When the surface of the division element is a tire surface, it is determined whether or not the surface of the division element is in contact with the ground, and the result is stored in, for example, the storage unit 54 (see FIG. 8) of the physical quantity prediction device 50 (Step 8). S303). When this procedure is executed by a computer, the determination result is stored on a RAM or a predetermined storage area of a hard disk. For example, a method based on an angle formed by the above-described reference vector and a normal vector of the surface of the dividing element can be used to determine whether or not the surface of the dividing element is grounded. It is determined whether or not the processing has been completed for all the divided element surfaces existing on the tire surface (step S304). When the processing has been completed for all the divided element surfaces, the contact determination procedure is completed.
[0065]
Next, referring back to FIG. When the ground contact determination procedure for the split element surface is completed, it is determined whether the node of the split element is touched (step S202). In this determination method, for example, when the vertical contact force between the node and the road surface exceeds 0, it is determined that the node is grounded. When it is determined that the determination target node is grounded, an area shared by the node (hereinafter referred to as a node area) is calculated (step S203). The procedure for calculating the nodal area will be described later.
[0066]
Next, the shear contact force Fs is calculated from the contact force component parallel to the road surface (step S204), and the difference between the node displacement Dn and the road surface displacement Dr between the time step and the immediately preceding time step. A slip amount L between the ground contact surface 11p and the road surface is calculated from Dn−Dr (step S205). The order of steps S203 to S205 does not matter. From the shear contact force Fs, the slip amount L, and the nodal area A, the friction energy E per unit area is calculated based on the equation (1) (step S206). This value is added to the total friction energy per unit area after the node starts grounding (step S207).
E = Fs × L / A (1)
Here, A is the nodal area.
[0067]
It is determined whether or not all nodes have been processed (step S208). If all nodes have not been processed, steps S202 to S207 are repeated. If all nodes are processed in this step, it is determined whether or not all steps have been processed (step S209). If all steps have not been processed, steps S202 to S208 are repeated to process all steps. Thus, the friction energy E prediction procedure is completed.
[0068]
Next, the procedure for calculating the nodal area will be described. FIG. 12 is a flowchart showing the procedure for calculating the nodal area. First, in order to calculate the nodal area, the division element surface information (node coordinates constituting the division element surface, the division shape of the division element surface, and the like) including the node whose area is to be calculated is acquired (step S401). In addition, since the node information inside a tread is unnecessary when calculating | requiring the friction energy E, only the node which exists in the surface will be the object of node area calculation. Next, the determination result as to whether or not the surface of the divided element to be grounded, which is stored in the storage means in the above-described procedure for determining the surface of the divided element, is acquired (step S402). In addition, the order of step S401 and S402 is not ask | required.
[0069]
In the procedure for calculating the node area, the surface of the split element to be grounded is the target, and therefore it is determined whether or not the surface of the split element including the node is grounded (step S403). This can be determined based on the determination result obtained in step S402 as to whether or not the surface of the split element is to be grounded.
[0070]
In the case of the surface of the divided element to be grounded, the node area of the node that is the area calculation target is calculated from the area of the surface of the divided element and the divided shape, and is added to the obtained node area (step S404). Then, the node area of the area calculation target is calculated for all the divided element surfaces including the node of the area calculation target that is in contact with the ground, and when it is added to the already obtained node area, the calculation of the area of the node ends. If the processing of all the divided element surfaces has not been completed, steps S403 and S404 are repeated to process all the divided element surfaces (step S405), and the calculation of the total nodal area is completed.
[0071]
In the physical quantity prediction method relating to tire wear according to the first embodiment, the contact of the surface of the dividing element is determined before calculating the nodal area based on the initial shape of the tire, and the frictional energy per unit area is calculated. Thereby, since it is not necessary to determine the grounding state of each division | segmentation element surface for every time step, calculation amount can be reduced. As a result, in computer simulation, the time required for calculation can be shortened, and effective evaluation can be performed even when hardware resources are limited.
[0072]
(Modification)
Next, a procedure of a physical quantity prediction method related to tire wear according to a modification of the first embodiment will be described. This prediction procedure is substantially the same as the above prediction procedure, but in the procedure for predicting the frictional energy per unit area, the split element surface that contacts the ground based on the deformed shape at the time of rolling of the tire 10 and the split element surface that does not touch the ground Is different for each time step. The other configuration is the same as the above prediction procedure, and thus the description thereof is omitted.
[0073]
In the physical quantity prediction method relating to tire wear according to this modification, the series of flow from the creation of the tire model to the prediction of the friction energy per unit area is the same as described in the first embodiment, so that the description thereof is omitted ( (See FIG. 9). FIG. 13 is a flowchart illustrating a procedure for predicting the friction energy per unit area according to the modification of the first embodiment. Note that steps S501 to S508 of the present prediction procedure are the same as steps S202 to S209 (see FIG. 10) of the procedure for predicting the friction energy per unit area according to the first embodiment, and thus description thereof is omitted.
[0074]
The procedure for calculating the nodal area (step S502) according to the present modification is as follows. FIG. 14 is a flowchart showing the procedure for calculating the nodal area. First, in order to calculate the nodal area, the division element surface information (node coordinates constituting the division element surface, the division shape of the division element surface, and the like) including the node whose area is to be calculated is acquired (step S601). Next, it is determined at each time step whether or not the surface of the dividing element including the node is grounded (step S602). This is because the surface of the divided element to be grounded is a target in the procedure for calculating the nodal area. For this determination method, the above-described method can be used. For example, a method based on an angle formed by the normal vector of the surface of the dividing element and the normal vector of the road surface and a predetermined ground contact reference angle can be used.
[0075]
In the case of the surface of the divided element to be grounded, the node area of the node that is the area calculation target is calculated from the area and the divided shape of the surface of the divided element, and is added to the already obtained node area (step S603). Then, the node area of the area calculation target is calculated for all the divided element surfaces including the node of the area calculation target that is in contact with the ground, and when it is added to the already obtained node area, the calculation of the area of the node ends. If the processing of all the divided element surfaces has not been completed, steps S602 and S603 are repeated to process the entire divided element surfaces (step S604), and the calculation of the total nodal area is completed. The frictional energy per unit area is calculated based on the nodal area at each time step determined in this way.
[0076]
In the physical quantity prediction method relating to tire wear according to this modification, the surface of a split element that contacts the ground is determined at each time step, and the frictional energy per unit area is calculated. Thereby, since a contact surface can be judged according to the deformation | transformation shape at the time of rolling of a tire, abrasion of a tire contact surface can be estimated with higher precision.
[0077]
In the physical quantity prediction method for tire wear according to the first embodiment of the present invention, a program prepared in advance is executed by a computer such as a personal computer or a workstation without using the physical quantity prediction apparatus 50 for the structure. Can be realized. This program can be distributed via a network such as the Internet. The program can also be executed by being recorded on a computer-readable recording medium such as a hard disk, a flexible disk (FD), a CD-ROM, an MO, and a DVD and being read from the recording medium by the computer. The same applies to the following embodiments.
[0078]
(Example)
Here, using the method for predicting wear energy described in the first embodiment, the vertical contact force and the friction energy among physical quantities related to tire wear were evaluated. A radial tire for a passenger car with a tire size of 215 / 45R17 was modeled, and the vertical contact force and friction energy were predicted by FEM. 15 is an explanatory diagram showing the vertical contact force distribution of Comparative Example 1, FIG. 16 is an explanatory diagram showing the vertical contact force distribution of Example 1, and FIG. 17 is an explanatory diagram showing the vertical contact force distribution of Example 2. It is. 18 is an explanatory diagram showing the friction energy distribution of Comparative Example 2, FIG. 19 is an explanatory diagram showing the friction energy distribution of Example 3, and FIG. 20 is an explanatory diagram showing the friction energy distribution of Example 4. is there.
[0079]
In Examples 1 and 2 and Comparative Example 1, the vertical contact force when the inner pressure was 230 kPa and the load was 3.8 kPa was evaluated. In Examples 3 and 4 and Comparative Example 2, the frictional energy was evaluated at an internal pressure of 230 kPa, a load of 3.8 kN, a speed of 36 km / h, a camber angle of 1 degree, and a turning force of 0.2 G. Examples 1 and 3 are based on vertical contact force or frictional energy per unit area, but are calculated including the groove wall 11s. Examples 2 and 4 are based on the vertical contact force or friction energy per unit area, and the groove wall 11s is excluded in the calculation of the nodal area. Comparative examples 1 and 2 are not converted into vertical contact force or friction energy per unit area.
[0080]
In Comparative Examples 1 and 2, since the vertical contact force or frictional energy per unit area is not used, it can be seen that unevenness occurs in the distribution of the vertical contact force or frictional energy (B in FIGS. 15 and 18). Area). In addition, as can be seen from FIGS. 15 to 17, in Examples 1 and 3, the vertical contact force per unit area or the frictional energy is calculated, so in the same place as Comparative Examples 1 and 2, from these It can also be seen that the vertical contact force or friction energy is uniformly distributed. However, since these are calculated including the groove wall 11s, the value of the block end portion in contact with the groove 18 is small.
[0081]
On the other hand, in Examples 2 and 4, the vertical contact force or friction energy per unit area is converted, and the groove wall 11s is excluded in the calculation of the nodal area. There is no unevenness in the distribution of contact force or friction energy, and an appropriate value is also shown at the block end in contact with the groove 18. As described above, according to the method for predicting frictional energy according to the present invention, it is possible to appropriately evaluate the frictional energy not only on the surface of the tire 10 in contact with the ground but also on the block end in contact with the groove 18. Thereby, since wear and uneven wear of the tire 10 can be predicted with high accuracy, the tire 10 can be effectively used for the actual design of the tire 10.
[0082]
(Embodiment 2)
In the second embodiment, a method for predicting uneven tire wear and tire wear using the friction energy prediction method described in the first embodiment will be described. FIG. 21 is a flowchart showing a procedure for predicting uneven wear of the tire according to Embodiment 2 of the present invention. First, friction energies E1, E2,... En are predicted under actual driving conditions such as free rolling, braking, driving, and turning (step S701). The procedure for predicting the frictional energy is as described in the first embodiment.
[0083]
Next, based on how often or how large these operating conditions occur, the respective conditions are weighted (step S702). In this weighting, the weighting coefficients C1, C2,... Cn are determined to have a predetermined size according to the frequency of occurrence of each operating condition. Then, each weighting coefficient is multiplied by the friction energy under each operating condition, and the sum C1 × E1 + C2 × E2 +... + Cn × En is obtained. This is the friction energy Ea during actual traveling (step S703). Based on this frictional energy Ea, uneven wear of the tire during actual vehicle traveling is predicted (step S704).
[0084]
FIG. 22 is a flowchart showing a procedure for predicting uneven wear and wear of the tire according to Embodiment 2 of the present invention. First, the rubber friction index G is obtained by a normal lambourne wear test specified in, for example, JIS K6246 (step S701 ′). Next, under actual driving conditions such as free rolling, braking, and driving, friction energies E1, E2,. (Step S703 ′), the friction energy Ea during actual vehicle travel is predicted (Step S704 ′). Since these procedures are as described above, description thereof is omitted. Based on the friction energy Ea thus determined, the partial wear and wear of the tire during actual vehicle travel are predicted (step S705 ′). For example, the product G × Ea of the rubber friction index G and Ea is obtained, and the reference tire G × Ea _b The wear rate can be predicted by comparing with.
[0085]
Thus, if the rubber friction index G is used, not only the friction energy distribution but also the friction speed can be predicted, so that it can be applied to a wide range of tire wear evaluation. Note that the physical quantity prediction method related to tire wear described in the second embodiment can be executed by a computer as in the first embodiment.
[0086]
【The invention's effect】
As described above, the physical quantity prediction method related to tire wear according to the present invention. ,or Is a physical quantity prediction program related to tire wear according to the present invention. In The tire wear is evaluated by converting the physical quantity related to tire wear into a physical quantity per unit area according to the area of the split element surface shared by the nodes included in the surfaces of the multiple split elements constituting the tire model. . For this reason, physical quantities such as shear contact force and frictional energy are not affected by element division or node position, so that wear on the entire contact surface can be accurately predicted. Further, since it is not necessary to measure the friction energy, it is possible to efficiently predict tire wear and improve the development efficiency.
[0087]
Also, a method for predicting physical quantities related to tire wear according to the present invention. Law, Or a prediction program for physical quantities related to tire wear according to the present invention. In Is to obtain the area shared by each node based on the division shape of the surface of the division element. This makes it possible to accurately calculate the area shared by the nodes regardless of the shape of the surface of the dividing element and accurately predict the physical quantity related to tire wear.
[0088]
Also, a method for predicting physical quantities related to tire wear according to the present invention. Law, Or a prediction program for physical quantities related to tire wear according to the present invention. In Includes a step of determining whether or not the surface of the dividing element is grounded. Thereby, for example, the surface of the dividing element that is not grounded, such as the surface of the dividing element in the groove wall or the like, can be excluded when the physical quantity is predicted. As a result, the frictional energy at the block end in contact with the groove can be appropriately evaluated, so that wear at such a portion can be accurately predicted. The frictional energy can be predicted with high accuracy by determining the ground contact state of the surface of the dividing element even at the end of the block where it has been extremely difficult to measure the frictional energy.
[0089]
Also, a method for predicting physical quantities related to tire wear according to the present invention. In law Is based on the initial shape of the tire, and before calculating the area shared by each node, it is determined in advance whether or not the surface of the dividing element is in contact with the ground. For this reason, whether or not the surface of the dividing element is in contact with the ground is determined only once based on the initial shape of the tire. In particular, when the prediction method according to the present invention is executed using a computer, there is an advantage that hardware resources can be saved.
[0090]
Also, a method for predicting physical quantities related to tire wear according to the present invention. In law Since it is determined at each moment whether or not the surface of the dividing element is in contact with the tire according to the deformed shape of the tire, it is possible to further improve the prediction accuracy of the physical quantity related to tire wear.
[0091]
Also, a method for predicting physical quantities related to tire wear according to the present invention. Law, Or a prediction program for physical quantities related to tire wear according to the present invention. In Since the wear of the tire is predicted using the friction energy per unit area, the wear can be predicted with high accuracy. This makes it possible to obtain effective information for the tire design. Further, since it is not necessary to measure frictional energy, the labor for predicting tire wear can be greatly reduced, and the tire development period can be shortened.
[0092]
Also, a method for predicting physical quantities related to tire wear according to the present invention. In law Since the rubber friction index is used, not only the wear distribution but also the actual wear rate can be predicted, and a wider evaluation of tire wear can be performed.
[0093]
Also, a method for predicting physical quantities related to tire wear according to the present invention. In law Is provided with processing means for executing the physical quantity prediction method relating to the tire wear. As a result, physical quantities such as shear contact force and friction energy can be obtained without being affected by element division and node position, so that wear on the entire contact surface can be accurately predicted. Further, since it is not necessary to measure the friction energy, it is possible to efficiently predict tire wear and improve the development efficiency.
[Brief description of the drawings]
FIG. 1 is a partial cross-sectional view showing a cross section of a tire to be predicted cut along a meridian plane including a rotation axis thereof.
FIG. 2 is a perspective view showing an example in which a tire to be predicted is divided into minute elements.
FIG. 3 is an explanatory diagram showing an example in which a ground surface on which a cap tread is in contact with the ground is divided into minute elements.
FIG. 4 is a plan view showing a ground plane divided into nine by rectangular elements.
FIG. 5 is a perspective view showing one of the dividing elements including a contact surface of the tire model.
FIG. 6 is an explanatory diagram showing an example of a method for determining the grounding of the surface of the dividing element.
FIG. 7 is an explanatory diagram showing another example of a method for determining the grounding of the surface of a split element.
FIG. 8 is an explanatory diagram showing an example of a physical quantity prediction device related to tire wear according to Embodiment 1 of the present invention;
FIG. 9 is a flowchart showing a procedure example of a physical quantity prediction method relating to tire wear according to Embodiment 1 of the present invention;
FIG. 10 is a flowchart showing a procedure for predicting friction energy per unit area.
FIG. 11 is a flowchart showing a procedure for determining whether or not the surface of a split element to be grounded.
FIG. 12 is a flowchart showing a procedure for calculating a nodal area.
FIG. 13 is a flowchart showing a procedure for predicting frictional energy per unit area according to a modification of the first embodiment.
FIG. 14 is a flowchart showing a procedure for calculating a nodal area.
15 is an explanatory diagram showing a vertical contact force distribution of Comparative Example 1. FIG.
16 is an explanatory diagram showing a vertical contact force distribution of Example 1. FIG.
17 is an explanatory diagram showing a vertical contact force distribution of Example 2. FIG.
18 is an explanatory diagram showing a friction energy distribution of Comparative Example 2. FIG.
FIG. 19 is an explanatory diagram showing the friction energy distribution of Example 3.
20 is an explanatory diagram showing a friction energy distribution of Example 4. FIG.
FIG. 21 is a flowchart showing a procedure for predicting uneven wear of a tire according to Embodiment 2 of the present invention.
FIG. 22 is a flowchart showing a procedure for predicting uneven wear and wear of a tire according to Embodiment 2 of the present invention.
[Explanation of symbols]
10 tires
11 Cap tread
11B block
11p Ground plane
11s groove wall
18 groove
20, 21 Split element surface
30 road surface
50 Physical quantity prediction device
52 processor
53 Input means
54 Memory

Claims

In predicting tire wear using the finite element method, boundary element method, and other analysis methods,
A procedure for creating a tire model by dividing a tire into a plurality of division elements;
Anda procedure for determining whether a node present on the surface of the parting element is grounded,
When the node is grounded, the area of the surface of the dividing element that is not grounded is excluded, the procedure for calculating the area of the surface of the dividing element shared by the node , and the calculated node is A procedure for converting a physical quantity related to tire wear into a physical quantity per unit area according to the shared area , and
When the nodes are not grounded, a procedure for determining whether or not the processing for all the nodes has been completed is executed.
Physical quantity of the prediction method for the wear of the tire characterized by a crotch.

In the procedure for calculating the area of the surface of the dividing element shared by the nodes,
A procedure for calculating the surface area of the dividing element including the node;
A procedure for determining a ratio of the nodes sharing the area of the surface of the dividing element based on the division shape of the surface of the dividing element including the node;
The area in which the nodes are sharing the area of the surface of the partition element, a step of adding the surface of all the dividing elements including the node,
The method for predicting a physical quantity related to wear of a tire according to claim 1, comprising:

3. The method according to claim 1 , wherein it is determined at each time step whether or not a node existing on the surface of the dividing element is grounded based on a deformed shape of the tire model when rolling. Of predicting physical quantities related to tire wear.

In the procedure determines whether the node that is present on the surface of the dividing element is grounded, based on the initial shape of the upper Symbol tires, that the surface of the dividing element to determine whether or not ground The method for predicting a physical quantity related to wear of a tire according to claim 1 or 2, characterized in that:

4. The physical quantity prediction method for tire wear according to claim 3, further comprising determining whether the surface of the dividing element is in contact with the ground according to a deformed shape of the tire.

According to the method for predicting a physical quantity related to tire wear according to any one of claims 1 to 5 under a free rolling, braking, driving, turning and other driving conditions, A procedure for predicting the friction energy of
Considering the frequency of occurrence of each operating condition, a procedure for weighting each operating condition;
From the friction energy per unit area under each driving condition and the weighting of each driving condition, a procedure for predicting the friction energy per unit area during actual vehicle travel;
A method for predicting a physical quantity related to wear of a tire, comprising:

The tire wear according to claim 6, further comprising a step of measuring a rubber wear index of rubber constituting the surface of the tire, and predicting the wear of the tire during running of the vehicle based on the rubber wear index. Prediction method for physical quantities.

Processing means for processing each procedure in the physical quantity prediction method for tire wear according to any one of claims 1 to 7,
Input means for giving physical properties of the material constituting the tire, operating conditions, boundary conditions and other data to the processing means;
Display means for displaying a prediction result by the processing means;
An apparatus for predicting a physical quantity related to tire wear.

In predicting physical quantities related to tire wear using the finite element method, boundary element method, and other analysis methods,
The procedure to create a tire model by dividing the tire into multiple elements,
Determining whether or not a node constituting the surface of the dividing element is grounded ,
When the node is grounded, the area of the surface of the dividing element that is not grounded is excluded, the procedure for calculating the area of the surface of the dividing element shared by the node , and the calculated node is A procedure for converting a physical quantity related to tire wear into a physical quantity per unit area according to the shared area , and
When the nodes are not grounded, a procedure for determining whether or not the processing for all the nodes has been completed is executed.
Physical quantity prediction program regarding tire wear, characterized in that to execute a call to the computer.

In the procedure for calculating the area of the surface of the dividing element shared by the nodes,
A procedure for calculating the surface area of the dividing element including the node;
A procedure for determining a ratio of the nodes sharing the area of the surface of the dividing element based on the division shape of the surface of the dividing element including the node;
Adding the area where the nodes share the area of the surface of the dividing element for the surfaces of all the dividing elements including the nodes;
10. The physical quantity prediction program relating to tire wear according to claim 9, wherein the computer is executed.

The computer is configured to determine at each time step whether or not a node existing on the surface of the dividing element is grounded based on a deformed shape of the tire model when rolling. Item 11. The physical quantity prediction program relating to tire wear according to Item 9 or 10.

According to the method for predicting a physical quantity related to tire wear according to any one of claims 1 to 5 under a free rolling, braking, driving, turning and other driving conditions, A procedure for predicting the friction energy of
Considering the frequency of occurrence of each operating condition, a procedure for weighting each operating condition,
From the friction energy per unit area under each driving condition and the weighting of each driving condition, a procedure for predicting the friction energy per unit area during actual vehicle travel;
A computer program for predicting physical quantities related to tire wear.