- •М.Е. Семенов, н.Н. Некрасова математическое моделирование физических процессов
- •Оглавление
- •Предисловие
- •Введение
- •Общие принципы математического моделирования
- •1.1. Математическое моделирование на основе фундаментальных законов природы
- •1.2. Вариационные принципы и математические модели
- •Контрольные вопросы и задания
- •Классификация уравнений математической физики как моделей атмосферных процессов. Уравнение волновых движений
- •2.1. Понятие об общем решении уравнений в частных производных
- •2.2. Классификация уравнений
- •2.3. Волновые уравнения
- •Продольные колебания стержня
- •Контрольные вопросы и задания
- •Методы решения волновых уравнений
- •3.1. Метод Даламбера
- •3.2. Собственные колебания
- •Контрольные вопросы и задания
- •3.3. Метод Фурье
- •Контрольные вопросы и задания
- •Уравнение теплопроводности
- •4.1. Вывод уравнения теплопроводности. Первая краевая задача
- •4.2. Решение уравнения теплопроводности методом Фурье
- •4.3. Охлаждение бесконечного стержня
- •Контрольные вопросы и задания
- •5. Хаос в дискретных моделях
- •5.1. Одномерные отображения
- •5.2. Теорема Шарковского
- •5.3. Двумерные отображения, сохраняющие площадь
- •Контрольные вопросы и задания
- •Система лоренца
- •Задача о конвекции в подогреваемом снизу слое
- •6.2. Вывод уравнений Лоренца
- •Контрольные вопросы и задания
- •Динамика системы лоренца
- •Результаты численного моделирования уравнений Лоренца
- •7.2. Ограниченность и диссипативность системы Лоренца
- •Контрольные вопросы и задания
- •Неподвижные точки. Устойчивость. Бифуркация
- •. Неподвижные точки
- •Устойчивость неподвижных точек
- •Бифуркации в системе Лоренца
- •Контрольные вопросы и задания
- •Обобщенные размерности
- •9.1. Информационные размерности
- •Величину d1 называют информационной размерностью.
- •9.2. Корреляционная размерность
- •Контрольные вопросы и задания
- •Обработка реализаций. Характеристики хаотической динамики
- •10.1. Реконструкция фазового пространства. Оценка корреляционной размерности по наблюдаемой
- •10.2. Вычисление ляпуновских показателей
- •Контрольные вопросы и задания
- •Заключение
- •Библиографический список
- •3 94006 Воронеж, ул. 20-летия Октября, 84
2.2. Классификация уравнений
С помощью замены переменных уравнение второго порядка
(2.6)
сводится к одному из простейших уравнений. Полагая, что коэффициент , введем новые независимые переменные: где произвольные, но различные числа (иначе не будут взаимно независимые функции). Так как
и
то имеет место соответствие поэтому
Умножим эти вторые производные соответственно на a, 2b и c, а затем их сложим, тогда левая часть уравнения (2.6) примет вид
где
Рассмотрим вспомогательное квадратное уравнение
Его корнями являются В зависимости от значений дискриминанта возможны три случая:
D > 0 – уравнение гиперболического типа
где
D = 0 – уравнение параболического типа
D < 0 – уравнение эллиптического типа
В общем случае вводятся новые переменные – дважды непрерывно дифференцируемые функции и .
Дифференциальное уравнение называется уравнением характеристик уравнения
2.3. Волновые уравнения
В математической физике под струной понимают гибкую, упругую нить. Натяжения, возникающие в струне в любой момент времени, направлены по касательной к ее профилю. Пусть струна длины в начальный момент направлена по отрезку оси от 0 до . Предположим, что концы струны закреплены в точках и . Если струну отклонить от ее первоначального положения, а потом предоставить самой себе или, не отклоняя положение, придать в начальный момент ее точкам некоторую скорость, или отклонить струну и придать ее точкам некоторую скорость, то точки струны будут совершать движения – говорят, что струна начнет колебаться. Задача заключается в определении формы струны в любой момент времени и определении закона движения каждой точки струны в зависимости от времени.
Рассмотрим малые отклонения точек струны от начального положения. В силу этого можно предполагать, что движение точек струны происходит перпендикулярно оси OX и в одной плоскости (рис. 2.1). При этом предположении процесс колебания струны описывается одной функцией , которая дает величину перемещения точек струны с абсциссой в момент .
Рис. 2.1. Малые отклонения струны в плоскости
Так как рассматриваются малые отклонения струны в плоскости , то будем предполагать, что длина элемента струны равняется ее проекции на ось , т.е. . Также предполагается, что натяжение Т во всех точках струны одинаковое. На концах этого элемента, по касательным к струне, действует сила .
Касательные образуют с осью 0Х углы и . Тогда проекция на ось сил, действующих на элемент , будет равна Так как угол мал, можно положить , тогда имеем, что
где 0 < θ < 1 – параметр согласно формуле приращений конечных разностей Лагранжа.
Чтобы получить уравнение движения, нужно внешние силы, приложенные к элементу, приравнять силе инерции. Пусть ρ – линейная плотность струны. Тогда масса элемента струны будет . Ускорение элемента равно . Следовательно, по принципу Даламбера будем иметь . Сокращая на и обозначая получаем уравнение движения
. (2.7)
Это и есть волновое уравнение – уравнение колебаний струны. Для полного определения движения струны одного уравнения (2.7) недостаточно. Искомая функция должна удовлетворять еще граничным условиям, указывающим, что делается на концах струны ( и ), и начальным условиям, описывающим состояние струны в начальный момент ( ).