Transient analysis of a piezoelectric strip with a permeable crack under anti-plane impact loads

The problem of a through permeable crack situated in the mid-plane of a piezoelectric strip is considered under anti-plane impact loads for two cases. The first is that the strip boundaries are free of stresses and of electric displacements, and the second is that the strip boundaries are clamped rigid electrodes. The method adopted is to reduce the mixed initial-boundary value problem, by using integral transform techniques, to dual integral equations, which are further transformed into a Fredholm integral equation of the second kind by introducing an auxiliary function. The dynamic stress intensity factor and energy release rate in the Laplace transform domain are obtained in explicit form in terms of the auxiliary function. Some numerical results for the dynamic stress intensity factor are presented graphically in the physical space by using numerical techniques for solving the resulting Fredholm integral equation and inverting Laplace transform.     

Transient analysis of dynamic crack propagation in piezoelectric materials

Abstract

In this paper, the transient analysis of semi‐infinite propagating cracks in piezoelectric materials subjected to dynamic anti‐plane concentrated body force is investigated. The crack surface is assumed to be covered with an infinitesimally thin, perfectly conducting electrode that is grounded. In analyzing this problem, it has characteristic lengths and a direct attempt towards solving this problem by transform and Wiener‐Hopf techniques (Noble, 1958) is not applicable. In order to solve this problem, a new fundamental solution for propagating cracks in piezoelectric materials is first established and the transient response of the propagating crack is obtained by superposition of the fundamental solution in the Laplace transform domain. The fundamental solution to be used is the responses of applying exponentially distributed traction in the Laplace transform domain on the propagating crack surface. Taking into account the quasi‐static approximation, exact analytical transient solutions for the dynamic stress intensity factor and the dynamic electric displacement intensity factor are obtained by using the Cagniard‐de Hoop method (Cagnard, 1939; de Hoop, 1960) of Laplace inversion and are expressed in explicit forms. Numerical calculations of dynamic intensity factors are evaluated and the results are discussed in detail. The transient solutions for stationary cracks have been shown to approach the corresponding static values after the shear wave of the piezoelectric material has passed the crack tip.     

压电复合材料中III型唇形裂纹问题的解析解

采用复变函数方法和保角映射技术,研究了压电复合材料中含唇形裂纹的无限大体远场受反平面机械载荷和面内电载荷作用下的反平面问题,利用复变函数中的留数定理和Cauchy积分公式,分别获得了电不可通和电可通两种边界条件下裂纹尖端场强度因子和机械应变能释放率的解析表达式。当唇形裂纹的高度趋于零时,可得到无限大压电复合材料中Griffith裂纹的解析解。若不考虑电场作用,所得解退化为经典材料的已知结果。数值算例显示了裂纹的几何尺寸和机电载荷对机械应变能释放率的影响规律。结果表明: 唇形裂纹高度的增加会阻碍裂纹的扩展;机械载荷总是促进裂纹的扩展;电载荷对裂纹扩展的影响与裂纹面电边界条件有关。     

Transient response of a center crack in a functionally graded piezoelectric strip under electromechanical impact

The dynamic fracture problem for a functionally graded piezoelectric strip containing a center crack parallel to the free boundaries is considered in this study. It is assumed that the electroelastic properties of the medium vary continuously in the thickness direction, and that the strip is under in-plane mechanical and electric impact loadings. Integral transform techniques and dislocation density functions are employed to reduce the problem to the solutions of a system of singular integral equations. The dynamic stress and electric displacement intensity factors versus time are presented for various values of dimensionless parameters representing the crack size, the material nonhomogeneity and the loading combination.     

A wedge crack in an anisotropic material under antiplane shear

A crack emanating from the apex of an infinite wedge in an anisotropic material under antiplane shear is investigated. An isotropic wedge crack subjected to concentrated forces is first solved by using the conformal mapping technique. The solution of an anisotropic wedge crack is obtained from that of the transformed isotropic wedge crack based on a linear transformation method. Expressions for the stress intensity factor for the anisotropic wedge crack with both concentrated and distributed loads are derived. The stress intensity factors are numerically calculated for generally orthotropic wedge cracks with various crack and wedge angles as well as anisotropic parameters.     

