压电体中裂纹与孤立电偶极子的相互作用

研究压电体裂纹与电偶极子的相互作用,得到问题的闭合解,包括应力-电位移场,裂纹张开位移和电势差,以及裂尖应力强度因子,结果表明,电偶极子的方向对裂纹场的影响可由压电体各向异性方向函数表示;当电偶极子位于裂尖附近时,在原点取在裂尖的局部极坐标系中电偶极子位置的极角对裂尖场的影响可由各向异性方向函数表示,电偶极子引起的裂尖应力强度因子与其距裂尖的距离的-3/2次幂成正比。     

Penny-shaped crack in transversely isotropic piezoelectric materials

Using a method of potential functions introduced successively to integrate the field equations of three-dimensional problems for transversely isotropic piezoelectric materials, we obtain the so-called general solution in which the displacement components and electric potential functions are represented by a singular function satisfying some special partial differential equations of 6th order. In order to analyse the mechanical-electric coupling behaviour of penny-shaped crack for above materials, another form of the general solution is obtained under cylindrical coordinate system by introducing three quasi-harmonic functions into the general equations obtained above. It is shown that both the two forms of the general solutions are complete. Furthermore, the mechanical-electric coupling behaviour of penny-shaped crack in transversely isotropic piezoelectric media is analysed under axisymmetric tensile loading case, and the crack-tip stress field and electric displacement field are obtained. The results show that the stress and the electric displacement components near the crack tip have (r ^−1/2) singularity. The project supported by the Natural Science Foundation of Shaanxi Province, China     

A crack in a piezoelectric layer bonded to a dissimilar elastic layer under transient load

Summary  A piezoelectric layer bonded to the surface of an elastic structure is considered. The piezoelectric and the elastic layers are infinite along the x-axis and have finite thickness in the y-direction. The polarization direction of the piezoelectric material is along the y-axis. By means of the method of singular integral equations, the solution in a Laplace transform plane is demonstrated. Laplace inversion yields the results in the time domain. Numerical values of the crack tip fields under in-plane transient electromechanical loading are obtained. The influence of layers thickness on stress and electric displacement intensity factors is investigated. Received 16 March 2000; accepted for publication 16 August 2000     

Stress intensity factor measurement of cracks using a piezoelectric element

We present a stress intensity factor (SIF) measurement method of cracks using a piezoelectric element and electrostatic voltmeter. In the method, an isotropic piezoelectric element is first adhered near the crack tip. Then, the surface electrodes are attached to the three different positions on the piezoelectric element. The electric potentials of the surface electrodes, which are proportional to the strain sum (ɛ_x+ɛ_y) on the structural member, are measured by an electrostatic voltmeter during load cycling. Mode I and mode II SIFs of the crack are estimated using the relationship between the SIF and (σ_x+σ_y). The applicability of the proposed method is examined through experiments and numerical analysis.     

Arbitrarily oriented crack near interface in piezoelectric bimaterials

Arbitrarily oriented crack near interface in piezoelectric bimaterials is considered. After deriving the fundamental solution for an edge dislocation near the interface, the present problem can be expressed as a system of singular integral equations by modeling the crack as continuously distributed edge dislocations. In the paper, the dislocations are described by a density function defined on the crack line. By solving the singular integral equations numerically, the dislocation density function is determined. Then, the stress intensity factors (SIFs) and the electric displacement intensity factor (EDIF) at the crack tips are evaluated. Subsequently, the influences of the interface on crack tip SIFs, EDIF, and the mechanical strain energy release rate (MSERR) are investigated. The J-integral analysis in piezoelectric bimaterals is also performed. It is found that the path-independent of J₁-integral and the path-dependent of J₂-integral found in no-piezoelectric bimaterials are still valid in piezoelectric bimaterials.     

Damage analysis of tetragonal perovskite structure ceramics implicated by asymptotic field solutions and boundary conditions

