^{1}

^{2}

^{2}

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.

Quantitative data and some qualitative analysis of the dispersion properties of normal waves in a cylinder with a liquid were carried out in [

An infinite length of deformable (viscoelastic) cylindrical shell of radius R with constant thickness h 0 , density ρ 0 , Poisson’s ratio ν 0 , filled with a viscous fluid with density at equilibrium. Fluctuations of a shell under a load, the density of which is denoted p_{1}, p_{2}, p_{n} respectively, can be described by following [

L u − L 0 ∫ 0 t R L ( t − τ ) L u ( τ ) d τ = ( 1 − ν 0 2 ) E 0 h 0 р + ρ 0 ( 1 − ν 0 2 ) E 0 ( ∂ 2 u ∂ t 2 ) , (1)

Here u = u ( u r , u θ , u z ) ―displacement vector points of the middle surface of the shell and membranes for Kirchhoff-Love it has a dimension equal to three ( u r = u ; u θ = v ; u z = w ) , 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 [_{n} and the sign “−” before the last component of the inertial member says that is considered positive motion towards the center of curvature); R L ( t − τ ) ―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 [

∂ ϑ ∂ t − ν ∗ Δ ϑ + 1 ρ 0 ∗ g r a d P − ν ∗ 3 g r a d d i v ϑ = 0

1 ρ 0 ∗ ∂ ρ ∗ ∂ t + d i v ϑ = 0 ; ∂ P ∂ ρ ∗ = a 0 2 , a 0 = c o n s t . u ˙ z = ϑ z , u ˙ r = ϑ r , u ˙ θ = ϑ θ , q z = − p r z , q r = − p r , q θ = − p r θ . (2)

p r z = μ ∗ ( ∂ ϑ z ∂ r + ∂ ϑ r ∂ z ) ; p r r = − p + λ ∗ ( ∂ ϑ r ∂ r + ∂ ϑ z ∂ z + ϑ r r ) + 2 μ ∗ ∂ ϑ r ∂ r ; p r θ = μ ∗ ( 1 r ∂ ϑ z ∂ θ + ∂ ϑ θ ∂ r − ϑ θ r ) .

Here, in the Equations (2) ϑ = ϑ ( ϑ r , ϑ θ , ϑ z ) ―the velocity vector of fluid particles; ρ ∗ and Р―disturbance density and fluid pressure; ρ 0 ∗ and а_{0}―density and sound velocity in the fluid at rest; ν ∗ , μ ∗ ―kinematic and dynamic viscosity;

for the second viscosity coefficient λ ∗ accepted ratio λ = − 2 3 μ • ; р r z , р r r , р r θ ―

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.

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

u ( r , t ) = u 1 ( r , t ) e − i ω R t ,

where ψ ( r , t ) ―slowly varying function of time, ω R ―real constant. Next, using the procedure of freezing [

L [ 1 − Γ C ( ω R ) − i Γ S ( ω R ) ] u 1 = ( 1 − ν 0 2 ) E 0 h 0 р + ρ 0 ( 1 − ν 0 2 ) E 0 ( ∂ 2 u 1 ∂ t 2 ) , (3)

where, for shell Kirchhoff-Love

L = ( ∂ 2 ∂ x 2 + 1 − ν 0 2 R 2 ∂ 2 ∂ φ 2 1 + ν 0 2 R ∂ 2 ∂ x ∂ φ ν 0 R ∂ ∂ x 1 + ν 0 2 R ∂ 2 ∂ x ∂ φ 1 + ν 0 2 ( 1 + 4 a ) ∂ 2 ∂ x 2 + ( 1 + a ) ∂ 2 ∂ φ 2 1 R 2 ∂ ∂ φ − a ( 2 − ν 0 ) ∂ 3 ∂ x 2 ∂ φ − a R 2 ∂ 3 ∂ φ 3 ν 0 R ∂ ∂ x 1 R 2 ∂ ∂ φ − a ( 2 − ν 0 ) ∂ 3 ∂ x 2 ∂ φ − a R 2 ∂ 3 ∂ φ 3 1 R 2 + a ( ∂ 2 ∂ x 2 + 1 R 2 ∂ 2 ∂ φ 2 ) 2 ) ,

Γ C ( ω R ) = ∫ 0 ∞ R L ( τ ) cos ω R τ d τ , Γ S ( ω R ) = ∫ 0 ∞ R L ( τ ) sin ω R τ d τ respectively, cosine

and sine Fourier transforms relaxation kernel material. As an example, the viscoelastic material take three parametric kernel relaxation R L ( t ) = A e − β t / t 1 − α , ρ―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 ξ = x / R and multiply by R^{2} system (3). The matrix of the resulting system will take the form

L = ( ∂ 2 ∂ ξ 2 + 1 − ν 0 2 ∂ 2 ∂ φ 2 1 + ν 0 2 ∂ 2 ∂ ξ ∂ φ ν 0 ∂ ∂ ξ 1 + ν 0 2 ∂ 2 ∂ ξ ∂ φ 1 − ν 0 2 ( 1 + 4 a ) ∂ 2 ∂ ξ 2 + ( 1 + a ) ∂ 2 ∂ φ 2 ∂ ∂ φ − a ( 2 − ν 0 ) ∂ 3 ∂ ξ 2 ∂ φ − a ∂ 3 ∂ φ 3 ν 0 ∂ ∂ ξ ∂ ∂ φ − a ( 2 − ν 0 ) ∂ 3 ∂ ξ 2 ∂ φ − a R 2 ∂ 3 ∂ φ 3 1 R 2 + a ( ∂ 2 ∂ ξ 2 + 1 R 2 ∂ 2 ∂ φ 2 ) 2 ) . (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:

∂ p r θ ∂ r + 2 p r θ r + ∂ p θ z ∂ z = ρ 0 ∗ ϑ ¨ θ , p r θ = μ ∗ ( ∂ ϑ θ ∂ r − ϑ θ r ) , p θ z = μ ∗ ∂ ϑ θ ∂ r , r = R 1 : G ¯ h 0 ∂ 2 u θ ∂ z 2 − ( ρ 0 h u ¨ θ ± σ φ r ) = 0 , G ¯ = E ¯ 2 ( 1 + v 0 ) , r = 0 : p r θ = 0 , E ¯ = Е 0 ( 1 − Γ C ( ω R ) − i Γ S ( ω R ) ) (5)

- Longitudinal transverse vibrations:

∂ р r r ∂ r + p r r − p θ θ r + ∂ p r z ∂ z = ρ 0 ∗ ∂ ϑ r ∂ t ∂ p r z ∂ r + p r z r + ∂ p z z ∂ z = ρ 0 ∗ ∂ ϑ z ∂ t

p r r = − p + λ ∗ k η d i v ϑ + 2 μ ∗ ∂ ϑ r ∂ r , p θ θ = − p + λ ∗ d i v ϑ + 2 μ ∗ ϑ r r p z z = − p + λ ∗ d i v ϑ + 2 μ ∗ ∂ ϑ z ∂ z (6)

р r z = μ ∗ ( ∂ ϑ z ∂ r + ∂ ϑ r ∂ z ) , р r θ = μ ∗ ( 1 r ∂ ϑ z ∂ θ + ∂ ϑ θ ∂ r + ϑ r r ) ∂ ρ ∗ ∂ t + ρ 0 d i v ϑ = 0 , d i v ϑ = ∂ ϑ r ∂ r + ϑ r r + ∂ ϑ z ∂ z , ∂ p ∂ ρ = а 0 2

r = R 1 : D ¯ ∂ 4 u r ∂ z 4 + C ¯ R 1 ( u r R 1 + ν 0 ∂ u z ∂ z ) + p r r + ρ 0 h 0 ∂ 2 u r ∂ t 2 = 0 , C ¯ ( ∂ 2 u r ∂ z 2 + v 0 R ∂ u r ∂ z ) − ( p r z ± ρ 0 h 0 ∂ 2 u z ∂ t 2 ) = 0 ,

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 [

( р r r , р r z , р r θ , u , ϑ ) Т = ∑ m ( σ r m ( ξ , φ , t ) , τ z m ( ξ , φ , t ) , τ φ m ( ξ , φ , t ) , u m ( ξ , φ , t ) , ϑ m ( ξ , φ , t ) ) Т , (7)

where u m ( ξ , φ , t ) = u m { U m , V m , W m } T , ϑ ( ξ , φ , t ) = ϑ m { ϑ r m , ϑ θ m , ϑ z m } T .

Expressions (7) in the form

( σ r m ( ξ , φ , t ) , τ z m ( ξ , φ , t ) , τ φ m ( ξ , φ , t ) ) Т = ( σ r cos ( m φ ) , τ z cos ( m φ ) , τ φ sin ( m φ ) ) T e i к z − i ω t ,

( u m ( ξ , ϕ , t ) , ϑ m ( ξ , φ , t ) ) Т = ( U m cos ( m φ ) , V m sin ( m φ ) , W m cos ( m φ ) , ϑ r cos ( m φ ) , ϑ θ cos ( m φ ) , ϑ z cos ( m φ ) ) T e i к z − i ω t ,

where σ r , τ z , τ φ , U m , V m , W m , ϑ r , ϑ θ , ϑ 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 , ⋯ . 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 θ = р r z = 0 , ϑ r = 0 . If the outer surface г = R assumed stationary, then u r = u z = u φ = 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 [

δ c = 2 π | Im ω | / Re ω (8)

Decrement is similar to the spatial

δ y = 2 π | Im k | / Re k .

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

c c = Re ω R , c y = ω 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 Y ¯ c . It can be in a superposition of waves

Y ¯ c = ( r , z , t ) = ∑ n ∫ − ∞ ∞ Y n c ( r , k ) exp [ t ( k z − ϖ n ( k ) t ) ] d k , (9)

where vectors Y ¯ n c are their own form of the problem of natural oscillations, normalized so that the spatial Fourier spectrum of the initial disturbance f ¯ ( r , z ) = Y ¯ c ( r , z , 0 ) forms a linear combination

f ¯ ( r , z ) = ∫ − ∞ ∞ F ( r , k ) e i k z d k , f ¯ ( r , k ) = ∑ n Y ¯ n c ( r , k ) . (10)

Similarly, the main vector of unknowns Y ¯ y boundary value problem is calculated according to the expression

Y ¯ y ( r , z , t ) = ∑ n ∫ − ϖ ϖ Y ¯ k y ( r , ω ) exp [ i k ( ω ) z − ω t ] d ω (11)

where Y ¯ k y forms steady-state oscillation, the linear combination of which should form a Fourier spectrum given boundary perturbation

q ¯ ( r , t ) = Y ¯ y ( r , 0 , t ) , q ¯ ( r , t ) = ∫ − ∞ ∞ q ( r , ω ) e − i ω t d ω , q ¯ ( r , ω ) = ∑ n ∫ − ∞ ∞ Y n y ( r , ω )

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 [

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

d τ φ d r = − ( ρ 0 ∗ ω 2 − i μ ∗ 2 ξ 2 ω ) ϑ θ − 2 τ φ r d v d r = ϑ θ r + i ω μ ∗ τ φ r = R 1 : h 0 ( G ¯ ξ 2 − ξ ρ 0 ω 2 ) ϑ θ ± τ φ = 0 r = 0 : τ φ = 0 (12)

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

d 2 ϑ θ d r 2 + d ϑ θ r d r + ( − ξ 2 + i ω v ∗ − 1 r 2 ) ϑ θ = 0 ; v ∗ = μ ∗ ρ 0 ∗ (13)

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

v = A 1 J 1 ( r − k 2 + i ω v ∗ ) = 0 . (14)

where J_{1}-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 ( R 1 − k 2 + i ω v ∗ ) = 0 (15)

from whence

ω n = − i ( v ∗ k 2 + Γ т 2 ) (16)

in the case of natural oscillations and

k n = − Γ n 2 + i ω v ∗ (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 ν ∗ . 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 + ω 2 a 2 + ω v ∗ a 3 p ˜ h ˜ R 2 + ( z J 0 ( z ) J 1 ( z ) − 2 ) = 0 (18)

which was first obtained in A. Guz [

p ˜ = ρ ∗ ρ 0 ; h ˜ = h R 1 ; z = R 1 − k 2 + i ω v ∗ ; a = G ρ 0

shear wave velocity shell: J_{0}-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 J_{0} and J_{1} power series

J 0 = 1 − z 2 4 = ⋯ ; J 1 ( z ) = z 2 ( 1 − z 2 8 + ⋯ ) ; (19)

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

− k 2 + ω 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 + ω 2 a 2 + i ω v ∗ 4 a 2 p ˜ h ˜ ( k 2 − i ω v ∗ ) = 0 (20)

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

k = ω a [ ( 1 + 1 4 p ˜ h ˜ ) / ( 1 − ω v ∗ 4 a 2 p ˜ h ˜ ) ] 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 ( z ) ≅ ( 2 π z ) 1 / 2 cos ( z − π 4 ) , J 1 ( z ) ≅ ( 2 π z ) 1 / 2 sin ( z − π 4 )

On the basis of (20) and (21) it is easy to show that for sufficiently large positive imaginary part z: J 0 ( z ) / J 1 ( z ) ≅ − i . Substituting (1) and further assuming

smallness ν ∗ in comparison with the ω k 2 , to obtain an approximate dispersion

equation, which is also contained in the [

− k 2 + ω 2 a 3 ( 1 + v ∗ ω р ˜ h ˜ R l + i 1.41 ) = 0 (22)

where, in the pursuit of the viscosity v ∗ to zero (and also tends ω to infinity),

we have a trivial result ω k → 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.

Consider the case of natural oscillations, when the shell is filled with liquid. In Figures 1-9 and

(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 [ 0.0120 , 0.0125 ] . 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

stresses linearly increase according to (6): c another, as shown in

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

∫ v h ( ∂ u φ ∂ z δ ∂ u φ ∂ z + ρ 1 ∂ u φ ∂ z δ u φ ) R 1 d φ d z − ∫ v ( σ r φ δ ε r φ + σ z φ δ ε z φ + ρ 0 ∂ 2 u φ ∂ z 2 δ u φ ) r d φ d r d z = 0 (23)

has the form. Choosing a linear function as the basis

u φ ( r , z , t ) = φ ( z , t ) r , (24)

and after substituting (24) into (23) and the standard procedure, we obtain where the parameters β and a_{0} are expressed in terms of the polar moments of inertia of the shell I_{1} and liquid I_{0} as follows

β = η I 0 G I 1 ; a 0 = a / ( 1 + I 0 ρ ˜ I 1 ) 1 2

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

( 1 + β ∂ ∂ t ) ∂ 2 φ ∂ z 2 = 1 a 0 2 ∂ 2 φ ∂ t 2 , (25)

The solution of (25) is represented in the form

φ ( z , t ) = φ 0 exp ( i ( k z − ω t ) ) .

where the following relations satisfy

а 0 2 к 2 ( 1 − i ω β ) − ω 2 = 0. (26)

Taking into account the relation I_{1}/I_{0} = 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

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

d ϑ r d r = − ϑ r r − i k ϑ z − p d ϑ z d r = i k ϑ r + 1 η ω τ ў d σ r d r = − ρ 0 ω 2 ϑ z + 2 i η ω ( d ϑ r d r − ϑ r r ) = i k τ z d τ z d r = − ρ 0 ω 2 ϑ z + 2 η ω k ( d ϑ r d r − i k ϑ z ) − i k σ r − τ z r (27)

With the boundary conditions

r = 0 : ϑ r = 0 , τ z = 0 ; r = R : D ∇ 4 u + C R ( u R + i v 0 k w ) + σ r − ρ 1 h ω 2 u = 0 ; C ( i v 0 k u R − ∇ 2 u ) − τ z + ρ 1 h ω 2 w = 0 ; C = E h 0 1 − ν 0 2 . (28)

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

p = − σ r + 2 i η ω ( i k u + ϑ r r ) ρ 0 C 0 2 − i ω ( k + 2 η ) (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.

The results of numerical study of natural oscillations. _{0} = ∞ − dot-dash line) and compressible (С_{0} = 0.1 − solid line) of the liquid. Shell parameters and coefficients of viscosity taken following: h_{0} = 0.05; р = 1.8; v_{0} = 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 ( ρ 0 E ) 1 2 , 1 ρ 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 _{1}, 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 [

two, and compressible-unlimited number of vibration modes. It is interesting to find out how this system can be shown a synergistic effect.

h = 0.05 ; ρ 3 = 80 ; v = 0.25 t ; η = 7.071 × 10 − 4 ; C 0 = ∞

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 [

d ϑ r d r = − ϑ r r − i k ϑ z − σ ρ 0 C 0 2 , d σ d r = − ρ 0 ω 2 ϑ r , σ = i ρ 0 ω 2 ϑ z (30)

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

ϑ z = A J 0 ( q r ) ; σ = i ρ 0 ω 2 k A I 0 ( q r ) ϑ z = i q R A I 1 ( q r ) ; q 2 = ω 2 C 0 2 − k 0 2 (31)

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

D ∇ 4 u + C R ( u R + i v 0 k w ) + σ r − ρ 1 h ω 2 u = 0 ; C ( i v 0 k u R − ∇ 2 u ) + ρ 1 h ω 2 w = 0 ; (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 U_{1}. The roots of the determinant of this system are the desired Eigen values, and its decision to define the relation between А and U_{1}.

For an incompressible fluid, there are two real own Bessel functions I_{0} and I_{1}

ω 1 = R ( E ρ 1 ) 1 2 ; ω 2 = [ Е R 1 ρ 1 ( 1 + h 2 R 1 к 4 ) / ( 1 + ρ n I ( k R 1 ) h ρ 1 I 1 ( k R 1 ) к ) ] 1 2 (33)

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

C R = ( E h 2 ρ 0 R 1 ) 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 ν = 0 and limiting the phase velocity of the transverse mode oscillation in the shell k → 0 is the velocity of waves Cortège Zhukov sky.

C k = C 0 C R ( C 0 2 + C R 2 ) 1 2 (35)

Numerical study showed that the critical value C_{k} 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 ( max i m ϖ ) / ( min i m ϖ ) → 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.

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 ( max i m ω ) / ( min i m ω ) → 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 р_{zφ} 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 р_{rφ}.

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 C_{y} 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 Rek_{2,1} ≥ Rek_{2,2}, and the phase velocity C_{y} with increasing ω tends to infinity.

Safarov, I.I., Teshaev, M.K. and Boltayev, Z.I. (2018) Properties of Wave Motion in a Cylindrical Shell, Interacting with Viscous Liquid. Open Access Library Journal, 5: e4563. https://doi.org/10.4236/oalib.1104563