Если одно из понятий незнакомо, вернитесь к предыдущему занятию и воспроизведите его численный пример. После этого переходите к новой постановке.
Место занятия в инженерной методике
Расчётная лабораторная работа · Построение карты отклика и локальных параметров
Инженерный вопрос. Достаточно ли точны производные карты около рабочей точки?
Потокосцепления, индуктивность и запас энергии превращают полевую модель в параметры электрической цепи. Они нужны для выбора числа витков, тока, магнитной связи и рабочего диапазона до расчёта динамики.
Локальный Якобиан потокосцеплений. Выполните действия практикума по его исходной таблице. Сравните столбцы Якобиана при двух приращениях и независимый малый токовый сдвиг.
Ответить на вопрос по своим результатам. Проверьте карту в независимых точках; сопоставьте интерполяцию или производную с прямым расчётом и повторите при меньшем шаге.
Что включить в отчёт
Результат этого занятия. Матрица с указанной точкой, шагом и ошибкой предсказания малого изменения.
Данные именно этого опыта. Сравните столбцы Якобиана при двух приращениях и независимый малый токовый сдвиг.
Оформление наблюдений. Координаты каждой точки, отклик, шаг производной и способ интерполяции; контрольные точки не должны совпадать с узлами определения параметра.
Как обосновать результат
Проверка из практикума. Для двух портов требуется четыре возмущённых решения; для N портов — 2N плюс базовый режим.
Независимый контроль. Сопоставьте прогноз с прямым расчётом в независимой точке и измените шаг производной или плотность карты отдельно от сетки.
Критерий выполнения. Все токи и угол обозначены; шаг проверен; матрица динамических индуктивностей не заменена отношениями потокосцеплений к токам; взаимность оценена численно.
Точность и расхождения. Разделите невязку полевого решения, ошибку интерполяции и чувствительность производной к приращению. Вблизи нулевого отклика сохраняйте абсолютную оценку с единицей.
Основание сравнения. Используйте определения карты, локальной производной и коэнергии из занятия. Вывод относится к проверенному диапазону; экстраполяция и переход через границу карты требуют отдельного обоснования.
Вывод для устройства
Объясните, в каком диапазоне карта может заменить прямые полевые расчёты в модели цепи или движения и почему за её пределами нужен новый расчёт.
Границы вывода. Параметр относится к указанной рабочей точке и определению индуктивности. Торцевые поля, сопротивление и нагрев не возникают автоматически из магнитостатической энергии.
Дополнительный расчётный метод
Зубцовый контур и взаимная проводимость при движении
Зубцовый контур служит элементом матрицы магнитной связи; серия положений показывает, как взаимная проводимость изменяется при движении.
Таблицы потокосцеплений описывают взаимное влияние нескольких обмоток в нелинейной магнитной системе. Их локальное дифференцирование позволяет определить матрицу индуктивностей в выбранной рабочей точке.
Соглашения об обозначениях и единицах
\(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(i + \Delta i) \approx \Psi(i) + L_{\text{дин}} \Delta i выводит эту матрицу непосредственно. Для каждого столбца меняется только один ток; остальные токи и положение фиксируются.
Центральная разность использует симметричные возмущения и обычно уменьшает ошибку усечения относительно односторонней. Но чрезмерно малый шаг усиливает влияние численного шума поля. Поэтому проверяют \Delta I и \Delta I/2, сохраняя одну сетку и одинаковые допуски решателя. Взаимность L_{jk} = L_{kj} относится к локальным производным одной консервативной точки. Покомпонентное деление \Psi_j/i_k не создаёт универсальную матрицу статических индуктивностей насыщенной многопортовой машины; присутствие постоянных магнитов дополнительно требует корректного учёта исходного потокосцепления.
В линейной системе \boldsymbol\Psi=\boldsymbol\Psi_0+\mathbf L\mathbf i одна постоянная матрица L описывает и конечные изменения относительно \boldsymbol\Psi_0, и локальные производные. При нелинейной B\text{–}H-характеристике именно матрица L_{\text{дин},jk}=\partial\Psi_j/\partial i_k однозначно задаёт малый отклик в фиксированной точке; универсальная статическая матрица не получается покомпонентным делением. Слова «касательная» и «дифференциальная» далее объясняют математический смысл динамической матрицы.
\Delta I — шаг тока, А; элементы локальной матрицы имеют единицу Гн.
Порядок проверки результата
Численный пример ниже является аналитическим учебным эталоном. Его значения вычислены по указанным формулам; они не обозначают результат решения новой модели методом конечных элементов. Практикум связывает эталон с полевой задачей: сначала воспроизведите уже проверенный результат базового раздела, затем измените только один исследуемый параметр. Именно такое сравнение позволяет отделить ошибку программы от влияния геометрии, сетки и физического приближения. Сохраните исходную модель, рабочую копию и таблицу с заголовками величин. В подписи каждой величины укажите единицу, положительное направление и способ получения; для гармонических величин дополнительно укажите амплитудную нормировку. Если условие задачи не выполнено, остановите серию и сохраните причину отказа вместе с последним корректным состоянием.
Пояснение программы
Функция linkage1 задаёт независимую нелинейную характеристику первой обмотки. Во всех возмущённых состояниях ток второй обмотки равен 1 \,\text{А}. Производная по первому току вычисляется центральной разностью с расстоянием 2\cdot\mathtt{step} между точками. Два значения шага позволяют проверить устойчивость производной. В системе нескольких связанных обмоток такую проверку выполняют для каждого столбца локальной матрицы индуктивностей.
Термин «дифференциальная индуктивность» в комментариях программы означает динамическую индуктивность L_{\text{дин}}: здесь вычисляется центральная оценка частной производной.
Таблица потокосцеплений как локальная модель
Каждая строка таблицы относится к одной геометрии, температуре и положению. Все компоненты потокосцепления снимаются для одного совместного токового состояния.
Постройте таблицу в пространстве независимых токов. При постоянных магнитах сохраните ненулевое смещение при нулевом токе.
В одной выбранной точке используйте парные возмущения каждого тока для столбцов Якобиана. Порядок токов и потокосцеплений должен совпадать во всех строках.
Проверьте предсказание в случайном малом смешанном направлении и затем уменьшите его масштаб. Ошибка линейного предсказания должна уменьшаться до уровня численного шума.
Границы применимости. Интерполяция допускается только внутри охваченного диапазона. Вне его программа должна явно сообщить об экстраполяции, особенно при насыщении.
Связь с общими определениями
Определения и границы применимости собраны в занятии Статическая и динамическая индуктивности. Используйте интегральную запись как исходную; алгебраические формулы применяйте после проверки перечисленных там условий.
Разобранный численный пример
Вычисленный аналитический эталон
Дано. Независимый математический эталон: \Psi_{1} = 0{,}1i_{1} + 0{,}02i_{2} + 0{,}002i_{1}^{3}; \Psi_{2} = 0{,}02i_{1} + 0{,}08i_{2}. Точка (i_{1},i_{2}) = (2{,}1) А; \Delta I = 0{,}1 \,\text{А}. Положительный кубический член даёт растущую динамическую индуктивность и не моделирует насыщение железа. Эталон проверяет только разностный алгоритм; насыщенная полевая карта исследуется отдельно.
В двухобмоточной модели задайте рабочие токи и сохраните базовый результат.
Для каждой цепи отдельно выполните два симметричных токовых возмущения и читайте \Psi обеих обмоток.
Соберите столбцы матрицы из центральных разностей. Между столбцами восстанавливайте базовые токи.
Повторите вычисления с половинным шагом. Отдельно оцените расхождение взаимных элементов; не заменяйте их средним значением до проверки.
Расширьте токовый набор до небольшой карты; вне рассчитанной области запретите интерполяцию.
Укажите наличие постоянных магнитов и тип материала; не выдавайте отношение потокосцепления к току за динамическую индуктивность.
Ожидаемый результат. Матрица в одной рабочей точке и контроль зависимости от шага; карта содержит все исходные \Psi.
Как проверить. Для двух портов требуется четыре возмущённых решения; для N портов — 2N плюс базовый режим.
Обработка и анализ результатов
Достаточно ли точны производные карты около рабочей точки?
Сравните столбцы Якобиана при двух приращениях и независимый малый токовый сдвиг.
Проверьте карту в независимых точках; сопоставьте интерполяцию или производную с прямым расчётом и повторите при меньшем шаге.
С чем сравнить. Используйте определения карты, локальной производной и коэнергии из занятия. Вывод относится к проверенному диапазону; экстраполяция и переход через границу карты требуют отдельного обоснования.
Итог работы. Матрица с указанной точкой, шагом и ошибкой предсказания малого изменения.
Этот эталон проверяет разностный алгоритм; полевые \Psi подставляются после независимого повторения GUI-модели.
Python · встроенный MicroPython FEMM 5
"""Аналитическая проверка центральной производной.Среда: встроенный MicroPython FEMM 5.Аналитический учебный пример; полевой решатель не вызывается.Параметры, единицы величин и порядок вывода пояснены по стадиям."""# Учебная нелинейная характеристика первой обмотки; токи в А, потокосцепление в Вб.deflinkage1(i1,i2):return0.1*i1+0.02*i2+0.002*i1**3# Уменьшаем возмущение тока вдвое для независимой проверки шага производной.forstepin[0.1,0.05]:# Центральная разность по первому току при неизменном втором токе 1 А.# Расстояние между двумя режимами равно 2·step; результат имеет размерность Гн.derivative=(linkage1(2+step,1)-linkage1(2-step,1))/(2*step)# Порядок вывода: шаг тока, А; дифференциальная индуктивность, Гн.print(step,derivative)
Расчётное исследование и разбор результатов
Эти материалы продолжают практикум: постановка, программа встроенного Python, сохранённые численные результаты и проверка их смысла. Графики строятся по результатам FEMM 5.
Магнитная система с несколькими обмотками описывается не одной величиной L, а отображением вектора токов в вектор потокосцеплений. При насыщении это отображение нелинейно; рабочая точка определяет отклик на малое возмущение.
Фрагмент таблицы Р4.2. Сохранённые результаты FEMM 5
Рисунок 17.5. Визуализация численных данных FEMM 5: Нелинейный магнитопровод. Это экспорт результатов поля, а не снимок окна программы.
Условная учебная кривая B(H), заданная в файле; это не характеристики промышленной стали. Статическая индуктивность (наклон секущей) и динамическая индуктивность (наклон касательной) при насыщении различаются.
Ориентиры в интерфейсеРисунок 17.6. Геометрия: двухобмоточный линейный сердечник, независимое первичное возбуждение Открыть в полном размереРисунок 17.7. Распределение магнитной индукции: двухобмоточный линейный сердечник, независимое первичное возбуждение Открыть в полном размереРисунок 17.8. Учебная нелинейная B–H-характеристика: подтверждение введённых точек и их единиц. Открыть в полном размереРисунок 17.9. Поле нелинейного двухобмоточного сердечника при первичном токе 3 А. Открыть в полном размереРисунок 17.10. Встроенная консоль FEMM: Таблицы потокосцеплений и локальная матрица индуктивностей. Аналитическая проверка программы занятия. Показана консоль исходной редакции программы, зафиксированной при съёмке. Текущая программа с русскими пояснениями проверена отдельно во встроенном Python FEMM; этот снимок не является новым запуском текущей редакции. Открыть в полном размере
Самостоятельная работа
Постройте матрицу в двух рабочих точках с разным насыщением. Сравните собственные и взаимные производные, сохранив исходные таблицы.
Критерий выполнения: Все токи и угол обозначены; шаг проверен; матрица динамических индуктивностей не заменена отношениями потокосцеплений к токам; взаимность оценена численно.
