On the fracture toughness of ferroelastic materials

The toughness enhancement due to domain switching near a steadily growing crack in a ferroelastic material is analyzed. The constitutive response of the material is taken to be characteristic of a polycrystalline sample assembled from randomly oriented tetragonal single crystal grains. The constitutive law accounts for the strain saturation, asymmetry in tension versus compression, Bauschinger effects, reverse switching, and strain reorientation that can occur in these materials due to the non-proportional loading that arises near a propagating crack. Crack growth is assumed to proceed at a critical level of the crack tip energy release rate. Detailed finite element calculations are carried out to determine the stress and strain fields near the growing tip, and the ratio of the far field applied energy release rate to the crack tip energy release rate. The results of the finite element calculations are then compared to analytical models that assume the linear isotropic K-field solution holds for either the near tip stress or strain field. Ultimately, the model is able to account for the experimentally observed toughness enhancement in ferroelastic ceramics.     

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.     

A phenomenological constitutive model for ferroelastic and ferroelectric hysteresis effects in ferroelectric ceramics

This paper is concerned with a macroscopic constitutive law for domain switching effects, which occur in ferroelectric ceramics. The three-dimensional model is thermodynamically consistent and is determined by two scalar valued functions: the Helmholtz free energy and a switching surface. In a kinematic hardening process the movement of the center of the switching surface is controlled by internal variables. In common usage, the remanent polarization and the irreversible strain are employed as internal variables. The novel aspect of the present work is to introduce an irreversible electric field, which serves instead of the remanent polarization as internal variable. The irreversible electric field has only theoretical meaning, but it makes the formulation very suitable for a finite element implementation, where displacements and the electric potential are the nodal degrees of freedom. The paper presents an appropriate implementation into a hexahedral finite brick element. The uni-axial constitutive model successfully reproduces the ferroelastic and the ferroelectric hysteresis as well as the butterfly hysteresis for ferroelectric ceramics. Furthermore it accounts for the mechanical depolarization effect, which occurs if the polarized ferroelectric ceramic is subjected to a compression stress.     

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.     

Fully coupled, multi-axial, symmetric constitutive laws for polycrystalline ferroelectric ceramics

In this paper, a general form for multi-axial constitutive laws for ferroelectric ceramics is constructed. The foundation of the theory is an assumed form for the Helmholtz free energy of the material. Switching surfaces and associated flow rules are postulated in a modified stress and electric field space such that a positive dissipation rate during switching is guaranteed. The resulting tangent moduli relating increments of stress and electric field to increments of strain and electric displacement are symmetric since changes in the linear elastic, dielectric and piezoelectric properties of the material are included in the switching surface. Finally, parameters of the model are determined for two uncoupled cases, namely non-remanent straining ferroelectrics and purely ferroelastic switching, and then for the fully coupled ferroelectric case.     

含随机分布缺陷的多晶铁电材料的有效电弹性能

针对铁电材料含有随机分布缺陷以及微结构在外场作用下发生变化的特点,建立起一个细观统计力学模型,考虑到缺陷和酶极化转动之间的相互影响,利用细观力学方法－Eshelby等效夹杂法和Mori-Tanaka的平均理论,具体分析了所含缺陷以及外场对铁电材料有效电弹性能和模量的影响。针对BaTiO3铁电陶瓷的有效电弹性能与常数的预报结果表明了缺陷的存在将增强材料的压电性能。     

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.     

Macroscopical non-linear material model for ferroelectric materials inside a hybrid finite element formulation

A new approach for modeling hysteretic non-linear ferroelectric ceramics is presented, based on a fully ferroelectric/ferroelastic coupled macroscopic material model. The material behavior is described by a set of yield functions and the history dependence is stored in internal state variables representing the remanent polarization and the remanent strain. For the solution of the electromechanical coupled boundary value problem, a hybrid finite element formulation is used. Inside this formulation the electric displacement is available as nodal quantity (i.e. degree of freedom) which is used instead of the electric field to determine the evolution of remanent polarization. This involves naturally the electromechanical coupling. A highly efficient integration technique of the constitutive equations, defining a system of ordinary differential equations, is obtained by a customized return mapping algorithm. Due to some simplifications of the algorithm, an analytical solution can be calculated. The automatic differentiation technique is used to obtain the consistent tangent operator. Altogether this has been implemented into the finite element code FEAP via a user element. Extensive verification tests are performed in this work to evaluate the behavior of the material model under pure electrical and mechanical as well as coupled and multi-axial loading conditions.     

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

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

Finite element simulation of a polycrystalline ferroelectric based on a multidomain single crystal switching model

In this paper, we compute the constitutive behavior of a ferroelectric ceramic by a plane strain finite element model, where each element represents a single grain in the polycrystal. The properties of a grain are described by the microscopic model for switching in multidomain single crystals of ferroelectric materials presented by Huber et al. [J. Mech. Phys. Solids 47 (1999) 1663]. The poling behavior of the polycrystal is obtained by employing the finite element formulation for electromechanical boundary value problems developed by Landis [Int. J. Numer. Meth. Eng. 55 (2002) 613]. In particular, we address the influence of the single grain properties and the interaction between grains, respectively.     

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.     

A unified framework for modeling hysteresis in ferroic materials

