In this paper, a kind of piezoelectric ceramics model based on hysteretic nonlinear theory has been developed. Van de Pol nonlinear difference item was introduced to interpret the hysteresis phenomenon of the voltage-strain curve of piezoelectric ceramics. The coupling relationship between voltage and stress was obtained in partial least-square regression method to describe the driftage phenomenon of the voltage-strain curve in different stress. Based on above, the final relationship among strain, stress and voltage was set up. The results of significance test showed that the new model could describe the hysteresis characteristics of piezoelectric ceramics in different stress well. The new piezoelectric ceramics model considers the effect of stress, and is easy to be analyzed in theory, which is helpful to vibration control.