Научная статья на тему 'Локальная оптимизация дискретных электромеханических систем'

Локальная оптимизация дискретных электромеханических систем Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

CC BY
94
26
i Надоели баннеры? Вы всегда можете отключить рекламу.

Аннотация научной статьи по электротехнике, электронной технике, информационным технологиям, автор научной работы — Букреев В. Г., Богданов А. А., Соснин С. К., Чащин А. К.

Рассматривается решение задачи синтеза оптимального управления электромеханическими системами с переменной структурой. Предлагается алгоритм локальной минимизации ошибки регулирования электромеханического объекта на каждом интервале дискретности импульсного преобразователя. В качестве примера рассмотрена задача стабилизации скорости электродвигателя постоянного тока с широтно-импульсным модулятором.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по электротехнике, электронной технике, информационным технологиям , автор научной работы — Букреев В. Г., Богданов А. А., Соснин С. К., Чащин А. К.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

The article considers a solution to an optimal control synthesis problem for electromechanical systems with variable structure. It gives the algorithm of local error control minimisation for electromechanical object in each sampling interval of pulse converter. It considers the problem of speed stabilisation for a DC-motor with pulse-width modulator as an example.

Текст научной работы на тему «Локальная оптимизация дискретных электромеханических систем»

По результатам расчетов дефектных опор можно сделать вывод:

• При ширине раскрытия трещин 1.. .3 мм, коррозия арматуры 15.30 % и фактическая прочность бетона 194.220 кгс/см2, опоры не выдерживает нагрузки и их необходимо усилить.

• При ширине раскрытия трещин 1 мм, коррозия арматуры 1.3 % и фактическая прочность бето-

на более 300 кгс/см2, опоры выдерживают нагрузку и для них необходима заделка трещин, прогноз срока службы и наблюдение при дальнейшей эксплуатации линий.

• При разрушении защитного слоя оголенную, покрытую коррозией арматуру (вид 3), необходимо усилить дополнительными элементами с заделкой разрушенной части.

СПИСОК ЛИТЕРАТУРЫ

1. Демин Ю.В., Демина Р.Ю., Горелов В.П. Обеспечение долговечности электросетевых материалов и конструкций в агрессивных средах. - Кн. 1. Теоретические основы / Под ред. д.т.н. проф. В.П. Горелова. - Новосибирск: НГАВТ, 1998. - 209 с.

2. Демин Ю.В., Демина Р.Ю., Горелов В.П. Обеспечение долговечности электросетевых материалов и конструкций в агрессивных средах. - Кн. 2. Практические рекомендации. / Под ред. д.т.н. проф. В.П. Горелова. - Новосибирск: НГАВТ, 1998. - 190 с.

3. Барг И.Г. Воздушные линии электропередачи. Вопросы эксплуатации и надежности. - М.: Энергоатомиздат, 1985. - 248 с.

УДК 68-83:681.513.68

ЛОКАЛЬНАЯ ОПТИМИЗАЦИЯ ДИСКРЕТНЫХ ЭЛЕКТРОМЕХАНИЧЕСКИХ СИСТЕМ

В.Г. Букреев, А.А. Богданов, С.К. Соснин, А.К. Чащин

Томский политехнический университет E-mail: vbuk@yandex.ru

Рассматривается решение задачи синтеза оптимального управления электромеханическими системами с переменной структурой. Предлагается алгоритм локальной минимизации ошибки регулирования электромеханического объекта на каждом интервале дискретности импульсного преобразователя. В качестве примера рассмотрена задача стабилизации скорости электродвигателя постоянного тока с широтно-импульсным модулятором.

Введение

Процедура оптимизации дискретных электромеханических систем предполагает стационарность параметров объекта управления на достаточно длительном интервале его функционирования [1]. Однако часто возникают ситуации, когда на интервалах дискретизации управляющего сигнала параметры непрерывного электромеханического объекта принимают различные значения. Данным свойством обладают, например широтно-импульсные системы и электромеханические системы с источником энергии ограниченной мощности. Для таких систем целесообразно организовать субоптимальный регулятор, обеспечивающий экстремальность назначенного критерия качества в некоторой ограниченной области. К таким регуляторам можно отнести локально-оптимальный регулятор, минимизирующий, например, ошибку регулирования в каждый момент времени [2].

1. Постановка задачи

Пусть длительность функционирования электромеханической системы представлена точками /0,/ь...,4,... на интервалы дискретности, значения которых равны Т=4+1-4. где к - целое действительное число. Внутри каждого интервала [4, 4+1] поведение электромеханического объекта (ЭМО) управления описывается уравнениями

