Mathematical Problems in Engineering
Volume 2020 (2020), Article ID 4797413, 15 pages
https://doi.org/10.1155/2020/4797413

Research Article

Exact Solutions for Interaction of Parallel Screw Dislocations with a Wedge Crack in One-Dimensional Hexagonal Quasicrystal with Piezoelectric Effects

Xin Lv1,2 and Guan-Ting Liu¹

¹College of Mathematics Science, Inner Mongolia Normal University, Hohhot 010022, China
²School of Mechanical Engineering, Tianjin University, Tianjin 300350, China

Correspondence should be addressed to Guan-Ting Liu; guantingliu@imnu.edu.cn

Received 7 October 2019; Revised 13 April 2020; Accepted 21 April 2020; Published 29 May 2020

Academic Editor: Francesco Aymerich

Copyright © 2020 Xin Lv and Guan-Ting Liu. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The purpose of this paper is to consider the interaction between many parallel dislocations and a wedge-shaped crack and their collective response to the external applied generalized stress in one-dimensional hexagonal piezoelectric quasicrystal, employing the complex variable function theory and the conformal transformation method; the problem for the crack is reduced to the solution of singular integral equations, which can be further reduced to solving Riemann–Hilbert boundary value problems. The analytical solutions of the generalized stress field are obtained. The dislocations are subjected to the phonon field line force, phason field line force, and line charge at the core. The positions of the dislocations are arbitrary, but the dislocation distribution is additive. The dislocation is not only subjected to the external stress and the internal stress generated by the crack, but also to the force exerted on it by other dislocations. The closed-form solutions are obtained for field intensity factors and the image force on a screw dislocation in the presence of a wedge-shaped crack and a collection of other dislocations. Numerical examples are provided to show the effects of wedge angle, dislocation position, dislocation distribution containing symmetric configurations and dislocation quantities on the field intensity factors, energy release rate, and image force acting on the dislocation. The principal new physical results obtained here are (1) the phonon stress, phason stress, and electric displacement singularity occur at the crack tip and dislocations cores, (2) the increasing number of dislocations always accelerates the crack propagation, (3) the effect of wedge angle on crack propagation is related to the distribution of dislocations, and (4) the results of the image force on the dislocation indicate that the dislocations can either be attracted or rejected and reach stable positions eventually.

1. Introduction

As a new structure of solid, quasicrystal (QC) was discovered by Shechtman et al. [1] in 1984. QC is a new solid structure with long-range orientation order and without translational symmetry which differs from crystals. The elastic problem of QC is related not only to the phonon field which describes lattice vibration, but also to the phason field which is about the quasiperiodic arrangement of atoms. As soon as the quasicrystal was discovered, it attracted wide attention because of its special symmetry and good performance. In the early discovery of the quasicrystal, the defects in the quasicrystal were observed [2]. Due to the unique properties of quasicrystal materials, many achievements have been made in the study of quasicrystal elasticity and fracture mechanics [3–9]. To name a few recent researches, Lou et al. [10] proposed a model of a thin elastic inclusion embedded in an infinite 1D hexagonal quasicrystal. Wang and Schiavone [11] investigated the elastic field near the tip of an anticrack in a homogeneous decagonal quasicrystal material subjected to plane strain deformations. Zhao et al. [12] studied the frictional contact problem and the adhesive contact problem in three-dimensional (3D) icosahedral quasicrystal. Li et al. [13] have given analytical solutions for arbitrarily shaped planar cracks in a two-dimensional (2D) hexagonal QC media.

