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

книги из ГПНТБ / Ковалев М.П. Динамическое и статическое уравновешивание гироскопических устройств

.pdf
Скачиваний:
46
Добавлен:
23.10.2023
Размер:
8.54 Mб
Скачать
Рис. 2.1. Схема действия си­ лы на точечную массу

Уравнение (2.8) выражает теорему о моменте количества движения тела. Согласно, этой теореме скорость конца вектора кинетического момента тела относительно неподвижной точки геометрически равна моменту действующих на тело внешних сил относительно той жё точки.

Если в твердом теле выбрать не­ которую точку, например центр масс, и совместить с этой точкой начало координат О системы £°т]0£0, оси ко­ торой параллельны соответствен­ ным осям инерциальной системы

то движение твердого тела можно рассматривать как совокуп­ ность двух движений: поступатель­ ного движения, определяемого дви­ жением выбранной точки тела, и движения (вращения) тела около этой точки относительно системы

S W

Уравнение (2.8) можно принять и при исследовании движения (вращения) твердого тела отно­

сительно системы І°г]0£0, которую будем в дальнейшем называть «неподвижной». При этом, как уже отмечалось, к моменту дей­ ствующих сил надо добавлять момент от сил инерции перенос­ ного (вместе с системой ё°Л°£0) движения тела. Общий момент М, а также кинетический момент тела G нужно брать теперь от­ носительно начала координат О неподвижной системы |°ц0?0.

2.2. ПОЛНАЯ ПРОИЗВОДНАЯ ВЕКТОРНОЙ ФУНКЦИИ ПО ВРЕМЕНИ

Предположим, что имеется некоторый вектор а, зависящий от времени t. Помимо неподвижной системы отсчета £°т]0(;0 вве­ дем еще систему координат XYZ (подвижная система) с тем же началом, вращающуюся относительно системы £°т]0£;0.

Будем

считать, что начало

вектора а совпадает с началом

координатных систем ^°т)0^0 и XYZ (рис.

2.2). Будем также счи­

тать,

что

подвижная система

вращается

относительно непод­

вижной с угловой скоростью со. Выразим

 

полную производную

da

т.

е. скорость изменения вектора

----- ,

а

относительно иепод-

dt

 

 

 

 

 

вижной системы, через величины, относящиеся __к подвижной

системе XYZ. В проекциях на оси подвижной системы вектор а можно записать следующим образом:

(2.9)

17

Дифференцируя (2.9) по времени, получим

da

,

' d a

I

dan

— ,

d a z 7-

dt

I 4------

J H------

* k

\

dt

'

d t

1

dt

d i s

d j ,

dk

( 2. 10)

« , ~ 77" "Ь a c/ ~TT 4" a z

dt

dt

dt

 

Выражение

da у —

— . rfaz

_da

(2. 11)

І-

dt ^ dt

dt

dt

 

представляет собой скорость изменения вектора а по отношению к подвижной системе XYZ. Эта ско­ рость называется локальной произ­

 

 

водной _вектора

а

и обозначается

 

 

как

d а

 

 

 

 

 

 

 

---- .

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

Вторая группа членов правой ча­

 

 

сти уравнения

(2.10)

при учете соот­

 

 

ношений

 

 

 

 

 

 

\ ° di

 

dj

— - dk

~ \у~й

 

 

———<« X I,—-

=ш X У, —— =«> X к

 

 

dt

 

dt

 

 

dt

 

 

 

и выражения

(2.9) может быть за­

 

 

писана в

следующем виде:

 

Рис. 2.2. Системы

координат

а,-- di

 

dJ +

az - ^ - = a x(uX i)+ -

dt

 

' dt

 

dt

 

-f <МШX у) + аг(м X Ä)=w X (avi- f ayj +

atk) = <»X

a.

Таким образом,

полная

производная вектора а по времени

определяется формулой

 

 

 

 

 

 

 

 

dt dt

 

X

 

 

 

(2.12)

где, как уже отмечалось, -4^- — скорость изменения вектора по.

отношению к подвижной системе координат XYZ, а со Ха — ско­

рость конца вектора а по отношению к неподвижной системе при вращении этого вектора вместе с подвижной системой XYZ.

2.3. ДИНАМИЧЕСКИЕ УРАВНЕНИЯ ДВИЖЕНИЯ ТВЕРДОГО ТЕЛА ОКОЛО НЕПОДВИЖНОЙ ТОЧКИ

