基于vb的窑墙温度场数值模拟程序的设计.doc

基于VB的窑墙温度场数值模拟程序的设计 李萍，曾令可，税安泽，金雪莉，王慧，刘艳春 （华南理工大学，广州 510640） 摘要：利用Visual Basic设计了窑墙温度场数值模拟系统，并利用所开发的系统程序进一步研究了窑墙各层材料物理性能参数及其厚度等对升温过程中窑墙内温度场、温度梯度、窑墙蓄热及散热的影响。 关键词：Visual Basic，温度场，窑墙，数值模拟 Design of Numerical Simulation Program for Thermal Field in Kiln Wall Based on Visual Basic LI Ping, ZENG Lingke, SHUI Anze, JIN Xueli, WANG Hui, LIU Yanchun (South China University of Technology, Guangzhou 510640) Abstract: Numerical simulation program for thermal field in kiln wall is developed by using Visual Basic. Further more, the influence of thermal properties and thickness of kiln wall materials on thermal field, thermal grads, heat accumulation and heat transfer is studied by using the program. Key words: Visual Basic, thermal field, kiln wall, numerical simulation 1 前言 陶瓷窑炉热工理论的研究主要涉及窑炉内速度场、温度场、浓度场以及制品本身的温度场、应力场等方面的内容，其基础理论为传热学，基本方程为能量、动量、质量三大守恒方程。这些方程一般都是非线性二阶偏微分方程，求解非常困难，须借助电子计算机才能完成。 窑炉在升温和降温过程中，窑体内部温度场是不断变化的。这种冷热循环很容易使窑体材料产生热疲劳和热损坏，从而缩短窑炉的寿命。根据窑体材料的物性参数及厚度等，利用计算机可模拟计算出窑体内部温度场。根据窑体材料所能承受的最大热应力和最大温差以及计算机模拟结果，制定出合理的升降温制度和窑墙材料的砌筑方式，为科学地进行热工操作提供依据。 2 数值计算方法 2.1 数学模型 对于窑墙的导热，由于在厚度方向上的温差远大于其他方向上的温差，可认为是发生在厚度方向上的一维、无内热源固体导热过程，导热微分方程可简化为： （） 式中：T——温度，℃； ——时间，s； ——导系数，m2/s它表征物体被加热和冷却时，物体内各部分温度趋向一致的能力 2.2 内节点方程组的建立 对于任一时刻时任一内节点，温度以表示。根据元体平衡法，可推导 （） 式中： 其中，；分别是节点1和n的导热系数；为节点i在时刻k时的比热；分别为区域i左右界面上的当量导热系数，可以用调和平均法计算： 2.3 边界节点方程的建立 内壁边界条件属于第一类边界条件，内壁温度与窑炉内温度一致。外壁边界条件属于第三类边界条件，根据元体平衡法，外壁节点方程可推导： （） 式中： ——自然对流换热系数，W/m2； ——辐射换热系数，W/m2； 为窑表面黑度，一般可取0.8C0＝5.669×10-8 W m-2 K-4 ——外壁温度，； ——外壁周围流体温度，。 节点方程组的求解法 所构成的n+1元方程组表示为： （） 此系数矩阵非零因素最多只有三个，第一行和最后一行只有两个非零因素。因此，从矩阵的第二行开始进行消元，把三元方程变为二元方程，直到最后一个方程变成一元方程，由此求得最后的一个节点的温度tn+1，然后逐个回代，直到求出第一个节点的温度t1为止，这种，称为TDMA）。 图1 温度场计算流程图 3.2 模拟程序系统运行流程图 图2 系统运行流程图 3.3程序特点 程序操作简单方便，图形显示直观易懂。 程序具有很强的灵活性和扩展性，可以方便地添加用户需要的功能。 最大可设5层不同窑墙材料。 在设置升温制度时，可以设置360个