Properties of Wave Motion in a Cylindrical Shell, Interacting with Viscous Liquid

Full-Text HTML XML Download as PDF (Size:368KB) PP. 1-23
The propagation of natural waves in a cylindrical shell (elastic or viscoelastic) that is in contact with a viscous liquid is considered. The problem reduces to solving spectral problems with a complex incoming parameter. The system of ordinary differential equations is solved numerically, using the method of orthogonal rotation of Godunov with a combination of the Muller method. The dissipative processes in the mechanical system are investigated. A mechanical effect is obtained that describes the intensive flow of mechanical energy.

1. Introduction

Quantitative data and some qualitative analysis of the dispersion properties of normal waves in a cylinder with a liquid were carried out in [1] [2] [3] [4] . In [5] [6] [7] [8] , data were obtained for cylinders filled with a liquid of solid and soft materials. It is studied in a rather wide range of frequencies and wave numbers. The effect of the wall thickness of the cylinder on the dispersion properties of normal waves in the range of variation of the relative internal radiation is from 0.3 to 0.99 [9] [10] . This article focuses on the dynamic behavior of a cylindrical shell (elastic or viscoelastic), which contains a viscous liquid. The problem of wave propagation in a filled or submerged liquid with a cylindrical shell is of great practical importance. The phenomenon of undulating fluid motion in elastic cylindrical shells has attracted the attention of many researchers [11] - [16] . In these papers, wave processes in an elastic cylindrical shell containing an ideal fluid are considered. Also, the refined classical shell equations, the average fluid flow (or gas) density, the influence of radial and longitudinal inertia forces are used. In [17] [18] [19] , the laws of wave processes in an elastic shell with a viscous fluid are analyzed in the model of linear equations of hydrodynamics of a viscous compressible fluid. Unlike other systems, a cylindrical shell (elastic or viscoelastic) and liquid (ideal or viscous) are considered as an inhomogeneous dissipative mechanical system [20] [21] [22] .

2. Statement of the Problem

An infinite length of deformable (viscoelastic) cylindrical shell of radius R with constant thickness ${h}_{0}$ , density ${\rho }_{0}$ , Poisson’s ratio ${\nu }_{0}$ , filled with a viscous fluid with density at equilibrium. Fluctuations of a shell under a load, the density of which is denoted p1, p2, pn respectively, can be described by following [11] [12] [14] , equations:

$Lu-{L}_{0}\underset{0}{\overset{t}{\int }}{R}_{L}\left(t-\tau \right)Lu\left(\tau \right)\text{d}\tau =\frac{\left(1-{\nu }_{0}^{2}\right)}{{E}_{0}{h}_{0}}р+{\rho }_{0}\frac{\left(1-{\nu }_{0}^{2}\right)}{{E}_{0}}\left(\frac{{\partial }^{2}u}{\partial {t}^{2}}\right),$ (1)

Here $u=u\left({u}_{r},{u}_{\theta },{u}_{z}\right)$ ―displacement vector points of the middle surface of the shell and membranes for Kirchhoff-Love it has a dimension equal to three $\left({u}_{r}=u;{u}_{\theta }=v;{u}_{z}=w\right)$ , and to membranes such as the dimension of Timoshenko u is five. Here, in addition to the axial, circumferential and normal movements added more angles of rotation normal to the middle surface in the axial and circumferential directions [22] ; ${\left\{u\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}v\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}w\right\}}^{\text{T}}$ ―the displacement vector with axial, radial and circumferential components, respectively (“+” sign in front of pn and the sign “−” before the last component of the inertial member says that is considered positive motion towards the center of curvature); ${R}_{L}\left(t-\tau \right)$ ―the core of relaxation; ${L}_{0}$ ―instantaneous modulus of elasticity.

The amplitudes of the oscillations are considered small, which allows you to record the basic relations in the framework of the linear theory. The system of linear equations of motion of a viscous barotropic liquid can be written as [22] :

$\frac{\partial \vartheta }{\partial t}-{\nu }^{\ast }\Delta \vartheta +\frac{1}{{\rho }_{0}^{\ast }}grad\text{\hspace{0.17em}}P-\frac{{\nu }^{\ast }}{3}grad\text{\hspace{0.17em}}div\text{\hspace{0.17em}}\vartheta =0$

$\begin{array}{l}\frac{1}{{\rho }_{0}^{\ast }}\frac{\partial {\rho }^{\ast }}{\partial t}+div\text{\hspace{0.17em}}\vartheta =0;\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\frac{\partial P}{\partial {\rho }^{\ast }}={a}_{0}^{2},\text{\hspace{0.17em}}{a}_{0}=const.\\ {\stackrel{˙}{u}}_{z}={\vartheta }_{z},\text{\hspace{0.17em}}{\stackrel{˙}{u}}_{r}={\vartheta }_{r},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\stackrel{˙}{u}}_{\theta }={\vartheta }_{\theta },\\ {q}_{z}=-{p}_{rz},\text{\hspace{0.17em}}{q}_{r}=-{p}_{r},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{q}_{\theta }=-{p}_{r\theta }.\end{array}$ (2)

$\begin{array}{l}{p}_{rz}={\mu }^{\ast }\left(\frac{\partial {\vartheta }_{z}}{\partial r}+\frac{\partial {\vartheta }_{r}}{\partial z}\right);\\ {p}_{rr}=-p+{\lambda }^{\ast }\left(\frac{\partial {\vartheta }_{r}}{\partial r}+\frac{\partial {\vartheta }_{z}}{\partial z}+\frac{{\vartheta }_{r}}{r}\right)+2{\mu }^{\ast }\frac{\partial {\vartheta }_{r}}{\partial r};\\ {p}_{r\theta }={\mu }^{\ast }\left(\frac{1}{r}\frac{\partial {\vartheta }_{z}}{\partial \theta }+\frac{\partial {\vartheta }_{\theta }}{\partial r}-\frac{{\vartheta }_{\theta }}{r}\right).\end{array}$

Here, in the Equations (2) $\vartheta =\vartheta \left({\vartheta }_{r},{\vartheta }_{\theta },{\vartheta }_{z}\right)$ ―the velocity vector of fluid particles; ${\rho }^{\ast }$ and Р―disturbance density and fluid pressure; ${\rho }_{0}^{\ast }$ and а0―density and sound velocity in the fluid at rest; ${\nu }^{\ast },{\mu }^{\ast }$ ―kinematic and dynamic viscosity;

for the second viscosity coefficient ${\lambda }^{\ast }$ accepted ratio $\lambda =-\frac{2}{3}{\mu }^{•}$ ; ${р}_{rz},{р}_{rr},{р}_{r\theta }$

components of the stress tensor in the fluid. Equation (1), respectively, kinematic and dynamic boundary conditions, which, because of the thin-walled shell, we will meet on the middle surface (r = R). Equations (1) and (2) is a closed system of relations hydro visco elastic cylindrical shell for containing a viscous compressible fluid. This is for shell obeying Kirchhoff-Love hypotheses. Be investigated joint shell and liquid fluctuations, harmonic of the axial coordinate z and decay exponentially over time, or time-harmonic and damped with respect to z.

3. Method of Solution

We accept the integral terms in (1) small, then the function

$u\left(r,t\right)={u}_{1}\left(r,t\right){\text{e}}^{-i{\omega }_{R}t}$ ,

where $\psi \left(r,t\right)$ ―slowly varying function of time, ${\omega }_{R}$ ―real constant. Next, using the procedure of freezing [23] , then the integral-differential Equation (1) takes the form

