Please wait a minute...
金属学报  2026, Vol. 62 Issue (3): 509-522    DOI: 10.11900/0412.1961.2024.00406
  研究论文 本期目录 | 过刊浏览 |
Al-4.7%Cu合金枝晶凝固动态收缩的数值模拟
朱宝丰, 李晨宇, 张士杰, 李日()
河北工业大学 材料科学与工程学院 天津 300401
Numerical Simulation of the Dynamic Contraction of Dendrite Solidification in the Al-4.7%Cu Alloy
ZHU Baofeng, LI Chenyu, ZHANG Shijie, LI Ri()
College of Materials Science and Engineering, Hebei University of Technology, Tianjin 300401, China
引用本文:

朱宝丰, 李晨宇, 张士杰, 李日. Al-4.7%Cu合金枝晶凝固动态收缩的数值模拟[J]. 金属学报, 2026, 62(3): 509-522.
Baofeng ZHU, Chenyu LI, Shijie ZHANG, Ri LI. Numerical Simulation of the Dynamic Contraction of Dendrite Solidification in the Al-4.7%Cu Alloy[J]. Acta Metall Sin, 2026, 62(3): 509-522.

全文: PDF(2758 KB)   HTML
摘要: 

合金铸件中微观缩孔的动态形成过程与枝晶凝固过程直接相关。为了模拟Al-4.7%Cu (质量分数)合金在凝固过程中发生的枝晶收缩现象,本工作建立了一个元胞自动机-格子Boltzmann方法(CA-LBM)耦合模型。该模型利用CA模拟枝晶生长过程中缩孔的形成,并通过LBM研究缩孔在液相中的扩散过程。首先用真空孔洞在液体中扩散均匀化的数值模拟验证了所建立的CA-LBM数值模型的正确性,然后对比计算了考虑和不考虑枝晶收缩的单枝晶凝固过程,最后计算了考虑收缩的多枝晶凝固收缩过程。对单枝晶凝固过程的计算结果表明,不考虑收缩的单枝晶在凝固结束后内部不存在孔洞,而考虑收缩的单枝晶在凝固结束后出现了均匀的微小缩孔;孔洞的出现对单枝晶的生长形貌有显著的影响,即会促进单枝晶二次分枝的出现。对多枝晶收缩模拟的结果也表明,微小孔洞出现在最后凝固的枝晶间交界处。缩孔总量的计算值与理论值的误差很小,证明所建立的数值模型是有效可靠的。

关键词 Al-Cu合金凝固枝晶收缩CA-LBM数值模型    
Abstract

The dynamic formation process of microshrinkage pole in alloy castings is directly related to the dendrite solidification process. To simulate dendrite shrinkage in the Al-4.7%Cu (mass fraction) alloy during solidification, we proposed a coupling model combining cellular automata (CA) and the lattice Boltzmann method (LBM), referred to as the CA-LBM model. In this model, CA was used to simulate the formation of shrinkage pores during dendrite growth, whereas LBM was applied to study the diffusion process of shrinkage pores in the liquid phase (fully liquid conditions). First, the accuracy of the proposed CA-LBM numerical model was verified through the numerical simulation of the diffusion homogenization of a vacuum cavity in liquid. Then, the solidification process of a single dendrite—with and without dendrite shrinkage were compared, followed by the calculations of the multi-dendrite solidification contraction process with and without dendrite shrinkage. Simulations of the single-dendrite solidification process indicated that no internal pores were formed in the single dendrites when shrinkage was not considered. However, when shrinkage was considered, uniform microshrinkage pores appeared in the single dendrites. Moreover, the presence of shrinkage pores notably influenced dendrite morphology by promoting secondary branching. The results of the multidendrite contraction simulation also showed that micro-shrinkage pores tended to form at the junctions of the last solidified dendrites. A comparison between the calculated number of shrinkage poles and the theoretical value showed a small error, indicating the effectiveness and reliability of the proposed numerical model.

