Energy principle of ferroelectric ceramics and single domain mechanical model

Many physical experiments have shown that the domain switching in a ferroelectric material is a complicated evolution process of the domain wall with the variation of stress and electric field. According to this mechanism, the volume fraction of the domain switching is introduced in the constitutive law of ferroelectric ceramic and used to study the nonlinear constitutive behavior of ferroelectric body in this paper. The principle of stationary total energy is put forward in which the basic unknown quantities are the displacement u _i , electric displacement D _i and volume fraction ρ _I of the domain switching for the variant I. Mechanical field equation and a new domain switching criterion are obtained from the principle of stationary total energy. The domain switching criterion proposed in this paper is an expansion and development of the energy criterion. On the basis of the domain switching criterion, a set of linear algebraic equations for the volume fraction ρ _I of domain switching is obtained, in which the coefficients of the linear algebraic equations only contain the unknown strain and electric fields. Then a single domain mechanical model is proposed in this paper. The poled ferroelectric specimen is considered as a transversely isotropic single domain. By using the partial experimental results, the hardening relation between the driving force of domain switching and the volume fraction of domain switching can be calibrated. Then the electromechanical response can be calculated on the basis of the calibrated hardening relation. The results involve the electric butterfly shaped curves of axial strain versus axial electric field, the hysteresis loops of electric displacement versus electric filed and the evolution process of the domain switching in the ferroelectric specimens under uniaxial coupled stress and electric field loading. The present theoretic prediction agrees reasonably with the experimental results given by Lynch. The project supported by the National Natural Science Foundation of China (10572138).     

Two-scale micromechanics-based probabilistic modeling of domain switching in ferroelectric ceramics

A two-scale micromechanics model is developed in this paper to analyze domain switching in ferroelectric ceramics, using a probabilistic domain switching criterion based on energetic analysis. The microstructure of ferroelectric ceramics at two distinct length scales, domains and grains, has been carefully analyzed. The interaction at domain level is accounted for by energy minimization theory, while the fluctuation at grain level is analyzed using ellipsoidal two-point correlation function. The model has been implemented by Monte Carlo method, and applied to simulate the electric poling and mechanical depoling of Pb(Zr_xTi_1-x)O₃ (PZT) ceramics across morphotropic phase boundary (MPB). The drastically different switching characteristics of PZT ceramics across MPB has been captured, and good agreement with experiments has been observed. The effects of the transformation strains and spontaneous polarizations are highlighted, confirming the proposition of Li et al. [2005. Domain switching in polycrystalline ferroelectric ceramics. Nature Materials 4, 776–781] that the strain compatibility plays a dominant role in domain switching in ferroelectric ceramics.     

A micromechanics method to study the effect of domain switching on fracture behavior of polycrystalline ferroelectric ceramics

ＩｎｔｒｏｄｕｃｔｉｏｎＲｅｃｅｎｔｌｙ ,ｔｈｅｆｅｒｒｏｅｌｅｃｔｒｉｃｃｅｒａｍｉｃｓｈａｓｓｕｃｈｅｘｃｅｌｌｅｎｔｃｈａｒａｃｔｅｒｉｓｔｉｃｓｏｆｐｉｅｚｏｅｌｅｃｔｒｉｃｉｔｙａｎｄｐｙｒｏｅｌｅｃｔｒｉｃｉｔｙｅｔｃ .ｔｈａｔｉｔｂｅｃｏｍｅｓｏｎｅｏｆｔｈｅｍｏｓｔｉｍｐｏｒｔａｎｔｆｕｎｃｔｉｏｎａｌｍａｔｅｒｉａｌｓ.Ｆｏｒｉｎｓｔａｎｃｅ ,ｔｈｅｗｉｄｅｌｙａｐｐｌｉｅｄｓｅｎｓｏｒｓ,ｔｒａｎｓｄｕｃｅｒｓａｎｄａｃｔｕａｔｏｒｓｅｔｃ .ａｒｅｍａｄｅｏｆｔｈｅｆｅｒｒｏｅ…     

