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

Работа выполнена на кафедре "Газотурбинные, паротурбинные установки и двигатели" Казанского государственного технического университета им. А.Н. Туполева и в ЗАО "НИИтурбокомпрессор им. В.Б. Шнеппа".

Научный руководитель - доктор технических наук, профессор Евгеньев Станислав Сергеевич

Официальные оппоненты: доктор технических наук, профессор Фафурин А.В., кандидат технических наук Сагадеев Р.Г.

Ведущая организация: ОАО "Авиамотор" (г. Казань)

Защита диссертации состоится " " июня 2006 г. в " " часов на заседании диссертационного совета Д 212.079.02 в Казанском государственном техническом университете им.А.Н. Туполева по адресу: 420111, г. Казань, ул. К. Маркса, 10.

С диссертацией можно ознакомиться в библиотеке Казанского государственного технического университета им. А.Н. Туполева.

Автореферат разослан " " ____________2006 г.

Ученый секретарь диссертационного совета

Общая характеристика работы

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

Опыт эксплуатации ЦК при высоких окружных скоростях вращения и давлениях рабочей среды выявил необходимость исследования нестационарных процессов в связи с усталостными поломками рабочих колес (РК), возникающих из-за наличия значительных переменных аэродинамических нагрузок. Среди известных нестационарных явлений в ЦК наиболее значительное отрицательное влияние на усталостную прочность колеса оказывает неравномерность распределения параметров потока по окружной координате около дисков колеса, связанная со сложной структурой потока на выходе из РК, обратным влиянием выходного устройства в случае безлопаточного диффузора (БЛД) и влиянием лопаточного диффузора (ЛД). При вращении ротора окружная неравномерность параметров потока, стационарная по отношению к корпусу, преобразуется в нестационарное поле давлений и скоростей по отношению к вращающемуся РК. Важно отметить, что этот вид нестационарности существует на всех режимах работы ЦК, допускаемых при эксплуатации. Кроме того, имеется возможность возникновения резонансного режима при совпадении частоты аэродинамической нагрузки с одной из собственных частот колебания какого-либо элемента РК, приводящего к опасным динамическим напряжениям и усталостным разрушениям.

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

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

Задачи исследования. В соответствии с поставленной целью решались следующие задачи:

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

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

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

Научная новизна работы. В процессе исследования получены новые научные результаты, которые выносятся на защиту:

1. В результате экспериментальных исследований трех характерных для ЦК ступеней концевого типа с коэффициентами расхода Ф>0>=0,025-0,07-0,09, содержащих входной аппарат, РК с цилиндрическими и пространственными лопатками, БЛД и ЛД и кольцевую камеру (КК), и анализа экспериментальных данных других авторов получены обобщающие зависимости неравномерности давлений за РК и около дисков РК от реактивности при разных режимах работы ступеней.

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

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

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

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

6. Выявлены резонансные режимы работы трех характерных РК, работающих с БЛД, ЛД и КК на основе построенных частотных диаграмм.

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

Практическая ценность. Разработанный метод расчета позволяет определять напряженно-деформированное состояние (НДС) РК при воздействии на него переменных аэродинамических нагрузок на стадии проектирования новых ЦК и заранее принять необходимые меры для повышения надежности их работы.

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

Реализация работы в промышленности. Методы расчета и программы внедрены в систему автоматического проектирования (САПР) ЦК в ЗАО "НИИтурбокомпрессор им. В.Б. Шнеппа".

