Заказать работу


Тип работы:
Предмет:
Язык работы:


МЕТОД СОПРЯЖЕННЫХ ГРАДИЕНТОВ С ПРЕДОБУСЛАВЛИВАНИЕМ В СХЕМАХ МКЭ

Работа №42073
Тип работыДипломные работы, ВКР
Предметматематика
Объем работы70
Год сдачи2019
Стоимость4900 руб.
ПУБЛИКУЕТСЯ ВПЕРВЫЕ
Просмотрено 39
Не подходит работа?

Узнай цену на написание

ВВЕДЕНИЕ 3
1 Метод конечных элементов 5
1.1 Схема метода конечных элементов для двумерной краевой задачи , , , , 5
2 Метод сопряженных градиентов для решения СЛАУ 7
3 Предобуелавливание для метода сопряженных градиентов 8
3.1 Метод Якоби 9
3.2 Метод SSOE 9
3.3 Метод декомпозиции области 9
3.4 Блочный метод Якоби 9
3.5 Блочный метод SSOR 9
4 Сравнение методов 11
4.1 Схема МКЭ для одномерной задачи е линейными конечными элементами 11
4.2 Схема МКЭ для одномерной задачи е квадратичными конечными элементами 14
4.3 Схема МКЭ для одномерной задачи е кубическими конечными элементами 18
4.4 Схема МКЭ для одномерной задачи е конечными элементами четвёртой
степени 21
4.5 Схема МКЭ для двумерной задачи 24
4.6 Анализ результатов 27
ЗАКЛЮЧЕНИЕ 29
СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ 30
ПРИЛОЖЕНИЕ 31

Многие явления можно описать с помощью станционарных процессов. Однако решение стационарных задач зачастую не может быть найдено аналитически, В этих случаях применяют различные численные методы, сводя задачу к конечномерной с помощью некоторого метода дискретизации. Одним из таких методов дискретизации является метод конечных элементов.
Однако при требовании высокой точности к решению исходной задачи размер итоговой системы линейных уравнений может неограниченно расти. Это приводит к проблеме решения больших систем линейных уравнений.
Методы решения систем линейных уравнений делятся на две основных группы: прямые и итерационные. Последние представляют из себя крайне широкую группу и выбор наиболее оптимального из этой группы является важной и в то же время не тривиальной задачей. Важно отметить, что у одних и тех же методов может быть несколько вариантов реализации, что ещё больше затрудняет выбор наиболее оптимального метода.
Данная работа посвящена реализациям нескольких вариантов предобуславливания для метода сопряженных градиентов в МКЭ для граничных задач,
С одной стороны, работа не претендует на универсальность - полученные в ней сравнительные эффективности разных методов подойдут не для любой СЛАУ, С другой стороны, она является не такой специализированной, чтобы рассматривать только один вид системы и потому полученные оценки можно считать применимыми к классу систем линейных уравнений, В частности, в данной работе рассматриваются системы, построенные по схеме МКЭ для одномерных и двумерных краевых задач. Основным свойством таких систем является сильная разреженность.
Кроме того, в этой работе производится сравнение итерационных методов не только между собой, но и с прямым методом. За счёт этого можно рассмотреть актуальность использования итерационных методов. Также в данной работе производится сравнение различных реализаций алгоритмов, одни из которых учитывают разреженную структуру матриц, а другие нет,
В рамках данной работы выполнены следующие задачи:
• Реализованы алгоритмы сборки системы метода конечных элементов с численным интегрированием для одномерной и двумерной краевых задач,
• Реализован метод сопряженных градиентов
• Реализованы различные методы предобуславливания для метода сопряженных градиентов
• Проведено сравнение различных видов предобуславливания для метода сопряженных градиентов между собой и сравнение метода сопряженных градиентов с прямым методом решения СЛАУ,Все реализованные алгоритмы были написаны на языке программирования python с использованием пакетов numpy и scipy.
В перспективе, хотелось бы расширить класс итерационных методов, участвующих в сравнении и исследовать более подробно возможности оптимизации методов, в частности распараллеливание вычислений, где возможно,
В рамках работы рассматривались две краевые задачи. Первая - одномерная краевая задача Дирихле вида:— ('p(x)U (x))' = f (x), Vx Е (a,b)
u(a) = Ua,
u(b) = ub.
Вторая - аналогичная двумерная задача Дирихле в квадратной области вида:
— div(p(x)Vu(x)) = f (х), Vx Е (0,1) х (0,1) u(x) = 0, x Е Г, где p(x) - скалярная функция векторного аргумента, Г - граница области (0,1) х (0,1),



Возникли сложности?

Нужна помощь преподавателя?

Помощь в написании студенческих
и аспирантских работ!


Из рассмотренных методов предобуславливания метод декомпозиции области проявил себя лучше всего. Данный метод оказался наиболее быстрым для большого набора тестов. Также было получено, что некоторые из итерационных методов сравнимы по времени решения задачи с прямым методом, В частности, к ним относятся блочные методы предобуелавливания.
При проведении работы для класса блочных предобуславливателей было выявлено, что учёт сильной разреженности системы положительно сказывается на времени решения задачи и может дать значительный прирост произодительноети.
Стоит отметить, что хоть ни один из методов не смог превозойти прямой метод по времени решения задачи, у некоторых методов есть ещё потенциал для увеличения производительности, В частности, метод декомпозиции может быть ускорен за счёт распараллеливания вычислений для независимых блоков. Кроме того, у итерационных методов есть важная особенность, которая в данной работе не была исследована - для итерационных методов может быть задано начальное приближение, что может дополнительно улучшить такой показатель, как время работы программы, К тому же исследованный список методов предобуславливания не является исчерпывающим. Все эти темы сподвигают на дальнейшее исследование использования метода сопряженных градиентов с предобуславливанием для систем МКЭ.


[1] Р, 3, Даутов, Программная реализация метода конечных элементов в MATLAB. Изд-во Казан, ун-та, 2014,
[2] В, Б, Андреев, Лекции по методу конечных элементов, Маке ПРЕСС, 2010,
[3] Д, Ортега, Введение в параллельные алгоритмы и векторные методы решения линейных систем,. Мир, 1991,
[4] А, А, Самарский и А, В, Гулин, Численные методы. Наука, 1989,


Работу высылаем на протяжении 30 минут после оплаты.

Пожалуйста, укажите откуда вы узнали о сайте!


Подобные работы


© 2008-2020 Cервис продажи образцов готовых курсовых работ, дипломных проектов, рефератов, контрольных и прочих студенческих работ.
.