Теорема о моменте количества движения твердого тела выра­ жается согласно (2.8) уравнением

— = 7Й,

(2.13)

18

в котором G обозначает кинетический момент тела относительно начала координат О неподвижной системы |°ті0£0, а М — момент сил (активных и инерционных от перекосного движения) относи­ тельно той же точки О.

Принимая в качестве вектора а вектор кинетического момен­ та G твердого тела, можно при учете формулы (2.12) записать уравнение (2.13) в виде

 

 

— 4- Ü X G —M.

 

(2.14)

Выразим векторное произведение через определитель

 

I» X о --

1

У

k

 

 

w.r

а и

Ш2

 

 

 

 

О,

Gy

Gz

 

После раскрытия этого определителя можно записать уравне­

ние (2.14) в проекциях подвижной системы XYZ

 

 

CtG; -f- ®yOziozGy= Mx;

 

 

dt

 

 

 

 

 

 

dGy

■“A

<*xGz = My\

(2.15)

 

dt

 

 

 

 

 

 

dG,

+ wxGy — %Gx= Mz

 

 

dt

 

 

 

I

 

Знак локальной производной здесь опущен, так как изменение

вектора G относительно подвижной системы

XYZ есть не что

иное, как скорость изменения проекций

этого

вектора на оси

X, Y, Z.

 

 

 

 

 

 

В дальнейшем в качестве твердого тела будем рассматривать

ротор (рис. 2.3).

Подвижную систему координат XYZ свяжем с

ротором следующим образом: ось Z направим по оси вращения

ротора, а оси X,

Y будем' считать не

 

 

 

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

 

 

 

ротор по отношению к этим осям

 

 

 

может поворачиваться. Угловую ско­

 

 

 

рость поворота связанных с телом

 

 

 

ротора осей Хи

относительно осей

 

 

 

X, Y обозначим через ср.

 

 

 

 

Можно доказать, что проекции

 

 

 

вектора G кинетического

момента

 

 

 

ротора на оси подвижной системы

 

 

 

XYZ определяются формулами

 

 

 

Ох= А ш х-

^

 

 

 

 

G y

>

 

(2.16) Рис. 2.3.

Ротор в системе

Gz —Cioz -ф-H ,

 

 

 

координат

19

Рис. 2.4. Схема гироскопа с перекошенными осями

где А, С — моменты инерции ротора 'соответственно относитель­ но экваториальных (X или У) и полярной (Z) осей; сщ-, соу, coz — проекции угловой скорости системы XYZ на оси этой же системы,

а Я = Сср— кинетический момент ротора при его вращении отно­

сительно системы XYZ (в дальнейшем

будем принимать ср =

= const).

следующие уравнения

Подставляя (2.16) в (1.15), получим

движения ротора относительно неподвижной системы £°г]0£;0, на­ зываемые динамическими уравнениями:

У Ц Ѵ~\~(С А) шу<лг Н (ау М х ;

Лму -} -(^ — С) со гш2— Н соЛ.= М у \

(2.17)

2.4. УРАВНЕНИЯ ГИРОСКОПА С ПЕРЕКОШЕННЫМИ ОСЯМИ ПОДВЕСА

Карданов подвес гироскопа ранее предполагался идеальным, т. е. кольца подвеса считались безынерционными, их оси враще­ ния взаимно перпендикулярными, а ось вращения ротора перперпендикулярной оси вращения внутреннего кольца. Однако вследст­ вие инерционности колец и возможного перекоса осей реальный трехстепен­ ный гироскоп в процессе собственных колебаний

няться от первоначально­ го положения даже при неподвижности его осно­ вания. Этот уход (дрейф) гироскопа при колебаниях основания возрастает.

Схема такого гироско­ па приведена на рис. 2.4. Она состоит из ротора с

осью симметрии Z\ и центром тяжести О, совпадающим с пере­ сечением осей карданова подвеса, а также из внутреннего и внешнего колец. Центры тяжести колец совпадают с точкой О, однако, несмотря на идеальность конфигурации колец, оси их вращения неточно перпендикулярны друг другу. Если Y' — ось вращения внутреннего кольца относительно внешнего, а У — ось в плоскости XY внешнего кольца, перпендикулярная оси X вра-

20

>

щения внешнего кольца, то угол е между Y' и Y представляет перекос внутреннего' кольца относительно внешнего.

