电磁搅拌作用下Al5％Cu凝固组织的数值模拟.docVIP

铸锭材料选用 Al-5%Cu 合金，此种合金是用纯铝(99.9%)和纯铜(99.9%)配制而成的，其 热物性参数如表 1 所示，并对其进行 DSC 分析，得到 DSC 曲线如图 1(b)，合金的固液相线 温度分别是 643.77℃和 664.72℃。电磁搅拌装置示意图及其几何参数如图 1(a)所示。该电磁 搅拌装置采用高磁导率的硅钢片制成搅拌器磁轭和铁芯系统，搅拌线圈采用中空的铜导线作 为绕组安装在铁芯上，形成一对极的电磁搅拌器。在搅拌过程中，电压的选择范围是 0~220 V，频率的选择范围是 0~50 HZ。为了使铸件横断面方向接近二维传热，铸模四壁材料选用 高纯石墨，铸模底部选用绝热性能较好的耐火材料。冷却过程采取空冷方式，采用 Ｋ 型热电 偶和 QuadTemp 温度记录仪连续测量中心点的温度变化，测温点布置如图 1(a)所示。　 表 1 铝合金热物性参数 Tab.1 Thermophysical properties of Al-5wt%Cu alloy 参数 参考值  　  (a) ? ? Cp ?o ? ? h 密度(Kg/m3) 热导率 W/(m·K) 比热 J/(Kg·K) 粘度系数 kg/(m·s) 相对磁导率 电阻率 ? ·m 2 (a) 电磁搅拌装置示意图  Thermocouple Graphite casting Electromagnetic stirrer 2780(Al-5%Cu合金熔体) 2100(石墨铸模) 192.5(Al-5%Cu合金熔体) 129(石墨铸模) 1086(Al-5%Cu合金熔体) 710(石墨铸模) 1.3?? 10-3 1(铜线圈) 2000(磁轭) 8(Al-5%Cu合金熔体) 1(石墨铸模) 1.5?? 10-7 (铜线圈) 0(磁轭) -7 8?? 10-6 (石墨铸模) 1000(熔体/铸模) 100(铸模/空气) (b) (b) Al-5%Cu 合金的 DSC 曲线  　 　 Fig.1 (a)Schematic diagram of EMS apparatus ; (b) DSC curve of Al-5wt%Cu alloyelectromagnetic force in melt　 2 数学模型的建立 为便于建立数学模型，假设铝合金熔体的流动是不可压缩的黏性流动过程。雷诺数小于 2100，可以采用层流模型。除熔体的黏度随温度发生变化外，假设磁轭、石墨坩埚、铜线圈 和铝合金熔体均为各向同性材料，磁导率均为常数；除焦耳热外无其他内热源，不考虑温度 变化对密度的影响，即忽略加热过程的热膨胀，密度始终设定为固定值；忽略位移电流的产 生；铝合金熔体流动对电磁场的影响不考虑，用时均电磁力代替时变电磁力与其它场耦合； 本研究主要考虑液态金属熔体的流动状态，不考虑固液两相流动的问题。为了计算方便，整 个电磁搅拌装置的计算模型可以简化成二维问题来处理。　 2.1　宏观模型的建立　 在以上假设条件下，控制铝合金熔体的流动方程是连续方程和纳维－斯托克斯方程，其 传热控制方程是能量方程。　 在熔体搅拌过程中，认为流体始终处于连续状态，其连续性方程：　 ? ? d u dt ?????0 e? A?T??T0??  (1) (2) (3) 式中，p 为压力，Pa ；ρ为铝合金熔体密度，kg/m3；μ为黏度；μ0 为黏度系数；A 为 常数，计算中取 0.5；T0 为平均形核温度。  能量方程：  ??? ? ? ?f s ?t  ???C p  dT dt  　　　　　　　　　　　　　　　　　　　　(4) 式中，T 是温度；Cp 为比热容；λ是导热系数；/ 是凝固潜热；VI是固相分数。 凝固潜热的释放是金属凝固过程区别于一般导热过程的显著特点，是联接宏微观现象、 进行宏微观耦合模拟的一个纽带，潜热的释放直接由微观范畴的晶粒形核、生长来确定，同 时枝晶生长释放的结晶潜热反馈到温度场计算中去，在本模型中对结晶潜热的处理采用温度 回升法。 2.2　微观形核生长模型的建立　 （1）形核模型　 形核模型采用的是 Rappaz 等[6,7]利用统计方法，由于模壁表面和熔体内部的形核条件有 很大差异，需要采用两种不同的形核条件分别处理模壁表面和熔体内部的形核。假设在异质 形核时，形核发生在模壁表面和熔体表面的一系列不同的形核位置上，而这些形核位置可以 由连续分布函数高斯分布来描述。对于一定的过冷度，晶粒密度由下式给出：  dn d (??)  ?  nmax 2????? ? 1??????????mn?? 2?? 2????????  (5)  式中，???mn 是平均形核过冷度；???? 过冷度的标准方差； nmax 最大晶粒密度