JP7056137B2 - 特定物質の解析方法及び特定物質の解析用コンピュータプログラム - Google Patents

Description

本発明は、特定物質の解析方法及び特定物質の解析用コンピュータプログラムに関し、例えば、２以上の物質を含む特定物質の解析方法及び特定物質の解析用コンピュータプログラムに関する。

従来、分子動力学を用いた複合材料のシミュレーション方法が提案されている（例えば、特許文献１参照）。特許文献１に記載の複合材料のシミュレーション方法では、モデル作成領域内にポリマーモデル及びフィラーモデルを作成した後、フィラーモデル表面の結合位置にポリマーモデルを結合する。これにより、特許文献１に記載の複合材料の解析方法では、フィラー表面におけるポリマー粒子の結合状態が、複合材料の材料特性に与える影響を解析することが可能となる。

特開２０１５－０６４２４２号公報

ところで、タイヤの耐摩耗性能を向上させるゴム材料の開発を加速させるためには、ゴム材料の変形に伴うナノ構造の破壊のメカニズムを明らかにすることが一助となる。変形前後のゴム材料のナノ構造の破壊を解析することにより、実際のタイヤに用いられるフィラー充填ゴムの破断強度向上の材料開発を加速させることができる。

従来の分子動力学による数値解析では、ゴム材料の変形に伴うナノ構造の破壊を解析するためには、ゴム材料中のポリマー粒子の粒子間結合を破断して消去する必要がある。しかしながら、従来の分子動力学による数値解析では、粒子間結合の破断による力学応答の再現性が必ずしも十分でなく、ゴム材料のナノ構造の破壊のメカニズムの解析が困難であった。

本発明は、このような実情に鑑みてなされたものであり、分子鎖の破断に伴う力学応答の再現性に優れた特定物質の解析方法及び特定物質の解析用コンピュータプログラムを提供することを目的とする。

本発明に係る特定物質の解析方法は、コンピュータを用いた分子動力学法による特定物質の解析方法であって、特定物質をモデル化した特定物質モデルを含む特定物質の解析用モデルを作成するステップと、解析対象となる特定物質モデルに属し、粒子間結合で結合された少なくとも一対の粒子の粒子間距離に第１閾値を設定するステップと、前記特定物質の粒子の質量を変化させるステップと、前記粒子間距離が前記第１閾値以上の場合に、前記粒子間結合を破断処理して前記解析用モデルの数値解析を実行するステップと、を含むことを特徴とする。

本発明に係る特定物質の解析方法によれば、粒子間距離の伸長によって粒子間結合が破断処理される特定物質の粒子の質量が変量するので、破断処理された特定物質の分子鎖の伸縮速度を変化させることでき、粒子間結合の破断に伴う力学応答の再現性が向上する。そして、破断処理を疑似的に行うこともできるので、粒子間距離が閾値以上の領域においても、粒子間結合が破断によって消去されずに解析用モデルを数値解析することも可能となる。これにより、特定物質の解析方法は、破断に伴う粒子間結合の物理的な消滅を防ぐことができるので、粒子間結合の疑似的な破断を再現することが可能となる。

本発明の特定物質の解析方法においては、さらに、前記特定物質モデルを架橋させるステップを含むことが好ましい。この方法により、特定物質の解析方法は、架橋反応を介して特定物質を予め架橋した状態で解析用モデルの数値解析をできるので、特定物質の架橋が粒子間結合の破断に及ぼす影響を解析可能となる。

本発明の特定物質の解析方法においては、前記破断処理される前記粒子間結合を含む分子鎖に属する前記特定物質の粒子の質量を変化させることが好ましい。これにより、質量を変化させる特定物質の粒子の数を削減できるので、計算時間をより短縮することが可能となる。

本発明の特定物質の解析方法においては、前記破断処理される前記分子鎖の末端部の前記特定物質の粒子の質量を変化させることが好ましい。これにより、質量を変化させる特定物質の粒子の数を更に削減できるので、計算時間をより一層短縮することが可能となる。

本発明の特定物質の解析方法においては、前記破断処理と共に前記特定物質の質量を変化させることが好ましい。これにより、力学応答の再現性への影響を低減しつつ特定物質の質量の変化に要する演算負荷を軽減することができる。

本発明の特定物質の解析方法においては、前記粒子間距離が前記第１閾値より小さい前記第２閾値以上となった際に、前記特定物質の粒子の質量を変化させることが好ましい。これにより、破断処理時の特定物質の分子鎖の急激な収縮を確実に防ぐことができるので、数値解析時の力学応答の再現性がより向上する。

本発明の特定物質の解析方法においては、前記破断処理された分子鎖の形状の変化が集束してから、質量を変化させた前記特定物質の粒子の質量の変化量を小さくすることが好ましい。これにより、分子鎖の運動の集束後に特定物質の質量の変化量が減少するので、分子鎖の運動が適切な状態となって数値解析時の力学応答の再現性がより向上する。

本発明の特定物質の解析方法においては、前記破断処理後の経過時間に応じて前記特定物質の粒子の質量の変化量を変更することが好ましい。これにより、粒子間結合の破断処理後の分子鎖の急激な収縮を防ぎつつ分子鎖の運動が適切な状態に調整できるので、数値解析時の力学応答の再現性がより向上する。

本発明の特定物質の解析方法においては、質量分布を設けて前記特定物質の粒子の質量を変化させることが好ましい。これにより、粒子間結合の破断処理後の特定物質の粒子の運動が質量分布に伴って連続的に変化するので、分子鎖の運動が適切な状態となり、数値解析時の力学応答の再現性がより向上する。

本発明の特定物質の解析方法においては、前記質量の変化が、質量の増大であることが好ましい。これにより、粒子間結合の破断処理後の特定物質の分子鎖の急激な収縮を防ぐことができるので、数値解析時の力学応答の再現性がより向上する。

本発明の特定物質の解析方法においては、前記破断処理が、破断結合演算用関数を用いた疑似破断処理であることが好ましい。この方法により、特定物質の解析方法は、解析用モデルの数値解析時における解析対象となる一対の粒子の疑似的な切断を再現できるので、例えば、フィラー充填ゴムなどの実際の複合材料におけるポリマー分子の破断などを精度良く再現することが可能となる。

本発明の特定物質の解析方法においては、前記特定物質が、ポリマー及びフィラーを含むことが好ましい。この方法により、特定物質の解析方法は、例えば、フィラー充填ゴムなどの実際の複合材料における解析結果を精度良く再現することが可能となる。

本発明の特定物質の解析用コンピュータプログラムは、上記特定物質の解析方法をコンピュータに実行させることを特徴とする。

本発明の特定物質の解析用コンピュータプログラムによれば、粒子間距離の伸長によって粒子間結合が破断処理された特定物質の粒子の質量が変量するので、破断処理された特定物質の分子鎖の伸縮速度を変化させることでき、粒子間結合の破断に伴う力学応答の再現性が向上する。そして、破断処理を疑似的に行うこともできるので、粒子間距離が閾値以上の領域においても、粒子間結合が破断によって消去されずに解析用モデルを数値解析することが可能となる。これにより、特定物質の解析方法は、破断に伴う粒子間結合の物理的な消滅を防ぐことができるので、粒子間結合の疑似的な破断を再現することが可能となる。

本発明によれば、分子鎖の破断に伴う力学応答の再現性に優れた特定物質の解析方法及び特定物質の解析用コンピュータプログラムを実現できる。

図１は、本発明の実施の形態に係る複合材料の解析方法の一例の概略を示すフロー図である。 図２は、本発明の実施の形態に係る複合材料の解析方法で用いられる解析用モデルの一例を示す概念図である。 図３Ａは、ポリマー粒子の質量の変化の説明図である。 図３Ｂは、ポリマー粒子の質量の変化の説明図である。 図３Ｃは、ポリマー粒子の質量の変化の説明図である。 図３Ｄは、ポリマー粒子の質量の変化の説明図である。 図４Ａは、ポリマー粒子の質量の変化の説明図である。 図４Ｂは、ポリマー粒子の質量の変化の説明図である。 図５は、一対のポリマー粒子の粒子間距離と粒子間結合の結合エネルギーとの関係を示す図である。 図６は、一対のポリマー粒子の粒子間距離と粒子間結合の結合力との関係を示す図である。 図７は、一対のポリマー粒子の粒子間結合の粒子間距離と結合エネルギーとの関係を示す図である。 図８は、一対のポリマー粒子の粒子間結合の粒子間距離と結合エネルギーとの関係を示す図である。 図９は、一対のポリマー粒子の粒子間結合の粒子間距離と結合エネルギーとの関係を示す図である。 図１０Ａは、モデル作成領域内での粒子間結合の破断位置の説明図である。 図１０Ｂは、モデル作成領域内での粒子間結合の破断位置の説明図である。 図１０Ｃは、モデル作成領域内での粒子間結合の破断位置の説明図である。 図１１は、フィラーモデル表面からの距離と破断座標の座標分布との関係の一例を示す図である。 図１２Ａは、第１解析時間Ｔ１及び第２解析時間Ｔ２の２つの解析時間での可視化の例の説明図である。 図１２Ｂは、第１解析時間Ｔ１及び第２解析時間Ｔ２の２つの解析時間での可視化の例の説明図である。 図１３Ａは、第１解析時間Ｔ１～第３解析時間Ｔ３の３つの解析時間での可視化の例の説明図である。 図１３Ｂは、第１解析時間Ｔ１～第３解析時間Ｔ３の３つの解析時間での可視化の例の説明図である。 図１３Ｃは、第１解析時間Ｔ１～第３解析時間Ｔ３の３つの解析時間での可視化の例の説明図である。 図１４Ａは、第１解析時間Ｔ１～第３解析時間Ｔ３の３つの解析時間での可視化の例の説明図である。 図１４Ｂは、第１解析時間Ｔ１～第３解析時間Ｔ３の３つの解析時間での可視化の例の説明図である。 図１４Ｃは、第１解析時間Ｔ１～第３解析時間Ｔ３の３つの解析時間での可視化の例の説明図である。 図１５Ａは、第１解析時間Ｔ１及び第２解析時間Ｔ２の２つの解析時間を集約して可視化する例の説明図である。 図１５Ｂは、第１解析時間Ｔ１及び第２解析時間Ｔ２の２つの解析時間を集約して可視化する例の説明図である。 図１６Ａは、第１解析時間Ｔ１及び第２解析時間Ｔ２の２つの解析時間を集約して可視化する例の説明図である。 図１６Ｂは、第１解析時間Ｔ１及び第２解析時間Ｔ２の２つの解析時間を集約して可視化する例の説明図である。 図１７は、本発明の本実施の形態に係る複合材料の解析方法の他の例の概略を示すフロー図である。 図１８は、本実施の形態に係る複合材料の解析方法及び複合材料の解析方法を実行する解析装置の機能ブロック図である。 図１９は、本発明の実施例に係る応力歪曲線を示す図である。

