61阅读

有限元法基础-边界元法:边界元法-简介,边界元法-边界元法的基础

发布时间:2017-10-06 所属栏目:数值计算方法

一 : 边界元法:边界元法-简介,边界元法-边界元法的基础

边界元法所属现代词,指的是一种继有限元法之后发展起来的一种新数值方法,与有限元法在连续体域内划分单元的基本思想不同,边界元法是公在定义域的边界上划分单元。

边界元法_边界元法 -简单介绍

边界元法是在有限元法之后发展起来的1种较精确有效的工程数值分析方法 。 又称边界积分方程-边界元法。它以定义在边界上的边界积分方程为控制方程,通过对边界分元插值离散,化为代数方程组求解。它与基于偏微分方程的区域解法相比,由于降低了问题的维数,而显著降低了自由度数,边界的离散也比区域的离散方便得多,可用较简单的单元准确地模拟边界形状,最终得到阶数较低的线性代数方程组。又由于它利用微分算子的解析的基本解作为边界积分方程的核函数 ,而具有解析与数值相结合的特点,通常具有较高的精度。特别是对于边界变量变化梯度较大的问题 ,如应力集中问题 ,或边界变量出现奇异性的裂纹问题,边界元法被公认为比有限元法更加精确高效。由于边界元法所利用的微分算子基本解能自动满足无限远处的条件,因而边界元法特别便于处理无限域以及半无限域问题。边界元法的主要缺点是它的应用范围以存在相应微分算子的基本解为前提,对于非均匀介质等问题难以应用,故其适用范围远不如有限元法广泛,而且通常由它建立的求解代数方程组的系数阵是非对称满阵,对解题规模产生较大限制。对一般的非线性问题,由于在方程中会出现域内积分项,从而部分抵消了边界元法只要离散边界的优点。

边界元法_边界元法 -边界元法的基础

边界元法是基于控制微分方程的基本解来建立相应的边界积分方程,再结合边界的剖分而得到的离散算式。
Jaswon和Symm于1963年用间接边界元法求解了位势问题;Rizzo[3]于1967年用直接边界元法求解了二维线弹性问题;Cruse[4]于1969年将此法推广到三维弹性力学问题。1978年,Brebbia用加权余量法推导出了边界积分方程,他指出加权余量法是最普遍的数值方法,如果以Kelvin解作为加权函数,从加权余量法中导出的将是边界积分方程——边界元法,从而初步形成了边界元法的理论体系,标志着边界元法进入系统性研究时期。

边界元法_边界元法 -边界元法的发展

经过近40年的研究和发展,边界元法已经成为1种精确高效的工程数值分析方法。在数学方面,不仅在一定程度上克服了由于积分奇异性造成的困难,同时又对收敛性、误差分析以及各种不同的边界元法形式进行了统一的数学分析,为边界元法的可行性和可靠性提供了理论基础。在方法与应用方面,现在,边界元法已应用到工程和科学的很多领域,对线性问题,边界元法的应用已经规范化;对非线性问题,其方法亦趋于成熟。在软件应用方面,边界元法应用软件已由原来的解决单一问题的计算程序向具有前后处理功能、可以解决多种问题的边界元法程序包发展。
我国约在1978年开始进行边界元法的研究,目前,我国的学者在求解各种问题的边界元法的研究方面做了很多的工作,并且发展了相应的计算软件,有些已经应用于工程实际问题,并收到了良好的效果。

边界元法_边界元法 -图书信息

书 名: 边界边界元法:边界元法-简介,边界元法-边界元法的基础_边界元法元法
作 者:姚振汉
出版社:高等教育出版社
出版时间: 2010年4月1日
ISBN: 9787040286090
开本: 16开
定价: 59.00元

边界元法_边界元法 -内容简单介绍

