МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ДВИЖЕНИЙ МАШУЩЕГО КРЫЛА ПРИ ВЗАИМОДЕЙСТВИИ С СОПРОТИВЛЯЮЩЕЙСЯ СРЕДОЙ
Аннотация и ключевые слова
Аннотация (русский):
Цель исследования. В статье рассматриваются вопросы моделирования биологически инспирированного движения беспилотной летающей платформы (БЛП) оснащенной машущими крыльями. Особое внимание уделено моделированию взаимодействия машущего крыла с сопротивляющейся воздушной средой. Поэтому, целью исследования является математическое моделирование приведенной силы, возникающей при взаимодействии крыла с сопротивляющейся средой. Задача, решению которой посвящена статья. Для выполнения поставленной цели необходимо провести математическое моделирование вращательного движения крыла БЛП с учетом кинематических и динамических особенностей двухкоординатного вращения для различных законов изменения углов поворота крыла относительно продольной и поперечной осей. Методы исследования. Для решения поставленной задачи использовались методы классической механики, теории устойчивости движения, оптимального управления, теории наблюдений, параметрической оптимизации. Новизна работы заключается в алгоритмах формирования законов изменения углов поворота крыльев относительно продольной и поперечной осей, представленных в виде гладких функций. Предложен метод определения приведенной силы, возникающей при движении крыла. Получены зависимости проекций приведенной силы взаимодействия крыла с воздухом, обеспечивающие создание подъемной и тяговой сил. Эти силы позволяют БЛП перемещаться в пространстве по заданной траектории. Результаты исследования. Проведено моделирование режимов движения машущего крыла и получены зависимости мгновенных значений проекций сил, создаваемых крылом при движении. Эти силы позволяют БЛП перемещаться в пространстве по заданной траектории. Показано, что управление движением БЛП осуществляется за счет изменения частоты и амплитуды колебаний крыла, а также изменением угла наклона нейтральной линии относительно которой происходит поворот крыла. Выводы: предложена модель приведенной силы, возникающей при движении машущего крыла. Построены графики зависимостей эффективной площади крыла (площади миделева сечения) для различных параметров законов изменения углов поворота от времени. Проведено моделирование машущего движения крыльев летательного аппарата. Определена область параметров, при которых проекции перемещения крыла создают горизонтальную F1x, и подъёмную F1z силы при различных значениях проекций площадей S1х, S1y, S1z.

Ключевые слова:
машущее крыло, подъёмная сила, угол поворота крыла, эффективная площадь крыла
Текст
Текст произведения (PDF): Читать Скачать

Введение

 

Развитие научных технологий, появление новых современных легких и прочных материалов и электронных компонентов приводит к возрождению интереса к моделированию движения летательных аппаратов с машущим крылом (ЛАМК), которые имитируют полет птиц или насекомых [11-13]. Такие механические устройства обладают возможностью вертикального взлета, планирования, зависания над объектами, а также они достаточно маневренны, как их биологические объекты [11-13, 17]. Интерес к идее копирования машущего полета природных аналогов - очень велик, однако его реализация связана с особенностями аэродинамической теории и отсутствием опыта разработки механизмов-приводов и их применения [1-12].

В природе крыло биологического объекта, например, птицы, реализует три вида вращения: вокруг оси, проходящей вдоль корпуса (взмахи вверх-вниз), поперечной оси корпуса (приводит к наклону крыла с целью изменения его угла атаки) и вокруг вертикальной оси, т.е. взмахи крыльев вперед и назад [11-14,17]. Помимо этих вращений, крыло птиц может изменять свою площадь во время взмаха для управления силой сопротивления в воздушной среде и развития подъёмной силы [1-6]. В отличие от птиц, у большинства насекомых площадь крыла не изменяется во время взмаха, поэтому, в данной статье остановимся на изучении движения недеформируемого крыла насекомого [13,14,16].

 

  

Материалы, модели, эксперименты и методы

 

1. Моделирование движения машущего крыла

При моделировании машущего крыла летательного аппарата будем рассматривать его движение в абсолютной системе координат OXYZ (рис. 1). Свяжем с корпусом робота относительную, подвижную систему координат C2X2Y2Z2, начало которой совпадает с центром тяжести корпуса С2. Ось C2X2 - направлена параллельно продольной оси корпуса, ось C2Y2 - направлена перпендикулярно плоскости C2X2Z2, а ось C2Z2 - перпендикулярно плоскости C2X2Y2. Плоскость C2X2Z2 является плоскостью симметрии аппарата [1-6, 13,17].

