- •Моделирование физических процессов Учебное пособие
- •Содержание
- •3. 8. Математическая модель массообмена ……………………………. 54
- •5.1.1. Формулы конечного дифференцирования …………………… 74
- •5.6.2. Метод Ньютона ………………………………………………. 95
- •Введение
- •VII этап
- •2. Структурные математические модели
- •Основы структурного моделирования
- •2.2. Типовые элементы
- •2.2.1. Способы математического описания звеньев и систем.
- •2.2.2. Типовые динамические звенья.
- •2.3. Схемы соединения типовых элементов
- •2.4. Структурные модели теплообменного оборудования
- •2.4.1. Составление математических моделей тепловых объектов.
- •3. Физические математические модели
- •3.1 Использование математических моделей
- •Под углом 57 к горизонту со скоростью 20 м/с
- •Изменение температуры при остывании тела
- •3.2 Математическая модель теплопроводности
- •. Математическая модель гидродинамики
- •3.4 Математическая модель теплообменника
- •3.5 Математическая модель парогенератора
- •С естественной циркуляцией
- •Математическая модель горения
- •Математическая модель образования токсичных
- •3. 8. Математическая модель массообмена
- •Дифференциальные уравнения диффузионного пограничного слоя
- •Моделирование теплоотдачи диффузией.
- •3.9. Математическая модель парового котла
- •3.9.1. Математическая модель динамических процессов.
- •3.9.2. Регулирование давления пара и тепловой нагрузки
- •3.10. Пример построения математической модели объекта
- •4. Планирование эксперимента
- •Полный трехфакторный эксперимент
- •Численные методы
- •5.1. Интерполяция
- •5.1.1. Формулы конечного дифференцирования
- •5.1.2. Сплайны
- •Способы получения краевых условий:
- •5.2. Решение систем алгебраических линейных уравнений
- •Метод прогонки. Метод прогонки используется для решения систем специального вида
- •5.3. Решение систем алгебраических нелинейных уравнений
- •5.3.1. Метод итераций
- •5.3.2. Метод Ньютона
- •5.4. Решение дифференциальных уравнений
- •5.4.1. Постановка задачи. Разностные схемы.
- •5.4.2. Разностные схемы. Метод Рунге – Кутта.
- •Пример решения дифференциального уравнения методом Рунге-Кутта.
- •Результаты решения уравнения методом Рунге-Кутта
- •5.5. Решение краевых задач методом прогонки
- •5.5.2. Распространение тепла в пластине радиатора.
- •Сравнение метода прогонки с точным решением
- •5.6. Решение краевых задач методом пристрелки
- •5.6.1. Метод пристрелки.
- •5.6.2. Метод Ньютона. Рассмотрим граничную задачу, определяемую дифференциальным уравнением второго порядка
- •5.7. Решение краевых задач методом конечных разностей
- •5.7.1. Метод конечных разностей.
- •5.7.2. Линейные дифференциальные уравнения.
- •Моделирующие программы
- •Смеси газов в реакторе
- •Перемещение свободной поверхности жидкости
- •Скоростей течения испаряющегося метана
- •Рекомендованный библиографический список
- •Пасько Петр Иванович «моделирование физических процессов»
- •346360, Г. Волгодонск, ул. Ленина 73/94
3. Физические математические модели
-
3.1 Использование математических моделей
С помощью функциональных математических моделей можно исследовать самые разнообразные физические, химические, социальные и др. состояния или объекты.
Рассмотрим свободное движение тела в поле тяготения Земли с учетом сил сопротивления.
По второму закону Ньютона ускорение тела равно
.
Ускорение свободного падения g постоянно и направлено вниз по вертикали, а силы вязкого трения направлены обратно направлению скорости и равны
,
где с – коэффициент лобового сопротивления; S – площадь мидельного сечения (наибольшее сечение тела в плоскости нормальной к скорости) тела; – плотность среды; u – скорость движения тела. Коэффициент лобового сопротивления равен:
-
с = 1,11
– для диска
;
с = 1,33
– для полусферы
;
с = 0,55
– для полусферы
;
с = 0,4
– для шара
;
с = 0,045
– для капли
;
с = 0,1
– для капли
.
В качестве примера движения тела в поле тяготения Земли с учетом сил сопротивления рассмотрим движение футбольного мяча, движущегося в начальный момент времени со скоростью u под углом к горизонту.
Уравнение движения мяча имеет вид
,
,
с начальными условиями
, , ,
и с конечными условиями , ,
где k – расчетный коэффициент, ; ; с – коэффициент лобового сопротивления, с = 0,4; S – площадь сечения шара по диаметру; – плотность воздуха; u – скорость мяча.
В качестве исходных данных выберем диаметр мяча d = 0,3 м, масса мяча m = 0,7 кг, начальная скорость u = 20 м/с, угол к горизонту = 57. Чтобы реализовать данную модель в MathCad можно решать систему дифференциальных уравнений первого порядка
, ,
,
с помощью команды rkfixed (метод Рунге-Кутта IV порядка с фиксированным шагом).
Исходные данные:
Решение системы дифференциальных уравнений движения мяча по горизонтали
Решение системы дифференциальных уравнений движения мяча по вертикали
Траектория движения мяча приведена на рис. 3.1.
Рис. 3.1. Траектория движения мяча брошенного