Моделирование динамики трехфазной асинхронной машины, также известной как индукционная машина, в единицах СИ или пу
Simscape / Электрический / Специализированные Энергосистемы / Электрические Машины

Блоки Asynchronous Machine SI Units и Asynchronous Machine pu Units реализуют трехфазную асинхронную машину (намотанный ротор, клетка с белкой или двойная клетка с белкой), смоделированную в выбираемой опорной системе dq (ротор, статор или синхронный). Обмотки статора и ротора соединены тросом с внутренней нейтральной точкой. Блок работает либо в генераторном, либо в моторном режиме. Режим работы диктуется признаком механического момента:
Если Tm положителен, машина действует как двигатель.
Если Tm отрицательный, машина действует как генератор.
Электрическая часть машины представлена моделью состояния-пространства четвёртого порядка (или шестого порядка для машины с двойным белым каркасом), а механическая часть - системой второго порядка. Все электрические переменные и параметры относятся к статору, который обозначается первичными знаками в следующих уравнениях машины. Все величины статора и ротора находятся в произвольной двухосной системе отсчета (dq frame). Используемые подстрочные индексы определены в этой таблице.
Нижний индекс | Определение |
|---|---|
d | Количество по оси d |
q | количество по оси q |
r | Количество ротора (намотанный ротор или одноконтактный) |
r1 | Количество ротора сепаратора 1 (двухконтактный) |
r2 | Количество ротора сепаратора 2 (двухконтактный) |
s | Количество статора |
l | Индуктивность утечки |
m | Намагничивающая индуктивность |

Vqs = Rsiqs + díqs/dt +
Vds = Rsids + d.ds/dt -
V 'qr = R' ri 'qr + d' qr/dt + (λ - λ r)
V 'dr = R' ri 'dr + dstart' dr/dt - (λ - startr)
Te = 1,5p (фдзики - фдзиды)

λ - Угловая скорость системы отсчета
startr - Электрическая угловая скорость
dwfqs = Lsiqs + Lmi 'qr
dwfds = Lsids + Lmi 'dr
start' qr = L 'ri' qr + Lmiqs
start' dr = L 'ri' dr + Lmids
Ls = Lls + Lm
L 'r = L' lr + Lm

Vqs = Rsiqs + díqs/dt +
Vds = Rsids + d.ds/dt -
0 = R 'r1i' qr1 + d
0 = R 'r1i' dr1 + dstart' dr1/dt - (λ -
0 = R 'r2i' qr2 + d
0 = R 'r2i' dr2 + d
Te = 1,5p (фдзики - фдзиды)

dwfqs = Lsiqs + Lm (i 'qr1 + i' qr2)
dwfds = Lsids + Lm (i 'dr1 + i' dr2)
start' qr1 = L 'r1i' qr1 + Lmiqs
start' dr1 = L 'r1i' dr1 + Lmids
start' qr2 = L 'r2i' qr2 + Lmiqs
start' dr2 = L 'r2i' dr2 + Lmids
Ls = Lls + Lm
L 'r1 = L' lr1 + Lm
L 'r2 = L' lr2 + Lm
ддтм м =
Параметры блока асинхронной машины определены в таблице. Все величины относятся к статору.
Параметры, общие для всех моделей | Определение |
|---|---|
Rs, Lls | Сопротивление статора и индуктивность утечки |
Lm | Намагничивающая индуктивность |
Ls | Полная индуктивность статора |
Vqs, iqs | напряжение и ток статора по оси q |
Vds, ids | Напряжение и ток статора d-оси |
ϕqs, ϕϕds | Потоки q-оси статора и d-оси |
ωm | Угловая скорость ротора |
Θm | Угловое положение ротора |
p | Количество пар полюсов |
ωr | Электрическая угловая скорость (startm × p) |
Θr | Угловое положение ротора электрооборудования (Startm × p) |
Те | Электромагнитный момент |
TM | Механический крутящий момент вала |
J | Комбинированный коэффициент инерции ротора и нагрузки. Для моделирования блокированного ротора установите значение «бесконечно». |
H | Комбинированный ротор и постоянная инерции нагрузки. Для моделирования блокированного ротора установите значение «бесконечно». |
F | Комбинированный коэффициент трения ротора и нагрузки |
Параметры, относящиеся к однооболочечному или намотанному ротору | Определение |
|---|---|
Л 'р | Суммарная индуктивность ротора |
R 'r, L' lr | Сопротивление ротора и индуктивность утечки |
V 'qr, i' qr | напряжение и ток ротора по оси q |
V 'dr, i' dr | напряжение и ток ротора по оси d |
«qr», «dr» | Потоки по оси q и оси d ротора |
Параметры, характерные для двухконтажного ротора | Определение |
|---|---|
R 'r1, L' lr1 | Сопротивление ротора и индуктивность утечки обоймы 1 |
R 'r2, L' lr2 | Сопротивление ротора и индуктивность утечки обоймы 2 |
L 'r1, L' r2 | Суммарные индуктивности ротора обоймы 1 и 2 |
i 'qr1, i' qr2 | q-осевой ток ротора обоймы 1 и 2 |
i 'dr1, i' dr2 | ток ротора d-оси сепаратора 1 и 2 |
«qr1», «dr1» | q-осевые и d-осевые потоки ротора сепаратора 1 |
«qr2», «dr2» | q-осевые и d-осевые потоки ротора сепаратора 2 |
Блоки асинхронной машины не содержат представления насыщения потоков утечки. Будьте осторожны при подключении идеальных источников к статору машины. Если выбрана подача питания на статор через трехфазный, Y-подключенный источник бесконечного напряжения, необходимо использовать три источника, соединенных в Y. Однако, если выбрано моделирование дельта-подключения источника, необходимо использовать только два источника, соединенных последовательно.

