
- •Г. М. Чечин
- •Оглавление
- •Часть 1. Введение в вычислительную физику...........................................................6
- •Часть 2. Элементы проблемного обучения при изучении темы
- •Часть 1. Введение в вычислительную физику модуль 1. Понятие о вычислительной физике
- •Содержание модуля 1
- •1.1. Некоторые исторические замечания
- •.1. История открытия Нептуна
- •1.1.2. Проблема Ферми-Пасты-Улама и открытие солитонов
- •1.2. Математическая модель
- •Тест рубежного контроля №1
- •Критерий оценки
- •Модуль 2. Простейшие дифференциальные уравнения
- •Содержание модуля 2
- •2.1. Движение тела под действием постоянной силы
- •2.2. Уравнение гармонического осциллятора
- •2.3. Математический маятник
- •2.4. Движение планет вокруг Солнца
- •Тест рубежного контроля №2
- •Критерий оценки
- •Модуль 3. Понятие о численных методах решения обыкновенных дифференциальных уравнений
- •Содержание модуля 3
- •. Метод Эйлера
- •3.2. О решении оду высших степеней и их систем
- •3.3. Недостатки метода Эйлера
- •3.4. Четырёхточечный метод Рунге-Кутты
- •3.5. Вычислительный эксперимент
- •Тест рубежного контроля №3
- •Критерий оценки
- •Модуль 4. Исследование периода колебаний математического маятника
- •Содержание модуля 4
- •4.1. Постановка задачи
- •4.2. Обезразмеривание задачи
- •4.3. Постановка прямого вычислительного эксперимента
- •4.4. Нахождение аналитической зависимости
- •Тест рубежного контроля №4
- •Критерий оценки
- •Часть 2. Элементы проблемного обучения при изучении темы «Динамика тел под влиянием сил гравитационного взаимодействия»
- •Модуль 5. Законы Кеплера
- •Содержание модуля 5
- •Тест рубежного контроля №5
- •Критерий оценки
- •Модуль 6. Закон всемирного тяготения Ньютона
- •Содержание модуля 6
- •Тест рубежного контроля №6
- •Критерий оценки
- •Модуль 7. Задача трех тел
- •Содержание модуля 7
- •Тест рубежного контроля №7
- •Критерий оценки
- •Модуль 8. Простые хореографии в задаче n тел Комплексная цель: Ознакомить студентов с необычными траекториями в задаче n тел с одинаковыми массами, которые получили название «простых хореографий».
- •Тест рубежного контроля №8
- •Критерий оценки
- •Литература
2.1. Движение тела под действием постоянной силы
Простейшее
дифференциальное уравнение, с которым
мы сталкиваемся в курсе физики, связано
с исследованием равноускоренного
движения материальной точки под действием
постоянной во времени силы. Очевидно,
эта ситуация описывается вторым законом
Ньютона
,
где
– постоянная сила. Поскольку ускорение
является второй производной от перемещения
по времени, это уравнение можно записать
в виде
(5)
Это уже есть дифференциальное уравнение второго порядка относительно переменной .
Поскольку
,
решение дифференциального уравнения
(5) можно получить с помощью элементарного
интегрирования. Действительно, мы имеем
,
и, стало быть, после однократного
интегрирования этого выражения получим
,
(6)
где
– некоторая произвольная постоянная.
Беря неопределенный интеграл от обеих частей последнего уравнения, находим
, (7)
где
– еще одна произвольная постоянная,
которая возникает за счет второго
интегрирования.
Формула (7)
приобретает совершенно иной вид, если
ввести более привычные обозначения для
входящих в нее постоянных величин. Пусть
есть постоянное во времени ускорение.
Константа
имеет смысл начальной координаты тела
(поскольку из (1) следует, что
)
и ее принято обозначать символом
.
Произвольная постоянная
имеет
смысл начальной скорости, поскольку из
формулы (6) ясно, что
.
В силу этого смысла,
принято обозначать символом
.
Тогда наше дифференциальное уравнение
принимает хорошо знакомый из школьного
курса физики вид
.
(8)
Это выражение есть формула для перемещения материальной точки при равноускоренном движении по прямой.
2.2. Уравнение гармонического осциллятора
Рис. 1
Рассмотрим колебания грузика массой m на пружинке с коэффициентом жесткости k, который лежит на плоском горизонтальном столе, предполагая, что трение грузика об поверхности стола отсутствует. Если грузик вывести из положения равновесия, он будет совершать колебания относительно этого положения. Эти колебания мы будем описываем зависящей от времени функцией , считая, что она определяет отклонение грузика из своего положения равновесия в момент времени t.
В горизонтальном
направлении на грузик действует только
одна сила – сила упругости пружинки,
,
определенная известным законом Гука
.
Деформация пружины является функцией времени, в силу чего, также является переменной.
Из второго закона Ньютона имеем
,
(9)
поскольку ускорение
является второй производной от смещения
:
.
Уравнение (9) можно переписать в форме
, (10)
где
.
Это уравнение получило название уравнение
гармонического осциллятора.
ЗАМЕЧАНИЕ. В математической литературе, при написании дифференциального уравнения обычно не указывают аргумент (t) около всех, зависящих от него функций. Такая зависимость предполагается по умолчанию. При использовании же математического пакета Maple в (10) необходимо указывать явную зависимость функции .
В отличие от
предыдущего примера движения тела под
действием постоянной силы в нашем случае
сила
изменяется с течением времени, и уравнение
(10) уже нельзя решить с помощью обычной
процедуры интегрирования. Попытаемся
угадать решение этого уравнения, зная,
что оно описывает некоторый колебательный
процесс. В качестве одного из возможных
решений уравнения (10) можно выбрать
следующую функцию:
.
(11)
Дифференцируя функцию (11), имеем
(12)
Подставляя выражение (12) в уравнение (10), убеждаемся, что оно удовлетворяется тождественно при любом значении t.
Однако, функция
(11) не является единственным решением
уравнения гармонического осциллятора.
Например, в качестве другого его решения
можно выбрать функцию
,
что также легко проверить аналогичным
образом. Более того, можно проверить,
что любая линейная комбинация этих двух
наугад названных решений
(13)
с постоянными коэффициентами A и B также является решениеv уравнения гармонического осциллятора.
Можно доказать, что зависящее от двух постоянных решение (13) является общим решением уравнения гармонического осциллятора (10). Это означает, что формула (13) исчерпывает все возможные решения этого уравнения. Иными словами, других частных решений, кроме тех, которые получаются из формулы (13) фиксацией произвольных постоянных А и В, уравнение гармонического осциллятора не имеет.
Заметим, что в физике наиболее часто приходится искать именно некоторые частные решения отдельных ОДУ или их систем. Рассмотрим этот вопрос более подробно.
Возбудить колебания в рассматриваемой нами системе грузика на пружинке можно разными способами. Пусть мы задали следующие начальные условия
(14)
Это значит, что в
начальный момент времени
грузик был отведен из положения равновесия
на величину a
и свободно отпущен (т.е. он начинает свое
движение с нулевой начальной скоростью).
Можно представить себе и много разных
других способов возбуждения, например,
грузику в положении равновесия «щелчком»
придается некоторая начальная скорость
и т.д. [общем случае,
].
Мы рассматриваем начальные условия (14) как некоторые дополнительные условия для выделения из общего решения (13) некоторого частного решения, соответствующего нашему способу возбуждения колебаний грузика.
Полагая t=0
в выражении (13), имеем
,
откуда следует, что B=a.
Таким образом, мы нашли одну из ранее
произвольных констант в решении (13).
Далее, дифференцируя
в формуле (13), имеем
.
(15)
Полагая в этом
выражении t=0
и учитывая второе начальное условие из
(14), получим
,
отсюда следует, что A=0
и, таким образом, исходное частное
решение имеет вид
.
(16)
Оно описывает колебательный режим рассматриваемой механической системы, который определяется условиями начального возбуждения (14).
Из школьного курса
физики известно, что в формуле (16) a
является амплитудой колебаний (она
задает максимальную величину отклонения
грузика от своего положения равновесия),
является
циклической частотой, а
– фазой колебаний (начальная фаза
оказывается при этом равной нулю).
Уравнение гармонического осциллятора (10) является примером линейного ОДУ. Это значит, что неизвестная функция и все ее производные входят в каждый член уравнения в первой степени 14. Линейные дифференциальные уравнения обладают чрезвычайно важным отличительным свойством: они удовлетворяют принципу суперпозиции. Это значит, что любая линейная комбинация двух каких либо решений линейного ОДУ также является его решением.
В рассматриваемом
нами примере уравнения гармонического
осциллятора, произвольная линейная
комбинация двух частных решений
и
является не просто каким-то новым
решением, но общим решением этого
уравнения (оно исчерпывает все возможные
его решения).
В общем случае,
это не так. Например, если бы мы имели
дело с линейным дифференциальным
уравнением третьего порядка, (т.е. если
бы в уравнение входила бы третья
производная
),
то линейная комбинация каких-либо двух
его частных решений также была бы
решением этого уравнения, но не
представляла бы собой его общее решение.
В курсе дифференциальных уравнений доказывается теорема о том, что общее решение ОДУ N-ого порядка (линейного или нелинейного) зависит от N произвольных постоянных. В случае нелинейного уравнения эти произвольные постоянные могут входить в общее решение (в отличие от (13)), нелинейным образом.
Принцип суперпозиции
играет в теории ОДУ исключительно важную
роль, поскольку с его помощью можно
построить общее решение дифференциального
уравнения в виде суперпозиции его
частных решений. Например, для случая
линейных ОДУ с постоянными коэффициентами
и их систем (уравнение гармонического
осциллятора относится именно к этому
типу уравнений) в теории дифференциальных
уравнений разработан общий метод
решения. Суть его заключается в следующем.
Ищется частное решение в виде
.
В результате его подстановки в исходное
уравнение, все зависящие от времени
множители сокращаются и мы приходим к
некоторому характеристическому
уравнению, которое для ОДУ N-ого
порядка представляет собой алгебраическое
уравнение N-ой
степени. Решая его, мы находим, тем самым,
все возможные частные решения, произвольная
линейная комбинация которых и дает
общее решение исходного ОДУ. Мы не будем
далее останавливаться на этом вопросе,
отсылая читателя к соответствующим
учебникам по теории дифференциальным
уравнениям, в которых можно найти
дальнейшие детали, в частности,
рассмотрение случая, когда характеристическое
уравнение содержит кратные корни.
Если рассматривается линейное ОДУ с переменными коэффициентами, (его коэффициенты зависят от времени), то принцип суперпозиции также справедлив, но построить в явном виде общее решение этого уравнение каким-либо стандартным методом, уже не представляется возможным. Мы вернемся к этому вопросу далее, обсуждая явление параметрического резонанса (см. раздел ??) и связанным с его исследованием уравненем Матье.