JP2007031541A - 高分子材料設計方法および装置 - Google Patents
高分子材料設計方法および装置 Download PDFInfo
- Publication number
- JP2007031541A JP2007031541A JP2005215695A JP2005215695A JP2007031541A JP 2007031541 A JP2007031541 A JP 2007031541A JP 2005215695 A JP2005215695 A JP 2005215695A JP 2005215695 A JP2005215695 A JP 2005215695A JP 2007031541 A JP2007031541 A JP 2007031541A
- Authority
- JP
- Japan
- Prior art keywords
- polymer
- physical property
- property value
- model
- design
- 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.)
- Granted
Links
Images
Landscapes
- Polymerisation Methods In General (AREA)
Abstract
【解決手段】 この高分子材料設計方法は、計算機を用いて所望の物性値条件を満たす高分子の構造を求めるものであって、高分子モデルの作成条件および物性値条件を入力するステップと、この作成条件の下でモノマーに基づく構成単位を高分子の設計単位とし、この設計単位を遺伝子として遺伝的アルゴリズムの個体となる高分子モデルを作成するステップと、作成された高分子モデルについて物性値を計算し、前記物性値条件に対する適合度を計算するステップとを有する。そして、遺伝的アルゴリズムを用いて、次世代以降の個体の作成および適合度計算のステップを繰り返し、適合度の大きい高分子の構造を探索することを特徴とする。
【選択図】 図1
Description
まず、以下に、高分子材料設計装置について図1を用いて説明する。
高分子材料設計装置は、高分子モデルの作成条件や物性値条件等を入力する入力装置10と、入力された高分子モデルの作成条件に基づいて、設計単位を構成単位とする高分子モデルを作成する高分子モデル作成装置20と、高分子モデルについて物性値を計算する物性値計算装置30と、物性値から適合度を計算し、遺伝的アルゴリズムを用いて適合度の大きい高分子の構造を探索する高分子構造探索装置40と、および高分子構造探索装置で探索した結果を表示する表示装置50から構成されている。
なお、これらの各装置やそれを構成する各部は、コンピュータがこれらの機能を行なう部分を仮想的に表すものであって、ハードウエアとしての具体的な構成が限定されるものではない。
まず、ステップ(1a)において、高分子モデルの作成条件を入力する。高分子モデルの作成条件としては、この実施の形態では、高分子モデルの原型、高分子モデル当たりに配列される設計単位の数、遺伝子として用いる設計単位の種類が挙げられる。
高分子モデルは、設計単位を構成単位として作成され、設計単位を複数個つなげることで高分子鎖を作成する。図3に示すように、原型として種々の骨組みが考えられるが、ここでは、2次元以上とする実益が無いので、1次元つまり直線構造を採用した。
高分子材料に求められる物性値を選択し、物性値条件を表現する適合度評価関数を定義し、入力する。物性値としては、密度、凝集エネルギー密度、溶解度パラメータ、ガラス転移温度、融点、体積弾性率、体積膨張率、定圧比熱、定積比熱、屈折率、固有複屈折率、配向複屈折率等が挙げられる。通常、これらの中から目的に応じた1つの物性値(X)を選択するが、場合によっては、複数の物性値(X,Y,・・)を選択する。そして、選択した物性値条件を表現する適合度評価関数F=f(X)を定義し、入力する。
F=X
と定義する。物性値Xが特定の値Aに近い程良いのであれば、
F=1/(X−A)2
と定義する。
物性値が複数の場合には、適合度評価関数Fもより複雑になる。例えば、
F=XaYb
F=1/{a(X−A)2+b(Y−B)2}
等となる。これらの式で、a,bは、複数の物性値X、Yに重み付けを与えるための定数である。
さらに、Head to Tail型構造やHead to Head型構造も考えるならば、選択する設計単位としては、エチレン、プロピレン−(1)、1−ブテン−(1)、1−ペンテン(1)、1−ヘキセン−(1)、プロピレン−(2)、1−ブテン−(2)、1−ペンテン(2)、1−ヘキセン−(2)、プロピレン−(3)、1−ブテン−(3)、1−ペンテン(3)、1−ヘキセン−(3)、プロピレン−(4)、1−ブテン−(4)、1−ペンテン(4)、および1−ヘキセン−(4)の17種類の設計単位が必要である。
高分子モデルの作成方法として、配列を使う。配列S(i)が高分子鎖であり、配列要素iが設計単位である。
図3にあるような直鎖型の高分子鎖ならば、1次元配列S(i)を使い、ステップ(2a)で選択した設計単位を1次元配列S(1),S(2),…,S(i),…,S(n)に入れる。これをN回繰り返し、図5のように、高分子モデルをN種類作成する。
ステップ(2b)で作成した高分子モデルを記述した構造ファイルを作成する。構造ファイルは、高分子モデルを構成する各原子の元素名、各原子の位置を表す3次元座標、および原子の結合のデータ等を有する。構造ファイルのフォーマットとして、図7に示したようなMDL−Mol形式をはじめとする汎用のフォーマットが挙げられる。構造ファイルのフォーマットの選び方としては、ステップ(3b)で使用する各汎用ログラムで読み込めるフォーマットでなければならない。
物性値としては、密度、凝集エネルギー密度、溶解度パラメータ、ガラス転移温度、融点、体積弾性率、体積膨張率、定圧比熱、定積比熱、屈折率、固有複屈折率、配向複屈折率等が挙げられる。これら物性値は、分子動力学計算、あるいは分子動力学計算に分子力学計算、分子軌道計算、モンテカルロ計算およびバンド計算のうち少なくとも1つを組み合わせて求められる。これらの計算は、例えば、財団法人化学技術戦略推進機構のOCTA、株式会社日本総合研究所のJ−OCTA、Accelrys Inc.のMS Modeling products等(いずれも商品名)の汎用プログラムを使って行うことができる。
ルーレット選択方式は、適合度評価関数の大きさに比例した割合で高分子モデルを選択する方式である。適合度評価関数の大きさに比例した面積でN個の領域に分割されているルーレットを作り、投げた玉の入った領域の高分子モデルを選択するという過程をM回繰り返して、M種類の高分子モデルを選択する方式である。
トーナメント選択方式は、高分子モデルをある数(トーナメントサイズ)だけランダムに選び出し、その中で適合度評価関数が最大の高分子モデルを選択するという過程をM回繰り返して、M種類の高分子モデルを選択する方式である。
エリート選択方式は、適合度評価関数が大きい高分子モデルを上位からM種類選択する方式である。
図9に例を示す。
20 高分子モデル作成装置
21 高分子モデル作成部
22 入力部
23 乱数発生部
24 設計単位データベース
30 物性値計算装置
31 データ変換部
32 物性値計算部
40 高分子構造探索装置
41 適合度計算部
42 高分子モデル選択部
43 交叉発生部
50 表示装置
Claims (6)
- 計算機を用いて、所望の物性値条件を満たす高分子の構造を求める高分子材料設計方法であって、
高分子モデルの作成条件および物性値条件を入力するステップと、
前記作成条件の下でモノマーに基づく構成単位を高分子の設計単位とし、該設計単位を遺伝子として遺伝的アルゴリズムの個体となる高分子モデルを作成するステップと、
作成された高分子モデルについて物性値を計算し、前記物性値条件に対する適合度を計算するステップとを有し、
遺伝的アルゴリズムを用いて、次世代以降の個体の作成および適合度計算のステップを繰り返し、適合度の大きい高分子の構造を探索することを特徴とする高分子材料設計方法。 - 前記高分子モデルの作成条件に、高分子モデルの原型の構造、高分子モデル当たりに配列される設計単位の数、高分子モデルを構成する設計単位の種類、の少なくとも1つを含むことを特徴とする請求項1に記載の高分子材料設計方法。
- 立体配置上の異性体を有する構成単位を、該異性体毎に異なる設計単位として高分子モデルを作成することを特徴とする請求項1または請求項2に記載の高分子材料設計方法。
- 前記遺伝的アルゴリズムにおける適合度は、所望の物性値条件を満たすほど大きな値を示す評価関数として定義されることを特徴とする請求項1ないし請求項3のいずれかに記載の高分子材料設計方法。
- 前記物性値は、前記高分子モデルについて、分子動力学計算、あるいは分子動力学計算に分子力学計算、分子軌道計算、モンテカルロ計算およびバンド計算のうち少なくとも1つまたは複数を用いて求められることを特徴とする請求項1ないし請求項4のいずれかに記載の高分子材料設計方法。
- 計算機を用いて、所望の物性値条件を満たす高分子の構造を求める高分子材料設計装置において、
高分子モデルの作成条件および物性値条件等を入力する装置と、
前記作成条件の下で設計単位を遺伝子として遺伝的アルゴリズムの個体となる高分子モデルを作成する装置と、
作成された高分子モデルについて物性値を計算する装置と、
物性値から適合度を計算し、遺伝的アルゴリズムを用いて適合度の大きい高分子の構造を探索する装置と、
探索結果を表示する装置と
を備えることを特徴とする高分子材料設計装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2005215695A JP4933755B2 (ja) | 2005-07-26 | 2005-07-26 | 高分子材料設計方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2005215695A JP4933755B2 (ja) | 2005-07-26 | 2005-07-26 | 高分子材料設計方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2007031541A true JP2007031541A (ja) | 2007-02-08 |
JP4933755B2 JP4933755B2 (ja) | 2012-05-16 |
Family
ID=37791173
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2005215695A Expired - Fee Related JP4933755B2 (ja) | 2005-07-26 | 2005-07-26 | 高分子材料設計方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP4933755B2 (ja) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008282389A (ja) * | 2007-03-30 | 2008-11-20 | Honda Research Inst Europe Gmbh | 間接的な設計コーディングを用いた設計の最適化 |
JP2010211761A (ja) * | 2009-03-12 | 2010-09-24 | Japan Agengy For Marine-Earth Science & Technology | 鎖状構造体情報の生成装置、生成方法、及び生成プログラム |
JP2023026082A (ja) * | 2021-08-12 | 2023-02-24 | 株式会社三菱総合研究所 | 情報処理装置、情報処理方法及びプログラム |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11281642A (ja) * | 1998-03-27 | 1999-10-15 | Toray Ind Inc | 分子集合体モデルの屈折率算出方法および装置 |
JP2003105090A (ja) * | 2001-09-28 | 2003-04-09 | Asahi Kasei Corp | 高分子材料の形態予測方法およびその装置 |
JP2004287812A (ja) * | 2003-03-20 | 2004-10-14 | Japan Research Institute Ltd | 物性値算出方法、物性値算出装置、画像作成装置、コンピュータプログラム及び記録媒体 |
JP2006282929A (ja) * | 2005-04-04 | 2006-10-19 | Taiyo Nippon Sanso Corp | 分子構造予測方法 |
-
2005
- 2005-07-26 JP JP2005215695A patent/JP4933755B2/ja not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11281642A (ja) * | 1998-03-27 | 1999-10-15 | Toray Ind Inc | 分子集合体モデルの屈折率算出方法および装置 |
JP2003105090A (ja) * | 2001-09-28 | 2003-04-09 | Asahi Kasei Corp | 高分子材料の形態予測方法およびその装置 |
JP2004287812A (ja) * | 2003-03-20 | 2004-10-14 | Japan Research Institute Ltd | 物性値算出方法、物性値算出装置、画像作成装置、コンピュータプログラム及び記録媒体 |
JP2006282929A (ja) * | 2005-04-04 | 2006-10-19 | Taiyo Nippon Sanso Corp | 分子構造予測方法 |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008282389A (ja) * | 2007-03-30 | 2008-11-20 | Honda Research Inst Europe Gmbh | 間接的な設計コーディングを用いた設計の最適化 |
JP2010211761A (ja) * | 2009-03-12 | 2010-09-24 | Japan Agengy For Marine-Earth Science & Technology | 鎖状構造体情報の生成装置、生成方法、及び生成プログラム |
JP2023026082A (ja) * | 2021-08-12 | 2023-02-24 | 株式会社三菱総合研究所 | 情報処理装置、情報処理方法及びプログラム |
JP7270005B2 (ja) | 2021-08-12 | 2023-05-09 | 株式会社三菱総合研究所 | 情報処理装置、情報処理方法及びプログラム |
Also Published As
Publication number | Publication date |
---|---|
JP4933755B2 (ja) | 2012-05-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Das et al. | Computational linear rheology of general branch-on-branch polymers | |
Geem | Optimal cost design of water distribution networks using harmony search | |
Lee et al. | A new meta-heuristic algorithm for continuous engineering optimization: harmony search theory and practice | |
Lee et al. | The harmony search heuristic algorithm for discrete structural optimization | |
CN106023015B (zh) | 课程学习路径推荐方法及装置 | |
Cruz et al. | Efficient differential evolution algorithms for multimodal optimal control problems | |
Cohen | Designing test suites for software interaction testing | |
Romero et al. | A modified genetic algorithm for reservoir characterisation | |
Paz et al. | Shapes of clusters and groups of galaxies: comparison of model predictions with observations | |
Chen et al. | Computer-aided design method for the components made of heterogeneous materials | |
US20070239411A1 (en) | Optimal Design Support System, Optical Design Support Method and Optimal Design Support Program | |
Read | From reactor to rheology in industrial polymers | |
JP4933755B2 (ja) | 高分子材料設計方法 | |
CN108563875A (zh) | 基于多目标优化的模拟电路测点和频率联合优选方法 | |
Wan et al. | On graph entropy measures based on the number of independent sets and matchings | |
Došlić et al. | Enumerative aspects of secondary structures | |
Zhai et al. | Automated selection of gate location for plastic injection molding processing | |
Sun et al. | A temperature field reconstruction method for spacecraft leading edge structure with optimized sensor array | |
Hallam et al. | Large barrier trees for studying search | |
Rosenthal et al. | Impact of population size and selection within a customized NSGA-II for biochemical optimization assessed on the basis of the average cuboid volume indicator | |
Ye et al. | Proposed strategy for the application of the modified harmony search algorithm to code-based selection and scaling of ground motions | |
Alvarado-Contreras et al. | Constitutive modeling of damage evolution in semicrystalline polyethylene | |
Hsu | Monte Carlo simulations of lattice models for single polymer systems | |
Zara et al. | Exact solution of a RNA-like polymer model on the Husimi lattice | |
JP7347148B2 (ja) | 高分子材料のシミュレーション方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20080724 |
|
RD05 | Notification of revocation of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7425 Effective date: 20090311 |
|
RD03 | Notification of appointment of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7423 Effective date: 20090316 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20101110 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20101116 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20110117 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20111206 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20120118 |
|
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: 20120214 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20120217 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20150224 Year of fee payment: 3 |
|
LAPS | Cancellation because of no payment of annual fees |