Учебная работа. Моделирование тепловой схемы ПГУ

1 Звезда2 Звезды3 Звезды4 Звезды5 Звезд (Пока оценок нет)
Загрузка...

Моделирование тепловой схемы ПГУ

Содержание

Введение

1. Технологический процесс

1.1 Применение газовой турбины в современной промышленности

1.2 краткое описание энергоблока ПГУ-190/220 Тюменской ТЭЦ-1

1.3 Выбор объекта моделирования и постановка задачи

1.4 Основные показатели работы камер сгорания

2. Программная среда разработки проекта

2.1 краткое описание метода конечных элементов

2.2 Основные принципы проектирования в ANSYS

предварительная подготовка и вход в программу

Основные стадии решения задач

приложение нагрузок и получение решения

Постпроцессорная обработка

Типы основных файлов, создаваемых программой

3. Проектирование модели исследуемого объекта

3.1 Построение графической модели

3.2 Наложение конечно-элементной сетки

3.3 Задание начальных условий и нагрузок

3.4 Подготовка и проведение расчетов

3.5 Определение оптимального коэффициента избытка воздуха

Выводы по главе

4. безопасность и экологичность проекта

4.1 Характеристика условий труда

4.2 Расчет искусственного освещения

Источники, виды и системы искусственного освещения

Расчет искусственного освещения с помощью коэффициента использования светового потока

4.3 Оценка экологичности проекта

4.4 Оценка чрезвычайных ситуаций

Выводы по разделу

5. кономическая эффективность проекта

5.1 Источники экономической эффективности

5.2 Расчет единовременных затрат

5.3 Оценка экономической целесообразности и выводы по разделу

Заключение

список использованных источников

Приложения

Введение

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

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

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

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

В представленном дипломном проекте разработана методика газодинамического анализа кольцевой камеры сгорания с использованием инженерного пакета ANSYS.

Методика включает в себя:

камера сгорание газовая турбина

разработку плоской графической модели камеры в виде радиального сечения в прикладном графическом пакете SolidWorks 2004;

наложение конечно-элементной сетки в ANSYS на импортированную в него графическую модель камеры;

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

представление результатов решения в виде анимации и диаграмм;

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

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

1. Технологический процесс

Благодаря повсеместному переходу в 90-е годы на использование природного газа в качестве основного топлива для электроэнергетики, газовые турбины заняли существенный сегмент рынка. Несмотря на то, что максимальная эффективность оборудования достигается на мощностях от 5 МВт и выше (до 300 МВт), некоторые производители выпускают модели в диапазоне 1-5 МВт.

значительная роль в техническом перевооружении должна принадлежать программе широкого внедрения в электроэнергетику газотурбинных и парогазовых технологий [14,18, 19]. В области создания газотурбинных установок (ГТУ) отечественная промышленность в основном находится на стадии разработок и испытаний опытных образцов и пока существенно отстает от мирового уровня.

.1 Применение газовой турбины в современной промышленности

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

Рисунок 1.1 — Газотурбинная энергетическая установка ГТЭ-110

Внедрение парогазовых и газотурбинных технологий обеспечит повышение КПД установок до 55 %, а в перспективе до 60 % и более. Это позволит существенно снизить прирост потребности ТЭС в топливе, улучшить экологические показатели электростанций, резко снизить капитальные затраты на новое строительство и техническое перевооружение.

продолжаются работы по созданию разных видов ГТУ [19]. В настоящее время РАО «ЕЭС России» совместно с предприятиями ГП НПКГ «Зоря — Машпроект» и ОАО «Сатурн» завершают на испытательном стенде Ивановской ГРЭС испытания опытного образца газотурбинной энергетической установки ГТЭ-110, изображенной на рисунке 1.1, которая станет основой современных отечественных парогазовых энергоблоков ПГУ-325.

В 2002 г. по заказу РАО «ЕЭС россии» ОАО ЛМЗ совместно с ОАО «Авиадвигатель» завершили разработку технического проекта газотурбинной установки ГТЭ-180. На базе ГТЭ-180 будут модернизироваться существующие газомазутные энергоблоки 300 МВт, а также разрабатываться пилотные проекты для ГРЭС-3 «Мосэнерго», Пермской ГРЭС, ТЭЦ-5 «Ленэнерго».

На Северо-Западной ТЭЦ «Ленэнерго» продолжается опытно-промышленная эксплуатация первого парогазового энергоблока мощностью 450 МВт. Начато строительство Калининградской ТЭЦ-2 с подобными парогазовыми установками. Принято решение о строительстве Сочинской ТЭЦ мощностью 76 МВт.

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

Для технического перевооружения КЭС и ТЭЦ имеются, в основном, следующие две альтернативы [17]:

применение самых современных технологий производства энергии — нового поколения газовых турбин (ГТ), паровых турбин (ПТ), парогазовых установок (ПГУ) и котлов — в первую очередь на имеющихся площадках ТЭС, что позволяет использовать существующие линии электропередачи для вывода мощности в энергосистему и использовать уже имеющиеся разрешения на размещение электростанции на существующей площадке;

использование имеющегося технического потенциала модернизации оборудования (паровых турбин, котлов, АСУТП и др.) в соответствии с современным уровнем техники.

Первый вариант означает последовательный вывод из эксплуатации существующих морально устаревших и неэкономичных установок с параллельным строительством на площадке энергоблоков с самыми современными техническими решениями. Примером такой стратегии может служить модернизация всей структуры электростанций бывшей ГДР в период с 1992 по 2000 гг.

Решающим для правильного выбора стратегии технического перевооружения и модернизации ТЭС является выбор топлива при условии, что будут действовать рыночные ценовые соотношения на газ, уголь и мазут, включая транспортные расходы. При этом решающую роль играют размеры страны. Так, например, в такой относительно небольшой стране, как Германия, с развитой инфраструктурой можно почти к любой площадке электростанции подвести газ, уголь и мазут и отвести от нее электроэнергию в сеть. А в такой большой стране, как Россия, как правило, доступен и экономически обоснован лишь один вид топлива — газ или уголь. Для ТЭС, работающих на газовом топливе в Европейской части РФ (доля газа 83%) переход на уголь по причине огромных расстояний от мест добычи угля экономически необоснован, учитывая при этом рыночные транспортные тарифы. Применение газа на городских ТЭЦ, прежде всего, обусловлено соображениями экологии.

При модернизации ТЭС на газе экономически целесообразно использовать самые современные ПГУ с КПД до 58% [18].

одной из особенностей электроэнергетики России является то, что на паросиловых электростанциях получила большое распространение комбинированная выработка тепла и электроэнергии [14]. В настоящий момент в РАО «ЕЭС России» на КЭС и ТЭЦ установлено свыше 132 000 МВт электрической мощности и более 187 000 МВт тепловой мощности. Здесь имеются возможности технического перевооружения ПСУ (паросиловая установка) — ТЭЦ в ПГУ (парогазовая установка) — ТЭЦ. При сохранении тепловой нагрузки электрическая мощность ТЭЦ с применением парогазовой технологии может увеличиться в 2-2,5 раза. Исходя из этого, концепцию технического перевооружения и модернизации ТЭС РАО «ЕЭС россии» целесообразно строить в строгом соответствии с топливным признаком, с учетом особенностей комбинированной выработки электроэнергии и тепла на различных типах турбин.

Принцип работы газовых турбин состоит в следующем: газ, нагнетаемый в камеру сгорания компрессором, смешивается с воздухом, формируя топливную смесь, и поджигается. Образующиеся продукты горения с высокой температурой (900°С-1200°С), проходя через несколько рядов лопаток, установленных на валу турбины, приводят к вращению турбины. Механическая энергия вала передается через понижающий редуктор электрическому генератору. Тепловая энергия выходящих из турбины газов поступает в теплоутилизатор. Вместо производства электричества, механическая энергия турбины может использоваться для работы насосов, компрессоров и т.п. Наиболее традиционным видом топлива для газовых турбин является природный газ, хотя это не исключает возможности использования других видов газообразного топлива. При этом газовые турбины предъявляют повышенные требования к качеству его подготовки (механические включения, влажность).

температура исходящих из турбины газов составляет 450°С-550°С. Количественное соотношение тепловой энергии к электрической у газовых турбин составляет от 1.5: 1 до 2.5: 1, что позволяет строить когенерационные системы, различающиеся по типу теплоносителя:

непосредственное (прямое) использование отходящих горячих газов;

Производство пара низкого или среднего давления (8-18 кг/см2) во внешнем котле;

производство горячей воды (лучше, когда требуемая температура превышает 140°С);

производство пара высокого давления.

КПД газовой турбины составляет 25% — 35%, в зависимости от параметров работы конкретной модели турбины и характеристик топлива. В составе когенерационных систем эффективность возрастает до 90% в расчете на условную единицу израсходованного топлива (по теплотворной способности). Газовые турбины обладают хорошими экологическими параметрами (эмиссия NOx на уровне 25 ppm).

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

.2 краткое описание энергоблока ПГУ-190/220 Тюменской ТЭЦ-1

Сочетание паротурбинных и газотурбинных установок, объединенных общим технологическим циклом, позволило снизить потери теплоты с уходящими газами ГТУ, использовать газы за газовыми турбинами в качестве подогретого окислителя при сжигании топлива в паровом котле, получить дополнительную тепловую и электрическую мощность за счет частичного вытеснения регенерации паровых турбоустановок и, в конечном итоге, повысить КПД электростанции [17, 20].

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

Технические характеристики основного тепломеханического оборудования блока № 1 при температуре охлаждающего воздуха — 7,5°С (среднеотопительная) приведены в таблице 1.

Таблица 1 — Основное тепломеханическое оборудование ПГУ-190/220

№№ ппНаименованиеЕдиницы измеренияХарактеристика1Паровой котлоагрегатЕ-500-13,8-560ГН (модель ТГЕ-435/ПГУ) 1.1Номинальная паропроизводительностьт/ч5001.2Номинальные параметры свежего пара: давление температура МПа (кгс/см2) 0С 13,8 (135) 5601.3температура питательной воды0С1681.4Температура уходящих газов0С1071.5КПД при номинальной паропроизводительности%94,171.6Массовая концентрация оксида азота NOх в дымовых азах при а = 1,4 и содержании NOх за ГТУ 50 мг/м3, при 02 = 15%. мг/м31251.7Массовая концентрация окиси углерода СО в дымовых газах при а =1,4 и выбросе СО за газовой турбиной 0 мг/м3мг/м33002Паровая турбинаТ-130/160-12,82.1максимальная мощность в конденсационном режимеМВт163,52.2Номинальная мощность в среднеотопительном режиме (-7,5°С) при тепловой нагрузке 140,5 Гкал/часМВт150,82.3Номинальный расход свежего парат/ч5002.4Номинальное абсолютное давление свежего пара перед топорными клапанами ЦВДМПа (кгс/см2) 12,8 (130) 2.5Номинальная температура свежего пара перед стопорными клапанами ЦВД0С5552.6Номинальный расход охлаждающей воды конденсаторам3 /час160002.7Удельный расход теплоты в конденсационном режимеккал/кВтч20373. Газовая турбина фирмы SiemensV 64. ЗА3.1Температура газов за газовой турбиной0С574,53.2Расход газов на выходе из газовой турбиныкг/с2003.3Расход топливакг/с (м3/ч) 4,17 (15,012) 3.4электрическая мощностьМВт70,853.5КПД ГТУ%35,123.6Концентрация NOх в газах за турбиноймг/м350

