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

Полный текст

Введение Сложные обводы и поверхности получили широкое распространение в лопаточных машинах (компрессоры и турбины). Проектирование и расчет таких обводов и поверхностей с применением современной вычислительной техники, а также воспроизведение на станках с программным управлением требует их корректного математического описания [1]. Одним из современных методов конструирования поверхностей является метод нелинейных преобразований, который позволяет получать искомую поверхность (обвод) как образ плоскости или какой-либо известной кривой. Данные и методы 1 1 В работе применено преобразование Гирста I2 c несобственным центром F¥ [1]. У рассматриваемой специализации кремоновых инволюций In с несобственным центром F¥ есть очень важное метрическое свойство: такое преобразование является эквиформным, т.е. сохраняются площади прообраза и образа. В использованной специализации квадратичной инволюции I2 инвариантной коникой является квадратичная парабола d2 в виде кривой Bezier 2-го порядка. Прообразом компрессорного изогнутого профиля является симметричный исходный компрессорный профиль С4 [2; 3], аэродинамические характеристики которого известны из экспериментальных данных. Все вычисления и построения выполнены в интегрированном математическом пакете MathCAD [5], для чего была разработана система автоматизированного проектирования (САПР) и подпрограмма для расчета координат опорных точек полиномов Bezier 2-5-го порядков (Pp - «Опорные точки»). Объектом построения был профиль С4 [2], симметричный, с максимальной относительной толщиной 10%, хордой профиля 100 мм, относительная абсцисса сечения с максимальной толщиной - 30%, радиус входной кромки r1 = 1,2 мм, радиус выходной кромки r2 = 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 VSvxx,VSvxy - векторы абсцисс и ординат спинки профиля, x1(Θ), y1(Θ) - параметрическое уравнение окружности входной кромки. Рис. 2. Область входной кромки профиля: ■ - верхний обвод профиля; · - нижняя половина обвода профиля; ▬ - входная кромка профиля [Fig. 2. Leading edge of the blade profile: ■ - top part contour; o - bottom part contour; ▬ - leading edge] Для начала аппроксимации вариантными расчетами были определены координаты начальной точки (абсцисса и ордината) А: 2 2 XA º 0.2YA(x) = r1 -(x - r1 ) YA = YA(XA)YA = 0.663325. Начальная точка А расположена на окружности входной кромки, и аппроксимационная кривая должна быть касательной к этой окружности. Эта точка является одной (левой) опорной точкой характеристического треугольника аппроксимационной кривой Bezier. Правая опорная точка С характеристического треугольника определялась методом последовательного приближения - перебором точек матрицы координат профиля. Для первого аппроксимационного участка правой опорной точкой стала точка вектора с индексом 9. Для расчета вектора опорных точек полиномов Bezier применялась разработанная подпрограмма Pp_basic_point(XA, xC, YA, yC, αA, αC, 2). Формат матрицы исходных данных для этой подпрограммы имеет вид Pp_basic_point(XA, xC, YA, yC, αA, αC, 2), где - XA, YA, xC, yC - координаты правой и левой (конечных) опорных точек кривой Bezier; αА, αС - углы наклона касательных в этих точках. Подпрограмма формирует матрицу координат опорных точек полиномов Bezier и параметрические уравнения соответствующих полиномов. Ниже приведена матрица координат опорных точек для профиля С4. Матрица координат опорных точек полиномов Bezier ⎛ 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 ⎟ ⎜⎝ 0 0 0 40 4.89 0 0 0 0 ⎟⎠ V1 = ⎜ 6 4 2 0 0 10 20 30 Рис. 3. Начальный участок аппроксимации обвода спинки: ▬ - аппроксимационная кривая Безье; § - исходные точки контура профиля С4 [Fig. 3. Initial region of approximation for blade back contour: ▬ - Bézier approximation curve; · - points of original C4 blade profile contour] На рисунке 3 показан первый (начальный) участок аппроксимации. Второй аппроксимационный участок начинается с элемента вектора с индексом 9 (последний элемент первого аппроксимационного участка) до 14. Крайний правый элемент второго участка определялся методом последовательного подключения крайнего элемента к вектору характеристического треугольника второго участка. Координаты опорных точек кривой Bezier 2-го участка рассчитываются по подпрограмме, которая ранее уже использовалась. На рисунке 4 показан второй участок аппроксимации первичного профиля. 6 y22bz(x22) 4 Vy 2 0 0 20 40 60 80 100 x22, Vx Рис. 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 в виде SC4 (x) = 1 1 r 2 -(x - r )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 r 2 - ⎢ x -(b - r )⎥2 if 99.47 < x £ b 2 ⎣ 2 ⎦ В силу симметрии уравнение корытца можно записать в виде WC4(x) = -SC4(x). На рисунке 5 показан симметричный исходный компрессорный профиль С4, аппроксимированный двумя окружностями входной и выходной кромок, тремя кривыми Bezier 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 Bézier curves] Построение изогнутого профиля методом нелинейных преобразований было выполнено по уравнениям вида SNL(x) = LmidC4(x) - SC4(x) WNL(x) = LmidC4(x) - WC4(x) где LmidC4(x) - инвариантная кривая 2-го порядка в виде кривой Bezier с углом на входе χ1 = 14°, на выходе χ2 = -10°, т.е. угол поворота потока в решетке ϑ = 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) = x ×cos⎛ 70 × p ⎞ + S (x)×sin⎛ 70 × p ⎞ ⎝⎜ 180⎟⎠ NL ⎝⎜ 180⎟⎠ y¢(x) = - x ×sin⎛ 70 × p ⎞ + S (x)×cos⎛ 70 × p ⎞ ⎝⎜ 180⎟⎠ NL ⎝⎜ 180⎟⎠ x ¢w(x) = x ×cos⎛ 70 × p ⎞ +W (x)×sin⎛ 70 × p ⎞ ⎝⎜ 180⎟⎠ NL ⎝⎜ 180⎟⎠ y¢w(x) = - x ×sin⎛ 70 × p ⎞ +W (x)×cos⎛ 70 × p ⎞ ⎜⎝ 180⎟⎠ NL ⎜⎝ 180⎟⎠ 20 0 y ¢(x) y ¢w(x) y ¢(x) y ¢w(x) -20 -40 -60 -80 -100 0 50 100 150 x ¢(x), x ¢w(x), x ¢(x) + t, x ¢w(x) + t Рис. 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 с главными центральными осями. Геометрические параметры изогнутого профиля С4 [Table. Geometric parameters of curved C4 compressor blade profile] Таблица № п/п Наименование параметров Имя Значение 1 Площадь профиля [Profile area] b F = ò (WNL (x) - SNL (x))dx 0 F F = 718.926 2 Статический момент инерции площади относительно оси x [First moment of area relative to x-axis] b WNL( x ) Sx = ò ò s ds dx 0 SNL( x ) Sx Sx = 2.805 · 103 3 Статический момент инерции площади относительно оси y [First moment of area relative to y-axis] b WNL( x ) Sy = ò ò x ds dx 0 SNL( x ) Sy Sy = 3.117 · 104 4 Момент инерции площади относительно оси x [Second moment of area relative to x-axis] b WNL( x ) Jx = ò ò s2 ds dx 0 SNL( x ) Jx Jx = 1.635 · 104 5 Момент инерции площади относительно оси y [Second moment of area relative to y-axis] b WNL( x ) Jy = ò ò x 2 ds dx 0 SNL( x ) Jy Jy = 1.778 · 106 6 Центробежный момент инерции [Product moment of area] b WNL( x ) Jxy = ò ò x × s ds dx 0 SNL( x ) Jxy Jxy = 1.24 · 105 7 Абсцисса центра тяжести профиля [x-coordinate of centre of gravity] XC = 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 - X 2 · F 1 C Jy 1 Jy = 4.263 · 105 1 10 Момент инерции площади относительно оси абсцисс центральных осей [Second moment of area relative to centroid x-axis] Jx = Jx - Y 2 · F 1 C Jx 1 Jx = 5.402 · 103 1 11 Центробежный момент инерции относительно центральных осей [Product moment of area relative to centroid axes] Jx y = Jxy - XC · YC · F 1 1 Jx y 1 1 Jx y = 2.36 · 103 1 1 12 Максимальный момент инерции относительно главных центральных осей [Maximum moment of area relative to centroid axes] Jmax = (Jx + Jy ) · 2-1 + Sqr 1 1 Jmax Jmax = 4.263 · 105 13 Минимальный момент инерции относительно главных центральных осей [Minimum moment of area relative to centroid axes] Jmin = (Jx + Jy ) · 2-1 - Sqr 1 1 Jmin Jmin = 5.389 · 103 14 Угол наклона главных центральных осей относительно исходной системы координат [Angle of inclination of centroid axes relative to the original coordinate system] j = a tan ⎡J × 1 ⎤ ⎢ x1y1 (J - J ) ⎥ ⎣ y1 min ⎦ φ 180 j× = 0.321 p 10 0 -10 0 20 40 60 80 Рис. 8. Фрагмент изогнутого компрессорного профиля С4: - главные центральные оси; ··· - средняя линия; ─ - обвод профиля [Fig. 8. Fragment of curved C4 compressor blade profile: - centroid axes; ··· - centreline; ─ - profile contour] Выводы Разработана в MathCAD система САПР для проектирования компрессорного профиля типа С4 и решетки на его базе. Получено аналитическое описание обводов исходного симметричного и изогнутого профилей, а также решетки на его базе. Аналитическая форма описания обводов профиля позволяет рассчитать геометрические параметры профиля, необходимые для расчетов на прочность, а также аналитически исследовать изменение тех характеристик обводов по длине хорды профиля или обвода, которые определяют газодинамические параметры течения и свойства пограничного слоя.

Об авторах

Леонид Валерьевич Виноградов

Российский университет дружбы народов

Автор, ответственный за переписку.
Email: vinogradov_lv@rudn.university

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

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

Владимир Константинович Мамаев

Российский университет дружбы народов

Email: mamaev_vk@rudn.university

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

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

Петр Платонович Ощепков

Российский университет дружбы народов

Email: oshchepkov_pp@rudn.university

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

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

Список литературы