Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Основы математического моделирования и численные методы

..pdf
Скачиваний:
8
Добавлен:
12.11.2023
Размер:
4.65 Mб
Скачать

Федеральное государственное бюджетное образовательное учреждение высшего образования «Пермский национальный исследовательский политехнический университет»

Н.М. Труфанова, А.Г. Щербинин, А.В. Казаков

ОСНОВЫ МАТЕМАТИЧЕСКОГО МОДЕЛИРОВАНИЯ

И ЧИСЛЕННЫЕ МЕТОДЫ

Утверждено Редакционно-издательским советом университета

в качестве учебного пособия

Издательство Пермского национального исследовательского

политехнического университета 2018

УДК 5196(075.8) Т80

Рецензенты:

канд. техн. наук, генеральный директор В.Г. Савченко (ООО «Богословский кабельный завод», г. Краснотурьинск), канд. техн. наук, доцент В.В. Черняев

(Пермский национальный исследовательский политехнический университет)

Труфанова, Н.М.

Т80 Основы математического моделирования и численные ме­ тоды учеб, пособие / Н.М. Труфанова, А.Г. Щербинин, А.В. Казаков. - Пермь Изд-во Перм. нац. исслед. политехи, ун-та, 2018.- 112 с.

ISBN 978-5-398-02123-3

Приведены основные термины и понятия теории математиче­ ского моделирования, рассмотрено применение метода конечных элементов и метода конечных разностей для решения задач матема­ тической физики, представлены основные методы решения систем алгебраических уравнений.

Предназначено для бакалавров и магистров высших учебных заведений направлений подготовки 13.03.02, 13.04.02 «Электроэнер­ гетика и электротехника», профилей «Конструирование и технологии

вэлектротехнике» и «Управление и информационные технологии

вэлектротехнике», а также для научных и инженерных работников, аспирантов, специалистов в области инженерных наук.

УДК 5196(075.8)

SBN 978-5-398-02123-3

©ПНИПУ, 2018

СОДЕРЖАНИЕ

 

1. Теория математического моделирования.......................................................

5

1.1. Понятие математической модели.......................................................

6

1.2. Этапы решения задач........................................

8

1.3. Классификация математических моделей и их свойства................

9

1.4. Примеры построения математических моделей.

 

Детерминированные физические модели..............................................

12

Пример 1. Падение тела...................................................................

12

Пример 2. Падение тела с учетом горизонтальной скорости...

13

Пример 3. Тело, брошенное под углом к горизонту....................

15

Пример 4. Падение тела с учетом сил сопротивления...............

17

1.5. Уравнения математической физики.................................................

20

1.6. Классификация уравнений математической физики....................

21

Пример 5. Колебание струны.........................................................

23

1.7. Виды погрешностей............................................................................

27

2. Метод конечных разностей.............................................................................

28

2.1. Основные понятия..............................................................................

28

2.2. Сетка и сеточные функции................................................................

29

2.3. Аппроксимация производных, порядок аппроксимации..............

35

2.4. Метод формальной замены производных.......................................

36

2.5. Построение разностных схем...........................................................

41

2.6. Интегроинтерполяционныи метод, метод баланса.......................

44

2.7. Устойчивость и сходимость разностных схем..............................

47

2.8. Экономичные схемы...........................................................................

50

3. Метод конечных элементов.............................................................................

52

3.1. Суть метода конечных элементов....................................................

52

3.2. Дискретизация области......................................................................

54

3.3. Типы конечных элементов.................................................................

54

3.4. Одномерный симплекс-элемент........................................................

56

3.5. Двухмерный симплекс-элемент........................................................

58

3.6. Локальная система координат одномерного

 

симплекс-элемента.....................................................................................

61

3.7. Локальная система координат для двухмерного

 

симплекс-элемента.....................................................................................

62

3.8. Интерполяционные полиномы для дискретизованной

 

области..........................................................................................................

64

3.9. Применение метода Галеркина в МКЭ......................................

67

3.10. Решение стационарной одномерной задачи

 

теплопроводности МКЭ...................................................................

69

3.11. Решение одномерной осесимметричной задачи

 

электростатики МКЭ........................................................................

75

3.12. Решение стационарной двухмерной задачи

 

теплопроводности МКЭ с использованием симплекс-элементов.....

80

3.13. Решение нестационарной одномерной задачи

 

теплопроводности МКЭ...................................................................

88

4. Методы решения систем алгебраических уравнений.............................

93

4.1. Метод переменных направлений................................................

93

4.2. Метод установления..................................................................

95

