  We consider a mixed finite element method for approximating the solution of
nearly incompressible elasticity and Stokes equations. The finite element
method is based on quadrilateral and hexahedral triangulation using primal and
dual meshes. We use the standard bilinear and trilinear finite element space
enriched with element-wise defined bubble functions with respect to the primal
mesh for the displacement or velocity, whereas the pressure space is
discretised by using a piecewise constant finite element space with respect to
the dual mesh.