Апробация работы. Результаты работы докладывались на XIV-XVI Всероссийских межвузовских научно-технических конференциях "Электромеханические и внутрикамерные процессы в энергетических установках, струйная акустика, диагностика технических систем, приборы и методы контроля природной среды, веществ, материалов и изделий" (г. Казань, Военный артиллеристский университет, 2002 г., 2003 г., 2004 г.), на VI и VII Международных научно-технических конференциях молодых специалистов "Исследование, конструирование и технология изготовления компрессорных машин" (г. Казань, ЗАО "НИИтурбокомпрессор им.В.Б. Шнеппа", 2002 г., 2004 г.), на XII Международной молодежной научной конференции "XII Туполевские чтения" (г. Казань, Казанский государственный технический университет им.А.Н. Туполева, 2004 г.), на XIII Международной научно-технической конференции по компрессоростроению "Компрессорная техника и пневматика в XXI веке", (г. Сумы. СНПО им. Фрунзе, 2004 г.), на Международной научно-технической конференции "Рабочие процессы и технология двигателей" (г. Казань, Казанский государственный технический университет им.А.Н. Туполева, 2005 г.), на II Международной научно-технической конференции "Авиадвигатели XXI века" (г. Москва, ЦИАМ, 2005 г.).

В целом работа докладывалась на кафедре "Газотурбинные, паротурбинные установки и двигатели" (г. Казань, Казанский государственный технический университет им.А.Н. Туполева, 2006 г.) и в ЗАО "НИИтурбокомпрессор им. В.Б. Шнеппа", 2004-2006 г. г.

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

Публикации. По теме диссертации опубликовано 14 печатных работ из них 9 статьи и 5 тезисов докладов.

Объем и структура работы. Диссертация состоит из введения, четырех глав, списка литературы и приложения. Работа содержит страниц машинописного текста, таблиц, рисунков. Список литературы включает 89 наименований. В приложении представлен акт внедрения.

Краткое содержание работы

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

центробежный компрессор рабочее колесо

Первая глава содержит анализ современного состояния проблемы усталостной прочности РК ЦК. Отражены исследования физической картины потока в областях, непосредственно примыкающих к РК, и нестационарности потока, определяемой окружной неравномерностью давления около дисков РК, а также приведены методы определения переменных аэродинамических нагрузок и динамических напряжений, действующих на РК. Отмечен значительный вклад Яновского М.И., Скубачевского Г.С., Жирицкого Г.С., Биргера И.А., Шорра Б.Ф., Раера Г.А., Риса В.Ф., Демьянушко И.В., Образцова И.А. и др. в развитие науки о прочности рабочих колес турбомашин. Отражены результаты проведенных ранее исследований Красильникова В.А., Алемасовой Н.А., Измайлова Р.А., Селезнева К.П., Евгеньева С.С., Ильина А.Л. и др. Все указанные работы способствовали улучшению прочностных и газодинамических характеристик. Однако задача получения метода определения переменных аэродинамических нагрузок и динамических напряжений применительно к центробежным ступеням с РК закрытого и полуоткрытого типа с пространственными лопатками в широком диапазоне коэффициентов расхода (Ф>0>=0,025-0,07-0,09) в них не ставилась. В литературе недостаточно отражены результаты экспериментальных исследований переменных аэродинамических нагрузок современных конструкций ступеней ЦК с РК закрытого и полуоткрытого типа с пространственными лопатками, БЛД, ЛД и КК. Отсутствует инженерный метод определения переменных аэродинамических нагрузок и динамических напряжений РК закрытого и полуоткрытого типа с пространственными лопатками.

В результате анализа сформулированы задачи настоящего исследования.

Во второй главе дано описание проведенного экспериментального исследования неравномерности потока по окружности и радиусу около дисков РК. В качестве объектов исследования приняты три характерные ступени концевого типа, каждая из которых включала: входной аппарат, РК, БЛД или ЛД и КК или внутреннюю улитку (ВУ). Параметры этих ступеней приведены в таблице 1.

Экспериментальные исследования проведены на стенде газодинамики ЗАО "НИИтурбокомпрессор им.В.Б. Шнеппа", аттестованном согласно ГОСТ 24555-81. Стенд оснащен всеми необходимыми системами, позволяющими безопасно эксплуатировать его в диапазоне частот вращения ротора исследуемой ступени n=5000-25000об/мин. Система измерения дает возможность определять величины температуры одновременно в 60 точках отбора, избыточного давления в 100 точках отбора и перепада давления на диафрагме с записью и обработкой их на ПЭВМ в процессе испытания.

