A Closed-Form Solution to the Mechanism of Interface Crack Formation with One Contact Area in Decagonal Quasicrystal Bi-Materials

Abstract

Cracks and crack-like defects in engineering structures have greatly reduced the structural strength. An interface crack with one contact area in a combined tension–shear field of decagonal quasicrystal bi-material is investigated. Based on the deformation compatibility equation and displacement potential function, the complex representation of stress and displacement is given. Using the mixed boundary conditions, the closed-form expressions for the stresses and the displacement jumps in the phonon field and phason field on the material interface are obtained. The results show that the stress intensity factor at the crack tip is zero for the phason field. The variation in the stress intensity factor and the length of the contact zone in the phonon field is given, and the result is consistent with the properties of the crystal. The design of safe engineering structures and the formulation of reasonable quality acceptance standards may benefit from the theoretical research carried out here.

1. Introduction

Quasicrystals, as a lightweight, high-strength material suitable for medium-temperature operation, are becoming a new functional and structural material. With the help of various special properties of quasicrystal materials, scientists have made new breakthroughs in various research fields, which has also promoted the development of quasicrystal material applications. Different from ordinary crystals [1,2,3], the two-dimensional decagonal quasicrystals are periodic in the direction of the decagonal rotational symmetry axis, while the arrangement of plane atoms perpendicular to the decagonal rotational symmetry axis is quasi-periodic, which leads to the additional elastic degrees of freedom that do not exist in ordinary crystals and increases the complexity of fracture mechanics research. The quasicrystal phase of the electron diffraction pattern with 10 rotationally symmetric axes was found by Bendersky [4]. At the same time, Feng et al. found the decagonal quasicrystal phase in rapidly cooled Al-Fe alloy [5]. In 1989, two types of dislocation in decagonal quasicrystals were confirmed by the comparative analysis of electron diffraction [6]. Based on the multiplication of the basis sites’ groups, Girzhon et al. proposed the model of the reciprocal lattice of decagonal quasicrystal [7]. By atomic resolution high-angle annular dark-field scanning transmission electron microscopy, Ma and He observed the largest decagonal subunits, which expanded to 5.2 nm in a decagonal shape [8].

There is a lot of the literature on defects in decagonal quasicrystals. By using the Eshelby method, the elastic field and energy of a decagonal quasicrystal with a special dislocation line are given [9]. Fan’s research team used the complex solution of classical elastic theory and introduced the displacement potential function and stress potential function to transform the final governing equation of the two-dimensional decagonal quasicrystal plane elastic problem into a quadruple harmonic equation [10]. Using this theory and conformal transformation of a complex function, the elliptical notch problem is solved [11]. Wang and Zhong studied the interaction between a semi-infinite crack and a line dislocation in a decagonal quasicrystal solid using the complex variable method [12]. Li constructed the complex potential theory of two-dimensional decagonal quasicrystals and further developed Muskhelishvili’s complex variable method [13]. Fan et al. studied the interface crack problem of two-dimensional decagonal quasicrystal bi-material using the propagation displacement discontinuity method [14]. Wang and Schiavone investigated the elastic field near the tip of an anti-crack in a homogeneous decagonal quasicrystal material and presented explicit and elegant expressions for the anti-crack contraction force [15]. The plane problem of a two-dimensional decagonal quasicrystal with a rigid circular inclusion under infinite tension and concentrated force is studied by Zhai et al. [16]. Based on the complex representation of stress and displacement of two-dimensional decagonal quasicrystal, the above problem is transformed into the Riemann boundary problem by using the analytic continuation principle of complex variable function. Zhao et al. extended the displacement and temperature discontinuity method to the two-dimensional decagonal quasicrystal coating structure and studied the mechanical behavior of the interface crack under thermal–mechanical loads [17]. The plane elastic problem of two asymmetric edge cracks in a two-dimensional decagonal quasicrystal elliptical hole under far-field tensile stress is considered by Yu [18]. Li et al. established a phase field framework to simulate the macroscopic brittle fracture of quasicrystal materials [19]. In this phase field model, the volume fraction parameter is introduced into the fracture toughness to reflect the phase wall effect for the first time.

The classical interface crack model [20] assumes that the crack is completely open, which leads to oscillating singularity at the crack tips. By assuming that there is a small contact zone near the crack tip, the unreal vibration singularity is eliminated [21]. Using a singular integral equation formulation, Qin and Mai investigated interface cracks with contact zones in thermo-piezoelectric materials [22]. An analytical solution for an interface crack with one contact zone in anisotropic material was studied by Herrmann and Loboda [23]. Kharun and Loboda studied the crack problem at the interface of two isotropic materials under mixed-mode loading [24]. The interface crack is assumed to be fully open, partially closed, friction-free contact zone, and fully closed. The problem is reduced to a homogeneous combined Dirichlet Riemann boundary value problem and solved in closed form. The problem of interface crack with a frictionless contact zone at the right crack tip between two semi-infinite piezoelectric/piezomagnetic spaces is considered by Herrmann et al. [25]. Saikia and Muthu studied the interface crack by using the non-intrinsic cohesive zone model to eliminate the stress singularity at the crack tip [26]. However, to the best of the authors’ knowledge, the problem of interface crack with contact zone at the crack tip in decagonal quasicrystals has not been studied.

In the present study, the interface crack theory in elastic fracture mechanics is extended to the elastic fracture mechanics of decagonal quasicrystal bi-material. The interface crack problem with a contact zone in decagonal quasicrystal bi-material is considered, which has a contact zone penetrating the solid along the quasi-periodic direction. By using the method of complex variable function, the mixed boundary value problem is transformed into the Dirichlet Riemann boundary value problem, and the closed solution of the problem is obtained.

2. Basic Equations

The stress and strain of decagonal quasicrystal satisfy generalized Hooke’s law [10,27,28]:

$$\begin{array}{l}{\sigma }_{xx}=C_{11}{\epsilon }_{xx}+C_{12}{\epsilon }_{yy}+R\left({\omega }_{xx}+{\omega }_{yy}\right)\\ {\sigma }_{yy}=C_{12}{\epsilon }_{xx}+C_{11}{\epsilon }_{yy}-R\left({\omega }_{xx}+{\omega }_{yy}\right)\\ {\sigma }_{xy}={\sigma }_{yx}=\left(C_{11}-C_{12}\right){\epsilon }_{xy}+R\left({\omega }_{yx}-{\omega }_{xy}\right)\\ H_{xx}=K_1{\omega }_{xx}+K_2{\omega }_{yy}+R\left({\epsilon }_{xx}-{\epsilon }_{yy}\right)\\ H_{yy}=K_2{\omega }_{xx}+K_1{\omega }_{yy}+R\left({\epsilon }_{xx}-{\epsilon }_{yy}\right)\\ H_{xy}=K_1{\omega }_{xy}-K_2{\omega }_{yx}-2R{\epsilon }_{xy}\\ H_{yx}=K_1{\omega }_{yx}-K_2{\omega }_{xy}+2R{\epsilon }_{xy}\end{array}$$

(1)

From Equation (1), the strain relation expressed by stress can be written as

$$\begin{array}{l}{\epsilon }_{xx}=\frac{1}{4\left(C_{12}+C_{66}\right)}\left({\sigma }_{xx}+{\sigma }_{yy}\right)+\frac{K_1+K_2}{4r}\left({\sigma }_{xx}-{\sigma }_{yy}\right)-\frac{R}{2r}\left(H_{xx}+H_{yy}\right)\\ {\epsilon }_{yy}=\frac{1}{4\left(C_{12}+C_{66}\right)}\left({\sigma }_{xx}+{\sigma }_{yy}\right)-\frac{K_1+K_2}{4r}\left({\sigma }_{xx}-{\sigma }_{yy}\right)+\frac{R}{2r}\left(H_{xx}+H_{yy}\right)\\ {\epsilon }_{xy}={\epsilon }_{yx}=\frac{K_1+K_2}{2r}{\sigma }_{xy}+\frac{R}{2r}\left(H_{xy}-H_{yx}\right)\end{array}$$

(2)

$$\begin{array}{l}{\omega }_{xx}=\frac{1}{2\left(K_1-K_2\right)}\left(H_{xx}-H_{yy}\right)+\frac{C_{66}}{2r}\left(H_{xx}+H_{yy}\right)-\frac{R}{2r}\left({\sigma }_{xx}-{\sigma }_{yy}\right)\\ {\omega }_{yy}=-\frac{1}{2\left(K_1-K_2\right)}\left(H_{xx}-H_{yy}\right)+\frac{C_{66}}{2r}\left(H_{xx}+H_{yy}\right)-\frac{R}{2r}\left({\sigma }_{xx}-{\sigma }_{yy}\right)\\ {\omega }_{xy}=\frac{R}{r}{\sigma }_{xy}+\frac{1}{2\left(K_1-K_2\right)}\left(H_{xy}+H_{yx}\right)+\frac{C_{66}}{2r}\left(H_{xy}-H_{yx}\right)\\ {\omega }_{yx}=-\frac{R}{r}{\sigma }_{xy}+\frac{1}{2\left(K_1-K_2\right)}\left(H_{xy}+H_{yx}\right)-\frac{C_{66}}{2r}\left(H_{xy}-H_{yx}\right)\end{array}$$

(3)

where ${\sigma }_{ks}$, ${\epsilon }_{ks}$ and $C_{ij}\left(k,s=x,y;i, j=1,2\right)$ are stresses, strains, and elastic constants in the phonon field, respectively; $H_{ks}$, ${\omega }_{ks}$, and $K_i$ are the stresses, strains, and elastic constants in the phason field, respectively; $R$ is the phonon–phason coupling elastic constant.

The strain relation of the plane elastic problem on the quasi-periodic plane of decagonal quasicrystal can be expressed by displacements as

$$\begin{array}{l}{\epsilon }_{xx}=\frac{\partial u_x}{\partial x}, {\epsilon }_{yy}=\frac{\partial u_y}{\partial y}, {\epsilon }_{xy}=\frac{1}{2}\left(\frac{\partial u_x}{\partial y}+\frac{\partial u_y}{\partial x}\right)\\ {\omega }_{xx}=\frac{\partial w_x}{\partial x}, {\omega }_{yy}=\frac{\partial w_y}{\partial y}, {\omega }_{xy}=\frac{\partial w_x}{\partial y}, {\omega }_{yx}=\frac{\partial w_y}{\partial x}\end{array}$$