本発明者は、従来の分子動力学を用いた複合材料の解析方法では、破断による力学応答の再現性が必ずしも十分でないことに着目した。そして、本発明者は、粒子間結合の破断処理後に特定物質の粒子の質量を変化させることにより、破断に伴う力学応答の再現性に優れた特定物質の解析方法を実現できることを見出した。さらに、本発明者は、分子結合を物理的に消滅させずに演算処理を行うことにより、分子結合の破断箇所の評価が可能となることも見出し、本発明を完成させるに至った。

以下、本発明の各実施の形態について、添付図面を参照して詳細に説明する。なお、本発明は、以下の各実施の形態に限定されるものではなく、適宜変更して実施可能である。なお、以下においては、解析対象となる特定物質が２以上の物質を含む複合材料である場合について説明しているが、本発明は、特定物質がポリマーなどの単一の物質である場合にも適用可能である。また、以下においては、複合材料がポリマー及びフィラーを含む例について説明しているが、本発明は、２種類の以上の物質を含有する複合材料にも適用可能である。また、本発明は、フィラー及びポリマー以外の物質を含有する複合材料にも適用可能である。

図１は、本実施の形態に係る複合材料の解析方法の概略を示すフロー図である。図１に示すように、本実施の形態に係る複合材料の解析方法は、コンピュータを用いた分子動力学法による複合材料の解析方法である。この特定物質の解析方法は、ポリマー（特定物質）をモデル化したポリマーモデル（特定物質モデル）及びフィラー（特定物質）をモデル化したフィラーモデル（特定物質モデル）を含む複合材料の解析用モデルを作成するステップＳＴ１１と、解析対象となるポリマーモデル又はフィラーモデルに属し、粒子間結合で結合された少なくとも一対の粒子の粒子間距離に閾値を設定するステップＳＴ１２と、粒子間距離が第１閾値以上の場合に、破断処理後にポリマー粒子の質量を変化させるステップＳＴ１３と、粒子間結合の破断処理をして解析用モデルの数値解析を実行するステップＳＴ１４と、を含む。

図２は、本実施の形態に係る複合材料の解析方法で用いられる解析用モデル１の一例を示す概念図である。図２に示すように、解析用モデル１は、例えば、一辺の長さが距離Ｌの略立方体形状の仮想空間であるモデル作成領域Ａ内でモデル化される。モデル作成領域Ａは、互いに直交するＸ軸、Ｙ軸及びＺ軸方向に広がる三次元空間となっている。解析用モデル１は、複数のフィラー粒子１１ａがモデル化された４つのフィラーモデル１１Ａ，１１Ｂ，１１Ｃ，１１Ｄと、複数のポリマー粒子２１ａ及び結合鎖２１ｂがモデル化された４つのポリマーモデル２１とを含む。なお、図２に示す例では、解析用モデル１が、４つのフィラーモデル１１Ａ，１１Ｂ，１１Ｃ，１１Ｄがモデル化された例について説明するが、モデル化されるフィラーモデルの数に制限はなく、解析用モデル１は、４つ未満のフィラーモデル１１を含んでいてもよく、４つを超えるフィラーモデル１１を含んでいてもよい。また、図２においては、４つのポリマーモデル２１のみを示しているが、解析用モデル１では、複数のポリマーモデル２１がモデル作成領域Ａ内の全域に亘って存在している。さらに、図２に示す例では、モデル作成領域Ａが、略直方体形状の仮想空間である例について示しているが、球状、楕円状、直方体形状、多面体形状など任意の形状であってもよい。

フィラーモデル１１は、複数のフィラー粒子１１ａがそれぞれ略球状体に集合した状態でモデル化される。また、フィラーモデル１１は、互いに所定間隔をとって離れた状態で配置されている。なお、フィラーモデル１１とは、相互に凝集した状態で外縁部が共有結合によって相互に連結されていてもよい。

フィラーとしては、例えば、カーボンブラック、シリカ、及びアルミナなどが含まれる。フィラー粒子１１ａは、複数のフィラーの原子が集合されてモデル化される。また、フィラー粒子１１ａは、複数のフィラー粒子１１ａが集合してフィラー粒子群を構成する。フィラー粒子１１ａは、複数のフィラー粒子１１ａ間の結合鎖（不図示）によって相対位置が特定されている。この結合鎖（不図示）は、フィラー粒子１１ａ間の結合距離である平衡長とばね定数とが定義されたバネとしての機能を有し、各フィラー粒子１１ａ間を拘束している。結合鎖は、フィラー粒子１１ａの相対位置及び捻り、曲げなどによって力が発生するポテンシャルが定義されているボンドである。フィラーモデル１１は、フィラーを分子動力学で取り扱うためのフィラー粒子１１ａの質量、体積、直径及び初期座標などを含む数値データである。フィラーモデル１１の数値データは、コンピュータに入力される。

ポリマーとしては、例えば、ゴム、樹脂、及びエラストマーなどが含まれる。ポリマー粒子２１ａは、複数のポリマーの原子が集合されてモデル化される。また、ポリマー粒子２１ａは、複数のポリマー粒子２１ａが集合してポリマー粒子群を構成する。ポリマーには、フィラーとの親和性を高める変性剤が必要に応じて配合される。この変性剤としては、例えば、水酸基、カルボニル基、及び原子団の官能基などが含まれる。ポリマーモデル２１は、複数のポリマー原子及び複数のポリマー原子の集合体であるポリマー粒子２１ａがモデル作成領域Ａ内に所定密度で充填されてモデル化される。ポリマー粒子２１ａは、複数のポリマー粒子２１ａ間の結合鎖２１ｂによって結合されて相対位置が特定されている。この結合鎖２１ｂは、ポリマー粒子２１ａ間の結合距離である平衡長とばね定数とが定義されたバネとしての機能を有し、各ポリマー粒子２１ａ間を拘束している。結合鎖２１ｂは、ポリマー粒子２１ａの相対位置及び捻り、曲げなどによって力が発生するポテンシャルが定義されているボンドである。また、結合鎖２１ｂは、複数のポリマー粒子２１ａが直列状に連結されてなるポリマーモデル２１間にも架橋結合（不図示）として結合されている。このポリマーモデル２１は、ポリマーを分子動力学で取り扱うための数値データ（ポリマー粒子２１ａの質量、体積、直径及び初期座標などを含む）である。ポリマーモデル２１の数値データは、コンピュータに入力される。

解析用モデル１は、分子動力学法による数値解析により各種物理量が取得される。数値解析としては、例えば、伸張解析、せん断解析などの変形解析及び緩和解析などの運動解析が挙げられる。これらの運動解析で取得する物理量は、運動解析の結果得られた変位などの値を用いてもよく、所定の演算処理を実行した歪みであってもよい。これらの中でも、運動解析としては、複合材料のコンパウンドの力学特性を解析可能となる観点から、変形解析が好ましい。

次に、本実施の形態に係る複合材料の解析方法について詳細に説明する。ステップＳＴ１１では、複数のフィラー粒子１１ａが集合してモデル化されたフィラーモデル１１及び複数のポリマー粒子２１ａが結合鎖２１ｂ（粒子間結合）を介して連結されてモデル化されたポリマーモデル２１を含む複合材料の解析用モデル１（図２参照）を作成する。

また、ステップＳＴ１１では、作成したフィラーモデル１１とポリマーモデル２１との間に相互作用を設定する。フィラーモデル１１とポリマーモデル２１との間の相互作用としては、例えば、分子間力及び水素結合などの引力及び斥力などの化学的な相互作用、及び共有結合などの物理的な相互作用が挙げられる。なお、フィラーモデル１１とポリマーモデル２１との間の相互作用は、フィラー粒子１１ａ間、ポリマー粒子２１ａ間及びフィラー粒子１１ａとポリマー粒子２１ａとの間に必要に応じて設定されるものであり、必ずしも全てのフィラー粒子１１ａ及びポリマー粒子２１ａに設定されるものではない。また、ポリマーモデル２１が複数の種類のポリマー粒子２１ａで構成されている場合には、複数の種類のポリマー粒子２１ａにそれぞれ相互作用を設定してもよい。また、複数の種類の各ポリマー粒子２１ａとフィラーモデル１１との相互作用は同一であってもよく、異なっていてもよい。例えば、ポリマー粒子Ａとフィラー粒子１１ａの相互作用とポリマー粒子Ｂとフィラー粒子１１ａの相互作用とは異なる相互作用を設定してもよい。

