- •Математический факультет Кафедра вычислительной математики и программирования
- •Содержание
- •2 Разностные методы решения задач математической физики….…..20
- •Разностные схемы для линейных дифференциальных уравнений
- •1 Математические модели физических систем
- •1.1 Физические системы и методы их исследования
- •Основные определения и понятия
- •Системы деформируемых твёрдых тел
- •1.1.3 Подходы к исследованию систем деформируемых твёрдых тел
- •1.1.4 Средства исследования состояния систем деформируемых твёрдых тел
- •Разностные схемы для основных уравнений математической физики
- •1.2.1 Основные уравнения математической физики и их классификация
- •1.2.2 Разностная аппроксимация основных производных. По определению производная функции одной переменной записывается в виде
- •Это значит, что аппроксимация
- •Общий принцип метода сеток
- •Сходимость, аппроксимация и устойчивость разностных схем
- •2 Разностные методы решения задач математической физики
- •2.1 Разностные схемы для линейных дифференциальных уравнений
- •2.1.1 Краевые задачи для линейных дифференциальных уравнений
- •2.1.2 О задаче Дирихле
- •2.1.3 Разностные аппроксимации уравнений эллиптического типа
- •2.1.4 Разностная аппроксимация граничных условий
- •2.2 Построение разностной аппроксимации задачи Дирихле для
- •2.2.1 Разностная схема задачи Дирихле для уравнения Пуассона.
- •2.2.2 Устойчивость и сходимость разностной задачи Дирихле
- •2.3. Разрешимость разностных уравнений и способы их решения. Правило Рунге
- •2.4. Разностные схемы для линейных дифференциальных
- •2.5. Разностные схемы для линейных дифференциальных уравнений
- •2.5.1. Общая постановка задачи.
- •3. Метод конечных элементов и суперэлементов
- •3.1. Физические предпосылки методов
- •3.2 Деформации твёрдых тел
- •3. 2.1 Линейно-упругие деформации твёрдых тел
- •3.2.2. Нелинейно-упругие деформации твёрдых тел
- •3.3. Прикладные вопросы метода конечных элементов
- •3.3.1. Основная концепция метода конечных элементов
- •3.3.2 Дискретизация и оптимальная нумерация узлов области
- •3.3.3. Метод конечных элементов в теории упругости.
- •3..4 Построение конечно - элементных соотношений
- •3.5. Аналитический алгоритм метода конечных элементов для исследования
- •Свойства матрицы жёсткости
- •3.7 Аналитический алгоритм метода суперэлементов
- •Методы исследования нелинейных математических
- •4.1 Общие предпосылки методов численного исследования нелинейных
- •4.2 Итерационные численные методы
- •4.3 Безитерационные методы исследования нелинейных математических моделей деформируемых твёрдых тел и их систем
- •4.4 Сравнительный анализ эффективности методов исследования
- •5 Основы методологии компьютерного объектно-
- •5.1 Основные понятия и определения
- •5.2 Построение виртуальной физической модели системы
- •5.3 Исходные данные при компьютерном моделировании
- •5.4 Деформации грунтовых оснований при устройстве отдельных
- •6 Технология компьютерного объектно-
- •6.1 Технологические этапы компьютерного объектно-ориентированного
- •6.2 Интерфейс визуального ввода данных
- •7 Программное обеспечение компьютерного
- •7.1 Структура программного обеспечения объектно-ориентированного моделирования системы деформируемых твёрдых тел
- •7.2 Технология объектно-ориентированного программирования
- •Объектно-ориентированного моделирования физической системы
- •7.3 Программный комплекс «Энергия - 2д » моделирования
- •7.4 Программный комплекс «Энергия - ос » моделирования
- •7.5 Программный комплекс «Энергия - 3д » моделирования
- •8 Верификация технологии и программного
- •8.1 Методологические аспекты верификации
- •8.2 Верификация технологии и программного обеспечения
- •8.3 Верификация технологии и программного обеспечения
- •8.4 Верификация технологии и программного обеспечения
- •Численные методы математической физики
- •1. Метод сеток для задачи дирихле
- •1.1. Основы метода сеток
- •1.3.Варианты задания
- •Лабораторная работа №2
- •Дополнение к лабораторной работе №2
- •Лабораторная работа №3 Метод сеток для уравнения параболического типа.
- •Лабораторная работа №4 Метод сеток для уравнения гиперболического типа
2 Разностные методы решения задач математической физики
2.1 Разностные схемы для линейных дифференциальных уравнений
эллиптического типа
2.1.1 Краевые задачи для линейных дифференциальных уравнений
эллиптического типа
Исследование математическими методами стационарных процессов различной физической природы приводит к краевым задачам для уравнений эллиптического типа. К дифференциальному уравнению (1.1) эллиптического типа присоединяются граничные условия и в целом они представляют краевую задачу. Граничные условия могут быть трёх видов:
Граничные условия первого рода:
и|г= φ1(М), (2.1)
Граничные условия второго рода:
|г = φ2(М), (2.2.)
где |г - производная по внешней нормали;
Граничные условия третьего рода:
[ + u] |г = φ3(М), (2.3.)
где , , φ1, φ2, φ3 - известные функции, М – точка границы Г.
Требуется определить функцию и(х,у), которая в области D удовлетворяла бы дифференциальному уравнению (1.1) эллиптического типа и одному из граничные условий. Эта задача совместно с граничным условием (2.1) называется задачей Дирихле, совместно с граничным условием (2.2) называется задачей Неймана, а совместно с граничным условием (2.3) называется смешанной граничной задачей, иначе задачей Дирихле – Неймана.
При рассмотрении некоторых вопросов численного решения граничных задач целесообразно воспользоваться понятием и свойствами гармонических функций. Функция и(х,у) называется гармонической, если она имеет непрерывные частные производные второго порядка в области D и удовлетворяет внутри этой области уравнению Лапласа
. (2.4)
Простейшими гармоническими функциями являются линейная функция u=ax+by+c. Исходя из определения гармонических функций, задача Дирихле может быть сформулирована следующим образом: найти функцию, гармоническую внутри данной области D и принимающую на её границе Г непрерывные заданные значения.
2.1.2 О задаче Дирихле
Единственность решения задачи Дирихле и непрерывная зависимость её от краевых условий (корректность краевой задачи) вытекает из следующих свойств гармонических функций.
Свойство 1 (принцип максимума).
Гармоническая в ограниченной области функция и непрерывная в замкнутой области = D+Г, не может принимать внутри этой области значений больших, чем максимальное её значение на границе Г, и меньших, чем меньшее её значение на границе Г.
Доказательство. Пусть М – максимум значений и(х,у) на границе Г. Допустим, что функция и(х,у) в некоторой точке Р0(x0,y0) внутри D принимает значение μ = и(х0 , у0 ), причём μ > M. Составим вспомогательную функцию
v(х,у) = и(х,у)+
где d – диаметр области D. Очевидно, имеем
v(х0 ,у0) = и(х0 ,у0) = μ,
причём, при (x,y) Г выполнено неравенство
v(х,у) ≤ М+ = < μ.
Следовательно, функция v(х,у) достигает своего максимума внутри области D в некоторой точке (х0,у0), причём в этой точке будут выполнены необходимые условия для максимума функции:
Из соотношения
= 0 +
вытекает, что по меньшей мере одна из производных положительна внутри D. Поэтому функция v(х,у) ни в какой внутренней точке области D не может иметь максимума, и тем самым приходим к противоречию, следовательно v(х,у) ≤ М.
Свойство 2 (единственность задачи Дирихле).
Задача Дирихле для замкнутой и ограниченной области может иметь лишь единственное решение, т.е. не существует двух непрерывных гармонических функций в замкнутой ограниченной области D, принимающих на границе одни и те же значения.
Доказательство. Допустим, что две функции и1(х,у) и и2(х,у) гармонические в области D, совпадают всюду на её границе. Рассмотрим функцию
и(х,у)= и1(х,у)- и2(х,у).
Очевидно, что и(х,у)- гармоническая функция, обращающаяся в нуль на границе. По свойству 1 эта функция не может принимать внутри D значений больше или меньше нуля, следовательно, и(х,у) 0 внутри D и и1(х,у) = и2(х,у).
Замечание. Из доказанного не следует, что задача Дирихле для ограниченной замкнутой области D имеет решение, это свойство только утверждает, что если решение существует, то оно единственно.
Свойство 3 (корректность задачи Дирихле). Решение задачи Дирихле для ограниченной и замкнутой области D непрерывно зависит от граничных данных.
Сущность этого свойства очевидна в силу свойства 1.