Cracking of ceramics with tetragonal perovskite grain structure is known to appear at different sites and scale level. The multiscale character of damage depends on the combined effects of electromechanical coupling, prevailing physical parameters and boundary conditions. These detail features are exhibited by application of the energy density criterion with judicious use of the mode I asymptotic and full field solution in the range of r/a=10⁻⁴ to 10⁻² where r and a are, respectively, the distance to the crack tip and half crack length. Very close to the stationary crack tip, bifurcation is predicted resembling the dislocation emission behavior invoked in the molecular dynamics model. At the macroscopic scale, crack growth is predicted to occur straight ahead with two yield zones to the sides. A multiscale feature of crack tip damage is provided for the first time. Numerical values of the relative distances and bifurcation angles are reported for the PZT-4 ceramic subjected to different electric field to applied stress ratio and boundary conditions that consist of the specification of electric field/mechanical stress, electric displacement/mechanical strain, and mixed conditions. To be emphasized is that the multiscale character of damage in piezoceramics does not appear in general. It occurs only for specific combinations of the external and internal field parameters, elastic/piezoelectric/dielectric constants and specified boundary conditions.     

Electromechanical impact of a crack in a functionally graded piezoelectric medium

The Mode-I transient response of a functionally graded piezoelectric medium is solved for a through crack under the in-plane mechanical and electric impact. Integral transforms and dislocation density functions are employed to reduce the problem to singular integral equations. Numerical results display the effects of the loading combination parameter λ and the material parameter βa on the dynamic stress intensity factor and electric displacement intensity factor. The energy density factor criterion is applied to obtain the maximum of the minimum energy density factor and the direction of crack initiation.     

Crack problems of piezoelectric materials

This paper presents an analysis of crack problems in homogeneous piezoelectrics or on the interfaces between two dissimilar piezoelectric materials based on the continuity of normal electric displacement and electric potential across the crack faces. The explicit analytic solutions are obtained for a single crack in piezoelectrics or on the interfaces of piezoelectric bimaterials. A class of boundary problems involving many cracks is also solved. For homogeneous materials it is found that the normal electric displacementD ₂ induced by the crack is constant along the crack faces which depends only on the applied remote stress field. Within the crack slit, the electric fields induced by the crack are also constant and not affected by the applied electric field. For the bimaterials with realH, the normal electric displacementD ₂ is constant along the crack faces and electric fieldE ₂ has the singularity ahead of the crack tip and a jump across the interface. The project is supported by the National Natural Science Foundation of China(No. 19704100) and the Natural Science Foundation of Chinese Academy of Sciences(No. KJ951-1-201).     

The Evolution of Stress Intensity Factors and the Propagation of Cracks in Elastic Media

When a crack Γ _s propagates in an elastic medium the stress intensity factors evolve with the tip x(s) of Γ _s . In this paper we derive formulae which describe the evolution of these stress intensity factors for a homogeneous isotropic elastic medium under plane strain conditions. Denoting by ψ=ψ(x,s) the stress potential (ψ is biharmonic and has zero traction along the crack Γ _s ) and by κ(s) the curvature of the crack at the tip x(s), we prove that the stress intensity factors A ₁(s), A ₂(s), as functions of s, satisfy:

Dynamic stress intensity factor K III and dynamic crack propagation characteristics of anisotropic materials

Based on the mechanics of anisotropic materials, the dynamic propagation problem of a mode Ⅲ crack in an infinite anisotropic body is investigated. Stress, strain and displacement around the crack tip are expressed as an analytical complex function, which can be represented in power series. Constant coefficients of series are determined by boundary conditions. Expressions of dynamic stress intensity factors for a mode Ⅲ crack are obtained. Components of dynamic stress, dynamic strain and dynamic displacement around the crack tip are derived. Crack propagation characteristics are represented by the mechanical properties of the anisotropic materials, i.e., crack propagation velocity M and the parameter ~. The faster the crack velocity is, the greater the maximums of stress components and dynamic displacement components around the crack tip are. In particular, the parameter α affects stress and dynamic displacement around the crack tip.     

Fracture-mechanical assessment of electrically permeable interface cracks in piezoelectric bimaterials by consideration of various contact zone models

