- •Фгбоу впо «Воронежский государственный технический университет» г.Е. Шунин с.А. Кострюков в.В. Пешков
- •Воронеж 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. Составить систему линейных алгебраических уравнений, получающуюся при аппроксимации по Галеркину решения задачи 2/x2+2/y2=0; x,y[–1, 1], =100 при y=1; /n= =–1– при x=1. Ввести три невязки: в области и на двух границах, соответствующих различным краевым условиям.
2. Рассмотреть дискретизацию уравнения 2/x2=–(x) методом взвешенных невязок, если в качестве весовых функций использовать -функции Дирака: Wj=(x–xj). Каким образом можно обеспечить выполнение краевых условий Дирихле и Неймана?
3. Провести дискретизацию уравнения
2/x2+2/y2+2/y=0
проекционным методом. Функцию на Г считать заданной.
4. Показать, что если исключить из рассмотрения невязку на Г2, то решение, полученное с помощью метода взвешенных невязок, будет удовлетворять условию /n=0.
5. Получить систему алгебраических уравнений для задачи (2.8), при условии, что =().
2.4 Конечные элементы и аппроксимация
Конечным элементом внутри рассматриваемой области называют некоторую подобласть e, геометрические размеры которой очень малы по сравнению с размерами области , но при этом остаются конечными. Элемент характеризуется числом геометрических узлов, типом аппроксимирующих функций и степенью аппроксимации неизвестной функции. Границы элементов могут быть как прямолинейными, так и криволинейными.
Решение трехмерной краевой задачи внутри элемента e можно представить в виде суммы
. (2.25)
Базисные функции Nm(e), относящиеся к элементу e, называются функциями формы этого элемента. Вне пределов e они тождественно равны нулю. Семейство {Nm(e)} должно обладать свойством полноты, т.е. при M (e)–(e) 0. Каждая функция формы обычно связывается с узлом m, причем Nm(e)(xn, yn, zn)=mn , т.е. она равна нулю во всех узлах, за исключением m-го.
На всей области определения решение можно представить в виде
. (2.26)
Последняя запись символическая; значение функции в точке (x, y, z) всецело определяется параметрами m, связанными с элементом e, которому принадлежит эта точка. Если обратиться к вариационной формулировке задачи, легко показать, что полный функционал задачи, относящийся ко всей области определения, равен сумме функционалов, вычисленных на каждом элементе. Поэтому необходимое условие экстремума можно применять для каждого конечного элемента, получать систему типа (2.15), причем интегрирование в (2.16 – 2.17) производится только по e и границе e, если последняя является частью . Сформированные таким образом локальные системы уравнений добавляются в одну общую (глобальную) систему по определенным правилам, рассмотренным ниже.
Поскольку при реализации МКЭ требуется вычислять интегралы от функций формы и их производных по области конечного элемента e, целесообразно ограничиться рассмотрением стандартного элемента e*, имеющего правильную форму. Поэтому, чтобы сделать процесс вычислений единообразным для всех элементов, делается замена
.
Здесь det J – якобиан преобразования. Теперь функции формы {Nm(e)} определяются раз и навсегда для стандартного элемента e*.
Рассмотрим МКЭ-процедуру формирования матрицы системы уравнений и вектора правых частей для задачи (2.8) на примере треугольника Лагранжа (двумерный случай). Такой конечный элемент обеспечивает на границах элементов непрерывность только функции, но не ее производных.
П
Рис. 2.1.
усть имеется треугольный элемент с вершинами (xi, yi), i=1, 2, 3. Введем произвольную точку Р с координатами (x, y). Определим ее положение с помощью симплексных координат по правилуi = Ai/A, i = 1, 2, 3, (2.27)
где А – площадь треугольника, Аi – площадь треугольника, образованного точкой Р и двумя вершинами исходного треугольника. Очевидно, что
1 + 2 + 3 = 1 и 0 i 1 внутри треугольника
Поскольку
то
. (2.28)
Формулы (28) преобразуют треугольный элемент произвольной формы, размера и ориентации в стандартный элемент (0,1), (1,0), (0,0) (рис. 2.1). Положим для простоты в (8) = const = 1, q = 0, а для возьмем аппроксимацию:
(2.29)
В случае линейной аппроксимации
Ni(e) = i, i = 1, 2, 3, (2.30)
. (2.31)
В результате локальная система уравнений запишется в виде:
S111 + S122 + S133 = T111 + T122 + T133,
S211 + S222 + S233 = T211 + T222 + T233, (2.32)
S311 + S322 + S233 = T311 + T322 + T333,
где
Подобным образом формируется система для каждого конечного элемента. Рассмотрим процесс объединения таких локальных систем в одну глобальную.
Пусть имеется два элемента с общей стороной 2–3 (рис. 2.2). Проводя дискретизацию на каждом конечном элементе, получим две системы вида (32). Переходя к глобальной нумерации узлов сетки и учитывая, что для обеспечения непрерывности функции ее значения i в совпадающих узлах должны быть равны, т.е.
2(1) = 2(2) 2 и 3(1) = 3(2) 3 ,
левые части уравнений для первого элемента примут вид:
, (2.33а)
, (2.33b)
, (2.33с)
Д
Рис. 2.2.
ля второго элемента, (2.34а)
, (2.34b)
. (2.34c)
Правые части уравнений получаются из левых заменой SijTij, ii. Складывая попарно уравнения для общих узлов ((2.33b) и (2.34a), (2.33c) и (2.34b)), получаем систему уравнений МКЭ для ансамбля из двух элементов:
,
(2.35)
.
Приведенный способ объединения элементов обобщается на случай любого числа элементов. Нетрудно заметить, что полученная матрица будет симметричной, положительно определенной и сильно разреженной (т.е. содержащей много нулей вне главной диагонали). При определенной нумерации можно добиться ленточного вида матрицы, при котором все ненулевые элементы группируются возле главной диагонали. Такой способ нумерации может иметь важное значение для компактного хранения матрицы. Симметричность и положительная определенность выгодны в вычислительном плане.
Можно показать, что сформированная в виде (32) и (35) система уравнений не будет иметь решения, поскольку ее определитель равен нулю. Необходимо модифицировать систему, чтобы учесть главное краевое условие (условие Дирихле) на Г1. Пусть в узлах 1 и 2 функция принимает фиксированные значения 1 и 2 соответственно. Тогда первое и второе уравнения в (35), соответствующие этим узлам, заменяются уравнениями
1 1 = 1 ,
1 2 = 2 ,
а в других уравнениях члены, содержащие 1 и 2, переносятся в правую часть с обратным знаком. Например, для 4-го уравнения в (35) имеем
.