$L\left[1-{\Gamma }^{C}\left({\omega }_{R}\right)-i{\Gamma }^{S}\left({\omega }_{R}\right)\right]{u}_{1}=\frac{\left(1-{\nu }_{0}^{2}\right)}{{E}_{0}{h}_{0}}р+{\rho }_{0}\frac{\left(1-{\nu }_{0}^{2}\right)}{{E}_{0}}\left(\frac{{\partial }^{2}{u}_{1}}{\partial {t}^{2}}\right),$ (3)

where, for shell Kirchhoff-Love

$L=\left(\begin{array}{ccc}\frac{{\partial }^{2}}{\partial {x}^{2}}+\frac{1-{\nu }_{0}}{2{R}^{2}}\frac{{\partial }^{2}}{\partial {\phi }^{2}}& \frac{1+{\nu }_{0}}{2R}\frac{{\partial }^{2}}{\partial x\partial \phi }& \frac{{\nu }_{0}}{R}\frac{\partial }{\partial x}\\ \frac{1+{\nu }_{0}}{2R}\frac{{\partial }^{2}}{\partial x\partial \phi }& \frac{1+\nu {}_{0}}{2}\left(1+4a\right)\frac{{\partial }^{2}}{\partial {x}^{2}}+\left(1+a\right)\frac{{\partial }^{2}}{\partial {\phi }^{2}}& \frac{1}{{R}^{2}}\frac{\partial }{\partial \phi }-a\left(2-{\nu }_{0}\right)\frac{{\partial }^{3}}{\partial {x}^{2}\partial \phi }-\frac{a}{{R}^{2}}\frac{{\partial }^{3}}{\partial {\phi }^{3}}\\ \frac{{\nu }_{0}}{R}\frac{\partial }{\partial x}& \frac{1}{{R}^{2}}\frac{\partial }{\partial \phi }-a\left(2-{\nu }_{0}\right)\frac{{\partial }^{3}}{\partial {x}^{2}\partial \phi }-\frac{a}{{R}^{2}}\frac{{\partial }^{3}}{\partial {\phi }^{3}}& \frac{1}{{R}^{2}}+a{\left(\frac{{\partial }^{2}}{\partial {x}^{2}}+\frac{1}{{R}^{2}}\frac{{\partial }^{2}}{\partial {\phi }^{2}}\right)}^{2}\end{array}\right),$

${\Gamma }^{C}\left({\omega }_{R}\right)=\underset{0}{\overset{\infty }{\int }}{R}_{L}\left(\tau \right)\mathrm{cos}{\omega }_{R}\tau \text{d}\tau$ , ${\Gamma }^{S}\left({\omega }_{R}\right)=\underset{0}{\overset{\infty }{\int }}{R}_{L}\left(\tau \right)\mathrm{sin}{\omega }_{R}\tau \text{d}\tau$ respectively, cosine

and sine Fourier transforms relaxation kernel material. As an example, the viscoelastic material take three parametric kernel relaxation ${R}_{L}\left(t\right)=A{\text{e}}^{-\beta t}/{t}^{1-\alpha }$ , ρ―material density shell; E―Young’s modulus; ν―Poisson’s ratio, $a={h}^{2}/12{R}^{2}$ Let’s move on to the dimensionless axial coordinate $\xi =x/R$ and multiply by R2 system (3). The matrix of the resulting system will take the form

$L=\left(\begin{array}{ccc}\frac{{\partial }^{2}}{\partial {\xi }^{2}}+\frac{1-{\nu }_{0}}{2}\frac{{\partial }^{2}}{\partial {\phi }^{2}}& \frac{1+{\nu }_{0}}{2}\frac{{\partial }^{2}}{\partial \xi \partial \phi }& {\nu }_{0}\frac{\partial }{\partial \xi }\\ \frac{1+{\nu }_{0}}{2}\frac{{\partial }^{2}}{\partial \xi \partial \phi }& \frac{1-{\nu }_{0}}{2}\left(1+4a\right)\frac{{\partial }^{2}}{\partial {\xi }^{2}}+\left(1+a\right)\frac{{\partial }^{2}}{\partial {\phi }^{2}}& \frac{\partial }{\partial \phi }-a\left(2-{\nu }_{0}\right)\frac{{\partial }^{3}}{\partial {\xi }^{2}\partial \phi }-a\frac{{\partial }^{3}}{\partial {\phi }^{3}}\\ {\nu }_{0}\frac{\partial }{\partial \xi }& \frac{\partial }{\partial \phi }-a\left(2-{\nu }_{0}\right)\frac{{\partial }^{3}}{\partial {\xi }^{2}\partial \phi }-\frac{a}{{R}^{2}}\frac{{\partial }^{3}}{\partial {\phi }^{3}}& \frac{1}{{R}^{2}}+a{\left(\frac{{\partial }^{2}}{\partial {\xi }^{2}}+\frac{1}{{R}^{2}}\frac{{\partial }^{2}}{\partial {\phi }^{2}}\right)}^{2}\end{array}\right).$ (4)

Expanding Equation (2) and (3) in coordinate form, it is easy to see that the relations (2) - (3) break up into independent boundary value problems:

- Torsional vibrations:

$\begin{array}{l}\frac{\partial {p}_{r\theta }}{\partial r}+\frac{2{p}_{r\theta }}{r}+\frac{\partial {p}_{\theta z}}{\partial z}={\rho }_{0}^{\ast }{\stackrel{¨}{\vartheta }}_{\theta },\\ {p}_{r\theta }={\mu }^{\ast }\left(\frac{\partial {\vartheta }_{\theta }}{\partial r}-\frac{{\vartheta }_{\theta }}{r}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{p}_{\theta z}={\mu }^{\ast }\frac{\partial {\vartheta }_{\theta }}{\partial r},\\ r={R}_{1}:\stackrel{¯}{G}{h}_{0}\frac{{\partial }^{2}{u}_{\theta }}{\partial {z}^{2}}-\left({\rho }_{0}h{\stackrel{¨}{u}}_{\theta }±{\sigma }_{\phi r}\right)=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\stackrel{¯}{G}=\frac{\stackrel{¯}{E}}{2\left(1+{v}_{0}\right)},\\ r=0:{p}_{r\theta }=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\stackrel{¯}{E}={Е}_{0}\left(1-{\Gamma }^{C}\left({\omega }_{R}\right)-i{\Gamma }^{S}\left({\omega }_{R}\right)\right)\end{array}$ (5)

- Longitudinal transverse vibrations:

$\begin{array}{l}\frac{\partial {р}_{rr}}{\partial r}+\frac{{p}_{rr}-{p}_{\theta \theta }}{r}+\frac{\partial {p}_{rz}}{\partial z}={\rho }_{0}^{\ast }\frac{\partial {\vartheta }_{r}}{\partial t}\\ \frac{\partial {p}_{rz}}{\partial r}+\frac{{p}_{rz}}{r}+\frac{\partial {p}_{zz}}{\partial z}={\rho }_{0}^{\ast }\frac{\partial {\vartheta }_{z}}{\partial t}\end{array}$

$\begin{array}{l}{p}_{rr}=-p+{\lambda }^{\ast }{k}_{\eta }div\text{\hspace{0.17em}}\vartheta +2{\mu }^{\ast }\frac{\partial {\vartheta }_{r}}{\partial r},\\ {p}_{\theta \theta }=-p+{\lambda }^{\ast }div\text{\hspace{0.17em}}\vartheta +2{\mu }^{\ast }\frac{{\vartheta }_{r}}{r}\\ {p}_{zz}=-p+{\lambda }^{\ast }div\text{\hspace{0.17em}}\vartheta +2{\mu }^{\ast }\frac{\partial {\vartheta }_{z}}{\partial z}\end{array}$ (6)

$\begin{array}{l}{р}_{rz}={\mu }^{\ast }\left(\frac{\partial {\vartheta }_{z}}{\partial r}+\frac{\partial {\vartheta }_{r}}{\partial z}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{р}_{r\theta }={\mu }^{\ast }\left(\frac{1}{r}\frac{\partial {\vartheta }_{z}}{\partial \theta }+\frac{\partial {\vartheta }_{\theta }}{\partial r}+\frac{{\vartheta }_{r}}{r}\right)\\ \frac{\partial {\rho }^{\ast }}{\partial t}+{\rho }_{0}div\vartheta =0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}div\text{\hspace{0.17em}}\vartheta =\frac{\partial {\vartheta }_{r}}{\partial r}+\frac{{\vartheta }_{r}}{r}+\frac{\partial {\vartheta }_{z}}{\partial z},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\frac{\partial p}{\partial \rho }={а}_{0}^{2}\end{array}$

$\begin{array}{l}r={R}_{1}:\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\stackrel{¯}{D}\frac{{\partial }^{4}{u}_{r}}{\partial {z}^{4}}+\frac{\stackrel{¯}{C}}{{R}_{1}}\left(\frac{{u}_{r}}{{R}_{1}}+{\nu }_{0}\frac{\partial {u}_{z}}{\partial z}\right)+{p}_{rr}+{\rho }_{0}{h}_{0}\frac{{\partial }^{2}{u}_{r}}{\partial {t}^{2}}=0,\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\stackrel{¯}{C}\left(\frac{{\partial }^{2}{u}_{r}}{\partial {z}^{2}}+\frac{{v}_{0}}{R}\frac{\partial {u}_{r}}{\partial z}\right)-\left({p}_{rz}±{\rho }_{0}{h}_{0}\frac{{\partial }^{2}{u}_{z}}{\partial {t}^{2}}\right)=0,\end{array}$

Let the wave process is periodic in z and fades over time, then is given a real wave number k, and the complex frequency is the desired characteristic value. Solution of (2) - (6) for the major unknowns satisfying constraints imposed above the dependence on time and coordinates z, should be sought in the form [24]

$\begin{array}{l}{\left({р}_{rr},{р}_{rz},{р}_{r\theta },u,\vartheta \right)}^{\text{Т}}\\ =\underset{m}{\sum }{\left({\sigma }_{rm}\left(\xi ,\phi ,t\right),{\tau }_{zm}\left(\xi ,\phi ,t\right),{\tau }_{\phi m}\left(\xi ,\phi ,t\right),{u}_{m}\left(\xi ,\phi ,t\right),{\vartheta }_{m}\left(\xi ,\phi ,t\right)\right)}^{\text{Т}},\end{array}$ (7)

where ${u}_{m}\left(\xi ,\phi ,t\right)={u}_{m}{\left\{{U}_{m},{V}_{m},{W}_{m}\right\}}^{\text{T}},\text{\hspace{0.17em}}\vartheta \left(\xi ,\phi ,t\right)={\vartheta }_{m}{\left\{{\vartheta }_{rm},{\vartheta }_{\theta m},{\vartheta }_{zm}\right\}}^{\text{T}}$ .

Expressions (7) in the form

$\begin{array}{l}{\left({\sigma }_{rm}\left(\xi ,\phi ,t\right),{\tau }_{zm}\left(\xi ,\phi ,t\right),{\tau }_{\phi m}\left(\xi ,\phi ,t\right)\right)}^{\text{Т}}\\ ={\left({\sigma }_{r}\mathrm{cos}\left(m\phi \right),{\tau }_{z}\mathrm{cos}\left(m\phi \right),{\tau }_{\phi }\mathrm{sin}\left(m\phi \right)\right)}^{\text{T}}{\text{e}}^{iкz-i\omega t},\end{array}$

$\begin{array}{l}{\left({u}_{m}\left(\xi ,\varphi ,t\right),{\vartheta }_{m}\left(\xi ,\phi ,t\right)\right)}^{\text{Т}}\\ =\left({U}_{m}\mathrm{cos}\left(m\phi \right),{V}_{m}\mathrm{sin}\left(m\phi \right),{W}_{m}\mathrm{cos}\left(m\phi \right),\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{{\vartheta }_{r}\mathrm{cos}\left(m\phi \right),{\vartheta }_{\theta }\mathrm{cos}\left(m\phi \right),{\vartheta }_{z}\mathrm{cos}\left(m\phi \right)\right)}^{\text{T}}{\text{e}}^{iкz-i\omega t},\end{array}$

where ${\sigma }_{r},{\tau }_{z},{\tau }_{\phi },{U}_{m},{V}_{m},{W}_{m},{\vartheta }_{r},{\vartheta }_{\theta },{\vartheta }_{z}$ ―Amplitude integrated vector-function; к―wavy number; С―phase velocity; ω―complex frequency; m―circumferential wave number (the number of district-wave), takes values $m=1,2,3,\cdots$ . When m = 0, happening Ax symmetrical vibrations. This approach allows you to seek a solution for every fixed value of the wave number of the district m independently.

In this way С, k, ω it is well-known real and complex spectral parameters of the type of problem.

To elucidate their physical meaning consider two cases:

1) $к={к}_{R}$ ; $С={С}_{R}+i{C}_{i}$ , Then the solution of (5) has the form of a sine wave х, whose amplitude decays over time;

2) $к={к}_{R}+i{к}_{I}$ ; $С={С}_{R}$ , Then at each point х fluctuations established, but х attenuate. In the case of axially symmetric on the axis r = 0 conditions must be satisfied conditions ${р}_{r\theta }={р}_{rz}=0$ , ${\vartheta }_{r}=0$ . If the outer surface г = R assumed stationary, then ${u}_{r}={u}_{z}={u}_{\phi }=0$ . The superposition of the solutions (8) forms an exponentially decaying over time the standing wave that describes the natural oscillations of a liquid and a cylindrical shell of finite length with boundary conditions. With infinite length sheath similarly specified type of movement (8) will be called private or free fluctuations. In the case of steady-state over time and fading coordinate the process, in contrast, is a well-known real rate of ω, as desired be a complex wave number k. In contrast to their own, these fluctuations will be called the established. Actual values of the ω in the first case, and k, second frequency have the physical meaning of the process in time and the coordinate, respectively. Imaginary part―the rate of decay of wave processes in time and Z, respectively [25] . The value of 1/Imk sometimes defined as the interval damped wave propagation. In the extreme case, the elastic range spread endless. The degree of attenuation of wave process in the time period is characterized by the logarithmic decrement

${\delta }_{c}=2\text{π}|\mathrm{Im}\omega |/\mathrm{Re}\omega$ (8)

Decrement is similar to the spatial

${\delta }_{y}=2\text{π}|\mathrm{Im}k|/\mathrm{Re}k$ .

You can also introduce the concept of phase velocity of its own and steady motions

${c}_{c}=\frac{\mathrm{Re}\omega }{R},{c}_{y}=\frac{\omega }{\mathrm{Re}k}$

The values Сс and Су have physical sense speeds of zero state at its own and steady oscillations, respectively, and, in contrast to the elastic (real) case, do not coincide with each other at the same frequencies. Two types of oscillations (and set their own), you can put two different formulations of the problem. And in the non-stationary case, namely the Cauchy problem for an infinite shell and boundary value problem for the semi-infinite interval changes Z. In either case, the solution is using the integral transformation of the decisions of the respective steady-state problems. For example, in the case of the Cauchy problem, the main vector of unknowns ${\stackrel{¯}{Y}}^{c}$ . It can be in a superposition of waves