Summary An interface crack with an artificial contact zone at the right-hand side crack tip between two piezoelectric semi-infinite half-planes is considered under remote mixed-mode loading. Assuming the stresses, strains and displacements are independent of the coordinate x ₂, the expression for the displacement jumps and stresses along the interface are found via a sectionally holomorphic vector function. For piezoceramics of the symmetry class 6 mm and for electrically permeable crack faces, the problem is reduced to a combined Dirichlet-Riemann boundary value problem which can be solved analytically. Further, analytical expressions for the stresses, electrical displacements, derivatives of elastic displacement jumps, stress and electrical intensity factors are found at the interface. Real contact zone lengths and the well-known oscillating solution are derived from the obtained solution as well. Analytical relationships between the fracture-mechanical parameters of various models are found, and recommendations are suggested concerning the application of numerical methods to the problem of an interface crack in the discontinuity area of a piezoelectric bimaterial. Received 16 March 1999; accepted for publication 31 May 1999     

功能梯度压电材料反平面裂纹问题

基于三维弹性理论和压电理论，导出了材料系数在横观各向同性平面内梯度分布的压电体的状态方程，进而对材料系数指数函数规律分布的半无限大压电体中的反平面裂纹问题进行了求解，利用Fourier变换给出了半无限大压电体中位移，应力，电势及电位移的解析表达式，并求得了裂纹尖端的应力强度因子和电位移强度因子，分析了不同的非均匀材料系数及几何尺寸对它们的影响。     

Antiplane internal crack normal to the edge of a functionally graded piezoelectric/piezomagnetic half plane

This paper deals with the antiplane magnetoelectroelastic problem of an internal crack normal to the edge of a functionally graded piezoelectric/piezomagnetic half plane. The properties of the material such as elastic modulus, piezoelectric constant, dielectric constant, piezomagnetic coefficient, magnetoelectric coefficient and magnetic permeability are assumed in exponential forms and vary along the crack direction. Fourier transforms are used to reduce the impermeable and permeable crack problems to a system of singular integral equations, which is solved numerically by using the Gauss-Chebyshev integration technique. The stress, electric displacement and magnetic induction intensity factors at the crack tips are determined numerically. The energy density theory is applied to study the effects of nonhomogeneous material parameter β, edge conditions, location of the crack and load ratios on the fracture behavior of the internal crack.     

压电材料平面问题的尖端场和应力强度因子的求解

利用Lekhnitskii理论和Stroh理论的相互联系,把已知的基于Lekhnitskii理论平面应变结果转化为Stroh理论形式的结果,直接获得Stroh公式中A,B的显式表达式,此方法可扩展到平面应力情况,然后导出压电材料平面应变问题的尖端场Williams形式的展开式,采用半权函数法计算有限大压电体平面问题应力和电位移强度因子。对无穷大板含中心裂纹的情况下本文结果和已有结果进行了比较,表明本文方法得到的结果精度可靠。本文方法的最大优点是可以求解有限压电体的应力强度因子,并且需要的单元少,精度高,实用性好。     

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.     

Investigation of the Dynamic Behavior of a Griffith Crack in a Piezoelectric Material Strip Subjected to the Harmonic Elastic Anti-plane Shear Waves by Use of the Non-local Theory

In this paper, the dynamic behavior of a Griffith crack in a piezoelectric material strip subjected to the harmonic anti-plane shear waves is investigated by use of the non-local theory for impermeable crack surface conditions. To overcome the mathematical difficulties, a one-dimensional non-local kernel is used instead of a two-dimensional one for the anti-plane dynamic problem to obtain the stress and the electric displacement near at the crack tip. By means of the Fourier transform, the problem can be solved with the help of two pairs of dual integral equations. These equations are solved using the Schmidt method. Contrary to the classical solution, it is found that no stress and electric displacement singularity is present near the crack tip. The non-local dynamic elastic solutions yield a finite hoop stress near the crack tip, thus allowing for a fracture criterion based on the maximum dynamic stress hypothesis. The finite hoop stress at the crack tip depends on the crack length, the thickness of the strip, the circular frequency of incident wave and the lattice parameter.     

Griffith crack moving in a piezoelectric strip

