二维三角形线性插值有限元法的原理与应用.docxVIP

电磁场数值计算结课大作业 电磁场数值计算结课大作业 PAGE PAGE # 二维三角形单元的有限元程序编写及应用 葛青宇1182201150 研电1803班 一、二维有限元理论介绍 1、电磁场边值问题 泊松方程是电磁场位函数⑴ 泊松方程是电磁场位函数 ⑴ (2)W 2u = f (2) .1 =U0 ；u a — :n =bu c -2,3 其中：(1)式称为泊松方程， 为电磁场的位函数， 为表征空间介质材料的常数， 表示 产生电磁场位函数的源；(2)式表示第一类边界条件， 为边界上已知的位函数；(3)式为 第三类边界条件， 和 为已知的边界条件系数。 2、二维伽辽金加权余量方程 以插值方式构造近似函数， 的近似解表示为 nn U f M nUn (4) 其中，基函数为M1 其中，基函数为M1,M2, 权函数与单元形状函数的关系, ■ ■ , Mnn ,相应待定常数为 U1,U2, ■■■ ,Unn根据基函数和 在有限元法中一个单元内基函数和权函数的值等于相应的形 状函数的值，表不为在单兀 e中，M n =Ni, M m =Nj。 e J e, j 采用伽辽金加权余量形式， 设一组与基函数相同的权函数 M1,M2, ■ ,Mnn ,将权函 数代入伽辽金加权余量方程，得一组方程 TOC \o 1-5 \h \z -.Mm(a—u)d,J = Mmfd 1 ⑸ Q Q 其中，m=1, 2,…，nn,对应节点的总体编号。对(4)式应用格林定理，并将第二、 三类边界条件代入得 a Mm u ud Mm(bu c)d- = Mmfd1 1 (6) Q T Q 其中m同上。 将不含未知函数u的项移至方程右端，得 av Mm ud「QbMmud】=Mmfd，「[cMmd] 将由基函数线性组合构成的近似函数代入方程，得. nn nn TOC \o 1-5 \h \z [atMm (MnUn)]d「J-RbMmV (MnUn)d： = (Mm f )d，J 「(CMm)d： (8) 「1 nl n鼻 m J 将近似函数的梯度运算表示出来，得 nn nn ⑶ MmrGMn)un]}d [bMm (Mnun)]d[ (9)(10)；-1 ng n 1 (9) (10) =(Mmf)dc [(CMm)d 权函数与求和无关，可以将其移到求和号之内，得 nn nn {a「GMmJMn)un]}d「-[[b (MmMn)un]d： TOC \o 1-5 \h \z n 1 n 4 =(Mmf)d? [(cMm)d： 将未知节点变量从积分中提出得： nn nn _ -(a^Mm ； Mn)d-un eMmMnWun = (Mmf)d「「一(cMm)d1 (11) n n」L 可以看出代数方程如下所示： Km,lUl - Km,2U2 U-Km,nUn L_l +Km,nnUnn = % (12) 对比式(10),可得： (13)(14)(15)(16)Km,n=⑻ Mm.、Mn)dC-D(bMmM (13) (14) (15) (16) Rm = . (Mmf)d，「[](cMm)d： 将(12)、(13)式化作单元上积分之和得形式，得： ne neb Km,n (a^ Mm ； Mn)dJ . (bMmMn)d- e^E eT Eb ne neb Rm C .(Mmf)d「「八.(cMm)d： e 丑 E e d Eb ne和neb分别是场域单元数和边界单元数， E和Eb分别是单个场域单元积分区域和单个 边界单元积分区域。 将基函数和权函数更换为单元形状函数，在单元上将基函数和权函数用双层下标标记： ne neb Km,n=￡ ((a^Mmei ^Mne,j)dQ-Z j (bM叫加「。」)d「 (17) e4 E e=1 Eb ne neb Km,n=￡ KaVMmei NMne,j)dG—￡ J(bM叫jMne,j)d「 (18) e 4 E ejEb 给定m和n的值，根据单元节点局部编号与整体编号的对应关系， 寻找单元上的i和j。 单元形状函数如下所示： ne neb Km,n 八 ⑻、Ni J Nj2-q (bNiNj)d： (19) e=1 E e』Eb ne neb Rm c .(Nif)d「「八 (cNi)d] (20) e=E e* 得到得方程组以矩阵形式表示为： Ku = R (21) 3、二维三角形有限元矩阵计算 在二维三角形单元中 (1)场域单元系数矩阵 ： Ke,i,j = .(a【Ni .\Nj)dj E y \ \ 11yj书 _yj^2 aL(yi 由- y42) (x 也一X4O」| xj — xj 书 (2 )2 (22) a ， —(yi1-yi2)(yj1-yj2) (^ ^-xi1)(Xj ^X