《电磁场与电磁波》课件_第3.5节.ppt

3.5有限差分法当所求问题的边界比较复杂，通常采用数值解法。目前，比较成熟的求解电磁场问题的数值解法很多，如主要有矩量法、有限差分法、有限元法、边界元法等等，采用计算机求数值解，理论上可以得到任意要求的精度。本节要点差分方程的建立简单迭代法超松弛法1.差分方程的建立1）一般来说，网格划分得愈细所能达到的精度愈高，但计算时间也愈长；2）网格的划分有不同的方法，这里仅介绍正方形网格划分。首先把求解区域划分成网格，把求解区域内连续的场分布用网格节点上的离散的数值解代替。有限差分法是把微分方程在给定点附近用差分代数方程代替而计算电位的一种近似方法。差分方程的建立（续）设二维平面场中每个正方形格子的边长为h?i-1,j和?i+1,j可以用在点(i,j)附近的泰勒级数展开为线与线的交点为节点设区域中某点(i,j)的电位为?i,j，则其上下左右四个点的电位分别为差分方程的建立（续）同理?i,j+1和?i,j-1在点(i,j)附近的泰勒级数展开为在h足够小的情况下，忽略4阶以上的高次项,将以上四式相加得泊松方程与拉普拉斯方程的差分形式设所研究区域中电荷密度为?V，点(i,j)电位满足泊松方程如果所研究的区域?V=0，则二维拉普拉斯方程的有限差分形式为在没有体电荷分布的区域，任意点的电位等于围绕它的四个点的电位的平均值。1.简单迭代法首先对待求节点设置初值；当初值给定后，利用拉普拉斯方程的有限差分形式，按一个固定的顺序（从左到右，从下到上）依次计算每点的电位，用围绕它的四个点的电位平均值作为新值，当所有的点计算完后，用它们的新值代替旧值就完成了一次迭代。然后再进行下一次迭代，直到每一点计算的新值和旧值之差小于指定的范围为止。特点是用前一次迭代得到的节点电位作为下一次迭代的初值。2.超松弛法与简单迭代法相比，超松弛法有两点重大改进：第一是计算每一节点电位时，把刚才得到的邻近点第二电位新值代入，即在计算(i,j)点的电位时，把它左边的点和下面的点的电位新值代入，即第二，再把上式写成增量的形式由于提前使用了新值，所以加快了收敛速度!增量超松弛法(续)为了加快收敛，我们引进松弛因子s，上式可以表达为松弛因子s的取值一般在1与2之间，它有一个最优值。如果松弛因子选择适当，收敛速度还将加快。当我们给予每点的增量超过使方程达到局部平衡所需要的值，将加快收敛速度。[例3-8]长方形截面的无限长导体槽，其上有一块与槽绝缘的盖板，槽的电位为零，盖板的电位为100V，求槽内的电位函数。分析：在直角坐标系中，矩形槽的电位满足拉普拉斯方程，利用Laplace方程的差分形式，采用超松弛法即可求得槽内的电位。取步长为1，将长方形截面划分成x方向网格数为16、y方向网格数为10的格子，共有16×10=160个网孔、17×11=187个节点，其中槽内节点为15×9=135个（待求），边界节点有187-135=52个（电位已知）。设迭代精度为10-6，利用MATLAB编制程序可以求得结果。导体槽内的电位分布、等位线及电力线这个问题采用简单迭代法需要的迭代次数为222，而采用超松弛法，当松弛因子取1.591时，在同样的精度下迭代次数只有40，显然超松弛法比简单迭代法收敛快。矩形槽内的电位分布是左右对称，说明计算的场域可以缩小一倍。[例3-9]长方形截面的无限长导体槽，边界条件如图所示，求槽内的电位函数。分析：在直角坐标系中，矩形槽的电位仍满足拉普拉斯方程，但这是一个含第二类边界条件的问题，仍然利用Laplace方程的差分形式，采用超松弛法求得槽内的电位，但须要保证每次迭代后保证第二类边界上各节点的取值相等。取步长为1，将长方形截面划分成x方向网格数为16、y方向网格数为10的格子，共有16×10=160个网孔、17×11=187个节点，其中槽内节点为15×9=135个（待求），边界节点有187-135=52个（电位已知）。导体槽内的电位分布、等位线及电力线矩形槽内的电位分布左右不对称，且其电位值全部大于零，这是由边界条件所决定的。有限差分法