An I-integral method for crack-tip intensity factor variation due to domain switching in ferroelectric single-crystals

In the present study, an I-integral method is established for solving the crack-tip intensity factors of ferroelectric single-crystals. The I-integral combined with the phase field model is successfully used to investigate crack-tip intensity factor variations due to domain switching in ferroelectricity subjected to electromechanical loadings, which exhibits several advantages over previous methods based on small-scale switching. First, the shape of the switching zone around a crack tip is predicted by the time-dependent Ginzburg–Landau equation, which does not require preset energy-based switching criterion. Second, the I-integral can directly solve the crack-tip intensity factors and decouple the crack-tip intensity factors of different modes based on superimposing an auxiliary state onto an actual state. Third, the I-integral is area-independent, namely, the I-integral is not affected by the integral area size, the polarization distributions, or domain walls. This makes the I-integral applicable to large-scale domain switching. To this end, the electro-elastic field intensity factors of an impermeable crack in PbTiO₃ ferroelectric single crystals are evaluated under electrical, mechanical, and combined loading. The intensity factors obtained by the I-integral agree well with those obtained by the extrapolation technique. From numerical results, the following conclusions can be drawn with respect to fracture behavior of ferroelectrics under large-scale switching. Under displacement controlled mechanical loading, the stress intensity factors (SIFs) decrease monotonically due to the domain switching process, which means a crack tip shielding or effective switching-induced toughening occurs. If an external electric field is applied, the electric displacement intensity factor (EDIF) increases in all cases, i.e., the formed domain patterns enhance the electric crack tip loading. The energy release rate, expressed by the crack-tip J-integral, is reduced by the domain switching in all examples, which underlines the switching-induced-toughening effect. In contrast, under stress controlled load, the SIF evolves due to large-scale switching to a stable value, which is higher than the non-switching initial value, i.e., fracture is promoted in this case.     

Effect of electric displacement saturation on the stress intensity factor for a crack in a ferroelectric ceramic

A crack in a ferroelectric ceramic with perfect saturation under electric loading is analyzed. The boundary of the electric displacement saturation zone ahead of the crack tip is assumed to be ellipse in shape. The shape and size of ferroelectric domain switching zone near a crack tip is determined based on the nonlinear electric theory. The stress intensity factor induced by ferroelectric domain switching under small-scale conditions is numerically obtained as a function of the electric saturation zone parameter and the ratio of the coercive electric field to the yield electric field. It is found that the stress intensity factor increases as the ratio of the semi-axes of the saturation ellipse increases.     

Numerical simulation of domain switching in multilayer ferroelectric actuators

Micromechanical finite element methods are developed based on a nonlinear constitutive model of ferroelectric polycrystals. Electromechanical behaviors ahead of an internal electrode tip are numerically simulated in multilayer ferroelectric actuators. Around the electrode edge, the nonuniform electric field generates a concentration of stress due to the incompatible strain as well as spontaneous strain. The preferred domain switching enhances the concentration of residual stress and may cause the actuators to crack. An electrically permeable crack emanating from an internal electrode is analyzed. A large scale domain switching zone is found in the vicinity of crack tips. The larger the actuating strain and electric field are, the larger the switching zone will be. The size of switching zone even reaches the scale of crack length with increasing electromechanical loading.     

Coalescence of Pore Columns by Domain Switching

The present paper studies the coalescence of pore columns in ferroelectric ceramics driven by back and forth domain switching under cyclic electric field. A finite element method that incorporates mass transfer capacity is formulated to simulate the evolution of point defects subjected to the kinetics of pore surface diffusion and domain wall migration. The merge of point defects provides a mechanism for the vacancy agglomeration that leads to the formation of large pores or microcracks.The project supported by the National Natural Science Foundation of China (10332020 and 10121202), and the National 973 Project (2004CB619304).     

