Numerical solution of polarization saturation/dielectric breakdown model in 2D finite piezoelectric media

The polarization saturation (PS) model [Gao, H., Barnett, D.M., 1996. An invariance property of local energy release rates in a strip saturation model of piezoelectric fracture. Int. J. Fract. 79, R25–R29; Gao, H., Zhang, T.Y., Tong, P., 1997. Local and global energy release rates for an electrically yielded crack in a piezoelectric ceramic. J. Mech. Phys. Solids 45, 491–510], and the dielectric breakdown (DB) model [Zhang, T.Y., Zhao, M.H., Cao, C.F., 2005. The strip dielectric breakdown model. Int. J. Fract. 132, 311–327] explain very well some experimental observations of fracture of piezoelectric ceramics. In this paper, the nonlinear hybrid extended displacement discontinuity-fundamental solution method (NLHEDD-FSM) is presented for numerical analysis of both the PS and DB models of two-dimensional (2D) finite piezoelectric media under impermeable and semi-permeable electric boundary conditions. In this NLHEDD-FSM, the solution is expressed approximately by a linear combination of fundamental solutions of the governing equations, which includes the extended point force fundamental solutions with sources placed at chosen points outside the domain of the problem under consideration, and the extended Crouch fundamental solutions with extended displacement discontinuities placed on the crack and the electric yielding zone. The coefficients of the fundamental solutions are determined by letting the approximated solution satisfy certain conditions on the boundary of the domain, on the crack face and the electric yielding zone. The zero electric displacement intensity factor in the PS model or the zero electric field strength intensity factor in the DB model at the outer tips of the electric yielding zone is used as a supplementary condition to determine the size of the electric yielding zone. Iteration approaches are adopted in the NLHEDD-FSM. The electric yielding zone is determined, and the extended intensity factors and the local J-integral are calculated for center cracks in piezoelectric strips. The effects of finite domain size, saturation property and different electric boundary conditions, as well as different models on the electric yielding zone and the local J-integral, are studied.     

Electrical nonlinear anti-plane shear crack in a functionally graded piezoelectric strip

The electrical nonlinear behavior of an anti-plane shear crack in a functionally graded piezoelectric strip is studied by using the strip saturation model within the framework of linear electroelasticity. The analysis is conducted on the electrically unified crack boundary condition with the introduction of the electric crack condition parameter that can describe all the electric crack boundary condition in accordance with the aspect ratio of an ellipsoidal crack and the permittivity inside the crack, in particular, including traditional permeable and impermeable crack boundary conditions. The resulting mixed boundary value problem is analysed and near tip field is obtained by using the integral transform techniques. Numerical results for the normalized five kinds of energy release rates under the small scale electric saturation condition are presented and compared to show the influences of the electric crack condition parameter with the variation of the ellipsoidal crack parameters, electric loads, functionally graded piezoelectric material gradation, crack length, electromechanical coupling coefficient, and crack location. It reveals that there are considerable differences between the results obtained from the traditional electric crack models and those obtained from the current unified crack model.     

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)     

Strip electric-magnetic breakdown model in a magnetoelectroelastic medium

Extending the polarization saturation model [Gao et al., 1997. Local and global energy release rates for an electrically yielded crack in a piezoelectric ceramic. J. Mech. Phys. Solids 45, 491-510] and the dielectric breakdown (DB) model [Zhang et al., 2005. The strip dielectric breakdown model. Int. J. Fract. 132, 311-327] in piezoelectric materials, the Strip Electric-Magnetic Breakdown (SEMB) model is proposed for electrically and magnetically impermeable crack in a magnetoelectroelastic medium to study the effect of the nonlinear character of electric field and magnetic field on fracture of magnetoelectroelastic materials. In the SEMB model, the electric field in the strip of the electric breakdown zone ahead of the crack tip is equal to the electric breakdown strength, while the magnetic filed in the strip of the magnetic breakdown zone is equal to the magnetic breakdown strength. By using the extended Stroh formalism and the extended dislocation modeling of a crack, the Griffith crack problem under the electrically and magnetically elastic-plastic condition in a magnetoelectroelastic medium is reduced to a set of dual integral equations. The sizes of the electric breakdown zone and the magnetic breakdown zone, the extended intensity factors and the local J-integral are obtained. The effect of the combined mechanical-electric-magnetic loadings on the local J-integral is studied.     