Quasicrystal is sensitive to force, electricity, magnetism, heat, and light, and these physical properties have been extensively studied. It is well known that QCs with piezoelectric effects have better performance in smart composite structures. Therefore, the study of quasicrystal materials considering the piezoelectric effect is of great significance. A lot of work on QC with piezoelectric effects has been published in the open literature. Li et al. [14] obtained three-dimensional general solutions to static problems of a 1D hexagonal piezoelectric QC. Yu et al. [15] considered an antiplane crack problem by the semi-inverse method in a 1D hexagonal piezoelectric QC. Yu et al. [16] also discussed the complex variable method for an antiplane elliptical cavity. Yang et al. [17] presented the fracture behavior of two asymmetrical limited permeable cracks emanating from an elliptical hole in a one-dimensional hexagonal quasicrystal with piezoelectric effects. Yang et al. [18] derived the closed-form solution for a multilayered two-dimensional decagonal QC plate. Guo et al. [19] analyzed an elliptical inclusion embedded in an infinite 1D hexagonal piezoelectric QC matrix. Li et al. [20] have given the exact solution for a layered piezoelectric quasicrystal nanoplate with nonlocal effect in cylindrical bending. Zhang et al. [21] addressed the static bending deformation of a functionally graded multilayered 1D hexagonal piezoelectric QC simply supported nanoplate. Zhou and Li [22] dealt with the problem of two collinear cracks in a one-dimensional hexagonal piezoelectric QC strip under mode III far-field loading, while they analyzed the electroelastic problem related to two collinear cracks of equal length and normal to the boundaries of a 1D hexagonal piezoelectric quasicrystal layer with the aid of a singular integral equation in ref. [23]. Besides, Zhou and Li [24] examined the electroelastic behavior induced by a penny-shaped dielectric crack of a 1D hexagonal piezoelectric quasicrystal plate of finite thickness, and they considered a Yoffe-type moving crack in a 1D hexagonal piezoelectric quasicrystal by using the Fourier transform technique in ref. [25]. Li and Liu [26] derived the electroelastic field of a straight dislocation in an icosahedral piezoelectric QC. Li et al. [27] displayed the interaction between a screw dislocation and an elliptical hole with two asymmetrical cracks in a 1D hexagonal quasicrystal with piezoelectric effects. Hu et al. [28] studied an antiplane interface crack between dissimilar 1D hexagonal piezoelectric QC.

Engineering materials always contain some defects, such as dislocations, inclusions, and cracks. The existence of these defects often seriously affects the safety and reliability of materials. There is not only one defect in the material, but usually different types of defects also exist at the same time, and the interaction between multiple defects has a very important impact on the mechanical properties and physical properties of the material. Dislocation activity is the cause of plastic deformation at the crack tip, which determines the crack propagation behavior, so the interaction between crack and dislocation has attracted more and more attention. Concerning the wedge-shaped crack which is an important geometric form of structures interacting with a single dislocation, Ohr et al. [29] performed the elastic analysis of the interaction between a wedge crack and a screw dislocation. Zhang et al. [30] obtained an analytical solution of the interaction of edge dislocation with the wedge crack by using the elastic technique of Muskhelishvili. Chen et al. [31] derived the closed solution of the electroelastic field for a screw dislocation situated near the tip two-bonded wedge-shaped dissimilar piezoelectric materials. Furthermore, they deduced the electroelastic stress field in the interaction of a piezoelectric screw dislocation with a wedge-shaped bimaterial interface in ref. [32]. Li and Zhao [33] coped with the problem of the interaction of a screw dislocation with the interface and wedge-shaped cracks in a one-dimensional hexagonal piezoelectric QC biomaterial using the conformal mapping method in conjunction with the image principle. Jiang and Liu [34] investigated the interaction between a screw dislocation and a wedge-shaped crack in the 1D hexagonal piezoelectric QC based on the conformal mapping method and the perturbation technique.

