工程电磁场第一章3概要.ppt

收敛速度与? 有明显关系： 收敛因子（ a ） 1.0 1.7 1.8 1.83 1.85 1.87 1.9 2.0 迭代次数（ N） 1000 269 174 143 122 133 171 发散 最佳收敛因子的经验公式（不唯一） （正方形场域、正方形网格） （矩形场域、正方形网格） 收敛速度与电位初始值、网格剖分粗细有关； 迭代次数与计算精度 ? 及收敛因子有关。 下 页 上 页 泊松方程的超松弛迭代格式 下 页 上 页 4. 边界条件离散化（Discrete Boundary Condition) 第一类边界条件 网格结点与边界重合 对称边界条件 3 1 2 0 4 对称边 下 页 上 页 网格结点与边界不重合 3 1 2 0 4 h1 h2 h 同理 设：h1=ph，h2=qh 第二类边界条件 分界面衔接条件 其中 介质分界面 下 页 上 页 3 1 2 0 4 h 边界节点赋已知电位值 赋节点电位初始值 累计迭代次数 N=0 N=N+1 按超松弛法进行一次迭代，求 打印 N Y 程序框图 下 页 上 页 上机作业要求： 试用超松弛迭代法求解接地金属槽内电位的分布。 给定边值：如图示； 给定初值： 误差范围： 下 页 上 页 接地金属槽半场域的网格剖分 按对称场差分格式求解电位的分布 或： 计算：1) 迭代次数 N = ? , 分布； 2) 按电位差 画出槽中等位线。 下 页 上 页 下 页 上 页 下 页 上 页 仿真所用模型的立体视图 下 页 上 页 面板电流密度分布 下 页 上 页 面板电流密度分布 50Hz电流走向和密度分布 下 页 上 页 面板电流密度分布 50Hz电流走向和密度分布 线脚电流强度 下 页 上 页 面板电流密度分布 50Hz电流走向和密度分布 线脚电流强度 面板磁场强度分布 50Hz 下 页 上 页 下 页 上 页 交流电流在触头上的分布 下 页 上 页 交流电流产生的纵向磁场在触头表面分布 下 页 上 页 交流电流产生的纵向磁场在触头间隙中心平面上分布 下 页 上 页 3.选做题 已知：无限长矩形屏蔽空腔中长直矩形导体的横截面如图示，且给定参数为 图1.6.8 无限长矩形屏蔽空 腔中长直矩形导体的横截面 要求 用超松弛选代法求解无限长矩形屏蔽空腔 中长直矩形导体周围的电位分布; 画出屏蔽腔中矩形导体周围等位线分布; 下 页 上 页 下 页 上 页 ? 磁轭 样品 d l1 a l2 磁轭 h US + - ? ，? a + - ? ，? US 最终解 下 页 上 页 均匀电场 x 均匀外电场中介质圆柱内外的电场 下 页 上 页 介质柱内电场均匀，并与外加电场E 平行； 表明 若?2?1 ，E2 E1； 若?2?1 ，E2 E1； 取圆柱坐标系，边值问题 长圆柱导体壳上部电压为-U下部电压为U，求柱内的电位分布。 下 页 上 页 例 解 ? -U U a 通解 因 奇函数 参考点 下 页 上 页 ? -U U a 等电位面 应用傅立叶级数展开： 下 页 上 页 下 页 上 页 1.6 有限差分法 Finite Difference Method 1. 数值计算的基本思想 下 页 上 页 当计算场域的边界几何形状复杂时，应用解析法分析较困难，这时可以采用数值计算的方法。 将电磁场连续域内的问题变换为离散系统的问题求解，用离散点的数值解逼近连续域内的真实解。 把求解连续函数的偏微分方程问题转换为求解离散点上的代数方程组的问题。 下 页 上 页 下 页 上 页 绝缘子的破坏情况 下 页 上 页 复合悬式绝缘子结构 下 页 上 页 下 页 不同憎水性等级的计算模型 RHC1 RHC2 RHC3 RHC4 RHC5 RHC1伞裙上表面电场分布 下 页 上 页 下 页 上 页 下 页 上 页 2. 二维差分方程的建立 下 页 上 页 有限差分的网格分割 场域的离散 不同的离散方式得到不同的差分方程 结点多，步长小，计算结果精确 网格法 h 注意 下 页 上 页 下 页 上 页 下 页 上 页 下 页 上 页 用差分代替微分 增量 一阶差分 一阶差商 二阶差商 中心差分 二维静电场边值问题 下 页 上 页 有限差分的网格分割 当场域中 或 或 下 页 上 页 用差分方程代替微分方程 五点差分格式 若场域离散为