Компрессорный профиль С4: автоматизация проектирования

Вестник РУДН. Серия: Инженерные исследования

RUDN Journal of Engineering researches

2018 Vol. 19 No. 2 155-164

http://journals.rudn.ru/engineering-researches

DOI 10.22363/2312-8143-2018-19-2-155-164 УДК 621.165.533

КОМПРЕССОРНЫЙ ПРОФИЛЬ С4: АВТОМАТИЗАЦИЯ ПРОЕКТИРОВАНИЯ

Л.В. Виноградов, В.К. Мамаев, П.П. Ощепков

Российский университет дружбы народов Российская Федерация, 117198, Москва, ул. Миклухо-Маклая, 6

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

Ключевые слова: компрессор, профиль С4, решетка, САПР, нелинейные преобразования

Сложные обводы и поверхности получили широкое распространение в лопаточных машинах (компрессоры и турбины). Проектирование и расчет таких обводов и поверхностей с применением современной вычислительной техники, а также воспроизведение на станках с программным управлением требует их корректного математического описания [1]. Одним из современных методов конструирования поверхностей является метод нелинейных преобразований, который позволяет получать искомую поверхность (обвод) как образ плоскости или какой-либо известной кривой.

Введение

Данные и методы

В работе применено преобразование Гирста !2 с несобственным центром F'J) [1]. У рассматриваемой специализации кремоновых инволюций ^ с несобственным центром F*1> есть очень важное метрическое свойство: такое преобразование является эквиформным, т.е. сохраняются площади прообраза и образа.

В использованной специализации квадратичной инволюции !2 инвариантной коникой является квадратичная парабола d2 в виде кривой Ве21ег 2-го порядка.

Прообразом компрессорного изогнутого профиля является симметричный исходный компрессорный профиль С4 [2; 3], аэродинамические характеристики которого известны из экспериментальных данных.

Все вычисления и построения выполнены в интегрированном математическом пакете MathCAD [5], для чего была разработана система автоматизированного проектирования (САПР) и подпрограмма для расчета координат опорных точек полиномов Ве21ег 2—5-го порядков (Pp — «Опорные точки»).

Объектом построения был профиль С4 [2], симметричный, с максимальной относительной толщиной 10%, хордой профиля 100 мм, относительная абсцисса сечения с максимальной толщиной — 30%, радиус входной кромки г1 = 1,2 мм, радиус выходной кромки г2 = 0,6 мм. Ниже приведена матрица координат обвода исходного профиля, где первая строка — абсциссы, вторая — ординаты.

Матрица координат компрессорного профиля С4

0 1,25 2,5 5,0 7,5 10 15 20 30 40 50 60 70 80 90 95 100

0 1,65 2,27 3,08 3,62 4,02 4,55 4,83 5,0 4,89 4,57 4,05 3,37 2,54 1,60 1,06 0

На рисунке 1 показан векторный график исходного профиля. График строится для получения первичного образа (прообраза) объекта проектирования и контроля корректности исходных данных и их ввода.

Рис. 1. Векторный график симметричного профиля С4:  — верхний обвод профиля;  — нижняя половина обвода профиля [Fig. 1. Vector plot of symmetric C4 blade profile:  — top part contour;  — bottom part contour]

Результаты вычислений

Для последующих вычислений векторы обводов исходного профиля были интерполированы кубическими сплайнами, для чего в пакете MathCAD [5] имеется стандартная функция, позволяющая работать как с аналитическими функциями.

На рисунке 2 VSvxxx, VSvxy — векторы абсцисс и ординат спинки профиля, x1(@), y1(@) — параметрическое уравнение окружности входной кромки.

Рис. 2. Область входной кромки профиля:  — верхний обвод профиля;

♦ — нижняя половина обвода профиля;--входная кромка профиля

[Fig. 2. Leading edge of the blade profile:  — top part contour; ♦ — bottom part contour;--leading edge]

Для начала аппроксимации вариантными расчетами были определены координаты начальной точки (абсцисса и ордината) А:

XA = 0.2YA(x) = ,Jr12 -(x -r1 )2 YA = YA(XA)YA = 0.663325.

Начальная точка А расположена на окружности входной кромки, и аппрокси-мационная кривая должна быть касательной к этой окружности. Эта точка является одной (левой) опорной точкой характеристического треугольника аппрок-симационной кривой Bezier. Правая опорная точка С характеристического треугольника определялась методом последовательного приближения — перебором точек матрицы координат профиля. Для первого аппроксимационного участка правой опорной точкой стала точка вектора с индексом 9.

Для расчета вектора опорных точек полиномов Bezier применялась разработанная подпрограмма

Pp_basic_point(XA, xC, YA, yC, aA, aC, 2).

Формат матрицы исходных данных для этой подпрограммы имеет вид

Pp_basic_point(XA, xC, YA, yC, aA, aC, 2),

где — XA, YA, xC, yC — координаты правой и левой (конечных) опорных точек кривой Bezier; аА, аС — углы наклона касательных в этих точках.

Подпрограмма формирует матрицу координат опорных точек полиномов Bezier и параметрические уравнения соответствующих полиномов.

Ниже приведена матрица координат опорных точек для профиля С4.

Матрица координат опорных точек полиномов Bezier

V1 =

' 0.2 3.528 40 0.2 0.663 1.032 1.918 6.423 3.073"

0.663 5.681 4.89 1.864 3.172 11.814 4.229 21.348 4.658

0 0 0 21.764 5.285 30.882 5.088 0 0

V 0 0 0 40 4.89 0 0 0 0

_......,- -—-—'

л

0

10

20

30

Рис. 3. Начальный участок аппроксимации обвода спинки:--аппроксимационная кривая Безье;

— исходные точки контура профиля С4

[Fig. 3. Initial region of approximation for blade back contour:--Bezier approximation curve;

— points of original C4 blade profile contour]

На рисунке 3 показан первый (начальный) участок аппроксимации.

Второй аппроксимационный участок начинается с элемента вектора с индексом 9 (последний элемент первого аппроксимационного участка) до 14. Крайний правый элемент второго участка определялся методом последовательного подключения крайнего элемента к вектору характеристического треугольника второго участка. Координаты опорных точек кривой Bezier 2-го участка рассчитываются по подпрограмме, которая ранее уже использовалась. На рисунке 4 показан второй участок аппроксимации первичного профиля.

y22bz

Vy ♦ ♦♦

♦ '

♦ ♦

♦ ♦

20

40 60

x22, Vx

80

100

Рис. 4. Второй аппроксимационный участок спинки профиля, построенный по уравнению y22bz(x22) [Fig. 4. Second approximation region for blade back contour constructed on the basis of y22bz(x22) equation]

Третий заключительный участок аппроксимации включал элемент с индексом 14 (последний элемент второго участка) и окружность выходной кромки. Определение параметров правой опорной точки характеристического треугольника было выполнено методом последовательного приближения, задавая абсциссу точки на окружности выходной кромки. Для профиля С4 абсцисса этой точки была определена как xC = 94,7.

По аналогии с предыдущими интервалами координаты опорных точек полиномов были определены по подпрограмме Pp — «Опорные точки».

После аппроксимации трех участков были записаны уравнения обвода спинки профиля С4 в виде

^C4(x)_

^/r2 - (х - r1)2 if 0 < x < 0.2 y2bz(x) if 0.2 < x < 40 y22bz (x) if 40 < x < 90 .

y222bz(x) if 90 < x < 99.47

^r22-Lx-(b-r2)J2 if 99.47< x< b

6

0

В силу симметрии уравнение корытца можно записать в виде

Же4(х) = -^(х).

На рисунке 5 показан симметричный исходный компрессорный профиль С4, аппроксимированный двумя окружностями входной и выходной кромок, тремя кривыми Ве21ег 2-го порядка и построенный по вышеприведенным уравнениям.

5 0 -5

0 20 40 60 80 100

