CC BY
12
B.C. Беднаржевский
К ВОПРОСУ СОЗДАНИЯ ЕДИНОЙ ИНФОРМАЦИОННОЙ СИСТЕМЫ В ТЕПЛОЭНЕРГЕТИКЕ
Предложена схема информационной технологии для моделирования энергетических паровых котлов как совокупность взаимодействия статических, параметрических и динамических моделей. Дано математическое описание моделей как функции определенных параметров проектирования. Практическая реализация этой схемы проведена на персональных ЭВМ IBM PC в ОАО «Сибэнергомаш», а также на кафедре котло- и реакторостроения Алтайского государственного технического университета.
Математическое моделирование энергетических котлов относится к классу сложных задач ввиду их единичного производства, и, как следствие, при алгоритмизации такого вида задач возникают определенные трудности. Эта проблема может быть решена на основе широкого использования средств вычислительной техники и компьютерных технологий. Основными целями математического моделирования энергетических котлоагрегатов являются:
- сокращение сроков разработки новых котлов и запуска их в производство вследствие более совершенной организации всего цикла проектирования и отладки проектных решений на моделях на ранних стадиях конструирования без изготовления дорогостоящего прототипа;
- повышение качества и конкурентоспособности изделий путем внедрения компьютерных технологий, включая геометрическое моделирование, математические методы анализа и оптимизации будущей конструкции.
Информационное обеспечение компьютерного моделирования паровых котлов включает в себя комплекс пакетов прикладных программ (ППП), состоящий из семи пакетов для проектирования котлоагрегата в режиме диалога: «Тепловой расчет котлоагрегата»; «Расчет температуры металла стенки трубы»; «Расчет на прочность цельносварных газоплотных конструкций»; «Аэродинамический расчет котлоагрегата»; «Расчет теплофизических свойств теплоносителей»; «Гидравлический расчет пароперегревателя»; «Расчет надежности работы котлоагрегата».
Во всех пакетах исходная информация в режиме диалога вводится с клавиатуры дисплея в последовательности, определяемой запросами ЭВМ. Режим диалога обеспечивают программы бесформатного ввода с дисплея числовых и символьных данных. В ходе расчетов на дисплей выводятся промежуточные результаты, на основании анализа которых пользователю предлагается выбрать стратегию поведения: продолжить расчет либо, изменив исходные данные, вернуться к началу программы.
Тепловой расчет - основа математического моделирования котлов. В работах [2, 3] реализован принцип расчета по поверхностям нагрева, несущим
основную тепловую нагрузку. К ним относятся: циркуляционный контур, ширмы (2-ая - 4-ая ступени), пароперегреватель (2-ая - 5-ая ступени), ступени воздухоподогревателя и водяного экономайзера. При компьютерном моделировании паровых котлов для быстрого анализа громоздкие цепочки схем теплообменников можно заменить несколькими теплоне-сущими поверхностями нагрева для упрощения алгоритма теплового расчета и в то же время сохранения качественного анализа конструкции.
При математическом моделировании парового котла последовательно прорабатываются определенные этапы [4]: расчет конструкции, изготовление чертежей, расчет динамики работы котла. Каждому этапу соответствует своя математическая модель: статическая, параметрическая, динамическая.
Статические модели представляют собой алгоритмы расчетов: конструктивного теплового, поверочного теплового, аэродинамического, гидравлического, перепада давления по паровому тракту, элементов котла на прочность, цельносварных экранов на прочность, системы пылеприготовления, трубопроводов на самокомпенсацию, цельносварных газоплотных потолков на прочность, каркаса котла, надежности работы поверхностей нагрева. Ниже каждый расчет укрупненно представлен как функция определенных значений:
- конструктивный тепловой НН=:Т(Р);
- поверочный тепловой ОБ=:Т(НН);
- аэродинамический Z=f(G);
- гидравлический G=f(HH,ТТ);
- температуры стенки трубы Т8=А(ТР, К);
- перепада давления по паровому тракту ZР=f(К,W);
- элементов котла на прочность 8М=:(К,ТТ, 81);
- цельносварных экранов на прочность 8Ю=А(К,ТТ);
- системы пылеприготовления ТК=А(КТ);
- цельносварных газоплотных потолков на прочность БЬ^К^);
- каркаса котла ИР^ДРЯ);
- трубопроводов на самокомпенсацию Ь=А(К, ТТ);
- надежности работы поверхностей нагрева WW=f(NS).
В уравнениях приняты следующие условные обозначения: НН - конструктивные характеристики поверхностей нагрева; Р - исходные данные для
проекта котла (паропроизводительность, вид топлива, температура и давление перегретого пара, температура питательной воды, уходящих газов и воздуха); GD - характеристики потока газов и диаметры труб; Z - сопротивление газовоздушных трактов; G - характеристики потока пароводяной смеси (истинный расход, действительная скорость, кратность циркуляции); ТТ - теплотехнические данные теплового расчета; TS - температура стенки трубы; TP - температура пара; К - конструктивные характеристики трубы (диаметр, толщина стенки); ZP - сопротивление тракта и давление среды; W - скорость пара; SM - минимальная толщина стенки трубы; SI - допустимые напряжения; SIG - суммарные напряжения; TN - температура горячего воздуха и производительность мельниц, сепараторов и т. д.; NT - вид и количество топлива; DL - прогибы потолков и эпюры напряжений; Q - нагрузки; UP - устойчивость и прочность балок каркаса; PR - прочностные характеристики элементов (площадь, момент сопротивления поперечного сечения, момент инерции, нагрузка); WW - параметр потока вынужденных остановов; NS -количество стыков, прямых участков сварных труб; L - удлинения трубопроводов.
Высокий уровень математического моделирования удается обеспечить за счет специализации моделей на определенный тип изделий, таких, например, как трубы, коллектора, цельносварные панели и т. п. Моделирование происходит по AutoLISP-программам в среде AutoCAD, что предполагает предварительную параметризацию объекта проектирования.
Под параметризацией понимается такое представление некоторого класса геометрических объектов, которое позволяет автоматически получать конкретный объект (экземпляр класса) путем задания значений параметров. В крайней точке этого спектра - полностью автоматический выпуск конструкторской и технологической документации, необходимой для подготовки производства и собственно производства котлоагрегата.
Поддержка технологии типового конструирования - один из примеров применения параметризации. Эта технология предполагает, что конструирование детали производится путем сопоставления ей некоторой типовой детали, определения значений размеров и задания дополнительных элементов - «элементов обогащения», которые отличают конструируемую деталь от типовой [5]. Для поддержки этой технологии нужно представить в параметризованной форме и элементы обогащения, и типовые детали. Использование параметризиро-ванных объектов способно радикально изменить стиль работы пользователя в AutoCAD-системах с применением AutoLISP-программ.
На чертеже обычно бывает различное количество сечений, видов, разрезов. Для их построе-
ния целесообразнее использовать параметризацию по размерам. В основе параметризации по размерам лежит известный факт, что геометрия чертежа может быть рассчитана по эскизу и значениям нанесенных на эскизе размеров. С точки зрения параметризации это означает, что любой чертеж с правильной и полной расстановкой размеров может служить описателем класса чертежей, различающихся значениями размеров. Таким образом, для создания параметризированного чертежа достаточно построить эталонный экземпляр-образец. Параметризация по размерам опирается на возможность расчета геометрии чертежа по приблизительному эскизу.
В нем определяется состав геометрических элементов, размерные отношения и наглядные геометрические отношения. Совокупность геометрических и размерных отношений при заданном наборе элементов образует систему ограничений, достаточную для расчета геометрии чертежа.
Все данные, которые используются для представления параметризированного чертежа, составляют его модель. В модели можно выделить две части: параметрическую и топологическую, по которым, соответственно, рассчитывается геометрия для заданных значений размеров и строится экземпляр чертежа. Параметрическая модель МР состоит из трех множеств: Е - геометрических элементов, О - геометрических отношений и Я - размерных отношений; МР={Е, О, Я}. Например, множество Е содержит точки, прямые и окружности. Включим в Е трубы и коллектора, образующие основу для чертежей цельносварных панелей топки, блоков пароперегревателя, ширм, водяного экономайзера, водоопускной системы. Отношение (размерное или геометрическое) определяется типом и набором участвующих в нем элементов. Используем следующие типы отношений.
Геометрические отношения: | | (Т1, Т2) - трубы Т1, Т2 параллельны; ± (Т1, Т2) - трубы Т1, Т2 перпендикулярны; « (Т,К) - торец трубы Т принадлежит коллектору К; « (Т1, Т2) - торец трубы Т1 принадлежит трубе Т2.
Размерные отношения: ф (Т1, Т1) - задан угол гиба трубы Т1; р (Т1, Т2) - задано расстояние между трубами Т1, Т2; р (К1, К2) - задано расстояние между коллекторами К1, К2.
Топологическая модель задает соответствие между элементами чертежа и геометрическими элементами параметрической модели. По топологической модели, используя полученные в результате расчета параметры геометрических элементов, строится экземпляр чертежа.
Многие задачи имеют по несколько решений, и, как правило, для одного набора значений суще-
B.C. Беднаржевский
К вопросу создания единой информационной системы в теплоэнергетике
ствует множество экземпляров, удовлетворяющих одной и той же модели.
Рассмотрим критерий, с помощью которого из этого множества выделяется единственный, в определенном смысле наиболее подходящий, экземпляр. Выбор должен делаться только на основании чертежа-образца и не должен требовать дополнительной информации от пользователя.
Перечислим ряд свойств, которым должен удовлетворять принцип выбора. Пусть по чертежу Е°, который будет использоваться в качестве образца параметризированного чертежа и которому соответствуют значения размеров r0 = (rf,..., rn0), построена модель MP={E0, G, R}. Для любого набора значений у1 процедура выбора должна определять из конечного множества экземпляров, удовлетворяющих модели MP, единственный экземпляр Е1, о котором будем говорить, что он похож на образец Е0, либо сообщать, что похожего экземпляра нет. Процедура выбора задает на множестве экземпляров, удовлетворяющих модели MP, бинарное отношение «похож», обозначаемое далее ». Пользователю было бы очень трудно работать в системе параметризации, в которой принцип выбора не удовлетворял бы следующим свойствам.
Свойство 1 (рефлексивность).
Е»Е.
Выбор для значений размеров, совпадающих со значениями размеров образца Е, должен дать сам образец Е.
Свойство 2 (симметричность).
Е0»Е1^Е1»Е0
Если при использовании образца Е0 с размерами у0 среди экземпляров с размерами у1 выбран экземпляр Е1, то при использовании образца Е1 среди экземпляров с размерами у0 должен быть выбран Е0.
Свойство 3 (транзитивность).
Е0»Е1, Е1»Е2^Е0»Е2
Свойство 4 (однозначность).
Е1»Е2, у (Е1)= у (Е2) ^Е1°Е2
Здесь у (Е) обозначает вектор значения размеров для экземпляра Е. Это свойство делает однозначным выбор при одинаковых значениях размеров. В частности, этим обеспечивается то, что начиная с двух различных, но похожих образцов при любых одинаковых значениях размеров мы получим один и тот же экземпляр.
Свойство 5 (непрерывность). Образец и выбранный экземпляр должны отличаться друг от друга незначительно при близких значениях размеров.
Описанный аппарат параметризации был реализован для построения рабочих чертежей коллекторов на языке AutoLISP для ПЭВМ IBM PC [6]. Предварительно была проведена классификация всех коллекторов. На полученном рабочем черте-
же показаны виды и разрезы торцов, сечения, а также элементы обогащения - штуцеры, транспортировочные ушки, торцы, отверстия. При программной реализации на ПЭВМ на укрупненную параметрическую модель коллектора MP={E, О, R} накладываются геометрические и размерные ограничения (до пяти рядов отверстий, до трех рядов штуцеров, длина коллектора не больше 4500 мм, 36 комбинаций левый-правый торцы). В зависимости от количества сечений и видов программа сама выбирает формат чертежа с основной надписью, проставляются размеры, допуски, пишется текст условий изготовления коллектора.
Динамическая модель котла необходима для настройки автоматических систем регулирования (САР) и выявления опасных отклонений параметров на стадии проектирования. Задача математического моделирования заключается в составлении системы линейных дифференциальных и алгебраических уравнений с постоянными коэффициентами, связывающей входные и выходные координаты расчетных участков и в решении этой системы уравнений относительно выходных координат с целью получения динамических характеристик на ЭВМ.
Математическая модель создается в предположении, что котельный агрегат является линейной детерминированной системой в условиях малых возмущений (сосредоточенные параметры) [7].
При составлении уравнений динамики применяется метод линеаризации, позволяющий использовать при исследованиях малые отклонения от стационарного режима.
Барабанный котел делится на ряд расчетных участков, каждый из которых неизменен в геометрическом отношении.
Необогреваемые трубопроводы и коллекторы присоединяются к соответствующим расчетным участкам, объемы и массы этих элементов прибавляются к объемам и массам соответствующих расчетных участков.
По характеру теплообмена участки барабанного котла разделяются на конвективные и радиационные, для радиационных участков принято, что тепловой поток не зависит от температуры рабочей среды.
За радиационный участок принимается циркуляционный контур барабанного котла, все остальные поверхности нагрева котла считаются конвективными, но изменение коэффициента теплоотдачи лучеиспусканием учитывается по соответствующим зависимостям.
При определении величины аккумулирования тепла в циркуляционном контуре в расчет вводится величина активного металла барабанов и коллекторов, изменение температуры которого следует за изменением температуры насыщенного пара.
Изменение температуры металла стенок труб и барабанов циркуляционного контура принимается равным изменению температуры насыщенного пара при данном давлении.
Объем пара в самом барабане под зеркалом испарения и опускных трубах не учитывается ввиду его малого численного значения по сравнению с объемом барабана и объемом пара в трубной системе циркуляционного контура.
Температура и энтальпия воды в циркуляционном контуре котла принимается равной температуре и энтальпии воды при давлении насыщенного пара в барабане.
Определение величины материального и теплового аккумулирования на всех расчетных участках (кроме циркуляционного контура) производится по плотности и энтальпии рабочей среды на выходе из участка.
Массовая и тепловая аккумуляция в газовоздушном тракте не учитывается.
Давление дымовых газов принято постоянным.
Коэффициенты теплоотдачи от газов к стенкам труб и от стенок труб к рабочей среде приняты пропорциональными изменению скоростей соответствующих сред.
Коэффициенты теплопроводности рабочей среды, критерий Прандтля и коэффициенты температуропроводности приняты постоянными.
Топка рассматривается как одно звено, зоны по высоте топки в расчет не принимаются, процессы горения в топке не моделируются, рассматриваются только процессы теплообмена, при этом принято, что процесс теплообмена в топке не оказывает влияния на параметры рабочей среды, протекающей внутри экранирующих труб, т. е. все поверхности в топке принимаются радиационными.
Температура горячего воздуха, поступающего в топку, принята постоянной, в случае необходимости выявления влияния изменения температуры горячего воздуха на динамические характеристики котла изменение температуры этого воздуха можно рассматривать как внешнее возмущение.
Линеаризация переменных производится относительно значений переменных величин в исходном стационарном состоянии. Параметры, характеризующие исходное состояние, а также конструктивные параметры принимаются постоянными для каждого режима.
Барабанный котел представляет собой динамическую систему, которая при моделировании разделяется на отдельные взаимосвязанные звенья, образующие поверхности нагрева и необогре-ваемые участки.
Барабанный котел при расчете динамических характеристик разбивается на 2 группы участков, к первой группе относится циркуляционный контур котла как участок с двухфазной средой. Этот контур является многомерной динамической системой, входными координатами которой являются изменения давления, температуры и расхода рабочей среды, а выходными координатами - изменения давления пара, уровня в барабане, расхода пара из барабана, температура газов на выходе из топки.
Ко второй группе относятся все участки котла с однофазной средой (пароперегреватель).
Теплообмен между средами и стенкой характеризуется коэффициентами теплоотдачи. Коэффициент теплоотдачи от внутренней стенки труб к рабочей среде не зависит от расхода, температуры и давления рабочей среды.
Коэффициент теплоотдачи от газов к наружной поверхности разделяющей стенки имеет конвективную и радиационную составляющие.
Конвективная составляющая зависит от расхода топлива и воздуха, средней температуры газов и коэффициента рециркуляции.
Радиационная составляющая зависит от температуры газов и температуры загрязнения стенки трубы.
Статические отклонения параметров определяют по специальной программе на персональном компьютере.
Описанная теория была применена для кот-лоагрегата БКЗ 420-140-9 производства ОАО «Си-бэнергомаш» [8].
Список использованной литературы:
1. Беднаржевский B.C. Комплекс пакетов прикладных программ для автоматизированного проектирования котлоагрегата на APM-M / / Тяжелое машиностроение, 1992, №12, с. 33-35.
2. Беднаржевский B.C. Автоматический тепловой расчет котлоагрегата на ЭВМ //Изв. вузов «Энергетика», 1995, №1-2, с. 54-57.
3. Беднаржевский В.С., Оскорбин Н.М. Основные положения теплового расчета паровых котлов на ЭВМ // Теплоэнергетика. 2002. №8. c.10-14.
4. Беднаржевский В.С., Оскорбин Н.М. Автоматизированное проектирование энергетических паровых котлов // Известия вузов. Проблемы энергетики. 2002. №1-2. С. 3-9.
5. Коваленко В.Н., Ревякин Ю.Г., Хухлаев Е.В. Параметризация машиностроительных чертежей, основанная на поэлементном расчете // Программирование. 1992. №2, С. 64-77
6. Беднаржевский B.C. Автоматизированное проектирование коллекторов энергетических котлов на ПЭВМ IBM PC // Тяжелое машиностроение, 1994, №11-12. С.13-14, 32.
7. РТМ 108.031.101-84. Котлы барабанные: Расчет динамических характеристик. - Л.: НПО ЦКТИ, 1986. - 80 с.
8. Беднаржевский В.С., Оскорбин Н.М. Динамическая математическая модель парового котла БКЗ 420-140-9 // Известия вузов. Проблемы энергетики. 2002. №3-4. С. 8-13.
