BLDC

Трехобмоточный бесщеточный двигатель постоянного тока с распределением трапеций

  • Библиотека:
  • Simscape/Электрический/Электромеханический/Постоянный магнит

  • BLDC block

Описание

Блок BLDC моделирует синхронную машину с постоянными магнитами с трехфазным статором с соединением обмоток в «звезду». Блок имеет четыре опции для определения распределения потока постоянных магнитов как функции от угла ротора. Две опции допускают простую параметризацию путем принятия идеальной трапеции для коэффициента противо-ЭДС. Для простой параметризации вы задаете редактирование потока или коэффициента противо-ЭДС, вызванную ротором. Другие две опции дают более точные результаты с помощью табличных данных, которые вы задаете. Для более точных результатов вы задаете либо частную производную редактирования, либо измеренный коэффициент противо-ЭДС, постоянная для заданной скорости ротора.

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

Моторная Конструкция

Этот рисунок показывает конструкцию мотора с одной парой полюсов на роторе.

Для соглашения о осях на предыдущем рисунке a-фазы и потоки постоянных магнитов выравниваются, когда угол θr ротора равен нулю. Блок поддерживает второе определение оси ротора. Для второго определения угол ротора является углом между магнитной осью a-фазы и осью q ротора.

Трапециевидная скорость изменения потока

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

Коэффициент противо-ЭДС - скорость изменения потока, заданная как

dΦdt=Φθdθdt=Φθω,

где:

  • Φ - редактирование потока постоянных магнитов.

  • θ - угол ротора.

  • ω - механическая скорость вращения.

Высота h трапеций скорости изменения профиля потока получают из пикового потока постоянных магнитов.

Интеграция Φθ по области значений 0,

Φmax=h2(θF+θW),

где:

  • Φmax - редактирование потока постоянных магнитов.

  • h - скорость изменения высоты профиля потока.

  • θF - угол области значений ротора, над которым коэффициент противо-ЭДС, который индуцирует поток постоянных магнитов в статоре, является постоянным.

  • θW - область значений углов ротора, в котором коэффициент противо-ЭДС увеличивается или уменьшается линейно, когда ротор движется с постоянной скоростью.

Переставляя предшествующее уравнение,

h=2Φmax/(θF+θW).

Электрические определяющие уравнения

Напряжения на обмотках статора определяются

[vavbvc]=[Rs000Rs000Rs][iaibic]+[dψadtdψbdtdψcdt],

где:

  • va, vb и vc являются внешними напряжениями, приложенными к трем электрическим соединениям мотора.

  • Rs - эквивалентное сопротивление каждой обмотки статора.

  • ia, ib и ic являются токами, текущими в обмотках статора.

  • dψadt,dψbdt, и dψcdt

    - скорости изменения магнитного потока в каждой обмотке статора.

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

[ψaψbψc]=[LaaLabLacLbaLbbLbcLcaLcbLcc][iaibic]+[ψamψbmψcm],

где:

  • ψa, ψb и ψc являются суммарными потоками, связывающими каждую обмотку статора.

  • Laa, Lbb и Lcc являются самоиндуктивностями обмоток статора.

  • Lab, Lac, Lba и т.д. являются взаимными индуктивностями обмоток статора.

  • ψam, ψbm и ψcm являются потоками постоянных магнитов, связывающими обмотки статора.

Индуктивности в обмотках статора являются функциями угла ротора, заданными как

Laa=Ls+Lmcos(2θr),

Lbb=Ls+Lmcos(2(θr2π/3)),

Lcc=Ls+Lmcos(2(θr+2π/3)),

Lab=Lba=MsLmcos(2(θr+π/6)),

Lbc=Lcb=MsLmcos(2(θr+π/62π/3)),

и

Lca=Lac=MsLmcos(2(θr+π/6+2π/3)),


