The hysteresis phenomena of ferroelectric/ferroelastic material in polarization procedure are investigated. Some assumptions are presented based on the published experimental data. The electrical yielding criterion, mechanical yielding criterion and isotropic hardening model are established. The flow theory in incremental forms in polarization procedure is presented. The nonlinear constitutive law for electrical-mechanical coupling is proposed phenomenologically. Finally, the nonlinear constitutive law expressed in a form of matrices and vectors, which is immediately associated with finite element analysis, is formulated. In the example problem of a rectangular specimen subjected to a uniaxial electric field, the procedure from virgin state to fully polarized state is simulated. Afterward, a uniaxial compressive loading is applied to depolarizing the specimen. Results are in agreement with the experimental data.ferroelectric/ferroelastic material, constitutive law, nonlinear, finite element method.Piezoelectric materials, due to their coupling behavior between electric and mechanical fields, are extensively used for electromechanical and electronic devices, such as actuators, sensors and transducers in intelligent systems. If such a material is polarized by an electric field above the coercive field at a temperature below the Curie point, its response to a small signal may be characterized as linear. Nowadays, application of piezoelectric materials involves rigorous loading and complicated geometry of the components so that the assumption of small signals is no longer justified. Instead, the nonlinear behavior of the materials may become dominant. Hence, it is significant to study the numerical method for simulating the nonlinear constitutive behavior of the electric devices made of piezoelectric materials.Curves characterizing piezoelectricity and ferroelectricity may be found in many textbooks and articles. These hystereses demonstrate ferroelastic and ferroelectric nonlinearity and the underlying physical mechanisms. Attempts for analysis and description of the ferroelectricity and