Файл: Колесников К.С. Упругий летательный аппарат как объект автоматического управления.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