- •ПРЕДИСЛОВИЕ
- •ЧАСТЬ I. ПРОСТРАНСТВА И КООРДИНАТЫ
- •Аффинные пространства
- •Аффинные евклидовы пространства
- •Аффинные координаты и преобразования
- •Криволинейные координаты
- •Криволинейные системы координат
- •Локальные базисы
- •ЧАСТЬ II. КИНЕМАТИКА
- •Кинематика точки
- •Коэффициенты Ламе. Проекции скорости точки на оси криволинейной системы координат
- •Проекции ускорения точки на оси ортогональной криволинейной системы координат
- •Описание движения точки в естественных координатах
- •Определение кривизны траектории точки по движению
- •Два примера движения точки
- •Кинематика твердого тела
- •Движение механической системы. Твердое тело. Число степеней свободы положения. Аффинное пространство и координаты, связанные с твердым телом
- •Группа движений аффинного евклидова пространства
- •Поступательное движение твердого тела
- •Вращение твердого тела вокруг неподвижной оси
- •Плоское движение твердого тела
- •Движение твердого тела вокруг неподвижной точки
- •Угловая скорость, формула Эйлера и движение твердого тела в общем случае
- •Сложное движение
- •Сложное движение точки. Основные понятия
- •Относительная производная
- •Теорема сложения скоростей в сложном движении точки
- •Теорема сложения ускорений в сложном движении точки
- •Теорема о сложении угловых скоростей в сложном движении твердого тела
- •ЧАСТЬ III. ДИНАМИКА
- •Уравнения движения и основные законы динамики механической системы
- •Принцип детерминированности и уравнение Ньютона
- •Инерциальные системы координат
- •Сила и масса. Второй и третий законы Ньютона
- •Законы сил
- •Две задачи динамики
- •Уравнения движения механической системы
- •Теорема об изменении главного вектора количества движения
- •Уравнение движения центра инерции
- •Кинетический момент относительно неподвижной точки и теорема о его изменении
- •Движение точки в центральном поле сил
- •Изменение кинетического момента, вычисляемого относительно подвижного полюса
- •Работа силы и изменение кинетической энергии материальной точки
- •Условия потенциальности силового поля
- •Кинетическая энергия системы и теорема Кенига
- •Теорема об изменении кинетической энергии системы
- •Движение точки в центральном поле сил (продолжение)
- •Динамика твердого тела
- •Масса и плотность. Геометрия масс
- •Основные законы динамики твердого тела
- •Вращение твердого тела вокруг неподвижной точки
- •Уравнения движения свободного твердого тела
- •Динамика точки с переменной массой
- •Уравнение Мещерского
- •Две задачи Циолковского
- •ЧАСТЬ IV. АНАЛИТИЧЕСКАЯ ДИНАМИКА
- •Общее уравнение механики
- •Cвязи, реакции. Обобщенные кооpдинаты
- •Изохронные ваpиации
- •Идеальные связи. Общее уравнение механики и принцип возможных перемещений
- •Общее уравнение механики в лагранжевых координатах
- •Уравнения Лагранжа
- •Уpавнения Лагpанжа II pода, их инваpиантность
- •Разрешимость уравнений Лагранжа II рода относительно старших производных
- •Обобщенный потенциал и уравнения Лагранжа II рода
- •Уравнения Лагранжа I рода и реакции идеальных связей
- •Канонические уравнения механики
- •Вывод канонических уравнений
- •Первые интегралы канонических уравнений
- •Метод Якоби решения канонических уравнений
- •Решение задачи о движении точки в центральном поле методом Якоби
- •Вариационные принципы механики
- •Дифференциальный принцип Даламбера-Лагранжа в декартовых переменных
- •Дифференциальный принцип Даламбера-Лагранжа в канонических переменных
- •Функционал и функция действия
- •Интегральный принцип наименьшего действия при изохронном варьировании (Принцип Гамильтона)
- •ЧАСТЬ V. ДОПОЛНИТЕЛЬНЫЕ СВЕДЕНИЯ
- •Структуры и пространства
- •Группы, кольца, поля
- •Векторные пространства
- •Метрические пространства
- •Банаховы, гильбертовы и евклидовы пространства
- •Пространства Rn
- •Тензоры
- •Сопряженные пространства
- •Два определения тензора и их эквивалентность
- •Примеры тензоров
- •Алгебраические операции над тензорами
- •ЧАСТЬ VI. МАТЕРИАЛЫ ДЛЯ ЗАНЯТИЙ
- •Проверочные вопросы
- •Экзаменационные вопросы
- •Тесты
- •Упражнения
- •Литература
- •Предметный указатель
- •Оглавление
Перейти к оглавлению на странице: 256
§4. Решение задачи о движении точки в центральном поле методом Якоби
Здесь, на примере этой задачи, имеющей многочисленные и важные приложения в таких, например, разделах науки, как астрономия, физика и биология, мы обсудим схему решения уравнений движения механической системы методом Якоби. Для этого мы последовательно выпишем для нее уравнения Лагранжа второго рода, канонические уравнения, уравнение Гамильтона-Якоби, а затем рассмотрим систему равенств метода Якоби, дающих общий интеграл канонических уравнений, то есть общее решение уравнений движения точки в центральном поле.
Уравнения Лагранжа второго рода
Различные аспекты задачи о движении точки в центральном поле мы рассмотрели в §§10, 13, 16 главы 6.
Сила, действующая на материальную точку M массы m в центральном поле сил в системе координат с началом в центре сил O задается формулой
~ |
= δ · Φ(r) · |
~r |
|
(4.1) |
F |
r |
, |
а соответствующее уравнение Ньютона имеет вид (см. (10.2) гл. 6)
¨ |
~r |
|
|
m~r = δ · Φ(r) · |
|
, |
(4.2) |
r |
−−→
где ~r = OM , r = |~r|, δ = ±1, а Φ(r) модуль силы, действующей на рассматриваемую материальную точку.
Напомним, что законы сил, дающих примеры центральных силовых полей, мы рассмотрели в §4 главы 6, а в §10 главы 6 мы выписали для них соответствующие величины δ · Φ(r).
Центральное поле сил является потенциальным, причем потенциал и потенциальная энергия задаются формулами (§13 гл. 6)
Z
U(x, y, z) = (r) = δ · Φ(r)dr, Π(x, y, z) = − (r), (4.3)
167
Перейти к оглавлению на странице: 256
а тогда функция Лагранжа вычисляется по формуле:
L = T − Π = |
m v2 |
+ (r). |
(4.4) |
2 |
Так как рассматриваемая механическая система не стеснена связями, то в качестве обобщенных координат можно использовать декартовы и любые криволинейные координаты. Мы далее используем декартовы и сферические координаты.
Уравнения Лагранжа в декартовых координатах
Если в качестве обобщенных координат использовать декартовы координаты q = (q1, q2, q3) = (x, y, z), то
p
r = |q| = q12 + q22 + q32,
(4.5)
v = |q˙| = pq˙12 + q˙22 + q˙32, L = m |q˙|2 /2 + (|q|).
Если в уравнения Лагранжа
|
d ∂L |
− |
∂L |
= 0 |
(4.6) |
||
|
|
|
|
|
|||
dt |
∂ q˙ |
∂ q |
подставить выражение L по формуле (4.5) и произвести все дифференцирования слева, то в результате получим уравнения Ньютона (4.2).
Уравнения Лагранжа в сферических координатах
Введем в рассмотрение в качестве обобщенных координат точки M ее сферические координаты q = (q1, q2, q3) = (r, ϕ, ϑ), задаваемые формулами x = r sin ϑ cos ϕ, y = r sin ϑ sin ϕ, z = r cos ϑ. Пользуясь формулами (1.2), (1.5) главы 3, последовательно получаем, что
Hr = 1, Hϕ = r sin ϑ, Hϑ = r, |
(4.7) |
|||||||
vr = r,˙ vϕ = r sin ϑ ϕ,˙ |
|
˙ |
||||||
vϑ = rϑ, |
|
|||||||
L = T − Π = m (r˙ |
2 |
+ r |
2 |
2 |
2 |
2 ˙2 |
)/2 + (r). |
(4.8) |
|
|
sin ϑ ϕ˙ |
|
+ r ϑ |
Упражнение 4.1. Выпишите уравнение Лагранжа в сферических координатах для рассматриваемой задачи, подставив функцию Лагранжа в (4.6).
168
Перейти к оглавлению на странице: 256
Канонические уравнения
Канонические уравнения относительно декартовых переменных
Для того, чтобы получить канонические уравнения для рассматриваемой задачи, введем импульсы p = (p1, p2, p3) и гамильтониан H по формулам (см. (1.1), (1.4), (1.5) и (4.4), (4.5)):
p = |
∂L |
= mq,˙ |
q˙(q, p, t) = |
p |
, |
(4.9) |
|
|
|
|
|||||
∂ q˙ |
m |
H(q, p, t) = pq˙(q, p, t) − L(q, q˙(q, p, t), t) =
(4.10)
= p2/m − L(q, p/m, t) = p2/2m − (|q|).
Если в канонические уравнения
q˙ = ∂H/∂p, p˙ = −∂H/∂q |
(4.11) |
подставить H по формуле (4.10) и произвести справа дифференцирования, то в результате получим уравнения:
q˙ = |
p |
|
p˙ = δ · Φ (|q|) · |
q |
(4.12) |
|
|
, |
|
. |
|||
m |
|q| |
Эти же уравнения получаются непосредственно из уравнений Ньютона (см. (4.2)), если использовать замену переменных q = (x, y, z), p = (mx,˙ my,˙ mz˙).
Канонические уравнения относительно сферических переменных
Используя формулу (4.8) для функции L в сферических переменных, введем импульсы p = (pr, pϕ, pϑ), соответствующие координатам q = (r, ϕ, ϑ), и гамильтониан H(q, p, t):
|
∂L |
|
|
|
∂L |
|
|
2 |
|
2 |
|
|
∂L |
|
˙ |
2 |
|
|
||
pr = |
|
= mr,˙ |
pϕ |
= |
|
= mϕr˙ |
|
sin |
|
ϑ, pϑ = |
|
= mϑr |
, |
(4.13) |
||||||
|
∂ ϕ˙ |
|
|
˙ |
||||||||||||||||
|
∂ r˙ |
|
|
|
|
|
|
|
|
|
|
∂ ϑ |
|
|
|
|
|
|
||
|
|
|
pr |
|
|
|
|
pϕ |
|
|
˙ |
|
|
pϑ |
|
|
|
|
||
r˙(q, p, t) = |
m |
, ϕ˙(q, p, t) = |
mr2 sin2 ϑ |
, ϑ(q, p, t) = |
mr2 |
, |
|
(4.14) |
169
Перейти к оглавлению на странице: 256
H(q, p, t) = pq˙(q, p, t) − L(q, q˙(q, p, t), t) =
|
|
|
1 |
2 |
|
pϕ2 |
|
pϑ2 |
! − (r). |
(4.15) |
||
= |
|
|
pr |
+ |
|
+ |
|
|
||||
2m |
r2 sin2 ϑ |
r2 |
|
|||||||||
|
|
|
Выпишите |
канонические |
уравнения в |
|||||||
|
Упражнение 4.2. |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
сферических переменных для рассматриваемой задачи, подставив функцию Гамильтона в равенства (4.11).
Уравнение Гамильтона-Якоби в декартовых и сферических переменных
Для рассматриваемой задачи о движении материальной точки в центральном поле уравнение Гамильтона-Якоби в декартовых переменных (см. (3.12) и (4.10)) имеет следующий вид:
|
|
|
∂ t + 2m |
k=1 (∂S/∂ qk)2 |
! |
− |
|
q12 + q22 + q32 |
= 0, (4.16) |
||||||||||||||||
|
|
|
∂S |
|
|
1 |
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
R |
X |
|
|
|
|
|
|
|
q |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
Φ(r) dr (см. (4.3)). |
|
|
|
|
|
|
|
|
|
|||||||||
где (r) = δ · |
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
Используя формулу (4.15), в случае сферических переменных |
||||||||||||||||||||||
получаем другое уравнение Гамильтона-Якоби: |
|
! − (r) = 0. |
|||||||||||||||||||||||
|
∂ t + |
2m |
|
∂ r |
|
|
+ r2 sin2 |
ϑ |
∂ϕ |
|
|
+ r2 |
∂ϑ |
||||||||||||
∂S |
1 |
|
|
|
∂S |
|
2 |
1 |
|
|
|
∂S |
|
2 |
1 |
∂S |
|
|
2 |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(4.17) |
Общий интеграл уравнений движения
Так как уравнение Гамильтона - Якоби (4.17) явно от t и ϕ не зависит, то его полный интеграл будем искать в виде:
|
|
|
|
|
|
|
S = a1t + a2ϕ + Z(r, ϑ). |
|
(4.18) |
|||||
|
Подставляя (4.18) в (4.17), получаем: |
|
! |
− (r) = 0. |
|
|||||||||
a1 |
+ 2m |
∂ r |
|
|
+ r2 sin2 ϑ |
+ r2 |
∂ϑ |
(4.19) |
||||||
|
1 |
|
|
∂Z |
|
2 |
|
a22 |
1 |
∂Z |
|
2 |
|
|
170
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Перейти к оглавлению на странице: 256 |
|||||||||||||||||||||
|
Функцию Z будем искать в виде |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
Z(r, ϑ) = Z1(r) + Z2(ϑ), |
|
|
(4.20) |
||||||||||||||||||||||||
что, после подстановки (4.20) в (4.19) дает равенство: |
|
|
|
||||||||||||||||||||||||||||||||||||
a1 |
+ 2m |
dr |
|
|
|
|
+ r2 |
|
|
|
sin2 |
ϑ + |
dϑ |
!! |
− (r) = 0. (4.21) |
||||||||||||||||||||||||
|
1 |
|
dZ |
1 |
|
2 |
|
|
|
|
1 |
|
|
|
|
a2 |
|
|
|
|
|
|
dZ |
2 |
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|||||||||
|
Полагая здесь |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a2 |
|
|
|
+ |
dZ |
|
2 |
2 |
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
2 |
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= a3, |
|
|
(4.22) |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
sin2 ϑ |
dϑ |
|
|
|
|||||||||||||||||||||
получаем: |
|
a1 |
|
+ 2m |
dr |
|
+ r23 |
! − (r) = 0. |
(4.23) |
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
dZ |
2 |
|
a2 |
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
Система (4.22), (4.23) решается в квадратурах: |
|
|
|
|||||||||||||||||||||||||||||||||||
|
|
Z1 = e1 |
Z r |
|
|
|
|
|
|
dr, |
e1 = ±1, |
(4.24) |
|||||||||||||||||||||||||||
|
|
2m ( (r) − a1) − r23 |
|||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a2 |
|
|
|
|
|
|
|
|||||
|
|
|
|
|
Z2 = e2 |
Z s |
|
|
|
|
|
|
e2 = ±1. |
(4.25) |
|||||||||||||||||||||||||
|
|
|
|
|
a32 − sin22 |
ϑ dϑ, |
|||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Согласно (4.18), (4.20), мы получили: |
|
|
|
|
|
|
||||||||||||||||||||||||||||||||
S = a1t+a2ϕ+e1 |
|
Z r |
|
|
dr+e2 |
Z s |
|
|
|
|
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a32 |
− sin22 ϑ dϑ. |
||||||||||||||||||
|
2m ( (r) − a1) − r23 |
||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a2 |
|
|
|
|
a2 |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(4.26) |
Упражнение 4.3. Покажите, что функция (4.26) удовлетво-
ряет условиям теоремы Якоби и примените эту теорему для нахождения общего интеграла уравнений движения материальной точки в центральном поле.
171