В работе использованы следующие методы исследования:

пневмометрический - для измерения поля статических давлений в зазорах между дисками РК и корпусом;

термоанемометрический - для выявления частот и амплитуд пульсаций скорости за РК с частотой опроса 2,5 кГц;

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

Таблица 1

№ ступени

Тип вых. устройства

Тип РК

>л2>

z>2>

z>3>

M>u>

>пол>

(расч)

(расч)

№1

КК

Закрытое с цилиндрическими лопатками

90

0,017

23

0,5-0,7

0,47

0,63

0,27

0,025

№2

КК

Закрытое с пространственными лопатками

45

0,068

13

Вариант с БЛД

-

0,5-0,7

0,67

0,82

0,3

0,07

Вариант с ЛД

10

0,5-0,7

0,67

0,84

0,23

0,062

№3

ВУ

Полу открытое с пространственными лопатками

50

0,059

17

19

1,1

0,65

0,81

0,34

0,09

Результаты экспериментальных исследований показали, что переменные безразмерные давления в зависимости от типа ступени, режимов работы, наличия БЛД или ЛД, покрывного или основного диска могут составлять 10-30% от среднего значения давления по окружности.

На основе собственных экспериментальных данных и данных других авторов построены обобщенные графики в виде зависимости неравномерности давления на (А - амплитуда неравномерности давления, - среднее давление по окружности на данном радиусе) от степени реактивности РК  (рис.1) и коэффициента затухания неравномерности давления , где - на любом радиусе, от радиуса для ступеней с БЛД и ЛД. Полученные зависимости позволяют найти неравномерность давления на любом радиусе для ступеней с подобной геометрией, исследованных нами и другими авторами. Для ступеней с БЛД и пространственными РК полученные данные дополняют приведенные в литературе данные для РК с цилиндрическими лопатками. Для ступеней с ЛД полученные обобщенные зависимости в литературе отсутствовали.

По обобщенным данным для ступеней с БЛД наибольшее влияние на относительную неравномерность давления на выходе из РК при оказывает реактивность ступени . С ее уменьшением, т.е. с ростом уровня скоростей за РК, величина заметно растет. Этот рост более заметен для течения около покрывного диска РК и менее - около основного, что связано с меньшей величиной протечки через уплотнение вала.

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

По полученным обобщенным данным для ступеней с ЛД с увеличением реактивности ступени  величина снижается менее заметно по сравнению с при использовании БЛД (рис.2).

Для основного диска величина К> заметно снижается на участке и не сильно зависит от режима работы ступени.

Проведен гармонический анализ распределения безразмерного статического давления по окружности за РК в области в программной среде Mathcad с использованием команды CFFT - стандартного преобразования Фурье. Выявлены гармоники, имеющие максимальные амплитуды: для ступени №1 с БЛД (, , ) k=1-6, для ступени №2 с БЛД (, , ) k=1-4 (рис.3).

Гармонический анализ распределения статических давлений по окружности за РК в области , полученных другими авторами, показал для ступени , , цилиндрические лопатки, (ЗАО "НИИтурбокомпрессор") k=1-2, для ступени , , цилиндрические лопатки, (АО "Невский завод") k=1-2. Для ступени , , прямые лопатки, (ЗАО "НИИтурбокомпрессор") k=1-5.

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

При исследовании пульсаций скорости термоанемометрическим методом измерительная аппаратура, включающая в себя проволочный датчик, термоанемометр DISA55M, линеаризатор, аналогово-цифровой преобразователь (АЦП) и ПЭВМ с программным пакетом для анализа измеренных сигналов, заимствована на кафедре 22 КГТУ им.А.Н. Туполева. С помощью программы, имеющейся в составе измерительной аппаратуры, исходный сигнал разложен в ряд Фурье и найдена энергия сигналов пульсаций скорости Е>k>. На основе построения графика зависимости Е>k> от кратности гармоник k выявлены k, имеющие максимальные амплитуды (рис.4). Для ступени №2 с БЛД k=1-4. Для ступени №1 с БЛД k=1. В связи с частотой опроса 2500 Гц выявить гармоники с кратностями k=z>2> и 2z>2> и выше при M>u>>2>=0,5 (f>0>=182Гц) не удалось. Однако с уменьшением частоты вращения эти гармоники также проявлялись.

