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

10. Аронович Г. В., Картвелишвили Н. А., Любимцев Я. К. Гидравлический удар и уравнительные резервуары. - М.: Наука, 1968. - 247 с.

11. Алышев В. М. Неустановившееся напорное движение реальной жидкости в трубопроводных системах: автореф. дис. ... д-ра техн. наук. - М.: МГМИ, 1987. -527 с.

12. Тарасевич В. В. О максимальном давлении при гидравлическом ударе, сопровождающемся разрывом сплошности потока // Гидротехническое строительство. - 1980. - № 8. - С. 15-18.

13. Рекомендации по расчету неустановившегося движения многофазной жидкости в напорных системах / Масс Е. И. [и др.] - М.: ВНИИ транспортного строительства, 1984. - 104 с.

14. Двухшерстов Г. И. Гидравлический удар в трубах некругового сечения и потоке жидкости между упругими стенками // Ученые записки МГУ. - 1948. - Т. 2. - С. 17-76.

15. Карамбиров С. Н. Совершенство-

вание методов расчета систем подачи и распределения воды в условиях многоре-жимности и неполной исходной информации: автореф. дис. ... д-ра техн. наук. -М.: МГУП, 2005. - 46 с.

16. Szumovski A. Pressure ave pattern in a liquid filling an elastic pipe // Arch. Mech., Stosow. - 1978. - Vol. 30. - №№ 4 and 5. - P. 645-656.

17. Егиазаров И. В., Картвелишвили Н. А., Первозванский А. А. К влиянию резинового шланга с воздухом при моделировании гидравлического удара // Известия АН СССР. Отделение технических наук. - 1957. - № 11. - С. 160-166.

18. Бегляров Д. С. Научное обоснование методов расчетов переходных процессов в напорных системах водоподачи с насосными станциями: автореф. дис. ... д-ра техн. наук. - М.: МГУП, 2007. - 59 с.

Материал поступил в редакцию 15.08.12. Картвелишвили Леонид Николаевич, доктор технических наук,главный научный сотрудник E-mail: l_n_k@mail.ru

УДК 502/504 : 532.5 Л. Н. КАРТВЕЛИШВИЛИ

Федеральное государственное бюджетное научное учреждение Центр научно-технической информации «Мелиоводинформ»

ПРИНЦИПЫ РАСЧЕТА ГИДРАВЛИЧЕСКОГО УДАРА И ИХ РАЗВИТИЕ

Дан краткий обзор современного состояния теории гидравлического удара. Отмече-но,что для консервативных (без учета потерь энергии) трубопроводов теория разработана детально. Выделены направления,касающиеся того,какие выражения принимать за исходные при составлении расчетных алгоритмов.

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

There is given a brief survey of the present state of the hydraulic impact theory. It is pointed out that for conservative (without consideration of energy losses) pipelines the theory is developed in detail. The directions are singled out concerning what expressions should be taken for the initial data when working out calculation algorithms.

Head pipeline, pressure pipeline, hydraulic impact, fluid flow velocity, losses of energy, system of pi pes, boundary conditions, princi ples of calculation, calculated algorithms.

Теоретически гидравлический удар рассматривали сначала как удар идеальной упругой жидкости в упругой трубе или в системе упругих труб, т. е. в гидравлических системах, без учета потерь энергии. Это консервативные системы. Исследования, в процессе которых определяли только скорость распространения ударных волн по трубам, начали проводить в конце XVIII века. Как указывается в [1], еще в 1775 году Л. Эйлер анализировал решения волновых уравнений применительно к системам кровообращения. В конце XIX столетия появилась фундаментальная работа Н. Е. Жуковского, в которой исследован весь процесс гидравлического удара в консервативной системе [2].

Уравнения гидравлического удара. Вывод дифференциальных уравнений гидравлического удара выполняется двумя путями. Гидравлический путь основан на следующих допущениях [2]:

инерция поперечных перемещений жидкости и оболочки трубопровода не учитывается;

неустановившееся движение жидкости в трубе рассматривается как одномерное;

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

Применение второго закона механики к массе жидкости, заключенной между двумя указанными сечениями, приводит к уравнению движения:

где g - ускорение свободного падения; и - средняя по сечению скорость потока; £ - время; х - расстояние по оси трубы; у = z + p/(pg) - пьезометрический напор; z - отметка оси трубы в сечении с абсциссой х; р - давление в этом сечении на оси; р - плотность жидкости; а1 - корректив количества движения; Ь - гидравлический уклон.

