(57)【要約】【課題】精度良く氷とゴム材との接触状態をシミュレーションする。【解決手段】ゴム材を数値解析が可能な複数個の要素を用いてモデル化したゴム材モデルを設定するステップＳ１と、氷を数値解析が可能な複数個の要素を用いてモデル化した氷モデルを設定するステップＳ２と、前記ゴム材モデルと氷モデルとを予め定めた境界条件で接触させ各々の変形計算を行ない必要な物理量を取得するシミュレーションステップＳ４〜Ｓ８とを含むとともに、氷モデルの変形計算は、各要素について、破壊の有無を判定する処理と、破壊と判定されたときに該要素を取り除いて破壊を表現する処理を含むことを特徴とするゴム材と氷との接触シミュレーション方法である。【選択図】　図２
2004025963.tif 000002
【特許請求の範囲】
【請求項１】
ゴム材を数値解析が可能な複数個の要素を用いてモデル化したゴム材モデルを設定するステップと、
氷を数値解析が可能な複数個の要素を用いてモデル化した氷モデルを設定するステップと、
前記ゴム材モデルと氷モデルとを予め定めた境界条件で接触させ各々の変形計算を行ない必要な物理量を取得するシミュレーションステップとを含むとともに、
前記氷モデルの変形計算は、該氷モデルの各要素について、破壊の有無を判定する処理と、破壊と判定されたときに該要素を取り除いて要素の破壊を表現する処理とを含むことを特徴とするゴム材と氷との接触シミュレーション方法。
【請求項２】
前記破壊の有無を判定する処理は、前記氷モデルの要素に作用する主応力、主ひずみ、せん断応力、ミーゼス応力又はひずみエネルギーの値を予め定めた値と比較することにより行われることを特徴とする請求項１記載のゴム材と氷との接触シミュレーション方法。
【請求項３】
前記氷モデルの変形計算は、破壊された要素を取り除いた後に新たな表面形状を認識するステップを含むことを特徴とする請求項１又は２記載のゴム材と氷との接触シミュレーション方法。
【請求項４】
前記ゴム材モデルは、ゴム部分をモデル化したゴムモデルと、このゴム部分よりも硬度が大きくかつ該ゴム部分から突出し前記氷モデルに接触する硬質材をモデル化した硬質材モデルとを含むとともに、
前記氷モデルの各々の要素が、前記硬質材モデルよりも小さいことを特徴とする請求項１乃至３のいずれかに記載のゴム材と氷との接触シミュレーション方法。
【請求項５】
前記硬質材モデルは、前記ゴムモデルの表面から略垂直に突出した軸状体をなすことを特徴とする請求項４記載のゴム材と氷との接触シミュレーション方法。
【請求項６】
前記軸状体が６ないし８角形状であることを特徴とする請求項５記載のゴム材と氷との接触シミュレーション方法。
【請求項７】
前記ゴム材モデルは、タイヤのトレッドゴムの転動時の挙動に合わせて移動しつつ前記氷モデルと接触する境界条件が与えられることを特徴とする請求項１乃至６のいずれかに記載のゴム材と氷との接触シミュレーション方法。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、例えば氷路面などを頻繁に走行するスタッドレスタイヤなどの開発に役立つゴム材と氷との接触シミュレーション方法に関する。
【０００２】
【従来の技術及び発明が解決しようとする課題】
スパイク夕イヤは、凍結路を走行する夕イヤとして使用されていたが、アスファルトより硬度が大きい金属ピンを使用しているため、走行時にアスファルトを損傷させてしまい粉塵を発生させる。このため、現在、殆どの地域で使用が禁止され、近年では、スパイク夕イヤに代えてスタッドレスタイヤが採用されつつある。スタッドレスタイヤは、前記金属ピンがないため、スパイクタイヤに比べると氷上性能に劣るが、近年では種々の改良が重ねられ、発売当初と比較して氷上性能が大きく向上しつつある。しかし、凍結路運転時の安全性をより一層高めるために、さらに高い氷上性能の実現が期待され、開発が続けられている。
【０００３】
例えばタイヤの氷上走行性能を種々改善するにあたっては、やはり現実の車両テストを数多く必要とする。とりわけ凍結路は人工的に作り出すのが困難であるため、限られた期間でしかテストできず、この種のタイヤの開発コストや開発期間を増大させる原因となっていた。
【０００４】
そこで、近年では有限要素法といった数値解析手法を用いたコンピューターシミュレーションにより、タイヤを試作しなくてもある程度の性能を予測・解析する方法が提案されている。しかしながら、従来の提案では、タイヤを舗装路面或いは水膜が存在する路面上を走行させるシミュレーションに止まる。水は、解析モデルでは一般に非圧縮性の完全流体として取り扱われる。一方、氷は、一定の応力内では実質的に形状が変化しないものの、ある一定の応力を超えると破壊し形状が変化する特性を有する。このような特性を有した氷とゴムとの接触シミュレーションには、上記従来の提案では対応することができない。
【０００５】
発明者らは、鋭意研究の結果、前述した氷を数値解析が可能な要素を用いてモデル化するとともに、その各要素に破壊応力を設定し、該破壊応力と、計算により得られた該要素に作用する応力とを比較して計算された応力の方が大きければ該要素を取り除くことにより破壊を表現でき、氷の物理現象をシミュレーション上に的確に取り込み得ること、また氷の表面形状を再構築することで、新たに生成された氷の形状が接触するゴムに与える影響といった相互作用を解析しうることを見出した。そして、これによって、例えばタイヤでの凍結路走行を行う際のゴムと氷の接触シミュレーションを精度良く行いうることを見出し本発明を完成させるに至った。
【０００６】
以上のように、本発明は、氷路面などを頻繁に走行するスタッドレスタイヤなどの開発に役立つゴム材と氷との接触シミュレーション方法を提供することを目的としている。
【０００７】
【課題を解決するための手段】
本発明のうち請求項１記載の発明は、ゴム材を数値解析が可能な複数個の要素を用いてモデル化したゴム材モデルを設定するステップと、氷を数値解析が可能な複数個の要素を用いてモデル化した氷モデルを設定するステップと、前記ゴム材モデルと氷モデルとを予め定めた境界条件で接触させ各々の変形計算を行ない必要な物理量を取得するシミュレーションステップとを含むとともに、前記氷モデルの変形計算は、該氷モデルの各要素について、破壊の有無を判定する処理と、破壊と判定されたときに該要素を取り除いて要素の破壊を表現する処理とを含むことを特徴とするゴム材と氷との接触シミュレーション方法である。
【０００８】
また請求項２の発明は、前記破壊の有無を判定する処理は、前記氷モデルの要素に作用する主応力、主ひずみ、せん断応力、ミーゼス応力又はひずみエネルギーの値を予め定めた値と比較することにより行われることを特徴とする請求項１記載のゴム材と氷との接触シミュレーション方法である。
【０００９】
また請求項３記載の発明は、前記氷モデルの変形計算は、破壊された要素を取り除いた後に新たな表面形状を認識するステップを含むことを特徴とする請求項１又は２記載のゴム材と氷との接触シミュレーション方法である。
【００１０】
また請求項４記載の発明は、前記ゴム材モデルは、ゴム部分をモデル化したゴムモデルと、このゴム部分よりも硬度が大きくかつ該ゴム部分から突出し前記氷モデルに接触する硬質材をモデル化した硬質材モデルとを含むとともに、前記氷モデルの各々の要素が、前記硬質材モデルよりも小さいことを特徴とする請求項１乃至３のいずれかに記載のゴム材と氷との接触シミュレーション方法である。
【００１１】
また請求項５記載の発明は、前記硬質材モデルは、前記ゴムモデルの表面から略垂直に突出した軸状体をなすことを特徴とする請求項４記載のゴム材と氷との接触シミュレーション方法である。
【００１２】
また請求項６記載の発明は、前記軸状体が６ないし８角形状であることを特徴とする請求項５記載のゴム材と氷との接触シミュレーション方法である。
【００１３】
また請求項７記載の発明は、前記ゴム材モデルは、タイヤのトレッドゴムの転動時の挙動に合わせて移動しつつ前記氷モデルと接触する境界条件が与えられることを特徴とする請求項１乃至６のいずれかに記載のゴム材と氷との接触シミュレーション方法である。
【００１４】
【発明の実施の形態】
以下本発明の実施の一形態を図面に基づき説明する。
図１には、本発明のシミュレーション方法を実施するためのシミュレーション装置１の一例が示されている。この装置１は、本体１ａと、入力手段としてのキーボード１ｂ、マウス１ｃと、出力手段としてのディスプレイ装置１ｄとから構成されている。本体１ａには、図示していないが、演算処理装置（ＣＰＵ）、ＲＯＭ、作業用メモリー、磁気ディスクなどの大容量記憶装置、ＣＤ−ＲＯＭやフレキシブルディスクのドライブ１ａ１、１ａ２などの記憶装置を適宜具えている。そして、前記大容量記憶装置は後述するシミュレーション方法を実行するための処理手順（プログラム）が記憶された記憶媒体を含んでいる。
【００１５】
図２には、本発明のシミュレーション方法の処理手順の一例が示されている。先ず本実施形態では、数値解析が可能な要素でゴム材をモデル化したゴム材モデルを設定する（ステップＳ１）。数値解析が可能とは、例えば有限要素法、有限体積法、差分法又は境界要素法といった少なくとも一つの数値解析法にて系全体の変位などを計算により求め得ることを意味し、本例ではこの数値解析法として有限要素法を採用する。
【００１６】
図３は、ゴム材モデル２の一例を３次元上に視覚化して表し、図４はその一部を取り出して拡大したものである。本実施形態では、ゴム材として、スタッドレスタイヤのトレッドに配されているゴムブロックの一部をモデル化したものを例示している。例えば図５に拡大して示すように、ゴムブロックＧは、ゴムからなるゴム部分Ｇ１と、このゴム部分Ｇ１に埋着された氷よりも硬くかつアスファルトよりも柔らかい繊維、例えばガラス繊維等などの短繊維からなる硬質材Ｇ２とから構成されている。前記硬質材Ｇ２は、ゴムブロックＧの路面と接地する接地面Ｐから略垂直に突出して形成されている。
【００１７】
本実施形態のゴム材モデル２は、このようなゴムブロックＧをモデル化している。即ち、ゴム材モデル２は、ゴム部分Ｇ１をモデル化したゴムモデル２Ａと、このゴムモデル２Ａよりも硬度が大きくかつ該ゴム部分Ｇ１から突出する硬質材Ｇ２をモデル化した硬質材モデル２Ｂとを含んでいる。
【００１８】
ゴムモデル２Ａ、硬質材モデル２Ｂは、いずれも有限個の小さな要素２ａ、２ｂ、２ｃ…に分割してモデル化されることにより、前記コンピュータ装置１にて取り扱い可能な数値データを構成する。具体的には、各要素２ａ、２ｂ、２ｃ…を特定する要素番号、節点座標値、要素形状などが記憶される。また例えばゴムモデル２Ａでは、モデル化の対象となったゴム部分Ｇ１の材料特性、例えば密度、ヤング率、減衰係数などに応じた材料特性が定義される。
【００１９】
同様に、硬質材モデル２Ａについては、モデル化の対象となった硬質材Ｇ２の材料特性を参考として、各要素にその特性が定義される。この材料特性の中には破断強度が含まれる。従って、この硬質材モデル２Ａに大きな曲げ応力などが作用すると折損が表現できる。このような材料特性は、各々の要素について変形計算を行う際に利用される。
【００２０】
特に限定はされないが、ゴムモデル２Ａについては、３次元要素、とりわけ４ないし６面体ソリッド要素が好ましく、本実施形態では６面体要素を主体的に使用している。また短繊維をなす硬質材モデル２Ｂについては、前記ゴムモデル２Ａの表面から略垂直に突出した軸状体でモデル化されており、図５に示したゴム内部に属する硬質材ないしその部分については省略している。なお実際の短繊維は、ほぼ円柱状をなすものであるが、これを要素分割するに際しては、前記軸状体として例えば全体が６ないし８角形状でモデル化することが好適である。端面が６角柱未満の軸状体では計算精度が低下する傾向があり、逆に８角柱よりも多くしても計算時間が増加し好ましくない。なおタイヤに使用される短繊維をモデル化した硬質材モデル２Ｂの外径は、概ね数十μｍ程度である。従って、ゴム材モデル２の全体の一辺の大きさも、せいぜい０．５〜１ｍｍ程度で設定するのが望ましい。
【００２１】
次に本実施形態では、氷をモデル化した氷モデル（路面形成物モデル）を設定する処理を行う（ステップＳ２）。氷モデルは、ゴム材モデル２と同様に数値解析が可能な要素を用いてモデル化される。
【００２２】
図６（Ａ）には氷モデル３の全体例を視覚化して示し、図６（Ｂ）にはその一部を取り出して拡大した図を示す。本実施形態の氷モデル３は、厚さＴ、巾Ｗ、長さＬの矩形状をなす氷板をモデル化したものを示す。該氷モデル３は、前記氷板を小さな６面体ソリッド要素３ａ、３ｂ、３ｃ…で分割することによりモデル化されたものを例示している。氷モデル３は、その一方の面が前記ゴム材モデル２を接触させる接触面Ｓａとなる。本実施形態では、初期状態において、この接触面Ｓａを水平かつ平坦な面として構成している。ただし、必要に応じて凹凸面としたり又勾配を設けたり適宜形状を変化させることもできる。
【００２３】
図７には、氷モデル３を構成する各要素の特性として、主応力と歪との関係を示す。前記各要素３、３ｂ、３ｃ…は、一定の応力値σａまではフックの法則に従って弾性変形するものとして取り扱われる。ヤング率は、例えば−１０℃において、９．０×１０１０（ｄｙｎ／ｃｍ２　）〜９．９×１０１０（ｄｙｎ／ｃｍ２　）程度、またポアソン比は０．３１〜０．３７程度に設定するのが望ましい。従って、前記応力値σａまでは実質的に弾性体としての特性を示す。一方、各要素にσａを超える応力が作用すると、その要素は破壊するように設定される。
【００２４】
図５に示したようなゴムブロックＧを有するスタッドレスタイヤで氷結路を走行した場合、ゴム表面から突出した短繊維状の硬質材Ｇ２は氷面に貫入して高いグリップ力を発揮する。これは放射光等を利用した種々の実験により確認されている。このような現象は、氷の弾性変形によるものではなく、硬質材Ｇ２が氷を破砕しかつ破砕により形成された穴に貫入しているためと考えられる。このため、氷との接触をコンピュータを用いて解析する場合、氷が破壊応力以上の力で破壊されることをシミュレーション上に的確に取り込まなければ精度の良い解析を行うことはできない。そこで、本実施形態では、上述のような応力値σａを氷モデル３の破壊応力として定義している。前記破壊応力σａは、特に限定されることなく種々定めることができるが、実測値などを参考として、例えば圧縮応力で７０ＭＰａ、せん断応力として０．１ＭＰａ程度に設定するのが好ましい。なお破壊した要素の取り扱い等については後述する。
【００２５】
本実施形態では、氷モデル３の要素の破壊の有無を判定するパラメータとして、主応力を用いているが、これに代えて例えば主ひずみ、せん断応力、ミーゼス応力又はひずみエネルギーといったパラメータを用いることも可能である。そして、採用する各パラメータ毎に要素の破壊を判定する閾値が予め設定され、この閾値よりも計算されたパラメータが大きければ要素の破壊と判定することができる。なお応力状態は、通常３つの直交する主応力（ベクトル）で表されるが、これを下記式により単軸引張応力σｖｍに置き換えたものがフォンミーゼス応力である。このフォンミーゼス応力σｖｍは、スカラー量となるため強度判定が比較的容易となる。
σｖｍ＝√［｛（σ１　−σ２　）２　＋（σ２　−σ３　）２　＋（σ３　−σ１　）２　｝／２］
ここで、σ１　、σ２　、σ３　は、主応力である。
【００２６】
また本実施形態の氷モデル３の各要素３ａ、３ｂ、ｃ…は、その全てが前記硬質材モデル２Ｂよりも小さいものを例示している。ここで、硬質材モデル２Ｂと氷モデル３の個々の要素とにおいて大きさを比較する場合、硬質材モデル２Ｂの大きさについてはその軸状体の最大外径とし、氷モデル３の各要素３ａ、３ｂ、３ｃ…の大きさはについては、図６（Ｂ）に示すように、接触面Ｓａと平行な面における一辺の長さｘ、ｙ又はｚの最も大きい値とし、これら対比することにより行うことができる。このように氷モデル３の要素３ａ、３ｂ、３ｃ…を硬質材モデル２Ｂよりも小さく設定することによって、氷モデル３に硬質材モデル２Ｂが突き刺さる際に形成される穴の輪郭形状、深さなどをより現実のものに近づけ、シミュレーション精度を高めるのに役立つ。
【００２７】
また本実施形態の氷モデル３は、図６（Ａ）に示したように、Ｘ軸に沿って小さな要素が並ぶ細分化領域Ｘａと、前記要素よりも大きな要素が並ぶ非細分化領域Ｘｂとを交互に設定するとともに、Ｙ軸側にも細分化領域Ｙａと、非細分化領域Ｙｂとを設定し細分化領域Ｘａ、Ｙａとが交差する交差部Ａに、最も小さな要素を区画している。このような交差部Ａは、例えばゴム材モデル２の硬質材モデル２Ｂとの接触部とすることにより、コンピュータによる計算時間の増大を抑制しつつ精度の向上に役立つ。なお以上のステップＳ１、Ｓ２は、実行順序を逆としても良く、また並列処理しても良い。
【００２８】
次に本実施形態では、境界条件等を設定する（ステップＳ３）。本実施形態のシミュレーションでは、図８に略示するように、ゴム材モデル２と氷モデル３とを接触させかつ接触面Ｓａに沿ってゴム材モデル２をスライドさせるシミュレーションを行う。このために設定される条件としては、氷モデル３の接触面Ｓａの初期状態の形状、ゴム材モデル２を氷モデル３に押し付ける垂直力Ｆａ、ゴム材モデル２と氷モデル３との間のせん断力Ｆｂ、ゴム材モデル２と氷モデル３との間の静ないし動摩擦係数μ、ゴム材モデル２、及び／又は氷モデル３の変形計算を行う際の時間増分などが挙げられる。前記摩擦係数としては、氷の状態に応じて種々設定することができ、例えば０．０１〜０．４程度が好適である。また変形開始からの計算を行う時間増分としては、特に限定はされないが、１０−９〜　１０−５ｓｅｃ程度とするのが望ましい。
【００２９】
次に本実施形態では、前記ゴム材モデル２と氷モデル３とを前記境界条件で接触させ各々の変形計算を行ない必要な物理量を取得するシミュレーションステップを行う（ステップＳ４ないしＳ８）。
【００３０】
図２において、ステップＳ４ないしＳ８から明らかなように、本実施形態では、ゴム材モデル２の変形計算と氷モデル３の変形計算とを個別に行う。そして、ゴム材モデル２の変形計算で得られた該ゴム材モデル２の形状や速度データを氷モデル３の変形計算時の境界条件として与えるとともに（ステップＳ８）、氷モデル３の変形計算で得られた接触面の形状をゴム材モデル２の変形計算時の境界条件として与える（ステップＳ７）ものを例示する。いわゆる「連成」である。以下、詳細に説明する。
【００３１】
ゴム材モデル２の変形計算の一例は図９に示される。
本例では先ず時間増分△ｔを経過した後のゴム材モデル２の変形計算を行う（ステップＳ４１）。ゴム材モデル２の変形計算においては、氷モデル３の接触面Ｓａは変位しない壁として取り扱われる。従ってゴム材モデル２の変位は該氷モデル３の接触面Ｓａにより規制される。
【００３２】
ゴム材モデル２の変形計算には、本例では有限要素法を採用し、例えば下記式で示される運動方程式が用いられる。この運動方程式からゴム材モデル２の時刻Δｔでの加速度、速度、変位などを求めることができる。
【数１】
2004025963.tif 000003
【００３３】
そして、所定時刻が経過したか否かを判断し（ステップＳ４２）、結果が真の場合（ステップＳ４２でＹ）、処理を終える。ステップＳ４２で偽の場合には、時間増分を加算して再度変形計算を行う（ステップＳ４１）。これにより、各時刻での変位が逐次計算できる。
【００３４】
図１０には、氷モデル３の変形計算の具体的な処理手順の一例を示す。
氷モデル３の変形計算においては、例えば氷モデル３の各要素に付された要素番号に対応した序数ｉを１に初期化し（ステップＳ５１）、要素番号ｉの要素に作用する応力を計算する（ステップＳ５２）。氷モデル３の変形計算には、ゴム材モデル２の形状、速度のデータが参照されて行われる。即ち、図１１（Ａ）に示すように、ゴム材モデル２（この例では硬質材モデル２Ｂ）の位置を参照し、、このゴム材モデル２を介して垂直力Ｆａが作用する状態で応力計算が行われる。
【００３５】
次に、要素番号ｉの要素について計算された応力が、予め設定された前記破壊応力σａよりも大きいか否かを判断する（ステップＳ５３）。計算された応力が前記破壊応力σａと同等ないしそれよりも小さい場合（ステップＳ５３でＮ）、前記序数ｉが最終要素の番号か否かを判断する（ステップＳ５６）。ステップＳ５６で最終要素を示す番号ではないと判断された場合、序数ｉに１を加算し（ステップＳ５７）、ステップＳ５２以降を繰り返す。また最終要素を示す番号である場合（ステップＳ５６でＹ）、即ち全ての要素について処理を終えた場合には、指定した時間が経過している場合（ステップＳ５８でＹ）には氷モデル３の変形計算の処理を終え、図２のステップＳ６以後を実行する。指定した時間が経過していない場合（ステップＳ５８でＮ）、ステップＳ５１以降を新たに時間増分を加算して繰り返す。
【００３６】
一方、計算された応力が前記破壊応力σａよりも大きいと判断された場合（ステップＳ５３でＹ）、当該要素番号ｉの要素を氷モデル３から取り除く（ステップＳ５４）。例えば図１１（Ａ）において、硬質材モデル２Ｂが接触している要素３ａ１〜３ａ６について、上記ステップＳ５３の結果がいずれもＹと判断された場合、図１１（Ｂ）に示すように、要素３ａ１〜３ａ６は氷モデル３から取り除かれる。これにより要素（氷モデル）の破壊が表現される。
【００３７】
次に、氷モデル３の接触面Ｓａの更新を行う（ステップＳ５５）。接触面Ｓａの更新は、取り除かれた要素番号ｉの要素が無いとした場合の氷モデル３の表面形状を新たに再認識し、このデータをシミュレーション装置１のメモリ等に記憶させることにより行われる。従って、図１１（Ｂ）の例では、このステップＳ５５により、氷モデル３の接触面Ｓａは、元の水平面を保持した水平面部Ｓａ１と、この水平面部Ｓａ１から１要素分凹んだ凹部Ｓａ２とを含むものとして更新されかつ記憶される。なお以後は、同様にステップＳ５６を繰り返す。
【００３８】
次に図２に示すように、ステップＳ７、Ｓ８では、それぞれ別々に独立させて計算されたゴム材モデル２と氷モデル３との変形計算結果から、お互いに必要なデータを受け渡しさせ両モデル２、３の連成を行う。例えば次ループでのゴム材モデル２の変形計算には、図１２（Ａ）に示すように、新たに更新された氷モデル３の接触面Ｓａの形状データが境界条件（壁）として与えられる。他方、次ループでの氷モデル３の変形計算には、図１２（Ｂ）に示すように、ゴム材モデル２の新たに計算された位置や速度などが条件として与えられる。なおこの連成は、同時刻におけるゴム材モデル２、氷モデル３の状態で行われる。図１３、図１４には、氷モデル３にゴム材モデル２を先ず接触させかつ水平方向に所定の速度でスライドさせたときの氷モデル３の破砕状況の一例を視覚化している。
【００３９】
上述のような処理を行うことにより、氷モデル３にはゴム材モデル２の位置の変化に伴う新たな応力を作用させることができ、他方、ゴム材モデル２については、氷モデル３の要素が破砕することによって形成される新たな接触面Ｓａが与えられる。そして、このような計算を微小時間毎に繰り返すことによって、ゴム材モデル２と氷モデル３との間の相互作用を考慮に入れつつ、ゴム材モデル２、氷モデル３の時々刻々と変化する変形状態を連成させて計算しうる。なおこれらの連成処理の計算手順は例えば一般に知られている有限要素法解析プログラム（例えばＬＳ−Ｄｙｎａ（ＬＳＴＣ社製）などを用いて自動計算することができる。
【００４０】
またステップＳ６では、計算終了となる予め指定した時間が経過したかを判断し、ステップＳ６でＹと判断された場合、計算結果を出力し（ステップＳ９）、処理を終える。ステップＳ６での計算を終える時間は、実行するシミュレーションに応じ安定した計算結果が得られるよう種々定めることができる。
【００４１】
計算結果の出力には種々の情報を含むことができる。例えば、ゴム材モデル２に水平力Ｆｂを与えた場合、そのときに氷モデル３へと伝えられる前後方向力を取り出すことにより、凍結路におけるスタッドレスタイヤの駆動性能（又は制動性能）を評価、改善するのに役立つ。そして、これらの出力結果から、必要なゴム材の材料特性の改良、硬質材の硬度、外径の改良などを行い、さらにはゴム材にサイピングや細溝などを追加してその形状、深さ、厚さなどを種々変化させることにより、好適なシミュレーション結果が得られたゴム材、ひいてはタイヤを実際に試作することができる。これにより、例えば冬用のタイヤ、とりわけ短繊維をゴム部に配合したスタッドレスタイヤの開発期間を大幅に短縮するとともに開発コストを低減できる。
【００４２】
またゴム材モデル２にスリップ角を与えて氷モデル３上を走行させた場合、ゴム材モデル２に生じる横力を出力することにより、氷上でのタイヤのコーナリング性能を評価、解析することができる。なお出力する情報は、これらの値に限定されず、必要に応じて種々のものを出力することができる。
【００４３】
またゴム材モデル２は、例えば図１５に示すように、タイヤ全体を小さな要素Ｔａ、Ｔｂ、Ｔｃ…でモデル化したタイヤモデルＴの一部として構成することもできる。この場合、タイヤの氷路走行時における制・駆動性能や、コーナリング性能を精度良く計算する事ができる点で望ましい。しかしながら、このような方法では、硬質材モデルの短繊維の外径が数十μｍ　と非常に小さいため、要素数が通常の有限要素モデルよりも大となり、また計算ステップもより細かくなるなど非常に多くの計算時間を必要とする。従って、計算時間を考慮すると、上述のようにタイヤのトレッドに含まれるゴムブロックの一部ないし全体を抜き出してたゴム材モデルでシミュレーションを行うことにより、より実用的な解析を行うことができる。
【００４４】
図１６には本発明の他の実施形態を示す。
この実施形態では、例えば図１６（Ａ）に示すように、タイヤモデルＴを用いて路面モデルＲで転動シミュレーションを行い、接地するトレッドゴムの任意の一要素ＳＴａについて、転動時の挙動を調べる。この挙動には、同図（Ｂ）に示すように、前記要素ＳＴａの転動時の軌跡、各時刻における垂直力ｆａ、水平力ｆｂなどを含む。そして、この挙動を前記ゴム材モデル２の境界条件として与えてシミュレーションを行うことが望ましい。これにより、少ない要素数でタイヤの転動により近づけてシミュレーションを行うことができ、計算時間の大幅な増大を招くことなく氷路でのタイヤの性能をシミュレーションすることができる。
【００４５】
図１７には、垂直圧力１００ｋＰａ、水平移動速度５５ｃｍ／ｓ、硬質材モデル２Ｂと氷モデル３との摩擦係数を０．０５、ゴムモデル２Ａと氷モデル３との摩擦係数を０．１としてシミュレーションを行い、氷モデル３からの抗力の大きさを時刻歴として出力したグラフを示す。なおタイプＡのゴム材モデルとして、硬質材モデル２Ｂの最大外径を３３μｍとしその配設密度は硬質材モデル間距離が３２０μｍとなるように設定した。またタイプＢのゴム材モデルとして、硬質材モデル２Ｂの最大外径を１１μｍとしその配設密度はタイプＡの４倍とした。
【００４６】
同じ速度でゴム材モデル２を氷モデル３上を移動させた場合、ゴム材モデル２と氷モデル３との間に働く抗力が大きいものは氷上グリップ力が高いと言える。図１７から、硬質材モデル２Ｂの配設密度は小さいが外径が大きいタイプＡの方が氷上グリップ力に優れることが判る。
【００４７】
以上本発明の一実施形態について説明したが、本発明はこのような形態に限定される訳ではない。例えばゴム材モデル２を固定された氷モデル３の上でスライドさせる際にゴム材モデル２を固定し該ゴム材モデル２と接触している氷モデル３を移動させることもできる。またゴム材モデル２は、硬質材モデル２Ｂを含まなくても良い。さらにゴム材モデル２は、解析対象をタイヤに限定するものではない。
【００４８】
【発明の効果】
上述したように、本発明のシミュレーション方法にあっては、ゴム材を実際に試作しなくとも、氷との接触時の性能を大凡知ることができる。従って氷と接触するゴム部を有する例えばスタッドレスタイヤなどの開発期間、開発コストを低減するのに役立つ。また本発明のシミュレーション方法では、氷モデルの変形計算において、その要素について計算された応力が予め設定された破壊応力よりも大きいときに該要素を取り除いて破壊を表現するステップを含むことにより、実際に生じる氷の破壊現象をシミュレーション上に的確に取り込むことができる。従って、より精度の高い、シミュレーションが可能となる。また請求項３記載の発明のように、氷モデルの変形計算において、破壊された要素を取り除いた後、新たな表面形状を認識するステップを含むときには、この形状の変化がゴム材に及ぼす相互作用をも考慮でき、さらに解析精度を高めうる。
【００４９】
また請求項４記載の発明のように、前記ゴム材モデルは、ゴムをモデル化したゴムモデルと、このゴムよりも硬度が大きくかつ該ゴム材から突出する硬質材をモデル化した硬質材モデルとを含むとともに、前記氷モデルの各々の要素を前記硬質材モデルよりも小さく設定したときには、硬質材が氷を貫く貫通孔などをより精度良く解析できる。
【００５０】
また請求項５記載の発明のように、前記硬質材モデルは、前記ゴムモデルの表面から略垂直に突出した軸状体をなすときには、スタッドレスタイヤのトレッド面と同様の構成をモデル化するのに特に有効となる。従って、スタッドレスタイヤの氷上走行性能の解析に役立つ。
【００５１】
また請求項７記載の発明のように、前記ゴム材モデルは、タイヤのトレッドゴムの転動時の挙動に合わせて移動しつつ前記氷モデルと接触する境界条件が与えられたときには、ゴム材モデルをタイヤの挙動と同様の条件で氷モデルに接触させることができるため、スタッドレスタイヤの氷上走行性能の解析に役立つ。またタイヤ全体をモデル化する必要がないため、モデルの作成時間や計算時間の短縮化にも役立つ。
【図面の簡単な説明】
【図１】本発明のシミュレーション方法を実施するためのコンピュータ装置の構成図である。
【図２】本発明のシミュレーション方法の処理手順の一例を示すフローチャートである。
【図３】本発明のゴム材モデルの斜視図である。
【図４】その要部拡大図である。
【図５】タイヤのトレッドのゴムブロックの断面図である。
【図６】（Ａ）は氷モデルの一例を示す斜視図、（Ｂ）はその一部拡大図である。
【図７】氷モデルの要素の応力と体積との関係を示すグラフである。
【図８】氷モデルと硬質材モデルとの接触を説明するための線図である。
【図９】ゴム材モデルの変形計算の一例を示すフローチャートである。
【図１０】氷モデルの変形計算の一例を示すフローチャートである。
【図１１】（Ａ）、（Ｂ）は氷モデルと硬質材モデルとの接触を説明するための線図である。
【図１２】（Ａ）、（Ｂ）は氷モデルと硬質材モデルとの接触を説明するための線図である。
【図１３】氷モデルのシミュレーションを視覚化して示す線図図である。
【図１４】氷モデルのシミュレーションを視覚化して示す線図図である。
【図１５】タイヤモデルの一例を示す斜視図である。
【図１６】（Ａ）は本発明の他の実施形態を示す他意やモデルの転動状態を示す線図、（Ｂ）はその一要素の軌跡を示す線図である。
【図１７】シミュレーションから得た抗力の時刻歴を示すグラフである。
【図１８】（Ａ）、（Ｂ）は実施例に用いた硬質材モデルの一例を示す斜視図である。
【符号の説明】
１　シミュレーション装置
２　ゴム材モデル
２Ａ　ゴムモデル
２Ｂ　硬質材モデル
３　氷モデル
【図１】
2004025963.tif 000004
【図２】
2004025963.tif 000005
【図３】
2004025963.tif 000006
【図４】
2004025963.tif 000007
【図５】
2004025963.tif 000008
【図６】
2004025963.tif 000009
【図７】
2004025963.tif 000010
【図８】
2004025963.tif 000011
【図９】
2004025963.tif 000012
【図１０】
2004025963.tif 000013
【図１１】
2004025963.tif 000014
【図１２】
2004025963.tif 000015
【図１３】
2004025963.tif 000016
【図１４】
2004025963.tif 000017
【図１５】
2004025963.tif 000018
【図１６】
2004025963.tif 000019
【図１７】
2004025963.tif 000020
【図１８】
2004025963.tif 000021