х(:) = Л1х(:)+Ь1 при : е х\1к, :к + тк ], хЦ) = Л2 х($)+Ь2 при : е \:к +т к, :к + Т ]. (1)

Здесь х(/) - «-мерный вектор переменных состояния, Ь1 и Ь2 - внешние воздействия, А1 и А2 -пхп - матрицы параметров. Управление объектом осуществляется за счет выбора последовательности значений т0,т1,.,тк, (интервал формирования управляющего воздействия). Примером подобных систем являются автоматические системы с широтно-импульсной модуляцией. Пусть 1^)=ст - желаемый скалярный выход объекта (ст - вектор-строка, т - символ транспонирования). Задачу управления можно поставить как задачу выбора такой последовательности значений т0,ть...,тк,...., при которых в установившемся режиме на каждом интервале времени Т величина г(0 незначительно отличается от некоторого заданного значения г, т.е. выполняется неравенство

) - г*| < е, (2)

где е - заданное малое число. Данное требование в другой формулировке можно рассматривать как задачу слежения за постоянным сигналом I на определенном интервале Т.

2. Решение уравнения динамического движения ЭМО

Предположим, что правые части ур. (1) не зависят от времени. Тогда решение данного уравнения на интервале времени [4,4+7] можно записать в виде

х(?к+тк^¿т) х(?к) +/1(тк), х(?к + Т) = и2(Т-тк) х(?к+тк)+^(Т-тк), где ^(/)=ехрА/ - матричная экспонента,

х(?) = |у1 (?) ЬЛ =(и, (?) -1) 42 Ьп 1 = 1, 2. (3)

* о

Здесь I - единичная матрица. Правая часть соотношения (3) записана в предположении, что матрица А невырожденная. В противном случае для вычисления интеграла нужно использовать разложение матрицы в ряд, т.е.

f(?) = Ь? + 4Ь?2¡2 +....

Отсюда можно получить дискретную модель объекта

х(к+1)= F(тк)х(к)+И(тк), к=0,1, 2 ... (4)

где

х(к )= х(?к ),

^ (т) = и 2 (Т-т) и ¡(г),

И (т) = и 2 (Т-т) т+аТ-т).

Таким образом, поставленная задача состоит в том, чтобы начиная с некоторого к-ого такта выполнялось неравенство (2).

Сформулировать достаточно строгие математически условия разрешимости данной задачи довольно трудно. Поэтому приведем некоторые общие соображения.

1. Значение г* должно быть достижимым, то есть должна существовать последовательность значений тт,...,^,...., при которой система из начального состояния переходит в окрестность значения г*. Здесь многое определяется величинами внешних воздействий Ь1 и Ь2. По-видимому, в большинстве практических случаев множество достижимых значений г=стх(к) можно оценить, если проинтегрировать систему (1) сначала при т=0, а затем т=Т.

2. Система (4) должен быть устойчивой, хотя бы для определенных значений параметра т в интервале [0,7].

3. Выход системы (1) должен быть регулируемым. Под этим понимается то, что при разных значениях т выход может увеличиваться или уменьшаться от такта к такту. Это свойство во многом определяется внешними воздействиями Ь1 и Ь2. Как видно из (4), для регулируемости системы необходимо, чтобы функция сТк(т) изменяла знак на интервале [0,7].

4. Предполагается, что величина интервала Тдостаточно мала, т.е. значительно меньше постоянных времени подсистем (1). Это необходимо, во-первых, для обоснования проводимой

ниже линеаризации уравнения (4) и, во-вторых, для обоснования отсутствия субгармонических колебаний внутри каждого интервала.

3. Локально-оптимальное управление

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

х(к+1)=Г0 х(к) + И0 +(Г1х(к) +И)тк, (5)

где

Е0=Е(0) = и2(Т), Е = и2(Т)ДА, ДА = А1 - А2,

Ио = f2(T )=(и 2 (Т) -1) А-%,

И 1= и2(Т) ДЬ, ДЬ = Ь - Ь2.

Из (5) следует, что ошибка слежения на к+1-ом такте равна

е(к +1) = г - сТх(к +1) = г - д0(х(к))- д1(х(к)) тк, где

д0(х) = сТ (Го х + И о) = сТи 2 (Т) х + сТИ о, д1(х) = сТ ^ х + И1) = сТи 2(Т)(ДАх + ДЬ).

Минимум модуля этой ошибки достигается при следующих условиях

о при р(х(к)) <о,

