Научная статья на тему 'Математическое моделирование динамических процессов в вибрационном кавитаторе'

Математическое моделирование динамических процессов в вибрационном кавитаторе Текст научной статьи по специальности «Физика»

CC BY
152
44
i Надоели баннеры? Вы всегда можете отключить рекламу.
Область наук
Ключевые слова
ГИДРОУПРУГИЕ КОЛЕБАНИЯ / ВЯЗКАЯ ЖИДКОСТЬ / ПЛАСТИНА / ОБОЛОЧКА / МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ / ДАВЛЕНИЕ ЖИДКОСТИ / АМПЛИТУДНО-ЧАСТОТНЫЕ ХАРАКТЕРИСТИКИ / HYDROELASTIC OSCILLATIONS / VISCOUS FLUID / PLATE / MATHEMATICAL MODELING / VIBRATION CAVITATIONS / AMPLITUDE-FREQUENCY CHARACTERISTICS

Аннотация научной статьи по физике, автор научной работы — Попов Виктор Сергеевич, Черненко Александр Викторович

В статье рассматривается механическая модель канала, стенки которого образованы упругой однородной пластиной, с различными видами опирания на концах, и жесткой пластиной. Между ними находится слой вязкой несжимаемой жидкости. Представленная механическая модель находится под воздействием гармонического виброускорения (канал установлен на вибрирующем основании). Используя метод приведенной массы, упругая пластина моделируется жесткой стенкой с упругой связью. Построена математическая модель рассматриваемого щелевого канала, которая представляет собой связанную систему уравнений динамики вязкой жидкости, уравнения динамики пластины, с соответствующими краевыми условиями. Решение поставленной математической модели производится методом возмущений в безразмерных переменных. В качестве малых параметров задачи выбраны относительная толщина слоя вязкой несжимаемой жидкости, и относительная амплитуда колебаний пластины как одномассовой системы. Использование метода возмущений дает возможность линеаризовать уравнения динамики жидкости. Найдены выражения для амплитудно-частотных характеристик пластины канала и давления в слое жидкости. Построены и исследованы амплитудные частотные характеристики коэффициента динамичности (отношение амплитудно-частотной характеристики динамического к статическому прогибу) и амплитудные частотные характеристики перемещения пластины и давления в жидкости на резонансных частотах. Произведены численные расчеты поведения амплитудно-частотных характеристик коэффициента динамичности и давления.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по физике , автор научной работы — Попов Виктор Сергеевич, Черненко Александр Викторович

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Mathematical modeling of dynamic processes in vibration cavitator

The article considers a mechanical model of the narrow channel which formed by elastic plate and a rigid one with various boundary conditions at the ends. Channel filled a viscous incompressible fluid and located on a vibration foundation. The dimensionless variables of the hydroelastic problem and its small parameters are suggested. The relative thickness of liquid layer and relative amplitude of channel wall oscillations are taken as small parameters. The elastic plate is modeled as one-mass system by equivalent mass method. The mathematical model of the narrow channel is built. It consist of fluid dynamics equation and plate dynamics equation, as well as appropriate boundary conditions. The hydroelasticity problem is linearized by perturbations method. The law of channel wall oscillations and hydrodynamic parameters distribution in liquid are defined. Frequency dependent function of dynamic pressure distribution along the channel and channel wall amplitude frequency characteristics are received. The calculations of channel wall amplitude frequency characteristics and liquid pressure at the channel center are done.

Текст научной работы на тему «Математическое моделирование динамических процессов в вибрационном кавитаторе»

Интернет-журнал «Науковедение» ISSN 2223-5167 http ://naukovedenie.ru/

Том 9, №1 (2017) http://naukovedenie.ru/vol9-1.php

URL статьи: http://naukovedenie.ru/PDF/15TVN117.pdf

Статья опубликована 20.02.2017

Ссылка для цитирования этой статьи:

Попов В.С., Черненко А.В. Математическое моделирование динамических процессов в вибрационном кавитаторе // Интернет-журнал «НАУКОВЕДЕНИЕ» Том 9, №1 (2017) http://naukovedenie.ru/PDF/15TVN117.pdf (доступ свободный). Загл. с экрана. Яз. рус., англ.