1 - -

2 - -

3 - -

4 - -

5 - -

6 - , -

7 - -

Рис. 2. Зависимость для покрывного диска РК ступеней с БЛД
1, 2, 3, 4 – данные других авторов; 5, 6, 7 – данные автора

Исследование пульсаций давлений на передней стенке за РК проведено с использованием аппаратуры, включающей датчик тензометрического типа, изготовленный в ЗАО "НИИтурбокомпрессор", АЦП и ПЭВМ с пакетом программ, позволяющих проводить гармонический анализ измеренных сигналов. Измерительная система имела частоту опроса 5000 Гц. В процессе измерения исходный сигнал разложен в ряд Фурье и построен график зависимости амплитуды пульсаций давления p>k> от кратности гармоник k (рис.5). На основе этого графика выявлены k, имеющие максимальные амплитуды. Для ступени №2 с БЛД k=1-4, k=z>2> и 2z>2>; для ступени №1 с БЛД k=1-2 и k=z>2>; для ступени №2 с ЛД k=1, k=z>2> и 2z>2>. Измерение пульсаций статических давлений за РК с частотой опроса 5000 Гц подтверждает данные, полученные с помощью термоанемометра, и дополняет их выявленными кратностями гармоник k=z>2> и 2z>2> имеющих максимальные амплитуды.

Рис. 3. Гармонический анализ зависимости за РК на

а) ступень №1 с БЛД (,)

б) ступень №2 с БЛД (, )

а) при , б) при , в) при

1 - датчик в центре (b>3>/2), 2 - датчик около покрывного диска.

Рис. 4. Зависимость Е>k>=f(k) на для ступени №2 с БЛД (, )

Рис. 5. Зависимость p>k>=f(k) на при

а) Ступень №1 с БЛД (, ), f>0>=182,4 Гц; 1 - ; 2 - ; 3 -

б) Ступень №2 с БЛД (, ), f>0>=182,4 Гц; 1 - , 2 - , 3 -

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

f>k>=kn,

где n - частота вращения РК (об/с), и использованы в дальнейшем для определения резонансов РК.

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

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

. (1)

и величина аэродинамической нагрузки с учетом амплитуды неравномерности по формуле

. (2)

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

Определение мгновенных изменений давлений на диски и лопатки в межлопаточных каналах РК, возникающих за счет разного противодавления из-за неравномерности давления на выходе РК за 1 оборот проводилось с помощью программного комплекса, в котором решаются задачи: расчета осесимметричного потока в меридиональном канале РК; обтекания решеток профилей на осесимметричных поверхностях тока. Среднее давление на выходе из РК при определялось на основе расчета потока на осесимметричных поверхностях тока при известных начальных параметрах (p>, Т>, М>u>, R, k, Ф>0>) и заданной геометрии ступени; оно сравнивалось со средним давлением p>2>, известным из опытной газодинамической характеристики ступени. Это отличие определялось коэффициентом коррекции К>, который меньше 1, т.к. при расчете не учитываются потери. К> рассчитывался по формуле

, (3)

С использованием величины А при , определяемой по методу, представленному выше, вычислялось максимальное и минимальное давления p>2> с учетом нестационарности по формуле

,. (4)

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

,. (5)

Заданные значения и достигали за счет подбора расчетом потока на осесимметричных поверхностях тока двух мгновенных за 1 оборот значений Ф>0> и соответствующих мгновенных эпюр давлений и скоростей, в канале РК. Последние сравнивали со средними значениями, полученными для среднего за 1 оборот давления по окружности за РК, соответствующего исходному Ф>0>. В результате находили предполагаемое мгновенное значение давления p> и относительной скорости W> по ширине канала на данном , которые, в последующем, использовались для расчета динамических напряжений в межлопаточных отсеках РК. Рассчитанные распределения относительных скоростей в каналах РК сравнивались с экспериментальными данными Красильникова В.А., Локшина И.Л., Тарасова А.Д. Получено удовлетворительное качественное совпадение, что позволяет использовать эту математическую модель.

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

