книги из ГПНТБ / Кашин Г.М. Автоматическое управление продольным движением упругого самолета
.pdfгде
\ \ 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
дгк) |
|
, 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 |