Stabilized Mixed Tetrahedrals with Volume and Area Bubble Functions at Large Deformations

In this contribution, stabilized mixed finite tetrahedral elements are presented in order to avoid volume locking and stress oscillations. Geometrically non-linear elastic problems are addressed. The mixed method of incompatible modes is considered. As a key idea, volume and area bubble functions are used for the method of incompatible modes [1], thus giving it the interpretation of a mixed finite element method with stabilization terms. Concerning non-linear problems these are non-linearly dependent on the current deformation state, however, linearly dependent stabilization terms are used. The approach becomes most attractive for the numerical implementation, since the use of quantities related to the previous Newton iteration step is completely avoided. The variational formulation for the standard two-field method, the method of incompatible modes in finite deformation problems is derived for a hyper elastic Neo-Hookean material. In the representative examples Cook's membrane problem and a block under central pressure illustrate the good performance of the presented approaches compared to existing finite element formulations. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Nonconforming tetrahedral finite elements for fourth order elliptic equations

This paper is devoted to the construction of nonconforming finite elements for the discretization of fourth order elliptic partial differential operators in three spatial dimensions. The newly constructed elements include two nonconforming tetrahedral finite elements and one quasi-conforming tetrahedral element. These elements are proved to be convergent for a model biharmonic equation in three dimensions. In particular, the quasi-conforming tetrahedron element is a modified Zienkiewicz element, while the nonmodified Zienkiewicz element (a tetrahedral element of Hermite type) is proved to be divergent on a special grid.

     

椭圆型方程四面体线元的超逼近与外推

重新讨论了三角线元的积分恒等式,使之适用于三维区域的拟一致四面体元,借此证明了椭圆型方程有限元解梯度有超逼近现象,函数值Richardson外推可以提高精度.     

A weak Galerkin mixed finite element method for the Helmholtz equation with large wave numbers

In this article, a new weak Galerkin mixed finite element method is introduced and analyzed for the Helmholtz equation with large wave numbers. The stability and well‐posedness of the method are established for any wave number k without mesh size constraint. Allowing the use of discontinuous approximating functions makes weak Galerkin mixed method highly flexible in term of little restrictions on approximations and meshes. In the weak Galerkin mixed finite element formulation, approximation functions can be piecewise polynomials with different degrees on different elements and meshes can consist elements with different shapes. Suboptimal order error estimates in both discrete H¹ and L² norms are established for the weak Galerkin mixed finite element solutions. Numerical examples are tested to support the theory.     

三维矩形域上泊松方程四面体线元的超逼近与外推

改进三角元的积分恒等式,使之适用于拟一致四面体元,借此证明了泊松方程四面体线元梯度有超逼近现象,函数值Richardson外推可以提高精度.     

On the stabilization of tetrahedral finite elements using volume and area bubble functions

In order to overcome the oscillatory effects of the bi-linear Galerkin formulation for tetrahedral elements the mixed method of incompatible modes and the mixed method of enhanced strains are reformulated, thus giving both the interpretation of a mixed finite element method with stabilization terms. For nonlinear problems these are nonlinearly dependent on the current deformation state, and therefore are replaced by linearly dependent stabilization terms. The approach becomes most attractive for the numerical implementation, since the use of quantities related to the previous Newton iteration step, typically arising for mixed enhanced elements, is completely avoided. The stabilization matrices for the mixed method of incompatible modes and the mixed method of enhanced strains are obtained with volume and area bubble functions. Cook's membrane problem illustrates successfully the stabilization effect for bi-linear tetrahedral elements. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

New mixed finite elements for plane elasticity and Stokes equations

We consider mixed finite elements for the plane elasticity system and the Stokes equation. For the unmodified Hellinger-Reissner formulation of elasticity in which the stress and displacement fields are the primary unknowns, we derive two new nonconforming mixed finite elements of triangle type. Both elements use piecewise rigid motions to approximate the displacement and piecewise polynomial functions to approximate the stress, where no vertex degrees of freedom are involved. The two stress finite element ...     

平面线弹性纯位移边值问题低阶非协调矩形元的Robust多重网格方法

1引 言 对于各向同性,均匀介质的平面线弹性问题,当Lamé常数λ→∞(泊松率v→0.5)时,即对于几乎不可压介质,通常的协调有限元格式的解往往不再收敛到原问题的解,或者达不到最优收敛阶,这就是所谓的闭锁现象(见[3],[7],[8]及[10]).究其原因,在通常的有限元分析中,其误差估计的系数与λ有关,当λ→∞时,该系数将趋于无穷大.因此为克服闭锁现象就需要构造特殊的有限元格式,使得当λ→∞时,有限元逼近解仍然收敛到原问题的解.     

Multiscale methods for elliptic homogenization problems

In this article we study two families of multiscale methods for numerically solving elliptic homogenization problems. The recently developed multiscale finite element method [Hou and Wu, J Comp Phys 134 (1997), 169–189] captures the effect of microscales on macroscales through modification of finite element basis functions. Here we reformulate this method that captures the same effect through modification of bilinear forms in the finite element formulation. This new formulation is a general approach that can handle a large variety of differential problems and numerical methods. It can be easily extended to nonlinear problems and mixed finite element methods, for example. The latter extension is carried out in this article. The recently introduced heterogeneous multiscale method [Engquist and Engquist, Comm Math Sci 1 (2003), 87–132] is designed for efficient numerical solution of problems with multiscales and multiphysics. In the second part of this article, we study this method in mixed form (we call it the mixed heterogeneous multiscale method). We present a detailed analysis for stability and convergence of this new method. Estimates are obtained for the error between the homogenized and numerical multiscale solutions. Strategies for retrieving the microstructural information from the numerical solution are provided and analyzed. Relationship between the multiscale finite element and heterogeneous multiscale methods is discussed. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Adaptive FEM with Stabilized Elements for Parameter Identification of Incompressible Hyperelastic Materials

This work is concerned with the identification of material parameters for isotropic, incompressible hyperelastic material models. Besides the principal stretch-based strain-energy function by Ogden an invariant-based strain-energy function by Rivlin/Saunders is considered for which parameter sensitivities are derived. The identification is formulated as a least-squares minimization problem based on the finite element method to account for inhomogeneous states of stresses and strains in the experimental data which is obtained by optical measurements. For the finite element method low-order tetrahedral elements in a mixed displacement-pressure formulation with stabilization are considered. Special attention is payed to an adaptive mesh-refinement based on a goal-oriented a posteriori error indicator to gain reliable material parameters. To approximate error terms an element-wise recovery technique based on enhanced gradients is introduced. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)