ee_generateIdealPMSMfluxData

Сгенерируйте сведенные в таблицу данные о потокосцеплении для идеального PMSM

Синтаксис

[F,T,dFdA,dFdB,dFdC,dFdX] = ee_generateIdealPMSMfluxData(PM,Ld,Lq,L0,A,B,C,X)
F = ee_generateIdealPMSMfluxData(PM,Ld,Lq,L0,A,B,C,X)
[F,T,dFdA,dFdB,dFdC,dFdX] = ee_generateIdealPMSMfluxData(PM,Ld,Lq,L0,D,Q,X)
F = ee_generateIdealPMSMfluxData(PM,Ld,Lq,L0,D,Q,X)

Описание

пример

[F,T,dFdA,dFdB,dFdC,dFdX] = ee_generateIdealPMSMfluxData(PM,Ld,Lq,L0,A,B,C,X) генерирует 4-D данные о потокосцеплении, включая крутящий момент и частные производные, для идеального постоянного магнита синхронного двигателя (PMSM).

Используйте эту функцию, чтобы создать тестовые данные для блока FEM-Parameterized PMSM, или в целях валидации или настроить модель, прежде чем фактические данные о потокосцеплении будут доступны.

пример

F = ee_generateIdealPMSMfluxData(PM,Ld,Lq,L0,A,B,C,X) генерирует 4-D матрицу потокосцепления F для идеального PMSM.

пример

[F,T,dFdA,dFdB,dFdC,dFdX] = ee_generateIdealPMSMfluxData(PM,Ld,Lq,L0,D,Q,X) генерирует 3-D данные о потокосцеплении, включая крутящий момент и частные производные, для идеального PMSM.

пример

F = ee_generateIdealPMSMfluxData(PM,Ld,Lq,L0,D,Q,X) генерирует 3-D матрицу потокосцепления F для идеального PMSM.

Примеры

Сгенерируйте 4-D данные о потокосцеплении

Задайте моторные параметры.

PM = 0.1;     % Permanent magnet flux
N = 6;        % Number of pole pairs
Ld = 0.0002;  % D-axis inductance
Lq = 0.0002;  % Q-axis inductance
L0 = 0.00018; % Zero-sequence inductance
Rs = 0.013;   % Stator resistance

Задайте фазу текущие векторы.

iA = linspace(-250,250,5);
iB = iA;
iC = iA;

Задайте угловой вектор ротора на основе количества пар полюса.

X = pi/180*linspace(0,360/N,180/N+1);

Плоские частные производные потокосцепления и крутящий момент с точки зрения A-, B-, C-токов и угла ротора

[F,T,dFdA,dFdB,dFdC,dFdX] = ee_generateIdealPMSMfluxData(PM,Ld,Lq,L0,iA,iB,iC,X);

Функция возвращает 4-D матрицу F потокосцепления, матрицу T крутящего момента 4-D и четыре 4-D матрицы для частных производных потокосцепления. Четыре матрицы частной производной соответствуют трем токам фазы и углу ротора, соответственно. Матричные размерности соответствуют трем токам фазы и углу ротора.

Сгенерируйте 4-D матрицу потокосцепления F

Задайте моторные параметры.

PM = 0.1;     % Permanent magnet flux
N = 6;        % Number of pole pairs
Ld = 0.0002;  % D-axis inductance
Lq = 0.0002;  % Q-axis inductance
L0 = 0.00018; % Zero-sequence inductance
Rs = 0.013;   % Stator resistance

Задайте фазу текущие векторы.

iA = linspace(-250,250,5);
iB = iA;
iC = iA;

Задайте угловой вектор ротора на основе количества пар полюса.

X = pi/180*linspace(0,360/N,180/N+1);

Плоские частные производные потокосцепления и крутящий момент с точки зрения A-, B-, C-токов и угла ротора

F = ee_generateIdealPMSMfluxData(PM,Ld,Lq,L0,iA,iB,iC,X);

Функция возвращает 4-D матрицу F потокосцепления. Матричные размерности соответствуют трем токам фазы и углу ротора.

Сгенерируйте 3-D данные о потокосцеплении

Задайте моторные параметры.

