FEMM 5

Часть VI · Электрические аппараты → § 39 · Рассеяние, потокосцепление и магнитная проводимость

Занятие 39.3

Распределённая МДС: от непрерывной обмотки к коэффициентам 1/2 и 1/3

Цель и задачи

Цель. Вывести коэффициенты приведения рассеяния и проверить их применимость.

Задачи

  • Вывести Um(x) из закона полного тока для равномерной обмотки.
  • Отдельно вычислить поток и потокосцепление рассеяния.
Что нужно вспомнить перед началом

Если одно из понятий незнакомо, вернитесь к предыдущему занятию и воспроизведите его численный пример. После этого переходите к новой постановке.

Место занятия в инженерной методике

Расчётная лабораторная работа · Выбор наблюдаемых величин и проверка интегрального закона

Инженерный вопрос. Почему один и тот же слой рассеяния даёт коэффициенты 1/2 и 1/3 для разных величин?

Магнитная цепь, потокосцепление, тяговая характеристика и потери образуют методику выбора зазора, обмотки и режима привода. Сила становится проектным критерием после сопоставления с нагрузкой механизма.

Последовательность решения, отчёт и критерии результата

Результат получают из наблюдаемых величин и проверяют в условиях данного практикума.

Последовательность решения

  1. Обосновать постановку и способ измерения. Определите ориентацию контура или поверхности, охватываемые источники и геометрический множитель. Закон связывает измеряемую величину с независимым источником.

  2. Два способа подсчитать потокосцепление. Проверьте знаки и единицы интегрального закона на разобранной геометрии.

  3. Соберите непрерывную и дискретную эпюры. Выполните действия практикума по его исходной таблице. Сведите профиль МДС, число сцепленных витков, поток, потокосцепление и энергию трёх сеток; отдельно измените распределение витков.

  4. Ответить на вопрос по своим результатам. Приведите результаты к одинаковым единицам и ориентации; вычислите невязку независимых способов измерения.

Что включить в отчёт

Результат этого занятия. Вывод обоих интегралов, таблица полевого эталона и объяснение изменения коэффициентов при неравномерной намотке.

Данные именно этого опыта. Сведите профиль МДС, число сцепленных витков, поток, потокосцепление и энергию трёх сеток; отдельно измените распределение витков.

Оформление наблюдений. Контур или область интегрирования, ориентация, точечные и интегральные величины раздельно, источник и невязка закона.

Как обосновать результат

Проверка из практикума. Есть оба независимых интеграла, размерности и проверка g=t^{2}; 1/2 и1/3 не объявлены универсальными.

Независимый контроль. Сопоставьте интеграл с источником или суммой потоков, а не с повторным чтением того же интеграла другим именем.

Критерий выполнения. Есть оба независимых интеграла, размерности и проверка g=t^{2}; 1/2 и1/3 не объявлены универсальными.

Точность и расхождения. Нормируйте невязку на физический масштаб источника или потока, объявленный до сравнения. Для баланса около нуля приводите также абсолютную невязку в единицах закона.

Основание сравнения. Используйте интегральный закон и ориентацию контура из теории занятия. Плотность, интеграл по сечению и полный объёмный результат сравнивают после учёта геометрического множителя.

Вывод для устройства

Объясните, даёт ли выбранный интеграл полный параметр устройства или только локальную характеристику его сечения.

Границы вывода. Магнитная сила отдельного привода не подтверждает работу всего аппарата: контактные усилия, пружины, трение, динамика и нагрев должны быть заданы и проверены по назначению.

Теория и вывод формул

Распределённая обмотка создаёт изменяющееся вдоль магнитопровода магнитное напряжение. Интегрирование этой зависимости позволяет получить два способа приведения рассеяния и установить границы соответствующих приближений.

