This paper deals with the approximation of the buckling coefficients and modes of a clamped plate modeled by the Reissner-Mindlin equations. These coefficients are related with the eigenvalues of a non-compact operator. We give a spectral characterization of this operator and show that the relevant buckling coefficients correspond to isolated nondefective eigenvalues. Then we consider the numerical solution of the buckling problem. For the finite element approximation of Reissner-Mindlin equations, it is well known that some kind of reduced integration or mixed interpolation has to be used to avoid locking. In particular we consider Duran-Liberman elements, which have been already proved to be locking-free for load and vibration problems. We adapt the classical approximation theory for non-compact operators to obtain optimal order error estimates for the eigenfunctions and a double order for the eigenvalues. These estimates are valid with constants independent of the plate thickness. W...