Выполнено при поддержке гранта РФФИ 15-01-01604-a УДК 532.517.2:539.3

Попов Виктор Сергеевич

ФГБОУ ВО «Саратовский государственный технический университет имени Гагарина Ю.А.», Россия, Саратов1

Профессор кафедры «Прикладная математика и системный анализ»

Доктор технических наук E-mail: vic_p@bk.ru ORCID: http://orcid.org/0000-0002-9582-7195 РИНЦ: http://elibrary.ru/author items.asp?id=13615 SCOPUS: https://www.scopus.com/authid/detail.uri?authorId=55627876864

Черненко Александр Викторович

ФГБОУ ВО «Саратовский государственный технический университет имени Гагарина Ю.А.», Россия, Саратов

Магистрант E-mail: 3chav@mail.ru РИНЦ: http://elibrary.ru/author items.asp?id=819753

Математическое моделирование динамических процессов в вибрационном кавитаторе

Аннотация. В статье рассматривается механическая модель канала, стенки которого образованы упругой однородной пластиной, с различными видами опирания на концах, и жесткой пластиной. Между ними находится слой вязкой несжимаемой жидкости. Представленная механическая модель находится под воздействием гармонического виброускорения (канал установлен на вибрирующем основании). Используя метод приведенной массы, упругая пластина моделируется жесткой стенкой с упругой связью. Построена математическая модель рассматриваемого щелевого канала, которая представляет собой связанную систему уравнений динамики вязкой жидкости, уравнения динамики пластины, с соответствующими краевыми условиями. Решение поставленной математической модели производится методом возмущений в безразмерных переменных. В качестве малых параметров задачи выбраны относительная толщина слоя вязкой несжимаемой жидкости, и относительная амплитуда колебаний пластины как одномассовой системы. Использование метода возмущений дает возможность линеаризовать уравнения динамики жидкости. Найдены выражения для амплитудно-частотных характеристик пластины канала и давления в слое жидкости. Построены и исследованы амплитудные частотные характеристики коэффициента динамичности (отношение амплитудно-частотной характеристики динамического к статическому прогибу) и амплитудные частотные характеристики перемещения пластины и

1 410054, Саратов, Саратовская обл., Политехническая ул., 77

давления в жидкости на резонансных частотах. Произведены численные расчеты поведения амплитудно-частотных характеристик коэффициента динамичности и давления.

Ключевые слова: гидроупругие колебания; вязкая жидкость; пластина; оболочка; математическое моделирование; давление жидкости; амплитудно-частотные характеристики

Работа ряда вибрационных машин происходит в условиях взаимодействия их рабочих элементов с жидкостью [1-6]. Вследствие этого, возникает необходимость изучения взаимодействия рабочих элементов машин с жидкостью, находящейся между этими элементами, при условии вибрации элементов (канала) и пульсации давления жидкости [7-11]. В [12] рассмотрена задача об изгибных колебаниях полосы, погруженной в идеальную несжимаемую жидкость, имеющей свободную поверхностью. Хаотические колебания упругой полосы при ее взаимодействии с потоком идеальной несжимаемой жидкости исследованы в [13], а колебания круглой пластины на свободной поверхности идеальной несжимаемой жидкости в [14]. В работе [15] поставлена и решена задача о колебаниях круглой пластины, погруженной в воду со свободной поверхностью. Задача об изгибных колебаниях стенки канала как балки-полоски, взаимодействующей с идеальной жидкостью, заполняющей канал, решена в [16]. В данной работе выполнено исследование причин возникновения вибрационной кавитации в охлаждающей жидкости, применительно к двигателю внутреннего сгорания.

С другой стороны, известны работы, в которых учитывается вязкость жидкости. В [17] выполнено исследование колебаний балки в потоке вязкой жидкости применительно к пьезопреобразователям. В [18] исследованы колебания консольнозакрепленной пластины, находящейся в потоке вязкой несжимаемой жидкости. Исследованию демпфирующих свойств слоя вязкой жидкости, на котором происходят гармонические колебания бесконечно длиной балки посвящена работа [19]. Вместе с тем, представляет интерес исследование колебаний упругой стенки канала, взаимодействующей с вязкой жидкостью применительно к вибрационному кавитатору.

