This paper is focused on the numerical simulation of piezoelectric effects, in view of their application in MEMS harvesters and actuators. In order to obtain a fast and accurate prediction of dynamic behavior, a simple 1-DOF model has been introduced for thin piezoelectric beams. The validation of the model has been carried out by comparison with full 3D simulations via commercial codes