где:

  • Ls является самоиндуктивностью статора на фазу - средняя самоиндуктивность каждой из обмоток статора.

  • Lm - индуктивность статора - колебание самоиндуктивности и взаимной индуктивности с изменением угла ротора.

  • Ms является взаимной индуктивностью статора - средняя взаимная индуктивность между обмотками статора.

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

Упрощенные уравнения

Определяющие уравнения напряжения и крутящего момента для блока

[vdvqv0]=P([vavbvc]Nω[ψamθrψbmθrψcmθr]),

vd=Rsid+LddiddtNωiqLq,

vq=Rsiq+Lqdiqdt+NωidLd,

v0=Rsi0+L0di0dt,

и

T=32N(iqidLdidiqLq)+[iaibic][ψamθrψbmθrψcmθr],

где:

  • vd, vq и v0 являются d осью, q осью и напряжениями нулевой последовательности.

  • P - трансформация Парка, заданная как

    P=2/3[cosθecos(θe2π/3)cos(θe+2π/3)sinθesin(θe2π/3)sin(θe+2π/3)0.50.50.5]

  • N - количество пар полюсов постоянных магнитов ротора.

  • ω - механическая скорость вращения ротора.

  • ψamθr,ψbmθr, и ψcmθr

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

  • id, iq и i0 являются d -осью, q -осью и токами нулевой последовательности, заданными как

    [idiqi0]=P[iaibic].

  • Ld = Ls + Ms + 3/2 <reservedrangesplaceholder2>. Ld - индуктивность d оси статора.

  • Lq = Ls + <reservedrangesplaceholder3> − 3/2 <reservedrangesplaceholder2>. Lq - индуктивность q оси статора.

  • L0 = Ls – 2 Ms. L0 - индуктивность нулевой последовательности статора.

  • T - крутящий момент ротора. Крутящий момент течет от корпуса мотора (блокирующий физический порт C) к ротору мотора (блокирующий физический порт R).

Вычисление потерь в железе

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

Термин, представляющий основной путь намагниченности, зависит от индуцированного напряжения статора RMS, Vmrms:

POC(Vmrms)=ahkVmrms+ajk2Vmrms2+aexk1.5Vmrms1.5

Это доминирующий термин во время работы без нагрузки. k является коэффициентом противо-ЭДС, постоянная, относящимся к RMS В на Гц. Он определяется как k=Vmrms/f, где f является электрической частотой. Первый член на правой стороне является потерей магнитного гистерезиса, второй - потеря токов Фуко, третий - избыточные потери. Три коэффициента, появляющиеся на числителях, получают из значений, которые вы обеспечиваете для гистерезиса без обратной связи, вихря и избыточных потерь.

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

PSC(Vdrms*)=bhkVdrms*+bjk2Vdrms*2+bexk1.5Vdrms*1.5

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

Тепловые порты

Блок имеет четыре дополнительных тепловых порта, по одному для каждой из трех обмоток и по одному для ротора. Эти порты по умолчанию скрыты. Чтобы открыть тепловые порты, щелкните правой кнопкой мыши блок в модели, выберите Simscape > Block choices, а затем выберите требуемый вариант блока с тепловыми портами: Composite three-phase ports | Show thermal port или Expanded three-phase ports | Show thermal port. Это действие отображает тепловые порты на значке блока и отображает параметры Temperature Dependence и Thermal Port. Эти параметры описаны далее на этой странице с описанием.

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

Переменные

Используйте настройки Variables, чтобы задать приоритет и начальные целевые значения для основных переменных перед симуляцией. Для получения дополнительной информации смотрите Задать приоритет и Начальный целевой объект для основных переменных.

Порты

Сохранение

расширить все

Расширяемый трехфазный порт.

Электрический порт сопоставлен с нейтральной фазой.

Зависимости

Чтобы включить этот порт, установите Winding Type равным Wye-wound и Zero sequence к Include.

Механический вращательный порт сопоставлен с ротором мотора.

Механический вращательный порт сопоставлен с корпусом мотора.

Тепловой порт, сопоставленный с обмоткой А. Для получения дополнительной информации смотрите Тепловые порты.