The research studies mentioned above were mainly focused on the single screw dislocation interacting with the various types of cracks. In engineering practice, structures with pairs of dislocations are very common. Dislocation emission and dislocation pile-up at the crack tip will result in crack blunting, inhibiting the growth of the crack and then improving the fracture toughness of the material. Concerning multiple dislocations interacting with cracks, the dislocation is not only subjected to the external stress and the internal stress generated by the crack, but also to the force exerted on it by other dislocations. Lin and Thomson [35] earlier investigated the elastic interaction between pairs of dislocations and a cleavage crack under external loads. Later, Juang and Lee [36] examined the interaction of the two-dislocation model and a surface crack by using the dislocation-modelling technique and the complex variable method. More recently, Liu and Yang [37] made a detailed analysis of the stress field of infinitely many parallel screw dislocations near a semi-infinite crack tip without applied loads, and dislocations considered here are located at particular points to make the total stress finite. Lv and Liu [38] addressed the analytical solution of the generalized stress field for parallel screw dislocations near a semi-infinite crack in magnetic-electric-mechanical coupling fields. However, for the piezoelectric QC, these problems become more complicated than the conventional material due to the introduction of the phason field and the electric field. To the best of the authors’ knowledge, the study on the interaction between the wedge-shaped crack and multiple dislocations in a 1D hexagonal QC with piezoelectric effect has not been reported yet. In this paper, through the singular integral technique, we consider the antiplane problem of the interaction between many parallel screw dislocations and a wedge-shaped crack under remote uniform loadings, and the dislocations are subjected to the phonon field line force, phason field line force, and line charge at the core. The positions of the dislocations are arbitrary. The advantage of the present model is demonstrated when dealing with many dislocations because the dislocation distribution is additive. The expressions for the generalized stress field, field intensity factors, energy release rate, and image force are derived in a closed form. The exact solutions are presented and compared with those already obtained by other authors. Numerical examples are presented for various parameters of the problem.

2. Basic Equations

The quasicrystal piezoelectric material considered here is transversely isotropic with a hexagonal symmetry, which has an isotropic basal plane of plane and an electric field poling direction and quasiperiodic direction of the axis. When the defects penetrate along the axis, the quasicrystal piezoelectric boundary value problem is treated as in the case of out-of-plane phonon field displacement , out-of- plane phason field displacement , and in-plane electric field . The basic equations can be written as follows [39].

In this case, the constitutive relations become

Generalized strain-displacement relations are given by

Equilibrium equations are

where a comma represents differentiation with respect to; the repeated indices denote that summation and are the phonon field stress and strain; and are the phason field stress and strain; and are the electric displacement and the electric field; , , and are the phonon field elastic, phason field elastic, and phonon-phason field coupling coefficients; and are the piezoelectric constants; is the dielectric constant.

Substituting equations (1) and (2) to (3), the governing equation becomes

in which

In general,

Equation (4) can be written as

3. Description of the Problem and Its Solutions

3.1. The Electroelastic field of Parallel Screw Dislocations

The physical problem is shown in Figure 1. There are parallel screw dislocations located at with the Burgers vectors in the vicinity of a wedge-shaped crack in the 1D hexagonal QC piezoelectric material which is under the far-field uniform loading . The dislocations are subjected to the phonon field line force , phason field line force , and line charge at the core. The wedge-shaped crack is infinitely long along the negative real axis, and denotes the wedge angle.

Figure 1: The interaction between a wedge-shaped crack and many parallel screw dislocations.

As shown in Figure 1, we have the following crack boundary conditions:

Proposing a conformal transformation function

which maps the wedge-shaped crack in the z plane to a semi-infinite crack along the negative real axis in the plane (see Figure 2). The problem turns into the interaction of parallel screw dislocations with a semi-infinite crack.

Figure 2: Conformal transformation.

In the plane, we first consider the analytical solution of the elastic field of a screw dislocation. Assuming that one of the dislocations is situated at the point with the Burgers vector , can be expressed by the imaginary part of an arbitrary analytic function following the theory of a complex variable function, and and

It is reasonable to assume that the potential function is composed of two parts; that is,

where is the complex potential for the dislocation in an infinite medium and is the additional term which assumes that the stress boundary conditions on the cleavage surfaces are satisfied. Based on the force conditions of the dislocation, can be chosen as

with

Then,

From equations (10)–(14), we have

Introducing the generalized strains and and the generalized stresses and , the generalized strain field and the stress field can be written as

According to equation (17), the generalized stress stimulated by dislocation which is working on is

Now, we study the electroelastic solution of parallel screw dislocations. Taking the form of the corresponding matrix as , the generalized stress generated by dislocation exerted on the arbitrary point is