铁电陶瓷微观电畴的成核率模型

大量的实验已经证实电畴翻转是铁电材料非线性和迟滞性本构曲线的根本原因。研究者已经对铁电陶瓷的微观电畴翻转行为进行了大量详细的研究。针对描述电畴成核的物理实验结果和经典的成核率实验数据,为了建立电畴翻转体积分数的演化方程提出了反应微观电畴翻转的成核率模型。针对铁电试样电畴随机分布的情况,应用该模型对铁电陶瓷的非线性本构行为进行了研究。理论结果与实验数据的比较表明,模型能较好的描述铁电材料的非线性本构行为。同时模型所揭示的微观反转的物理本质可进一步的指导宏观唯象模型的改进。     

The influence of mechanical constraint upon the switching of a ferroelectric memory capacitor

The role of mechanical constraint upon the switching response of a ferroelectric thin film memory capacitor is explored. The memory capacitor is represented by a two dimensional ferroelectric island whose non-linear behaviour is modelled by a crystal plasticity constitutive law within the finite element method. The switching response of the device, in terms of remnant charge storage, is determined as a function of geometry and constraint. Various types of constraint on the ferroelectric capacitor are considered, including the presence of a silicon dioxide passivation layer, a silicon substrate and metallic electrodes. The effect of the relative resistance to 90 degree switching and 180 degree switching is also explored in a tetragonal ferroelectric device. Throughout the study, the finite element calculations are compared with the behaviour of a material element subjected to various degrees of mechanical constraint.     

Fatigue crack growth induced by domain switching under electromechanical load in ferroelectrics

Reliability calls for a better understanding of the failure of ferroelectric ceramics. The fracture and fatigue of ferroelectric ceramics under an electric field or a combined electric and mechanical loading are investigated. The small-scale domain-switching model is modified to analyze failure due to fracture and fatigue. Effects of anisotropy and electromechanical load coupling are taken into account. Analytical expressions are obtained for domain-switching regions near the crack tip such that of 90° domain switching can be distinguished from 180° domain switching in addition to different initial poling directions. The crack tip stress intensity variation of ferroelectric ceramics due to the domain switching is analyzed. A positive electric field tends to enhance the propagation of an insulating crack perpendicular to the poling direction, while a negative field impedes it. Fatigue crack growth under various coupling loads and effects of the stress field and electric field on near field stress intensity variation are analyzed. Predicted crack growth versus cyclic electric field agrees well with experiment.     

铁电陶瓷微观电畴演化的成核率模型

大量的实验已经证实电畴翻转是铁电材料非线性和迟滞性本构曲线的根本原因.根据这些描述电畴运动的物理实验结果,提出了成核率模型来模拟铁电材料的非线性本构行为.进一步由模型建立了基于微观电畴运动的电畴翻转体积分数演化方程.应用该模型得到的理论计算结果与实验数据的比较表明,模型能够描述铁电材料的非线性本构行为.同时模型所揭示的微观电畴成核的演化行为可进一步的指导宏观唯象模型的建立与改进.     

Rate dependent properties of perovskite type tetragonal piezoelectric materials using micromechanical model

In this paper, a three-dimensional micromechanical model is presented for simulation of the rate dependent properties of certain perovskite type tetragonal piezoelectric materials. The model is based on linear constitutive, nonlinear domain switching, and linear kinetics theories. The simulation starts with a virgin bulk material of randomly oriented grains. Then the material is electrically loaded with an alternating voltage of various frequencies, which are in the order of 0.01 Hz to 1 Hz. An energy equation in combination with a probability function is used to determine the onset of the domain switching inside the grains. Such a probability function leads to a better phenomenological model for the domain switching even for electrical loadings, which are in a range far below the coercive fields. The propagation of the domain wall during the domain switching process in grains is modeled by means of linear kinetics relations after domain nucleation. The response of the bulk ceramic is predicted by averaging the response of individual grains using Euler angles for the transformation from local coordinates of the grains to global coordinate. Electric displacement hysteresis loops for different frequencies and amplitudes of the alternating electric fields are simulated. A simple micromechanical model without the probabilistic approach is compared with the one that takes it into account. Both models give important insights into the rate dependency of piezoelectric materials, which was observed in some experiments reported in the literature.     