Тепловой порт сопоставлен с обмоткой B. Для получения дополнительной информации смотрите Тепловые порты.

Тепловой порт сопоставлен с обмоткой C. Для получения дополнительной информации смотрите Тепловые порты.

Тепловой порт сопоставлен с ротором. Для получения дополнительной информации см. Раздел «Тепловые порты»

Параметры

расширить все

Ротор

Выберите строение обмотки:

  • Wye-wound - Обмотки - рана Уая.

  • Delta-wound - Обмотки размотаны. a -фаза соединяется между портами a и b, b фаза между портами b и c и c фаза между портами c и a.

Параметризация для определения распределения потока постоянных магнитов как функции от угла ротора. Выбирать:

  • Perfect trapezoid - specify maximum flux linkage для определения максимального редактирования потока для постоянного магнита и угла ротора, где коэффициент противо-ЭДС является постоянным. Блок принимает идеальную трапецию для коэффициента противо-ЭДС. Это значение по умолчанию.

  • Perfect trapezoid - specify maximum rotor-induced back emf для определения максимальной индуцированной ротором коэффициента противо-ЭДС и соответствующей скорости ротора. Блок принимает идеальную трапецию для коэффициента противо-ЭДС.

  • Tabulated - specify flux partial derivative with respect to rotor angle для задания значений для частного производного редактирования потока и соответствующих углов ротора.

  • Tabulated - specify rotor-induced back emf as a function of rotor angle для определения измеренного коэффициента противо-ЭДС, постоянная и соответствующей скорости и углов вращения ротора.

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

Зависимости

Чтобы включить этот параметр, установите Back EMF profile равным Perfect trapezoid - specify maximum flux linkage.

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

Зависимости

Чтобы включить этот параметр, установите Back EMF profile равным Perfect trapezoid - specify maximum flux linkage или Perfect trapezoid - specify maximum rotor-induced back emf.

Пик индуцированного ротором коэффициента противо-ЭДС в обмотки статора.

Зависимости

Чтобы включить этот параметр, установите Back EMF profile равным Perfect trapezoid - specify maximum rotor-induced back emf.

Вектор значений индуцированного ротором заднего ЭДС как функции от угла ротора. Первое и последнее значения должны быть одинаковыми и обычно оба равны нулю. Для получения дополнительной информации смотрите параметр Corresponding rotor angles. Первое и последнее значения одинаковы, потому что поток цикличен с периодом 2π/N, где N - количество пар полюсов постоянных магнитов.

Зависимости

Чтобы включить этот параметр, установите Back EMF profile равным Tabulated - specify rotor-induced back emf as a function of rotor angle.

Вектор значений для частного производного редактирования потока (где редактирование потока - это количество витков обмотки) относительно угла ротора. Первое и последнее значения должны быть одинаковыми и обычно оба равны нулю. Для получения дополнительной информации смотрите параметр Corresponding rotor angles. Первое и последнее значения одинаковы, потому что поток цикличен с периодом 2π/N, где N - количество пар полюсов постоянных магнитов.

Вектор углов ротора, где задан поток редактирования частной производной или вызванной ротором коэффициента противо-ЭДС. Угол ротора определяется как угол между магнитной осью a-фазы и осью d-фазы. То есть, когда угол равен нулю, магнитные поля, обусловленные ротором и a-фазовой обмоткой, выравниваются. Это определение используется независимо от настройки блока для определения угла ротора. Первое значение равняется нулю, а последнее значение равняется2π/N, где N - количество пар полюсов постоянных магнитов.

Зависимости

Чтобы включить этот параметр, установите Back EMF profile равным Tabulated - specify flux partial derivative with respect to rotor angle или Tabulated - specify rotor-induced back emf as a function of rotor angle.

Задайте скорость ротора, соответствующую максимальной индуцированной ротором коэффициенту противо-ЭДС.

Зависимости

Чтобы включить этот параметр, установите Back EMF profile равным Perfect trapezoid - specify maximum rotor-induced back emf или Tabulated - specify rotor-induced back emf as a function of rotor angle.

