- •Модуль 1 Основы теории разностных схем
- •Комплексная цель модуля
- •Введение
- •1.3. Сетки и сеточные функции. Вводные понятия
- •1.4. Построение разностных схем
- •1.4.1. Интегро-интерполяционный метод
- •1.4.2. Консервативные и не консервативные разностные схемы
- •1.4.3. Исследование аппроксимации дифференциального уравнения
- •1.4.4 . Исследование погрешности аппроксимации граничного условия
- •1.4.5. Понятие о погрешности аппроксимации, устойчивости и сходимости разностной схемы на примере краевой задачи для оду второго порядка
- •1.4.6. Разностные операторы и их свойства
- •1. Разностные тождества
- •1.1. Формулы разностного дифференцирования произведения
- •1.2. Первая разностная формула Грина
- •1.3. Вторая разностная формула Грина
- •2. Оператор разностной производной второго порядка и его собственные значения
- •3. Обобщение оценок для скалярных произведений, содержащих оператор второй разностной производной с переменными коэффициентами
- •1.4.7. Построение разностных схем для уравнения теплопроводности
- •1.4.8. Связь аппроксимации, устойчивости, сходимости для линейных разностных задач
- •1.4.9. Исследование устойчивости разностных задач методом гармоник
- •1.4.10. Неявная схема для уравнения теплопроводности
- •1.4.11. Схема с весами
- •1.4.12. Трехслойные разностные схемы
- •Исследование устойчивости методом гармоник однородного уравнения с нулевыми граничными условиями. Необходимое условие устойчивости
- •1.4.13. Примеры трехслойных разностные схем для уравнения теплопроводности
- •1.5. Принцип максимума и следствия из него
- •1.5.1. Аппроксимация задачи Дирихле для уравнения Пуассона в прямоугольнике
- •1.5.2. Каноническая форма сеточных уравнений
- •1.5.3. Каноническая форма общего вида для сеточного уравнения
- •Разностная схема, аппроксимирующая задачу Дирихле для уравнения Пуассона в прямоугольнике:
- •Принцип максимума
- •1.5.4. Следствия из принципа максимума
- •1.5.5. Теорема сравнения . Мажоранта
- •1.5.6. Следствие из теоремы сравнения
- •1.5.7. Оценка решения неоднородного уравнения
- •1.5.8. Исследование устойчивости разностных схем на основе принципа максимума и следствий из него
- •1.5.9. Примеры решения задач на принцип максимума и следствий из него
- •1.6. Энергетический метод исследования устойчивости
- •1.6.1. Некоторые сведения из теории операторов
- •Знакоопределенные операторы
- •1.6.2. Свойства оператора симметричной разностной производной второго порядка в гильбертовом пространстве
- •Свойства разностных эллиптических операторов
- •Построение энергетической эквивалентности оператора второй разностной производной с переменными коэффициентами
- •Корректность операторно-разностных уравнений
- •1.7.1. Общие сведения
- •1.7.2. Краткие сведения об обратимости оператора в линейных нормированных пространствах
- •1.6.3. Свойства операторов разностных производных первого порядка
- •Каноническая форма двухслойных схем и ее применение к исследованию устойчивости
- •1.7.4. Необходимое и достаточное условие устойчивости в на
- •Энергетическое пространство
- •1.7.5. Исследование устойчивости в случае несамосопряженного оператора
- •Экономичные схемы для нестационарных задач математической физики
- •Схемы переменных направлений
- •Факторизованные схемы
- •Понятие о монотонных схемах
- •Понятие асимптотической устойчивости.
- •Понятие о трехслойных разностных схемах.
- •Введение в метод конечных элементов
- •Кусочно-линейные восполнения сеточных функций.
- •Понятие о методе конечных элементов.
- •10. Общее описание метода.
- •20. Пример построения и исследования конечно-элементой схемы.
- •Свойства приближенного решения мнк.
- •10. Вводные сведения.
- •20. Вариационные свойства приближенного решения.
- •30. Формулы вложения.
- •Сходимость приближенного решения к точному решению.
- •1.9. Заключение по модулю
- •Вопросы для самоконтроля
- •1.10. Проектное задание
- •Контрольная работа №1.
- •Контрольная работа №2.
- •1.11. Тест рубежного контроля
- •Бланк правильных ответов
Введение в метод конечных элементов
Метод конечных элементов носит название вариационно-разностного метода. Реже встречается его другое название как проекционно-сеточный метод.
Этот метод применяют в случае сложной геометрии области, где ищется решение задачи, также является наиболее естественным, когда исходная задача имеет вариационную формулировку.
§1
Кусочно-линейные восполнения сеточных функций.
Фактически с кусочно-полиномеальными восполнениями мы имели дело при рассмотрении интерполирования, когда разговор шел о сплайнах.
Поставим задачу.
Пусть - непрерывная функция на отрезке . Построим, в общем случае, не равномерную сетку . Пусть в узлах сетки заданы значения функции .
Необходимо построить кусочно-линейную, т.е. линейную на каждом отрезке, функцию , которая в узлах сетки совпадает с исходной функцией , т.е. .
Рассмотрим два способа построения неизвестной функции.
1 -й способ построения функции .
.
(1) ,
,
.
Функция не имеет производных в точках , но в точках имеются односторонние производные и можно разумным способом доопределить полную производную в этих точках.
2-й способ построения функции .
Этот способ связан с введением базисных функций.
П остроим базис.
(2)
(3) .
___________________________________________________________________________________________________
‼ Доказать, что построенная система функций линейно не зависима, т.е. любая линейная комбинация равна нулю только тогда когда все коэффициенты равны нулю на отрезке [a,b]. Показать, что все функции базиса квадратично интегрируемы, т.е. .
___________________________________________________________________________________________________
Построенная система функций является базисом в классе кусочно-линейных функций, определенных на отрезке [a,b] у которых углы наклона могут меняться только в узлах сетки.
Очевидно, что
(4) .
Особо выделим случай финитных функций: , тогда:
(5) .
В дальнейшем нам пригодятся производные, вычислим их.
(6)
Удобство формул (2), задающих базисные функции состоит в том, что носитель функции составляет существенно малую часть от области определения функции.
Для дальнейшего нам потребуются леммы, которые позволяют оценить среднеквадратичное уклонение функции от в зависимости от гладкости функции и шагов сетки.
Лемма 1
Пусть функция , т.е. . Тогда имеет место неравенство: (7) .
Доказательство:
Рассмотрим на интервале выражение вида .
.
Итак, получено равенство:
(8) .
Возведем равенство (8) в квадрат и проинтегрируем.
(9) .
Оценим следующее выражение:
применим неравенство Коши – Буняковского:
.
Итак, была получена оценка
.
Вернемся к равенству (9).
.
Из последнего неравенства получаем,
.
Взяв , затем суммируя полученное равенства по , мы получаем неравенство (7).
Лемма 2 Оценка среднеквадратичного отклонения функции от .
Пусть функция , т.е. . Тогда имеет место оценка:
(10) .
Доказательство:
Имеет место тождество
откуда следует равенство (10).
___________________________________________________________________________________________________
‼ Получить аналог оценок (7) и (10) в равномерной метрике .
___________________________________________________________________________________________________
§2