Функция Лагранжа и обобщённые координаты: вывод движения

Ньютоновская механика требует расписать все силы, включая силы реакции связей, и работать с векторами в декартовых осях. Для системы с нитями, стержнями и поверхностями это быстро превращается в громоздкие проекции. Лагранжев подход меняет правила игры: вместо сил вы оперируете одной скалярной величиной - функцией Лагранжа, а движение описываете в удобных обобщённых координатах, которые сами учитывают связи. Силы реакции при этом выпадают автоматически. Ниже разберём, что такое лагранжиан, как выбрать обобщённые координаты и как из них механически получить уравнения движения. Чтобы сразу почувствовать связь между координатой и траекторией, ниже стоит интерактивный вывод для маятника.

Что такое функция Лагранжа

Функция Лагранжа (лагранжиан) механической системы - это разность кинетической и потенциальной энергий:

$$L = T - U.$$

Здесь $T$ - суммарная кинетическая энергия всех тел системы, $U$ - её потенциальная энергия. Лагранжиан - скалярная функция: у него нет направления, и записать его обычно проще, чем расписать векторный баланс сил. Вся информация о динамике системы упакована в одну формулу $L$, зависящую от координат, скоростей и, возможно, времени.

Ключевая идея в том, что движение реальной системы доставляет экстремум интегралу от лагранжиана по времени - это принцип наименьшего действия. Из него и выводится центральное уравнение метода, но на практике пользоваться им можно, даже не погружаясь в вариационное исчисление: достаточно рецепта, который мы сейчас и соберём.

Обобщённые координаты: что это и зачем

Обобщённые координаты - это любой набор независимых параметров $q_1, q_2, \dots, q_s$, который однозначно задаёт конфигурацию системы. Число $s$ равно числу степеней свободы. Главное преимущество: координаты можно выбирать под геометрию задачи, а не привязываться к декартовым осям.

Рассчитать для своей задачи

Несколько примеров, как обобщённая координата экономит силы:

• Маятник на нити. В декартовых осях нужны $x$ и $y$ плюс уравнение связи $x^2 + y^2 = \ell^2$. В обобщённой координате достаточно одного угла $\theta$ - связь уже «зашита» в выбор.
• Брусок на наклонной плоскости. Удобная координата - путь $s$ вдоль плоскости. Нормальная реакция в уравнения вообще не входит.
• Двойной маятник. Две координаты - углы $\theta_1$ и $\theta_2$. Декартовых координат было бы четыре с двумя уравнениями связей.

Связи, которые можно выразить как уравнения на координаты (нить, стержень, поверхность), называются голономными - именно для них метод работает в простейшей форме. Подробнее о подсчёте независимых параметров - в заметке про степени свободы механической системы.

Обобщённые скорости и импульсы

Производные обобщённых координат по времени $\dot q_i = dq_i/dt$ называются обобщёнными скоростями. Кинетическая энергия выражается через них, поэтому лагранжиан - это функция $L(q_i, \dot q_i, t)$.

По аналогии с обычным импульсом вводят обобщённый импульс:

$$p_i = \frac{\partial L}{\partial \dot q_i}.$$

Для угла $\theta$ это будет момент импульса, для линейной координаты - обычный импульс. Соответствующая ему «обобщённая сила» - это частная производная $\partial L / \partial q_i$. Если лагранжиан не зависит явно от какой-то координаты $q_i$ (такая координата называется циклической), то её обобщённый импульс сохраняется - это прямой путь к законам сохранения.

Уравнение Эйлера-Лагранжа

Сердце метода - уравнение Эйлера-Лагранжа. Для каждой обобщённой координаты $q_i$ оно записывается одинаково:

$$\frac{d}{dt}\left(\frac{\partial L}{\partial \dot q_i}\right) - \frac{\partial L}{\partial q_i} = 0.$$

Сколько степеней свободы - столько и таких уравнений. Это и есть уравнения движения системы. Обратите внимание: ни одной силы реакции в выводе не появляется - для голономных связей они выпадают сами. Сравните с ньютоновским разбором, где нормальную реакцию или натяжение нити приходится вводить явно и потом исключать.

Рассчитать для своей задачи

Рабочий алгоритм метода Лагранжа сводится к четырём шагам:

1. Выбрать обобщённые координаты $q_i$ по числу степеней свободы.
2. Выразить через них $T$ и $U$ и составить лагранжиан $L = T - U$.
3. Для каждой $q_i$ вычислить $\partial L/\partial \dot q_i$, $\partial L/\partial q_i$ и подставить в уравнение Эйлера-Лагранжа.
4. Решить полученные дифференциальные уравнения движения.

Полный пример: математический маятник

Покажем рецепт на классике. Возьмём грузик массой $m$ на невесомой нити длины $\ell$. Степень свободы одна, поэтому обобщённая координата - угол отклонения $\theta$ от вертикали.

