ПАРАМЕТРИЧЕСКАЯ ИДЕНТИФИКАЦИЯ ВЕНТИЛЬНОГО ЭЛЕКТРОПРИВОДА АЗИМУТАЛЬНОЙ ОСИ ТЕЛЕСКОПА ТРАЕКТОРНЫХ ИЗМЕРЕНИЙ
А.Г. Ильина
УДК 681.587.73
ПАРАМЕТРИЧЕСКАЯ ИДЕНТИФИКАЦИЯ ВЕНТИЛЬНОГО ЭЛЕКТРОПРИВОДА АЗИМУТАЛЬНОЙ ОСИ ТЕЛЕСКОПА ТРАЕКТОРНЫХ
ИЗМЕРЕНИЙ
А.Г. Ильина
Рассмотрена процедура параметрической идентификации безредукторного высокомоментного прецизионного электропривода опорно-поворотного устройства телескопа траекторных измерений, проводимая в целях синтеза параметров регуляторов его подчиненной САР. В результате проведенной идентификации вычислены параметры теоретической модели электропривода, позволившие в достаточной для решаемой задачи степени приблизить свойства математической модели к свойствам объекта управления. Ключевые слова: прецизионный электропривод, параметрическая идентификация.
Введение
Синтез прецизионной системы управления телескопом траекторных измерений (ТТИ) невозможен без адекватной объекту математической модели моментного электропривода опорно-поворотного устройства (ОПУ). Развитие машинно-ориентированных методов проектирования систем автоматического регулирования (САР) позволяет использовать для получения таких моделей методы параметрической и непараметрической идентификаций.
Основными особенностями моментных электроприводов ОПУ ТТИ являются [1, 2]: широкий диапазон скоростей слежения при малых значениях среднеквадратичных ошибок наведения, инфранизкие скорости движения оптической оси ТТИ (единицы угловых секунд за секунду), использование уникальных прецизионных датчиков координат, нежесткость конструкции осей ОПУ ТТИ и, как следствие, возникновение механического резонанса на частотах от единиц Гц. Поэтому в качестве электромеханических преобразователей в исполнительных осях ОПУ рассматриваемых ТТИ (рис. 1, а) с целью максимального упрощения кинематических передач выбираются встраиваемые в ОПУ высокомоментные синхронные машины с постоянными магнитами на роторе, используемые в режиме вентильного двигателя (ВД) [3, 4]. Эти машины отличаются малым энергопотреблением на инфранизких скоростях и большим отношением вращающего момента к моменту инерции ротора, определяющим предельное быстродействие машины. Термин «моментный» означает, что двигатели спроектированы специально для режима вращения с малой скоростью (строго говоря, для пускового режима), когда можно не учитывать все потери, связанные с быстрым вращением ротора.
В настоящей статье ставится задача получения адекватной математической модели моментного электропривода азимутальной оси совместно с его нагрузкой – ОПУ ТТИ, описывающей этот сложный механический объект управления в виде двух сосредоточенных масс, соединенных нежесткой связью, т.е. системы, обладающей одним механическим резонансом.
Решение поставленной задачи обеспечивается посредством определения части коэффициентов уравнений модели с помощью параметрической идентификации. В век повсеместного внедрения информационных технологий сложнейший математический аппарат параметрической и непараметрической идентификации, превращенный Mathworks при сотрудничестве пионера современной идентификации и всемирно известного специалиста в области теории автоматического управления Л. Льюнга в доступный и удобный инструмент, не является более серьезным препятствием к ее практическому применению в задачах синтеза систем управления. Однако успешное применение параметрической идентификации требует задания априорных значений искомых параметров, для чего, учитывая сложность конструкции рассматриваемого объекта управления, целесообразно подвергнуть анализу 3-d модель конструкции (рис. 1, б) для нахождения значений параметров неэлектрической природы. 3-d модель может быть реализована в любой из множества сред разработки, таких как SolidWorks или AutoCAD, а анализ конструкции на частотные и другие свойства легко реализуется в пакетах типа ANSYS.
В статье представлена последовательность операций, проводимая с применением упомянутых пакетов и сред и позволяющая получить адекватную ОПУ ТТИ двухмассовую математическую модель.
Рассмотрим математическую модель вентильного двигателя (ВД), работающего на двухмассовую нагрузку с нежесткой связью. Известно, что математическая модель ВД может быть без существенной потери информативности упрощена и представлена моделью двигателя постоянного тока (ДПТ) [1, 3]. Тогда система уравнений электропривода, работающего на нагрузку с нежестким валом, будет иметь вид
Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2010, № 4(68)
33
ПАРАМЕТРИЧЕСКАЯ ИДЕНТИФИКАЦИЯ ВЕНТИЛЬНОГО ЭЛЕКТРОПРИВОДА ...
L
di dt
i
c e
1
u;
J1
d 1 dt
cmi
d 1 dt
1;
db dt
c 12
(1
M c1 2 );
M12 M каб.пер.; M каб.пер.
M 12
b m(1
2 );
K
1;
J
2
d2 dt
M12 M c2 ;
d
2
dt
2 ,
(1)
где L , R – индуктивность и сопротивление фазы ВД, i , u – ток в фазе ВД и напряжение, приложенное к
фазе ВД, ce – конструктивная постоянная по ЭДС, cm – постоянная по моменту, 1 , 2 – скорости вра-
щения первой и второй массы, M , M c1 , M c2 , M12 – электромагнитный момент двигателя, моменты сопротивления вращению первой и второй масс и момент, создаваемый нежесткой связью, Mкаб.пер. – момент возмущения, создаваемый кабельным переходом, связанный с углом поворота первой массы коэффициен-
том жесткости K , c12 – коэффициент жесткости упругой связи, m – коэффициент демпфирования, J1 , J2
– приведенные к оси вращения моменты инерции первой и второй масс, 1 , 2 – углы поворота первой и
второй масс,
zp
– число пар полюсов двигателя,
cos arctg z p1Tf
R2 z p1RTf 2
1 R
– коэффициент приведения
ВД
к
ДПТ
(для
моментных
двигателей,
работающих
на
инфранизких
скоростях
1 R
)[3].
Параметры,
требующие определения для приведения математической модели (1) к реальному объекту – J1 , J2 , ce ,
c12 , cm , m , тогда как известными можно считать подлежащие измерению R и Tf .
При проведении параметрической идентификации требуется определить примерные (так называемые начальные или априорные) значения искомых параметров, так как при большом количестве неизвестных возможно множество решений поставленной задачи, и полученная таким образом модель может
оказаться недостоверной. Априорные значения параметров J1 , J2 и c12 могут быть получены по 3-d
модели (рис. 1, б, в) ОПУ ТТИ и его реальной частотной характеристике (рис. 2). На рис. 1 представлены фотография опытного образца ОПУ ТТИ (а) и его 3-d модель (б–г), реали-
зованная в SolidWorks. САПР Solid Works позволяет провести расчет приведенных к азимутальной оси моментов инерции
«сборок» (рис. 1, в, и 1, г) и, таким образом, установить приближенные значения моментов инерции J1 и
J2 обеих сосредоточенных масс. Из выражения, определяющего резонансную частоту двухмассового
механизма, 0р
c12 (J1 J2 ) J1 J 2
4 (J1 J2 ) / J1 3 [6], и значения нижней частоты резонанса на экспери-
ментальной АЧХ опытного образца ОПУ ТТИ (рис. 2) находим значение жесткости упругой связи с12.
Априорные значения ce и cm , связанные соотношением cm 1,5ce [3], можно грубо оценить как отно-
шение выходного и входного сигналов ОПУ ТТИ в установившемся режиме работы, так как при жесткой
связи ( c12 ) обеих масс передаточная функция вентильного электропривода без учета влияния момента сопротивления кабельного перехода примет вид
W p
p u p
Tf JsR cecm
1 / ce
p2
JsR cecm
p
, 1
где Js J1 J2 – суммарный момент инерции.
34 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2010, № 4(68)
А.Г. Ильина
Рис. 1. Опытный образец ОПУ ТТИ (а), 3-d модель ОПУ-ТТИ (продольный разрез)(б), 3-d модель второй
массы, характеризующейся моментом инерции J2 (в), 3-d модель первой массы, характеризующейся моментом инерции J1 (г)
Параметрическая идентификация ОПУ ТТИ проводится методом Estimating ODE parameters [5]. Этот метод может применяться для проведения параметрической идентификации в нашем случае, так как система дифференциальных уравнений, описывающих объект, известна. В качестве входного сигнала из ряда стандартных тестовых сигналов для решения поставленной задачи выбран меандр напряжения. При других стандартных тестовых сигналах (импульсная характеристика или белый шум в объекте) могут возникнуть колебания во всем частотном спектре, что, учитывая наше представление объекта двухмассовой моделью с одним резонансом, является некорректным. При использовании линейно нарастающего входного задания частотная характеристика будет искажена демпфирующим воздействием кабельного перехода.
Рассмотрим процесс идентификаций методом Estimating ODE parameters. Для проведения иденти-
фикации модель (1) следует представить в матричной записи (здесь e t – белый шум):
x t Ax t Bu t Ke t ,
y t Cx t Du t e t ,
(2)
x 0 x0.
Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2010, № 4(68)
35
ПАРАМЕТРИЧЕСКАЯ ИДЕНТИФИКАЦИЯ ВЕНТИЛЬНОГО ЭЛЕКТРОПРИВОДА ...
Рис. 2. Амплитудно-частотная характеристика опытного образца ОПУ ТТИ (рис. 1, а)
Так как на рис. 3 наблюдается низкая зашумленность сигнала скорости, то влиянием шума измерений и шума объекта можно пренебречь. Матрицы исследуемой системы имеют следующий вид:
y i 1 1 2 2 T ; u U b 1_ отн b ;
1 Tf cm J1
A 0
0
0
ce L
m J1 1
m J2 0
0 c12 J1
0 c12 J 2
0
00
m J1
c12 J1
0 0 ;
m J2 c12 J 2
1 0
1 / L
0
0
0
0
0
0
0
B 0
0
Mc1 / J1 K 00
0 ;
0
0 0 0 Mc2 / J1
C 0 1 0 0 0, K 0 0 0 0 0, D 0 0 0 0 ;
x0 0 0 0 0 0 .
Для вычисления значений параметров записанной в матричном виде системы уравнений в рабочей
области Matlab задаются следующие команды:
model=idgrey(‘OPU’,parameter,’c’,0.001,0);model=pem(dannie,model);,
(3)
где “OPU” представляет собой функцию, задающую структуру модели:
function [A,B,C,D,K,x0] = OPU(par,T,CDmfile,aux)
A=[(-1/Tf), -par(1)/L, 0, 0, 0; 1.5*par(1)/par(2), -par(3)/ par(2), -par(4)/ par(2), par(3)/ par(2),
par(4)/par(2); 0, 1, 0, 0, 0;0, par(3)/par(5), par(4)/ par(5), - par(3)/ par(5), - par(4)/ par(5); 0, 0, 0, 1, 0];
B=[1/L 0 0 0;0 –par(6)/J1 -par(7) 0;0 0 0 0;0 0 0 –par(8)/J2; 0 0 0 0];
C=[0, 1, 0, 0, 0];
D=[0 0 0 0];
K=zeros(5,1);
x0=[0; 0; 0; 0; 0];
“parameter” – файл, содержащий априорные значения искомых параметров; dannie=iddata(s,u,Td) –
входные данные идентификации, сформированные из сигналов: s – массив данных скорости объекта, u –
массив данных управляющего напряжения, Td – шаг измерений в секундах.
36 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2010, № 4(68)
А.Г. Ильина
В результате выполнения команд (3) с учетом принятых допущений и при условии, что искомое решение существует, формируется модель заданной структуры с параметрами, близкими по значениям к указанным в файле parameter.m, и обладающая динамическими свойствами, сходными с реальным объектом. На рис. 3 представлены графики изменения угловой скорости при подаче входного меандра напряжения, полученные на опытном образце ОПУ ТТИ и при моделировании системы уравнений (1) где в качестве значений параметров используются полученные параметрической идентификацией методом Estimating ODE parameters.
Рис. 3. Сигналы изменения скорости разомкнутого вентильного электропривода при реакции на входной сигнал управления типа меандр: 1 – опытный образец ОПУ ТТИ (рис. 1, а), 2 – найденная математическая
модель
Как уже отмечалось, конечной целью получения адекватной математической модели объекта является синтез его системы управления. Продемонстрируем работу замкнутого контура скорости опытного образца ОПУ ТТИ (рис. 1, а), настроенного на технический оптимум (ТО), причем параметры ПИ регулятора контура рассчитаны на основании динамических свойств найденной с помощью идентификации модели. На рис. 4 представлены графики эталонной кривой отработки настроенной на ТО САР скорости входного ступенчатого воздействия и угловая скорость в САР скорости вращения оптической оси ОПУ ТТИ при реакции на тот же сигнал.
Рис. 4. Эталонная кривая (кривая 1) переходного процесса в замкнутом контуре скорости математической модели (1), настроенном на технический оптимум при отработке входного ступенчатого задающего сигнала, и аналогичный процесс в контуре скорости опытного образца ОПУ ТТИ (кривая 2)
Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2010, № 4(68)
37
ДИНАМИЧЕСКОЕ НАБЛЮДЕНИЕ НЕЛИНЕЙНЫХ ДВОИЧНЫХ...
Выводы
Проведена параметрическая идентификация, в результате которой получена математическая модель, удовлетворительно описывающая поведение реального объекта – опытного образца ОПУ ТТИ.
Применение разработанной модели при синтезе САР опытного образца ОПУ ТТИ приводит к гарантированному обеспечению желаемых качественных показателей проектируемой системы, а также значительно сокращает время на ее настройку.
Литература
1. Кротенко В.В., Толмачев В.А., Томасов В.С., Синицын В.А. Синтез микропроцессорной системы управления электропривода опорно-поворотного устройства // Изв. вузов. Приборостроение. – 2004. – Т. 47. – № 11. – С. 23–30.
2. Гурьянов А.В., Кротенко В.В. Параметрический синтез цифровых систем управления с широтноимпульсными преобразователями и эталонными моделями // Изв. вузов. Приборостроение. – 2006. – Т. 49. – № 3. – С. 17–24.
3. Микеров А.Г. Управляемые вентильные двигатели малой мощности: Учеб. пособие. – СПб: Издво ГЭТУ, 1997. – 64 с.
4. Takeshita T., Matsui N. Sensorless Brushless DC Motor Drive with EMF Constant Identifier // IECON-94. –1994. – V.1. – P. 8–13.
5. Ljung L. System Identification Toolbox 7 User’s Guide [Электронный ресурс]. – Режим доступа: http://www.mathworks.com/access/helpdesk/help/pdf_doc/ident/ident.pdf, своб.
6. Шрейнер Р.Т. Системы подчиненного регулирования электроприводов: учеб. пособие. – Екатеринбург: Изд-во ГОУ ВПО «Рос. гос. проф.-пед. ун-т», 2008. – 279 с.
Ильина Аглая Геннадьевна
– Санкт-Петербургский государственный университет информационных технологий, механики и оптики, научный сотрудник, dvanoska@mail.ru
38 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2010, № 4(68)
УДК 681.587.73
ПАРАМЕТРИЧЕСКАЯ ИДЕНТИФИКАЦИЯ ВЕНТИЛЬНОГО ЭЛЕКТРОПРИВОДА АЗИМУТАЛЬНОЙ ОСИ ТЕЛЕСКОПА ТРАЕКТОРНЫХ
ИЗМЕРЕНИЙ
А.Г. Ильина
Рассмотрена процедура параметрической идентификации безредукторного высокомоментного прецизионного электропривода опорно-поворотного устройства телескопа траекторных измерений, проводимая в целях синтеза параметров регуляторов его подчиненной САР. В результате проведенной идентификации вычислены параметры теоретической модели электропривода, позволившие в достаточной для решаемой задачи степени приблизить свойства математической модели к свойствам объекта управления. Ключевые слова: прецизионный электропривод, параметрическая идентификация.
Введение
Синтез прецизионной системы управления телескопом траекторных измерений (ТТИ) невозможен без адекватной объекту математической модели моментного электропривода опорно-поворотного устройства (ОПУ). Развитие машинно-ориентированных методов проектирования систем автоматического регулирования (САР) позволяет использовать для получения таких моделей методы параметрической и непараметрической идентификаций.
Основными особенностями моментных электроприводов ОПУ ТТИ являются [1, 2]: широкий диапазон скоростей слежения при малых значениях среднеквадратичных ошибок наведения, инфранизкие скорости движения оптической оси ТТИ (единицы угловых секунд за секунду), использование уникальных прецизионных датчиков координат, нежесткость конструкции осей ОПУ ТТИ и, как следствие, возникновение механического резонанса на частотах от единиц Гц. Поэтому в качестве электромеханических преобразователей в исполнительных осях ОПУ рассматриваемых ТТИ (рис. 1, а) с целью максимального упрощения кинематических передач выбираются встраиваемые в ОПУ высокомоментные синхронные машины с постоянными магнитами на роторе, используемые в режиме вентильного двигателя (ВД) [3, 4]. Эти машины отличаются малым энергопотреблением на инфранизких скоростях и большим отношением вращающего момента к моменту инерции ротора, определяющим предельное быстродействие машины. Термин «моментный» означает, что двигатели спроектированы специально для режима вращения с малой скоростью (строго говоря, для пускового режима), когда можно не учитывать все потери, связанные с быстрым вращением ротора.
В настоящей статье ставится задача получения адекватной математической модели моментного электропривода азимутальной оси совместно с его нагрузкой – ОПУ ТТИ, описывающей этот сложный механический объект управления в виде двух сосредоточенных масс, соединенных нежесткой связью, т.е. системы, обладающей одним механическим резонансом.
Решение поставленной задачи обеспечивается посредством определения части коэффициентов уравнений модели с помощью параметрической идентификации. В век повсеместного внедрения информационных технологий сложнейший математический аппарат параметрической и непараметрической идентификации, превращенный Mathworks при сотрудничестве пионера современной идентификации и всемирно известного специалиста в области теории автоматического управления Л. Льюнга в доступный и удобный инструмент, не является более серьезным препятствием к ее практическому применению в задачах синтеза систем управления. Однако успешное применение параметрической идентификации требует задания априорных значений искомых параметров, для чего, учитывая сложность конструкции рассматриваемого объекта управления, целесообразно подвергнуть анализу 3-d модель конструкции (рис. 1, б) для нахождения значений параметров неэлектрической природы. 3-d модель может быть реализована в любой из множества сред разработки, таких как SolidWorks или AutoCAD, а анализ конструкции на частотные и другие свойства легко реализуется в пакетах типа ANSYS.
В статье представлена последовательность операций, проводимая с применением упомянутых пакетов и сред и позволяющая получить адекватную ОПУ ТТИ двухмассовую математическую модель.
Рассмотрим математическую модель вентильного двигателя (ВД), работающего на двухмассовую нагрузку с нежесткой связью. Известно, что математическая модель ВД может быть без существенной потери информативности упрощена и представлена моделью двигателя постоянного тока (ДПТ) [1, 3]. Тогда система уравнений электропривода, работающего на нагрузку с нежестким валом, будет иметь вид
Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2010, № 4(68)
33
ПАРАМЕТРИЧЕСКАЯ ИДЕНТИФИКАЦИЯ ВЕНТИЛЬНОГО ЭЛЕКТРОПРИВОДА ...
L
di dt
i
c e
1
u;
J1
d 1 dt
cmi
d 1 dt
1;
db dt
c 12
(1
M c1 2 );
M12 M каб.пер.; M каб.пер.
M 12
b m(1
2 );
K
1;
J
2
d2 dt
M12 M c2 ;
d
2
dt
2 ,
(1)
где L , R – индуктивность и сопротивление фазы ВД, i , u – ток в фазе ВД и напряжение, приложенное к
фазе ВД, ce – конструктивная постоянная по ЭДС, cm – постоянная по моменту, 1 , 2 – скорости вра-
щения первой и второй массы, M , M c1 , M c2 , M12 – электромагнитный момент двигателя, моменты сопротивления вращению первой и второй масс и момент, создаваемый нежесткой связью, Mкаб.пер. – момент возмущения, создаваемый кабельным переходом, связанный с углом поворота первой массы коэффициен-
том жесткости K , c12 – коэффициент жесткости упругой связи, m – коэффициент демпфирования, J1 , J2
– приведенные к оси вращения моменты инерции первой и второй масс, 1 , 2 – углы поворота первой и
второй масс,
zp
– число пар полюсов двигателя,
cos arctg z p1Tf
R2 z p1RTf 2
1 R
– коэффициент приведения
ВД
к
ДПТ
(для
моментных
двигателей,
работающих
на
инфранизких
скоростях
1 R
)[3].
Параметры,
требующие определения для приведения математической модели (1) к реальному объекту – J1 , J2 , ce ,
c12 , cm , m , тогда как известными можно считать подлежащие измерению R и Tf .
При проведении параметрической идентификации требуется определить примерные (так называемые начальные или априорные) значения искомых параметров, так как при большом количестве неизвестных возможно множество решений поставленной задачи, и полученная таким образом модель может
оказаться недостоверной. Априорные значения параметров J1 , J2 и c12 могут быть получены по 3-d
модели (рис. 1, б, в) ОПУ ТТИ и его реальной частотной характеристике (рис. 2). На рис. 1 представлены фотография опытного образца ОПУ ТТИ (а) и его 3-d модель (б–г), реали-
зованная в SolidWorks. САПР Solid Works позволяет провести расчет приведенных к азимутальной оси моментов инерции
«сборок» (рис. 1, в, и 1, г) и, таким образом, установить приближенные значения моментов инерции J1 и
J2 обеих сосредоточенных масс. Из выражения, определяющего резонансную частоту двухмассового
механизма, 0р
c12 (J1 J2 ) J1 J 2
4 (J1 J2 ) / J1 3 [6], и значения нижней частоты резонанса на экспери-
ментальной АЧХ опытного образца ОПУ ТТИ (рис. 2) находим значение жесткости упругой связи с12.
Априорные значения ce и cm , связанные соотношением cm 1,5ce [3], можно грубо оценить как отно-
шение выходного и входного сигналов ОПУ ТТИ в установившемся режиме работы, так как при жесткой
связи ( c12 ) обеих масс передаточная функция вентильного электропривода без учета влияния момента сопротивления кабельного перехода примет вид
W p
p u p
Tf JsR cecm
1 / ce
p2
JsR cecm
p
, 1
где Js J1 J2 – суммарный момент инерции.
34 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2010, № 4(68)
А.Г. Ильина
Рис. 1. Опытный образец ОПУ ТТИ (а), 3-d модель ОПУ-ТТИ (продольный разрез)(б), 3-d модель второй
массы, характеризующейся моментом инерции J2 (в), 3-d модель первой массы, характеризующейся моментом инерции J1 (г)
Параметрическая идентификация ОПУ ТТИ проводится методом Estimating ODE parameters [5]. Этот метод может применяться для проведения параметрической идентификации в нашем случае, так как система дифференциальных уравнений, описывающих объект, известна. В качестве входного сигнала из ряда стандартных тестовых сигналов для решения поставленной задачи выбран меандр напряжения. При других стандартных тестовых сигналах (импульсная характеристика или белый шум в объекте) могут возникнуть колебания во всем частотном спектре, что, учитывая наше представление объекта двухмассовой моделью с одним резонансом, является некорректным. При использовании линейно нарастающего входного задания частотная характеристика будет искажена демпфирующим воздействием кабельного перехода.
Рассмотрим процесс идентификаций методом Estimating ODE parameters. Для проведения иденти-
фикации модель (1) следует представить в матричной записи (здесь e t – белый шум):
x t Ax t Bu t Ke t ,
y t Cx t Du t e t ,
(2)
x 0 x0.
Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2010, № 4(68)
35
ПАРАМЕТРИЧЕСКАЯ ИДЕНТИФИКАЦИЯ ВЕНТИЛЬНОГО ЭЛЕКТРОПРИВОДА ...
Рис. 2. Амплитудно-частотная характеристика опытного образца ОПУ ТТИ (рис. 1, а)
Так как на рис. 3 наблюдается низкая зашумленность сигнала скорости, то влиянием шума измерений и шума объекта можно пренебречь. Матрицы исследуемой системы имеют следующий вид:
y i 1 1 2 2 T ; u U b 1_ отн b ;
1 Tf cm J1
A 0
0
0
ce L
m J1 1
m J2 0
0 c12 J1
0 c12 J 2
0
00
m J1
c12 J1
0 0 ;
m J2 c12 J 2
1 0
1 / L
0
0
0
0
0
0
0
B 0
0
Mc1 / J1 K 00
0 ;
0
0 0 0 Mc2 / J1
C 0 1 0 0 0, K 0 0 0 0 0, D 0 0 0 0 ;
x0 0 0 0 0 0 .
Для вычисления значений параметров записанной в матричном виде системы уравнений в рабочей
области Matlab задаются следующие команды:
model=idgrey(‘OPU’,parameter,’c’,0.001,0);model=pem(dannie,model);,
(3)
где “OPU” представляет собой функцию, задающую структуру модели:
function [A,B,C,D,K,x0] = OPU(par,T,CDmfile,aux)
A=[(-1/Tf), -par(1)/L, 0, 0, 0; 1.5*par(1)/par(2), -par(3)/ par(2), -par(4)/ par(2), par(3)/ par(2),
par(4)/par(2); 0, 1, 0, 0, 0;0, par(3)/par(5), par(4)/ par(5), - par(3)/ par(5), - par(4)/ par(5); 0, 0, 0, 1, 0];
B=[1/L 0 0 0;0 –par(6)/J1 -par(7) 0;0 0 0 0;0 0 0 –par(8)/J2; 0 0 0 0];
C=[0, 1, 0, 0, 0];
D=[0 0 0 0];
K=zeros(5,1);
x0=[0; 0; 0; 0; 0];
“parameter” – файл, содержащий априорные значения искомых параметров; dannie=iddata(s,u,Td) –
входные данные идентификации, сформированные из сигналов: s – массив данных скорости объекта, u –
массив данных управляющего напряжения, Td – шаг измерений в секундах.
36 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2010, № 4(68)
А.Г. Ильина
В результате выполнения команд (3) с учетом принятых допущений и при условии, что искомое решение существует, формируется модель заданной структуры с параметрами, близкими по значениям к указанным в файле parameter.m, и обладающая динамическими свойствами, сходными с реальным объектом. На рис. 3 представлены графики изменения угловой скорости при подаче входного меандра напряжения, полученные на опытном образце ОПУ ТТИ и при моделировании системы уравнений (1) где в качестве значений параметров используются полученные параметрической идентификацией методом Estimating ODE parameters.
Рис. 3. Сигналы изменения скорости разомкнутого вентильного электропривода при реакции на входной сигнал управления типа меандр: 1 – опытный образец ОПУ ТТИ (рис. 1, а), 2 – найденная математическая
модель
Как уже отмечалось, конечной целью получения адекватной математической модели объекта является синтез его системы управления. Продемонстрируем работу замкнутого контура скорости опытного образца ОПУ ТТИ (рис. 1, а), настроенного на технический оптимум (ТО), причем параметры ПИ регулятора контура рассчитаны на основании динамических свойств найденной с помощью идентификации модели. На рис. 4 представлены графики эталонной кривой отработки настроенной на ТО САР скорости входного ступенчатого воздействия и угловая скорость в САР скорости вращения оптической оси ОПУ ТТИ при реакции на тот же сигнал.
Рис. 4. Эталонная кривая (кривая 1) переходного процесса в замкнутом контуре скорости математической модели (1), настроенном на технический оптимум при отработке входного ступенчатого задающего сигнала, и аналогичный процесс в контуре скорости опытного образца ОПУ ТТИ (кривая 2)
Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2010, № 4(68)
37
ДИНАМИЧЕСКОЕ НАБЛЮДЕНИЕ НЕЛИНЕЙНЫХ ДВОИЧНЫХ...
Выводы
Проведена параметрическая идентификация, в результате которой получена математическая модель, удовлетворительно описывающая поведение реального объекта – опытного образца ОПУ ТТИ.
Применение разработанной модели при синтезе САР опытного образца ОПУ ТТИ приводит к гарантированному обеспечению желаемых качественных показателей проектируемой системы, а также значительно сокращает время на ее настройку.
Литература
1. Кротенко В.В., Толмачев В.А., Томасов В.С., Синицын В.А. Синтез микропроцессорной системы управления электропривода опорно-поворотного устройства // Изв. вузов. Приборостроение. – 2004. – Т. 47. – № 11. – С. 23–30.
2. Гурьянов А.В., Кротенко В.В. Параметрический синтез цифровых систем управления с широтноимпульсными преобразователями и эталонными моделями // Изв. вузов. Приборостроение. – 2006. – Т. 49. – № 3. – С. 17–24.
3. Микеров А.Г. Управляемые вентильные двигатели малой мощности: Учеб. пособие. – СПб: Издво ГЭТУ, 1997. – 64 с.
4. Takeshita T., Matsui N. Sensorless Brushless DC Motor Drive with EMF Constant Identifier // IECON-94. –1994. – V.1. – P. 8–13.
5. Ljung L. System Identification Toolbox 7 User’s Guide [Электронный ресурс]. – Режим доступа: http://www.mathworks.com/access/helpdesk/help/pdf_doc/ident/ident.pdf, своб.
6. Шрейнер Р.Т. Системы подчиненного регулирования электроприводов: учеб. пособие. – Екатеринбург: Изд-во ГОУ ВПО «Рос. гос. проф.-пед. ун-т», 2008. – 279 с.
Ильина Аглая Геннадьевна
– Санкт-Петербургский государственный университет информационных технологий, механики и оптики, научный сотрудник, dvanoska@mail.ru
38 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2010, № 4(68)