Паровой котел типа Е-500-13,8-560ГН (модель ТГЕ-435/ПГУ) предназначен для работы в составе парогазовой установки с газовой турбиной фирмы «Siemens» V 64.3А. Котел выполняется с естественной циркуляцией, однобарабанным, в газоплотном исполнении для работы под наддувом. Котел снабжен установкой производства собственного конденсата для регулирования температуры пара впрыском. На пусковых режимах впрыск осуществляется питательной водой. Барабан оборудован устройством обогрева и расхолаживания. Котел имеет П-образную компоновку поверхностей нагрева и состоит из топочной камеры и опускного газохода. Топочная камера призматическая с 8-ю вихревыми горелками, расположенными в два яруса на боковых стенах топки. Над горелками установлены сопла для подвода части окислителя из общего короба перед горелками. В опускном газоходе расположены последовательно по ходу газов три ступени конвективного подогревателя, водяной экономайзер, совмещенный с газоводоподогревателем высокого давления (ГВПВД), а также газоводоподогреватель низкого давления (ГВПНД). ГВПНД предназначен для подогрева конденсата паровой турбины перед деаэратором 0,7 МПа при номинальном режиме.

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

Турбина паровая теплофикационная Т-130/160-12,8 предназначена для работы в составе парогазового энергоблока с газовой турбиной фирмы «Siemens» V 64.3А. Турбина сопрягается с генератором переменного тока ТЗФП-160-2УЗ ОАО «Электросила», с воздушным охлаждением, монтируемым на общем с турбиной фундаменте. Турбина представляет собой одновальный двухцилиндровый агрегат, состоящий из однопоточного цилиндра высокого давления и двухпоточного цилиндра низкого давления. подогрев основного конденсата осуществляется последовательно в охладителях основных эжекторов, конденсаторе пара лабиринтовых уплотнений, газоводяном подогревателе низкого давления, деаэраторе. Подогреватели низкого давления №№ 1 и 2 конструктивно выполнены в одном корпусе с соответствующим подогревателем сетевой воды. В режиме работы ПГУ ПНД байпасируются. Турбина имеет, два теплофикационных отбора пара — верхний и нижний, предназначенные для подогрева сетевой воды. Система маслоснабжения турбины питает маслом подшипники турбины и генератора, систему автоматического регулирования турбины и систему гидроподъема роторов. Конденсатор турбины 130КП-9000-1 предназначен для работы на пресной воде и имеет общую поверхность охлаждения 9000м с номинальным расходом охлаждающей воды 16000м3. Конденсатор двухкорпусной с подключение корпусов к разным выхлопам ЦНД. Конденсатор имеет по два встроенных устройства для приема редуцированного пара от БРОУ и горячей воды. С турбиной поставляется устройство шариковой очистки трубок конденсатора и фильтры предочистки.

Фирма «Siemens» проводит фундаментальные исследования, разработку компонентов и процессов, направленные на совершенствование всей концепции электростанции. Эти работы позволяют ожидать повышения КПД паровых электростанций до 50%. Для газовых турбин в ближайшие годы будет превышен показатель 40%, для ПГУ — 60%. повышение КПД приводит к уменьшению расхода топлива и снижению выбросов вредных веществ.

Ставшие возможными вследствие новых разработок увеличение температуры на входе в турбину и повышение КПД компонентов привели к существенному увеличению мощности ГТУ и их КПД. В настоящий момент в странах с частотой тока в сети 50 Гц Фирмаодновременно с развитием газовых турбин фирма «Siemens» разрабатывала и поставляла парогазовые установки с возрастающим КПД. Высокоэффективные ГТУ для рынка с частотой тока 50 Гц номинальной мощностью от 67 до 265 МВт обеспечивают широкий диапазон мощностей парогазовых электростанций одновальных и многовальных схем (даже с отборами тепла). обычная станция комбинированного цикла в настоящий момент может иметь КПД (нетто) до 57,3%. среди электростанций, находящихся в эксплуатации в настоящее время, одновальным блоком ПГУ с газовой турбиной V94.3A в Otahuhu (новая Зеландия) уже достигнут КПД (нетто) 58%.

Таблица 2 — Газовые турбины фирмы «Siemens»

МаркаМощность, МВтКПД, %Расход уходящих газов, кг/сТемпература уходящих газов,°СV64.3A6734,7191589V94.215734,4509537V94.2A19035,2520585V94.3A26538,5656584

Газовая турбина фирмы «Siemens» V 64.3А представляет собой одновальную машину стационарной конструкции, заключенную в один общий корпус. генератор ТЗФГ-63-2УЗ ОАО «Электросила» присоединяется к газовой турбине со стороны компрессора через редуктор.

Герметичный внешний корпус компрессора и турбины состоит из четырех частей:

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

корпус компрессора, служащего держателем направляющих аппаратов компрессора;

центрального корпуса с кольцом камеры сгорания;

корпуса турбины.

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

Сжигание топлива происходит в кольцевой камере сгорания с гибридными горелками, которые обеспечивают низкий уровень выбросов окислов азота при сжигании газового топлива [5]. Кольцевая камера сгорания гарантирует высокую равномерность поля температур горячих газов перед лопаточным аппаратом турбины. конструкция системы сжигания отличается простотой и высокой степенью надежности. горячий газ, поступающий из камеры сгорания, подвергается расширению в четырехступенчатой турбине. Лопатки ротора турбины изготовлены из высокотемпературных сплавов, выдерживающих большие механические и тепловые нагрузки. Ротор объединяет секции компрессора и газовой турбины на одном валу. Газы от газовой турбины направляются к котлу по осевому диффузору.

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

Все направляющие и рабочие лопатки компрессора и турбины допускают возможность своей индивидуальной замены непосредственно на месте установки турбины. Турбина снабжена валоповоротным механизмом, предназначенным для поворота ротора газовой турбины после её останова для предотвращения неравномерного охлаждения, которое может привести к деформации вала. Газотурбинная установка укомплектована маслосистемой, топливной системой, системой пожаротушения, оборудованием промывки ГТУ (поставка фирмы «Siemens»).

1.3 Выбор объекта моделирования и постановка задачи

В качестве объекта исследования в данном дипломном проекте выбрана камера сгорания газовой турбины промышленного типа, на примере агрегата фирмы «Siemens» модели V64.3A [20]. Данная ГТУ относится к классу средних по мощности промышленных турбин с непрерывной кольцевой камерой сгорания.

Камеры сгорания в ГТУ играют главенствующую роль. Остановимся подробнее на её основных функциях и конструктивных особенностях.

В камере сгорания ГТУ осуществляется нагрев рабочего тела до заданной температуры за счет сжигания топлива в потоке сжатого воздуха [5].

К камерам сгорания предъявляются следующие требования:

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

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

для увеличения срока службы они должны иметь надежное охлаждение, особенно наиболее нагретых частей;

высокая экономичность на всех режимах работы ГТУ;

возможно меньшее гидравлическое сопротивление;

надежный запуск;

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

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

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

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

На рисунке 1.2 (а) показана схема простейшей камеры сгорания — прямой цилиндрический канал, соединяющий компрессор с турбиной. К сожалению, такое простое устройство непригодно из-за недопустимо больших потерь давления.

Рисунок 1.2 — Стадии развития схемы традиционной камеры сгорания газотурбинного двигателя

Потери давления пропорциональны квадрату скорости воздушного потока. поскольку скорость воздуха на выходе из компрессора близка к 150 (м/с), потери давления при этом могут достигать четвертой части общего повышения давления в компрессоре. Для снижения потерь давления до приемлемого уровня используют, как показано на рисунке 1.2 (б) диффузор, с помощью которого скорость воздуха уменьшают приблизительно в 5 раз.

Однако этого недостаточно, так как для предотвращения срыва пламени и поддержания устойчивого процесса горения необходимо с помощью обратных токов создать зону малых скоростей. На рисунке 1.2 (в) показано, как этого можно достичь посредством простой пластины. Такое устройство имеет, однако, один недостаток, который заключается в том, что необходимое для получения заданной величины повышения температуры отношение топливо-воздух существенно превышает предел воспламеняемости смесей углеводородов с воздухом. В идеальном случае коэффициент избытка воздуха a близок к 1,25, хотя, например, при желании снизить выбросы окислов азота, эта величина может быть увеличена до 1,6. указанный недостаток может быть устранен, если простой стабилизатор заменить, как показано на рисунке 1.2 (г), перфорированной жаровой трубой. В жаровой трубе создается зона малых скоростей, в которой процесс горения поддерживается циркуляционным потоком продуктов сгорания, непрерывно поджигающим поступающую в камеру свежую топливовоздушную смесь. избыточная (ненужная для горения) часть воздуха вводится в жаровую трубу за зоной горения, где она перемешивается с горячими продуктами сгорания, понижая, таким образом, их температуру до приемлемого для турбины уровня.

Существующие камеры сгорания можно разделить на следующие основные типы:

а)индивидуальные;

б)секционные (многотрубчатые);

в)кольцевые;

г)трубчато-кольцевые.

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

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

существуют две разновидности индивидуальных камер сгорания: цилиндрические и угловые. В цилиндрической камере сгорания, изображенной на рисунке 1.3, воздух разделяется на два потока: первичный и вторичный. Первичный воздух поступает через воздухо-направляющее устройство 1 в пламенную трубу 4, куда через форсунку 2 (или горелку) подается топливо. Расход первичного воздуха регулируется в зависимости от расхода топлива поворотом лопаток воздухо-направляющего устройства 1, что осуществляется посредством специальных рычагов управления. Вторичный (охлаждающий) воздух пропускается через кольцевое пространство между пламенной трубой 4 и корпусом 3 камеры сгорания. При движении он интенсивно охлаждает стенки трубы и корпуса. Выходя из кольцевого пространства, вторичный воздух попадает в объем А, где он смешивается с продуктами сгорания, понижая тем самым их температуру до заданного значения.

Рисунок 1.3 — Схема цилиндрической камеры сгорания

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

В цилиндрических камерах можно установить не одну, а несколько форсунок, что увеличивает надёжность работы и позволяет регулировать тепловую мощность камеры сгорания изменением числа работающих форсунок. Объемная теплонапряженность этих камер составляет 20000-30000 кВт/м3 при давлении 0,4-0,45 МПа, а тепловая мощность камеры сгорания достигает 3000 кДж/ч, расход воздуха — 2,5·105 м3/ч,

К преимуществам индивидуальных цилиндрических камер сгорания относятся простота конструкции и сравнительно малые потери давления, достигающие 1,5-3,0 %. основными недостатками этих камер являются большие масса и габариты.

Секционные (многотрубчатые) камеры сгорания представляют собой конструкцию, в которой объединено несколько (6-16) параллельно работающих цилиндрических камер (секций), часто связанных между собой пламяпередающими патрубками.

Секция многотрубчатой камеры сгорания, в соответствии с рисунком 1.4, состоит из пламенной трубы и кожуха 8. пламенная труба включает в себя головку, состоящую из лопаточного завихрителя 3, тарелки 2 и конуса 4, и корпус, состоящий из цилиндрической части 5 и двух конических участков, соединенных между собой конусным кольцом 6.

Первичный воздух поступает через входной кожух 1 в головку пламенной трубы. Часть его направляется в зону горения через лопаточный завихритель 3, а оставшаяся часть идет туда через многочисленные отверстия в тарелке 2 и конусе 4. кроме того, на цилиндрической части пламенной трубы 5 имеется еще два ряда отверстий, через которые дополнительно поступает воздух, необходимый для горения при полной нагрузке ГТУ. Вторичный воздух идет по кольцевому пространству между пламенной трубой и кожухом 8 и затем поступает в зону смешения через четыре ряда отверстий в конической части пламенной трубы 7. Наибольшая часть охлаждаемого воздуха входит внутрь пламенной трубы через большое число отверстий малого диаметра в конусном кольце 6.

Рисунок 1.4 — Секция многотрубчатой камеры сгорания