This paper addresses the development of a unified framework for quantifying hysteresis and constitutive nonlinearities inherent to ferroelectric, ferromagnetic and ferroelastic materials. Because the mechanisms which produce hysteresis vary substantially at the microscopic level, it is more natural to initiate model development at the mesoscopic, or lattice, level where the materials share common energy properties along with analogous domain structures. In the first step of the model development, Helmholtz and Gibbs energy relations are combined with Boltzmann theory to construct mesoscopic models which quantify the local average polarization, magnetization and strains in ferroelectric, ferromagnetic and ferroelastic materials. In the second step of the development, stochastic homogenization techniques are invoked to construct unified macroscopic models for nonhomogeneous, polycrystalline compounds exhibiting nonuniform effective fields. The combination of energy analysis and homogenization techniques produces low-order models in which a number of parameters can be correlated with physical attributes of measured data. Furthermore, the development of a unified modeling framework applicable to a broad range of ferroic compounds facilitates material characterization, transducer development, and model-based control design. Attributes of the models are illustrated through comparison with piezoceramic, magnetostrictive and shape memory alloy data and prediction of material behavior.     

铁电材料的疲劳失效行为

在航空航天、核能发电等重大装备技术领域, 作为高温传感/驱动/能量收集器件的敏感材料——铋层状结构铁电(BLSF)陶瓷在复杂载荷环境下的疲劳失效问题严重限制着器件寿命和可靠性的提高. 本文以BLSF陶瓷的应用需求为背景, 围绕铁电材料的疲劳裂纹扩展与电畴极化翻转及其相互作用机制等关键问题, 综述了铁电材料在热、力、电三种载荷及其耦合作用下疲劳失效行为的研究现状, 并根据当前铁电材料的一些新发展、新应用对其未来研究方向进行了展望, 旨在为高性能、长寿命铁电/压电器件设计提供参考.      

多孔极化PZT95/5铁电陶瓷单轴压缩力学响应与放电特性

采用添加造孔剂的方法制备了四种不同孔隙率PZT95/5铁电陶瓷,对其进行电场极化,随后开展了准静态单轴压缩实验,讨论了畴变、相变以及孔隙率对极化PZT95/5铁电陶瓷的力学响应与放电特性的影响. 研究结果表明：(1)多孔极化PZT95/5铁电陶瓷非线性力学响应行为主要归因于畴变和相变的共同作用,与内部孔洞变形和坍塌基本无关;(2) 在准静态单轴压缩下极化PZT95/5铁电陶瓷的去极化机制是畴变和相变的共同作用;(3) 孔隙率对极化PZT95/5铁电陶瓷的弹性模量、压缩强度有明显的影响,而对断裂应变的影响较小;(4)极化PZT95/5铁电陶瓷畴变和相变开始的临界应力都随着孔隙率的增大而线性衰减,但相变开始的临界体积应变却不依赖孔隙率;(5)极化PZT95/5铁电陶瓷电荷饱和释放量随着孔隙率呈线性减小,但孔隙率对电荷释放速率基本没有影响。     

Computational modeling of rate-dependent domain switching in piezoelectric materials

In this contribution a micromechanically motivated model for rate-dependent switching effects in piezoelectric materials is developed. The proposed framework is embedded into a three-dimensional finite element setting whereby each element is assumed to represent an individual grain. Related dipole (polarization) directions are thereby initially randomly oriented at the element level to realistically capture the originally un-poled state of grains in the bulk ceramics. The onset of domain switching processes is based on a representative energy criterion and combined with a linear kinetics theory accounting for time-dependent propagation of domain walls during switching processes. In addition, grain boundary effects are incorporated by making use of a macromechanically motivated probabilistic approach. Standard volume-averaging techniques with respect to the response on individual grains in the bulk ceramics are later on applied to obtain representative hysteresis and butterfly curves under macroscopically uniaxial loading conditions at different loading frequencies. It turns out that the simulations based on the developed finite element formulation nicely match experimental data reported in the literature.     

A Micromechanics-based Finite Element Model for the Constitutive Behavior of Polycrystalline Ferromagnets

A micromechanics-based finite element model for the constitutive behavior of polycrystalline ferromagnets is developed. In the model, the polycrystalline solid is assumed to comprise numerous single crystals with randomly distributed crystallographic orientations, and the single crystals, in turn, consist of ferromagnetic domains, each of which is represented by a cubic element. The dipole directions of the domains are randomly assigned to simulate the crystallographic nature of ferromagnetic polycrystals. A switching criterion for the domains is specified at the microscopic level. The macroscopic constitutive behavior is obtained by averaging the microscopic/local behavior of each domain. The developed model has been applied to the simulation of a ferromagnetic material. With appropriate material parameters adopted, hysteresis loops of the predicted magnetic induction versus magnetic field and those of the strain versus magnetic field are shown to agree well with experimental observations.The project supported by the National Natural Science Foundation of China (90205030, 10472088, 10425210), the National Basic Research Program of China (2006CB601202) and the State Administration of the Foreign Experts Affairs Through the “111” Project (B06024) The English text was polished by Yunming Chen.     

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)     

Multi-grain analysis versus self-consistent estimates of ferroelectric polycrystals

Ferroelectrics are polycrystalline materials consisting of intragranular regions with different polarization directions, called domains. The domains can be switched into different states by the application of an electric field or mechanical stress. We study the influence of grain-to-grain interactions on the overall and local switching behavior. The behavior inside each grain is represented by the micromechanics model of Huber et al. [1999. A constitutive model for ferroelectric polycrystals. J. Mech. Phys. Solids 47 (8), 1663-1697]. The predictions of a self-consistent model of the polycrystal response are compared with those of a multi-grain model in which grains are represented individually. In one flavor of the multi-grain model, each grain is represented by a single finite element, while in the other the fields inside each grain are captured in more detail through a fine discretization. Different electrical and mechanical loading situations are investigated. It is found that the overall response is only mildly dependent on the accuracy with which grain-to-grain interactions are captured, while the distribution of grain-average stresses is quite sensitive to the resolution of the intragranular fields.