Если одно из понятий незнакомо, вернитесь к предыдущему занятию и воспроизведите его численный пример. После этого переходите к новой постановке.
Место занятия в инженерной методике
Расчётная лабораторная работа · Выбор наблюдаемых величин и проверка интегрального закона
Инженерный вопрос. Почему точка, профиль, контур и область дают разные результаты?
Эти действия входят в подготовку расчёта любого устройства: они устанавливают связь источника, геометрии и наблюдаемой величины. Их результат позволяет выбрать допустимую постановку и способ извлечения будущей характеристики.
Последовательность решения, отчёт и критерии результата
Результат получают из наблюдаемых величин и проверяют в условиях данного практикума.
Последовательность решения
Обосновать постановку и способ измерения. Определите ориентацию контура или поверхности, охватываемые источники и геометрический множитель. Закон связывает измеряемую величину с независимым источником.
Поток и потокосцепление. Проверьте знаки и единицы интегрального закона на разобранной геометрии.
Четыре способа измерения. Выполните действия практикума по его исходной таблице. Извлеките предусмотренные практикумом точечные, контурные и объёмные величины.
Результат этого занятия. Таблица величин, единиц и областей измерения без смешения плотности и интеграла.
Данные именно этого опыта. Извлеките предусмотренные практикумом точечные, контурные и объёмные величины.
Оформление наблюдений. Контур или область интегрирования, ориентация, точечные и интегральные величины раздельно, источник и невязка закона.
Как обосновать результат
Проверка из практикума. Каждая величина имеет единицу и определение области; для полной энергии не пропущен воздух.
Независимый контроль. Сопоставьте интеграл с источником или суммой потоков, а не с повторным чтением того же интеграла другим именем.
Критерий выполнения. Другой учащийся может воспроизвести измерения по таблице без угадывания выбранных областей.
Точность и расхождения. Нормируйте невязку на физический масштаб источника или потока, объявленный до сравнения. Для баланса около нуля приводите также абсолютную невязку в единицах закона.
Основание сравнения. Используйте интегральный закон и ориентацию контура из теории занятия. Плотность, интеграл по сечению и полный объёмный результат сравнивают после учёта геометрического множителя.
Вывод для устройства
Объясните, даёт ли выбранный интеграл полный параметр устройства или только локальную характеристику его сечения.
Границы вывода. Простой эталон проверяет конкретное действие или физический закон; перенос на устройство требует заново согласовать геометрию, материалы, режим и объёмный множитель.
Дополнительный расчётный метод
Векторная алгебра и векторный магнитный потенциал
Провод с током даёт количественный эталон B(r), A(r), циркуляции и потока. Через него проверяют калибровку потенциала, локальные производные и допустимость скалярного магнитного потенциала.
Получаемые величины.B_r и B\tau прямого провода; разности A и поток на единицу длины; циркуляция H; ошибка численной производной и сеточной аппроксимации
Различить скалярное и векторное поле и вычисляемые операции.
Сопоставить аналитическое поле внутри и снаружи провода с точками и контурами FEMM.
Изменить отсчёт A, сохранив B и разности потенциала.
Проверить скачок скалярного магнитного потенциала вокруг тока и объяснить необходимость разреза области.
Страницы источника. PDF 1–31 полностью: 1–12 операции и потенциалы; 13–23 поле прямого провода; 24–28 поток, Биo–Савар, уравнения потенциала; 29–31 разностная схема.
Переход от распределения поля к инженерным характеристикам требует выбора измеряемой величины и соответствующего интеграла. Значения в точках, интегралы по контуру и объёмные интегралы рассматриваются с учётом их физического смысла и единиц.
Соглашения об обозначениях и единицах
\(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» без координаты, способа сглаживания и материала не является достаточным описанием: максимум может относиться к искусственному острому углу.
Точечное поле и сглаживание
На линейном треугольном элементе потенциал интерполируется линейно, а его производные кусочно постоянны. Показанное B может быть дополнительно сглажено для удобства просмотра. Сглаживание не добавляет новых физических данных и иногда скрывает скачок около границы материала. Для количественного сравнения используется один режим обработки и фиксированная координата. В точке на самой границе следует уточнять, с какой стороны берётся значение; у идеального острого угла локальное поле может иметь сингулярный характер.
Поток через поверхность
В плоскопараллельной модели отрезок сечения, вытянутый вдоль Depth, образует поверхность. Поток через неё равен интегралу нормальной B по линии, умноженному на глубину. Знак задаётся нормалью, связанной с направлением обхода. Если окно обработки результатов уже учитывает Depth, повторно умножать результат не нужно. В осесимметрии поверхность образуется вращением линии и содержит вес 2\pi r. Интеграл через стержень измеряет поток этого стержня, но не автоматически потокосцепление всех витков.
Для плоской задачи d в м, \mathrm dl в м, B в Тл, итог \Phi в Вб; нормаль n единичная.
Потокосцепление и результаты цепи
Для тонких замкнутых витков \Psi=\sum_k\int_{S_k}\mathbf B\cdot\mathbf n_k\,\mathrm dS=\sum_k\oint_{C_k}\mathbf A\cdot\mathrm d\mathbf l. Контуры ориентированы по положительному току, а нормали — правилом правой руки. Это сумма интегралов отдельных витков; один произвольный контур по сердечнику не заменяет геометрию всей обмотки.
Потокосцепление суммирует потоки, сцепляющиеся с отдельными витками. В распределённой обмотке витки могут сцепляться с разными потоками; для неё полезен расчёт через потенциал и определённую схему намотки. В FEMM значение \Psi получают в Circuit Properties. Оно зависит от витков и их знаков, тогда как локальное поле при отключённой измерительной обмотке может быть неизменным. Поэтому ошибку \Psi часто нужно искать в соединении обмоток, а не в сетке железа.
\Psi=\sum_{q=1}^{N}\Phi_q=\sum_{q=1}^{N}\oint_{C_q}\mathbf A\cdot\mathrm d\mathbf l
Для тонких замкнутых витков \Phi_q — поток через поверхность с границей C_q; ориентация согласована с положительным током. Единица потокосцепления — Вб.
Интегралы по области
Для энергии выделяется вся область, где хранится учитываемое поле, включая воздух. Выделение только железа пропускает энергию рассеяния и зазора. Для силы по взвешенному тензору выделяется тело, на которое действует сила, а не весь магнитопровод без разбора. Для потерь выбирают материал с соответствующим механизмом и корректной частотой. Прежде чем сравнивать два интеграла, проверьте, относятся ли они к одному телу и одной глубине. Название интеграла не отменяет условий применимости его формулы.
W_{\text{м}} в Дж; w_{\text{м}} в \text{Дж}/\text{м}^3; формула плотности ограничена линейной изотропной обратимой связью B=\mu H.
Минимальная таблица результата
В отчёт включаются параметры источника, частота, глубина или осесимметрия, число элементов и определения измерений. Для B указывают координату и компоненту, для потока — поверхность и знак, для \Psi — цепь и витки, для энергии — выбранные области. Можно добавить карту и профиль B вдоль линии, но они должны пояснять числовые ответы. Такой формат позволяет повторить измерение и отличить изменение физики от изменения визуальной шкалы. Для каждого вопроса укажите нужный тип измерения: точка, контур или интеграл по области.
Ориентация измерения сохраняется в таблице
Когда контур проводится слева направо, его нормаль выбирается правилом ориентации, и поток получает соответствующий знак. Разворот того же контура меняет знак, но не модуль интеграла. Аналогично, смена начала обмотки меняет знак её потокосцепления. Для сравнения двух методов нужно согласовать именно эти ориентации, иначе числа будут противоположными при физически одинаковом поле. Энергия не зависит от направления обхода контура, а знак компоненты силы определяется положительным направлением выбранной оси. В таблице результатов следует указывать не только значения, но и определения величин. Например: поток через поперечное сечение левого стержня, нормаль вверх; потокосцепление цепи первой обмотки с указанным началом; сила на отдельный якорь по оси y.
От локального поля к наблюдаемой величине
Каждая интегральная величина требует собственной поверхности, объёма или электрической цепи. Цветовая шкала индукции не является измерением потокосцепления.
Для потока выберите ориентированную поверхность. В плоской постановке её элемент равен длине участка контура, умноженной на глубину; знак определяется выбранной нормалью.
Для потокосцепления сложите потоки всех витков. У распределённой обмотки разные витки могут охватывать разные потоки, поэтому произведение одного локального потока на число витков требует обоснования.
Для энергии интегрируйте плотность по всем областям поля, включая воздух. Уточните, является результат мгновенной статической энергией или величиной, определённой для гармонической задачи.
Границы применимости. При сравнении двух моделей должны совпадать область интегрирования, число витков, глубина и амплитудное соглашение. Максимум поля у идеального острого угла может не иметь конечного сеточного предела.
Пояснения к официальному руководству
Подробности, которые следует учитывать при настройке и проверке модели. Все учебные пояснения.
Сглаживание поля: что восстанавливается и что остаётся приближением
В плоской задаче линейный потенциал на треугольнике даёт постоянное по элементу производное поле. На общей вершине несколько соседних элементов поэтому предлагают разные значения B. Карта без сглаживания показывает этот кусочно-постоянный результат. Сглаживание создаёт узловое представление производного поля для интерполяции и отображения; исходные узловые значения потенциала и сама сетка от этого не становятся решениями более высокого порядка.
Общий метод восстановления по группе соседних элементов мог бы аппроксимировать значения поля в центрах треугольников плоскостью по наименьшим квадратам. Однако на нерегулярной сетке такая задача может оказаться плохо обусловленной. Для магнитного поля FEMM использует более простое среднее: вклад каждого соседнего элемента взвешивается обратно пропорционально расстоянию от его центра тяжести до рассматриваемого узла. Поэтому большие и маленькие элементы влияют не одинаково, а восстановление не является простым средним всех ближайших значений.
Каждому соседнему треугольнику приписывается вес, обратный расстоянию от его центра до рассматриваемого узла. Сглаживание делает поле непрерывнее на рисунке, но не создаёт новую сетку.
Обозначения и единицы
\mathbf B_e — индукция в центре элемента, Тл; \mathbf x_e и \mathbf x_{\mathrm{node}} — координаты центра и узла, м; w_e — вес, м⁻¹; восстановленная \mathbf B_{\mathrm{node}} имеет единицу Тл.
Условия применения. На интерфейсах действуют специальные правила; нельзя усреднять несовместимые предельные значения через физический разрыв.
В PDF правило задано словами «inversely proportional to the distance»; нормированная формула раскрывает это описание.
Для электростатики процедура другая: в окрестности узла аппроксимируют узловые напряжения плоскостью и используют её градиент для восстановления электрического поля и смещения с учётом материала. Специальные случаи на границах и поверхностях обрабатываются отдельно. Магнитную формулу усреднения B поэтому нельзя выдавать за точное описание любого поля любого модуля.
На границе материалов физически непрерывны определённые компоненты, а не весь вектор. Без поверхностного заряда непрерывна нормальная D; при разных \varepsilon нормальная E обычно разрывна. Аналогично для границы раздела магнитных сред без поверхностного тока непрерывны нормальная B и тангенциальная H. Сглаживание изображения не должно скрывать эти физически обусловленные скачки. Сравнивайте значения с одной выбранной стороны границы раздела и фиксируйте режим сглаживания.
У острых углов поле может иметь сингулярность. Сглаживание способно изменить показанный максимум, но не создать конечное физическое значение для идеальной острой кромки. Для оценки изоляции указывают реальный радиус скругления; для интегральных сил и энергии отдельно проверяют сходимость нужного интеграла. Сравнение сглаженного и несглаженного профиля полезно как диагностический опыт, но не является самостоятельной строгой оценкой погрешности.
Дано. Однородная B=0{,}5 \,\text{Тл} нормально пересекает стержень шириной 10 \,\text{мм} и глубиной 20 \,\text{мм}. Все 100 витков сцепляются с этим потоком.
Создайте плоскопараллельную модель в миллиметрах с Depth=20 мм и Frequency=0. Постройте два концентрических квадрата: наружный x,y=\pm 30 \,\text{мм}, внутренний x,y=\pm 20 \,\text{мм}. Область между ними — один линейный магнитопровод с \mu_r=1000 и \sigma=0.
Обмотка 1 охватывает левый стержень: прямоугольники x=-40\ldots -30 и x=-20\ldots -10 \,\text{мм}, в обоих y=-10\ldots 10 \,\text{мм}. Назначьте цепь «Обмотка 1» типа series и соответственно +100 и -100 витков.
Обмотка 2 охватывает правый стержень: прямоугольники x=10\ldots 20 и x=30\ldots 40 \,\text{мм}, y=-10\ldots 10 \,\text{мм}. Для положительного взаимного потокосцепления назначьте цепь «Обмотка 2» и соответственно -50 и +50 витков. Подпишите на схеме выбранный положительный обход потока.
Замкните внешнюю квадратную область x,y=-150\ldots 150 \,\text{мм} и назначьте Prescribed A=0. Метки воздуха необходимы во всех независимых замкнутых воздушных областях, включая окно. Материал обмоток: \mu_r=1, \sigma=0, без дополнительного источника J.
Задайте размер элемента в железе и обмотках не более 2 \,\text{мм}. Сначала возбудите только одну цепь, затем другую; каждый раз сохраняйте отдельную копию задачи и оба потокосцепления из View → Circuit Properties.
При I_{1}=1 \,\text{А},I_{2}=0 снимите B в середине левого стержня. Проведите поперёк этого стержня контур и измерьте нормальный поток; подпишите направление контура.
Из результатов цепи снимите \Psi_{1} и \Psi_{2}. Для энергии выделите всё железо, обмотки и воздух, а затем отдельно только железо; сравните значения.
Постройте профиль B вдоль линии через окно и стержни. Сравните сглаженный и несглаженный вид, не меняя исходную модель.
Ожидаемый результат. Поток стержня и \Psi/N близки при небольшом рассеянии, но не обязаны совпадать точно. Полная энергия не меньше энергии только железа.
Как проверить. Каждая величина имеет единицу и определение области; для полной энергии не пропущен воздух.
Обработка и анализ результатов
Почему точка, профиль, контур и область дают разные результаты?
Извлеките предусмотренные практикумом точечные, контурные и объёмные величины.
Приведите результаты к одинаковым единицам и ориентации; вычислите невязку независимых способов измерения.
С чем сравнить. Используйте интегральный закон и ориентацию контура из теории занятия. Плотность, интеграл по сечению и полный объёмный результат сравнивают после учёта геометрического множителя.
Итог работы. Таблица величин, единиц и областей измерения без смешения плотности и интеграла.
Пошаговая настройка расчётного примера
Откройте отдельную страницу с последовательными действиями, точными параметрами, накопительными файлами этапов и программами с русскими комментариями.
Рисунок 4.1. Визуализация численных данных FEMM 5: Двухобмоточный магнитопровод: первичное возбуждение. Это экспорт результатов поля, а не снимок окна программы.
Двумерная модель активной части. Модель стали с постоянной магнитной проницаемостью не описывает насыщение. Идеализированные обмотки с нулевой проводимостью не учитывают все потери. Сходимость силы и момента требует отдельной проверки.
Ориентиры в интерфейсеРисунок 4.2. Интерфейс FEMM 5: осесимметричная катушка, геометрические объекты, метки и воздушная область. Открыть в полном размереРисунок 4.3. Поле воздушной катушки: осесимметричный магнитостатический эталон FEMM 5. Открыть в полном размереРисунок 4.4. Результаты цепи Coil: заданный ток и рассчитанное потокосцепление. Открыть в полном размереРисунок 4.5. Интегральная энергия всей магнитной системы; выделение включает катушку и окружающий воздух. Открыть в полном размереРисунок 4.6. Измерение B и H в точке (−25; 15) мм внутри левого стержня. Открыть в полном размереРисунок 4.7. Ориентированный контур поперёк стержня и интеграл нормальной индукции. Открыть в полном размереРисунок 4.8. Профиль |B| через стержни и окно: сглаживание выключено. Контур от (−45; 15) до (45; 15) мм, 400 точек. Position — расстояние от начала контура в мм; |B| — магнитная индукция в Тл. Единицы восстановлены по модели и выбранной величине: само окно FEMM их не печатает. Открыть в полном размереРисунок 4.9. Профиль |B| через стержни и окно: сглаживание включено. Контур от (−45; 15) до (45; 15) мм, 400 точек. Position — расстояние от начала контура в мм; |B| — магнитная индукция в Тл. Единицы восстановлены по модели и выбранной величине: само окно FEMM их не печатает. Открыть в полном размере
Самостоятельная работа
Подготовьте таблицу B, \Phi, \Psi, W_{\text{м}} для одной модели с точным описанием всех измерений.
Критерий выполнения: Другой учащийся может воспроизвести измерения по таблице без угадывания выбранных областей.