При использовании блоков асинхронной машины в дискретных системах может потребоваться подключение небольшой паразитной резистивной нагрузки на клеммах машины, чтобы избежать цифровых колебаний. Большие времена выборки требуют больших нагрузок. Оптимальная резистивная нагрузка пропорциональна времени выборки. С шагом 25 мкс в системе с частотой 60 Гц минимальная нагрузка составляет приблизительно 2,5% номинальной мощности машины. Например, асинхронная машина 200 МВА в энергосистеме, дискретизированная с временем выборки 50 мкс, требует приблизительно 5% резистивной нагрузки или 10 МВт. Если время выборки уменьшено до 20 мкс, достаточна резистивная нагрузка 4 МВт.
Клеммы статора блоков асинхронной машины обозначены буквами A, B и C. Клеммы ротора обозначены буквами a, b и C. Нейтральные соединения обмоток статора и ротора недоступны. Предполагается трехпроводное Y-соединение.
power_pwm В примере используется блок асинхронной машины в режиме двигателя. Пример состоит из асинхронной машины в системе управления скоростью с разомкнутым контуром.
Ротор машины имеет короткое замыкание, а статор питается ШИМ-инвертором, построенным из блоков Simulink и сопряженным с блоком асинхронной машины через блок источника управляемого напряжения. Инвертор использует синусоидальную широтно-импульсную модуляцию. Базовую частоту синусоидальной опорной волны устанавливают на 60 Гц, а частоту треугольной несущей - на 1980 Гц. Эта частота соответствует коэффициенту модуляции частоты mf 33 (60 Гц x 33 = 1980).
Машина 3 ВД подключена к постоянной нагрузке номинального значения (11,9 Н· м). Запускается и достигает заданной скорости 1,0 pu при t = 0,9 секунды.
Параметры машины такие же, как и у блока асинхронной машины SI, за исключением индуктивности утечки статора, которая установлена в два раза выше нормального значения для моделирования сглаживающего индуктора, размещенного между инвертором и машиной. Также для получения показанных результатов использовали стационарный опорный кадр.
Откройте окно power_pwm пример. В параметрах моделирования требуется небольшой относительный допуск из-за высокой скорости переключения инвертора.

Выполните моделирование и проверьте скорость и крутящий момент машины.

Первый график показывает скорость машины от 0 до 1725 об/мин (1,0 pu). На втором графике показан электромагнитный крутящий момент, создаваемый машиной. Поскольку статор питается ШИМ-инвертором, наблюдается шумный крутящий момент.
Однако этот шум не виден в скорости, потому что он отфильтрован инерцией машины, но его можно видеть в токах статора и ротора.

Посмотрите на выходной сигнал ШИМ-инвертора. Поскольку на шкале времени моделирования ничего интересного не видно, график концентрируется на последних моментах моделирования.

power_asm_sat пример иллюстрирует эффект насыщения блока асинхронной машины.

Два одинаковых трехфазных двигателя (50 л.с., 460 В и 1800 об/мин) моделируются с насыщением и без него для наблюдения воздействия насыщения на токи статора. В примере реализованы два различных моделирования.
Первое моделирование - это стационарное испытание без нагрузки. Эта таблица содержит значения параметров насыщения и измерений, полученных путем моделирования различных рабочих точек на насыщенном двигателе (без нагрузки и в установившемся состоянии).
Параметры насыщения | Измерения | ||
|---|---|---|---|
Всат (Vrms L-L) | Изат (пик А) | Vrms L-L | Is_A (пик А) |
- | - | 120 | 7.322 |
230 | 14.04 | 230 | 14.03 |
- | - | 250 | 16.86 |
- | - | 300 | 24.04 |
322 | 27.81 | 322 | 28.39 |
- | - | 351 | 35.22 |
- | - | 382 | 43.83 |
414 | 53.79 | 414 | 54.21 |
- | - | 426 | 58.58 |
- | - | 449 | 67.94 |
460 | 72.69 | 460 | 73.01 |
- | - | 472 | 79.12 |
- | - | 488 | 88.43 |
506 | 97.98 | 506 | 100.9 |
- | - | 519 | 111.6 |
- | - | 535 | 126.9 |
- | - | 546 | 139.1 |
552 | 148.68 | 552 | 146.3 |
- | - | 569 | 169.1 |
- | - | 581 | 187.4 |
598 | 215.74 | 598 | 216.5 |
- | - | 620 | 259.6 |
- | - | 633 | 287.8 |
644 | 302.98 | 644 | 313.2 |
- | - | 659 | 350 |
- | - | 672 | 383.7 |
- | - | 681 | 407.9 |
690 | 428.78 | 690 | 432.9 |
Следующий график иллюстрирует эти результаты и показывает точность модели насыщения. Измеренные рабочие точки соответствуют кривой, построенной на основе данных параметров насыщения.

Можно наблюдать другие эффекты насыщения токов статора, выполняя моделирование с заблокированным ротором или с различными значениями крутящего момента нагрузки.
[1] Краузе, П.К., О. Васынчук и С.Д. Sudhoff, Анализ электрических машин, IEEE ® Press, 2002.
[2] Мохан, Н., Т. М. Унделанд и В. П. Роббинс, Power Electronics: Converters, Applications, and Design, John Wiley & Sons, Inc., Нью-Йорк, 1995 год, раздел 8.4.1.