^{1}

^{*}

^{1}

^{1}

In this paper, on the basis of Liu’s complex function and conformal mapping methods, supplemented by local coordinate system method, e-type piezoelectric material and elastic wave scattering and dynamic stress concentrations problems with double holes question are studied, and an analytical solution is given to the problems. On the basis of multiple scattering of elastic wave theory, put forward the study about microscopic dynamics model to dynamic stress in the structure of piezoelectric composites as well as dynamic playing field. As an example, the numerical results of the dynamic stress distribution around the hole in case double equal diameter holes are given in the paper, and the influence of incident wave number and hole-spacing parameters on the dynamic stress concentration factor is analyzed.

In engineering practice, various holes of different shapes are often needed in piezoelectric materials and structures. The hole spacing will have an influence on the stress state. If the hole spacing is approximately solved according to a single hole, a large error is bound to occur.

In order to solve this problem, Hu [

If adopting elastic dynamics theory and wave function expansion method, approximately solving the row of holes problem according to a single hole will produce larger errors. On the basis of single arbitrary hole, further solve the elastic wave scattering and dynamic stress concentrations problems in piezoelectric media with double holes.

The governing equation and constitutive relation of the steady-state inverse plane dynamics problem of piezoelectrics used in this chapter are shown in (1) and (2).

4 ∂ 2 w ∂ ζ ∂ ζ ¯ + k 2 w = 0 4 ∂ 2 w ∂ ζ ∂ ζ ¯ = 0 (1)

τ x z = c 44 ( ∂ w ∂ ζ + ∂ w ∂ ζ ¯ ) + e 15 ( ∂ ϕ ∂ ζ + ∂ ϕ ∂ ζ ¯ )

τ y z = i c 44 ( ∂ w ∂ ζ − ∂ w ∂ ζ ¯ ) + i e 15 ( ∂ ϕ ∂ ζ − ∂ ϕ ∂ ζ ¯ )

D x = e 15 ( ∂ w ∂ ζ + ∂ w ∂ ζ ¯ ) − κ 11 ( ∂ ϕ ∂ ζ + ∂ ϕ ∂ ζ ¯ ) (2)

D y = i e 15 ( ∂ w ∂ ζ − ∂ w ∂ ζ ¯ ) − i κ 11 ( ∂ ϕ ∂ ζ − ∂ ϕ ∂ ζ ¯ )

τ x z + i τ y z = 2 c 44 ∂ w ∂ ζ ¯ + 2 e 15 ∂ ϕ ∂ ζ ¯

D x + i D y = 2 e 15 ∂ w ∂ ζ ¯ − 2 κ 11 ∂ ϕ ∂ ζ ¯

Considering the case of an infinite piezoelectric material with two holes, a steady-state SH electroacoustic wave is incident along × axis, ignoring the time factor, the corresponding out-of-plane displacement field w ( i ) and in-plane potential field ϕ ( i ) can be indicated as:

w ( i ) = w 0 ∑ n = − ∞ ∞ i n J n ( k | Ω ( η ) | ) { Ω ( η ) | Ω ( η ) | } n ϕ ( i ) = e 15 κ 11 w ( i ) } (3)

In the analysis of calculation, the local coordinate method can be used to convert the internal force component in the local polar coordinate system to the polar coordinate system to be calculated. Considering the multiple scattering between each hole, the elastic wave scattering field generated by the MTH hole in polar coordinate system ( r m , θ m ) can be described as follows:

w ( s ) = ∑ m = 1 2 ∑ n = − ∞ ∞ A n m H n ( 1 ) ( k r m ) e i n θ m ϕ ( s ) = e 15 k 11 w ( s ) + ∑ m = 1 2 ∑ n = 0 ∞ B n m ( k r m ) − n e i n θ m } (4)

Of which, A n m , B n m ( m = 1 , 2 ) are respectively scattering wave mode coefficient produced by the m-th hole, determined by boundary conditions.

The total wave field of the anti-plane shear wave should be superposed by the incident field and the scattering field when solving the open-hole boundary value problem of the dielectric, then the total field of the piezoelectric material with double holes is:

w ( t ) = w ( i ) + w ( s ) ϕ ( t ) = ϕ ( i ) + ϕ ( s ) } (5)

