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

книги из ГПНТБ / Павловский М.А. Влияние погрешностей изготовления и сборки гироприборов на их точность

.pdf
Скачиваний:
10
Добавлен:
23.10.2023
Размер:
6.09 Mб
Скачать

от односторонней функции

ßi

 

 

1ГГ = і І г -

» ( t - % - ) +

ß 1 ( T ) o ( T - т о - ) .

 

Воспользовавшись известным свойством дельта-функции

б (т — т 0 - )

 

 

 

 

Рх (т) 6 (т -

то-) -

ßi (то-) Ô (т - то-),

(ІѴ.55)

d p ;

запишем в

виде

 

выражение для - ^ -

 

Вторая производная с учетом формулы (IV.55) будет равна

dx2

I

dpi

Ô(T —То-)

 

 

dô (т — тп—)

+

ß x

( т 0 - ) -

dx

 

~*~

rfx

 

х = т 0 .

 

 

 

 

 

 

d2 ßl

1,

dp*

1

.

1 \

Умножив

теперь - ^ - н а

 

~ и

ßi н а

(1 — ^ г ) ,

сложив полученные выражения и сравнив найденный ре­

зультат с уравнением (IV.54), будем

иметь

R(t, то-) =

Рх(то-) +

1

; б (т — то-) +

 

 

т=т 0 _

+

ßl C*>-) d ô ( x

— х 0 - )

(IV.56)

 

 

dx

'

В дальнейшем для краткости будем писать т 0 - = т0 , ßi =

=h и / (т) 1 (т + ^о-) = / (т).

Решение уравнения (IV.54) запишем в виде

^ ( т ) = С і Ф і ( т ) + С 2 ф 2 (т),

(ІѴ.57)

где ф 1 і 2 — линейно независимые частные решения однород­ ного уравнения; С 1 ( С 2 — величины, которые должны быть определены так, чтобы при подстановке выражения (IV.57) в (IV.54) получилось тождество. Величины Сх и С 2 найдем, используя метод вариации произвольных постоянных. Вы­ ражения для них должны удовлетворять системе уравнений

С і ф 1 + С 2 ф 2 = 0, С^ф, + С 2 ф 2 = F ( T ) .

141

Отсюда

где

 

 

А(т) .

Фі

Фг

 

 

 

 

 

 

фі

ф 2

 

 

 

 

 

 

 

 

 

 

 

является определителем Вронского уравнения (IV.54).

Подставив значения

Сх и С2 в выражения

(IV.57),

полу­

чим

 

 

 

 

 

 

 

 

 

 

ßj

=

j * t—Фі W Фа (ц ) +

Фг W Фі f

(») d u

(IV.58)

 

 

То

 

 

 

 

 

 

 

где переменная интегрирования обозначается буквой

и для

отличия

от

переменной

т.

 

 

 

 

 

Если

ввести

обозначение

 

 

 

 

 

 

 

,

.

1 І Ф і М

Ф г ( т )

(ІѴ.59)

 

 

q{%, и) =

Л (и)

Ф і ( « )

Ф г ( " )

 

 

 

 

 

 

 

 

 

 

 

 

 

го формула (IV.58) примет следующий окончательный вид:

т

 

Рг (т) = J ^ ( T , и) F {и) du.

(IV.60)

r„

Д л я выяснения физической сущности q (т, и) рассмот­ рим случай подачи на вход системы сигнала в виде дельтафункции, действующей в момент т = т0 :

F(u) = ô(u — То). Применив формулу (IV.60), получим

j q (т, u)ô(u~

т 0 ) du = q (т, т0 ).

(IV.61)

То

Таким образом, q (т, т0 ) есть импульсная переходная функ­ ция системы, описываемой уравнением (IV.54).

Построим теперь решение уравнения (IV.54), которое является хорошо известным уравнением Бесселя. Линейно независимые частные решения уравнения (IV.54) будут равны

ф 1 ( т ) = ./1 (т), ф 2 ( т ) = УѴ1(т),

(ІѴ.62)

где J у (т), Ni (т) — функции Бесселя первого и второго рода первого порядка.

142

Используя дифференциальную формулу [20, 62] для

функций Бесселя dJ г

-1 х

„ „, - g -

—„ 1 — „ 0

и элементарное функциональное

уравнение JxN0— J0NX =

2

 

 

 

= - ^ - . найдем выражение для

определителя

Вронского

А(т)

=

 

(IV. 63)

Подставив формулы (IV.62) и (IV.63) в выражение (IV.61), получим импульсную переходную функцию движения внут­ ренней рамки ГВ в режиме запуска гиромотора на вираже объекта

Яи (т. т0 ) = -

т о lJi (т) Л/х (т„) -

(г) Л (т.,)].

(IV.64)

Из этой формулы видно, что импульсная

переходная

функция системы

с

переменными параметрами,

в

отличие

от системы с постоянными параметрами, зависит не только

от времени т,

но и от момента т 0

приложения входного им­

пульса.

 

 

 

 

 

 

 

Перед тем,

как перейти к построению реакции системы

на входной сигнал R (т, т0 ), для сокращения записи введем

вспомогательные

функции

 

 

 

 

Чп (т> и)

=

п2

и \3Х

(т) Nx

(и)

(т)Л(«)];

 

