第三章不稳定导热的有限元.ppt

第三章 不稳定导热的有限元解法 魏艳红 南京航空航天大学 第三章 不稳定导热的有限元解法 一、数学基础 二、有限元法的解题思想和步骤 三、二维热传导方程的有限元解法 一、数学基础 差分法：以差分代替微分，对基本方程离散，建立以节点参数为未知量的线性方程组，而求得近似解。 优点：线性方程组的计算格式比较简单 不足：差分格式大多采用正方形、矩形和正三角形 有限元法：对连续体本身进行离散，根据变分原理求解问题 优点：适合于各种复杂形状和复杂边界条件的数值计算 不足：计算过程复杂 一、数学基础1、行列式 代数余子式 一、数学基础2、矩阵 矩阵A可逆的充要条件 定理1 若方阵A可逆，则|A| ? 0 。 证明 A可逆，即有A的逆矩阵A?1，使得AA?1 = E。故|A|A?1|=|A||A?1| = |E| = 1，所以|A| ? 0。 定理2 若|A| ? 0，则方阵A可逆，且 其中A*为方阵A的伴随阵: 一、数学基础 一、数学基础 3、变分方法 研究泛函的极大值和极小值的方法 1）泛函定义 给定两点1和2，连接这两点曲线的长度: 这样就建立了一个函数关系： I=I[y(x)]，称I是y(x)的泛函。自变量是个函数，因变量是普通变量。 一、数学基础 2）、泛函和函数 一、数学基础 3）、泛函和变分 一、数学基础 4、差值函数 线性差值：求过曲线y(x)上已知点A(xi,yi)、B(xi+1,yi+1)的直线方程： 一、数学基础 5、形函数 1）一维不稳定导热 求解区间[0，L]划分为有限个互不重叠的小区间。 构造的差值函数： 一、数学基础 对于三角形单元，通常假设单元e上的温度是x,y的线性函数。 一、数学基础 一、数学基础 一、数学基础 形函数(Ni, Nj, Nm)的特点： Ni, Nj, Nm是x, y的线性函数，与插值函数具有同样的类型 Ni（xi,yi）=1 , Ni（xj,yj）= Ni（xm,ym）=0 一、数学基础 二、有限元法的解题思想和步骤 1、有限元解题思想 从数学角度讲，某一泛函取极值所需要的充要条件等价于求解相应的微分方程式加边界条件。从而可利用泛函取极值的变分计算来代替微分方程及边界条件的求解。 2、有限元解题步骤 1）找到导热微分方程对应的泛函 二、有限元法的解题思想和步骤 2、有限元解题步骤 2）单元划分 二、有限元法的解题思想和步骤 2、有限元解题步骤 3）建立温度的差值函数 对于三角形单元：T=f(Ti,Tj,Tm) T=NiTi+NjTj+NmTm 4）单元变分计算 二、有限元法的解题思想和步骤 2、有限元解题步骤 5) 总体合成 得到线性方程组。 6) 求解线性方程组 三、二维热传导方程的有限元解法 1、数学模型 无内热源、假定热物理性能为常数。 三、二维热传导方程的有限元解法 3、区域离散化 三、二维热传导方程的有限元解法 4、温度差值函数的建立 对于三角形单元 T=f(Ti,Tj,Tm) T=NiTi+NjTj+NmTm 三、二维热传导方程的有限元解法 三、二维热传导方程的有限元解法 三、二维热传导方程的有限元解法 三、二维热传导方程的有限元解法 三、二维热传导方程的有限元解法 三、二维热传导方程的有限元解法 三、二维热传导方程的有限元解法 三、二维热传导方程的有限元解法 三、二维热传导方程的有限元解法 三、二维热传导方程的有限元解法 三、二维热传导方程的有限元解法 三、二维热传导方程的有限元解法 三、二维热传导方程的有限元解法 三、二维热传导方程的有限元解法 三、二维热传导方程的有限元解法 三、二维热传导方程的有限元解法 三、二维热传导方程的有限元解法 三、二维热传导方程的有限元解法 解题方法和步骤 把传热问题转化为成变分问题； 对物体进行有限元分割，把变分问题近似的表达成线性方程组； 求解线性方程组，将所得的解作为热传导问题解的近似值。 I[T(x,y)]的变分问题则转化为多元函数求极值问题： 温度场T(x,y)离散成T1,T2……Tn等n个节点温度，则泛函I[T(x,y)]实际上是一个多元函数：I(T1,T2,……,Tn)： 将区域划分成有限个三角形单元 例如，分成E个单元，n个节点 2、泛函 对应的泛函： 目标：寻找温度场T，使?I=0，即：寻找是泛函达 到极值的函数。 将一个矩形区域，划分成多个直角三角形。 设直角边长为h，(?x =?y=h) 节点x=rh,y=sh (r, s为正整数) 此节点记为(r,s)，（相当于（x,y）点） 5、单元变分的计算 将求解区域分成有限个单元后，泛函I(T)变成各个单元内泛函的积分。 5、单元变分 （