有限元法基础热传导和热应力讲课文档.pptVIP

有限元法基础热传导和热应力

第1页，共53页。

11传热分析与热应力

本章参考书

孔祥谦.有限单元法在传热学中的应用.第三版，科学出版社，1998

王勖成.有限单元法.清华大学出版社，2003.第12章

曾攀.有限元分析及应用.清华大学出版社，2004.第8章

2

第2页，共53页。

11.传热分析与热应力

传热是广泛存在的自然现象，只要有温度差存在，就会有热量的传递，只要有热量的输入和输出，就会引起温度的变化。

传热分析的对象是固体、液体和气体，其应用包括热量交换、化学反应、材料相变、能量转换等。

温度的变化和不均匀分布，引起结构出现应力变化，称为热应力。

当以应力分析为目的时，为确定温度场，需要对固体进行传热计算，以便确定相关的热应力。

3

第3页，共53页。

11传热分析与热应力

11.1传热问题的基本方程

固体热传导的现象

4

第4页，共53页。

11传热分析与热应力

术语和单位

在国际标准单位制中，传热分析的术语和单位

c

比热容[J/(kg·K)]

T

温度（K或°C）

Q

热流（W/m2）

h

对流换热系数[W/(m2·K)]

t

时间（s）

k

导热系数[W/(m·K)]]

质量密度（kg/m3）

q

单位体积热生成率(W/m3)

Stefan-Boltzman常数[=5.67×10-8W/(m2·K4)]

5

第5页，共53页。

11传热分析与热应力

控制方程

对于微元dxdydz，生成的热量为

微元体内的净流出热流量为

由于热量的储存使内能增加，即

由能量守恒定律，在微元内有

6

第6页，共53页。

11传热分析与热应力

由于方程不封闭，需补充条件。各向同性均匀材料的Fourier定律

于是，有

第一类边界条件：边界上已知温度T，即

第二类边界条件：边界上已知热流密度，即

第三类边界条件：已知与物体相接触的流体介质的温度和换热系数为已知，即

7

第7页，共53页。

11传热分析与热应力

11.2变分原理与有限元

瞬态热传导问题变分泛函为

稳态问题，温度不随时间变化

泛函的变分取驻值，可得控制方程和第二类和第三类边界条件

第一类边界条件应强制满足，称为本质边界条件；

第二、第三类边界条件是自然边界条件。

8

第8页，共53页。

11传热分析与热应力

有限元法

将物体离散为n个单元体，即，将单元内的温度场用节点上的温度插值，有

代入泛函，泛函称为

9

第9页，共53页。

11传热分析与热应力

将单元矩阵组合后，得到物体的总体矩阵，对泛函变分取极值，即

得有限元方程组

[KT]称为热传导矩阵，[C]为热容矩阵，{RT}为等效节点温度载荷列阵。

对稳态问题，上述方程是线性的，热传导方程是对称正定的，求解类方法似于结构分析。

对瞬态问题，需采