Abstract
The article considers dynamic processes mathematical modeling in a mechanical system, consisting of an elastic hollow cylinder, surrounded by an elastic medium and containing viscous liquid and vibrating coaxial rigid cylinder. The amplitude frequency characteristic for investigating bending cylinder oscillations as onemass system is defined. It is shown, that the constructed amplitude characteristic makes it possible to define the considered system resonance frequencies oscillations. The calculations demonstrated the significance of taking into account viscous liquid inertia and the surrounding elastic medium.
1. Introduction
The problems of mathematical modeling of elastic construction elements interaction with liquid are of theoretical and practical interest [1]. For example, [28] deal with hydroelasticity problems of homogeneous plates, interacting with viscous liquid layer in various settings. Reference [9] considers the bending oscillations of a flat channel wall, interacting with pulsating viscous liquid. References [1012] investigate the hydroelastic oscillations of threelayered plates, interacting with viscous pulsating liquid, as well as under inertial excitation. References [1318] research in the axissymmetric problem of geometrically regular and ribbed cylinder shells, as well as the ones, forming an annular channel. The nonlinear longitudinal deformation wave propagation in cylinder shells, filled with viscous liquid and in axissymmetric setting is presented in [19, 20]. Reference [21] considers membrane oscillations on an elastic foundation on the pool bottom, the pool being filled with an ideal incompressible liquid with a free surface. References [22, 23] are devoted to the plate oscillation investigation, the plate standing on the elastic foundation and interacting with viscous liquid layer.
However, mathematical modeling of hydroelastic oscillations of a hollow cylinder, surrounded by an elastic medium and containing a viscous liquid and an oscillator, remained beyond the framework of the sources, mentioned above. Therefore, the problem under consideration is of undoubted interest for the research.
2. Mathematical formulation
Let us consider an elastic hollow cylinder, surrounded by an elastic medium (Fig. 1). The cylinder length is $l$, its thickness ${h}_{0}$ is significantly less than inner radius ${R}_{1}$. The elastic cylinder is simply supported at its edges. The rigid cylinder with radius $R$ is inside the elastic cylinder, mentioned above. The inner and outer cylinders surfaces are parallel to each other. The radial clearance between the cylinders is filled with viscous incompressible liquid. The seals, hindering the liquid leakage out of radial clearance, are supposed to exist on the edges. The rigid cylinder produces harmonic oscillations according to the set law in the plane, being normal to its axis. The hollow cylinder, surrounded by elastic medium, makes bending oscillations, caused by inner cylinder vibrations. Cylinders oscillations amplitudes are considerably less, than liquid layer thickness $\delta $. The elastic medium, surrounding the inner cylinder, can be considered as Winkler foundation.
Fig. 1A schematic diagram of an elastic hollow cylinder, surrounded by an elastic medium and containing viscous liquid and vibrating coaxial rigid cylinder
Let us introduce Cartesian coordinate system $xyz$ and cylindrical coordinate system $r\theta y$. We place coordinate systems centers in inner cylinder geometric centre in unperturbed state.
The absolute rigid cylinder displacement law has the form of:
where ${z}_{m}$ is the rigid cylinder oscillations amplitude, $\omega $ is the frequency, $t$ is the time.
The damping at the expense of viscous incompressible liquid is present in the considered oscillation system. As result, transient processes will go out in the course of time. That is why, we will restrict to the consideration the stationary harmonic oscillations study [24].
In accordance with the approach, suggested in [18], we study the deflection of a hollow cylinder, surrounded by an elastic medium, by means of onemass system. In so doing we consider it in the form of a hollow rigid cylinder with equivalent mass ${m}^{*}$ and an elastic connection with equivalent stiffness coefficient $n$. We find the equivalent mass ${m}^{*}$ under the condition of distributed and onemass systems kinematic energy coincidence, the equivalent stiffness coefficient $n$ being defined by means of hollow cylinder maximal deflection ${w}_{max}$ under uniform loading, equivalent to one Newton force (i.e. ${F}^{*}=$ 1 N):
Here $w$ is the hollow cylinder deflection form function, ${w}_{a}$ is the hollow cylinder deflection function value at the mass concentration point of onemass system, ${m}_{c}$ is the distributed mass per hollow cylinder length unit.
Let us choose the bending oscillations form to be sufficiently close to the expected one and satisfying the boundary conditions at the edges. The deflection of an annular cross section beam, based on Winkler foundation under the uniformly distributed load is chosen as such a form. To define this deflection, we present it in the form of a trigonometric series by longitudinal coordinate, limited by its first member. By taking into account that the equivalent mass is concentrated in the centre (i.e. at the point $y=$ 0), we obtain:
where ${\rho}_{0}$ is the hollow cylinder material density, $J$ is the moment of inertia of hollow cylinderbeam crosssection, $E$ is the Young’s modulus of outer cylinder material, $k$ is the elastic coefficient of Winkler foundation.
The transition to onemass system allows to write down the hollow cylinder motion equation in the form of:
where $N$ is the force, acting on the hollow cylinder from liquid; ${z}_{1}={z}_{1m}{f}_{1}\left(\omega t\right)$ is the displacement law of hollow cylinder, surrounded by elastic medium (within the frameworks of onemass model); $p$ is the liquid pressure.
Viscous incompressible liquid dynamic equations, presented by NavierStocks equations and continuity equation, as well [25] have the form of:
$=\frac{1}{\rho}\frac{\partial p}{\partial r}+\nu \left(\frac{{\partial}^{2}{V}_{r}}{\partial {r}^{2}}\right.+\frac{1}{r}\frac{\partial {V}_{r}}{\partial r}+\left.\frac{1}{{r}^{2}}\frac{{\partial}^{2}{V}_{r}}{\partial {\theta}^{2}}+\frac{{\partial}^{2}{V}_{r}}{\partial {y}^{2}}\frac{2}{{r}^{2}}\frac{\partial {V}_{\theta}}{\partial \theta}\frac{{V}_{r}}{{r}^{2}}\right),$
$\begin{array}{l}\frac{\partial {V}_{\theta}}{\partial t}+{V}_{r}\frac{\partial {V}_{\theta}}{\partial r}+\frac{{V}_{\theta}}{r}\frac{\partial {V}_{\theta}}{\partial \theta}+{V}_{y}\frac{\partial {V}_{\theta}}{\partial y}+{V}_{r}\frac{{V}_{\theta}}{r}\\ =\frac{1}{\rho r}\frac{\partial p}{\partial \theta}+\nu \left(\frac{{\partial}^{2}{V}_{\theta}}{\partial {r}^{2}}\right.+\frac{1}{r}\frac{\partial {V}_{\theta}}{\partial r}+\left.\frac{1}{{r}^{2}}\frac{{\partial}^{2}{V}_{\theta}}{\partial {\theta}^{2}}+\frac{{\partial}^{2}{V}_{\theta}}{\partial {y}^{2}}+\frac{2}{{r}^{2}}\frac{\partial {V}_{r}}{\partial \theta}\frac{{V}_{\theta}}{{r}^{2}}\right),\end{array}$
$\frac{\partial {V}_{y}}{\partial t}+{V}_{r}\frac{\partial {V}_{y}}{\partial r}+\frac{{V}_{\theta}}{r}\frac{\partial {V}_{y}}{\partial \theta}+{V}_{y}\frac{\partial {V}_{y}}{\partial y}=\frac{1}{\rho}\frac{\partial p}{\partial y}+\nu \left(\frac{{\partial}^{2}{V}_{y}}{\partial {r}^{2}}\right.+\frac{1}{r}\frac{\partial {V}_{y}}{\partial r}+\left.\frac{1}{{r}^{2}}\frac{{\partial}^{2}{V}_{y}}{\partial {\theta}^{2}}+\frac{{\partial}^{2}{V}_{y}}{\partial {y}^{2}}\right),$
$\frac{\partial {V}_{y}}{\partial r}+\frac{{V}_{r}}{r}+\frac{1}{r}\frac{\partial {V}_{\theta}}{\partial \theta}+\frac{\partial {V}_{y}}{\partial y}=0.$
Here $v$ is the kinematical coefficient of the liquid viscosity; $\rho $ is the liquid density, ${V}_{r}$, ${V}_{y}$, ${V}_{\theta}$ is the liquid velocity vector projection on coordinate axis.
The boundary conditions of Eq. (5) are the noslip conditions for viscous liquid and the ones for the pressure in the absence of leakage at the edges [18, 25]:
${V}_{r}=\mathrm{c}\mathrm{o}\mathrm{s}\theta {z}_{m}\frac{df}{dt},{V}_{\theta}=\mathrm{s}\mathrm{i}\mathrm{n}\theta {z}_{m}\frac{df}{dt},{V}_{y}=0,r=R,$
$\frac{\partial p}{\partial y}=0,y=\pm \frac{l}{2}.$
3. Solutions and discussions
Let us introduce dimensionless variables and the small parameters:
Here $\psi $, $\lambda $ are the small parameters.
As in the considered setting $\psi =o\left(1\right)$, $\lambda =\u043e\left(1\right)$ and ${z}_{1m}/{z}_{m}=O\left(1\right)$, the member at $\psi $ and $\lambda $ can be omitted in the liquid dynamics Eq. (5) and in the corresponding boundary conditions Eq. (6), written down in variables Eq. (7). As a result, we obtain liquid dynamics problem:
$\frac{\partial {U}_{\xi}}{\partial \xi}+\frac{\partial {U}_{\theta}}{\partial \theta}+\frac{\partial {U}_{\zeta}}{\partial \zeta}=0.$
With boundary conditions:
${U}_{\xi}=\left(\frac{df}{d\tau}\right)\mathrm{c}\mathrm{o}\mathrm{s}\theta ,\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}{U}_{\theta}=0,{U}_{\zeta}=0,\xi =0,\frac{\partial P}{\partial \zeta}=0,\zeta =\pm 1.$
The solution of the problem Eq. (8) with boundary conditions Eq. (9) for the established harmonic oscillations regime takes the form of:
${U}_{\xi}=\frac{1}{2{\epsilon}^{2}}\left[\frac{\partial}{\partial \tau}\left(\frac{{\partial}^{2}P}{\partial {\theta}^{2}}+\frac{1}{{\sigma}^{2}}\frac{{\partial}^{2}P}{\partial {\zeta}^{2}}\right)\left({\stackrel{}{\mathrm{\Psi}}}_{1}\left(\xi \right)\xi \right)+\left(\frac{{\partial}^{2}P}{\partial {\theta}^{2}}+\frac{1}{{\sigma}^{2}}\frac{{\partial}^{2}P}{\partial {\zeta}^{2}}\right){\stackrel{}{\mathrm{\Phi}}}_{1}\left(\zeta \right)\right]+\frac{df}{d\tau}\mathrm{c}\mathrm{o}\mathrm{s}\theta ,$
$P=2{\epsilon}^{2}\alpha \left(\frac{{d}^{2}f}{d{\tau}^{2}}\frac{{z}_{1m}}{{z}_{m}}\frac{{d}^{2}{f}_{1}}{d{\tau}^{2}}\right)\mathrm{c}\mathrm{o}\mathrm{s}\theta +12\gamma \left(\frac{df}{d\tau}\frac{{z}_{1m}}{{z}_{m}}\frac{d{f}_{1}}{d\tau}\right)\mathrm{c}\mathrm{o}\mathrm{s}\theta .$
Here we introduce the symbols:
$\alpha \left(\omega \right)=\frac{\mathrm{\epsilon}\left(\mathrm{\epsilon}\left(\mathrm{c}\mathrm{h}\epsilon +\mathrm{c}\mathrm{o}\mathrm{s}\epsilon \right)\left(\mathrm{s}\mathrm{h}\epsilon +\mathrm{s}\mathrm{i}\mathrm{n}\epsilon \right)\right)}{{\epsilon}^{2}\left(\mathrm{c}\mathrm{h}\epsilon +\mathrm{c}\mathrm{o}\mathrm{s}\epsilon \right)2\epsilon \left(\mathrm{s}\mathrm{h}\epsilon +\mathrm{s}\mathrm{i}\mathrm{n}\epsilon \right)+2\left(\mathrm{c}\mathrm{h}\epsilon \mathrm{c}\mathrm{o}\mathrm{s}\epsilon \right)},$
$\stackrel{}{\mathrm{\Psi}}\left(\xi \right)={F}_{2}\left(\epsilon \xi \right){D}_{1}{F}_{1}\left(\epsilon \xi \right)2{F}_{4}\left(\epsilon \xi \right){D}_{2},\stackrel{}{\mathrm{\Phi}}\left(\xi \right)=2{F}_{3}\left(\epsilon \xi \right){F}_{2}\left(\epsilon \xi \right){D}_{2}2{F}_{4}\left(\epsilon \xi \right){D}_{1},$
${F}_{1}\left(\mathrm{\epsilon}\mathrm{\xi}\right)=\mathrm{c}\mathrm{h}\epsilon \xi \mathrm{c}\mathrm{o}\mathrm{s}\epsilon \xi ,\mathrm{}\mathrm{}\mathrm{}{F}_{2}\left(\epsilon \xi \right)=\frac{1}{2}\left[\mathrm{c}\mathrm{h}\epsilon \xi \mathrm{s}\mathrm{i}\mathrm{n}\epsilon \xi +\mathrm{s}\mathrm{h}\epsilon \xi \mathrm{c}\mathrm{o}\mathrm{s}\epsilon \xi \right],{F}_{3}\left(\epsilon \xi \right)=\frac{1}{2}\mathrm{s}\mathrm{h}\epsilon \xi \mathrm{s}\mathrm{i}\mathrm{n}\epsilon \xi ,$
${F}_{4}\left(\epsilon \xi \right)=\frac{1}{4}\left[\mathrm{c}\mathrm{h}\epsilon \xi \mathrm{s}\mathrm{i}\mathrm{n}\epsilon \xi \mathrm{s}\mathrm{h}\epsilon \xi \mathrm{c}\mathrm{o}\mathrm{s}\epsilon \xi \right],{D}_{1}=\frac{\mathrm{s}\mathrm{h}\epsilon \mathrm{s}\mathrm{i}\mathrm{n}\epsilon}{\mathrm{c}\mathrm{o}\mathrm{s}\epsilon +\mathrm{c}\mathrm{h}\epsilon},{D}_{2}=\frac{\mathrm{s}\mathrm{i}\mathrm{n}\epsilon +\mathrm{s}\mathrm{h}\epsilon}{\mathrm{c}\mathrm{o}\mathrm{s}\epsilon +\mathrm{c}\mathrm{h}\epsilon},$
$\epsilon \left(\omega \right)=\sqrt{\frac{{\delta}^{2}\omega}{2\nu}},{\stackrel{}{\mathrm{\Psi}}}_{1}\left(\xi \right)=\frac{1}{2\epsilon}\left[\frac{\mathrm{s}\mathrm{h}\epsilon \xi \mathrm{s}\mathrm{i}\mathrm{n}\epsilon \xi (\mathrm{s}\mathrm{i}\mathrm{n}\epsilon \mathrm{s}\mathrm{h}\epsilon )\mathrm{c}\mathrm{h}\epsilon \xi \mathrm{c}\mathrm{o}\mathrm{s}\epsilon \xi (\mathrm{s}\mathrm{i}\mathrm{n}\epsilon +\mathrm{s}\mathrm{h}\epsilon )}{\mathrm{c}\mathrm{h}\epsilon +\mathrm{c}\mathrm{o}\mathrm{s}\epsilon}\right.$
$\left.+\frac{\mathrm{s}\mathrm{i}\mathrm{n}\epsilon +\mathrm{s}\mathrm{h}\epsilon}{\mathrm{c}\mathrm{h}\epsilon +\mathrm{c}\mathrm{o}\mathrm{s}\epsilon}+\mathrm{c}\mathrm{h}\epsilon \xi \mathrm{s}\mathrm{i}\mathrm{n}\epsilon \xi +\mathrm{s}\mathrm{h}\epsilon \xi \mathrm{c}\mathrm{o}\mathrm{s}\epsilon \xi \right],$
$\left.\frac{\mathrm{s}\mathrm{i}\mathrm{n}\epsilon \mathrm{s}\mathrm{h}\epsilon}{\mathrm{c}\mathrm{h}\epsilon +\mathrm{c}\mathrm{o}\mathrm{s}\epsilon}\mathrm{c}\mathrm{h}\epsilon \xi \mathrm{s}\mathrm{i}\mathrm{n}\epsilon \xi +\mathrm{s}\mathrm{h}\epsilon \xi \mathrm{c}\mathrm{o}\mathrm{s}\epsilon \xi \right].$
In special cases we have ${\stackrel{}{\mathrm{\Psi}}}_{1}\left(0\right)=$ 0, ${\stackrel{}{\mathrm{\Phi}}}_{1}\left(0\right)=$ 0, $\stackrel{}{\mathrm{\Psi}}\left(1\right)=$ –1, $\stackrel{}{\mathrm{\Psi}}\left(0\right)=$ –1, $\stackrel{}{\mathrm{\Phi}}\left(0\right)=$ 0, $\stackrel{}{\mathrm{\Phi}}\left(1\right)=$ 0.
Let us note that in the case of strongly viscous liquid, we can omit inertia members in the Eq. (8). In this case we must take $\alpha =$ 0 and $\gamma =$ 1 at the pressure expression. As a result, the expressions for vector velocity liquid movement component have the form:
${U}_{\theta}=\frac{1}{2}\frac{\partial P}{\partial \theta}\left({\xi}^{2}\xi \right),{U}_{\zeta}=\frac{1}{2{\sigma}^{2}}\frac{\partial P}{\partial \zeta}\left({\xi}^{2}\xi \right).$
By substituting Eq. (10) into Eq. (7), we obtain the expression for dimensional pressure. With taking into account the obtained expression for dimensional pressure, Eq. (4) will take the form of:
where $M=m\nu {R}_{1}(\omega {\delta}^{3}{)}^{1}2{\epsilon}^{2}\alpha $, $K=m\nu {R}_{1}{\delta}^{3}12\gamma $, $m=\pi {R}^{2}l\rho $.
We write down Eq. (11) solution for the stationary harmonic oscillations regime in the form of:
Here $A\left(\omega \right)=\sqrt{\left({M}^{2}{\omega}^{4}+{K}^{2}{\omega}^{2}\right)/{\left(n\left(M+{m}^{*}\right){\omega}^{2}\right)}^{2}+{K}^{2}{\omega}^{2}}$ is the amplitude frequency characteristic of a hollow cylinder, surrounded by an elastic medium and containing viscous liquid, as well as vibrating cylinder, $\phi =\mathrm{a}\mathrm{r}\mathrm{c}\mathrm{t}\mathrm{g}\left(K\omega /\left((M+{m}^{*}){\omega}^{2}n\right)\right)\mathrm{a}\mathrm{r}\mathrm{c}\mathrm{t}\mathrm{g}\left(K/M\omega \right)$.
We introduce the result of the hollow cylinder amplitude frequency characteristics with the following parameters, for example: ${R}_{1}=$ 5∙10^{–2} m, $l=$ 6∙10^{–1} m, ${h}_{0}=$ 4∙10^{3} m, $\rho =$1.84∙10^{3} kg/m^{3}, $v=$ 2.5·10^{–4} m^{2}/s, $E=$ 1.96∙10^{11} Pa, ${\rho}_{0}=$ 7.87∙10^{3} kg/m^{3}. In the cases of taking into account surrounding elastic medium we adopt elastic coefficient of Winkler foundation $k=$ 10^{8} N/m^{3}. While in the cases of elastic medium neglecting we adopt $k=$ 0. In the course of calculation, the cases of neglecting inertia members in liquid dynamics equations we take $\alpha =$ 0 and $\gamma =$ 1 consequently $M=0$, $K=12m\nu {R}_{1}{\delta}^{3}\gamma $. The calculation results are presented in Fig. 2.
The solid blue line is the model for $k=$ 10^{8} N/m^{3 }with inertia of viscous liquid being taken into account; the dash blue line is the model for $k=$ 0 with inertia of viscous liquid being taken into account; the solid black line is the model for $k=$ 10^{8} N/m^{3} and without inertia of viscous liquid being taken into account; the dash black line is the model for $k=$ 0 and without inertia of viscous liquid being taken into account.
4. Comparison with the experimental data
Reference [26] presents the data of experimental research of the amplitudes of oscillation rigid cylinder interacting with viscous liquid. This cylinder has an elastic suspension and is inside the elastic cylinder filled with a viscous fluid. In other words, the rigid cylinder forms an inner wall of the annular channel, while an elastic one forms its outer wall. The annular channel formed by the cylinders is set on the foundation vibrating in the vertical plane. By using our approach, we can work out the mathematical model and make the comparison of model calculation and experimental data. To do this, we consider outer elastic cylinder in the form of a rigid one with equivalent mass and an elastic connection with equivalent stiffness coefficient. The equivalent mass and equivalent stiffness coefficient are defined by Eqs. (2), (3). Then the scheme of mechanical system which was investigated in [26] can be presented as in Fig. 3.
Fig. 2The calculation results
Fig. 3Scheme of mechanical system
Within this calculation scheme the equations of the annular channel walls movement take the form of:
Here ${\ddot{z}}_{0}$ is the vibrating foundation acceleration; ${z}_{1}$ is the outer cylinder movement law; ${z}_{2}$ is the inner cylinder movement law; ${m}_{1}$ is equivalent mass of outer cylinder, calculated according to Eq. (3); ${m}_{2}$ is the inner cylinder mass; ${n}_{1}$ is the equivalent stiffness coefficient of outer cylinder suspension, calculated according to Eq. (3) where $k=$ 0; ${n}_{2}$ is the stiffness coefficient of inner cylinder suspension (electronic spring); $N$ is the viscous liquid layer between cylinder reaction defined by Eq. (4).
Taking into consideration foundation acceleration ${\ddot{z}}_{0}$ we write down the expression for pressure in Eq. (7) as $p={p}_{0}+\rho \nu \lambda \omega {\psi}^{2}P\rho {R}_{2}{\ddot{z}}_{0}\mathrm{c}\mathrm{o}\mathrm{s}\theta $. Further, we find expression for $N$ out of the solutions of hydrodynamics Eq. (10) in the form of:
Here the expressions for $m$, $M$, $K$ coincide with Eq. (11).
By substituting Eq. (14) into Eq. (13) and solving them for stationary harmonic oscillations regime under ${\ddot{z}}_{0}={z}_{0m}{\omega}^{2}\mathrm{s}\mathrm{i}\mathrm{n}\omega t$, we find the inner cylinder movement law:
$a=\left({n}_{1}{m}_{1}{\omega}^{2}\right)\left({n}_{2}\left({m}_{2}+M\right){\omega}^{2}\right)\left({n}_{2}{m}_{2}{\omega}^{2}\right)M{\omega}^{2},$
$b=K\omega \left({n}_{1}+{n}_{2}\left({m}_{1}+{m}_{2}\right){\omega}^{2}\right),$
${a}_{2}=\left({m}_{2}m\right)\left({n}_{1}{m}_{1}{\omega}^{2}\right)M{\omega}^{2}\left({m}_{1}+{m}_{2}\right),{b}_{2}=K\omega \left({m}_{2}+{m}_{1}\right).$
Here ${z}_{0m}{\omega}^{2}$ is the acceleration amplitude of vibrating foundation.
Fig. 4 presents the amplitude oscillations of the inner cylinder under various frequencies of the foundation vibrations; the latter being defined in experimentally in [26] and calculated by Eq. (15) under the foundation vibroacceleration amplitude ${z}_{0m}{\omega}^{2}=$ 1 g. In the course of calculating the following parameters we defined according to [26]: ${\rho}_{0}=$ 8.4·10^{3} kg/m^{3}, $\rho =$ 2·10^{3} kg/m^{3}, $\nu =$ 3.5·10^{6} m^{2}/s, $E=$ 9.81·10^{10} Pа, $l=$ 0.08 m, ${\delta}_{0}/l=$ 6.25·10^{3}, ${h}_{0}/l=$ 1.25·10^{1}, ${R}_{2}/l=$ 3.125·10^{1}, ${R}_{1}/l=$ 31.875·10^{2}, ${m}_{2}=$ 225·10^{3} kg, ${n}_{2}=$ 1.75·10^{4}^{}kg/s^{2}.
The comparison of experimental data with calculation ones by the mathematical model worked out by authors presents a good coincidence. This fact proves that the suggested approach to the consideration of the cylinder elastic properties on the basis of onemass model with equivalent mass and stiffness coefficient is quite legitimate.
Fig. 4The curve 1 is calculation data by Eq. (15); the curve 2 is experimental data from [26]
5. Conclusions
Thus, mathematical model for investigating bending oscillations of a hollow cylinder, surrounded by an elastic medium, filled with viscous liquid and containing vibrating rigid cylinder is suggested. The calculations showed that taking into account only the surrounding elastic medium leads to increasing oscillation resonance frequencies. The resonance frequencies reduce, and amplitude oscillations increase at resonance frequencies are caused by taking into account only liquid movement inertia. Taking into account medium elasticity and viscous incompressible liquid movement inertia simultaneously, leads to further increase resonance frequencies values and the corresponding oscillations amplitudes. Therefore, in the course of modeling real mechanical systems it is necessary to take into consideration viscous incompressible liquid inertia inside hollow cylinder, as well as surrounding elastic medium presence.
References

