- •Фгбоу впо «Воронежский государственный технический университет» г.Е. Шунин с.А. Кострюков в.В. Пешков
- •Воронеж 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
Задания
11. Для изображенных областей и указанных граничных условий найти распределение электростатического поля, построить эквипотенциальные кривые и графики изменения напряженности поля вдоль некоторых линий. Для проводящих границ использовать различные варианты условий (u=u0, u=const или .
I. Задачи с плоской геометрией
1) Плоская пластина во внешнем поле
(отношение сторон 41, 81, 121); |
2) Две пластины во внешнем поле
a:b:c = 8:1:5; 6:1:2; 4:1:4 |
3) Бесконечно длинный экран с поперечным сечением вида:
a:b=2:1; 1:1; c:a=1:4; 1:3; 1:2; d:a = 1:12; 1:8; 1:4 |
4) Система электродов во внешнем поле
|
5) Цилиндрический конденсатор с электродами внутри (показана четверть области). r2/r1 = 2; 3; 4; = 15; 30; d = 0.1(r2–r1); 0.2(r2–r1); 0.4(r2–r1) |
|
6) Заряженные проводники во внешнем поле
a:b = 1:1; c:a = 1:5; d:c = 1:3; 1:5.
|
|
II. Осесимметричные задачи
7) Цилиндрический экран
|
8) Полый шар с отверстием
|
|
9) Две полусферы во внешнем поле
|
10)Кольцо П-образного сечения
|
|
10) Цилиндрический конденсатор с двумя кольцами внутри а) кольца проводящие (u=const) б) кольца из диэлектрика ( 1) |
|
12. Для односвязных сверхпроводниковых систем, представленных в заданиях 23-26 (стр. 34-36) с соответствующими конфигурациями, найти распределение магнитного поля через векторный магнитный потенциал A.
Указание. В случае плоской геометрии данная задача сводится к уравнению Лапласа A=0 с условием A=const на поверхности односвязного сверхпроводника и граничными условиями 1-го и 2-го рода для учета внешнего поля.
Лабораторная работа № 4 Решение краевой задачи для уравнения Пуассона
Уравнение Пуассона имеет вид
.
Ему подчиняется, например, стационарное распределение температуры или потенциала в области с распределенными источниками. Оно, как правило, дополняется граничными условиями 1-го, 2-го и 3-го рода. Могут ставиться и дополнительные условия по типу рассмотренных в лабораторной работе №3.
Пример 6. Найти стационарное распределение температуры в пластине L-образной формы, рассмотренной в примерах 2 и 3, все стороны которой поддерживаются при нулевой температуре. Пластина нагревается постоянным током, выделяющим в единице объема тепло Q. Данная задача сводится к решению уравнения Пуассона u=–Q/k (k – коэффициент теплопроводности) в соответствующей области с условием 1-го рода u|Г = 0.
Порядок решения:
1) Выполнить пункты 1–4 из примера 3.
2) Удалить граничные условия на всех линиях («Удал» «Все ГУ»).
3) Ввести на всех граничных линиях условие Дирихле u=0.
4) Задать уравнение u=–1 («Файл» «Уравнение», отметить пробелом «Пуассона», перейти на «More», определить функцию F(u,r,t) как const0 и ввести число 1, нажать Enter и Esc).
5) Осуществить проверку введенных параметров («Файл» «Информация» …)
6) Выйти из программы pre2d. Дальнейшие действия – решение и анализ результатов – как в рассмотренных выше примерах.
Пример 7. Найти распределение температуры в прямоугольной пластине, внутри которой имеется источник тепла круглой формы с постоянной плотностью (расчетная область показана на рис. П13). На сторонах пластины поддерживаются условия обмена тепла с внешней средой по закону .
Задача сводится к решению уравнения Пуассона только внутри подобласти с источником (т.е. круга), в остальной части области решается уравнение Лапласа. Из-за симметрии задачи ограничимся рассмотрением только половины области.
Порядок решения:
1) В препроцессоре загрузить геометрию задачи, решенной в примере 4 (вместе с сеткой).
2) Отменить все граничные и дополнительные условия, выбирая соответствующий пункт в меню «Удал».
3) Разбить зоны с номерами 1, 5, 9, 12, 16, 20, если они еще не разбиты, на конечные элементы.
4) Поскольку согласно параметрам загруженной задачи для всех зон установлено уравнение Лапласа, необходимо заменить его на уравнение Пуассона, даже если последнее будет действовать лишь для некоторых зон («Файл» «Уравнение», отметить пробелом «Пуассона», перейти на «More», нажать Enter и Esc).
Для указанных в п. 3 зон ввести отличную от нуля плотность источника тепла, например, Q=3 («Файл» «Неоднородности среды» [1, 5, 9, 12, 16, 20] <Enter> «Вид уравнения» «F(u,r,t)» «const0» [3] <Enter> <Esc> <Esc>).
5) Задать граничные условия на левой, правой и верхней сторонах прямоугольной области: («Файл» «Граничные условия» «Смешанные», затем по запросу ввести номера линий, на которых задаются условия: 38, 42, 54, 53, 51 и т.д., далее, выбрав для функции B(r, t) пункт «const0», ввести число 1, после перейти с помощью Tab на функцию «Q(r,t)», определив для нее «0», и выйти из режима ввода граничных условий, с помощью клавиш Tab и Enter).
На нижней стороне прямоугольника выполняется однородное условие Неймана (его специально задавать не нужно).
6) Проверить все введенные данные («Файл» «Информация») и выйти из препроцессора.
7) Провести расчет распределения функции u (программы appl_fem, difeqt) и войти в постпроцессор (post2d), в котором вывести на экран картинку распределения поля (линии уровня, матрица стрелок, цветовая карта), построить графики функции u, ее градиента вдоль линии y=0, проверить выполнение условий на границах области.
Модификации. Провести аналогичный расчет поля для других граничных условий:
а) u|x=a=0, u|x=b=0, , ;
б) , , , ;
в) , , ,
(здесь предполагается axb, 0yd).
Решить данную задачу как осесимметричную, считая, что геометрия области задана в координатах (r, z) цилиндрической системы, принимая z=x, r=y. В этом случае подобласть источников поля будет представлять собой шар, а вся область задачи – цилиндр.