纯应力边界条件下的平面弹性问题的一个二阶非协调矩形元

1引言对于几乎不可压材料,当材料的Lame常数λ→∞时,通常的低阶协调有限元的解不再收敛到原问题的解或达不到最优收敛阶,这就是弹性材料的Locking现象[1,2,3,5,4].为了克服Locking问题,目前已提出许多有效的方法.例如:采用高次元[6,7,8]、非协调元[4,9,10,11]和混合有限元方法[1,2,12,13]等.     

各向异性网格下Bogner-Fox-Schmit元的超收敛

1引言 有限元超收敛的研究已有三十多年的历史,至今为止已取得了丰富的成果,可见[3],[18],[10],[6],[5]以及[17].1981年,陈传淼(见[2]345-372页)考虑了四阶板问题有限元解的超收敛性,得到了高一阶的超收敛结果.1995年,林群和罗平[8]用积分恒等式技巧再次研究这个问题,在均匀矩形网格的条件下,得到了更好的结论,有限元解与有限元插值函数之间的误差在H2范数下,有高二阶的超收敛.     

无限元多重网格算法

在求偏微分方程数值解时,往往需要解一个规模很大的代数方程组,而多重网格是一种十分有效的迭代方法.大量数值试验证明,它具有很高的收敛速度.理论分析表明,这种迭代法的收敛速度并不随网格的加密而降低,这一突出优点是其它迭代方法望尘莫及的. 在使用有限元多重网格算法时,如果区域边界的角点使解具有奇性,理论分析会遇到     

解一类非线性方程组的球形迭代法

对一类非线性方程组,本文给出一种球形迭代解法。它的基本思想是:以空间某一固定点为球心,给出一球形区域为方程组解的初估计范围,以球中任一点为迭代初值,按某一格式迭代,当迭代解超出这一区域,则将这球形区域的半径扩大,同时重新从迭代初值出发迭代,我们将证明,当球形区域最终包含方程组解时,迭代解至多有限次超出球形区域,直至收敛到方程组的解。这种解法具有大范围收敛性,同时允许迭代格式中带有误差项,适合这种解法的非线性方程组较[1],[2],[3]中的更广,而迭代格式中的误差项又比[4]中更一般。     

线弹性问题在交错网格上的有限差分算法

1引言对线弹性材料的形变分析在工程实践中有着广泛应用,而数值方法是求解该问题的一种重要方法[1].其中,低阶有限元方法常用来处理可压缩线弹性问题[2,3].但是当材料几乎不可压缩时(拉梅常数λ→∞),普通有限元方法的收敛性会被破坏[4].这种现象被称为”locking”.基于交错网格的有限差分方法是科学计算和数值分析领域中的一种有效方法,被广泛应用于流体问题、多孔介质流问题以及麦克斯韦方程组等.     

基于网格变动的油水驱动半正定问题迎风混合元方法

在网格随时间变动的有限元空间上研究了不可压缩的两相渗流驱动问题.分别对饱和度方程扩散矩阵正定和半正定的情形,提出了基于网格变动的迎风混合元方法混合元逼进压力方程,饱和度方程的对流项采用迎风格式来处理,扩散项则采用推广的混合元来逼进.在网格任意变动的情形下得到几乎最优的误差估计;对正定问题的格式进行改进,即在两个网格之间投影变化时采取近似解的线性构造,可以得到与固定网格时相同的最优收敛阶.     

非协调Wilson有限元的多重网格方法

§1.引言 非协调Wilson有限元[1—3]对解弹性力学方程有实用价值,在工程上有用。本文分析Wilson元的多重网格法,给出用多重网格方法求得的近似解按L~2模和能量模的最佳收敛阶误差估计。对于W-循环,可以证明其计算量与离散空间的维数为同一量级O(N_k)。 考虑二阶椭圆Dirchlet边值问题:     

