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

Трухан. Динамика твердого тела

.pdf
Скачиваний:
50
Добавлен:
03.06.2015
Размер:
733.35 Кб
Скачать
OXY,

11

центральные моменты инерции A, B, C. Построить главные

оси

инерции

тела

для

точки O

с координатами

ξO = a, ηO = b, ζO = O.

 

 

 

Решение. В точке O построим систему координат

OXYZ с осями,

параллельными главным центральным осям.

Так

как точка

O

лежит

в плоскости,

содержащей две

главные центральные оси, то по выводу 4 из задачи 1.1 ось OZ - главная для точки O. Значит нужно найти только две главные оси.

Пусть некоторая ось u лежит в плоскости проходит через точку O и образует с осью OX угол α. Тогда направляющие косинусы оси u в системе OXYZ будут равны соответственно cosα, sinα, 0. Момент инерции тела относительно оси u приобретает вид

Ju = J x cos2 α + J y sin2 α 2J xy sinα cosα. (1.19)

Так как в данном случае

cos2 α +cos2 β +cos2 γ = cos2 α +sin2 α =1

независимо от значения угла α, то для отыскания

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

можно найти из условия стационарности функции Ju :

 

 

dJu

= 0.

 

(1.20)

Что дает

 

dα

 

 

 

 

 

 

 

 

 

 

 

 

 

J x sin 2α + J y sin 2α 2J xy cos 2α = 0.

Откуда

1

 

 

2Jxy

 

α =

arctg

.

2

J y

Jx

 

 

 

 

Но из задачи 1.1 следует

J x = A +mb2 , J y = B + ma2 , Jxy = mab.

 

 

 

12

 

 

И, значит,

1

 

2mab

 

α =

 

 

 

arctg

B A +m(a2 b2 )

.

(1.21)

2

Таким образом, повернув ось OX вокруг точки O на угол α, определяемый равенством (1.21), мы совместим ее с

главной осью. Момент инерции относительно построенной главной оси может быть найден в результате подстановки α из равенства (1.21) в выражение (1.19).

§ 2. Подсчет основных динамических величин для твердого тела с одной неподвижной точкой

2.1. Кинетический момент. В каждый момент времени твердое тело с неподвижной точкой О имеет такое

же распределение скоростей Vi точек, как при вращении

вокруг неподвижной оси, т.е.

 

 

 

 

 

 

 

 

 

 

Vi

=

ω

×

ri

, где

ω

- вектор

мгновенной угловой скорости,

ri

- радиус-вектор i -й точки,

проведенный из неподвижной точки. Поэтому кинетический момент тела относительно точки О будет

KO = miri ×Vi = miri ×(ω ×ri )

или с помощью тензора инерции J KO = Jω.

Если тензор инерции J записан в системе осей, главных для точки О, то он имеет диагональный вид. Вектор угловой

скорости в этих осях

 

= pi

 

 

 

 

 

 

 

 

 

ω

 

+qj

+rk ,

и вектор

кинетического момента

 

 

 

O = Api

 

+ Bqj

 

+Crk

 

 

(2.1)

 

K

 

 

 

( A, B, C - осевые моменты инерции). Векторы ω и KO , вообще говоря, не коллинеарны. Они, очевидно, совпадают

13

по направлению лишь в том случае, если ось вращения совпадает с главной осью. При этом KO = Jω = λω , где λ является коэффициентом пропорциональности. И, значит, J λEω = 0. Условие нетривиальности решения этого

уравнения: detJik λδik = 0. Отсюда следует, что

собственные векторы матрицы J определяют направления главных осей инерции, а собственные значения этой матрицы

моменты инерции относительно главных осей.

2.2. Кинетическая энергия.

Кинетическая энергия

T =

1

miVi

2

=

1

mihi

2ω2 =

1

JOωω2. (2.2)

2

 

2

2

 

 

 

 

 

 

 

 

Здесь JOω = mihi2 - момент инерции тела относительно мгновенной оси вращения, hi - расстояние i-й точки от этой оси.

