CC BY
149
E. A. Novikov
DEVELOPMENT OF (4,2)-METHOD OF ORDER 3 FOR SOLVING STIFF PROBLEMS
We obtain the coefficients of (4,2)-method for solving stiff systems in which both the main and the intermediate numerical schemes are L-stable. We construct an accuracy control inequality. Numerical tests confirming efficiency and workability of the constructed variable step algorithm were held.
Keywords: stiff problem, (m,k)-methods, accuracy control.
© Новиков Е. А., 2011
УДК 621.318.562.5
А. Н. Пахомов, М. Ф. Коротков, А. А. Федоренко
МОДАЛЬНОЕ УПРАВЛЕНИЕ ЭЛЕКТРОПРИВОДОМ ПЕРЕМЕННОГО ТОКА
Приведена методика синтеза модальных регуляторов координат векторной системы «преобразователь частоты-асинхронный двигатель» методом стандартных уравнений. Дана оценка качества процессов регулирования координат путем анализа результатов имитационного моделирования системы в среде МаіЬаЬ.
Ключевые слова: модальный регулятор, электропривод переменного тока, векторная система.
Теория систем векторного управления частотно регулируемого асинхронного электропривода разработана достаточно полно [1; 2]. Контуры регулирования в системе векторного управления выполняются, как правило, в соответствии с принципами подчиненного регулирования координат, что ограничивает их быстродействие и, как следствие, точность в динамических режимах. Обеспечить предельное быстродействие и точность в динамических режимах возможно, снабдив систему так называемым модальным регулятором, построенным на основе суммирования обратных связей по вектору состояния. Вопросы построения таких регуляторов применительно к векторным системам асинхронного электропривода в литературе отражены недостаточно, в связи с чем в данной работе ставится задача разработки методики проектирования модальных регуляторов векторных систем частотно управляемого асинхронного электропривода и проверка ее эффективности с помощью имитационного моделирования в среде МаіЬаЬ.
В качестве объекта управления принята получившая наибольшее распространение система «преобразователь частоты с автономным инвертором напряжения с широтно-импульсной модуляцией - асинхронный двигатель». Поскольку у таких преобразователей автономный инвертор формирует не только частоту, но и амплитуду выходного напряжения, влияние звена постоянного тока на динамические свойства системы при синтезе можно не учитывать. Кроме того, частота модуляции современных преобразователей весьма высока, что позволяет пренебречь также его дискретными свойствами. Изложенное дает возможность представить в первом приближении преобразо-
ватель частоты безынерционным линейным звеном с коэффициентом передачи кп.
Математическая модель асинхронного двигателя (АД) с учетом общепринятых допущений [2; 3], в декартовой системе координат и-у, вращающейся с произвольной скоростью юк, в форме Коши имеет следующий вид [3]:
di
1
dt RT
+ k2 Pn RT
шШ 2v +rokilv;
T2
1
-u,„---------ilv +-
dilv
dt RX T
k2 Pn
T2
RT
2u -“kilu;
= R2k2ilu - ^ Ш2u + К - Pnш)Ш2v ; dt T 2
= R2k2hv - T Ш2v - (“k - Pn®)Y2u ;
dt T
^ = Шх i Ш 2u - ^ ilu Ш 2v - 1M c
dt
2J
2 J lu 2v J
(l)
рующих векторов напряжения статора и1 , тока статора іі, потокосцепления ротора Т2 соответственно, на оси и и у декартовой системы координат; ю = юэл / рп - угловая скорость вращения ротора;
где ulu , ulv , ilu
lv
юэл - угловая скорость вращения электромагнитного
поля; Мс - момент статического сопротивления на
валу двигателя [1-3].
Уравнения динамики АД в векторно-матричной форме имеют вид
XT = F(X, G).
(2)
(
X =
Л
*lv
Т
1 2u
Т
и G =
Си ^
lu
(3)
Линеаризованную модель АД будем искать в виде стандартной векторно-матричной формы уравнений состояния
AX = AAX + BAG,
(4)
A =
dF
dX
я dF
и B =--------
g=g0 dG
X = X,
(5)
ления ротора Т 2 = Т 2
Т 2v = 0. Переменная
A =
где X и G - векторы переменных состояния и входных воздействий, причем
Т O k2 T2 RэTэ O
O l k2 pn k2 pn
Тэ RэTэ~0 RэTэ
R2 k2 O -l T2 O
O 3pn k2 Т л T 1 2u 0 3pn k2 . л t lvO O
Т
2 J
(б)
Г kn O
RэTэ
0 kn RэTэ
0 O
V 0 O
где А и В - матрицы динамики и входа соответственно, которые необходимо определить; АХ и АО -приращения векторов переменных состояния и входных воздействий соответственно.
Применим метод линеаризации на основе разложения уравнений (1) в ряд Тейлора. Тогда матрицы А и В определятся в виде
2 J
и B =
Расчет модальных регуляторов производится отдельно для двух каналов управления: канала регулирования потокосцепления ротора АД и канала регулирования скорости АД. Таким образом, строится два контура регулирования аналогично построению систем двухзонного регулирования скорости электропривода постоянного тока [1].
Разобьем матрицы динамики и входа (6) на две пары для синтеза модальных регуляторов по каждому из каналов регулирования.
Матрицы канала регулирования потокосцепления ротора:
( _1_
a, =
где Х0 и 00 - переменные состояния и входные воздействия в точке установившегося режима (центра разложения), значения которых определяются по уравнениям статики.
Взяв производные от уравнений (1) по соответствующим переменным, видим, что матрицы А и В имеют высокий порядок. Для их упрощения необходимо выбрать угловую скорость вращения юк системы координат и-у и сориентировать ось и по определенной векторной переменной АД.
Целесообразно использовать систему координат и-у с угловой скоростью вращения магнитного поля юк = ю0эл, ориентированной по вектору потокосцеп-
R2 k2
b, =
bn
Г К Л
R3t3
O
Матрицы канала регулирования скорости вращения ротора АД:
A2 =
T
k2 ph RэTэ
Т
1 о
3?п k2 2 J
O
Г k Л
П
и B2 =
0
V У
юк = ю0эл в этом случае определяется зависимостью га0эл = /Оь, ^2и,ю), т. е. она уже не является входным воздействием.
Синтез модального регулятора проведем для управляющего воздействия, т. е. положим Мс = 0 .
Наличие на входе двигателя преобразователя частоты учтем в матрице входа коэффициентом кп. Если скомпенсировать влияние перекрестных обратных связей по току статора, получим матрицы
Характеристический полином замкнутого по вектору состояния объекта управления можно записать в виде
G (p ) = |p •E-(At.-Bt-K)| = D (p), (7)
где K - коэффициент модального регулятора; E - единичная матрица; D(p)- желаемый характеристический полином.
Поскольку корни характеристического уравнения D(p)=0 определяют динамические свойства системы, замкнутой через модальный регулятор, то для получения желаемых динамических показателей требуется,
7l
чтобы замкнутая обратная связь по вектору состояния обеспечивала желаемые корни характеристического полинома. Таким образом, необходимо подобрать матрицу К модального регулятора так, чтобы получить в комплексной плоскости желаемое распределение корней характеристического полинома.
Синтез модального регулятора выполняется методом стандартных уравнений в нормированной форме [1]. Для получения нормированной формы исходное уравнение п-го порядка делится на коэффициент при старшей степени Б0 и приобретает корневую форму:
рп + Сх -Юб • рп-1 +... + С,. Ч • рп-‘ +... + Сп_1 •< = 0, (8) где юб - базовая частота (среднегеометрический корень), юб = :
С=
Di
D0 -Юб
Делением уравнения (8) на осуществляется пе-
реход к нормированному уравнению:
Со • 5п + С! • 5п-1 +... + С, • +... + Сп-1 • 5 +1 = 0 (9)
с относительным оператором S = -
d d т
и относи-
(
a,, + a22 + Cq
Л
a^2 + С, • a22 + С, + a,2
(11)
Tn.n = 4,8 Юб-
В схеме приняты следующие обозначения.
Модель АД, реализованная по уравнениям состояния (1) - AD.
Коэффициенты усиления, обеспечивающие требуемые потокосцепление и скорость ротора АД при наличии ОС - k1u, k2u.
Блоки преобразователя координат - m1, m3, m5, реализованные по уравнениям
ix = ia • cos9 + ip • sin9;j
iy = ia • sin9-ip • cos9;J
ua = ux • cos9-uy • sin9;1
Mp = ux • sin9 + uy • cos9; J
F„
cos 0 =
cos 0 =
jFxr 2 + Fyr 2
Fr
jFxr 2 + Fyr 2
F
F = J F.2 + Fy.2.
Блок компенсаций перекрестных связей лизованный по уравнениям
BK, реа-
тельным временем Х = Юб • ґ.
Рассчитаем коэффициенты обратных связей (ОС) при настройке системы на модульный оптимум: перерегулирование ст = 4,5 %, а относительное время переходного процесса Тп.п = 2,8 ®б"‘. Нормированное желаемое уравнение имеет вид [1]
£2 + ^2 • £ +1. (10)
После перехода к абсолютным единицам подставим желаемый характеристический полином в выражение (7) и приравняем значения коэффициентов при одинаковых степенях р, получим вектор коэффициентов модального регулятора для канала регулирования потокосцепления ротора:
U у* = U y,x - iy'
U уу = Uу,у + .х
• RlITl
kn
Я21 • Ь11
Аналогично производится расчет второго канала регулирования угловой скорости вращения ротора АД. В качестве примера принята настройка по биномиальному разложению (5 + 1)п: перерегулирование ст = 0 %, а время переходного процесса
Приведем расчет переходных процессов в замкнутой системе асинхронного электропривода при помощи пакета имитационного моделирования Matlab фирмы The MathWorks, Inc. Воспользуемся соответствующей структурной схемой замкнутой системы асинхронного электропривода в неподвижной системе координат (рис. 1), в которой переменные ОЭМ характеризуют поведение реальных координат АД.
где иух, иуу - управляющие воздействия на входе системы.
Переходные процессы в представленной модели изображены на рис. 2. Они совпадают с переходными процессами во вращающейся системе координат, что подтверждает правильность преобразования координат. Сначала подается задание на потокосцепление ротора АД, которое успешно отрабатывает канал регулирования потокосцепления, после чего производится скачкообразное изменение задания на скорость вращения ротора и включается в работу второй канал. И, наконец, на третьем участке осуществляется моделирование наброса нагрузки на валу АД.
Характер переходных процессов соответствует заданным настройкам. Следовательно, регуляторы позволяют обеспечить требуемый уровень скорости и потокосцепления как в неподвижной, так и во вращающейся системах координат.
Библиографические ссылки
1. Терехов В. М., Осипов О. И. Системы управления электроприводов : учеб. для студ. высш. учеб. заведений ; под ред. В. М. Терехова. М. : Академия, 2005.
2. Соколовский Г. Г. Электроприводы переменного тока с частотным регулированием : учеб. для студ. высш. учеб. заведений. М. : Академия, 2006.
3. Карагодин М. С., Федоренко А. А. Уравнения динамики частотно-управляемых электроприводов : учеб. пособие / Краснояр. политехн. ин-т. Красноярск, 1985.
п
Рис. 1. Модель электропривода в неподвижной системе координат
Модуль потокосцепления ротора
Рис. 2. Графики переходных процессов в неподвижной системе координат
A. N. Pakhomov, M. F. Korotkov, A. A. Fedorenko ALTERNATING CURRENT ELECTRIC DRIVE MODAL CONTROL
The technique of synthesis of modal regulators of coordinates of vector system “the converter of frequency-asynchronous engine " is resulted by a method of the standard equations. The estimation of quality ofprocesses of regulation of coordinates by the analysis of results of simulation modeling of the system in the environment of MatLab is given.
Keywords: modal regulator, alternating current electric drive, vector system.
© Пахомов А. Н., Коротков М. Ф., Федоренко А. А., 2011
УДК 004.912
К. В. Полянский
ПОСТРОЕНИЕ ФРЕЙМОВОЙ МОДЕЛИ ПЕРЕВОДА С ИСПОЛЬЗОВАНИЕМ КЛАСТЕРИЗАЦИИ ТЕРМОВ
Рассмотрена фреймовая модель представления знаний в 1Р-системах машинного перевода. Предложен алгоритм сегментации исходного и целевого текста через связь. Проанализированы различные методы кластеризации термов, определены наиболее эффективные из них для разбиения текста на кластеры.
Ключевые слова: машинный перевод, сегментация текста, кластеризация термов, фреймовая модель.
Важным этапом в 1Р-переводе (машинном переводе, использующем ресурсы информационно-поисковых систем) на стадии синтеза является сопоставление исходного текста (ИЯ-текста) и релевантных текстов на целевом языке (ЦЯ-текстов), выявление в них схожих сегментов. Процесс такого сопоставления выполняется в несколько шагов:
1) сегментация текста;
2) кластеризация сегментов;
3) построение фреймовой модели структуры текста. Рассмотрим каждый шаг подробнее.
Сегментация текста. Для анализа структуры
предложений ИЯ- и ЦЯ-текстов необходимо поделить эти предложения на логические сегменты, где каждый сегмент будет семантически самостоятельной единицей. Сегментом назовем непрерывный фрагмент тек-
