非协调元解椭圆主程的一种新格式.pdfVIP

引言 科学、技术、工程中许多数学建模问题需要用微分方程来描述，许多近代自然科学都 把微分方程作为基本方程．微分方程研究的来源极广，历史久远Ⅳe伽￡帆和Lei轨iz创造微分 和积分运算时就指出了它们的互逆性．当人们用微积分学去研究天文学，几何学，力学，物理 学等等问题时提出大量微分方程和微分方程组的问题．进入20世纪以后，大量的边缘科学 诸如电磁流体力学，动力气象学，海洋动力学，化学流体力学，地下水动力学等等蓬勃发展，提 出更多的新型的微分方程问题．微分方程的求解问题成为当时数学家们的重要任务．初期 人们把求通解作为微分方程的主要目标，如果我们知道微分方程通解的表达式，问题对应的 特解就迎刃而解．我们还可以了解对某些参数的依赖情况，对解进行其他研究．后来人们逐 渐放弃了求通解．首先．能求得通解的方程显然是很少的．目前人们已经对许多微分方程如 线性方程，恰当方程，常系数二阶线性方程和方程组等求出通解，但是大多数微分方程无法 求出通解，如黎卡提方程： 筹=p(z)可2+g(z)∥+7’(z)． 时，人们需要的不是方程的任一解，而是求得满足某些补充条件的特解．Cauchy认为这是人 们放弃求通解的最重要的和决定性的原因．但是绝大多数微分方程(特别是偏微分方程)的 定解问题很难求出解析解，因此由数学和计算机紧密结合的一门工具性、方法性、边缘交 叉性的新学科应运而生一科学与工程计算．如今大量从事力学、物理学、天文学的研究人 员．从事电子、机械、动力、航空、地质勘探等的工程技术人员也要把这么学科作为自己领 域的主要研究手段． 有限元作为一种求微分方程数值解的主要方法，对数值方法的理论和实践带了巨大的 冲击．有限元方法求解微分方程定界问题可以归结为求一个等价的变分问题，有限元法的步 骤可以归纳为： (1)把定解问题转化为变分形式； (2)按求解区域选择单元形状； (3)构造基函数：形成有限元方程并求解； (4)分析收敛性及误差． wi2sDn元是工程计算中常用的一种有限元，与双线性协调元相比它不仅可求出问题 解的函数值，还可以求出解的二阶导数值，这一点在多尺度渐近展开理论分析中显得非常 重要【lJ．wi2sD扎元用的是分片二次多项式，用传统方法求解二阶椭圆问题，收敛阶只能达 到D(危)．文献f31给出一个反例说明这个收敛阶是最优的，不能被改进．借鉴文献[6]本文改 进了离散双线性型并在双线性型上添加了稳定项后使收敛阶提高到D(九2)，另一方面本文 使用的范数要比传统方法中的范数大，更显示出本文方法的高收敛性，数值实验进一步说 明了理论分析的正确性． 四阶椭圆边值问题在弹性力学等领域中有着非常重要的应用．冯康，石钟慈等在f91中 将线性弹性力学的基础，组合弹性结构的数学理论和弹性结构问题的有限元方法有机结 合起来，着重弹性结构问题的数学提法的精确性，详细介绍了其有限元方法．众所周知，四 阶问题通常的有限元离散要求形函数空间是H2(Q)的一个子空间，这就意味着我们需要整 体C·的分片光滑函数．但是，构造这样的C1函数要么非常困难，要么会使得单元相当复杂．最 常用的方法是采用非协调元，此时有限元空间中的函数不需要满足C1连续性．关于四阶板弯 人在f101中也分析了其收敛性，得到了常规的收敛阶0(^)，本文在双线性型后加了稳定项后 使误差阶在较大范数下依然能达到0(^)． 2 第一章弟一早 耿亩刘以预备知识 1．1泛函分析与sobolev空间相关知识【7][8】 1．1．1 L播Milgram引理 且满足： (1)存在Mo，使得In(z，可)J≤M恻川训，Vz，可∈日； (2)存在6o，使得ln(z，z)I≥洲z悒比∈日． 则存在唯一的z∈日，使得 n(z：y)=2(z)，V可∈日． 1．1．2线性赋范空间的模等价 定义1．1．2．1：设¨lI，和”l|z为线性空间x上的两个范数：我们称”112比¨II·强是指： j IIz。112-÷o号IIznIllo：(佗-÷oo)． 如果”IIz比¨II。强且”Il。比”II。强，则称