Соглашения об обозначениях и единицах
  • \(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\). При постоянных магнитах и взаимном возбуждении сначала учитывают ненулевое смещение потокосцепления.
  • Положительные направления, оси и специальные индексы определяются в конкретной задаче. Единицу берите из условия или заголовка таблицы; перед подстановкой в формулу СИ переводите миллиметры в метры, миллигенри в генри и градусы в радианы.
  • Численные результаты обычно показаны с четырьмя значащими цифрами; для сравнения близких значений в таблицах сходимости сохранено больше цифр. Полные значения остаются в исходных расчётных данных. Малые остатки вычисления координат на осях показаны как ноль; малые физические величины сохраняются в показательной записи.
Указатель формул этого занятия
  1. (39.8) · Что распределено и по какой координате
  2. (39.9) · Вывод коэффициента по потоку
  3. (39.10) · Вывод коэффициента по потоку
  4. (39.11) · Почему в потокосцеплении возникает квадрат
  5. (39.12) · Почему в потокосцеплении возникает квадрат
  6. (39.13) · Поток по длине стержня и связь со схемой
  7. (39.14) · Обобщение и пределы
  8. (39.15) · Связь двух нормировок рассеяния

Что распределено и по какой координате

Рассмотрим идеализированную П-образную магнитную систему, сопротивлением стали в которой пренебрегают. Катушка содержит N витков, равномерно распределённых на длине \ell. Координата x отсчитывается от края, где замкнутый контур рассеяния не охватывает витков, до края полного охвата обмотки. В выбранном слое рассеяния элементу dx соответствует геометрическая проводимость a dx; постоянный коэффициент a измеряется в Гн/м. Это допущение относится к принятой модели слоя и не предполагает постоянства поля в произвольном реальном аппарате.

Контур через идеальную сталь и выбранный воздушный слой охватывает Nx/\ell витков. Поэтому его магнитное напряжение равно части полной МДС. Именно геометрия охвата создаёт линейную эпюру; распределение J по площади и конечная \mu могут её изменить.

F=NI,\quad t=x/\ell,\quad N_x=Nt,\quad U_m(x)=Ft,\quad d\Lambda_g=a\,dx,\quad\Lambda_g=a\ell
(39.8)
\Lambda_g — геометрическая проводимость всего слоя, ещё не приведённая к полной МДС.

Вывод коэффициента по потоку

Поток каждого тонкого слоя равен его проводимости, умноженной на местное магнитное напряжение. Суммируем именно потоки, поэтому дополнительного множителя числа охваченных витков пока нет.

  1. Интегрируем линейную эпюру магнитного напряжения.

    \Phi_\sigma=\int_0^\ell U_m(x)a\,dx=Fa\ell\int_0^1t\,dt=F\Lambda_g/2
    (39.9)
  2. Определяем эквивалентную проводимость, включённую под полную МДС.

    \Lambda_{\sigma\Phi}=\Phi_\sigma/F=\Lambda_g/2
    (39.10)

Почему в потокосцеплении возникает квадрат

Тот же элементарный поток сцеплен только с Nx витками. Нельзя сначала сложить все потоки рассеяния, а потом умножить результат на полное N. Интеграл содержит произведение двух одинаковых долей t: одна относится к возбуждению потока, другая — к охвату витков.

  1. Суммируем потокосцепления тонких слоёв.

    \Psi_\sigma=\int_0^\ell N_x\,d\Phi_\sigma=NF\Lambda_g\int_0^1t^2\,dt=NF\Lambda_g/3
    (39.11)
  2. Получаем проводимость для потокосцепления, индуктивность и энергию линейного рассеяния.

    \Lambda_{\sigma\Psi}=\Lambda_g/3,\quad L_\sigma=N^2\Lambda_g/3,\quad W_\sigma=W_\sigma^\prime=\tfrac12L_\sigma I^2
    (39.12)

Поток по длине стержня и связь со схемой

Пусть рабочий поток \Phi_{\delta} проходит через весь основной путь. Выбранное сечение с координатой x пересекают только те потоки рассеяния, чьи трубки входят в стержень дальше этого сечения. В объявленной ориентации это интеграл от x до \ell. Эпюра потока параболическая, хотя эпюра магнитного напряжения линейна.

