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.     

Multiple cracks in thermoelectroelastic bimaterials

The formulation for thermal stress and electric displacement in an infinite thermopiezoelectric plate with an interface and multiple cracks is presented. Using Green's function approach and the principle of superposition, a system of singular integral equations for the unknown temperature discontinuity defined on each crack face is developed and solved numerically. The formulation can then be used to calculate some fracture parameters such as the stress–electric displacement and strain energy density factor. The direction of crack growth for many cracks in thermopiezoelectric bimaterials is predicted by way of the strain energy density theory. Numerical results for stress–electric displacement factors and crack growth direction at a particular crack tip in two crack system of bimaterials are presented to illustrate the application of the proposed formulation.     

On a semi-infinite crack penetrating a piezoelectric circular inhomogeneity with a viscous interface

We investigate a semi-infinite crack penetrating a piezoelectric circular inhomogeneity bonded to an infinite piezoelectric matrix through a linear viscous interface. The tip of the crack is at the center of the circular inhomogeneity. By means of the complex variable and conformal mapping methods, exact closed-form solutions in terms of elementary functions are derived for the following three loading cases: (i) nominal Mode-III stress and electric displacement intensity factors at infinity; (ii) a piezoelectric screw dislocation located in the unbounded matrix; and (iii) a piezoelectric screw dislocation located in the inhomogeneity. The time-dependent electroelastic field in the cracked composite system is obtained. Particularly the time-dependent stress and electric displacement intensity factors at the crack tip, jumps in the displacement and electric potential across the crack surfaces, displacement jump across the viscous interface, and image force acting on the piezoelectric screw dislocation are all derived. It is found that the value of the relaxation (or characteristic) time for this cracked composite system is just twice as that for the same fibrous composite system without crack. Finally, we extend the methods to the more general scenario where a semi-infinite wedge crack is within the inhomogeneity/matrix composite system with a viscous interface.     

Dynamic interaction of multiple arbitrarily oriented planar cracks in a piezoelectric space: A semi-analytic solution

The problem of determining the electro-elastic fields around arbitrarily oriented planar cracks in an infinite piezoelectric space is considered. The cracks which are acted upon by a transient load are either electrically impermeable or permeable. A semi-analytic method based on the theory of exponential Fourier transformation is proposed for solving the problem in the Laplace transform domain. The Laplace transforms of the jumps in the displacements and electric potential across opposite crack faces are determined by solving a system of hypersingular integral equations. Once these displacement and electric potential jumps are obtained, the displacements and electric potential and other physical quantities of interest, such as the crack tip stress and electric displacement intensity factors, can be computed with the help of a suitable algorithm for inverting Laplace transforms. The stress and electric displacement intensity factors are computed for some specific cases of the problem.     

压电、压磁材料球对称问题的通解

研究压电、压磁材料在球坐标系下,不计体力、体电荷和体电流的情况下,由平衡方程、梯度方程、压电和压磁的本构方程导出应力、应变、位移、电位移、电场强度、电位势、磁感强度和磁位势各未知量的通解.考虑不同的边界条件,将其通解应用到应力、电学短路以及位移、电学开路和磁场分布的边界条件中,得到不同边界条件下问题的解.     

A linearized theory of elastic dielectrics which couples electrically to flexural vibrations

Assuming that the free energy depends on the deformation gradient and the spatial electric field, we derive the expressions for the Cauchy stress tensor and the spatial electric displacement from an observer invariant quadratic form of the free energy via the strict definitions of these quantities. Specific forms of the Piola-Kirchhoff stress tensor and the material electric displacement are then deduced and linearized in a particular sense. As an application of the resulting theory, we formulate the problem of an electrically driven disc within the context of the classical bending theory of thin plates. The material of the disc is assumed to have at most the symmetry of a hexagonal system of classC _6v.The resulting coupled differential equations for the axial mechanical displacement of the middle surface and the material electric potential indicate that the problem is not empty. This result is of particular interest in view of the fact that it is generally held that the classical theory of piezoelectricity does not permit such couplings to occur.     

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     

压电薄膜表面贯穿裂纹的有限元分析

基于有限元软件ANSYS数值模拟,计算了激光作用下的压电薄膜表面贯穿裂纹外场应力强度因子和电位移强度因子,并且研究了90°畴变所诱致的畴变增韧行为。首先,求解无裂纹压电薄膜在激光作用下的热-力-电响应,将求得的应力和电位移场反向作用于裂纹面,求解裂纹尖端处的外场应力和电位移强度因子,然后基于小范围畴变理论求解了90°畴变所致的屏蔽应力强度因子。讨论了薄膜表面裂纹的外场应力强度因子、电位移强度因子及屏蔽应力强度因子随激光作用时间和裂纹位置的变化关系,从而预测压电薄膜体系在加热工作状况下的裂纹扩展和断裂行为。     