Ротор ’предполагается перекошенным относительно внутрен­ него кольца. Перекос определяется углами уі и уг, характеризу­ ющими отклонение оси Z\ вращения ротора от главной оси инер­ ции Z внутреннего кольца, у і— угол отклонения в плоскости симметрии Y'Z' внутреннего кольца, а у2 — угол отклонения оси

ротора от этой плоскости симметрии. Последний угол можно рассматривать как отклонение оси Z\ вращения ротора относи­ тельно диаметральной оси Y' этого ротора, лежащей в плоскости симметрии Y'Z' внутреннего кольца.

Все упомянутые оси приведены на рис. 2.5. Система XYZ (см. рис. 2.5, а) связана с внешним кольцом, причем ось X направле­ на по оси вращения этого кольца. — система, связанная с ос­ нованием. Ось £ этой системы совпадает с осью X вращения внешнего кольца. Основание предполагается подвижным в инер­

циальном пространстве, так что систему

следует рассматри­

вать как инерциальную, а — угол поворота

внешнего кольца

вокруг оси X относительно инерциальной системы |т]£. X'Y'Z'

система, получающаяся в результате поворота системы. ZTZ во­

круг оси Z на угол е (где е — упомянутый

выше перекос оси

вращения внутреннего кольца), xyz — система, связанная с внут­ ренним кольцом (см. рис. 2.5, б). Она получается из системы

21

поворотом на угол р вокруг осп Y'. Углы а и ß определяют поло­

жение колец подвеса в инерциальной системе отсчета

Оси

Ц'і и уі определяются положением оси z x вращения ротора,

но не

участвуют в собственном вращении

ротора вокруг

оси Z\

(оси

Резаля). Система xxy xz x получается

из системы xyz

поворотом

вокруг оси хна угол ум, а затем поворотом вокруг оси у х на угол \ 2 (см. рис. 2.5, в). Углы уі и у2 представляют составляющие пе­

рекоса от оси Z\ ротора относительно внутреннего кольца. Перед тем как переходить к составлению уравнений движе­

ния, сделаем

ряд предположений. Будем предполагать оси х, у,

z и Ху, у I, z 1

главными осями инерции соответственно внутренне­

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

Если на ротор не действуют моменты вокруг осп вращения Z\, то согласно закону о моменте количества движения составля­ ющая абсолютной угловой скорости п ротора по его оси враще­ ния z\ постоянна. Для определения углов а, ß необходимы два уравнения. Чтобы найти одно из этих уравнений, рассмотрим систему, состоящую из внутреннего кольца и ротора. Использу­ ем тот факт, что скорость изменения составляющей момента ко­ личества движения этой системы по оси у вследствие отсутствия моментов вокруг этой оси равна нулю. Основываясь на уравне­ нии

— ф ь Х О = 1 ,

(2.18)'

где G — момент количества движения рассматриваемой

систе­

мы;

 

М — момент, действующий на систему.

 

Разложим это векторное уравнение на составляющие по осям системы xyz, принимаемой в качестве подвижной системы коор­

динат. В уравнении (2.18) со обозначает абсолютную угловую скорость подвижной системы xyz. Принимая во внимание М = 0 и представляя уравнение (2.18) в виде

clG

і

j

k

= 0,

(2.19)

dt

wx

шу

wz

 

 

 

gx

Sy

gz

 

 

где Г, /, к — орты по осям х, у, z, а gx, g v, ^ — проекции кинети­ ческого момента G на эти оси, для составляющей уравнения (2.19) по оси у вращения внутреннего кольца получим

gym-

+ £.cü,*= °.

(2-20)

