Abstract: A multivariable wavelet finite element method (FEM) is presented, which is based on Hellingger-Reissner variational principle. The interpolating wavelet bases are constructed in order to deal with boundary conditions conveniently, and two-dimensional interpolating wavelet bases in product form are used to construct the generalized field functions of thick plate. In calculating variables the stress-strain relations and the differential calculation are bypassed, resulting in high variable accuracy.