Чоі (т> и)

=

п2

u[J0

(т)Л^! (и) - І Ѵ 0 ( т ) / , ( « ) ] ;

 

 

 

л

 

 

(")•- Л ^

(т) J0(u)Y,

(IV. 65)

9 ю ( т . ")

=

2~ " [ • М ^ о

 

loo (т. ")

=

я

u[J0

(т) УѴ0

(и) ~N0(x)

J0(u)].

 

Y

 

Подставив выражения для R (т, т0 ) в соотношение (IV.60), представим, следуя [54], реакцию системы на входной сиг­ нал R (т, т0 ), обусловленный начальными условиями (IV.53), как общую импульсную переходную функцию

W(x, т0 ) =•

_!_

R

I

Éh.

 

х0

р

ю +

dx

_ L

R

 

I Éh.

<7іі(т До) — ßioVioC^o)- (IV.66)

т 0

P l 0

+

dx

x—xa

 

 

 

143

Легко убедиться, что для ГВ в установившемся

режиме

импульсная, общая импульсная и ступенчатая

переходная

функции по ß 1

будут равны

соответственно

 

 

 

 

 

 

 

 

<7со(т,т0) =

sin(T — т 0 ) ;

 

 

 

 

 

 

 

 

 

 

sin (т — т0 ) —ß0 cos (т — т0 );

(IV.67)

 

 

 

 

/ICO(T,T0 ) = _[ <7C» (T, и)

X

 

 

 

 

 

 

 

 

То

 

 

 

 

 

 

 

 

 

X

du

=

1 — cos (т —

т 0 ) .

 

 

 

 

 

Д л я

 

сравнения

на

 

 

 

 

рис. 34, а, б представле­

 

 

 

 

ны

импульсные

и

сту­

 

 

 

 

пенчатые

 

переходные

 

 

 

 

функции

для

ГВ в пус­

 

 

 

 

ковом и установившемся

 

 

 

 

режимах

на вираже объ­

 

 

 

 

екта. Отметим, что об­

 

 

 

 

щая

импульсная

 

пере­

2\

 

 

 

ходная

функция

систе­

 

 

 

мы

по

ß

(по а)

описыва­

 

 

 

 

 

 

 

 

ет

поведение

ГВ

с

 

вы­

 

 

 

 

ключенной системой кор­

 

 

 

 

рекции

на вираже

 

под

 

 

 

 

действием начальных ус­

 

 

 

 

ловий,

т.

е.

описывает

 

 

 

 

свободные

колебания

 

 

 

 

гироскопа.

Вынужден­

 

 

 

 

ные

колебания

ГВ

 

под

 

 

 

 

действием моментов кор­

 

 

 

 

рекции

найдем

как

 

ре­

 

 

 

 

акцию

системы

на

вход­

 

 

 

 

ной сигнал / (т) при по­

мощи интеграла

^ ап(%,

u)f(u)du.

 

 

 

 

 

 

 

 

 

 

 

То

 

 

 

 

 

 

 

 

 

 

 

Таким образом, с учетом (IV.59) решение

уравнения

(IV.54) будет

таким:

 

X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h=W

(г,

г,,) +

J <7н (т, и) f (и) du.

 

 

(IV.68)

144

Выражения для ai найдем, зная ßx , из уравнения (IV.60). После ряда элементарных преобразований, используя тож­ дества

J qu

(т, и) du =

1 +

 

д10

(т, т0 ) -

 

2*Щ.

+

 

j "

T

 

 

 

 

 

 

 

 

 

 

 

T

 

 

j q01

(t,

a) d« =

q00

(T,

T„) + - f -

^ Ц р ^ - +

J

du,

То

 

 

 

 

 

 

 

 

 

 

 

T0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(IV.69)

выражения

для ai,

ßi представим в виде

 

 

 

 

« 1

«1у =

( » 1 0

a 2 y )

Çoi — (ßio — ß2y) ÇoO' (IV.70)

 

 

ßl ßly =

(«и — <*2y) <7ll — (ßlO — ß2y)

O10,

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r,

~a

 

M

Ш г

ш *

4

-

ТпШц Sign ß,

+

 

 

«.y - «le. +

 

^

 

 

 

 

/ x n M | f t s i g n P l

 

 

 

 

 

 

 

 

 

\

 

 

 

W.v

 

tùy /-J

"

 

 

 

 

 

 

 

 

 

 

 

 

 

 

to

 

 

a2 y = a l o o

4- —

+ •

 

 

 

 

 

 

ßio —

ßi« —

 

 

 

"1/

 

 

 

ха>і/

 

т0

 

 

 

 

 

 

 

Ply =

Pi«, — —

 

 

Н

 

 

 

 

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

П

\

 

 

Cûy

 

CD, / J

и

 

 

 

 

 

 

 

 

 

 

 

 

 

То

 

 

 

 

 

 

о

 

о

 

 

w.t

t n 2 f t sign a.].

 

 

 

 

 

P2y =

Pico

 

 

 

 

 

Аналогичные выражения дл я а х и ß b описывающие дви­ жение главной оси ГВ в установившемся режиме при вира­ же объекта, имеют вид

ßi — ßiy =

( a io —a i y ) s i n ( T T o) +

(ßi0 — ßiy) cos (T

т0 );

a i «iy =

К о a i y ) cos (T To)

1 0 — ßiy ) sin (T

т„),

 

 

 

