очень сложными [6], поскольку прошли путь интуитивного обобщения и согласования с опытом. Соответственным образом должны строиться формальные модели поиска решений. В данной работе показана роль выделения и использования стационарных инвариантов графов - деревьев кратчайших путей.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Иванов Д.С. Управление цепями поставок. - СПб.: Изд-во политехи. ун-та, 2009.
2. Розенберг И.Н., Беляков С.Л. Программные интеллектуальные оболочки геоинформаци-онных систем. - М.: Научный мир, 2010.
3. Берштейн Л.С., Беляков С.Л., Боженюк А.В. Использование нечетких темпоральных графов для моделирования в ГИС // Известия ЮФУ. Технические науки. - 2012. - № 1 (126). - С. 121-127.
4. Кристофидес Н. Теория графов. Алгоритмический подход. - М.: Мир, 1978.
5. Ерусалимский Я.М., Скороходов В.А. Графы с вентильной достижимостью. Марковские процессы и потоки в сетях // Известия вузов. Северо-Кавказский регион. Естесственные науки. - 2003. - № 3. - С. 3-5.
6. Люггер Д.Ф. Искусственный интеллект: стратегии и методы решения сложных проблем: Пер. с англ. - 4-е изд. - М.: Вильямс, 2005
Статью рекомендовал к опубликованию д.т.н., профессор В.Е. Золотовский.
Беляков Станислав Леонидович - Федеральное государственное автономное образовательное учреждение высшего профессионального образования «Южный федеральный университет»; e-mail: beliacov@yandex.ru; 347928, г. Таганрог, пер. Некрасовский, 44; тел.: 88634371743; кафедра прикладной информатики; д.т.н.; профессор.
Боженюк Александр Витальевич - e-mail: avb002@yandex.ru; кафедра прикладной информатики; д.т.н.; профессор.
Beliacov Stanislav Leonidovich - Federal State-Owned Autonomy Educational Establishment of Higher Vocational Education “Southern Federal University”; e-mail: beliacov@yandex.ru; 44, Nekrasovskiy, Taganrog, 347928, Russia; phone: +78634371743; the department of applied information science; dr. of eng. sc.; professor.
Bozhenyuk Alexandr Vitalievich - e-mail: avb002@yandex.ru; the department of applied information science; dr. of eng. sc.; professor.
УДК 681.513
В.Х. Пшихопов, А.Р. Гайдук, М.Ю. Медведев, В.Е. Беляев, Н.К. Полуянович,
Ю.П. Волощенко
ЭНЕРГОСБЕРЕГАЮЩЕЕ УПРАВЛЕНИЕ ТЯГОВЫМИ ПРИВОДАМИ ЭЛЕКТРОПОДВИЖНОГО СОСТАВА*
Рассматривается задача оптимизации движения электроподвижного состава по критерию энергопотребления с использованием комплексной математической модели, включающей уравнения механического движения и уравнения электромеханического и электрических преобразователей энергии. Приводится методика планирования движения состава вдоль, пути, заданного прямолинейными отрезками различных наклонов. Методика демонстрируется численным примером. На базе комплексной модели строятся алгоритмы управления, реализующие заданную оптимального траекторию движения. Погрешности модели учитываются в системе управления путем оценивания возмущений в процессе движения.
Электроподвижный состав; энергосбережение оценивание возмущений.
* Работа поддержана грантом РФФИ № 12-08-13112-офи_м_РЖД. 192
V.Kh. Pshikhopov, A.R. Gaiduk, М.Уи. Medvedev, V.E. Belyaev, N.K. Poluyanovich,
Yu.P. Voloshchenko
POWER SAVING CONTROL FOR ELECTRICAL TRAIN
This paper presents power saving control for electrical trains. The problem is solved on base of a complex mathematical model. The model describes mechanical movement of train, electromechanical drives, electrical convertors, and power grid. Method of planning of optimal train movement is proposed. The numerical example is presented. On base of the complex model the optimal control is designed. Adaptation of the control system is based on asymptotical robust estimators of disturbances. The results of computer simulation is presented.
Electrical train; estimation of disturbances; power saving.
Введение. Актуальность проблемы управления электроподвижным составом обусловлена интенсивным развитием силовой электроники, приводной техники, повышением скоростей движения и необходимостью энергосбережения [1]. Развитие технических средств позволяет повысить уровень автоматизации, в связи с чем в настоящее время имеется возможность использовать для управления тяговыми приводами электропоездов современные методы и алгоритмы управления [2-4], включая современные методы построения оптимального управления [5-8].
При построении системы автоматического управления движением, в том числе, электропоезда целесообразно применять иерархическую структуру управления, включающую уровень планирования движения, тактический уровень управления и уровень управления исполнительными механизмами [3].
С точки зрения энергосбережения эффект может быть достигнут на любом из перечисленных уровней системы управления. Так, на уровне исполнительных механизмов достижение энергосберегающего эффекта возможно как за счет применения более эффективных преобразователей, так и за счет управления, обеспечивающего соответствующие режимы работы. Однако наиболее существенный эффект может быть получен на уровне планирования движения при учете не только механических, но и электрических процессов.
В данной работе рассматривается решение задачи движения электропоезда при движении вдоль пути с заданным профилем.
Экономичное управление подвижным составом. Предположим, профиль пути из пункта А в пункт В известен и имеет вид приведенный на рис. 1. Он имеет возвышение на h 1 м и понижение на h2 м, но понижение меньше повышения, т.е.
h 2< h1. Предположим также, что известна средняя скорость V0, при которой состав (электровоз + вагон) преодолевает путь из А в В за заданное время с учетом времени на разгон и остановку поезда.
При движении по возвышенности с постоянной скоростью V, очевидно при
подъёме придется увеличивать тягу, т.е. затрачивать дополнительную энергию, а при спуске тормозить, тоже затрачивая дополнительную энергию. При этом исключается возможность преобразования высокой потенциальной энергии состава, полученной им на подъёме, в его кинетическую энергию.
Такое преобразование станет возможным, если только на вершине скорость состава будет меньше V0. Таким образом, если состав будет двигаться с переменной скоростью, появляется возможность экономии энергии, так как в этом случае указанную потенциальную энергию состава можно без потерь преобразовать в его кинетическую энергию движения. Однако, из-за того, что в данном случае h 2< h 1
после проезда возвышенности скорость состава окажется меньше требуемой V0 ,
что приведет к опозданию. С тем, чтобы исключить это, необходимо сообщить составу дополнительную кинетическую энергию. Наиболее целесообразным представляется это сделать на участке пути до вершины подъема.
Рис. 1. Графики профиля пути, скорости и тяги
Итак, возникает задача определения дополнительного ускорения и длины пути, на котором составу можно сообщить такое количество энергии, при котором не только его скорость в конце спуска будет равна требуемой V0, но и средняя скорость прохождения возвышенности будет тоже равна V.
Если тягу не менять на всем участке 5і - S3, то скорость состава будет изменяться так, как показано на рис. 1 жирной линией. Как видно на участке 5^ S2
она будет падать до значения V2 , а затем возрастать до значения У3. Если считать
массу состава сосредоточенной в точке, то кинетическая энергия в точке 8 2, из-за перехода её в потенциальную, будет равна Е2 = тУ02 - mgh 1= mV^. Следовательно, скорость состава равна
V? =Л/п2-?^ . (1)
Примем, что выполняется условие
gh1< к2. (2)
На участке 8 2 - 83 происходит обратное преобразование энергий, поэтому в точке 83 скорость
vз =\І V2 + g ф— 2). (3)
Или
Vз =а/У0?Г^ . (4)
Таким образом, скорость состава в точке 83, как и средняя скорость участке 81- 83, будут меньше требуемой У0.
В данной работе предлагается сохранять эти скорости равными У0, за счет придания дополнительного постоянного ускорения И составу на участке, предшествующем точке 32. Обозначим Зу расстояние, на котором действует это ускорение. Дополнительное возрастание скорости состава Удт показано на рис. 1. Значение этой величины в конце пути 3 у определяется известным равенством
^т =А/^г+2И3у - V (5)
Для того чтобы выполнялись указанные выше условия необходимо выбрать 3у и И так, чтобы, согласно рис. 1, выполнялись следующие условия:
а) приращение скорости состава за счет ускорения, должно быть равно потери скорости из-за разности высот Н1 и Н 2, т.е. согласно (5)
V2 + 2а-У0 = V -У3; (6)
б) площадь Ру (на рис. 1 она обозначена цифрой 1) должна быть равна пло-
щади Рп (на рис. 1 она обозначена цифрой 2):
Ру = Рп . (7)
Согласно рис. 1 площадь Ру определяется выражением
Ру=2 (>У02+2азу - Уо) . (8)
Из равенства (6) получаем 2а3у = 3Уо2 - 4УоУз + У32 и, подставляя в (8), при-
ходим к выражению
ру =
2 (у1 (2Уо - Уз)2 - Уо) 3у = ±(Уо - Уз)3у . (9)
Аналогично, согласно рис. 1, имеем
р. = ^(32 -31) + (3з -32).
Отсюда
Р. = 1 [Уо(2Бз-Я1 -32)-У2(3з -50-Уз3 -^У]. (1о)
Наконец, подставляя равенства (9) и (1о) в (7), получим
2(Уо -Уз)3у = 2[Уо(23^ 31 - 32) -У23 - 31) -Уз(3з -32)]. Отсюда
= [Уо(23^^ -Б2)-У^ -^)-Уз(3з -Я2)]. (11)
(Уо Уз)
С учетом выражения (11) из равенства (6) выводим
а = -^[зУо2-4УоУз +Уз2]. (12)
у
Полученные выражения (11) и (12) дают решение поставленной задачи. Приведем численный пример применения этих соотношений.
Уо3 = °3 °оа = _ _ = 20,84 м/с,
Пусть профиль пути соответствует рис. 1, причем Б1 = 2000 м, 82 = 2300 м, Б3 = 2500 м, h 1= 30 м, h 2 = 15 м, V 0= 20 м/с, g = 10 м / с2. Необходимо найти Б у (м), И (м/с2), при которых скорость в конце спуска, а также средняя скорость прохождения возвышенности будут равны V). Оценить расход энергии.
Согласно (1) и (4) У2 = ^202 -10 -30 = 10 м/с, У3 = ^202-10 -15 = 15,81 м/с, а У2 - У3 = -5,81 м/с. При этом по формуле (11) имеем Бу = 1393 м, а по формуле
(12) находим а = 0,06645 м/с2. Следовательно, если ускорение состава заканчивается в точке Б2, то точка Оа начала ускорения состава имеет координату Боа = Б2 - Бу = 2300 -1393 = 907 м.
Легко установить, что средние скорости на участках и времена их прохождения равны: Боа - Б1 ^ У01 = 21,65 м/с, гт = 50,427 с; Б1 - Б2 ^ У12 = 15,056 м/с,
г12 = 15,056 с; Б2 -Б3 ^У23 = 18,25 м/с, г23 = 10,959 с. Поэтому средняя скорость на этих участках равна
5, - ^ 1593
76,436
т.е. погрешность не превышает 5 %.
Таким образом, за счет движения состава с переменной скоростью и дополнительных затрат энергии на придание составу местного ускорения, обеспечено поддержание заданной средней скорости движения состава на всем пути, исключая участки разгона и торможения.
На придание составу местного ускорения дополнительно затрачено энергии в объеме
Ед = mgh2=150 - т.
Если бы скорость поддерживалась постоянной и при прохождении возвышенности, то на преодоление подъема на h 1= 30 м со скоростью У 0= 20 м/с была бы затрачена энергия
Е12 = mgh 1= 300 - т.
Кроме того, на торможение состава на спуске между 82 и 83 потребовалась бы энергия равная изменению потенциальной энергии на этом участке, т.е. Е23 = mg(^- h 2) = 10 - (30 - 15)m = 150m. Допустим, половина этой энергии (75m)
обеспечивается работой сил сопротивления движению, а вторая половина тоже (75m) обеспечивается активным торможением. Тогда полные затраты энергии на поддержание постоянной скорости движения состава будут равны 375m. Следовательно, при движении с переменной скоростью экономия составит
375m - 150m 225 . ,
Л =---------------=-----= 0,6,
375m 375
т.е. 60 %. С учетом конкретных условий, ряда не учтенных потерь и того, что состав имеет некоторую длину, эффект, конечно, будет ниже, но все же довольно существенным.
Система управления движением. Для реализации описанного выше движения состава можно применить адаптивную систему управления, схема которой приведена на рис. 2.
Рис. 2. Система управления движением
Предполагается, что тяговый электродвигатель постоянного тока (ТЭД), представленный на рис. 3, питается от управляемого преобразователя (УПр), обеспечивающего такие напряжения на двигателе, при которых ток якоря I изменяется плавно и не превышает допустимых значений.
и
Яя
ия
и„
]—7
Рис. 3. Схема тягового двигателя ТЛ-2К1
Уравнения двигателя ТЛ-2К1 с последовательным возбуждением с учетом управляемого преобразователя УПр имеют вид [1]:
(13)
(14)
КК = -(К + К + ^вт)Іа -ІрСеФ(/ц )+ Явт/ц + КупЛ
К ^= (/а - /ц ) Квт .
Кц = 2
ЛФ(/ц )
‘ц
, К = сеФ(1ц )1а1р / К,
(15)
где Яа = 0,0317 + 0,0356 = 0,0673 Ом, Ьа = 0,11 - индуктивность якорной обмотки; Кв = 0,025 Ом,; Квт = 0,005 Ом; Се = 2,31; функция Ф (/ц) Вхс/рад определена в работе [1]; ір = 3,826 - передаточное отношение редуктора; р = 2 -число пар полюсов ТЭД; ^ = 350 - число витков обмотки возбуждения; V -скорость состава м/с, Лд - напряжение на двигателе; Як = 0,625 м - радиус колес электровоза, Кт - сила тяги электровоза.
Система имеет два регулятора скорости ПИ и тяги ПИР, а также фильтр Ф, который предназначен для плавного ввода изменения заданной скорости движения состава. Согласно [9] это звено второго порядка с передаточной функцией
1
Т1Т2 р + (Т1 + Т2) р + 1
(16)
где Т1 = 2 с, Т2 = 1,5 с.
В схеме имеется переключатель П регуляторов. Так как профиль пути и масса состава m известны, то расчет величин Sy и tt а также дополнительной тяги
Fjd = a m выполняется перед стартом или во время движения до точки Soa. Уравнения движения состава имеют вид:
x = x3, (17)
m1x3 = -c(x1 - x2) -|u(x3 - x4) + F + W01(x3) + W31(x1), (18)
x2 = x4 , (19)
m2X4 = c(x1 - x2) + ц(x3 - x4) + W02(x4) + W31(x2), (20)
где m1 = 120 т и m2 = 71 т - массы электровоза и вагона;
W01(x3) = (1,9 + 0,01x3 + 0,0003x^)G1, Wa1(x1) = m1 g sin 5(x1), (21)
, . _ 8 + 0,1x4 + 0,0025x4 ґ ■ с-ґ ■
W02K) = (0,7 +----------17-75-±)Q2, Wd 2( x2) = m2 g sin 5( x2) , (22)
5 J,-|5(x¡)|, при движении вверх по склону,
' I |5(x¡ )|, придвижениивнизпо склону.
Здесь |5(x¡)| - острый угол наклона профиля пути в точке с координатой x¡, по горизонтальной оси; x3 , x4 - скорости электровоза и вагона соответственно в км/ч. G1 = m1g = 120 -10 = 1200 кН. Q2 = m2 g = 71-10 = 710 кН.
Значения коэффициента с и |і определяются интерполяцией по таблицам или по выражения приведенным в [9].
При разгоне состава и движении с постоянной скоростью (участок до точки Soa) переключатель П находится в положении 1. Система работает в режиме стабилизации скорости. На этом интервале фиксируется сила тяги F0, которая определяется весом состава и сопротивлением движению данного состава по горизонтальному пути с заданной скоростью V0. При достижении точки Soa переключатель П переводится в состояние 2 в результате чего в систему включается регулятор силы тяги, на его вход поступает сумма Fjd + F0. В этом положении переключатель П остается до прохождении точки S2 , после чего он переводится снова в положение 1. В результате в положении 2 переключателя П система управления на расстоянии Sy будет сообщать составу дополнительное ускорение tt, благодаря
чему к моменту прохождения составом точки S3 , его средняя скорость прохождения участка S0a - S3 и текущая скорость будут равны заданной V0. Это обеспечит
прибытие состава в пункт В в заданное время.
Расчет коэффициента передачи регулятора по скорости. В результате линеаризации в окрестности V0 = 20 м/с получены значения нелинейных функций:
W01 = 3977,76 Н; W02 = 1324,6 Н; F^ = 5302,4 Н; Ia 0 = 90 А; Ua 0 = 613 В.
Уравнения в отклонениях тягового двигателя имеют вид:
(p + 63,9) AF1, (p) = -1308,32AV (p) +16,7AIR (p) + 3339,12AUa (p),
(p + 3,33)AI^ (p) = 3,33AIa (p), ДFт (p) = 36,73AIa (p),
или
(0,3р2 + 64,2р + 63,45)А^т (р) = -1308,32АУ(р) + 3339,12АЦ, (р). (23)
Уравнения в отклонениях для состава:
(р2 + 5,5р)АУ (р) = 8,33 • 10-6 А^т (р) + / (р), (23)
где / (р) - изображение всех сил сопротивления движению, зависящих от профиля пути.
Объединив уравнения (23) и (24) и промоделировав замкнутую систему, найдем, что при ПИ-регуляторе в виде А и = 3500(У0 - У) переходной процесс имеет длительность около 8 с и не имеет перерегулирования.
Коэффициент передачи регулятора по тяге ПИР (см. рис. 2) подбирается в режиме управления по моменту.
Заключение. Предложенный метод позволяет решать задачу энергосберегающего управления движением состава, исходя из известного профиля пути. На основе предложенного подхода можно сформировать программную траекторию движения, которая может корректироваться в соответствии с текущим состоянием системы. Например, если в результате действия возмущений и погрешностей средняя скорость состава отличается от требуемой, то на следующие участи движения производится коррекция скорости движения.
Для повышения точности управления, учета нелинейностей, входящих в уравнения движения состава, и адаптации синтезируемой системы управления к неопределенным параметрам, возмущениям можно также использовать нелинейные методы управления движением и различные методы оценивания [10].
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Бахвалов Ю.А., Зарифьян А.А., Кашников В.Н., Колпахчьян П.Г., Плохов Е.М., Янов В.П. Моделирование электромеханической системы электровоза с асинхронным тяговым приводом. - М.: Транспорт, 2001.
2. Пшихопов В.Х., Медведев М.Ю. Структурный синтез автопилотов подвижных объектов с оцениванием возмущений // Информационно-измерительные и управляющие системы.
- 2006. - № 1. - С. 103-109.
3. Пшихопов В. Х. Позиционно-траекторное управление подвижными объектами. - Таганрог: Изд-во ТТИ ЮФУ, 2009. - 183 с.
4. Пшихопов В.Х., Медведев М.Ю. Блочный синтез робастных систем при ограничениях на управления и координаты состояния // Мехатроника, автоматизация и управление.
- 2011. - № 1. - С. 2-8.
5. Пшихопов В.Х. Оптимально по быстродействию траекторное управление электромеханическими манипуляционными роботами // Известия высших учебных заведений. Электромеханика. - 2007. - № 1. - С. 51-57.
6. Медведев М.Ю. Синтез замкнутых оптимальных по быстродействию управлений каскадными нелинейными динамическими системами с ограничениями на координаты // Мехатроника, автоматизация и управление. - 2009. - № 7. - С. 2-6.
7. Медведев М.Ю. Синтез субоптимальных управлений нелинейными многосвязными динамическими системами // Мехатроника, автоматизация и управление. - 2009. - № 12. - С. 2-8.
8. Гайдук А.Р., Капустян С.Г., Шаповалов И.О. Оптимальное перемещение тела интеллектуальным роботом // Мехатроника, автоматизация, управление. - 2009. - № 7. - С. 43-46.
9. Пудовиков О.Е. Управление длинносоставными тяжеловесными грузовыми поездами // Управление большими системами. - 2011. - Вып. 29. - С. 214-231.
10. Медведев М.Ю. Структура и алгоритмическое обеспечение нелинейного наблюдателя производных в условиях действия случайных шумов // Известия ЮФУ. Технические науки. - 2008. - № 12 (89). - С. 20-25.
Статью рекомендовал к опубликованию д.т.н., профессор Р. А. Нейдорф.
Пшихопов Вячеслав Хасанович - Федеральное государственное автономное образовательное учреждение высшего профессионального образования «Южный федеральный университет»; e-mail: pshichop@rambler.ru; 347928, г. Таганрог, пер. Некрасовский, 44; тел.: 88634371694; кафедра электротехники и мехатроники; зав. кафедрой; д.т.н.
Медведев Михаил Юрьевич - e-mail: medvmihal@gmail.com; кафедра электротехники и мехатроники; д.т.н.; профессор.
Беляев Виктор Егорович - e-mail: toe@tsure.ru; кафедра электротехники и мехатроники; к.т.н.; старший научный сотрудник.
Полуянович Николай Константинович - кафедра электротехники и мехатроники; к.т.н., доцент.
Волощенко Юрий Петрович - кафедра электротехники и мехатроники; к.т.н., доцент.
Гайдук Анатолий Романович - e-mail: gaiduk_2003@mail.ru; тел.: 88634371689; кафедра систем автоматического управления; д.т.н.; профессор.
Pshikhopov Vyacheslav Khasanjvich - Federal State-Owned Autonomy Educational Establishment of Higher Vocational Education “Southern Federal University”; e-mail: pshichop@rambler.ru; 44, Nekrasovskiy, Taganrog, 347928, Russia; phone: +78634371694; the department of electrical engineering and mechatronics; head the department; dr. of eng. sc.
Medvedev Мikhail Yur’evich - e-mail: medvmihal@gmail.com; the department of electrical engineering and mechatronics; dr. of eng. sc.; professor.
Belyaev Viktor Egorovich - e-mail: toe@tsure.ru; the department of electrical engineering and mechatronics; cand. of eng. sc.; senior scientist.
Poluyanovich Nikolay Kjnstantinovich - the department of electrical engineering and mechatronics; cand. of eng. sc.; associate professor.
Voloshchenko Yurij Petrovich - the department of electrical engineering and mechatronics; cand. of eng. sc.; associate professor.
Gaiduk Anatoly Romanovich - e-mail: gaiduk_2003@mail.ru; phone: +78634371689; the department of automatic control systems; dr. of eng. sc.; professor.
УДК 681.513
В.Х. Пшихопов, М.Ю. Медведев
АЛГОРИТМЫ ОЦЕНИВАНИЯ В СИСТЕМЕ УПРАВЛЕНИЯ АВТОНОМНОГО РОБОТИЗИРОВАННОГО ДИРИЖАБЛЯ*
Предлагаются алгоритмы оценивания возмущений в системе управления автономного дирижабля. Рассматриваемые алгоритмы предназначены для использования в рамках непрямого адаптивного позиционно-траекторного управления и реализуются непосредственно в процессе функционирования системы управления. Предлагаются алгоритмы оценивания аддитивных внешних возмущений, параметров дирижабля, а также его ускорений. Оценки строятся на базе нелинейной многосвязной модели дирижабля. В работе приводятся результаты численных экспериментов, подтверждающих работоспособность предложенных алгоритмов.
Управление дирижаблями; оценивание возмущений; оценивание параметров.
* Работа поддержана грантом Президента Российской Федерации № НШ-1557.2012.10 и грантом Евросоюза ЕР 7 «МАТ - многоуровневая транспортная система».