По теореме Гюйгенса-Штейнера JOω = JM ω +md 2

(d - расстояние между мгновенной осью вращения и осью, проходящей через центр масс параллельно мгновенной,

JM ω - момент инерции тела относительно последней).

Выражение кинетической энергии можно представить в форме

T = 12 md 2ω2 + 12 JM ωω2 = 12 mVM 2 + 12 J M ωω2 , (2.2)

где VM - скорость центра масс тела. Формула (2.2)

представляет собой аналитическое выражение теоремы Кенига для твердого тела. Обратим внимание, что кинетическая энергия относительного движения

T '=

1

JM ωω2

есть кинетическая энергия вращательного

2

 

 

 

14

движения вокруг оси, проходящей через центр масс

параллельно мгновенной.

 

 

 

 

 

 

 

Выражение

(2.2)

кинетической

энергии

в

осях

Oξηζ , главных для

неподвижной

точки,

записывается в

виде

 

 

1

(Ap2 + Bq2

+Cr 2 ).

 

 

 

 

 

T =

 

 

(2.3)

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

Аналогично формула (2.2) в главных центральных

осях может быть представлена так:

 

 

 

2 ).

 

T =

1

mVM

2

+

1

(AM pM

2 + BM qM

2

+CM rM

(2.4)

 

 

 

 

2

 

 

2

 

 

 

 

 

 

 

Задача 2.1. Однородный круговой конус массы m с

высотой h и углом при вершине 60D катится без проскальзывания по неподвижному конусу с углом при

вершине 120D (см. рис. 2.1). Найти величину вектора кинетического момента KO и угол β, который образует этот

вектор с вектором ω. Подсчитать кинетическую энергию T конуса, если центр его основания движется с постоянной по величине скоростью V.

Решение. Величину вектора KO найдем с помощью равенства

KO =

O 60o O1

120o

A2 p2 + B2q2 +C 2r 2 . (2.5)

Так как центральный эллипсоид инерции эллипсоид вращения, а точка O лежит на главной центральной

оси, то главные оси для этой точки получаются

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

Рис. 2.1

 

 

 

 

 

15

 

 

 

 

 

 

 

 

центральных осей MξMηMζM .

Центр масс однородного

кругового конуса лежит на расстоянии

 

3

h от его вершины,

4

 

3

 

 

 

 

 

 

 

 

 

т.е. OM =

h. Главные

центральные моменты

инерции

4

 

 

 

 

 

 

 

 

 

 

 

 

конуса равны

 

 

3

 

mR2 , где R = h

 

AM = BM =

 

3

m(4R2 +h2 ), CM

=

 

-

80

 

 

 

 

10

 

 

 

 

3

 

радиус

основания

конуса.

 

Таким

образом,

A = B = AM

_

ω

ζ

+mOM 2

=

 

13

mh2

,

C =CM =

 

 

1

mh2.

Вектор

20

10

 

 

 

 

 

 

 

 

ω

 

_

 

 

 

 

 

угловой

 

скорости

 

 

 

 

 

 

 

направлен

 

по

Ko η

 

 

 

 

 

 

 

 

 

ηM

 

мгновенной

 

оси

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

вращения,

 

 

ζM M

 

совпадающей с линией

O

 

 

 

 

 

касания

подвижного и

 

 

 

 

 

неподвижного конусов.

ξ

ξM

 

Центр

 

 

 

основания

 

 

 

 

 

 

конуса

 

 

движется

со

Рис. 2.2

 

 

 

 

 

скоростью

V

 

и

 

 

 

 

 

находится

 

 

от

 

 

 

 

 

 

 

 

мгновенной оси на расстоянии h2 , поэтому угловая скорость

ω =

V

=

2V

. Проекции вектора

ω

на оси системы Oξηζ :

h 2

h

 

 

 

 

 

 

 

 

 

p = 0,

q =

V

,

r =

3 V

. Подставляя в (2.5) найденные

 

 

 

 

 

 

h

 

 

h

 

 

 

значения моментов инерции и компонентов вектора угловой

скорости, получим KO =

181 mVh. Вектор

 

O образует с

K

 

20

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

16

 

 

 

 

 

 

 

 

осью

 

Oζ

 

угол

 

 

α = arctg

Bq

= arctg

13 .

Откуда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Cr

 

 

 

2 3

 

β =

30D α

.

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для подсчета

воспользуемся

формулой (2.3).

Используя найденные при подсчете

KO значения

A, B, C и

p, q, r, получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

13

 

 

2 V

2

 

1

 

2 3V

2

 

 

19

 

2

 

 

T =

 

mh

 

 

mh

 

 

 

mV

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

=

 

 

.

 

2

20

 

h

2

10

 

h

2

 

 

40

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

§ 3. Уравнения движения тела с неподвижной точкой. Случай Эйлера и случай Лагранжа

3.1. Уравнения движения. Твердое тело с одной неподвижной точкой имеет три степени свободы.

Часто в качестве координат, однозначно задающих положение твердого тела с неподвижной точкой, берут углы Эйлера: ψ - угол

ζ

Z

 

прецессии,

 

θ

-

угол

_.

_._

η

нутации,

ϕ

-

угол

ϕ

θ ψ

 

собственного

вращения

 

_

Y

(см. рис. 3.1). Эти углы

 

и

их

производные

ψ

θ.

 

связаны с

 

проекциями

ϕ

 

вектора

 

 

угловой

X

ξ

 

скорости

ω

на главные

 

оси

кинематическими

 

 

 

Рис. 3.1

 

соотношениями Эйлера:

 

p =ψ sinθ sin ϕ+θ cosϕ

 

 

 

 

 

 

 

q =ψ sinθ cos ϕ−θ sin ϕ

 

 

 

 

 

(3.1)

r =ψ cosθ .

17

Уравнения динамики тела с одной неподвижной точкой получим с помощью теоремы об изменении кинетического момента

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

K

O

=

 

 

 

 

 

 

 

O.

(3.2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Запишем вектор

 

 

 

O в проекциях на

главные

оси,

K

построенные для неподвижной точки

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

O = Api

 

+ Bqj

 

 

+Crk

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как система главных осей подвижная и жестко

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dk

=

ω

 

 

 

 

 

 

 

 

 

di

=

,

 

 

dj

 

=

 

×

 

,

 

 

 

 

связана с телом, то

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

×k

.

 

 

 

ω

×i

 

 

 

ω

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d'

 

 

 

 

dt

 

 

dt

 

dt

 

 

 

Обозначим

K

O

 

 

 

 

 

 

 

 

 

 

 

= Api

+ Bqj +Crk

. Тогда

 

теорему

об

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

изменении кинетического момента можно записать в форме

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dK

O

=

d' KO

+

 

×

 

O =

 

 

O .

(3.3)

 

 

 

 

 

ω

K

M

 

 

 

 

 

dt

 

 

 

 

 

 

 

dt

 

 

d'

 

O

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь

- изменение вектора KO

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

dt

выбранной подвижной системы координат (иначе относительная скорость конца вектора KO ), а ω ×KO -

изменение вектора KO за счет движения подвижной системы координат (переносная скорость конца вектора KO ). Таким образом, движение конца вектора KO рассматривается как

сложное: вместе с подвижной системой координат, жестко связанной с телом и относительно нее.

Проецируя равенство (3.3) на главные оси, построенные для неподвижной точки, получаем динамические уравнения Эйлера:

Ap +(C B)qr = Mξ

18

 

Bq +(A C)pr = Mη

(3.4)

Cr +(B A)qp = Mζ .

 

Уравнения (3.4) совместно с уравнениями (3.1) описывают движение тела с неподвижной точкой и записываются в подвижных осях, жестко связанных с телом, главных для неподвижной точки.

Общее решение системы (3.4) и (3.1) позволяет в каждый момент найти значения углов Эйлера, однозначно задающих положение тела с неподвижной точкой, т.е.

ψ =ψ(t,C1,...,C6 ),

θ =θ(t,C1,...,C6 ), (3.5) ϕ = ϕ(t,C1 ,...,C6 ).

Известно, что лишь в трех случаях при произвольных начальных условиях интегрирование системы уравнений Эйлера может быть доведено до квадратур. Это случаи Эйлера, Лагранжа и Ковалевской. Здесь будет рассмотрено решение задач в случаях Эйлера и Лагранжа.

3.2. Случай Эйлера. Этот случай соответствует движению по инерции, когда MO = 0.

Тогда при произвольной форме тела сохраняется вектор кинетического момента KO и кинетическая энергия

T. Интегрирование системы динамических уравнений Эйлера может быть доведено до квадратур, которые получаются в виде эллиптических интегралов. При некоторых начальных условиях эти квадратуры вырождаются в обыкновенные интегралы, и тогда p, q, r

являются гиперболическими функциями времени.

Задача решена до конца, если найдена зависимость эйлеровых углов от времени.

Так как KO = const, возьмем неподвижную ось OZ

вдоль вектора кинетического момента. Тогда

19

Ap = KO sinθ sin ϕ,

Bq = KO sinθ cosϕ,

Cr = KO cosθ.

Из (3.6) находим

Cr(t),

cosθ =

 

KO

tg ϕ =

Ap(t)

Bq(t).

(3.6)

(3.7)

(3.8)

Чтобы найти ψ, умножим

первое из уравнений (3.1) на

sin ϕ, а второе на cos ϕ и сложим

 

 

 

 

 

 

psin ϕ+ q cosϕ =ψ sinθ.

 

Откуда

 

 

 

psin ϕ+ q cosϕ

 

 

 

 

ψ =

.

 

 

 

 

 

 

 

 

 

 

sinθ

 

 

 

 

 

Используя (3.7) и (3.8), найдем выражения

 

 

sinθ =

1

A2 p2 + B2q2 ,

 

 

Ap

 

KO

 

 

 

 

Bq

 

sin ϕ =

 

 

cosϕ =

 

 

A2 p2 + B2q2 ,

 

A2 p2 + B2q2 .

 

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

ψ =

 

 

Ap2

+ Bq2

K

O

(3.9)

 

 

A2 p2

+ B2q2

 

 

 

 

 

и

 

 

 

 

 

Ap2 (t)+ Bq2 (t)

 

 

 

 

 

 

t

 

 

ψ =ψO + KO 0

 

dt. (3.10)

 

A2 p2 (t)+ B2q2 (t)

Таким образом, углы Эйлера θ, ϕ и ψ могут быть найдены

из равенств (3.7), (3.8) и (3.10), если определены p(t), q(t), r(t).

20

Задача 3.1. Однородный эллипсоид с моментами инерции A, 32 A, 2A движется вокруг своего неподвижно

закрепленного центра масс. В начальный момент эллипсоиду сообщена угловая скорость Ω =ω 2i +ωk , где i , j, k -

единичные орты главных центральных осей инерции. Получить уравнения движения эллипсоида в квадратурах.

Решение. Динамические уравнения Эйлера в этом

случае

 

 

 

 

 

 

 

dp

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

+ qr = 0,

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3 dq

 

2 pr = 0,

 

(3.11)

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4 dr

+qp = 0.

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

Умножив первое уравнение на p, третье на r

и приравняв

затем левые части этих уравнений, получим

 

 

p

dp

= 2r

 

dr

 

 

p2 p

2

= r2 r

2

.

 

 

 

 

 

O

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

dt

 

 

 

 

 

2

 

O

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как pO = rO

2, то p = r

 

2 в любой момент времени.

Разделив первое уравнение (3.11) на второе, получим

уравнение с разделяющимися переменными

 

 

 

 

 

 

 

 

4 p dp = −3q dq.

 

 

Откуда

p2 =

1

(4 pO

2

3q2 ),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

3q2 ).

 

 

 

 

и, кроме того,

r 2

=

1

(4 pO

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Связь между p, q и r можно было получить и

непосредственно, используя интеграл энергии и закон сохранения кинетического момента: