JP5919384B2 - 材料組織予測装置、製品製造方法及び材料組織予測方法 - Google Patents

材料組織予測装置、製品製造方法及び材料組織予測方法 Download PDF

Info

Publication number
JP5919384B2
JP5919384B2 JP2014532693A JP2014532693A JP5919384B2 JP 5919384 B2 JP5919384 B2 JP 5919384B2 JP 2014532693 A JP2014532693 A JP 2014532693A JP 2014532693 A JP2014532693 A JP 2014532693A JP 5919384 B2 JP5919384 B2 JP 5919384B2
Authority
JP
Japan
Prior art keywords
calculation
phase
nucleation
calculated
temperature
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
Application number
JP2014532693A
Other languages
English (en)
Other versions
JPWO2014033928A1 (ja
Inventor
佐野　光彦
光彦 佐野
小原　一浩
一浩 小原
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Toshiba Mitsubishi Electric Industrial Systems Corp
Original Assignee
Toshiba Mitsubishi Electric Industrial Systems Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Toshiba Mitsubishi Electric Industrial Systems Corp filed Critical Toshiba Mitsubishi Electric Industrial Systems Corp
Application granted granted Critical
Publication of JP5919384B2 publication Critical patent/JP5919384B2/ja
Publication of JPWO2014033928A1 publication Critical patent/JPWO2014033928A1/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B21—MECHANICAL METAL-WORKING WITHOUT ESSENTIALLY REMOVING MATERIAL; PUNCHING METAL
    • B21B—ROLLING OF METAL
    • B21B37/00—Control devices or methods specially adapted for metal-rolling mills or the work produced thereby
    • C—CHEMISTRY; METALLURGY
    • C21—METALLURGY OF IRON
    • C21D—MODIFYING THE PHYSICAL STRUCTURE OF FERROUS METALS; GENERAL DEVICES FOR HEAT TREATMENT OF FERROUS OR NON-FERROUS METALS OR ALLOYS; MAKING METAL MALLEABLE, e.g. BY DECARBURISATION OR TEMPERING
    • C21D11/00—Process control or regulation for heat treatments
    • C—CHEMISTRY; METALLURGY
    • C21—METALLURGY OF IRON
    • C21D—MODIFYING THE PHYSICAL STRUCTURE OF FERROUS METALS; GENERAL DEVICES FOR HEAT TREATMENT OF FERROUS OR NON-FERROUS METALS OR ALLOYS; MAKING METAL MALLEABLE, e.g. BY DECARBURISATION OR TEMPERING
    • C21D11/00—Process control or regulation for heat treatments
    • C21D11/005—Process control or regulation for heat treatments for cooling
    • C—CHEMISTRY; METALLURGY
    • C21—METALLURGY OF IRON
    • C21D—MODIFYING THE PHYSICAL STRUCTURE OF FERROUS METALS; GENERAL DEVICES FOR HEAT TREATMENT OF FERROUS OR NON-FERROUS METALS OR ALLOYS; MAKING METAL MALLEABLE, e.g. BY DECARBURISATION OR TEMPERING
    • C21D9/00—Heat treatment, e.g. annealing, hardening, quenching or tempering, adapted for particular articles; Furnaces therefor
    • C21D9/52—Heat treatment, e.g. annealing, hardening, quenching or tempering, adapted for particular articles; Furnaces therefor for wires; for strips ; for rods of unlimited length
    • C21D9/54—Furnaces for treating strips or wire
    • C21D9/56—Continuous furnaces for strip or wire
    • C21D9/573—Continuous furnaces for strip or wire with cooling
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/20—Metals
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/20—Metals
    • G01N33/204—Structure thereof, e.g. crystal structure
    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B15/00—Systems controlled by a computer
    • G05B15/02—Systems controlled by a computer electric
    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D23/00—Control of temperature
    • G05D23/19—Control of temperature characterised by the use of electric means
    • G05D23/1917—Control of temperature characterised by the use of electric means using digital means
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • 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
    • C—CHEMISTRY; METALLURGY
    • C21—METALLURGY OF IRON
    • C21D—MODIFYING THE PHYSICAL STRUCTURE OF FERROUS METALS; GENERAL DEVICES FOR HEAT TREATMENT OF FERROUS OR NON-FERROUS METALS OR ALLOYS; MAKING METAL MALLEABLE, e.g. BY DECARBURISATION OR TEMPERING
    • C21D2211/00—Microstructure comprising significant phases
    • C21D2211/004—Dispersions; Precipitations
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N25/00—Investigating or analyzing materials by the use of thermal means
    • G01N25/02—Investigating or analyzing materials by the use of thermal means by investigating changes of state or changes of phase; by investigating sintering
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/24—Sheet material

Landscapes

  • Engineering & Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mechanical Engineering (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Health & Medical Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Materials Engineering (AREA)
  • Metallurgy (AREA)
  • Organic Chemistry (AREA)
  • Thermal Sciences (AREA)
  • General Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • General Health & Medical Sciences (AREA)
  • Biochemistry (AREA)
  • Analytical Chemistry (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Automation & Control Theory (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Medicinal Chemistry (AREA)
  • Food Science & Technology (AREA)
  • Computer Hardware Design (AREA)
  • Geometry (AREA)
  • Evolutionary Computation (AREA)
  • Software Systems (AREA)
  • Operations Research (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Probability & Statistics with Applications (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Investigating And Analyzing Materials By Characteristic Methods (AREA)

Description

本発明は、温度変化により相変態をする材料の組織を予測することができる材料組織予測装置、製品製造方法及び材料組織予測方法に関する。
近年、環境問題及び省エネルギーの観点から、工業製品に用いる材料の材質の向上が望まれている。また、材料においては、引張強度、降伏応力、又は伸びなどで代表される材質特性は、顕微鏡レベルの材料組織と密接に関係していることが知られている。このため、熱間加工後条件、又は、その後の冷却条件を操作して、材料に与える温度変化の指標となる温度変化曲線を変化させることにより、製品の材質特性の創り込みがされている。
また、材料を冷却後、研磨して、操作型電子顕微鏡を応用したＥＢＳＤ(electron back scattering diffraction)装置で観察すると、析出相（新相）の各結晶粒の結晶方位が得られる。母相の結晶方位と析出相の結晶方位には一定の関係がある。これを利用して母相の結晶方位と粒界を再構築できることが知られている（非特許文献１参照）。また、冷却中の拡散変態の進行の定式化に用いる核生成速度及び成長速度に関する研究がされている（非特許文献２参照）。
温度制御装置としては、圧延材が相変態を起こすことにより発生する変態発熱の量を変態発熱モデルを用いて予測して、その変態発熱の量を補償しながら、巻取温度を所定の温度目標値に一致させるように制御することが開示されている（特許文献１参照）。
所望の材質特性を得るための適正な製造条件を確立するためには、試作を繰り返す必要がある。しかし、このような試作を繰り返す作業は、膨大な費用と手間が掛かる。そこで、鉄鋼材料では、熱間圧延プロセスを対象として、製品である材料の材料組織及び材質を予測する組織材質予測技術が開発されている。組織材質予測技術は、所望の材質特性を得るために温度変化曲線などの製造条件を適正化するために利用されている。
ところが、既存の組織材質予測モデルは、材料組織内を均一であるものとして取り扱っている。従って、材料組織内での結晶粒の大きさ及び形状についての違い又はばらつきについては考慮されていない。しかし、実際には、製造条件によっては、材料組織内で、結晶粒毎に大きさ及び形状が異なり、かつこれらのばらつきが大きくなる場合がある。これらは、材質が変動する一因となる。
一方、材料の熱間状態からの冷却過程における固相変態は、(i)結晶の核生成と、(ii) 生成した結晶の粒成長の２つの現象からなる。後者（結晶粒の成長）に対し、フェーズフィールド(MPF, multi-phase-field)法とよばれる解析手法が知られている。例えば、単結晶成長において、フェーズフィールド法により粒成長によるグレイン形状の変化を計算するシミュレーション方法が開示されている（特許文献２参照）。しかし、前者（核生成）については、実際の核生成の位置とタイミングは母相組織と温度変化により決まるにも関わらず、従来のフェーズフィールド法では、計算開始時の初期核配置を手動で与えるなどごく簡易的な考慮に留まっている。特許文献２には、単結晶において、転位によるひずみエネルギーにより引き起こされる核生成が考慮されているものの、鉄鋼を初めとする多結晶金属においては、母相組織の粒界付近で核生成が生じやすいことなどが知られているように、特許文献２の方法の多結晶金属の核生成への適用には限界がある。
したがって、従来のフェーズフィールド法では、材料組織内での結晶粒の大きさ及び形状の違い及びこれらのばらつきについて、粒成長による定性的な変化の傾向を把握することができたものの、核生成挙動の考慮が十分でないために、定量的な予測精度は不十分である。従って、冷却後の材質の予測など、定量的な評価が必要となる用途への、この解析手法の適用は困難である。
このように、材質などの定量的な評価に用いるために、従来のフェーズフィールド法によって、温度変化により相変態をする材料の組織を予測することは困難である。
特開２００５−２９７０１５号公報 特開２００７−０９１５４３号公報
森本 敬治、外４名，「ＥＢＳＤデータによるバリアント解析プログラムの開発」，「鉄と鋼」，社団法人日本鉄鋼協会，２００７年，Ｖｏｌ．９３，Ｎｏ．９，ｐ．２７−３５ 末広 正芳、外４名，「低炭素普通鋼の冷却中の変態進行の定式化」，「鉄と鋼」，社団法人日本鉄鋼協会，１９８７年，第７３年，第８号，ｐ．１１０−１１７
本発明が解決しようとする課題は、温度変化により相変態をする材料の組織を定量的に予測することができる材料組織予測装置、製品製造方法及び材料組織予測方法を提供することにある。
本発明の観点に従った材料組織予測装置は、温度変化により母相から析出相に相変態する材料の組織を予測する材料組織予測装置であって、前記材料に前記温度変化を与えるための温度条件に基づいて、前記材料の計算対象領域内の複数の計算点の温度を計算する温度計算手段と、前記温度計算手段により計算された前記複数の計算点の温度に基づいて、前記計算対象領域内の核生成回数を計算する核生成回数計算手段と、前記核生成回数計算手段により計算された前記核生成回数及び予め決められた前記複数の計算点の核生成の起こり易さに基づいて、前記複数の計算点から前記析出相の核を生成する析出相生成点を決定する析出相生成点決定手段と、前記析出相生成点決定手段により前記析出相生成点が繰り返し決定される度に、前記析出相生成点における前記析出相の粒成長を計算する粒成長計算手段と、前記粒成長計算手段により計算された前記析出相の粒成長に基づいて、前記材料の組織を予測する材料組織予測手段とを備える。
図１は、本発明の実施形態に係る材料組織予測装置の構成を示す構成図である。 図２は、実施形態に係る材料組織予測装置が適用された薄板熱間圧延ラインの構成を示す構成図である。 図３は、実施形態に係る薄板熱間圧延ラインにおける材料の計算対象領域の位置を示す模式図である。 図４は、実施形態に係る核生成順位決定部により各計算点をレベル分けした金属材サンプルの模式図である。 図５は、実施形態に係る核生成順位決定部により各計算点の核生成順位を決定した金属材サンプルの模式図である。
以下、図面を参照して、本発明の実施形態を説明する。
（実施形態）
図１は、本発明の実施形態に係る材料組織予測装置１の構成を示す構成図である。なお、図中における同一部分には同一符号を付してその詳しい説明を省略し、異なる部分について主に述べる。
まず、材料組織予測装置１が適用される一般的な金属材加工ラインについて説明する。
金属材加工ラインでは、金属材の加熱又は冷却の少なくとも一方を行う。
加熱装置は、誘導加熱装置又はガス加熱装置などである。冷却装置は、水冷装置、油冷装置又は強制空冷装置などである。これらの装置は、複数あってもよい。また、加熱装置と冷却装置の配置順序は任意である。
ここでは、金属材加工ラインの一例として薄板熱間圧延ライン２に材料組織予測装置１を適用した構成について説明する。また、金属材ＭＴを主に冷却する場合について説明するが、金属材ＭＴを加熱する場合も同様である。
図２は、実施形態に係る材料組織予測装置１が適用された薄板熱間圧延ライン２の構成を示す構成図である。図２では、薄板熱間圧延ライン２のランアウトテーブル付近の構成を示している。
薄板熱間圧延ライン２は、金属材ＭＴから所望の材質の製品を製造するための設備である。薄板熱間圧延ライン２には、仕上圧延機１１、水冷装置１２、巻取機１３、搬送ローラ１４、複数の放射温度計１５、及び光学式金属検知器１６が設けられている。
金属材ＭＴは、加熱された金属材ＭＴが仕上圧延機１１により圧延された後、搬送ローラ１４により搬送されつつ、水冷装置１２より冷却される。水冷装置１２は、例えば、ラミナー水冷装置である。水冷装置１２より冷却された金属材ＭＴは、巻取機１３に巻き取られる。放射温度計１５及び光学式金属検知器１６により検出された金属材ＭＴの各種データは、温度制御装置５に送信される。
温度制御装置５は、金属材ＭＴの温度変化を予め設定された冷却曲線（温度変化曲線）に一致させるように、放射温度計１５により計測された温度及び材料組織予測装置１から受信した金属材ＭＴの変態に伴う発熱の予測値に基づいて、水冷装置１２を制御する。
材料組織予測装置１は、主にコンピュータで構成されている。材料組織予測装置１は、薄板熱間圧延ライン２から受信した金属材ＭＴの各種データに基づいて、金属材ＭＴの組織の予測及び金属材ＭＴの材質評価などをするための演算処理をする。材料組織予測装置１は、演算処理結果を温度制御装置５に送信する。演算処理結果には、変態に伴う発熱の予測値などの金属材ＭＴについての情報が含まれている。
材料組織予測装置１は、金属材ＭＴの組織の予測及び材質の評価の計算を、金属材ＭＴが薄板熱間圧延ライン２の所定の位置に到達した場合、又は、材料組織予測装置１の操作画面から計算開始の指示があった場合などに開始する。所定の位置とは、例えば、水冷装置１２の上流側に設置された放射温度計１５の直下である。所定の位置への金属材ＭＴの到着は、ライン上の要所に設置された光学式金属検知器１６による検出信号又は搬送ローラ１４の回転計による回転数に基づいて計算される搬送距離実績値などにより検知される。
金属材ＭＴの組織の予測をするための計算対象となる計算対象領域は、予め決められている。計算対象領域は、金属材ＭＴの温度条件又は冷却条件（加熱条件）などの材料組織に影響がある処理条件を代表する位置となる。計算対象領域の数は任意であるが、多過ぎると計算量が増大する。ここでは、薄板熱間圧延ライン２における金属材ＭＴの計算対象領域は、図３に示すように、圧延後の金属材ＭＴの長手方向の先端から５ｍの位置Ｐ１、長手方向の中央の位置Ｐ２、及び長手方向の尾端から５ｍの位置Ｐ３の３点とする。また、これらの位置Ｐ１〜Ｐ３は、板表面から内部方向に板厚の１／４の部分とする。
また、各計算対象領域の大きさが小さ過ぎると、結晶粒の形状特性値を正しく決定できない。一方、計算対象領域の大きさがある一定以上あれば、形状特性値の統計量の精度はあまり変わらない。但し、計算対象領域の大きさが大き過ぎると、計算量が増大する。従って、結晶粒の形状特性値を評価する場合には、想定される析出相の結晶粒径の５倍〜５０００倍とするのが望ましい。
図３を参照して、計算点Ｐｃの配置について説明する。材料組織予測装置１は、各計算対象領域内に多数の計算点Ｐｃをメッシュ状に配置する。これらの計算点Ｐｃの間隔は、大き過ぎると結晶粒の形状特性値を正しく決定できず、小さ過ぎると計算量が増大する。従って、計算対象領域内の計算点Ｐｃの間隔は、想定される析出相の結晶粒径の１／５００〜１／５とするのが望ましい。
材料組織予測装置１は、各計算点の情報を保持するため、温度計算値、核生成順位、変態に伴い拡散する元素（例えば、オーステナイトからフェライトに変化する場合は炭素）の濃度ｃ、及びスカラー量φについての記憶領域を確保する。
温度計算値、核生成順位、及び拡散元素濃度の記憶領域の容量は、(計算対象領域の数)×(計算対象領域の大きさ)／（計算点の間隔）となる。
スカラー量φの記憶領域の容量は、（結晶の総数）×(計算対象領域の数)×(計算対象領域の大きさ)／（計算点の間隔）となる。
スカラー量φは、各計算点における結晶粒の存在確率を表す。計算点が結晶粒の内部の場合、スカラー量φは１である。計算点が結晶粒の外部の場合、スカラー量φは０である。計算点が粒界の場合、スカラー量φは０よりも大きく１よりも小さい。
図１は、本実施形態に係る材料組織予測装置１の構成を示すブロック図である。図１中の矢印は、主に情報（又はデータ）の流れを示している。
材料組織予測装置１は、核生成順位決定部１０１、初期温度計算部１０２、温度計算部１０３、核生成回数計算部１０４、核生成部１０５、粒成長計算部１０６、変態発熱計算部１０７、時間ステップ更新部１０８、最終ステップ判定部１０９、形状特性値計算部１１０、統計量計算部１１１、材質計算部１１２、許容範囲判定部１１３、判定結果出力部１１４、下工程対処部１１５、温度変化曲線設定部１１６、温度変化曲線修正部１１７、統計量入力部１１８、核生成頻度式修正部１１９、及び温度変更部１２０を備えている。
まず、作業者は、金属材サンプルから母相の材料組織に関する情報（以下、「母相組織情報」という。）を取得する。母相組織情報の取得方法の具体的な例について、以下に説明する。
金属材サンプルは、材料組織を予測する金属材ＭＴと合金成分及び上工程のプロセス（温度条件又は加工率）がほぼ同一条件の金属である。作業者は、この金属材サンプルを母相のまま急冷（組織凍結）し、研磨し、エッチングする。このように処理された金属材サンプルを、顕微鏡で観察すると、母相の粒界が際立った画像が得られる。金属材ＭＴが鉄鋼材料の場合には、ピクリン酸溶液で腐食することにより、母相であるオーステナイト組織の粒界を明瞭に観察することができる。一般にエッチング後の画像では、傾角（隣接する結晶間の結晶方位の差）の小さい粒界ほど細く薄く現れる。目視により又は画像解析ソフトウェアを用いて、この画像に現れる粒界の太さ及び濃さを検出することで傾角を知ることができる。
また、金属材サンプルを冷却後、研磨して、走査型電子顕微鏡を応用したＥＢＳＤ(electron back scattering diffraction)装置で観察した画像から、析出相の各結晶粒の結晶方位が得られる。母相の結晶方位と析出相の結晶方位には一定の関係がある。この関係を利用して母相の結晶方位と粒界を再構築できる。従って、ＥＢＳＤ装置で観察した画像から粒界の傾角を把握してもよい。
このようにして、作業者は、母相組織情報を取得する。このように画像から母相組織情報を取得する行為は、作業者が人為的に行ってもよいし、材料組織予測装置１がプログラムにより実行してもよい。
核生成順位決定部１０１は、母相組織情報に基づいて金属材ＭＴ内の各計算対象領域の各計算点の核生成順位を決定する。核生成順位は、母相からの析出相の核生成の順序を表す。ここでは、核生成順位が小さい程、先に核生成されるものとする。多結晶金属では、複数の粒界が交差する部分又は傾角の大きい粒界の部分で核生成が起こり易い。従って、このような部分の計算点の核生成順位を高くするように、各計算点の核生成順位を決定する。
図４及び図５を参照して、核生成順位決定部１０１による核生成順位の決定方法について説明する。図４は、核生成順位決定部１０１により各計算点をレベル分けした金属材サンプルの模式図である。図５は、核生成順位決定部１０１により各計算点の核生成順位を決定した金属材サンプルの模式図である。図４及び図５中の破線Ｌ１，Ｌ２，Ｌ３，Ｌ４は、母相の粒界の位置を示している。粒界の各部分は、最も近い計算点に属するものとする。
核生成順位決定部１０１は、各粒界Ｌ１〜Ｌ４についてレベルを決定する。粒界Ｌ１〜Ｌ４は、金属材サンプルの画像における太さ及び濃さに応じて、レベル１からレベル４の４段階に分ける。粒界Ｌ１〜Ｌ４は、画像における粒界が太く濃いほど、核生成の優先度の高いレベルが付けられる。太さ又は濃さを判断する基準は、予め決められた絶対的な基準でもよいし、全ての粒界の中での位置付けをするための相対的な基準でもよい。ここでは、粒界Ｌ１〜Ｌ４は、太く濃い粒界をレベル４、やや太く濃い粒界をレベル３、通常の粒界をレベル２、細く薄い粒界をレベル１のようにレベル分けする。
次に、レベル分けした粒界Ｌ１〜Ｌ４に基づいて、各計算点のレベル分けをする。図４の升目（計算点）の中央の数値は、その計算点のレベルを示している。各計算点は、次のようにレベル分けする。
レベル４の粒界Ｌ４が交差する計算点をレベル４、レベル３の粒界Ｌ３が交差する計算点をレベル３、レベル２の粒界Ｌ２が交差する計算点をレベル２、レベル１の粒界Ｌ１が交差する計算点をレベル１とする。粒界Ｌ１〜Ｌ４が交差しない計算点を、核生成順位が最も低くなるレベル０とする。複数の粒界Ｌ１〜Ｌ４が交差する計算点（即ち、トリプルジャンクションに相当する計算点）を核生成順位が最も高くなるレベル５とする。
核生成順位決定部１０１は、各計算点のレベル分け後、各計算点の核生成順位を決定する。図５の升目の中央の数値は、その計算点の各生成順位を示している。各計算点は、次のように核生成順位を決定する。
核生成順位決定部１０１は、レベルの高い計算点から順に核生成順位を振る。従って、レベル５の計算点には、最も小さい番号が振られる。即ち、これらの計算点では、早期に核生成がされる。また、レベル０の計算点には、最も大きい番号が振られる。また、同一のレベルの複数の計算点に対しては、作業者の意図しない偏りのある選択又は画像解析ソフトウェアによる走査順による影響を避けるために、乱数で核生成順位を決める。
初期温度計算部１０２は、金属材ＭＴの各計算対象領域の各計算点の初期温度を計算する。初期温度計算部１０２は、放射温度計１５による測定値を用いて、初期温度を計算する。なお、初期温度の計算は、接触式温度計による測定値又は上工程の実績値などに基づく予測値を用いてもよい。初期温度計算部１０２は、計算した初期温度を補正してもよい。例えば、初期温度計算部１０２は、金属材ＭＴの温度測定位置と計算対象領域が離れている場合、又は、放射温度計１５による測定後時間が経過した場合には、金属材ＭＴの表面からの熱放射、空気の対流、又は材料内部の熱伝導などを考慮して補正する。なお、計算対象領域の大きさが温度勾配に比べて十分に小さく、同一の計算対象領域内の複数の計算点の温度差が殆ど無視できる場合には、その計算対象領域内の全ての計算点の初期温度を同一値としてもよい。このように初期温度を決定しても、計算精度の低下は限定的であるため、実用上は、この計算精度の低下を考慮して運用すればよい。
材料組織予測装置１は、核生成順位決定部１０１及び初期温度計算部１０２による計算後、時間ステップ毎に、温度計算部１０３、核生成回数計算部１０４、核生成部１０５、粒成長計算部１０６、変態発熱計算部１０７、時間ステップ更新部１０８、及び最終ステップ判定部１０９による計算を行う。材料組織予測装置１は、時間ステップ毎に実行されるこれらの一連の計算を、金属材ＭＴの冷却（又は加熱）の開始から終了まで計算する。
温度計算部１０３は、金属材ＭＴの各計算対象領域の各計算点の温度を時間ステップ毎に計算する。金属材ＭＴの冷却（又は加熱）が実際に開始される前は、温度計算部１０３は、目標とする温度変化曲線を用いて、各計算点の温度を計算する。金属材ＭＴの冷却（又は加熱）が実際に行われている間は、温度計算部１０３は、計算時点までの実績の温度変化曲線及び計算時点以降の目標とする温度変化曲線を用いて、各計算点の温度を計算する。金属材ＭＴの冷却（又は加熱）が実際に行われた後は、温度計算部１０３は、実績の温度変化曲線を用いて、各計算点の温度を計算する。例えば、実績の温度変化曲線は、冷却水流量又は誘導加熱装置の電流などの実績値に基づいて、差分法を用いた温度モデルにより計算される。
核生成回数計算部１０４は、母相と析出相の自由エネルギー差及び拡散係数に基づいて、計算対象領域内での核生成回数ΔＮαを計算する。まず、核生成回数計算部１０４は、母相からの析出相の核生成頻度Ｉを計算する。核生成頻度Ｉは、例えば、次式で表される。
Figure 0005919384
ここで、Ｒはガス定数である。Ｔは絶対温度である。Ｄγは変態に伴い拡散する元素の母相における拡散係数であり、元素濃度と温度の関数である。ΔＧｖは、母相と析出相の自由エネルギー差であり、熱力学データベースによる値を利用した値である。ｋ１は界面エネルギー及び転位密度に依存する係数である。ｋ２は活性化エネルギーに依存する係数である。
核生成回数計算部１０４は、計算した核生成頻度Ｉを用いて、ある時間ステップにおいて計算対象領域内で核生成させる回数ΔＮαを、次式により計算する。回数ΔＮαは、整数であるため、計算結果の少数部を残差ｒｅｓとして記憶する。記憶した残差ｒｅｓは、次の時間ステップの計算に加える。
Figure 0005919384
Figure 0005919384
ここで、Δtは時間増分である。関数ＴＲＵＮＣは、小数点以下を切り捨てる関数である。Ｖγは、計算対象領域内の未変態部の体積（２次元計算の場合は面積）である。未変態部か否かは、各計算点のスカラー量φｉが全ての結晶粒（ｉ＝１〜ｎ）について０になっているか否かで判断する。
核生成部１０５は、時間ステップ毎に、計算対象領域内に核生成回数計算部１０４により計算された核生成回数ΔＮαに従って、ΔＮα個の新たな析出相の核を生成する。まず、核生成部１０５は、核生成順位の高い方から計算点を探索する。核生成部１０５は、探索した計算点が析出相に含まれていないことを確認する。計算点が析出相に含まれていないことを確認するには、計算点のスカラー量φiが全ての結晶粒（ｉ＝１〜ｎ）について０になっていることを確認すればよい。探索した計算点が、他の析出相の結晶粒に既に含まれていた場合は、核生成順位が次に高い計算点を探索する。これを繰り返して、核生成部１０５は、析出相に含まれていない計算点を決定する。核生成部１０５は、決定した計算点に新たな結晶核を生成する。
具体的には、次のような演算処理をする。核生成部１０５は、結晶粒の総数ｎを１つ増やし、増やした結晶粒についてのスカラー量φｎの記憶領域を新たに割り付ける。核生成部１０５は、結晶核が生成された計算点の新に生成された結晶粒についてのスカラー量φｎを１にし、それ以外の計算点の新に生成された結晶粒についてのスカラー量φｎを０とする。
粒成長計算部１０６は、フェーズフィールド法を用いて、析出相の粒成長による析出相の界面の移動を計算する。具体的には、粒成長計算部１０６は、全ての計算点のスカラー量φｉを次式により更新する。
Figure 0005919384
Figure 0005919384
ここで、ｎは結晶粒の総数である。ｉ，ｊ，ｋは、結晶粒の番号である。結晶粒の番号が１のときは、母相を示している。ηは計算点間隔である。Ｇｉ,Ｇｊは、熱力学データベースによる値を利用した結晶粒の自由エネルギーを表している。μは界面の易動度（モビリティ）である。σは界面エネルギーである。μ，σは材料の合金組成による定数である。
また、粒成長計算部１０６は、変態に伴い拡散する元素の濃度cを次の拡散方程式により更新する。
Figure 0005919384
Figure 0005919384
ここで、γは母相、αは析出相をそれぞれ示す。Ｄαは変態に伴い拡散する元素の析出相における拡散係数である。∇はベクトル微分演算子である。
変態発熱計算部１０７は、ギブス−ヘルムホルツの式により、結晶粒の自由エネルギーから温度Ｔにおける母相のエンタルピＨα及び析出相のエンタルピＨγを次式により計算する。
Figure 0005919384
Figure 0005919384
次に、変態発熱計算部１０７は、粒成長による析出相界面の移動の計算結果に基づいて、単位重量あたりの変態発熱量Ｑを次式により計算する。
Figure 0005919384
ここで、ｗｉは計算点ｉの占める部分の重量である。Ｍはモル質量である。
変態発熱計算部１０７は、計算した単位重量あたりの変態発熱量Ｑを温度制御装置５に伝送する。温度制御装置５は、変態発熱計算部１０７により計算された変態発熱量Ｑに基づいて、目標とする温度変化曲線に従って金属材ＭＴを冷却するように、薄板熱間圧延ライン２の水冷装置１２を制御する。
変態発熱計算部１０７による演算処理終了後、時間ステップ更新部１０８は、時間ステップを更新する。更新による時間ステップの増分は、上述の拡散方程式による界面の移動量が、計算点の間隔を超えない範囲で行う。
最終ステップ判定部１０９は、終了条件を満たしているか否かの判定を行う。終了条件を満たしていない場合、材料組織予測装置１は、温度計算部１０３による演算処理に戻って、上述の一連の演算処理を繰り返す。これにより、金属材ＭＴの計算対象領域内の材料組織の変化がシミュレーションされる。
最終ステップ判定部１０９による終了条件を満たした場合、材料組織予測装置１は、以降の演算処理を実行する。なお、以降に説明する演算処理についても、上述した繰り返し行われる一連の演算処理に含めてもよい。
形状特性値計算部１１０は、各結晶粒について、スカラー量φに基づいて、各種形状特性値を計算する。形状特性値は、粒面積、面積率、粒度、等価円直径、等価楕円の長軸長さ、等価楕円の短軸長さ、等価楕円の長短軸比、周囲長、真円度、粒体積、体積率、等価球直径、等価楕円体の長軸長さ、等価楕円体の短軸長さ、等価楕円体の長短軸比、又は表面積などである。形状特性値計算部１１０は、例えば、番号がｉの結晶粒のスカラー量φiを計算対象領域全体に亘って積算し、粒面積Ａｉ（２次元計算の場合）又は粒体積Ｖｉ（３次元計算の場合）を計算する。
Figure 0005919384
Figure 0005919384
形状特性値計算部１１０は、上式により計算した粒面積Ａｉ又は粒体積Ｖｉに基づいて、次式により等価円直径ｄｉを計算する。
Figure 0005919384
Figure 0005919384
統計量計算部１１１は、形状特性値計算部１１０により計算された複数の結晶粒の形状特性値の統計量を計算する。統計量は、形状特性値の度数、最大値、最小値、分散値、平均値、粒面積に基づく加重平均値、粒体積に基づく加重平均値、標準偏差値、ガンマ分布係数、ワイブル分布係数、切断法に基づく平均粒径、比較法に基づく平均粒径、又はＡＳＴＭ(American Society for Testing and Materials)規格に基づく粒度番号などである。
２次元計算では、統計量計算部１１１は、粒面積に基づく重み付き平均粒径を次式により計算する。
Figure 0005919384
ここで、ｎは結晶粒の数である。
３次元計算では、統計量計算部１１１は、粒体積に基づく重み付き平均粒径を次式により計算する。
Figure 0005919384
２次元計算では、統計量計算部１１１は、粒面積に基づく重み付き粒径標準偏差を次式により計算する。
Figure 0005919384
３次元計算では、統計量計算部１１１は、粒体積に基づく重み付き粒径標準偏差を次式により計算する。
Figure 0005919384
材質計算部１１２は、統計量計算部１１１により計算された形状特性値の統計量に基づいて、金属材ＭＴの材質を計算する。計算された材質は、冷却後（又は加熱後）の金属材ＭＴの予測される材質となる。材質は、降伏応力（ＹＳ）、引張強度（ＴＳ）、伸び（ＥＬ）、硬さ、靭性、穴拡げ性、又は磁性などである。
例えば、材質計算部１１２は、Ｈａｌｌ−Ｐｅｔｃｈ則により、降伏応力ＹＳを次式により計算する。
Figure 0005919384
ここで、ｋは実験などで予め決められた定数である。
許容範囲判定部１１３は、材質計算部１１２により計算された金属材ＭＴの材質が予め決められた許容範囲内か否かを判定する。許容範囲判定部１１３は、判定結果を判定結果出力部１１４に出力する。
判定結果出力部１１４は、許容範囲判定部１１３による判定結果を表示画面又は記憶媒体などに出力する。判定結果出力部１１４は、金属材ＭＴの材質が許容範囲を外れた場合、例えば、次のような演算処理をする。判定結果出力部１１４は、許容範囲を外れた計算対象領域の金属材ＭＴにおける位置、又は許容範囲を外れた結晶粒の計算対象領域における位置などを表示画面又は記憶媒体に出力する。判定結果出力部１１４は、金属材ＭＴの材質が許容範囲を外れたことを作業者に知らせるための警報（画面表示又は警報音）を出力する。判定結果出力部１１４は、金属材ＭＴの許容範囲を外れた部分について、下工程で対処するために、下工程対処部１１５に必要な情報を出力する。
下工程対処部１１５は、許容範囲判定部１１３により許容範囲外と判定された金属材ＭＴについて、下工程で対処するための演算処理をする。ここで、下工程とは、例えば、熱間薄板圧延ラインの場合、切断、溶接、冷間圧延、酸洗、レベラー、又はスキンパスなどである。下工程対処部１１５は、薄板熱間圧延ライン２に設置された光学式金属検知器１６又は搬送ローラ１４の回転計などにより、金属材ＭＴの許容範囲外と判定された計算対象領域を追跡する。下工程対処部１１５は、金属材ＭＴの許容範囲外の部分を処理させるために、処理を担当する下工程に追跡情報などを伝送する。例えば、下工程が切断装置による切断の場合、切断装置は、金属材ＭＴの許容範囲外の部分が到達すると、該当部分の切断をする。これにより、金属材ＭＴから材質が許容範囲外の部分が取り除かれる。
次に、許容範囲判定部１１３により金属材ＭＴの材質が許容範囲を外れたと判定された場合に、金属材ＭＴに与える温度を変更する方法について説明する。
温度変更部１２０は、許容範囲判定部１１３により金属材ＭＴの材質が許容範囲外と判定された場合、統計量計算部１１１により計算された析出相の結晶粒の形状特性値を統計処理した結果が許容範囲内になるように、許容範囲外とされた計算対象領域の温度変化曲線又は入出熱量を修正するための計算をする。
例えば、統計量計算部１１１により計算された粒径標準偏差が許容範囲を外れた場合、温度変更部１２０は、温度制御装置５の制御に用いる温度変化曲線を折れ線で表した場合のある節点ｐの温度Ｔｐを次式により変更する。
Figure 0005919384
ここで、σｄＡＩＭは、粒径標準偏差の目標値である。例えば、σｄＡＩＭは、許容範囲の中央値である。
温度変更部１２０は、温度Ｔｐに微小値±ΔＴを加えた場合の粒径標準偏差の２つの計算値σｄ（Ｔｐ＋ΔＴ），σｄ（Ｔｐ−ΔＴ）を試算する。温度変更部１２０は、試算結果に基づいて影響係数を次式により計算する。
Figure 0005919384
温度Ｔｐが操業上の上限値Ｔｐｍａｘを超える場合、温度変更部１２０は、温度Ｔｐを上限値Ｔｐｍａｘにする。温度Ｔｐが操業上の下限値Ｔｐｍｉｎを下回る場合、温度変更部１２０は、温度Ｔｐを下限値Ｔｐｍｉｎにする。温度変更部１２０は、計算した節点pの温度Ｔｐを温度変化曲線修正部１１７に出力する。
温度変化曲線修正部１１７は、温度変化曲線設定部１１６に設定されている温度変化曲線の節点ｐの温度を、温度変更部１２０により計算された温度Ｔｐに修正する。これにより、粒径標準偏差が許容範囲内となるように温度変化曲線が修正される。
次に、核生成頻度式の係数を修正する方法について説明する。
まず、作業者は、冷却（又は加熱）が完了後の金属材ＭＴの形状特性値の統計量の実績値を次のように測定する。
作業者は、冷却（又は加熱）が完了後の金属材ＭＴの顕微鏡写真又は電子顕微鏡写真を撮影する。作業者は、撮影した写真を画像解析して、形状特性値の統計量の実績値を得る。例えば、鉄鋼の試験片を研磨後、フェライト組織をナイタール液で腐食すると、フェライト結晶粒の粒界を明瞭に観察することができる。作業者は、このように処理された鉄鋼の試験片を光学顕微鏡又は走査型電子顕微で写真に撮る。撮影した画像を画像処理装置に読み込ませて、輝度調整後、濃淡の適切な閾値で二値化することで粒界が抽出される。その後、作業者は、各結晶粒の画素数をカウントして、面積に換算した後集計する。これにより、平均粒径又は粒径標準偏差などの統計量が得られる。その他に、平均粒径又は粒径標準偏差などの統計量は、ＥＢＳＤ装置により得られる。また、作業者は、金属材ＭＴに縦波超音波を繰り返し反射させる。反射させた検出波形の減衰特性から平均粒径などの統計量が得られる。
作業者は、このようにして得られた形状特性値の統計量の実績値を、材料組織予測装置１の統計量入力部１１８に入力する。材料組織予測装置１への実績値の入力方法は、画像処理装置から伝送してもよいし、記憶媒体を介して入力してもよい。
統計量入力部１１８には、統計量計算部１１１により計算された統計量の予測値及び測定された統計量の実績値が入力される。統計量入力部１１８は、統計量の予測値及び統計量の実績値を核生成頻度式修正部１１９に出力する。
核生成頻度式修正部１１９は、統計量の予測値及び統計量の実績値に基づいて、次のように核生成頻度式の係数を修正する。
核生成頻度式修正部１１９は、統計量の予測値σｄと統計量σｄＡＣＴの実績値の偏差を計算する。核生成頻度式修正部１１９は、この偏差を小さくするように、核生成回数計算部１０４で用いる式（１）に示す核生成頻度式の係数ｋ１，ｋ２を修正する。
係数ｋ１の修正は、次式により行う。
Figure 0005919384
核生成頻度式修正部１１９は、係数ｋ１に微小値±Δｋ１を加えた場合の粒径標準偏差σｄの２つの計算値σｄ（ｋ１＋Δｋ１），σｄ（ｋ１−Δｋ１）を試算する。核生成頻度式修正部１１９は、試算結果に基づいて影響係数を次式により計算する。
Figure 0005919384
係数ｋ１が操業上の上限値ｋ１ｍａｘを超える場合、核生成頻度式修正部１１９は、係数ｋ１を上限値ｋ１ｍａｘにする。係数ｋ１が操業上の下限値ｋ１ｍｉｎを下回る場合、核生成頻度式修正部１１９は、係数ｋ１を下限値ｋ１ｍｉｎにする。
このように修正された係数ｋ１，ｋ２を材料の合金組成の違いに応じて層別に記憶させてもよい。例えば、鉄鋼材料の場合であれば、炭素含有量、ニオブ含有量、クロム含有量などに応じて、極低炭素鋼、低炭素鋼、中炭素鋼、高炭素鋼、ニオブ鋼、クロム鋼のように層別に区分する。このようにすることで、合金組成の違いによる影響を補正することができる。
ここで、結晶粒径のバラツキ度合いは、核生成挙動との相関が強い。例えば、初期に核生成した結晶粒は十分に粒成長し、結晶粒径が大きく成り易い。このため、核生成頻度式の係数の算定には、形状特性値として、結晶粒径の標準偏差、ガンマ分布係数、ワイブル分布係数、又はＡＳＴＭ粒度番号などの結晶粒径の分布に関する統計量を用いるのが適している。
ここでは、係数ｋ１を修正する例を示したが、係数ｋ２を修正する場合も同様である。また、核生成頻度式として式（１）を用いたが、他の式でも同様の手順で核生成頻度式の係数を修正することができる。
本実施形態によれば、以下の作用効果を得ることができる。
従来のフェーズフィールド法では、冷却過程（又は加熱過程など）のシミュレーションで、妥当なタイミングで妥当な数量の固相変態の起点となる核を生成させることができていないために、定量的な予測精度が不十分であった。
これに対して、本実施形態に係る材料組織予測装置１では、金属材ＭＴの核生成及び粒成長を組み合わせて計算することで、母相組織の粒界付近などの適正な位置に、適正な数量の核を適正なタイミングで生成させることができ、温度変化により相変態をする金属材ＭＴの組織を定量的に予測することができる。
従って、材料組織予測装置１は、金属材ＭＴの組織の定量的な予測結果に基づいて、金属材ＭＴの熱処理（冷却、加熱、又は等温過程など）による材質を予測することで、精度の高い材質の予測をすることができる。
また、材料組織予測装置１は、熱処理後の金属材ＭＴの材質の予測結果を評価することで、評価結果を製品の製造プロセスに反映させることができる。例えば、製品を製造するシステム（例えば、薄板熱間圧延ライン２）において、複数の評価結果を収集することで、以降に処理される金属材ＭＴの材質の品質を向上させることができる。また、製品を製造するシステムにおいて、熱処理された金属材ＭＴの材質の評価結果が許容範囲外であった場合、金属材ＭＴの許容範囲外と判定された部分を下工程で対処することで、その金属材ＭＴが加工された製品の材質を許容範囲内に修正することができる。
さらに、材料組織予測装置１は、粒成長の計算結果に基づいて、変態発熱量を計算する。薄板熱間圧延ライン２は、材料組織予測装置１により計算された変態発熱量に基づいて、金属材ＭＴに与える温度を制御することで、固相変態の挙動を反映した温度制御をすることができる。
また、材料組織予測装置１は、計算した析出相の結晶粒の形状特性値の統計量の予測値が、熱処理後の析出相の結晶粒の形状特性値の統計量の実績値に一致するように、核生成頻度式の係数を修正する。核生成は確率的な現象であり、理論的に十分に解明されていない。このため、核生成頻度式にはある程度の誤差がある。従って、核生成頻度式の係数を実績値に基づいて修正することで、材料組織予測装置１による核生成挙動の計算精度を高めることができる。これにより、材料組織予測装置１による金属材ＭＴの組織の予測精度を高めることができる。
なお、実施形態において、予測の対象とする材料は、説明したものに限らない。例えば、相変態をする多結晶材料として、材料料、鉄鋼材料（溶質に炭素、マンガン、シリコン、ニオブ、ニッケル、アルミニウム、窒素などを含む）、磁性体材料（銅、鉄、白金、パラジウム、コバルト、クロム、ニッケル、ネオジウム、バリウム、ビスマス、サマリウムなどからなる合金や酸化物）、誘電体、又は、ポリスチレン−ポリメタクリル酸メチルといったブロック共重合体のような高分子有機材料などに代表される機能性材料などでもよい。例えば、磁性体材料では、加熱又は冷却などの熱処理（等温過程、アニール及び時効処理なども含む）によって、核形成もしくは相分離を起こし、結果さまざまに微細組織が形成される。また、誘電体及び高分子材料も、同様に、加熱又は冷却による熱処理によってさまざまな組織形態を形成する。従って、実施形態と同様の構成で、これらの材料の熱処理中の相変態又は相分離によって生成される材料組織を定量的に予測してもよい。例えば、材料組織予測装置による予測方法は、有機材料の結晶化プロセスに適用してもよい。
実施形態において、予測する相変態の母相及び析出相は、任意の組合せを選択してよい。例えば、鉄鋼材料では、フェライト相、オーステナイト相、パーライト相、マルテンサイト相、ベイナイト相、セメンタイト相の任意の組合せの相変態による材料組織の予測をしてもよい。また、鉄鋼材料に限らず、他の材料においても同様に、任意の組合せの相変態による材料組織の予測をしてもよい。
実施形態において、核生成順位の決定方法は、説明したものに限らない。例えば、材料組織予測装置１に、典型的な事例について核生成順位が決定されたサンプル情報を記憶媒体に予め記憶させておいてもよい。材料組織予測装置１は、核生成順位の決定時に、予測対象の金属材ＭＴと合金成分及び上工程のプロセス（温度又は加工率など）が類似している事例のサンプル情報を検索して、サンプル情報に振られた核生成順位に従って、核生成順位を決定してもよい。このように、核生成順位を決定することで、材料組織予測装置１による計算量又は作業者による作業量を軽減することができる。
実施形態では、析出相の組織の時間発展を演算することで、材料組織の予測をしたが、析出相と同様に、母相の組織の時間発展を演算してもよい。
本実施形態に係る材料組織の予測方法は、製造システム（例えば、薄板熱間圧延ライン）に限らず、どのようなシステム又は装置に組み込まれていてもよい。また、材料組織予測装置１は、実施形態で説明した構成に限らない。他の装置にある任意の組合せの機能が材料組織予測装置１にあってもよいし、材料組織予測装置１の一部の機能が他の装置にあってもよい。
材料組織予測装置を構成するコンピュータはどのように構成されていてもよい。例えば、コンピュータは、各種プロセッサなどの演算部、メモリ又はハードディスクなどの記憶部を備える。この他に、コンピュータは、入出力部及び表示部などを備えていてもよい。また、コンピュータは、マイクロコンピュータ又はパーソナルコンピュータなどの種別を問わない。
実施形態において、材料組織予測装置（コンピュータ）で演算していた処理の任意の一部の処理を、人が行ってもよい。実施形態に係る材料組織の予測方法に基づく原理及び理論は、自然法則を利用するものである。従って、材料組織の予測方法の一部の処理を人が行ったとしても、実施形態に係る材料組織の予測方法は、全体として自然法則を利用するものである。
なお、本発明は上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態にわたる構成要素を適宜組み合わせてもよい。

Claims (20)

  1. 温度変化により母相から析出相に相変態する材料の組織を予測する材料組織予測装置であって、
    前記材料に前記温度変化を与えるための温度条件に基づいて、前記材料の計算対象領域内の複数の計算点の温度を計算する温度計算手段と、
    前記温度計算手段により計算された前記複数の計算点の温度に基づいて、前記計算対象領域内の核生成回数を計算する核生成回数計算手段と、
    前記核生成回数計算手段により計算された前記核生成回数及び予め決められた前記複数の計算点の核生成の起こり易さに基づいて、前記複数の計算点から前記析出相の核を生成する析出相生成点を決定する析出相生成点決定手段と、
    前記析出相生成点決定手段により前記析出相生成点が繰り返し決定される度に、前記析出相生成点における前記析出相の粒成長を計算する粒成長計算手段と、
    前記粒成長計算手段により計算された前記析出相の粒成長に基づいて、前記材料の組織を予測する材料組織予測手段と
    を備えることを特徴とする材料組織予測装置。
  2. 前記粒成長計算手段により計算された前記析出相の粒成長に基づいて、前記析出相の結晶粒の形状特性値を計算する形状特性値計算手段と、
    前記形状特性値計算手段により計算された前記形状特性値に基づいて、複数の前記形状特性値の統計量を計算する統計量計算手段と、
    前記統計量計算手段により計算された前記統計量に基づいて、前記材料の材質を予測する材質予測手段と
    を備えることを特徴とする請求項１に記載の材料組織予測装置。
  3. 冷却後又は加熱後の前記材料から得られる前記形状特性値の統計量である実績統計量を取得する実績統計量取得手段と、
    前記統計量計算手段により計算された前記統計量及び前記実績統計量取得手段により取得された前記実績統計量に基づいて、前記核生成回数計算手段により前記核生成回数を計算するための計算式を修正する核生成回数計算修正手段と
    を備えることを特徴とする請求項２に記載の材料組織予測装置。
  4. 前記材料組織予測手段により予測された前記材料の組織に基づいて、前記材料の材質を評価する材質評価手段
    を備えることを特徴とする請求項１から請求項３のいずれか１項に記載の材料組織予測装置。
  5. 前記核生成回数計算手段は、前記母相と前記析出相の自由エネルギー差及び前記母相の拡散係数に基づいて、前記核生成回数を計算すること
    を特徴とする請求項１から請求項４のいずれか１項に記載の材料組織予測装置。
  6. 前記複数の計算点に前記母相からの前記析出相の核生成の順位を決定する核生成順位決定手段を備え、
    前記析出相生成点決定手段は、前記核生成順位決定手段により決定された前記析出相の核生成の順位を用いて、前記析出相生成点を決定すること
    を特徴とする請求項１から請求項５のいずれか１項に記載の材料組織予測装置。
  7. 前記温度計算手段は、前記材料に既に加えられた前記温度条件の実績値に基づいて、前記複数の計算点の温度を計算すること
    を特徴とする請求項１から請求項６のいずれか１項に記載の材料組織予測装置。
  8. 温度変化により母相から析出相に相変態する材料から製品を製造するための製造方法であって、
    前記材料に前記温度変化を与えるための温度を制御し、
    前記温度を制御するための温度条件に基づいて、前記材料の計算対象領域内の複数の計算点の温度を計算し、
    計算した前記複数の計算点の温度に基づいて、前記計算対象領域内の核生成回数を計算し、
    計算した前記核生成回数及び予め決められた前記複数の計算点の核生成の起こり易さに基づいて、前記複数の計算点から前記析出相の核を生成する析出相生成点を決定し、
    前記析出相生成点が繰り返し決定される度に、前記析出相生成点における前記析出相の粒成長を計算し、
    計算した前記析出相の粒成長に基づいて、前記材料の組織を予測し、
    予測した前記材料の組織に基づいて、前記析出相の結晶粒の形状特性値を計算し、
    計算した前記形状特性値に基づいて、複数の前記形状特性値の統計量を計算し、
    計算した前記統計量を利用して、前記製品を製造すること
    を含むことを特徴とする製品製造方法。
  9. 計算した前記析出相の粒成長の計算結果に基づいて、前記析出相の前記相変態に伴う変態発熱量を計算することを含み、
    前記温度の制御は、前記温度条件及び計算した前記変態発熱量を用いること
    を特徴とする請求項８に記載の製品製造方法。
  10. 計算した前記統計量に基づいて、前記材料の材質を予測し、
    前記材質の予測結果に基づいて、前記製品を製造すること
    を含むことを特徴とする請求項８に記載の製品製造方法。
  11. 予測した前記材料の組織に基づいて、前記材料の材質を評価すること
    を含むことを特徴とする請求項８に記載の製品製造方法。
  12. 前記材料の材質の評価結果に基づいて、下工程で前記材料を処理すること
    を含むことを特徴とする請求項１１に記載の製品製造方法。
  13. 前記材料の材質の評価結果に基づいて、前記温度条件を変化させること
    を特徴とする請求項１１に記載の製品製造方法。
  14. 温度変化により母相から析出相に相変態する材料の組織を予測する材料組織予測方法であって、
    前記材料に前記温度変化を与えるための温度条件に基づいて、前記材料の計算対象領域内の複数の計算点の温度を計算し、
    計算した前記複数の計算点の温度を用いて前記計算対象領域内の核生成回数を計算し、
    計算した前記核生成回数及び予め決められた前記複数の計算点の核生成の起こり易さに基づいて、前記複数の計算点から前記析出相の核を生成する析出相生成点を決定し、
    前記析出相生成点が繰り返し決定される度に、前記析出相生成点における前記析出相の粒成長を計算し、
    計算した前記析出相の粒成長に基づいて、前記材料の組織を予測することを含むことを特徴とする材料組織予測方法。
  15. 計算した前記析出相の粒成長に基づいて、前記析出相の結晶粒の形状特性値を計算し、
    計算した前記形状特性値に基づいて、複数の前記形状特性値の統計量を計算し、
    計算した前記統計量に基づいて、前記材料の材質を予測すること
    を含むことを特徴とする請求項１４に記載の材料組織予測方法。
  16. 冷却後又は加熱後の前記材料から得られる前記形状特性値の統計量である実績統計量を取得し、
    計算した前記統計量及び取得した前記実績統計量に基づいて、前記核生成回数を計算するための計算式を修正すること
    を含むことを特徴とする請求項１５に記載の材料組織予測方法。
  17. 予測した前記材料の組織に基づいて、前記材料の材質を評価すること
    を含むことを特徴とする請求項１４から請求項１６のいずれか１項に記載の材料組織予測方法。
  18. 前記核生成回数の計算は、前記母相と前記析出相の自由エネルギー差及び前記母相の拡散係数を用いること
    を含むことを特徴とする請求項１４から請求項１７のいずれか１項に記載の材料組織予測方法。
  19. 前記複数の計算点に前記母相からの前記析出相の核生成の順位を決定し、
    前記析出相生成点の決定は、決定した前記順位を用いること
    を特徴とする請求項１４から請求項１８のいずれか１項に記載の材料組織予測方法。
  20. 前記複数の計算点の温度の計算は、前記材料に既に加えられた前記温度条件の実績値を用いること
    を特徴とする請求項１４から請求項１９のいずれか１項に記載の材料組織予測方法。
JP2014532693A 2012-08-31 2012-08-31 材料組織予測装置、製品製造方法及び材料組織予測方法 Active JP5919384B2 (ja)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2012/072234 WO2014033928A1 (ja) 2012-08-31 2012-08-31 材料組織予測装置、製品製造方法及び材料組織予測方法

Publications (2)

Publication Number Publication Date
JP5919384B2 true JP5919384B2 (ja) 2016-05-18
JPWO2014033928A1 JPWO2014033928A1 (ja) 2016-08-08

Family

ID=50182774

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2014532693A Active JP5919384B2 (ja) 2012-08-31 2012-08-31 材料組織予測装置、製品製造方法及び材料組織予測方法

Country Status (7)

Country Link
US (1) US9767227B2 (ja)
JP (1) JP5919384B2 (ja)
KR (1) KR101659702B1 (ja)
CN (1) CN104602830B (ja)
AU (1) AU2012388949B2 (ja)
IN (1) IN2015DN01600A (ja)
WO (1) WO2014033928A1 (ja)

Families Citing this family (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106794499B (zh) 2014-10-10 2018-10-12 杰富意钢铁株式会社 材料特性值推定方法、材料特性值推定装置、及钢带的制造方法
CN104820873B (zh) * 2015-05-13 2017-12-26 中国环境科学研究院 一种基于金属定量构效关系的淡水急性基准预测方法
JP6601762B2 (ja) * 2015-09-11 2019-11-06 株式会社日本製鋼所 鋼の熱処理シミュレーション方法および鋼の熱処理シミュレーションプログラム
DE102016100811A1 (de) 2015-09-25 2017-03-30 Sms Group Gmbh Verfahren und Ermittlung der Gefügebestandteile in einer Glühlinie
CN106053471A (zh) * 2016-08-04 2016-10-26 南京钢铁股份有限公司 一种真空热腐蚀条件下显现微合金钢奥氏体晶界的方法
JP6770236B2 (ja) * 2017-02-23 2020-10-14 富士通株式会社 熱伝導性の推定方法、熱伝導性の推定装置、及びプログラム
JP6747391B2 (ja) * 2017-06-30 2020-08-26 Ｊｆｅスチール株式会社 材料特性推定装置及び材料特性推定方法
JP7146377B2 (ja) * 2017-09-12 2022-10-04 株式会社Uacj 材料組織予測方法および材料組織予測装置
CN108460190B (zh) * 2018-02-06 2022-03-04 北京卫星制造厂有限公司 一种预测蒙皮法兰结构焊缝预置应力焊接的方法
CN109187610B (zh) * 2018-07-20 2022-03-18 中冶连铸技术工程有限责任公司 铸坯微观组织模拟方法
CN108890529B (zh) * 2018-07-25 2023-06-23 浙江工业大学 光催化钴基合金加工控制系统及控制方法
CN108959832B (zh) * 2018-09-25 2022-12-02 武汉理工大学 用于优化M50NiL钢渗碳淬火工艺的晶粒长大预测方法
JP7147645B2 (ja) * 2019-03-18 2022-10-05 日本製鉄株式会社 結晶粒解析装置、結晶粒解析方法、及びプログラム
CN111553032B (zh) * 2020-04-27 2023-09-15 湖北文理学院 整体叶盘铣削温度预测方法、装置、设备及存储介质
WO2022059186A1 (ja) * 2020-09-18 2022-03-24 株式会社東芝 粒度推定装置、粒度推定方法、粒度推定プログラム、粒度推定システム。
CN113063650B (zh) * 2021-03-16 2023-02-28 天津市新天钢钢铁集团有限公司 在20℃～-40℃温度条件下钢板扩孔试验方法
CN114235864B (zh) * 2021-11-11 2024-05-10 泛锐云智科技(郑州)有限公司 一种用于ebsd测试的铋合金试样的制备方法
CN117219213B (zh) * 2023-11-07 2024-01-23 中北大学 相场法定量预测含铌奥氏体不锈钢铌化物析出规律的方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0587800A (ja) * 1991-05-23 1993-04-06 Nippon Steel Corp 鋼板の材質予測方法
JPH05142126A (ja) * 1991-11-21 1993-06-08 Nippon Steel Corp 鋼板の材質予測方法
JPH06264137A (ja) * 1993-03-11 1994-09-20 Nippon Steel Corp 多結晶材料の材料組織制御条件決定法
JP2012032404A (ja) * 2006-06-20 2012-02-16 Toshiba Mitsubishi-Electric Industrial System Corp 組織材質測定システム

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0541825A4 (en) * 1991-06-04 1995-10-11 Nippon Steel Corp Method of estimating material of steel product
CN1201880C (zh) * 2002-01-11 2005-05-18 中国科学院金属研究所 一种热轧过程带钢组织演变与性能预测的方法
JP4402502B2 (ja) 2004-04-13 2010-01-20 東芝三菱電機産業システム株式会社 巻取温度制御装置
JP2007091543A (ja) 2005-09-29 2007-04-12 Canon Inc グレイン形状計算方法
CN101591729B (zh) * 2009-06-19 2011-01-26 东北大学 预测板带钢热变形中奥氏体动态再结晶组织演变的方法
CN102234716A (zh) * 2011-08-03 2011-11-09 天津大学 纯铜热变形过程中动态再结晶组织演变的预测方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0587800A (ja) * 1991-05-23 1993-04-06 Nippon Steel Corp 鋼板の材質予測方法
JPH05142126A (ja) * 1991-11-21 1993-06-08 Nippon Steel Corp 鋼板の材質予測方法
JPH06264137A (ja) * 1993-03-11 1994-09-20 Nippon Steel Corp 多結晶材料の材料組織制御条件決定法
JP2012032404A (ja) * 2006-06-20 2012-02-16 Toshiba Mitsubishi-Electric Industrial System Corp 組織材質測定システム

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
JPN6015049836; 村松眞由、味岡秀恭、岡田敏、只野裕一、志澤一之: '大角粒界を有するサブグレインからの核生成と再結晶に関する多結晶Ｐｈａｓｅ-Ｆｉｅｌｄシミュレーション' 理論応用力学講演会　講演論文集 第５６回, 20070718, ２３０, 日本学術会議「機械工学委員会・土木工学・建築学委 *
JPN6015049837; 山本隆史、山中晃徳、高木知弘、冨田佳宏: 'Ｍｕｌｔｉ-Ｐｈａｓｅ-Ｆｉｅｌｄ法による鉄鋼材料の拡散相変態シミュレーション' 計算力学講演会講演論文集 第２２回, 20091010, ｐ４９-５０, 日本機械学会 *

Also Published As

Publication number Publication date
CN104602830A (zh) 2015-05-06
IN2015DN01600A (ja) 2015-07-03
AU2012388949B2 (en) 2016-11-10
KR20150038317A (ko) 2015-04-08
US9767227B2 (en) 2017-09-19
AU2012388949A1 (en) 2015-03-19
WO2014033928A1 (ja) 2014-03-06
KR101659702B1 (ko) 2016-09-26
CN104602830B (zh) 2017-03-01
US20150178415A1 (en) 2015-06-25
JPWO2014033928A1 (ja) 2016-08-08

Similar Documents

Publication Publication Date Title
JP5919384B2 (ja) 材料組織予測装置、製品製造方法及び材料組織予測方法
Samanta et al. Formation of bainite below the MS temperature: Kinetics and crystallography
Bok et al. Thermo-mechanical-metallurgical modeling for hot-press forming in consideration of the prior austenite deformation effect
Bernier et al. An alternative to the crystallographic reconstruction of austenite in steels
EP1608472B1 (en) A system for on-line property prediction for hot rolled coil in a hot strip mill
WO2014034964A2 (en) Method and system for predicting material structure
Pietrzyk et al. Conventional and multiscale modeling of microstructure evolution during laminar cooling of DP steel strips
Wang et al. A micromechanics-based machine learning model for evaluating the microstructure-dependent rolling contact fatigue performance of a martensitic steel
Kubota et al. Analysis of recrystallization behavior of hot-deformed austenite reconstructed from electron backscattering diffraction orientation maps of lath martensite
Perlade et al. Development of 3rd generation Medium Mn duplex steels for automotive applications
Pernach et al. Numerical modeling of phase transformation in dual phase (DP) steel after hot rolling and laminar cooling
Szeliga et al. Stochastic model describing evolution of microstructural parameters during hot rolling of steel plates and strips
Asgharzadeh et al. Cellular automata modeling of the kinetics of static recrystallization during the post-hydroforming annealing of steel tube
Bombač et al. Amelioration of surface cracking during hot rolling of AISI D2 tool steel
Razali et al. Practical approach for determining material parameters when predicting grain size after static recrystallization
Narayanan et al. The Causes and Effects of Preannealed Radiative Property Variations Across Full‐Hard Advanced High‐Strength Steel Coils
Ahmels et al. Manufacturing Induced Properties: Determination, Understanding, and Beneficial Use
Zhu Phase-field modeling of microstructure evolution in low-carbon steels during intercritical annealing
Bachniak et al. Selection of the optimization method for identification of phase transformation models for steels
Dehghani et al. Effect of thermomechanical processing on forming limit diagrams predicted by neural networks
Hartl et al. Machine learning driven prediction of mechanical properties of rolled aluminum and development of an in-situ quality control method based on electrical resistivity measurement
Cabus et al. The influence of grain morphology on texture measured after phase transformation in multiphase steels
Ayaz et al. Improvement of grain size and dome height of microalloyed steels using Taguchi method based on grey relational grade in controlled rolling
Hwang Inhomogeneity of microstructure and mechanical properties in compression-type bulk formed specimens
JP4144476B2 (ja) 鍛造焼ならしシミュレーション方法

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20160210

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: 20160315

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20160411

R150 Certificate of patent or registration of utility model

Ref document number: 5919384

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

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250