с тяжелым наносным и гидрологическим режимами: дис. ... канд. техн. наук. - Новочеркасск, НИМИ, 1975. - 184 с.
Материал поступил в редакцию 12.05.2015. Кловский Алексей Викторович, аспирант, E-mail: Alexey.Klovskiy@yandex.ru Тел. 8 (903) 541-07-85
Румянцев Игорь Семенович,¡доктор тех-
нических наук,профессор,заслуженный деятель науки РФ, почетный член РААСН Козлов Дмитрий Вячеславович, доктор технических наук, профессор, проректор по инновационному развитию E-mail: kozlovdv@mail.ru Тел. 8 (499) 976-29-62
УДК 502/504 : 627.8 : 532.5 В. А. ФАРТУКОВ
ЗАО «Бюро сервиса и эксплуатации» BSM, г. Москва
М. В. ЗЕМЛЯННИКОВА
Федеральное государственное бюджетное образовательное учреждение высшего образования «Российский государственный аграрный университет — МСХА имени К. А. Тимирязева», г. Москва
СТАЦИОНАРНЫЕ НЕЛИНЕЙНЫЕ КОЛЕБАНИЯ НИЖНЕГО БЬЕФА
В статье рассмотрен случай нелинейных колебаний, возникающих в нижнем бьефе при установившемся режиме течения открытого водного потока. Приводится нелинейное уравнение, в котором энергия водного потока диссипирует при больших амплитудах колебания водной поверхности и генерируется при малых значениях амплитуд колебания. Эта система обладает предельными циклами, колеблющихся около состояния, при котором приток и диссипация энергии сбалансированы и определено наличие бифуркаций векторных полей течения воды. Установлено, что в рассмотренной колебательной системе незатухающие колебания практически могут существовать при наличии некоторого источника энергии, который компенсирует расход энергии, возникший за счет присутствия диссипативных сил. Получено нелинейное уравнение для описания протекающего колебательного процесса в открытом водном потоке нижнего бьефа гидротехнического сооружения при сопряжении бьефов в виде гидравлического прыжка. Уравнение позволяет определять параметры нелинейных колебаний (амплитуда волны, частота колебаний, длина волны) возникающих в нижнем бьефе при установившемся режиме течения потока. Результаты вычислений позволяют осуществить коррекцию размеров гасителей энергии водного потока при решении многих практических задач гидротехнических сооружений.
Структура потока, установившийся режим, диссипация энергии, бифуркация,
автоколебательный процесс.
The article considers a case of non-linear vibrations arising downstream under a stable regime of open water flow. There is given a non-linear equation in which the energy of water flow dissipates at high amplitudes of water surface vibrations and generates at small values of vibration amplitudes. This system possesses limit cycles vibrating near the state under which the inflow and energy dissipation are balanced and availability of bifurcations of vector fields of water flow is determined. It is established that in the considered vibration system undamped vibrations can practically exist at the availability of some source of energy which compensates power consumption arising due to present dis-sipative forces. The non-linear equation is obtained for describing a flowing process in the open water current of the lower pond of the hydraulic structure at ponds conjugation in the way of hydraulic jump. The equation allows determining parameters of non-linear vibrations (wave amplitude, frequency of vibrations,wave length) arising downstream under the stable regime of flow current. The results of calculations make it possible to carry out correction of sizes of energy dissipaters of the water flow when deciding many practical tasks of hydraulic structures.
Flow structure, steady regime, energy dissipation, bifurcation, astable process.
Одной из характерных особенностей режима течения водного потока в нижнем бьефе (для случая сопряжения бьефов в виде гидравлического прыжка) является наличие стоячих волн на поверхности потока. Данный режим течения воды в зоне сопряжения сопровождается волнами с характерной амплитудой и длиной.
Природа их образования и развития нуждается в изучении физики протекающего процесса с последующей количественной оценкой его параметров, которые оказывают существенное влияние как на работу гасителя водной энергии, так и на его (гасителя) размеры.
В основе модели процесса образования стационарных нелинейных колебаний в нижнем бьефе находится система уравнений мелкой воды в приближении Сен-Венана, причем система уравнений Сен-Венана представлена в дивергентной форме, позволяющей охватить так называемые разрывные решения.
Как известно, система уравнений Сен-Венана является следствием уравнений Рейнольдса [1, 2], при этом сглаживание турбулентных пульсаций нормального уровня (то есть тех пульсаций, которые порождены трением на границе «жидкость - омываемая твердая поверхность», а не турбулентных пульсаций возникающих в зонах отрывных течений) возможно при осреднении уравнений На-вье-Стокса.
В зоне гидравлического прыжка, как и во всех других случаях отрывных течений, образуется макротурбулентность, имеющая характерный временной масштаб Тт, который существенно больше временного масштаба Тп, отвечающего обычному уравнению турбулентности без отрывных течений с зависимыми от времени t (Тт > t > Тп) характеристиками прыжкового потока. Так же необходимо заметить, что уравнение сопряженных глубин может быть применено только в случае увеличения периода временного сглаживания.
Тогда используя систему дифференциальных уравнений, описывающих нестационарный режим течения водного потока в зоне прыжкового сопряжения, и, предварительно произведя необходимые подстановки, получили нелинейное уравнение локальной нестационарности (1).
3' 2015
1 -О)
ё
до
Л
<Л
^2 -
дt
дt
= 41
ё - h2 ё - К
= 0;
42'
где к1, й2 - соответственно величина первой и второй сопряженных глубин; 4 , д2 - соответственно величина удельного расхода в зонах первой и второй сопряженных глубин; ис - скорость центра массы объема; А - площадь боковой поверхности гидравлического прыжка при аппроксимации продольного профиля в виде квадратичной параболы.
Необходимо заметить, что в случае увеличения периода временного сглаживания до значения, при котором производные д(vcQ.)/дt и дQ./дt обратятся в ноль, будет иметь место уравнение сопряженных глубин Беланже-Бресса.
Авторы статьи рассмотрели основную систему дифференциальных уравнений нестационарного режима движения водного потока в зоне прыжкового сопряжения, приведенную в [3, 4], представив ее в безразмерном виде, разделив на длину прыжка I и у/1 / ё :
d
Ъ2йЪ2 dt
3ё - Н'
dt
+
42
ё ■ К д ■ \
+
+
(К - V ) = 0;
21 &К
3 dt
= 41 - ^
(1) (2)
где 4 - удельный расход; к1 и к2 - соответственно, первая и вторая сопряженные глубины; К -осредненное значение второй сопряженной глубины; I - длина гидравлического прыжка.
Выразили
2/ 7 ¿Й2
4=4- % ■ '"¿Г-
подставили это выражение в уравнение (1), раскрыв скобки, произведя сокращения подобных членов уравнения, получили:
Г К ак ^
I2
а
2 2 dt
3ёй '2
dt
- 43' - */(ёй2) ак-
■49 '2/( ^)[ а I
(К22 - К)
= о (3)
И
2
1
или
?к а2н
4/1 • *1
а^ + ^ (ёК)-/з /(ёК2)1Г
dh0
+4912 (^ )|dh2
9 ^ | аг ) gh1
Я\1(К - V)
ёК
* ' ,(V - V ) = 0.
+
0.
(4)
Разделив почленно уравнение (4) на (г2^) / (3ёК '2), получили:
а2к2 '2 *12 - 4 ёК '2 31*1 ак2
г к2 ёК2 з г к2 ёК2 аг
4 ёК \ зг2( ак2 ! 3ёК '2 *
9 г к2 ёК2 V аг) I2к2 ёК1 3ёК' (н22 - К)
г2К
= 0.
(5)
или
аг2
3ё • К' • к
2 кр
г2 •
( I
К М
V V
'2
1 +
2 У
г • к '22| 1 +
4Н ■
ас
аг
ас!2 3ё• к3кР
зн'
1 +
Н'
г2 • к
2 )
3ё •К'
К \2 +2К '2 •С + С2 -К!2
= 0.
г2 • к '
1 +
2)
Проверка
на
размерность
полученного выражения: т т • т • т3 т2 • т • т
2 + 2 2 2 2 + с с- т • т с • т • т- с
После проведения необходимых сокращений в уравнении (5), получили:
а Н2 3Н '2 2
аг2 г2 к
*1
4Н ан
*1
гк
аг
т • т
2 2 т • с
т • т
т • т • т
2 2 2 2 ' с • т • т с • т • т
которая показала, что размерность не нарушена. Тогда:
4 н \ ( ан2
3Н
з к22 V аг
г к2к1
2 *12
3ёК '2 (н22 - н12)
г2 к '
= 0.
(6)
В уравнении (6) приняли, что
К3кр = *12/ ё :
а2к2 3ёК'
н* -
4*1к '2 ак2 4 к '21 ак
аг2 г2К22 кр
г2к аг + з к2 V аг
3ёК '2 ,3 + 4 К \ ( аК2
" Н "т _
г2К2К1 кр з к22 V аг
3ёК -2 (н22 - н12) = 0 г2К22 2 .
3ёК '2
г2К
кр
Так как Н2 = Н '2 + ^ , то
2 2 а2 (н '2+с) 3ёК '2 н3кР
аг2 г2 (н '2+с)
4?1к '2 ^ (К '2+ С) +
г ( к\+с)2
4Н '„
3 (Н\+с)2
а (к\+с)! 3ё (к\+с) ^ +
аг
(к\+с)(г2к1)
3ё • к '2 ((н '2+С)- V)
г2 •( к '2+с)
= 0
ас
аг2
3ё • К
кр
4*1
'•11 + н'
•К '•|1 + н
а
2 аг
зн '«•|1 + н'
(
2 V аг,
3ё • К
кр
3ё
'2 •И + С
н \
н '22 + 2н ^2•C + C2 -н21
= 0.
(8)
В уравнении (8) приняли длину гидравлического прыжка равной трем значениям второй сопряженной глубины
(7) г = '2, тогда:
А?
3ё • К
кр
4?х
9н '22 • н '2 1 +
н'
зн '2 •! 1 +
н'
а
^ аг
'ас!2 зё • кз
зн \ 1 +
н'
зё
чаг) 9н '22 • н1
(н '22+2н '2 •с + с2 -н21)
9н1 +
н '
= 0.
(9)
Проверка размерности полученного уравнения (9):
+
2
+
т т • т3 т2 • т т2
2 2 2 2 2 с с •т •т с т • с т• с
т • т3 т • т2 т
- + ■
2 2 2 2 2 с •т •т с • т с
кр
dt2
9Л Л V | 1 +
Л'
4?1
3Л + Л"
2' dt
¿о2
кр
3Л '2„ • I 1
а2с+. л3кР
((
м2 3h '3
1 -
л
3Л'
С2
л
2 У
24h
2 у
з/ 4Л32
кр
3Л ,320
1-
3Л'
ас
2 У
Л:
кр
3Л 'V Л1 3Л ,22
24Л 2 dt
((г Г2 ^
1 - —+ —— Л '2 24Л '22 у
УЛ
( Л '2 л2
^ + Л V С + — - Л21 2 2 1
3Л '„
1 + ^
/у
2
3Л '2 24Л ,22 у
£1 =».
Раскрыв скобки, получили:
dí2 3Л,3
(
9Л '4
72Л!
4Л7
кр
4С • Л/2
2Ь/2
кр
С2л
л
кр
3Л'
,5/ 9Л /2„
18Л'
dí
4С
3Л '2 9Л ,22 18Л ,32
Л '2
Л V С С2
Л21
3Л ,22 Л1 +-2- + 6Л '2а 3Л '22 6Л ,22 6Л'
Л '2а С2 Л \ С3
18Л '32 9Л '32 18Л '32 ' 18Л '32 '
С2 л '22 , С3 Л \ С2л21
12Л '42 6Л '4а ' 12Л '42 12Л '42
Исходя из условий не возмущений в уравнении (12) исключили члены:
Л
кр
Л '
кр
Л '2
Л2
Уравнение (9) привели к безразмерному виду, разделив почленно на Л'2 и .-УЛ^^А^ :
Л г V dt У 3Л '„• Л,
3Л VI1+Л; , 2 1
1 ? г ^+Л ■• - £1=0- (10)
(12)
3' 2015
3Л '32 3Л 'V Л1 6 Л '2 6 Л '22
Одновременно решение полученного уравнения проводили методом малого параметра [5-9]. Приняли за малый параметр | = 1/ Л '2 [3], и уравнение (12) приняло следующий вид:
Л3
^ 4 dt2 ^
9С
кр ■ |
(Л!кр'
V 72 у
с2 -
-и
3-2 4Л\
dt
// 4Л32кр ¿С
и/2--кр • Г • — -
^ 9 dt
Л'
Избавились от рациональности, произведя необходимые разложения:
3/
-и.7'2 • ^+ 1. С2 -
18 ^^ 6
2 1 ^ 2 1 ^2 3 1 >-3 ,
-|и---д - и — <; -и---(, +
18 9 18
Л2 1 1
+и3' 1Т^ + ^? + И3-^«3 +
+и4' 1Т-и4- + 13- С = 0. (13)
После алгебраических преобразований получили:
,2 ^ ■
dt2
3
dt 3 V dt
-и2' (d5T+и2• С5 -
36
-И2 • ^ + = 0. 18 3
или
(11) ^
4С ( ¿с
dt2
9 V dt
-И' ^ • 1^1 +
?
-и
3/
3/ 4Л/2
2__
3 dt
36 18
(с+^ (2 ^
= 0.
(14)
Полученное уравнение описывает стационарные нелинейные колебания в нижнем бьефе.
Выводы
Полученное нелинейное уравнение описывает протекающий колебательный процесс в открытом водном потоке нижнего бьефа гидротехнического сооружения для случая сопряжения бьефов в виде гидравлического прыжка. Это выражение позволяет определять параметры нелинейных колебаний (амплитуда волны, частота колебаний,
(бб|
2
1
3
длина волны) возникающих в нижнем бьефе при установившемся режиме течения потока. Результаты вычислений применимы для коррекции размеров гасителей энергии водного потока при решении многих практических задач гидротехнических сооружений.
1. Грушевский М. С. Неустановившееся движение воды в реках и каналах. — Л.: Гидрометеоиздат, 1982. — 288 с.
2. Кюнж Ж. А., Холли Ф. М., Вервей А. Численные методы в задачах речной гидравлики. — М.: Энергоатомиздат, 1985. — 256 с.
3. Землянникова М. В. Фартуков В. А. Обобщенные нелинейные уравнения локальной нестационарности // Экологическая устойчивость природных систем и роль природообустройства в ее обеспечении: сб. материалов Всероссийской научно-технической конференции. — М.: МГУП, 200з. — С. 1зб-1з7.
4. Землянникова М. В. Фартуков В. А. Уравнения локальной нестационарности при прыжковых сопряжениях. / Сборник материалов Всероссийской научно-технической конференции «Экологическая устойчивость природных систем и роль природообустройства в ее обеспечении». — М.: МГУП, 200з. — С.1з7-1з8.
5. Найфэ А. Ю. Методы возмущений.
- М.: Мир, 1976. - 454 с.
6. Моисеев Н. Н. Асимптотические методы нелинейной механики. — М.: Наука, 1969. — 379 с.
7. Кузьмина Р. П. Асимптотические методы для обыкновенных дифференциальных уравнений. — М.: Едиториал УРСС, 2003. — 336 с.
8. Гукенхеймер Дж., Холмс Ф. Нелинейные колебания, динамические системы и бифуркации векторных полей. — Москва-Ижевск: Институт компьютерных исследований, 2002. — 560 с.
9. Боголюбов Н. Н., Митропольский Ю. А. Асимптотические методы в теории нелинейных колебаний. — М.: Государственное издательство технико-теоретической литературы, 1955. — 407 с.
Материал поступил в редакцию 16.06.2014. Фартуков Василий Александрович, кандидат технических наук, доцент E-mail: vasfar@mail.ru Тел.: 8 (916) 653-17-59 Землянникова Марина Владимировна, кандидат технических наук, профессор кафедры «Гидрологии, гидрогеологии и регулирования стока» Тел.: 8 (499) 976-22-27