WO2018221353A1 - 材料組織計算装置および制御プログラム - Google Patents

Abstract

従来よりも正確に、工業プロセスにおけるＡｌの材料組織を予測する。情報処理装置（１）において、工程間情報統合部は、工程計算部に含まれる第ｉ工程計算部にＰＣ（ｉ）およびＭＳ（ｉ，０）を供給する。第ｉ工程計算部は、材料組織計算部にＭＳ（ｉ，ｔ）およびＴＭＰ（ｉ，ｔ）を供給することで、当該材料組織計算部にＭＳ（ｉ，ｔｆｉ）を算出させ、かつ、当該ＭＳ（ｉ，ｔｆｉ）を工程間情報統合部（１１）に供給する。工程間情報統合部（１１）は、第ｉ工程計算部から取得したＭＳ（ｉ，ｔｆｉ）を、ＭＳ（ｉ＋１，０）として設定する。

Description

材料組織計算装置および制御プログラム

　本発明は、工業プロセスによって製造されるＡｌ（アルミニウム）の材料組織を計算する材料組織計算装置に関する。

　近年、工業プロセスによって製造される金属材料の製造条件を検討するために、各製造工程における当該金属材料の特性の変化を予測する技術が提案されている。一例として、特許文献１には、Ａｌ合金板の製造工程において、当該Ａｌ合金板の材質を予測する方法が開示されている。

日本国公開特許公報「特開２００２－２２４７１号公報（２００２年８月１３日公開）」

　但し、以下に述べるように、工業プロセスにおいてＡｌの材料組織を予測（計算）するための方法については、改善の余地がある。本発明の一態様は、従来よりも正確に、工業プロセスにおけるＡｌの材料組織を予測することを目的とする。

　上記の課題を解決するために、本発明の一態様に係る材料組織計算装置は、第１工程から第Ｎ工程までを含む、アルミニウムを製造する工業プロセスについて、ｉを、１≦ｉ≦Ｎを満たす自然数、第ｉ工程の時刻ｔについて、ｔ＝０を当該第ｉ工程の開始時刻、ｔ＝ｔｆｉを当該第ｉ工程の終了時刻、第ｉ工程における上記アルミニウムの製造条件を示す情報をＰＣ（ｉ）、第ｉ工程の時刻ｔにおける上記アルミニウムの加工熱処理条件を示す情報をＴＭＰ（ｉ，ｔ）、第ｉ工程の時刻ｔにおける上記アルミニウムの材料組織を示す情報をＭＳ（ｉ，ｔ）、として、ＰＣ（１）～ＰＣ（Ｎ）およびＭＳ（１，０）を、あらかじめ設定された情報として取得する工程間情報統合部と、第１工程計算部から第Ｎ工程計算部までを含む工程計算部と、ＭＳ（ｉ，ｔ）およびＴＭＰ（ｉ，ｔ）に基づいて、ＭＳ（ｉ，ｔ）の時間変化を算出する材料組織計算部と、を備えており、上記材料組織計算部は、複数の冶金的現象のそれぞれの時間変化を算出する、複数の計算モジュールを含み、上記工程計算部における第ｉ工程計算部は、ＰＣ（ｉ）に基づいてＴＭＰ（ｉ，ｔ）を算出し、複数の上記計算モジュールのそれぞれには、共通のデータ構造を有するＭＳ（ｉ，ｔ）が供給され、ｉ＝１からｉ＝Ｎまでの順に、上記工程間情報統合部は、上記第ｉ工程計算部にＰＣ（ｉ）およびＭＳ（ｉ，０）を供給し、第ｉ工程において、ｔ＝０からｔ＝ｔｆｉまでの順に、上記第ｉ工程計算部は、（ｉ）上記材料組織計算部にＭＳ（ｉ，ｔ）およびＴＭＰ（ｉ，ｔ）を供給することで、当該材料組織計算部にＭＳ（ｉ，ｔｆｉ）を算出させ、かつ、（ｉｉ）当該ＭＳ（ｉ，ｔｆｉ）を上記工程間情報統合部に供給し、上記工程間情報統合部は、上記第ｉ工程計算部から取得したＭＳ（ｉ，ｔｆｉ）を、ＭＳ（ｉ＋１，０）として設定する。

　本発明の一態様に係る材料組織計算装置によれば、従来よりも正確に、工業プロセスにおけるＡｌの材料組織を予測することが可能となる。

実施形態１に係る情報処理装置の要部の構成を示す機能ブロック図である。 図１の情報処理装置における全体的な処理の流れを例示する図である。 図１の情報処理装置における工程間情報統合部の処理の流れを例示する図である。 図１の情報処理装置における第ｉ工程計算部の処理の流れを例示する図である。 図１の情報処理装置における材料組織計算部の処理の流れを例示する図である。 実施例におけるＰＣの項目を示す図である。 実施例におけるＭＳの項目を示す図である。 実施例における全体的な処理の流れを例示する図である。 実施例における計算結果を示す図である。 実施例、比較例、および各参考例のそれぞれについての、計算の処理の違いを示す図である。 実施例、比較例、および各参考例のそれぞれについての比較結果を示す図である。 比較例におけるＭＳの項目を示す図である。 比較例における全体的な処理の流れを例示する図である。 比較例における計算結果を示す図である。 参考例１における全体的な処理の流れを例示する図である。 参考例１における計算結果を示す図である。 参考例２における全体的な処理の流れを例示する図である。 参考例２における計算結果を示す図である。

　〔実施形態１〕
　以下、本発明の実施形態１について、図１～図１８に基づいて詳細に説明する。実施形態１の情報処理装置１（材料組織計算装置）は、Ａｌの最も初期的な材料組織を示す情報（後述するＭＳ（１，０））と、Ａｌの各製造工程における製造条件を示す情報（後述するＰＣ（１）～ＰＣ（Ｎ））とに基づいて、Ａｌの材料組織の時間的な変化を計算（予測）する。

　（用語の説明）
　はじめに、情報処理装置１の具体的な構成の説明に先立ち、本明細書において用いる各用語について述べる。

　「製造工程」とは、Ａｌ（Ａｌ製品）を製造するための任意の工程を含む。製造工程の一例としては、Ａｌ展伸材の加工および熱処理に関する工程（例：均質化処理、熱間圧延、冷間圧延、熱間押出、引抜き、熱間鍛造、冷間鍛造、溶体化処理、時効処理、および焼鈍等）が挙げられる。

　「製造条件」（Processing Condition，以下「ＰＣ」とも呼称する）とは、各製造工程において、Ａｌの品質（特性）に関連する任意の条件を含む。ＰＣには、各製造工程における設備の設定値が含まれる。一例として、冷間圧延の場合には、ＰＣには、開始板厚、終了板厚、加工度、圧延速度、圧延ワークロール径等のそれぞれの設定値が含まれる。ＰＣについては、後により詳細に述べる。

　「加工熱処理条件」（Thermo-Mechanical Processing，以下「ＴＭＰ」とも呼称する）とは、所定のＰＣのもとで、各時刻においてＡｌに与えられる温度、歪量、および歪速度等を意味する。本発明の一態様に係る材料組織計算装置では、ＰＣのデータセット（以下述べるＰＣ（ｉ））は、ＴＭＰを算出できるように設定されていることが必要である。ＴＭＰについても、後により詳細に述べる。

　ＴＭＰの計算には、ＦＥＭ（Finite Element Method，有限要素法）等の公知のシミュレーション技術が用いられてよい。一例として、熱間押出の場合には、Ａｌの合金成分、コンテナの温度、ビレット温度、押出形状，および押出速度等の設定値を用いて、ＦＥＭによりＴＭＰを算出できる。

　「材料組織」（Micro-Structure，以下「ＭＳ」とも呼称する）とは、微視的な材料組織（金属材料組織）を意味する。ＭＳには、例えば固溶元素、析出物、晶出物、転位、亜結晶粒、および結晶粒等が含まれる。ＭＳについても、後により詳細に述べる。

　「材料組織の変化」とは、ＭＳの変化を意味する。ＭＳの変化は、固溶、析出、回復、および再結晶等の現象（冶金的現象）に伴って生じる。ＭＳの変化には、固溶元素の変化、第二相粒子の変化、亜結晶粒の変化、および結晶粒の変化等が含まれる。

　以下に述べるように、情報処理装置１によれば、ある製造工程の開始時刻から終了時刻までのＭＳの時間的な変化を計算できる。そして、当該製造工程の終了時刻におけるＭＳの計算結果を、次に製造工程における初期的なＭＳとして用いることができる。これにより、複数の異なる製造工程によって製造されるＡｌについて、最終的なＭＳ（最後の製造工程が完了した時点のＭＳ）を算出できる。つまり、最終製品としてのＡｌの材質を予測できる。

　その結果、予測されたＡｌの材質に基づいてＰＣを設定できるので、最終製品の特性を効果的に制御することが可能となる。例えば、最終製品の強度バラツキの抑制、および、熱間押出における材料表面のひび割れの抑制等が可能となる。加えて、各製造工程において、ＭＳを要因とする不具合を低減できる。それゆえ、より高品質なＡｌを製造することが可能となる。

　情報処理装置１によって最終的なＭＳ（つまり材質）を予測することが可能なアルミ製品としては、アルミ板材（圧延材）、アルミ箔材、アルミ押出材、およびアルミ鍛造材等が挙げられる。一例として、これらのアルミ製品の製造工程には、以下の１～４に示す各工程が含まれる。

　＜１．アルミ板材（圧延材）の場合の製造工程の例＞
　溶解工程、脱ガス工程、鋳造工程、連続鋳造工程、均質化処理工程、熱間粗圧延工程、熱間仕上圧延工程、冷間圧延工程、溶体化処理工程、時効処理工程、矯正工程、焼鈍工程、および表面処理工程。

　＜２．アルミ箔材の場合の製造工程の例＞
　溶解工程、脱ガス工程、鋳造工程、連続鋳造工程、均質化処理工程、熱間粗圧延工程、熱間仕上圧延工程、冷間圧延工程、溶体化処理工程、時効処理工程、矯正工程、焼鈍工程、表面処理工程、および箔圧延工程。

　＜３．アルミ押出材の場合の製造工程の例＞
　溶解工程、脱ガス工程、鋳造工程、均質化処理工程、熱間押出工程、引き抜き工程、溶体化処理工程、時効処理工程、矯正工程、焼鈍工程、表面処理工程、および切断工程。

　＜４．アルミ鍛造材の場合の製造工程の例＞
　熱間鍛造工程（アルミ鋳物材、アルミ圧延材、またはアルミ押出材を素材とする）、冷間鍛造工程、溶体化処理工程、時効処理工程、および焼鈍工程。

　（情報処理装置１の目的）
　所望の特性を有するＡｌ製品を製造するためには、材料内部の微視的な金属組織（つまりＭＳ）を最適に制御する必要がある。ＭＳのサイズは、通常は５００μｍ以下のサイズである。また、上述の通り、ＭＳは、Ａｌ合金への元素の添加、および、製造工程中の加工歪または熱処理によって、材料内部において生じる冶金的現象（例：固溶、析出、晶出、塑性変形、回復、および再結晶等）により変化する。

