Table Of ContentCopyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
МИНОБРНАУКИ РОССИИ
ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ ОБРАЗОВАТЕЛЬНОЕ
УЧРЕЖДЕНИЕ ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ
«САМАРСКИЙ ГОСУДАРСТВЕННЫЙ АЭРОКОСМИЧЕСКИЙ
УНИВЕРСИТЕТ имени академика С.П. КОРОЛЕВА
(НАЦИОНАЛЬНЫЙ ИССЛЕДОВАТЕЛЬСКИЙ УНИВЕРСИТЕТ)»
О.В. Батурин, Д.А. Колмакова, В.Н. Матвеев,
Г.М. Попов, Л.С. Шаблий
Расчет пространственной структуры потока в ступени
осевого компрессора в программном комплексе ANSYS CFX
Электронное учебное пособие
САМАРА
2011
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
УДКСГАУ: 621.438.001.2
ББК СГАУ: 39.55
Авторы: Батурин Олег Витальевич,
Колмакова Дарья Алексеевна,
Матвеев Валерий Николаевич,
Попов Григорий Михайлович,
Шаблий Леонид Сергеевич
Рецензент: д-р техн. наук, профессор С.В. Фалалеев
Батурин, Расчет пространственной структуры потока в ступени
осевого компрессора в программном комплексе ANSYS CFX [Электронный
ресурс] : электрон, учеб. пособие / О.В. Батурин, Д.А. Колмакова, В.Н. Матвеев,
Г.М. Попов, Л.С. Шаблий; Минобрнауки России, Самар, гос. аэрокосм. ун-т.
С.П. Королева (нац. исслед. ун-т). - Электрон, текстовые и граф. дан.
(3,1 Мбайт). - Самара, 2011. - 1 эл. опт. диск (CD-ROM). - Систем, требования:
ПК Pentium; Windows 98 или выше.
Учебное пособие содержит краткие сведения о расчетном исследовании
рабочего процесса в ступени осевого компрессора с использованием
современного CAD-программного комплекса Ansys CFX, основанного на
решении уравнений Навье - Стокса. Подробно описаны процесс создания
расчетной модели, создания сетки конечных элементов в автоматизированном
режиме, наложение граничных условий, проведения расчета и обработки его
результатов.
Пособие разработано на кафедре теории двигателей летательных аппаратов
СГАУ и предназначено для студентов, обучающихся по курсам «Теория, расчет
и проектирование авиационных двигателей и энергетических установок»,
«Теория и расчет лопаточных машин», «Агрегаты наддува двигателей», а также
для самостоятельной работы студентов-дипломников, аспирантов, магистрантов
и научных работников, чьи исследования связаны с турбомашиностроением.
О Самарский государственный
аэрокосмический университет, 2011
2
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
ОГЛАВЛЕНИЕ
ОСНОВНЫЕ УСЛОВНЫЕ ОБОЗНАЧЕНИЯ 6
ВВЕДЕНИЕ 9
1. ОСНОВНЫЕ ЭТАПЫ СОЗДАНИЯ РАСЧЕТНЫХ МОДЕЛЕЙ
ПОТОКА РАБОЧЕЕО ТЕЛА В ПРОТОЧНОЙ ЧАСТИ ТРЕХМЕРНОЙ
12
СТУПЕНИ ОСЕВОЕО КОМПРЕССОРА В ПРОЕРАММНОМ
КОМПЛЕКСЕ ANSYS CFX
2. ОСНОВНЫЕ СВЕДЕНИЯ О ПОСТРОЕНИИ РАСЧЕТНЫХ
МОДЕЛЕЙ ПОТОКА В МЕЖЛОПАТОЧНЫХ КАНАЛАХ ВЕНЦОВ
14
ОСЕВЫХ КОМПРЕССОРОВ В ПРОЕРАММЕ TurboGrid ПАКЕТА
ANSYSCFX
2.1. Типовой алгоритм построения расчетных моделей потока в
14
межлопаточных каналах в программе TurboGrid
2.2. Независимый запуск программы TurboGrid 16
2.3. Базовые сведения о работе в программе TurboGrid 16
3. ПОДГОТОВКА ИСХОДНЫХ ДАННЫХ И ПОСТРОЕНИЕ
ЕЕОМЕТРИИ МОДЕЛИ ПОТОКА РАБОЧЕЕО ТЕЛА В ПРОТОЧНОЙ 28
ЧАСТИ СТУПЕНИ ОСЕВОЕО КОМПРЕССОРА
3.1. Описание геометрии межлопаточного канала в программа Tur
28
boGrid
3.2. Подготовка файлов с описанием геометрии расчетной области
31
для программы TurboGrid
3.3. Загрузка в программу TurboGrid геометрии межлопаточного ка
33
нала венца осевого компрессора
3.4. Задание величины и геометрии радиального зазора 40
3.5. Корректировка размеров входного и выходного участков 42
4. ПОСТРОЕНИЕ СЕТКИ КОНЕЧНЫХ ЭЛЕМЕНТОВ В СОЗДАННОЙ
45
РАСЧЕТНОЙ ОБЛАСТИ В ПРОЕРАММЕ TurboGrid
4.1. Создание блочной топологии сетки конечных элементов в меж-
45
лопаточном канале
4.2. Настройка параметров сетки конечных элементов 47
4.3. Оценка качества конечноэлементной сетки 51
4.4. Еенерация сетки конечных элементов в расчетной области 55
4.5. Завершение работы с моделью в TurboGrid и передача создан
56
ной расчетной модели потока в программный комплекс Ansys CFX
3
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
5. ПОДГОТОВКА РАСЧЕТНОЙ МОДЕЛИ ТЕЧЕНИЯ ВОЗДУХА В
59
СТУПЕНИ КОМПРЕССОРА В ПРОГРАММЕ ANSYS CFX - PRE
5.1. Запуск программы ANSYS CFX-Pre 59
5.2. Создание проекта в программе ANSYS CFX-Pre 60
5.3. Определение основных параметров турбомашины (Basic Set
61
tings)
5.4. Загрузка компонентов, составляющих ступень компрессора
62
(Component Definition)
5.5. Определение граничных условий турбомашины и настройка ре
65
шателя (Physics Definition)
5.6. Проверка созданных интерфейсов (Interface Defimition)n гранич- ^
ных условий {Boundary Definition)
5.7. Завершение работы с шаблонами Turbomachinery 69
5.8. Редактирование свойств рабочего тела 70
5.9. Запись созданной модели с граничными условиями и передача
ее на решение
6. РАСЧЕТ ТЕЧЕНИЯ В СТУПЕНИ ОСЕВОГО КОМПРЕССОРА В
ПРОГРАММНОМ КОМПЛЕКСЕ С ИСПОЛЬЗОВАНИЕМ 75
ПРОГРАММНОГО КОМПЛЕКСА ANSYS CFX
1. ОБРАБОТКА РЕЗУЛЬТАТОВ РАСЧЕТА ТЕЧЕНИЯ В СТУПЕНИ
ОСЕВОГО КОМПРЕССОРА В ПРОГРАММНОМ КОМПЛЕКСЕ ANSYS 80
CFX
7.1. Запуск программы ANSYS CFX - Post 80
7.2. Загрузка результатов решения 81
7.3. Определение результатов расчета как результатов расчета тур- ^
бомашины
7.4. Отображение на экране нескольких межлопаточных каналов 84
7.5. Постройка графиков изменения параметров потока вдоль по
верхности лопаток осевого компрессора в произвольном по высоте 85
сечении
7.6. Управление отображением расчетной области 87
7.7. Создание турбоповерхностей 88
7.8. Построение картины распределения параметров трехмерного
потока в ступени осевого компрессора на созданной турбоповерхно- 90
сти
7.9. Построение линий тока в ступени осевого компрессора 96
7.10. Построения поля векторов скорости в ступени осевого ком
97
прессора на созданной ранее турбоповерхности
4
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
7.11. Определение интегральных параметров потока в ступени осе- ^
вого компрессора
СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ 100
5
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
ОСНОВНЫЕ УСЛОВНЫЕ ОБОЗНАЧЕНИЯ
высота лопатки, хорда профиля;
скорость потока в абсолютной системе координат, теплоем
кость;
изобарная теплоемкость;
D - диаметр;
F - площадь сечения;
G - массовый расход;
i - энтальпия;
1рю 1на " 3™ атаки;
к - показатель адиабаты;
L - удельная работа;
I - длина;
М - число Маха;
N - мощность;
и - частота вращения ротора;
р - давление;
Р - сила, тяга;
г - радиус;
R - газовая постоянная;
Т - температура;
t - шаг решетки;
t - относительный шаг;
и - окружная скорость;
w - относительная скорость;
число лопаток в ЛВ, число ступеней, число гребешков лаби-
z
ринтного уплотнения;
т(А); тг(А);
газодинамические функции;
е(А); q(A)
угол между фронтом решетки и вектором абсолютной скоро
а
сти;
угол между фронтом решетки и вектором относительной ско
р рости;
88 - величина зазора;
угол отставания потока в абсолютной и относительной систе
8а, 8р
мах координат;
у - угол установки профиля;
jj - коэффициент полезного действия;
6
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Л - приведенная скорость;
ц - коэффициент расхода;
як' - степень повышения давления компрессора;
р - плотность;
рст - степень реактивности ступени;
а - коэффициент восстановления полного давления;
т - угловая скорость вращения ротора, угол заострения профиля;
£ - коэффициент потерь
Индексы
* — параметр торможения;
1 — сечение на входе в рабочее колесо;
2 — сечение на выходе из рабочего колеса;
3 — сечение на выходе из направляющего аппарата;
е — воздух;
рз — относящийся к радиальному зазору;
ср относящийся к среднему сечению;
—
сот — относящийся к ступени;
а — осевое направление;
opt оптимальный;
—
р — изобарный процесс (p=const);
г — радиальное направление;
s — изоэнтропический (идеальный) процесс;
и — окружное направление;
параметр, рассматриваемый в относительном движении, т.е. в
w — системе координат, связанной с вращающимся рабочим коле
сом;
х, у, z — проекции на координатные оси;
А — изменение величины;
Л — связанный с влиянием скорости потока;
Е — суммарная величина.
7
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Сокращения
ВМ — верхнее меню;
ВНА — входной направляющий аппарат;
ГДФ — газодинамические функции;
ГМ — главное меню;
ГТД — газотурбинный двигатель;
ГТУ — газотурбинная установка;
ГУ — граничное условие;
КПД — коэффициент полезного действия;
ЛВ — лопаточный венец;
JIKM — левая кнопка мыши;
НА — направляющий аппарат;
ОС — операционная система;
ПЧ — проточная часть;
ПКМ — правая кнопка мыши;
РЗ — радиальный зазор;
РК — рабочее колесо;
СКМ — средняя кнопка мыши;
СК — система координат;
ТРД — турбореактивный двигатель;
(англ. аббр. Computational Fluid Dynamics) -«вычислительная га
зовая динамика».
Остальные обозначения, индексы и условные сокращения объяснены в тексте.
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
ВВЕДЕНИЕ
Газодинамическое совершенствование узлов газотурбинных
двигателей как авиационного, так и наземного назначения является
основным направлением их развития, поскольку они непосредст
венно влияют на топливную экономичность и конкурентоспособ
ность изделия в целом. Например, повышение КПД турбины высо
кого давления ТРД с умеренными параметрами цикла на 7 % приво
дит к снижению удельного расхода топлива на 0,7%. Аналогичное
влияние оказывают и другие турбомашины на совершенство ГТД.
В настоящее время турбомашины достигли высокого уровня со
вершенства. Современные ГТД характеризуются высокими пара
метрами цикла, что приводит к существенному уменьшению высот
лопаток, а тенденция сокращения числа ступеней требует примене
ния высоконагруженных сверхзвуковых венцов. Поэтому достиже
ние заложенного в техническом задании значения КПД турбома
шины, не говоря уже о его повышении, является сложной научно-
технической задачей. В случае проектирования компрессоров за
дача осложняется еще и требованием обеспечения достаточных за
пасов газодинамической устойчивости, а они, как известно, снижа
ются с ростом параметров цикла. Причем за последние 20...30 лет
время на решение задачи сократилось от десяти до 3...5 лет, что
требует резко сократить число доводочных испытаний и, по сути,
получать конструкцию близкую к окончательной уже «с листа».
Сложившаяся ситуация привела к полному пересмотру про
цесса проектирования и газодинамической доводки узлов ГТД и
ГТУ. В него внедрены новые методики, содержащие наиболее со
вершенные физические модели. Широко применяются методы вы
числительной газовой динамики, основанные на численном реше
нии системы уравнений Навье - Стокса, и способные заменить
большую часть дорогих натурных экспериментов. В настоящий
момент времени решение уравнений Навье-Стокса является самым
9
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
совершенным методом расчета газодинамических процессов. При
менение методов вычислительной газовой динамики позволяет на
этапе первоначального проектирования выявить действительную
картину течения, определить необходимые изменения элементов
ПЧ, ведущие к устранению обнаруженных недостатков, а также от
крывает возможности для поиска новых прогрессивных решений.
Это позволяет проводить значительную часть проектировочных и
доводочных работ на компьютере, сократив число экспериментов
до минимума.
В мировой практике уже однозначно установлено, что переход
к новому процессу проектирования позволяет снизить количество
средств на доводку авиационного (газотурбинного) двигателя бо
лее, чем на порядок.
В данной работе рассматривается последовательность форми
рования виртуальных моделей потока воздуха в ступени осевого
компрессора газотурбинного двигателя для исследования методами
вычислительной газовой динамики. Созданная подобным образом
расчетная модель в дальнейшем может быть использована для ана
лиза рабочего процесса спроектированной ступени компрессора,
выявления ошибок проектирования и поиска путей увеличения ее
газодинамической эффективности.
С/'/) - м о д с л и р о в а н и с любого течения условно можно разделить
на следующие этапы:
- выбор расчетной схемы, граничных условий и основных
допущений;
- построение расчетной модели: включает в себя определение
границ рассматриваемой области, создание расчетной сетки,
назначение граничных условий, задание параметров счета;
- непосредственно расчет;
- визуализация результатов расчета и определение основных
расчетных параметров (интегральных значений скорости,
давления, температуры и т.п.) и анализ выполненных расче-
ю
