フィラーが配合されたゴム材料の変形をシミュレートするゴム材料のシミュレーション方法であって、前記ゴム材料を数値解析が可能な要素でモデル化したゴム材料モデルを設定するステップと、前記ゴム材料モデルに条件を設定して変形計算を行うステップと、前記変形計算から必要な物理量を取得するステップとを含むとともに、前記ゴム材料モデルは、厚さを有した三次元のモデルからなり、かつ、ゴムマトリックスをモデル化したマトリックスモデルと、フィラーをモデル化したフィラーモデルとを少なくとも含み、しかも前記ゴム材料モデルは、前記厚さ方向に同一断面形状が連続することを特徴とするゴム材料のシミュレーション方法。前記ゴム材料モデルを設定するステップは、x-y二次元座標上に、前記マトリックスモデル及び前記フィラーモデルを設定した二次元モデルを設定するステップと、前記二次元モデルの各節点を厚さ方向であるz方向に少なくとも1つ複写するステップと、前記複写された節点と前記二次元モデルの節点とを繋いで前記三次元のゴム材料モデルを設定するステップとを含む請求項1記載のゴム材料のシミュレーション方法。前記複写するステップは、前記二次元モデルの各節点をz方向に少なくとも2つ複写するステップを有することにより、前記三次元のゴム材料モデルは、前記厚さ方向の一方の面及び他方の面に表れる表面節点と、これらの表面節点に挟まれる内部節点とを有する請求項2記載のゴム材料のシミュレーション方法。前記表面節点に、厚さ方向に移動不能な条件を定義し、かつ、前記内部節点に、厚さ方向に移動可能な条件を定義するステップをさらに含む請求項3記載のゴム材料のシミュレーション方法。ところで、上述のようなゴム材料モデルa1又はa2を用いた変形シミュレーションにおいて、例えば、フィラーの充填率の増加やフィラーの分散状態の関係で、ゴム材料モデルの一部に局所的に大きな変形が生じると、計算が途中で中断されることがある。特に変形の自由度が小さい二次元のゴム材料モデルa1を使用した場合、変形の進行に伴い、要素e1の歪が大きくなって変形計算の続行が不可能になることがしばしばある。このように、二次元シミュレーションは、三次元シミュレーションに比べると、モデル作成時間は短縮できるが、変形に対して計算が脆弱である。したがって、ゴム材料モデルa1、a2において、フィラーの体積率が同一の場合、明らかに三次元モデルの方が大きな変形をシミュレートすることが可能である。この原因について、詳細には分かっていないが、二次元平面歪の場合、要素の節点の移動が平面内に限られることに関連していると考えられる。本発明は、以上のような問題点に鑑み案出なされたもので、ゴム材料モデルの作成時間の短縮化を図りつつ、大変形のシミュレーション計算が可能なゴム材料のシミュレーション方法を提供することを課題としている。本発明のうち請求項1記載の発明は、フィラーが配合されたゴム材料の変形をシミュレートするゴム材料のシミュレーション方法であって、前記ゴム材料を数値解析が可能な要素でモデル化したゴム材料モデルを設定するステップと、前記ゴム材料モデルに条件を設定して変形計算を行うステップと、前記変形計算から必要な物理量を取得するステップとを含むとともに、前記ゴム材料モデルは、厚さを有した三次元のモデルからなり、かつ、ゴムマトリックスをモデル化したマトリックスモデルと、フィラーをモデル化したフィラーモデルとを少なくとも含み、しかも前記ゴム材料モデルは、前記厚さ方向に同一断面形状が連続することを特徴とする。また請求項2記載の発明は、前記ゴム材料モデルを設定するステップは、x-y二次元座標上に、前記マトリックスモデル及び前記フィラーモデルを設定した二次元モデルを設定するステップと、前記二次元モデルの各節点を厚さ方向であるz方向に少なくとも1つ複写するステップと、前記複写された節点と前記二次元モデルの節点とを繋いで前記三次元のゴム材料モデルを設定するステップとを含む請求項1記載のゴム材料のシミュレーション方法である。また請求項3記載の発明は、前記複写するステップは、前記二次元モデルの各節点をz方向に少なくとも2つ複写するステップを有することにより、前記三次元のゴム材料モデルは、前記厚さ方向の一方の面及び他方の面に表れる表面節点と、これらの表面節点に挟まれる内部節点とを有する請求項2記載のゴム材料のシミュレーション方法である。また請求項4記載の発明は、前記表面節点に、厚さ方向に移動不能な条件を定義し、かつ、前記内部節点に、厚さ方向に移動可能な条件を定義するステップをさらに含む請求項3記載のゴム材料のシミュレーション方法である。本発明のシミュレーション方法で用いられるゴム材料モデルは、厚さを有した三次元のモデルからなり、かつ、ゴムマトリックスをモデル化したマトリックスモデルと、フィラーをモデル化したフィラーモデルとを少なくとも含んでいる。このような三次元のゴム材料モデルは、要素の節点の変形自由度が二次元のそれに比べて大きいため、大きな変形計算が可能になる。また、本発明で用いられるゴム材料モデルは、厚さ方向に同一断面形状が連続している。このようなゴム材料モデルは、二次元座標上に、前記マトリックスモデル及び前記フィラーモデルを設定した二次元モデルを設定し、これを厚さ方向に複写して容易に作成しうるため、三次元モデルでありながらもモデルの作成時間を短縮できる。従って、本発明のシミュレーション方法によれば、ゴム材料モデルの作成時間の短縮化を図りつつ、大変形のシミュレーション計算が可能になる。以下、本発明の実施の一形態が図面に基づき説明される。本発明は、フィラーが配合されたゴム材料の変形をシミュレートするゴム材料のシミュレーション方法である。前記フィラーとしては、例えば、例えば、カーボンブラック、シリカ、クレー、タルク、炭酸マグネシウム又は水酸化マグネシウム等が挙げられる。前記ゴム材料としては、例えば、天然ゴムNR、イソプレンゴムIR、ブチルゴムIIR、ブタジエンゴムBR、スチレンブタジエンゴムSBR、スチレンイソプレンブタジエンゴムSIBR、エチレンプロピレンジエンゴムEPDM、クロロプレンゴムCR又はアクリロニトリルブタジエンゴムNBRなどが挙げられる。図1には、本発明のシミュレーション方法を実施するためのコンピュータ装置1が示されている。このコンピュータ装置1は、本体1aと、入力手段としてのキーボード1b、マウス1cと、出力手段としてのディスプレイ装置1dとから構成されている。本体1aには、図示していないが、演算処理装置CPU、ROM、作業用メモリー、磁気ディスクなどの大容量記憶装置、CD-ROMやフレキシブルディスクのドライブ1a1、1a2が設けられている。そして、前記大容量記憶装置には後述する本発明のシミュレーション方法を実行するための処理手順プログラムが記憶されている。図2には、上記コンピュータ装置1を用いたシミュレーション方法の処理手順の一例が示される。本実施形態では、先ずゴム材料モデルが設定され、コンピュータ装置1に入力されるステップS1。図3には、このようなゴム材料モデルを設定する手順がさらに詳しく示される。本実施形態のゴム材料モデルの設定方法は、x-y二次元座標上に、マトリックスモデル及びフィラーモデルを設定した二次元モデルが設定されるステップS11。図4には、このような二次元モデルM2の一例が示されている。該二次元モデルM2は、解析しようとするゴム材料実在するか否かを問わないの微小領域が、有限個の小さな要素2a、2b、2c...を用いて離散化されたものである。各要素2a、2b、2c...には、例えば、2次元平面としての三角形ないし四辺形の要素が好ましく用いられる。各要素は、いずれも節点と、該節点間を繋ぐ辺とから構成される。また、この二次元モデルM2は、ゴムの任意の領域の微視構造を表しており、その一辺の長さは、例えば縦横それぞれ300nm×300nmとされている。また、二次元モデルM2には、ゴムマトリックス部分がモデル化されたマトリックスモデル3グレーの部分と、このマトリックスモデル3の中に分散して配されかつフィラー充填材がモデル化されたフィラーモデル4白い部分とを少なくとも含んでいる。前記マトリックスモデル3は、二次元モデルM2の主要部を構成している。また、フィラーモデル4は、マトリックスモデル3中に、予め定められた充填率に従い、分散して配置されている。本実施形態のフィラーモデル4は、球形の粒子が複数個連結したような形状で定義されているが、このような形状に限定されるものではない。また、必要に応じて、前記マトリックスモデル3と前記フィラーモデル4との間に介在しかつ両者の間の界面を形成する厚さが小さい界面モデル図示省略などを設けることもできる。次に、本実施形態では、前記二次元モデルM2の各節点を、厚さ方向であるz方向に少なくとも1つ複写するステップが行われるステップS12。このようなステップは、図5a乃至cに時系列的に示されている。先ず、図5aには、前記二次元モデルM2が簡略化して示される。該二次元モデルM2には複数個の節点e1が含まれている。次に、図5bには、前記二次元モデルM2の各節点e1が、厚さ方向であるz方向に少なくとも1つ、この例では、それぞれ距離を異ならせて2つ複写された状態が示される。これにより、二次元モデルM2からz方向にそれぞれ異なる距離を隔てて新たな節点e2及びe3が形成される。なお、節点e2同士、及び、節点e3同士は、それぞれ二次元モデルM2と同様に辺にて連結されているつまり、新たな2つの二次元モデルが、元の二次元モデルM2からz方向にオフセットされて形成される。。次に、図5cに示されるように、次に、前記複写された節点e2、e3と、これに対応する各々の前記二次元モデルM2の節点とを繋いで前記三次元のゴム材料モデルM3を設定するステップが行われるステップS13。即ち、z方向で隣り合う節点e1、e1及びe3をそれぞれ繋ぐことにより、多数の三次元の要素で構成された三次元のゴム材料モデルM3が形成される。このような三次元のゴム材料モデルM3は、前記厚さ方向z方向に同一断面形状が連続するものとして構成される。これにより、本実施形態の前記三次元のゴム材料モデルM3は、前記厚さ方向z方向の一方の面P1及び他方の面P2に表れる表面節点e1、e3と、これらの表面節点e1、e3に挟まれる内部節点e2とを有する。さらに、本実施形態では、前記表面節点e1、e3に、厚さ方向z方向に移動不能な条件を定義するただし、x及びy方向は移動可能とする一方、前記内部節点e2に、厚さ方向に移動可能即ち、三次元に移動可能な条件を定義する処理が含まれるステップS14、S15。従って、このような三次元のゴム材料モデルM3は、厚さ方向の要素数が増えることになるため、三次元に動くことが可能な節点が増えることにより、収束解が得られやすくなり、ひいては大きな変形まで計算することができる。なお、図4の二次元モデルM2を用いて上記手順によって得られた詳細な三次元のゴム材料モデルM3は、図6に示される。次に、新たに生成された三次元の各要素には、物性値等が入力されるステップS16。即ち、マトリックスモデル3を構成する要素についてはマトリックスゴムの物性値が、またフィラーモデル4を構成する要素についてはフィラーの物性値がそれぞれ入力されるステップS16。次に、変形シミュレーションを実行するために必要な各種の条件等が設定されるステップS2。設定される条件としては、上記三次元のゴム材料モデルM3を変形させるための変形条件が挙げられる。本実施形態では、図6のy方向に任意の平均ひずみ速度を加えてゴム材料モデルM3に引張変形を与える条件が定義される。ただし、変形条件は種々定めうるのは言うまでもない。次に本実施形態のシミュレーション方法では、上述のように設定された三次元のゴム材料モデルM3と、設定された条件に基づいて、変形シミュレーションが行われるステップS3。変形シミュレーション変形計算は、例えば有限要素法を用いたエンジニアリング系の解析アプリケーションソフトウエア例えば米国リバモアソフトウェアテクノロジー社で開発改良されたLS-DYNA等を用いて行うことができる。また、本シミュレーションは、均質化法漸近展開均質化法に基づいて行われる。均質化法は、図6に示した微視構造であるゴム材料モデルM3均質化法ではユニットセルとも呼ばれるを周期的に持っているゴム材料全体を表現する変数xiと、前記微視構造を表現するyiとの独立した2変数が用いられ、微視的スケールと巨視的スケールという異なる尺度の場におけるそれぞれ独立した変数を漸近展開することにより、図6に示した微視構造のモデル構造を反映させたゴム材料全体の平均的な力学応答を近似的に求めることができる。前記変形計算が行われると、その結果から必要な物理量を取得することができるステップS4。本実施形態では、物理量として、応力と歪の関係が出力される。また、前記ゴム材料モデルM3の各要素の時系列的な変形状態や物理量の分布を視覚化して表示することもできる。図7には、このような三次元のゴム材料モデルM3を用いて引っ張り試験を行った結果が示される。実施例には、図6に示した三次元のゴム材料モデルM3が用いられ、比較例には、図4に示した二次元のゴム材料モデルM2が用いられた。いずれも、均質化法を用い、同一の歪速度100mm/minでy軸方向の引張変形を与えた。なお、均質化法では、特開2010-205165号と同様、巨視領域20mm×20mmとし、最大変形量を2mm最大ひずみ10とした。図7から明らかなように、比較例では、約2の歪で計算がストップしたが、実施例では約10の歪までストップすることなく変形計算が可能であった。この理由としては、三次元のゴム材料モデルM3では、内部節点がz方向に移動しうる自由度を有するため。計算の安定性が増しているためと推察される。なお、実施例のモデルの作成時間については、比較例のモデルの作成時間と同様であった。