Рассмотрим схему вибрационного кавитатора, который может быть использован для обеззараживания жидкости [20]. Схема кавитатора представлена на рис. 1. Стенки 1, 2 образуют канал, который заполнен вязкой несжимаемой жидкостью 3. Стенка 2 обладает упругой податливостью и на основе метода приведенной массы моделируются одномассовой моделью, т.е. жесткой стенкой с упругим подвесом [4, 11, 21]. Ширина канала 2£ значительно b, его

длины, а толщина слоя жидкости ^0 значительно больше Zlm амплитуды колебания стенки. Введем декартовую систему координат x, у, z. Будем учитывать, что геометрический размер канала b во много раз превосходит размер £, и считать, геометрический размер в направлении оси у неограниченным, то есть перейдем к рассмотрению плоской задачи. Канал располагается на виброосновании, которое совершает гармонические колебания. На торцах канала

наблюдается свободное истечение в полости с давлением слева и справа Ро.

Закон движения вибрирующего основания представим в виде:

z0 = E, sin at

(1)

здесь: Ez - амплитуда колебаний основания; ю - частота; t - время.

.г, г zn = -Eya2 sin at = -kg sin at

Тогда виброускорение основания можно представить как 0 z Л

Здесь Ez = kgla , k - коэффициент виброперегрузки, g - ускорение свободного падения.

Динамика вязкой жидкости, находящейся в канале между стенками плоского канала, описывается системой уравнений Навье-Стокса и уравнением неразрывности [4, 22, 23]:

dux dt

■ + u„

dux

dx

■ + u„

dux

dz

1 dp f d 2ux d2 u

p dx

■ + v

2

dx2

+

dz2

duz dt

- + u

du, dx

■ + u„

duz_

dz

= -zn---— + v

1 dp f d 2uz d 2u

p dz

2

dx2

+

dz7

dur du„

dx dz

= 0

(2)

u , u

где: x' z - проекции вектора скорости жидкости на оси координат; p - давление;

р, v - плотность и коэффициент кинематической вязкости жидкости.

Рисунок 1. Расчетная схема вибрационного кавитатора (составлено авторами)

Краевые условия - условия прилипания жидкости к стенкам и условия свободного истечения на левом и правом торце [2, 4, 6, 21]:

их = 0 , uz = z2 при z = ^0 ; ux = 0 , uz = 0 при z = 0

p = P0-pz0(z-50 -z)

здесь z2 = z2mf2(Pt)

при

x

= -£ P = P0-pz0(z-50 -z2)

при

x = i

(3)

(4)

закон движения стенки 2. Уравнения движения стенки 2 имеет следующий вид:

m2( z 2 + z0) + n2z 2 = F

(5)

здесь: m2 - приведенная масса стенки 2; n2 - приведенная жесткость упругого подвеса стенки 2, определяются по [3, 22, 23]. F - сила, действующая со стороны сдавливаемого слоя жидкости на пластину.

Выражения для силы F запишется как:

>. n

z

ь е

F=-Я q*dxdy

о -е

(6)

здесь

q22 =-p + 2pv(duJdz) -

нормальное напряжение, действующее со стороны

жидкости на стенку [23].

Введем безразмерные переменные

¥

= 5о/е << 1 Л = z2mlSo << 1 Re = 8l®/v т = at Г = x/е С = zjS^. uz = z

'0 . uz = z 2maUC .

= z 2maUd¥ ■ P = Po +PVz 2ma(S¥2) ^ P - Pz 0( z - S0 - z 2)

(7)

Здесь у, X, Re - параметры, характеризующие задачу.

Подставляя (7) в (3)-(6) получаем задачу упругогидродинамики в безразмерном виде, включающую в себя уравнения динамики жидкости

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Re

¥ Re

~sT

-+A

- + U„

dU^ ~dC

DP 2 д2 U^ д2U — + ¥-r+" ■>

дГ дС

r

д Uc

- + А