Секционные камеры сгорания выполняют обычно в виде единого моноблока, в котором все секции заключены в общий корпус. каждая секция имеет одну форсунку, впрыскивающую топливо по направлению потока. Секционные камеры сгорания отличаются компактностью, обеспечивают высокую полноту сгорания топлива и устойчиво работают в различных эксплуатационных условиях. недостатком их является сравнительно большие потери давления (2,5-7,5%). Тепловая мощность отдельной секции составляет в среднем (0,7-1,7) ·103 кВт, а иногда достигает 3,5·103 кВт. Объемная теплонапряженность у камер этого типа высокая — (100-160) ·103 кВт/м3.

В кольцевых камерах сгорания, в соответствии с рисунком 1.5, зона горения I имеет форму кольцевой полости обычно шириной 150-200 м, которая образуется цилиндрами 1 и 2. Два других соосно-расположенных цилиндра (9 и 8) составляют кожух камеры. Первичный воздух через воздухопроводящее устройство 4 поступает в зону горения I. Вторичный воздух направляется по кольцевым зазорам 6 и 7 к смесительным насадкам 5, через которые поступает в зону II, где смешивается с продуктами сгорания, понижая тем самым их температуру. В воздухоподводящем устройстве 4, на входе в зону горения I по всей окружности расположены форсунки 3. За счет этого обеспечивается хорошее перемешивание топлива с воздухом и горение по всему кольцевому пространству. Число форсунок может достигать 10-20, но иногда это бывает одна вращающаяся форсунка.

рисунок 1.5 — Схема кольцевой камеры сгорания

объемная теплонапряженность у кольцевых камер примерно такая же, как и у секционных, но потери давления несколько больше (до 10 %). По сравнению с секционными камерами они имеют меньший рабочий объем и более равномерное поле температур газа на выходе. Зато кольцевые камеры сложнее в изготовлении и доводке, труднодоступны для осмотра в ходе эксплуатации.

Трубчато-кольцевая камера сгорания представляет собой конструктивное совмещение элементов секционной и кольцевой камер. Так же, как и у кольцевой камеры, кожух её образуется наружным и внутренним соосно-расположенными цилиндрами. А в кольцевом пространстве между этими цилиндрами размещается ряд отдельных пламенных труб, снабженных форсунками. Трубы соединяются друг с другом пламяпередающими патрубками, которые предназначены для передачи пламени, зажигания и выравнивания давления между трубами. Трубчато-кольцевые камеры имеют теплонапряженность и потери давления приблизительно такие же, как секционные камеры. Они компактнее кольцевых камер и более просты в доводке. Небольшие размеры пламенных труб упрощают их изготовление и разборку.

Для работы на жидком топливе в камерах сгорания обычно применяют центробежные форсунки, изображенные на рисунке 1.6. Они просты по конструкции, надежны в работе и обеспечивают хорошее распыливание топлива. К форсунке топливо подаётся насосом 5 под давлением не менее 1,0-1,5 МПа. Поступает оно сначала в кольцевую полость 1, а затем через ряд тангенциально расположенных каналов 2 направляется в вихревую камеру 3, в которой приобретает вращательно-поступательное движение. При выходе из форсунки топливо распыляется под действием центробежных сил.

В центробежных форсунках регулировать расход топлива за счет изменения его давления можно не более чем в 2-2,5 раза. Для обеспечения более широкого диапазона регулирования применяют двухступенчатые форсунки и форсунки с перепуском топлива. У двухступенчатых (двухконтурных) форсунок на малых расходах работает лишь одна первая ступень. Для увеличения расхода топлива к ней подключается вторая ступень. У форсунок с перепуском топлива вихревая камера 3 соединена с регулируемым клапаном 4, который перепускает часть топлива обратно в подводящий трубопровод или же в расходный бак 6.

рисунок 1.6 — Центробежная форсунка с перепуском топлива

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

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

На сегодняшний день время подготовки производства и себестоимость продукции пристально изучаются с целью их сокращения для поддержания конкурентоспособности изделий и лидерства предприятий на рынке, поскольку требования, предъявляемые к предприятиям, постоянно возрастают. Наиболее короткий цикл подготовки производства в сочетании с низкими производственными затратами может быть реализован путем сокращения традиционных методов моделирования и испытаний. С внедрением методов компьютерного моделирования и анализа исключается необходимость повторного создания опытных образцов и проведения их испытаний вследствие корректировки конструкции по результатам предыдущих испытаний, так как корректировка конструкции осуществляется после компьютерного анализа, а испытание опытных образцов носит верификационный характер, как показано на рисунке 1.7 [15].

Рисунок 1.7 — жизненные циклы изделий, основанные на традиционных методах моделирования и компьютерном анализе

Для камер сгорания основными исследуемыми параметрами являются:

а)геометрия камеры сгорания (акустический анализ);

б)тепловой поток (термодинамический анализ);

в)режимы работы горелочного устройства (сопряженный анализ).

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

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

.4 основные показатели работы камер сгорания

Тепловая мощность камеры (кВт) выражается количеством тепла, которое выделяется в единицу времени при полном сгорания топлива [5]:

, (1.1)

где

В — расход сжигаемого топлива, кг/с;

— низшая теплота сгорания топлива, кДж/кг.

объемная теплонапряженность q (кВт/м3) характеризует компактность, а значит, и эффективность использования объема:

, (1.2)

где — объем камеры сгорания, м3 (принимается обычно равным объему пламенной (жаровой) трубы).

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

, (1.3)

где — давление на входе в камеру, МПа.

Потери энергия в камере сгорания состоят из тепловых потерь и потерь давления.

Тепловой к. п. д. камеры сгорания учитывает все тепловые потери:

, (1.4)

где — потери теплоты от неполного сгорания топлива (химический и физический недожог). У современных камер сгорания эти потери не должны превышать 1-5 % общего расхода теплоты при работе на всем диапазоне рабочих нагрузок и 1-3 % при работе на расчетной нагрузке;

— потери за счет отдачи теплоты, в окружающее пространство нагретой поверхностью камеры и примыкающих к ней трубопроводов. Эти потери обычно бывают не более 0,5 % расхода теплоты.

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

Полные потери давления в камере сгорания складываются из гидравлических потерь и дополнительных.

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

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

Полные потери давления принято выражать в долях или процентах по отношению к давлению полного торможения воздуха на входе в камеру:

, (1.5)

где — полная потеря давления в камере;

— давление полного торможения газов на выходе из камеры.

Давление полного торможения воздуха в газе:

,, (1.6)

где — статическое давление воздуха на входе и газа на выходе, МПа; rВ rГ — плотность воздуха и газа, кг/м3; wВ wГ — средние скорости воздуха во входном сечении и газа в выходном сечения камеры, м/с.

В современных конструкциях камер сгорания полные потери давления обычно бывают в пределах = 1-3 %.

Потери давления в камере сгорания снижают к. п. д. ГТУ. Это можно учесть с помощью аэродинамического к. п. д. камеры сгорания , который обычно составляет: . Общий к. п. д. камеры сгорания можно выразить в виде произведения:

. (1.7)

У современных камер сгорания .

Общий коэффициент избытка воздуха в камере сгорания:

, (1.8)

где и — соответственно расход воздуха и топлива в камере, кг/с;

— теоретически необходимое количество воздуха для сжигания 1 кг топлива рассчитывается элементарному составу топлива.

Как уже отмечалось в современных ГТУ a = 4-8.

непосредственно в процессе сгорания участвует не все количество воздуха GB, а только часть его — первичный воздух, поступающий в зону горения:

, (1.9)

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

Обычно .

2. Программная среда разработки проекта

.1 краткое описание метода конечных элементов

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

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

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

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

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

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

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

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

В настоящее время существует достаточно много программных продуктов для решения отдельных классов задач, основанных на методе конечных элементов [13]. Можно подобрать программный продукт практически для любой задачи. Следует отметить, что многие коммерческие программы чрезвычайно дороги (речь идет о десятках тысяч евро), но в отличие от дешевых и бесплатных программных пакетов, они способны представить более высокое качество и скорость решения задач.

один из наиболее мощных коммерческих программных продуктов — это ANSYS [13,15]. ANSYS известен на рынке уже более двадцати лет и является наиболее распространенным средством для научных и инженерных расчетов. особенностью ANSYS является чрезвычайно широкий спектр задач, которые он в состоянии решать. Сюда входят задачи: расчетов на прочность (как линейные, так и нелинейные), теплообмена, гидродинамики, смешанные и даже акустики.

Многоцелевой конечно-элементный пакет инженерного анализа ANSYS позволяет производить расчет в таких дисциплинах как прочность, теплофизика, динамика жидкостей и ряд других. программа ANSYS представляет собой, компьютерную систему для проектирования и выполнения связанного междисциплинарного анализа методом конечных элементов. Программный комплекс ANSYS сертифицирован согласно серии стандартов ISO 9000, 9001, а также бессрочно аттестован в Госатомнадзоре РФ на проведение расчетов для атомной промышленности. ANSYS обеспечивает двустороннюю связь практически со всеми известными CAD-системами [2].

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

За рубежом для расчета с помощью метода конечных элементов широко применяются такие программные пакеты, как NASTRAN, ANSYS, ABAQUS, SAP, ADINA, COSMOS.

.2 основные принципы проектирования в ANSYS

Программный комплекс ANSYS представляет собой многоцелевой пакет для решения сложных проблем физики и механики.

Программа располагает широким перечнем расчетных средств, которые могут учесть разнообразные конструктивные нелинейности. Многоцелевые функции комплекса ANSYS обеспечиваются наличием в нем многочисленного семейства отдельных специализированных программ, имеющих много общих функций, однако математическое обеспечение которых рассчитано на решение отдельных классов задач [13]./Multiphysics — программа для широкого круга инженерных дисциплин, которая позволяет проводить расчеты в области прочности, распространения тепла, механики жидкостей и газов, электромагнетизма, а также решать связанные задачи./Mechanical — программа для решения сложных задач прочности конструкций, теплопередачи, акустики. Эта программа позволяет определять перемещения, напряжения, усилия, температуры давления и другие параметры, важные для оценки механического поведения материалов и прочности конструкции. Данная программа является подмножеством ANSYS Multiphysics./Structural — выполняет сложный прочностной анализ конструкций с учетом разнообразных нелинейностей, среди которых геометрическая и физическая нелинейности, нелинейное точного моделирования поведения больших и сложных расчетных моделей. Данная программа является подмножеством ANSYS/Mechanical./Thermal — это отдельная программа, выделенная из пакем ANSYS/Mechanical, для решения тепловых стационарных и нестационарных задач./LS-DYNA — программа, предназначенная для решения прочностных задач динамики при больших нелинейностях. Эта программа может использоваться для численного моделирования процессов формообразования материалов, анализа аварийных столкновений и ударов при конечных деформациях, включая пробивание, нелинейное программу, обладающую возможностями ANSYS/Multiphysics, но имеющую ограничения по размерам расчетной модели. Эта программа предназначена, в основном, для учебных целей.

предварительная подготовка и вход в программу

Прежде, чем входить в программу, необходимо выполнить следующее:

установить разрешение дисплея не хуже 1024 х 765 пиксел;

цветовая палитра должна включать не менее 256 цветов.

Вход в интерактивный режим осуществляется в следующей последовательности:

настройка параметров;

запуск сессии Run Interactive Now.

после входа в интерактивный режим диалог пользователя с программой осуществляется через многооконный «Графический интерфейс пользователя (GUI)», изображенный на рисунке 2.1, в котором меню представляет собой перечень нескольких логически связанных процедур управления программой, расположенных в отдельных окнах. Кратко опишем некоторые из окон GUI.

Рисунок 2.1 — Графический интерфейс пользователя программы ANSYS

Верхнее горизонтальное окно 1 представляет собой «Меню утилит (Utility menu)». Меню утилит содержит набор часто используемых процедур, которые отображены здесь для доступа в любой момент работы программы.

Окно 2 представляет собой «Окно ввода (ANSYS input)» — область для набора команд и вывод сообщений в «Output Window». Имеется возможность обратиться к списку введенных ранее команд. Команды можно извлекать из файла регистрации (log-файла) введенных ранее команд и (или) входных файлов для последующего ввода.

