Abstract
We conduct numerical experiments of a recently proposed elasto-plastic constitutive model for nonlinear piezoelectric problems, which is based on thermo mechanics. This model is in the same framework of classical elasto-plasticity for metals and hence, computational return mapping algorithm can be constructed for nonlinear electro-mechanical finite element analyses. The numerical experiments under electrical loadings show not only its numerical robustness for large increments but also its capability for the high representation of phenomenological behaviors.