В четвертой главе приведено построение КЭ модели РК. Качество ее построения проверено согласованием расчетов по МКЭ (рис.6) с опытными данными по собственным частотам и формам колебаний исследованных РК, полученными методом голографической интерферометрии (рис.7).

Рассчитанные без учета вращения с помощью МКЭ собственные частоты f> и формы колебаний РК ступени №2, представленные в таблице 2, хорошо согласуются с экспериментальными частотами f> и формами колебаний. Сравнение экспериментальных и расчетных частот и форм колебаний показывает хорошее совпадение, что позволяет говорить о применимости, в дальнейшем, расчетного метода для определения собственных частот и форм колебаний.

Таблица 2

Число узловых диаметров

РК ступени №2 ()

f>, Гц

f>, Гц

2

1856

1842

0,7

3

3671

3643

0,7

4

5080

4815

5,2

5

6010

5615

6,5

МКЭ, реализованный в ANSYS, позволил рассчитать также собственные частоты f>' и формы колебания РК с учетом его вращения. При вращении жесткость РК увеличивается, что приводит к росту собственных частот.

На основе определенных в данной работе f>k> f> f>'=f (n) построены частотные диаграммы (диаграммы Кэмпбелла) для РК трех исследованных ступеней. В точке пересечения этих кривых определены частоты вращения РК, при которых возникнет резонанс.

Проведен анализ динамической и статической прочности РК трех исследованных ступеней, работающих в реальных компрессорах. В результате этого анализа обнаружено, что максимальная интенсивность динамических напряжений > наблюдается на периферии РК в месте стыка лопатки с покрывным диском (рис.8 РК ступени №2). Также определены запасы усталостной прочности для покрывного диска исследованных РК по формуле

, (6)

где - предельное амплитудное напряжение при асимметричных циклах нагружения для покрывного диска РК с лопаткой из стали 07Х16Н6.

определяется из диаграммы выносливости по рассчитанной величине статического напряжения >m>. Запасы прочности по динамическим и статическим напряжениям для трех исследованных РК приведены в таблице 3.

Таблица 3

РК ступени №1

РК ступени №2

РК ступени №3

Запасы по динамическим напряжениям n>a>

13,8

1,9

12,4

Запасы по статическим напряжениям n>m>

1,54

1,67

1,36

Из таблицы 3 видно, что для РК ступени №2 запас усталостной прочности покрывного диска , что меньше нормированного значения [n>] =3 и прочность покрывного диска не выполняется. В связи с этим возможно разрушение покрывного диска в периферийной зоне РК. Для увеличения усталостной прочности РК необходимо изменить его конструкцию.

Выводы

1. Экспериментально исследована неравномерность поля давлений и скоростей по окружности около дисков РК и на его периферии с помощью пневмометрических и безынерционных измерительных приборов для трех, широко применяемых компактных центробежных ступеней концевого типа с БЛД и ЛД, ВУ и КК, охватывающих три характерных для ЦК значения коэффициента расхода Ф>0> (0,025-0,07-0,09), в диапазоне M>u>=0,5-0,7 и M>u>=1,1 для ступени №3 и имеющих закрытые и полуоткрытые РК с цилиндрическими и пространственными лопатками. В результате полученных экспериментальных данных показано, что переменная составляющая давления для исследованных ступеней может достигать 10-30% от средних значений по окружности за РК.

2. В результате анализа собственных экспериментальных данных и данных других авторов для закрытых и полуоткрытых РК с углами >л2>=90, 60, 50, 45, 32 получены значения амплитуд неравномерности давления по окружности (от влияния ВУ и ЛД) на участках от наружного радиуса РК до радиуса уплотнения. Выявлено, что неравномерность давления является максимальной около покрывного диска РК при и использовании БЛД.

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

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