(4)

After eliminating the displacement, the deformation coordinate equation expressed by strain is

$$\begin{array}{l}\frac{{\partial }^2{\epsilon }_{xx}}{\partial y^2}+\frac{{\partial }^2{\epsilon }_{yy}}{\partial x^2}=2\frac{{\partial }^2{\epsilon }_{xy}}{\partial x\partial y}\\ \frac{\partial {\omega }_{xy}}{\partial x}=\frac{\partial {\omega }_{xx}}{\partial y}\\ \frac{\partial {\omega }_{yx}}{\partial y}=\frac{\partial {\omega }_{yy}}{\partial x}\end{array}$$

(5)

Introducing stress potential functions $\varphi (x,y)$, ${\psi }_1(x,y)$, ${\psi }_2(x,y)$, the stress–strain relationship can be expressed as

$$\begin{array}{l}{\sigma }_{xx}=\frac{{\partial }^2\varphi }{\partial y^2}, {\sigma }_{yy}=\frac{{\partial }^2\varphi }{\partial x^2}, {\sigma }_{xy}={\sigma }_{yx}=-\frac{{\partial }^2\varphi }{\partial x\partial y}, \\ H_{xx}=\frac{\partial {\psi }_1}{\partial y}, H_{xy}=-\frac{\partial {\psi }_1}{\partial x}, H_{yx}=-\frac{\partial {\psi }_2}{\partial y}, H_{yy}=\frac{\partial {\psi }_2}{\partial x}\end{array}$$

(6)

Substituting Equation (5) into Equations (2) and (3), one obtains

$$\begin{array}{l}{\epsilon }_{xx}=\frac{1}{4\left(C_{12}+C_{66}\right)}{\nabla }^2\varphi +\frac{K_1+K_2}{4r}\left(\frac{{\partial }^2\varphi }{\partial y^2}-\frac{{\partial }^2\varphi }{\partial x^2}\right)-\frac{R}{2r}\left(\frac{\partial {\psi }_1}{\partial y}+\frac{\partial {\psi }_2}{\partial x}\right)\\ {\epsilon }_{yy}=\frac{1}{4\left(C_{12}+C_{66}\right)}{\nabla }^2\varphi -\frac{K_1+K_2}{4r}\left(\frac{{\partial }^2\varphi }{\partial y^2}-\frac{{\partial }^2\varphi }{\partial x^2}\right)+\frac{R}{2r}\left(\frac{\partial {\psi }_1}{\partial y}+\frac{\partial {\psi }_2}{\partial x}\right)\\ {\epsilon }_{xy}={\epsilon }_{yx}=-\frac{K_1+K_2}{2r}\frac{{\partial }^2\varphi }{\partial x\partial y}+\frac{R}{2r}\left(\frac{\partial {\psi }_2}{\partial y}-\frac{\partial {\psi }_1}{\partial x}\right)\end{array}$$

(7)

$$\begin{array}{l}{\omega }_{xx}=\frac{1}{2\left(K_1-K_2\right)}\left(\frac{\partial {\psi }_1}{\partial y}-\frac{\partial {\psi }_2}{\partial x}\right)+\frac{C_{66}}{2r}\left(\frac{\partial {\psi }_1}{\partial y}+\frac{\partial {\psi }_2}{\partial x}\right)-\frac{R}{2r}\left(\frac{{\partial }^2\varphi }{\partial y^2}-\frac{{\partial }^2\varphi }{\partial x^2}\right)\\ {\omega }_{yy}=-\frac{1}{2\left(K_1-K_2\right)}\left(\frac{\partial {\psi }_1}{\partial y}-\frac{\partial {\psi }_2}{\partial x}\right)+\frac{C_{66}}{2r}\left(\frac{\partial {\psi }_1}{\partial y}+\frac{\partial {\psi }_2}{\partial x}\right)-\frac{R}{2r}\left(\frac{{\partial }^2\varphi }{\partial y^2}-\frac{{\partial }^2\varphi }{\partial x^2}\right)\\ {\omega }_{xy}=-\frac{R}{r}\frac{{\partial }^2\varphi }{\partial x\partial y}-\frac{1}{2\left(K_1-K_2\right)}\left(\frac{\partial {\psi }_1}{\partial x}+\frac{\partial {\psi }_2}{\partial y}\right)+\frac{C_{66}}{2r}\left(\frac{\partial {\psi }_2}{\partial y}-\frac{\partial {\psi }_1}{\partial x}\right)\\ {\omega }_{yx}=\frac{R}{r}\frac{{\partial }^2\varphi }{\partial x\partial y}-\frac{1}{2\left(K_1-K_2\right)}\left(\frac{\partial {\psi }_1}{\partial x}+\frac{\partial {\psi }_2}{\partial y}\right)+\frac{C_{66}}{2r}\left(\frac{\partial {\psi }_1}{\partial x}-\frac{\partial {\psi }_2}{\partial y}\right)\end{array}$$

(8)

Substituting Equations (7) and (8) into Equation (5) yields

$$\begin{array}{l}\left(\frac{1}{2\left(C_{12}+C_{66}\right)}+\frac{K_1+K_2}{2r}\right){\nabla }^2{\nabla }^2\varphi +\frac{R}{r}\left(\frac{\partial }{\partial y}{\Pi }_1{\psi }_1-\frac{\partial }{\partial x}{\Pi }_2{\psi }_2\right)=0\\ \left(\frac{r}{K_1-K_2}+C_{66}\right){\nabla }^2{\psi }_1+R\frac{\partial }{\partial y}{\Pi }_1\varphi =0\\ \left(\frac{r}{K_1-K_2}+C_{66}\right){\nabla }^2{\psi }_2-R\frac{\partial }{\partial x}{\Pi }_2\varphi =0\end{array}$$

(9)

where

${\nabla }^2=\frac{{\partial }^2}{\partial x^2}+\frac{{\partial }^2}{\partial y^2}, {\Pi }_1=3\frac{{\partial }^2}{\partial x^2}-\frac{{\partial }^2}{\partial y^2}, {\Pi }_2=3\frac{{\partial }^2}{\partial y^2}-\frac{{\partial }^2}{\partial x^2},$$r=C_{66}\left(K_1+K_2\right)-2R^2$.

Introducing a new function $G(x,y)$, one obtains

$$\varphi =r_1{\nabla }^2{\nabla }^{2}G,{\psi }_1=-R\frac{\partial }{\partial y}{\Pi }_1{\nabla }^{2}G,{\psi }_2=R\frac{\partial }{\partial x}{\Pi }_2{\nabla }^{2}G$$

(10)

where $r_1=\frac{r}{K_1-K_2}+C_{66}$.

Equation (9) is automatically satisfied, so $G$ is called the stress function. One has [29]

$${\nabla }^2{\nabla }^2{\nabla }^2{\nabla }^{2}G=0$$

(11)

Therefore, the final control function based on stress potential is a quadruple harmonic equation. Substituting Equation (10) into Equation (6), one obtains

$$\begin{array}{l}{\sigma }_{xx}=r_1\frac{{\partial }^2}{\partial y^2}{\nabla }^2{\nabla }^{2}G, {\sigma }_{yy}=r_1\frac{{\partial }^2}{\partial x^2}{\nabla }^2{\nabla }^{2}G, \\ {\sigma }_{xy}={\sigma }_{yx}=-r_1\frac{{\partial }^2}{\partial x\partial y}{\nabla }^2{\nabla }^{2}G,\end{array}$$

(12)

$$\begin{array}{l}{H}_{xx}=-R\frac{{\partial }^2}{\partial y^2}{\Pi }_1{\nabla }^{2}G, H_{xy}=R\frac{{\partial }^2}{\partial x\partial y}{\Pi }_1{\nabla }^{2}G, \\ H_{yx}=-R\frac{{\partial }^2}{\partial x\partial y}{\Pi }_2{\nabla }^{2}G, H_{yy}=R\frac{{\partial }^2}{\partial x^2}{\Pi }_2{\nabla }^{2}G\end{array}$$

(13)

Based on the method of stress potential function, Li and Fan [30] developed the complex variable function solution of the quartic harmonic equation. The fundamental solution of Equation (11) is

$$G=2\mathrm{Re}\left(g_1\left(z\right)+\overline{z}{g}_2\left(z\right)+\frac{1}{2}{\overline{z}}^2{g}_3\left(z\right)+\frac{1}{6}{\overline{z}}^3{g}_4\left(z\right)\right)$$

(14)

in which $g_j(z)(j=1,2,3,4)$ is about four analytic functions of complex variables; $z=x+iy$, and $i=\sqrt{-1}$. Superscript “$-$” indicates complex conjugation, i.e. $z=x-iy$.

Substituting Equation (14) into Equations (12) and (13), the complex expression of stress function is obtained as

$$\begin{array}{l}{\sigma }_{xx}=-32r_1\mathrm{Re}\left(g_3^{\left(4\right)}\left(z\right)+\overline{z}{g}_4^{\left(4\right)}\left(z\right)-2g_4^{‴}\left(z\right)\right)\\ {\sigma }_{yy}=32r_1\mathrm{Re}\left(g_3^{\left(4\right)}\left(z\right)+\overline{z}{g}_4^{\left(4\right)}\left(z\right)+2g_4^{‴}\left(z\right)\right)\\ {\sigma }_{xy}={\sigma }_{yx}=32r_1\mathrm{Im}\left(g_3^{\left(4\right)}\left(z\right)+\overline{z}{g}_4^{\left(4\right)}\left(z\right)\right)\end{array}$$

(15)

