Closed form solution for two unequal collinear semi-permeable straight cracks in a piezoelectric media

The problem of two unequal collinear straight cracks weakening a poled transversely isotropic piezoelectric ceramic is addressed under semi-permeable electric boundary conditions on the crack faces. The plate has been subjected to combined in-plane normal(to the faces of the cracks) mechanical and electric loads. Problem is formulated employing Stroh formalism and solved using complex variable technique. The elastic field, electric field and energy release rate are obtained in closed analytic form. A case study is presented for poled PZT-5H cracked plate to study the effect of prescribed mechanical load, electric load, inter-crack distance and crack lengths on crack arrest parameters stress intensity factor (SIF), electric displacement intensity factor (EDIF) and mechanical and total energy release rates (ERR). Moreover a comparative study is done of impermeable and semi-permeable crack face boundary conditions on SIF, EDIF and ERR, and results obtained is presented graphically. It is observed that the effect of dielectric medium in the crack gap cannot be ignored.     

A numerical analysis of cracks emanating from an elliptical hole in a 2-D elasticity plate

This paper is concerned with the stress intensity factors (SIFs) of cracks emanating from an elliptical hole in an infinite or a finite plate under biaxial loads by using a boundary element method, which consists of the non-singular displacement discontinuity element presented by Crouch and Starfield and the crack-tip displacement discontinuity elements due to the author. In the boundary element implementation the left or the right crack-tip element is placed locally at the corresponding left or right crack tip on top of the constant displacement discontinuity elements that cover the entire crack surface and the other boundaries. A few numerical examples are included to show that the present approach is very efficient and accurate for the calculating the SIFs of crack problems in an infinite or a finite plate. The present numerical results of cracks emanating from an elliptical hole under biaxial loads can reveal the effect of the elliptical aspect ratio and the transverse load on the SIFs.     

梯度材料中矩形裂纹的对偶边界元方法分析

采用对偶边界元方法分析了梯度材料中的矩形裂纹. 该方法基于层状材料基本解,以非裂纹边界的位移和面力以及裂纹面的间断位移作为未知量. 位移边界积分方程的源点配置在非裂纹边界上,面力边界积分方程的源点配置在裂纹面上. 发展了边界积分方程中不同类型奇异积分的数值方法. 借助层状材料基本解,采用分层方法逼近梯度材料夹层沿厚度方向力学参数的变化. 与均匀介质中矩形裂纹的数值解对比,建议方法可以获得高精度的计算结果. 最后,分析了梯度材料中均匀张应力作用下矩形裂纹的应力强度因子,讨论了梯度材料非均匀参数、夹层厚度和裂纹与夹层之间相对位置对应力强度因子的影响.      

Complex variable approach in studying modified polarization saturation model in two-dimensional semipermeable piezoelectric media

A modified polarization saturation model is proposed and addressed mathematically using a complex variable approach in two-dimensional(2 D) semipermeable piezoelectric media. In this model, an existing polarization saturation(PS) model in 2D piezoelectric media is modified by considering a linearly varying saturated normal electric displacement load in place of a constant normal electric displacement load, applied on a saturated electric zone. A centre cracked infinite 2D piezoelectric domain subject to an arbitrary poling direction and in-plane electromechanical loadings is considered for the analytical and numerical studies. Here, the problem is mathematically modeled as a non-homogeneous Riemann-Hilbert problem in terms of unknown complex potential functions representing electric displacement and stress components. Having solved the Hilbert problem, the solutions to the saturated zone length, the crack opening displacement(COD), the crack opening potential(COP), and the local stress intensity factors(SIFs) are obtained in explicit forms. A numerical study is also presented for the proposed modified model, showing the effects of the saturation condition on the applied electrical loading, the saturation zone length, and the COP. The results of fracture parameters obtained from the proposed model are compared with the existing PS model subject to electrical loading, crack face conditions, and polarization angles.     

