Уравнение движения из функции лагранжа

Принцип наименьшего действия в аналитической механике

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

Не все так просто

На самом деле я немного обманул, сказав, что тела всегда двигаются так, чтобы минимизировать действие. Хотя в очень многих случаях это действительно так, можно придумать ситуации, в которых действие явно не минимально.

Например, возьмем шарик и поместим его в пустое пространство. На некотором отдалении от него поставим упругую стенку. Допустим, мы хотим, чтобы через некоторое время шарик оказался в том же самом месте. При таких заданных условиях шарик может двигаться двумя разными способами. Во-первых, он может просто оставаться на месте. Во-вторых, можно его толкнуть по направлению к стенке. Шарик долетит до стенки, отскочит от нее и вернется обратно. Понятно, что можно толкнуть его с такой скоростью, чтобы он вернулся в точно нужное время.

Оба варианта движения шарика возможны, но действие во втором случае получится больше, потому что все это время шарик будет двигаться с ненулевой кинетической энергией.

Как же спасти принцип наименьшего действия, чтобы он был справедлив и в таких ситуациях? Об этом мы поговорим в следующий раз.

Так в чем же, с моей точки зрения, проблема?

Проблема в том, что автор, приводя данный пример допустил ряд фундаментальных ошибок. Она усугубляется тем, что планируемая вторая часть, со слов автора, будет опираться на эти ошибки. Руководствуясь принципом наполнения ресурса достоверной информацией я вынужден выступить с разъяснением своей позиции по данному вопросу более развернуто, и формат комментариев для этого маловат.

Данная статья расскажет о том, как строится механика на базе ПНД, и постарается объяснить читателю, что проблема, которую ставит автор цитируемой публикации отсутствует.

1. Определение действия по Гамильтону. Принцип наименьшего действия

— функция Лагранжа, для некоторой механической системы, в которой (опуская аргументы в дальнейшем) T — кинетическая энергия системы; П — потенциальная её энергия; q(t) — вектор обобщенных координат этой системы, являющийся функцией времени. при этом полагают, что моменты времени t₁ и t₂ — фиксированы.

Почему функционал, а не функция? Потому, что функция, по определению есть правило, по которому одному числу из области определения (аргументу функции) ставится в соответствие другое число из области значений. Функционал отличается тем, что качестве его аргумента выступает не число, а целая функция. В данном случае это закон движения механической системы q(t), определенный по крайней мере на промежутке времени между t₁ и t₂.

Многолетние (и это мягко сказано!) труды ученых-механиков (включая умопомрачительного Леонарда Эйлера), позволили сформулировать

Принцип наименьшего действия:

Механическая система, для которой задана функция Лагранжа , движется таким образом, что закон её движения q(t) доставляет минимум функционалу

называемому действием по Гамильтону.

Уже из самого определения ПНД следует тот факт, что данный принцип приводит к уравнениям движения лишь для ограниченного класса механических систем. Для каких? А давайте разберемся.

2. Границы применимости принципа наименьшего действия. Некоторые определения для самых маленьких

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

Рассмотрим движущуюся по траектории AB точку, к которой приложена сила . Бесконечно малое перемещение точки по траектории определяется вектором , направленным по касательной к траектории.

Элементарной работой силы на перемещении называют скалярную величину, равную

Тогда, полная работа силы на перемещении точки по траектории AB есть криволинейный интеграл

Кинетической энергией точки T называют работу, которую должны совершить приложенные к точке массой m силы, для того чтобы из состояния покоя перевести точку в движение со скоростью

В соответствии со вторым законом Ньютона

Вычислим строго стоящее под знаком интеграла скалярное произведение, для чего продифференцируем по времени скалярное произведение вектора скорости самого на себя

С другой стороны,

Дифференцируя это равенство по времени, имеем

Сравнивая (1) и (2) приходим к выводу, что

Тогда, спокойно вычисляем работу, раскрывая криволинейный интеграл через определенный, взяв в качестве пределов модуль скорости точки в начале и в конце траектории

