- •Математический факультет Кафедра вычислительной математики и программирования
- •Содержание
- •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 Метод сеток для уравнения гиперболического типа
1.3.Варианты задания
Граничные условия
Вари- анты |
u/AB |
u/BC |
u/CD |
uAD |
1 2 3 4 5 6 7 8 9 10 |
30y 30y 50y(1-y2) 20y 0 30 sin( πy) 30(1- y) 50 sin( πy) 40y2 50y
|
30(1-x2) 30cos(πx/2) 0 20 50x(1-x) 20x 20√x 30√x 40 50(1-x)
|
0 30cos(πy/2) 0 20y2 50y(1-y2) 20y 20y 30y2 40 0
|
0 0 50 sin( πx) 50x(1-x) 50x(1-x) 30x(1-x) 30x(1-x) 50 sin( πx) 40 sin( πx/2) 60, 0 ≤ x < 1/2 60(1-x), 1/2≤ x < 1
|
Лабораторная работа №2
Используя метод сеток, составить решение дифференциального уравнения Лапласа
, удовлетворяющее на окружности заданным начальными условиями; шаг дискретизации h=1. Уточнение решения производить до сотых долей при помощи процесса Либмана.
Пример решения задачи.
U (x,y) |Г = 0,5 (|x| + |y|).
Используя симметрию заданных начальных условий, построим решение только в 1 четверти (рис.1). Возьмем шаг h=1 и составим таблицу значений x и y:
-
Х
0
1
2
3
4
У
3
2,90
2,60
1,98
0
Y
2
U4
U8
U9
Рис. 1
На рисунке A, B, C, D, E, F – граничные узлы, Ui , i = 1,2,3,4,5,6,7,8,9,10,11 – внутренние.
Вычислим значения функции U (x,y) на границе:
A(0;3); U(A) = 0,5*(0+3) = 1,5,
B(1;2,90); U(B) = 0,5*(1+2,9) = 1,95,
C(2;2,60); U(C) = 0,5*(2+2,6) = 2,3,
D(3;1,98); U(D) = 0,5*(3+1,98) = 2,49,
E(3,77;1); U(E) = 0,5*(3,77+1) = 2,39,
F(4;0); U(F) = 0,5*(4+0) = 2.
Для определения начальных значений функции U(x,y) во внутренних точках составим систему уравнений, содержащих эти значения. Каждое уравнение получается приравниванием значения функции во внутренней точке среднему арифметическому четырех значений функции в соседних точках:
U1 = (1,5 + U4 +2U2)/4, U2 = (1,95 + U1 +U3 +U5)/4, U3 = (4,79 + U2 + U6)/4,
U4 = (U1 + U8 + 2U5)/4, U5 = (U2 +U4 +U6 +U9)/4, U6 = (U3 +U5 + U7 + U10)/4,
U7 = (4,88 + U6 + U11)/4, U8 = 4U5 /4, U9 = (U8 + U10 + 2U5)/4,
U10 = (U9 +U11 + 2U6)/4, U11 = (4,78+2U6)/4.
Решая эту систему, получим U1 = 1,91, U2 = 2,05, U3 = 2,10, U4 = 2,05, U5 = 2,
U6 = 2,18, U7 = 2,34, U8 = 2,11, U9 = 2,13, U10 = 2,19, U11 = 2,28.
Найденные значения функции U(x,y) позволяют составить шаблон №1, в котором внутренние значения соответствуют найденным, а граничные получаются в результате уточнения предыдущих граничных значений по формуле линейной интерполяции
Где Ah – узловая граничная точка, А – ближайшая к Ah точка, лежащая на границе;
Bh – ближайшая к Ah узловая точка, лежащая внутри области; - расстояние между точками А и Ah , взятое со знаком плюс, если точка Ah лежит внутри области, и со знаком минус, если она лежит вне области.
В данной задаче имеем:
№ 1.
1,5 |
1,94 |
2,43 |
|
|
1,91 |
2,05 |
2,10 |
2,49 |
|
2,05 |
2,11 |
2,18 |
2,34 |
2,40 |
2,11 |
2,13 |
2,19 |
2,28 |
2 |
Процесс Либмана заключается в уточнении значений, входящих в шаблон №1. Каждый следующий шаблон получается следующим образом из предыдущего: значения функции во внутренних точках равны среднему арифметическому четырех соседних значений предыдущего шаблона, а значения функции в граничных точках находятся по формуле линейной интерполяции, уже использованной при получении шаблона № 1. Это уточнение производится до тех пор, пока два последовательных шаблона не совпадут с заданной степенью точности. В результате вычислений получим следующую последовательность шаблонов:
№ 2. № 3.
1,5 |
1,94 |
2,31 |
|
|
|
1,5 |
1,94 |
2,33 |
|
|
1,91 |
2,02 |
2,29 |
2,49 |
|
|
1,90 |
2,06 |
2,25 |
2,49 |
|
2,06 |
2,10 |
2,18 |
2,34 |
2,40 |
|
2,05 |
2,10 |
2,23 |
2,32 |
2,41 |
2,09 |
2,13 |
2,19 |
2,22 |
2 |
|
2,10 |
2,12 |
2,18 |
2,22 |
2 |
№ 4 № 5.
1,5 |
1,94 |
2,31 |
|
|
|
1,5 |
1,94 |
2,33 |
|
|
1,92 |
2,05 |
2,28 |
2,49 |
|
|
1,91 |
2,06 |
2,26 |
2,49 |
|
2,05 |
2,12 |
2,21 |
2,34 |
2,40 |
|
2,06 |
2,11 |
2,23 |
2,33 |
2,41 |
2,09 |
2,12 |
2,20 |
2,20 |
2 |
|
2,09 |
2,14 |
2,9 |
2,22 |
2 |
№ 6. № 7.
1,5 |
1,94 |
2,31 |
|
|
|
1,5 |
1,94 |
2,32 |
|
|
1,92 |
2,06 |
2,28 |
2,49 |
|
|
1,92 |
2,06 |
2,27 |
2,49 |
|
2,062 |
2,12 |
2,22 |
2,34 |
2,40 |
|
2,06 |
2,12 |
2,23 |
2,33 |
2,41 |
2,20 |
2,13 |
2,20 |
2,22 |
2 |
|
2,10 |
2,20 |
2,22 |
2,22 |
2 |
№ 8
1,5 |
1,94 |
2,32 |
|
|
1,92 |
2,06 |
2,27 |
2,49 |
|
2,06 |
2,12 |
2,23 |
2,33 |
2,41 |
2,10 |
2,13 |
2,20 |
2,22 |
2 |
Шаблон №8 является ответом.
1.3. Задача. Используя метод сеток, составить решение дифференциального уравнения Лапласа
с заданным начальными условиями; шаг дискретизации h=1. Уточнение решения производить до сотых долей при помощи процесса Либмана.