Рис. 5. Профиль С4, построенный двумя окружностями входной и выходной кромок и тремя кривыми Безье 2-го порядка [Fig. 5. C4 blade profile constructed by two circles of leading and trailing edges and three second-order Bezier curves]

Построение изогнутого профиля методом нелинейных преобразований было выполнено по уравнениям вида

Snl

WNL

где LmidC4(x) — инвариантная кривая 2-го порядка в виде кривой Bezier с углом на входе х1 = 14°, на выходе х2 = -10°, т.е. угол поворота потока в решетке 0 = 24°, что соответствует практически максимальному углу поворота в компрессорной решетке [2].

Следует отметить, что при изгибе симметричного профиля инвариантная кривая является средней линией изогнутого профиля.

Параметры характеристического треугольника для инвариантной кривой были рассчитаны по подпрограмме расчета опорных точек. На рисунке 6 представлены изогнутый профиль и инвариантная кривая — средняя линия профиля.

5 0 -5

0 20 40 60 80

Рис. 6. Фрагмент изогнутого компрессорного профиля:

... — инвариантная кривая (средняя линия);--обвод профиля

[Fig. 6. Fragment of curved compressor blade profile: ... — invariant curve (centreline);--profile contour]

На рисунке 7 показана решетка компрессорных профилей типа С4 с углом установки профиля в решетке у = 70°, относительным шагом t = 1.

Параметрические уравнения координат изогнутого профиля решетки имеют

х'(x) = x- cosI 70+ SNL(x)• sinI 70--П-

1 180J NL I 180

y' (x) = -x-sinI 70| + SNL(x)-cosI 70

y I 180 J NLX ' | 180

x'w(x) = x• cosI 70+ WNL(x)-sinI 70•-П-

V ; 1 180 J I 180

y'w(x) = -x-sinI 70+ WNL(x)-cosI 70--П-l 180 J ^ 180

у '

у 'w

20

-20

-40

-60

-80

-100

50 100

x' (x), x' w(x), x' (x) + f, x'w(x) + t

150

Рис. 7. Решетка изогнутых компрессорных профилей типа С4

с углом установки у = 70° и относительным шагом t = 1 [Fig. 7. Cascade of curved C4 compressor blade profiles with angle of incidence у = 70° and pitch ratio t = 1]

Аналитические зависимости обводов профиля удобны для расчета геометрических характеристик профиля, которые используются в прочностных расчетах [2]. В таблице представлены наименования, формулы, имена по программам и численные значения геометрических параметров изогнутого профиля С4 в исходной системе координат.

На рисунке 8 представлен изогнутый компрессорный профиль С4 с главными центральными осями.

0

0

Таблица

Геометрические параметры изогнутого профиля С4 [Table. Geometric parameters of curved C4 compressor blade profile]

№ п/п Наименование параметров Имя Значение

1 Площадь профиля [Profile area] b F = j(WNL (x) - SNL (x))dx 0 F F = 718.926

2 Статический момент инерции площади относительно оси x [First moment of area relative to x-axis] b WNL{x) Sx = J J sdsdx 0 SNL(x) Sx Sx = 2.805  103

3 Статический момент инерции площади относительно оси y [First moment of area relative to y-axis] b WNL(x) Sy = J J xdsdx 0 SNL(x) Sy Sy = 3.117  104

4 Момент инерции площади относительно оси x [Second moment of area relative to x-axis] b WNL(x) Jx =J J s2 dsdx 0 SNL( x) Jx Jx = 1.635  104

5 Момент инерции площади относительно оси y [Second moment of area relative to y-axis] b WNL(x) Jy = J J x2 ds dx 0 SNL( x) Jy Jy = 1.778  106

6 Центробежный момент инерции [Product moment of area] b WNLix) Jxy =J J x sdsdx 0 SNL(x) Jxy Jxy = 1.24  105

7 Абсцисса центра тяжести профиля [x-coordinate of centre of gravity] = Sy  F-1 XC XC = 43.356

8 Ордината центра тяжести [y-coordinate of centre of gravity] Yc = Sx  F-1 Yc YC = 3.902