3.2. The Interaction of Parallel Screw Dislocations with Crack

It can be seen from the above discussion that the distribution force induced by dislocations on a semi-infinite crack surface is

Here, the condition of the free crack surface is considered; that is,

where Re denotes the real part of an analytical function vector.

The physical meaning of the boundary condition equation (21) is that the additional distribution force represents the cancelling force exerted on the boundary equal to the negative of . With the aid of equation (20), the boundary condition equation (21) leads to a system of the Riemann–Hilbert problem as

Using the standard methods given by Muskhelishvili [40], the general solution of equation (21) arrives at

where

and we have considered the single-valued condition of generalized displacements, i.e., .

The solution of equation (20) together with equations (23) and (24) gives

Considering the stress field of the crack due to the far-field uniform loadings, equation (20) combined with equation (25) yields

in which and

Therefore, the analytic solution of electroelastic interaction of many parallel screw dislocations with the crack is

If and , equation (28) reduces to the result in ref. [41] under some restrictions.

4. Field Intensity Factors and Energy Release Rate

4.1. Field Intensity Factors

In the plane, according to ref. [31], the phonon field stress, the phason field stress, and the electric displacement intensity factor at the wedge-shaped crack tip are defined as

Actually, only , , , , , and are singular at the crack tip. Substituting equation (28) to (29) and setting , we obtain

Equations (30)–(32) agree with Chen’s results [31] if there is only one dislocation and the phason field and the applied generalized stress are not considered. If the applied loading is and the matrix only contains one dislocation whose core are not under the line force, the above equations can reduce to Li’s results [3]. On the basis of the above restrictions, if the phason field is neglected, the solutions agree with those in ref. [34]. Moreover, equations (30)–(32) suggest that is independent of and , is independent of p and , and is independent of and , if all , , , , , and are supposed to be positive.

4.2. Energy Release Rate

The energy release rate can be calculated as [42]

where .

Thus, the energy release rate for the interaction of the wedge-shaped crack with screw dislocations is derived as

5. Image Force

The image force is expressed as [43]

where stresses and strains with superscript “” are the perturbation stresses and strains at the dislocation point.

Therefore, the image force exerted on the dislocation at ( is excluded from ) in the presence of a wedge-shaped crack and a collection of other dislocations at can be expressed as

where

Noting that the above solution can be reduced to Lin and Thomson [35] in the case of and the material is pure elastic. Furthermore, equation (36) is comprised of three terms. The first term containing the quadratic of coming from the image dislocation. The second term containing and which are the interaction forces from the other dislocations acting on the dislocation . The last one containing which is induced by the applied generalized stress.

When all the dislocations are located at the real axis, the solution can be simplified as

i.e.,

6. Numerical Example

Although the field intensity factors, energy release rate, and the image force acting on the dislocation are expressed explicitly in the previous sections, it is not straightforward that how the dislocation location, dislocation distribution, the wedge angle, and the number of dislocations affect them. The numerical examples are therefore performed in this section. The material constants we choose are listed in Table 1 [14].

Table 1: Material constants in the 1D hexagonal piezoelectric quasicrystal.

The Burgers vectors, line forces, and line charge at the dislocation core are

The critical energy release rate is taken as and the wedge angle in the following numerical discussions is . For the sake of discussion, the magnitude of all Burgers vectors is regarded as and all the dislocation have the same value of .

Figures 3–5, respectively, show the phonon field stress, phason field stress, and electric displacement component contours due to the dislocations. When the dislocations are located at the points and , one may note that as the lines of , , and get closer to the crack surface, they become parallel to the surface and the generalized stress becomes smaller. The phonon field stress vanishes when it reaches the crack surface, as well as the phason field stress and electric displacement , so that the boundary condition is satisfied. The lines of , , and are antisymmetric about the straight line , but the lines of , , and are symmetric about the straight line . These results also reveal that the phonon field stress, phason field stress, and the electric displacement field show square root singularity at the crack tip while having one order singularity at the location of dislocation. The phonon field stress contours have the same trend as the phason field stress contours but is opposite to the electric field. In addition, the magnitude of the electric displacement is less than the magnitude of the phason field stress, but the magnitude of the phason field stress is less than the magnitude of the phonon field stress.

