# Construction and Application of Experimental Formula for Nonlinear Behavior of Ferroelectric Ceramics Switched by Electric Field at Room Temperature during Temperature Rise

## Article information

## Abstract

A poled lead zirconate titanate (PZT) cube specimen that is switched by an electric field at room temperature is subject to temperature increase. Changes in polarization and thermal expansion coefficients are measured during temperature rise. The measured data are analyzed to obtain changes in pyroelectric coefficient and strain during temperature change. Empirical formulae are developed using linear or quadratic curve fitting to the data. The nonlinear behavior of the materials during temperature increase is predicted using the developed formulae. It is shown that the calculation results can be compared successfully with the measured values, which proves the accuracy and reliability of the developed formulae for the nonlinear behavior of the materials during temperature changes.

**Keywords:**Thermal expansion coefficient; Temperature; Pyroelectric coefficient; Switching; Empirical formula

## 1. Introduction

Piezoelectric ceramics are widely used in devices and equipment such as ferroelectric memory (FRAM), ultrasonic motors, sonar detectors, and micro electro mechanical systems (MEMS). Research on the application and development of piezoelectric ceramics is also widely pursued. The application of a significant load or electric field to a piezoelectric system can result in excessive concentration of stress or electric fields in the ceramic material due to the complexity of the system or the presence of defects, which causes unexpected domain switching and rapid temperature changes. As a result, the internal structure of the ceramic material often changes to produce macroscopic nonlinear behavior and variations in the material properties. Degradation in system performance is observed due to such unexpected domain switching. In order to prevent system instability due to nonlinear behavior, experimental observation and analysis of nonlinear behavior of piezoelectric ceramic materials under various load conditions must be made. A wide range of related studies have been carried out. Zhou and Kamlah,1) Liu and Huber,2) and Kim3) observed and analyzed the nonlinear behavior of the materials when electric field and stress are applied at room temperature. Recently, research on the nonlinear behavior of the materials at room temperature has been extended to high temperatures. Lee and Kim4) and Selten *et al*.5) measured and analyzed the nonlinear behavior of the materials caused by electric fields at room and high temperatures, while Weber *et al.*6) measured and analyzed the nonlinear behavior caused by stresses at room and high temperatures. Grunbichler *et al.*7) used the finite element method to predict and understand the nonlinear behavior of piezoelectric actuators due to stress, electric fields, and temperature changes. Also, to develop an empirical formula for predicting experimental observations of piezoelectric materials, Ji and Kim8,9,10) observed and analyzed the nonlinear behavior of the materials when they are subject to stress and electric fields at room and high temperatures.

In the present work, relatively simple empirical formulae for the experiment results of Ji and Kim11) are presented. Moreover, the presented empirical formulae were used to calculate the polarization and strain responses of piezoelectric materials during temperature increases, and comparisons were made with measurements.

## 2. Experimental Procedure

In the present work, the experimental data of Ji and Kim11) were used to present new formulae for the prediction of ferroelectric ceramic behavior at high temperatures. The experimental method of Ji and Kim11) is presented simply here. For ferroelectric specimens, PZT specimens of regular hexahedron shape with 10 mm sides were used. The followings are the properties of the specimen material, provided by the manufacturer (PZT5H1, Morgan Technical Ceramics, UK): the Curie point is 200°C, coupling factor *k** _{p}* = 0.60, piezoelectric coefficients

*d*

_{31}= −250 × 10

^{−12}mV

^{−1}and

*d*

_{33}= 620 × 10

^{−12}mV

^{−1}, elastic compliance coefficients

*s*

_{33}= 21.9 × 10

^{−12}m

^{2}N

^{−1}and

*s*

_{11}= 17.7 × 10

^{−12}m

^{2}N

^{−1}. Here, the subscript 3 refers to the polarization direction; the subscript 1 refers to the transverse polarization direction. Initially, the specimen is poled in the negative direction and a total of 14 electric fields of different magnitude were applied at the reference temperature of 20°C. The range of reference remnant polarization

^{−2}. After poling, the leakage current density was measured over 5,000 s and the specimen was stabilized, then followed by temperature increases to 110°C at a rate of 1.2°Cmin

^{−1}along with the invar specimen (Product No. 318-0285-3, Danyang, China). During temperature increases, remnant polarization

*ɛ*

_{G/R}and

*ɛ*

_{G/S}of the invar and ceramic specimens were measured. Here, the thermal output changes refer to the variations of strain gage outputs value during temperature rises. Then, thermal expansion coefficients in the longitudinal and transverse directions

*α*

_{3}and

*α*

_{1}of the ceramic specimen were calculated using the thermal expansion coefficient of the invar specimen

*α*

_{R}.12) Polarization was measured indirectly using the Sawyer-Tower circuit; a Keithley 6514 was used to measure the voltage of the capacitor serially connected to the specimen. Temperature was controlled by placing the specimen in an oil bath (2408 PID controller, EUROTHERM, UK), shown in Fig. 1, filled with an insulating oil (MICTRANS Class1-No2, MICHANG OIL IND. CO., Pusan, Korea), and using a heat coil installed on the bottom of the oil bath. All data were collected at a sampling rate of 100 Hz and the data were processed using the LABVIEW program through a DAQ board (PCI 6221, National Instruments, TX, USA).

## 3. Results and Discussion

### 3.1. Pyroelectric Coefficient and Strain Calculations During Temperature Increases