PM = 0.1;     % Permanent magnet flux
N = 6;        % Number of pole pairs
Ld = 0.0002;  % D-axis inductance
Lq = 0.0002;  % Q-axis inductance
L0 = 0.00018; % Zero-sequence inductance
Rs = 0.013;   % Stator resistance

Задайте d-ось и q-ось текущие векторы.

iD = linspace(-250,250,5);
iQ = iD;

Задайте угловой вектор ротора на основе количества пар полюса.

X = pi/180*linspace(0,360/N,180/N+1);

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

[F,T,dFdA,dFdB,dFdC,dFdX] = ee_generateIdealPMSMfluxData(PM,Ld,Lq,L0,iD,iQ,X);

Функция возвращает 3-D матрицу F потокосцепления, 3-D матрицу T крутящего момента и четыре 3-D матрицы для частных производных потокосцепления. Четыре матрицы частной производной соответствуют трем токам фазы и углу ротора, соответственно. Матричные размерности соответствуют d-оси и токам q-оси и углу ротора.

Сгенерируйте 3-D матрицу потокосцепления F

Задайте моторные параметры.

PM = 0.1;     % Permanent magnet flux
N = 6;        % Number of pole pairs
Ld = 0.0002;  % D-axis inductance
Lq = 0.0002;  % Q-axis inductance
L0 = 0.00018; % Zero-sequence inductance
Rs = 0.013;   % Stator resistance

Задайте d-ось и q-ось текущие векторы.

iD = linspace(-250,250,5);
iQ = iD;

Задайте угловой вектор ротора на основе количества пар полюса.

X = pi/180*linspace(0,360/N,180/N+1);

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

F = ee_generateIdealPMSMfluxData(PM,Ld,Lq,L0,iD,iQ,X);

Функция возвращает 3-D матрицу потокосцепления F. Матричные размерности соответствуют d-оси и токам q-оси и углу ротора.

Входные параметры

свернуть все

Пиковое потокосцепление постоянного магнита, в weber-поворотах, заданных как скаляр.

Типы данных: double

D- индуктивность оси, в henries, заданном как скаляр.

Типы данных: double

Q- индуктивность оси, в henries, заданном как скаляр.

Типы данных: double

Индуктивность нулевой последовательности, в henries, заданном как скаляр.

Типы данных: double

Текущая A-фаза, в амперах, заданных как вектор. Вектор должен монотонно увеличиваться, и двухсторонний (содержите и положительные и отрицательные величины). Лучшая практика состоит в том, чтобы включать нулевой ток как одну из точек. Используйте этот входной параметр, чтобы сгенерировать 4-D данные о потокосцеплении.

Типы данных: double

Текущая B-фаза, в амперах, заданных как вектор. Вектор должен монотонно увеличиваться, и двухсторонний (содержите и положительные и отрицательные величины). Лучшая практика состоит в том, чтобы включать нулевой ток как одну из точек. Используйте этот входной параметр, чтобы сгенерировать 4-D данные о потокосцеплении.

Типы данных: double

Текущая C-фаза, в амперах, заданных как вектор. Вектор должен монотонно увеличиваться, и двухсторонний (содержите и положительные и отрицательные величины). Лучшая практика состоит в том, чтобы включать нулевой ток как одну из точек. Используйте этот входной параметр, чтобы сгенерировать 4-D данные о потокосцеплении.

Типы данных: double

Текущая D-ось, в амперах, заданных как вектор. Вектор должен монотонно увеличиваться, и двухсторонний (содержите и положительные и отрицательные величины). Лучшая практика состоит в том, чтобы включать нулевой ток как одну из точек. Используйте этот входной параметр, чтобы сгенерировать 3-D данные о потокосцеплении.

Типы данных: double

Текущая Q-ось, в амперах, заданных как вектор. Вектор должен монотонно увеличиваться, и двухсторонний (содержите и положительные и отрицательные величины). Лучшая практика состоит в том, чтобы включать нулевой ток как одну из точек. Используйте этот входной параметр, чтобы сгенерировать 3-D данные о потокосцеплении.

Типы данных: double

Угол ротора, в радианах, заданных как вектор. Значения должны быть в диапазоне от нуля до 2π/N, где N является количеством пар полюса.

