空気入りタイヤを有限個の要素でモデル化したタイヤモデルを設定するステップと、前記タイヤモデルと接触する接触面を有するリムをモデル化したリムモデルを設定するステップと、前記タイヤモデルを前記リムモデルに装着する条件を与えて前記リムモデルに前記タイヤモデルを嵌合させたタイヤリム組立体モデルを設定するステップとを含むタイヤのシミュレーション方法であって、前記リムモデルの前記接触面は、関数によって定義されることを特徴とするタイヤのシミュレーション方法。前記タイヤモデルは、タイヤ回転軸を含む断面における二次元形状で設定され、前記リムモデルの前記接触面は、リムの回転軸を含む断面における二次元形状で設定され、かつ、前記接触面は、直線を表す関数と、円弧を表す関数とを滑らかに繋いで定義される請求項1記載のタイヤのシミュレーション方法。前記二次元形状のタイヤモデルと、前記二次元形状のリムモデルとを嵌合したタイヤリム組立体モデルをタイヤ周方向に展開して三次元のタイヤリム組立体モデルを設定するステップをさらに含む請求項2記載のタイヤのシミュレーション方法。しかしながら、上記のシミュレーション方法では、いずれもリムモデルが、実際のリムを微小要素で分割して作られる。このような作業は、リムを多数の要素で分割するための節点の座標を決定し、かつ、全ての座標をコンピュータに記憶させる必要があり、モデルの作成に多くの時間を要するという問題があった。また、図11に示されるように、リムaは、通常、空気入りタイヤのビード部bが着座するリムシート面a1と、ほぼ垂直にのびビード部bの外側面を支えるリムフランジ面a2とが交わりるリムヒール部a3等が滑らかな円弧からなる曲面で構成されている。一方、有限要素法で使用される微小要素は、二次元の場合には直線、三次元の場合には平面で構成されるため、リムヒール部a3等の曲面を精度良く表現することができない。このため、従来のシミュレーション方法では、タイヤモデルをリムモデルに嵌合させるステップにおいて、予期しない振動が生じたり、計算精度が悪化する等の問題があった。本発明は、以上のような問題点に鑑み案出なされたもので、リムモデルのタイヤモデルと接触する接触面を、関数によって定義することを基本として、モデルの作成時間を短縮するとともに、計算精度をも向上しうるタイヤのシミュレーション方法を提供することを主たる目的としている。本発明のうち請求項1記載の発明は、空気入りタイヤを有限個の要素でモデル化したタイヤモデルを設定するステップと、前記タイヤモデルと接触する接触面を有するリムをモデル化したリムモデルを設定するステップと、前記タイヤモデルを前記リムモデルに装着する条件を与えて前記リムモデルに前記タイヤモデルを嵌合させたタイヤリム組立体モデルを設定するステップとを含むタイヤのシミュレーション方法であって、前記リムモデルの前記接触面は、関数によって定義されることを特徴とする。また請求項2記載の発明は、前記タイヤモデルは、タイヤ回転軸を含む断面における二次元形状で設定され、前記リムモデルの前記接触面は、リムの回転軸を含む断面における二次元形状で設定され、かつ、前記接触面は、直線を表す関数と、円弧を表す関数とを滑らかに繋いで定義される請求項1記載のタイヤのシミュレーション方法である。また請求項3記載の発明は、前記二次元形状のタイヤモデルと、前記二次元形状のリムモデルとを嵌合したタイヤリム組立体モデルをタイヤ周方向に展開して三次元のタイヤリム組立体モデルを設定するステップをさらに含む請求項2記載のタイヤのシミュレーション方法である。本発明は、空気入りタイヤを有限個の要素でモデル化したタイヤモデルを設定するステップと、前記タイヤモデルと接触する接触面を有するリムをモデル化したリムモデルを設定するステップと、前記タイヤモデルを前記リムモデルに装着する条件を与えて前記リムモデルに前記タイヤモデルを嵌合させたタイヤリム組立体モデルを設定するステップとを含むタイヤのシミュレーション方法において、リムモデルの前記接触面を、関数によって定義することを特徴とする。このように、本発明で用いられるリムモデルは、関数によって接触面を特定しうる結果、従来のリムモデルのように、多数の微小要素多数の節点を用いてリムを分割する必要が無い。従って、リムモデルの作成時間が低減される。また、本発明では、リムモデルの接触面が、関数を用いて定義されるため、平面はもとより、曲面についても二次関数等を用いて定義できる。従って、リムモデルの接触面を、滑らかに表現することもできる。これにより、タイヤモデルをリムモデルに嵌合させるステップにおいて、予期しない振動の発生や計算精度の悪化が抑制される。以下、本発明の実施の一形態が図面に基づき説明される。図1には、本実施形態のタイヤのシミュレーション方法の処理手順の一例が示されており、この処理は、コンピュータ図示省略を用いて実行される。本実施形態では、先ず、空気入りタイヤを有限個の要素でモデル化したタイヤモデルが設定されるステップS1。図2に示されるように、本実施形態のタイヤモデル2は、タイヤ回転軸を含む断面形状を有する二次元形状で設定される。該タイヤモデル2は、解析対象となる空気入りタイヤ図示省略を有限個の微小な要素2a、2b、2c...に分割して離散化することにより、コンピュータにて取り扱い可能な数値データとして設定される。具体的には、各要素2a、2b、2c...の節点座標値、要素形状、材料特性例えば密度、ヤング率、減衰係数などが定義され、コンピュータに記憶される。各要素2a、2b、2c...には、例えば2次元平面としては四辺形要素が好適である。また、3次元のタイヤモデルを作成する場合には、複雑形状を表現するのに適した4乃至6面体のソリッド要素などが用いられる。本実施形態のタイヤモデル2には、トレッド部の表面のパターン形状縦溝も忠実に再現されているが、これらの形状を省略することもできる。また、空気入りタイヤは、カーカスやベルト層等の繊維コード材で補強されるが、これらの繊維コード材は、例えば強度異方性が定義されたシェル要素2Sなどを用いてタイヤモデル2の中に組み入れられることが望ましい。さらに、空気入りタイヤには、リムに装着されるビード部に、非伸張性材料からなる環状のビードコアが埋設されているが、これについても、例えば剛要素を用いてタイヤモデル2の中に組み入れられるのが望ましい。次に、本実施形態では、リムをモデル化したリムモデルが設定されるステップS2。図3には、本実施形態のリムモデル3が視覚化して示されている。該リムモデル3は、リムの回転中心を含む断面形状を有する二次元形状で定義されており、図3には、右側のみが示されている。リムモデル3は、タイヤモデル2のビード部と接触する接触面4を有し、本実施形態のリムモデル3は、この接触面4のみで構成されている。本実施形態において、前記リムモデル3の接触面4は、基点Oからタイヤ軸方向外側にのびる直線からなるリムシート面4aと、このリムシート面4aのタイヤ軸方向の外端Aに滑らかに接続されかつ円弧からなるヒール円弧面4bと、このヒール円弧面4bのタイヤ軸方向外側の外端Bに滑らかに接続されかつ直線でのびるフランジ内方面4cと、このフランジ内方面4cの外端Cに滑らかに接続されかつ円弧からなるフランジ外円弧面4dとから構成されている。このフランジ外円弧面4dは、外端Dで終端している。前記接触面4において、直線からなるリムシート面4aは、例えば一次関数y=f1x{ただし、区間は点O~点A}で、また、直線からなるフランジ内方面4cは、例えば一次関数y=f3x{ただし、区間は点B~点C}でそれぞれ定義される。本実施形態のリムシート面4aは、タイヤ軸方向に対して5 ゚程度で傾斜している。同様に、フランジ内方面4cもタイヤ半径方向に対して5 ゚程度で傾斜している。また、接触面4において、円弧からなるヒール円弧面4bは二次関数y=f2x{ただし、区間は点A~B}で、また、円弧からなるフランジ外円弧面4dは、二次関数y=f4x{ただし、区間は点C~D}でそれぞれ定義される。そして、上記各関数は、区間を示す点O、A乃至Dの座標とともに、コンピュータに記憶される。また、本実施形態では、リムモデル3の前記接触面4は、実際のリムの変形が微小であることに鑑み、変化しない剛体表面として条件付けられる。ただし、このような条件に限定されるものではない。このように、本実施形態のリムモデル3の接触面4は、一次関数で定義される直線及び二次関数で定義される円弧だけを用いて、基点Oから外端Dまで連続してかつ凹凸なく滑らかに構成される。このような接触面4を有するリムモデル3は、本来のリムが有する輪郭形状に基づいて関数を定義するだけで足り、従来のように、多数の微小要素でリムを分割離散化する必要がない。従って、リムモデル3の作成時間が大幅に短縮される。図3には、タイヤモデル2の右側のビード部と接触する接触面4Rを有する右側のリムモデル3Rが示されているが、図4aに示されるように、タイヤモデル2の左側のビード部と接触する接触面4Lを有する左側のリムモデル3Lも、右側のリムモデル3Rと左右対称形状で形成することができる。なお、左右のリムモデル3L、3Rは、図4aに示したように、離間して互いに独立したもの以下、分離タイプということがある。でも良いし、図4bに示されるように、互いに連結されて一体化されたもの以下、一体タイプということがある。でも良い。次に、本実施形態では、タイヤモデル2をリムモデル3に装着する条件を与え、リムモデル3にタイヤモデル2を嵌合させたタイヤリム組立体モデルが設定されるステップS3。上記処理は、種々の方法で実行させることができる。例えば、リムモデル3が、図4aに示した分離タイプの場合、図5に示されるように、先ず、タイヤモデル2及びリムモデル3の互いの回転軸を揃えるとともに、左右のリムモデル3L、3Rのリム幅Rwaを、タイヤモデル2のビード部の幅Bwよりも大きく設定することにより、接触面4がタイヤモデル2のビード部と干渉しないように位置させる。次に、図6に示されるように、左右のリムモデル3L、3Rのリム幅を、正規のリム幅Rwになるようにタイヤ軸方向内側に移動させる。また、移動の後、タイヤモデル2の内腔面2iに空気圧に相当する等分布荷重が適用される。さらに、この際、両モデル2及び3の接触が判定され、接触していると判定された場合、一方のモデルは他方のモデルに対して壁として働き、互いの表面を通して力のやりとりが行われるとともに、タイヤモデル2の変形が計算される。そして、力の釣り合いが得られるまでこのような計算が繰り返される。この方法は、例えば上記特許文献2に記載されている。他方、リムモデル3が、図4bに示した一体タイプの場合、そのリム幅は、予め正規のリム幅Rwで定義される。このため、図7に示されるように、タイヤモデル2及びリムモデル3の互いの回転軸を揃えるとともに、例えばタイヤモデル2のビード部の幅Bwをリム幅Rwよりも小さくなるよう変形させ、互いに干渉しないように位置させる。次に、タイヤモデル2の内腔面2iに空気圧に相当する等分布荷重が適用することにより、タイヤモデル2を膨張変形させる。この際、上記の場合と同様、両モデル2及び3の接触が判定され、接触していると判定された場合、一方のモデルは他方のモデルに対して壁として働き、互いの表面を通して力のやりとりが行われるとともに、タイヤモデル2の変形が計算される。なお、この方法は、例えば上記特許文献1に記載されている。これらの処理を行うことにより、例えば、図6に示したように、リムモデル3にタイヤモデル2を所定の空気圧で嵌合させたタイヤリム組立体モデル5を設定することができる。なお、上記ステップS3において、タイヤモデル2と、リムモデル3の接触面4との接触判定がなされるが、本実施形態のような関数で定義された接触面4を有するリムモデル3を使用した場合、微小要素で分割されたリムモデルを使用する場合に比して、上記接触判定を高速に計算できる利点もある。即ち、従来では、タイヤモデル2の各節点が、リムモデル3の個々の要素毎に内側又は外側にあるかを計算して判断しなければならない。これに対して、本実施形態の場合、接触面4が関数で定義されているので、タイヤモデル2側の各節点が、この関数で表される面の内側にあるか又は外側にあるかだけを調べれば良い。また、関数の数種類は、要素の数よりも少ないため、本実施形態では、上記判定処理をより高速に行うことができる他、要素分割のように、分割方法にも依存しない利点がある。次に、本実施形態では、上記タイヤリム組立体モデル5を用いて、各種の変形シミュレーション例えば、転動シミュレーション、接地形状シミュレーション、耐久シミュレーション等が行われステップS4、その結果から各種の物理量が取得されるステップS5。なお、リムモデル3が分離タイプの場合、左右のリムモデル3L、3Rの並進自由度を一致させるのが望ましい。なお、上記実施形態では、二次元のタイヤモデル2及びリムモデル3を例に挙げて説明したが、本発明は、このような実施形態に限定されるものではない。即ち、上述のように、二次元形状のタイヤモデル2とリムモデル3とを嵌合させたタイヤリム組立体モデル5をタイヤ周方向に展開し、三次元のタイヤリム組立体モデルを設定しても良いのは言うまでもない。以上本発明の実施形態について詳細に説明したが、本発明は、上記実施形態に限定されるものではなく、本発明の要旨を逸脱しない範囲において種々変形することができるのは言うまでもない。サイズ195/65R15の解析対象のタイヤと、15×6.0JのリムJATMA規格準拠について、それぞれ二次元のタイヤモデル、リムモデルが下記の仕様にて設定され、モデルの作成時間や嵌合時の精度等が評価された。[実施例のタイヤリム組立体の仕様]タイヤモデル:図2に準拠し、要素数は1926個とした。リムモデル:図3に準拠した分離タイプとし、左右各々4種類の関数で定義された。[比較例のタイヤリム組立体の仕様]タイヤモデル:上記と同様リムモデル:図9に示されるように、分離タイプとし、左右各々17個計34要素の微小要素eで定義された。なお、図9中、丸印は節点を示す。先ず、リムモデルの作成時間については、実施例を100とした場合、比較例は150であり、50も時間が多くかかっていることが確認できた。また、嵌合ステップの計算に要した時間は、実施例を100とした場合、比較例は131であった。さらに、図10には、嵌合ステップを完了した状態が視覚化して示されており、図10bの比較例では、タイヤモデルのビード部がリムモデルからはみ出しており、嵌合状態が精度良く再現されていないのに対して、図10aの実施例では、精度の良い嵌合が得られていることが確認できた。