Effective Computational Methods for the Modeling of Ferroelectroelastic Hysteresis Behavior