где со,;, (Oy, cor — проекции угловой скорости системы на оси этой системы.

22

Проекции абсолютной угловой скорости на оси X, у, z имеют вид

(і)г = а cos $ cos 3;

іо = 3 — а sin в; ioz = а cos г sin 3.

внутреннего кольца

(2.21)

Если /Г, В', С —-главные моменты инерции внутреннего коль­ ца, то составляющие момента количества движения этого кольца по указанным осям определяются формулами

g ’r= A 'a cos s cos ß;

gy = B' (ß — а sin s);

( 2.22)

g'z= C'a cos s sin 3.

Проекции угловой скорости системы луг/iZi, связанной с осью Z\ ротора, на оси этой системы имеют вид

(Од.,== a cos г cos у, cos ß + (и — а sin è) X

X sin Yi sin Y2— cx cos £ cos Yi sin Yo sin ß;

totfl= (^ — а sin e) cos Yi + а cos s sin y, sin [3;

(2.23)

(ог1 = а cos s cos Yi cos Y2 sin ß-j-а cos s X

X sin Y2 cos 3 — (3— а sin s) sin Yi cos-Уг-

Зная эти проекции, легко найти составляющие момента коли­ чества движения ротора по осям лу, у и z \

g.vі=

gyi= Amyn gzi = Cn,

(2.24)

где А, С — главные моменты инерции ротора, а п — проекция уг­ ловой скорости ротора на ось zb являющаяся, как уже упомина­ лось, постоянной.

Таким образом, составляющие момента количества движения системы «внутреннее кольцо — ротор» по осям х, у, z имеют вид

 

+

cos y2 + Szi sin y2;

 

 

-

+

sin Yi sin Yo + g'ai cos \i — g sin у, cos у2:

(2.25)

 

gz= gz — gxi cos Y, sin y2-j"gу 1 Sin Yi +

cos yxcos Y;

 

23

Подставляя (2.21) — (2.25) в уравнение (2.20), получим пер­ вое из двух необходимых уравнений движения

[В' (ß — а sin s)-j- A sin Yi sin y2[а cos e cos y2 cos ß-j-

+ (fi — a sin e) sin Yi sin y2— а cos's cos Yi sin y2 sin p]-|-

+ A cos Yi [(P — а sin e) cos Y]+ a cos e sin Yi sin ß]—

Cn. sin Yi cos y2) — а cos s cos ß [C'a cos e sin ß —

A' cos Yi sin y2[a cos e cos y2cos ß — (ß— a sin s) sin Yi sin y2

— a cos s sin ß cos Yi sin y2]+ A sin Yi [(ß— a sin s) cos Yi +

- fa cos e sin Yi sin 8]-\-Cu cos y(. cos y2) -fa cos £ sin ß X

X j A'a cos e cos ß -f A cos y2[a cos s cos y2 cos ß — (ß — a sin e) X

X sin Yi sin y2 —a cos s cos y1sin y2sin ß] -\-Cn sin y2} = 0 . (2.26)

Второе уравнение движения можно получить, рассматривая скорость изменения составляющей момента 'количества движе­ ния всей системы «ротор — внутреннее кольцо — внешнее коль­ цо» по оси вращения внешнего кольца. Если, как предполагалось, трение отсутствует, эта скорость равна нулю. Поскольку ось X является осью инерциальной системы, составляющая момента количества движения всей системы в направлении оси X постоян­ на. Эта составляющая определяется выражением

A"a-|-g-A.cos вcos 3 — g y sin

cos s sin ß = const,

(2.27)

где А"а — момент количества движения внешнего кольца

(А"

момент инерции этого кольца), а gx, g u, gz— составляющие мо­ мента количества движения системы «ротор— внутренее коль­ цо» [см. (2.22) — (2.25)].

Подставляя в (2.26) выражения этих составляющих, получим

Â"a-)-cos s cos ß {А'а cös г cos ß-f-A cos y2[a cos s cos y2 cos ß —

— (ß— a sin г) sin Yi sin y2 — a cos s cos Yi sin y.2sin ß]-j-

-fC /г sin y2) — sin E — a sin s)-j- A sin Yi sin y2X

X [“ cos e cos y2 cos .3 —(ß — a sin e) sin Yi sin y2 —

— a cos s cos Yi sin y2 sin ß]-|- A cos Yi [(ß' —a sin s) cos Yi-f

-j- a cos Esin Yi sin p]— Cti sin Yi cos y2) -f cos e sin ß X

X (C'a cos s sin ß — A cos Yi sin y2 [a cos в cos y2cosß —

— (ß— a sin s) sin Yi sin y2 — a cos s cos Yi sin y2 sin ß]-j-

-j- A sin Yi [(ß— a sin s) cos Yi + a cos в sin Yi sin ß]-j-

 

-f Cn cos Yi cos y2] = const.

(2.28)

24

Если через gx обозначить левую часть

(2.28), то искомое

уравнение движения примет вид

 

- ^ = 0.

(2.29)

dt

 

Уравнения движения (2.26), (2.28) представляют собой сис­ тему двух нелинейных дифференциальных уравнений с .неизвест­ ными а и ß. Точного метода решения таких уравнений не сущест­ вует. Однако решение можно получить в виде функционального ряда по степеням малого параметра. Опишем этот метод реше­ ния для начальных условий

t = 0, а = а0, а = й,

(2.30)

?= &>, £= 0.

Эти начальные условия получаются при сообщении внешнему

