Piezoelectric bimorph actuators are used in a variety of applications, including micro positioning, vibration control, and micro robotics. The nature of the aforementioned applications calls for the dynamic characteristics identification of actuator at the embodiment design stage. For decades, many linear models have been presented to describe the dynamic behavior of this type of actuators; however, in many situations, such as resonant actuation, the piezoelectric actuators exhibit a softening nonlinear behavior; hence, an accurate dynamic model is demanded to properly predict the nonlinearity. In this study, first, the nonlinear stress-strain relationship of a piezoelectric material at high frequencies is modified. Then, based on the obtained constitutive equations and Euler-Bernoulli beam theory, a continuous nonlinear dynamic model for a piezoelectric bending actuator is presented. Next, the method of multiple scales is used to solve the discretized nonlinear differential equations. Finally, the results are compared with the ones obtained experimentally and nonlinear parameters are identified considering frequency response and phase response simultaneously. Also, in order to evaluate the accuracy of the proposed model, it is tested out of the identification range as well.