(a)

(b)

(c)

(d)

Figure 3: Contours of phonon field stress components are induced by dislocations situated at and , (a) and (b) , and at and , (c) and (d) .

(a)

(b)

(c)

(d)

Figure 4: Contours of phason field stress components are induced by dislocations situated at and , (a) and (b) , and at and , (c) and (d) .

(a)

(b)

(c)

(d)

Figure 5: Contours of electric displacement components are induced by dislocations situated at and , (a) and (b) , and at and , (c) and (d) .

Supposing there are already two dislocations located at and , respectively, the relationship between the field intensity factors and the positions of the other dislocation is exhibited in Figures 6–8 without remote uniform loading. Let us define the normalized phonon field intensity factor, phason field intensity factor, and electric displacement intensity factor as , , and , respectively. It can be seen that the value of and is always negative, but the value of is always positive. We find that the effects of the dislocation position on the three intensity factors are different. Figures 6(a) and Figure 6(b) show that the shielding effect in the phonon field quickly decreases with the value of the normalized distance increasing, while decreases with the value of increasing. From Figures 7(a) and Figure 7(b), it is interesting to note that increases with the increasing , while decreases with the increasing . When the dislocation is located at the real axis, the stress intensity factor is positive. Figures 8(a) and Figure 8(b) indicate that the shielding effect on the electric field decreases as increases, while it decreases as decreases. Taking the stress intensity factor of the phonon field as an example, Figure 9(a) illustrates that increases with , and the larger the wedge angle is, the stronger the shielding effect will become. The variations in the three field intensity factors with the wedge angle are plotted in Figure 9(b). When the two dislocations are located at and , it can be found that the phonon field stress intensity factor increases as the wedge angle increases, but has little effect on the phason field stress intensity factor and the electric displacement intensity factor since the material constants related to the phonon field are far greater than those related to the phason field and the electric field.

(a)

(b)

Figure 6: (a) varies with . (b) varies with .

(a)

(b)

Figure 7: (a) varies with . (b) varies with .

(a)

(b)

Figure 8: (a) varies with . (b) varies with .

(a)

(b)

Figure 9: (a) varies with of different . (b) Field intensity factors vary with .

Similarity, there exist two dislocations located at and , respectively. In this case, Figures 10(a) and 10(b) more clearly show the retarding effects, respectively, of and on the energy release rate without far-field loading. The normalized energy release rate increases as d θ decreases. In contrast, increases as increases. Figure 10(a) implies that contributes none to the energy release rate when its value increases to a certain extent. Figure 10(b) illustrates that does not affect the energy rate when all the dislocations are situated at the real axis. Figures 11(a) and Figure 11(b) depict the effect of the wedge angle on the energy release rate when the dislocations are in different positions. The energy release rate is positive whether the dislocations are located at the real axis or symmetrically distributed at the imaginary axis, which means the wedge angle accelerates the crack propagation. increases as the number of dislocations increases. Also, increases with the increasing when all the dislocations are located at the real axis, while it increases first and then decreases when the dislocations are situated at the imaginary axis symmetrically.

(a)

(b)

Figure 10: (a) varies with . (b) varies with .

(a)

(b)

Figure 11: (a) Energy release rate varies with when all the dislocations are located at the real axis. (b) Energy release rate varies with when all the dislocations are located at the imaginary axis symmetrically.

Making all dislocations located at the real axis, and and , in the absence of the remote applied generalized stress, Figure 12(a) shows that is always negative when the dislocation is between the crack tip and dislocation , which means dislocation is always attracted. When the dislocation from moves to , it is first rejected and then attracted, and the attractive force is gradually increased and the repulsive force is gradually reduced. When the dislocation is farther away from , it is always rejected and the repulsive force is reduced. Figure 12(b) implies that the force on the dislocation is always positive and the magnitude decreases when increasing the wedge angle.

(a)

(b)

