当前课程知识点：有限元分析及应用 > 第13讲 有限元分析的应用领域引论(2) > 13.1 传热问题的有限元分析：基本原理 > Video 13.1
同学们好
首先回顾一下上一讲的主要内容
上一讲我们讲解了
结构振动有限元分析的基本原理
结构振动的有限元分析实例
弹塑性问题有限元分析的基本原理
弹塑性问题有限元分析的非线性方程求解
这一讲我们将讲解
传热问题的有限元分析的基本原理
传热问题的有限元分析实例
热应力问题的有限元分析基本原理
最后我们将介绍热应力问题的有限元分析实例
传热是日常生活和工程实际中广泛存在的自然现象
它包含这几个方面
首先是热量的输入和输出
它会引起热量的传递
进而引进温度的变化
从而产生热应力
因此进行热问题的分析一般包含两个部分
第一是通过传热分析来确定温度场
第二是在获得温度场的基础上
计算所产生的热应力
这是涡轮发动机在热工作环境下的传热分析
这是涡轮盘
这是进行传热分析得到的温度场
进行根据温度场来计算它的热应力
通过傅里叶传热定律和能量守恒定律
我们可以得到传热问题的控制方程
它的基本变量是温度场T
它是坐标点x,y,z的函数
同时也是时间的函数
在这个控制方程中间
我们可以看见
ρ是材料的密度
cT是材料的比热
分别沿x,y,z方向的热传导系数
Q是物体内部的热源强度
传热问题有三类边界条件
第一类，我们叫S1边界条件
也叫Dirichlet条件
它是给定温度的分布T上面一横
第二类边界条件我们叫S2
它是给定热流密度的Neumann条件
它是分别沿着偏T偏x
在x方向的热传导
然后在边界法线的方向余弦
三个方向加起来等于在S2边界上
给定的热流密度qf
第三类边界条件是S3
也叫对流换热的Neumann条件
它的基本方程的表达跟第二类边界条件类似
只是右端是与环境温度有一个热交换
其中，对热的换热系数是hc
整个物体的表面应该是包含S1+S2+S3
传热问题的初始条件
也就是说Initial Condition
它是在T0时刻整个物体的温度场分布
传热问题的变分提法就是在
满足三类边界条件及初始条件的许可温度场中
真实的温度场使得以下的泛函取极小值
这就是我们的泛函
这里的кx,кy,кz是热传导系数
Q是内部热源
cT是比热
它是关于温度的二阶导数
还有关于时间的一阶导数
对整个的物体求一个积分
在实际问题的处理过程中
第二类和第三类边界条件
由于比较难以事先满足
所以将这两个条件耦合进泛函里面去
所以我们要定义一个新的泛函
这个新的泛函除了刚才那几项以外
再加上给定热流密度qf
以及对流换热方面的一个表达
这么两项
传热问题我们分两类
一类就是稳态情况
就是温度不随时间变化
就是温度对时间的偏导为0
进行有限元分析首先要进行几何离散化
在离散化的基础上进行单元的温度场的函数表达
它的函数模式就是表达成节点的温度列阵
乘上前面的形状函数
那么节点的温度列阵就是
在前面我们变分表达的基础上
我们来求极值
这样我们就可以得到
当然它们都是关于温度的表达
我们叫作单元的传热矩阵
我们叫等效温度载荷矩阵
单元热传导矩阵的具体表达是这样的
当然前面要乘上一个x方向的热传导系数
同样在y,z方向作相应的表达
对整个单元作积分
由于我们考虑到S3的边界
所以这一项就是S3边界上由对流热交换条件
对传热矩阵所做的贡献
在相应的等效温度载荷矩阵里面
这项是内部热源的等效温度载荷
这部分是给定热流的等效温度载荷
这部分是对流热交换的等效温度载荷
我们可以看出，传热问题实际上是C0问题
温度场是一个标量场
因此所构造的有限元分析列式是比较简单的
对于瞬态的传热问题的分析
它的问题场对于时间是不等于0的
因此我们的节点温度是关于时间的函数
它的插值关系和前面的稳态实际上是相同的
只是节点的温度是
同样我们由传热问题的泛函极值原理
我们可以得到瞬态传热问题的有限元分析的列式
它叫单元的比热矩阵，它等于
它是节点温度与时间的导数，它等于
瞬态传热问题可以转化为一组
以时间t为独立变量的线性常微分方程组
它同样也是，在几何域上要进行离散
在时间域上也要进行离散
也就是说要把时间分成若干个步长
并且进行时间函数的插值
通常有两点插值，还有三点插值的公式
当然要进行一个循环的计算
这里面的关键是怎么取步长
它会影响到整个问题的收敛性
必须根据解的稳定性理论
来给出一个最大收敛步长的条件
这样才能得到稳定的解
并且控制好计算的误差
下面我们构造一个最简单的平面3节点三角形传热单元
3个节点i,j,m
其中在mj边界上有BC(S2)或BC(S3)
这两个热交换的传热条件
那么节点描述，同样我们有节点坐标(xp,yp)
p是分别等于i,j,m这3个节点
节点的温度是Ti,Tj,Tm
单元场的描述分别是
同样这个N是我们的形状函数矩阵
它的插值关系和我们静态问题的3节点三角形单元
的插值关系是一样的
这是Ni的关系的计算和3个节点的几何坐标相关的
我们考虑3种情况
第一个，这个单元没有传热边界
即为完全的内部单元
第二种情况是
如果该单元的jm边界为S2时
我们来推导这个相应的表达
第三种情况
如果这个单元的jm边为第三类传热边界
也就是S3
这个时候，我们来推导相应的单元的表达
第一种情况是完全内部的单元
没有传热边界
那么单元的传热矩阵是这么两项
кx和кy，分别乘上相应的3X3的矩阵
相应的等效温度载荷矩阵
它在Ti,Tj,Tm三个自由度上的等效温度载荷分别是
Q是内部热源的强度
A是单元的面积
ρ是材料的密度
对于jm边为S2的边界的时候
得到的单元传热矩阵前面两项也是一样
得到的等效温度载荷矩阵和前面类似
但是是在Tj,Tm两个自由度上面加上两项
因为它主要是影响到jm这条边的这两个节点
所以要加上这个边界上由热流密度引起的等效载荷
对于jm为传热边界S3的时候
单元的传热矩阵前两项和前面一样
第三项是由对流换热条件引起的3X3的矩阵
相应的等效温度载荷的矩阵
同样，除了内部热源强度引起的等效载荷以外
还由于j和m这两个节点上引起的
由于和环境温度引起的对流热交换的
这么一个等效载荷
-有限的单元 无限的能力
--Video
-课程大纲
--课程大纲
-第一章第一节测试题
-1.1 力学的分类：质点、刚体、变形体的力学
-第一章第二节测试题
-1.2 变形体力学的要点
-第一章第三节测试题
-1.3 微分方程求解的方法
--勘误
-第一章第四节测试题
-1.4 关于函数逼近的方式
-第一章第五节测试题
-1.5 针对复杂几何域上的函数表征及逼近
-第一章第六节测试题
-1.6 有限元的核心：针对复杂几何域的分片函数逼近
-第一章第七节测试题
-1.7 有限元发展的历史和软件
-课后讨论
--讨论题
-课后练习
-第二章第一节测试题
-2.1 弹簧的力学分析原理
-第二章第二节测试题
-2.2 弹簧单元与杆单元的比较
-第二章第三节测试题
-2.3 杆单元的坐标变换
-第二章第四节测试题
-2.4 一个四杆结构的实例分析
-2.5 四杆结构的ANSYS实例分析
--ANSYS
-课后讨论
--讨论题
-第三章第一节测试题
-3.1 力学描述的基本思路及关于变形体材料的基本假设
-第三章第二节测试题
-3.2 指标记法
-第三章第三节测试题
-3.3 关于三大变量及三大方程的思路
-第三章第四节测试题
-3.4 平面问题的平衡方程构建
-第三章第五节测试题
-3.5 平面问题的几何方程构建
--勘误
-第三章第六节测试题
-3.6 平面问题的物理方程构建
-第三章第七节测试题
-3.7 两类边界条件
-课后讨论
--讨论题
-第四章第一节测试题
-4.1 几种特殊情况的讨论
--勘误
-第四章第二节测试题
-4.2 简单拉杆问题的完整弹性力学求解
-第四章第三节测试题
-4.3 平面纯弯梁的描述及求解
-第四章第四节测试题
-4.4 空间弹性问题的完整描述
--勘误
-第四章第五节测试题
-4.5 关于张量的描述及理解
-课后讨论
--讨论题
-第五章第一节测试题
-5.1 变形体力学方程求解的主要方法分类及试函数方法
-第五章第二节测试题
-5.2 平面弯曲梁求解的试函数方法-残值处理法
-第五章第三节测试题
-5.3 如何降低对试函数的高阶导数的要求
-第五章第四节测试题
-5.4 平面弯曲梁求解的虚功原理
-第五章第五节测试题
-5.5 平面弯曲梁求解的最小势能原理的变分基础
-第五章第六节测试题
-5.6 一般弹性问题的能量原理
-课后讨论
--讨论题
-第六章第一节测试题
-6.1 基于试函数的经典方法与有限元方法
-第六章第二节测试题
-6.2 有限元方法中的自然离散与逼近离散
-第六章第三节测试题
-6.3 有限元方法中的基本步骤
-第六章第四节测试题
-6.4 经典方法及有限元方法的比较
-课后讨论
--讨论题
-第七章第一节测试题
-7.1 局部坐标系中的杆单元构建及MATLAB编程
-第七章第二节测试题
-7.2 局部坐标系中的平面纯弯梁单元构建及MATLAB编程
-第七章第三节测试题
-7.3 局部坐标系中的一般梁单元构建(组装)
-第七章第四节测试题
-7.4 梁单元的坐标变换
-第七章第五节测试题
-7.5 分布力的处理
-7.6 门型框架结构的实例分析及MATLAB编程
--【知识点7.6】 门型框架结构的实例分析及MATLAB编程(PDF)
-7.7 门型框架结构的ANSYS实例分析
--【知识点7.7】 ANSYS算例-门型框架结构有限元GUI操作与命令流(PDF)
-课后讨论
--讨论题
-第八章第一节测试题
-8.1 平面3节点三角形单元及MATLAB编程
-第八章第二节测试题
-8.2 平面4节点矩形单元及MATLAB编程
-第八章第三节测试题
-8.3 轴对称单元
-第八章第四节测试题
-8.4 分布力的处理
-8.5 平面矩形薄板分析的MATLAB编程
--【知识点8.5】 平面矩形薄板分析的MATLAB编程(PDF)
-8.6 平面矩形薄板的ANSYS实例分析
--【知识点8.6】 ANSYS算例-平面矩形薄板有限元GUI操作与命令流(PDF)
-课后讨论
--讨论题
-第九章第一节测试题
-9.1 空间4节点四面体单元及MATLAB编程
-第九章第二节测试题
-9.2 空间8节点正六面体单元及MATLAB编程
-第九章第三节测试题
-9.3 参数单元的原理
-第九章第四节测试题
-9.4 数值积分
-9.5 典型空间问题的MATLAB编程
--【知识点9.5】 典型空间问题的MATLAB编程(PDF)
-9.6 典型空间问题的ANSYS分析实例
--【知识点9.6】 ANSYS算例-典型空间问题有限元GUI操作与命令流(PDF)
-课后讨论
--讨论题
-第十章第一节测试题
-10.1 节点编号与存储带宽
-第十章第二节测试题
-10.2 形状函数矩阵与刚度矩阵的性质
-第十章第三节测试题
-10.3 边界条件的处理与支反力的计算
-第十章第四节测试题
-10.4 位移函数构造与收敛性要求
-第十章第五节测试题
-10.5 C0单元与C1单元
-第十章第六节测试题
-10.6 单元的拼片试验
-第十章第七节测试题
-10.7 有限元分析数值解的精度与性质
-第十章第八节测试题
-10.8 单元应力计算结果的误差与平均处理
-第十章第九节测试题
-10.9 控制误差和提高精度的h方法和p方法
-课后讨论
--讨论题
-第十一章第一节测试题
-11.1 1D高阶单元
-第十一章第二节测试题
-11.2 2D高阶单元
-第十一章第三节测试题
-11.3 3D高阶单元
-第十一章第四节测试题
-11.4 基于薄板理论的弯曲板单元
-第十一章第五节测试题
-11.5 子结构与超级单元
-课后讨论
--讨论题
-第十二章第一节测试题
-12.1 结构振动的有限元分析：基本原理
-第十二章第二节测试题
-12.2 结构振动的有限元分析实例
-第十二章第三节测试题
-12.3 弹塑性问题的有限元分析：基本原理
-第十二章第四节测试题
-12.4 弹塑性问题的有限元分析：非线性方程求解
-课后讨论
--讨论题
-第十三章第一节测试题
-13.1 传热问题的有限元分析：基本原理
-第十三章第二节测试题
-13.2 传热问题的有限元分析实例
-第十三章第三节测试题
-13.3 热应力问题的有限元分析：基本原理
-第十三章第四节测试题
-13.4 热应力问题的有限元分析实例
-课后讨论
--讨论题
-【基本建模Project1】2D问题：带孔平板的有限元分析
--Doc I-1
-【基本建模Project2】3D问题：花型卡盘网格划分的控制
--Doc I-2
-【应用建模Project3】振动模态分析：斜拉桥的模态分析
--Doc I-3
-【应用建模Project4】弹塑性分析：厚壁圆筒受内压的弹塑性分析
--Doc I-4
-【应用建模Project5】传热分析：钢制圆柱冷却过程温度场的瞬态问题
--Doc I-5
-【应用建模Project6】热应力分析：桁架结构的温度及装配应力分析
--Doc I-6
-【高级建模Project7】结构的概率：大型液压机机架的概率设计分析
--Doc I-7
-【高级建模Project8】p方法的建模与应用：平面问题的p型单元建模与分析
--Doc I-8
