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

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.)eir dependence upon the polar angle in the region close to the tip of a crack is deduced, and the field intensity factors then follow. )ese exhibit that the phenomenon of crack shielding is dependent upon the relative spacing of the cracks. )e analogous analyses, that have not been given previously, involving nonpiezoelectric 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.


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 to nonuniform mechanical and electric loads in onedimensional 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.e corresponding results when the material is simply quasicrystalline or solely piezoelectric or purely isotropic elastic can be deduced as special cases.e discovery of some unusual complex metallic alloys having quasiperiodic structures was first widely announced by Shechtman et al. [3].ey are now known universally as quasicrystals.e 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.us, 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.
e underlying equations that govern the linear elastic theory of quasicrystals within the framework of the mechanics of continua are now widely established.
eir 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 threedimensional equations of elasticity of quasicrystals to include quasi-static electric effects.
is has enabled piezoelectric quasicrystal boundary value problems to be addressed.
e 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 onedimensional 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 threedimensional 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.
e 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.
ey provided a useful quite comprehensive review of previous related works of this type.
In Section 2, 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.e relevant constitutive equations for one-dimensional hexagonal piezoelectric quasicrystals with a point group 6 mm are displayed.e 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.
ese lead to closed-form representations for the resulting phonon, phason, and electric fields and intensity factors created.
eir 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.

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 threedimensional general equations for the deformation components of the phonon, phason, and electric fields within 2 Advances in Materials Science and Engineering 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 (x 1 , x 2 , x 3 ), 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: where p after a comma, for p � i, j, k, l, denotes partial differentiation with respect to x p .u i and σ ij are the components of the phonon displacement vector and stress tensor, w i and H ij are the components of the phason displacement vector and stress tensor, and D i and E i are the components of the electric displacement and field.Respectively, c ijkl , R ijkl , K ijkl , e ijk , e ijk ′ , and c ijkl 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.
e material contains an infinite stack of equally spaced, Griffith-type shear cracks parallel to the x-axis, each of width 2a, that is subjected at infinity to nonuniform phonon, phason, and electrical loads which create a mode III deformation.
e constitutive equations describing the interrelationships between the components u X , ε XY , σ XY , w X , w zX , H zX , D X , and E X , 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, or z, for deformations of the piezoelectric quasicrystal can be elegantly written in the matrix notation as follows: with Here, the material constants involved have been abbreviated, using the standard contracted notation of Voigt with i and j having integer values, to c ij , the phonon elastic moduli, K i , the phason elastic moduli, R i , the phonon-phason coupling Advances in Materials Science and Engineering elastic moduli, e ij and e ij ′ , 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 eld 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.e study here is for piezoelectric quasicrystals, which are the most general material currently available for which the analysis is amenable.
e results can be reduced to those for hexagonal piezoelectric crystals with no phason eld involved by throughout putting R K e 15 ′ 0 and then those for purely isotropic elastic materials with no piezoelectric e ects by also putting e 15 0, with c 44 μ, the second Lame constant.
e cracks are considered to be parallel to each other and a distance h apart, so that they are centred at x 0 on the planes y nh for n 0, ±1, ±2, . . ., as shown in Figure 1.
It is convenient to write the electric eld vector, E, in terms of an electric potential, ϕ, so that E −∇ϕ. ( All the eld variables in the antiplane deformation which is created are independent of z. us, from Equation ( 4), the nonzero phonon and phason strain components are given in terms of the components, u z and w z , of the phonon and phason displacements by In Equations ( 7)-( 9) and subsequently, for the sake of brevity, the material constants R 3 and K 2 are simply written as R and K. e nonuniform phonon, phason, and electrical loads, T(x), H(x), and D(x), applied at in nity are symmetrical, so that with T(−x) T(x), H(−x) H(x), and D(−x) D(x) 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 , H zy , w zy , D y , or E y at in nity.

A Piezoelectric Quasicrystal Screw Dislocation
It is convenient to rst summarize the fundamental features of a "piezoelectric quasicrystal screw dislocation" upon which the subsequent analysis here is founded.is 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, u z and w z , of the phonon and phason displacements and 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 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: 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  Likewise, if required, the components of the phonon and phason strain and electric field can be deduced using Equations ( 6), (11), and (5).

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.e 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.
e discontinuities in the phonon and phason displacement components and the electric potential must therefore have density functions, f(x), g(x), and j(x), 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 at y � 0 is free of phonon and phason stresses and electric displacement.
From Equations ( 13)-( 15), at a general point on the xaxis, these dislocation arrays produce phonon and phason stresses and electric displacement having components: Advances in Materials Science and Engineering e equations of equilibrium that necessarily must be fulfilled to satisfy the imposed boundary conditions (10) are that 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: e 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 and the property that the density functions are odd functions.is enables the left-hand side of Equation (20), for example, to be rewritten as With the introduction of new variables given by the integral equation (20) for f(x) thus becomes where 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 Likewise, it is found that the phason and electric densities that fulfill Equations ( 21) and ( 22) are where g 2 (x 2 ) � g 2 (sinh (πx/h)) ≡ g(x) and j 2 (x 2 ) � j 2 (sinh (πx/h)) ≡ j(x).It is of note that, generally, these above density functions all depend upon the loads T(x), Advances in Materials Science and Engineering H(x), and D(x) that are applied, the values of the geometric constants a and h, and the data for the piezoelectric quasicrystal moduli.
To illustrate the kind of analysis which is involved, it can be seen from Equations ( 10) and ( 13), for example, that 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 where the functions F F (θ), R(x, y), and θ(x, y) are defined for F(x) � T(x), H(x), and D(x) by with e branches of the square root function in Equation ( 34) are chosen to specify that θ is zero for −a < x < a and y � 0+ and extended elsewhere by analytic continuation.
Similarly utilizing Equations ( 14) and ( 15) yields 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, T(x), H(x), and D(x), 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 coordinates r and ψ, with r ≪ a, so that en, as r ⟶ 0, Equation (34) shows that approximately and correspondingly, Equations ( 32)-(38) yield Advances in Materials Science and Engineering where the phonon and phason stress and electric displacement intensity factors, K S T , K S H , and K S D , at the tip of a crack within a stack of cracks, are de ned for F(x) T(x), H(x), and D(x) by It can be observed that K S F ⟶ K F in the limit as h ⟶ ∞, where which is in agreement with the intensity factors derived by Tupholme [25] for a single, nonuniformly loaded Gri thtype shear crack.Furthermore, the polar components of the phonon stress, phason stress, and electric displacement, σ ψz , H zψ , and D ψ , respectively, near a crack tip are deduced from Equations ( 41) and (42) to be with their variation with the nonuniform loads contained entirely in the intensity factors given by Equation ( 43). e distribution around a crack tip with the angle ψ of the scaled phonon stress component, r √ σ ψz /K S T , as evaluated from Equation (45), is depicted in Figure 2.
Clearly from the above, the phonon and phason stress and electric displacement components exhibit the same 1/ 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 T(x), H(x), and D(x) 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 eld intensity factors solely depend on the respective loads T(x), H(x), and D(x), as well as the geometric parameters, a and h. is 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 are uniform, so that F(x) F constant, the integral in Equation (43) can be calculated exactly and then K S F reduces to is 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) us, with uniform loads T(x) T constant, H(x) H constant, and D(x) D constant on the stack of cracks, for F T, H, and D, so in this situation each of the eld intensity factors at a crack tip in the stack is of smaller magnitude than that near a corresponding single uniformly loaded crack. is so-called crack shielding behaviour is illustrated in Figure 3, where the variation with the relative spacing h/a of the normalized eld intensity factors, K S F /K F , of a uniformly loaded stack is plotted.It is seen to occur decreasingly as h/a increases from 0 to approximately 15. is phenomenon was exhibited similarly, for example, in the numerical investigation of multiple cracks by Zhou et al. [30].