$$\begin{array}{l}{H}_{xx}=32R\mathrm{Re}\left(g_2^{\left(5\right)}\left(z\right)+\overline{z}{g}_3^{\left(5\right)}\left(z\right)+\frac{1}{2}{\overline{z}}^2{g}_4^{\left(5\right)}\left(z\right)-g_3^{\left(4\right)}\left(z\right)-\overline{z}{g}_4^{\left(4\right)}\left(z\right)\right)\\ H_{yy}=-32R\mathrm{Re}\left(g_2^{\left(5\right)}\left(z\right)+\overline{z}{g}_3^{\left(5\right)}\left(z\right)+\frac{1}{2}{\overline{z}}^2{g}_4^{\left(5\right)}\left(z\right)+g_3^{\left(4\right)}\left(z\right)+\overline{z}{g}_4^{\left(4\right)}\left(z\right)\right)\\ H_{xy}=-32R\mathrm{Im}\left(g_2^{\left(5\right)}\left(z\right)+\overline{z}{g}_3^{\left(5\right)}\left(z\right)+\frac{1}{2}{\overline{z}}^2{g}_4^{\left(5\right)}\left(z\right)+g_3^{\left(4\right)}\left(z\right)+\overline{z}{g}_4^{\left(4\right)}\left(z\right)\right)\\ H_{yx}=-32R\mathrm{Im}\left(g_2^{\left(5\right)}\left(z\right)+\overline{z}{g}_3^{\left(5\right)}\left(z\right)+\frac{1}{2}{\overline{z}}^2{g}_4^{\left(5\right)}\left(z\right)-g_3^{\left(4\right)}\left(z\right)-\overline{z}{g}_4^{\left(4\right)}\left(z\right)\right)\end{array}$$

(16)

New analytic functions $f_j(z)(j=2,3,4)$ are introduced to make

$$f_2\left(z\right)=g_2^{\left(4\right)}\left(z\right), f_3\left(z\right)=g_3^{‴}\left(z\right), f_4\left(z\right)=g_4^{″}\left(z\right)$$

(17)

Then, Equations (15) and (16) can be rewritten as

$$\begin{array}{l}{\sigma }_{xx}=-32r_1\mathrm{Re}\left(f_3^{\prime }\left(z\right)+\overline{z}{f}_4^{″}\left(z\right)-2f_4^{\prime }\left(z\right)\right)\\ {\sigma }_{yy}=32r_1\mathrm{Re}\left(f_3^{\prime }\left(z\right)+\overline{z}{f}_4^{″}\left(z\right)+2f_4^{\prime }\left(z\right)\right)\\ {\sigma }_{xy}={\sigma }_{yx}=32r_1\mathrm{Im}\left(f_3^{\prime }\left(z\right)+\overline{z}{f}_4^{″}\left(z\right)\right)\end{array}$$

(18)

$$\begin{array}{l}{H}_{xx}=32R\mathrm{Re}\left(f_2^{\prime }\left(z\right)+\overline{z}{f}_3^{″}\left(z\right)+\frac{1}{2}{\overline{z}}^2{f}_4^{‴}\left(z\right)-f_3^{\prime }\left(z\right)-\overline{z}{f}_4^{″}\left(z\right)\right)\\ H_{yy}=-32R\mathrm{Re}\left(f_2^{\prime }\left(z\right)+\overline{z}{f}_3^{″}\left(z\right)+\frac{1}{2}{\overline{z}}^2{f}_4^{‴}\left(z\right)+f_3^{\prime }\left(z\right)+\overline{z}{f}_4^{″}\left(z\right)\right)\\ H_{xy}=-32R\mathrm{Im}\left(f_2^{\prime }\left(z\right)+\overline{z}{f}_3^{″}\left(z\right)+\frac{1}{2}{\overline{z}}^2{f}_4^{‴}\left(z\right)+f_3^{\prime }\left(z\right)+\overline{z}{f}_4^{″}\left(z\right)\right)\\ H_{yx}=-32R\mathrm{Im}\left(f_2^{\prime }\left(z\right)+\overline{z}{f}_3^{″}\left(z\right)+\frac{1}{2}{\overline{z}}^2{f}_4^{‴}\left(z\right)-f_3^{\prime }\left(z\right)-\overline{z}{f}_4^{″}\left(z\right)\right)\end{array}$$

(19)

3. Statement of the Problem

Consider an interface crack with one contact area between two bonded semi-infinite decagonal quasicrystals. $\sigma$ and $\tau$ are uniformly loaded in the phonon field, and ${\sigma }_{xx}^{\left({\rm I}\right)\infty }$, ${\sigma }_{xx}^{\left({\rm I}{\rm I}\right)\infty }$, $H_{xx}^{\left({\rm I}\right)\infty }$, $H_{xx}^{\left({\rm I}{\rm I}\right)\infty }$ are applied at infinity, as shown in Figure 1.

The continuity conditions on the interface can be written as follows