Dynamic stress intensity factor due to concentrated loads on a propagating semi-infinite crack in orthotropic materials

The elastodynamic response of an infinite orthotropic material with a semi-infinite crack propagating at constant speed under the action of concentrated loads on the crack faces is examined. Solution for the stress intensity factor history around the crack tip is found for the loading modes I and II. Laplace and Fourier transforms along with the Wiener-Hopf technique are employed to solve the equations of motion. The asymptotic expression for the stress near the crack tip is analyzed which lead to a closed-form solution of the dynamic stress intensity factor. It is found that the stress intensity factor for the propagating crack is proportional to the stress intensity factor for a stationary crack by a factor similar to the universal function k(v) from the isotropic case. Results are presented for orthotropic materials as well as for the isotropic case.     

A line dislocation interacting with a semi-infinite crack in piezoelectric solid

Closed-form analytical solutions are presented for the physical problem of a semi-infinite crack interacting with a line dislocation under the loading of a line force and a line charge in two-dimensional infinite anisotropic piezoelectric medium. The crack can be a conventional Griffith crack or an anti-crack (a rigid line inhomogeneity). Using the extended Stroh formalism and perturbation technique, the explicit expressions of the field intensity factors and the image force on the dislocation are computed as functions of dislocation location and material constants. The results are discussed and compared with those from special cases existed in the literature. The analytical solutions obtained can be applied to studying interacting cracks and crack branching problems in piezoelectric solids.     

The mode I crack problem for layered piezoelectric plates

The plane strain singular stress problem for piezoelectric composite plates having a central crack is considered. For the case of the crack which is normal to and ends at the interface between the piezoelectric plate and the elastic layer, the order of stress singularity around the tip of the crack is obtained. The Fourier transform technique is used to formulate the problem in terms of a singular integral equation. The singular integral equation is solved by using the Gaus–Jacobi integration formula. Numerical calculations are carried out, and the main results presented are the variation of the stress intensity factor as functions of the geometric parameters, the piezoelectric material properties and the electrical boundary conditions of the layered composites.     

Scattering of the harmonic anti-plane shear waves by a crack in functionally graded piezoelectric materials

In the present paper, the dynamic behavior of a Griffth crack in the functionally graded piezoelectric material (FGPM) is investigated. It is assumed that the elastic stiffness, piezoelectric constant, dielectric permittivity and mass density of the FGPM vary continuously as an exponential function, and that FGPM is under the anti-plane mechanical loading and in-plane electrical loading. By using the Fourier transform and defining the jumps of displacement and electric potential components across the crack surface as the unknown functions, two pairs of dual integral equations are derived. To solve the dual integral equations, the jumps of the displacement and electric potential components across the crack surface are expanded in a series of Jacobi polynomial. Numerical examples are provided to show the effects of material properties on the stress and the electric displacement intensity factors.     

The interface crack problem for bonded piezoelectric and orthotropic layers under antiplane shear loading

The primary objective of this paper is to study the influence of the electroelastic interactions on the stress intensity factor in bonded layers of piezoelectric and orthotropic materials containing a crack along the interface under antiplane shear. Attention is given to a two-layer hybrid laminate formed by adding a layer of piezoelectric ceramic to a unidirectional graphite/epoxy composite or an aluminum layer. Electric displacement or electric field is prescribed on the surfaces of the piezoelectric layer. The problem is formulated in terms of a singular integral equation which is solved by using a relatively simple and efficient technique. A number of examples are given for various material combinations. The results show that the effect of the electroelastic interactions on the stress intensity factor and the energy release rate can be highly significant. This revised version was published online in July 2006 with corrections to the Cover Date.     

Transient thermal stress analysis of an edge crack in a functionally graded material