Окно 3 представляет собой «главное меню (ANSYS Main Menu)». Оно содержит основные функции и этапы выполнения программы, которые группируются в располагаемые сбоку всплывающие (динамические) меню, вид которых зависит от продвижения по программе.

Окно 4 является «Графическим окном». Оно представляет собой область для вывода такой графической информации, как конечно-элементная модель или графики результатов анализа.

Окно 5 — это «Линейка инструментов (Toolbar)». Линейка инструментов позволяет пользователю создавать кнопки и иметь быстрый доступ к часто исполняемым командам.

Окно 6 не представлено на рисунке 2.1 Это «Окно вывода (Output Window)», предназначенное для показа текстовых сообщений программы.

Основные стадии решения задач

Решение задач с помощью программы ANSYS состоит из трех этапов: препроцессорная (предварительная) подготовка (Preprocessing), получение решения (Solving the Equations) и постпроцессорная обработка результатов (Postprocessing). Приведем краткое содержание основных шагов при выполнении каждого из этапов. Более подробное описание с примерами приведено ниже [10].

На стадии препроцессорной подготовки выполняется выбор типа расчета, построение модели и приложение нагрузок (включая и граничные условия). Здесь задаются необходимые для решения исходные данные. Пользователь выбирает координатные системы и типы конечных элементов, указывает упругие постоянные и физико-механические свойства материала, строит твердотельную модель и сетку конечных элементов, выполняет необходимые действия с узлами и элементами сетки, задает уравнения связи и ограничения. Можно также использовать модуль статистического учета для оценки ожидаемых размеров файлов и затрат ресурсов памяти.

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

Исходные данные, введенные при препроцессорной подготовке, становятся частью центральной базы данных программы. Эта база данных разделена на таблицы координатных систем, типов элементов, свойств материала, ключевых точек, узлов сетки, нагрузок и т.д. Как только в таблице появляются некоторые данные, на них становится возможным ссылаться по входному номеру таблицы. например, могут быть определены несколько координатных систем, которые активизируются простой ссылкой на соответствующий номер системы (входной номер таблицы). Кроме того, существует набор команд управления базой данных, чтобы выделить некоторую ее часть для определенных операций. Выделение необходимых данных можно проводить по местоположению геометрических объектов, графическим примитивам твердой модели, типам конечных элементов, видам материалов, номерам узлов и элементов и т.п. Так, например, сложные граничные условия можно легко указать или изменить, используя геометрическое элементов.

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

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

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

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

В программе ANSYS доступны следующие два способа моделирования: нисходящий и восходящий. В первом случае пользователь указываем только самый высокий порядок сложности объектов модели. Используемые обычно объекты (такие, как сферы и призмы, т.е. формы, которые называются геометрическими примитивами) могут быть созданы за одно обращение к меню. Например, пользователь определяет объемный примитив, а программа автоматически находит связанные с ним поверхности, линии и ключевые точки. Примитивы позволяют непосредственно указывать геометрические формы. В программе ANSYS можно легко и быстро определить в двумерном случае такие формы, как окружности и прямоугольники, или параллелепипеды, сферы, конусы и цилиндры в трехмерном. после того как геометрические объекты указаны (с помощью примитивов, считыванием данных из файлов формата. IGES или непосредственным построением), к ним можно применять операции булевой алгебры. При импортировании геометрии в формате. IGES пользователь имеет возможность управлять значениями допусков на слияние объектов модели, выявлять «проблемные» области и возможные ошибки.

В случае непосредственного создания модели в интерактивном режиме работы чаще всего применяется так называемое «восходящее моделирование». При восходящем моделировании пользователь строит модель, начиная с объектов самого низкого порядка. сначала задаются ключевые точки, затем связанные с ними линии, поверхности и объёмы — именно в таком порядке.

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

Еще одним эффективным методом построения модели в программе ANSYS является построение некоторой поверхности с помощью так называемого метода «обтягивания каркаса». С помощью этого метода можно задать некоторый набор поперечных сечений, а затем дать программе команду построить поверхность, которая будет точно соответствовать указанным сечениям.

После того, как построена модель, строится ее конечно-элементный аналог (т.е. сетка узлов и элементов).

библиотека конечных элементов программы ANSYS содержит более 80 типов, каждый из которых определяет, среди прочего, применимость элемента к той или иной области расчетов (прочностной, тепловой, магнитный и электрический анализы, движение жидкости или связанные задачи), характерную форму элемента (линейную, плоскую, в виде бруска и т.д.), а также двухмерность (2-D) или трехмерность (3-D) элемента.

После выбора типа элементов необходимо задать их константы. Константы элемента — это свойства, специфичные для данного типа элемента. например, для элемента ВЕАМЗ — балочного 2-D элемента — константами являются площадь поперечного сечения, момент инерции, высота и др.

Свойства материала требуются для большинства типов элементов. В зависимости от области приложения свойства могут быть линейными, нелинейными и (или) анизотропными.

Линейные свойства могут зависеть или не зависеть от температуры, быть изотропными или ортотропными. Зависимость свойств от температуры имеет форму полинома (вплоть до четвертой степени) или задается таблично.

Нелинейные соотношения, такие как кривые деформирования, кривые намагничивания материала, кривые ползучести, обычно задаются в виде таблицы.

В программе ANSYS предусмотрено четыре способа генерации сетки: использование метода экструзии, создание упорядоченной сетки, создание произвольной сетки (автоматически) и адаптивное построение.

Метод экструзии (выдавливания) используется для превращения областей двумерной сетки в трехмерные объекты, состоящие из параллелепипедов, клиновидных элементов или их комбинации. Процесс экструзии осуществляется с помощью процедур смещения из плоскости, буксировки, поступательного и вращательного перемещений.

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

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

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

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

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

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

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

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

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

Приложение нагрузок и получение решения

после того, как при препроцессорной подготовке построена расчетная модель, можно переходить к стадии решения задачи. Этот этап включает в себя задание вида анализа и его опций, нагрузок, шага решения и заканчивается запуском на счет конечно-элементной задачи [10,12].

Программа ANSYS предусматривает два метода решения задач: h-метод и р-метод. Первый из перечисленных методов может применяться при любом типе расчетов (статический, динамический, тепловой и т.п.), в то время как второй метод может использоваться только в линейном статическом анализе. При прочих равных условиях, h-метод требует более частой сетки, чем р — метод.

Под нагрузками понимаются как внешние и внутренние усилия, так и граничные условия в виде ограничений на перемещения. В программе ANSYS нагрузки разделены на следующие категории:

ограничения степеней свободы;

сосредоточенные силы;

поверхностные нагрузки;

объемные силы;

инерционные нагрузки.

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

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

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

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

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

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

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

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

Для того чтобы получить решение за минимальное время, программа ANSYS переупорядочивает расположение элементов и узлов.

Постпроцессорная обработка

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

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

Итогом работы программы на постпроцессорной стадии является графическое и (или) табличное представление результатов. Графическое изображение может быть выведено на монитор в интерактивном режиме во время постпроцессорной обработки или преобразования в твердую копию.

На стадии получения решения результаты записываются в базу данных программы ANSYS и в так называемый «файл результатов». Результаты, полученные на каждом дополнительном шаге решения, накапливаются как наборы данных.

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

Типы основных файлов, создаваемых программой

В процессе моделирования программа ANSYS может создавать следующие типы файлов [10]:

а). log — протокольный (журнальный) файл, хранящий историю работы в виде ANSYSкоманд, т.е. полная запись программы (текстовый файл) на всех этапах решения задачи (препроцессор, решение, постпроцессор);

б).inp или .dаt — файл ввода программы при пакетном (batch) режиме работы;

в).db или .dbb — база данных программы (бинарный, в кодировке ANSYS), сохраняет геометрию модели, граничные условия, решение;

г).emat — файл матрицы элементов;

д).grph — файл графики;

е).err — файл ошибок, содержит все ошибки и предупреждения, выданные программой в процессе решения;

ж).mac — файл с макрокомандами;

з).nod — файл с узлами сетки;

и).elm — файл с элементами сетки;

к).rst — бинарный файл, содержащий результаты прочностного анализа;

л).rth — бинарный файл, содержащий результаты теплового анализа;

м).rmg — бинарный файл, содержащий результаты электромагнитного анализа;

н).rfl — бинарный файл, содержащий результаты гидродинамического анализа;

о).SI_MPL — файл — библиотека, в котором хранятся заданные пользователем свойства материала.

3. Проектирование модели исследуемого объекта

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

Программа FLOTRAN и модуль FLOTRAN CFD (расчетная гидрогазодинамика) программы ANSYS предлагают всесторонние инструментальные средства для анализа двумерных и трехмерных полей потока жидкости или газа.

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

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

проектирование графической модели объекта;

выбор вида конечно-элементной сетки и её наложение на графическую модель;

определение граничных условий и нагрузок;

определение вида и параметров решения;

проведение решения и вывод результатов;

анализ полученных результатов решения.

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

.1 Построение графической модели

Графическая модель камеры сгорания была разработана в прикладном графическом пакете SolidWorks 2004. Она представляет собой двумерную конструкцию линий, являющихся границами плоскостей радиального сечения форсунки горелочного устройства и рабочей области камеры сгорания. Для того чтобы иметь возможность применять модель в пакете ANSYS, она была сохранена в форме IGES — объекта.

после импортирования файла модели она будет выглядеть, как показано на рисунке 3.1 необходимо доработать графическую модель средствами пакета. Применение графических возможностей препроцессора ANSYS делает построенную модель более «дружелюбной» к другим компонентам инженерного пакета, в том числе, исключает возможные проблемы при импортировании плоскостей из других графических процессоров (SolidWorks, 3D Max Studio, AutoCAD и других).

Доработка включает в себя:

масштабирование импортированных линий, для соответствия реальных размеров камеры сгорания, размерам графической модели, выраженных в единицах СИ;

создание плоскостей, необходимых для последующего нанесения конечно-элементной сетки;

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

Создаются две плоскости, в соответствии с рисунком 3.2: одна из них статическая, не участвует в анализе, представляет собой разрез форсунки горелочного устройства; другая — рабочая плоскость, радиальный разрез внутренней полости кольцевой камеры сгорания. Именно на рабочую область будет в последующем наложена конечно-элементная сетка.

Рисунок 3.1 — Импортированная SolidWorks в ANSYS графическая модель камеры сгорания

рисунок 3.2 — Плоскости, созданные средствами препроцессора ANSYS

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

.2 Наложение конечно-элементной сетки

Библиотека элементов ANSYS содержит более 150 различных типов элементов [10,12]. каждый тип элемента имеет уникальный номер и имя, которое идентифицирует категорию элемента: BEAM4, PLANE77, SOLID96, и так далее. Как правило, в различных вариантах программы доступны не все виды элементов. Это зависит от того, какие виды расчета могут использоваться в данном варианте программы.

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

Для двумерного газодинамического анализа во FLOTRAN используем элементы типа FLUID141.

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

Для элементов FLOTRAN CFD скорости определяются в соответствии с законом сохранения импульса, давление определяется в соответствии с законом сохранения массы, температура определяется в соответствии с законом сохранения энергии.

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

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

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

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

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

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

каждая из областей должна быть разбита на соответствующее сложности и требуемой точности количество элементов.

Для того, чтобы плоскость была разбита на элементы, сначала нужно разбить линии, являющиеся внешними границами областей. У первой области линии разбиваются с расстоянием между узлами 0.03 (м), у второй — 0.04 (м) и у третьей — 0.03 (м). Выбор расстояния между узлами должен быть основан на необходимой точности решения, а также времени выполнения расчетов. далее осуществляется выбор формы элементов и вида разбиения. В данном случае применим треугольную форму, по причинам, обозначенным выше, и выберем свободное расположение элементов.

