浮法玻璃成型过程数值模拟研究进展研究.pdf

浮法玻璃成型过程数值模拟研究进展1 张勤陈泽敬胡桅林李志信 (清华大学航天航空学院北京 100084) 摘要：本文在介绍浮法玻璃成形工艺及建立数学模型基础上，分别对数值模拟技术在锡槽中 保护气体，玻璃带、锡液流动模拟的研究进展进行了概述．针对目前浮法玻璃成型数 值模拟技术存在精度低，完整性差，无法适应实际生产等缺陷，提出了完整的模拟模 型，且对成型工艺中的节能与产品质量提高进行了探讨． 关键词：浮法玻璃成型数值模拟锡槽 1 引言 浮法工艺为世界平板玻璃生产的主导工艺。而成型则是浮法技术的核心，在基本 理论研究和生产反复实践中，国外浮法技术从开始的“强冷重热拉薄法”到“徐冷拉 薄法”发展到“双向拉薄工艺”，使成型工艺不断改进【M3。 随着建筑、汽车、电子等行业对浮法玻璃日益增长的需求下，我国在2007年平板 玻璃产量达到了5-3亿重量箱，其中浮法玻璃为83％。在浮法玻璃产品产量快速增长 的同时，新品种不断增加。但与此同时，国内大多数企业仍停留在国外80年代成型工 艺水平。少数企业虽可生产品种范围扩大到1．0--25mm厚度范围，但加工参数依赖经 验调节，无法实现工艺参数重复的一致性，影响了玻璃成型质量。以直觉知识为理论 基础，以经验设计为技术手段的成型工艺，已无法适应浮法玻璃市场的发展与竞争【3棚。 数值模拟技术的发展为浮法玻璃成型提供了新的设计方法与性能预测手段，其能减少 中试，优化工艺工程，缩短开发周期，提高产品质量，产品变更时，也能缩短参数调 试时间。 20世纪80年代以来，美、日、荷等发达国家的大型电子产品企业联合一些科研 机构、高校大力开展了玻璃成型过程中物理及力学分析、温度场分析，开发了玻璃成 型有限元分析软件。国内的玻璃成型数值模拟起步晚，而且主要利用通用成型分析软 件或有限元计算软件进行玻璃局部问题进行分析计算【5】。 的拉制模型，对浮法成型的数值模拟进行了初步探索[61。但鉴于浮法玻璃的工艺复杂， 其数值模拟涉及到材料学、流变学、传热学、计算力学等多学科交叉渗透，且现有文 献集中于浮法工艺的各分体模型模拟，尚无法对浮法玻璃成型过程进行全面模拟[71。 1本文研究得到“十一五”科技支撑计划支持 一5l一 鉴于市场对高质量、新产品的迫切需求，建立完整浮法玻璃成型模型，开发适应实际 生产的浮法玻璃成型模拟软件已经称为我国浮法技术可持续发展的保障。 2浮法工艺及数学模型 浮法玻璃成型工艺过程为来自池窑的经熔化、澄清、冷却的优质玻璃液，在锡槽 中漂浮于熔融锡液表面，完成摊平、拉薄(增厚)、抛光、冷却、固化等过程，成为高 质量平板玻璃18】。熔融玻璃液从玻璃窑炉冷却部经流槽进入锡槽中，到形成玻璃带从 锡槽中出来的整个过程，必须正确控制玻璃冷却速度以适应玻璃液摊平抛光、拉薄成 型或增厚定型、冷却硬化的工艺过程【91。根据玻璃液带在锡槽中的温度可以根据成型 工艺要求加以自动调节控制，浮法玻璃工艺如图l所示【101。 玻璃成型依赖于外力、热传递与本身流变性质的影响作用。玻璃流变性质主要指 玻璃粘度、表面张力和弹性等。玻璃带在锡槽中的冷却和硬化，取决于在成型过程中 连续地同周围介质进行热传递所产生的温度场，该温度场由传热过程制约，与玻璃液 本身及其周围介质的热物理性质，即比热、导热率等有关【11。 ‰ 嗡 图1 浮法玻璃成型工艺示意图 整个浮法玻璃成型系统可以看作多个分体系统的组合，包含玻璃液、锡液、保护 气体、耐火材料和钢壳等。尽管这些分体系统相互关联，但是由于分体系统内的作用 规律相差较大，数学物理模型与计算方法也不同。其中，玻璃带温度与速度分布研究 是成型过程数值模拟的难点与重点，其数学模型如下所示。 连续方程 等+V(以；)=。 (1) 动量方程 宴：-7+土v； (2) gt Pg