Пусть точка перемещается в пространстве по произвольной траектории AB. Вычислим, какую работу при этом совершит сила (3)

Так как проекции силы на оси координат зависят исключительно от этих самых координат, всегда можно найти функцию

Тогда, выражение для работы преобразуется к виду

где — значения функции U(x, y, z) в точках A и B соответственно. Таким образом работа рассматриваемой нами силы не зависит от траектории точки, а определяется только значениями функции U в начале и в конце траектории. Такая сила называется консервативной силой, а соответствующая ей функция U(x, y, z) — силовой функцией. Очевидно, что , а так же равенство нулю работы консервативной силы при движении по замкнутой траектории. Говорят так же, что функция U(x, y, z) задает в пространстве силовое поле.

Потенциальной энергией точки, в пространстве с заданным силовым полем, называют работу внешних сил, приложенных к ней, которую они совершают при перемещении точки в заданное координатами (x, y, z) положение в пространстве из некоторого произвольного положения, выбранного в качестве начала отсчета уровня потенциальной энергии.

— потенциальная энергия точки в положении A, а

— потенциальная энергия точки в положении B. Учитывая всё вышесказанное снова вычислим работу потенциальных сил на перемещении из точки A в точку B

Таким образом, работа консервативных сил равна изменению потенциальной энергии точки, взятому с обратным знаком

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

3. Понятие о вариациях обобщенных координат. Постановка вариационной задачи

где s — число степеней свободы данной системы.

Действительный, но неизвестный пока нам, закон движения данной системы определяется зависимость обобщенных координат (4) от времени. Рассмотрим одну из обобщенных координат , полагая аналогичные рассуждения и для всех остальных координат.

Рисунок 1. Действительное и окольное движение механической системы

На рисунке зависимость изображена красной кривой. Выберем два произвольных фиксированных момента времени t₁ и t₂, полагая t₂ > t₁. Положение системы договоримся называть начальным положением системы, а — конечным положением системы.

Однако, я ещё раз настаиваю на том, чтобы нижеследующий текст был прочтен внимательно! Несмотря на то, что мы задаемся начальным и конечным положением системы, ни первое положение, ни второе, нам заранее неизвестны! Равно как и неизвестен закон движения системы! Эти положения рассматриваются именно как начальное и конечное положение, безотносительно конкретных значений.

Далее мы полагаем, что из начального положения в конечное система может придти разными путями, то есть зависимость может быть любой кинематически возможной. Действительное движение системы будет существовать в единственном варианте (красная кривая), остальные кинематически возможные варианты будем называть окольными движениями (синяя кривая на рисунке). Разность между действительным и окольным движением

будем называть изохронными вариациями обобщенных координат

В данном контексте вариации (5) следует понимать как бесконечно малые функции, выражающие отклонение окольного движения от действительного. Малая «дельта» для обозначения выбрана не случайно и подчеркивает принципиальное отличие вариации от дифференциала функции. Дифференциал — главная линейная часть приращения функции, вызванного приращением аргумента. В случае с вариацией изменение значения функции при постоянном значении аргумента вызвано изменением вида самой функции! Мы не варьируем аргумент, в роли которого выступает время, поэтому вариация называется изохронной. Мы варьируем правило по которому каждому значению времени приводится в соответствие некоторое значение обобщенных координат!

По сути, мы варьируем закон движения, по которому система из начального состояния перемещается в конечное состояние. Начальное и конечное состояние определяются действительным законом движения, но я ещё раз подчеркиваю — их конкретные значения нам не известны и могут быть любыми кинематически возможными, мы лишь полагаем, что они существуют и система гарантированно перемещается из одного положения в другое! В начальном и конечном положении системы мы не варьируем закон движения, поэтому вариации обобщенных координат в начальном и конечном положении равны нулю