Key wordsAl-Cu alloy    solidification    dendrite shrinkage    CA-LBM numerical model
收稿日期: 2024-11-28     
ZTFLH:  TG111.4  
基金资助:国家自然科学基金项目(51975182)
通讯作者: 李 日,hbcllr@hebut.edu.cn,主要从事晶体凝固过程及组织控制的多相多物理场多尺度的宏微观耦合数值模拟研究
Corresponding author: LI Ri, professor, Tel: (022)60202006, E-mail: hbcllr@hebut.edu.cn
作者简介: 朱宝丰,男,1997年生,硕士生
图1  微观缩松算法模型
ParameterVariableUnitValue
Solidus temperatureTmK917
Liquidus temperatureTlK821.7
Density of liquid phaseρlkg·m-32606
Solid densityρskg·m-32743
Solute diffusion coefficientDLm-2·s-13.0 × 10-9
Latent heat of crystallizationLJ·kg-1397000
Liquidus slopemlm·K·%-1-3.44
Liquid specific heat capacitycpJ·kg-1·K-11100
Thermal diffusion coefficientαm-2·s-12.7 × 10-7
Equilibrium distribution coefficientk-0.145
Liquid viscosityνm-2·s-11.2 × 10-6
Anisotropy coefficientε-0.0467
Gibbs-Thomson coefficientΓm·K2.4 × 10-7
表1  Al-4.7%Cu的热物性参数
图2  模型流程图
图3  标准反弹格式
图4  九网格下的真空孔洞的扩散模型
图5  真空孔洞扩散过程的九网格验证
图6  真空孔洞扩散过程的二十五网格和一百网格验证
图7  真空孔洞扩散过程的多网格流场的计算域
图8  纯液相中的缩孔扩散(前3000步)
图9  纯液相中的缩孔扩散(3000步后)
图10  多网格流场缩孔扩散示意图(含障碍物)
图11  两点扩散过程示意图
图12  单枝晶生长的模拟结果(不考虑收缩、考虑缩孔产生及扩散)
图13  单枝晶计算中得到的缩孔总量与理论值的对比
图14  多枝晶模拟结果(25个形核点位置及不同位向的多晶粒)
图15  多枝晶生长示意图
图16  多枝晶计算中得到的缩孔总量与理论值的对比
[1] Zhu J D, Cockcroft S L, Maijer D M. Modeling of microporosity formation in A356 aluminum alloy casting [J]. Metall. Mater. Trans., 2006, 37A: 1075
[2] Lee P D, Chirazi A, See D. Modeling microporosity in aluminum-silicon alloys: A review [J]. J. Light Met., 2001, 1: 15
doi: 10.1016/S1471-5317(00)00003-1
[3] Dong S Y, Zhuang Z, Xu Q Y, et al. Preliminary study on simulation of microporosity evolution and fatigue life of aluminum alloy casting [J]. Int. J. Cast Met. Res., 2003, 15: 309
doi: 10.1080/13640461.2003.11819503
[4] Nabavizadeh S A, Eshraghi M, Felicelli S D. A comparative study of multiphase lattice Boltzmann methods for bubble-dendrite interaction during solidification of alloys [J]. Appl. Sci., 2019, 9: 57
doi: 10.3390/app9010057
[5] Khalajzadeh V, Beckermann C. Simulation of shrinkage porosity formation during alloy solidification [J]. Metall. Mater. Trans., 2020, 51A: 2239
[6] Kannan P, Balasubramanian K, Rajeswari N. Process optimization for pore defect in casting—A review [J]. Appl. Mech. Mater., 2015, 813-814: 581
doi: 10.4028/www.scientific.net/AMM
[7] Lee P D, Hunt J D. Hydrogen porosity in directional solidified aluminium-copper alloys: In situ observation [J]. Acta Mater., 1997, 45: 4155
doi: 10.1016/S1359-6454(97)00081-5
[8] Lee P D, Hunt J D. Measuring the nucleation of hydrogen porosity during the solidification of aluminium-copper alloys [J]. Scr. Mater., 1997, 36: 399
doi: 10.1016/S1359-6462(96)00411-3
[9] Atwood R C, Sridhar S, Lee P D. Equations for nucleation of hydrogen gas pores during solidification of aluminium seven weight percent silicon alloy [J]. Scr. Mater., 1999, 41: 1255
doi: 10.1016/S1359-6462(99)00297-3
[10] Atwood R C, Sridhar S, Zhang W, et al. Diffusion-controlled growth of hydrogen pores in aluminium-silicon castings: In situ observation and modelling [J]. Acta Mater., 2000, 48: 405
doi: 10.1016/S1359-6454(99)00363-8
[11] See D, Atwood R C, Lee P D. A comparison of three modeling approaches for the prediction of microporosity in aluminum-silicon alloys [J]. J. Mater. Sci., 2001, 36: 3423
doi: 10.1023/A:1017955927601
[12] Atwood R C, Lee P D. A three-phase model of hydrogen pore formation during the equiaxed dendritic solidification of aluminum-silicon alloys [J]. Metall. Mater. Trans., 2002, 33B: 209
[13] Lee P D, Atwood R C, Dashwood R J, et al. Modeling of porosity formation in direct chill cast aluminum-magnesium alloys [J]. Mater. Sci. Eng., 2002, A328: 213
[14] Hamilton R W, See D, Butler S, et al. Multiscale modeling for the prediction of casting defects in investment cast aluminum alloys [J]. Mater. Sci. Eng., 2003, A343: 290
[15] Atwood R C, Lee P D. Simulation of the three-dimensional morphology of solidification porosity in an aluminium-silicon alloy [J]. Acta Mater., 2003, 51: 5447
doi: 10.1016/S1359-6454(03)00411-7
[16] Chaijaruwanich A, Lee P D, Dashwood R J, et al. Evolution of pore morphology and distribution during the homogenization of direct chill cast Al-Mg alloys [J]. Acta Mater., 2007, 55: 285
doi: 10.1016/j.actamat.2006.08.023
[17] Kareh K M, Lee P D, Atwood R C, et al. Pore behaviour during semi-solid alloy compression: Insights into defect creation under pressure [J]. Scr. Mater., 2014, 89: 73
doi: 10.1016/j.scriptamat.2014.06.033
[18] Puncreobutr C, Lee P D, Kareh K M, et al. Influence of Fe-rich intermetallics on solidification defects in Al-Si-Cu alloys [J]. Acta Mater., 2014, 68: 42
doi: 10.1016/j.actamat.2014.01.007
[19] Pequet C, Rappaz M, Gremaud M. Modeling of microporosity, macroporosity, and pipe-shrinkage formation during the solidification of alloys using a mushy-zone refinement method: Applications to aluminum alloys [J]. Metall. Mater. Trans., 2002, 33A: 2095
[20] Stefanescu D M. Computer simulation of shrinkage related defects in metal castings—A review [J]. Int. J. Cast Met. Res., 2005, 18: 129
doi: 10.1179/136404605225023018
[21] Carlson K D, Lin Z P, Beckermann C. Modeling the effect of finite-rate hydrogen diffusion on porosity formation in aluminum alloys [J]. Metall. Mater. Trans., 2007, 38B: 541
[22] Yao L, Cockcroft S, Zhu J D, et al. Modeling of microporosity size distribution in aluminum alloy A356 [J]. Metall. Mater. Trans., 2011, 42A: 4137
[23] Wu W, Sun D K, Dai T, et al. Modeling of dendritic growth and bubble formation [J]. Acta Phys. Sin., 2012, 61(15): 39
[23] 吴 伟, 孙东科, 戴 挺 等. 枝晶生长和气泡形成的数值模拟 [J]. 物理学报, 2012, 61(15): 39
[24] Zhu M F, Li Z Y, An D, et al. Cellular automaton modeling of microporosity formation during solidification of aluminum alloys [J]. ISIJ Int., 2014, 54: 384
doi: 10.2355/isijinternational.54.384
[25] Niyama E, Uchida T, Morikawa M, et al. A method of shrinkage prediction and its application to steel casting practice [J]. J. Jpn. Foundry. Soc., 1982, 54: 507
[26] Polyakov S, Korotchenko A, Bast J. Use of the Niyama criterion to predict porosity of the mushy zone with deformation [J]. Arch. Foundry Eng., 2011, 11: 131
[27] Liu C C, He J S. Numerical prediction of shrinkage defect in aluminum alloy electron beam weld based on Niyama criterion program [J]. Adv. Mater. Res., 2014, 1030-1032: 130
[28] Ignaszak Z. Discussion on usability of the Niyama criterion for porosity predicting in cast iron castings [J]. Arch. Foundry Eng., 2017, 17: 196
doi: 10.1515/afe-2017-0115
[29] Ivanina E S, Monastyrskiy V P, Ershov M Y. Quantitative estimation of formation of shrinkage porosity by the niyama criterion [J]. Inorg. Mater. Appl. Res., 2022, 13: 100
doi: 10.1134/S2075113322010130
[30] Carlson K D, Beckermann C. Prediction of shrinkage pore volume fraction using a dimensionless Niyama criterion [J]. Metall. Mater. Trans., 2009, 40A: 163
[31] Monroe C, Beckermann C. Prediction of hot tearing using a dimensionless Niyama criterion [J]. JOM, 2014, 66: 1439
doi: 10.1007/s11837-014-0999-7
[32] Guo J, Beckermann C, Carlson K, et al. Microporosity prediction and validation for Ni-based superalloy castings [J]. IOP Conf. Ser., 2015, 84: 012003
[33] Khalajzadeh V, Carlson K D, Backman D G, et al. A pore-centric model for combined shrinkage and gas porosity in alloy solidification [J]. Metall. Mater. Trans., 2017, 48A: 1797
[34] Zhang A, Guo Z P, Jiang B, et al. Phase-field modeling of microstructure and gas porosity evolution during solidification of alloys: A review [J]. Chin. J. Nonferrous Met., 2021, 31: 2976
[34] 张 昂, 郭志鹏, 蒋 斌 等. 合金凝固组织和气孔演变相场模拟研究进展 [J]. 中国有色金属学报, 2021, 31: 2976
[35] Backer G, Wang Q G. Microporosity simulation in aluminum castings using an integrated pore growth and interdendritic flow model [J]. Metall. Mater. Trans., 2007, 38B: 533
[36] Liu Y Q. Investigation on solidification process of aluminum alloys and casting technologies of complicated structural components [D]. Xi'an: Northwestern Polytechnical University, 2015
[36] 刘永勤. 铝合金凝固过程及其结构件的铸造工艺研究 [D]. 西安: 西北工业大学, 2015
[37] Lee P D, Hunt J D. Hydrogen porosity in directionally solidified aluminium-copper alloys: A mathematical model [J]. Acta Mater., 2001, 49: 1383
doi: 10.1016/S1359-6454(01)00043-X
[38] Huang J L, Conley J G, Mori T. Simulation of microporosity formation in modified and unmodified A356 alloy castings [J]. Metall. Mater. Trans., 1998, 29B: 1249
[39] Gui Z L, Dai T, Zhu M F. Present status of numerical simulation of micro-porosity in aluminum alloys [J]. Spec. Cast. Nonferrous Alloys, 2007, 10: 766
[39] 桂仲林, 戴 挺, 朱鸣芳. 铝合金显微疏松数值模拟的研究现状 [J]. 特种铸造及有色合金, 2007, 10: 766
[40] Su M, Sun D K, Zhu M F. Application of lattice boltzmann method to dendrite growth simulation [J]. Mater Mec. Eng., 2007, 31(2): 75
[40] 苏 敏, 孙东科, 朱鸣芳. Lattice Boltzmann方法在枝晶生长模拟中的应用 [J]. 机械工程材料, 2007, 31(2): 75
[41] Li Z Y, Zhu M F, Dai T. Modeling of microporosity formation in an Al-7%Si Alloy [J]. Acta Metall. Sin., 2013, 49: 1032
doi: 10.3724/SP.J.1037.2013.00128
[41] 李正扬, 朱鸣芳, 戴 挺. Al-7%Si合金显微气孔形成的模拟研究 [J]. 金属学报, 2013, 49: 1032
doi: 10.3724/SP.J.1037.2013.00128
[42] Meidani H, Jacot A. Phase-field simulation of micropores constrained by the dendritic network during solidification [J]. Acta Mater., 2011, 59: 3032
doi: 10.1016/j.actamat.2011.01.043
[43] Zhang A, Guo Z P, Jiang B, et al. Multiphase and multiphysics modeling of dendrite growth and gas porosity evolution during solidification [J]. Acta Mater., 2021, 214: 117005
doi: 10.1016/j.actamat.2021.117005
[44] Du L F, Wang L L, Zheng B, et al. Phase-field modeling of microstructure evolution during solidification in presence of gas bubble [J]. Comput. Mater. Sci., 2016, 114: 94
doi: 10.1016/j.commatsci.2015.12.018
[45] Liu L, Pian S, Zhang Z, et al. A cellular automaton-lattice Boltzmann method for modeling growth and settlement of the dendrites for Al-4.7%Cu solidification [J]. Comput. Mater. Sci., 2018, 146: 9
doi: 10.1016/j.commatsci.2018.01.015
[46] Wang L. Simulation of flow around cylinder with lattice Boltzmann method [J]. Acta Sci. Nat. Univ. Pekin., 2002, 38: 647
[46] 王 龙. 圆柱绕流的LBM模拟 [J]. 北京大学学报(自然科学版), 2002, 38: 647
[47] Zhu M F, Stefanescu D M. Virtual front tracking model for the quantitative modeling of dendritic growth in solidification of alloys [J]. Acta Mater., 2007, 55: 1741
doi: 10.1016/j.actamat.2006.10.037
[48] Jelinek B, Eshraghi M, Felicelli S, et al. Large-scale parallel lattice Boltzmann-cellular automaton model of two-dimensional dendritic growth [J]. Comput. Phys. Commun., 2014, 185: 939
doi: 10.1016/j.cpc.2013.09.013
[1] 杨林洁, 张丽丽, 赵九洲, 江鸿翔, 何杰. Zn-Pb偏晶合金微量界面活性元素选取[J]. 金属学报, 2026, 62(3): 532-540.
[2] 贾玉亮, 张勇佳, 史泽楷, 沈旭, 殷亚军, 施长坤, 周建新, 吕志刚. 高温合金定向凝固雀斑形成机理及控制技术研究进展[J]. 金属学报, 2026, 62(2): 309-327.
[3] 王飞翔, 陈忠奉, 尹晓宇, 熊良华, 谢红兰, 邓彪, 肖体乔. 基于X射线体视成像实现高温合金熔体凝固三维显微结构的原位观测[J]. 金属学报, 2025, 61(7): 1109-1118.
[4] 王强, 李小兵, 郝俊杰, 陈波, 张滨, 张二林, 刘奎. 一种新型Ti-Al-Mn-Nb合金的固态相变行为[J]. 金属学报, 2025, 61(7): 1060-1070.
[5] 赵广迪, 李阳, 姚晓雨, 王亮, 李渭滨, 潘玉华, 李维娟, 王兆宇. BFe-Cr-B-C合金凝固行为、强韧性及耐磨性的影响[J]. 金属学报, 2025, 61(5): 699-716.
[6] 杨明辉, 李星吾, 孙崇昊, 阮莹. 定向凝固与固态相变双联协控下Monel K-500合金的组织和力学性能[J]. 金属学报, 2025, 61(4): 561-571.
[7] 梁炫, 侯廷平, 张东, 谭昕暘, 吴开明. 中碳Nb合金化钢液析碳化物的析出行为[J]. 金属学报, 2025, 61(4): 653-664.
[8] 徐涛, 邓安元, 李阳, 王恩刚. 连铸结晶器内Ar气泡破碎-聚合和捕获行为的数值模拟[J]. 金属学报, 2025, 61(12): 1895-1910.
[9] 孙昊, 江鸿翔, 赵九洲, 张丽丽, 何杰. 微重力条件下Gd-Co-Ti合金原位粒子复合凝固组织的形成[J]. 金属学报, 2025, 61(12): 1933-1944.
[10] 吴宇轩, 唐子渊, 张保泽, 郭晓玉, 罗颖, 刘铁, 王强. 强磁场下Al-Si共晶合金定向凝固组织演变及溶质迁移行为[J]. 金属学报, 2025, 61(11): 1615-1624.
[11] 林美, 郭博静, 王志军, 李俊杰, 王雷, 王锦程, 何峰. 选区激光熔化Ni58Cr23Fe10W5Ti2Ta1Nb1 多主元合金的裂纹形成机理及抑制[J]. 金属学报, 2025, 61(11): 1703-1714.
[12] 李俊杰, 李盼悦, 黄立清, 郭杰, 吴京洋, 樊凯, 王锦程. 真空自耗电弧熔炼铸锭凝固行为多尺度模拟研究进展[J]. 金属学报, 2025, 61(1): 12-28.
[13] 杨林, 马长松, 刘连杰, 李金富. 过冷(Fe1 -x Co x)79.3B20.7 合金的凝固[J]. 金属学报, 2025, 61(1): 99-108.
[14] 余东, 马威龙, 王亚莉, 王锦程. Au-Pt合金凝固-固态相变微观组织演化相场法模拟[J]. 金属学报, 2025, 61(1): 109-116.
[15] 王叶青, 付珂, 赵永柱, 苏礼季, 陈正. Fe7(CoNiMn)80B13 共晶高熵合金的深过冷非平衡凝固行为及微观组织演变[J]. 金属学报, 2025, 61(1): 143-153.