Распределенное имитационное моделирование системы теплоснабжения зданий в среде VisSim

УДК 517.958:532.5; 697.1

РАСПРЕДЕЛЕННОЕ ИМИТАЦИОННОЕ МОДЕЛИРОВАНИЕ СИСТЕМЫ ТЕПЛОСНАБЖЕНИЯ ЗДАНИЙ В СРЕДЕ VISSIM

А.А. Басалаев

В статье рассматривается подход к построению распределенной имитационной модели централизованной системы теплоснабжения зданий с использованием ЭВМ в среде имитационного моделирования VisSim. Представленный подход основан на разделении модели на взаимосвязанные подсистемы, моделирование которых может производиться на узлах вычислительной сети одновременно. Предложены способы синхронизации подсистем, полученных в результате разделения исходной модели. Описаны преимущества и недостатки предложенных способов синхронизации. Представленные результаты моделирования систем теплоснабжения с различным числом потребителей подтверждают, что скорость расчета модели с большим количеством элементарных вычислительных блоков и нелинейных уравнений увеличивается при разделении ее на взаимосвязанные подсистемы.

Ключевые слова: распределенные вычисления, имитационное моделирование, система централизованного теплоснабжения, VisSш.

Введение

Сегодня для управления большими сложными системами все чаще используются имитационные модели. Такие модели зачастую представляют собой сложные системы нелинейных уравнений с большим количеством неизвестных, для нахождения которых требуется выполнение большого объема вычислений. При этом вычислительных мощностей процессора и оперативной памяти отдельных ЭВМ уже недостаточно для обеспечения процесса расчета. Поэтому для решения таких систем становятся все более актуальными разработка и применение методов распределенных и параллельных вычислений. При этом развиваются преимущественно два направления распределенного имитационного моделирования:

1. Разработка решателей систем нелинейных уравнений, позволяющих распределить сам процесс решения единой сложной системы уравнений между ЭВМ. При развитии этого направления наряду с математическими методами активно применяются методы параллельного программирования.

2. Разработка методов модификации моделей, направленная на изменение структуры и параметров модели таким образом, чтобы ее было возможно рассчитывать параллельно на нескольких ЭВМ или ядрах процессора.

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

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

Метод разделения на подсистемы не зависит от решателя, используемого той или иной средой моделирования. Поскольку подсистемы моделируются в отдельных процессах, то распарал-

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

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

Принцип построения распределенной имитационной модели

Модель теплоснабжения зданий может быть достаточно большой или отдельные ее блоки могут быть достаточно сложными, но в общем случае в системе есть источники, соединительные элементы (трубопроводы, трубопроводные тройники и т. д.) и потребители. Рассмотрим принцип разбиения модели, построенной, например, в среде имитационного моделирования VisSim согласно методике, описанной в [4, 5] (рис. 1).

| Источник И Труба 7 ^

Исходная модель системы

гНІРУбаТ}

Делитель потока 1

точка разделения -

-^ГПГуба6 |

Делитель потока 2 Труба 1 и Здание 1

Труба 2 Здание 2

Делитель потока 3 Труба 3 Здание 3

Труба 4 Здание 4

| Источник!-^ Труба 7 ^

Делитель потока 1

Модель подсистемы №1 гНЗйбЦ

Делитель Труба 1 Здание 1

потока 2 Труба 2 Здание 2

Ос Вход на трубу 6

Модель подсистемы №2

узел коммуникации с моделью №2

узел коммуникации с

моделью №1 — | Выход с делителя потоков 1