Типы данных: double

Выходные аргументы

свернуть все

Потокосцепление, в weber-поворотах, возвратилось как матрица. Матрица может быть четырехмерной или 3D, в зависимости от синтаксиса, используемого, чтобы вызвать функцию. В четырехмерной матрице первые три измерения соответствуют токам фазы, и четвертая размерность соответствует углу ротора. В 3D матрице первые две размерности соответствуют d - оси и q - токи оси, и третья размерность соответствует углу ротора.

Закрутите, в N*m, возвращенном как матрица. Матрица может быть четырехмерной или 3D, в зависимости от синтаксиса, используемого, чтобы вызвать функцию. В четырехмерной матрице первые три измерения соответствуют токам фазы, и четвертая размерность соответствует углу ротора. В 3D матрице первые две размерности соответствуют d - оси и q - токи оси, и третья размерность соответствует углу ротора.

Частная производная потокосцепления относительно текущей A-фазы, возвращенной как матрица. Матрица может быть четырехмерной или 3D, в зависимости от синтаксиса, используемого, чтобы вызвать функцию. В четырехмерной матрице первые три измерения соответствуют токам фазы, и четвертая размерность соответствует углу ротора. В 3D матрице первые две размерности соответствуют d - оси и q - токи оси, и третья размерность соответствует углу ротора.

Частная производная потокосцепления относительно текущей B-фазы, возвращенной как матрица. Матрица может быть четырехмерной или 3D, в зависимости от синтаксиса, используемого, чтобы вызвать функцию. В четырехмерной матрице первые три измерения соответствуют токам фазы, и четвертая размерность соответствует углу ротора. В 3D матрице первые две размерности соответствуют d - оси и q - токи оси, и третья размерность соответствует углу ротора.

Частная производная потокосцепления относительно текущей C-фазы, возвращенной как матрица. Матрица может быть четырехмерной или 3D, в зависимости от синтаксиса, используемого, чтобы вызвать функцию. В четырехмерной матрице первые три измерения соответствуют токам фазы, и четвертая размерность соответствует углу ротора. В 3D матрице первые две размерности соответствуют d - оси и q - токи оси, и третья размерность соответствует углу ротора.

Частная производная потокосцепления относительно угла ротора, возвращенного как матрица. Матрица может быть четырехмерной или 3D, в зависимости от синтаксиса, используемого, чтобы вызвать функцию. В четырехмерной матрице первые три измерения соответствуют токам фазы, и четвертая размерность соответствует углу ротора. В 3D матрице первые две размерности соответствуют d - оси и q - токи оси, и третья размерность соответствует углу ротора.

Алгоритмы

Поток, соединяющий каждую обмотку, имеет вклады от постоянного магнита плюс эти три обмотки. Поэтому общий поток дан [1]:

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

Laa=Ls+Lmпотому что(2θr)Lbb=Ls+Lmпотому что(2(θr2π/3))Lcc=Ls+Lmпотому что(2(θr+2π/3))Lab=Lba=MsLmпотому что(θr+π/6)Lbc=Lcb=MsLmпотому что(θr+π/62π/3)Lca=Lac=MsLmпотому что(θr+π/6+2π/3)ψam=ψmпотому чтоθeψbm=ψmпотому что(θe2π/3)ψbm=ψmпотому что(θe+2π/3)

Здесь, Θ e является электрическим углом, который связан с углом ротора Θ r Θ e = N ·Θ r. Функция принимает, что поток постоянного магнита, соединяющий обмотку A-фазы, в максимуме для Θ e = 0.

Функциональный вывод F соответствует ψ сведенный в таблицу как функция текущей A-фазы, текущая B-фаза, текущая C-фаза, и угол ротора.

Ls, Lm и Ms связаны с входными параметрами Ld, Lq и L0:

Ls=L03+Ld3+Lq3Ms=Ld6L03+Lq6Lm=Ld3Lq3

Ссылки

[1] Андерсон, пополудни Анализ Неработающих Энергосистем. 1-й Выпуск. Нажатие Wiley-IEEE, июль 1995, p.187.

Введенный в R2017a