ТУРБИННЫЙ ПРОФИЛЬ A3K7 - NACA: АВТОМАТИЗАЦИЯ ПРОЕКТИРОВАНИЯ

Полный текст

В NACA был разработан турбинный реактивный профиль (серия A3K7), для построения которого задана координатным способом форма средней линии изогнутого профиля и исходный симметричный профиль, толщина которого наращивается на среднюю линию. (Это один из возможных вариантов проектирования турбинных профилей осевых машин [5].) Толщина профиля приводится в зависимости от расстояния вдоль хорды, однако, если средняя линия имеет значительную кривизну (при большом угле поворота потока), т.е. если длина средней линии на 30% превышает длину хорды, распределение толщины дается обычно вдоль средней линии [2]. В работе применено преобразование Гирста I 2 c несобственным центром F 1  [3]. У рассматриваемой специализации кремоновых инволюций I n с несобственным центром F 1  есть очень важное метрическое свойство: такое преобразование является эквиформным, т.е. сохраняют площади прообраза и образа. В использованной в данной работе специализации квадратичной инволюции I 2 , инвариантной коникой является серия квадратичных парабол d 2 в виде кривых Bezier 2-го порядка. Прообразом турбинного изогнутого профиля является симметричный исходный турбинный профиль A3K7 [2], аэродинамические характеристики которого известны из экспериментальных данных. Все вычисления и построения выполнены в интегрированном математическом пакете Mathcad, для чего была разработана система автоматизированного проектирования (САПР), включающая ряд ранее разработанных программный модулей. Объектом конструирования был профиль серии A3K7 [1] симметричный, с максимальной относительной толщиной 20%, хордой профиля 100 мм, относительная абсцисса сечения с максимальной толщиной - 20%, радиус входной кромки r 1 = 4,407% от хорды, радиус выходной кромки r 2 = 1% от хорды. Важно отметить, что NACA определяет входной и выходной углы (тангенсы углов) установки лопатки не на входной и выходной кромках, а на расстояниях 0,5 и 95% от длины хорды [2]. (Как указывалось ранее, заданы координаты средней линии в процентах от хорды.) В таблице 1 приведена матрица координат средней линии профиля, где первая, третья, пятая строки - абсциссы, вторая, четвертая, шестая - ординаты. (Матрица координат параллельно записывается в текстовый структурированный файл для возможности последующего ее использования в других средах.) Таблица 1 Матрица координат средней линии турбинного реактивного профиля A3K7 На рисунке 1 показан фрагмент векторного графика средней линии в исходной системе координат. График строится для получения первичного образа объекта проектирования и контроля корректности исходных данных и их ввода. Рис. 1. Векторный график (фрагмент) средней линии профиля турбинного A3K7 (NACA): _ - интерполяция средней линии кубическим сплайном; • • • - точки средней линии, координаты которых заданы (векторный график средней линии) Для последующих вычислений и построений средняя линия, заданная координатным способом, аппроксимируется двумя отрезками прямых в области входной и выходной кромок с заданными углами наклона и тремя криволинейными участками в виде дуг кривых Безье 2-го порядка. Координаты точек сопряжения соседних аппроксимируемых участков определяются итерационным методом с оптимизацией индекса корреляции (максимальное значение) [4]. Для вычисления индекса корреляции в Mathcad имеется специальная функция, что существенно сокращает время вычислений. Следует отметить, что полученное уравнение аппроксимированной средней линии имеет аналитический формат, что дает возможность применять цифровую технологию не зависимо от языка программирования и соответствует современным технологиям при проектировании инженерных объектов. Интегральное уравнение средней линии имеет вид (формат Mathcad) На рисунке 2 показана исходная средняя линия в виде интерполяционного графика и средняя линия, построенная по интегральному уравнению (1). Рис. 2. Средняя линия профиля A3K7: • • • - исходные точки, заданные координатным способом; _ - аппроксимированная средняя линии, построенная по формуле (1) Рассчитанный индекс корреляции между исходным вектором значений и аппроксимированными равен corr = 0,999723. С практической точки точность аппроксимации близка к точности округления исходных данных. (В сложных случаях для лучшего приближения исходных и аппроксимированных данных можно увеличить количество участков аппроксимации.) Следующим этапом проектирования было построение исходного симметричного профиля, заданного распределением толщины профиля по длине хорды. В таблице 2 представлена матрица координат обводов симметричного исходного профиля A3K7. Для оценки корректности исходных данных был построен векторный график исходного профиля и интерполяционный вариант (рис. 3). Рис. 3. Векторный график исходного турбинного профиля A3K7: • • • - исходные точки обводов профиля; _ - интерполяция исходных данных кубическим сплайном Как отмечалось ранее входная и выходная кромки профиля оформлены в виде дуг окружностей соответствующих радиусов. Как можно видеть из графика (см. рис. 3), в области выходной кромки имеется особенность обвода профиля - перегиб контура. Эту особенность следует учесть при аппроксимации профиля кривыми Безье 2-го порядка. Для определения координат точки перегиба использовалось условие равенства нулю второй производной интерполяционного обвода профиля. На рисунке 4 в формате Mathcad показан фрагмент алгоритма решения для точки перегиба контура. Определение координат точки перегиба обвода спинки профиля Рис. 4. Фрагмент программы для определения абсциссы точки перегиба контура (формат Mathcad) Следующим этапом проектирования является определение координат точки сопряжения входной кромки в виде дуги окружности и серии точек обвода контура, аппроксимированной кривой Безье 2-го порядка. Итерационным способом, варьируя абсциссу начальной точки на дуге окружности и задавая конечную точку вектора аппроксимации, при условии максимума индекса корреляции corr max были найдены координаты начальной опорной точки A кривой Безье для профиля A3K7. Продолжение этого этапа - процедура максимального увеличения интервала аппроксимации для сокращения участков аппроксимации обводов при условии максимума индекса корреляции. Для первого участка аппроксимации абсцисса начальной опорной точки равна X AC = 0,3 мм. Конечной точкой первого участка аппроксимации была точка с абсциссой равной 20 мм. Определив начальную и конечную опорные точки параболы Безье 2-го порядка, определялись координаты промежуточной опорной точки на пересечении касательных, исходящих из начальной и конечной опорных точек. (Следует напомнить, что координаты опорных точек: начальной, средней и конечной, являются коэффициентами параметрических уравнений кривой Безье.) Расчеты для первого аппроксимационного участка показали, что индекс корреляции равен corr = 0,999931. Для решения второго участка аппроксимации характеристика начальной опорной точки известна, так как эта точка есть конечная опорная точка предыдущего участка. Поэтому следует только тем же методом найти конечную опорную точку второго (последующего) участка. Для второго участка были получены следующие данные: абсцисса начальной опорной точки 20 мм, конечной опорной точки - 60 мм, индекс корреляции - corr = 0,99969. Третий участок аппроксимации. Особенность аппроксимации третьего участка заключается в том, что необходимо определить характеристики конечной опорной точки. Для заданного профиля конечной опорной точкой были координаты экстремальной точки дуги окружности выходной кромки с абсциссой равной Х с = (b - r 2 ) и ординатой r 2 . Дополнительным условие для конечной точки было условие параллельности оси абсцисс исходной системы координат касательной, исходящей из конечной точки. При этих условиях можно построить характеристический треугольник для аппроксимирующей кривой Безье 2-го порядка. Для третьего участка аппроксимации индекс корреляции был получен равным corr = 0,99981. После аппроксимации трех участков кривыми Безье 2-го порядка было записано интегральное уравнение обвода спинки (выпуклой части контура) профиля A3K7 с учетом того, что входная и выходная кромки выполнены в виде дуг окружностей. Спинка симметричного исходного профиля в аналитическом виде в формате Mathcad описывается интегральным уравнением В силу симметричности обводов профиля, уравнение вогнутой части профиля (корытца) может быть записано тем же уравнением что и спинка, но со знаком минус. На рисунке 5 показан симметричный исходной профиль A3K7 в исходной системе координат, построенный по аппроксимационному интегральному уравнению (2). Рис. 5. Симметричный исходный профиль A3K7, построенный по аппроксимационному аналитическому уравнению Как указано ранее, изогнутый турбинный реактивный профиль на базе профиля A3K7, проектируется методом нелинейных преобразований, где прообразом служит симметричный профиль, инвариантной кривой средняя линия, образом - изогнутый профиль. Для построения изогнутого профиля были получена система уравнений преобразования: На рисунке 6 показан изогнутый турбинный реактивный профиль серии A3K7 NACA с углом поворота потока 50°, построенный на базе симметричного исходного профиля и средней линии методом нелинейных преобразований Гирста с несобственным центром. Полученная аналитическая форма описания обводов симметричного и изогнутого профилей серии A3K7 позволяет аналитически рассчитать все геометрические характеристики, необходимые при расчете на прочность лопатки, формируемых профилями каналов решетки и др. Как указывалось ранее, примененный нелинейный метод обладает свойством эквиформности, т.е. площади исходного симметричного профиля и изогнутого профиля равны. Для сконструированного профиля площади оказались равными F = 1,199 · 10 3 мм 2 . Проведенная работа позволяет сделать следующие выводы. 1. Разработана система автоматизированного проектирования, которая позволяет проектировать турбинный реактивный профиль серии A3K7. 2. Предложенная форма описания обводов профиля дает возможность аналитически исследовать весь комплекс геометрических параметров профиля, строить решетки в других программах и на других языках. 4. Разработанная САПР позволяет автоматически варьировать параметры профиля и решетки: геометрический масштаб, угол изгиба профиля, угол установки профиля в решетке, относительный шаг решетки и др. 5. САПР выполнена по модульной схеме, что делает систему модернизационно пригодной: совершенствование системы, дополнение ее новыми модулями встраивание системы в другие САПР более высокого уровня, например, в САПР «Турбина» и др.

Об авторах

Л В Виноградов

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

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

Виноградов Леонид Валерьевич, кандидат технических наук, доцент, ведущий инженер департамента машиностроения и приборостроения инженерной академии РУДН. Сфера научных интересов: системы автоматизированного проектирования, турбины, агрегаты турбонаддува.

ул. Миклухо-Маклая, 6, Москва, Россия, 117198

В К Мамаев

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

Email: mamaev_vk@rudn.university

Мамаев Владимир Константинович, старший преподаватель департамента машиностроения и приборостроения инженерной академии РУДН. Сфера научных интересов: системы автоматизированного проектирования, турбины, акустические и волновые процессы.

ул. Миклухо-Маклая, 6, Москва, Россия, 117198

П П Ощепков

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

Email: oshchepkov_pp@rudn.university

Ощепков Петр Платонович, кандидат технических наук, доцент департамента машиностроения и приборостроения инженерной академии РУДН. Сфера научных интересов: системы автоматизированного проектирования, эксплуатация двигателей внутреннего сгорания, альтернативные топлива.

ул. Миклухо-Маклая, 6, Москва, Россия, 117198

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