В характеристиках дронов эти производные даются обычно по относительным угловым скоростям вида
x 2lV x ;
y 2lV y.
Для перехода от величин x и y следует воспользоваться соотношениями
x 2 laV x ;
y 2 laV y .
Аналогично имеем:
|
m y m |
|
|
|
y |
l |
|
|
|
; |
|
|
|
|
2 V |
|
x |
x |
|
|
|
|
|
|
|
|
|
a |
|
|
my x m |
|
|
x |
|
l |
|
|
y |
; |
|
|
2 aV |
|
|
|
|
|
|
|
|
|
|
my y m |
|
y |
l |
|
|
y |
. |
|
2 aV |
|
|
|
|
|
|
|
|
|
Для возмущений f1, f2 и f3 |
имеем выражения |
|
|
|
|
|
|
|
|
' |
|
|
f1 p z f1 |
|
|
|
|
G |
|
|
l2 |
' |
|
f2 |
|
|
|
|
|
|
|
f2 |
; |
ρ0SV02 |
|
|
l |
|
|
|
|
|
|
|
|
|
P |
|
|
l3 |
|
|
|
f3 |
|
|
|
|
f3' |
, |
|
2 |
|
|
|
|
|
|
|
l |
|
|
|
|
ρ0SV0 |
|
|
|
|
где f1' , f2' , f3' – возмущения, вызванные, например, ударными вол-
нами вблизи БПЛА, взрывами и т. д; l2 – плечо момента крена;
l3 – расстояние между двигателями, имеющими разные тяги.
Уравнения (13.5) и (13.6), устанавливающие связь между регулируемыми величинами β, γ, φ, ψ, x , y и регулирующими факто-
рами δэ и δн и характеризующие динамические свойства БПЛА
в боковом движении, описывают линейную математическую модель бокового движения летательного аппарата.
В случае горизонтального полета 0 с учетом малости угла ψ из уравнений получаем
ωx γ; ωy ψ
или в относительной форме
x p ;
Уравнения (13.5) принимают вид
p n11 β n12 p n14 n21β p n22 p n23 n31β n32 pγ p n33
y p .
n13 pψ f1; |
|
|
|
|
pψ nэδэ f2; |
(13.7) |
pψ n δ |
н |
f |
3 |
. |
|
н |
|
|
|
В табл. 13.1 даны ориентировочные значения коэффициентов nik
уравнений бокового движения легкого самолета для двух режимов полета.
|
|
|
Таблица 13.1 |
Ориентировочные значения коэффициентов nik |
|
|
|
|
|
|
|
БПЛА |
|
Коэффициент |
H 6 км, |
|
H 10 км, |
|
M 0,77, |
|
M 0,8, |
|
|
|
|
|
a 2,5 c |
|
a 3,8 c |
|
n11 |
0,156 |
0,097 |
|
n12 |
0 |
0 |
|
n13 |
–1 |
|
–1 |
|
n14 |
–0,039 |
|
–0,039 |
|
n21 |
15,8 |
9,5 |
|
n22 |
6,7 |
4,82 |
|
n23 |
0,43 |
0,41 |
|
n31 |
5,76 |
4,3 |
|
n32 |
0,037 |
0,0058 |
|
n33 |
0,22 |
0,16 |
|
nэ |
30,7 |
19 |
|
nн |
3,18 |
2,26 |
|
13.2. Передаточные функции конкретных моделей бокового движения беспилотного летательного аппарата
Рассмотрим частные случаи уравнений бокового движения. Простейшим боковым движением БПЛА является движение рыскания без крена, когда в силу большой инерции можно пренебречь движением центра масс под действием боковых сил. При таком движении продольная ось БПЛА совершает колебания относительно вектора
скорости, поворот которого не учитывается. Опуская первое и второе уравнения (13.7) и полагая в последнем уравнении 0 , получим
n31β p2 n33 p ψ nнδн f3. |
(13.8) |
В этом уравнении при указанных выше предположениях угол
есть угол между вектором скорости и продольной осью БПЛА. Следовательно, этот угол равен углу скольжения, т. е. ψ = β, и уравне-
ние (13.8) принимает вид
p2 n33 p n31 ψ nнδн f3.
Передаточная функция будет
Следующим частным видом бокового движения является плоское движение со скольжением при неизменном угле крена. Полагая в уравнениях (13.7) γ 0 и опуская второе уравнение, получим
p n |
β pψ f ; |
|
|
11 |
1 |
|
(13.9) |
p2 n33 p ψ n31β nнδн f3 |
|
. |
|
|
|
|
|
При выводе этих уравнений предполагается, |
что угол атаки |
близок к нулю, поэтому sin α 0 |
и cos α 1. |
|
|
Из уравнений (13.9) получаем выражения для передаточных
функций: |
|
|
|
|
|
|
|
|
|
|
W |
p |
|
|
|
n |
|
; |
|
|
|
2 |
2 |
|
|
|
β |
|
|
p |
|
|
|
|
|
|
|
2d0 0 p 0 |
|
|
(13.10) |
|
p |
|
|
|
nн p n11 |
|
|
|
W |
|
|
|
|
|
, |
|
p2 2d0 0 p 02 |
p |
|
ψ |
|
|
|
|
где
2d0 0 n11 n33,
02 n31 n11n33.
Из выражений (13.10) видно, что БПЛА является нейтральным по отношению к углу рыскания ψ.
Рассмотрим случай бокового движения, реализуемого приближенно в начальный момент крена БПЛА, при этом пренебрегаем
изменением курса ψ 0 . В этом случае уравнения (13.7) принимают вид (третье уравнение опускаем)
p n |
β n |
p n |
ψ f ; |
|
11 |
12 |
14 |
1 |
|
n21β p2 n22 p γ nэδэ f2 |
|
. |
|
|
|
|
|
Отсюда находим выражения для передаточных функций, полагая
f1 f2 |
0 : |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
W |
|
p |
nэ p n11 |
; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
γ |
|
|
|
|
|
|
p |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n |
n |
|
p n |
|
|
|
|
|
|
|
W |
p |
э |
|
12 |
|
14 |
|
|
, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ψ |
|
|
|
|
|
|
|
p |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где |
|
|
p p3 c' |
|
|
|
|
|
|
|
|
|
|
|
|
|
p2 |
c' |
p c' |
; |
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
2 |
|
|
|
3 |
|
|
|
|
c' |
n |
n ; |
|
c' |
n n |
|
n n ; |
c' |
n n . |
|
1 |
11 |
22 |
|
2 |
11 |
22 |
|
12 |
21 |
|
|
|
3 |
21 |
14 |
Рассмотрим также случай движения по крену без скольжения. В этом случае из системы (13.7) при пренебрежении моментом рыскания находим
p2 n22 p γ n23 pψ nэδэ.
13.3. Автоматическое управление угловыми движениями летательного аппарата
Вектор скорости центра масс БПЛА как твердого тела в каждый момент времени направлен по касательной к траектории полета. Следовательно, для изменения направления траектории полета в пространстве необходимо менять направления вектора скорости. Величина и направление вектора скорости определяются величиной и направлением вектора тяги двигателя и аэродинамическими характеристиками БПЛА. Вектор тяги обычно ориентируется по отношению к корпусу БПЛА, и на большинстве режимов полета его направление остается неизменным или меняется в ограниченных пределах.
Аэродинамические характеристики БПЛА (подъемная сила, боковая сила, сила лобового сопротивления и др.) зависят от углового положения вектора тяги по отношению к вектору скорости. Так, например, для изменения подъемной или боковой сил необходимо изменять углы между векторами тяги и скорости полета в вертикальной (угол атаки) или горизонтальной (угол скольжения) плоскостях.
Таким образом, вектор тяги на данном режиме полета должен занимать определенное положение по отношению к вектору скорости. Поскольку вектор тяги связан с осями БПЛА, то БПЛА должен занимать определенное угловое положение по отношению к вектору скорости, различное в разных условиях полета. Задание определенных угловых положений, осуществляемое управляющей системой, составляет основную задачу управления угловыми движениями БПЛА.
При этом угловые координаты и их производные, характеризующие движение БПЛА по отношению к центру масс, остаются неизменными или меняются по определенным законам. Обычно производится управление углами тангажа, крена, рыскания, атаки, скольжения и их производными.
При построении САУ (систем автоматического управления) угловыми движениями БПЛА необходимо изучать динамику переходных процессов, влияние структуры и параметров систем управления на качество процессов, а также рассматривать вопросы синтеза систем управления для получения заданного качества процесса управления.
13.3.1.Управление угловой скоростью дрона
Всистемах автоматического управления полетом имеются контуры управления угловыми скоростями БПЛА, служащие для формирования демпфирующих моментов и, следовательно, для улучшения качества переходного процесса. Для этих же целей применяются самостоятельные системы управления угловыми скоростями, назы-
ваемые демпферами.
Рассмотрим систему управления угловой скоростью тангажа, представленную структурной схемой на рис. 13.4. В этой схеме закон управления принят в виде
δв k 3 ,
где k – передаточное число;
3 – заданная угловая скорость тангажа, а передаточная функция БПЛА по угловой скорости представлена в прямоугольнике.
Рис. 13.4. Схема управления угловой скоростью тангажа
Передаточная функция системы, какследуетизрис. 13.4, имеетвид
W p |
k&n |
p n |
|
|
|
|
в |
22 |
|
. |
(13.20) |
p |
2 |
|
|
2 |
|
|
2d p |
|
Здесь
2d 2d0 0 nвk ;
(13.21)
2 02 nвn22k ,
где 0 и d0 – соответственно собственная частота и коэффициент
затухания движений БПЛА.
Видно, что передаточная функция системы управления угловой скоростью тангажа равна произведению передаточных функций колебательного и форсирующего звеньев. Для выбора передаточно-
го числа k заметим, что наилучшее качество процесса в колеба-
тельном звене получается при d 0,7 / 1, причем верхнее значение следует брать при наличии форсирования. Исключая из уравнений (13.12) частоту , получим выражение для k :
|
|
|
k |
2d 0 |
|
2d |
|
d |
n22 |
|
|
|
|
|
|
|
|
|
|
|
1 |
0 |
. |
|
|
(13.13) |
|
|
|
|
|
|
|
|
|
|
|
|
nв |
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
В табл. 13.2 приведены значения передаточного числа k |
и соб- |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ственной частоты 0 |
для контуров управления угловой скоростью |
дронов № 1 и 2, параметры которых взяты из табл. 13.2. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблица 13.2 |
|
|
|
|
Значения k |
|
и |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2d |
|
|
|
|
|
|
|
|
k |
|
|
|
|
|
|
|
, |
|
|
nв |
n11 |
|
|
|
a , |
|
|
|
|
k & |
–1 |
|
|
с |
Самолет |
|
|
|
|
|
|
|
|
с |
|
|
|
|
, |
|
|
|
0 |
|
|
|
|
|
|
|
|
|
град/с |
с |
|
р |
|
|
№ 1 |
49 |
2,4 |
|
5,25 |
|
43,9 |
|
3,8 |
|
0,154 |
|
0,585 |
2,07 |
|
3 |
|
№ 2 |
24,5 |
2,6 |
|
4,94 |
|
15,13 |
|
2,1 |
|
0,128 |
|
0,28 |
2,3 |
2,7 |
|
Заметим, что относительные значения передаточного числа k
и собственной частоты 0 вычислены по формулам (13.12) и (13.13), а размерные – определены из соотношений
k' k a ,
0 .a
Время регулирования в обоих случаях не превышает 3 с.
При рассмотрении системы управления угловой скоростью рыскания возьмем структурную схему (рис. 13.5), подобную схеме, показанной на рис. 13.4.
Рис. 13.5. Схема управления угловой скоростью рыскания
Передаточная функция этой системы может быть представлена в виде
W p |
k |
n |
|
p n |
|
|
|
|
|
н |
11 |
|
, |
(13.23) |
p |
2 |
|
|
|
2 |
|
|
2d p |
|
|
где
2d 2d0 0 k nн;
2 02 k nнn11.
Поскольку структура передаточной функции (13.14) аналогична структуре передаточной функции (13.11), то передаточное число
k можно определить из формулы (13.14), в которой параметры
продольного движения дрона следует заменить параметрами бокового движения.
В табл. 13.3 приведены значения передаточного числа k и соб-
ственной частоты 0 для контуров управления угловой скоростью рыскания № 1 и 2; коэффициент относительного затухания d 1 .
Таблица 13.3
|
|
|
Значения k и |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2d |
|
2 |
|
|
|
k |
|
|
|
|
|
|
|
, с |
Самолет |
nн |
n11 |
τa , с |
|
|
|
|
k |
–1 |
р |
|
0 0 |
0 |
|
|
|
|
град/с |
, |
с |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
№ 1 |
3,18 |
0,156 |
0,376 |
5,8 |
2,5 |
|
1,37 |
|
3,42 |
1,02 |
6,15 |
№ 2 |
22,5 |
0,26 |
1,05 |
54 |
2,1 |
|
0,29 |
|
0,61 |
3,34 |
1,78 |
Из сравнения передаточных чисел k |
и |
k |
|
видно, что требуемое |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
значение передаточного числа в канале рыскания значительно больше, чем в канале тангажа. Объясняется это малым естественным демпфированием дрона по углу рыскания.
При рассмотрении систем управления угловой скоростью тангажа и рыскания предполагалось, что законы управления реализованы точно. Рассмотрим теперь систему управления угловой скоростью крена, причем будем полагать, что динамические погрешности сервопривода значительны.
На рис. 13.6 динамика сервопривода отражена внутренним контуром, а динамика управления угловой скоростью – внешним контуром. Передаточная функция БПЛА по угловой скорости крена взята в виде инерционного звена.
Рис. 13.6. Схема управления угловой скоростью крена
На рис. 13.7 показана одноконтурная структурная схема, преобразованная из схемы, приведенной на рис 13.6. В этой схеме параметры передаточной функции сервопривода определяются выражениями