${\stackrel{¯}{Y}}^{c}=\left(r,z,t\right)=\underset{n}{\sum }\underset{-\infty }{\overset{\infty }{\int }}{Y}_{n}^{c}\left(r,k\right)\mathrm{exp}\left[t\left(kz-{\varpi }_{n}\left(k\right)t\right)\right]\text{d}k$ , (9)

where vectors ${\stackrel{¯}{Y}}_{n}^{c}$ are their own form of the problem of natural oscillations, normalized so that the spatial Fourier spectrum of the initial disturbance $\stackrel{¯}{f}\left(r,z\right)={\stackrel{¯}{Y}}^{c}\left(r,z,0\right)$ forms a linear combination

$\stackrel{¯}{f}\left(r,z\right)=\underset{-\infty }{\overset{\infty }{\int }}F\left(r,k\right){\text{e}}^{ikz}\text{d}k,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\stackrel{¯}{f}\left(r,k\right)=\underset{n}{\sum }{\stackrel{¯}{Y}}_{n}^{c}\left(r,k\right)$ . (10)

Similarly, the main vector of unknowns ${\stackrel{¯}{Y}}^{y}$ boundary value problem is calculated according to the expression

${\stackrel{¯}{Y}}^{y}\left(r,z,t\right)=\underset{n}{\sum }\underset{-\varpi }{\overset{\varpi }{\int }}{\stackrel{¯}{Y}}_{k}^{y}\left(r,\omega \right)\mathrm{exp}\left[ik\left(\omega \right)z-\omega t\right]\text{d}\omega$ (11)

where ${\stackrel{¯}{Y}}_{k}^{y}$ forms steady-state oscillation, the linear combination of which should form a Fourier spectrum given boundary perturbation

$\stackrel{¯}{q}\left(r,t\right)={\stackrel{¯}{Y}}^{y}\left(r,0,t\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}\stackrel{¯}{q}\left(r,t\right)=\underset{-\infty }{\overset{\infty }{\int }}q\left(r,\omega \right){\text{e}}^{-i\omega t}\text{d}\omega ,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\stackrel{¯}{q}\left(r,\omega \right)=\underset{n}{\sum }\underset{-\infty }{\overset{\infty }{\int }}{Y}_{n}^{y}\left(r,\omega \right)$

Obviously, the solutions (8) and (9) have a meaning only when there are (10) and (11). So, there are four possible variants of steady motions, which are discussed below, and established their own systems fluctuations shell―fluid inside and outside the sheath liquid [25] . Substituting the solution (7) in the system of differential Equations (2) - (6) we obtain a system of ordinary differential equations with complex coefficients, which is solved by Godunov’s orthogonal sweep method with a combination of method of Muller [26] in the complex arithmetic.

4. Torsional Vibrations

After performing in (5) the change of variables (7) permitting relations describing stationary torsional vibrations of the shell liquid, formulated in the form of the spectral boundary value problem for a system of two ordinary differential equations

$\begin{array}{l}\frac{\text{d}{\tau }_{\phi }}{\text{d}r}=-\left({\rho }_{0}^{\ast }{\omega }^{2}-i{\mu }^{\ast 2}{\xi }^{2}\omega \right){\vartheta }_{\theta }-\frac{2{\tau }_{\phi }}{r}\\ \frac{\text{d}v}{\text{d}r}=\frac{{\vartheta }_{\theta }}{r}+\frac{i}{\omega {\mu }^{\ast }}{\tau }_{\phi }\\ r={R}_{1}:{h}_{0}\left(\stackrel{¯}{G}{\xi }^{2}-\xi {\rho }_{0}{\omega }^{2}\right){\vartheta }_{\theta }±{\tau }_{\phi }=0\\ r=0:{\tau }_{\phi }=0\end{array}$ (12)

First investigate fluctuations of fluid in the walls. Equations (12) can be converted to a single equation for the displacement v

$\frac{{\text{d}}^{2}{\vartheta }_{\theta }}{\text{d}{r}^{2}}+\frac{\text{d}{\vartheta }_{\theta }}{r\text{d}r}+\left(-{\xi }^{2}+i\frac{\omega }{{v}^{\ast }}-\frac{1}{{r}^{2}}\right){\vartheta }_{\theta }=0;\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{v}^{\ast }=\frac{{\mu }^{\ast }}{{\rho }_{0}^{\ast }}$ (13)

The solution of Equation (13) is limited at r = 0 has the form

$v={A}_{1}{J}_{1}\left(r\sqrt{-{k}^{2}+i\frac{\omega }{{v}^{\ast }}}\right)=0$ . (14)

where J1-Bessel function of the first order, and A is an arbitrary constant. Given the immobility of the shell, we obtain the dispersion equation

${J}_{1}\left({R}_{1}\sqrt{-{k}^{2}+i\frac{\omega }{{v}^{\ast }}}\right)=0$ (15)

from whence

${\omega }_{n}=-i\left({v}^{\ast }{k}^{2}+{\Gamma }_{т}^{2}\right)$ (16)

in the case of natural oscillations and

${k}_{n}=\sqrt{-{\Gamma }_{n}^{2}+i\frac{\omega }{{v}^{\ast }}}$ (17)

in the case of steady-state oscillations. Here, through the Гn marked the roots of Bessel functions assigned to R. As it can be seen from (15), (16) own motion aperiodicity always on time, with the anchor points are fixed (the phase velocity С0 = 0), while the steady motion are oscillatory in nature, as the nodal point move at the speed of Су, a monotonically increasing from zero to indefinitely with an increase in viscosity or ${\nu }^{\ast }$ . These characteristic features of the motion of a viscous medium will appear in the following more complex example.

Let us now consider the relation (12) in the case of the internal arrangement of the liquid. This problem can be solved in the same way using special features and have a dispersion equation

$-{k}^{2}+\frac{{\omega }^{2}}{{a}^{2}}+\frac{\omega {v}^{\ast }}{{a}^{3}\stackrel{˜}{p}\stackrel{˜}{h}{R}^{2}}+\left(z\frac{{J}_{0}\left(z\right)}{{J}_{1}\left(z\right)}-2\right)=0$ (18)

which was first obtained in A. Guz [22] . Here we have introduced new designations

$\stackrel{˜}{p}=\frac{{\rho }^{\ast }}{{\rho }_{0}};\text{\hspace{0.17em}}\text{\hspace{0.17em}}\stackrel{˜}{h}=\frac{h}{{R}_{1}};\text{\hspace{0.17em}}\text{\hspace{0.17em}}z={R}_{1}\sqrt{-{k}^{2}+i\frac{\omega }{{v}^{\ast }}};\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}a=\sqrt{\frac{G}{{\rho }_{0}}}$

shear wave velocity shell: J0-Bessel function of zero order.

The direct solution of the Equation (18) comes up against certain difficulties caused by the need to calculate the Bessel functions of complex argument. Therefore we examine (18) by means of asymptotic representations of these functions at small and large arguments z. The smallness of z occurs in the low-frequency vibrations. According to the known expansion J0 and J1 power series

${J}_{0}=1-\frac{{z}^{2}}{4}=\cdots ;\text{\hspace{0.17em}}\text{\hspace{0.17em}}{J}_{1}\left(z\right)=\frac{z}{2}\left(1-\frac{{z}^{2}}{8}+\cdots \right);$ (19)

Hold the expansions (19) only the first term, we obtain

$-{k}^{2}+\frac{\omega }{{a}^{2}}=0$

dispersion equation of torsional vibrations or dry shell filled with an ideal liquid, keeping in (19) on the first two terms, we have the equation

$-{k}^{2}+\frac{{\omega }^{2}}{{a}^{2}}+i\frac{\omega {v}^{\ast }}{4{a}^{2}\stackrel{˜}{p}\stackrel{˜}{h}}\left({k}^{2}-i\frac{\omega }{{v}^{\ast }}\right)=0$ (20)