Highly reliable bipolar resistive switching in sol-gel derived lanthanum-doped PbTiO3 thin film： Coupling with ferroelectricity？

Nanoscale PbxLa1-,Ti1-x/4O3 （PLT） thin film has been fabricated on Pt／Ti／SiO2／Si substrates by chemical solution deposition （CSD） method. Ferroelectricity of the fresh-made PLT thin film has been clearly detected through piezoelectric force microscopy （PFM） by writing reversible ferroelectric domains. However, PLT thin film also shows off-standard ferroelectric hysteresis loops highly dependent on frequency, indicating large amount of mobile space charges in the film. Subsequent current-voltage （C-V） studies show that sandwich-like Pt／PLT／Pt structure exhibits notable bipolar resistive switching （BRS） characteristics with high stability （〉 103 switching cycles）. It is found that the C-V curves of both high- and low-resistance states have the feature of space-charge-limited current （SCLC） conduction, indicating important roles of defects in the conduction. X-ray photoelectron spectroscopy measurement further verifies that oxygen vacancies based conductive filament mechanism is likely responsible for the observed RS effect. Our demonstration of stable RS effect in the PLT thin film and its possible coupling with ferroelectricity is promising in device development and applications, such as development of ferroelectric-tunable RS memories.     

Electrical body forces and electrical tractions in the nonlinear response of ferroelectric actuators

The influence of the electrical body forces and electrical tractions on the nonlinear response of ferroelectric stack actuators is analytically investigated. While the role of the electrical body forces and tractions in the response of piezoelectric actuators is well documented (and in many cases is not significant), the questions of their effect on ferroelectric active materials is still of interest. To examine this influence, the analytical model for the electro-mechanical behavior of a ferroelectric stack actuator is augmented to account for the electrical body forces along the actuator and the electrical tractions at the material–electrode interfaces. Focusing on the effect of the electrical forces and tractions on the ferroelectric domain switching phenomenon, the model is used for the numerical analysis of a ferroelectric layer and for the comparison with the case that neglects the electrical body forces and traction. The comparison theoretically designates cases in which the effect of the electrical body forces and tractions may be prominent and other cases where the classical approach that neglects these effects can be adopted.     

循环电载下热效应对铁电材料畴变的影响

铁电材料在力、电载荷作用尤其是循环载荷作用下有明显的热效应.热效应会在铁电材料中引起应力场和电场.另外,在外载作用下铁电材料发生约束畴变时也会产生附加的应力场,这些附加的应力场和电场都会对铁电材料的畴变产生影响.而且在循环载荷作用的情况下,这种影响会逐渐累积.但在以往的研究中,很少涉及这种影响.该文就循环电载下热效应引起的力、电场和约束畴变时产生的附加应力场对铁电材料畴变的影响进行了初步研究.     

Time dependence of effective properties of polycrystalline ferroelectric ceramics

In this paper, based on Merz[7] experimental results and classical nucleation theory, a micromechanics statistical model is proposed to describe the relation between the special microstructure-level evolution phenomena-domain switching and macro-response. The polycrystalline ferroelectric ceramics treated as a composition of switched domain and unswitched domain, the approaches of Eshelby's equivalent inclusion and Mori-Tanaka's mean field theory are used to analyze and predict its effective electroelastic properties. The model can incorporate the effects of time dependence of domain switching and shape of individual crystalline. To the BaTiO₃ polycrystalline ceramics, the analytical results are in good agreement with the experimental results.     

Analysis of the electromechanical behavior of ferroelectric ceramics based on a nonlinear finite element model