　従来、ＭＳの制御に関しては多くの理論的研究が行われている。このため、上述の冶金的現象の多くは、所定の数式（理論式または実験式）により定式化されている。従って、実験室レベルでの単純な合金成分および製造工程を考慮する場合には、ＭＳ制御の方針（換言すれば、製造条件（ＰＣ））を決定しやすい。

　しかしながら、工業プロセス（工業製品）の場合には、Ａｌの合金成分および製造工程はいずれも複雑である。さらに、材料において複数の冶金的現象が同時に生じ、ＭＳの変化に影響を及ぼす。また、複数の冶金的現象は、相互作用的である。このため、工業プロセスの場合には、ＭＳ制御の方針を決定することが容易ではなかった。

　さらに、各製造工程におけるＭＳの変化は、前の製造工程におけるＭＳの影響を強く受ける。このため、最終製品のＭＳを制御するためには、全製造工程におけるＭＳの変化を制御することが必要である。

　しかしながら、従来の冶金理論では、単一の現象（冶金的現象）または単一の製造工程のみを考慮して、上述の定式化が行われている。このため、従来の冶金理論では、単一の現象または単一の製造工程におけるＭＳの変化しか計算することができず、工業プロセスにおけるＭＳを計算するには不十分であった。

　このため、工業プロセスでのＰＣを決定する場合には、トライアンドエラー（試作検討）によりＰＣを調整していく方法が主流である。しかしながら、試作検討によりＰＣを調整する場合には、多大なコスト（例：時間、労力、および試作費用）を要する。

　試作検討のコスト（回数）を減らすために、工業プロセスにおいて理論的な方法によりＭＳを予測する方法も提案されている。しかしながら、従来の方法では、ＭＳを変化させる複数の冶金的現象が同時に生じた場合に、当該冶金的現象の相互作用がＭＳの変化に及ぼす影響については十分に考慮されていない。

　また、従来の技術では、複数の製造工程に亘ってＭＳを予測するための方法も十分ではない。例えば、上述の特許文献１の技術では、Ａｌ合金板の材質を予測するために使用される材料組織因子の数が少ない。このため、特許文献１の技術は、複数の製造工程に亘るＭＳの変化の複雑性を取り扱うには不十分である。

　以上の点を鑑み、情報処理装置１は、工業プロセスにおける試作検討のコストを低減することを目的として、本願の発明者ら（以下、発明者ら）によって想到された。具体的には、以下に述べるように、情報処理装置１は、工業プロセスにおいて、複数の製造工程に亘るＭＳの変化を計算できるように構成されている。

　（ＰＣについて）
　以下、情報処理装置１において用いられる各データについて述べる。実施形態１では、工業プロセスにおけるＡｌの製造工程が、第１工程から第Ｎ工程までのＮ個の工程を含む場合を考える。Ｎは、Ｎ≧２を満たす自然数である。また、上記製造工程におけるｉ番目の工程を、第ｉ工程とする。ｉは、１≦ｉ≦Ｎを満たす自然数である。

　実施形態１では、第ｉ工程における製造条件（ＰＣ）を示すデータセットを、ＰＣ（ｉ）として表す。ＰＣ（ｉ）は、製造条件情報と称されてもよい。ＰＣ（ｉ）（つまり、ＰＣ（１）～ＰＣ（Ｎ））は、ユーザ入力によって、情報処理装置１にあらかじめ設定されている。

　ＰＣ（ｉ）に含まれる情報（データ、設定値）の例は、以下の１～５の通りである。但し、ＰＣ（ｉ）に含まれる情報は、以下の１～５の例に限定されない。上述のように、ＰＣ（ｉ）は、ＴＭＰ（より具体的には、以下に述べるＴＭＰ（ｉ，ｔ））を算出できるように設定されていればよい。

　＜１．熱間粗圧延工程においてＰＣ（ｉ）に含まれる情報の例＞
　被加工物の材質、スラブ寸法（幅、厚さ、長さ）、開始温度、終了温度、パス数、各パスの材料寸法（幅、厚さ、長さ）、各パスの圧下量、各パスの圧延時間、各パス間の保持時間、各パスの圧延速度、ワークロール径、およびクーラント量等。

　＜２．熱間仕上圧延工程においてＰＣ（ｉ）に含まれる情報の例＞
　被加工物の材質、パス数、開始温度、終了温度、パス数、各パスの材料寸法（幅、厚さ、長さ）、各パスの圧下量、各パスの圧延時間、各パス間の保持時間、各パスの圧延速度、各パスのワークロール径、巻取りコイルサイズ、巻取り後の冷却速度、およびスプール径等。

　＜３．冷間圧延工程においてＰＣ（ｉ）に含まれる情報の例＞
　被加工物の材質、パス数、開始温度、終了温度、パス数、各パスの材料寸法（幅、厚さ、長さ）、各パスの圧下量、各パスの圧延時間、各パス間の保持時間、各パスの圧延速度、各パスのワークロール径、巻取りコイルサイズ、巻取り後の冷却速度、およびスプール径等。

　＜４．熱間押出工程においてＰＣ（ｉ）に含まれる情報の例＞
　被加工物の材質、コンテナ形状、製品形状、開始コンテナ温度、開始ダイス温度、開始ビレット温度、ラム速度、および押出速度等。

　＜５．均質化処理工程においてＰＣ（ｉ）に含まれる情報の例＞
　鋳塊寸法（幅、厚さ、長さ）、昇温速度、保持温度、保持時間、冷却速度、および冷却時間等。

　（ＴＭＰについて）
　第ｉ工程において、当該工程が開始されてから経過した時間をｔとして表す。ｔ＝０は、第ｉ工程の開始時刻を表す。なお、各工程の区別のため、第ｉ工程における時刻ｔは、ｔｉと表記されてもよい。例えば、ｔ１は、第１工程における時刻ｔである。

　ＴＭＰ（ｉ，ｔ）は、第ｉ工程の時刻ｔにおいて、Ａｌに与えられる温度（加工温度）Ｔ、歪量（歪み）ε、および歪速度ｄεのそれぞれを示すデータセットである。ＴＭＰ（ｉ，ｔ）は、加工熱処理条件情報と称されてもよい。

　以下、ＴＭＰ（ｉ，ｔ）に含まれる温度Ｔ、歪量ε、および歪速度ｄεをそれぞれ、温度Ｔ（ｉ，ｔ）、歪量ε（ｉ，ｔ）、および歪速度ｄε（ｉ，ｔ）とも表す。ＴＭＰ（ｉ，ｔ）は、第ｉ工程の時刻ｔにおいて、Ａｌに与えられる熱、および当該Ａｌの加工状態（加工条件）を示す指標となる。ＴＭＰ（ｉ，ｔ）は、時刻ｔの関数、または、時刻ｔに対応するデータセットとして規定されてよい。

　以下に述べるように、温度Ｔ、歪量ε、および歪速度ｄεはいずれも、ＡｌにおけるＭＳの変化を高精度に予測（算出）するために必須のパラメータである。このため、ＴＭＰには、温度Ｔ、歪量ε、および歪速度ｄεの３種類のデータが含まれている必要がある。

　温度Ｔは、析出、回復、および再結晶等の現象（製造工程において生じ得る冶金的現象）に大きく影響を与えるパラメータである。

　歪量εは、材料に塑性加工が加えられる工程において、当該材料内に形成される加工組織を決定する支配的なパラメータである。

　また、歪速度ｄεは、熱間加工時に重要なパラメータである。Ａｌは、動的な回復または再結晶が特に生じやすい金属であるためである。このため、ＡｌにおけるＭＳの変化を高精度に算出するためには、温度Ｔおよび歪量εのみなならず、歪速度ｄεをも考慮することが必要である。

　上述のように、ユーザ入力によって設定されたＰＣ（ｉ）は、ＴＭＰ（ｉ，ｔ）を算出できるように設定されている。このため、情報処理装置１は、ＰＣ（ｉ）に基づいて、ＴＭＰ（ｉ、ｔ）を算出する。

　ＴＭＰを算出することが可能な製造工程であれば、任意の製造工程に対して実施形態１の方法を適用できる。このため、情報処理装置１は、任意の製造工程に対する机上検討が可能であるように構成できる。ＰＣ（ｉ）に基づいてＴＭＰ（ｉ，ｔ）を算出する方法は、製造工程に応じて任意に設定されてよい。以下の１～３に、その一例を示す。

　＜１．均質化処理工程におけるＴＭＰ（ｉ，ｔ）の算出方法＞
　均質化処理工程は、熱処理工程の１つである。このため、第ｉ工程が均質化処理工程である場合には、Ａｌに加工のための外力は印加されない。従って、歪量εおよび歪速度ｄεは、全ての時刻ｔにおいて０であると見なすことができる。つまり、ε（ｉ，ｔ）＝０およびｄε（ｉ，ｔ）＝０とすることができる。

　それゆえ、均質化処理工程では、温度Ｔの時間変化のみを考慮すればよい。つまり、温度Ｔ（ｉ，ｔ）＝Ｔ（ｔ）のみを算出すればよい。一例として、均質化処理工程における、開始温度Ｔｓ（℃）、昇温速度Ｈｕ（℃／ｓ）、保持温度Ｔｈ（℃）、保持時間ｔｈ（ｓ）、冷却速度Ｈｄ（℃／ｓ）、および終了温度Ｔｅ（℃）のそれぞれのデータが、ＰＣ（ｉ）に含まれる場合を考える（後述の図６も参照）。

　この場合、均質化処理工程における昇温→保温（保持）→冷却の各ステップにおける温度Ｔ（ｔ）を、以下の式（１）、

によって近似できる。

　式（１）では、均質化処理工程における温度Ｔ（ｔ）が、３つの線形関数（一次関数）によって近似されている。式（１）において、１段目の式は昇温時における温度Ｔ（ｔ）を、２段目の式は保温時における温度Ｔ（ｔ）を、３段目の式は冷却時における温度Ｔ（ｔ）を、それぞれ表す。式（１）に所定の時刻ｔを代入することにより、当該時刻ｔにおける温度Ｔ（ｔ）を得ることができる。

　なお、均質化処理工程における温度Ｔ（ｔ）を近似する式は、式（１）に限定されなくともよい。例えば、所定の時間ごとにＡｌの温度が直接的に測定されている場合には、各測定時刻における温度の測定値が、情報処理装置１（より具体的には、後述する第ｉ工程計算部）に供給されてよい。

　この場合、情報処理装置１は、各測定点のデータに対して公知の補間処理を用いることにより、所望の形式の関数（例：線形関数または非線形関数）によって温度Ｔ（ｔ）を表現してよい。

　＜２．冷間圧延工程におけるＴＭＰ（ｉ，ｔ）の算出方法＞
　一例として、冷間圧延工程における、ワークロール径、板厚、および圧延速度のデータが、ＰＣ（ｉ）に含まれる場合を考える。この場合、圧延時の歪量εおよび歪速度ｄεをそれぞれ、圧延理論に基づいて（圧延理論に基づく理論式を用いて）算出できる。

