ElectromagneticFieldFEA说课.ppt

* 电工电磁场理论与应用 i(t) 电路/磁路分析方法的局限性 假設: 磁通是集中在磁路中, 而且在磁路中是均勻分布的. 電流在導体(線圈)中是均勻分布的--直流電阻. 實際: 磁通是呈現空間分佈的特性, 而且在磁路中也是不均勻分布的. 電流在導體中是不均勻分布的, 而且在高頻下, 存在明顯涡流效應. Flux distribution Flux density Flux distribution with gap 濾波電感周圍的磁場 轉折銅皮的直流電流分布 導体的渦流效應 电磁场基本方程难以直接求解，因此一般通过引入中间变量(位值)。 Electric scalar potential Magnetic vector potential 电磁场场控制方程边界条件： 定解问题 = 场域控制方程 + 场域边界条件 + 时域初始条件 第一类边界条件： 第二类边界条件： 第三类边界条件： 电磁场方程的边值问题与定解条件 电磁场场控制方程： 电磁场的数值分析技术 实际电气装置复杂，根本无法得到理论解 科学技术, 研究开发以及產品設計的客观需求 计算机技术的飞速发展（速度，容量） 计算机软件技术的发展（人机界面，图形处理，多媒体） 数值计算方法的发展（高效算法，矩阵存储技术，网格剖分技术） 电磁场数值分析的必要性与可能性： 电磁场数值分析： 以计算机为计算工具，以数值分析为计算方法，求解电磁场方程的数值解。 电磁场数值计算的根本就是偏微分方程的数值分析。 Electrostatic field j: Electric scalar potential Magnetostatic field A: Magnetic vector potential Eddy current field A: Magnetic vector potential Conduction field jc: Electric scalar potential Electrostatic field Magnetostatic field Conduction field i 工程电磁场问题的种类 Eddy current field 电磁场数值分析的有限元素法(FEM) y=y(x) F=F[f(x)] Function (函数) Functional Function (泛函) a b y(x) Timea-b =F[y(x)] F=F[f(x)]=min. Variational Calculus (泛函变分) 有限元素法(Finite Element Method)是求解偏微分方程的最常用的方法 泛函 偏微分方程 泛函变分问题,即泛函极值的问题 泛函变分 x y 0 e i j m e S 3) Region discrete (Meshing) 4) Potential interpolation (Linear or higher order) 5) Element assemble (Summing) 6) Extremization (Matrix Eqs.) 1) Field governing equations 2) Equivalent functional function @ 7) Solving matrix Eqs. 有限元分析的基本过程 二维拉普拉斯场的有限元分析举例 等价变分问题 将场域离散剖分为e0个单元，则泛函可用各个单元积分的总和表示 设单元剖分为三角形，则单元内的位函数可表示为 e 1 2 3 将单元三个节点上的位函数值j1, j2, j3 代入，可求得 a, b 和 c。 表示为 所以 所以 单元系数矩阵 在三角形单元下，单元系数矩阵的计算 单元系数矩阵相加（依各单元矩阵的元素在总体矩阵中的位置），形成总体系数矩阵 e 1 2 3 1 2 3 4 5 6 ① ② ③ ④ 总体系数矩阵的各个元素的一般形式为 求极值 得到有限元方程 边界条件的强加 边界条件强加方法 I 边界条件强加方法 II 1 2 3 4 5 6 1 2 3 4 5 6 有限元方程总体系数矩阵的特点 对称矩阵 矩阵主对角元素占优（是正定矩阵） 稀疏矩阵（每一行上的非零元素的个数等于与该行节点有关联的节点总数） 矩阵元素具有带宽分布的性质（带宽大小取决于节点的编号） X X O X X X X X O X X O X X X O O X X X X X O X X X X X X D