Summary  The dynamic problem of an impermeable crack of constant length 2a propagating along a piezoelectric ceramic strip is considered under the action of uniform anti-plane shear stress and uniform electric field. The integral transform technique is employed to reduce the mixed-boundary-value problem to a singular integral equation. For the case of a crack moving in the mid-plane, explicit analytic expressions for the electroelastic field and the field intensity factors are obtained, while for an eccentric crack moving along a piezoelectric strip, numerical results are determined via the Lobatto–Chebyshev collocation method for solving a resulting singular integral equation. The results reveal that the electric-displacement intensity factor is independent of the crack velocity, while other field intensity factors depend on the crack velocity when referred to the moving coordinate system. If the crack velocity vanishes, the present results reduce to those for a stationary crack in a piezoelectric strip. In contrast to the results for a stationary crack, applied stress gives rise to a singular electric field and applied electric field results in a singular stress for a moving crack in a piezoelectric strip. Received 14 August 2001; accepted for publication 24 September 2002 The author is indebted to the AAM Reviewers for their helpful suggestions for improving this paper. The work was supported by the National Natural Science Foundation of China under Grant 70272043.     

Investigation of the dynamic behavior of two parallel symmetric cracks in piezoelectric materials use of non-local theory

In this paper, the dynamic behavior of two parallel symmetric cracks in piezoelectric materials under harmonic anti-plane shear waves is investigated by use of the non-local theory for permeable crack surface conditions. To overcome the mathematical difficulties, a one-dimensional non-local kernel is used instead of a two-dimensional one for the problem to obtain the stress occurs near the crack tips. By means of the Fourier transform, the problem can be solved with the help of two pairs of dual integral equations that the unknown variables are the jumps of the displacement along the crack surfaces. These equations are solved using the Schmidt method. Numerical examples are provided. Contrary to the previous results, it is found that no stress and electric displacement singularity is present near the crack tip. The non-local elastic solutions yield a finite hoop stress near the crack tip, thus allowing for a fracture criterion based on the maximum stress hypothesis. The finite hoop stress at the crack tip depends on the crack length, the frequency of the incident wave, the distance between two cracks and the lattice parameter of the materials, respectively. Contrary to the impermeable crack surface condition solution, it is found that the dynamic electric displacement for the permeable crack surface conditions is much smaller than the results for the impermeable crack surface conditions. The results show that the dynamic field will impede or enhance crack propagation in the piezoelectric materials at different stages of the dynamic load.     

功能梯度压电板条中电绝缘型运动裂纹的电弹性场

基于三维弹性理论和压电理论，对材料系数按指数函数规律分布的功能梯度压电板条中的反平面运动裂纹问题进行了求解。利用Fourier积分变换方法将电绝缘型运动裂纹问题化为对偶积分方程，并进一步归结为易于求解的第二类Fredholm积分方程。通过渐近分析，获得了裂纹尖端应力、应变、电位移和电场的解析解，给出了裂纹尖端场各个变量的角分布函数，并求得了裂纹尖端场的强度因子，分析了压电材料物性梯度参数、几何尺寸及裂纹运动速度对它们的影响。结果表明，对于电绝缘型裂纹，功能梯度压电板条中运动裂纹尖端附近的各个场变量都具有-1／2阶的奇异性；当裂纹运动速度增大时，裂纹扩展的方向会偏离裂纹面。     

A method for measuring mode I crack tip constraint under static and dynamic loading conditions

A novel experimental technique for measuring crack tipT-stress, and hence in-plane crack tip constraint, in elastic materials has been developed. The method exploits optimal positioning of stacked strain gage rosette near a mode I crack tip such that the influence of dominant singular strains is negated in order to determineT-stress accurately. The method is demonstrated for quasi-static and low-velocity impact loading conditions and two values of crack length to plate width ratios (a/W). By coupling this new method with the Dally-Sanford single strain gage method for measuring the mode I stress intensity factorK _I , the crack tip biaxiality parameter is also measured experimentally. Complementary small strain, static and dynamic finite element simulations are carried out under plane stress conditions. Time histories ofK _I andT-stress are computed by regression analysis of the displacement and stress fields, respectively. The experimental results are in good agreement with those obtained from numerical simulations. Preliminary data for critical values ofK _I and β for dynamic experiments involving epoxy specimens are reported. Dynamic crack initiation toughness shows an increasing trend as β becomes more negative at higher impact velocities.