　例えば、入側の板厚をｈｉｎ（ｍ）、出側の板厚ｈｏｕｔ（ｍ）、ロールを半径Ｒｒ（ｍ）、圧延速度をＶｒ（ｍ／ｓ）とすると、以下の式（２）、

によって、歪量εおよび歪速度ｄεを算出できる。

　材料がロールを通過する時刻ｔについてのデータが、ＰＣ（ｉ）にさらに含まれていれば、式（２）の歪量εおよび歪速度ｄεをそれぞれ、ε（ｉ，ｔ）およびｄε（ｉ，ｔ）として表現できる。つまり、歪量εおよび歪速度ｄεをそれぞれ、時刻ｔの関数として表現できる。

　なお、冷間圧延工程において、単純な計算を行う場合には、温度Ｔとして室温を用いればよい。但し、一部の製品の場合には、加工による発熱によって生じるＭＳの変化を考慮することが好ましい。この場合、温度Ｔも、時刻ｔの関数Ｔ（ｔ）として表現することが好ましい。

　但し、理論式を用いて、冷間圧延工程における温度Ｔ（ｔ）を高精度に算出することは必ずしも容易ではない。そこで、ＰＣ（ｉ）には、冷間圧延工程の前後の温度（冷間圧延工程の開始時の温度、冷間圧延工程の終了時の温度）が含まれていることが好ましい。この場合、冷間圧延工程前後の温度を用いて、冷間圧延工程における温度Ｔの時間的な変化を、近似式を用いて概略的に表現できる。

　＜３．熱間押出工程におけるＴＭＰ（ｉ，ｔ）の算出方法＞
　熱間押出工程においては、材料の幾何学的形状が複雑であるため、単純な計算式のみからＴＭＰを算出することは必ずしも容易ではない。但し、コンテナ径、押出速度、および製品形状等のデータがＰＣ（ｉ）に含まれていれば、ＦＥＭを用いてＴＭＰを算出できる。

　（ＭＳについて）
　ＭＳ（ｉ，ｔ）は、第ｉ工程の時刻ｔにおけるＭＳを示すデータセットである。ＭＳ（ｉ，ｔ）は、材料組織情報と称されてもよい。一例として、ＭＳ（ｉ，ｔ）には、主要な添加元素の含有量、主要な添加元素の固溶量、代表的な第二相粒子の平均半径および数密度、再結晶率、結晶粒径、および歪エネルギーのデータが含まれていてよい。

　ここで、歪エネルギーとは、材料内の亜結晶粒が有する自由エネルギーの総量を意味する。亜結晶粒は、加工により生ずる転位によって、または、熱による転位の再配列によって生じる。

　ＭＳ（ｉ，ｔ）には、代表的な第二相粒子の粒子サイズの分布を示すデータがさらに含まれていることが好ましい。第二相粒子がＭＳの変化に及ぼす影響は、当該第二相粒子の粒子サイズによって異なるためである。それゆえ、第二相粒子の平均半径および数密度に加えて、当該第二相粒子の粒子サイズの分布をさらに用いることにより、ＭＳを高精度に算出できる。

　また、材料（Ａｌ）の再結晶部と未再結晶部とで、ＭＳ（ｉ，ｔ）が個別に取り扱われることがさらに好ましい。つまり、ＭＳ（ｉ，ｔ）を、材料の再結晶状態と未再結晶状態とで、別の情報（別のデータセット）として取り扱うことが好ましい。

　複数の製造工程を連結する計算においては、ある製造工程（第ｉ工程）における計算の完了時に、材料の再結晶部と未再結晶部が混在した状態のままで、計算結果が次の製造工程（第ｉ＋１工程）における計算を行う工程計算部（第ｉ＋１工程計算部）に与えられる場合が多いためである。

　この場合、ＭＳ（ｉ，ｔ）を再結晶部と未再結晶部との平均値（平均情報）として取り扱うよりも、ＭＳ（ｉ，ｔ）を再結晶部と未再結晶部とで個別の情報として取り扱う方が、計算精度を向上させることができる。

　一例として、材料の再結晶の有無によって、当該材料における析出速度および加工組織の形成具合は大きく変化する。すなわち、材料の再結晶部と未再結晶部とで、複数の冶金的現象のそれぞれの時間変化の態様は大きく異なりうる。それゆえ、ＭＳ（ｉ，ｔ）を再結晶部と未再結晶部とで個別に取り扱うことにより、材料組織の予測をより高精度に行うことができる。

　また、ＭＳ（ｉ，ｔ）は、全ての製造工程（第１工程～第Ｎ工程）において、共通の項目（換言すれば、同じ種類のフォーマットまたはデータ構造）によって定義されている。これにより、情報処理装置１（より具体的には、工程間情報統合部１１）において、第ｉ工程計算部（後述）から取得したＭＳ（ｉ，ｔｆｉ）を、第ｉ＋１工程計算部（後述）へと容易に伝達できる。例えば、以下に述べる第１工程計算部１２ａから第２工程計算部１２ｂへと、ＭＳ（１，ｔｆ１）＝ＭＳ（２，０）を容易に伝達できる。

　（情報処理装置１の構成）
　図１は、情報処理装置１の要部の構成を示す機能ブロック図である。情報処理装置１は、制御部１０、入力部７０、表示部７５、および記憶部９０を備える。

　制御部１０は、情報処理装置１の各部を統括的に制御する。制御部１０の機能は、記憶部９０に記憶されたプログラムを、ＣＰＵ（Central Processing Unit）が実行することで実現されてよい。記憶部９０は、制御部１０が実行する各種のプログラム、および当該プログラムによって使用されるデータを格納する。

　制御部１０は、工程間情報統合部１１、工程計算部１２、および材料組織計算部１３を備える。工程計算部１２は、第１工程計算部１２ａから第Ｎ工程計算部１２Ｎまでの、Ｎ個の下位的な工程計算部（サブ工程計算部）を備える。

　工程計算部１２において、第ｉ工程における計算を行うサブ工程計算部を、第ｉ工程計算部と称する。図１には、便宜上、第１工程計算部１２ａ、第２工程計算部１２ｂ、および第Ｎ工程計算部１２Ｎが図示されている。なお、後述の実施例では、Ｎ＝３の場合を例示する。

　材料組織計算部１３は、組織計算モジュール群１３０を備える。組織計算モジュール群１３０は、複数の冶金的現象のそれぞれの時間変化を計算するモジュール（計算モジュール）である。

　組織計算モジュール群１３０は、例えば、計算モジュールとして、再結晶モジュール１３０ａ、固溶析出モジュール１３０ｂ、加工組織モジュール１３０ｃ、および回復モジュール１３０ｄを備える。

　再結晶モジュール１３０ａは、材料（Ａｌ）の再結晶の時間変化を算出する。固溶析出モジュール１３０ｂは、材料の固溶析出の時間変化を算出する。加工組織モジュール１３０ｃは、材料の加工組織の時間変化を算出する。回復モジュール１３０ｄは、材料の回復の時間変化を算出する。

　Ａｌは、各種の金属材料のうち、回復が生じやすい材料である。このため、回復モジュール１３０ｄを設けることにより、Ａｌの回復（Ａｌにおいて生じやすい冶金的現象）を考慮して、複数の相互作用的な冶金的現象の時間変化を算出できる。それゆえ、材料組織をより正確に予測することが可能となる。

　但し、組織計算モジュール群１３０に、さらなるモジュールが設けられてもよい。組織計算モジュール群１３０によって、複数の冶金的現象によって生じる、ＭＳの変化を計算できる。

　複数のモジュールのそれぞれ（例：再結晶モジュール１３０ａ、固溶析出モジュール１３０ｂ、加工組織モジュール１３０ｃ、および回復モジュール１３０ｄ）には、共通のデータ構造を有するＭＳ（ｉ，ｔ）が供給される。例えば、複数のモジュールのそれぞれには、図７に示される項目を有するＭＳ（ｉ，ｔ）が供給される。

　以下に述べるように、工程間情報統合部１１、工程計算部１２（第１工程計算部１２ａ～第Ｎ工程計算部１２Ｎ）、および材料組織計算部１３（組織計算モジュール群１３０）が設けられることにより、Ａｌの材料組織の変化を予測できる（以下に述べる図２も参照）。

　入力部７０は、ユーザの入力（ユーザ入力）を受け付ける。上述のように、ＰＣ（１）～ＰＣ（Ｎ）は、ユーザ入力によって設定される。また、ＭＳ（１，０）（つまり、第１工程における材料組織情報の初期値）も、ユーザ入力によって設定される（後述の図３も参照）。

　表示部７５は、各種のデータ（特に、制御部１０の計算結果）をグラフとして表示する。また、表示部７５は、ユーザ入力のための入力画面を表示してもよい。表示部７５が設けられることにより、ユーザにＧＵＩ（Graphic User Interface）を提供できるので、ユーザの利便が向上する。なお、入力部７０および表示部７５としてタッチパネルを用いた場合、入力部７０と表示部７５とを一体の部材として実現できる。

　（工程間情報統合部１１）
　図２は、情報処理装置１における全体的な処理の流れを例示する図である。図２には、工程間情報統合部１１と工程計算部１２と材料組織計算部１３との間のデータの流れ、および、処理の流れが概略的に示されている。図３は、工程間情報統合部１１の処理Ｓ１～Ｓ６の流れを例示するフローチャートである。以下、図２を参照しつつ、図３に基づいて、工程間情報統合部１１の動作について述べる。

　なお、図２では、各工程の区別のために、第ｉ工程におけるＴＭＰ（ｉ，ｔ）およびＭＳ（ｉ，ｔ）がそれぞれ、ＴＭＰ（ｉ，ｔｉ）およびＭＳ（ｉ，ｔｉ）として表記されている。

　まず、工程間情報統合部１１は、ユーザによって設定されたＭＳ（１，０）およびＰＣ（１）～ＰＣ（Ｎ）を取得する（Ｓ１）。そして、工程間情報統合部１１は、ｉ＝１に設定する（Ｓ２）。ｉは、図３の処理におけるループカウントである。

　続いて、工程間情報統合部１１は、第ｉ工程計算部に、ＭＳ（ｉ，０）およびＰＣ（ｉ）を供給する（Ｓ３）。第ｉ工程計算部がＭＳ（ｉ，ｔｆｉ）を取得する処理については、後述の図４に示されている。ｔｆｉは、第ｉ工程の終了時刻である。ｔｆｉは、あらかじめ設定されていてもよいし、あるいは第ｉ工程計算部によって変更可能に設定されてもよい。

　一例として、ｉ＝１の場合には、工程間情報統合部１１は、第１工程計算部１２ａに、ＭＳ（１，０）およびＰＣ（１）を供給する（図２における、工程間情報統合部から第１工程計算部に向かう矢印を参照）。この場合、第１工程計算部１２ａは、ＭＳ（１，０）およびＰＣ（１）を取得して、ＭＳ（ｉ，ｔｆ１）を算出する。