Figure 12: (a) The force acting on the dislocation varies with the distance from the dislocation to the crack tip when the other dislocations are located at and . (b) The force acting on the dislocation varies with when the other dislocations are located at and.

7. Conclusions

In this paper, the fracture mechanic model of the two-dimensional crack and parallel dislocations is established. The interaction between parallel screw dislocations and a wedge-shaped crack in a one-dimensional hexagonal quasicrystal with piezoelectric effects is studied, and the closed-form solutions of the electroelastic field, field intensity factors, energy release rate, and generalized Peach-Koehler force are obtained. The analytical solutions can reduce to the results from the existing literature. The disadvantage of the present solution is that the dislocations are considered to distribute uniformly in the numerical examples which do not reflect the microstructure of the material. The solutions are limited in numbers of dislocations, and only the same sign dislocations are considered here. From the numerical illustrations, some important results are summarized in the following outline:

(1)The trend of contours of the phonon field stress is similar to the phason field stress but is different from the electric displacement.(2)The stress concentration occurs at the crack tip and dislocations cores. The effects of the location of dislocations on the three field intensity factors are different, and the wedge angle always increases the shielding effect.(3)The increasing number of dislocations always accelerates the crack propagation. The wedge angle accelerates the crack growth when dislocations are located at the real axis, but it first promotes crack growth and then demotes crack propagation when dislocations are located at the imaginary axis symmetrically.(4)The image force acting on the dislocation in the presence of a crack and other dislocations decreases when increasing the wedge angle. The dislocation is either attracted or rejected in the process of moving from left to right along the real axis.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This project was supported by the National Natural Key R&D Program of China (Grant no. 2017YFC1405605) and the Natural Science Foundation of Inner Mongolia Autonomous Region, China (Grant no. 2018MS01005).

References

D. Shechtman, I. Blech, D. Gratias, and J. W. Cahn, “Metallic phase with long-range orientational order and no translational symmetry,” Physical Review Letters, vol. 53, no. 20, pp. 1951–1953, 1984.

Z. Zhang and K. Urban, “Transmission electron microscope observations of dislocations and stacking faults in a decagonal Al-Cu-Co alloy,” Philosophical Magazine Letters, vol. 60, no. 3, pp. 97–102, 1989.

L. H. Li and G. T. Liu, “A screw dislocation interacting with a wedge-shaped cracks in one-dimensional hexagonal quasicrystal,” Acta Physica Sinica, vol. 61, no. 8, Article ID 086103, 2012.

L. H. Li and G. T. Liu, “The interaction between the dislocation and elliptic notch in one-dimensional hexagonal quasi-crystals,” Modern Physics Letters B, vol. 23, no. 28, pp. 3397–3407, 2009.

Y. Gao, B. X. Xu, B. S. Zhao, and T. C. Chang, “New general solutions of plane elasticity of one- dimensional quasicrystal,” Physica Status Solidi (B), vol. 245, no. 1, pp. 20–27, 2008.

X. F. Li and T. Y. Fan, “Elastic analysis of a mode II crack in a decagonal quasi-crystal,” Chinese Physics, vol. 11, no. 3, pp. 266–271, 2002.

J.-H. Guo and Z.-X. Lu, “Exact solution of four cracks originating from an elliptical hole in one-dimensional hexagonal quasicrystals,” Applied Mathematics and Computation, vol. 217, no. 22, pp. 9397–9403, 2011.

J. H. Guo and G. T. Liu, “Exact analytic solutions for an elliptic hole with asymmetric collinear cracks in a one-dimensional hexagonal quasi-crystal,” Chinese Physics B, vol. 17, no. 7, pp. 2610–2620, 2008.

G. T. Liu, R. P. Guo, and T. Y. Fan, “On the interaction between dislocations and cracks in one- dimensional hexagonal quasi-crystals,” Chinese Physics, vol. 12, no. 10, pp. 1149–1155, 2003.

F. Lou, T. Cao, T. Qin, and C. Xu, “Plane analysis for an inclusion in 1D hexagonal quasicrystal using the hypersingular integral equation method,” Acta Mechanica Solida Sinica, vol. 32, no. 2, pp. 249–260, 2019.

