Ученые из МФТИ рассчитали аэродинамику машущего крыла робоптицы

Исследователи из МФТИ с коллегами рассчитали аэродинамику орнитоптера — робоптицы, которая держится в воздухе за счет машущих крыльев. Моделирование показало, что при частоте взмахов 8 Гц даже плоское крыло нулевой толщины создает не сопротивление, а тягу, причем средняя подъемная сила остается такой же, как у неподвижной модели.

Численное моделирование аэродинамического обтекания орнитоптера с использованием деформируемых расчетных сеток выполнено в рамках государственного задания Минобрнауки России №075-00558-26-00 от 12.01.2026.

Подавляющее большинство современных беспилотников построено по самолетной или вертолетной схеме. У аппаратов, построенных по этим схемам, одинаковые недостатки: сильный шум, невысокая энергоэффективность и ограниченная маневренность.

Птицы всех этих недостатков лишены. Этим объясняется устойчивый интерес инженеров к орнитоптерам — летательным аппаратам, которые используют для своей работы механику движения крыльев птиц.

Крыло птицы одновременно опускается и поднимается, поворачивается вокруг собственной оси, складывается и меняет форму. Подобная динамика выходит далеко за область применимости традиционных аэродинамических методов расчета.

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

Исследователи из МФТИ и ЦНИИ РТК пошли другим путем и применили деформируемые расчетные сетки. Сетка не перестраивается заново, а растягивается и сжимается вместе с крылом, сохраняя постоянное число ячеек. В их подходе физические законы сохранения выполняются точно, а не приближенно.

Прототипом для математической модели послужила механическая модель орнитоптера, разработанная в МФТИ. Крылья в ней сделаны из непроницаемой ткани, натянутой на направляющие спицы, размах составляет 0,908 м.

Такая конструкция позволила заменить крылья в расчете поверхностями нулевой толщины, а корпус аппарата не рассматривать вовсе: в модели остались только плоскости крыльев и хвостового оперения.

Нулевая толщина заодно исключила эффект Кноллера—Бетса, при котором тягу создает сам колеблющийся профиль конечной толщины. Значит, любая тяга, полученная в расчете, обязана исключительно кинематике взмаха.

Рисунок 1. Механическая модель орнитоптера, разработанная в МФТИ — (а), и ее математическое представление (б). Крылья размахом 0,908 м заданы поверхностями нулевой толщины, корпус аппарата в расчете не учитывается. Источник: Вестник ЮУрГУ. Серия «Математическое моделирование и программирование»

Сложное маховое движение ученые свели к суперпозиции двух простых периодических вращений. Первое из них называется взмахом: колебание плоскости крыла вокруг продольной оси аппарата с амплитудой 15 градусов.

Второе называется круткой: вращение крыла вокруг его собственной оси с амплитудой 5–7 градусов. Сдвиг фаз между двумя колебаниями зафиксировали так, чтобы максимум угла крутки приходился на максимум вертикальной скорости крыла.

Скорость набегающего потока задали равной 20 м/с — примерно так летят мигрирующие птицы вроде уток или гусей. Угол атаки, под которым крыло встречает поток, перебирали в трех вариантах: 0, 3 и 6 градусов.

Дальше последовали три серии расчетов. В первой меняли частоту взмахов, во второй — угол атаки, в третьей — амплитуду крутки.

При частоте 4 Гц подъемная сила заметно колебалась, но ее среднее значение, 6,8 Н, совпадало со значением для неподвижной модели. Сила сопротивления в среднем составляли 0,26 Н против 0,34 Н у неподвижного крыла и лишь на коротком интервале периода становились отрицательными, превращаясь в слабую тягу.

При частоте 8 Гц картина качественно изменилась. Амплитуды колебаний выросли, среднее значение подъемной силы осталось прежним — 6,8 Н, а продольная сила изменила знак: вместо сопротивления за период появилась тяга 0,206 Н.

Любопытным оказалось также и то, что максимум горизонтальной тяги по времени совпал с максимумом подъемной силы, и оба экстремума пришлись на опускание крыла.

Работа опубликована в журнале «Вестник ЮУрГУ. Серия „Математическое моделирование и программирование”». Исследование поддержано Минобрнауки России, соглашение №075-03-2024-117. Работа проводилась при финансовой поддержке гранта Российского научного фонда №24-19-00420 (разработка механического прототипа орнитоптера и его цифровой модели).

Изображения: МФТИ

Источник: Минобрнауки России