Stress state of a piezoelectric ceramic body with a plane crack under antisymmetric loads

The static equilibrium of an electroelastic transversely isotropic space with a plane crack under antisymmetric mechanical loads is studied. The crack is located in the plane of isotropy. Relationships are established between the stress intensity factors (SIFs) for an infinite piezoceramic body and the SIFs for a purely elastic body with a crack of the same form under the same loads. This makes it possible to find the SIFs for an electroelastic body without the need to solve specific electroelasitc problems. As an example, the SIFs are determined for a piezoelastic body with penny-shaped and elliptic cracks under shear __________ Translated from Prikladnaya Mekhanika, Vol. 42, No. 2, pp. 32–42, February 2006.     

Limited permeable crack in an interlayer between piezoelectric materials with different zones of electrical saturation and mechanical yielding

A plane problem for two identical piezoelectric semi-infinite spaces adhered by means of a thin isotropic interlayer is considered. It is assumed that a crack of a limited electric permeability occurs in the interlayer parallel to its faces. Combined electromechanical loading is prescribed at infinity. It is assumed that the interlayer is softer than the adherent materials. To avoid the singularities, which are typical for the Griffith crack model, two distinct zones – a zone of mechanical yielding and a zone of electrical saturation – of unknown lengths are introduced as crack continuations. These lengths can be essentially different, with the zone of mechanical yielding significantly longer or shorter than the zone of electrical saturation. Assuming that the interlayer thickness tends to zero, a constant normal stress is prescribed in the zone of mechanical yielding and a saturated electrical displacement is prescribed in the zone of electrical saturation. Outside of these zones, the semi-infinite spaces are assumed to be perfectly bonded. This formulation results in a linear fracture mechanics problem with unknown pre-fracture zone lengths. The problem, formulated mathematically by a system of two equations of linear relationship, is solved exactly. The unknown yield and saturated zones lengths are found from the conditions of finiteness of stress and electrical displacement at the ends of these zones for the both cases when the electrical saturated zone is longer and shorter than the zone of mechanical yielding. It is shown that the same equation as for the Griffith crack model can be used for the determination of the electrical displacement in the crack region. The main results of the paper are obtained in the form of simple analytical equations which are convenient for engineering applications. Some numerical illustrations in graphical and tabular form show dependencies of the pre-fracture zone lengths, the energy release rate, the mechanical displacement and electrical potential jumps on the electromechanical loading and the electrical permeability of the crack medium.     

SELF-SIMILAR CRACK EXPANSION METHOD FOR THE STRESS INTENSITY FACTOR ANALYSIS

The Self-Similar Crack Expansion (SSCE) method is proposed to evaluate stress intensi-ty factors at crack tips, whereby stress intensity factors of a crack can be determined by the crackopening displacement over the crack, not just by the local displacement around the crack tip. The crackexpansion rate is estimated by taking advantage of the crack self-similarity. Therefore, the accuracy ofthe calculation is improved. The singular integrals on crack tip elements are also analyzed and are pre-cisely evaluated in terms of a special integral analysis. Combination of these two techniques greatly in-creases the accuracy in estimating the stress distribution around the crack tip. A variety of two-dimen-sional cracks, such as subsurface cracks, edge cracks, and their interactions are calculated in terms ofthe self-similar expansion rate. Solutions are satisfied with errors less than 0.5% as compared with theanalytical solutions. Based on the calculations of the crack interactions, a theory for crack interactionsis proposed such that for a group of aligned cracks the summation of the square of SIFs at the right tipsof cracks is always equal to that at the left tips of cracks. This theory was proved by the mehtod ofSelf-Similar Crack Expansion in this paper.     

Application of hypersingular integral equation method to three-dimensional crack in electromagnetothermoelastic multiphase composites

