[2018年最新整理]4-1导热问题的数值解法(课件1).pptVIP

第四章 导热问题的数值解法 数值法：在很大程度上弥补了分析法的缺点，适应性强，特别对于复杂问题更显其优越性；与实验法相比成本低。 数值求解的步骤 中心差分 例2 ：某方形物体，导热系数为常数，已知各边界温度如图所示，试用高斯-塞德尔迭代法求其内部节点1、2、3、4点的温度。 1．首先选择空间坐标间隔 Δx ，即距离步长。对二维问题一般使 Δy=Δx　。 2．对显式格式差分方程，根据方程的稳定性条件选择允许的最大时间步长 Δτ ； 在稳定性条件允许范围内，Δτ 越大，计算工作量越小，但精度较差；对一维问题，一般取　　　　 ，即可满足工程精度要求；　　　 表面温度在不同 Δτ 时的计算值 4-3 非稳态导热问题的数值解法 4-3 非稳态导热问题的数值解法 对于隐式差分方程，Δx、Δτ 可任意选取，不必进行稳定性条件校核； 4-3 非稳态导热问题的数值解法 五、节点方程组的求解 四、设立迭代初场 §4-2 稳态导热问题的数值解法 五、求解代数方程组 写出所有内节点和边界节点的温度差分方程式： 直接解法 迭代解法 直接解法：矩阵求逆、高斯消元法等 缺点：所需内存较大、方程数目多时不便、不适用于非线性问题（若物性为温度的函数，节点温度差分方程中的导热系数不再是常数，而是温度的函数。这些系数在计算过程中要相应地不断更新） §4-2 稳态导热问题的数值解法 五、求解代数方程组 迭代解法：简单迭代（Jacobi迭代）、高斯-赛德尔（Gauss-Seidel）迭代、块迭代、交替方向迭代等 先对要计算的场作出假设（给定初始值）、在迭代计算过程中不断予以改进、直到计算结果与假定值的结果相差小于允许值。称迭代计算已经收敛。 Gauss-Seidel迭代法：每次迭代时总是使用节点温度的最新值 §4-2 稳态导热问题的数值解法 五、求解代数方程组 在计算后面的节点温度时应采用最新值： 根据第 k 次迭代的数值： Gauss-Seidel迭代法 判断迭代是否收敛的准则： or or ε 为允许的偏差，一般取10-3～10-6 为k次迭代得到的计算域温度最大值 计算域温度存在近于0的值时采用 §4-2 稳态导热问题的数值解法 五、求解代数方程组 §4-2 稳态导热问题的数值解法 五、求解代数方程组 如何避免迭代发散？ 必须满足对角占优原则：每个迭代变量的系数总大于/等于该式中其它变量系数绝对值的代数和 (参考教材例题4-1) 该条件可表示为： 例4-1：利用高斯-赛德尔迭代法求解下列方程： 先将上式改写成迭代形式： 假设t1、t2、t3的初始值均取零，迭代值为： 如此经过七次迭代后，在四位有效数字内得到了与精确解一致的结果。 0 0 0 3.625 5.675 3.769 1.735 4.545 4.996 1.864 4.038 5.058 1.983 3.980 5.013 2.003 3.994 5.000 2.001 4.000 5.000 2.000 4.000 5.000 0 1 2 3 4 5 6 7 t1 t2 t3 迭代次数 假设t1、t2、t3的初始值均取零，迭代值为（经三次迭代）： 0 0 0 32 -36 -155 522 -396 -3355 8722 -3996 -61755 0 1 2 3 t1 t2 t3 迭代次数 200℃ 600℃ 400℃ 800℃ 2 3 4 解：1）列节点方程（内节点） 1： 2： 3： 4： 200℃ 600℃ 400℃ 800℃ 2 3 4 设: t10=500℃,t20=650℃, t30=650℃,t40=750℃ 500 650