Dielectric crack problem for a magnetoelectroelastic strip with functionally graded properties

The paper presents a fracture analysis for an electromagnetically dielectric crack in a functionally graded magnetoelectroelastic strip. It is considered that the material properties are varying exponentially along the width direction. Under the assumption of the in-plane magneto-electro-mechanical loadings, the dielectric crack is simulated by using the semi-permeable crack-face boundary conditions. The Fourier transform technique is applied to solve the boundary-value problem and four coupling singular integral equations are determined. A nonlinear system of algebraic equations is further derived and solved numerically to determine the electromagnetic field inside the crack. Then the field intensity factors of stress, electric displacement, and magnetic induction are given. Through the numerical computations, the effects of the material non-homogeneity and the permeability of crack interior on the electric displacement and the magnetic induction at the crack faces are studied. The variations of the intensity factors of stress, electric displacement, and magnetic induction versus the geometry of the crack, the strip width, and the material non-homogeneity are presented in graphics respectively.     

Green's function for thermopiezoelectric plates with holes of various shapes

Summary Thermoelectroelastic problems for holes of various shapes embedded in an infinite matrix are considered in this paper. Based on Lekhnitskii's formalism, the technique of conformal mapping and the exact electric boundary conditions on the hole boundary, the thermoelectroelastic Green's function has been obtained analytically in terms of a complex potential. As an application of the proposed function, the problem of an infinite plate containing a crack and a hole is analysed. A system of singular integral equations for the unknown temperature discontinuity and the discontinuity of elastic displacement and electric potential (EDEP) defined on crack faces is developed and solved numerically. Numerical results for stress and electric displacement (SED) intensity factors of the crack-hole system are presented to illustrate the application of the proposed formulation. Received 7 October 1998; accepted for publication 26 January 1999     

Analysis of a Griffith crack at the interface of two piezoelectric materials under anti-plane loading

The main objective of this work is the contribution to the study of the piezoelectric structures which contain preexisting defect (crack). For that, we consider a Griffith crack located at the interface of two piezoelectric materials in a semi-infinite plane structure. The structure is subjected to an anti-plane shearing combined with an in-plane electric displacement. Using integral Fourier transforms, the equations of piezoelectricity are converted analytically to a system of singular integral equations. The singular integral equations are further reduced to a system of algebraic equations and solved numerically by using Chebyshev polynomials. The stress intensity factor and the electric displacement intensity factor are calculated and used for the determination of the energy release rate which will be taken as fracture criterion. At the end, numerical results are presented for various parameters of the problem; they are also presented for an infinite plane structure.     

Subinterface crack in an anisotropic piezoelectric bimaterial

In this paper, the problem of a subinterface crack in an anisotropic piezoelectric bimaterial is analyzed. A system of singular integral equations is formulated for general anisotropic piezoelectric bimaterial with kernel functions expressed in complex form. For commonly used transversely isotropic piezoelectric materials, the kernel functions are given in real forms. By considering special properties of one of the bimaterial, various real kernel functions for half-plane problems with mechanical traction-free or displacement-fixed boundary conditions combined with different electric boundary conditions are obtained. Investigations of half-plane piezoelectric solids show that, particularly for the mechanical traction-free problem, the evaluations of the mechanical stress intensity factors (electric displacement intensity factor) under mechanical loadings (electric displacement loading) for coupled mechanical and electric problems may be evaluated directly by considering the corresponding decoupled elastic (electric) problem irrespective of what electric boundary condition is applied on the boundary. However, for the piezoelectric bimaterial problem, purely elastic bimaterial analysis or purely electric bimaterial analysis is inadequate for the determination of the generalized stress intensity factors. Instead, both elastic and electric properties of the bimaterial’s constants should be simultaneously taken into account for better accuracy of the generalized stress intensity factors.     

Anti-plane shear crack in a functionally gradient piezoelectric material

The main objective of this paper is to study the singular nature of the crack-tip stress and electric displacement field in a functionally gradient piezoelectric medium having material coefficients with a discontinuous derivative. The problem is considered for the simplest possible loading and geometry, namely, the anti-plane shear stress and electric displacement in-plane of two bonded half spaces in which the crack is parallel to the interface. It is shown that the square-root singularity of the crack-tip stress field and electric displacement field is unaffected by the discontinuity in the derivative of the material coefficients. The problem is solved for the case of a finite crack and extensive results are given for the stress intensity factors, electric displacement intensity factors, and the energy release rate. Project supported by the National Natural Science Foundation of China (No. 10072041), the National Excellent Young Scholar Fund, of China (No. 10125209) and the Teaching and Research Award Program for Outstanding Young Teachers in Higher Education Institutions of MOE, P. R. C..     

Transient response of a cracked piezoelectric strip under arbitrary electro-mechanical impact

