Файл: Свешников А.А. Вероятностные методы в прикладной теории гироскопов.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) стационарными с нулевыми математиче скими ожиданиями. Предположим, что эти функции, кроме того,