In this paper, coupled electromechanical behavior of a vibrational energy harvesting system composed of a unimorph piezoelectric laminated beam with a large attached tip mass is investigated. To achieve this goal, first the electromechanically coupled partial differential equations governing the lateral displacement and output voltage of the harvester are extracted through exploiting the Hamilton’s principle. Considering vibration damping due to mechanical to electrical energy conversion, a complex modal analysis is performed to extract the complex eigenfrequencies and eigenfunctions of the system. Furthermore, an exact analytical solution is presented for the system response to the harmonic base excitations, including output voltage and harvested power. To validate the analytical results, at the next step a finite element simulation is conducted through ABAQUS software. To perform a fully-coupled analysis which brings into account the effect of harvesting circuit, user subroutine User-defined Amplitude (UAMP) is utilized to calculate the voltage–current relation and impose the correct electrical charge on the electrodes in each step by monitoring the output voltage of the system at previous time increments. Results of both analytical and numerical simulations are compared for a Micro-Electro-Mechanical Systems (MEMS) harvester as a case study, where a very good agreement is observed between them. © 2018, IMechE 2018.