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

книги / Системы управления летательными аппаратами и их силовыми установками

..pdf
Скачиваний:
8
Добавлен:
12.11.2023
Размер:
7.92 Mб
Скачать

Отсюда можно найти угол атаки, обеспечивающий реализацию рикошетирующего режима полета:

а >

2G

(7.11)

 

К*Рv X

7.3.Алгоритмы реализации требуемого угла атаки

Взаключение рассмотрим принципы реализации расчетного значения угла атаки а р.

Как известно, непосредственно регулировать угол атаки невоз­ можно. Управлять можно только углом тангажа. Поэтому алгоритм управления углом атаки можно представить в следующем виде:

1.Измерение составляющих вектора скорости vxi,vyi.

2.Определение угла наклона вектора скорости 0, :

0,- = arctg = — .

(7.12)

3. Определение угла атаки ЛА

 

а , = К \ ~ 6,-

(7.13)

4. Определение расчетного значения угла атаки

 

 

 

(7.14)

либо

 

 

%

= С-

(7.15)

5. Определение угла тангажа

 

 

0,-=

( а , - а р/).

(7.16)

Глава 8 ОПРЕДЕЛЕНИЕ ПАРАМЕТРОВ ДВИЖЕНИЯ АБЛА

8.1. Определение вертикальной составляющей действительной скорости и высоты полета АБЛА

Как видно из вышеизложенного, управление полетом АБЛА в вертикальной плоскости осуществляется с помощью регулирования угла атаки, для определения которого необходимо наличие информа­ ции о действительных параметрах движения АБЛА, в частности со­ ставляющих вектора действительной скорости vx,v y.

При определении горизонтальной скорости vA. проблемы не воз­

никает, так как

 

vx = wx,

(8.1)

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

Значительно сложнее реша­ ется задача по определению вер­ тикальной составляющей вектора действительной скорости. Рас­ смотрим данный вопрос подробнее (рис. 8.1).

Представим зависимость для вертикальной составляющей ско­ рости:

j( Wy-g)(lt

(8.2)

Ускорение силы тяжести в полете можно вычислять с помощью известных зависимостей, например, с помощью зависимости (3.21). Для упрощения процедуры исследования при решении рассматри­ ваемой задачи используем более простую зависимость для опредоления g :

g = g<>

J l

(8.3)

R3+H y

 

здесь g0 - ускорение силы тяжести на поверхности Земли. Предста­ вим зависимость для вычисления высоты полета:

Il = Я 0 +vy0t + '\\(wv - g)clt2

(8.4)

0 0

 

Будем считать, что Я 0 = О и vy0 = 0, тогда

 

H = ])(wy - g ) d t2

(8.5)

00

 

Анализируя зависимости (8.3) и (8.5), можно

отметить, что

g = / (Я), а Я = f ( g ) , т.е. в системе определения ускорения силы тяжести имеется положительная обратная связь. Поэтому погрешно­ сти вычисления g (методические и инструментальные) с течением времени будут возрастать, что приведет к возрастанию ошибок при вычислении параметров вертикального канала: действительной ско­ рости и высоты полета. Все это обусловит возрастание погрешностей при вычислении значений угла атаки, а значит, приведет к ошибкам при управлении полетом АБЛА.

В заключение можно отметить, что у баллистических ЛА дли­ тельность управления полетом на активном участке траектории не­ велика (5-6 мин), поэтому погрешности вычисления параметров вер­ тикального канала не достигают недопустимо больших значений. Однако управление движением аэробаллистических ЛА осуществля­ ется в течение всего полета (30-40 мин), поэтому погрешности вы­ числения параметров вертикального канала могуг существенно пре­ высить допустимые пределы.

8.2. Определение ускорения силы тяжести при различных режимах полета АБЛА

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

I. Режим планирования. При реализации режима планирования выполняется условие v —0. Представим зависимость для верти­ кального ускорения

= vy + g-

(8-6)

Учитывая приведенное условие, получим

 

£ = |й'г|-

(8-7)

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

связь, благодаря чему не будет происходить накопление ошибок при определении g.

2. Режим рикошетирования.

Рассмотрим методику определения ускорения силы тяжести при реали­ зации режима рикошетирования (рис. 8.2).

Представим зависимость для определения вертикальной сос­ тавляющей кажущейся скорости АБЛА

wy = vy + \ë dl

(8.8)

0

 

В уравнении (8.8) два неизвестных: v,, и g. Если исключить vy,

то можно будег вычислить значение g. Рассмотрим данный вопрос подробнее.

Запишем зависимость для высоты полета АБЛА (см. рис. 8.2)

 

H - Н() - vw,t -I АН sinco,, /

(8.9)

где Я 0 - начальная высота полета;

 

v

- постоянная составляющая вертикальной скорости;

 

АН - амплитуда переменной составляющей высоты полета;

()}р

частота рикошета.

 

Продифференцируем данное уравнение, получим зависимость

для определения vy :

 

 

vy = -v ^-i-AZ/cOpcos a y .

(8.10)

Анализируя (8.12), можно отметить, что в момент времени, ко­

гда cos со ,/ =0,

 

 

vy = -vy0.

(8.11)

Задача состоит в нахождении некоторой математической функ­ ции, с помощью которой можно будет определить данный момент

времени.

 

Из выражения (8.10) получим зависимости для

vy и v,,.

vy = -АЯсОр sin со,,/,

(8.12)

v,, = -ДЯсОр cos со,,/

(8.13)

Примем во внимание, что в выражении для кажущегося уско­

рения

 

ùy^V y+ g,

(8.14)

£ = £o + £r*l A gsm ay,

(8.15)

где g0 - начальное значение ускорения силы тяжести;

gc постоянная составляющая скорости изменения ускорения

силы тяжести;

 

