Уравнение лагранжа 2 рода пример решения

iSopromat.ru

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

Для такой системы можно записать s уравнений, которые называются уравнениями Лагранжа второго рода или дифференциальными уравнениями движения в обобщенных координатах:

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

Уравнения Лагранжа второго рода представляют собой дифференциальные уравнения второго порядка относительно обобщенных координат q₁, q₂,…q_s.

Дважды интегрируя эти уравнения и определяя по начальным условиям постоянные интегрирования, получим систему уравнений движения в обобщенных координатах:

Уважаемые студенты!
На нашем сайте можно получить помощь по техническим и другим предметам:
✔ Решение задач и контрольных
✔ Выполнение учебных работ
✔ Помощь на экзаменах

Примеры решений дифференциальных уравнений второго порядка методом Лагранжа

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

Пример 1

Решить дифференциальное уравнение второго порядка с постоянными коэффициентами методом вариации постоянных Лагранжа:
(1)

Шаг 1. Решение однородного уравнения

Вначале мы решаем однородное дифференциальное уравнение:
(2)
Ищем решение в виде . Составляем характеристическое уравнение:

Это уравнение второго порядка.

Решаем квадратное уравнение:
.
Корни кратные: . Фундаментальная система решений уравнения (2) имеет вид:
(3) .
Отсюда получаем общее решение однородного уравнения (2):
(4) .

Шаг 2. Вариация постоянных – замена постоянных функциями

Варьируем постоянные C ₁ и C ₂ . То есть заменим в (4) постоянные и на функции:
.
Ищем решение исходного уравнения (1) в виде:
(5) .

Находим вторую производную:
.
Подставляем в исходное уравнение (1):
(1) ;

.
Поскольку и удовлетворяют однородному уравнению (2), то сумма членов в каждом столбце последних трех строк дает нуль и предыдущее уравнение приобретает вид:
(7) .
Здесь .

Вместе с уравнением (6) мы получаем систему уравнений для определения функций и :
(6) :
(7) .

Решение системы уравнений

Решаем систему уравнений (6-7). Выпишем выражения для функций и :
.
Находим их производные:
;
.

Решаем систему уравнений (6-7) методом Крамера. Вычисляем определитель матрицы системы:

.
По формулам Крамера находим:
;
.

Итак, мы нашли производные функций:
;
.
Интегрируем (см. Методы интегрирования корней). Делаем подстановку
; ; ; .

Общее решение исходного уравнения:

;
.

Пример 2

Решить дифференциальное уравнение методом вариации постоянных Лагранжа:
(8)

Шаг 1. Решение однородного уравнения

Решаем однородное дифференциальное уравнение:

(9)
Ищем решение в виде . Составляем характеристическое уравнение:

Это уравнение имеет комплексные корни:
.
Фундаментальная система решений, соответствующая этим корням, имеет вид:
(10) .
Общее решение однородного уравнения (9):
(11) .

Шаг 2. Вариация постоянных – замена постоянных функциями

Теперь варьируем постоянные C ₁ и C ₂ . То есть заменим в (11) постоянные на функции:
.
Ищем решение исходного уравнения (8) в виде:
(12) .

Далее ход решения получается таким же, как в примере 1. Мы приходим к следующей системе уравнений для определения функций и :
(13) :
(14) .
Здесь .

Решение системы уравнений

Решаем эту систему. Выпишем выражения функций и :
.
Из таблицы производных находим:
;
.

Решаем систему уравнений (13-14) методом Крамера. Определитель матрицы системы:

.
По формулам Крамера находим:
;
.

Первый интеграл немного сложней (см. Интегрирование тригонометрических рациональных функций). Делаем подстановку :

.
Поскольку , то знак модуля под знаком логарифма можно опустить. Умножим числитель и знаменатель на :
.
Тогда
.

Общее решение исходного уравнения:

.

Автор: Олег Одинцов . Опубликовано: 05-08-2013 Изменено: 19-06-2017

Maple: составление уравнений Лагранжа 2 рода и метод избыточных координат

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

С Maple (на кафедре была 6-я версия, а у лоточников домой была куплена 8-я) познакомился ещё студентом, когда начинал работать над будущей кандидатской под крылом моего первого (ныше покойного) научного руководителя. Были и добрые люди, что помогли на самом первом этапе разобраться с пакетом и начать работать.

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

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

Груз, висящий на нити длины L = 1,1 м, привязанной к гвоздю, толкнули так, что он поднялся, а затем ударился в гвоздь. Какова его скорость в момент удара о гвоздь? Ускорение свободного падения g = 10 м/с 2 .

Если не придираться к некоторонной туманности условия, то задача достаточно проста, а её решение, полученное путем довольно громоздких для школьника выкладок, в общем виде дает результат

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

