# Nonuniformly Loaded Stack of Antiplane Shear Cracks in One-Dimensional Piezoelectric Quasicrystals

Authors: G. E. Tupholme
Journal: Advances in Materials Science and Engineering
 (2018)

---

## Abstract

Representations in a closed form are derived, using an extension to the method of dislocation layers, for the phonon and phason stress and electric displacement components in the deformation of one-dimensional piezoelectric quasicrystals by a nonuniformly loaded stack of parallel antiplane shear cracks. Their dependence upon the polar angle in the region close to the tip of a crack is deduced, and the field intensity factors then follow. These exhibit that the phenomenon of crack shielding is dependent upon the relative spacing of the cracks. The analogous analyses, that have not been given previously, involving non-piezoelectric or non-quasicrystalline or simply elastic materials can be straightforwardly considered as special cases. Even when the loading is uniform and the crack is embedded in a purely elastic isotropic solid, no explicit representations have been available before for the components of the field at points other than directly ahead of a crack. Typical numerical results are graphically displayed.


---

## Body

## 1. Introduction

Since the ground-breaking investigation of strip cracks by Griffith [1], almost a century ago, there has been a flourishing technological need for analyzing the increasingly complex behaviour of defects.For elastic isotropic solids, Sneddon and Lowengrub [2] obtained the stress-intensity factor in the vicinity of the tip of a strip crack within a row of cracks under the influence of a constant load. Analogously numerous other crack configurations in various materials have been presented. However, invariably most of them do not give representations in a closed form for the field components at a general point, but only directly ahead of the crack.Here, closed-form expressions are derived for the components of the phonon and phason stresses and electric displacement and field intensity factors created around a stack of parallel shear cracks subjected tononuniform mechanical and electric loads in one-dimensional piezoelectric quasicrystals, using an extension of the method of dislocation arrays. Currently, these are apparently the most comprehensive materials for which it is readily feasible to apply this technique. The corresponding results when the material is simply quasicrystalline or solely piezoelectric or purely isotropic elastic can be deduced as special cases.The discovery of some unusual complex metallic alloys having quasiperiodic structures was first widely announced by Shechtman et al. [3]. They are now known universally as quasicrystals. The continuing development of the performance and design of smart devices that exploit the effects of piezoelectric coupling of quasicrystals is aided by analyzing the influence of cracks within them.Experimentally, quasicrystalline materials have been demonstrated to be extremely brittle and often to contain very high densities of multiple cracks. Thus, it is desirable to have a breadth of basic exact solutions readily available for fundamental defect problems, against which more general numerical and experimental studies of analytically intractable practical crack situations can be compared and verified as particular cases.The underlying equations that govern the linear elastic theory of quasicrystals within the framework of the mechanics of continua are now widely established. Their solutions continue to be extensively developed for studying a wide variety of boundary value problems and especially those related to cracks. A fundamental analysis that extended the traditional theory of cracks in linear elastic solids to quasicrystals was presented by Li et al. [4]. Subsequently, Li [5, 6] formulated a general explicit solution of elastohydrodynamic problems in decagonal quasicrystals which enabled wave characteristics in them to be analyzed. Comprehensive accounts and discussions of, and references to, the literature on quasicrystals are conveniently provided for an interested reader by, for example, Ding et al. [7], Guo et al. [8], Fan et al. [9], and Fan [10, 11]. Yadav [12] gave an interesting summary of the 56 papers that were presented in 2016 at the most recent 13th International Conference on Quasicrystals (ICQ13).Henceforth, Altay and Dökmeci [13] extended the three-dimensional equations of elasticity of quasicrystals to include quasi-static electric effects. This has enabled piezoelectric quasicrystal boundary value problems to be addressed.The underpinning results of Li and Liu [14] on the physical property tensors enabled Wang and Pan [15] to give detailed closed-form representations of the field components of a screw dislocation moving uniformly in one-dimensional hexagonal piezoelectric quasicrystals. Yang et al. [16] obtained and investigated numerically the electric/elastic field of a straight stationary dislocation based upon the generalized Stroh formalism.To aid in providing general static solutions to three-dimensional boundary value problems in hexagonal piezoelectric quasicrystals, Li et al. [17] applied the rigorous operator theory. Yu et al. [18] and Yu et al. [19] used techniques of complex variable theory for studying a stationary uniformly loaded antiplane Griffith crack and an elliptical cavity, respectively, in a one-dimensional hexagonal piezoelectric quasicrystal.Within the last two years, the burgeoning depth of research into the effects of piezoelectricity in quasicrystals has continued. The internal and interfacial Green’s functions of piezoelectric quasicrystalline bi-materials were obtained by Zhang et al. [20]. Guo and Pan [21] proposed a three-phase cylinder model for composites of piezoelectric quasicrystals. Fundamental solutions for three-dimensional cracks were established by Fan et al. [22] using unit point extended displacement discontinuities and the boundary integral equation method. By means of conformal mappings and complex variables methods, Yang and Li [23] found analytical solutions for a crack emanating from a circular inclusion and Guo et al. [24] studied a composite matrix containing an elliptical inclusion. Both works give the corresponding solutions for a uniformly loaded Griffith crack as special cases. Explicit expressions for the field components of a moving nonconstantly loaded antiplane single crack embedded in an infinite region and within a half-space of one-dimensional hexagonal piezoelectric quasicrystals were determine by Tupholme [25, 26]. Most recently, Fourier transforms and triple integral equations were used by Zhou and Li [27] to derive solutions in terms of complete elliptical integrals of the first and third kinds for two collinear, constantly loaded, shear cracks in a strip of the piezoelectric quasicrystal.Within purely piezoelectric media, in which no phason fields are present, many analyses of various crack configurations have been undertaken. It is of note that, in particular, Han and Wang [28] considered a periodic row of parallel, mode I, uniformly loaded cracks situated perpendicular to the surfaces of a piezoelectric strip, by reducing this mixed boundary value problem to the numerical solution of a system of hypersingular integral equations. A periodic distribution of constantly loaded mode III slant cracks was analyzed by Li and Lee [29] using a collocation method and quadrature formula to solve numerically the governing singular integral equations. Zhou et al. [30] studied multiple equally spaced, parallel antiplane shear cracks in infinite piezoelectric materials subjected to constant shear stress loading. Fourier transforms led to dual integral equations which were expanded as a series of Jacobi polynomials and then solved numerically using the Schmidt method. Recently, Monfared and Ayatollahi [31] investigated piezoelectric solids weakened by multiple constantly loaded cracks using Fourier transforms that yielded a set of integral equations for a numerical solution. They provided a useful quite comprehensive review of previous related works of this type.In Section2, the general equations in three dimensions that govern deformed piezoelectric quasicrystals are formulated and the physical stack of cracks problem being considered here is stipulated. The relevant constitutive equations for one-dimensional hexagonal piezoelectric quasicrystals with a point group 6 mm are displayed. The properties of a “piezoelectric quasicrystal screw dislocation” which underpin the current investigation are then outlined in Section 3. A necessary adaptation of the dislocation layer technique is exploited in Section 4 to generate the singular integral equations satisfied by the density functions of the dislocation arrays. These lead to closed-form representations for the resulting phonon, phason, and electric fields and intensity factors created. Their variations with the polar angle near a crack tip are deduced, and graphical numerical results are displayed. In the concluding Section 5, the principal features of this analysis are summarized.


