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

За материальную точку принимают тело пренебрежительно малых размеров,
атакже бесконечно малые частицы тела, на которые условно разделяют тело.
Вследующем разделе покажем, что каждому телу можно сопоставить материальную точку, посредством условной концентрации масс тела в его центре масс и параллельного переноса всех приложенных сил в центр масс.
2.4. Пример
Подпружиненный груз массы m движется без трения по горизонтальной плоскости вдоль оси Ox. Жесткость горизонтальной пружины равна c. В начальное мгновение t0 = 0 груз располагался на отрицательной ветви оси Ox на расстоянии b и имел модуль скорости v0, причем начальная скорость была направлена против направления оси (рис.1).
Найти:
•динамическое уравнение
•сформулировать задачу Коши и посредством её решения найти кинематическое уравнение движения точки
•Найти реакцию опоры N.
Решения найти в символьном виде, а также - при численных значениях b
= 0,2, v0 = 4, c = 8, m = 2, x&0 = −4. .
Решение.
1. Выпишем все данные задачи в буквенной и численной форме. По этим данным оформим рисунок, на который нанесем условия задачи. На ри-
сунке покажем начальное состояние объекта и его произвольное «положи-
тельное» состояние с произвольными положительными значениями координаты и скорости [v > 0, x > 0]. На объекте в положительном состоянии покажем приложенные силы (рис. 1), которые условно перенесем в центр масс (геометрический центр) тележки.
|
y |
N |
|
|
|
v0 |
|
F |
|
G |
|
|
j |
|
v |
|
|
|
|
|
x0 |
O |
i |
x |
x |
Заданное началь- |
Текущее «положительное» |
ное состояние |
состояние объекта |
[v0x <0, x0 <0]. |
[x > 0, x > 0]. |
|
& |
Рис. 1. Колебания тела вдоль горизонтальной оси |
|
m = 2 кг, с = 8 Н/м, x0= – b = |
– 0,2 м, x0 = −v0 = −4 м/с, F = −cxi , |
&
R = F +G + N , G = −mg = −mg j , N = N j
17

