铸钢铸铁件缩孔缩松缺陷预测的数值模拟研究及软件开发
|
摘要
铸件充型凝固过程数值模拟是当前材料加工领域的研究热点。许多发达国家都拥有了自己的商品化软件,已经应用大生产实践当中,国内也出现了几款商品化软件,但与国外水平相比,差距较大。基于此,本文采用较新的技术建立数值模拟的基础框架,总体发展目标是开发出具有自主知识产权的具有国内外先进水平的铸件成型数值模拟软件。本论文研究开发了铸钢、铸铁砂型铸件凝固过程温度场计算模块及缩孔缩松缺陷预测模块,材料数据库模块、后处理数据可视化等模块,进一步完善了前处理模块,基本形成了完整的数值模拟软件框架。
在温度场计算方面,离散模型采用显式直接差分法,采用温度回升法和等效比热容法分区段应用来处理潜热,界面换热系数采用等效界面换热系数,并编制了完善的温度场计算程序。在此基础上完成了铸件凝固过程孤立域的实时快速准确计算。对于一般铸钢件,采用孤立域即缩即补法来预测集中缩孔,用Niyama准则来预测缩松缺陷。对于铸铁件,采用凝固过程收缩、膨胀动态叠加法和最终区域的收缩膨胀法来预测缩孔缩松缺陷。模拟结果与实际铸件浇注结果对比基本吻合。
在材料数据库方面,在全面分析比较当前各种数值模拟软件数据库,采用paradox进行数据库开发,完成了具有增加、删除、修改等功能的简易实用软件材料数据库,并将其成功嵌入到数值模拟软件中。
在数据可视化方面,采用目前最先进的OpenGL技术,在C++Builder开发环境下,编制了完整的前处理、中间计算及后处理图形显示程序,使得图形能够实时快速平移、旋转、缩放、剖切等变换,实现了对铸件的全方位立体观察。并且实现了对铸造凝固过程的凝固状态的动态彩色化显示。
在软件界面开发方面,本文利用C++Builder丰富的构件,编制了一个图形化友好的操作界面,并整合了前处理、可视化模拟计算、后处理计算结果显示和材料数据库等程序,使之成为一个操作简单,功能相对完整的铸造过程数值模拟软件。
Numerical simulation for foundry process is one of the hot research topics in material processing. Many kinds of commercial casting simulation software built by several developed countries have been put into industrial application now, and in contrast, there are also several software poor of strong capabilities created by a few Chinese units. In this paper, by adapting development to the numerical simulation, a numerical profile software including modulus of temperature field simulation, shrinkage defects prediction, material database, data visualization, and etc has been developed. We have independent intellectual property rights for it.
In temperature field solving module, the direct difference scheme was used for numerical discretization, and the temperature compensation method and equivalent specific heat method were used together to deal with the latent heat, and the equivalent coefficient of heat transfer was used to deal with the interface coefficient, so a temperature field simulation module was set up, and also the isolate hot spots during the solidification of casting can be traced rapidly and accurately just in time. For steel castings, ELSDM (Equivalent Liquid Surface Descending Method) and Niyama Method were applied to search the shrinkage and porosity defects; however, for iron castings, DECAM (Dynamic Expansion Contraction Accumulation Method) was employed. Then a few calculation instances verified the validity of them.
For the database module in our software, by investigating several kinds of software from domestic or abroad, a thermal-physical database management system was developed on the base of PARADOX, with functions of adding, deleting, editing and querying data records in lists of the software database.
In data visualization module, OpenGL, an advanced graphics technology, was adopted to set up the data visualization program under the developing environment of C++Builder. At last the functions such as the graphic’s real-time translation, rotation, and scaling were carried out, which makes it possible to perceive the whole casting from all aspects. In addition, the dynamic process of casting solidification in color display was achieved.
In order to construct software interface, the abundant components of C++Builder was used, and then a friendly interface was built up, including the pre-process module, the temperature calculation module, the database system, and the post-process module. In the end, the easy operating casting simulation software was finished.
在温度场计算方面,离散模型采用显式直接差分法,采用温度回升法和等效比热容法分区段应用来处理潜热,界面换热系数采用等效界面换热系数,并编制了完善的温度场计算程序。在此基础上完成了铸件凝固过程孤立域的实时快速准确计算。对于一般铸钢件,采用孤立域即缩即补法来预测集中缩孔,用Niyama准则来预测缩松缺陷。对于铸铁件,采用凝固过程收缩、膨胀动态叠加法和最终区域的收缩膨胀法来预测缩孔缩松缺陷。模拟结果与实际铸件浇注结果对比基本吻合。
在材料数据库方面,在全面分析比较当前各种数值模拟软件数据库,采用paradox进行数据库开发,完成了具有增加、删除、修改等功能的简易实用软件材料数据库,并将其成功嵌入到数值模拟软件中。
在数据可视化方面,采用目前最先进的OpenGL技术,在C++Builder开发环境下,编制了完整的前处理、中间计算及后处理图形显示程序,使得图形能够实时快速平移、旋转、缩放、剖切等变换,实现了对铸件的全方位立体观察。并且实现了对铸造凝固过程的凝固状态的动态彩色化显示。
在软件界面开发方面,本文利用C++Builder丰富的构件,编制了一个图形化友好的操作界面,并整合了前处理、可视化模拟计算、后处理计算结果显示和材料数据库等程序,使之成为一个操作简单,功能相对完整的铸造过程数值模拟软件。
Numerical simulation for foundry process is one of the hot research topics in material processing. Many kinds of commercial casting simulation software built by several developed countries have been put into industrial application now, and in contrast, there are also several software poor of strong capabilities created by a few Chinese units. In this paper, by adapting development to the numerical simulation, a numerical profile software including modulus of temperature field simulation, shrinkage defects prediction, material database, data visualization, and etc has been developed. We have independent intellectual property rights for it.
In temperature field solving module, the direct difference scheme was used for numerical discretization, and the temperature compensation method and equivalent specific heat method were used together to deal with the latent heat, and the equivalent coefficient of heat transfer was used to deal with the interface coefficient, so a temperature field simulation module was set up, and also the isolate hot spots during the solidification of casting can be traced rapidly and accurately just in time. For steel castings, ELSDM (Equivalent Liquid Surface Descending Method) and Niyama Method were applied to search the shrinkage and porosity defects; however, for iron castings, DECAM (Dynamic Expansion Contraction Accumulation Method) was employed. Then a few calculation instances verified the validity of them.
For the database module in our software, by investigating several kinds of software from domestic or abroad, a thermal-physical database management system was developed on the base of PARADOX, with functions of adding, deleting, editing and querying data records in lists of the software database.
In data visualization module, OpenGL, an advanced graphics technology, was adopted to set up the data visualization program under the developing environment of C++Builder. At last the functions such as the graphic’s real-time translation, rotation, and scaling were carried out, which makes it possible to perceive the whole casting from all aspects. In addition, the dynamic process of casting solidification in color display was achieved.
In order to construct software interface, the abundant components of C++Builder was used, and then a friendly interface was built up, including the pre-process module, the temperature calculation module, the database system, and the post-process module. In the end, the easy operating casting simulation software was finished.
引文
[1] 李维容,高维钢.从《铸造》50 年看我国铸造业的发展,铸造,2002(2).
[2] 邱宗文. 铸件凝固过程温度场三维数值模拟软件开发. 硕士学位论文,东南大学,2004.3
[3] 陈海清 李华基 曹阳等 铸件凝固过程数值模拟 重庆大学出版社 1991 年 4 月
[4] 熊守美 许庆彦 康进武 铸造过程模拟仿真技术 北京 机械工业出版社 2004.9
[5] 曹洪吉,宋延沛,王文焱. 铸造过程计算机模拟研究应用现状与发展. 河南科技大学学报, 2006(1):5～8
[6] 崔吉顺.大型铸件缩孔缩松形成多种预测判据的研究.硕士学位论文,清华大学工学机械工程,2001.5
[7] 曹辉. 铸件凝固进程的数值模拟方法研究. 硕士学位论文,大连理工大学,2003.6
[8] 王君卿,孙逊. MCSP622004 国际会议概况和铸造过程数值模拟技术发展动向. 铸造, 2005(4):307～313
[9] 朱鸣芳,于金,戴挺. 金属凝固过程数值模拟的最新进展——第10 届MCWASP 国际学术会92议论文述评. 铸造,2005(12):115～120
[10] 大中逸雄. 铸物,1981(53):673～679.
[11] 千千岩健儿. 铸物,1983(55):271～278.
[12] 李日. 专家系统及数值模拟在铸造工艺设计中应用的研究,博士学位论文,西北工业大学,1998 年
[13] 吴孟怀.铸造工艺计算机辅助设计基础,西安:西北工业大学,1993
[14] Kubo etal. Mathematical Modelling of prosity Formation in Solidification Metallurgical Transactions,1985,16B :359～366
[15] Davies V del. Feeding Range Determination by Numerically Computed Heat distribution. AFS Cast Metals Research J. 1975,11(2) :33～44
[16] Niyama E,etal. Predicting Shrinkage in Large Steel Castings from Temperature Gradiert Calculation. AFS International Cast Meetals J. 1981. 6(2):PP10～22
[17] 贾宝仟等. 逐渐缩松缩孔判据 G/ T 的理论基础及应用.铸造,1996(4):13～15.
[18] 贾宝仟等.逐渐缩松缩孔判据 G/ T 临界值的实验研究.铸造技术,1997(1):41-42.
[19] 李文珍等.铸钢件缩孔缩松预测的实验研究.铸造技术,1995(4): 37-40.
[20] 王春乐. 铸钢件缩孔缩松预测方法及判据浅析. 山西机械,2003(4):8～10
[21] 田学雷,李光友,李成栋. 铸钢件缩松判据的改进. 铸造,2001(6):346～348
[22] S.I.Karsay 著.白天申,林家骝,李延龄译,球墨铸铁浇口和冒口.清华大学出版社.1983
[23] 李嘉荣,柳百成,刘睿. G / R法在球铁件缩孔缩松预测中的应用研究.铸造.1998(7): 9-13
[24] Y.-F.Chen etal.Improvements on Solidification Modeling for A356.2 Aluminum Alloy.AFS.Trans.96-180
[25] 李邦盛,狂飞,吴士平. TiAl合金熔模精铸件缩孔形成的数值模拟及预测. 特种铸造及有色 合金,2005(5):280～283
[26] 徐宏,程军. 铝合金铸件缩孔缩松预测新技术. 中国有色金属学报,2001(1):34～37
[27] 林首位,徐宏,侯华. 铸造CAE工程数据库开发技术研究. 华北工学院学报, 2001(1):20～24
[28] 杨秉俭, 王理勤. 铸造材料热物性值数据库管理系统的创建. 铸造技术, 1998(6):9～11
[29] 王向群,刘瑞祥. 铸造CAE热物性参数库的设计与开发. 中国铸造装备与技术,2004(01)
[30] 邱宗文, 仇春荔, 潘冶. 数据库技术在铸造业中的应用. 铸造,2003(1):45～48
[31] 王喜军,安阁英.铸造热物性参数数据库系统的开发与研究.中国铸造装备与技术,1998(1):11～14
[32] 刘伟涛,李万,杨秉俭. 铸件凝固过程数值模拟后处理的研究. 西安交通大学学报, 1999(10):52～56
[33] 怀刚,李梅娥. 基于有限差分网格的凝固模拟后处理研究. 铸造,2003(5):350～352
[34] 刘晶峰,刘瑞祥. 铸造CAE系统后置处理模块中“无锯齿化”算法. 铸造,2002(11):706～708
[35] 胡志勇,李殿中. 铸件凝固过程模拟结果的三维OpenGL显示. 铸造,1999(10):19～21
[36] 曹敏.等值面的高质量显示算法实现三维物体表面的重建.电脑开发与应用,2002(14): 14～15
[37] 吴树森. 材料加工冶金传输原理.北京:机械工业出版社,2001.5
[38] 张波, 于九明, 由宝财. 铸造过程温度场的数值模拟. 中国铸造装备与技术, 2004(2)
[39] 陈国权. 铝合金铸件/金属型界面间隙及换热系数通用化的研究. 铸造,1993(8):1～6
[40] 安晓卫,王承志,宋广胜. 铸件凝固温度场有限元分析中界面热阻的处理.计算力学学报,2005(1)
[41] 宋广胜,路明. 铸件凝固温度场有限元分析中界面热阻处理. 沈阳航空工业学院学报, 2002(4):20～22
[42] 殷刚,章志超. 相变蓄冷材料的选择与相变潜热的测定,哈尔滨商业大学学报,2004(4)
[43] 周建兴,刘瑞祥,陈立亮. 凝固过程数值模拟中的潜热处理方法. 铸造,2001(7):404～407
[44] 白云峰,徐达鸣,郭景杰. 采用温度回升法对任意结晶区间的铸件凝固结晶潜热的数值计算. 金属学报,2003(6).623～629
[45] 刘永刚,陈 光,孙国雄. 一种新型潜热处理模型显微偏析数值模拟的应用. 特种铸造及有色合金, 2005(5)
[46] 李东辉,高云宝,辛启斌. 铸件凝固潜热的处理方法与应用研究. 铸造,2004(12)
[47] 刘增林. 铸铁件热节区的凝固与补缩. 铸造,1993(3):29～31
[48] 康进武,熊守美,柳百成. 采用多热节和即缩即补方法预测铸钢件缩孔的研究. 铸造. 2000(8):478～481
[49] 贾宝仟,熊守美,柳百成. 铸件凝固过程孤立域动态划分及缩孔缩松数值模拟. 铸造技术,1996(5):15～17
[50] 徐宏,侯华,赵宇宏. 铸件凝固模拟液态孤立多熔池动态判定及缩孔预测. 铸造技术, 2003(4):319～320
[51] 田学雷,李成栋. 铸钢件凝固模拟中对多个缩孔的预测技术. 中国铸造装备与技术,2001(5):16～17
[52] 陆文华,李隆盛.铸造合金及其熔炼.北京:机械工业出版社,1997(10)
[53] 刘伟涛,杨秉俭,朱静林. 球墨铸铁件收缩缺陷判据的研究. 西安交通大学学报,2000(3):78～82
[54] 周亘. 球墨铸铁件冒口补缩失败原因分析--球墨铸铁缩松、缩孔问题探讨.现代铸铁,2004(4):7～14
[55] 杨杰,董怀宇,熊守美. 灰铁凝固过程中缩孔缩松的预测. 金属学报,2005(9):929～932
[56] 宋玉华,任善之,冯毓青. 球铁件三维缩孔预测的实验研究及其在计算机上的实现. 材料科学与工艺,1997(1):114～116
[57] Richard S.Wright,Jr. Benjamin Lipchak 著,徐波 译. OpenGL 超级宝典.北京:人民邮电出版社 (第三版), 2005.9
[58] Edward Angel. OpenGL 程序设计指南(第 2 版).北京: 清华大学出版社,2005
[59] Mike Bailey.Teaching OpenGl shaders:Hands-on, interactive, andimmediate feedback, Computers & Graphics .2007.02
[60] 李胜睿. 计算机图形学实验教程OpenGL版. 北京:机械工业出版社,2004
[61] 饶文碧,陈旭辉,徐锐. 基于C++Builder的图形图像处理方法. 交通与计算机,2002(2):61～63
[62] 万方,熊家军. 基于C++Builder6.0实现OpenGL程序框架. 电脑学习,2006(2):34～36
[63] Bo Chen,Harry H.Cheng. Interpretive OpenGL for computer graphics. Computers and Graphics.2005.29: PP331～339
[64] 孟飞. 基于OpenGL 的应用与探讨. 科技广场,2006(4):63～64
[65] 龚安,宫法明. 用OpenGL实现交互式的三维图形显示. 科学技术与工程,2006(6):1428～1430
[66] YIN Aijun,QIN Shuren,etal.Visualization of three-dimensional data field and its application in machine testing.Chinese Journal of Mechanical Engineering. 2006. 1:PP81～84
[67] 焦藏桢,董晟全. OpenGL 技术在金属凝固数值模拟图形中的应用. 铸造设备研究,2004(6):32～34
[68] 李文强,华祖耀. 基于OpenGL的仿真可视化研究. 计算机仿真,2005(6):158～160
[69] 欧训勇. OpenGL程序移植技术浅析. 计算机应用研究,2002(9):158～160
[70] 励映群,王润杰. OpenGL中实现飞行仿真观察变换的几种方法. 计算机仿真,2003(5):71～73
[71] 王胜全,范燮昌,谈士力. OpenGL中的坐标变换分析. 机械与电子,2004(2):65～67
[72] 林茂松,尹显明,夏季. OpenGL中矩阵和矩阵堆栈的使用方法. 机械,2000(27):97
[73] 邱宗文,孙国雄. 金属凝固数值模拟可视化中OpenGL的应用研究. 铸造,2003(12):1173～1175
[74] 崔洪斌,唐俊源. 利用OpenGL实现动态观察. 现代制造工程,2004(7):44～45
[75] 王凤林,陈立亮. OpenGL技术在铸造CAE后处理模块中的应用. 铸造,2002(12):775～777
[76] 陶永劲,陈立亮. OpenGL技术在铸造CAE前处理模块中的应用. 特种铸造及有色合金,2001(6):23～24
[77] 王福生,齐国清. 二值图像中目标物体轮廓的边界跟踪算法. 大连海事大学学报,2006(32): 62～67
[78] 孟飞. 基于OpenGL 的计算机动画技术. 科技广场,2005(1):109～110
[79] 冯渊. 利用OpenGL实现RGBA模式下的光照处理. 福建电脑,2004(4):23～24
[80] 宋永恒,荆涛,周舰. 三维凝固过程温度场可视化研究. 新技术新工艺,2001(12):8～10
[81] 陈小全. 基于OpenGL的钢坯温度场可视化研究及其应用. 硕士学位论文,昆明理工大学2005.7
[82] 牛连强,王金东,杨放. OpenGL图形消隐质量的改善. 计算机应用. 2001(4):36～38
[83] 赵卫东,卫 刚. 在OpenGL下面消隐和线消隐的实现. 计算机工程,2002(6):233～234
[84] 张仁忠,王志勋. 基于OpenGL 双缓存动画仿真软件的开发. 应用科技,2003(7):38～40
[85] 贺再红,杨鼎强. 基于OpenGL双缓存的多媒体动画实现. 计算技术与自动化,2001(1):4～38
[86] 余明兴,吴明哲. Borland C++ Builder6 程序设计经典.北京:科学出版社,2004.1
[87] 孙江涛.例说 C++Builder.北京:北京大学出版社,2000.4
[2] 邱宗文. 铸件凝固过程温度场三维数值模拟软件开发. 硕士学位论文,东南大学,2004.3
[3] 陈海清 李华基 曹阳等 铸件凝固过程数值模拟 重庆大学出版社 1991 年 4 月
[4] 熊守美 许庆彦 康进武 铸造过程模拟仿真技术 北京 机械工业出版社 2004.9
[5] 曹洪吉,宋延沛,王文焱. 铸造过程计算机模拟研究应用现状与发展. 河南科技大学学报, 2006(1):5～8
[6] 崔吉顺.大型铸件缩孔缩松形成多种预测判据的研究.硕士学位论文,清华大学工学机械工程,2001.5
[7] 曹辉. 铸件凝固进程的数值模拟方法研究. 硕士学位论文,大连理工大学,2003.6
[8] 王君卿,孙逊. MCSP622004 国际会议概况和铸造过程数值模拟技术发展动向. 铸造, 2005(4):307～313
[9] 朱鸣芳,于金,戴挺. 金属凝固过程数值模拟的最新进展——第10 届MCWASP 国际学术会92议论文述评. 铸造,2005(12):115～120
[10] 大中逸雄. 铸物,1981(53):673～679.
[11] 千千岩健儿. 铸物,1983(55):271～278.
[12] 李日. 专家系统及数值模拟在铸造工艺设计中应用的研究,博士学位论文,西北工业大学,1998 年
[13] 吴孟怀.铸造工艺计算机辅助设计基础,西安:西北工业大学,1993
[14] Kubo etal. Mathematical Modelling of prosity Formation in Solidification Metallurgical Transactions,1985,16B :359～366
[15] Davies V del. Feeding Range Determination by Numerically Computed Heat distribution. AFS Cast Metals Research J. 1975,11(2) :33～44
[16] Niyama E,etal. Predicting Shrinkage in Large Steel Castings from Temperature Gradiert Calculation. AFS International Cast Meetals J. 1981. 6(2):PP10～22
[17] 贾宝仟等. 逐渐缩松缩孔判据 G/ T 的理论基础及应用.铸造,1996(4):13～15.
[18] 贾宝仟等.逐渐缩松缩孔判据 G/ T 临界值的实验研究.铸造技术,1997(1):41-42.
[19] 李文珍等.铸钢件缩孔缩松预测的实验研究.铸造技术,1995(4): 37-40.
[20] 王春乐. 铸钢件缩孔缩松预测方法及判据浅析. 山西机械,2003(4):8～10
[21] 田学雷,李光友,李成栋. 铸钢件缩松判据的改进. 铸造,2001(6):346～348
[22] S.I.Karsay 著.白天申,林家骝,李延龄译,球墨铸铁浇口和冒口.清华大学出版社.1983
[23] 李嘉荣,柳百成,刘睿. G / R法在球铁件缩孔缩松预测中的应用研究.铸造.1998(7): 9-13
[24] Y.-F.Chen etal.Improvements on Solidification Modeling for A356.2 Aluminum Alloy.AFS.Trans.96-180
[25] 李邦盛,狂飞,吴士平. TiAl合金熔模精铸件缩孔形成的数值模拟及预测. 特种铸造及有色 合金,2005(5):280～283
[26] 徐宏,程军. 铝合金铸件缩孔缩松预测新技术. 中国有色金属学报,2001(1):34～37
[27] 林首位,徐宏,侯华. 铸造CAE工程数据库开发技术研究. 华北工学院学报, 2001(1):20～24
[28] 杨秉俭, 王理勤. 铸造材料热物性值数据库管理系统的创建. 铸造技术, 1998(6):9～11
[29] 王向群,刘瑞祥. 铸造CAE热物性参数库的设计与开发. 中国铸造装备与技术,2004(01)
[30] 邱宗文, 仇春荔, 潘冶. 数据库技术在铸造业中的应用. 铸造,2003(1):45～48
[31] 王喜军,安阁英.铸造热物性参数数据库系统的开发与研究.中国铸造装备与技术,1998(1):11～14
[32] 刘伟涛,李万,杨秉俭. 铸件凝固过程数值模拟后处理的研究. 西安交通大学学报, 1999(10):52～56
[33] 怀刚,李梅娥. 基于有限差分网格的凝固模拟后处理研究. 铸造,2003(5):350～352
[34] 刘晶峰,刘瑞祥. 铸造CAE系统后置处理模块中“无锯齿化”算法. 铸造,2002(11):706～708
[35] 胡志勇,李殿中. 铸件凝固过程模拟结果的三维OpenGL显示. 铸造,1999(10):19～21
[36] 曹敏.等值面的高质量显示算法实现三维物体表面的重建.电脑开发与应用,2002(14): 14～15
[37] 吴树森. 材料加工冶金传输原理.北京:机械工业出版社,2001.5
[38] 张波, 于九明, 由宝财. 铸造过程温度场的数值模拟. 中国铸造装备与技术, 2004(2)
[39] 陈国权. 铝合金铸件/金属型界面间隙及换热系数通用化的研究. 铸造,1993(8):1～6
[40] 安晓卫,王承志,宋广胜. 铸件凝固温度场有限元分析中界面热阻的处理.计算力学学报,2005(1)
[41] 宋广胜,路明. 铸件凝固温度场有限元分析中界面热阻处理. 沈阳航空工业学院学报, 2002(4):20～22
[42] 殷刚,章志超. 相变蓄冷材料的选择与相变潜热的测定,哈尔滨商业大学学报,2004(4)
[43] 周建兴,刘瑞祥,陈立亮. 凝固过程数值模拟中的潜热处理方法. 铸造,2001(7):404～407
[44] 白云峰,徐达鸣,郭景杰. 采用温度回升法对任意结晶区间的铸件凝固结晶潜热的数值计算. 金属学报,2003(6).623～629
[45] 刘永刚,陈 光,孙国雄. 一种新型潜热处理模型显微偏析数值模拟的应用. 特种铸造及有色合金, 2005(5)
[46] 李东辉,高云宝,辛启斌. 铸件凝固潜热的处理方法与应用研究. 铸造,2004(12)
[47] 刘增林. 铸铁件热节区的凝固与补缩. 铸造,1993(3):29～31
[48] 康进武,熊守美,柳百成. 采用多热节和即缩即补方法预测铸钢件缩孔的研究. 铸造. 2000(8):478～481
[49] 贾宝仟,熊守美,柳百成. 铸件凝固过程孤立域动态划分及缩孔缩松数值模拟. 铸造技术,1996(5):15～17
[50] 徐宏,侯华,赵宇宏. 铸件凝固模拟液态孤立多熔池动态判定及缩孔预测. 铸造技术, 2003(4):319～320
[51] 田学雷,李成栋. 铸钢件凝固模拟中对多个缩孔的预测技术. 中国铸造装备与技术,2001(5):16～17
[52] 陆文华,李隆盛.铸造合金及其熔炼.北京:机械工业出版社,1997(10)
[53] 刘伟涛,杨秉俭,朱静林. 球墨铸铁件收缩缺陷判据的研究. 西安交通大学学报,2000(3):78～82
[54] 周亘. 球墨铸铁件冒口补缩失败原因分析--球墨铸铁缩松、缩孔问题探讨.现代铸铁,2004(4):7～14
[55] 杨杰,董怀宇,熊守美. 灰铁凝固过程中缩孔缩松的预测. 金属学报,2005(9):929～932
[56] 宋玉华,任善之,冯毓青. 球铁件三维缩孔预测的实验研究及其在计算机上的实现. 材料科学与工艺,1997(1):114～116
[57] Richard S.Wright,Jr. Benjamin Lipchak 著,徐波 译. OpenGL 超级宝典.北京:人民邮电出版社 (第三版), 2005.9
[58] Edward Angel. OpenGL 程序设计指南(第 2 版).北京: 清华大学出版社,2005
[59] Mike Bailey.Teaching OpenGl shaders:Hands-on, interactive, andimmediate feedback, Computers & Graphics .2007.02
[60] 李胜睿. 计算机图形学实验教程OpenGL版. 北京:机械工业出版社,2004
[61] 饶文碧,陈旭辉,徐锐. 基于C++Builder的图形图像处理方法. 交通与计算机,2002(2):61～63
[62] 万方,熊家军. 基于C++Builder6.0实现OpenGL程序框架. 电脑学习,2006(2):34～36
[63] Bo Chen,Harry H.Cheng. Interpretive OpenGL for computer graphics. Computers and Graphics.2005.29: PP331～339
[64] 孟飞. 基于OpenGL 的应用与探讨. 科技广场,2006(4):63～64
[65] 龚安,宫法明. 用OpenGL实现交互式的三维图形显示. 科学技术与工程,2006(6):1428～1430
[66] YIN Aijun,QIN Shuren,etal.Visualization of three-dimensional data field and its application in machine testing.Chinese Journal of Mechanical Engineering. 2006. 1:PP81～84
[67] 焦藏桢,董晟全. OpenGL 技术在金属凝固数值模拟图形中的应用. 铸造设备研究,2004(6):32～34
[68] 李文强,华祖耀. 基于OpenGL的仿真可视化研究. 计算机仿真,2005(6):158～160
[69] 欧训勇. OpenGL程序移植技术浅析. 计算机应用研究,2002(9):158～160
[70] 励映群,王润杰. OpenGL中实现飞行仿真观察变换的几种方法. 计算机仿真,2003(5):71～73
[71] 王胜全,范燮昌,谈士力. OpenGL中的坐标变换分析. 机械与电子,2004(2):65～67
[72] 林茂松,尹显明,夏季. OpenGL中矩阵和矩阵堆栈的使用方法. 机械,2000(27):97
[73] 邱宗文,孙国雄. 金属凝固数值模拟可视化中OpenGL的应用研究. 铸造,2003(12):1173～1175
[74] 崔洪斌,唐俊源. 利用OpenGL实现动态观察. 现代制造工程,2004(7):44～45
[75] 王凤林,陈立亮. OpenGL技术在铸造CAE后处理模块中的应用. 铸造,2002(12):775～777
[76] 陶永劲,陈立亮. OpenGL技术在铸造CAE前处理模块中的应用. 特种铸造及有色合金,2001(6):23～24
[77] 王福生,齐国清. 二值图像中目标物体轮廓的边界跟踪算法. 大连海事大学学报,2006(32): 62～67
[78] 孟飞. 基于OpenGL 的计算机动画技术. 科技广场,2005(1):109～110
[79] 冯渊. 利用OpenGL实现RGBA模式下的光照处理. 福建电脑,2004(4):23～24
[80] 宋永恒,荆涛,周舰. 三维凝固过程温度场可视化研究. 新技术新工艺,2001(12):8～10
[81] 陈小全. 基于OpenGL的钢坯温度场可视化研究及其应用. 硕士学位论文,昆明理工大学2005.7
[82] 牛连强,王金东,杨放. OpenGL图形消隐质量的改善. 计算机应用. 2001(4):36～38
[83] 赵卫东,卫 刚. 在OpenGL下面消隐和线消隐的实现. 计算机工程,2002(6):233～234
[84] 张仁忠,王志勋. 基于OpenGL 双缓存动画仿真软件的开发. 应用科技,2003(7):38～40
[85] 贺再红,杨鼎强. 基于OpenGL双缓存的多媒体动画实现. 计算技术与自动化,2001(1):4～38
[86] 余明兴,吴明哲. Borland C++ Builder6 程序设计经典.北京:科学出版社,2004.1
[87] 孙江涛.例说 C++Builder.北京:北京大学出版社,2000.4