Summary: There have been many efforts, dating back four decades, to develop stable mixed finite elements for the stress-displacement formulation of the plane elasticity system. This requires the development of a compatible pair of finite element spaces, one to discretize the space of symmetric tensors in which the stress field is sought, and one to discretize the space of vector fields in which the displacement is sought. Although there are number of well-known mixed finite element pairs known for the analogous problem involving vector fields and scalar fields, the symmetry of the stress field is a substantial additional difficulty, and the elements presented here are the first ones using polynomial shape functions which are known to be stable. We present a family of such pairs of finite element spaces, one for each polynomial degree, beginning with degree two for the stress and degree one for the displacement, and show stability and optimal order approximation. We also analyze some obstructions to the construction of such finite element spaces, which account for the paucity of elements available.

##### MSC:

74S05 | Finite element methods in solid mechanics |

74G15 | Numerical approximation of solutions for equilibrium problems in solid mechanics |

65N30 | Finite elements, Rayleigh-Ritz and Galerkin methods, finite methods (BVP of PDE) |

65N12 | Stability and convergence of numerical methods (BVP of PDE) |

65N15 | Error bounds (BVP of PDE) |