Conclusions
e elds of deformation generated around a stack of nonuniformly loaded parallel mode III shear cracks in onedimensional piezoelectric quasicrystals are studied by appropriately extending the technique of dislocation layers for simulating the cracks.
e system of three coupled singular integral equations which must be ful lled by the density functions of the required dislocation arrays in order to satisfy the imposed boundary conditions is derived.ese 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-in nite series that are involved as summations of the contributions from all the cracks.Explicit expressions for the components of the phonon, phason, and electric elds are derived in a closed form, and their angular variations around a crack tip are deduced.Elegant representations for the eld 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 eld 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 spacing h/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 simpli ed results for non-quasicrystals or non-piezoelectric or purely elastic media follow.
Representative numerical data are displayed graphically.

Appendix
e components σ yz and σ xz of the phonon stress are given by the expressions ( 30) and (31), respectively, for example, and involve integrals of the types for r(x) f(x), g(x), and j(x).
It is particularly convenient to simultaneously calculate these by considering the complex combination K(x, y) of them de ned by (A.4) In rewriting Equation (A.3) in the form (A.4), the doubly-infinite summation has been evaluated using the identities and the property that all the density functions are odd, in a manipulation and with a notation which is reminiscent of that in which the identity (23) was used to rewrite Equation (20) as Equation (26).
When the representations for r 2 (x ″ 2 ) in Equations ( 27)-( 29) are incorporated into Equation (A.4), K(x, y) is seen to involve relevant constant multiples of double integrals of the type: where F(x) � T(x), H(x), or D(x).e inner integral in the above expression can be calculated using complex variable contour integration methods by applying Cauchy's residue theorem to evaluate the integral around a contour in the complex plane enclosing the singularities at x 2 ′ and sinh π(x + iy)/h   along the branch cuts from −a 2 to a 2 and a large circle around the origin to give where the square root function's branches are chosen as in Equation (34).e real and imaginary parts of the representation (A.7) can then straightforwardly be deduced to give the integrals I(x, y) and J(x, y) in terms of a function F F (θ), defined for F(x) � T(x), H(x), and D(x) by If preferred, alternatively, partial fractions, complex conjugate functions, and contour integration methods can used to calculate I(x, y) and J(x, y).

FFigure 3 : 1 x 2 􏼒
Figure 3: Variation of the normalized eld intensity factors with the relative spacing in a stack of uniformly loaded cracks.