計算機上で、フィラー充填ポリマーの物性についてシミュレーションを行うため、ポリマー成分及びフィラーの平衡状態のモデルを作成する方法であって、前記ポリマー成分及び前記フィラーのモデルをそれぞれ用意するステップと、前記ポリマー成分モデルと前記フィラーモデルとを混合するステップと、前記フィラーモデルに前記ポリマー成分モデルが結合した状態を作成するステップとを含むことを特徴とするフィラー充填ポリマーの平衡モデルの作成方法。前記フィラー表面モデルと前記ポリマー成分モデルとの結合は、前記ポリマー成分モデルと前記フィラーモデルとの距離が設定した閾値以下となったときに行われることを特徴とする求項1に記載のフィラー充填ポリマーの平衡モデルの作成方法。前記平衡モデルは、散逸粒子動力学法DPD法を用いて作成されることを特徴とする請求項1に記載のフィラー充填ポリマーの平衡モデルの作成方法。前記ポリマー成分は、ゴム成分であることを特徴とする請求項1に記載のフィラー充填ポリマーの平衡モデルの作成方法。本発明の目的は、フィラーとポリマー成分との結合を考慮することによって、従来技術の問題を解決し、ポリマー成分及びフィラーの平衡状態を高精度に再現できるフィラー充填ポリマーの平衡モデルの作成方法を提供することにある。本発明者は、計算機上で、フィラー充填ポリマーの物性についてシミュレーションを行うため、ポリマー成分及びフィラーの平衡状態のモデルを作成する方法について、上記目的を達成するべく鋭意検討した結果、前記ポリマー成分のモデル及び前記フィラーのモデルをそれぞれ用意するステップ、並びに、前記ポリマー成分モデルと前記フィラーモデルとを混合するステップを含むことに加えて、前記フィラーモデルに前記ポリマー成分モデルが結合した状態を作成するステップを新たに含むことで、ポリマー成分及びフィラーの平衡状態を高精度に再現できることを見出した。本発明は、このような知見に基づきなされたもので、その要旨は以下の通りである。1計算機上で、フィラー充填ポリマーの物性についてシミュレーションを行うため、ポリマー成分及びフィラーの平衡状態のモデルを作成する方法であって、前記ポリマー成分及び前記フィラーのモデルをそれぞれ用意するステップと、前記ポリマー成分モデルと前記フィラーモデルとを混合するステップと、前記フィラーモデルに前記ポリマー成分モデルが結合した状態を作成するステップとを含むことを特徴とするフィラー充填ポリマーの平衡モデルの作成方法。2前記フィラー表面モデルと前記ポリマー成分モデルとの結合は、前記ポリマー成分モデルと前記フィラーモデルとの距離が設定した閾値以下となったときに行われることを特徴とする上記1に記載のフィラー充填ポリマーの平衡モデルの作成方法。3前記平衡状態の再現は、散逸粒子動力学法DPD法を用いることを特徴とする上記1に記載のフィラー充填ポリマーの平衡モデルの作成方法。4前記ポリマー成分は、ゴム成分であることを特徴とする上記1に記載のフィラー充填ポリマーの平衡モデルの作成方法。本発明によれば、従来の方法に比べて、ポリマー成分及びフィラーの平衡状態を高精度に再現できるフィラー充填ポリマーの平衡モデルの作成方法を提供することができる。以下に、図を参照しながら、本発明のフィラー充填ポリマーの平衡モデル作成方法を具体的に説明する。本発明によるフィラー充填ポリマーの平衡モデルの作成方法は、計算機上で、フィラー充填ポリマーの物性についてシミュレーションを行うにあたって、ポリマー成分及びフィラーの平衡状態のモデルを作成する方法である。そして本発明は、図1に示すように、前記ポリマー成分及び前記フィラーのモデルをそれぞれ用意するステップ図1と、前記ポリマー成分モデルと前記フィラーモデルとを混合するステップ図1と、前記フィラーモデルに前記ポリマー成分モデルが結合した状態を作成するステップ図1とを含むことを特徴とする。上記構成を具えることで、前記フィラーモデルに前記ポリマー成分モデルが結合した状態を加味した上で、フィラー充填ポリマーの平衡状態を再現することができる結果、従来のフィラー表面とポリマー成分との結合を考慮しないシミュレーション用モデルに比べて、より実物に近く優れた精度のシミュレーションが行える。さらに、本発明では、従来の方法に、1つの工程前記フィラーモデルに前記ポリマー成分モデルが結合した状態を作成するステップ図1を加えるだけで実現できることから、時間やコストが増大するという問題もない。なお、前記計算機とは、計算のための機械、器具のことをいい、通常、電子計算機コンピュータのことをいう。前記フィラー充填ポリマーとは、フィラーを含有するポリマーのことをいう。また、前記ポリマー成分及びフィラーの平衡状態とは、ポリマー成分及びフィラーが熱平衡状態最安定な状態であることをいい、前記フィラー充填ポリマーの平衡状態を再現することができれば、その物性について、高い精度でシミュレーションを行うことが可能となる。前記ポリマー成分としては、特に限定はされず、ゴム成分、樹脂成分、低分子液体成分等が挙げられ、これらを1種又は複数含むことができるが、従来技術に比べてポリマー成分及びフィラーの平衡状態を高精度に再現できる点から、前記ポリマー成分はゴム成分であることが好ましい。前記フィラーとしては、特に限定はされないが、カーボンブラック、シリカ、グラファイト化カーボン、ハイジライト等が挙げられる。前記フィラー充填ポリマーの物性については、計算機上でポリマー成分及びフィラーの平衡状態を再現することで、シミュレートできる物性であれば特に限定はされない。例えば、ヒステリシスロス、弾性率、摩耗特性等の物性を挙げることができる。<ステップ1>本発明では、まず、図1に示すように、前記ポリマー成分及び前記フィラーのモデルをそれぞれ用意するステップステップ1が行われる。前記ポリマー成分モデル及び前記フィラーのモデルの作成方法については、特に限定はされず、通常採用される方法に従って作成することができる。前記フィラーモデルについては、例えば図2~に示すような3次元のモデルを作成することができる。図2に示す単位フィラーを連結することで、図2及びに示すフィラーモデルを作成する。前記ポリマー成分モデルについては、例えば図3~に示すような3次元のモデルを作成することができる。図3に示す単位ポリマーを連結することで、図3及びに示すポリマー成分モデルを作成する。<ステップ2>本発明では、前記ステップ1の後、図1に示すように、前記ポリマー成分モデルと前記フィラーモデルとを混合するステップステップ2が行われる。前記ポリマー成分モデルと前記フィラーモデルとの混合状態の作成については、特に限定はされず、通常採用される方法に従って作成することができる。例えば、図4~に示すように、前記ポリマー成分モデル及び前記フィラーモデルを同じボックスに入れた後図4、ボックスサイズを小さくしながら図4、所定の運動方程式に従い各粒子の時間的発展を計算することで、前記ポリマー成分モデルと前記フィラーモデルとの混合状態を形成することができる図4。<ステップ3>本発明では、前記ステップ2の後、図1に示すように、前記フィラーモデルに前記ポリマー成分モデルが結合した状態を作成するステップステップ3が行われる。従来のフィラー充填ポリマーの解析モデルでは、図5に示すように、前記フィラーと前記ポリマー成分とが結合されていないが、本発明では、図5及びに示すように、前記フィラーモデルに前記ポリマー成分モデルを結合することで、解析モデルのポリマー成分及びフィラーの平衡状態が、より実物の状態に近づくことになる。前記フィラーモデルに前記ポリマー成分モデルを結合する方法については、特に限定はされないが、前記ポリマー成分モデルと前記フィラーモデルとの距離が定められた閾値以下となったときに結合するように設定することで、実物に近いポリマー成分及びフィラーの平衡状態を再現できる点で好ましい。一律に結合させた場合には、必要以上に前記ポリマー成分モデルと前記フィラーモデルとの結合が多くなり、所望の平衡状態を得られないおそれがある。前記閾値については、前記フィラー及び前記ポリマー成分の種類によって、適宜数値を変更することができる。また、前記フィラー及び前記ポリマーのそれぞれの結合点の位置についても指定することができる。前記ステップ3の後、ポリマー成分及びフィラーの平衡状態が再現された解析モデルがえられる図1。その後、平衡状態が再現された解析モデルを用いて、物性を予測するためのシミュレーションが行われる。また、前記ポリマー成分及びフィラーの平衡状態については、一定時間が経過した後、各パラメータ圧力、エネルギー及び構造などの変化がなくなることで判断されることが好ましい。ここで、前記平衡状態の再現は、計算に要する時間やコストを節約しながらフィラー充填ポリマーの初期状態を再現できる手法であるという点から、散逸粒子動力学法DPD法を用いることが好ましい。なお、シミュレーションの方法については特に限定はされず、従来のシミュレーション方法を採用すればよい。例えば、分子動力学法、有限要素法等が挙げられる。以下に、実施例を挙げて本発明を更に詳しく説明するが、本発明は下記の実施例に何ら限定されるものではない。実施例1実施例1では、計算機上で、前記フィラーモデルに前記ゴム成分モデルの片末端の主鎖が結合している状況で得られる平衡状態を散逸粒子動力学法DPD法によって再現した。具体的には、まず前記フィラー成分及び前記ポリマー成分のモデルをそれぞれ作成するステップ1を行った。前記フィラー成分のモデルとして、図2で示すように、DPD法による単位粒子以下、粒子と呼ぶを中心粒子から上下左右前後6方向に3粒子ずつ並べ、その隣接する粒子に下記式Aにより記述される調和振動子型のポテンシャルエネルギーUを発生させることで単位フィラーを作成した。作成した単位フィラーを図2に示すように、周期的に400個配置し、フィラー成分のモデルを作成した。なお、式A中のパラメータCは、C=2000とした。U/kBT=Cr2Aここで、rは隣接する粒子間距離である。また本実施例では、kBT=10とした。また、前記ポリマー成分のモデルとして、図3で示すように、粒子を直線上に40個並べ、隣接する粒子間に上記式Aで与えられるポテンシャルが発生するように単位ポリマーを作成し、作製した単位ポリマーを、図3に示すように空間的にランダムに3200個配置することで、モデルを作成した。この場合、式A中のパラメータCは、C=20とした。ステップ1の後、前記ゴム成分モデルと前記フィラーモデルとを混合するステップ2を行った。ステップ2では、ステップ1で作成したフィラーモデル及びポリマーモデルを、図4に示すように、ボックス領域内に配置した後、全単位粒子をDPD運動方程式下記で詳細説明に従って20000ステップ計算させ、それと同時にボックスのサイズを縮小を行った図4を参照。。ここで、前記ボックスのサイズ縮小スピードは、20000ステップ計算後の系におけるボックス内の単位粒子の単位体積当たりの数密度が、30になるように調整した。以上の処理を行うことで、最終的に図4に示されるるような前記ゴム成分と前記フィラーモデルの混合状態を得た。DPD運動方程式の詳細DPD運動方程式及びその積分方法は、非特許文献:RDGrootandPBWarrenJChemPhys1074423の、式1~式9を利用し、3次元周期境界条件の下で計算を実施した。必要になるパラメータとして、非特許文献の式3中の相互作用パラメータaijは、粒子種によらず常にaij=250、式5中のパラメータとして、kBT=10及びγ=45、式9中のパラメータとして、デルタt=001及びλ065、と選択した。ステップ2の後、前記フィラーモデルに前記ゴム成分モデルの末端が結合する状態を作成するステップ3を行った。ステップ3では、各ポリマーの構成粒子中の片末端において、構成する粒子とフィラー表面を構成する粒子フィラー表面粒子の間の距離が、粒子の直径の12倍以下に近づいたときに結合できるという命令、及び、結合した際にはそれら粒子間には式Aの調和振動子型ポテンシャルが発生するという命令を与えることで、全粒子をDPD方程式に従って10000ステップ計算する処理を行った。ゴム成分モデルの末端が結合する状態末端結合モデルのイメージを図5に示す。ステップ3の後、上記非特許文献の式3中のパラメータaijを、フィラー構成粒子とポリマー構成粒子の間のパラメータだけaij=400になるように変更し、全粒子をDPD運動方程式に従って100000ステップだけ計算することで、サンプルとなるゴム材料の平衡状態を再現した解析モデルを得た。実施例2実施例2では、前記フィラーモデルに前記ゴム成分モデルの主鎖が結合する状態を作成するステップ3において、各ポリマーの構成粒子のうち中心に位置する粒子に対して、構成する粒子とフィラー表面を構成する粒子フィラー表面粒子の間の距離が、粒子の直径の12倍以下に近づいたときに結合できるという命令、及び、結合した際にはそれら粒子間には式Aの調和振動子型ポテンシャルが発生するという命令を与えることで、全粒子をDPD方程式に従って10000ステップ計算する処理を行ったこと以外は、実施例1と同様の条件によって、サンプルとなるゴム材料の平衡状態を再現した解析モデルを得た。ゴム成分モデルの主鎖が結合する状態主鎖結合モデルのイメージを図5に示す。比較例比較例では、ステップ3として、全粒子をDPD運動方程式に従って10000ステップだけ計算する処理を行ったこと以外、つまり、前記フィラーモデルの表面と前記ゴム成分モデルとが結合される状態を作成しないこと以外は、実施例1と同様の条件によって、サンプルとなるゴム材料の平衡状態を再現した解析モデルを得た。評価各実施例及び比較例において得られたゴム材料の解析モデルについて、1解析モデルの断面の観察を行うことで評価を行った。具体的には、フリーソフトOCTA内のイメージングツールを用いて、解析モデルの断面画像を得た。得られた解析モデルの断面画像については、図6に示す図6が実施例1、図6が実施例2、図6が比較例である。。各実施例及び比較例において得られたゴム材料の解析モデルについて、2分散指数の算出を行った。分散指数については、得られたモデルに対して、異なる単位フィラーはお互いの表面距離が粒子の直径の105倍より近い時は接触していると仮定し、その下で接触している単位フィラーから構成されるフィラークラスターが何個の単位フィラーから構成しているかを評価し、その個数平均を計算することにより算出した。各解析モデルの分散指数について表1に示す。なお、表1の分散指数については、比較例の分散指数を100に規格化した場合の指標値として示し、数値が小さいほど、ゴム材料中にフィラーが均一に分散されており、良好な結果を示す。図6及び表1の結果から、実施例1及び2の分析モデルについては、比較例のモデルに比べて、ゴム材料中のフィラーが均一に分散されており、より実物のゴム材料に近い結果を示していることがわかった。