Файл: Свешников А.А. Вероятностные методы в прикладной теории гироскопов.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 10.04.2024

Просмотров: 257

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

164

ОСНОВНЫЕ У РАВНЕНИ Я ПРИКЛАДНОЙ ГИРОСКОПИИ

[ГЛ. 3

К уравнению типа (367) можно отнести, например, второе уравнение (24) ГН с коррекцией, уравнение (58) ГВ с маятниковой коррекцией, уравнение (141) ГТ, в котором используется трех­ степенной астатический гироскоп, уравнение (173) поплавкового интегрирующего гироскопа, уравнение (207) гироскопического интегратора линейных ускорений объекта.

3) Уравнение второго порядка

ä -f- -f- а2а = X (t), (3.368)

где ах и а2 — постоянные коэффициенты.

К уравнениям типа (368) можно отнести, например, уравне­ ния (42), (55), (56) колебаний физического маятника, уравнение (6 6 ) ГВ при наличии статической неуравновешенности гироскопа, уравнения (109), (110) ИВ с демпфированием, уравнение (122) ГТ, уравнение (149) ВГ, уравнение (270) колебаний корабля с гиро­

скопическим успокоителем

качки активного типа, уравнение

(288) ГК Фуко, уравнение

(320) ГК с косвенной коррекцией,

уравнение (352) акселерометра. Уравнение (368) является наи­ более распространенным типом дифференциальных уравнений, характеризующих движение ГУ.

Возможны случаи, когда в уравнении (368) коэффициенты

или а2могут быть равны нулю. Так при

ах= 0 имеем уравнение

ä -f- a2a = X (t).

(3.369)

Куравнениям этого типа относятся, например, уравнения

(81)ГМ без демпфирования, уравнения (95) ИВ при отсутствии устройства для обеспечения демпфирования колебаний гироскопа.

Если в (368) положить а2 =0, то получим

&+ a1i = X(t).

(3.370)

Этому типу соответствует, например, уравнение (164) поплав­ кового интегрирующего гироскопа.

4) Системы уравнений с постоянными коэффициентами. Вид этих уравнений может быть различным. Здесь в качестве примеров можно указать уравнения (6 ) астатического трехстепенного гиро­ скопа, уравнения (77) ГМ, уравнения (102) ИВ с учетом стати­ ческой неуравновешенности гироскопа, уравнения (230) одно­ осного ГС на качке, уравнения (261) гироскопического успокои­ теля качки пассивного типа, уравнения (293) однороторного ГК.

В некоторых случаях при анализе ГУ приходится учитывать случайный характер коэффициентов линейных уравнений, опре­ деляющих поведение ГУ. При этом коэффициенты могут быть случайными величинами или случайными функциями времени. К уравнениям первого типа можно, например, отнести уравнения (28) ГН с учетом статической неуравновешенности ротора, урав­


§ 3.7]

ОСНОВНЫЕ ТИПЫ УРАВНЕНИЙ ГУ

165

нения (64) ГВ при наличии также статической неуравновешенности гироскопа, уравнения (105) ИВ с учетом изменения ее параметров.

Примерами второго типа уравнений, когда коэффициенты уравнений являются случайными функциями времени, могут служить уравнения (34) ГН на качке, уравнение (57) колебаний физического маятника, уравнения (85) ГМ, уравнение (114) ГТ, уравнение (163) поплавкового интегрирующего гироскопа, урав­ нения (248) трехосного ГС, уравнения (303) двухроторного ГК.

При исследовании ряда ГУ приходится учитывать нелиней­ ность уравнений движения вследствие того, что искомые пере­ менные являются аргументами нелинейных выражений, входящих в уравнения, характеризующих поведение ГУ. Такими уравне­ ниями являются, например, уравнения (1 ) трехстепенного аста­ тического гироскопа, уравнения (34) ГН на качке, уравнения (74) ГВ с нелинейной характеристикой коррекции, уравнение (ИЗ) ГТ, уравнение (162) поплавкового интегрирующего гироскопа.