Указав нужную нам (рабочую) плоскость, в соответствии с рисунком 3.3, получаем область с наложенной конечно-элементной сеткой.

Рисунок 3.3 — Нанесенная на рабочую площадь конечно-элементная сетка

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

.3 Задание начальных условий и нагрузок

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

К граничным условиям относится задание стенок камеры сгорания. Скорость узлов у них приравнена к нулю (Vx = 0 и Vy = 0). Допущением является то, что модель не учитывает излучение, принимаемое стенками от пламени. Вместо этого, на них наложена термальная нагрузка в виде функциональной зависимости от значения нагрузок на входящие потоки воздуха и топливного газа.

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

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

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

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

К сожалению, версия пакета, в которой проводилась разработка дипломного проекта, не обладает возможностью моделирования реакции горения [8]. Из современных инженерных пакетов серии ANSYS такими возможностями обладает модуль CFX, который введен в состав полнофункциональной версии ANSYS 10.0. Исходя из этого, добиться необходимой точности решения, можно только имея в наличии аналитический тепловой расчет камеры сгорания. В рамках данного дипломного проекта предполагается создать эмуляцию пламени фиксированного размера и с варьируемой температурой. кроме того, чтобы добиться приемлемых результатов и для расчета переноса вещества, было решено создать поле обратных скоростей. Поле обратных скоростей имитирует участие в реакции горения кислорода атмосферного воздуха, который подается компрессором в камеру сгорания. параметры этих имитаторов находятся в зависимости от параметров входящих газов. Зависимость определена опытным путем.

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

На этом построение конечно-элементной модели камеры сгорания газовой турбины можно считать завершенным. результат представлен на рисунке 3.4 Красным цветом показаны нагрузки типа скорость (Velocity), желтым — температурные (Temperature), зеленым — давления (Pressure).

Рисунок 3.4 — Конечно-элементная модель камеры сгорания газовой турбины V64.3A

далее следует выбор стратегии решения и выполнение программой расчетных итераций.

.4 Подготовка и проведение расчетов

Расчет термодинамики газов внутри камеры сгорания необходимо проводить в несколько этапов. Этап в инженерном пакете ANSYS именуется шагом решения.

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

На первом шаге компоненты переноса не выделяются. Уравнения переноса решаются с воздухом в системе СИ. Для решения на данном шаге необходимо небольшое количество итераций и в данном случае оно равно 50. Командный код этого шага представлен в листинге 1.

Листинг 1 — первый шаг решения

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

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

Таблица 3 — основные свойства компонентов газовой смеси

ПараметрТопливный газВоздух из компрессораНазваниеgasvozМолекулярная масса16.0428.96Плотность при нормальных условиях0.71681.2928Варьирование плотностидадаВязкость при нормальных условиях1.786E-0051.2067E-005Варьирование вязкостидадаТеплопроводность0.025980.02674Коэффициент массовой диффузии2.601E-0052.149E-005

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

после задания компонентных нагрузок все настройки многокомпонентного решения сделаны. Далее необходимо указать число итераций, в нашем случае их 100. На этом очередной шаг решения завершен. Его командный код представлен в листинге 2.

Листинг 2 — второй шаг решения

На следующем этапе необходимо решить уравнение энергии для получения распределения температур по рабочей области. Для этого отключаем расчет переноса вещества и включаем опцию решения уравнения энергии.

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

Листинг 3 — третий шаг решения

! Подготовка к решению уравнений энергии,solu,temp,t! Включить решатель уравнения энергии вещества,solu,flow,f! заморозить поле потока,meth,temp,3! Активировать решатель PCCR,conv,temp,1. e-8! критерий сходимости для PCCR, iter,exec,50! Для решения необходимо несколько итераций,relx,temp,1.0! Отсутствует следующим шагом будет совместное решение уравнений энергии и переноса вещества. При этом включается опция расчета потоков и устанавливается необходимое число итераций. Для первичной обработки нужно около 50 итераций. Полученный результат для улучшения сходимости решения обрабатываем еще раз с числом итераций 100. Командный код этого шага представлен в листинге 4.

после уведомления об успешном решении, предварительный расчет камеры сгорания завершен. Данные расчетов можно просмотреть в интерактивном режиме в основном постпроцессоре (General Post). Но для получения более точного с физической точки зрения решения, необходимо учесть то, что поток в камере сгорания является сжимаемым. Это можно определить и из геометрии объекта, которая имеет вид сопла, и из того, что воздух компрессора нагнетается под давлением, поступая с очень большой скоростью.

Листинг 4 — Четвертый шаг решения

О том, что поток сжимаемый, можно судить и по выходной скорости отработавших газов. Если она близка к скорости звука, то поток в сопле считается сжимаемым. Скорость звука в газе равна примерно 514 (м/с). В нашем расчете скорость выходного потока составляет примерно 600-750 (м/с). Следовательно, поток можно считать сжимаемым. Для решения уравнений переноса в этом случае ANSYS позволяет включить дополнительную опцию расчета сжимаемого потока.

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

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

Листинг 5 — пятый, шестой и седьмой шаги решения

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

.5 Определение оптимального коэффициента избытка воздуха

В последнее десятилетие были введены новые нормативные требования к уровню выбросов вредных веществ наземными газотурбинными установками (ГТУ) [2]. Самое существенное ограничение касается уровня выбросов оксидов азота, который не должен превышать 50 мг/нм3 при 15 % O2. Применение концепции предварительной подготовки топливовоздушной смеси позволяет снизить эмиссию оксидов азота до уровня, требуемого нормативными документами. однако при этом возникает проблема обеспечения устойчивости камеры сгорания к пульсациям давления, возникающим вследствие случайных турбулентных возмущений или из-за изменения режима работы камеры сгорания [4].

Ведущие зарубежные фирмы, такие как Rolls-Royce, АBB, Siemens, решают данную проблему либо за счет системы управления, поддерживающей постоянное соотношение между расходами воздуха и топлива, поступающими в камеру сгорания, либо используя системы подавления колебаний давления, включающих в себя исполнительные механизмы, содержащие подвижные части в топливном и воздушном трактах камеры сгорания. Применение первого метода уменьшает КПД ГТУ, а использование систем подавления колебаний значительно увеличивает стоимость камеры сгорания. кроме того, применение подвижных частей, работающих с высокой частотой, снижает надежность камеры сгорания. поэтому более перспективным путем решения данной проблемы является применение систем пассивного подавления автоколебаний, то есть при создании камер сгорания данного типа более предпочтительным является применение пассивных систем подавления колебаний давления, в которых не предполагается высокочастотного управления работой камеры сгорания.

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

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

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

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

При проектировании ГТУ коэффициент принимается согласно установленным нормам, в условиях эксплуатации — устанавливаются экспериментально. Исходя из опыта применения газовых турбин в промышленности, a принимает значение от 4 до 8.

Разработанная методика решения задачи термодинамики газов в камере сгорания турбины позволяет произвести оценку оптимального значения коэффициента избытка воздуха. Это достигается путем проведения ряда экспериментов с варьированием входных параметров воздуха компрессора и топливного газа.

Было проведено пять экспериментов для a = 2, 4, 6, 8,10. Результаты экспериментов представлены в приложении Б.

Камера сгорания этого класса имеет следующие рабочие показатели [19]:

температура на выхлопе сопла Твых = 1200-1300 (°С);

скорость потока на выходе сопла Vвых = 750 (м/с).

На диаграммах шкала температур имеет размерность Кельвин, а скорость задана в условных единицах. одна условная единица 1у. е. = 33 (м/с)

Из приведенных диаграмм можно сделать вывод, что стабильная работа турбины будет при a = 6. снижение a приводит к чрезмерному выбросу топлива, а увеличение — к увеличению температуры выхлопа и возможным пульсациям пламени.

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

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

Выводы по главе

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

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

4. безопасность и экологичность проекта

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

Установленная электрическая мощность — 420 МВт, установленная тепловая мощность — 1463 Гкал/час.

Основное топливо — природный газ, поступающий с Уренгойского месторождения.

резервное топливо — мазут (резервная подача мазута рассчитана на 3 дня эксплуатации станции).

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

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

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

.1 Характеристика условий труда

Автоматизированное рабочее место инженера I категории представляет собой операторский ПЭВМ, обеспечивающий работу специалиста в операторной энергоблока. Условия труда описаны в таблице 4.

Таблица 4 — характеристика условий труда

ПоказательЗначениеОбъем помещения, м36х10х3,5=210характеристика тяжести работыСредняяПлощадь помещения, м260Тип системы вентиляцииобщеобменная

Фактическое состояние условий труда оператора по степени вредности и опасности на рабочем месте представлены в таблице 5.

Таблица 5 — Степень вредности и опасности на рабочем месте

№НаименованиеЕд. изм. Фактическое значениеДопустимое значениеФизические факторы1ШумдБА58-60602Температура воздуха°С24,022,0-24,03Скорость движения воздухам/с<0.10.14Естественное освещениеКЕО, %1,51,55Относительная влажность%4040-60%6Освещённостьлк260-520300-500химические факторы1Азбестосодержащая пыльмг/м3522Серная кислотамг/м3Не выявлен13Аммиакмг/м3Не выявлен20

Также имеют место и другие вредные производственные факторы:

ионизирующее излучение;

высокая напряженность электромагнитного поля;

повышенный уровень вибрации.

Также в зоне работы могут иметь место опасные производственные факторы:

вращающиеся и движущиеся машины, и механизмы, подвижные части оборудования;

повышенная температура поверхностей оборудования;

расположение на территории производства сосудов высокого давления;

напряжение в электроустановках и электрической цепи;

расположение устройств на высоте более 1,3м относительно поверхности пола, земли;

газоопасные места;

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

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

при повышенном уровне шума нужно применять противошумные наушники, вкладыши «Беруши»;

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

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

при выполнении работ на участках с температурой воздуха выше 32°С необходимо применять режим труда с интервалами времени для отдыха и охлаждения;

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

при нахождении в помещениях с действующим технологическим оборудованием (за исключением щитов управления ТГ) необходимо носить защитную каску для защиты головы от ударов случайными предметами;

при работе на высоте более 1,3 м над уровнем пола (земли) площадки — предохранительный пояс и страховочный канат или работать с лесов;

при недостаточной освещенности рабочей зоны следует применять дополнительное местное освещение (фонари, переносные электрические светильники);

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

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

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

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

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

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

Оператору должны выдавать бесплатно согласно отраслевым нормам следующие средства индивидуальной защиты:

костюм хлопчатобумажный — на 12 мес.;

каска защитная — на 2 года;

очки защитные — до износа;

«беруши» — до износа;

куртка х/б на утепляющей прокладке — на 2 года;

ботинки кожаные — на 12 мес.

.2 Расчет искусственного освещения

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

Освещение производственных помещений характеризуется количественными и качественными показателями.

К количественным показателям относятся: лучистая энергия и лучистый поток, световой поток, сила света, яркость и освещенность.

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

При подборе источников света и размещении рабочих мест в помещении следует ограничивать прямую блесткость от источников освещения, при этом яркость светящихся поверхностей (окна, светильники и др.), находящихся в поле зрения, не должна быть более 200 Кд/м2.

Для освещения помещений с мониторами и ПЭВМ рекомендуется применять светильники серии ЛПО36 с зеркализованными решетками, укомплектованные высокочастотными пускорегулирующими аппаратами (ВЧ ПРА). Применение светильников без рассеивателей и экранирующих решеток не допускается.

Яркость светильников общего освещения в зоне углов излучения от 50° до 90° с вертикалью в продольной и поперечной плоскостях должна составлять не более 200 Кд/м2, защитный угол светильников должен быть не менее 40°. светильники местного освещения должны иметь не просвечивающий отражатель с защитным углом не менее 40 градусов.

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

Источники, виды и системы искусственного освещения

Для искусственного освещения производственных помещений обычно применяют газоразрядные лампы (люминесцентные низкого давления, дуговые ртутные, металлогалогенные и ксеноновые — высокого давления).

