УДК 533.6.011
DOI: 10.18698/2308-6033-2019-1-1838
Влияние стреловидности стабилизаторов на аэродинамические характеристики беспилотного летательного аппарата
© В.О. Москаленко, А.А. Косырев МГТУ им. Н.Э. Баумана, Москва, 105005, Россия
Беспилотным летательным аппаратам уделяется сегодня большое внимание в армиях многих стран мира. Установлено, что улучшить аэродинамические характеристики беспилотного летательного аппарата можно за счет изменения стреловидности стабилизаторов. Для исследования влияния угла стреловидности стабилизаторов на изменение аэродинамических характеристик беспилотного летательного аппарата Predator проведено математическое моделирование обтекания модели беспилотного летательного аппарата с различной стреловидностью стабилизаторов при дозвуковой скорости набегающего потока. В результате численного расчета в пакете SolidWorks Flow Simulation получены аэродинамические характеристики беспилотного летательного аппарата с различной стреловидностью стабилизаторов. Построены графики, отображающие зависимость аэродинамических коэффициентов от угла атаки, получены поля давления в вертикальных плоскостях течения. Проведено исследование полей давления. Проанализировано влияние картины течения на изменение аэродинамических характеристик беспилотного летательного аппарата со стабилизаторами прямой и обратной стреловидности. Выявлено преимущество беспилотного летательного аппарата со стабилизаторами обратной стреловидности перед аппаратами со стабилизаторами прямой стреловидности и без стреловидности.
Ключевые слова: беспилотный летательный аппарат, аэродинамические характеристики, программный комплекс SolidWorks, стабилизаторы обратной стреловидности, стабилизаторы прямой стреловидности
Введение. Одним из актуальных и перспективных путей совершенствования боевого потенциала современных вооруженных сил является создание и внедрение в процессы боевой подготовки и возможных боевых действий беспилотных летательных аппаратов (БПЛА) [1]. В армиях многих стран мира уделяется достаточно большое внимание боевым комплексам с БПЛА, ведется активная разработка принципиально новых моделей и конструкций планера-беспилотника, который был бы рассчитан на повышенные нагрузки в сравнении с пилотируемой авиацией. Так, успешное применение БПЛА Predator активизировало работы в США по созданию БПЛА нового поколения, которые будут способны не только поражать наземные цели, но и вести воздушный бой [2-6].
Улучшение тактико-технических характеристик БПЛА, в том числе и аэродинамических характеристик, является важной и актуальной задачей. Повлиять на аэродинамические характеристики БПЛА можно посредством изменения угла стреловидности х крыла (стабилизато-
ров). Различают крылья (стабилизаторы) прямой и обратной стреловидности соответственно с положительным и отрицательным углом стреловидности х [7]. Модели крыла с различной стреловидностью показаны на рис. 1.
Рис. 1. Модели крыла с различной стреловидностью (вид сверху)
К применению стабилизаторов (крыльев) с обратной стреловидностью (ОС) на БПЛА, особенности обтекания которых еще мало изучены, проявляется значительный интерес [8-10]. Примером летательного аппарата с крылом ОС может послужить экспериментальный самолет С-37 («Беркут»), созданный для реализации новейших перспективных решений, которые могут быть использованы на боевых самолетах будущего [11-14]. Работы по улучшению характеристик крыльев летательных аппаратов не прекращаются по настоящее время [15, 16].
Эта статья посвящена исследованию коэффициентов подъемной силы суа, лобового сопротивления сха, моментов тангажа mz и рыскания ту, а также аэродинамического качества К БПЛА Predator со стабилизаторами различной стреловидности для углов атаки а = = 0.. .14°. Проведен анализ обтекания БПЛА Predator со стабилизаторами прямой и обратной стреловидности (соответственно ПС и ОС).
Таким образом, одна из основных целей исследования — изучение влияния стреловидности стабилизаторов, включая и обратную, на аэродинамические характеристики БПЛА Predator с помощью программного комплекса SolidWorks Flow Simulation [17].
Решение тестовой задачи. Перед выполнением численного эксперимента по изучению влияния угла стреловидности стабилизаторов на аэродинамические характеристики БПЛА была решена тестовая задача по расчету аэродинамических характеристик самолета МИГ-29 [18]. На примере этой задачи проведена оценка размеров расчетной области, выбрана расчетная сетка и определены погрешности расчета аэродинамических характеристик по сравнению с известными аэродинамическими характеристиками МИГ-29 [19].
С использованием пакета SolidWorks Flow Simulation и модели турбулентности к—8 были определены аэродинамические коэффициенты: сопротивления сха и подъемной силы суа при различных углах атаки для М« = 0,6. При этом проводилась табличная адаптация сетки (в качестве единицы измерения были выбраны итерации). Таким образом ячейки сетки дробились в областях с особенностями течения для обеспечения установленной точности решения (рис. 2). В итоге количество ячеек сетки после адаптации составило порядка 1,8 млн.
Рис. 2. Расчетная сетка после неоднократного уплотнения
По результатам эксперимента [19] и численного расчета построены зависимости коэффициентов подъемной силы суа (рис. 3, а) и силы лобового сопротивления сха (рис. 3, б) от угла атаки а для МИГ-29.
На рисунке видно, что результаты численного расчета и эксперимента совпадают с ошибкой по сха примерно 15 %, по суа примерно 6 % (окрестность нуля в расчет не бралась).
Были предприняты попытки уменьшения осредненных ошибок по сха и суа. Увеличивались размеры расчетной области, а также количество ячеек с 1,8 до 2,5 млн. При этом осредненные ошибки оставались примерно в тех же пределах, в связи с чем расчеты в дальнейшем были проведены для количества ячеек, равного 1,8 млн.
Полученные систематические ошибки в определении суммарных коэффициентов сха и суа не повлияли на результаты анализа и выводы, так как анализировались не абсолютные значения коэффициентов, а их приращения А по отношению к стабилизаторам ОС и стабилизаторам ПС.
Суа 1,2 1,0 0,8 0,6 0,4 0,2
0 2 4 6 8 10 12 14 16 18 а, град
а
Сха
0,5 0,4 0,3 0,2 ОД
О 2 4 6 8 10 12 14 16 18 а, град
б
Рис. 3. Зависимость коэффициентов подъемной силы суа (а) и силы лобового сопротивления сха (б) от угла атаки а: 1, 4 — эксперимент [19]; 2, 3 — расчет
Расчет аэродинамических характеристик БПЛА со стабилизаторами различной стреловидности. После решения тестовой задачи был проведен численный эксперимент по оценке влияния угла стреловидности х стабилизаторов и угла атаки а на аэродинамические характеристики БПЛА. Согласно данным [7], оптимальным углом стреловидности, при котором достигается максимальный прирост аэродинамического качества крыла, является х = -30°. Был рассмотрен БПЛА Predator длиной L = 7,5 м и размахом крыла l = 14,84 м.
Модель БПЛА Predator показана на рис. 4. Стабилизаторы 1 и 3 имеют удлинение X = 3,4, а стабилизатор 2 — X = 1,6. При численном эксперименте изменялась стреловидность стабилизаторов 1 и 3. Стреловидность стабилизатора 2 не изменялась ввиду его малого удлинения. Рассматривались стабилизаторы как ПС при х = 30°, так и при х = -30°.
1 2
Рис. 4. Модель БПЛА Predator: 1-3 — порядковые номера стабилизаторов
Область расчета представляла собой параллелепипед, в центре которого был расположен БПЛА. Габариты расчетной области выбирались из соображений обеспечения обтекания БПЛА невозмущенным потоком, включая максимальные углы атаки, и составляли в высоту, длину и ширину соответственно 158 м х 150 м х 150 м.
Результаты численного расчета аэродинамических характеристик БПЛА с различной стреловидностью стабилизаторов. Были рассчитаны аэродинамические коэффициенты суа, сха, mz, ту, а также аэродинамическое качество К для БПЛА со стабилизаторами ОС при х = -30°, без стреловидности (БС) при х = 0 и ПС при х = 30° и Мда = 0,3 (рис. 5-9), а также зависимости приращения коэффициентов подъемной силы Асуа и аэродинамического качества ДК БПЛА со стабилизаторами ОС над стабилизаторами ПС и БС от угла атаки, которые представлены на рис. 10. Приращения коэффициентов вычислены по следующим формулам:
Асуа (ОС-ПС) = суа (ОС) — суа (ПС); (1)
АК( ОС-ПС ) = К( ОС ) — К( ПС ); (2)
Асуа (ОС—БС) = суа (ОС) — суа (БС); (3)
АК(ОС—БС) = К( ОС) — К( БС), (4)
где суа (ос) — коэффициент подъемной силы БПЛА со стабилизаторами ОС; суа (пс) — коэффициент подъемной силы БПЛА со стабилизаторами ПС; суа (бс) — коэффициент подъемной силы БПЛА со стабилизаторами БС; К(ос), К(ПС), К(БС) — аэродинамическое качество БПЛА со стабилизаторами ОС, ПС, БС соответственно.
Суа
0,7 0,6 0,5 0,4 0,3 0,2 ОД
0 2 4 6 8 10 12 а, град
Рис. 5. Зависимость коэффициента подъемной силы суа БПЛА от угла атаки а при различных углах стреловидности стабилизаторов:
1 — х = -30°; 2 — х = 0; 3 — х = 30°
Сха 0,16 0,14 0,12 0,10 0,08 0,06 0,04 0,02
0 2 4 6 8 10 12 а, град
Рис. 6. Зависимость коэффициента силы лобового сопротивления сха БПЛА от угла атаки а при различных углах стреловидности стабилизаторов:
1 — х = -30°; 2 — х = 0; 3 — х = 30°
а, град
Рис. 7. Зависимость аэродинамического качества К БПЛА от угла атаки а при различных углах стреловидности стабилизаторов:
1 — х = -30°; 2 — х = 0; 3 — х = 30°
а, град
Рис. 8. Зависимость коэффициентов момента тангажа ш2 БПЛА от угла атаки а при различных углах стреловидности стабилизаторов:
1 — х = -30°; 2 — х = 0; 3 — х = 30°
Р, град
Рис. 9. Зависимость коэффициентов момента рыскания шу БПЛА от угла скольжения р при различных углах стреловидности стабилизаторов:
1 — х = -30°; 2 — х = 0; 3 — х = 30°
ЬСуа
0,025 0,020
0,015'
0,010 0,005
0
-0,005 -0,010
О 2 4 6 8 10 12 а, град
а
Д К 0,8 0,6 0,4 0,2
0 -0,2
О 2 4 6 8 10 12 а, град
б
Рис. 10. Зависимость Асуа (а) и АК (б) БПЛА со стабилизаторами ОС, ПС
и БС от угла атаки а:
1 -АСуа (ОС-ПС); 2-АСуа (ОС-БС> 3-АК(ОС-ПС); 4-АК(ОС-БС)
На рис. 10 наблюдается положительный прирост аэродинамического качества БПЛА АК со стабилизаторами ОС по сравнению со стабилизаторами ПС во всем исследуемом диапазоне углов атаки (а = = 0.. .14°) и стабилизаторов ОС по сравнению со стабилизаторами БС при а = 0.11°.
Положительный прирост коэффициента подъемной силы БПЛА Асуа со стабилизаторами ОС по сравнению со стабилизаторами ПС и стабилизаторов ОС по сравнению со стабилизаторами БС наблюдается во всем исследуемом диапазоне углов атаки и имеет два максимума на углах атаки порядка 4 и 12°.
Наибольший прирост аэродинамического качества стабилизаторов ОС по сравнению со стабилизаторами БС составляет 0,75 при а = 4°, а стабилизаторов ОС по сравнению со стабилизаторами ПС — 0,91 при а = 6°. Этот прирост обеспечивается за счет некоторого увеличения коэффициентов суа и снижения сха для стабилизаторов ОС (см. рис. 5, 6). При этом зависимости коэффициентов продольного момента mz(a) практически не отличаются для стабилизаторов ОС, БС и ПС (см. рис. 8), что свидетельствует об имеющемся запасе про-
дольной устойчивости для всех рассмотренных стабилизаторов, включая ОС. Зависимости для коэффициентов момента рыскания также имеют отрицательный градиент, что свидетельствует о положительном запасе путевой устойчивости. Для стабилизаторов ОС эти зависимости имеют несколько большие значения, что должно благоприятно сказаться на маневренности.
В результате численного расчета были получены картины полей давления и линий тока (картин обтекания) БПЛА для различных значений углов стреловидности стабилизаторов и углов атаки. Так, на рис. 11 представлены характерные картины обтекания БПЛА со
б
Рис. 11. Картины обтекания БПЛА со стабилизаторами ПС (х = 30°) и ОС (х = -30°) под углом а = 12°
стабилизаторами ПС и ОС при а = 12°. Видна разница между обтеканием стабилизаторов БПЛА ПС и ОС. Так, при обтекании стабилизатора ПС набегающий поток перемещается от бортовой хорды к концевой (рис. 11, а), а при обтекании стабилизатора ОС поток перемещается в обратном направлении (рис. 11, б). Из-за смены направления перетекания у стабилизатора ОС снижается интенсивность концевого вихря, что уменьшает индуктивное сопротивление и повышает общее аэродинамическое качество крыла.
Для БПЛА со стабилизаторами различной стреловидности формируются разные картины полей давления (рис. 12). Зоны изменения давления находятся внутри расчетной области, что подтверждает правильность выбора ее границ. Вблизи передней кромки стабилизатора
а
б
Рис. 12. Поле давления на расстоянии г = 1 м от сечения ХОУ при обтекании дозвуковым потоком БПЛА под углом атаки а = 12° (в сечении: крыло 4 и стабилизатор 5) при х = 30° (а) и х = -30° (б)
давление по сравнению с давлением набегающего потока увеличивается (области 1 красного цвета), на верхней части профиля создается область пониженного давления (области 2 голубого и области 3 зеленого цвета).
Заключение. Таким образом, в результате анализа приведенных численных расчетов установлено, что применение стабилизаторов ОС с углом стреловидности х = -30° по сравнению со стабилизаторами БС позволяет увеличить коэффициент подъемной силы БПЛА во всем исследуемом диапазоне углов атаки (а = 0...14°), а также аэродинамическое качество при а = 0.11°. Наибольший прирост аэродинамического качества БПЛА составляет 0,75 при а = 4°, а наибольший прирост коэффициента подъемной силы — 0,012 при а = 6°.
Применение стабилизаторов ОС с углом стреловидности х = -30° по сравнению со стабилизаторами ПС с х = 30° позволяет увеличить коэффициент подъемной силы БПЛА и аэродинамическое качество во всем исследуемом диапазоне углов атаки. Наибольший прирост аэродинамического качества БПЛА составляет 0,91 при а = 6°, а наибольший прирост коэффициента подъемной силы — 0,027 при а = 4°. При этом в исследованном диапазоне углов атаки а и скольжения в для всех типов стабилизаторов БПЛА имеет запасы продольной и путевой устойчивости.
Необходимо отметить, что БПЛА со стабилизаторами ОС в исследованном диапазоне углов атаки имеет преимущество в аэродинамическом качестве перед БПЛА со стабилизаторами ПС и БС.
ЛИТЕРАТУРА
[1] Полтавский А.В., Жумабаева А.С., Бикеев Р.Р. Многофункциональные комплексы беспилотных летательных аппаратов: развитие в системе вооружения. Надежность и качество сложных систем, 2016, № 1 (13), с. 39-46.
[2] Балыко Ю. Беспилотники НАТО в локальных конфликтах. Военный парад, 2008, № 1-2, с. 38-39.
[3] Мужичек С.М. Беспилотные летательные аппараты: мода или насущная необходимость. Фазотрон, 2005, № 3, с. 52-55.
[4] Семенов С.С. Американская стратегия развития боевых беспилотных самолетов. Аэрокосмическое обозрение, 2008, № 3, с. 21-23.
[5] Романов И.В. Беспилотные летательные аппараты в XXI веке. Наука среди нас, 2018, № 4 (8), с. 1-2.
[6] Полтавский А.В., Бородуля В.М. Развитие беспилотной авиации в армиях зарубежных стран. Стратегическая стабильность, 2007, № 1, с. 45-53.
[7] Москаленко В.О., Косырев А.А. Влияние угла стреловидности на аэродинамические характеристики крыла. Инженерный журнал: наука и инновации, 2018, вып. 2, с. 1-2. DOI: 10.18698/2308-6033-2018-2-1735
[8] Chambers J.R. Modeling Flight. Washington, US National Aeronautics and Space, 2010, 202 p.
[9] Frederick A.J. Sweeping Forward. The National Aeronautics and Space Administration, 2013, 328 p.
[10] Hepperle M. MDO of Forward Swept Wings. Braunschweig, DLR Institute of Aerodynamics and Flow Technology, 2008, pp. 14-19.
[11] Бауэрс П. Летательные аппараты нетрадиционных схем. Москва, Мир, 1991, 320 с.
[12] Ильин В. А. Экспериментальный самолет С-37 «Беркут». Вестник авиации и космонавтики, 1999, № 3, с. 31-33.
[13] Ильин В.А. Полет «Беркута». Авиация и космонавтика, 1998, № 1, с. 1-3.
[14] Докучаев А. «Беркут» раздвигает горизонты. Красная звезда, 25 сентября 2000.
[15] Коробкова Ю.П., Москаленко В.О. Численный расчет взаимодействия плоской струи с сносящим дозвуковым потоком. Аэрокосмический научный журнал, 2017, № 3, с. 53-63.
[16] Москаленко В.О., Красников И.Ю. Численный расчет взаимодействия плоской струи с сносящим дозвуковым потоком. Аэрокосмический научный журнал, 2016, № 2, с. 30-40.
[17] Алямовский А.А. SolidWorks Simulation. Как решать практические задачи. Санкт-Петербург, БХВ-Петербург, 2012, 445 с.
[18] Гордон Е., Фомин А., Михеев А. МИГ-29. Любимая книга, 1998, 256 с.
[19] Практическая аэродинамика самолета МИГ-29. Министерство обороны СССР, 1987, с. 14-35.
Статья поступила в редакцию 09.10.2018
Ссылку на эту статью просим оформлять следующим образом: Москаленко В.О., Косырев А.А. Влияние стреловидности стабилизаторов на аэродинамические характеристики беспилотного летательного аппарата. Инженерный журнал: наука и инновации, 2019, вып. 1.
http://dx.doi.org/10.18698/2308-6033-2019-1-1838
Москаленко Валерий Осипович — канд. техн. наук, доцент кафедры «Динамика и управление полетом ракет и космических аппаратов» МГТУ им. Н.Э. Баумана. Автор более 150 научных работ в области аэрогазодинамики и управления движением летательных аппаратов. e-mail: [email protected]
Косырев Антон Александрович — студент 6-го курса кафедры «Динамика и управление полетом ракет и космических аппаратов» МГТУ им. Н.Э. Баумана. e-mail: [email protected]
The effect of stabilizer sweep on aerodynamic characteristics of UVA
© V.O. Moskalenko, A.A. Kosyrev Bauman Moscow State Technical University, Moscow, 105005, Russia
Nowadays armed forces of many countries of the world pay a lot of attention to unmanned aerial vehicles (UAV). It's stated that estimation of stabilizers sweep can improve aerodynamic characteristics of UAV. We conducted mathematical simulation of UAV streamlining with different stabilizer sweep in conditions of approaching windflow at subsonic speed. It helped to study the impact of stabilizer sweep angle on changes of aerodynamic characteristics of UAV Predator. As a result of numeral computations with SolidWorks Flow Simulation, we received aerodynamic characteristics of UAV with different stabilizer sweep. We plotted graphs showing the dependence of aerodynamic coefficient from attack angle. Pressure fields in vertical plane of the stream were received. Research of pressure fields was conducted. We analysed the impact of the stream on changes of aerodynamic characteristics of UAV with back sweep and forward sweep stabilizers. We detected advantages of UAV with forward sweep stabilizers over UAV with back sweep stabilizers or without any sweep.
Keywords: unmanned aerial vehicles, aerodynamic characteristics, program complex SolidWorks, forward sweep stabilizers, back sweep stabilizers
REFERENCES
[1] Poltavskiy A.V., Zhumabaeva A.S., Bikeev R.R. Nadezhnost i kachestvo slozhnykh sistem — Reliability and Quality of Complex Systems, 2016, no. 1 (13), pp. 39-46.
[2] Balyko Y. Voenny parad — Military parade, 2008, no. 1-2, pp. 38-39.
[3] Muzhichek S.M. Fazotron — Phazotron, 2005, no. 3, pp. 52-55.
[4] Semenov S.S. Aehrokosmicheskoe obozrenie — Aerospace review, 2008, no. 3, pp. 21-23.
[5] Romanov I.V. Nauka sredi nas (Science is among us), 2018, no. 4 (8), pp. 1-2.
[6] Poltavskiy A.V., Borodulya V.M. Strategicheskaya stabilnost (Strategic stability), 2007, no. 1, pp. 45-53.
[7] Moskalenko V.O, Kosyrev A.A. Inzhenernyy zhurnal: nauka i innovatsii — Engineering Journal: Science and Innovation, 2018, issue 2, pp. 1-2. DOI: 10.18698/2308-6033-2018-2-1735
[8] Chambers J.R. Modeling Flight. Washington, US National Aeronautics and Space, 2010, 202 p.
[9] Frederick A.J. Sweeping Forward. The National Aeronautics and Space Administration, 2013, 328 p.
[10] Hepperle M. MDO of Forward Swept Wings. Braunschweig, DLR Institute of Aerodynamics and Flow Technology, 2008, pp. 14-19.
[11] Bowers P.M. Boeing Aircraft Since 1916. Naval Institute Press, 1989, 668 p. [In Russ.: Bowers P. Letatelnye apparaty netraditsionnykh skhem. Moscow, Mir Publ., 1991, 320 p.].
[12] Il'in V.A. Vestnik aviatsii i kosmonavtiki — Aerospace herald, 1999, no. 3, pp. 31-33.
[13] Il'in V.A. Aviatsiya i kosmonavtika (Aviation and Astronautics), 1998, no. 1, pp. 1-3.
[14] Dokuchaev A. Krasnaia zvezda (Red Star), 25 sentyabrya 2000.
[15] Korobkova Yu.P., Moskalenko V.O. Aerokosmicheskiy nauchnyy zhurnal — Aerospace Scientific Journal, 2017, no. 3, pp. 53-63.
[16] Moskalenko V.O., Krasnikov I.Y. Aerokosmicheskiy nauchnyy zhurnal — Aerospace Scientific Journal, 2016, no. 2, pp. 30-40.
[17] AlyamovskiyA.A. SolidWorks Simulation. Kak reshat prakticheskie zadachi [SolidWorks Simulation. How to solve practical problems]. St. Petersburg, BKhV-Peterburg Publ., 2012, 445 p.
[18] Gordon E., Fomin A., Mikheev A. MIG-29 [MiG-29]. Lyubimaya kniga Publ., 1998, 256 p.
[19] Prakticheskaya aehrodinamika samolyota MIG-29 [MiG-29 practical aerodynamics]. Ministerstvo oborony SSSR Publ., 1987, pp.14-35.
Moskalenko V.O., Cand. Sc. (Eng.), Assoc. Professor, Department of Dynamics and Flight Control of Rockets and Spacecraft, Bauman Moscow State Technical University. Author of over 150 research papers in the field of aerogasdynamics and aircraft motion control. e-mail: [email protected]
Kosyrev A.A., 6th year student, Department of Dynamics and Flight Control of Rockets and Spacecraft, Bauman Moscow State Technical University. e-mail: [email protected]