Количество пар полюсов постоянных магнитов на роторе.

Контрольная точка для измерения угла ротора. Значение по умолчанию Angle between the a-phase magnetic axis and the d-axis. Это определение показано на рисунке Конструкции. Когда вы выбираете это значение, потоки ротора и a-фазы выравниваются, когда угол ротора равен нулю.

Другое значение, которое можно выбрать для этого параметра Angle between the a-phase magnetic axis and the q-axis. Когда вы выбираете это значение, ток a-фазы генерирует максимальный крутящий момент, когда угол ротора равен нулю.

Статор

Выберите точность моделирования:

  • Constant Ld and Lq - Ld и Lq значения являются постоянными и определяются их соответствующими параметрами.

  • Tabulated Ld and Lq - Ld и Lq значения вычисляются в оперативном режиме из интерполяционных таблиц токов DQ следующим образом:

    Ld=f1(id,iq)

    Ld=f2(id,iq)

Зависимости

Чтобы включить этот параметр, установите Stator parameterization равным Specify Ld, Lq, and L0.

Выберите Specify Ld, Lq, and L0 или Specify Ls, Lm, and Ms.

Индуктивность оси D.

Зависимости

Чтобы включить этот параметр, установите Stator parameterization равным Specify Ld, Lq, and L0 и Modeling fidelity к Constant Ld and Lq.

Q-составляющая индуктивность.

Зависимости

Чтобы включить этот параметр, установите Stator parameterization равным Specify Ld, Lq, and L0 и Modeling fidelity к Constant Ld and Lq.

Вектор тока прямой оси, iD.

Зависимости

Чтобы включить этот параметр, установите Stator parameterization равным Specify Ld, Lq, and L0 и Modeling fidelity к Tabulated Ld and Lq.

Вектор тока с квадратурной осью, iQ.

Зависимости

Чтобы включить этот параметр, установите Stator parameterization равным Specify Ld, Lq, and L0 и Modeling fidelity к Tabulated Ld and Lq.

Матрица Ld.

Зависимости

Чтобы включить этот параметр, установите Stator parameterization равным Specify Ld, Lq, and L0 и Modeling fidelity к Tabulated Ld and Lq.

Матрица Lq.

Зависимости

Чтобы включить этот параметр, установите Stator parameterization равным Specify Ld, Lq, and L0 и Modeling fidelity к Tabulated Ld and Lq.

Индуктивность нулевой последовательности.

Зависимости

Чтобы включить этот параметр, либо:

  • Установите Winding Type значение Wye-wound, Zero sequence к Include, и Stator parameterization к Specify Ld, Lq, and L0.

  • Установите Winding Type значение Delta-wound и Stator parameterization к Specify Ld, Lq, and L0.

Средняя индуктивность каждой из трех обмоток статора.

Зависимости

Чтобы включить этот параметр, установите Stator parameterization равным Specify Ls, Lm, and Ms.

Колебания самоиндуктивности и взаимной индуктивности обмоток статора с углом ротора.

Зависимости

Чтобы включить этот параметр, установите Stator parameterization равным Specify Ls, Lm, and Ms.

Средняя взаимная индуктивность между обмотками статора.

Зависимости

Чтобы включить этот параметр, установите Stator parameterization равным Specify Ls, Lm, and Ms.

Сопротивление каждой из обмоток статора.

Опция включения или исключения членов нулевой последовательности.

  • Include - Включите условия нулевой последовательности. Чтобы расставить приоритеты для точности модели, используйте эту настройку по умолчанию. Использование этой опции:

  • Exclude - Исключить условия нулевой последовательности. Чтобы определить приоритеты скорости симуляции для симуляции рабочего стола или развертывания приложений, выберите эту опцию.

Зависимости

Этот параметр видим только, когда вы устанавливаете параметр Winding Type равным Wye-wound.

Потери в железе

Задайте вычислительную модель потерь в железе.

