Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

RESEARCH OF OCSILLATIONS OF QUADROCOPTER UNDER INFLUENCE OF EXTERNAL PERIODIC DISTURBANCE

Popov N.I. 1 Emelyanova O.V. 1 Jatsun S.F. 1 Savin А.I. 1
1 Southwest State University
The mathematical modeling of movement of a quadrocopter which describes the interconnected electrical and mechanical processes in electromechanical system of drives of screws of a quadrocopter is shown in this paper. The movement of quadrocopter body o the center of mass, around an axis is described using the theorem of angular momentum relatively. The analytical solution of the nonlinear differential equations describing movements of quadrocopter body is found. Modeling of roll angle controlled movement of the device in a hanging mode considering the proportional regulator is described. Influence of dynamic parameters on oscillations character of quadrocopter under external periodic disturbances is investigated. It’s also established that the value of proper frequency of oscillations of a quadrocopter depends on gain of the proportional regulator, geometrical sizes of a quadrocopter, level of operating voltage, and aerodynamic characteristics of screw drives of a quadrocopter. It is established that there is an area of parameters of regulator coefficient of proportionality, which are providing the minimum value of amplitude of the forced ocsillations.
quadrocopter
linear regulator
resonance regimes
1. Biderman V.L. Teorija mehanicheskih kolebanij (Theory of mechanical oscillations). M.: The higher school, 1980. рp. 408.
2. Braverman A.S., Vajntrub A.P. Dinamika vertoleta. Predel’nye rezhimy poleta («Dinamika of the helicopter». Limit modes of flight). M.: Mechanical engineering, 1988. рp. 280.
3. Emelianova O.V., Popov N.I., Jatsun S.F. Modelirovanie dvizhenija kvadrorotacionnogo letajushhego robota. Aktual’nye voprosy nauki. Materialy VIII Mezhdunarodnoj nauchno-prakticheskoj konferencii. («Modeling of movement of the quadrotor helicopter flying robot». Topical issues of science. Materials VIII of the International scientific and practical conference.) Moscow, Satellite + . 2013. pp 6–8.
4. Emelianova O.V., Popov N.I., Jatsun S.F. Modelirovanie dvizhenija kvadrokoptera v prostranstve. Aviakosmicheskie tehnologii (AKT-2013). Trudy XIV Vserossijskoj nauchno-tehnicheskoj konferencii i shkoly molodyh uchenyh, aspirantov i studentov. («Modeling of movement of a quadrotor helicopter in space». Aerospace technologies (AKT-2013). Works XIV of the All-Russian scientific and technical conference and school of young scientists, graduate students and students). Voronezh: JSC Firm Elist, 2013. pp. 131–138.
5. Zagordan A.M. Jelementarnaja teorija vertoljota. («Elementary theory of the helicopter»). Military publishing house of the Ministry of Defence of USSR. Moscow, 1955. рp. 216.
6. Jatsun S.F., Emelianova O.V., Popov N.I. Izuchenie dvizhenija kvadrokoptera v vertikal’noj ploskosti. Aktual’nye voprosy tehnicheskih nauk (II): materialy mezhdunarodnoj zaoch. nauch. konf. («Studying of movement of a quadrotor helicopter in the vertical plane». Topical issues of technical science (II): materials of the international correspondence scientific conference). Perm, Mercury, 2013. pp. 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.