　続いて、工程間情報統合部１１は、第ｉ工程計算部からＭＳ（ｉ，ｔｆｉ）を取得する（Ｓ４）。一例として、ｉ＝１の場合には、工程間情報統合部１１は、第１工程計算部１２ａから、ＭＳ（１，ｔｆ１）を取得する（図２における、第１工程計算部から工程間情報統合部に向かう矢印を参照）。

　続いて、工程間情報統合部１１は、ＭＳ（ｉ，ｔｆｉ）＝ＭＳ（ｉ＋１，０）として、ＭＳ（ｉ＋１，０）を設定する（Ｓ５）。つまり、工程間情報統合部１１は、第ｉ工程の完了時における材料組織情報を、第ｉ＋１工程（次工程）の開始時における材料組織情報として設定する。

　そして、工程間情報統合部１１は、ｉ＝Ｎであるか否かを判定する（Ｓ６）。つまり、工程間情報統合部１１は、現時点の計算対象である第ｉ工程が、最終工程（第Ｎ工程）であるか否かを判定する。

　ｉ＝Ｎでない場合（Ｓ６でＮＯ）、工程間情報統合部１１は、ｉ＝ｉ＋１に設定する（Ｓ７）。つまり、工程間情報統合部１１は、ｉの値を１つカウントアップする。そして、Ｓ３に戻り、同様の処理を繰り返す。

　一例として、ｉ＝１の場合を考える。この場合、Ｓ５において、工程間情報統合部１１は、ＭＳ（１，ｔｆ１）＝ＭＳ（２，０）として、ＭＳ（２，０）を設定する。続いて、工程間情報統合部１１は、ｉ＝２としてｉを設定し、Ｓ３に戻る。

　そして、Ｓ３において、工程間情報統合部１１は、第２工程計算部１２ｂに、ＭＳ（２，０）（＝ＭＳ（１，ｔｆ１））およびＰＣ（２）を供給する（図２における、工程間情報統合部から第２工程計算部に向かう矢印を参照）。

　続いて、工程間情報統合部１１は、Ｓ４において、第２工程計算部１２ｂからＭＳ（２，ｔｆ２）を取得する（図２における、第２工程計算部から工程間情報統合部に向かう矢印を参照）。以下、第３工程以降についても同様である。

　そして、ｉの値が順次カウントアップされ、ｉ＝Ｎとなった場合（Ｓ６でＹＥＳ）に、工程間情報統合部１１は、処理を終了する。工程間情報統合部１１は、第１工程から第Ｎ工程までの計算において得られた全ての材料組織情報を、最終的な計算結果として出力してよい。

　また、工程間情報統合部１１は、第１工程から第Ｎ工程までの計算において得られた全ての加工熱処理条件情報を、最終的な計算結果として出力してよい。以下に述べるように、加工熱処理条件情報（ＴＭＰ（ｉ，ｔ））は、第ｉ工程計算部によって算出される。

　工程間情報統合部１１は、ユーザ入力に応じて、最終的な計算結果の一部（例：ユーザが閲覧を所望する計算結果）のみを、表示部７５に選択的に表示させてよい。また、工程間情報統合部１１は、上述の各処理の途中において、記憶部９０に計算結果を格納してもよい。

　なお、記憶部９０には、ユーザによってあらかじめ測定された、実際の材料におけるＭＳを示すデータ（以下、ＭＳ実測データ）（実測情報）が記憶されていてよい。工程間情報統合部１１は、ＭＳ実測データを用いて、当該ＭＳ実測データの測定対象となった材料と同じ製造条件（ＰＣ（ｉ））のもとで算出されたＭＳを評価してよい。

　より具体的には、工程間情報統合部１１は、ＭＳ実測データに対するＭＳ（ｉ，ｔ）の誤差を算出することにより、当該ＭＳ（ｉ，ｔ）の誤差（以下、計算誤差）を評価してよい。さらに、工程間情報統合部１１は、計算誤差が小さくなるよう、工程計算部および材料組織計算部における計算に用いられるパラメータ（例：係数）を調整してよい。このように、工程間情報統合部１１は、ＭＳ実測データを用いて、計算精度を向上させる機能を有していることが好ましい。

　なお、ＭＳ実測データには、ＭＳの全項目（後述の図７を参照）に対応する実測値が含まれている必要はない。ＭＳ実測データには、ＭＳの項目の少なくとも一部分に対応する実測値が含まれていればよい。

　（工程計算部１２）
　図４は、工程計算部１２に含まれる第ｉ工程計算部の処理Ｓ１１～Ｓ１８の流れを例示するフローチャートである。以下、図４を参照し、第ｉ工程計算部の動作について述べる。

　なお、図４におけるｉの値は、上述の処理Ｓ２において、工程間情報統合部１１によってあらかじめ設定されているものとする。この点については、後述する図５についても同様である。つまり、図４および図５はいずれも、所定の第ｉ工程における処理を示す。

　まず、第ｉ工程計算部は、工程間情報統合部１１から、ＭＳ（ｉ，０）および製造条件ＰＣ（ｉ）を取得する（Ｓ１１）。なお、上述のように、ＭＳ（ｉ，０）は、ＭＳ（ｉ－１，ｔｆｉ）に等しい。

　続いて、第ｉ工程計算部は、第ｉ工程における時刻ｔ（時刻ｔｉ）を、ｔ＝０（開始時刻，初期時刻）に設定する（Ｓ１２）。そして、第ｉ工程計算部は、ＰＣ（ｉ）に基づいてＴＭＰ（ｉ，ｔ）を算出する（Ｓ１３）。一例として、ｉ＝１かつｔ＝ｔ１＝０の場合を考える。この場合、第１工程計算部１２ａは、ＭＳ（１，０）を材料組織計算部１３に供給する（図２も参照）。

　なお、上述のとおり、第ｉ工程計算部においてＴＭＰ（ｉ，ｔ）を算出する方法は、製造工程に応じて異なる手法が用いられてよい。例えば、焼鈍等の工程であれば、炉の温度設定プログラムから、温度履歴を単純に模擬した数式を設定できる。また、圧延または引抜き等の、比較的単純な塑性加工を伴う工程であれば、圧延理論または材料力学の理論に基づいて数式を設定できる。この場合、第ｉ工程計算部は、当該数式を用いてＴＭＰ（ｉ，ｔ）を算出できる。

　あるいは、鍛造または押出等の複雑な工程については、ＦＥＭなどのシミュレーション技術を用いることにより、第ｉ工程計算部においてＴＭＰ（ｉ，ｔ）を算出すればよい。

　そして、第ｉ工程計算部は、ＭＳ（ｉ，ｔ）およびＴＭＰ（ｉ，ｔ）を材料組織計算部１３に供給する（Ｓ１４）。一例として、ｉ＝１かつｔ＝ｔ１＝０の場合には、第１工程計算部１２ａは、ＭＳ（１，０）およびＴＭＰ（１，０）を材料組織計算部１３に供給する（図２も参照）。

　続いて、第ｉ工程計算部は、材料組織計算部１３から、ＭＳ（ｉ，ｔ＋Δｔ）を取得する（Ｓ１５）。一例として、ｉ＝１かつｔ＝ｔ１＝０の場合、第１工程計算部１２ａは、材料組織計算部１３から、ＭＳ（１，Δｔ）を取得する（図２も参照）。なお、Δｔは微小時間である。Δｔの値は、材料組織計算部１３において設定されてよい。

　なお、Δｔの値は、一定であってもよいし、一定でなくともよい。一例として、異なるｉに対して、異なるΔｔの値が設定されてよい。あるいは、同じｉに対しても、計算の途中において、Δｔの値が変更されてもよい。

　続いて、第ｉ工程計算部は、ｔ＝ｔｆｉであるか否かを判定する（Ｓ１６）。つまり、第ｉ工程計算部は、時刻ｔが第ｉ工程における終了時刻まで進行したか否かを判定する。

　ｔ＝ｔｆｉでない場合（Ｓ１６でＮＯ）、第ｉ工程計算部は、ｔ＝ｔ＋Δｔに設定する（Ｓ１７）。つまり、第ｉ工程計算部は、第ｉ工程における時刻ｔを、Δｔだけカウントアップ（進行）させる。そして、Ｓ１３に戻り、同様の処理を繰り返す。

　他方、第ｉ工程計算部は、ｔ＝ｔｆｉである場合（Ｓ１６でＹＥＳ）、第ｉ工程計算部は、工程間情報統合部１１に、ＭＳ（ｉ，ｔｆｉ）を供給する（Ｓ１８）（上述の図３のＳ４・Ｓ５も参照）。一例として、ｉ＝１かつｔ＝ｔ１＝０の場合、第１工程計算部１２ａは、工程間情報統合部１１に、ＭＳ（１，ｔｆ１）を供給する。

　以上のように、第ｉ工程計算部は、ＰＣ（ｉ）に基づいてＴＭＰ（ｉ，ｔ）を算出する。そして、第ｉ工程計算部は、材料組織計算部１３にＭＳ（ｉ，ｔ＋Δｔ）を順次算出させてゆき、材料組織計算部１３からＭＳ（ｉ，ｔｆｉ）を取得する。

　そして、上述のように、工程間情報統合部１１は、第ｉ工程計算部からＭＳ（ｉ，ｔｆｉ）を取得し、当該ＭＳ（ｉ，ｔｆｉ）をＭＳ（ｉ＋１，０）として第ｉ＋１工程計算部に供給する。このように、工程間情報統合部１１は、第ｉ工程計算部における第ｉ工程の計算が完了すると、第ｉ＋１工程計算部に第ｉ＋１工程（次の製造工程）の計算を開始させる。

　（材料組織計算部１３）
　図５は、材料組織計算部１３の処理Ｓ２１～Ｓ２３の流れを例示するフローチャートである。以下、図５を参照し、材料組織計算部１３の動作について述べる。

　まず、材料組織計算部１３は、第ｉ工程計算部からＭＳ（ｉ，ｔ）およびＴＭＰ（ｉ，ｔ）を取得する（Ｓ２１）。続いて、材料組織計算部１３は、組織計算モジュール群１３０によって、ＭＳ（ｉ，ｔ）およびＴＭＰ（ｉ，ｔ）に基づいて、ＭＳ（ｉ，ｔ＋Δｔ）を算出する（Ｓ２２）。

　そして、材料組織計算部１３は、算出したＭＳ（ｉ，ｔ＋Δｔ）を、第ｉ工程計算部に供給する（Ｓ２３）。

　組織計算モジュール群１３０は、現時点における加工熱処理条件（ＴＭＰ（ｉ，ｔ））を用いて、現時点の材料組織（ＭＳ（ｉ，ｔ））から、微小時間Δｔ後の材料組織（ＭＳ（ｉ，ｔ＋Δｔ））を算出（予測）する機能を有する。

