Статическая и динамическая индуктивности: определения, энергия и расчёт в FEMM
Цель и задачи
Цель. Различать секущую и дифференциальную индуктивности по данным модели.
Задачи
Определять статическую и динамическую индуктивности через \Psi (I), объяснять секущую и касательную и доказывать совпадение в линейной системе без смещения.
Получать потокосцепление обмотки в FEMM без повторного умножения на число витков; выбирать допустимый способ расчёта индуктивности.
Если одно из понятий незнакомо, вернитесь к предыдущему занятию и воспроизведите его численный пример. После этого переходите к новой постановке.
Место занятия в инженерной методике
Расчётная лабораторная работа · Определение потокосцеплений, индуктивности и запаса энергии
Инженерный вопрос. Когда два определения совпадают и почему расходятся при насыщении?
Потокосцепления, индуктивность и запас энергии превращают полевую модель в параметры электрической цепи. Они нужны для выбора числа витков, тока, магнитной связи и рабочего диапазона до расчёта динамики.
Последовательность решения, отчёт и критерии результата
Результат получают из наблюдаемых величин и проверяют в условиях данного практикума.
Последовательность решения
Обосновать постановку и способ измерения. Выберите секущее, дифференциальное или взаимное определение параметра; запишите условия применимости энергетической формулы.
Ответить на вопрос по своим результатам. Рассчитайте требуемый параметр по определениям занятия и сопоставьте предусмотренные независимые оценки при одном рабочем состоянии.
Что включить в отчёт
Результат этого занятия. Две зависимости индуктивности с рабочими точками, шагом и энергетическими ограничениями.
Данные именно этого опыта. Сведите ток, потокосцепление и локальные приращения воздушной катушки и насыщенного магнитопровода.
Оформление наблюдений. Токи, направления обмоток, потокосцепления и энергия; для матрицы — независимое возбуждение каждого столбца и проверочный режим.
Как обосновать результат
Проверка из практикума. Сопоставьте секущую индуктивность, локальную производную и энергетическое отношение в одних рабочих точках. Уточните приращение тока и сетку раздельно; сформулируйте границы магнитостатической оценки.
Независимый контроль. Используйте независимую энергетическую оценку, взаимность либо совместное возбуждение в соответствии с линейностью и процедурой занятия.
Критерий выполнения. Определения L_{\mathrm{ст}} и L_{\mathrm{дин}} записаны с фиксированными параметрами; секущая и касательная объяснены геометрически, линейный случай доказан отдельно. Выполнены полные решения при пяти токах на каждой из двух сеток; сохранены .fem, .ans, токи, потокосцепления, энергия и размеры сетки. Для исходных учебных моделей отклонение L_{\mathrm{дин}} от опубликованной оценки на соответствующей сетке не превышает 0,5%; изменение при уменьшении шага менее 0,1%, при данном сгущении сетки менее 0,5%. Это критерии воспроизведения этих моделей, не универсальные допуски FEMM. Для нелинейного магнитопровода отдельно вычислены L_{\mathrm{ст}}, L_{\mathrm{дин}} и 2W/I^{2} и объяснено их различие; для линейной катушки подтверждено совпадение двух определений. Для двух обмоток сохранены оба столбца якобиана, исходная асимметрия и условия взаимности. Аналитический пример, нативные результаты и непроверенные дополнительные режимы разделены.
Точность и расхождения. Уточняйте токовый шаг отдельно от сетки; при нулевом токе не вычисляйте секущее отношение. Для малых взаимных коэффициентов показывайте абсолютное расхождение и знак.
Основание сравнения. Используйте определения секущей, дифференциальной и взаимной индуктивности, приведённые в занятии. Энергетическая формула линейной системы не переносится автоматически на насыщенную систему.
Вывод для устройства
Свяжите найденный параметр с моделью электрической цепи, выбором зазора или допустимым запасом энергии, сохраняя зависимость от рабочего тока.
Границы вывода. Параметр относится к указанной рабочей точке и определению индуктивности. Торцевые поля, сопротивление и нагрев не возникают автоматически из магнитостатической энергии.
Теория и вывод формул
Отношение \Psi /I и производная \partial \Psi /\partial I совпадают в линейной системе без смещения, но имеют разный смысл при насыщении. Геометрия секущей и касательной связывается с матрицей индуктивностей и проверенными расчётами FEMM на двух сетках.
Соглашения об обозначениях и единицах
\(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\). При постоянных магнитах и взаимном возбуждении сначала учитывают ненулевое смещение потокосцепления.
Положительные направления, оси и специальные индексы определяются в конкретной задаче. Единицу берите из условия или заголовка таблицы; перед подстановкой в формулу СИ переводите миллиметры в метры, миллигенри в генри и градусы в радианы.
Численные результаты обычно показаны с четырьмя значащими цифрами; для сравнения близких значений в таблицах сходимости сохранено больше цифр. Полные значения остаются в исходных расчётных данных. Малые остатки вычисления координат на осях показаны как ноль; малые физические величины сохраняются в показательной записи.
Исходная величина — потокосцепление обмотки \Psi, а индуктивность описывает его зависимость от тока. Для тонких последовательных витков \Psi равно алгебраической сумме потоков через ориентированные поверхности витков. По теореме Стокса каждый поток можно вычислить циркуляцией векторного потенциала A по замкнутому контуру соответствующего витка. Контуры и поверхности должны иметь согласованную ориентацию. Отдельный незамкнутый отрезок обмотки не определяет калибровочно независимое потокосцепление.
Частная формула \Psi =N\Phi допустима, когда все N витков сцеплены с одним и тем же ориентированным потоком. При рассеянии потоки отдельных витков различаются. Для распределённой обмотки применяют усреднение с функцией распределения витков; FEMM возвращает уже полное потокосцепление Circuit Properties с учётом назначенных витков и их знаков. Умножать этот результат на N повторно нельзя.
Зафиксируем геометрию x, температуру, материальные характеристики и токи остальных обмоток. В квазистатической однозначной модели получается функция \Psi (I;x). Только после задания этих условий можно говорить об отношении \Psi /I или о производной по I.
\Psi=\sum_{k=1}^{N}\Phi_k=\sum_{k=1}^{N}\int_{S_k}\mathbf B\cdot d\mathbf S=\sum_{k=1}^{N}\oint_{C_k}\mathbf A\cdot d\mathbf l
Равенство потоков — дополнительное условие, а не определение потокосцепления.
Статическая и динамическая: физическое назначение и геометрия графика
Для одной обмотки без независимого смещения потока статическая индуктивность L_{\mathrm{ст}}(I_{0})=\Psi (I_{0})/I_{0} связывает полное потокосцепление с данным током. На графике \Psi (I) это наклон секущей из начала координат в рабочую точку. Динамическая индуктивность L_{\mathrm{дин}}(I_{0})=\partial \Psi /\partial I описывает малое приращение потокосцепления около той же точки; геометрически это наклон касательной. Названия «секущая» и «касательная» сохраняются как математические пояснения к основным названиям.
Терминологическое соглашение курса: основные названия — «статическая» и «динамическая» индуктивности. В литературе для L_{\mathrm{дин}} употребляются «дифференциальная», «инкрементальная», differential, incremental, tangent. В этом курсе слово «динамическая» означает именно локальную производную квазистатической характеристики. Оно само по себе не означает, что выполнен переходный расчёт, учтены частотные потери или получена динамическая петля гистерезиса.
При I_{0}=0 отношение \Psi /I не вычисляют делением на ноль. Если \Psi (0)=0 и функция дифференцируема, его предел равен начальной динамической индуктивности. Если есть ненулевое смещение \Psi (0), такого конечного предела у полного отношения может не быть.
Следствие \Psi =I L_{\mathrm{ст}} при фиксированной геометрии. Подстановка L_{\mathrm{ст}} вместо L_{\mathrm{дин}} теряет второй член.
Почему в линейной системе индуктивность одна
Если материалы линейны, геометрия фиксирована, постоянных магнитов и других независимых источников поля нет, то \Psi =LI. Прямая проходит через начало координат, поэтому её секущая и касательная совпадают во всех точках. Записывать два разных численных параметра для такой катушки нет физического основания. Небольшое расхождение вычисленных значений требует проверки дискретизации и численного дифференцирования.
Различие L_{\mathrm{ст}} и L_{\mathrm{дин}} существенно для нелинейного магнитопровода с характеристикой B–H: материал меняет локальный отклик, а вместе с ним меняется функция \Psi (I). Насыщение часто уменьшает динамическую индуктивность быстрее статической. Это следствие формы кривой в данном диапазоне; универсальное неравенство для произвольных характеристик без дополнительных условий утверждать нельзя.
Перемещение ротора или якоря меняет L(x) даже в полностью линейных материалах. Это зависимость от геометрии, а не нелинейность по току: при каждом фиксированном x сохраняется L_{\mathrm{ст}}=L_{\mathrm{дин}}=L(x). Производная d\Psi /dt при движении дополнительно содержит член, связанный с изменением x.
Даже при одной линейной индуктивности движение создаёт дополнительный вклад в напряжение.
Напряжение, постоянные магниты и смещение другими обмотками
В пассивном соглашении знаков напряжение на квазистатической обмотке равно RI+d\Psi /dt. При неподвижной геометрии и фиксированных остальных токах коэффициент при dI/dt — динамическая индуктивность. ЭДС самоиндукции при согласованном обходе имеет противоположный знак: e=-d\Psi /dt. Указание направления тока и полярности выводов обязательно.
Постоянный магнит или фиксированный ток другой обмотки может создать \Psi см\ne 0. Даже при линейных материалах тогда \Psi =\Psi см+LI. Полное отношение \Psi /I=L+\Psi см/I уже не является собственной индуктивностью, хотя материалы линейны. Чтобы сравнивать статическую и динамическую индуктивности именно токового отклика, для линейной системы сначала отделяют \Psi см=\Psi (0), затем вычисляют [\Psi (I)-\Psi (0)]/I.
В насыщенной системе вычитание отдельного расчёта «только магниты» из общего расчёта не создаёт точной суперпозиции: рабочая проницаемость в этих решениях различается. Конечная разность [\Psi (I)-\Psi (0)]/I описывает средний отклик на выбранном интервале, а не локальную динамическую индуктивность. Для малого возмущения нужны производные при сохранённых магнитах и остальных источниках.
Отделение смещения точно для линейной системы при одних и тех же геометрии и независимых источниках.
Энергия и коэнергия: почему 2W/I^{2} не универсальная индуктивность
Для неподвижной обратимой однозначной катушки без независимых источников и с нулевым начальным потоком энергию получают интегрированием I по \Psi, коэнергию — интегрированием \Psi по I. На графике эти площади дополняют друг друга до прямоугольника I\Psi. При линейной зависимости обе площади треугольные и одинаковые.
В нелинейном магнитопроводе энергия W=\int _{\Omega }\int _{0}^{B} H(b)\cdot db dV, тогда как коэнергия W'=\int _{\Omega }\int _{0}^{H} B(h)\cdot dh dV. Выражение B\cdot H/2 допустимо для линейного обратимого материала без смещения; его нельзя безусловно подставлять при насыщении. Соответственно, 2W/I^{2} и 2W'/I^{2} в общем случае не равны ни L_{\mathrm{ст}}, ни L_{\mathrm{дин}}.
Формула 2W/I^{2} для линейной катушки требует энергии всего магнитного поля. Если выбран только сердечник, воздух исключён или условие открытой границы описывает энергию за пределами расчётной области, проверка неполна. Основным способом получения индуктивности обмотки в FEMM служит потокосцепление цепи, а энергетический способ — независимая проверка при выполненных предпосылках.
Локальная кривизна коэнергии; второе дифференцирование энергии численно чувствительнее разности потокосцеплений.
Несколько обмоток: динамическая индуктивность становится матрицей
В системе из m обмоток изменение тока j меняет потокосцепления всех обмоток k. Поэтому динамическая индуктивность — якобиан карты \Psi (I_{1},\ldots ,I_{m},x). Столбец j получают, изменяя только I_j и считывая все \Psi_k. Остальные токи, включая нулевые токи разомкнутых вторичных обмоток, остаются фиксированными. Изменение тока вторичной цепи из-за нагрузки задаёт другую задачу и другой входной отклик.
В линейном случае после отделения независимого потока получается постоянная матрица L(x), связывающая вектор токов с вектором потокосцеплений. Для её столбцов используют отдельные возбуждения. Деление \Psi _{1} на I_{1} при одновременно ненулевом I_{2} включает взаимный вклад и не выделяет L_{11}.
Для обратимой магнитостатической модели с гладкой коэнергией матрица динамических индуктивностей симметрична, поскольку смешанные производные коэнергии совпадают. Конечный шаг и сетка дают малую асимметрию; проверяют её уменьшение, а не исправляют результаты принудительным усреднением до выяснения причины. У многомерной нелинейной системы один вектор \Psi и один вектор I не определяют единственную «статическую матрицу»: требуется явно указанный способ разложения.
(\mathbf L_{\mathrm{дин}})_{kj}=\left.\frac{\partial\Psi_k}{\partial I_j}\right|_{x,I_{\ell\ne j}},\qquad\delta\boldsymbol\Psi=\mathbf L_{\mathrm{дин}}\,\delta\mathbf I
Для движущейся системы отдельно учитывают изменение потокосцепления с положением; другие изменяемые параметры требуют своих членов.
Расчёт в FEMM: линейная катушка
Задайте Frequency=0, фиксированную геометрию и один ненулевой ток цепи. Убедитесь, что в назначенных областях нет нелинейных B–H-кривых и постоянных магнитов; наличие неиспользованного материала в библиотеке модели само по себе не меняет поле. Сохраните рабочую копию, выполните анализ и откройте соответствующее решение.
В Circuit Properties считайте ток и Flux Linkage. Во встроенном Python FEMM команда current, voltage, linkage = mo_getcircuitproperties("Coil") возвращает именно эти три величины. Для линейной модели без смещения L=linkage/current. Число витков уже учтено. В осесимметричной задаче результат относится к полному телу вращения; дополнительно умножать на 2\pi r нельзя. В плоской задаче проверьте заданную глубину и число одинаковых секций.
Для энергетической проверки выделите все блоки магнитного поля командами mo_clearblock(); mo_groupselectblock(), затем получите W=mo_blockintegral(2) и W'=mo_blockintegral(17). Сравните \Psi /I, 2W/I^{2} и 2W'/I^{2}. Повторите расчёт с удвоенным током: \Psi должно удвоиться, энергия увеличиться вчетверо, индуктивность остаться прежней. Точность согласования оцените вместе с сеткой и открытой границей.
Расчёт в FEMM: нелинейная катушка и выбор шага
Сначала решите полную нелинейную задачу при рабочем токе I_{0}. Получите L_{\mathrm{ст}}=\Psi (I_{0})/I_{0} при отсутствии независимого смещения. Затем решите две новые полные нелинейные задачи с токами I_{0}-\Delta I и I_{0}+\Delta I. Все остальные источники, B–H-кривые, геометрия и граничные условия сохраняются. Эти два решения должны заново учитывать изменение насыщения; простое изменение тока без анализа не меняет результат.
Повторите пару при \Delta I/2. Разность центральных оценок характеризует влияние шага. Для гладкой характеристики ошибка центральной разности имеет порядок \Delta I^{2}; чрезмерно малый шаг усиливает влияние допуска нелинейного решателя и вычитания близких чисел. Начните, например, с 1–5% характерного тока, а окончательный шаг выбирайте по наблюдаемой устойчивости. Около нулевого тока используйте абсолютный масштаб, а не процент от нуля.
Отдельно сгущайте сетку и сравнивайте L_{\mathrm{ст}} и L_{\mathrm{дин}}. Внутри серии по току желательно сохранять одну и ту же сетку: смена узлов может загрязнять разность. Одинаковое число элементов не доказывает совпадение сеток. В приведённой ниже проверке выполнения в FEMM дополнительно подтверждено совпадение SHA-256 координат узлов и связности элементов для пяти токов каждой серии. Два шага и две сетки — проверка чувствительности; сами по себе они не устанавливают строгую верхнюю границу ошибки.
Загружаемый скрипт для встроенного Python выполняет пять полных решений для каждой заданной сетки, сохраняет рабочие копии и сырые значения \Psi. Для второго столбца матрицы смените изменяемую цепь на Secondary и задайте её рабочий ток I_{0}=0, сохранив Primary=3 А в исходной модели. При I_{0}=0 скрипт не вычисляет \Psi /I, но центральная производная остаётся определённой.
Для матрицы изменяют один ток, но считывают потокосцепления всех цепей.
Incremental и Frozen: разные способы использования рабочего состояния
В FEMM 5 поле Previous Solution задаёт сохранённый файл .ans базового решения. Режим Incremental использует рабочую точку для малосигнальной задачи с дифференциальным откликом материала. Режим Frozen сохраняет проницаемость базового состояния. Поэтому слова «сохранённая проницаемость» недостаточно: необходимо назвать режим и величину, которая была сохранена.
Для изотропного нелинейного материала даже малосигнальный векторный отклик зависит от направления возмущения относительно рабочего поля; он не сводится к одной скалярной величине B/H во всех направлениях. Incremental учитывает локальную линеаризацию, тогда как Frozen удобен для разложения вкладов источников при зафиксированном распределении проницаемости. Коэффициенты Frozen не следует автоматически подписывать как L_{\mathrm{дин}}.
Во встроенном интерфейсе команда mi_setprevious("base.ans",1) выбирает Incremental, значение 2 — Frozen. Файл base.ans должен соответствовать указанным рабочим токам, положению и сетке; геометрия и сетка для связанной задачи должны быть согласованы с ним. Решение с единичным возмущающим током после линеаризации означает единицу малого сигнала, а не полный новый режим насыщенной машины.
Основной воспроизводимый опыт этой страницы использует полные нелинейные центральные разности. Дополнительное задание — получить тот же якобиан через Incremental и сопоставить результаты при выбранной частоте. Режим Frozen исследуется отдельно как разложение источников; равенство его коэффициентов динамическим здесь не заявляется. Работа конкретных настроек Incremental/Frozen в скачиваемом опыте не проверялась.
Гистерезис и переменный ток: границы термина «динамическая»
Обычная заданная в магнитостатической модели B–H-кривая — однозначная материальная зависимость. В реальном материале с гистерезисом \Psi зависит от предыстории и состояния; центральная разность по независимым магнитостатическим решениям не воспроизводит малую петлю. Для такой задачи нужны явно выбранная ветвь, модель состояния и подходящее экспериментальное или динамическое обоснование. Начальную кривую намагничивания нельзя автоматически считать малосигнальной характеристикой при любой истории.
При Frequency≠0 ток, напряжение и потокосцепление обычно комплексны; вихревые токи, скин-эффект и магнитные потери меняют фазу и распределение поля. Отношение комплексных амплитуд \widehat{\Psi }/\widehat I может быть частотно зависимым комплексным коэффициентом. Это не вещественная производная DC-характеристики. В протокол включают частоту, рабочее подмагничивание, амплитуду сигнала и соглашение о фазорах.
Для эквивалентной последовательной цепи при соглашении \exp (j\omega t) можно определить L_{\mathrm{сер}}=\operatorname{Im}(Z)/\omega, где Z=\widehat U/\widehat I, если \omega >0 и модель действительно допускает такое представление. Сопротивление \operatorname{Re}(Z) включает потери, предусмотренные постановкой. При насыщении конечный синусоидальный сигнал и малый сигнал около DC-состояния — разные задачи; результат гармонического приближения не доказывает форму переходного процесса.
Эквивалентная последовательная индуктивность на заданной частоте и при заданном соглашении о фазорах; не универсальная замена L_{\mathrm{дин}} DC.
Проверенный линейный опыт FEMM 5: одна индуктивность
16 октября 2026 г. заново решены рабочие копии учебной модели solenoid.fem: осесимметричная воздушная катушка, I_{0}=1\,\mathrm{А}, \Delta I=0{,}1\,\mathrm{А} и \Delta I/2=0{,}05\,\mathrm{А}. На исходной сетке 7081 элемент получены \Psi =95{,}851415\,\mathrm{мкВб} и L_{\mathrm{ст}}=L_{\mathrm{дин}}=95{,}851415\,\mathrm{мкГн}. Совпадение двух способов здесь следует из линейности.
При уменьшении заданных размеров элементов у меток вдвое получилось 18000 элементов. На этой сетке \Psi =96{,}138188\,\mathrm{мкВб}, L_{\mathrm{ст}}=L_{\mathrm{дин}}=96{,}138188\,\mathrm{мкГн}; изменение с сеткой составило 0,2983%. Полная полевая энергия W=48{,}072919 мкДж даёт 2W/I^{2}=96{,}145838\,\mathrm{мкГн}, что отличается от \Psi /I на 0,00796%. Это численное расхождение энергетического интеграла и потокосцепления, а не две физические индуктивности линейной катушки.
Исходные .fem и .ans сохранены без изменений. Это фактические новые решения FEMM headless engine v_{5}.0 06Oct2026 на двух конечных сетках, а не аналитический результат бесконечно длинного соленоида и не полная физическая валидация катушки.
Проверенный нелинейный опыт FEMM 5: различие более чем в пять раз
В модели nonlinear-3.fem задан плоский двухобмоточный магнитопровод глубиной 20 мм, N_{1}=100, N_{2}=50, I_{1}=3\,\mathrm{А} и I_{2}=0. Назначенный сердечнику материал NonlinearSteel содержит восемь точек B–H: (0 Тл;0 А/м), (0,25;100), (0,5;200), (0,9;400), (1,2;800), (1,45;2000), (1,6;5000), (1,7;10000). В расчётных областях постоянные магниты не назначены. Геометрию, полярности витков и граничные условия берите из файла, не восстанавливайте по одной строке параметров.
На исходной сетке 11476 элементов: \Psi _{1}=0{,}02865061625\,\mathrm{Вб}, L_{\mathrm{ст}}=9{,}550205\,\mathrm{мГн}, центральная оценка L_{\mathrm{дин}}=1{,}767301\,\mathrm{мГн} при \Delta I=0{,}1\,\mathrm{А} и 1,766 мГн при \Delta I=0{,}05\,\mathrm{А}. Изменение шага составляет 0,04971%. Отношение 2W/I^{2}=4{,}857707\,\mathrm{мГн} не совпадает ни с L_{\mathrm{ст}}, ни с L_{\mathrm{дин}}, как и предсказывает нелинейная теория.
На уточнённой сетке 22602 элемента: \Psi _{1}=0{,}02868537280\,\mathrm{Вб}, L_{\mathrm{ст}}=9{,}561791\,\mathrm{мГн}, L_{\mathrm{дин}}=1{,}769631\,\mathrm{мГн} при \Delta I=0{,}1\,\mathrm{А} и 1,769 мГн при \Delta I=0{,}05\,\mathrm{А}. Изменение шага — 0,04948%, изменение L_{\mathrm{дин}} между сетками — 0,13188%. Статическая индуктивность примерно в 5,406 раза больше динамической. Эти две оценки отвечают на разные вопросы: полное потокосцепление и малый дополнительный отклик.
Дополнительная серия с изменением только I_{2} вокруг нуля даёт второй столбец якобиана. На уточнённой сетке получена приведённая ниже матрица. Взаимные элементы отличаются приблизительно на 0,0144% от их среднего значения: это остаточная численная асимметрия при конечном шаге, а не нарушение взаимности обратимой модели. Скрипт и квитанция сохраняют исходные несреднённые результаты.
Коэнергия W'=\int _{0}^{I}\Psi (i)di=L_{0}I_{s}^{2}\ln (\cosh (I/I_{s}))=0{,}468483297\,\mathrm{Дж}. Энергия W=I\Psi -W'=0{,}354038392\,\mathrm{Дж}.
Энергетическое отношение 2W/I_{0}^{2}=78{,}6752\,\mathrm{мГн}. Оно не равно ни 91,3913 мГн, ни 50,3969 мГн.
Центральная разность даёт 50,41 мГн при \Delta I=0{,}1\,\mathrm{А} и 50,4 мГн при \Delta I=0{,}05\,\mathrm{А}. Ошибка относительно точной производной уменьшается примерно вчетверо.
Проверка:W+W'=I\Psi =0{,}822521688\,\mathrm{Дж}. При I\to 0 обе индуктивности стремятся к L_{0}=120\,\mathrm{мГн}. Числа этого примера относятся к аналитической \tanh-кривой виджета, а не к B–H-материалу модели nonlinear-3.fem.
Практикум в FEMM 5
Сравните две индуктивности на воздушной катушке и насыщенном магнитопроводе
Скачайте модели solenoid.fem и nonlinear-3.fem по ссылкам в источниках. Сохраните оригиналы отдельно; для сгущения создайте новые копии. Для каждой метки уменьшите положительный заданный размер элемента вдвое, сохранив геометрию, материалы, витки и границы.
Скачайте inductance_embedded-ru.py с русскими пояснениями. Укажите абсолютные пути MODELS и новый каталог OUTPUT. Запустите текст скрипта во встроенной Python-консоли FEMM 5. Для первой нелинейной серии оставьте VARY="Primary", I0=3, DELTA=0.1, READ_CIRCUITS=["Primary","Secondary"].
Скрипт сохранит отдельную рабочую .fem до изменения тока, решит пять состояний на каждой сетке и создаст samples.json и summary.json. Проверьте, что исходная модель не содержит Previous Solution: эта работа использует самостоятельные полные нелинейные решения.
Вычислите L_{\mathrm{ст}}=\Psi _{1}/3 и L_{\mathrm{дин}} по центральным разностям при двух шагах. Сопоставьте с опубликованными значениями на соответствующей сетке; сравните отдельно изменение шага и изменение сетки.
Повторите с VARY="Secondary", I0=0 и новым OUTPUT. В исходной копии Primary должен остаться равным 3 А. Соберите матрицу: строки — считываемые потокосцепления, столбцы — изменяемые токи.
Для линейной серии укажите исходную и уточнённую solenoid.fem, VARY="Coil", READ_CIRCUITS=["Coil"], I0=1. Сравните L_{\mathrm{ст}}, L_{\mathrm{дин}} и 2W/I^{2}; объясните остаточное численное расхождение энергии.
Для дополнительного исследования подготовьте связанную Incremental-задачу от сохранённого рабочего .ans и сравните с центральной производной. Укажите частоту и амплитуду возмущения. Frozen исследуйте как отдельное разложение; эту дополнительную часть нельзя отметить выполненной только на основании основной квитанции.
Ожидаемый результат. Уточнённая линейная модель: L_{\mathrm{ст}}=L_{\mathrm{дин}}\approx 96{,}1382\,\mathrm{мкГн}. Уточнённая нелинейная модель: L_{\mathrm{ст}}\approx 9{,}56179\,\mathrm{мГн}, L_{\mathrm{дин},11}\approx 1{,}76876\,\mathrm{мГн}; два определения отличаются существенно.
Как проверить. Для основной нативной проверки выполнены 30 решений: три серии по десять состояний; отдельно ещё 10 решений воспроизвели скачиваемый встроенный скрипт. Исходные файлы сохранены без изменений. Сравнение относится только к этим моделям, сеткам и DC-постановкам.
Обработка и анализ результатов
Когда два определения совпадают и почему расходятся при насыщении?
Сведите ток, потокосцепление и локальные приращения воздушной катушки и насыщенного магнитопровода.
Рассчитайте требуемый параметр по определениям занятия и сопоставьте предусмотренные независимые оценки при одном рабочем состоянии.
С чем сравнить. Используйте определения секущей, дифференциальной и взаимной индуктивности, приведённые в занятии. Энергетическая формула линейной системы не переносится автоматически на насыщенную систему.
Итог работы. Две зависимости индуктивности с рабочими точками, шагом и энергетическими ограничениями.
Расчётные исследования по теме
Изучите физическую постановку, выполните программу во встроенном Python FEMM 5 и проверьте результаты по сохранённым моделям и независимым критериям.
Насыщение изменяет связь тока и потокосцепления. Постоянная индуктивность перестаёт описывать одновременно установившийся поток, запасённую энергию и малое приращение напряжения.
Снимки интерфейса FEMM показывают геометрию, настройки и работу с результатами в FEMM 5. Откройте изображение, чтобы рассмотреть подписи и значения в полном размере.
Насыщение дросселя: статическая индуктивность и энергия. Геометрия учебной модели; подписи скрыты для читаемости. Открыть в полном размереСохранённое решение «Насыщение дросселя: статическая индуктивность и энергия». Модуль магнитной индукции |B|, Тл, и линии статического магнитного поля; рабочая область модели. Открыть в полном размере
Самостоятельная работа
Выполните восемь самостоятельных расчётов из блока вопросов и основную работу FEMM для обеих моделей. Представьте таблицу L_{\mathrm{ст}}, L_{\mathrm{дин}}(\Delta I), L_{\mathrm{дин}}(\Delta I/2), 2W/I^{2}, количества элементов и относительных изменений; для нелинейной модели добавьте полный якобиан. Объясните, какая индуктивность входит в u-RI при неподвижной геометрии, и отдельно разберите смещение магнитом и движение.
Критерий выполнения: Определения L_{\mathrm{ст}} и L_{\mathrm{дин}} записаны с фиксированными параметрами; секущая и касательная объяснены геометрически, линейный случай доказан отдельно. Выполнены полные решения при пяти токах на каждой из двух сеток; сохранены .fem, .ans, токи, потокосцепления, энергия и размеры сетки. Для исходных учебных моделей отклонение L_{\mathrm{дин}} от опубликованной оценки на соответствующей сетке не превышает 0,5%; изменение при уменьшении шага менее 0,1%, при данном сгущении сетки менее 0,5%. Это критерии воспроизведения этих моделей, не универсальные допуски FEMM. Для нелинейного магнитопровода отдельно вычислены L_{\mathrm{ст}}, L_{\mathrm{дин}} и 2W/I^{2} и объяснено их различие; для линейной катушки подтверждено совпадение двух определений. Для двух обмоток сохранены оба столбца якобиана, исходная асимметрия и условия взаимности. Аналитический пример, нативные результаты и непроверенные дополнительные режимы разделены.
Первоисточники и дополнительное чтение по занятию
[33] — Разделы 2.2, 2.3, 3.3.4–3.3.5; поиск mo_getcircuitproperties и mi_setprevious
Пояснения и типичные ошибкиСамостоятельно: \Psi =0{,}02I Вб, I=4\,\mathrm{А}. Найдите L_{\mathrm{ст}}, L_{\mathrm{дин}}, W и W'.
L_{\mathrm{ст}}=L_{\mathrm{дин}}=20\,\mathrm{мГн}; \Psi =0{,}08\,\mathrm{Вб}; W=W'=0{,}16\,\mathrm{Дж}. Линейная характеристика без смещения: секущая и касательная совпадают.
Самостоятельно: \Psi =0{,}03+0{,}01I Вб из-за постоянного магнита, I=2\,\mathrm{А}. Почему \Psi /I не равно собственной индуктивности?
\Psi =0{,}05\,\mathrm{Вб}, \Psi /I=25\,\mathrm{мГн}, но токовая индуктивность равна 10 мГн. После отделения смещения [\Psi (I)-\Psi (0)]/I=10\,\mathrm{мГн} и d\Psi /dI=10\,\mathrm{мГн}; различие полного отношения вызвано смещением, а не нелинейностью.
Самостоятельно: \Psi (2{,}9\,\mathrm{А})=0{,}02847\,\mathrm{Вб} и \Psi (3{,}1\,\mathrm{А})=0{,}02883\,\mathrm{Вб}, \Psi (3\,\mathrm{А})=0{,}02865\,\mathrm{Вб}. Найдите две индуктивности.
L_{\mathrm{ст}}=0{,}02865/3=9{,}55\,\mathrm{мГн}. L_{\mathrm{дин}}\approx (0{,}02883-0{,}02847)/0{,}2=1{,}80\,\mathrm{мГн}. Это учебная округлённая таблица; для контроля шага нужны дополнительные решения, а не только эти три точки.
Самостоятельно: для предыдущей округлённой таблицы при dI/dt=100\,\mathrm{А/с} и неподвижном якоре найдите u-RI.
u-RI\approx L_{\mathrm{дин}}\cdot dI/dt=0{,}0018\cdot 100=0{,}18\,\mathrm{В}. Подстановка L_{\mathrm{ст}} дала бы 0,955 В и неверный малосигнальный отклик.
Самостоятельно: линейная катушка имеет L(x)=0{,}02+2x Гн, где x в метрах; I=3\,\mathrm{А}, dI/dt=0, dx/dt=0{,}1\,\mathrm{м/с}. Каково u-RI?
При любом фиксированном x статическая и динамическая индуктивности равны. Напряжение движения u-RI=I(dL/dx)(dx/dt)=3\cdot 2\cdot 0{,}1=0{,}6\,\mathrm{В}. Переменная L(x) не означает нелинейность по току.
Самостоятельно: в линейной системе L_{11}=10\,\mathrm{мГн}, M_{12}=4\,\mathrm{мГн}, I_{1}=2\,\mathrm{А}, I_{2}=3\,\mathrm{А}. Что даёт \Psi _{1}/I_{1}?
\Psi _{1}=0{,}01\cdot 2+0{,}004\cdot 3=0{,}032\,\mathrm{Вб}; \Psi _{1}/I_{1}=16\,\mathrm{мГн}. Собственная индуктивность по-прежнему 10 мГн. Чтобы выделить её, задают I_{2}=0 в линейной модели или берут \partial \Psi _{1}/\partial I_{1} при фиксированном I_{2}.
Самостоятельно: 100 витков сцеплены с 0,2 мВб каждый, а ещё 100 — с 0,1 мВб каждый. Найдите \Psi и L_{\mathrm{ст}} при 2 А.
\Psi =100\cdot 0{,}2+100\cdot 0{,}1=30\,\mathrm{мВб}; L_{\mathrm{ст}}=15\,\mathrm{мГн}. Формула 200·0,2 мВб ошибочно приписала бы всем виткам один поток. При выводе Circuit Properties 30 мВб повторное умножение на 200 запрещено.
Самостоятельно: Z=2+j6{,}283185\,\mathrm{Ом} при 50 Гц. Найдите эквивалентную последовательную индуктивность и поясните её статус.
L_{\mathrm{сер}}=6{,}283185/(2\pi \cdot 50)\approx 0{,}02\,\mathrm{Гн}=20 мГн. Это параметр последовательной эквивалентной цепи на 50 Гц при заданных условиях; его нельзя автоматически объявить производной DC-кривой.
Почему при насыщении нельзя вычислять любую индуктивность как 2W/I^{2}?
В нелинейном случае W — интеграл I по \Psi, а не треугольник LI^{2}/2. Энергетическое отношение не совпадает в общем случае ни с \Psi /I, ни с d\Psi /dI. В выполненном нелинейном опыте три числа различаются.
Почему один и тот же I_{0} ещё не задаёт динамическую индуктивность реального материала с гистерезисом?
Нужны состояние и предыстория намагничивания, выбранная ветвь или малая петля и модель её отклика. Однозначная B–H-кривая FEMM не содержит произвольную историю реальной петли.
Какая проверка различает влияние шага и влияние сетки?
На одной сетке сравнить центральные разности при \Delta I и \Delta I/2; затем при том же рабочем состоянии повторить обе оценки на уточнённой сетке. Сохранять сырые \Psi, допуск решателя и состав фиксированных источников.
Почему коэффициенты Frozen нельзя без проверки назвать динамическими индуктивностями?
Frozen фиксирует проницаемость базового состояния для выбранного разложения источников. Динамическая индуктивность определяется производной полного отклика; Incremental использует локальную линеаризацию. Эти операции различаются в нелинейной модели.
На что обратить внимание
Называть \Psi /I и d\Psi /dI двумя различными индуктивностями линейной катушки без смещения.
Повторно умножать Flux Linkage из Circuit Properties на N или на 2\pi r.
Подставлять 2W/I^{2} при насыщении, независимом поле магнита или неполном объёме интегрирования.
Изменять положение ротора или ток другой обмотки одновременно с токовым возмущением и называть результат частной производной.
Считать конечный ненулевой ток возмущения «маленьким» без проверки \Delta I/2 или игнорировать смену сетки.
Отождествлять Frozen и Incremental или считать магнитостатическую B–H-кривую полной моделью гистерезиса.
Сравнивать аналитический виджет \tanh с нативной моделью другого магнитопровода как один и тот же эталон.