Electric fields of various magnitudes were applied in the direction opposite to the polarization of PZT specimens, causing poling in the negative longitudinal direction at reference temperature 20°C. Switching is induced by the applied electric field and an internal state of specific reference remnant polarization is obtained. Then, the electric field was removed, and temperature is increased to 110°C. During temperature increases, remnant polarization
*α*_{3} and *α*_{1} were measured. Fig. 2(a) shows measured remnant polarization
*α*_{3} and *α*_{1} during temperature increases. Only five states of reference remnant polarization
*p*_{3} is obtained. The values of pyroelectric coefficients, calculated for the temperature range between 20°C to 110°C, for five chosen reference remnant polarization states, are shown as symbols in Fig. 3(a). In the figure, pyroelectric coefficients are plotted with respect to temperature for five constant values of reference remnant polarization
*p*_{3} fit well with straight line, given by Equation (1) below:

where, *a** _{P}* refers to the slope of the fitting straight line and

*b*

*refers to the intercept of the line at zero temperature. The values of*

_{P}*a*

*and*

_{P}*b*

*are plotted verse*

_{P}^{−2}and +0.263 Cm

^{−2}in Fig. 3(b). In the figure,

*a*

*is shown to fit with a straight line and*

_{P}*b*

*with a parabolic curve, both expressed as*

_{P}Substitution of Eq. (2) into Eq. (1) results in Eq. (3).

Using Eq. (3), one can estimate pyroelectric coefficient *p*_{3} when reference remnant polarization
*θ* are given.

Now turn to the calculation of longitudinal and transverse strains. Thermal expansion coefficient is the rate of change in strain with respect to temperature. Thus, Eq. (4) below can be used to calculate strains from measured thermal expansion coefficients.

Among fourteen measured states, longitudinal and transverse strains are calculated for five specific states of remnant strains and plotted with temperature in Figs. 4(a) and (b), respectively. In Fig. 4(a), as the magnitude of electric field increases,
^{−6}, it increases with temperature. This can be explained by microscopic process of domain switching. When the magnitude of electric field is small, switching is not sufficient and most domains remain parallel to the longitudinal *x*_{3} direction. When electric field is sufficiently large, the ratio of domains with polarizations parallel to the transverse *x*_{1} axis increases due to electric field-induced switching, leading to increases in
*x*_{1} axis, thus resulting in lower increases in

Let us apply the equations in Eqs. (1) – (3) to calculate strains. The same procedure is used for both longitudinal and transverse strains, so only the calculation of longitudinal strain will be dealt with here. The variations of longitudinal thermal expansion coefficient *α*_{3} with temperature for five chosen states of

where *a*_{S3} refers to the slope of fitting straight line and *b*_{S3} the intercept of the line at zero the temperature. The distributions of *a*_{S3} and *b*_{S3} over
^{−6} to −458.2 × 10^{−6} are plotted in Fig. 5(b). The data of *a*_{S3} and *b*_{S3}, marked with symbols, fit well with straight lines, expressed as

Substitution of Eq. (6) into Eq. (5) results in Eq. (7) below.

Using Eq. (7), one can estimate longitudinal thermal expansion coefficient *α*_{3} when reference remnant longitudinal strain
*θ* are given. In a similar manner, one can estimate transverse thermal expansion coefficient when
*θ* are given.

### 3.2 Comparison of the Proposed Formula and Experimental Results

Pyroelectric coefficient *p*_{3} is the rate of change in remnant polarization
*θ*, so it can be expressed as

Changes in remnant polarization

Substituting Eq. (2) into Eq. (9) gives Eq. (10) below, in which the variation of remnant polarization
*θ* is determined for given values
*θ*_{0} are given.

In this work, reference temperature *θ*_{0} is 20°C. The values of all coefficients in Eq. (10) are listed in Table 1. Applying Eq. (10) and Table 1 to the data in Fig. 2(a) gives the prediction results shown in Fig. 6. The symbols in Fig. 6 represent the measured data and the curves represent the calculation results by Eq. (10). Fig. 6(a) shows the plots of changes in remnant polarization
*θ* and Fig. 6(b) the plots of pyroelectric coefficient *p*_{3} with *θ*. In Fig. 6(a), the measurements and predictions of remnant polarization

Experimental formulae for longitudinal and transverse remnant strains
_{2} and (6)_{2}. The values of all coefficients for the strain formulae are found in Table 1. Fig. 7 shows the calculation results for strains. Fig. 7(a) shows the plots for longitudinal remnant strain
*θ*.
*α*_{3} and Fig. 7(d) transverse thermal expansion coefficient *α*_{1} plotted versus temperature *θ*. The agreement of calculated thermal expansion coefficients with measured ones is relatively good.

## 4. Conclusions

In the present work, fourteen different electric fields were applied to poled piezoelectric ceramic specimens in the direction opposite to polarization in order to reach a specific state of polarization and strains. Then electric field was removed, and the temperature of specimen was increased, along with that of the invar specimens. Pyroelectric coefficients were calculated using remnant polarizations measured during temperature increase, and strains using thermal expansion coefficient. Then, experimental formulae were developed and used to predict pyroelectric coefficient and strains during temperature rises. The predicted results were compared with the experimental results. For all fourteen states investigated, polarization and strain variations were predicted accurately. The predictions for pyroelectric and thermal expansion coefficients were in relatively good agreement with experimental results. Thus this study shows that the suggested experiment-based formulae can be used to predict the macroscopic behavior, i.e., polarization and strain changes, of ferroelectric materials during temperature increases, including changes in thermal properties such as pyroelectric and thermal expansion coefficients.

## Acknowledgments

This research was supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF), funded by the Ministry of Science, ICT and Future Planning (2015R1A2A2A01005067).