　組織予測モジュール群に含まれる各モジュールは、公知の方法（例：ＭＳの予測のための公知の理論式または実験式）を用いた計算手法によって、各冶金的現象に伴うＭＳの時間的な変化を計算（予測）してよい。つまり、各モジュールは、所定のモデル（計算モデル、物理モデル）を用いて、ＭＳの変化を予測してよい。但し、各モジュールにおけるＭＳの予測方法は特に限定されない。組織予測モジュール群は、ＭＳの変化を予測できれば、任意の計算手法を用いてよい。

　一例として、再結晶モジュール１３０ａは「Acta Materialia 51 (2003) 1453-1468」に記載の方法を、固溶析出モジュール１３０ｂは、「Acta Materialia 45 (1997) 4231-4240」に記載の方法を、加工組織モジュール１３０ｃは、「Journal of Materials Processing Technology 117 (2001) 333-340」に記載の方法を、それぞれ用いて実装されてよい。

　各モジュールは、各モジュールの計算において必要な材料組織パラメータを相互に参照する。例えば、再結晶モジュール１３０ａが再結晶を計算する場合には、第二相粒子の情報が不可欠となる。このため、再結晶モジュール１３０ａは、固溶析出モジュール１３０ｂによって算出された固溶析出状態の計算結果を参照し、当該計算結果を第二相粒子の情報として用いる。

　固溶析出モジュール１３０ｂでの固溶析出計算において、粒界析出を考慮するためには、再結晶モジュール１３０ａにおいて算出された、結晶粒組織の計算結果が必要である。このため、固溶析出モジュール１３０ｂは、再結晶モジュール１３０ａの計算結果を参照し、固溶析出を計算する。

　各モジュールにおける計算を、微小時間Δｔごとに行うことで、複数の冶金的現象の相互作用を考慮して、ＭＳの時間変化を計算することが可能となる。つまり、複数の相互作用的な冶金的現象の時間変化を算出して、ＭＳ（ｉ，ｔ＋Δｔ）を得ることができる。なお、Δｔが大きい場合、材料組織の時間変化を精度よく計算することができない。このため、Δｔはある程度小さい値として設定（選択）される必要がある。

　ここで、組織計算モジュール群１３０での計算における、所定の材料組織パラメータ（材料組織の状態を示す変数）をＺ（ｔ）と表す。Δｔは、時刻ｔから時刻ｔ＋Δｔまでの間におけるＺ（ｔ）の変化量に着目して設定されてよい。

　一例として、Δｔは、Ｚ（ｔ）からＺ（ｔ＋Δｔ）への変化量が、０．５％以下となるように設定されてよい。つまり、Δｔは、Ｚ（ｔ）からＺ（ｔ＋Δｔ）への時間的な変化が、ある程度小さくなるように設定されればよい。

　なお、Δｔは、所定の１つのＺ（単一の材料組織パラメータ）または単一のモジュールによって設定されない。Δｔは、組織計算モジュール群１３０（全てのモジュール）において用いられる全てのＺの中で、時間領域における変化量が最も大きいＺにより決定される。このように、変化量が最も大きいＺに対してΔｔを設定することにより、複数の現象の相互作用を十分に考慮できる。

　（情報処理装置１の効果）
　情報処理装置１によれば、種々のＰＣ（ｉ）（製造条件）に対して、ＭＳ（ｉ，ｔ）を算出できる。つまり、各工程の進行に伴うＭＳの変化を予測できる。それゆえ、製品（Ａｌ）の製造時に得られるＭＳ（すなわち、ＭＳ（Ｎ，ｔｆＮ））を従来よりも高精度に予測できる。すなわち、工業プロセスにおけるＡｌの材料組織を従来よりも高精度に予測できる。

　このため、例えば、ＰＣ（ｉ）に含まれるパラメータ（製造パラメータ）を所定の条件範囲中で変化させて計算を行うことで、当該製造パラメータがＭＳに及ぼす影響を系統的に知ることができる。

　また、複数の製造パラメータのそれぞれに条件範囲を指定し、当該範囲内で各製造パラメータをランダムに変化させて繰り返し計算させることで、所望のＭＳを得るためのＰＣ（ｉ）を探索することもできる。

　このように、情報処理装置１によれば、ＭＳに起因する特性を安定化させるように、ＰＣ（ｉ）を設定することが可能となる。それゆえ、工業プロセスにおける試作検討のコストの回数を大幅に減少させることができる。

　さらに、情報処理装置１によれば、以下の実施例において詳述するように、共通のデータ構造を有するＭＳ（ｉ，ｔ）を用いて、材料組織計算部１３に計算を行わせることができる。以下、当該構成の利点について、実施例を参照し説明する。

　（実施例）
　「鋳造工程→均質化処理工程→熱間粗圧延工程→熱間仕上圧延工程→冷間圧延工程」という製造工程からなるＨ１ｎ調質のＡｌ板材（製品）を考える。冷間圧延後の製品強度の安定化を目的として、情報処理装置１による計算を利用した例を示す。

　Ｈ１ｎ調質の製品の場合、冷間圧延後の製品強度は、合金成分および冷間圧延の加工度に強く依存する。但し、冷間圧延後の製品強度は、熱間仕上圧延後の再結晶率にも影響されることが知られている。

　合金成分および冷間圧延の加工度は、ＰＣとして容易に制御できる。他方、熱間仕上圧延後の再結晶率は、安定した制御が難しく、製品強度のバラツキが生じる一因となっている。そこで、発明者らは、情報処理装置１を用いて、「均質化処理工程→熱間粗圧延工程→熱間仕上圧延工程」の後のＭＳの予測を行い、製品強度を安定化させる対策を検討した。

　発明者らは、「均質化処理工程→熱間粗圧延工程→熱間仕上圧延工程」という連続した３つの工程を、計算対象とした。つまり、発明者らは、均質化処理工程を第１工程、熱間粗圧延を第２工程、熱間仕上圧延工程を第３工程として（Ｎ＝３として）、情報処理装置１による計算を行った。

　図６には、本実施例におけるＰＣの項目が示されている。具体的には、図６には、ＰＣ（１）（均質化処理工程におけるＰＣのデータセット）、ＰＣ（２）（熱間粗圧延工程におけるＰＣのデータセット）、およびＰＣ（３）（熱間仕上圧延工程におけるＰＣのデータセット）が、それぞれ示されている。

　なお、以下に述べるように、ＭＳ（１，０）としては、工場での製造指示の条件および製造実績から、均質化処理に供されるＡｌ鋳塊のＭＳが与えられている。

　これらのＰＣを用いて、各工程におけるＴＭＰを算出した。まず、均質化処理工程（第１工程）におけるＴＭＰについて、温度Ｔ（ｔ）＝Ｔ（１，ｔ）を上述の式（１）によって計算した。上述の通り、歪量ε（ｔ）＝ε（１，ｔ）および歪速度ｄε（ｔ）＝ｄε（１，ｔ）はいずれも０とした。

　熱間粗圧延工程（第２工程）および熱間仕上圧延工程（第３工程）では、圧延中の温度変化は不明である。このため、圧延による温度変化は一瞬で生じると仮定して、温度Ｔ（ｔ）（温度Ｔ（２，ｔ）およびＴ（３，ｔ））を算出した。

　また、歪量ε（ｔ）（歪量ε（２，ｔ）およびε（３，ｔ））および歪速度ｄε（ｔ）（歪速度ｄε（２，ｔ）およびε（３，ｔ））については、圧延ロールを通る瞬間のみ非零の数値が生じるとした。圧延ロールを通る瞬間以外の時刻では、歪量ε（ｔ）および歪速度ｄε（ｔ）はいずれも０に設定した。

　歪量εおよび歪速度ｄεの大きさは、圧下量、圧延速度、およびロール径を用いて、圧延理論により算出した。熱間仕上圧延工程では、材料がコイルに巻き取られた後の冷却中にもＭＳの変化が起こる。熱間仕上圧延工程におけるＴＭＰには、加工は含まれないので、単純に冷却速度の実測値から温度変化Ｔ（ｔ）を算出した。

　図７には、本実施例におけるＭＳの項目が示されている。図７のＭＳは、第１工程～第３工程に共通して用いられている。鋳塊をサンプリングして材料組織観察を行い、第１工程（均質化処理工程）におけるＭＳの初期値（すなわち、ＭＳ（１，０））を実測した。なお、第１工程における材料の初期状態は鋳塊であるため、再結晶率は１であり、歪エネルギーは０である。

　以上のデータに基づき、情報処理装置１を用いてＭＳの計算を行った。図８には、本実施例における全体的な処理の流れが例示されている。図８は、上述の図３において、Ｎ＝３とした場合に相当する。

　ＭＳを算出するモデル（物理モデル、計算モデル）を調整するための実測データを取得するために、１０ロット分の工場材を抽出した。当該工場材について、均質化処理工程後および熱間仕上圧延工程のそれぞれの後においてサンプリングを行い、ＭＳを実測した。続いて、（ｉ）均質化処理工程後の析出物の平均粒子半径Ｒｐと、（ｉｉ）熱間仕上圧延工程後の再結晶率Ｙおよび再結晶部分の平均粒径Ｄと、を定量化した。

　これらの実測データを用いて、モデルの調整を行った。ここで、モデル調整用の材料組織因子をＸとする。Ｘは、上述のＲｐ、Ｄ、およびＹのうちの任意の１つである。また、ＸｅｘｐをＸの実測値、ＸｃａｌをＸの計算値とする。

　本実施例では、情報処理装置１（工程間情報統合部１１）において、以下の式（３）、

によって、Ｘの誤差評価値Ｅを算出した。そして、Ｒｐ、Ｄ、およびＹの３つについての誤差評価値の平均値Ｅｍを最小化するよう、モデル（より具体的には、モデルに含まれるパラメータ）を調整した。

　実施形態１では、ＭＣＭＣ法（マルコフ連鎖モンテカルロ法，Markov chain Monte Carlo method）を用いた学習によって、モデルの調整を行った。但し、モデルの調整には、その他の任意の公知の手法が用いられてよい。例えば、学習を利用しない方法によって、モデルの調整を行ってもよい。

　図９には、上述のように調整したモデルを用いた計算結果が示されている。図９における各ロットに対する記号は、
　　Ｒｅｘｐ：Ｒｐの実測値；
　　Ｒｃａｌ：Ｒｐの計算値；
　　ＥＲｐ：Ｒｐの誤差評価値；
　　Ｄｅｘｐ：Ｄの実測値；
　　Ｄｃａｌ：Ｄの計算値；
　　ＥＤ：Ｄの誤差評価値；
　　Ｙｅｘｐ：Ｙの実測値；
　　Ｙｃａｌ：Ｙの計算値；
　　ＥＹ：Ｙの誤差評価値；
　　Ｅｍ：誤差評価値の平均値；
の通りである。なお、Ｅｍ＝（ＥＲｐ＋ＥＤ＋ＥＹ）／３である。また、１０ロット間におけるＥｍの平均値を、全体誤差平均値Ｅｍｍと称する。これらの記号は、後述の比較例および各変形例においても同様である。全体誤差平均値Ｅｍｍが０．２を下回れば、従来の経験則に依存した試作検討（トライアンドエラー）に代替しうる予測精度が得られることが期待される。

