预览加载中,请您耐心等待几秒...
1/10
2/10
3/10
4/10
5/10
6/10
7/10
8/10
9/10
10/10

亲,该文档总共57页,到这已经超出免费预览范围,如果喜欢就直接下载吧~

如果您无法下载资料,请参考说明:

1、部分资料下载需要金币,请确保您的账户上有足够的金币

2、已购买过的文档,再次下载不重复扣费

3、资料包下载后请先用软件解压,在使用对应软件打开

有限元方法有限元法是求解偏微分方程问题的一种重要数值方法,它的基础分两个方面:一是变分原理,二是剖分插值.从第一方面看,有限元法是Ritz-Galerkin方法的一种变形.它提供了一种选取“局部基函数”的新技巧,从而克服了Ritz-Galerkin方法选取基函数的固有困难.从第二方面看,它是差分方法的一种变形.差分法是点近似,它只考虑在有限个离散点上函数值,而不考虑在点的邻域函数值如何变化;有限元方法考虑的是分段(块)的近似.因此有限元方法是这两类方法相结合,取长补短而进一步发展了的结果.在几何和物理条件比较复杂的问题中,有限元方法比差分方法有更广泛的适应性.§7.两点边值问题的有限元方法1.写出Ritz形式的变分问题2.区域剖分设为的有限维子空间,它的元素为.要构造,只需构造单元基函数.构造单元基函数所遵循的原则是:3.确定单元基函数2021/5/182021/5/182021/5/18(7.6)(7.8)(7.10)(7.11)(7.16)于是有(7.17)从总刚度矩阵和总荷载向量的形成过程可以看出,的计算,实际上是把中四个元素在适当的位置上“对号入座”地叠加,的计算也是如此.我们引入,只是为了叙述方便,实际上,在编制程序时并不需要.显然,方程组(7.18)的系数矩阵是对称正定的三对角矩阵,因此可采用追赶法求出在节点上的近似值.§7.两点边值问题的有限元方法与边值问题(7.1)、(7.2)等价的Galerkin变分问题是:求,使得(7.19)其中仍用分段线性函数构成的试探函数空间替代,将代入(7.19),则得到所满足的线性代数方程组(7.20)这和方程组(7.6)是完全一样的.与容易看出,方程组(7.20)的系数矩阵就是总刚度矩阵.在总刚度矩阵形成的过程中,注意到(7.21)而从而有即故有这就是有限元方程(7.18).由上述看出,按Galerkin法推导有限元方程更加直接方便.尤其重要的是.按这一观点推导的有限元方程,不仅适用于定常的微分方程定解问题,而且也适用于不定常的微分方程定解问题,因此具有广泛的适应性.例7.1用有限元方法解边值问题将区间[0,1]等分成4个单元.解利用上述分析结果,我们只需构造出单元刚度矩阵和单元荷载向量,然后合成为总刚度矩阵和总荷载向量.注意到(7.13)和(7.15),并将形成单元上的中点值则不难得到其中,,单元的中点为于是有如果把单元刚度矩阵和单元荷载向量“扩大”,便得到和为类似地,可写出和.然后进行叠加,便得到总刚度矩阵和总荷载向量:依边界条件即在中划去首末两行和首末两列,在中划去首末两行,便得到如下线性代数方程组:解之,得§8.二维椭圆边值问题的有限元方法与边值问题(8.1)、(8.2)等价的Galerkin变分问题是:求,使得(8.3)其中8.1区域剖分正如前章所言,对高维区域的剖分与对一维区域的剖分有很大不同.对一维区域无论作哪一种剖分,其单元仍然是一个区间,对不同的剖分只是区间长度不同而已.对高维区域而言,不同的剖分其单元的形状各异,如对二维区域,剖分后的子区域可以是三角形、矩形或四边形.限于篇幅,本书只讨论剖分后所得的子区域是三角形的情况,这种剖分称为三角形剖分.将区域划分成有限个三角形单元,剖分方法见前章,那里曾假定剖分的单元应是锐角三角形.现在我们去掉这一限制,只假定不同的单元是无重叠的内部,且单元的顶点不是其它单元边的内点.当然还要尽量避免出现大钝角的三角形.在物理量变化剧烈的地方,单元要划分得细密一些,变化缓和的地方,划分得稀一些.划分好单元之后,要对单元和节点进行编号.设是区域中的单元总数,将全区域中的单元统一编号,单元号记为.全区域中的节点也要按一定的顺序统一编号,记全区域中共有个节点,节点号记为节点编号的一般原则是尽可能使同一单元内的节点号比较接近.以后可以看到,单元内节点序号的差值决定了总体系数矩阵的带宽.8.2确定单元基函数与一维情形一样,为了构造试探函数空间我们只需在每个单元上构造插值基函数.这里,我们仅考虑三角形单元上的线性插值函数.为了便于后面积分的计算,我们先将直角坐标转换为面积坐标.1.面积坐标及有关公式(1)面积坐标的定义设是以为顶点的任意三角形单元,面积为,我们规定的次序按逆时针方向排列.在中(图8.1)任意一点的位置,可用它在直角坐标系中的两个坐标值来确定.如果我们过点作与三个顶点的连线,形成三个小三角形那么一旦的值确定后,这三个三角形的面积也就有了确定的值;反之,这三个小三角形的面积确定之后,点也就有了确定的位置,由此可见,三角形单元中任意一点的位置,除了可用直角坐标来确定外,还可以用连接点与与三个节点所形成的小三角形的面积来确定.用分别表示这三个小三角形的面积,显然令(8.4)则称这三个比值为点的面积坐标.由定义可知,所以,并不是互相独立的,