次に、ステップＳＴ１２では、ポリマー粒子２１ａの粒子間距離に所定の第１閾値を設定する。粒子間距離としては、ポリマー粒子２１ａを連結する結合鎖２１ｂの長さを用いてもよく、一対のポリマー粒子２１ａ間の直線距離を用いてもよい。なお、本実施の形態においては、ステップＳＴ１２において、解析対象となる一対のポリマー粒子２１ａの粒子間距離に所定の第１閾値を設定する例について説明するが、粒子間距離の第１閾値は、解析対象となる複合材料に応じて一対のフィラー粒子１１ａ間の粒子間距離に設定してもよい。

次に、ステップＳＴ１３では、ポリマーモデル２１のポリマー粒子２１ａの質量を変化させる。図３Ａ～図３Ｄは、ポリマー粒子２１ａの質量の変化の説明図である。なお、図３Ａ～図３Ｄにおいては、複数のポリマーモデル２１による三次元網目構造の一部を模式的に示し、一部のポリマー粒子２１ａのみを示している。

図３Ａ及び図３Ｂに示すように、ステップＳＴ１３では、結合鎖２１ｂ（粒子間結合）で結合された一対のポリマー粒子２１ａ－１，２１ａ－２を、質量変換粒子２２ａ－１，２２ａ－２に変換して質量を増大させる。この結果、図３Ｃに示すように、ステップＳＴ１４において、結合鎖２１ｂが破断処理されて破断結合鎖２１ｂｘとなった場合であっても、質量変換粒子２２ａ－１，２２ａ－２が属するポリマーモデル２１の分子鎖の質量が増大するので、一対の質量変換粒子２２ａ－１，２２ａ－２に作用する加速度が減少して間の粒子間距離の急激な増大を防ぐことができる。これに対して、図３Ｄに示すように、一対のポリマー粒子２１ａ－１，２１ａ－２を質量変換粒子２２ａ－１，２２ａ－２に変換しない場合には、結合鎖２１ｂが破断処理に伴い一対のポリマー粒子２１ａ－１，２１ａ－２の粒子間距離が急激に増大する。このように、ステップＳＴ１３においては、一対のポリマー粒子２１ａ－１，２１ａ－２を、質量変換粒子２２ａ－１，２２ａ－２に変換して質量を変化させることにより、ポリマーモデル２１の分子鎖の急激な収縮を防ぐことができるので、数値解析時の力学応答の再現性が向上する。

上述した実施の形態では、一対のポリマー粒子２１ａ－１，２１ａ－２の質量を増大させる例について説明したが、ポリマー粒子２１ａは、質量変換粒子２２ａ－１，２２ａ－２に変換して質量を減少させてもよい。この場合であっても、解析対象となる複合材料を構成する物質に応じて、結合鎖２１ｂの破断に伴うポリマーモデル２１の分子鎖の運動を変化させることができるので、数値解析時の力学応答の再現性が向上する。

上述した実施の形態では、結合鎖２１ｂで結合された一対のポリマー粒子２１ａ－１，２１ａ－２を質量変換粒子２２ａ－１，２２ａ－２に変換する例について説明したが、質量を変化させるポリマー粒子２１ａは適宜変更可能である。例えば、ステップＳＴ１３では、ポリマーモデル２１に属するポリマー粒子２１ａの全てを質量変換粒子２２ａに変換してもよく、ポリマーモデル２１に属するポリマー粒子２１ａの一部を質量変換粒子２２ａに変換してもよい。

図４Ａ及び図４Ｂは、ポリマー粒子２１ａの質量の変化の説明図である。なお、図４Ａ～図４Ｂにおいては、複数のポリマーモデル２１による三次元網目構造の一部を模式的に示し、一部のポリマー粒子２１ａのみを示している。

ステップＳＴ１３では、ステップＳＴ１４で破断される結合鎖２１ｂを含むポリマーモデル２１の分子鎖に属するポリマー粒子２１ａの質量を変化させることが好ましい。図４に示す例では、ポリマーモデル２１の三次元網目構造の中で、破断結合鎖２１ｘを含む分子鎖２１Ｓに属するポリマー粒子２１ａを変換する例を示している。これにより、ポリマーモデル２１に属する全てのポリマー粒子２１ａを質量変換粒子２２ａに変換する場合と比較して、質量変換粒子２２ａに変換するポリマー粒子２１ａの数を削減できるので、計算時間をより短縮することが可能となる。また、この場合には、ポリマーモデル２１の分子鎖における破断される末端部のポリマー粒子２１ａの質量を変化させることが好ましい。この場合、例えば、破断処理により破断結合鎖２１ｂｘとされる分子鎖２１Ｓの末端及び末端近傍のポリマー粒子２１ａを質量変換粒子２２ａとする。これにより、ポリマーモデル２１に属する全てのポリマー粒子２１ａを質量変換粒子２２ａに変換する場合と比較して、質量変換粒子２２ａに変換するポリマー粒子２１ａの数を削減できるので、計算時間をより一層短縮することが可能となる。

また、ステップＳＴ１３では、質量分布を設けてポリマー粒子２１ａの質量を変化させることが好ましい。図４Ｂに示す例では、破断結合鎖２１ｘとされる末端のポリマー粒子２１ａを質量変換粒子２２ａ－１，２２ａ－２とし、末端から離れた末端部のポリマー粒子２１ａを質量変換粒子２２ａ－３，２２ａ－４とする。そして、量変換粒子２２ａ－１，２２ａ－２の質量の変化量に対して、質量変換粒子２２ａ－３，２２ａ－４の質量の変化量を変化させて質量分布を設けている。この場合には、末端の質量変換粒子２２ａ－１，２２ａ－２の質量変化量に対して、末端から距離又は空間的に離れた質量変換粒子２２ａ－３，２２ａ－４の質量変化量を相対的に小さくすることが好ましい。これにより、結合鎖２１ｂの破断処理後のポリマー粒子２１ａの運動が質量分布に伴って連続的に変化するので、ポリマーモデル２１の分子鎖の運動が適切な状態となり、数値解析時の力学応答の再現性がより向上する。

なお、上述した実施の形態では、ポリマー粒子２１ａの結合鎖２１ｂを破断処理した後にポリマー粒子２１ａの質量を変化させる例について説明したが、ステップＳＴ１３は解析用モデル１の作成後の任意の時点で実行可能である。ステップＳＴ１３は、例えば、ステップＳＴ１１で解析用モデル１を作成した後に実行してもよく、ステップＳＴ１２で第１閾値Ｓを設定した後に実行してもよく、ステップＳＴ１４における結合鎖２１ｂの破断処理の直前に実行してもよい。

ステップＳＴ１３では、ステップＳＴ１４における結合鎖２１ｂの破断処理と共に実行してポリマー粒子２１ａの質量を変化させることが好ましい。これにより、力学応答の再現性への影響を低減しつつポリマー粒子２１ａの質量の変化に要する演算負荷を軽減することができる。

また、ステップＳＴ１３では、ポリマー粒子２１ａの粒子間距離に第１閾値より小さい第２閾値を設定し、粒子間距離が第２閾値以上となった際に実行することが好ましい。これにより、ポリマー粒子２１ａの質量の変化が破断処理前に確実に実現できるので、破断処理時のポリマー粒子２１ａの分子鎖の急激な収縮を確実に防ぐことができ、数値解析時の力学応答の再現性がより向上する。第２閾値は、ポリマー粒子２１ａの質量の変化に伴う計算負荷を低減する観点及び破断処理前にポリマー粒子２１ａの質量の変化を確実に実現する観点から、例えば、第１閾値に対して、０．９５倍以上１．０倍未満が好ましい。

また、ステップＳＴ１３では、破断処理によって破断されるポリマーモデル２１の分子鎖の形状の変化が集束してから、ポリマー粒子２１ａの質量と質量変換粒子２２ａの質量と変化量を小さくすることが好ましい。これにより、ポリマーモデル２１の分子鎖の運動の集束後にポリマー粒子２１ａの質量の変化量が減少するので、ポリマーモデル２１の分子鎖の運動が適切な状態となって数値解析時の力学応答の再現性がより向上する。ポリマーモデル２１の分子鎖の形状の変化の集束は、例えば、ポリマーモデル２１の末端間距離と慣性半径との関係に基づいて判断する。この場合、ポリマーモデル２１の分子鎖は、例えば、ポリマーモデル２１の末端間距離：慣性半径＝１：６となった際に形状の変化が集束したと判断できる。

また、ステップＳＴ１３では、ステップＳＴ１４での破断処理後の経過時間に応じて質量変換粒子２２ａの質量を変化させて、ポリマーモデル２１ａからの質量の変化量を変更することが好ましい。これにより、結合鎖２１ｂの破断処理後のポリマーモデル２１の分子鎖の急激な収縮を防ぎつつ、分子鎖の運動が適切な状態となるので、数値解析時の力学応答の再現性がより向上する。この場合、例えば、質量変換粒子２２ａの質量をポリマー粒子２１ａに対して増大させた場合には、経過時間に応じて質量変換粒子２２ａの質量を小さくしてもよい。これにより、分子鎖の破断に伴う分子鎖の収縮が集束しつつポリマーモデル２１の分子量が質量変化前の状態に徐々に戻るので、力学応答の再現性がより向上する。また、質量変換粒子２２ａの質量をポリマー粒子２１ａに対して増大させる場合には、例えば、指数関数などを用いて短時間で連続的に大きくしてもよい。この場合には、ポリマーモデル２１の分子鎖の質量が滑らかに増大するので、分子鎖の収縮を精度良く再現することができる。上述した質量の変化量は、瞬間的に質量を瞬間的に変化させるように段階的に変更してもよく、質量を徐々に変化させるように連続的に変更してもよい。また、ステップＳＴ１３では、ステップＳＴ１４での破断処理後に質量変換粒子２２ａをポリマー粒子２１ａに戻してもよい。