一阶拟线性非齐次偏微分方程存在整体光滑解的条件

在文献[1],[2]中讨论了一阶拟线性齐次偏微分方程 Cauchy 问题(1)(2)关于整体光滑解的存在性问题.文献[1]得到了λ_i=λ_i(u)时 Cauchy 问题(1)、(2)存在整体光滑解的充要条件;文献[2]进而得到了λ_i=λ_i(t,x,u)时 Cauchy 问题(1)、(2)存在整体光滑解的充要条件。本文将用[1]、[2]的思想方法,讨论一阶拟线性非齐次偏微分方程 Cauchy 问题     

n维矩形域上椭圆问题有限元单方向外推

1 引言 Richardson外推应用于椭圆偏微方程边值问题有限元法始于1978年(见[1],并于1983年在理论研究方面取得突破性进展(见[2]).自那以后有限元外推得到迅速发展,成为一个富于竞争的国际性研究课题(见[3],[4],[5]及其所列参考文献).但是通常的有限元外推需同时在每一个方向上分半加密网格,因此,对n维问题,细网格的结点数是粗网格的2~n倍,结果当n较大时(高维问题),细网格上的计算工作量十分庞大.为了克服这个缺点,发展了有限元单方向外推.对Poisson方程边值问题,[6]研究了2维矩形域上双线性有限元解的单方向外推,[7]研究了3维矩形域上三线性有限元单方向外推必须的插值渐近展开式,[8]研究了n维矩形域上n线性有限元解的区域分裂外推.本文旨在研究n维矩形域上Poisson方程边值问题及其对应的本征值问题n线性有限元解的单方向外推.始终假设本文出现的函数u是连续的.     

半线性椭圆问题Petrov-Galerkin逼近及亏量迭代

本文考虑了二阶半线性椭圆问题的Petrov-Galerkin逼近格式,用双二次多项式空间作为形函数空间,用双线性多项式空间作为试探函数空间,证明了此逼近格式与标准的二次有限元逼近格式有同样的收敛阶.并且根据插值算子的逼近性质,进一步证明了半线性有限元解的亏量迭代序列收敛到Petrov-Galerkin解.     

Stabilized Nonconforming Mixed Finite Element Method for Linear Elasticity on Rectangular or Cubic Meshes

Based on the primal mixed variational formulation, a stabilized nonconforming mixed finite element method is proposed for the linear elasticity on rectangular and cubic meshes. Two kinds of penalty terms are introduced in the stabilized mixed formulation, which are the jump penalty term for the displacement and the divergence penalty term for the stress. We use the classical nonconforming rectangular and cubic elements for the displacement and the discontinuous piecewise polynomial space for the stress, where the discrete space for stress are carefully chosen to guarantee the well-posedness of discrete formulation. The stabilized mixed method is locking-free. The optimal convergence order is derived in the $L^2$-norm for stress and in the broken $H^1$-norm and $L^2$-norm for displacement. A numerical test is carried out to verify the optimal convergence of the stabilized method.     

MULTIGRID METHODS FOR MORLEY ELEMENT ON NONNESTED MESHES

1.IntroductionWeconsidersomemultigridalgorithmsforthebiharmonicequationdiscretizedbyMoneyelementonnonnestedmeshes.TOdefineamultigridalgorithm,certainintergridtransferoperatorhastobeconstructed.Throughtakingtheaveragesofthenodalvariables,weconstructanintergridtransferoperatorforMoneyelementonnonnestedmeshesthatsatisfiesacertainstableapproximationpropertywhichplaysakeyroleinmultigridmethodsfornonconformingplateelementsonnonnestedmeshes.Theso--calledregularity-approximaticnassurnptionisestablis…     

A NEW NONCONFORMING MIXED FINITE ELEMENT SCHEME FOR THE STATIONARY NAVIER-STOKES EQUATIONS

In this article, a new stable nonconforming mixed finite element scheme is proposed for the stationary Navier-Stokes equations, in which a new low order CrouzeixRaviart type nonconforming rectangular element is taken for approximating space for the velocity and the piecewise constant element for the pressure. The optimal order error estimates for the approximation of both the velocity and the pressure in L2-norm are established, as well as one in broken H1-norm for the velocity. Numerical experiments are given which are consistent with our theoretical analysis.     

PRECONDITIONING MIXED ELEMENT METHOD FOR NONSELFADJOINT AND INDEFINITE SECOND ORDER ELLIPTIC PROBLEMS*

In this paper, an equivalence between mixed element method and nonconforming element method for nonselfad joint and indefinite second order elliptic problems is established without using any bubble functions. It is proved that the H~1-condition number of preconditioned operator B_h~(-1)A_h is uniformly bounded and its B_h-singular values cluster in a positive finite interval, where A_h is the equivalent nonconforming element discretization of nonselfad joint and indefinite second order elliptic operator A, B_h is usual noncon forming element discretization of selfadjoint and positive definite second order elliptic operator B. Finally a simple V-cycle multigrid implementation of B_h~(-1) is given.     

一个改进的Reissner-Mindlin矩形元

本文提出了一个改进的Reissner-Mindlin矩形非协调元方法：旋度用连续双线性元逼近，横向位移用旋转矩形非协调元逼近，而作为中间变量的剪切力用增广的分片常数元逼近，我们证明：该方法具有关于板厚一致稳定性和一致最优收敛性。     

非嵌套有限元的W循环多重网格方法

本文讨论了下述情形：１非嵌套网格；２曲边有限元；３非协调元；４拟协调元；５有限元的型函数有特殊性质，都能导致非嵌套的有限元空间．对一个包括上述情形的问题给出了非嵌套有限元的Ｗ循环多重网格方法，并证明了它的收敛性。     

EXPLICIT ERROR ESTIMATES FOR MIXED AND NONCONFORMING FINITE ELEMENTS

In this paper, we study the explicit expressions of the constants in the error estimates of the lowest order mixed and nonconforming finite element methods. We start with an explicit relation between the error constant of the lowest order Raviart-Thomas interpolation error and the geometric characters of the triangle. This gives an explicit error constant of the lowest order mixed finite element method. Furthermore, similar results can be ex- tended to the nonconforming P1 scheme based on its close connection with the lowest order Raviart-Thomas method. Meanwhile, such explicit a priori error estimates can be used as computable error bounds, which are also consistent with the maximal angle condition for the optimal error estimates of mixed and nonconforming finite element methods.     

MULTIGRID FOR THE MORTAR ELEMENT METHOD WITH LOCALLY P_1 NONCONFORMING ELEMENTS

In this paper we study the theoretical properties of multigrid algorithm for discretization of the Poisson equation in 2D using a mortar element method under the assumption that the triangulations on every subdomain are uniform. We prove the convergence of the W-cycle with a sufficiently large number of smoothing steps. The     

Multigrid for the mortar element method for P1 nonconforming element

In this paper, a multigrid algorithm is presented for the mortar element method for P1 nonconforming element. Based on the theory developed by Bramble, Pasciak, Xu in [5], we prove that the W-cycle multigrid is optimal, i.e. the convergence rate is independent of the mesh size and mesh level. Meanwhile, a variable V-cycle multigrid preconditioner is constructed, which results in a preconditioned system with uniformly bounded condition number. Received May 11, 1999 / Revised version received April 1, 2000 / Published online October 16, 2000     

MULTIGRID METHODS FOR THE GENERALIZED STOKES EQUATIONS BASED ON MIXED FINITE ELEMENT METHODS

1. IntroductionIn this paper, we consider the fOllowing generalized stationary Stokes equations:where fl is a bounded convex domain in R', u represents the velocity of fluid, p its pressure; Fand G are external fOrce and source terms. Note that the source…