Б01: https://doi.org/10.33408/2519-237X.2020.4-3.328
УДК 532.59+627.8
ОБ ИНТЕГРИРОВАНИИ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ НЕУСТАНОВИВШЕГОСЯ ПОСТЕПЕННО ИЗМЕНЯЮЩЕГОСЯ ДВИЖЕНИЯ ПОТОКА В ОТКРЫТОМ РУСЛЕ В УСЛОВИЯХ ВЫСОКОГОРЬЯ ПРИ ПРОРЫВЕ ПЛОТИНЫ
Стриганова М.Ю., Шаталов И.М., Самедов С.А., Щербакова М.К., Недашковская И.В., Рабченя В.С.
Цель. Выбор научно-технически обоснованного метода интегрирования уравнений неустановившегося постепенно изменяющегося движения потока в открытом русле в условиях высокогорья при прорыве плотины.
Методы. Методы приближенного интегрирования дифференциальных уравнений.
Результаты. Предложено использование конечно-разностного метода интегрирования дифференциальных уравнений (метода характеристических уравнений С.А. Хри-стиановича) для неустановившегося постепенно или плавно изменяющегося движения в условиях высокогорья при прорыве плотины.
Область применения исследований. Предложенное решение уравнений неустановившегося постепенно изменяющегося движения в условиях высокогорья при прорыве плотин можно использовать при решении практических задач по определению основных параметров волны перемещения и зон вредного воздействия на прилегающую к высокогорью территорию.
Ключевые слова: плотина, высокогорье, неустановившееся движение, волна перемещения, параметры волны, скорость потока (волны), глубина потока (волны), время перемещения, длина распространения.
(Поступила в редакцию 2 июня 2020 г.)
Введение
Решение задачи о неустановившемся движении потока жидкости в открытом русле, как правило, заключается в интегрировании системы двух уравнений: уравнения баланса расхода (т. е. уравнения неразрывности постепенно или плавно изменяющегося неустановившегося движения потока жидкости в открытом русле) и уравнения динамического равновесия (т. е. дифференциального уравнения постепенно или плавно изменяющегося неустановившегося движения потока жидкости в открытом русле).
В работе [1] эта система уравнений представлена в следующем виде:
8 б 8 ю л —- +-= 0,
81 8г (1)
О - Аб2 )g
е 8 ю 8 V 8 V
•10- '
а 0--ь а V-
В 81 0 8 г 81
где б - объемный расход воды, м3/с; ш - площадь живого сечения, м2; V - средняя скорость, м/с; I - длина рассматриваемого участка, м; г - момент времени, с; г - уклон свободной поверхности воды; А - удельное сопротивление русла, с2/м2; А б2 - уклон трения /тр; В - ширина русла по поверхности потока, м; а,0 - коэффициент Буссинеска; а - коэффициент Кориолиса.
Учитывая, что б = га, и введя обозначение (/ - Аб2)е = Е, перепишем систему уравнений в следующем виде:
д ю дv д ю
-+ ю — + v-= 0,
д t dl dl (2)
„ g дю дv дv
E =---+ a 0--+ a v—.
B д1 0 дt дl
В результате решения системы уравнений определяются две основные функции: Q = fi(l,t) и о = f2(l,t). Зная эти функции, в практических расчетах можно определить основные параметры (или характеристики) потока: среднюю скорость - v = fi(l,t) и глубину - h = f4(l,t) в любом створе потока и далее построить мгновенный профиль свободной поверхности потока или волны перемещения.
Приведенная система уравнений является системой нелинейных дифференциальных уравнений гиперболического типа, обладающих двумя совокупностями характеристик.
Интегрирование системы уравнений (2) в общем случае представляет достаточно большие трудности, поэтому в инженерной практике широкое применение нашли методы приближенного интегрирования этих уравнений.
Рассмотрим один из методов решения уравнений (2) с учетом возможности использования компьютерных технологий.
Так как система уравнений (2) относится к классу уравнений гиперболического типа с двумя совокупностями характеристик, эту систему уравнений можно заменить уравнениями соответствующих им характеристик, которые могут быть решены методами приближенного интегрирования. Исследование системы уравнений неустановившегося движения потока в открытом русле и решения их методом характеристических уравнений впервые было проведено С.А. Христиановичем [2].
И.М. Бернадский, В. А. Архангельский и другие разработали методику решения дифференциальных уравнений неустановившегося движения жидкости в открытых руслах способом конечных разностей для определенных значений временных отрезков (мгновений при t = const). Такой метод в расчетной практике получил название метода мгновенных режимов, хотя более уместно его назвать методом мгновенных фрагментов (или ситуаций).
Основная часть
В реальных условиях высокогорья при прорыве плотины происходит относительно постепенное опорожнение водохранилища, при котором наблюдается падение уровня воды в водохранилище, уменьшение расхода воды в начальном створе и увеличение расхода в конечном сечении прямой отрицательной волны перемещения [1].
Движение воды в теле такой волны перемещения хорошо описывается двумя дифференциальными уравнениями баланса расхода и уравнением динамического равновесия [2].
Для численного (с использованием компьютерных технологий) решения системы уравнений (2) для условий высокогорья при прорыве плотины предлагается использовать метод характеристических уравнений С. А. Христиановича [2], при этом движение воды в теле волны перемещения будем считать постепенно или плавно изменяющимся.
Рассмотрим более подробно метод характеристических уравнений С. А. Христиановича, дополнив его моделированием волн на мелкой воде методом частиц [3].
На схеме движения потока (рис. 1) позиция 1 соответствует свободной поверхности в момент времени ti, а позиция 2 - в момент времени t2. Рассмотрим неустановившееся движение потока воды между сечениями I-I и II-II на бесконечно малом расстоянии А/.
Представим неустановившееся постепенно или плавно изменяющееся движение в виде параллельных прямолинейных отрезков линий тока (см. рис. 1). Допустим, что в области решения системы уравнений (2) задан отрезок некоторой кривой линии тока функцией / = /(t) и значения функций v = v(/,t) и ш = ш (/,t).
1 - начальное положение поверхности волны перемещения в момент времени 2 - конечное положение поверхности волны перемещения в момент времени t + Рисунок 1. - Отрезок некоторой кривой линии тока
Тогда для каждой точки этой кривой (точки представляют собой частицы движущейся жидкости) / = /(^ можно записать:
й ю dt
д V й/ д V
д/ ^ д t '
д ю й/ д ю
д / ^ д t
или
д V
д ю
~д7
dv dt й ю dt
д V й/ д/ dt' д ю й/ д/ dt
(3)
Подставив значения и д- в систему (2), после простейших преобразований бу-
дем иметь:
д V
ю--ъ
д /
V - ■
ё/ Л д ю
dt) д /
й ю dt
й/ Лд V я д ю
„„ й — „ dv
а V - а 0 — I--ъ —--= Е - а 0 —.
dt ) д /В д / 0 dt
(4)
Откуда легко получить
dt dt
/
а V - а
а V - а
V
dt)
ач_Л
dt
gю В
Яю В
дv
дю "д/
Е - а,
dv dt
/
/
а V - а
V
й/ ^ йю
dt
+
Я йю В dt
dt
dt
ю
dvЛ Е - а 0 —
<# у
(5а)
В условиях высокогорья в водотоках всегда наблюдается быстрое течение воды, называемое бурным. Бурный поток движется в развитом турбулентном режиме, для которого коэффициенты Кориолиса и Буссинеска можно принимать равными 1,0. В этом случае систему уравнений (9) можно записать в более компактном виде, удобном для анализа и последующего решения.
ш л2 dt у
Ш ^2 dt у
go B
go B
дv дю
эТ
dv
E--
dt
d[ dt
+
g Шю B dt '
dt_ dt
Шю dt
/
ю
dv E--
v dt j
(5)
Системы уравнений (5а) и (5) можно применять для расчета неустановившегося постепенно или плавно изменяющегося движения в открытых руслах произвольной формы поперечного сечения. Такие русла называются непризматическими, и для них должно соблюдаться условие ш = _Д1; И), где ш - площадь живого сечения потока; I - длина потока; И -глубина потока.
Однако в условиях высокогорья при растекании бурного потока его глубина значительно меньше его ширины В (т. е. В>>И) и поперечное или живое сечение такого потока близко к прямоугольной форме. Прямоугольное сечение, равно как и любое другое сечение правильной формы (например, трапецеидальное, треугольное, овальное, параболическое и т. д.) относится к призматическим руслам, для которых ш = _ДИ) и дojд1 = 0 . Учитывая это для условий высокогорья при прорыве плотины, можем записать:
df
av - a0 —
0 dt
v —
d__
dt
Шю dt
Шю dt
- ю
- ю
dv
E - a0 —
0 dt
0
dv
E--
dt
0.
(6a)
(6)
С физической точки зрения производная dl|dt является скоростью распространения фронта (тела) волны перемещения, причем эта волна перемещения имеет два направления (прямое и обратное), для определения которых можно записать следующие уравнения, приведенные в [2]:
(dljdt)1 = v + yjgo / B и (dljdt)2 = v — ^go / B .
(7)
Подставив выражения (7) в уравнение (6) получим формулы для определения скоростей распространения фронта (тела) прямой и обратной волны перемещения:
dv = —.
g
Bo
d o + Edt и dv =+.
g
Bo
d о + Edt
В обозначениях
(7) перепишется в виде
yjg о / B = a1 и v — yjg о / B =
(dljdt)1 = a1 и (dljdt)2 = a2 .
(8)
(9)
(10)
Для любого элементарного участка Л1 постепенного или плавно изменяющегося движения в конечных разностях систему уравнений (10) запишем в виде:
A/j = a1At1, А/2 = a2 At2.
(11)
Уравнения системы (11) для мгновенных режимов или фрагментов времени tl, t2, ... являются уравнениями прямых линий на элементарном участке кривой А/ с постоянными угловыми коэффициентами а1 и а2, т. к. в любой зафиксированный момент времени общее неустановившееся движение является установившимся, при котором значения б, ш, V являются постоянными величинами.
При интегрировании в конечных разностях в координатной плоскости (/, ^ для любого заданного отрезка кривой / = /(^ можно записать расчетные уравнения (11) в следующем виде:
где /а, ta, /ь, и - координаты точек (частиц жидкости) в пределах кривой / = /({), в которых известны (или заданы) функции (или значения) V и ш; /т, ^ - неизвестные координаты точек (частиц жидкости) кривой / = представляющей собой линию тока неустановившегося движения потока воды.
Представим уравнения (8) определения скоростей распространения фронта (тела) волны перемещения (прямой и обратной) в конечных разностях:
где Va, vь и ша, шь - начальные скорости распространения и площади живого сечения волны перемещения; Vm и шт - текущая скорость распространения и площадь живого сечения волны перемещения (прямой и обратной).
Для решения практических задач и компьютерного моделирования неустановившегося постепенно или плавно изменяющегося движения воды в условиях высокогорья в виде волны перемещения прямой или обратной, положительной или отрицательной наиболее применим метод конечных приращений [2]. Подобный метод был использован Томпсоном для расчета прямоугольных русел, который с некоторыми дополнениями и изменениями можно распространить на русла произвольной формы поперечного сечения.
Рассмотрим русло произвольной формы поперечного сечения (см. рис. 1). Разделим это русло на элементарные участки А/, в пределах которых площадь живого сечения Аш будет изменяться постепенно (или плавно). Рассмотрим конкретный элементарный участок, в начальном сечении которого, как и в последующих сечениях, параметры неустановившегося потока (глубины И, скорости V, площади живых сечений ш и т. д.) известны в данный момент времени t и в последующие отрезки времени Аt. Предположим, что на рисунке 1 линия 1 соответствует положению свободной поверхности волны перемещения в начальный момент времени ^ а линия 2 - это положение свободной поверхности той же волны по истечении отрезка времени Аt, т. е. в момент времени t+Аt.
Определим средние значения параметров неустановившегося потока в любом его сечении для отрезка времени А1:
/ - / = а1(t - t ),
т а 1V т а/'
/ - = а2(t - tь),
т Ь 2\ т Ь
(12)
(13)
ю
= (ю1 + ю 2 + ю 3 + ю 4) /4,
(14)
В = (В1 + В2 + Вз + В4 )/4, Л = ( + Я2 + Яз + Я4)/ 4, V = ( + у2 + У3 + У4 )/ 4,
С = (С1 + С 2 + С 3 + С 4) / 4,
где ш - площадь живого сечения, м2; V - средняя скорость, м/с; В - ширина русла по поверхности потока, м; Я - гидравлический радиус, м; С - коэффициент Шези, м0,5/с. В системе уравнений (1) неустановившегося постепенно или плавно изменяющегося движения уклон трения /тр = Л02 на элементарном участке потока А/ можно выразить из урав-
72 - дк
v
дю
нения Шези /тр = —2— . С учетом того что для призматического русла — = 0 и / = /0 ——,
С Я д/ д/
где /о - уклон дна водотока, уравнение движения системы (1) представимо в виде:
г2 дк а0 ду ау ду
v
lo=^= +
+
+
C2R dl g dt g dl
v2
lo=J^ +
dh 1 dv v dv +--+
С2 Я д/ я дг g д/ •
Частные производные в конечных приращениях представим в виде:
(15а) (15)
dh 1 ( h 3 - h1 3 1 i h 4 - h 2 1 = h i+h 2 - h 3 - h 4
dl 2 t Al ' Al J 2 Al '
dv 1 ( v 3-v 1 + v 4 - v 2 '1 = v 1 + v 2 - v 3 - v 4
dl 2 t Al Al J 2 Al '
dQ d( wv) = 1 ( w3v 3-w1v 1 + W4v 4 Ш2v 2 1 =
dl dl 2 t Al A l J
dh 1 (h2-h1 + h 4 - h 3 1 = h 1 - h 2+h 3 - h 4
dt 2 t A t ' At J 2 At '
dv 1 ( v 2-v 1 + v 4 - v 3 1 = v 1 - v 2 + v 3 - v 4
dt 2 t A t At J 2 At '
dw 1 (ш2-w 1 + Ш4-Ш31= ш 1 -ш2+w3-w4
d t 2 t A t At J 2 A t
Qxv 1+Ш2v 2-Ш3У 3-Ш4v 4
2А l
(16)
где к;, к2, кз, к4 - глубина потока в рассматриваемых сечениях за отрезок времени Аг; У1, У2, уз, У4 - средние скорости за отрезок времени Аг; Ш1 ,Ш2 ,шз ,Ш4 - площади живых сечений за отрезок времени Аг.
Согласно уравнениям (16) уравнения (15а) и (15) в конечных разностях примут вид:
v2
l =__h 1 +h 2-h 3-h 4 -а v 1-v 2 +v 3-v 4 avv 1 +v 2-v 3-v 4
If) - ^ry --Ыа \J*V
0 C2R 2А1 0 2 gAt 2 g Al
(17а)
. _ _ h 1 +h 2-h 3-h 4 _ У i-V 2 +У 3-V 4 _ 1 +v 2 -v 3 -v 4
0 C2R 2 А/ 2 g At 2 g Al ' (17)
Уравнение баланса расхода (или уравнение неразрывности) постепенно или плавно изменяющегося неустановившегося движения потока жидкости в открытом русле системы (1) в конечных приращениях согласно (16) принимает вид
c1v 1 +c2v 2-c3v3-c4v4 со1 -с2 +с3-с4 _ о
2А/ 2At " ' (18)
Уравнения (17а), (17) и (18) позволяют найти параметры h и v неустановившегося потока в любой отрезок времени At и в любых сечениях этого потока, а также построить кривую свободной поверхности волны перемещения (прямой и обратной) в условиях высокогорья при прорыве плотины.
Заключение
1. Представленный конечно-разностный метод интегрирования дифференциальных уравнений неустановившегося постепенно или плавно изменяющегося движения для определенных отрезков времени t = const (метод мгновенных режимов или фрагментов) является достаточно приближенным. Однако этот метод наиболее полно отвечает требованиям реальной инженерной практики и позволяет осуществить компьютерное моделирование процесса распространения волны перемещения (как прямой, так и обратной) в условиях высокогорья при прорыве плотины.
2. Предложенный метод интегрирования дифференциальных уравнений неустановившегося плавно изменяющегося движения по способу конечных разностей будет проверен экспериментально в лабораторных и натурных условиях.
ЛИТЕРАТУРА
1. Стриганова, М.Ю. Математическая модель пространственно изменяющегося неустановившегося движения потока при прорыве напорных и гидротехнических сооружений в условиях высокогорья / М.Ю. Стриганова [и др.] / Вестник Университета гражданской защиты МЧС Беларуси. - 2020. - Т. 4, № 1. - С. 48-58. DOI: 10.33408/2519-237X.2020.4-1.48.
2. Богомолов, А.И. Гидравлика / А.И. Богомолов, А.И. Михайлов; 2-е изд., перераб. и доп. - М.: Стройиздат, 1972. - 648 с.
3. Богомолов, С.В. Моделирование волн на мелкой воде методом частиц / С.В. Богомолов, Е.В. Захаров, С.В. Зеркаль // Математическое моделирование. - 2002. - Т. 14, № 3. - С. 103-116.
Об интегрировании дифференциальных уравнений неустановившегося постепенно изменяющегося движения потока в открытом русле в условиях высокогорья при прорыве плотины
On the integration of differential equations of unsteady gradually changing flow in an open channel in highland conditions when a dam breaks
Стриганова Марина Юрьевна
кандидат технических наук, доцент
Государственное учреждение образования «Университет гражданской защиты Министерства по чрезвычайным ситуациям Республики Беларусь», кафедра автоматических систем безопасности, доцент
Адрес: ул. Машиностроителей, 25, 220118, г. Минск, Беларусь e-mail: striganovam@tut.by ORCID: 0000-0002-8100-733X
Шаталов Игорь Михайлович
Белорусский национальный технический университет, кафедра гидротехнического и энергетического строительства, водного транспорта и гидравлики, старший преподаватель
Адрес: пр-т Независимости, 65,
220013, г. Минск, Беларусь e-mail: shatl19@yandex.by ORCID: 0000-0002-5348-5318
Самедов Самедага Абзар оглы Министерство по чрезвычайным ситуациям Республики Азербайджан, главное оперативное управление, заместитель начальника Адрес: ул. М. Мушвига, 501,
AZ1073, г. Баку, Азербайджан e-mail: ssamedaga@yandex.by ORCID: 0000-0002-4241-8080
Щербакова Мария Константиновна
Белорусский национальный технический университет, кафедра гидротехнического и энергетического строительства, водного транспорта и гидравлики, старший преподаватель
Адрес: пр-т Независимости, 65,
220013, г. Минск, Беларусь e-mail: mk_shcherb@tut.by ORCID: 0000-0002-8864-0517
Marina Yu. Striganova
PhD in Technical Sciences, Associate Professor
State Educational Establishment «University of Civil Protection of the Ministry for Emergency Situations of the Republic of Belarus», Chair of Automatic System Security, Associate Professor
Address: ul. Mashinostroiteley, 25,
220118, Minsk, Belarus e-mail: striganovam@tut.by ORCID: 0000-0002-8100-733X
Igor'M. Shatalov
Belarusian National Technical University, Chair of Hydrotechnical and Power Engineering, Water Transport and Hydraulics, Senior Lecturer
Address: pr-t Nezavisimosti, 65, 220013, Minsk, Belarus e-mail: shatl19@yandex.by ORCID: 0000-0002-5348-5318
Abzar S. Samedov
Ministry of Emergency Situations of the Republic of Azerbaijan, Main Operational Department, Deputy Head
Address: ul. M. Mushviga, 501,
AZ1073, Baku, Azerbaijan e-mail: ssamedaga@yandex.by ORCID: 0000-0002-4241-8080
Mariya K. Shcherbakova
Belarusian National Technical University, Chair of Hydrotechnical and Power Engineering, Water Transport and Hydraulics, Senior Lecturer
Address: pr-t Nezavisimosti, 65, 220013, Minsk, Belarus e-mail: mk_shcherb@tut.by ORCID: 0000-0002-8864-0517
Недашковская Ирина Васильевна
Белорусский национальный технический университет, управление подготовки научных кадров высшей квалификации, методист Адрес: пр-т Независимости, 65,
220013, г. Минск, Беларусь e-mail: nedash19@yandex.by ORCID: 0000-0002-7059-9122
Рабченя Виктория Сергеевна
Белорусский национальный технический университет, факультет энергетического строительства, студент Адрес: пр-т Независимости, 65,
220013, г. Минск, Беларусь e-mail: rabash19@yandex.by ORCID: 0000-0002-9350-1693
Irina V. Nedashkovskaya
Belarusian National Technical University, Department for the Training of Scientific Personnel of Higher Qualification, Methodologist
Address: pr-t Nezavisimosti, 65, 220013, Minsk, Belarus e-mail: nedash19@yandex.by ORCID: 0000-0002-7059-9122
Viktoriya S. Rabchenya
Belarusian National Technical University, Department of Energy Construction, Student
Address: pr-t Nezavisimosti, 65, 220013, Minsk, Belarus e-mail: rabash19@yandex.by ORCID: 0000-0002-9350-1693
DOI: https://doi.org/10.33408/2519-237X.2020.4-3.328
ON THE INTEGRATION OF DIFFERENTIAL EQUATIONS OF UNSTEADY
GRADUALLY CHANGING FLOW IN AN OPEN CHANNEL IN HIGHLAND CONDITIONS WHEN A DAM BREAKS
Striganova M.Yu., Shatalov I.M., Samedov S.A., Shcherbakova M.K., Nedashkovskaya I.V., Rabchenya V.S.
Purpose. The choice of a scientifically and technically substantiated method for integrating the equations of unsteady gradually changing flow movement in an open channel in highland conditions when a dam breaks.
Methods. Approximate integration methods for differential equations.
Findings. The use of the finite-difference method for integrating differential equations (method of characteristic equations of S.A. Khristianovich) is proposed for unsteady gradually changing flow in highland conditions when a dam breaks.
Application field of research. The proposed solution of the equations of unsteady gradually changing movement in highland conditions when breaking dams can be used to solve practical problems in determining the main parameters of the displacement wave and the zones of harmful effects on the territory adjacent to the high mountains.
Keywords: dam, highland, unsteady movement, displacement wave, wave parameters, flow velocity (waves), flow depth (waves), displacement time, propagation length.
(The date of submitting: June 2, 2020) REFERENCES
1. Striganova M.Yu., Shatalov I.M., Samedov S.A., Nedashkovskaya I.V., Rabchenya V.S. Matematich-eskaya model' prostranstvenno izmenyayushchegosya neustanovivshegosya dvizheniya potoka pri proryve napornykh i gidrotekhnicheskikh sooruzheniy v usloviyakh vysokogor'ya [Mathematical model of a spatially variable unstable flow motion at the breakthrough of hydrotechnical structures under the conditions of highland]. Journal of Civil Protection, 2020. Vol. 4. No. 1. Pp. 48-58. (rus). DOI: 10.33408/2519-237X.2020.4-1.48.
2. Bogomolov A.I., Mikhaylov K.A. Gidravlika [Hydraulics]. Moscow: Stroyizdat, 1972. 648 p. (rus)
3. Bogomolov S.V., Zakharov E.V., Zerkal' S.V. Modelirovanie voln na melkoy vode metodom chastits [The shallow water wave simulating by particle method]. Matematicheskoe modelirovanie, 2002. Vol. 14. No. 3. Pp. 103-116. (rus)