次に、ステップＳＴ１４では、ポリマー粒子２１ａの粒子間距離が第１閾値以上の場合に、ポリマー粒子２１ａの結合鎖２１ｂを破断処理して解析用モデル１の数値解析を実行する。ここでの数値解析としては、変温解析及び変圧解析などの数値解析、伸張解析、せん断解析などの変形解析及び緩和解析などの運動解析などの各種数値解析が挙げられる。また、破断処理は、結合鎖２１ｂ（粒子間結合）を実際に消去する破断処理をしてもよく、結合鎖２１ｂの結合エネルギー及び結合力の少なくとも一方を低下させる疑似破断処理をしてもよい。以下、疑似破断処理の場合について説明する。

ここで、ポリマー粒子２１ａの結合鎖２１ｂの結合エネルギー及び結合力と粒子間距離との関係について説明する。図５は、一対のポリマー粒子２１ａの粒子間距離と結合鎖２１ｂの結合エネルギーとの関係を示す図であり、図６は、一対のポリマー粒子２１ａの粒子間距離と結合鎖２１ｂの結合力との関係を示す図である。なお、図５においては、縦軸に一対のポリマー粒子２１ａの結合鎖２１ｂの結合エネルギーを示し、横軸に一対のポリマー粒子２１ａの粒子間距離を示している。図６においては、縦軸に一対のポリマー粒子２１ａの結合鎖２１ｂの結合力を示し、横軸に一対のポリマー粒子２１ａの粒子間距離を示している。

図５に示すように、解析用モデル１を用いた数値解析では、結合鎖２１ｂによって結合された一対のポリマーモデル粒子２１ａは、粒子間距離が変化することにより結合鎖２１ｂの結合エネルギーが増減する。結合鎖２１ｂの結合エネルギーは、一対のポリマー粒子２１ａの粒子間距離が増大するにつれて減少して最小値（図５の点Ｐ１参照）をとった後、更に粒子間距離が増大するにつれて増大してゆく（図５の破線Ｌ１参照）。

また、図６に示すように、解析用モデル１を用いた数値解析では、結合鎖２１ｂによって結合鎖２１ｂされた一対のポリマーモデル粒子２１ａは、粒子間距離が変化することにより結合鎖２１ｂの結合力が増減する。結合鎖２１ｂの結合力は、ポリマー粒子２１ａの粒子間距離が増大するにつれて増大し、変曲点（図６の点Ｐ２参照）を過ぎた後、更に増大してゆく（図６の破線Ｌ３参照）。

ところで、一般的な解析用モデル１を用いた数値解析では、一対のポリマー粒子２１ａ間の結合鎖２１ｂの破断を解析するためには、図５及び図６に示したような一対のポリマー粒子２１ａ間の結合鎖２１ｂの結合エネルギー及び結合力の増大を防ぐために、粒子間距離が一定以上となった際に結合鎖２１ｂなどを消去して破断する必要がある。このため、結合鎖２１ｂが破断された一対のポリマー粒子２１ａ間では、ポリマー粒子２１ａの連結関係が失われるので、破断された結合鎖２１ｂの再現及び特定は困難であった。

本発明者は、図５及び図６に示した一対のポリマー粒子２１ａの結合鎖２１ｂの結合エネルギー及び結合力と粒子間距離との関係に着目した。そして、本発明者は、一対のポリマー粒子２１ａ間の粒子間距離に所定の閾値Ｓを設けることを着想した。さらに、本発明者は、粒子間距離が第１閾値Ｓ以上となった場合にポリマー粒子２１ａ間の結合鎖２１ｂの結合エネルギー及び結合力の少なくとも一方を低下させて解析することにより、一対のポリマー粒子２１ａの結合鎖２１ｂを維持した状態で結合鎖２１ｂの破断を再現できることを見出し、本発明を完成させるに至った。

本実施の形態では、図５の実線Ｌ２に示すように、粒子間距離が所定の第１閾値Ｓ以上となった際に、一対のポリマー粒子２１ａ間の結合鎖２１ｂの結合エネルギーを低下させる。また、図６の実線Ｌ４に示すように、粒子間距離に設けた第１閾値Ｓ以上となった際に、一対のポリマー粒子２１ａ間の結合鎖２１ｂの結合力を低下させる。これにより、一対のポリマー粒子２１ａ間の粒子間距離が長くなり第１閾値Ｓ以上となっても、結合鎖２１ｂの結合エネルギー及び結合力が低下するので、第１閾値Ｓ以上でも結合鎖２１ｂが維持された状態で解析用モデル１の数値解析を継続することができる。この結果、一対のポリマー粒子２１ａ間の結合鎖２１ｂを維持しつつ結合鎖２１ｂの破断を再現することが可能となる。

破断処理が疑似破断処理である場合には、一対のポリマー粒子２１ａの粒子間距離に所定の第１閾値Ｓを設定し、粒子間距離が第１閾値Ｓ以上の場合には、粒子間距離が第１閾値Ｓ未満の場合に対して、結合鎖２１ｂの結合エネルギー及び結合力を低下させる破断結合演算用関数を用いて結合鎖２１ｂを演算する。当該破断結合演算用関数としては、例えば、下記式（１）に示すものが挙げられる。この破断結合演算用関数を用いることにより、図５及び図６に示した例と同様に、粒子間距離ｒが所定の第１閾値Ｓ未満の場合には、一対のポリマー粒子２１ａ間の粒子間距離ｒに応じて結合エネルギー及び結合力が増減し、粒子間距離ｒが所定の第１閾値Ｓ以上の場合には、結合エネルギー及び結合力がゼロとなる。なお、下記式（１）については、本発明の効果を奏する範囲で適宜変更可能である。なお、以下においては、結合エネルギーを低下させ例について説明するが、結合力を低下させる場合にも同様に実施可能である。

このように、破断処理が疑似破断処理である場合には、粒子間距離が第１閾値Ｓ以上の領域においては、結合鎖２１ｂの結合エネルギー及び結合力の少なくとも一方が低下するので、結合鎖２１ｂを破断することなく解析用モデル１の数値解析を継続することが可能となる。これにより、破断に伴う結合鎖２１ｂの物理的な消滅を防ぐことができるので、結合鎖２１ｂを物理的に消滅させずに疑似的な切断を再現することが可能となる。したがって、数値解析時に破断した結合鎖２１ｂの破断箇所を特定することが可能となる複合材料の解析方法を実現できる。

なお、図５に示した例では、第１閾値Ｓ以上で結合エネルギーをゼロにする例について説明したが、結合エネルギーは、解析用モデル１の運動解析が可能であれば、必ずしもゼロまで低下させる必要はない。第１閾値Ｓ以上での結合エネルギーは、例えば、結合エネルギーの最小点Ｐ１以下に低下させることが好ましい。同様に、図６に示した例では、第１閾値Ｓ以上で結合力をゼロにする例について説明したが、結合力は、解析用モデル１の数値解析が可能であれば、必ずしもゼロまで低下させる必要はない。第１閾値Ｓ以上での結合力は、例えば、結合力の変曲点Ｐ２以下に低下させることが好ましく、結合力を実質的に消失させる略ゼロまで低下させることがより好ましい。なお、ここでの略ゼロとは、実質的にゼロに近い数値範囲のことであり、若干のゼロ以上ゼロ以下の数値範囲を含むものとする。

また、破断結合演算用関数としては、第１閾値Ｓ以上で結合鎖２１ｂの結合エネルギー及び結合力をゼロとする破断結合演算用関数を用いて解析用モデル１の数値解析を実行することが好ましい。これにより、解析用モデル１の数値解析時における解析対象となる一対の粒子の疑似的な切断を精度良く再現できるので、例えば、フィラー充填ゴムなどの実際の複合材料におけるポリマー分子の破断などを精度良く再現することが可能となる。

また、上述した実施の形態では、上記式（１）に示す１つの破断結合演算用関数を用いて結合鎖２１ｂの結合エネルギー又は結合力を低下させる例について説明したが、複数の破断結合演算用関数を組み合わせて結合鎖２１ｂの結合エネルギー又は結合力の少なくとも一方を低下させてもよい。図７は、一対のポリマー粒子２１ａの結合鎖２１ｂの粒子間距離と結合エネルギーとの関係を示す図である。図７においては、縦軸に一対のポリマー粒子２１ａの結合鎖２１ｂの結合エネルギーを示し、横軸に一対のポリマー粒子２１ａの粒子間距離を示している。

図７に示す例では、粒子間距離に応じて結合鎖２１ｂの演算に用いる破断結合演算用関数を第１破断結合演算用関数と第２破断結合演算用関数との間で切替えて演算処理する。粒子間距離が所定の第１閾値Ｓ未満の場合には（図７の実線Ｌ５参照）、第１破断結合演算用関数としての下記式（２）に示す式を用いて結合鎖２１ｂの結合エネルギーを演算処理し、粒子間距離が所定の第１閾値Ｓ以上の場合には（図７の実線Ｌ６参照）、第２破断結合演算用関数としての下記式（３）に示す式を用いて結合鎖２１ｂの結合エネルギーを演算処理する。このように、第１破断結合演算用関数及び第２破断結合演算用関数を第１閾値Ｓを境に切り替えて用いることにより、図５及び図６に示した例と同様に、粒子間距離ｒが所定の第１閾値Ｓ未満の場合には、一対のポリマー粒子２１ａ間の粒子間距離ｒに応じて結合エネルギーが増減し、粒子間距離ｒが所定の第１閾値Ｓ以上の場合には、結合エネルギーがゼロとなる。なお、下記式（２）及び下記式（３）については、本発明の効果を奏する範囲で適宜変更可能である。