Труба 6 [

Делитель Труба 3 Здание 3

потока 3 Труба 4 Здание 4

Рис. 1. Структурное разделение исходной модели на подсистемы

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

Связи между блоками в исходной модели

Модель №1

Модель №2

Рис. 2. Функциональное разделение исходной модели на подсистемы

Изначально имеем информационные связи между частями модели из 6 переменных. Это температуры, давления и расходы в подающем и обратном трубопроводах.

Для первой модели узел коммуникации представляет собой потребителя. На узел коммуникации как на потребителя подается теплоноситель с определенной температурой T1 и известным перепадом давлений P1 и P2. Эти данные передаются узлом коммуникации первой модели в узел коммуникации второй модели, например, через ЛВС с использованием стека протокола TCP/IP. Во второй модели узел коммуникации является источником с уже известными характеристиками (принятыми из первой модели P1, P2 и T1), на основании которых во второй модели производится расчет расходов теплоносителя G1 и G2 и температуры теплоносителя в возвращающем трубопроводе T2, что и характеризует вторую модель как потребителя. В свою очередь, рассчитанные параметры передаются узлом коммуникации второй модели в узел коммуникации первой. И на основе информации о второй модели как о потребителе происходит пересчет параметров тепловой сети в первой модели.

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

Математически система нелинейных уравнений была разделена на 2 части. При этом модифицировались все уравнения, в которых участвовали 6 переменных исходной связи. Остальные уравнения были разделены между частями модели.

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

Синхронизация процессов моделирования подсистем

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

Отметим преимущества данного подхода:

- простота реализации (только настройка параметров проекта в среде моделирования);

- нет задержек при моделировании частей системы и системы в целом, связанных с ожиданием данных;

- время передачи данных не влияет на время расчета модели.

Недостатки этого подхода заключаются в следующем:

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

- необходимо подбирать временной запас шага моделирования из-за возможного увеличения длительности расчета отдельного шага при переходных процессах;

- расчет на каждом шаге моделей с малой реальной длительностью шага растягивается до величины глобального шага.

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

Преимущества данного подхода - это то, что:

- каждая модель рассчитывается в собственном модельном времени;

- ввиду этого при наличии синхронизатора нет необходимости расчета и ввода временного запаса;

- процесс синхронизации не зависит от длительности расчета шагов частей модели.

Недостатки этого подхода заключаются в том, что:

- модель усложняется дополнительными блоками синхронизации;

- время расчета модели зависит от времени передачи данных;

- части модели с малой длительностью шага моделирования простаивают в ожидании такта синхронизации.

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

Результаты моделирования

Описанный подход был исследован на базе модели системы теплоснабжения ФГБОУ ВПО «Южно-Уральский государственный университет» (НИУ) [5]. В ходе экспериментального моделирования на ЭВМ с использованием 2 ядер с частотой 2,3 ГГц процессора AMD А10-4600М при разбиении статической модели, построенной в среде VisSim, на 2 подсистемы с применением подсистемы синхронизации были получены результаты, представленные в таблице и на графике (рис. 3). Модуль синхронизации в модели подсистемы № 1 состоит из 400 элементарных вычислительных блоков. Модуль синхронизации в модели подсистемы № 2 - из 120 блоков. В таблице число блоков в подсистемах указано без учета блоков модулей синхронизации. В ходе каждого эксперимента было произведено 10000 шагов моделирования. Для увеличения количества блоков в экспериментах система дополнялась однотипными группами объектов, состоящих из потребителей, трубопроводов и трубопроводных тройников.

Результаты моделирования

Число функциональных блоков в исходной системе Время выполнения исходной системы, с Число блоков в подсистеме № 1 Число блоков в подсистеме № 2 Время выполнения при разделении, с

2500 14,7 1300 1200 83,4

5000 58,7 2600 2400 98,5

7500 207,2 3800 3700 133,4

10000 583,0 5000 5000 143,7

12500 1229,8 6200 6300 314,4

Согласно графику (рис. 3), время расчета исходной модели с увеличением количества блоков возрастает быстрее, чем время расчета модели, разделенной на подсистемы. Это связано не только с тем, что каждая подсистема рассчитывается в отдельном процессе, но и с тем, что в каждой подсистеме меньше и элементарных вычислительных блоков, и решаемых на каждом шаге моделирования уравнений.

Рис. 3. Результаты моделирования

Заметим, что для исходных систем с количеством блоков 2500 и 5000 разбиение на подсистемы нецелесообразно, поскольку в этих случаях на каждом шаге моделирования время обмена данными между подсистемами больше времени расчета каждой подсистемы.

Заключение

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

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

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

Литература

1. Mukai, H. Parallel algorithms for solving systems of nonlinear equations / H. Mukai // Computers & Mathematics with Applications. - 1981. - Vol. 7, no. 3. - P. 235-250.

2. Фельдман, Л.П. Эффективные методы распараллеливания численного решения задачи Коши для обыкновенных дифференциальных уравнений / Л.П. Фельдман, О.А. Дмитриева //Математическое моделирование. - 2001. - Т. 13, № 7. - С. 66-72.

3. Веников, В.А. Кибернетические модели электрических систем: учеб. пособие для вузов / В.А. Веников, О.А. Суханов. - М. : Энергоиздат, 1982. - 328 с.

4. Шишкин, М.В. Моделирование теплогидравлических систем в среде VisSim / М.В. Шишкин, Д.А. Шнайдер // Вестник ЮУрГУ. Серия «Компьютерные технологии, управление, радиоэлектроника». - 2004. - Вып. 3. - № 9(38). - С. 120-123.

5. Басалаев, А.А. Модель системы энергоэффективного теплоснабжения зданий / А.А. Басалаев // Вестник ЮУрГУ. Серия «Компьютерные технологии, управление, радиоэлектроника».

- 2012. - Вып. 16. - № 23(282). - С. 216-218.

Басалаев Александр Анатольевич, аспирант кафедры автоматики и управления, ЮжноУральский государственный университет (г. Челябинск), sasha-basalaev@mail.ru

Bulletin of the South Ural State University Series “Computer Technologies, Automatic Control, Radio Electronics”

2013, vol. 13, no. 2, pp. 82-87

DISTRIBUTED SIMULATION MODELING OF DISTRICT HEAT SUPPLY SYSTEM IN VISSIM

A.A. Basalaev, South Ural State University, Chelyabinsk, Russian Federation, sasha-basalaev@mail. ru

This paper describes an approach to constructing distributed simulation model of district heat supply system using a computer in VisSim simulation environment. The approach is based on division a model into interconnected subsystems, simulation of which can be processed by area network nodes simultaneously. Ways to synchronize model subsystems obtained by the division of initial system are supposed. Advantages and disadvantages of these ways are described. The presented results of simulation of district heat supply systems

with different number of consumers confirm that calculation speed of a model with large number of basic computational units and nonlinear equations increases in case of the division this model into interconnected subsystems.

Keywords: distributed computing, simulation modeling, district heat supply system, VisSim.

References

1. Mukai H. Parallel algorithms for solving systems of nonlinear equations. Computers & Mathematics with Applications, 1981, vol. 7, no. 3, pp. 235-250.

2. Feldman L.P., Dmitriewa O.A. Effective Methods of Multisequencing of Caushi Problem's Numeral Decision for Ordinary Differential Equations [Effektivnye metody rasparallelivanija chislen-nogo reshenija zadachi Koshi dlja obyknovennykh differentsial'nykh uravneniy]. Matematicheskoe modelirovanie [MathematicalModelling], 2001, vol. 13, no. 7, pp. 66-72.