(IV.7I)

145

где

 

 

 

 

Юг

+ colft

sign

ßi

 

 

 

 

 

 

 

ßiy = ßI m

x

CU

Sign «!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из формул (IV.70) и (IV.71) следует, что в пределах

каждой четверти фазовой

плоскости

Oidißi

(рис.

35) при

 

 

 

 

 

вираже

 

объекта

 

главная

 

 

 

 

 

ось

ГВ

 

в

установившемся

 

 

 

 

 

режиме

 

будет двигаться по

 

 

 

 

 

окружности,

а

в

 

режиме

 

 

 

 

 

запуска — по

 

сходящейся

 

 

 

 

 

спирали. Д л я

больших зна­

 

 

 

 

 

чений т,

воспользовавшись

 

 

 

 

 

асимптотическими

 

пред­

 

 

 

 

 

ставлениями

функций Бес­

 

 

 

 

 

селя,

будем

иметь

 

 

 

 

 

 

<7u =

 

V -

sin( t — т0 ),

 

 

 

 

 

 

 

 

 

 

 

 

 

8 ѵ-т0

 

 

 

 

 

 

То

cos

 

(т — т0 );

 

 

Рис. 35.

 

 

?10

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

] / ^ С 0 8 ( Т

— Т0 ),

?00

j / ^ s i n (

 

 

т0 ).

(IV.72)

 

Подставив соотношения (IV.63) в выражения

(IV.70) и

исключив параметр (т — т0 ),

получим

 

уравнение

 

фазовой

траектории для ГВ в режиме запуска

 

 

 

 

 

 

 

К

-

а1 у )« + (ßi -

ßl y )2 =

-f-

[(а1 0 -

а2 у )2

+ (ß10

 

ß2y)2].

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(IV.73)

Аналогичное уравнение д л я ГВ в

установившемся

 