　図９に示されるように、本実施例では、各ロットの誤差評価値について、十分に小さい平均値Ｅｍが得られた。このため、全体誤差平均値Ｅｍｍについても、０．２を十分に下回る小さい値（０．０８）が得られた（図９のハッチング部を参照）。

　製造パラメータを様々に変更させてＭＳを算出することで、当該製造パラメータに応じて、仕上熱間圧延工程後の再結晶率の計算値Ｙｃａｌを得ることができる。それゆえ、情報処理装置１によれば、仕上熱間圧延工程後のＹｃａｌの値によって、冷間圧延後の製品強度をシミュレーションによって検討できる（机上検討できる）。

　図１０には、実施例、比較例（後述）、および各参考例（後述）のそれぞれについて、計算処理の違いが示されている。図１０に示されるように、本実施例では、各工程間の情報連携（データの受け渡し）は、工程間情報統合部１１によって行われる。また、各工程におけるＴＭＰの計算は、工程計算部１２によって行われる。また、各工程におけるＭＳの計算は、材料組織計算部１３によって行われる。

　従来、冷間圧延後の製品強度を検討する場合には、工場で試験材を１０ロット程度作成し、仕上熱延工程後の再結晶率が１００％に近くなる製造条件を、試作検討によって模索せざるを得なかった。例えば、Ａｌ圧延材の１０トン鋳塊１本分のコストをＣ（単位：円）とすると、Ｋ回の試作検討を行った場合には、概ね「Ｃ×Ｋ」の費用が発生する。

　しかしながら、上述のように、情報処理装置１によれば、従来よりも正確に、工業プロセスにおけるＡｌのＭＳを予測できる。このため、精度の高い机上検討により、試作検討の回数を、例えば半分以下に削減することが期待できる。従って、例えば「Ｃ×Ｋ／２」以上の費用低減が可能となる。

　図１１には、実施例、比較例、および各参考例のそれぞれについての比較結果が示されている。本実施例では、工程間情報統合部１１によって、工程計算部１２に対して、計算に必要な情報（ＰＣ（ｉ）およびＭＳ（ｉ，０））を自動的に供給できる。このため、全工程（第１工程～第Ｎ工程）の計算を、工程間情報統合部１１が工程計算部１２に一括して行わせることができる。それゆえ、ユーザの工数（計算のための工数）を低減できる。つまり、ユーザに対して利便性に優れた計算システム（計算プログラム）を提供できる。

　さらに、本実施例では、工程間情報統合部１１によって全工程に亘るモデル調整を行うことができるため、工程全体での計算の最適化を図ることができる。それゆえ、高い計算精度を得ることも可能となる。

　また、本実施例では、全工程において、共通の項目（データ構造）を有する材料組織情報（ＭＳ（ｉ，ｔ））を用いて、材料組織計算部１３に計算を行わせることが可能である。つまり、共通の（１つの）材料組織計算部１３において、ＭＳを計算することが可能である。これにより、工程に応じて個別の（異なる）材料組織計算モジュールを準備する必要がないため、計算システムの開発工数を大幅に低減できる。

　本実施例では、「均質化処理工程（第１工程）→熱間粗圧延処理工程（第２工程）→熱間仕上圧延処理工程（第３工程）」の３つの工程を検討対象とした。情報処理装置１の構成によれば、これらの工程とは異なる別の工程（第４工程）を検討する場合にも、当該別の工程に対応する工程計算部（第４工程計算部）を開発するのみでよい。

　このように、本発明の一態様に係る計算システムは、検討すべき工程に応じてカスタマイズ可能なパッケージ（ＭＳ予測パッケージ）であるので、他の工程に適用するためのプログラムの開発が簡便となる。それゆえ、様々な工程についての検討を行う場合にも、計算システムの開発に要する人工数を大幅に低減できる。

　（比較例）
　従来技術によるＭＳの予測計算の一例を、比較例として示す。比較例では、上述の実施例と同じＰＣおよびＴＭＰを用いて、実施例と同様の計算を行った。図１２は、比較例におけるＭＳの項目を示す図である。図１３は、比較例における全体的な処理の流れを例示する図である。

　図１３に示されるように、比較例における情報処理装置１Ｘは、第１工程（均質化処理工程）の材料組織計算モジュール（第１材料組織計算モジュール）と、第２工程（熱間粗圧延工程）の材料組織計算モジュール（第２材料組織計算モジュール）と、第３工程（熱間仕上圧延工程）の材料組織計算モジュール（第３材料組織計算モジュール）と、を備える。

　つまり、情報処理装置１Ｘは、３つの工程のそれぞれに対して、個別の材料組織計算モジュールを備える。この点において、比較例における情報処理装置は、情報処理装置１と異なる。各材料組織計算モジュール（第１材料組織計算モジュール～第３材料組織計算モジュール）には、各工程のＴＭＰ計算のアルゴリズムとＭＳ計算のアルゴリズムとが実装されている。各材料組織計算モジュールは、実施例と同様の方法によって、ＰＣからＴＭＰおよびＭＳを算出する。

　また、図１３に示されるように、情報処理装置１Ｘは、工程間情報統合部１１を備えていない。この点においても、情報処理装置１Ｘは、情報処理装置１と異なる。情報処理装置１Ｘでは、現在の第ｉ工程（例：第１工程）から次の第ｉ＋１工程（例：第２工程）へと計算を進める場合におけるデータの受け渡しは、ユーザがマニュアル操作によって行う必要がある（図１０および図１１も参照）。

　情報処理装置１Ｘでは、工程ごとに独立した材料組織計算モジュールによって計算が行われる。このため、各工程の計算において使用される材料組織情報の種類が異なる。具体的には、第１材料組織計算モジュール（均質化処理の材料組織計算モジュール）における計算では、再結晶率、結晶粒径、および歪エネルギーの値は使用されない。

　このため、図１２のハッチング部に示されるように、第１材料組織計算モジュールでは、再結晶率、結晶粒径、および歪エネルギーは、材料組織情報の項目に含まれていない。それゆえ、第１材料組織計算モジュールには、材料組織情報として、再結晶率、結晶粒径、および歪エネルギーの値を入出力する機能が設けられていない。

　従って、第２工程（熱間粗圧延工程）の計算のための材料組織情報（ＭＳ（２，０））を第２材料組織計算モジュールに受け渡す場合には、再結晶率、結晶粒径、および歪エネルギーの値をユーザが手動で選択する処理が必要となる。このため、情報処理装置１Ｘでは、ＭＳの計算に要するユーザの工数が、情報処理装置１に比べて増大する。

　さらに、情報処理装置１Ｘでは、各工程間の情報が個別に取り扱われているため、個々の個別の材料組織計算モジュール内でしかモデル調整を行うことができない。つまり、全工程に亘るモデル調整を行うことができない。それゆえ、比較例では、実施例に比べて計算精度が低下する。

　図１４には、比較例における計算結果が示されている。図１４に示されるように、比較例では、各ロットにおける誤差評価値の平均値Ｅｍは、実施例の場合に比べて十分に大きい。さらに、全体誤差平均値Ｅｍｍについても、０．２を十分に上回る大きい値（０．３５）となった（図１４のハッチング部を参照）。すなわち、比較例における計算では、従来の試作検討に代替しうる予測精度が得られないことが確認された。

　また、情報処理装置１Ｘでは、工程ごとに材料組織計算モジュール（プログラム）を設計する必要があるため、プログラム設計に要する人工数が、実施例と比較して大となる。例えば、情報処理装置１Ｘを別の工程（例：第４工程）の検討に適用する場合には、第４工程に応じた材料組織計算モジュールを設計する必要がある。このことから、情報処理装置１Ｘは、別の工程への検討に適用するためのカスタマイズ性が、情報処理装置１に比べて十分に低いと言える。

　（参考例１）
　参考例１では、上述の実施例と同じＰＣ、ＴＭＰ、およびＭＳを用いて、実施例と同様の計算を行った。図１５は、参考例１における全体的な処理の流れを例示する図である。図１５に示されるように、参考例における情報処理装置１Ｐは、比較例の情報処理装置１Ｘに対して、情報処理装置１の工程間情報統合部１１を付加した構成である。なお、参考例１におけるＭＳの項目は、比較例と同様である（図１２を参照）。

　上述の比較例と同様に、第１材料組織計算モジュール（均質化処理の材料組織計算モジュール）における計算では、再結晶率、結晶粒径、および歪エネルギーの値は使用されない。但し、情報処理装置１Ｐでは、工程間情報統合部１１によって、第１材料組織計算モジュールにおける再結晶率、結晶粒径、および歪エネルギーの値が、そのまま第２材料組織計算モジュールに受け渡される。

　参考例１においても、工程間情報統合部１１は、実施例と同様に、各工程のモデル調整を行う。但し、各工程の材料組織計算部は個別であるため（各工程において共有されていないため）、ある工程でのモデル調整の結果は、別の工程には反映されない。

　例えば、参考例１において、第１工程（均質化処理工程）では、実測された平均粒子半径（Ｒｅｘｐ）との誤差を低減するように、モデル調整が行われるが、当該モデル調整の結果は、第２工程（熱間粗圧延工程）および第３工程（熱間仕上圧延工程）には反映されない（図１０および図１１も参照）。

　参考例１によれば、工程間情報統合部１１によって、各工程間の各工程の情報連携が自動化されるため、ＭＳの計算に要するユーザの工数を、比較例に比べて低減することができる。しかしながら、参考例１においても、比較例と同様に、工程ごとに独立した材料組織計算モジュールによって計算が行われるため、全行程に亘るモデル調整ができない。それゆえ、参考例１においても、実施例と同程度の計算精度を得ることはできない。

　図１６には、参考例における計算結果が示されている。図１６に示されるように、参考例１においても、各ロットにおける誤差評価値の平均値Ｅｍは、実施例の場合に比べて十分に大きい。さらに、全体誤差平均値Ｅｍｍについても、０．２を十分に上回る大きい値（０．２９）となった（図１６のハッチング部を参照）。すなわち、参考例１における計算においても、従来の試作検討に代替しうる予測精度が得られないことが確認された。

　また、情報処理装置１Ｐにおいても、情報処理装置１Ｘと同様に、工程ごとに材料組織計算モジュールを設計する必要があるため、プログラム設計に要する人工数が、実施例と比較して大となる。つまり、情報処理装置１Ｐも、情報処理装置１Ｘと同様に、カスタマイズ性が情報処理装置１に比べて十分に低いと言える。

　（参考例２）
　参考例２では、上述の実施例と同じＰＣ、ＴＭＰ、およびＭＳを用いて、実施例と同様の計算を行った。図１７は、参考例１における全体的な処理の流れを例示する図である。図１７に示されるように、参考例における情報処理装置１Ｑは、情報処理装置１から工程間情報統合部１１を取り除いた構成である。なお、参考例２におけるＭＳの項目は、実施例と同様である（図７を参照）。