Газоразрядные лампы имеют высокую светоотдачу (до 100 лм/Вт) и большой срок службы. Серийно выпускаются люминесцентные лампы низкого давления, различающиеся спектральным составом света:

ЛД (дневного света);

ЛДЦ (дневного света с улучшенной цветопередачей);

ЛБ (белого света);

ЛТБ (теплого белого света);

ЛХБ (холодного белого света).

Расчет искусственного освещения с помощью коэффициента использования светового потока

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

Для расчета лучше подходит определение освещенности с помощью коэффициента использования светового потока, так как в данном случае соблюдается два условия: равномерное освещение горизонтальных поверхностей и отсутствие крупных затеняющих предметов [1].

Основная формула для расчетов имеет вид:

, (4.1)

где E — нормируемая освещенность, лк, при классе и подклассе зрительной работы 4В (E=200 лк);- коэффициент запаса, для газоразрядных светильников (k3=1.5);- освещаемая площадь, м2 (S = 60 м2);- коэффициент неравномерности освещения (для люминесцентных ламп z=1.1);- число ламп в светильнике (n=2);- число светильников (N=6);

h — коэффициент использования светового потока.

Находим индекс помещения:

, (4.2)

где A,B,h — соответственно длина, ширина, высота подвеса светильников над рабочей поверхностью.

В соответствии с формулой (4.2) получаем индекс помещения I = 1,5.

Коэффициент отражения потолка rп принимаем 70 %.

Коэффициент отражения стен rст принимаем 60 %.

В соответствии с полученными значениями выбираем коэффициент использования светового потока (для светильников, снабженных люминесцентными лампами) h=58 %.

Подставляя данные в формулу (4.1) получаем

.

Исходя из полученного результата, можно сделать вывод, что в данном случае подходят люминесцентные лампы ЛБ 40.

основные характеристики ламп следующие:

мощность 40 Вт;

световой поток 3000 лм;

световая отдача 75 лм/Вт.

точное число светильников рассчитывается по формуле:

. (4.3)

Подставляя точное число светильников:

.

Следовательно, для данного помещения достаточно 6 светильников ЛБ 40 (по 2 лампы в светильнике).

Схема расположения светильников приведена на рисунке 4.1.

рисунок 4.1 — Схема расположения светильников

Светотехнические расчеты являются основополагающими при проектировании осветительных установок.

.3 Оценка экологичности проекта

выхлоп газовой турбины, когда энергоблок работает по схеме ПГУ, поступает в топку энергетического котла.

основными контролируемыми параметрами выхлопа являются в частности концентрация оксида азота.

Расчет выбросов загрязняющих веществ при сжигании топлива производиться по формуле [11]:

, (4.4)

где B — расход топлива, (т/год);

— теплота сгорания натурального топлива, (МДж/кг);

— коэффициент, зависящий от мощности котлоагрегата;

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

Вышеперечисленные коэффициенты для природного газа: (МДж/кг), , .

Расход газа на котел в год .

Рассчитаем массу выбросов в соответствии с формулой (4.4):

Концентрация равна = 0.077 ().

Предельно допустимая концентрация равна 0.085 ().

Плата за выброс рассчитывается по следующей формуле:

, (4.5)

где — масса выброса в тоннах;

— норматив в рублях за тонну;

— коэффициент экологичности;

— коэффициент индексации.

необходимые параметры имеют следующие значения: , , , .

Подставляя эти значения в формулу (4.5) получаем расходы за выброс загрязняющих веществ:

.

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

.4 Оценка чрезвычайных ситуаций

Рассматриваемый теплоэнергетический объект связан с высокой экологической опасностью, а в частности:

выброс вредных и химически активных веществ в окружающую среду в случае аварии;

взрывы и разрушение прилегающих к производственной территории объектов и строений;

в случае масштабной аварии возникает реальная опасность прекращения подачи тепла в теплосети, также перебои с подачей электроэнергии;

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

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

первый энергоблок ПГУ-190/220;

турбинное отделение;

котельное отделение.

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

краткая характеристика пожарной безопасности этого здания приведена в таблице 6.

Таблица 6 — Оценка пожаробезопасности основного здания

Наименование объектаКатегория зданияКраткая характеристика зданияСтепень огнестойкости зданияСгораемые материалыКласс пожараЦехГЖелезобетонное здание со стальным каркасом, крышей из листовой профилированной стали. Котлы, турбины, Электроустановки высокой мощности. Электрические кабели. IМазут, Газ — пропан, Изоляционная лентаЕ

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

Углекислый газ является бесцветным, не имеющим запаха, электрически не проводящим, инертным газом, пригодным для тушения пожара. Для того чтобы гасить пожар посредством СО2, необходимо увеличить концентрацию по объёму до 34% и более. Требуемое количество СО2 передаётся на защищённый участок по системе трубопроводов. специальные форсунки для СО2 распространяют СО2 (внутри трубопроводов имеет жидкую форму) по защищённому участку. Выпуск СО2 может создать опасность для персонала в виде кислородного голодания и пониженной видимости.

системы пожаротушения на СО2 используются для защиты дорогостоящего оборудования, которое может быть повреждено системой пожаротушения на воде. Газ СО2 из системы противопожарной защиты выпускается внутри герметизированного корпуса с оборудованием с интенсивностью, создающей расчётную концентрацию (изначально 37%), чтобы обеспечить полное погашение огня.

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

Автоматическое управление системы пожаротушения на СО2 действует следующим образом.

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

включение звукового аварийного сигнала (сирена с первичным предупреждающим звуковым сигналом), мигающие лампы внутри корпуса и панели с подсветкой на входах;

общая звуковая и индивидуальная визуальная аварийная сигнализация на щите пожарной сигнализации.

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

включение звуковой предупреждающей сигнализации выпуска газа (вторичный звуковой сигнал или звук сирены);

включение общей звуковой и индивидуальной визуальной аварийной сигнализации на щите пожарной сигнализации.

Срабатывание временной задержки выпуска газа позволяет персоналу в течение приблизительно 30 секунд покинуть защищённый участок (помещение). При ручном пуске времени меньше, около 10 секунд.

Кнопки аварийного прерывания предусмотрены на каждом входе на участки, защищённые СО2. Данные кнопки остановят запуск системы только в том случае, если они нажаты до того, как сработал датчик на второй линии. Если процесс пожаротушения начат, его нельзя остановить. Только постоянное нажатие кнопки аварийного прерывания предотвратит приведение системы в действие.

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

Выводы по разделу

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

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

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

5. кономическая эффективность проекта

5.1 Источники экономической эффективности

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

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

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

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

Так же нужно учитывать экономию средств при доводке и эксплуатации объекта моделирования на предприятии заказчике.

5.2 Расчет единовременных затрат

В общем случае единовременные затраты предприятия изготовителя определяются как затраты на разработку и программирование задачи, и представлены в следующей формуле [7]:

, (5.1)

где Краз — затраты на разработку системы, руб.;

Кпрог — затраты на программирование, руб.

перечень исходных данных предприятия разработчика для расчета единовременных затрат представлен в таблице 7.

Таблица 7 — Данные для расчета единовременных затрат предприятия разработчика

№ПоказательЗначение1Заработная плата разработчика, руб. 150002Коэффициент доплат к заработной плате, доли ед. 0,23Районный коэффициент, доли ед. 0,154Коэффициент отчисления в социальные фонды, доли ед. 0,265Время разработки системы, месяцы36время использования ЭВМ для разработки программы, час4507Коэффициент накладных расходов, доли ед. 0,048Годовой фонд работы ЭВМ, час18009Годовой фонд оплаты труда персонала обслуживающего ЭВМ, р. 7200010Норма амортизационных отчислений ЭВМ, доли ед. 0,1211Норма амортизационных отчислений здания, доли ед. 0,0212Площадь занимаемая ЭВМ, м2613Стоимость одного м2 здания, руб. 550014Стоимость ЭВМ, руб. 3500015Коэффициент накладных расходов на экспл. ЭВМ, доли ед. 0,0416Потребляемая мощность ЭВМ, Вт35017Стоимость кВт/часа, руб. 0,9218Коэффициент затрат на ремонт ЭВМ (от стоимости), доли ед. 0,0519Затраты на материалы для эксплуатации ЭВМ, руб. 114020Коэффициент затрат на содержание ЭВМ, доли ед. 0,1521Ставка НДС, доли ед. 0,1822Коэффициент перевода единиц времени16823Коэффициент интенсивного использования мощности ЭВМ0,75

Произведем расчет затрат на постановку задачи моделирования, представленные в следующем виде:

, (5.2)

где

Зо — месячный оклад разработчика, р.;

Траз — трудоемкость разработки проекта и проектной документации, ч×мес.;

Кд, Кр — соответственно коэффициенты доплат к заработной плате и районный, доли ед.;

Ксн — коэффициент отчислений на социальные нужды, доли ед.;

Кн. раз — коэффициент накладных расходов, доли ед.

Перечень элементов трудоемкости разработки представлен в таблице 8.

Таблица 8 — трудоемкость отдельных видов работ

Стадии разработкиТрудоемкость, чел. месяц1. изучение патентов0.22. Изучение литературных источников0.53. Разработка технического задания0.34. Разработка эскизного проекта0.25. Разработка технического проекта0.36. Разработка рабочего проекта17. Внедрение проекта0.3

Рассчитаем Краз по формуле (5.2):

далее произведем расчет затрат на программирование модели и её решение. затраты на разработку программной модели и её решение можно рассчитать следующим образом:

, (5.3)

где Тпрог — время на создание программы, мес.;

Кн. прогр — коэффициент накладных расходов, доли ед.;мч — стоимость машино-часа ЭВМ, р.;

Кч — коэффициент перевода единиц времени.

Стоимость машино-часа ЭВМ рассчитывается по формуле:

, (5.4)

где Sэкс — годовые эксплуатационные расходы, связанные с обслуживанием ЭВМ, р.;

Тпол — годовой фонд работы ЭВМ, час.

Эксплуатационные расходы рассчитываются по формуле:

экс = 12×ЗП × (1 + Кд) × (1 + Кр) × (1 + Ксн) + А + Тр + Э + М + Нрэкс, (5.5)

где ЗП — месячная оплата труда всего обслуживающего персонала в сумме, р.;

А — амортизационные отчисления от стоимости ЭВМ и здания, р. /год;

Тр — затраты на ремонт, р/год;

Э — затраты на электроэнергию, р/год;

М — затраты на материалы, р.;

Нрэкс — накладные расходы, связанные с эксплуатацией ЭВМ, р/год.

затраты на амортизацию вычисляются по формуле:

А = Кэвм×Нэвм + Сзд×Sзд×Нзд, (5.6)

где Кэвм — балансовая стоимость ЭВМ, р.;

Нэвм — норма амортизационных отчислений от стоимости ЭВМ, доли ед.;

Сзд — стоимость 1 м2 здания, р/м2;зд — площадь, занимаемая ЭВМ, м2;

Нзд — норма амортизационных отчислений от стоимости здания, доли ед.

Затраты на ремонт вычислим по формуле:

Тр = Кэвм × Ктрэвм, (5.7)

где

Ктрэвм — коэффициент, учитывающий затраты на ремонт ЭВМ.

Затраты на электроэнергию, потребляемую ЭВМ за год эксплуатации определяем по формуле:

Э = Ц×Тпол×N×Км, (5.8)

где Ц — цена за один кВт×ч электроэнергии, р.; N — потребляемая мощность, кВт; Км — коэффициент интенсивного использования мощности вычислительной техники.

В годовые эксплуатационные затраты по обслуживанию ЭВМ входят также накладные расходы, которые рассчитываются по формуле:

Нрэкс = 12×Зо× (1 + Кд) × (1 + Кр) ×Кнэкс, (5.9)

