CC BY
39
УДК 622.45
A.B. Шалимов, 8912-48-58-977, shalimovav@mail.ru (Россия, Пермь, Горный институт УрО РАН)
АДАПТАЦИЯ МЕТОДА УЗЛОВЫХ ДАВЛЕНИЙ К РАСЧЁТАМ ВОЗДУХОРАСИРЕДЕЛЕНИЯ В РУДНИЧНЫХ ВЕНТИЛЯЦИОННЫХ СЕТЯХ
Проведён аналитический обзор численных методов моделирования воздухорас-пределения в вентиляционных сетях. Исследованы на сходимость в различных условиях алгоритмы методов последовательных приближении, контурных расходов и узловых давлений. Установлены причины плохой сходимости метода узловых давлений, на основании анализа которых произведена корректировка метода и устранены выявленные в нём недостатки.
Ключевые слова: воздухораспределение, аэродинамические сопротивления, естественная тяга, тепловые депрессии, напорные характеристики, невязка давлений, сходимость.
Современные рудники и шахты представляют собой сложные разветвлённые многоуровневые сети, состоящие из большого количества горизонтальных, вертикальных и наклонных выработок. Проветривание таких вентиляционных сетей осуществляется, как правило, несколькими источниками тяги, - наземными и подземными вентиляторами главного проветривания и вспомогательными вентиляторами местного проветривания. Определённый вклад в движение воздуха вносят также естественная тяга и тепловые депрессии, возникающие в негоризонтальных выработках при значительной разнице температур воздуха в них. Потребность в моделировании проветривания вентиляционных сетей подобной сложности вызывает необходимость в разработке специализированных математических методов расчёта, позволяющих с достаточной точностью и за реальное время определять воздухораспределение в различных условиях.
Вентиляционная сеть задаётся ориентированным графом, аэродинамическими сопротивлениями всех ветвей и напорами всех источников тяги. Математическая постановка задачи расчёта воздухораспределения заключается в составлении уравнений 1-го и 2-го законов сетей. Пусть сеть состоит из п ветвей и V узлов. Первый закон сетей записывается в виде у-1 уравнений в узлах:
где - номер узла (от 1 до V-1); у - номера выработок, инцидентных узлу с номером s, Q(s) - расходы воздуха.
Поскольку общее число уравнений должно быть равно числу неизвестных п, то оставшееся число уравнений ш=п-^-1). Это уравнения, отражающие второй закон сетей:
(1)
^(±д*Ч))=Е(±У' оЧ) , (2)
3 3
где ' - номер (от 1 до т) уравнения (независимого контура); у - номера выработок, содержащихся в контуре с номером I; Аж(} - напор источника тяги
в выработке с номером у контура с номером '. Схема сети задаётся в виде ориентированного графа, на котором все выработки имеют направление предполагаемого движения воздуха (ветви со стрелкой). Стрелки расставляются без нарушения первого закона сетей во всех узлах, т.е. нет узлов, в которые стрелки только входили бы или только выходили. В подсистеме уравнений (1) для определённости можно принять, что знаки перед расходами положительны, если стрелки направлены в узел, и отрицательны, если направлены они из узла. В подсистеме уравнений (2) знаки зависят от направления обхода контура: по часовой или против часовой стрелки (выбор направления может быть любой). Если направление напора источника тяги Ал(,г) совпадает с направлением обхода, то ставится знак “+”, если нет
- то “-“ (знак же самого напора, так же как и знак расхода, определяется совпадением с направлением стрелки выработки). Если направление стрелки выработки совпадает с направлением обхода, то перед Я(') ставится знак “+”, если наоборот, то “-“. Знак расхода Q(1) в (2) от направления
обхода не зависит и определяется только направлением стрелки выработки: “+”- по стрелке, “-“- против стрелки.
Система уравнений (1) - (2) является нелинейной, и для её решения, вообще говоря, могут быть использованы стандартные математические методы решения нелинейных уравнений, например, метод Ньютона. Однако, как показывает практика решения вентиляционных задач, метод сходится только при достаточно «удачном» выборе начального приближения, подбор которого для систем большой размерности является весьма длительной и трудоёмкой процедурой. Следует заметить, что метод Ньютона предназначен для решения уравнений и систем уравнений произвольной нелинейности, а система (1) - (2) таковой не является. Она состоит из двух подсистем, одна из которых (1) является линейной, а вторая (2) - квадратичной, т.е. слабо нелинейной. Это обстоятельство послужило поводом к разработке других, более эффективных специализированных методов расчёта систем уравнений такого типа. Первым из них был метод последовательных приближений (МПП), разработанный в первой половине прошлого века. Суть метода заключается в последовательной поконтурной увязке расходов с постепенным уменьшением невязки давлений в контурах. Метод достаточно прост в реализации, изначально был предназначен для ручных вычислений (электронно-вычислительной техники в то время ещё не было) и имеет широкую область сходимости по начальным приближениям в отличие от метода Ньютона. С появлением и развитием ЭВМ метод совершенствовался на предмет ускорения сходимости (путём увязки расхо-
дов сразу по нескольким контурам) и на сегодняшний день является наиболее широко используемым при решении сетевых вентиляционных задач.
С наибольшей точностью и математической корректностью процедура уменьшения невязки давлений в контурах была реализована в методе контурных расходов (МКР), в котором увязка расходов производится сразу во всех независимых по расходам ветвях сети. Метод достаточно сложен в реализации, занимает значительно больше оперативной памяти, чем МПП, сходится на порядок быстрее и может по праву считаться наилучшим вариантом метода последовательных приближений. К сожалению, ввиду того, что необходимое для программной реализации метода развитие вычислительной техники было достигнуто не так давно, МКР для решения вентиляционных задач распространения пока не получил.
В работах описан ещё один специализированный метод решения системы уравнений (1) - (2) того же уровня сложности и математической корректности, что и МКР, но принципиально от него отличающийся. Это метод межузловых депрессий (ММД) ([4]) (другое название - метод узловых давлений). Суть метода заключается в уменьшении невязки расходов воздуха в узлах путём увязки падений давлений во всех независимых по давлениям выработках (либо увязки давлений во всех узлах). Таким образом, ММД симметричен МКР относительно порядка решения подсистем (1) и (2). Преимуществом ММД перед МКР является тот факт, что уменьшаемая невязка расходов в узлах позволяет говорить о точности их вычисления, в то время как результат решения уравнений с помощью МКР позволяет оценивать лишь точность вычисления давлений. По результатам проведённых численных экспериментов установлено, что ММД имеет гораздо худшую сходимость, чем МКР. Один метод (МКР) имеет отличную сходимость, другой (ММД) - практически не сходится. Выяснить причину расходимости ММД и найти возможность её устранения можно на основе сравнения и анализа математических алгоритмов обоих методов.
Алгоритм метода контурных расходов имеет следующую структуру.
1. Задаются начальные приближения расходов воздуха Qk в хордах графа (независимых по расходам выработках сети).
2. Из решения линейной подсистемы уравнений (1) находятся значения остальных расходов воздуха в ветвях дерева графа (зависимых по расходам выработках сети).
3. Подстановкой значений расходов воздуха в правую часть подсистемы уравнений (2) находятся невязки давлений во всех независимых контурах.
4. Поскольку падения давлений во всех выработках могут быть представлены как функции независимых расходов, то увязка этих расходов может быть связана с невязкой давлений в независимых контурах с помощью стандартной процедуры разложения функции многих переменных в
ряд Тейлора с сохранением слагаемых только первого порядка (по аналогии с методом Ньютона):
¿р, = £
і ,к
дО 4
о,
дО,
' АОк . (3)
-7
где і - номер независимого контура от 1 до т;, - номер ветви от 1 до п; к -номер хорды 1 до т; дР, - невязка давлений в і-м контуре;; - расход воз-
духа в,-й ветви; Ок - расход воздуха в к-ой хорде; Ал,= Дя^і) - напорная характеристика источника тяги в ,-й ветви (если нет, то 0); Я. - сопротивление і'-й ветви в і-м контуре (-Я,, 0 или +Я, в зависимости от наличия ветви в контуре и её ориентации относительно обхода контура). Производные
О К Яр.)
находятся из решения подсистемы уравнений (1), производные
определяются аналитически или численно в зависимости от слож-
10,
ности напорных характеристик. Из решения системы уравнений (3) определяются значения увязочных расходов А0к.
5. Определяется очередное приближение расходов воздуха в хордах Окг+1) = Ок + А0к (г -номер итерации), и процесс повторяется с пункта 2 расчёт прекращается, когда значения невязок давлений в контурах становятся меньше заданного значения.
Применительно к ММД система уравнений (1) - (2) записывается в
виде:
0 = £(±р<->) , (4)
Р 1в1в1 \ = Р, + д*, (в,) , (5)
0 = Х (±в *)) (6)
У
где * - номер узла; ' - номер независимого контура;, - индекс суммирования по номерам ветвей, инцидентных узлу с номером * в подсистеме (6), и по номерам ветвей входящих в контур с номером ' в подсистеме (4); Р, -
падение давления вместе с напором источника тяги Дл| в ветви с номером, (с индексом * - для ветвей, инцидентных узлу *, с индексом ' - для ветвей
контура '). Знаки выбираются так же как и в системе (1) - (2). Алгоритм
метода узловых давлений имеет следующий вид.
1. Задаются начальные приближения падений давлений Р* в ветвях дерева графа (независимых по давлениям выработках сети). Либо задаются давления во всех узлах графа, а Р* находятся.
2. Из решения линейной подсистемы уравнений (4) находятся значения остальных падений давления на хордах графа (зависимых по давлениям выработках сети).
3. По формуле (5) определяются все расходы воздуха.
4. Подстановкой значений найденных расходов в правую часть подсистему уравнений (6) находятся невязки расходов во всех узлах.
5. Поскольку расходы воздуха во всех выработках могут быть представлены как функции независимых падений давлений, то увязка этих падений может быть связана с невязкой расходов в узлах по аналогии с
где * - номер узла от 1 до у;, - номер ветви от 1 до п; ? - номер ветви дерева 1 до у; §в* - невязка расходов в *-ом узле; Р, - падение давления в ,-й ветви (с напором источника тяги); Рг - падение давления в ¿-й ветви дерева (с напором источника тяги); Я*, - сопротивление ,-й ветви (-Я,, да или +Я, в зависимости от инцидентности ветви Б-му узлу и её направленности - в узел (+) или из узла (-)). Если Я*,<0, то знак «-» выносится из под корня. дР,
Производные —- находятся из решения подсистемы уравнений (4). Из
дР(
решения системы уравнений (7) определяются значения увязочных давлений ДРг.
6. Определяется очередное приближение перепадов давлений в ветвях дерева Р/г+1) = Р/ + ДР? (г -номер итерации), и процесс повторяется с пункта 2 расчёт прекращается, когда значения невязок расходов в узлах становятся меньше заданного значения.
Алгоритмы МПП, МКР и ММД были реализованы численно и проанализированы на предмет сходимости при расчёте воздухораспределения для вентиляционных сетей реальных рудников. Результат - МПП сходится медленно, МКР сходится на порядок быстрее и ММД не сходится вообще. Оказалось, что причина плохой сходимости ММД заключается вовсе не в преобладании недиагональных элементов в узловой матрице (Максвелла), как это представлено в [2], а в большой разнице аэродинамических сопротивлений выработок реальных рудников (от 10-10 кмюрг для сбоек до 1000кмюрг для перемычек). Численный эксперимент показал, что ММД начинает сходиться при разнице сопротивлений не более трёх порядков, а при разнице в два порядка сходится не хуже МКР. Причина плохой сходимости ММД состоит в том, что при большой разнице сопротивлений двух соседних ветвей на одной из них с большим сопротивлением падает весь напор, а на другой ветви с меньшим сопротивлением падение напора практически равно нулю. Ноль оказывается в знаменателе выражения, что и не даёт итерационному процессу ММД сойтись, т.к. он осциллирует с огромным шагом вокруг этих «нулевых» решений. Выравнивание же сопротивлений «уничтожает» нулевые решения и позволяет ММД сойтись. Исходя
МКР:
(7)
из этого, можно сделать вывод, что метод узловых давлении в классическом виде не пригоден для расчёта воздухораспределения в реальных вентиляционных сетях, поскольку область сходимости метода в три порядка на разброс сопротивлений выработок слитком ограничена.
Несмотря на указанные недостатки, ММД может быть модернизирован на предмет улучшения сходимости следующим образом. Чтобы избежать появления нулевых величин в знаменателе, в качестве независимого базиса следует взять не перепады давлений, а расходы воздуха в ветвях
дерева, т.е. не Pj, a Qj =
І
Pl
R,
и раскладывать невязки расходов в узлах по
AQj, а не по ДРу. В результате этих преобразований изменяются пункты 5 и 6 алгоритма ММД (пункты 1 - 4 остаются прежними):
7. Поскольку расходы воздуха во всех выработках могут быть представлены как функции расходов в ветвях дерева, то увязка этих расходов может быть связана с невязкой расходов в узлах:
SQs ~ÇS AQ'
(8)
где £ - номер узла от 1 до у; і - номера инцидентных узлу б ветвей; ґ - номер ветви дерева 1 до V; дQs - невязка расходов в £-м узле; Q¡ - расход в і-
дР
ой ветви; Qt - расход в (-й ветви дерева. Производные
дР
находятся из
решения подсистемы уравнений (4), после чего определяются производные
dQt
2 Rt\Qt
d (A^) 3Qt
2 R.
Ql
d {Аж.)
dQ
dP
Производные
^(A^-) dQ
определяются анали-
j y
тически или численно в зависимости от сложности напорных характеристик источников тяги Ащ = Акі(Оіі) . Из решения системы уравнений (8) определяются значения увязочных расходов ДQt.
8. Определяется очередное приближение перепадов давлений в ветвях дерева Р{(г+1) = Я; (Q(( г) + AQt) г) + AQt - Аж{ (Qt(г) + Ад{) (г -номер
итерации), и процесс повторяется с пункта 2. Расчёт прекращается, когда значения невязок расходов в узлах становятся меньше заданного значения. Проведённая корректировка алгоритма себя полностью оправдала -
ММД.
A.V. Shalimov
ADAPTATION JUNCTION PRESSURES METHOD TO CALCULATING AIR-DISTRIBUTION IN MINING VENTILATION NETWORKS
Analytical survey of modeling air-distribution in ventilation networks numerical methods was made. Algorithm convergence at different conditions of successive approximations method and methods of outline flows and junction pressures were researched. Reasons of bad junction pressures method survey were ascertained and correcting this method was made and ascertained shortcomings were eliminated.
Key words: air-distribution, aerodynamic resistance, natural draft, thermal barometric depression, pressure characteristic, pressures residual, survey.
Получено 20.04. 11