　情報処理装置１Ｑは、工程間情報統合部１１を備えていないので、比較例の情報処理装置１Ｘと同様に、工程間の情報を連携する機能を有していない。このため、第ｉ工程から第ｉ＋１工程へと計算を進める場合におけるデータの受け渡しは、ユーザがマニュアル操作によって行う必要がある（図１０および図１１も参照）。それゆえ、情報処理装置１Ｑにおいても、情報処理装置１Ｘと同様に、ＭＳの計算に要するユーザの工数が、情報処理装置１に比べて増大する。

　また、情報処理装置１Ｑは、各工程計算部に共通の材料組織計算部１３を有しているものの、工程間情報統合部１１を備えていないので、全工程に亘るモデル調整を行うことができない。それゆえ、参考例２においても、実施例と同程度の計算精度を得ることはできない。

　図１８には、情報処理装置１Ｑにおける計算結果が示されている。図１８に示されるように、参考例２においても、各ロットにおける誤差評価値の平均値Ｅｍは、実施例の場合に比べて十分に大きい。さらに、全体誤差平均値Ｅｍｍについても、０．２を十分に上回る大きい値（０．４０）となった（図１８のハッチング部を参照）。すなわち、参考例２における計算においても、従来の試作検討に代替しうる予測精度が得られないことが確認された。

　但し、情報処理装置１Ｑは、各工程計算部に共通の材料組織計算部１３を有しているため、情報処理装置１と同様に、高いカスタマイズ性を有している。

　〔変形例〕
　実施形態１では、ＴＭＰおよびＭＳを示す情報（加工熱処理条件情報および材料組織情報）がそれぞれ、ＴＭＰ（ｉ，ｔ）およびＭＳ（ｉ，ｔ）として表される場合を例示した。つまり、ＡｌのＴＭＰおよびＭＳが、Ａｌの位置（マクロ的な位置）によらず一定であるとして取り扱う場合を例示した。

　但し、ＴＭＰ（ｉ，ｔ）およびＭＳ（ｉ，ｔ）を、Ａｌを空間的に分割することによって得られた複数の部分領域（エリア）ごとに、別の情報として取り扱うこともできる。つまり、ＡｌのＴＭＰおよびＭＳが、Ａｌの位置（以下、ｘ）により異なるものとして取り扱うこともできる。

　具体的には、加工熱処理条件情報をＴＭＰ（ｉ，ｘ，ｔ）、材料組織情報をＭＳ（ｉ，ｘ，ｔ）として表現できる。すなわち、ｘを示す情報（位置情報）をさらに付加してＴＭＰおよびＭＳを表現してよい。当該構成によれば、部分領域ごとにＡｌのＴＭＰおよびＭＳを考慮できるので、ＡｌのＭＳをさらに高精度に予測することが可能となる。

　一例として、Ａｌの板材において、板の表面から中心までの領域を複数のエリアに分割し、それぞれエリアにおけるＭＳの変化を算出できる。具体的には、既知のＰＣ（ｉ）を用いて、「ＰＣ（ｉ）→ＴＭＰ（ｉ，ｘ，ｔ）→ＭＳ（ｉ，ｘ，ｔ）」の順に、ＭＳ（ｉ，ｘ，ｔ）を算出できる。つまり、複数のエリアごとに、実施形態１と同様に、ＴＭＰおよびＭＳを算出すればよい。

　一例として、複数のエリアごとにＴＭＰを算出することにより、ＰＣ（ｉ）からＴＭＰ（ｉ，ｘ，ｔ）を算出できる。この場合、上述のように、圧延理論またはＦＥＭを用いて、複数のエリアごとにＴＭＰを算出すればよい。

　〔実施形態２〕
　実施形態１では、工業プロセスにおけるＡｌの材料組織を予測する手法について説明を行った。但し、本発明の一態様に係る情報処理装置における計算手法の適用対象は、Ａｌのみに限定されない。

　例えば、当該計算手法は、工業プロセスにおける所定の金属製品の材料組織を予測するために適用されてもよい。一例として、当該計算手法は、工業プロセスにおける鉄（Ｆｅ）または銅（Ｃｕ）の材料組織を予測するために適用されてもよい。

　〔ソフトウェアによる実現例〕
　情報処理装置１の制御ブロック（特に制御部１０）は、集積回路（ＩＣチップ）等に形成された論理回路（ハードウェア）によって実現してもよいし、ＣＰＵ（Central Processing Unit）を用いてソフトウェアによって実現してもよい。

　後者の場合、情報処理装置１は、各機能を実現するソフトウェアであるプログラムの命令を実行するＣＰＵ、上記プログラムおよび各種データがコンピュータ（またはＣＰＵ）で読み取り可能に記録されたＲＯＭ（Read Only Memory）または記憶装置（これらを「記録媒体」と称する）、上記プログラムを展開するＲＡＭ（Random Access Memory）などを備えている。そして、コンピュータ（またはＣＰＵ）が上記プログラムを上記記録媒体から読み取って実行することにより、本発明の目的が達成される。上記記録媒体としては、「一時的でない有形の媒体」、例えば、テープ、ディスク、カード、半導体メモリ、プログラマブルな論理回路などを用いることができる。また、上記プログラムは、該プログラムを伝送可能な任意の伝送媒体（通信ネットワークや放送波等）を介して上記コンピュータに供給されてもよい。なお、本発明の一態様は、上記プログラムが電子的な伝送によって具現化された、搬送波に埋め込まれたデータ信号の形態でも実現され得る。

　〔まとめ〕
　本発明の態様１に係る材料組織計算装置は、第１工程から第Ｎ工程までを含む、アルミニウムを製造する工業プロセスについて、ｉを、１≦ｉ≦Ｎを満たす自然数、第ｉ工程の時刻ｔについて、ｔ＝０を当該第ｉ工程の開始時刻、ｔ＝ｔｆｉを当該第ｉ工程の終了時刻、第ｉ工程における上記アルミニウムの製造条件を示す情報をＰＣ（ｉ）、第ｉ工程の時刻ｔにおける上記アルミニウムの加工熱処理条件を示す情報をＴＭＰ（ｉ，ｔ）、第ｉ工程の時刻ｔにおける上記アルミニウムの材料組織を示す情報をＭＳ（ｉ，ｔ）、として、ＰＣ（１）～ＰＣ（Ｎ）およびＭＳ（１，０）を、あらかじめ設定された情報として取得する工程間情報統合部と、第１工程計算部から第Ｎ工程計算部までを含む工程計算部と、ＭＳ（ｉ，ｔ）およびＴＭＰ（ｉ，ｔ）に基づいて、ＭＳ（ｉ，ｔ）の時間変化を算出する材料組織計算部と、を備えており、上記材料組織計算部は、複数の冶金的現象のそれぞれの時間変化を算出する、複数の計算モジュールを含み、上記工程計算部における第ｉ工程計算部は、ＰＣ（ｉ）に基づいてＴＭＰ（ｉ，ｔ）を算出し、複数の上記計算モジュールのそれぞれには、共通のデータ構造を有するＭＳ（ｉ，ｔ）が供給され、ｉ＝１からｉ＝Ｎまでの順に、上記工程間情報統合部は、上記第ｉ工程計算部にＰＣ（ｉ）およびＭＳ（ｉ，０）を供給し、第ｉ工程において、ｔ＝０からｔ＝ｔｆｉまでの順に、上記第ｉ工程計算部は、（ｉ）上記材料組織計算部にＭＳ（ｉ，ｔ）およびＴＭＰ（ｉ，ｔ）を供給することで、当該材料組織計算部にＭＳ（ｉ，ｔｆｉ）を算出させ、かつ、（ｉｉ）当該ＭＳ（ｉ，ｔｆｉ）を上記工程間情報統合部に供給し、上記工程間情報統合部は、上記第ｉ工程計算部から取得したＭＳ（ｉ，ｔｆｉ）を、ＭＳ（ｉ＋１，０）として設定する。

　上述のように、従来の技術では、複数の製造工程を含む工業プロセスにおける、Ａｌ（アルミニウム）のＭＳ（材料組織）の変化を取り扱うことが困難であった。しかしながら、上記の構成によれば、種々のＰＣ（ｉ）（製造条件）に対して、ＭＳ（ｉ，ｔ）を算出できる。つまり、各工程の進行に伴うＭＳの変化を予測できる。

　それゆえ、Ａｌの製造時に得られるＭＳ（つまり、ＭＳ（Ｎ，ｔｆＮ））を従来よりも高精度に予測することが可能となる。すなわち、工業プロセスにおけるＡｌのＭＳを従来よりも高精度に予測できる。その結果、本発明の一態様に係る材料組織計算装置によれば、工業プロセスにおける試作検討のコストの回数を大幅に減少させることができる。

　また、本発明の一態様に係る材料組織計算装置によれば、第１工程計算部から第Ｎ工程計算部までのそれぞれ（第１工程から第Ｎ工程までのそれぞれ）において、共通のデータ構造を有するＭＳ（ｉ，ｔ）を用いて、材料組織計算部に計算を行わせることができる。つまり、共通の材料組織計算部を用いて計算を行うことができる。それゆえ、工程に応じて個別の材料組織計算部を準備する必要がない。

　このように、本発明の一態様に係る計算システムによれば、検討すべき工程に応じたカスタマイズが可能であるので、計算システムの開発工数を大幅に低減できるという利点も得られる。

　本発明の態様２に係る材料組織計算装置は、上記態様１において、上記第ｉ工程計算部は、ＰＣ（ｉ）に基づいて、第ｉ工程の時刻ｔにおいて、上記アルミニウムに加えられる温度、歪量、および歪速度のそれぞれを示す情報として、ＴＭＰ（ｉ，ｔ）を算出してよい。

　上述のように、温度、歪量、および歪速度はいずれも、ＡｌにおけるＭＳの変化を高精度に予測するために重要なパラメータである。それゆえ、上記の構成によれば、工業プロセスにおけるＡｌのＭＳをさらに高精度に予測できる。

　本発明の態様３に係る材料組織計算装置は、上記態様１または２において、Δｔを微小時間として、上記材料組織計算部は、ＭＳ（ｉ，ｔ）およびＴＭＰ（ｉ，ｔ）に基づいて、複数の相互作用的な冶金的現象の時間変化を算出して、ＭＳ（ｉ，ｔ＋Δｔ）を得ることを特徴とする請求項１または２に記載の材料組織計算装置。

　上記の構成によれば、複数の相互作用的な冶金的現象に伴うＭＳの時間変化を考慮できるので、工業プロセスにおけるＡｌのＭＳをさらに高精度に予測できる。

　本発明の態様４に係る材料組織計算装置は、上記態様１から３のいずれか１つにおいて、ＭＳ（ｉ，ｔ）には、第ｉ工程の時刻ｔにおける、添加元素の添加量、添加元素の固溶量、第二相粒子の粒子サイズ分布、再結晶率、および歪エネルギーのそれぞれを示す情報が含まれている。

　上記の構成によれば、ＭＳとして多様なパラメータを考慮できるので、工業プロセスにおけるＡｌのＭＳをさらに高精度に予測できる。

