两类变分不等式问题的有限元方法.docxVIP

第一章 概论 1．1常用记号及主要引理 我们在所考虑的有限元逼近问题中，剖分单元一般用T表示，T可为三角形 单元(或矩形单元)，d，为T的顶点，，。为其相应的三条边(或四条边)， i=1,2，3(或扛1,2，3，4)。设_，s．分别表示边，。上单位外法向导数及单位切向导 数。如图所示： ·q=v(q)表示单元顶点q处的函数值，卜；‰为，，边上外法向导数积分 ’‘j rm 平均值。 ·对于三角形单元，设^．：生型△．．掣表示三角形的面积坐标，z：1，2，3 其中 也 Ⅵ  = q = 一 一 b 睁坼 X = h盼 匕述 三 式中 ．u ％引雎 ● Z 3 撒％ ， 由解析 几何 知 kH 1l △l=ltJ，：1l l弓儿1l 2wl 4-W2+w3 =口lb2一“2b{ =a2b3一a3b2 2口3bl一口1b 并记A。=lal ‘对于矩形单元不妨设中点在(o，o)点，边长分别为 =吲=2a，II：I=川=2b 1．2二阶与四阶变分不等式问题的已有结果 对于二阶变分不等式问题的逼近，自70年代起有以下的研究结果， 文献[13】利用协调线形元分别在齐次与非齐次的边界条件下获得了D(向)的 误差阶。 [14]利用协调二次元在自由边界具有有限长度的假设以及z=0的条件下， 获得了误差阶o(^i)v巧0。 [41，【5】]引入非协调元对二阶变分不等式问题进行逼近，先后采用非协调线 形元以及Wilson元，于[2】中得到了与[1 3]中相同的结果。同时，[1]$fl用协调二 次元，在没有自由边界具有有限长度的假设下，获得了与[14]相同的误差阶 !一J o(h：1VJ0。 关于四阶变分不等式问题的有限元方法，到目前为止研究工作还相当少 d [23]对固支情形曲率障碍问题构造了Morley元逼近，给出了收敛性分析， 但没给出误差估计。 【9]对简支情形曲率障碍问题构造了Clough Tother协调元逼近，得到了 !一土 o(h 2，)的误差估计。 【3]对简支情形曲率障碍问题，在解满足一定的光滑性的条件下给出了 Morley非协调元逼近的误差估计c矗(H湘+19l。+≠,lt／tf。) 1．3本文的主要结果 x,t于-次元逼近，本文第三章在保持[1】的误差阶的前提下，将适用范围由 [1]中的凸多边形区域扩展至具有光滑边界的凸区域，主要结论为定理2．2 关于四阶变分不等式问题，，本文对[3]的逼近格式进行修正，给出了新的误 差估计结果，同时利用能量模的正交性，证明所得结果对另一个8一自由度的矩形 元同样适用。主要结论为定理3 2 1及定理3 3 1 第二章二阶变分不等式的逼近问题 2．1基本知识与已有结果 设Q是R 2内的凸区域，f∈r，z∈日2(n)，g是H 2(Q)内的一个函数在边 界上的迹，z≤g在aQ上。 考虑下述二阶问题： 勘毗 足q 懈叫 Ⅳ K (2．1．1) K=p∈JVj(n)：v≥zaP在Q中，v=g在弛上} (2 I 2) 其中 订(址V)=廖”Vvdx， (／，V)=L，‘vdr  (213) 由文章【14]知问题(2 1 1)等价于下述微分问题 一△“=，，在n+=仁∈n：(r)’z(r)冲， 一Au≥，，在Q。=扛∈Q：¨(x)=z(x))[扣， (2 l 4) “2 z，在Q中，n=g在aQ上 关于问题(2 1 1)的正则性，我们有如下引理 引理2．I r圳耶61如果f∈r nB矿(Q)，(晶z)∈c3每)，n的边界足够光滑且在Q 的边界上z≤g，则问题(1．1)的解“满足 zf∈W s-’(Q)，1p∞，2+三一万52+土 (2．1 5) P P 考虑问题(2 1 1)的协调有限元逼近，令L是QNiIiNN#，％=ut；LT，T为 6 剖分单元，吒cL2(n。)是范数为 。的有限元空间．K。c K是世的有限元逼近且 为吒的闭凸集。那么问题(2 1 1)的有限元逼近如下 J求％∈K。，使得 1口(zf。，v。一z，。)≥(，，v。——甜。)vv。∈足。  (216’ 其中 口@一，V一一“一)=VⅣ。·Vvha6c 关于问题(2 I 1)，[1】利用二次元逼近在凸多边形区域己得定理2 1 定理2．1设问题(2．1 1)的解具有正则性：”∈W”(n)v1‘p。，52+上， P l抄∈H‘(Q)n￡●(Q)，z∈W”n∥2。(Q)，“。为(2 1 6)解，则 II,,-．。ll，。≤西_31 6-，vJo (2 l 7) 为以后的证明需要，引入下面引理 引理2．1‘”1 若对所有T∈Th，v∈B(T)，]x。∈T使得v(r1)：0。则下面不等式成立 Ilql：-(hTlvl。 (2 1．8) 引理2．2设f，z，。分别为(2 1 1)与(2 1 6)的解，则vv。∈K。，