кольцу в момент / = 0 угловой скорости а = й, что соответстңует нанесению удара по внешнему кольцу в момент ^ = 0.

При небольшом значении начальной угловой скорости Q ре­ шение можно искать в виде функционального ряда по степеням малого параметра й. Следуя этому методу, можно записать

a = a°+a /2 +

a',Ss+ . . . ,

(2.31)

ß = ß" + ß'Q +

ß"£2+ • • (

2 .

3 2 )

где а0, а', а", ... ß°, ß', ß"— неизвестные функции времени. Эти функции можно определить, подставляя ряды (2.31), (2.32) в уравнения движения (2.26), (2.28). Членами, содержащими Й в степени выше второй, можно пренебречь. Учет членов с более высокими степенями приводит к усложнению выкладок. Вместе с тем основные свойства решения выявляются при учете лишь членов до второго порядка малости.

При использовании рядов (2.31), (2.32) составляются при­ ближенные выражения для sin ß, cos ß, sin ß cos ß:

sin ß ~ sin (ß°-f ß'Q + ^ 22)= sin ß° cos (ß'2 -f ß"S2)-{-

+ cosß0sin(P/Q-f ß"22) ~ s in ß o |l -- L ß '2g^ -)-

+

cosß°(ß'2 +

ß"22);

'

(2.33) '

cos В~ cos ß° ^ 1

- -L ß'2Q2j _

sin ßo (ß'g + ß"Q2).

(2.34)

Приближенные выражения для sin ß cos ß и cos2 ß можно по­ лучить, перемножая (2.33), (2.34) и возводя в квадрат (2.34). При этом в окончательных выражениях должны быть сохранены только члены с параметром й в степени не выше второй. Эти выражения, а также принятая форма решения (2.31), (2.32) подставляются в уравнения движения (2.26), (2.28). Левые час-, ти уравнений превращаются в степенные ряды в отношении па­

25

раметра й. Чтобы уравнения удовлетворялись, коэффициенты при различных степенях й, зависящие от а и ß, далжны равнять­ ся нулю. Принимая во внимание в каждом пз уравнений лишь члены с нулевой степенью Й, получим два дифференциальных уравнения, переменными в которых являются а0 и ß°. Выбирая для а0 II ß° начальные условия

(2.35)

совместимые с начальными условиями (2.30) для а, ß и приня­ той формой решения (2.31), (2.32), путем решения этих уравне­ ний получим а0, ß° как функции времени.

Принимая далее во внимание члены с первой степенью й, по­ лучим два дифференциальных уравнения в отношении перемен­ ных a', ß'. В этих уравнениях помимо других фигурируют также члены, содержащие теперь уже известные функции а0, ß°. Пере­ носим эти члены в правые части уравнений. Решая эти уравне­ ния при начальных условиях

^ = 0, а' = 0,

а '= 1 ,

ß'= 0,

ß' = 0,

(2.36)

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

условиями

(2.30), (2.35),

(2.36) и

принятой формой решения

(2.31),

(2.32),

получим а'

и ß' как

функции времени.

Принимая во внимание члены со второй степенью й, получим два дифференциальных уравнения (их правые части содержат члены с теперь уже известными а' и ß'), решая которые при ну­ левых начальных условиях, находим а", ß" как функции време­ ни. В данном случае нулевые начальные условия совместимы с начальными условиями (2.30), (2.35), (2.36) и принятой формой решения (2.31), (2.32).

Подставляя, наконец, найденные значения а0, а', а", ß°, ß', ß" в (2.31), (2.32), получим искомые приближенные решения урав­ нений движения (2.26), (2.28).

Из-за сложности общих уравнений движения (2.26), (2.28) наиболее целесообразно установить отдельно влияние каждого из указанных параметров карданова подвеса на дрейф гироскопа.

Если положить -уI = Y2— 0 и учитывать только перекос е оси

вращения внутреннего кольца, в этом случае уравнения движе­ ния принимают вид

/г3—(/, sin s) а — {Си cos s) а cos ß-f-

-\-[{А-\- А' — C') cos2 е] а- sin 3 cos ß — 0;

— (/2sin s)ß-|-(Cß cos s)ß cos ß-j-

+ [Л"+ + А' - C') cos2 e cos2 ß +

-j- C' cos2 s — /, sin2 e] а —

— 2 [(Л + А' C') cos2 e] aß sin ß cos ß = 0,

где і 2 = А+В'.

26

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