и любой другой точке плоскости Oxy , т.е. равны нулю два центробежных момента инерции Jxz=Jyz=0. В связи с этим можно утверждать, что центробежные моменты инерции (вместе со статическими моментами массы тела) составляют меру материальной несимметрии тела относительно координатных плоскостей.

oЕсли тело имеет ось материальной симметрии, (тело вращения) то для любой точки О, взятой на оси симметрии, главными осями инерции являются ось симметрии Oz и любые два перпендикуляра к

этой оси Ox, Oy, т.е. Jxy=Jxz= Jyz=0 Таким образом, главные оси инерции точно изготовленного волчка, маховика, снаряда заранее известны для их центров масс и других точек, взятых на оси материальной симметрии.

3.6.Перепроецирование вектора с применением матрицы поворота

Рассмотрим две системы координат:

o исходную n-систему Oxnynzn с ортами in, jn, kn

o s-систему с ортами is, js, ks, повернутую от исходной

zs

 

 

zn

r

M

 

 

 

 

 

kn

 

 

ys

 

ks

 

 

 

 

 

js

 

 

 

 

 

 

 

 

O

 

 

 

 

in

 

 

 

 

yn

is

jn

 

xn

 

 

 

 

 

 

 

 

 

 

 

 

 

xs

 

 

 

Рис. 5 Радиус-вектор точки М в исходной и повернутой системах отсчета

Радиус-вектор r любой точки M можно разложить на составляющие, выразить через орты и координаты точки в той и другой системе отсчета:

 

r = xs

 

+ ys

 

 

 

 

 

 

 

+ yn

 

 

 

 

 

 

 

 

is

js + zs ks = xn

in

jn + zn kn

 

 

 

Умножим скалярно это равенство на орт

 

 

 

 

 

. Аналогичную операцию

 

is

 

 

выполним с ортом

 

 

 

 

, затем с ортом

 

 

 

 

. В результате получаем три мат-

 

 

js

 

ks

ричных выражения:

i i

 

 

 

 

 

 

 

 

 

 

i j

 

 

 

 

 

 

 

 

 

 

 

i k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

=[x y

z

n

]

jn

is

 

y

s

=[x y

z

n

] jn

js

 

z

s

=[x y

z

n

]

jn

ks

 

s

n n

 

 

 

n

s

 

 

n n

 

 

n

 

s

 

 

 

 

n n

 

 

 

n

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

js

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

kn

is

 

 

 

 

 

 

 

 

 

kn

 

 

 

 

 

 

 

 

 

 

 

kn

ks

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

35

Объединим три выражения в одну вектор-строку, т.е. выполним горизонтальную конкатенацию (сцепление) трех выражений. Получим одно матричное выражение:

rs = rnCsn при rs =[xs ys zs ], rn =[xn yn zn ]

(3.16)

in iS

Csn = jn iS

kn iS

in jS jn jS kn jS

in

kS

 

in

 

(3.17)

 

 

 

 

 

 

[is js ks ]

jn kS

jn

 

k

n

k

 

 

k

 

 

 

 

 

S

 

 

n

 

 

Здесь rs и rn – вектор-строки вектора r в двух системах отсчета, Csn - матрица поворота s-системы от n-системы отсчета.

Скалярные произведения ортов по существу являются направляющими косинусами осей s -системы по отношению к осям n-системы, например in js = cos(xn , ys ) . В выражении (3.17) формально введено понятие скалярно-

го произведения матрицы-столбца ортов на матрицу-строку ортов. Выражения (3.16) и (3.17) можно запомнить с помощью следующего мнемонического правила. В формуле перепроецирования вектора r из n-системы в s- систему координат (3.16) индексы повторяются в неизменном порядке, а в матрице косинусов (3.17)- в обратном порядке, т.е. орты со вторым индексом выстроены в столбец, а с первым индексом - в строку. Матрица Csn обладает

 

1

свойством

Csn

=Csn =Cns , т.е. обратная матрица равна транспонированной

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

Пример Пусть в роли повернутой s-системы выступает “нулевая” система Oxyz ,а в роли исходной системы – “первая” система Ox1y1z1. Тогда на основании общей формулы (3.16) получим:

r = rC ,

r = rC ,

(3.18)

1

01

1

10

 

где r =[x y z], r1 =[x1 y1 z1 ], C01 =[i1

j1 k1 ]' [i j k], C10 =[i j k]' [i1

j1 k1 ]

C01 - матрица поворота нулевой системы от первой, C10 -первой системы от

нулевой.

Замечание. Не только радиус-вектор r , но и любой вектор a перепроецируется из системы Ox1y1z1 на систему Oxyz и обратно по матричным формулам

a1 = aC10 = a C01, a = a1C01 при a1 =[ax1 ay1 az1] , a =[ax ay az ] (3.19)

Эти формулы перепроецирования, содержащие первые степени направляющих косинусов, выражают основное тензорное свойство любого вектора a. В связи с этим любой вектор называют также тензором первого ранга в пространстве трех измерений.

3.7.Тензор инерции в точке твердого тела

