- •Фгбоу впо «Воронежский государственный технический университет» г.Е. Шунин с.А. Кострюков в.В. Пешков
- •Воронеж 2014
- • Шунин г.Е., Кострюков с.А., Пешков в.В., 2014
- •1 Современное состояние разработок сверхпроводящих электромагнитных подвесов
- •1.1 Основные типы сверхпроводящих электромагнитных подвесов
- •1.2 Методы расчёта и компьютерного моделирования сверхпроводящих подвесов
- •1.3 Компьютерные системы конечно-элементного анализа
- •2 Основные положения метода
- •2.1 Сущность метода конечных элементов
- •2.2 Вариационные методы дискретизации
- •Упражнения
- •2.3 Проекционные методы дискретизации
- •Упражнения
- •2.4 Конечные элементы и аппроксимация
- •Упражнения
- •2.5 Решение обыкновенных дифференциальных уравнений
- •Упражнения
- •2.6 Решение дифференциальных уравнений с частными производными
- •Упражнения
- •3 Физико-математическое моделирование конструктивных элементов сверхпроводящих электромагнитных подвесов
- •3.1 Физико-математическая модель
- •3.2 Конечно-элементная дискретизация уравнений
- •Упражнения
- •3.3 Особенности решения задач для открытых многосвязных систем
- •Упражнение
- •3.4 Моделирование экранов
- •Упражнения
- •Упражнения
- •Упражнение
- •4 Конечно-элементный комплекс программ fempdesolver
- •4.1 Структура и возможности комплекса программ fempdeSolver
- •4.2 Препроцессор
- •4.3 Процессор
- •4.4 Постпроцессор
- •5 Моделирование сверхпроводникового гравиинерциального датчика
- •5.1 Геометрическая модель датчика
- •А) Цилиндрический подвес с плоской катушкой
- •Б) Цилиндрический подвес с катушкой квадратного сечения
- •5.2 Моделирование распределения магнитного поля в рабочем объеме датчика
- •А) Цилиндрический подвес с плоской катушкой
- •Б) Цилиндрический подвес с катушкой квадратного сечения
- •5.3 Моделирование распределения электростатического поля в емкостном датчике смещений пробного тела
- •Лабораторная работа № 1 Решение краевой задачи для обыкновенного дифференциального уравнения
- •Задания
- •Лабораторная работа № 2 Решение краевой задачи для уравнения Лапласа
- •Задания
- •Лабораторная работа № 3 Решение краевой задачи для уравнения Лапласа с дополнительными условиями
- •Задания
- •I. Задачи с плоской геометрией
- •II. Осесимметричные задачи
- •Лабораторная работа № 4 Решение краевой задачи для уравнения Пуассона
- •Задания
- •Лабораторная работа № 5 Решение краевой задачи при наличии физически неоднородных сред
- •Задания
- •Лабораторная работа № 6 Решение уравнения Лапласа в области с разрезами
- •Задания
- •Лабораторная работа № 7 Решение краевой задачи для уравнения Лондонов
- •Задания
- •Лабораторная работа № 8 Решение смешанной задачи для уравнения теплопроводности
- •Задания
- •Лабораторная работа № 9 Сверхпроводниковые подвесы
- •Описание интерфейса препроцессора
- •«Выход»
- •Горячие клавиши препроцессора
- •Программа appl_fem
- •Процессор
- •Описание интерфейса постпроцессора
- •Меню «Файл»
- •Меню «Вид»
- •Меню «Поле»
- •Меню «График»
- •Меню «Таблица»
- •Меню «Печать»
- •Меню «Вычислить»
- •Меню «Опции»
- •Меню «Помощь»
- •«Выход»
- •Учебное издание
- •394026 Воронеж, Московский просп., 14
Упражнения
1. Найти функции формы Ni = i (i=1, 2, 3, 4) для тетраэдрального элемента Лагранжа 1-го порядка, используя преобразование
x = Xii, y = Yii, z = Zii, 1 = i.
Показать, что внутри тетраэдра 0 1.
2. Найти квадратичные функции формы для одномерного элемента Лагранжа [x1,x2]. Указание. Использовать аппроксимацию неизвестной функции =C12+ C2+C3 , где =(x–x1)/(x2–x1) – симплексная координата, {Сi} – коэффициенты, подлежащие определению. Записав это выражение для каждого из трех узлов (=0; 0.5; 1), получим систему уравнений относительно {Сi}. Решив ее, найдем разложение = i Ni().
3. Определить квадратичные функции формы для стандартного треугольного элемента Лагранжа (см. рис.1). См. указание к упр. 2; рассмотреть аппроксимацию =C112+ C222+ C312+ C41+ C52+С6 для точек (0,1), (0,0), (1,0), (0,0.5), (0.5,0), (0.5,0.5).
4. Определить функции формы для билинейного прямоугольного элемента –11, –11. См. указание к упр. 2; рассмотреть аппроксимацию =C1 + C2 + C3 + C4 для точек (–1,–1), (–1,1), (1,–1), (1,1).
5. Показать, что функции формы должны удовлетворять условию Ni = 1.
6. Записать выражение для grad для линейного и квадратичного треугольника Лагранжа. Найти скачок градиента на границах элементов.
7. Показать, что полный функционал для ансамбля конечных элементов не равен сумме элементных функционалов, если для аппроксимации неизвестной функции используются базисные функции класса С0. Какую ошибку при этом имеет конечно-элементное решение?
8. Для уравнения Лапласа =0 получить матрицу и вектор правых частей локальной системы линейных алгебраических уравнений, соответствующей квадратичному одномерному элементу, билинейному четырехугольнику и тетраэдральному элементу Лагранжа.
2.5 Решение обыкновенных дифференциальных уравнений
Теоретические сведения. Обыкновенными дифференциальными уравнениями называются такие уравнения, которые содержат одну или несколько производных от искомой функции = (x). Их можно записать в виде
,
где x – независимая переменная.
Общее решение обыкновенного дифференциального уравнения n-го порядка содержит n произвольных постоянных C1, C2, ... , Cn, т. е. имеет вид
= (x, C1, C2, ... , Cn).
Частное решение обыкновенного дифференциального уравнения получается из общего, если произвольным постоянным придать определенные значения. В зависимости от способа задания дополнительных условий для получения частного решения дифференциального уравнения существуют два различных типа задач: задача Коши и краевая задача. Если эти условия задаются в одной точке, то такая задача называется задачей Коши. Дополнительные условия в задаче Коши называются начальными условиями, а точка x = x0, в которой они задаются, – начальной точкой. Если же дополнительные условия задаются в более чем одной точке, т.е. при разных значениях независимой переменной, то такая задача называется краевой. На практике обычно граничные условия задаются в двух точках x = a и x = b, являющихся границами области решения дифференциального уравнения.
МКЭ обычно используется для решения краевых задач. Задачи Коши чаще всего решаются методами Эйлера, Рунге–Кутта, Адамса, прогноза-коррекции и т. д.
Пример. Рассмотреть конечно-элементную формулировку для уравнения
при условии, что = 0 при x = 0 и /x = 0 при x = Lx.
Покажем, что стационарное значение функционала
определяет решение данной краевой задачи, если пробные функции удовлетворяют краевому условию = 0 при x = 0.
Найдем первую вариацию F:
.
Для первого слагаемого используем интегрирование по частям
.
Поскольку множество допустимых функций удовлетворяет краевому условию при x = 0, то x=0 = +x=0 = 0, первое слагаемое равно нулю.
Из условия стационарности F = 0 и в силу произвольности имеем:
.
Краевое условие при x=Lx выполняет роль естественного краевого условия, так как оно выполняется для любой функции , обеспечивающей стационарное значение функционала. Краевое условие при x=0 – главное условие.
Рис. 2.3.
Следуя методу конечных элементов, разобьем отрезок [0,Lx] на M одномерных элементов (рис. 3). На каждом элементе будем использовать линейную аппроксимацию. В таком элементе неизвестная функция определяется ее значениями в двух узлах x1 и x2 :
(x) = 1N1 + 2N2 , (2.36)
где функции формы N1, N2 имеют вид
.
Подставляя аппроксимацию (2.36) в функционал F() и удовлетворяя условию стационарности F/i = 0, приходим к системе линейных алгебраических уравнений
S111 + S122 = F1,
S211 + S222 = F2,
где
.
Здесь he x2 – x1.
Аналогичным образом строятся системы уравнений для всех элементов, а затем добавляются в глобальную систему способом, описанным в разд. 3. Чтобы окончательно сформировать систему, необходимо учесть главное граничное условие при x = 0 (см. там же).