《(完整版)大连理工大学高等数值分析有限元简述》由会员分享，可在线阅读，更多相关《(完整版)大连理工大学高等数值分析有限元简述（9页珍藏版）》请在装配图网上搜索。
1、椭圆与抛物微分方程的有限元法有限元法是与差分法并驾齐驱的一套求解偏微分方程的方法。它的基本想法是，首先把微分方程转化成一种变分方程（微分积分方程），从而降低了对解的光滑性和边值条件的要求；然后，把求解区域划分成有限个单元（有限元）,构造分片光滑函数，这个光滑函数由其在单元顶点上的函数值决定；最后，把这个分片光滑函数带入到上述微分积分方程中去，就得到关于单元顶点函数值的一个线性方程组，解之即得有限元解。与差分法相比，有限元法易于处理边界条件，易于利用分片高次多项式等等来提高逼近精度。空间Hm作为例子，我们将考虑区间I=（。,1）上的微分方程。用L2（I）表示在I上勒贝格平方可积函数的集合，Hm（
2、I）表示本身以及直到m阶的导数都属于L2（I）的函数的集合。我们下面用到的主要是h1（/）。这里所说的导数准确地说是应该是广义导数，对此我们不予详细说明，只需知道比如说，连续的分片线性函数（折线函数）就属于H1（I），其广义导数是分片常数函数。另外，我们还用到空间H1（I）二vGH1（I）,v（0）二0。E（空间=函数集合。）微分方程考虑两点边值问题一(pu)+qu=f,xe(0,1)1)u(0)=02)u(1)=03)其中p,q,f都是区间(0,1)上的光滑函数，q0,并且pp，p是00一个正常数。用h1(I)中任一函数v乘(1)式两端，并在0,1E上积分，得f1-(pu)v+quv-fvd
3、x=00利用分部积分，并注意u(1)=0和v(0)=0，得4)以此代入到(4)得到f1(pu)vdx=-puv|1+J*1puvdx=000f1puvdx0f1(puv+quv-fv)dv=005)为了方便，定义w-vdxa(w,v)=(pw,v)+(qw,v)7)8)则相应于微分方程(1)-(3)的变分方程为：求ueH1(I)满E足a(u,v)=(f,v)VveH1(I)(9)E注意在(9)中不出现二阶导数。可以证明，满足微分方程(1)-(3)的光滑解一定满足变分方程(9)。(9)的解称之为(1)-(3)的广义解，它可能只有一阶导数，因此可能不是(1)-(3)的解；但是如果它在通常意义下二阶
4、可微，则一定也是(1)-(3)的解。另外注意，在变分方程(9)中，我们强制要求广义解u满足边值条件u(0)=0，因而称之为强制(或本质)边界条件；而对边值条件u，(1)=0，则不加要求。但是可以证明，如果广义解u在通常意义下二阶可微，则一定有u，(1)=0，即这个边界条件自然满足。这类边界条件称之为自然边界条件。总之，变分方程(9)不但降低了对解的光滑性的要求，也降低了对边值条件的要求。有限元空间构造有限元法的第一步与差分法一样，也是对求解区间作网格剖分0二xxLx二1。相邻节点x,x之间01ni-1i的小区间I=x,x称为第i个单元，其长度为h二x-x。记ii-1iiii-1h=maxhoi
5、在空间h1(I)中，按如下原则选取有限元空间V:它的元Eh素u(x)满足所谓本质边界条件u(0)二0，在每一单元上是m次hh多项式，并且在每个节点上都是连续的。当m二1时，就得到最简单的线性元，这时每个ueV可表为i=1,2,L,n10)hhx-xx-x，u(x)=iu+u,xeIhhi-1hii其中u=u(x),hiu=u(0)=00h线性元的另外一种表示方法是利用以下具有局部支集的基函数：x-xi,xxxhi-1ii存x-x=1i,xhii+1，在别处xxi+1i=1,2,L,n-111)存x-x1+n,xhn-1n、,在别处xxn12)图2.线性元的基函数显然，任一ueV可以表为hhu(
6、x)二工u申(x)hiii=1有限元方程将变分方程（9）局限在有限元空间上考虑，就得到有限元方程：求有限元解ueV满足hha（u,v）=（f,v）VveV（14）hhhhh注意到u和v都可以表示成（13）形式，容易看出（14）等hh价于如下的线性方程组：求节点上的近似解u,L,u满足15)1n工a（甲,甲）u=（f,甲）,j=1,L,nijiji=1这个线性方程组是三对角的，可以用追赶法求解。可以把微分方程（1）、变分方程（9）和有限元方程（15）比喻为确定“好人”的三种标准：他每一时刻表现都好；每一个人都说他好；一个遴选委员会说他好。误差估计可以证明，微分方程（1）-（3）的解u和有限元方程
7、（14）或（15）的解u之间的误差满足hIIu-uII+hIIuf-ufllChIIuII（16）hh其中c是一个常数；IIII表示gi）范数，定义为|v|=J（v,v）=Jb|v|2dx2,VveU（I）（17）a二维椭圆方程有限元法以二维区域上的Poisson方程第18)一边值问题为例：学学二f(x，y)，(x,y)eG19)Qx2Qy2uI=0r其中g是以r为边界的一个二维区域。利用Green公式，容易推出相应的变分方程：求ueHi(G)满足020)a(u,v)=(f,v)，VveHi(G)0由在边界r上为零且广义偏导数在区域g上勒21)22)其中空间Hi(G)0贝格可积的所有函数组成(
8、w,v)三口wvdxdyGdwdvdwQv、a(w,v)JJ(+)dxdyQxQxdydyG二维区域上最常用的剖分是形如下图的三角剖分我们可以相应地构造三角剖分上的线性元。对内点集合Gh(例如上图中3,6,5这三个点)中每个节点i，定义其基函数甲(x,y)为一个分片线性函数，它在节点i取值为1而在i所有其他节点为0。这样，有限元空间V中任一元素就可以h表示成u(x)丄ucp(x)。把它带入到变分方程(20)便得有限hiiieGh元方程：求g上的近似解u满足hi工a(p,p)u二(f,p),VjeG(23)ijijhieGh高次元可以从两个途径来提高有限元法的精度,一个是加密网格,另一个是利用高
9、次元。例如对于一维问题,可以使用所谓Hermite三次元，它在每一个单元/=卜,x上是一ii-1i24)个三次多项式,由两个端点上的函数值和导数值总共4个参数确定。这时,相应于(16)我们有误差估计IIu-uII+hIIu-ufllCh4fIIu(k)|hhk=0其中u(k)表示k阶导数。对于二维问题也可以使用高次元，但是其定义要稍微复杂一点。抛物方程有限元法考虑一维抛物方程u(x,0)=u(x),0x10u(0,t)=0,lu(1,t)=0,0tT(27)lx其中系数p,q,f都是x和t的已知光滑函数，初值u(x)是x的已0知光滑函数。它的变分方程为：求讥兀t)使得对每一个固定的都有te0,
10、T，其中u(x,t)eH1(I)，并且E28)lu29)(,v)+a(u,v)=(f,v),VveH1(I)ltE(w,v)三J1wvdxdy30)0lwlva(w,v)=(p-，亍)+(qw,v)lxlx抛物方程有限元法的通常做法是在时间方向用差分法在空间方向用有限元法。象在(10)中那样，可以关于变量x构造线性有限元空间v。令时间方向步长为工。若时间方向h用向前差商，空间方向用线性有限元，并记fk=f(x,加)，则有限元方程为：对k=1,L,K=T/t，逐层求uk=uk申(x)eV满足hiihi=1uk+1uk(htv)+a(uk,v)=(fk,v),VveVhhhhhh31)这相当于在每
11、一层要解一个线性方程组：工(q,q)(ut)+Ya(qijtii=1i=1,q)uk=(fk,q),Vj=1,Ljij,n或者稍微整理一下：工(q,q)uk+1iji工(q,q)ukTa(q,q)uk+t(fk,q),Vj=1,Lijiijij如果在时间方向用梯形公式，则类似于（31）得到所谓Crank-Nicolson格式：zUk+1-Uk、Uk+1+Uk、zfk+1+fk（33）（hh,v）+a（hh,v）=（,v）,VveVph2h2hhh习题1设边值条件为u（0）=o，u，（i）=o，步长为h=05。写出相应的线性元的各个基函数，并图示。习题2假设如习题1，并设p=1，q=o，f=1，具体写出线性元有限元方程相应的线性方程组。习题3仿照（32），将Crank-Nicolson格式（33）写成线性方程组形式。习题4将边界条件（3）换成u（1）=o，试推出相应于（14）的有限元方程。
(完整版)大连理工大学高等数值分析有限元简述
