Рисунок 1
Автор: Б. Ф. Лямаев
Источник: Политерм
История изучения нестационарных гидродинамических процессов, включающих в частности гидравлический удар, насчитывает уже более ста лет и начинается со ставшей в настоящее время уже классической работы Н.Е. Жуковского [1]. В прошедшие годы явление гидравлического удара неизменно привлекало внимание ученых и инженеров во многих странах в связи с научным и практическим аспектами. К началу шестидесятых годов прошлого столетия основные физические и гидродинамические процессы, протекающие в элементах трубопроводных систем при гидравлическом ударе, были в основном изучены. В то же время было установлено, что в реальных (сложных) гидросистемах специфика переходных процессов в основном определяется многократным наложением отраженных от конструктивных неоднородностей системы волн давления и трансформацией этих волн при их прохождении по длине трубопроводов. В сложных гидросистемах, как и во многих других инженерных задачах, решению препятствует «проклятие многомерности», усугубляемое нелинейностью процессов.
На переходные процессы оказывают влияние большое количество факторов, например, наличие в трубах растворенного в жидкости воздуха, деформационная податливость стенок труб и другие явления. В частности, было установлено, что на характер нестационарных процессов в трубопроводных системах кардинальное влияние оказывают кавитационные разрывы сплошности, которые могут возникать в любой точке системы при падении давления там ниже величины давления насыщенных паров жидкости. Поэтому не удивительно, что в конце шестидесятых годов прошлого века система, состоящая из двух труб разного диаметра, считалась у специалистов «сложной»и трудно поддающейся расчету.
Важным вкладом в методику расчета гидравлического удара было создание графического метода, изложенного лучше всего в монографии Л. Бержерона. Однако, даже специалистами в области гидравлического удара, расчеты могли выполняться лишь для достаточно простых гидросистем, состоящих из одной, двух труб, и без учета нелинейных и кавитационных процессов. Все это не позволяло гарантировать достаточную, даже для инженерных систем, точность расчета и возможность массового использования метода в повседневной практике.
Очередной (решающий) шаг в расчетах гидравлического удара - отказ от аналитических и графических методов в пользу численных, позволило сделать начало массового использования электронно-вычислительных машин. Их внедрение и массовое использование позволило проводить расчеты практически без ограничения сложности сетей.
Определенный итог выполненных в этом направлении работ и исследований был подведен в монографиях В. Л. Стритера и Е. Б. Уайли (США), Б. Ф. Лямаева, Г. П. Небольсина и В. А. Нелюбова (СССР), опубликованных в шестидесятых - семидесятых годах двадцатого века. В этих работах приведены методы и результаты решения основных проблем при создании компьютерных моделей для проведения численного эксперимента по изучению нестационарных процессов в реальных (сложных) трубопроводных системах на ЭВМ. Это позволило существенно повысить точность и достоверность расчетов, в том числе за счет привлечения более точных нелинейных моделей для описания распространения волн в трубах и поведения конструктивных элементов системы, учета кавитационных явлений.
Современным подходом при разработке программ для расчета режимов работы инженерных сетей, в частности и для расчета нестационарных гидродинамических процессов в трубопроводных системах, является создание программно-расчетных комплексов для компьютеров на основе использования геоинформационных технологий (ГИС). Это позволяет решать весь перечень задач, связанных с разработкой, проектированием и эксплуатацией гидросистем различного назначения для конкретных условий.
Применение геоинформационных технологий для создания моделей водопроводных сетей и, в частности, исследования нестационарных процессов дает существенное преимущество по сравнению с использованием других программных продуктов, не имеющих в своей основе геоинформационной системы. Это преимущество определяется следующим:
Возможностью использования графической подосновы (карты города, района, населенного пункта) для нанесения на нее водопроводной сети.
Простотой нанесения на карту города схемы водопроводной сети с ее привязкой к существующим зданиям и сооружениям, наглядностью представления информации, быстрым вводом исходных данных, необходимых для выполнения расчетов и удобством анализа полученных результатов расчета.
Процесс нанесение водопроводной сети на карту города максимально автоматизирован. Для этого достаточно при помощи мыши в графическом редакторе изобразить схему водопроводной сети. При этом автоматически создается топологически связный граф сети с привязкой соответствующих баз данных к каждому объекту.
В последнее время в связи с частыми отключеньями электроэнергии участились случаи возникновения гидроударов, поэтому более острой становится потребность в защите от гидроудара. При выборе защитных устройств существенную помощь может оказать компьютерная программа, позволяющая сразу оценить эффективность планируемых мероприятий.
Компания «Политерм», имеющая богатый опыт в разработке программного обеспечения для гидравлических и тепловых стационарных расчетов, предлагает новый программный продукт, предназначенный для расчета переходных процессов (в том числе и при возникновении гидравлического удара) в трубопроводных системах. Программный продукт создан с использованием опыта, накопленного научными школами одного из старейших вузов России Военного инженерно-технического университета. Программа, выполняющая расчеты гидравлического удара, оформлена в виде модуля расширения (plug-in) геоинформационной системы Zulu и входит в программно-расчетный комплекс ZuluHydro.
Кратко опишем возможности программы и используемую математическую модель. Известно, что переходные процессы в водопроводных сетях возникают чаще всего при включении и выключении насосов, а также при открытии и закрытии задвижек. Гидравлический удар является следствием быстрого выполнения одного или нескольких перечисленных действий и может усугубляться наличием воздушных полостей в трубах. В программе предусмотрена возможность выбора пользователем источника возмущения стационарного процесса. В процессе расчета и по его завершению пользователю для анализа переходных процессов программа предоставляет следующую информацию:
В процессе расчета пользователь имеет возможность наблюдать в реальном времени распространение бегущих волн давления и скорости вдоль любого маршрута.
По результатам расчета производится построение графиков наибольшего и наименьшего давлений в каждой точке вдоль этого маршрута.
Приводятся графики изменения давления в зависимости от времени для ряда выбранных точек наблюдения.
В базы данных заносятся значения наибольшего и наименьшего давлений для каждого участка и узла сети с указанием времени возникновения этих давлений, а для участка указываются и соответствующие точки.
ГИС Zulu позволяет с помощью запросов найти по базе данных наиболее опасные участки и узлы и на основании этой информации можно при необходимости выбрать (или добавить) другие точки наблюдения и изменить маршрут для построения графиков.
Отметим, что наряду с зонами повышенного давления определенную опасность представляют и зоны пониженного давления. В этом случае возможен подсос загрязненных грунтовых вод в систему водоснабжения, что может привести к возникновению эпидемий, кроме того, могут проявляться эффекты, связанные с кавитацией.
Остановимся подробнее на основных моментах описания математической модели, использованной при разработке программы расчета гидравлического удара.
Трубопроводная сеть (с математической точки зрения - взвешенный ориентированный граф) состоит из участков (труб) и узлов (соединение труб, потребитель, воздушный колпак, резервуар и так далее). С математической точки зрения насос и задвижка представляют собой участки, хотя на схеме они изображаются узлами. Пользователю о таких вещах беспокоиться не приходится - кодировка сети выполняется программой автоматически.
Распространение волн сжатия и разрежения в трубах описывается квазилинейной системой дифференциальных уравнений первого порядка с частными производными
В одном программном комплексе можно проводить расчеты, как стационарных режимов, так и нестационарных для любых по сложности систем водоснабжения.Приведем пример расчета переходных процессов для простейшей сети, состоящей из резервуара, колодца, задвижки и потребителя. На рисунках 1 и 2 показаны окно управления расчетами и окно результатов расчета.
Рисунок 1
Рисунок 1
На верхнем графике наибольшее и наименьшее давление в каждой точке вдоль трубы за все время численного эксперимента.
На среднем графике остановлена бегущая волна в последний момент расчета.
На нижнем графике показано изменение давления в зависимости от времени в трех выбранных точках наблюдения.
Вывод: Организацией «Политерм» создан современный программный комплекс для расчета стационарных и переходных процессов в трубопроводных системах любой сложности.
1. Жуковский Н. Е. О гидравлическом ударе в водопроводных трубах. М.-Л., Гостехиздат, 1949. 103 с.
2. Бержерон Л. От гидравлического удара в трубах до разряда в электрической сети. Общий графический метод расчета. (Перевод с франц.) М., Машгиз, 1962. 348 с.
3. Streeter V. L., Wylie E. B. Hydraulic Transients. New York, Mc Craw - Hill, 1967. 317 p.
4. Лямаев Б Ф., Небольсин Г. П., Нелюбов В. А. Стационарные и переходные процессы в сложных гидросистемах. Методы расчета на ЭВМ. Под ред. Б. Ф. Лямаева.- Л., Машиностроение. Ленинградское отделение, 1978. - 192 с.
5. Фокс Д. А. Гидравлический анализ неустановившегося течения в трубопроводах. М, 1981.