Нелинейная динамика. 2012. Т. 8. № 2. С. 219-229. Полнотекстовая версия в свободном доступе http://nd.ics.org.ru
УДК: 531.36,531.38,531.66
М8С 2010: 70Е55,70Е35,70К42, 70Е50
К динамике твердого тела, несущего материальную точку
А. П. Маркеев
Рассматривается система, образованная «несущим» твердым телом и «несомой» материальной точкой, которая движется по заданному закону вдоль кривой, жестко скрепленной с телом. Движение происходит в однородном поле тяжести над неподвижной абсолютно гладкой горизонтальной плоскостью. При движении несущее тело может соударяться с плоскостью. Величина коэффициента восстановления при ударе считается произвольной. Получены уравнения, описывающие как свободное движение системы над плоскостью, так и моменты соударений. Указано несколько частных решений уравнений движения, и в некоторых случаях исследована их устойчивость.
В случае, когда несущее тело динамически симметрично, а материальная точка движется вдоль оси симметрии по произвольному закону, найдено общее решение уравнений свободного движения тела в квадратурах, обобщающее решение, соответствующее классической регулярной прецессии в случае Эйлера.
Показано, что для существования поступательного движения несущего тела в режиме его свободного полета над плоскостью в общем случае необходимо, чтобы материальная точка двигалась относительно тела в соответствии с законом площадей.
Ключевые слова: динамика систем твердых тел, соударение, периодическое движение, устойчивость
Получено 19 марта 2012 года После доработки 25 марта 2012 года
Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (11-01-00322), Программы поддержки ведущих научных школ (НШ - 4149.2012.1) и гранта Правительства РФ для государственной поддержки научных исследований, проводимых под руководством ведущих ученых в российских образовательных учреждениях ВПО (дог. №11.034.31.0039).
Маркеев Анатолий Павлович [email protected]
Институт проблем механики им. А. Ю. Ишлинского РАН 119526, Россия, г. Москва, пр. Вернадского, д. 101, корп. 11
1. Уравнения движения
Рассмотрим движение твердого тела массы М в однородном поле тяжести над неподвижной горизонтальной плоскостью. Движение отнесем к неподвижной системе координат ОаХаУаЕа с началом в точке Оа этой плоскости, ось Оа2а направим вертикально вверх (рис. 1). Пусть Охуг — жестко связанная с телом система координат, ее начало поместим в центр масс О тела, а оси Ох, Оу, Ог направим вдоль главных центральных осей инерции тела. Соответствующие главные моменты инерции обозначим через ’]х,3у,^г. Ориентацию тела зададим при помощи матрицы А косинусов углов между осями неподвижной и связанной систем координат. В частности, единичный вектор п вертикальной оси OaZa в связанной системе координат Охуг задается компонентами азі, аз2, аээ-
Рис. 1. Системы координат для описания движения тела, несущего материальную точку.
Предположим, что вдоль некоторой кривой, жестко скрепленной с телом, движется материальная точка Р массы т. Ее положение в системе координат Охуг задается вектором ОР, компоненты х,у,г которого — заданные функции времени: х = х(Ь), у = у(Ь), г = г(Ь). Положение центра масс С рассматриваемой системы — твердое тело, несущее материальную точку, задается вектором, который в системе Охуг определяется равенством ОС = т/(т +
Через 5 будем обозначать точку поверхности «несущего» твердого тела, ближайшую к неподвижной горизонтальной плоскости ОаХаУа. Ее положение в системе Охуг задается вектором ОБ, компоненты которого х*, у*, г* определяются формой поверхности тела и его ориентацией в системе OaXaYaZa.
Во время движения «несущее» тело может соударяться с плоскостью Оа,XaYa. Считаем, что плоскость абсолютно гладкая, а коэффициент восстановления е при ударе произволен
Уравнения движения получим из теоремы об изменении количества движения системы
0 и теоремы об изменении ее кинетического момента К относительно центра масс С.
Уравнения свободного движения тела над плоскостью. В случае свободного движения имеем два векторных уравнения
п
+ М) ОР.
(1.1)
где д — ускорение свободного падения. Если Хс, Yc, Zc — компоненты радиус-вектора ОаС центра масс С в неподвижной системе координат ОаХаYaZa, то первое из векторных уравнений движения (1.1) запишется в виде трех скалярных уравнений
(1.2)
Получим скалярную форму второго из уравнений (1.1). Для этого перепишем его в виде
| 1^хК = 0, аі
(1.3)
где первое слагаемое представляет собой вектор, компоненты которого равны производным по времени от проекций вектора К на оси жестко связанной с несущим телом системы координат Охуг. Вектор ш в (1.3) — мгновенная угловая скорость несущего тела, ее проекции на оси Ох, Оу и Ог будем обозначать, соответственно, через шх, шу и .
Воспользовавшись хорошо известными правилами [1], можно получить следующее выражение для вектора кинетического момента рассматриваемой системы относительно ее центра масс:
тМ
К — Зш + ц ОР х Уг
ц
т + М'
(1.4)
Здесь через Л обозначена симметрическая матрица вида
Зх + ц(у2 + г2)
цху цхг
— і
цху Зу + ц(г2 + х2)
цуг
-цхг -цуг Зх + ц(х2 + у2)
— I
(1.5)
\г — вектор скорость точки Р относительно несущего тела, имеющий в системе координат Охуг, жестко связанной с телом, компоненты х, у, г.
Замечание. Можно сказать, что Л — матрица тензора инерции для точки О некоторого воображаемого твердого тела, представляющего собой рассматриваемое несущее твердое тело, к которому в точке с координатами х, у, г жестко прикреплена материальная точка массы ц. Очевидно, что элементы этой матрицы могут принимать только конечные значения. Будем предполагать, что <1еЛ отличен от нуля, а элементы матрицы Л-1, обратной к матрице Л, — ограниченные функции времени. Последнее предположение исключает, например, такие экзотические случаи, когда при каких-либо значениях £ геометрия масс упомянутого воображаемого тела отвечает геометрии масс бесконечно тонкого стержня.
Компоненты вектора К в системе координат Охуг вычисляются по формулам Кх — [Зх + ц(у2 + г2)]шх - цхушу - цхгшх + ц(уг - гу) (х, у, г).
(1.6)
Здесь и в дальнейшем невыписанные соотношения получаются из выписанных круговой перестановкой символов, указанных в круглых скобках.
Векторное уравнение (1.3) запишется в виде следующих трех скалярных уравнений:
[Зх + ц(у2 + г2)]их - цхуйу - цхгшх + (Зх - Зу)шу+ +ц[(х^х + ушу + гшх)(ушх - гшу) + 2(уу + гг)шх-
(1.7)
-2х(ушу + гшх) +уг - гу] — 0
(х,у,г).
З
Согласно (1.2), в случае свободного движения центр масс системы C движется по параболе или вертикальной прямой. Движение же несущего тела относительно его центра масс O описывается уравнениями (1.7). К этим уравнениям следует еще добавить кинематические уравнения Пуассона:
ац = uzai2 - Uyai3 (1,2, 3; x,y,z) i = 1,2,3. (1.8)
Если решение уравнений (1.2) и (1.7), (1.8) найдено, то положение центра масс O несущего тела находится из соотношения
ОйО = ОйС — А —2Ч7 OP.
m + M
Система динамических уравнений (1.7) может быть исследована независимо от кинематических уравнений (1.8). Из (1.1) следует, что эта система допускает интеграл
K% + Ky + KZ = K2 = const, (1.9)
где величины Kx,Ky,Kz задаются равенствами (1.6).
Уравнения, соответствующие моментам с о уд а р е н и й. При помощи классической теории удара [2, 3] получим связь кинематических величин, характеризующих движение системы до и после соударения несущего тела с плоскостью. Считаем, что положение центра масс системы в абсолютном пространстве, положение точки P в несущем теле и ориентация несущего твердого тела в абсолютной системе координат OaXaYaZa во время удара не изменяются; следовательно, во время удара неизменны координаты x*,y*, z* точки S поверхности несущего тела, которой оно соударяется с плоскостью. Кроме того, предполагаем, что во время удара не изменяется вектор Vr относительной скорости точки P, что обеспечивается соответствующим управлением, задающим режим x = x(t),y = = y(t),z = z(t) движения точки P относительно несущего тела.
Так как плоскость абсолютно гладкая, то ударный импульс I вертикален и классические уравнения теории удара имеют вид
Q+ - Q- = In, K+ - K- = Id. (1.10)
Здесь верхними значками « —» и «+» обозначены значения величин до и после удара. Через d в (1.10) обозначен вектор CS х n; его компоненты dx ,dy, dz не изменяются во время удара:
mm
Лх = \У* ~ т + МУ) 033 “ V"* “ т + М*) °'32 ^Х’У’ г; а'31’ а'32’ а'33^
Принимая во внимание, что О — (т + М)\с, где \а — вектор абсолютной скорости центра масс, и учитывая выражение (1.4) для вектора кинетического момента, из (1.10) находим
Х+=Л~ Ї'+ = Ї'СГ, 2ї = 2с+^м: (1-11)
ш + — ш_ +1 з-1 а. (1.12)
Для вычисления величины модуля ударного импульса I воспользуемся соотношением
— -еу-, (1.13)
где e — коэффициент восстановления, а vn = (Vs, n) — проекция на вертикальную ось OaZa вектора абсолютной скорости точки S поверхности несущего тела. Можно показать, что
Vn = ZC-----(v.r, n) + (ш, d). (1.14)
m + M
Из (1.11)—(1.14) находим
1 = [Zc - т (Уг,п) + (u;“,d)]. (1.15)
Здесь введено обозначение
к = —Ц-р +(J-4d). (1.16)
т + М v v '
Величина Vr в (1.15) вычисляется в момент соударения.
Соотношения (1.2), (1.7), (1.8), (1.11), (1.12) с учетом равенств (1.15), (1.16) образуют замкнутую систему уравнений, описывающую движение рассматриваемой материальной системы «несущее тело-материальная точка» на промежутках времени, включающих в себя как интервалы свободного полета, так и моменты соударений несущего тела с плоскостью.
2. О поступательном движении несущего тела
Для приложений (например, в динамике роботов) представляет интерес вопрос о возможности задания такого движения точки P относительно несущего тела, чтобы оно могло совершать поступательное движение.
Условия существования поступательного движения. Положив в уравнениях (1.7) их = 0, uy = 0, uz = 0, получим три условия
yz — zy = 0, zx — xz = 0, xy — yx = 0, (2.1)
необходимые для существования поступательного движения несущего тела на этапе его свободного движения над плоскостью.
Эти условия, очевидно, выполнены в том частном случае, когда точка P движется относительно несущего тела по произвольной прямой с постоянной, произвольной по величине, скоростью:
x = xo + x о t, y = yo + y ot, z = zo + zot. (2.2)
В общем же случае условия (2.1) выполнены только тогда, когда движение точки P относительно несущего тела происходит в соответствии с законом площадей [4], то есть когда выполняются равенства
yz - zy = Сх, zx - xz = Cy, xy - yx = Cz, (2.3)
где cx,cy,cz — произвольные постоянные. Следовательно, траектория точки P относительно несущего тела должна быть кривой, лежащей в фиксированной относительно тела плоскости, которая перпендикулярна вектору с с компонентами cx,cy,cz; при этом радиус-вектор OP за равные промежутки времени должен заметать равные площади. В частном случае, когда с = 0, движение точки P должно происходить вдоль прямой, проходящей через центр масс O несущего тела, но закон ее движения вдоль этой прямой произвольный (в отличие
от случая (2.2), когда прямая не обязательно проходит через центр масс тела, но движение вдоль прямой равномерное).
Поступательное движение несущего тела может существовать и при наличии его соударений с плоскостью. Действительно, пусть выполнены условия (2.1) и несущее тело движется над плоскостью поступательно. Предположим, что при этом материальная точка Р движется в несущем теле вдоль вертикальной прямой £*, которая проходит через центр масс О тела и пересекает поверхность тела в ее наинизшей точке 5, причем в этой точке касательная плоскость к поверхности тела перпендикулярна £*. В момент соударения вектор СБ вертикален, поэтому d = СБ х п = 0 и, согласно равенству (1.12), вектор угловой скорости несущего тела остается неизменным. Таким образом, в рассмотренном примере поступательное движение несущего тела существует не только при свободном движении рассматриваемой системы над плоскостью, но и на промежутках времени, включающих и моменты соударений несущего тела и плоскости.
Об устойчивости поступательного движения. Задача об устойчивости поступательного движения несущего тела является довольно сложной. В частном случае, когда несущее тело динамически симметрично, а точка Р движется вдоль оси симметрии по гармоническому закону, эта задача изучалась в статье [5].
Здесь рассмотрим устойчивость поступательного движения несущего тела по отношению к возмущениям величин Шх, Шу, шг на этапе свободного движения тела над плоскостью, причем ограничимся случаем, когда точка Р движется вдоль прямой, проходящей через центр масс несущего тела. Закон движения точки Р вдоль этой прямой не конкретизируется.
Для рассматриваемого относительного движения точки Р постоянные сх, су, сх в соотношениях (2.3) равны нулю и выражения (1.6) для проекций вектора кинетического момента на оси системы координат Охуг будут линейными однородными функциями величин
Шx, Шу, :
Кх = и + Му2 + г2)]шх - цхуШу - цхгш? (х,у,г). (2.4)
Покажем, что по отношению к возмущениям величин шх,шу,шг поступательное движение несущего тела на этапе его свободного движения над плоскостью является устойчивым. Это можно сделать при помощи теоремы Ляпунова об устойчивости движения [6].
Действительно, по предположению матрица (1.5) имеет ограниченную обратную, поэтому из (2.4) следует, что задачи об устойчивости по отношению к возмущениям величин шх,шу,шг и Кх,Ку,Кг эквивалентны.
В качестве функции Ляпунова V примем функцию из левой части интеграла (1.9):
V = К2 + Ку2 + К2. (2.5)
Из сказанного выше и очевидной знакоопределенности функции (2.5) относительно переменных Кх, Ку, Кх сразу следует устойчивость поступательного движения тела относительно возмущений величин Шх,Шу,Шг.
3. О свободном движении динамически симметричного несущего тела
Предположим, что Зх = ,1у, то есть считаем, что ось Ог является осью динамической симметрии несущего тела. И пусть точка Р движется вдоль оси симметрии по заданному
закону: x = 0,y = 0,z — z(t). В этом случае компоненты (1.6) вектора кинетического момента K в системе координат Oxyz будут иметь вид
Kx = (Jx + №z )wx, Ky = (Jx + №z )wy, Kz = Jz Wz, (3.1)
а система динамических уравнений (1.7), описывающая свободный полет тела над плоскостью, имеет вид
^z2)lox] + {Jz — Jx — liz2)u)yUJz = 0, (3-2)
+ fiz2)ujy] — (Jz — Jx — i~iz2)luzlux = 0, (3-3)
= 0- (3-4)
Эти уравнения имеют интеграл (1.9), который с учетом равенств (3.1) можно записать в виде
(Jx + ц-z2)2 (и2 + W2) + J2w2 = K2 = const. (3.5)
Кроме того, имеется интеграл
wz — = const. (3.6)
Стационарное вращение вокруг оси симметрии и его устойчивость. Пусть и — произвольная ось фиксированного направления в неподвижной системе координат OaXaYaZa. Направляющие косинусы этой оси во вращающейся вместе
с телом системе координат Oxyz обозначим через ax, ay, az. Очевидно геометрическое
тождество
o?x + Q*y + а2 — 1, (3.7)
а также соотношение
(Jx + /j.z2 )(wxax + wy ay) + Jz wz az — Ku — const. (3.8)
Последнее соотношение означает, что проекция кинетического момента K системы на лю-
бую ось, фиксированную в неподвижной системе координат, не изменяется со временем (см. второе из уравнений (1.1)).
Существует движение несущего тела, когда его ось симметрии имеет фиксированное направление в пространстве (и, например, совпадает с направлением оси и), а само тело вращается вокруг оси симметрии с постоянной угловой скоростью w°. Для такого стационарного вращения
wx ° wy ° wz wz, ax 0 ay 0, az 1.
При помощи теоремы Ляпунова об устойчивости покажем устойчивость этого вращения относительно возмущений величин wx, wy, wz, ax, ay, az. Будем полагать, что — 0. (Если - 0, то движение тела поступательно и, согласно предыдущему разделу, является
устойчивым по отношению к возмущениям величин wx, uy, uz.)
Введем возмущения Ci (i — 1, 2,..., 6) по формулам
wx — Съ wy — С2, wz — wz + Cз, ax — ay — С5, az — 1 + Сб.
Из (3.5)—(3.8) следует, что в возмущенном движении имеют место интегралы
V1 — K2 — J— (Jx + ^z2)2 (С2 + С2) = const, V2 — C3 — const,
V3 = ax + O? + a2 — 1 = 2С6 + С4 + c2 + c2 = const,
V4 = Ku — Jz Wz = (Jx + ^z2)(ClC4 + С2С5) + Jz (w° + Сз)Сб = const.
Функцию Ляпунова возьмем в виде такой суммы:
V = Vi + V,2 + V32 + V42. (3.9)
Она будет определенно положительной, если функция V4 будет знакоопределенной при тех значениях Ci (i = 1,2,..., 6), для которых Vi — 0, V2 — 0, V3 — 0. Из последних трех равенств получаем
6 = 0, 6 = 0, 6 = 0, 6 = -|(62 + 62) + ---- (3-10)
В выражении для 6 не выписаны члены выше второй степени относительно 6, 6.
При выполнении условий (3.10) функция V4 будет функцией двух переменных С4, С5:
V4 = ) + . . . ,
где многоточие обозначает совокупность членов выше второй степени относительно С4, С5. Эта функция, очевидно, является знакоопределенной. Тем самым устойчивость стационарного вращения доказана.
Общее решение уравнений движения. В рассматриваемом случае, когда Jx — Jy, x = 0, y = 0, общее решение уравнений движения несущего тела при произвольном законе движения z — z(t) материальной точки P вдоль его оси динамической симметрии можно получить в квадратурах. Это решение обобщает классическое решение задачи о движении динамически симметричного твердого тела вокруг неподвижной точки в случае Эйлера, которое отвечает регулярной прецессии тела.
Удобно рассматривать движение несущего тела в поступательно движущейся системе координат OXYZ, начало которой находится в центре масс O тела, а направление оси OZ совпадает с неизменным направлением вектора кинетического момента K в неподвижной системе координат OaXaYaZa. Ориентацию несущего тела относительно системы координат OXYZ зададим при помощи углов Эйлера ф,в,ф (рис. 2).
Из интеграла (3.6) и соотношения Jzwz — K cos в, где K — неизменный по величине модуль кинетического момента (см. рис. 2 и формулы (3.1)), следует, что угол нутации в постоянен во все время движения: в — в0 — const, причем
т , ,0
cos (Р = —~г^~ ■ (3.11)
K
При в = в0 имеем такие выражения для проекций вектора угловой скорости тела на оси Ox, Oy, Oz:
wx — sin в0 sin >£, wy — sin в0 cos ф, w° — t^cos в0 + ф. (3.12)
А проекции Kx,Ky кинетического момента на оси Ox,Oy, согласно (3.1), задаются равен-
ствами
(Jx + /j,z2)wx — K sin в0 sin ф, (Jx + ^z2)wy — K sin в0 cos ф. (3.13)
Рис. 2. К случаю свободного движения при наличии динамической симметрии.
Из соотношений (3.12) и (3.13) находим
Ф = TT~i' * = ш°(1 --rz-^)- (ЗЛ4)
Jx + l^z2 \ Jx + l^z2/
Если бы величина Jx + /iz2 не зависела от времени, то производные ф и ф были бы постоянными и движение несущего тела было бы регулярной прецессией. В рассматриваемом же случае движение тела является более сложным, нежели регулярная прецессия. Для него угол нутации постоянен (в — в0), а углы прецессии и собственного вращения определяются из равенств (3.14) при помощи квадратур.
4. О некоторых других частных случаях свободного движения несущего тела
Укажем еще несколько точных решений динамических уравнений (1.7), существующих при специальном выборе закона движения точки P относительно несущего тела.
П лоские движения. Пусть z = 0, то есть материальная точка P во все время движения остается в главной плоскости инерции Oxy несущего тела. Тогда система уравнений (1.7) имеет решение, в котором wx —0,wy — 0, а wz (t) удовлетворяет дифференциальному уравнению
[Jz + i(x2 + y2)]W z + i[xy — yx + 2(xx + yy)wz ] = 0. (4.1)
Это уравнение допускает интеграл
[Jz + i(x2 + y2)]wz + i(xy — yx) — const, (4.2)
и при заданном движении точки P в плоскости Oxy задача нахождения wz как функции времени сводится к квадратуре.
Стационарные вращения. Несущее тело может совершать стационарные вращения вокруг его главных центральных осей инерции, например, вокруг оси Oz:
wx — 0, wy — 0, wz — w — const. (4.3)
Подставив эти величины в уравнения (1.7), получим условия, которым должно удовлетворять относительное движение точки P для существования стационарного вращения (4.3):
zyw2 + yz — zy — 2xzw — 0, (4.4)
—xzw2 + zx — xz — 2yzw — 0, (4.5)
xy — yx + 2(xx + yy)w — 0. (4.6)
Не останавливаясь на подробном анализе этих условий, укажем только некоторые частные случаи существования стационарного вращения (4.3).
1. Случай, когда точка P движется вдоль оси Oz, то есть когда x = 0,y = 0, а z — z(t) — произвольная функция времени. Для динамически симметричного тела это стационарное вращение рассмотрено в предыдущем разделе.
2. Случай, когда точка P движется по прямой, параллельной главной оси инерции Oz, по гармоническому закону с частотой w, равной угловой скорости вращения тела:
x — x0 — const, y — y0 — const, z + w2z — 0. (4.7)
3. Случай, когда точка P совершает периодическое движение, в котором
2
х = R cos ^t, у = —R sin ^ t, z + z = 0. (4.8)
Определяемая этими соотношениями траектория точки P является замкнутой кривой, лежащей на поверхности прямого кругового цилиндра произвольного радиуса R, ось же цилиндра совпадает с осью Oz; при этом угловая скорость стационарного вращения (4.3) равна удвоенной частоте периодического движения точки P относительно несущего тела.
Представляет интерес непростая задача об устойчивости упомянутых стационарных вращений.
Список литературы
[1] Маркеев А. П. Теоретическая механика. М.-Ижевск: НИЦ «Регулярная и хаотическая динамика», Институт компьютерных исследований, 2007. 592 с.
[2] Иванов А. П. Динамика систем с механическими соударениями. М.: Международная программа образования, 1997. 336 с.
[3] Леви-Чивита Т., Амальди У. Курс теоретической механики: Т. 2, Ч. 2: Динамика систем с конечным числом степеней свободы. М.: ГИИЛ, 1951. 556с.
[4] Дубошин Г. Н. Небесная механика: Основные задачи и методы. М.: Физматлит, 1963. 587 с.
[5] Маркеев А. П. О периодическом движении твердого тела, несущего материальную точку, при наличии его соударений с горизонтальной плоскостью // Нелинейная динамика, 2012, т. 8 , № 1, с. 71-81.
[6] Малкин И. Г. Теория устойчивости движения. М.: Наука, 1966. 532 с.
On the dynamics of a rigid body carrying a material point
Anatoly P. Markeev
A.Ishlinsky Institite for Problems in Mechanics, Russian Academy of Sciences
pr. Vernadskogo 101-1, Moscow, 119526, Russia
In this paper we consider a system consisting of an outer rigid body (a shell) and an inner body (a material point) which moves according to a given law along a curve rigidly attached to the body. The motion occurs in a uniform field of gravity over a fixed absolutely smooth horizontal plane. During motion the shell may collide with the plane. The coefficient of restitution for an impact is supposed to be arbitrary. We present a derivation of equations describing both the free motion of the system over the plane and the instances where collisions with the plane occur. Several special solutions to the equations of motion are found, and their stability is investigated in some cases. In the case of a dynamically symmetric body and a point moving along the symmetry axis according to an arbitrary law, a general solution to the equations of free motion of the body is found by quadratures. It generalizes the solution corresponding to the classical regular precession in Euler’s case. It is shown that the translational motion of the shell in the free flight regime exists in a general case if the material point moves relative to the body according to the law of areas.
MSC 2010: 70E55,70F35,70K42, 70E50
Keywords: rigid body dynamics, collision, periodic motion, stability
Received March 19, 2012, accepted March 25, 2012
Citation: Rus. J. Nonlin. Dyn., 2012, vol. 8, no. 2, pp. 219-229 (Russian)