Sobolev方程的两类非协调元高精度分析.pdf

Abstract Inthis orderWilsonelementand dementare to paper，thehigher Carey applied Sobolev ofthetwoelementsarediscussedunder equation．Higheraccuracyanalysis of of the structureelementsand80menewtech- anisotropicmeshes．Byusing particular errorestimateand area8sameasresultsderivedOil optimal niques，the superclose，which meshesare isusedto regular derived．Then，theinterpolation technique postprocessing obtainthe globalsuperconvergeilce． Words：Sobolev finite meshes Key equation；nonconformingelement；anisotropie superconvergence；postprocessingtechnique． 11 前言 基于三角形网格剖分的Dirichlet问题的分片线性逼近，在我国，计算数学家冯 康先生首先独立于西方发明了这种方法【2，3】，有限元方法的实质是根据变分原 理用有限维空间的离散解去逼近无穷维空间的连续解，有限维空间的构造， Sobolev空间插值理论以及微分方程正则性理论都是这种方法能够实现的理论 前提．近年来，随着计算机的蓬勃发展，有限元方法的研究也蒸蒸日上．从20 世纪60年代至今，经过40多年的研究和发展，有限元方法已成为一门理论完 善、应用广泛的数值计算方法． 在传统有限元方法的研究中，剖分的正则性假设或拟一致性假设【4，5】(即 警≤c或=h≤c，这里hK为剖分单元K的直径，pK为K的最大内切圆直径， h=maxⅣhⅣ，元=minKh耳，e是一个常数)是传统有限元方法进行理论分析的前 提条件．但最近的一些研究成果f6-17]表明，这种假设对一些窄边有限元和各 向异性元并非是必要的．同时有些问题定义在窄边区域，如果用正则性剖分， 计算量将非常大；另一方面有些问题的解呈各向异性，即沿某个方向解变化非 常剧烈，而沿另外方向解变化平缓，这时若采用各向异性剖分可以用较少的自 理在插值误差分析中已不再适用(见[7-9])．而且对于非协调元来说，其传统边 界估计技巧也不再适用，因为此时对单元的长边，其误差估计中出现的嚣项 在专著中研究了各向异性下Lagrange型协调元的误差分析，并给出了判别一 个单元是否具有各向异性特征的判别定理．从而利用具有各向异性特征的单元 使我们得以用较小的自由度得到同样的最优误差估计效果．但这种方法有时难 以验证．陈绍春，石东洋和赵永成【10】对它进行了改进，给出了一个更为一般 的各向异性判别定理，并将之应用到实际问题中(包括各向异性非协调元)，取 得了许多有意义和价值的成果(见【10，14，15，16，17】)． 首先我们先利用积分恒等式技巧，在各向异性矩形网格之下对Wilson元进 行了分析，通过引入一些技巧得到了与正则网格下相同的最优误差估计及超逼 近性质，进而通过插值后处理的技巧，得到了相应