тк = \р(х(к)) при о < р(х(к)) < Т, (6)

Т при Т < р(х(к)).

где

Р(х) = (г*- Яо (х)) /Я1 (х). (7)

Таким образом, соотношение (6) определяет локально-оптимальный закон управления. Результаты моделирования поведения различных систем при приведенном выше алгоритме управления показывают, что система сначала при постоянных значениях т или Т переходит в окрестность значения г*, а затем начинается собственно процесс регулирования.

4. Локально-оптимальное управление последовательно включенными динамическими звеньями

Рассмотрим задачу управления двумя последовательными динамическими звеньями, между которыми имеется ключ, который в течение каждого такта включается в течение интервала времени [¡¡„¡+т] и выключается в течение интервала [4+т,4+7]. Поведение такой системы можно описать уравнением

" х4?)'

_ х2(? )_

I • А,.

I • А

х4?) + 81

х2(? ) 8 2

(8)

где х() и х2(/) - фазовые переменные первого и второго звена, g1и g2- внешние воздействия, действующие на звенья. Строку ст, определяющую требуемый выход

объекта, зададим в виде сТ=[0Т,с7], т.е. нас интересует выход второго звена. Параметр I, определяющий работу ключа, определяется следующим выражением:

1 [1 при tk <t < ^ +тк,

10 при tk +тк < t < tk + T. (9)

Таким образом, матрицы А1 и А2, фигурирующие в (1), отличаются только недиагональными блоками. При этом матрица А2 является блочно-диаго-нальной, и поэтому

~ип(Г) 0 " _ 0 и22(Г) _

и2 (Г )=

где и11(Т)=ехрА11Т и и22(7)=ехрА22Т. Кроме того,

А Л2

0

Л2,

Л12 0

К

(ии(Т) -1) Лп^1 (и22 (Т) - I) Л^g2

К = 0.

В результате выражения для д0и дь которые входят в (7), принимают вид

% (к?1= ¡^(к), где параметры определяются как

¿■Т =сТи22 (Г), ¡Т =СТ2 и22 (Г) Л21,

^0 = сТ (и22 (Г ) - I)

Следует отметить, что выражение для д0, которое входит в числитель выражения (7), зависит от переменных х2(/) - фазовых координат второго звена, а выражение для дь которое входит в знаменатель выражения (7), зависит от х() - фазовых координат первого звена. Кроме того, эти выражения не зависят от g1 - воздействия на первое звено.

5. Пример решения задачи стабилизации скорости

вращения двигателя постоянного тока

Электромеханический объект состоит из последовательно включенных силового фильтра и собственно двигателя. Уравнения работы силового фильтра имеют вид:

Ь/1/^) + Я,/^) + ^) = пм ^),

С/у/ ^) = у ^) - ИЦ). (10)

Уравнения движения двигателя имеют вид [1]: )+К I ^) + С^) = IVг ^), ^) = С^ ($) - М ($). (11)

Здесь Щ, Ь, С, $ - ток, напряжение, индуктивность, емкость и сопротивление фильтра, /(/) - сила тока в якоре, а>(() - угловая частота вращения двигателя, «,„(/) - входное напряжение фильтра (напряжение источника питания), М(/) - момент нагрузки двигателя, Ь и $ - индуктивность и сопротивление обмотки якоря, С - конструктивный коэффициент двигателя, / - момент инерции на валу двигателя.

Параметр I определяет работу ключа, согласно (9). В этой системе два внешних воздействия: напряжение

и;„(0, которое через фильтр поступает на двигатель на каждом такте Т в течение интервала [0,т], и нагрузка М(/), которая приложена к валу двигателя. Для того, чтобы скорость двигателя приняла заданное значение 6, необходимо, чтобы нагрузка М(/) была меньше некоторого значения. Ур. (10) и (11) можно представить в векторно-матричной форме, аналогичной (8),

'/^)

V/(t) )

о^)

0 " 0

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

где

аи = - Я/ /Ь/ .

^2 0 0 " у (t) V,

0 1а23 0 V/(t) 0

07} а" а34 ) V

0 | а43 0 6) _ _ то _

= -1/Ь1 > а21 = 1/Су , а23 = =-1/С/

а32 = 1/> аз3 = / Ь > а34 = -Се / Ь :

а 43 = С / 3,

V) = ип/ Ь >

т0 =-М / 3.

