《art1浙江大学电气学院电磁场数值计算课件》由会员分享,可在线阅读,更多相关《art1浙江大学电气学院电磁场数值计算课件(91页珍藏版)》请在金锄头文库上搜索。
1、电磁场数值分析 需要课件的同学,请发邮件至:,或者下课后把你的邮件地址告诉我。 如果方便,也请把手机号告诉我(非电工理论与新技术专业同学),以便有事联系,谢谢!December 2, 20091art1浙江大学电气学院电磁场数值计算1 教材或参考书 倪光正:工程电磁场数值计算,机械工业出版社胡之光:电机电磁场的分析与计算,机械工业出版社 (修订版)谢德馨:工程电磁场数值分析与综合,机械工业出版社汤蕴璆:电机内的电磁场,科学出版社(Second Edition)电磁场数值分析December 2, 20092art1浙江大学电气学院电磁场数值计算2 课程的主要内容 电磁场正问题的分析、计算方法电
2、磁场数值分析n概述 n二维电磁场的有限元方法(线性、非线性)n三维电磁场的有限元法n前、后处理技术December 2, 20093art1浙江大学电气学院电磁场数值计算2 课程的主要内容 电磁场逆问题的分析、计算方法电磁场数值分析n概述 n全局优化算法 n表面响应模型n矢量(多目标)优化算法n鲁棒优化问题December 2, 20094art1浙江大学电气学院电磁场数值计算2 课程的主要内容 上机作业电磁场数值分析n利用有限元法计算典型工程电磁场问题n编制典型优化算法源程序December 2, 20095art1浙江大学电气学院电磁场数值计算3 总体安排 电磁场正问题 4次课 电磁场逆问
3、题 3次课上机作业 1次课 电磁场数值分析December 2, 20096art1浙江大学电气学院电磁场数值计算4 考核方法 上机作业50%+期末考试 50%电磁场数值分析n利用有限元法计算典型工程电磁场问题n编制典型优化算法源程序December 2, 20097art1浙江大学电气学院电磁场数值计算一、 概述 1.电磁场数值分析电磁场正问题数值分析对于描述电磁现象普遍规律的电磁场微分(积分)方程, 通过某种方法和运算,将原微分(积分)方程转化为离散的代数方程,然后利用高性能的计算机直接求解代数方程,得到原问题在一系列离散点上的数值解。差分/差商代替微分/微商有限求和代替积分Decembe
4、r 2, 20098art1浙江大学电气学院电磁场数值计算一、 概述 2.主要内容电磁场正问题数值分析离散化前处理Preprocessing-Preprocessor 将计算区域划分为一系列单元(节点),建立节点 方程求解求解代数方程组Solvingsolver后处理由 基本场量进一步计算电磁量和参数 Postprocessing-PostprocessorDecember 2, 20099art1浙江大学电气学院电磁场数值计算一、 概述 3.优、缺点电磁场正问题数值分析普适性强、用户拥有弹性大一个特定问题的边界条件、电气结构、激励等输入特性可以不必编入基本方程,而由用户输入现已发展到用户通过
5、图形界面(UGI)以非常直观的方式输入 December 2, 200910art1浙江大学电气学院电磁场数值计算一、 概述 3.优、缺点电磁场正问题数值分析不依赖领域知识用户不必具有高度专业化的电磁磁场理论,数学以及数值计算等方面的知识,就能应用商业软件解决具体的工程电磁场问题。December 2, 200911art1浙江大学电气学院电磁场数值计算一、 概述 3.优、缺点电磁场正问题数值分析理论上任何问题都可以解决December 2, 200912art1浙江大学电气学院电磁场数值计算一、 概述 3.优、缺点电磁场正问题数值分析受计算机硬件、软件发展水平制约从线性到非线性、从二维到三维
6、,再到电磁场逆问题和多物理耦合场的数值分析和计算December 2, 200913art1浙江大学电气学院电磁场数值计算一、 概述 3.优、缺点电磁场正问题数值分析No close-form expression输入和输出之间一般没有显式关系表达式,不利于分析参数变化对性能指标的影响 December 2, 200914art1浙江大学电气学院电磁场数值计算一、 概述 4.数学模型电磁场正问题数值分析1) Maxwells eqution积分形式全电流定律电磁感应定律December 2, 200915art1浙江大学电气学院电磁场数值计算一、 概述 4.数学模型电磁场正问题数值分析1) M
7、axwells equation积分形式磁通连续性原理/磁场中的高斯定理电场中的高斯定理December 2, 200916art1浙江大学电气学院电磁场数值计算一、 概述 4.数学模型电磁场正问题数值分析1) Maxwells eqution微分形式December 2, 200917art1浙江大学电气学院电磁场数值计算一、 概述 4.数学模型电磁场正问题数值分析1) Maxwells equation相量形式(复数形式) December 2, 200918art1浙江大学电气学院电磁场数值计算一、 概述 4.数学模型电磁场正问题数值分析1) Maxwells equationconst
8、itutive relations高频条件下上述媒质特性参数将为复数 December 2, 200919art1浙江大学电气学院电磁场数值计算一、 概述 4.数学模型电磁场正问题数值分析1) Maxwells equation不同媒质交界条件December 2, 200920art1浙江大学电气学院电磁场数值计算一、 概述 4.数学模型电磁场正问题数值分析2) 场矢量的微分方程December 2, 200921art1浙江大学电气学院电磁场数值计算一、 概述 4.数学模型电磁场正问题数值分析3) 时变电磁场中的电磁位December 2, 200922art1浙江大学电气学院电磁场数值计
9、算一、 概述 4.数学模型电磁场正问题数值分析4) 动态位的波动方程与洛仑兹规范引入洛仑兹规范 Lorentz GaugeDecember 2, 200923art1浙江大学电气学院电磁场数值计算一、 概述 4.数学模型电磁场正问题数值分析5) 电磁位的涡流方程与电导率规范时变电磁场中,如果求解区内存导电媒质,则会在其中感应涡流。在低频电磁场中(小于1010Hz),可忽略位移电流。此时Js:源电流密度,Je涡流密度December 2, 200924art1浙江大学电气学院电磁场数值计算一、 概述 4.数学模型电磁场正问题数值分析5) 电磁位的涡流方程与电导率规范电导率规范December 2
10、, 200925art1浙江大学电气学院电磁场数值计算一、 概述 4.数学模型电磁场正问题数值分析6) 静态场中的位函数及其满足的基本方程(A) 静电场(有源、无旋)引入标量电位函数 ,并令其与电场强度E的关系为December 2, 200926art1浙江大学电气学院电磁场数值计算一、 概述 4.数学模型电磁场正问题数值分析6) 静态场中的位函数及其满足的基本方程(A) 静电场(有源、无旋)无源区Poissons equationDecember 2, 200927art1浙江大学电气学院电磁场数值计算一、 概述 4.数学模型电磁场正问题数值分析6) 静态场中的位函数及其满足的基本方程(A
11、) 静电场(有源、无旋)拉普拉斯方程(Laplaces equation)December 2, 200928art1浙江大学电气学院电磁场数值计算一、 概述 4.数学模型电磁场正问题数值分析6) 静态场中的位函数及其满足的基本方程(B) 恒定磁场(有旋、无源)对于求解区域内包含源的恒定磁场问题,引入矢磁位A,其与磁通密度B的关系为 December 2, 200929art1浙江大学电气学院电磁场数值计算一、 概述 4.数学模型电磁场正问题数值分析7) 静态场中的位函数及其满足的基本方程(B) 恒定磁场(有旋、无源)引入库仑规范(Coulomb Gage)December 2, 200930
12、art1浙江大学电气学院电磁场数值计算一、 概述 4.数学模型电磁场正问题数值分析6) 静态场中的位函数及其满足的基本方程(B) 恒定磁场(有旋、无源)矢量泊松方程 矢量拉氏方程 December 2, 200931art1浙江大学电气学院电磁场数值计算一、 概述 4.数学模型电磁场正问题数值分析6) 静态场中的位函数及其满足的基本方程(B) 恒定磁场(有旋、无源)对于求解区内不含电流源的恒定磁场问题,既可引入矢磁量位函数,也可以引入标量位函数。对于标量磁位函数 December 2, 200932art1浙江大学电气学院电磁场数值计算一、 概述 4.数学模型电磁场正问题数值分析6) 静态场中
13、的位函数及其满足的基本方程(B) 恒定磁场(有旋、无源)December 2, 200933art1浙江大学电气学院电磁场数值计算一、 概述 4.数学模型电磁场正问题数值分析7) 动态电磁场中的位函数对及其方程矢量磁位与标量电位A-V-ADecember 2, 200934art1浙江大学电气学院电磁场数值计算一、 概述 4.数学模型电磁场正问题数值分析2) 动态电磁场中的位函数及其满足的基本方程在涡流区(不包含源区) 没引入库仑规范 December 2, 200935art1浙江大学电气学院电磁场数值计算一、 概述 4.数学模型电磁场正问题数值分析2) 动态电磁场中的位函数及其满足的基本方
14、程在涡流区(不包含源区) 引入库仑规范 December 2, 200936art1浙江大学电气学院电磁场数值计算一、 概述 4.数学模型电磁场正问题数值分析7) 动态电磁场中的位函数对及其方程在非涡流区(包含源区) 引入库仑规范 December 2, 200937art1浙江大学电气学院电磁场数值计算一、 概述 5.二维稳态场分析的统一数学模型电磁场正问题数值分析December 2, 200938art1浙江大学电气学院电磁场数值计算一、 概述 6.定解条件电磁场正问题数值分析1) 边界条件(A) 第一类边界条件(Dirichlet Boundary Condition)给定的是整个场域
15、边界S上的场函数值 December 2, 200939art1浙江大学电气学院电磁场数值计算一、 概述 6.定解条件电磁场正问题数值分析1) 边界条件(B) 第二类边界条件(Neumann Boundary Condition)给定的是场函数在边界S上的法向导数值 December 2, 200940art1浙江大学电气学院电磁场数值计算一、 概述 6.定解条件电磁场正问题数值分析1) 边界条件(C) 第三类边界条件(Cauchy Boundary Condition)给定的是边界S上的场函数与其法向导数的线性组合December 2, 200941art1浙江大学电气学院电磁场数值计算一、
16、 概述 6.定解条件电磁场正问题数值分析1) 边界条件(C) 第三类边界条件(Cauchy Boundary Condition)经过变换,可转化为第一类边界条件。如在表面阻抗法(Surface Impedance Method)的应用中。December 2, 200942art1浙江大学电气学院电磁场数值计算一、 概述 6.定解条件电磁场正问题数值分析2) 初始条件(Initial Condition)与时间坐标t相联系,给出初始瞬间待求场函数 u 在场域各处的值December 2, 200943art1浙江大学电气学院电磁场数值计算一、 概述 7.主要方法电磁场正问题数值分析1) 有限
17、元法Finite element method2) 有限差分法Finite difference methodDecember 2, 200944art1浙江大学电气学院电磁场数值计算一、 概述 7.主要方法电磁场正问题数值分析3) 边界单元法Boundary element method4) 棱边元法Edge element methodDecember 2, 200945art1浙江大学电气学院电磁场数值计算一、 概述 7.主要方法电磁场正问题数值分析5) 积分方程法Integral equation method6) 模拟电荷法Charge simulation methodDecemb
18、er 2, 200946art1浙江大学电气学院电磁场数值计算一、 概述 7.主要方法电磁场正问题数值分析7) 矩量法Method of moment8) 传输线法Transmission line methodDecember 2, 200947art1浙江大学电气学院电磁场数值计算一、 概述 7.主要方法电磁场正问题数值分析9) 几何光学法GO10) 无网格法Element free or meshless methodDecember 2, 200948art1浙江大学电气学院电磁场数值计算一、 概述 8.有关源的处理问题电磁场正问题数值分析尽量采用易于处理的源,如体、面源模型等。Dec
19、ember 2, 200949art1浙江大学电气学院电磁场数值计算二、 有限元法 1. 引言1) 历史60多年历史,40多年的电气工程应用历史December 2, 200950art1浙江大学电气学院电磁场数值计算二、 有限元法 1. 引言2) 主要特点 (A)离散化过程保持了明显的物理意义 变分原理描述了支配物理现象的物理学中的最小作用原理(如力学中的最小势能原理、静电学中的汤姆逊定理) December 2, 200951art1浙江大学电气学院电磁场数值计算二、 有限元法 1. 引言2) 主要特点 (A)离散化过程保持了明显的物理意义 基于问题固有的物理特性而予以离散化处理,列出计算
20、公式,当可保证方法的正确性、数值解的存在与稳定性等前提要素December 2, 200952art1浙江大学电气学院电磁场数值计算二、 有限元法 1. 引言2) 主要特点 (B)优异的解题和适应能力 适合场域边界几何形状以及媒质物理性质变异情况复杂的问题的求解、计算 不同媒质分界面上的边界条件是自动满足 December 2, 200953art1浙江大学电气学院电磁场数值计算二、 有限元法 1. 引言2) 主要特点 (B)优异的解题和适应能力 第二、三类边界条件不必作单独的处理December 2, 200954art1浙江大学电气学院电磁场数值计算二、 有限元法 1. 引言2) 主要特点
21、 (C)适宜编程 方便地编写通用计算程序,使之构成模块化的子程序集合,适应计算功能延拓的需要,从而即可构成各种高效能的计算软件包December 2, 200955art1浙江大学电气学院电磁场数值计算二、 有限元法 1. 引言2) 主要特点 (D)影响巨大、辐射效果明显 从数学理论意义上讲,有限元法作为应用数学的一个重要分支,很少有其他方法应用得这样广泛。它使微分方程的解法与理论面目一新,推动了泛函分析与计算方法的发展December 2, 200956art1浙江大学电气学院电磁场数值计算二、 有限元法 1. 引言3) 计算实例 中心思想:将微分方程代数方程计算(求解) 以矩形槽中的磁场为
22、例说明有限元法如何将微分方程转化为代数,进而求解December 2, 200957art1浙江大学电气学院电磁场数值计算二、 有限元法 1. 引言3) 计算实例 December 2, 200958art1浙江大学电气学院电磁场数值计算二、 有限元法 1. 引言3) 计算实例 齿、槽中心线为磁力线AE为一根磁力线FG为一根磁力线CB近似为磁力线求解区域为ABCDEA围成的矩形域December 2, 200959art1浙江大学电气学院电磁场数值计算二、 有限元法 1. 引言3) 计算实例 引入标量磁位u泛定方程 December 2, 200960art1浙江大学电气学院电磁场数值计算二、
23、 有限元法 1. 引言3) 计算实例 引入标量磁位u-定界条件 December 2, 200961art1浙江大学电气学院电磁场数值计算二、 有限元法 1. 引言3) 计算实例 边值问题December 2, 200962art1浙江大学电气学院电磁场数值计算二、 有限元法 1. 引言3) 计算实例 (A)等价的条件变分问题在应用有限元法时,首先将与边值问题等价位条件变分问题,或建立与其对应的余量表达式 December 2, 200963art1浙江大学电气学院电磁场数值计算二、 有限元法 1. 引言3) 计算实例 (A)等价的条件变分问题December 2, 200964art1浙江大
24、学电气学院电磁场数值计算二、 有限元法 1. 引言3) 计算实例 (A)等价的条件变分问题W(u):能量泛函W(u)=min表示通过u的选取使得能量泛函达到极小值。December 2, 200965art1浙江大学电气学院电磁场数值计算二、 有限元法 1. 引言3) 计算实例 (A)等价的条件变分问题由能量泛函取极值来确定函数u值的问题,叫做变分问题,带有条件的变分问题称为条件变分问题。 December 2, 200966art1浙江大学电气学院电磁场数值计算二、 有限元法 1. 引言3) 计算实例 (B)单元分析为方便计算,设矩形边AB、CB、CD的长度皆为单位长度1December 2
25、, 200967art1浙江大学电气学院电磁场数值计算二、 有限元法 1. 引言3) 计算实例 (B)将区域作单元剖分,并在区域中构造插值函数i)剖分采用三角形单元对求解域离散4个单元(element)6个节点(node)December 2, 200968art1浙江大学电气学院电磁场数值计算二、 有限元法 1. 引言3) 计算实例 (B)单元分析i)剖分将组成每个单元的节点依次定义为 i、j、m,并按逆时针方向编号e i j m1 1 2 32 4 3 23 4 6 34 5 3 6December 2, 200969art1浙江大学电气学院电磁场数值计算二、 有限元法 1. 引言3) 计
26、算实例 (B)单元分析ii)构造插值函数(shape function)在每个三角形单元中,近似地认为u值在3个节点的u值ui、uj、um线性变化December 2, 200970art1浙江大学电气学院电磁场数值计算二、 有限元法 1. 引言3) 计算实例 (B)单元分析ii)构造插值函数(shape function)单元面积December 2, 200971art1浙江大学电气学院电磁场数值计算二、 有限元法 1. 引言3) 计算实例 (B)单元分析ii)构造插值函数(shape function)December 2, 200972art1浙江大学电气学院电磁场数值计算二、 有限元
27、法 1. 引言3) 计算实例 (B)单元分析ii)构造插值函数(shape function)December 2, 200973art1浙江大学电气学院电磁场数值计算二、 有限元法 1. 引言3) 计算实例 (B)单元分析ii)构造插值函数(shape function)December 2, 200974art1浙江大学电气学院电磁场数值计算二、 有限元法 1. 引言3) 计算实例 (B)单元分析ii)构造插值函数(shape function)December 2, 200975art1浙江大学电气学院电磁场数值计算二、 有限元法 1. 引言3) 计算实例 (B)单元分析ii)构造插值函
28、数(shape function)December 2, 200976art1浙江大学电气学院电磁场数值计算二、 有限元法 1. 引言3) 计算实例 (B)单元分析iii)建立单元代数方程组将插值函数及其导数代入,将能量泛函极值问题转化为能量函数(选定基函数)的极值问题建立单元代数方程组December 2, 200977art1浙江大学电气学院电磁场数值计算二、 有限元法 1. 引言3) 计算实例 (B)单元分析iii)建立单元代数方程组December 2, 200978art1浙江大学电气学院电磁场数值计算二、 有限元法 1. 引言3) 计算实例 (B)单元分析iii)建立单元代数方程组
29、December 2, 200979art1浙江大学电气学院电磁场数值计算二、 有限元法 1. 引言3) 计算实例 (B)单元分析iii)建立单元代数方程组December 2, 200980art1浙江大学电气学院电磁场数值计算二、 有限元法 1. 引言3) 计算实例 (B)单元分析iii)建立单元代数方程组单元1December 2, 200981art1浙江大学电气学院电磁场数值计算二、 有限元法 1. 引言3) 计算实例 (B)单元分析iii)建立单元代数方程组单元2December 2, 200982art1浙江大学电气学院电磁场数值计算二、 有限元法 1. 引言3) 计算实例 (B
30、)单元分析iii)建立单元代数方程组单元3December 2, 200983art1浙江大学电气学院电磁场数值计算二、 有限元法 1. 引言3) 计算实例 (B)单元分析iii)建立单元代数方程组单元4December 2, 200984art1浙江大学电气学院电磁场数值计算二、 有限元法 1. 引言3) 计算实例 (C)总体合成i)总体合成December 2, 200985art1浙江大学电气学院电磁场数值计算二、 有限元法 1. 引言3) 计算实例 (C)总体合成ii)根据能量函数取极值条件得到代数方程组December 2, 200986art1浙江大学电气学院电磁场数值计算二、 有
31、限元法 1. 引言3) 计算实例 (C)总体合成iii)边界条件处理将u1=u2=u0、u4=u5=u6=0带入有December 2, 200987art1浙江大学电气学院电磁场数值计算二、 有限元法 1. 引言3) 计算实例 (C)求解December 2, 200988art1浙江大学电气学院电磁场数值计算二、 有限元法 1. 引言3) 计算实例 (D)后处理计算磁场强度或者磁通密度December 2, 200989art1浙江大学电气学院电磁场数值计算二、 有限元法 1. 引言4) 有限元法求解偏微分方程定解问题的全过程 (1)建立与边值问题等价位条件变分问题,或建立与其对应的余量表达式 (2)单元分析 (3)总体合成 (4)方程求解 December 2, 200990art1浙江大学电气学院电磁场数值计算二、 有限元法 1. 引言4) 有限元法求解偏微分方程定解问题的全过程 (5)后处理 December 2, 200991art1浙江大学电气学院电磁场数值计算