《边界元法》内容简单介绍:边界元法是在有限元法之后发展起来的1种精确高效的工程分析数值方法。经过近五十年的发展,它不仅在固体与结构分析领域成为有限元法最重要的1种补充,而且在微机电系统电磁场分析和大型结构电磁波散射分析等领域也得到广泛应用。
《边界元法》分为传统边界元法的基本内容和近年发展的快速多极边界元法等新进展2大部分。前七章包含了传统边界元法的基本内容,分为3个单元:前三章为数学力学基础部分,介绍各种问题边界积分方程的建立;第四、第五章为基本数值方法部分,包括分元离散,数值积分和方程求解,并结合二维问题介绍其程序实现;第六、第七章为几类应用专题,主要是含时间问题、几种非线性问题和反问题。
第八、第九章介绍快速多极边界元法和大规模快速多极边界元并行算法,第十二章介绍与边界积分方程相关的边界型无网格法。另外在第十、第十一两章简要介绍国际上边界元法比较成功的应用,包括在机械、结构工程中的应用,和声场、电磁场分析设计中的应用。
书中的内容多于48学时或32学时的课程能够讲授的内容,便于不同学校、不同专业的老师根据需要选讲部分内容,同时为研究生提供课外的补充学习材料。《边界元法》附带光盘,提供弹性力学平面问题的边界元法C++和Fonran源程序、1个三维位势问题的常规和快速边界元分析程序的执行文件,以及相应的考题和算例,供读者试用。
《边界元法》也可以作为有关教师和工程技术人员学习边界元法的参考书。

边界元法_边界元法 -图书目录

引言
1 边界元法的数学基础
2 边界元法的发展历史
3 我国边界元法研究概况
4 边界元法研究的最新进展
5 边界元法的应用举例
6 边界元法的优缺点
7 本书的内容安排
参考文献
第一章 位势问题的边界积分方程与边界元法
1 调和方程的基本定解问题
2 Green等式、基本解及解的积分表达式
3 边界积分方程的建立
4 对于一般问题的推广
5 位势问题的边界元法简单介绍
习题
附录1指标符号与笛卡儿张量简单介绍
A1.1 指标符号
A1.2 矢量
A1.3 张量和张量场
参考文献
第二章 线弹性静力学问题的边界积分方程
1 线弹性静力学定解问题的微分提法
2 Betti定理、Kelvin解及Somigliana等式
3 线弹性静力学的边界积分方程
4 建立基本解的1种一般方法
习题
参考文献
第三章 几种常见的直接法和间接法边界积分方程
1 核函数的扩充
2 回转体问题
2.1 变截面轴的扭转问题
2.2 轴对称问题
2.3 回转体的弯曲问题
3 弹性薄板弯曲问题
3.1 弹性薄板弯曲问题的微分提法
3.2 弹性薄板弯曲问题的基本边界积分方程
3.3 弹性薄板弯曲问题的补充边界积分方程
4 弹性裂纹问题的对偶边界积分方程
4.1 位移边界积分方程
4.2 面力边界积分方程
5 半空间、半平面问题
5.1 半空间问题
5.2 半平面问题
6 位势问题的间接法边界积分方程
7 虚应力法建立的边界积分方程
8 位移间断法建立的边界积分方程
9 域外回线虚载荷法建立的回线积分方程
10 域外奇点法建立的边界积分方程
11 边界积分方程的正则化和基本解的恒等式
习题
参考文献
第四章 二维问题的边界元数值方法与程序实现
1 边界的离散化
1.1 二维域边界线的几何描述及单元自动划分
1.2 二维域的边界线元单元描述
2 边界积分方程的离散化
2.1 由加权余量法配点格式将边界积分方程化为线性代数方程组
2.2 核函数与形函数乘积的等精度Gauss积分
2.3 奇异积分的处理
2.3.1 弱奇异积分的处理
2.3.2 Cauchy主值积分和超奇异积分的简单特解法
2.3.3 Cauchy主值积分和超奇异积分的有限部分积分的一般性处理方法
3 方程的求解以及边界应力、内点位移和应力的确定
3.1 离散化的边界积分方程的求解
3.2 边界应力的确定
3.3 内点位移和应力的确定
4 边界元法计算误差的1种直接估计
4.1 内点变量趋于边界点极限的确定
4.2 边界元解误差的1种直接估计
4.3 基于边界元解误差直接估计的边界元自适应计算简例
5 边界元子域法
5.1 边界元链状子域法
5.2 边界元重复相似子域法
5.3 边界元行列子域法
习题
附录4 弹性力学平面问题边界元分析软件部分源程序
A4.1 FORTRAN程序BIEBE2说明
A4.2 C++源程序
参考文献
第五章 三维问题的边界元数值方法
1 边界的离散化
1.1 用I-J映射法自动划分单元
1.2 三维域的边界面元单元描述
2 边界积分方程的离散化
2.1 核函数与形函数乘积的等精度Gauss积分
2.2 弱奇异积分的处理
2.3 奇异积分和近奇异积分的简单特解法
2.4 Cauchy主值积分的直接计算法
2.5 超奇异积分的有限部分积分
3 线性代数方程组的求解
4 裂纹问题对偶边界元法
4.1 裂纹面的分元离散
4.2 确定应力强度因子的方法
5 边界元一有限元耦合方法
习题
附录5 三维位势问题边界元分析软件使用说明
A5.1 FMBEM-LAPLACE3D简单介绍
A5.2 使用方法
A5.2.1 前处理
A5.2.2 运行程序
A5.2.3 结果信息
A5.3 前处理转换软件PAq、一TRANSI。ATOR使用说明
A5.3.1 生成MSC.PATRAN前处理文件
A5.3.2 生成BEM_INP.DAq、前处理文件
参考文献
第六章 与时间有关问题的边界元法
1 瞬态热传导问题
1.1 Laplace变换法
1.2 边界元一时间差分耦合法
1.3 与时间有关的基本解
2 弹性动力学问题
2.1 基于与时间有关基本解的边界积分方程与边界元法
2.1.1 与时间有关的基本解
2.1.2 时间一空间域的边界积分方程
2.1.3 时间一空间域的弹性动力学边界元法
2.2 Laplace变换法
2.2.1 在Laplace变换空间的边界积分方程
2.2.2 边界积分方程的离散
2.2.3 Laplace反演方法
2.3 双重互易法
习题
附录6 弹性动力学边界元法补充公式
A6.1 弹性动力学二维问题的时间一空间域基本解
A6.2 弹性动力学三维问题时间一空间域基本解的时间积分
A6.3 弹性动力学三维问题的1种新的时空域边界积分方程
A6.4 弹性动力学三维问题的1种新的高效的时空域边界元法
参考文献