2. Рассматривая на рисунке «положительное» состояние объекта, по формулам (2.1), (2.2) составим векторное динамическое уравнение:
ma = F +G + N .
Проецируем его на оси x, y: |
my = −bmg + N, при |
y ≡ 0. |
mx = −cx, |
||
&& |
&& |
&& |
Движение вдоль оси Oy отсутствует, поэтому второе уравнение системы фактически является алгебраическим уравнением. Из него находим реакцию связи (опоры) N = mg. Остается одно дифференциальное (динамическое) уравнение m&x&+ cx = 0 .
3. Упрощаем вид динамического уравнения с целью уменьшения количества буквенных констант и присоединяем начальные условия, показанные на рисунке, начальный фазовый вектор. Получаем задачу Коши, записанную в символьном виде с тремя символьными константами: k,b,v0
&& |
+k |
2 |
x = 0 |
при k = с/ m ; x(0) = x0 = – b, x(0) = vx0 = – v0. |
x |
|
4. Решаем задачу Коши в символьном виде. (В сложных задачах ищем численное решение задачи Коши на компьютере.) В данном случае существует символьное решение вида:
x = x0 cos (kt) + (vx0/ k) sin(kt) = –b cos(kt) – (v0/ k) sin(kt).
Оно соответствует гармоническим колебаниям с амплитудой A, со смещен-
ной фазой (kt–α), с частотой k : |
x = A cos(kt–α) |
|||||
при A = |
2 |
+(vxo / k) |
2 |
|
& |
/ A) = −arccos(−b / A), |
xo |
|
, α = sgn(x0 )arccos(x0 |
||||
где sgn(x0 ) = −1, поскольку |
x0 = −v0 < 0. |
|
||||
|
& |
|
|
|
& |
|
5. Подставляя численные значения буквенных параметров, получаем окончательное кинематическое уравнение движения груза:
k = 2, A = 2.01,α = −1.67, x = −0,2cos(2t) −2sin(2t) = 2,01cos(2t +1,67).
Замечание. Динамическое уравнение, полученное для «положительного» состояния, как правило верно и для прочих состояний. Но в некоторых задачах целесообразна проверка правильности полученного дифференциального уравнения для различных «неположительных» динамических состояний. А именно, в задачах с сухим или квадратичным относительно скорости трением проверка приведет к появлению в динамическом уравнении двойного знака (±), поскольку в выражении силы имеется двойной знак. В задачах с вязким трением, выражаемым линейной нечетной функцией, такая проверка не требуется.
2.5. Второй закон динамики в кинетостатической форме Даламбера
Векторное уравнение (2.1) в результате перенесения всех векторов в левую часть и введения нового обозначения записывается в следующей даламберовой форме:
n |
|
|
|
|
|
|
|
|
(2.8) |
|
|||||||||
∑Fs + S = 0 при обозначении S = −ma . |
|
||||||||
s=1 |
|
||||||||
|
|
|
|
|
|
|
|
|
18 |

Это уравнение в проекциях на оси декартовой системы координат представляет собой систему трех скалярных уравнений
n |
n |
n |
(2.9) |
∑X s + Sx = 0, ∑Ys + S y = 0 , ∑Zs + Sz = 0. |
|
||
s=1 |
s=1 |
s=1 |
|
Вектор S = −ma имеет размерность силы, его называют силой инерции материальной точки. Уравнения (2.8) трактуются следующим образом: Любая материальная точка в ИСО движется так, что в каждый момент времени система приложенных к ней сил уравновешивается даламберовой силой инерции. Иными словами — система сил, включая силу инерции, удовлетворяет векторному уравнению “статического” равновесия (2.8) или — системе трех скалярных уравнений (2.9), эквивалентных векторному уравнению.
Силу инерции обычно называют фиктивной силой, поскольку она не приложена к материальной точке, все реальные приложенные силы уже содержатся в векторной сумме, стоящей в уравнении (2.8). Но можно условно применять и другую трактовку даламберовой силы инерции.
Будем применять упрощенную, формальную трактовку сил инерции. В случае, когда второй закон динамики записан в основной форме (2.1), считаем, что сила инерции не существует. Если закон записан в кинетостатической форме (2.8), то будем формально считать, что сила инерции существует и приложена к материальной точке. В том и другом случае для материальной точки записывается одно и то же векторное уравнение, но в разных обозначениях. Кинетостатическая форма удобна при определении реакций связей и решении прочностных задач в дисциплине «сопротивление материалов».
Условно представим себе, что материальная точка имеет вид бесконечно малой связи-капсулы, заполненной материалом. Действующие силы приложенны к внешней поверхности связи-капсулы, материал оказывает давление D изнутри. По закону Ньютона векторная сумма сил, действующих на безмассовую связь-капсулу равна нулю, поскольку mкапсa = 0 , отсюда
n
∑Fs + D = 0 . Сравнивая эти уравнения с (2.8), заключаем что даламберову
s=1
силу инерции можно трактовать как силу давления материальной точки на связи. Отметим, что в такой трактовке вторая массовая сила - сила тяжести материала также приложена изнутри. Согласно с замечанием выдающегося французского ученого П.Аппеля (1855-1930) при бросании камня поверхность ладони (связь) ощущает силу инерции камня в виде давления на руку. Но эта реальная сила приложена к связи, не приложена к камню и в этом смысле она считается фиктивной. Уравнение Даламбера (2.8) означает, что векторная сумма сил, приложенных к материальной точке и инерционных противодействующих сил, создаваемых этой точкой, равна нулю, удовлетворяют закону равенства действия и противодействия.
Геометрические (голономные) безмассовые связи объектов имеют разнообразные представления. Например, во вращающимся маховике, состоящем из множества материальных точек, можно мысленно выделить множест-
19

во связей, образующих безмассовую пористую структуру, к которой приложено множество внутренних сил давления: массовых сил инерции, массовых сил тяжести и приложенных внешних сил. Объединенная система сил вызывает деформации связей маховика и даже его разрушения при большой угловой скорости. Названная система сил удовлетворяет уравнениям статического равновесия.
2.6. Закон динамики относительного движения
Рассмотрим любую неинерциальную (подвижную) систему отсчета Охуz, которая произвольно движется в некоторой инерциальной (неподвижной) системе отсчета O1x1y1z1 Считаем, что подвижная система в общем слу-
чае имеет некоторую скорость vO и ускорение aO полюса О. При этом она вращается вокруг полюса с некоторой угловой скоростью ω и угловым ускорением ε, где векторы ω и ε могут изменяться по величине и направлению, образуя между собой переменный угол.
Радиус-вектор материальной точки в подвижной системе отсчета и его вектор-строку обозначим r(x, y, z) и r =[x y z]. Ускорение a1 точки в непод-
вижной системе О1 х1 у1 z1 назовем абсолютным ускорением, ускорение a1 точки в подвижной системе Охуz – относительным ускорением. Разложим a1
по теореме Кориолиса на три составляющие: относительное, кориолисово и переносное ускорения
a1 |
= |
a |
+ |
ac |
+ |
ae |
= &&+ ω |
× & +ω |
× ω |
× |
r ) |
+ε |
× |
r |
+ |
aO . |
(2.10) |
|
|
|
|
r 2 r |
( |
|
|
|
|
|
Здесь переносное ускорение представлено суммой трех составляющих осестремительного ускорения, вращательного ускорения и ускорения полюса, производные от r вычисляет “относительный” наблюдатель, считающий
“свою” систему Охуz неподвижной, её орты i , j ,k – постоянными по направ-
лению.
Подставляя выражение (2.10) в уравнение (2.1), получим векторное динамическое уравнение относительного движения материальной точки в об-
щем виде |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
&& |
|
|
& |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
+ω×(ω×r ) +ε ×r + aO − R / m = 0 , |
|
||||||||||||||||||||||
r |
+ 2ω×r |
(2.11) |
||||||||||||||||||||||
|
|
|
|
|
|
|
n |
|
|
n |
||||||||||||||
|
где |
|
= ∑ |
|
|
|
|
|
|
|
|
|
||||||||||||
|
R |
Fs = ∑( X si |
+Ys j + Zs k ) . |
|
||||||||||||||||||||
|
|
|
|
|
|
s=1 |
|
|
s=1 |
|
Таким образом, закон (2.11) движения материальной точки в подвижной системе отсчета отличается от закона Ньютона в неподвижной системе &&r1 − R / m = 0 наличием в общем случае четырех дополнительных слагаемых.
Векторное уравнение (2.11), эквивалентно следующему матричному дифференциальному уравнению относительного движения материальной точки:
&& |
& % |
%2 |
% |
(2.12) |
r |
+ 2rω + r(ω |
+ε) + aO − R / m =[0 0 0], |
|
20

где в круглых скобках поставлена квадратная матрица, образованная из про-
екций векторов ω и ε по правилам (1.4)-(1.6). Уравнение (2.12) раскрывается в виде системы трех дифференциальных уравнений, записанных в строку.
Напомним, что ω% и ε% - кососимметрические (3x3) матрицы, а ω%2 симметрическая матрица.
В обозначениях системы MATLAB уравнения (2.12) записываются следующим образом:
D2r+2*Dr*w+r*[w^2+e]+a0-R/m=zeros(1,3) (2.13)
Здесь введены латинские обозначения ω → w, ε → e; вектор-строка равно-
n
действующей системы приложенных внешних сил R = ∑Fs =[Rx , Ry , Rz ] ;
s=1
вектор-строка ускорения материальной точки относительно подвижной системы отсчета D2r = [D2x, D2y, D2z]; строка относительной скорости Dr = [Dx Dy Dz] ; строка абсолютного ускорения начала координат подвижной
системы Oxyz a0 = [aox,aoy,aoz].
Пусть силы, приложенные к материальной точке, разделены на две группы: силы, заданные в проекциях на подвижную систему Oxyz с равнодействующей вектор-строкой R0 = [Rx0 , Ry0 , Rz0 ] и силы, заданные в проекциях на
ИСО Ox1y1z1 с равнодействующей вектор-строкой R =[R1 |
, R1 |
, R1 |
]. В таком |
||||||||||||||||||||||||||||||||||||||||||||||||||
1 |
|
|
|
x1 |
|
y1 |
|
|
z1 |
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||
случае следует перепроецировать R1 и в формуле (2.13) полагать |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
R = R0 + R1 C01 |
|
|
|
|
|
|
|
|
|
|
|
|
(2.14) |
||||||||||||||||||||||||||||
|
|
|
|
|
|
′ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
k ; k1 i , k1 |
j ,k1 k ] |
|||||||||||||||||||||||||||||||||||||
C01 =[i1 , j1, k1] [i , j ,k ] =[i1 i , i1 j , i1 k ; j1 i , j1 j , j1 |
|||||||||||||||||||||||||||||||||||||||||||||||||||||
Левую часть матричного динамического уравнения (2.13) запишем сле- |
|||||||||||||||||||||||||||||||||||||||||||||||||||||
дующим образом: |
|
|
|
|
|
|
|
|
|
|
|
|
(2.15) |
|
|||||||||||||||||||||||||||||||||||||||
V = D2r+2*Dr*s33(w) + r*sp33(w, e)+a0 – R/m , |
|
|
|
|
где применены M-функции s33 и sp33 вида (1.7), (1.9), а все векторыстроки считаем представленными в проекциях на Oxyz. Система MATLAB раскроет матричное выражение (2.15), возвратив его в виде вектор-строки W из трех выражений. Сохранение только первого выражения достигается функцией W =W(:,1), сохранение первого и второго выражения достигается функцией W(:,3)=[ ]. Все символьные обозначения должны быть предварительно представлены функцией syms...real. Затем выражение (2.15) (вместе с начальными условиями r(0) =r0, Dr(0) = v0 вносится в операцию Dsolve() , возвращающую символьное или численное решение.
Уравнение (2.12) можно представить системой в форме Коши для век- тор-строки V = [v,r] , где v = r& . Имеем с добавлением нулевых слагаемых
21