( DU, DU, Л

U—l+u,—£ r дГ С дС

дР 2

=--+ ¥

дС +

2 д2Uc ¥ -С +

д U

дГ дС

(8)

r +-С = о

дГ дС

и уравнения движения стенок канала

Ыр^2ma f

m2( z 2 + z 0) + n2 z 2 = 2bePo +■

So¥2

П p - 2¥

du, Л

С

дС

dr

(9)

m = 2pSn be где 0 - масса жидкости в канале.

При этом граничные условия (4), (5) запишутся в виде

n = 0, n = dfL

dT при С = ^Шт) ; U r

Uf = 0, U, = 0 с=0

С при ь ;

(10)

р=o r=+i

Р 0 при r - f.

Для тонкого слоя жидкости у << 1. В нулевом приближении по у уравнения (8) и (9)

2 — п

упрощаются, так как в них можно положить ¥ = .

Учитывая, что перемещение стенки 2 значительно меньше толщины слоя жидкости,

можно утверждать, что X = o(1), zfm/z2m =0(1. Тогда в нулевом приближении по X, рассматривая асимптотическое разложение

р=р +яр +..., ^ =Uo + ,,UC =исо +яи0 +...

динамики жидкости в виде уравнений [2, 3, 4]

получим линеаризованную задачу

i

ReUi=_ spL+ dP^=Q SUJL+ ди,0 _ o дт dfc dC2 ' dC ' dfc dC

и граничных условий

j г

US0 =0 UC0 = -f При C = 1. U(0 = 0, UC0 = 0 при C = 0.

? ?

p = 0 при fc = ±1 и уравнения движения абсолютно жесткой стенки канала

1

m2(z2 + z0) + П2z 2 = 2biP0 + ЫР^2т®(дъ¥2)^ j P0dfc

(11)

(12)

(13)

Решение задачи (11) с учетом граничных условий (12) при установившемся гармоническом законе вибрации позволяет записать давление в виде

1 fc

P = -1)

2s2a dLfL + 12у f

dT dT

(14)

где введены обозначения [3, 4]:

r 1 r r — r r + r

s»=Re/2, a(p)=^4", № = -1 s»^^, r =1 + -_-^,r = r__l r + r2 6 r + r2 s(v) s(a)

r3 =-shs(a>)Kchs(v)+coss(v)), r4 = sins(a>)Kchs(v)+coss(v))

Уравнение (13) с учетом решения (14) примет вид

(m2 + M)z2 + 2Kz2 + n2z2 -Mzl -2Kzl = 2bip0 -m2'z0

(15)

здесь введены обозначения по [4, 13] 2K = 8ibPv(^2) 7 ; M = 4ibPv(v^2)Sa/3. Решение уравнения (5) при гармоническом виброускорении основания (1) имеет вид

z2 = 2bip<) + + arctg(2Kv/((m2 + M)o2 -n2)))]

П2 (16)

лг \ m2

А(а) = ■

J\(n2 -(m + M )rn2 ]2 + [2Kv]2 где v L 2 2 J L J - амплитудные частотные характеристики.

m n u

Приведенная масса 2 и приведенная жесткость 2 определены по методу приведенной

массы [24, 25]. Рассмотрим шарнирное опирание пластины на торцах [26] в этом случае по

m = 17 m n2 = 48 J

m2 - - m0 2 i3

методу приведенной массы будем иметь 35 и 2 .

jJK_

m = 2bh i p j? io

Здесь 0 0 2H0 - масса пластины; E - модуль Юнга материала пластины; 12

- момент инерции пластины; p° - плотность материала пластины.

Для расчета было принято, что стенка канала выполнена из стали с параметрами р0 = 7,87 *103кг/м3 h0 = 0,005м Е = 1,96*10пПа М0 = 0,3 b = 0,27м l = 0,09м

р = 1*103кг/м3 v = 1,006 10-6, м2/с

Расчет амплитудно-частотных характеристик производился в безразмерных величинах

А(а)

в виде А(0) в программе Maple [27], на основе полученных амплитудно-частотных характеристик был построен график в зависимости от a при шарнирном опирании.

Рисунок 2. Зависимость амплитудных частотных характеристик от a

(составлено авторами)

Для сравнения также рассмотрим жесткое защемление пластины на торцах [26] в этом

192 iQnEJ

m =-т mn n7 = 192——

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

- * 2 4 734 0 2 е3

случае по методу приведенной массы будем иметь и 2 , которое

