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

книги из ГПНТБ / Кашин Г.М. Автоматическое управление продольным движением упругого самолета

.pdf
Скачиваний:
9
Добавлен:
24.10.2023
Размер:
10.78 Mб
Скачать

где

\ \ j

« V / =

™ z i j ^

d f k d X j

d r k .

d x j

1

1-f A

1

1-j- A

1

1 +

A

/ * ( -

i

Гк

 

4

d s k

r k

d X j

d f l

 

d x j

 

d f i

l

!

i

 

^

 

d f

i

d z f

d r k

,

_

~?

d x j

'-‘k

 

1

d f 2

2 i

 

 

 

i

-

d X j

 

 

1

d f i ~

i 2 i ;

 

 

i

-ь / 2

d ^ j

z. =z, .

 

1

d f i

 

 

ii + f t

d z j

T ‘u

i

d s k

1

dRk \

 

d x j

Rk

дх; I

( n — A ) \ H

Л*Гу

dR.

+

1)M2] ^ ^*- +

Cbty

 

d X j

d/fe <^7- t^Os

dyj

dsk

+ 7 r ( ^ + sH ^ [ ? 2+ ( ^ - i ) M 2+ tg aXi]};

Rk

=/*(— % +

dst

1

. соэф,-

) :

 

s k difj

Rk dyj

 

 

rk dyj

 

 

tg Хг (<7 c o s —e* sin ф(.);

 

(2.

170)

гц

 

 

 

 

 

 

 

 

:2q cos ф, tg X/+

(P - z ’i) sin ф, tg ул -

 

 

dy j

 

- ч {

sin ф;

dZ:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dy j

 

 

,3

 

_ * . =

[=Р + (Л - 1)M2]

Sk _dsk,

 

гк

 

dR,

 

 

 

 

 

 

 

 

 

 

(J(/y

 

 

 

 

 

 

dy}

 

Rk

 

 

+

{k -

1)M2 + tg2 x/] cos <|>,. - e* sin ф,- tg x/} — -

 

- ^ - H + < i W + ( k - m 2+№xi}) t

;

 

 

 

 

 

 

 

 

 

dtfj

d f u

_

f

l

1

d r k ,

1

^

1_

! sin^,-\ .

dzj

 

 

 

rk

dzj

sk

d zj

Rk

d zj

Я

d r k

_

-

t g

l i

(4 sLn

+

4 cos ф,);

 

 

dzj

=

JU ^

 

 

 

 

 

rh

 

 

 

 

 

 

ds

:= -

2<7 sin ф, tg X i +

 

 

C O S ф, tg x, -

( / > —

« ; )

dz

dz,

E* (cos <|»,

dz.

120

 

 

 

 

 

 

 

 

 

 

I

 

 

^ -

= [? 4 (* -

l)Ma]i5-

.Ё£*

Rk

+

l

 

 

dzj

 

 

 

Rk

dzj

 

 

 

 

 

■'r(k— l)M 2 + tg2x;] sin^.-f

cos 6,. tg y;

 

 

 

 

\zl-2r ^ [ f + {k~

1)M2 +

 

. (2.170)

 

 

 

 

 

 

 

 

 

+*■*]>

 

dzj

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

(k=\,2).

 

 

 

 

 

 

 

 

 

 

)

Пределы интегрирования в (2. 168)

и (2. 170)

имеют следую­

щие значения:

 

 

 

 

 

 

 

'li

Z2i

 

-*

—*

 

" i Zij,

Z2j

 

 

 

П р и

Zij,

Z 2]‘

 

” !

 

 

 

 

 

 

 

 

2bd

 

26a

 

 

2n-

 

S

Z2i= z\]

при

—- i L ^ z l j ^ d L

 

26

 

 

 

2b0

2b0

 

26n

Z n =

 

 

Z J | = A

при г ^ <

h

-* \

//

 

260

— , z2/ > — ;

 

 

 

260

 

 

260

J ^ 2 6 0

 

Z \ l =

Z*i],

Z 2i = Z * 2j

При

— - ^ - O iy ,

2b0

 

CM

 

-*

 

 

 

 

2,o0

 

l:

Z \ l =

 

//

ПРИ

 

/; ^ -*

, /;

-* -

 

Z \ } ,

*2/= —

 

26q

 

, 22 ; > - f ;

 

 

 

 

26о

 

 

26q

26q

 

z*/ =