$$\left\{\begin{array}{l}\left({\sigma }_{yy}^{\left(\mathrm{I}\right)}\left(x,0\right)-\mathrm{i}{\sigma }_{xy}^{\left(\mathrm{I}\right)}\left(x,0\right)\right)-\left({\sigma }_{yy}^{\left(\mathrm{I}\mathrm{I}\right)}\left(x,0\right)-\mathrm{i}{\sigma }_{xy}^{\left(\mathrm{I}\mathrm{I}\right)}\left(x,0\right)\right)=0\\ \left(H_{yy}^{\left(\mathrm{I}\right)}\left(x,0\right)-\mathrm{i}{H}_{xy}^{\left(\mathrm{I}\right)}\left(x,0\right)\right)-\left(H_{yy}^{\left(\mathrm{I}\mathrm{I}\right)}\left(x,0\right)-\mathrm{i}{H}_{xy}^{\left(\mathrm{I}\mathrm{I}\right)}\left(x,0\right)\right)=0\end{array}\right.\hspace{1em}x\in \left(-\infty ,\infty \right)$$

(20)

$$\left\{\begin{array}{l}\left(u_x^{\left(\mathrm{I}\right)}\left(x,0\right)+\mathrm{i}{u}_y^{\left(\mathrm{I}\right)}\left(x,0\right)\right)-\left(u_x^{\left(\mathrm{I}\mathrm{I}\right)}\left(x,0\right)+\mathrm{i}{u}_y^{\left(\mathrm{I}\mathrm{I}\right)}\left(x,0\right)\right)=0\\ \left(w_x^{\left(\mathrm{I}\right)}\left(x,0\right)+\mathrm{i}{w}_y^{\left(\mathrm{I}\right)}\left(x,0\right)\right)-\left(w_x^{\left(\mathrm{I}\mathrm{I}\right)}\left(x,0\right)+\mathrm{i}{w}_y^{\left(\mathrm{I}\mathrm{I}\right)}\left(x,0\right)\right)=0\end{array}\right.\hspace{1em}x\in \left(-\infty ,c\right)\cup \left(a,\infty \right)$$

(21)

The boundary conditions on the crack face can be expressed as follows:

$$\left\{\begin{array}{l}{\sigma }_{xy}^{\left(\mathrm{I}\right)}\left(x,0\right)=0\\ H_{xy}^{\left(\mathrm{I}\right)}\left(x,0\right)=0\\ u_y^{\left(\mathrm{I}\right)}\left(x,0\right)-u_y^{\left(\mathrm{I}\mathrm{I}\right)}\left(x,0\right)=0\\ w_y^{\left(\mathrm{I}\right)}\left(x,0\right)-w_y^{\left(\mathrm{I}\mathrm{I}\right)}\left(x,0\right)=0\end{array}\right.\hspace{1em}x\in \left(b,a\right)$$

(22)

$$\left\{\begin{array}{l}{\sigma }_{yy}^{\left(\mathrm{I}\right)}\left(x,0\right)-\mathrm{i}{\sigma }_{xy}^{\left(\mathrm{I}\right)}\left(x,0\right)=0\\ H_{yy}^{\left(\mathrm{I}\right)}\left(x,0\right)-\mathrm{i}{H}_{xy}^{\left(\mathrm{I}\right)}\left(x,0\right)=0\end{array}\right.\hspace{1em}x\in \left(c,b\right)$$

(23)

where ${\sigma }_{12}^{(k)}\left(x,y\right)$; ${\sigma }_{22}^{(k)}\left(x,y\right)$, $u_1^{(k)}\left(x,y\right)$, and $u_2^{(k)}\left(x,y\right)$ are the phonon field of shear stresses, normal stresses, and displacements along $x$- and $y$-axes, respectively. $H_{12}^{(k)}\left(x,y\right)$, $H_{22}^{(k)}\left(x,y\right)$, $w_1^{(k)}\left(x,y\right)$, and $w_2^{(k)}\left(x,y\right)$ are the phason field of shear stresses, normal stresses, and displacements along $x$- and $y$-axes, respectively. Subscripts $k={\rm I}$ and $k=\mathrm{I}\mathrm{I}$ mean, respectively, the upper and lower half-planes. Intervals $(-\infty ,c)\cup (a,\infty )$, $\left[c,b\right]$, and $(b,a)$ denote the bond, open part of the crack, and contact zone, respectively.

4. Theoretical Derivation of Interface Stresses and Displacement Jump

For the plane elastic problems, the stresses and displacements in the decagonal quasicrystal of point groups, 10 can be expressed in terms of the complex potentials as follows:

$$\left\{\begin{array}{l}{\sigma }_{yy}^{\left(k\right)}\left(x,y\right)-\mathrm{i}{\sigma }_{xy}^{\left(k\right)}\left(x,y\right)=32r_1^{\left(k\right)}\left(f_{4k}^{\prime }\left(z\right)+\overline{f_{3k}^{\prime }\left(z\right)}+z\overline{f_{4k}^{″}\left(z\right)}+\overline{f_{4k}^{\prime }\left(z\right)}\right)\\ u_x^{\left(k\right)}\left(x,y\right)+\mathrm{i}{u}_y^{\left(k\right)}\left(x,y\right)=32r_6^{\left(k\right)}{f}_{4k}\left(z\right)-32r_5^{\left(k\right)}\left(\overline{f_{3k}\left(z\right)}+z\overline{f_{4k}^{\prime }\left(z\right)}\right)\end{array}\right.$$

(24)

$$\left\{\begin{array}{l}{H}_{yy}^{\left(k\right)}\left(x,y\right)-\mathrm{i}{H}_{xy}^{\left(k\right)}\left(x,y\right)=-32R^{\left(k\right)}\left(\overline{f_{2k}^{\prime }\left(z\right)}+z\overline{f_{3k}^{″}\left(z\right)}+\frac{1}{2}{z}^2\overline{f_{4k}^{‴}\left(z\right)}+\overline{f_{3k}^{\prime }\left(z\right)}+z\overline{f_{4k}^{″}\left(z\right)}\right)\\ w_x^{\left(k\right)}\left(x,y\right)+\mathrm{i}{w}_y^{\left(k\right)}\left(x,y\right)=32r_7^{\left(k\right)}\left(\overline{f_{2k}\left(z\right)}+z\overline{f_{3k}^{\prime }\left(z\right)}+\frac{1}{2}{z}^2\overline{f_{4k}^{″}\left(z\right)}\right)\end{array}\right.$$

(25)

$$\left\{\begin{array}{l}{\sigma }_{xx}^{\left(k\right)}\left(x,y\right)+{\sigma }_{yy}^{\left(k\right)}\left(x,y\right)=64r_1^{\left(k\right)}\left(f_{4k}^{\prime }\left(z\right)+\overline{f_{4k}^{\prime }\left(z\right)}\right)\\ H_{xx}^{\left(k\right)}\left(x,y\right)+H_{yy}^{\left(k\right)}\left(x,y\right)=-32R^{\left(k\right)}\left(f_{3k}^{\prime }\left(z\right)+\overline{z}{f}_{4k}^{″}\left(z\right)+\overline{f_{3k}^{\prime }\left(z\right)}+z\overline{f_{4k}^{″}\left(z\right)}\right)\end{array}\right.$$

(26)

in which

$r^{\left(k\right)}=C_{66}^{\left(k\right)}\left(K_1^{\left(k\right)}+K_2^{\left(k\right)}\right)-2R^{\left(k\right)2}$, $r_1^{\left(k\right)}=\frac{r^{\left(k\right)}}{K_1^{\left(k\right)}-K_2^{\left(k\right)}}+C_{66}^{\left(k\right)}$,

$r_2^{\left(k\right)}=\frac{1}{2\left(C_{12}^{\left(k\right)}+C_{66}^{\left(k\right)}\right)}+\frac{K_1^{\left(k\right)}+K_2^{\left(k\right)}}{2r^{\left(k\right)}}$, $r_3^{\left(k\right)}=\frac{R^{\left(k\right)2}}{r^{\left(k\right)}}$, $r_4^{\left(k\right)}=\frac{K_1^{\left(k\right)}+K_2^{\left(k\right)}}{r^{\left(k\right)}}$,

$r_5^{\left(k\right)}=r_1^{\left(k\right)}{r}_4^{\left(k\right)}-r_3^{\left(k\right)}, r_6^{\left(k\right)}=4r_1^{\left(k\right)}{r}_2^{\left(k\right)}-r_3^{\left(k\right)}-r_1^{\left(k\right)}{r}_4^{\left(k\right)}, r_7^{\left(k\right)}=\frac{R^{\left(k\right)}}{K_1^{\left(k\right)}-K_2^{\left(k\right)}}$

Introducing two new functions

$$s_k\left(z\right)={\overline{f}}_{3k}\left(z\right)+z{\overline{f}}_{4k}^{\prime }\left(z\right)$$

(27)

$$t_k\left(z\right)={\overline{f}}_{2k}\left(z\right)+z{\overline{f}}_{3k}^{\prime }\left(z\right)+\frac{1}{2}{z}^2{\overline{f}}_{4k}^{″}\left(z\right)$$

(28)

Replacing $z$ with $\overline{z}$ in Equations (27) and (28) to obtain the following expressions

$$\overline{f_{3k}\left(z\right)}=s_k\left(\overline{z}\right)-\overline{z}\overline{f_{4k}^{\prime }\left(z\right)}$$

(29)

$$\overline{f_{2k}\left(z\right)}=t_k\left(\overline{z}\right)-\overline{z}\overline{f_{3k}^{\prime }\left(z\right)}-\frac{1}{2}{\overline{z}}^2\overline{f_{4k}^{″}\left(z\right)}$$

(30)

The differential of the above formula is

$$\overline{f_{3k}^{\prime }\left(z\right)}+\overline{f_{4k}^{\prime }\left(z\right)}=s_k^{\prime }\left(\overline{z}\right)-\overline{z}\overline{f_{4k}^{″}\left(z\right)}$$

(31)

$$\overline{f_{3k}^{″}\left(z\right)}+2\overline{f_{4k}^{″}\left(z\right)}=s_k^{″}\left(\overline{z}\right)-\overline{z}\overline{f_{4k}^{‴}\left(z\right)}$$

(32)

$$\overline{f_{2k}^{\prime }\left(z\right)}+\overline{f_{3k}^{\prime }\left(z\right)}=t_k^{\prime }\left(\overline{z}\right)-\overline{z}\overline{f_{4k}^{″}\left(z\right)}-\overline{z}\overline{f_{3k}^{″}\left(z\right)}-\frac{1}{2}{\overline{z}}^2\overline{f_{4k}^{‴}\left(z\right)}$$

(33)

Substituting Equations (27), (29) and (31) into Equation (24) yields

$$\left\{\begin{array}{l}{\sigma }_{yy}^{\left(k\right)}\left(x,y\right)-\mathrm{i}{\sigma }_{xy}^{\left(k\right)}\left(x,y\right)=32r_1^{\left(k\right)}\left({{\rm I}}_k\left(z\right)+{\Gamma }_k\left(\overline{z}\right)+\left(z-\overline{z}\right)\overline{{{\rm I}}_k^{\prime }\left(z\right)}\right)\\ u_x^{\left(k\right)}\left(x,y\right)+\mathrm{i}{u}_y^{\left(k\right)}\left(x,y\right)=32r_6^{\left(k\right)}{f}_{4k}\left(z\right)-32r_5^{\left(k\right)}\left(s_k\left(\overline{z}\right)+\left(z-\overline{z}\right)\overline{{{\rm I}}_k\left(z\right)}\right)\end{array}\right.$$

(34)

Substituting Equations (28), (30) and (32) into Equation (25), one obtains

$$\left\{\begin{array}{l}{H}_{yy}^{\left(k\right)}\left(x,y\right)-\mathrm{i}{H}_{xy}^{\left(k\right)}\left(x,y\right)=-32R^{\left(k\right)}\left(H_k\left(\overline{z}\right)+\left(z-\overline{z}\right)\left({\Gamma }_k^{\prime }\left(\overline{z}\right)-\overline{{{{\rm I}}^{\prime }}_k\left(z\right)}\right)+\frac{1}{2}{\left(z-\overline{z}\right)}^2\overline{{{\rm I}}_k^{″}\left(z\right)}\right)\\ w_x^{\left(k\right)}\left(x,y\right)+\mathrm{i}{w}_y^{\left(k\right)}\left(x,y\right)=32r_7^{\left(k\right)}\left(t_k\left(\overline{z}\right)+\left(z-\overline{z}\right)\overline{f_{3k}^{\prime }\left(z\right)}+\frac{1}{2}\left(z^2-{\overline{z}}^2\right)\overline{f_{4k}^{″}\left(z\right)}\right)\end{array}\right.$$

(35)

Introduce the following notation

$${{\rm I}}_k\left(z\right)=f_{4k}^{\prime }\left(z\right), {\Gamma }_k\left(\overline{z}\right)=s_k^{\prime }\left(\overline{z}\right), {\mathrm{H}}_k\left(\overline{z}\right)=t_k^{\prime }\left(\overline{z}\right)$$

(36)

Thus,

$$\begin{array}{l}\overline{f_{3k}^{\prime }\left(z\right)}={\Gamma }_k\left(\overline{z}\right)-\overline{{{\rm I}}_k\left(z\right)}-\overline{z}\overline{{{\rm I}}_k^{\prime }\left(z\right)},\\ f_{3k}^{\prime }\left(z\right)=\overline{{\Gamma }_k\left(\overline{z}\right)}-{{\rm I}}_k\left(z\right)-z{{\rm I}}_k^{\prime }\left(z\right),\end{array}$$

(37)

$\overline{f_{3k}^{″}\left(z\right)}={\Gamma }_k^{\prime }\left(\overline{z}\right)-2\overline{{{\rm I}}_k^{\prime }\left(z\right)}-\overline{z}\overline{{{\rm I}}_k^{″}\left(z\right)}$

Substituting Equation (37) into Equations (34) and (35), one obtains

$$\left\{\begin{array}{l}{\sigma }_{yy}^{\left(k\right)}\left(x,y\right)-\mathrm{i}{\sigma }_{xy}^{\left(k\right)}\left(x,y\right)=32r_1^{\left(k\right)}\left({{\rm I}}_k\left(z\right)+{\Gamma }_k\left(\overline{z}\right)+\left(z-\overline{z}\right)\overline{{{\rm I}}_k^{\prime }\left(z\right)}\right)\\ \frac{\partial \left(u_x^{\left(k\right)}\left(x,y\right)+\mathrm{i}{u}_y^{\left(k\right)}\left(x,y\right)\right)}{\partial x}=32r_6^{\left(k\right)}{{\rm I}}_k\left(z\right)-32r_5^{\left(k\right)}\left({\Gamma }_k\left(\overline{z}\right)+\left(z-\overline{z}\right)\overline{{{\rm I}}_k^{\prime }\left(z\right)}\right)\end{array}\right.$$

(38)

$$\left\{\begin{array}{l}{H}_{yy}^{\left(k\right)}\left(x,y\right)-\mathrm{i}{H}_{xy}^{\left(k\right)}\left(x,y\right)=-32R^{\left(k\right)}\left({\mathrm{H}}_k\left(\overline{z}\right)+\left(z-\overline{z}\right)\left({\Gamma }_k^{\prime }\left(\overline{z}\right)-\overline{{{{\rm I}}^{\prime }}_k\left(z\right)}\right)+\frac{1}{2}{\left(z-\overline{z}\right)}^2\overline{{{\rm I}}_k^{″}\left(z\right)}\right)\\ \frac{\partial \left(w_x^{\left(k\right)}\left(x,y\right)+\mathrm{i}{w}_y^{\left(k\right)}\left(x,y\right)\right)}{\partial x}=32r_7^{\left(k\right)}\left({\mathrm{H}}_k\left(\overline{z}\right)+\left(z-\overline{z}\right)\left({\Gamma }_k^{\prime }\left(\overline{z}\right)-\overline{{\Gamma }_k^{\prime }\left(z\right)}\right)+\frac{1}{2}{\left(z-\overline{z}\right)}^2\overline{{{\rm I}}_k^{″}\left(z\right)}\right)\end{array}\right.$$

(39)

Equation (26) can be represented as

$$\left\{\begin{array}{l}{\sigma }_{xx}^{\left(k\right)}\left(x,y\right)+{\sigma }_{yy}^{\left(k\right)}\left(x,y\right)=128r_1^{\left(k\right)}\mathrm{Re}{{\rm I}}_k\left(z\right)\\ H_{xx}^{\left(k\right)}\left(x,y\right)+H_{yy}^{\left(k\right)}\left(x,y\right)=-64R^{\left(k\right)}\mathrm{Re}\left(\overline{{\Gamma }_k\left(\overline{z}\right)}-{{\rm I}}_k\left(z\right)+\left(\overline{z}-z\right){{\rm I}}_k^{\prime }\left(z\right)\right)\end{array}\right.$$

(40)

Combined with continuity conditions Equations (20) and (21), the following equations are given:

$$\left\{\begin{array}{l}{r}_1^{\left({\rm I}\right)}{{\rm I}}_{{\rm I}}^{+}\left(x\right)-r_1^{\left(\mathrm{I}\mathrm{I}\right)}{\Gamma }_{\mathrm{I}\mathrm{I}}^{+}\left(x\right)=r_1^{\left(\mathrm{I}\mathrm{I}\right)}{{\rm I}}_{\mathrm{I}\mathrm{I}}^{-}\left(x\right)-r_1^{\left(\mathrm{I}\right)}{\Gamma }_{\mathrm{I}}^{-}\left(x\right)\hspace{1em}x\in \left(-\infty ,\infty \right)\\ r_6^{\left(\mathrm{I}\right)}{{\rm I}}_{\mathrm{I}}^{+}\left(x\right)+r_5^{\left(\mathrm{I}\mathrm{I}\right)}{\Gamma }_{\mathrm{I}\mathrm{I}}^{+}\left(x\right)=r_6^{\left(\mathrm{I}\mathrm{I}\right)}{{\rm I}}_{\mathrm{I}\mathrm{I}}^{-}\left(x\right)+r_5^{\left(\mathrm{I}\right)}{\Gamma }_{\mathrm{I}}^{-}\left(x\right)\hspace{1em}x\in \left(-\infty ,c\right)\cup \left(a,\infty \right)\end{array}\right.$$

(41)

$$\left\{\begin{array}{l}{R}^{\left(\mathrm{I}\mathrm{I}\right)}{\mathrm{H}}_{\mathrm{I}\mathrm{I}}^{+}\left(x\right)=R^{\left(\mathrm{I}\right)}{\mathrm{H}}_{\mathrm{I}}^{-}\left(x\right)\hspace{1em}x\in \left(-\infty ,\infty \right)\\ r_7^{\left({\rm I}{\rm I}\right)}{\mathrm{H}}_{\mathrm{I}\mathrm{I}}^{+}\left(x\right)=r_7^{\left({\rm I}\right)}{\mathrm{H}}_{\mathrm{I}}^{-}\left(x\right) x\in \left(-\infty ,c\right)\cup \left(a,\infty \right)\end{array}\right.$$

(42)

where superscripts “$+$” and “$-$” denote the limit values of the analytical functions, when $z\to x+i0$ and $z\to x-i0$, respectively.

Both sides of Equation (41) represent the boundary values of two analytical functions in their respective half-planes, and the two functions can be analytically extended into the whole plane. For the phonon field, we can introduce the following functions:

$$A_1\left(z\right)=\left\{\begin{array}{l}{r}_1^{\left(\mathrm{I}\right)}{{\rm I}}_{\mathrm{I}}\left(z\right)-r_1^{\left(\mathrm{I}\mathrm{I}\right)}{\Gamma }_{\mathrm{I}\mathrm{I}}\left(z\right)\hspace{1em}y>0\\ r_1^{\left(\mathrm{I}\mathrm{I}\right)}{{\rm I}}_{\mathrm{I}\mathrm{I}}\left(z\right)-r_1^{\left(\mathrm{I}\right)}{\Gamma }_{\mathrm{I}}\left(z\right)\hspace{1em}y<0\end{array}\right.$$

(43)

and

$$A_2\left(z\right)=\left\{\begin{array}{l}{r}_6^{\left(\mathrm{I}\right)}{{\rm I}}_{\mathrm{I}}\left(z\right)+r_5^{\left(\mathrm{I}\mathrm{I}\right)}{\Gamma }_{\mathrm{I}\mathrm{I}}\left(z\right)\hspace{1em}y>0\\ r_6^{\left(\mathrm{I}\mathrm{I}\right)}{{\rm I}}_{\mathrm{I}\mathrm{I}}\left(z\right)+r_5^{\left(\mathrm{I}\right)}{\Gamma }_{\mathrm{I}}\left(z\right)\hspace{1em}y<0\end{array}\right.$$

(44)

Corresponding, for the phason field, introducing new functions $B_1\left(z\right)$ and $B_2\left(z\right)$

$$B_1\left(z\right)=\left\{\begin{array}{l}{R}^{\left({\rm I}{\rm I}\right)}{\mathrm{H}}_{{\rm I}{\rm I}}\left(z\right)\hspace{1em}y>0\\ R^{\left({\rm I}\right)}{\mathrm{H}}_{{\rm I}}\left(z\right)\hspace{1em}y<0\end{array}\right.$$

(45)

and

$$B_2\left(z\right)=\left\{\begin{array}{l}{r}_7^{\left({\rm I}{\rm I}\right)}{\mathrm{H}}_{{\rm I}{\rm I}}\left(z\right)\hspace{1em}y>0\\ r_7^{\left({\rm I}\right)}{\mathrm{H}}_{{\rm I}}\left(z\right)\hspace{1em}y<0\end{array}\right.$$

(46)

$A_1\left(z\right)$ and $B_1\left(z\right)$ are two analytical functions in the whole plane; $A_2\left(z\right)$ and $B_2\left(z\right)$ are analytical in the region $\left(-\infty ,c\right)\cup \left(a,\infty \right)$, and the functions ${{\rm I}}_k\left(z\right)$, ${\Gamma }_k\left(z\right)$ and ${\mathrm{H}}_k\left(z\right)$ are bounded at infinity. Based on Liouvill’s theorem, we can conclude that $A_1\left(z\right)$ and $B_1\left(z\right)$ are constants in the whole plane. That is

$$A_1\left(z\right)\equiv A$$

(47)

$$B_1\left(z\right)\equiv B$$

(48)

After algebraic manipulations, Equations (43) and (44) can be represented in the form

$$\left\{\begin{array}{l}{r}_1^{\left(\mathrm{I}\right)}{{\rm I}}_{\mathrm{I}}\left(z\right)=g\left(A_2\left(z\right)+r_8^{\left(\mathrm{I}\mathrm{I}\right)}A\right)\hspace{1em}y>0\\ r_1^{\left(\mathrm{I}\mathrm{I}\right)}{{\rm I}}_{\mathrm{I}\mathrm{I}}\left(z\right)=g\gamma \left(A_2\left(z\right)+r_8^{\left(\mathrm{I}\right)}A\right)\hspace{1em}y<0\end{array}\right.$$

(49)

and

$$\left\{\begin{array}{l}{r}_1^{\left(\mathrm{I}\mathrm{I}\right)}{\Gamma }_{\mathrm{I}\mathrm{I}}\left(z\right)=g\left(A_2\left(z\right)+r_8^{\left(\mathrm{I}\mathrm{I}\right)}A\right)-A\hspace{1em}y>0\\ r_1^{\left(\mathrm{I}\right)}{\Gamma }_{\mathrm{I}}\left(z\right)=g\gamma \left(A_2\left(z\right)+r_8^{\left(\mathrm{I}\right)}A\right)-A\hspace{1em}y<0\end{array}\right.$$

(50)

$$\left\{\begin{array}{l}{\mathrm{H}}_{{\rm I}{\rm I}}\left(z\right)=\frac{B}{R^{\left({\rm I}{\rm I}\right)}}\hspace{1em}y>0\\ {\mathrm{H}}_{{\rm I}}\left(z\right)=\frac{B}{R^{\left({\rm I}\right)}}\hspace{1em}y<0\end{array}\right.$$

(51)

where

$g=\frac{r_1^{\left(\mathrm{I}\right)}{r}_1^{\left(\mathrm{I}\mathrm{I}\right)}}{r_1^{\left(\mathrm{I}\right)}{r}_5^{\left(\mathrm{I}\mathrm{I}\right)}+r_1^{\left(\mathrm{I}\mathrm{I}\right)}{r}_6^{\left(\mathrm{I}\right)}}, \gamma =\frac{r_1^{\left(\mathrm{I}\right)}{r}_5^{\left(\mathrm{I}\mathrm{I}\right)}+r_1^{\left(\mathrm{I}\mathrm{I}\right)}{r}_6^{\left(\mathrm{I}\right)}}{r_1^{\left(\mathrm{I}\mathrm{I}\right)}{r}_5^{\left(\mathrm{I}\right)}+r_1^{\left(\mathrm{I}\right)}{r}_6^{\left(\mathrm{I}\mathrm{I}\right)}}, r_8^{\left(k\right)}=\frac{r_5^{\left(k\right)}}{r_1^{\left(k\right)}}$.

The expressions for stresses in the phonon field can be rewritten as

$${\sigma }_{yy}\left(x,y\right)-\mathrm{i}{\sigma }_{xy}\left(x,y\right)=\left\{\begin{array}{l}32\left(g\left(A_2\left(z\right)+r_8^{\left(\mathrm{I}\mathrm{I}\right)}A\right)+g\gamma \left(A_2\left(z\right)+r_8^{\left(\mathrm{I}\right)}A\right)-A+g\left(z-\overline{z}\right)\overline{A_2^{\prime }\left(z\right)}\right)\hspace{1em}y>0\\ 32\left(g\gamma \left(A_2\left(z\right)+r_8^{\left(\mathrm{I}\right)}A\right)+g\left(A_2\left(z\right)+r_8^{\left(\mathrm{I}\mathrm{I}\right)}A\right)-A+g\gamma \left(z-\overline{z}\right)\overline{A_2^{\prime }\left(z\right)}\right)\hspace{1em}y<0\end{array}\right.$$

(52)

and

$${\sigma }_{xx}\left(x,y\right)+{\sigma }_{yy}\left(x,y\right)=128\left\{\begin{array}{l}g\mathrm{Re}\left(A_2\left(z\right)+r_8^{\left(\mathrm{I}\mathrm{I}\right)}A\right),\hspace{1em}y>0\\ g\gamma \mathrm{Re}\left(A_2\left(z\right)+r_8^{\left(\mathrm{I}\right)}A\right),\hspace{1em}y<0\end{array}\right.$$

(53)

The expressions for stresses in the phason field can be rewritten as

$$H_{yy}\left(x,y\right)-\mathrm{i}{H}_{xy}\left(x,y\right)=-32\left\{\begin{array}{l}B+R^{\left({\rm I}\right)}\left(z-\overline{z}\right)\frac{g}{r_1^{\left(\mathrm{I}\right)}}\left(\gamma A_2^{\prime }\left(\overline{z}\right)-\overline{A_2^{\prime }\left(z\right)}\right)+\frac{1}{2}{\left(z-\overline{z}\right)}^2\overline{A_2^{″}\left(z\right)},\hspace{1em}y>0\\ B+R^{\left({\rm I}{\rm I}\right)}\left(z-\overline{z}\right)\frac{g}{r_1^{\left(\mathrm{I}\mathrm{I}\right)}}\left(A_2^{\prime }\left(\overline{z}\right)-\gamma \overline{A_2^{\prime }\left(z\right)}\right)+\frac{1}{2}\gamma {\left(z-\overline{z}\right)}^2\overline{A_2^{″}\left(z\right)},\hspace{1em}y<0\end{array}\right.$$

(54)

and

$$H_{xx}\left(x,y\right)+H_{yy}\left(x,y\right)=-64\left\{\begin{array}{l}\frac{R^{\left({\rm I}\right)}}{r_1^{\left(\mathrm{I}\right)}}\mathrm{Re}\left(g\left(\gamma {\overline{A}}_2\left(z\right)-A_2\left(z\right)+\left(\overline{z}-z\right)A_2^{\prime }\left(z\right)\right)+\left(g\gamma r_8^{\left(\mathrm{I}\right)}-g{r}_8^{\left(\mathrm{I}\mathrm{I}\right)}-1\right)A\right),\hspace{1em}y>0\\ \frac{R^{\left({\rm I}{\rm I}\right)}}{r_1^{\left(\mathrm{I}\mathrm{I}\right)}}\mathrm{Re}\left(g\left({\overline{A}}_2\left(z\right)-\gamma A_2\left(z\right)+\left(\overline{z}-z\right)\gamma A_2^{\prime }\left(z\right)\right)+\left(g{r}_8^{\left(\mathrm{I}\mathrm{I}\right)}-g\gamma r_8^{\left(\mathrm{I}\right)}-1\right)A\right),\hspace{1em}y<0\end{array}\right.$$

(55)

The expressions for displacements in the phonon and phason fields can be rewritten in the form

$$\frac{\partial \left(u_x\left(x,y\right)+\mathrm{i}{u}_y\left(x,y\right)\right)}{\partial x}=32\left\{\begin{array}{l}{r}_6^{\left({\rm I}\right)}g\left(A_2\left(z\right)+r_8^{\left({\rm I}{\rm I}\right)}A\right)\\ -r_5^{\left({\rm I}\right)}\left(g\gamma \left(A_2\left(z\right)+r_8^{\left({\rm I}\right)}A\right)-A+g\left(z-\overline{z}\right)\overline{A_2^{\prime }\left(z\right)}\right),\hspace{1em}y>0\\ r_6^{\left({\rm I}{\rm I}\right)}g\gamma \left(A_2\left(z\right)+r_8^{\left({\rm I}\right)}A\right)\\ -r_5^{\left({\rm I}{\rm I}\right)}\left(g\left(A_2\left(z\right)+r_8^{\left({\rm I}{\rm I}\right)}A\right)-A+g\gamma \left(z-\overline{z}\right)\overline{A_2^{\prime }\left(z\right)}\right),\hspace{1em}y<0\end{array}\right.$$

(56)

and

$$\frac{\partial \left(w_x\left(x,y\right)+i{w}_y\left(x,y\right)\right)}{\partial x}=\left\{\begin{array}{l}32r_7^{\left({\rm I}\right)}\left(\frac{B}{R^{\left({\rm I}\right)}}+\frac{g}{r_1^{\left(\mathrm{I}\right)}}\left(z-\overline{z}\right)\left(\gamma A_2^{\prime }\left(\overline{z}\right)-\overline{{A^{\prime }}_2\left(z\right)}+\frac{1}{2}\left(z-\overline{z}\right)\overline{A_2^{″}\left(z\right)}\right)\right),\hspace{1em}y>0\\ 32r_7^{\left({\rm I}{\rm I}\right)}\left(\frac{B}{R^{\left({\rm I}{\rm I}\right)}}+\frac{g}{r_1^{\left(\mathrm{I}\mathrm{I}\right)}}\left(z-\overline{z}\right)\left(A_2^{\prime }\left(\overline{z}\right)-\gamma \overline{{A^{\prime }}_2\left(z\right)}+\frac{1}{2}\gamma \left(z-\overline{z}\right)\overline{A_2^{″}\left(z\right)}\right)\right),\hspace{1em}y<0\end{array}\right.$$

(57)

Using Equations (49)–(54), we can introduce the following functions:

$$F\left(z\right)=A_2\left(z\right)+pA$$

(58)

$$H\left(z\right)=B_2\left(z\right)$$

(59)

Substituting Equations (58) and (59) into Equations (54) and (55), one obtains

$$\frac{{\sigma }_{yy}\left(x,y\right)-\mathrm{i}{\sigma }_{xy}\left(x,y\right)}{32g}=\left\{\begin{array}{l}F\left(z\right)+\gamma F\left(\overline{z}\right)+\left(z-\overline{z}\right){\overline{F}}^{\prime }\left(\overline{z}\right)\hspace{1em}y>0\\ \gamma F\left(z\right)+F\left(\overline{z}\right)+\gamma \left(z-\overline{z}\right){\overline{F}}^{\prime }\left(\overline{z}\right)\hspace{1em}y<0\end{array}\right.$$

(60)

and

$${\sigma }_{xx}\left(x,y\right)+{\sigma }_{yy}\left(x,y\right)=128\left\{\begin{array}{l}g\mathrm{Re}\left(F\left(z\right)+\left(r_8^{\left(\mathrm{I}\mathrm{I}\right)}-p\right)A\right),\hspace{1em}y>0\\ g\mathrm{Re}\left(\gamma F\left(z\right)+\gamma \left(r_8^{\left(\mathrm{I}\right)}-p\right)A\right),\hspace{1em}y<0\end{array}\right.$$

(61)

$$\frac{H_{yy}\left(x,y\right)-\mathrm{i}{H}_{xy}\left(x,y\right)}{-32g}=\left\{\begin{array}{l}{R}^{\left({\rm I}\right)}\left({\mathrm{H}}_{{\rm I}}\left(\overline{z}\right)+\frac{\left(z-\overline{z}\right)}{r_1^{\left({\rm I}\right)}}\left(\left(\gamma F^{\prime }\left(\overline{z}\right)-{\overline{F}}^{\prime }\left(\overline{z}\right)\right)+\frac{1}{2}\left(z-\overline{z}\right){\overline{F}}^{″}\left(\overline{z}\right)\right)\right),\hspace{1em}y>0\\ R^{\left({\rm I}{\rm I}\right)}\left({\mathrm{H}}_{{\rm I}{\rm I}}\left(\overline{z}\right)+\frac{\left(z-\overline{z}\right)}{r_1^{\left({\rm I}{\rm I}\right)}}\left(F^{\prime }\left(\overline{z}\right)-\gamma {\overline{F}}^{\prime }\left(\overline{z}\right)\right)+\frac{\gamma }{2}\left(z-\overline{z}\right){\overline{F}}^{″}\left(\overline{z}\right)\right),\hspace{1em}y<0\end{array}\right.$$

(62)

and

$$H_{xx}\left(x,y\right)+H_{yy}\left(x,y\right)=-64\left\{\begin{array}{l}\frac{R^{\left({\rm I}\right)}}{r_1^{\left(\mathrm{I}\right)}}\mathrm{Re}\left(g\left(\gamma \overline{F}\left(z\right)-F\left(z\right)+\left(\overline{z}-z\right)F^{\prime }\left(z\right)\right)+gp\left(1-\gamma \right)A+\left(g\gamma r_8^{\left(\mathrm{I}\right)}-g{r}_8^{\left(\mathrm{I}\mathrm{I}\right)}-1\right)A\right),\hspace{1em}y>0\\ \frac{R^{\left({\rm I}{\rm I}\right)}}{r_1^{\left(\mathrm{I}\mathrm{I}\right)}}\mathrm{Re}\left(g\left(\overline{F}\left(z\right)-\gamma F\left(z\right)+\left(\overline{z}-z\right)\gamma F^{\prime }\left(z\right)\right)+gp\left(\gamma -1\right)A+\left(g{r}_8^{\left(\mathrm{I}\mathrm{I}\right)}-g\gamma r_8^{\left(\mathrm{I}\right)}-1\right)A\right),\hspace{1em}y<0\end{array}\right.$$

(63)

$$\frac{\partial \left(u_x+\mathrm{i}{u}_y\right)}{\partial x}=32g\left\{\begin{array}{l}\frac{r_6^{\left({\rm I}\right)}}{r_1^{\left({\rm I}\right)}}\left(F\left(z\right)+r_8^{\left({\rm I}{\rm I}\right)}A-pA\right)\\ -\frac{r_5^{\left({\rm I}\right)}}{r_1^{\left({\rm I}\right)}}\left(\gamma \left(F\left(\overline{z}\right)-pA+r_8^{\left({\rm I}\right)}A\right)-\frac{A}{g}+\left(z-\overline{z}\right){\overline{F}}^{\prime }\left(\overline{z}\right)\right)\hspace{1em}y>0\\ \frac{r_6^{\left({\rm I}{\rm I}\right)}}{r_1^{\left({\rm I}{\rm I}\right)}}\gamma \left(F\left(z\right)+r_8^{\left({\rm I}\right)}A-pA\right)\\ -\frac{r_5^{\left({\rm I}{\rm I}\right)}}{r_1^{\left({\rm I}{\rm I}\right)}}\left(\left(F\left(\overline{z}\right)+r_8^{\left({\rm I}{\rm I}\right)}A-pA\right)-\frac{A}{g}+\left(z-\overline{z}\right)\gamma {\overline{F}}^{\prime }\left(\overline{z}\right)\right)\hspace{1em}y<0\end{array}\right.$$

(64)

$$\frac{\partial }{\partial x}\left(w_x\left(x,y\right)+\mathrm{i}{w}_y\left(x,y\right)\right)=\left\{\begin{array}{l}32r_7^{\left({\rm I}\right)}\left({\mathrm{H}}_{{\rm I}}\left(\overline{z}\right)+\frac{g}{r_1^{\left({\rm I}\right)}}\left(z-\overline{z}\right)\left(\left(\gamma F^{\prime }\left(\overline{z}\right)-{\overline{F}}^{\prime }\left(\overline{z}\right)\right)+\frac{1}{2}\left(z-\overline{z}\right){\overline{F}}^{″}\left(\overline{z}\right)\right)\right),\hspace{1em}y>0\\ 32r_7^{\left({\rm I}{\rm I}\right)}\left({\mathrm{H}}_{{\rm I}{\rm I}}\left(\overline{z}\right)+\frac{g}{r_1^{\left({\rm I}{\rm I}\right)}}\left(z-\overline{z}\right)\left(\left(F^{\prime }\left(\overline{z}\right)-\gamma {\overline{F}}^{\prime }\left(\overline{z}\right)\right)+\frac{\gamma }{2}\left(z-\overline{z}\right){\overline{F}}^{″}\left(\overline{z}\right)\right)\right),\hspace{1em}y<0\end{array}\right.$$

(65)

Thus, the complex function expressions of the stresses and displacement jump derivatives on the interface are written as

$$\left\{\begin{array}{l}{\sigma }_{yy}^{\left({\rm I}\right)}\left(x,0\right)-\mathrm{i}{\sigma }_{xy}^{\left({\rm I}\right)}\left(x,0\right)=32g\left(F^{+}\left(x\right)+\gamma F^{-}\left(x\right)\right)\\ \frac{\partial }{\partial x}\left(\left(u_x^{\left({\rm I}\right)}\left(x\right)+\mathrm{i}{u}_y^{\left({\rm I}\right)}\left(x\right)\right)-\left(u_x^{\left({\rm I}{\rm I}\right)}\left(x\right)+\mathrm{i}{u}_y^{\left({\rm I}{\rm I}\right)}\left(x\right)\right)\right)=32\left(F^{+}\left(x\right)-F^{-}\left(x\right)\right)\end{array}\right.$$

(66)

$$\left\{\begin{array}{l}{H}_{yy}^{\left({\rm I}\right)}\left(x,0\right)-\mathrm{i}{H}_{xy}^{\left({\rm I}\right)}\left(x,0\right)=-32\left(K_1^{\left({\rm I}\right)}-K_2^{\left({\rm I}\right)}\right){\mathrm{H}}^{-}\left(x\right)\\ \frac{\partial }{\partial x}\left(\left(w_x^{\left({\rm I}\right)}\left(x,0\right)+\mathrm{i}{w}_y^{\left({\rm I}\right)}\left(x,0\right)\right)-\left(w_x^{\left({\rm I}{\rm I}\right)}\left(x,0\right)+\mathrm{i}{w}_y^{\left({\rm I}{\rm I}\right)}\left(x,0\right)\right)\right)=32\left(r_7^{\left(\mathrm{I}\right)}{\mathrm{H}}^{-}\left(x\right)-r_7^{\left(\mathrm{II}\right)}{\mathrm{H}}^{+}\left(x\right)\right)\end{array}\right.$$

(67)

5. Complex Potential Solution of the Problem

From the derivation in the previous section, the problem is transformed into the homogeneous Dirichlet–Riemann boundary value problem

$$\left\{\begin{array}{l}{F}^{+}\left(x\right)+\gamma F^{-}\left(x\right)=0,\hspace{1em}x\in \left(c,b\right)\\ \mathrm{Im}{F}^{\pm }\left(x\right)=0,\hspace{1em}x\in \left(b,a\right)\end{array}\right.$$

(68)

$$\left\{\begin{array}{l}{H}^{-}\left(x\right)=0,\hspace{1em}x\in \left(c,b\right)\\ \mathrm{Im}{H}^{\mp }\left(x\right)=0,\hspace{1em}x\in \left(b,a\right)\end{array}\right.$$

(69)

From the second equation of the above equations, one obtains

$$\mathrm{H}\left(z\right)=0$$

(70)

Thus, the stresses and displacements of the phason field can be expressed as

$$\frac{H_{yy}\left(x,y\right)-\mathrm{i}{H}_{xy}\left(x,y\right)}{-32g}=\left\{\begin{array}{l}\frac{R^{\left({\rm I}\right)}}{r_1^{\left({\rm I}\right)}}\left(z-\overline{z}\right)\left(\left(\gamma F^{\prime }\left(\overline{z}\right)-{\overline{F}}^{\prime }\left(\overline{z}\right)\right)+\frac{1}{2}\left(z-\overline{z}\right){\overline{F}}^{″}\left(\overline{z}\right)\right),\hspace{1em}y>0\\ \frac{R^{\left({\rm I}{\rm I}\right)}}{r_1^{\left({\rm I}{\rm I}\right)}}\left(z-\overline{z}\right)\left(\left(F^{\prime }\left(\overline{z}\right)-\gamma {\overline{F}}^{\prime }\left(\overline{z}\right)\right)+\frac{\gamma }{2}\left(z-\overline{z}\right){\overline{F}}^{″}\left(\overline{z}\right)\right),\hspace{1em}y<0\end{array}\right.$$

(71)

$$\frac{\partial }{\partial x}\left(w_x\left(x,y\right)+\mathrm{i}{w}_y\left(x,y\right)\right)=\left\{\begin{array}{l}\frac{32g{r}_7^{\left(1\right)}}{r_1^{\left({\rm I}\right)}}\left(z-\overline{z}\right)\left(\left(\gamma F^{\prime }\left(\overline{z}\right)-{\overline{F}}^{\prime }\left(\overline{z}\right)\right)+\frac{1}{2}\left(z-\overline{z}\right){\overline{F}}^{″}\left(\overline{z}\right)\right),\hspace{1em}y>0\\ \frac{32g{r}_7^{\left(1{\rm I}\right)}}{r_1^{\left({\rm I}{\rm I}\right)}}\left(z-\overline{z}\right)\left(\left(F^{\prime }\left(\overline{z}\right)-\gamma {\overline{F}}^{\prime }\left(\overline{z}\right)\right)+\frac{\gamma }{2}\left(z-\overline{z}\right){\overline{F}}^{″}\left(\overline{z}\right)\right),\hspace{1em}y<0\end{array}\right.$$

(72)

when $x\in \left(-\infty ,c\right)\cup \left(a,\infty \right)$, $F^{+}\left(x\right)=F^{-}\left(x\right)=F\left(x\right)$ is valid.

Using the conditions at infinity, one obtains

$$32g\left(1+\gamma \right)F\left(x\right)=\sigma -\mathrm{i}\tau ,\hspace{1em}x\in \left(-\infty ,c\right)\cup \left(a,\infty \right)$$

(73)

Function $F\left(z\right)$ is analytic at infinity; one obtains

$${F\left(z\right)|}_{z\to \infty }=\frac{\sigma -\mathrm{i}\tau }{32g\left(1+\gamma \right)}$$

(74)

The general solution of Equation (68) of the combined Dirichlet–Riemann boundary value problem from [30] is unbounded at all points $a$, $b$, $c$ and can be written as

$$F\left(z\right)=\mathrm{i}P\left(z\right)X_1\left(z\right)+Q\left(z\right)X_2\left(z\right)$$

(75)

where

$P\left(z\right)=C_{1}z+C_2, Q\left(z\right)=D_{1}z+D_2,$

$X_1\left(z\right)=\frac{e^{\mathrm{i}\phi \left(z\right)}}{\sqrt{\left(z-a\right)}\sqrt{\left(z-c\right)}}, X_2\left(z\right)=\frac{e^{\mathrm{i}\phi \left(z\right)}}{\sqrt{\left(z-b\right)}\sqrt{\left(z-c\right)}}$

$\phi \left(z\right)=2\epsilon \ln \frac{\sqrt{\left(a-b\right)\left(z-c\right)}}{\sqrt{\left(a-c\right)\left(z-b\right)}+\sqrt{\left(b-c\right)\left(z-a\right)}}$

$\epsilon =\frac{1}{2\pi }\ln \gamma$

Replacing $z$ with $\overline{z}$ in Equations (31)–(33), one has

$$F\left(\overline{z}\right)=\mathrm{i}P\left(\overline{z}\right)X_1\left(\overline{z}\right)+Q\left(\overline{z}\right)X_2\left(\overline{z}\right)$$

(76)

derived by differentiation

$$\begin{array}{l}{F}^{\prime }\left(z\right)=\mathrm{i}{C}_1{X}_1\left(z\right)+\mathrm{i}P\left(z\right)X_1^{\prime }\left(z\right)+D_1{X}_2\left(z\right)+Q\left(z\right)X_2^{\prime }\left(z\right)\\ F^{″}\left(z\right)=2\mathrm{i}{C}_1{X}_1^{\prime }\left(z\right)+\mathrm{i}P\left(z\right)X_1^{″}\left(z\right)+2D_1{X}_2^{\prime }\left(z\right)+Q\left(z\right)X_2^{″}\left(z\right)\\ \overline{F^{\prime }\left(z\right)}=\mathrm{i}{C}_1\overline{X_1\left(z\right)}+\mathrm{i}\overline{P\left(z\right)}\overline{X_1^{\prime }\left(z\right)}+D_1\overline{X_2\left(z\right)}+\overline{Q\left(z\right)}\overline{X_2^{\prime }\left(z\right)}\\ \overline{F^{″}\left(z\right)}=2\mathrm{i}{C}_1\overline{X_1^{\prime }\left(z\right)}+\mathrm{i}\overline{P\left(z\right)}\overline{X_1^{″}\left(z\right)}+2D_1\overline{X_2^{\prime }\left(z\right)}+\overline{Q\left(z\right)}\overline{X_2^{″}\left(z\right)}\end{array}$$

(77)

$X_1\left(z\right)$ and $X_2\left(z\right)$ have an oscillating singularity at the point $z=c+\mathrm{i}0$ and square-root singularities at the points $z=a+\mathrm{i}0$ and $z=b+\mathrm{i}0$, respectively. $X_1\left(z\right)$ and $X_2\left(z\right)$ at infinity can be written as

$$\begin{array}{l}{X}_1\left(z\right)=z^{-2}{e}^{\mathrm{i}\beta }\left(z+\mathrm{i}{\beta }_1+\frac{c+a}{2}\right)+O\left(z^{-3}\right)\\ X_2\left(z\right)=z^{-2}{e}^{\mathrm{i}\beta }\left(z+\mathrm{i}{\beta }_1+\frac{c+b}{2}\right)+O\left(z^{-3}\right)\end{array}$$

(78)

where

$\beta =\epsilon \ln \frac{\sqrt{a-c}-\sqrt{b-c}}{\sqrt{a-c}+\sqrt{b-c}},{\beta }_1=\epsilon \sqrt{\left(a-c\right)\left(b-c\right)},$

The arbitrary constants $C_1$, $C_2$, $D_1$, $D_2$ have the following forms

$$\begin{array}{l}{C}_1=-\tau \cos \beta -\sigma \sin \beta ,C_2=-\frac{c+a}{2}{C}_1-{\beta }_1{D}_1,\\ D_1=\sigma \cos \beta -\tau \sin \beta ,D_2={\beta }_1{C}_1-\frac{c+b}{2}{D}_1,\end{array}$$

(79)

The stresses and the derivatives of the displacement jumps for $z=x+\mathrm{i}0$ can be expressed as follows:

for $x>a$:

$$\begin{array}{l}{\sigma }_{yy}^{\left({\rm I}\right)}\left(x,0\right)-\mathrm{i}{\sigma }_{xy}^{\left({\rm I}\right)}\left(x,0\right)=32g\left(1+\gamma \right)\left(\frac{Q\left(x\right)}{\sqrt{x-b}}+\mathrm{i}\frac{P\left(x\right)}{\sqrt{x-a}}\right)\frac{\exp \left[\mathrm{i}\phi \left(x\right)\right]}{\sqrt{x-c}}\\ H_{yy}^{\left({\rm I}\right)}\left(x,0\right)-\mathrm{i}{H}_{xy}^{\left({\rm I}\right)}\left(x,0\right)=0\end{array}$$

(80)

for $x\in \left(b,a\right)$:

$$\begin{array}{l}{\sigma }_{yy}^{\left(1\right)}=\frac{32gP\left(x\right)}{\sqrt{\left(x-c\right)\left(a-x\right)}}\left(\left(1-\gamma \right)ch{\phi }_0\left(x\right)+\left(1+\gamma \right)sh{\phi }_0\left(x\right)\right)\\ \hspace{1em}\hspace{1em}+\frac{32gQ\left(x\right)}{\sqrt{\left(x-c\right)\left(x-b\right)}}\left(\left(1+\gamma \right)ch{\phi }_0\left(x\right)+\left(1-\gamma \right)sh{\phi }_0\left(x\right)\right)\end{array}$$

(81)

$${\left(u_x^{\left({\rm I}\right)}-u_x^{\left({\rm I}{\rm I}\right)}\right)}^{\prime }=\frac{2}{\sqrt{x-c}}\left(\frac{P\left(x\right)}{\sqrt{a-x}}ch{\phi }_0\left(x\right)+\frac{Q\left(x\right)}{\sqrt{x-b}}sh{\phi }_0\left(x\right)\right)$$

(82)

for $x\in \left(c,b\right)$:

$$\begin{array}{l}{\left(u_x^{\left({\rm I}\right)}\left(x\right)+\mathrm{i}{u}_y^{\left({\rm I}\right)}\left(x\right)\right)}^{\prime }-{\left(u_x^{\left({\rm I}{\rm I}\right)}\left(x\right)+\mathrm{i}{u}_y^{\left({\rm I}{\rm I}\right)}\left(x\right)\right)}^{\prime }=\frac{32\left(1+\gamma \right)}{\sqrt{\gamma }}\left(\frac{P\left(x\right)}{\sqrt{a-x}}-\mathrm{i}\frac{Q\left(x\right)}{\sqrt{b-x}}\right)\frac{\exp \left(\mathrm{i}{\phi }^{*}\left(x\right)\right)}{\sqrt{x-c}}\\ {\left(w_x^{\left({\rm I}\right)}\left(x\right)+\mathrm{i}{w}_y^{\left({\rm I}\right)}\left(x\right)\right)}^{\prime }-{\left(w_x^{\left({\rm I}{\rm I}\right)}\left(x\right)+\mathrm{i}{w}_y^{\left({\rm I}{\rm I}\right)}\left(x\right)\right)}^{\prime }=0\end{array}$$

(83)

where

$\begin{array}{l}{\phi }^{*}\left(x\right)=2\epsilon \ln \frac{\sqrt{\left(a-b\right)\left(x-c\right)}}{\sqrt{\left(a-c\right)\left(b-x\right)}+\sqrt{\left(b-c\right)\left(a-x\right)}},\\ {\phi }_0\left(x\right)=2\epsilon \arctan \sqrt{\frac{\left(b-c\right)\left(a-x\right)}{\left(a-c\right)\left(x-b\right)}}\end{array}$

From Equations (80) and (82), it can be seen that the normal stress in the phonon field has a square-root singularity for $x\to b+0$, and the shear stress has the same singularity for $x\to a+0$. The relevant stress intensity factors in the phonon field can be defined as

$$K_{S1}=\underset{x\to b+0}{\lim }\sqrt{2\left(x-b\right)}{\sigma }_{yy}\left(x,0\right), K_{S2}=\underset{x\to a+0}{\lim }\sqrt{2\left(x-a\right)}{\sigma }_{xy}\left(x,0\right),$$

(84)

and can be written as

$$K_{S1}=64g\sqrt{\pi \gamma }\frac{Q\left(b\right)}{\sqrt{b-c}}, K_{S2}=-32g\left(1+\gamma \right)\sqrt{\frac{2\pi }{a-c}}P\left(a\right)$$

(85)

Further, the stress intensity factors can be expressed in the form

$$\begin{array}{l}{K}_{S1}=\frac{\sqrt{\pi \gamma }}{1+\gamma }\left(\sqrt{b-c}\left(\sigma \cos \beta -\tau \sin \beta \right)-2\epsilon \sqrt{a-c}\left(\sigma \sin \beta +\tau \cos \beta \right)\right)\\ K_{S2}=\sqrt{\frac{\pi }{2}}\left(\sqrt{a-c}\left(\sigma \sin \beta +\tau \cos \beta \right)+2\epsilon \sqrt{b-c}\left(\sigma \cos \beta -\tau \sin \beta \right)\right)\end{array}$$

(86)

The asymptotic behavior of the stresses and the displacement jumps in the phonon field at the points $a$ and $b$ can be written as

$$\begin{array}{l}{{\sigma }_{yy}^{\left({\rm I}\right)}\left(x,0\right)|}_{x\to b+0}=\frac{K_{S1}}{\sqrt{2\pi \left(x-b\right)}}, {{\sigma }_{xy}^{\left({\rm I}\right)}\left(x,0\right)|}_{x\to a+0}=\frac{K_{S2}}{\sqrt{2\pi \left(x-a\right)}}, \\ {\left(u_y^{\left({\rm I}\right)}\left(x,0\right)-u_y^{\left({\rm I}{\rm I}\right)}\left(x,0\right)\right)|}_{x\to b-0}=\frac{K_{S1}}{16g\sqrt{2\pi \gamma }}\sqrt{b-x}, \\ {\left(u_x^{\left({\rm I}\right)}\left(x,0\right)-u_x^{\left({\rm I}{\rm I}\right)}\left(x,0\right)\right)|}_{x\to a-0}=\frac{K_{S2}}{8g\left(1+\gamma \right)\sqrt{2\pi }}\sqrt{a-x}\end{array}$$

(87)

6. Numerical Results and Discussion

The elastic constants of Al Ni Co quasicrystal alloy are taken as the elastic constants of the phonon field, phason field, and coupling constants of the phonon–phason field [31], as shown in

Table 1. Material I constants of decagonal quasicrystal.

Elastic Constants	The Value of Elastic Constant
Phonon field elastic constant/GPa	$C_{11}=234.33$$, C_{66}=88.46$
Phason field elastic constant/GPa	$K_1=122$$, K_2=24$
Coupling constant/GPa	$R=0.8846$

and

Table 2. Material II constants of decagonal quasicrystal.

Elastic Constants	The Value of Elastic Constant
Phonon field elastic constant/GPa	$C_{11}=214.33$$, C_{66}=68.46$
Phason field elastic constant/GPa	$K_1=102$$, K_2=22$
Coupling constant/GPa	$R=0.8646$

. In order to avoid matrix ill condition caused by material parameters in different orders of magnitude, the material constants are dimensionless, and the dimensionless stress intensity factors in the phonon field are obtained.

It is clear from Equations (84)–(86) that the length of the contact zone depends only on point $b$. Figure 2 shows the change in the stress intensity factor with the crack contact zone, where the relative contact zone length of a crack with a right contact zone is given. It can be found that for any point, the model framework that only considers the contact area can clearly define the area with a large contact area. The smaller the contact length, the greater the normal stress intensity factor. This is consistent with the trend of classical elasticity.

7. Conclusions

The interface crack contact zone of decagonal quasicrystal bi-materials under far-field mixed loading is studied. Based on the theory of complex variable function, the problem is transformed into a Dirichlet–Riemann problem for analytical solution. The expressions of stress, stress intensity factor, and displacement jump along the material interface are obtained by using the closed analytical formula of the interface crack in the single contact zone of the decagonal quasicrystal bi-materials, and the relationship between the fracture mechanics parameters of the interface crack is given. The analytical expression obtained can be used to verify some numerical analysis and can also accurately show the physical nature of the crack problem in the contact zone of decagonal quasicrystal materials.