the root of which, for example, in the case of steady-state oscillations is given by

$k=\frac{\omega }{a}{\left[\left(1+\frac{1}{4\stackrel{˜}{p}\stackrel{˜}{h}}\right)/\left(1-\frac{\omega {v}^{\ast }}{4{a}^{2}\stackrel{˜}{p}\stackrel{˜}{h}}\right)\right]}^{1/2}$ . (21)

The physical interpretation of (18) is provided below. Consider now the situation when z is large enough, which corresponds to a high-frequency vibrations and low viscosity. In this case the asymptotic formulas for the Bessel functions have the form

${J}_{0}\left(z\right)\cong {\left(\frac{2}{\text{π}z}\right)}^{1/2}\mathrm{cos}\left(z-\frac{\text{π}}{4}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}{J}_{1}\left(z\right)\cong {\left(\frac{2}{\text{π}z}\right)}^{1/2}\mathrm{sin}\left(z-\frac{\text{π}}{4}\right)$

On the basis of (20) and (21) it is easy to show that for sufficiently large positive imaginary part z: ${J}_{0}\left(z\right)/{J}_{1}\left(z\right)\cong -i$ . Substituting (1) and further assuming

smallness ${\nu }^{\ast }$ in comparison with the $\frac{\omega }{{k}^{2}}$ , to obtain an approximate dispersion

equation, which is also contained in the [17]

$-{k}^{2}+\frac{{\omega }^{2}}{{a}^{3}}\left(1+\sqrt{\frac{{v}^{\ast }}{\omega }}\frac{\stackrel{˜}{р}}{\stackrel{˜}{h}R}\frac{l+i}{1.41}\right)=0$ (22)

where, in the pursuit of the viscosity ${v}^{\ast }$ to zero (and also tends ω to infinity),

we have a trivial result $\frac{\omega }{k}\to 0$ , which was obtained at low ω from Equation

(20). Equation (22) when an unacceptably high viscosities. In this case, the phase velocity C unlimited increases with ω. This example shows inconsistencies of various asymptotic estimates in the mid-frequency vibrations. Thus, the analysis of wave processes asymptotic methods in the first approximation is not possible to establish the limits of applicability of formulas and calculations to estimate the error. In this paper for solving spectral problems using a direct numerical integration of permitting relations of the type (12) by the method of orthogonal shooting in complex arithmetic. This approach avoids the above difficulties associated with the calculation of Bessel functions of complex argument. Another advantage is due to the specificity of the orthogonal sweep method, which is due to the procedure orthonormality can solve highly rigid system with a boundary layer. As a result of a numerical study has found that the problem of natural oscillations (12) admits no more than one complex value ω, corresponding vibrations of the shell together with the adjacent liquid layers to it. The rest found the Eigen values appeared purely imaginary. They correspond to a periodic motion of a fluid with almost stationary shell. Proper form corresponding complex values also are complex, that is, the phase of joint oscillations of the shell and liquid is not the same along the radius. In the case of steady-state oscillations all the calculated Eigen values k and their own forms are complex.

5. Numerical Results

Consider the case of natural oscillations, when the shell is filled with liquid. In Figures 1-9 and Figure 1, Figure 2, Figure 4, Figure 5 show, respectively, depending on the dispersion curves Reω, Imω, δ the wave number, k―the first mode, in which the damping coefficients of the smallest, and the Eigen values are complex Bat. In accordance with the numbering of graphs asked four different values of the coefficient η 1) 0.0009: 2) 0.0018 3) 0.15 4) 0.018 (а = 0.6199; $\stackrel{˜}{\rho }=0.0529$ ; $\stackrel{˜}{h}=0.0101$ ; $R=1$ ; ${\nu }_{0}=0.25$ .) for the remaining parameters according to (1) In Figure 3, Figure 6, Figures 7-9 to show their own forms Rev for values k equal to 1 and 8, respectively. It is easy to notice the difference in the behavior characteristic of the dispersion curves 1.2 and 3.4. In the last two cases, there is a wave number since a variable with only takes purely imaginary values corresponding to a periodic motion of the system. For curves 1.2 with less viscosity real part of the Eigen values Reω nonzero at all wave numbers and the damping rate has a finite limit at infinity. The greater the viscosity, the earlier start a periodic traffic (curves 3, 4) and the higher limit of the damping rate

Figure 1. Dependence of the real part of the complex frequencies (Reω) to wave numbers (k) for different values of η. 1-0.0009; 2-0.0018; 3-0.18, 4-0.19, 5-according to the formula (20); By the formula (22).

Figure 2. Dependence of the logarithmic decrement (δc) on the wave numbers (k) for different values of η. 1-0.0009; 2-0.0018; 3-0.18, 4-0.19, 5-0.20; 6-0.22.

Figure 3. Dependence of ϑ on the wave number r, for different values of the viscosity of a liquid. 1-0.0009; 2-0.0018; 3-0.18; 4-0.19.

(curves 1, 2). It follows that where is a minimum critical viscosity ηk, above which a zone of high wave numbers of the first mode, there are a periodic wave number. As a result of numerical experiment, it was found that the critical values of the coefficient of viscosity ηk, is in the range $\left[0.0120,0.0125\right]$ . Analyzing the dependence of energy dissipation on the wave number, two opposite tendencies should be noted. As the wave number increases, at a fixed amplitude, tangential

Figure 4. Dependence of the imaginary part of the complex frequencies (Im) to wave numbers (k) for different values of η:1-0.0009; 2-0.0018; 3-0.18; 4-0.19; by the formula (22).

Figure 5. Dependence of the spatial decrement on the wave number k for different values of η: 1-0.0009; 2-0.0018; 3-0.18; 4-0.19.

stresses linearly increase according to (6): c another, as shown in Figure 3, localization of the fluid motion amplitudes near the shell simultaneously results, which leads to a decrease in the mass of fluid involved in the motion, as well as tangential stresses. The difference in the behavior of curves 1, 2 and 3, 4 is due to which of the two tendencies prevails. At small wave numbers, a linear dependence of the Eigen function v on the radius is observed, that is, the entire mass of the liquid is involved in the motion. Ask increases, the central part of the liquid begins to “not keep up” with the vibrations of the shell, which leads to the localization of the amplitudes. The rate of localization depends on the viscosity of the liquid. If the localization occurs slowly, then starting from some k (owing to the growth of stresses), the self-motions become aperiodic (curves 3, 4). If, on the other hand, the average amplitude of the fluid oscillation decreases rapidly

Figure 6. Dependence of ϑ on the wave number r. When $\omega =8,\eta =0.0018$ .

Figure 7. Dependence of ϑ on the wave number. When $\omega =8,\eta =0.018$ .

Figure 8. Dependence of ϑ on the wave number r. When $\omega =16,\eta =0.0018$ .

Figure 9. Dependence of ϑ on the wave number r. When $\omega =16,\eta =0.018$ .

enough, the motions will always remain oscillatory (curves 1, 2). In this case, large voltage wave numbers prevail over voltages, and increase with increasing localization. In view of the latter circumstance, the damping coefficient always increases with increasing k. The linear dependence of the shape on the radius at small k also indicates the fulfillment of the flat-section hypothesis on which the elementary theory of viscoelastic rods is based. Using the Ritz method one can find the parameters of the Fought core model and determine the limits of applicability of this model in the framework of the hydrodynamic theory, but for a narrower class of straight rods of circular cross section. Variation equation of the principle of possible displacements, equivalent to the relations