Вектор-строка, длиной 3, потерь в железе из-за гистерезиса, Эдди и лишних потерь, соответственно, на частоте, заданной Electrical frequency at which losses determined.

Зависимости

Чтобы включить этот параметр, установите Iron-loss равным Empirical.

Вектор-строка, длиной 3, потерь в железе с короткой схемой из-за гистерезиса, Эдди и лишних потерь, соответственно, на частоте, заданной Electrical frequency at which losses determined.

Зависимости

Чтобы включить этот параметр, установите Iron-loss равным Empirical.

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

Зависимости

Чтобы включить этот параметр, установите Iron-loss равным Empirical.

Результат короткой схемы ток фазы RMS при измерении потерь короткой схемы.

Зависимости

Чтобы включить этот параметр, установите Iron-loss равным Empirical.

Механический

Инерция ротора, прикрепленного к механическому переводному порту R. Значение может быть нулем.

Вращательное демпфирование.

Температурная зависимость

Эти параметры появляются только для блоков с открытыми тепловыми портами. Для получения дополнительной информации см. Раздел «Тепловые порты»

Температура, для которой приведены параметры двигателя.

Коэффициент α в уравнении, относящем сопротивление температуре, как описано в Тепловой модели для блоков привода. Значение по умолчанию для меди.

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

Тепловой порт

Эти параметры появляются только для блоков с открытыми тепловыми портами. Для получения дополнительной информации см. Раздел «Тепловые порты»

Значение тепловой массы для обмоток A, B и C. Тепловая масса является энергией, необходимой для повышения температуры на одну степень.

Тепловая масса ротора, то есть энергия, необходимая для повышения температуры ротора на одну степень.

Процент потерь в железе основного пути потока сопоставлен с магнитным путем через ротор. Он определяет, какую часть нагрева потерь в железе относят к тепловому порту HR ротора и какую часть относят к трем обмоткам тепловых портов HA, HB и HC.

Зависимости

Чтобы включить этот параметр, установите Iron-loss равным Empirical.

Процент потерь в железе пути поперечного зуба, сопоставленных с магнитным путем через ротор. Он определяет, какую часть нагрева потерь в железе относят к тепловому порту HR ротора и какую часть относят к трем обмоткам тепловых портов HA, HB и HC.

Зависимости

Чтобы включить этот параметр, установите Iron-loss равным Empirical.

Примеры моделей

BLDC Position Control with Thermal Model

Управление положением BLDC с помощью тепловой модели

Управляйте углом ротора в электроприводе на базе BLDC. BLDC включает тепловую модель и эмпирические потери в железе. Идеальный источник крутящего момента обеспечивает нагрузку. Подсистема Control использует основанную на ПИ структуру каскадного регулирования с тремя циклами управления: внешний цикл управления положения, цикл управления скорости и внутренний цикл управления током. BLDC питается управляемым трехфазным инвертором. Сигналы управления ключами для инвертора получаются из сигналов зала. В симуляции используются ссылки на шаги. Начальная температура обмоток статора и ротора установлена на 25 степени Цельсия. Температура окружающей среды составляет 27 степени Цельсия. Подсистема Возможностей содержит возможности, которые позволяют вам видеть результаты симуляции.

Ссылки

[1] Кундур, П. Устойчивость системы Степени и Управление. Нью-Йорк, Нью-Йорк: McGraw Hill, 1993.

[2] Андерсон, П. М. Анализ неисправных степеней. Hoboken, NJ: Wiley-IEEE Press, 1995.

[3] Меллор, П.Х., Р. Вробель и Д. Холлидей. «Вычислительно эффективная модель потери в железе для бесщеточных машин переменного тока, которая удовлетворяет номинальным потокам и полевым ослабленным операциям». Конференция по электрическим машинам и приводам IEEE. Май 2009.

Расширенные возможности

Генерация кода C/C + +
Сгенерируйте код C и C++ с помощью Coder™ Simulink ®

.
Введенный в R2013b