欢迎来到天天文库
浏览记录
ID:9851958
大小:654.51 KB
页数:69页
时间:2018-05-11
《ansoft 高级培训班教材 ansoft hfss的有限元理论基础》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、Ansoft高级培训班教材AnsoftHFSS的有限元理论基础谢拥军编著西安电子科技大学Ansoft培训中心第一章概述第二章有限元的基本理论及三维有限元分析2.1电磁场边值问题及其变分原理2.2有限元方法的原理――从一维简单例子来看其建模过程2.3三维时谐场有限元问题2.4有限元方程组的求解第三章电磁内问题和散射问题的有限元分析方法3.1电磁内问题3.2电磁散射问题第一章概述AnsoftHFSS软件是应用有限元方法的原理来编制的,深入的了解有限元方法的理论基础,及其在电磁场与微波技术领域的应用原理,对于我们灵活、准确地使用AnsoftHFSS软件来解决实际工
2、程问题能够提供帮助。这一部分教材的内容就是在结合AnsoftHFSS软件中涉及到的有限元技术,力争在最小的篇幅和最短的时间里为学员建立理论结合实际的有限元方法的基本概念。第二章有限元的基本理论及三维有限元分析有限元方法是近似求解数理边值问题的一种数值技术,大约有40年的历史。他首先在本世纪40年代被提出,在50年用于飞机的设计。在六七十年代被引进到电磁场问题的求解中。2.1电磁场边值问题及其变分原理电磁场的边值问题和很多的物理系统中的数学模型中的边值问题一样,都可以用区域Ω内的控制微分方程(电磁场问题中可以是泊松方程、标量波动方程和矢量波动方程等)和包围区域
3、的边界Γ上的边界条件(可以是第一类的Dirichlet条件和第二类的Neumann条件,或者是阻抗和辐射边界条件等)来定义。微分方程可表示为:(2.1)式中,是微分算符,是激励函数,是未知量。对于电磁场边值问题,只有少数情况可以得到解析解。很多的时候我们采用基于变分原理的数值方法去求其近似解,比如伽辽金方法。在伽辽金方法中,我们首先定义非零的残数:(2.2)的最佳近似应能满足:(2.3)这里表示残数加权积分(也可称为误差泛函),是所选择的加权函数。进一步地,我们可以将近似解展开为:(2.4)式中,是定义在区域Ω内的展开函数,是待定的展开系数。并且我们将加权函
4、数选为:(2.5)这时,式(2.3)变为:(2.6)这样问题的求解就转化为能够使上式最小化的展开系数的线性问题的求解,将(2.6)式写为矩阵形式:(2.7)的元素为:(2.8)的元素为:(2.9)2.2有限元方法的原理――从一维的例子来看其建模的过程从上一小节的内容我们可以看到电磁场边值问题变分解法的这样的两个特点:(1)变分问题已经将原来电磁场边值问题的严格求解变为求解在泛函意思下的弱解,这个解可以和原来的解式不一样的。(2)在电磁场边值问题的变分方法中,展开函数(也可成为试探函数)是由定义在全域上的一组基函数组成,这种组合必须能够表示真实解,也必须满足
5、适当的边界条件,这对于二维、三维问题是非常困难的。很自然的,人们认为如果采用组成全域的子域上的一组基函数能够提高近似解对于真实解的逼近精度。这就是有限元方法。下面我们通过一个简单的一维例子来看看有限元方法的建模过程和其方法的特点。考虑一个均匀充填介电常数为ε的平板电容器,如图2.1所示:如果我们假设电场只有x方向的分量,问题就可以简化为一维问题。问题的支配方程为:(2.10)其边界条件为:(2.11)利用(2.10)式与权函数构成内积,仿照(2.3)式的方法我们可以给出这里的误差泛函:(2.12)如图2.2所示,我们可以将一维区域离散化为N段(单元),每一小
6、段又有编号为“1”和“2”的两个端点(结点),也称为“本地”序号,当然,与单元一样每个结点还有相应的全域序号。如果我们假设在单元内部电位函数按照线性规律变化,也就是对于单元内部的函数进行一阶插值:(2.13)特别的,在两个结点和处我们令其值分别为和,则(2.13)式可以重新写为(实际上和就成为了这一子域上的待求的系数):(2.14)其中:,,,,那么这时候在离散化的意义下,泛函(2.12)式可以写为:(2.15)其中,k是结点的全域序号,K是所有结点的总数,是第k个结点的子域。由于结点和单元的关系,我们可以在单元内选取(i=1,2)做为权函数,在利用一些矢量
7、运算恒等式,我们可以得到:(2.16)式中,n为单元的序号,N为总的单元数。注意到在离散化子域上有:(2.17)(2.18)实际问题中,应该是域内无源,所以为零。则在每个单元内(2.16)式的左边可以写为线性表达式:(2.19)(2.19)具体的我们可以用图2.1所示的例子来进行数值实现。在图2.1的离散化情况下我们有3个未知数,即对应结点全域序号的,和(而其中的和又有边界条件给定)。首先将(2.19)式对应单元1中的线性表达式的值带入到求解全部3个未知数的全域矩阵中:(2.20)再将(2.19)式对应单元2中的线性表达式的值带入到求解全部3个未知数的全域矩
8、阵中,构成全域矩阵方程:(2.21)再在(2.21)
此文档下载收益归作者所有