PZT-4紧凑拉伸试样的断裂分析

基于线性压电材料的复势理论,通过解析分析,导出了一种分析有限压电板裂纹问题的解析数值方法. 首先,计算了含中心裂纹有限板的断裂参数,与Woo和Wang的解析数值法(Int J Fract, 1993, 62: 203$\sim$218)相比较,表明该方法具有很高的精度和很好的计算效率. 随后,采用该方法和有限元法计算了PZT-4紧凑拉伸试样在绝缘裂纹面边界条件下断裂时的断裂参数,发现各断裂参数的临界值分散性很大,不能作为压电材料的单参数断裂准则. 进而,针对试样真实的裂隙形状,采用有限元法计算了裂隙尖端的应力、电位移场,比较了裂隙内介质的介电性能对裂隙尖端场的影响,计算了带微裂纹的真实裂隙模型的断裂参数并进行了理论分析.      

Unified electrical boundary conditions for a crack interacting with a dislocation in piezoelectric media

This paper investigates the interaction problem between a dislocation and a finite crack in piezoelectric media. Analytical solutions for the generalized two-dimensional problem of a dislocation that is interacting with a finite crack in piezoelectric media are formulated via Stroh formalism. The analysis is conducted on the unified electrical crack boundary condition with the introduction of the electric crack condition parameter that can describe all the electric crack boundary conditions. The two ideal crack boundary conditions, namely, the electrically impermeable and permeable crack assumptions are obtained as two special cases for the current solutions. Based on the complex variable method and the perturbation technique, closed form solutions are obtained. The field intensity factors at the crack tip and the image forces on the dislocation due to the crack are computed and discussed.     

Investigation of the crack problem in non-local piezoelectric materials under combined electromechanical loadings

The present investigation of the crack problem in piezoelectric materials is performed based on the non-local theory. After some manipulations, the impermeable crack, the permeable crack （the crack gap is full of NaCI solution）, and the semi-permeable crack （the crack gap is full of air or silicon oil） are reduced to a uniform formulation by assuming the normal electric displacement on the crack surfaces to be an unknown variable. Thus, a triple integral equation with the unknown normal electric displacement is established. By using the Newton iterative method and solving the triple integral equation, it is found that the normal electric displacement on the crack surfaces is no longer a constant as determined by previous studies, rather, it depends upon the remote combined electromechanical loadings. Numerical results of the stresses and electric displacement fields show that there are no singularities at the crack tips so that the stresses remain finite. It is of great significance that the concrete electric boundary condition on the crack surfaces exerts significant influence on the near-tip fields and in this way plays an important role in evaluating the crack stability in the non-local piezoelectric materials. More specifically, the impermeable crack model always overestimates the finite stresses at the crack tips, whereas the permeable crack model always underestimates them.     

The effective electroelastic property of piezoelectric media with parallel dielectric cracks

Existing studies on the coupled electroelastic behaviour of cracked piezoelectric media have been based mostly on the electrically impermeable and permeable crack models. The current paper presents a study of the effective electroelastic property of piezoelectric media weakened by parallel cracks using a dielectric crack model with the electric boundary condition along the crack surfaces being governed by the opening displacement. The theoretical formulation is obtained using the dilute model of distributed cracks and the solution of a single dielectric crack problem. It is observed that the effective electroelastic property of cracked piezoelectric media is nonlinear and sensitive to loading conditions. Different modes of crack deformation are predicted and discussed. Attention is paid to the transition between electrically permeable and impermeable crack models.     

