1、电磁场数值计算上机作业报告一、 有限差分法及原理有限差分法基本思想是把连续的定解区域用有限个离散点构成的网格来代替,这些离散点称作网格的节点；把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似；把原方程和定解条件中的微商用差商来近似,积分用积分和来近似,于是原微分方程和定解条件就近似地代之以代数方程组,即有限差分方程组 ,解此方程组就可以得到原问题在离散点上的近似解.然后再利用插值方法便可以从离散解得到定解问题在整个区域上的近似解.在采用数值计算方法求解偏微分方程时,若将每一处导数由有限差分近似公式替代,从而把求解偏微分方程的问题转换成求解代数方程的问题。求解拉普拉斯方程：20
2、xy为简单起见，将场域分成足够小的正方形网格，网格线之间的距离为 h， 0。节点 0、1、2 、3 、4 上的电位分别用 0、 1、 2、 3和 4表示。点 1、点 3 在 x0 处可微，沿 x 方向在 x0 处的泰勒级数展开式为 231000011()()().!hhhxx2330000()()().!x点 2、点 4 在 y0 处可微，沿 y 方向在 y0 处的泰勒级数展开式为232000011()()().!hhy2340000()()().!yy忽略高次项 221234000()()4hxy 稍作变化得到拉普拉斯方程的五点差分格式： 12340利用超松弛迭代法求解以上差分方程，二维场拉
3、普拉斯方程等距剖分差分格式公式为：+1, =,+4(+1,+,+1+1-1,+1,-1-4,)其中 为超松驰因子， . 1 2计算流程如图 1 所示：图 1 超松驰迭代法计算流程利用上述方法求解两题。二、 第一题计算长直接地金属槽中的电场分布。金属槽横截面如图所示，其侧壁与底面电位均为零，顶盖电位相对值为 10。槽内电位函数满足拉普拉斯方程。计算槽内电位分布。要求：（1）先用正方形网格粗分，每边取 4 个网格计算，取不同的松弛因子，比较其收敛速度。取计算精度为千分之一。 （2）划分网格加倍，计算电位分布，并与上面计算结果比较。（1） 用正方形网格，每边取 4 个网格计算。图 2 题 1 示意图
4、网格分割如图 3 所示：取 =1，计算得结果如图 4：图 3 每边 4 个网格分割图 4 4 个网格分割计算结果具体数值表示结果如图 5 所示：取不同的超松驰因子 ，得不同的收敛速度如下表：超松驰因子 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9循环次数14 11 8 10 12 16 19 27 44 82可见随着 从 1 到 2 的增加，收敛速度先变快后变慢。在 =1.2 时收敛速度最快。图 5 4 个网格分割计算结果数值表示表 1 不同超松驰因子 下的迭代次数(2) 划分网格加倍，计算结果如下图所示图 6 8 个网格分割计算结果具体数值表示为：与划分为 4
5、个网格的结果做比较有下表：0.00 0.25a 0.5a 0.75a aX轴Y 轴8 格 4 格 8 格 4 格 8 格 4 格 8 格 4 格 8 格 4 格a10.00 10.00 10.00 10.00 10.00 10.00 10.00 10.00 10.00 10.00 0.75a 0.00 0.00 4.31 4.29 5.36 5.27 4.31 4.29 0.00 0.00 0.5a 0.00 0.00 1.84 1.87 2.50 2.50 1.84 1.87 0.00 0.00 0.25a 0.00 0.00 0.69 0.71 0.96 0.98 0.69 0.71 0.
6、00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 可见在拟合过程中，划分网格较少时，距离 y=a 较近处，数据偏小，而距离 y=0 较近处，数据偏大。可见划分网格较多时可以使得到的电势分布趋于平滑。表 2 相同节点位置不同网格划分计算数值对比图 7 8 个网格分割计算结果数值表示三、 第二题计算电机电枢槽气隙磁位分布：如图所示，忽略曲度效应，采用直角坐标系，计算二维恒定场。令定子和转子磁铁中 =，EF 是对称线，考虑对称性有：1|=1002|=0|=0 |=0求 G1（ABCDEF）内部和边界上全部节点磁位，计算精度取
7、 1/1000。设g=8，h=66 ，BC=6 ，DE=5 。已知磁位 也满足拉普拉斯方程，在所求区域内部，利用超松弛迭代法给出二维场拉普拉斯方程等距剖分差分格式公式为：+1, =,+4(+1,+,+1+1-1,+1,-1-4,)在 AB 处，有 ，转化为差分形式可得 ，所以在 AB 边界|=0 +1,=+1-1, 处的差分形式方程为：+1, =,+4(2+1,+,+1+1,-1-4,)图 8 题 2 示意图同理在 EF 处满足的差分形式方程为：+1, =,+4(2-1,+,+1+1,-1-4,)按正方形网格划分，网格大小为 11，根据对称性，得整个区域的磁位如图 9：在 G1 区域中，在每个节点用数值表示为（放大可看具体数值）：图 9 题 2 计算结果