представлено на рисунке 3.

Рисунок 3. Зависимость амплитудных частотных характеристик от a

(составлено авторами)

Принимая во внимание безразмерные переменные (7), давление в размерном виде представляется как

.. , с- ч 1 pva ,,x. 2 2 l

- 2s az 2mf 2+-z

12r df2(at)

a

" 2m

dt

Или

.. , с- ч 1 pva ,,x. 2 14

p = p0 -pz0 (z - 5) ) + 2 (^)2 - 1)

= p0 -pz0(z -50) +

2 д.щ1 l

^((f )2 -1)

2 50щ l

2s2a .. 12/ .

a

22

z +

a

2s 2a

a

z 2 + 12/2

p0 -pz0(z-50) + ((y)2 - 1)[Mpz2 + ^2]

1 pv 2s2a pvas2 1 pv

M =—-= ———, K =—12/

p О .,.2 с- ,,2 ' p ~ ~ 2 /

(17)

где _ a a 2

2 a a

С учетом вида решения (16), давление можно представить в виде:

2

p = p0 -pz0(z-S0) + kgА(a)((X)2 -1)[-Mpa2sin(at + p2) + Kpacos(at + p2)]. Или

x

p = po -pz0 (z - S0 ) + kgAp (a)(^)2 - 1)sin(at + P2 +Ф)

Ap (a) = J{Mp®2 )2 +(2Kpa)2 А(а) atjv где p vv p ' v p ' - АЧХ давления,

, Kpa Kp tgP = - -

Mp a Mp a

В итоге

p = po -pz0 (z - S0 ) + kgAp (a)((T)2 - 1)sin(at + Pi + p)

l

Далее будем строить АЧХ давления для определения численного значения величин.

(18)

Ap (a) = yl(Mpa2 )2 +(2Kpa)2 А(a)

А{о) =

m.

