Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Мехатроника и динамика мини-роботов

.pdf
Скачиваний:
2
Добавлен:
24.11.2025
Размер:
4.17 Mб
Скачать

с2 ik ;

2dc c 1 ; k k .

Рис. 13.7. Схема преобразования из схемы рис. 13.6

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

Передаточная функция системы (см. рис. 13.6 и 13.7) будет

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n k k

 

 

 

 

 

 

 

 

 

 

 

 

э

 

 

 

 

 

 

.

(13.15)

 

 

 

2

2d

 

p ω

2

 

p n

 

 

 

 

3

 

p

 

 

 

 

 

 

 

 

 

 

c

c

 

c

22

э

 

 

 

Из структуры выражения (13.15) видно, что для улучшения качества процесса необходимо увеличивать собственную частоту серво-

привода c . В пределе, при достаточно большой частоте c , ди-

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

увеличение передаточного числа kγ приводит к сокращению вре-

мени переходного процесса.

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

350

13.3.2.Управление ускорением центра масс дрона

Вряде случаев системы управления полетом БПЛА включают контуры управления перегрузками, т. е. ускорениями центра масс.

Такие контуры управления являются внешними по отношению к контурам управления угловыми скоростями.

Поскольку нормальная и боковая перегрузки пропорциональны соответственно углам атаки и скольжения, то контуры управления перегрузкой эквивалентны контурам управления углами α и β .

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

Рассмотрим нормальные и боковые ускорения, которые определяются скоростью изменения углов наклона траектории в вертикальной и горизонтальной плоскостях:

J&y Vθ&,

J&z Vθ&,

 

б

где θ α; θб ψ β.

Если воспользоваться первыми из уравнений (13.7), то (при

условии f2 0 и

f1 0 ) получим

 

 

 

 

 

J

y

n V ;

J

z

Vn .

(13.16)

 

 

22

 

11

 

Передаточные функции самолета по ускорениям J y и Jz при

входах δв и δэ можно получить из передаточных функций и соот-

ношений (13.16):

 

J y

 

 

n n

V

 

 

 

 

 

в

22

 

;

 

δв

 

p2 2d

 

p 2

 

 

 

0

0

0

 

 

 

J&

 

 

n n V

.

 

z

 

э

11

 

 

 

δэ

 

 

p2 2d0ω0 p ω02

 

 

Рассмотрим теперь схему системы стабилизации нормального ускорения (рис. 13.8). Для измерения ускорения J y и Jz можно

351

применить акселерометры, оси чувствительности которых совпадают соответственно с нормальной Oy и поперечной Oz осями аппарата.

Рис. 13.8. Схема системы стабилизации нормального ускорения

Запишем закон управления в соответствии со схемой рис. 13.8 в виде

δв k jy jy0 .

(13.17)

Уравнение движения системы будет

p2 2d0 0 p 02 kn22nвV jy kn22nвVjy0.

Отсюда видно, что при принятом законе управления (13.17) обеспечивается повышение запаса статической устойчивости БПЛА. Такой же эффект будет при введении в закон управления угла ата-

ки вместо нормального ускорения &jy . Демпфирование же БПЛА

будет неудовлетворительным. Для получения приемлемого переходного процесса контур управления ускорением необходимо применять совместно с контуром управления угловой скоростью. Сказанное здесь об управлении нормальным ускорением полностью относится и к управлению боковым ускорением.

13.3.4. Управление углом тангажа посредством статического автопилота

Рассмотрим статическую систему автоматического управления углом тангажа (рис. 13.9), включающую контур управления угловой

352

скоростью и контур управления углом тангажа. Передаточная функция БПЛА взята в предположении постоянства скорости самолета. На структурной схеме не показаны внешние возмущения f2 и f3 , дей-

ствующиена БПЛА. Закон управления системы берем в виде

δ

в

k

 

 

k p ,

(13.18)

 

 

3

 

 

где 3 – заданное значение угла тангажа.

Рис. 13.9. Схема управления углом тангажа

Решая уравнение (13.18) совместно суравнениями (13.17), получим

