Ferroelectric materials have been wide deployed in many applications due to their capability of converting energy between the electrical and mechanical types automatically.But the most challenging task is nonlinear (hysteretic) behaviors in both the mechanical and electrical fields caused by polarization.In order to model and control the ferroelectric materials and devices, it is essential to construct a suitable model for its hysteretic dynamics.On the other side.it is always beneficial if the hysteretic dynamics can be described accurately by Ordinary Differential Equations (ODE) since theories for system analysis and control were well established for ODE.