An exact solution of crack problems in piezoelectric materials

IntroductionInrecentyearscrackproblemsinpiezoelectricmaterialhavereceivedmuchattention.Manytheoreticalanalyseshavebeengivenby[1～16].Itshouldbe,however,notedthatalltheaboveanalysesarebasedonaso-calledimpermeablecrackassumphon,i.e.thecrackfacesareassumedtobeimpermeabletoelectricfield,sotheelectricdisplacementvanishesinsidethecrack.Usingthisassumption,onewillobtainthefollowingresultS[2'3'5,6'9'16]=whentheelectricloadsaresolelyaPPliedatinLfinity,theelectricdisplacementissquare-rootsingularatthe…     

Study of a propagating finite crack in functionally graded piezoelectric materials considering dielectric medium effect

This paper provides a study of the problem of a propagating finite crack under in-plane loading in functionally graded piezoelectric materials (FGPMs). The analytical formulations are developed by Fourier transforms and the resulting singular integral equations are solved by using Chebyshev polynomials. By using a dielectric crack model with deformation-dependent electric boundary condition, numerical simulations are made to show the effects of the dielectric medium, the gradient of material properties and the speed of crack propagation on the fracture parameters, such as the stress, electric displacement and crack opening displacement intensity factors. A critical state for the electromechanical loading applied to the FGPMs is observed, which determines whether the traditionally impermeable (or permeable) crack model serves as the upper or lower bound for the dielectric model. The validity of this dielectric crack model is also examined by comparing the results of different existing crack models.     

A mixed electric boundary value problem for a two-dimensional piezoelectric crack

In this paper, a mixed electric boundary value problem for a two-dimensional piezoelectric crack problem is presented, in the sense that the crack face is partly conducting and partly impermeable. By the analytical continuation method, the unknown electric charge distributions on the upper and lower conducting crack faces are reduced to two decoupled singular integral equations and then these two equations are converted into algebraic equations to find the full field solution. Though the results suggest that the stress intensity factors at the crack tip are identical to those of conventional piezoelectric materials, but the electric field and electric displacement are related to the electric boundary conditions on the crack faces. The electric field and electric displacement are singular not only at crack tips but also at the junctures between the impermeable part and conducting parts. Numerical results for the variations of the electric field, electric displacement field and J-integral with respect to the normalized impermeable crack length are shown. Some discussions for the energy release rate and the J-integral are made.     

Generalized 2D problem of piezoelectric media containing collinear cracks

The generalized 2D problem in piezoelectric media with collinear cracks is addressed based on Stroh's formulation and the exact electric boundary conditions on the crack faces. Exact solutions are obtained, respectively, for two special cases: one is that a piezoelectric solid withN collinear cracks is subjected to uniform loads at infinity, and the other is that a piezoelectric solid containing a single crack is subjected to a line load at an arbitrary point. It is shown when uniform loads are applied at infinity or on the crack faces that, the stress intensity factors are the same as those of isotropic materials, while the intensity factor of electric displacement is dependent on the material constants and the applied mechanical loads, but not on the applied electric loads. Moreover, it is found that the electric field inside any crack is not equal to zero, which is related to the material properties and applied mechanical-electric loads. The project supported by the National Natural Science Foundation of China (19772004)     

Coulomb traction on a penny-shaped crack in a three dimensional piezoelectric body

The axisymmetric problem of a penny-shaped crack embedded in an infinite three-dimensional (3D) piezoelectric body is considered. A general formulation of Coulomb traction on the crack surfaces can be obtained based on thermodynamical considerations of electromechanical systems. Three-dimensional electroelastic solutions are derived by the classical complex potential theory when Coulomb traction is taken into account and the poling direction of piezoelectric body is perpendicular to the crack surfaces. Numerical results show that the magnitude of Coulomb tractions can be large, especially when a large electric field in connection with a small mechanical load is applied. Unlike the traditional traction-free crack model, Coulomb tractions induced by an applied electric field influence the Mode I stress intensity factor for a penny-shaped crack in 3D piezoelectric body. Moreover, compared to the current model, the traditional traction-free crack model always overestimates the effect of the applied electric load on the field intensity factors and energy release rates, which has consequences for 3D piezoelectric fracture mechanics.     