Это послужило катализатором для того, чтобы взять да и откопать свои старые задумки, накопленные ещё со времен работы в оргкомитете Всероссийской Олимпиады студентов по теоретической механике — три года подряд занимался там подготовкой задач компьютерного конкурса. Задумки касались автоматизации построения уравнений движений для механических систем с неудерживающими связями и трением, используя известные всем уравнения Лагранжа 2 рода

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

Что касается Maple, то его библиотека для решения задач вариационного исчисления дает возможность быстро получить уравнения Эйлера-Лагранжа, решение которых минимизирует действие по Гамильтону, что применимо для консервативных систем

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

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

1. Метод избыточных координат

Рассматриваем механическую систему, имеющую s степеней свободы, положение которой описывается вектором обобщенных координат . Пусть также имеется r неудерживающих связей, к числу реакций которых можно причислить и трение покоя, при превышении предельного значения переходящее в активную силу трения скольжения, направление которой противоположно направлению относительной скорости скольжения.

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

Составим s + r уравнений движения в форме уравнений Лагранжа 2 рода

содержащие s+r неизвестных координат и r неизвестных реакций связей. Считая связи удерживающими, дополняем данную систему уравнениями связей (для простоты рассматривая геометрические связи) в виде

получаем замкнутую систему уравнений, из которой находятся значения реакций

являющиеся функциями первых s (независимых) обобщенных координат и скоростей и они могут быть расчитаны на любом шаге интегрирования уравнений движения (1). Для удерживающих связей типа «нить/поверхность» уравнения (1) и (2) надо дополнить условием освобождения от связи

а для связей с сухим трением вида

где F_j и N_j соответственно касательная и нормальная составляющая реакции; v_j — проекция скорости относительного проскальзывани точки приложения реакции.

Таким образом, уравнения (1) — (4) представляют собой полную математическую модель движения рассматриваемой механической системы.

Засим с теорией можно покончить и перейти к практике

2. Maple-функции построения и анализа уравнений Лагранжа

Для решения этой задачи была написана Maple-библиотека lagrange, содержащая четыре функции

LagrangeEQs — построение уравнений движения в форме Лагранжа 2 рода

В качестве входных параметров функция принимает выражение кинетической энергии T как функцию обобщенных координат и обобщенных скоростей; массив обобщенных координат q; массив радиус-векторов точек приложения сил r и массив векторов сил F.

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

Здесь надо отметить, что система уравнений геометрических связей eqs должна содержать избыточные координаты в явном виде, то есть иметь вид

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

ReduceSystem — преобразование уравнений движения с учетом уравнений связей

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

SolveAccelsReacts — решение уравнений движения относительно реакций и обобщенных ускорений

Данная функция принимает на вход систему уравнений движения eqs, преобразованную с учетом уравнений связей. Она линейна относительно вторых производных независимых координат и реакций связей. Другие входные параметры: q — вектор независимых координат; R — массив реакций, относительно которых необходимо разрешить уравнения движения.

Теперь проиллюстрируем, как применять описанное «хозяйство» в деле

3. Задача о маятнике на тонкой нерастяжимой нити

Расчетная схема будет такой. В качестве обобщенной координаты выбираем угол наклона нити к вертикали.

Поскольку нить — неудерживающая связь, нас будет интересовать её реакция, а значит введем дополнительную, избыточную координату r(t).

Приступаем. Чистим память и подключаем библиотеку линейной алгебры

Подключаем библиотеку lagrange

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

На выходе получаем выражение для кинетической энергии (для вставки сюда использована функция latex(), генерирующая результат в LaTeX-нотации)

Формируем массив сил и массив координат точек их приложения

Скармливаем всё функции LagrangeEQs()

получая на выходе уравнения движения

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

Далее задаем уравнение связи — пока нить натянута, справедливо условие

преобразуем систему с учетом этого условия и находим реакцию связи

Сила натяжения нити равна

Система (5) — (7) является полной системой уравнений движения груза, с учетом возможности провисания нити. Теперь подготовим её к численному интегрированию. Для начала разрешим её относительно ускорений, передав в SolveAccelsReacts() уравнения (5) и (6), вектор обобщенных координат и пустой массив реакций

получая на выходе

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

Готовим исходные данные и систему уравнений движения

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

Теперь проверяем «школьное» решение задачи

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

Замечу, что погрешности попадания в гвоздь — вынужденная мера: в полярных координатах, которые были использованы, задача имеет особенность, понятную из уравнения (8). Поэтому r(t) сравнивалось не с нулем, а с величиной eps достаточно малой, чтобы получить решение, и достаточно большой, чтобы численный решатель fsolve() не сходил с ума. Однако это нисколько не умаляет практической ценности изложенных результатов.

Вместо заключения

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

Тестовую версию библиотеки можно качнуть тут