В схеме потоков используется \Lambda_{\sigma\Phi}; в энергетической схеме и цепи обмотки используется \Lambda_{\sigma\Psi}. Для рабочего потока, сцепленного со всеми витками, \Psi_{\delta}=N\Phi_{\delta}. У броневого аппарата с внутренним рабочим зазором эпюры и координаты меняются; переносить эту формулу без новой схемы нельзя.

\Phi_c(x)=\Phi_\delta+\tfrac12F\Lambda_g(1-t^2),\quad\Psi=N\Phi_\delta+NF\Lambda_g/3
(39.13)
Поток сечения и \Psi /N отличаются, даже когда сопротивление стали равно нулю.

Обобщение и пределы

Для неравномерной катушки введём g(t)=N охв(t)/N и непостоянную плотность проводимости a(t). При идеальной стали Um=Fg, поэтому два коэффициента становятся взвешенными интегралами g и g^{2}. При конечной \mu Um уже включает падения в стали; тогда его нужно находить из связанной магнитной цепи или поля.

Коэффициенты приведения c_{\Phi}=1/2 и c_{\Psi}=1/3 действуют при g=t, постоянной a и выбранной одинаковой длине слоя. Смещение катушки, широкое сечение, насыщение, внешнее и лобовое рассеяние требуют нового расчёта. Замена широкой катушки тонким токовым слоем сохраняет NI, но может изменить g, поле окна и \Psi.

c_\Phi={\int g(t)a(t)\,dt\over\int a(t)\,dt},\quad c_\Psi={\int g^2(t)a(t)\,dt\over\int a(t)\,dt}
(39.14)
c_{\Phi} и c_{\Psi} — коэффициенты приведения одного геометрического слоя: \Lambda_{\sigma\Phi}=c_{\Phi}\Lambda_g и \Lambda_{\sigma\Psi}=c_{\Psi}\Lambda_g. Они отличаются от коэффициентов рассеяния k_{\Phi}=\Phi_{\text{ст}}/\Phi_{\delta} и k_{\Psi}=\Psi /(N\Phi_{\delta}) из занятия о потокосцеплении.

Связь двух нормировок рассеяния

В занятии о потокосцеплении k_{\Phi}=\Phi_{\text{ст}}/\Phi_{\delta} и k_{\Psi}=\Psi /(N\Phi_{\delta}) сравнивают полный поток или полное потокосцепление с рабочим путём. Здесь c_{\Phi}=\Lambda_{\sigma\Phi}/\Lambda_g и c_{\Psi}=\Lambda_{\sigma\Psi}/\Lambda_g приводят только выбранный слой рассеяния к полной МДС. Поэтому численно сравнивать k с 1/2 или 1/3 нельзя: знаменатели и состав учитываемого потока различны.

Для данной идеализированной схемы \Phi_{\delta}=F\Lambda_{\delta}, \Phi_{\text{ст}}=F\Lambda_{\delta}+F\Lambda_g c_{\Phi}, а \Psi /N=F\Lambda_{\delta}+F\Lambda_g c_{\Psi}. После деления на рабочий поток получается явная связь двух нормировок. Она требует одного и того же слоя и отсутствия неучтённых ветвей; при конечной проницаемости и иной картине потоков определения сохраняются, но эта простая связь может не выполняться.

Числовая сверка разобранного примера: \Lambda_g=6 мкГн, \Lambda_{\delta}=12 мкГн, c_{\Phi}=1/2, c_{\Psi}=1/3. Получаем k_{\Phi}=1.25 и k_{\Psi}=7/6\approx 1.166667. Действительно, поток у ярма 3 мВб, рабочий поток 2.4 мВб, \Psi =1.12 Вб и N\Phi_{\delta}=0.96 Вб: 3/2.4=1.25, 1.12/0.96=7/6. Прежний секционный пример с 50, 55 и 60 мкВб имеет другие исходные потоки и даёт свои k_{\Phi}=1.20, k_{\Psi}=1.10; это не противоречие.