9 Момент инерции площади относительно оси ординат центральных осей [Second moment of area relative to centroid y-axis] Jy, = Jy - X2  F Jy1 Jy = 4.263  105 y1

10 Момент инерции площади относительно оси абсцисс центральных осей [Second moment of area relative to centroid x-axis] Jx, = Jx - Yc  F JX1 Jx = 5.402  103 x1

11 Центробежный момент инерции относительно центральных осей [Product moment of area relative to centroid axes] Jx,y, = Jxy - Xc  Yc  F Jx,y, Jx,y, = 2.36  103

12 Максимальный момент инерции относительно главных центральных осей [Maximum moment of area relative to centroid axes] Jmax = (Jx, + Jy,) ' 2 + Sqr Jmax Jmax = 4.263  105

13 Минимальный момент инерции относительно главных центральных осей [Minimum moment of area relative to centroid axes] Jmin = (Jx, + Jy,) ' 2 - Sqr Jmin Jmin = 5.389  103

14 Угол наклона главных центральных осей относительно исходной системы координат [Angle of inclination of centroid axes relative to the original coordinate system] 9 = atan Jx1y1--1- _ (Jy1 - Jmin) _ ф Ф-180=0.321 П

10

0 s

-10  0

Рис. 8. Фрагмент изогнутого компрессорного профиля С4:

----главные центральные оси;  — средняя линия;--обвод профиля

[Fig. 8. Fragment of curved C4 compressor blade profile: ----centroid axes;----centreline;--profile contour]

Выводы

Разработана в MathCAD система САПР для проектирования компрессорного профиля типа С4 и решетки на его базе.

Получено аналитическое описание обводов исходного симметричного и изогнутого профилей, а также решетки на его базе.

Аналитическая форма описания обводов профиля позволяет рассчитать геометрические параметры профиля, необходимые для расчетов на прочность, а также аналитически исследовать изменение тех характеристик обводов по длине хорды профиля или обвода, которые определяют газодинамические параметры течения и свойства пограничного слоя.

СПИСОК ЛИТЕРАТУРЫ

[1] Иванов Г.С. Конструирование технических поверхностей (математическое моделирование на основе нелинейных преобразований). М.: Машиностроение, 1987. 192 с.

[2] ХолщевниковК.В. Теория и расчет авиационных лопаточных машин. М.: Машиностроение, 1970. 610 с.

[3] Селезнев К.П., Галеркин Ю.Б., Анисимов С.А., Митрофанов В.П., Подобуев Ю.С. Теория и расчет турбокомпрессоров Л.: Машиностроение (Ленинградское отделение), 1987. 392 с.

[4] Скубачевский Г.С. Авиационные газотурбинные двигатели. Конструкция и расчет деталей. 5-е изд. перераб. и доп. М.: Машиностроение, 1981. 550 с.

[5] МакаровЕ. Инженерные расчеты в MathCAD: учебный курс. СПб.: Питер, 203. 448 с.

© Виноградов Л.В., Мамаев В.К., Ощепков П.П., 2018

История статьи:

Дата поступления в редакцию: 02 февраля 2018 Дата принятия к печати: 17 апреля 2018

Для цитирования:

Виноградов Л.В., Мамаев В.К., Ощепков П.П. Компрессорный профиль С4: автоматизация проектирования // Вестник Российского университета дружбы народов. Серия: Инженерные исследования. 2018. Т. 19. № 2. С. 155—164. DOI 10.22363/2312-8143-2018-19-2155-164

Сведения об авторах:

Виноградов Леонид Валерьевич — кандидат технических наук, доцент, учебный мастер департамента машиностроения и приборостроения Инженерной академии, Российский

университет дружбы народов. Область научных интересов: системы автоматизированного проектирования, лопаточные машины, газодинамика лопаточных машин. Контактная информация: e-mail: vinogradov_lv@rudn.university