p3 a1 p2 a2 p a3 b0 p a3 3 n0 p n32 f2 p n22 f2 ,

где

a1 2d0 0 nвk ;

a2 02 nв k n22k ;

a3 nвn22k ;

b0 nвk .

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

3 при слабом реагировании на возмущения f2 и f3 . Если передаточные числа k и k выбрать достаточно большими, то реакция системы на возмущения f2 и f3 будет слабой.

353

Выбор передаточных чисел k и k осуществляется в два этапа. Сначала выберем значение передаточного числа k из условия за-

данного переходного процесса во внутреннем контуре (см. рис. 13.9), передаточная функция для которого имеет вид

 

 

nв p n22

 

 

 

,

y

 

p2 2d p 2

где

2 2

n n k ;

2d 2d n k .

0

в 22

0 0 в

Затем выберем такое значение передаточного числа k , чтобы коэффициент затухания был оптимальным, например d 1. Находим

 

1

 

 

 

 

2d0ω0

 

2

 

 

 

 

 

 

k&

 

2d 2n

1

1

 

ω0

 

2d

 

ω

 

.

 

d 2n

d 2n2

 

 

 

nв

22

 

 

 

 

 

0

 

0

 

 

 

 

 

 

 

22

 

22

 

 

 

 

 

 

Для внешнего замкнутогоконтура(см. рис. 13.9) можно написать

 

 

nвk p n22

,

3

p3 A p2

A 2 p 3

 

1

2

 

где

 

 

nвk

 

3

 

 

A 2d;

A 1

 

;

n k

 

n .

 

1

2

2

 

в

22

 

 

 

 

 

 

354

Глава 14. ДИНАМИКА БЕСПИЛОТНОГО ЛЕТАТЕЛЬНОГО АППАРАТА ПРИ УЧЕТЕ УПРУГОСТИ

ЭЛЕМЕНТОВ КОНСТРУКЦИИ

14.1. Балочные модели изгибных колебаний корпуса дрона

Ранее в главах 12–13 предполагалось, что БПЛА является абсолютно жестким телом. Такое предположение справедливо до тех пор, пока собственные частоты замкнутой системы «автопилот– БПЛА» значительно ниже наименьшей частоты собственных изгибных колебаний корпуса БПЛА. При больших скоростных напорах воздуха могут возникать аэроупругие колебания. Аэроупругие колебания, возникающие в результате взаимодействия упругих

иаэродинамических сил, совершаются с частотами, близкими к собственным изгибным колебаниям конструкции (фюзеляжа, крыльев).

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

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

Рассмотрим влияние упругих колебаний БПЛА на его передаточную функцию в продольном движении.

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

икрыльев (рис. 14.1). На рисунке показаны φ1 x и φ2 x – упру-

гие линии фюзеляжа; M1, , M6 – аэродинамические моменты,

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

гать, что его упругие линии φ1 x , φ2 x , подобны упругим линиям балки со свободными концами.

355

Рис. 14.1. Схема действия сил и моментов на упругий летательный аппарат

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

 

2

 

 

2

y

 

 

2

y

F x y

 

 

 

EI x

 

 

m x

 

 

 

 

2

 

 

 

2

 

x

 

x

2

t

x

 

 

 

 

 

 

(14.1)

D x 2 y T t x xi .

x t

Граничные условия для балки со свободными концами сводятся к следующим граничным условиям:

2 y 0,t

 

3 y 0,t

0;

2 y l,t

 

3 y l,t

0.

(14.2)

x2

x3

x2

x3

 

 

 

 

 

В выражениях (14.1) и (14.2) приняты обозначения:

x– текущая координата вдоль продольной оси БПЛА; l – длина БПЛА;

t – время;

yx, t – прогиб в точке x ;

EI x – жесткость на изгиб;

mx – закон распределения массы по продольной оси;

F x – коэффициент пропорциональности восстанавливающей аэродинамической силы, направленной нормально к продольной оси;

356

D x – коэффициент аэродинамического демпфирования;

T t – сила, создаваемая управляющим органом нормально

кпродольной оси;

– угол отклонения управляющего органа (например, руля

высоты);

