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)     

退化椭圆问题的最小二乘混合有限元方法

本文研究了电磁场中关于共振现象的一类退化的椭圆问题 ,提出了最小二乘混合有限元方法 .这一方法的好处是可以去掉传统混合元空间的LBB条件所得到的系数矩阵是对称正定的 ,使得法语解更加方便 .本文得到了最小二乘混合有限元方法的L2 和H1估计 .     

Numerical electroseismic modeling: A finite element approach

Electroseismics is a procedure that uses the conversion of electromagnetic to seismic waves in a fluid-saturated porous rock due to the electrokinetic phenomenon. This work presents a collection of continuous and discrete time finite element procedures for electroseismic modeling in poroelastic fluid-saturated media. The model involves the simultaneous solution of Biot’s equations of motion and Maxwell’s equations in a bounded domain, coupled via an electrokinetic coefficient, with appropriate initial conditions and employing absorbing boundary conditions at the artificial boundaries. The 3D case is formulated and analyzed in detail including results on the existence and uniqueness of the solution of the initial boundary value problem. Apriori error estimates for a continuous-time finite element procedure based on parallelepiped elements are derived, with Maxwell’s equations discretized in space using the lowest order mixed finite element spaces of Nédélec, while for Biot’s equations a nonconforming element for each component of the solid displacement vector and the vector part of the Raviart-Thomas-Nédélec of zero order for the fluid displacement vector are employed. A fully implicit discrete-time finite element method is also defined and its stability is demonstrated. The results are also extended to the case of tetrahedral elements. The 2D cases of compressional and vertically polarized shear waves coupled with the transverse magnetic polarization (PSVTM-mode) and horizontally polarized shear waves coupled with the transverse electric polarization (SHTE-mode) are also formulated and the corresponding finite element spaces are defined. The 1D SHTE initial boundary value problem is also formulated and approximately solved using a discrete-time finite element procedure, which was implemented to obtain the numerical examples presented.     

On the use of physical spline finite element method for acoustic scattering

This paper presents the comparison of physical spline finite element method (PSFEM), in which differential equations are incorporated into interpolations of basic elements, with least-squares finite element method (LSFEM) and mixed Galerkin finite element method (MGFEM) on the numerical solution of one dimensional Helmholtz equation applied to an acoustic scattering problem. Firstly, all three methods are explained in detail and then it is shown that PSFEM reaches higher precision in a shorter time with fewer nodes than the other methods. It is also observed that this method is well suited for high frequency acoustic problems. Consequently, the results of PSFEM point out better efficiency in terms of number of unknowns and accuracy level.     

Radiation and Scattering by Complex Conformal Antennas on a Circular Cylinder

A technique for characterizing and designing complex conformal antennas flush-mounted on a singly-curved surface is presented. This approach is based on the hybrid finite element–boundary integral (FE–BI) method. A related method was proposed in the past utilizing cylindrical-shell finite element and roof-top rectangular basis functions for the boundary integral. Although that method proved very powerful for analyzing cylindrical–rectangular patch arrays flush-mounted to a circular cylinder, the requirement for uniform meshing in the aperture ultimately limited its usefulness. In this present formulation, tetrahedral elements are used to expand the volumetric electric fields while similar basis functions are used for the boundary integral. The curvature of the aperture is explicitly included via the use of the circular cylinder dyadic Green's function. After presentation of the formulation and validation using several well-understood examples, an example is presented that illustrates the capabilities of this method for modeling complex conformal antennas heretofore not examined by rigorous methods due to inherent limitations of the various published methods.     

伪双曲方程的新混合有限元方法

构造分析一类二阶伪双曲方程的H1-Galerkin扩展混合有限元方法,该方法采用了扩展混合有限元方法和H1-Galerkin混合有限元方法相结合的技巧.新的格式同时保持了扩展混合有限元方法和H1-Galerkin混合有限元方法的优点.该混合格式与标准的混合格式相比能同时逼近三个变量:未知函数、梯度和流量(系数乘以梯度),并且不必满足LBB相容性条件.     

关于Ciarlet-Raviart混合有限元法的一个注记

本文推广解双调合方程的Ciarlet-Raviart混合有限元方案：用二次元逼近流函数φ．一次元逼近涡度-Δφ．在拟一致三角形剖分的条件下,证明了推广方案具有φ和-Δφ都用二次元逼近的标准Ciarlet－Raviart方案同样的精度阶．     

双调和方程混合元的一种新格式

本文介绍了双调和方程混合元的一种新格式,用双二次多项式逼近流函数,双一次多项式逼近涡函数.在拟一致矩形剖分的条件下,证明了此格式具有与C-R格式中分别用双二次多项式逼近相同的收敛阶.     

Normal and tangential continuous elements for least‐squares mixed finite element methods

We consider a finite element discretization of the primal first‐order least‐squares mixed formulation of the second‐order elliptic problem. The unknown variables are displacement and flux, which are approximated by equal‐order elements of the usual continuous element and the normal continuous element, respectively. We show that the error bounds for all variables are optimal. In addition, a field‐based least‐squares finite element method is proposed for the 3D‐magnetostatic problem, where both magnetic field and magnetic flux are taken as two independent variables which are approximated by the tangential continuous and the normal continuous elements, respectively. Coerciveness and optimal error bounds are obtained. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2004.     

三维Poisson方程特征值的四种有限元解及比较

应用三维EQ1rot元、三维Crouzeix-Raviart元、八节点等参数元、四面体线性元计算三维Poisson方程的近似特征值.计算结果表明:三维EQ1rot元和三维Crouzeix-Raviart元特征值下逼近准确特征值,八节点等参数元、四面体线性元特征值上逼近准确特征值,三维EQr1ot元和三维Crouzeix-Raviart元外推特征值下逼近准确特征值.计算结果还表明三维Crouzeix-Raviart元是一种计算效率较高的非协调元.     

粘弹性方程一种新的分裂正定混合元法

本文用分裂正定混合有限元方法研究二阶粘弹性方程. 首先构造一种新的分裂正定混合变分形式和基于这种分裂正定混合变分形式关于时间的半离散格式, 然后绕开关于空间变量的半离散化格式, 直接从时间半离散出发构造出全离散化的分裂正定混合有限元格式, 并给出这种分裂正定混合有限元解的误差估计. 这种研究思路使得理论论证变得更简单,这是处理二阶粘弹性方程的一种新的尝试.