A New Family of Efficient Conforming Mixed Finite Elements on Both Rectangular and Cuboid Meshes for Linear Elasticity in the Symmetric Formulation

A new family of mixed finite elements is proposed for solving the classical Hellinger--Reissner mixed problem of the elasticity equations. For two dimensions, the normal stress of the matrix-valued stress field is approximated by an enriched Brezzi--Douglas--Fortin--Marini element of order $k$, the shear stress by the serendipity element of order $k$, and the displacement field by an enriched discontinuous vector-valued $P_{k-1}$ element. The degrees of freedom on each element of the lowest order element, which is of first order, are 10 plus 4. For three dimensions, the normal stress is approximated by an enriched Raviart--Thomas element of order $k$, each component of the shear stress by a product space of the serendipity element space of two variables and the space of polynomials of degree $\leq k-1$ with respect to the rest variable, and the displacement field by an enriched discontinuous vector-valued $Q_{k-1}$ element. The degrees of freedom on each element of the lowest order element, which is of fi...