JP6963119B2 - 熱力学的平衡状態の予測装置、予測方法、及び予測プログラム - Google Patents
熱力学的平衡状態の予測装置、予測方法、及び予測プログラム Download PDFInfo
- Publication number
- JP6963119B2 JP6963119B2 JP2020553834A JP2020553834A JP6963119B2 JP 6963119 B2 JP6963119 B2 JP 6963119B2 JP 2020553834 A JP2020553834 A JP 2020553834A JP 2020553834 A JP2020553834 A JP 2020553834A JP 6963119 B2 JP6963119 B2 JP 6963119B2
- Authority
- JP
- Japan
- Prior art keywords
- model
- prediction
- learning
- thermodynamic equilibrium
- input
- 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 25
- 239000013077 target material Substances 0.000 claims description 35
- 238000013461 design Methods 0.000 claims description 30
- 238000010586 diagram Methods 0.000 claims description 27
- 239000000203 mixture Substances 0.000 claims description 24
- 230000006870 function Effects 0.000 claims description 23
- 238000012549 training Methods 0.000 claims description 21
- 238000010801 machine learning Methods 0.000 claims description 19
- 229910000838 Al alloy Inorganic materials 0.000 claims description 13
- 238000004519 manufacturing process Methods 0.000 claims description 12
- 238000013528 artificial neural network Methods 0.000 claims description 10
- 210000002569 neuron Anatomy 0.000 claims description 10
- 238000013459 approach Methods 0.000 claims description 9
- 239000000463 material Substances 0.000 claims description 7
- 238000013135 deep learning Methods 0.000 claims description 6
- 238000013277 forecasting method Methods 0.000 claims description 2
- 238000010587 phase diagram Methods 0.000 description 11
- 238000004364 calculation method Methods 0.000 description 10
- 150000001875 compounds Chemical class 0.000 description 10
- 239000000956 alloy Substances 0.000 description 8
- 229910002059 quaternary alloy Inorganic materials 0.000 description 8
- 238000004088 simulation Methods 0.000 description 7
- 229910045601 alloy Inorganic materials 0.000 description 5
- 238000004891 communication Methods 0.000 description 5
- 229910052802 copper Inorganic materials 0.000 description 5
- 229910052749 magnesium Inorganic materials 0.000 description 5
- 229910052710 silicon Inorganic materials 0.000 description 5
- 239000000654 additive Substances 0.000 description 4
- 230000000996 additive effect Effects 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 239000007864 aqueous solution Substances 0.000 description 2
- 239000000919 ceramic Substances 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000000724 energy-dispersive X-ray spectrum Methods 0.000 description 2
- 238000009434 installation Methods 0.000 description 2
- 229910000531 Co alloy Inorganic materials 0.000 description 1
- 229910000881 Cu alloy Inorganic materials 0.000 description 1
- 229910000640 Fe alloy Inorganic materials 0.000 description 1
- 229910000861 Mg alloy Inorganic materials 0.000 description 1
- 229910000914 Mn alloy Inorganic materials 0.000 description 1
- 229910000990 Ni alloy Inorganic materials 0.000 description 1
- 229910001069 Ti alloy Inorganic materials 0.000 description 1
- 238000002441 X-ray diffraction Methods 0.000 description 1
- 229910001297 Zn alloy Inorganic materials 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 229910052804 chromium Inorganic materials 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 229910052742 iron Inorganic materials 0.000 description 1
- 229910052748 manganese Inorganic materials 0.000 description 1
- 229910052759 nickel Inorganic materials 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000007637 random forest analysis Methods 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 229910052725 zinc Inorganic materials 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/30—Prediction of properties of chemical compounds, compositions or mixtures
-
- 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
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- C—CHEMISTRY; METALLURGY
- C22—METALLURGY; FERROUS OR NON-FERROUS ALLOYS; TREATMENT OF ALLOYS OR NON-FERROUS METALS
- C22C—ALLOYS
- C22C1/00—Making non-ferrous alloys
- C22C1/02—Making non-ferrous alloys by melting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C60/00—Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/26—Composites
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/70—Machine learning, data mining or chemometrics
Description
本発明は、熱力学的平衡状態の予測装置、予測方法、及び予測プログラムに関する。
合金、セラミックス、水溶液、化学反応等の熱力学的平衡状態を計算機シミュレーションにより予測する手法としてＣＡＬＰＨＡＤ（ＣＡＬｃｕｌａｔｉｏｎ ｏｆ ＰＨＡｓｅ Ｄｉａｇｒａｍ ａｎｄ Ｔｈｅｒｍｏｄｙｎａｍｉｃｓ）法が知られている（例えば特許文献１）。
ＣＡＬＰＨＡＤ法は平衡状態として最小のギブズエネルギーを与える状態を探索するシミュレーション手法であり、計算に時間がかかる。１つの材料組成に対して熱力学的平衡計算を行って状態図を予測するには、ＣＡＬＰＨＡＤ法は現実的な時間内に実施できるが、多数の材料組成に対して状態図を予測してスクリーニングを実施するには時間がかかり過ぎ、現実的ではない。例えば合金材料について１００００通りの説明変数（組成や製造条件の組）をスクリーニングすべく全合金組成の状態図を算出しようとすると、各組の計算に約９０秒かかるので、合計で約２５０時間の時間を要する。
本発明は、熱力学的平衡状態を短時間で計算できる熱力学的平衡状態の予測装置、予測方法、及び予測プログラムを提供することを目的とする。
本発明の実施態様は、以下に示す構成を備える。
［１］ 対象材料の熱力学的平衡状態の予測装置であって、
前記対象材料の設計条件に係る説明変数の入力に基づいて、前記熱力学的平衡状態に係る目的変数を出力するモデルと、
所定の設計条件の入力と、該設計条件に基づき生じ得る前記熱力学的平衡状態の出力とを含む学習用データを作成する学習用データ作成部と、
前記学習用データ作成部により作成された前記学習用データを用いて、前記モデルの入出力関係が前記学習用データの入出力関係に近づくように機械学習を行うモデル学習部と、
前記対象材料の熱力学的平衡状態の予測に用いる予測用説明変数を設定する説明変数設定部と、
前記モデル学習部による機械学習が済んだ前記モデルへの前記予測用説明変数の入力に基づき、前記熱力学的平衡状態の予測結果である予測用目的変数を前記モデルから出力する予測部と、
を備える予測装置。
前記対象材料の設計条件に係る説明変数の入力に基づいて、前記熱力学的平衡状態に係る目的変数を出力するモデルと、
所定の設計条件の入力と、該設計条件に基づき生じ得る前記熱力学的平衡状態の出力とを含む学習用データを作成する学習用データ作成部と、
前記学習用データ作成部により作成された前記学習用データを用いて、前記モデルの入出力関係が前記学習用データの入出力関係に近づくように機械学習を行うモデル学習部と、
前記対象材料の熱力学的平衡状態の予測に用いる予測用説明変数を設定する説明変数設定部と、
前記モデル学習部による機械学習が済んだ前記モデルへの前記予測用説明変数の入力に基づき、前記熱力学的平衡状態の予測結果である予測用目的変数を前記モデルから出力する予測部と、
を備える予測装置。
［２］ 前記予測部により出力された前記予測用目的変数に基づき前記熱力学的平衡状態の状態図を作成して表示する状態図表示部を備える、［１］に記載の予測装置。
［３］ 前記モデルは多層ニューラルネットワークであり、
前記モデル学習部はディープラーニングを用いて前記モデルの学習を行う、
［１］または［２］に記載の予測装置。
前記モデル学習部はディープラーニングを用いて前記モデルの学習を行う、
［１］または［２］に記載の予測装置。
［４］ 前記モデルの出力である前記目的変数は、前記対象材料の熱力学的平衡状態における相分率であり、
前記多層ニューラルネットワークの出力層のニューロンにはソフトマックス関数が用いられる、
［３］に記載の予測装置。
前記多層ニューラルネットワークの出力層のニューロンにはソフトマックス関数が用いられる、
［３］に記載の予測装置。
［５］ 前記学習用データ作成部は、前記設計条件の所定範囲の組合せから成る前記説明変数を作成し、ＣＡＬＰＨＡＤ法を用いて前記目的変数を算出して、前記作成した前記説明変数と前記算出した前記目的変数とを含む前記学習用データを作成する、
［１］〜［４］のいずれか１項に記載の予測装置。
［１］〜［４］のいずれか１項に記載の予測装置。
［６］ 前記対象材料がアルミニウム合金であり、
前記説明変数が前記アルミニウム合金の組成及び製造条件を含み、
前記目的変数が前記アルミニウム合金の熱力学的平衡状態での相分率を含む、
［１］〜［５］のいずれか１項に記載の予測装置。
前記説明変数が前記アルミニウム合金の組成及び製造条件を含み、
前記目的変数が前記アルミニウム合金の熱力学的平衡状態での相分率を含む、
［１］〜［５］のいずれか１項に記載の予測装置。
［７］
対象材料の熱力学的平衡状態の予測方法であって、
前記対象材料の設計条件に係る説明変数の入力に基づいて、前記熱力学的平衡状態に係る目的変数を出力するモデルに関して、所定の設計条件の入力と、該設計条件に基づき生じ得る前記熱力学的平衡状態の出力とを含む学習用データを作成する学習用データ作成ステップと、
前記学習用データ作成ステップにて作成された前記学習用データを用いて、前記モデルの入出力関係が前記学習用データの入出力関係に近づくように機械学習を行うモデル学習ステップと、
前記対象材料の熱力学的平衡状態の予測に用いる予測用説明変数を設定する説明変数設定ステップと、
前記説明変数設定ステップにて設定された前記予測用説明変数の、前記モデル学習ステップにて機械学習が済んだ前記モデルへの入力に基づき、前記熱力学的平衡状態の予測結果である予測用目的変数を前記モデルから出力する予測ステップと、
を含む予測方法。
対象材料の熱力学的平衡状態の予測方法であって、
前記対象材料の設計条件に係る説明変数の入力に基づいて、前記熱力学的平衡状態に係る目的変数を出力するモデルに関して、所定の設計条件の入力と、該設計条件に基づき生じ得る前記熱力学的平衡状態の出力とを含む学習用データを作成する学習用データ作成ステップと、
前記学習用データ作成ステップにて作成された前記学習用データを用いて、前記モデルの入出力関係が前記学習用データの入出力関係に近づくように機械学習を行うモデル学習ステップと、
前記対象材料の熱力学的平衡状態の予測に用いる予測用説明変数を設定する説明変数設定ステップと、
前記説明変数設定ステップにて設定された前記予測用説明変数の、前記モデル学習ステップにて機械学習が済んだ前記モデルへの入力に基づき、前記熱力学的平衡状態の予測結果である予測用目的変数を前記モデルから出力する予測ステップと、
を含む予測方法。
［８］ 対象材料の熱力学的平衡状態の予測プログラムであって、
前記対象材料の設計条件に係る説明変数の入力に基づいて、前記熱力学的平衡状態に係る目的変数を出力するモデルに関して、所定の設計条件の入力と、該設計条件に基づき生じ得る前記熱力学的平衡状態の出力とを含む学習用データを作成する学習用データ作成機能と、
前記学習用データ作成機能により作成された前記学習用データを用いて、前記モデルの入出力関係が前記学習用データの入出力関係に近づくように機械学習を行うモデル学習機能と、
前記対象材料の熱力学的平衡状態の予測に用いる予測用説明変数を設定する説明変数設定機能と、
前記説明変数設定機能により設定された前記予測用説明変数の、前記モデル学習機能により機械学習が済んだ前記モデルへの入力に基づき、前記熱力学的平衡状態の予測結果である予測用目的変数を前記モデルから出力する予測機能と、
をコンピュータに実現させるための予測プログラム。
前記対象材料の設計条件に係る説明変数の入力に基づいて、前記熱力学的平衡状態に係る目的変数を出力するモデルに関して、所定の設計条件の入力と、該設計条件に基づき生じ得る前記熱力学的平衡状態の出力とを含む学習用データを作成する学習用データ作成機能と、
前記学習用データ作成機能により作成された前記学習用データを用いて、前記モデルの入出力関係が前記学習用データの入出力関係に近づくように機械学習を行うモデル学習機能と、
前記対象材料の熱力学的平衡状態の予測に用いる予測用説明変数を設定する説明変数設定機能と、
前記説明変数設定機能により設定された前記予測用説明変数の、前記モデル学習機能により機械学習が済んだ前記モデルへの入力に基づき、前記熱力学的平衡状態の予測結果である予測用目的変数を前記モデルから出力する予測機能と、
をコンピュータに実現させるための予測プログラム。
本発明によれば、熱力学的平衡状態を短時間で計算できる熱力学的平衡状態の予測装置、予測方法、及び予測プログラムを提供することができる。
以下、添付図面を参照しながら実施形態について説明する。説明の理解を容易にするため、各図面において同一の構成要素に対しては可能な限り同一の符号を付して、重複する説明は省略する。
図１〜図５を参照して実施形態に係る熱力学的平衡状態の予測装置１（以下では単に「予測装置１」とも表記する）の構成を説明する。図１は、実施形態に係る予測装置１の概略構成を示すブロック図である。予測装置１は、複数の組成からなる材料、または、複数の製造条件の組合せにより製造される材料を含む対象材料の熱力学系平衡状態や状態図を予測するための装置である。本実施形態では、予測対象材料の一例として、アルミニウム合金を挙げて説明する。
図１に示すように、予測装置１は、モデル２と、学習用データ作成部３と、モデル学習部４と、説明変数設定部５と、予測部６と、状態図表示部７とを備える。
モデル２は、対象材料の設計条件に係る説明変数（アルミニウム合金の組成、製造条件）の入力に基づいて、熱力学的平衡状態に係る目的変数（熱力学的平衡状態における化合物の相分率）を出力する。モデル２は、教師あり学習モデルであり、熱力学的平衡状態の予測を行う前段階として、モデル学習部４によって機械学習が行われて、説明変数と目的変数との対応関係、すなわちモデル２の入出力関係が学習される。
本実施形態では、モデル２は図１に示すように入力層と、複数の中間層と、出力層とを有する多層ニューラルネットワークである。モデル２の入力層の各ニューロンは、説明変数の各項目と同数設けられ、各項目の数値が入力される。図１では、製造条件として温度Ｔ（℃）と、原料組成として添加元素３種のＳｉ、Ｃｕ、Ｍｇの重量百分率（ｗｔ％）の数値が入力層の各ニューロンに入力される。モデル２の出力層の各ニューロンは、目的変数の各項目と同数設けられ、各項目の数値が出力される。
また、モデル２の多層ニューラルネットワークの出力層にはソフトマックス関数が用いられる。つまり、モデル２の複数の出力のそれぞれが０から１の範囲であり、かつ、複数の出力の総和が１となる。上述のように本実施形態ではモデル２の出力は熱力学的平衡状態における各化合物の相分率なので、ソフトマックス関数を出力層に用いることにより、モデル２の出力値に追加演算を行わずに相分率として利用でき、計算コストを低減できる。
学習用データ作成部３は、モデル２の学習用データを作成する。学習用データは、所定の設計条件の入力と、該設計条件に基づき生じ得る熱力学的平衡状態の出力とを含む。学習用データ作成部３は、アルミニウム合金の設計条件（組成、製造条件）の所定範囲の組合せから成る説明変数を作成し、ＣＡＬＰＨＡＤ法を用いて目的変数を算出して、作成した説明変数と、算出した目的変数とを含む学習用データを作成する。
図２は、４元系の学習用データの一例を示す図である。図２の学習用データは、説明変数として、アルミニウム合金の組成のうち添加元素３種のＳｉ、Ｃｕ、Ｍｇの重量百分率（ｗｔ％）、及び製造条件のうち温度（℃）の４個の項目を含み、目的変数として、入力の３種の組成に基づく９種の化合物の相分率を含む、４元系の学習用データである。４元系の学習用データを採用する場合、図１に例示するように、モデル２の入力層のニューロンの数は４個、出力層のニューロンの数は９個となる。
なお、説明変数の組成Ｓｉ、Ｃｕ、Ｍｇの重量百分率は、例えば各元素の数値を所定範囲内で振った数値群を作成し、各元素の数値群の全組合せである。また、説明変数の温度の数値は、所定の温度範囲（例えば０〜１０００℃）内で振った数値群である。これらの組成の全組合せのそれぞれに、温度の数値群が組み合わされて説明変数の組が生成される。目的変数の化合物の数は、説明変数に含む組成の内容に応じて決まる。
図４は、９元系の学習用データの一例を示す図である。図４の学習用データは、説明変数として、アルミニウム合金の組成として添加元素８種Ｓｉ、Ｆｅ、Ｃｕ、Ｍｎ、Ｍｇ、Ｃｒ、Ｎｉ、Ｚｎの重量百分率、及び製造条件の温度（℃）の９個の項目を含み、目的変数として入力の８種の組成に基づく３５種の化合物の相分率を含む、９元系の学習用データである。９元系の学習用データを採用する場合、モデル２の入力層のニューロンの数は９個、出力層のニューロンの数は３５個となる。
なお、説明変数に含める組成の種類としては上記の添加元素３種、８種等に限られず、任意の種類と数を設定できる。また、説明変数に含める製造条件の項目は温度以外のものでもよく、ガス雰囲気などを含めてもよい。
モデル学習部４は、学習用データ作成部３により作成された学習用データを用いて、モデル２の入出力関係が学習用データの入出力関係に近づくように機械学習を行う。本実施形態では、モデル学習部４はディープラーニングを用いてモデル２の機械学習を行う。
図３は、図２に示した４元系の学習データを用いた場合の学習結果を示す図である。図３の（ａ）は学習用データの目的変数から作成した状態図を示し、（ｂ）はその拡大図である。図３の（ｃ）は学習用データの説明変数の入力に基づく学習済みのモデル２の出力から作成した状態図を示し、（ｄ）はその拡大図である。状態図は、目的変数のＳｉ、Ｃｕ、Ｍｇの重量百分率の各組合せごとに、目的変数の温度範囲に亘って作成される。なお図３の例では目的変数は実際には９項目であるが、紙面の都合上代表的なものを図示している。図３に示す各状態図の横軸は温度（℃）を示し、縦軸は熱力学的平衡状態における各化合物の相分率を示す。図３の（ａ）、（ｂ）と（ｃ）、（ｄ）とを比較すれば、４元系では学習済みのモデル２の出力が学習用データに則していることがわかる。
図５は、図４に示した９元系の学習データを用いた場合の学習結果を示す図である。図５の（ａ）は学習用データの目的変数から作成した状態図を示し、（ｂ）は学習用データの説明変数の入力に基づく学習済みのモデル２の出力から作成した状態図を示す。なお図５の例では目的変数は実際には３５項目であるが、紙面の都合上代表的なものを図示している。図５の（ａ）と（ｂ）とを比較すれば、９元系の学習データを用いた場合でも４元系の学習データを用いた場合と同様に学習済みのモデル２の出力が学習用データに則していることがわかる。
なお、モデル学習部４によるモデル２の機械学習は、学習用データの各データセットの出力誤差を解消するよう学習する構成でもよいし、図３（ｃ）や図５（ｂ）に例示したモデル２出力に基づく状態図が、図３（ａ）や図５（ａ）に例示した学習用データの状態図に近づくように学習する構成でもよい。
説明変数設定部５は、対象材料の熱力学的平衡状態の予測に用いる予測用の説明変数を設定する。説明変数設定部５は、例えば表示装置上のＧＵＩなどに各種設計条件の入力画面を表示して設計者に入力を促すことによって、説明変数を設定できる。
予測部６は、説明変数設定部５により設定された予測用説明変数を、モデル学習部４による機械学習が済んだモデル２へ入力することによって、熱力学的平衡状態の予測結果である予測用の目的変数をモデル２から出力する。
状態図表示部７は、予測部６により出力された予測用の目的変数に基づき、対象材料の熱力学的平衡状態の状態図を作成して表示する。
図６は、予測装置１のハードウェア構成を示すブロック図である。図６に示すように、予測装置１は、物理的には、ＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）１０１、ＧＰＵ（Ｇｒａｐｈｉｃｓ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）１０８、主記憶装置であるＲＡＭ（Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ）１０２およびＲＯＭ（Ｒｅａｄ Ｏｎｌｙ Ｍｅｍｏｒｙ）１０３、入力デバイスであるキーボード及びマウス等の入力装置１０４、ディスプレイ等の出力装置１０５、ネットワークカード等のデータ送受信デバイスである通信モジュール１０６、ハードディスク等の補助記憶装置１０７、などを含むコンピュータシステムとして構成することができる。
図１に示す予測装置１の各機能は、ＣＰＵ１０１、ＲＡＭ１０２等のハードウェア上に所定のコンピュータソフトウェア（予測プログラム）を読み込ませることにより、ＣＰＵ１０１の制御のもとで通信モジュール１０６、入力装置１０４、出力装置１０５を動作させるとともに、ＲＡＭ１０２や補助記憶装置１０７におけるデータの読み出し及び書き込みを行うことで実現される。すなわち、本実施形態の予測プログラムをコンピュータ上で実行させることで、予測装置１は、図１のモデル２と、学習用データ作成部３と、モデル学習部４と、説明変数設定部５と、予測部６と、状態図表示部７として機能する。
本実施形態の予測プログラムは、例えばコンピュータが備える記憶装置内に格納される。なお、予測プログラムは、その一部又は全部が、通信回線等の伝送媒体を介して伝送され、コンピュータが備える通信モジュール１０６等により受信されて記録（インストールを含む）される構成としてもよい。また、予測プログラムは、その一部又は全部が、ＣＤ−ＲＯＭ、ＤＶＤ−ＲＯＭ、フラッシュメモリなどの持ち運び可能な記憶媒体に格納された状態から、コンピュータ内に記録（インストールを含む）される構成としてもよい。
図７は、実施形態に係る予測装置１により実施される熱力学的平衡状態の予測処理のフローチャートである。
ステップＳ１では、学習用データ作成部３により、モデル２の学習用データが作成される（学習用データ作成ステップ）。学習用データ作成部３は、設計者による組成や製造条件（温度など）の指定範囲の入力に応じて、この指定範囲内を網羅するように説明変数の組を作成し、作成したすべての説明変数の組を使ってＣＡＬＰＨＡＤ法により目的変数を算出し、説明変数と紐づけて、例えば図２、図４に示す４元系や９元系の学習用データを作成する。
ステップＳ２では、モデル学習部４により、ステップＳ１にて作成された学習用データを用いて、モデル２の機械学習が行われる（モデル学習ステップ）。モデル学習部４は、学習用データの説明変数の入力に応じて、この説明変数に紐づけられた目的変数と合致した出力となるように、多層ニューラルネットワークの各層間のウェイトを調整して、モデル２の学習を行う。モデル学習部４は、例えばディープラーニングを用いてモデル２の学習を行う。
ステップＳ３では、説明変数設定部５により、対象材料の熱力学的平衡状態の予測に用いる予測用の説明変数が設定される（説明変数設定ステップ）。
ステップＳ４では、予測部６により、学習済みのモデル２を用いて対象材料の熱力学的平衡状態が予測される（予測ステップ）。予測部６は、ステップＳ３にて設定された説明変数をステップＳ２にて機械学習が済んだモデル２へ入力し、モデル２から出力される熱力学的平衡状態の予測結果である予測用の目的変数を取得する。
ステップＳ５では、状態図表示部７により、予測部６により出力された予測用の目的変数に基づき、対象材料の熱力学的平衡状態の状態図が作成されて表示される。
本実施形態の効果を説明する。本実施形態の予測装置１は、対象材料の設計条件に係る説明変数の入力に基づいて、熱力学的平衡状態に係る目的変数を出力するモデル２と、所定の設計条件の入力と、該設計条件に基づき生じ得る熱力学的平衡状態の出力とを含む学習用データを作成する学習用データ作成部３と、学習用データ作成部３により作成された学習用データを用いて、モデル２の入出力関係が学習用データの入出力関係に近づくように機械学習を行うモデル学習部４と、対象材料の熱力学的平衡状態の予測に用いる予測用説明変数を設定する説明変数設定部５と、モデル学習部４による機械学習が済んだモデル２への予測用説明変数の入力に基づき、熱力学的平衡状態の予測結果である予測用目的変数をモデル２から出力する予測部６と、を備える。
この構成により、学習済みのモデル２に予測用の目的変数を入力するだけで、この目的変数に応じた熱力学的平衡状態の予測結果である予測用目的変数がモデル２から出力されるので、従来のＣＡＬＰＨＡＤ法などのシミュレーションに比べて目的変数を算出するまでの計算コストを格段に低減でき、熱力学的平衡状態を短時間で計算できる。例えば１００００通りの説明変数（組成や製造条件の組）をスクリーニングすべく全組成の状態図を算出しようとすると、シミュレーションの場合、各組の計算に約９０秒かかるので、合計で約２５０時間の時間を要するが、本実施形態の場合、各組の計算時間は約３ミリ秒なので、合計で約３０秒で済む。
また、この構成により、学習済みのモデル２の汎化能力を利用して、学習用データの説明変数とは異なる未知の入力に対しても適切な予測用の目的変数を出力できるので、熱力学的平衡状態の予測を高精度に行うことができる。つまり、従来のシミュレーション手法では未知の入力の場合には再度シミュレーション演算を行う必要があったのに対して、本実施形態では一度学習済みのモデル２を獲得すれば、モデル２の追加学習を行わなくても未知の入力に対して適切な出力を得ることができる。これにより、組成の網羅的な解析を短時間かつ高精度で実施できるので、従来のシミュレーションでは現実的ではなかった大量の組成のスクリーニングも可能となり、より最適な設計条件の抽出が可能となる。
また、本実施形態の予測装置１は、予測部６により出力された予測用目的変数に基づき熱力学的平衡状態の状態図を作成して表示する状態図表示部７を備える。この構成により、設計者に対して対象材料の熱力学的平衡状態の予測結果を視覚的に提示できるので、設計者が予測結果を容易に理解できる。
また、本実施形態の予測装置１では、モデル２が多層ニューラルネットワークであり、モデル学習部４はディープラーニングを用いてモデルの学習を行うので、学習を高速かつ高精度に行うことができ、熱力学的平衡状態の予測をより高精度にできる。
また、本実施形態の予測装置１では、モデル２の出力である目的変数は、対象材料の熱力学的平衡状態における相分率であり、モデル２の多層ニューラルネットワークの出力層のニューロンはソフトマックス関数に基づき出力値を算出する。この構成により、出力の総和が常に１となるように維持できるので、相分率をモデル２の出力とする本実施形態では、各出力をそのまま相分率として使用でき、計算コストをさらに低減できる。
また、本実施形態の予測装置１では、学習用データ作成部３は、説明変数を作成し、ＣＡＬＰＨＡＤ法を用いて目的変数を算出して学習用データを作成するので、目的変数を精度良く算出でき、学習用データの精度を向上できる。
以上、具体例を参照しつつ本実施形態について説明した。しかし、本開示はこれらの具体例に限定されるものではない。これら具体例に、当業者が適宜設計変更を加えたものも、本開示の特徴を備えている限り、本開示の範囲に包含される。前述した各具体例が備える各要素およびその配置、条件、形状などは、例示したものに限定されるわけではなく適宜変更することができる。前述した各具体例が備える各要素は、技術的な矛盾が生じない限り、適宜組み合わせを変えることができる。
上記実施形態では、予測装置１が熱力学系平衡状態を予測する対象材料としてアルミニウム合金を挙げたが、アルミニウム合金以外の合金系を対象材料としてもよい。このような合金系としては、Ｆｅ合金、Ｃｕ合金、Ｎｉ合金、Ｃｏ合金、Ｔｉ合金、Ｍｇ合金、Ｍｎ合金、Ｚｎ合金などが挙げられる。また、セラミックス、水溶液、化学反応等の合金以外を対象材料としてもよい。
上記実施形態では、モデル２として多層ニューラルネットワークを例示し、モデル２の機械学習手法としてディープラーニングを例示したが、モデル２や学習手法はこれに限られず、遺伝的アルゴリズムなどの他の教師あり学習モデルや、ランダムフォレスト回帰やカーネルリッジ回帰などの他の機械学習手法を用いることもできる。
上記実施形態では、学習用データ作成部３が、モデル２の学習用データをＣＡＬＰＨＡＤ法で計算して作成する構成を例示したが、熱力学計算等を行ったシミュレーション結果ではなく実験結果を用いて学習用データを作成する構成でもよい。例えば、組成を変えて製造した種々の合金材料について、種々の温度に保持して平衡状態に達した後、この平衡状態を凍結するため低温へ急冷して、各種分析によって合金材料中における化合物の種類と相分率を算出する方法が挙げられる。分析の方法としては、例えばＸ線回折測定により得られたピーク位置から化合物の種類を特定し、ピーク強度比から相分率を算出する方法や、走査型電子顕微（Ｓｃａｎｎｉｎｇ Ｅｌｅｃｔｒｏｎ Ｍｉｃｒｏｓｃｏｐｅ（ＳＥＭ））や透過型電子顕微鏡（Ｔｒａｎｓｍｉｓｓｉｏｎ Ｅｌｅｃｔｒｏｎ Ｍｉｃｒｏｓｃｏｐｅ（ＴＥＭ））といった電子顕微鏡観察において、観察した対象視野に対してエネルギー分散型Ｘ線分析（Ｅｎｅｒｇｙ ｄｉｓｐｅｒｓｉｖｅ Ｘ−ｒａｙ ｓｐｅｃｔｒｏｍｅｔｒｙ（ＥＤＳ））を行い、相の種類を特定するとともに、観察像における各相の粒子のそれぞれの合計面積率から相分率を算出する方法が考えられる。
上記実施形態では、予測装置１が状態図表示部７によって、最終的に熱力学的平衡状態の状態図を出力する構成を例示したが、予測装置１が予測部６から出力される予測結果の数値を出力し、状態図を出力しない構成でもよい。
上記実施形態では、モデル２の出力が熱力学的平衡状態における化合物の相分率である構成を例示したが、モデル２の出力は熱力学的平衡状態に係る情報であればよく、相分率以外でもよい。
本国際出願は２０１８年１０月３１日に出願された日本国特許出願２０１８−２０６０１７号に基づく優先権を主張するものであり、２０１８−２０６０１７号の全内容をここに本国際出願に援用する。
１ 熱力学的平衡状態の予測装置
２ モデル
３ 学習用データ作成部
４ モデル学習部
５ 説明変数設定部
６ 予測部
７ 状態図表示部
２ モデル
３ 学習用データ作成部
４ モデル学習部
５ 説明変数設定部
６ 予測部
７ 状態図表示部
Claims (8)
- 対象材料の熱力学的平衡状態の予測装置であって、
前記対象材料の設計条件に係る説明変数の入力に基づいて、前記熱力学的平衡状態に係る目的変数を出力するモデルと、
所定の設計条件の入力と、該設計条件に基づき生じ得る前記熱力学的平衡状態の出力とを含む学習用データを作成する学習用データ作成部と、
前記学習用データ作成部により作成された前記学習用データを用いて、前記モデルの入出力関係が前記学習用データの入出力関係に近づくように機械学習を行うモデル学習部と、
前記対象材料の熱力学的平衡状態の予測に用いる予測用説明変数を設定する説明変数設定部と、
前記モデル学習部による機械学習が済んだ前記モデルへの前記予測用説明変数の入力に基づき、前記熱力学的平衡状態の予測結果である予測用目的変数を前記モデルから出力する予測部と、
を備える予測装置。 - 前記予測部により出力された前記予測用目的変数に基づき前記熱力学的平衡状態の状態図を作成して表示する状態図表示部を備える、請求項１に記載の予測装置。
- 前記モデルは多層ニューラルネットワークであり、
前記モデル学習部はディープラーニングを用いて前記モデルの学習を行う、
請求項１または２に記載の予測装置。 - 前記モデルの出力である前記目的変数は、前記対象材料の熱力学的平衡状態における相分率であり、
前記多層ニューラルネットワークの出力層のニューロンにはソフトマックス関数が用いられる、
請求項３に記載の予測装置。 - 前記学習用データ作成部は、前記設計条件の所定範囲の組合せから成る前記説明変数を作成し、ＣＡＬＰＨＡＤ法を用いて前記目的変数を算出して、前記作成した前記説明変数と前記算出した前記目的変数とを含む前記学習用データを作成する、
請求項１〜４のいずれか１項に記載の予測装置。 - 前記対象材料がアルミニウム合金であり、
前記説明変数が前記アルミニウム合金の組成及び製造条件を含み、
前記目的変数が前記アルミニウム合金の熱力学的平衡状態での相分率を含む、
請求項１〜５のいずれか１項に記載の予測装置。 - コンピュータが実行する、対象材料の熱力学的平衡状態の予測方法であって、
前記対象材料の設計条件に係る説明変数の入力に基づいて、前記熱力学的平衡状態に係る目的変数を出力するモデルに関して、所定の設計条件の入力と、該設計条件に基づき生じ得る前記熱力学的平衡状態の出力とを含む学習用データを作成する学習用データ作成ステップと、
前記学習用データ作成ステップにて作成された前記学習用データを用いて、前記モデルの入出力関係が前記学習用データの入出力関係に近づくように機械学習を行うモデル学習ステップと、
前記対象材料の熱力学的平衡状態の予測に用いる予測用説明変数を設定する説明変数設定ステップと、
前記説明変数設定ステップにて設定された前記予測用説明変数の、前記モデル学習ステップにて機械学習が済んだ前記モデルへの入力に基づき、前記熱力学的平衡状態の予測結果である予測用目的変数を前記モデルから出力する予測ステップと、
を含む予測方法。 - 対象材料の熱力学的平衡状態の予測プログラムであって、
前記対象材料の設計条件に係る説明変数の入力に基づいて、前記熱力学的平衡状態に係る目的変数を出力するモデルに関して、所定の設計条件の入力と、該設計条件に基づき生じ得る前記熱力学的平衡状態の出力とを含む学習用データを作成する学習用データ作成機能と、
前記学習用データ作成機能により作成された前記学習用データを用いて、前記モデルの入出力関係が前記学習用データの入出力関係に近づくように機械学習を行うモデル学習機能と、
前記対象材料の熱力学的平衡状態の予測に用いる予測用説明変数を設定する説明変数設定機能と、
前記説明変数設定機能により設定された前記予測用説明変数の、前記モデル学習機能により機械学習が済んだ前記モデルへの入力に基づき、前記熱力学的平衡状態の予測結果である予測用目的変数を前記モデルから出力する予測機能と、
をコンピュータに実現させるための予測プログラム。
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2018206017 | 2018-10-31 | ||
| JP2018206017 | 2018-10-31 | ||
| PCT/JP2019/041716 WO2020090617A1 (ja) | 2018-10-31 | 2019-10-24 | 熱力学的平衡状態の予測装置、予測方法、及び予測プログラム |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPWO2020090617A1 JPWO2020090617A1 (ja) | 2021-09-02 |
| JP6963119B2 true JP6963119B2 (ja) | 2021-11-05 |
Family
ID=70462444
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2020553834A Active JP6963119B2 (ja) | 2018-10-31 | 2019-10-24 | 熱力学的平衡状態の予測装置、予測方法、及び予測プログラム |
Country Status (5)
| Country | Link |
|---|---|
| US (1) | US20210406433A1 (ja) |
| EP (1) | EP3876241A4 (ja) |
| JP (1) | JP6963119B2 (ja) |
| CN (1) | CN112912969A (ja) |
| WO (1) | WO2020090617A1 (ja) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111681714B (zh) * | 2020-07-02 | 2023-06-20 | 兰州大学 | 一种在定向凝固包晶合金中生长非典型三次枝晶的方法 |
| JP2022132895A (ja) * | 2021-03-01 | 2022-09-13 | 株式会社Uacj | 合金材料の特性を予測する製造支援システム、予測モデルを生成する方法およびコンピュータプログラム |
| JP2022190226A (ja) * | 2021-06-14 | 2022-12-26 | トヨタ自動車株式会社 | 圧力容器のひずみ解析装置および圧力容器の製造方法 |
| WO2023095680A1 (ja) | 2021-11-24 | 2023-06-01 | 株式会社レゾナック | 予測装置、学習装置、予測方法、学習方法、予測プログラム及び学習プログラム |
| WO2024085071A1 (ja) * | 2022-10-17 | 2024-04-25 | 株式会社レゾナック | 組成最適化装置、組成最適化方法及び組成最適化プログラム |
Family Cites Families (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH1055348A (ja) * | 1996-08-08 | 1998-02-24 | Bridgestone Corp | 多成分系材料最適化解析装置および方法 |
| CN102663498B (zh) * | 2012-04-28 | 2014-06-18 | 武汉大学 | 一种9%Cr马氏体耐热钢焊缝金属Ac1点的预测方法 |
| JP5749233B2 (ja) | 2012-08-31 | 2015-07-15 | 株式会社東芝 | 材料組織の予測方法及び予測システム |
| GB201302743D0 (en) * | 2013-02-18 | 2013-04-03 | Rolls Royce Plc | Method and system for designing a material |
| WO2015054637A1 (en) * | 2013-10-10 | 2015-04-16 | Scoperta, Inc. | Methods of selecting material compositions and designing materials having a target property |
| US20160034614A1 (en) * | 2014-08-01 | 2016-02-04 | GM Global Technology Operations LLC | Materials property predictor for cast aluminum alloys |
| JP6630640B2 (ja) * | 2016-07-12 | 2020-01-15 | 株式会社日立製作所 | 材料創成装置、および材料創成方法 |
| US10830747B2 (en) * | 2016-10-26 | 2020-11-10 | Northwestern University | System and method for predicting fatigue strength of alloys |
| JP6730216B2 (ja) * | 2017-03-23 | 2020-07-29 | 株式会社日立製作所 | エレベータ管理システム、及び、エレベータの管理方法 |
| JP6711784B2 (ja) | 2017-06-02 | 2020-06-17 | 小澤 隆 | 建築図面作成システム及び建築図面作成方法 |
-
2019
- 2019-10-24 CN CN201980070275.3A patent/CN112912969A/zh active Pending
- 2019-10-24 US US17/288,185 patent/US20210406433A1/en active Pending
- 2019-10-24 EP EP19880330.6A patent/EP3876241A4/en active Pending
- 2019-10-24 WO PCT/JP2019/041716 patent/WO2020090617A1/ja unknown
- 2019-10-24 JP JP2020553834A patent/JP6963119B2/ja active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN112912969A (zh) | 2021-06-04 |
| WO2020090617A1 (ja) | 2020-05-07 |
| EP3876241A4 (en) | 2022-01-19 |
| US20210406433A1 (en) | 2021-12-30 |
| EP3876241A1 (en) | 2021-09-08 |
| JPWO2020090617A1 (ja) | 2021-09-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6963119B2 (ja) | 熱力学的平衡状態の予測装置、予測方法、及び予測プログラム | |
| Dimiduk et al. | Perspectives on the impact of machine learning, deep learning, and artificial intelligence on materials, processes, and structures engineering | |
| Ramprasad et al. | Machine learning in materials informatics: recent applications and prospects | |
| Johnson et al. | Invited review: Machine learning for materials developments in metals additive manufacturing | |
| Zeng et al. | Revealing high-fidelity phase selection rules for high entropy alloys: A combined CALPHAD and machine learning study | |
| Balachandran et al. | Adaptive strategies for materials design using uncertainties | |
| JP6976456B2 (ja) | 材料設計装置、材料設計方法、及び材料設計プログラム | |
| WO2020203922A1 (ja) | 結晶形予測装置、結晶形予測方法、ニューラルネットワークモデルの製造方法、及びプログラム | |
| Nguyen et al. | A novel binary particle swarm optimization algorithm and its applications on knapsack and feature selection problems | |
| Kaufmann et al. | Efficient few-shot machine learning for classification of EBSD patterns | |
| Del Vecchio et al. | Support Vector Representation Machine for superalloy investment casting optimization | |
| CN114207729A (zh) | 材料特性预测系统以及材料特性预测方法 | |
| WO2006126647A1 (ja) | 石油精製プラントのパラメータ予測装置及びパラメータ予測方法 | |
| JP2022132895A (ja) | 合金材料の特性を予測する製造支援システム、予測モデルを生成する方法およびコンピュータプログラム | |
| Steingrimsson et al. | Machine learning and data analytics for design and manufacturing of high-entropy materials exhibiting mechanical or fatigue properties of interest | |
| Jarosz et al. | Metamodeling and optimization of a blister copper two-stage production process | |
| LeSar | Materials informatics: an emerging technology for materials development | |
| Gao et al. | Data-driven machine learning for alloy research: Recent applications and prospects | |
| JPWO2016203757A1 (ja) | 制御装置、それを使用する情報処理装置、制御方法、並びにコンピュータ・プログラム | |
| Shen et al. | Predicting phases and hardness of high entropy alloys based on machine learning | |
| JP7456273B2 (ja) | データ解析システム、データ解析方法及びデータ解析プログラム | |
| Virdi et al. | Feature selection with LASSO and VSURF to model mechanical properties for investment casting | |
| Vazquez et al. | A deep neural network regressor for phase constitution estimation in the high entropy alloy system Al-Co-Cr-Fe-Mn-Nb-Ni | |
| Abu Doush et al. | A sensitivity analysis for harmony search with multi-parent crossover algorithm | |
| Singh et al. | Machine learning for weather forecasting: XGBoost vs SVM vs random forest in predicting temperature for visakhapatnam |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20210310 |
|
| A871 | Explanation of circumstances concerning accelerated examination |
Free format text: JAPANESE INTERMEDIATE CODE: A871 Effective date: 20210310 |
|
| A975 | Report on accelerated examination |
Free format text: JAPANESE INTERMEDIATE CODE: A971005 Effective date: 20210409 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20210608 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20210720 |
|
| 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: 20211012 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20211014 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 6963119 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |