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.     

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.     

Study on fracture behavior of ferroelectric ceramics under combined electromechanical loading by using a moiré interferometry technique

This paper discusses an in situ observation of fracture behavior around a crack tip in ferroelectric ceramics under combined electromechanical loading by use of a moiré interferometry technique. The deformation field induced by the electric field and the stress concentration near the crack tip in three-points bending experiments was measured. By analysis of the moiré images it is found that under a constant mechanical load, the electric field almost has no effect on the crack extension in the case that the directions of the poling, electric field and crack extension are perpendicular to each other. When the poling direction is parallel to the crack extension direction and perpendicular to the electric field, the strain decreases faster than that calculated by FEM with and without electrical loading as one goes away from the crack tip. In addition, as the electric field intensity increases, the strain near the crack tip increases, and the strain concentration becomes more significant. The project supported by the National Natural Science Foundation of China (10132010, 10025209, 10232023)     

Multiscale behavior of crack initiation and growth in piezoelectric ceramics

The multiscale nature of cracking in ferroelectric ceramics is explored in relation to the crack growth enhancement and retardation behavior when the direction of applied electric field is reversed with reference to that of poling. An a priori knowledge of the prevailing fracture behavior is invoked for the energy dissipated in exchange of the macro- and micro-crack surface. To avoid the formalism of developing a two-scale level model, a single dominant crack is considered where the effect of microcracking could be reflected by stable crack growth prior to macro-crack instability. This is accounted for via a length ratio parameter λ. Micro- and macro-crack damage region would necessarily overlap in the simplified approach of applying equilibrium mechanics solutions to different scale ranges that are connected only on the average over space and time. The strain energy density theory is applied to determine the crack growth segments for conditions of positive, negative and zero electric field. The largest and smallest crack segments were found to correspond, respectively, to the positive and negative field. All of the three piezoceramics PZT-4, PZT-5H and P-7 followed such a trend. This removes the present-day controversy arising from the use of the energy release rate concept that yields results independent of the sign of the electric field. Interaction of non-similar crack growth with the direction of electric field is also discussed in relation to Mode II cracking. The crack initiation angle plays a dominant role when the growth segment is sufficiently small. Otherwise, a more complex situation prevails where consideration should also be given to the growth segment length. Failure stresses of Modes I and II cracking are also obtained and they are found to depend not only on the electric field density but also on crack length and the extent of slow crack growth damage. These findings suggest a series of new experiments.     

Analysis of electric boundary condition effects on crack propagation in piezoelectric ceramics

There are three types of cracks: impermeable crack, permeable crack and conducting crack, with different electric boundary conditions on faces of cracks in piezoelectric ceramics, which poses difficulties in the analysis of piezoelectric fracture problems. In this paper, in contrast to our previous FEM formulation, the numerical analysis is based on the used of exact electric boundary conditions at the crack faces, thus the common assumption of electric impermeability in the FEM analysis is avoided. The crack behavior and elasto-electric fields near a crack tip in a PZT-5 piezoelectric ceramic under mechanical, electrical and coupled mechanical-electrical loads with different electric boundary conditions on crack faces are investigated. It is found that the dielectric medium between the crack faces will reduce the singularity of stress and electric displacement. Furthermore, when the permittivity of the dielectric medium in the crack gap is of the same order as that of the piezoelectric ceramic, the crack becomes a conducting crack, the applied electric field has no effect on the crack propagation. The project supported by the National Natural Science Foundation of China (19672026, 19891180)     

Unconventional domain band structure at the crack tip in ferroelectric ceramics

Domain polarization switch near the tip of a crack or an electrode plays a critical role in the fracture or toughening of ferroelectric ceramics. The intensive electric field near a crack tip stimulates local domain switching. Experiment indicates that the domain band structure in front of an indentation crack under lateral electric loading is unconventional, attributed to the highly localizing crack tip electric field. The partially switched ferroelectric grain resembles a banded Eshelby inclusion embedded in a polycrystalline ferroelectric matrix. The domain wall energy for unconventional domain structures is estimated via arrays of misfit dislocations. Mesomechanics analysis quantifies the unconventional domain band structures. The predicted parameters include the volume fraction, the thickness, and the orientation of switched domain bands.     

The interaction between crack and electric dipole of piezoelectricity

Discrete dipoles located near the crack tip play an important role in nonlinear electric field induced fracture of piezoelectric ceramics. A physico-mathematical model of dipole is constructed of two generalized concentrated piezoelectric forces with equal density and opposite sign. The interaction between crack and electric dipole in piezoelectricity is analyzed. The closed form solutions, including those for stress and electric displacement, crack opening displacement and electric potential, are obtained. The function of piezoelectric anisotropic direction,p _a (θ)=cosθ+p _a sinθ, can be used to express the influence of a dipole's direction. In the case that a dipole locates near crack tip, the piezoelectric stress intensity factor is a power function with −3/2 index of the distance between dipole and crack tip. Supported by National Natural Science Foundation of China(No. 10072033)     

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

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

Effect of electric boundary conditions on crack propagation in ferroelectric ceramics

In this paper, the effect of electric boundary conditions on Mode I crack propagation in ferroelectric ceramics is studied by using both linear and nonlinear piezoelectric fracture mechanics. In linear analysis, impermeable cracks under open circuit and short circuit are analyzed using the Stroh formalism and a rescaling method. It is shown that the energy release rate in short circuit is larger than that in open circuit. In nonlinear analysis, permeable crack conditions are used and the nonlinear effect of domain switching near a crack tip is considered using an energy-based switching criterion proposed by Hwang et al.（Acta Metal. Mater.,1995）. In open circuit, a large depolarization field induced by domain switching makes switching much more diffcult than that in short circuit. Analysis shows that the energy release rate in short circuit is still larger than that in open circuit, and is also larger than the linear result. Consequently,whether using linear or nonlinear fracture analysis, a crack is found easier to propagate in short circuit than in open circuit, which is consistent with the experimental observations of Kounga Njiwa et al.（Eng. Fract. Mech., 2006）.     

Numerical analysis of quasi-static crack branching in brittle solids by a modified displacement discontinuity method

Mechanism of quasi-static crack branching in brittle solids has been analyzed by a modified displacement discontinuity method. It has been assumed that the pre-existing cracks in brittle solids may propagate at the crack tips due to the initiation and propagation of the kink (or wing) cracks. The originated wing cracks will act as new cracks and can be further propagated from their tips according to the linear elastic fracture mechanics (LEFM) theory. The kink displacement discontinuity formulations (considering the linear and quadratic interpolation functions) are specially developed to calculate the displacement discontinuities for the left and right sides of a kink point so that the first and second mode kink stress intensity factors can be estimated. The crack tips are also treated by boundary displacement collocation technique considering the singularity variation of the displacements and stresses near the crack tip. The propagating direction of the secondary cracks can be predicted by using the maximum tangential stress criterion. An iterative algorithm is used to predict the crack propagating path assuming an incremental increase of the crack length in the predicted direction (straight and curved cracks have been treated). The same approach has been used for estimating the crack propagating direction and path of the original and wing cracks considering the special crack tip elements. Some example problems are numerically solved assuming quasi-static conditions. These results are compared with the corresponding experimental and numerical results given in the literature. This comparison validates the accuracy and applicability of the proposed method.     

Phase-field modeling of crack propagation in piezoelectric and ferroelectric materials with different electromechanical crack conditions

We present a family of phase-field models for fracture in piezoelectric and ferroelectric materials. These models couple a variational formulation of brittle fracture with, respectively, (1) the linear theory of piezoelectricity, and (2) a Ginzburg–Landau model of the ferroelectric microstructure to address the full complexity of the fracture phenomenon in these materials. In these models, both the cracks and the ferroelectric domain walls are represented in a diffuse way by phase-fields. The main challenge addressed here is encoding various electromechanical crack models (introduced as crack-face boundary conditions in sharp models) into the phase-field framework. The proposed models are verified through comparisons with the corresponding sharp-crack models. We also perform two dimensional finite element simulations to demonstrate the effect of the different crack-face conditions, the electromechanical loading and the media filling the crack gap on the crack propagation and the microstructure evolution. Salient features of the results are compared with experiments.     

Coupled solution for anisotropic piezoelectric materials with cracks

The coupled elastic and electric fields for anisotropic piezoelectric materials with electrically permeable cracks are analyzed by using Stroh formula in anisotropic elasticity. It is shown from the solution that the tangent component of the electric field strength and the normal component of the electric displacement along the faces of cracks are all constants, and the electric field intensity and electric displacement have the singularity of type (1/2) at the crack tip. The energy release rate for crack propagation depends on both the stress intensity factor and material constants. The electric field intensity and electric displacement inside electrically permeable cracks are all constants.     

Electric-field-induced fatigue crack growth in ferroelectric ceramics

Electric-field-induced fatigue crack growth in pre-cracked PZT ferroelectric ceramics is experimentally investigated in this work. It is found that the crack open and close under an alternating electric field is a major mechanism of crack propagation. The experimental results also show that the frequency, waveform, as well as the amplitude ratio, of the electric loading, play important roles in electric-field-induced fatigue cracking. Empirical formulations of fatigue crack propagation rates are obtained based on the experimental results. It is revealed that the crack grows at a nearly constant rate when the loading frequency is below 100 Hz. However, with the increase of the loading frequency over 125 Hz, the crack propagation rate diminishes rapidly.     