4.3. Методы решения сеточных уравнений......................................

95

4.4. Прямые методы.........................................................................

97

4.4.1. Метод Гаусса (метод исключения)................................

97

4.4.2. Метод прогонки...........................................................

102

4.5. Итерационные методы..............................................................

104

4.5.1. Метод Зейделя..............................................................

104

4.5.2. Метод Гаусса - Зейделя...............................................

106

4.6. Сравнение методов..................................................................

110

Список литературы...................................................................................

111

1. ТЕОРИЯ МАТЕМАТИЧЕСКОГО МОДЕЛИРОВАНИЯ

Во второй половине XX столетия широкое распространение получила новая методология исследования сложных процессов и яв­ лений, которая основывается на методе математического моделиро­ вания и реализованных на его основе вычислительных эксперимен­ тах. Использование математических моделей для описания не оче­ видных, а предполагаемых процессов и явлений происходило

иранее. Например:

-анализ движения Урана с учетом возмущений от известных планет позволил в 1846 г. Леверье предположить, что существует планета Нептун. Он определил массу, блеск и собственное движение планеты, что через несколько месяцев было подтверждено Галле в Берлинской обсерватории;

-К.Э. Циолковский определил, что для преодоления земного притяжения достаточно первой космической скорости, а не скорости света.

Тем не менее методы математического моделирования не ис­ пользовались для исследования сложных систем и явлений (в техни­ ке, экономике, биологии и т.д.). В технической области отсутствие адекватных математических моделей приводило к развитию много­ численных инженерных методов расчета, носивших частный харак­ тер, и к господству натурных экспериментальных исследований, когда невозможность проведения полной обработки полученных ре­ зультатов приводила к неоднозначным решениям.

Положение дел существенно изменилось во второй половине XX в., когда развитие средств вычислительной техники позволило исследователям использовать современные ПЭВМ в качестве нового эффективного средства моделирования сложных технических и иных систем. Сегодня невозможно представить объекты или явле­ ния, при изучении которых не использовались бы методы математи­ ческого моделирования и вычислительного эксперимента. Разраба­ тываются и активно применяются сложные математические модели технических объектов, различных технологических, геологических

1.1. Понятие математической модели

Методы математического моделирования позволяют до разра­ ботки реального объекта (системы) рассматривать различные режи­ мы работы, вероятные процессы и выбирать оптимальные парамет­ ры воздействия для реализации необходимых свойств объекта или состояния системы. Использование вычислительных экспериментов как способа реализации математических моделей позволяет всесто­ ронне исследовать объекты, процессы или явления в целом.

Кроме того, на основе математического моделирования интен­ сивно разрабатываются системы автоматизированного проектирова­ ния, управления и обработки данных.

Основными задачами математического моделирования явля­ ются: определение законов в различных системах (природе, обще­ стве, технике) и запись их на математическом языке.

Одними из простейших математических моделей являются сле­ дующие:

- второй закон Ньютона (F = та) определяет зависимость меж­ ду массой тела т, силой, действующей на него F, и ускорением а;

- закона Ома (/= U/R) определяет зависимость между электри­ ческим напряжением цепи U, силой тока/и сопротивлением R.

Таким образом, можно сказать, что математическая модель процесса, объекта или явления - это формулировка его свойств с помощью математического языка для получения новых знаний (свойств) об изучаемом объекте с помощью формальных методов, а математическое моделирование - процесс построения и изучения математических моделей реальных процессов и явлений. Все естест­ венные и общественные науки, использующие математический ап­ парат, по сути, занимаются математическим моделированием: заме­ няют реальный объект его моделью и затем изучают последнюю.

Как и в случае любого моделирования, математическая модель не описывает полностью изучаемое явление, и возникают вопросы о применимости полученных таким образом результатов. Нужно от­ метить, что математическое моделирование - процесс не формали­ зованный, т.е. нет четких алгоритмов или рекомендаций по созда­ нию моделей.

Иным подходом в изучении процессов и явлений является эм­ пирический, недостатками которого являются зачастую высокая стоимость натурного эксперимента и получение частного решения, не позволяющего оценить явление в целом.

Пусть зависимость некоторой характеристики объекта у от па­ раметра х имеет экстремальный характер (рис. 1.1). Необходимо оп­ ределить значение функции в экстремуме. На рис. 1.1, а приведены экспериментально определенные значения функции при конкретных значениях параметра х. Видно, что наибольшее значение функция у достигает в двух точках х» и х5. Для уточнения истинного значения максимума функции у необходимо проводить дополнительные изме­ рения.