x xi – дельта-функция, определяемая соотношениями

 

x xi 0 при x xi ,

 

xi

 

x xi dx 1 при

 

 

 

x xi.

0

 

 

 

Для упрощения задачи предположим, что величины восстанавливающей и демпфирующей аэродинамических сил, действующих нормально к продольной оси, зависят только от вращательного эквивалентного жесткого БПЛА, имеющего одинаковую аэродинамическую форму и обладающего одинаковым распределением масс:

F x

y

F x

y

R ;

 

 

 

x

 

 

 

 

 

 

 

x

 

 

(14.3)

 

2 y

 

2 yR

 

 

D x

D x

,

 

 

t x

t x

 

 

 

 

 

 

 

где yR x,t – длина дуги, описываемой точкой недемпфированной

оси в данном плоском вращении, происходящем под действием аэродинамических сил.

Для вращательного движения эквивалентного жесткого тела имеем

yR x,t x x0 t ,

(14.4)

где x0 – координата центра масс;

– угол поворота тела около центра масс (этот угол в зависимости от направления соответствует углу тангажа или рыскания).

357

Используя выражение (14.4), вместо зависимостей (14.3) можем написать

F x

y

F x t ;

 

x

 

 

 

 

 

 

 

2 y

 

t

 

D x

D x

.

 

 

 

t x

 

t

 

 

 

 

Учитывая, что масса эквивалентного

жесткого тела равна

m x x x0 , а коэффициенты демпфирующей и восстанавливающей сил соответственно D x иF x , дифференциальное уравнение вращательного движения тела запишем в виде

m x x

x

 

d2θ

 

D x

F x θ TδΔ x x

.

dt2

 

 

 

 

 

 

0

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

i

Умножим обе части этого уравнения на x x0

и проинтегрируем

по x в пределах от 0 до l :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d2θ

lm x x x

2 dx l D x x x

dx

 

 

 

 

 

 

 

dt2 0

 

 

 

 

 

 

0

 

 

 

 

 

dt 0

 

 

0

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

θ F x x x0 dx Tδ l x0 .

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вводя обозначения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I

0

lm x x x

 

2 dx; d

0

 

1

l D

x x x

dx;

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I0 0

 

 

 

 

 

 

 

 

 

 

 

1

l F x

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

μ

 

 

x

x

dx, μ

 

 

T l x

,

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

1

 

I0

0

 

 

 

 

 

 

I0 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

можем написать

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d2θ

 

+ d

 

μ

 

θ = μ δ.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 dt

 

 

 

 

 

 

 

 

 

 

 

 

dt2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

1

 

 

 

 

358

Преобразовав уравнение по Лапласу, получим

 

θ p

μ1

δ p .

(14.5)

p2 d0 p μ0

Преобразуя уравнение (14.1) и граничные условия (14.2) по

Лапласу, найдем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d2

EI

x

 

dy x, p

 

p

2m x y x, p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

2

 

dx

2

 

 

(14.6)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Tδ p x l pD x F x p ;

 

 

 

 

 

 

d

2 y

0, p

 

d3 y 0, p

0;

 

 

 

 

 

 

 

 

 

 

 

dx2

 

 

 

 

dx3

 

 

 

(14.7)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d2 y l, p

 

 

d3 y l, p

 

 

 

 

 

 

 

 

 

0.

 

 

 

 

 

 

 

 

 

 

 

 

dx2

 

 

 

 

dx3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя выражение (14.5) в уравнение (14.6), получим

 

 

 

 

d

2

 

 

 

 

2

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

EI d

 

 

p2m

x y

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

dx

2

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(14.8)

 

 

 

 

 

 

 

 

 

 

 

μ1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T x x

 

 

 

 

 

 

 

 

 

 

 

pD x F x

δ

p .

 

 

 

2

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

p

d0 p μ0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для решения этого неоднородного уравнения, сначала рассмотрим решение однородного уравнения

d2

 

d2η

 

 

EI x

 

 

p2m x 0

dx2

dx2

 

 

 

 

 

 

 

 

при граничных условиях (14.7), в которых следует положить y .

359