Исходя из принципа наименьшего действия, действительное движение системы должно быть таким, чтобы доставлять минимум функционалу действия. Варьирование координат вызывает изменение функционала действия. Необходимым условием достижения функционалом действия экстремального значения является равенство нулю его вариации

4. Решение вариационной задачи. Уравнения Лагранжа 2-го рода

Загоним всё под один интеграл, и так как для вариаций справедливы все операции над бесконечно малыми величинами, преобразуем этот крокодил к виду

Исходя из определения обобщенной скорости

Тогда выражение (8) преобразуется к виду

Второе слагаемое интегрируется по частям

Исходя из условия (7), имеем

тогда, получаем уравнение

При произвольных пределах интегрирования равенство нулю определенного интеграла обеспечивается равенством нулю подынтегральной функции

С учетом того, что вариации обобщенных координат независимы, (11) справедливо только в случае равенства нулю всех коэффициентов при вариациях, то есть

Никто не мешает нам умножить каждое из уравнений на (-1) и получить более привычную запись

Уравнения (12) и есть решение задачи. И вот на этом моменте ещё раз внимание — решение вариационной задачи по принципу наименьшего действия, это не функция, доставляющая минимум действию по Гамильтону, а система дифференциальных уравнений, решая которое таковую функцию можно найти. В данном случае это дифференциальное уравнение Лагранжа 2-го рода, записанное через функцию Лагранжа, то есть в формулировке для консервативных механических систем.

И всё, на этом принцип наименьшего действия заканчивается, а начинается теория обыкновенных дифференциальных уравнений, которая, в частности, гласит, что решением уравнения (12) является вектор-функция вида

где C₁. C_2s — произвольные константы интегрирования.

ПНД — фундаментальный принцип, позволяющий получить уравнения движения системы, для которой определена функция Лагранжа

Точка! В задачах аналитической механики вышеперечисленные выкладки больше не нужно проделывать, достаточно использовать их результат (12). Функция, удовлетворяющая уравнению (12) есть закон движения системы, удовлетворяющей ПНД.

5. Задача с шариком и стенкой

Теперь вернемся к той задаче, с которой всё началось — об одномерном движении шарика около абсолютно упругой стенки. Разумеется, для данной задачи можно получить дифференциальные уравнения движения. Так как это дифференциальные уравнения движения, то любое, я подчеркиваю это, любое их решение доставляет минимум функционалу действия, а значит ПНД выполняется! Общее решение уравнений движения шарика можно изобразить в виде так называемого фазового портрета рассматриваемой механической системы. Вот этот фазовый портрет

Рисунок 2. Фазовый портрет системы в задаче с шариком

По горизонтальной оси откладывается координата шарика, по вертикальной — проекция скорости на ось x. Может это покажется странным, но данный чертеж отражает все возможные фазовые траектории движения шарика, при любых начальных, или если вам так хочется, краевых условиях. На самом деле параллельных прямых на графике бесконечно много, на чертеже показаны некоторые из них и направление движения по фазовой траектории.

Это — общее решение уравнения движения шарика. Каждая из этих фазовых траекторий доставляет минимум функционалу действия, что непосредственно следует из выкладок, проделанных выше.

Что делает автор задачи? Он говорит: вот шарик покоится, и за промежуток времени от t_A до t_B действие равно нулю. Если шарик толкнуть к стенке, то за тот же промежуток времени действие будет больше, так как у шарика отличная от нуля и неизменная кинетическая энергия. Но почему шарик движется к стенке, ведь в покое действие будет меньше? Значит ПНД испытвает проблемы и не работает! Но мы обязательно решим это в следующей статье.

То что говорит автор — бред. Почему? Да потому, что он сравнивает действия на различных ветвях одной и той же действительной фазовой траектории! Между тем, при применении ПНД, сравнивается действие на действительной траектории и на множестве окольных траекторий. То есть происходит сравнение действия на реальной траектории с действием на тех траекториях, которых нет в природе, и никогда не будет!