## 2. Fundamental Governing Equations of Piezoelectric Quasicrystals and the Stack of Cracks Problem Formulation

In differential and variationally invariant forms, Altay and Dökmeci [13] presented and studied the governing three-dimensional general equations for the deformation components of the phonon, phason, and electric fields within a linearized theory of piezoelectric quasicrystalline continua. With repeated suffices i,j,k,l=1,2,3 summed according to the traditional Einstein convention, in relation to a rectangular system of fixed Cartesian coordinates x1,x2,x3, the appropriate quasi-static equilibrium equations and constitutive equations of such media which are devoid of any applied body forces or electric charge densities can be expressed neatly as follows:(1)σij,i=0,Hij,i=0,Di,i=0,σij=cijkluk,l+ul,k2+Rijklwk,l−ekijEk,Hij=Rklijuk,l+ul,k2+Kijklwk,l−ekij′Ek,Di=ekijuj,k+uk,j2+ekij′wj,k−εijEj,where p after a comma, for p=i,j,k,l, denotes partial differentiation with respect to xp. ui and σij are the components of the phonon displacement vector and stress tensor, wi and Hij are the components of the phason displacement vector and stress tensor, and Di and Ei are the components of the electric displacement and field. Respectively, cijkl, Rijkl, Kijkl, eijk, eijk′, and cijkl are the phonon elastic, the phonon-phason coupling, the phason elastic, the phonon and phason piezoelectric and the dielectric constants.A one-dimensional hexagonal piezoelectric quasicrystal medium of a point group 6 mm which is initially in a reference undisturbed stress-free state is studied here. With respect to a fixed Cartesian system of coordinates (x, y, z), the positive z-axis coincides with its quasiperiodic direction and it is periodic in the x-y plane.The material contains an infinite stack of equally spaced, Griffith-type shear cracks parallel to thex-axis, each of width 2a, that is subjected at infinity to nonuniform phonon, phason, and electrical loads which create a mode III deformation.The constitutive equations describing the interrelationships between the componentsuX, εXY, σXY, wX, wzX, HzX, DX, and EX, respectively, of the phonon displacement vector, phonon strain tensor, phonon stress tensor, phason displacement vector, phason strain tensor, phason stress tensor, electric displacement vector, and electric field vector, for X and Y=x,y,orz, for deformations of the piezoelectric quasicrystal can be elegantly written in the matrix notation as follows:(2)σxxσyyσzzσyzσxzσxyHzzHzxHzy=c11c12c13000R100c12c11c13000R100c13c13c33000R2000002c440000R300002c4400R3000000c11−c12000R1R1R2000K10000002R300K200002R30000K2εxxεyyεzzεyzεxzεxywzzwzxwzy−00e3100e3100e330e150e150000000e33′e15′000e15′0ExEyEz,(3)DxDyDz=00002e1500e15′00002e150000e15′e31e31e33000e′3300εxxεyyεzzεyzεxzεxywzzwzxwzy+ε11000ε11000ε33ExEyEz,with(4)εXY=12∂uX∂Y+∂uY∂X,wzX=∂wz∂X,Here, the material constants involved have been abbreviated, using the standard contracted notation of Voigt with i and j having integer values, to cij, the phonon elastic moduli, Ki, the phason elastic moduli, Ri, the phonon-phason coupling elastic moduli, eij and eij′, the piezoelectric moduli, and εij, the dielectric moduli.It is noteworthy to mention at this stage that no previous investigations have been presented deriving explicit expressions for the deformed field components, or their angular variations near the tip of a crack, within a stack of nonuniformly loaded parallel cracks, even in hexagonal piezoelectric crystals or purely isotropic elastic, nor even when the loads are uniform. The study here is for piezoelectric quasicrystals, which are the most general material currently available for which the analysis is amenable. The results can be reduced to those for hexagonal piezoelectric crystals with no phason field involved by throughout puttingR=K=e15′=0 and then those for purely isotropic elastic materials with no piezoelectric effects by also putting e15=0, with c44=μ, the second Lame constant.The cracks are considered to be parallel to each other and a distanceh apart, so that they are centred at x=0 on the planes y=nh for n=0,±1,±2,…, as shown in Figure 1.Figure 1
Stack of parallel Griffith-type shear cracks.It is convenient to write the electric field vector,E, in terms of an electric potential, ϕ, so that(5)E=−∇ϕ.All the field variables in the antiplane deformation which is created are independent ofz. Thus, from Equation (4), the nonzero phonon and phason strain components are given in terms of the components, uz and wz, of the phonon and phason displacements by(6)εxz=12∂uz∂x,εyz=12∂uz∂y,wzx=∂wz∂x,wzy=∂wz∂y,and, from the constitutive equations (2) and (3), the corresponding nonzero phonon and phason stresses and displacement components are(7)σyz=2c44εyz+Rwzy−e15Ey,σxz=2c44εxz+Rwzx−e15Ex,(8)Hzy=2Rεyz+Kwzy−e15′Ey,Hzx=2Rεxz+Kwzx−e15′Ex,(9)Dy=2e15εyz+e15′wzy+ε11Ey,Dx=2e15εxz+e15′wzx+ε11Ex.In Equations (7)–(9) and subsequently, for the sake of brevity, the material constants R3 and K2 are simply written as R and K.The nonuniform phonon, phason, and electrical loads,Tx,Hx, and Dx, applied at infinity are symmetrical, so that(10)σyzx,y⟶Tx,Hzyx,y⟶Hx,Dyx,y⟶Dx,asx2+y2⟶∞,with T−x=Tx, H−x=Hx, and D−x=Dx for all x.Likewise, an interested reader could analyze analogously the mode III deformations which arise by specifying instead any other selection of three of the componentsσyz,εyz,Hzy,wzy,Dy, or Ey at infinity.