Ag

- амплитуда переменной составляющей ускорения силы

тяжести.

 

 

Из выражения (8.14) получим

 

 

% = v y+g,

(8.16)

где

 

 

 

g ^ g c + Ag<ùpcos(ùpt.

(8.17)

Подставим в выражение (8.16) зависимости (8.13) и (8.17).

Тогда

 

 

 

wy = cûp (- АЯсОр+ Ag)cosoy +gr .

(8.18)

Необходимо учесть, что ускорение силы тяжести с течением

времени

изменяется незначительно, поэтому можно

считать, что

gr = 0. При учете данного допущения зависимость (8.18) примет вид

wy = а>р(~ Д/7(0р + Ag)coscop*.

(8.19)

Как видно из выражения (8.19), функция wy = /( /)

обращается

в нуль в момент времени, когда coscop =0. В этом случае, как указы­

валось выше, выполняется условие (8.11).

Необходимо отметить, что в отличие от функции (8.10), кото ­ рую невозможно реализовать на борту ЛА ввиду отсутствия инфор­

мации о

vv0, АН, cùp, функция (8.19) может быть достаточно просто

реализована. Действительно

 

 

Wj+i-Wj

( 8.20)

 

At

 

 

где At -

шаг интегрирования (/ = 1,2,3...).

 

Информация о кажущемся ускорении (й^) поступает с акселе­ рометров. Таким образом, на борту может быть сравнительно легко определен момент времени /, (/ = 1,2,3...), когда функция wy обра­

щается в нуль.

Используя выражение (8.8), представим зависимости для верти­ кальной составляющей кажущейся скорости АБЛА для /-го и (/ + 1)-го моментов времени (/ = 1,2,3,...) :

U

w v i = ~ v y c +

\ g d t ,

(8.21)

О

 

 

 

*1+1

(8.22)

w y ( M ) = - v y c +

\ s d t -

 

о

 

 

Вычтем (8.21) из (8.22). Полагая, что за время, равное

/(+1

,

g - const, тогда получим

 

 

 

WAM) ~ ия = S (*M “ О •

(8.23)

Из выражения (8.23) найдем зависимость для вычисления g:

 

W y i M )

^ yi

(8.24)

Si

 

*м - h

здесь vv;,w(+, - значения кажущейся скорости АБЛА в момент вре­ мени, когда wv = 0.

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

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

Глава 9 ПРИНЦИПЫ УПРАВЛЕНИЯ ПОЛЕТОМ АВЛА

ПРИ ОРГАНИЗАЦИИ МАНЕВРА

9.1.Общие сведении

оманевре аэробаллистического ЛА

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

Чаще всего траектория движения маневрирующего АБЛА пред­ ставляет собой колебания относительно некоторой базовой траекто­ рии, соединяющей по кратчайшему пути центр масс АБЛА с точкой приземления. Колебания АБЛА при маневре в общем случае можно описать с помощью следующей математической зависимости:

2Н(0 = 2(Г + ГИ),

(9.1)

здесь zMкоордината маневра;

Ти - период маневра.

На первом этапе рассмотрим принципы организации бокового маневра, т.е. маневра - в плоскости местного горизонта. Данный ма­ невр АБЛА осуществляется за счет боковой составляющей аэроди­ намической силы, которая появляется при наличии угла крена АБЛА.

9.2. Выбор закона управлении углом крена АБЛА при осуществлении бокового маневра

Представим уравнение движения АБЛА по оси г в нормальной системе координат при учете основной силы, действующей в напр.т лении, перпендикулярном вертикальной плоскости (рис. 9.1):

m z ~ R z,

(9.2)

где R. боковая составляющая аэродинамической силы,

Rt = R s m y M,

(9.3)

где фм - угол крена. Будем считать, что при осуществлении маневра

угол крена не превысит 30°.

 

На основании данного

поло­

жения можно принять допущение

Rz =Яфм-

(9.4)

Подставим (9.4) в (9.2), по­

лучим

 

mz = R<pM.

(9.5)

Для организации колебатель­ ного движения относительно оси z необходимо создать консерватив­

ную систему управления маневром. С этой целью примем закон управления углом крена АБЛЛ в виде

Фм = ~K zz-

(9-6)

При учете (9.6) уравнение движения АБЛА преобразуется сле­ дующим образом:

mz + KzRz = 0.

(9.7)

Перейдем в область оператора Лапласа и получим характери­ стическое уравнение системы

р * + Е Л = о.

(9.8)

т

 

Найдем корни данного уравнения

p , , 2 = ± j M ^ .

(9.9)

V т

Тогда решение уравнения движения АБЛА можно представить в виде

z = y4Msin<oM/.

(9.Ю)

Анализ зависимости (9.9) показывает, что движение АБЛА име­ ет колебательный характер.

Найдем выражения для частоты и амплитуды маневра (сом, Аы). Частота маневра

(9.11)

Для получения зависимости для амплитуды маневра необходи­ мо учесть, что при колебательном движении АБЛА относительно оси z, угловое движение АБЛА по крену также представляет собой гар­ монические колебания

Фм = Фм„т sincùj,

(9.12)

где <рМ1гах - максимальное значение угла крена.

 

Подставив (9.12) в (9.5), получим

 

g _ ^ Ф м шах sino>M/.

(9.13)

т

Проинтегрируем данное выражение дважды и при учете начальных условий (z0 = 0, z„ = О) найдем зависимость для координаты маневра:

z

* Ф м Пгах sincoMf.

(9.14)

 

 

rrml

 

Тогда можно записать, что

 

 

 

 

А

*|

 

—*Фмшах

(9.15)

 

 

2

 

 

 

/И(0М

 

или при учете (9.11)

 

 

 

 

АМ

Фм

(9.16)

к .