(57)【要約】【課題】鋳造条件が変動しても鋳物の品質を保証することができる鋳造方案設計方法を提供する。【解決手段】複数の要素から形成される解析モデルを用いた鋳物の伝熱凝固解析に基づく鋳造方案設計方法であって、前記解析モデルに対して複数水準の鋳物の鋳造条件を設定する段階（Ｓ３００）と、前記複数水準の鋳物の鋳造条件にそれぞれ対応する前記鋳物の熱物性値を算出する段階（Ｓ３１０）と、前記算出した熱物性値を用いて、前記複数水準の鋳造条件にそれぞれ対応する前記鋳物の伝熱凝固解析を実行する段階（Ｓ３２０）と、前記複数水準の鋳造条件にそれぞれ対応する伝熱凝固解析から得られた情報が、所定の品質許容範囲内であるか否かを判定する段階（Ｓ３４０）と、を含んでいる。【選択図】図３
2008155230.tif 000002
【特許請求の範囲】
【請求項１】
複数の要素から形成される解析モデルを用いた鋳物の伝熱凝固解析に基づく鋳造方案設計方法であって、
前記解析モデルに対して複数水準の鋳物の鋳造条件を設定する段階と、
前記複数水準の鋳物の鋳造条件にそれぞれ対応する前記鋳物の熱物性値を算出する段階と、
前記算出した熱物性値を用いて、前記複数水準の鋳造条件にそれぞれ対応する前記鋳物の伝熱凝固解析を実行する段階と、
前記複数水準の鋳造条件にそれぞれ対応する伝熱凝固解析から得られた情報が、所定の品質許容範囲内であるか否かを判定する段階と、
を含むことを特徴とする鋳造方案設計方法。
【請求項２】
前記所定の品質許容範囲に含まれるか否かを判定する段階の後に、
前記複数水準の鋳造条件にそれぞれ対応する伝熱凝固解析から得られた情報の少なくとも一つの情報が前記品質許容範囲外である場合、前記解析モデルを変更する段階をさらに含むことを特徴とする請求項１に記載の鋳造方案設計方法。
【請求項３】
前記解析モデルを変更する段階の後に、
前記変更した解析モデルに対して請求項１に記載の前記鋳物の熱物性値を算出する段階から前記所定の品質許容範囲内であるか否かを判定する段階を実行することを特徴とする請求項２に記載の鋳造方案設計方法。
【請求項４】
前記複数水準の鋳造条件にそれぞれ対応する前記鋳物の伝熱凝固解析を実行する段階は、
前記鋳物の溶湯の凝固前の熱物性値を用いて相互に隣接する前記要素間の伝熱計算を実行する段階と、
前記伝熱計算に基づき前記要素の冷却速度を算出する段階と、
算出した冷却速度に基づき、前記鋳物の溶湯が共晶凝固を開始するときの凝固時の熱物性値と前記鋳物の溶湯が共晶凝固を完了した後の凝固後の熱物性値とを算出する段階と、
前記冷却速度から算出した前記凝固時の熱物性値と前記凝固後の熱物性値とを用いて伝熱凝固解析を実行する段階と、
を含むことを特徴とする請求項１に記載の鋳造方案設計方法。
【請求項５】
前記伝熱計算に基づき前記要素の冷却速度を算出する段階は、
前記要素の温度が所定の温度に達したときに実行されることを特徴とする請求項４に記載の鋳造方案設計方法。
【請求項６】
前記所定の温度は、前記鋳物の溶湯がある一定以上の時間をかけて平衡凝固したときの共晶凝固開始温度であることを特徴とする請求項５に記載の鋳造方案設計方法。
【請求項７】
前記凝固時の熱物性値は、前記鋳物の共晶凝固開始温度と当該鋳物の凝固潜熱量とを含む熱物性値であることを特徴とする請求項４に記載の鋳造方案設計方法。
【請求項８】
前記凝固後の熱物性値は、前記算出した冷却速度と前記鋳物の凝固組織との関係から算出される熱物性値であることを特徴とする請求項４に記載の鋳造方案設計方法。
【請求項９】
前記冷却速度から算出した前記凝固時の熱物性値と前記凝固後の熱物性値とを用いて伝熱凝固解析を実行する段階は、
共晶凝固を開始した要素に対して前記凝固時の熱物性値を設定する段階と、
共晶凝固が完了した要素に対して前記凝固後の熱物性値を設定する段階と、
を含むことを特徴とする請求項４に記載の鋳造方案設計方法。
【請求項１０】
前記凝固時の熱物性値は、前記冷却速度から算出される共晶凝固開始温度を含み、
前記共晶凝固を開始した要素に対して前記凝固時の熱物性値を設定する段階は、
前記共晶凝固開始温度に達した要素に対して、前記凝固時の熱物性値を設定することを特徴とする請求項９に記載の鋳造方案設計方法。
【請求項１１】
複数水準の鋳物の鋳造条件を設定する段階は、
前記鋳物を形成する材料の成分、前記鋳物を形成する材料の溶湯温度、接種材の種類、接種材の添加量、およびフェーディング時間のうちの少なくともいずれか一つを前記鋳造条件として設定することを特徴とする請求項１に記載の鋳造方案設計方法。
【請求項１２】
前記所定の品質許容範囲内か否かを判定する段階の後に、
前記伝熱凝固解析によって得られた情報を加工して表示する段階をさらに含むことを特徴とする請求項１に記載の鋳造方案設計方法。
【請求項１３】
複数の要素から形成される解析モデルを用いた鋳物の伝熱凝固計算に基づく鋳造方案設計装置であって、
前記解析モデルに対して複数水準の鋳物の鋳造条件を設定する鋳造条件設定手段と、
前記複数水準の鋳物の鋳造条件に対応する前記鋳物の熱物性値を算出する熱物性値算出手段と、
前記算出した熱物性値を用いて、前記複数水準の鋳造条件にそれぞれ対応する前記鋳物の伝熱凝固解析を実行する伝熱凝固解析手段と、
前記複数水準の鋳造条件にそれぞれ対応する伝熱凝固解析から得られた情報が、所定の品質許容範囲内か否かを判定する判定手段と、
を備えることを特徴とする鋳造方案設計装置。
【請求項１４】
複数の要素から形成される解析モデルを用いた鋳物の伝熱凝固計算に基づく鋳造方案設計プログラムであって、
前記解析モデルに対して複数水準の鋳物の鋳造条件を設定する手順と、
前記複数水準の鋳物の鋳造条件に対応する前記鋳物の熱物性値を算出する手順と、
前記算出した熱物性値を用いて、前記複数水準の鋳造条件にそれぞれ対応する前記鋳物の伝熱凝固解析を実行する手順と、
前記複数水準の鋳造条件にそれぞれ対応する伝熱凝固解析から得られた情報が、所定の品質許容範囲内か否かを判定する手順と、
をコンピュータに実行させることを特徴とする鋳造方案設計プログラム。
【請求項１５】
請求項１４に記載の鋳造方案設計プログラムを記録したコンピュータ読み取り可能な記録媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、鋳造方案設計方法に係わり、詳細には、電子計算機によるシミュレーションを利用した鋳造方案設計方法に関する。
【背景技術】
【０００２】
健全で安価な鋳物を製造するためには、鋳造品の形状や鋳造の設計方針を事前に検討することが必要である。そのための手段の一つとして、電子計算機（以下、「コンピュータ」と称する）を用いた鋳造解析が広く用いられている。
【０００３】
鋳造解析の種類は、流動解析、変形解析、伝熱凝固解析など多岐に渡り、特に、伝熱凝固解析は、鋳物の機械的特性の予測（例えば、鋳物の硬度や強度）、引け巣の予測、または金属組織の不良による鋳造欠陥の予測をするための重要な解析方法である（例えば、特許文献１および特許文献２）。
【０００４】
一般的な伝熱凝固解析は、解析モデル（解析対象物となる鋳物や鋳型領域を複数の要素（以下、「セル」と称する）に分割して構成されるモデル）において、相互に隣接するセル間の微少時間における熱流束を計算し、これを順次、微少時間間隔ごとに繰り返し、解析モデルの温度変化を求める解析方法である。
【０００５】
各セル間の熱流束の計算は、セルごとに設定された熱物性値、例えば、凝固開始温度、凝固潜熱量、比熱、熱伝導率などを用いて行われる。
【０００６】
ところで、日常の操業において鋳造条件が変動する場合、これに応じて鋳物の熱物性値も変動することが考えられる。熱物性値の変動は、鋳物の凝固形態に直接影響を与えるので、鋳物の品質に重大な影響を及ぼす可能性がある。
【０００７】
しかしながら、従来の伝熱凝固解析方法では、鋳造条件の変動にかかわらず、鋳物の熱物性値として１種類の鋳物の金属材種につき１水準の値のみを用い、解析結果が目標とする鋳物の品質を満足するまで解析モデルの形状を変化させて伝熱凝固解析を繰り返し実行していた。
【特許文献１】特開平７−１１３７７１号公報
【特許文献２】特開平８−３１３４６４号公報
【発明の開示】
【発明が解決しようとする課題】
【０００８】
したがって、上記従来の伝熱凝固解析方法では、鋳造条件の変動に係わらず解析に用いた１水準の鋳物の熱物性値に対しての鋳物の品質が保証されるに留まり、鋳造条件の変動に伴い変化する熱物性値に対しての鋳物の品質を保証することができないという問題点がある。
【０００９】
本発明は、上記問題点を解決するために成されたものであり、鋳造条件が変動しても鋳物の品質を保証することができる鋳造方案設計方法を提供することを目的とする。
【課題を解決するための手段】
【００１０】
上記目的を達成するための本発明に係る鋳造方案設計方法は、複数の要素から形成される解析モデルを用いた鋳物の伝熱凝固解析に基づく鋳造方案設計方法であって、前記解析モデルに対して複数水準の鋳物の鋳造条件を設定する段階と、前記複数水準の鋳物の鋳造条件にそれぞれ対応する前記鋳物の熱物性値を算出する段階と、前記算出した熱物性値を用いて、前記複数水準の鋳造条件にそれぞれ対応する前記鋳物の伝熱凝固解析を実行する段階と、前記複数水準の鋳造条件にそれぞれ対応する伝熱凝固解析から得られた情報が、所定の品質許容範囲内であるか否かを判定する段階と、を含むことを特徴とする。
【００１１】
また、上記の目的を達成するために本発明に係る鋳物の鋳造方案設計装置は、複数の要素から形成される解析モデルを用いた鋳物の伝熱凝固計算に基づく鋳造方案設計装置であって、前記解析モデルに対して複数水準の鋳物の鋳造条件を設定する鋳造条件設定手段と、前記複数水準の鋳物の鋳造条件に対応する前記鋳物の熱物性値を算出する熱物性値算出手段と、前記算出した熱物性値を用いて、前記複数水準の鋳造条件にそれぞれ対応する前記鋳物の伝熱凝固解析を実行する伝熱凝固解析手段と、前記複数水準の鋳造条件にそれぞれ対応する伝熱凝固解析から得られた情報が、所定の品質許容範囲内か否かを判定する判定手段と、を備えることを特徴とする。
【００１２】
また、上記の目的を達成するために本発明に係る鋳物の鋳造方案設計プログラムは、複数の要素から形成される解析モデルを用いた鋳物の伝熱凝固計算に基づく鋳造方案設計プログラムであって、前記解析モデルに対して複数水準の鋳物の鋳造条件を設定する手順と、前記複数水準の鋳物の鋳造条件に対応する前記鋳物の熱物性値を算出する手順と、前記算出した熱物性値を用いて、前記複数水準の鋳造条件にそれぞれ対応する前記鋳物の伝熱凝固解析を実行する手順と、前記複数水準の鋳造条件にそれぞれ対応する伝熱凝固解析から得られた情報が、所定の品質許容範囲内か否かを判定する手順と、をコンピュータに実行させることを特徴とする。
【発明の効果】
【００１３】
以上のように構成された本発明に係る鋳物の鋳造方案設計方法、鋳造方案設計装置、および鋳造方案設計プログラムによれば、複数水準の鋳造条件にそれぞれ対応する鋳物の熱物性値を算出して伝熱凝固解析を実行している。この結果、鋳造条件の変動に伴い変化する鋳物の熱物性値に対応した鋳物の品質を保証することができ、精度の高い鋳造方案設計方法を提供することができる
【発明を実施するための最良の形態】
【００１４】
本発明者は、鋳物の熱物性値の変化が鋳物の溶湯が凝固するときの温度変化を示す冷却曲線から読み取れることに着目した。そして、この冷却曲線から、鋳造条件と鋳物の熱物性値との関係式を導出し、この関係式を伝熱凝固解析に適用することによって、複数水準の鋳造条件に対応した鋳物の品質を予測することが可能となることを見出した。
【００１５】
以下に、本発明に係る鋳造方案設計方法を第１実施形態および第２実施形態に分けて、図面を参照しながら詳細に説明する。なお、同じ構成要素には、同じ参照番号を付してある。
【００１６】
まず、本発明に係る鋳造方案設計方法を実施するための鋳造方案設計装置について詳細に説明しておく。
【００１７】
図１は、本発明に係る鋳造方案設計方法を実施するための鋳造方案設計装置の概略構成を示すブロック図である。図１に示すように、本実施の形態における鋳造方案設計装置（以下、「コンピュータ」と称する）１００は、ＣＰＵ１１１、ＲＡＭ１１３、ＲＯＭ１１５、ハードディスク１１７、ディスプレイ１２１、および入力部１３１を備えており、これらは信号をやり取りするためのバス１４１を介して相互に接続されている。
【００１８】
コンピュータ１００は、例えばパーソナルコンピュータなどの電子計算機であり、伝熱凝固解析プログラム（以下、「シミュレーションプログラム」と称する）に基づいて伝熱凝固解析を実行したり、解析結果から得られた情報を加工して表示したりするものである。解析結果から得られる情報は、例えばセルの温度の経時的変化、共晶凝固開始温度、凝固潜熱量、冷却速度、または凝固時間などであり、解析結果から得られる各種の情報（以下、これらの情報を「欠陥判定パラメータ」と称する）が該当する。そして、欠陥判定パラメータに基づき鋳物の比較・評価、例えば製品として合格品であるか否か（鋳物の品質が保証されているか否か）などの判断が行われる。
【００１９】
ＣＰＵ１１１は、シミュレーションプログラムに基づいて上記各部の制御や伝熱凝固解析に必要な各種演算処理を実行するものである。ＣＰＵ１１１は、鋳造条件設定部（鋳造条件設定手段）、熱物性値算出部（熱物性値算出手段）、伝熱凝固解析部（および伝熱計算手段および伝熱凝固解析手段）、判定部（判定手段）、冷却速度算出部（冷却速度算出手段）として機能する。
【００２０】
ＲＡＭ１１３は、作業領域として一時的にプログラムやデータを記憶するものである。
【００２１】
ＲＯＭ１１５は、コンピュータ１００の基本動作を制御する各種プログラムやパラメータなどをあらかじめ格納しているものである。
【００２２】
ハードディスク１１７は、ＯＳ（オペレーティングシステム）やコンピュータ１００の所定の動作を制御するためのプログラムやパラメータを格納するものである。ハードディスク１１７には、解析モデルの作成および変更、各種熱物性値の算出、解析結果から得られた情報の加工または表示、欠陥判定パラメータに基づく完成品の比較・評価を実行するプログラムなど、その他、一般的な伝熱凝固解析に必要なプログラムがあらかじめ記憶されている。さらに、ハードディスク１１７は、解析結果を記憶する記憶領域としても機能する。なお、伝熱凝固解析に必要なプログラムは、記録媒体（例えば、ＣＤ−ＲＯＭ、ＤＶＤ−ＲＯＭなどの記録媒体）にあらかじめ記憶されているものでも良く、この記録媒体からプログラムを直接読み取って伝熱凝固解析をコンピュータ１００に実行させても良い。
【００２３】
ディスプレイ１２１は、例えばＣＲＴディスプレイや液晶ディスプレイなどであり、解析結果によって得られた各種の情報を表示するものである。
【００２４】
入力装置１３１は、マウス、キーボード、またはタッチパネルなどのポインティングデバイスであり、ユーザーからの入力を受け付けるものである。
【００２５】
以上のように構成されたコンピュータ１００を用いて、本発明に係る鋳造方案設計方法が実行される。
【００２６】
以下に、本発明の第１実施形態に係る鋳造方案設計方法について図面を参照しながら詳細に説明する。
【００２７】
［第１実施形態］
図２から図９は、本発明の第１実施形態に係る鋳物の鋳造方案設計方法および鋳造方案設計装置の説明に供する図である。図２は、本実施の形態の鋳造方案設計方法によって設計される鋳物を例示したものを、図３は図１に示した鋳造方案設計装置の動作フローチャートを示し、このフローチャートは、本実施の形態に係る鋳造方案設計方法の手順に相当するものである。図４〜図９は本実施の形態に係る鋳造方案設計方法の処理内容を説明するために用いられる図である。
【００２８】
図２を参照しつつ、本実施の形態の鋳造方案設計方法で設計される鋳物について説明する。図２は、本実施の形態の鋳造方案設計方法によって設計される鋳物２００およびチラー２１０を示したものである。図２（Ａ）に示すように、鋳物２００は、車両用鋳鉄部品であって、ベース部２００ａとノーズ部２００ｂとを備えるカム形状の部材である。図２（Ｂ）に示すように、鋳物２００のベース部２００ａの肉厚は、ノーズ部２００ｂの肉厚よりも薄く設計されている。したがって、ベース部２００ａの部位Ａは軸孔加工されるため比較的低い硬度が求められる部位であり、ノーズ部２００ｂの部位Ｂは耐摩耗性から比較的高い硬度が求められる部位である。本実施の形態では、鋳物２００およびチラー２１０を解析モデルとしている。
【００２９】
また、本実施の形態では、鋳造条件として、鋳物の溶湯が溶解炉から注湯取鍋に出湯された時刻からの経過時間を示すフェーディング時間を用い、２水準のフェーディング時間に対して鋳造される鋳物２００が所定の品質保証範囲内に含まれるか否かを判定する場合を例に挙げて説明する。以下、本実施の形態に係る鋳造方案設計方法を詳細に説明する。
【００３０】
図３に示すように、まず、解析モデルに対して複数水準の鋳造条件を設定する（Ｓ３００）。
）。具体的には、解析対象物となる鋳物の複数の要素から形成される解析モデルを作成し、この解析モデルに対して、鋳造に関する複数水準の鋳造条件を設定する。
【００３１】
解析モデルの作成は、解析モデルの作成に必要なデータ、例えば形状データ（解析対象となる鋳物の形状、鋳物の設計形状、鋳型の形状など、鋳物の凝固解析を行うために必要な形状データ）、要素分割数、鋳造時間、鋳物を形成する材料、鋳物の溶湯の温度など、一般的なシミュレーションに必要な各種データを入力装置１３１からコンピュータ１００へ入力することによって開始され、入力された各種データに基づきコンピュータ１００が演算することによって終了する。解析モデルはメッシュモデルとして作成される。初期条件の設定は、入力装置１３１から各種データをコンピュータ１００へ入力することによって行なわれる。なお、解析モデルの作成に必要なデータや初期条件、鋳造条件などは、あらかじめハードディスク１１７などに記憶させておいても良いし、作業者が手入力しても良い。また、要素分割数はメッシュモデルのセル数（要素数）に等しく、要素分割はシミュレーションを行う際のメッシュモデルに対して行われる。セルとはメッシュモデルの各要素を指す。
【００３２】
前記鋳造条件の設定は、鋳物を形成する材料の成分、鋳物を形成する材料の溶湯温度、接種材の種類、接種材の添加量、またはフェーディング時間のうちの少なくともいずれか一つを鋳造条件として設定する。
【００３３】
本実施の形態では、複数水準の鋳造条件として、フェーディング時間が０分（以下、「フェーディング時間０分」と称する）およびフェーディング時間が１５分（以下、「フェーディング時間１５分」と称する）の２つの水準の鋳造条件を設定している。既に説明したように、フェーディング時間とは、鋳物の溶湯が溶解炉から注湯取鍋に出湯した時刻を起点とした経過時間を意味し、このフェーディング時間が長いほど鋳物の溶湯の黒鉛化能は低下する。通常、生産工場におけるフェーディング時間は最長１５分であることから、鋳造条件としてフェーディング時間０分およびフェーディング時間１５分の２つの水準を採用している。
【００３４】
次に、設定された複数水準の鋳造条件のうち、一の鋳造条件を選択する（Ｓ３１０）。具体的には、鋳造条件であるフェーディング時間０分およびフェーディング時間１５分のうちから一の鋳造条件を選択、例えば、フェーディング時間０分を選択する。
【００３５】
次に、選択した鋳物の鋳造条件に対応する熱物性値を算出する（Ｓ３２０）。熱物性値の算出は、上述したように、ハードディスク１１７にあらかじめ記憶されている鋳造条件と鋳物の熱物性値との関係式（例えば、鋳造条件に対応した熱物性値を算出する変換テーブル）を用いて、選択した鋳造条件に対応する鋳物の熱物性値を算出する。本実施の形態では、ステップＳ３１０の処理で選択された鋳造条件であるフェーディング時間０分に対応する熱物性値が算出される。
【００３６】
以下に、鋳造条件に対応する熱物性値の算出方法を詳細に説明する。本実施の形態では、鋳造条件としてフェーディング時間を採用しているので、フェーディング時間と熱物性値の一つである共晶凝固開始温度との関係を例に挙げて説明する。
【００３７】
図４は、フェーディング時間による鋳物の冷却曲線の変化を示したものである。横軸は時間、縦軸は温度を示している。冷却曲線Ａ０はフェーディング時間０分の冷却曲線であり、このときの共晶凝固開始温度はＴ０である。冷却曲線Ａ１５はフェーディング時間１５分の冷却曲線であり、このときの共晶凝固開始温度はＴ１５である。図４に示すように、共晶凝固開始温度は、フェーディング時間が長くなるほど低下していく。これは、鋳造条件を変更した場合、鋳造条件の変更に伴い冷却曲線が変化するので、冷却曲線に対応する熱物性値も変化することを示している。
【００３８】
そして、これらの冷却曲線を数値化して処理することにより、フェーディング時間と共晶凝固開始温度との関係式を導出することができる。他の熱物性値も、同様にして導出することができる。
【００３９】
次に、ステップＳ３２０の処理で算出した熱物性値を用いて、伝熱凝固解析を実行する（Ｓ３３０）。伝熱凝固解析の基本原理は、解析モデル（解析対象物となる鋳物や鋳型領域などを複数のセルに分割して構成されるモデル）において、セルの特徴に応じた熱物性値を設定し、相互に隣接するセル間の微少時間間隔における熱流束を計算する。そして、順次、微少時間間隔ごとに熱流束の計算を繰り返して解析モデルの温度変化を求める。
【００４０】
熱流束の計算は、一般的な伝熱計算の手法を用いて実行される。ここで、セルの特徴に応じた設定とは、鋳物（例えば、図２に示す鋳物２００）に属するセルに対しては所定の熱物性値が設定され、それ以外のセル、すなわち、鋳型領域（例えば、図２に示すチラー３００）に属するセルに対してはそれぞれが属する材料の熱物性値が設定されることを意味する。熱物性値としては、例えば凝固開始温度、液相線温度、比熱、熱伝導率などの伝熱計算に必要な物性値である。
【００４１】
なお、解析開始時においては、総てのセルに所定の初期温度が与えられているが、伝熱計算により、微少時間間隔後には初期温度は変化していく。解析開始時では、鋳物に属するセルの温度は、鋳型領域に属するセルの温度よりも高く設定されている。このため、伝熱計算の結果として鋳物に属するセルの温度は、当然、低下していくことになる。
【００４２】
本実施の形態では、解析結果から得られた温度変化のデータを数値化して処理し、セルごとの冷却速度を算出している。冷却速度は、以下のように算出される。
【００４３】
セルが共晶凝固開始温度ＴＬ達した時点から微少時間間隔である‘Δｔ’時間だけ前の当該セルの温度をＴとすると、冷却速度‘ｄＴ／ｄｔ’は、下記の式（１）で表せる。
【００４４】
【数１】
2008155230.tif 000003
【００４５】
このように、冷却速度を算出する理由は、冷却速度は鋳物の品質の指針となる鋳物の機械的特性、例えば、鋳物の硬度に影響を及ぼす重要な因子だからである（冷却速度と鋳物の硬度との関係の詳細な説明は後述する）。
【００４６】
次に、総ての鋳造条件に対して伝熱凝固解析を実行したか否かを判定する（Ｓ３４０）。総ての鋳造条件に対して伝熱凝固解析を実行していない場合（Ｓ３４０：ＮＯ）、他の鋳造条件が選択され（Ｓ３７０）、他の鋳造条件に基づきステップＳ３２０以下の処理を再び実行する。具体的には、フェーディング時間１５分を選択し、ステップＳ３１０以下の処理を再び実行する。
【００４７】
そして、総ての鋳造条件に対して伝熱凝固解析を実行した場合（Ｓ３４０：ＹＥＳ）、当該鋳造条件に対応する伝熱凝固解析から得られた情報が、所定の品質許容範囲内であるか否かを判定する（Ｓ３５０）。具体的には、ステップＳ３３０の処理で算出されたセルごとの冷却速度に基づいて鋳物の硬度が算出され、算出された硬度が所定の品質許容範囲内（以下、「規格値」と称する）であるか否かを判定している。以下に、冷却速度と鋳物の硬度との関係およびステップＳ３５０の処理内容について詳細に説明する。
【００４８】
図５は、冷却速度と鋳物の硬度との関係を示したものである。横軸は、図２に示す鋳物２００の部位Ｃからの距離、縦軸は図２に示す鋳物２００の硬度を示している。
【００４９】
図５に示すように、冷却速度と鋳物の硬度とはほぼ比例関係を有している。
【００５０】
前記規格値は、例えば、部位Ｃからの距離ＬＡにおける部位Ａ（図２（Ａ）参照）の硬度の規格値はＨＡ以下であり、部位Ｃからの距離ＬＢにおける部位Ｂ（図２（Ａ）参照）の硬度の規格値はＨＢ以上であると定める。なお、既に説明したように、部位Ａは比較的低い硬度が要求される部位であり、部位Ｂは比較的高い高度が要求される部位である。
【００５１】
この場合、図５に示す関係から、部位Ａの冷却速度はＶＡ以下であり、部位Ｂの冷却速度はＶＢ以上を規格値として要求されていることが分かる。なお、冷却速度は、鋳物の硬度に影響を及ぼす因子であるため、本実施の形態では、冷却速度と鋳物の硬度との関係から、規格値を定めている。
【００５２】
図６は、鋳物の先端部（図２における部位Ｃ）からの距離と冷却速度との関係を示したものである。横軸は、図２に示す部位Ｃからの距離、縦軸は冷却速度を示している。また、曲線Ｃ０はフェーディング時間０分、曲線Ｃ１５はフェーディング時間１５分の場合を示している。上述したように、部位Ａの冷却速度はＶＡ以下、部位Ｂの冷却速度はＶＢ以上が規格値として要求されているが、図６に示すように、フェーディング時間０分の場合、部位Ａの冷却速度は冷却速度ＶＡ以下の規格値を満たす一方、部位Ｂの冷却速度は冷却速度ＶＢ以上の規格値を満たしていない。
【００５３】
また、フェーディング時間１５分の場合、部位Ｂの冷却速度は冷却速度ＶＢ以上の規格値を満たす一方、部位Ａの冷却速度は冷却速度ＶＡ以下の規格値を満たしていない。ステップＳ３５０の処理では、このようにして、冷却速度と鋳物の硬度との関係から規格値内であるか否かを判定するようにしている。
【００５４】
再び、ステップＳ３５０の処理内容の説明に戻る。上述のように、少なくとも一つの解析結果が規格値内でない場合（Ｓ３５０：ＮＯ）、解析モデルを変更する（Ｓ３６０）。解析モデルの変更は、例えば、図７に示すように、鋳鉄部品２００のベース部２００ａの肉厚が大きくすると共に、チラー２１０の底部の容量が拡大するように変更する。なお、解析モデルの変更は、ユーザーが手動で実行しても良いし、ハードディスク１１７にあらかじめ記憶された解析モデルを自動的に読み込んで変更しても良い。なお、図７において図２と同じ構成要素には、同じ参照番号を付してある。そして、変更後の解析モデルを用いて、ステップＳ３１０以下の処理を繰り返す。
【００５５】
そして、総ての鋳造条件に対応する伝熱凝固解析から得られた結果が、規格値内であると判定した場合、処理を終了する（Ｓ３５０：ＹＥＳ）。具体的には、以下のような処理が実行されている。
【００５６】
図８は解析モデルの変更後の鋳物の先端部（図７における部位Ｃ）からの距離と冷却速度との関係を示したものである。横軸は、図７に示す部位Ｃからの距離、縦軸は冷却速度を示している。また、曲線Ｃ’０はフェーディング時間０分、曲線Ｃ’１５はフェーディング時間１５分の場合を示している。
【００５７】
図８に示すように、フェーディング時間０分の場合、部位Ａの冷却速度は冷却速度ＶＡ以下の規格値を満たし、部位Ｂの冷却速度は冷却速度ＶＢ以上の規格値を満たしている。
【００５８】
また、図７に破線で示されるフェーディング時間１５分の場合、部位Ｂの冷却速度は冷却速度ＶＢ以上の規格値を満たし、部位Ａの冷却速度は冷却速度ＶＡ以下の規格値を満たしている。
【００５９】
このように、規格値を満たすまで解析モデルの形状を変更し、鋳造条件のフェーディング時間０分とフェーディング時間１５分の２水準に対する解析結果が規格値を満たせば処理を終了する。
【００６０】
以上のように構成された本発明に係る鋳物の鋳造方案設計方法およびその方案設計装置によれば、複数水準の鋳造条件に対応した熱物性値が算出され、算出された熱物性値を用いて伝熱凝固解析が実行される。そして、解析結果に基づいて鋳物の品質が判定され、鋳物の品質が所定の規格値内であると判定されるまで解析モデルが変更され、伝熱凝固解析が繰り返し実行される。この結果、鋳造条件の変動に伴い変化する鋳物の熱物性値に対応した鋳物の品質を保証することができ、精度の高い鋳造方案設計方法を提供することができる。
【００６１】
［第２実施形態］
第１実施形態では、鋳物の熱物性値が鋳造条件によって変化することを考慮することで、従来の方法よりも、精度の高い解析結果を得ることができることを説明した。鋳造条件は、主として鋳型に注湯されるまでの因子であるが、鋳型に注湯された後、鋳物の熱物性値に影響を及ぼす因子もある。鋳型に注湯された後、鋳物の熱物性値に影響を及ぼす因子として鋳物の溶湯の冷却速度が挙げられる。そこで、本実施の形態では、鋳造条件に応じて熱物性値が変化することに加えて、鋳物の溶湯の冷却速度に応じて鋳物の熱物性値が変化する場合について述べる。
【００６２】
図９から図１２は、本発明の実施の形態に係る鋳物の伝熱凝固解析方法の説明に供する図である。図９および図１０は図１に示した鋳造方案設計装置の動作フローチャートを示し、このフローチャートは、本実施の形態に係る鋳造方案設計方法の手順に相当するものである。図１１から図１２は、本実施の形態に係る鋳造方案設計方法における熱物性値の算出方法を説明するために用いられる図である。なお、本実施の形態において第１実施形態と同じ鋳造条件で処理を実行している。したがって、第１実施形態と同じ処理内容については、重複記載を避けるため詳細な説明は省略する。
【００６３】
図９に示すように、まず、解析モデルに対して複数水準の鋳造条件を設定する（Ｓ４００）。
【００６４】
次に、設定された複数水準の鋳造条件のうち一の鋳造条件を選択する（Ｓ４１０）。
【００６５】
次に、選択した鋳造条件に対応する熱物性値を算出する（Ｓ４２０）。ステップ４２０の処理において算出される熱物性値は、図３のステップＳ３２０の処理と同様に算出される熱物性値であるが、これは、鋳造条件であるフェーディング時間のみを考慮した熱物性値であるので、本実施の形態では、鋳物の溶湯の凝固前の熱物性値（以下、「凝固前の熱物性値」と称する）を意味する。
【００６６】
次に、ステップＳ４２０の処理で算出された熱物性値を用いて、伝熱凝固解析処理が実行される（Ｓ４３０）。ここで、ステップＳ４３０の処理内容について詳細に説明する。
【００６７】
図１０は、ステップＳ４３０のサブルーチンの処理内容を示すフローチャートである。
【００６８】
図１０に示すように、まず、解析モデルに対してステップＳ４２０の処理で算出された凝固前の熱物性が設定され、相互に隣接するセル間の伝熱計算を実行する（Ｓ４３１）。凝固前の熱物性値の設定はセルの特徴に応じて設定し、伝熱計算を実行する。
【００６９】
次に、伝熱計算の結果に基づきセルの温度が所定の温度（以下、「冷却速度算出温度」と称する）に達したか否かを判定する（Ｓ４３２）。セルが冷却速度算出温度に達しない場合、ステップＳ４３１以下の処理を繰り返し（Ｓ４３２：ＮＯ）、セルの温度が冷却速度算出温度に達した場合（Ｓ４３２：ＹＥＳ）、ステップＳ４３３の処理に進む。ステップＳ４３２の処理では、セルごとに冷却速度算出温度に達したか否かを判定している。
【００７０】
次に、冷却速度算出温度に達したセルに対して、当該セルの冷却速度を算出する（Ｓ４３３）。以下に、ステップＳ４３３における冷却速度の算出方法を詳細に説明する。
【００７１】
セルの温度が所定の温度に達した時点におけるセルの温度をＴｎ、当該セルが所定の温度に達した時点から微少時間間隔である‘Δｔ’時間だけ前の当該セルの温度をＴｎ−１とすると、冷却速度‘ｄＴ／ｄｔ’は、下記の式（２）で表される。
【００７２】
【数２】
2008155230.tif 000004
【００７３】
なお、冷却速度算出温度は任意の温度に設定することができるが、設定する温度は、実際の鋳造における共晶凝固開始温度以上であることが好ましい。しかしながら、現時点における解析上では、実際の鋳造における共晶凝固開始温度は未知である。そこで、冷却速度算出温度としては、鋳物の溶湯が無限遠の時間をかけて平衡凝固したときの共晶凝固開始温度とすることが好ましく、特に鋳物を形成する材料が鋳鉄の場合には、鋳物溶湯中の炭素が完全黒鉛化したときの共晶凝固開始温度とすることが好ましい。このような温度を設定すれば、解析モデルのうち、どのセルにおいても冷却速度算出温度が実際の鋳造における共晶凝固開始温度より低くなることが避けられるからである。
【００７４】
次に、算出した冷却速度に基づき、前記鋳物の溶湯の共晶凝固が開始されるときの凝固時の熱物性値（以下、「凝固時の熱物性値」と称する）と前記鋳物の溶湯の共晶凝固が完了した後の凝固後の熱物性値（以下、「凝固後の熱物性値」と称する）とを算出する（Ｓ４３４）。凝固時の熱物性値および凝固後の熱物性値は、ステップＳ４３３の処理から求めた冷却速度を用いてセルごとに算出される。
【００７５】
以下に、凝固時の熱物性値および凝固後の熱物性値の算出方法を詳細に説明する。凝固時の熱物性値として、共晶凝固開始温度、凝固潜熱量の算出方法を例に挙げて説明する。
【００７６】
凝固時の熱物性値のひとつである共晶凝固開始温度ＴＬおよび凝固潜熱量ｑは、冷却速度‘ｄＴ／ｄｔ’の関数として、以下の式で表すことができる。
【００７７】
【数３】
2008155230.tif 000005
【００７８】
【数４】
2008155230.tif 000006
【００７９】
ここで、ＴＬｍａｘは鋳物溶湯中の炭素元素が完全黒鉛化した場合の共晶凝固開始温度である。
【００８０】
以下に、本実施の形態で用いられる式（３）、（４）の算出方法を詳細に説明する。
【００８１】
図１１は、人為的に冷却速度を変化させて鋳物を鋳造した場合における冷却曲線の概略図である。横軸は時間、縦軸は鋳物の温度を示している。この冷却曲線は、実際に鋳物の鋳造を繰り返して測定した結果を示したものである。このように冷却速度を変化させて測定する理由は、既に説明したように、鋳造条件を変更すれば冷却曲線も変化し、結果的にこれに対応する熱物性値も変化するからである。
【００８２】
冷却曲線Ｂ１は鋳物を徐冷して炭素元素が完全黒鉛化したときの冷却曲線であり、このときの共晶凝固開始温度はＴ３である。共晶凝固開始温度Ｔ３はＴＬｍａｘと等しい。冷却曲線Ａ２は人為的に冷却速度を変えなかったときの冷却曲線であり、このときの共晶凝固開始温度はＴ２である。冷却曲線Ａ３は鋳物を急冷して完全レデブライト化したときの冷却曲線であり、このときの共晶凝固開始温度はＴ１である。図に示すように、冷却曲線Ｂ２は冷却曲線Ｂ１〜Ｂ３により囲まれる領域に存在している。このとき、共晶凝固開始温度Ｔ５は共晶凝固開始温度Ｔ４〜Ｔ６の範囲に存在している。
【００８３】
図１２は、図１１に示した冷却曲線Ｂ１〜Ｂ３の測定結果から共晶凝固開始温度と共晶凝固開始時の冷却速度との関係を求めた結果を示したものである。横軸は共晶凝固開始時の冷却速度、縦軸は共晶凝固開始温度を示している。
【００８４】
図１２に示すように、共晶凝固開始温度は共晶凝固開始時の冷却速度に反比例して低下している（曲線Ｂ’１〜Ｂ’３参照）。一の鋳造条件における鋳物の冷却速度と鋳物の共晶凝固開始温度の関係を示す曲線がＢ’１である場合、鋳造条件を変化させれば、鋳物の冷却速度と鋳物の共晶凝固開始温度の関係は、例えば曲線Ｂ’２、または、曲線Ｂ’３に変化する。具体的には、例えば鋳物の冷却速度がＶＣである場合、曲線Ｂ’１における共晶凝固開始温度はＴ７であるので、鋳造条件を変更すれば、曲線Ｂ’２における共晶凝固開始温度はＴ６に変化することが分かる。また、曲線Ｂ’２において、鋳物の冷却速度がＶＣからＶＤに変化すると、共晶凝固開始温度は、Ｔ６からＴ７に変化することが分かる。
【００８５】
なお、それぞれの曲線は、特定の共晶凝固開始温度および特定の冷却速度で変曲点を示している。この変曲点は、チル臨界冷却速度と称され、チル臨界冷却速度よりも共晶凝固開始時の冷却速度が大きい場合には、鉄−セメンタイト系共晶凝固が起こり、チル臨界冷却速度よりも冷却速度が小さい場合には、鉄−黒鉛系共晶凝固が起こることが知られている。
【００８６】
ここで、本実施の形態では、凝固時の熱物性値および凝固後の熱物性値を算出するために用いる冷却速度は、ステップＳ４３３の処理で説明したように、実際の溶湯の共晶凝固開始温度ではなく、鋳物の溶湯が無限遠の時間をかけて平衡凝固したときの共晶凝固開始温度（鋳物を形成する材料が鋳鉄の場合には、鋳物溶湯中の炭素が完全黒鉛化したときの共晶凝固開始温度）を用いて算出している。しかしながら、実際の鋳造における溶湯の共晶凝固開始温度と、鋳物の溶湯中の炭素元素が完全黒鉛化した場合の共晶凝固開始温度との差は実用上無視できる範囲であるため、図１２に示す曲線を数式化することによって、数式（３）として表すことができる。
【００８７】
また、数式（４）も同様に、図２に示した冷却曲線Ｂ１〜Ｂ３の測定結果から算出している。数式（４）は、凝固潜熱量を変数とした一般的な反復解析（詳細な説明は省略する）により算出する。
【００８８】
ところで、凝固後の熱物性値は、上述のように冷却速度から算出することができるが、冷却速度からそのまま算出した熱物性値を用いるのではなく、凝固組織（例えば、パーライト／フェライトの面積比率）との関係から算出される熱物性値を用いることが好ましい。
【００８９】
このように、冷却速度から算出される凝固後の熱物性値をそのまま用いず、凝固組織との関係から求められる熱物性値を用いる理由は、次の通りである。
【００９０】
鋳物は、冷却速度によって、凝固後の金属組織が異なる。例えばダクタイルは、急冷するとチルを生じ、徐冷するとフェライトを晶出することが知られている。このように、鋳物を急冷した場合と徐冷した場合とでは凝固後の金属組織が異なるので、凝固後の熱物性値は異なることになる。
【００９１】
したがって、冷却速度から算出される値をそのまま凝固後の熱物性値として用いるよりも、上記のように冷却速度と凝固組織との関係から、凝固後の熱物性値を算出することが好ましい。この結果、凝固組織を考慮した精度の高い解析結果が得られることになる。
【００９２】
なお、鋳物を形成する材料が鋳鉄の場合、冷却速度が所定の速度より大きい場合には、鉄−セメンタイト系共晶凝固における冷却速度に基づき凝固時の熱物性値および凝固後の熱物性値を算出し、当該冷却速度が所定の速度以下の場合には、鉄−黒鉛系共晶凝固における冷却速度に基づき凝固時の熱物性値および凝固後の熱物性値を算出するようにすることが好ましい。
【００９３】
このように、所定の冷却速度を基準として熱物性値を算出する理由は、チル臨界冷却速度を挟んで、冷却速度と、凝固時および凝固後との熱物性値の関係式が変化するからである。
【００９４】
したがって、所定の速度としてチル臨界冷却速度を設定し、冷却速度がチル臨界冷却速度より大きい場合には、鉄−セメンタイト共晶凝固における冷却速度との関係式を用いて凝固時および凝固後の熱物性値を算出し、冷却速度がチル臨界冷却速度以下の場合には、鉄−黒鉛共晶凝固における冷却速度との関係式を用いて凝固時および凝固後の熱物性値を算出することが好ましい。
【００９５】
次に、セルがステップＳ４３４の処理で算出した共晶凝固開始温度に達したか否かを判定する（Ｓ４３５）。セルの温度が前記算出した共晶凝固開始温度に達していない場合、ステップＳ４３４以下の処理を繰り返し（Ｓ４３５：ＮＯ）、当該セルの温度が前記算出した共晶凝固開始温度に達した場合（Ｓ４３５：ＹＥＳ）、ステップＳ４３６の処理へ進む。
【００９６】
次に、前記共晶凝固開始温度に達したセルに対して、前記凝固時の熱物性値を設定して伝熱凝固解析を実行する（Ｓ４３６）。
【００９７】
ここで、鋳物を形成する材料が鋳鉄の場合、凝固時の熱物性値として、共晶凝固開始温度と凝固潜熱量とを用いることが好ましい。鋳物の場合、冷却速度が大きいと鉄−セメンタイト共晶凝固が起こりやすく、冷却速度が小さいと鉄−黒鉛共晶凝固が起こりやすい。鉄−セメンタイト共晶凝固と鉄−黒鉛共晶凝固では熱物性値のなかでも共晶凝固開始温度と凝固潜熱量に顕著な違いがあるからである。このように、凝固時の熱物性値のうち、共晶凝固開始温度と凝固潜熱量とをセルごとに付与することで、より解析精度が向上する。
【００９８】
次に、前記伝熱凝固解析の結果に基づき、前記セルが凝固を完了したか否かを判定する（Ｓ４３７）。セルの凝固が完了していない場合、ステップＳ４３６以下の処理を繰り返し（Ｓ４３７：ＮＯ）、当該セルの凝固が完了した場合、ステップＳ４３８の処理へ進む（Ｓ４３７：ＹＥＳ）。
【００９９】
次に、凝固が完了したセルに対して前記凝固後の熱物性値を設定し、伝熱凝固解析を実行する（Ｓ４３８）。そして、すべてのセルの凝固が完了すれば処理を終了し、ステップＳ４４０（図９参照）へ進む。
【０１００】
次に、すべての鋳造条件に対して伝熱凝固解析を実行したか否かを判定する（Ｓ４４０）。すべての鋳造条件に対して伝熱凝固解析を実行していない場合（Ｓ４４０：ＮＯ）、他の鋳造条件が選択され（Ｓ４７０）、ステップＳ４１０以下の処理を繰り返す。そして、すべての鋳造条件に対して伝熱凝固解析を実行した場合（Ｓ４４０：ＹＥＳ）、ステップＳ４５０の処理へ進む。
【０１０１】
次に、それぞれの前記鋳造条件に対応する伝熱凝固解析から得られた情報が、規格値内であるか否かを判定する（Ｓ４５０）。鋳造条件のうちの少なくとも一つの解析結果が規格値内でないと判定した場合（Ｓ４５０：ＮＯ）、解析モデルを変更する（Ｓ４６０）。そして、総ての鋳造条件に対応する伝熱凝固解析から得られた結果が、規格値内である判定した場合（Ｓ４５０：ＹＥＳ）、処理を終了する。
【０１０２】
以上のように構成された本発明に係る鋳造方案設計方法およびその設計装置によれば、複数水準の鋳造条件に対応した熱物性値が算出され、さらに、鋳物の冷却速度から算出される熱物性値を用いて伝熱凝固解析が実行される。この結果、鋳造条件の変動に伴い変化する鋳物の熱物性値に対応した解析結果が得られることになり、精度の高い鋳造方案を提供することができる。
【０１０３】
また、本発明に係る鋳造方案設計方法およびその設計装置では、伝熱凝固解析結果から得られる情報を加工することによって、解析モデルの全部または一部における欠陥判定パラメータを数値、色彩、グラフ、または図形などを用いて表示することができる。
【０１０４】
また、本発明に係る鋳造方案設計方法をコンピュータが読み取り可能なようにプログラム化したり、このプログラム化したデータを、コンピュータ読み取り可能な記録媒体に記録したりすることもできる。
【産業上の利用可能性】
【０１０５】
本発明は、鋳物の伝熱凝固解析、または、鋳物の伝熱凝固解析に基づく鋳造方案に関する技術分野に有用である。
【図面の簡単な説明】
【０１０６】
【図１】本発明に係る鋳造方案設計方法を実施するための鋳造方案設計装置の概略構成を示すブロック図である。
【図２】本発明に係る鋳造方案設計方法によって設計される鋳物を例示したものである。
【０１０７】
【図３】本発明の第１実施形態に係る鋳造方案設計装置の動作フローチャートである。
【図４】フェーディング時間による鋳物の冷却曲線の変化を示す図である。
【図５】冷却速度と鋳物の硬度との関係を示す図である。
【図６】鋳物の先端部からの距離と冷却速度との関係を示す図である。
【図７】図２に示した解析モデルを変更した解析モデルを示す図である。
【図８】鋳物の先端部からの距離と冷却速度との関係を示す図である。
【図９】本発明の第２実施形態に係る鋳造方案設計装置の動作フローチャートである。
【図１０】本発明の第２実施形態に係る鋳造方案設計装置の動作フローチャートである。
【図１１】本発明の第２実施形態に係る鋳造方案設計方法における熱物性値の算出方法を説明する図である。
【図１２】本発明の第２実施形態に係る鋳造方案設計方法における熱物性値の算出方法を説明する図である。
【符号の説明】
【０１０８】
１００    鋳造方案設計装置、
１１１    ＣＰＵ、
１１３    ＲＡＭ、
１１５    ＲＯＭ、
１１７    ハードディスク、
１２１    ディスプレイ、
１３１    入力部、
２００    鋳物
２１０    チラー
Ａ０、Ａ１５    冷却曲線、
Ｂ１〜Ｂ３    冷却曲線、
Ｃ０、Ｃ１５、Ｃ’０、Ｃ’１５    鋳物先端部からの距離と冷却速度との関係を示す曲線、
Ｂ’１〜Ｂ’３    共晶凝固開始時の冷却速度の変化曲線
Ｔ０、Ｔ１５、Ｔ１〜Ｔ７    共晶凝固開始温度、
ＶＡ〜ＶＤ    冷却速度。
【図１】
2008155230.tif 000007
【図２】
2008155230.tif 000008
【図３】
2008155230.tif 000009
【図４】
2008155230.tif 000010
【図５】
2008155230.tif 000011
【図６】
2008155230.tif 000012
【図７】
2008155230.tif 000013
【図８】
2008155230.tif 000014
【図９】
2008155230.tif 000015
【図１０】
2008155230.tif 000016
【図１１】
2008155230.tif 000017
【図１２】
2008155230.tif 000018