Crack spacing effect for a piezoelectric cylinder under electro-mechanical loading or transient heating

This paper investigates the fracture problem of a piezoelectric cylinder with a periodic array of embedded circular cracks. An electro-mechanical fracture mechanics model is established first. The model is further used to the thermal fracture analysis of a piezoelectric cylinder subjected to a sudden heating on its outer surface. The temperature field and the associated thermal stresses and electric displacements are obtained and are added to the crack surface to form a mixed-mode boundary value problem for the electro-mechanical coupling fracture. The stress and stress intensities are investigated for the effect of crack spacing. Strength evaluation of piezoelectric materials under the transient thermal environment is made and thermal shock resistance of the medium is given.     

Crack energy density in arbitrary direction for piezoelectrics and its characteristic features in electromechanical mixed mode crack

The concepts of crack energy density (CED) and its derivatives in arbitrary direction were established for piezoelectric material and, keeping their application to mixed mode fracture in mind, the characteristic features of them as fracture parameters were investigated based on the approximate equations for CED and its derivatives. That is, CED and its derivatives in arbitrary direction are defined first and separation into their each mode contribution is made. Subsequently, path independent integral expressions of them are derived, and then using them, approximate equations of each mode contribution of CED are obtained concretely for the case where linear singular solution is known. The resulting equations are then used to investigate the effects of electric field and electrical boundary condition on CED and its derivatives. An infinite piezoelectric plane with a crack inclined with respect to the poling direction is considered as a numerical example. Mode I contribution of mechanical CED is mainly employed as a possible fracture parameter for the study and it was shown that applied electric field significantly influences on fracture parameters especially for the impermeable crack perpendicular to the poling direction. The effect of electric field has the tendency to decrease as crack inclination angle increases. It was also found that, even for the impermeable crack perpendicular to the poling direction, crack propagation could be deviated from self-similar direction under a strong negative electric field, and this fact is qualitatively consistent with an existing experimental observation. For the ideally sharp crack with no width, impermeable and Hao and Shen type boundary conditions are admissible showing qualitative agreement with experimental results, but exact boundary condition is not suitable and finally consistent with permeable boundary condition.     

基于新型裂尖杂交元的压电材料断裂力学研究

提出了一种裂尖邻域杂交元模型，将其与标准杂交应力元结合来求解压电材料裂纹尖 端的奇性电弹场和断裂参数的数值解．裂纹尖端杂交元的建立步骤为：1) 利用高次内插有限元特征法求解特征问题，得到反映裂尖奇异性电弹场状况的特 征值和特征角分布函数；2) 利用广义Hellinger-Reissner变分泛函以及特征问题的解来建立裂尖邻域杂交元模型．该 方法求解电弹场时，摒弃了传统有限元方法中裂尖奇异性场需要借助解析解的做法，也避免 了单纯有限元方法中需要在裂尖端部进行高密度单元划分．采用PZT5板中心裂纹问题 作为考核例，数值结果显示了良好的精确性．作为进一步应用，求解了含中心界面裂纹 的PZT4-PZT5两相压电材料的应力强度因子和电位移强度因子．所有的算例都考虑 了3种裂纹面电边界条件．     

Two parallel limited-permeable mode-I cracks or four parallel limited-permeable mode-I cracks in the piezoelectric materials