なお、上述した実施の形態では、一対のポリマー粒子２１ａの粒子間距離に特定数値の第１閾値Ｓを設定する例について説明したが、第１閾値としては、粒子間距離に所定の数値範囲を有する閾値範囲として設定してもよい。図８及び図９は、一対のポリマー粒子２１ａにおける粒子間距離と結合鎖２１ｂの結合エネルギーとの関係を示す図である。なお、図８及び図９においては、縦軸に結合鎖２１ｂの結合エネルギーを示し、横軸に一対のポリマー粒子２１ａの粒子間距離を示している。

図８に示す例では、粒子間距離の所定数値範囲に閾値範囲ＳＲを設定している。この場合には、上述した図５～図７に示した例と同様に、閾値範囲ＳＲのいずれかで結合エネルギーがゼロとなるようにする（図８の実線Ｌ７参照）。この場合には、図５に示したように、特定の第１閾値Ｓのみで結合エネルギーが低下するのではなく、閾値範囲ＳＲ内のいずれかの値で結合エネルギーが低下するので、結合鎖２１ｂの破断の発生確率を調整することが可能となる。この結果、必ずしも特定の結合長で分子鎖が切断されない実際の複合材料における結合鎖２１ｂの破断を精度よく再現することが可能となる。なお、閾値範囲ＳＲとしては、例えば、ポリマー粒子２１ａの平衡結合長の１．２倍以上に設定することが好ましく、１．５倍以上に設定することがより好ましい。また、結合エネルギーは、必ずしもゼロまで低下させる必要はなく、また低下後の結合エネルギーは、一定値となるようにしてもよく、所定範囲内で変動するようにしてもよい。

図９に示す例では、ステップＳＴ１３において、閾値範囲ＳＲの下限値から上限値に向けて結合鎖２１ｂの結合エネルギー及び結合力の少なくとも一方を減少させる破断結合演算用関数を用いて数値解析を実行する。この場合には、例えば、図９の実線Ｌ８～実線Ｌ１０に示すように、閾値範囲ＳＲ内での結合エネルギーの減少率を任意に調整した破断結合演算用関数を使用してもよい。これにより、閾値範囲ＳＲの下限値から上限値の範囲で徐々に結合エネルギーが減少するので、結合鎖２１ｂの破断の発生確率を精度よく調整することが可能となり、必ずしも特定の結合長で分子鎖が切断されない実際の複合材料における結合鎖２１ｂの破断をより高精度で再現することが可能となる。図９に示す例で用いられる破断結合演算用関数としては、例えば、実線Ｌ８～実線Ｌ１０に対応したｎ次曲線を示す関数及び指数関数などが挙げられる。

また、上述した実施の形態では、ポリマー粒子２１ａの粒子間距離が第１閾値Ｓ以上となると所定の破断結合演算用関数を適用して結合鎖２１ｂを演算する例について説明したが、粒子間距離の時間平均値が第１閾値以上となった場合に、結合鎖２１ｂの結合エネルギー及び結合力を低下させる破断結合演算用関数を用いて解析用モデル１の数値解析を実行してもよい。これにより、例えば、ブラウン運動などによって、一時的に粒子間距離が第１閾値以上となった一対のポリマー粒子２１ａにおける結合鎖２１ｂの結合エネルギー及び結合力を低下させる破断結合演算用関数を用いた演算を除外することができる。この結果、実際の複合材料における結合鎖２１ｂの破断を精度よく再現することが可能となる。時間平均値の第１閾値としては、評価時間間隔の１／１０００が好ましく、１／１００がより好ましい。

ところで、結合鎖２１ｂによって連結された一対のポリマー粒子２１ａが複数存在する場合には、解析用モデル１の数値解析後、複数の一対のポリマー粒子２１ａの粒子間距離が順次閾値Ｓ以上となり、結合鎖２１ｂの結合エネルギー及び結合力の少なくとも一方が順次低下した状態となる。このため、所定解析時間後の特定の解析時間に複数の一対のポリマー粒子２１ａの粒子間距離が閾値Ｓ以上となっていた場合には、結合鎖２１ｂの結合エネルギー及び結合力の少なくとも一方が低下したポリマーモデル２１の座標をそれぞれ特定しても、必ずしも十分にそれぞれの正確な破断位置を評価することができない場合がある。

そこで、上記実施の形態においては、粒子間距離が閾値Ｓ以上となった時点の結合鎖２１ｂの座標を破断位置として特定して評価してもよい。図１０Ａ～図１０Ｃは、モデル作成領域Ａ内での結合鎖２１ｂの破断位置の説明図である。なお、図１０Ａにおいては、第１解析時間Ｔ１の状態を示し、図１０Ｂにおいては、第２解析時間Ｔ２の状態を示し、図１０Ｃにおいては、第３解析時間Ｔ３の状態を示している。

本実施の形態では、複数の解析時間について、それぞれ結合鎖２１ｂの結合エネルギー及び結合力が低下した結合鎖２１ｂの位置を特定する。図１０Ａに示す例では、第１解析時間Ｔ１では、ポリマーモデル２１Ａは、フィラーモデル１１の近傍で粒子間距離が閾値Ｓ以上となり、結合エネルギー又は結合力が低下した破断結合鎖２１ｂｘが生じている。フィラーモデル１１の近傍に存在するポリマーモデル２１Ｂ，２１Ｃは、ポリマー粒子２１ａの粒子間距離が閾値Ｓ未満となり結合鎖２１ｂが残存している。この第１解析時間Ｔ１では、ポリマーモデル２１Ａの座標を破断座標Ｘ_１として特定する。

次に、所定時間経過後の第２解析時間Ｔ２では、ポリマーモデル２１Ｂは、フィラーモデル１１の近傍で粒子間距離が閾値Ｓ以上となり、結合エネルギー又は結合力が低下した破断結合鎖２１ｂｘが生じている。フィラーモデル１１の近傍に存在するポリマーモデル２１Ｃは、ポリマー粒子２１ａの粒子間距離が閾値Ｓ未満となり結合鎖２１ｂが残存している。一方、移動によりフィラーモデル１１から離れたポリマーモデル２１Ａは、結合エネルギー又は結合力が低下した破断結合鎖２１ｂｘが維持されている。この第２解析時間Ｔ２では、ポリマーモデル２１Ｂの座標を新たな破断座標Ｘ_２として特定し、第１解析時間Ｔ１で既に破断結合鎖２１ｂｘが生じたポリマーモデル２１Ａの現座標は破断座標として新たに特定しない。

さらに、所定時間経過後の第３解析時間Ｔ３では、ポリマーモデル２１Ｃは、フィラーモデル１１の近傍で粒子間距離が閾値Ｓ以上となり、結合エネルギー又は結合力が低下した破断結合鎖２１ｂｘが生じている。移動によりフィラーモデル１１表面から離れたポリマーモデル２１Ａ，２１Ｂは、ポリマー粒子２１ａの結合鎖２１ｂが閾値Ｓ以上となり破断結合２１ｂｘが維持されている。この第３解析時間Ｔ３では、ポリマーモデル２１Ｃの座標を破断座標Ｘ_３として新たに特定し、第１解析時間Ｔ１で既に破断結合鎖２１ｂｘが生じたポリマーモデル２１Ａの現座標及び第２解析時間Ｔ２で既に破断結合鎖２１ｂｘが生じたポリマーモデル２１Ｂの現座標は破断座標として特定しない。

このように、連続する第１解析時間Ｔ１～第３解析時間Ｔ３中に粒子間距離が閾値Ｓ以上となった破断座標Ｘ_１～Ｘ_３を順次特定することにより、例えば、図１１に示すように、フィラーモデル１１表面からの距離と破断座標の座標分布とが得られる。これにより、図１０Ａ～図１０Ｃに示した解析用モデル１の数値解析では、フィラーモデル１１からの距離が近くなるにつれて破断座標Ｘ_１～Ｘ_３の座標分布が増大することが分かる。この結果から、図１１に示すように、結合鎖２１ｂの破断されやすい場所を評価することができるので、フィラーモデル１１表面からの距離と破断確率の関係などを評価することが可能となる。

なお、図１０Ａ～図１０Ｃに示した例では、結合鎖２１ｂに代表点を設定して破断位置（破断座標）を特定してもよい。例えば、図１０Ａに示した例では、結合鎖２１ｂである結合鎖２１ｂにおける一対のポリマー粒子２１ａとの重心（中点）又はポリマー粒子２１ａの座標と重なる端点などを代表点として破断位置として特定する。これにより、長さが増大した結合鎖２１ｂの代表点の座標を破断位置として評価することができるので、破断位置の評価が容易となる。

また、図１０Ａ～図１０Ｃに示した例では、粒子間距離が所定値以上となった結合鎖２１ｂを可視化してもよい。これにより、疑似的に破断した一対のポリマー粒子２１ａ間の結合鎖２１ｂである破断結合鎖２１ｂｘを目視で確認することができるので、数値解析時に破断した結合鎖２１ｂの破断箇所の特定が容易となる。同様に、結合鎖２１ｂが破断して破断結合鎖２１ｂｘが生じた破断座標Ｘ_１～Ｘ_３を可視化してもよい。これにより、破断しやすい場所を目視で評価できるので、破断しやすい場所の評価が容易となる。また、代表点を可視化してもよい。これにより、長さが増大した結合鎖２１ｂの全体を可視化せずに代表点を可視化するので、破断した結合鎖２１ｂの確認が容易となる。結合鎖２１ｂの可視化は、例えば、破断結合鎖２１ｂｘ以外の結合鎖２１ｂを非表示としてもよく、結合鎖２１ｂの透明度を高めてもよく、破断結合鎖２１ｂｘの色及び太さを結合鎖２１ｂと変更して表示してもよい。これにより、破断結合鎖２１ｂｘを強調することができ、目視で容易に確認することが可能となる。

