This article presents a numerical analysis to investigate the natural frequencies and harmonic response of a perforated cantilever beam attached to two layers of piezoelectric materials by using the finite element method for the first time. The bimorph piezoelectric is composed of 3 layers; two of them at the outer are piezoelectric, and the inner isotropic material. A higher order 3-D 20-node solid element that exhibits quadratic displacement behavior is exploited to discretize the isotropic layer, and coupled piezoelectric 3D element with twenty nodes is used to mesh the top and bottom layers. CIRCU94 element is added to act as a resistor part of the model. The proposed model is validated with previous works. The numerical parametric studies are presented to illustrate the effects of perforation geometry, the number of rows, the resistance on the natural frequencies, frequency response, and power. It is found that the thickness has a positive relationship with the natural frequency. Perforations help in producing higher voltage, and the best shape is rectangular perforations, and to produce higher voltage, two rows of rectangular perforations should be applied.