An edge crack in a strip of a functionally graded material (FGM) is studied under transient thermal loading conditions. The FGM is assumed having constant Young's modulus and Poisson's ratio, but the thermal properties of the material vary along the thickness direction of the strip. Thus the material is elastically homogeneous but thermally nonhomogeneous. This kind of FGMs include some ceramic/ceramic FGMs such as TiC/SiC, MoSi₂/Al₂O₃ and MoSi₂/SiC, and also some ceramic/metal FGMs such as zirconia/nickel and zirconia/steel. A multi-layered material model is used to solve the temperature field. By using the Laplace transform and an asymptotic analysis, an analytical first order temperature solution for short times is obtained. Thermal stress intensity factors (TSIFs) are calculated for a TiC/SiC FGM with various volume fraction profiles of the constituent materials. It is found that the TSIF could be reduced if the thermally shocked cracked edge of the FGM strip is pure TiC, whereas the TSIF is increased if the thermally shocked edge is pure SiC.     

含垂直极化方向穿透的共线半无限裂纹的狭长压电体的平面问题

通过构造新的广义保角映射, 利用广义复变函数方法研究了在电不可渗透边界条件下含有沿垂直于极化方向穿透的共线半无限裂纹的狭长压电体的平面问题, 给出了裂纹尖端处的各场强度因子的解析解。此外, 当狭长体高度和裂纹尺寸按一定的趋势变化时, 还可以得到压电复合材料另外几种新构型的平面问题的解析解。通过数值算例分析得到了两裂纹之间的距离、 狭长体的高度及裂纹面上的受载长度对场强度因子的影响规律。     

Transient response of a cracked magnetoelectric material under the action of in-plane sudden impacts

Dynamic analysis of a crack embedded in a magnetoelectric material is made when subjected to in-plane mechanical, electric and magnetic impacts. The Laplace and Fourier transforms are applied to reduce the associated initial- and mixed-boundary value problem to dual integral equations, and then to singular integral equations with Cauchy kernel. By numerically solving the resulting equation, the dynamic field intensity factors as well as CODs, and energy release rates near the crack tip are evaluated and presented graphically. The effects of applied magnetic and electric impacts on crack growth are discussed. Obtained results show that, different from the static results, applied magnetic and electric impacts can strongly affect dynamic stress intensity factors.     

Transient thermoelectroelastic response of a functionally graded piezoelectric strip with a penny-shaped crack

Transient response of a penny-shaped crack in a plate of a functionally graded piezoelectric material (FGPM) is studied under thermal shock loading conditions. It is assumed that the thermoelectroelastic properties of the strip vary continuously along the thickness of the strip, and that the crack faces are completely insulated. By using both the Laplace and Hankel transforms, the thermal and electromechanical problems are reduced to a singular integral equation and a system of singular integral equations which are solved numerically. The intensity factors vs. time for various crack size, crack position and material nonhomogeneity are obtained.     

引入减基法的压电层合板瞬态响应分析

摘要：用减基法（RBM）结合有限元法、傅里叶变换和Newmark直接积分法,研究了压电层合板在机电耦合载荷下的瞬态响应。用层单元将层合板沿厚度方向进行离散,得到时间域内的运动方程,通过傅里叶变换得到波数域内的控制方程。应用Newmark直接积分法求解波数域内的位移和电势,并在Newmark法求解过程嵌入减基法,构造减基空间,把结构的等效刚度矩阵、质量矩阵和载荷列向量映射到减基空间降阶,得到减缩的Newmark增量式,从而快速求解得到原结构波数域响应,通过傅里叶逆变换得到时域内的响应。以PZT-5A/0°PVDF铺层两相材料复合压电层合板为算例,分析了机电耦合线载荷激励下,位移场和电势场的瞬态响应情况。计算结果表明,求解过程引入减基法能更快得到结构的瞬态响应,并保证了精度。     

Dynamic stress intensity factors of two 3D rectangular permeable cracks in a transversely isotropic piezoelectric material under a time‐harmonic elastic P‐wave

