电磁场5分离变量法、有限差分法、有限元法简析.pptVIP

（4）求解方程 思考：（1）有限元的解跟有限差分法的解有何根本不同？ （2）有限元的系数阵总是对称的吗？ x 0 0 0 0.2 0.0361 0.0360 0.4 0.0628 0.0625 0.6 0.0710 0.0708 0.8 0.0525 0.0523 1.0 0 0 与有限差分法（FDM）相比，有限差分法是对点的离散，得到一系列离散点上的解；而有限元（FEM）是对区域的离散（单元），尽管所求的是节点上的自由度，但它的解在场域中每一个点上都有定义。 所以，即是有限元节点上的解是精确的，有限元的整个解仍然是近似的。好的数据处理技术可以从该近似解中提取更精确的分析结果。 线性单元中，如果所求的自由度是电位j，单元中的电场 E是场量；节点上的 E 取邻近单元的平均。 一些补充说明：关于有限元的解 以二维静电场泊松方程的求解为例。 2. 有限元法基本原理与实施步骤：二维问题 目标：依据加权余量法，利用分域基，建立离散的代数方程组，即确定系数{Kij} 和{bi}。 场域离散 二维问题常使用三角形单元离散，便于处理复杂的场域形状，容易实现。 单元：互不重叠，覆盖全部场域；每个单元内介质是 单一、均匀的。 节点：网格的交点，待求变量的设置点。 该步骤需要记录的信息： 节点编号、节点坐标 节点属性(激励源、是否边界等) 单元编号 单元节点编号 单元介质 * 欢迎学习 《电磁场》教案 公共邮箱—文件中心—网盘： 账号：scu_yingwei@163.com 密码：scu2015142536 泊松方程 所有静电场问题的求解都可归结为在一定条件下寻求泊松方程或拉普拉斯方程的解的过程。(解二阶偏微分方程) 微 分 方 程 边 界 条 件 外边界条件 内分界条件 环路定律 高斯定律 静电场定解问题 1.4 静电场定解问题（边值问题） 静电场定解问题 静电场定解问题 1.5 分离变量法 分离变量法采用正交坐标系，将变量分离后得到微分方程的通解， 当场域边界与正交坐标面重合或平行时，才可确定积分常数，得到边值问题的解。 1.5.1 解题的一般步骤： 2）分离变量，将偏微分方程分离成几个常微分方程； 3）解常微分方程，并叠加得到通解； 1）写出边值问题（微分方程和边界条件）； 4）利用边界条件确定积分常数，最终得到电位的特解。 只含有一个变量的微分方程，采用积分法求解。含有两个变量的微分方程，可以采用分量变量法求解。 例1.5.1 试求长直接地金属槽内电位的分布。 解: 1）确定边值问题 1.5.2 应用实例 1. 直角坐标系中的分离变量法（二维场） （D 域内） 图1.5.1 接地金属槽的截面 y 2）分离变量试探解 则 -分离常数, 设 代入微分方程得 电位方程为 ——二阶常系数齐次方程 ——拉普拉斯方程 双曲函数 即 kn 为实数时， 若 ， 若 ， 若 ， 3）解常微分方程，将各特解线性叠加得通解。 4）利用给定边界条件确定积分常数，最终得到电位函数的解。 图1.5.1 接地金属槽的截面 y 通解 沿 x方向作正弦变化，知 题设 比较系数求常数 当 时， 当 时， 等式无法成立！ 若金属槽盖电位 ，再求槽内电位分布？ 通解 等式两端同乘以 ，然后从 积分 左式 当 时， 右式 ＝ 代入式（1） 代入通解 n＝奇数 图1.5.3 接地金属槽内 的等位线分布 解：1）取圆柱坐标系，边值问题 根据对称性 例1.5.2 垂直于均匀电场 E 放置 一根无限长均匀介质圆柱棒 ， 试求 圆柱内外 和 E 的分布。 均匀电场中的介质圆柱棒 自然边界条件 当 时， 当 时， 代入方程整理 分离变量,设 3）通解 拆分为两个方程 2） 根据 （自然边界条件），得 当 时， 根据 4）利用给定边界条件确定积分常数 当 时， 通解 根据 得到 比较系数 当n=1时， 当 时，An=Bn= 0， 则最终解 由分界面 的衔接条件，得 介质圆柱内外的电场 求电场强度E 1.6 有限差分法 1.6.1 二维泊松方程的差分格式 （1） 二维静电场边值问题 基本思想：将