A finite element model for rate-dependent behavior of ferroelectric ceramics

In this article, materials within a crystallite are modeled by continuum particles consisting of various types of ferroelectric variants which are characterized by their mass fractions. The constitutive behavior of each type of variant is characterized by a proposed Helmholtz free energy potential. Polarization switching is modeled by continuous changes of mass fractions which are governed by a onset criterion and a kinetic relation. A finite element algorithm is developed using the virtual work principle. The simulated results on the rate dependence in the polarization and strain responses to applied alternating electric field of different frequencies are in qualitative consistence with experimental observations. The rate-dependent behavior is explained in terms of changes of mass fractions of the variants that polarization switching involves, in response to the loading programs of different loading rates.