ANSYS有限元分析软件在热分析中的应用

收稿日期:2003-12-16

张建峰(1978～　),硕士研究生;271019　山东省泰安市。

ANSYS 有限元分析软件在热分析中的应用

张建峰　王翠玲　吴玉萍(山东科技大学机电学院材料系)

顾　明

(济南钢铁集团第一炼钢厂)

摘　要　热分析是广泛应用于各个领域的一种分析工具,ANSYS 作为有限元分析软件在热分析方面具有强大的功能。本文介绍了ANSYS 热分析的基本原理、方法,综述了ANSYS 有限元软件在热分析中的应用现状,及应用ANSYS 进行热分析的发展趋势。关键词　ANSYS 有限元　热分析　应用

Application of ANSYS in H eat 2analysis

Zhang Jianfeng Wang Cuiling Wu Yuping (Shandong University of Science and Technology )

Gu Ming

(Jinan Iron and Steel Group Corporation )

Abstract Heat 2analysis is an analytical implement widely used in many areas ,and ANSYS ,FEA software ,has mighty function for heat 2analysis.Fundamental principles and methods of use are in 2troduced.And present state on application of ANSYS in heat 2analysis and the prospect are generalized in this paper.

K ey w ords ANSYS FEA heat 2analysis application

1　引言

热分析是广泛应用于各个领域的一种分析工具。国际热分析协会(简称ICTA )的命名委员会于1977年给的定义是:热分析是在程序控制温度下测量物质的物理性质与温度关系的一类技术。定义中的程序控制温度是指按某种规律加热或冷却,通常是线性升温和线性降温。在实际生产过程中,常常会遇到各种各样的热量传递问题:如计算某个系统或部件的温度分布、热量的获取或损失、热梯度、热流密度、热应力、相变等等。所涉及的部门包括:能源、化工、冶金、建筑、电子、航空航天、制冷、农业、船舶等。比如机械加工,往往需要估算和控制工件的温度场,分析不同条件下,不同材料及几何形状对温

度场变化的影响,以及防止加工过程中工件缺陷

的产生。因此,热分析在工业生产及科学研究中具有重要的作用。

早在1887年,热分析实验首先由Le Chate 2lier 开始,当时是用热曲线方法分析陶瓷材料,

热曲线是用电流计、照相底片和切光器自动记录下来。1915年,出现了另一范畴的热分析,即

热重分析(TG )。1925年,日本电器工程师Ku 2jirai 和Akahira 首次用热重分析数据预测电绝缘

性能材料的使用寿命。而第三种热分析方法热膨胀仪出现于第二次世界大战之前。所以,在那个时期,这三大分析技术已经使用起来了。但是,仪器中全自动控制和记录还谈不上。

第二次世界大战之后,特别是上世纪50年代以来,自动控制和记录的技术得到了发展。在50年代中期,日本的全自动D TA 仪器已经有3

台。以后,自动化热天平和功率补偿式DSC 亦成为市售商品。1965年,在英国的阿伯丁举行

了首届国际热分析会议。

在1992年,另一划时代改革就是把温度调节引入到DSC。突破了由于先前的分析是在线性条件下,只能够局限于热焓测定的限制,同时较好地解决了由于热容而引起的相漂移问题。

自上个世纪后半叶以来,计算机技术突飞猛进,应用计算机技术进行热分析成为热分析发展中的一个飞跃性进步。

有限元法是以电子计算机为手段的“电算”方法,它以大型问题为对象,未知数的个数可以成千上万,因而为解决复杂的力学问题提供了一个有效的工具。由于有限元法强有力的、广泛的分析功能,以及固体力学的数学物理方程与很多其它领域的相应方程可归于同一类方程,因此很自然的被推广应用于分析其它领域问题,尤其是热分析中的场问题,甚至成了这一领域主要的分析方法。应用计算机这一先进手段,以有限元理论为基础进行数值模拟,则可以提高产品加工质量,省时省力,降低成本。ANSYS作为有效的有限元分析软件,应运而生。

2　ANSYS热分析简介

1970年,John Swanson博士创建了ANSYS 公司,致力于设计分析软件的开发,ANSYS程序的第一个版本仅提供了热分析及线结构分析功能,像当时的大多数程序一样,它只是一个批处理程序,且只能在大型计算机上运行。70年代初,ANSYS程序发生了很大变化,非线性、子结构以及更多的单元类型被加入了程序。70年代末,交互方式的加入是该程序最为显著的变化,它大大地简化了模型生成和结果评价(前处理和后处理)。在进行分析之前,可用交互式图形来验证模型的几何形状、材料及边界条件;在分析完成之后计算结果的图形显示立即可用于分析检验。今天该程序的功能更加强大,使用更加便利。ANSYS作为新颖的有限元分析软件在处理热分析问题方面具有强大的功能,热分析问题一直是ANSYS重要的应用领域。而且其界面友好,易于掌握,可以随心所欲地选择图形用户界面方式(GU I)或命令流方式进行计算,同时具有强大的网格划分功能及强大的结果后处理功能。

在ANSYS公司推出的众多产品中,有五种产品能够进行热分析,包括ANSYS/Multi2 physics、ANSYS/Mechanical、ANSYS/Thermal、ANSYS/ELO TRAN、ANSYS/ED。

3　ANSYS热分析方法理论

311　有限元基础

用有限元法计算温度场并没有什么特别之处,稳态温度场与弹性力学振动问题相似,所不同的是场变量与物理参数。在弹性力学中,场变量是位移,是向量场,其物理参数是刚度系数。而在热传导问题中,场变量是温度,是标量场,其物理参数是热传导系数,从场论的角度讲,位移场、应力场也都属于场论的范畴,它们是矢量场。因此,求解弹性力学问题所用的划分单元、求解节点阵列、求单元刚度矩阵并构成总刚度矩阵、引入边界条件、求解方程组等一系列方法,都能用于求解温度场。

在ANSYS所进行的结构力学分析中,主要使用的是能量方程式,通过Rayleigh2Ritz方法,导出有限元分析的刚性矩阵。在Rayleigh2Ritz 方法中,使用能量方程式〔1〕:

U=∫v

1

2EI

(v″)2-pv dv(1)然而,并不是所有问题都适合用能量方程式来处理。对于纯量场问题,比如热分析等问题,由于微分方程式比能量方程式更容易获得,因此比较适合采用G alerkin方法,直接生成系统的刚性矩阵。方程式如下〔1〕:

U=∫l0(c i EI(N″i)2-N i P)dx(2)程序可处理热传递的三种基本类型:传导、对流和辐射。并可进行稳态和瞬态、线性和非线性分析。热分析还具有可以模拟材料固化和熔解过程的相变分析能力以及模拟热与结构应力之间的热-结构耦合分析能力。

312　ANSYS如何进行热分析

ANSYS进行热分析的基本原理是先将所处理的对象划分成有限个单元(包含若干节点),然后根据能量守恒原理求解一定边界条件和初始条件下每一节点处的热平衡方程,由此计算出各节点温度,继而进一步求解出其他相关量。

以PLAN E55单元为例,该单元为四节点四