- •Лекции по математическому моделированию
- •Математическое моделирование. Математическая модель в задачах оптимизации. Элементарные математические модели
- •Элементарные математические модели
- •Примеры моделей, получаемых из фундаментальных законов природы
- •4. Движение шара, присоединенного к пружине
- •Вариационные принципы и математические модели
- •Общая схема принципа Гамильтона.
- •Третий способ получения модели системы «шарик – пружина».
- •3. Колебания маятника в поле сил тяжести
- •4. Заключение
- •Универсальность математических моделей
- •1. Жидкость в u – образном сосуде.
- •2. Колебательный электрический контур.
- •3. Малые колебания при взаимодействии двух биологических популяций.
- •4. Заключение.
- •Сохранение массы вещества
- •1. Поток частиц в трубе.
- •2. Основные предположения о гравитационном режиме течения грунтовых вод.
- •3. Баланс массы в элементе грунта.
- •4. Замыкание закона сохранения массы.
- •5. О некоторых свойствах уравнения Буссинеска.
- •6. Основные выводы.
- •Сохранение энергии
- •1. Предварительные сведения о процессах теплопередачи.
- •2. Вывод закона Фурье из молекулярно-кинетических представлений.
- •3. Уравнение баланса тепла.
- •4. Постановка типичных краевых условий для уравнения теплопроводности.
- •5. Об особенностях моделей теплопередачи.
- •Совместное применение нескольких фундаментальных законов
- •1. Предварительные понятия газовой динамики.
- •2. Уравнение неразрывности для сжимаемого газа.
- •3. Уравнения движения газа.
- •4. Уравнение энергии.
- •Фильтрация смеси нефти и воды в пористой среде
- •Математическая модель фильтрации
- •Модель переноса примеси при однокомпонентной фильтрации
- •Модель переноса примеси при многокомпонентной фильтрации
- •Математическое моделирование физических процессов
- •1. Изменение атмосферного давления с изменением расстояния от поверхности Земли.
- •2. Задача об остывании тела.
- •3. Падение тел у земной поверхности.
- •4. Режимы течения. Вязкость. Число Рейнольдса.
- •5. Формула Стокса.
- •6. Сила гидравлического сопротивления.
- •Математическое программирование. Понятие линейного программирования. Виды задач линейного программирования. Геометрическая интерпретация задач линейного программирования
- •1. Понятие математического программирования
- •2. Понятие линейного программирования. Виды задач линейного программирования
- •3. Геометрическая интерпретация задач линейного программирования
- •1. Понятие нелинейного программирования
- •2. Классификация методов нелинейного программирования
- •2.1. Задача нелинейного программирования при ограничениях – неравенствах
- •4. Геометрическая интерпретация задач нелинейного программирования
Вариационные принципы и математические модели
Дадим упрощенную формулировку вариационного принципа Гамильтона для механической системы. На его основе выведем уравнения движения на пружине и маятника в поле сил тяжести. Сопоставим результаты получения моделей из фундаментальных законов и вариационного принципа.
-
Общая схема принципа Гамильтона.
Пусть имеется механическая система, формального или строгого определения которой пока давать не будем, имея в виду, что все взаимодействия между элементами такой системы определяются законами механики (один из простейших примеров, система «шарик – пружина»). Ведем понятие обобщенных координат , полностью определяющих положение механической системы в пространстве. Величина может быть декартовой координатой (например, координата в системе «шарик – пружина»), радиусом-вектором, угловой координатой, набором координат материальных точек, составляющих систему и т.д. Величину естественно называть обобщенной скоростью механической системы в момент времени . Набор величин и определяет состояние механической системы во все моменты времени.
Для описания механической системы вводится функция Лагранжа, которая, в простейших случаях, имеет явный смысл и записывается в следующем виде
, (1)
где и - кинетическая и потенциальная энергии соответственно.
Введем величину , называемой действием:
. (2)
Интеграл (2), очевидно, является функционалом от обобщенной координаты , т.е. функции , заданной на отрезке , он ставит в соответствие некоторое число (действие).
Принцип Гамильтона для механической системы гласит: если система движется по законам механики, то - стационарная функция для, или
. (3)
Фигурирующая в принципе наименьшего действия (3) функция - некоторая пробная функция, обращающаяся в ноль, в моменты и удовлетворяющая тому условию, что - возможная координата данной системы (в остальном произвольна).
Смысл принципа (3) в том, что из всех априори мыслимых (допускаемых) траекторий (движений) системы между моментами выбирается (реализуется) движение, доставляющее минимум функционалу действие (отсюда и происходит название принципа). Функция называется вариацией величины .
Итак, схема применения принципы Гамильтона (3) для построения моделей механических систем состоит в следующем: определяются обобщенные координаты и обобщенные скорости системы, строятся функция Лагранжа и функционал действия , минимизация которого на вариациях координаты и дает искомую модель.
-
Третий способ получения модели системы «шарик – пружина».
Воспользуемся принципов Гамильтона для построения модели движения шарика, соединенного с пружиной. В качестве обобщенной координаты системы естественно выбрать обычную эйлерову координату шарика . Тогда обобщенная скорость - обычная скорость шарика. Функция Лагранжа (1), равная , записывается через значения кинетической и потенциальной энергии системы:
.
Для величины действия получаем выражение
.
Теперь в соответствии со схемой, вычислим действия на вариациях координаты:
.
Последнюю формулу необходимо продифференцировать по (учитывая, что функции от не зависят):
,
и положить в нее :
.
Правая часть этого выражения (равного нулю в согласии с принципом Гамильтона) с помощью интегрирования ее первого члена по частям и с учетом того, что в моменты , преобразуется к виду:
.
Поскольку пробная функция , фигурирующая в формулировке принципа наименьшего действия, произвольна, то часть выражения, стоящая под знаком интеграла в квадратных скобках, должна быть равна нулю во все моменты времени :
,
т.е. движение системы должно описываться уравнением, полученным из закона Ньютона (первый способ) и закона сохранения энергии (второй способ). Все три подхода эквивалентны.