режиме

имеет

вид

 

 

 

 

 

 

 

 

 

 

 

 

К

-

а1 у )2 + (ßi -

ßl y )2 = ( а 1 0

- а1 у )2

+

1 0

-

ßl y )2 .

(IV.74)

 

Из

выражений

(IV.70),

(IV.71),

(IV.73),

(IV.74)

 

следует,

что центр спирали и центр окружности скачкообразно ме­ няют свое положение на фазовой плоскости при изменении знаков ai и ßi-

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

146

Следовательно, радиус спирали или окружности при изме нении знака aj или ßj скачкообразно получает определенное приращение, причем приращение радиуса спирали может быть значительно больше приращения радиуса окружности. Это значит, что главная ось ГВ в режиме запуска, при про­ чих равных условиях, быстрее, чем в установившемся режиме, будет приходить в положение равновесия. Преиму­ щество ГВ в режиме запуска при вираже объекта тем ощу­ тимее, чем больше переключений системы коррекции про­ изойдет за время разгона гиромотора. Характер движения гировертикали на вираже показан на рис. 35.

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

статора на наружную ось \—г sin ß] , который, действуя

аналогично коррекции по рамке, сильно демпфирует колеба­ ния ГВ.

§6. Влияние моментов сил сухого трения

иинерционности маятников на автоколебания ГВ

Во многих случаях при тщательных лабораторных ис­ пытаниях наблюдаются автоколебания гировертикалей и курсовых гироскопов с коррекцией от ЖМП . Отличитель­ ной особенностью автоколебаний является малая амплитуда (несколько угловых минут) и большой период (десятки се­ кунд) [7] . В случае использования гироприборов в режиме разгона ротора амплитуда автоколебаний на начальном этапе разгона достигает нескольких градусов, а период уменьшается до нескольких секунд. Это существенно влияет на точность показаний гировертикалей (см. рис. 36). Как будет показано ниже, инерционность Ж М П при наличии моментов сил сухого трения в опорах карданового подвеса является основной причиной возникновения в некоторых случаях устойчивых автоколебаний гировертикали относи­ тельно положения равновесия. Впервые одна из постоянных времени Ж М П была определена в работе [55], в которой ЖМП представлен в виде апериодического звена первого порядка. В работе [7] Ж М П представляется в виде звена второго порядка.

147

При исследовании движения ГВ в режиме разгона ро­ тора будем пользоваться методом «замораживания» коэф­ фициентов.

Уравнения движения ГВ вблизи положения равновесия запишем в виде

Ясс — /СхѴаН-Мат sign Р = 0;

(IV.75)

— Щ — Яр + /C2 vß - f M I T sign а = 0.

К этим уравнениям необходимо добавить уравнения дви-

Л

град

2

4

S

8

/0

12

/4 %сек

Рис. 36.

жения жидкости Ж М П [7] в пределах зоны пропорциональ­ ности

(ТУ

+

2g7> +

1) ѵа

=

- (ТУ

+

1) а;

 

(ТУ

+

2£Тр +

l ) v ß

=

(ТУ

+

1) р,

(IV. 76)

где 7\ і — постоянная времени и относительный коэффи­ циент демпфирования Ж М П , Т = 0,4 сек, I = 6 [8] .

Системы уравнений (IV.75), (IV.76) содержат существен­ но нелинейные члены, обусловленные трением, наличие ко­ торого в данном случае может стать одной из причин авто­ колебаний гировертикалей. Покажем это. Используя ме­ тод гармонической линеаризации, выражения, содержащие

148

трение,

представим в

виде

 

 

 

 

 

 

соіт sign а = q±pa;

со2

т sign ß = ^2 pß,

(IV.77)

 

4Л*,Т

 

4 М 2 т

 

 

 

 

 

г д е Ii =

-ШІГ

1* =

ü t o F

> a > b

-

амплитуда

автоколе-

баний соответственно

наружной и внутренней рамок;

со —

частота

автоколебаний.

 

 

 

 

 

Подставив

(IV.77)

в исходные

уравнения,

после

ряда

элементарных преобразований получим следующее характе­ ристическое уравнение системы:

а0рв

+ ахръ + а2р4

 

