Автоматика. Экономика



бет1/5
Дата04.07.2016
өлшемі403.07 Kb.
  1   2   3   4   5





Раздел 5



Автоматика. Экономика




УДК (519.9+518.5):532.54




Н.И. КАРАСЕВ
Н.И. ТОМИЛОВА

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





Метод узловых давлений (МД), развитый в теории гидравлических цепей [1], обобщает метод узловых потенциалов из теории электрических цепей, сочетая преобразования Максвелла к узловым напорам с методом Ньютона, и базируется на процедуре линеаризации системы уравнений материального баланса в узлах с замыкающими соотношениями, описывающими стационарное изотермическое течение жидкости на пассивных и активных ветвях гидравлических цепей теплоснабжающих систем (ТСС). На основе исследований, выполненных в [1], сделан вывод о большем быстродействии и лучшей сходимости метода контурных расходов (МКР) по сравнению с МД, что на многие годы исключило попытки привлечения МД для конструирования автоматизированных решателей задач потокораспределения для информационных технологий, обеспечивающих разработку теплогидравлических режимов ТСС.

Однако в задачах машинного анализа сложных электрических цепей метод узловых потенциалов получил преимущественное развитие из-за большей вычислительной эффективности. Именно это обстоятельство побудило нас обратить внимание на другой подход к формированию системы уравнений узловых напоров для сложных гидравлических цепей, который опирается на инверсные характеристики пассивных и активных элементов трубопроводных сетей, замыкающих системы уравнений стационарного изотермического режима ТСС. Вычислительная эффективность этого подхода, как показали последующие исследования его программной реализации на реальных гидравлических цепях большой размерности [2], намного превосходят показатели эффективности метода МД, описанного в [1]. Модель стационарного гидравлического режима большой теплоснабжающей системы (ТСС) может быть выражена следующей нелинейной системой уравнений, записанной относительно узловых напоров исходя из закона о материальном балансе в линейно-независимых узлах гидравлической цепи



(1)

где hi, qi, Qjсоответственно, напор в i-м узле, расход на i-м участке, узловой расход в j-м узле;


aji, n, m — соответственно, элемент матрицы соединений, число участков и узлов в гидравлической цепи.

(2)

где R0i, R1i, R2i, n1, ri, μi — соответственно, параметры энергетической характеристики i-го насосного агрегата, число насосных агрегатов, гидравлическое сопротивление i-го участка и признак принадлежности участка к активной либо пассивной ветви цепи.

Если предположить существование решения векторного уравнения узловых напоров f(h) = 0 и заменить в нём h* = h0 + δh, где — некоторое начальное приближение к решению h*, то после разложения векторного уравнения f(h0 + δh) = 0 в двучленный ряд Тейлора по степеням малых приращений напоров в узлах δh получим векторное уравнение

(3)

покоординатная запись которого приводит к следующей системе линейных алгебраических уравнений относительно малых приращений напоров в линейно-независимых узлах:



(4)

где коэффициенты при δhi есть элементы матрицы Якоби для нелинейной вектор-функции f(h), вычисляемые в точке h0.

Решая линейную систему (4) относительно приращений или поправок δhj, , получим лишь их приближённые значения , с помощью которых можно улучшить исходные значения всех напоров в узлах. Представленные рассуждения привели к известному итерационному методу Ньютона — в данном случае для системы нелинейных уравнений узловых напоров. Элементы Якобиана для линеаризованной системы уравнений узловых напоров могут быть выражены посредством следующих соотношений:

(5)

(6)

(7)

Введем векторно-матричные обозначения для линеаризованной системы узловых напоров и получим компактные соотношения для диагональных и недиагональных элементов Якобиана



(8)

где — вектор поправок узловых напоров;


вектор правых частей;
— Якобиан вектор-функции узловых напоров.

Для пассивных элементов диагональный элемент Якобиана вычисляется из следующего соотношения:



а для активных ветвей

В самом общем случае диагональные элементы определяются следующим образом:

(9)

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

и (10)

Недиагональные элементы матрицы D могут получать числовые значения из следующих соотношений. Если узлы j и t связаны i-м пассивным элементом, то



Для тепловых сетей αi = 2 для а поэтому



(11)

Если узлы j и t связаны активным i-м элементом, то



(12)

Следующие свойства матрицы коэффициентов D в линеаризованной системе уравнений узловых напоров позволяют выбрать эффективные методы решения большеразмерных линейных систем:

- матрица D является симметрической, так как для неё справедливы отношения или djt = dtj. Настоящее свойство объективно обеспечено топологией технологической структуры теплоснабжающих систем, в которых каждая ветвь соединена парой узлов, а поэтому при табличном представлении структуры сети тождественность элементов, симметрично расположенных относительно главной диагонали, неизбежна;

- для всех узлов трубопроводной сети, непосредственно не связанных с опорным узлом, обеспечено строгое равенство диагональных элементов сумме (по модулю) остальных элементов в соответствующей строке матрицы D, как это следует из полученных соотношений (11)-(12);

