Stokes问题基于泡函数的简化的稳定化混合元格式的收敛性

利用泡函数导出Stokes问题的两种新的、简化的稳定化混合有限元格式.并证明这些格式与通常带泡函数的稳定化格式具有相同的收敛性,但是自由度可以大大减少.     

二阶椭圆问题的混合有限元法的泡函数稳定性及其后验误差估计

本文讨论二阶椭圆问题的混合有限元逼近的一种泡函数稳定性,并给出其基于简化的稳定化格式的先验误差估计和后验误差估计.该方法较通常的格式(例如,Raviaxt-Thomas方法的同阶格式)节省大量的自由度.     

平面弹性力学问题的混合元法的泡函数稳定性及其后验误差估计

该文讨论平面弹性力学问题的混合元法的泡函数稳定性,并导出基于简化的稳定化格式的一种先验误差估计和后验误差估计．这种简化的稳定化格式较通常的格式节省自由度．     

抛物型方程基于POD方法的时间二阶中心差的时间二阶精度简化有限元格式

本文将特征正交分解(Proper Orthogonal Decomposition,简记为POD)方法应用于抛物型方程通常时间二阶中心差的时间二阶精度有限元格式(简称为通常格式),简化其为一个自由度极少但具有时间二阶精度的有限元格式,并给出简化的时间二阶中心差的时间二阶精度有限元格式(简称为简化格式)解的误差分析.数值...     

二阶椭圆问题的各向异性混合元简化格式和后验误差估计

对二阶椭圆问题构造了一个非常规各向异性Hermite型矩形单元.并基于泡函数对其构造了一种简化的稳定化混合元格式.同时给出了格式的收敛性分析和后验误差估计.     

抛物型方程基于POD方法的时间二阶精度CN有限元降维格式

将特征正交分解(proper orthogonal decomposition, 简记为POD) 方法应用于抛物型方程通常的时间二阶精度Crank-Nicolson (简记为CN) 有限元格式, 简化其为一个自由度极少的时间二阶精度CN 有限元降维格式, 并给出简化的时间二阶精度CN 有限元解的误差分析. 数值例子表明在简化的时间二阶精度CN 有限元解和通常的时间二阶精度CN 有限元解之间的误差足够小的情况下, 简化的时间二阶精度CN 有限元格式能大大地节省自由度, 而且时间步长可以比时间一阶精度的格式取大10 倍, 以至能更快计算到所要时刻数值解, 减少计算机计算过程的截断误差, 提高计算速度和计算精度,从而验证降维时间二阶精度CN 有限元格式用于解类似于抛物型方程的时间依赖方程是很有效的.     

二阶椭圆问题的混合有限元估计

关于二阶椭圆问题，[1]中给出了一种混合有限格式，但所用的自由度太多，给实     

Cahn-Hilliard方程的新型二阶稳定化半隐有限元方法

基于Crank-Nicolson/Adams-Bashforth离散,一种新型的二阶稳定化半隐有限元格式被建立用来求解Cahn-Hilliard方程.在此格式中,通过添加一个新型的二阶稳定项,得到一个满足离散的能量耗散定律的线性系统.空间离散考虑Galerkin有限元方法,从而获得时空全离散格式.算法的稳定性被考虑,同时给出相应的误差估计.理论结果表明,所提出的算法具有二阶精度.最后,数值算例验证所提算法的有效性.     

导算子范数的一个估计