A three-dimensional crack problem in electromagnetothermoelastic multiphase composites (EMTE-MCs) under extended loads is investigated in this paper. Using Green’s functions, the extended general displacement solutions are obtained by the boundary element method. This crack problem is reduced to solving a set of hypersingular integral equations coupled with boundary integral equations, in which the unknown functions are the extended displacement discontinuities. Then, the behavior of the extended displacement discontinuities around the crack front terminating at the interface is analyzed by the main-part analysis method of hypersingular integral equations. Analytical solutions for the extended singular stresses, the extended stress intensity factors (SIFs) and the extended energy release rate near the crack front in EMTE-MCs are provided. Also, a numerical method of the hypersingular integral equations for a rectangular crack subjected to extended loads is put forward with the extended displacement discontinuities approximated by the product of basic density functions and polynomials. In addition, distributions of extended SIFs varying with the shape of the crack are presented. The results show that the present method accurately yields smooth variations of extended SIFs along the crack front.     

An anti-plane propagating crack in a functionally graded piezoelectric strip bonded to a homogeneous piezoelectric strip

Summary A finite crack propagating at constant speed in a functionally graded piezoelectric strip (FGPS) bonded to a homogeneous piezoelectric strip is considered. It is assumed that the electroelastic material properties of the FGPS vary exponentially across the thickness of the strip, and that the bimaterial strip is under combined anti-plane mechanical shear and in-plane electrical loads. The analysis is conducted for the electrically unified crack boundary condition, which includes both the traditional permeable and the impermeable ones. By using the Fourier transform, the problem is reduced to the solution of Fredholm integral equations of the second kind. Numerical results for the stress intensity factor and the crack sliding displacement are presented to show the influences of the crack propagation speed, electric loads, FGPS gradation, crack length, electromechanical coupling coefficient, properties of the bonded homogeneous piezoelectric strip and crack location.     

Energy density and release rate study of an elliptic-arc interface crack in piezoelectric solid under anti-plane shear

The anti-plane problem of an elliptical inhomogeneity with an interfacial crack in piezoelectric materials is investigated. The system is subjected to arbitrary singularity loads (point charge and anti-plane concentrated force) and remote anti-plane mechanical and in-plane electrical loads. Using the complex variable method, the explicit series form solutions for the complex potentials in the matrix and the inclusion regions are derived. The electroelastic field intensity factors, the corresponding energy release rates and the generalized strain energy density at the cracks tips are then provided. The influence of the aspect ratio of the ellipse, the crack geometry and the electromechanical coupling coefficient on the energy release rate and the strain energy density is discussed and shown in graphs. The results indicate that the energy release rate increases with increment of the aspect ratio of the ellipse and the influence of electromechanical coupling coefficient on the energy release rate is significant. The strain energy density decreases with increment of the aspect radio of the ellipse and it is always positive for the cases discussed. The energy release rate, however, can be negative when both mechanical and fields are applied.     

Effect of electric displacement saturation on the hysteretic behavior of ferroelectric ceramics and the initiation and propagation of cracks in piezoelectric ceramics

This paper presents a computational investigation of a proposed simplified account for electric displacement saturation on the hysteretic behavior of initially unpoled ferroelectric ceramics as well as on the initiation and propagation of cracks in poled ferroelectric ceramics within the linear regime of piezoelectricity. For the latter case, experimental observations suggest an odd dependency of the onset of crack initiation in these brittle materials on the orientation of the applied electric field with respect to their poling direction which contradicts theoretical results which propose an even dependency of the energy release rate on the applied electric field within the framework of anisotropic linear piezoelectricity. Electric non-linearities arising at regions of inhomogeneities such as inclusions or at the crack tip are proposed in the literature to avoid this discrepancy. Electric displacement saturation is one such non-linear effect which is investigated in this work. A simplified account of this effect is proposed based on an exponential saturation model of the identified material parameters which can be related to this non-linearity. Its advantage over the superposition of a complex function onto the singular solution of a crack within the framework of linear piezoelectricity lies in the straightforward extension of the proposed approach to problems where no analytical solutions exist. This is outlined based on its incorporation into a rate-dependent ferroelectric model accounting for polarization switching as well as based on its incorporation into a finite element framework capable of simulating the initiation and propagation of cracks in piezoelectric ceramics through strong discontinuities in the displacement field and the electric potential. It is shown that besides the determination of the crack initiation onset also the crack propagation direction is influenced by the appearance of saturation zones arising at the crack tip normal to the polarization direction. The numerically obtained crack paths are found to be close to the experimentally reported results.     

