A Weighted Least-Squares Mixed Finite Element Formulation for Quasi-Incompressible Elastodynamics

The purpose of this work is the application of the least-squares finite element method to an elastodynamic, quasi-incompressible problem under small strain assumptions. Therefore a mixed finite element based on a weighted least-squares formulation is developed. The L₂-norm minimization of the time-discretized residuals of the given first-order system of partial differential equations leads to a functional depending on displacements and stresses. In the numerical example the proposed mixed element is compared to an alternative approach, which is based on a least-squares mixed finite element with improved momentum balance, see [1]. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A Finite Element Formulation for the Nonlinear Analysis of Piezoelectric Three-Dimensional Beam Structures

A finite element formulation for a three-dimensional piezoelectric beam which includes geometrical and material nonlinearities is presented. To account for the piezoelectric effect, the coupling between the mechanical stress and the electrical displacement is considered. Based on the Timoshenko theory, an eccentric beam formulation is introduced which provides an efficient model to analyze piezoelectric structures. The geometrically nonlinear assumption allows the calculation of large deformations including buckling analysis. A quadratic approximation of the electric potential through the cross section of the beam ensures the fulfilment of the charge conservation law exactly. This assumption leads to a finite element formulation with six mechanical and five electrical degrees of freedom per node. To take into account the typical ferroelectric hysteresis phenomena, a nonlinear material model is essential. For this purpose, the phenomenological Preisach model is implemented into the beam formulation which provides an efficient determination of the remanent part of the polarization. The applicability of the introduced beam formulation is discussed with respect to available data from literature. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A New Finite Element Space for Expanded Mixed Finite Element Method

In this article, we propose a new finite element space Λ$_h$ for the expanded mixed finite element method (EMFEM) for second-order elliptic problems to guarantee its computing capability and reduce the computation cost. The new finite element space Λ$_h$ is designed in such a way that the strong requirement V$_h\subset$Λ$_h$ in [9] is weakened to {v$_h\in$V$_h$; divv$_h$=0}$\subset$Λ$_h$ so that it needs fewer degrees of freedom than its classical counterpart. Furthermore, the new Λ$_h$ coupled with the Raviart-Thomas space satisfies the inf-sup condition, which is crucial to the computation of mixed methods for its close relation to the behavior of the smallest nonzero eigenvalue of the stiff matrix, and thus the existence, uniqueness and optimal approximate capability of the EMFEM solution are proved for rectangular partitions in $\mathbb{R}^d, d=2,3$ and for triangular partitions in $\mathbb{R}^2$. Also, the solvability of the EMFEM for triangular partition in $\mathbb{R}^3$ can be directly proved without the inf-sup condition. Numerical experiments are conducted to confirm these theoretical findings.     

拟线性抛物型积分微分方程的一个新最低阶混合元格式

对一类拟线性抛物型积分微分方程构造了一个新的最低阶三角形协调混合元格式,并直接利用单元插值的性质,给出了相应的收敛性分析和H~1-模及L~2-模意义下的最优误差估计.     

A Least-Squares Mixed Finite Element for Quasi-Incompressible Elastodynamics

In the present work, a mixed finite element based on a modified least-squares formulation is proposed. Here, we consider the time-dependent equations for quasi-incompressible elastodynamics under small strain assumptions. The main goal is to obtain an accurate approximation of both displacements and stresses in particular for the lowest-order element. Basis for the element formulation is a weak form resulting from a least-squares method. The L₂-norm minimization of the time-discretized residuals of the given first-order system leads to a functional depending on approximations for displacements and stresses. By introducing a time-independent displacement test function, a weak form is derived. A numerical example concerning quasi-incompressible elasticity shows the performance of the approach for the lowest-order element RT₀P₁as. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

一个最低阶新混合元格式的超收敛分析

对一类拟线性抛物型积分微分方程构造一个新的最低阶三角形协调混合元格式,在抛弃传统有限元分析中不可缺少的工具Ritz-Vblterra投影的前提下,直接利用单元插值的性质及积分恒等式技巧,给出了相应变量的超逼近及超收敛结果,弥补了已有文献的不足.     

A Mixed Stiffness Finite Element Method for Navier-Stokes Equation

The paper is devoted to the study and analysis of the mixed stiffness finite element method for the Navier-Stokes equations, based on a formulation of velocity-pressure-stress deviatorics. The method used low order Lagrange elements and leads to optimal error order of convergence for velocity, pressure, and stress deviatorics by means of mesh-dependent norms defined in this paper. The main advantage of the MSFEM is that the stream function can not only be employed to satisfy the divergence constraint but stress deviatorics can also be eliminated at the element level so that it is unnecessary to solve a larger algebraic system containing stress multipliers, or to develop a special code for computing the MSFE solutions of Navier-Stokes equations because we can use the computing codes used in solving the Navier-Stokes equations with the velocity-pressure formulation, or even the computing codes used in solving the problems of solid mechanics.     

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

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

A Four-node Finite Element for Piezoelectric Shell Structures in Convective Co-ordinates

