With the increasing world population and the limited availability of renewable energy sources, all countries face the challenge of energy production. The need to pay attention to renewable energy sources is useful in removable piezoaeroelastic harvesters in order to prevent the loss of this type of energy. In this type of harvester, energy is harvested through interaction between the structure and the fluid, and then transferred into electrical energy at the scale of milliwatt or microwatt. In the present study, it is tried to increase the efficiency of energy harvesting models. Cantilevered beam is considered, in which there are two piezoelectric layers on the top and bottom of the beam, which are underlying its base excitation, and at the free end of the beam, the airfoil is hinged and exposed to the air flow. A similar model is also modeled, with the difference that the circular cylinder is attached at the free end of the beam. The purpose is to investigate the start speed of the flutter and the nonlinear response of the systems in both linear and nonlinear geometric strains afterwards. In order to investigate the speed of the start of the flutter, for the models (airfoil and circular cylinder), the Peters linear model was used for simplicity and better application, and the linearized model of the Van der Pol equation was used around the origin coordinates, respectively. To investigate the nonlinear response of the airfoil fluctuation, an ONERA model describing dynamic stall is used and for the circular cylinder, the use of the Van der Pol model is used. It is also assumed that the base excitation is harmonic. At the first, the nonlinear differential equation governing the transverse vibration of the beam was converted to a number of ordinary differential equations using the Galerkin method and then dissolved by a Rang-Kuta numerical method. Finally, by studying the response of the governing equations, the influence of system parameters on the vibrational behavior of the beam and the output voltage is studied. The results show that applying the base excitation increases the amplitude of the vibration and output voltage for both models. In addition, the displacement of the ends of the beam, the flap rotational airfoil and the output voltage for both models in linear and nonlinear states are investigated and compared. The results show that, in the nonlinear state, the vibrational amplitude and the output voltage are increased compared to the linear state. Keywords : Energy harvesting, Airfoil, Aeroelasticity, Flutter, ONERA model, Peters model, Dynamic stall, Van der Pol