
- •Учебное пособие
- •Предисловие
- •Список основных обозначений
- •Выдающиеся ученые, создатели основ механики
- •ВВЕДЕНИЕ
- •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. Теорема об изменении кинетического момента механической системы
- •Приложение
- •Равновесие плоской системы сил
- •Равновесие пространственной системы сил
- •Динамика относительного движения материальной точки
- •Литература

1.МАТЕМАТИЧЕСКИЕ СВЕДЕНИЯ
1.1.Матрицы и векторы
Матрицей называется квадратная или прямоугольная таблица, запол-
ненная числами, либо - постоянными или переменными символами.
В основном мы будем рассматривать квадратные матрицы и
матрицы-строки |
вида |
[a,b,c,d ] =[a b c d ], а также матрицы-столбцы, |
|
T |
|
′ |
|
{a b c d} =[a b c d] |
|
=[a b c d] , получаемые в результате операции транспо- |
нирования матриц-строк и отмечаемые фигурными скобками. Например, ра- диус-вектору r (x, y, z) и любым другим трехкомпонентным векторам, имеющим геометрический образ и правило проектирования на различные
оси, сопоставляется вектор-строка (матрица-строка) |
r = [x y z] |
или вектор- |
||||||||||||||
|
|
|
|
|
′ |
|
|
|
|
|
|
|
|
|
|
|
столбец {r} ={x y z} = [x y z] проекций. |
|
|
|
|
|
|
|
|||||||||
|
Скалярное произведение двух векторов a и |
|
|
вычисляется как произ- |
||||||||||||
|
b |
|||||||||||||||
ведение строки a на столбец {b} =bT : |
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
c = a |
|
= a{b} = ax bx |
+ ay by + az bz . |
|
|
|
|
|
|
(1.2) |
|||
|
|
|
b |
|
|
|
|
|
|
|||||||
|
Вектор-строки векторного произведения c = a × |
|
= − |
|
×a |
|
|
|||||||||
|
b |
b |
|
|
||||||||||||
определяется формулой |
% |
|
|
|
|
|
|
|
|
|
(1.3) |
|||||
|
|
|
% |
|
|
|
|
|
|
|
|
|
||||
|
~ |
|
|
|
c =ba = −ab . |
|
|
|
|
|
|
|
|
|
|
|
Здесь a — кососимметрическая матрица третьего порядка, вида |
|
|
||||||||||||||
~ |
0 |
az |
− ay |
|
|
|
|
|
|
|
|
|
|
(1.4) |
||
|
0 |
|
− ay ; |
− az 0 ax ; ay − ax |
0] |
|||||||||||
a |
= − az |
ax = [0 az |
|
|||||||||||||
|
|
− ax |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
ay |
0 |
|
|
|
|
|
|
|
|
|
|
|
Итак, вектор-строка векторного произведения равна строке b второго сомножителя, умноженной на кососимметрическую матрицу a% первого сомножителя, либо взятому со знаком минус произведению строки a первого со-
множителя и кососимметрической матрицы b% второго сомножителя. В MATLAB вектор-строка c векторного произведения a ×b возвращается функ-
цией cross вида c=cross(a,b).
В результате транспонирования выражений (1.3) и учета свойства ко-
~ |
~ |
находим выражения |
для векторного |
сосимметрической матрицы a |
′ = −a , |
||
произведения в виде вектор-столбца |
~ |
|
|
|
~ |
(1.5) |
|
{c} = −a{b} = b{a}. |
Квадрат кососимметрической матрицы есть симметрическая матрица
7

% |
|
|
−a2y −az2 |
axay |
axaz |
|
|
(1.6) |
|||||
2 |
|
ayax |
2 |
|
2 |
ay az |
|
= |
|||||
a |
|
= |
−ax |
−az |
|
|
|||||||
|
|
|
a |
a |
x |
a |
a |
y |
−a2 |
−a2 |
|
|
|
|
|
|
z |
|
z |
|
x |
y |
|
|
|
=[−a2y −az2 ,axay ,axaz ;ay ax , −ax2 −az2 ,ayaz ;az ax ,az ay , −ax2 −a2y ]
Вэтих выражениях применен формат записи матриц в пакете МatLab, где строки отделяются знаком (;), а элементы в строке отделены пробелами или запятыми.
Библиотека M-функций в механике.
Кососимметрическую матрицу a размера (3 3) можно сформировать в
~ ×
среде МatLab на основе функции cross. А именно, для строки a =[ax ,ay ,az ] создадим M-функцию s33(a), возвращающую (3×3) матрицу вида
function s = s33(a) |
(1.7) |
s=[ cross(a,[1 0 0]) ; cross(a,[0 1 0]) ; |
|
cross(a,[0 0 1])]; |
|
~2 |
|
М-функция ss33(a) формирует симметрическую (3×3)-матрицу a |
|
function s = ss33(a) |
(1.8) |
s = s33(a) s33(a); |
|
М-функция sp33(a,b) функция двух векторов-строк a , b (угловые скорость и ускорение вращающегося тела: a =ω, b =ε )
function s = sp33(a, b) |
(1.9) |
s = ss33(a) + s33(b); |
|
M-функции для пространственной системы сил.
В кинетостатике и статике при вычислении бивектора W силы
F=[X Y Z] можно применять (3×6)-матрицу радиус-вектора r =[x, y, z], опре-
деляемую М-функцией s36(r) , W=F*s36(r) |
|
|
function s = |
s36(r) |
(1.10) |
s = [s33(r), |
eye(3)]; |
|
где eye(3) = diag[1 1 1] — единичная диагональная матрица третьего порядка.
M-функция координат r = [x y z] точки приложения для силы F = nF e , заданной её модулем и ортом e , где nF=| F | =norm(F), e = [ex ey ez]
function s |
= s16(r, e) |
(1.11) |
s=e*s36(r) |
|
|
M-функция координат для силы F = [0 Y Z], вводимой без нуля как, F = [Y Z]
function s = sx26(r) |
(1.12) |
s=s36(r)
s=s(1, :) = [ ];
8

M-функция координат r = [x y z] силы F = [X 0 Z], вводим как F = [X Z]
function s = sy26(r) |
(1.13) |
s=s36(r)
s=s(2, :) = [ ];
M-функция координат силы F = [X Y 0], вводим как F = [X Y]]
function s = sz26(r) |
(1.14) |
s=s36(r)
s=s(3, :) = [ ];
M-функция координат для силы F = [X 0 0] ,вводим как F = X , без нулей
function s = syz16(r) |
(1.15) |
s=s36(r) s=s(1, :) ;
M-функция координат для силы F = [0 Y 0], вводим как F = Y
function s = sxz16(r) |
(1.16) |
s=s36(r) s=s(2, :) ;
M-функция координат для силы F = [0 0 Z], вводим как F = Z
function s = sxy16(r) |
(1.17) |
s=s36(r) s=s(3, :);
M-функция для момента пары сил с моментом M = M e = nMe , e =[ex ey ez ]– вектор-строка орта e . M = nM*e, бивектор пары W = nM*sp16(e)
function s |
= sp16(e) |
(1.18) |
|||
s |
= [ex |
ey |
ez |
0 0 0] |
|
% |
или s |
= [e, |
zeros(1, 3)]; |
|
M-функции для плоской системы сил, расположенной в плоскости Oxy
Вычисление (2×3)-матрицы координат для случая плоской системы сил, когда r=[x y], F = [X Y], W = F*s23(r) - бивектор силы F
function s |
= |
s23(r) |
(1.19) |
||
s = [0 -1; |
1 |
0]*[r;[0 –1]; [1 0]]’; |
|
||
Бивектор скалярного момента пары сил M Z = ± | M |, W = nM [±1,0,0] |
|||||
function s |
= |
sp13(M) |
(1.20) |
||
s = M * [1 |
0 |
0]; |
|
Бивектор единичной силы e , заданной строкой e = [ex ey], при r = [x y] function s = s13(e)
s = e * s23(r);
9
1.2.Конкатенация матриц и линейных функций
Вертикальная конкатенация (сцепление) двух матричных выражений, содержащих один и тот же столбцевой множитель {X}, вида A1{X} и A2{X}, где A1, A2 — прямоугольные матрицы с одинаковым числом столбцов, приводит к сцепленному по вертикали выражению A{X}, где матрица A = [A1; A2]
— результат вертикальной конкатенации матриц A1 и A2.
Горизонтальная конкатенация двух выражений XB1 и XB2, где X есть n-строка, B1 и B2 — прямоугольные матрицы с n-строками, n1 и n2-столбцами, приводит к объединенному выражению XB, где B = [B1, B2] — nx(n1+n2)- матрица, результат горизонтальной конкатенации матриц.
Любую сумму произведений вида ax +by +cz + ku можно представить как произведение вектора-строки (матрицы-строки) [a,b,c,k] образованной из первых сомножителей и вектора-столбца [x, y, z,u] ' вторых сомножителей
или наоборот, — строки вторых сомножителей, умноженной на столбец первых сомножителей. Например, возьмем три суммы произведений и каждую из них представим в виде произведения строки на столбец:
a1x +b1 y +c1z + k1u =[x, y, z,u][a1,b1,c1,k1 ]' |
(1.21) |
a2 x +b2 y +c2 z + k2u =[x, y, z,u][a2 ,b2 ,c2 ,k2 ]'
a3 x +b3 y + c3 z + k3u =[x, y, z,u][a3 ,b3 ,c3 ,k3 ]'
Три линейных алгебраических выражения объединяем по горизонтали в одно матричное строчное выражение VA, где V=[ x y z u], A=[a1 a2 a3 ; b1 b2 b3 ; c1 c2 c3 ; k1 k2 k3]. Либо сцепляем вертикально в столбцевое выражение A΄{V}
при {V}=V΄.
Сумма произведений строчных матриц и согласованных с ними по умножению прямоугольных матриц объединяется в произведение удлиненной вектор строки и блочной матрицы, полученной в результате вертикальной конкатенации прямоугольных матриц. Например 1x3 – вектор-строка вида:
[x, y, z]A33 +[u, v]B23 + wC13 =[x, y, z, u, v, w][A33; B23; C13 ], |
(1.22) |
где вторые матричные сомножители имеют размеры (3x3), (2x3), (1x3).
1.3.Перепроецирование векторов и тензоров второго ранга
Пусть имеется две декартовы системы координат: Oxnynzn и Oxsyszs, первую назовем — n-системой или исходной системой, вторую — повернутой системой или — s-системой, с отсчетом углов от исходной системы. Расмотрим в них один и тот же вектор r , он определяется строками проекций, которые обозначим соответственно rn=[xn yn zn] и rs=[xs ys zs].
Верна матричная формула перепроецирования вектора r из n-системы в s-систему, связывающая его проекции на оси s-системы, n-системы:
10
rs = rnCsn |
(1.23) |
Здесь Csn — матрица поворота s-системы от n-системы, элементами которой являются направляющие косинусы осей повернутой системы, т.е. скалярные произведения соответствующих ортов осей,
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
is in |
|
|
js in |
|
ks in |
in |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(1.24) |
|||||
|
|
|
|
|
|
|
|
js jn |
|
|
|
|
|
|
|
|
[is |
js ks ] |
||||||||||||||||||||||||
Csn = is jn |
|
ks jn |
≡ jn |
|
||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
js kn |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
is kn |
ks kn |
kn |
|
|
|
|
|
|
|
|
Здесь последнее выражение представляет собой скалярное произведение столбца ортов исходной системы и строчной матрицы ортов повернутой систем, которое раскрывается формально по правилу перемножения матриц.
Мнемоническое правило: в формуле (1.23) индексы s, n повторяются в том же порядке, а в произведении матриц (1.24) – в обратном порядке.
Общая формула перепроектирования любого вектора a |
из n-системы на |
||||||||||||
s-систему и обратно |
(1.25) |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
′ |
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|||||||||
as = anCsn , an = asCns при Сns =[is js ks ]' [in jn kn ] = Csn |
|
Примем, например, за исходную систему Ox1y1z1 и за повернутую — “нулевую” систему Oxyz. Тогда по формуле (1.23) находим строку r=[x y z]
проекций вектора r |
на “нулевую” систему: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
′ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
r = rC |
при |
С |
=[i |
|
j |
|
|
|
|
|
|
|
|
|
j k ] |
|
|||||||||||||
|
k ] [i |
|
|||||||||||||||||||||||||||
1 |
01 |
|
01 |
1 |
|
1 |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
И наоборот, перепроецирование r |
из системы Oxyz в Ox1y1z1 имеет вид: |
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
′ |
|
|
|
|
|
|
|
|
|
|
|
|
|
(1.26) |
r1 = rC10 |
при |
С10 =[i |
|
j |
|
|
|
|
[i1 |
|
j1 k1] |
||||||||||||||||||
|
|
k ] |
|
|
Пусть например, имеется четыре системы координат, из них вторая Ox2y2z2 повернута от первой Ox1y1z1 на угол ψ вокруг Oz1 , третья от второй повернута на угол θ вокруг Ox2 , четвертая от третьей на угол ϕ вокруг Oz3.
Векторы-строки вектора r на эти оси обозначим rS = [xS yS zS ] , s =1,2,3,4 . По формулам (1.24), (1.23) находим
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
cψ −sψ |
0 |
|
|
|
|||
|
|
|
i1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
C |
|
(ψ, z) = |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= sψ |
cψ |
0 |
, |
r |
= rC |
|
21 |
|
j |
|
|
i |
|
j |
|
k |
2 |
|||||||||||||||
|
|
|
1 |
|
2 |
2 |
|
|
|
|
|
|
|
2 |
1 21 |
||||||||||
|
|
|
|
k |
|
|
|
|
|
|
|
|
|
|
|
|
0 |
0 |
1 |
|
|
|
|||
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
i |
|
|
|
|||||||
|
|
|
2 |
|
|
|
|
|
|
|||
|
|
|
|
|
||||||||
C32 (θ, x) = j2 |
|
i3 |
||||||||||
|
|
|
|
|
|
|
|
|
||||
|
|
k2 |
|
|
|
|
||||||
|
|
cϕ −sϕ |
||||||||||
C |
43 |
(ϕ, z) = sϕ cϕ |
||||||||||
|
|
|
0 |
|||||||||
|
|
0 |
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
0 |
0 |
|
|
|
|
(1.27) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
j3 |
|
|
= |
|
cθ −sθ |
, r3 = r2C32 |
|||||||||
k3 |
0 |
|
|
||||||||||||
|
|
|
|
|
|
0 |
|
sθ |
cθ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
, r |
= r C |
43 |
= rC C |
32 |
C |
43 |
,cψ ≡ cosψ, sψ ≡ sinψ |
||||||
|
|
4 |
|
3 |
1 |
21 |
|
|
|
||||||
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Отсюда находим матрицу поворота C41 четвертой системы от первой системы и матрицу C14 обратного поворота.
11