边界元法_边界元法 -边界元法在传热学中的应用

前言

数值计算方法及计算机仿真是目前解决现代工程技术问题的
1个重要手段,在航空、航天、化工、机械、冶金、电力、能源、安全、
仪器仪表、计量测试等各个领域都有着广泛的应用。特别是近几
年来计算机硬件和软件技术的高速发展,为数值计算与仿真研究
开辟了广阔的应用前景。
常见的数值计算的离散方法主要有有限差分法、有限元法和
边界元法等。有限差分法是历史上最早采用的数值方法,这种方
法是用网格线的交点为节点,在每个节点上控制方程中的微分用
差分代替,从而在每个节点形成1个代数方程,求解这些代数方程
即可获得所需的数值解。有限差分法的优点是最容易实施,缺
点是数值解稳定性难以保证,并且对复杂区域的适应性较差。有
限元法则是把区域划分成一系列单元,在每个单元上取数个点作
为节点,然后通过对控制方程积分来获得离散方程。有限元法的
优点是对不规则区域的适应性好,但计算量较大,需要划分巨大数
量的单元。这2种方法的共同缺点是存在解的不稳定性问题。有
限差分法和有限元法的研究和应用已经比较成熟,并且都已有了
商业化的通用软件。边界元法是应用格林公式选择适当的权函数
把空间求解域上的偏微分方程转换成为其边界上的积分方程。边
界元法的最大优点是使求解空间的维数降低了一阶,从而大大减
少了计算时间和存储容量。另外,由于边界元法不需要对区域内
进行离散,因此温度的测量点位置即可任意选择。更进1步的
是,不像其他数值方法,边界元法可以直接给出未知表面的温度和
热流,而表面的热流的精确测量要比温度的测量困难得多。这些
优点极大地提高了边界元法在实际工程中的应用。
由于边界元法的研究相对较晚,目前在国际上还没有很好的
通用软件。在国内有关边界元法解决实际工程问题的论文和论著
都比较少,特别是边界元法数值计算的实用程序也难得一见。本
书是根据作者多年来在边界元法和数值传热学研究的基础上完成
的。在书中作者提供了许多边界元法实际应用的c语言源程序,
包含了许多实用的编程技巧。这些程序都经过Vc++严格调试,
每个子程序都是独立模块,并且自成1个完整的体系,不需要任何
其他环境,可以独立运行计算。希望能为使用边界元法解决实际
工程问题的科技工作者提供有用的参考。
本书的第1章介绍了边界元法的理论基础,其中的公式都进
行了严格的推导。第2章介绍了边界元法在稳态导热问题中的应
用。第3章介绍了边界元法在非稳态导热问题中的应用。第4章
介绍了边界元法在导热反问题中的应用。第5章给出一些数值计
算中要用到的相关的数值计算方法。
作者十分感谢中国计量学院的崔志尚教授、李希靖教授、袁昌
明教授在科研和学术上的指导和帮助,感谢国防工业出版社江洪
湖编辑的大力支持。由于作者水平有限,书中难免有缺点和错误,
敬请读者批评指正。
吴洪潭
于中国计量学院