+ а3ря + а^р2

+ аьр + ав

=

О, (IV.78)

где

 

 

 

 

 

 

 

 

 

 

 

а0

= Т%

а±

=

4СГ*с + (kt + k2) 7 м +

k0T\

а2

=

2 £ Г 3 (kx

+ k2) + T%k2

+

2T2c(l

+ 2g2) +

 

 

 

 

+

k0kxT* +

4t,T%,

 

 

 

a3 = 4£Tc + 2(k,+

 

k2) T2 + 4t?T%

+ 2T%

+

kak%T\

at = с 4- 2kxk2T2

 

+ 21T (kx + k2)

+ 4СГк0

+

 

2k0k1T2;

«5 = k l + К + К + 4 T k 0 k V

a0 = k A +

k0kl>

 

 

« i , 2 = -jf-

 

, k0 = -yj- ; с = 1 - f qiq2.

 

Д л я отыскания параметров автоколебаний подставим в ха­ рактеристическое уравнение (IV.78) р = /со и приравняем вещественную и мнимую его части нулю

ßjCo5 — а3 со3 -(- а5со = О,

 

— а0со° 4- а2 со4 — а4 соа 4- а 0 = 0.

(IV.79)

Из первого уравнения (IV.79) получим

®І2

2а,

Подставляя в это выражение значение коэффициентов ха­ рактеристического уравнения и учитывая, что

* о О і - 2 ; (kx +

k2)T*^Tc,

после ряда элементарных преобразований получаем сле­ дующие приближенные формулы для определения частоты

149

автоколебаний:

 

 

 

 

 

 

г,2

~

1 .

г? ~

/ г і + К + К + ilTkJi^

 

.

,тѵ

я т

wi

« - j T j - ,

2 Ä ;

^рт^

 

(IV.ou;

Подставляя

 

значение частот (IV.80) во второе

уравнение

(IV.79), легко убедиться, что для частоты ах

с

= 0,

что

не

имеет физического

смысла.

Д л я частоты

со2

получается

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

 

 

 

V 2

+

V

+

ö 2

=

0,

 

 

(IV.81)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ь0 = я с г 2 -

г, Ъ, =

2 Г 2

( 1 + 2g2) -

г (а., -

с),

 

 

 

Ь2 = а 2

- 2 Г с ( 1 + 2 ? ) - ^ ,

 

г = ^ - .

 

 

Очевидно,

что

при

указанных

значениях

параметров

Ж М П

(или близких к

ним)

Ьх

>

0,

Ь2 <

0.

 

 

 

Необходимым и достаточным условием существования

автоколебаний

при

с =

1 +

 

 

будет неравенство

с >

1,

которое может выполняться как при Ьи >

0, так и при Ьи

<

< 0. Можно

показать,

что

условие

с >

1 сводится

к

не­

равенствам

 

 

 

 

 

 

 

 

 

 

 

 

 

+ \Ь2\>Фо

 

+ ьі)>0

 

 

при

Ь 0 > 0 ;

 

 

 

+

I b21 >

I Ь 0 — Ъг I

 

при £>0 <

0,

 

 

которые выполняются в довольно широком диапазоне

пара­

метров

гировертикалей.

 

 

 

 

 

 

 

 

 

Отметим, что автоколебания будут наблюдаться и в том

случае,

когда одна

из коррекций

отключена, но частота

их

будет

меньше.

 

 

 

 

 

 

 

 

 

 

 

Приведенные формулы

для

частоты и

амплитуды

авто­

колебаний являются довольно сложными. В некоторых слу­ чаях определенное упрощение дает следующее представле-

ление уравнений движения Ж М П

[55]:

 

 

 

 

( 7 > +

1) ѵ« =

-

а,

2р

+

1) v ß

=

- ß,

(IV.82)

где Т2 =

2\Т.

 

 

 

 

 

 

 

 

 

 

Тогда

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

уравнение

примет

вид

 

а 0 р 4

+

аіР*

+

а2р2

+

а3р

+ а, =

0,

(ІѴ.83)

где

 

 

 

 

 

 

 

 

 

 

 

 

а0

=

Т\с;

 

=

2 7 V

+

kjl

 

 

150

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