Lev I. Mogilevich1 , Victor S. Popov2 , Dmitry V. Kondratov3 , Lev N. Rabinskiy4
1, 2Yuri Gagarin State Technical University of Saratov, Saratov, Russia
3Russian Presidential Academy of National Economy and Public Administration, Saratov, Russia
4Moscow Aviation Institute (National Research University), Moscow, Russia
Journal of Vibroengineering, Vol. 19, Issue 8, 2017, p. 5758-5766.
Received 14 January 2017; received in revised form 3 October 2017; accepted 13 October 2017; published 31 December 2017
Copyright © 2017 JVE International Ltd. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
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 one-mass 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.
Keywords: hydroelasticity, viscous liquid, elastic medium, hollow cylinder, oscillations, amplitude frequency characteristic, one-mass system.
The problems of mathematical modeling of elastic construction elements interaction with liquid are of theoretical and practical interest . For example, [2-8] deal with hydroelasticity problems of homogeneous plates, interacting with viscous liquid layer in various settings. Reference  considers the bending oscillations of a flat channel wall, interacting with pulsating viscous liquid. References [10-12] investigate the hydroelastic oscillations of three-layered plates, interacting with viscous pulsating liquid, as well as under inertial excitation. References [13-18] research in the axis-symmetric problem of geometrically regular and ribbed cylinder shells, as well as the ones, forming an annular channel. The non-linear longitudinal deformation wave propagation in cylinder shells, filled with viscous liquid and in axis-symmetric setting is presented in [19, 20]. Reference  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.
Let us consider an elastic hollow cylinder, surrounded by an elastic medium (Fig. 1). The cylinder length is , its thickness is significantly less than inner radius . The elastic cylinder is simply supported at its edges. The rigid cylinder with radius 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 . The elastic medium, surrounding the inner cylinder, can be considered as Winkler foundation.
Fig. 1. A 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 and cylindrical coordinate system . We place coordinate systems centers in inner cylinder geometric centre in unperturbed state.
The absolute rigid cylinder displacement law has the form of:
where is the rigid cylinder oscillations amplitude, is the frequency, 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 .
In accordance with the approach, suggested in , we study the deflection of a hollow cylinder, surrounded by an elastic medium, by means of one-mass system. In so doing we consider it in the form of a hollow rigid cylinder with equivalent mass and an elastic connection with equivalent stiffness coefficient . We find the equivalent mass under the condition of distributed and one-mass systems kinematic energy coincidence, the equivalent stiffness coefficient being defined by means of hollow cylinder maximal deflection under uniform loading, equivalent to one Newton force (i.e. 1 N):
Here is the hollow cylinder deflection form function, is the hollow cylinder deflection function value at the mass concentration point of one-mass system, 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 0), we obtain:
where is the hollow cylinder material density, is the moment of inertia of hollow cylinder-beam cross-section, is the Young’s modulus of outer cylinder material, is the elastic coefficient of Winkler foundation.
The transition to one-mass system allows to write down the hollow cylinder motion equation in the form of:
where is the force, acting on the hollow cylinder from liquid; is the displacement law of hollow cylinder, surrounded by elastic medium (within the frameworks of one-mass model); is the liquid pressure.
Viscous incompressible liquid dynamic equations, presented by Navier-Stocks equations and continuity equation, as well  have the form of:
Here is the kinematical coefficient of the liquid viscosity; is the liquid density, , , is the liquid velocity vector projection on coordinate axis.
The boundary conditions of Eq. (5) are the no-slip conditions for viscous liquid and the ones for the pressure in the absence of leakage at the edges [18, 25]:
Let us introduce dimensionless variables and the small parameters:
Here , are the small parameters.
As in the considered setting , and , the member at and 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:
With boundary conditions:
The solution of the problem Eq. (8) with boundary conditions Eq. (9) for the established harmonic oscillations regime takes the form of:
Here we introduce the symbols:
In special cases we have 0, 0, –1, –1, 0, 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 0 and 1 at the pressure expression. As a result, the expressions for vector velocity liquid movement component have the form:
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 , , .
We write down Eq. (11) solution for the stationary harmonic oscillations regime in the form of:
Here is the amplitude frequency characteristic of a hollow cylinder, surrounded by an elastic medium and containing viscous liquid, as well as vibrating cylinder, .
We introduce the result of the hollow cylinder amplitude frequency characteristics with the following parameters, for example: 5∙10–2 m, 6∙10–1 m, 4∙10-3 m, 1.84∙103 kg/m3, 2.5·10–4 m2/s, 1.96∙1011 Pa, 7.87∙103 kg/m3. In the cases of taking into account surrounding elastic medium we adopt elastic coefficient of Winkler foundation 108 N/m3. While in the cases of elastic medium neglecting we adopt 0. In the course of calculation, the cases of neglecting inertia members in liquid dynamics equations we take 0 and 1 consequently , . The calculation results are presented in Fig. 2.
The solid blue line is the model for 108 N/m3 with inertia of viscous liquid being taken into account; the dash blue line is the model for 0 with inertia of viscous liquid being taken into account; the solid black line is the model for 108 N/m3 and without inertia of viscous liquid being taken into account; the dash black line is the model for 0 and without inertia of viscous liquid being taken into account.
Reference  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  can be presented as in Fig. 3.
Fig. 2. The calculation results
Fig. 3. Scheme of mechanical system
Within this calculation scheme the equations of the annular channel walls movement take the form of:
Here is the vibrating foundation acceleration; is the outer cylinder movement law; is the inner cylinder movement law; is equivalent mass of outer cylinder, calculated according to Eq. (3); is the inner cylinder mass; is the equivalent stiffness coefficient of outer cylinder suspension, calculated according to Eq. (3) where 0; is the stiffness coefficient of inner cylinder suspension (electronic spring); is the viscous liquid layer between cylinder reaction defined by Eq. (4).
Taking into consideration foundation acceleration we write down the expression for pressure in Eq. (7) as . Further, we find expression for out of the solutions of hydrodynamics Eq. (10) in the form of:
Here the expressions for , , coincide with Eq. (11).
By substituting Eq. (14) into Eq. (13) and solving them for stationary harmonic oscillations regime under , we find the inner cylinder movement law:
Here 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  and calculated by Eq. (15) under the foundation vibroacceleration amplitude 1 g. In the course of calculating the following parameters we defined according to : 8.4·103 kg/m3, 2·103 kg/m3, 3.5·10-6 m2/s, 9.81·1010 Pа, 0.08 m, 6.25·10-3, 1.25·10-1, 3.125·10-1, 31.875·10-2, 225·10-3 kg, 1.75·104kg/s2.
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 one-mass model with equivalent mass and stiffness coefficient is quite legitimate.
Fig. 4. The curve 1 is calculation data by Eq. (15); the curve 2 is experimental data from 
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.
The reported study was funded by Russian Foundation for Basic Research (RFBR) according to the Research Project No 15-01-01604-a and President of Russian Federation Grant No. MD-6012.2016.8.