目 录

第l章边界元法的理论基础……………………………
1.1边界元法概述……………………………………
1.2加权余量法………………………………………
1.3 函数……………………………………………
1.4基本解…………………………………………·
1.5边界积分方程…………………………………·
l_6格林公式………………………………………·
1.7常单元…………………………………………·
1.8线性单元………………………………………·
1.9 次单元………………………………………·
lI 10角点处理………………………………………·
参考文献………………………………………………·
第2章边界元法在稳态导热问题中的应用…………·
2.1 --维平面稳态导热问题的边界元分析………-
2.1.1 混合边界条件下的边界积分方程…·
2.1.2混合边界条件下的边界离散方程…·
2.1.3 二维平面稳态导热问题的通用边界元
C程序…………………………………一
2.2轴对称稳态导热问题的边界元分析…………·
2.2.1边界积分方程………………………一
2.2.2基本解……………………··
2.2.3边界元方程及单元插值…一
参考文献……………………………………一
第3章边界元法在非稳态导热问题中的应用
3.1边界积分方程………………………一
3.2边界元方程…………………………一
3.3时间步长划分………………………-.
3.4区域积分项…………………………一
3.5边界单元的离散……………………··
3.5.1 非对角元素的计算………··
3.5.2 对角元素的计算…………··
3.6 非稳态导热问题通用边界元C程序··
3.6.1 域内划分十六个四边形……··
3.6.2 域内划分八个三角形……一
3.6.3 域内划分160个三角形…··
参考文献………………………………………一
第4章边界元法在导热反问题中的应用……
4.1 导热反问题概述……………………一
4.2 边界形状导热反问题………………
4.3共轭梯度法…………………………
4.4边界元数值计算……………………
4.5 数值模拟实验………………………
4.6 导热反问题的边界元C程序………
参考文献………………………………………
第5章相关的数值方法………………………·
5.1 高斯求积公式………………………-
5.2 辛普生求积公式……·
5.3 线性方程组的数值解法
5.4 随机数的产生………·
参考文献………………………·

二 : 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?中的任何二个发生共线的情况。

本文标题:有限元法基础-边界元法:边界元法-简介,边界元法-边界元法的基础
本文地址: http://www.61k.com/1081736.html

61阅读| 精彩专题| 最新文章| 热门文章| 苏ICP备13036349号-1