X. Wang and P. Schiavone, “Elastic field near the tip of an anticrack in a decagonal quasicrys-talline material,” Applied Mathematics and Mechanics, vol. 41, no. 3, pp. 401–408, 2020.

X. Zhao, X. Li, and S. Ding, “Two kinds of contact problems in three-dimensional icosahedral quasicrystals,” Applied Mathematics and Mechanics, vol. 36, no. 12, pp. 1569–1580, 2015.

Y. Li, M. Zhao, Q.-H. Qin, and C. Fan, “Analysis solution method for 3D planar crack problems of two-dimensional hexagonal quasicrystals with thermal effects,” Applied Mathematical Modelling, vol. 69, pp. 648–664, 2019.

X. Y. Li, P. D. Li, T. H. Wu, M. X. Shi, and Z. W. Zhu, “Three-dimensional fundamental solutions for one-dimensional hexagonal quasicrystal with piezoelectric effect,” Physics Letters A, vol. 378, no. 10, pp. 826–834, 2014.

J. Yu, J. Guo, E. Pan, and Y. Xing, “General solutions of plane problem in one-dimensional quasicrystal piezoelectric materials and its application on fracture mechanics,” Applied Mathematics and Mechanics, vol. 36, no. 6, pp. 793–814, 2015.

J. Yu, J. Guo, and Y. Xing, “Complex variable method for an anti-plane elliptical cavity of one-dimensional hexagonal piezoelectric quasicrystals,” Chinese Journal of Aeronautics, vol. 28, no. 4, pp. 1287–1295, 2015.

J. Yang, Y.-T. Zhou, H.-L. Ma, S.-H. Ding, and X. Li, “The fracture behavior of two asymmetrical limited permeable cracks emanating from an elliptical hole in one-dimensional hexagonal quasicrystals with piezoelectric effect,” International Journal of Solids and Structures, vol. 108, pp. 175–185, 2017.

L.-Z. Yang, Y. Gao, E. Pan, and N. Waksmanski, “An exact solution for a multilayered two-dimensional decagonal quasicrystal plate,” International Journal of Solids and Structures, vol. 51, no. 9, pp. 1737–1749, 2014.

J. Guo, Z. Zhang, and Y. Xing, “Antiplane analysis for an elliptical inclusion in 1D hexagonal piezoelectric quasicrystal composites,” Philosophical Magazine, vol. 96, no. 4, pp. 349–369, 2016.

Y. Li, L. Yang, Y. Gao, and E. Pan, “Cylindrical bending analysis of a layered two-dimensional piezoelectric quasicrystal nanoplate,” Journal of Intelligent Material Systems and Structures, vol. 29, no. 12, pp. 2660–2676, 2018.

L. Zhang, J. Guo, and Y. Xing, “Nonlocal analytical solution of functionally graded multilayered one-dimensional hexagonal piezoelectric quasicrystal nanoplates,” Acta Mechanica, vol. 230, no. 5, pp. 1781–1810, 2019.

Y.-B. Zhou and X.-F. Li, “Two collinear mode-III cracks in one-dimensional hexagonal piezoelectric quasicrystal strip,” Engineering Fracture Mechanics, vol. 189, no. 15, pp. 133–147, 2018.

Y.-B. Zhou and X.-F. Li, “Exact solution of two collinear cracks normal to the boundaries of a 1D layered hexagonal piezoelectric quasicrystal,” Philosophical Magazine, vol. 98, no. 19, pp. 1780–1798, 2018.

Y.-B. Zhou and X.-F. Li, “Fracture analysis of an infinite 1D hexagonal piezoelectric quasicrystal plate with a penny-shaped dielectric crack,” European Journal of Mechanics—A/Solids, vol. 76, pp. 224–234, 2019.

Y.-B. Zhou and X.-F. Li, “A Yoffe-type moving crack in one-dimensional hexagonal piezoelectric quasicrystals,” Applied Mathematical Modelling, vol. 65, pp. 148–163, 2019.

