Scattering and stress concentration of SH guided waves by a semicircular hole on the boundary of an infinite piezoelectric ceramic plate

In this paper, the dynamic issue related to a semicircular depression in a piezoelectric ceramic plate is investigated. The fundamental mechanical problems of dynamic stress concentration factor and electric field intensity concentration factor are addressed using the complex function approach, wave function expansion method, and repeated mirror method, and the analytical expressions for the relevant stress concentration are presented. The influences of the dimensionless parameter wave number, strip thickness, and guided wave order on the stress concentration of the SH guided wave incident within the strip are explored.The results indicate that, with a fixed plate thickness, the 0 order high frequency guided wave exerts a more pronounced influence on the stress concentration near the hole, resulting in more severe damage. Under the electromechanical coupling effect, the detrimental impact of the electric field far exceeds that of the stress field. Consequently, it is essential to focus on the adverse effects imposed by the electric field.

1 Introduction

Piezoelectric ceramics possess extraordinary electromechanical coupling properties and hold an irreplaceable value both in scientific research and practical applications. In terms of research significance, they provide a crucial path for exploring the inherent relationship between the micro - structure and macroscopic properties of materials. This helps to deepen the understanding of the physical properties of crystals and promote the development of basic sciences. In the energy field, the energy - harvesting potential of piezoelectric ceramics is expected to relieve the energy pressure and make self - powered devices possible. At the application level, piezoelectric ceramics can accurately measure parameters such as pressure and acceleration, serving key scenarios like aerospace and automotive safety. In ultrasonic technology, they enable efficient cleaning and precision machining. In the medical field, ultrasonic diagnostic and therapeutic equipment, by leveraging piezoelectric ceramics, achieve accurate disease detection and treatment, greatly enhancing the medical level and efficiency. However, during the design and use of piezoelectric ceramic structures, a series of problems such as void defects often occur. The existence of defects at the boundaries can lead to brittle fracture and failure of the structure. This paper conducts research on the problem of open holes at the boundaries of strip - shaped piezoelectric ceramic structures, taking advantage of the SH guided wave, which features fast propagation speed, high detection efficiency, low attenuation, and strong applicability.

In recent years, the application of equipment made of piezoelectric ceramics has been widely concerned in structural damage detection.Miao et al. [1] introduced a double-layer metasurface aimed at isolating the fundamental shear horizontal wave (SH0 wave). Dutta et al. [2] presented a novel study comparing the effects of single and dual porosity on the propagation of horizontally polarized shear waves (SH waves) in corrugated elastic void materials. Hu et al. [3] studied the propagation of SH waves in elastic semi-infinite media with semi-elliptical craters on the surface. Hemalatha et al. [4] analyzed and studied SH waves in a rotating functionally graded magnetoelectroelastic system composed of an elastic substrate and a linearly changing functionally graded magnetoelectroelastic layer. Wang et al. [5] improved the damage imaging accuracy of anisotropic plates with uniform thickness and variable thickness by combining a new compensation term based on the theoretical excitation directivity of quasi-SH0 modes in anisotropic materials. The MSBPT proposed by Miao et al. [6] is composed of a meta-substrate two-column thickness shear mode PZT wafer, which is achieved by developing a unidirectional SH0 deflection of a substrate-based piezoelectric transducer under a single driving source. Qu et al. [7] studied the dynamic stress concentration coefficient of a semi-cylindrical depression on the surface of an infinite elastic plate under the action of SH-waves and gave an analytical solution. Li et al. [8] used the complex function method to study wave propagation in a non-uniform half-space with semi-cylindrical surface bumps. Guo et al. [9] investigated theoretical methods, function expansion method, complex mirror image method and other methods to study the wave scattering problem in an elastic half-space anisotropic medium containing a semicircular canyon and a movable cylindrical cavity. Qi et al. [10–13] conducted a systematic and comprehensive study on the theoretical models of half space, infinite space, and semi-infinite space that contain holes and inclusions, thereby providing a valuable theoretical foundation for understanding the wave field. Sun et al. [14] examined the displacement amplitude at the interface of an elastic half-space with a semicircular depression covered by a viscous liquid under the action of shear horizontal waves. Zhang et al. [15] considered the exponential distribution of material parameters along the coordinate axis and studied the dynamic characteristics of a circular inclusion in an inhomogeneous piezomagnetic half-space with anti-plane shear wave propagation. Using the theory of complex variables in plane elasticity, Liu et al. [16] proposed a wave function expansion method to analyze the external scattering of SH wave surfaces caused by two symmetrical circular cavities situated in two bonded exponentially graded half-spaces. Li et al. [17] discussed in detail the effects of magnetic field and compressive stress on the mechanical displacement, dynamic stress and magnetic potential of the SH wave around the circular hole. Kuo et al. [18] proposed an exact analysis for the scattering of anti-plane shear wave by a piezoelectric circular cylinder in piezomagnetic matrix with imperfect interfaces. Yang et al. [19] used complex functions and Green’s method to study the dynamic anti-plane characteristics of collinear periodic cracks at Type III interfaces in a bi-material half-space. Singh et al. [20] examined the scattering of Love waves in interface cracks in thin lossy media (viscoelastic materials).The fracture and scattering problems of piezoelectric quasicrystals [21–27] have also received extensive attention. Research in this area has significantly enriched the applications of piezoelectric materials.

The research on boundary perforations in strip - shaped piezoelectric ceramic structures holds immense significance and remarkable value. In the realm of fundamental research, the presence of boundary perforations modifies the stress distribution and vibration characteristics within piezoelectric ceramics. This not only facilitates a profound comprehension of the microscopic mechanisms underlying the piezoelectric effect but also furnishes a theoretical foundation for the optimization of material properties. In the field of acoustics, the insights derived from the study of this model can serve as crucial theoretical references. They play a pivotal role in optimizing the performance of ultrasonic transducers, enhancing the precision of medical imaging techniques, and improving the efficacy of non - destructive testing methods.

2 Steady state governing equation

Assuming that the Z-axis is the direction of electrical polarization, the steady-state control Equation of the infinitely long ceramic piezoelectric medium is Equation 1.

