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

第一章 概论 1．1常用记号及主要引理 我们在所考虑的有限元逼近问题中，剖分单元一般用T表示，T可为三角形 单元(或矩形单元)，d，为T的顶点，，。为其相应的三条边(或四条边)， i=1,2，3(或扛1,2，3，4)。设_，s．分别表示边，。上单位外法向导数及单位切向导 数。如图所示： ·q=v(q)表示单元顶点q处的函数值，卜；‰为，，边上外法向导数积分 ’‘j rm 平均值。 ’ 7 7 ·对于三角形单元，设^．：生型掣表示三角形的面积坐标，z：1，2，3 △．． 其中 = = 一 也 b q 一 坼 X 睁 Ⅵ = h盼 匕 三 ● 3 ， 述 式中 ．u ％引雎 Z撒％ 由解析几何知 kH1l △l=ltJ，：1l l弓儿1l 2wl4-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】中得到了与[13]中相同的结果。同时，[1]$fl用协调二 次元，在没有自由边界具有有限长度的假设下，获得了与[14]相同的误差阶 !一J o(h：1VJ0。 关于四阶变分不等式问题的有限元方法，到目前为止研究工作还相当少 d [23]对固支情形曲率障碍问题构造了Morley元逼近，给出了收敛性分析， 但没给出误差估计。 Tother协调元逼近，得到了 【9]对简支情形曲率障碍问题构造了Clough !一土 2，)的误差估计。 o(h 【3]对简支情形曲率障碍问题，在解满足一定的光滑性的条件下给出了 Morley非协调元逼近的误差估计c矗(H湘+19l。+≠,lt／tf。) 1．3本文的主要结果 x,t于-次元逼近，本文第三章在保持[1】的误差阶的前提下，将适用范围由 [1]中的凸多边形区域扩展至具有光滑边界的凸区域，主要结论为定理2．2 关于四阶变分不等式问题，，本文对[3]的逼近格式进行修正，给出了新的误 差估计结果，同时利用能量模的正交性，证明所得结果对另一个8一自由度的矩形 2 31 元同样适用。主要结论为定理3 1及定理3 第二章二阶变分不等式的逼近问题 2．1基本知识与已有结果 设Q是R 2(Q)内的一个函数在边 2内的凸区域，f∈r，z∈日2(n)，g是H 界上的迹，z≤g在aQ上。 考虑下述二阶问题： (2．1．1)