Flux-Based PMSM

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

  • Библиотека:
  • Блоксет силовых агрегатов/Двигатели/Электродвигатели и инверторы

  • Flux-Based PMSM block

Описание

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

Модели электродвигателя на основе потока учитывают магнитное насыщение и потери в железе. Чтобы вычислить магнитное насыщение и потерю в железе, блок Flux-Based PMSM использует обратные редактирования потока. Чтобы получить параметры блоков, можно использовать конечноэлементный анализ (FEA) или измерить напряжения фазы с помощью динамометра.

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

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

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

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

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

ВычислениеУравнение
q- и d-составляющие напряженияvd=dψddt+Rsidωeψqvq=dψqdt+Rsiq+ωeψd
ток q- и d-осейid=f(ψd,ψq)iq=g(ψd,ψq)
Электромеханический крутящий моментTe=1.5P[ψdiqψqid]

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

ωm

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

ωe

Электрическая скорость ротора

Θda

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

Rs, Rr

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

iq, id

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

vq, vd

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

Ψq, Ψd

q- и d-составляющие потоки магнитов, соответственно

P

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

Te

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

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

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

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

ωe=Pωmdθedt= ωe

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

Кларк

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

xα= 23xa 13xb 13xcxβ= 32xb 32xc

Парк

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

xd= xαcosθe+ xβsinθexq= xαsinθe+ xβcosθe

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

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

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

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

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

xα= xdcosθe xqsinθexβ= xdsinθe+ xqcosθe

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

ωm

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

P

Пары шестов двигателей

ωe

Электрическая скорость ротора

Θe

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

x

Ток фазы или напряжение

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

Скорость вращения ротора задается:

ddtωm=1J(TeTfFωmTm)dθmdt=ωm

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

J

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

F

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

θm

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

Tm

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

Te

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

Tf

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

ωm

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

Учет степени

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

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

PwrInfo

PwrTrnsfrd - Степень между блоками

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

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

PwrMtr

Механическая степень

Pmot

Pmot= ωmTe
PwrBus

Электрическая степень

Pbus

Pbus= vania+ vbnib+vcnic

PwrNotTrnsfrd - Степень через контур блока, но не переданный

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

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

PwrElecLoss

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

Pelec

Pelec= 32(Rsisd2+Rsisq2)
PwrMechLoss

Механические потери степени

Pmech

Когда Port Configuration установлено на Torque:

Pmech= (ωm2F+ |ωm|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-осей статора

van, vbn, vcn

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

ωm

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

F

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

Te

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

Tf

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

Оптимизация памяти интерполяционной таблицы

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

Чтобы включить текущие вычисления, подходящие для целей генерации кода, которые ограничивают память, выберите Enable memory optimized 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

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

ia

A

IbStator

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

ib

A

IcStator

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

ic

A

IdSync

Ток прямой оси

id

A

IqSync

Квадратурный ток по оси

iq

A

VdSync

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

vd

V

VqSync

Квадратурная ось напряжение

vq

V

MtrSpd

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

ωm

рад/с

MtrPos

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

θm

рад

MtrTrq

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

Te

Н· м

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 configuration.

Угловая скорость двигателя, ωmtr, в рад/с.

Зависимости

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

Параметры

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

Опции блока

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

Зависимости

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

Интегрирование шага расчета для дискретной симуляции, в с.

Зависимости

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

В этой таблице представлены строения портов.

Строение портаСоздает Input portСоздает выходной порт

Torque

LdTrq

MtrSpd

Speed

Spd

MtrTrq

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

d-составляющая потока, Ψd, точек останова, в Wb.

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

Зависимости

Чтобы создать этот параметр, выберите Enable memory optimized 2D LUT.

q-составляющая потока, Ψq, точек останова, в Wb.

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

Зависимости

Чтобы создать этот параметр, выберите Enable memory optimized 2D LUT.

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

Если вы задаете Enable memory optimized 2D LUT, блок преобразует данные в одинарную точность.

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

Если вы задаете Enable memory optimized 2D LUT, блок преобразует данные в одинарную точность.

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

Зависимости

Чтобы создать этот параметр, выберите Enable memory optimized 2D LUT и Specify Extrapolation.

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

Зависимости

Чтобы создать этот параметр, выберите Enable memory optimized 2D LUT и Specify Extrapolation.

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

Зависимости

Чтобы создать этот параметр, выберите Enable memory optimized 2D LUT и Specify Extrapolation.

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

Зависимости

Чтобы создать этот параметр, выберите Enable memory optimized 2D LUT и Specify Extrapolation.

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

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

Начальный поток, Ψq0 и Ψd0 оси D и Q в Wb.

Начальное угловое положение ротора, θm0, в рад.

Начальная скорость вращения ротора, ωm0, в рад/с.

Зависимости

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

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

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

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

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

Зависимости

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

Ссылки

[1] Hu, Dakai, Yazan Alsmadi, and Longya Xu. «Высокоточное нелинейное моделирование IPM на основе измеренных редактирований потока обмотки статора». IEEE® Сделки по отраслевым приложениям, том 51, № 4, июль/август 2015 года.

[2] Чен, Сяо, Джиэбин Ван, Bhaskar сенатор, Панайотис Лазари, Tianfu Солнце. «Высокочастотная и в вычислительном отношении эффективная модель для внутренних машин постоянного магнита, рассматривая магнитное насыщение, пространственную гармонику и эффект потери в железе». Сделки IEEE по промышленной электронике, том 62, № 7, июль 2015 года.

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

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

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

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