$\begin{array}{l}\underset{v}{\int }h\left(\frac{\partial {u}_{\phi }}{\partial z}\delta \frac{\partial {u}_{\phi }}{\partial z}+{\rho }_{1}\frac{\partial {u}_{\phi }}{\partial z}\delta {u}_{\phi }\right){R}_{1}\text{d}\phi \text{d}z\\ \text{ }-\underset{v}{\int }\left({\sigma }_{r\phi }\delta {\epsilon }_{r\phi }+{\sigma }_{z\phi }\delta {\epsilon }_{z\phi }+{\rho }_{0}\frac{{\partial }^{2}{u}_{\phi }}{\partial {z}^{2}}\delta {u}_{\phi }\right)r\text{d}\phi \text{d}r\text{d}z=0\end{array}$ (23)

has the form. Choosing a linear function as the basis

${u}_{\phi }\left(r,z,t\right)=\phi \left(z,t\right)r,$ (24)

and after substituting (24) into (23) and the standard procedure, we obtain where the parameters β and a0 are expressed in terms of the polar moments of inertia of the shell I1 and liquid I0 as follows

$\beta =\frac{\eta {I}_{0}}{G{I}_{1}};\text{\hspace{0.17em}}{a}_{0}=a/{\left(1+\frac{{I}_{0}}{\stackrel{˜}{\rho }{I}_{1}}\right)}^{\frac{1}{2}}$

Equation (23) describes the torsional vibrations of a viscoelastic Feucht rod according to the relations

$\left(1+\beta \frac{\partial }{\partial t}\right)\frac{{\partial }^{2}\phi }{\partial {z}^{2}}=\frac{1}{{a}_{0}^{2}}\frac{{\partial }^{2}\phi }{\partial {t}^{2}},$ (25)

The solution of (25) is represented in the form

$\phi \left(z,t\right)={\phi }_{0}\mathrm{exp}\left(i\left(kz-\omega t\right)\right)$ .

where the following relations satisfy

${а}_{0}^{2}{к}^{2}\left(1-i\omega \beta \right)-{\omega }^{2}=0.$ (26)

Taking into account the relation I1/I0 = 4h, it is easy to see that Equation (26) coincides with Equation (22), which was obtained for the asymptotic solution of problem (16) for small oscillation frequencies. In Figure 1, Figure 4, the dotted lines show the dispersion curves of natural oscillations found from Equation (22). As follows from the figures, a satisfactory coincidence of dotted and continuous lines is observed in the region of small wave numbers whose upper bound exceeds unity in this case and increases with increasing viscosity of the liquid. In the short-wavelength range, there is a discrepancy due to the localization of the oscillation amplitudes near the shell. Small wave numbers correspond to the natural vibrations of long finite tubes. We now turn to an analysis of the steady-state oscillations of a shell filled with a liquid. Figures 1-9 shows the dispersion curves and waveforms for two values of the viscosity coefficient (below and above the critical value) 1) 0.0018, 2) 0.018 and the same values of the remaining parameters as in (22). In the first case of relatively low viscosity, the results of the calculation are in good agreement with the asymptotic solutions of the Goose Equation (18) at high frequencies.

6. Longitudinal-Transverse Vibrations

This section analyzes the stationary longitudinal-transverse vibrations of a shell filled with fluid, which according to (6) can be described by a system of four ordinary differential equations

$\begin{array}{l}\frac{\text{d}{\vartheta }_{r}}{\text{d}r}=-\frac{{\vartheta }_{r}}{r}-ik{\vartheta }_{z}-p\\ \frac{\text{d}{\vartheta }_{z}}{\text{d}r}=ik{\vartheta }_{r}+\frac{1}{\eta \omega }{\tau }_{ў}\\ \frac{\text{d}{\sigma }_{r}}{\text{d}r}=-{\rho }_{0}{\omega }^{2}{\vartheta }_{z}+2i\eta \omega \left(\frac{\text{d}{\vartheta }_{r}}{\text{d}r}-\frac{{\vartheta }_{r}}{r}\right)=ik{\tau }_{z}\\ \frac{\text{d}{\tau }_{z}}{\text{d}r}=-{\rho }_{0}{\omega }^{2}{\vartheta }_{z}+2\eta \omega k\left(\frac{\text{d}{\vartheta }_{r}}{\text{d}r}-ik{\vartheta }_{z}\right)-ik{\sigma }_{r}-\frac{{\tau }_{z}}{r}\end{array}$ (27)

With the boundary conditions

$\begin{array}{l}r=0:{\vartheta }_{r}=0,\text{\hspace{0.17em}}{\tau }_{z}=0;\\ r=R:D{\nabla }^{4}u+\frac{C}{R}\left(\frac{u}{R}+i{v}_{0}kw\right)+{\sigma }_{r}-{\rho }_{1}h{\omega }^{2}u=0;\\ C\left(i{v}_{0}k\frac{u}{R}-{\nabla }^{2}u\right)-{\tau }_{z}+{\rho }_{1}h{\omega }^{2}w=0;\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}C=\frac{E{h}_{0}}{1-{\nu }_{0}^{2}}.\end{array}$ (28)

The value of р in the first equation of system (27) is defined through the main unknowns according to the expression

$p=\frac{-{\sigma }_{r}+2i\eta \omega \left(iku+\frac{{\vartheta }_{r}}{r}\right)}{{\rho }_{0}{C}_{0}^{2}-i\omega \left(k+2\eta \right)}$ (29)

The spectral problems (27), (28), as in the case of longitudinal-transverse vibrations were solved by orthogonal shooting. To find the roots of the characteristic equation method were used Mueller.

7. Numerical Results

The results of numerical study of natural oscillations. Figure 10 shows the dispersion curves Reω the wave number k for the case of an incompressible (С0 = ∞ − dot-dash line) and compressible (С0 = 0.1 − solid line) of the liquid. Shell parameters and coefficients of viscosity taken following: h0 = 0.05; р = 1.8; v0 = 0.25; h = 6.011 × 10−4; к = −2η/3. Here and henceforth given dimensionless quantities

for which the units of length and mass density are $R,R{\left(\frac{{\rho }_{0}}{E}\right)}^{\frac{1}{2}},\frac{1}{{\rho }_{0}}$ . For an

incompressible fluid, there are two modes, corresponding mainly longitudinal (curve 1) and preferably a cross (curve 2) fluctuations in the shell, with complex Eigen values. All other traffic have their own imaginary Eigen values, that is a periodic in time. The dashed lines in Figure 11 are designated the dispersion curves corresponding to the vibrations of a shell with an ideal incompressible fluid. The solution of the latter problem is given below. It should be noted that, unlike the dry shell joint oscillations transverse vibrations of said sheath fluid density p1, it takes place on a smaller compared to the frequency of longitudinal vibrations in the entire range of the wave number. When administered viscosity oscillation frequency of the first mode decreases, apparently due to the involvement of additional masses in movement of fluid in the boundary layer and in the second mode appears critical wave number restricting oscillatory motions bottom region. In [25] , who investigated the steady oscillations, noted the desire for zero phase velocity of the lowest mode with decreasing frequency. Proper motion of the shell and the viscous compressible fluid has an infinite number of modes. The paper S. Vasin et al. [26] using asymptotic methods of solving, the latter effect could not be found. Figure 11 shows the dispersion curves for the first four events with a minimum of vibration frequencies (curves 3, 4, 5, 6) in ascending order of magnitude Reω. Comparing curves 1.2 and 3.4 together, we can see that the second worse than the first few vibration modes of the shell-compressible fluid to the selected parameters are satisfactorily described by a model of an incompressible fluid in the region of wave numbers k < 1. This gives grounds for the study of the said system in the first approximation neglect the compressibility of the fluid. System elastic shell is a viscous liquid dissipations-inhomogeneous viscoelastic body at a radial coordinate. Moreover, in contrast to the earlier torsional vibrations here for an incompressible fluid, there are

