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.