Научный журнал
Фундаментальные исследования
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,074

ИССЛЕДОВАНИЕ КОЛЕБАНИЙ КВАДРОКОПТЕРА ПРИ ВНЕШНИХ ПЕРИОДИЧЕСКИХ ВОЗДЕЙСТВИЯХ

Попов Н.И. 1 Емельянова О.В. 1 Яцун С.Ф. 1 Савин А.И. 1
1 ФГБОУ ВПО «Юго-Западный государственный университет»
В работе рассмотрена математическая модель, описывающая периодические колебания квадрокоптера, в которой учтены взаимосвязанные электрические и механические процессы в электромеханической системе приводов винтов квадрокоптера. С использованием теоремы динамики об изменении момента количества движения описано вращение корпуса квадрокоптера относительно центра масс. Аналитически найдено решение нелинейных дифференциальных уравнений, описывающих движения квадрокоптера. Проведено моделирование управляемого движения аппарата по углу крена в режиме висения с учетом свойств пропорционального регулятора. Исследовано влияние динамических параметров на характер колебаний квадрокоптера при внешних периодических воздействиях и установлено, что величина собственной частоты колебаний коптера зависит от коэффициента усиления пропорционального регулятора, геометрических размеров коптера, величины управляющего напряжения и аэродинамических характеристик винтовых приводов коптера. Установлено, что существует область параметров коэффициента пропорциональности регулятора, обеспечивающих минимальное значение резонансной амплитуды вынужденных колебаний.
квадрокоптер
линейный регулятор
резонансные режимы
1. Бидерман В.Л. Теория механических колебаний. – М.: Высшая школа, 1980. – 408 с.
2. Браверман А.С., Вайнтруб А.П. Динамика вертолета. Предельные режимы полета. – М.: Машиностроение, 1988. – 280 с.
3. Емельянова О.В., Попов Н.И., Яцун С.Ф. Моделирование движения квадроротационного летающего робота // Актуальные вопросы науки. Материалы VIII Международной научно-практической конференции. – М., Спутник + . 2013. –C. 6–8.
4. Емельянова О.В., Попов Н.И., Яцун С.Ф. Моделирование движения квадрокоптера в пространстве // Авиакосмические технологии (АКТ-2013). Труды XIV Всероссийской научно-технической конференции и школы молодых ученых, аспирантов и студентов. – Воронеж: ООО Фирма «Элист», 2013. – C. 131–138.
5. Загордан А.М. Элементарная теория вертолёта. Военное издательство Министерства Обороны Союза ССР. – М., 1955. – 216 с.
6. Яцун С.Ф., Емельянова О.В., Попов Н.И. Изучение движения квадрокоптера в вертикальной плоскости // Актуальные вопросы технических наук (II): материалы международной заоч. науч. конф. – Пермь: Меркурий, 2013.– С. 66–69.
7. Bresciani Т. Modeling, identification and control of a quadrotor helicopter. Master’s thesis, Department of Automatic control, Lund University, October 2008, p. 170.
8. Tahar М., Zemalache K.M., Omari A. Control of under-actuated X4-flyer using indegral Backstepping controller. Przeglad elektrotechniczny (Electrical review), ISSN 0033-2097, R.87 NR 10/2011, Р. 251–256.
9. Bemporad А., Rocchi C. Decentralized linear time-varying model predictive control of a formayion of unmanned aerial vehicles/ 50th IEEE Conference on decision and control and European control conference, Orlando, FL, USA,1011, Р. 7488–7493.

Последнее время началось бурное развитие нового класса небольших беспилотных квадроротационных летательных аппаратов (квадрокоптеров), способных нести полезную нагрузку в виде контролирующих приборов, видеоаппаратуры. В связи с этим расширяется область применения беспилотных летательных аппаратов (БЛА), увеличивается роль стабилизации их полета. Устойчивость – свойство аппарата восстанавливать без вмешательства оператора кинематические параметры невозмущенного движения и возвращаться к исходному режиму полета после прекращения действия возмущений [1–4, 6]. При отсутствии управления на большинстве режимов полета БЛА являются неустойчивыми, при малейшем внешнем воздействии начинается «неуправляемое снижение», и квадрокоптер падает. Причиной их неустойчивости является отсутствие восстанавливающего момента по отклонению относительно центра масс и малое демпфирование этого движения. Под режимом висения будем понимать способность квадрокоптера держаться в воздухе на одном месте, т.е. висеть [5]. Этот режим используется при необходимости вести какие-либо наблюдения с воздуха, находясь на одном месте. Кратковременное висение БЛА у земли может производиться перед его вертикальной посадкой или после отрыва от земли.