Сложные гироскопические системы, например, гирокомпасы, гирогоризонткомпасы, трехосные силовые гиростабилизаторы, инерциальные навигационные системы характеризуются систе­ мами уравнений высокого порядка (см. гл. 7), содержащими, как правило, нелинейные члены, случайные величины в коэффициен­ тах уравнений и отдельные коэффициенты уравнений, имеющие характер случайных функций. Исследование таких сложных урав­ нений часто удается свести различными приемами к исследованию уравнений более простого типа, упомянутых ранее. Поэтому мы начинаем с рассмотрения простейших уравнений ГУ, постепенно увеличивая их сложность.

Г Л А В А 4

ИССЛЕДОВАНИЕ ГИРОСКОПИЧЕСКИХ УСТРОЙСТВ, ОПИСЫВАЕМЫХ ЛИНЕЙНЫМИ ДИФФЕРЕНЦИАЛЬНЫМИ УРАВНЕНИЯМИ

§ 4.1. Линейное уравнение первого порядка, не содержащее зависимой переменной

1. Вывод общих формул. Начнем исследование дифферен­ циальных уравнений, описывающих поведение ГУ, с простейшего уравнения первого порядка, не содержащего зависимой перемен­ ной. Рассмотрим для определенности уравнение движения гиро­ скопа направления (ГН) в рамках прецессионной теории, т. е. уравнение (3.24)

(4.1)

где а — угловая скорость ухода гироскопа; М — возмущающий момент, являющийся случайной функцией времени, а Н — кине­ тический момент гироскопа, обычно принимаемый постоянной величиной. При исследовании инструментальных погрешностей гироскопа H(t) нужно считать случайной функцией времени. Предположим, как это обычно и имеет место, что непостоянство кинетического момента вызвано небольшими отклонениями в ско­ рости вращения ротора гироскопа и, следовательно, отклонение H(t) от номинального значения кинетического момента Н0 можно считать малым. Поэтому положим

Я(Й = Я0 + Я 1 (Й>

(4.2)

где H^t) будем считать малой величиной в том смысле, что вто­ рыми и более высокими степенями Н^і) будем пренебрегать, т. е.,

1

H^ t)

1

Hi(t)

В качестве возмущаю-

например, считать Ho +

 

-----Щ~

щего момента М будем учитывать момент сил трения М тх в оси вращения внутреннего карданова кольца подвеса ГН. В этом

случае уравнение (1 ) можно

представить в

виде

а ~ Н п

' Н~ Н г (t)Mti

(4.3)

Предположим, что математические ожидания и корреляцион­ ные функции момента сил трения m(t) и Km (tv іъ) и случайного



§ 4.11 УРАВНЕНИЕ 1-ГО ПОРЯДКА Б Е З ЗАВИСИМОЙ ПЕРЕМЕННОЙ Ш

отклонения кинетического момента H x{t), hy(t) и Kh(tv t2) известны,

а начальное значение угла а,

H^t)

и М тх(t)

взаимно независимы.

В этом случае для первых

двух

моментов

случайной функции

а(і) могут быть получены простые расчетные формулы. Действи­

тельно, интегрируя

обе части

уравнения

(3), будем иметь

 

t

t

 

