
- •Учебное пособие
- •Предисловие
- •Список основных обозначений
- •Выдающиеся ученые, создатели основ механики
- •ВВЕДЕНИЕ
- •1. МАТЕМАТИЧЕСКИЕ СВЕДЕНИЯ
- •1.1. Матрицы и векторы
- •1.2. Конкатенация матриц и линейных функций
- •1.3. Перепроецирование векторов и тензоров второго ранга
- •1.4. Решение системы алгебраических линейных уравнений MATLAB
- •2. ДИНАМИКА МАТЕРИАЛЬНОЙ ТОЧКИ
- •2.1. Инерциальные системы отсчета
- •2.2. Поступательные (невращающиеся) система отсчета
- •2.4. Пример
- •2.5. Второй закон динамики в кинетостатической форме Даламбера
- •2.6. Закон динамики относительного движения
- •2.7. Уравнения относительного движения в форме Даламбера
- •2.8. Основной закон динамики относительного движения
- •2.9. Случай поступательного движения системы отсчета
- •2.10. Случай поступательной системы отсчета с началом на поверхности Земли
- •2.11. Прямолинейное движение материальной точки во вращающейся системе отсчета
- •3. ДИНАМИКА ТВЕРДОГО ТЕЛА
- •3.1. Геометрия масс твердого тела
- •3.2. Статические моменты массы и координаты центра масс тела
- •3.3. Осевые моменты инерции и радиусы инерции тела
- •3.4. Осевые моменты инерции площади поперечного сечения тела
- •3.5. Центробежные моменты инерции твердого тела
- •3.6. Перепроецирование вектора с применением матрицы поворота
- •3.7. Тензор инерции в точке твердого тела
- •3.8. Угловая скорость сферического движения тела
- •3.9. Кинетический момент сферического движения тела
- •3.10. Динамические уравнения Эйлера
- •3.11. Элементарная теория гироскопа
- •3.12. Динамика произвольного движения твердого тела
- •4.1. Сведения из статики
- •4.2. Момент силы и пары сил, расположенных в плоскости
- •4.3 Бивектор плоской системы сил
- •4.4. Уравнения равновесия плоской системы сил
- •4.5 Бивектор и M-функции пространственной системы сил
- •4.6. Матричное уравнение равновесия пространственной системы сил
- •4.7. Бивекторы реакций опор
- •4.8. Пример конкатенации матриц в матричном уравнении равновесия
- •4.9. Кинетостатическое матричное уравнение произвольного движения тела
- •5. ДИНАМИКА ГОЛОНОМНЫХ МЕХАНИЧЕСКИХ СИСТЕМ
- •5.1. Фазовый вектор системы (вектор состояния системы)
- •5.2. Кинетическая энергия твердого тела и голономной стационарной одностепенной системы
- •5.3. Мощность силы
- •5.4. Работа и потенциальная энергия
- •5.5. Обобщенная сила одностепенной голономной стационарной системы
- •5.6. Обобщенная сила нестационарной одностепенной системы
- •5.7. Работа обобщенной силы одностепенной стационарной системы
- •5.8. Идеальные связи и реакции связей
- •5.9. Уравнение Лагранжа для одностепенной стационарной голономной системы
- •5.10. Уравнение Лагранжа для голономной нестационарной одностепенной системы
- •6. ОБЩИЕ ТЕОРЕМЫ ДИНАМИКИ МЕХАНИЧЕСКОЙ СИСТЕМЫ
- •6.1. Уравнения Лагранжа для голономной стационарной системы
- •6.2.Уравнения равновесия стационарной системы в обобщенных координатах
- •6.3. Матричная форма уравнений Лагранжа
- •6.4. Матричное уравнение Гамильтона
- •6.5. Теорема о движении центра масс механической системы
- •6.6. Теорема об изменении количества движения механической системы
- •6.7. Теорема импульсов
- •6.8. Теорема об изменении кинетической энергии
- •6.9. Теорема об изменении кинетического момента механической системы
- •Приложение
- •Равновесие плоской системы сил
- •Равновесие пространственной системы сил
- •Динамика относительного движения материальной точки
- •Литература

и любой другой точке плоскости 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