Математическое моделирование и компьютерные технологии в задачах проектирования энергетических паровых котлов
Автор: В.С. Беднажевский
Источник: Математическое моделирование и компьютерные технологии
в задачах проектирования энергетических паровых котлов. Том 7, №6, 2002
с. 13-20
Математическое моделирование является одним из наиболее быстро развивающихся направлений научно-технического прогресса, современные темпы которого обусловливают необходимость повышения качества и сокращения продолжительности проектирования новых видов продукции.
Математическое моделирование энергетических котлов относится к классу сложных задач ввиду их единичного производства, и, как следствие, при алгоритмизации такого вида задач возникают определенные трудности. Эта проблема может быть решена на основе широкого использования средств вычислительной техники и компьютерных технологий. Основными целями математического моделирования энергетических котлоагрегатов являются: сокращение сроков разработки новых котлов и запуска их в производство вследствие более совершенной организации всего цикла проектирования и отладки проектных решений на моделях на ранних стадиях конструирования без изготовления дорогостоящего прототипа; повышение качества и конкурентоспособности изделий путем внедрения компьютерных технологий, включая геометрическое моделирование, математические методы анализа и оптимизации будущей конструкции.
Комплекс пакетов прикладных программ для моделирования котлоагрегатов
Информационное обеспечение компьютерного моделирования паровых котлов включает в себя комплекс пакетов прикладных программ, состоящий из семи пакетов прикладных программ (ППП) для проектирования котлоагрегата в режиме диалога: «Тепловой расчет котлоагрегата», «Расчет температуры металла стенки трубы», «Расчет на прочность цельносварных газоплотных конструкций», «Аэродинамический расчет котлоагрегата», «Расчет теплофизических свойств теплоносителей», «Гидравлический расчет пароперегревателя», «Расчет надежности работы котлоагрегата».
В ППП теплового расчета котлоагрегата входят вычисления: объемов газообразных продуктов сгорания, объемных долей трехатомных газов, концентрации золы, таблицы энтальпий теоретического расхода воздуха и газообразных продуктов сгорания, теплового баланса котлоагрегата, топки, ширмовых пароперегревателей, конвективных пароперегревателей, хвостовых поверхностей нагрева.
Алгоритмы расчета теплового режима объединены с инвертированной базой данных «Котлоагрегат», позволяющей по известному значению ключевого признака определить адреса всех записей, которые его содержат.
Пакет прикладных программ «Расчет температуры металла стенки трубы» предназначен для определения средней температуры металла стенки неравномерно обогреваемой трубы ширмового и конвективного пароперегревателей.
Пакет прикладных программ «Расчет на прочность цельносварных газоплотных конструкций» предназначен для определения теплонапряжений в цельносварных газоплотных экранах котла в зоне максимального уровня тепловых потоков с учетом дополнительных нагрузок. Пакет производит расчет напряжений от внутреннего давления, весовых нагрузок, хлопка в топке, неравномерного распределения тепловых потоков по ширине экранов, разности температур экрана и балки жесткости, находит суммарные напряжения и оценивает прочность
Пакет прикладных программ «Аэродинамический расчет котлоагрегата» предназначен для вычисления аэродинамических сопротивлений газового и воздушного трактов котельного агрегата, он позволяет для газового тракта котлоагрегата рассчитывать следующие поверхности: пароперегреватель, фестон, водяной экономайзер, воздухоподогреватель. Поток газов может смывать шахматные, коридорные, продольные пучки. Предусмотрен расчет участков от воздухоподогревателя до дымовой трубы, дымовая труба.
Для всех участков рассчитываются повороты в канале, изменение сечения, раздающие и собирающие короба, диффузоры. Для воздушного тракта котельной установки можно учесть следующие виды сопротивлений: измерительная шайба, поворот на N градусов, карман Рихтера, трение на участке, диффузор, изменение сечения, калорифер, воздухонагреватель, конфузор, сопротивление входа в трубу, шибер, сопла, горелки, собирающий тройник.
Пакет прикладных программ «Расчет теплофизических свойств теплоносителей» предназначен для расчета термодинамических и теплофизических свойств воды, водяного пара и дымовых газов. В состав пакета входят 17 программ. Программы расчета энтальпий газов по температуре, критерия Прандтля для дымовых газов, теплопроводности дымовых газов, кинематической вязкости дымовых газов применимы для температуры 0 ... 3000 °С. Программы расчета температуры пара по энтальпии и давлению, энтальпии пара по давлению и температуре применяются при температуре 330 ... 560 °С и давлении 10 ... 20 МПа. Программы расчета энтальпии воды по температуре и давлению, температуры воды по энтальпии и давлению, удельного объема воды используются при температуре 215 ... 320 °С и давлении 10 ... 20 МПа. Расчет удельного объема воды применим для всей области до критического давления. Программы расчета удельного объема пара на кривой насыщения, энтальпии воды на линии насыщения, энтальпии пара на линии насыщения, температуры воды и водяного пара на линии насыщения, удельного объема воды на линии насыщения применяются при температуре насыщения. Программы расчета динамической вязкости воды и водяного пара, теплопроводности воды и водяного пара применяются для температуры 0 ... 800 °С и давления до 100 МПа.
Пакет прикладных программ «Гидравлический расчет пароперегревателя» предназначен для расчета перепада давления по паровому тракту пароперегревателя котлоагрегата. Основные программы пакета: расчет удельного объема пара, внутреннего сечения труб, перепада давления, скорости пара, гидравлических сопротивлений входа, выхода, поворота, трения по длине труб, вывода на печать результатов расчета.
На этапе моделирования необходимо кроме основных расчетов (теплового, гидравлического, прочностного и т.д.) проводить вероятностную оценку надежности каждого варианта конструкции котлоагрегата. Пакет прикладных программ «Расчет надежности котлоагрегата» предназначен для расчета надежности работы поверхностей нагрева котлоагрегата. Параметр потока отказов по каждой поверхности рассчитывается как сумма произведений числа элементов на интенсивность отказов. Элемент складывается из прямых участков, бесштуцерной приварки змеевиков к камерам, гибов труб, контактных сварных стыков. Также определяется характеристика потока вынужденных остановов на 105 ч работы.
Во всех пакетах исходная информация в режиме диалога вводится с клавиатуры в последовательности, определяемой запросами ЭВМ. Режим диалога обеспечивают программы бесформатного ввода с дисплея числовых и символьных данных. В ходе расчетов на дисплей выводятся промежуточные результаты, на основании анализа которых пользователю предлагается выбрать стратегию поведения: продолжить расчет либо, изменив исходные данные, вернуться к началу программы. Тепловой расчет – основа математического моделирования котлов. В работе реализован принцип расчета по поверхностям нагрева, несущим основную тепловую нагрузку. К ним относятся: циркуляционный контур, ширмы (2 – 4-я ступени), пароперегреватель (2 – 5-я ступени), ступени воздухоподогревателя и водяного экономайзера. При компьютерном моделировании паровых котлов для быстрого анализа громоздкие цепочки схем теплообменников можно заменить несколькими теплонесущими поверхностями нагрева для упрощения алгоритма теплового расчета и в то же время сохранения качественного анализа конструкции.
При математическом моделировании парового котла последовательно прорабатываются определенные этапы: расчет конструкции, изготовление чертежей, расчет динамики работы котла. Каждому этапу соответствует своя математическая модель: статическая, параметрическая, динамическая.
Статические модели
Статические модели представляют собой алгоритмы расчетов: конструктивного теплового, поверочного теплового, аэродинамического, гидравлического, перепада давления по паровому тракту, элементов котла на прочность, цельносварных экранов на прочность, системы пылеприготовления, трубопроводов на самокомпенсацию, цельносварных газоплотных потолков на прочность, каркаса котла, надежности работы поверхностей нагрева. Ниже каждый расчет укрупненно представлен как функция определенных значений:
В уравнениях приняты следующие условные обозначения: НН − конструктивные характеристики поверхностей нагрева; Р − исходные данные для проекта котла (паропроиз водительность, вид топлива, температура и давление перегретого пара, температура питательной воды, уходящих газов и воздуха); СО − характеристики потока газов и диаметры труб; Z − сопротивление газовоздушных трактов; С − характеристики потока пароводяной смеси (истинный расход, действительная скорость, кратность циркуляции); ТТ − теплотехнические данные теплового расчета; ТБ − температура стенки трубы; ТР − температура пара; К − конструктивные характеристики трубы (диаметр, толщина стенки); ZP − сопротивление тракта и давление среды; Ш − скорость пара; БМ − минимальная толщина стенки трубы; Б 1 − допустимые напряжения; Б 1С − суммарные напряжения; TN − температура горячего воздуха и производительность мельниц, сепараторов и т.д. NT − вид и количество топлива; DL − прогибы потолков и эпюры напряжений; Q − нагрузки; UP − устойчивость и прочность балок каркаса; PR − прочностные характеристики элементов (площадь, момент сопротивления поперечного сечения, момент инерции, нагрузка); WW − параметр потока вынужденных остановов; NS − количество стыков, прямых участков сварных труб; L − удлинения трубопроводов.
Динамическая модель
Динамическая модель котла необходима для настройки автоматических систем регулирования (САР) и выявления опасных отклонений параметров на стадии проектирования. Задача математического моделирования заключается: а) в составлении системы линейных дифференциальных и алгебраических уравнений с постоянными коэффициентами, связывающей входные и выходные координаты расчетных участков; б) в решении этой системы уравнений относительно выходных координат с целью получения динамических характеристик на ЭВМ.
Математическая модель создается в предположении, что котельный агрегат является линейной детерминированной системой в условиях малых возмущений (сосредоточенные параметры). При составлении уравнений динамики применяется метод линеаризации, позволяющий использовать при исследованиях малые отклонения от стационарного режима.
Барабанный котел делится на ряд расчетных участков, каждый из которых неизменен в геометрическом отношении.
Необогреваемые трубопроводы и коллекторы присоединяются к соответствующим расчетным участкам, объемы и массы этих элементов прибавляются к объемам и массам соответствующих расчетных участков.
По характеру теплообмена участки барабанного котла разделяются на конвективные и радиационные, для радиационных участков принято, что тепловой поток не зависит от температуры рабочей среды.
За радиационный участок принимается циркуляционный контур барабанного котла, все остальные поверхности нагрева котла считаются конвективными, но изменение коэффициента теплоотдачи лучеиспусканием учитывается по соответствующим зависимостям.
При определении величины аккумулирования тепла в циркуляционном контуре в расчет вводится величина активного металла барабанов и коллекторов, изменение температуру которого следует за изменением температуры насыщенного пара. Изменение температуры металла стенок труб и барабанов циркуляционного контура принимается равным изменению температуры насыщенного пара при данном давлении.
Объем пара в самом барабане под зеркалом испарения и опускных трубах не учитывается ввиду его малого численного значения по сравнению с объемом барабана и объемом пара в трубной системе циркуляционного контура. Температура и энтальпия воды в циркуляционном контуре котла принимаются равными температуре и энтальпии воды при давлении насыщенного пара в барабане.
Определение величины материального и теплового аккумулирования на всех расчетных участках (кроме циркуляционного контура) производится по плотности и энтальпии рабочей среды на выходе из участка. Массовая и тепловая аккумуляция в газовоздушном тракте не учитывается. Давление дымовых газов принято постоянным. Коэффициенты теплоотдачи от газов к стенкам труб и от стенок труб к рабочей среде приняты пропорциональными изменению скоростей соответствующих сред.
Коэффициенты теплопроводности рабочей среды, критерий Прандтля и коэффициенты температуропроводности приняты постоянными.
Топка рассматривается как одно звено, зоны по высоте топки в расчет не принимаются, процессы горения в топке не моделируются, рассматриваются только процессы теплообмена, при этом принято, что процесс теплообмена в топке не оказывает влияния на параметры рабочей среды, протекающей внутри экранирующих труб, т.е. все поверхности в топке принимаются радиационными.
Температура горячего воздуха, поступающего в топку, принята постоянной, в случае необходимости выявления влияния изменения температуры горячего воздуха на динамические характеристики котла изменение температуры этого воздуха можно рассматривать как внешнее возмущени.
Выводы
Предложена схема компьютерной технологии для моделирования энергетических паровых котлов как совокупность взаимодействия статических, параметрических и динамических моделей. Дано математическое описание моделей как функции определенных параметров проектирования. Практическая реализация этой схемы проведена на кафедре котло и реакторостроения Алтайского государственного технического университета.
Список использованной литературы
1. Беднаржевский В.С. Комплекс пакетов прикладных программ для автоматизированного проектирования котлоагрегата на АРМ-М // Тяжелое машиностроение. 1992. №12. С. 33-35.
2. Беднаржевский В.С. Автоматический тепловой расчет котлоагрегата на ЭВМ // Изв. вузов. Энергетика. 1995. №1-2. С. 54-57.
3. Беднаржевский В.С., Оскорбин Н.М. Основные положения теплового расчета паровых котлов на ЭВМ // Теплоэнергетика. 2002. №7. С. 10-14.
4. РТМ 108.031.101-84. Котлы барабанные: Расчет динамических характеристик. Л.: НПО ЦКТИ, 1986. 80 с.