## 3. A Piezoelectric Quasicrystal Screw Dislocation

It is convenient to first summarize the fundamental features of a “piezoelectric quasicrystal screw dislocation” upon which the subsequent analysis here is founded. This extends to a piezoelectric quasicrystalline medium the traditional notion of a screw dislocation in a linear solely elastic solid material by the introduction of the concept of an extended Burgers vector that has discontinuities in the components,uzandwz, of the phonon and phason displacements and ϕ, the electric potential, of magnitudes b,d, and b4, respectively, across a slip plane.Wang and Pan [15] presented explicit general expressions for the components of the fields of this type of straight dislocation situated parallel to the z-axis at the origin in a one-dimensional hexagonal piezoelectric quasicrystal with a point group 6 mm. From these, it can be deduced that, for such a stationary screw dislocation, the phonon and phason displacement components and the electric potential can be written as(11)uzIIIx,y=b2πtan−1yx,wzIIIx,y=d2πtan−1yx,φIIIξ,y=b42πtan−1yx,where a superscript III is attached to the field variables to indicate that they relate to an antiplane mode III deformation.For convenience of presentation, it is appropriate to introduce the piezoelectrically stiffened phonon and phason elastic moduli,c¯44 and K¯, and the phonon-phason coupling elastic constant, R¯, as defined in the following equations:(12)c¯44=c44+e152ε11,K¯=K+e15′2ε11,R¯=R+e15e15′ε11.Recalling the relations (5) and (6), it follows by substituting the expressions (11) into (7)–(9) that the necessary components of the phonon and phason stresses and electric displacement are(13)σxzIIIx,y=−c44b+Rd+e15b42πyx2+y2,σyzIIIx,y=c44b+Rd+e15b42πxx2+y2,(14)HzxIIIx,y=−Rb+Kd+e15′b42πyx2+y2,HzyIIIx,y=Rb+Kd+e15′b42πxx2+y2,(15)DxIIIx,y=−e15b+e15′d−ε11b42πyx2+y2,DyIIIx,y=e15b+e15′d−ε11b42πxx2+y2.Likewise, if required, the components of the phonon and phason strain and electric field can be deduced using Equations (6), (11), and (5).


