工程电磁场第一章3选编.pptVIP

1.5 分离变量法 Separation Variable Method 1. 分离变量法的思想 把一个包含多个自变量的函数用含单个自变量的函数的乘积表示，从而将偏微分方程分离为几个常微分方程分别求解，最后利用边界条件确定积分常数，得到级数形式的解。 分离变量法解题的一般步骤： 分离变量法采用正交坐标系，当场域边界与正交坐标面重合或平行时，分离变量法是一种有效的方法。 2. 直角坐标系中的分离变量法（二维场） 分离变量，设解答为： 代入微分方程 分离常数的取值有三种情况： 除以12 分离常数 特解1 指数函数 特解2 特解3 通解为所有特解的叠加 对于具体问题，根据边界条件确定积分常数， 积分常数的确定一般有： 比较系数法 傅立叶级数展开法 试求长直接地金属槽内电位的分布。 边值问题（D 域内） 例 解 代入边界条件，确定积分常数 比较系数 傅立叶级数 代入通解 n＝奇数 2. 圆柱坐标系中的分离变量法（二维场） 设电位只是 和 的函数，拉普拉斯方程为： 分离变量, 设 代入微分方程 分离常数 令=n2 欧拉方程 因 周期函数  = 正整数 不是周期函数 欧拉方程 取圆柱坐标系，边值问题 垂直于均匀电场E 放置一根无限长均匀介质圆柱棒，试求圆柱内外 和 E 的分布。 例 解 通解 利用给定边界条件确定积分常数 参考点 对称性 n=1 比较系数 由分界面的衔接条件 最终解 均匀电场 表明 取圆柱坐标系，边值问题 长圆柱导体壳上部电压为-U下部电压为U，求柱内的电位分布。 例 解 通解 因 奇函数 参考点 等电位面 应用傅立叶级数展开： 1.6 有限差分法 Finite Difference Method 1. 数值计算的基本思想 当计算场域的边界几何形状复杂时，应用解析法分析较困难，这时可以采用数值计算的方法。 将电磁场连续域内的问题变换为离散系统的问题求解，用离散点的数值解逼近连续域内的真实解。 把求解连续函数的偏微分方程问题转换为求解离散点上的代数方程组的问题。 绝缘子的破坏情况 复合悬式绝缘子结构 不同憎水性等级的计算模型 RHC1伞裙上表面电场分布 2. 二维差分方程的建立 场域的离散 不同的离散方式得到不同的差分方程 结点多，步长小，计算结果精确 网格法 用差分代替微分 增量 一阶差分 一阶差商 二阶差商 中心差分 二维静电场边值问题 用差分方程代替微分方程 五点差分格式 若场域离散为矩形网格，写出差分格式 例 解 结点1， x = x0 + h ，结点3 x = x0 -h 略去高次项 同理 3.差分方程组的求解 ( Solution Method ) 差分方程的特点 当步长h 减小，结点增加，方程数很大； 方程组的系数是有规律的； 各方程的项数只有5项。 采用逐次近似的方法求解，常用的方法为超松弛迭代法 1）松弛法 假定结点电位的初值，代入差分方程，计算各结点余数； 修正余数最大点的电位，减小该点余数，再重新计算各结点的余数； 重复减小最大余数的过程，直至各余数都达很小； 松弛法的步骤 为达到精度，细分网格，重复以上过程。 在接地方形导体管中有一圆形导线（很细），电压为100V，求管线间的电位分布。 例 解 对称性，只需求八分之一区域 划分网格 细分网格，重复以上步骤，提高精度。 松弛法计算简单； 不论初值如何，必收敛于最后解答； 收敛速度慢。 2）迭代法 网格编号，假定结点电位的初值，作为解的零次近似值，代入差分方程得一次近似值； 判断误差 同步迭代 结束计算 继续计算 特点 同步迭代法计算简单；但用计算机解时需要两套存储单元； 收敛速度较慢。 高斯—赛德尔迭代法（异步迭代） 特点 用计算机解时只需要一套存储单元； 迭代时有一半用了迭代的新值，收敛速度较快。 3）超松弛迭代法 a —加速收敛因子（1 a 2） 将异步迭代和松弛法结合，采用了松弛法中余数表示形式。 收敛速度与 有明显关系： 收敛因子（ 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) 第一类边界条件 网格结点与边界重合 对称边界条件 对称边 网格结点与边界不重合 同理 设：h1=ph，h2=qh 第二类边界条件 分界面衔接条件 N=N+1 N