JP7159809B2 - Rubber material simulation method and rubber material manufacturing method - Google Patents
Rubber material simulation method and rubber material manufacturing method Download PDFInfo
- Publication number
- JP7159809B2 JP7159809B2 JP2018219327A JP2018219327A JP7159809B2 JP 7159809 B2 JP7159809 B2 JP 7159809B2 JP 2018219327 A JP2018219327 A JP 2018219327A JP 2018219327 A JP2018219327 A JP 2018219327A JP 7159809 B2 JP7159809 B2 JP 7159809B2
- Authority
- JP
- Japan
- Prior art keywords
- rubber material
- model
- interaction potential
- polymer chain
- coarse
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Description
本発明は、物性値が異なる複数の高分子成分を含むゴム材料の破壊現象を、コンピュータを用いてシミュレーションするための方法などに関する。 TECHNICAL FIELD The present invention relates to a method for simulating, using a computer, a fracture phenomenon of a rubber material containing a plurality of polymer components having different physical property values.
下記非特許文献１は、ホモポリマーＡ、Ｂを含むポリマーブレンドについて、それらの界面での破壊特性を評価するためのシミュレーション方法を提案している。下記非特許文献１の方法では、先ず、ホモポリマーＡを、ビーズＡを用いてモデル化したビーズスプリングモデルＡが設定される。次に、ホモポリマーＢを、ビーズＢを用いてモデル化したビーズスプリングモデルＢが設定される。ビーズスプリングモデルＡのビーズＡ、Ａ間、及び、ビーズスプリングモデルＢのビーズＢ、Ｂ間には、相互作用ポテンシャルがそれぞれ定義されている。
The following
次に、下記非特許文献１の方法では、予め定められた単位セルに、複数のビーズスプリングモデルＡ、Ｂが配置される。単位セルには、ビーズスプリングモデルＡとビーズスプリングモデルＢとがそれぞれ別の領域に配置され、それらの領域間で界面が形成されている。そして、下記非特許文献１の方法では、単位セルの変形が計算され、ポリマーブレンドの破壊エネルギー等が計算される。
Next, in the method of
上記非特許文献１の方法では、ポリマーモデルについて引張試験を行った場合、主として、ビーズスプリングモデルＡが配置された領域と、複数のビーズスプリングモデルＢが配置された領域との界面で破壊が生じる現象が計算される。しかしながら、実際のポリマーブレンドを変形させると、ホモポリマーＡ、Ｂのうち、耐変形性を示す物性値（例えば、弾性率）が相対的に小さいポリマー内で破壊が観察される。したがって、上記非特許文献１の方法では、ポリマーブレンドの破壊特性を精度よく評価するには、さらなる改善の余地があった。
In the method of Non-Patent
発明者らは、鋭意研究を重ねた結果、ホモポリマーＡ、Ｂの物性値の差異を考慮して、ビーズＡ、Ａ間の相互作用ポテンシャル、及び、ビーズＢ、Ｂ間の相互作用ポテンシャルを互いに異なるものとして定義することが、ポリマーブレンドの破壊特性を評価する上で有効であることを知見した。 As a result of intensive research, the inventors have found that the interaction potential between beads A and A and the interaction potential between beads B and B are mutually calculated in consideration of the difference in physical property values of homopolymers A and B. We have found that defining as different is useful in assessing the fracture properties of polymer blends.
本発明は、以上のような実状に鑑み案出されたもので、ゴム材料の破壊特性をより精度良く評価することが可能なゴム材料のシミュレーション方法などを提供することを主たる目的としている。 The present invention has been devised in view of the actual situation as described above, and a main object thereof is to provide a rubber material simulation method and the like capable of more accurately evaluating the fracture characteristics of the rubber material.
本発明は、予め定められた物性値について差異を有する第１高分子成分と第２高分子成分とを少なくとも含むゴム材料の破壊現象を、コンピュータを用いてシミュレーションするための方法であって、前記第１高分子成分の分子鎖構造を複数の第１粗視化粒子を用いて表現した第１高分子鎖モデルを前記コンピュータに入力する工程と、前記第２高分子成分の分子鎖構造を複数の第２粗視化粒子を用いて表現した第２高分子鎖モデルを前記コンピュータに入力する工程と、前記第１高分子鎖モデルと前記第２高分子鎖モデルとの間に界面が形成されるように、前記第１高分子鎖モデル及び前記第２高分子鎖モデルを、前記ゴム材料の一部に対応する仮想空間であるセルに配置したゴム材料モデルを定義する工程と、前記第１粗視化粒子間の相互作用を定義するための第１相互作用ポテンシャルを定義する工程と、前記第２粗視化粒子間の相互作用を定義するための第２相互作用ポテンシャルを、前記物性値の差異に基づいて、前記第１相互作用ポテンシャルとは異なるものとして定義する工程と、前記コンピュータが、前記ゴム材料モデルの変形を計算して、前記第１相互作用ポテンシャル、及び、前記第２相互作用ポテンシャルに基づく前記界面での破壊現象を計算する工程とを含むことを特徴とする。 The present invention provides a method for simulating, using a computer, the destruction phenomenon of a rubber material containing at least a first polymer component and a second polymer component having predetermined physical property values different from each other, the method comprising: inputting into the computer a first polymer chain model representing the molecular chain structure of the first polymer component using a plurality of first coarse-grained particles; A step of inputting a second polymer chain model expressed using the second coarse-grained particles into the computer, and forming an interface between the first polymer chain model and the second polymer chain model defining a rubber material model in which the first polymer chain model and the second polymer chain model are arranged in a cell that is a virtual space corresponding to a part of the rubber material; a step of defining a first interaction potential for defining interaction between coarse-grained particles; and calculating a deformation of the rubber material model to generate the first interaction potential and the second interaction potential based on the difference between the first interaction potential and the second interaction potential. and calculating a fracture phenomenon at the interface based on the action potential.
本発明に係る前記ゴム材料のシミュレーション方法において、前記第１相互作用ポテンシャル及び前記第２相互作用ポテンシャルは、相互作用の強さの定義するための強度パラメータを含み、前記第２相互作用ポテンシャルを定義する工程は、前記第２相互作用ポテンシャルに、前記第１相互作用ポテンシャルとは異なる前記強度パラメータを定義する工程を含んでもよい。 In the rubber material simulation method according to the present invention, the first interaction potential and the second interaction potential include a strength parameter for defining the strength of interaction, and the second interaction potential is defined The step of determining may comprise defining a different strength parameter for the second interaction potential than for the first interaction potential.
本発明に係る前記ゴム材料のシミュレーション方法において、前記第１相互作用ポテンシャル及び前記第２相互作用ポテンシャルは、以下の式で定義されてもよい。
ｒij：粒子間の距離
ｒij,cut：カットオフ距離
εij：強度パラメータ
σij：粒子の直径に相当
In the rubber material simulation method according to the present invention, the first interaction potential and the second interaction potential may be defined by the following equations.
r ij : Distance between particles r ij , cut : Cutoff distance ε ij : Strength parameter σ ij : Corresponding to particle diameter
本発明に係る前記ゴム材料のシミュレーション方法において、前記強度パラメータを定義する工程は、予め定められた複数種類の前記強度パラメータεij毎に、前記物性値を計算する工程と、前記物性値の差異に基づいて、前記複数種類の前記強度パラメータεijから、前記第２相互作用ポテンシャルの前記強度パラメータεijを選択する工程とを含んでもよい。 In the rubber material simulation method according to the present invention, the step of defining the strength parameter includes the step of calculating the physical property value for each of a plurality of predetermined strength parameters ε ij , and calculating the difference between the physical property values. and selecting the strength parameter ε ij of the second interaction potential from the plurality of types of strength parameters ε ij based on.
本発明に係る前記ゴム材料のシミュレーション方法において、前記物性値は、弾性率を含んでもよい。 In the rubber material simulation method according to the present invention, the physical property value may include an elastic modulus.
本発明の製造方法は、請求項１ないし５のいずれかに記載の前記シミュレーション方法において、前記性能が良好であると評価された前記ゴム材料モデルに基づいて、前記ゴム材料を製造する工程を含むことを特徴とする。
The production method of the present invention includes the step of producing the rubber material based on the rubber material model evaluated as having good performance in the simulation method according to any one of
本発明のゴム材料のシミュレーション方法は、第１高分子成分の分子鎖構造を複数の第１粗視化粒子を用いて表現した第１高分子鎖モデルをコンピュータに入力する工程と、第２高分子成分の分子鎖構造を複数の第２粗視化粒子を用いて表現した第２高分子鎖モデルを前記コンピュータに入力する工程とを含んでいる。さらに、本発明の前記シミュレーションは、前記第１高分子鎖モデルと前記第２高分子鎖モデルとの間に界面が形成されるように、前記第１高分子鎖モデル及び前記第２高分子鎖モデルを、前記ゴム材料の一部に対応する仮想空間であるセルに配置したゴム材料モデルを定義する工程を含んでいる。 The rubber material simulation method of the present invention comprises the steps of: inputting into a computer a first polymer chain model representing the molecular chain structure of a first polymer component using a plurality of first coarse-grained particles; inputting into the computer a second polymer chain model representing the molecular chain structure of the molecular component using a plurality of second coarse-grained particles. Furthermore, the simulation of the present invention includes the first polymer chain model and the second polymer chain model such that an interface is formed between the first polymer chain model and the second polymer chain model. It includes a step of defining a rubber material model in which the model is arranged in a cell, which is a virtual space corresponding to a portion of the rubber material.
本発明の前記シミュレーション方法は、前記第１粗視化粒子間の相互作用を定義するための第１相互作用ポテンシャルを定義する工程と、前記第２粗視化粒子間の相互作用を定義するための第２相互作用ポテンシャルを、前記第１高分子成分と前記第２高分子成分との物性値の差異に基づいて、前記第１相互作用ポテンシャルとは異なるものとして定義する工程とを含んでいる。これにより、本発明の前記シミュレーション方法は、前記物性値の差異を考慮することが可能な前記ゴム材料モデルを作成することができる。 The simulation method of the present invention comprises a step of defining a first interaction potential for defining interaction between the first coarse-grained particles, and a step of defining an interaction between the second coarse-grained particles. defining a second interaction potential of as being different from the first interaction potential based on the difference in physical property values between the first polymer component and the second polymer component . Thereby, the simulation method of the present invention can create the rubber material model that can consider the difference in the physical property values.
本発明の前記シミュレーション方法は、前記コンピュータが、前記ゴム材料モデルの変形を計算して、前記第１相互作用ポテンシャル、及び、前記第２相互作用ポテンシャルに基づく前記界面での破壊現象を計算する工程を含んでいる。これにより、本発明の前記シミュレーション方法は、前記第１高分子成分と前記第２高分子成分とを含むゴム材料の破壊特性を、精度良く評価することができる。 In the simulation method of the present invention, the computer calculates deformation of the rubber material model to calculate a fracture phenomenon at the interface based on the first interaction potential and the second interaction potential. contains. Thereby, the simulation method of the present invention can accurately evaluate the fracture characteristics of the rubber material containing the first polymer component and the second polymer component.
以下、本発明の実施の一形態が図面に基づき説明される。
本実施形態のゴム材料のシミュレーション方法（以下、単に「シミュレーション方法」ということがある。）では、ゴム材料の破壊現象が、コンピュータを用いてシミュレーション（計算）される。
An embodiment of the present invention will be described below with reference to the drawings.
In the rubber material simulation method of the present embodiment (hereinafter sometimes simply referred to as "simulation method"), the destruction phenomenon of the rubber material is simulated (calculated) using a computer.
図１は、ゴム材料のシミュレーション方法を実行するためのコンピュータの一例を示す斜視図である。コンピュータ１は、本体１ａ、キーボード１ｂ、マウス１ｃ、及び、ディスプレイ装置１ｄを含んでいる。本体１ａには、例えば、演算処理装置（ＣＰＵ）、ＲＯＭ、作業用メモリ、磁気ディスクなどの記憶装置、及び、ディスクドライブ装置１ａ１、１ａ２が設けられている。記憶装置には、本実施形態のシミュレーション方法を実行するためのソフトウェア等が予め記憶されている。
FIG. 1 is a perspective view showing an example of a computer for executing a rubber material simulation method. A
本実施形態のゴム材料は、第１高分子成分と第２高分子成分とを少なくとも含んでいる。第１高分子成分及び第２高分子成分としては、例えば、天然ゴム、ブタジエンゴム、ポリイソプレンゴム、又は、スチレンブタジエンゴム等が例示されるが、これらに限定されるわけではない。 The rubber material of this embodiment contains at least a first polymer component and a second polymer component. Examples of the first polymer component and the second polymer component include, but are not limited to, natural rubber, butadiene rubber, polyisoprene rubber, styrene-butadiene rubber, and the like.
第１高分子成分及び第２高分子成分は、例えば、上記の高分子成分の中から、互いに異なる高分子成分がそれぞれ選択される。本実施形態では、予め定められた物性値について差異を有するように、第１高分子成分及び第２高分子成分がそれぞれ選択される。物性値については、特に限定されない。本実施形態の物性値には、耐変形性を示す弾性率（ヤング率）が含まれる。本実施形態において、第１高分子成分がブタジエンゴムであり、第２高分子成分がスチレンブタジエンゴムである場合が例示される。 For the first polymer component and the second polymer component, for example, polymer components different from each other are selected from among the polymer components described above. In the present embodiment, the first polymer component and the second polymer component are each selected so that they have different predetermined physical property values. Physical properties are not particularly limited. The physical property values of this embodiment include an elastic modulus (Young's modulus) that indicates deformation resistance. In this embodiment, the case where the first polymer component is butadiene rubber and the second polymer component is styrene-butadiene rubber is exemplified.
本実施形態のゴム材料には、架橋剤が含まれている。架橋剤としては、硫黄である場合が例示される。 The rubber material of this embodiment contains a cross-linking agent. A cross-linking agent is exemplified by sulfur.
図２は、ゴム材料のシミュレーション方法及び製造方法の処理手順の一例を示すフローチャートである。本実施形態のシミュレーション方法では、先ず、コンピュータ１に、ゴム材料モデルが入力される（ゴム材料モデル入力工程Ｓ１）。図３は、ゴム材料モデル入力工程Ｓ１の処理手順の一例を示すフローチャートである。図４は、第１高分子鎖モデル３、第２高分子鎖モデル４、及び、架橋剤粒子モデル５が配置されたセル１５の一例を示す概念図である。図５は、第１高分子鎖モデル３及び架橋剤粒子モデル５の一例を示す概念図である。図４は、一つの第１高分子鎖モデル３、一つの第２高分子鎖モデル４、及び、２つの架橋剤粒子モデル５を代表して示している。
FIG. 2 is a flow chart showing an example of the procedure of the rubber material simulation method and manufacturing method. In the simulation method of this embodiment, first, a rubber material model is input to the computer 1 (rubber material model input step S1). FIG. 3 is a flow chart showing an example of the procedure of the rubber material model input step S1. FIG. 4 is a conceptual diagram showing an example of the
本実施形態のゴム材料モデル入力工程Ｓ１では、先ず、図４及び図５に示されるように、第１高分子成分（本実施形態では、ブタジエンゴム）の分子鎖構造を、複数の第１粗視化粒子７を用いて表現した第１高分子鎖モデル３が、コンピュータ１に入力される（工程Ｓ１１）。本実施形態の第１高分子鎖モデル３は、粗視化モデル（本実施形態では、Kremer-Grestモデル）として定義されている。本実施形態の第１高分子鎖モデル３は、複数の第１粗視化粒子７と、隣接する第１粗視化粒子７、７間を結合する第１結合鎖モデル８とを含んで構成されている。
In the rubber material model input step S1 of the present embodiment, first, as shown in FIGS. The first
第１粗視化粒子７は、第１高分子成分（本実施形態では、ブタジエンゴム）の分子鎖のモノマー（図示省略）又はモノマーの一部分をなす構造単位を置換したものである。これにより、第１高分子鎖モデル３には、複数個（例えば、１０～５０００個）の第１粗視化粒子７が設定される。第１粗視化粒子７は、後述の分子動力学計算において、運動方程式の質点として取り扱われる。このため、第１粗視化粒子７には、例えば、質量、体積、粒子径Ｄ１、又は、電荷などのパラメータが定義される。
The first coarse-
本実施形態の第１結合鎖モデル８は、第１粗視化粒子７、７間に、伸びきり長が設定された第３相互作用ポテンシャルＰ３によって定義される。本実施形態の第３相互作用ポテンシャルＰ３は、下記式（２）の非調和ポテンシャルＵch（ｒ）によって定義される。
The first
ｒ：粒子モデル間の距離
ｋ：粒子モデル間のばね定数
Ｒ0：伸びきり長
なお、距離ｒ、及び、伸びきり長Ｒ0は、第１粗視化粒子７の中心７ｃ（図５に示す）の座標に基づいて設定される。
r: distance between particle models k : spring constant between particle models R 0 : stretched length ) is set based on the coordinates of
非調和ポテンシャルＵch（ｒ）の各定数及び各変数については、例えば、論文１（ 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）に基づいて適宜設定される。これにより、工程Ｓ１１では、第１粗視化粒子７が伸縮自在に拘束された直鎖状の第１高分子鎖モデル３を定義することができる。第１高分子鎖モデル３は、コンピュータ１に記憶される。
For each constant and each variable of the anharmonic potential U ch (r), see, for example, Article 1 (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). As a result, in step S11, the linear first
次に、本実施形態のゴム材料モデル入力工程Ｓ１では、図４に示されるように、第２高分子成分（本実施形態では、スチレンブタジエンゴム）の分子鎖構造を、複数の第２粗視化粒子１１を用いて表現した第２高分子鎖モデル４が、コンピュータ１に入力される（工程Ｓ１２）。本実施形態の第２高分子鎖モデル４は、第１高分子鎖モデル３と同様に、粗視化モデルとして定義されている。本実施形態の第２高分子鎖モデル４は、複数の第２粗視化粒子１１と、隣接する第２粗視化粒子１１、１１間を結合する第２結合鎖モデル１２とを含んで構成されている。
Next, in the rubber material model input step S1 of the present embodiment, as shown in FIG. The second
第２粗視化粒子１１は、第２高分子成分（本実施形態では、スチレンブタジエンゴム）の分子鎖のモノマー（図示省略）又はモノマーの一部分をなす構造単位を置換したものである。これにより、第２高分子鎖モデル４には、複数個（例えば、１０～５０００個）の第２粗視化粒子１１が設定される。第２粗視化粒子１１は、後述の分子動力学計算において、運動方程式の質点として取り扱われる。このため、第２粗視化粒子１１には、第１粗視化粒子７と同様のパラメータが定義される。
The second coarse-
第２結合鎖モデル１２は、第２粗視化粒子１１、１１間に、伸びきり長が設定された第４相互作用ポテンシャルＰ４によって定義される。本実施形態の第４相互作用ポテンシャルＰ４は、上記式（２）の非調和ポテンシャルＵch（ｒ）によって定義される。非調和ポテンシャルＵch（ｒ）の各定数及び各変数については、第１高分子鎖モデル３の各定数及び各変数と同様に設定される。これにより、工程Ｓ１２では、第２粗視化粒子１１が伸縮自在に拘束された直鎖状の第２高分子鎖モデル４を定義することができる。第２高分子鎖モデル４は、コンピュータ１に記憶される。
The second
次に、本実施形態のゴム材料モデル入力工程Ｓ１は、図４及び図５に示されるように、架橋剤をモデル化した架橋剤粒子モデル５が、コンピュータ１に入力される（工程Ｓ１３）。本実施形態の架橋剤粒子モデル５は、一つの独立した粒子モデルとして設定される。架橋剤粒子モデル５は、分子動力学計算において、運動方程式の質点として取り扱われる。このため、架橋剤粒子モデル５には、例えば、質量、体積、粒子径Ｄ２、又は、電荷などのパラメータが定義される。
Next, in the rubber material model input step S1 of the present embodiment, as shown in FIGS. 4 and 5, a cross-linking
架橋剤粒子モデル５の粒子径Ｄ２については、適宜設定することができる。本実施形態の粒子径Ｄ２は、第１粗視化粒子７及び第２粗視化粒子１１の粒子径Ｄ１と同一に設定されている。架橋剤粒子モデル５は、コンピュータ１に記憶される。
The particle diameter D2 of the cross-linking
次に、本実施形態のゴム材料モデル入力工程Ｓ１では、図４に示されるように、セル１５が、コンピュータ１に入力される（工程Ｓ１４）。本実施形態のセル１５は、解析対象のゴム材料の一部に対応する仮想空間である。セル１５は、互いに向き合う三対の平面１５ａ、１５ｂを有する直方体として定義されている。各平面１５ａ、１５ｂには、周期境界条件が定義されている。このようなセル１５は、例えば、一方の平面１５ａから出て行った第１高分子鎖モデル３、第２高分子鎖モデル４、及び、架橋剤粒子モデル５の一部が、反対側の平面１５ｂから入ってくるように計算することができる。したがって、一方の平面１５ａと、反対側の平面１５ｂとが連続している（繋がっている）ものとして取り扱うことができる。
Next, in the rubber material model input step S1 of this embodiment, as shown in FIG. 4, the
セル１５の一辺の各長さＬ１ａ、Ｌ１ｂ及びＬ１ｃは、適宜設定することができる。本実施形態の長さＬ１ａ、Ｌ１ｂ及びＬ１ｃは、第１高分子鎖モデル３、及び、第２高分子鎖モデル４の拡がりを示す量である慣性半径（図示省略）の２倍以上が望ましい。これにより、セル１５は、分子動力学計算において、周期境界条件による自己のイメージとの衝突の発生を防ぐことができるため、第１高分子鎖モデル３、及び、第２高分子鎖モデル４の空間的拡がりを適切に計算することができる。また、セル１５の大きさは、例えば１気圧で安定な体積に設定される。これにより、セル１５は、解析対象のゴム材料の少なくとも一部の体積を定義することができる。セル１５は、コンピュータ１に記憶される。
Each length L1a, L1b, and L1c of one side of the
次に、本実施形態のゴム材料モデル入力工程Ｓ１では、コンピュータ１が、解析対象のゴム材料に含まれる成分（第１高分子成分、第２高分子成分、及び、架橋剤）の分布を求める（分布計算工程Ｓ１５）。図６は、分布計算工程Ｓ１５の処理手順の一例を示すフローチャートである。
Next, in the rubber material model input step S1 of the present embodiment, the
本実施形態の分布計算工程Ｓ１５では、先ず、第１高分子成分、第２高分子成分、及び、架橋剤のχパラメータが求められる（工程Ｓ６１）。χパラメータは、Flory-Huggins理論に基づくものであり、第１高分子成分、第２高分子成分、及び、架橋剤の相溶性を表すものである。χパラメータの値が小さいほど、良溶媒であることを示している。第１高分子成分と第２高分子成分とによって形成される界面は、相溶性が高い（即ち、χパラメータが小さい）ほど、互いに混じり合った混合状態になりやすい。逆に、相溶性が低い（χパラメータが大きい）ほど、第１高分子成分及び第２高分子成分が互いに混じり合い難くなり、相分離した界面が形成される。 In the distribution calculation step S15 of this embodiment, first, the χ parameters of the first polymer component, the second polymer component, and the cross-linking agent are obtained (step S61). The χ parameter is based on the Flory-Huggins theory and represents the compatibility of the first polymer component, the second polymer component, and the cross-linking agent. A smaller value of the χ parameter indicates a better solvent. At the interface formed by the first polymer component and the second polymer component, the higher the compatibility (that is, the smaller the χ parameter), the more likely they are to be in a mixed state. Conversely, the lower the compatibility (the larger the χ parameter), the more difficult it is for the first polymer component and the second polymer component to mix with each other, forming a phase-separated interface.
本実施形態の工程Ｓ６１では、ブタジエンゴム間の第１χパラメータ、スチレンブタジエンゴム間の第２χパラメータ、ブタジエンゴムとスチレンブタジエンゴムとの間の第３χパラメータが求められる。さらに、工程Ｓ６１では、架橋剤間の第４χパラメータ、ブタジエンゴムと架橋剤との間の第５χパラメータ、及び、スチレンブタジエンゴムと架橋剤との間の第６χパラメータが求められる。各χパラメータは、従来と同様に、（株）ＪＳＯＬ社製のソフトマテリアル総合シミュレーター（Ｊ－ＯＣＴＡ）に含まれるＳＵＳＨＩモデラを用いて、原子団寄与法、モンテカルロ法、分子動力学法、又は、ＱＳＰＲ法に基づいて計算することができる。各χパラメータは、コンピュータ１に記憶される。
In step S61 of the present embodiment, a first χ parameter between butadiene rubbers, a second χ parameter between styrene-butadiene rubbers, and a third χ parameter between butadiene rubbers and styrene-butadiene rubbers are determined. Further, in step S61, a fourth χ parameter between cross-linking agents, a fifth χ parameter between butadiene rubber and cross-linking agent, and a sixth χ parameter between styrene-butadiene rubber and cross-linking agent are determined. As in the past, each χ parameter is calculated using the SUSHI modeler included in JSOL Co., Ltd.'s soft material general simulator (J-OCTA), atomic group contribution method, Monte Carlo method, molecular dynamics method, or It can be calculated based on the QSPR method. Each χ parameter is stored in
次に、本実施形態の分布計算工程Ｓ１５では、解析対象のゴム材料に含まれる成分（第１高分子成分、第２高分子成分、及び、架橋剤）の分布が求められる（工程Ｓ６２）。本実施形態では、第１高分子成分、第２高分子成分、及び、架橋剤の配置のエントロピーを考慮した静的なＳＣＦ法（ Self-Consistent Field Method ）を用いて、ゴム材料に含まれる第１高分子成分（ブタジエンゴム）の分布、第２高分子成分（スチレンブタジエンゴム）の分布、及び、架橋剤の分布がそれぞれ求められる。ＳＣＦ法は、例えば、（株）ＪＳＯＬ社製のソフトマテリアル総合シミュレーター（Ｊ－ＯＣＴＡ）に含まれるＳＵＳＨＩを用いて計算することができる。 Next, in the distribution calculation step S15 of this embodiment, the distribution of the components (the first polymer component, the second polymer component, and the cross-linking agent) contained in the rubber material to be analyzed is obtained (step S62). In the present embodiment, a static SCF method (Self-Consistent Field Method) that considers the entropy of the arrangement of the first polymer component, the second polymer component, and the cross-linking agent is used to The distribution of the first polymer component (butadiene rubber), the distribution of the second polymer component (styrene-butadiene rubber), and the distribution of the cross-linking agent are obtained. The SCF method can be calculated using, for example, SUSHI included in JSOL Co., Ltd. Soft Material General Simulator (J-OCTA).
ＳＣＦ法の入力値としては、第１高分子成分（ブタジエンゴム）の分子量、第２高分子成分（スチレンブタジエンゴム）の分子量、架橋剤の分子量、及び、χパラメータ（本実施形態では、第１χパラメータ～第６χパラメータ）が用いられる。これにより、工程Ｓ６２では、解析対象のゴム材料に含まれる第１高分子成分（ブタジエンゴム）の分布（濃度分布）、第２高分子成分（スチレンブタジエンゴム）の分布（濃度分布）、及び、架橋剤の分布（濃度分布）がそれぞれ求められる。図７は、第１高分子成分の分布、第２高分子成分の分布、及び、架橋剤の分布の一例を示すグラフである。このグラフでは、第１高分子成分、第２高分子成分及び架橋剤の体積分率と、それらのゴム材料での空間位置との関係を示している。ゴム材料に含まれる成分の分布は、コンピュータ１に入力される。
Input values for the SCF method include the molecular weight of the first polymer component (butadiene rubber), the molecular weight of the second polymer component (styrene-butadiene rubber), the molecular weight of the cross-linking agent, and the χ parameter (in this embodiment, the first χ parameter to sixth χ parameter) are used. As a result, in step S62, the distribution (concentration distribution) of the first polymer component (butadiene rubber), the distribution (concentration distribution) of the second polymer component (styrene-butadiene rubber), and The distribution (concentration distribution) of the cross-linking agent is determined respectively. FIG. 7 is a graph showing an example of the distribution of the first polymer component, the distribution of the second polymer component, and the distribution of the cross-linking agent. This graph shows the relationship between the volume fractions of the first polymeric component, the second polymeric component and the cross-linking agent and their spatial positions in the rubber material. The distribution of components contained in the rubber material is input to
次に、本実施形態のゴム材料モデル入力工程Ｓ１では、第１高分子成分、第２高分子成分及び架橋剤の分布に基づいて、コンピュータ１が、第１高分子鎖モデル３、第２高分子鎖モデル４、及び、架橋剤粒子モデル５をセル１５内に配置する（工程Ｓ１６）。本実施形態の工程Ｓ１６では、ＤＢＭＣ法（Density Biased Monte Carlo）を用いて、複数の第１高分子鎖モデル３、複数の第２高分子鎖モデル４、及び、複数の架橋剤粒子モデル５が、セル１５内に配置される。ＤＢＭＣ法は、例えば、上記したソフトマテリアル総合シミュレーター（Ｊ－ＯＣＴＡ）に含まれるＣＯＧＮＡＣを用いて計算することができる。
Next, in the rubber material model input step S1 of the present embodiment, the
ＤＢＭＣ法の入力値としては、工程Ｓ６２で求められた第１高分子成分（ブタジエンゴム）の分布（濃度分布）、第２高分子成分（スチレンブタジエンゴム）の分布（濃度分布）、及び、架橋剤の分布（濃度分布）が用いられる。これにより、本実施形態の工程Ｓ１６では、解析対象のゴム材料に含まれる第１高分子成分の分布、第２高分子成分の分布、及び、架橋剤の分布（図７に示す）に基づいて、複数の第１高分子鎖モデル３、複数の第２高分子鎖モデル４、及び、複数の架橋剤粒子モデル５を、セル１５内に配置することができる。図８は、高分子成分の分布に基づいて第１高分子鎖モデル３、及び、第２高分子鎖モデル４が配置されたセル１５を部分的に拡大した概念図である。図８において、図４に示した架橋剤粒子モデル５、第１結合鎖モデル８及び第２結合鎖モデル１２を省略するとともに、第２高分子鎖モデル４を着色して表示している。
The input values for the DBMC method include the distribution (concentration distribution) of the first polymer component (butadiene rubber), the distribution (concentration distribution) of the second polymer component (styrene-butadiene rubber), and the cross-linking The agent distribution (concentration distribution) is used. As a result, in step S16 of the present embodiment, based on the distribution of the first polymer component, the distribution of the second polymer component, and the distribution of the cross-linking agent (shown in FIG. 7) contained in the rubber material to be analyzed, , a plurality of first
図８に示されるように、工程Ｓ１６では、図７に示した各高分子成分の分布に基づいて、複数の第１高分子鎖モデル３（ブタジエンゴム）と、複数の第２高分子鎖モデル４（スチレンブタジエンゴム）との界面構造が作成されたゴム材料モデル１６が定義される。 As shown in FIG. 8, in step S16, a plurality of first polymer chain models 3 (butadiene rubber) and a plurality of second polymer chain models are based on the distribution of each polymer component shown in FIG. 4 (styrene-butadiene rubber) is defined.
次に、本実施形態のゴム材料モデル入力工程Ｓ１では、図４に示されるように、第１高分子鎖モデル３の第１粗視化粒子７、及び、第２高分子鎖モデル４の第２粗視化粒子１１に、引力及び斥力が作用する相互作用ポテンシャルが定義される（工程Ｓ１７）。本実施形態の相互作用ポテンシャルは、一対の第１高分子鎖モデル３、３の第１粗視化粒子７、７間に定義される第５相互作用ポテンシャルＰ５（図示省略）、一対の第２高分子鎖モデル４、４の第２粗視化粒子１１、１１間に定義される第６相互作用ポテンシャルＰ６（図示省略）、及び、第１高分子鎖モデル３の第１粗視化粒子７と第２高分子鎖モデル４の第２粗視化粒子１１との間に定義される第７相互作用ポテンシャルＰ７を含んでいる。第５相互作用ポテンシャルＰ５～第７相互作用ポテンシャルＰ７は、例えば、下記式（１）のＬＪポテンシャルＵLJ
ij（ｒ）で定義することができる。
Next, in the rubber material model input step S1 of the present embodiment, as shown in FIG. An interaction potential is defined in which an attractive force and a repulsive force act on the coarse-grained particles 11 (step S17). The interaction potentials of this embodiment are a fifth interaction potential P5 (not shown) defined between the first coarse-
ｒij：粒子間の距離
ｒij,cut：カットオフ距離
εij：強度パラメータ
σij：粒子の直径に相当
なお、距離ｒij、カットオフ距離ｒij,cut、及び、粒子の直径に相当するσijは、第１粗視化粒子７の中心７ｃ（図５に示す）の座標に基づいて定義される。
r ij : Distance between particles r ij , cut : Cut-off distance ε ij : Strength parameter σ ij : Corresponds to diameter of particle Note that distance r ij , cut-off distance r ij , cut , and diameter of particle σ ij is defined based on the coordinates of the
第５相互作用ポテンシャルＰ５～第７相互作用ポテンシャルＰ７の各定数及び各変数については、例えば、論文２（Ali Makke, Michel Perez, Olivier Lame and Jean-Louis Barrat著、 「Mechanical testing of glassy and rubbery polymers in numerical simulations: Role of boundary conditions in tensile stress experiments 」J. Chem. Phys. vol.131, 014904 (2009)）に基づいて適宜設定される。なお、第５相互作用ポテンシャルＰ５～第７相互作用ポテンシャルＰ７の各強度パラメータεijについては、適宜設定することができる。強度パラメータεijは、ブタジエンゴム及びスチレンブタジエンゴムのχパラメータと相関がある。したがって、第５相互作用ポテンシャルＰ５～第７相互作用ポテンシャルＰ７の各強度パラメータεijは、ブタジエンゴム及びスチレンブタジエンゴムのχパラメータに基づいて同定されるのが望ましい。本実施形態の第５相互作用ポテンシャルＰ５～第７相互作用ポテンシャルＰ７の各強度パラメータεijには、例えば、１．０に設定される。第５相互作用ポテンシャルＰ５～第７相互作用ポテンシャルＰ７は、コンピュータ１に記憶される。
For each constant and each variable of the fifth interaction potential P5 to the seventh interaction potential P7, see, for example, paper 2 (Ali Makke, Michel Perez, Olivier Lame and Jean-Louis Barrat, "Mechanical testing of glassy and rubbery polymers in numerical simulations: Role of boundary conditions in tensile stress experiments” J. Chem. Phys. vol.131, 014904 (2009)). Note that each intensity parameter ε ij of the fifth interaction potential P5 to the seventh interaction potential P7 can be set as appropriate. The strength parameter ε ij correlates with the χ parameter of butadiene rubber and styrene-butadiene rubber. Therefore, each strength parameter ε ij of the fifth interaction potential P5 to the seventh interaction potential P7 is desirably identified based on the χ parameter of butadiene rubber and styrene-butadiene rubber. Each intensity parameter ε ij of the fifth interaction potential P5 to the seventh interaction potential P7 in this embodiment is set to 1.0, for example. The fifth interaction potential P5 to the seventh interaction potential P7 are stored in the
次に、本実施形態のゴム材料モデル入力工程Ｓ１では、架橋剤粒子モデル５に、引力及び斥力が作用する相互作用ポテンシャルが定義される（工程Ｓ１８）。本実施形態の相互作用ポテンシャルは、一対の架橋剤粒子モデル５、５間に定義される第８相互作用ポテンシャルＰ８、架橋剤粒子モデル５と第１高分子鎖モデル３の第１粗視化粒子７との間に定義される第９相互作用ポテンシャルＰ９、及び、架橋剤粒子モデル５と第２高分子鎖モデル４の第２粗視化粒子１１との間に定義される第１０相互作用ポテンシャルＰ１０を含んでいる。第８相互作用ポテンシャルＰ８～第１０相互作用ポテンシャルＰ１０は、例えば、上記式（１）のＬＪポテンシャルＵLJ
ij（ｒ）を用いて定義することができる。
Next, in the rubber material model input step S1 of the present embodiment, interaction potentials where attractive force and repulsive force act on the cross-linking
第８相互作用ポテンシャルＰ８～第１０相互作用ポテンシャルＰ１０の強度パラメータεijについては、適宜設定することができる。強度パラメータεijは、架橋剤のχパラメータと相関がある。従って、強度パラメータεijは、架橋剤のχパラメータに基づいて同定されるのが望ましい。第８相互作用ポテンシャルＰ８～第１０相互作用ポテンシャルＰ１０は、コンピュータ１に記憶される。
The intensity parameters ε ij of the eighth interaction potential P8 to the tenth interaction potential P10 can be set appropriately. The strength parameter ε ij correlates with the χ parameter of the crosslinker. Therefore, the strength parameter ε ij is desirably identified based on the χ parameter of the crosslinker. The eighth interaction potential P8 to the tenth interaction potential P10 are stored in the
次に、本実施形態のゴム材料モデル入力工程Ｓ１では、コンピュータ１が、分子動力学計算に基づいて、第１高分子鎖モデル３及び第２高分子鎖モデル４を対象に構造緩和を計算する（工程Ｓ１９）。本実施形態の分子動力学計算では、例えば、セル１５について所定の時間、第１高分子鎖モデル３の第１粗視化粒子７、及び、第２高分子鎖モデル４の第２粗視化粒子１１が古典力学に従うものとして、ニュートンの運動方程式が適用される。そして、各時刻での第１粗視化粒子７及び第２粗視化粒子１１の動きが、単位時間ステップ毎に追跡される。
Next, in the rubber material model input step S1 of the present embodiment, the
分子動力学計算では、セル１５において、圧力及び温度が一定、又は体積及び温度が一定に保たれる。これにより、工程Ｓ１９では、実際のゴム材料の分子運動に近似させて、第１高分子鎖モデル３、及び、第２高分子鎖モデル４の初期配置を精度よく緩和することができる。本実施形態では、第１高分子鎖モデル３、及び、第２高分子鎖モデル４の初期配置が十分に緩和されるまで行われる。構造緩和の計算は、例えば、上記Ｊ－ＯＣＴＡに含まれるＣＯＧＮＡＣを用いて処理することができる。
In molecular dynamics calculations, the pressure and temperature, or the volume and temperature, are kept constant in the
本実施形態の工程Ｓ１９では、第１粗視化粒子７、７間に定義される第３相互作用ポテンシャルＰ３（図５に示す）、及び、第２粗視化粒子１１、１１間に定義される第４相互作用ポテンシャルＰ４（図５に示す）が、非調和ポテンシャルＵch（ｒ）として定義されている。これにより、工程Ｓ１９では、第１粗視化粒子７、７間、及び、第２粗視化粒子１１、１１間で、柔軟に伸縮することができるため、構造緩和をスムーズに行うことができる。さらに、本実施形態の第３相互作用ポテンシャルＰ３及び第４相互作用ポテンシャルＰ４の各定数及び各変数には、同一のものが採用されている。このため、セル１５において、第１高分子鎖モデル３の配置、及び、第２高分子鎖モデル４の配置に偏りが生じるのを防ぐことができる。
In step S19 of the present embodiment, a third interaction potential P3 (shown in FIG. 5) defined between the first coarse-
本実施形態の工程Ｓ１９では、架橋剤粒子モデル５を対象とした構造緩和が計算されていない。これにより、架橋剤粒子モデル５は、セル１５において、実際の架橋剤の分布を維持することができる。
Structural relaxation for the cross-linking
次に、本実施形態のゴム材料モデル入力工程Ｓ１は、コンピュータ１が、隣接する分子鎖モデル（第１高分子鎖モデル３又は第２高分子鎖モデル４）について、一対の粗視化粒子（第１粗視化粒子７又は第２粗視化粒子１１）を連結する（工程Ｓ２０）。本実施形態の工程Ｓ２０では、一対の粗視化粒子（第１粗視化粒子７又は第２粗視化粒子１１）を、架橋剤粒子モデル５を介して連結している。図９（ａ）、（ｂ）は、一対の粗視化粒子を連結する工程の一例を説明する概念図である。
Next, in the rubber material model input step S1 of the present embodiment, the
図９（ａ）、（ｂ）に示されるように、工程Ｓ２０では、架橋剤粒子モデル５の位置から予め定められた範囲Ｒ（図９（ａ）に示す）内に存在する一対の粗視化粒子（本例では、第１粗視化粒子７、７）を連結している。範囲Ｒについては、適宜設定することができる。本実施形態の範囲Ｒは、例えば、架橋剤粒子モデル５の中心５ｃからの半径ｒｄが、０．８σ以内の領域に設定されている。
As shown in FIGS. 9A and 9B, in step S20, a pair of coarse grains existing within a predetermined range R (shown in FIG. 9A) from the position of the cross-linking
工程Ｓ２０では、図９（ｂ）に示されるように、一対の粗視化粒子（本例では、第１粗視化粒子７、７）が、架橋剤粒子モデル５を介して連結される。本実施形態では、一方の粗視化粒子（第１粗視化粒子７又は第２粗視化粒子１１）と架橋剤粒子モデル５との間、及び、他方の粗視化粒子（第１粗視化粒子７又は第２粗視化粒子１１）と架橋剤粒子モデル５との間が、第３結合鎖モデル１８を介して接続される。
In step S20, a pair of coarse-grained particles (first coarse-
第３結合鎖モデル１８は、第１１相互作用ポテンシャルＰ１１によって設定される。本実施形態の第１１相互作用ポテンシャルＰ１１は、上記式（２）の非調和ポテンシャルＵch（ｒ）によって定義される。また、上記式（２）において、第１１相互作用ポテンシャルＰ１１の定数及び変数は、例えば、図４に示した第３相互作用ポテンシャルＰ３及び第４相互作用ポテンシャルＰ４の定数又は変数と同様に定義されうる。
The third
このように、工程Ｓ２０では、隣接する高分子鎖モデル（本実施形態では、第１高分子鎖モデル３又は第２高分子鎖モデル４）の一対の粗視化粒子（本実施形態では、第１粗視化粒子７又は第２粗視化粒子１１）が、架橋剤粒子モデル５を介して連結される。これにより、解析対象の架橋されたゴム材料を再現したゴム材料モデル１６（図４に示す）を作成することができる。ゴム材料モデル１６は、コンピュータ１に記憶される。
Thus, in step S20, a pair of coarse-grained particles (in this embodiment, the second One coarse-
図４に示されるように、第１高分子鎖モデル３及び第２高分子鎖モデル４について、第１粗視化粒子７、７間に定義される第３相互作用ポテンシャルＰ３、及び、第２粗視化粒子１１、１１間に定義される第４相互作用ポテンシャルＰ４は、非調和ポテンシャルＵch（ｒ）によって定義されている。さらに、第３相互作用ポテンシャルＰ３及び第４相互作用ポテンシャルＰ４の各定数及び各変数には、同一の値が設定されている。このため、第１高分子鎖モデル３及び第２高分子鎖モデル４は、第１高分子成分及び第２高分子成分の物性値の差異が考慮されていない。したがって、工程Ｓ１１～工程Ｓ２０を経て作成された相分離界面（即ち、第１高分子成分及び第２高分子成分が互いに分離した界面）のゴム材料モデル１６は、例えば、後述の変形（引張試験）の計算に用いられると、上記非特許文献１と同様に、第１高分子鎖モデル３が配置された領域と、第２高分子鎖モデル４が配置された領域との界面で破壊が生じる現象が計算される。
As shown in FIG. 4, for the first
しかしながら、第１高分子成分及び第２高分子成分を含む相分離した実際のゴム材料を変形させると、第１高分子成分及び第２高分子成分のうち、耐変形性を示す物性値（例えば、弾性率）が相対的に小さい第１高分子成分又は第２高分子成分の一方（本実施形態では、第１高分子成分）内で破壊が観察される。 However, when the phase-separated actual rubber material containing the first polymer component and the second polymer component is deformed, the physical property value (for example, , elastic modulus) is relatively small (in this embodiment, the first polymer component).
発明者らは、鋭意研究を重ねた結果、第１高分子成分及び第２高分子成分の物性値の差異を考慮して、第１粗視化粒子７、７間の相互作用、及び、第２粗視化粒子１１、１１間の相互作用を互いに異なるものとして定義することが、第１高分子成分及び第２高分子成分を含む実際のゴム材料の破壊特性を評価する上で有効であることを知見した。
As a result of intensive studies, the inventors have found that the interaction between the first coarse-
このような知見に基づいて、本実施形態のシミュレーション方法（ゴム材料モデル入力工程Ｓ１）では、第１高分子成分と第２高分子成分との物性値の差異に基づいて、第２粗視化粒子１１、１１間の相互作用を定義するための第２相互作用ポテンシャルＰ２を、第１粗視化粒子７、７間の相互作用を定義するための第１相互作用ポテンシャルＰ１とは異なるものとして定義している。
Based on such knowledge, in the simulation method (rubber material model input step S1) of the present embodiment, the second coarse-grained Assuming that the second interaction potential P2 for defining the interaction between the
本実施形態のゴム材料モデル入力工程Ｓ１では、第１粗視化粒子７、７間の相互作用を定義するための第１相互作用ポテンシャルＰ１（図４に示す）が定義される（工程Ｓ２１）。第１相互作用ポテンシャルＰ１は、相互作用の強さの定義するための強度パラメータεijを含み、上記式（１）のＬＪポテンシャルによって定義される。本実施形態の工程Ｓ２１では、上記式（２）で定義された第１粗視化粒子７、７間の第３相互作用ポテンシャルＰ３（図４に示す）を無効にし、上記式（１）で定義される第１相互作用ポテンシャルＰ１（図４に示す）が定義される。
In the rubber material model input step S1 of the present embodiment, a first interaction potential P1 (shown in FIG. 4) for defining the interaction between the first coarse-
第１相互作用ポテンシャルＰ１の各定数及び各変数については、例えば、上記論文２に基づいて適宜設定される。本実施形態の第１相互作用ポテンシャルＰ１の強度パラメータεijは、１．０に設定される。第１相互作用ポテンシャルＰ１は、コンピュータ１に記憶される。
Each constant and each variable of the first interaction potential P1 are appropriately set based on the above-mentioned article 2, for example. The strength parameter ε ij of the first interaction potential P1 in this embodiment is set to 1.0. The first interaction potential P1 is stored in
次に、本実施形態のゴム材料モデル入力工程Ｓ１では、第２粗視化粒子１１、１１間の相互作用を定義するための第２相互作用ポテンシャルＰ２（図４に示す）が定義される（第２相互作用定義工程Ｓ２２）。第２相互作用ポテンシャルＰ２は、相互作用の強さの定義するための強度パラメータεijを含み、上記式（１）のＬＪポテンシャルによって定義される。第２相互作用定義工程Ｓ２２では、第１高分子成分と第２高分子成分との物性値の差異に基づいて、第１相互作用ポテンシャルＰ１（図４に示す）とは異なるものとして、第２相互作用ポテンシャルＰ２を定義している。図１０は、第２相互作用定義工程Ｓ２２の処理手順の一例を示すフローチャートである。
Next, in the rubber material model input step S1 of the present embodiment, a second interaction potential P2 (shown in FIG. 4) for defining the interaction between the second coarse-
本実施形態の第２相互作用定義工程Ｓ２２では、先ず、上記式（２）で定義された第２粗視化粒子１１、１１間の第４相互作用ポテンシャルＰ４（図４に示す）を無効にする工程Ｓ３１、及び、第２粗視化粒子１１、１１間に、上記式（１）で定義される第２相互作用ポテンシャルＰ２（図４に示す）を定義する工程Ｓ３２が実施される。
In the second interaction definition step S22 of the present embodiment, first, the fourth interaction potential P4 (shown in FIG. 4) between the second coarse-
本実施形態の工程Ｓ３２において、第２相互作用ポテンシャルＰ２の各定数及び各変数については、例えば、上記論文２に基づいて適宜設定される。なお、強度パラメータεijは、後述の強度パラメータ定義工程Ｓ３３において、第１相互作用ポテンシャルＰ１とは異なる強度パラメータεijが定義される。 In the step S32 of the present embodiment, each constant and each variable of the second interaction potential P2 are appropriately set based on the above-mentioned article 2, for example. Note that the strength parameter ε ij different from the first interaction potential P1 is defined in the strength parameter definition step S33, which will be described later.
次に、本実施形態の第２相互作用定義工程Ｓ２２では、第２相互作用ポテンシャルＰ２に、第１相互作用ポテンシャルＰ１とは異なる強度パラメータεijが定義される（強度パラメータ定義工程Ｓ３３）。図１１は、強度パラメータ定義工程Ｓ３３の処理手順の一例を示すフローチャートである。 Next, in the second interaction definition step S22 of the present embodiment, a strength parameter ε ij different from the first interaction potential P1 is defined for the second interaction potential P2 (strength parameter definition step S33). FIG. 11 is a flow chart showing an example of the processing procedure of the intensity parameter definition step S33.
本実施形態の強度パラメータ定義工程Ｓ３３では、先ず、予め定められた複数種類の強度パラメータεij毎に、物性値を計算する（物性値計算工程Ｓ４１）。図１２は、物性値計算工程Ｓ４１の処理手順の一例を示すフローチャートである。図１３は、高分子鎖モデル１７が配置されたセル２１の一例を示す概念図である。
In the strength parameter definition step S33 of the present embodiment, first, a physical property value is calculated for each of a plurality of types of predetermined strength parameters ε ij (physical property value calculation step S41). FIG. 12 is a flow chart showing an example of the processing procedure of the physical property value calculation step S41. FIG. 13 is a conceptual diagram showing an example of the
本実施形態の物性値計算工程Ｓ４１では、先ず、図１３に示されるように、物性値を計算するための高分子鎖モデル１７が、コンピュータ１に入力される（工程Ｓ５１）。高分子鎖モデル１７は、図４に示した第１高分子鎖モデル３及び第２高分子鎖モデル４と同様に、複数の粗視化粒子１９と、隣接する粗視化粒子１９、１９間を結合する結合鎖モデル２０とを含んで構成されている。高分子鎖モデル１７には、図４に示した第１高分子鎖モデル３又は第２高分子鎖モデル４がそのまま用いられてもよいし、第１高分子鎖モデル３及び第２高分子鎖モデル４とは独立して定義されてもよい。
In the physical property value calculation step S41 of the present embodiment, first, as shown in FIG. 13, a
粗視化粒子１９は、後述の分子動力学計算において、運動方程式の質点として取り扱われる。このため、粗視化粒子１９には、図４に示した第１粗視化粒子７及び第２粗視化粒子１１と同様のパラメータが定義される。
The coarse-
結合鎖モデル２０は、粗視化粒子１９、１９間に設定される第１５相互作用ポテンシャルＰ１５によって定義される。第１５相互作用ポテンシャルＰ１５は、図４に示した第１相互作用ポテンシャルＰ１及び第２相互作用ポテンシャルＰ２と同様に、上記式（１）のＬＪポテンシャルによって定義される。第１５相互作用ポテンシャルＰ１５の各定数及び各変数については、例えば、強度パラメータεijを除いて、例えば、第１相互作用ポテンシャルＰ１及び第２相互作用ポテンシャルＰ２と同一のものを設定することができる。高分子鎖モデル１７は、コンピュータ１に記憶される。
次に、本実施形態の物性値計算工程Ｓ４１では、セル２１が、コンピュータ１に入力される（工程Ｓ５２）。本実施形態のセル２１は、図４に示したゴム材料モデル１６のセル１５とは独立して定義される。セル２１は、ゴム材料モデル１６のセル１５と同一の大きさに設定されても良いし、異なる大きさに設定されてもよい。セル２１は、コンピュータ１に記憶される。
Next, in the physical property value calculation step S41 of this embodiment, the
次に、本実施形態の物性値計算工程Ｓ４１では、複数の高分子鎖モデル１７がセル２１内に配置され（工程Ｓ５３）、高分子鎖モデル１７の粗視化粒子１９に、引力及び斥力が作用する第１６相互作用ポテンシャルＰ１６が定義される（工程Ｓ５４）。本実施形態の第１６相互作用ポテンシャルＰ１６は、一対の高分子鎖モデル１７、１７の粗視化粒子１９、１９間にそれぞれ定義される。第１６相互作用ポテンシャルＰ１６は、第５相互作用ポテンシャルＰ５～第７相互作用ポテンシャルＰ７（図４に示す）と同様に、上記式（１）のＬＪポテンシャルＵLJ
ij（ｒ）で定義することができる。
Next, in the physical property value calculation step S41 of the present embodiment, a plurality of
第１６相互作用ポテンシャルＰ１６の各定数及び各変数については、例えば、上記論文２に基づいて適宜設定される。本実施形態の第１６相互作用ポテンシャルＰ１６の強度パラメータεijは、１．０に設定される。第１６相互作用ポテンシャルＰ１６は、コンピュータ１に記憶される。
Each constant and each variable of the sixteenth interaction potential P16 are appropriately set based on, for example, Paper 2 above. The strength parameter ε ij of the sixteenth interaction potential P16 in this embodiment is set to 1.0. The sixteenth interaction potential P16 is stored in
次に、本実施形態の物性値計算工程Ｓ４１では、コンピュータ１が、分子動力学計算に基づいて、高分子鎖モデル１７を対象に構造緩和を計算する（工程Ｓ５５）。分子動力学計算については、上述のとおりである。本実施形態では、高分子鎖モデル１７の初期配置が十分に緩和されるまで行われる。これにより、単一の高分子成分を含むゴム材料をモデル化したゴム材料モデル（以下、単に「単一ゴム材料モデル」ということがある。）２２が定義される
Next, in the physical property value calculation step S41 of the present embodiment, the
次に、本実施形態の物性値計算工程Ｓ４１は、コンピュータ１が、予め定められた複数種類の強度パラメータεij毎に、物性値を計算する（工程Ｓ５６）。工程Ｓ５６では、高分子鎖モデル１７の第１５相互作用ポテンシャルＰ１５に設定される複数種類の強度パラメータεij毎に、単一ゴム材料モデル２２の伸長計算をして、物性値が計算される。
Next, in the physical property value calculation step S41 of the present embodiment, the
複数種類の強度パラメータεijについては、例えば、第１高分子成分と第２高分子成分との物性値の差異に基づいて、適宜設定することができる。強度パラメータεijの一例としては、第１高分子成分と第２高分子成分との物性値に差異が無い状態の強度パラメータεij（例えば、１．０）から、第１高分子成分と第２高分子成分との物性値の差異が大きくなる強度パラメータεij（例えば、２．４）の間で複数設定（例えば、εij＝１．０、１．２、１．４、１．６、１．８、２．０、２．２、及び、２．４）されるのが望ましい。なお、強度パラメータεijは、このような態様に限定されない。 A plurality of types of intensity parameters ε ij can be appropriately set, for example, based on the difference in physical property values between the first polymer component and the second polymer component. As an example of the strength parameter ε ij , the strength parameter ε ij (for example, 1.0) when there is no difference in physical property values between the first polymer component and the second polymer component is Multiple settings (e.g., ε ij =1.0, 1.2, 1.4, 1.6) among strength parameters ε ij (e.g., 2.4) that increase the difference in physical property values from the two polymer components , 1.8, 2.0, 2.2 and 2.4). Note that the strength parameter ε ij is not limited to such a mode.
工程Ｓ５６では、先ず、複数種類の強度パラメータεijから選択された一つの強度パラメータεijが、高分子鎖モデル１７の第１５相互作用ポテンシャルＰ１５に定義される。次に、工程Ｓ５６では、単軸引張り試験に基づいて、単一ゴム材料モデル２２を一方向に（例えば、Ｙ軸方向に０％～５０％）伸長させる変形シミュレーションが実施され、その物性値（本実施形態では、弾性率）が計算される。変形シミュレーションについては、後述のゴム材料モデル１６を用いたゴム材料のシミュレーション方法と同様の手順で行うことができる。
In step S56, first, one strength parameter ε ij selected from a plurality of types of strength parameters ε ij is defined as the fifteenth interaction potential P15 of the
工程Ｓ５６において、物性値の計算は、設定された複数の強度パラメータεij毎に実施される。図１４は、強度パラメータεij毎に計算された応力－歪み曲線の一例を示すグラフである。図１４では、強度パラメータεijが１．０、１．６、及び、２．４のみを代表して示している。図１５は、物性値の差異（ヤング率の比Ｅ／Ｅ1.0）と、強度パラメータεijとの関係を示すグラフである。ヤング率の比Ｅ／Ｅ1.0は、強度パラメータεij＝１．０のヤング率Ｅ1.0と他の強度パラメータεijのヤング率Ｅとの割合を示している。複数種類の強度パラメータεij毎に計算された物性値は、コンピュータ１に記憶される。
In step S56, physical property values are calculated for each of the plurality of set strength parameters ε ij . FIG. 14 is a graph showing an example of a stress-strain curve calculated for each strength parameter ε ij . FIG. 14 representatively shows only intensity parameters ε ij of 1.0, 1.6, and 2.4. FIG. 15 is a graph showing the relationship between the physical property value difference (Young's modulus ratio E/E 1.0 ) and the strength parameter ε ij . The Young's modulus ratio E/E 1.0 indicates the ratio of the Young's modulus E 1.0 for the strength parameter ε ij =1.0 to the Young's modulus E for the other strength parameter ε ij . Physical property values calculated for each of the plurality of types of intensity parameters ε ij are stored in the
次に、図１１に示されるように、本実施形態の強度パラメータ定義工程Ｓ３３は、物性値の差異に基づいて、複数種類の強度パラメータεij（例えば、εij＝１．０、１．２、１．４、１．６、１．８、２．０、２．２、及び、２．４）から、第２相互作用ポテンシャルＰ２の強度パラメータεijを選択する（工程Ｓ４２）。工程Ｓ４２では、先ず、実際の第１高分子成分と第２高分子成分との物性値の差異が求められる。物性値（弾性率）の差異としては、第１高分子成分のヤング率Ｅ1と第２高分子成分のヤング率Ｅ2との比Ｅ2／Ｅ1が求められる。第１高分子成分のヤング率Ｅ1と第２高分子成分のヤング率Ｅ2は、例えば、ＪＩＳＫ６２５１「加硫ゴム及び熱可塑性ゴム－引張特性の求め方」に記載の試験方法に準拠して引張特性を求め、歪みに対して応力が直線的に変化する範囲（即ち、線形応答が得られた範囲）内の引張応力と、これに対応する歪みとの比から求めることができる。 Next, as shown in FIG. 11, in the strength parameter definition step S33 of the present embodiment, a plurality of types of strength parameters ε ij (eg, ε ij =1.0, 1.2) are calculated based on differences in physical property values. , 1.4, 1.6, 1.8, 2.0, 2.2, and 2.4), select the strength parameter ε ij of the second interaction potential P2 (step S42). In step S42, first, the actual difference in physical property values between the first polymer component and the second polymer component is obtained. As the difference in the physical property value ( elastic modulus), the ratio E2 / E1 between the Young's modulus E1 of the first polymer component and the Young's modulus E2 of the second polymer component is obtained. The Young's modulus E 1 of the first polymer component and the Young's modulus E 2 of the second polymer component are, for example, in accordance with the test method described in JISK6251 "Vulcanized rubber and thermoplastic rubber - Determination of tensile properties". Tensile properties are determined and can be determined from the ratio of the tensile stress within the range in which the stress varies linearly with strain (ie, the range in which a linear response is obtained) to the corresponding strain.
次に、本実施形態の工程Ｓ４２では、図１５のグラフにおいて、実際の第１高分子成分と第２高分子成分との物性値の差異（ヤング率の比Ｅ2／Ｅ1）に対応する強度パラメータεijが選択される。図１５のグラフにおいて、例えば、物性値の差異（ヤング率の比Ｅ2／Ｅ1）が２．５の場合、強度パラメータεij＝２．０が選択される。 Next, in step S42 of the present embodiment, in the graph of FIG . 15 , there is An intensity parameter ε ij is chosen. In the graph of FIG. 15, for example, when the difference in physical property values (Young's modulus ratio E 2 /E 1 ) is 2.5, the strength parameter ε ij =2.0 is selected.
次に、本実施形態の強度パラメータ定義工程Ｓ３３では、第２相互作用ポテンシャルＰ２（図４に示す）に、選択された強度パラメータ（本例では、εij＝２．０）が定義される（工程Ｓ４３）。強度パラメータεij＝２．０は、第１相互作用ポテンシャルＰ１（図４に示す）の強度パラメータεij＝１．０とは異なるものである。これにより、第２相互作用定義工程Ｓ２２では、第２相互作用ポテンシャルＰ２を、第１高分子成分と第２高分子成分との物性値の差異（本実施形態では、ヤング率の比Ｅ2／Ｅ1）に基づいて、第１相互作用ポテンシャルＰ１とは異なるものとして定義することができる。したがって、本実施形態のシミュレーション方法（ゴム材料モデル入力工程Ｓ１）は、第１高分子成分と第２高分子成分との物性値の差異を考慮することが可能なゴム材料モデル１６（図４に示す）を作成することができる。ゴム材料モデル１６は、コンピュータ１に記憶される。
Next, in the strength parameter definition step S33 of the present embodiment, the selected strength parameter (ε ij =2.0 in this example) is defined for the second interaction potential P2 (shown in FIG. 4) ( step S43). The strength parameter ε ij =2.0 is different from the strength parameter ε ij =1.0 of the first interaction potential P1 (shown in FIG. 4). Thus, in the second interaction definition step S22, the second interaction potential P2 is defined as the difference in physical property values between the first polymer component and the second polymer component (in this embodiment, the ratio of Young's modulus E 2 / E 1 ) can be defined as different from the first interaction potential P1. Therefore, the simulation method (rubber material model input step S1) of the present embodiment uses the rubber material model 16 (Fig. 4 shown) can be created. A
本実施形態のゴム材料モデル入力工程Ｓ１では、第２相互作用定義工程Ｓ２２の後に、コンピュータ１が、分子動力学計算に基づいて、図４に示した第１高分子鎖モデル３及び第２高分子鎖モデル４を対象に構造緩和を計算するのが望ましい。これにより、第２相互作用定義工程Ｓ２２で設定される第２相互作用ポテンシャルＰ２の変化によって、分子同士の凝集度合い、セル２１のサイズ、セル２１内の密度（即ち、セル２１内の粒子数／セル２１のサイズ）が変化しても、新たな第２相互作用ポテンシャルＰ２に基づいて平衡化された構造を作成することができるため、セル２１を適切な密度にすることができる。本実施形態では、NPT（N：粒子数、P：圧力、T：温度）一定の条件下で構造緩和が行われるのが望ましく、また、セル１５の圧力がゼロ（又は、ゼロに近い任意の値）になるように、構造緩和が計算されるのが望ましい。これにより、ポテンシャルの異なる構造間で初期応力が異なることを避けることができ、界面状態の異なる構造間での応力比較が可能になる。
In the rubber material model input step S1 of the present embodiment, after the second interaction definition step S22, the
本実施形態のゴム材料モデル入力工程Ｓ１では、架橋剤をモデル化した架橋剤粒子モデル５が、セル１５に配置されたが、このような態様に限定されない。本実施形態のゴム材料モデル入力工程Ｓ１では、架橋剤粒子モデル５（図４及び図５に示す）を省略してもよい。この場合、工程Ｓ２０では、架橋剤粒子モデル５を介さずに、一対の粗視化粒子（第１粗視化粒子７又は第２粗視化粒子１１）を連結するのが望ましい。また、ゴム材料に配合される充填剤（フィラー）をモデル化したフィラーモデル（図示省略）が、セル１５に配置されてもよい。これにより、充填剤が配合されたゴム材料を再現したゴム材料モデル１６を作成することができる。
In the rubber material model input step S1 of the present embodiment, the cross-linking
次に、本実施形態のシミュレーション方法では、コンピュータ１が、図４に示したゴム材料モデル１６の変形を計算（以下、単に、「変形計算」ということがある。）する（工程Ｓ２）。ゴム材料モデル１６の変形計算としては、例えば、引張変形や圧縮変形など適宜選択することができる。本実施形態の変形計算は、引張変形が行われる。
Next, in the simulation method of this embodiment, the
工程Ｓ２では、Ｙ軸方向において、ゴム材料モデル１６の一端１６ａ及び他端１６ｂを離間させて、ゴム材料モデル１６の引張変形が計算される。本実施形態では、例えば、特開２０１６－８１２９７号公報に記載の方法と同様に、分子動力学計算の単位時間（ＭＤステップ）あたりの歪が１０－５以上の速度Ｖ１で、ゴム材料モデル１６をアフィン変形に基づいて引張変形させている。また、工程Ｓ２では、第１高分子鎖モデル３、第２高分子鎖モデル４、及び、架橋剤粒子モデル５の熱運動（分子運動）が計算される。
In step S2, the tensile deformation of the
次に、本実施形態のシミュレーション方法では、ゴム材料モデル１６の応力が計算される（工程Ｓ３）。工程Ｓ３では、工程Ｓ２で与えられた歪に基づいて、ゴム材料モデル１６の応力が単位時間毎に計算される。さらに、工程Ｓ３では、ゴム材料モデル１６の引張変形によって計算される第１粗視化粒子７、第２粗視化粒子１１、及び、架橋剤粒子モデル５の動きが、単位時間毎に追跡される。図１６（Ａ）は、ゴム材料モデル（歪：０．０）を示す図である。図１６（Ｂ）は、ゴム材料モデル（歪：０．１２）を示す図である。図１６（Ｃ）は、ゴム材料モデル（歪：０．３０）を示す図である。図１６（Ｄ）は、ゴム材料モデル（歪：０．６０）を示す図である。
Next, in the simulation method of this embodiment, the stress of the
工程Ｓ３では、ゴム材料モデル１６に与えられた歪、各相互作用ポテンシャル及び運動方程式に基づいて、第１高分子鎖モデル３、第２高分子鎖モデル４、及び、架橋剤粒子モデル５の熱運動が計算される。これらの熱運動が計算されることにより、セル１５には、第１高分子鎖モデル３、第２高分子鎖モデル４、及び、架橋剤粒子モデル５（図示省略）が存在しない小領域が形成される。このような小領域は、ゴム材料モデル１６に形成された空孔（ボイド）２４として定義される。
In step S3, the heat of the first
図１６（Ａ）～（Ｄ）に示されるように、本実施形態の変形計算では、第２高分子鎖モデル４よりも小さい強度パラメータεijが定義された第１高分子鎖モデル３において、破壊が生じる（空孔２４が生じる）現象が計算された。これは、実際のゴム材料の破壊現象に合致している。したがって、本実施形態のシミュレーション方法では、ゴム材料モデル入力工程Ｓ１で作成されたゴム材料モデル１６が、変形計算を行うシミュレーションに用いられることにより、第１高分子成分と第２高分子成分とを含むゴム材料の破壊特性を、精度良く評価することができる。ゴム材料モデル１６の応力、第１高分子鎖モデル３、第２高分子鎖モデル４、及び、架橋剤粒子モデル５の各座標値は、コンピュータ１に記憶される。
As shown in FIGS. 16A to 16D, in the deformation calculation of this embodiment, in the first
次に、本実施形態のシミュレーション方法では、ゴム材料モデル１６の変形の計算結果に基づいて、ゴム材料の性能が評価される（工程Ｓ４）。本実施形態の工程Ｓ４では、例えば、ゴム材料モデル１６の変形結果に基づいて計算される応力－歪み曲線や、空孔（ボイド）２４の成長曲線等に基づいて、ゴム材料の強度が評価される。
Next, in the simulation method of this embodiment, the performance of the rubber material is evaluated based on the deformation calculation result of the rubber material model 16 (step S4). In step S4 of the present embodiment, for example, the strength of the rubber material is evaluated based on the stress-strain curve calculated based on the deformation result of the
本実施形態では、ゴム材料の性能（変形の計算結果）が良好であると評価された場合（工程Ｓ４において、「Ｙ」）、ゴム材料モデル１６に基づいて、ゴム材料が製造される（工程Ｓ５）。一方、ゴム材料の性能が良好でないと判断された場合（工程Ｓ４において、「Ｎ」）、ゴム材料の諸条件を変更し（工程Ｓ６）、かつ、図３に示した手順に基づいて、ゴム材料モデル１６を再作成して（工程Ｓ７）、及び、工程Ｓ１～工程Ｓ４が再度実施される。このように、本実施形態のシミュレーション方法、及び、製造方法では、ゴム材料の性能が良好になるまで、ゴム材料の諸条件が変更されるため、良好な強度（破壊特性）を有するゴム材料を、効率よく製造することができる。
In the present embodiment, when the performance of the rubber material (deformation calculation result) is evaluated as good ("Y" in step S4), the rubber material is manufactured based on the rubber material model 16 (step S5). On the other hand, if it is determined that the performance of the rubber material is not good (“N” in step S4), various conditions of the rubber material are changed (step S6), and based on the procedure shown in FIG. The
以上、本発明の特に好ましい実施形態について詳述したが、本発明は図示の実施形態に限定されることなく、種々の態様に変形して実施しうる。 Although the particularly preferred embodiments of the present invention have been described in detail above, the present invention is not limited to the illustrated embodiments and can be modified in various ways.
予め定められた物性値について差異を有する第１高分子成分と第２高分子成分とを少なくとも含むゴム材料の数値解析用のゴム材料モデルが作成された（実施例、比較例）。 A rubber material model for numerical analysis of a rubber material containing at least a first polymer component and a second polymer component having predetermined physical property values different from each other was created (Example, Comparative Example).
実施例及び比較例では、第１高分子成分の分子鎖構造を複数の第１粗視化粒子を用いて表現した第１高分子鎖モデル、及び、第２高分子成分の分子鎖構造を複数の第２粗視化粒子を用いて表現した第２高分子鎖モデルが、コンピュータに入力された。第１高分子鎖モデルには、第１粗視化粒子間の相互作用を定義するための第１相互作用ポテンシャルが定義された。第２高分子鎖モデルには、第２粗視化粒子間の相互作用を定義するための第２相互作用ポテンシャルが定義された。さらに、実施例及び比較例では、架橋剤をモデル化した架橋剤粒子モデルが、コンピュータに入力された。 In Examples and Comparative Examples, a first polymer chain model expressing the molecular chain structure of the first polymer component using a plurality of first coarse-grained particles, and a plurality of molecular chain structures of the second polymer component were used. was input into the computer. A first interaction potential for defining the interaction between the first coarse-grained particles was defined in the first polymer chain model. A second interaction potential was defined in the second polymer chain model to define the interaction between the second coarse-grained particles. Furthermore, in Examples and Comparative Examples, a cross-linking agent particle model modeling a cross-linking agent was input into the computer.
そして、実施例及び比較例では、図７に示した第１高分子成分の分布、第２高分子成分の分布、及び、架橋剤の分布に基づいて、第１高分子鎖モデル、第２高分子鎖モデル、及び、架橋剤粒子モデルがセル内に配置され、ゴム材料モデルが作成された。 Then, in the examples and comparative examples, the first polymer chain model, the second polymer chain model, the second polymer chain model, and the second polymer chain model based on the distribution of the first polymer component, the distribution of the second polymer component, and the distribution of the cross-linking agent shown in FIG. A molecular chain model and a cross-linking agent particle model were placed in a cell to create a rubber material model.
実施例の第２相互作用ポテンシャルは、図１１に示した強度パラメータ定義工程の処理手順にしたがい、第１高分子成分と第２高分子成分との物性値の差異に基づいて、第１相互作用ポテンシャルとは異なるものとして定義された。一方、比較例の第２相互作用ポテンシャルは、物性値の差異を考慮せずに、第１相互作用ポテンシャルと同一のものが採用された。そして、実施例及び比較例のゴム材料モデルの引張変形が計算された。共通仕様は、次のとおりである。 The second interaction potential of the example is calculated based on the difference in physical property values between the first polymer component and the second polymer component according to the processing procedure of the strength parameter definition step shown in FIG. defined as different from potential. On the other hand, the second interaction potential of the comparative example was the same as the first interaction potential without considering the difference in physical property values. Then, the tensile deformation of the rubber material models of Examples and Comparative Examples was calculated. Common specifications are as follows.
第１、第２高分子鎖モデル：
第１、第２粗視化粒子の個数（鎖長）：２００個
各本数：４０本
第３、第４相互作用ポテンシャル：
粒子モデル間のばね定数ｋ：３０［ε／σ2］
伸びきり長Ｒ0：１［σ］
架橋剤粒子モデル：２２９個
セル：
体積：１６２２９［σ3］
実施例の密度：０．９４［／σ3］
比較例の密度：０．８７［／σ3］
架橋密度：０．０１［σ3］
第５～第１０相互作用ポテンシャル：
カットオフ距離ｒij,cut：２．５［σ］
強度パラメータεij：１．０［ε］
第１１相互作用ポテンシャル：
粒子モデル間のばね定数ｋ：３０［ε／σ2］
伸びきり長Ｒ0：１［σ］
強度パラメータ定義工程：
高分子鎖モデル：
粗視化粒子の個数（鎖長）：２００個
分子数：４０本
第１５相互作用ポテンシャル：
粒子モデル間のばね定数ｋ：３０［ε／σ2］
伸びきり長Ｒ0：１［σ］
強度パラメータε：１．０、１．２、１．４、１．６、
１．８、２．０、２．２、２．４［ε］
第１６相互作用ポテンシャル：
カットオフ距離ｒij,cut：２．５［σ］
強度パラメータεij：１．０［ε］
粒子の直径に相当σij：１［σ］
第１相互作用ポテンシャル：
強度パラメータεij：１［ε］
カットオフ距離ｒij,cut：２．５［σ］
粒子の直径に相当σij：１［σ］
第２相互作用ポテンシャル：
強度パラメータεij（実施例）：２［ε］
強度パラメータεij（比較例）：１［ε］
カットオフ距離ｒij,cut：２．５［σ］
粒子の直径に相当σij：１［σ］
計算条件：
時間刻み（単位時間）：０．００６［τ］
合計時間：６０００［τ］
ひずみ速度：０．０００１［１／τ］
最大歪み：０．６
ポアソン比：０
First and second polymer chain models:
Number of first and second coarse-grained particles (chain length): 200
Each number: 40
Third and fourth interaction potentials:
Spring constant k between particle models: 30 [ε/σ 2 ]
Stretched length R 0 : 1 [σ]
Crosslinker particle model: 229 cells:
Volume: 16229 [σ 3 ]
Example density: 0.94 [/σ 3 ]
Density of Comparative Example: 0.87 [/σ 3 ]
Crosslink density: 0.01 [σ 3 ]
5th to 10th interaction potentials:
Cutoff distance r ij , cut : 2.5 [σ]
Strength parameter ε ij : 1.0 [ε]
Eleventh interaction potential:
Spring constant k between particle models: 30 [ε/σ 2 ]
Stretched length R 0 : 1 [σ]
Strength parameter definition step:
Polymer chain model:
Number of coarse-grained particles (chain length): 200
Number of molecules: 40
Fifteenth interaction potential:
Spring constant k between particle models: 30 [ε/σ 2 ]
Stretched length R 0 : 1 [σ]
strength parameter ε: 1.0, 1.2, 1.4, 1.6,
1.8, 2.0, 2.2, 2.4 [ε]
Sixteenth interaction potential:
Cutoff distance r ij , cut : 2.5 [σ]
Strength parameter ε ij : 1.0 [ε]
Corresponding to particle diameter σ ij : 1 [σ]
First interaction potential:
Intensity parameter ε ij : 1 [ε]
Cutoff distance r ij , cut : 2.5 [σ]
Corresponding to particle diameter σ ij : 1 [σ]
Second interaction potential:
Strength parameter ε ij (Example): 2 [ε]
Strength parameter ε ij (comparative example): 1 [ε]
Cutoff distance r ij , cut : 2.5 [σ]
Corresponding to particle diameter σ ij : 1 [σ]
Calculation condition:
Time step (unit time): 0.006 [τ]
Total time: 6000 [τ]
Strain rate: 0.0001 [1/τ]
Maximum strain: 0.6
Poisson's ratio: 0
図１７（Ａ）は、比較例のゴム材料モデル（歪：０．０）を示す図である。図１７（Ｂ）は、比較例のゴム材料モデル（歪：０．１２）を示す図である。図１７（Ｃ）は、比較例のゴム材料モデル（歪：０．３０）を示す図である。図１７（Ｄ）は、比較例のゴム材料モデル（歪：０．６０）を示す図である。比較例では、主として、第１高分子鎖モデル３が配置された領域と、第２高分子鎖モデル４が配置された領域との界面で破壊が生じる現象が計算され、実際のゴム材料の破壊現象に合致しなかった。
FIG. 17A is a diagram showing a rubber material model (strain: 0.0) of a comparative example. FIG. 17B is a diagram showing a rubber material model (strain: 0.12) of a comparative example. FIG. 17C is a diagram showing a rubber material model (strain: 0.30) of a comparative example. FIG. 17D is a diagram showing a rubber material model (strain: 0.60) of a comparative example. In the comparative example, a phenomenon in which fracture occurs mainly at the interface between the region in which the first
一方、実施例では、図１６（Ａ）～（Ｄ）に示されるように、第２高分子鎖モデル４よりも小さい強度パラメータεijが定義された第１高分子鎖モデル３が配置された領域において、破壊が生じる（空孔２４が生じる）現象が計算されており、実際のゴム材料の破壊現象に合致した。
On the other hand, in the example, as shown in FIGS. 16(A) to 16(D), the first
図１８は、実施例のゴム材料モデル、及び、比較例のゴム材料モデルの応力－歪み曲線を示すグラフである。実施例は、比較例に比べて、ヤング率が大きく計算された。これは、実際のゴム材料のヤング率の傾向に近似している。 FIG. 18 is a graph showing stress-strain curves of the rubber material model of Example and the rubber material model of Comparative Example. The calculated Young's modulus of the example is larger than that of the comparative example. This approximates the trend of Young's modulus of actual rubber materials.
このように、実施例のゴム材料モデルは、変形計算を行うシミュレーションに用いられることにより、第１高分子成分と第２高分子成分とを含むゴム材料の破壊特性を、精度良く評価できることを確認できた。 Thus, it was confirmed that the rubber material model of the example can accurately evaluate the fracture characteristics of the rubber material containing the first polymer component and the second polymer component by using the simulation for performing deformation calculation. did it.
Ｓ１１ 第１高分子鎖モデルを入力する工程
Ｓ１２ 第２高分子鎖モデルを入力する工程
Ｓ２１ 第１相互作用ポテンシャルを定義する工程
Ｓ２２ 第２相互作用ポテンシャルを定義する工程
S11 Step of inputting the first polymer chain model S12 Step of inputting the second polymer chain model S21 Step of defining the first interaction potential S22 Step of defining the second interaction potential
Claims (6)
前記第１高分子成分の分子鎖構造を複数の第１粗視化粒子を用いて表現した第１高分子鎖モデルを前記コンピュータに入力する工程と、
前記第２高分子成分の分子鎖構造を複数の第２粗視化粒子を用いて表現した第２高分子鎖モデルを前記コンピュータに入力する工程と、
前記第１高分子鎖モデルと前記第２高分子鎖モデルとの間に界面が形成されるように、前記第１高分子鎖モデル及び前記第２高分子鎖モデルを、前記ゴム材料の一部に対応する仮想空間であるセルに配置したゴム材料モデルを定義する工程と、
前記第１粗視化粒子間の相互作用を定義するための第１相互作用ポテンシャルを定義する工程と、
前記第２粗視化粒子間の相互作用を定義するための第２相互作用ポテンシャルを、前記物性値の差異に基づいて、前記第１相互作用ポテンシャルとは異なるものとして定義する工程と、
前記コンピュータが、前記ゴム材料モデルの変形を計算して、前記第１相互作用ポテンシャル、及び、前記第２相互作用ポテンシャルに基づく前記界面での破壊現象を計算する工程とを含む、
ゴム材料のシミュレーション方法。 A method for simulating, using a computer, the destruction phenomenon of a rubber material containing at least a first polymer component and a second polymer component having predetermined differences in physical property values, comprising:
inputting into the computer a first polymer chain model representing the molecular chain structure of the first polymer component using a plurality of first coarse-grained particles;
inputting into the computer a second polymer chain model representing the molecular chain structure of the second polymer component using a plurality of second coarse-grained particles;
The first polymer chain model and the second polymer chain model are part of the rubber material so that an interface is formed between the first polymer chain model and the second polymer chain model. A step of defining a rubber material model placed in a cell, which is a virtual space corresponding to
defining a first interaction potential for defining interactions between the first coarse-grained grains;
defining a second interaction potential for defining an interaction between the second coarse-grained particles as being different from the first interaction potential based on the difference in the physical property value;
the computer calculating deformation of the rubber material model to calculate a fracture phenomenon at the interface based on the first interaction potential and the second interaction potential;
How to simulate rubber materials.
前記第２相互作用ポテンシャルを定義する工程は、前記第２相互作用ポテンシャルに、前記第１相互作用ポテンシャルとは異なる前記強度パラメータを定義する工程を含む、請求項１記載のゴム材料のシミュレーション方法。 said first interaction potential and said second interaction potential comprising a strength parameter for defining the strength of the interaction;
2. The method of simulating a rubber material according to claim 1, wherein the step of defining said second interaction potential comprises defining said strength parameter for said second interaction potential that is different from said first interaction potential.
ｒij：粒子間の距離
ｒij,cut：カットオフ距離
εij：強度パラメータ
σij：粒子の直径に相当 3. The rubber material simulation method according to claim 2, wherein said first interaction potential and said second interaction potential are defined by the following equations.
r ij : Distance between particles r ij , cut : Cutoff distance ε ij : Strength parameter σ ij : Corresponding to particle diameter
前記物性値の差異に基づいて、前記複数種類の前記強度パラメータεijから、前記第２相互作用ポテンシャルの前記強度パラメータεijを選択する工程とを含む、請求項３記載のゴム材料のシミュレーション方法。 The step of defining the strength parameter is a step of calculating the physical property value for each of a plurality of predetermined types of the strength parameter εij ;
4. The rubber material simulation method according to claim 3, further comprising the step of selecting said strength parameter ε ij of said second interaction potential from said plurality of types of said strength parameter ε ij based on the difference in said physical property value. .
ゴム材料の製造方法。 6. The rubber material simulation method according to any one of claims 1 to 5, comprising the step of manufacturing the rubber material based on the rubber material model evaluated as having good performance,
A method of manufacturing a rubber material.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2018219327A JP7159809B2 (en) | 2018-11-22 | 2018-11-22 | Rubber material simulation method and rubber material manufacturing method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2018219327A JP7159809B2 (en) | 2018-11-22 | 2018-11-22 | Rubber material simulation method and rubber material manufacturing method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2020086835A JP2020086835A (en) | 2020-06-04 |
| JP7159809B2 true JP7159809B2 (en) | 2022-10-25 |
Family
ID=70908217
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2018219327A Active JP7159809B2 (en) | 2018-11-22 | 2018-11-22 | Rubber material simulation method and rubber material manufacturing method |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP7159809B2 (en) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003318363A (en) | 2002-04-25 | 2003-11-07 | Sharp Corp | Salient electrode bonded type semiconductor device and method of manufacturing the same |
| JP6254325B1 (en) | 2017-08-21 | 2017-12-27 | 住友ゴム工業株式会社 | Coarse-grained molecular dynamics simulation method for polymer materials |
-
2018
- 2018-11-22 JP JP2018219327A patent/JP7159809B2/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003318363A (en) | 2002-04-25 | 2003-11-07 | Sharp Corp | Salient electrode bonded type semiconductor device and method of manufacturing the same |
| JP6254325B1 (en) | 2017-08-21 | 2017-12-27 | 住友ゴム工業株式会社 | Coarse-grained molecular dynamics simulation method for polymer materials |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2020086835A (en) | 2020-06-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP5923069B2 (en) | Method for simulating polymer materials | |
| US9990450B2 (en) | Simulation method for high polymer material | |
| JP6408856B2 (en) | Method for simulating polymer materials | |
| JP6254325B1 (en) | Coarse-grained molecular dynamics simulation method for polymer materials | |
| JP2005121535A (en) | Method for simulating rubber material | |
| JP5530480B2 (en) | Method for simulating polymer materials | |
| JP6776876B2 (en) | Simulation method for polymer materials | |
| JP5592921B2 (en) | Method for simulating polymer materials | |
| JP6414929B2 (en) | How to create an all-atom model | |
| JP5001883B2 (en) | Rubber material simulation method | |
| JP7159809B2 (en) | Rubber material simulation method and rubber material manufacturing method | |
| JP4697870B2 (en) | Simulation method for viscoelastic materials | |
| JP6794820B2 (en) | How to create a polymer material model | |
| JP5749973B2 (en) | Rubber material simulation method | |
| JP6554995B2 (en) | Method for simulating polymer materials | |
| JP6575062B2 (en) | Method for simulating polymer materials | |
| JP6101159B2 (en) | Calculation method of energy loss of polymer materials | |
| JP6962160B2 (en) | Coarse-grained molecular dynamics simulation method for polymer materials | |
| JP6711186B2 (en) | Polymer material simulation method | |
| JP2020085692A (en) | Rubber material simulation method | |
| JP7290037B2 (en) | Rubber material simulation method and rubber material manufacturing method | |
| JP6408874B2 (en) | Method for simulating modified polymer materials | |
| JP7107111B2 (en) | Method for determining interaction potential of polymeric material, method for creating polymeric material model, method for simulating polymeric material, and method for manufacturing polymeric material | |
| JP2024007219A (en) | Polymer material analyzing method | |
| JP2022139140A (en) | Filler model creation method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20210928 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20220817 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20220823 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20220829 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20220913 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20220926 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 7159809 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |