Mathematical modeling of hydroelastic walls oscillations of the channel on Winkler foundation under vibrations
Lev I. Mogilevich^{1} , Victor S. Popov^{2} , Anna A. Popova^{3} , Aleftina V. Christoforova^{4}
^{1, 2, 3}Yuri Gagarin State Technical University of Saratov, Saratov, Russia
^{4}Saratov State University, Saratov, Russia
^{2}Corresponding author
Vibroengineering PROCEDIA, Vol. 8, 2016, p. 294-299.
Received 29 May 2016; accepted 24 August 2016; published 7 October 2016
The bending oscillations of a narrow slit channel walls with highly viscous liquid inside and put on a vibrating Winkler foundation are investigated. The channel walls bending oscillations laws are discovered on the basis of hydroelasticity problem solution, as well as pressure in the liquid ones. The deflections amplitudes distribution and liquid pressure along the channel functions are constructed. The obtained results allow investigating dynamic processes, conditioned by constructions elastic elements and viscous liquid interaction in lubrication system, damping and various devices and units.
Keywords: hydroelastic oscillations, viscous liquid, plate, vibrations, Winkler foundation.
1. Introduction
The elastic elements of constructions with liquid oscillations are of great importance for modern machines, units and devices analysis [1]. In reference [2] the rectangular plate oscillations plunged into an ideal incompressible liquid with a free surface are investigated. Reference [3] studies the plate chaotic oscillations under interaction with ideal incompressible flow. The problem of round plate oscillations on incompressible liquid fee surface is considered in reference [4], whereas reference [5] investigates round plate oscillations in water with free surface. The problem of bending oscillations of channel wall interacting with ideal liquid, which fills the channel, is solved in [6], and vibration cavitations in cooling liquid, surrounding internal-combustion engine shell having been solved on its basis. The investigation of the beam hydroelastic oscillation in viscous flow with application to piezoelectric element which can be used for receiving energy from the flow is made in [7]. Reference [8] studies the cantilevering fixed beam oscillations in an unlimited volume of viscous incompressible flow. Reference [9] investigates interminable plate oscillations on viscous liquid layer. The investigation of terminable length vibrating plates interaction with viscous incompressible liquid layer between them is made in [10]. References [11, 12] consider analogues problems for round disks with one of them being three-layered. The plate oscillations, caused by pressure pulsating in viscous liquid layer are considered in [13]. References [14, 15] deal with three-layered plates interaction with viscous liquid layer. The oscillations of membrane with elastic foundation, situated on the channel bottom, with ideal incompressible liquid with a free surface inside the channel are investigated in [16]. However, for practical purposes, the consideration of thin liquid layer viscosity impact on elastic channel walls oscillations, for the case of the channel having vibrating elastic foundation, is of great interest.
2. Statement of the problem
Let us consider a narrow slit channel, filled with highly viscous incompressible liquid (see Fig. 1). The channel walls make two elastic plates 1 and 2, being of similar geometrical sizes in the plan and simply supported at the butt ends. The channel size $b\gg \mathcal{l}$. The plates thicknesses are ${h}_{1}$, ${h}_{2}$, the liquid layer thickness is ${\delta}_{0}\ll \mathcal{l}$. The plate 1 is put on Winkler elastic foundation. The considered mechanical system performs harmonic oscillations with assigned vibro-acceleration ${\ddot{z}}_{0}={E}_{z}{d}^{2}{f}_{0}\left(\omega t\right)/d{t}^{2}=-{E}_{z}{\omega}^{2}{f}_{0}\left(\omega t\right)$, ${f}_{0}\left(\omega t\right)=\mathrm{s}\mathrm{i}\mathrm{n}\omega t$. Here ${E}_{z}{\omega}^{2}$ is the vibro-acceleration amplitude, $\omega $ is the oscillation frequency, $t$ is the time. The vibro-acceleration amplitude can be assigned in gravity acceleration units $g$, i.e. ${E}_{z}{\omega}^{2}=kg$, here $k$ is the vibro-overloading coefficient. Plates 1, 2 perform bending oscillations with amplitudes ${w}_{1m}\ll {\delta}_{0}$, ${w}_{2m}\ll {\delta}_{0}$. In the considered assignment, the liquid flow in the channel is viewed as a creeping one and liquid leakage at the butt ends being free in the same liquid with constant pressure ${p}_{0}$.
Fig. 1. A schematic diagram of a narrow slit channel on Winkler foundation
Let us put Cartesian coordinate system in the middle plate 1 surface center. Taking into account, that $b\gg \mathcal{l}$ we will consider the problem in a flat setting. Strong damping, conditioned by liquid viscosity, takes place in the considered oscillating system. Therefore, transition processes in the system will get out rapidly and the stationary oscillations will emerge. Further, according to reference [17], we will focus on studying stationary forced harmonic oscillations regime.
3. The theory
The channel wall bending oscillations equations take the form of:
where ${w}_{i}$ is the deflection of $i$th plate, ${D}_{i}$ is the bending stiffness of $i$th plate, ${\rho}_{i}$ is the density of $i$th plate, $\chi $ is the elastic coefficient of Winkler foundation, $q=-p+2\rho \nu \partial {u}_{z}/\partial z$ is the normal stress, acting from the liquid layer on the plate [18], $p$ is the pressure, ρ is the density of the liquid, $\nu $ is the kinematical coefficient of the liquid viscosity, ${u}_{z}$ is the liquid velocity vector project on coordinate axis $z$.
Boundary conditions for Eq. (1) are the conditions of simply supported edges:
Highly viscous incompressible liquid dynamics in a narrow slit channel is described by Navier-Stokes equations and continuity equation, in which local and convective inertia members are omitted, according to hydrodynamic lubrication theory [18]:
where ${u}_{x}$ is the liquid velocity vector project on coordinate axis x.
Boundary conditions for Eq. (3) are the conditions of liquid adhesion to the walls and the conditions for pressure at the both ends:
${u}_{x}=\frac{\partial {u}_{1}}{\partial t},{u}_{z}=\frac{\partial {w}_{1}}{\partial t},\text{at}z=\frac{h{}_{1}{}^{\mathrm{}}\mathrm{}}{2}+{w}_{1},$
where ${u}_{i}$ is the longitudinal shift of $i$th plate, $i=$ 1, 2.
Let us introduce dimensionless variables:
${w}_{i}={w}_{im}{W}_{i},\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}{u}_{i}={u}_{im}{U}_{i},\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}p={p}_{0}+P\frac{\nu \rho {w}_{2m}^{}\omega {\psi}^{-2}}{{\delta}_{0}}-\rho {\ddot{z}}_{0}\left(z-{\delta}_{0}-{w}_{2}-\frac{{h}_{1}}{2}\right),$
${\lambda}_{i}=\frac{{w}_{im}^{}}{{\delta}_{0}}\ll 1,\psi =\frac{{\delta}_{0}}{\mathcal{l}}\ll 1,$
where $\psi $, ${\lambda}_{i}$ are the small parameters, characterizing the problem, $i=$1, 2.
By substituting Eq. (6) in liquid dynamic Eq. (3) and the corresponding boundary Eq. (4) and (5), we get:
${U}_{\xi}=\psi \left(\frac{{u}_{1m}}{{w}_{2m}}\right)\frac{\partial {U}_{1}}{\partial \tau},{U}_{\zeta}=\left(\frac{{w}_{1m}}{{w}_{2m}}\right)\frac{\partial {W}_{1}}{\partial \tau}\text{at}\zeta ={\lambda}_{1}{W}_{1},$
The normal stress, acting from the liquid layer on the plate, will be written down, as:
Further, we note that in the considered setting $\psi =o\left(1\right)$, ${u}_{im}/{w}_{2m}=O\left(1\right)$, ${\lambda}_{1}={\lambda}_{2}\left({w}_{1m}/{w}_{2m}\right)=o\left(1\right)$, and therefore, the member at $\psi $, ${\psi}^{2}$, ${\lambda}_{1}$, ${\lambda}_{2}$ in Eq. (7), boundary Eq. (8), (9) and normal stress Eq. (10) can be omitted. As a result, we get liquid dynamic problem in the form of the Eq. (11) and boundary Eq. (12):
$P=0\text{at}\xi =\pm 1$ .
The solutions of the Eq. (11) with boundary Eq. (12) has the form of:
$P=12\underset{1}{\overset{\xi}{\int}}\underset{}{\overset{}{\int}}\frac{\partial \stackrel{~}{W}}{\partial \tau}d\xi d\xi -6\left(\xi -1\right)\underset{-1}{\overset{1}{\int}}\underset{}{\overset{}{\int}}\frac{\partial \stackrel{~}{W}}{\partial \tau}d\xi d\xi ,\stackrel{~}{W}={W}_{2}-\frac{{w}_{1m}}{{w}_{2m}}{W}_{1}.$
By substituting Eq. (13) in the expression from the normal stress Eq. (10), written down in zero approximation at $\psi $ and ${\lambda}_{2}$, and considering the normal stress in channel walls dynamics Eq. (1), we get integro-differential equations of channel walls hydroelastic oscillations:
$=(-1{)}^{i}\left({p}_{0}+12\nu \rho {w}_{2m}^{}\omega ({\psi}^{2}{\delta}_{0}{)}^{-1}\left(\underset{1}{\overset{\xi}{\int}}\underset{}{\overset{}{\int}}\frac{\partial \stackrel{~}{W}}{\partial \tau}d\xi d\xi -\left(\xi -1\right)/2\underset{-1}{\overset{1}{\int}}\underset{}{\overset{}{\int}}\frac{\partial \stackrel{~}{W}}{\partial \tau}d\xi d\xi \right)\right).$
With consideration of boundary Eq. (2), (5) we will look for the assigned equations solution in the form of:
here the upper index 0 means the solution, corresponding to static pressure ${p}_{0}$.
By substituting Eq. (15) into Eq. (14) and setting the remaining members of Eq. (14) in series of trigonometric functions of longitudinal coordinate $\xi $ and then, equating the coefficients at similar trigonometric functions, we obtain the expression for ${R}_{ik}^{0}$:
and the system of ordinary differential equations for ${R}_{ik}\left(\tau \right)$:
${w}_{1m}{a}_{22}^{}\frac{d{R}_{1k}^{}}{d\tau}+{w}_{2m}{a}_{23}^{}{R}_{2k}^{}+{w}_{2m}{a}_{24}^{}\frac{d{R}_{2k}^{}}{d\tau}=-{d}_{2}{E}_{z}{\omega}^{2}\frac{{d}^{2}{f}_{0}}{d{\tau}^{2}}.$
Here, it is taken into account that for the stationary harmonic oscillations regime we have ${d}^{2}{R}_{ik}/d{\tau}^{2}=\u2013{R}_{ik}$ and the following symbols having been introduced:
${a}_{22}^{}={a}_{14}^{}=-{a}_{12}^{},{a}_{23}^{}={D}_{2}{\mathcal{l}}^{-4}(\left(2k-1\right)\pi /2{)}^{4}-{\rho}_{2}^{}{h}_{2}^{}{\omega}^{2},{a}_{24}^{}={a}_{12}^{}=-{a}_{14}^{},$
$2{K}_{k}=12\rho \nu {\psi}_{}^{-2}{\delta}_{0}^{-1}(2/\left(\left(2k-1\right)\pi \right){)}^{2},{d}_{i}=-4(-1{)}^{k}(\left(2k-1\right)\pi {)}^{-1}\left({\rho}_{i}{h}_{i}+\left(2-i\right)\rho {\delta}_{0}\right).$
By solving the system Eq. (17) for the stationary harmonic oscillations regime we find:
here $i=$1, 2 for plates 1 and 2, the following symbols having been introduced:
${F}_{2k}^{}={d}_{2}^{}{a}_{11}^{},{C}_{1k}^{}={C}_{2k}^{}={a}_{12}^{}({d}_{1}^{}+{d}_{2}^{}).$
Considering Eq. (16) and (18), we can present the final expressions for plate deflections as:
$+{E}_{z}{\omega}^{2}{A}_{i}(x,\omega )\mathrm{s}\mathrm{i}\mathrm{n}(\omega t+{\varphi}_{i}(x,\omega \left)\right),$
${A}_{i}=\sqrt{{Q}_{i}^{2}+{N}_{i}^{2}},\mathrm{}\mathrm{}{\varphi}_{i}=\mathrm{a}\mathrm{r}\mathrm{c}\mathrm{t}\mathrm{g}\left(\frac{{N}_{i}}{{Q}_{i}}\right),$
${Q}_{i}=\sum _{k=1}^{\mathrm{\infty}}\frac{{F}_{k}{F}_{ik}+{C}_{k}{C}_{ik}}{{F}_{k}^{2}+{C}_{k}^{2}}\mathrm{c}\mathrm{o}\mathrm{s}\frac{(2k-1)\pi x}{2\mathcal{l}},\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}{N}_{i}=\sum _{k=1}^{\mathrm{\infty}}\frac{{F}_{k}{C}_{ik}-{C}_{k}{F}_{ik}}{{F}_{k}^{2}+{C}_{k}^{2}}\mathrm{c}\mathrm{o}\mathrm{s}\frac{(2k-1)\pi x}{2\mathcal{l}}.$
Taking into account Eqs. (18), (15) in the expression for pressure Eq. (13) we finally obtain the law of dimensional pressure distribution in the channel in the form:
$\mathrm{\Pi}\left(x,\omega \right)=\sqrt{{B}_{k}^{2}+{S}_{k}^{2}},\mathrm{}\mathrm{}\mathrm{}\mathrm{}{\varphi}_{p}(x,\mathrm{\omega})=\mathrm{a}\mathrm{r}\mathrm{c}\mathrm{t}\mathrm{g}\left(\frac{{S}_{k}}{{B}_{k}}\right),$
${B}_{k}=\sum _{k=1}^{\infty}2{K}_{k}\omega {C}_{k}\frac{{F}_{1k}^{}-{F}_{2k}^{}}{{F}_{k}^{2}+{C}_{k}^{2}}\mathrm{c}\mathrm{o}\mathrm{s}\frac{(2k-1)\pi x}{2\mathcal{l}},$
${S}_{k}=\sum _{k=1}^{\infty}2{K}_{k}\omega {F}_{k}\frac{{F}_{1k}^{}-{F}_{2k}^{}}{{F}_{k}^{2}+{C}_{k}^{2}}\mathrm{c}\mathrm{o}\mathrm{s}\frac{(2k-1)\pi x}{2\mathcal{l}}.$
4. Conclusions
The first member of the Eq. (19) is the channel walls static deflection conditioned by static pressure in the liquid ${p}_{0}$, the second member is the deflection, conditioned by foundation vibration. The first component in pressure distribution law Eq. (20) represents static liquid pressure in the channel. The second component is dynamic liquid pressure, conditioned by foundation vibro-acceleration (an analogue of pressure due to the gravity). The third component is pressure, conditioned by channel walls oscillations under foundation vibro-acceleration.
It must be noted that the first and the second pressure components cannot exceed the assigned values. The third component pressure value is defined by the function $\mathrm{\Pi}\left(x,\omega \right)$, which represents pressure distribution frequency dependent function along the channel. Under longitudinal coordinate fixed value $x$, the assigned function is the pressure amplitude characteristic in the assigned channel cross-section. The same can be stated in relation to the deflections amplitudes distribution functions ${A}_{i}\left(x,\omega \right)$, being frequency dependent ones. The functions ${\phi}_{i}\left(x,\omega \right)\text{,}$${\phi}_{p}\left(x,\omega \right)$ are frequency dependent ones of the plates deflections phase shifts distribution and the pressure phase shift along the channel, correspondingly.
The investigation of the above mentioned functions behavior, depending on oscillations frequency, makes it possible to study dynamic processes in the channel under consideration. The obtained expressions for channel walls deflections analysis allows to state that the elastic coefficient of Winkler foundation in a static state causes deflection reduction of the first plate, deflections amplitude dynamics reducing in both plates in a dynamic state. Thus, the paper results can be used for calculation and analysis hydroelastic oscillations of channel walls on elastic foundation and for calculation of hydraulic dampers, pumps elements, hydrodrives, lubrication, cooling and fuel supply systems.
Acknowledgements
The research is made under the financial support of the Russian Foundation for Basic Research (RFFI) Grant 15-01-01604-a.
References
- Gorshkov A. G., Morozov V. I., Ponomarev A. T., Shklyarchuk F. N. Aerouprugost konstruktsii (Aerogidroelasticity of designs), Fizmathlit, Moscow, 2000, (in Russian). [Search CrossRef]
- Haddara M. R., Cao S. A study of the dynamic response of submerged rectangular flat plates. Marine Structures, Vol. 9, Issue 10, 1996, p. 913-933. [Search CrossRef]
- 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. 303-309. [Search CrossRef]
- 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. 261-283. [Search CrossRef]
- Askari E., Jeong K.-H., Amabili M. Hydroelastic vibration of circular plates immersed in a liquid-filled container with free surface. Journal of Sound and Vibration, Vol. 332, Issue 12, 2013, p. 3064-3085. [Search CrossRef]
- Indeitsev D. A., Polypanov I. S., Sokolov S. K. Calculation of cavitation life-time of ship engine liner. Problemy Mashinostraeniya i Nadezhnos’ti Mashin, Vol. 4, 1994, p. 59-64. [Search CrossRef]
- Akcabay D. T., Young Y. L. Hydroelastic response and energy harvesting potential of flexible piezoelectric beams in viscous flow. Physics of Fluids, Vol. 24, Issue 5, 2012. [Search CrossRef]
- Faria Cassio T., Inman Daniel J. Modeling energy transport in a cantilevered Euler-Bernoulli beam actively vibrating in Newtonian fluid. Mechanical Systems and Signal Processing, Vol. 45, 2014, p. 317-329. [Search CrossRef]
- Önsay T. Effects of layer thickness on the vibration response of a plate-fluid layer system. Journal of Sound and Vibration, Vol. 163, 1993, p. 231-259. [Search CrossRef]
- 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. 322-331. [Search CrossRef]
- 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. 375-388. [Search CrossRef]
- Ageev R. V., Mogilevich L. I., Popov V. S. Vibrations of the walls of a slot channel with a viscous fluid formed by three-layer and solid disks. Journal of Machinery Manufacture and Reliability, Vol. 43, Issue 1, 2014, p. 1-8. [Search CrossRef]
- 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. 17-35. [Search CrossRef]
- Mogilevich L. I. Dynamics of interaction between a squeezed layer of a viscous incompressible fluid and an elastic three-layer plate. Mechanics of Solids, Vol. 43, Issue 5, 2008, p. 785-792. [Search CrossRef]
- Grushenkova E. D., Mogilevich L. I., Popov V. S., Rabinsky L. N., Kuznetsova E. L. Mathematical model of three-layer plate interaction with viscous incompressible liquid layer under foundation vibration. Applied Mathematical Sciences, Vol. 9, Issue 2015, 112, p. 5551-5559. [Search CrossRef]
- 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. 903-907. [Search CrossRef]
- Panovko Y. G., Gubanova I. I. Stability and Oscillations of Elastic Systems. Consultants Bureau Enterprises, Inc., New York, N. Y., 1965. [Search CrossRef]
- Loitsyanskii L. G. Mechanics of Liquids and Gases. Pergamon Press, Oxford, 1966. [Search CrossRef]