Если одно из понятий незнакомо, вернитесь к предыдущему занятию и воспроизведите его численный пример. После этого переходите к новой постановке.
Место занятия в инженерной методике
Расчётная лабораторная работа · Выбор расчётной области, симметрии и граничных условий
Инженерный вопрос. Какая величина задана разностью A, а какая определяется решением смешанной краевой задачи?
Магнитная цепь, потокосцепление, тяговая характеристика и потери образуют методику выбора зазора, обмотки и режима привода. Сила становится проектным критерием после сопоставления с нагрузкой механизма.
Последовательность решения, отчёт и критерии результата
Результат получают из наблюдаемых величин и проверяют в условиях данного практикума.
Последовательность решения
Обосновать постановку и способ измерения. Проверьте совместимость границы с геометрией, источниками и искомым режимом. При секторе заранее выведите восстановление полного результата.
Определите проводимость собственной области. Выполните действия практикума по его исходной таблице. Сведите разность потенциала, глубину, поток, магнитное напряжение двух допустимых путей и проводимость сеточной серии.
Результат этого занятия. Таблица нормировок и сходимости с объяснением идеальных полюсных границ и отличия этой задачи от полного аппарата.
Данные именно этого опыта. Сведите разность потенциала, глубину, поток, магнитное напряжение двух допустимых путей и проводимость сеточной серии.
Оформление наблюдений. Размер области, тип границы, размер сектора и множитель; внутренняя контрольная величина и полный восстановленный интеграл.
Как обосновать результат
Проверка из практикума. В \Phi нет лишнего \mu _{0}; знаки и единицы объявлены; \lambda относится к точной геометрии и подтверждён сходимостью.
Независимый контроль. Сравните с большей областью либо полной геометрией, предусмотренной занятием; внешнюю границу меняйте отдельно от сетки.
Критерий выполнения. В \Phi нет лишнего \mu _{0}; знаки и единицы объявлены; \lambda относится к точной геометрии и подтверждён сходимостью.
Точность и расхождения. Оцените изменение одной контрольной величины при удалении границы или восстановлении полной модели. При близком к нулю результате укажите абсолютную разность и физический масштаб.
Основание сравнения. Используйте граничные условия и симметрию из теории занятия. Удалённая граница — приближение; симметрия допустима лишь при согласовании геометрии, источников и искомого режима.
Вывод для устройства
Обоснуйте сокращение вычислений без потери именно той характеристики, ради которой рассчитывается устройство.
Границы вывода. Магнитная сила отдельного привода не подтверждает работу всего аппарата: контактные усилия, пружины, трение, динамика и нагрев должны быть заданы и проверены по назначению.
Теория и вывод формул
Воздушная краевая задача позволяет определить магнитную проводимость без явного моделирования катушки. Два интегральных способа расчёта используются для сопоставления результатов и отделения масштаба возбуждения от геометрических и сеточных факторов.
Соглашения об обозначениях и единицах
\(i(t)\), \(u(t)\) — мгновенные значения; \(\widehat I\), \(\widehat U\) — комплексные амплитуды; \(\underline I_{\text{д}}\), \(\underline U_{\text{д}}\) — комплексные действующие значения. Для синусоиды \(\widehat I=\sqrt{2}\,\underline I_{\text{д}}\); действующее значение тока равно \(|\underline I_{\text{д}}|\). Нормировка обязательно учитывается в формулах мощности.
\(\Phi=\int_S\mathbf B\cdot\mathrm d\mathbf S\) — поток через ориентированную поверхность; \(\Psi=\sum_k\oint_{C_k}\mathbf A\cdot\mathrm d\mathbf l\) — потокосцепление последовательных витков с согласованными замкнутыми контурами. Только при одинаковом сцепленном потоке каждого витка \(\Psi=N\Phi\). Обе величины измеряются в веберах: число витков безразмерно. Запись «Вб·виток» подчёркивает суммирование.
\(\theta_{\text{м}}\), \(\Omega\) (также \(\omega_{\text{мех}}\)) — механические угол и угловая скорость. Для основной гармоники \(\theta_{\text{э}}=p\theta_{\text{м}}\), \(\omega_{\text{э}}=p\Omega\), где \(p\) — число пар полюсов. Геометрический период зубцов задаётся отдельно. В производных по углу используйте радианы.
Статическая индуктивность \(L_{\text{ст}}=\Psi/i\) — наклон секущей из начала координат; это определение используется при \(\Psi(0)=0\) и \(i\ne0\). Динамическая индуктивность \(L_{\text{дин}}=\partial\Psi/\partial i\) — наклон касательной (дифференциальная, инкрементальная); положение и остальные токи фиксированы. Для линейной характеристики без смещения обе равны одной индуктивности \(L\). При постоянных магнитах и взаимном возбуждении сначала учитывают ненулевое смещение потокосцепления.
Положительные направления, оси и специальные индексы определяются в конкретной задаче. Единицу берите из условия или заголовка таблицы; перед подстановкой в формулу СИ переводите миллиметры в метры, миллигенри в генри и градусы в радианы.
Численные результаты обычно показаны с четырьмя значащими цифрами; для сравнения близких значений в таблицах сходимости сохранено больше цифр. Полные значения остаются в исходных расчётных данных. Малые остатки вычисления координат на осях показаны как ноль; малые физические величины сохраняются в показательной записи.
П4 — воздушная краевая задача идеальных полюсов. Катушка явно не моделируется, сопротивление железа принято равным нулю. На поверхности идеальной стали B направлена нормально, поэтому для A_z задаётся естественное условие \partial A_z /\partial n=0. На известных линиях потока устанавливают A_z =\mathrm{const}. В воздухе J=0 и \Delta Az=0.
Самостоятельная рассчитанная геометрия курса: C(-50{,}0), D(0{,}50), E(0{,}25), M(40{,}25), N(40{,}0), мм. CD — четверть окружности радиуса 50 мм с центром (0,0). Область ограничена C–N–M–E–D–дуга–C. На CD задано A_z =0.002 Вб/м, на MN — A_z =0; на CN, DE и EM — нулевой нормальный градиент. Глубина равна 20 мм. Это конкретная смешанная краевая задача, а не аппроксимация открытой трёхмерной системы произвольной формы.
Поток из векторного потенциала
Для плоского поля B_x =\partial A_z /\partial y и B_y =-\partial A_z /\partial x. Если нормаль лежит слева от направления обхода, ориентированный интеграл Bn dl равен A_z в начале минус A_z в конце контура. При развороте контура знак потока меняется. Модуль полного потока между CD и MN равен b|\Delta Az|.
В формулах 3.33–3.34 книги на с. 134 напечатан лишний \mu _{0}. Это противоречит B=\operatorname{rot} A и числам самого П4. Здесь A_z имеет единицу Вб/м, поэтому \Phi =b\Delta Az, а \mu _{0} появляется при вычислении H=B/\mu _{0} и нормировке \lambda. При заданных входах модуль потока должен быть 40 мкВб.
Длины в формулах переводятся в метры; интеграл FEMM автоматически использует глубину, заданную в постановке.
Интеграл H и независимость от пути
В односвязной воздушной области без свободного тока интеграл H между двумя фиксированными точками не зависит от пути. Определите Um вдоль MN и по альтернативному воздушному пути с теми же концами. Не интегрируйте по несуществующему железу, которое заменено идеальным граничным условием.
Сравните интеграл Bn через полное сечение с b\Delta Az. Короткий отрезок только через середину торца может не охватывать выпучивание и даст другой поток. Проводимость относится к объявленным полюсным поверхностям и их магнитному напряжению, а не к локальному максимуму B.
Сходимость и масштаб
Размер внутренних элементов и максимальный угол сегмента дуги меняйте отдельно. Грубая дуга меняет геометрию: сильное сгущение воздуха не исправляет её форму. Сравнивайте \lambda, Um, баланс потока и повторное открытие точных FEM/ANS, а не только гладкость линий.
В линейной задаче удвоение \Delta Az удваивает B, \Phi и Um, сохраняя \Lambda и \lambda. Удвоение глубины b удваивает \Phi и \Lambda, сохраняя B, Um и \lambda. При подобном изменении всех размеров сечения \lambda также сохраняется; однако масштаб B зависит от того, удерживается ли \Delta Az или Um.
Книжная таблица 5.2 уточняет \lambda от 2.464 до 2.492 для её размеров. У собственной геометрии курса мелкая сетка дала \Phi \approx 39.9992 мкВб, Um≈631.451 А и \lambda \approx 2.520414. Эти результаты относятся к отдельной нативной серии fields; сравнение с 2.492 не является проверкой её точности.
Разобранный численный пример
Точный прямоугольный калибровочный участок
Дано. Две идеальные плоскости разделены \delta =1 мм; ширина потока a=10 мм, глубина b=20 мм; Um=300 А.
Откройте отдельную программу fields встроенным Python FEMM 5; проверьте координаты CDEMN и русские имена свойств.
Нанесите на схему \Delta Az=0.002 Вб/м, b=20 мм, заданный поток 40 мкВб и неизвестное Um. Катушки и конечной \mu стали в этой задаче нет.
Сопоставьте три сетки и отдельно аппроксимации дуги; повторно откройте сохранённую модель и решение.
Получите \Phi из Bn и b\Delta Az, Um из H; вычислите \Lambda, \lambda и относительные расхождения методов.
Удвойте отдельно \Delta Az и b. Обработайте отношения результатов и таблицу сходимости; каждой изменённой геометрии соответствует собственная квитанция.
Ожидаемый результат. Протокол краевой задачи, таблица \Phi/Um/\Lambda /\lambda и два опыта масштабирования.
Как проверить. В \Phi нет лишнего \mu _{0}; знаки и единицы объявлены; \lambda относится к точной геометрии и подтверждён сходимостью.
Обработка и анализ результатов
Какая величина задана разностью A, а какая определяется решением смешанной краевой задачи?
Сведите разность потенциала, глубину, поток, магнитное напряжение двух допустимых путей и проводимость сеточной серии.
Сравните контрольную величину в неизменной внутренней области и оцените влияние изменённой границы.
С чем сравнить. Используйте граничные условия и симметрию из теории занятия. Удалённая граница — приближение; симметрия допустима лишь при согласовании геометрии, источников и искомого режима.
Итог работы. Таблица нормировок и сходимости с объяснением идеальных полюсных границ и отличия этой задачи от полного аппарата.
Программа проверяет аналитические тождества через assert. Она не строит модель и не заменяет результаты полевого расчёта этого занятия.
Python · встроенный MicroPython FEMM 5
"""Самостоятельный аналитический контроль занятия.Среда: встроенный MicroPython FEMM 5; используется только стандартная арифметика.Этот код проверяет формулы и единицы; нативные поля имеют отдельные квитанции."""importmathmu0=4e-7*math.pi;a=0.01;gap=0.001;depth=0.02;um=300# Прямоугольный участок — отдельный точный калибратор, не ответ области CDEMN.b=mu0*um/gap;flux=b*a*depth;da=flux/depthlam=flux/(mu0*depth*um)assertabs(lam-10)<1e-12assertabs(flux-75.3982236862e-6)<1e-12# В собственной области П4 полный поток задаётся краевой разностью Az.assertabs(depth*0.002-40e-6)<1e-15assertabs(flux-depth*da)<1e-15print("Калибровочный поток, Вб",flux,"Лямбда",lam)
Место занятия в расчёте магнитной системы
Воздушная краевая задача полюс–плоскость. Получите Φ из ΔA, Uм из интеграла H и безразмерную λ; независимо проверьте сетку, пути интегрирования и масштабирование.
Самостоятельные воздушные краевые задачи и полный П-магнит с поворотным якорем. Координаты, ось поворота, зазоры и нормировки задаются программой и протоколом геометрии. Откройте модель, воспроизведите указанный режим и сопоставьте свои измерения с таблицей результатов.
Рисунок 39.5. Собственная учебная схема геометрии и измеряемых величин. Размеры, токовые пути, источники и группы заданы программой встроенного Python.
Собственные 2D магнитостатические адаптации методик книги. Проверенные поля FEMM 5, не 3D/физические испытания аппарата.
Самостоятельная работа
Обоснуйте сетку по интегральным величинам и покажите влияние отдельного уточнения дуги.
Критерий выполнения: В \Phi нет лишнего \mu _{0}; знаки и единицы объявлены; \lambda относится к точной геометрии и подтверждён сходимостью.
Первоисточники и дополнительное чтение по занятию
[100] — С. 227–229, 287–296. Формулы 3.33–3.34 на с. 134 исправлены по B=rot A на с. 174 и собственному примеру П4. Печатная нумерация; ссылка ведёт на фрагмент.
Пояснения и типичные ошибкиПочему здесь нет катушки?
Возбуждение задаётся значениями A_z на границах; задача определяет проводимость идеализированного воздушного поля.
Означает ли A_z =\mathrm{const}, что H=0?
Нет: постоянство вдоль границы ограничивает направление B, но нормальная производная A_z может быть ненулевой.
На что обратить внимание
Не используйте книжное \lambda как универсальную константу.
Идеальная \mu \to \infty не позволяет судить о насыщении стали.
