CC BY
136
УДК 629.733:004.942:629.7.035.5
КОМПЬЮТЕРНОЕ МОДЕЛИРОВАНИЕ ВОЗДУШНОГО ВИНТА ПРИ ИССЛЕДОВАНИИ
*
ЕГО ТЯГОВО-МОЩНОСТНЫХ ХАРАКТЕРИСТИК
Р. А. Нейдорф, Ю. Л. Сигида
Кафедра «Программное обеспечение вычислительной техники и автоматизированных систем», ФГБОУ ВПО «Донской государственный технический университет», г. Ростов-на-Дону; [email protected]
Ключевые слова и фразы: аэродинамические свойства; аэростат; воздушный винт; компьютерное моделирование; математическая модель; стратосфера.
Аннотация: Рассмотрены возможности имитационного моделирования воздушно-винтового движителя. Поставлена и решена задача исследования тяговых и мощностных возможностей воздушно-винтового движителя типичной для дирижаблей конфигурации. Исследование проведено путем построения математической модели, связывающей установившиеся значения мощности на валу и тягового усилия на оси винта со скоростью его вращения и фактической скоростью движения воздушного судна. Обоснован выбор ориентированной на решаемую задачу расчетной схемы компьютерного моделирования аэро- и термодинамических процессов, протекающих при вращении винта. Дано описание разработанного для задачи плана эксперимента, предусматривающего варьирование основных параметров работы винта, движения и среды полета. Приведены результаты проверки.
Введение. При математическом моделировании динамики движения летательных аппаратов возникает задача идентификации ряда констант и функций, которые не могут быть рассчитаны аналитически. К таким функциям относятся зависимости аэродинамических коэффициентов взаимодействия корпуса аппарата с воздушным потоком [1 - 3], а также тягово-мощностные рабочие характеристики активных движителей, используемых в полете воздушным судном. В настоящее время на дирижаблях и других воздушных кораблях с аэростатическим принципом удержания в воздухе применяются исключительно винтовые воздушные движители. Несмотря на внешнюю простоту конструкции винта, процессы, протекающие при его вращении и взаимодействии с частицами воздуха, чрезвычайно разнообразны и сложны [4]. Поэтому как в научной, так и в технической литературе практически отсутствуют приемлемые уравнения связи «форма лопасти -расход топлива - мощность на валу - обороты винта - тяговое усилие - скорость продольного движения винта с буксируемым корпусом относительно потока».
В связи с этим приходится решать задачи построения математической модели воздушно-винтового движителя и оптимизации его геометрических параметров с использованием натурных экспериментов, либо компьютерного моделирования аэродинамических процессов. В связи с огромной сложностью аналитического решения нелинейных систем дифференциальных уравнений в частных производных в настоящее время все больше предпочтения отдается последнему подходу.
* По материалам доклада на конференции ММТТ-2014.
Рассмотрим задачу, в которой ставились методологические цели: освоение и демонстрация технологии эксперимента по компьютерному моделированию воздушного винта и математической обработке его результатов с получением рабочих характеристик винта.
Постановка задачи. Для исследования выбран пятилопастной винт диаметром 2 м, предназначенный для установки на дирижаблях [5]. Необходимо с использованием методов и средств компьютерного моделирования получить данные по его мощностным, тяговым и скоростным свойствам. На основании этих данных требуется графически построить и исследовать рабочие характеристики воздушно-винтового движителя при движении воздушного судна вблизи земной поверхности. Исследованию подлежат наиболее информативные характеристики, связывающие момент и полную мощность на валу винта, тягу на его оси, полезную мощность, преобразованную в движение, и коэффициент полезного действия винта с частотой его вращения и фактической скоростью движения относительно воздушного потока. Максимальное значение данной скорости для дирижаблей составляет ~ 200 км/ч, то есть около 60 м/с.
Методы и инструменты решения задачи. На рисунке 1, a изображена геометрическая модель выбранного для исследования воздушного винта. Для компьютерного моделирования применен метод SMM (Sliding Mesh Model) [6], основанный на вращении расчетной сетки, который является наиболее точным, однако требует больших вычислительных мощностей компьютера. Расчетная область для исследования винта изображена на рис. 1, б.
Параметры имитационного моделирования винта:
Предполагаемая скорость движения дирижабля, м/с .. 0... 70 (шаг 10 м/с)
Давление, Па.................................................................... 101325
Скорость вращения пропеллера, об/мин....................... 1000 - 3000
(шаг 500 об/мин)
Температура окружающей среды, °С ........................... 10
Временной шаг моделирования, с ................................ 1,11 • 10-3
Угол поворота за один шаг............................................ 10°
Число шагов.....................................................................360
Модель турбулентности ................................................k-g realizable
На поверхностях фронтального, бокового и заднего интерфейсов (см. рис. 1, б) измерялись при компьютерном моделировании давление Р, скорость V и массовый расход О. Кроме того, измерялись механические и энергетические параметры взаимодействия винта с воздушным потоком: момент на валу Мх и осевая тяга Рх.
а) б)
Рис. 1. Геометрическая модель винта (а) и расчетная область (б)
Кроме данных эксперимента вычислены параметры, не поддающиеся прямому измерению или моделированию, к которым относятся следующие величины:
(1)
■shaft = Мхюх;
Ndrag FxV;
K = ^use _
N
drag
N
shaft
K = — ^loss _
Ns
shaft
N,
(2)
(3)
(4)
drag
где А^ьаД; - мощность, развиваемая двигателем на валу винта; - мощность,
затрачиваемая на движение воздушного судна с преодолением сопротивления воздуха; Кше, - коэффициенты полезного действия и потери мощности
соответственно. Коэффициент введен как виртуальная характеристика
свойств винта и соответствует случаю, когда скорость потока превышает максимальную скорость потока, которую может создать винт при заданной скорости вращения. В этом случае поток создает на валу винта противомомент и противо-силу, так что эти переменные меняют знак. Знак минус в уравнении (4) отражает потерю мощности.
Анализ полученных данных. На рисунке 2 показано семейство зависимостей тяги ¥х, развиваемой винтом, от скорости движения винта в воздушной среде при различных скоростях вращения винта.
С увеличением скорости движения тяга уменьшается. При увеличении числа оборотов винта в единицу времени тяга увеличивается. Уменьшение тяги с ростом скорости движения при постоянной скорости вращения винта приводит к тому, что при некоторой скорости vcr она становится нулевой. При дальнейшем росте скорости начинается торможение винта, и сила становится отрицательной.
Величина V,,,. является пределом скоростных возможностей винта для создания тяги и показателем его скоростной эффективности. Результаты эксперимента позволяют построить и использовать при выборе типа винта графики зависимости скоростной эффективности различных конструкций от скорости вращения. Для исследуемого винта данная зависимость оказалась практически линейной.
Важной эксплуатационной характеристикой воздушного винта является зависимость развиваемой им при неподвижном положении тяги от скорости вращения (рис. 3). Она показывает стартовые характеристики движителя и хорошо описывается близкой к параболической зависимостью вида
F-10-10" 9>
8 7 ■ 6 54 3 2 1
, Н
6 9f
8 S
4
3 и ч < 7>f
2
1
Fdrag (ю) = -31,39ю - 0,003014ю2
+
+0,000001272ю3
(5)
0 10 20 30 40 50 60 V, м/с
Рис. 2. Зависимость тяги от скорости движения воздушного судна при различных скоростях вращения винта, об/мин:
1 - 500; 2 - 1000; 3 - 1500; 4 - 2000; 5 - 2500; 6 - 3000; 7- 9 - сила лобового сопротивления при фронтальной площади 6; 11,6 и 23 м2
2 2,5 3 <в-10-3, об/мин
Рис. 3. Тяга, развиваемая в неподвижном положении:
• - экспериментальные данные; компьютерное моделирование
Семейство графиков тяговых зависимостей винта от скорости его вращения удобно использовать в сочетании с аэродинамическими характеристиками лобового сопротивления воздушного судна. Совместно они позволяют построить более информативную рабочую характеристику системы «аппарат -движитель», которая определяет маршевые возможности воздушного судна.
На рисунке 2 вместе с тяговыми характеристиками винта построены зависимости лобового сопротивления ^ (5", V) макета классического дирижабля типа «парсеваль» с масштабируемой характеристической площадью 5 максимального фронтального сечения. Абсциссы точек пересечения характеристик являются значениями равновесной скорости воздушного судна моделируемой формы с маршевым двигателем, оснащенным исследуемым винтом.
Таким образом, возникает возможность оценивать и корректировать свойства системы «дирижабль - винт» (Д-В система) [7] в следующих аспектах.
Во-первых, можно с высокой точностью оценивать и моделировать свойства построенной Д-В системы. Во-вторых - корректировать геометрические параметры винта, добиваясь требуемых тактико-технических данных проектируемой Д-В системы. В-третьих, появляется возможность уточнять геометрические и аэродинамические характеристики аппарата, для которых выбран конкретный маршевый движитель, изменять конструкцию нельзя или нецелесообразно.
Наиболее эффективным способом решения задачи построения Д-В системы является ее оптимизация. При этом наиболее эффективным критерием оптимизации эксплуатационных характеристик дирижабля является максимизация -Ки8е. В связи с этим представляет интерес исследование полученных в эксперименте данных в целях определения оптимального скоростного режима работы рассматриваемого в качестве примера винта.
На рисунке 4, а построены зависимости полученных из эксперимента значений ^Ьай и Кривые зависимости мощности, отдаваемой двигателем на валу
^Жай (жирные линии со сплошными маркерами) от скорости движения воздушного судна расположены выше кривых зависимости мощности, затрачиваемой на движение этого судна (тонкие линии с контурными маркерами). Превы-
шение мощности на валу винта над отдаваемой им полезной мощностью увеличивается с уменьшением скорости. При V = 0 мощность движения тоже равна 0.
На рисунке 4, б построены кривые зависимости Ки8е, переходящие при пересечении нуля в К[088, для обособленно работающего в движущемся потоке винта. Следует отметить особенность, что максимум Ки8е оказывается практически независимым (с учетом погрешности компьютерного эксперимента) от скорости вращения винта. Их значения для исследуемой конструкции находятся в пределах 0,655...0682, а максимумы соответствуют скоростям движения, которые при заданных оборотах винта минимизируют энергозатраты на полет дирижабля. Это открывает перспективы эффективной оптимизации двигательных систем воздушных судов с винтовыми приводами.
АЧ0
6 5,5 5 4,5 4 3,5 3 2,5 2 1,5 1
0,5 0
- 0,5 - 1
- 1,5 -2
, Вт
3- i-1
10 \
5, \
,71 \
8
9 \
7 \
1—
---- \ !
1 \ ^ \ s
\ \
kk
use
1
0,9 0,8 0,7 0,6 0,5 0,4 0,3 0,2 0,1
1 3 1 4 , 6 1 71 8
1
0 10 20 30 40 50 60 V, м/с 0 10 20 30 40 50 60 V, м/с а) б)
Рис. 4. Графики зависимостей мощностей движения (1, 3, 5, 7, 9,11)
и на валу (2, 4, 6, 8,10,12) (а) и КПД (б) от скорости движения воздушного судна при различных скоростях вращения винта, об/мин:
1, 2,13 - 500; 3, 4, 14 - 1000; 5, 6,15 - 1500; 7, 8,16 - 2000; 9,10,17 - 2500; 11, 12, 18 - 3000
Анализ построенных характеристик (см. рис. 2) показывает, что при юх е (1000,..., 3000) оптимальные по Ки8е скорости дирижабля приблизительно отвечают ряду V е (20,..., 60), следовательно можно вывод о том, какие обороты винта необходимо устанавливать при различных заданиях на скорость полета. Так, при установленном техническими требованиями значении крейсерской скорости дирижабля V = 50 м/с, что соответствует 180 км/ч, оптимальные обороты винта составляют 2500 об/мин. В этой точке тяговую характеристику винта пересекает характеристика дирижабля с профилем «парсеваль» и площадью сечения 11,6 м2, которая представляет график зависимости лобового сопротивления аппарата от скорости полета. Следовательно, Д-В система дирижабля исследованного винта образует оптимальную пару. Примечательно, что такая Д-В система близка к оптимуму и на других оборотах с соответствующими скоростями.
Выводы. Для задачи исследования воздушно-винтового движителя компьютерное моделирование позволило оценить тяговые и мощностные характеристики выбранного винта, найти оптимальные режимы его работы и показать практическую независимость оптимального значения КПД от скорости движения.
Список литературы
1. Impact of the Feeder Aerodynamics Characteristics on the Power of Control Actions in Steady and Transient Regimes : SAE Technical Paper : 2012-01-2112 / Pshikhopov V. [et al.]. - 2012.
2. Stability Analysis of the MAAT Feeder Airship During Ascent and Descent with Wind Disturbances : SAE Technical Paper : 2013-01-2111 / Neydorf R. [et al.]. -2013.
3. Aerodynamic Characteristics Study and Possible Improvements of MAAT Feeder Airships : SAE Technical Paper : 2013-01-2112 / Voloshin V. [et al.]. - 2013.
4. Нейдорф, Р. А. Специфика построения расчетной области при компьютерном моделировании взаимодействия вращающихся тел и воздушного потока / Р. А. Нейдорф, Ю. Л. Сигида // Сборник материалов Девятой Всероссийской научно-практической конференции «Перспективные системы и задачи управления» и Четвертой молодежной школы-семинара «Управление и обработка информации в технических системах» / Координац. науч.-техн. центр систем упр. Южного Фе-дер. округа при Южном федер. университете. - Таганрог, 2014. - C. 23 - 30.
5. GrabCAD [Электронный ресурс]. - Режим доступа : https://grabcad.com/ home (дата обращения: 10.09.2013).
6. Modeling Flows Using Sliding and Deforming Meshes [Электронный ресурс] / Fluent Inc. - Режим доступа : http://aerojet.engr.ucdavis.edu/fluenthelp/html/ug/ node436.htm (дата обращения: 11.08.2013).
7. Propeller for Airships [Электронный ресурс] // MT-Propeller. - Режим доступа : http://www.mt-propeller.com/en/entw/pro_air.htm (дата обращения: 8.08.2013).
Computer Simulation of Airscrew in the Study of its Traction Power Capacity
R. A. Neudorf, Yu. L. Sigida
Department "Software Computer Technology and Automated Systems", Don State Technical University, Rostov-on-Don; [email protected]
Key words and phrases: aerodynamic properties; airscrew; airship; computer simulation; mathematical model; stratosphere.
Abstract: The possibilities of simulation of airscrew propeller have been studied. The authors stated and solved the problem of studying traction and power capabilities of airscrew propeller of a typical airship configuration. The authors constructed a mathematical model relating the steady-state values of the shaft power and traction effort on the screw axis with a speed of rotation and the actual speed of the aircraft. The choice of the design scheme for computer simulation of aerodynamic and thermodynamic processes occurring during the rotation of the screw was justified. The authors described the experiment plan, permitting the variation of the basic parameters of the screw, movement and flight environment. The test results were reported.
References
1. Pshikhopov V., Medvedev M., Neydorf R., Krukhmalev V. et al. Impact of the Feeder Aerodynamics Characteristics on the Power of Control Actions in Steady and Transient Regimes, SAE Technical Paper 2012-01-2112, 2012, doi: 10.4271/2012-012112.
2. Neydorf R., Sigida Y., Voloshin V., Chen Y. Stability Analysis of the MAAT Feeder Airship During Ascent and Descent with Wind Disturbances, SAE Technical Paper 2013-01-2111, 2013, doi:10.4271/2013-01-2111.
3. Voloshin V., Chen Y., Neydorf R., Boldyreva A. Aerodynamic Characteristics Study and Possible Im-provements of MAAT Feeder Airships, SAE Technical Paper 2013-01-2112, 2013, doi:10.4271/2013-01-2112.
4. Neidorf R.A., Sigida Yu.L. Sbornik materialov Devyatoi Vserossiiskoi nauchno-prakticheskoi konferentsii "Perspektivnye sistemy i zadachi upravleniya" i Chetvertoi molodezhnoi shkoly-seminara "Upravlenie i obrabotka informatsii v tekhnicheskikh sistemakh" (Proceedings of the Ninth All-Russian scientific-practical conference "Advanced Systems and Control Problems" and the Fourth Youth School-Seminar "Management and information processing in technical systems"), Taganrog, 2014, pp. 23-30.
5. https://grabcad.com/home (дата обращения: 10.09.2013).
6. http://aerojet.engr.ucdavis.edu/fluenthelp/html/ug/node436.htm (accessed 11.08.2013).
7. http://www.mt-propeller.com/en/entw/pro_air.htm (accessed 8.08.2013).
Computermodellierung der Luftschraube bei der Forschung ihrer Zugleistungscharakteristiken
Zusammenfassung: Es werden die Möglichkeiten der Imitationmodellierung des Luftschraubefahrwerks betrachtet. Es ist die Aufgabe der Forschung der Zug- und Leistungsmöglichkeiten des Luftschraubefahrwerks der für die Luftschiffe typischen Konfiguration gestellt und entschieden. Die Forschung wird durch die Konstruktion des mathematischen Modells, das die bestimmten Bedeutungen der Macht auf der Welle und der Zugbemühungen auf der Achse der Schraube mit der Geschwindigkeit ihres Drehens und mit der tatsächlichen Geschwindigkeit der Bewegung des Luftfahrzeugs verbindet, durchgeführt. Es wird die Auswahl des auf die gelöste Aufgabe ausgerichteten Rechenschemas der Computermodellierung der luft- und thermodynamischen Prozesse rechtfertigt, die beim Drehen der Schraube verlaufen. Es wird der für die Aufgabe entwickelte Plan des Experimentes, der das Variieren der Hauptparameter der Arbeit der Schraube, der Bewegung und der Umgebung des Flugsvorsieht, beschrieben. Es werden die Ergebnisse der Prüfung gebracht.
Modélage informatique d'hélice d'avion lors de l'étude de ses caractéristiques de traction et de puissance
Résumé: Sont examinées les possibilités du modélage d'imitation du propulseur d'hélice d'avion. Est mis et résolu le problème de l'étude des possibilités de traction et de puissance du propulseur d'hélice d'avion typique pour la configuration des dirigeables. L'étude est accompagnée de la construction du modèle mathématique. Est argumenté le choix du schéma du modélage informatique des processus aéro- et thermodynamiques. Est décrit le plan de l'expériment. Sont cités les résultats du contrôle.
Авторы: Нейдорф Рудольф Анатольевич - доктор технических наук, профессор, заведующий кафедрой «Программное обеспечение вычислительной техники и автоматизированных систем»; Сигида Юрий Леонидович - аспирант кафедры «Программное обеспечение вычислительной техники и автоматизированных систем», ФГБОУ ВПО «Донской государственный технический университет», г. Ростов-на-Дону.
Рецензент: Туголуков Евгений Николаевич - доктор технических наук, профессор кафедры «Техника и технологии производства нанопродуктов», ФГБОУ ВПО «ТГТУ».
