Файл: Колесников К.С. Упругий летательный аппарат как объект автоматического управления.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 27.06.2024
Просмотров: 153
Скачиваний: 1
3 / Л £ /
u «
d^(x,t) \ . . . |
дЧ(х, t) |
|
a y * - ) + mW |
0 ; |
|
dt2 |
||
tni\i И- m1Л1~b ^lfii — 0; |
||
d2l{x,t) |
|
d2Ux,t) |
dx2 |
■} = °, EJ---- = 0 |
|
|
dx2 |
при x = 0;
* ( £ / — Ь г Ч = - |
|
« - й Н - = » |
||||||
dx2 |
|
|
|
|
|
dx2 |
||
|
|
|
|
|
|
|
|
при x = /; |
|
|
|
|
|
|
|
|
(4.36) |
U h - 0 ) = |
|
l(h + |
|
0), |
|
d£(*,Q |
||
|
|
|
|
|
|
|
dx |
Л -0 |
dl(x,t) |
|
, |
EJ |
d2l{x,t) |
||||
dx |
|
|
л+о |
|
, |
|
||
|
|
|
|
|
dx2 |
/1-0 |
||
|
|
|
£ / |
d2j(x,t) |
|
|||
|
|
|
|
dx2 |
1h+0 |
|||
|
|
|
|
|
||||
L |
( |
' |
e i * |
|
M |
|
h- 0 |
|
dx |
|
|
dx2 |
|
||||
± ( EJ |
|
|
dx2 |
|
). |
|
= |
"M £ 1+ 111). |
d x ' |
|
|
|
Л+0 |
|
|||
В точке с координатой х = хд.у находится датчик угло |
||||||||
вого положения аппарата, |
входной сигнал которого |
|||||||
|
|
|
М О = |
— |
|
дЦх, |
t) |
|
|
|
|
|
dx |
XД.у |
|||
|
|
|
|
|
|
|
Рассмотрим вынужденные колебания упругого аппарата, если угол б изменяется по закону
б = б 0е/“ *.
Реакцию конструкции на входной сигнал можно предста вить в виде
£ (* , 0 = £ {х) е/и (, т]1 ( 0 — тц е /“ *.
196
Уравнения и граничные условия (4.36) преобразуются к выражениям, аналогичным (2.46), за исключением гра ничного условия
d |
d2Ux) |
= — Pybo- |
|
dx |
dx2 |
||
|
Общий вид функции g(x), характеризующий деформации конструкции при' вынужденных колебаниях, согласно (2.47) — (2.53) представим в виде
I (х) = С щ и ( х ) + С 2« 1 2 ( х ) .
Удовлетворив граничным условиям при х = 1, получим следующие уравнения для определения постоянных С\ и С2:
C1U31(/)-(- CzUiz(l) = 0;
Cl«41 (/) + |
С2«42(0 |
- Рубо. |
|
Находим |
|
|
|
Ci = Ру«зг(0 |
бо, С2 = |
Ру«31 (0 |
|
А(со) |
|
А (со) |
|
где |
«31 (0 |
«32 (0 |
|
А(со) = |
|||
«41(0 |
«42(0 |
||
|
Как следует из метода решения, величины |(х), гц, так же как «з1(/ ), и41 (/) и т. д., являются функциями частоты вынужденных колебаний, поэтому и определитель А (со) есть функция со. Равенство А(со)=0 определяет частоты собственных колебаний системы. Частотная характери стика летательного аппарата будет определяться из вы ражения
д%(х)
(4.37)
dx Xя.У
Поскольку демпфирование конструкции не учитывалось, то сдвиг фаз может быть равным нулю или я, что опре деляется знаком выражения (4.37).
Так как
dU*) |
--- Ci«21 (Хд .у)-)- С2«22(-^д.у) , |
|
dx |
||
ду |
197
т о
Щ (/“ ) = |
ЯуАо((р) |
|||
А (со) |
||||
|
||||
где |
|
«22(-^д.у) |
||
«21 (-^Д.у) |
|
|||
Ао(со) — |
|
«32 ( 0 |
||
«31(0 |
|
|||
Определитель Ао(со) обращается |
в |
нуль на тех частотах |
||
колебаний, которые соответствуют |
нулям передаточной |
функции W'y(s), расположенным на мнимой оси плоско сти 5. Покажем, что задачу о нахождении нулей пере даточной функции U?y(s), расположенных на мнимой оси, можно сформулировать как задачу о нахождении собст венных частот некоторой эквивалентной системы.
Рассмотрим краевую задачу, определяемую следую щим дифференциальным уравнением и граничными ус
ловиями, |
|
|
1 |
|
|
(Р I |
d2t{x) |
\ |
|
d_ |
E J d \ W |
= 0, |
EJ d\ (x) |
|
dx |
d,x2 |
|
dx2 |
|
|
|
|
при x |
0; |
|
|
|
|
i" (4.38) |
|
EJ |
dx2 |
= 0 при x = l, |
|
|
|
|
|
dl(x)
--------- —■0 цри x — Хду. dx
Остальные граничные условия в точке крепления упруго подвешенных масс имеют вид, аналогичный (4.36).
Общее решение этой задачи можно искать в виде
|(х) = СцМц (х) -р C2«i2(x).
Постоянные С\ и С2 определяются из следующей системы однородных уравнений:
C\U%i{X^.y) -р С2Ц22(Хд.у) — 0;
Ci«31 (/) -р С2«зг(/) = 0.
198
Нетривиальные решения возможны только при таких и, когда
| « 2 i(* A .y ) |
Ы22(* д .у ) |
| __ |
i Мз1 (/) |
Из2(/) |
I |
что совпадает с условием А0(со)=0.
Таким образом, нули передаточной функции W'y(s), расположенные на мнимой оси, совпадают с собственны ми частотами краевой задачи (4.38).
Используя метод прогонки, можно определять также нули передаточной функции, расположенные на действи тельной оси комплексной плоскости, если в уравнении (4.38) полагать величину ю2<0.
4.4. ПЕРЕДАТОЧНЫЕ ФУНКЦИИ УПРУГОЙ РАКЕТЫ С УЧЕТОМ КОЛЕБАНИЙ ЖИДКОСТИ И ИХ СВОЙСТВА
Общий вид передаточных функций упругого летатель ного аппарата представлен соотношениями (4.23), (4.32), (4.33). Распределение нулей и полюсов передаточной функции на комплексной плоскости 5 полностью опреде ляет свойства упругого летательного аппарата как объек та регулирования. Разумеется, в общем случае ввиду сложности задачи, значения нулей и полюсов и их зави симость от конструктивных параметров летательного аппарата могут быть определены только с помощью вы числительных машин. Однако некоторые качественные зависимости параметров передаточных функций могут быть получены аналитически для упрощенной постанов ки задачи. Эти упрощения имеют следующий характер: во-первых, из всего многообразия степеней свободы, ха рактеризующих возмущенное движение в целом, выбира ются две-три степени свободы , определяющие тот или иной вид возмущенного движения; во-вторых, в уравне ниях движения такой упрощенной системы опускаются второстепенные члены. Здесь приводятся упрощенные пе редаточные функции упругого летательного аппарата, которые позволяют более конкретно представить зависи мость нулей и полюсов от конструктивных параметров.
Вначале рассмотрим передаточную функцию ракеты как твердого тела. Уравнения возмущенного движения, если пренебречь инерционными и демпфирующими со-
199
с т а в л я ю щ и м и у п р а в л я ю щ и х о р г а н о в , и м е ю т в и д
0-1—iZ + d-i-iz -(- (d-io -f- е_ю) ф +
+ |
(fr-io + § -io ) Ф = — § -ie 6 ; |
(4.39), |
||
do—iZ -f- #-ooi|) -(- ( doo |
Coo) ф “b ^ооф — — §об6. |
|||
Здесь принято |
q~i(t) = z{t) = z — перемещение центра |
|||
тяжести ракеты, |
a qo(t) = ф(/) = ф — угол |
поворота ра |
||
кеты. |
сигнал |
датчика углового |
положения |
|
Входной |
Ру(О = Ф, следовательно, передаточная функция
Wy.T(s) |
Ф($) |
|
6(s) |
||
|
Если пренебречь движением центра тяжести, то из вто рого уравнения (4.39) получим
______ §06____________
Wy.T ( s )
OooS2 + {doo + Cqo) s + doo
Коэффициенты уравнений (4.39) имеют следующий физи ческий смысл: a_i_i = m — масса ракеты; a00 = J — мо мент инерции относительно оси Oz, проходящей через центр тяжести:
§ —16 = |
Ру, §06 = Ру (-^Д -Кц.т) > |
где Хд — координата точки приложения управляющей силы. Передаточная функция теперь может быть пред ставлена в виде
1Ку.т( s ) — |
Ру { хд |
*ц.т) |
|
1 |
|
|
/ |
|
s2+ |
(4.40) |
|
|
|
|
2h*Ts + ю2 |
||
|
|
|
|
|
*Т |
где |
|
|
|
|
|
п и |
___ |
Д)0 + |
С00 |
2 |
&00 |
Z/7*t |
— |
~ |
) |
С0*т |
— - ■ |
Для динамически неустойчивого аппарата со*т2<0, а для устойчивого о)*т2>0; обычно величина /г*т> 0. Рас пределение полюсов передаточной функции (4.40) и со-
200