\Lambda_{\sigma\Phi}=c_\Phi\Lambda_g,\quad\Lambda_{\sigma\Psi}=c_\Psi\Lambda_g,\qquad k_\Phi=1+{\Lambda_g\over\Lambda_\delta}c_\Phi,\quad k_\Psi=1+{\Lambda_g\over\Lambda_\delta}c_\Psi
(39.15)
c описывает приведение геометрического слоя, k — полный поток или потокосцепление относительно рабочего пути. Единичное слагаемое относится к рабочему потоку.

Разобранный численный пример

Два способа подсчитать потокосцепление

Дано. N=400, I=0.5 А, \Lambda_g=6 мкГн, \Lambda_{\delta}=12 мкГн; равномерная катушка и идеальная сталь.

  1. МДС равна 200 А; \Phi_{\delta}=2.4 мВб и \Phi \sigma =0.6 мВб.

  2. \Psi_{\delta}=0.96 Вб, \Psi \sigma =0.16 Вб. Проводимость для \Psi равна 12+6/3=14 мкГн.

  3. Полный поток у ярма 3 мВб, но N\Phi ярма=1.2 Вб превышает реальное \Psi.

Ответ: \Psi =1.12 Вб, L=2.24 Гн; L_{\sigma}=0.32 Гн, W_{\sigma}=0.04 Дж. Подстановка N\Phi ярма завышает \Psi на 7.14%.

Проверка: Проверка \Psi /I=N^{2}(\Lambda_{\delta}+\Lambda_g/3) даёт тот же результат; числа являются аналитическим примером.

Практикум в FEMM 5

Соберите непрерывную и дискретную эпюры

Таблица 39.5. Исходные данные практикума
Исходная величинаЗначениеЕдиница
Слои10
N400
I0,5А
\Lambda_g6мкГн
  1. Разбейте t∈[0,1] на 10, 20 и 40 слоёв; используйте середины слоёв для d\Phi и d\Psi.
  2. Постройте таблицу t, Nx, Um, d\Phi, d\Psi; отдельно суммируйте два последних столбца.
  3. Измените g(t) на t^{2} и получите c_{\Phi}=1/3, c_{\Psi}=1/5 при той же a; объясните изменение охвата витков.
  4. В собственной полевой модели сравните \Phi двух сечений и \Psi из цепных свойств. Подпишите, что поле не обязано иметь постоянную a.
  5. Обработайте данные: абсолютные значения, нормированные коэффициенты и изменение с числом слоёв; храните источники и допущения отдельно от протокола выполнения в FEMM.

Ожидаемый результат. Две эпюры, таблица интегрирования и объяснение отличия \Phi ярма от \Psi /N.

Как проверить. Есть оба независимых интеграла, размерности и проверка g=t^{2}; 1/2 и1/3 не объявлены универсальными.

Обработка и анализ результатов

Почему один и тот же слой рассеяния даёт коэффициенты 1/2 и 1/3 для разных величин?

  1. Сведите профиль МДС, число сцепленных витков, поток, потокосцепление и энергию трёх сеток; отдельно измените распределение витков.
  2. Приведите результаты к одинаковым единицам и ориентации; вычислите невязку независимых способов измерения.

С чем сравнить. Используйте интегральный закон и ориентацию контура из теории занятия. Плотность, интеграл по сечению и полный объёмный результат сравнивают после учёта геометрического множителя.

Итог работы. Вывод обоих интегралов, таблица полевого эталона и объяснение изменения коэффициентов при неравномерной намотке.

Встроенный Python: Распределённая МДС: от непрерывной обмотки к коэффициентам 1/2 и 1/3 · Порядок запуска

Проверить вывод и единицы встроенным Python

Программа проверяет аналитические тождества через assert. Она не строит модель и не заменяет результаты полевого расчёта этого занятия.