また、上述した実施の形態においては、可視化を数値解析中の複数の解析時間において実行してもよい。図１２Ａ及び図１２Ｂは、第１解析時間Ｔ１及び第２解析時間Ｔ２の２つの解析時間での可視化の例の説明図であり、図１３Ａ～図１３Ｃ及び図１４Ａ～図１４Ｃは、第１解析時間Ｔ１～第３解析時間Ｔ３の３つの解析時間での可視化の例の説明図である。なお、図１２Ａ～図１４Ｃにおいては、一対のポリマー粒子２１ａの粒子間距離が閾値Ｓ以上となった破断座標Ｘ_１～Ｘ３を略球状に模式的に示している。

図１２Ａ及び図１２Ｂに示す例では、図１２Ａに示すように、第１解析時間Ｔ１でフィラーモデル１１の近傍に破断座標Ｘ_１が発生して可視化された後、図１２Ｂに示すように、第２解析時間Ｔ２では、フィラーモデル１１の移動に伴って破断座標Ｘ_１がフィラーモデル１１との相対位置が保たれた状態で表示される。すなわち、第１解析時間Ｔ１で発生した破断座標Ｘ_１は、第１解析時間Ｔ１での絶対座標Ｘ_１１からフィラーモデル１１との相対位置が保たれた新たな破断座標Ｘ_１として表示される。これにより、第１解析時間Ｔ１から第２解析時間Ｔ２に解析時間が進展しても、フィラーモデル１１の周囲で発生した破断座標Ｘ_１の相対位置は維持されて表示される。

図１３Ａ～図１３Ｃに示す例では、図１３Ａに示すように、第１解析時間Ｔ１でフィラーモデル１１の近傍に破断座標Ｘ_１が発生して可視化される。そして、図１３Ｂに示すように、第２解析時間Ｔ２では、第１解析時間Ｔ１で発生した破断座標Ｘ_１がフィラーモデル１１との相対位置が保たれた状態でフィラーモデル１１と共に移動すると共に、第２解析時間Ｔ２で発生した破断座標Ｘ_２がフィラーモデル１１の近傍で新たに可視化される。さらに、第３解析時間Ｔ３では、第１解析時間Ｔ１で発生した破断座標Ｘ_１及び第２解析時間Ｔ２で発生した破断座標Ｘ_２がフィラーモデル１１との相対位置が保たれた状態でフィラーモデル１１と共に移動すると共に、第３解析時間Ｔ３で発生した破断座標Ｘ_３がフィラーモデル１１の近傍で新たに可視化される。このように可視化することにより、フィラーモデル１１の位置を第１解析時間Ｔ１～第３解析時間Ｔ３での座標に応じて変化させた場合であっても、フィラーモデル１１の周囲で発生した破断座標Ｘ_１～Ｘ_３の相対座標を保って表示することが可能となる。この結果、フィラーモデル１１の周囲のポリマーモデル２１の結合鎖２１ｂの破断を再現することが可能となる。

図１４Ａ～図１４Ｃに示す例では、モデル作成領域Ａ内でのフィラーモデル１１の位置は固定して表示される。まず、図１４Ａに示すように、第１解析時間Ｔ１でフィラーモデル１１の近傍に破断座標Ｘ_１が発生して可視化される。そして、図１４Ｂに示すように、第２解析時間Ｔ２では、フィラーモデル１１及び第１解析時間Ｔ１で発生した破断座標Ｘ_１の表示位置が維持された状態で、第２解析時間Ｔ２で発生した破断座標Ｘ_２が新たに可視化される。さらに、図１４Ｃに示すように、第３解析時間Ｔ３では、フィラーモデル１１、第１解析時間Ｔ１で発生した破断座標Ｘ_１及び第２解析時間Ｔ２で発生した破断座標Ｘ_２の表示位置が維持された状態で、第３解析時間Ｔ３で発生した破断座標Ｘ_３が新たに可視化される。このようにモデル作成領域Ａ内でフィラーモデル１１の位置を固定した状態で可視化することによっても、フィラーモデル１１の周囲で発生した破断座標Ｘ_１～Ｘ_３の相対座標を保って表示することが可能となる。この結果、フィラーモデル１１の周囲のポリマーモデル２１の結合鎖２１ｂの破断を再現することが可能となる。これらにより、複数の解析時間毎の結合鎖２１ｂの破断の進展を確認することが可能となる。

また、上記実施の形態においては、数値解析を解析用モデル１中に含まれる複数のポリマーモデル２１又はフィラーモデル１１についてそれぞれ実行し、得られた複数の数値解析の結果を集約して可視化して評価してもよい。図１５Ａ及び図１５Ｂは、第１解析時間Ｔ１及び第２解析時間Ｔ２の２つの解析時間を集約して可視化する例の説明図である。本実施の形態では、モデル作成領域Ａ内には、第１フィラーモデル１１Ａ及び第２フィラーモデル１１Ｂが存在し、第１フィラーモデル１１Ａ及び第２フィラーモデル１１Ｂのそれぞれに対して数値解析を実行する。

図１５Ａに示すように、第１解析時間Ｔ１では、第１フィラーモデル１１Ａの近傍でポリマーモデル２１の結合鎖２１ｂの破断が発生し、第１フィラーモデル１１Ａの近傍に１つの破断座標Ｘ_１が可視化される。また、図１５Ｂに示すように、第２解析時間Ｔ２では、第１フィラーモデル１１Ａ及び第２フィラーモデル１１Ｂがそれぞれモデル作成領域Ａ内で移動する。また、第２解析時間Ｔ２では、第１フィラーモデル１１Ａの近傍及び第２フィラーモデル１１Ｂの近傍でそれぞれ１つのポリマーモデル２１の結合鎖２１ｂの破断が発生し、第１フィラーモデル１１Ａの近傍に破断座標Ｘ_２が可視化されて新たに追加され、第２フィラーモデル１１Ｂの近傍に破断座標Ｘ_３が新たに可視化される。この結果、第２解析時間Ｔ２では、第１フィラーモデル１１Ａの近傍で発生した２つの破断座標Ｘ_１，Ｘ_２及び第２フィラーモデル１１Ｂの近傍で発生した１つの破断座標Ｘ_３の３つの破断座標Ｘ_１～Ｘ_３が集約して表示される。

このように、本実施の形態では、解析用モデル１中の２つの第１フィラーモデル１１Ａ及び第２フィラーモデル１１Ｂを１つのモデル作成領域Ａ内に集約して表示できるので、解析結果が迅速に得られると共に解析結果の理解が容易となる。

また、図１５Ａ及び図１５Ｂに示す例においては、数値解析を解析用モデル１中に含まれる複数のポリマーモデル２１又はフィラーモデル１１についてそれぞれ実行し、得られた複数の数値解析の結果を解析用モデル１中に指定した特定の第１フィラーモデル１１Ａに集約して可視化して評価してもよい。図１６Ａ及び図１６Ｂは、第１解析時間Ｔ１及び第２解析時間Ｔ２の２つの解析時間を集約して可視化する例の説明図である。

図１６Ａ及び図１６Ｂに示すように、本実施の形態では、モデル作成領域Ａ内には、図１５Ａ及び図１５Ｂに示した第１フィラーモデル１１Ａ（図１６Ａ及び図１６Ｂにおいて不図示）及び第２フィラーモデル１１Ｂ（図１６Ａ及び図１６Ｂにおいて不図示）の計算結果を集約して表示する代表フィラーモデル１１Ｅが存在する。この代表フィラーモデル１１Ｅは、図１５Ａ及び図１５Ｂに示した第１フィラーモデル１１Ａ及び第２フィラーモデル１１Ｂの近傍で発生する破断座標Ｘ_１～Ｘ_３を集約して表示するものである。なお、図１６Ａ及び図１６Ｂに示す例では、代表フィラーモデル１１Ｅは、モデル作成領域Ａ内の座標が固定されて表示される。

図１６Ａに示すように、第１解析時間Ｔ１では、第１フィラーモデル１１Ａの近傍にのみ破断座標Ｘ_１が可視化されるので、第１フィラーモデル１１Ａの近傍に発生した破断座標Ｘ_１が代表フィラーモデル１１Ｅに対応する座標に表示される。図１６Ｂに示すように、第２解析時間Ｔ２では、第１フィラーモデル１１Ａの近傍で発生した破断座標Ｘ_２及び第２フィラーモデル１１Ｂの近傍で発生した破断座標Ｘ_３が代表フィラーモデル１１Ｅに新たに追加される。この結果、代表フィラーモデル１１Ｅの周囲には、第１フィラーモデル１１Ａ及び第２フィラーモデル１１Ｂを含めたフィラーモデル１１全体の周囲に発生した破断座標Ｘ_１～Ｘ_３が集約して表示される。

このように、本実施の形態では、解析用モデル１中の２つの第１フィラーモデル１１Ａ及び第２フィラーモデル１１Ｂを特定の１つの代表フィラーモデル１１Ｅとしてモデル作成領域Ａ内に集約して表示できるので、解析結果が迅速に得られると共に解析結果の理解が容易となる。代表フィラーモデル１１Ｅとしては、解析用モデル１に含まれる任意のフィラーモデル１１を特定してもよく、また解析用モデル１とは別個に作成した新たなフィラーモデル１１などのその他のモデルを特定してもよい。また、代表フィラーモデル１１Ｅに破断座標の情報を集約するフィラーモデル１１としては、解析用モデル１中の全フィラーモデル１１としてもよく、解析用モデル１中に指定した複数のフィラーモデル群としてもよい。複数のフィラーモデル群としては、例えば、特定のフィラーモデル１１を基準として他のフィラーモデル１１との間の距離が所定の閾値以上離れた分散フィラーモデル群としてもよく、特定のフィラーモデル１１を基準として他のフィラーモデル１１との間の距離が所定の閾値未満の凝集フィラーモデル群としてもよい。

