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


РЕАЛИЗАЦИЯ МКЭ ПРОИЗВОЛЬНОГО ПОРЯДКА ТОЧНОСТИ ДЛЯ РЕШЕНИЯ КРАЕВЫХ ЗАДАЧ ДЛЯ НЕЛИНЕЙНЫХ ОДУ ВТОРОГО ПОРЯДКА В MATLAB

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

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

ВВЕДЕНИЕ 3
1 Смешанная краевая задача 6
1.1 Интегральное тождество 6
1.2 Определение схемы МКЭ 7
1.3 Система уравнений МКЭ 8
1.4 Метод Ньютона 10
2 Алгоритм вычисления Fи F0 для схемы МКЭ 11
2.1 Алгоритм вычисления F 11
2.1.1 Алгоритм сборки вектора F (вклад элементов) 12
2.1.2 Алгоритм вычисления F' 15
2.2 Алгоритм вычисления F0 15
2.2.1 Алгоритм сборки матрицы K(вклад элементов) 17
2.2.2 Алгоритм вычисления K' 18
3 Ортогональные полиномы, квадратуры типа Гаусса и Лобатто 20
3.1 Ортогональные полиномы 20
3.2 Квадратурные формулы типа Гаусса и Лобатто 22
4 Результаты тестовых расчетов 25
4.1 Тестовая задача 1 25
4.2 Тестовая задача 2 28
4.3 Тестовая задача 3 31
ЗАКЛЮЧЕНИЕ 33
СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ 35
ПРИЛОЖЕНИЕ. Код программы


Дифференциальные уравнения представляют собой математические модели различных процессов и явлений. Они обширно применяются на практике, в частности с целью описания переходных процессов.
Появлению нелинейных дифференциальных уравнений способствовали нелинейные задачи в механике. В этих уравнениях участвовали координаты тел, их скорости и ускорения, которые рассматриваются как функции от времени. На данный момент нелинейные дифференциальные уравнения изучены намного меньше, чем линейные.
Нелинейные ОДУ - раздел математики, изучающий теорию и способы решения нелинейных уравнений, содержащих искомую функцию и ее производные различных порядков одного аргумента (обычные нелинейные дифференциальные) или нескольких аргументов (нелинейные дифференциальные уравнения в частных производных).
Данные уравнения возникают во многих физических задачах: в механике, термоупругости, оптике, гидродинамике, аэродинамике, радиотехнике, электротехнике и др., но они ограничиваются не только физикой, но и другими областями науки: биологии, химии, экономики. Редко удается получить точное решение данных уравнений, поэтому решения находятся приближенно. Для поиска приближенного решения используют различные численные методы, например: метод конечных разностей, метод Галеркина и метод конечных элементов.
Целью выпускной работы является программная реализация в среде Mat Lab метода конечных элементов произвольного порядка точности для решения на конечном отрезке (a, b) уравнения
--d k(x,u, ) + k0(x,u, d--) = 0
dx dx/ dx/
при произвольном сочетании краевых условий первого рода
u(a) = ga; u(b) = дь,
или третьего рода
, / . . du(a)
—к a,u(a), —-—
X dx )
ko(b,u(b),^^u(b^- + аь (u(b))
dx
в граничных точках отрезка [a, b]. Здесь функции k = k(x,u,z), k0= = k0(x,u, z), aa = aa(u) и ab = ab(u) и чпс.1 a ga, gbявляются заданными.
Для выполнения этого задания необходимо:
1) На примере смешанных краевых условий получить слабую формулировку задачи; дать определение пространства решений и тестовых функций МКЭ; дать определение схемы МКЭ с численным интегрированием на основе квадратурных формул Гаусса и Лобатто.
2) Получить алгоритм формирования системы МКЭ.
3) Получить алгоритм приближенного определения производной Гато оператора схемы МКЭ.
4) Определить способ представления исходных данных о решаемой задаче и схеме МКЭ и программно реализовать пункты 4,5.
5) Написать и отладить программу решения схемы МКЭ, пользуясь необходимыми функциями MatLab.
6) Написать и отладить функцию решения тестовых задач и определения порядка точности метода.
7) Написать функцию решения задачи, точное решение которой неизвестно, представляя решение и его производную в графическом виде.
Также на основе написанной программы необходимо:
а) Решить ряд тестовых задач с известным точным решением и определить на основе вычислений порядок точности найденного решения и его производной в зависимости от размера элемента и степени полинома в норме Lfa, b).
б) Решить ряд задач, точное решение которых неизвестно.
Результатом работы является комплекс программ в среде MatLab, направленный на решение краевых задач для нелинейных обыкновенных дифференциальных уравнений второго порядка вышеуказанным методом, изученным в книгах [1,2,3]. Все написанные функции были тщательно про-тестированы.

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

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

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


В ходе выполнения выпускной квалификационной работы был реализован отлаженный комплекс программ в среде MatLab метода конечных элементов произвольного порядка точности для решения на конечном отрезке [a, b] уравнения
2k(x,u,^+N,UP = 0 dx dx/ dx/
при произвольном сочетании краевых условий первого
u(a) = ga, u(b) = gb,
или третьего рода
k0^b, u(b)
в граничных точка отрезка [а, Ь]. Также были написаны и отлажены функции:
1) Функция решения тестовых задач с известным точным решением и определения порядка точности метода на основе вычислений в зависимости от размера элемента и степени полинома в норме L2(a, b).
2) Функция решения задач, точное решение которых неизвестно и представления полученного решения и его производную в графическом виде.
Во время проведения работы было изучено большое количество теоретического материала, научных статей. На основе решений на различных сетках ряда тестовых задач, приведенных выше, можно сказать что метод конечных элементов является простым, гибким и оптимальным методом и позволяет наиболее эффективно решать нелинейные краевые задачи. Все цели были успешно достигнуты, тесты показали точность формул и доказали, что написанные программы являются надежными и эффективными.


[1] Сьярле Ф. Метод конечных элементов для эллиптических задач М.: Мир, 1980. 512 с.
[2] Р.З. Даутов, М.М. Карчевский Введение в теорию метода конечных элементов. Изд. 2. исправленное. - Казань: Казанский университет. - 2011. - 238 с, 2012, 238 с..
[3] Р. 3. Даутов, Програмная реализация метода конечных элементов в MATLAB — Казань, Изд-во КФУ, 2014.
[4] Глазырина Л.Л., Карчевский М.М., Введение в численные методы, КФУ, 2012
[5] Бахвалов И.С., Жидков Н.П., Кобельков Г.М., Численные методы, М.: Наука, 1987.
[6] Сегё Г., Ортогональные многочлены, М.: Физматгиз, 1962.


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




© 2008-2022 Cервис помощи студентам в выполнении работ