A comprehensive theoretical analysis of a dynamic thermoferroelectric pre-stressed bimorph energy harvester is performed. The analysis takes into account pyroelectric and thermal expansion effects and general analytical expressions for the energy conversation coefficients are presented for a bilayer. These derived coefficients (transformation coefficients) are for situation when mechanical, electrical, and thermal fields are present and they are important for sensors, actuators, and energy harvesters. They are a function of material properties, location of boundary conditions, vibration frequency, and in plane compressive/tensile follower force. Numerical simulations of the analytical results are presented. Effects of volume fraction, material properties, applied mechanical loads, and boundary conditions on the harvesting coefficients are shown in the figures. The results for a cantilever and a simply-supported beam are obtained as particular cases. The result for a low frequency (static) system is obtained as a particular case by approaching the vibration frequency to zero. It is shown that volume fraction, material properties, plain compressive/tensile follower force, the location of the boundary conditions, and the vibrational frequency of the bimorph strongly influence the strain distribution, and this in turn influences the charge coefficient and the generation of energy.