Python · встроенный MicroPython FEMM 5
"""Самостоятельный аналитический контроль занятия.
Среда: встроенный MicroPython FEMM 5; используется только стандартная арифметика.
Этот код проверяет формулы и единицы; нативные поля имеют отдельные квитанции.
"""
import math
n=400; current=0.5; lg=6e-6; ld=12e-6
phi=0.0; psi=0.0; cells=400
# Один множитель t относится к возбуждению потока, второй — к охвату витков.
for k in range(cells):
    t=(k+0.5)/cells
    dp=n*current*t*lg/cells
    phi+=dp; psi+=n*t*dp
# Проверяем разные приведения по потоку и потокосцеплению.
assert abs(phi/(n*current*lg)-0.5)<1e-12
assert abs(psi/(n*n*current*lg)-1/3)<1e-6
assert abs((n*n*ld+psi/current)-2.24)<1e-5
print("Поток утечки, Вб",phi,"Потокосцепление, Вб",psi)

Место занятия в расчёте магнитной системы

Непрерывный вывод коэффициентов приведения 1/2 и 1/3. Проинтегрируйте распределённые МДС и охват витков; объясните разные веса потока и потокосцепления и их связь с энергией.

Опорная теория: 8.7 · Определение Ψ как суммы сцеплений; 39.1 · Рассеяние аппарата и приведение к полной МДС.

Применение: 44.4 · Учесть напряжение боковых трубок и пересечение обмотки; 44.1 · Связать неполное сцепление с взаимной индуктивностью.

Маршрут магнитных систем в основном оглавлении · Каталог рассчитанных конструкций

Расчётный образец FEMM 5

Локальный плоский эталон равномерного токового слоя с обратным поверхностным током на границе A=0; N=300, ширина 50 мм, высота 40 мм, глубина 20 мм. Откройте модель, воспроизведите указанный режим и сопоставьте свои измерения с таблицей результатов.

Таблица 39.6. Результаты расчётного образца
ВеличинаРезультатЕдиница
Коэффициент приведения c_{\Phi}0,51
Коэффициент приведения c_{\Psi}0,33331
Потокосцепление603,2 × 10⁻⁶Вб
Энергия поля301,6 × 10⁻⁶Дж
Коэнергия301,6 × 10⁻⁶Дж
Конечные элементы725591
Собственная схема вывода и сходимость рассчитанных коэффициентов. Иллюстрация построена по формулам и результатам, не является снимком интерфейса.
Рисунок 39.3. Собственная схема вывода и сходимость рассчитанных коэффициентов. Иллюстрация построена по формулам и результатам, не является снимком интерфейса.

Точный локальный плоский эталон равномерно распределённой МДС. Обратный поверхностный ток на A=0. Скалярные коэффициенты 1/2 и 1/3 относятся к равномерной намотке и принятому полю; не универсальны.

Самостоятельная работа

Сравните равномерную и квадратичную функции охвата, а затем объясните, какой геометрический эксперимент способен нарушить принятое приближение.

Критерий выполнения: Есть оба независимых интеграла, размерности и проверка g=t^{2}; 1/2 и1/3 не объявлены универсальными.
Первоисточники и дополнительное чтение по занятию
  1. [100] — С. 15–29, 39–43, 136–137. Печатная нумерация. Издательская ссылка ведёт на фрагмент; полная предоставленная копия изучена отдельно. Выводы и численные примеры курса самостоятельные.
Пояснения и типичные ошибки
Откуда в \Psi возникает вторая степень t?

Местный поток пропорционален охваченной МДС, а его сцепление дополнительно взвешивается числом охваченных витков.

Почему коэффициент энергии не равен 1/2 по потоку?

Энергия линейной катушки равна I\Psi /2, поэтому требуется коэффициент 1/3 по потокосцеплению, а не 1/2 по потоку.

На что обратить внимание

  • Не умножайте весь поток утечки на N.
  • Не считайте линейность Um доказательством линейности потока вдоль стержня.