Для описания движения крыльев 1,3 свяжем плоскостью крыла c системой координат С2XiYiZi (i=1,3), где С2 -  центр масс фюзеляжа (корпуса). Оси С2Xi совпадают с осями поперечного вращения, а ось СiYi  принадлежит плоскостям крыльев, совпадают с её кромкой и с продольной осью вращения крыла. Таким образом, системы координат С2XiYiZi (i=1,3) могут вращаться относительно корпуса вместе с крыльями, совершая двухкоординатное вращение,  определяемое углами ai , относительно осей CiYi и ji, относительно осей СiXi .

Будем рассматривать крылья как два твердых тела, которые сочетают движения по углам φi и αi, вносящие значительный вклад при формировании подъемной и тяговой сил.

Пусть  ,  или                               .                     (1)

При моделировании необходимо учесть важные кинематические параметра крыла, такие как амплитуда взмаха, частота и углы поворота ai и ji, в зависимости от положения крыла [13-17].

 

Рис.1. Схема систем координат, связанных с крыльями: 1,3 – крылья, 2 – фюзеляж (корпус)

Fig. 1. Diagram of the coordinate systems associated with wings: 1,3 - wings, 2 - fuselage (hull)

 

 

2. Кинематическая модель машущего крыла

Для удобства, будем рассматривать только крыло 1, положение которого определено поворотом на угол a1 и последующим поворотом на угол j12 (рис. 2) [13, 14]. Вращение крыла относительно оси С2Y1 происходит с угловой скоростью , а относительно оси C2X1 -с   угловой скоростью .

 

 

Рис. 2. Схема поворота системы координат крыла 1: ` - скорость центра

масс С1   крыла; а, б – геометрические параметры

Fig. 2. The rotation scheme of the wing coordinate system 1:  is  speed

of the center of masses C1 of a wing; a, b - geometric parameters

 

 

Расстояние между точками С2 и С1 вдоль оси C2X2 равно а, вдоль оси C2Z2 - равно нулю, а вдоль оси C2Y2b; а и b – переменные величины, зависящие от положения точки С1. Для приближенных расчетов примем их постоянными величинами.

Положение центра масс С1, во второй (2) системе координат определяется радиус-вектором:

.                       (2)

Матрицы, определяющие векторы, заданные в первой (1) системе координат, с учетом двух вращений, имеют вид:

 

.                       (3)

На данном этапе рассмотрим законы изменения углов поворота  и  при движении крыльев в виде гладких функций:

;                ,                      (4)

;         ;           ,                            (5)

 

где  - отклонение от нейтральной линии крыла (ось Х2)  и амплитуда  колебаний крыла, относительно продольной оси (рис. 3); g - угол сдвига фазы; - отклонение от нейтральной линии крыла (ось Y2) и амплитуда колебаний крыльев, относительно поперечной оси (рис. 4); w- угловая скорость.

 

а)                                                             б)                                         в)                                                           г)

Рис. 3. Схема расположения крыла относительно продольной оси при следующих параметрах:

а- б -  в - г -  

Fig. 3. The layout of the wing relative to the longitudinal axis with the following parameters: а- б -  в - г -  

 

а)                                                                            б)

Рис. 4. Схема возможного расположения начального положения

 крыла 1 относительно поперечной оси: 

а - при ;  б - при  

Fig. 4. Diagram of the possible location of the initial position

of the wing 1 relative to the transverse axis:

а - при ;  б - при  

 

 

На рис. 5, 6 представлены графики, иллюстрирующие изменение углов  и  во времени для различных схем расположения крыла.

 

Рис. 5. Зависимости углов поворота  крыльев от времени:

1,2,3,4 соответствуют схемам а,б,в,г  расположения

крыла рис. 3 соответственно

Fig. 5. The dependences of the angles of rotation of the wings on time:

1,2,3,4 correspond to the schemes а,б,в,г of the wing arrangement

of figure 3, respectively

 

а)                                                           б)

Рис. 6. Зависимости углов поворота  крыльев от времени:

 а - для схемы на  рис. 4,а;  б - для схемы на рис. 4,б

Fig. 6. The dependence of the angles of rotation of the wings on time:

a - for the scheme in figure 4,a; b- for the scheme in figure 4,b

 

 

Анализ графиков на рис. 5, 6 показывает, что существует сдвиг по фазе углов поворота крыльев в зависимости от их начального (нулевого) и последующего отклонения, относительно продольной и поперечной осей. Выявленные факторы можно использовать, для моделирования управления движением машущего крыла.  

При движении крыла относительно корпуса вверх, его плоскость может поворачиваться на угол  и становиться практически вертикальной (рис.3, в). В этот момент площадь миделевого сечения S близка к нулю. При движении крыла вниз плоскость крыла поворачивается и может становиться горизонтальной (рис. 3, а,б,г), при этом площадь миделева сечения достигает своего максимального значения. Соответственно, это приводит к увеличению силы сопротивления. Скорость изменения силы сопротивления зависит от того какой закон изменения площади миделева сечения принят для управления движением.

Площади проекций крыла 1 на плоскости, перпендикулярные осям C2X2, C2Y2, C2Z2, представим как (рис. 7):

 

                               ;                            (6)

где S0 - площадь сечения крыла.

 

                                                               а)                                                                           б)

Рис. 7. Расчетные зависимости проекций  площади крыла S, S, S1z  при периодическом движении крыла вверх-вниз от времени при частоте w=25,12 рад/с; S0=0,36´0,4=0,48 м2 и следующих параметрах:

а - ; ;   б - ;

Fig. 7. Calculated dependences of the wing area projections  S1х, S1х, S1z  with periodic up-and-down wing movement on time at a frequency of w=25.12 rad/s; S0=0.36×0.4=0.48 m2 and the following parameters:

 а - ; ;   б - ;

 

 

3. Определение силы взаимодействия крыла с воздухом

Рассмотрим схемы внешней распределенной силы взаимодействия с сопротивляющейся средой`F, приложенной  в точке С1, при движении крыльев вверх (рис. 8, а) и вниз (рис. 8, б).

 

 

                            

                                               а)                                                                                          б)