36

Рассмотрим твердое тело. Возьмем декартову систему координат Охуz c ортами i , j , k , скрепленную с телом. Из осевых и центробежных моментов

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

J x

JO = J yx

J zx

J xy

J xz

 

(3.20)

J

y

J

 

, где Jxy= Jyx , Jxz=Jzx, Jyz=Jzy = −∫∫∫ ρyzdV (кг*м2)

 

 

yz

 

 

J zy

J z

 

 

 

 

 

 

 

 

Матрицу (3.20) называют матрицей тензора инерции тела в точке О,

или кратко — тензором инерции тела. Величины элементов матрицы (3.20) зависят от направления осей в теле, изменяются при изменении этих направлений и зависят от выбора полюса О.

Замечание. Во многих учебниках не вводится знак минус в определяющих формулах (3.14), но при этом входят со знаком минус все центробежные моменты инерции в матрице (3.20). Данное отличие не является существенными, но предпочтительнее иметь в матрице только положительные знаки.

Свойство тензора инерции, тензора второго ранга в пространстве трех измерений, излагается в следующей теореме.

Теорема1: Матрица инерции JО1=J1, отнесенная к декартовой системе Ох1у1z1 определяется через исходную матрицу JО , заданную в системе Oxyz , перемножением трех матриц: транспонированной матрицы поворота “первой” системы от исходной С10= C01´, матрицы JO и матрицы поворота. А именно,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3.21)

J

 

= C '

J

 

C

при

 

 

 

 

 

 

[i

 

j

k ]

1

O

C =[i j k ]

 

 

 

10

 

10

 

10

 

 

 

 

 

 

1

1

1

 

И наоборот, тензор инерции заданный в системе Ох1у1z1 , проецируется на систему Oxyz согласно формуле: JO = C01´ J1 C01

Формула (3.21) напоминает формулу перепроецирования вектора a1=a C10 , но отличается добавлением множителя C10´ = C01 , поставленного на первом месте.

Теорема 2 Для каждого полюса О тела существует такая декартова система координат Оxyz, в которой матрица инерции J0 имеет диагональный вид

JO =diag [ Jx , J y , Jz ] diag [ A, B, C ]

Диагональные элементы A, B, C этой матрицы являются главными мо-

ментами инерции в точке О тела, а оси Оx, Оy, Оz - главными осями инерции. Здесь все элементы вне главной диагонали, т.е. все центробежные моменты инерции, равны нулю.

Главные моменты инерции Jx=А, Jy=В, Jzможно находить как три корня характеристического уравнения det (λE - J1 ) = 0 , где Е — единичная

диагональная матрица третьего порядка, J1 -матрица, вычисленная в исходной системе Ox1 y1z1 . В пакете MATLAB они находятся как три собственных

действительных значения квадратной симметрической матрицы третьего порядка J1 .

37

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

Определение Тензором инерции в точке О тела называется физическая величина, задаваемая матрицей вида (3.20) в конкретно выбранной системе Охуz и вычисляемая в любой повернутой системе по формуле (3.21).

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

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

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

полуосями, равными 1/ A , 1/ B , 1/ C . Положение в теле такого трехпараметрического эллипсоида в окрестности точки O определяется тремя угловыми параметрами. Таким образом, эллипсоид в целом определяется шестью параметрами, также как и матрица тензора инерции.

Геометрическое свойство эллипсоида инерции, построенного для выбранной точки O тела: Расстояние ОМ от точки O до любой точки М эллипсоида инерции обратно пропорционально корню квадратному из момента инерции тела относительно линии ОМ.

Момент инерции Je относительно произвольной оси ОL, заданной ортом e в системе Оxyz, можно находить через тензор инерции Jo по формуле

Jе = е Jo е´,

(3.22)

где Jo имеет вид (3.20), e =[ex ey ez ] и e′=[ex ey ez ]T

есть строчная и

столбцевая матрицы прекции орта, называемых направляющими косинусами оси ОL. В результате перемножения строки e на (3х3)-матрицу Jo и далее ум-

ножения полученной строки на стоблец eполучаем искомую скалярную ве-

личину — осевой момент инерции Jе. Выражение (3.22) в развернутом виде:

Je = Jxex2+Jyey2+Jzez2+2Jxyexey+2Jxzexez+2Jyzeyez.

(3.23)

В качестве орта e можно, например, выбрать орты осей повернутой от Охуz системы О1х1у1z1 . Последовательно принимая в качестве e орты i1 , j1 и k1 ,

получаем по формуле (3.23) выражения осевых моментов инерции J x1 , J y1, J z1. Напомним, что по общей матричной формуле (3.22) находятся

не только осевые, но и центробежные моменты инерции J x1 y1 , J x1 z1, J y1 z1. .

Все эти выражения возвращаются, раскрываются автоматически, в системах MATLAB, Maple и др. в символьном или численном виде.

Рассмотрим тело массы m в двух системах координат Cхуz , О1х1у1z1 со взаимно параллельными осями, причем пусть система Cхуz имеет начало в

38