FEMM 5

Часть V · Электрические машины → § 34 · Потери и неопределённость

Занятие 34.2

Чувствительность к исходным данным и неопределённость

Цель и задачи

Цель. Оценить устойчивость вывода к исходным данным.

Задачи

  • Выполнять контролируемые возмущения входов и объяснять статистический либо интервальный смысл оценки.
  • Ранжировать входы по контролируемым допустимым возмущениям.
Что нужно вспомнить перед началом

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

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

Расчётная лабораторная работа · Оценка достоверности инженерного вывода

Инженерный вопрос. Какой исходный параметр сильнее всего меняет полезный результат?

Свяжите требование по потоку, ЭДС, моменту, мощности, потерям или нагреву с наблюдаемой величиной и сравнимыми режимами. Такой переход позволяет выбрать необходимую сложность модели машины.

Лабораторная работа:

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

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

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

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

  2. Вычисленный аналитический эталон. Проверьте чувствительность по одному изменяемому фактору и правила распространения для принятой модели.

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

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

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

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

Данные именно этого опыта. Сохраните базовый результат и отдельные изменения геометрии, материала и режима.

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

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

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

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

Критерий выполнения. Каждая оценка имеет источник, предпосылки и единицу; возмущения и цепь проверены; лишние значащие цифры исключены.

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

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

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

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

Границы вывода. Проверенный двумерный вопрос является частью проектирования. Неучтённые торцевые, тепловые, механические и динамические эффекты перечисляются по назначению задачи.

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

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

Соглашения об обозначениях и единицах
  • \(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. (34.5) · Физическая и вычислительная постановка
  2. (34.6) · Чувствительность и распространение неопределённости
  3. (34.7) · Численный пример

Физическая и вычислительная постановка

Относительная чувствительность S_x = (x/q)\partial q/\partial x показывает, на сколько процентов меняется результат при малом процентном изменении входа. Её выводят из линейного разложения q(x + \Delta x) \approx q(x) + (\partial q/\partial x)\Delta x. Для индуктивности магнитной цепи L \propto N^2 S/\delta получаются S_N = 2, S_S = 1, S_\delta = -1. Это удобный независимый эталон для численного прохода.

Неопределённость входов и ошибка дискретизации имеют разные причины. При независимых входных стандартных неопределённостях первого порядка сумма квадратов взвешенных чувствительностей даёт дисперсию результата. Если известны только пределы входов, сумма модулей линейных вкладов даёт локальную оценку первого порядка, а не гарантированную границу нелинейного результата. Гарантированный интервал требует анализа исходной функции на всём заданном диапазоне; он не является доверительным интервалом. Корреляции размеров и материалов требуют перекрёстных членов. При сильной нелинейности дополнительно проверяют конечные отклонения и границы карты; одну производную нельзя переносить на большой диапазон насыщения.

S_x=\frac{x}{q}\frac{\partial q}{\partial x},\quad\frac{u(q)}{|q|}\approx\sqrt{\sum_kS_{x_k}^2\left[\frac{u(x_k)}{x_k}\right]^2}
(34.5)
Последняя формула предполагает независимые стандартные неопределённости и локальную линейность.

Порядок проверки результата

Численный пример ниже является аналитическим учебным эталоном. Его значения вычислены по указанным формулам; они не обозначают результат решения новой модели методом конечных элементов. Практикум связывает эталон с полевой задачей: сначала воспроизведите уже проверенный результат базового раздела, затем измените только один исследуемый параметр. Именно такое сравнение позволяет отделить ошибку программы от влияния геометрии, сетки и физического приближения. Сохраните исходную модель, рабочую копию и таблицу с заголовками величин. В подписи каждой величины укажите единицу, положительное направление и способ получения; для гармонических величин дополнительно укажите амплитудную нормировку. Если условие задачи не выполнено, остановите серию и сохраните причину отказа вместе с последним корректным состоянием.

Пояснение программы

При возведении чувствительности зазора в квадрат знак исчезает. Корень суммы квадратов применим к независимым стандартным неопределённостям. Для локального интервального предела используются сумма модулей и иной смысл исходных оценок. Умножение на 100 выводит проценты. Независимость данных и локальная линейность проверяются отдельно; программа не определяет вероятностную модель исходных данных.

Чувствительность и распространение неопределённости

Линейное распространение применимо при малых неопределённостях и гладкой зависимости результата. Корреляции входов должны быть явно рассмотрены.

  1. Изменяйте каждый параметр в обе стороны при повторном согласовании рабочего режима. Производная при фиксированном токе отличается от производной при фиксированном напряжении.

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

  3. Для связанных входов добавьте ковариационные члены. Корень суммы квадратов отдельных вкладов получается только при отсутствии корреляций.

    u^2(q)\approx\sum_{j,k}\frac{\partial q}{\partial x_j}\frac{\partial q}{\partial x_k}\operatorname{cov}(x_j,x_k)
    (34.6)
Таблица 34.4. Обозначения и единицы
ВеличинаФизический смыслЕдиница СИ
qПроверяемый результат; единица определяется выбором результатапо результату
x_kНезависимый входной параметрпо параметру
u(q)Стандартная неопределённость результатакак у результата

Границы применимости. Численная сходимость и неопределённость исходных данных — разные величины. Систематическую ошибку двумерного приближения нельзя автоматически включить в статистическую сумму как случайный независимый вклад.

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

Вычисленный аналитический эталон

Дано. L \propto S/\delta; относительные стандартные неопределённости S и \delta равны 2 \,\% и 5 \,\%, независимы; N фиксировано.

  1. S_S = 1; S_\delta = -1.

    \frac{u(L)}{L}\approx\sqrt{0{,}02^2+0{,}05^2}=0{,}0538516=5{,}38516\,\%
    (34.7)
  2. u(L)/L = \sqrt{0{,}02^{2} + 0{,}05^{2}} = 0{,}053851648 = 5{,}38516 \,\%.

  3. Если вместо стандартных неопределённостей заданы именно пределы \pm2\,\% для площади и \pm5\,\% для зазора, сумма 2\,\%+5\,\%=7\,\% — лишь оценка первого порядка. Для положительных входов точная монотонность L\propto S/\delta даёт интервал 0{,}98/1{,}05\le L/L_0\le1{,}02/0{,}95, то есть от -6{,}667 \,\% до +7{,}368 \,\%. Стандартные неопределённости из исходного условия не являются такими гарантированными пределами.

Ответ: Стандартная относительная неопределённость 5{,}385 \,\% при заявленных предпосылках; интервальный предел имеет другой смысл.

Проверка: Сеточная ошибка не включена молча в стандартную неопределённость исходных размеров.

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

Чувствительность и происхождение неопределённости

Таблица 34.5. Исходные данные практикума
Исходная величинаЗначениеЕдиница
Входызазор, размеры, B\text{–}H, B_{\text{ост}}, R
Возмущения\pm 1; \pm 2%
Масштаб результатазадан
  1. Выберите два наиболее существенных непрерывных параметра с источниками неопределённости.
  2. Выполните базовую точку и симметричные возмущения каждого входа при неизменных остальных параметрах.
  3. Получите чувствительности и повторите с половинным возмущением, используя проверенную сетку.
  4. Укажите тип неопределённости: стандартная, интервальная или дискретный вариант; не смешивайте их.
  5. Для питания от источника напряжения пересогласуйте цепь после каждого изменения; иначе исследуется другая задача.
  6. Представьте результат с единицей и обоснованным числом значащих цифр, отдельно указав численную сходимость.

Ожидаемый результат. Ключевые входы ранжированы по влиянию; оценка результата имеет определённый математический смысл.

Как проверить. Проверена локальная линейность, независимость либо корреляция входов и повторное согласование режима.

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

Какой исходный параметр сильнее всего меняет полезный результат?

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

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

Итог работы. Таблица чувствительности с происхождением диапазонов и ограниченным выводом.

Программа с русскими пояснениями: Численная проверка выведенного эталона · Порядок запуска

Численная проверка выведенного эталона

Код выполняется во встроенном Python FEMM 5. Он проверяет приведённый аналитический вывод; получение полевых данных выполняется отдельным практикумом с указанной моделью.

Python · встроенный MicroPython FEMM 5
"""Чувствительность к исходным данным и неопределённость.

Среда: встроенный MicroPython FEMM 5.
Аналитический учебный пример; полевой решатель не вызывается.
Параметры, единицы величин и порядок вывода пояснены по стадиям.
"""

import math

# Относительные стандартные неопределённости площади и зазора, доли единицы.
area_relative, gap_relative = 0.02, 0.05
# Локальные относительные чувствительности для L ∝ S/g: +1 по площади и −1 по зазору.
area_sensitivity, gap_sensitivity = 1.0, -1.0
# Корень суммы квадратов применим при независимых исходных неопределённостях.
standard_relative = math.sqrt(
    (area_sensitivity * area_relative) ** 2 + (gap_sensitivity * gap_relative) ** 2
)
# Сумма модулей даёт иной результат — локальный интервальный предел.
# Обе оценки предполагают малые изменения и локальную линейность.
interval_relative = (
    abs(area_sensitivity) * area_relative + abs(gap_sensitivity) * gap_relative
)
print("Стандартная неопределенность, %:", 100 * standard_relative)
print("Локальный интервальный предел, %:", 100 * interval_relative)

Расчётное исследование и разбор результатов

Эти материалы продолжают практикум: постановка, программа встроенного Python, сохранённые численные результаты и проверка их смысла. Графики строятся по результатам FEMM 5.

Пошаговая настройка расчётного примера

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

Все расчётные примеры курса

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

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

Таблица 34.6. Результаты расчётного образца
ВеличинаРезультатЕдиница
Частота0Гц
Активная длина80мм
Ток цепи Фаза А1А
Потокосцепление цепи Фаза А0,001816Вб
Ток цепи Фаза С-2А
Потокосцепление цепи Фаза С-0,2922Вб
Ток цепи Фаза В1А
Потокосцепление цепи Фаза В0,2904Вб
Момент на исходной учебной сетке-1,156Н·м
Энергия по локальному направлению намагниченности35,15Дж
Модуль B в контрольной точке0,8503Тл
Координаты контрольной точки (x, y)56,42; 17,25мм
Момент на проверенной уточнённой сетке-1,152Н·м
Визуализация численных данных FEMM 5: Машина с поверхностными постоянными магнитами. Это экспорт результатов поля, а не снимок окна программы.
Рисунок 34.6. Визуализация численных данных FEMM 5: Машина с поверхностными постоянными магнитами. Это экспорт результатов поля, а не снимок окна программы.

Двумерная модель активной части. Модель стали с постоянной магнитной проницаемостью не описывает насыщение. Идеализированные обмотки с нулевой проводимостью не учитывают все потери. Сеточная проверка момента этой рабочей точки выполнена отдельно; её критерии, методы и границы приведены в протоколе. Здесь энергия получена отдельным интегрированием с фактическим локальным направлением намагниченности. Способ интегрирования приведён в программе шага 11. У острых кромок модели стали с постоянной проницаемостью возможны большие локальные значения B. Они не описывают реальный материал и сами по себе не доказывают насыщение или повреждение. Проверяйте нелинейную характеристику материала и сходимость поля в физически заданной области.

Ориентиры в интерфейсе
Геометрия: поверхностные постоянные магниты и проверенная трёхфазная обмотка
Рисунок 34.7. Геометрия: поверхностные постоянные магниты и проверенная трёхфазная обмотка Открыть в полном размере
Распределение магнитной индукции: поверхностные постоянные магниты и проверенная трёхфазная обмотка
Рисунок 34.8. Распределение магнитной индукции: поверхностные постоянные магниты и проверенная трёхфазная обмотка Открыть в полном размере
Машина с поверхностными постоянными магнитами: выбранная сетка воздушного зазора.
Рисунок 34.9. Машина с поверхностными постоянными магнитами: выбранная сетка воздушного зазора. Открыть в полном размере
Машина с поверхностными постоянными магнитами: момент ротора на выбранной сетке.
Рисунок 34.10. Машина с поверхностными постоянными магнитами: момент ротора на выбранной сетке. Открыть в полном размере
Встроенная консоль FEMM: Чувствительность к исходным данным и неопределённость. Аналитическая проверка программы занятия. Показана консоль исходной редакции программы, зафиксированной при съёмке. Текущая программа с русскими пояснениями проверена отдельно во встроенном Python FEMM; этот снимок не является новым запуском текущей редакции.
Рисунок 34.11. Встроенная консоль FEMM: Чувствительность к исходным данным и неопределённость. Аналитическая проверка программы занятия. Показана консоль исходной редакции программы, зафиксированной при съёмке. Текущая программа с русскими пояснениями проверена отдельно во встроенном Python FEMM; этот снимок не является новым запуском текущей редакции. Открыть в полном размере

Расчётные исследования по теме

Изучите физическую постановку, выполните программу во встроенном Python FEMM 5 и проверьте результаты по сохранённым моделям и независимым критериям.

Все темы учебного практикума

Окна FEMM 5

Снимки интерфейса FEMM показывают геометрию, настройки и работу с результатами в FEMM 5. Откройте изображение, чтобы рассмотреть подписи и значения в полном размере.

Машина с локальным ослаблением одного постоянного магнита до 60% исходной остаточной индукции: геометрия выбранного положения ротора 15°. Подписи скрыты; размеры магнитов сохранены.
Машина с локальным ослаблением одного постоянного магнита до 60% исходной остаточной индукции: геометрия выбранного положения ротора 15°. Подписи скрыты; размеры магнитов сохранены. Открыть в полном размере
Локальное ослабление одного магнита до 60% исходной остаточной индукции; положение ротора 15°. Статическая карта |B|, Т, и линии поля. Влияние дефекта оценивают сравнением потокосцеплений с исправной моделью при тех же условиях.
Локальное ослабление одного магнита до 60% исходной остаточной индукции; положение ротора 15°. Статическая карта |B|, Т, и линии поля. Влияние дефекта оценивают сравнением потокосцеплений с исправной моделью при тех же условиях. Открыть в полном размере

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

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

Критерий выполнения: Каждая оценка имеет источник, предпосылки и единицу; возмущения и цепь проверены; лишние значащие цифры исключены.
Первоисточники и дополнительное чтение по занятию
  1. [33] — 10–15, 39–52, 90–94
  2. [81] — Раздел Packages · Локальная копия
  3. [5]
Пояснения и типичные ошибки
Равна ли сеточная разность измерительной неопределённости?

Нет. Это разные источники; их объединение требует выбранной методики и предпосылок.

Можно ли число витков возмутить на 0,1 витка?

Для обычной дискретной обмотки нет; число витков анализируется допустимыми целыми конструктивными вариантами.

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

  • Фиксированные токи при изменённом сопротивлении не описывают чувствительность машины с питанием от источника напряжения.
  • Нельзя заявлять 95 \,\% доверие без вероятностной модели и соответствующего коэффициента охвата.