Stabilized dual-mixed method for the problem of linear elasticity with mixed boundary conditions

We extend the applicability of the augmented dual-mixed method introduced recently in [4, 5] to the problem of linear elasticity with mixed boundary conditions. The method is based on the Hellinger–Reissner principle and the symmetry of the stress tensor is imposed in a weak sense. The Neuman boundary condition is prescribed in the finite element space. Then, suitable Galerkin least-squares type terms are added in order to obtain an augmented variational formulation which is coercive in the whole space. This allows to use any finite element subspaces to approximate the displacement, the Cauchy stress tensor and the rotation.