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

К1 – 1-ая ступень компрессора.
К2 – 2-ая ступень компрессора.
ОВ – охладитель воздуха.
Т1–1-ая камера сгорания.
Т2 – 2-ая камера сгорания.
ГТ1 – 1-ая ступень газовой турбины.
ГТ2 – 2-ая ступень газовой турбины.
ЭГ – электрогенератор.
В – воздух.
Вд – вода.
Т – топливо.
д.г. – дымовые газы.
1. Формализованное описание схемы.
Для описания схемы составим граф вида:

Которую опишем с помощью структурной матрицы (матрица связи).
| I | II | III | IV | V | VI | VII | VIII |
| |
| -1 | |||||||||
| -1 | -1 | ||||||||
| -1 | |||||||||
| -1 | |||||||||
| -1 | |||||||||
| -1 | |||||||||
| -1 | |||||||||
| -1 | -1 | ||||||||
| -1 | |||||||||
| -1 | |||||||||
| -1 | |||||||||
| -1 | |||||||||
| -1 | -1 |
Входящую связь обозначаем - +1, выходящую – -1.
Матрица видов связи:
| Взд | вода | т-во | ПС | МЭ | ЭЭ | |
ПС – продукты сгорания.
МЭ – механическая энергия.
ЭЭ – электрическая энергия.
2. Дополнительно составляем матрицы: смежности, пр-сов, контуров.
Матрица смежности:
| I | II | III | IV | V | VI | VII | VIII | |
| I | ||||||||
| II | ||||||||
| III | ||||||||
| IV | ||||||||
| V | ||||||||
| VI | ||||||||
| VII | ||||||||
| VIII | ------ | ------ | ------ | ------ | ----- | ------- | ------ | ------ |
Имеется пустая VIII строка, это значит, что связь 8 не входит в контур, и ее можно рассчитать отдельно.
Матрица пр-сов: Матрица контуров:
| I | -2 | ||||
| II | -4 | -5 | |||
| III | -6 | ||||
| IV | -8 | ||||
| V | -9 | -13 | -14 | -15 | |
| VI | |||||
| VII | -12 | -16 | |||
| VIII | -17 |
| |||||||||||
3. Сис-ма балансовых у-ний.
Составим упрощенную сис-му у-ний, состоящую только из материальных и энергетических балансов.
Сис-му сводим в таблицу:
| № эл-та | Графическое у-ние | Балансовое у-ние | примечание |
| I |
|
|
(1)
|
| II |
|
|
(2)
|
| III |
|
|
(3)
|
| IV |
|
| (4) (5) |
| V |
|
|
(6)
|
| VI |
|
| (7) (8) |
| VII |
|
|
(9)
|
| VIII |
|
| (10) |
Параллельно в записи у-ний рекомендуется составлять таблицу переменных, в следующей форме:
| № эл-та | наименование | единицы измерения | обозначения | идентификатор | диапазон измерений | примечание |
4. сис-ма ограничения.
В данном разделе указываются диапазоны измерения основных переменных, и тем самым определяется область исследования.
5. ф-ции качества.
В качестве критерия эффективности принимаем КПД установки, и в этом случае ф-ция качества будет представлена у-нием:
.
6. анализ переменных мат. модели и выбор зависимых и независимых пар-ров.
Расходы:
.
Энтальпия:
.
Мощности:
.
Коэф. потерь:
У нас 8 эл-ов, поэтому будет 8 значений потерь.
Общее число переменных:
.
Степень свободы:
.
В 1-ую очередь задаются конструктивные пар-ры, к которым в данном случае можно отнести коэф. потерь, далее задаются теплофизические пар-ры (т-ра и давление) по которым находим энтальпию. Часто приходиться для этого вводить дополнительное у-ние.
В частности, в данном случае, т-ра на выходе из компрессора или газовой турбины м.б. однозначно определена, если заданная т-ра на входе и соответственно степень сжатия или расширения. Из группы расходов и мощностей задается минимальное число переменных, обычно это или расход т-ва или мощ-ть.
В данном случае мы можем определить электрическую мощ-ть установки, и тогда расходы т-ва будут зависимыми, или можем произвести расчет в предельных единицах, задав расход воздуха равного единице. В рез-те получим удельную электрическую мощ-ть и удельные расходы т-ва отнесенный на единичный расход т-ва.
В табличных переменных в примечании указывается, какие переменные явл-ся зависимыми, какие нет, какие регламентированные какие нет.
Если G1 = 1 кг/с, то остается 10 у-ний, все энтальпии необходимо задавать.
Для определения, из какого у-ния находим, неизвестный пар-р, находит сокращенную матрицу функциональных связей, где строки считаются по числу у-ний, столбцы – по неизвестным пар-рам.
|
|
|
|
|
|
|
|
|
|
| |
| |||||||||||
| |||||||||||
| |||||||||||
| |||||||||||
| |||||||||||
| |||||||||||
| |||||||||||
| |||||||||||
| |||||||||||
| |||||||||||
1-ца в i-j эл-те матрицы означает, что j-ая переменная входит в i-ое у-ние. Присматриваем все столбцы, в которых только одна единица. Далее просматриваем строки, и ищем, где есть только одна 1-ца.
Если в строке находится только одна 1-ца, то данное у-ние м.б. рассчитано отдельно от всей сис-мы.
Алгоритм повторяется, но уже по 2 переменным (по двум единицам).
(1)
(2)
(3)
(6)
(9)






