(57)【要約】【課題】粘弾性材料を含む回転体の転動状態における耐久性を予測できるようにすること。【解決手段】粘弾性材料を含む回転体の回転体モデルの幾何形状や応力等の解析を実行する（ステップＳ１０１）。次に、前記解析によって得られた数値データから、発熱エネルギ関連の特性を評価するために必要な物理量を抽出する（ステップＳ１０２）。次に、抽出した物理量と、発熱を予測するのに必要変数を用いて、回転体モデルの発熱エネルギ関連の特性を評価する（ステップＳ１０３）。そして、この発熱エネルギ関連の特性と、回転体モデルの転動の状態を表す変数等を用いて回転体モデルの耐久性を予測する（ステップＳ１０４）。【選択図】        図１
2005306113.tif 000002
【特許請求の範囲】
【請求項１】
回転体の形状と当該回転体を形成する材料の特性とにより回転体モデルを作成するモデル作成手順と、
前記回転体モデルから解析可能な数値データを作成する解析手順と、
前記回転体モデルの転動時における物理量を前記解析手順によって得られた結果から抽出する抽出手順と、
前記抽出手順で抽出した前記物理量と発熱を予測するための変数とにより前記回転体モデルの発熱エネルギ関連の特性を評価する評価手順と、
前記評価手順で得られた発熱エネルギ関連の特性を用いて回転体の耐久性を予測する耐久性予測手順と、
を含むことを特徴とする回転体の耐久性予測方法。
【請求項２】
前記物理量は応力及びひずみを含んでおり、
前記応力及び前記ひずみから局所座標を参照した応力、ひずみに変換する演算を行い前記回転体モデルの一点における応力、ひずみの成分を演算し、
前記回転体モデルの円周方向に隣接する諸点の応力、ひずみを順次演算し、当該回転体モデルの一周分の応力、ひずみを得ることにより応力、ひずみの変化特性の曲線を導出し、一周分の応力、ひずみの有限次数のフーリエ級数展開の演算を行い、各フーリエ次数ごとに曲線の振幅、位相を演算し、
前記回転体モデルを形成する前記材料の損失係数に応じた位相遅れをひずみ値に与えての、各フーリエ次数ごとのヒステリシスループの面積の演算に基づきフーリエ次数とヒステリシスループ面積の積の総和を導出し、
以上の一連の演算過程を応力、ひずみの全成分について反復実行し各成分ごとの該総和を演算することによりそれぞれの位置における発熱エネルギ密度を導出し、
当該位置の発熱エネルギ密度にて発熱エネルギ密度特性が近似される、当該位置を含む領域の体積と発熱エネルギ密度の積を演算することにより、当該領域における発熱エネルギを導出し、
以上の一連の演算および導出の過程を回転体の全体について反復実行し回転体全体の発熱エネルギを導出する方法を用いて、前記発熱エネルギ関連の特性を評価することを特徴とする請求項１に記載の回転体の耐久性予測方法。
【請求項３】
前記発熱エネルギ関連の特性を、前記回転体モデルのひずみエネルギ密度の振幅から算出する方法を用いて評価することを特徴とする請求項１に記載の回転体の耐久性予測方法。
【請求項４】
前記解析手順は、前記回転体モデルを数値解析法で取り扱い可能な有数個の要素として分割することを特徴とする請求項１〜３のいずれか１項に記載の回転体の耐久性予測方法。
【請求項５】
前記回転体モデルが転動している状態、或いは、前記回転体モデルの進行方向と当該回転体モデルの回転の中心軸とのなす角度が保持された状態で当該回転体モデルが転動している状態を数値解析することを特徴とする請求項１〜４のいずれか１項に記載の回転体の耐久性予測方法。
【請求項６】
前記解析手順は、前記回転体モデルの回転の中心軸と当該回転体モデルとが一体となって連動する設定を行い、当該中心軸に速度を定義することで、前記回転体モデルの転動する速度を設定することを特徴とする請求項１〜５のいずれか１項に記載の回転体の耐久性予測方法。
【請求項７】
前記解析手順は、前記回転体モデルに接触している接触体を定義し、当該接触体の移動速度を定義することにより、前記回転体モデルの転動する速度を設定することを特徴とする請求項１〜５のいずれか１項に記載の回転体の耐久性予測方法。
【請求項８】
前記解析手順は、前記回転体モデルの移動方向を定義し、当該移動方向と前記中心軸とのなす角度に変化を定義することで、当該回転体モデルの移動方向の角度を設定することを特徴とする請求項１〜７のいずれか１項に記載の回転体の耐久性予測方法。
【請求項９】
前記解析手順は、前記回転体モデルの移動方向を定義し、当該移動方向と前記中心軸とのなす角度に変化を定義することで、当該回転体モデルの接触体に対する傾き角度を設定することを特徴とする請求項１〜８のいずれか１項に記載の回転体の耐久性予測方法。
【請求項１０】
前記解析手順は、前記接触体に角度変化を定義することで、前記回転体モデルの移動方向の角度、又は前記回転体モデルの接触体に対する傾き角度の少なくとも一方を設定することを特徴とする請求項１〜９のいずれか１項に記載の回転体の耐久性予測方法。
【請求項１１】
前記回転体はタイヤであることを特徴とする請求項１〜１０のいずれか１項に記載の回転体の耐久性予測方法。
【請求項１２】
請求項１〜１１のいずれか１項に記載の回転体の耐久性予測方法をコンピュータに実行させることを特徴とする耐久性予測用コンピュータプログラム。
【請求項１３】
回転体の形状と当該回転体を形成する材料の特性とにより作成される回転体モデルから解析可能な数値データを作成する回転体解析部と、
前記回転体モデルの物理量を前記回転体解析部によって得られた結果から抽出する物理量抽出部と、
前記物理量抽出部で抽出された前記物理量と発熱を予測するための変数とにより前記回転体モデルの発熱エネルギ関連の特性を評価する発熱エネルギ評価部と、
前記発熱エネルギ評価部で得られた発熱エネルギ関連の特性を用いて回転体の耐久性を予測する耐久性予測部と、
を備えることを特徴とする回転体の耐久性予測装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、構成材料に粘弾性材料を含む回転体の耐久性を評価する技術に関し、さらに詳しくは、このような回転体の転動状態における耐久性を評価できる回転体の耐久性予測方法及び回転体の耐久性予測用コンピュータプログラム、並びに回転体の耐久性予測装置に関する。
【背景技術】
【０００２】
従来の回転体の耐久性能の評価は、回転体を実際に試作し、試験を実施する等の方法で行われていたため、多くの時間、費用、及び労力が必要であったが、近年におけるコンピュータ技術の発達により、実際に回転体を試作せずに、コンピュータ上の数値シミュレーションで回転体の性能や耐久性を予測する技術が用いられるようになってきている。このようなシミュレーションでは、例えば回転体がタイヤである場合には、有限要素法などの数値解析によってタイヤのひずみエネルギを求め、ひずみエネルギからタイヤの発熱を予測し、このタイヤの発熱と走行距離とにより、転動抵抗を予測して性能を評価している。
【０００３】
しかし、タイヤは複数の異なる材料からなる部材によって形成されているため、このような手法では正確なひずみエネルギを得ることができず、また、タイヤの各部位ごとに異なる変形周波数を考慮することができない、などの理由により忠実度の高い性能の評価を行うことが困難であった。そこで、かかる問題を解決するため、特許文献１には、応力、ひずみのヒステリシスループの面積から発熱エネルギ及びタイヤ転動抵抗を推定する方法及び装置が開示されている。この特許文献１では、応力及びひずみの成分ごとに材料の損失係数を変えて取り扱うことができ、また、タイヤの各変形周波数ごとにヒステリシスループが計算されるので、タイヤの各部位ごとに異なる変形周波数を考慮することができるので、忠実度の高い性能の評価を行うことができる。
【０００４】
【特許文献１】特開平１１−２３７３３２号公報
【発明の開示】
【発明が解決しようとする課題】
【０００５】
しかし、上記の方法及び装置は、静止している状態の数値解析結果から発熱を予測していたため、コーナーリングや制駆動時等の動的な状態の発熱分布を再現していなかった。さらに、回転体の耐久性を予測するためには様々な動的条件（例えばタイヤであれば、スリップ角をパラメトリックに変更する場合など）を考慮に入れた数値解析を行う必要があるが、上述の方法及び装置では不可能であった。
【０００６】
本発明は、上記に鑑みてなされたものであって、構成材料に粘弾性材料を含む回転体の転動状態における耐久性を評価できる回転体の耐久性予測方法及び回転体の耐久性予測用コンピュータプログラム、並びに回転体の耐久性予測装置を提供することを目的とする。
【課題を解決するための手段】
【０００７】
上述した課題を解決し、目的を達成するために、本発明に係る回転体の耐久性予測方法は、回転体の形状と当該回転体を形成する材料の特性とにより回転体モデルを作成するモデル作成手順と、前記回転体モデルから解析可能な数値データを作成する解析手順と、前記回転体モデルの転動時における物理量を前記解析手順によって得られた結果から抽出する抽出手順と、前記抽出手順で抽出した前記物理量と発熱を予測するための変数とにより前記回転体モデルの発熱エネルギ関連の特性を評価する評価手順と、前記評価手順で得られた発熱エネルギ関連の特性を用いて回転体の耐久性を予測する耐久性予測手順と、を含むことを特徴とする。
【０００８】
この回転体の耐久性予測方法は、粘弾性材料を含む回転体の耐久性を評価するにあたり、転動状態、即ち、定常回転も含む動的な状態の回転体モデルの幾何形状、材料特性、及び解析目的に対して必要な条件等を解析可能な数値データとして解析手順で作成している。抽出手段で前記数値データから転動状態の前記回転体モデルの物理量、つまり、応力及びひずみを抽出する。この物理量と前記回転体を形成する粘弾性材料の特性等から、前記回転体モデルの発熱エネルギ関連の特性を評価手順によって評価する。この評価手順によって得られた発熱エネルギ関連の特性と回転体の転動の状態を表す変数等を用いて、発熱と耐久性との相関関係より耐久性を予測する耐久性予測手段によって、回転体モデルの耐久性を予測する。この結果、構成材料に粘弾性材料を含む回転体の転動状態における耐久性を評価できる。また、本発明に係る回転体の耐久性予測用コンピュータプログラムによれば、前述の回転体の耐久性予測方法がコンピュータを利用して実現できる。
【０００９】
また、次の発明に係る回転体の耐久性予測方法は、前記回転体の耐久性予測方法において、前記物理量は応力及びひずみを含んでおり、前記応力及び前記ひずみから局所座標を参照した応力、ひずみに変換する演算を行い前記回転体モデルの一点における応力、ひずみの成分を演算し、前記回転体モデルの円周方向に隣接する諸点の応力、ひずみを順次演算し、当該回転体モデルの一周分の応力、ひずみを得ることにより応力、ひずみの変化特性の曲線を導出し、一周分の応力、ひずみの有限次数のフーリエ級数展開の演算を行い、各フーリエ次数ごとに曲線の振幅、位相を演算し、前記回転体モデルを形成する前記材料の損失係数に応じた位相遅れをひずみ値に与えての、各フーリエ次数ごとのヒステリシスループの面積の演算に基づきフーリエ次数とヒステリシスループ面積の積の総和を導出し、以上の一連の演算過程を応力、ひずみの全成分について反復実行し各成分ごとの該総和を演算することによりそれぞれの位置における発熱エネルギ密度を導出し、当該位置の発熱エネルギ密度にて発熱エネルギ密度特性が近似される、当該位置を含む領域の体積と発熱エネルギ密度の積を演算することにより、当該領域における発熱エネルギを導出し、以上の一連の演算および導出の過程を回転体の全体について反復実行し回転体全体の発熱エネルギを導出する方法を用いて、前記発熱エネルギ関連の特性を評価することを特徴とする。
【００１０】
この回転体の耐久性予測方法は、前記物理量の応力とひずみから算出されるヒステリシスループの面積の総和を、応力とひずみの全成分について反復実行して各成分ごとの当該総和を演算し、それぞれの位置における発熱エネルギ密度を導出している。さらに、この発熱エネルギ密度と当該位置の体積から発熱エネルギを導出する。この発熱エネルギから前記耐久性予測手順によって前記回転体モデルの耐久性を予測する。この結果、構成材料に粘弾性材料を含む回転体の転動状態における耐久性を評価できる。また、本発明に係る回転体の耐久性予測用コンピュータプログラムによれば、前述の回転体の耐久性予測方法がコンピュータを利用して実現できる。
【００１１】
また、次の発明に係る回転体の耐久性予測方法は、前記回転体の耐久性予測方法において、前記発熱エネルギ関連の特性を、前記回転体モデルのひずみエネルギ密度の振幅から算出する方法を用いて評価することを特徴とする。
【００１２】
この回転体の耐久性予測方法は、前記回転体モデルを転動させたときのひずみの変化によるひずみエネルギ密度の振幅から発熱エネルギを予測している。この発熱エネルギから前記耐久性予測手順によって前記回転体モデルの耐久性を予測する。この結果、構成材料に粘弾性材料を含む回転体の転動状態における耐久性を評価できる。また、本発明に係る回転体の耐久性予測用コンピュータプログラムによれば、前述の回転体の耐久性予測方法がコンピュータを利用して実現できる。
【００１３】
また、次の発明に係る回転体の耐久性予測方法は、前記回転体の耐久性予測方法において、前記解析手順は、前記回転体モデルを数値解析法で取り扱い可能な有数個の要素として分割することを特徴とする。
【００１４】
この回転体の耐久性予測方法は、前記回転体モデルを有限要素法、有限体積法、差分法、または境界要素法など数値解析法で取り扱い可能な有数個の要素として分割する解析手段によって数値データとして解析している。これにより、より正確に回転体モデルの解析を行うことができる。また、本発明に係る回転体の耐久性予測用コンピュータプログラムによれば、前述の回転体の耐久性予測方法がコンピュータを利用して実現できる。
【００１５】
また、次の発明に係る回転体の耐久性予測方法は、前記回転体の耐久性予測方法において、前記回転体モデルが転動している状態、或いは、前記回転体モデルの進行方向と当該回転体モデルの回転の中心軸とのなす角度が保持された状態で当該回転体モデルが転動している状態を数値解析することを特徴とする。
【００１６】
この回転体の耐久性予測方法は、回転体モデルの進行方向と中心軸とのなす角度がいずれの状態においても回転体モデルが転動している状態の数値解析をすることにより、確実に回転体モデルの転動時の数値解析を行なうことができる。また、本発明に係る回転体の耐久性予測用コンピュータプログラムによれば、前述の回転体の耐久性予測方法がコンピュータを利用して実現できる。
【００１７】
また、次の発明に係る回転体の耐久性予測方法は、前記回転体の耐久性予測方法において、前記解析手順は、前記回転体モデルの回転の中心軸と当該回転体モデルとが一体となって連動する設定を行い、当該中心軸に速度を定義することで、前記回転体モデルの転動する速度を設定することを特徴とする。
【００１８】
この回転体の耐久性予測方法は、上記のように回転体モデルの中心軸と当該回転体モデルとが一体となって連動するように設定することにより、当該中心軸の転動速度を定義するだけで容易に回転体モデルの転動速度を設定できる。また、本発明に係る回転体の耐久性予測用コンピュータプログラムによれば、前述の回転体の耐久性予測方法がコンピュータを利用して実現できる。
【００１９】
また、次の発明に係る回転体の耐久性予測方法は、前記回転体の耐久性予測方法において、前記解析手順は、前記回転体モデルに接触している接触体を定義し、当該接触体の移動速度を定義することにより、前記回転体モデルの転動する速度を設定することを特徴とする。
【００２０】
この回転体の耐久性予測方法は、回転体モデルに接触する接触体を定義し、当該回転体モデルに対する接触体の移動速度を定義することにより、回転体モデルの転動速度を設定しなくても、容易に回転体モデルの転動速度を設定できる。また、本発明に係る回転体の耐久性予測用コンピュータプログラムによれば、前述の回転体の耐久性予測方法がコンピュータを利用して実現できる。
【００２１】
また、次の発明に係る回転体の耐久性予測方法は、前記回転体の耐久性予測方法において、前記解析手順は、前記回転体モデルの移動方向を定義し、当該移動方向と前記中心軸とのなす角度に変化を定義することで、当該回転体モデルの移動方向の角度を設定することを特徴とする。
【００２２】
この回転体の耐久性予測方法は、定義した移動方向と前記回転体モデルの回転の中心軸とがなす角度に変化を定義することにより、回転体モデルの移動方向の角度を設定することができる。回転体は、例えば接触面に対して当該回転体の回転方向と異なる方向に移動することがあるが、上記のように中心軸と移動方向とがなす角度に変化を与えることにより、回転体モデルが転動しながら、転動の方向と異なる方向に移動する状態の評価をすることができる。この結果、回転体の耐久性予測方法を、より実際の回転体の使用状態に近付けることができ、忠実度の高い回転体の耐久性予測方法とすることができる。また、本発明に係る回転体の耐久性予測用コンピュータプログラムによれば、前述の回転体の耐久性予測方法がコンピュータを利用して実現できる。
【００２３】
また、次の発明に係る回転体の耐久性予測方法は、前記回転体の耐久性予測方法において、前記解析手順は、前記回転体モデルの移動方向を定義し、当該移動方向と前記中心軸とのなす角度に変化を定義することで、当該回転体モデルの接触体に対する傾き角度を設定することを特徴とする。
【００２４】
これにより、回転体モデルが転動しながら、転動の方向と異なる方向に移動する状態の評価をすることができる。この結果、回転体の耐久性予測方法を、より実際の回転体の使用状態に近付けることができ、忠実度の高い回転体の耐久性予測方法とすることができる。回転体モデルの接触体に対する傾き角度には、例えばキャンバー角度やトー角度がある。
【００２５】
また、次の発明に係る回転体の耐久性予測方法は、前記回転体の耐久性予測方法において、前記解析手順は、前記接触体に角度変化を定義することで、前記回転体モデルの移動方向の角度、又は前記回転体モデルの接触体に対する傾き角度の少なくとも一方を設定することを特徴とする。
【００２６】
この回転体の耐久性予測方法は、回転体モデルに接触する接触体の、当該回転モデルの回転方向に対する角度変化を定義することにより、回転体モデルの前記接触体に対する移動方向の角度、又は前記回転体モデルの接触体に対する傾き角度の少なくとも一方を設定することができる。この結果、上記と同様に忠実度の高い回転体の耐久性予測方法とすることができる。また、本発明に係る回転体の耐久性予測用コンピュータプログラムによれば、前述の回転体の耐久性予測方法がコンピュータを利用して実現できる。
【００２７】
また、次の発明に係る回転体の耐久性予測方法は、前記回転体の耐久性予測方法において、前記回転体はタイヤであることを特徴とする。
【００２８】
この回転体の耐久性予測方法は、上記のように粘弾性材料を含む回転体の具体例としては、例えば自動車の車輪に使用されるタイヤが挙げられ、タイヤの耐久性を評価することができる。このタイヤは、荷重を受けつつ高速で転動するため、転動状態での耐久性の評価が重要であり、本発明は、タイヤの評価に適した評価を行うことができる。また、本発明に係る回転体の耐久性予測用コンピュータプログラムによれば、前述の回転体の耐久性予測方法がコンピュータを利用して実現できる。
【００２９】
また、次の発明に係る回転体の耐久性予測装置は、回転体の形状と当該回転体を形成する材料の特性とにより作成される回転体モデルから解析可能な数値データを作成する回転体解析部と、前記回転体モデルの物理量を前記回転体解析部によって得られた結果から抽出する物理量抽出部と、前記物理量抽出部で抽出された前記物理量と発熱を予測するための変数とにより前記回転体モデルの発熱エネルギ関連の特性を評価する発熱エネルギ評価部と、前記発熱エネルギ評価部で得られた発熱エネルギ関連の特性を用いて回転体の耐久性を予測する耐久性予測部と、を備えることを特徴とする。
【００３０】
この回転体の耐久性予測装置は、粘弾性材料を含む回転体の耐久性を評価するにあたり、転動状態、即ち、定常回転も含む動的な状態の回転体モデルの幾何形状、材料特性、及び解析目的に対して必要な条件等を、解析可能な数値データとして作成している。この数値データからは、転動状態の前記回転体モデルの物理量、つまり、応力及びひずみを抽出する。この物理量と前記回転体を形成する粘弾性材料の特性等から、前記回転体モデルの発熱エネルギ関連の特性を評価する。この発熱エネルギ関連の特性と回転体の転動の状態を表す変数等を用いて、発熱と耐久性との相関関係より回転体モデルの耐久性を予測する。この結果、構成材料に粘弾性材料を含む回転体の転動状態における耐久性を評価することができる。
【発明の効果】
【００３１】
以上説明したように、本発明にかかる回転体の耐久性予測方法では、構成材料に粘弾性材料を含む回転体の転動状態における耐久性を評価できる、という効果を奏する。また、本発明に係る回転体の耐久性予測用コンピュータプログラムによれば、前述の回転体の耐久性予測方法がコンピュータを利用して実現できる、という効果を奏する。また、本発明に係る回転体の耐久性予測装置では、構成材料に粘弾性材料を含む回転体の転動状態における耐久性を評価できる、という効果を奏する。
【発明を実施するための最良の形態】
【００３２】
以下、この発明につき図面を参照しつつ詳細に説明する。なお、この発明を実施するための最良の形態によりこの発明が限定されるものではない。また、下記発明を実施するための最良の形態における構成要素には、当業者が容易に想定できるもの、あるいは実質的に同一のものが含まれる。以下の説明では、粘弾性材料を含む回転体としてタイヤを例とするが、本発明の適用対象はタイヤに限られるものではない。また、タイヤの種類も空気入りタイヤに限定されるものではなく、本発明の適用できるタイヤの種類は問わない。
【実施例】
【００３３】
本発明に係る回転体の耐久性予測方法及び評価装置は、粘弾性材料を含む回転体の耐久性を評価するにあたって、定常回転を含む動的な状態、即ち転動状態の回転体を解析し、発熱エネルギを求めることによって、前記回転体の耐久性を評価する点に特徴がある。
【００３４】
本発明は、動的な状態の回転体を解析し、発熱エネルギを求めることによって、発熱と耐久性との相関関係より、回転体の実際の使用状況での耐久性を評価できることに着目したものである。例えば、回転体が自動車の車輪に使用するタイヤの場合、自動車のコーナーリング時や、制動時、加速をする際の駆動時など、実際の使用状況での発熱エネルギを求めることができる。この結果、回転体の耐久性予測方法及び評価装置での耐久性の評価を、前記回転体を様々な状況で実際に使用した際の実際の耐久性に近付けることができる。次に、本発明に係る回転体の耐久性予測方法及び耐久性予測装置の具体的な構成について説明する。
【００３５】
図１は、本発明に係る回転体の耐久性予測方法の処理手順を示すフローチャートである。まず、本発明に係る回転体の耐久性予測装置の構成について説明する。図２は、本発明に係る回転体の耐久性予測装置を示す説明図である。回転体の耐久性予測装置（以下耐久性予測装置という）５０は、処理部５０ｐと記憶部５０ｍとで構成される。処理部５０ｐは、さらに回転体解析部５１と、物理量抽出部５２と、発熱エネルギ評価部５３と、耐久性予測部５４とで構成される。また、この耐久性予測装置５０には入出力ポート５６が備えられており、この入出力ポート５６を介して入出力装置６０が接続されている。そして、入出力装置６０に備えられた入力手段６１で耐久性予測対象の回転体であるタイヤのゴムや補強コードの物性値、あるいは境界条件等を処理部５０ｐや記憶部５０ｍへ入力する。ここで、入力手段６１には、キーボード、マウス等の入力デバイスを使用することができる。
【００３６】
処理部５０ｐを構成する回転体解析部５１と、物理量抽出部５２と、発熱エネルギ評価部５３と、耐久性予測部５４とは、入出力ポート５６によって相互に接続されており、相互に評価条件や演算結果その他の情報をやり取りできるように構成される。また、記憶部５０ｍも入出力ポート５６に接続されており、処理部５０ｐの演算結果を一時的に格納する。記憶部５０ｍに対しては、処理部５０ｐを構成する前記回転体解析部５１等が入出力ポート５６を介してアクセスできるようになっている。
【００３７】
記憶部５０ｍには、本発明に係る回転体の耐久性予測方法を実現するコンピュータプログラムが格納されている。ここで、記憶部５０ｍは、ハードディスク装置や光磁気ディスク装置、又はフラッシュメモリ等の不揮発性のメモリ（ＣＤ−ＲＯＭ等のような読み出しのみが可能な記憶媒体）や、ＲＡＭ（Random Access Memory）のような揮発性のメモリ、あるいはこれらの組み合わせにより構成することができる。
【００３８】
また、上記コンピュータプログラムは、コンピュータシステムにすでに記録されているコンピュータプログラムとの組み合わせによって、本発明に係るタイヤの摩耗に関する物理量の予測方法を実現できるものであってもよい。また、図２における処理部５０ｐの機能を実現するための上記コンピュータプログラムをコンピュータ読み取り可能な記録媒体に記録して、この記録媒体に記録されたコンピュータプログラムをコンピュータシステムに読み込ませ、実行することにより本発明に係る予測方法を実行してもよい。なお、ここでいう「コンピュータシステム」とは、ＯＳや周辺機器などのハードウェアを含むものとする。
【００３９】
処理部５０ｐは、メモリ及びＣＰＵにより構成されている。回転体の耐久性を評価する際には、設定されたタイヤモデル及び入力データに基づいて、処理部５０ｐが前記コンピュータプログラムを当該処理部５０ｐに組み込まれたメモリに読み込んで演算する。その際に処理部５０ｐは、適宜記憶部５０ｍへ演算途中の数値を格納し、また格納した数値を取り出して演算を実行する。なお、この処理部５０ｐは、前記コンピュータプログラムの代わりに専用のハードウェアにより実現されるものであってもよい。
【００４０】
耐久性の評価結果は、入出力装置６０の表示手段６３に表示される。ここで、表示手段６３には、ＣＲＴ（Cathode Ray Tube）や液晶表示装置等を使用することができる。また、記憶部５０ｍは、処理部５０ｐに内蔵されるものであっても、他の装置（例えばデータベースサーバ）内にあってもよい。さらに、上記耐久性予測装置５０は、入出力装置６０を備えた端末装置から、通信により処理部５０ｐや記憶部５０ｍにアクセスするものであってもよい。次に、この耐久性予測装置５０を用いて本発明に係る回転体の耐久性予測方法を実現する場合の処理手順について説明する。
【００４１】
本発明に係る回転体の耐久性予測方法を実現するにあたって、まず、本発明に係る耐久性予測装置５０が備える回転体解析部５１は、耐久性予測対象の回転体であるタイヤの幾何形状、材料特性及び解析目的に対し必要な条件等を解析可能な数値データにするよう解析を実行する（ステップＳ１０１、図１）：転動解析。ここで、本発明では、前記解析を実行する際に、ＦＥＭ（Finite Element Method：有限要素法）を使用する。なお、本発明に係るタイヤの物理量の予測方法に適用できる解析手法はＦＥＭに限られず、ＢＥＭ（Boundary Element Method：境界要素法）、ＦＤＭ（Finite Differences Method：有限差分法）、ＦＶＭ（Finite Volumes Method：有限体積法）等も使用できる。タイヤの種類や境界条件等によって最も適当な解析手法を選択し、又は複数の解析手法を組み合わせて使用することが好ましい。
【００４２】
図３−１は、耐久性予測対象の回転体であるタイヤを示す側面図である。図３−２は、耐久性予測対象の回転体であるタイヤを微小要素に分割した例を示す一部斜視図である。図３−１、２に示すように、耐久性予測対象の回転体であるタイヤ１は、有限要素法に基づき、有限個の微小要素１０1、１０2・・・１０n等に分割されて、３次元有限要素モデル（３次元の回転体モデル、すなわち３次元のタイヤモデル１０）が作成される。有限要素法に基づく微小要素とは、例えば２次元平面においては三角形要素、四角形要素、３次元体としては四面体ソリッド要素、五面体ソリッド要素、六面体ソリッド要素等のソリッド要素や、三角形シェル要素、四角形シェル要素等のシェル要素等、コンピュータで用いうる要素とすることが望ましい。このようにして分割された微小要素は、解析の過程において３次元座標を用いて逐一特定される。
【００４３】
なお、本発明で回転体モデルあるいはタイヤモデル１０というときには、本発明に適用できる有限要素法、境界要素法その他の解析手法に基づき作成された解析用のモデルをいう。この実施例においては、有限要素法に基づいて微小要素に分割された有限要素モデルが、回転体モデルあるいはタイヤモデル１０となる。
【００４４】
前記幾何形状等の解析では、前記タイヤモデル１０の動的な状態、即ち、転動状態における動的解析によって、解析可能な数値データを作成している。前記タイヤモデル１０の転動状態を定義する際には、まず、前記タイヤモデル１０の回転の中心軸と当該タイヤモデル１０とが連動する設定を行う。その後、当該回転軸に回転速度を定義することにより、タイヤモデル１０全体の転動する速度を定義することができる。また、前記タイヤモデル１０の移動方向を定義し、当該移動方向に対する前記中心軸の角度を定義することにより、当該タイヤモデル１０の移動方向に対するタイヤモデル１０の角度を定義することができる。これにより、転動状態のタイヤモデル１０のスリップ角を定義することができる。
【００４５】
このように回転体解析部５１では、タイヤモデル１０の転動時の様々な状態を詳細に定義し、それらを動的解析することによって、忠実度の高い数値データとして解析している。この転動時の定義は、例えば、物理量が要素内を移動する速度を定義することで、回転体の転動状態を定義することができる。また、この解析では、上記のようにタイヤの転動を考慮し、準静的解析を含む動的解析によって、タイヤモデル１０の転動時の物理量である各微小要素の応力σｉやひずみεｉを求めている。このように、モデル１０を転動させることによって、より実際の現象に近い状態で応力σｉやひずみεｉを求めることができる。
【００４６】
次に、この数値データから、発熱エネルギ密度等の発熱エネルギ関連の特性を評価するために必要な物理量を抽出する（ステップＳ１０２、図１）。即ち、耐久性を評価する発熱エネルギ密度を計算するために必要な物理量を抽出する。この抽出する物理量は、具体的には、上記のように微小要素に分割された前記タイヤモデル１０の、転動時に刻々と変化する各微小要素の応力σｉと、ひずみεｉとを抽出する。
【００４７】
次に、抽出した物理量と、発熱を予測するのに必要変数を用いて、タイヤモデル１０の発熱エネルギ関連の特性を評価する（ステップＳ１０３、図１）。この発熱エネルギを評価するに際して、前記発熱エネルギ評価部５３が前記物理量より発熱エネルギを導出する手順について説明する。この発熱エネルギ導出手順は、リサージュ法を用いて、前記各微小要素の応力σｉと、ひずみεｉとにより、発熱エネルギ密度を求める。なお、ここでいう発熱は、タイヤを構成する粘弾性材料（例えばゴム）の変形に起因するタイヤの内部発熱である。図４は、本発明に係る回転体の耐久性予測方法における発熱エネルギ導出手順を示すフローチャートである。図５−１は、タイヤモデルを構成する微小要素の回転を示す概念図である。図５−２は、タイヤモデルの微小要素が１回転した場合における微小要素の応力変化を示す説明図である。図５−１に示すタイヤモデル１０の微小要素１０ｎが１回転（２π）すると、微小要素１０ｎの応力σｉは、回転角ωｔに対して、図５−２に示すように変化する。
【００４８】
タイヤの発熱エネルギ密度を求めるにあたっては、前記発熱エネルギ評価部５３は、タイヤモデル（３次元のタイヤモデル）１０を１回転させた場合における各微小要素の応力σｉ及びひずみεｉをフーリエ級数展開する（ステップＳ２０１）。例えば、図５−２に示した微小要素１０ｎの応力σｉでは、回転角ωｔが０〜２πまで、すなわちタイヤモデル１０が１回転した場合における応力σｉの変化をフーリエ級数展開する。
【００４９】
このように、タイヤモデル１０の１回転あたりにおける発熱エネルギを用いれば、当該発熱エネルギとタイヤモデル１０の回転数との積を求めることにより、所定回転数回転後におけるタイヤモデル１０の総発熱エネルギを簡易に求めることができる。本実施例ではタイヤモデル１０の１回転あたりにおける発熱エネルギを求めるが、タイヤモデル１０の単位回転数あたりにおける発熱エネルギを求めれば、同様の作用・効果を得ることができる。
【００５０】
次に、発熱エネルギ評価部５３は、応力σｉのフーリエ級数の各フーリエ次数について、材料の損失係数に応じた位相遅れδをひずみεｉに与え、各フーリエ次数における応力σｉとひずみεｉとのヒステリシスループの面積を演算する（ステップＳ２０２）。この内容について次に説明する。
【００５１】
タイヤを構成するゴムは粘弾性材料である。タイヤの転動時において粘弾性材料であるゴムの繰り返し変形によって発生するエネルギの一部が熱エネルギ（ヒステリシスロス）となり、このエネルギロスが転がり抵抗となる。図６−１は、粘弾性体を変形させた場合における応力とひずみとの関係を示す説明図である。図６−２は、応力とひずみとのヒステリシスループを示す説明図である。図６−１に示すように、粘弾性材料では応力σｉに対してひずみεｉの位相がδだけ遅れる。
【００５２】
これにより、図６−２に示すように、１周期（図６−１においては、回転角ωｔ1）の変形に対して応力σｉとひずみεｉとのヒステリシスループ５を描くことができる。このヒステリシスループ５の面積が、粘弾性材料の１周期における変形で発生した発熱エネルギであり、粘弾性材料の１周期における変形で損失したエネルギ、すなわち転動抵抗に相当する。したがって、各微小要素のヒステリシスループ５の面積から、タイヤの転動時における発熱エネルギ（転動抵抗）を求めることができる。なお、前記位相遅れδは、粘弾性材料の種類によって異なる損失係数に応じて決定される。なお、粘弾性材料は周波数依存性があるため、得られたひずみの振幅に近い条件下で求められた損失係数を用いるとさらに発熱エネルギの予測精度が向上するため好ましい。
【００５３】
発熱エネルギ評価部５３は、各フーリエ次数における応力σｉとひずみεｉとのヒステリシスループ５の面積を求めたら、フーリエ次数とヒステリシスループ面積との積の総和を演算する（ステップＳ２０３）。この一連の演算過程を応力、ひずみの全成分（３次元の場合は計６成分）について発熱エネルギ評価部５３が反復実行することにより（ステップＳ２０４；Ｎｏ）、応力、ひずみの各成分毎に前記総和を演算する（ステップＳ２０４；Ｙｅｓ）。さらに、発熱エネルギ評価部５３は、これらの過程を全微小要素について反復実行し（ステップＳ２０５；Ｎｏ、ステップＳ２０６）、タイヤモデル１０の全微小要素の発熱エネルギを導出して（ステップＳ２０５；Ｙｅｓ、ステップＳ２０７）、モデル化したタイヤ全体の発熱エネルギ導出手順が終了する。
【００５４】
このように発熱エネルギは、タイヤモデル１０を形成する粘弾性材料によっても異なるため、発熱を予測するのに必要な変数は、例えば、タイヤモデル１０を形成するゴムに力が作用した際の位相遅れの度合いなど、ゴムや補強コードの物性が挙げられる。発熱エネルギ評価部５３は、この発熱を予測するのに必要な変数と、上記の物理量とにより、タイヤモデル１０の発熱エネルギ関連の特性を評価する。なお、この発熱エネルギ関連の特性は、タイヤモデル１０全体でみる場合と、前記微小要素ごとに評価する場合とがあるが、必要に応じてどちらか一方を使い分ける、或いは双方を使用してもよい。
【００５５】
次に、この発熱エネルギ関連の特性と、タイヤモデル１０の転動の状態を表す変数等を用いてタイヤモデル１０の耐久性を予測する（ステップＳ１０４、図１）。この耐久性の予測は、発熱と耐久性との相関関係、及びタイヤモデル１０の転動の状態、例えば、荷重の加わり方や、安全率なども考慮して予測をする。タイヤモデル１０の転動の状態を表す変数は、例えば、所定速度で走行させたときにおいてタイヤが破壊するまでの距離がある。また、発熱エネルギ関連の特性には、例えば、タイヤの回転軸を通る子午面断面内におけるすべての微小要素中における、発熱エネルギ密度の最大値等を用いることができる。
【００５６】
以上、本発明によれば、粘弾性材料を含む回転体の動的な状態における耐久性に影響を与えるパラメータとして、原則として粘弾性材料の変形に起因する発熱を採用する。その結果、粘弾性材料を含む回転体の実際の使用状況に即した耐久性を予測することができる。なお、前記回転体解析部５１での解析時の、タイヤモデル１０の転動速度の定義は、タイヤモデル１０に接触している接触体を定義し、タイヤモデル１０に対する当該接触体の移動速度を定義することにより、当該タイヤモデル１０の転動速度を定義してもよい。これにより、実際のタイヤの転動状態と同様に、接触体がタイヤモデル１０に接触した状態での転動状態を再現できる。
【００５７】
また、タイヤモデル１０のスリップ角の定義は、前記接触体の移動方向と、前記タイヤモデル１０の前記接触体と接触する部分の転動方向との角度とを定義することによって、前記タイヤモデル１０のスリップ角を定義することができる。これにより、実際のタイヤの転動状態と同様に、接触体がタイヤモデル１０に接触した状態での転動状態を再現できる。
【００５８】
また、前記発熱エネルギ評価部５３での、タイヤモデル１０の発熱エネルギ関連の特性を評価は、当該タイヤモデル１０のひずみエネルギ密度の振幅から算出する方法を用いて評価してもよい。
【００５９】
（評価例）
本発明に係る回転体の耐久性予測方法及び評価装置によって、タイヤの耐久性を評価した。耐久性予測対象であるタイヤは、タイヤサイズ１８５／６５Ｒ１４のタイヤとし、このタイヤと同一サイズの６種類のタイヤを試作して、本発明に係る回転体の耐久性予測方法によりタイヤの耐久性を予測した。試験条件は、前記サイズのタイヤを車両に装着して１００ｋｍ／ｈで走行させた場合において、前記タイヤが破壊されるまでの走行距離を予測した。この予測結果と、本発明に係る回転体の耐久性予測方法で得られた結果との関係を、図７に示す。
【００６０】
図７は、走行距離に対する発熱エネルギ密度の最大値を示す説明図である。なお、図７の横軸に示す発熱エネルギ密度の最大値は、タイヤの回転軸を通る子午面断面内におけるすべての微小要素中の最大値である。そして、図７の右側に向かうに従って発熱エネルギ密度の最大値は大きくなっている。また、図７の縦軸は、上記の実走行によってタイヤが破壊するまでの走行距離を示しており、上方に向かうに従って破壊されるまでの走行距離が長くなっている。
【００６１】
図７から明らかなように、発熱エネルギ密度の最大値が小さいほど、タイヤが破壊するまでの走行距離は長くなり、発熱エネルギ密度の最大値が大きくなるに従って、タイヤが破壊されるまでの走行距離は短くなる。これは、本発明に係る回転体の耐久性予測方法で得られた発熱エネルギ密度の最大値が低いほど、タイヤの転動時の耐久性が高く、発熱エネルギ密度の最大値が高いほど、タイヤの転動時の耐久性が低いことを示している。実際の走行においては、タイヤの発熱エネルギが大きくなる程、タイヤの耐久性は低下することが分かっている。従って、この実走行による評価結果から、本発明に係る回転体の耐久性予測方法は、実際の回転体の挙動を十分に再現して、回転体の耐久性を高い精度で評価できるといえる。
【産業上の利用可能性】
【００６２】
以上のように、本発明に係る回転体の耐久性予測方法及び回転体の耐久性予測用コンピュータプログラム、並びに回転体の耐久性予測装置は、回転体の耐久性予測に有用であり、特に、タイヤのように、粘弾性材料を構成材料に含み、且つ、回転体の動的な状態における耐久性を予測する場合に適している。
【図面の簡単な説明】
【００６３】
【図１】本発明に係る回転体の耐久性予測方法の処理手順を示すフローチャートである。
【図２】本発明に係る回転体の耐久性予測装置を示す説明図である。
【図３−１】耐久性予測対象の回転体であるタイヤを示す側面図である。
【図３−２】耐久性予測対象の回転体であるタイヤを微小要素に分割した例を示す一部斜視図である。
【図４】本発明に係る回転体の耐久性予測方法における発熱エネルギ演算手順を示すフローチャートである。
【図５−１】タイヤモデルを構成する微小要素の回転を示す概念図である。
【図５−２】タイヤモデルの微小要素が１回転した場合における微小要素の応力変化を示す説明図である。
【図６−１】粘弾性体を変形させた場合における応力とひずみとの関係を示す説明図である。
【図６−２】応力とひずみとのヒステリシスループを示す説明図である。
【図７】走行距離に対する発熱エネルギ密度の最大値を示す説明図である。
【符号の説明】
【００６４】
１  タイヤ
５  ヒステリシスループ
１０  タイヤモデル
１０ｎ  微小要素
５０  耐久性予測装置
５０ｐ  処理部
５０ｍ  記憶部
５１  回転体解析部
５２  物理量抽出部
５３  発熱エネルギ評価部
５４  耐久性予測部
５６  入出力ポート
６０  入出力装置
６１  入力手段
６３  表示手段
【図１】
2005306113.tif 000003
【図２】
2005306113.tif 000004
【図３−１】
2005306113.tif 000005
【図３−２】
2005306113.tif 000006
【図４】
2005306113.tif 000007
【図５−１】
2005306113.tif 000008
【図５−２】
2005306113.tif 000009
【図６−１】
2005306113.tif 000010
【図６−２】
2005306113.tif 000011
【図７】
2005306113.tif 000012