工程有限元方法温度场有限元分析PPT(31页)_详细.ppt

稳态热传导分析有限元列式 矩阵K和P的元素分别表示如下 改写成单元集成的形式 瞬态热传导分析有限元列式 瞬态温度场与稳态温度场主要的差别是瞬态温度场的场函数温度不仅是空间域的函数，而且还是时间域的函数。但是时间和空间两种域并不耦合，因此建立有限元格式时可以采用部分离散的方法。 将空间域离散为有限个单元体，在典型单元内温度仍可以近似的用节点温度插值得到，但要注意此时节点温度是时间的函数，即 瞬态热传导分析的插值形式 瞬态热传导分析有限元列式 等效积分形式 瞬态热传导分析有限元列式 代入上式并对W内的积分进行分部积分可得 令（伽辽金形式） 瞬态热传导分析有限元列式 由离散格式 则瞬态热传导问题离散方程 C热容矩阵 平面稳态热传导问题实例 假定单元上的温度呈线性变化，单元内任一点(x, y)的温度可表示为 解方程可以得到 三角形单元 平面稳态热传导问题实例 记 三角形ijm的面积 第一行、第二行、第三行各元素的代数余子式 代入热传导矩阵，得到 平面稳态热传导问题实例 平面稳态问题实例 平面稳态热传导问题实例 平面稳态热传导问题实例 取k=2 平面稳态热传导问题实例 平面稳态热传导问题实例 代入后得到 平面稳态热传导问题实例 采用ABAQUS进行算例验证 操作演示 平面稳态热传导问题实例 操作演示 有限元方法与应用温度场有限元分析 张有为 工程力学系 内容概要 热传导问题概述 热传导问题基本方程 稳态热传导分析有限元列式 瞬态热传导分析有限元列式 平面稳态热传导分析实例 热传导问题概述 热传导问题概述 FSW演示 Laser演示 热传导问题概述 热传导基本方程 在分析工程问题时，经常要了解工件内部的温度分布情况，例如发动机的工作温度、金属工件在热处理过程中的温度变化、流体温度分布等。物体内部的温度分布取决于物体内部的热量交换，以及物体与外部介质之间的热量交换，一般认为是与时间相关的。在一般三维问题中，瞬态温度场的场变量 在直角坐标中应满足下述热传导（Fourier热传导）微分方程 导热系数 物体内部的热源密度 比热容 密度 升温需要的热量 由x, y, z方向传入的热量 内部热源产生的热量 温度场问题（热传导）基本方程 热传导基本方程 （1）给定物体边界上的温度，称为第一类边界条件 热流密度 对流换热系数 热传导问题的边界条件 域的 边界条件 在 边界上 在 边界上 在 边界上 （2）给定物体边界上的流量（输入或输出），称为第二类边界条件 （3）给定物体边界上的对流换热，称为第三类边界条件 当q=0时为绝热边界条件 热传导基本方程 当在一个方向上，例如z方向温度变化为零时，瞬态温度场微分方程就退化为二维问题的热传导微分方程 (在 内) 这时场变量不再是z的函数。场变量同时应满足的边界条件是 二维热传导问题基本方程及边界条件 热传导基本方程 对于轴对称问题，在柱坐标中场函数应满足的微分方程是 (在 内) 边界条件是 轴对称热传导问题基本方程及边界条件 热传导基本方程 瞬态热传导问题分析的初始条件 (当t=0) 瞬态热传导方程的解即为满足控制方程、初始及边界条件的场函数 ，其应是坐标和时间的函数。 初始条件及三维稳态热传导问题基本方程 如果边界上 及内部的Q不随时间变化，则经过一定时间的热交换后，物体内个点温度也将不再随时间变化，即 ，这时瞬态热传导方程就退化为稳态热传导方程。三维问题的稳态热传导方程为 稳态热传导问题基本方程 热传导基本方程 二维问题的稳态热传导方程为 轴对称问题的稳态热传导方程为 (在 内) (在 内) 二维及轴对称稳态热传导问题基本方程 稳态热传导分析有限元列式 稳态热传导问题的变分原理（考虑各向同性） 加权残量法 在 边界上 在 边界上 在 边界上 稳态热传导分析有限元列式 分部积分 稳态热传导分析有限元列式 分部积分 稳态热传导分析有限元列式 稳态热传导问题的变分原理（考虑各向同性） 由上式得出稳态热传导问题的变分原理如下 是插值函数，它是 型插值函数，它亦具有性质 稳态热传导分析有限元列式 将空间域 离散为有限个单元体，在典型单元内各点的温度 可以近似的用单元的节点温度 插值得到 温度插值 稳态热传导分析有限元列式 离散控制方程 K热传导矩阵，P温度载荷列阵