CC BY
201
УДК 621.438
ПРОФИЛИРОВАНИЕ РАБОЧИХ КОЛЕС РАДИАЛЬНО-ОСЕВЫХ ТУРБИН
С ПОМОЩЬЮ КРИВЫХ БЕЗЬЕ
©2011 О.В.Батурин
Самарский государственный аэрокосмический университет имени академика С.П. Королёва
(национальный исследовательский университет)
В работе описывается алгоритм построения профилей рабочих лопаток радиально-осевых турбин с помощью кривых Безье. Поверхность лопаток, полученных по данному алгоритму, имеет монотонно меняющуюся кривизну, что благоприятно скажется на КПД. Кроме того, лопатки радиальны во всех сечениях, перпендикулярных оси вращения, что позволяет существенно снизить изгибные напряжения в пере лопатки.
Радиально-осевая турбина, профилирование, меридиональное сечение, рабочая лопатки, кривая Безье.
Общий уровень эффективности турбины определяется выбранным сочетанием основных режимных параметров, таких как параметр нагруженности Ут, степень реактивности рст и т.д. Кроме того, КПД турбины зависит от потерь в проточной части, которые определяются аэродинамическим совершенством межлопаточного канала. Оно обеспечивается на этапе профилирования проточной части (ПЧ). Задача этого этапа заключатся в том, чтобы найти такую форму лопаток, которая обеспечивает требуемый поворот потока с минимальными потерями при сохранении приемлемой прочности и технологичности изготовления. Грамотно проведенное профилирование позволяет достигнуть, а в ряде случаев и превысить, значение КПД, заложенное при газодинамическом проектировании.
Процесс профилирования радиально -осевой турбины условно можно разделить на несколько этапов:
1) построение меридионального обвода проточной части;
2) построение средних линий профилей в контрольных сечениях;
3) построение профилей лопаток в контрольных сечениях.
Профилирование лопаток РК радиальной турбины целесообразно проводить в цилиндрической системе координат (сшг) (рис. 1), центр которой находится на оси вращения, осевое направление оа - совпадает с осью вращения РК, а радиальная ось ог проходит через точку на внутреннем обводе проточной части с максимальным значением радиуса. Окружное направление он перпендикулярно двум другим осям.
используемая при профилировании
На первом этапе проектирования РК радиально-осевой турбины формируются меридиональные обводы проточной части -обводы в плоскости, проходящей через ось вращения ротора {гоа, рис. 1). Они во многом определяют распределение скоростей по длине межлопаточного канала и в конечном итоге - уровень потерь в РК.
Меридиональное сечение состоит из втулочного (внутреннего) и периферийного (внешнего) обводов, ограничивающих лопатку со стороны диска и сверху, а также промежуточных сечений.
Для обеспечения благоприятного распределения скоростей и снижения потерь изменение как внешних, так и внутренних обводов меридионального сечения должно быть плавным, без резкого изменения кривизны. Радиусы кривизны должны быть максимальными, а на границах обводов (входе и выходе РК) - стремиться к бесконечности (в этом случае контур обводов приближается к прямой линии).
В доступной технической литературе существует множество подходов к профилированию меридиональных обводов. Широко распространен подход, согласно которому они формируются дугами окружности. Внутренний обвод может строиться независимо от наружного в виде дуги, величина радиуса которой задается некой полуэмпирической формулой, но чаще рекомендуется находить внутренний обвод как огибающую окружностей, касающихся наружного обвода, диаметр которых плавно уменьшается от выхода к входу. Построение меридиональных обводов в виде дуг является самым простым, но не самым лучшим с точки зрения минимизации потерь в проточной части. В ряде работ контур наружного обвода мери-
дионального сечения рекомендуется описывать полиномом третьей степени.
В некоторых зарубежных публикациях есть упоминание о построении элементов РК с помощью кривых Безье. Данный подход видится перспективным, поскольку кривая Безье имеет важное свойство - кривизна меняется плавно, непрерывно и без разрывов на всем протяжении.
Кривая Безье строится на основе опорных точек. Двумя из них являются концы кривой, а другими могут быть произвольные точки в любом количестве. Для профилирования элементов проточной части радиаль-но-осевых турбин целесообразно применять кривые с тремя или четырьмя опорными точками (рис. 2).
а б
Рис. 2. Примеры построения кривых Безье с тремя (а) и четырьмя (б) опорными точками
В данной статье приводится общий подход к построению меридиональных обводов радиально-осевой турбины с помощью кривых Безье на основе трех опорных точек. Две точки являются концами кривой, а третья опорная точка находится в месте пересечения касательных, проведенных к кривой в концах.
Построение втулочного и периферийного меридионального обвода РК производится на основе данных, полученных в результате газодинамического расчета. Исходными данными для построения меридионального профиля являются: наружный диаметр РК на входе Щ], втулочный, средний и внутренний диаметры на выходе из РК И^вт, Г)2ср-, 1-Ьк-, его ширина В и высота лопатки на входе Ъ]. Эти данные определяют положение точек А, В, С, Д являющихся крайними (опорными) для кривых, описывающих втулочный и периферийный обводы (рис. 3).
В случае, если кривая имеет три опорных точки, то для каждого обвода следует
найти по одной недостающей точке (/• - для внутреннего обвода и Е - для наружного контура). Они лежат на пересечении касательных, проведенных к кривой в крайних точках (А, I) - для внутреннего обвода и В,С - для наружного контура).
Для того, что найти положение точки /•, через точку А проводится прямая, наклоненная под углом в]6Ш к вертикали, а через точку В - прямая под углом 9дт к горизонтали. Аналогично точка Е находится в месте пересечения прямой, проходящей через точку В под углом в]К к вертикали, и прямой, проходящей через точку С под углом 0\: к горизонтали (рис. 3).
Примечание. Положительное направление угла Оц соответствует повороту вертикальной прямой против часовой стрелки. Положительное направление угла 0;: соответствует повороту горизонтальной прямой по часовой стрелке.
а
Рис. 3. Построение меридиональных обводов с помощью кривых Безъе с тремя опорными точками
Величина угла Оц обычно находится в интервале Оц =0... 10а величина угла Ö2i находится в диапазоне 02i =10... 15°.
Находятся координаты недостающей базовой точки F верхнего обвода. Для этого на основе выбранных ранее углов наклона касательных в крайних точках записываются уравнения прямых, проходящих через точки AFnFD:
уравнение прямой AF
r = tg(eie„,+90)tf + A,M;
Уравнение прямой FD:
= tgИ:„„)а + /Л„,„ -Вtg{-в2т1), м.
Зная, что эти прямые имеют общую точку F, приравняв их, можно найти значение осевой координаты, соответствующей этой точке:
Аналогичным образом находятся координаты дополнительной базовой точки верхнего меридионального обвода.
Записываются уравнения прямых ВЕ:
г = Щ (0„. + 90). а Щ - Щ (0„. + 90) Ъх, м и ЕС.
г = (~е2„ )а + 1К-В-Щ (~е2вт), м.
Находится радиальная координата точ-
ки Е:
аЕ =
м.
аF =
м.
Найденное значение осевой координаты подставляется в уравнение одной из прямых АР или ЕР, и находится значение радиальной координаты точки Е:
р2к-в-гё(-в2к)-р^гё(01к + 90)
1ё(01к+90)-1ё(-02к) Затем найденное значение осевой координаты подставляется в уравнение одной из прямых ВЕ или ЕС и находится значение радиальной координаты точки Е:
гЕ = Хё(еы +90)аЕ +Д -1§(б1к +90)6,,м.
Таким образом, на данном этапе известны радиальные гг и осевые аг координаты всех трех опорных точек периферийного и втулочного обводов. Это позволяет записать уравнения кривых Безье.
Для наружного (периферийного) обвода:
а.
= (1 - t2)aB + 2f(l - t)aE + ¿2ас;
г1к =(1-Г)Гв + 2({\-()Гб + ГГс.
Для внутреннего (втулочного) обвода: ат, = С1 -12)аА +2/(1-/>/,,+ 12ав;
В этих уравнениях параметр t меняется в интервале 0 <1 < 1.
Для построения кривых следует задаться п значениями параметра / из интервала 0 < 1 < 1. Удобнее всего будет, если интервал [0,1] будет разбит на одинаковые отрезки. Для корректного построения число п должно быть достаточно большим (например, 20). Число точек на верхнем и нижнем обводе должно совпадать.
По аналогичному принципу строятся меридиональные сечения корпуса над лопатками и произвольное число сечений лопатки по высоте ПЧ. Промежуточные сечения могут строиться двумя путями. Первый полностью повторяет тот, по которому строились наружные обводы. При этом величины наклона касательных 0, необходимые для построения кривых Безье, находятся по значениям ввт и 9К, принимая линейное изменение угла по высоте ПЧ. Второй путь основан на делении отрезков, соединяющих одноименные точки (с одинаковой величиной 1) наружного и внутреннего обвода в пропорции, соответствующей высоте требуемого сечения.
На рис. 4 приводится пример спрофилированного по приведенному выше алгоритму меридионального сечения проточной части радиальной турбины регистровой системы наддува дизеля Д21-26ДГ,
Качество проектирования меридионального обвода РК предлагается контролировать на основе зависимости высоты лопатки по длине проточной части. Высота лопаток в одноименных точках наружного и внутреннего обвода (в точках с одинаковыми значениями ) находится с помощью следующей формулы:
/ * (а -а. У + (/• — г V,л/.
1 у V ¡к ют / \ 1К хш /
Зная длину средней линии канала, строится упомянутый выше график (рис. 5). Критерием качества проектируемого канала является плавный и монотонный рост высоты лопатки без резких скачков.
Л , лыв
120 -100 -ВО -60 -40 -20 -О -
О 20 40 60 й, ЛЫИ
Рис.4. Спроектированный меридиональный профиль
РКрадиально - осевой турбины для регистровой системы наддува дизеля Д21-26ДГ
ВЫХОД
ВХОД
О 5 10 15 20
Рис. 5. Зависимость высоты лопатки по длине проточной части, используемая для оценки качества профилирования меридиональных обводов
Следующий шаг по профилированию рабочих лопаток заключается в отыскании формы их средних линий в контрольных сечениях. Средняя линия задается в виде набора точек в цилиндрической системе координат (а, г, <р). Меридиональная проекция средних линий была построена в предыдущем разделе и, строго говоря, координаты а, г точек, лежащих на контрольных сечениях, уже известны.
Важным требованием, предъявляемым к турбине - прочность лопаток. Для уменьшения изгибных напряжений рабочие лопатки будут проектироваться таким образом, что лопатка в любом осевом сечении будет радиальной. Это достигается за счет того, что угловая координата средних линий <р будет одинакова для всех сечений.
Кроме того, рабочую точку турбины предпочтительно выбирать таким образом,
чтобы, рабочие лопатки во входном сечении были направлены радиально (/3/л=90°>.
Таким образом, задача построения средней линии в любом сечении сводится к отысканию ее угловой координаты в контрольном сечении. В качестве него было принято среднее сечение, параметры в котором были рассчитаны в проектном расчете.
Рис. 6. К поиску угловой координаты средней линии профиля
Для нахождения угловой координаты ср строится проекция средней линии на окружное направление. Оно представляет собой изогнутую кривую (рис. 6). Ее окружная протяженность равна шагу между лопатками в контрольном сечении (С]Ъ а осевая длина равна ширине рабочего колеса В. Для обеспечения высокого аэродинамического качества лопатки окружная проекция средней линии строится с помощью кривой Безье, поскольку последняя имеет плавное изменение кривизны на всем протяжении. Как и при построении меридиональных обводов РК для формирования кривой Безье использовались три опорных точки: точка начала кривой К, лежащая на осевой линии, точка В и точка М пересечения осевой линии и прямой проведенной через точку В под углом А?лсР от вертикали (рис. 6). Расстояние, на которое точка К отстоит от точки Т, обычно находится в диапазоне 0,25...0,3В.
Поиск координат данных точек не составляет труда и осуществляется по алгоритму, описанному в предыдущем разделе. Таким образом находятся значения окружных координат и, которые, опираясь на значение радиуса контрольного сечения гср, пе-
ресчитываются в значение угловой координаты:
На следующем этапе рассчитываются координаты точек, лежащих на поверхности спинки и корытца лопатки в контрольных сечениях.
Построение происходит в следующей последовательности:
Шаг 1. Выбирается толщина лопатки во втулочном сечении на входе 8] и угол заострения лопатки на входе // (рис. 7).
Шаг 2. Выбирается толщина лопатки во втулочном сечении на выходе <5? и угол заострения лопатки на выходе у? (рис. 7).
Рис. 7. Задание толщины лопатки
Величина угла заострения обычно находится в диапазоне <5=0,5...2°, а толщина лопатки во втулочном сечении в первом приближении может быть принята равной 8=2... 10 мм в зависимости от размерности турбины. При этом следует следить, чтобы толщина профиля на периферии по технологическим соображениям не была меньше 1..1,5мм.
Шаг 3. По величине углов заострений и толщины профиля на втулке находится толщина лопаток в контрольных сечениях, пологая, что она меняется линейно.
Шаг 4. Полагая, что толщина лопатки меняется линейно от входа к выходу, находится величина полутолщины 8/2 в разных точках по длине контрольных сечений.
Шаг 5. Зная радиальную координату точки, которой соответствует значение полутолщины, находится разница угловых координат поверхности лопатки и средней линии:
А ф: =-град.
л -т;.
Рис. 8. Внешний вид спроектированного рабочего колеса радиально-осевой турбины для регистровой системы наддува дизеля Д21-26ДГ
Шаг 6. Находятся угловые координаты поверхности спинки и корытца:
<Рс, =(Рсрп + Ь(р„град;
Задняя поверхность диска также строится на основе кривой Безье с тремя опорными точками.
По описанной выше методике и на основе проведенного проектного расчета была
построена полная геометрия РК спроектированной радиально-осевой турбины для регистровой системы наддува дизеля Д21-26ДГ. Его внешний вид показан на рис.8.
Работа выполнена при финансовой поддержке Правительства Российской Федерации (Минобрнауки) на основании постановления Правительства РФ №218 от 09.04.2010.
BLADING OF RADIAL-AXIAL TURBINE WHEEL USING BEZIER CURVE
©2011 О. V. Baturin
Samara State Aerospace University named after academician S.P. Korolyov
(National Research University)
The paper describes the constructing algorithm for radial-axial turbine rotor blades using Bezier curves. The surface of the blade obtained by this algorithm is inonotonically changing curvature and it impact positively on efficiency. In addition, blades are radial in all section perpendicular to the axis of rotation; it can significantly reduce the bending stresses in theblade body.
Radial turbines, blading, meridional section, rotor blade, Bezier curve.
Информация об авторах
Батурин Олег Витальевич, кандидат технических наук, доцент кафедры теории двигателей летательных аппаратов Самарский государственный аэрокосмический университет имени академика С.П. Королёва (национальный исследовательский университет). Е - mail: oleg.v.baturin@gmail.com. Область научных интересов: рабочие процессы в лопаточных машинах, вычислительная газовая динамика, агрегаты наддува ДВС.
Baturin Oleg Vitalievich, Cand idate of Technical Science, associate professor of the chair of theory of engine for flying vehicle, Samara State Aerospace University named after academician S.P. Korolyov (National Research University). E - mail: oleg. v.baturin@gmail.com. Area of research: workflows in turbomachines, computational fluid dynamics, turbocharger.