Рис. 8. Схема крыла робота при относительном движении: а – вниз; б – вверх;

1 - нейтральная линия крыла относительно корпуса; 2 – крыло; 3 – корпус

Fig. 8. The diagram of the robot wing with relative motion: a - down; b - up;

1 - the neutral line of the wing relative to the body; 2 - the wing; 3 - the body

 

 

Здесь приняты следующие обозначения  - проекции скорости центра масс Сна оси связанной с корпусом системы координат, `F1 приведенная сила взаимодействия крыла с воздухом.

 

 В общем случае: , или .

Силу   представим в виде вектора:

,                                                          (7)

где F1x, F1y, F1z – проекции силы крыла 1 на оси подвижной второй (2) системы координат:

;

;                                              (8)

,

 

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

Анализ формул (8) показывает, что продольная , поперечная   и подъемная  силы крыла изменяются в зависимости от скорости движения крыла и модуля его площади [11-14].

Уравнения проекции скоростей центра масс крыла .на оси связанной, первой (1)  системы координат, полученные на основании кинематических уравнений, изложенных в [14,17],   будут иметь вид:

 

 

;                                                    (9)

 

Результаты

 

На основании уравнений (8) с учетом (6) и (9) проведем моделирование проекций силы в зависимости от времени при следующих параметрах a=0,2 м; b=0,6 м; w=25,12 рад/с, С=0,5; r=1,29 кг/м3.

Определим средние за период значения продольной , поперечной  и подъемной   сил крыла для графиков рис. 9 (таблица):

 

;       ;       .                  (10)

                                                                                                Таблица

Средние значения проекции силы

Рис. 9

а

0,216

0

0

б

0,317

1,133

0,407

в

0,039

1,694

-0,918

 

 

Анализ графиков на рис. 9, а так же средние за период значения силы  крыла показывают, что характер изменения проекций силы существенно зависит от параметров, задающих движение крыла ; Изменяя эти постоянные, можно получить силы, обеспечивающие подъем летательного аппарата (рис. 9, б), движение по горизонтальной прямой (рис. 9, а) и спуск (рис. 9,в).

 

а)

б)

в)

Рис. 9. Зависимости проекций силы  взаимодействия крыла с воздухом от времени        при следующих параметрах: а - ; ; б - ; ;    

в - ;

Fig. 9. The dependences of the projections of the force  of the interaction of the wing with the air on time with the following parameters: а - ; ; б - ; ;    

в - ;

Заключение


В результате исследований предложена модель приведенной силы, возникающей при движении машущего крыла. Построены графики зависимостей эффективной площади крыла (площади миделева сечения) для различных параметров законов изменения углов поворота от времени. Проведено моделирование машущего движения крыльев летательного аппарата. Определена область параметров, при которых проекции перемещения крыла создают горизонтальную F1x, и подъёмную F1zx силы при различных значениях проекций площадей S, S1y, S1z.

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

1. Jatsun S., Efimov S., Emelyanova O., Leon A.S.M., Davalos P.J.C. Modeling and control architecture of an autonomous mobile aerial platform for environmental monitoring. In 2019 International Conference on Information Systems and Computer Science (INCISCOS). IEEE, 2019. Р. 177-182. DOI:https://doi.org/10.1109/INCISCOS49368.2019.00036.

2. Jatsun S.F., Korenevskiy N.A., Efimov S.V., Korovin E.N. An Automated System for Monitoring the Environment and Assessing People's Status in Extreme Situations Using a Flying Robot. Biomedical Engineering. 2018;52(4). DOI:https://doi.org/10.1007/s10527-018-9832-x.

3. Karimian S., Jahanbin Z. Bond graph modeling of a typical flapping wing micro-air-vehicle with the elastic articulated wings. Meccanica. 2020;55(6): 1263-1294. DOI:https://doi.org/10.1007/s11012-020-01162-w

4. Nguyen K., Au L.T.K., Phan H.V., Park S.H., Park H.C. Effects of wing kinematics, corrugation, and clap-and-fling on aerodynamic efficiency of a hovering insect-inspired flapping-wing micro air vehicle. Aerospace Science and Technology. 2021;118: 106990. DOI:https://doi.org/10.1016/j.ast.2021.106990.

5. Smart bird. URL: https://ru.wikipedia.org/wiki/SmartBird.

6. SmartBird. URL: https://www.festo.com/group/ru/cms/10238.htm.

7. Vorochaeva L.Y., Efimov S.V., Loktionova O.G., Yatsun S.F. Motion Study of the Ornithopter with Periodic Wing Oscillations. Journal of Computer and Systems Sciences International. 2018;57(4):672-687. DOI:https://doi.org/10.1134/S1064230718040147.

8. Word Science. URL: http://wordscience.org/uchenye-garvarda-sozdali-letayushhego-robota-robo-fly.html.

9. Yoon S.H., Cho H., Lee J., Kim C., Shin S.J. Effects of camber angle on aerodynamic performance of flapping-wing micro air vehicle. Journal of Fluids and Structures. 2020;97:103101. DOI:https://doi.org/10.1016/j.jfluidstructs.2020.103101.

10. Емельянова О.В., Поляков Р.Ю., Ефимов С.В., Яцун С.Ф. Мобильный летательный комплекс для раннего обнаружения очагов возгорания. Фундаментальные и прикладные проблемы техники и технологии, 2018;3(329):136-141. ISBN: 978-5-9929-0514-4.

11. Ефимов С.В. Коршунов Е.В., Поляков Р.Ю., Тарасов О.С., Яцун С.Ф. Экспериментальные исследования параметров машущего крыла для прыгающе-летающего робота. Фундаментальные исследования. 2015;2(11):252-256. ISSN 1812-7339.

12. Поляков Р.Ю., Яцун С.Ф. Моделирование полета инсектоптера для исследования параметров машущего крыла. Системы управления, сложные системы: моделирование, устойчивость, стабилизация, интеллектуальные технологии. 2020. С. 431-436. ISBN: 978-5-00151-153-3.

13. Поляков Р.Ю., Яцун С.Ф. Математическое моделирование полета стрекозоподобного миниробота. Системы управления, сложные системы: моделирование, устойчивость, стабилизация, интеллектуальные технологии. 2020. С. 467-473. ISBN: 978-5-00151-153-3.

14. Поляков Р.Ю. Мобильная приборная платформа для системы экологического мониторинга загрязнения токсичными газами атмосферного воздуха: специальность 05.11.13 «Приборы и методы контроля природной среды, веществ, материалов и изделий»: дис. на соискание научной степени канд. тех. наук / Поляков Роман Юрьевич; Юго-Западный гос. ун-т. Курск, 2019. 165 с. Библиогр.: с.154-165. https://www.dissercat.com/

15. Поляков Р. Ю., Ефимов С. В., Яцун С.Ф. Робот-инсекоптер для мониторинга окружающей среды Современные проблемы гражданской защиты. 2015. №. 3 (16). ISSN 2658-6223.

16. Поляков Р.Ю., Ефимов С.В., Валуйский В.Е., Мозговой Н.В. Проведение исследования движения многозвенного летающего робота для мониторинга окружающей среды. Пожарная безопасность: проблемы и перспективы. 2014;1(1 (5)). ISSN 2411-3778.

17. Яцун С.Ф., Локтионова О.Г., Ворочаева Л.Ю., Емельянова О.В. Конструкция и система управления робота-орнитоптера, оснащенного крыльями и хвостом. Известия Юго-Западного государственного университета. 2018: 22(2):18-26. DOIhttps://doi.org/10.21869/2223-1560-2018-22-2-18-26.

Войти или Создать
* Забыли пароль?