工程电磁场-第7章-二维泊松方程的有限元法.ppt_第1页
工程电磁场-第7章-二维泊松方程的有限元法.ppt_第2页
工程电磁场-第7章-二维泊松方程的有限元法.ppt_第3页
工程电磁场-第7章-二维泊松方程的有限元法.ppt_第4页
工程电磁场-第7章-二维泊松方程的有限元法.ppt_第5页
已阅读5页，还剩36页未读， 继续免费阅读

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

7.3二维泊松方程的有限元法,以二维静电场泊松方程的求解为例。,目标：依据加权余量法，利用分域基，建立离散的代数方程组，即确定系数Kij和bi。,场域离散二维问题常使用三角形单元离散，便于处理复杂的场域形状，容易实现。,单元：互不重叠，覆盖全部场域；每个单元内介质是单一、均匀的。节点：网格的交点，待求变量的设置点。,该步骤需要记录的信息：节点编号、节点坐标节点属性(激励源、是否边界等)单元编号单元节点编号单元介质,三角形单元内的基函数设三角形三个顶点处待求函数值分别为u1,u2,u3。如果单元足够小，可以采用线性近似，将单元内任意p点的u(x,y)表示为,代入三个顶点的坐标和函数值，可以解出a、b、c。得到,单元节点的编号按逆时针方向排列！,其中，,记住我们的任务寻找基函数,对比,可得,基函数Ni常被称为插值函数或者形状函数，具有以下性质：（1）是插值的；（2）（3）在相邻单元的公共边界上，Ni是连续的，从而通过Ni构造的逼近函数也是连续的。,在积分中，对于确定的i，j的有效取值为i本身以及与节点i相联的周围节点，积分的有效区域为以i、j为公共节点的所有三角形单元，在这些单元中Ni、Nj才有交叠。,计算系数阵,这些积分可以分单元进行。例如对右图所示的局部编码，K01、K00以及b0的计算公式为：,计算系数阵,以下把单元e的贡献记为,这样，就有,每个或的计算都在具体的单元内单独考虑（称为单元分析）。,由于单元很小，做单元分析时通常可以取f(e)为常数值（可以认为等于三个顶点上的平均值）。因此,右端项元素：,上述以节点为序的分析过程对于有限元原理的说明是易于理解的。而在实际编程中，更有效率的是以单元为序，逐个计算单元系数阵K(e)，然后合成整体系数阵K。单元系数阵K(e)定义为设i,j,m是节点的整体编号，元素Kij在整体矩阵中的实际位置是第i行、j列；因此必须合成到整体矩阵的第i行、j列元素上。,单元矩阵：,整体矩阵合成：,第一类边界条件（强加边界条件）,第一类边界节点是指边界上函数值已知。因此处理方法是，合成整体系数阵之后，将该节点所在行的主元素
人人文库> 全部分类> 教育资料 > 课件下载

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论