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

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

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

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

Добавлен: 10.04.2024

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

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

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

g 4.1]

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

І75

Формулы (45) и (47) отличаются от формул (23) и (25), полу­ ченных для сухого трения при отсутствии вращения подшипни­ ков тем, что теперь под знаком интеграла явно входит время и, следовательно, функция MTX(t) не является стационарной.

Умножение характеристической функции на множитель вида еіиЬ эквивалентно прибавлению к математическому ожиданию случайной величины (функции) неслучайной величины Ь. По­ этому формулы (45) и (47) можно переписать в виде

 

т (t) =

^ Е і (и) ~ + т°,

(4.48)

1 ’ У --

Ql

 

 

2

 

 

 

Tl

 

 

 

—СО

— [тй(У — fn°][m{t2) ™Р],

(4.49)

 

 

где штрихом отмечены характеристические функции ординат

случайной функции [9 (£)+avcos о>0£ ]. Интегралы в формулах (45) и (47) по виду не отличаются от выражений (23) и (25), и следова­ тельно, для вычисления m(t) и Km(tv у можно воспользоваться формулами, аналогичными (31) и (36):

ln(t) = m° + QJl-2F'i>(0)],

(4.50)

Кт(У У - 40* [Fi і, (0, 0) - Ft, (0) Fi, (0)J +

 

+ D [M°rx][in (У — m°] [m (t2) — m°\,

(4.51)

где штрихом отмечены функции распределения ординат случай­ ной функции [9 (£)-[-avcos w0t ].

Если, например, 0 (t) — стационарный нормальный процесс, то

^ ( 0 ) = 4 1 +

COS

С0|)1

(4.52)

 

 

Гп (t) = ,п° + 0 ЖФ (

.

 

(4.53)

где Ф (ж) — обозначение функции Лапласа. Для Fe„è,(0, 0 ) в этом случае имеем

 

ѵП - к* ' 2er? (l-fc2)

 

(ac-h2)

 

 

 

х/2

 

іЛа(і-кг)

X

 

г і,л ( М ) = $

2ъа

2оga\/2 ua е ѳ

 

 

 

X

1 — б»( ___ ь-

___ )]

d-s

 

 

 

VafiV/aV/l

к* )

 


17(i ГУ , О П И С Ы В А Е М Ы Е Л И Н Е Й Н Ы М И У Р А В Н Е Н И Я М И ІГЛ. 4

где

 

 

 

а — 1 — Л: sin 2ср,

Ъ= — у | (cos <»+ — к cos <»0t2) cos f

+

 

+

(cos Ш0t2к cos «>+) sin tp],

c = cos2 i»0tj +

cos2c»0£2 2/c cos (»+

cos со0t„, к =

К j (x)/o?.

Дальнейшие вычисления требуют применения численного инте­ грирования.

Аналогичным образом вычисляется m(t) и Km(t\, t2) в том слу­ чае, когда подшипники вращаются в противоположных направле­ ниях с постоянной скоростью, изменяя направление вращения через каждые п оборотов (подшипник типа «роторейс»). При этом формулы (50) и (51) остаются в силе, но в них нужно заменить т°, D [М®х ], Qx значениями параметров, соответствующих харак­ теру движения опор в подшипнике типа «роторейс», а частоту со0 на частоту шх [см. (2. 109)].

Сравнение формул (50) и (51) с формулами (31) и (36) доказы­ вает, что в случае вращения подшипников выражения для rn(t) и Ат(т) приобретают колебательный характер, так что при инте­ грировании по формулам (5) и (6 ) можно ожидать уменьшения значений соответствующих интегралов сравнительно с тем слу­ чаем, когда нет вращения подшипников.

Аналогичным образом могут быть подсчитаны первые моменты случайной функции a{t) при колебательном движении опор на оси внутреннего карданова кольца в противофазе, применяемого иногда для снижения отрицательного влияния моментов сил су­ хого трения в подшипниках. Момент сил трения в этом случае

имеет вид [см.

(2. 105)J

 

 

М „ (*) = ML +

Ql sign (t) +

a, C0S <Vl +

 

 

+

Q2sign [Ѳ(t) — av cos ш0£],

(4.54)

где M^x — знакопостоянная составляющая момента трения, а зна­ копеременные составляющие моментов трения и Q2 могут быть различными из-за неидентичности подшипников.

Применяя общие формулы (50) и (51), получим

т (t) =

іи» - Q

x\1 -

2F-(0)] -

Q2[1 - 2F" (0)],

(4.55)

A„(/„ У = а д / '’ііЛ (0 ,

0 ) _

2 /^ ( 0 ) /^

(0 ) 1 +

 

+

Ш

П

л ( ° ’

° ) -

2 F l

(°) (°)l +

 

+ WiQ* \n iA (0, 0) -

2F l (0) F l (0)1 +

 

 

+

D \M°rx][m + ) — m°\\tn (£,,) ■— m°],

(4.56)

где штрихом отмечены функции распределения случайных вели­



§ 4.1]

У Р А В Н Е Н И Е 1-ГО П О Р Я Д К А

Б Е З ЗАВИСИМ ОЙ П Е Р Е М Е Н Н О Й

177

чин

[0 (<x)+avcos «Vj 1 и [Ѳ(£2)+ а ,cos u)0i2 l;

двумя

штрихами —

случайных величин [О(tj)—avcos u)0£ j и [0 (t2)—avcos ш ^],

а тремя

штрихами — функции распределения случайных величин

[Ö(ft)-f-

- ha,cos Шді?!] и

[6 (f2)—a,cos м0і2].

для ä(t)

и D[a(<)l,

Получены

приближенные

формулы

показывающие, что при применении разновращения подшипни­ ков дисперсия ухода оси гироскопа направления уменьшается.

5. Уход астатического гироскопа, вызванный неуравнове­ шенностью и другими факторами. Выше мы рассмотрели диф­ ференциальное уравнение (1 ), характеризующее уход астати­ ческого гироскопа при учете момента трения M rx(t) в оси подвеса. Аналогичное уравнение возникает и в некоторых других случаях, при которых уход гироскопа связан не с наличием момента трения в осях подвеса, а с другими возмущающими моментами.

Рассмотрим два наиболее интересных случая: уход ГН вслед­ ствие статической неуравновешенности гироскопа и уход астати­ ческого гироскопа вследствие упругой деформации ротора [б3]. В первом случае в соответствии с формулой (3.29) будем иметь

(при

ß0 = 0 )

 

 

 

 

 

(4.57)

где

Р — вес гироскопа (ротор+ кожух);

Іг — расстояние

вдоль

оси ротора от центра тяжести гироскопа до точки подвеса;

пт. —

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

 

 

Уравнение (57) отличается от уравнения (1) формально только

тем,

что вместо момента трения М Тх (t)

здесь стоит выражение

 

 

 

(4.58)

Следовательно, все формулы этого параграфа останутся в силе, если входящую в них корреляционную функцию Km(tv t2) за­ менить на K v {tx, t2), а m(t) на у (t). Будем считать дебаланс Іг случайной величиной с математическим ожиданием, равным нулю, и заданной дисперсией а|ж, а вертикальное ускорение (t)

стационарной случайной функцией времени, не зависящей от Іг и имеющей нулевое математическое ожидание. В этом случае, применяя к обеим частям равенства (58) операцию нахождения математического ожидания, получим

(4.59)

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

12 А. А. Свешников, С. G. С’ивкин


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

будет сходиться к нулю при увеличении числа приборов. Исходя из определения корреляционной функции (1.60), имеем

K y {tx, g

= p ^ [ i

+ l / c , , w ] = A^(x)

(X= * ,-* ,).

(4.60)

Таким

образом,

корреляционная функция К у {і) будет

зави­

сеть только от разности моментов времени и, следовательно, функция Y (t) будет стационарной.

Корреляционная функция (60), а также вычисленная по этой формуле дисперсия ухода а(t) характеризуют разброс показаний различных гироскопов.

В некоторых задачах представляет интерес исследование ошибки конкретного, взятого для испытания гироскопа. В этом случае величина Іг не меняется в течение всего опыта и, следовательно, задача сводится к определению условного математического ожи­ дания ошибки гироскопа и ее условной дисперсии при конкрет­

ной реализации случайной величины lz.

В соответствии с этим

вместо математического ожидания и корреляционной функции

У (і), определяемых формулами (59) и (60),

необходимо вычислить

условное математическое ожидание y^(t) и условную

корреляцион­

ную функцию Ky\is (tv t2) (при фиксированном

значении Іг),

которые в этом случае будут иметь вид

 

У1ж«) = Р1„

(4.61)

і2) = Р Ч І ^ К ш^ ) = К уѴ^ ) .

(4.62)

Во втором из упомянутых случаев, т. е. при рассмотрении ухода астатического гироскопа вследствие упругих деформаций ротора, в соответствии с формулой (3. 30), имеем (при ßo=0)

®W = ТПГ7 (с- — “V« W w*(*)>

(4-63)

где сг и сУі — коэффициенты жесткости ротора вдоль соответст­ вующих осей, а wSl (t) и wz (t) — проекции вектора ускорения точки

подвеса на экваториальную ось и ось собственного вращения ротора. Следовательно, обозначив

z (*)= г г (с* ” су») “V. (*) w*(*)>

(4-64)

мы снова можем использовать все формулы настоящего параграфа, заменив т (t) на z(t) и Km(tv t2) на t2). Будем считать коэф­ фициенты жесткости неслучайными величинами, а случайные функции wVl (t) и wz (t) стационарными с нулевыми математиче­ скими ожиданиями. Предположим, что эти функции, кроме того,