• / 41
  • 下载费用:10 金币  

有限元法0099.ppt

关 键 词:
有限元法0099.ppt
资源描述:
1. 有限元法,上一讲,利用加权余数法和变分法将偏微分方程转化为代数方程组求解,该矩阵方程包括系数矩阵、激励源矩阵和边界矩阵,而计算这些矩阵的元素时,常常用到分部积分法。如果为了计算精度而选取很多个尝试函数,那么计算这些为数众多的分部积分既十分复杂又很费时间,并且很难用计算机进行数值计算。 因此,我们需要寻找一个改进的方法来简化计算,并设法利用计算机进行处理,有限元法就是其中的一种。,通过尝试函数的选取,近似解满足1类边界条件,,基本思想:,1. 有限元法,在有限元方法中,场域被分割成许多很小的子区域,通常称为“单元”或“有限元”。对所有子区域进行独立的处理和运算,便对一个整体问题进行局部化处理。通过选取恰当的尝试函数,使每个单元的计算都变得非常简单,经过对每个单元重复而简单的计算,再将其结果总和起来,便可以得到用整体矩阵表达的整个区域的解,这一整体矩阵又常常是稀疏短阵,可以更进一步简化和加快求解过程。由于计算机非常适合于重复性的计算和处理过程,所以整体矩阵的形成过程很容易使用计算机来实现,,以满足帕松方程的有源静电场为例:,2. 一维有限元法,由对称性,将上述问题化为一维问题,用有限元方法求解,2. 一维有限元法,场域单元化 (剖分) 场域(0,1)分割成四个“单元”,e1、e2、e3、e4。单元大小与精度和计算量有关。处理复杂场域和激励源带来极大方便。,作为一种数值计算方法,有限元法并非用来寻求问题的解析解。实际上,很多工程问题目前都无法找到解析解。有限元的作用就在于求解分布场的势函数在每个节点上的近似值,而势函数在单元的其它位置的值可用插值来表示:如果采用线性捅值的方法来表示分布势函数,则称为一阶有限元法,如果采用高阶插值法表示分布势函数,则称为高阶有限元法,通常一阶插值即可。,2. 一维有限元法,本例,场域分割成4个单元,5个节点,求场域内电势分布,转化为求5个节点的电位即可。 场域内其它点(各单元内)的电位,由5个节点电位来插值表示。(一阶插值、高阶插值),对一 维场域来说,单元就是一个线段; 对二维场域,有限元单元形状可为二角形、矩形等,单元形状对有限元的简化有影响,通常为三角形,2. 一维有限元法,单元化后,整体场域积分转化为局部单元求和每个单元可独立处理,把两个要求解的量联系起来,有限元中令待定系数就是节点电位,当然尝试函数要重新确定,2. 一维有限元法,总体系数矩阵化为了局部系数矩阵的求和这个场域的求解变为求解剖分后节点上的电位,,,2. 一维有限元法,重新确定尝试函数的思想,尝试函数本质上是在场域内对真解的近似。以前的尝试函数是针对整个场域来选取的(随意性大,有技巧经验因素)现在整个场域被剖分为N个单元,每个单元内的真解由该单元尝试函数线性组合来近似,(即尝试函数也局部化,)当单元足够多(剖分足够细)时每单元内的尝试函数可以非常简单(线段)尝试函数:一维一阶有限元尝试函数为直线段;二维一阶为平面 (形函数)任意单元内的形函数只与该单元节点坐标有关,与其它单元无关,2. 一维有限元法,某单元内尝试函数确定,,,,,,,,,,,,,,,,,,,,,,,2. 一维有限元法,某单元内尝试函数的表达(局部坐标),,,,,,,,,,,,,,,,,,,,,任意单元内尝试函数的形式相同,单元内才不为0,构成的局部系数矩阵对整体矩阵有贡献,在考查的单元之外,尝试函数为0,即i, j不属于同一单元时,构成的局部系数矩阵为0,只有当i, j=i+1属于同一单元时,构成的局部系数矩阵才不为0,才对整体矩阵有贡献。这样每个单元的计算分开、独立了,且每个单元的计算有是相同的、重复的,便于计算机求解。,,,2. 一维有限元法,局部系数矩阵的计算,2. 一维有限元法,局部系数矩阵元素的计算,q是问题的激励,2. 一维有限元法,局部系数矩阵的计算,任意单元e的局部系数矩阵,只与该单元的长度(坐标)有关,与其它单元无关,该矩阵的i,i+1的行列上有值,其它位置为0。,其它所有单元的系数矩阵结构形式完全相同,只是不同单元的长度le不同而已。 这样极便于计算机作重复计算。,局部单元系数矩阵相加就可构成整个问题的整体系数矩阵(称为封装),2. 一维有限元法,本例中局部系数矩阵,2. 一维有限元法,整体矩阵的封装,各单元系数矩阵按照节点编号次序,在整体矩阵的相应行列填充并相加即可,2. 一维有限元法,整个问题的代数方程组:,称为有限元刚度矩阵,但不能直接求解,需要消去1行、1列。,2. 一维有限元法,由边界条件对整个问题的代数方程组消元:,由问题的边界条件,第5 个节点电位为0.5V,已知,故消去该节点的方程:5行5列。必有这一步,实际上原K矩阵行列式的值为0,本质上是找参考电位,由于参考电位为0.5V,所以1-4号节点计算值+0.5,有限元计算与真
展开阅读全文
  微传网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
0条评论

还可以输入200字符

暂无评论,赶快抢占沙发吧。

关于本文
本文标题:有限元法0099.ppt
链接地址:https://www.weizhuannet.com/p-9811947.html
微传网是一个办公文档、学习资料下载的在线文档分享平台!

微传网博客

网站资源均来自网络,如有侵权,请联系客服删除!

 网站客服QQ:80879498  会员QQ群:727456886

copyright@ 2018-2028 微传网络工作室版权所有

     经营许可证编号:冀ICP备18006529号-1 ,公安局备案号:13028102000124

收起
展开