
- •М.Е. Семенов, н.Н. Некрасова математическое моделирование физических процессов
- •Оглавление
- •Предисловие
- •Введение
- •Общие принципы математического моделирования
- •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)
недостаточно. Искомая функция
должна
удовлетворять еще граничным условиям,
указывающим, что делается на концах
струны (
и
),
и начальным условиям, описывающим
состояние струны в начальный момент
(
).