Abstract: In this paper, based on the multi-scale finite element method, the meso-heterogeneous information is introduced into the macroscopic response by constructing the displacement multi-scale basis functions and the potential multi-scale basis functions. Coupling additional terms are also incorporated for the coupling between the displacement fields in each coordinate direction. The constructed multi-scale basis functions can directly and efficiently transfer the information between coarse scale and fine-scale, addressing problems on the macro scale and reducing the computational costs. The numerical simulation results show the multi-scale finite element method not only guarantees calculation accuracy but also achieves higher computational efficiency compared with the traditional finite element method, and thus it provides an effective means for the numerical simulation of the dynamic problems of heterogeneous piezoelectric structures.