5. Выполнен расчет по МКЭ статических и динамических напряжений и соответствующих запасов прочности на основе действия центробежных сил и поля переменных давлений, действующих с внешней и внутренней стороны дисков трех исследованных типов РК, работающих в составе реальных ЦК, и даны рекомендации по совершенствованию конструкций. В частности, для РК ступени №2 с БЛД запас по динамическим напряжениям составляет 1,9, что требует изменения конструкции.

6. С помощью пневмометрических и безынерционных систем измерений давлений и скоростей в области выхода из РК и анализа спектра частот их колебаний по ширине канала определены номера гармоник k частот воздействия f=kn переменных давлений, соответствующих высоким амплитудам, для исследованных типов РК на разных режимах работы. Для РК, работающих с БЛД и ВУ (или КК) они находятся в диапазоне от k=1 до k=6 в зависимости от степени реактивности  и k=z>2> (z>2> - число лопаток РК), а работающих с ЛД составляют k=1, k=z>2> и k=z>3> (z>3> - число лопаток ЛД).

7. Методом голографической интерферометрии определены собственные частоты и формы колебаний исследованных трех типов закрытых (>л2>=90 и 45) и полуоткрытого ( (>л2>=50) РК и проведено их сравнение с результатами расчета собственных частот и форм колебаний по МКЭ с использованием пакета ANSYS (лицензия № 24281/101217 ЗАО "НИИтурбокомпрессор им.В.Б. Шнеппа"). Совпадение опытных и расчетных данных в пределах 3-5% показало приемлемость построенной КЭ-модели и она использована для расчета собственных частот и форм колебаний с учетом вращения РК.

8. В результате анализа построенных частотных диаграмм выявлены резонансные режимы работы трех исследованных типов РК с >л2>=90, 45 и 50. Например, для ступени с закрытым РК (>л2>=90, , цилиндрические лопатки, БЛД и КК) в диапазоне n=11000-15000об/мин определены три резонансные частоты вращения n=13811об/мин (k=5, два узловых диаметра), n=11854об/мин (k=z>2>=23, 5 узловых диаметров) и n=14251об/мин (k= z>2>=23, 6 узловых диаметров). Широкое закрытое РК (>л2>=45, , пространственные лопатки, ЛД и КК) имело только одну резонансную частоту вращения n=11240об/мин (k=z>3>=10, 2 узловых диаметра). В полуоткрытом РК (>л2>=50, , ЛД, пространственные лопатки) в интервале n=15000-25000об/мин резонансные частоты вращения проявлялись при n=21000об/мин (k=z>2>=17, зонтичные колебания) и n=15000об/мин (k=z>2>=17, 6 узловых диаметров; k=z>3>=19, 8 узловых диаметров). Для этого РК выявлены также характерные формы колебаний лопаток входного вращающегося аппарата.

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

По теме диссертации опубликованы следующие работы

  1. Футин В.А. Анализ погрешностей измерений при газодинамических испытаниях сменных проточных частей для ГПА. - Материалы докладов ХIV Всероссийской межвузовской научно-технической конференции "Внутрикамерные процессы в энергетических установках, акустика, диагностика, экология".Ч.I. Казань, 2002, с.89-92.

  2. Футин В.А. Изучение нестационарных процессов за рабочим колесом центробежного компрессора. - Тезисы докладов VI научно-технической конференции молодых специалистов "Исследование, конструирование и технология изготовления компрессорных машин". Казань, 2002, с.18-19.

  3. Футин В.А. Экспериментальное исследование распределения давления за рабочим колесом центробежного компрессора. - Сборник материалов XV Всероссийской межвузовской научно-технической конференции "Электромеханические и внутрикамерные процессы в энергетических установках, струйная акустика и диагностика, приборы и методы контроля природной среды, веществ, материалов и изделий". Ч. II. Казань, 2003, с.36-37.

  4. Футин В.А. Экспериментальное определение неравномерности распределение давления по окружности около дисков закрытых рабочих колес центробежных ступеней концевого типа. - Сборник материалов XVI всероссийской межвузовской научно-технической конференции "Электромеханические и внутрикамерные процессы в энергетических установках, струйная акустика, диагностика технических систем, приборы и методы контроля природной среды, веществ, материалов и изделий".Ч.I. Казань, 2004, с.231-232.

  5. Евгеньев С.С., Футин В.А. Экспериментальное определение пульсаций давлений и скоростей за рабочим колесом центробежной ступени концевого типа. - Сборник материалов XVI всероссийской межвузовской научно-технической конференции "Электромеханические и внутрикамерные процессы в энергетических установках, струйная акустика, диагностика технических систем, приборы и методы контроля природной среды, веществ, материалов и изделий".Ч.I. Казань, 2004, с.229-231.

  6. Евгеньев С.С., Футин В.А. Экспериментальное исследование пульсаций скоростей и статических давлений за рабочим колесом центробежного компрессора. - Тезисы докладов VII международной научно-технической конференции молодых специалистов "Исследование, конструирование и технология изготовления компрессорных машин" посвященная 80-летию со дня рождения В.Б. Шнеппа, Казань, 2004, с.

  7. . Евгеньев С.С., Футин В.А. К определению неравномерности распределения давления по окружности около дисков закрытых рабочих колес центробежных компрессоров. - Тезисы докладов VII международной научно-технической конференции молодых специалистов "Исследование, конструирование и технология изготовления компрессорных машин" посвященная 80-летию со дня рождения В.Б. Шнеппа, Казань, 2004, с.

  8. Евгеньев С.С., Футин В.А. Исследование неравномерности потока за рабочим колесом центробежного компрессора. - Сборник научных трудов под редакцией доктора техн. наук И.Г. Хисамеева "Проектирование и исследование компрессорных машин". Вып.5, Казань, 2004, с.124-139.

  9. Евгеньев С.С., Футин В.А. Исследование пульсаций давлений и скоростей за рабочим колесом ступени концевого типа центробежного компрессора. - Труды XIII Международной научно-технической конференции по компрессоростроению "Компрессорная техника и пневматика в XXI веке".Ч.I., Сумы, 2004, с.83-90.

  10. Евгеньев С.С., Футин В.А. Распределение давления по окружности около дисков закрытых рабочих колес ступеней концевого типа центробежных компрессоров. - Труды XIII Международной научно-технической конференции по компрессоростроению "Компрессорная техника и пневматика в XXI веке".Ч.I., Сумы, 2004, с.125-131.

  11. Евгеньев С.С., Футин В.А. Определение резонансных режимов рабочих колес центробежных компрессоров, испытывающих воздействие переменных давлений. - Материалы Международной молодежной научной конференции "XII Туполевские чтения". Том I., Казань, 2004, с.221-223.

  12. Евгеньев С.С., Футин В.А. Распределение давления по окружности около дисков закрытых рабочих колес ступеней концевого типа центробежных компрессоров. // - Компрессорная техника и пневматика., № 8, 2004, с.28-30.

  13. Евгеньев С.С., Футин В.А., Каримов А.Х., Макаева Р.Х., Царева А.М. Определение резонансных частот вращения закрытых рабочих колес центробежных компрессоров. - Тезисы докладов Международной научно-технической конференции "Рабочие процессы и технология двигателей", 23-27 мая 2005 г. КГТУ им.А.Н. Туполева, Казань. Казань: Изд-во Казан. гос. техн. ун-та, 2005, с. 198-200.

  14. Евгеньев С.С., Футин В.А. Определение переменных аэродинамических нагрузок и динамических напряжений, действующих на рабочие колеса центробежных компрессоров. - II Международная научно-техническая конференция "Авиадвигатели XXI века". Сб. тезисов, том I. - М.: ЦИАМ, 2005, с.323-325.