CC BY
18УДК 62-762.62
Технические науки
Ильичев Владимир Юрьевич, к.т.н., доцент кафедры «Тепловые двигатели и
гидромашины »
Калужский филиал ФГОУ ВО «Московский государственный технический университет имени Н.Э. Баумана (национальный исследовательский
университет), г. Калуга, Россия Качурин Алексей Витальевич, магистрант кафедры «Тепловые двигатели и
гидромашины»
Калужский филиал ФГОУ ВО «Московский государственный технический университет имени Н.Э. Баумана (национальный исследовательский
университет), г. Калуга, Россия
СОЗДАНИЕ И АПРОБАЦИЯ ПРОГРАММЫ ДЛЯ РАСЧЁТА
ХАРАКТЕРИСТИК ВОЗДУШНЫХ ЛАБИРИНТНЫХ УПЛОТНЕНИЙ
Аннотация: В статье приведено описание процесса создания программы для расчёта воздушных лабиринтных уплотнений, применяемых, например, в газотурбинных двигателях (ГТД). Описаны функции языка Python, использованные для написания кода программы (в том числе мощная функция нахождения теплофизических параметров CoolProp), приведены примеры результатов расчёта для одного из уплотнений ГТД при варьировании некоторых исходных параметров. По результатам расчётов сделаны выводы, намечены направления и определены методы дальнейших исследований воздушных лабиринтных уплотнений.
Ключевые слова: лабиринтное уплотнение, газотурбинный двигатель, автоматизация расчётов, язык Python, библиотека CoolProp.
Annotation: The article describes the process of creating a program for calculating air labyrinth seals used, for example, in gas turbine engines (GTE).
Python language functions used to write the program code (including the powerful function of finding the thermophysical parameters CoolProp) are described, examples of calculation results for one of the GTE seals are given when varying some initial parameters. Based on the results of the calculations, conclusions were drawn, directions were outlined and methods for further research of air labyrinth seals were determined.
Keywords: labyrinth seal, gas turbine engine, calculation automation, Python language, library CoolProp.
Введение
Актуальной проблемой при проектировании и эксплуатации газотурбинных двигателей (ГТД) является повышение их экономичности [9]. Важнейшим из способов, позволяющих повысить КПД, является уменьшение утечек воздуха через уплотнения.
В ГТД наибольшее распространение получили воздушные лабиринтные уплотнения, представляющие из себя совокупность последовательно расположенных зазоров, в которых поток воздуха ускоряется, и камер, в которых воздух замедляется с одновременным возникновением завихрений и потерь энергии (а в конечном итоге, давления). Указанные зазоры расположены между роторными и статорными деталями в разных узлах двигателя. Лабиринтные уплотнения стали настолько популярным конструктивным решением благодаря простоте конструкции, высокому ресурсу при расчётном режиме работы двигателя, возможности работы при высоких скоростях и температурах. Существенным недостатком лабиринтных уплотнений является то, что их работа основана на истечении, поэтому в любом случае наблюдаются потери рабочей газообразной среды [1].
Так как во время работы ГТД ротор обычно сильно перемещается в осевом направлении относительно статора, единственной надёжной конструкцией лабиринтных уплотнений в данном случае является её прямоточный вариант (рис. 1).
Корпус
Вал
Рис. 1. Конструктивная схема прямоточного лабиринтного уплотнения: s - шаг гребней, 5 -
радиальный зазор
Эффективность работы таких уплотнений зависит главным образом от конструктивных особенностей и шага их гребней, а также от величины радиального зазора [10]. Минимальная величина радиального зазора при проектировании задаётся для расчётных условий работы двигателя, однако она может существенно увеличиваться, например, при задевании и частичном разрушении гребней уплотнения об ответную часть при возникновении вибраций. Поэтому очень важна разработка методики автоматизированного расчёта уплотнений (компьютерной программы), позволяющей оценить утечки воздуха в зависимости от радиального зазора, а также от других факторов, подбор которых возможен при проектировании, и изменение которых может наблюдаться во время работы двигателя.
В связи с этим, целью данной работы явилась разработка кода программы на универсальном широко распространённом и бесплатном языке Python [12], с применением ряда специальных библиотек функций. Результатами расчётов являются: диаграмма распределения давлений по гребням (использующаяся, например, для расчёта их на прочность), а также зависимость утечки рабочей среды (воздуха) от величины зазора в уплотнении, с учётом срабатывания (уменьшения длины) гребней во время работы и соответствующего увеличения зазора между роторной и статорной деталями. Также в случае необходимости можно построить зависимость утечки от любого другого исходного параметра (например, от шага гребней).
Материал и методы исследования
В качестве основной методики расчёта лабиринтных уплотнений приняты хорошо зарекомендовавшие себя в предыдущих исследованиях [1, 5] руководящие технические материалы [8].
Код созданной программы для расчёта по данной методики включает в себя следующие блоки:
1. Импорт библиотеки математических функций Math, модуля вывода графической информации Matplotlib.pyplot [2] и библиотеки определения теплофизических параметров веществ CoolProp, ранее использованной в работах авторов [3].
2. Определение переменных, обозначающих следующие исходные данные: начальное и конечное давление воздуха, геометрические параметры гребней и др.
3. С помощью математических функций, описанных в руководящих технических материалах, определяется: рекомендуемое число гребней, критическое давление воздуха (от которого зависит выбор расчётного конечного давления), среднегеометрический диаметр зазора уплотнения.
4. Вычисляется давление за каждым гребнем, заполняются массивы номеров гребней и соответствующих давлений и по данным массивам строится диаграмма распределения давления по гребням.
5. При заданным температуре и давлении воздуха перед уплотнением определяется его удельный объём. Здесь очень полезной оказывается функция PropSI [4] модуля CoolProp, позволяющая автоматизировать данный процесс (в отличие от традиционного способа определения свойств газов по таблицам и графикам).
6. Вычисляется ещё ряд параметров, необходимых для нахождения протечек воздуха через уплотнение. Для одного из параметров, характеризующего поправку, учитывающую ширину кромки гребня, в программе Excel построена аппроксимирующая зависимость, представленная на рис. 2.
в
1,4 1,35 1,3 1,25 1,2 1,15 1,1 1,05 1
у = 1,1418х°,(
0,1
1,1 2,1 3,1 4,1 5,1 6,1
Ь/б
Рис. 2. Определение аппроксимирующей зависимости для поправки в к расчетной величине зазора, учитывающей влияние ширины кромки гребня
7. Определяются протечки при разных значениях исходных параметров (в частности, при разной величине зазора и высоте гребня), заполняются массивы значений зазоров и соответствующих протечек и по данным массивам строится диаграмма их взаимозависимости.
Пример расчёта
В качестве примера рассмотрим расчёт уплотнения 5-й опоры универсального газотурбинного модуля МСУ - 800 производства ПАО «Кадви», г. Калуга с применением разработанной программы. Исходные данные к расчёту приведены в табл. 1.
Таблица 1. Исходные данные для расчёта лабиринтного уплотнения.
Параметр Значение
Начальное давление Р0, Па 239282,26
Конечное давление Р, Па 119641,13
Высота гребня Н, м 0,0015
Ширина гребня Ь, м 0,0003
Номинальный зазор б, м 0,00005
Шаг 1, м 0,002
Осевой разбег с, м 0,002
Длина уплотнения Д, м 0,0292
Радиус скругления кромки г, м 0,00001
При номинальном зазоре в уплотнении расчётная протечка составляет 0,0065 кг/с. На рис. 3 представлен построенный в п. 4 рассмотренного кода программы график зависимости давления за каждым гребнем от номера гребня.
Рис. 3. Рассчитанное давление за каждым гребнем лабиринтного уплотнения
Из графика следует, что от гребня к гребню перепад давления увеличивается. Поэтому, например, при дальнейшем расчёте на прочность достаточно рассчитать только последний (15-й) гребень, на который действует максимальный перепад давлений.
Также программа позволяет проследить, каким образом будет изменяться протечка, например, при снашивании гребней уплотнения - то есть, при одновременном увеличении зазора и уменьшении высоты гребней. На рис. 4. показано изменение протечки начиная от номинальной высоты гребней, и заканчивая их полным износом.
■-1-1-1-1-1-1-1-1-1
0.0000 0.0002 0.0004 0.0006 0.0008 0.0010 0.0012 0.0014 0.0016
Зазор в уплотнении,м
Рис. 4. Изменение протечки воздуха через уплотнение при увеличении зазора и уменьшении
высоты гребней (при их износе)
По графику можно сделать вывод, что при увеличении зазора в уплотнении протечки стремительно увеличиваются, а при полном износе гребней (из-за очень большого зазора - около 1,5 мм) уплотнение полностью перестаёт выполнять свою функцию, а протечка по сравнению с расчётным режимом увеличивается в 77 раз. Данная ситуация является недопустимой, -при этом, согласно документации, двигатель нуждается в немедленной остановке, выявлении и устранении причин задевания уплотнений. После этого узлы уплотнений должны быть заменены на новые.
Заключение
Таким образом, цель описываемой научно-исследовательской работы выполнена в полном объёме: изучена методика расчёта лабиринтных уплотнений, согласно которой создана программа, позволяющая определять любые зависимости характеристик уплотнения от исходных параметров как в процессе их проектирования, так и в ходе их эксплуатации [11]. В программе использована выведенная аппроксимационная зависимость, а главное, уникальные средства языка Python, позволяющие полностью автоматизировать
расчётный процесс.
На следующем этапе исследований планируется верификация получаемых с помощью описанной программы данных с помощью проведения широкого ряда газодинамических расчётов в программе Ansys [6]. Разработанную в данной работе методику следует в дальнейшем дополнить расчётом лабиринтных уплотнений не только прямоточного типа, но и разного рода ступенчатых конструкций, с целью сравнения их характеристик.
Разработанную программу можно с успехом применять для обучения как основам расчётов деталей машин, так и современным программным средствам автоматизации проектирования [7].
В результате хотелось бы отметить, что язык программирования Python с огромным количеством библиотек функций, разработанных для него на данный момент, позволяет наглядным и достаточно простым путём решить не только поставленную в данном исследовании задачу, но и множество других проблем в различных областях науки и техники.
Библиографический список:
1. Ильичев В.Ю., Ибрагимов С.В. Влияние угла наклона гребней лабиринтных уплотнений турбин на протечки пара. // Электронный журнал: наука, техника и образование. 2017. № СВ2 (13). С. 15-22.
2. Ильичев В.Ю. Использование рекурсивных функций для создания фрактальной графики средствами языка Python. // Системный администратор. 2021. № 3 (220). С. 92-95.
3. Ильичев В.Ю., Лужецкий А.А. Методика технико-экономического обоснования применения пластинчатых рекуператоров в малоразмерных ГТУ. // Научное обозрение. Технические науки. 2021. № 1. С. 40-45.
4. Ильичев В.Ю. Разработка программных продуктов с использованием модуля Python CoolProp для исследования эффективности утилизации тепла продуктов сгорания газообразных топлив. // Системный администратор. 2020. № 11 (216). С. 80-83.
5. Ильичев В.Ю., Юрик Е.А., Ибрагимов С.В. Разработка методики определения осевых усилий, возникающих в лабиринтных уплотнениях паровых турбин. // СтройМного. 2018. № 1 (10). С. 3.
6. Насонов Д.А. Комбинированный подход к моделированию динамики зубчатых передач. // Известия Тульского государственного университета. Технические науки. 2011. № 5-3. С. 53-58.
7. Насонов Д.А., Леонтьев М.Ю., Раевский В.А., Смоловик А.Е. Моделирование и расчет амплитудно-частотных характеристик корпуса двухступенчатого планетарного редуктора. // Известия Тульского государственного университета. Технические науки. 2018. № 7. С. 218-224.
8. РТМ 108.020.33-86. Уплотнения лабиринтные стационарных паровых и газовых турбин и компрессоров. Проектирование и расчет. - М.: НПО ЦКТИ, 1988.
9. Трошенькин Б.А., Трошенькин В.Б. Совершенствование термодинамических циклов газотурбинных установок. // Современная наука: исследования, идеи, результаты, технологии. 2014. № 2 (15). С. 45-48.
10. Шмаков А.Ф., Модорский В.Я. Исследование зазоров в лабиринтных уплотнениях при действии газодинамических и тепловых нагрузок в FSI постановке. // Научно-технический вестник Поволжья. 2016. № 5. С. 216-218.
11. Greshta V., Tkach D., Sotnikov E., Pavlenko D., Klymov O. Studying and designing improved coatings for labyrinth seals of gasturbine engine turbines. // Eastern-European Journal of Enterprise Technologies. 2018. Т. 4. № 12 (94). С. 5663.
12. Ilyichev V.Yu. Development of a program for Lorentz attractor research and its use. // The Complex Systems. 2021. № 1 (11). С. 58-64.