また、上述した実施の形態では、１つの解析モデル１中に含まれる複数のフィラーモデル１１の周囲に発生した破断座標Ｘ_１～Ｘ_３の数値解析の解析結果を集約して表示する例について説明したが、複数の解析用モデル１（例えば、１０個の解析用モデル）を用いて別途演算した解析結果を１つの解析用モデル１の解析結果に投影したアンサンブル結果を可視化してもよい。これにより、多数のフィラーモデル１１の周囲で発生する多数の破断座標の解析結果を集約して表示できるので、効率良く演算結果を解析することが可能となる。

なお、上述した実施の形態においては、図１７に示すように、ポリマーをモデル化したポリマーモデル２１及びフィラーをモデル化したフィラーモデル１１を含む複合材料の解析用モデルを作成するステップＳＴ２１と、ポリマーモデル２１を架橋させるステップＳＴ２２と、解析対象となるフィラーモデル１１に属する少なくとも一対の粒子の結合鎖２１ｂの粒子間距離に閾値を設定するステップＳＴ２３と、ポリマー粒子２１ａの質量を変化させるステップＳＴ２４と、粒子間距離が閾値以上の場合に、粒子間結合の破断処理をして解析用モデル１の数値解析を実行するステップＳＴ２５とを含むようにしてもよい。これにより、複合材料の解析方法は、架橋反応を介してポリマーモデル２１を予め架橋した状態で解析用モデルの数値解析をできるので、ポリマーモデル２１の架橋が第１結合鎖２１ｂの破断に及ぼす影響をより正確に解析可能となる。

次に、本実施の形態に係る複合材料の解析方法、及び複合材料の解析用コンピュータプログラムについてより詳細に説明する。図１８は、本実施の形態に係る複合材料の解析方法及び複合材料の解析方法を実行する解析装置の機能ブロック図である。

図１８に示すように、本実施の形態に係る複合材料の解析方法は、処理部５２と記憶部５４とを含むコンピュータである解析装置５０が実現する。この解析装置５０は、入力手段５３を備えた入出力装置５１と電気的に接続されている。入力手段５３は、複合材料の解析用モデルの作成対象であるポリマー及びフィラーの各種物性値、ポリマー及びフィラーを含有する複合材料を用いた伸張試験結果の実測結果、及び解析における境界条件などを処理部５２又は記憶部５４へ入力する。入力手段５３としては、例えば、キーボード、マウスなどの入力デバイスが用いられる。

処理部５２は、例えば、中央演算装置（ＣＰＵ：Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）及びメモリを含む。処理部５２は、各種処理を実行する際にコンピュータプログラムを記憶部５４から読み込んでメモリに展開する。メモリに展開されたコンピュータプログラムは、各種処理を実行する。例えば、処理部５２は、記憶部５４から予め記憶された各種処理に係るデータを必要に応じて適宜メモリ上の自身に割り当てられた領域に展開し、展開したデータに基づいて複合材料の解析用モデルの作成及び複合材料の解析用モデルを用いた複合材料の解析に関する各種処理を実行する。

処理部５２は、モデル作成部５２ａと、条件設定部５２ｂと、解析部５２ｃとを含む。モデル作成部５２ａは、予め記憶部５４に記憶されたデータに基づき、分子動力学法により複合材料の解析用モデル１を作成する際のフィラー及びポリマーなどの複合材料の粒子数、分子数、分子量、分子鎖長、分子鎖数、分岐、形状、大きさ、反応時間、反応条件及び作成する解析用モデルに含まれる分子数である目標分子数などの構成要素の配置、設定及び計算ステップ数などの粗視化モデルの設定を行う。また、モデル作成部５２ａは、フィラー粒子１１ａ間、ポリマー粒子２１ａ間及びフィラー・ポリマー粒子の水素結合、分子間力などの相互作用などの各種計算パラメーターの初期条件の設定を行う。また、モデル作成部５２ａは、必要に応じてポリマーモデル２１の架橋による架橋結合の作成などの架橋解析などを作成してもよい。

フィラー粒子１１ａ間の相互作用及びポリマー粒子２１ａ間の相互作用を調整する計算パラメーターとしては、下記式（４）で表されるレナード・ジョーンズポテンシャルのσ、εを用い、これらが調整される。ポテンシャルを計算する上限距離（カットオフ距離）を大きくすることで、遠距離まで働いた引力、斥力を調整できる。なお、フィラー粒子１１ａ間の相互作用及びポリマー粒子２１ａ間の相互作用が一定値になるまで順次、フィラー粒子１１ａ間の相互作用及びポリマー粒子２１ａ間の相互作用のパラメーターを小さくすることが好ましい。レナード・ジョーンズポテンシャルのσ、εを大きな値から徐々に本来の値に近づけることにより、分子を不自然な状態に導かない穏やかな速度で粒子の接近を行うことができる。また、カットオフ距離も徐々に小さくすることにより、適正な範囲で引力、斥力を調整できる。

条件設定部５２ｂは、変温解析及び変圧解析などの数値解析、伸張解析、せん断解析などの変形解析及び緩和解析などの運動解析などの各種数値解析条件を設定する。

解析部５２ｃは、条件設定部５２ｂによって設定された解析条件に基づいて解析用モデル１の各種数値解析を実行する。また、解析部５２ｃは、モデル作成部５２ａによって作成された複合材料の解析用モデル１を用いて分子動力学法による数値解析を実行して物理量を取得する。ここでは、解析部５２ｃは、数値解析として、伸張解析、せん断解析などの変形解析及び緩和解析などの運動解析などを実行する。また、解析部５２ｃは、数値解析の結果得られた変位などの値又は得られた値に所定の演算処理を実行した歪みなどの物理量を取得する。

記憶部５４は、ハードディスク装置、光磁気ディスク装置、フラッシュメモリ及びＣＤ－ＲＯＭなどの読み出しのみが可能な記録媒体である不揮発性のメモリ、並びに、ＲＡＭ（Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ）のような読み出し及び書き込みが可能な記録媒体である揮発性のメモリが適宜組み合わせられる。

記憶部５４には、入力手段５３を介して解析対象となる複合材料の解析用モデルを作成するためのデータであるゴムカーボンブラック、シリカ、及びアルミナなどのフィラーのデータ、ゴム、樹脂、及びエラストマーなどのポリマーのデータ、予め設定した物理量履歴である応力歪み曲線及び本実施の形態に係る複合材料の解析方法、複合材料の解析方法を実現するためのコンピュータプログラムなどが格納されている。このコンピュータプログラムは、コンピュータ又はコンピュータシステムに既に記録されているコンピュータプログラムとの組み合わせによって、本実施の形態に係る複合材料の解析方法を実現できるものであってもよい。ここでいう「コンピュータシステム」とは、ＯＳ（Ｏｐｅｒａｔｉｎｇ Ｓｙｓｔｅｍ）及び周辺機器などのハードウェアを含むものとする。

表示手段５５は、例えば、液晶表示装置等の表示用デバイスである。なお、記憶部５４は、データベースサーバなどの他の装置内にあってもよい。例えば、解析装置５０は、入出力装置５１を備えた端末装置から通信により処理部５２及び記憶部５４にアクセスするものであってもよい。

次に、再び図１を参照して、本実施の形態に係る複合材料の解析方法についてより詳細に説明する。

まず、図１に示すように、モデル作成部５２ａが、所定のモデル作成領域Ａ内にポリマー粒子２１ａ及び結合鎖２１ｂを含む未架橋の複数のポリマーモデル２１を作成すると共にフィラー粒子１１ａを含む複数のフィラーモデル１１を含む解析用モデル１を作成する（ステップＳＴ１１）。未架橋のポリマーモデル２１は、図２に示したように、複数のポリマー粒子２１ａが結合鎖２１ｂによって連結されてなるものである。ここでは、モデル作成部５２ａは、作成したフィラーモデル１１中に未架橋のポリマーモデル２１を配置する。次に、モデル作成部５２ａは、初期条件の設定の後、平衡化計算を行う。平衡化計算では、所定の温度、密度及び圧力で、初期設定後の各種構成要素が平衡状態に到達する所定の時間、分子動力学計算を行う。そして、モデル作成部５２ａは、初期条件の設定及び平衡化の計算処理後に、計算領域内に設定した複合材料の解析用モデルを作成するモデル作成領域Ａ内に、ポリマーモデル２１及びフィラーモデル１１を含む解析用モデル１を配置する。また、モデル作成部５２ａは、必要に応じてポリマーにフィラーとの親和性を高める水酸基、カルボニル基、及び原子団の官能基などの変性剤を配合してもよい。また、モデル作成部５２ａは、作成したポリマーモデル２１に架橋解析により架橋結合２１ｅを導入してもよい。なお、モデル作成部５２ａは、必要に応じて、複合材料の解析用モデル１に、分子間力及び水素結合などの引力及び斥力などの化学的な相互作用、及び共有結合などの物理的な相互作用を設定してもよい。

次に、条件設定部５２ｂが、モデル作成部５２ａで作成した複合材料の解析用モデル１を用いた分子動力学法による架橋解析、数値解析及び運動解析（シミュレーション）を実行するための各種条件を設定する。条件設定部５２ｂは、入力手段５３からの入力及び記憶部５４に記憶されている情報に基づいて各種条件を設定する。各種条件としては、解析を実行するフィラーモデル１１の位置及び数、フィラー原子、フィラー原子団、フィラー粒子１１ａ及びフィラー粒子群の位置及び数、フィラー粒子番号、ポリマーの分子鎖の位置及び数、ポリマー原子、ポリマー原子団、ポリマー粒子２１ａ及びポリマー粒子群の位置及び数、ポリマー粒子番号、結合鎖２１ｂの位置及び数、結合鎖２１ｂの番号、予め設定した物理量履歴である応力歪み曲線及び条件を変更しない固定値などが含まれる。