Strip-electromechanical model solution for piezoelectric plate cut along two semi-permeable collinear cracks

A strip electric saturation and mechanical yielding model solution is proposed for a piezoelectric plate cut along two equal collinear semi-permeable mode-I cracks with electrical polarization reaching a saturation limit and normal stress reaching a yield stress along a line segment in front of the cracks. By using Stroh formalism and complex variable technique, we derived the analytical solution for the field quantities. Three different situations are investigated when developed electrical saturation zone is bigger/smaller or equal to the developed mechanical yield zone. Numerical results show that the effect of different electric boundary conditions on the crack opening displacement and crack opening potential drop is significant and should not be ignored. The influence of electric load displacement on the energy release rate is also investigated for PZT-4, PZT-5H and BaTiO₃ ceramics, and it may assists for the correct choosing of ceramic for specific job.     

Effect of electric fields on fracture behavior of ferroelectric ceramics

The asymptotic problem of a semi-infinite crack perpendicular to the poling direction in a ferroelectric ceramic subjected to combined electric and mechanical loading is analyzed to investigate effect of electric fields on fracture behavior. Electromechanical coupling induced by the piezoelectric effect is neglected in this paper. The shape and size of the switching zone is shown to depend strongly on the relative magnitude between the applied electric field and stress field as well as on the ratio of the coercive electric field to the yield electric field. A universal relation between the crack tip stress intensity factor and the applied intensity factors of stress and electric field under small-scale conditions is obtained from the solution of the switching zone. It is found that the ratio of the coercive electric field to the yield electric field plays a significant role in determining the enhancement or reduction of the crack tip stress intensity factor. The fracture toughness variation of ferroelectrics under combined electric and mechanical loading is also discussed.     

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)     

各向异性压电材料平面裂纹的耦合场分析

用Stroh方法分析了各向异性压电材料电导通型裂纹问题的耦合场。结果表明,裂纹面上的切向电场强度和法向电位移均为常数,在裂纹尖端有由弹性场的耦事作用产生的奇异电导通裂纹模型中的静电场对裂纹尖端扩展的能量释放率不作贡献。     

Electrically driven cracks in piezoelectric ceramics: experiments and fracture mechanics analysis

Piezoelectric systems like multilayer actuators are susceptible to damage by crack propagation induced by strain incompatibilities. These can arise under electric fields for example between the electroded and external regions. Such incompatibilities have been realised in thin rectangular model specimens from PZT-piezoelectric ceramics with top and bottom electrodes only close to one edge. Under an electric field, controlled crack propagation has been observed in situ in an optical microscope. The crack paths are reproducible with very high accuracy. Small electrode widths lead to straight cracks with two transitions between stable and unstable crack growth regions, while large electrode widths result in curved cracks with four transitions. Fracture mechanics analysis is able to explain the different crack paths. An iteration method is developed to simulate the curved crack propagation also for strong curvature of the crack paths using the finite element method. The computed crack contours exhibit excellent quantitative agreement with the experiment with respect to their shape, the stages of stable and unstable crack propagation and the transitions between them. Finally, also the crack length as a function of the electric field can be predicted.     

Energy density theory formulation and interpretation of cracking behavior for piezoelectric ceramics

Any attempt made to separate energy into electrical and mechanical parts may lead to inconsistencies as they do not necessarily decouple. This is illustrated by application of the energy density function in the linear theory of piezoelasticity. By assuming that a critical energy density function prevails at the onset of crack initiation, it is possible to establish the relative size of an inner and outer damage zone around the crack tip; they correspond to the ligaments at failure caused by pure electric field and pure mechanical load. On physical grounds, the relative size of these zones must depend on the relative magnitude of the mechanical and electrical load. Hence, they can vary in size depending on the electromechanical material and damage resistance properties. Numerical results are obtained for the PZT-4, PZT-5H, and P-7 piezoelectric ceramics. These two ligaments for the two damage zones may coincide for appropriate values of the applied electrical field and mechanical load.Explicit expression of the energy density factor S is derived showing the mixed mode electromechanical coupling effects. The factor S can increase or decrease depending on the direction of the applied electric field with reference to the poling direction. This is in contrast to the result obtained from the energy release rate quantity, which remains unchanged for electric field in the direction of poling or against it.     

Fracture prediction for piezoelectric ceramics based on the electric field saturation concept

The electrical field saturation model is applied to the fracture prediction of piezoelectric materials containing electrically impermeable cracks. This model is analogously similar to the electric displacement saturation model that available in the literature. An electrical field saturation strip near the crack front is introduced in the analytical model. The stress intensity factor K and the energy release rate G are obtained in closed-form. It is found that fracture predictions based on K and G criteria are identical. Fracture predictions based on the electric field saturation model and the electric displacement model are also found to be the same.