Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Численное моделирование колебательных 2FSI-процессов в компрессорах газоперекачивающих агрегатов

..pdf
Скачиваний:
1
Добавлен:
20.11.2023
Размер:
7.11 Mб
Скачать

Рис. 2.6. ГУ на входе, выходе и стенках при учете газодинамики ЛУ

Таблица 2.3

ГУ на входе, выходе и стенках при учете газодинамики ЛУ

 

 

Скоростьвращения

 

Цветовое

Тип ГУ

Параметр

 

, об/мин

 

 

п/п

 

 

 

 

 

 

 

обозначение

 

 

4200

 

5160

 

6000

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

Вход

Полное давление P*, МПа

5,585

 

5,933

 

6,166

 

 

2

Выход

Статическое давление P, МПа

5,014

 

5,012

 

5,019

 

 

3

Стенка

Скорость потока V, м/с

0

 

0

 

0

 

 

 

 

 

 

 

 

 

 

 

 

ГУ на входе, выходе и стенках трехмерных областей 2FSIсистемы при моделировании одной ступени компрессора с учетом газодинамики РУ представлены на рис. 2.8 и в табл. 2.5.

ГУ на контактных поверхностях взаимодействия областей «газ – газ», «газ – конструкция» и «конструкция – конструкция» 2FSI-системы при моделировании одной ступени компрессора с учетом газодинамики РУ представлены на рис. 2.9 и в табл. 2.6.

51

Рис. 2.7. ГУ на контактных поверхностях взаимодействия областей «газ – газ», «газ – конструкция» и «конструкция – конструкция» при учете газодинамики ЛУ

Таблица 2.4

ГУ на контактных поверхностях взаимодействия областей «газ – конструкция» и «конструкция – конструкция» при учете газодинамики ЛУ

Поверх-

Поверх-

Тип взаимо-

Математическая

 

Цветовое

п/п

ность 1

ность 2

действия

формулировка

 

обозначение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P S

i

F1(1)

F S2(1)

 

 

 

 

 

 

Проточная

 

 

i

 

 

 

 

i

 

 

 

(1)

 

1

ЛУ

Fluid – Solid

 

 

 

 

 

 

(1)

 

 

 

 

 

 

 

 

 

 

 

 

частьЛУ

UiS2(1)

 

(1)

UiF1(1)

 

