“…Another framework to construct stable weakly symmetric mixed finite elements was presented in [11], where two approaches were particularly proposed with the first one based on Stokes's problems and the second one based on interpolation operators. To keep the symmetry of discrete stresses, a second way is to relax the continuity of the normal components of discrete stresses across the internal edges or faces of grids, which leads to nonconforming mixed FEMs with strong symmetric stress tensor [4,9,10,25,32,33,37,[47][48][49]. In 2002, based on the elasticity complexes, the first family of symmetric conforming mixed elements with polynomial shape functions was proposed for two-dimensional cases in [8], which was extended to threedimensional cases in [2].…”