$\left\{\begin{array}{l}{c}_{44}{\nabla }^{2}w+{\mathrm{e}}_{15}{\nabla }^2\varphi +\rho {\omega }^{2}w=0\\ e_{15}{\nabla }^{2}w-{\kappa }_{11}{\nabla }^2\varphi =0\end{array}\right.(1)$

$w$, $\varphi$ and $\omega$ are respectively expressed as the out-of-plane displacement, potential and circular frequency of the piezoelectric material. The specific expression of potential is Equation 2.

$\varphi =\frac{e_{15}}{{\kappa }_{11}}\left(w+\phi \right)(2)$

Substitute Equation 2 into Equation 1 and further simplify it to Equation 3.

$\left\{\begin{array}{l}{\nabla }^{2}w+k^{2}w=0\\ {\nabla }^2\phi =0\end{array}\right.(3)$

There is a relationship in the above Equation: $k^2=\frac{\rho {\omega }^2}{c_{44}+\frac{e_{15}^2}{{\kappa }_{11}}}$.

In the polar coordinate system $\left(r,\theta \right)$, the governing Equation can transform the relationship with r and $\theta$ as independent variables into Equation 4.

$\left\{\begin{array}{l}\frac{{\mathrm{\partial }}^{2}w}{\mathrm{\partial }{r}^2}+\frac{1}{r}\frac{\mathrm{\partial }w}{\mathrm{\partial }r}+\frac{1}{r^2}\frac{{\mathrm{\partial }}^{2}w}{\mathrm{\partial }{\theta }^2}+k^{2}w=0\\ \frac{{\mathrm{\partial }}^2\phi }{\mathrm{\partial }{r}^2}+\frac{1}{r}\frac{\mathrm{\partial }\phi }{\mathrm{\partial }r}+\frac{1}{r^2}\frac{{\mathrm{\partial }}^2\phi }{\mathrm{\partial }{\theta }^2}=0\end{array}\right.(4)$

In the polar coordinate system $\left(r,\theta \right)$, the constitutive Equation can be expressed as Equation 5.

$\left\{\begin{array}{l}{\tau }_{rz}=c_{44}\frac{\mathrm{\partial }w}{\mathrm{\partial }r}+e_{15}\frac{\mathrm{\partial }\varphi }{\mathrm{\partial }r}\\ {\tau }_{\theta z}=c_{44}\frac{1}{r}\frac{\mathrm{\partial }w}{\mathrm{\partial }\theta }+e_{15}\frac{1}{r}\frac{\mathrm{\partial }\varphi }{\mathrm{\partial }\theta }\\ D_r=e_{15}\frac{\mathrm{\partial }w}{\mathrm{\partial }r}-{\kappa }_{11}\frac{\mathrm{\partial }\varphi }{\mathrm{\partial }r}\\ D_{\theta }=e_{15}\frac{1}{r}\frac{\mathrm{\partial }w}{\mathrm{\partial }\theta }-{\kappa }_{11}\frac{1}{r}\frac{\mathrm{\partial }\varphi }{\mathrm{\partial }\theta }\end{array}\right.(5)$

In the complex plane $\left(\eta ,\overline{\eta }\right)$, complex variables $\eta =x+\mathrm{i}y$ and $\overline{\eta }=x-\mathrm{i}y$ are introduced. Let $\eta =r{e}^{\mathrm{i}\theta }$ and $\overline{\eta }=r{e}^{-\mathrm{i}\theta }$, the constitutive Equation is further organized into Equation 6.

$\left\{\begin{array}{l}{\tau }_{rz}=c_{44}\left(\frac{\mathrm{\partial }w}{\mathrm{\partial }\eta }{e}^{\mathrm{i}\theta }+\frac{\mathrm{\partial }w}{\mathrm{\partial }\overline{\eta }}{e}^{-\mathrm{i}\theta }\right)+e_{15}\left(\frac{\mathrm{\partial }\varphi }{\mathrm{\partial }\eta }{e}^{\mathrm{i}\theta }+\frac{\mathrm{\partial }\varphi }{\mathrm{\partial }\overline{\eta }}{e}^{-\mathrm{i}\theta }\right)\\ {\tau }_{\theta z}=\mathrm{i}{c}_{44}\left(\frac{\mathrm{\partial }w}{\mathrm{\partial }\eta }{e}^{\mathrm{i}\theta }-\frac{\mathrm{\partial }w}{\mathrm{\partial }\overline{\eta }}{e}^{-\mathrm{i}\theta }\right)+\mathrm{i}{e}_{15}\left(\frac{\mathrm{\partial }\varphi }{\mathrm{\partial }\eta }{e}^{\mathrm{i}\theta }-\frac{\mathrm{\partial }\varphi }{\mathrm{\partial }\overline{\eta }}{e}^{-\mathrm{i}\theta }\right)\\ D_r=e_{15}\left(\frac{\mathrm{\partial }w}{\mathrm{\partial }\eta }{e}^{\mathrm{i}\theta }+\frac{\mathrm{\partial }w}{\mathrm{\partial }\overline{\eta }}{e}^{-\mathrm{i}\theta }\right)-{\kappa }_{11}\left(\frac{\mathrm{\partial }\varphi }{\mathrm{\partial }\eta }{e}^{\mathrm{i}\theta }+\frac{\mathrm{\partial }\varphi }{\mathrm{\partial }\overline{\eta }}{e}^{-\mathrm{i}\theta }\right)\\ D_{\theta }=\mathrm{i}{e}_{15}\left(\frac{\mathrm{\partial }w}{\mathrm{\partial }\eta }{e}^{\mathrm{i}\theta }-\frac{\mathrm{\partial }w}{\mathrm{\partial }\overline{\eta }}{e}^{-\mathrm{i}\theta }\right)-\mathrm{i}{\kappa }_{11}\left(\frac{\mathrm{\partial }\varphi }{\mathrm{\partial }\eta }{e}^{\mathrm{i}\theta }-\frac{\mathrm{\partial }\varphi }{\mathrm{\partial }\overline{\eta }}{e}^{-\mathrm{i}\theta }\right)\end{array}\right.(6)$

3 Structure model

The initial defects of the piezoelectric ceramic plate strip structure do not always occur inside the strip structure. Defects at the edge positions also happen from time to time.To effectively study the damage of piezoelectric ceramic plates, this paper analyzes the problem of a semi - circular hole being drilled at the edges of the piezoelectric strip.The structural model of an infinitely long piezoelectric ceramic plate containing a semi-cylindrical depression is shown in Figure 1. The radius of semicircular recess is $a$. The mass density of piezoelectric ceramic is $\rho$. The elastic constant is $c_{44}$. The piezoelectric constant is $e_{15}$. The dielectric constant is ${\kappa }_{11}$.The upper and lower boundaries of piezoelectric ceramic plate are $B_U$ and $B_L$ respectively. An electric field can be formed inside the depression and the dielectric constant of air is ${\kappa }_0$.Establish a plane coordinate system $\left(o,x,y\right)$ with the center of semicircular depression as coordinate origin.

Figure 1

Figure 1. An infinite length piezoelectric ceramic plate structure model with a semi-cylindrical depression.

4 Wave field in a slab structure

4.1 Incident wave field in a piezoelectric ceramic plate structure

The SH guided wave expression is shown in Equation 7.

$\left\{\begin{array}{l}{w}_m=f_m\left(y\right)\exp \left[\mathrm{i}\left(k_{m}x-\omega t\right)\right]\\ {\varphi }_m=\frac{e_{15}}{{\kappa }_{11}}{f}_m\left(y\right)\exp \left[\mathrm{i}\left(k_{m}x-\omega t\right)\right]\end{array}\right.(7)$

The upper and lower boundaries of the piezoelectric ceramic lath structure satisfy stress freedom and electrical insulation conditions as shown in Equation 8.

$\left\{\begin{array}{l}{\tau }_{yz}={\left.c_{44}\frac{\mathrm{\partial }w}{\mathrm{\partial }y}+e_{15}\frac{\mathrm{\partial }\varphi }{\mathrm{\partial }y}\right|}_{y=0,-h}=0\\ {\left.D_y=e_{15}\frac{\mathrm{\partial }w}{\mathrm{\partial }y}-{\kappa }_{11}\frac{\mathrm{\partial }\varphi }{\mathrm{\partial }y}\right|}_{y=0,-h}=0\end{array}\right.(8)$

$f_m\left(y\right)$ is the interference term in the y direction satisfying Equation 9. m is the order of guided waves. The SH-guided wave shape is shown in Figure 2. $w_m^1$ and $w_m^2$ are the amplitudes of the corresponding propagating guided waves. When m is an even number, $w_m^1=0$. When m is an odd number, $w_m^2=0$. Only when $k_m$ is a real number, $\exp \left[\mathrm{i}\left(k_{m}x-\omega t\right)\right]$ can represent a traveling wave propagating in the x-axis direction. When the m-order SH guided wave is incident, the required wave number should meet $k>m\pi /h$.

$f_m\left(y\right)=w_m^1\sin \left[q_m\left(y+\frac{h}{2}\right)\right]+w_m^2\cos \left[q_m\left(y+\frac{h}{2}\right)\right](9)$

Figure 2

Figure 2. Mode of SH guided wave.

The condition that $q_m$ satisfies is Equation 10.

$q_m=\frac{m\pi }{h}(10)$

$k_m$ and $q_m$ satisfy Equation 11.

$q_m^2=k^2-k_m^2(11)$

By superimposing the guided waves, all displacement waves that comply with the stress-free and electrically insulated conditions of the upper and lower boundaries of the strip piezoelectric material medium can be derived, as illustrated in Equation 12.

$\left\{\begin{array}{l}{w}^{\left(i\right)}=w_0{\displaystyle \sum _{m=0}^{+\infty }}{w}_m=w_0{\displaystyle \sum _{m=0}^{+\infty }}{f}_m\left(y\right)\exp \left[\mathrm{i}\left(k_{m}x-\omega t\right)\right]\\ {\varphi }^{\left(i\right)}=\frac{e_{15}}{{\kappa }_{11}}{w}^{\left(i\right)}=\frac{e_{15}}{{\kappa }_{11}}{w}_0{\displaystyle \sum _{m=0}^{+\infty }}{f}_m\left(y\right)\exp \left[\mathrm{i}\left(k_{m}x-\omega t\right)\right]\end{array}\right.(12)$

The displacement is expressed as Equation 13.

$\left\{\begin{array}{l}{w}^{\left(i\right)}=w_0\left\{w_m^{\left(1\right)}\sin \left[q_m\left(y+\frac{h}{2}\right)\right]+w_m^{\left(2\right)}\cos \left[q_m\left(y+\frac{h}{2}\right)\right]\right\}{e}^{\mathrm{i}{k}_{m}x}\\ {\varphi }^{\left(\mathrm{i}\right)}=\frac{e_{15}}{{\kappa }_{11}}{w}_0\left\{w_m^{\left(1\right)}\sin \left[q_m\left(y+\frac{h}{2}\right)\right]+w_m^{\left(2\right)}\cos \left[q_m\left(y+\frac{h}{2}\right)\right]\right\}{e}^{\mathrm{i}{k}_{m}x}\end{array}\right.(13)$

The stress is expressed as Equation 14.

$\left\{\begin{array}{l}{\tau }_{xz}^{\left(i\right)}=w_0{c}_m\mathrm{i}{k}_{\mathrm{m}}{e}^{\mathrm{i}{k}_{m}x}\left\{w_m^{\left(1\right)}\sin \left[q_m\left(y+\frac{h}{2}\right)\right]+w_m^{\left(2\right)}\cos \left[q_m\left(y+\frac{h}{2}\right)\right]\right\}\\ {\tau }_{yz}^{\left(i\right)}=w_0{c}_m{q}_m{e}^{\mathrm{i}{k}_{m}x}\left\{w_m^{\left(1\right)}\cos \left[q_m\left(y+\frac{h}{2}\right)\right]-w_m^{\left(2\right)}\sin \left[q_m\left(y+\frac{h}{2}\right)\right]\right\}\end{array}\right.(14)$

The polar coordinate expression of stress is Equation 15.

$\left\{\begin{array}{l}{\tau }_{rz}^{\left(i\right)}={\tau }_{xz}^{\left(i\right)}\cos \theta +{\tau }_{yz}^{\left(i\right)}\sin \theta \\ {\tau }_{\theta z}^{\left(i\right)}=-{\tau }_{xz}^{\left(i\right)}\sin \theta +{\tau }_{yz}^{\left(i\right)}\cos \theta \end{array}\right.(15)$

4.2 Scattered wave field in a piezoelectric ceramic plate structure

The expression for scattered waves in the elastic wave field is Equation 16.

$w_0^s=w_0{\displaystyle \sum _{n=-\infty }^{+\infty }}{A}_n{H}_n^{\left(1\right)}\left(k\left|\eta \right|\right){\left(\frac{\eta }{\left|\eta \right|}\right)}^n(16)$

Using the mirror image method, the scattered guided wave is $w_0^s$ and the excited potential function is ${\varphi }_0^s$. The stress freedom and electrical insulation conditions are satisfied on the upper and lower boundaries. Use $p$ to represent the number of mirrors, let $L_1^p=h_1+d_1^P$, $L_2^p=-\left(h_2+d_2^P\right)$. When the number of mirrors $p$ is an odd number, $w_p^{s1}$ (mirroring along the upper boundary) and $w_p^{s2}$ (mirroring along the lower boundary) can be expressed as Equation 17.

$\left\{\begin{array}{l}{w}_p^{s1}=w_0{\displaystyle \sum _{n=-\infty }^{+\infty }}{A}_n{H}_n^{\left(1\right)}\left(k\left|\eta -\mathrm{i}{L}_1^p\right|\right){\left(\frac{\eta -\mathrm{i}{L}_1^p}{\left|\eta -\mathrm{i}{L}_1^p\right|}\right)}^{-n}\\ w_p^{s2}=w_0{\displaystyle \sum _{n=-\infty }^{+\infty }}{A}_n{H}_n^{\left(1\right)}\left(k\left|\eta -\mathrm{i}{L}_2^p\right|\right){\left(\frac{\eta -\mathrm{i}{L}_2^p}{\left|\eta -\mathrm{i}{L}_2^p\right|}\right)}^{-n}\end{array}\right.(17)$

When $p$ is an even number, $w_p^{s1}$ (mirroring along the upper boundary) and $w_p^{s2}$ (mirroring along the lower boundary) can be expressed as Equation 18.

$\left\{\begin{array}{l}{w}_p^{s1}=w_0{\displaystyle \sum _{n=-\infty }^{+\infty }}{A}_n{H}_n^{\left(1\right)}\left(k\left|\eta -\mathrm{i}{L}_1^p\right|\right){\left(\frac{\eta -\mathrm{i}{L}_1^p}{\left|\eta -\mathrm{i}{L}_1^p\right|}\right)}^n\\ w_p^{s2}=w_0{\displaystyle \sum _{n=-\infty }^{+\infty }}{A}_n{H}_n^{\left(1\right)}\left(k\left|\eta -\mathrm{i}{L}_2^p\right|\right){\left(\frac{\eta -\mathrm{i}{L}_2^p}{\left|\eta -\mathrm{i}{L}_2^p\right|}\right)}^n\end{array}\right.(18)$

When $p$ is an odd number, $d_1^p$ and $d_2^p$ can be expressed as Equation 19.

$\left\{\begin{array}{l}{d}_1^p=h_1+\left(p-1\right)\left(h_1+h_2\right)\\ d_2^p=h_2+\left(p-1\right)\left(h_1+h_2\right)\end{array}\right.(19)$

When $p$ is an even number, $d_1^p$ and $d_2^p$ can be expressed as Equation 20.

$\left\{\begin{array}{l}{d}_1^p=h_2+\left(p-1\right)\left(h_1+h_2\right)\\ d_2^p=h_1+\left(p-1\right)\left(h_1+h_2\right)\end{array}\right.(20)$

The scattered wave expression of the elastic wave field is Equation 21.

$w^s=w_0^s+{\displaystyle \sum _{p=1}^{+\infty }}\left(w_p^{s1}+w_p^{s2}\right)(21)$

The potential expression for the scattered wave and its excitation is Equation 22.

${\varphi }^s=\frac{e_{15}}{{\kappa }_{11}}\left(w_0^s+{\phi }_0^s\right)+\frac{e_{15}}{{\kappa }_{11}}{\displaystyle \sum _{p=1}^{+\infty }}\left(w_p^{s1}+{\phi }_p^{s1}+w_p^{s2}+{\phi }_p^{s2}\right)(22)$

The specific expression of ${\phi }_0^s$ is Equation 23.

${\phi }_0^s=w_0{\displaystyle \sum _{n=1}^{\infty }}\left[B_n{\left|\eta \right|}^{-n}{\left(\frac{\eta }{\left|\eta \right|}\right)}^{-n}+C_n{\left|\overline{\eta }\right|}^{-n}{\left(\frac{\overline{\eta }}{\left|\overline{\eta }\right|}\right)}^{-n}\right](23)$

When $p$ is an odd number, ${\phi }_p^{s1}$ and ${\phi }_p^{s2}$ can be expressed as Equation 24.

$\left\{\begin{array}{l}{\phi }_p^{s1}=w_0{\displaystyle \sum _{n=1}^{\infty }}\left[B_n{\left|\overline{\eta }+\mathrm{i}{L}_1^p\right|}^{-n}{\left(\frac{\overline{\eta }+\mathrm{i}{L}_1^p}{\left|\overline{\eta }+\mathrm{i}{L}_1^p\right|}\right)}^{-n}+C_n{\left|\eta -\mathrm{i}{L}_1^p\right|}^{-n}{\left(\frac{\eta -\mathrm{i}{L}_1^p}{\left|\eta -\mathrm{i}{L}_1^p\right|}\right)}^{-n}\right]\\ {\phi }_p^{s2}=w_0{\displaystyle \sum _{n=1}^{\infty }}\left[B_n{\left|\overline{\eta }+\mathrm{i}{L}_2^p\right|}^{-n}{\left(\frac{\overline{\eta }+\mathrm{i}{L}_2^p}{\left|\overline{\eta }+\mathrm{i}{L}_2^p\right|}\right)}^{-n}+C_n{\left|\eta -\mathrm{i}{L}_2^p\right|}^{-n}{\left(\frac{\eta -\mathrm{i}{L}_2^p}{\left|\eta -\mathrm{i}{L}_2^p\right|}\right)}^{-n}\right]\end{array}\right.(24)$

When $p$ is an even number, ${\phi }_p^{s1}$ and ${\phi }_p^{s2}$ can be expressed as Equation 25.

$\left\{\begin{array}{l}{\phi }_p^{s1}=w_0{\displaystyle \sum _{n=1}^{\infty }}\left[B_n{\left|\eta +\mathrm{i}{L}_1^p\right|}^{-n}{\left(\frac{\eta +\mathrm{i}{L}_1^p}{\left|\eta +\mathrm{i}{L}_1^p\right|}\right)}^{-n}+C_n{\left|\overline{\eta }-\mathrm{i}{L}_1^p\right|}^{-n}{\left(\frac{\overline{\eta }-\mathrm{i}{L}_1^p}{\left|\overline{\eta }-\mathrm{i}{L}_1^p\right|}\right)}^{-n}\right]\\ {\phi }_p^{s2}=w_0{\displaystyle \sum _{n=1}^{\infty }}\left[B_n{\left|\eta +\mathrm{i}{L}_2^p\right|}^{-n}{\left(\frac{\eta +\mathrm{i}{L}_2^p}{\left|\eta +\mathrm{i}{L}_2^p\right|}\right)}^{-n}+C_n{\left|\overline{\eta }-\mathrm{i}{L}_2^p\right|}^{-n}{\left(\frac{\overline{\eta }-\mathrm{i}{L}_2^p}{\left|\overline{\eta }-\mathrm{i}{L}_2^p\right|}\right)}^{-n}\right]\end{array}\right.(25)$

The displacement and potential expression of medium I in the infinitely long piezoelectric ceramic plate structure is Equation 26.

$w=w^i+w^s,\varphi ={\varphi }^i+{\varphi }^s(26)$

4.3 Wavefield in the depression

There is only electric field but no elastic field in the depression. The potential inside the hole satisfies the Laplace Equation ${\nabla }^2{\varphi }^{\mathrm{c}}=0$, that is, Equation 27.

${\varphi }^c=\frac{e_{15}}{{\kappa }_{11}}{w}_0{D}_0+\frac{e_{15}}{{\kappa }_{11}}{w}_0{\displaystyle \sum _{n=1}^{+\infty }}\left[D_n{k}^n{\left|\eta \right|}^n{\left(\frac{\eta }{\left|\eta \right|}\right)}^n+E_n{k}^n{\left|\overline{\eta }\right|}^n{\left(\frac{\overline{\eta }}{\left|\overline{\eta }\right|}\right)}^n\right](27)$

5 Boundary conditions and definite solutions to the system of integral equations

The boundary conditions for stress freedom, potential and normal electric displacement continuity are Equation 28.

${\tau }_{rz}={\tau }_{rz}^i+{\tau }_{rz}^s=0,D_{\mathrm{r}}=D_{\mathrm{r}}^{\mathrm{c}},\varphi ={\varphi }^i+{\varphi }^s={\varphi }^{\mathrm{c}}(28)$

The definite solution integral Equation set established from Equation 28 is as follows Equation 29.

$\left\{\begin{array}{l}{\displaystyle \sum _{n=-\infty }^{+\infty }}{A}_n{\xi }_n^{\left(11\right)}+{\displaystyle \sum _{n=1}^{+\infty }}{B}_n{\xi }_n^{\left(12\right)}+{\displaystyle \sum _{n=1}^{+\infty }}{C}_n{\xi }_n^{\left(13\right)}={\zeta }^{\left(1\right)}\\ {\displaystyle \sum _{n=1}^{+\infty }}{B}_n{\xi }_n^{\left(22\right)}+{\displaystyle \sum _{n=1}^{+\infty }}{C}_n{\xi }_n^{\left(23\right)}+{\displaystyle \sum _{n=1}^{+\infty }}{D}_n{\xi }_n^{\left(24\right)}+{\displaystyle \sum _{n=1}^{+\infty }}{E}_n{\xi }_n^{\left(25\right)}={\zeta }^{\left(2\right)}\\ {\displaystyle \sum _{n=-\infty }^{+\infty }}{A}_n{\xi }_n^{\left(31\right)}+{\displaystyle \sum _{n=1}^{+\infty }}{B}_n{\xi }_n^{\left(32\right)}+{\displaystyle \sum _{n=1}^{+\infty }}{C}_n{\xi }_n^{\left(33\right)}+{\displaystyle \sum _{n=0}^{+\infty }}{D}_n{\xi }_n^{\left(34\right)}+{\displaystyle \sum _{n=1}^{+\infty }}{E}_n{\xi }_n^{\left(35\right)}={\zeta }^{\left(3\right)}\end{array}\right.(29)$

Multiply the left and right sides of Equation system (29) by $e^{-\mathrm{i}f\theta }\left(f=0,\pm 1,\pm 2,...,\pm N\right)$. This Equation is integrated at $\left(-\pi ,\pi \right)$ to obtain Equation 30.

$\left\{\begin{array}{l}{\displaystyle \sum _{n=-\infty }^{+\infty }}{A}_n{\xi }_{fn}^{\left(11\right)}+{\displaystyle \sum _{n=1}^{+\infty }}{B}_n{\xi }_{fn}^{\left(12\right)}+{\displaystyle \sum _{n=1}^{+\infty }}{C}_n{\xi }_{fn}^{\left(13\right)}={\zeta }_f^{\left(1\right)}\\ {\displaystyle \sum _{n=1}^{+\infty }}{B}_n{\xi }_{fn}^{\left(22\right)}+{\displaystyle \sum _{n=1}^{+\infty }}{C}_n{\xi }_{fn}^{\left(23\right)}+{\displaystyle \sum _{n=1}^{+\infty }}{D}_n{\xi }_{fn}^{\left(24\right)}+{\displaystyle \sum _{n=1}^{+\infty }}{E}_n{\xi }_{fn}^{\left(25\right)}={\zeta }_f^{\left(2\right)}\\ {\displaystyle \sum _{n=-\infty }^{+\infty }}{A}_n{\xi }_{fn}^{\left(31\right)}+{\displaystyle \sum _{n=1}^{+\infty }}{B}_n{\xi }_{fn}^{\left(32\right)}+{\displaystyle \sum _{n=1}^{+\infty }}{C}_n{\xi }_{fn}^{\left(33\right)}+{\displaystyle \sum _{n=0}^{+\infty }}{D}_n{\xi }_{fn}^{\left(34\right)}+{\displaystyle \sum _{n=1}^{+\infty }}{E}_n{\xi }_{fn}^{\left(35\right)}={\zeta }_f^{\left(3\right)}\end{array}\right.(30)$

Where $\left\{\begin{array}{l}{\xi }_{fn}^{\left(jl\right)}=\frac{1}{2\pi }{\int }_{-\pi }^{\pi }{\xi }_n^{\left(jl\right)}{e}^{-\mathrm{i}f\theta }d\theta \\ {\zeta }_f^{\left(j\right)}=\frac{1}{2\pi }{\int }_{-\pi }^{\pi }{\zeta }^{\left(j\right)}{e}^{-\mathrm{i}f\theta }d\theta \end{array}\right.$, $\left(j=1,2,3;l=1,2,3,4,5\right)$.

Equation 30 is an infinite system of Equations for solving the coefficients $A_n$, $B_n$, $C_n$, $D_n$ and $E_n$. The cylindrical function has convergence. If the largest positive integer of the sum is N, the number of unknown coefficients is 2N + 1, N, N, N+1, N, N, N+1, N, a total of 6N + 2 unknowns. The number of Equations is (2N + 1) +2N +(2N + 1) = 6N + 2 Equations. The problem is solvable with the same number of unknowns as the number of Equations.

6 Stress concentration

The dynamic stress concentration factor can be expressed as Equation 31.

${\tau }_{\theta \mathrm{z}}^{*}=\left|\frac{{\tau }_{\theta \mathrm{z}}^{\mathrm{{\rm I}}}}{{\tau }_0}\right|(31)$

Where ${\tau }_{\theta \mathrm{z}}^{\mathrm{{\rm I}}}={\tau }_{\theta \mathrm{z}}^i+{\tau }_{\theta z}^s$; ${\tau }_0=w_0{c}_{m}k$ is the amplitude of incident wave shear stress.

The electric field intensity coefficient can be expressed as Equation 32.

$E_{\theta }^{*}=\left|\frac{E_{\theta }^{\mathrm{{\rm I}}}}{E_0}\right|(32)$

Among them $E_{\theta }^{\mathrm{{\rm I}}}=E_{\theta }^i+E_{\theta }^s$, its specific expression is $E_{\theta }=-\mathrm{i}\left(\frac{\mathrm{\partial }\varphi {\rm I}}{\mathrm{\partial }\eta }{\mathrm{e}}^{\mathrm{i}\theta }-\frac{\mathrm{\partial }\varphi {\rm I}}{\mathrm{\partial }\overline{\eta }}{\mathrm{e}}^{-\mathrm{i}\theta }\right)$. $E_0=\frac{e_{15}}{{\kappa }_{11}}{w}_{0}k$ is the amplitude of the electric field intensity of the incident wave.

7 Results and analysis

This paper uses dimensionless parameter analysis. Since the dielectric constant of most piezoelectric media is three orders of magnitude greater than that of a cavity (vacuum or air), ${\kappa }_{11}/{\kappa }_0=1000$ is selected. h* = h/a represents the ratio of the thickness of the ceramic piezoelectric plate to the radius of the recess. m is the order of guided waves. ka is the wave numb.

Figure 3 discusses the stress concentration at the hole edge as a function of h* for ka = 0.5 and m = 1. It can be seen from the figure that the electric field intensity concentration coefficient at the hole edge is obviously greater than the dynamic stress concentration coefficient. When h* = 3, the maximum value of the electric field intensity coefficient at the edge of the recess is about 2.5 times larger than the maximum value of the dynamic stress coefficient. With the increase of dimensionless h*, the stress concentration coefficient at the edge of the hole tends to weaken. The value corresponding to h* = 3 is 5.18 times higher than that corresponding to h* = 10.The value corresponding to h* = 3 is 5.02 times higher than that corresponding to h* = 10. It is not difficult to see that the distance of the lower boundary has a significant reduction in the stress concentration at the hole edge. When the thickness of the ceramic piezoelectric plate reaches a certain value, the stress concentration distribution gradually becomes a regular circle, which is basically close to about 270°. When SH guided waves impinge upon piezoelectric ceramic structures harboring defects, due attention must be accorded to the harm inflicted by electric fields. In the event of structural damage, appropriately regulating the plate thickness can effectively mitigate stress concentration, thereby safeguarding the piezoelectric ceramic plate from being damaged.

Figure 3

Figure 3. Distribution of stress around the depression with h* when ka = 0.5: (A) Dynamic stress concentration coefficient (B) Electric field intensity concentration coefficient.

Figure 4 shows the change of stress concentration at the hole edge with h* when m = 1 and ka = 1.0. It is different from the physical phenomenon in Figure 3. When h* = 3, the stress concentration at the edge of the depression decreases significantly and the stress distribution diagram is no longer a single petal. As h* increases, the stress concentration does not always decrease but reaches a maximum at h* = 5. As h* increases, the stress concentration factor gradually approaches from the facing wave surface to the back wave surface. When it increases to a certain value, it is basically close to the position directly below the concave edge. It can be concluded that when designing this kind of ceramic piezoelectric plate. If there is a depression on the boundary, it is not always the most dangerous right below the hole edge. The electric field strength at the hole edge is still greater than the dynamic stress. Similarly, in engineering practice, we should still focus on the adverse effects of the electric field.

Figure 4

Figure 4. Distribution of stress around the depression with h* when ka = 1.0: (A) Dynamic stress concentration coefficient (B) Electric field intensity concentration coefficient.

Figure 5 shows the change of stress concentration at the hole edge with h* when m = 1 and ka = 1.5. Figure 5 can see that as h* increases. Both ${\tau }_{\theta \mathrm{z}}*$ and $E_{\theta }^{*}$ on the concave edge tend to increase. When h* = 10, the maximum values of hole edges ${\tau }_{\theta \mathrm{z}}*$ and $E_{\theta }^{*}$ are respectively 1.64 times and 1.65 times larger than the maximum values of ${\tau }_{\theta \mathrm{z}}*$ and $E_{\theta }^{*}$ corresponding to h* = 3. With the increase of thickness, the stress map shifts gradually from the front wave surface to the back wave surface. At this wave number, it is essential to maintain precise control over the thickness of the ceramic piezoelectric plate during the design process. In the presence of defects at the boundary, the thickness of the strip should be adjusted accordingly to mitigate stress concentration and prevent damage.

Figure 5

Figure 5. Distribution of stress around the depression with h* when ka = 1.5: (A) Dynamic stress concentration coefficient (B) Electric field intensity concentration coefficient.

Figure 6 shows the change of stress concentration at the hole edge with h* when m = 1 and ka = 2.0. It can be seen from the figure that the electric field intensity concentration coefficient at the edge of the recess increases little compared with the dynamic stress concentration coefficient. The impact and changes of the two are basically the same. As h* increases, the values of ${\tau }_{\theta \mathrm{z}}*$ and $E_{\theta }^{*}$ increase significantly. When the wave number reaches a certain value, the change of stress concentration is not very obvious. The effect of electromechanical coupling is essentially similar to that of the influence. Appropriately adjusting the incident wave number serves as an effective approach to mitigate the mechanical and electrical stress concentrations at the hole’s edge. This, in turn, diminishes the structural damage resulting from stress concentration.

Figure 6

Figure 6. Distribution of stress around the depression with h* when ka = 2.0: (A) Dynamic stress concentration coefficient (B) Electric field intensity concentration coefficient.

Figure 7 shows the change of stress concentration at the hole edge with the guided wave order m when h* = 5 and ka = 0.5. Figure 7A shows that the dynamic stress concentration coefficient at the recessed edge increases as the number of guided waves increases. The maximum value ${\tau }_{\theta \mathrm{z}}*$ corresponding to the guided wave order m = 2 is 6.07 times and 5.38 times larger than the values corresponding to m = 0 and m = 1, respectively. As the number of guided waves increases, the peak of ${\tau }_{\theta \mathrm{z}}*$ moves forward in turn. When the guided wave order m = 2, the wave peak appears at a position of 214°. Similarly, Figure 7B has a similar change pattern. When the guided wave orders are m = 0, 1, and 2, the electric field intensity concentration coefficient at the hole edge is nearly 2–3 times higher than the dynamic stress concentration coefficient. This manifests that the influence exerted by the electric field assumes a preeminent role. In the event of damage to ceramic plates and belts, extraordinary attention ought to be dedicated to the deleterious impacts imposed by the electric field.

Figure 7

Figure 7. Distribution of stress around the depression with guided wave order m when ka = 0.5: (A) Dynamic stress concentration coefficient (B) Electric field intensity concentration coefficient.

Figure 8 shows the change of stress concentration at the hole edge with the guided wave order m when h* = 5 and ka = 1.0. The ${\tau }_{\theta \mathrm{z}}*$ and $E_{\theta }^{*}$ of recessed edge do not always increase with the increase of guided wave order. When m = 2, the corresponding maximum value for high-order guided waves is 3–4 times higher than that for m = 0 and m = 1. When the guided wave order m = 2, there are two peaks in the trend of first increasing and then oscillating attenuation. The stress concentration diagrams of m = 0 and m = 1 first increase to reach the peak and then decrease. At this wave number, the $E_{\theta }*$ at the hole edge is greater than ${\tau }_{\theta \mathrm{z}}*$ in three different orders. When SH guided waves propagate at a frequency with ka = 1.0, it is essential to concentrate on the hazards brought about by the incidence of high - order guided waves.

Figure 8

Figure 8. Distribution of stress around the depression with guided wave order m when ka = 1.0. (A) Dynamic stress concentration coefficient (B) Electric field intensity concentration coefficient.

Figure 9 shows the change of stress concentration at the hole edge with the guided wave order m when h* = 5 and ka = 1.5. The dynamic stress and electric field strength at the edge of the hole have basically the same trend with the number of guided waves. When the guided wave order m = 0, the peaks of ${\tau }_{\theta \mathrm{z}}*$ and $E_{\theta }^{*}$ are larger than the values corresponding to the other two orders. When m = 2, the waveform changes more significantly. When this form of SH guided wave is incident, the adverse effects of relatively low-order electric fields should be paid attention to.

Figure 9

Figure 9. Distribution of stress around the depression with waveguide order m when ka = 1.5: (A) Dynamic stress concentration coefficient (B) Electric field intensity concentration coefficient.

Figure 10 shows how the stress concentration at the hole edge changes with the guided wave order m when h* = 5 and ka = 2.0. Under the action of mechano-electric coupling, the changing trends of the dynamic stress coefficient and the electric field intensity coefficient at the hole edge are similar. The effect of electric field is not significantly greater than that of force. The stress concentration caused by the incidence of 0th-order guided waves is significantly higher than that of first-order and second-order guided waves and shows a weakening trend. The order in which stress concentration reaches its peak is second order, first order and 0th order, all showing the trend of oscillation attenuation. Under the influence of this frequency band, the impact of the 0th-order guided wave incidence is more serious.

Figure 10

Figure 10. Distribution of stress around the depression with waveguide order m when ka = 2.0: (A) Dynamic stress concentration coefficient (B) Electric field intensity concentration coefficient.

8 Conclusion

In this paper, based on the wave theory of SH guided waves, the complex variable function method, repeated image method and guided wave expansion method are used to study the stress concentration problem at the hole edge. The effects of dimensionless parameters ka, h* and order m on the dynamic stress coefficient and electric field intensity concentration coefficient at the hole edge are discussed. Relevant conclusions are drawn to protect the safe use of piezoelectric ceramic plate structures.

(1) When the first-order guided wave is incident, the wave number ka = 0.5, and the plate thickness is relatively thick, the stress concentration at the edge of the hole can be effectively reduced, paying special attention to the influence under the action of the electric field. When the wave number ka > 0.5 and the plate thickness is relatively thin, the reduction of stress at the hole edge will cause less damage to the ceramic piezoelectric plate. Therefore, under different wave numbers, the plate thickness needs to be reasonably controlled to protect the piezoelectric ceramic element from damage.

(2) When the plate thickness is constant, the wave numbers are ka = 0.5 and ka = 1.0 respectively, and the impact of the incident high-order guided waves is significantly higher than that of the low-order guided waves. When ka = 1.5 and ka = 2.0, the impact of the zero-order guided wave incidence on the stress concentration at the hole edge is more obvious, resulting in more serious damage. For different wave numbers SH guided waves incident, reasonable selection of the guided wave order can reduce the stress concentration to a certain extent.

(3) The laws governing the influence of plate thickness, wave number, and guided wave order are not static. The dimensionless parameters under discussion must be deliberated in light of specific conditions. Generally, under the force electricity coupling, the damage to ceramic piezoelectric plates induced by electric fields is comparatively more significant. Therefore, greater emphasis should be placed on the impact of electric fields.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author contributions

EQ: Funding acquisition, Software, Writing–original draft. HQ: Data curation, Writing–review and editing. JG: Validation, Writing–review and editing. SL: Methodology, Software, Writing–review and editing. YL: Resources, Investigation, Writing–original draft.

Funding

The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported by the National Natural Cultivation General Project of Basic Scientific Research Operating Expenses for Universities in Heilongjiang Province (145409328).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

1. Miao H, Cao X, Fu M. Double-layer metasurface for blocking the fundamental SH wave. Smart Mater Structures (2024) 33(9):095044. doi:10.1088/1361-665X/ad7215

CrossRef Full Text | Google Scholar

2. Dutta R, Das S, Ahmad H, Gupta V. Comparative analysis of double and single porosity effects on SH-wave induced vibrations in periodic porous lattices. Soil Dyn Earthquake Eng (2024) 186:108919. doi:10.1016/j.soildyn.2024.108919

CrossRef Full Text | Google Scholar

3. Hu H, Dai M, Gao CF. Scattering of SH wave in an elastic half-space by a semi-elliptical crater with surface elasticity. Appl Math Model (2024) 135:759–71. doi:10.1016/j.apm.2024.07.017

CrossRef Full Text | Google Scholar

4. Hemalatha K, Kumar S. Propagation of SH-wave in a rotating functionally graded magneto-electro-elastic structure with corrugated interface. Mech Sol (2024) 59(3):1635–58. doi:10.1007/s42417-024-01365-5

CrossRef Full Text | Google Scholar

5. Wang J, Li B. SH-wave based defect imaging method for CFRP plates with variable thickness. Nondestructive Test Eval (2024) 1-27:1–27. doi:10.1080/10589759.2024.2358384

CrossRef Full Text | Google Scholar

6. Miao H, Du Y. Metasubstrate-based SH guided wave piezoelectric transducer for unidirectional beam deflection without time delay. Smart Mater Structures (2024) 33(1):015038. doi:10.1088/1361-665X/ad1890

CrossRef Full Text | Google Scholar

7. Qu E, Qi H, Guo J, Wang L, Yang J, Liu S. Dynamic response analysis of SH-guided waves in a strip-shaped elastic medium for a semi-cylindrical depression. Archive Appl Mech (2023) 93(3):1241–58. doi:10.1007/s00419-022-02325-9

CrossRef Full Text | Google Scholar

8. Li X, Song Y, Yang Z, Carrera E. Scattering of SH waves by a semi-cylindrical bump in an inhomogeneous half-space. Int J Appl Mech (2023) 15(04):2250034. doi:10.1142/S175882512250034X

CrossRef Full Text | Google Scholar

9. Guo D, Bian J, Song Y, Yang Y, Yang Z. Response of semicircular canyons and movable cylindrical cavities to SH waves in anisotropic half-space geology. Geophysics (2024) 89(5):C197–C209. doi:10.1190/geo2023-0598.1

CrossRef Full Text | Google Scholar

10. Qi H, Chu F, Guo J, Yang R. Surface motion of a half-space containing an elliptical-arc canyon under incident SH waves. Mathematics (2020) 8(11):1884. doi:10.3390/math8111884

CrossRef Full Text | Google Scholar

11. Hui Q, Fuqing C, Jing G. Dynamic stress concentration of an elliptical cavity in a semi-elliptical hill under SH-waves. Earthquake Eng Eng Vibration (2021) 20:347–59. doi:10.1007/s11803-021-2024-9

CrossRef Full Text | Google Scholar

12. Qi H, Xiang M, Guo J. The dynamic stress analysis of an infinite piezoelectric material strip with a circular cavity. Mech Adv Mater Structures (2021) 28(17):1818–26. doi:10.1080/15376494.2019.1709676

CrossRef Full Text | Google Scholar

13. Fan Z, Qi H, Zhang Y, Wu H, Chu F, Guo J. Anti-plane dynamic response characteristics of a semi-infinite plate with cylindrical hole defect. Thin-Walled Structures (2024) 202:112038. doi:10.1016/j.tws.2024.112038

CrossRef Full Text | Google Scholar

14. Sun B, Guo J, Zhang M, Zhang G. Dynamic response of a half-space with a depression covered by viscous liquid under SH waves. Soil Dyn Earthquake Eng (2022) 154:107139. doi:10.1016/j.soildyn.2021.107139

CrossRef Full Text | Google Scholar

15. Zhang X, Qi H. Propagation of SH-waves in inhomogeneous piezoelectric/piezomagnetic half-space with circular inclusion. Acta Mechanica (2022) 233(9):3829–52. doi:10.1007/s00707-022-03315-2

CrossRef Full Text | Google Scholar

16. Liu Q, Zhao M, Liu Z. Wave function expansion method for the scattering of SH waves by two symmetrical circular cavities in two bonded exponentially graded half spaces. Eng Anal Boundary Elem (2019) 106:389–96. doi:10.1016/j.enganabound.2019.05.015

CrossRef Full Text | Google Scholar

17. Li Q, Gu C, Lei D, Ou Z. Scattering of shear horizontal (SH) waves by a circular hole in an infinite piezomagnetic material. Acta Mechanica Solida Sinica (2024) 1-12:891–902. doi:10.1007/s10338-024-00508-1

CrossRef Full Text | Google Scholar

18. Kuo HY, Yu SH. Effect of the imperfect interface on the scattering of SH wave in a piezoelectric cylinder in a piezomagnetic matrix. Int J Eng Sci (2014) 85:186–202. doi:10.1016/j.ijengsci.2014.08.006

CrossRef Full Text | Google Scholar

19. Yang J, Liu S, Liu Y, Yang B. Dynamic stress intensity factor for periodic interfacial cracks in a bi-material half space disturbed by SH waves. J Electromagn Waves Appl (2022) 36(18):2769–84. doi:10.1080/09205071.2022.2106901

CrossRef Full Text | Google Scholar

20. Singh S, Saw GR, Chakraborty G. Scattering of Love wave from an interface crack in a viscoelastic waveguide layer bonded to a piezoelectric substrate: an analytical estimate and numerical validation. Mech Adv Mater Structures (2024) 31(16):3875–88. doi:10.1080/15376494.2023.2186547

CrossRef Full Text | Google Scholar

21. Yang J, Xu Y, Ding S, Li X. Dynamic fracture of a partially permeable crack in a functionally graded one-dimensional hexagonal piezoelectric quasicrystal under a time-harmonic elastic SH-wave. Mathematics Mech Sol (2023) 28(9):1939–58. doi:10.1177/10812865221138838

CrossRef Full Text | Google Scholar

22. Yang Z, Yu X, Xu W, Xu C, Zhou Z. Accurate fracture analysis of electrically permeable/impermeable cracks in one-dimensional hexagonal piezoelectric quasicrystal junction. Mathematics Mech Sol (2019) 24(12):4032–50. doi:10.1177/1081286519865002

CrossRef Full Text | Google Scholar

23. Ma Y, Zhao X, Lu S, Yang J, Zhou Y, Ding S. Dynamic behavior of interface cracks in 1D hexagonal piezoelectric quasicrystal coating–substrate structures subjected to plane waves. J Eng Mech (2024) 150(12):04024090. doi:10.1061/JENMDT.EMENG-7808

CrossRef Full Text | Google Scholar

24. Liu Y, Tang X, Duan P, Wang T. Three-dimensional thermo-electro-elastic field in one-dimensional hexagonal piezoelectric quasi-crystal weakened by an elliptical crack. Mathematics Mech Sol (2022) 27(7):1233–54. doi:10.1177/10812865211059219

CrossRef Full Text | Google Scholar

25. Ma Y, Zhou Y, Yang J, Lu S, Zhao X, Ding S. Scattering of plane waves from an interface crack between the 1D hexagonal quasicrystals coating and the elastic substrate. Acta Mechanica (2024) 236:289–303. doi:10.1007/s00707-024-04153-0

CrossRef Full Text | Google Scholar

26. Hu K, Yang W, Fu J, Chen Z, Gao CF. Analysis of an anti-plane crack in a one-dimensional orthorhombic quasicrystal strip. Mathematics Mech Sol (2022) 27(11):2467–79. doi:10.1177/10812865211073814

CrossRef Full Text | Google Scholar

27. Hu K, Yang W, Fu J, Chen Z, Gao CF. Analysis of a mode-I crack in a one-dimensional orthorhombic quasicrystal strip. Mathematics Mech Sol (2023) 28(3):635–52. doi:10.1177/10812865221091748

CrossRef Full Text | Google Scholar