次に、解析用モデル１に相互作用を設定して変温解析、変圧解析などの各種数値解析を実行する。解析部５２ｃは、必要に応じて、例えば、フィラー粒子１１ａ間、ポリマー粒子２１ａ間、フィラー粒子１１ａとポリマー粒子２１ａとの間の相互作用及びフィラー粒子１１ａとポリマー粒子２１ａとが結合鎖で結合した状態の相互作用などを設定する。

解析部５２ｃは、解析対象とする結合鎖２１ｂによって粒子間結合された一対のポリマーモデル２１の粒子間距離に第１閾値を設定する（ステップＳＴ１２）。ここでは、解析部５２ｃは、閾値としては、粒子間距離に特定の第１閾値を設定してもよく、所定の数値範囲を有する閾値範囲を第１閾値として設定してもよい。また、解析部５２ｃは、ポリマーモデル２１のポリマー粒子２１ａを質量変換粒子２２ａに変換して質量を変化させる（ステップＳＴ１３）。ここでは、解析部５２ｃは、質量変換粒子２２ａの質量をポリマー粒子２１ａに対して増大させてもよく、減少させてもよい。また、解析部５２ｃは、ポリマーモデル２１の全てのポリマー粒子２１ａを質量変換粒子２２ａに変換してもよく、一部のポリマー粒子２１ａを質量変換粒子２２ａに変換してもよい。また、解析部５２ｃは、第１閾値より小さい第２閾値に基づいてポリマー粒子２１ａを質量変換粒子２２ａに変換してもよい。また、解析部５２ｃは、破断処理後の経過時間に応じて質量変換粒子２２ａの質量を変化させてもよい。

次に、解析部５２ｃは、複合材料の解析用モデル１を用いた分子動力学法による緩和解析、伸張解析、せん断解析及び変形解析などの運動解析などの各種数値解析を実行する（ステップＳＴ１４）。ここでは、解析部５２ｃは、解析対象となる一対のポリマーモデル２１ａの結合鎖２１ｂが第１閾値以上となった際に、結合鎖２１ｂを破断処理して数値解析を実行する。また、解析部５２ｃは、解析対象となる一対のポリマーモデル２１ａの粒子間距離が第１閾値未満の場合に第１破断結合演算用関数を使用し、当該粒子間距離が第１閾値以上となった際に、結合鎖２１ｂの結合エネルギー及び結合力が低下する第２破断結合演算用関数を用いて解析用モデル１の各種数値解析を実行してもよい。さらに、解析部５２ｃは、結合鎖２１ｃの作成後、質量変換粒子２２ａをポリマー粒子２１ａに戻してもよい。

また、解析部５２ｃは、数値解析による運動解析の結果得られる運動変位及び公称応力又は運動変位を演算して得られる公称歪みなどの各種物理量を取得する。このような数値解析及び運動解析により、解析時間毎に変化する解析用モデル全体のポリマー分子の結合長及びポリマー粒子速度、架橋点間と自由末端の速度又は結合長、配向などの物理量などのセグメントの状態変化を表す数値と歪みとの関係、解析時間毎に変化するポリマー分子の結合長及びポリマー粒子速度などのセグメントの状態変化を表す数値と圧力又は解析時間との関係、解析時間毎に変化するポリマー分子の結合長及びポリマー粒子速度などのセグメントの状態変化を表す数値と温度又は解析時間との関係などを評価できるので、ポリマー分子の局所的な分子状態変化のより詳細な解析が可能となる。

また、解析部５２ｃは、数値解析によって得られたポリマーモデル２１の破断座標を特定し、特定した破断座標を評価する。ここでは、解析部５２ｃは、破断した結合鎖２１ｂを可視化して評価してもよく、破断座標を可視化して評価してもよい。さらに、解析部５２ｃは、複数のフィラーモデル１１の周囲に発生した破断座標を集約して評価してもよく、複数のフィラーモデル１１の周囲に発生した破断座標を１つの代表フィラーモデル１１Ｅに集約して評価してもよい。また、解析部５２ｃは、複数の解析用モデル１を用いて別途解析した解析結果を集約して評価してもよい。次に、解析部５２ｃは、解析した複合材料の解析結果を記憶部５４に格納する。

（実施例）
次に、本発明の効果を明確にするために行った実施例について説明する。なお、本発明は、以下の実施例によって何ら限定されるものではない。

本発明者らは、上記実施の形態に係る複合材料の解析方法を用いてポリマー粒子２１ａの質量を増大させた解析用モデル１を数値解析した場合（実施例１）と、実際の複合材料を用いて応力歪曲線を測定した場合（参考例１）と、ポリマー粒子２１ａの質量を変化させずに解析用モデル１を数値解析した場合（比較例１）と、結合鎖２１ｂの破断を考慮せずに解析用モデル１を数値解析した場合（比較例２）とを対比して評価した。

図１９は、本実施例に係る応力歪曲線を示す図である。図１９に示すように、上記実施の形態に係る複合材料の解析方法を用いた場合（実施例１：実線Ｌ１１参照）には、実際の複合材料を用いて応力歪曲線を測定した場合（参考例１：一点鎖線Ｌ１２参照）と同様に、応力と共に歪みが一定値まで上昇した後、応力は略一定に維持される結果となった。これに対して、ポリマー粒子２１ａの質量を変化させずに解析用モデル１を数値解析した場合（比較例１：二点鎖線Ｌ１３参照）には、応力が上昇して歪みが一定値まで上昇した後、応力が著しく減少する結果が得られた。また、結合鎖２１ｂの破断を考慮しない場合（比較例２：破線Ｌ１４参照）は、応力の上昇に伴い歪みが連続的に上昇した。この結果は、実施例１では、ポリマー粒子２１ａの質量の増大により破断したポリマーモデル２１の分子鎖の収縮速度を低減できたために、参考例１と同様の結果が得られたものと考えられる。

このように、上記実施例によれば、ポリマー粒子２１ａの質量を変化させると共に粒子間距離が第１閾値以上となった場合に、結合鎖２１ｂを破断処理して解析用モデル１を数値解析するので、破断に伴う力学応答の再現性に優れた複合材料の解析方法を実現できることが分かる。

１ 解析用モデル
１１，１１Ａ，１１Ｂ，１１Ｃ，１１Ｄ フィラーモデル
１１Ｅ 代表フィラーモデル
１１ａ フィラー粒子
２１，２１Ａ，２１Ｂ，２１Ｃ ポリマーモデル
２１ａ ポリマー粒子
２１ｂ，２１ｃ 結合鎖
２１ｂｘ 破断結合鎖
２２ａ 質量変換粒子
５０ 解析装置
５１ 入出力装置
５２ 処理部
５２ａ モデル作成部
５２ｂ 条件設定部
５２ｃ 解析部
５３ 入力手段
５４ 記憶部
５５ 表示手段
Ａ モデル作成領域
Ｘ_１，Ｘ_２，Ｘ_３ 破断座標

Claims (13)

1. コンピュータを用いた分子動力学法による特定物質の解析方法であって、
   ２以上の物質を含むゴム材料である特定物質をモデル化した特定物質モデルを含む特定物質の解析用モデルを作成するステップと、
   解析対象となる特定物質モデルに属し、粒子間結合で結合された少なくとも一対の粒子の粒子間距離に第１閾値を設定するステップと、
   前記特定物質の粒子の質量を変化させるステップと、
   前記粒子間距離が前記第１閾値以上の場合に、前記粒子間結合を破断処理して前記解析用モデルの変形解析又は運動解析である数値解析を実行するステップと、
   を含むことを特徴とする、特定物質の解析方法。
2. さらに、前記特定物質モデルを架橋させるステップを含む、請求項１に記載の特定物質の解析方法。
3. 前記破断処理される前記粒子間結合を含む分子鎖に属する前記特定物質の粒子の質量を変化させる、請求項１又は請求項２に記載の特定物質の解析方法。
4. 前記破断処理される前記分子鎖の末端部の前記特定物質の粒子の質量を変化させる、請求項３に記載の特定物質の解析方法。
5. 前記破断処理と共に前記特定物質の質量を変化させる、請求項１から請求項４のいずれか１項に記載の特定物質の解析方法。
6. 前記粒子間距離が前記第１閾値より小さい第２閾値以上となった際に、前記特定物質の粒子の質量を変化させる、請求項１から請求項５のいずれか１項に記載の特定物質の解析方法。
7. 前記破断処理された分子鎖の形状の変化が集束してから、質量を変化させた前記特定物質の粒子の質量の変化量を小さくする、請求項１から請求項６のいずれか１項に記載の特定物質の解析方法。
8. 前記破断処理後の経過時間に応じて前記特定物質の粒子の質量の変化量を変更する、請求項１から請求項７のいずれか１項に記載の特定物質の解析方法。
9. 質量分布を設けて前記特定物質の粒子の質量を変化させる、請求項１から請求項８のいずれか１項に記載の特定物質の解析方法。
10. 前記質量の変化が、質量の増大である、請求項１から請求項９のいずれか１項に記載の特定物質の解析方法。
11. 前記破断処理が、破断結合演算用関数を用いた疑似破断処理である、請求項１から請求項１０のいずれか１項に記載の特定物質の解析方法。
12. 前記特定物質が、ポリマー及びフィラーを含む、請求項１から請求項１１のいずれか１項に記載の特定物質の解析方法。
13. 請求項１から請求項１２のいずれか１項に記載の特定物質の解析方法をコンピュータに実行させることを特徴とする、特定物質の解析用コンピュータプログラム。

Images