- •ПРЕДИСЛОВИЕ
- •ОПРЕДЕЛЕНИЕ И НАЗНАЧЕНИЕ МОДЕЛИРОВАНИЯ
- •1.1. ЧТО ТАКОЕ МОДЕЛЬ?
- •Место моделирования среди методов познания
- •Определение модели
- •Определение модели
- •Цели моделирования
- •1.2. КЛАССИФИКАЦИЯ МОДЕЛЕЙ
- •Материальное моделирование
- •Идеальное моделирование
- •Когнитивные, концептуальные и формальные модели
- •1.3. Классификация математических моделей
- •Классификационные признаки
- •Классификация математических моделей в зависимости от оператора модели
- •Классификация математических моделей в зависимости от целей моделирования (рис. 1.11)
- •Классификация математических моделей в зависимости от методов реализации (рис. 1.12)
- •Задания для самостоятельного выполнения
- •Глава 2
- •ЭТАПЫ ПОСТРОЕНИЯ МАТЕМАТИЧЕСКОЙ МОДЕЛИ
- •2.1. ОБСЛЕДОВАНИЕ ОБЪЕКТА МОДЕЛИРОВАНИЯ
- •2.2. КОНЦЕПТУАЛЬНАЯ ПОСТАНОВКА ЗАДАЧИ МОДЕЛИРОВАНИЯ
- •2.3. МАТЕМАТИЧЕСКАЯ ПОСТАНОВКА ЗАДАЧИ МОДЕЛИРОВАНИЯ
- •2.4. ВЫБОР И ОБОСНОВАНИЕ ВЫБОРА МЕТОДА РЕШЕНИЯ ЗАДАЧИ
- •stop
- •2.5. РЕАЛИЗАЦИЯ МАТЕМАТИЧЕСКОЙ МОДЕЛИ В ВИДЕ ПРОГРАММЫ ДЛЯ ЭВМ
- •2.6. ПРОВЕРКА АДЕКВАТНОСТИ МОДЕЛИ
- •2.7. ПРАКТИЧЕСКОЕ ИСПОЛЬЗОВАНИЕ ПОСТРОЕННОЙ МОДЕЛИ И АНАЛИЗ РЕЗУЛЬТАТОВ МОДЕЛИРОВАНИЯ
- •Задания для самостоятельного выполнения
- •ПРИМЕРЫ МАТЕМАТИЧЕСКИХ МОДЕЛЕЙ
- •3.1. СТАТИЧЕСКИЙ АНАЛИЗ КОНСТРУКЦИЙ
- •Содержательная постановка задачи
- •Концептуальная постановка задачи
- •Математическая постановка задачи
- •Решение задачи
- •Анализ результатов
- •Концептуальная постановка
- •Математическая постановка задачи
- •Методика решения задачи
- •Анализ результатов
- •3.3. ДИНАМИКА ПОПУЛЯЦИИ
- •Содержательная постановка задачи
- •Концептуальная постановка задачи
- •Решение задачи
- •Анализ результатов
- •Математическая постановка задачи для модели Ферхюльста
- •Решение задачи
- •Анализ результатов
- •Численное исследование модели Ферхюльста
- •3.4. МОДЕЛЬ КОНКУРЕНЦИИ ДВУХ ПОПУЛЯЦИЙ
- •Математическая постановка задачи
- •Качественный анализ задачи
- •Численное исследование модели конкуренции популяций
- •3.5. ГАРМОНИЧЕСКИЙ ОСЦИЛЛЯТОР
- •Содержательная постановка задачи
- •Концептуальная постановка задачи
- •Математическая постановка задачи
- •Решение задачи
- •Качественный анализ задачи
- •Численное исследование модели
- •Качественный анализ задачи
- •Решение задачи
- •Задания для самостоятельного выполнения
- •Глава 4
- •СТРУКТУРНЫЕ МОДЕЛИ
- •4.1. ЧТО ТАКОЕ СТРУКТУРНАЯ МОДЕЛЬ?
- •4.2. СПОСОБЫ ПОСТРОЕНИЯ СТРУКТУРНЫХ МОДЕЛЕЙ
- •4.3. ПРИМЕРЫ СТРУКТУРНЫХ МОДЕЛЕЙ
- •Задания для самостоятельного выполнения
- •МОДЕЛИРОВАНИЕ В УСЛОВИЯХ НЕОПРЕДЕЛЕННОСТИ
- •5.1. ПРИЧИНЫ ПОЯВЛЕНИЯ НЕОПРЕДЕЛЕННОСТЕЙ И ИХ ВИДЫ
- •5.2. МОДЕЛИРОВАНИЕ В УСЛОВИЯХ НЕОПРЕДЕЛЕННОСТИ, ОПИСЫВАЕМОЙ С ПОЗИЦИЙ ТЕОРИИ НЕЧЕТКИХ МНОЖЕСТВ
- •5.3. МОДЕЛИРОВАНИЕ В УСЛОВИЯХ СТОХАСТИЧЕСКОЙ НЕОПРЕДЕЛЕННОСТИ
- •5.4. МОДЕЛИРОВАНИЕ МАРКОВСКИХ СЛУЧАЙНЫХ ПРОЦЕССОВ
- •Задания для самостоятельного выполнения
- •ЛИНЕЙНЫЕ И НЕЛИНЕЙНЫЕ ЯВЛЕНИЯ
- •6.1. О ЗАКОНЕ ГУКА И ГРАНИЦАХ ЛИНЕЙНОСТИ
- •6.3. О ПОСТРОЕНИИ СПЛОШНЫХ МОДЕЛЕЙ. ВЫВОД ВОЛНОВОГО УРАВНЕНИЯ ИЗ ЗАКОНОВ МЕХАНИКИ
- •6.4. РЕШЕНИЕ ВОЛНОВОГО УРАВНЕНИЯ МЕТОДОМ ФУРЬЕ
- •6.6. УРАВНЕНИЯ МАКСВЕЛЛА
- •6.7. О КЛАССИФИКАЦИИ КВАЗИЛИНЕЙНЫХ СИСТЕМ
- •6.8. СВЯЗЬ НЕПРЕРЫВНОГО И ДИСКРЕТНОГО НА ПРИМЕРАХ УРАВНЕНИЯ КОЛЕБАНИЙ СТРУНЫ И УРАВНЕНИЯ ШРЕДИНГЕРА
- •6.9. О ПОЛЬЗЕ ФЕНОМЕНОЛОГИИ ПРИ ПОСТРОЕНИИ МАТЕМАТИЧЕСКИХ МОДЕЛЕЙ
- •6.10. АНАЛИЗ ПОДОБИЯ И РАЗМЕРНОСТИ
- •6.11. АВТОМОДЕЛЬНОСТЬ
- •6.12. САМООРГАНИЗАЦИЯ И СТРУКТУРЫ В НЕЛИНЕЙНЫХ СРЕДАХ
- •6.13. О НЕЛИНЕЙНЫХ ВОЛНАХ В СПЛОШНЫХ СРЕДАХ
- •6.14. ИЕРАРХИЧЕСКИЕ МОДЕЛИ ТУРБУЛЕНТНОСТИ И МНОГОМАСШТАБНЫЕ ФУНКЦИОНАЛЬНЫЕ БАЗИСЫ
- •Иерархический базис для турбулентных полей
- •Одномерный иерархический базис
- •Двумерный базис
- •6.15. ВЕЙВЛЕТЫ
- •Непрерывное вейвлет-преобразование
- •6.16. ВЕЙВЛЕТ-АНАЛИЗ ВРЕМЕННЫХ КОЛЕБАНИЙ
- •6.17. О ФРАКТАЛАХ И ИХ ПРИМЕНЕНИИ
- •Примеры фракталов
- •Подобие и скейлинг
- •Множества Мандельброта и Жюлиа
- •Фрактальная размерность кластеров
- •Экспериментальные методы определения фрактальной размерности
- •Модель случайных фракталов для описания растущих дендритных структур
- •Результаты применения модели случайных фракталов
- •6.18. НЕЛИНЕЙНЫЕ МОДЕЛИ ДНК
- •Структура и физические свойства ДНК
- •Модель Пейрара-Бишопа
- •Задания для самостоятельного выполнения
- •Глава 7
- •7.2. ИМИТАТОР СИСТЕМЫ МАССОВОГО ОБСЛУЖИВАНИЯ
- •7.3. КЛЕТОЧНЫЕ АВТОМАТЫ
- •7.4. МОДЕЛИРОВАНИЕ ДИСЛОКАЦИЙ В МЕТАЛЛЕ
- •Самоорганизация дислокаций в модели клеточных автоматов
- •Задания для самостоятельного выполнения
- •ЯЗЫК ФОРМАЛЬНОГО ОПИСАНИЯ АЛГОРИТМОВ
- •====== Приложение 2
- •П2.1. Решение уравнений высоких степеней и трансцендентных уравнений с одним неизвестным
- •П2.2. Решение систем линейных уравнений
- •БИБЛИОГРАФИЧЕСКИЙ СПИСОК
- •ПРЕДМЕТНЫЙ УКАЗАТЕЛЬ
- •ОГЛАВЛЕНИЕ
- •Введение в математическое моделирование
= = = = = = = = |
Глава 3 |
ПРИМЕРЫ МАТЕМАТИЧЕСКИХ МОДЕЛЕЙ
—Должна же я стерпеть двух-трех гусениц, — сказала Роза, — если хочу познакомиться с бабоч ками. Они, должно быть, прелестны.
А. де Сент-Экзюпери
Настоящая глава посвящена иллюстрации рассмотренных выше этапов разработки математических моделей на простых примерах, часть которых известна из различных предметов школьного курса. Подобный выбор примеров обусловлен стремлением с достаточной (даже для школьников старших классов и студентов младших кур сов) прозрачностью показать именно «технологию» создания мате матических моделей и в то же время обратить внимание читателей на некоторые нюансы, присущие большинству из них. Практичес ки все приведенные здесь результаты могут быть получены само стоятельно теми, на кого ориентировано предлагаемое пособие. При этом следует отметить, что разработка математических моделей ре альных процессов и явлений существенно сложнее, чем в приме рах данной главы, и требует весьма глубоких знаний, интуиции, догадки.
Как уже отмечалось, аналитические модели позволяют получить искомое решение в виде конечной (или счетной) комбинации эле ментарных функций, что дает возможность более полно изучить свойства объекта моделирования, его качественное поведение, при меняя хорошо развитые методы исследования аналитических фун кций.
В тех случаях, когда по тем или иным причинам приходится ограничиваться некоторым конечным «набором» аналитических выражений (например, несколькими членами ряда, если решения
представляют собой бесконечные ряды), не являющимся точным ре шением задачи, модель получается приближенной, допускающей определенные погрешности в результатах. Однако и в этом случае сохраняется возможность качественного анализа ее поведения.
Следует отметить, что в настоящее время аналитические и при ближенные решения имеются для весьма небольшого числа задач, представляющих практический интерес. Как правило, прикладные научные и технические проблемы требуют учета такого громадно го числа параметров и связей, в том числе нелинейных, что полу чение для них аналитических или приближенных решений на со временном уровне развития математики не представляется возмож ным.
При численном подходе математические соотношения модели заменяются некоторым приближенным аналогом, что достигается дискретизацией принятых уравнений и в итоге заменой непрерыв ного решения конечным для ограниченного числа значений аргу ментов.
Благодаря развитию средств вычислительной техники числен ные методы решения сложных прикладных задач получили широ кое распространение. Задачи теплопроводности, механики жидко стей и газов, механики деформируемого твердого тела и многие другие были решены в основном благодаря широкому использова нию сеточных* методов. Вместе с тем следует иметь в виду, что неквалифицированное применение разностных схем к решению уравнений, описывающих эволюционные процессы, приводит к получению решений, далеких от истинных. Поэтому понятен ин терес к теории численных методов, позволяющей еще на стадии разработки алгоритма решения сложной инженерной задачи выяс нить условия успешной реализации вычислительной процедуры модели. Основная идея решения задач численными методами на ЭВМ заключается в сведении исходной задачи к решению системы алгебраических (линейных или нелинейных) уравнений. При этом, естественно, возникает вопрос о разрешимости этой системы. Чис
* В более широком смысле следует говорить о проекционных мето дах отыскания решения операторного уравнения в заданном подпростран стве, основанных на проектировании искомого решения уравнения на некоторое (вообще говоря, другое, как правило, —конечномерное) под пространство. Проекционные методы являются основой построения различ ных вычислительных схем решения краевых задач. К их числу относятся такие широко известные методы, как методы Галеркина, Ритца, конеч ных элементов, наименьших квадратов, конечно-разностные, коллокаций и ряд других.
ленное решение сходится к точному, если при неограниченном увеличении числа алгебраических уравнений решение дискретизи рованной задачи стремится к решению исходной задачи. Посколь ку решить систему с бесконечным числом алгебраических уравне ний невозможно, весьма актуальным представляется вопрос об оценке погрешности получаемых численных решений исходной задачи.
3.1. СТАТИЧЕСКИЙ АНАЛИЗ КОНСТРУКЦИЙ
В качестве первого примера рассмотрим порядок построения простейшей математической модели, описывающей состояние кон струкции в условиях равновесия. Необходимость в построении по добных моделей возникает как при анализе прочности и устойчи вости конструкции в целом, так и при анализе усилий (или на пряжений), возникающих в отдельных ее элементах. Если моделирование исследуемого объекта ведется в условиях равнове сия, то время в качестве отдельного параметра явно в модели не учитывается. Задачи подобного типа относят к задачам статики, а анализ состояния объекта, проведенный с использованием моделей данного типа, называют статическим. Рассмотрим статический анализ однопролетного арочного моста (рис. 3.1).
Рис. 3.1. Конструкция однопролетного арочного моста
Содержательная постановка задачи
Арка моста состоит из двух одинаковых половинок, соединен ных между собой с помощью одного шарового шарнира Е, рас положенного на средней линии моста (рис. 3.1,6). К опорам арка крепится с помощью четырех цилиндрических шарниров А, В, С и D. Подобный способ крепления позволяет половинкам моста по ворачиваться друг относительно друга, компенсируя изменение своих размеров при сезонных колебаниях температуры.
Известны геометрические размеры элементов конструкции арки моста: длина моста (расстояние между опорами А и В) L = 14 м; ширина моста w = 4 м; высота арки моста h = 6 м. Обе половинки моста имеют одинаковый вес G= GY= G2= 120 кН. Центр тяжести половинок моста расположен в продольной вертикальной плоско сти симметрии моста. Линия действия силы тяжести G половинки арки моста пересекает плоскость X Z (рис. 3.1,6) на расстоянии а = 1,5 м от оси Z
Модель должна позволять определять усилия в узлах крепле ния арки моста к опорам.
Концептуальная постановка задачи
Статический анализ конструкции проведем при следующих до пущениях:
>объектом исследования является арка моста, состоящая из двух одинаковых половинок;
>арка крепится к опорам моста с помощью четырех цилинд рических шарниров, оси которых параллельны оси Z и ле жат в плоскости XZ;
>половинки арки моста соединены между собой шаровым шар ниром, расположенным в точке Е вертикальной плоскости продольной симметрии моста;
>все элементы конструкции считаются абсолютно твердыми телами, т.е. пренебрегаем их деформациями при нагружении внешними усилиями;
>пренебрегаем наличием люфта и трения в шарнирах А, В, С, D и Е.
Сучетом сформулированных допущений требуется определить величины реакции опор RA, RB, Rc, RD арки моста.
Тела, препятствующие свободному перемещению исследуемой конструкции, называются связями. В данном случае в качестве свя зей выступают опоры арки моста. При построении модели стати-