01数值计算基础.pptVIP

电磁场数值计算 1.What? 1.1数值计算 ——自然界中的每一现象都可借助物理定律，按照与各种主要量相联系的代数方程、微分方程或积分方程来描述。 ——在科学研究及工程应用中，人们主要关心的量便是某个数学物理方程的解，包括解析解和数值解。 ——数值计算是研究各种数学问题的数值方法设计、分析、有关的数学理论和具体实现的一门学科。 1.2 电磁场理论 Maxwell方程组由4个方程构成的偏微分方程组，加上介质中的本构方程，以及两种媒质交界面的边界条件，还有Lorentz力公式，这些简洁的公式几乎可以解释所有的电磁现象，我们的任务就是在各类工程电磁问题中尽可能精确地求解这些方程： 1.3工程应用 电磁场数值计算在多个工程领域中都得到应用，例如： ——电力系统：高压(高压输电线、绝缘子)、电机、变压器、电缆等； ——电子与微波：高速PCB、波导、谐振腔、辐射、天线等； ——相关领域：感应加热、无损检测、电磁成形、电磁生物效应等； 2.How? 2.1有限元法(Finite Element Method) 现代FEM第一个成功的尝试，是将刚架位移法推广应用于弹性力学平面问题，这是Turner，Clough等人在分析飞机结构时于1956年得到的成果。 1960年Clough首次提出了“有限元法”的名称； 60年代，科学家证明了FEM是基于变分原理的Ritz法的另一种形式；并进一步利用加权余量法来确定单元特性和建立有限元方程，主要利用的是Galerkin法； 直到1968年，FEM才开始应用于电磁问题。 FEM的基本思想是分片插值，即： ——将连续的求解区域离散为一组有限个、且按一定方式相互联结在一起的单元的组合体； ——利用每个单元内假设的近似函数来分片表示全求解域上待求的未知场函数； ——单元内的近似函数通常由未知场函数及其导数在单元的各个节点的数值及其插值函数来表达； 这样未知函数从一个连续的无限自由度问题变成离散的有限自由度问题。随着单元数目的增加，单元尺寸的缩小，或单元自由度的增加及插值函数阶数的提高，近似解将收敛于精确解。 FEM相比其它数值方法的优点在于： ——理论基础成熟； ——计算格式规范统一，利于编程； ——适应性高，适合各种复杂形状的区域； ——求解精度高； 由于这些优异的特性，在短短几十年时间里，FEM成为了绝大多数物理和工程问题中(机械、航空、汽车、船舶、土木、海洋工程、电气电子、压力容器等)应用最广泛的一种计算机辅助分析方法。 在电磁分析领域，除了FEM以外，也有其它有效的数值方法，例如：矩量法(MOM)、边界元法(BEM)、时域有限差分法(FDTD)等等。 2.2 有限元软件 前处理程序 (几何建模、赋材料属性、网格剖分) 单元计算程序 (包括单元刚度矩阵的计算和总体刚度矩阵的组装) 代数方程组求解程序 (得到各个离散单元内的未知量值) 后处理程序 (通过插值得到区域每点的值，将结果数据可视化以及进一步处理) 数学基础 有限元法 可以从变分原理导出， 也可以从加权余量法导出。 前者需要补充泛函、变分法、欧拉方程、泛函 极值等数学知识，推导过程比较复杂； 后者相对比较直观，而且应用范围更广，推导 过程简单。 边界元法可以直接应用加权余量法（格林定理，积分方程） 通过了解二维静电场，三角形线性插值有限元法的基本原理，实施过程 重点掌握：加权余量概念，有限元法（和边界元法）的基本原理和实施方法。 1 加权余量法 1. 加权余量概念 假定边值问题方程 Lu=f 式中，u为未知函数。 L是算符（算子），表示对u的一种运算， f为已知函数。 为求u，设有一组完备、线性无关的函数 u1, u2, …, uk, …,（近似函数） 取其前m项的线性组合作为u的近似解ū。 若当m→∞时，有ū→ u， 则称u1, u2, …, uk, …为基函数序列， uk为基函数， c1, c2, …, cm为待定系数。 u的近似解表示为 ū= ∑mj=1 cj uj 将近似解代入方程，得余量 R=L ū - f= ∑mj=1 cj L uj – f 如果余量为零，说明已经满足方程。 即ū是方程的精确解。 一般情况下余量不为零。 只能放松约束， 强制余量的加权积分为零。 即 ∫? wi Rd ?=0 （i=1,2,…,n） 式中wi为权函数， w1，w2，…，wk，…为权函数序列， 权函数之间要求线性无关。 权函数的不同选择导致不同的近似方法。 2. 几种加权余量法 （1）配点法 在求解