The Arnold–Winther mixed FEM in linear elasticity. Part I: Implementation and numerical verification
نویسندگان
چکیده
Following some de Rham complex, Arnold and Winther have recently proposed a symmetric mixed finite element method (MFEM) in linear elasticity. This paper describes the implementation of the symmetric MFEM and its 30 30 local stress stiffness matrices and studies the implementation of the lowest-order scheme for general boundary conditions. Numerical experiments in model examples in computational mechanics illustrate the robust locking-free convergence behavior and support the theoretically predicted experimental convergence rates. 2008 Elsevier B.V. All rights reserved.
منابع مشابه
A Multigrid Preconditioner for the Mixed Formulation of Linear Plane Elasticity
In this paper, we develop a multigrid preconditioner for the discrete system of linear equations which results from the mixed formulation of the linear plane elasticity problem using the Arnold-Winther elements. This, in turn, can be reduced to the problem of finding a multigrid preconditioner for the form (·, ·)+(div ·,div ·) in the symmetric tensor space resulting from Arnold-Winther elements...
متن کاملRectangular mixed elements for elasticity with weakly imposed symmetry condition
We present new rectangular mixed finite elements for linear elasticity. The approach is based on a modification of the Hellinger-Reissner functional in which the symmetry of the stress field is enforced weakly through the introduction of a Lagrange multiplier. The elements are analogues of the lowest order elements described in Arnold, Falk and Winther [ Mixed finite element methods for linear ...
متن کاملLower Order Finite Element Approximations of Symmetric Tensors on Simplicial Grids in R
In this paper, we construct, in a unified fashion, lower order finite element subspaces of spaces of symmetric tensors with square-integrable divergence on a domain in any dimension. These subspaces are essentially the symmetric H(div) − Pk (1 ≤ k ≤ n) tensor spaces, enriched, for each n − 1 dimensional simplex, by (n+1)n 2 H(div) − Pn+1 bubble functions when 1 ≤ k ≤ n − 1, and by (n−1)n 2 H(di...
متن کاملModeling quasi-static crack growth with the extended finite element method Part II: Numerical applications
In Part I [Int. J. Solids Struct., 2003], we described the implementation of the extended finite element method (XFEM) within Dynaflowe, a standard finite element package. In our implementation, we focused on two-dimensional crack modeling in linear elasticity. For crack modeling in the X-FEM, a discontinuous function and the near-tip asymptotic functions are added to the finite element approxi...
متن کاملSemilinear Mixed Problems on Hilbert Complexes and Their Numerical Approximation
Arnold, Falk, and Winther recently showed [Bull. Amer. Math. Soc. 47 (2010), 281–354] that linear, mixed variational problems, and their numerical approximation by mixed finite element methods, can be studied using the powerful, abstract language of Hilbert complexes. In another recent article [arXiv:1005.4455], we extended the Arnold–Falk–Winther framework by analyzing variational crimes (a la...
متن کامل