Скорость грузика по дуге равна $\ell\dot\theta$, поэтому кинетическая энергия:

$$T = \frac{1}{2} m \ell^2 \dot\theta^2.$$

Высота подъёма над нижней точкой - $\ell(1 - \cos\theta)$, значит потенциальная энергия:

$$U = m g \ell (1 - \cos\theta).$$

Лагранжиан получается вычитанием:

$$L = \frac{1}{2} m \ell^2 \dot\theta^2 - m g \ell (1 - \cos\theta).$$

Теперь два нужных слагаемых уравнения Эйлера-Лагранжа:

$$\begin{aligned} \frac{\partial L}{\partial \dot\theta} &= m \ell^2 \dot\theta, \\ \frac{d}{dt}\frac{\partial L}{\partial \dot\theta} &= m \ell^2 \ddot\theta, \\ \frac{\partial L}{\partial \theta} &= - m g \ell \sin\theta. \end{aligned}$$

Подставляем в $\dfrac{d}{dt}\dfrac{\partial L}{\partial \dot\theta} - \dfrac{\partial L}{\partial \theta} = 0$ и сокращаем на $m\ell^2$:

$$\ddot\theta = -\frac{g}{\ell}\sin\theta.$$

Это точное уравнение движения маятника, полученное без единого упоминания силы натяжения нити. При малых углах $\sin\theta \approx \theta$, и оно превращается в уравнение гармонических колебаний с периодом $T_0 = 2\pi\sqrt{\ell/g}$. Калькулятор выше показывает, как при большой амплитуде реальный период уходит от этой оценки.

Рассчитать для своей задачи

Связь с энергией и законами сохранения

Метод Лагранжа естественно ведёт к законам сохранения. Если лагранжиан не зависит явно от времени, сохраняется полная энергия системы $E = T + U$ - на графике энергий в калькуляторе видно, как кинетическая и потенциальная перекачиваются друг в друга, а их сумма остаётся постоянной.

Если же $L$ не зависит от какой-то координаты $q_i$, сохраняется соответствующий обобщённый импульс $p_i$. Например, для частицы в центральном поле лагранжиан не зависит от полярного угла - отсюда сразу следует сохранение момента импульса, без отдельного вывода. Эта связь симметрий и законов сохранения формализуется теоремой Нётер, но даже на школьном уровне циклические координаты дают готовые интегралы движения.

Частые ошибки

• Путают знак в потенциальной энергии. $U$ растёт с высотой; для маятника это $+mg\ell(1-\cos\theta)$, а в лагранжиан она входит со знаком минус ($L = T - U$). Перепутанный знак ломает всё уравнение.
• Берут зависимые координаты. Обобщённые координаты должны быть независимыми. Если взять и $x$, и $y$ маятника, придётся тащить уравнение связи - теряется весь выигрыш метода.
• Забывают про полную производную по времени. В члене $\frac{d}{dt}\frac{\partial L}{\partial \dot q}$ дифференцируют по $t$ всё выражение, включая координаты внутри. Частую ошибку даёт остановка на частной производной.
• Пытаются добавить силы реакции вручную. Для голономных связей они уже исключены. Дописывать натяжение нити «для надёжности» - значит вводить лишнее.
• Считают $\partial L/\partial \dot\theta$ и $\partial L/\partial \theta$ одинаково. Это разные производные: по скорости и по координате. Их нельзя путать местами.

FAQ

Чем функция Лагранжа отличается от полной энергии?

Полная энергия - это сумма $E = T + U$, а лагранжиан - разность $L = T - U$. Энергия сохраняется (если $L$ не зависит от времени), а лагранжиан - рабочая функция, из которой выводят уравнения движения. Это разные объекты, хотя оба строятся из $T$ и $U$.

Можно ли любую механическую задачу решить через лагранжиан?

Метод универсален для систем с голономными связями и потенциальными силами. Для сил трения или неголономных связей нужны обобщения (диссипативная функция Рэлея, множители Лагранжа). Но огромный класс учебных задач - маятники, бруски, блоки, колебания - решается прямо по базовому рецепту.

Как выбрать обобщённые координаты правильно?

Берите столько независимых параметров, сколько у системы степеней свободы, и выбирайте те, в которых геометрия проще всего: углы для вращений, путь вдоль поверхности для скольжения, расстояния для радиального движения. Удачный выбор координат - половина решения: он сразу убирает связи и сокращает алгебру.

Коротко

Функция Лагранжа $L = T - U$ и обобщённые координаты - это скалярный способ описать механику без возни с силами реакции. Выбираете координаты по числу степеней свободы, составляете лагранжиан, прогоняете его через уравнение Эйлера-Лагранжа $\frac{d}{dt}\frac{\partial L}{\partial \dot q} - \frac{\partial L}{\partial q} = 0$ - и получаете уравнения движения. Циклические координаты сразу дают законы сохранения, а вся громоздкость ньютоновских проекций исчезает.