L. Li and G. Liu, “Study on a straight dislocation in an icosahedral quasicrystal with piezoelectric effects,” Applied Mathematics and Mechanics, vol. 39, no. 9, pp. 1259–1266, 2018.

L. H. Li, X. W. Cui, and J. H. Guo, “Interaction between a screw dislocation and an elliptical hole with two asymmetrical cracks in a one-dimensional hexagonal quasicrystal with piezoelectric effect,” Applied Mathematics and Mechanics, 2020.

K. Hu, H. Jin, Z. Yang, and X. Chen, “Interface crack between dissimilar one-dimensional hexagonal quasicrystals with piezoelectric effect,” Acta Mechanica, vol. 230, no. 7, pp. 2455–2474, 2019.

S. M. Ohr, S. J. Chang, and R. Thomson, “Elastic interaction of a wedge crack with a screw dislocation,” Journal of Applied Physics, vol. 57, no. 6, pp. 1839–1843, 1998.

T. Y. Zhang, P. Tong, H. Ouyang, and S. Lee, “Interaction of an edge dislocation with a wedge crack,” Journal of Applied Physics, vol. 78, no. 8, pp. 4873–4880, 1995.

B. J. Chen, Z. M. Xiao, and K. M. Liew, “A screw dislocation in a piezoelectric bi-material wedge,” International Journal of Engineering Science, vol. 40, no. 15, pp. 1665–1685, 2002.

B. J. Chen, Z. M. Xiao, and K. M. Liew, “A piezoelectric screw dislocation near a wedge-shaped bi-material interface,” International Journal of Solids and Structures, vol. 40, no. 9, pp. 2041–2056, 2003.

L. H. Li and Y. Zhao, “Interaction of a screw dislocation with interface and wedge-shaped cracks in one-dimensional hexagonal piezoelectric quasicrystal biomaterial,” Mathematical Problems in Engineering, vol. 2019, Article ID 1037297, 7 pages, 2019.

L. J. Jiang and G. T. Liu, “The interaction between a screw dislocation and wedge-shaped crack in one-dimensional hexagonal piezoelectric quasicrystal,” Chinese Physics B, vol. 26, no. 4, Article ID 044601, 2017.

I.-H. Lin and R. Thomson, “Cleavage, dislocation emission, and shielding for cracks under general loading,” Acta Metallurgica, vol. 34, no. 2, pp. 187–206, 1986.

R. R. Juang and S. Lee, “Elastic interaction between general parallel screw dislocations and a surface crack,” Journal of Applied Physics, vol. 59, no. 10, pp. 3421–3429, 1986.

G. T. Liu and L. Y. Yang, “Interaction between infinitely many dislocations and a semi-infinite crack in one-dimensional hexagonal quasicrystal,” Chinese Physics B, vol. 26, no. 9, Article ID 094601, 2017.

X. Lv and G. T. Liu, “Interaction between many parallel screw dislocations and a semi-infinite crack in a magnetoelectroelastic solid,” Chinese Physics B, vol. 27, no. 7, Article ID 074601, 2018.

G. Altay and M. C. Dökmeci, “On the fundamental equations of piezoelasticity of quasicrystal media,” International Journal of Solids and Structures, vol. 49, no. 23-24, pp. 3255–3262, 2012.

N. I. Muskhelishvili, Singular Integral Equations, Noordhoff, Groningen, Holland, 1953.

R. J. Hao and J. X. Liu, “Interaction of a screw dislocation with a semi-infinite interfacial crack in a magneto-electro-elastic bi-material,” Mechanics Research Communications, vol. 33, no. 3, pp. 415–424, 2006.

K. Y. Lee, W. G. Lee, and Y. E. Pak, “Interaction between a semi-infinite crack and a screw dislocation in a piezoelectric material,” Journal of Applied Mechanics, vol. 67, no. 1, pp. 165–170, 2000.

Y. E. Pak, “Crack extension force in a piezoelectric material,” Journal of Applied Mechanics, vol. 57, no. 3, pp. 647–653, 1990.