Рассматривая баланс вещества в отсеке между теми же сечениями, можно получить уравнение неразрывности для жидкости, подчиняющейся закону Гука: ад | № | ^ Ар | Ам> _ Ах Аг Е0 с!г й '

где Q = Fu - расход жидкости; F - площадь сечения трубы; Е0 - ее объемный адиабатический модуль упругости; d - диаметр трубы; w - расход жидкости через единицу площади внутренней поверхности стенок трубы (для проницаемых стенок).

Уравнения (1) и (2) должны быть дополнены реологическим уравнением для материала трубы и выражениями для Ь и w. Реологическое уравнение для «упругой» трубы при справедливости гипотезы Жуковского есть

где К - коэффициент упругого отпора стенок, т. е. давление, вызывающее приращение внутреннего радиуса на единицу длины.

Формулы расчета К для труб различных конструкций из упругих материалов приводятся в работе [3].

Из (2) и (3) в случае непроницаемых стенок, т. е. при w = 0 можно получить уравнение

йи

¿У.+

= О,

(4)

& g сЬс

где с - скорость распространения упругих возмущений по трубопроводу, т. е. скорость распростра нения ударных волн; С

С =

(5)

(2)

КА

где с0 - скорость распространения звука в жидкости.

Основные дифференциальные уравнения гидравлического удара (1) и (4) справедливы и для труб некруглого сечения, но скорость распространения ударных волн в таких трубах определяется по другим формулам [2]. При этом, например, небольшая эллиптичность трубы (отношение длин осей эллипса 0,9...0,95), как указано в [4], практически не влияет на эту скорость, но при отношении длин осей около 0,6 скорость распространения ударных волн снижается в несколько раз. Вообще значительные отклонения сечения трубы от формы круга приводят к значительному снижению скорости ударных волн [5].

Второй путь вывода уравнений гидравлического удара - из уравнений гидродинамики, уравнений теории оболочек и условий контакта жидкости и оболочки - является более строгим. Таким путем получены уравнения только для консервативных систем с тонкостенными трубами [4].

Отношение скоростного напора a1u2/(2g), входящего в уравнение (1), к изменениям пьезометрического напора при гидравлическом ударе, имеет порядок и/с, как было показано еще Жуковским. Величина и/с очень мала, так как значения

с в среднем около 1000 м/с, а и обычно не превышает 7 м/с. В связи с этим скоростным напором в теории гидравлического удара пренебрегают. Тогда уравнение (1) принимает следующий вид:

+ ^ +1 = 0. (6)

g а ах

Для идеальной жидкости г = 0. В этом случае, дифференцируя уравнение (4) по t, а уравнение (6) по х и исключая переменную и, а затем дифференцируя уравнение (4) по х , а (6) по t и исключая переменную у, можно получить выражения

2 2 2 2 у 2 д у д и 2 о и (7)

= с

= с

д дх д дх

совпадающие с уравнениями колебаний струны.

Существует также другая аналогия, которой посвящена монография [6]. Если считать, что гидравлический уклон i линейно зависит от скорости ^ = ти) и что имеется соответствие: у - напряжение; и

- сила тока; m,1/g, g/c2 - соответственно омическое сопротивление, самоиндукция и емкость на единицу длины проводника

- то уравнения (4) и (6) совпадают с уравнениями течения электрического тока по линейному проводнику. Если труба пористая и Ру - утечка жидкости на единицу длины, а проводимостью изоляции проводника О нельзя пренебречь, то Р соответствует О. Значение данной аналогии уравнений ограничено из-за различия граничных условий, от которых зависит решение - эти условия совпадают лишь для узлов гидравлических и электрических сетей и не совпадают для устройств на концах ветвей.

Получив уравнения (7) первым из описанных путей, Жуковский указал их общий интеграл:

Ау=у-у=ф--)+

с с

= (8)

с с с

где ф и 1 - произвольные функции, зависящие от краевых условий; у0 и и0 - начальные значения у и и.

Используя выражения (8), он сделал физический анализ гидравлического удара, получил широко известную и широко применяемую формулу прямого удара - Ду = сД u/g, рассмотрел приложения этих результатов. Во всех последующих исследованиях гидравлического удара в

консервативных системах использовались выражения (8) или выражения, полученные из них путем тождественных преобразований.