## 4. Stack of Shear Cracks Analysis Using an Extension of the Dislocation Layers Technique

Initially for isotropic elastic solids, the “dislocation layer technique” was devised for modelling loaded strip cracks by equivalent continuous distributions of elastic dislocations which replicate the displacement jump across the crack faces, as usefully described, for example, by Bilby and Eshelby [32] and Lardner [33]. Appropriate arrays of piezoelectric quasicrystal screw dislocations are determined here to extend this concept to the stack of cracks in piezoelectric quasicrystals under consideration.The screws on the portions to the right of the origin on each crack are positive and those to the left of the origin are negative. The discontinuities in the phonon and phason displacement components and the electric potential must therefore have density functions,fx, gx, and jx, respectively, which are odd functions of x.Because of the symmetry of the problem, the density functions will be identical on each of the cracks, and moreover, it is thus only necessary to ensure that the “central” crack of the stack aty=0 is free of phonon and phason stresses and electric displacement.From Equations (13)–(15), at a general point on the x-axis, these dislocation arrays produce phonon and phason stresses and electric displacement having components:(16)σyzIIIx,0=bc442π∫−aa∑n=−∞∞x−x′x−x′2+n2h2fx′dx′+dR2π∫−aa∑n=−∞∞x−x′x−x′2+n2h2gx′dx′+b4e152π∫−aa∑n=−∞∞x−x′x−x′2+n2h2jx′dx′,(17)HzyIIIx,0=bR2π∫−aa∑n=−∞∞x−x′x−x′2+n2h2fx′dx′+dK2π∫−aa∑n=−∞∞x−x′x−x′2+n2h2gx′dx′+b4e′152π∫−aa∑n=−∞∞x−x′x−x′2+n2h2jx′dx′,(18)DyIIIx,0=be152π∫−aa∑n=−∞∞x−x′x−x′2+n2h2fx′dx′+de′152π∫−aa∑n=−∞∞x−x′x−x′2+n2h2gx′dx′−b4ε112π∫−aa∑n=−∞∞x−x′x−x′2+n2h2jx′dx′.The equations of equilibrium that necessarily must be fulfilled to satisfy the imposed boundary conditions (10) are that(19)σyzIIIx,0=−Tx,HzyIIIx,0=−Hx,DyIIIx,0=−Dx,forx<a.By substituting the expressions (16)–(18) into the conditions (19), a system of three coupled simultaneous equations is obtained. After suitable algebraic manipulation, its solutions lead to three singular integral equations to be satisfied by the density functions, fξ, gξ, and jξ, in the following forms:(20)∫−aa∑n=−∞∞x−x′x−x′2+n2h2fx′dx′=−2πbε11c¯44K¯−R¯2×ε11K¯Tx−ε11R¯Hx−e15′R¯−e15K¯Dx,(21)∫−aa∑n=−∞∞x−x′x−x′2+n2h2gx′dx′=−2πdε11c¯44K¯−R¯2×−ε11R¯Tx+ε11c¯44Hx+e15′c¯44−e15R¯Dx,(22)∫−aa∑n=−∞∞x−x′x−x′2+n2h2jx′dx′=−2πb4ε11c¯44K¯−R¯2×−e′15R¯−e15K¯Tx+e′15c¯44−e15R¯Hx+c44K−R2Dx.The combined contributions to the integrands in Equations (20)–(22) from each of the individual cracks in the stack, as represented by the sum of the doubly-infinite series, can be neatly evaluated using the identity(23)∑n=−∞∞zz2+n2h2=πhcothπzh,and the property that the density functions are odd functions. This enables the left-hand side of Equation (20), for example, to be rewritten as(24)πh∫−aacoshπx′/hsinhπx/h−sinhπx′/hfx′dx′.With the introduction of new variables given by(25)x2=sinhπxh,a2=sinhπah,x2′=sinhπx′h,the integral equation (20) for fx thus becomes(26)∫−a2a2f2x2′x2−x2′dx2′=−2πbε11c¯44K¯−R¯2×ε11K¯Tx−ε11R¯Hx−e15′R¯−e15K¯Dx,forx2<a2,where f2x2=f2sinhπx/h≡fx. In the interpretation of this improper integral, its Cauchy principal value is needed in accordance with the Plemelj formulae. Classical techniques, as expounded by, for example, Muskhelishvili [34] and Gakhov [35], readily yield its suitable solution in the form(27)f2x2=−2πbε11c¯44K¯−R¯21a22−x221/2∫−a2a2a22−x′221/2x′2−x2×ε11K¯Tx2′−ε11R¯Hx2′−e15′R¯−e15K¯Dx2′dx2′.Likewise, it is found that the phason and electric densities that fulfill Equations (21) and (22) are(28)g2x2=−2πdε11c¯44K¯−R¯21a22−x221/2∫−a2a2a22−x′221/2x2′−x2×−ε11R¯Tx2′+ε11c¯44Hx2′+e15′c¯44−e15R¯Dx2′dx2′,(29)j2x2=−2πb4ε11c¯44K¯−R¯21a22−x221/2∫−a2a2a22−x′221/2x2′−x2×−e15′R¯−e15K¯Tx2′+e15′c¯44−e15R¯Hx2′+c44K−R2Dx2′dx2′,where g2x2=g2sinhπx/h≡gx and j2x2=j2sinhπx/h≡jx. It is of note that, generally, these above density functions all depend upon the loads Tx, Hx, and Dx that are applied, the values of the geometric constants a and h, and the data for the piezoelectric quasicrystal moduli.Having determined representations for the densities required for the boundary conditions to be satisfied, closed-form expressions can be derived from Equations (6), (10), (11), (13)–(15), and (27)–(29) for all the phonon, phason, and electric field components of interest.To illustrate the kind of analysis which is involved, it can be seen from Equations (10) and (13), for example, that(30)σyzx,y=Tx+bc442π∫−aa∑n=−∞∞x−x″x−x″2+y−nh2fx″dx″+dR2π∫−aa∑n=−∞∞x−x″x−x″2+y−nh2gx″dx″+b4e152π∫−aa∑n=−∞∞x−x″x−x″2+y−nh2jx″dx″,(31)σxzx,y=−bc442π∫−aa∑n=−∞∞y−nhx−x″2+y−nh2fx″dx″−dR2π∫−aa∑n=−∞∞y−nhx−x″2+y−nh2gx″dx″−b4e152π∫−aa∑n=−∞∞y−nhx−x″2+y−nh2jx″dx″.At first sight, substitution directly of the expressions for the densities from Equations (27)–(29) into the above leads to very unwieldy double integral representations for these components of stress. But in fact, as summarized in Equations (A.1)–(A.11) of the appendix, after the x″ integration has been evaluated first and much rearrangement and cancellation algebraically, it is found that concisely(32)σyzx,y=Tx+FTθ,σxzx,y=−FTθ−π2,where the functions FFθ, Rx,y, and θx,y are defined for Fx=Tx,Hx, and Dx by(33)FFθ=−1πR∫−a2a2M1x,x2′,θM2x,x2′a22−x2′21/2Fx2′dx2′,(34)Reiθ=a22−sinh2πx+iyh1/2,with(35)M1x,x2′,θ=sinhπxhcosπyh−x2′sinθ+coshπxhsinπyhcosθ,(36)M2x,x2′=sinhπxhcosπyh−x2′2+cosh2πxhsin2πyh.The branches of the square root function in Equation (34) are chosen to specify that θ is zero for −a<x<aandy=0+ and extended elsewhere by analytic continuation.Similarly utilizing Equations (14) and (15) yields(37)Hzyx,y=Hx+FHθ,Hzxx,y=−FHθ−π2,(38)Dyx,y=Dx+FDθ,Dxx,y=−FDθ−π2.With the imposed boundary conditions (10), it can be seen from Equations (32)–(38) and recalling Equation (25) that the phonon stress components, phason stress components, and the electric displacement components individually are decoupled and depend, respectively, upon the applied loads, Tx,Hx, and Dx, but that they are all dependent upon the geometric parameters a and h.Close to a crack tip, the practically interesting variation of these field components with the polar angle can be studied by defining polar coordinatesr and ψ, with r≪a, so that(39)x=a+rcosψ,y=rsinψ.Then, asr⟶0, Equation (34) shows that approximately(40)R∼r2πhsinhπahcoshπah1/2,θ∼−π−ψ2,and correspondingly, Equations (32)–(38) yield(41)σyzr,ψ∼KTSrcosψ2,Hzyr,ψ∼KHSrcosψ2,Dyr,ψ∼KDSrcosψ2,(42)σxzr,ψ∼−KTSrsinψ2,Hzxr,ψ∼−KHSrsinψ2,Dxr,ψ∼−KDSrsinψ2,where the phonon and phason stress and electric displacement intensity factors, KTS, KHS, and KDS, at the tip of a crack within a stack of cracks, are defined for Fx=Tx,Hx, and Dx by(43)KFS=12πhsinhπa/hcoshπa/h1/2∫−aasinhπa/h+sinhπx′/hsinhπa/h−sinhπx′/h1/2coshπx′hFx′dx′.It can be observed thatKFS⟶KF in the limit as h⟶∞, where(44)KF=1π2a∫−aaa+x′a−x′1/2Fx′dx′,which is in agreement with the intensity factors derived by Tupholme [25] for a single, nonuniformly loaded Griffith-type shear crack.Furthermore, the polar components of the phonon stress, phason stress, and electric displacement,σψz, Hzψ, and Dψ, respectively, near a crack tip are deduced from Equations (41) and (42) to be(45)σψzr,ψ∼KTSrcosψ2,Hzψr,ψ∼KHSrcosψ2,Dψr,ψ∼KDSrcosψ2,with their variation with the nonuniform loads contained entirely in the intensity factors given by Equation (43).The distribution around a crack tip with the angleψ of the scaled phonon stress component, rσψz/KTS, as evaluated from Equation (45), is depicted in Figure 2.Figure 2
Variation of the scaled phonon stress component,rσψz/KTS, around a crack tip.Clearly from the above, the phonon and phason stress and electric displacement components exhibit the same1/r crack-tip behaviour as that close to the tip of a loaded crack embedded in a classical purely elastic isotropic medium.It is also observed from Equations (27)–(29) that all the density functions are dependent not only upon the geometrical parameters a and h and the imposed loads Tx, Hx, and Dx but also upon the piezoelectric quasicrystal material constants. But on the other hand, it is seen from Equations (41)–(43) and (45) that the components of the phonon stress, phason stress, and electric displacement are decoupled and that they and the field intensity factors solely depend on the respective loads Tx, Hx, and Dx, as well as the geometric parameters, a and h. This phenomenon of decoupling and independence of material constants was also exhibited in the study of parallel constantly loaded periodic cracks in piezoelectric media by Li and Lee [29].In contrast, however, if there is an interest in the phonon and phason displacement components and the electric potential in the vicinity of a crack in the stack, it is clear that if they are deduced using Equations (11) and (27)–(29), they will be heavily related additionally to the material moduli.When all the loads applied to the stack of cracks areuniform, so that Fx=F=constant, the integral in Equation (43) can be calculated exactly and then KFS reduces to(46)KFS=h2πtanhπah1/2F.This is in agreement with the results of Li and Lee [29] for piezoelectric media.Correspondingly, for a single uniformly loaded crack, the expression (44) can be calculated exactly as(47)KT=a2T,KH=a2H,KD=a2D.Thus, withuniform loads Tx=T=constant, Hx=H=constant, and Dx=D=constant on the stack of cracks,(48)KFS=tanhπa/hπa/h1/2KF,for F=T,H, and D, so in this situation each of the field intensity factors at a crack tip in the stack is of smaller magnitude than that near a corresponding single uniformly loaded crack. This so-called crack shielding behaviour is illustrated in Figure 3, where the variation with the relative spacing h/a of the normalized field intensity factors, KFS/KF, of a uniformly loaded stack is plotted. It is seen to occur decreasingly as h/a increases from 0 to approximately 15. This phenomenon was exhibited similarly, for example, in the numerical investigation of multiple cracks by Zhou et al. [30].Figure 3
Variation of the normalized field intensity factors with the relative spacing in a stack of uniformly loaded cracks.It is plain from Equations (27)–(29) that the analyses above would be invalid when c¯44K¯−R¯2=0. But as yet, for one-dimensional hexagonal piezoelectric quasicrystals, the values presented of the moduli data are not entirely reliable. However, Li et al. [17] typically reported that c44=5.0×1010Nm−2, R=1.2×109Nm−2, K=3.0×108Nm−2, e15=−0.138Cm−2, e′15=−0.160Cm−2, and ε11=82.6×10−12C2N−1m−2. Then, clearly c¯44K¯−R¯2 is nonzero, since the value of c¯44K¯ is much greater than that of R¯2.