Мамаев Владимир Константинович — старший преподаватель департамента машиностроения и приборостроения Инженерной академии, Российский университет дружбы народов. Область научных интересов: системы автоматизированного проектирования, тур-бомашины, акустические и волновые процессы. Контактная информация: e-mail: mamaev_ vk@rudn.university

Ощепков Петр Платонович — кандидат технических наук, доцент департамента машиностроения и приборостроения Инженерной академии, Российский университет дружбы народов. Область научных интересов, системы автоматизированного проектирования, эксплуатация двигателей внутреннего сгорания, альтернативные топлива. Контактная информация: e-mail: oshchepkov_pp@rudn.university

COMPUTER-AIDED DESIGN OF TYPE C4: COMPRESSOR BLADE PROFILE

L.V. Vinogradov, V.K. Mamaev P.P. Oschepkov

Peoples' Friendship University of Russia (RUDN University) 6, Miklukho-Maklaya str., Moscow, 117198, Russian Federation

Abstract. The paper presents the results of designing a compressor blade profile of type C4 using CAD, which implements the method of nonlinear transformation. Girst transformation was used by the authors. All calculations and graphical representations were made in the computer algebra system Mathcad. For the calculation and graphical representation the centreline is approximated with three nonlinear segments in the form of second order Bezier curve. After determining the first and last point of the second order Bezier curve, the coordinates of the middle point on the intersection of tangential lines that start at the first and last point of the curve are determined. After the approximation of three sections with second order Bezier curves the integral function of the blade back (convex side of the contour) of the profile C4 considering the leading and trailing edges are made as arcs. According to the analytical formulas, the geometric parameters of the curved profile were calculated, which are necessary for calculating the strength of the compressor blade.

Key words: compressor, computer-aided design, profile C4, nonlinear transformation

REFERENCES

[1] Ivanov G.S. Konstruirovanie tekhnicheskikh poverkhnostei

[2] Holshevnikov K.V. Teoriya i raschet aviatsionnykh lopatochnykh mashin [Theory and analysis of aircraft impeller machines]. Moscow: Mashinostroenie Publ., 1970. 610 p.

[3] Seleznev K.P., Galyorkin J.B., Anisimov S.A., Mitrofanov VP., Podobuev U.S. Teoriya i raschet turbokompressorov [Theory and analysis of turbocompressors]. L.: Mashinostroenie Publ.

[4] Skubachvsky G.S. Aviatsionnye gazoturbinnye dvigateli. Konstruktsiya i raschet detalei [Aircraft gas turbine engines. Construction and analysis of parts]. Moscow: Mashinostroenie Publ., 1981. 550 p.

[5] Makarov E. Inzhenernye raschety v MathCAD: uchebnyi kurs [Engineering calculations in Mathcad. Educational course]. SPb.: Piter, 2003. 448 p.

Article history:

Received: February 02, 2018 Date accepted: April 17, 2018

For citation:

Vinogradov L.V., Mamaev V.K., Oschepkov P.P. (2018). Computer-aided design of type C4: compressor blade profile. RUDN Journal of Engineering Researches, 19(2), 155—164. DOI 10.22363/2312-8143-2018-19-2-155-164

Bio Note:

Leonid V. Vinogradov — Candidate of Technical Sciences, education master of the Department of Mechanical and Tool Engineering at the Academy of Engineering of Peoples' Friendship University of Russia. Research interests: computer aided design systems, turbines, gas dynamic of turbomachines. Contact information: e-mail: vinogradov_lv@rudn.university

Vladimir K. Mamaev — senior professor of the Department of Mechanical and Tool Engineering at the Academy of Engineering of Peoples' Friendship University of Russia. Research interests: computer aided design systems, turbines, acoustic and wave processes. Contact information: e-mail: mamaev_vk@rudn.university

Petr P. Oshcepkov — Candidate of Technical Sciences, Assistant Professor of the Department of Mechanical and Tool Engineering at the Academy of Engineering of Peoples' Friendship University of Russia. Research interests: computer aided design systems, internal combustion engines operation, alternative fuels. Contact information: e-mail: oshchepkov_pp@rudn.university