- в реальных технологических схемах ТСС с одним независимым узлом могут быть связаны две-пять ветвей, а поэтому матрица D имеет много незаполненных элементов, т.е. относится к классу слабо заполненных матриц, а линеаризованная система DδH = F — к классу больших разреженных систем уравнений. Для ТСС любой технологической схемы заполненность матрицы D является своеобразной структурной константой, так как зависит только от числа узлов, ветвей и их соединений. Свойство слабой заполненности матрицы D определяет необходимость применения специальных методов нахождения корней линеаризованной системы DδH = F, исключающих операции с нулевыми элементами;

- система линеаризованных уравнений узловых напоров DδH = F хорошо обусловлена, так как приведенные выше свойства матрицы D для реальных сетей с физически правдоподобными значениями параметров элементов обеспечивают выполнимость условия существования и единственности решения detD ≠ 0.

Однако тщательные исследования алгоритмов прямых методов на базе гауссовских исключений с эффективными процедурами учета заполненности матрицы коэффициентов сделали эти методы более предпочтительными по сравнению с итеративными [3-4]. В этой связи однозначно решается вопрос выбора численного метода для решения большеразмерной системы линеаризованных уравнений узловых напоров — это эффективные алгоритмы метода исключения Гаусса с эффективными методами порядка исключения неизвестных. Реализация этих алгоритмов связана с использованием сложных схем хранения и доступа к данным в памяти реализующей вычислительной системы, а поэтому эффективность алгоритмов для разреженных матриц более чем в какой-либо другой области вычислений зависит от качества их машинной реализации.

Ньютоновский процесс сходится со скоростью квадратичной сходимости в малой окрестности корня h* + δh нелинейной системы уравнений узловых напоров (1) при условии Последнее условие есть не что иное, как условие хорошей обусловленности линеаризованной системы уравнений узловых напоров, которое выполнимо всегда при соответствующем выборе начального приближения.

Достаточно эффективным критерием завершения итерационного процесса Ньютона для систем уравнений общего вида является условие . Действительно, вблизи корня ньютоновские итерации сходятся квадратично, а поэтому, если этот критерий выполнен, то . Выбирая ε≈10–5–10–6, можно получить решение с не менее чем десятком значащих цифр в среде машинной арифметики двойной точности. Однако исходя из необходимости обеспечения равномерного приближения к решению по множеству узлов, достаточно эффективным условием выхода из ньютоновского итерационного процесса для систем нелинейных уравнений узловых напоров оказывается критерий

(13)

где Е1 — заданная допустимая максимальная невязка материального баланса по множеству линейно-независимых узлов моделируемой ТСС.

Удачное завершение ньютоновского итерационного процесса приводит к решению задачи определения напоров в узлах ТСС, т.е. к нахождению значений всех компонентов векторов (h*, Δh*, Р*).

Тогда решение задачи о потокораспределении на всех активных и пассивных ветвях ТСС сводится к вычислению компонентов вектора расходов с помощью следующих инверсных (расходно-напорных) характеристик:



(14)

СПИСОК ЛИТЕРАТУРЫ

1. Меренков А.П., Хасилев В.Я. Теория гидравлических цепей. М.: Наука, 1985. 277 с.

2. Карасев Н.И. и др. Моделирование установившихся гидравлических режимов больших теплоснабжающих систем с дроссельным управлением // Электронное моделирование. 1990. №1. С. 72-76.

3. Влах И., Сингхал К. Машинные методы анализа и проектирования электронных схем / Пер. с англ. М.: Радио и связь, 1988. 560 с.

4. Джордж А., Лю Дж. Численное решение больших разряженных систем уравнений / Пер. с англ. М.: Мир, 1984. 333 с.


УДК 621.0892




Т.Л. ТЕН


Каталог: wp-content -> uploads -> docs -> trudi%20univer
trudi%20univer -> Научные сообщения Әож 62-523=512. 122
trudi%20univer -> Пак ю. Н., Шильникова и. О., Пак д. Ю. Методологические аспекты организации самостоятельной образовательной деятельности студентов в контексте госо нового поколения
trudi%20univer -> Машиностроение. Металлургия Әож 621. 91. 02
trudi%20univer -> Проблемы высшей школы
trudi%20univer -> Машиностроение. Металлургия Әож 669. 779. 052: 553. 322 МҰхтар а. А
trudi%20univer -> Геотехнологии. Безопасность жизнедеятельности Әож 622. 271 СӘбденбекұлы ө
trudi%20univer -> Проблемы высшей школы
trudi%20univer -> Машиностроение. Металлургия Әож 621. 735. 34=512. 122 Ішкі беттерді өңдеуге арналған жайғыш бастиектерінің тозуға төзімділігін арттыру К. Т. Шеров
trudi%20univer -> Машиностроение. Металлургия
trudi%20univer -> Пак ю. Н., Нарбекова б. М., Пак д. Ю. Компетентностный подход в госо нового поколения и качество образования


Достарыңызбен бөлісу:
  1   2   3   4   5


©dereksiz.org 2019
әкімшілігінің қараңыз

    Басты бет