单元刚度方程的建立.doc

位移法在平面稳态传热有限元计算格式中的应用 刘大全1 孙利民2 （1 郑州大学工程力学系，郑州，450001 2 郑州大学工程力学系，郑州，450001） 摘 要 本文以弹性力学平面问题中的单元节点力平衡理论为基础，避开了在热传导方程上的数学操作，直接利用稳态传热时，单元与单元之间热传导的热流量和外界热系统通过单元边界传入的热流量的在节点处的等效热量转化关系，建立了计算域上内部单元和边界单元的单元刚度方程，再应用节点热量平衡理论，组装后的总体刚度方程是一组节点热流量平衡方程，物理概念明确，推导过程简单。 关键词 平面稳态传热，热流量，等效热量，节点热量平衡 一﹑前言 弹性力学平面问题的有限元法是有限元学习中的基础内容。连续弹性体离散化后，单元与单元之间节点取为铰链连接，单元刚度方程中的单元节点力表示铰链约束中的销钉对单元节点的作用力，它也是单元与单元之间作用力在节点处等效简化的结果，刚度方程是一组节点力平衡方程。 平面传热问题有限元计算格式通常采用泛函分析法和伽辽金加权余量法两种方法建立。它们是求解满足一定边界条件的热传导方程的数值计算方法，因数理分析过程规范而得到广泛的应用。 二﹑单元刚度方程的建立 平面稳态热传导问题，整个计算域离散化为三节点三角形单元，单元与单元之间只在节点处连接。计算域分内部和边界两类单元计算，单元坐标及单元边界上外法线方向如图1所示。 图1 单元坐标及边界外法线方向 1．内部单元刚度方程 定义三角形单元温度函数用节点温度值，和的插值形式是： (1) 其中 ，， 是三角形单元的面积。 四﹑结论 与参考文献[1]采用泛函分析法和伽辽金加权余量法比较，本文推导的平面稳态热传导有限元计算格式是正确的。可以发现，借鉴位移法导出方程时，避免了应用较深奥的变分原理和伽辽金加权余量法，而直接应用了物理概念比较明确，数学形式又比较简单的热量平衡原理。同时，单元刚度方程左端的“力项”就是单元与单元之间或外界介质与单元之间热流量在节点处的等效值，只要单元的节点温度值确定后，单元与单元之间的热流量就可以求出，这是与弹性力学平面分析单元节点位移和单元与单元之间的作用力计算是一样的。 本文推导方法，同样适用于平面矩形单元和等参元的计算，对第一、二类边界条件可以采用本文第三类边界条件同样的处理方法。 参考文献 [1] David V.Hutton．Fundamentals of Finite Element Analysis．The McGraw-hill Companies，2004． [2] 孔祥谦．有限单元法在传热学中的应用．北京：科学出版社，1998． [3] 王秀喜．计算力学基础．合肥：中国科学技术大学出版社，2009． 作者简介 刘大全，女，1976年8月生，讲师，电子邮箱： ；研究领域：转子动力学及流固耦合运动分析。 项目基金 2009年度河南省高等教育改革研究项目。 参考文献格式: 1. 专著 专著是指以单行本或多卷册形式在限定期限内出版的非连续出版物，包括图书、古籍、学位论文、技术报告、会议文集、汇编、多卷书、丛书等。其著录格式为： [序号]主要责任者.题名:其他题名信息[文献类型标志(电子文献必备,其他文献任选)].其他责任者(任选).版本项.出版地:出版者,出版年:引文页码 [引用日期](联机文献必备,其他电子文献任选).获取和访问路径(联机文献必备). 示例： [1] 广西壮族自治区林业厅.广西自然保护区[M].北京：中国林业出版社，1993. 2. 期刊 [序号]析出文献主要责任者.析出文献题名[文献类型标志].刊名:其他刊名信息,年,卷(期):页码[引用日期]. 示例： [1] 张旭,张通和,易钟珍,等.采用磁过滤MEVVA源制类金刚石膜的研究[J].北京师范大学学报：自然科学版，2002，38(4)：478-481. [2] 周桂莲，许育彬，杨智全，等.认清市场形势化解“学报情结” ：我国农业学报的现状与发展趋势分析[J]. 编辑学报, 2005, 17(3)：209-211.