Gorshkov A. G., Morozov V. I., Ponomarev A. T., Shklyarchuk F. N. Aerogidroelasticity of Designs. Fizmathlit, Moscow, 2000, (in Russian).

Avramov K. V., Strel’nikova E. A. Chaotic oscillations of plates interacting on both sides with a fluid flow. International Applied Mechanics, Vol. 50, Issue 3, 2014, p. 303309.

Amabili M. Vibrations of circular plates resting on a sloshing liquid: solution of the fully coupled problem. Journal of Sound and Vibration, Vol. 245, Issue 2, 2001, p. 261283.

Askari E., Jeong K.H., Amabili M. Hydroelastic vibration of circular plates immersed in a liquidfilled container with free surface. Journal of Sound and Vibration, Vol. 332, Issue 12, 2013, p. 30643085.

Akcabay D. T., Young Y. L. Hydroelastic response and energy harvesting potential of flexible piezoelectric beams in viscous flow. Physics of Fluids, 24, p. 52012.

Faria Cassio T., Inman Daniel J. Modeling energy transport in a cantilevered EulerBernoulli beam actively vibrating in Newtonian fluid. Mechanical Systems and Signal Processing, Vol. 45, 2014, p. 317329.

Mogilevich L. I., Popov V. S., Popova A. A. Dynamics of interaction of elastic elements of a vibrating machine with the compressed liquid layer lying between them. Journal of Machinery Manufacture and Reliability, Vol. 39, Issue 4, 2010, p. 322331.

Mogilevich L. I., Popov V. S. Investigation of the interaction between a viscous incompressible fluid layer and walls of a channel formed by coaxial vibrating discs. Fluid Dynamics, Vol. 46, Issue 3, 2011, p. 375388.

Ageev R. V., Kuznetsova E. L., Kulikov N. I., Mogilevich L. I., Popov V. S. Mathematical model of movement of a pulsing layer of viscous liquid in the channel with an elastic wall. PNRPU Mechanics Bulletin, Vol. 3, 2014, p. 1735.

Mogilevich L. I. Dynamics of interaction between a squeezed layer of a viscous incompressible fluid and an elastic threelayer plate. Mechanics of Solids, Vol. 43, Issue 5, 2008, p. 785792.

Ageev R. V., Mogilevich L. I., Popov V. S. Vibrations of the walls of a slot channel with a viscous fluid formed by threelayer and solid disks. Journal of Machinery Manufacture and Reliability, Vol. 43, Issue 1, 2014, p. 18.

Grushenkova E. D., Mogilevich L. I., Kuznetsova E. L., Popov V. S., Rabinsky L. N. Mathematical model of threelayer plate interaction with viscous incompressible liquid layer under foundation vibration. Applied Mathematical Sciences, Vol. 9, Issue 112, 2015, p. 55515559.

