Abstract: A nonlinear constitutive model for piezoelectric ceramics is proposed, in which the polarization switching and electric saturation are taken into account. Based on the model, the nonlinear finite element analysis is implemented. Using localJ-integral as a fracture criterion, a relation curve of fracture loads against electric fields is obtained. Qualitatively, the curve is in agreement with the experimental observations reported by Park and Sun.