Введение в уравнение Лагранжа

Аналитическая механика выводит уравнения движения из одной функции — Лагранжиана.

Уравнение Лагранжа второго рода — уравнение движения через функцию Лагранжа $L = T - \Pi$: $\frac{d}{dt}\frac{\partial L}{\partial \dot q} - \frac{\partial L}{\partial q} = 0$.

Метод Ньютона требует расписывать все силы и реакции связей — для сложных систем это громоздко. Аналитическая механика Лагранжа предлагает иной путь: описать систему обобщёнными координатами $q$ (любыми удобными параметрами — углами, длинами) и составить одну функцию — Лагранжиан:

$$L = T - \Pi,$$

разность кинетической и потенциальной энергии. Уравнение движения тогда получается механически, по единой формуле (уравнение Лагранжа второго рода):

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

Гениальность подхода в том, что реакции идеальных связей выпадают сами — не нужно вводить неизвестные реакции опор.

Выведем уравнение осциллятора

Проверим метод на знакомом пружинном маятнике. Координата $q = x$, кинетическая энергия $T = \frac{m\dot x^2}{2}$, потенциальная $\Pi = \frac{c x^2}{2}$. Посчитаем производные численно и убедимся, что получается $m\ddot x + c x = 0$.

m = 1.0
c = 100.0
h = 1e-6

# L(x, v) = T - П = 0.5*m*v^2 - 0.5*c*x^2
def Lagr(x, v):
    return 0.5 * m * v**2 - 0.5 * c * x**2

x0, v0 = 0.03, 0.4
# dL/dv (частная производная по скорости)
dL_dv = (Lagr(x0, v0 + h) - Lagr(x0, v0 - h)) / (2 * h)
# dL/dx (частная производная по координате)
dL_dx = (Lagr(x0 + h, v0) - Lagr(x0 - h, v0)) / (2 * h)

print(f"dL/dv = {dL_dv:.4f}  (теория m*v = {m*v0:.4f})")
print(f"dL/dx = {dL_dx:.4f}  (теория -c*x = {-c*x0:.4f})")
print("уравнение: d/dt(m*v) - (-c*x) = m*a + c*x = 0")

Вывод:

dL/dv = 0.4000  (теория m*v = 0.4000)
dL/dx = -3.0000  (теория -c*x = -3.0000)
уравнение: d/dt(m*v) - (-c*x) = m*a + c*x = 0

$\frac{\partial L}{\partial \dot x} = m\dot x$, его производная по времени $m\ddot x$; $\frac{\partial L}{\partial x} = -c x$. Уравнение Лагранжа даёт $m\ddot x + c x = 0$ — в точности уравнение осциллятора, полученное теперь без единого упоминания сил.

Вершина классической механики

Аналитическая механика Лагранжа (и развивший её формализм Гамильтона) — это вершина и завершение классической механики, плод работы математиков XVIII–XIX веков. Её красота в универсальности: один и тот же рецепт — составить $L = T - \Pi$ и применить уравнение Лагранжа — решает задачи любой сложности, от маятника до движения планет, не требуя изобретать подход заново для каждой системы. Но значение её гораздо шире механики. Принцип наименьшего действия, лежащий в основе метода, оказался настолько фундаментальным, что на нём построена вся современная физика: квантовая механика (интегралы по траекториям Фейнмана), теория поля, общая теория относительности — все формулируются через лагранжианы и действие. Изучая скромное уравнение Лагранжа на примере пружинки, вы прикасаетесь к идее, которая описывает Вселенную на самом глубоком уровне. Это достойное завершение курса теоретической механики: от простого вектора силы в первом уроке мы пришли к принципу, из которого выводится движение всего сущего.

Почему это важно

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

Как работает под капотом

Откуда берётся такая «магическая» формула? Она следует из принципа наименьшего действия: реальное движение системы между двумя моментами времени таково, что интеграл $S = \int L\,dt$ (действие) минимален. Условие минимума функционала (уравнение Эйлера—Лагранжа из вариационного исчисления) и есть $\frac{d}{dt}\frac{\partial L}{\partial \dot q} = \frac{\partial L}{\partial q}$. Член $\frac{\partial L}{\partial \dot q} = m\dot q$ — это обобщённый импульс, а $\frac{\partial L}{\partial q} = -\frac{\partial \Pi}{\partial q}$ — обобщённая сила. Так уравнение Лагранжа оказывается переформулировкой $\frac{dp}{dt} = F$ — того же второго закона Ньютона, но в координатах, удобных для конкретной задачи, и со скрытыми реакциями связей.

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

Брать Лагранжиан как $T + \Pi$ вместо $T - \Pi$.
Путать частную производную по $q$ и по $\dot q$ — это разные операции.
Забывать брать полную производную по времени от $\partial L/\partial \dot q$.
Включать в потенциальную энергию работу реакций идеальных связей — они не входят.

Итог

Аналитическая механика описывает систему обобщёнными координатами и Лагранжианом $L = T - \Pi$.
Уравнение Лагранжа: $\frac{d}{dt}\frac{\partial L}{\partial \dot q} - \frac{\partial L}{\partial q} = 0$.
Реакции идеальных связей выпадают автоматически — огромная экономия.
Метод вытекает из принципа наименьшего действия и эквивалентен законам Ньютона.