Womersley J. R. Oscillatory motion of a viscous liquid in a thinwalled elastic tubeI: The linear approximation for long waves. Philosophical Magazine, Vol. 46, 1955, p. 199221.

Kondratov D. V., Kondratova J. N., Mogilevich L. I. Oscillating laminar fluid flow in a cylindrical elastic pipe of annular crosssection. Fluid Dynamics, Vol. 44, Issue 4, 2009, p. 528539.

Kondratov D. V., Kondratova J. N., Mogilevich L. I. Studies of the amplitude frequency characteristics of oscillations of the tube elastic walls of a circular profile during pulsed motion of a viscous fluid under the conditions of rigid jamming on the buttends. Journal of Machinery Manufacture and Reliability, Vol. 38, Issue 3, 2009, p. 229234.

Kondratov D. V., Kondratova Y. N., Popov V. S., Plaksina I. V. The problem of a hydroelasticity for a tube ringtype a profile with elastic, geometrically irregular outer shell at pressure influence. Izvestiya of Saratov University. New Series. Series: Mathematics. Mechanics. Informatics, Vol. 13, Issue 3, 2013, p. 7076, (in Russian).

Kondratov D. V., Kondratova J. N., Kuznetsova E. L., Mogilevich L. I., Rabinsky L. N. Mathematical model of elastic ribbed shell dynamics interaction with viscous liquid pulsating layer. Applied Mathematical Sciences, Vol. 9, Issue 71, 2015, p. 35253531.

Mogilevich L. I., Popov V. S. Dynamics of interaction of an elastic cylinder with a layer of viscous incompressible fluid. Mechanics of Solids, Vol. 39, Issue 5, 2004, p. 146155.

Blinkov Yu A., Blinkova Yu A., Mogilevich L. I. Nonlinear waves in coaxial cylinder shells containing viscous liquid inside with consideration for energy dispersion. Computational Continuum Mechanics, Vol. 6, Issue 3, 2013, p. 336345, (in Russian).

Blinkov Yu A., Kuznetsova E. L., Mogilevich L. I., Rabinskiy L. N. Nonlinear waves in viscoelastic physically nonlinear cylindrical shell with a viscous incompressible fluid and surrounded by an elastic medium. Journal on Composite Mechanics and Design, Vol. 21, Issue 2, 2015, p. 251261, (in Russian).

Alekseev V. V., Indeitsev D. A., Mochalova Yu A. Resonant oscillations of an elastic membrane on the bottom of a tank containing a heavy liquid. Technical Physics, The Russian Journal of Applied Physics, Vol. 44, Issue 8, 1999, p. 903907.

Kuznetsova E. L., Mogilevich L. I., Popov V. S., Rabinsky L. N. Mathematical model of the plate on elastic foundation interacting with pulsating viscous liquid layer. Applied Mathematical Sciences, Vol. 10, Issue 23, 2016, p. 11011109.

Mogilevich L. I., Popov V. S., Popova A. A., Christoforova A. V. Mathematical modeling of hydroelastic walls oscillations of the channel on winkler foundation under vibrations. Vibroengineering Procedia, Vol. 8, 2016, p. 294299.

Panovko Y. G., Gubanova I. I. Stability and Oscillations of Elastic Systems. Consultants Bureau Enterprises, New York, 1965.

Loitsyanskii L. G. Mechanics of Liquids and Gases. Pergamon Press, Oxford, 1966.

Konovalov S. F. Theory of Vibration Resistance Accelerometers. Mashinostroenie, Moscow, 1991, (in Russian).
Cited by
About this article
The reported study was funded by Russian Foundation for Basic Research (RFBR) according to the Research Project No 150101604a and President of Russian Federation Grant No. MD6012.2016.8.