林贵平《传热学》第四章_静态场边值问题解法15大作业.pdfVIP

五、数值法 1、电磁场边值问题的数值解法 由于在许多实际应用中，电磁场的边界条件过于复杂而使得 无法采用精确的解析解。因此，往往采用数值解, 常用的有： 有限差分法 （Finite-diffence method ） 有限单元法 (finite-element method) 矩量法 (method of moment 简称:MOM) 每种方法都是将一个连续域分成有限个离散点，然后求解一 系列代数方程（需要计算机编程）而不是微分或积分方程。 其中有限差分法是数值解的一种较简单的解法。作为数值解 法的入门，本课程只对该内容进行介绍 2、有限差分法的基本思路 在用有限差分法求解经典问题时： 首先，需要把求解的区域划分成网格(如正方网格) ，把 求解区域内连续的场分布，用求网节点上的离散的数值接代 替。只要将网格划分得充分地细，就能达到足够的精确度。 当然，网格划分越细，精度也就越高。但同时，所花费的计 算时间也就越长。 其次，就是将静场问题的微分方程化为差分方程，然后 进行求解。 关于网格的划分，有不同的方法。下面以正方形网格 （最 简单情况下的一种处理方法）为例，说明有限差分求解静态 场问题的思路。 设X轴上邻近O点的一点的电位为Φx，用泰勒公式展开为： 故1点的电位为： 3点的电位为： 当h很小时，4阶以上的高次项都可以略去，得到： 同样地，由2点和4点的电位，可以得到： 将上两式再相加，得到： 2 2 考虑到（在二维情况下 ∂ Φ ∂ Φ ρ + − 有）泊松方程为： ∂x 2 ∂y 2 ε 0 得到： ρ 其中 F 时称为泊松方程的有限差分方程。 ε 0 ρ 0 时，就得到拉普拉斯方程的有限差分方程。即 上式表明：任一点的电位等于围绕它的四个点的电位的平均值。 当用网络将区域划分后，对于每一个网络写出类似的式子，就 得到方程数与未知电位的网格点数目相等的线形方程组，已知 的边界条件在离散化后成为边界上的节点的已知电位值。 以一顶部电位为1000V，两侧及底部的电位为零的正方形 截面的无限长金属盒为例，为说明问题方便，在正方形区 域内的水平和垂直方向各画三条平行的等间距直线，如图 所示。 从而，有： 或 或 …………………… 用矩阵可表示为： 上面的线性方程可由计算机编程求解。 在实际的工程问题中，为了满足足够的计算精度，网格的划分 是非常细的，因此即使采用计算机编程求解，其工作量也是很 大的。 3、迭代法 3.1、简单迭代法 其步骤是： （1）先对每—格点设一初值。 这个初值完全可以任意给定，但较好的估计初值、可以较快 地求得解。 右图是初步确定初值的步骤及结果 （2）按固定顺序依次计算每点新值 n+1 1( n n n n ) 即:ϕ ϕ +ϕ +ϕ +ϕ i,j 4 i+1,j i,j +1 i−1,j i,j −1 用有限差分法求解金属盒内电位初值 （3 ）迭代终止