* (0 = * (0) +

$ м тх (tj) dt1~щ\

Я 1 (*i) Мтх {tx) dtv (4.4)

 

о

о

 

Находя математическое ожидание обеих частей последнего равенства, пользуясь формулой (1.91) для корреляционной функции суммы независимых слагаемых и используя формулу (1.87) для корреляционной функции интеграла от случайной функции, получим

 

 

 

 

 

t

t

 

ä (*) =

ä (0 ) +

_i_ J fn (fj) dtx— —- j h, (tv) in (fj) dtv

(4.5)

 

 

 

t21

0

0

 

 

 

 

 

 

 

K , ( * L h ) =

Щ 5 5 K m ( Xl - h ) d x l d x 2 +

 

 

 

f2

0

0

 

 

 

 

 

 

 

 

 

+ Щ

\

\

{K h (Xl>

S ) K m (Tl. T2) +

К ( T j ) h 2 (v2) K m ( T j , T 2)

+

0

0

0

 

 

 

 

 

 

к к

+ m ( T i ) m ( t 2 ) K h ( x l5 T 2) j d x , d x 2 —

 

 

 

 

x 2) dxl(7x2 + D [a (0)],

 

— щ

S

S VK (Xj) +

(X2)] K m (xlf

(4.6)

0

 

о

0

 

 

 

 

где при определении корреляционной функции от интегралов при­ няты во внимание соотношения (обозначим для краткости Z (t) =

КЛХѴТ2) = Kh(Xl, \)Km(Xl, т2) +

+ К (xj) h Y(x2) K m (xlt x2) + m (хД rh (x2) K h {zv x2)

(4.7)

Д»ж(т1 . T2) = Âi (X2) Km( X j, x2),

получающиеся, если воспользоваться формулами (1.62) и (1.63) для корреляционной функции и для взаимной корреляционной функции.


168

ГУ, ОПИСЫВАЕМЫЕ ЛИНЕЙНЫ М И УРАВНЕНИЯМИ

[ГЛ. '<

Положив в (6 ) t^= £2= tj находим формулу для дисперсии угла a.(t) в момент времени t:

t t

D [а (0] = £ r j 5 K m(tlt t2) dx,dx2+

00

it

5 (^A (ti. *2)я„(*і. Ts) + M ii)M T2) ^ h , ^ ) +

0 о 0

 

 

 

 

+ nUx1)m(x2) K h(xv

TjJJdxjdtj,—

 

t

t

 

 

— ■Щ S

S [^I (xi) + h (xa)] K rn (xi>

x2) dxi dx2 + D [ а (0)].

(4 .8 )

о0

Втом случае, когда случайные функции H1(t) и MTX(t) ста­ ционарны, формулы (5) и (8 ) упрощаются и принимают вид

где при выводе последней формулы введены новые переменные интегрирования £= т2 и т= т2— и выполнено интегрирова­ ние по £.

Таким образом, первые два момента ухода гироскопа полно­ стью определяются двумя первыми моментами случайных функ­ ций H x{t) и М (t).

Формулы (8 ) и (9) показывают, что даже при стационарных случайных функциях H^t) и M TX(t) уход гироскопа а(t) не явля­ ется стационарной случайной функцией.

Математическое ожидание отклонения кинетического момента гироскопа H(t) обычно может быть принято равным номинальному значению кинетического момента Н0, а Dla(0)] можно не учи­ тывать, так как нас интересует ошибка, накопившаяся за время t,

§ 4.1] У РАВНЕНИ Е 1-ГО ПОРЯДКА Б Е З ЗАВИСИМОЙ ПЕРЕМЕННОЙ

169

В этом случае hx= 0,

а формула (9) еще больше упрощается и при­

нимает вид

 

 

 

 

t

 

 

 

 

D[cc(f)] = ^ S (t — x) Кт{х) dx +

 

 

 

О

t

 

 

 

 

 

 

 

+

X) [Kh(x) Km(X) +

т*К„ (X)] dx.

(4.10)

 

0

 

 

 

Корреляционная

функция Kh( x) (или

Kh [tx,

t2) — если H(i)

не обладает свойством стационарности) должна

быть определена

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

Рассмотрим три основных случая.

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

MTx = nb(t),

(4.11)

где Ѳ( 0

составляющая угловой скорости вращения объекта

вдоль оси,

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

Используя формулу (1.74), находим

 

 

 

К т W = п * К 6 (X ) =

.

(4.12)

Подставляя последнее выражение в (10), учитывая в этой формуле только первое слагаемое и производя интегрирование по частям, окончательно получим

о „о

(4.13)

D г* (01 = ^ - [ * 9 (О ) - Я . (01-

Если считать время t большим времени корреляции случай­ ной функции Ѳ(£) (т. е. большим времени, в течение которого кор­ реляционная функция Кц (t) практически затухает до нуля), то формула (13) даст

D M O ]= -? rD [0 (f)].

(4.14)

Следовательно, при наличии вязкого трения дисперсия ошибки гироскопа направления после окончания переходного процесса бу­ дет величиной постоянной, а средняя квадратическая ошибка