(57)【要約】【課題】  タイヤの解析等においてビード部分の解析を高精度で実施する。【解決手段】  実際のタイヤからビード部分を切り出してねじったり引っ張ったり等のビード剛性を試験する（ステップ１２０，１２２）。ビード部分をその断面が１又は数個の有限要素からなる簡略化した有限要素モデルに定めかつ異方性材料で定義し、解析を行ったときのねじりや引っ張りの反力が実測の反力と一致するように、異方性材料のヤング率やポアソン比を調整する（ステップ１２４）。この結果の値（各方向のヤング率やポアソン比）を材料定数として、次のステップ１２６において入力する。このように、ヤング率やポアソン比で定義した異方性材料のビード部分は、簡略した有限要素モデルでも実際のビードと等価な材料物性を実現できる。【選択図】        図６
2006062554.tif 000002
【特許請求の範囲】
【請求項１】
使用状態でタイヤの挙動を模擬的に解析するために、タイヤを数値計算モデルに対応させて計算するタイヤモデルの一部を構成するタイヤのビード部分の解析モデルであって、
１または複数の金属線及び該金属線周辺部位のタイヤ材料から構成されるタイヤのビード部分について、異方性材料で定義したことを特徴とする
タイヤのビード部分の解析モデル。
【請求項２】
使用状態でタイヤの挙動を模擬的に解析するために、タイヤを数値計算モデルに対応させて計算するタイヤモデルの一部を構成するタイヤのビード部分の解析モデルであって、
１または複数の金属線及び該金属線周辺部位のタイヤ材料から構成されるタイヤのビード部分について、該ビード部分の断面領域を複数の要素領域に定め、少なくとも一つの要素領域に対応する部位を等方性材料で定義すると共に、それ以外の要素領域に対応する部位の少なくとも１つを異方性材料で定義したことを特徴とする
タイヤのビード部分の解析モデル。
【請求項３】
使用状態でタイヤの挙動を模擬的に解析するために、タイヤを数値計算モデルに対応させて計算するタイヤモデルであって、
請求項１または請求項２のタイヤのビード部分の解析モデルを含むことを特徴とする
タイヤモデル。
【請求項４】
車両の挙動を模擬的に解析するために、車両の各部の構造を数値計算モデルに対応させて計算する車体モデルであって、
請求項１または請求項２のタイヤのビード部分の解析モデル、或いは請求項３のタイヤモデル、を含むことを特徴とする
車体モデル。
【請求項５】
タイヤを数値計算モデルに対応させて計算するタイヤモデルの一部を構成するタイヤのビード部分の解析方法であって、
１または複数の金属線及び該金属線周辺部位のタイヤ材料から構成されるタイヤのビード部分を、異方性材料と定義し、
前記ビード部分を前記タイヤから取り出し、
前記取り出したビード部分について、複数の方向の引っ張り試験及びねじり試験の少なくとも一方の試験を予め行い、
試験結果の各方向の物理量に基づいて前記異方性材料の材料物性を定める
ことを特徴とするタイヤのビード部分の解析方法。
【請求項６】
タイヤを数値計算モデルに対応させて計算するタイヤモデルの一部を構成するタイヤのビード部分の解析方法であって、
１または複数の金属線及び該金属線周辺部位のタイヤ材料から構成されるタイヤのビード部分を、多数要素に分割して形成したビード部分の解析モデルを定め、
定めたビード部分の解析モデルを用いて、複数の方向の引っ張り解析及びねじり解析の少なくとも一方の解析を行い、
前記解析結果に基づいて前記ビード部分の引っ張り剛性及びねじり剛性の少なくとも一方の剛性を求め、
前記ビード部分を異方性材料モデルに置き換えると共に、該異方性材料モデルの引っ張り剛性及びねじり剛性の少なくとも一方の剛性が、前記求めた剛性となるように、該異方性材料モデルの材料物性を決定する
ことを特徴とするタイヤのビード部分の解析方法。
【請求項７】
前記異方性材料モデルの材料物性は、異方性部材のヤング率またはポアソン比であることを特徴とする請求項６に記載のタイヤのビード部分の解析方法。
【請求項８】
次の各ステップを含むタイヤの挙動シミュレーション方法。
（ａ）数値計算モデルとして接地及び転動により変形を与えることが可能なタイヤモデルとして、請求項１または請求項２に記載のタイヤのビード部分の解析モデルを含むタイヤモデルを定めるステップ。
（ｂ）タイヤ性能を使用状態で解析するために、前記タイヤモデルに使用条件を付与するステップ。
（ｃ）前記タイヤモデルの変形計算を実行するステップ。
（ｄ）前記ステップ（ｃ）におけるタイヤモデルに生じる物理量を求めるステップ。
（ｅ）前記物理量によりタイヤの性能を評価するステップ。
【請求項９】
コンピュータによってタイヤの挙動を解析するために、次の各ステップを含むことを特徴とするタイヤの挙動解析プログラム。
（Ａ）数値計算モデルとして接地及び転動により変形を与えることが可能なタイヤモデルとして、請求項１または請求項２に記載のタイヤのビード部分の解析モデルを含むタイヤモデルを定めるステップ。
（Ｂ）タイヤ性能を使用状態で解析するために、前記タイヤモデルに使用条件を付与するステップ。
（Ｃ）前記タイヤモデルの変形計算を実行するステップ。
（Ｄ）前記ステップ（Ｃ）におけるタイヤモデルに生じる物理量を求めるステップ。
（Ｅ）前記物理量によりタイヤの性能を評価するステップ。
【請求項１０】
コンピュータによってタイヤの挙動を解析するためのタイヤの挙動解析プログラムを記録した記録媒体であって、次の各ステップを含むことを特徴とするタイヤの挙動解析プログラムを記録した記録媒体。
（１）数値計算モデルとして接地及び転動により変形を与えることが可能なタイヤモデルとして、請求項１または請求項２に記載のタイヤのビード部分の解析モデルを含むタイヤモデルを定めるステップ。
（２）タイヤ性能を使用状態で解析するために、前記タイヤモデルに使用条件を付与するステップ。
（３）前記タイヤモデルの変形計算を実行するステップ。
（４）前記ステップ（３）におけるタイヤモデルに生じる物理量を求めるステップ。
（５）前記物理量によりタイヤの性能を評価するステップ。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、タイヤのビード部分の解析モデル、タイヤモデル、車体モデル、タイヤのビード部分の解析方法、タイヤの挙動シミュレーション方法、タイヤの挙動解析プログラム及びタイヤの挙動解析プログラムを記録した記録媒体にかかり、自動車等に使用される空気入りタイヤの性能を解析するときに用いられるタイヤのビード部分の解析モデル、タイヤモデル、車体モデル、そのモデルを用いた、、タイヤのビード部分の解析方法、タイヤの挙動シミュレーション方法、タイヤの挙動解析プログラム及びタイヤの挙動解析プログラムを記録した記録媒体に関する。
【背景技術】
【０００２】
タイヤ挙動についての解析は、実際に設計・製造したタイヤを計測したり自動車に装着して得た性能試験結果を用いたりしたものから、計算機（コンピュータ）環境の発達に伴って、計算機上でシミュレーションによって実現できるようになってきている。このタイヤ挙動をシミュレーションによって解析する主要な方法としては、有限要素法（ＦＥＭ）等の数値解析手法が主に用いられている。ＦＥＭは、構造体を有限個の要素でモデル化して、コンピュータを用いて構造体の挙動を解析する手法であり、その特徴から構造体を有限個の要素に分割する（以下、ＭＥＳＨ分割、または要素分割という。）ことが必要である。予測精度の高いタイヤ挙動をシミュレーションするためには、有限個の要素で構成されるシミュレーション用のタイヤモデル（数値データから構成されている）を如何に実際のタイヤ形状と同じように製作するかが重要である。
【０００３】
また、タイヤは、ホイルに組み付けられて用いられる。そこで、タイヤをＦＥＭによりシミュレーションするときにはホイルの部分を同時にモデル化する技術が知られている（例えば、特許文献１参照）。この場合は、タイヤが路面に押し付けられたときに、タイヤのビード部分がホイルにセットされた位置からずれようとする。通常のタイヤも同様な傾向があり、例えば６０ｋｍ／ｈで直進走行している車両が、急旋回しＪの字を描くようなターンをする場合、突然大きな力がタイヤに働き、タイヤがホイルからズレる現象が生じ、極端な場合にはズレ落ちて脱輪する場合がある。同様の現象は、その状態の規模に違いはあるが通常の走行条件でも起こり得る。
【特許文献１】特開２００２−３５０２９４号公報
【発明の開示】
【発明が解決しようとする課題】
【０００４】
しかしながら、タイヤをモデル化する場合、ビード部分を断面において１つの要素でモデル化することが一般的である。また、タイヤのモデル化手法に基づくタイヤーホイール組立体の有限要素モデル（ＦＥＭモデル）でも、ビード部分を断面において１つの要素でモデル化することが一般的である。タイヤのビード部分の断面において１つの要素でしかモデル化しない場合、実際のビードのような適切な動きを再現することは困難である。このため、コーナリング時に実際には、ビード部分が微小にホイルからずれる変形をするが、この変形量がシミュレーションと合致せず、高精度のシミュレーションができないときがある。
【０００５】
また、ビードは、数本の針金が複数束なって構成されていることが普通であり、この１本１本の針金を丁寧に有限要素モデル化することは複雑であるのに加えて、有限要素の要素数を増大させることになり、結果的に計算コスト高（計算時間がかかる）となってしまう。
【０００６】
すなわち、有限要素法を用いたタイヤーリムの解析においてはビード部分を簡略化して有限要素モデル化するのが一般的であり、その場合、ビードを構成する針金１本１本を正確にモデル化していないので、解析において実際のビードが変形している状態と異なった状態がシミュレーションされる場合がある。特に、簡略化したビード部分に、そのビードを通常構成する金属の物性値を設定した有限要素を使用すると、実際のビードよりも硬くモデル化される傾向にある。そのため、実際のタイヤではビード部分がリムに装着された位置からずれ動いているにも関わらず、シミュレーションの結果ではそれほど大きくずれ動かない結果を得る場合があった。
【０００７】
本発明は、上記事実を考慮して、有限要素法（ＦＥＭ）等の数値解析手法によるタイヤの解析等において、計算上でビード部分の解析を高精度で実施することを可能とする、タイヤのビード部分の解析モデル、タイヤモデル、車体モデル、タイヤのビード部分の解析方法、タイヤの挙動シミュレーション方法、タイヤの挙動解析プログラム及びタイヤの挙動解析プログラムを記録した記録媒体を得ることが目的である。
【課題を解決するための手段】
【０００８】
上記目的を達成するために本発明は、タイヤのビード部分について実際に則した挙動解析を可能としたものである。
【０００９】
詳細には、請求項１に記載の発明は、使用状態でタイヤの挙動を模擬的に解析するために、タイヤを数値計算モデルに対応させて計算するタイヤモデルの一部を構成するタイヤのビード部分の解析モデルであって、１または複数の金属線及び該金属線周辺部位のタイヤ材料から構成されるタイヤのビード部分について、異方性材料で定義したことを特徴とする。
【００１０】
本発明では、タイヤを数値計算モデルに対応させて計算するとき、１または複数の金属線及び該金属線周辺部位のタイヤ材料から構成されるタイヤのビード部分について、異方性材料で定義する。このように、ビード部分を異方性材料で定義することにより、実際のビードとほぼ等価なねじり剛性や引っ張り剛性を有するものとして扱うことができる。
【００１１】
請求項２に記載の発明は、使用状態でタイヤの挙動を模擬的に解析するために、タイヤを数値計算モデルに対応させて計算するタイヤモデルの一部を構成するタイヤのビード部分の解析モデルであって、１または複数の金属線及び該金属線周辺部位のタイヤ材料から構成されるタイヤのビード部分について、該ビード部分の断面領域を複数の要素領域に定め、少なくとも一つの要素領域に対応する部位を等方性材料で定義すると共に、それ以外の要素領域に対応する部位の少なくとも１つを異方性材料で定義したことを特徴とする。
【００１２】
ビード部分は、１または複数の金属線及び該金属線周辺部位のタイヤ材料から構成されることが多く、それらの部位によっては、等方性を有したり異方性を有したりする。そこで、ビード部分の断面領域を複数の要素領域に定めて、少なくとも一つの要素領域に対応する部位を等方性材料で定義すると共に、それ以外の要素領域に対応する部位の少なくとも１つを異方性材料で定義する。これによって、ビード部分に混在する等方性及び異方性を考慮した材料として扱うことができる。
【００１３】
請求項３に記載の発明は、使用状態でタイヤの挙動を模擬的に解析するために、タイヤを数値計算モデルに対応させて計算するタイヤモデルであって、請求項１または請求項２のタイヤのビード部分の解析モデルを含むことを特徴とする。
【００１４】
前記のようにモデル化したビード部分を含んでタイヤをモデル化することにより、実際のビードとほぼ等価なねじり剛性や引っ張り剛性を有するビード部分を含んだタイヤモデルを構築できる。
【００１５】
請求項４に記載の発明は、車両の挙動を模擬的に解析するために、車両の各部の構造を数値計算モデルに対応させて計算する車体モデルであって、請求項１または請求項２のタイヤのビード部分の解析モデル、或いは請求項３のタイヤモデル、を含むことを特徴とする。
【００１６】
車両は、その挙動を解析する場合、タイヤを装着した状態で解析することが一般的である。この場合、前記タイヤのビード部分を解析モデルとして採用したり、前記ビード部分を含むタイヤモデルを解析モデルとして採用したりすることによって、より柔軟な解析が可能となる。例えば、実際の現象では、車両が突然大きく旋回した場合には、タイヤのビード部分が、ホイルからずれる場合があるが、このときの現象を異方性材料でビード部分を定義することにより、従来手法よりも正確に解析することが可能である。
【００１７】
請求項５に記載の発明は、タイヤを数値計算モデルに対応させて計算するタイヤモデルの一部を構成するタイヤのビード部分の解析方法であって、１または複数の金属線及び該金属線周辺部位のタイヤ材料から構成されるタイヤのビード部分を、異方性材料と定義し、前記ビード部分を前記タイヤから取り出し、前記取り出したビード部分について、複数の方向の引っ張り試験及びねじり試験の少なくとも一方の試験を予め行い、試験結果の各方向の物理量に基づいて前記異方性材料の材料物性を定める。
【００１８】
タイヤのビード部分を解析する場合、そのビード部分の材料物性を特定することが重要である。そこで、実際のビード部分をタイヤから取り出して引っ張り試験やねじり試験を行う。その結果に等価となるように、異方性材料の材料物性を定めることで、タイヤのビード部分を実際のビードに合致する解析を行うことができる。
【００１９】
請求項６に記載の発明は、タイヤを数値計算モデルに対応させて計算するタイヤモデルの一部を構成するタイヤのビード部分の解析方法であって、１または複数の金属線及び該金属線周辺部位のタイヤ材料から構成されるタイヤのビード部分を、多数要素に分割して形成したビード部分の解析モデルを定め、定めたビード部分の解析モデルを用いて、複数の方向の引っ張り解析及びねじり解析の少なくとも一方の解析を行い、前記解析結果に基づいて前記ビード部分の引っ張り剛性及びねじり剛性の少なくとも一方の剛性を求め、前記ビード部分を異方性材料モデルに置き換えると共に、該異方性材料モデルの引っ張り剛性及びねじり剛性の少なくとも一方の剛性が、前記求めた剛性となるように、該異方性材料モデルの材料物性を決定する。
【００２０】
まず、金属などの１本１本を含む詳細なビード部分の一部または全部を有限要素モデルで忠実にモデル化する。次に、この詳細なビード部分の有限要素モデルを、ねじったり引っ張ったりする解析シミュレーションを実行する。そして、詳細なビード部分の有限要素モデルと同じ部分について、材質を異方性材料とする簡略化した有限要素モデルを作成し、この簡略化したビード部分の有限要素モデルを同じようにねじったり引っ張ったりしたときのねじり反力引っ張り反力が詳細なモデルを用いた時と同等になるように、その異方性材料モデルの材料物性を決定する。これにより、ビード部分のモデル化を簡略化しつつ実際の正確なビード部分の解析が可能となる。
【００２１】
請求項７に記載の発明は、前記異方性材料モデルの材料物性は、異方性部材のヤング率またはポアソン比であることを特徴とする。
【００２２】
このように、異方性材料モデルの材料物性を、異方性部材のヤング率やポアソン比で定めることによって、異方性材料モデルを容易かつ簡単に操作することができる。
【００２３】
前記モデルを用いることで、ビード部分の挙動を実際に近い状態で予測でき、タイヤの挙動をシミュレーションすることができる。詳細には、（ａ）数値計算モデルとして接地及び転動により変形を与えることが可能なタイヤモデルとして、請求項１または請求項２に記載のタイヤのビード部分の解析モデルを含むタイヤモデルを定めるステップ、（ｂ）タイヤ性能を使用状態で解析するために、前記タイヤモデルに使用条件を付与するステップ、（ｃ）前記タイヤモデルの変形計算を実行するステップ、（ｄ）前記ステップ（ｃ）におけるタイヤモデルに生じる物理量を求めるステップ、（ｅ）前記物理量によりタイヤの性能を評価するステップ、を含むタイヤの挙動シミュレーション方法によって、使用状態でタイヤの挙動を模擬的に解析することができる。
【００２４】
また、コンピュータによってタイヤの挙動を解析する場合、（Ａ）数値計算モデルとして接地及び転動により変形を与えることが可能なタイヤモデルとして、請求項１または請求項２に記載のタイヤのビード部分の解析モデルを含むタイヤモデルを定めるステップ、（Ｂ）タイヤ性能を使用状態で解析するために、前記タイヤモデルに使用条件を付与するステップ、（Ｃ）前記タイヤモデルの変形計算を実行するステップ、（Ｄ）前記ステップ（Ｃ）におけるタイヤモデルに生じる物理量を求めるステップ、（Ｅ）前記物理量によりタイヤの性能を評価するステップ、を含むタイヤの挙動解析プログラムによって、精度良くタイヤの挙動を解析することができる。
【００２５】
さらに、コンピュータによってタイヤの挙動を解析する場合、（１）数値計算モデルとして接地及び転動により変形を与えることが可能なタイヤモデルとして、請求項１または請求項２に記載のタイヤのビード部分の解析モデルを含むタイヤモデルを定めるステップ、（２）タイヤ性能を使用状態で解析するために、前記タイヤモデルに使用条件を付与するステップ、（３）前記タイヤモデルの変形計算を実行するステップ、（４）前記ステップ（３）におけるタイヤモデルに生じる物理量を求めるステップ、（５）前記物理量によりタイヤの性能を評価するステップ、を含むタイヤの挙動解析プログラムを記憶媒体に記憶するようにし実行させ、データ収集するようにすれば、精度良くタイヤの挙動解析が予測可能となる。
【発明の効果】
【００２６】
以上説明したように本発明によれば、ビードとリムの接触部分の挙動を考慮して精度良くタイヤの挙動解析を可能とするビード部分のモデルを提供でき、効率的なタイヤ開発を実現できる、という効果がある。
【発明を実施するための最良の形態】
【００２７】
以下、図面を参照して本発明の実施の形態を詳細に説明する。本実施の形態は、タイヤの挙動解析に本発明を適用したものである。
【００２８】
〔第１実施の形態〕
図１には本発明のタイヤの挙動シミュレーション方法を実施するためのパーソナルコンピュータの概略が示されている。このパーソナルコンピュータは、データ等を入力するためのキーボード１０、予め記憶された処理プログラムに従ってタイヤの性能を予測するコンピュータ本体１２、及びコンピュータ本体１２の演算結果等を表示するＣＲＴ１４から構成されている。
【００２９】
なお、コンピュータ本体１２には、記録媒体としてのフレキシブルディスク（ＦＤ）が挿抜可能なフレキシブルディスクユニット（ＦＤＵ）を備えている。なお、後述する処理ルーチン等は、ＦＤＵを用いてフレキシブルディスクＦＤに対して読み書き可能である。従って、後述する処理ルーチンは、予めＦＤに記録しておき、ＦＤＵを介してＦＤに記録された処理プログラムを実行してもよい。また、コンピュータ本体１２にハードディスク装置等の大容量記憶装置（図示省略）を接続し、ＦＤに記録された処理プログラムを大容量記憶装置（図示省略）へ格納（インストール）して実行するようにしてもよい。また、記録媒体としては、ＣＤ−ＲＯＭやＤＶＤ等の光ディスクや、ＭＤ，ＭＯ等の光磁気ディスクがあり、これらを用いるときには、上記ＦＤＵに代えてまたはさらに、対応する装置を用いればよい。また、パーソナルコンピュータの他に、ワークステーションやスーパーコンピュータをタイヤ解析に用いてもよいことは勿論である。
【００３０】
（挙動シミュレーション）
図２は、本実施の形態にかかるタイヤの挙動解析プログラムの処理ルーチンを示すものである。ステップ１００では、挙動解析の対象となるタイヤの設計案（タイヤ形状、構造、材料など）を定める。なお、ステップ１００における設定はタイヤ設計案に限定されるものではなく、現存するタイヤを解析する場合を含む。すなわち、現存するタイヤそのものを対象のモデルとして設定してもよい。次のステップ１０２では、タイヤ設計案を数値解析上のモデルに落とし込むためのタイヤのタイヤモデルを作成する。このタイヤモデルの作成は、用いる数値解析手法により若干異なる。本実施の形態では数値解析手法として有限要素法（ＦＥＭ）を用いるものとする。従って、上記ステップ１０２で作成するタイヤモデルは、有限要素法（ＦＥＭ）に対応した要素分割、例えば、メッシュ分割によって複数の要素に分割され、タイヤを数値的・解析的手法に基づいて作成されたコンピュータプログラムヘのインプットデータ形式に数値化したものをいう。この要素分割とはタイヤ、及び路面（後述）等の対象物を小さな幾つかの（有限の）小部分に分割することをいう。この小部分ごとに計算を行い全ての小部分について計算した後、全部の小部分を足し合わせることにより全体の応答を得ることができる。
【００３１】
上記ステップ１０２のタイヤモデルの作成では、図３に示すタイヤモデル作成ルーチンが実行される。まず、ステップ１１２では、タイヤ径方向断面のモデル（タイヤ断面モデル、すなわちタイヤ断面データ）を作成する。また、タイヤ断面内のゴム、補教材（ベルト、プライ等、鉄・有機繊維等でできた補強コードをシート状に束ねたもの）をそれぞれ有限要素法のモデル化手法に応じてモデル化する。次のステップ１１４では、２次元データであるタイヤ断面データ（タイヤ径方向断面のモデル）を周方向に一周分（３６０度）展開し、タイヤの３次元（３Ｄ）モデルを作成する。
【００３２】
図４はタイヤ断面モデルの一例を示すもので、複数のゴム部材毎に分割されたカーカス２２を有する空気入りタイヤ２０を示している。このカーカス２２はビード２６により折り返されている。このカーカス２２の内側はインナーライナー２４とされ、インナーライナー２４に延長上にはビードゴム３６が配置している。また、折り返されたカーカス２２により形成される略三角形状の領域はビードフィラー２８とされている。カーカス２２の上方には、ベルト３０が配置しており、このベルト３０の半径方向外側には溝が形成されたトレッドゴム３２が配置し、カーカス２２の軸方向外側にはサイドゴム３４が配置している。なお、タイヤ断面モデルをゴム部材毎に複数分割した例を挙げたが、設計目的によって３角形等の任意の形状に分割してもよい。
【００３３】
次に、図３のステップ１１６では、タイヤ各部のゴムの構成材料を設定する。このステップでは、タイヤの各部位に応じた剛性などの材料特性の構成材料を選択する。次のステップ１１８では、ビード部分をモデル化する。本実施の形態では、ビード部分を簡略化して有限要素モデル化する。これは、数本の針金が複数束なって構成されていることが多いビードについて、１本１本の針金毎に有限要素モデル化することが複雑でかつ有限要素の要素数を増大させることが計算コスト高（計算時間がかかる）となることを解消するためである。そして、そのビード部分の要素に異方性材料を設定している。
【００３４】
ここで、本実施の形態におけるビード部分の剛性について図５を参照して説明する。図５には、複数の金属線からなるビードを長方形材料（異方性材料）に置き換えたときの剛性の説明図である。図５の例では、断面積を変えずに、１０本の金属製のビード線の束を、断面が長方形の異方性材に置き換えた場合を示す。図５（Ａ）は１０本の金属製のビード線の束、（Ｂ）は置き換えた断面長方形の異方性材料、（Ｃ）はｘ軸周りの曲げ剛性、（Ｄ）はｙ軸周りの曲げ剛性、（Ｅ）はｚ軸周りの曲げ剛性、（Ｆ）はせん断剛性、について示している。
【００３５】
まず、１０本の金属線の束からなるビードを、同じ断面積を持つ等方性の金属としてＦＥＭでモデル化した場合、曲げ剛性Ｉｘ、Ｉｙ、ねじり剛性Ｉｚが実際より過大に評価されていた。これは、複数のビード線からなる金属の間にはゴムが介在し、ゴムは金属よりはるかに柔らかいため、同じ断面積を持つ１つの金属と置き換えた場合よりも、実際のビードは曲がりやすく、ねじりやすいからである。
【００３６】
例えば、各種実験を行った結果、引っ張り剛性Ｅｚは、同じ断面積を持つ等方性金属に置き換え可能だが、曲げ剛性Ｉｘは等方性金属に置き換えた場合よりも３０％も実際は柔らかく、曲げ剛性Ｉｙは３５％柔らかくなった。また、ねじり剛性Ｉｚは、等方性の金属に置き換えた場合よりも、３８％柔らかくなったことが分かった。また、ビードコアを断面についてせん断する時のせん断剛性τｘｙも、５５％柔らかいことがわかった。
【００３７】
このように、ビード部分を等方性金属に置き換えると、引っ張りや曲げ、ねじりなどのようにビードを変形させる方向によって、ビード硬さが違うので、等方性材料に置き換えることは正確でなかった。このため、ビードの線１本１本を全てＦＥＭモデル化して解析するべきであるが、これはＦＥＭ解析の要素数を膨大として、工業的には時に非効率な場合がある。そこで、ビード部分に、ビード部の曲げ、ねじり、引っ張り強度を考慮した、異方性材料を採用することにした。すなわち、本実施の形態では、金属の束を、同じ断面積を持つ異方性材料で置き換える。
【００３８】
なお、置き換えた後の断面積は、置き換える前の断面積と必ずしも同じでなくても良い。置き換えた部材が、同じねじり剛性、引っ張り剛性、せん断剛性を持てば、断面積の差は解析結果に悪影響を与えないからである。また、置き換える断面の形状は、長方形だけではなく、円形や、そのほかの多角形や楕円でも良い。
【００３９】
また、異方性材料の定義は、（歪）＝（異方性材料成分）・（応力）で構成された次式で表すことができる。
【００４０】
【数１】
2006062554.tif 000003
【００４１】
但し、Ｅはヤング率、Ｇはせん断弾性率、νはポアソン比であり、添え字は方向を表している。
【００４２】
次に、上述の剛性を考慮したビード部分の解析モデルの作成は（ステップ１１８）、タイヤのビード部分について異方性材料で材料物性を定義した有限要素モデルでモデル化することでなされる。この処理では図６に示す異方性材料の材料定数の設定プロセスにより得られた値を入力することによって達成される。このプロセスでは、まず、図６のステップ１２０で、実際のタイヤからビード部分（円環状）を切り出す。その切り出したビード部分単体について、次のステップ１２２ではビード剛性を試験する。例えばねじったり引っ張ったり曲げたりする動きを試験装置を用いて再現し、そのときのねじり反力あるいは引っ張り反力を計測する。図７には、ビード剛性の試験方法の概要を示した。図７（Ａ）はビード環、（Ｂ）はせん断剛性測定、（Ｃ）は引っ張り剛性測定、（Ｄ）は第１の曲げ剛性測定、（Ｅ）は第２の曲げ剛性測定、（Ｆ）はねじり剛性測定を示している。
【００４３】
次のステップ１２４では、上記ステップ１２２の試験結果からビード部分の異方性材料を決定する。具体的には、解析モデル（複数の針金とその周辺ゴムからなるビードをその断面が１又は数個の有限要素からなる簡略化した有限要素モデル）について、上記計測と同様の動きでねじったり引っ張ったり曲げたりするシミュレーションを実行する。ここで、解析モデル（簡略化されたビード部分の有限要素モデル）の材質は異方性材料で定義し（上記の式による）、ねじりや引っ張りを行ったときの反力が実測の反力と一致するように、異方性材料のヤング率やポアソン比を調整する。ヤング率やポアソン比の値を適切に定義することにより異方性材料として設定したビード部分は、簡略したビード円環形状の有限要素モデルでも実際のビードと等価なねじり剛性、引っ張り剛性を実現できる。この結果の値（各方向のヤング率やポアソン比）を材料定数として、次のステップ１２６において入力する。
【００４４】
上記のようにして簡略したビード円環形状の有限要素モデル（解析モデル）を含むタイヤモデルを作成した後には、図２のステップ１０４へ進み、路面の設定すなわち路面モデルの作成と共に路面状態の入力がなされる。このステップ１０４では、路面をモデル化し、そのモデル化した路面を実際の路面状態に設定するために入力するものである。路面のモデル化は、路面形状を要素分割してモデル化し、路面の摩擦係数μを選択設定することで路面状態を入力する。例えば、路面状態により乾燥（ＤＲＹ）、濡れ（ＷＥＴ）、氷上、雪上、非舗装等に対応する路面の摩擦係数μが存在するので、摩擦係数μについて適正な値を選択することで、実際の路面状態を再現させることができる。
【００４５】
なお、流体モデルを作成して、路面とタイヤモデルの間に設けても良い。流体モデルは、タイヤの一部（または全部）および接地面、タイヤが移動・変形する領域を含む流体領域を分割し、モデル化するものである。
【００４６】
このようにして、路面状態の入力がなされると、次のステップ１０６において、境界条件の設定がなされる。この境界条件とは、タイヤモデルに解析上すなわちタイヤの挙動をシミュレートする上で必要なものであり、タイヤモデルに付与する各種条件である。このステップ１０６の境界条件の設定では、まず、タイヤモデルには内圧を与えて、タイヤモデルに回転変位及び直進変位（変位は力、速度でも良い）の少なくとも一方と、予め定めた負荷荷重とを与える。なお、路面との摩擦を考慮する場合は、回転変位（または力、速度でもよい）もしくは直進変位（または力、速度でもよい）のどちらか一方のみでよい。
【００４７】
次に、ステップ１０６までに作成されたり設定されたりした数値モデルをもとに、解析としてのタイヤモデルの変形計算を行う。すなわち、上記ステップ１０６で境界条件の設定が終了すると、ステップ１０８へ進み、タイヤモデルの変形計算を行う。このステップ１０８では、タイヤモデルおよび与えた境界条件より、有限要素法に基づいてタイヤモデルの変形計算を行う。この変形計算は、タイヤ転動時の状態を得るために（過渡的な状態を得るために）、タイヤモデルの変形計算を繰り返し（例えば１msec以内の計算を繰り返して行い）、その度に境界条件を更新するようにしてもよい。また、変形計算は、タイヤ変形が定常状態となることを想定した予め定めた計算時間を採用することができる。次のステップ１１０では、上述の計算結果を出力する。この計算結果とは、タイヤ変形時の物理量を採用する。具体的には、サイドのたわみ量や接地形状、接地圧分布、タイヤ中心に作用する横力等である。
【００４８】
なお、計算結果の出力は、タイヤの接地部の形状や接地圧の分布、タイヤ中心に作用する力等の値または分布を可視化することを採用してもよい。これらは計算結果の値や変化量または変化率、力の向き（ベクトル）そして分布から導出することができ、それらをタイヤモデル周辺やパターン周辺とを共に線図等で表せば、把握しやすく提示可能な可視化をすることができる。
【００４９】
このように、本実施の形態では、タイヤのビード部分を、その断面が１又は数個の有限要素からなる簡略化した有限要素モデルでモデル化して、そのビードモデル（解析モデル）の材質を異方性材料で定義することによって、実際のビードに等価な材料特性となるようなねじり剛性、引っ張り剛性を実現できる。従って、ビード部分は簡略化した要素により実際形状のように針金の１本１本をモデル化する手間が省け、また、タイヤ全体の有限要素の数を減らすことができ、計算コストを低減すなわち計算時間の短縮することができる。
【００５０】
〔第２実施の形態〕
本実施の形態は、上記実施の形態と同様の構成のため、同一部分には同一符号を付して詳細な説明を省略する。上記では、タイヤからビード部分を実際に取り出して、実験によってビードの材料特性を特定したが、本実施の形態では実際の試験による材料特性の測定に代えて計算により求めるものである。
【００５１】
本実施の形態では、図６のプロセスを図８の処理ルーチンに代えて実行する。このプロセスでは、まず、図８のステップ１３０で、タイヤモデルのビード部分（円環状）について詳細な解析モデルを作成する。すなわち、針金などの金属線の１本１本を含む詳細なビード部分の一部または全部を有限要素モデルで忠実にモデル化する。図９には、ビード部分をモデル化することの説明図を示した。図９（Ａ）は６本の細い金属線がねじられながら太い金属線を取り囲んでいる実際のケーブルビードの断面図であり、（Ｂ）は要素分割したモデルである。図９（Ｂ）では、要素（ＭＥＳＨ）が２次元で表されているが、実際には図の厚み方向にも広がる３次元のものである。
【００５２】
次のステップ１３２では、ステップ１３０で作成した詳細なビード部分の有限要素モデルを、ねじったり引っ張ったり曲げたりする解析シミュレーション（ＦＥＭ）を実行する。なお、このモデル（詳細なビード部分の有限要素モデル）の材質は金属線の部分は金属でその間はゴムでモデル化する。このモデルに対して、図７で示したビード剛性試験と同様の変形を与えるシミュレーションを実施する。例えばねじったり引っ張ったり曲げたりする動きをシミュレーションで再現して、そのときのねじり反力あるいは引っ張り反力を計算する。
【００５３】
次のステップ１３４では、上記ステップ１３２のシミュレーション結果からビード部分の異方性材料を決定する。具体的には、解析モデル（複数の針金とその周辺ゴムからなるビードをその断面が１又は数個の有限要素からなる簡略化した有限要素モデル）について、上記ねじったり引っ張ったり曲げたりするシミュレーションを実行する。この場合、詳細なビード部分の有限要素モデルと同じ部分について、材質を異方性材料とする簡略化した有限要素モデルを作成し、この簡略化したビード部分の有限要素モデルを同じようにねじったりひっぱたりしたときのねじり反力引っ張り反力が詳細なモデルを用いた時と同等になるように、その異方性材料を定義する。すなわち、解析モデル（簡略化されたビード部分の有限要素モデル）の材質は異方性材料で定義し（前記の式による）、ねじりや引っ張りを行ったときの反力が実測の反力と一致するように、異方性材料のヤング率やポアソン比を調整する。ヤング率やポアソン比の値を適切に定義することにより異方性材料として設定したビード部分は、簡略したビード円環形状の有限要素モデルでも実際のビードと等価なねじり剛性、引っ張り剛性を実現できる。この結果の値（各方向のヤング率やポアソン比）を材料定数として、次のステップ１３６において入力する。
【００５４】
このように、本実施の形態では、実際のタイヤを用いることなくＦＥＭによる解析によってビード部分をモデル化できるので、解析処理の負荷を軽減することができる。
【００５５】
〔第３実施の形態〕
本実施の形態は、上記実施の形態と同様の構成のため、同一部分には同一符号を付して詳細な説明を省略する。また、本実施の形態では、解析対象としてタイヤとリムの組み合わせた組立体モデルを採用する。
【００５６】
タイヤを様々なリム（ホイール）と組み合わせて利用することを考えた場合、リムとタイヤを１体の構造物として解析を行なうことが重要である。本実施の形態では、図２の処理に代えて図１０の処理が実行される。図２と図１０の差異は、ステップ１０３の組み立てモデルの作成処理が追加されたことと、図２のステップ１０８をステップ１０７の組立体モデルの変形計算に代えたことである。ステップ１０３ではリムをモデル化すると共に、タイヤモデルにリムモデルを組み付けた組立体モデルを作成する。この組立体モデルについてステップ１０７で変形計算を行う。
【００５７】
このように、本実施の形態では、タイヤモデルとリムモデル（ホイール）からなる組立体モデルを用いてＦＥＭによる解析を行うので、タイヤ単体のみならずリム付きタイヤの解析を容易にすることができる。
【００５８】
〔第４実施の形態〕
本実施の形態は、上記実施の形態と同様の構成のため、同一部分には同一符号を付して詳細な説明を省略する。また、本実施の形態では、解析対象としてタイヤを取り付けた車両モデルを採用する。
【００５９】
タイヤは装着された車両において性能が発揮される。そこで、タイヤ単体の解析以外に、タイヤを車両モデルに取り付けて、旋回するような解析を実施する場合にも、本発明は有効である。例えば、実際の現象では、車両が突然大きく旋回した場合には、タイヤのビード部分が、ホイルからずれる場合があるが、このときの現象を異方性材料でビード部分を定義することにより、従来手法よりも正確に解析することが可能である。
【００６０】
本実施の形態では、図２の処理に代えて図１１の処理が実行される。図２と図１１の差異は、ステップ１０３と同一の組み立てモデルの作成処理（ステップ１０３Ａ）と車両モデルの作成処理（ステップ１３０Ｂ）が追加されたことと、図２のステップ１０８をステップ１０９の走行模擬計算に代えたことである。ステップ１０３Ａではリムをモデル化すると共に、タイヤモデルにリムモデルを組み付けた組立体モデルを作成し、次にその組立体モデルを取り付けるための車両モデルを作成しかる組立体モデルを取り付ける。この車両モデルについてステップ１０９で走行模擬計算を行う。
【００６１】
このように、本実施の形態では、タイヤを装着する車両を含めたモデルを用いたＦＥＭによる解析を行ので、車両の挙動に沿ったタイヤの解析を容易に行うことができる。
【実施例】
【００６２】
次に、本発明の実施例を詳細に説明する。
【００６３】
〔第１実施例〕
本実施例としてモデル化・試作したタイヤは、タイヤサイズは１８５／６０Ｒ１４であり、内圧１．８ｋｇｆ／ｃｍ2（＝１７６ｋＰａ）、荷重４００ｋｇ（＝３９．２ｋＮ）で、タイヤを路面に押しつけた後に横方向に３０ｍｍせん断する解析を実施した結果を示す（図１２（Ａ）参照）。このとき、ＦＥＭモデルは、リムとビード部分のタイヤ表面は接触要素を配置して接触させている（図１２（Ｂ）参照）。また、図１２（Ｃ）、（Ｄ）には、横力付与前後のビード部分の周辺を示した。
【００６４】
横力を付与したとき、リムの表面とビードとの距離Ｄを比較した結果を以下の表１に示す。このとき、解析試験に併せて、同様の試験を実際に行い、断面形状をＸ線を用いて撮影し、リムとビードとの距離を実験的に求めた結果と、等方性金属で定義したときの解析結果、異方性金属で定義したときの解析結果を示した。
【００６５】
【表１】
2006062554.tif 000004
【００６６】
表１から理解されるように、ビード部分を単に金属（例えば鉄）でモデル化した場合、ビードが硬く振る舞い、リムからの移動量（ズレ量）が少ない解析結果になっている。これに対して、ビードの曲げやねじりを考慮し、これを等価に表す異方性材料で定義した場合、実験と程同様に変形する結果を得ることができ、正確な解析が実行されたことが理解できる。
【００６７】
〔第２実施例〕
本実施例は第１実施例と同一のタイヤをモデル化・試作し、速度時速５０ｋｍでタイヤを転動させる解析を行い、スリップ角ＳＡを０度から１０度まで変化させる解析を行った。その結果を図１３に示す。
【００６８】
図から理解されるように、実測値ではスリップ角ＳＡが高くなりと、横力が小さくなる傾向にある。ビード部分を等方性の金属（等方性材料）でモデル化すると、ビード部分が硬く、変形が微小であり、特性は平坦になり、上記傾向が現れない。これに対して、ビードの曲げやねじりを考慮し、これを等価に表す異方性材料で定義した場合、実験と程同様に変形する結果を得ることができ、正確な解析が実行されたことが理解できる。
【図面の簡単な説明】
【００６９】
【図１】本発明の実施の形態にかかる、タイヤの挙動シミュレーション方法を実施するためのパーソナルコンピュータの概略図である。
【図２】本発明の第１実施の形態にかかる、タイヤの挙動解析プログラムの処理の流れを示すフローチャートである。
【図３】タイヤモデル作成処理の流れを示すフローチャートである。
【図４】タイヤ断面モデルを示す線図である。
【図５】ビード部分の剛性の説明図であり、（Ａ）はビード線の束、（Ｂ）は断面長方形の異方性材料、（Ｃ）はｘ軸周りの曲げ剛性、（Ｄ）はｙ軸周りの曲げ剛性、（Ｅ）はｚ軸周りの曲げ剛性、（Ｆ）はせん断剛性、について示している。
【図６】第１実施形態にかかる異方性材料の材料定数の設定プロセスの流れを示すフローチャートである。
【図７】ビード剛性試験概要の説明図であり、（Ａ）はビード環、（Ｂ）はせん断剛性測定、（Ｃ）は引っ張り剛性測定、（Ｄ）は第１の曲げ剛性測定、（Ｅ）は第２の曲げ剛性測定、（Ｆ）はねじり剛性測定を示している。
【図８】第２実施形態にかかる異方性材料の材料定数の設定プロセスの流れを示すフローチャートである。
【図９】ビード部分をモデル化の説明図で、（Ａ）は実際のケーブルビードの断面図、（Ｂ）は要素分割したモデル図である。
【図１０】本発明の第３実施の形態にかかる、組立体の挙動解析プログラムの処理の流れを示すフローチャートである。
【図１１】本発明の第４実施の形態にかかる、車両の挙動解析プログラムの処理の流れを示すフローチャートである。
【図１２】第１実施例の解析結果を示し、（Ａ）は解析条件、（Ｂ）はリムとビード部分の配置（Ｃ）は横力付与前、（Ｄ）は、横力付与後を示す線図である。
【図１３】第２実施例の解析結果を示す特性図である。
【符号の説明】
【００７０】
１０    キーボード
１２    コンピュータ本体
１４    ＣＲＴ
３０    タイヤモデル
ＦＤ    フレキシブルディスク（記録媒体）
【図１】
2006062554.tif 000005
【図２】
2006062554.tif 000006
【図３】
2006062554.tif 000007
【図４】
2006062554.tif 000008
【図５】
2006062554.tif 000009
【図６】
2006062554.tif 000010
【図７】
2006062554.tif 000011
【図８】
2006062554.tif 000012
【図９】
2006062554.tif 000013
【図１０】
2006062554.tif 000014
【図１１】
2006062554.tif 000015
【図１２】
2006062554.tif 000016
【図１３】
2006062554.tif 000017