The weight function for cracks in piezoelectrics

The weight function in fracture mechanics is the stress intensity factor at the tip of a crack in an elastic material due to a point load at an arbitrary location in the body containing the crack. For a piezoelectric material, this definition is extended to include the effect of point charges and the presence of an electric displacement intensity factor at the tip of the crack. Thus, the weight function permits the calculation of the crack tip intensity factors for an arbitrary distribution of applied loads and imposed electric charges. In this paper, the weight function for calculating the stress and electric displacement intensity factors for cracks in piezoelectric materials is formulated from Maxwell relationships among the energy release rate, the physical displacements and the electric potential as dependent variables and the applied loads and electric charges as independent variables. These Maxwell relationships arise as a result of an electric enthalpy for the body that can be formulated in terms of the applied loads and imposed electric charges. An electric enthalpy for a body containing an electrically impermeable crack can then be stated that accounts for the presence of loads and charges for a problem that has been solved previously plus the loads and charges associated with an unsolved problem for which the stress and electric displacement intensity factors are to be found. Differentiation of the electric enthalpy twice with respect to the applied loads (or imposed charges) and with respect to the crack length gives rise to Maxwell relationships for the derivative of the crack tip energy release rate with respect to the applied loads (or imposed charges) of the unsolved problem equal to the derivative of the physical displacements (or the electric potential) of the solved problem with respect to the crack length. The Irwin relationship for the crack tip energy release rate in terms of the crack tip intensity factors then allows the intensity factors for the unsolved problem to be formulated, thereby giving the desired weight function. The results are used to derive the weight function for an electrically impermeable Griffith crack in an infinite piezoelectric body, thereby giving the stress intensity factors and the electric displacement intensity factor due to a point load and a point charge anywhere in an infinite piezoelectric body. The use of the weight function to compute the electric displacement factor for an electrically permeable crack is then presented. Explicit results based on a previous analysis are given for a Griffith crack in an infinite body of PZT-5H poled orthogonally to the crack surfaces.     

Transient response of a piezoelectric ceramic with two coplanar cracks under electromechanical impact

The transient response of two coplanar cracks in a piezoelectric ceramic under antiplane mechanical and inplane electric impacting loads is investigated in the present paper. Laplace and Fourier transforms are used to reduce the mixed boundary value problems to Cauchy-type singular integral equations in Laplace transform domain, which are solved numerically. The dynamic stress and electric displacement factors are obtained as the functions of time and geometry parameters. The present study shows that the presence of the dynamic electric field will impede or enhance the propagation of the crack in piezoelectric ceramics at different stages of the dynamic electromechanical load. Moreover, the electromechanical response is greatly affected by the ratio of the space of the cracks and the crack length.     

拉伸载荷作用下单边缺陷-边裂纹应力强度因子研究

利用杂交位移不连续法研究拉伸载荷作用下矩形板中单边缺陷-边裂纹(半圆孔裂纹和半方孔裂纹)问题,给出了这三种平面弹性裂纹问题的应力强度因子的详细数值解。通过半圆孔裂纹问题和半方孔裂纹问题与单边裂纹问题的应力强度因子的比较,发现半圆孔和半方孔对单边裂纹有屏蔽影响。此外,本文的研究结果表明,杂交位移不连续法用于分析平面弹性有限体中复杂裂纹问题的应力强度因子简单且又准确。     

压电体中考虑表面效应时开裂孔洞的断裂特征