By using the well-developed integral transform methodology, the dynamic response of stress and electric displacement around a finite crack in an infinite piezoelectric strip are investigated under arbitrary dynamic anti-plane loads. The dynamic stress intensity factors and electric displacement are obtained analytically. It is shown that the dynamic crack-tip stress and electric field still have a square-root singularity. Numerical computations for the dynamic stress intensity factor show that the electric load has a significant influence on the dynamic response of stress field. The higher the ratio of the crack length to the width of the strip, the higher the peak value of the dynamic stress intensity factor is. On the other hand, the dynamic response of the electric field is determined solely by the applied electric load. The electric field will promote or retard the propagation of the crack depending on the time elapse since the application of the external electro-mechanical loads. The project supported by the National Natural Science Foundation of China and the Post-Doctor Science Foundation of China     

ANALYSIS OF THE DYNAMIC BEHAVIOR OF A GRIFFITH PERMEABLE CRACK IN PIEZOELECTRIC MATERIALS WITH THE NON-LOCAL THEORY

The dynamic behavior ofa Griffith permeable crack under harmonic anti-plane shearwaves in the piezoelectric materials is investigated by use of the non-local theory.To overcome themathematical difficulties,a one-dimensional non-local kernel is used instead of a two-dimensionalone for the anti-plane dynamic problem to obtain the stress and the electric displacement near thecrack tips.By means of Fourier transform,the problem can be solved with a pair of dual integralequations that the unknown variable is the jump of the displacement across the crack surfaces.These equations are solved with the Schmidt method and numerical examples are provided.Con-trary to the previous results,it is found that no stress and electric displacement singularities arepresent at the crack tip.The finite hoop stress and the electric displacement depend on the cracklength,the lattice parameter of the materials and the circle frequency of the incident waves.Thisenables us to employ the maximum stress hypothesis to deal with fracture problems in a naturalway.     

Stress intensity factors for an interface crack with impermeable surfaces between dissimilar electrostrictive materials

The asymptotic problem of a semi-infinite interface crack between dissimilar electrostrictive materials that are subjected to electric loading is numerically analyzed by using the finite element method. Numerical results of electric displacement fields are obtained on the basis of the mathematical equivalence of the mode III problem and an electrostatic problem. The shape and the size of saturation zones are explored as a function of the ratio of the saturated electric displacements of dissimilar electrostrictive materials. In contrast with conventional wisdom, the ratio of the permittivities is shown to exert a negligibly small influence on electric displacement fields. For various combinations of the material properties of dissimilar electrostrictive materials, stress fields and stress intensity factors are systemically calculated by using the numerical results of electric displacement fields. The effects of the electric, elastic, and electrostrictive properties on stress intensity factors are demonstrated.     

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.     

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

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

Numerical analysis for a crack in piezoelectric material under impact

In this paper, a numerical analysis of impact interfacial fracture for a piezoelectric bimaterial is provided. Starting from the basic equilibrium equation, a dynamic electro-mechanical FEM formulation is briefly presented. Then, the path-independent separated dynamic J integral is extended to piezoelectric bimaterials. Based on the relationship of the path-independent dynamic J integral and the stress and electric displacement intensity factors, the component separation method is used to calculate the stress and electric displacement intensity factors for piezoelectric bimaterials in this finite-element analysis. The response curves of the dynamic J integral, the stress and electric displacement intensity factors are obtained for both homogeneous material (PZT-4 and CdSe) and CdSe/PZT-4 bimaterial. The influences of the piezoelectricity and the electro-mechanical coupling factor on these responses are discussed. The effects of an applied electric field are also discussed.     

Dynamic analysis of a penny-shaped dielectric crack in a magnetoelectroelastic solid under impacts

The transient response of a magneto-electro-elastic material with a penny-shaped dielectric crack subjected to in-plane magneto-electro-mechanical impacts is made. To simulate an opening crack with a dielectric interior, the crack-face electromagnetic boundary conditions are supposed to depend on the crack opening displacement and the jumps of electric and magnetic potentials across the crack. Four ideal crack-face electromagnetic boundary conditions involving a combination of electrically permeable or impermeable and magnetically permeable or impermeable assumptions can be reduced. The Laplace and Hankel transform techniques are further utilized to solve the mixed initial-boundary-value problem. Three coupling Fredholm integral equations are obtained and solved by the composite Simpson's rule. Dynamic field intensity factors of stress, electric displacement, magnetic induction, crack opening displacement (COD), electric potential and magnetic potential are given in the Laplace transform domain. By means of a numerical inversion of the Laplace transform, numerical results are calculated to show the variations of the physical parameters of concern versus the normalized time in graphics. The effects of applied electric and magnetic loads on the dynamic intensity factors of stress and COD, and the dynamic energy release rate for a BaTiO₃-CoFe₂O₄ composite with a penny-shaped vacuum crack are discussed in detail.