Гидравлический удар в простых трубопроводах. Л. Аллиеви исследовал гидравлический удар в простом трубопроводе [6, 7]. (Простым называют трубопровод, имеющий постоянный диаметр и постоянную скорость распространения ударных волн, начинающийся от большого резервуара, давление во входном сечении которого можно считать постоянным, и заканчивающийся запорным устройством, через которое жидкость вытекает в атмосферу). В работе [6] Аллиеви употреблял выражения (8), а в [7], используя граничное условие в начале трубопровода у — уо = 0, выразил функцию 1 через ф и, исключив ф из (8), получил уравнение в конечных разностях:

(9)

имеющее название цепного уравнения Аллиеви. Индексы указывают моменты времени, к которым относятся значения у и и на расстоянии х от начала трубопровода, где т = 2х/с - двойное время пробега волной удара расстояния х. Если х = I -длина трубопровода, то т называют фазой гидравлического удара. Закон истечения, связывающий расход запорного устройства или скорость в конце трубопровода с напором, под которым происходит истечение, Аллиеви принимал в следующей форме:

и

И,

= а

У

У.

(10)

где а-А/4, - гидравлическое открытие, характеризующее положение запорного устройства; А -некоторая функция открытия запорного устройства, зависящая от его конструкции; А0 - значение А при полностью открытом запорном устройстве; и* и у* -значения и и у при установившемся режиме с полностью открытым запорным устройством (а = 1).

В работе также принято, что гидравлическое открытие изменяется от начального значения а0 до конечного а к за некоторое время Т по линейному закону. Из выражений (9) при х = I и (10) получается нелинейное уравнение в конечных разностях:

У- «У- 1 У■ \У■

Аллиеви, анализируя данное уравнение, определил, что при линейном изменении а величина \у — у\ достигает максимального значения в конце первой фазы (£ = |) или в конце движения запорного устройства (в последнем случае удар называется предельным: давление в трубе постепенно приближается к некоторому предельному значению).

Используя теорию Аллиеви, М. Га-риель рассмотрел вопрос о наибольшем возможном гидравлическом ударе и установил, что если полностью (ак = 0) закрывать запорное устройство от различных начальных открытий а0, во всех случаях с одной и той же скоростью 1/Т , то у — у, будет иметь максимальное значение не при а0 =1, а при а0 = |/Тв, т. е. при таком начальном открытии, когда время закрытия равно фазе удара [8].

Для этого случая Гариель получил 2 1и,

у* ' у- = иг- (12)

Такой удар называется коротким ударом.

Решения для сложных трубопроводов. Следующим этапом в развитии теории гидравлического удара был переход к трубопроводам, которые характеризуются сложной интерференцией ударных волн в точках разветвления (узлах): волна удара, подойдя к узлу по какой-либо ветви трубопровода, частично проходит в другие ветви - преломляется, а частично отражается обратно в ту же ветвь. (Частным случаем сложного трубопровода является телескопический трубопровод, у которого к каждому узлу примыкают только две ветви). В работе «О гидравлическом ударе в водопроводных трубах» Н. Е. Жуковский рассмотрел сложные трубопроводы с небольшим количеством ветвей при мгновенном закрытии запорных устройств. Решения для более сложных случаев, найденные из анализа многократного отражения и преломления ударных волн (представленные, например, Ч. Егером в 1933 году), громоздки и не совсем удобны для использования даже с применением компьютера. Значительно более удобным является применение уравнений Шниде-ра-Бержерона, получающихся путем исключения произвольных функций ф и f из уравнений (8):

У?

с , - - —(и

8 '

Vм -

2

- «,"), (13)

У? = -Ми; 2 8* 2 где М и N - сечения; индексы - моменты времени, к которым относятся значения у и и; |1 - двойное время пробега волной удара расстояния между сечениями М и N.

Сложение уравнений (13) и замена в них £ на £ - |/2 дает [9]:

К +

- 2у\ = - <). (14)

-I 8

В ряде случаев выражением (14) пользоваться проще, чем уравнениями Шнидера-Бержерона, из которых оно получено. Если сечение N совпадает с началом простого трубопровода, то = уо

2

и (14) обращается в цепное уравнение Ал-лиеви (9).

На основе выражений (13) Р. Леви (1928), О. Шнидер и Л. Бержерон разработали способы графического расчета гидравлического удара в сложных трубопроводах [10]. Уравнения (13) описывают сразу суммарное действие всех волн, учтенных в функциях ф и поэтому при их использовании не надо рассматривать отдельные волны, их отражение и преломление в местах разветвления трубопровода. Но сегодня и эти методы не имеют широкого применения в связи с неудобством их реализации на компьютере.

