材料有限元分析课件.ppt

* * 第二章 传热问题的有限元法 2.1. 热传导问题的基本方程 微体升温所需热量 由 x, y, z 方向传入微体的热量 微体内热源产生的热量 (2-1)式表明：微体升温所需热量应与传入微体的热量和微体内部产生的热量平衡。 式(2-2)：第一类边界条件（强制边界条件），边界上给定温度； 式(2-3)：第二类边界条件（自然边界条件），边界上给定热流量； 当 q = 0 时，绝热边界条件。 式(2-4)：第三类边界条件（自然边界条件），边界上给定对流换热。 特例： (1) 一维传热 (2) 二维传热 (3) 轴对称传热 求解瞬态温度场问题的实质 求解初始条件 下，满足瞬态热传导方程(2-1)及边界条件的场函数 ?(?,t)。 稳态传热问题 当 时，得稳态热传导方程 2.2. 稳态热传导问题 2.2.1. 稳态热传导有限元的一般格式（以二维问题为例） (2-13)的含义：微分方程(2-12)满足边界条件的有限元解之前提是其余量的加权积分等于零。 将(2-12)代入(2-13)进行部分积分 将空间域 ? 离散为有限个单元体，在典型的单元内各点的温度可近似用单元的结点温度插值得到，即 利用伽辽金法选择权函数 在? 边界上 (2-19)即为求解n个结点温度场的方程组。写成一般式 K? = P 其中 K 热传导矩阵， ? 结点温度场阵列，P 温度载荷阵列。 各单元对热传导矩阵的贡献，用 表示。 第三类边界条件对热传导矩阵的修正，用 表示。 给定热流产生的温度载荷 热交换产生的温度载荷 单元热源产生的温度 在边界上 P16，P17 2.2.2. 稳态热传导问题举例（平面3结点三角形单元） 求解如图所示的平面（二维）温度场 返回17 由第一章的内容得知，平面3结点三角形单元中的结点插值函数为 对结点插值函数求导后代入(2-21)的第一项得热传导矩阵元素 和单元热传导矩阵 对于具有第三类边界条件（存在对流换热）的边界单元，还应考虑边界条件引起的对热传导矩阵的修正。即将插值函数Ni代入(2-21)的第三项，得 式中 L = rs 对流边界的长度。 如果单元rsp中只有rs边存在对流换热，则对该单元热传导矩阵的修正为： 假设，求解问题的热源密度Q和给定热源q均为常数，于是可得单元热源产生的温度载荷、对流换热产生的温度载荷和给定热流产生的温度载荷： 最后将(2-23)～(2-27)的结果代入 K? = P 中，便可解得用3结点三角形单元离散的平面稳态热传导问题的温度场。 2.3. 瞬态热传导问题 2.3.1. 瞬态热传导有限元的一般格式（以二维问题为例） 此时的结点温度既是空间（用坐标表示）的函数，又是时间的函数。 即： 假设：近似函数 满足 上的边界条件，则 令余量的加权积分为零 用伽辽金法选择权函数 在边界上 代入到平面问题的热传导方程中，得 矩阵K,C和P中的元素分别为： 经过上述离散处理，将原来求解平面瞬态温度场的偏微分方程转变成了求解常微分方程(2-33)的初值问题，即 返回23 *