CC BY
67
2008 НАУЧНЫЙ ВЕСТНИК МГТУ ГА № 135
серия Эксплуатация воздушного транспорта и ремонт авиационной техники. Безопасность полетов
УДК 629.735
ОПТИМИЗАЦИЯ КОЛИЧЕСТВА ЛОПАТОК КЛАССИЧЕСКОЙ СТУПЕНИ ВЕНТИЛЯТОРА, С ТОЧКИ ЗРЕНИЯ УМЕНЬШЕНИЯ ШУМА В САМОМ ИСТОЧНИКЕ
В.В. СГАДЛЕВ
Статья представлена доктором технических наук, профессором Елисовым Л.Н.
В статье рассмотрен пример оптимизации классической ступени вентилятора, с точки зрения уменьшения шума в самом источнике.
Интенсивное развитие авиационной техники, создание нового поколения пассажирских и транспортных самолетов потребовали решения ряда сложных научно-технических задач в области аэродинамики, прочности, двигателестроения, приборостроения, материаловедения и технологии. Проблема снижения шума самолетов, являющаяся частью общей проблемы охраны окружающей среды, приобрела важное значение в авиации. Вот почему акустические характеристики пассажирских самолетов стали одним из показателей, определяющих их конкурентоспособность. С целью ограничения роста авиационного шума разработаны стандарты, соблюдение которых является необходимым условием эксплуатации пассажирских самолетов. Поэтому очень важной практической проблемой является проблема расчета акустических характеристик современных двигателей.
Реактивный двигатель является сложным источником шума, поскольку шум образуется во всех его узлах: компрессоре, камере сгорания, турбине, реактивном сопле, реактивной струе газов. Кроме того, значительным источником шума является вентилятор.
В ЦИАМ разработан аналитический метод, предназначенный для анализа спектрального состава тонального шума лопаточных машин и выводящий результаты анализа в графическом представлении, основанный на теории Тайлера-Софрина. Метод реализован в виде программы, на языке Visual С++, работающей под Microsoft Windows. Программа имеет удобный пользовательский интерфейс и обеспечивает графическую визуализацию полученных результатов.
В статье рассмотрен пример классической ступени вентилятора, ступени с двумя вращающимися в разные стороны рабочими колесами.
Математическая модель Моды в канале
Возмущения, распространяющиеся по прямолинейному цилиндрическому каналу, как простому (трубе), так и коаксиальному, делятся на пять типов.
Энтропийное возмущение плотности.
Два типа вихревых возмущений скорости.
Два типа звуковых возмущений.
В свою очередь, звуковые возмущения, в зависимости от характеризующих их параметров, могут быть разделены на две категории:
Затухающие.
Незатухающие.
Из этих возмущений энтропийные и вихревые возмущения распространяются только в направлении потока. Звуковые возмущения могут распространяться как по потоку, так и против него.
Возмущения в канале описываются в виде конечной суммы частных решений уравнения Эйлера (мод) в цилиндрической системе координат, характеризуемых следующими параметрами. т - ёбоаТ аау +апоТ оа т - аде! ооаёш1 а -^ёПёТ п - бааёаёш 1 а ^епёТ а - а! 1 ёёооаа -оадт аие пааёа
Приведем выражение для тех частных решений линеаризованных уравнений Эйлера в цилиндрической системе координат {х, г,в}, которые описывают распространение звуковых возмущений. В качестве параметров пульсаций выберем плотность р, скорость в направлении оси канала - и , радиальную компоненту скорости у'г и угловую скорость м>'в . Пусть канал имеет радиус Е00, среднее стационарное поле течения в нем однородно, и скорость течения и0 направлена вдоль оси X. Тогда можно показать, что в этом случае имеем следующее выражение для частного решения, определяемого указанным набором параметров:
( 1 Л
~ Ут (1тпГ )
(р л /
и
у:
ч
V р J
= аЫт
к,.
пр0
1
Ут (1тпГ )
/
Ут (1тп
^ ттпр0 т Ут 1тпГ )
^ютпр0 Г
Ут (1тпГ )
где Ут (2) - функция Бесселя (2 - произвольный аргумент). В данном выражении использованы сокращения
к,
т -М ± А,
-т
1 - М2 -1 ± I А'
________,
1 - М2
У
где с - скорость звука.
В вышеприведенных уравнениях плюс надо взять для возмущения, распространяющегося в положительном направлении оси X, а минус - для возмущения, распространяющегося в отрицательном направлении. Нормировочный фактор Мттп в программном комплексе определяется
)
Г
с
таким образом, чтобы максимальное значение давления было равно амплитуде. Параметр 1тп
/
определяется граничным условием Ут (1тпК0 ) = 0 . Это уравнение имеет бесконечное число корней, поэтому радиальное число п определяет номер корня уравнения, которому соответствует данное конкретное значение 1тп.
Отметим, что выражение (1) описывает распространение незатухающей звуковой волны только при действительном А,тп , т. е. при:
1 -(1 -м2)V,] >0.
В случае, если данное условие распространения не выполняется, волна в канале экспоненциально затухает. Расстояние, на котором волна затухает в два раза, равно
ln (2) с 1 - M2
exp О (і m2)f1mnC"] 1
V V О У
Предположим теперь, что волна распространяется по коаксиальному каналу с внешним радиусом Rmax и внутренним радиусом Rmln. В этом случае, вместо функций Бесселя в выражении
для частного решения следует использовать функции следующего вида:
/ /
Gm ( Rmln ,1mn , r ) = Nm ( KmnRmn ) Jm (Kj )- Jm ( KmnRmn ) Nm (1mnr ) ,
где Nm (z) - функция Неймана.
Функция Gm (Rmln, 1mn, r) автоматически удовлетворяет граничному условию на внутреннем радиусе. Граничное условие на внешнем радиусе - Gm (Rmln, 1тп, Rmax) = 0 как и в случае цилиндрического канала, позволяет определить значение ln .
Таким образом, в случае, если граница расчетной области соответствует сечению простого цилиндрического или коаксиального каналов, среднее поле течения на границе однородно и отсутствует завихренность, можно задать условие на границе расчетной области в виде набора мод, описываемых выражением (І).
Соотношения Тайлера-Софрина
Подробный анализ мод, возникающих при ротор-статор взаимодействии, был проведен в бО-х годах Тайлером и Софриным. Для однорядного вентилятора они показали, что возникают только моды, азимутальные числа которых связаны с числами лопаток ротора В и статора V следующим соотношением:
m = nB + sV,
где n номер гармоники, n>0 , а s = -»L« , где s целое число.
При этом w = nBW.
Заметим, что угловая скорость вращения ротора Q может быть как положительной, так и отрицательной величиной, поэтому и частота моды может иметь произвольный знак. Для удобства перепишем данные соотношения таким образом, чтобы частоты мод лежали в положительной полуоси, а направления их вращения определялись бы знаком m. Нетрудно показать, что в этом случае
(О = nB |Q| m = ( nB + sV ) sign ( W ) .
Для биротативного вентилятора необходимо потребовать, чтобы соотношения Тайлера-Софрина выполнялись как в системе координат, связанной с первым ротором, так и в системе координат, связанной со вторым. После несложных преобразований легко показать, что
О = \пВ10.1 + sB2W2 \;
т = (пБ1 + sB2) бій g (пБ10.1 + &В202); п = 0;
где изменение п ограничено неотрицательными числами, чтобы избежать дублирования мод. Здесь Б1 и В2 числа лопаток, а 01 и ^2 - окружные скорости для первого и второго роторов соответственно.
Анализ полученных данных позволяет, при правильном выборе режимных и геометрических параметров вентилятора, организовать отсечку неблагоприятных звуковых частот, при сохранении высоких аэродинамических характеристик вентилятора. Метод может быть использован на начальном этапе проектирования, при выборе основных режимных и геометрических параметров вентилятора.
Методика оценки зависимости шума вентилятора от числа лопаток ротора и статора
Так как одним из основных источников шума двигателя является вентилятор, то снижение шума вентилятора очень важная проблема. Этот шум зависит от режимных и геометрических параметров вентилятора, а также числа лопаток ротора и статора для классической схемы вентилятора, или числа лопаток обоих роторов для биротативной схемы вентилятора.
Обычно режимные и геометрические параметры проточной части заданы заранее и их изменение не рекомендовано. Поэтому мы будем исследовать зависимость шума вентилятора от соотношения количества лопаток ротора и статора для классической схемы вентилятора или соотношения числа лопаток роторов для биротативной схемы вентилятора.
Чтобы исследовать эту зависимость, мы предполагаем, что звуковая мощность распределена равномерно по группам мод с одинаковыми азимутальными числами. Следовательно, в этом приближении шум в заданном диапазоне частот определяется количеством мод в данном диапазоне. Для усовершенствования методики планируется ввести весовые коэффициенты моды для лучшего учета человеческой физиологии слуха.
Результаты расчета
Мы предполагаем, что звуковая мощность распределена равномерно по группам мод с одинаковыми азимутальными числами. Следовательно, в этом приближении шум в заданном диапазоне частот определяется количеством мод в данном диапазон.
Начальными данными для работы по оптимизации количества лопаток классической ступени вентилятора были выбраны параметры одной из модельных ступеней, исследовавшихся в ЦИАМе:
Число Маха = 0.339;
Скорость звука = 336м/с;
Частота вращения ротора = 1400.69рад/с;
Внешний радиус канала = 0.2м;
Внутренний радиус канала = 0.03м;
Число лопаток рабочего колеса = 18;
Число лопаток статора = 30.
Общий вид классической ступени вентилятора и схемы расположения венцов приведены на рис. 1 и рис. 2 соответственно.
Оптимизация проводилась в диапазоне частот от 0 до 18 кГц потому, что эта модель является модельной и диапазон частот от 0 до 18 кГц примерно соответствует диапазону частот от 0 до 9 кГц полноразмерного вентилятора. Взят диапазон азимутальных чисел от -115 до 115, так как этот диапазон удовлетворяет всем модам, которые входят в диапазон частот от 0 до 18 кГц.
Для начала работы над ступенью надо рассмотреть, что же излучает не оптимизированная ступень вентилятора (рис. 3). Как мы видим на рисунке, в данных диапазонах частот и азимутальных чисел присутствует 10 излучаемых мод.
Процесс оптимизации проходил следующим образом, мы изменяли количество лопаток рабочего колеса и находили оптимальное для него количество лопаток статора. Проводили несколько итераций и сравнивали найденные соотношения
Далее представлены наиболее интересные нам варианты.
Первый из них имеет следующую конфигурацию:
Количество лопаток первого ротора - 18.
Количество лопаток второго ротора - 34.
Количество излучаемых мод представлено - 6.
Излучаемые моды представлены на рис. 4.
Второй вариант получен увеличением количества лопаток рабочего колеса на 3 и уменьшением лопаток статора на 1, при этих изменениях получили то же самое количество излучаемых мод равное 6, как и в предыдущем варианте. Но этот вариант интересен тем, что предполагает снижение количества лопаток статора (рис. 5).
Последний вариант был получен из предыдущего путем добавления десяти лопаток к статору, и получили снижение количества излучаемых мод до 3 (рис. 6).
В результате оптимизации количество лопаток на первом рабочем колесе увеличено на 3, а на втором рабочем колесе - на 9, что позволило снизить количество излучаемых мод с 10 до 3.
Дальнейшая оптимизация была прекращена из-за соображений большого увеличения веса ступени вентилятора.
Расчет спектрального состава тонального шума вентилятора проводился с помощью программы, написанной автором статьи на языке УІБиаІ С++, работающей под Microsoft Windows. Результат был получен путем анализа накопленных аналитических решений, основанных на теории Тайлера-Софрина и теории распространения акустических волн в канале.
Выводы
В работе представлен метод оценки акустических характеристик вентилятора. Метод основан на получении аналитических решений теории Тайлера-Софрина и теории распространения возмущений в канале. Разработана программа для получения этих решений и представления их в удобной для пользователя форме.
После проведенной верификации программы была проведена оптимизация числа лопаток модельной ступени вентилятора, исследовавшаяся в ЦИАМе, после ее оптимизации количество мод уменьшилось на 30%, что предполагает соответственное снижение уровней шума.
Таким образом, предложенный метод оптимизации числа лопаток с нашей точки зрения является перспективным.
Стоит отметить, что данная методика будет более совершенной, если в нее ввести весовые коэффициенты моды для лучшего учета человеческой физиологии слуха и более оптимальной настройки характеристик ЗПК.
Рис. 1. Трехмерная математическая модель вентилятора
X
Рис. 2. Схема расположения венцов
18000
17067
Ь ЗЬ
14^/0
240Ь
1147*(Нг)
10540
Рис. 3. Моды, излучаемые не оптимизированной ступенью вентилятора
Рис. 4. Моды, излучаемые после первого шага оптимизации
Рис. 5. Моды, излучаемые после второго шага оптимизации
Рис. 6. Моды излучаемые оптимизированной ступенью вентилятора
ЛИТЕРАТУРА
1. Tyler J.M., Sofrin J.E. Axlal Flow Compressor Nolse Studles. - SAE Trans., vol. 70, 19б2, pp. 309-332.
2. Нюхтиков, М.А., Россихин, А.А., Сгадлев, В.В. Численное моделирование тонального шума вентиляторов современных двигателей // Научный Вестник МГТУ ГА, серия Эксплуатация воздушного транспорта и ремонт авиационной техники. Безопасность полетов, №10Б(11), 200б. - С. 23-29.
BLADES NUMBER OPTIMIZATION OF THE CLASSICAL FAN STAGE, FROM THE POINT OF VIEW OF NOISE REDUCTION IN THE SOURCE
Sgadlev V.V.
In the article we presented the example of optimization for the classical fan stage, from the point of view of noise reduction in the source.
Сведения об авторе
Сгадлев Василий Викторович, 1983 г.р., окончил МГТУ ГА (2005), аспирант МГТУ ГА, автор 3 научных работ, область научных интересов - численное моделирование акустических характеристик.
