exponenta event banner

PMSM на основе потока

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

  • Библиотека:
  • Блок силового агрегата/движитель/электродвигатели и инверторы

  • Flux-Based PMSM block

Описание

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

Модели двигателей на основе потока учитывают магнитное насыщение и потери железа. Для расчета магнитного насыщения и потерь железа блок PMSM на основе потока использует обратную связь потока. Для получения параметров блока можно использовать конечноэлементный анализ (МКЭ) или измерять фазовые напряжения с помощью динамометра.

По умолчанию блок устанавливает для параметра «Тип моделирования» значение Continuous для использования непрерывного времени выборки во время моделирования. Если требуется создать код для целей с двойной и одинарной точностью с фиксированным шагом, рекомендуется установить для параметра значение Discrete. Затем укажите параметр Sample Time, Ts.

Чтобы включить вычисления потерь мощности, подходящие для целей генерации кода, ограничивающих объем памяти, выберите Включить оптимизированную память 2D LUT.

Трехфазная синусоидальная модель электрической системы

Блок реализует уравнения, которые выражены в неподвижном опорном (dq) кадре ротора. Ось d выравнивается по оси A. Все величины в системе отсчета ротора относятся к статору.

Блок использует эти уравнения.

ВычислениеУравнение
q- и d-осевое напряжениеvd = dü ddt + Rsid
q- и d-осевой токid = f , «) iq = g («d», «q»)
Электромеханический крутящий моментTe = 1,5P [

Уравнения используют эти переменные.

ωm

Механическая частота вращения ротора

ωe

Электрическая частота вращения ротора

Θda

dq электрический угол статора относительно оси а ротора

Rs, Rr

Сопротивление обмоток статора и ротора соответственно

iq, id

q- и d-осевой ток соответственно

vq, vd

q- и d-осевое напряжение соответственно

Λ q, Λ d

магнитный поток по оси q- и d соответственно

P

Количество пар полюсов

Те

Электромагнитный момент

Преобразовывает

Для вычисления напряжений и токов в сбалансированных трехфазных (a, b) величинах, квадратурных двухфазных (α, β) величинах и вращающихся (d, q) опорных кадрах блок использует преобразования Кларка и Парка.

В уравнениях преобразования.

starte =  Pstartmdstartedt = starte

ПреобразоватьОписаниеУравнения

Кларк

Преобразует сбалансированные трехфазные величины (a, b) в сбалансированные двухфазные квадратурные величины (α, β).

 = 23xa  13xb  13xxxβ =  32xb 32xc

Парк

Преобразует сбалансированные двухфазные ортогональные стационарные величины (α, β) в ортогональный вращающийся опорный кадр (d, q).

xd  = xαcosstarte +  xβsinü exq =  xαsinstarte + xβcosstarte

Обратный Кларк

Преобразует сбалансированные двухфазные квадратурные величины (α, β) в сбалансированные трехфазные величины (a, b).

xa  = xaxb =  12xα +  32xβxc =  12xα − 32xβ

Обратный парк

Преобразует ортогональный вращающийся опорный кадр (d, q) в сбалансированные двухфазные ортогональные стационарные величины (α, β).

 = xdcosstarte  xqsinü exβ =  xdsinstarte + xqcosstarte

Преобразования используют эти переменные.

ωm

Механическая частота вращения ротора

P

Пары полюсов двигателя

ωe

Электрическая частота вращения ротора

Θe

Электрический угол ротора

x

Фазный ток или напряжение

Механическая система

Угловая скорость ротора определяется:

ddtü m = 1J (Te Tf Fü m Tm) d

Уравнения используют эти переменные.

J

Комбинированная инерция ротора и нагрузки

F

Комбинированное вязкое трение ротора и нагрузки

θm

Механическое угловое положение ротора

TM

Крутящий момент вала ротора

Те

Электромагнитный момент

Tf

Комбинированный момент трения ротора и нагрузки

ωm

Механическая частота вращения ротора

Учет мощности

Для учета мощности блок реализует эти уравнения.

Сигнал шины ОписаниеПеременнаяУравнения

PwrInfo

PwrTrnsfrd - Мощность, передаваемая между блоками

  • Положительные сигналы указывают на поток в блок

  • Отрицательные сигналы указывают на выход потока из блока

PwrMtr

Механическая мощность

Pmot

Pmot  = startmTe
PwrBus

Электроэнергия

Pbus

Pbus  = vania + vbnib + vcnic

PwrNotTrnsfrd - Мощность, пересекающая границу блока, но не передаваемая

  • Положительные сигналы указывают на вход

  • Отрицательные сигналы указывают на потерю

PwrElecLoss

Потеря резистивной мощности

Pelec

Пелек  = 32 (Rsisd2 + Rsisq2)
PwrMechLoss

Потеря механической мощности

Pmech

Если для параметра Port Configuration установлено значение Torque:

Pmech  = ( startm2F + | startm 'Tf)

Если для параметра Port Configuration установлено значение Speed:

Pmech  =  0

PwrStored - Скорость изменения накопленной энергии

  • Положительные сигналы указывают на увеличение

  • Отрицательные сигналы указывают на снижение

PwrMtrStored

Накопленная мощность двигателя

Pstr

Pstr  = Pbus +  Pmot +   Pelec + Pmech

Уравнения используют эти переменные.

RS

Сопротивление статора

ia, ib, ic

Статорный ток фаз a, b и c

isq, isd

Токи q- и d-оси статора

фургон, vbn, vcn

Напряжение фазы a, b и c статора

ωm

Угловая механическая скорость ротора

F

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

Те

Электромагнитный момент

Tf

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

Оптимизация памяти таблицы подстановки

Данные для Соответствующих токов d-оси, id и Соответствующих токов q-оси, таблиц поиска iq являются функциями потока d-оси и q-оси.

Чтобы включить текущие вычисления, подходящие для целей создания кода, ограничивающих объем памяти, выберите Включить оптимизированную память 2D LUT. Блок использует линейную интерполяцию для оптимизации текущих значений таблицы поиска для генерации кода. В этой таблице приводится сводная информация о реализации оптимизации.

Сценарий использованияВнедрение

Поток d- и q-оси выравнивается со значениями точки останова таблицы поиска.

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

Поток по осям d- и q не выравнивается со значениями точки останова таблицы поиска, но находится в пределах диапазона.

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

Поток по осям d- и q не выравнивается со значениями точки останова таблицы поиска и выходит за пределы диапазона.

Не удается вычислить ток, оптимизированный для памяти. Блок использует экстраполированные данные.

Экстраполяция

Таблицы поиска, оптимизированные для создания кода, не поддерживают экстраполяцию данных, находящихся вне допустимого диапазона. Однако предварительно рассчитанные значения экстраполяции можно включить в таблицу поиска потерь мощности, выбрав команду Указать экстраполяцию (Specify Extrapolation).

Блок использует параметры конечной точки для изменения размера табличных данных.

Ввод данных пользователемЭкстраполяция

Порты

Вход

развернуть все

Входной крутящий момент вала ротора, Tm, в Н· м.

Зависимости

Чтобы создать этот порт, выберите Torque для параметра Port Configuration.

Угловая скорость ротора, λ м, в рад/с.

Зависимости

Чтобы создать этот порт, выберите Speed для параметра Port Configuration.

Напряжение на клеммах статора, Va, Vb и Vc, в В.

Зависимости

Чтобы создать этот порт, выберите Speed или Torque для параметра Port Configuration.

Продукция

развернуть все

Сигнал шины содержит эти блочные вычисления.

Сигнал ОписаниеПеременнаяЕдиницы

IaStator

Ток фазы статора А

ia

A

IbStator

Ток фазы статора B

ib

A

IcStator

Фазный ток статора C

ic

A

IdSync

Постоянный ток оси

я бы

A

IqSync

Ток квадратурной оси

IQ

A

VdSync

Напряжение по прямой оси

vd

V

VqSync

Напряжение квадратурной оси

vq

V

MtrSpd

Угловая механическая скорость ротора

ωm

рад/с

MtrPos

Механическое угловое положение ротора

θm

радиус

MtrTrq

Электромагнитный момент

Те

Н· м

PwrInfo

PwrTrnsfrd

PwrMtr

Механическая мощность

Pmot

W
PwrBus

Электроэнергия

Pbus

W

PwrNotTrnsfrd

PwrElecLoss

Потеря резистивной мощности

Pelec

W
PwrMechLoss

Потеря механической мощности

Pmech

W
PwrStoredPwrMtrStored

Накопленная мощность двигателя

Pstr

W

Фазы a, b, c тока, ia, ib и ic, в А.

Крутящий момент двигателя, Tmtr, в Н· м.

Зависимости

Чтобы создать этот порт, выберите Speed для параметра конфигурации Port.

Угловая скорость мотора, startmtr, в рад/с.

Зависимости

Чтобы создать этот порт, выберите Torque для параметра конфигурации Port.

Параметры

развернуть все

Параметры блока

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

Зависимости

Установка для типа моделирования значения Discrete создает параметр Sample Time, Ts.

Время интеграционной выборки для дискретного моделирования, в с.

Зависимости

Установка для типа моделирования значения Discrete создает параметр Sample Time, Ts.

В этой таблице представлены конфигурации портов.

Конфигурация портаСоздание входного портаСоздание выходного порта

Torque

LdTrq

MtrSpd

Speed

Spd

MtrTrq

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

d-осевой флюс, Λ d, точки останова, в Wb.

Размер хранилища точек останова потока, n1, безразмерный. Блок производит повторную выборку соответствующего тока d-оси, id и соответствующего тока q-оси, iq данных на основе размера хранилища.

Зависимости

Чтобы создать этот параметр, выберите Включить оптимизированную память 2D LUT.

q-осевой поток, Λ q, точки останова, в Wb.

Размер хранилища точек останова потока, n2, безразмерный. Блок производит повторную выборку соответствующего тока d-оси, id и соответствующего тока q-оси, iq данных на основе размера хранилища.

Зависимости

Чтобы создать этот параметр, выберите Включить оптимизированную память 2D LUT.

Массив значений для тока d-оси, id, как функции M d-потоки, N q-fluxes, Λ q, in A. Каждое значение задает ток для конкретной комбинации потока по оси d- и q. Размер массива должен соответствовать размерам, определенным векторами потока.

Если установлен флажок «Включить память, оптимизированную 2D LUT», блок преобразует данные в единую точность.

Массив значений тока по оси q, id, как функция M d-потоки, N q-fluxes, Λ q, in A. Каждое значение задает ток для конкретной комбинации потока по оси d- и q. Размер массива должен соответствовать размерам, определенным векторами потока.

Если установлен флажок «Включить память, оптимизированную 2D LUT», блок преобразует данные в единую точность.

Конечная точка максимальной экстраполяции точки останова потока u1max в Wb.

Зависимости

Чтобы создать этот параметр, выберите Включить оптимизированную память 2D LUT и Указать экстраполяцию.

Минимальная конечная точка экстраполяции точки останова потока, u1min, в Wb.

Зависимости

Чтобы создать этот параметр, выберите Включить оптимизированную память 2D LUT и Указать экстраполяцию.

Конечная точка максимальной экстраполяции точки останова потока u2max в Wb.

Зависимости

Чтобы создать этот параметр, выберите Включить оптимизированную память 2D LUT и Указать экстраполяцию.

Минимальная конечная точка экстраполяции точки останова потока, u2min, в Wb.

Зависимости

Чтобы создать этот параметр, выберите Включить оптимизированную память 2D LUT и Указать экстраполяцию.

Сопротивление фазы статора, Rs, в Ом.

Пары полюсов двигателей, П.

Начальный поток по осям d- и q, Ψq0 и Ψd0, в Wb.

Исходная угловая позиция ротора, startm0, в рад.

Исходная угловая скорость ротора, λ m0, в рад/с.

Зависимости

Чтобы включить этот параметр, выберите Torque параметр конфигурации.

Механические свойства ротора:

  • Инерция, Дж, в кгм ^ 2

  • Вязкое демпфирование, F, в Н· м/( рад/с)

  • Статическое трение, Tf, в Н· м

Зависимости

Чтобы включить этот параметр, выберите Torque параметр конфигурации.

Ссылки

[1] Ху, Дакай, Язан Алсмади и Лунъя Сюй. «Высокоточное нелинейное моделирование IPM на основе измеренной связи потока обмотки статора». Транзакции IEEE ® по отраслевым приложениям, том 51, № 4, июль/август 2015 г.

[2] Чэнь, Сяо, Цзябинь Ван, Бхаскар Сен, Панагиотис Ласари, Тяньфу Сунь. «Модель высокой точности и вычислительной эффективности для внутренних машин с постоянными магнитами, учитывающая магнитное насыщение, пространственную гармонику и эффект потерь железа». IEEE Transactions on Industrial Electronics, том 62, № 7, июль 2015 года.

[3] Оттоссон, Дж., М. Алакула. «Компактная реализация контроллера ослабления поля». Международный симпозиум по силовой электронике, электрическим приводам, автоматизации и движению, июль 2006 года.

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

Создание кода C/C + +
Создайте код C и C++ с помощью Simulink ® Coder™

.
Представлен в R2017b