铝合金阶梯试块温度场的数值模拟.docVIP

铝合金阶梯试块温度场的数值模拟 华 中 科 技 大 学 　 熊小青 Ξ　 吴志超 　 叶升平 摘 　 要 　 探讨了通用分析软件 ANSY S 在铸件凝固温度场数值模拟上的应用 。选取阶梯试样作为标准试块 , 分别进行了 数值模拟和实际浇注试验 。 结果表明 , 通过建立合理的数理模型和设置合理的边界条件 ,ANSY S 获得的数值模拟结果和实 际浇注试验基本符合 。 关键词 :ANSYS 　 数值模拟 　 铝合金 中图分类号 :TP39　　 文献标识码 :A　 文章编号 :1001-2249(2004 02-0038-03 　　 ANSY S 软件是目前应用最为广泛的大型通用有限 元分析软件 , 可用于结构 、 力学 、 热应力等多物理场的分 析 , 能够实现结构 、 温度场 、 流场之间的耦合分析 , 并且 具有强大的前后处理功能 [1,2]。本课题探讨了 。 1　 , 用有限元分析软 件 ANSY S 选取阶梯试块进行分析的目 的是为了能够直观地看到铸件不同厚度部位的凝固情 况 , 以及判断缩松缩孔可能出现的位置 。 由于铸型及铸 件在结构上是对称的 , 因此其所受载荷也是对称的 , 利 用其对称性 , 建立了包括铸型 、 铸件和浇道砂型的统一 对称有限元数值分析模型 。 1. 1　 控制方程 在铸件和铸型的不稳定导热过程中 , 温度与时间和 空间的关系可用傅立叶 (F ourier 导热微分方程来描述 , 在直角坐标系中 , 模型的温度场函数 T (x , y , z , t 应满 足以下导热微分方程 。 ρc t =x (K x x +y (K y y +z (K z z +ρQ (1 式中 , 温度 T 为坐标 (x , y , z 及时间 t 的函数 ; ρ为材 料密度 ; c 为材料的比热容 ; K x , K y , K z 为材料沿 x , y , z 方向的传热系数 (在此 K x =K y =K z ; Q 为单位体积的 内部生成热量 。 1. 2　 铸件 2铸型界面上的传热边界条件 在许多情况下 , 界面传热往往是整个系统换热的关 键 , 在砂型铸造铸件生产过程中 , 由于在砂型内表面涂 刷了一层涂料 , 因此在铸件和铸型的界面上形成热阻 , 降低了铸型的内表面在金属液接触下的温度值 。在模 型的建立和计算中 , 在铸型的界面上建立了气体的导热 模型 。 界面模型的当量导热系数为 : =g (2 , k d k , h g 。 1. 体系的对流边界条件主要包括铸型与外界空气的 自然对流换热边界和阶梯试块与空气的对流换热边界。 根据相似准则 , 采用 Nu , R e , P r 等无量纲数表示求对流 换热系数 a 的整理式 , 即 a = Nuk L (3 式中 , N u 为努塞尔数 , K f 为流体导热系数 , L 为特征长 度 (此时 L 取砂型和金属型对流表面的周长 1. 4　 凝固潜热的处理 处理潜热的方法有温度回升法 (或称温度补偿法 、 有效比热法 (又称当量比热法或等价比热法 、 比热焓 法 、 假想热流法等 。 考虑到比热焓在克服控制方程高度 非线性方面的优越性 , 本文凝固潜热采用比热焓进行处 理 。 比热焓 H 定义为 : H =∫ T T 0(c -Q i T d T (4 式中 , c 是比热容 , ξi 为新相生成量 , Q i 为 ξi 相全转变 发生时的潜热量 , T 和 T 0为温度区间的值。 2　 ANSY S 模拟过程 2. 1　 前处理及 ANSY S 求解 有限元模型是真实系统的数值模型 。确定单元类 型 、 材料特性 、 几何模型以及网格划分都是温度场分析 的关键环节 。 (1 有限元模型选择单元类型 , 最终可把选用单元 的范围缩小为 ANSY S 单元类型库中的 8节点三维热实 体单元 s olid70和 20节点的三维热实体单元 s olid90, 在 本例中选用的单元是 s olid90。 83 计算机应用技术 特种铸造及有色合金 　 2004年第 2期 Ξ熊小青 , 女 ,1978年出生 , 硕士 , 华中科技大学材料科学与工程学院 , 武汉 (430074 , 电话 :027E -mail :xxqnihao@vip. sina. com 收稿日期 :2003-09-15 (2 对材料的特性进行定义 , 选定有关参数。 该分析 中所需要的参数主要有导热系数 λ和比热容 c 及密度 ρ。 经过查表 [3]可知 , 砂型的密度为 1660kg/m 3, 砂型的初始温 度为 289K, 比热容为 0. 25619J/(g ? K , 导热系数为 0. 016621W/(m ? K , 铸件的密度为 2660kg/m 3。 (