一 : 边界元法:边界元法-简介,边界元法-边界元法的基础
边界元法所属现代词,指的是一种继有限元法之后发展起来的一种新数值方法,与有限元法在连续体域内划分单元的基本思想不同,边界元法是公在定义域的边界上划分单元。
边界元法_边界元法 -简单介绍
边界元法是在有限元法之后发展起来的1种较精确有效的工程数值分析方法 。 又称边界积分方程-边界元法。它以定义在边界上的边界积分方程为控制方程,通过对边界分元插值离散,化为代数方程组求解。它与基于偏微分方程的区域解法相比,由于降低了问题的维数,而显著降低了自由度数,边界的离散也比区域的离散方便得多,可用较简单的单元准确地模拟边界形状,最终得到阶数较低的线性代数方程组。又由于它利用微分算子的解析的基本解作为边界积分方程的核函数 ,而具有解析与数值相结合的特点,通常具有较高的精度。特别是对于边界变量变化梯度较大的问题 ,如应力集中问题 ,或边界变量出现奇异性的裂纹问题,边界元法被公认为比有限元法更加精确高效。由于边界元法所利用的微分算子基本解能自动满足无限远处的条件,因而边界元法特别便于处理无限域以及半无限域问题。边界元法的主要缺点是它的应用范围以存在相应微分算子的基本解为前提,对于非均匀介质等问题难以应用,故其适用范围远不如有限元法广泛,而且通常由它建立的求解代数方程组的系数阵是非对称满阵,对解题规模产生较大限制。对一般的非线性问题,由于在方程中会出现域内积分项,从而部分抵消了边界元法只要离散边界的优点。边界元法_边界元法 -边界元法的基础
边界元法是基于控制微分方程的基本解来建立相应的边界积分方程,再结合边界的剖分而得到的离散算式。边界元法_边界元法 -边界元法的发展
经过近40年的研究和发展,边界元法已经成为1种精确高效的工程数值分析方法。在数学方面,不仅在一定程度上克服了由于积分奇异性造成的困难,同时又对收敛性、误差分析以及各种不同的边界元法形式进行了统一的数学分析,为边界元法的可行性和可靠性提供了理论基础。在方法与应用方面,现在,边界元法已应用到工程和科学的很多领域,对线性问题,边界元法的应用已经规范化;对非线性问题,其方法亦趋于成熟。在软件应用方面,边界元法应用软件已由原来的解决单一问题的计算程序向具有前后处理功能、可以解决多种问题的边界元法程序包发展。边界元法_边界元法 -图书信息
书 名: 边界元法边界元法_边界元法 -内容简单介绍
《边界元法》内容简单介绍:边界元法是在有限元法之后发展起来的1种精确高效的工程分析数值方法。经过近五十年的发展,它不仅在固体与结构分析领域成为有限元法最重要的1种补充,而且在微机电系统电磁场分析和大型结构电磁波散射分析等领域也得到广泛应用。边界元法_边界元法 -图书目录
引言边界元法_边界元法 -边界元法在传热学中的应用
前言
数值计算方法及计算机仿真是目前解决现代工程技术问题的目 录
第l章边界元法的理论基础……………………………二 : 93有限元法的基础理论
61阅读/ www.61k.net
一、里兹法与迦辽金法(摘自电磁场有限元方法 金建铭)
1. 里兹法
里兹法是一种变分方法,其中边值问题用变分表达式(也称泛函)表示,泛函的极小值
对应于给定边界条件下的控制微分方程。通过求泛函相对于其变量的极小值可得到近似解。
2. 伽辽金法
伽辽金法属于残数加权方法类型,它通过对微分方程的残数求加权的方法得到方程的
解。
?是方程的近似解,将u?代入方程可得到非零的残数: 若u
??f r?Lu
?的最佳近似应能使残数r在?内所有点上有最小值。残数加权方法要求: u
Ri???ird??0 ?
这里Ri表示残数的加权积分,?i是所选的加权函数。
在伽辽金法中,加权函数与近似解展开中所用的函数相同。通常,这样可得到最精确的
解。
二、有限元方法
里兹法和伽辽金法中,在整个解域内找出能表示或至少近似表示问题真实解的试探函数
是非常重要的。然而对于许多问题,这个步骤是十分困难的,对二维和三维问题尤其如此。
为此,我们可将整个区域划分成小子域,并应用定义在每个子域上的试探函数。因为子域是
小区域,因而在每一子域内函数的变化不大,所以定义在子域上的试探函数通常比较简单。
这正是有限元法的基本思想。应用里兹法的过程通常称为里兹有限元法或变分有限元法,而
应用伽辽金方法的过程通常称为伽辽金有限元方法。
有限元法与经典里兹法和伽辽金法的不同之处是在试探函数的公式上。在经典里兹法和
伽辽金法中,试探函数由定义在全域上的一组基函数组成。这种组合必须能够(至少近似)
表示真实解,也必须满足适当的边界条件。在有限元法中,试探函数是由定义在组成全域的
子域上的一组基函数构成。因为子域很小,所以定义在子域上的基函数能够十分简单。
三、关于形函数(摘自有限元法在电磁计算中的应用 张榴晨)
对于一个待求的微分方程,用一组线性独立的尝试函数?i和待定系数Ci来表示方程的
近似解,并用加权余数法(迦辽金法)来求解这些待定系数。求解待定系数的代数方程组为:
?[?i?1n???i??jd?]Ci??q?jd??j?1,2,?,n
这里?j为所选择的加权函数,应用迦辽金法时,所选取的加权函数即为尝试函数。
有限元中应用的尝试函数代表了单元上近似解的一种插值关系,它决定了近似解在单元
上的形状。因此尝试函数在有限元法中又称为形函数。对于一维有限元来说,形函数为一个
直线段;对一维高阶有限元来说,形函数为一个曲线段;对二维一阶有限元来说,形函数为
一个平面;对二维高阶有限元来说,形函数为一个曲面;三维有限元来说,形函数为多维平
面或曲面。选择形函数时可以使一个任意元上的函数只与该元所对应的节点势函数值有关,
而与其它各点的值无关。
1. 一维有限元
对于一维有限元来说,形函数分段线性。该形函数?i在节点i上的值为1,并在与节点
i相邻的两个单元上线性减小,直到在相邻节点i?1和i?1上分别减小为零。其形函数的形
式为:
?i??ix??i
根据形函数的性质即可得到两个线性方程组,解这两个线性方程组即得到系数?i和?i,
进而得到?i。应用同样的方法可以得到每一个节点上的形函数表达式。
2. 二维有限元
以三角形单元为例,单元的顶点分别为i,j,k。每个顶点都对应于一个形函数。考虑一
阶有限元法,即用线性插值的方法表示一个有限元上的势函数分布。这样,三角形作为一个
单元所对应的三个形函数都是分段线性的函数,并可以看作由几个平面组成。对应于节点i
的形函数?i在节点i上取单位值(即为1),并由此在单元e上直线下降,直到在其它另外
两个节点j和k上降为零,在单元之外的区域则一直保持零值。于是这三个节点值就决定了
该形函数的形状。形函数?i可由一个线性函数表达为:
?i??i??ix??iy
四、势函数分布
有限元的作用在于求解分布场的势函数在每个节点上的近似值,而势函数在单元的其它
位置的值,可用插值的原理来表示。如果采用线性插值的方法来表示分布势函数,则称为一
阶有限元法,如果有限元法采用高阶插值法表示分布势函数,则称为高阶有限元法。
任意一个单元上的势函数分布由这个单元上的节点势函数值及相应的形函数表示,对
于一个一维单元有:
?e????i?1 ii??i?1
整个区域的势函数分布则由每一个单元的势函数分布相加得到。
五、关于对称性的利用
利用对称性可以减少节点和单元的数目,从而节省用于建立模型和计算近似解的时间。
从另外一个角度来讲,如果维持节点和单元的数目不变,则利用对称性可以对这四分之一区
域作更为详细地划分,即单元的尺寸可大为减小,从而提高近似解的精度。
case 1(a):一个两极电容器的静电场分布问题。假设两个极板都接在0.5v电源端,极板间
的距离为2,极板间充有密度为?的自由电荷。这样,电容器的激励和几何形状都关于y轴
对称。由于这种关系的存在,我们只要求解电场在整个区域的一半的分布即可,而另一半的
电场分布则能从对称关系而得到。从边界条件的观点出发,这种对称的结构导致电力线垂直
穿过y轴,使电势在该对称轴上沿x方向的变化率为零。
描述这一问题的微分方程和边界条件为:
?2????x?(0,1)
?x?1?0.5
?0??
?nx?0
见原文page62
case 1(b):当考虑极板端部的边缘效应时(考虑一个方向的边缘效应),这个问题则为一个
二维边值问题。假设电势分布在与极板垂直平等的截面上分布相同,那么只需计算一个截面
的电势分布即可。在y轴两侧,几何结构及施加电压都对称,因此y轴为对称轴,即沿y
轴的电势对法线方向(x方向)的变化率为零(是否可以理解为等势面与y轴垂直而与x轴
平行),这样y轴便构成一个齐次诺伊曼边界条件。从另外一个角度看,该问题的几何结构
x轴对称,而施加电压则x轴反对称,因此沿x轴的电势应取两极板施加电压的中值,即零
值。因此x轴构成一个狄利克莱边界条件。于是计算区域被减为原区域的四分之一。
为了把边缘效应考虑进去,计算区域还应考虑电容器周围空间,并应将无限远处设为
零电位参考点。但从实际意义来讲,假设电容器外某一定距离的空间处为零电位参考点即可
满足实际需要。设沿x轴方向的2m处电位为零。
描述这一问题的微分方程和边界条件为:
?2??0
?
??1
?2?10??
?4?3?0 ???n?0
case 2:一个同轴传输线,两个同芯长方形导体间充满了线性介质,其介电常数为?。假设
两导体间加有直流电压10V,导体间贮有密度为?的自由电荷,传输线的长度远远大于其
截面的长和宽,那么可以认为电场在传输线各个截面上的分布都相同,因此只需求解电场在
某个截面的分布,从而检查绝缘材料的工况。该问题的几何形状、介质及激励都x轴和y
轴对称,因此只需求解整个截面的四分之一即可。在对称轴的两侧,电势对于该轴线的法向
变化率为零。从电力线的观点出发,也可以说电力线垂直对称轴线。
描述这一问题的微分方程和边界条件为:
?2???qx,y??(q??)?
???1
?3?0?10
?2???n???
?n?4?0
见原文page68。
case 3:一个简单的变压器的静电场分布。变压器结构及激励关于x轴对称,因此沿x轴有齐次的诺依曼条件成立;变压器的几何结构y轴对称,且激励y轴反对称,因此沿y轴存在着磁势为零的狄利克莱条件。这样,只需选取变压器的四分之一就能完全解出磁场分布。
这一静磁场的描述方程为帕松方程及拉普拉斯方程,由于只考虑二维结构,磁势只在沿z轴的方向不为零,而沿其他方向的分量均为零。对于这样一个二维问题的矢量磁势A实际上被简化为一个标量(或为沿一个方向的矢量)Az:
1
?r
??0
??
?n?2Az???0J ?0
四、高阶有限元
一阶有限元的形函数是一个线性函数,应用一阶有限元形成的系数矩阵为稀疏矩阵,其求解简单但精度较低。解决精度问题的方法之一为采用高阶有限元法。
对于一维问题用曲线来逼近,即采用高阶插值的方法通常可以得到更为精确的结果。例如采用二阶插值:
?i?ai?bix?cix2(i?1,2,3)
二阶有限元由三个点构成,其中的两个点为构成该有限元区域的端点,而第三个点通常选在这两点中间,这三个点分别对应三个二阶形函数,也就是说在一个单元中,每个节点都对应一个二阶形函数。
一般来说高阶有限元允许我们我用尺寸较大的单元来描述具体工程问题,特别是几何形状较复杂的问题。从另外一个角度来说,如果保持单元的尺寸不变,相对于一阶有限元来说,高阶有限元能够提高近似解的精确度。但高阶有限元也增加了计算量,因此更耗时。
Page106给出了一个结构分析中所用的几种三维有限元法的比较,其结果很能说明有限元形状和疏密对计算精度的影响。
五、单元与插值函数(有限单元法基本原理和数值方法 王勖成)
一个函数在域内其本身连续,它的一阶导数具有有限个不连续点但在域内可积,这样的函数称之为具有C0连续性的函数。类似地,如果微分算子A出现的最高阶导数是n阶,则要求函数u必须有连续的n?1阶导数,即函数应具有Cn?1阶连续性。
关于单元插值函数的形式,有限单元法中几乎全部采用不同阶次幂函数的多项式。这是因为它们具有便于运算和易于滑坡路收敛性要求的优点。如果采用幂函数多项式作为单元的插值函数,对于只满足C0连续性的单元(称C0型单元),单元内的未知场函数的线性变化能够仅用角(或端)结点的参数表示。对于它的二次变化则必须在角(或端)结点之间的边界上适当配置一个边内结点(二次单元)。它的三次变化则必须在每个边界上配置二个边内结点(三次单元)。配置边内结点的另一原因是常常要求单元的边界是曲线的,沿边界配置适当的边内结点可以构成二次或更高次多项式来描述它们。
1. 一维单元
对于n 个结点的一维单元,Ni可以采用n?1次Lagrange插值多项式li(n?1)(x),即令:
Ni(x)?li(n?1)(x)??nx?xj
xi?xjj?1,j?i?(x?x1)(x?x2)?(x?xi?1)(x?xi?1)?(x?xn)(xi?x1)(xi?x2)?(xi?xi?1)(xi?xi?1)?(xi?xn)
(i?1,2,?,n)
构造插值函数时,一般用无量纲表示,一维单元中的无量纲为长度坐标,更一般化的可称为自然坐标。
为构造形式的Lagrange单元方便,还可把插值多项式写为:
Ni(x)?li(n?1)(x)??nfj(?)
fj(?i)
j?1,j?i
其中fj(?)????j表示任一点?至点?j的距离,也是j点坐标???j表示成方程形式fj(?)????j?0的左端项。显然,fj(?j)?0。Ni?li(n?1)的展开式中包含了除fi(?)外的所有fj(?)(j?1,2,?,i?1,i?1,?,n)的因子,从而保证了Ni(?j)?0(j?i)这一要求。
这一因子引入li(n?1)的分母是为了保证满足fj(?i)??i??j是点I的坐标代入fj(?)的结果,
Ni(?i)?1的要求。
2. 二维单元
1) 三角形单元
I. 一次单元
对于3结点三角形单元,引入面积坐标
Li?Ai/A
则单元插值函数可表示为: (i?1,2,3)
Ni?Li(i?1,2,3)
II. 二次单元
二次单元有六个结点,各结点的面积坐标分别标注在图中。(page98)需要构造的插值函数表示成:
Ni(x)??
(i)2fj(i)(L1,L2,L3)fj(i)(L1i,L2i,L3i)(i) j?1其中的fj(L1,L2,L3),fj(L1i,L2i,L3i)赋与了和三角形单元相对应的几何意义。
fj(L1,L2,L3)是通过除结点i以外所有结点的二根直线方程fj(L1,L2,L3)?0的左端项。如当i?1时,fj(1)(i)(i)分别是通过结点4,6的直线方程f1(L1,L2,L3)?L1?1/2?0和通过
(1)(i)(1)结点2,5,3的直线方程f2(L1,L2,L3)?L1?0的左端项。fj(L1i,L2i,L3i)中的L1i,L2i,L3i是结
点i的面积坐标。所以得到:
93有限元法的基础理论_有限单元法基本原理和数值方法
N1?
L1?1/2L1
1/21
?(2L1?1)L1
类似的方法可以得到其它的形函数。
III. 三次单元
为保证二维域三次多项式的完备性,三次单元应有10个结点,按与二次单元相同的步骤按画线法构造它的插值函数。Page98
2) 矩形单元
如果所研究问题的总体域是矩形的,采用矩形单元将比三角形单元更有效。
I. Lagrange矩形单元
构造任意的Lagrange矩形单元插值函数的一个简便而系统的方法是利用二个坐标方向适当方次Lagrange多项式的乘积。在第J行和第I列的结点上的插值函数是:
Ni?NIJ?lI(r)(?)lJ(p)(?)
虽然构造Lagrange插值函数很容易,但是这一类型的单元存在一定缺点,主要是出现了随插值函数方次增高而增加的内结点,从而增加了单元的自由度数,而这些自由度的增加通常并不能提高单元的精度。
II. Serendipity单元
通常我们希望结点仅配置在单元的边界上,并且在实际应用中常希望同一单元的不同边界有不同数目的结点,这样可以实现不同阶次单元之间的过渡,从而可能在求解的不同区域采用不同精度的单元。假定开始只有四处角结点,对应这些结点的插值函数可利用双一次Lagrange多项式构造,即:
Ni?
1
(1??0)(1??0)4
(i?1,2,3?)
其中,?0??i?,?0???i
如果增加边内结点,则与它对应的插值函数可以按划线法构造或直接表示成?(或?)方向二次和?(或?)方向一次Lagrange多项式的乘积。例如增加结点5(位于结点1,2之间),则:
N5?
1
(1??2)(1??)2
需要指出的是,对于现在5个结点的情况,N5满足N5j??5j(j?1,2,?,5)要求,而
??
(是否理解为插值函数在该结点上为Ni(i?1,2,3,4)不再满足Ni5??i5(j?1,2)的要求了。
1而在其它结点上为零,在1,2结点间增加了结点5后若仍然采用线性插值函数则对于N1将不满足插值函数的特性,因为N1在结点5处不为零,为了保证在结点5处N1为零则必须使插值函数曲线在该点过零。如图page103图3.15所示,这时的插值函数是一条曲线。)为满足要求,N1,N2需要修正为:
?1?1N1?N1?N5,N2?N2?N5
22
类似地,可以讨论增加边内结点6,7,8的情况,
?11
N1?N1?N5?N8,
22?11
N3?N3?N6?N7,
221
N5?(1??2)(1??),
21
N7?(1??2)(1??),
2?11N2?N2?N5?N6
22?11N4?N4?N7?N8
22 1
N6?(1??2)(1??)
21
N8?(1??2)(1??)
2
若5,6,7,8结点中任一个不存在,则对应的插值函数为0。
对于高次的Serendipity单元,可用同样的方法构造它的插值函数。例如对于p次单元的边内结点,它的插值函数可表示成?(或?)方向p次和和?(或?)方向一次Lagrange多项式的乘积,而对于角结点则其插值函数可表示成一双线性函数和用适当分数分别乘以相邻两个边界上的各个边界结点插值函数之和,以保证它在边内结点上的值也为0。因为角结点插值函数中双线性函数之附加项可随相邻边界上的边内结点的增减而变化,所以边内结点数是灵活的,这对于构造变结点数的过渡单元是很适合的。 3. 三维单元
1) 四面体单元
四面体单元与二维情况的三角形单元类似,插值函数是在三维坐标内的各次完全多项式。在各个面上的结点配置与同次的二维三角形单元相同,函数是相应二维的完全多项式。
根据三维四面体单元的几何特点,引进的自然坐标是体积坐标。 2) Serendipity单元 3) Lagrange单元
以上内容详见原文page106。 六、等参单元
对于一个给定问题的求解域,预期用较少的单元即可获得需要精度的解答,但是用较少的形状规则的单元离散几何形状比较复杂的求解域常会遇到困难,因此需要寻找适当的方法将规则形状的单元转化为其边界为曲线或曲面的相应单元。在有限元法中最普遍采用的变换方法是等参变换,即单元几何形状的变换和单元的场函数采用相同数目的结点参数及相同的插值函数进行变换。采用等参变换的单元称为等参单元。 1. 等参变换
为将局部(自然)坐标中几何形状规则的单元转换成总体(笛卡尔)坐标中几何形状扭曲的单元以满足对一般形状求解域进行离散化的需要,要建立一个坐标变换。
最方便的方法是将这种坐标变换也表示成插值函数的形式:
x??Ni?xi
i?1
m
y??Ni?yi
i?1
m
z??Ni?zi
i?1
m
其中m是用以进行坐标变换的单元结点数,xi,yi,zi是这些结点在总体(笛卡尔)坐标内的坐标值,N?为形函数,实际上它也是用局部(自然)坐标表示的插值函数。通过上式
建立起两个坐标系间的变换,从而将自然坐标内的形状规则的单元(母单元)变换为笛卡尔坐标内的形状扭曲的单元(子单元)。
坐标变换关系式和函数的插值表示式在形式上是相同的。如果坐标变换和函数插值采用相同的结点并采用相同的插值函数则称这种变换为等参变换。 2. 导数间的变换
函数Ni对?的偏导数可表示成为:
?Ni?Ni?x?Ni?y?Ni?????x????y???
z
?z??
对于其它两个坐标(?,?),可写出类似的表达式,将它们集合成矩阵形式则有:
???Ni???x?y
?z????????????????Ni???Ni?
?x???x????N??i??y?z???????????x??????????????Ni???N?
i?
y??J?????y? ?????Ni???x?y?z????
?Ni?????????
??
?????Ni?
???z??????z??
式中J为Jacobi矩阵,可记作?(x,y,z)
?(?,?,?)
,根据上述的插值关系,J可以显式地表示为
自然坐标的函数:
?m?N'
??m
?N'
m
i??ixii?i??
??N'i???N'
?N'
'2??
yiz1i??Nm?i?ii?i????????
?????
??xy1J??(x,y,z)??m?N'
?(?,?,?)????i??
xi
i?i?m
?N'?m
iy?N''i?N'2??N'?1
m???
x2y2iz?????N1
i??i?i
??i?i??????????????
????
m?????N'ixii?i??
?m
?N'iy?N'i?m???xm
ym
i?i??
?m
?N'i?????N'
z1?N'2ii?i????????
??
???
??
d?,d?,d?在笛卡儿坐标系内形成的体积微元是:
dV?d??(d??d?)
而:
d??
?x?y???d?i???d?j?z
??d?kd???x??d?i??y??d?j??z??d?k
d???x??d?i??y??d?j??z??d?k
其中,i,j,k是笛卡儿坐标x,y,z方向的单位向量。将d?,d?,d?代入dV可得:
z1?z?2????zm??
?x???xdV?
???x???y???y???y???z???z
d?d?d??Jd?d?d? ???z??
3. 等参变换的条件
两个坐标间一对一变换的条件是Jacobi行列式J不为0,等参变换作为一种坐标变换也必须服从这个条件。
如果J?0则表明笛卡尔坐标中体积微元(或面积微元)为0,即在自然坐标中的体积微元d?d?d?(或面积微元d?d?)对应笛卡尔坐标中的一个点,这种变换显然不是一一对应的。
笛卡尔坐标中的面积微元可直接表示成:
dA?d??d??d?d?sin(d?,d?)
而: dA?Jd?d?
于是有:J?
d?d?sin(d?,d?)
d?d?
可见,只在以下三种情况之一成立,即
d??0,或d??0,或sin(d?,d?)?0就会出现J?0。因此,为了保证变换的一一
对应,应防止因任意的二个结点退化为一个结点而导致d?,d?,d?中的任意一个为0,还应防止因单元过分歪曲而导致的d?,d?,d?中的任何二个发生共线的情况。
本文标题:有限元法基础-边界元法:边界元法-简介,边界元法-边界元法的基础61阅读| 精彩专题| 最新文章| 热门文章| 苏ICP备13036349号-1