In this paper, the quadratic nonconforming brick element （MSLK element） intro- duced in [10] is used for the 3D Stokes equations. The instability for the mixed element pair MSLK-P1 is analyzed, where the vector-valued MSLK element approximates the velocity and the piecewise P1 element approximates the pressure. As a cure, we adopt the piecewise P1 macroelement to discretize the pressure instead of the standard piecewise P1 element on cuboid meshes. This new pair is stable and the optimal error estimate is achieved. Numerical examples verify our theoretical analysis.
