CN100458801C - 一种铸件凝固模拟中预测缩松的方法 - Google Patents
一种铸件凝固模拟中预测缩松的方法 Download PDFInfo
- Publication number
- CN100458801C CN100458801C CNB2006101662298A CN200610166229A CN100458801C CN 100458801 C CN100458801 C CN 100458801C CN B2006101662298 A CNB2006101662298 A CN B2006101662298A CN 200610166229 A CN200610166229 A CN 200610166229A CN 100458801 C CN100458801 C CN 100458801C
- Authority
- CN
- China
- Prior art keywords
- active cell
- unit
- adjacent
- calculating
- prediction
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
Description
一、技术领域
本发明涉及一种铸件凝固模拟中预测缩松的方法,属于计算机应用技术领域。
二、技术背景
铸件凝固模拟软件主要是通过模拟铸件的凝固过程来预测铸造缺陷的形成,所预测的铸造缺陷主要是缩孔和缩松,预测缩松是铸件凝固模拟软件的一项重要的任务。
目前主要应用的预测缩松的判据是 其中G是温度梯度(K/cm),R是冷却速度(K/min),C是临界值。
所有这些仅仅是考虑了凝固条件对缩松形成的影响,而与缩松形成有着直接影响的几个因素却没有考虑。如与补缩条件有关的金属液的压力P,补缩通道上影响压力损失的固相率梯度Gs等。
三、发明内容
为克服现有技术的不足,本发明提供一种铸件凝固模拟中预测缩松的方法。
本发明方法的预测缩松判据公式为: 根据该公式对缩松进行预测,并取得了良好的效果。
一种铸件凝固模拟中预测缩松的方法,根据直接影响缩松形成的各个因素,将铸件凝固模拟软件运行于计算机上,在预测缩松判据公式 中,P为金属液的压力,Gs固相率梯度,R是冷却速度(K/min),C是临界值,具体确定步骤为:
1)计算程序开始,利用通用的差分方法计算出铸件的温度场,根据温度场,对铸件液态区进行金属液压力P的计算;
2)根据温度场计算结果,确定铸件是否有即将凝固单元即当前单元?如果有,对当前单元进行缩松判据的计算;如果没有,计算程序返回步骤1),继续计算铸件的温度场;对于不是即将凝固时刻的单元不进行缩松判据的计算;
3)根据步骤2)确定的当前单元,计算当前单元与相邻26个单元的距离;
4)根据当前单元与相邻26个单元的距离以及当前单元和相邻单元的固相率差,计算当前单元在相邻26个单元方向上的固相率梯度Gs,并选取固相率梯度最大者为当前单元的固相率梯度Gs;
5)计算当前单元的等价冷却速度R,计算方法为,将凝固模拟过程中一个时间步长内释放的潜热换算成当前单元温度的降低,此温度的降低与该单元实际温度的降低之和除以一个时间步长作为当前单元的冷却速度,称之为等价冷却速度;
所述的与当前单元相邻的26个单元是,面相邻的6个,棱相邻的12个和顶点相邻的8个,当前单元与相邻26个单元的距离是当前单元中心点与相邻26个单元中心点之间的距离。
在铸件凝固模拟中,首先计算出铸件的温度场,温度场的计算采通用的有限差分方法,结晶潜热、边界条件、固相率等有关凝固问题的处理都在这里完成。根据温度场的情况,对液态区域进行金属液压力的计算,压力计算可以采用流场模拟的方式计算,也可以直接利用计算液体静压力的方式计算金属液的静压力。压力计算结束以后,保留至后续的缩松判据使用。某单元缩松判据的计算要在该单元即将凝固的时刻进行,该时刻就是:所计算单元的固相率刚刚达到1的前一个时刻,即刚刚凝固的前一个时刻。对于计算缩松判据的单元来说,需要计算该单元(称当前单元)与相邻26个单元的固相率梯度和它本身的等价冷却速度。首先要计算当前单元与相邻26个单元的距离,与当前单元相邻的26个单元是,面相邻的6个,棱相邻的12个和顶点相邻的8个,当前单元与相邻26个单元的距离是两个单元中心点的距离。固相率梯度是当前单元与相邻26个单元的固相率差除以它们之间的距离,取其中最大者而得。等价冷却速度就是把当前单元一个时间步长内释放的凝固潜热折算为该单元的温度降低的数值与该单元实际温度降低的数值叠加在一起,然后除以一个时间步长即可得到当前单元的等价冷却速度。按照判据公式计算判据值将此值与临界值C相比较,判据值小于等于临界值的单元设置为缩松单元,否则不为缩松单元。对于不同的合金临界值也不同,要由实验确定。
四、附图说明
图1是本发明中确定铸件即将凝固单元的数值模拟方法流程图。
图2是本发明中计算预测缩松判据的方法流程图。
五、具体实施方式
实施例:一种铸件凝固模拟中预测缩松的方法,根据直接影响缩松形成的各个因素,将铸件凝固模拟软件运行于计算机上,在预测缩松判据公式 中,P为金属液的压力,Gs固相率梯度,R是冷却速度(K/min),C是临界值,具体确定步骤为:
1)计算程序开始,利用通用的差分方法计算出铸件的温度场,根据温度场,对铸件液态区进行金属液压力P的计算;
2)根据温度场计算结果,确定铸件是否有即将凝固单元即当前单元?如果有,对当前单元进行缩松判据的计算;如果没有,计算程序返回步骤1),继续计算铸件的温度场;对于不是即将凝固时刻的单元不进行缩松判据的计算;
3)根据步骤2)确定的当前单元,计算当前单元与相邻26个单元的距离;
4)根据当前单元与相邻26个单元的距离以及当前单元和相邻单元的固相率差,计算当前单元在相邻26个单元方向上的固相率梯度Gs,并选取固相率梯度最大者为当前单元的固相率梯度Gs;
5)计算当前单元的等价冷却速度R,计算方法为,将凝固模拟过程中一个时间步长内释放的潜热换算成当前单元温度的降低,此温度的降低与该单元实际温度的降低之和除以一个时间步长作为当前单元的冷却速度,称之为等价冷却速度;
所述的与当前单元相邻的26个单元是,面相邻的6个,棱相邻的12个和顶点相邻的8个,当前单元与相邻26个单元的距离是当前单元中心点与与相邻26个单元中心点之间的距离。
在铸件凝固模拟中,首先计算出铸件的温度场,温度场的计算采通用的有限差分方法,结晶潜热、边界条件、固相率等有关凝固问题的处理都在这里完成。根据温度场的情况,对液态区域进行金属液压力的计算,压力计算可以采用流场模拟的方式计算,也可以直接利用计算液体静压力的方式计算金属液的静压力。压力计算结束以后,保留至后续的缩松判据使用。某单元缩松判据的计算要在该单元即将凝固的时刻进行,该时刻就是:所计算单元的固相率刚刚达到1的前一个时刻,即刚刚凝固的前一个时刻。对于计算缩松判据的单元来说,需要计算该单元(称当前单元)与相邻26个单元的固相率梯度和它本身的等价冷却速度。首先要计算当前单元与相邻26个单元的距离,与当前单元相邻的26个单元是,面相邻的6个,棱相邻的12个和顶点相邻的8个,当前单元与相邻26个单元的距离是两个单元中心点的距离。固相率梯度是当前单元与相邻26个单元的固相率差除以它们之间的距离,取其中最大者而得。等价冷却速度就是把当前单元在一个时间步长内释放的凝固潜热折算为该单元的温度降低的数值与该单元实际温度降低的数值叠加在一起,然后除以一个时间步长即可得到当前单元的等价冷却速度。按照判据公式计算判据值将此值与临界值C相比较,判据值小于等于临界值的单元设置为缩松单元,否则不为缩松单元。
Claims (2)
1、一种铸件凝固模拟中预测缩松的方法,其特征在于,根据预测缩松判据公式 利用计算机,按如下步骤进行:
1)计算程序开始,利用通用的差分方法计算出铸件的温度场,根据温度场,对铸件液态区进行金属液压力P的计算;
2)根据温度场计算结果,确定铸件是否有即将凝固单元,即当前单元,如果有,对当前单元进行缩松判据的计算;如果没有,计算程序返回步骤1),继续计算铸件的温度场;对于不是即将凝固单元的不进行缩松判据的计算;
3)根据步骤2)确定的当前单元,计算当前单元与相邻26个单元的距离;
4)根据当前单元与相邻26个单元的距离以及当前单元和相邻单元的固相率差,计算当前单元在相邻26个单元方向上的固相率梯度Gs,并选取固相率梯度最大者为当前单元的固相率梯度Gs;
5)计算当前单元的等价冷却速度R,计算方法为,将凝固模拟过程中一个时间步长内释放的潜热换算成当前单元温度的降低,此温度的降低与该单元实际温度的降低之和除以一个时间步长作为当前单元的冷却速度,称之为等价冷却速度;
6)根据步骤3)、4)、5)的计算结果,利用公式计算判据值,如果判据值与临界值C的关系满足预测缩松判据公式 即确定当前单元为缩松。
2、如权利要求1所述的一种铸件凝固模拟中预测缩松的方法,其特征在于,所述的与当前单元相邻的26个单元是,面相邻的6个,棱相邻的12个和顶点相邻的8个,当前单元与相邻26个单元的距离是当前单元中心点与相邻26个单元中心点之间的距离。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2006101662298A CN100458801C (zh) | 2006-12-22 | 2006-12-22 | 一种铸件凝固模拟中预测缩松的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2006101662298A CN100458801C (zh) | 2006-12-22 | 2006-12-22 | 一种铸件凝固模拟中预测缩松的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1996319A CN1996319A (zh) | 2007-07-11 |
CN100458801C true CN100458801C (zh) | 2009-02-04 |
Family
ID=38251407
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2006101662298A Expired - Fee Related CN100458801C (zh) | 2006-12-22 | 2006-12-22 | 一种铸件凝固模拟中预测缩松的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100458801C (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101963608A (zh) * | 2009-07-25 | 2011-02-02 | 比亚迪股份有限公司 | 金属液凝固过程中预测缩松的方法及缩松连续预测方法 |
CN102274947A (zh) * | 2011-08-16 | 2011-12-14 | 中北大学 | 一种铝合金低压铸造的缩孔缩松预测方法 |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102072914B (zh) * | 2009-11-23 | 2012-11-21 | 比亚迪股份有限公司 | 金属液凝固过程中预测缩孔的方法 |
CN102236726B (zh) * | 2010-04-23 | 2013-01-30 | 比亚迪股份有限公司 | 一种金属液凝固过程中预测缩松的方法及缩松连续预测方法 |
CN104749194A (zh) * | 2013-12-25 | 2015-07-01 | 沈阳铸造研究所 | 一种预测定向凝固铸件中显微缩松的方法 |
CN107844852B (zh) * | 2017-10-20 | 2021-10-01 | 太原理工大学 | 一种模拟铸钢件砂型铸造过程的缩松缺陷预测方法 |
WO2019087435A1 (ja) * | 2017-11-06 | 2019-05-09 | 株式会社I2C技研 | 鋳造時における凝固解析方法、鋳造方法及び電子プログラム |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006159221A (ja) * | 2004-12-03 | 2006-06-22 | Mitsubishi Materials Corp | マイクロポロシティの発生予測方法、マイクロポロシティ発生予測プログラム及び記録媒体 |
CN1873401A (zh) * | 2006-06-28 | 2006-12-06 | 中国船舶重工集团公司第十二研究所 | 一种基于有限元的球墨铸铁件缩孔缺陷预测方法 |
-
2006
- 2006-12-22 CN CNB2006101662298A patent/CN100458801C/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006159221A (ja) * | 2004-12-03 | 2006-06-22 | Mitsubishi Materials Corp | マイクロポロシティの発生予測方法、マイクロポロシティ発生予測プログラム及び記録媒体 |
CN1873401A (zh) * | 2006-06-28 | 2006-12-06 | 中国船舶重工集团公司第十二研究所 | 一种基于有限元的球墨铸铁件缩孔缺陷预测方法 |
Non-Patent Citations (2)
Title |
---|
铸钢件缩松判据的改进. 田学雷,李光友,李成栋,王桂进.计算机应用,第50卷第6期. 2001 |
铸钢件缩松判据的改进. 田学雷,李光友,李成栋,王桂进. 计算机应用,第50卷第6期. 2001 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101963608A (zh) * | 2009-07-25 | 2011-02-02 | 比亚迪股份有限公司 | 金属液凝固过程中预测缩松的方法及缩松连续预测方法 |
CN102274947A (zh) * | 2011-08-16 | 2011-12-14 | 中北大学 | 一种铝合金低压铸造的缩孔缩松预测方法 |
CN102274947B (zh) * | 2011-08-16 | 2013-04-03 | 中北大学 | 一种铝合金低压铸造的缩孔缩松预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN1996319A (zh) | 2007-07-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100458801C (zh) | 一种铸件凝固模拟中预测缩松的方法 | |
CN109500371A (zh) | 一种板坯动态二冷和轻压下控制系统 | |
CN104493121B (zh) | 一种大方坯连铸生产过程的凝固末端位置在线检测方法 | |
Lyu et al. | Mold simulator study on the initial solidification of molten steel near the corner of continuous casting mold | |
CN104331629A (zh) | 一种连铸结晶器保护渣液态、固态渣膜与气隙厚度非均匀分布的计算方法 | |
CN107052292A (zh) | 一种基于热物性参数分布计算的连铸坯热跟踪计算方法 | |
El-Bealy | On the formation of macrosegregation and interdendritic cracks during dendritic solidification of continuous casting of steel | |
CN109711113B (zh) | 一种大方坯凸形辊连铸重压下偏析行为预测方法 | |
JP6835297B1 (ja) | 鋳型内凝固シェル厚推定装置及び鋳型内凝固シェル厚推定方法 | |
CN102228972B (zh) | 用于连铸结晶器凝固传热过程的计算方法 | |
CN111950107B (zh) | 一种多流连铸中间包内流场流动特性的分析方法 | |
Smyrnov et al. | Application of numerical model of continuous cast bloom crystallization to improve the efficiency of mechanical soft reduction technology | |
CN115455760A (zh) | 铸坯在线热状态跟踪方法 | |
CN101168188A (zh) | 动态轻压下控制方法 | |
CN102078947A (zh) | 用于连铸结晶器凝固传热过程热流密度的计算方法 | |
CN102228974B (zh) | 连铸结晶器内钢液凝固传热过程仿真方法 | |
CN102228971B (zh) | 在线模拟连铸结晶器内钢液凝固传热过程的仿真方法 | |
Ji et al. | Online dynamic control of secondary cooling for the continuous casting process | |
JP2664572B2 (ja) | 連続鋳造における鋳片未凝固部分の温度予測方法 | |
Dautov et al. | On 3D dynamic control of secondary cooling in continuous casting process | |
Xia et al. | Simulation of thermomechanical behavior during continuous casting process based on MiLE method | |
RU2787109C1 (ru) | Устройство для оценки толщины затвердевшей корочки в кристаллизаторе и способ оценки толщины затвердевшей корочки в кристаллизаторе | |
JP5742601B2 (ja) | 鋳片の凝固完了位置算出方法および凝固完了位置算出装置 | |
SARIOĞLU et al. | Simulation of solidification process used in twin roll casting | |
Chen | Control of constrained moving-boundary process with application to steel continuous casting |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20090204 Termination date: 20111222 |