No abstract

Abstract. We present a new stabilized mixed finite element method for the linear elasticity problem in R 2 . The approach is based on the introduction of Galerkin least-squares terms arising from the constitutive and equilibrium equations, and from the relation defining the rotation in terms of the displacement. We show that the resulting augmented variational formulation and the associated Galerkin scheme are well posed, and that the latter becomes locking-free and asymptotically locking-free for Dirichlet and mixed boundary conditions, respectively. In particular, the discrete scheme allows the utilization of Raviart-Thomas spaces of lowest order for the stress tensor, piecewise linear elements for the displacement, and piecewise constants for the rotation. In the case of mixed boundary conditions, the essential one (Neumann) is imposed weakly, which yields the introduction of the trace of the displacement as a suitable Lagrange multiplier. This trace is then approximated by piecewise linear elements on an independent partition of the Neumann boundary whose mesh size needs to satisfy a compatibility condition with the mesh size associated to the triangulation of the domain. Several numerical results illustrating the good performance of the augmented mixed finite element scheme in the case of Dirichlet boundary conditions are also reported.Mathematics Subject Classification. 65N12, 65N15, 65N30, 74B05.

In this note we analyze a modified mixed finite element method for second-order elliptic equations in divergence form. As a model we consider the Poisson problem with mixed boundary conditions in a polygonal domain of R 2 . The Neumann (essential) condition is imposed here in a weak sense, which yields the introduction of a Lagrange multiplier given by the trace of the solution on the corresponding boundary. This approach allows to handle nonhomogeneous Neumann boundary conditions, theoretically and computationally, in an alternative and usually easier way. Then we utilize the classical Babuška-Brezzi theory to show that the resulting mixed variational formulation is well posed. In addition, we use Raviart-Thomas spaces to define the associated finite element method and, applying some elliptic regularity results, we prove the stability, unique solvability, and convergence of this discrete scheme, under appropriate assumptions on the mesh sizes. Finally, we provide numerical results illustrating the performance of the algorithm for smooth and singular problems.

No abstract

No abstract