где Кнэкс — коэффициент накладных расходов, связанных с эксплуатацией ЭВМ.

Затраты на материалы определяем по формуле:

, (5.10)

где i — вид материала; Цi — цена i-того материала, р.; Мi — количество i-го материала.

Расчет может быть представлен в виде таблице 9.

Таблица 9 — перечень и стоимость материалов используемых ЭВМ

Наименование материалаКоличество в годЦена за ед., р. Стоимость, р. Гибкие магнитные диски, штук61590Тонер2150300Бумага, пачек. 4150600Салфетки обтирочные1015150Итого1140

Подставив данные из таблицы 7 в (5.6) получаем затраты на амортизацию (А):

А = 35000×0,12 + 5500×6×0,02 = 4860 (руб).

Подставив данные из таблицы 7 в (5.7) получаем затраты на ремонт:

Тр = 35000×0,05 = 1750 (руб).

Подставив данные из таблицы 7 в (5.8) получаем затраты на электроэнергию:

Э = 0,92×1800×0,350×0,75 = 434,7 (руб).

Подставив данные из таблицы 7 в (5.9) получаем накладные расходы:

Нрэкс = 12×15000× (1 + 0,2) × (1 + 0,15) ×0,04 = 9936 (руб).

Расчет затрат на материалы представлен в виде таблицы 9: М = 1140 (руб). Подставив результаты расчетов (5.6), (5.7), (5.8), (5.9), (5.10) в (5.5) получим эксплуатационные расходы Sэкс:

экс = 72000× (1 + 0,2) × (1 + 0,15) × (1 + 0,26) + 4860 + 1750 + 434,7 + 1140

+ 9936 = 134518,30 (руб).

Подставив данные из таблицы 7 и результат вычислений (5.5) в формулу (5.4) получим стоимость одного машино-часа (Смч):

Смч = 134518,30/1800 = 74,73 (руб).

Исходя из полученных результатов для (5.4) и исходных данных таблицы 7 находим капитальные затраты на разработку программной модели и её решение (Кпрог):

Кпрог = 15000×1,3× (1 + 0,2) × (1 + 0,15) × (1 + 0,26) × (1 + 0,04) +

+ 74,73 ×1,3×168 = 45706,75 (руб).

В итоге единовременные затраты составляют:

К = Краз + Кпрог = 40687,92 + 45706,75 = 86394,67 (руб).

5.3 Оценка экономической целесообразности и выводы по разделу

Проектирование какого либо промышленного агрегата связано, прежде всего, с аналитическими расчетами и испытаниями. Для камеры сгорания ГТУ основными значимыми параметрами, которые влияют на стоимость проектирования и изготовления, являются:

определение оптимальной геометрии стенок камеры сгорания;

оптимальный выбор материала стенок камеры сгорания;

оптимальный выбор горелочного устройства.

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

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

Определение геометрии стенок камеры весьма трудоемкий процесс. От неё, так или иначе, зависят такие показатели работы ГТУ как экологичность, КПД и надежность агрегата в целом. методика проведения газодинамического расчета, предложенная в данном проекте позволяет, меняя геометрию, проследить за акустическими свойствами, тепловым распределением и другими параметрами камер сгорания.

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

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

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

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

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

Заключение

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

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

Был проведен ряд экспериментов на примере ГТУ производства фирмы «Siemens» серии V64.3A [20], которые показали состоятельность данной методики при анализе оптимальности коэффициента избытка воздуха. Точность модели небольшая, но и она позволяет сделать оценку критических режимов работы по основным показателям работы камеры сгорания, это скорость потока на выходе из камеры, температура на выходе и общее состояние пламени (по турбулентной диаграмме).

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

список использованных источников

1Безопасность жизнедеятельности и промышленная безопасность: Учебное пособие / Под ред. проф. В.Д. Шантарина — Тюмень: ТюмГНГУ, 2001. — 308с.

2Вестник» Самарского государственного аэрокосмического университета имени академика С.П. Королева, №2 — Материалы IV Всероссийской научно-практической конференции «Процессы горения, теплообмена и экология тепловых двигателей». Самара: СГАУ, 2002. — 111с.

Зенкевич О. метод конечных элементов в технике: Пер. с англ. — М.: мир, 1975. — 541с., ил.

Кашапов Р.С., Максимов Д.А., Скиба Д.В., Куликов С.В., Баштанников М.Н. исследование автоколебаний давления в камере сгорания с предварительным смешением топлива // Газотурбинные технологии №4 (13). Рыбинск, 2001. С.34-37.

Латыпов Р.Ш. вопросы рациональной эксплуатации газотурбинных установок: Учебное пособие. Уфа: УГНТУ, 2000. — 100с.

Методические указания по дипломному проектированию для студентов очного и заочного обучения специальности 220201 — Управление и информатика в технических системах. — Тюмень: ТюмГНГУ, 2006. — 28с.

Методические указания к оценке экономической эффективности автоматизированных систем в курсовом и дипломном проектировании для студентов специальности АСОиУ, АТП, УИТС дневного и заочного обучения — Тюмень: ТюмГНГУ, 2003. — 33с.

«Моделирование процессов горения в ANSYS CFX. Хитрых Д., «ANSYS Solutions. русская редакция» № 3 (май) ‘2006.

Норри Д., де Фриз Ж. Введение в метод конечных элементов: Пер. с англ. — М.: мир, 1981. — 304с., ил.

Справочная система CAE-пакета ANSYS.

Старикова Г.В., Столбова Н.В., Дорофеева Э.С. Охрана окружающей природной среды: Учебное пособие. — Тюмень: ТюмГНГУ, 2000. — 94с.

Чигарев А.В., Кравчук А.С., Смалюк А.Ф. ANSYS для инженеров: Справ. пособие. М.: Машиностроение-1, 2004, 512с.

13www.ansys. spb.ru <#"center">Приложения

приложение А

Листинг программы

!***************************************************************************

!Программа газодинамического расчета камеры сгорания газовой турбины v64.3A *

!автор: студент 5-го курса УИТС, Панышев Евгений Александрович *

!***************************************************************************

!Начало

!Импортирование модели формата IGES в ANSYS

/AUX15,IGES,NODEFEAT,MERGE,YES,SOLID,YES,SMALL,YES,GTOLER, DEFA,’chambre’,’IGS’,’.GrafModel’

/REPLOT,RESIZE

/AUTO, 1

/REP

/VIEW, 1 ,,,1

/ANG, 1

/REP,FAST

/REPLOT,RESIZE

!Выбор единиц измерения

/units,si0020

!Задание констант моделирования

/PREP7=1.54 !Скорость воздуха=450 !температура воздуха=1.57 !Скорость газа=350 !температура газа=(vvoz*tvoz+vgas*tgas)/vgas-vvoz*vvoz*15=tflame*0.6 !температура стенок(т.к. не учитываем излучения пламени)

!Сжатие номеров,KP,LINE

!Постройка площадей,2,26,4,2,21 !Выбор линий первой (рабочей) плоскости,2,20,2,19,2,18,2,17,2,10,2,9,2,8,2,28,2,4,2,1,2,3,2,27,2,5,2,6,2,7,2,11,2,12,2,13,2,14,2,15,2,24,2,23,2,25,2,26,2,22,P51X !Постройка рабочей плоскости,2,6,4,2,24 !Выбор линий второй плоскости,2,23,2,25,2,16,2,26,2,22,P51X !Постройка второй плоскости

!Подключение фильтров графического меню

/NOPR

/PMETH,OFF,0,PR_SET,1,PR_STRUC,0,PR_THERM,1,PR_FLUID,1,PR_ELMAG,0,MAGNOD,0,MAGEDG,0,MAGHFE,0,MAGELC,0,PR_MULTI,2,PR_CFD,1

/GO

!Выбор типа элементов,1,FLUID141

!!!!!!!!Разбитие на элементы рабочей площади

!Разметка линий,2,26,4,ORDE,5,2,1,2,3,2,-15,2,17,2,-28,P51X, , ,-1, ,1,5,15,4,ORDE,4,5,11,5,-15,5,17,5,-26,S, , ,P51X,5,15,4,ORDE,4,5,11,5,-15,5,17,5,-26,_Y,LINE, , , ,P51X,_Y1,LINE,,_Y,_Y1,0.03, , , , , , ,1 !Разбивка линий по расстоянию между узлами (0.03),5,6,4,ORDE,2,5,5,5,-10,S, , ,P51X,5,6,4,ORDE,2,5,5,5,-10,_Y,LINE, , , ,P51X,_Y1,LINE,,_Y,_Y1,0.04, , , , , , ,1 !Разбивка линий по расстоянию между узлами (0.04),5,5,4,ORDE,5,5,1,5,3,5,-4,5,27,5,-28,S, , ,P51X,5,5,4,ORDE,5,5,1,5,3,5,-4,5,27,5,-28,_Y,LINE, , , ,P51X,_Y1,LINE,,_Y,_Y1,0.03, , , , , , ,1 !Разбивка линий по расстоянию между узлами (0.03)

!Выбор вида элементов (треугольные),1,2D,0,_Y,AREA, , , , 1 !Выбор площади для разбиения,_Y1,AREA,S,_Y,_Y1 !Наложение конечно-элементной сетки на выбранную площадь,_Y,_Y1,_Y2

!!!!!!!!!!!!!!!!Задание начальных условий и нагружений

! Задание параметров стенок КС

/UI,MESH,OFF,ALL,5,22,4,ORDE,6 !Выбор линий,5,3,5,-14,5,18,5,-23,5,25,5,-28,S, , ,P51X,ALL,S,1,5,239,1,ORDE,10,5,1,5,-100,5,112,5,-218,5,224,5,-233,5,240,5,-249,5,255,5,-266,_Y,NODE,R, , ,P51X,_Y1,NODE,S,_Y,_Y,_Y1,VX,0 !Задние скорости на границы КС,_Y1,VY,0,_Y1,temp,tsten !Задание температуры стенок КС,_Y1

! Задание скорости течения воздуха на левый патрубок,S, , , 17,ALL,S,1,5,7,1,ORDE,3,5,245,5,249,5,-254,_Y,NODE,R, , ,P51X,_Y1,NODE,S,_Y,_Y,_Y1,VX,vvoz!Задание нагрузки,_Y1,VY,0,_Y1

!Задание температуры потока на левый патрубок,2,7,1,ORDE,3,2,245,2,249,2,-254

/GO,P51X,TEMP,tvoz !Задание нагрузки

!Задание скорости потока на правый патрубок,S, , , 15 !Выбор линии,ALL!Выбор всех узлов на линии,S,1,5,7,1,ORDE,3,5,212,5,218,5,-223,_Y,NODE,R, , ,P51X,_Y1,NODE,S,_Y,_Y,_Y1,VX,vvoz*sin(0.2618) !Задание нагрузки,_Y1,VY,-vvoz*sin(1.309)!в соответствии с глобальной системой координат,_Y1

!Задание темепературы потока на правый патрубок,2,7,1,ORDE,3,2,212,2,218,2,-223

/GO,P51X,TEMP,tvoz !Задание значения нагрузки

!Задание скорости потока газа,S, , , 24 !Выбор линий,ALL!Выбор всех узлов линии,S,1,5,8,1,ORDE,3,5,228,5,233,5,-239,_Y,NODE,R, , ,P51X,_Y1,NODE,S,_Y,_Y,_Y1,VX,vgas*sin(0.8727) !Задание значения нагрузки,_Y1,VY,-vgas*sin(0.6981)!в соответствии с глобальной системой координат,_Y1

!Задание температуры потока газа,2,8,1,ORDE,3,2,228,2,233,2,-239

/GO,P51X,TEMP,tgas !Задание значения нагрузки

!Задание граничного условия на выход КС,S, , , 1 !Выбор линии,ALL!Выбол всех узлов линии,S,1,5,13,1,ORDE,3,5,91,5,100,5,-111,_Y,NODE,R, , ,P51X,_Y1,NODE,S,_Y,_Y