Не понятно? Объясню ещё более доходчиво. Рассмотрим состояние покоя. Оно описывается ветвью фазового портрета, совпадающего с осью абсцисс. Координата не меняется с течением времени. Это действительное движение. А какое же движение будет окольным. Любое другое кинематически возможное. Например малые колебания шарика около рассматриваемого нами положения покоя. Задача допускает колебания шарика вдоль оси х? Допускает, значит такое движение кинематически возможно и может рассматриваться как одно их окольных

Почему же шарик таки покоится? Да потому, что действие в состоянии покоя, вычисленное на фиксированном промежутке времени от t_A до t_B, будет меньше действия, при малых колебаниях на том же промежутке времени. Значит колебаниям и любому другому «шевелению» шарика природа предпочитает покой. В полном соответствии с ПНД.

Допустим мы толкнули шарик в сторону стенки. Пусть мы толкнули его как хочет автор, со скоростью, подобранной из краевых условий, так чтобы в момент времени t_B шарик оказался в том же положении, откуда стартовал. Шарик, с постоянной скоростью долетит до стенки, упруго отскочит и вернется в начальное положение в момент времени t_B, опять таки с постоянной скоростью. Ок, это действительное движение. Какое движение будет одним из окольных? Например, если шарик будет двигаться к стенке и от стенки со скоростью, меняющейся со временем. Такое движение возможно кинематически? Возможно. Почему же модуль скорости шарика не меняется? Да потому, что действие на такой фазовой траектории будет иметь минимальное значение, в сравнении с любом другим вариантом, где скорость зависит от времени.

Вот и всё. Ничего такого волшебного тут не происходит. ПНД работает безо всяких проблем.

Выводы и пожелания

ПНД — фундаментальный закон природы. Из него, в частности, вытекают законы механики, например дифференциальные уравнения движения (12). ПНД говорит нам о том, что природа устроена так, что уравнение движения консервативной механической системы выглядит именно как выражение (12) и никак иначе. Большего от него и не требуется.

Не нужно придумывать проблем там где их нет.

ЛАГРА́НЖА УРАВНЕ́НИЯ

ЛАГРА́НЖА УРАВНЕ́НИЯ ме­ха­ни­ки, обык­но­вен­ные диф­фе­рен­ци­аль­ные урав­не­ния вто­ро­го по­ряд­ка, опи­сы­ваю­щие дви­же­ние ме­ха­нич. сис­тем под воз­дей­ст­ви­ем при­ло­жен­ных к ним сил. Вы­ве­де­ны Ж. Ла­гран­жем в 1788 в двух фор­мах: Л. у. 1-го ро­да – урав­не­ния в де­кар­то­вых ко­ор­ди­на­тах с не­оп­ре­де­лён­ны­ми мно­жи­теля­ми Ла­гран­жа, и Л. у. 2-го ро­да – урав­не­ния в обоб­щён­ных ла­гран­же­вых ко­ор­ди­на­тах.

Теоретическая механика. Уравнения Лагранжа

В этой статье мы попробуем разобраться с такой темой, как «Уравнения Лагранжа». Вообще, уравнения Лагранжа довольно полезная штука, например, на их основе решаются задачи на малые колебания. В МГТУ им. Баумана в третьем семестре предлагается самостоятельное домашнее задание, в котором нужно записать уравнения Лагранжа для системы с двумя степенями свободы.

Итак, типовое задание выглядит так.

Прежде чем броситься решать эту задачу, посмотрим на задание и проанализируем его. Есть призма 3, которая движется поступательно по горизонтальной плоскости без трения. В призме сделан паз 2, по которому движется шарик 1. Если вы помните темы прошлого семестра, то легко увидите, что шарик совершает сложное движение — переносное поступательное вместе с призмой 3 и относительное поступательное по пазу 2. Далее есть стержень 4, который соединяет призму и каток 5. Очевидно, что скорость центра катка С равна скорость призмы. Каток движется без скольжения, это важный момент. Движение системы описывается двумя обобщенными координатами, которые любезно выбрал для нас составитель задания.