В 1948 и 1952 годах были опубликованы два алгоритма, не имеющие подобного недостатка [11, 12]. Первый из них заключается в следующем. В (14) £ заменяется сначала на t — \У, затем на I — 2ц, на ^ — Зц и так далее до тех пор, пока все индексы в этом выражении не станут такими, при которых соответствующие у и и будут относиться к установившемуся режиму, т. е. иметь значения у0 и и. Затем все строки суммируются, и после перехода к безразмерным величинам получается:

IIм = им -

где и =

Я

У-у.

), (15)

Р /-1 ¡=о

, = Т*

Я, н. т 8рн. 1 *

- одинаковые для всей исследуемой системы (а в остальном произвольные) расход, напор и время соответственно.

Если считать М за точку разветвления, а N1,..,Nn - за противоположные концы отходящих от М ветвей, то из (15) и уравнения неразрывности

+ ... + я?» = 0 (16)

получается:

фМ +5. (У, Iм -УФ*1

X 1 4 Лш^ Т-1У1 ^ VI .

¡=0 2

1 I

) = 0,

м

Т-1УЯ

(17)

где 5 =—/У— (к = 1, ..., п) - коэффициенты

1 Р>/ Р. преломления ударных волн в точке М.

Величина Т* выбирается так, чтобы для всей системы каждое А было целым числом. Получающаяся система уравнений, дополненная граничными условиями, связывающими Ф и и в указанных вершинах, позволяет последовательно определить Ф во всех вершинах графа сети для целых значений т.

Второй метод [12] основан на записи уравнений (13) в следующей форме:

= о:,

где Я =

1

+ и,

гг

п

= п\

1

о,

р ' р

Из (19) следует: 2

■А - и.

+ П = — Р

а - п = 2и.

(18)

(19)

(20)

После подстановки в уравнение неразрывности (16) значения и из (19) получаем:

фм =

у у Пм

^ X Лш^ X

(21)

Пользуясь этим выражением, гидравлический удар в точке разветвления следует выражать через функции П, если скорости течения считаются положительными по направлению от точки разветвления, и через функции О, если скорости течения считаются положительными по направлению к точке разветвления. В 1981 году Д. Авакумовичем был предложен алгоритм, основанный на операторном преобразовании уравнения (14).

При проведении расчетов для сложных трубопроводов необходимо, чтобы Ф были целыми числами. Для этого в качестве Т* принимают наибольший общий делитель фаз т всех ветвей системы, причем Т* может оказаться очень малым числом. В подобных случаях изменяют длины ветвей системы так, чтобы увеличить Т* и соответственно уменьшить максимальное значение Ф. Описанные решения получены в предположении, что гидравлические сопротивления по длине отсутствуют ^ = 0).

Граничные условия, создаваемые запорными устройствами на трубопро-

водах. Граничными условиями для уравнений гидравлического удара обычно называют условия, характеризующие запорные устройства на трубопроводах. Например, условие (10) соответствует обычным запорным устройствам (задвижки, краны и др.) при истечении жидкости в атмосферу, а, например, для реактивных турбин связь расхода с напором сложнее и может быть выражена полуэмпирической формулой [9]:

и

= а

У-

(23)

где к - постоянная, зависящая от типа турбины (а для поворотно-лопастных турбин - и от угла установки лопастей); п, п0 - соответственно фактическая и синхронная частота вращения турбины.

Для активных турбин к = 0 и (23) переходит в (10). Для насосов и насосотур-бин формул, аналогичных (23), нет, и связь между расходом, напором, частотой вращения и открытием регулирующих устройств дается в виде характеристик насосов и на-сосотурбин. Часто и для турбин используют не формулу (23), а их универсальные характеристики.

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

Правильное задание краевых условий очень значимо. В [4] показано, что процессы гидравлического удара при граничном условии типа (10)

У.

(24)

( у1 и у2 - пьезометрические напоры, соответственно выше и ниже по течению относительно запорного устройства) и при упрощенном условии и/и* = а, часто встречающемся в теоретических работах, с одинаковой функцией а(1) протекают совершенно по-разному.

