Если одно из понятий незнакомо, вернитесь к предыдущему занятию и воспроизведите его численный пример. После этого переходите к новой постановке.
Место занятия в инженерной методике
Расчётная лабораторная работа · Построение карты отклика и локальных параметров
Инженерный вопрос. Как выбор поля H внутри участка меняет потребную МДС и результат дробления?
Магнитная цепь, потокосцепление, тяговая характеристика и потери образуют методику выбора зазора, обмотки и режима привода. Сила становится проектным критерием после сопоставления с нагрузкой механизма.
Последовательность решения, отчёт и критерии результата
Результат получают из наблюдаемых величин и проверяют в условиях данного практикума.
Сравните два алгоритма и причину ошибки. Выполните действия практикума по его исходной таблице. Сохраните потоки и индукции на границах участков, оба расчёта магнитных напряжений и независимые полевые сечения.
Ответить на вопрос по своим результатам. Проверьте карту в независимых точках; сопоставьте интерполяцию или производную с прямым расчётом и повторите при меньшем шаге.
Что включить в отчёт
Результат этого занятия. Две согласованные таблицы алгоритмов и сравнение с FEMM при одной B(H), геометрии и выбранном рабочем потоке.
Данные именно этого опыта. Сохраните потоки и индукции на границах участков, оба расчёта магнитных напряжений и независимые полевые сечения.
Оформление наблюдений. Координаты каждой точки, отклик, шаг производной и способ интерполяции; контрольные точки не должны совпадать с узлами определения параметра.
Как обосновать результат
Проверка из практикума. Оба алгоритма используют одни данные; полушаг и полный баланс показаны; сходимость не объявлена доказательством правильной карты поля.
Независимый контроль. Сопоставьте прогноз с прямым расчётом в независимой точке и измените шаг производной или плотность карты отдельно от сетки.
Критерий выполнения. Оба алгоритма используют одни данные; полушаг и полный баланс показаны; сходимость не объявлена доказательством правильной карты поля.
Точность и расхождения. Разделите невязку полевого решения, ошибку интерполяции и чувствительность производной к приращению. Вблизи нулевого отклика сохраняйте абсолютную оценку с единицей.
Основание сравнения. Используйте определения карты, локальной производной и коэнергии из занятия. Вывод относится к проверенному диапазону; экстраполяция и переход через границу карты требуют отдельного обоснования.
Вывод для устройства
Объясните, в каком диапазоне карта может заменить прямые полевые расчёты в модели цепи или движения и почему за её пределами нужен новый расчёт.
Границы вывода. Магнитная сила отдельного привода не подтверждает работу всего аппарата: контактные усилия, пружины, трение, динамика и нагрев должны быть заданы и проверены по назначению.
Теория и вывод формул
Расчёт по участкам восстанавливает требуемую МДС по полю зазора и распределению потока в полюсе. Сравниваются приближение постоянной напряжённости на полушаге и уточнение по промежуточному сечению.
Соглашения об обозначениях и единицах
\(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\). При постоянных магнитах и взаимном возбуждении сначала учитывают ненулевое смещение потокосцепления.
Положительные направления, оси и специальные индексы определяются в конкретной задаче. Единицу берите из условия или заголовка таблицы; перед подстановкой в формулу СИ переводите миллиметры в метры, миллигенри в генри и градусы в радианы.
Численные результаты обычно показаны с четырьмя значащими цифрами; для сравнения близких значений в таблицах сходимости сохранено больше цифр. Полные значения остаются в исходных расчётных данных. Малые остатки вычисления координат на осях показаны как ноль; малые физические величины сохраняются в показательной записи.
Заданы желаемая B_{\delta} однородной части зазора, геометрия, безгистерезисная H(B) и проводимости боковых трубок. Искомая величина — NI. В половинной симметричной схеме \delta означает половину полного промежутка. Множитель 2 возвращают только при сборке полного контура. Границы трубок согласуют с поперечными сечениями полюса.
По мере удаления от узкого торца боковое рассеяние увеличивает поток, так что \Phi k+1>\Phi k. При этом индукция может как возрастать, так и уменьшаться, поскольку одновременно изменяется площадь Sk. Такое распределение рабочих точек не позволяет описывать сильно насыщенный полюс одной средней проницаемостью \mu. В пошаговом расчёте Uk обозначает магнитное напряжение, накопленное от плоскости симметрии зазора до начала рассматриваемого участка.
Кривая H(B) данного П6-образца отличается от кривой прежней серии П-магнита 38.2: H(1,45 Тл)=1000 А/м здесь и 600 А/м там. Для сравнения двух алгоритмов П6 используется одна и та же кривая именно П6; межсерийные различия геометрии и материала не выдают за ошибку одного метода.
\Phi 0 может включать торцевое выпучивание и потому не обязано совпадать с B_{\delta}, умноженным на геометрическую площадь торца.
Алгоритм 1: боковой поток сосредоточен в середине
На входе участка известны Uk, \Phi k, Sk, \ell k и \Lambda k бокового потока, условно привязанного к середине. Сначала на первой половине длины H принимается равной входной Hk. Полученное напряжение середины возбуждает боковую трубку. После добавления потока вычисляются B и H на следующем сечении, а полное падение на участке определяется трапецией по H концов.
Следующий участок использует обновлённые \Phi k+1 и Uk+1. Это дискретная адаптация метода П6, а не конечный элемент. Участки особенно важно сгущать там, где H быстро меняется. Все величины должны вычисляться по одной и той же H(B).
Для предварительного расчёта середины разделим \Lambda k на две равные части. Первую питаем напряжением Uk, получая промежуточный поток \Phi k*. По нему и промежуточной площади находим Hk*. Падение на первой половине вычисляем средним Hk и Hk*, а не постоянным входным Hk. Затем всю \Lambda k питаем исправленным напряжением середины и обновляем конец участка теми же формулами, что в методе 1.
Если боковая проводимость распределена неравномерно, \Lambda k/2 заменяют проводимостью фактической первой полудлины. Промежуточную площадь берут по геометрии. Равное деление всех трубок без анализа их расположения может ухудшить результат.
У трубки, пересекающей катушку, Um определяется охваченным током и падениями в стали. В первом приближении при равномерной намотке и слабонасыщенном ярме геометрическая \Lambda i приводится по потоку с долей охваченных витков \alpha i. Приведение относится к напряжению своего контура, а не автоматически к U_{0} зазора.
Сложите падения на полюсе, зазорах, ярме и обратном тракте одного замкнутого пути. При заданном NI вместо B_{\delta} потребуется внешняя итерация по потоку до нулевой невязки. Сравнение с полем проводится на одинаковых контурах: распределение отдельных интегралов H по деталям может зависеть от выбранного пути.
Все члены относятся к одному физическому замкнутому пути.
Проверка применимости
У автора первый способ обычно завышал, второй несколько занижал МДС рассмотренных систем. Это наблюдение не является теоремой: в другой дискретной схеме знак разности может измениться. Оба метода сравнивают с одинаковыми H(B), геометрией и картой трубок. Уменьшение длины участка не исправляет неверную картину рассеяния.
При сильном насыщении ярма нельзя считать его почти эквипотенциальным. При значительном заполнении окна катушкой меняется охват её тока трубками. Требуется уточнение связанной схемы или поля. Указанные в книге 5 Т относятся к специальным конструкциям и сопоставлениям автора, а не к общему пределу точности.
В новой серии курса при целевой B_{\delta}=1.25 Тл половинная МДС поля равна примерно 583.889 А, а оба метода с 64 участками дают около 599.87 А: общая ошибка 2.74%, но ошибка отдельно выделенного падения на полюсе достигает 31.25%. При 1.45 Тл методы дают примерно 1344.5 А против полевых 973.28 А, то есть около 38.14%. Сходимость дробления не устранила ошибку проводимостей, извлечённых из идеального поля.
Объявите одну таблицу Sk, \ell k, \Lambda k и согласуйте трубки с сечениями. Начните с приведённого одиночного участка.
Выполните построчный расчёт обоих методов: Uk, \Phi k, Bk, Hk, Uc, \Delta \Phi, \Delta U. Метод 2 не заимствует конечное H из метода 1.
Раздробите участки вдвое, распределяя \Lambda по фактическим областям. Сравните полную МДС и максимальную B.
Сопоставьте новую нативную серию при 1.25 и 1.45 Тл: отдельно полную МДС, падение на полюсе и поле ярма.
Обработайте невязку NI, разность методов и сходимость дробления; отдельно оцените ошибку исходной карты рассеяния.
Ожидаемый результат. Таблицы двух алгоритмов, карта участков и сравнение общего и локального падения МДС.
Как проверить. Оба алгоритма используют одни данные; полушаг и полный баланс показаны; сходимость не объявлена доказательством правильной карты поля.
Обработка и анализ результатов
Как выбор поля H внутри участка меняет потребную МДС и результат дробления?
Сохраните потоки и индукции на границах участков, оба расчёта магнитных напряжений и независимые полевые сечения.
Проверьте карту в независимых точках; сопоставьте интерполяцию или производную с прямым расчётом и повторите при меньшем шаге.
С чем сравнить. Используйте определения карты, локальной производной и коэнергии из занятия. Вывод относится к проверенному диапазону; экстраполяция и переход через границу карты требуют отдельного обоснования.
Итог работы. Две согласованные таблицы алгоритмов и сравнение с FEMM при одной B(H), геометрии и выбранном рабочем потоке.
Программа проверяет аналитические тождества через assert. Она не строит модель и не заменяет результаты полевого расчёта этого занятия.
Python · встроенный MicroPython FEMM 5
"""Самостоятельный аналитический контроль занятия.Среда: встроенный MicroPython FEMM 5; используется только стандартная арифметика.Этот код проверяет формулы и единицы; нативные поля имеют отдельные квитанции."""importmatharea=1e-4;phi=80e-6;u=100.0;length=0.01;side=0.2e-6defh(b):return100*b+1000*b**3# Метод 1 использует входное H на первой половине длины.h0=h(phi/area);uc1=u+h0*length/2b1=(phi+side*uc1)/areadu1=length*(h0+h(b1))/2# Метод 2 сначала предсказывает промежуточный поток и его H.hm=h((phi+side*u/2)/area);uc2=u+length*(h0+hm)/4b2=(phi+side*uc2)/area;du2=length*(h0+h(b2))/2assertabs(du1-8.55228673337)<1e-9assertabs(du2-8.57010085112)<1e-9assertb1>phi/areaandb2>phi/areaprint("Падения на участке, А",du1,du2)
Место занятия в расчёте магнитной системы
Два поучастковых алгоритма насыщенного полюса. Проследите присоединение бокового потока и уточнение H на полушаге; сравните общую МДС и её отдельные составляющие с независимым полем.
Поучастковая цепная оценка и независимые поля FEMM; переменное питание с внешним замыканием вторичного контура. Конструкция и допущения указаны отдельно для каждого опыта. Откройте модель, воспроизведите указанный режим и сопоставьте свои измерения с таблицей результатов.
Рисунок 44.4. Собственная учебная схема геометрии и измеряемых величин. Размеры, токовые пути, источники и группы заданы программой встроенного Python.
Поучастковая цепная оценка и независимые поля FEMM; переменное питание с внешним замыканием вторичного контура. Конструкция и допущения указаны отдельно для каждого опыта.
Самостоятельная работа
Покажите случай, когда суммарная МДС выглядит точной из-за малого вклада ошибочного участка; оцените ошибку самого участка.
Критерий выполнения: Оба алгоритма используют одни данные; полушаг и полный баланс показаны; сходимость не объявлена доказательством правильной карты поля.
Первоисточники и дополнительное чтение по занятию
[100] — С. 30–43, 319–324, приложение П6. Печатные 319–324 соответствуют PDF-страницам 320–325. Изложение и одиночный численный пример курса самостоятельные.
Пояснения и типичные ошибкиВ чём разница методов?
В оценке H на первой половине участка: входное значение либо среднее с предварительным промежуточным сечением.
Устраняют ли 1000 участков неучтённое рассеяние?
Нет: это уточнение дискретизации прежней схемы, а не восстановление исключённых потоков.
На что обратить внимание
Не путайте полный промежуток и его половину.
Разность методов не является гарантированной оценкой ошибки.
