用宏微观耦合模型模拟铝合金凝固过程 2004.docVIP

用宏微观耦合模型模拟铝合金凝固过程 3Ξ 南昌大学　许　林ΞΞ　杨湘杰　郭洪民 摘　要　用一种宏微观耦合模型对Al 22. 5S i 在金属型中的凝固情况进行了模拟。该模型结合改进的元胞自动机模型与 有限差分法, 考虑了温度场、浓度场和微观生长过程。该改进的元胞自动机模型与经典元胞自动机相比较, 不仅考虑了温度场扩散, 而且考虑了溶质在液相中的扩散、界面曲率和溶质富集对枝晶尖端过冷度的影响。在宏微观耦合模型中采用大小两套网格分别进行宏观温度场和微观组织演变的计算。在3种不同的浇注温度下得到了3种典型的柱状晶向等轴晶转变的图形。考虑边界条件后的模拟结果发现等轴晶的形核情况与理论分析有所差别。 关键词:元胞自动机　宏微观耦合　微观组织模拟 中图分类号:TG316, TP311　　文献标识码:A　文章编号:1001-2249(2004 03-0037-03 　　预测凝固过程中微观组织的演变是控制微观组织 形成的关键因素[1]。Old field 在1966年把传热和形核及晶粒生长结合起来, 创建了宏微观模型[2]。Brown 等采用M onte 2Carlo 方法研究了合金构成、模块温度及液相线斜率等因素对晶粒结构的影响[3]。进入20世纪90年代,Rappaz 等综合确定性模型和随机性模型提出了元胞自动机(C A 模型[4]。由于数值模拟的计算量很大, 所以模拟空间一直受到限制。为解决这个问题, 提出很多方法。如采用并行计算提高计算速度、合方法减小计算量[1]等一些工作[1,5]。元胞自动机[6]。 1　数学模型 1. 1　元胞自动机模型 微观组织演变层面是一个随机元胞自动机模型, 有 4个组成要素[6]:元胞几何形状; 元胞状态(在本模型中元胞有生长能力 ; 邻居结构(本模型采用M oor 邻居 ; 演变规则(晶粒的形核、长大控制元胞4种状态的转化 。 模型采用M onte Carlo 方法处理概率均匀形核, 以G auss 分布描述熔体内与模壁表面的均匀形核。由于型壁表面和熔体内部的形核条件相差很大, 所以本模型采用两种不同的数据分别处理型壁表面和熔体内部的形核, 即两种不同的高斯分布[7]。 一旦元胞形核, 该元胞就会沿着最优生长方向以一定的速度长大, 且生长速度受尖端过冷度和溶质浓度的影响。模型中尖端生长速度的计算采用Rappaz 和K urz 提出的一种简化关系式[5]。详细算法可参考文献[7]。 1. 2　溶质扩散模型 假定固液界面溶质局部平衡 C s (t =kC L (t (1 式中, C s (t 为界面固相的平衡浓度; C L (t 为液相的平衡浓度; k , , 则溶质扩散控制方t =x (D L x +y ( D L y +C L (1-k t (2 式中, D L 为溶质在液相中的扩散系数; D s 为溶质在固相中的扩散系数; f s 为固相分数。1. 3　温度场计算 宏观温度场的模拟是进行微观组织模拟的前提。模型采用单元热平衡法求解热传导微分方程, 金属和铸型间的边界条件用复合热阻法处理。二维宏观温度场传导控制方程的有限差分法模型如下。 ρ(i C (i V (i Δt = ∑4j =1 W (i , j [ T (t , j -T (t , i ](3 其中 W (i , j = h (i , j +K (i + K (j (4 式中, ρ(i 为i 单元密度; C (i 为i 单元比热容; V (i 为i 单元体积; j 为i 单元邻居; h (i , j 为界面换热系 数; K (i , K (j 为i , j 单元导热系数; D (i , j 为i 单元到j 单元界面的距离; D (j , i 为j 单元到i 单元界面的距离; S (i , j 为i , j 单元面积; L 为网格步长。 7 3计算机应用技术特种铸造及有色合金　2004年第3期 ΞΞΞ许林, 男,1980年出生, 硕士生, 南昌大学机电工程学院, 南昌(330029 , 电话:0791-8145387,E -mail :x13424@sohu. com 收稿日期:2003-12-15 江西省自然科学基金资助项目(0250006 , 江西省教育厅科技项目, 南昌大学基础理论基金资助项目(ncuj -14 2　宏微观耦合模型 模型以顺序集成为主线, 采用不同空间尺度和时间 尺度计算宏观传输现象和微观组织演变, 即双重网格划分(粗网格和细网格 和双重时间步(宏观时间步和微观时间步 。在模型中称粗网格为单元, 细网格为元胞。在宏观单元中用有限差分法计算温度场, 在微观元胞中用改进的元胞自动机处理枝晶生长过程。并借助潜热和固相分数将两个尺度的信息紧密结合起来。2. 1　宏观温度场和微观组织演变的耦合 首先把模拟空间划分成一个个宏