- •Фгбоу впо «Воронежский государственный технический университет» г.Е. Шунин с.А. Кострюков в.В. Пешков
- •Воронеж 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
5 Моделирование сверхпроводникового гравиинерциального датчика
Целью вычислительного эксперимента по моделированию сверхпроводникового гравиинерциального датчика было определение его электромеханических характеристик, необходимых для построения уравнений движения ПТ и анализа его отклика на гравиинерциальные поля.
В режиме постоянного тока обобщенная сила, действующая вдоль координаты qi , определяется выражением [42]
, i = 1, ..., 6, (5.1)
где L – индуктивность системы, координатам q1 , q2 , q3 соответствуют линейные степени свободы пробного тела, а координатам q4 , q5 , q6 – угловые степени свободы. При i = 1, 2, 3 в (5.1) понимаются силы Fq, действующие на пробное тело по координатам q1 = x, q2 = y, q3 = z, а при i = 4, 5, 6 – моменты сил Мq при поворотах пробного тела на углы q4 = , q5 = , q6 = вокруг осей q1 , q2 , q3 соответственно.
Матрица жесткостей подвеса, характеризующая отклик системы на возмущения по всем степеням свободы, определяется как
, i, j =1, ..., 6, (5.2)
а собственные частоты колебаний соответствующих степеней свободы –
, i = 1, 2, 3;
, i = 4, 5, 6. (5.3)
Таким образом, было необходимо получить ряд зависимостей: L=L(z), L=L(y), L=L(), Fz=Fz(z), Fz=Fz(y), Fz=Fz(), Fy=Fy(), Mx=Mx(), Mx=Mx(y). Для определения зависимостей L=L(z), Fz=Fz(z) достаточно решить серию осесимметричных задач, а для остальных, когда расчетная область не обладает осевой симметрией (угловые или боковые смещения ПТ), задачи становятся трехмерными. Для уменьшения числа степеней свободы и облегчения процесса разбиения витки плоской и цилиндрической катушек подвеса были заменены кольцами прямоугольного сечения, плотно прилегающими друг к другу. При этом линия разреза задавалась одна со скачком, равным суммарному току витков. Такая аппроксимация соответствует случаю пренебрежения неоднородностью поля в небольшой окрестности витков. В диссертации [19] показано, что при этом интегральные характеристики подвеса меняются не слишком сильно, причем L и Fz имеют значения на ~10 % меньше, чем в случае учета геометрии структуры витков. Это обеспечивает сокращение числа степеней свободы в осесимметричных задачах почти в 3 раза, а решить трехмерную задачу на ПК без такой замены вообще не представляется возможным.
5.1 Геометрическая модель датчика
Были рассмотрены два варианта конструкции цилиндрического подвеса, различающихся формой одной из катушек. В первом из них использовалась плоская однослойная катушка в виде диска, расположенная в непосредственной близости от пробного тела, во втором подъемная сила создавалась катушкой квадратного сечения, значительно удаленной от ПТ (рис. 5.1–5.3).
А) Цилиндрический подвес с плоской катушкой
Принципиальная схема подвеса данной конструкции представлена на рис. 5.1. Он состоит из сверхпроводникового пробного тела (ПТ) 1, выполненного в виде полого цилиндра с крышкой, плоской однослойной сверхпроводниковой катушки 2 и цилиндрической однослойной сверхпроводниковой катушки 3, намотанных на каркасе 4 из гиперпроводника или сверхпроводника. Гиперпроводящий каркас выполняет также роль электромагнитного демпфера колебаний ПТ вблизи положения равновесия. Катушки 2 и 3 соединены последовательно и образуют единую короткозамкнутую сверхпроводящую цепь с тепловым ключом и контактами для запитки током. При запитке их током I создается неоднородное магнитное поле, в котором, вследствие эффекта Мейсснера, осуществляется левитация сверхпроводникового ПТ в некотором равновесном положении, соответствующем рабочему зазору d. Требуемое значение зазора достигается за счет подбора величины тока запитки I. Подвес может работать в двух режимах: постоянного тока и постоянного магнитного потока в короткозамкнутой сверхпроводящей цепи катушек. Информация о смещении ПТ под действием гравиинерциальных сил считывается с помощью емкостного датчика 5, представляющего собой два электрода в виде кольца и диска с равными площадями поверхности.
Параметры подвеса: плоская катушка – 35 витков, радиус внешний и внутренний 7.5 и 3.15 мм соответственно; цилиндрическая катушка – 40 витков, радиус 9.05 мм; длина боковой стенки ПТ 10.6 мм, радиус 10 мм, радиус отверстия 2 мм; толщи-
Рис. 5.1. Конструктивная схема датчика (вариант а).
Р ис. 5.2. Геометрическая модель датчика: 1 – сверхпроводящий экран; 2 – диэлектрическая подложка; 3 – емкостной датчик смещений ПТ; 4 – сверхпроводящее пробное тело; 5-6 – сверхпроводящие катушки; 7 – каркас.
на стенок 0.4мм; диаметр проволоки 0.1мм. Масса ПТ m=3.85г, момент инерции J относительно осей x, y равен 2.1гcм2.
В качестве сверхпроводникового конструкционного материала используется ниобий, обладающий первым критическим полем Bс1 ~ 0.14 Тл при T~4K. Соответствующее значение магнитного давления Fm ~ 8103 Нм–2.