/GO,_Y1,PRES,0!Задание значения нагрузки,_Y1!Сохранение конечно-элементной модели с нагружениями

!Создание поля температур для эмуляции пламени,5,1032,1,ORDE,84 !Выбор элементов поля,5,443,5,453,5,459,5,468,5,-470,5,475,5,480,5,697,5,700,5,-703,5,705,5,-706,5,708,5,721,5,725,5,-732,5,925,5,928,5,-931,5,933,5,946,5,949,5,-958,5,1142,5,-1150,5,1161,5,-1162,5,1165,5,-1172,5,1341,5,1344,5,-1345,5,1348,5,-1356,5,1365,5,-1367,5,1369,5,-1374,5,1539,5,-1540,5,1543,5,-1550,5,1566,5,-1573,5,1645,5,-1654,5,1725,5,1728,5,-1735,5,1749,5,-1758,5,1813,5,-1836,5,1896,5,-1903,5,1911,5,-1921,5,1943,5,-1944,5,1949,5,-1951,5,1957,5,-1985,5,2039,5,-2045,5,2050,5,2053,5,-2062,5,2075,5,-2117,5,2123,5,-2144,5,2163,5,-2172,5,2175,5,2177,5,-2184,5,2187,5,-2266,5,2268,5,-2269,5,2271,5,2274,5,-2914,S, , ,P51X,2,1032,1,ORDE,84,2,443,2,453,2,459,2,468,2,-470,2,475,2,480,2,697,2,700,2,-703,2,705,2,-706,2,708,2,721,2,725,2,-732,2,925,2,928,2,-931,2,933,2,946,2,949,2,-958,2,1142,2,-1150,2,1161,2,-1162,2,1165,2,-1172,2,1341,2,1344,2,-1345,2,1348,2,-1356,2,1365,2,-1367,2,1369,2,-1374,2,1539,2,-1540,2,1543,2,-1550,2,1566,2,-1573,2,1645,2,-1654,2,1725,2,1728,2,-1735,2,1749,2,-1758,2,1813,2,-1836,2,1896,2,-1903,2,1911,2,-1921,2,1943,2,-1944,2,1949,2,-1951,2,1957,2,-1985,2,2039,2,-2045,2,2050,2,2053,2,-2062,2,2075,2,-2117,2,2123,2,-2144,2,2163,2,-2172,2,2175,2,2177,2,-2184,2,2187,2,-2266,2,2268,2,-2269,2,2271,2,2274,2,-2914

/GO,P51X,TEMP,tflame!Задание нагрузки

!Задние поля обратного потока в пламени,5,255,1,ORDE,106 !Выбор элементов,5,728,5,-730,5,925,5,931,5,946,5,951,5,-956,5,1144,5,-1149,5,1165,5,1167,5,-1170,5,1349,5,-1354,5,1370,5,-1374,5,1545,5,-1550,5,1571,5,-1573,5,1730,5,-1735,5,1755,5,-1758,5,1896,5,-1903,5,1912,5,1919,5,-1921,5,2040,5,-2045,5,2054,5,2060,5,-2062,5,2163,5,2165,5,-2170,5,2178,5,2183,5,-2184,5,2279,5,2281,5,-2285,5,2290,5,2295,5,-2298,5,2387,5,-2390,5,2394,5,-2396,5,2398,5,-2402,5,2482,5,2484,5,-2487,5,2492,5,-2493,5,2496,5,-2499,5,2507,5,-2510,5,2576,5,2578,5,-2582,5,2584,5,-2589,5,2594,5,-2600,5,2624,5,-2627,5,2658,5,-2667,5,2671,5,-2678,5,2698,5,-2702,5,2728,5,-2730,5,2733,5,-2736,5,2739,5,-2747,5,2764,5,-2769,5,2790,5,-2795,5,2798,5,-2806,5,2820,5,-2825,5,2841,5,-2843,5,2846,5,-2854,5,2862,5,-2868,5,2878,5,2880,5,-2881,5,2884,5,-2891,5,2894,5,-2899,5,2906,5,2909,5,-2914,S, , ,P51X,5,255,1,ORDE,106,5,728,5,-730,5,925,5,931,5,946,5,951,5,-956,5,1144,5,-1149,5,1165,5,1167,5,-1170,5,1349,5,-1354,5,1370,5,-1374,5,1545,5,-1550,5,1571,5,-1573,5,1730,5,-1735,5,1755,5,-1758,5,1896,5,-1903,5,1912,5,1919,5,-1921,5,2040,5,-2045,5,2054,5,2060,5,-2062,5,2163,5,2165,5,-2170,5,2178,5,2183,5,-2184,5,2279,5,2281,5,-2285,5,2290,5,2295,5,-2298,5,2387,5,-2390,5,2394,5,-2396,5,2398,5,-2402,5,2482,5,2484,5,-2487,5,2492,5,-2493,5,2496,5,-2499,5,2507,5,-2510,5,2576,5,2578,5,-2582,5,2584,5,-2589,5,2594,5,-2600,5,2624,5,-2627,5,2658,5,-2667,5,2671,5,-2678,5,2698,5,-2702,5,2728,5,-2730,5,2733,5,-2736,5,2739,5,-2747,5,2764,5,-2769,5,2790,5,-2795,5,2798,5,-2806,5,2820,5,-2825,5,2841,5,-2843,5,2846,5,-2854,5,2862,5,-2868,5,2878,5,2880,5,-2881,5,2884,5,-2891,5,2894,5,-2899,5,2906,5,2909,5,-2914,_Y,NODE,R, , ,P51X,_Y1,NODE,S,_Y,_Y,_Y1,VX,-(0.01*exp(tflame/1000)-0.04)*sin(0.8727) !Скорость противотока газов,_Y1,VY,(0.01*exp(tflame/1000)-0.04)*sin(0.6981),_Y1,ENKE,-1,_Y1!Сохранение конечно-элементной модели со всеми нагружениями

!Настройка параметров решателя,SOLU,TRAN,0 !Используем стационарный анализ,SOLU,FLOW,1 !Используем расчет потоков,SOLU,TEMP,0 !Не используем температурный расчет,SOLU,TURB,1 !Используем модель турбулентности (по умолчанию),SOLU,COMP,0,SOLU,VOF,0,SOLU,SFTS,0,SOLU,IVSH,0,SOLU,SWRL,0,SOLU,SPEC,0,SOLU,ALE,0 !В расчете участвуют все элементы

/prep7,conv,iter,5!Задание параметров сходимости,CONV,PRES,1e-008,

!Задание параметров среды

/solu,iter,exec,50!необходимо несколько итераций для начального расчета,temp,nomi,350 !начальная температура 350K,prot,dens,air-si !Начальное решение для воздуха,vary,dens,true,prot,visc,air-si,vary,visc,true,prot,cond,air-si,vary,cond,true,prot,spht,air-si!Провести 50 итераций для воздуха при 350K!Подготовка к многокомпонентному переносу

!Задание параметров многокомпонентного переноса

/prep7,prot,dens,cmix !Растворенная смесь для плотности,prot,visc,cmix !Растворенная смесь для вязкости,1,1,2 !установка двух компонент смеси,solu,spec,t !Включение опции многокомпонентного переноса,1 !Задание компонента 1 (GAS) в качестве алгебраического

!!!!!!! Задание свойств каждого из компонентов

!Для топливного газа,1,gas,16.04!Молекулярная масса,1,DENS,GAS,0.7168,300,1.01325E+5!Плотность,1,dens,t!Варьирование плотности,1,.3!Массовая доля вещества,1,VISC,GAS,3.786E-4,300,1.01325E+5!Вязкость,1,mdif,GAS,1.601E-6!Массовая диффузия,1,visc,t!Варьирование вязкости,1,cond,GAS,.01598!теплопроводность,1,1,2,voz,28.96!Молекулярная масса,2,DENS,GAS,1.2928,300,1.01325E+5!Плотность,2,dens,t!Варьирование плотности,2,.7!Массовая доля вещества,2,VISC,GAS,2.2067E-5!Вязкость,2,mdif,GAS,2.149E-5!Массовая диффузия,2,visc,t!Варьирование вязкости,2,cond,GAS,.02674!теплопроводность,2,1

!Степень стабилизации параметров,1,1,,2,1,

!Настройка решателя многокомпонентного переноса вещества,2,1.0,2,3,2,,,2,1.e-9,1,3,1,,,2,1.e-9

!!!!!!! установка граничных условий смесей,s,,,15,17,,1,all,gas,0.0 !Для подвода воздуха,all,voz,1.0,s,,,24,,1,all,gas,1.0!Для подвода топливного газа,all,voz,0.0

!!!!!!! конец установки граничных условий

/prep7,iter,exec,100 !Необходимо выполнить 100 глобальных итераций,conv,outp,land !Настроить монитор сходимости

/solu!Провести 100 итераций с активированными компонентами

!Подготовка к решению уравнений энергии,solu,temp,t !Включить решатель уравнения энергии вещества,solu,flow,f !заморозить поле потока,meth,temp,3 !Активировать решатель PCCR,conv,temp,1.e-8 !критерий сходимости для PCCR,iter,exec,50 !Для решения необходимо несколько итераций,relx,temp,1.0 !Отсутствует необходимо выполнить 100 дополнительных итераций

!для улучшения результатов расчета!Решение

!!Улучшение модели турбулентности и задание её параметров

/solu,TURB,MODL,4!Выбор модели турбулентности,NKET,C1MX,0.43,!Настройка параметров,NKET,C2,1.9,,NKET,SCTK,1,,NKET,SCTD,1.2,!Решение с новой моделью турбулентности

!!!!!!!Включение опции сжимаемого потока

!Расчет уравнения энергии

/SOLU,SOLU,TRAN,0,SOLU,FLOW,0 !Отключение опции переноса вещества,SOLU,TEMP,1,SOLU,TURB,1,SOLU,COMP,1 !Сжимаемый поток,SOLU,VOF,0,SOLU,SFTS,0,SOLU,IVSH,0,SOLU,SWRL,0,SOLU,SPEC,1,SOLU,ALE,0

!Совместное решение,SOLU,TRAN,0,SOLU,FLOW,1 !Включение опции переноса вещества,SOLU,TEMP,1,SOLU,TURB,1,SOLU,COMP,1 !Сжимаемый поток,SOLU,VOF,0,SOLU,SFTS,0,SOLU,IVSH,0,SOLU,SWRL,0,SOLU,SPEC,1,SOLU,ALE,0

!************************************************************************

!Расчет модели камеры сгорания газовой турбины закончен*

!************************************************************************

!Подготовка к просмотру результатов решения

/POST1,LAST

!Просмотр результатов ведется в интерактивном режиме

!Конец

приложение Б

Диаграммы результатов экспериментов

Рисунок Б.4 — Диаграмма поля распределения температур в камере сгорания газовой турбины V64.3A (a = 8)

рисунок Б.6 — Диаграмма векторного поля распределения скоростей (a = 8)

рисунок Б.7 — Диаграмма поля распределения температур в камере сгорания газовой турбины V64.3A (a = 10)

рисунок Б.8 — Диаграмма кинетической энергии турбулентности (a = 10)

рисунок Б.9 — Диаграмма векторного поля распределения скоростей (a = 10)

рисунок Б.10 — Диаграмма поля распределения температур в камере сгорания газовой турбины V64.3A (a = 4)

рисунок Б.11 — Диаграмма кинетической энергии турбулентности (a = 4)

рисунок Б.12 — Диаграмма векторного поля распределения скоростей (a = 4)

рисунок Б.13 — Диаграмма поля распределения температур в камере сгорания газовой турбины V64.3A (a = 2)

рисунок Б.15 — Диаграмма векторного поля распределения скоростей (a = 2)

Учебная работа. Моделирование тепловой схемы ПГУ