Неустойчивость БЛА вынуждает систему автоматического управления (САУ) непрерывно вмешиваться в управление, чтобы обеспечить устойчивость на режиме висения [4–7]. Поэтому актуальной задачей является исследование динамических явлений, возникающих при внешних периодических воздействиях, а также выбор стратегии управления и параметров регулятора.

Математическая модель квадрокоптера в режиме висения

Режим висения квадрокоптера будем рассматривать в плоскости zСy (рис. 1). Считаем, что т. С, являющаяся центром квадрокоптера, остаётся неподвижной в режиме висения. С этой точкой связана система координат Схyz, тогда подвижная система координат Сх1 y1 z1 связана с корпусом, который поворачивается относительно оси х на угол φ.

pic_6.tif

Рис. 1. Расчётная схема квадрокоптера при для изучения движений относительно оси Сx

За обобщенную координату примем угол поворота φ. Колебания квадрокоптера происходят под действием сил Т2, Т4, момента сопротивления popov01.wmf и момента внешнего периодического воздейсвия М. Винты приводятся во вращение электродвигателями, генерирующими моменты МЭДi, а вращению винтов с угловыми скоростями wi препятствует момент внешних нагрузок МHi (рис. 2).

pic_7.wmf

Рис. 2. Упрощенная схема электропривода винта: 1 – винт; 2 – редуктор; 3 – электродвигатель; u – управляющее напряжение; МЭД – момент, создаваемый электродвигателем; МHi – момент внешних нагрузок винта относительно оси вращения; МВi – крутящий момент винта; JP, JM – момент инерции ротора вокруг оси винта и оси мотора соответственно; wi, Wi – угловые скорости электродвигателя и пропеллера соответственно

Пусть момент внешнего периодического воздействия представлен в виде M = M0sin(w0t), где М0, w0 – амплитуда и частота внешнего возмущения.

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

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

popov02.wmf (1)

где Li, Ri – индуктивность и активное сопротивление обмотки ротора; ui – управляющее напряжение, i – ток в обмотке электромагнита, wi, Wi – угловые скорости электродвигателя и пропеллера соответственно; сЕ – коэффициент пропорциональности, называемый постоянной ЭДС двигателя; Ji – приведенный осевой момент инерции.

Момент, создаваемый электродвигателями:

popov03.wmf (2)

где cMi – коэффициент пропорциональности момента.

Связь между угловой скоростью электродвигателя wi и винта Wi:

wi/Wi = N, (3)

где N – передаточное отношение редуктора.

Момент внешних нагрузок МНi будем определять как

popov04.wmf i = 2, 4. (4)

где popov05.wmf – момент сопротивления вращению винта; m – коэффициент нагрузки, приведенный к валу двигателя; h – кпд электропривода.

Движение корпуса квадрокоптера, относительно оси Сх описывается дифференциальным уравнением

popov06.wmf (5)

где Jx – осевой момент инерции; popov07.wmf – главный момент внешних сил относительно оси Ох:

popov08.wmf (6)

где l – расстояние между центром квадрокоптера и центром пропеллера; bz – аэродинамическая постоянная.

С учетом (6) уравнение (5) будет иметь вид

popov09.wmf (7)

Система дифференциальных уравнений (1), (7) описывает колебания квадрокоптера относительно оси Сх при внешнем периодическом воздействии. В общем случае система является нелинейной и её решение удобно выполнить с помощью численных методов. Однако, сделав некоторые предположения, можно получить аналитическое решение. Пусть электродвигатели, используемые в квадрокоптере, обладают небольшой индуктивностью, т.е. popov10.wmf popov11.wmf, popov12.wmf, popov13.wmf popov14.wmf, тогда уравнения (1) примут вид:

popov15.wmf (8)

где

popov16.wmf popov17.wmf

popov18.wmf popov19.wmf

Откуда

popov20.wmf popov21.wmf (9)

Обозначим:

popov22.wmf popov23.wmf

тогда:

popov24.wmf (10)

С учетом (10) уравнение (7) примет вид:

popov25.wmf (11)

где popov26.wmf.

Пусть управляющее воздействие определяется по следующему алгоритму:

popov27.wmf popov28.wmf (12)

где popov29.wmf – постоянное напряжение питания; k – пропорциональный коэффициент.

Тогда уравнение (11) примет вид:

popov30.wmf (13)

или

popov31.wmf (14)

Получили уравнение вынужденных колебаний, где popov32.wmf popov33.wmf – коэффициент демпфирования; popov34.wmf – циклическая частота свободных колебаний.

Решение дифференциального уравнения (14) имеет вид [1]:

popov35.wmf (15)

Здесь А, В – амплитуды свободных и вынужденных колебаний; g, b – соответственно начальная фаза свободных колебаний и сдвиг фаз вынужденных колебаний по отношению к возмущающей силе.

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

popov36.wmf

popov37.wmf (16)

Рассмотрим предельные случаи:

1. Если l << w0, тогда

popov38.wmf β = 180°. (17)

2. Если l ≈ w0, тогда

popov39.wmf β = 90°. (18)

3. Если l >> w0, тогда

popov40.wmf tgβ = 0. (19)

Анализ уравнений показывает, что при высокочастотных возмущениях (17) в режиме l << w0 на величину амплитуды вынужденных колебаний оказывает влияние коэффициент сопротивления m и частота внешнего возмущения w0, причем с возрастанием последних амплитуда колебаний падает.

При низкочастотных возмущениях (19) в режиме l >> w0 коэффициент сопротивления m и частота внешнего возмущения w0 не оказывают воздействия на величину амплитуды вынужденных колебаний, а влияет осевой момент инерции Jx, повышающий стабильность квадрокоптера и частота свободных колебаний l, которая зависит от геометрических размеров квадрокоптера, аэродинамических коэффициентов, свойств электроприводов, управляющего напряжения, передаточного числа редуктора, коэффициента усиления пропорционального регулятора системы автоматического управления.

При частоте возмущения, равной частоте собственных колебаний системы (l ≈ w0), амплитуда вынужденных колебаний значительно возрастает (18), диссипативные силы уравновешиваются силами инерции при любом значении амплитуды колебаний. Внешняя возмущающая сила оказывается неуравновешенной. Поэтому амплитуда колебаний квадрокоптера, определяемая коэффициентом сопротивления m и частотой свободных колебаний l, в режимах, близких к резонансному, может быть значительной.

Исследование влияния частоты внешнего периодического воздействия на амплитуду вынужденных колебаний для различных коэффициентов усиления пропорционального регулятора системы автоматического управления (САУ) (рис. 3) показывает, что присутствует область частот, при которых амплитуда резонансных вынужденных колебаний корпуса квадрокоптера минимальна.

pic_8.tif

Рис. 3. Зависимость амплитуды вынужденных колебаний от величины коэффициента усиления пропорционального регулятора

Выводы

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

Рецензенты:

Локтионова О.Г., д.т.н., доцент, проректор по учебной работе, ФГБОУ ВПО «Юго-Западный государственный университет», г. Курск;

Кобелев Н.С., д.т.н., профессор, заведующий кафедрой теплогазоснабжения и вентиляции, ФГБОУ ВПО «Юго-Западный государственный университет», г. Курск.

Работа поступила в редакцию 27.01.2014.


Библиографическая ссылка

Попов Н.И., Емельянова О.В., Яцун С.Ф., Савин А.И. ИССЛЕДОВАНИЕ КОЛЕБАНИЙ КВАДРОКОПТЕРА ПРИ ВНЕШНИХ ПЕРИОДИЧЕСКИХ ВОЗДЕЙСТВИЯХ // Фундаментальные исследования. – 2014. – № 1. – С. 28-32;
URL: http://www.fundamental-research.ru/ru/article/view?id=33497 (дата обращения: 18.10.2019).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.074