In this paper, a derivation δ_AB mapping into a ideal I of B(H) is considered, when A,B ∈B(H) and I is a norm ideal. If Ran(δ_AB)?I, let δ_AB:B(H)→I denote the induced operator and let λ be the scalar such that A- λ∈I, B-λ∈I, we estimate the norm of δ_AB as follows‖A-λ‖+‖B-λ‖≥‖δ_AB‖≥‖A- λ‖+‖B-λ‖ when W_N(A-λ)∩W_N(λ - B)≠?, where W_N(A- λ) denotes the normalized maximal numerical range and ‖A-λ‖ denotes the norm of A-λ∈I. In particular when I=C_p(lp, we prove that ‖δ_AB‖_p=‖A-λ‖_p+‖B-λ‖_p if and only if ‖A-λ‖=‖A-λ‖_p and W_N(A-λ)∩W_N(λ-B)≠?. At last, some examples show that the estimate as above is exact.     

二阶非线性差分方程有界解振动的充分必要条件

在本文中,我们给出了非线性二阶差分方程△（p_n△y_n）+q_nf（y_n-r_n）=0有界解振动的充分必要条件和比较定理,所得结果推广了文[3,6,7]的相应定理。     

AN ERROR ANALYSIS METHOD SPP-BEAM AND A CONSTRUCTION GUIDELINE OF NONCONFORMING FINITE ELEMENTS FOR FOURTH ORDER ELLIPTIC PROBLEMS

Under two hypotheses of nonconforming finite elements of fourth order elliptic problems,we present a side-patchwise projection based error analysis method(SPP-BEAM for short).Such a method is able to avoid both the regularity condition of exact solutions in the classical error analysis method and the complicated bubble function technique in the recent medius error analysis method.In addition,it is universal enough to admit generalizations.Then,we propose a sufficient condition for these hypotheses by imposing a set of in some sense necessary degrees of freedom of the shape function spaces.As an application,we use the theory to design a P3 second order triangular H2 non-conforming element by enriching two P4 bubble functions and,another P4 second order triangular H2 nonconforming finite element,and a P3 second order tetrahedral H2 non-conforming element by enriching eight P4 bubble functions,adding some more degrees of freedom.     

A new quadratic nonconforming finite element on rectangles

A new quadratic nonconforming finite element on rectangles (or parallelograms) is introduced. The nonconforming element consists of P₂ ⊕ Span{x²y,xy²} on a rectangle and eight degrees of freedom. Our element is essentially of seven degrees of freedom since the degree of freedom associated with the integration on rectangle is essentially of bubble‐function nature. Global basis functions are constructed for both Dirichlet and Neumann type of problems; accordingly the corresponding dimensions are counted. The local and global interpolation operators are defined. Error estimates of optimal order are derived in both broken energy and L²(Ω) norms for second‐order of elliptic problems. Brief numerical results are also shown to confirm the optimality of the presented quadratic nonconforming element. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Advanced finite element formulations for modeling thin piezoelectric structures

This contribution is concerned with mixed finite element formulations for modeling piezoelectric beam and shell structures. Due to the electromechanical coupling, specific deformation modes are joined with electric field components. In bending dominated problems incompatible approximation functions of these fields cause incorrect results. These effects occur in standard finite element formulations, where interpolation functions of lowest order are used. A mixed variational approach is introduced to overcome these problems. The mixed formulation allows for a consistent approximation of the electromechanical coupled problem. It utilizes six independent fields and could be derived from a Hu-Washizu variational principle. Displacements, rotations and the electric potential are employed as nodal degrees of freedom. According to the Timoshenko theory (beam) and the Reissner-Mindlin theory (shell), the formulations account for constant transversal shear strains. To incorporate three dimensional constitutive relations all transversal components of the electric field and the strain field are enriched by mixed finite element interpolations. Thus the complete piezoelectric coupling is appropriately captured. The common assumption of vanishing transversal stress and dielectric displacement components is enforced in an integral sense. Some numerical examples will demonstrate the capability of the presented finite element formulation. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

An asymptotically optimal Schur complement reduction for the Stokes equation

Summary. In this paper we develop an efficient Schur complement method for solving the 2D Stokes equation. As a basic algorithm, we apply a decomposition approach with respect to the trace of the pressure. The alternative stream function-vorticity reduction is also discussed. The original problem is reduced to solving the equivalent boundary (interface) equation with symmetric and positive definite operator in the appropriate trace space. We apply a mixed finite element approximation to the interface operator by iso triangular elements and prove the optimal error estimates in the presence of stabilizing bubble functions. The norm equivalences for the corresponding discrete operators are established. Then we propose an asymptotically optimal compression technique for the related stiffness matrix (in the absence of bubble functions) providing a sparse factorized approximation to the Schur complement. In this case, the algorithm is shown to have an optimal complexity of the order , q = 2 or q = 3, depending on the geometry, where N is the number of degrees of freedom on the interface. In the presence of bubble functions, our method has the complexity arithmetical operations. The Schur complement interface equation is resolved by the PCG iterations with an optimal preconditioner. Received March 20, 1996 / Revised version received October 28, 1997     

Optimal second order rectangular elasticity elements with weakly symmetric stress

We present new second order rectangular mixed finite elements for linear elasticity where the symmetry condition on the stress is imposed weakly with a Lagrange multiplier. The key idea in constructing the new finite elements is enhancing the stress space of the Awanou’s rectangular elements (rectangular Arnold–Falk–Winther elements) using bubble functions. The proposed elements have only 18 and 63 degrees of freedom for the stress in two and three dimensions, respectively, and they achieve the optimal second order convergence of errors for all the unknowns. We also present a new simple a priori error analysis and provide numerical results illustrating our analysis.     

Computation of the Three-Dimensional Stress State in Composite Shell Structures with Mixed Finite Elements

Being able to compute the complete three-dimensional stress state in layered composite shell structures is essential in order to examine complicated interlaminar failure modes such as delamination. We lay out a mixed finite element formulation with independent displacements, rotations, stress resultants and shell strains. A mixed hybrid shell element with 4 nodes and 5 or 6 nodal degrees of freedom is developed, so that the element formulation can also be used for problems with shell intersections. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

p-version finite elements and finite cells for finite strain elastoplastic problems

In this paper, the p-version finite element method and its fictitious domain extension, the finite cell method, are extended to the finite strain J₂ plasticity. High-order shape functions are used for the finite element approximation of volume-preserving plastic dominated deformations. The accuracy and efficiency of p-version elements and cells in the finite plastic strain range are evaluated by the computation of two benchmark problems. It is shown that they provide locking free behavior and simplified meshing. These results are verified in comparison with the results of h-version elements in F-bar formulation. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)