A nonlinear finite element (FE) model based on domain switching was proposed to study the electromechanical behavior of ferroelectric ceramics. The incremental FE formulation was improved to avoid any calculation instability. The problems of mesh sensitivity and convergence, and the efficiency of the proposed nonlinear FE technique have been assessed to illustrate the versatility and potential accuracy of the said technique. The nonlinear electromechanical behavior, such as the hysteresis loops and butterfly curves, of ferroelectric ceramics subjected to both a uniform electric field and a point electric potential has been studied numerically. The results obtained are in good agreement with those of the corresponding theoretical and experimental analyses. Furthermore, the electromechanical coupling fields near (a) the boundary of a circular hole, (b) the boundary of an elliptic hole and (c) the tip of a crack, have been analyzed using the proposed nonlinear finite element method (FEM). The proposed nonlinear electromechanically coupled FEM is useful for the analysis of domain switching, deformation and fracture of ferroelectric ceramics.The project supported by the National Natural Science Foundation of China (10025209, 10132010 and 90208002), the Research Grants of the Council of the Hong Kong Special Administrative Region, China (HKU7086/02E) and the Key Grant Project of the Chinese Ministry of Education (0306)     

A ferroelectric and ferroelastic 3D hexahedral curvilinear finite element

An isoparametric 3D electromechanical hexahedral finite element integrating a 3D phenomenological ferroelectric and ferroelastic constitutive law for domain switching effects is proposed. The model presents two internal variables which are the ferroelectric polarization (related to the electric field) and the ferroelastic strain (related to the mechanical stress). An implicit integration technique of the constitutive equations based on the return-mapping algorithm is developed. The mechanical strain tensor and the electric field vector are expressed in a curvilinear coordinate system in order to handle the transverse isotropy behavior of ferroelectric ceramics. The hexahedral finite element is implemented into the commercial finite element code Abaqus® via the subroutine user element. Some linear (piezoelectric) and non linear (ferroelectric and ferroelastic) benchmarks are considered as validation tests.     

Three-dimensional finite element simulations of ferroelectric polycrystals under electrical and mechanical loading

Complex, non-linear, irreversible, hysteretic behavior of polycrystalline ferroelectric materials under a combined electro-mechanical loading is a result of domain wall motion, causing simultaneous expansion and contraction of unlike domains, grain sub-divisions that have distinct spontaneous polarization and strain. In this paper, a 3-dimensional finite element method is used to simulate such a polycrystalline ferroelectric under electrical and mechanical loading. A constitutive law due to Huber et al. [1999. A constitutive model for ferroelectric polycrystals. J. Mech. Phys. Solids 47, 1663-1697] for switching by domain wall motion in multidomain ferroelectric single crystals is employed in our model to represent each grain, and the finite element method is used to solve the governing conditions of mechanical equilibrium and Gauss's law. The results provide the average behavior for the polycrystalline ceramic. We compare the outcomes predicted by this model with the available experimental data for various electromechanical loading conditions. The qualitative features of ferroelectric switching are predicted well, including hysteresis and butterfly loops, the effect on them of mechanical compression, and the response of the polycrystal to non-proportional electrical loading.     

A constitutive model for thermo-electro-mechanical behavior of ferroelectric polycrystals near room temperature

A constitutive model that can be used to predict thermo-electro-mechanical linear and nonlinear behavior of ferroelectric polycrystals near room temperature is proposed. A ferroelectric polycrystal is modeled by an agglomerate of 210 single crystallites that are distributed regularly over all directions. A variant in a single crystallite is characterized by a Gibbs free energy function whose coefficients have linear dependency on temperature. A dissipation inequality for domain switching is derived from the restriction of the second law of thermodynamics. Domain switching process is governed by a viscoplastic switching law with temperature-dependent switching parameters. The responses of the proposed model to electric field and mechanical stress loading at room and elevated temperatures are calculated and compared qualitatively with experimental observations available in literature.