книги / Численное моделирование колебательных 2FSI-процессов в компрессорах газоперекачивающих агрегатов
..pdfТаблица 2.17
Значения максимальных амплитуд перемещений конструкции при учете газодинамики РУ (2FSI)
|
|
|
ω=4200об/мин |
|
|
||
Точка |
Переме- |
Частота, |
Переме- |
Частота, |
Переме- |
Частота, |
|
щенияUx, |
щенияUy, |
щенияUz, |
|||||
|
|||||||
|
мкм |
Гц |
мкм |
Гц |
мкм |
Гц |
|
|
|
|
|
||||
На границе |
|
|
|
|
|
|
|
контакта РУ |
|
|
|
|
|
|
|
с зазором лаби- |
4,938 |
100 |
0,116 |
500 |
5,554 |
400 |
|
ринтного уплот- |
|
|
|
|
|
|
|
нения РУ |
|
|
|
|
|
|
|
В центре |
|
|
|
|
|
|
|
торцевого |
3,890 |
100 |
0,127 |
500 |
4,417 |
500 |
|
сечения РУ |
|
|
|
|
|
|
|
На оси |
0,078 |
100 |
0,043 |
100 |
3,827 |
500 |
|
вращения вала |
|||||||
|
|
|
|
|
|
||
|
|
|
ω=5160об |
/мин |
|
|
|
На границе |
|
|
|
|
|
|
|
контакта РУ |
|
|
|
|
|
|
|
с зазором лаби- |
6,567 |
600 |
0,21 |
500 |
8,454 |
500 |
|
ринтного уплот- |
|
|
|
|
|
|
|
нения РУ |
|
|
|
|
|
|
|
В центре |
|
|
|
|
|
|
|
торцевого |
5,172 |
600 |
0,231 |
500 |
7,704 |
500 |
|
сечения РУ |
|
|
|
|
|
|
|
На оси |
0,104 |
100 |
0,069 |
100 |
6,821 |
500 |
|
вращения вала |
|||||||
|
|
|
|
|
|
||
|
|
|
ω=6000об |
/мин |
|
|
|
На границе |
|
|
|
|
|
|
|
контакта РУ |
8,747 |
600 |
0,248 |
500 |
9,974 |
500 |
|
с зазором лаби- |
|||||||
ринтного уплот- |
|
|
|
|
|
|
|
нения РУ |
|
|
|
|
|
|
|
В центре |
|
|
|
|
|
|
|
торцевого |
6,888 |
600 |
0,277 |
500 |
9,119 |
500 |
|
сечения РУ |
|
|
|
|
|
|
|
На оси |
0,132 |
100 |
0,093 |
100 |
8,095 |
500 |
|
вращения вала |
|||||||
|
|
|
|
|
|
81
Таблица 2.18
Сравнение работы газодинамических сил от времени в 1FSI- и 2FSI-постановках при учете газодинамики РУ
1FSI-постановка при учете |
2FSI-постановка при учете |
||
газодинамики РУ |
|
газодинамики РУ |
|
График зависимости работы |
Возрастание (+) затухание (–) колебаний |
График зависимости рабо- |
Возрастание (+) затухание (–) колебаний |
газодинамических |
ты газодинамических сил |
||
сил от времени |
от времени |
||
|
|
||
|
|
|
|
|
+ |
|
+ |
ω = 4200 об/мин |
|
ω = 4200 об/мин |
|
|
+ |
|
+ |
ω = 5160 об/мин |
|
ω = 5160 об/мин |
|
|
+ |
|
+ |
ω = 6000 об/мин |
|
ω = 6000 об/мин |
|
Анализ аэроупругих колебаний ротора компрессора ГПА (учет всей конструкции и газодинамики РУ) показал, что по результатам проведения 2FSI-расчетов на всех режимах работа газодинамических сил не изменила знак. Вместе с тем, в 2FSI-подходе исчезла высокочастотная составляющая характера изменения работы от времени,аамплитудаколебанийуменьшилисьв1,5...2раза.
82
По максимальным значениям амплитуд перемещений ротора на различных режимах вращения (см. табл. 2.17) построены критические области работы системы с учетом газодинамики РУ.
Критическая область работы системы с учетом газодинамики РУ ограничена допустимыми с точки зрения условий эксплуатации (до 40 мкм) [40–46, 75] значениями перемещений, и представлена на рис. 2.29.
а |
б |
в |
Рис. 2.29. Гистограммы максимальных перемещений при различных скоростях вращения ротора ГПА при учете газодинамики РУ
в 2FSI-постановке (Ux, Uy – радиальные направления, Uz – осевое направление): а – на границе контакте РУ с зазором лабиринтного уплотнения РУ, б – в центре торцевого сечения РУ, в – на оси вращения вала
На оси вращения вала осевые перемещения значительно больше радиальных и составляют до 8 мкм. При учете газодинамики РК без газодинамики РУ осевые перемещения составляли до 19 мкм, т.е. даже отдельный учет РУ при относительно небольших массово-габаритных размерах вносит существенный вклад в колебания вала ротора компрессора ГПА. Характер перемещений вблизи лабиринтного уплотнения РУ и в центре торцевого сечения РУ близок, а максимальное значение не превышает 10 мкм, что меньше предельного 40 мкм, обусловленного эксплуатационными особенностями ротора компрессора ГПА [40–46, 75].
83
2.9. ИССЛЕДОВАНИЕ АЭРОУПРУГИХ КОЛЕБАТЕЛЬНЫХ ПРОЦЕССОВ В ОДНОСТУПЕНЧАТОМ ВАРИАНТЕ КОНСТРУКЦИИ КОМПРЕССОРА С СОВМЕСТНЫМ УЧЕТОМ ГАЗОДИНАМИКИ РАБОЧЕГО КОЛЕСА,ЛАБИРИНТНОГО УПЛОТНЕНИЯ
И РАЗГРУЗОЧНОГО УСТРОЙСТВА.ВЗАИМОВЛИЯНИЕ 2FSI-ПОДСИСТЕМ ПРИ РАЗДЕЛЬНОМ И СОВМЕСТНОМ УЧЕТЕ ГАЗОДИНАМИКИ РК,ЛУ,РУ
При исследовании колебательных процессов ротора с учетом газодинамики рабочего колеса (РК), лабиринтного уплотнения (ЛУ) и разгрузочного устройства (РУ) учитывались следующие элементы:
–вал;
–опоры;
–рабочее колесо;
–лабиринтное уплотнение;
–разгрузочное устройство;
–газодинамика рабочего колеса;
–газодинамика лабиринтного уплотнения;
–газодинамика разгрузочного устройства.
Расчеты проведены в 2FSI-постановке для режимов 4200, 5160, 6000 об/мин. Моделируемый отрезок времени составил 0,01 с с шагом 5·10–5 с.
На рис. 2.30, 2.31 представлены поля скорости, давления и линии тока в газодинамической области ротора на скорости вращения компрессора 5160 об/мин.
На следующем этапе был проведен спектральный анализ следующих динамически изменяющихся внутренних переменных, отражающих поведение 2FSI-системы:
– реакции опор в подшипниках (АЧХ представлены
втабл. 2.19, 2.20);
–перемещения конструкции ротора в точках (рис. 2.32): на лопатке РК – точка 1, на границе контакта ЛУ с зазором ЛУ – точка 2, на границе контакта покрывного диска РК с зазором ЛУ –
84
точка 3, на границе контакте РУ с зазором лабиринтного уплотнения РУ – точка 4, в центре торцевого сечения РУ – точка 5, и на оси вращения вала – точка 6 (максимальные значения амплитуд перемещений при раздельном и совместном учете РК, ЛУ, РУ представлены в табл. 2.21);
–работа газодинамических сил, определяемая суперпозицией по направлениям произведений газодинамических сил, действующих на РК, ЛУ, РУ и перемещений ротора в точке на оси вращения вала. Приращение работы за период показывает возрастание или затухание колебательных процессов, что важно для оценки виброустойчивости роторов компрессора ГПА (графики при раздельном и совместномучетеРК,ЛУ, РУпредставленывтабл.2.22).
Граничные условия представлены в табл. 2.7–2.8.
а |
б |
|
Рис. 2.30. Распределение газодинамических параметров
вгазодинамической области ротора: а – поле скорости,
б– поле давления (ANSYS CFX-Post 18.0)
Рис. 2.31. Распределение 3D-линий тока в газодинамической области ротора (ANSYS CFX-Post 18.0)
85
Таблица 2.19
АЧХ реакций опоры подшипника со стороны РУ при учете газодинамики РК, ЛУ, РУ (2FSI)
ω = 4200 об/мин |
ω = 5160 об/мин |
ω = 6000 об/мин |
Fx
Fy
Fz
Таблица 2.20
АЧХ реакций опоры подшипника со стороны трансмиссии при учете газодинамики РК, ЛУ, РУ (2FSI)
ω = 4200 об/мин |
ω = 5160 об/мин |
ω = 6000 об/мин |
Fx
Fy
Fz
86
Точка 2 |
Точка 1 |
|
|
|
Точка 4 |
Точка 3 |
Точка 5 |
Точка 6
Рис. 2.32. Расположение точек в твердотельной модели ротора при учете газодинамики РК, ЛУ и РУ (ANSYS Mechanical 18.0)
Таблица 2.21
Значения максимальных амплитуд перемещений конструкции при раздельном (подразд. 2.6–2.8) и совместном (настоящий раздел) учете газодинамики РК, ЛУ, РУ (2FSI)
Точка |
Тип |
Переме- |
Частота, |
Переме- |
Частота, |
Переме- |
Частота, |
|
учета |
щения |
Гц |
щения |
Гц |
щения |
Гц |
||
|
Ux, мкм |
Uy, мкм |
Uz, мкм |
|||||
|
|
|
|
|
||||
|
|
ω = |
4200 об/ |
мин |
|
|
|
|
|
Раздельно |
Режим |
50 |
16,944 |
50 |
11,342 |
50 |
|
На лопатке РК |
резонанса |
|||||||
Совмест- |
Режим |
100 |
17,150 |
100 |
17,160 |
100 |
||
|
||||||||
|
но |
резонанса |
||||||
На границе |
Раздельно |
0,00067 |
3200 |
0,00197 |
1800 |
0 |
– |
|
контакта ЛУ |
Совмест- |
0 |
– |
0 |
– |
0 |
– |
|
с зазором ЛУ |
но |
|||||||
|
|
|
|
|
|
|||
На границе |
Раздельно |
5,002 |
1700 |
1,937 |
3200 |
13,476 |
100 |
|
контакта по- |
Совмест- |
Режим |
|
|
|
|
|
|
крывного диска |
100 |
10,89 |
100 |
12,159 |
100 |
|||
РК с зазором ЛУ |
но |
резонанса |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
На границе |
Раздельно |
4,938 |
100 |
0,116 |
500 |
5,554 |
400 |
|
контакта РУ |
|
|
|
|
|
|
|
|
с зазором лаби- |
Совмест- |
Режим |
100 |
12,945 |
100 |
8,861 |
100 |
|
ринтного уплот- |
но |
резонанса |
||||||
нения РУ |
|
|
|
|
|
|
|
87
Продолжение табл. 2.21
Точка |
Тип |
Переме- |
Частота, |
Переме- |
Частота, |
Переме- |
Частота, |
|
учета |
щения |
Гц |
щения |
Гц |
щения |
Гц |
||
|
Ux, мкм |
Uy, мкм |
Uz, мкм |
|||||
|
|
|
|
|
||||
В центре |
Раздельно |
3,890 |
100 |
0,127 |
500 |
4,417 |
500 |
|
торцевого |
Совмест- |
Режим |
100 |
10,191 |
100 |
8,824 |
100 |
|
сечения РУ |
но |
резонанса |
||||||
На оси |
Раздельно |
2,882 |
100 |
1,768 |
150 |
7,872 |
50 |
|
вращения |
Совмест- |
0,580 |
100 |
0,316 |
100 |
8,663 |
100 |
|
вала |
но |
|
|
|
|
|
|
|
|
|
|
ω = |
5160 об/мин |
|
|
||
|
Раздельно |
Режим |
50 |
19,617 |
50 |
11,342 |
600 |
|
На лопатке РК |
резонанса |
|||||||
Совмест- |
Режим |
100 |
22,437 |
100 |
26,552 |
100 |
||
|
||||||||
|
но |
резонанса |
||||||
На границе |
Раздельно |
0,00111 |
1800 |
0,00325 |
1800 |
0 |
– |
|
контакта ЛУ |
Совмест- |
0 |
– |
0 |
– |
0 |
– |
|
с зазором ЛУ |
но |
|||||||
|
|
|
|
|
|
|||
На границе |
Раздельно |
7,708 |
1700 |
2,799 |
3200 |
21,493 |
100 |
|
контакта по- |
Совмест- |
Режим |
|
|
|
|
|
|
крывного диска |
100 |
14,398 |
100 |
17,382 |
100 |
|||
РК с зазором ЛУ |
но |
резонанса |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
На границе |
Раздельно |
6,567 |
600 |
0,210 |
500 |
8,454 |
500 |
|
контакта РУ |
|
|
|
|
|
|
|
|
с зазором лаби- |
Совмест- |
Режим |
100 |
16,957 |
100 |
11,189 |
100 |
|
ринтного уплот- |
но |
резонанса |
||||||
нения РУ |
|
|
|
|
|
|
|
|
В центре |
Раздельно |
5,172 |
600 |
0,231 |
500 |
7,704 |
500 |
|
торцевого |
Совмест- |
Режим |
100 |
13,366 |
100 |
11,345 |
100 |
|
сечения РУ |
но |
резонанса |
||||||
На оси |
Раздельно |
0,458 |
50 |
0,444 |
50 |
9,270 |
50 |
|
вращения |
Совмест- |
0,104 |
100 |
0,33 |
100 |
11,195 |
100 |
|
вала |
но |
|
|
|
|
|
|
|
|
|
ω = |
6000 об/ |
мин |
|
|
|
|
|
Раздельно |
Режим |
75 |
22,3 |
75 |
18,63 |
600 |
|
На лопатке РК |
резонанса |
1 |
||||||
Совмест- |
Режим |
400 |
9,072 |
400 |
46,27 |
400 |
||
|
||||||||
|
но |
резонанса |
5 |
|||||
На границе |
Раздельно |
0,00362 |
3300 |
0,00988 |
3300 |
0,00001 |
1800 |
|
контакта ЛУ |
Совмест- |
|
|
|
|
|
|
|
с зазором |
0 |
– |
0 |
– |
0 |
– |
||
ЛУ |
но |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
88
Окончание табл. 2.21
Точка |
Тип |
Переме- |
Частота, |
Переме- |
Частота, |
Переме- |
Частота, |
учета |
щения |
Гц |
щения |
Гц |
щения |
Гц |
|
|
Ux, мкм |
Uy, мкм |
Uz, мкм |
||||
|
|
|
|
|
|||
На границе |
Раздельно |
25,429 |
1800 |
8,704 |
3300 |
14,349 |
1800 |
контакта по- |
Совмест- |
Режим |
|
|
|
|
|
крывного диска |
400 |
3,722 |
400 |
32,003 |
400 |
||
РК с зазором ЛУ |
но |
резонанса |
|
|
|
|
|
|
|
|
|
|
|
|
|
На границе |
Раздельно |
4,938 |
100 |
0,116 |
500 |
5,554 |
400 |
контакта РУ |
|
|
|
|
|
|
|
с зазором лаби- |
Совмест- |
Режим |
400 |
0,761 |
400 |
14,206 |
1500 |
ринтного уплот- |
но |
резонанса |
|||||
нения РУ |
|
|
|
|
|
|
|
В центре |
Раздельно |
6,888 |
600 |
0,277 |
500 |
9,119 |
500 |
торцевого |
Совмест- |
Режим |
400 |
0,709 |
400 |
12,900 |
400 |
сечения РУ |
но |
резонанса |
|||||
На оси |
Раздельно |
Режим |
75 |
22,300 |
75 |
18,631 |
600 |
резонанса |
|||||||
вращения |
|
|
|
|
|
|
|
Совмест- |
|
|
|
|
|
|
|
вала |
0,110 |
400 |
0,009 |
400 |
13,466 |
400 |
|
|
но |
|
|
|
|
|
|
Таблица 2.22
Сравнение работы газодинамических сил от времени при раздельном (подразд. 2.6–2.8) и совместном (настоящий раздел) учете газодинамики РК, ЛУ, РУ (2FSI)
2FSI-постановка при раздельном |
2FSI-постановка при совместном |
|||
учете газодинамики РК / ЛУ / РУ |
учете газодинамики РК, ЛУ, РУ |
|||
График зависимости работы |
Возрастание (+) затухание (–) колебаний |
График зависимости рабо- |
Возрастание (+) затухание (–) колебаний |
|
газодинамических |
ты газодинамических сил |
|||
сил от времени |
от времени |
|||
|
|
|||
|
ω = 4200 |
об/мин |
|
|
Работа газодинамических сил на РК |
Работа газодинамических сил на РК, |
|||
выделенная из совместной модели |
||||
|
|
|||
|
– |
|
– |
|
|
|
|
|
89
Продолжение табл. 2.22
2FSI-постановка при раздельном |
2FSI-постановка при совместном |
|||
учете газодинамики РК / ЛУ / РУ |
учете газодинамики РК, ЛУ, РУ |
|||
График зависимости работы |
Возрастание(+) затухание(–) колебаний |
График зависимости рабо- |
Возрастание(+) затухание(–) колебаний |
|
газодинамических |
|
ты газодинамических сил |
|
|
сил от времени |
|
от времени |
|
|
|
|
|
|
|
Работа газодинамических сил на ЛУ |
Работа газодинамических сил на ЛУ, |
|||
выделенная из совместной модели |
||||
|
|
|||
|
+ |
|
+ |
|
|
|
|
|
|
Работа газодинамических сил на РУ |
Работа газодинамических сил |
на РУ, |
||
выделенная из совместной модели |
||||
|
|
|||
|
+ |
|
– |
|
|
|
|
|
|
|
ω = 5160 |
об/мин |
|
|
Работа газодинамических сил на РК |
Работа газодинамических сил на РК, |
|||
выделенная из совместной модели |
||||
|
|
|||
|
+ |
|
– |
|
|
|
|
|
|
Работа газодинамических сил на ЛУ |
Работа газодинамических сил на ЛУ, |
|||
выделенная из совместной модели |
||||
|
|
|||
|
+ |
|
+ |
|
|
|
|
|
90