3. Venikov V.A., O.A. Sukhanov. Kiberneticheskiye modeli elektricheskikh sistem: Ucheb. poso-biye dlya vuzov [Cybernetic Models of Electrical Systems: University Textbook]. Moscow, Energoizdat Publ., 1982. 328 p.

4. Shishkin M.V., Shnaider D.A. Simulation of Thermal-Hydraulic Systems in VisSim Environment [Modelirovanie teplogidravlicheskikh sistem v srede Vissim]. Vestnik Juzhno-Ural'skogo Gosu-darstvennogo Universiteta. Serija «Komp'juternye tekhnologii, upravlenie, radioelektronika» [The Bulletin of the South Ural State University. Seria «Computer Technologies, Automatic Control, Radio Electronics»], 2004, no. 9(38), issue 3, pp. 120-123.

5. Basalaev A.A. Model of Energy Efficient Heating for Buildings [Model' sistemy energoeffek-tivnogo teplosnabzhenija zdaniy]. Vestnik Juzhno-Ural'skogo Gosudarstvennogo Universiteta. Serija «Komp'juternye tekhnologii, upravlenie, radioelektronika» [The Bulletin of the South Ural State University. Seria «Computer Technologies, Automatic Control, Radio Electronics»], 2012, no. 23(282), issue 16, pp. 216-218.

Поступила в редакцию 12 марта 2013 г.