Figure 10. Addiction Reω with the wave number k in the case of an incompressible fluid.

Figure 11. Addiction Im a the wave number in the case of a compressible fluid.

two, and compressible-unlimited number of vibration modes. It is interesting to find out how this system can be shown a synergistic effect. Figure 11 shows the dispersion curves (2) for the following parameters of the shell and liquid:

$h=0.05;\text{\hspace{0.17em}}{\rho }_{3}=80;\text{\hspace{0.17em}}v=0.25t;\text{\hspace{0.17em}}\eta =7.071×{10}^{-4};\text{\hspace{0.17em}}{C}_{0}=\infty$

Dash-dotted lines correspond to fluctuations in the dry shell. The dashed lines show the frequency dependence for the case of an ideal fluid v = 0. In contrast with the previously discussed embodiment, the density p = 8, in this case partial frequency (v = 0) of the longitudinal and transverse vibrations of the shell with a perfect fluid intersect. It is natural to expect that the v near the intersection of partial frequencies will be a strong connectedness of both modes, leading to increased energy, resulting in a synergistic effect. Indeed, the presence of events demonstrates the effect of the conversion of Vina-longitudinal mode in transverse and longitudinal cross-section in a change of the wave number in the vicinity of the intersection of partial frequencies. Violation of the monotony of growth and synergies. Compared to the previous description of this effect there are two features. Firstly, the effect is far from the place of approximation curves of two modes, secondly, damping factor curves do not intersect. In [27] investigated the coherence of joint oscillations of ideal compressible gas and the shell with the help of diagrams wines. As he examined the frequency of partial oscillations of gas in rigid walls and an empty shell. Returning to Figure 12, Figure 13 we note a similar manifestation of the effect of wines in places of convergence curves 4.5 and 5.6. In these areas in Figure 14 there is a synergistic effect for the curves. It is interesting to trace the influence of fluid viscosity on connectivity modes. 3.4 Curves in Figure 15 correspond to the value of the viscosity coefficient η = 0.11 at constant other parameters. In this case, fashion predominantly transverse vibrations are defined on a finite interval of the wave of change, and the effect of guilt is not observed, indicating a loose coupling modes. Another large increase in viscosity (η = 0.13, curve 5) leads to the fact that fashion is everywhere transverse vibrations becomes a periodic and у longitudinal oscillations appear critical wave numbers, limiting the scope of the vibration motions of the top. The physical nature of the observed effect is revealed when analyzing the vibrations of a shell filled with a perfect fluid. The equations of harmonic oscillations of an ideal liquid is easy to deduce from (27), formally putting viscosity coefficients equal to zero.

$\frac{\text{d}{\vartheta }_{r}}{\text{d}r}=-\frac{{\vartheta }_{r}}{r}-ik{\vartheta }_{z}-\frac{\sigma }{{\rho }_{0}{C}_{0}^{2}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\frac{\text{d}\sigma }{\text{d}r}=-{\rho }_{0}{\omega }^{2}{\vartheta }_{r},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\sigma =i{\rho }_{0}{\omega }^{2}{\vartheta }_{z}$ (30)

General solution of (26) satisfying the finiteness condition unknown at zero, has the form

$\begin{array}{l}{\vartheta }_{z}=A{J}_{0}\left(qr\right);\text{\hspace{0.17em}}\text{\hspace{0.17em}}\sigma =i{\rho }_{0}\frac{{\omega }^{2}}{k}A{I}_{0}\left(qr\right)\\ {\vartheta }_{z}=i\frac{q}{R}A{I}_{1}\left(qr\right);\text{\hspace{0.17em}}\text{\hspace{0.17em}}{q}^{2}=\frac{{\omega }^{2}}{{C}_{0}^{2}}-{k}_{0}^{2}\end{array}$ (31)

Figure 12. The dependence of ReW on the wave number r.

Figure 13. The dependence of ReW on the wave number r.

where А arbitrary constant: J0, Jl,-Bessel functions of zero and first order, respectively. The boundary conditions at the r = R similarly written conditions (28)

$\begin{array}{l}D{\nabla }^{4}u+\frac{C}{R}\left(\frac{u}{R}+i{v}_{0}kw\right)+{\sigma }_{r}-{\rho }_{1}h{\omega }^{2}u=0;\\ C\left(i{v}_{0}k\frac{u}{R}-{\nabla }^{2}u\right)+{\rho }_{1}h{\omega }^{2}w=0;\end{array}$ (32)

where w―axial movement of the shell, which is not now coincides with the axial movement of the liquid. After substitution of the solutions (22) of (23) there is a system of homogeneous linear algebraic equations in the unknown А and U1. The roots of the determinant of this system are the desired Eigen values, and its decision to define the relation between А and U1.

For an incompressible fluid, there are two real own Bessel functions I0 and I1

${\omega }_{1}=R{\left(\frac{E}{{\rho }_{1}}\right)}^{\frac{1}{2}};\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\omega }^{2}={\left[\frac{Е}{{R}_{1}{\rho }_{1}}\left(1+{h}^{2}{R}_{1}{к}^{4}\right)/\left(1+\frac{{\rho }_{n}I\left(k{R}_{1}\right)}{h{\rho }_{1}{I}_{1}\left(k{R}_{1}\right)к}\right)\right]}^{\frac{1}{2}}$ (33)

Figure 14. The dependence Reω on the wave number k (for the case of an incompressible fluid).

Figure 15. The dependence of Imω on the wave number k.

Unlike dry shell here second frequency locking is absent and the phase speed at low k equal to the

${C}_{R}={\left(\frac{Eh}{2{\rho }_{0}{R}_{1}}\right)}^{\frac{1}{2}}$ (34)

which coincides with the speed of the wave Ressala (see. the review at the beginning of this chapter). In the case of a compressible fluid $\nu =0$ and limiting the phase velocity of the transverse mode oscillation in the shell $k\to 0$ is the velocity of waves Cortège Zhukov sky.

${C}_{k}=\frac{{C}_{0}{C}_{R}}{{\left({C}_{0}^{2}+{C}_{R}^{2}\right)}^{\frac{1}{2}}}$ (35)

Numerical study showed that the critical value Ck does not depend on the viscosity of the liquid, but with increasing η weakening the dependence of oscillations of Poisson’s ratio, so that the ratio $\left(\mathrm{max}im\varpi \right)/\left(\mathrm{min}im\varpi \right)\to 1$ and own form U it becomes flat. As follows from the above results, generally within the engineering problem statement, we cannot adequately describe the longitudinal vibrations of the cylindrical shell filled with a viscous fluid via rod theory.

8. Conclusions

1) Numerical investigation has shown that the critical value Vk does not depend on the viscosity of the liquid, but as the value of η increases, the dependence of the Poisson’s coefficient is weakened, so that the ratio $\left(\mathrm{max}im\omega \right)/\left(\mathrm{min}im\omega \right)\to 1$ and the proper form U becomes flat. However, in some particular cases, namely: for a large viscosity or at a critical value of the Poisson ratio, it is possible to propose a method for estimating the natural frequency based on a bar model of the type (26).

2) Analyzing the dependence of energy dissipation on the number of waves, it should be noted that there are two opposite tendencies, since the number of waves and tangential stresses р increases for a fixed amplitude v. And also, as shown by numerical results, localization of fluid motion amplitudes near the shell occurs simultaneously, which results in a decrease in the mass of fluid involved in the motion, as well as tangential stresses р.