w ( t ) = ∑ n = − ∞ ∞ [ w 0 i n J n ( k r ) + ∑ m = 1 2 A n m H n ( 1 ) ( k r m ) ] e i n θ m (6)

ϕ ( t ) = e 15 κ 11 ∑ n = − ∞ ∞ ∑ m = 1 2 [ w 0 i n J n ( k r ) + ∑ m = 1 2 A n m H n ( 1 ) ( k r m ) ] e i n θ m + ∑ n = 0 ∞ ∑ m = 1 2 B n ( k r m ) − n e i n θ m (7)

There is no elastic displacement field in a circular hole, but only the electric potential field φ c , and the charge density is zero. Therefore, the solution of this equation should satisfy the Laplace equation ∇ 2 φ 2 = 0 . Considering that the electric potential in a circular hole cannot be infinite, it should be finite, so the equation can be written as:

ϕ c = ∑ n = 0 ∞ ∑ m = 1 2 C n ( k r m ) n e i n θ m (8)

Then the corresponding stress can be expressed as:

T ρ z = k c 44 2 ( 1 + λ ) ∑ n = − ∞ ∞ ∑ m = 1 2 { w 0 i n [ η m ρ m Ω ′ ( η m ) | Ω ′ ( η m ) | J n − 1 ( k | Ω ( η m ) | ) { Ω ( η m ) | Ω ( η m ) | } n − 1 − η m ¯ ρ m Ω ′ ( η m ) ¯ | Ω ′ ( η m ) | J n + 1 ( k | Ω ( η m ) | ) { Ω ( η m ) | Ω ( η m ) | } n + 1 ] + A n [ η m ρ m Ω ′ ( η m ) | Ω ′ ( η m ) | H n − 1 ( 1 ) ( k | Ω ( η m ) | ) { Ω ( η m ) | Ω ( η m ) | } n − 1 − η m ¯ ρ m Ω ′ ( η m ) ¯ | Ω ′ ( η m ) | H n + 1 ( 1 ) ( k | Ω ( η m ) | ) { Ω ( η m ) | Ω ( η m ) | } n + 1 ] } − e 15 ∑ n = 0 ∞ ∑ m = 1 2 B n n k − n η m ¯ ρ m Ω ′ ( η m ) ¯ | Ω ′ ( η m ) | ( Ω ( η m ) ¯ ) − n − 1 (9)

D ρ = e 15 ∂ w ∂ r − κ 11 ∂ ϕ ∂ r = κ 11 ∑ n = 0 ∞ ∑ m = 1 2 B n n k − n η m ¯ ρ m Ω ′ ( η m ) ¯ | Ω ′ ( η m ) | ( Ω ( η m ) ¯ ) − n − 1 (10)

D ρ c = − κ 0 ϕ c = − κ 0 ∑ n = 0 ∞ ∑ m = 1 2 C n n k n η m ρ m Ω ′ ( η m ) ¯ | Ω ′ ( η m ) | ( Ω ( η m ) ¯ ) − n − 1

Of which, λ = e 15 2 c 44 κ 11 is dimensionless piezoelectric constant, κ 0 is dielectric constant in vacuum.

Piezoelectric medium with a single arbitrary hole is studied. On η plane, the open hole, stress free, potential and normal electric displacement are set as the free boundary conditions, 6 boundary conditions can be given.

τ ρ m z | ρ 1 = a 1 = 0 D ρ m | ρ 1 = a 1 = D ρ m c | ρ 1 = a 1 ϕ | ρ m = a 1 = ϕ c τ ρ m z | ρ 1 = a 2 = 0 D ρ m | ρ 2 = a 2 = D ρ m c | ρ 2 = a 2 ϕ | ρ m = a 2 = ϕ c } (11)

Of which, a 1 , a 2 are the radii of the double holes.

Substitute Formulas (6), (7) into the opening boundary condition Formula (11), According to the orthogonality of the function system, the six mode coefficients A n 1 , B n 1 , C n 1 , A n 2 , B n 2 , C n 2 to be solved can be determined by:

∑ j = 1 6 ∑ n = − ∞ ∞ E n i X n = E i ( i = 1 , 2 , 3 , 4 , 5 , 6 ) (12)

E n = [ E 11 n E 12 n E 13 n E 14 n E 15 n E 16 n E 21 n E 22 n E 23 n E 24 n E 25 n E 26 n E 31 n E 32 n E 33 n E 34 n E 35 n E 36 n E 41 n E 42 n E 43 n E 44 n E 45 n E 46 n E 51 n E 52 n E 53 n E 54 n E 55 n E 56 n E 61 n E 62 n E 63 n E 64 n E 65 n E 66 n ] , X n = [ A n 1 B n 1 C n 1 A n 2 B n 2 C n 2 , E i = [ E 1 E 2 E 3 E 4 E 5 E 6

Multiply exp ( − i s θ m ) by both ends of Formula (12), and integrate on interval ( − π , π ) , get infinite algebraic equations as follows:

∑ n = − ∞ ∞ E n s X n = E s (13)

Of which, E n s = 1 2 π ∫ − π π E n exp ( − i s θ j ) d θ j , E s = 1 2 π ∫ − π π E i exp ( − i s θ j ) d θ j .

According to the definition of dynamic stress concentration in the opening, the dynamic stress concentration coefficient is the ratio between the annular dynamic stress around the opening and the annular stress amplitude of the incident wave in the direction of incidence, i.e.

T θ z = i k 2 ( 1 + λ ) ∑ n = 0 ∞ { [ ε n η m ρ m Ω ′ ( η m ) | Ω ′ ( η m ) | J n − 1 ( k | Ω ( η m ) | ) { Ω ( η m ) | Ω ( η m ) | } n − 1 − η m ¯ ρ m Ω ′ ( η m ) ¯ | Ω ′ ( η m ) | J n + 1 ( k | Ω ( η m ) | ) { Ω ( η m ) | Ω ( η m ) | } n + 1 ] + ∑ m = 1 2 A n m [ η m ρ m Ω ′ ( η m ) | Ω ′ ( η m ) | H n − 1 ( 1 ) ( k | Ω ( η m ) | ) { Ω ( η m ) | Ω ( η m ) | } n − 1 + η m ¯ ρ m Ω ′ ( η m ) ¯ | Ω ′ ( η m ) | H n + 1 ( 1 ) ( k | Ω ( η m ) | ) { Ω ( η m ) | Ω ( η m ) | } n + 1 ] } + i e 15 ∑ n = 0 ∞ ∑ m = 1 2 B n n k − n η m ¯ ρ m Ω ′ ( η m ) ¯ | Ω ′ ( η m ) | ( Ω ( η m ) ¯ ) − n − 1 (14)

Among it, the comparison between maximum stress factor and static in the formula is the largest dynamic stress concentration factor.

With steady wave w ( i ) incidence along x axis, For Circular aperture mapping function with radius a 1 = a 2 = a , it can be taken as:

Ω = a η (15)

According to the formula of elastic wave scattering and dynamic stress concentration in piezoelectric materials with two holes, taking circular holes as an example, the corresponding calculation program is worked out, and n = 15, Poisson/v = 0.3, and dimensionless wave number Ka = 2.0 - 5.0.

This paper, based on elastic dynamics theory, uses Liu’s complex function and conformal mapping methods, and studies e-type piezoelectric material and elastic wave scattering and dynamic stress concentrations problems with double holes question. An analytical solution to the problems and numerical calculation results is given. By analysis and calculation results: 1) when the wave number 𝑘𝑎 is a constant, the maximum of the dynamic stress concentration factor varies as the hole spacing d/a varies, with different parameter 𝜆. Small hole spacing corresponds to greater impact on the dynamic stress concentration factor, and large hole spacing corresponds to less impact on the dynamic stress concentration factor. When the hole spacing d/a = 12.0, the dynamic stress concentration factor is no longer influenced by the adjacent holes, it is almost the same as the impact of a single object; 2) when the hole spacing d/a is a constant, as the wave numberkaincreases,the impact of the hole spacing on the dynamic stress concentration factor grows.

The authors declare no conflicts of interest regarding the publication of this paper.

Li, Z.H., Liu, H.Y. and Zhen, W.Z. (2020) Elastic Wave Scattering and Dynamic Stress Concentrations around Double Holes in Piezoelectric Media. Journal of Applied Mathematics and Physics, 8, 3060-3069. https://doi.org/10.4236/jamp.2020.812224