^[(n2 -(m + M )a2 ]2 +[2Ka]2

T^ A/IVf'> Н'Ь I j^j /It/ I I U--14W ~

1де vlv^v^ ^jl j- была найдена ранее.

Т.е. получаем АЧХ давления для построения и исследования в виде

г>/(Mpa2 )2 +(2Kpa)2

Ap (a) =

m^ (Mpa) +(2Kpa) 4(n2 -(m + M )a2 ]2 +[2Ka]2

Рассмотрим шарнирное опирание пластины на торцах [26] в этом случае по методу

17 n = 48 EJ

m = — m0 n2 = 48^3 приведенной массы будем иметь 35 и 2 .

j=bhL

m = 2bh е р j? 1о

Здесь 0 0 2Н0 - масса пластины; E - модуль Юнга материала пластины; 12

- момент инерции пластины; Po - плотность материала пластины.

г> Ар (a)

Расчет амплитудно-частотных характеристик давления производился в виде р , для

определения численных величин данного параметра в программе Maple [27], на основе

полученных амплитудно-частотных характеристик давления были построены графики в

зависимости от a при шарнирном опирании, которое представлено на рисунке 4.

Рисунок 4. Зависимость амплитудных частотных характеристик от a

(составлено авторами)

Для сравнения также рассмотрим жесткое защемление пластины на торцах [26] в этом

192

случае по методу приведенной массы будем иметь представлено на рисунке 5.

EJ

m =-4 m0 щ =192

и 2 , которое

4,73

Лр(а>)

1400D

12DDQ

10000

в о о о

■6DQQ

40DD

О 2DD 400 ЙОО SOO 1000 12DD 1400 16D0 1S00 2D DO 22DD 2-100

Рисунок 5. Зависимость амплитудных частотных характеристик от a

(составлено авторами)

Решения и расчеты АЧХ перемещений стенки канала и давления в жидкости с помощью программы Maple позволяют сделать вывод, что на значение резонансных частот и амплитуд колебаний оказывает большое влияние вид опирания пластины. Расчеты АЧХ давления в жидкости показывают, что численное значение давления в центре канале превышает давления насыщенного пара на два-три порядка. Таким образом, расчетным путем показана возможность возбуждения путем вибрации кавитационного поля для обеззараживания жидкости, и возможность найти частоту вибрации - резонансную частоту, на которой возможен этот процесс, а также определить условия для проведения натурного эксперимента.

По полученным результатам при моделировании колебаний стенки и давления в центре канала вибрационного кавитатора на базе одномассовой модели можно сделать вывод, что разработанная математическая модель имеет практическое значение. В частности, предложенная математическая модель упругих колебаний стенок, может быть актуальна для создания безреагентных технологий очистки на базе эффекта вибрационной кавитации.

2D00--

D-

со, радi с

ЛИТЕРАТУРА

1. Башта Т.М. Машиностроительная гидравлика. М.: Машиностроение, 1971. 672 с.

2. Попова А.А. Математическое моделирование динамических процессов в виброопоре с упругими элементами конструкции // Вестник Саратовского государственного технического университета. 2007. Т. 1. №4. С. 25-31.

3. Могилевич Л.И., Попов В.С. Прикладная гидроупругость в машино- и приборостроение. Саратов: Саратовский ГАУ, 2003. 156 с.

4. Могилевич Л.И., Попов В.С., Попова А.А. Динамика взаимодействия упругих элементов вибромашины со сдавливаемым слоем жидкости, находящимся между ними // Проблемы машиностроения и надежности машин. 2010. №4. С. 23-32.

5. Могилевич Л.И., Кондратов Д.В. Упругогидродинамика машин и приборов на транспорте. М.: РГОТУПС, 2007. 168 с.

6. Попова А.А. Математическое моделирование процессов взаимодействия вязкой жидкости с тонкостенными ребристыми элементами гидродинамических демпферов и трубопроводов. Диссертация на соискание ученой степени кандидата технических наук / Саратовский государственный технический университет. Саратов, 2008. 174 с.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

7. Агеев Р.В., Могилевич Л.И., Попов В.С., Попова А.А. Движение вязкой жидкости в плоском канале, образованном вибрирующим штампом и шарнирно опертой пластиной // Труды МАИ. 2014. №78.

8. Агеев Р.В., Кузнецова Е.Л., Куликов Н.И., Могилевич Л.И., Попов В.С. Математическая модель движения пульсирующего слоя вязкой жидкости в канале с упругой стенкой // Вестник Пермского национального исследовательского политехнического университета. Механика. 2014. №3. С. 1735.

9. Агеев Р.В., Могилевич Л.И., Попов В.С. Колебания стенок щелевого канала с вязкой жидкостью, образованного трехслойным и твердым дисками // Проблемы машиностроения и надежности машин. 2014. №1. С. 3-11.

10. 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. 2015. Т. 9. №112. С. 5551-5559.

11. Могилевич Л.И., Попов В.С. Исследование взаимодействия слоя вязкой несжимаемой жидкости со стенками канала, образованного соосными вибрирующими дисками // Известия Российской академии наук. Механика жидкости и газа. 2011. №3. С. 42-55.

12. Haddara M.R. and Cao S. A Study of the Dynamic Response of Submerged Rectangular Flat Plates // Marine Structures. Vol. 9. №10. 1996. P. 913-933.

13. Аврамов К.В., Стрельникова Е.А. Хаотические колебания пластинок при их двустороннем взаимодействии с потоком движущейся жидкости // Прикладная механика. 2014. Т. 50. №3. С. 86-93.

14. Amabili M. Vibrations of Circular Plates Resting on a Sloshing Liquid: Solution of the Fully Coupled Problem // Journal of Sound and Vibration. 2001. Vol. 245. №2. P. 261283.

15. 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. 2013. Vol. 332. №12. P. 3064-3085.

16. Индейцев Д.А., Полипанов И.С., Соколов С.К. Расчет кавитационного ресурса втулки судовых двигателей // Проблемы машиностроения и надежности машин. 1994. №4. С. 59-64.

17. Akcabay D.T., Young Y.L. Hydroelastic Response and Energy Harvesting Potential of Flexible Piezoelectric Beams in Viscous Flow // Physics of Fluids. 2012, Vol. 24. №5.

18. 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. №2. 2014. P. 317-329.

19. Onsay T. Effects of layer thickness on the vibration response of a plate-fluid layer system // Journal of Sound and Vibration. 1993. Vol. 163. №2. P. 231-259.

20. Черненко А.В. Обеззараживание жидкости на основе эффекта вибрационной кавитации // Совершенствование методов гидравлических расчетов водопропускных и очистных сооружений. 2014. Т. 1. №1 (41). С. 88-92.

21. Ageev R.V., Mogilevich L.I., Popov V.S., Popova A.A., Kondratov D.V. Mathematical model of pulsating viscous liquid layer movement in a flat channel with elastically fixed wall // Applied Mathematical Sciences. 2014. Т. 8. №157-160. С. 7899-7908.

22. Козлов В.Г. Устойчивость периодического движения жидкости в плоском канале // Изв. АН СССР. МЖГ. 1979. №6. С. 24-32.

23. Лойцянский Л.Г. механика жидкости и газа. М.: Дрофа, 2003. 840 с.

24. Курс сопротивления материалов // М.М. Филоненко-Бородич, С.И. Изюмов, В.А. Олисов, И.Н. Кудрявцев, Л.И. Мальгинов. М.: Гостехиздат. Ч. 1, 1955. 644 с., Ч. 2, 1956. 539 с.

25. Могилевич Л.И., Попов В.С. Динамика взаимодействия упругого цилиндра со слоем вязкой несжимаемой жидкости // Известия Российской академии наук. Механика твердого тела. 2004. №5. С. 179.

26. Попова А.А. Взаимодействие шарнирно закрепленной пластины со слоем вязкой жидкости // Научные труды SWorld. 2014. Т. 29. №4. С. 3-5.

27. Говорухин В.Н., Цибулин В.Г. Введение в Maple V. Математический пакет для всех. М.: Мир, 1997.

Popov Victor Sergeevich

Yuri Gagarin state technical university of Saratov, Russia, Saratov

E-mail: vic_p@bk.ru

Chernenko Aleksandr Viktorovich

Yuri Gagarin state technical university of Saratov, Russia, Saratov

E-mail: 3chav@mail.ru

Mathematical modeling of dynamic processes in vibration cavitator

Abstract. The article considers a mechanical model of the narrow channel which formed by elastic plate and a rigid one with various boundary conditions at the ends. Channel filled a viscous incompressible fluid and located on a vibration foundation. The dimensionless variables of the hydroelastic problem and its small parameters are suggested. The relative thickness of liquid layer and relative amplitude of channel wall oscillations are taken as small parameters. The elastic plate is modeled as one-mass system by equivalent mass method. The mathematical model of the narrow channel is built. It consist of fluid dynamics equation and plate dynamics equation, as well as appropriate boundary conditions. The hydroelasticity problem is linearized by perturbations method. The law of channel wall oscillations and hydrodynamic parameters distribution in liquid are defined. Frequency dependent function of dynamic pressure distribution along the channel and channel wall amplitude frequency characteristics are received. The calculations of channel wall amplitude frequency characteristics and liquid pressure at the channel center are done.

Keywords: hydroelastic oscillations; a viscous fluid; a plate; mathematical modeling; vibration cavitations; amplitude-frequency characteristics

REFERENCES

1. Basta T.M. Engineering hydraulics. Moscow: Mashinostroenie, 1971. 672.

2. Popov A.A. Mathematical modeling of dynamic processes in vibration support with elastic parts of constraction // Bulletin of Saratov state technical University. 2007. Vol. 1. No. 4. S. 25-31.

3. Mogilevich L.I., Popov V.S. Applied Hydroelasticity of Mashino-and Instrument Making: Saratov: Saratov State Agrarian University, 2003. p. 156.

4. Mogilevich L.I., Popov V.S., Popova A.A. Dynamics of the interaction of elastic elements with the vibrator to compress a layer of fluid between them // Journal of Machinery Manufacture and Reliability. 2010, Volume 39, Issue 4, pp. 322-331.

5. Mogilevich L.I., Kondratov D.V. Oppositionality machines and devices for transport. M.: RGOTUPS, 2007. p. 168.

6. Popov A.A. Mathematical modeling of the interaction of a viscous liquid with thin-walled ribbed elements of hydrodynamic dampers and pipelines. The dissertation on competition of a scientific degree of candidate of technical Sciences, Saratov state technical University. Saratov, 2008. 174 P.

7. Ageev, R.V., Mogilevich L.I., Popov V.S., Popova A.A. Viscous liquid movement in a flat channel, formed by a vibrating stamp and simply supported plate // Journal "Trudy MAI" Moscow Aviation institute (National Research University). 2014. No. 78.

8. 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. 2014. No. 3. Pp. 17-35. DOI: 10.15593/perm.mech/2014.3.02.

9. Ageev, R.V., Mogilevich L.I., Popov V.S. the Vibrations of the walls of a slot channel with a viscous fluid formad by three-layer and solid disks // Journal of Machinery Manufacture and Reliability. 2014, Volume 43, Issue 1, pp. 1-8.

10. Grushenkova E.D., Mogilevich L.I., Popov V.S., L.N. Rabinsky, Kuznetsova E.L. Mathematical model of three-layer plate interaction with viscous incompressible liquid layer under foundation vibration // Applied Mathematical Sciences. 2015. T. 9. No. 112. S. 5551-5559.

11. Mogilevich L.I., Popov V.S. Investigation of the interaction between a viscous incompressible fkuid layer and waals of a channel formad by coaxial vibrating discs // Fluid Dynamics. 2011, Volume 46, Issue 3, pp 375-388.

12. Haddara M.R. and Cao S. A Study of the Dynamic Response of Submerged Rectangular Flat Plates // Marine Structures. Vol. 9. №10. 1996. P. 913-933.

13. Avramov K.V., Strelnikova E.A. Chaotic Vibrations of Plates Two-Sided Interacting with Flux of Moving Fluid // Applied mechanics. 2014. T. 50. №3. p. 86-93.

14. Amabili M. Vibrations of Circular Plates Resting on a Sloshing Liquid: Solution of the Fully Coupled Problem // Journal of Sound and Vibration. 2001. Vol. 245. №2. P. 261283.

15. 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. 2013. Vol. 332. №12. P. 3064-3085.

16. Indeitsev D.A., Polipanov I.S., Sokolov S.K. Problemy mashinostroeniya i nadezhnosti mashin, 1994, no.4, pp. 59-64.

17. Akcabay D.T., Young Y.L. Hydroelastic Response and Energy Harvesting Potential of Flexible Piezoelectric Beams in Viscous Flow // Physics of Fluids. 2012, Vol. 24. №5.

18. 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. №2. 2014. P. 317-329.

19. Onsay T. Effects of layer thickness on the vibration response of a plate-fluid layer system // Journal of Sound and Vibration. 1993. Vol. 163. №2. P. 231-259.

20. Chernenko V.A. Disinfection liquid based on the effect of vibratory cavitation Improvement of methods of hydraulic calculations of culverts and water treatment facilities. 2014. Vol. 1. No. 1 (41). S. 88-92.

21. Ageev R.V., Mogilevich L.I., Popov V.S., Popova A.A., Kondratov D.V. Mathematical model of pulsating viscous liquid layer movement is in a flat channel with elastically fixed wall // Applied Mathematical Sciences. 2014. Vol. 8. No. 157-160. S. 7899-7908.

22. Kozlov V.G. Stability of periodic motion of fluid in a planar channel // Fluid Dynamics. 1979, Volume 14, Issue 6, pp. 904-908.

23. Loytzyansky L.G. Mechanics of liquid and gas. M.: Drofa, 2003. 840 p.

24. The course of strength of materials // M.M. Filonenko-Borodich, S.I. Izyumov, V.A. Olisov, I.N. Kudryavtsev, L.I. Malginov. M: Gostekhizdat. H 1, 1955. 644 p., Part 2, 1956. 539 S.

25. Mogilevich L.I., Popov V.S. Dynamics of interaction of an elastic cylinder with a layer of viscous incompressible fluid // Mech. Solids. 39 (5), 146-155 (2004).

26. Popov A.A. The Interaction of a hinged plate with a layer of a viscous liquid // proceedings SWorld. 2014. Vol. 29. No. 4. S. 3-5.

27. Govorukhin, V.N., Cibulin V.G. Introduction to the Maple Mathematical package for all. M.: Mir, 1997.

i Надоели баннеры? Вы всегда можете отключить рекламу.