3) For low viscosity, the frequency Rek of both modes is close to one in the low-frequency region, and at high frequencies the phase velocity Cy corresponds to the tendency to speed in the dry shell. Damping coefficients grow approximately linearly, and in the second mode, this coefficient is always greater than that of the first. In the case of a higher viscosity over the entire range of variation, the frequency ω of the real part of the wave numbers Rek2,1 ≥ Rek2,2, and the phase velocity Cy with increasing ω tends to infinity.

Conflicts of Interest

The authors declare no conflicts of interest.

Cite this paper

Ibrahimovich Safarov, I. , Khudoyberdiyevich Teshaev, M. and Ixtiyorovich Boltayev, Z. (2018) Properties of Wave Motion in a Cylindrical Shell, Interacting with Viscous Liquid. Open Access Library Journal, 5, 1-23. doi: 10.4236/oalib.1104563.

 [1] Padilla, F., de Billy, M. and Quentin, G. (1999) Theoretical and Experimental Studies of Surface Waves on Solid-Fluid Interfaces When the Value of the Fluid Sound Velocity Is Located between the Shear and the Longitudinal Ones in the Solid. The Journal of the Acoustical Society of America, 106, 666-673. https://doi.org/10.1121/1.427084 [2] Glorieux, C., Van de Rostyne, K., Nelson, K., Gao, W., Lauriks, W. and Thoen, J. (2001) On the Character of Acoustic Waves at the Interface between Hard and Soft Solids and Liquids. The Journal of the Acoustical Society of America, 110, 1299-1306. https://doi.org/10.1121/1.1396333 [3] Zhu, J., Popovics, J.S. and Schubert, F. (2004) Leaky Rayleigh and Scholte Waves at the Fluid-Solid Interface Subjected to Transient Point Loading. The Journal of the Acoustical Society of America, 116, 2101-2110. https://doi.org/10.1121/1.1791718 [4] Bayon, A., Gascon, F. and Nieves, F.J. (2005) Estimation of Dynamic Elastic Constans from the Amplitude and Velocity of Rayleigh Waves. The Journal of the Acoustical Society of America, 117, 3469-3477. https://doi.org/10.1121/1.1898663 [5] Safarov, I.I., Teshaev, M.Kh. and Akhmedov, M.Sh. (2018) Free Oscillations of a Toroidal Viscoelastic Shell with a Flowing Liquid. American Journal of Mechanics and Applications, 6, 37-49. http://www.sciencepublishinggroup.com/j/ajma [6] Safarov, I.I. and Boltayev, Z.I. (2018) Methods for Assessing the Seismic Resistance of Subterranean Hydro Structures under the Influence of Seismic Waves. American Journal of Physics and Applications, 6, 51-62. http://www.sciencepublishinggroup.com/j/ajpa [7] Safarov, I.I., Teshayev, M.Kh., Boltayev, Z.I. and Akhmedov, M.Sh. (2017) Damping Properties of Vibrations of Three-Layer VIscoelastic Plate. International Journal of Theoretical and Applied Mathematics, 3, 191-198. [8] Gazis, D.C. (1959) Three-Dimension Investigation of the Propagation of Waves in Hollow Circular Cylinders. II. Numerical Results. The Journal of the Acoustical Society of America, 31, 573-578. https://doi.org/10.1121/1.1907754 [9] Rosenberg, R.L. and Thurston, R.N. (1977) Relationship between Plate and Surface Modes of a Tube. The Journal of the Acoustical Society of America, 61, 1499. https://doi.org/10.1121/1.381450 [10] Ter-Hakobyan, G.L. (2013) Refinement of the Results of the Effect of the Fluid on the Propagation of Waves in an Elastic Cylindrical Shell. Journal Basic Research, Engineering Sciences, 10, 516-520. [11] Sorokin, S.V. (1997) Fluid-Structure Interaction and Structural Acoustics. Book of Lecture Notes, Technical University of Denmark, Lyngby, 188 p. [12] Safarov, I.I., Akhmedov, M.Sh. and Boltayev, Z.I. (2016) The Actual Waves in Layered Media. Lambert Academic Publishing, Saarbrücken, 192s. [13] Safarov, I.I., Boltayev, Z.I. and Akhmedov, M.Sh. (2016) Properties of Wave Motion in a Fluid-Filled Cylindrical Shell. Lambert Academic Publishing, Saarbrücken, 105 р. [14] Guz A.N. ,et al. (2016)Wave Propagation in a Cylindrical Shell with a Viscous Compressible Fluid Prikladnaia Mekhanika 16, 10-20. [15] Safarov, I.I., Akhmedov, M.Sh. and Boltayev, Z.I. (2015) Dissemination Sinusoidal Waves in of a Viscoelastic Strip. Global Journal of Science Frontier Research: F Mathematics & Decision Sciences, 15, 39-60. [16] Safarov, I.I., Teshaev, M.X., Akhmedov, M.Sh. and Rajabov, O. (2017) Distribution Natural Waves on the Viscoelastic Cylindrical Body in Plane Strain State. Case Studies Journal, 6, 1-8. http://www.casestudiesjournal.com [17] Safarov, I.I., Teshaev, M.X., Akhmedov, M.Sh. and Boltaev, Z.I. (2017) Properties of Wave Motion in a Cylindrical Shell Is Contact with a Viscous Fluid. Case Studies Journal, 6, 9-35. http://www.casestudiesjournal.com [18] Safarov, I.I., Teshaev, M.Kh., Boltaev, Z.I. and Nuriddinov, B.Z. (2017) Of Own and Forced Vibrations of Dissipative Inhomogeneous Mechanical Systems. Applied Mathematics, 8, 1001-1015. http://www.scirp.org/journal/am https://doi.org/10.4236/am.2017.87078 [19] Safarov, I.I., Teshaev, M.X., Akhmedov, M.Sh. and Boltaev, Z.I. (2017) Distribution Free Waves in Viscoelastic Wedge with and Arbitrary Angle Tops. Applied Mathematics, 8, 736-745. https://doi.org/10.4236/am.2017.85058 [20] Safarov, I.I., Akhmedov, M.Sh. and Boltayev, Z.I. (2016) Ducting in Extended Plates of Variable Thickness. Global Journal of Science Frontier Research: F Mathematics & Decision Sciences, 16, 33-66. [21] Safarov, I.I., Teshaev, M.Kh. and Boltayev, Z.I. (2016) Propagation of Linear Waves in Extended Lamellar Bodies. Lambert Academic Publishing, Saarbrücken, 315 p. [22] Safarov, I.I., Teshaev, M.H. and Boltaev, Z.I. (2012) Wave Processes in Mechanical Waveguide. Lambert Academic Publishing, Saarbrücken, 217 p. [23] Bozorov, M.B., Safarov, I.I. and Shokin, Y.I. (1966) Numerical Simulation of Vibrations Dissipative Homogeneous and Heterogeneous Mechanical Systems. Moscow, 188 p. [24] Kayumov, S.S. and Safarov, I.I. (2002) Propagation and Diffraction of Waves in Dissipative—Inhomogeneous Cylindrical Deformable Mechanical Systems. Tashkent, 214 p. [25] Grinchenko, V.T. and Myaleshka, V.V. (1981) Harmonic Waves in Elastic Bodies. Kiev, 284. [26] Vasin, S.V. and Mikolyuk, V.V. (1983) Free Oscillations Tolerable Cylindrical Shells Separated by a Viscous Fluid Hydro Aeromechanics. And the Theory of Elasticity. No. 3, 108-116. [27] Beginners, Y.N. (1996) The Study of the Spectra of Natural Frequencies of Cylindrical Shells Containing a Compressible Fluid. Science, 74 p.