JP2014138950A - 鋳造シミュレーション装置及び鋳造シミュレーション方法 - Google Patents
鋳造シミュレーション装置及び鋳造シミュレーション方法 Download PDFInfo
- Publication number
- JP2014138950A JP2014138950A JP2013008639A JP2013008639A JP2014138950A JP 2014138950 A JP2014138950 A JP 2014138950A JP 2013008639 A JP2013008639 A JP 2013008639A JP 2013008639 A JP2013008639 A JP 2013008639A JP 2014138950 A JP2014138950 A JP 2014138950A
- Authority
- JP
- Japan
- Prior art keywords
- temperature
- solidus
- solidus temperature
- history
- elements
- 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.)
- Pending
Links
Images
Landscapes
- Investigating And Analyzing Materials By Characteristic Methods (AREA)
Abstract
【課題】内部凝固欠陥の発生を精度良く予測可能な鋳造シミュレーション技術を提供する。
【解決手段】実施形態によれば、鋳造シミュレーション装置10は、鋳型に充填された湯を構成する各要素について、湯が充填された開始時点から湯が凝固し、冷えて鋳物が完成するまでの温度の時間変化である温度履歴を算出する温度履歴算出部22と、算出された各要素の温度履歴に基づいて、鋳物の内部に内部凝固欠陥が生じるか否かを判定する凝固欠陥判定部24とを有している。温度履歴算出部22は、予め設定された第1の固相線温度に達する時点が遅くなるに従って、第1の固相線温度に比べて低い温度に設定された固相線温度まで、凝固前の物性値を用いて温度履歴を算出する。
【選択図】図3
【解決手段】実施形態によれば、鋳造シミュレーション装置10は、鋳型に充填された湯を構成する各要素について、湯が充填された開始時点から湯が凝固し、冷えて鋳物が完成するまでの温度の時間変化である温度履歴を算出する温度履歴算出部22と、算出された各要素の温度履歴に基づいて、鋳物の内部に内部凝固欠陥が生じるか否かを判定する凝固欠陥判定部24とを有している。温度履歴算出部22は、予め設定された第1の固相線温度に達する時点が遅くなるに従って、第1の固相線温度に比べて低い温度に設定された固相線温度まで、凝固前の物性値を用いて温度履歴を算出する。
【選択図】図3
Description
本発明の実施形態は、鋳造シミュレーション技術に関し、特に、鋳物の内部に凝固欠陥が生じるか否かを予測するためのシミュレーション技術に関する。
鋳物の鋳造を行う場合、鋳型内に充填された溶融した金属(以下、単に「湯」と記す)のうち、最後に凝固する部分には、空隙(以下、凝固欠陥と記す)が生じることがある。特に、細長い形状の鋳物に鋳造する際、湯の内部に鋳型の長手方向における温度勾配[℃/mm]が小さい部分が生じることがあり、当該部分のうち最後に凝固する部分には、凝固欠陥(以下、内部凝固欠陥と記す)が生じることがある。
上述のような鋳造の過程について、コンピュータによるシミュレーション(以下、単に「鋳造シミュレーション」と記す)を行うことが知られている。このような鋳造シミュレーションにおいては、内部凝固欠陥の発生の有無を判定するパラメータとして、経験則に基づき見出されたパラメータであるG/√R(G:温度勾配[℃/mm]、R:凝固速度[℃/sec]、以下、「凝固パラメータ」と記す)が用いられることが多い。この凝固パラメータが、内部凝固欠陥が生じるか否かを判定する所定の閾値(ニイヤマ・クライテリアと称される)以下である場合に、鋳物に内部凝固欠陥が発生すると判定される。
この凝固パラメータは、鋳型に充填された湯が、凝固し、冷えて鋳物となるまでの温度分布の時間変化を算出することにより、求めることができる。すなわち、鋳造シミュレーションにおいては、鋳型に充填された湯を多数の要素に分割し、各要素について、鋳型に湯が充填された時点から、湯が凝固し、冷えて鋳物となるまでの温度の時間変化(以下、温度履歴と記す)を算出することにより、上述の凝固パラメータを算出することができる。
ところで、上述した湯の凝固が完了する温度である「固相線温度」は、平衡状態の場合と、所定の凝固速度を以って凝固した場合で異なる。湯が冷やされて凝固する過程(以下、単に「凝固過程」と記す)における固相線温度は、その凝固速度により変化し、且つ、平衡状態の固相線温度に比べて低くなる。特に、鋳物の材料をニッケル基合金や高合金鋼(high-alloy steel)とした場合、上述のような凝固過程の固相線温度は、平衡状態の固相線温度に比べて大幅に低くなることがある。
上述のような鋳造シミュレーションにおいて、固相線温度として、平衡状態の固相線温度を用い、各要素が当該温度に達した時点から凝固前の物性値から凝固後の物性値に切り替えて各要素の温度履歴を算出した場合、当該温度履歴に基づいて凝固パラメータを算出して内部凝固欠陥の判定を行うと、実際に鋳造を行った場合には内部凝固欠陥が生じるにもかかわらず、内部凝固欠陥が生じないと判定される場合があった。
なお、特許文献1においては、流動限界固相率を用いて凝固欠陥の判定を行うシミュレーション手法が提案されている。しかし、当該手法は、引け巣等、鋳物の表面に生じる凝固欠陥の発生を予測するのに適しているものの、上述した「内部凝固欠陥」の発生の判定を行うためには、十分な実験データがないと、凝固欠陥が生じるか否かを判定する閾値を決定することができないという問題がある。
したがって、上述した内部凝固欠陥の発生を、鋳造シミュレーションにより予測する場合、従来の凝固パラメータG/√Rと、その閾値を用いつつ、より精度良く内部凝固欠陥の発生を予測する技術が要望されている。
そこで、本発明の実施形態は、内部凝固欠陥の発生を精度良く予測可能な鋳造シミュレーション技術を提供することを目的とする。
上述の目的を達成するため、本発明の実施形態の鋳造シミュレーション装置は、鋳型に充填された湯を構成する各要素について、湯が充填された開始時点から湯が凝固し、冷えて鋳物が完成するまでの温度の時間変化である温度履歴を算出する温度履歴算出部と、算出された各要素の温度履歴に基づいて、鋳物の内部に内部凝固欠陥が生じるか否かを判定する凝固欠陥判定部と、を有し、前記温度履歴算出部は、予め設定された第1の固相線温度に達する時点が遅くなるに従って、第1の固相線温度に比べて低い温度に設定された固相線温度まで、凝固前の物性値を用いて温度履歴を算出することを特徴とする。
また、本発明の実施形態の鋳造シミュレーション方法は、鋳型に充填された湯を構成する各要素について、湯が充填された開始時点から湯が凝固し、冷えて鋳物が完成するまでの温度の時間変化である温度履歴を温度履歴算出部が算出するステップと、算出された温度履歴に基づいて、鋳物の内部に凝固欠陥が生じるか否かを凝固欠陥判定部が判定するステップと、を有し、前記温度履歴を温度履歴算出部が算出するステップは、予め設定された第1の固相線温度に達する時点が遅くなるに従って、第1の固相線温度に比べて低い温度に設定された固相線温度まで、凝固前の物性値を用いて温度履歴を算出することを特徴とする。
本発明の実施形態によれば、湯が液相から固相になる過程において多様な析出物が生じて、残りの湯の成分の融点が下がっても、第1の固相線温度に達する時点が遅いものについては、より低い温度に設定された第2の固相線温度まで、凝固前の物性値を用いて算出するため、より高い精度で各要素の温度履歴を算出し、この温度履歴に基づいて内部凝固欠陥の発生を、より精度良く予測することができる。
以下に、本発明の実施形態について図面を参照して説明する。なお、以下に説明する実施形態により、本発明が限定されるものではなく、その要旨を逸脱しない範囲において種々の変更が可能である。
〔第1の実施形態〕
まず、本実施形態の鋳造シミュレーション技術が適用される鋳造の一例について図1及び図2を用いて説明する。図1は、鋳物の形状と、鋳物の内部に生じる内部凝固欠陥の一例を示す説明図である。図2は、内部凝固欠陥が生じる範囲における温度勾配の一例を示すグラフである。
まず、本実施形態の鋳造シミュレーション技術が適用される鋳造の一例について図1及び図2を用いて説明する。図1は、鋳物の形状と、鋳物の内部に生じる内部凝固欠陥の一例を示す説明図である。図2は、内部凝固欠陥が生じる範囲における温度勾配の一例を示すグラフである。
図1及び図2の例において、鋳型1は、軸心が水平方向に延びている略円筒状をなしており、その軸心を図1に一点鎖線Aで示す。鋳型1には、例えば、砂型が用いられる。鋳型1には、溶融した金属である「湯」が、鋳型1の軸方向の一方側の上部に設けられた湯溜り(又は押湯)3から供給される。湯溜り(又は押湯)3から供給された湯は、図に矢印Bで示すように、鋳型1内を軸方向に流れて、鋳型1内に充填される。鋳型1に充填された湯は、時間が経過するに従って温度が低下することにより、液相から徐々に固相に変化する。湯は、ある温度に達すると完全に凝固する。この「完全に凝固する温度」が、いわゆる「固相線温度」であり、図2に一点鎖線Cで示す。固相線温度は、主に、湯を構成する金属の組成と、湯の冷却速度に応じて定まる。
鋳型1内における湯の温度分布の一例を、図2に直線Dで示す。図1及び図2に示すように、細長い形状の鋳型1の場合、湯の内部には、鋳型1の軸方向(長手方向)においては、固相線温度Cの近傍において温度勾配[℃/mm]が小さい領域(図1に斜線で示す)が生じることがある。当該領域の温度勾配を図2に二点鎖線Eで示す。この領域のうち、最後に凝固する部分には、図1に示すように、内部凝固欠陥5が生じることがある。このような内部凝固欠陥5の発生を予測するための、本実施形態の鋳造シミュレーション装置及び鋳造シミュレーション方法の一例について以下に説明する。
次に、本実施形態の鋳造シミュレーション装置について説明する。図3は、本実施形態の鋳造シミュレーション装置のハードウェアの構成を模式的に示すブロック図である。
図3に示すように、本実施形態において、鋳造シミュレーション装置10は、一般的なコンピュータとして構成されており、各種の演算処理を行う中央演算処理装置(以下、演算部と記す)20と、各種の情報を記憶する補助記憶装置(以下、記憶部と記す)30とを有している。なお、記憶部30には、ハードディスク及びEPROMが用いられている。本実施形態の鋳造シミュレーション方法において用いられるプログラム及びデータベースは、記憶部30に予め記憶されている。一方、演算部20においては、各種のプログラムの演算処理が行われる。なお、演算部20には、主記憶装置(RAM)が含まれている。加えて、鋳造シミュレーション装置10には、操作者が各種の情報を入力する入力装置(以下、入力部と記す)40と、演算部20により演算処理された各種の情報を出力して表示する表示装置(以下、表示部と記す)50とを有している。鋳造シミュレーション装置10において、演算部20、記憶部30、入力部40、及び表示部50は、情報通信線12を介して接続されており、相互に各種の情報を送受信可能に構成されている。なお、演算部20における演算処理、及び記憶部30が記憶する情報の詳細については後述する。
本実施形態において、演算部20は、鋳型1に充填された湯を構成する各要素について、温度の時間変化(以下、温度履歴と記す)を算出する機能(以下に「温度履歴算出部」と記す)22と、算出された各要素の温度履歴に基づいて、上述した内部凝固欠陥が生じるか否かを判定する機能(以下、「凝固欠陥判定部」と記す)24とを有している。以下に、鋳造シミュレーション装置10が実行する情報処理、すなわち本実施形態の鋳造シミュレーション方法について、図1〜図5を用いて説明する。図4は、本実施形態の鋳造シミュレーション方法の全体的な情報処理を示すフローチャートである。図5は、本実施形態の鋳造シミュレーション方法のうち、図4にステップS06で示す凝固シミュレーションを示すサブルーチンである。
なお、図5に示す凝固シミュレーションのルーチンは、鋳造シミュレーション装置10の演算部20において温度履歴を算出する温度履歴算出部22により実行される。なお、この凝固シミュレーションには、一例として有限要素法が用いられる。上述した湯は、予め多数の要素に分割されている。各要素に関する情報は、予め記憶部30に記憶されている。
なお、上述した実施形態においては、シミュレーションには、有限要素法を用いるものとしたが、シミュレーションの態様は、これに限定されるものではない。シミュレーションに差分法を用いることもでき、この場合、上述した「要素」は、「セル」と呼称される。
図3及び図4に示すように、鋳造シミュレーション装置10の演算部20は、ステップS02において、鋳造シミュレーションの条件に関する各種の情報を取得する。当該情報には、鋳型(鋳物)や押湯の形状及び寸法のデータ、鋳物の材料の物性値に関するデータ、鋳型の温度、鋳型と鋳物(又は湯)との間との熱伝達率等、境界条件に関するデータなどが含まれている。これらの情報は、入力部40から入力され、データベースとして記憶部30に予め記憶されている。
なお、鋳物の材料の物性値に関するデータには、湯の凝固が完了する前に用いられる各種の物性値(以下、単に「凝固前の物性値」と記す)と、湯の凝固が完了した後に用いられる各種の物性値(以下、単に「凝固後の物性値」と記す)が含まれている。また、鋳物の材料の物性値に関するデータには、鋳物の材料に応じて設定された複数の固相線温度が含まれている。固相線温度の設定手法については、後述する。
そして、ステップS04において、演算部20は「湯流れシミュレーション」を実行する。「湯流れシミュレーション」は、鋳型に湯が注ぎ込まれてから、鋳型に湯が充填されるまでの湯の流動と温度分布を解析するものである。これにより、演算部20は、鋳型に湯が充填された状態、すなわち後述する「凝固シミュレーション」の開始時点における湯の温度分布を決定する。
そして、ステップS06において、演算部20は、図5に示す「凝固シミュレーション」のルーチンを実行する。「凝固シミュレーション」においては、鋳型に湯が充填された時点(以下、開始時点と記す)から、液相にあった湯が、凝固して固相となり、さらに冷えて鋳物が完成するまでの温度履歴を算出する。なお、演算部20のうち、当該温度履歴を算出する部分を、以下の説明において「温度履歴算出部」と記し、図3に符号22で示す。
図5に示すように、ステップS102において、算出部は、各種の情報を取得する。この情報には、上述した凝固前の物性値、凝固後の物性値、第1の固相線温度、第2の固相線温度が含まれている。なお、凝固前の物性値、及び凝固後の物性値は、予め実験等により求められており、予め記憶部30に記憶されている。
ここで、第2の固相線温度の設定手法について、図6を用いて説明する。図6は、ある金属材料の固相率に対する液相線温度及び固相線温度を示すグラフである。なお、図6において、平衡状態を実線で示し、凝固過程の一例を破線で示している。
図6に示すように、固相率が0%であるときの温度すなわち液相線温度TLは、液相にある金属において凝固が開始される温度であり、平衡状態(実線)と凝固過程(破線)でほぼ同一である。これに対して、固相率が100%であるときの温度すなわち固相線温度は、平衡状態(実線)と凝固過程(破線)が大きく異なり、凝固過程の固相線温度TS2は、平衡状態の固相線温度TS1に比べて低い値となっている。
なお、鋳物の材料が低合金鋼(low-alloy steel)である場合、平衡状態の固相線温度TS1と凝固過程の固相線温度TS2との差異は小さい。しかし、材料がNi基合金や高合金鋼(high-alloy steel)である場合は、平衡状態の固相線温度TS1と凝固過程の固相線温度TS2との差異は、図6に示すように大きくなる。これは、液相から固相になる過程において多様な析出物が生じ、残りの湯の成分の融点が下がるためである。
このような金属で構成される湯を、ある程度の凝固速度[℃/sec]を以って凝固させる場合(破線参照)、図に点Fで示すように、温度が、平衡状態の固相線温度TS1に達しても、完全に固相とならない、すなわち完全に凝固していない。当該材料が完全に凝固するのは、平衡状態の固相線温度TS1に比べて低い値である凝固過程の固相線温度TS2に達した時点となる。
そこで、本実施形態の鋳造シミュレーション装置10においては、第1の固相線温度として、平衡状態の固相線温度TS1が予め設定されている。平衡状態の固相線温度は、予め実験により求められており、記憶部30に予め記憶されている。加えて、第2の固相線温度として、第1の固相線温度に比べて低い温度である凝固過程の固相線温度TS2が予め設定されている。なお、凝固過程の固相線温度は、平衡状態の固相線温度と、凝固速度[℃/sec]、及び材料の組成に基づいて、熱力学平衡計算を行うことにより求めることが可能である。
図5に戻り、ステップS104において、各要素について、温度が第1の固相線温度に達するまで、凝固前の物性値を用いて温度履歴の計算を、温度履歴算出部22が行う。この計算は、所定の時間ステップごとに、各要素について並行して行われる。このとき、他の要素との熱伝導等を計算する場合は、凝固前の物性値が用いられる。
そして、ステップS106において、全ての要素(以下、単に「全要素」と記す)のうち、第1の固相線温度に達した順番が、閾値として予め設定されている「所定の順番」に比べて遅いものであるか否かを判定する。すなわち、全要素のうち温度が第1の固相線温度に達している要素の割合が、所定の閾値に比べて高いか否かを判定する。なお、当該「所定の順番」は、予め演算部20に設定されており、記憶部30に予め記憶されている。
第1の固相線温度に達した順番が、所定の順番に比べて早い場合(ステップS106、No)、すなわち、当該要素が第1の固相線温度に達した時点は比較的早いものであり、当該時点において、全要素のうち既に温度が第1の固相線温度に達している要素の割合が低いと判定した場合には、当該要素は、完全に凝固しているものとして、ステップS110に進む。
ステップS110において、当該要素については、凝固前の物性値から凝固後の物性値に切り替える。凝固後の物性値を用いて、他の要素との間の熱伝導解析を行なうことにより、温度履歴算出部22は当該要素の温度履歴を算出する。
一方、ステップS106において、第1の固相線温度に達した順番が、所定の順番に比べて遅い場合(ステップS106、Yes)、すなわち、当該要素が第1の固相線温度に達した時点は比較的遅いものであり、当該時点において、全要素のうち既に温度が第1の固相線温度に達している要素の割合が高いと判定した場合には、当該要素は、未だ完全に凝固していないものとして、ステップS108に進む。
ステップS108において、当該要素については、未だ完全に凝固していないため、温度が第2の固相線温度に達するまで、凝固前の物性値を用いて温度履歴を算出する。当該要素の温度が、第2の固相線温度に達すると、当該要素も完全に凝固したものと判定する。そして、ステップS110において、凝固後の物性値を用いて熱伝導解析を行うことにより当該要素が冷えて鋳物が完成するまでの温度履歴を、温度履歴算出部22が算出する。
以上、ステップS102〜S110に示す計算(凝固シミュレーション)を、各要素について行うことにより、湯が充填された開始時点から、湯が凝固し、冷えて鋳物が完成するまでの温度履歴を算出する。
上述した凝固シミュレーションが完了すると、図4に示すステップS08において、凝固パラメータである、G/√Rの計算を行う。なお、Gは、温度勾配[℃/mm]であり、Rは、凝固速度[℃/sec]である。温度勾配G及び凝固速度Rは、上述した凝固シミュレーションにおいて算出された各要素の温度履歴から求めることができる。
そして、ステップS10において、鋳物の内部に「内部凝固欠陥」が生じるか否かを各要素について判定する。具体的には、凝固パラメータ:G/√Rを、予め設定されている閾値(例えば、ニイヤマ・クライテリア)と比較することにより、内部凝固欠陥が生じるか否かを判定することができる。
以上に説明したように本実施形態の鋳造シミュレーション方法においては、鋳型に充填された湯を構成する各要素について、湯が充填された開始時点から湯が凝固し、冷えて鋳物が完成するまでの温度の時間変化である温度履歴を温度履歴算出部22が算出するステップS06(S102〜S110)と、算出された各要素の温度履歴に基づいて、鋳物の内部に内部凝固欠陥が生じるか否かを凝固欠陥判定部24が判定するステップS08,S10とを有している。温度履歴を温度履歴算出部22が算出するステップS06(S102〜S110)は、予め設定された第1の固相線温度に達する時点が遅くなるに従って、第1の固相線温度に比べて低い温度に設定された固相線温度(第2の固相線温度)まで、凝固前の物性値を用いて温度履歴を算出するものとした。これにより、湯が液相から固相になる過程において多様な析出物が生じて、残りの湯の成分の融点が下がっても、第1の固相線温度に達する時点が遅いものについては、より低い温度に設定された第2の固相線温度まで、凝固前の物性値を用いて計算するため、より高い精度で各要素の温度履歴を算出し、この温度履歴に基づいて内部凝固欠陥の発生を、より精度良く予測することができる。特に、鋳物の材料が、Ni基合金である場合や高合金鋼である場合に、本実施形態の鋳造シミュレーション方法は、内部凝固欠陥の発生を精度良く予測することができる。
また、本実施形態の鋳造シミュレーション方法において、温度履歴を温度履歴算出部22が算出するステップS06(S102〜S110)は、全要素のうち第1の固相線温度に達する時点が早い要素については、第1の固相線温度に達するまで凝固前の物性値を用いて温度履歴を算出し(ステップS104)、全要素のうち第1の固相線温度に達する時点が遅い要素については、第1の固相線温度に比べて低い温度に設定された第2の固相線温度に達するまで凝固前の物性値を用いて温度履歴を計算する(ステップS108)ものとした。第1の固相線温度に達する時点が遅い要素については、凝固前の物性値を用いて計算する温度の下限値を、第1の固相線温度から第2の固相線温度に変更するだけで、より精度良く各要素の温度履歴を算出することができる。
また、本実施形態の鋳造シミュレーション方法において、温度履歴を温度履歴算出部22が算出するステップS06(S102〜S110)は、温度が第1の固相線温度に達するまでは凝固前の物性値を用いて温度履歴を算出し(ステップS104)、温度が第1の固相線温度に達したときに、全要素のうち温度が第1の固相線温度に達している要素の割合が所定の閾値に比べて高い場合には、第2の固相線温度に達するまで凝固前の物性値を用いて温度履歴を算出する(ステップS108)ものとした。第1の固相線温度に達している要素の割合が高い場合、すなわち実際には、析出物が既に析出して、融点が低下していると推定される要素については、より低い温度に設定された第2の固相線温度に達するまで凝固前の物性値を用いて温度履歴を算出するため、より高い精度で各要素の温度履歴を算出することができる。
なお、本実施形態の鋳造シミュレーション方法において、温度履歴を温度履歴算出部22が算出するステップS06(S102〜S110)は、温度が第1の固相線温度に達したとき全要素のうち第1の固相線温度に達した順番が所定の順番に比べて遅い(ステップS106、Yes)場合に、全要素のうち温度が第1の固相線温度に達している要素の割合が所定の閾値に比べて高いと判定するものとした。全要素のうち第1の固相線温度に達している要素の割合が高いか否かを、簡素な計算手法で実現することができる。
なお、上述した実施形態においては、上述した凝固シミュレーションにおいて、各要素について、温度履歴を計算している最中に、第1の固相線温度から、より低い温度である第2の固相線温度に変更するものとしたが、第1の固相線温度に達する時点が遅くなるに従って、第1の固相線温度に比べて低い温度に設定する態様は、これに限定されるものではない。以下に、他の態様について説明する。
〔第2の実施形態〕
本実施形態の鋳造シミュレーション方法について図3、図4、及び図7を用いて説明する。図7は、図4に示す鋳造シミュレーション方法のうち、本実施形態の凝固シミュレーションを示すサブルーチンである。本実施形態の鋳造シミュレーション方法は、凝固シミュレーションのみが、上述した第1の実施形態と異なっており、以下に詳細を説明する。なお、第1の実施形態と略共通の構成については、同一の符号を付して説明を省略する。
本実施形態の鋳造シミュレーション方法について図3、図4、及び図7を用いて説明する。図7は、図4に示す鋳造シミュレーション方法のうち、本実施形態の凝固シミュレーションを示すサブルーチンである。本実施形態の鋳造シミュレーション方法は、凝固シミュレーションのみが、上述した第1の実施形態と異なっており、以下に詳細を説明する。なお、第1の実施形態と略共通の構成については、同一の符号を付して説明を省略する。
図3及び図7に示すように、ステップS202において、温度履歴算出部22は、各種の情報を取得する。この情報には、上述した凝固前の物性値、第1の固相線温度が含まれている。
そして、ステップS204において、各要素について、温度が第1の固相線温度に達するまで、凝固前の物性値を用いて温度履歴の算出を、温度履歴算出部22が行う。この計算は、所定の時間ステップごとに、各要素について並行して行われる。
そして、ステップS206において、第1の固相線温度に到達した順番に応じて、要素を分別する。具体的には、第1の固相線温度に達する時点が早い要素と、遅い要素に、温度履歴算出部22が分別する。なお、この分別は、全要素のうち、第1の固相線温度に達した順番が、所定の順番(閾値)に比べて遅いものであるか否かにより、決定することができる。
そして、ステップS208において、物性値の再設定を行う。全要素のうち第1の固相線温度に達する時点が早いと分別された要素については、所定の固相線温度として、そのまま第1の固相線温度に設定される。一方、第1の固相線温度に達する時点が遅いと分別された要素については、所定の固相線温度として、第1の固相線温度に比べて低い温度である第2の固相線温度に設定される。つまり、第1の固相線温度に達する時点が遅い要素については、第1の固相線温度から第2の固相線温度に変更される。なお、これら固相線温度以外の「凝固前の物性値」も、第1の固相線温度に達する時点が早い要素と、遅い要素で異ならせることが好ましい。
そして、ステップS210において、ステップS208で設定された固相線温度まで温度履歴の算出を行う。すなわち、第1の固相線温度に達する時点が早いと分別された要素については、上述した開始時点から第1の固相線温度に達するまで、凝固前の物性値を用いて温度履歴を算出する。これと共に、第1の固相線温度に達する時点が遅いと分別された要素については、開始時点から第1の固相線温度に比べて低い温度である第2の固相線温度に達するまで、凝固前の物性値を用いて温度履歴を算出する。
そして、ステップS212において、各要素について設定された所定の固相線温度に温度が到達すると、当該要素については、凝固前の物性値から凝固後の物性値に切り替える。第1の固相線温度に達するのが早い要素については、第1の固相線温度に達した時点で、凝固前の物性値から凝固後の物性値に切り替える。一方、第1の固相線温度に達するのが遅い要素については、第2の固相線温度に達した時点で、凝固前の物性値から凝固後の物性値に切り替える。設定された固相線温度に達した後は、凝固後の物性値を用いて、他の要素との間の熱伝導解析を行なうことにより、温度履歴算出部22は当該要素が冷えて鋳物が完成するまでの温度履歴を算出する。
以上、ステップS202〜S212に示す計算(凝固シミュレーション)を、各要素について行うことにより、湯が充填された開始時点から、湯が凝固し、冷えて鋳物が完成するまでの温度履歴を算出する。凝固シミュレーションが完了すると、図4に示すステップS08において、凝固パラメータである、G/√Rの計算を行う。そして、ステップS10において、鋳物の内部に「内部凝固欠陥」が生じるか否かを各要素について判定する。凝固パラメータ:G/√Rを、予め設定されている閾値と比較することにより、内部凝固欠陥が生じるか否かを判定することができる。
以上に説明したように本実施形態の鋳造シミュレーション方法においては、温度履歴を温度履歴算出部22が算出するステップ(S202〜S212)は、全要素について、温度が第1の固相線温度に達するまで、凝固前の物性値を用いて温度履歴を算出し(S204)、全要素を、第1の固相線温度に達する時点が早い要素と、遅い要素とに分別する(S206)。そして、物性値を再設定し(S208)、凝固シミュレーションを開始時点からやり直す。第1の固相線温度に達する時点が早いと分別された要素については、前記開始時点から第1の固相線温度まで、凝固前の物性値を用いて温度履歴を算出すると共に、第1の固相線温度に達する時点が遅いと分別された要素については、前記開始時点から第2の固相線温度に達するまで凝固前の物性値を用いて温度履歴を算出する(S210)。湯が充填された開始時点から湯が凝固し、冷えて鋳物が完成するまでの凝固シミュレーションの途中において、要素に応じて固相線温度を変更する必要がなくなる。各要素の温度履歴の算出を、より簡素な計算手法で実現することができる。
なお、上述した実施形態において、全要素のうち第1の固相線温度に達する時点が遅い要素については、第1の固相線温度に比べて低い温度に設定された第2の固相線温度に達するまで凝固前の物性値を用いて温度履歴を計算するものとした。しかし、本発明に係る第1の固相線温度に比べて低い温度に設定された「所定の固相線温度」の態様は、この態様に限定されるものではない。例えば、上述した第1及び第2の固相線温度に加え、各要素の温度が第1の固相線温度に達する時点が遅くなるに従って、さらに低い温度に設定された第3の固相線温度まで、凝固前の物性値を用いて温度履歴を算出するものとしても良い。また、各要素の温度が第1の固相線温度に達する時点が遅くなるに従って、設定される固相線温度を徐々に低下させるものとしても良い。
〔その他の実施形態〕
なお、上述した実施形態の鋳造シミュレーション装置10(図3参照)は、本発明を分かりやすく説明するために詳細に記載したものであり、必ずしも説明した全ての構成を備えるものに限定されるものではない。また、上述した実施形態の鋳造シミュレーション装置10の構成の一部を他の構成に置き換えることも可能である。また、各実施形態の構成の一部について、他の構成の追加、削除、置換をすることが可能である。
なお、上述した実施形態の鋳造シミュレーション装置10(図3参照)は、本発明を分かりやすく説明するために詳細に記載したものであり、必ずしも説明した全ての構成を備えるものに限定されるものではない。また、上述した実施形態の鋳造シミュレーション装置10の構成の一部を他の構成に置き換えることも可能である。また、各実施形態の構成の一部について、他の構成の追加、削除、置換をすることが可能である。
また、上述した鋳造シミュレーション装置10の各構成、機能、例えば演算部20、温度履歴算出部22、凝固欠陥判定部24、記憶部30等は、それらの一部又は全部を、例えば、集積回路で設計する等によりハードウェアで実現しても良い。また、上記の各構成、機能等は、プロセッサがそれぞれの機能を実現するプログラムを解釈し、実行することにより、ソフトウェアで実現しても良い。各機能を実現するプログラム、テーブル、ファイル等の情報は、メモリやハードディスク等の記憶装置、または、SDカードやDVD等の記録媒体におくことができる。
また、上述した各実施形態においては、鋳造シミュレーションの計算手法として、有限要素法を用いた場合について説明したが、本発明に係る鋳造シミュレーションの計算手法として差分法を用いるものとしても良い。なお、差分法を用いた場合、上述した「要素」については、「セル」と呼称される。
以上、本発明のいくつかの実施形態について説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これらの実施形態はその他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これらの実施形態やその変形例は、発明の要旨に含まれると同様に、特許請求の範囲に記載された発明と均等の範囲に含まれる。
1 鋳型
10 鋳造シミュレーション装置
20 演算部
22 温度履歴算出部
24 凝固欠陥判定部
30 記憶部
40 入力部
50 表示部
10 鋳造シミュレーション装置
20 演算部
22 温度履歴算出部
24 凝固欠陥判定部
30 記憶部
40 入力部
50 表示部
Claims (12)
- 鋳型に充填された湯を構成する各要素について、湯が充填された開始時点から湯が凝固し、冷えて鋳物が完成するまでの温度の時間変化である温度履歴を算出する温度履歴算出部と、
算出された各要素の温度履歴に基づいて、鋳物の内部に内部凝固欠陥が生じるか否かを判定する凝固欠陥判定部と、
を有し、
前記温度履歴算出部は、予め設定された第1の固相線温度に達する時点が遅くなるに従って、第1の固相線温度に比べて低い温度に設定された固相線温度まで、凝固前の物性値を用いて温度履歴を算出する
ことを特徴とする鋳造シミュレーション装置。 - 前記温度履歴算出部は、
全要素のうち第1の固相線温度に達する時点が早い要素については、第1の固相線温度に達するまで凝固前の物性値を用いて温度履歴を算出し、
全要素のうち第1の固相線温度に達する時点が遅い要素については、第1の固相線温度に比べて低い温度に設定された第2の固相線温度に達するまで凝固前の物性値を用いて温度履歴を算出する
ことを特徴とする請求項1に記載の鋳造シミュレーション装置。 - 前記温度履歴算出部は、
温度が第1の固相線温度に達するまでは凝固前の物性値を用いて温度履歴を算出し、
温度が第1の固相線温度に達したとき、全要素のうち温度が第1の固相線温度に達している要素の割合が所定の閾値に比べて高い場合には、第2の固相線温度に達するまで凝固前の物性値を用いて温度履歴を算出する
ことを特徴とする請求項2に記載の鋳造シミュレーション装置。 - 前記温度履歴算出部は、温度が第1の固相線温度に達したとき全要素のうち第1の固相線温度に達した順番が所定の順番に比べて遅い場合に、全要素のうち温度が第1の固相線温度に達している要素の割合が所定の閾値に比べて高いと判定する
ことを特徴とする請求項3に記載の鋳造シミュレーション装置。 - 前記温度履歴算出部は、
全要素について、温度が第1の固相線温度に達するまで、凝固前の物性値を用いて温度履歴を算出し、
全要素を、第1の固相線温度に達する時点が早い要素と、遅い要素とに分別し、
第1の固相線温度に達する時点が早いと分別された要素については、前記開始時点から第1の固相線温度まで、凝固前の物性値を用いて温度履歴を算出すると共に、
第1の固相線温度に達する時点が遅いと分別された要素については、前記開始時点から第2の固相線温度に達するまで凝固前の物性値を用いて温度履歴を算出する
ことを特徴とする請求項2に記載の鋳造シミュレーション装置。 - 前記第1の固相線温度は、平衡状態の固相線温度である、ことを特徴とする請求項1ないし請求項5に記載の鋳造シミュレーション装置。
- 湯を構成する材料に対する、平衡状態の固相線温度および第2の固相線温度が、予め記憶されている記憶部を、さらに備えることを特徴とする請求項1ないし請求項6に記載の鋳造シミュレーション装置。
- 鋳型に充填された湯を構成する各要素について、湯が充填された開始時点から湯が凝固し、冷えて鋳物が完成するまでの温度の時間変化である温度履歴を温度履歴算出部が算出するステップと、
算出された温度履歴に基づいて、鋳物の内部に凝固欠陥が生じるか否かを凝固欠陥判定部が判定するステップと、
を有し、
前記温度履歴を温度履歴算出部が算出するステップは、
予め設定された第1の固相線温度に達する時点が遅くなるに従って、第1の固相線温度に比べて低い温度に設定された固相線温度まで、凝固前の物性値を用いて温度履歴を算出する
ことを特徴とする鋳造シミュレーション方法。 - 前記温度履歴を温度履歴算出部が算出するステップは、
全要素のうち第1の固相線温度に達する時点が早い要素については、第1の固相線温度に達するまで凝固前の物性値を用いて温度履歴を算出し、
全要素のうち第1の固相線温度に達する時点が遅い要素については、第1の固相線温度に比べて低い温度に設定された第2の固相線温度に達するまで凝固前の物性値を用いて温度履歴を算出する
ことを特徴とする請求項8に記載の鋳造シミュレーション方法。 - 前記温度履歴を温度履歴算出部が算出するステップは、
温度が第1の固相線温度に達するまでは凝固前の物性値を用いて温度履歴を算出し、
温度が第1の固相線温度に達したときに、全要素のうち温度が第1の固相線温度に達している要素の割合が所定の閾値に比べて高い場合には、第2の固相線温度に達するまで凝固前の物性値を用いて温度履歴を算出する
ことを特徴とする請求項9に記載の鋳造シミュレーション方法。 - 前記温度履歴を温度履歴算出部が算出するステップは、
温度が第1の固相線温度に達したとき全要素のうち第1の固相線温度に達した順番が所定の順番に比べて遅い場合に、全要素のうち温度が第1の固相線温度に達している要素の割合が所定の閾値に比べて高いと判定する
ことを特徴とする請求項10に記載の鋳造シミュレーション方法。 - 前記温度履歴を温度履歴算出部が算出するステップは、
全要素について、温度が第1の固相線温度に達するまで、凝固前の物性値を用いて温度履歴を算出し、
全要素を、第1の固相線温度に達する時点が早い要素と、遅い要素とに分別し、
第1の固相線温度に達する時点が早いと分別された要素については、前記開始時点から第1の固相線温度まで、凝固前の物性値を用いて温度履歴を算出すると共に、
第1の固相線温度に達する時点が遅いと分別された要素については、前記開始時点から第2の固相線温度に達するまで凝固前の物性値を用いて温度履歴を算出する
ことを特徴とする請求項9に記載の鋳造シミュレーション方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2013008639A JP2014138950A (ja) | 2013-01-21 | 2013-01-21 | 鋳造シミュレーション装置及び鋳造シミュレーション方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2013008639A JP2014138950A (ja) | 2013-01-21 | 2013-01-21 | 鋳造シミュレーション装置及び鋳造シミュレーション方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2014138950A true JP2014138950A (ja) | 2014-07-31 |
Family
ID=51416007
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2013008639A Pending JP2014138950A (ja) | 2013-01-21 | 2013-01-21 | 鋳造シミュレーション装置及び鋳造シミュレーション方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2014138950A (ja) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111940680A (zh) * | 2020-07-14 | 2020-11-17 | 陕西柴油机重工有限公司 | 一种中高速大功率柴油机球墨铸铁飞轮的成型方法 |
CN117034176A (zh) * | 2023-10-08 | 2023-11-10 | 西安石油大学 | 一种定向凝固过程铸件异常预测系统 |
-
2013
- 2013-01-21 JP JP2013008639A patent/JP2014138950A/ja active Pending
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111940680A (zh) * | 2020-07-14 | 2020-11-17 | 陕西柴油机重工有限公司 | 一种中高速大功率柴油机球墨铸铁飞轮的成型方法 |
CN111940680B (zh) * | 2020-07-14 | 2024-05-24 | 陕西柴油机重工有限公司 | 一种中高速大功率柴油机球墨铸铁飞轮的成型方法 |
CN117034176A (zh) * | 2023-10-08 | 2023-11-10 | 西安石油大学 | 一种定向凝固过程铸件异常预测系统 |
CN117034176B (zh) * | 2023-10-08 | 2024-01-09 | 西安石油大学 | 一种定向凝固过程铸件异常预测系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Kavoosi et al. | Influence of cooling rate on the solidification behavior and microstructure of IN738LC superalloy | |
Monroe et al. | Prediction of hot tearing using a dimensionless Niyama criterion | |
Bazhenov et al. | Simulation of fluidity and misrun prediction for the casting of 356.0 aluminum alloy into sand molds | |
Reilly et al. | Process modeling of low-pressure die casting of aluminum alloy automotive wheels | |
JP2014138950A (ja) | 鋳造シミュレーション装置及び鋳造シミュレーション方法 | |
JP6665849B2 (ja) | 鋳造品の機械的特性予測方法、機械的特性予測システム、機械的特性予測プログラム及び機械的特性予測プログラムを記録したコンピュータ読取り可能な記録媒体 | |
KR101312474B1 (ko) | 용융 합금의 응고 해석 방법 및 그 응고 해석 프로그램을 저장한 컴퓨터-판독가능 저장 매체 | |
KR100986923B1 (ko) | 응고 해석 방법 및 장치 | |
Wang et al. | A submerged-gate casting method | |
JP5351575B2 (ja) | 凝固割れ予測方法、これを用いた鋳造方法、凝固割れ予測装置、及び凝固割れ予測プログラム | |
JP4501845B2 (ja) | 鋳造解析装置及び鋳造解析方法 | |
Li et al. | Influence of mold temperature on microstructure and shrinkage porosity of the A357 alloys in gravity die casting | |
JP2007125589A (ja) | 鋳造解析装置及び鋳造解析方法 | |
WO2009133602A1 (ja) | 鋳鉄鋳物の引け巣の予測および防止方法 | |
JP2009262235A (ja) | 凝固割れ予測方法、これを用いた鋳造方法、凝固割れ予測装置、及び凝固割れ予測プログラム | |
Grilo et al. | Thermal analysis strategy for Al-7Si-Mg alloy sand-casting modelling | |
JP2020168636A (ja) | 鋳造解析方法 | |
Manjabacas et al. | A comparison between traditional criteria and FEM analysis results for gravity casting feeding and risering systems | |
Fardi-Ilkhchy et al. | Effect of Pressure on Microstructure and Cooling Curves of A356 Aluminum Alloy During Solidification | |
Pulisheru et al. | Prediction of Hot Spot and Hot Tear of the Al–Cu Cast Alloy by Casting Simulation Software | |
Mahesh et al. | Simulation of sand casting to predict porosity for Ductile Iron | |
JP2008155248A (ja) | 鋳物の伝熱凝固解析方法 | |
JP2005111544A (ja) | 鋳物の凝固解析方法、凝固解析プログラム、記録媒体、マイクロポロシティの発生予測方法、マイクロポロシティ発生予測プログラム及び記録媒体 | |
JP2006159221A (ja) | マイクロポロシティの発生予測方法、マイクロポロシティ発生予測プログラム及び記録媒体 | |
KR20230081466A (ko) | 주조품 미세조직 특성값 예측 시스템 및 방법 |