Estimates for the combined effect of boundaryapproximation and numerical integration on the approximation of(simple) eigenvalues and eigenvectors of 4th order eigenvalue problems with variable/constant coefficientsin convex domains with curved boundary by an isoparametric mixed finite element method, which,in the particular case of bending problems ofaniso-/ortho-/isotropic plates with variable/constant thickness, gives a simultaneous approximation to bending momenttensor field $\Psi= (\psi_{ij})_{1 \le i,j \le 2}$ anddisplacement field `u', have been developed.