研究了反平面机械载荷和面内电载荷作用下压电体中考虑表面效应时孔边双裂纹问题的断裂特征。基于Gurtin-Murdoch表面理论模型，通过构造映射函数，利用复势电弹理论获得了应力场和电位移场的闭合解答。给出了裂纹尖端应力强度因子、电位移场强因子和能量释放率的解析解。讨论了开裂孔洞几何参数和施加力电载荷对电弹场强因子和能量释放率的影响。     

Closed form solution of stress intensity factors for cracks emanating from surface semi-spherical cavity in finite body with energy release rate method

In this paper, a new semi-analytical and semi-engineering method of the closed form solution of stress intensity factors (SIFs) of cracks emanating from a surface semi-spherical cavity in a finite body is derived using the energy release rate theory. A mode of crack opening displacements of a normal slice is established, and the normal slice relevant functions are introduced. The proposed method is both effective and accurate for the problem of three-dimensional cracks emanating from a surface cavity. A series of useful results of SIFs are obtained.     

压电效应下一维六方准晶双材料孔边裂纹的反平面问题研究

利用复变函数方法和保角变换技术研究了压电效应下一维六方准晶双材料中圆孔边单裂纹的反平面问题．考虑电不可渗透型边界条件,运用保角变换和Stroh公式得到了弹性体受远场剪切力和面内电载荷作用下裂纹尖端应力强度因子和能量释放率的解析解. 数值算例分析了几何参数、远场受力、电位移载荷对能量释放率的影响．结果表明：裂纹长度、耦合系数和远场剪切力的减小可以抑制裂纹的扩展．不考虑电场时,声子场应力对能量释放率的影响较小．本文的研究结果可作为研究一维六方压电准晶双材料孔边裂纹问题的理论基础,同时为压电准晶及其复合材料的设计、制备、优化和性能评估提供理论依据．     

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.     

直接计算应力强度因子的扩展有限元法

系统地给出了直接计算应力强度因子的扩展有限元法。该方法以常规有限元法为基础,利用单位分解法思想,通过在近似位移表达式中增加能够反映裂纹面的不连续函数及反映裂尖局部特性的裂尖渐进位移场函数,间接体现裂纹面的存在,从而无需使裂纹面与有限元网格一致,无需在裂尖布置高密度网格,也不需要后处理就可以直接计算出应力强度因子,并且大大简化了前后处理工作。最后通过两个简单算例验证了该方法的精度,分析了影响计算结果的因素,并与采用J积分计算的应力强度因子作了对比,得出了两种方法计算精度相当的结论。     

General treatment of mode III interfacial crack problems in piezoelectric materials

Summary  The anti-plane problem of N collinear interfacial cracks between dissimilar transversely isotropic piezoelectric media, which are subjected to piecewise uniform out-of-plane mechanical loading combined with in-plane electric loading at infinity, and also a line loading at an arbitrary point, is addressed by using the complex function method. In comparison with other relevant works, the present study has two features: one is that the analysis is based on the permeable crack model, i.e. the cracks are considered as permeable thin slits, and, thus, both the normal component of electric displacement and the tangential component of electric field are assumed to be continuous across these slits. The other feature is that explicit closed-form solutions are given not only in piezoelectric media, but also inside cracks when the media are subjected to the most general loading. It is shown that the singularities of electric displacement and electric field in the media are always dependent on that of stress for the general case of loading, and all the singularities of field variables are independent of the applied uniform electric loads at infinity. For the interfacial cracks the electric field is square-root singular at the crack tips and shows jumps across the interface, while the normal component of the electric field is linearly variable inside the crack, but the tangential component is square-root singular. However, for a homogeneous medium with collinear cracks, the electric field is always nonsingular in the medium while the electric displacement exhibits square-root singularity. Moreover, in this case, the electric field inside any crack is equal to a constant when uniform loads are applied at infinity. Received 22 November 1999; accepted for publication 20 July 2000