Numerical analysis of a locking-free mixed finite element method for a bending moment formulation of Reissner-Mindlin plate model