JP6097130B2 - 高分子材料のシミュレーション方法 - Google Patents
高分子材料のシミュレーション方法 Download PDFInfo
- Publication number
- JP6097130B2 JP6097130B2 JP2013085077A JP2013085077A JP6097130B2 JP 6097130 B2 JP6097130 B2 JP 6097130B2 JP 2013085077 A JP2013085077 A JP 2013085077A JP 2013085077 A JP2013085077 A JP 2013085077A JP 6097130 B2 JP6097130 B2 JP 6097130B2
- Authority
- JP
- Japan
- Prior art keywords
- model
- coarse
- parameter
- atom
- length
- 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
- 238000000034 method Methods 0.000 title claims description 72
- 239000002861 polymer material Substances 0.000 title claims description 38
- 238000004364 calculation method Methods 0.000 claims description 113
- 238000000329 molecular dynamics simulation Methods 0.000 claims description 66
- 229920000642 polymer Polymers 0.000 claims description 58
- 238000004088 simulation Methods 0.000 claims description 53
- 239000011324 bead Substances 0.000 claims description 45
- 125000004432 carbon atom Chemical group C* 0.000 claims description 26
- 125000004429 atom Chemical group 0.000 claims description 21
- 238000006243 chemical reaction Methods 0.000 claims description 18
- 229910052799 carbon Inorganic materials 0.000 claims description 17
- 238000005307 time correlation function Methods 0.000 claims description 14
- 239000000463 material Substances 0.000 claims description 4
- 230000000087 stabilizing effect Effects 0.000 claims description 3
- 239000002245 particle Substances 0.000 description 14
- 230000008569 process Effects 0.000 description 11
- 238000012545 processing Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 7
- 238000002474 experimental method Methods 0.000 description 7
- 238000000324 molecular mechanic Methods 0.000 description 7
- 230000000052 comparative effect Effects 0.000 description 5
- 239000005062 Polybutadiene Substances 0.000 description 4
- 125000004435 hydrogen atom Chemical group [H]* 0.000 description 4
- 229920002857 polybutadiene Polymers 0.000 description 4
- 238000000342 Monte Carlo simulation Methods 0.000 description 2
- 229920003193 cis-1,4-polybutadiene polymer Polymers 0.000 description 2
- 229920001971 elastomer Polymers 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000005452 bending Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 235000020057 cognac Nutrition 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000000806 elastomer Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000000155 melt Substances 0.000 description 1
- 238000006116 polymerization reaction Methods 0.000 description 1
- 238000005381 potential energy Methods 0.000 description 1
- 239000011347 resin Substances 0.000 description 1
- 229920005989 resin Polymers 0.000 description 1
- 239000005060 rubber Substances 0.000 description 1
- 238000010058 rubber compounding Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000003936 working memory Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C10/00—Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computing Systems (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Description
ここで、
<Rc 2>:粗視化モデルの末端間ベクトルの長さの2乗のアンサンブル平均
Lc:粗視化モデルの全長
Ac:ビーズ間の平衡長
Nc:1つの粗視化モデルあたりのRouseビーズの個数
ここで、
bc:粗視化モデルの長さのパラメータ
τc:粗視化モデルの時間のパラメータ
Nc:1つの粗視化モデルあたりのRouseビーズの個数
kB:ボルツマン定数
T:絶対温度
ここで、
<Rf 2>:全原子モデルの末端間ベクトルの長さの2乗のアンサンブル平均
Lf:全原子モデルの全長
ここで、
bf:全原子モデルの長さのパラメータ
τf:全原子モデルの時間のパラメータ
Nf:1つの全原子モデルあたりのRouseビーズの個数
Lf:全原子モデルの全長
kB:ボルツマン定数
T:絶対温度
図1に示されるように、本実施形態の高分子材料のシミュレーション方法(以下、単に「シミュレーション方法」ということがある)は、高分子材料の反応を、コンピュータを用いて評価するためのものである。高分子材料としては、例えば、ゴム、樹脂又はエラストマー等が含まれる。
ここで、各定数及び変数は、次のとおりである。
rij:各ビーズ間距離
rc:カットオフ距離
ε:各ビーズ間に定義される相互ポテンシャルの強度
σ:各ビーズ間に定義される相互ポテンシャルの距離に関するパラメータ(分子動力学の分野では、LJ球の直径と呼ばれる)
なお、ビーズ間の距離rij及びカットオフ距離rcは、各ビーズ3の中心3a(図3に示す)を基準として定義される。
強度ε:1.0
距離σ:1.0
カットオフ距離rc:21/6
ここで、
<Rc 2>:粗視化モデルの末端間ベクトルの長さの2乗のアンサンブル平均
Lc:粗視化モデルの全長
Ac:ビーズ間の平衡長
Nc:1つの粗視化モデルあたりのRouseビーズの個数
Cc=A×exp(−t/τc)…(6)
ここで、
Cc:粗視化モデルの末端間ベクトルの時間相関関数
A:フィッティングパラメータ
t:時間
τc:粗視化モデルの時間のパラメータ
ここで、
bc:粗視化モデルの長さのパラメータ
τc:粗視化モデルの時間のパラメータ
Nc:1つの粗視化モデルあたりのRouseビーズの個数
kB:ボルツマン定数
T:絶対温度
ここで、
<Rf 2>:全原子モデルの末端間ベクトルの長さの2乗のアンサンブル平均
Lf:全原子モデルの全長
Cf=A×exp(−t/τf)…(6)
ここで、
Cf:全原子モデルの末端間ベクトルの長さの時間相関関数
A:フィッティングパラメータ
t:時間
τf:粗視化モデルの時間のパラメータ
ここで、
bf:全原子モデルの長さのパラメータ
τf:全原子モデルの時間のパラメータ
Nf:1つの全原子モデルあたりのRouseビーズの個数
Lf:全原子モデルの全長
kB:ボルツマン定数
T:絶対温度
高分子鎖:
構造:cis-1,4ポリブタジエン
重合度:2070
仮想空間:
1辺の長さLs:18σ
6 粗視化モデル
11 全原子モデル
Claims (8)
- コンピュータに、高分子鎖を、複数のビーズでモデル化した粗視化モデルを入力する工程、
前記コンピュータが、予め定めた仮想空間に配置された前記粗視化モデルを用いて分子動力学計算を行う工程、
前記コンピュータが、前記粗視化モデルのRouseのパラメータを求める第1計算工程、
前記コンピュータに、任意の高分子鎖を、複数の炭素を含む原子でモデル化した全原子モデルを入力する工程、
前記コンピュータが、予め定めた仮想空間に配置された前記全原子モデルを用いて分子動力学計算を行う工程、
前記コンピュータが、前記全原子モデルのRouseのパラメータを求める第2計算工程、及び、
前記コンピュータが、前記粗視化モデルのRouseのパラメータと、前記全原子モデルのRouseのパラメータとに基づいて、前記粗視化モデルを用いた分子動力学計算での単位系を、前記任意の高分子鎖の単位系に換算する換算工程を含み、
前記第2計算工程は、前記全原子モデルの全長Lfを求める全長計算工程を含み、
前記全長計算工程は、前記コンピュータが、前記全原子モデルの前記炭素原子の一端の固定する工程と、
前記コンピュータが、分子動力学計算を行いながら、前記全原子モデルの前記炭素原子の他端を、前記一端から離れる方向へ強制移動させる工程と、
前記コンピュータが、前記分子動力学計算によって強制的に引き伸ばされた前記全原子モデルを分子力学法による計算を行うことにより、前記全原子モデルの構造を安定化する工程と、
前記コンピュータが、前記全原子モデルの前記一端と前記他端との間の距離から前記全原子モデルの全長Lfを求める工程とを含むことを特徴とする高分子材料のシミュレーション方法。 - 前記粗視化モデルのRouseのパラメータは、長さのパラメータbcと、時間のパラメータτcと、摩擦のパラメータζcとを含み、
前記全原子モデルのRouseのパラメータは、長さのパラメータbfと、時間のパラメータτfと、摩擦のパラメータζfとを含み、
前記換算工程は、前記粗視化モデルを用いた分子動力学計算での単位長さ、単位時間、又は、単位質量を、前記任意の高分子鎖の長さ、時間、又は、質量にそれぞれ換算する請求項1に記載の高分子材料のシミュレーション方法。 - 前記粗視化モデルの前記長さのパラメータbcは、下記式(1)で定義される請求項2に記載の高分子材料のシミュレーション方法。
ここで、
<Rc 2>:粗視化モデルの末端間ベクトルの長さの2乗のアンサンブル平均
Lc:粗視化モデルの全長
Ac:ビーズ間の平衡長
Nc:1つの粗視化モデルあたりのRouseビーズの個数 - 前記粗視化モデルの前記時間のパラメータτcは、前記粗視化モデルの末端間ベクトルの時間相関関数が、1/eに減少するのに要する時間である請求項2又は3に記載の高分子材料のシミュレーション方法。
- 前記粗視化モデルの前記摩擦のパラメータζcは、下記式(2)で定義される請求項2乃至4のいずれかに記載の高分子材料のシミュレーション方法。
ここで、
bc:粗視化モデルの長さのパラメータ
τc:粗視化モデルの時間のパラメータ
Nc:1つの粗視化モデルあたりのRouseビーズの個数
kB:ボルツマン定数
T:絶対温度 - 前記全原子モデルの前記長さのパラメータbfは、下記式(3)で定義される請求項2乃至5のいずれかに記載の高分子材料のシミュレーション方法。
ここで、
<Rf 2>:全原子モデルの末端間ベクトルの長さの2乗のアンサンブル平均
Lf:全原子モデルの全長 - 前記時間のパラメータτfは、前記全原子モデルの末端間ベクトルの時間相関関数が、1/eに減少するのに要する時間である請求項2乃至6のいずれかに記載の高分子材料のシミュレーション方法。
- 前記摩擦のパラメータζfは、下記式(4)で定義される請求項2乃至7のいずれかに記載の高分子材料のシミュレーション方法。
ここで、
bf:全原子モデルの長さのパラメータ
τf:全原子モデルの時間のパラメータ
Nf:1つの全原子モデルあたりのRouseビーズの個数
Lf:全原子モデルの全長
kB:ボルツマン定数
T:絶対温度
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2013085077A JP6097130B2 (ja) | 2013-04-15 | 2013-04-15 | 高分子材料のシミュレーション方法 |
US14/219,699 US9864845B2 (en) | 2013-04-15 | 2014-03-19 | Simulation method for macromolecular material |
EP20140161985 EP2793150A1 (en) | 2013-04-15 | 2014-03-27 | Simulation method for macromolecular material |
KR1020140044485A KR102359411B1 (ko) | 2013-04-15 | 2014-04-14 | 고분자 재료의 시뮬레이션 방법 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2013085077A JP6097130B2 (ja) | 2013-04-15 | 2013-04-15 | 高分子材料のシミュレーション方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2014206913A JP2014206913A (ja) | 2014-10-30 |
JP6097130B2 true JP6097130B2 (ja) | 2017-03-15 |
Family
ID=50391039
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2013085077A Active JP6097130B2 (ja) | 2013-04-15 | 2013-04-15 | 高分子材料のシミュレーション方法 |
Country Status (4)
Country | Link |
---|---|
US (1) | US9864845B2 (ja) |
EP (1) | EP2793150A1 (ja) |
JP (1) | JP6097130B2 (ja) |
KR (1) | KR102359411B1 (ja) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10019827B2 (en) | 2013-12-31 | 2018-07-10 | Disney Enterprises, Inc. | Material point method for simulation of granular materials |
US10210287B2 (en) * | 2013-12-31 | 2019-02-19 | Disney Enterprises, Inc. | Augmented material point method for simulating phase changes and varied materials |
CN104915473B (zh) * | 2015-05-15 | 2019-04-30 | 华东理工大学 | 一种测量晶体表面结构之间距离的方法 |
JP6050903B1 (ja) * | 2016-01-22 | 2016-12-21 | 住友ゴム工業株式会社 | 高分子材料のシミュレーション方法 |
JP6651254B2 (ja) * | 2016-03-29 | 2020-02-19 | 住友重機械工業株式会社 | シミュレーション方法、シミュレーションプログラム、及びシミュレーション装置 |
JP6711186B2 (ja) * | 2016-07-14 | 2020-06-17 | 住友ゴム工業株式会社 | 高分子材料のシミュレーション方法 |
EP3447717A1 (en) * | 2017-08-24 | 2019-02-27 | Tata Consultancy Services Limited | Systems and methods for determining properties of composite materials for predicting behaviour of structures |
JP6844930B2 (ja) * | 2017-09-05 | 2021-03-17 | 住友重機械工業株式会社 | シミュレーション方法及びシミュレーション装置 |
CN110442887B (zh) * | 2018-05-03 | 2021-04-06 | 中国科学院化学研究所 | 一种研究重复聚酰亚胺体系热力学性质的模型构建方法 |
CN110070918B (zh) * | 2019-04-02 | 2022-12-27 | 重庆邮电大学 | 基于分子间相互作用的粗粒化方法 |
CN114974446A (zh) * | 2022-06-24 | 2022-08-30 | 中国科学院兰州化学物理研究所 | 一种基于分子动力学模拟的Stribeck曲线计算方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012177609A (ja) * | 2011-02-25 | 2012-09-13 | Bridgestone Corp | 高分子材料のモデル作成方法 |
-
2013
- 2013-04-15 JP JP2013085077A patent/JP6097130B2/ja active Active
-
2014
- 2014-03-19 US US14/219,699 patent/US9864845B2/en active Active
- 2014-03-27 EP EP20140161985 patent/EP2793150A1/en not_active Ceased
- 2014-04-14 KR KR1020140044485A patent/KR102359411B1/ko active IP Right Grant
Also Published As
Publication number | Publication date |
---|---|
KR20140123915A (ko) | 2014-10-23 |
KR102359411B1 (ko) | 2022-02-07 |
JP2014206913A (ja) | 2014-10-30 |
EP2793150A1 (en) | 2014-10-22 |
US20140309971A1 (en) | 2014-10-16 |
US9864845B2 (en) | 2018-01-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6097130B2 (ja) | 高分子材料のシミュレーション方法 | |
JP6085224B2 (ja) | フィラー間の相互作用ポテンシャルの計算方法 | |
JP5923069B2 (ja) | 高分子材料のシミュレーション方法 | |
JP6408856B2 (ja) | 高分子材料のシミュレーション方法 | |
JP6254325B1 (ja) | 高分子材料の粗視化分子動力学シミュレーション方法 | |
JP6353290B2 (ja) | 高分子材料モデル作成方法 | |
JP6082303B2 (ja) | 高分子材料のシミュレーション方法 | |
JP6414929B2 (ja) | 全原子モデルの作成方法 | |
JP2016081183A (ja) | ゴム材料のシミュレーション方法 | |
JP7040152B2 (ja) | 高分子材料のシミュレーション方法 | |
JP5592921B2 (ja) | 高分子材料のシミュレーション方法 | |
JP6200193B2 (ja) | 高分子材料のシミュレーション方法 | |
JP6554995B2 (ja) | 高分子材料のシミュレーション方法 | |
JP6055359B2 (ja) | 高分子材料のシミュレーション方法 | |
JP6575062B2 (ja) | 高分子材料のシミュレーション方法 | |
JP6101159B2 (ja) | 高分子材料のエネルギーロスの計算方法 | |
JP2018112525A (ja) | 高分子材料のシミュレーション方法 | |
JP6711186B2 (ja) | 高分子材料のシミュレーション方法 | |
JP2017182231A5 (ja) | ||
JP6050903B1 (ja) | 高分子材料のシミュレーション方法 | |
JP2015102972A (ja) | 粒子集団の座標を定義する方法 | |
JP2016070889A (ja) | 高分子材料のシミュレーション方法 | |
JP6434805B2 (ja) | 高分子材料のシミュレーション方法 | |
JP6962160B2 (ja) | 高分子材料の粗視化分子動力学シミュレーション方法 | |
JP7159809B2 (ja) | ゴム材料のシミュレーション方法、及び、ゴム材料の製造方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20160122 |
|
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: 20170207 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20170217 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6097130 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |