SPATIAL PROBLEM OF INTERACTION P (LONGITUDINAL) -WAVES ON A CYLINDRICAL CAVITY IN AN ELASTIC
MEDIUM
1N.U. Kuldashov, 2M. Choriev, 3U.A.Urolov
1,2,3 Tashkent Institute of Chemical Technology https://doi.org/10.5281/zenodo.7965947
Abstract. The spatial problem of diffraction of elastic waves, the angle of incidence relative to the longitudinal axis of the cavity, is considered. The aim of the study is the diffraction of elastic harmonic waves in a cavity depending on the angle of incidence. An analytical solution of the problem of diffraction of elastic waves in a cavity is obtained, which makes it possible to determine the force factors of the stress-strain state of the cavityfrom the angle of incidence. The environment is described by the Helmholtz equation and the solution of which is described by the Bessel and Neumann functions with complex arguments. The stress and permeability of the cavity point take on a maximum value in the region of long waves. A technique and algorithm for solving the spatial problem of diffraction of hormone waves in a cavity has been developed.
Keywords: cavity, wave, wave equation, wave incidence angle, Bessel and Neumann function.
1. Introduction. Many scientists dealt with the problems of diffraction of plane harmonic waves in homogeneous isotropic elastic bodies [1,2,3]. When elastic waves interact with inhomogeneities of the medium or with obstacles, the wave field changes to some extent. Due to the interaction of the wave with the interface, reflected and transmitted plane waves appear, depending on the angle of incidence. The incident elastic wave is considered to be flat, the angle of incidence of the wave front relative to the longitudinal axis is considered in [5,7]. The diffraction of elastic waves in inhomogeneous cylindrical elastic bodies placed in a deformable medium is the subject of works [4, 8].
The papers [8, 9] consider the problem of elastic wave diffraction in a viscoelastic cylindrical body. At the contact of the cylindrical body with the medium, the conditions of rigid contact are set. At infinity, the Somerfeld radiation conditions are set.
Scattering of elastic waves in an ideal infinite circular cavity is the simplest diffraction problem with an exact solution. In addition, problems that have exact solutions are the problems of scattering of acoustic waves in a sphere, an ellipsoid, and in other bodies whose surfaces are the coordinate surfaces of the corresponding curvilinear coordinate systems.
In this paper, we consider the diffraction of elastic harmonic waves in a cavity, which depends on the angle of incidence. Numerical results are obtained
2. Method
2.1. Problem statement and solution technique
In a cylindrical coordinate r,d, z system, the propagation of longitudinal P (incident) and PP (reflected), transverse PS (reflected along r,6) and PZ (reflected along z ) waves is considered. The design scheme is shown in (ras1)
Fig.1.
When the longitudinal P waves fall. The linear equation of motion of the medium is determined by a system of differential equations in partial derivatives, expressed in terms of displacements, by equations (Lame equations) [10]:
d2U
jlN2u + (/I + 2ju)graddiSu = p-
dt
2 '
(1)
À =
vE
vE
(1 + v)(l - 2v ) 2(1 + v )'
(2)
Here u{ux,u _) is the displacement vector of the medium points; p- medium density;
V - Poisson's ratio; , X, / - Lame coefficient.
We represent the medium displacement vector as a potential
u = gradcp + rotij/,!//(!//. ,<//0,i//z) (3)
(3) substituting (1) we obtained the following differential of the medium equation by the Helmholtz equations [9]:
* c2P d t2
0;
V Vz
l a2Wz
c2 at¿
0;
(4)
V2We
Voo + 2 dV
r
vw-h +
r- dO 2 dWo
1 d!wo
C dt
1 d2Wr
= 0 ;
r2 dO c2 dt2
= 0;
V2 =
d2 1 d 1 d2
+
+
+
d2
dr2 r dr r2 dO2 dz2 C , C - longitudinal and transverse velocity of wave propagation in the medium.
The boundary condition is set on the free surface of the cavity
r
+ Gípp) = G(rp) + G(pp) = G(P) + G(Pp) = 0
rr rr
rz rz
(5)
Solutions to equation (4) are sought in the form
Pr = e
Vz =e
, „ wcosO0 i (nO--0 z-wt )
„ _ wcosO0 i (nd+-0 z - wt )
/ \ . _ wcosft,
V ' i ( n°+——0 z-wt )
V
= e
P( r ) V( r) V(r )
(6)
r J
(6) we substitute into (4), as a result we obtain the equations of the medium in the form of the Bessel equation [6]:
ô p(r)+1 dp(r) +
dr r dr
C 2 2
& /1 2 n —(i - cos2 ^o) 2
vS r
p(r ) = 0
0 Vz(r)+i ôVz(r) +
^ 2
& 2 X n
ôr r ôr a2Vr,g(r) i ôVr,o(r)
J
2
—(i - cos Yo) —2
v<2
Vz (r) = 0
J
ôr2
r ôr
2
&(i - cos2 Y0) - -?(n2 + i) v CS r2
v ,o (r) = 0
J
Where
r
a =
, \
& 2 ^ n"
—(i - cos2 $0) - —
v S r J
ß=
Í 2
& 2 . n
—(i-cos Y0) —2 r
v<2
2 2
Q =
C 2 &
i
\
— (i - cos2 Y0) j(n + i)
vC2
Corresponding wave numbers of longitudinal and transverse waves. The general solution (7) in potentials has the following form.
m i (nO--°z-œt )
(ar )
cp{pp) = AsmO,- ^
n=0
.. _ &cosO0 X i ( nO+-°z-fflt )
Hi^)(ar )
n=0
Vps) ^ f^
vV J
( Pz )
v C J
cos
, &cosY0 . / „(2) X i ( nO+-/J0z-&t )! \
n=0
Hi2\Qr )
H:\ßr )
V n
<?(p>- incident expansion wave, <p(pp>- reflected expansion wave, Ips>,lpz> - reflected shear wave along the Z, 9, r axis;
à pp)
J
(ps) pz)
(7)
(8)
(9)
c
c
(pQ - amplitude of longitudinal incident waves; Oc, y0 - angle of incident and reflected waves; Jn (ar ), H (n2\ar ), H<2)(Q r ), H <2)( ßr ) - Bessel and Hankel function of the second kind of the n-th order; A, B, C - unknown coefficients which is determined from the boundary conditions (5). When the value of the argument (7) increases more than Pi, then the asymptotic value is used to calculate the value of the Bessel and Hankel function. At infinity, the Sommerfeld radiation conditions for the potentials of reflected waves must be satisfied [8]. Voltage vector.
1
2ß
G
d2 ç
dç
1 -v d2 ç v 1 - 2v dr2 r2(1 - 2v) dO2 + r(1 - 2v) ~dr
+
+
v
1 - 2v dz2 1 r 2 d2ç
d Ç v 1 d V 1 dVz
r drdO 2 dç
r1 dO
d Vo
drdz
i2
2^
d2Vr dr dz 1
r dr dO 1 dVr.
1 d vz d vz 1 dvz 1 d vo
r2 dO r2 dO1 dr2 r dr r dOdz
(10)
r dr
2ß
G,
v
d2ç
v
v d2ç --- + -
1 - 2v dr2 r2(1 - 2v) dO2 r(1 - 2v) dr drdz
dÇ + d 2Vo
1 dVo 1 d 2Vr
r dz 1
2^
+
goo =
r dOdz v d2 ç
1 - v d2ç 1 - v dç
1 - 2v dr2
,2 „ 1 -i2
r2(1 - 2v) dO2 r(1 - 2v) dr
+
v d ç 1 d v 1 dv2 1 dye 1 d v .
1 - 2v dz2 r drdQ r2 36 r dz r dzdQ' Since the spatial problem is considered, the angle of incidence, reflection of waves according to the Fermat, Hugenz and Snell's laws is the law of geolatric seismology. Can be generalized to rays vertically inhomogeneous media [5,7]
sinO0 _ sin^0
Cpo c\
Cs
/0 = arcsin(—^sinO )
CPo :
Substituting (9) into (10) we obtain the value
1 -Grr) =çcsinO(^ 1 v
2ß
a2j'i(a r)
\
va
r (1 - 2v)
Jn(ar ) -
1 - 2v y
v co2 cos2 O
+
v(-n2)
1 - 2v c
r 2(1 - 2v) 0 Jn (ar)
J„ (a r) +
Po
/
2j
-G (pp) = A ■ sinft
1 -v 1 - 2v
va
a
2- Hl»" (ar ) + •
<-"2) (1 - 2v)
H„<2) (ar)
+
V
r
(1 - 2v)
H»2" (ar) -„ V ^f^" H„<2) (ar)
(1 - 2v)
c„
Po
-Bc°s2r0 ■ — Hn{2)(Qr) + c■ cosYq 1 ■ i■ n ■ß■ Hn(2)(ßr)-1 ■ i ■ n ■ Hn(2)(ßr)
1
1
1 r
—G„a(p) = a0 sin6>0
2j
1 (2 1 Gftpp) = A ■ sin& 2
r 1 1
— ■ i ■ n -a- Jn' (ar) —- ■ i ■ n --Jn (ar)
Vr r
(11)
2j
i ■ n ■ a ■ Hn(2)/ (ar) - — ■ i ■ n ■ Hn(2) (ar)
-B ■ cos2 r0 ■ n ■— ■ Hn(2) (Qr) + C ■ cosr
f
V r
1- ■ (-n) ■ Hn^ {ßr) - p2 ■ Hn^ {ßr) +1 ■ ß ■ H :n(2)/ (ßr) 1
G (p) = %sinft
2ju
V
V
1 - 2v
r
.2 t //
a2 ■ Jn" (ar) +
(1 -2v)
\-n2) ■ Jn(ar )
•(1 -2v)
a
Jj (ar)
2j
Grz (pp) = A ■ sin 6>0
^H™(ar) + pV)-(-n2) ■
H(2) (ar )■
+B ■ cos Y0 ■ 1 -
V
"0 - 2v)
a H„(2)// (ar)
+
' —coY ■ Q ■ Hn ^ (Qr) +1 ■ —coY ■ Hn ^ (Qr)
V CSO
2 j
Gee = %sinfto
V a2 ■ Jn" (ar) +
r CSo
1 - v
V
1 - 2v
v 1 -v
+—;-ra^ Jn (ar) +
c 2 2 ^ \ — ■cos ft
•(1 - 2v)
(1 - 2v)
(1 - 2v) JÄar )
(-n2 )■ Jn(ar ) +
po J
(12)
1
so
2ß
GJpp) = A • sinJ
1 iœ
1 -v
1 -v
-^^a2 • H (2)//(ar) + 2; , • (-n2)• H (2)(ar) + 1-2v H V } r2 (l-2v) n V }
1 - v
1-^-a-H (2)/(ar) + • (l - 2v) n V } 1-2v
f 2 2 A œ cos Yo
Hn (ar )
so y
- B • cos2 Yo-----Hn (2) (ar ) + C • cos y0
r c
i • n
ß^ Hn12" (ßr )
Substituting stress (11) into boundary conditions (5), we obtain a system of complex algebraic equations [3x3] in the form
vT r
(13)
[C](ABC)'=[P] ;
[C] is a quadratic complex matrix whose elements have the following form:
Cu = sin J
—a2 • Hn(2)// (ar) + V( n \ H(2) (ar) + , a x Hn(2)/ (ar) 1 - 2v n V } r2 (1 - 2v) n V } r(1 - 2v) n V r
'H )
v (O2 cos2 J
(1 - 2v)
■H, <2> (ar )
Po
Cl2 = cos2 ro • Hfï (ar );
c
r
C13 = cos Yo
\
C21 = sinJ0
- • i • n •ß• Hn(2)/ (ßr)--• i • n • Hj2)(ßr)
Vr r y
/ 2 2 A
- • i • n-a-Hn (2)/ (ar )—- • i • n • Hn (2) (ar )
V
r
r
y
C- = cos2 ro • n ( H™ (ar);
c
C23 = cos Yo • f 4r • (-n) • (ßr ) - p2 • Hn^ (ßr ) +1 ß. 2)/ (ßr )
V r
M-v
C31 = sinJo
———a2 • Hn(2)// (ar) + -- • (-n2) • Hn(2) (ar) +
1 - 2v " V ' r2 (1 - 2v) K ) n y >
v
a Hn^ {ar)
v r.
C32 = cosYo
(1 - 2v)
i(cosYo • a • Hn (2)/ (ar )+1 • iœcosYo • Hn (2) (ar )
V Cso
; C33 = o;
[ P] - Vector column of external loads
1
P =çosinOo
1 -v 1 - 2v
aJ'l(ar )
v(-n2) „ ,
+ ^-— J, (ar ) +
r2(1 - 2v) n
va
r (1 - 2v)
J'n(ar )■
v <c2 cos2 O
1 - 2v c 2
Jn (ar );
Po
( 1 1 P = çsinO0 --i-n-a-Jn'{ar) —-•i-n-Jn{ar)
P3 =çosinOo
v 2 r // / \ v a2- Jj {ar) +
1 - 2v
{1 - 2v)
{-n2 )- Jn {ar )
v
•f1 - 2v)
a
Jj {ar)
Solving the complex algebraic equations (13) we find the unknown complex coefficients A, B, C, substituting these coefficients into (11.12) we find the corresponding stresses. 3. Results and analysis
Numerical results and graphs obtained using C++ software
0.59 1.09 1.59 2.09 2.59
Fig 2.
Fig 3.
Fig 6.
Fig 7.
In fig 1, 2, 3 Figure 3 shows the hoop stresses of the medium cavity when soft soil is considered; the maximum hoop stresses at pi/3 reaches 6.54 dimensionless units, at pi/4 it reaches 5.53 dimensionless units, at pi/6 it reaches 6.26 dimensionless units. When sandstone is considered, the maximum hoop stresses at pi/3 reaches 10.27 dimensionless units, at pi/4 it reaches 6.39 dimensionless units, at pi/6 it reaches 3.49 dimensionless units. In three cases, it can be seen that when a longitudinal wave falls with an increase in the angle of incidence, the hoop stress of the cavity also increases
4. Conclusions
1. A technique and algorithm have been developed for solving the spatial problem of diffraction of harmonic waves in a cavity, an elastic medium.
2. when a longitudinal wave falls with an increase in the angle of incidence relative to the longitudinal axis of the cavity, the hoop stress increases.
REFERENCES
1. A. N. Guz, V. D. Kubenko, and M. A. Cherevko, Diffraction of elastic waves. Kyiv: Naukova Dumka, 1978. 308 p.
2. A. N. Guz, Elastic waves in bodies with initial stresses. T. 2. Patterns of distribution. Kyiv: Naukova Dumka, 1986. 536 p.
3. Baron M.L. Parnas R., Disspacements is and velocities produced by the diffraction of a pressure wave by a cylindrical cavity in an elastic medium. Trans ASME. J. Appl. Mech., 1962.29 No. 2 385-395
4. Safarov I.I., Kuldashov N.U., Boltaev Z.I. Oscillations and Waves in a Layered Homogeneous Viscoelastic Medium. International Journal of Emerging Engineering Research and Technology. Volume 6, Issue 3, 2018, PP 27-32.
5. Landau L.D., Lifshitz E.M. The Theory of Elasticity. - Moscow: Science. Ch. ed. Phys.-Math. literature, 1965. - 204 p.
6. Gray, E. Bessel functions and their applications to physics and mechanics. - M .: Publishing house of foreign. literature, 1953. - 372 p.125
7. K. Aki and P.G. Richards, 1980. Quantitative Seismology, Theory and Methods. Volume I: 557 pages
8. Grinchenko V.T., Meleshko V.V. Harmonic oscillations and waves in elastic bodies. - K .: Nauk. Dumka, 1981. - 284 p.
9. Safarov I.I., Teshaev M., Toshmatov E., Boltaev Z. and Homidov.F.F. Torsional vibrations of a cylindrical shell in a linear viscoelastic medium / IOP Conference Series: Materials Science and Engineering 883 (1) 012190 (2020).
10. Teshaev M.K, Safarov I.I, Kuldashov N.U, Ishmamatov M.R and Ruziev T.R, On the Distribution of Free Waves on the Surface of a Viscoelastic Cylindrical, Cavity / Journal of Vibrational Engineering and Technologies 8 (4) pp 579-585. (2020)