/? +

£2 sm X,

 

 

T

\2

 

 

 

 

 

M

— S2c0s X ;- ? “ cos X/;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(£=1,2),

 

 

 

где k — размах вихря г вдоль оси Oz'.

Если подкоренное выражение в (2. 171) является отрицатель­ ным или равно нулю, то в этом случае составляющие скоростей в (2. 169) следует полагать равными нулю. Если пределы интег­

рирования zu и 22/ в (2. 171) являются постоянными величинами, не зависящими от Xj, уj и zj, то их производные в (2.170) сле­ дует полагать равными нулю. В противном случае эти производ­ ные определяются выражениями:

д*ы

дгк)

дгы

дгМ

д*ы

дгN

 

 

дXj

дх

дуj

dyj

0*1

дг1'

 

сч

 

s i n x / ~ ( - l ) ft —

е 2

C O S 2 XI

COS X;,

см

дх 4

 

 

 

 

] /

( i

^2 COS

2 V,— дг2

 

 

 

 

 

 

 

 

г ъ

 

121

Рис. 2.28. Пересечение характеристического конуса с вершиной в точке х, у, г с плоскостью у = 0

дгк)

 

, so sin

 

I

cos /i sin'fy + ( - l

x; cos l i sin O'; — q cos yf

ду,

 

 

X

 

 

 

гп2cos2у.,- — q2

 

 

 

 

 

X COS/j,

 

 

CSI

 

 

 

 

d4j

 

! y* £2 sin

CM

cos /j cos 6r

7,- COS 7,; COS Цj -1- q S in

 

 

 

 

 

 

I

 

_ T _ \ 2

 

 

 

( ТГ)" ~ £2 c°s2 7-‘ __

 

 

 

 

M

 

X cosZ(.,

 

 

(£= 1,2).

Таким образом, имея приведенные выражения в координатах, связанных с вихрем г, легко переписать их в основной системе координат, связанной с «рылом, и получить компоненты скоростей в точке / от всей вихревой системы

крыла.

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

щенных . скоростей ср, построить

' функцию <р(/, т) по заданному направлению / И продифференцировать ее. Результат дифференцирования и даст искомое значение скоса потока по этому направлению. Таким обра­ зом, задача сводится к отысканию

ф(т) в точке х, у, z пространства, проекция которой на плоскость xOz показана на рис. 2.28. Аналогичные задачи при мгновенных смещениях «рыла, постепенном его входе в резкоограниченный порыв и при отклонениях руля для точки, лежащей в плоскости xOz, были рассмотрены выше. В данном случае метод решения по существу ничем не отличается от рассмотренных, только в некоторые формулы необходимо ввести координату у. В частно­ сти, формула (2.163) для расчета безразмерного времени Xi и t 2 будет иметь вид

хг= х

mih

т i + п ]

/

{т — i) (п j) + у?/№

 

2b0k

2

' +

 

М

— ДТ,

 

 

 

 

 

( г =

1,

2)

 

122

где у ~ ~ , а / в зависимости от вида задачи является либо раз­

махом крыла, либо размахом руля.

 

 

 

Входящий во все приведенные выше формулы

параметр ф,

в данном случае записывается следующим образом:

 

ri —

— ^i) —

(2. 173)

В соответствии с (2. 173) интегралы

(2.154) принимают вид

J\mn(т) —4Л(т) \Y [т i-\- 1 ){п — у-)- 1)— y2//z2 —

 

V(m — i)(n — y-f 1)— f/2/A2-j-

 

+ -----, l

= ( a r c t g l / ,(/!-»

+ В ( " - У+

1)

— arctg ]|/ /

^ / Л2

/

 

[V(m — i-\- !)(«. — у) — г/2//г2—

 

i){n — j) — уг1кг)

H e

] / &

 

+ 1)(П — j)

1 —

 

y2i№

J2mn{l) = mhJ\mn(*) — h2(x) ( К (m—/ + 1) (n—y+1)—г/2/Л2)3

«5

n — j + 1

_ ( l / ( m — t) (n — у +

\ 3

1) — i/2/ft2)3

Л — У+ 1

