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