第一章 有限差分法.ppt

******************************************超松弛迭代法求解具有稀疏系数矩阵的大型差分方程组，其中最优的就是超松弛迭代法（SuccessiveOverRelaxation，SOR）。为了说明SOR方法，首先介绍雅可比法和高斯-赛德尔法雅可比法（Jacobi）就是要使迭代值能精确的满足前一次各点的电位值所能表示的差分方程第30页,共42页，星期六，2024年，5月超松弛迭代法高斯-赛德尔法是雅可比法的改进方法，主要针对减少内存消耗，只需存储一组完整的数组。它采取的措施是对每一次迭代尽量采用最新计算的值来替换上一次迭代的旧值。结果收敛速度比雅可比法快一倍。第31页,共42页，星期六，2024年，5月超松弛迭代法逐次超松弛法是对高斯-赛德尔法的改进，该方法的核心是借助于一收敛因子w作用到高斯-赛德尔迭代公式。当时w＝1，就回到高斯-赛德尔法。当w2时，迭代过程变得及其不稳定。只有1w2，才能提高收敛速度。第32页,共42页，星期六，2024年，5月超松弛迭代法正方形第一类边界条件时长方形第一类边界条件时第33页,共42页，星期六，2024年，5月场强与电、磁积分量的计算通过上述差分方程组的求解，在获得场域内各结点上待求位函数后，往往还需求场中的场强分布，以及其他有关的积分特性（如磁通量和磁导、电导、电容等磁路及电路参数等）。第34页,共42页，星期六，2024年，5月场强与电、磁积分量的计算无论是静电场、恒定电流场或恒定磁场，其通量可一般地表示为所分析的静电场中的电容C、恒定电流场中的电导G或恒定磁场中的磁导等电路或磁路参数P就可按下式计算第35页,共42页，星期六，2024年，5月典型算例分析设长直接地金属槽的横截面如图所示，其侧壁与底面电位均为零，顶盖电位的相对值为10。试求槽中间电位分布第36页,共42页，星期六，2024年，5月典型算例分析※、场问题分析。直角坐标系，槽内电位函数满足Laplace方程，构成如下的第一类边值问题第37页,共42页，星期六，2024年，5月典型算例分析※、离散场域。用简洁的正方形网格对场域D各方向进行等分剖分p，q※、场域内差分格式。采用Laplace五点差分格式第38页,共42页，星期六，2024年，5月典型算例分析※、超松弛迭代计算。用超松弛迭代法计算差分方程※、边界条件。本例给定为第一类边值，边界条件的差分离散化应直接赋值方式第39页,共42页，星期六，2024年，5月典型算例分析※、初值。取零初值※、收敛条件指标。规定当网格内点的相邻两次迭代近似值的绝对误差的绝对值小于给定精度时，终止迭代。第40页,共42页，星期六，2024年，5月典型算例分析计算流程第41页,共42页，星期六，2024年，5月*感谢大家观看第42页,共42页，星期六，2024年，5月****************************************关于第一章有限差分法介绍有限差分方法是一种微分方法，自上世纪五十年代以来得到了广泛的应用，该方法概念清晰，方法简单，直观。虽然其与变分法相结合所形成的有限元法更有效，但有限差分还是以其固有特点在数值计算中有其重要地位，是应用最多的一种数值方法。为求解由偏微分方程定解问题所构造的数学模型，有限差分法是将定解区域（场区）离散化为网格离散节点的集合。并以各离散点上函数的差商来近似该点的偏导数，使待求的偏微分方程定解问题转化为一组相应的差分方程。根据差分方程组解出各离散点处的待求函数值——离散解。第2页,共42页，星期六，2024年，5月1、差分与差商用差分代替微分，是有限差分法的基本出发点。这一点由微分原理保证的，当自变量的差分趋于零时，差分变成微分第3页,共42页，星期六，2024年，5月差分与差商前向差分后向差分中心差分第4页,共42页，星期六，2024年，5月差分与差商通过泰勒公式分析上面差分精度，在点上的一阶导数的逼近度可由泰勒公式展开两式相减第5页,共42页，星期六，2024年，5月差分与差商前向、后向差分截断于，具有h的一阶精度，而中心差分法截断于，具有h的二阶精度，中心差分的精度比较高。函数f（x）的二阶导数前向差分前向差分第6页,共42