(1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P S

 

F1(2)

F S2(2)

 

 

 

 

Проточная

Покрывной

 

i i

 

 

 

 

 

i

 

(2)

 

2

Fluid – Solid

 

 

 

 

 

 

(2)

 

 

 

 

 

 

 

 

 

 

частьЛУ

диск РК

UiS2(2)

 

(2)

UiF1(2)

 

(2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

Основной

Вал

Solid – Solid

UiS1(3)

 

UiS2(3)

 

 

 

 

 

 

 

 

 

диск РК

 

 

 

 

 

 

 

(3)

 

 

 

(3)

 

 

 

 

 

 

 

 

 

 

 

 

 

52

Окончание табл. 2.4

Поверх-

Поверх-

Тип взаимо-

Математическая

Цветовое

п/п

ность 1

ность 2

действия

формулировка

обозначение

4

Основной

РУ

Solid – Solid

UiS1(4)

 

UiS2(4)

 

 

 

 

диск РК

 

 

 

(4)

 

 

(4)

 

 

 

 

 

 

 

5

РУ

Вал

Solid – Solid

UiS1(5)

(5)

UiS2(5)

 

(5)

 

 

 

 

 

 

 

 

 

Рис. 2.8. ГУ на входе, выходе и стенках при учете газодинамики РУ

Таблица 2.5

ГУ на входе, выходе и стенках при учете газодинамики РУ

 

 

Скоростьвращения ,

 

Цветовое

Тип ГУ

Параметр

 

об/мин

 

 

п/п

 

 

 

 

 

обозначение

 

 

4200

5160

6000

 

 

 

 

 

 

 

 

 

 

 

 

 

1

Вход

Полное давление P*, МПа

5,570

5,843

6,128

 

 

2

Выход

Статическое давление P, МПа

5,066

5,066

5,066

 

 

 

 

 

 

 

 

 

 

3

Стенка

Скорость потока V, м/с

0

0

0

 

 

 

 

 

 

 

 

 

 

53

Рис. 2.9. ГУ на контактных поверхностях взаимодействия областей «газ – газ», «газ – конструкция» и «конструкция – конструкция» при учете газодинамики РУ

Таблица 2.6

ГУ на контактных поверхностях взаимодействия областей «газ – газ», «газ – конструкция» и «конструкция – конструкция» при учете газодинамики РУ

Поверх-

Поверх-

Тип взаимо-

Математическая

 

Цветовое

п/п

ность 1

ность 2

действия

формулировка

 

обозначение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P S

i

F1(1)

F S2(1)

 

 

 

 

 

 

Проточная

 

 

i

 

 

 

 

 

i

 

 

 

(1)

 

1

РУ

Fluid – Solid

 

 

 

 

 

 

 

(1)

 

 

 

 

 

 

 

 

 

 

 

 

 

часть РУ

UiS2(1)

 

(1)

UiF1(1)

 

(1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P S

i

 

F1(2)

F S2(2)

 

 

 

 

Проточная

Основной

 

i

 

 

 

 

 

i

 

(2)

 

2

Fluid – Solid

 

 

 

 

 

 

 

(2)

 

 

 

 

 

 

 

 

 

 

 

часть РУ

диск РК

UiS2(2)

 

UiF1(2)

 

(2)

 

 

 

(2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

54

Окончание табл. 2.6

Поверх-

Поверх-

Тип взаимо-

Математическая

 

Цветовое

п/п

ность 1

ность 2

действия

 

формулировка

 

обозначение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P S

i

F1(3)

 

F S 2(3)

 

 

 

 

Проточная

 

 

i

 

 

 

 

 

i

 

(3)

 

3

Вал

Fluid – Solid

 

 

 

 

 

 

(3)

 

 

 

 

 

 

 

 

 

 

часть РУ

U iS2(3)

 

(3)

U iF1(3)

 

(3)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

Основной

Вал

Solid – Solid

U iS1(4)

 

 

 

 

U iS2(4)

 

 

 

 

 

 

диск РК

 

 

 

 

 

 

(4)

 

 

 

 

(4)

 

 

 

 

 

 

 

 

 

 

 

 

 

5

Основной

РУ

Solid – Solid

UiS1(5)

 

 

 

 

UiS2(5)

 

 

 

 

 

 

 

 

диск РК

 

 

 

 

 

 

 

(5)

 

 

 

(5)

 

 

 

 

 

 

 

 

 

 

 

 

 

6

РУ

Вал

Solid – Solid

U iS1(6)

 

(6)

U iS2(6)

 

(6)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ГУ на входе, выходе и стенках трехмерных областей 2FSI-системы при моделировании одной ступени компрессора с совместным учетом газодинамики РК, ЛУ, РУ представлены на рис. 2.10 и в табл. 2.7.

Рис. 2.10. ГУ на входе, выходе и стенках при учете газодинамики РУ

55

Таблица 2.7

ГУ на входе, выходе и стенках при совместном учете газодинамики РК, ЛУ, РУ

 

 

Скорость вращения

Цветовое

Тип ГУ

Параметр

, об/мин

п/п

 

 

 

 

 

обозначение

 

 

4200

5160

6000

 

 

 

 

 

 

 

1

Вход

Полное давление P*, МПа

5,066

5,066

5,066

 

2

Выход

Массовый расход G, кг/с

200

200

200

 

из РК

 

3

Выход

Статическое давление P,

5,066

5,066

5,066

 

из РУ

МПа

 

4

Стенка

Скорость потока V, м/с

0

0

0

 

 

 

 

 

 

 

 

ГУ на контактных поверхностях взаимодействия областей «газ – газ», «газ – конструкция» и «конструкция – конструкция» 2FSI-системы при моделировании одной ступени компрессора с совместным учетом газодинамики РК, ЛУ, РУ представлены на рис. 2.11 и в табл. 2.8.

Рис. 2.11. ГУ на контактных поверхностях взаимодействия областей «газ – газ», «газ – конструкция» и «конструкция – конструкция» при совместном учете газодинамики РК, ЛУ, РУ

56

Таблица 2.8

ГУ на контактных поверхностях взаимодействия областей «газ – газ», «газ – конструкция» и «конструкция – конструкция» при совместном учете газодинамики РК, ЛУ, РУ

Поверх-

Поверх-

Тип взаимо-

Математическая

Цветовое

п/п

ность 1

ность 2

действия

формулировка

обозначение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

Проточная

Входное

Fluid – Fluid

PF1(1)

 

 

 

PF2(1)

 

 

 

 

 

 

 

 

 

 

 

 

часть РК

устройство РК

i

 

 

 

 

 

(1)

i

(1)

 

 

 

 

 

 

 

 

 

 

 

2

Проточная

Выходное

Fluid – Fluid

PF1(2)

 

 

 

PF2(2)

 

 

 

 

 

 

 

 

 

 

 

 

часть РК

устройство РК

i

 

 

 

 

(2)

i

 

 

 

(2)

 

 

 

 

 

 

 

 

 

 

 

 

 

3

Проточная

Проточная

Fluid – Fluid

P F1(3)

 

 

 

PF2(3)

 

 

 

 

 

 

 

 

 

 

 

 

часть РК

частьЛУ

i

 

 

 

 

(3)

i

 

 

 

 

(3)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

Проточная

Проточная

Fluid – Fluid

P F1(4)

 

 

 

P F2(4)

 

 

 

 

 

 

 

 

 

 

 

 

часть РК

частьЛР

i

 

 

 

 

(4)

i

 

 

 

 

(4)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P S

i

F1(5)

F S2(5)

 

 

 

 

 

Проточная

 

 

i

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

(5)

 

5

Лопатки РК

Fluid – Solid

 

 

 

 

 

 

 

(5)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

часть РК

UiS2(5)

 

(5)

UiF1(5)

 

 

 

 

(5)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P S

i

F1(6)

F S2(6)

 

 

 

 

Проточная

Основной

 

i

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

(6)

 

6

Fluid – Solid

 

 

 

 

 

 

 

(6)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

часть РК

диск РК

UiS2(6)

 

(6)

UiF1(6)

 

 

 

 

(6)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P S

i

F1(7)

F S2(7)

 

 

 

 

Проточная

 

 

i

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

(7)

 

7

Вал

Fluid – Solid

 

 

 

 

 

 

 

(7)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

часть РК

UiS2(7)

 

(7)

UiF1(7)

 

 

 

 

(7)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P S

i

F1(8)

F S2(8)

 

 

 

Проточная

Покрывной

 

i

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

(8)

 

8

Fluid – Solid

 

 

 

 

 

 

 

(8)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

часть РК

диск РК

UiS2(8)

 

(8)

UiF1(8)

 

 

 

 

(8)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P S

i

F1(9)

F S2(9)

 

 

 

Проточная

 

 

i

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

(9)

 

9

ЛУ

Fluid – Solid

 

 

 

 

 

 

 

(9)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

частьЛУ

UiS2(9)

 

(9)

UiF1(9)

 

 

 

 

(9)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P S

 

F1(10)

F S2(10)

 

 

 

Проточная

Покрывной

 

i i

 

 

 

 

 

i

 

 

 

 

 

 

 

 

(10)

 

10

Fluid – Solid

 

 

 

 

 

 

 

(10)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

частьЛУ

диск РК

UiS2(10)

 

(10)

UiF1(10)

 

(10)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

57

Окончание табл. 2.8

Поверх-

Поверх-

Тип взаимо-

Математическая

 

 

 

Цветовое

п/п

ность 1

ность 2

действия

формулировка

 

 

 

обозначение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P S

i

F1(11)

 

F S2(11)

 

 

 

 

 

Проточная

 

 

i

 

 

 

 

 

i

 

 

(11)

 

11

РУ

Fluid – Solid

 

 

 

 

 

 

(11)

 

 

 

 

 

 

 

 

 

 

 

часть РУ

UiS2(11)

 

(11)

UiF1(11)

 

(11)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P S

i

F1(12)

 

F S2(12)

 

 

 

 

 

Проточная

Основной

 

i

 

 

 

 

 

i

 

 

(12)

 

12

Fluid – Solid

 

 

 

 

 

 

(12)

 

 

 

 

 

 

 

 

 

 

 

 

часть РУ

диск РК

UiS2(12)

 

 

(12)

UiF1(12)

 

 

 

 

 

 

 

 

 

 

 

 

 

(12)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P S

i

F1(13)

 

F S2(13)

 

 

 

 

 

 

Проточная

 

 

i

 

 

 

 

 

i

 

 

 

(13)

 

13

Вал

Fluid – Solid

 

 

 

 

 

 

(13)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

часть РУ

UiS2(13)

 

(13)

UiF1(13)

 

 

 

 

 

 

 

 

 

(13)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

14

Основной

Вал

Solid – Solid

UiS1(14)

 

 

 

UiS2(14)

 

 

 

 

 

 

 

 

 

диск РК

 

 

 

 

 

 

 

(14)

 

 

 

 

(14)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

15

Основной

РУ

Solid – Solid

UiS1(15)

 

 

 

UiS2(15)

 

 

 

 

 

 

 

 

 

 

диск РК

 

 

 

 

 

 

 

(15)

 

 

 

 

(15)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

16

РУ

Вал

Solid – Solid

UiS1(16)

 

(16)

UiS2(16)

 

 

(16)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2.3. РАЗРАБОТКА СЕТОЧНЫХ КОНЕЧНО-ЭЛЕМЕНТНЫХ МОДЕЛЕЙ ДЛЯ ГАЗОДИНАМИЧЕСКИХ И ТВЕРДОТЕЛЬНЫХ ОБЛАСТЕЙ РОТОРА ГПА

Для всех 11 элементов конструкции ротора и газодинамических полостей при учете одной ступени компрессора (см. рис. 2.1) построены отдельные сеточные модели. В табл. 2.9 приводятся основные сведения об используемых модулях для разработки моделей с указанием количества конечных элементов (см. рис. 2.1).

Использование специализированных модулей ANSYS и проведение сеточной оптимизации позволило разработать сеточные модели для проведения численного моделирования 2FSI-процессов в компрессоре ГПА с учетом одной ступени при раздельном и совместном рассмотрении газодинамики РК, ЛУ, РУ.

58

Таблица 2.9

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

 

Тип

Модуль

Модуль

Кол-во

Модель

разработки

разработки

конечных

п/п

модели

 

геометрии

сетки

элементов

1

Газодинамическая полость

Fluid

DesignModeler

ICEM CFD

83 520

входного устройства РК

2

Газодинамическая полость

Fluid

DesignModeler

ICEM CFD

222 460

выходного устройства РК

3

Газодинамическая полость

Fluid

BladeGen

TurboGrid

1 573 726

проточной части РК

 

 

 

 

 

4

Газодинамическая

Fluid

DesignModeler

ICEM CFD

694 051

полость РУ

 

 

 

 

 

5

Газодинамическая

Fluid

DesignModeler

ICEM CFD

180 376

полость ЛУ

 

 

 

 

 

6

Основной диск РК

Solid

BladeGen +

ANSYS Mesher

37 701

DesignModeler

7

Покрывной диск РК

Solid

BladeGen +

ANSYS Mesher

38 706

 

 

 

DesignModeler

ANSYS Mesher

 

8

Лопатки РК

Solid

BladeGen

41 582

9

Вал ротора

Solid

DesignModeler

ANSYS Mesher

58 328

10

РУ

Solid

DesignModeler

ANSYS Mesher

132 322

11

ЛУ

Solid

DesignModeler

ANSYS Mesher

22 264

 

Общее количество

конечных элементов

3 085 036

Количество конечных элементов относительно небольшое – ≈3,1 млн при совместном учете, что критически важно для последующих 2FSI-расчетов, так как сама по себе задача совместного динамического моделирования газа и конструкции в 3D-поста- новке крайне ресурсоемка даже с таким количеством ячеек. Предельно ресурсоемкая задача запускалась на 32 восьмиядерных процессорах «Intel Xeon E5-2680» ВВК ПНИПУ [36, 100], при этом использовалось 2 Тб оперативной памяти. Результаты задачи заняли около 2,5 Тб на распределенных жестких дисках.

Меньшее количество элементов в твердых телах (Solid) по сравнению с газодинамическими объемами (Fluid) обусловлено малыми деформациями в конструкции и использованием упру-

59

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

Использование модулей ICEM CFD и TurboGrid для моделирования газодинамических полостей ротора компрессора ГПА позволило разработать гексагональную структурированную сетку с более высоким качеством, где в явном виде прорабатывалась блочная O-grid топология, корректно описывающая особенности течения в пограничном слое при использовании SSTмодели (см. подразд. 2.2).

Использование SST-модели турбулентности и применение специализированных модулей ANSYS, позволивших уменьшить количество конечных элементов без потери качества моделей, позволило более корректно описать поведение 2FSI-системы в сравнении с [74], а также значительно (на 2 порядка) уменьшить потребное время расчета. При этом продолжительность моделируемого отрезка времени увеличилась в 20 раз, что отразилось на разрешающей способности проведения спектрального анализа, который необходим для поиска возможных областей возникновения аэроупругих колебаний ротора компрессора ГПА.

Ниже представлены сеточные модели отдельных элементов конструкции и газодинамических полостей ротора компрессора ГПА при учете одной ступени компрессора.

С целью более качественного и удобного построения сеточной модели геометрия газодинамической полости РК была построена в специализированном программном модуле для турбомашин BladeGen. Построена проточная часть одной лопатки рабочего колеса, которая в дальнейшем размножена и развернута в газодинамическом модуле ANSYS CFX для полноразмерного 3D-моделирования (360°) проточной части РК (рис. 2.12).

Сеточная модель газодинамической полости проточной части РК построена для одного объемного сектора 32,73° в моду-

ле TurboGrid (рис. 2.13) и размножена до 360°.

60

Соседние файлы в папке книги