Моделирование проводилось при следующих числовых значениях (использовался двигатель ДК-1-2.3): 7=0,1 мс, ю*=10 рад/с, 1=2 мГн, С=4 мФ, Ь =2,24 мГн, $,=0,95 Ом, Сг=0,0231 Вс, /=0,00221 кгм2, и;=48 В, М=0,1 Н.м и нулевых начальных условиях. Оптимизация производилась по методике, изложенной в п. 4. ш, рад/с

12

ю

1 //2 / 3

I с

0.02

0.04

0.06

0.08

0.1

Рис. 1. Переходные процессы частоты вращения двигателя при разных моментах нагрузки: М: 1) 0,1; 2) 0,3; 3) 0,5 Нм

¡;, ¡, А; У;, В; ш, рад/с

2

... 3 11

"У-:..... 4

[ \ У-':

1

', с

0.02

0.04

0.06

0.08

0.1

Рис. 2. Переходные процессы в двигателе и силовом фильтре широтно-импульсного преобразователя: 1) ¡(), 2)уМ, 3) ¡(1), 4) о()

На рис. 1 представлены переходные процессы частоты вращения двигателя при варьировании нагрузки. На рис. 2 приведены переходные процессы для всех переменных состояния. Моделирование процессов произведено в системе МаНаЬ 7.

Выводы

Анализ результатов моделирования позволяет утверждать, что переходный процесс скорости дви-

гателя носит апериодический характер, и стабилизация заданной скорости осуществляется практически без перерегулирования примерно за 100 интервалов дискретности управляющего сигнала при варьировании момента нагрузки в широких пределах. Можно также заключить, что локально-оптимальный регулятор является в некотором смысле адаптивным по отношению к возмущениям и обеспечивает заданную точность регулирования.

СПИСОК ЛИТЕРАТУРЫ

1. Борцов Ю.А., Поляхов Н.Д., Путов В.В. Электромеханические системы с адаптивным и модальным управлением. - Л.: Энер-гоатомиздат, 1984. - 216 с.

2. Букреев В.Г., Параев Ю.И. Адаптивные регуляторы в дискретных системах управления сложными электромеханическими объектами. - Томск: Изд-во Том. ун-та, 2000. - 278 с.

3. Черноруцкий Г.С., Сибрин А.П., Жабреев В.С. Следящие системы автоматических манипуляторов. - М.: Наука, 1987. - 272 с.

УДК 621.312

КОММУТИРУЮЩИЕ СВОЙСТВА ОБМОТОК ЯКОРЕЙ КОЛЛЕКТОРНЫХ МАШИН ЭЛЕКТРОПРИВОДОВ НА НАЧАЛЬНОМ ЭТАПЕ КОММУТАЦИИ

С.И. Качин, Ю.С. Боровиков, О.С. Качин, Е.Н. Клыжко

Томский политехнический университет E-mail: borovikov@tpu.ru

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

Наличие непосредственной связи качества коммутационных процессов с конструкциями и параметрами якорных обмоток общеизвестно и отражено в ряде работ различных авторов [1-4].

Основными факторами, определяющими характер коммутационного процесса и зависящими от свойств обмотки якоря, являются средняя реактивная ЭДС секции, ее чувствительность к влиянию основного поля и поля дополнительных полюсов, степень взаимного демпфирования секций на протяжении времени их коммутации сопротивлением скользящего контакта, а также в течение времени горения дуги под сбегающим краем щетки.

Влияние указанных выше факторов на коммутационный процесс в той или иной степени исследовалось применительно к различным конструкциям якорных обмоток [1]. Наиболее подробно изучены коммутирующие свойства простых петлевых и волновых обмоток, которые широко применяются в машинах малой мощности.

Причем основное внимание при этом уделялось особенностям коммутации в течение промежутка времени, когда щетка непосредственно контактирует с набегающей и сбегающей ламелями исследуемой секции. Поэтому основополагающие реко-

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

Вместе с тем более детальное изучение особенностей процесса коммутации на начальном его этапе, когда коммутируемая секция электрически замкнута электрощеткой (так называемый процесс собственно коммутации [1]) позволяет выработать рекомендации по конструированию новых обмоток электрических машин, обладающих лучшими коммутирующими свойствами в сравнении с известными техническими решениями [2].

Частный случай процесса собственно коммутации при щеточном перекрытии, равном единице, был достаточно подробно рассмотрен в [1], где выявлено неоднозначное влияние величины индуктивности секции на коммутирующие свойства обмотки якоря в зависимости от соотношения других параметров коммутируемого контура и сформулированы условия достижения оптимальной величины индуктивности секции.

Значительно более сложными являются процессы коммутации при щеточном перекрытии,

i Надоели баннеры? Вы всегда можете отключить рекламу.