The dynamic behavior of two 3D rectangular permeable cracks in a transversely isotropic piezoelectric material is investigated under an incident harmonic stress wave by using the generalized Almansi's theorem and the Schmidt method. The problem is formulated through double Fourier transform into three pairs of dual integral equations with the displacement jumps across the crack surfaces as the unknown variables. To solve the dual integral equations, the displacement jumps across the crack surfaces are directly expanded as a series of Jacobi polynomials. Finally, the relations among the dynamic stress field and the dynamic electric displacement filed near the crack edges are obtained, and the effects of the shape of the rectangular crack, the characteristics of the harmonic wave, and the distance between two rectangular cracks on the stress and the electric intensity factors in a piezoelectric composite material are analyzed. Copyright © 2015 John Wiley & Sons, Ltd.     

Transient response of a surface crack subjected to dynamic anti-plane concentrated loadings

In this study, the transient response of a surface crack in an elastic solid subjected to dynamic anti-plane concentrated loadings is investigated. The angles of the surface crack and the half-plane are 60° and 90°. In analyzing this problem, an infinite number of diffracted and reflected waves generated by the crack tip and edge boundaries must be taken into account and it will make the analysis extremely difficult. The solutions are determined by superposition of the proposed fundamental solution in the Laplace transform domain and by using the method of image. The fundamental solution to be used is the problem for applying exponentially distributed traction on the crack faces. The exact transient solutions of dynamic stress intensity factor are obtained and expressed in formulations of series form. The solutions are valid for an infinite length of time and have accounted for the contribution of an infinite number of diffracted waves. The explicit value of the dynamic overshot for the perpendicular surface crack is obtained from the analysis. Numerical results are evaluated which indicate that the dynamic stress intensity factors will oscillate near the correspondent static values after the first three or six waves have passed the crack tip.     

Investigation of anti-plane shear behavior of a Griffith crack in a piezoelectric material by using the non-local theory

In this paper, the behavior of a Griffith crack in a piezoelectric material under anti-plane shear loading is investigated by using the non-local theory for impermeable crack surface conditions. By using the Fourier transform, the problem can be solved with two pairs of dual integral equations. These equations are solved using Schmidt method. Numerical examples are provided. Contrary to the previous results, it is found that no stress and electric displacement singularity is present at the crack tip.     

On contact zone models for an electrically impermeable interface crack in a piezoelectric bimaterial

An electrically impermeable interface crack between two semi-infinite piezoelectric planes under remote mechanical tension-shear and electrical loading is studied. Assuming the stresses, strains and displacements are independent on the coordinate x ₂ the expressions for the elastic displacement and potential jumps as well as for the stresses and electrical displacement along the interface via a sectionally holomorphic vector function are found. Introducing an artificial contact zone at the right crack tip and assuming the materials possess the symmetry class 6 mm the problem is reduced for a wide range of bimaterial compounds to a combination of combined Dirichlet–Riemann and Hilbert boundary value problems which are solved analytically. From these solutions clear analytical expressions for characteristic mechanical and electrical parameters are derived. As particular cases of the above mentioned solution the classical (oscillating) and contact zone solutions are obtained. Further, a comparison with an associated solution for an electrically permeable crack has been performed. The fracture mechanical parameters for all models via the remote loads are found analytically and important relationships between these parameters are obtained. Due to these relationships an important algorithm of a numerical method applicable for the investigation of an interface crack in a finite sized piezoelectric bimaterial is suggested.     

Multiple interfacial cracks in dissimilar piezoelectric layers under time harmonic loadings

In this study, the dislocation‐based model is developed to study the interaction between time‐harmonic elastic waves and multiple interface cracks in 2 bonded dissimilar piezoelectric layers. In this model, cracks are represented by a distribution of so‐called electro‐elastic dislocations whose density is to be determined by satisfying the boundary conditions. Using the Fourier transform, this formulation leads to 2 singular integral equations, which can be solved numerically for the densities of electro‐elastic dislocations on a crack surface. The formulation is used to determine dynamic field intensity factors for multiple interface cracks without limitation of number of cracks. The dynamic field intensity factors are then calculated for both permeable and impermeable crack, and finally, numerical results are presented to illustrate the variation of these quantities with the electromechanical coupling, crack spacing, and the frequency of loading.