In this contribution, an isoparametric piezoelectric shell element is presented which is based on convective coordinates and which allows for the analysis of arbitrary shell geometries. A two-field variation formulation [1, 2] is used in which the displacements and the electric potentials serve as independent variables. Especially, for thin-walled structures under certain boundary conditions and load cases, the displacement based element tend to shear and membrane locking. In order to avoid this poor behaviour, the Assumed Natural Strain (ANS) method [3] is introduced into the piezoelectric shell element. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

蒸汽沉淀化学反应方程混合元法的离散格式研究

蒸汽沉淀化学反应过程有着极其广泛的应用,其数学模型归结为一个包含流速场,温度场,压力场和气体溶质场的非线性偏微分方程组．用混合有限元方法研究蒸汽沉淀化学反应方程组,导出其半离散化和全离散化的混合元格式,并证明这些格式的解的存在性和收敛性(误差估计)．用混合元法处理究蒸汽沉淀化学反应方程组,可以同时求出流速场,温度场,压力场和气体溶质场的数值解. 因此该研究既具有重要的理论意义,又具有广泛的应用前景．     

Raviart-Thomas混合元的超收敛

考虑二阶椭圆方程Dirichlet边值问题在正则矩形网格上k阶RaviartThomas混合有限元的超收敛.对有限元解经插值处理后,与通常的有限元最优误差估计相比,收敛速度提高了两阶.     

A Mixed Finite Element Method for Fourth Order Eigenvalue Problems

A new mixed finite element method for approximating eigenpairs of IV order elliptic eigenvalue problems with Dirichlet boundary conditions has been given. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

A Piezoelectric Finite Shell Element for the Geometrically Nonlinear Analysis of Sensors

A geometrically nonlinear finite element formulation to analyze piezoelectric shell structures is presented. The formulation is based on the mixed field variational functional of Hu–Washizu. Within this variational principle the independent fields are displacements, electric potential, strains, electric field, stresses and dielectric displacements. The mixed formulation allows an interpolation of the strains and the electric field through the shell thickness, which is an essential advantage when using a three dimensional material law. It is remarked that no simplification regarding the constitutive relation is assumed. The normal zero stress condition and the normal zero dielectric displacement condition are enforced by the independent resultant stress and resultant dielectric displacement fields. The shell structure is modeled by a reference surface with a four node element. Each node possesses six mechanical degrees of freedom, three displacements and three rotations, and one electrical degree of freedom, which is the difference of the electric potential through the shell thickness. The developed mixed hybrid shell element fulfills the in–plane, bending and shear patch tests, which have been adopted for coupled field problems. A numerical investigation of a smart antenna demonstrates the applicability of the piezoelectric shell element under the consideration of geometrical nonlinearity. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

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

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

广义神经传播方程的一种修正混合有限元方法的误差分析

利用修正的H~1-Galerkin混合有限元方法研究了广义神经传播方程,论证了其半离散解的存在唯一性,得到了半离散解的最优阶误差估计,该方法的优点是不需验证LBB相容性条件.     

周期结构区域中压电问题的双尺度有限元方法

近年来纤维压电复合材料的力电性能预测已发展为一个重要的研究领域．对力电耦合周期结构的复合材料问题,通过引入匹配的边界层得到了电势与位移解的新型双尺度有限元计算方法,建立了电势与位移的双尺度耦合关系,分析了双尺度有限元解的误差.数值算例验证了方法的有效性．     

A Characteristic Mixed Finite Element Two-Grid Method for Compressible Miscible Displacement Problem

A nonlinear parabolic system is derived to describe compressible miscible displacement in a porous medium. The concentration equation is treated by a mixed finite element method with characteristics (CMFEM) and the pressure equation is treated by a parabolic mixed finite element method (PMFEM). Two-grid algorithm is considered to linearize nonlinear coupled system of two parabolic partial differential equations. Moreover, the $L^q$ error estimates are conducted for the pressure, Darcy velocity and concentration variables in the two-grid solutions. Both theoretical analysis and numerical experiments are presented to show that the two-grid algorithm is very effective.     

非线性抛物问题的二重网格混合元算法

针对一类非线性抛物方程的混合元形式,本文提出了二重网格算法.该算法是在网格大小为H的粗网格上求解—个非线性系统,再在网格大小为h的细网格上进行两次线性计算.算法第二步和第三步的误差分别为O(△_t~2 h~(k 1) H~(2K 2)),O(△_t~2 h~(k 1) h~(-d/2)H~(4k 4)),其中k为逼近空间的多项式的次数,d为空间维数.该估计对H的选取起了很大的作用.对于粗网格上的非线性计算,本文给出了L~p(2≤p<∞)模误差估计.     

Plate Contact问题的混合有限元逼近

论文考虑了Plate Contact问题的混合有限元逼近,其变分问题为第二类四阶椭圆变分不等问题.首先,根据正则化方法,得到原问题的正则化问题.再根据网格依赖范数技巧,考虑了正则化问题的Ciarlet-Raviart混合有限元逼近,并证明了真解与逼近解之间的误差估计.最后通过数值算例验证了理论分析的结果.