Итак, приступим к решению.

Поскольку обобщенных координат две (две степени свободы), система уравнений Лагранжа будет выглядеть так:

Расчет начинаем с записи уравнений связи — выражаем скорости всех ключевых точек и тел, имеющих массу, через обобщенные координаты. Из сказанного ранее понятно, что нам понадобится линейная скорость призмы 3, линейная скорость катка 5, угловая скорость катка 5 и скорость шарика 1. С поступательным движением все просто

С угловой скоростью катка тоже все понятно. Так как проскальзывание отсутствует

Самое трудное — выразить скорость шарика 1. Как мы уже говорили, он совершает сложное движение, значит, его скорость складывается из относительной и переносной. Переносная — это скорость поступательного движения призмы 3. Относительное — скольжение вдоль паза 2, которое описано координатой S. Значит

Векторно складываем эти две скорости

Второе выражение здесь — это теорема косинусов. Если нанести все векторы на рисунок, станет понятно, почему так.

Определившись со скоростями, записываем выражение для кинетической энергии системы Т. Полная кинетическая энергия складывается из кинетических энергий всех тел, обладающих массой. То есть в нашем случае, тел 1, 3, 5.

Шарик 1 обладает энергией

Призма 3 движется поступательно

Каток 5 совершает плоское движение, так что его кинетическая энергия складывается из энергии поступательного и вращательного движений

Полная кинетическая энергия системы

Для записи уравнений Лагранжа это выражение нужно несколько раз продифференцировать.

Сначала по координате x. Частные производные

Производную по x с точкой дифференцируем по времени

Теперь то же самое по координате S. Частные производные

Производная по времени

Левая часть уравнений Лагранжа готова. Займемся правой частью. Для нее нужно посчитать обобщенные силы по каждой координате. Есть несколько способов это сделать, мы предпочитаем делать это через элементарную работу на малом приращении координаты. В общем случае формула выглядит так

На практике это применяется следующим образом. Сначала нанесем на рисунок все действующие силы. В нашем случае это сила упругости пружины и силы тяжести.

Сначала считаем обобщенную силу по координате x. Для этого мысленно «замораживаем» координату S, и позволяем системе свободно двигаться по координате x. То есть шарик «приклеивается» к пазу 2, и внутри него никуда не движется. Все перемещение происходит по координате x. Очевидно, что сила упругости работу не совершает, так как ее длина не меняется. Очевидно, что силы тяжести работу не совершают, так как движение происходит горизонтально. Официальным языком это записывается так

Теперь обобщенная сила по координате S. Мысленно «замораживаем» координату x. Получается, что призма 3 вместе с пазом 2 и катком 5 стоит на месте, а внутри неподвижного паза движется шарик. Сила упругости совершает работу, также как и сила тяжести шарика 1. Пружина была растянута на величину статической деформации δ и дополнительно растянута на S в произвольный момент времени, то есть сила упругости равна с·(δ+S). Работа силы упругости отрицательна, так как пружина растягивается. Работа силы тяжести шарика 1 положительна, так как шарик движется вниз. Силы тяжести призмы 3 и катка 5 работу не совершают, так как эти тела покоятся. Получаем

Собственно, все. Собираем все посчитанные величины в уравнения Лагранжа и получаем систему дифференциальных уравнений, описывающих движение системы.

Для проверки можно посмотреть размерности, в обеих частях выражения размерности должны совпадать (обычно это ньютоны).

Конечно, разные задачи немного отличаются в ходе решения, но алгоритм всех задач примерно такой.

1) Определить число степеней свободы и выбрать обобщенные координаты

2) Записать уравнения связей

3) Записать выражение для кинетической энергии

4) Взять необходимые производные

5) Записать обобщенные силы по каждой координате

6) Записать уравнения Лагранжа

Если что-то не получается, не отчаивайтесь, мы всегда рады помочь.