　本発明の態様５に係る材料組織計算装置は、上記態様１から４のいずれか１つにおいて、上記材料組織計算部には、上記材料組織の時間変化を算出するための計算モデルがあらかじめ設定されており、上記工程間情報統合部は、あらかじめ取得された上記材料組織の実測データとの誤差を低減するように、上記計算モデルに含まれるパラメータを調整してよい。

　上記の構成によれば、実測データに基づいて計算モデルに含まれるパラメータを調整できるので、工業プロセスにおけるＡｌのＭＳをさらに高精度に予測できる。

　本発明の態様６に係る材料組織計算装置は、上記態様１から５のいずれか１つにおいて、ＭＳ（ｉ，ｔ）は、上記アルミニウムの再結晶状態と未再結晶状態とで、別の情報として取り扱われていることが好ましい。

　上述のように、Ａｌの再結晶の有無によって、複数の冶金的現象のそれぞれの時間変化の態様は大きく異なりうる。それゆえ、上記の構成によれば、工業プロセスにおけるＡｌのＭＳをさらに高精度に予測できる。

　本発明の態様７に係る材料組織計算装置は、上記態様１から６のいずれか１つにおいて、上記材料組織計算部は、複数の上記計算モジュールとして、上記アルミニウムの回復の時間変化を算出する回復モジュールを、少なくとも備えていることが好ましい。

　上述のように、Ａｌは回復が生じやすい材料である。それゆえ、上記の構成によれば、回復モジュールによってＡｌの回復（Ａｌにおいて生じやすい冶金的現象）を考慮できるので、工業プロセスにおけるＡｌのＭＳをさらに高精度に予測できる。

　本発明の態様８に係る材料組織計算装置は、上記態様１から７のいずれか１つにおいて、ＴＭＰ（ｉ，ｔ）およびＭＳ（ｉ，ｔ）は、上記アルミニウムを空間的に分割することによって得られた複数の部分領域ごとに、別の情報として取り扱われていることが好ましい。

　上記の構成によれば、部分領域ごとにＴＭＰおよびＭＳを考慮できるので、工業プロセスにおけるＡｌのＭＳをさらに高精度に予測できる。

　本発明の各態様に係る材料組織計算装置は、コンピュータによって実現してもよく、この場合には、コンピュータを上記材料組織計算装置が備える各部（ソフトウェア要素）として動作させることにより上記材料組織計算装置をコンピュータにて実現させる材料組織計算装置の制御プログラム、およびそれを記録したコンピュータ読み取り可能な記録媒体も、本発明の範疇に入る。

　〔付記事項〕
　本発明は上述した各実施形態に限定されるものではなく、請求項に示した範囲で種々の変更が可能であり、異なる実施形態にそれぞれ開示された技術的手段を適宜組み合わせて得られる実施形態についても本発明の技術的範囲に含まれる。さらに、各実施形態にそれぞれ開示された技術的手段を組み合わせることにより、新しい技術的特徴を形成できる。

　〔本発明の別の表現〕
　なお、本発明の一態様は、以下のようにも表現できる。

　本発明の一態様に係る材料組織計算装置は、アルミニウム製品の製造工程において、（ｉ）初期組織に対して、材料に加えられた加工熱処理条件から材料組織情報の時間変化を計算する材料組織計算部と、（ｉｉ）前工程の材料組織情報と対象工程の製造方法および製造条件から、対象工程において生ずる材料組織情報の時間変化を、材料組織計算部を呼び出して計算し、工程終了後の材料組織情報を計算する工程計算部と、（ｉｉｉ）工程計算部が計算した材料組織情報を次工程に渡す工程間情報統合部と、を備え、上工程から順に生ずる材料組織変化を引き継いで、各工程におけるアルミニウムの材料組織情報を計算し、各工程および最終製品の材料組織を計算する。

　また、本発明の一態様に係る材料組織計算装置は、アルミニウム製品の製造工程において、（ｉ）対象工程の製造条件から、加工熱処理条件すなわち材料に加えられる温度、歪量、歪速度を計算し、計算された温度、歪量、歪速度と、前工程の材料組織情報から、対象工程において生ずる材料組織変化を計算する工程計算部と、（ｉｉ）材料組織を、材料組織情報すなわち主要添加元素の添加量、主要添加元素の固溶量、主要な第二相粒子のサイズ分布関数、再結晶率、歪エネルギーを含む材料組織パラメータ群で表現し、前工程の材料組織情報と対象工程の加工熱処理条件から、対象工程において、複数の冶金的現象が並列に進行することで生ずる材料組織変化を、複数の物理式あるいは経験式を組み合わせて計算する材料組織計算部と、をさらに備える。

　また、本発明の一態様に係る材料組織計算装置において、工程間情報統合部は、部分的に実測された各工程における材料組織因子を教師データとして、工程計算部および／あるいは材料組織計算部に含まれる係数を調整して予測精度を向上させる機能を有する。

　１　情報処理装置（材料組織計算装置）
　１１　工程間情報統合部
　１２　工程計算部
　１２ａ　第１工程計算部
　１２Ｎ　第Ｎ工程計算部
　１３　材料組織計算部
　１３０ａ　再結晶モジュール（計算モジュール）
　１３０ｂ　固溶析出モジュール（計算モジュール）
　１３０ｃ　加工組織モジュール（計算モジュール）
　１３０ｄ　回復モジュール（計算モジュール）
　ＰＣ（ｉ）　製造条件情報（第ｉ工程におけるＡｌの製造条件を示す情報）
　ＴＭＰ（ｉ，ｔ）　加工熱処理条件情報（第ｉ工程の時刻ｔにおけるＡｌの加工熱処理条件を示す情報）
　ＭＳ（ｉ，ｔ）　材料組織情報（第ｉ工程の時刻ｔにおけるＡｌの材料組織を示す情報）
　ｉ　工程の番号
　Ｎ　工程の数
　ｔ　時刻
　Δｔ　微小時間
　ｔｉ　第ｉ工程の時刻
　ｔｆｉ　第ｉ工程の終了時刻

Claims (9)

1. 　第１工程から第Ｎ工程までを含む、アルミニウムを製造する工業プロセスについて、
   　　ｉを、１≦ｉ≦Ｎを満たす自然数、
   　　第ｉ工程の時刻ｔについて、ｔ＝０を当該第ｉ工程の開始時刻、ｔ＝ｔｆｉを当該第ｉ工程の終了時刻、
   　　第ｉ工程における上記アルミニウムの製造条件を示す情報をＰＣ（ｉ）、
   　　第ｉ工程の時刻ｔにおける上記アルミニウムの加工熱処理条件を示す情報をＴＭＰ（ｉ，ｔ）、
   　　第ｉ工程の時刻ｔにおける上記アルミニウムの材料組織を示す情報をＭＳ（ｉ，ｔ）、
   　として、
   　ＰＣ（１）～ＰＣ（Ｎ）およびＭＳ（１，０）を、あらかじめ設定された情報として取得する工程間情報統合部と、
   　第１工程計算部から第Ｎ工程計算部までを含む工程計算部と
   　ＭＳ（ｉ，ｔ）およびＴＭＰ（ｉ，ｔ）に基づいて、ＭＳ（ｉ，ｔ）の時間変化を算出する材料組織計算部と、を備えており、
   　上記材料組織計算部は、複数の冶金的現象のそれぞれの時間変化を算出する、複数の計算モジュールを含み、
   　上記工程計算部における第ｉ工程計算部は、ＰＣ（ｉ）に基づいてＴＭＰ（ｉ，ｔ）を算出し、
   　複数の上記計算モジュールのそれぞれには、共通のデータ構造を有するＭＳ（ｉ，ｔ）が供給され、
   　　ｉ＝１からｉ＝Ｎまでの順に、
   　　上記工程間情報統合部は、上記第ｉ工程計算部にＰＣ（ｉ）およびＭＳ（ｉ，０）を供給し、
   　　　第ｉ工程において、ｔ＝０からｔ＝ｔｆｉまでの順に、
   　　　上記第ｉ工程計算部は、（ｉ）上記材料組織計算部にＭＳ（ｉ，ｔ）およびＴＭＰ（ｉ，ｔ）を供給することで、当該材料組織計算部にＭＳ（ｉ，ｔｆｉ）を算出させ、かつ、（ｉｉ）当該ＭＳ（ｉ，ｔｆｉ）を上記工程間情報統合部に供給し、
   　　上記工程間情報統合部は、上記第ｉ工程計算部から取得したＭＳ（ｉ，ｔｆｉ）を、ＭＳ（ｉ＋１，０）として設定することを特徴とする材料組織計算装置。
2. 　上記第ｉ工程計算部は、ＰＣ（ｉ）に基づいて、
   　第ｉ工程の時刻ｔにおいて、上記アルミニウムに加えられる温度、歪量、および歪速度のそれぞれを示す情報として、ＴＭＰ（ｉ，ｔ）を算出することを特徴とする請求項１に記載の材料組織計算装置。
3. 　Δｔを微小時間として、
   　上記材料組織計算部は、ＭＳ（ｉ，ｔ）およびＴＭＰ（ｉ，ｔ）に基づいて、
   　複数の相互作用的な冶金的現象の時間変化を算出して、ＭＳ（ｉ，ｔ＋Δｔ）を得ることを特徴とする請求項１または２に記載の材料組織計算装置。
4. 　ＭＳ（ｉ，ｔ）には、
   　第ｉ工程の時刻ｔにおける、添加元素の添加量、添加元素の固溶量、第二相粒子の粒子サイズ分布、再結晶率、および歪エネルギーのそれぞれを示す情報が含まれていることを特徴とする請求項１から３のいずれか１項に記載の材料組織計算装置。
5. 　上記材料組織計算部には、上記材料組織の時間変化を算出するための計算モデルがあらかじめ設定されており、
   　上記工程間情報統合部は、あらかじめ取得された上記材料組織の実測データとの誤差を低減するように、上記計算モデルに含まれるパラメータを調整することを特徴とする請求項１から４のいずれか１項に記載の材料組織計算装置。
6. 　ＭＳ（ｉ，ｔ）は、上記アルミニウムの再結晶状態と未再結晶状態とで、別の情報として取り扱われていることを特徴とする請求項１から５のいずれか１項に記載の材料組織計算装置。
7. 　上記材料組織計算部は、複数の上記計算モジュールとして、
   　上記アルミニウムの回復の時間変化を算出する回復モジュールを、少なくとも備えていることを特徴とする請求項１から６のいずれか１項に記載の材料組織計算装置。
8. 　ＴＭＰ（ｉ，ｔ）およびＭＳ（ｉ，ｔ）は、上記アルミニウムを空間的に分割することによって得られた複数の部分領域ごとに、別の情報として取り扱われていることを特徴とする請求項１から７のいずれか１項に記載の材料組織計算装置。
9. 　請求項１に記載の材料組織計算装置としてコンピュータを機能させるための制御プログラムであって、上記工程計算部、上記工程計算部、および上記材料組織計算部としてコンピュータを機能させるための制御プログラム。

Images