## 5. Conclusions

The fields of deformation generated around a stack of nonuniformly loaded parallel mode III shear cracks in one-dimensional piezoelectric quasicrystals are studied by appropriately extending the technique of dislocation layers for simulating the cracks.The system of three coupled singular integral equations which must be fulfilled by the density functions of the required dislocation arrays in order to satisfy the imposed boundary conditions is derived. These are neatly manipulated into singular integral equations that are amenable to solution exactly, using an appropriate explicit hyperbolic trigonometric function identity representation of the doubly-infinite series that are involved as summations of the contributions from all the cracks. Explicit expressions for the components of the phonon, phason, and electric fields are derived in a closed form, and their angular variations around a crack tip are deduced. Elegant representations for the field intensity factors are deduced.Even for a stack of uniformly loaded cracks, general representations such as these have not been presented before for any materials whatsoever. For example, in the special case of a stack of constantly loaded cracks in just piezoelectric media, only the field components along the crack lines, and thus the intensity factors, have been derived. Comparisons are drawn between the nonuniformly loaded and uniformly loaded situations.It is demonstrated that the phenomenon of crack shielding occurs for small values of the relative spacingh/a, by illustrating graphically the variation of the ratio of the intensity factors for a stack of uniformly loaded cracks to those of a corresponding single crack.As particular cases, of course, the previously unavailable corresponding simplified results for non-quasicrystals or non-piezoelectric or purely elastic media follow.Representative numerical data are displayed graphically.



---

*Source: 4847837-2018-09-23.xml*