УРАВНЕНИЯ ДВИЖЕНИЯ
1. УРАВНЕНИЯ ДВИЖЕНИЯ
1.1. Динамические и кинематические уравнения Эйлера
Рассмотрим движение свободного твердого тела относительно прямоугольной неподвижной системы координат . Согласно известной теореме Шаля, любое движение твердого тела можно рассматривать как совокупность поступательного движения, определяемого движением произвольной точки тела (полюса), и движением тела вокруг этой точки как неподвижной. Из основных теорем динамики следует, что за полюс удобно взять центр масс, поскольку в этом случае движение определяется наиболее просто. Действительно, согласно теореме о движении центра масс, последний движется как материальная точка, к которой приложены все внешние силы системы, а теоремы об изменении кинетического момента и кинетической энергии для движения относительно центра масс (относительно поступательно движущейся прямоугольной системы координат с началом в центре масс - кениговой системы координат ) формулируется точно так же, как и для движения вокруг неподвижной точки.
Пусть - масса тела, - скорость центра масс, - кинетический момент тела в его движении относительно центра масс, и - главный вектор и главный момент внешних сил относительно точки . Тогда из теоремы о движении центра масс и теоремы об изменении кинетического момента имеем два векторных дифференциальных уравнения
, (1.1)
. (1.2)
Если , , - координаты центра масс тела в неподвижной системе координат , а , , - проекции вектора на оси , , , то уравнение (1.1) запишется в виде трех скалярных уравнений
, , . (1.3)
Кинетический момент тела относительно центра масс, точки , определяется по формуле
, (1.4)
где векторы и обозначают соответственно радиус-вектор и скорость произвольной точки твердого тела с массой ( - вектор мгновенной угловой скорости). Преобразуя правую часть формулы (1.4), можно представить кинетический момент в виде произведения тензора инерции на вектор мгновенной угловой скорости /1/
. (1.5)
Введем прямоугольную систему координат , оси которой жестко свяжем с движущимся телом и расположим по его главным осям инерции для точки . Абсолютная производная связана с локальной производной (вычисляемой в подвижной системе координат ) формулой
. (1.6)
Учитывая формулу (1.6), запишем уравнение (1.2) в виде
. (1.7)
В подвижных осях координат тензор , вектор и вектор можно записать в следующем виде
,
, (1.8)
,
где постоянные , , - главные моменты инерции тела для центра ;
, , - единичные орты осей , , .
Проектируя теперь уравнение (1.7) на подвижные оси , получим три скалярных уравнения
,
, (1.9)
.
Уравнения (1.9) называются динамическими уравнениями Эйлера.
В общем случае правые части уравнений (1.3) и (1.9) могут зависеть от ориентации твердого тела в пространстве. Ориентация твердого тела в пространстве в каждый момент времени определяется положением подвижной системы координат относительно кениговой системы координат и может быть задана тремя углами Эйлера , , . На рис.1 показан угол прецессии как угол между осью и линией узлов , угол нутации - как угол между осями и , угол собственного вращения - как угол между линией узлов и осью .
Получим выражения проекций мгновенной угловой скорости твердого тела через углы Эйлера и их производные. Введем в дополнение к ортам , , осей , , орт , направленный по оси , и орт , направленный по линии узлов . Спроектируем вектор угловой скорости
(1.10)
на оси и получим следующие соотношения
,
, (1.11)
.
Из ортогональности осей следует, что
, .
Проектируя орт на ось и на плоскость , а затем на оси и , получим
, , . (1.12)
Кроме того, заметим, что
, .
После подстановки найденных выражений в соотношения (1.11) получим следующие кинематические уравнения Эйлера
,
, (1.13)
.
Систему уравнений (1.13) можно разрешить относительно производных углов Эйлера. Умножим первое уравнение на , второе на , сложим их и получим уравнение для угловой скорости прецессии . Уравнение для угловой скорости нутации получим, умножая первое уравнение на , второе на и вычитая второе уравнение из первого. Уравнение для угловой скорости собственного вращения получим из третьего уравнения, учитывая полученное выражение для . Таким образом
,
, (1.14)
.
Уравнения (1.3), (1.9), (1.13) или (1.14) образуют замкнутую систему дифференциальных уравнений, описывающую движение свободного твердого тела.
1.2. Уравнения движения в обобщенном случае Лагранжа
Рассмотрим движение свободного твердого тела, близкого к динамически симметричному, под действием медленно меняющихся во времени внешнего момента , действующего в плоскости угла нутации, и малых возмущающих моментов , направленных соответственно по осям , , (см. рис. 1). В этом случае уравнения движения (1.9) и (1.14) можно записать в виде
,
,
,
, (1.15)
,
,
.
Здесь - вектор медленных переменных (к их числу могут относиться, например, параметры движения центра масс); - отношение осевого и экваториального () моментов инерции тела относительно центра масс; - отношение внешнего момента, действующего в плоскости угла нутации, к экваториальному моменту инерции; - проекции вектора возмущающего момента на главные оси инерции тела, отнесенные к экваториальному моменту инерции; - малый параметр, характеризующий величину возмущений.
Движение твердого тела под действием нутационного момента , который имеет вид нечетного ряда Фурье по углу нутации с медленно меняющимися во времени коэффициентами, и малого возмущающего момента можно определить как движение твердого тела в обобщенном случае Лагранжа.
При и имеет место движение твердого тела в случае Эйлера, а при , и () имеет место движение тяжелого твердого тела в случае Лагранжа в классической постановке /1/.
Для ряда практических задач, например, для задачи о входе неуправляемого твердого тела в атмосферу планеты /6,7/, характерен случай, когда величина нутационного момента определяется бигармонической зависимостью от угла нутации
, (1.16)
где , - медленно меняющиеся коэффициенты, на знаки которых никаких ограничений не накладывается.
1.3. Канонические переменные Эйлера в обобщенном случае Лагранжа
Рассмотрим невозмущенное д