导热问题的数值解法.pptVIP

节点温度迭代次数05.005.0020.0020.0017.759.0021.4419.6129.449.1521.6719.7139.559.2021.7419.7449.609.2221.7719.754.4非稳态导热问题的有限差分法空间离散和时间离散稳态导热：边界值问题非稳态导热：初值问题、边界值问题非稳态导热过程的有限差分法：也是将导热微分方程变成差分方程。节点温度不仅和周围节点温度有关，而且还和时间有关。将时间用等间隔划分为若干等分，则当需要确定（k＋l）时刻的温度分布时必须知道前一时间间隔k时刻的温度分布。二维非稳态导热微分方程式的有限差分表达式常物性、无热源二维非稳态导热问题的微分方程式：将求解区域离散成有限个网格，从过程开始起计算，以划分时间间隔，以k表示其序号。节点在时刻的温度可表示为：(4-23)ooo第四章导热问题的数值解法Numericalmethodforheatconduction理论分析法；(2)数值计算法；(3)实验法1求解导热问题的三种基本方法：01理论分析方法:直接对微分方程在给定的定解条件下进行积分，获得解析解(closesolution)数值计算法:把在时间和空间连续的物理量的场，用有限个离散点上的值的集合来代替，通过求解按一定方法建立起来的关于这些点上的物理量值的代数方程，从而获得离散点上被求物理量的值；并称之为数值解(numericalsolution)实验法:在传热学基本理论的指导下，采用实验的方法对所研究对象的传热过程进行实验研究，从而求得所求量的方法2.三种方法的基本求解过程:023三种方法的特点3三种方法的特点(1)分析法a能获得所研究问题的精确解，可以为实验和数值计算提供比较依据；b局限性很大，对复杂的问题无法求解；c分析解具有普遍性，各种情况的影响清晰可见。(2)数值法在很大程度上弥补了分析法的缺点，适应性强，特别对于复杂问题更显其优越性；与实验法相比成本低。(3)实验法传热学的基本研究方法，a适应性不好；b费用昂贵。数值解法：有限差分法(finite-differencemethod)、有限元法(finiteelementmethod)、有限体积法(finitevolumemethod)、边界元法(boundary-elementmethod)、离散元法(discreteelementmethod)······2区域离散化3建立节点物理量的代数方程1建立控制方程及定解条件5获得数值解并分析结果4求解代数方程组传热（物理）问题的数值求解过程将求解区域离散、以节点网格代替物体，以每个节点的温度作为未知量01在节点上用差分代替微分，将微分方程式近似地变成差分方程式——线性的代数方程组02解此代数方程组，得到节点上温度的近似值031、基本思想：4.1有限差分法的基本原理函数在点的泰勒级数展开形式为：函数、、、在点x的泰勒级数展开式分别为：2、函数f（x）在点x的导数的有限差分表达式：由式(b)得：01由式(a)与式(b)相减得：02由式(a)和式(c)消去03得：04由式(b)和式(d)消去05得：060708091011由式(a)得：由式(a)和式(b)消去得：由(e)式～(j)式分别略去以上各项得一阶、二阶导数向前、向后及中心差分公式为：、及(j)、一阶导数向前差分：一阶导数向后差分：一阶导数中心差分：二阶导数向后差分：二阶导数中心差分：二阶导数向前差分：一阶导数的差分及其误差：引入下列关系式：函数f（x）在点x的一、二阶导数的有限差分表达式分别为：图4-2有限差分表达式的几何意义式中：向前和向后差分的误差比中心差分的误差高，中心差分应用较广。4.2稳态导热问题的差分表达式1。内部节点的差分方程式物理性质参数为常数的具有内热源的二维稳态导热方程：图4-3间距为Δx、Δy的矩形网格（a）将整个区域划分步长为的矩形有限差分网格，节点