论文部分内容阅读
四边形单元与三角形单元相比,具有网格简单、刚度矩阵带宽小、计算量小等优点,而矩形单元对求解区域的边界的近似有一定的局限性,因此将任意四边形单元用于求解问题更有理论意义和实际价值。
G.Acosta,RDurSn在文献中指出弱角条件RDP(N,ψ)似乎是最弱的网格条件,并证明了在此条件下凸四边形单元K上Q1等参元在HI模下的最优插值误差估计,即RDP(N,ψ)是凸四边形单元K上Q1等参元在H1模下取得最优插值误差估计的充分条件.但是否是必要条件呢?他们将它作为一个悬而未决的问题.Z.Zhang在f691中推测它是必要的.然而,本文通过一个反例证明了弱角条件RDP(N,ψ)不是凸四边形单元K上Q1等参元在H1模下取得最优插值误差估计的必要条件,从而圆满解决了这一问题.
对于Stokes问题,本文构造了两个新的低阶非协调任意四边形单元.新单元与以前用于讨论相同问题的单元相比兼有构造简单;使用方便;形函数空间的基函数的次数低;能较好的逼近求解区域边界等优势.由于新单元不含任何协调部分,误差估计比较困难,通过采用更加简洁的证明方法得到了最优误差估计.特别指出的是,其中一个单元在矩形网格下,还是一个Locking-free元,可用于平面弹性问题.另外,本文通过引入变网格思想又将该单元用于非定常Stokes问题得到了最优误差估计.
众所周知,在构造有限元时,单元自由度和形函数空间应当匹配,一方面使有限元空间中的函数在跨越单元边界时,具有某种积分意义下的连续性,同时自由度又是有限元离散方程的未知量,应取得简单,使整体自由度尽量的少.同时具备这两个性质的自由度的选取往往是很困难的.陈绍春教授和石钟慈院士提出了构造有限元的双参数法恰好解决了这一矛盾,并由此构造了许多有价值的双参数有限元(见[54,70,71]等).但有关双参数有限元方法用于变分不等式问题的研究还不多.因此本文将一个12参双参数矩形板元用于四阶位移障碍变分不等式问题,同时通过引入与已有文献不同的新技巧给出了与常规有限元相同的最优误差估计.值得一提的是,我们的结论对几乎所有已知的非协调元都成立。
各向异性有限元方法的研究是目前有限元领域内的亮点和难点之一。由于在各向异性网格剖分下,传统的Sobolev插值理论不能直接利用.虽然T.Apel[38]给出了一个检验单元能否应用于各向异性网格的标准,但应用起来很不方便.陈绍春教授和石东洋教授概括了T.Apel的结果并提出了一个便于操作的插值定理并将其广泛应用于二阶问题和四阶问题中的许多著名的单元(见[20,64,65,66,67]等)。在实际工程计算中,有限元的超逼近和超收敛分析占有非常重要的地位,她一直是数值分析家们研究的热点问题之一.以林群院士为首的科研小组利用积分恒等式技巧对许多单元作了研究并得到了很好的结果。但几乎所有关于这些单元的超收敛性结果都要求网格满足正则性假设和拟一致假设㈣.有关各向异性网格下非协调元的超逼近和整体超收敛分析的文献还不多见。
本文首先将一个低阶的具有各向异性特征的非协调四边形单元(仅有三个独立的自由度)用于求解Stokes方程并给出了一个组合加罚方法,它不仅可以在罚因子较大时得到与传统有限元加罚法中罚因子较小时同样的精度,而且对同一级别的罚因子其收敛阶是传统方法的二倍,这是Falk等人提出的加罚外推法所不能得到的结果,该方法还可以应用于Navier-Stokes方程.
其次,我们将各向异性有限元用于耦合问题,利用具有各向异性特征的双线性元和双二次元构造了一个非协调mortar元,利用积分恒等式技巧得到了与传统方法完全相同的超逼近结果.需要指出的是,本文所提供的思想具有一般性,可用于其它类型有限元的研究.
最后作者将改进的五节点矩形元用于二阶问题,利用其自身独有的特点,在各向异性网格下研究了它的超逼近和超收敛性质.数值检验的结果与研究的理论分析是相吻合的.这对进一步发展二阶问题数值解的后验误差估计方法和设计自适应算法是很有帮助的.