This work presents multiphysics numerical analysis of piezoelectric actuators realized using the finite element method (FEM) and their performances to analyze the structure-electric interaction in three-dimensional (3D) piezoelectric continua. Here, we choose the piezoelectric bimorph actuator without the metal shim and with the metal shim as low-frequency problems and a surface acoustic wave device as a high-frequency problem. More attention is given to low-frequency problems because in our application micro air vehicle’s wings are actuated by piezoelectric bimorph actuators at low frequency. We employed the Newmark’s time integration and the central difference time integration to study the dynamic response of piezoelectric actuators. Monolithic coupling, noniterative partitioned coupling and partitioned iterative coupling schemes are presented. In partitioned iterative coupling schemes, the block Jacobi and the block Gauss–Seidel methods are employed. Resonance characteristics are very important in micro-electro-mechanical system (MEMS) applications. Therefore, using our proposed coupled algorithms, the resonance characteristics of bimorph actuator is analyzed. Comparison of the accuracy and computational efficiency of the proposed numerical finite element coupled algorithms have been carried out for 3D structure–electric interaction problems of a piezoelectric actuator. The numerical results obtained by the proposed algorithms are in good agreement with the theoretical solutions.