+ — ^2(t) ( ^ ( w —г + 1) (Я—у)—у2//г2)3

( У (m—г) (и—у)—(/2/Л2) 3

п — у

И— У

(У(~я-у + 1) (<я— 1+ 1) - ^ 2/Й2):

Jidnn(X)= nhJ\ H X)----Г А2(Т)

т — / + 1

(У (п —у) — i + 1)

з

г/2//г2)3

т i + 1

 

( К(/г — у + 1) (от — г) — у2//г2)3

/ К —

/

( у" (п — у) (от — о г/2//г2)3

/гг — г

123

Все остальные формулы остаются без изменения. Имея в виду изложенное, легко получить аэродинамические характеристики системы «механизированное крыло — механизированный стаби­ лизатор» при произвольных во времени смещениях ее точек, а также при действии порыва произвольного профиля. По ана­ логии с рассмотренным выше установившемся движением можно также представить себе самолет, схематизированный системой тонких пластин (см. рис. 2. 14). Тогда, используя приведенные формулы, можно получить приближенные выражения для опреде­ ления аэродинамических характеристик системы «механизиро­ ванное крыло — корпус— механизированный стабилизатор».

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

Г л а в а

III

УПРУГИЙ САМОЛЕТ КАК ОБЪЕКТ АВТОМАТИЧЕСКОГО УПРАВЛЕНИЯ

Система автоматического управления самолета включает контуры управления угловым короткопериодическим движением и контуры управления движением центра масс.

Ниже будут рассмотрены короткопериодическое движение самолета и система управления этим движением.

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

При разработке САУ, выполняющей указанные функции, приходится решать ряд проблем, основные из которых следую­ щие:

выбор состава и места установки датчиков;

синтез корректирующих устройств и фильтров, обеспечи­ вающих необходимое качество регулирования;

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

разработка принципов построения самонастраивающихся контуров САУ упругого самолета;

разработка способов фильтрации помех САУ от упругих колебаний и т. д.

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

кпродольному движению сверхзвукового тяжелого самолета

(СТС) и сверхзвукового легкого самолета (СЛС).

В этой главе основное внимание будет уделено динамической модели упругого самолета, анализу характеристик сигналов дат­

чиков в точках их измерения, а также

рассмотрению возмуще­

ний, действующих на самолет в полете.

 

3 . 1 . УРАВНЕНИ Я ДВИЖ ЕНИ Я

УПРУГОГО САМ ОЛЕТА

В первой главе подробно были рассмотрены уравнения коле­ баний упругого самолета при различных видах схематизации

125

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

любых

видов

продольного движения летательного аппарата

и при действий на него воздуш­

Г

ного порыва как

в дозвуковом,

так и

сверхзвуковом полетах.

 

Таким

образом,

подставляя

 

аэродинамические

характери­

 

стики, полученные по форму­

 

лам второй главы, в уравнения,

 

приведенные в

первой главе,

 

получим в общем случае систе­ мы дифференциальных уравне-

Рис.

3. 1. Балансировоч­

Рис. 3.2. Характеристики управляе­

ные

кривые тяжелого са­

мости тяжелого самолета:

 

молета:

-------------ж е с т к и й с а м о л е т ;

ж е с т к и й

с а м о л е т ;

— • — у п р у г и й с а м о л е т ;

л е т н ы е и с п ы т а н и я

— • — у п р у г и й

с а м о л е т ;

 

------------л е т н ы е и с п ы т а н и я

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

В зависимости от характера исследования по материалам первой и второй глав можно получить системы уравнений как для расчета характеристик самолета с учетом статической аэро­ упругости, так и для расчета его характеристик с учетом дина­ мической аэроупругости. Первая система представляет собой хорошо известные уравнения продольного движения жесткого самолета [28], но с аэродинамическими коэффициентами, рассчи­ танными с учетом статической аэроупругости по формулам раз­ делов 1.2 и 2.2. Как показывают результаты расчетов (рис. 3. 1 и 3.2), эта система может с успехом применяться для расчета балансировочных характеристик и маневров упругого самолета при полете в спокойной атмосфере. Вторая система представляет

126

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

*уравнениям (1.61) с применением аэродинамических характе­ ристик из раздела 2.3 (рис. 3.3). Если же определяются и рас-

Ми тм

Рис. 3. 3. Амплитудно-фазовые частотные

Рис. 3.4. Динамическое

нагру­

характеристики

тяжелого

самолета при

жение

фюзеляжа

по

методу

действии вертикального порыва:

сосредоточенных

масс;

------------ р а с ч е т ;

-------------л е т н ы е

и с п ы т а н и я

--------- р а с ч е т ;

-------------- э к с п е р и м е н т

пределенные

по конструкции характеристики

самолета,

то

еле-

дует пользоваться формулами для расчета динамических дефор­ маций или по методу форм собственных свободных колебаний или по методу сосредоточенных масс. При надлежащим образом выбранном числе тонов упругих колебаний и числе сосредото­ ченных масс оба метода дают практически одинаковые резуль­ таты. Сходимость с экспериментальными данными также удов­ летворительная [8] (рис. 3. 4).

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

127

дальнейшего анализа и синтеза САУ будем использовать аэро­ динамические характеристики в квазиустановившемся движении без учета нестационарное™ обтекания. Заметим здесь, что при учете нестационарное™ обтекания общий подход при исследова­ нии САУ сохраняется.

Как было показано в разд. 1.3, уравнения для исследования динамики системы с дискретно распределенными параметрами по методу форм имеют вид

t

 

(i = 1, 2, . . . , п'к

'(1.67)

где qi — обобщенная координата г'-го тона колебаний; toг — собственная частота колебаний г-ro тона;

Напомним, что обобщенная масса Mi и обобщенная сила Qt определяются соответственно по формулам:

и =1 и=1

где

Ри — сосредоточенная сила, приложенная в и-й точке; т ц — и-я сосредоточенная масса системы;

riu — смещение м-ой массы по i-й форме колебаний; k — число сосредоточенных масс.

Сосредоточенная сила Ри представляет собой аэродинамиче­ скую нагрузку на участке отсека конструкции, включающего сосредоточенную массу. Эта нагрузка рассчитывается по данным разд. 2. 2 в виде

1 . р=1

где 6р, бр—соответственно угол и скорость отклонения р-го руля; д — число рулей;

WE— вертикальная скорость воздушного порыва.

Производные Ру

в

(3. 1)

представляют собой интегралы по

площадям

отсеков от

соответствующих

производных давления

АР. Имея это в виду и подставляя (3. 1)

в уравнения движения,

получим

 

 

 

 

 

4i — 2

\djflj +

bjiQil + 2

\пЪрЧ^р+ ni

qbp\ - г nw qW в , ( 3 . 2 )

7 = 1

 

Р = 1

 

Р

 

 

 

{ i ~ U

2

 

128

г д е

— U \ Р чЦ х , z )r i (x, z) dxdz\

 

М; 1.1

l

c c

APq‘(x, z )r t(x, zjdxdz,

:—

\\

M;

.).)

(3. 3)

 

S*

 

 

APbp (x, z)ri(x, z)d^dz',

Здесь 5* — площадь интегрирования по поверхности самолета;

v— скорость полета.

В(3. 2) предполагается мгновенный охват самолета порывом WB(t). Для учета конструкционного демпфирования к коэффи­ циентам dji необходимо добавить, например, по формуле (1.83)

члены вида

Ada=

-----Loy

при г = 3, 4 ,..., п. Здесь

6 ,— лога­

 

 

я

 

 

рифмический декремент затухания.

ином виде,

Уравнения

(3.2)

часто

записывают в несколько

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

g1 = H = v ( b — а ) , <7а = » ,

(3.4)

где и — путевой угол атаки корневого сечения поверхности;

4— угол тангажа;

Я— смещение самолета по высоте.

Подставляя (3.4) в (3. 2), окончательно получим

П

а — dn a - ф d2 1 &

- ) -

^

[d^qj - ф

bjXq; ] - j -

 

 

 

 

П

 

 

 

d 12a - ф

^ 2

2 ^

+ 2

j

^ j

]

 

 

 

 

 

nwbWB',

 

(3. 5)

 

 

 

 

 

 

 

 

 

 

 

П

 

 

 

qt

d ria - ф

d 2ib - ф

^

\djtqj-\~bjtqj\ - ф

-ф Ягэ?Д 4- «8к?(.8к-ф tlwqW B\

(i = 3, 4, . . . , n). I

5

3819

129

Соседние файлы в папке книги из ГПНТБ