(57)【要約】【課題】  変形に伴う物理量を精度よく比較する。【解決手段】  コンピュータを用いて、フィラーとポリマーとの相互作用が互いに異なる少なくとも２種類の高分子材料の性能を比較するための方法である。このシミュレーション方法は、各高分子材料に基づいて、フィラーとポリマーとの間の相互作用ポテンシャルの強度が互いに異なる高分子材料モデルを、コンピュータにそれぞれ設定する工程Ｓ１と、コンピュータが、相互作用ポテンシャルに基づいて、各高分子材料モデルの構造緩和を計算する工程Ｓ２と、構造緩和後の各高分子材料モデルをそれぞれ変形させて、各高分子材料モデルの変形に伴う物理量を計算する工程Ｓ３とを含んでいる。構造緩和を計算する工程Ｓ３は、各高分子材料モデルの圧力、温度、並びに、フィラー及びポリマーを構成する粒子の数が互いに同一となる条件の下で、構造緩和を計算する。【選択図】図３
2019036253.tif 000002
【特許請求の範囲】
【請求項１】
コンピュータを用いて、フィラーとポリマーとの相互作用が互いに異なる少なくとも２種類の高分子材料の性能を比較するための方法であって、
前記各高分子材料に基づいて、前記フィラーと前記ポリマーとの間の相互作用ポテンシャルの強度が互いに異なる高分子材料モデルを、前記コンピュータにそれぞれ設定する工程と、
前記コンピュータが、前記相互作用ポテンシャルに基づいて、前記各高分子材料モデルの構造緩和を計算する工程と、
前記コンピュータが、構造緩和後の前記各高分子材料モデルをそれぞれ変形させて、前記各高分子材料モデルの前記変形に伴う物理量を計算する工程とを含み、
前記構造緩和を計算する工程は、前記各高分子材料モデルの圧力、温度、並びに、前記フィラー及び前記ポリマーを構成する粒子の数が互いに同一となる条件の下で、前記構造緩和を計算する、
高分子材料の粗視化分子動力学シミュレーション方法。
【請求項２】
前記物理量は、前記各高分子材料モデルの応力を含む請求項１記載の高分子材料の粗視化分子動力学シミュレーション方法。
【請求項３】
前記物理量は、前記各高分子材料モデルの空孔に関するパラメータを含む請求項１又は２記載の高分子材料の粗視化分子動力学シミュレーション方法。
【請求項４】
前記空孔に関するパラメータは、前記空孔の合計体積、又は、前記各高分子材料モデルの体積に対する前記空孔の前記合計体積の割合である請求項３記載の高分子材料の粗視化分子動力学シミュレーション方法。
【請求項５】
前記コンピュータが、前記各高分子材料モデルの前記物理量を比較する工程をさらに含む請求項１乃至４のいずれかに記載の高分子材料の粗視化分子動力学シミュレーション方法。
【請求項６】
前記相互作用ポテンシャルは、レナードジョーンズポテンシャルであり、
前記強度は、前記レナードジョーンズポテンシャルのポテンシャルの深さεである請求項１乃至５のいずれかに記載の高分子材料の粗視化分子動力学シミュレーション方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、高分子材料の粗視化分子動力学シミュレーション方法に関し、詳しくは、コンピュータを用いて、フィラーとポリマーとの相互作用が互いに異なる少なくとも２種類の高分子材料の性能を比較するための方法に関する。
【背景技術】
【０００２】
下記特許文献１は、コンピュータを用いて、フィラーとポリマーとを含む高分子材料のシミュレーション方法を提案している。下記特許文献１のシミュレーション方法では、先ず、高分子材料に基づいて、数値計算用の高分子材料モデルがコンピュータに設定される。フィラーとポリマーとの間には、斥力又は引力を作用させるための相互作用ポテンシャルが定義される。そして、下記特許文献１のシミュレーション方法では、相互作用ポテンシャルに基づいて、高分子材料モデルの構造緩和が計算される。
【０００３】
相互作用ポテンシャルには、その強度を定義するためのパラメータが含まれている。相互作用ポテンシャルの強度を小さくすると、フィラーとポリマーとの間の物理吸着が弱くなる。また、相互作用ポテンシャルの強度を大きくすると、フィラーとポリマーとの間の物理吸着が強くなる。
【先行技術文献】
【特許文献】
【０００４】
【非特許文献１】K.Hagita, H.Morita, H.Takano著、「Molecular dynamics simulation study of a fracture of filler-filled polymer nanocomposites」、Poymer、ELSEVIER、2016年9月2日、第99巻、p.368-375
【発明の概要】
【発明が解決しようとする課題】
【０００５】
上記非特許文献１では、相互作用ポテンシャルの強度が互いに異なる高分子材料モデルを設定して、高分子材料モデルの破壊に関するシミュレーションが行われている。しかしながら、上記非特許文献１のシミュレーションでは、相互作用ポテンシャルの強度が互いに異なる高分子材料モデルを対象とする構造緩和の計算方法が明示されていないため、各高分子材料モデルの変形に伴う物理量を、精度よく比較することが難しいという問題があった。
【０００６】
発明者らは、鋭意研究を重ねた結果、相互作用ポテンシャルの強度が互いに異なる高分子材料モデルの構造緩和をそれぞれ適切に計算しないと、フィラーに対するポリマーの分布がそれぞれ変化し、緩和計算後の各高分子材料モデルの内部応力が互いに大きく異なることを知見した。そして、内部応力が互いに大きく異なる各高分子材料モデルでは、変形に伴う物理量を精度良く比較できないことを見出した。
【０００７】
本発明は、以上のような実状に鑑み案出されたもので、変形に伴う物理量を精度よく比較することができる高分子材料の粗視化分子動力学シミュレーション方法を提供することを主たる目的としている。
【課題を解決するための手段】
【０００８】
本発明は、コンピュータを用いて、フィラーとポリマーとの相互作用が互いに異なる少なくとも２種類の高分子材料の性能を比較するための方法であって、前記各高分子材料に基づいて、前記フィラーと前記ポリマーとの間の相互作用ポテンシャルの強度が互いに異なる高分子材料モデルを、前記コンピュータにそれぞれ設定する工程と、前記コンピュータが、前記相互作用ポテンシャルに基づいて、前記各高分子材料モデルの構造緩和を計算する工程と、前記コンピュータが、構造緩和後の前記各高分子材料モデルをそれぞれ変形させて、前記各高分子材料モデルの前記変形に伴う物理量を計算する工程とを含み、前記構造緩和を計算する工程は、前記各高分子材料モデルの圧力、温度、並びに、前記フィラー及び前記ポリマーを構成する粒子の数が互いに同一となる条件の下で、前記構造緩和を計算することを特徴とする。
【０００９】
本発明に係る前記高分子材料の粗視化分子動力学シミュレーション方法において、前記物理量は、前記各高分子材料モデルの応力を含んでもよい。
【００１０】
本発明に係る前記高分子材料の粗視化分子動力学シミュレーション方法において、前記物理量は、前記各高分子材料モデルの空孔に関するパラメータを含んでもよい。
【００１１】
本発明に係る前記高分子材料の粗視化分子動力学シミュレーション方法において、前記空孔に関するパラメータは、前記空孔の合計体積、又は、前記各高分子材料モデルの体積に対する前記空孔の前記合計体積の割合であってもよい。
【００１２】
本発明に係る前記高分子材料の粗視化分子動力学シミュレーション方法において、前記コンピュータが、前記各高分子材料モデルの前記物理量を比較する工程をさらに含んでもよい。
【００１３】
本発明に係る前記高分子材料の粗視化分子動力学シミュレーション方法において、前記相互作用ポテンシャルは、レナードジョーンズポテンシャルであり、前記強度は、前記レナードジョーンズポテンシャルのポテンシャルの深さεであってもよい。
【発明の効果】
【００１４】
本発明の高分子材料の粗視化分子動力学シミュレーション方法は、前記各高分子材料に基づいて、前記フィラーと前記ポリマーとの間の相互作用ポテンシャルの強度が互いに異なる高分子材料モデルを、前記コンピュータにそれぞれ設定する工程を含んでいる。さらに、本発明の高分子材料の粗視化分子動力学シミュレーション方法は、前記コンピュータが、前記相互作用ポテンシャルに基づいて、前記各高分子材料モデルの構造緩和を計算する工程と、構造緩和後の前記各高分子材料モデルをそれぞれ変形させて、前記各高分子材料モデルの前記変形に伴う物理量を計算する工程とを含んでいる。
【００１５】
前記構造緩和を計算する工程は、前記各高分子材料モデルの圧力、温度、並びに、前記フィラー及び前記ポリマーを構成する粒子の数が互いに同一となる条件の下で、前記構造緩和を計算している。これにより、本発明の高分子材料の粗視化分子動力学シミュレーション方法は、構造緩和後の前記各高分子材料モデルの内部応力を互いに近似させることができるため、各高分子材料モデルの変形に伴う物理量を精度よく比較できる。
【図面の簡単な説明】
【００１６】
【図１】高分子材料の粗視化分子動力学シミュレーション方法を実行するコンピュータの一例を示す斜視図である。
【図２】ポリマーの一例を示す構造式である。
【図３】高分子材料の粗視化分子動力学シミュレーション方法の処理手順の一例を示すフローチャートである。
【図４】モデル設定工程の処理手順の一例を示すフローチャートである。
【図５】フィラーモデル及びポリマーモデルが配置されたセルの一例を示す概念図である。
【図６】図５のＡ部拡大図である。
【図７】フィラーモデルのフィラー粒子モデルの拡大図である。
【図８】ポリマーモデルの一例を示す概念図である。
【図９】フィラーモデル及びポリマーモデルの相互作用ポテンシャルの一例を説明する概念図である。
【図１０】緩和工程の処理手順の一例を示すフローチャートである。
【図１１】第２高分子材料モデルの体積分率と、時刻との関係を示すグラフである。
【図１２】第１高分子材料モデルの平衡内部応力及び第２高分子材料モデルの内部応力と、時刻との関係を示すグラフである。
【図１３】（ａ）は、変形計算前の高分子材料モデルの一部を示す概念図、（ｂ）は、変形計算後の高分子材料モデルの一部を示す概念図である。
【図１４】実施例の応力と歪みとの関係を示すグラフである。
【図１５】実施例の空孔の体積分率と歪みとの関係を示すグラフである。
【図１６】比較例の各高分子材料モデルの内部応力と時刻との関係を示すグラフである。
【図１７】比較例の空孔の体積分率と時刻との関係を示すグラフである。
【発明を実施するための形態】
【００１７】
以下、本発明の実施の一形態が図面に基づき説明される。
本実施形態の高分子材料の粗視化分子動力学シミュレーション方法（以下、単に「シミュレーション方法」ということがある。）では、コンピュータを用いて、フィラーとポリマーとの相互作用が互いに異なる少なくとも２種類（本実施形態では、２種類）の高分子材料の性能が比較される。
【００１８】
図１は、シミュレーション方法を実行するコンピュータの一例を示す斜視図である。コンピュータ１は、本体１ａ、キーボード１ｂ、マウス１ｃ及びディスプレイ装置１ｄを含んでいる。この本体１ａには、例えば、演算処理装置（ＣＰＵ）、ＲＯＭ、作業用メモリ、磁気ディスクなどの記憶装置、及び、ディスクドライブ装置１ａ１、１ａ２が設けられている。また、記憶装置には、本実施形態のシミュレーション方法を実行するためのソフトウェア等が予め記憶されている。
【００１９】
フィラーとしては、例えば、シリカ、カーボンブラック、又は、アルミナ等が採用される。ポリマーとしては、例えば、天然ゴム、合成ゴム、又は、樹脂等が採用される。
【００２０】
図２は、ポリマーの一例を示す構造式である。本実施形態のポリマーとしては、図２に示されるように、cis-１，４ポリイソプレン（以下、単に「ポリイソプレン」ということがある。）が例示される。ポリイソプレンを構成するポリマーは、メチン基等（例えば、−ＣＨ＝、＞Ｃ＝）、メチレン基（−ＣＨ２−）、及び、メチル基（−ＣＨ３）によって構成されるイソプレンのモノマー（イソプレン分子）３が、重合度ｎで連結されて構成されている。なお、ポリマーには、ポリイソプレン以外のものが用いられてもよい。
【００２１】
本実施形態のシミュレーション方法において、性能が比較される高分子材料モデルには、第１高分子材料及び第２高分子材料が含まれている。これらの第１高分子材料及び第２高分子材料は、フィラーとポリマーとの相互作用が互いに異なっている。本実施形態の第１高分子材料は、第２高分子材料よりも相互作用が弱いものとして区分されるが、第２高分子材料よりも相互作用が強いものであってもよい。また、本実施形態の相互作用としては、物理吸着である場合が例示される。
【００２２】
フィラーとポリマーとの間の物理吸着は、主にフィラーとポリマーとの間に働くファンデルワールス相互作用及び静電相互作用によって、ポリマーがフィラーの表面に吸着される現象である。フィラーとポリマーとの間の物理吸着が強いほど、高分子材料の耐破壊性が向上すると期待される。このような物理吸着の強弱は、実際の高分子材料において、例えば、分子の化学構造や分子密度等を調節することで適宜設定することができる。一方、後述のシミュレーションに用いられる高分子材料モデルでは、例えば、フィラーとポリマーとの間の相互作用ポテンシャルの強度（ポテンシャルの深さ）に設定される値を調節することで、物理吸着の強弱を定義することができる。
【００２３】
図３は、シミュレーション方法の処理手順の一例を示すフローチャートである。本実施形態のシミュレーション方法は、先ず、高分子材料に基づいて、相互作用ポテンシャルの強度が互いに異なる高分子材料モデルが、コンピュータ１に設定される（モデル設定工程Ｓ１）。本実施形態のモデル設定工程Ｓ１では、第１高分子材料をモデル化した第１高分子材料モデルと、第２高分子材料をモデル化した第２高分子材料モデルとがそれぞれ設定される。図４は、モデル設定工程Ｓ１の処理手順の一例を示すフローチャートである。
【００２４】
本実施形態のモデル設定工程Ｓ１では、先ず、高分子材料の一部に対応する仮想空間であるセルが、コンピュータ１に定義される（工程Ｓ１１）。セル４は、各高分子材料モデル１０（本例では、第１高分子材料モデル１０Ａ及び第２高分子材料モデル（図示省略））毎に定義される。図５は、フィラーモデル６、及び、ポリマーモデル７が配置されたセル４の一例を示す概念図である。
【００２５】
セル４は、少なくとも互いに向き合う一対の面５、５、本実施形態では、互いに向き合う三対の面５、５を有しており、直方体又は立方体（本例では、立方体）として定義されている。各面５、５には、周期境界条件が定義されている。このようなセル４が用いられることにより、後述の粗視化分子動力学計算において、例えば、ポリマー（図２に示す）をモデル化した後述のポリマーモデル７について、一方側の面５ａから出て行ったポリマーモデル７の一部が、他方側の面５ｂから入ってくるように計算することができる。従って、セル４は、一方側の面５ａと、他方側の面５ｂとが連続している（繋がっている）ものとして取り扱うことができる。
【００２６】
セル４の一辺の各長さＬ１は、適宜設定することができる。本実施形態の長さＬ１は、後述のポリマーモデル７の拡がりを示す量である慣性半径（図示省略）の３倍以上が望ましい。これにより、後述の粗視化分子動力学計算では、周期境界条件による自己のイメージとの衝突の発生を防いで、ポリマーモデル７の空間的拡がりを適切に計算することができる。また、セル４の大きさは、例えば１気圧で安定な体積に設定される。これにより、セル４は、解析対象の高分子材料の少なくとも一部の体積を定義することができる。
【００２７】
図６は、図５のＡ部拡大図である。セル４の内部には、例えば、立方体状に区分された複数の小領域９が定義されている。各小領域９には、節点９ｔが設定されている。このような小領域９は、後述する分子動力学計算において、フィラーモデル６の小粒子１２（図９に示す）、及び、ポリマーモデル７の粗視化粒子１５（図９に示す）の追跡等に用いられる。小領域９の１辺の長さＬ２は、例えば、粗視化粒子１５の直径に対して０．１〜５倍程度に設定されるのが望ましい。セル４は、コンピュータ１に記憶される。
【００２８】
次に、本実施形態のモデル設定工程Ｓ１では、図５に示したセル４の内部に、フィラーをモデル化したフィラーモデル６が定義される（工程Ｓ１２）。本実施形態の工程Ｓ１２では、各高分子材料モデル１０（本例では、第１高分子材料モデル１０Ａ及び第２高分子材料モデル（図示省略））のセル４毎に、フィラーモデル６がそれぞれ定義される。
【００２９】
本実施形態のフィラーモデル６は、セル４の内部で凝集した複数のフィラー粒子モデル１１によって定義されている。本実施形態では、実際の高分子材料の電子線透過画像のフィラーの一次粒子の位置に基づいて、フィラー粒子モデル１１が配置されている。これにより、フィラーモデル６は、実際のフィラーの形状を精度よく表現することができる。図７は、フィラーモデル６のフィラー粒子モデル１１の拡大図である。
【００３０】
各フィラー粒子モデル１１は、複数の小粒子１２を含んで構成されている。小粒子１２は、後述の粗視化分子動力学計算において、運動方程式の質点として取り扱われる。即ち、小粒子１２には、質量、直径、電荷又は初期座標などのパラメータが定義される。
【００３１】
各フィラー粒子モデル１１には、隣接する小粒子１２、１２間の相対位置を固定する拘束条件が定義されてもよいし、隣接する小粒子１２、１２間を拘束する結合鎖モデル（図示省略）が定義されても良い。これにより、フィラーモデル６は、後述の粗視化分子動力学計算において、フィラー粒子モデル１１の形状が維持され、フィラーモデル６の形状を実際のフィラーの形状に近似させることができる。
【００３２】
小粒子１２、１２間を拘束する結合鎖モデル（図示省略）は、例えば、論文（ Kurt Kremer & Gary S. Grest 著、「Dynamics of entangled linear polymer melts: A molecular-dynamics simulation」、J. Chem Phys. vol.92, No.8, 15 April 1990、p5057-5086）に記載されているレナードジョーンズポテンシャルとＦＥＮＥポテンシャルとの和で定義することができる。また、各ポテンシャルに定義される定数については、上記論文に基づいて、適宜設定することができる。
【００３３】
フィラー粒子モデル１１を構成する複数の小粒子１２のうち、フィラー粒子モデル１１の外面を構成する小粒子１２には、例えば、官能基をモデル化した官能基モデル（図示省略）が設けられてもよい。これにより、後述の粗視化分子動力学計算、及び、高分子材料の変形計算において、官能基によって変化するフィラーとポリマーとの相互作用を考慮することができる。フィラーモデル６は、コンピュータ１に記憶される。
【００３４】
次に、本実施形態のモデル設定工程Ｓ１では、図５に示したセル４の内部に、ポリマーをモデリングしたポリマーモデル７が定義される（工程Ｓ１３）。本実施形態の工程Ｓ１３では、各高分子材料モデル１０（本例では、第１高分子材料モデル１０Ａ及び第２高分子材料モデル（図示省略））のセル４毎に、ポリマーモデル７がそれぞれ定義される。
【００３５】
本実施形態の工程Ｓ１３では、セル４の内部において、フィラーモデル６が配置されてない領域に、少なくとも一つ（例えば、１０個〜１，０００，０００個）のポリマーモデル７が配置される。これにより、工程Ｓ１３では、フィラーモデル６との重なりを回避しながら、ポリマーモデル７を定義することができる。
【００３６】
図８は、ポリマーモデル７の一例を示す概念図である。本実施形態のポリマーモデル７は、図２に示したポリマーの分子構造を、複数の粗視化粒子１５でモデリングしたものである。隣接する粗視化粒子１５、１５の間には、結合鎖モデル１６で連結されている。
【００３７】
粗視化粒子１５は、ポリマーのモノマー又はモノマーの一部分をなす構造単位を置換したものである。ポリマーがポリイソプレンである場合には、上記論文に基づいて、例えば１．７３個分のモノマー３（図２に示す）を構造単位として、１個の粗視化粒子１５に置換される。これにより、各ポリマーモデル７には、複数（例えば、１０〜５０００個）の粗視化粒子１５が設定される。
【００３８】
粗視化粒子１５は、後述の粗視化分子動力学計算において、運動方程式の質点として取り扱われる。即ち、粗視化粒子１５には、例えば、質量、直径、電荷又は初期座標などのパラメータが定義される。
【００３９】
図９は、フィラーモデル６、及び、ポリマーモデル７の相互作用ポテンシャルの一例を説明する概念図である。結合鎖モデル１６は、粗視化粒子１５、１５間に、伸びきり長が設定された相互作用ポテンシャルＰ１によって定義される。相互作用ポテンシャルＰ１については、適宜定義することができる。相互作用ポテンシャルＰ１には、例えば、従来と同様に、レナードジョーンズポテンシャルとＦＥＮＥポテンシャルとの和で定義することができる。各ポテンシャルに定義される定数については、上記論文に基づいて、適宜設定することができる。これにより、粗視化粒子１５が伸縮自在に拘束された直鎖状のポリマーモデル７を定義することができる。ポリマーモデル７は、コンピュータ１に記憶される。
【００４０】
次に、本実施形態のモデル設定工程Ｓ１では、隣接するフィラーモデル６及びポリマーモデル７に、相互作用ポテンシャルが定義される（工程Ｓ１４）。本実施形態の工程Ｓ１４では、各高分子材料モデル１０（本例では、第１高分子材料モデル１０Ａ及び第２高分子材料モデル（図示省略））において、相互作用ポテンシャルがそれぞれ定義される。
【００４１】
本実施形態の工程Ｓ１４では、図９に示されるように、フィラーモデル６の小粒子１２、又は、ポリマーモデル７の粗視化粒子１５に、下記の相互作用ポテンシャルＰ２〜Ｐ４が定義される。
相互作用ポテンシャルＰ２：フィラーモデル６の小粒子１２と
フィラーモデル６の小粒子１２との間
相互作用ポテンシャルＰ３：ポリマーモデル７の粗視化粒子１５と
ポリマーモデル７の粗視化粒子１５との間
相互作用ポテンシャルＰ４：フィラーモデル６の小粒子１２と
ポリマーモデル７の粗視化粒子１５との間
【００４２】
上記相互作用ポテンシャルＰ２〜Ｐ４は、従来と同様に、レナードジョーンズポテンシャルで定義することができる。相互作用ポテンシャルＰ２〜Ｐ４の定数については、上記論文に基づいて、適宜設定することができる。これらの相互作用ポテンシャルＰ２〜Ｐ４は、コンピュータ１に記憶される。
【００４３】
次に、本実施形態のモデル設定工程Ｓ１では、高分子材料モデル１０（本例では、第１高分子材料モデル１０Ａ及び第２高分子材料モデル（図示省略））毎に、フィラーとポリマーとの間の相互作用ポテンシャルの強度を互いに異ならせる（工程Ｓ１５）。工程Ｓ１５では、フィラーモデル６の小粒子１２とポリマーモデル７の粗視化粒子１５との間に定義された相互作用ポテンシャルＰ４の強度を、高分子材料モデル１０毎に互いに異ならせている。
【００４４】
本実施形態の工程Ｓ１５では、高分子材料モデル１０（本例では、第１高分子材料モデル１０Ａ及び第２高分子材料モデル（図示省略））毎に、相互作用ポテンシャルＰ４として定義されているレナードジョーンズポテンシャルのポテンシャルの深さεを互いに異ならせている。上述したように、本実施形態の第１高分子材料は、第２高分子材料よりも物理吸着が弱いものとして区分されるため、第１高分子材料モデル１０Ａの相互作用ポテンシャルＰ４（即ち、レナードジョーンズポテンシャル）のポテンシャルの深さεが、第２高分子材料モデルの相互作用ポテンシャルＰ４のポテンシャルの深さεよりも小さく設定される。これにより、モデル設定工程Ｓ１では、フィラーとポリマーとの間の相互作用ポテンシャルの強度が互いに異なる高分子材料モデル１０をそれぞれ設定することができる。各高分子材料モデル１０は、コンピュータ１に記憶される。なお、第１高分子材料モデル１０Ａ、及び、第２高分子材料モデルの相互作用ポテンシャルＰ４のポテンシャルの深さεについては、適宜設定することができ、本実施形態では、以下のように設定される。
第１高分子材料モデルのポテンシャルの深さε：０．５
第２高分子材料モデルのポテンシャルの深さε：２．０
【００４５】
次に、本実施形態のシミュレーション方法は、コンピュータ１が、相互作用ポテンシャルに基づいて、各高分子材料モデル１０（本例では、第１高分子材料モデル１０Ａ及び第２高分子材料モデル（図示省略））の構造緩和を計算する（緩和工程Ｓ２）。本実施形態の緩和工程Ｓ２では、粗視化分子動力学計算に基づいて、各高分子材料モデル１０の構造緩和がそれぞれ計算される。
【００４６】
本実施形態の粗視化分子動力学計算では、例えば、セル４について所定の時間、フィラーモデル６、及び、ポリマーモデル７が古典力学に従うものとして、ニュートンの運動方程式が適用される。そして、各時刻でのフィラーモデル６、及び、ポリマーモデル７の動きが、シミュレーションの単位時間毎に追跡される。構造緩和の計算は、例えば（株）ＪＳＯＬ社製のソフトマテリアル総合シミュレーター（Ｊ−ＯＣＴＡ）に含まれるＣＯＧＮＡＣ、又は、ＶＳＯＰを用いて処理することができる。
【００４７】
ところで、従来のシミュレーション方法では、上記論文等と同様に、各高分子材料モデル１０の体積、温度、並びに、フィラー及びポリマーを構成する粒子の数（フィラーモデル６の小粒子１２、及び、ポリマーモデル７の粗視化粒子１５）が一定となる条件（即ち、ＮＶＴ）の下で、各高分子材料モデル１０の構造緩和が計算されていた。なお、各高分子材料モデル１０は、体積、温度及び粒子数が同一に設定されている。このため、従来のシミュレーション方法の構造緩和計算では、各高分子材料モデル１０において、フィラーに対するポリマーの分布（即ち、フィラーモデル６に対するポリマーモデル７の粗視化粒子１５の分布）がそれぞれ変化し、各高分子材料モデル１０の内部応力がそれぞれ異なってしまう。このような内部応力がそれぞれ異なる各高分子材料モデル１０を用いて変形計算が行われた場合、変形に伴う物理量を精度良く比較できないという問題がある。
【００４８】
本実施形態の緩和工程Ｓ２は、各高分子材料モデル１０（本例では、第１高分子材料モデル１０Ａ及び第２高分子材料モデル（図示省略））の圧力、温度、並びに、フィラー及びポリマーを構成する粒子（フィラーモデル６の小粒子１２、及び、ポリマーモデル７の粗視化粒子１５）の数が互いに同一となる条件の下で、構造緩和が計算される。図１０は、緩和工程Ｓ２の処理手順の一例を示すフローチャートである。
【００４９】
本実施形態の緩和工程Ｓ２では、先ず、モデル設定工程Ｓ１で設定された高分子材料モデル１０（本例では、第１高分子材料モデル１０Ａ及び第２高分子材料モデル（図示省略））から、一つの高分子材料モデル１０（以下、単に、「基準高分子材料モデル」ということがある。）が選択される（工程Ｓ２１）。本実施形態の基準高分子材料モデルとしては、第１高分子材料モデル１０Ａが選択される場合が例示されるが、このような態様に限定されるわけではなく、例えば、第２高分子材料モデルが選択されてもよい。
【００５０】
次に、本実施形態の緩和工程Ｓ２では、選択された一つの高分子材料モデル（基準高分子材料モデル）１０について、従来のシミュレーション方法と同様に、構造緩和が計算される（工程Ｓ２２）。工程Ｓ２２では、基準高分子材料モデル（第１高分子材料モデル１０Ａ）について、高分子材料モデル１０の体積、温度、並びに、フィラー及びポリマーを構成する粒子の数（フィラーモデル６の小粒子１２の数、及び、ポリマーモデル７の粗視化粒子１５の数）が一定となる条件（即ち、ＮＶＴアンサンブル）の下で、構造緩和が計算される。
【００５１】
工程Ｓ２２では、従来のシミュレーション方法と同様の手順に基づいて、基準高分子材料モデル（本例では、第１高分子材料モデル１０Ａ）において、フィラーモデル６、及び、ポリマーモデル７の初期配置が十分に緩和されるまで計算される。これにより、工程Ｓ２２では、フィラーモデル６、及び、ポリマーモデル７の平衡状態（構造が緩和した状態）を、確実に計算することができる。構造緩和された基準高分子材料モデルは、コンピュータ１に記憶される。
【００５２】
次に、本実施形態の緩和工程Ｓ２では、構造緩和が計算されていない他の高分子材料モデル１０（本例では、第２高分子材料モデル（図示省略））が選択され（工程Ｓ２３）、選択された高分子材料モデル１０の構造緩和が計算される（工程Ｓ２４）。
【００５３】
工程Ｓ２４では、選択された高分子材料モデル１０（本例では、第２高分子材料モデル（図示省略））について、基準高分子材料モデル（本例では、第１高分子材料モデル１０Ａ）の圧力、温度、並びに、フィラー及びポリマーを構成する粒子の数（フィラーモデル６の小粒子１２の数、及び、ポリマーモデル７の粗視化粒子１５の数）と同一となる条件の下で、構造緩和が計算される。
【００５４】
上述したように、モデル設定工程Ｓ１において、各高分子材料モデル１０（本例では、第１高分子材料モデル１０Ａ及び第２高分子材料モデル（図示省略））は、体積、温度及び粒子数が同一に設定されている。このため、工程Ｓ２４では、選択された高分子材料モデル１０（本例では、第２高分子材料モデル）について、基準高分子材料モデル（本例では、第１高分子材料モデル１０Ａ）とは異なる体積に変化させながら構造緩和が計算されることで、基準高分子材料モデルの圧力、温度、並びに、フィラー及びポリマーを構成する粒子の数と同一にしている。
【００５５】
図１１は、第１高分子材料モデル１０Ａの体積に対する第２高分子材料モデル（図示省略）の体積の比率（体積分率）と、時刻との関係を示すグラフである。図１２は、第１高分子材料モデルの平衡内部応力及び第２高分子材料モデル（図示省略）の内部応力と、時刻との関係を示すグラフである。本実施形態では、選択された高分子材料モデル１０（本例では、第２高分子材料モデル）の相互作用ポテンシャルＰ４の強度が、基準高分子材料モデル（本例では、第１高分子材料モデル１０Ａ）の相互作用ポテンシャルＰ４の強度に比べて大きい。このような場合、工程Ｓ２４では、図１１に示されるように、選択された高分子材料モデル（第２高分子材料モデル）の体積を、基準高分子材料モデル（第１高分子材料モデル１０Ａ）の体積よりも小さくして、構造緩和が計算される。これにより、選択された高分子材料モデル１０（第２高分子材料モデル）のフィラーに対するポリマーの分布と、基準高分子材料モデル（第１高分子材料モデル１０Ａ）のフィラーに対するポリマーの分布とが異なっていても、図１２に示されるように、選択された高分子材料モデル（第２高分子材料モデル（図示省略））の内部応力を、基準高分子材料モデル（第１高分子材料モデル１０Ａ）の内部応力に近似させることができる。
【００５６】
工程Ｓ２４では、選択された高分子材料モデル１０（本例では、第２高分子材料モデル（図示省略））において、フィラーモデル６、及び、ポリマーモデル７の初期配置が十分に緩和されるまで計算される。選択された高分子材料モデル１０の圧力の制御は、例えば、アンダーソン（Andersen）の圧力制御や、パリネロ・ラーマン（Parrinello-Rahman）等の圧力制御に基づいて実施される。これにより、工程Ｓ２４では、上記条件に基づいて、フィラーモデル６、及び、ポリマーモデル７の平衡状態（構造が緩和した状態）を、確実に計算することができる。構造緩和された高分子材料モデル１０（第２高分子材料モデル）は、コンピュータ１に記憶される。
【００５７】
次に、本実施形態の緩和工程Ｓ２では、全ての高分子材料モデル１０（本例では、第１高分子材料モデル１０Ａ及び第２高分子材料モデル（図示省略））の構造緩和が計算されたか否かが判断される（工程Ｓ２５）。工程Ｓ２５において、全ての高分子材料モデル１０の構造緩和が計算された場合（工程Ｓ２５において、「Ｙ」）、次の工程Ｓ３（図３に示す）が実施される。他方、工程Ｓ２５において、全ての高分子材料モデル１０の構造緩和が計算されていないと判断された場合（工程Ｓ２５において、「Ｎ」）、工程Ｓ２３〜工程Ｓ２５が再度実施される。これにより、緩和工程Ｓ２では、各高分子材料モデル１０の圧力、温度、並びに、フィラー及びポリマーを構成する粒子の数（フィラーモデル６の小粒子１２の数、及び、ポリマーモデル７の粗視化粒子１５の数）が互いに同一となる条件の下で、構造緩和を計算することができる。緩和計算後の各高分子材料モデル１０の内部応力を近似させることができる。
【００５８】
次に、本実施形態のシミュレーション方法では、コンピュータ１が、構造緩和後の各高分子材料モデル１０（本例では、第１高分子材料モデル１０Ａ及び第２高分子材料モデル（図示省略））をそれぞれ変形させて、各高分子材料モデル１０の変形に伴う物理量を計算する（工程Ｓ３）。
【００５９】
工程Ｓ３では、図５に示されるように、各高分子材料モデル１０（本例では、第１高分子材料モデル１０Ａ及び第２高分子材料モデル（図示省略））を、予め定められた方向に引っ張る単軸引張試験が計算される。高分子材料モデル１０を引っ張る方向については、適宜選択することができ、本例では、ｚ軸方向に引っ張っている。
【００６０】
高分子材料モデル１０の変形計算は、例えば、特許文献（特開２０１６−８１２９７号公報）に記載された内容の手順に従い、ｚ軸方向において、高分子材料モデル１０の一端（図５に示したセル４の一方側の面５ａ）、及び、高分子材料モデル１０の他端（図５に示したセル４の他方側の面５ｂ）が互いに離間するように、高分子材料モデル１０の伸長が計算される。高分子材料モデル１０の伸長計算は、分子動力学計算に基づいて計算される。
【００６１】
図１３（ａ）は、変形計算前の高分子材料モデル１０の一部を示す概念図である。図１３（ｂ）は、変形計算後の高分子材料モデル１０の一部を示す概念図である。高分子材料モデル１０の変形計算前において、セル４の内部の各小領域９には、上述した構造緩和計算により、フィラーモデル６（図示省略）、及び、ポリマーモデル７のいずれかが配置されている。
【００６２】
工程Ｓ３では、高分子材料モデル１０の伸長計算により、フィラーモデル６、及び、ポリマーモデル７の熱運動が計算される。このような熱運動は、高分子材料モデル１０に与えられた歪み、図９に示した上記相互作用ポテンシャルＰ１〜Ｐ４、及び、運動方程式に基づいて計算される。これにより、図１３（ｂ）に示されるように、セル４には、フィラーモデル６、ポリマーモデル７、及び、カップリング剤モデル８が配置されない小領域９が形成される。このような小領域９は、高分子材料モデル１０に形成された空孔（ボイド）２６として定義される。このような空孔２６により、高分子材料モデル１０の破壊が再現される。
【００６３】
変形条件としては、ポリマーモデル７が配置される領域の少なくとも一部に空孔２６を形成することができれば、適宜設定することができる。変形条件の一例としては、高分子材料モデル１０に与えられる歪みが０．１〜０．３程度であり、また、１歪み当たりの変形速度Ｖ１（図５に示す）が１０００〜１００００τ程度である。
【００６４】
工程Ｓ３では、各高分子材料モデル１０（本例では、第１高分子材料モデル１０Ａ及び第２高分子材料モデル（図示省略））の変形開始から変形終了までの間、分子動力学計算の単位ステップ（ＭＤステップ）毎に、各高分子材料モデル１０の変形に伴う物理量が計算される。物理量としては、適宜採用することができる。本実施形態の物理量は、高分子材料モデル１０の応力が含まれる。工程Ｓ３では、高分子材料モデル１０に与えられる歪毎に、高分子材料モデル１０の伸長方向の応力成分が計算される。これにより、本実施形態のシミュレーション方法では、各高分子材料モデル１０の応力と歪との関係を示す応力−歪曲線を取得することができる。
【００６５】
さらに、本実施形態の物理量としては、各高分子材料モデル１０（本例では、第１高分子材料モデル１０Ａ及び第２高分子材料モデル（図示省略））の空孔２６に関するパラメータを含んでいる。空孔２６に関するパラメータについては、例えば、高分子材料モデル１０の系全体に形成された空孔２６の合計体積や、各高分子材料モデル１０の体積に対する空孔２６の合計体積の割合（以下、単に「空孔の体積分率」ということがある。）など、適宜選択することができる。工程Ｓ３では、高分子材料モデル１０に与えられる歪毎に、高分子材料モデル１０の系全体に形成される空孔２６の合計体積が計算される。これにより、本実施形態のシミュレーション方法では、各高分子材料モデル１０の空孔２６に関するパラメータ（本例では、空孔の体積分率）と、歪との関係を示すグラフを取得することができる。空孔２６の体積分率は、上記変形計算によって生じたひずみによる高分子材料モデル１０（セル４）の体積変化量に対する空孔２６の体積の割合であってもよい。
【００６６】
本実施形態の工程Ｓ３では、構造緩和後の内部応力を互いに近似させた各高分子材料モデル１０（本例では、第１高分子材料モデル１０Ａ及び第２高分子材料モデル（図示省略））の変形が計算されるため、各高分子材料モデル１０の物理量の尺度を互いに揃えることができる。従って、本実施形態のシミュレーション方法では、フィラーとポリマーとの相互作用を考慮した各高分子材料モデル１０の変形を、物理量の尺度を揃えて計算することができるため、各高分子材料モデル１０の変形に伴う物理量を精度良く比較することができる。各高分子材料モデル１０の物理量は、コンピュータ１に記憶される。
【００６７】
次に、本実施形態のシミュレーション方法では、コンピュータ１が、各高分子材料モデル１０（本例では、第１高分子材料モデル１０Ａ及び第２高分子材料モデル（図示省略））の物理量を比較する（工程Ｓ４）。工程Ｓ４では、各高分子材料モデル１０について、応力及び空孔に関するパラメータがそれぞれ比較される。これにより、本実施形態のシミュレーション方法では、フィラーとポリマーとの相互作用（本例では、物理吸着）が互いに異なる高分子材料について、それらの性能（例えば、破壊特性等）を評価することができるため、高分子材料の開発に役立つ。また、本実施形態のシミュレーション方法では、様々な相互作用を設定して、高分子材料モデル１０の物理量を評価するパラメータスタディにより、高分子材料の開発の方向性を提示することが可能となる。
【００６８】
次に、本実施形態のシミュレーション方法では、コンピュータ１が、各高分子材料モデル１０（本例では、第１高分子材料モデル１０Ａ及び第２高分子材料モデル（図示省略））について、物理量が良好な高分子材料モデル１０が存在するか否かを判断する（工程Ｓ５）。工程Ｓ５では、予め定められた閾値に基づいて、各高分子材料モデル１０の物理量の良否が判断される。本実施形態では、各高分子材料モデル１０の応力が、閾値よりも大きい場合に、耐破壊性能に優れていると判断される。さらに、本実施形態では、各高分子材料モデル１０の空孔２６に関するパラメータが閾値よりも小さい場合に、耐破壊性能に優れていると判断される。なお、閾値については、例えば、高分子材料に求められる耐破壊性能等に応じて適宜設定することができる。
【００６９】
工程Ｓ５において、物理量が良好な高分子材料モデル１０が存在すると判断された場合（工程Ｓ５において、「Ｙ」）、物理量が良好な高分子材料モデル１０に設定された諸条件に基づいて、高分子材料が製造される（工程Ｓ６）。他方、工程Ｓ５において、物理量が良好な高分子材料モデル１０が存在しないと判断された場合（工程Ｓ５において、「Ｎ」）、フィラーとポリマーとの相互作用（本例では、物理吸着）が変更され（工程Ｓ７）、工程Ｓ２〜工程Ｓ５が再度実施される。これにより、本実施形態のシミュレーション方法は、フィラーとポリマーとの相互作用に基づいて、耐破壊性能に優れる高分子材料を開発することができる。
【００７０】
本実施形態のシミュレーション方法では、図５に示したセル４の内部に、フィラーモデル６及びポリマーモデル７が定義されたが、このような態様に限定されない。例えば、セル４の内部には、ポリマーにフィラーを結合させるためのカップリング剤をモデリングしたカップリング剤モデル（図示省略）がさらに定義されてもよい。
【００７１】
高分子材料モデル１０には、隣接するポリマーモデル７を連結するための架橋モデル（図示省略）が定義されてもよい。この場合、予め定められた架橋点に基づいて、ポリマーモデル７の粗視化粒子１５、１５間を架橋モデルで連結させるのが望ましい。これにより、この実施形態のシミュレーション方法では、架橋された高分子材料モデル１０の変形に伴う物理量を、精度よく比較することができる。
【００７２】
さらに、カップリング材料モデルは、例えば、ポリマーモデル７と同様に、複数の粗視化粒子と、隣接する粗視化粒子を結合する結合鎖モデルとで定義することができる。また、少なくとも一部のフィラーモデル６及びポリマーモデル７は、カップリング剤モデルを介して連結される。これにより、この実施形態のシミュレーション方法では、カップリング反応後の高分子材料モデル１０について、変形に伴う物理量を精度よく比較することができる。
【００７３】
以上、本発明の特に好ましい実施形態について詳述したが、本発明は図示の実施形態に限定されることなく、種々の態様に変形して実施しうる。
【実施例】
【００７４】
図３に示した処理手順に基づいて、フィラーとポリマーとの間の相互作用ポテンシャルの強度が互いに異なる高分子材料モデルが設定された（実施例、比較例）。高分子材料モデルとしては、上述した第１高分子材料モデル及び第２高分子材料モデルが設定された。そして、実施例及び比較例では、各高分子材料モデルの構造緩和が計算され、構造緩和後の各高分子材料モデルの変形に伴う物理量が計算された。
【００７５】
実施例の構造緩和を計算する工程は、各高分子材料モデルの圧力、温度、並びに、フィラー及びポリマーを構成する粒子の数が互いに同一となる条件の下で、各高分子材料モデルの構造緩和が計算された。一方、比較例では、特許文献（特開２０１６−８１２９７号公報）のシミュレーション方法と同様に、各高分子材料モデルの体積、温度、並びに、フィラー及びポリマーを構成する粒子の数が一定となる条件の下で、各高分子材料モデルの構造緩和が計算された。
【００７６】
実施例の変形計算は、下記の摩擦係数及び変形速度に基づいて、ｚ軸方向の引張計算が実施された。比較例の変形計算は、下記の初期摩擦係数及び初期変形速度に基づいて、ｚ軸方向の引張計算が実施された後に、摩擦係数を上記論文の摩擦係数と同一に設定し、かつ、高分子材料モデルの体積を固定して、９００τの緩和計算が実施された。実施例及び比較例の共通仕様については、次のとおりである。
セル：
一辺の長さＬ１：３５０σ
小領域の一辺の長さＬ２：１．５σ
フィラーモデル：
体積分率：１６．２％
フィラー粒子モデルの個数：１９７３個
フィラー粒子モデルの直径：２８．２σ
１つのフィラー粒子モデルを構成する小粒子：１６、５８９個
ポリマーモデル：
個数：１００、０００本
１つのポリマーモデルを構成する粗視化粒子：１０００個
架橋密度：０．００７２４５（１／σ3）
カップリング剤モデル：
１つのカップリング剤モデルを構成する粗視化粒子：１個
１つのフィラー粒子モデルあたりのカップリング剤モデル：８０個
変形計算：
歪み：０．１５
実施例：
摩擦係数：上記論文の摩擦係数と同一
変形速度：６３００τ／歪み
比較例：
初期摩擦係数：上記論文の１０倍
初期変形速度：６．３τ／歪み
【００７７】
図１２に示したように、実施例では、構造緩和後の各高分子材料モデルの内部応力を互いに近似させることができた。図１４は、実施例の応力と歪みとの関係を示すグラフである。図１５は、実施例の空孔の体積分率と歪みとの関係を示すグラフである。なお、体積分率は、ひずみ０．１５によるセルの体積変形量に対する空孔の体積の割合として計算された。図１４に示されるように、第２高分子材料モデルは、第１高分子材料モデルに比べて、応力が大きく計算された。図１５に示されるように、第２高分子材料モデルは、第１高分子材料モデルに比べて、５〜１０％歪において、空孔の体積分率が小さく計算された。このように、実施例は、各高分子材料モデルの物理量の尺度を互いに揃えることができたため、フィラーとポリマーとの相互作用（物理吸着）を考慮した各高分子材料モデルの変形を計算することができた。従って、実施例は、各高分子材料モデルの変形に伴う物理量を、精度良く比較することができた。
【００７８】
図１６は、比較例の各高分子材料モデルの内部応力と時刻との関係を示すグラフである。図１７は、比較例の空孔の体積分率と時刻との関係を示すグラフである。図１６に示されるように、比較例では、構造緩和後の各高分子材料モデルの内部応力が互いに大きく異なった。従って、比較例は、各高分子材料モデルの物理量の尺度を互いに揃えることができなかった。これにより、比較例は、フィラーとポリマーとの相互作用（物理吸着）を考慮した各高分子材料モデルの変形を計算できなかったため、図１７に示されるように、第１高分子材料モデルの物理量（空孔の体積分率）と、第２高分子材料モデルの物理量（空孔の体積分率）とが略同一となった。従って、比較例は、各高分子材料モデルの変形に伴う物理量を、精度良く比較することができなかった。
【符号の説明】
【００７９】
Ｓ１    相互作用ポテンシャルの強度が互いに異なる高分子材料モデルを設定する工程
Ｓ２    各高分子材料モデルの構造緩和を計算する工程
Ｓ３    各高分子材料モデルの変形に伴う物理量を計算する工程
【図１】
2019036253.tif 000003
【図２】
2019036253.tif 000004
【図３】
2019036253.tif 000005
【図４】
2019036253.tif 000006
【図５】
2019036253.tif 000007
【図６】
2019036253.tif 000008
【図７】
2019036253.tif 000009
【図８】
2019036253.tif 000010
【図９】
2019036253.tif 000011
【図１０】
2019036253.tif 000012
【図１１】
2019036253.tif 000013
【図１２】
2019036253.tif 000014
【図１３】
2019036253.tif 000015
【図１４】
2019036253.tif 000016
【図１５】
2019036253.tif 000017
【図１６】
2019036253.tif 000018
【図１７】
2019036253.tif 000019