Рис. 1.1. Сравнение эмпирического и математического подхода

С другой стороны, построив математическую модель объекта (см. рис. 1.1, б) и определяя функциональную зависимость некой ха­ рактеристики объекта у от параметрах в виде у =Дх), можно опреде­ лить единственно возможное значение максимума функции, прирав­ няв к нулю ее производную.

1.2.Этапы решения задач

Внастоящее время невозможно представить себе решение крупной естественно-научной или промышленной задачи без ис­ пользования вычислительной техники. Технология исследования сложных проблем, основанная на построении и анализе с помощью ЭВМ математических моделей изучаемого объекта, называется вы­ числительным экспериментом [1].

Этапы вычислительного эксперимента могут быть представле­ ны в виде следующей схемы:

Объект

Математи­

Численный метод

Программа

Проведение

(дискрет, модель

исследо­

ческая

для ЭВМ,

вычислений и анализ

и вычислительный

вания

модель

отладка

результатов

алгоритм)

 

 

 

 

1. Производится физическая постановка задачи: описывается объект исследования, формулируются законы, которым он подчиня­ ется, отбрасываются несущественные факторы, мало влияющие на объект, обосновываются гипотезы.

2. Создается математическая модель объекта, чаще всего осно­ ванная на законах сохранения и представляющая собой систему дифференциальных (алгебраических, интегральных) уравнений.

3. Выбирается метод реализации математической модели. На этом этапе возникает необходимость использования численных методов и применения ЭВМ. Результатом использования численно­ го метода, позволяющего перейти от решения систем дифференци­ альных или интегральных уравнений к системе алгебраических, является дискретный аналог математической модели (дискретная модель), которая реализуется на ЭВМ.

4. Для реализации численного метода разрабатывается алго­ ритм решения, составляется программа для ЭВМ, производится от­ ладка программы - проверка соответствия результатов математиче­ ской модели тестовым задачам.

5. Проведение вычислений и анализ результатов. Прежде всего полученные результаты анализируются с целью подтверждения со­ ответствия математической модели исследуемому явлению, сопос­ тавляются с известными экспериментальными или теоретическими результатами. В ряде случаев при необходимости корректируется математическая модель или вычислительный процесс.

1.3. Классификация математических моделей и их свойства

Существует ряд подходов к классификации математических моделей, которые являются достаточно условными.

1.Средства построения определяют следующие возможные ви­ ды моделей:

-текстовые или словесные, которые представляют собой неко­ торое упорядоченное описание, подчиняющееся некоторой схеме (полицейский протокол, поэма Пушкина, инструкция по использо­ ванию чего-либо);

-натурные модели (модель Вселенной, модель самолета);

-абстрактные модели включают в себя математические, ком­ пьютерные, дискретные и ряд других.

2.Модели могут быть классифицированы по области приме­

нения:

-физические модели;

-биологические;

-социологические;

-экономические;

-климатические;

-экологические и т.д.

3.По уровню моделирования модели подразделяют следующим образом:

-эмпирическая модель основывается на результатах экспери­

мента;

-теоретическая модель основывается на математическом опи­ сании явлений;

-смешанная, или полуэмпирическая модель сочетает в себе первые два способа описания.

4.В зависимости от используемого математического аппарата:

-модели, включающие в себя обыкновенные дифференциаль­ ные уравнения;

-модели, записанные посредством уравнений в частных про­

изводных;

-вероятностные модели и т.д.

5.По цели моделирования:

-описательные модели используются для описания процессов или объектов, изменить или повлиять на которые человек не может,

алишь фиксирует особенности и закономерности, например модель движения астероида, кометы;

-оптимизационные модели создаются для отыскания наилуч­ шего решения задач при выполнении ряда условий и ограничений. Примером может быть известная задача о волке, козе и капусте, ко­ торые необходимо переправить через реку. В оптимизационной мо­ дели всегда присутствует ряд параметров, с помощью которых воз­ можно влиять на процесс решения и конечный результат;

-модели игровые (различные компьютерные игры);

-тренажерные (обучающие) модели;

-имитационная модель, предназначенная для определения возможных вариантов развития событий путем варьирования неко­ торых или всех параметров модели.

6. Если параметры модели зависят от времени, то различают:

- стационарные (статические) модели -

параметры модели

не зависят от времени;

 

- нестационарные (динамические) модели -

параметры модели

зависят от времени.

7. В зависимости от поведения системы модели подразделяют на дискретные, когда поведение системы описывается только в от­ дельные моменты времени, и непрерывные - непрерывно в любой момент времени.