В то же время, как показал Д. Спарр (1915), условие (24) можно упростить, приняв У\~ У^ = У' ДУ , разлагая с учетом этого правую часть (24) в ряд по степеням Ду и отбрасывая высшие степени Ду. Это приводит к замене нелинейной зависимости (24) линейной:

1 + -Ах.

и

и.

= а

2 у.

(25)

использование которой дает ошибку менее 1% в максимальной величине удара, если тах |Ду| < У./2.

Выводы

Таким образом, состояние теории гидравлического удара характеризуется следующими главными положениями.

Для консервативных (без учета потерь энергии) трубопроводов теория разработана детально. Существует общий интеграл дифференциальных уравнений гидравлического удара. Представленный в виде уравнений Шнидера-Бержерона

(13) или обобщенного уравнения Аллиеви

(14), он применяется для расчетов удара в простых трубопроводах. На основе этих уравнений получены алгоритмы расчета для сложных трубопроводных сетей.

Для трубопроводных систем с диссипацией (с учетом потерь) энергии общего способа расчета гидравлического удара нет.

Среди различных направлений развития теории гидравлического удара можно выделить два направления, касающиеся того, какие выражения принимать за исходные при составлении расчетных алгоритмов. Многие исследования основаны на решении дифференциальных уравнений гидравлического удара как без учета сопротивлений, так и для линейного или линеаризованного законов сопротивлений при определенных граничных условиях. Такой подход позволяет проводить детальные исследования и расчеты, но не всегда удобен из-за недостаточной гибкости. Другое направление связано с использованием в качестве исходных не дифференциальных уравнений удара, а их общего интеграла в виде уравнений Шнидера-Бержерона (13), уравнения Аллиеви или обобщенного уравнения Аллиеви (14). Как уже отмечалось, на основе этих зависимостей получены удобные для реализации на компьютере методы расчета гидравлического удара в сложных трубопроводах (без учета сопротивлений по длине) [11, 12]. М. Маршал, Г. Флеш и П. Сьютер (1968) использовали решения Аллиеви для разработки программы расчета гидравлического удара в трубопроводах сложных гидроэнергетических систем. В [4] дано обобщение (14) для малых линеаризованных гидравлических сопротивлений. Сделаны попытки обобщить уравнение Аллиеви на случай сложных

№ 4' 2012

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

1. Paynter H. M. Fluid Transients in Engineering Systems: Handbook of Fluid Dynamics; V.L. Streeter, Ed. - New York: McGraw-Hill, 1961. - Chapter 20.

2. Жуковский Н. Е. О гидравлическом ударе в водопроводных трубах: полн. собр. соч. - M. - Л.: Гл. ред. авиац. литературы, 1937. - Т. 7. - С. 5S-157.

3. Мостков М. А. Гидравлический справочник. - M.-Л.: Госстройиздат, 1954. -532 с.

4. Картвелишвили Н. А. Динамика напорных трубопроводов. - M.: Энергия, 1979. - 224 с.

5. Daily J. W. and Deemer K. The Unsteady-Flow Water Tunnel at the Massachusetts Institute of Technology // Trans. Asme. - 1954. - № 1. - P. S1-SS.

6. Allievi L. Teoria generale del moto perturbato dell' acqua nei tubi in pressio-ne. - Annali della Societa' degli Ingeneri ed Architetti, 1903.

7. Allievi L. Teoria del colpo d'ariete. -Milano, 1913. - 210 p.

S. Gariel M. Etude sur les maxima de surpression dans les phe'nomenes de coup de b'elier // Revue generale de l' electricite.

- 191S. - № 23. - P. 403-410.

9. Картвелишвили Н. А. Неустановившиеся режимы в силовых узлах гидроэлектрических станций. - M.-Л.: Госэнер-гоиздат, 1951. — 136 с.

10. Бержерон Л. От гидравлического удара в трубопроводах до разряда в электрической сети (общий графический метод расчета); пер. с франц. (с приложением И. А. Чарного и Г. Д. Розенберга). -M.: Mашгиз, 1962. - 34S с.

11. Картвелишвили Н. А. Расчет гидравлического удара в сложных системах // Гидротехническое строительство.

- 194S. - № 3. - С. 15-20.

12. Мостков М. А., Башкиров А. А. Расчеты гидравлического удара. - M. -Л.: Госэнергоиздат, 1952. - 24S с.

Ыатериал поступил в редакцию 15.0Б.12. Картвелишвили Леонид Николаевич, доктор технических наук. главный научный сотрудник E-mail: l_n_k@mail.ru

II"