In this paper, the basic solutions of two parallel mode-I cracks or four parallel mode-I cracks in the piezoelectric materials were investigated by means of the Schmidt method for the limited-permeable electric boundary conditions. The electric permittivity of air in the crack was considered. Through the Fourier transform, the problems can be solved with the help of two pairs of dual integral equations, in which the unknown variables were the jumps of the displacements across the crack surfaces, not the dislocation density functions. To solve the dual integral equations, the jumps of the displacements across the crack surfaces were directly expanded in a series of Jacobi polynomials. Finally, the effects of the distance between two parallel cracks, the distance between two collinear cracks and the electric boundary conditions on the stress and the electric intensity factors in the piezoelectric materials are analyzed. These results can be used for the strength evaluation of the piezoelectric materials with multi-cracks. The crack shielding effect is also present in the piezoelectric materials.     

Isolated crack in three-dimensional piezoelectric solid. Part II: Stress intensity factors for circular crack

This is part II of the work concerned with finding the stress intensity factors for a circular crack in a solid with piezoelectric behavior. The method of solution involves reducing the problem to a system of hypersingular integral equations by application of the unit concentrated displacement discontinuity and the unit concentrated electric potential discontinuity derived in part I [1]. The near crack border elastic displacement, electric potential, stress and electric displacement are obtained. Stress and electric displacement intensity factors can be expressed in terms of the displacement and the potential discontinuity on the crack surface. Analogy is established between the boundary integral equations for arbitrary shaped cracks in a piezoelectric and elastic medium such that once the stress intensity factors in the piezoelectric medium can be determined directly from that of the elastic medium. Results for the penny-shaped crack are obtained as an example.     

2D dynamic analysis of cracks and interface cracks in piezoelectric composites using the SBFEM

The dynamic stress and electric displacement intensity factors of impermeable cracks in homogeneous piezoelectric materials and interface cracks in piezoelectric bimaterials are evaluated by extending the scaled boundary finite element method (SBFEM). In this method, a piezoelectric plate is divided into polygons. Each polygon is treated as a scaled boundary finite element subdomain. Only the boundaries of the subdomains need to be discretized with line elements. The dynamic properties of a subdomain are represented by the high order stiffness and mass matrices obtained from a continued fraction solution, which is able to represent the high frequency response with only 3–4 terms per wavelength. The semi-analytical solutions model singular stress and electric displacement fields in the vicinity of crack tips accurately and efficiently. The dynamic stress and electric displacement intensity factors are evaluated directly from the scaled boundary finite element solutions. No asymptotic solution, local mesh refinement or other special treatments around a crack tip are required. Numerical examples are presented to verify the proposed technique with the analytical solutions and the results from the literature. The present results highlight the accuracy, simplicity and efficiency of the proposed technique.     

Effects of electric field on crack growth for a penny-shaped dielectric crack in a piezoelectric layer

The assumptions of impermeable and permeable cracks give rise to significant errors in determining electro-elastic behavior of a cracked piezoelectric material. The former simply imposes that the permittivity or electric displacement of the crack interior vanishes, and the latter neglects also the effects of the dielectric of an opening crack interior. Considering the presence of the dielectric of an opening crack interior and the permeability of the crack surfaces for electric field, this paper analyzes electro-elastic behavior induced by a penny-shaped dielectric crack in a piezoelectric ceramic layer. In the cases of prescribed displacement or prescribed stress at the layer surfaces, the Hankel transform technique is employed to reduce the problem to Fredholm integral equations with a parameter dependent nonlinearly on the unknown functions. For an infinite piezoelectric space, a closed-form solution can be derived explicitly, while for a piezoelectric layer, an iterative technique is suggested to solve the resulting nonlinear equations. Field intensity factors are obtained in terms of the solution of the equations. Numerical results of the crack opening displacement intensity factors are presented for a cracked PZT-5H layer and the effect of applied electric field on crack growth are examined for both cases. The results indicate that the fracture toughness of a piezoelectric ceramic is affected by the direction of applied electric fields, dependent on the elastic boundary conditions.