СИСТЕМЫ УПРАВЛЕНИЯ НА ТРАНСПОРТЕ
УДК 629.42.077
ПРИНЦИПЫ ПОСТРОЕНИЯ И МОДЕЛИ СИСТЕМЫ АВТОМАТИЧЕСКОГО УПРАВЛЕНИЯ ВТОРЫМ ЛОКОМОТИВОМ ПРИ ВИРТУАЛЬНОЙ СЦЕПКЕ
БАРАНОВ Леонид Аврамович, д-р техн. наук, профессор, профессор кафедры «Управление и защита информации»; e-mail: baranov.miit@gmail.com
БЕСТЕМЬЯНОВ Петр Филимонович, д-р техн. наук, профессор, профессор кафедры «Автоматика, телемеханика и связь на железнодорожном транспорте»; e-mail: bestpetr4@gmail.com БАЛАКИНА Екатерина Петровна, канд. техн. наук, доцент кафедры «Управление и защита информации»; e-mail balakina_e@list.ru
ПУДОВИКОВ Олег Евгеньевич, д-р техн. наук, доцент, профессор кафедры «Электропоезда и локомотивы»; e-mail: olegep@mail.ru
Российский университет транспорта (МИИТ), Москва
(М (М
о
(М
о
£ с
ш
га
0
S
<в
I-
о 0) L0
LQ
©
Рассмотрены модель системы автоматического управления локомотивом второго поезда при виртуальной сцепке и способы анализа качества управления. Управляемой величиной является рассчитываемое в процессе движения поездов расстояние между «хвостом» первого и «головой» второго поезда, обеспечивающее безопасность движения при экстренном торможении первого поезда. Это расстояние названо «длиной виртуальной сцепки». Введена формализация, позволяющая учесть отсутствие синхронизации при передаче информации о координате и скорости первого локомотива и такта работы цифровой системы управления вторым локомотивом. Рассмотрены погрешности измерений скорости и координат, связанные с запаздыванием при передаче информации, возможностью отказа от декодирования получаемой информации по радиоканалу приемником второго локомотива. При определении длины виртуальной сцепки рассмотрены различные способы ее верхней оценки, выбран способ, учитывающий измерение скорости движения и координаты обоих поездов, путь экстренного торможения первого поезда, возможность начала экстренного торможения первого поезда в начале цикла передачи информации, методические и инструментальные погрешности измерения. В качестве оператора, преобразующего рассогласования при выборе управления вторым локомотивом, рассмотрен цифровой пропорциональный интегральный дифференцирующий закон управления. Разработана модель системы автоматического управления, относящаяся к классу импульсных систем, особенностью которой является расчет управляющего воздействия в процессе движения поездов. Параметрический синтез системы автоматического управления рекомендуется проводить методами имитационного моделирования. Выбраны критерии качества управления.
Ключевые слова: виртуальная сцепка; экстренное торможение; служебное торможение; система автоматического управления; радиоканал; длина виртуальной сцепки; модель системы; имитационное моделирование; закон управления; погрешности измерений; вычисление длины виртуальной сцепки; критерии качества управления.
DOI: 10.20295/2412-9186-2022-8-04-377-388
▼ Введение
Одним из возможных направлений повышения пропускной способности железнодорожных линий является использование «виртуальной сцепки», при которой по радиоканалу первым ведущим поездом передается второму поезду информация о своей координате и скорости. Второй поезд выбирает свой режим движения таким образом, чтобы при экстренном торможении первого поезда не произошло аварийного столкновения. Эта технология
является развитием направления, связанного с движением соединенных поездов, имеющих физическое объединение, реализованное сцепкой локомотива второго поезда с последним вагоном первого.
13 сентября 1972 г. первый соединенный поезд, в СССР оборудованный телемеханической системой «Братск 1», передающей команды управления локомотиву второго поезда, был успешно проведен на участке Вихоревка — Чуна Восточно-Сибирской железной дороги.
Эта система разработана в Московском институте инженеров железнодорожного транспорта (ныне Российский университет транспорта — РУТ (МИИТ)), использовала для обмена информацией между локомотивами высокочастотный канал связи по контактной сети и была реализована еще на навесных полупроводниковых элементах [1]. В 1973 г. уже на микросхемах модернизированный вариант системы (Братск 1 м.), выполненный на интегральных микросхемах, прошел успешные испытания на участке Россошь — Лихая Юго-Восточной железной дороги. Это был первый опыт в нашей стране использования интегральной микросхемотехники на подвижном составе [2]. В дальнейшем технологией и системой управления соединенных поездов занимались ООО «АВП Технология», АО «ВНИИЖТ», АО «Московский тормозной завод ТРАНСМАШ», ООО «НПО САУТ».
Ограниченность длины станционных путей создает трудности при формировании и разъединении соединенных поездов, что приводит к занятию горловин и перегонов. Этих недостатков лишена технология виртуальной сцепки. Информация при виртуальной сцепке о координате и скорости ведущего поезда передается по радиоканалу. Работы в этом направлении ведутся как в нашей стране, так и за рубежом [3—7]. При внедрении виртуальной сцепки необходимо решить ряд задач, связанных с обеспечением безопасности движения. В частности, необходимо разработать систему контроля целостности состава [8, 9], контроля исправности рельсов, реализуемой в современных системах в рамках «контрольного режима» [10], использовать радиоканал с высоким уровнем достоверности и защиты при обмене цифровой информацией, реализовать систему точного определения координат движущихся поездов, разработать систему автоматического управления движением второго поезда с определенными характеристиками качества управления и парированием с нормированной высокой вероятностью опасных отказов. Должны быть решены технологические проблемы, определяемые согласованием числа поездов с виртуальной сцепкой с количеством станционных путей [3].
В рамках данной статьи р ешается круг вопр о -сов, связанных с принципами построения
системы автоматического управления движением второго из виртуально сцепленных поездов, разработки модели системы и методики анализа качества управления.
Целью системы управления второго поезда при виртуальной сцепке является обеспечение во время движения такого расстояния между «головой» второго и «хвостом» первого поезда, при котором экстренная остановка первого поезда не приводит к столкновению поездов. В дальнейшем для упрощения описания будем называть это расстояние длиной виртуальной сцепки — L. При создании системы управления необходимо выбрать оператора, преобразующий рассогласования А заданной и измеренной длиной виртуальной сцепки в сигнал управления режимами локомотива второго поезда. При тяговом приводе постоянного тока управление — это номер позиции управления N (тяговой, тормозной, выбег), при тяговом приводе на асинхронных двигателем — это сила тяги , торможения Кторм либо выбег.
Модель системы автоматического управления
Для описания модели управления используем следующие обозначения: пусть время t = (т + £)Т для первого локомотива. Здесь т = 0,1,2,...,0<£<1 (рис. 1, а), Т — шаг временной дискретизации при измерении координаты «головы» 51* и скорости У1 первого поезда. Аналогично для второго поезда t = (п + £)Т, п = 0,1,2,.,0<£<1 (рис. 1, б), 52 — координата «головы» и скорость У2 второго поезда.
Такая формализация позволяет учесть отсутствие синхронизации при измерении координаты пути и скорости поездов. Шаг временной дискретизации Т должен быть согласован с временем доставки информации по радиоканалу от первого локомотива ко второму. Так, если время передачи пакета незначительно меньше Т, то координата пути и скорость первого локомотива, соответствующая моменту времени тТ , получены вторым локомотивом в момент (т + 1)Т = (п + £ )Т и, следовательно, запаздывают на Т , что вносит дополнительную погрешность от запаздывания при определении рассогласования.
{п-4)Т (п-З)Т (п-2)Т (п-1 )Т пТ {п-1)Т {п-2)Т (п-З)Т
Рис. 1. Введение дискретного времени
Пусть I — длина первого поезда. Обозначим координату его «хвоста» как 51 = -1. Тогда в момент времени пТ в системе управления второго поезда известны значение координаты «хвоста» первого поезда, измеренной в момент времени (т - 1)Т = (п -2 + е1 )Т, и погрешность от запаздывания равна:
А [пТ ] = 51 [пТ ]-51 [(п -2 + е1 )Т ]. (1)
Эта величина соответствует ситуации, когда все последовательно передаваемые пакеты с первого локомотива принимаются вторым локомотивом и отсутствуют случаи отказа приемника от декодирования при обнаружении ошибки. Если первая передаваемая кодовая серия не принимается приемником, то
А [пТ ] = 51 [пТ ]-51 [(п -3 + е1 )Т ]. (2)
В общем случае, когда непринятыми могут быть подряд к кодовых серий:
А51>к [пТ] = [пТ] [(п - 2 - к + е1 )Т] . (3)
Аналогично погрешность от запаздывания при учете времени передачи информации составляет:
Ау [пТ ] = У1 [пТ ] -У1 [(п - 2 -к + е1 )Т ]. (4)
Максимальные значения погрешностей соответствуют е1 = 0 и к = ктах. Пусть Р — вероятность отказа от декодирования кодовых серий. События отказа от декодирования статистически независимы. Тогда вероятность того, что будет принята (к +1) -я кодовая серия,
составляет Рк (1 - Р). Величина ктах , следовательно, и максимальная оценка погрешности от запаздывания могут быть приняты с определенной вероятностью. Данное рассмотрение соответствует ситуации, когда вероятность трансформации кодовой серии пренебрежимо мала. Современные методы помехозащи-щенного кодирования позволяют решить эту задачу.
Кроме погрешностей А51 и Ау, определяемых запаздыванием, связанных с передачей информации с первого локомотива, имеют место методические и инструментальные погрешности измерения скорости и координаты. Пусть А51а, А52а, АУ1а, АУ2а — соответственно значения этих погрешностей измерения скорости и координаты первого и второго поездов. Погрешности Ауа и А51а суммируются с погрешностями А51 и Ау . В частности, методическая погрешность измерения скорости определяется квантованием по уровню и погрешностью от запаздывания [11, 12].
Последняя объясняется тем, что измеряется средняя скорость поезда за время измерения Ту < Т [11, 12].
Анализ величины погрешностей измерительных тактов особенно существенен, так как они определяют показатели качества системы управления и не парируются обратной связью. Модель измерительных тактов системы управления длиной виртуальной сцепки приведена на рис. 2. Так как информация о координате и скорости первого поезда используется в системе управления вторым поездом с запаздыванием, определяемом временем передачи, в модели включены соответствующие звенья с запаздыванием на (к + 2 - е1 )Т.
Рис. 2. Модель измерительных тактов
В измерительных тактах скорости звенья с
1 - е" рТ
передаточными функциями ——— определяют усреднение за время Ту при измерении скорости [11, 12], импульсный элемент моделирует временную дискретизацию — преобразование функции непрерывного аргумента t в функцию дискретного аргумента пТ, п = 0,1,2,... .Аналогичны тракты измерения координат. Отдельно рассмотрим составляющие погрешности измерения скорости. Сумму методической погрешности квантования по уровню и инструментальными погрешностями обозначим как АК|МИ или АКМаи . Вместе с
1 - е"рТу
тем наличие звена-, частотная характера
ристика которого
1 - e
wTv
-jwTv sm—- . Tv 2 „-]wу
jwT-
wT-2
e
гово-
T
рит о запаздывании на в тракте измерения
и, следовательно, дополнительной погрешности от запаздывания. Эту погрешность
Т
можно оценить сверху [11, 12] как атах, где атах — максимальное ускорение или замед-
T
ление поезда. Тогда AK* < AV,™ + omav — и
i a i a max
AV* < AV™ + . При аналитических оцен-
ках длины виртуальной сцепки из условий обеспечения безопасности будем считать, что
сумма А^Ми + атах у и Г2МГ +равны
соответственно АР^ и АР^. При имитационном моделировании функционирования
системы отдельно моделируется звено
1-e
- рт-
PT
а составляющие AV1™ и AV2™ задаются своими законами распределения плотности вероятности.
В качестве оператора управления выберем цифровой ПИД (пропорциональный, интегральный, дифференцирующий) закон управления. Управление и в этом случае определяется выражением:
п
и = к1 А [пТ] + к2 ^А [тТ] + к3 А [пТ] -
т=0
-А[(п-1)Т], (5)
где к1, к2, к3 — коэффициенты;
А[пТ ] — величина рассогласования; А[пТ ] - А[(п - 1)Т ] — обратная первая разность рассогласования;
п
^А[тТ ] — эквивалент интегрирования
т=0
рассогласования А[тТ ]. Величина рассогласования:
A[nT] = L[nT]-La3U [nT],
(6)
где Ь [пТ] — заданная в момент пТ требуемая длина виртуальной сцепки; Ьизм [пТ ] — измеренная на втором локомотиве в момент времени пТ длина виртуальной сцепки.
LH3M [nT] = ^2изм [nT] - StoM [nT]
(7)
где 52изм [пТ] — известная в момент пТ координата «головы» второго локомотива; 51изм [пТ ] — полученное значение устройством управления второго локомотива из канала в момент времени пТ значения координаты «хвоста» первого поезда, измеренное на первом локомотиве в момент времени (п -к -1 -е1 )Т .
Способы преобразования величины и при дискретном управлении тягой и торможением в номер N позиции управления является самостоятельной задачей. В этом случае появляется дополнительное требование — уменьшение числа переключений позиций управления тяговым приводом при выполнении заданных показателей качества функционирования системы автоматического управления.
Выбор длины виртуальной сцепки
Длина виртуальной сцепки Ь[пТ] — выбирается из условий обеспечения безопасности движения при экстренном торможении первого поезда. Для идеального случая, когда отсутствуют погрешности измерений скорости и координаты поездов:
L [nT] = S* [nT] -1 - S2 [nT] +
+S
1эторм
(V S2CJtopm V ),
(8)
где S1эторм (V) — путь экстренного торможения первого поезда, движущегося со скоростью V1; S2cmopM (V2) — путь служебного торможения второго поезда, движущегося со скоростью V2.
Учитывая наличие погрешностей от запаздывания, методических, аппаратных погрешностей при измерении скоростей и координат поездов, разброс в задании путей экстренного и служебного торможения, величина L выбирается для «худших условий» — максимальных величин погрешностей.
Максимальная величина погрешностей запаздывания AS1 к [nT ] определяется из (3) при е1 = 0 и величине к, при которой Рк (1 - P) < 10-8 согласно SIL 4.
Оттуда:
к >■
8 + lg (1 - Р) ' log Р '
(9)
Оценку сверху погрешности от запаздывания вычислим как путь, пройденный первым поездом при максимальной скорости за максимальное время запаздывания:
AST =(к + 2)TVm
(10)
Тогда оценка максимальной погрешности координаты первого поезда, используемой для управления вторым поездом, определяется выражением:
ASmax = ASmах + AS™* < ASmax + (к + 2)TVm
(11)
где AS1max — максимум модуля погрешности AS1a.
Оценка сверху погрешности измерения координаты второго поезда определяется только методическими и инструментальными составляющими:
Л2тах = ^гТ*,
где — максимум модуля погрешности Л^.
Тогда:
Ь ^торм ^ах Мэ (У )■+ ЛТ +
+Л^2таХ +Л1 + (к + 2)ТУ тах, (12)
где Л! — погрешность в задании длины первого поезда.
Для верхней оценки величины Ь можно не учитывать путь экстренного торможения первого поезда и принять:
Ь > ^торм (У2тах ) + ЛТ + Л^Г + Л!. (13) ЛУ2
Величина £2Торм (У2тах) может определяться следующим образом:
^2торм (У2тах )
У2
' От
2
(14)
где ^тах — максимальное замедление второго поезда при служебном торможении.
В выражениях (12) и (13) погрешность от запаздывания при передаче информации о координате первого поезда учитывается оценкой сверху, что соответствует управлению вторым поездом с запасом по условиям безопасности. Вместе с тем эту погрешность нужно учитывать при имитационном моделировании всей системы управления для анализа изменений длины виртуальной сцепки в процессе движения.
Полученная верхняя оценка длины виртуальной сцепки не учитывает погрешностей, возникающих в системе автоматического управления вторым поездом, анализ которых проводится методами имитационного моделирования. После разработки адекватной математической модели системы, верифицированного программного обеспечения модели проводится параметрический синтез системы. Задачей параметрического синтеза является выбор значений
коэффициентов к1,к2,к3 закона управления, при которых обеспечиваются заданные показатели качества управления. Учитывая, что при функционировании системы автоматического управления длины виртуальной сцепки имеются динамические и статические погрешности, требуемая величина Ь должна быть увеличена на максимальное значение суммарной погрешности системы.
Описанный выше способ определения минимальной длины виртуальной сцепки и управления вторым поездом по этой длине не единственен. Рассмотрим вычисления длины Ь с учетом известной скорости У2 второго поезда и управления без учета пути экстренного торможения первого поезда. В этом случае на каждом такте управления может быть вычислена величина Ь по выражению (13) после замены на ^Торм |>2изм К] + АУ2тах ] , где
^2торм (Утах )
Г — максимальная величина погрешности измерения скорости второго поезда, и 52 вычисляется по (14) после замены У21 на
2 тах
У2изм [пТ] + ЛУ2тах ] . При этом допустимая величина Ь уменьшается. Следующий способ уменьшения длины виртуальной сцепки является учет пути экстренного торможения первого поезда, движущегося со скоростью У1 . Длина пути возможного экстренного торможения первого поезда при У1 изм [пТ] - ЛУ^ вычисляется на борту этого поезда и передается в общем пакете вместе с координатой и скоростью второму поезду, где и учитывается при вычислении Ь согласно (12). При этом расчет первого слагаемого в (12) зависит от выбора одного из двух способов, описанных выше.
Вместе с тем расчет пути экстренного торможения первого поезда должен учитывать возможность уменьшения его скорости за время (к + 2)Т передачи информации второму поезду (рассматриваем длительность передачи для худшего случая). За время (к + 2)Т скорость первого поезда перед экстренным торможением может уменьшится с У до У1 - gmax (к + 2)Т , где &тах — максимальное замедление поезда. Следовательно, путь экстренного торможения первого поезда, движущего со скоростью У1 , для гарантии безопасности должен быть уменьшен с учетом возможного уменьшения его скорости и определяться выражением:
^ [У У -
Ауут ах
(к + 2)Т ] =
тах (к + 2)Т]2
2 Я
(15)
где ятах — максимальное замедление, которое
можно выбрать равное яэ.
Приведем анализ различных способов вычисления Ь — длины виртуальной сцепки. Для дальнейшего рассмотрения введем обозначения этих способов:
1-й способ — для расчета Ь используется путь служебного торможения второго поезда от скорости У2тах , путь экстренного торможения первого поезда принимается равным 0;
2-й способ — для расчета Ь используется путь служебного торможения второго поезда с учетом начала торможения от скорости У2изм [пТ] + АУ2аах , путь экстренного торможения первого поезда принимается равным 0;
3-й способ — путь служебного торможения второго поезда определяется от скорости начала торможения У2тах (так же как и в первом способе), учитывается путь экстренного торможения от скорости уизм [пТ ] - Ау^ах первого поезда;
4-й способ — путь служебного торможения 2-го поезда определяется от скорости У2изм [пТ] + АУ™ах (так же как и во втором способе), учитывается путь экстренного торможения от скорости У1изм [пТ ] - Ау^ах первого поезда.
Выражение для определения длины виртуальной сцепки для рассмотренных четырех способов сведены в табл. 1.
Из качественных соображений очевидно, что минимальную из рассмотренных способов определения длину виртуальной сцепки обеспечивает 4-й способ при обеспечении условий безопасности одинаковых максимальных погрешностях трактов измерения и методах передачи информации.
Для получения численных результатов зададим следующие исходные данные: А^™* = 10 м;
АС* = 10м; ДУГ = 1 км/ч; А У™ = 1 км/ч, к = 10,
Т = 140 мС; Ясл = 0,6 м/с2; & = 0,7м/с2; ^ = 0,7м/с2; У2тах = 80 км/ч. В выбранных данных Т = 140 мС при передаче сообщения длиной 24 байта по радиоканалу на частоте 160 МГц. При использовании радиоканала на частоте 2130 МГц время передачи 24 информационных байта составляет 70 мС. Результаты расчета сведены в табл. 2.
Обратим внимание, что допустимая длина виртуальной сцепки для второго способа, когда У2изм = У2тах = 80 км/ч превышает соответствующую величину для первого способа. Это вызвано тем, что У2изм = 80 км/ч известно с погрешностью, равной 1 км/ч. При выборе L в соответствии с первым и третьим способами величина У2тах известна точно.
Анализ результатов третьего способа определения L говорит о том, что его использование
Таблица 1. Выражения для расчета длины виртуальной сцепки
Способ £ > 0
1 У 2 Ь > 2тах + Д™ах + Д2таах + А1 2Ясл (16)
2 Ь > (У2изм [пТ] + ДУ2 а ) + дтах + дтах +Д 1а 2а 2 Ясл (17)
У2
Ь > ^пах + д тах + а тах + д -
1а 2 а
2Ясл
[Уизм [пТ]-АУаах -атах (к + 2)Т]2 2 £э
(18)
(У [пТ] + АУтах)2 [Уизм [пТ]-АУ"ах -атах (к + 2)Т]
Ь > 2изм [ 1 2а ' + Дт ах + Дтах + д [ 1изм Ь 1 1а т ах V ) ]
2 Я с
2ёэ
(19)
3
Таблица 2. Допустимая длина виртуальной сцепки
Способ 1
^2тах = 80 км/ч ^2изм = 80 км/ч ^2изм = 60 км/ч уп = 40 км/ч 2изм
1 441,4 м — — —
2 — 451,75 м 269,13 м 108,11 м
3 140,05 м пРи ^1изм = 80 км/ч 150,52 м пРи ^1изм = 80 км/ч 281,11 м пРи ^1изм = 60 км/ч 377,88 м пРи ^1изм = 40 км/ч
4 — 150,4 м пРи ^1изм = 80 км/ч 108,84 м пРи ^1изм = 60 км/ч 44,59 м пРи ^1изм = 40 км/ч
приводит к увеличению допустимой величины длины виртуальной сцепки при низких скоростях движения за счет уменьшения пути экстренного торможения первого поезда при завышенной верхней оценке пути служебного торможения второго поезда. Перспективным для использования в автоматической системе управления вторым поездом можно считать четвертый способ, учитывающий значения измеренных координат и скоростей первого и второго поездов, при выполнении:
ГУ +ЛУтах I2
1 2из^ 2 ^ + Лт ах + Л2таах + Л! >
1а 2а
2 gсл
>[У^изм -ЛУ^ -gэ (к + 2)Т]2
> ^э .
(20)
где Ь0 =Л5тах + ЛСх +Л!;
ЛУтах + е
глу1а Т ¿»тах
(к + 2)Т
ЛУ тах
Т = У _+ У
М ' 2изм ' ' 1изм
Величина Ь0 определяется только максимальными погрешностями измерений и известна заранее.
Величина Ь1 зависит как от максимальных погрешностей измерений, так и от измеренных значений скорости и возможного максимального уменьшения скорости первого поезда во время передачи информации. Ь1 можно оценить сверху, приняв У2изм = У2тах ,У1изм = ^тах .
Тогда верхняя оценка этой составляющей
ЛУтах Ь = У + У
М ' 2тах ' ' 1тах
+ gmax (к + 2)Т
Условие (20) выполняется, так как в противном случае скорость второго поезда настолько мала, что его путь служебного торможения меньше пути экстренного торможения первого поезда. Следовательно, система управления может увеличить скорость второго поезда и приблизить его к первому. При равенстве скоростей движения величина Ь приведена в табл. 2. Результаты расчета количественно показывают, как изменяется длина виртуальной сцепки при учете пути экстренного торможения первого поезда.
Проанализируем дополнительно выражение (19) для четвертого способа вычисления длины виртуальной сцепки. После несложных преобразований получаем:
У2 у 2
Ь > 2изм 1изм
^сл ^э
+ Ьо + Ь
(21)
и
у 2 У 2
Ь > 2изм 1изм
^сл 2 gэ
+ Ьо + Ь .
(22)
Сумму величин Ь0 + Ь1 по существу можно определить как длину защитного промежутка при сближении поездов. Для приведенных выше исходных данных Ь0 + Ь* = 80,37 м.
Таким образом, длина виртуальной сцепки в соответствии с выражением (22) определяется как сумма, разность путей служебного торможения второго и экстренного торможения первого поездов при измеренных значениях их скоростей и длины защитного промежутка.
3. Блок-схема модели
Блок-схема модели для анализа функционирования системы управления вторым поездом в условиях применения виртуальной сцепки
Рис 3. Блок-схема модели
при заданном управлении и1 ф первого поезда приведена на рис. 3.
Модель поезда задается системами известных дифференциальных уравнений [13]. В частности, поезд может рассматриваться как материальная точка, «нерастяжимая нить» [14]. Возможно, использование более сложных моделей, если исследователя интересует продольная динамика поезда [15—19]. Модели измерительных трактов координат и скоростей поездов с учетом запаздывания при передаче информации соответствуют описанным выше и приведенным на рис 2. Закон управления реализуется блоком ПИД, моделирующим преобразованием рассогласование Д[пТ] в управление и [пТ]. Идеальный импульсный элемент ИЭт преобразует функцию дискретного аргумента и [пТ] в решетчатую функцию — последовательность 5-функций с весами и [пТ] . Решетчатая функция экстраполятором нулевого порядка ЭНП, передаточная функция
которого 1
Р
преобразуется в ступенча-
тый сигнал управления и, поступающий на вход объекта — вход управления локомотивом. Задаваемая системе управления величина длины виртуальной сцепки вычисляется в «Блоке задания» при известных измеренных
величинах Уlу3иЗifгrГ}, Уизм[(п к 1-Ё1 )Т] и
заданных Ь0 и ¿1 в соответствии с выражением (22). Измеренное в момент времени значение длины виртуальной сцепки Ьизм [пТ вычисляется как разность 51изм [(п -к -1 -е)Т и 52изм [пТ]. Величина рассогласования А[пТ вычисляется как разность Ь [пТ] - Ьизм [пТ].
Качество управления системы характеризуют:
— Ьг [пТ] = 51 [пТ] - 52 [пТ] — зависимость действительной длины виртуальной сцепки от времени получаемой при различных временах хода первого поезда по заданному участку, где 51 [пТ] — зависимость координаты пути «хвоста» первого поезда от времени;
- РТ
е
52 [пТ] — зависимость координаты «головы» второго поезда от времени;
— Л [пТ] = Ь [пТ] - Ьg [пТ] — зависимость погрешности реализации заданной длины виртуальной сцепки от действительной, получаемой при различных временах хода первого поезда по заданному участку;
— максимальные, средние и среднеквадра-тические оценки Ь [пТ] и Л[пТ] по времени при движении по заданному участку;
— усредненные по времени и по множеству средние и среднеквадратические оценки Ьg [пТ] и Л[пТ] при движении по различным участкам;
— число переключений управления при релейно-контакторном управлении привода локомотива;
— обязательное отсутствие случаев Л[пТ ] < 0, что определяет удовлетворение условий безопасности.
Таким образом, разработанная модель позволяет создать алгоритмическое и программное обеспечение для проведения имитационных экспериментов функционирования системы автоматического управления вторым локомотивом при виртуальной сцепке. По результатам имитационных экспериментов при выбранных показателях качества управления проводится параметрический синтез системы и вычисляются характеристики ее функционирования.
Заключение
Повышение использования пропускной способности участков железных дорог возможно при использовании технологии виртуальной сцепки. Реализация этой технологии связана с решением сложных технических и технологических задач.
Существенной проблемой при внедрении технологии виртуальной сцепки является, в частности, создание системы автоматического управления движением второго локомотива при передаче информации по радиоканалу по скорости и координате первого локомотива.
Предложенные принципы построения системы автоматического управления движением второго локомотива при виртуальной сцепке, разработка математической модели
системы, выбор показателей качества управления являются фундаментальным для решения задачи параметрического синтеза этой системы, определения характеристик ее функционирования.
Описанные в данной статье результаты исследований являются определенным этапом, предшествующим разработке алгоритмического программного и информационного обеспечения цифрового двойника технологии виртуальной сцепки. ▲
Библиографический список
1. Баранов Л. А. Система автоматического и телемеханического управления электроподвижным составом / Л. А. Баранов, Е. В. Ерофеев, В. И. Астрахан и др.; под ред. Л. А. Баранова. — М.: Транспорт, 1984. — 311 с.
2. Баранов Л. А. 35 лет кафедре «Автоматика и телемеханика» / Л. А. Баранов // Автоматическое управление технологическими процессами на транспорте. Юбилейный сборник научных трудов. — М.: МИИТ, 1996. — Вып. 892. — С. 3-8.
3. Бушуев С. В. Повышение пропускной способности участка железной дороги с применением технологии виртуальной сцепки / С. В. Бушуев, К. В. Гундырев, Н. С. Голочалов // Автоматика на транспорте. — 2021. — № 1. — С. 7-20. — DOI: 10.20295/2412-9186-2021 -7-1-7-20.
4. Климова Е. Е. Технология виртуальной сцепки поездов как инструмент, повышения пропускной и провозной способности линии / Е. Е. Климова, Л. Е. Пилинушка,
B. С. Рябов // Транспортная инфраструктура сибирского региона: материалы десятой Международной научно-практической конференции. — 2019. — С. 60-64.
5. Flammini F. Towards Railway Virtual Coupling / F. Flammini, S. Marrone, R. Nardone et al. // International Conference of Electrical Systems for Aircraft, Raiway, Ship Propulsion and Road Vehicles and International Transportation Electrification Conference. — 2019.
6. Mitchell I. Train Convoys or Virtual Coupling / I. Mitchell, E. Goddard, F. Montes // IRSE news. Institution of railway signal engineers. — 2016. — Iss. 219.
7. Goikoetxea J. Roadmap Towards the Wireless Virtual Coupling of Trains / J. Goikoetxea // Springer International Publishing Switzerland. — 2016. — Pp. 3-9.
8. Никитин А. Б. Результаты исследования технических средств контроля целостности поездов / А. Б. Никитин, И. В. Кушпиль // Автоматика на транспорте. — 2020. — № 4. C. 411-434. — DOI: 10.20295/2412-9186-2020-6-4411-434.
9. Quante F. Контроль полносоставности грузовых поездов / F. Quante // Железные дороги мира. — 2005. — № 2. — С. 43-48.
10. Bestemfyanov P. F. Energy efficient algorithms assessment of the Rail circuit operation / P. F. Bestemfyanov // Russian Electrical Engineering. — 2017. — Vol. 88. — № 9. —
C. 557-562.
11. Баранов Л. А. Квантование по уровню и временная дискретизация в цифровых системах управления / Л. А. Баранов. — М.: Энергоатомиздат, 1990. С. 306.
12. Баранов Л. А. Микропроцессорные системы автоведения поездов / Л. А. Баранов, Я. М. Головичер, Е. В. Ерофеев и др.; под ред. Л. А. Баранова. — М.: Транспорт, 1990. — С. 272.
13. Розенфельд В. Е. Теория электрической тяги / В. Е. Розен-фельд, И. П. Исаев, Н. Н Сидоров; под ред. И. П. Исаева. — М.: Транспорт, 1995. — С. 294.
14. Баранов Л. А. Оптимизация управления движением поездов / Л. А. Баранов, Е. В. Ерофеев, И. С. Милешин и др.; под ред. Л. А. Баранова. — М.: МИИТ, 2011. — С. 164.
15. Лисицын А. Л. Нестационарные режимы тяги (Тяговое обеспечение перевозочного процесса) / А. Л. Лисицын, Л. А. Мугинштейн. — М.: Интекст, 1996. — 159 с.
16. Gao Gj. of longitudinal forces of coupler devices in emergency braking process for heavy haul trains / Gj. Gao, W. Chen, J. Zhang et al. // J. Cent. South Univ. — 2017. — № 24. Pp. 2449-2457. — DOI: https://doi.org/10.1007/ s11771-017-3656-9
17. Pugi L. Modelling the longitudinal dynamics of long freight trains during the braking phase. / L. Pugi, D. Fioravanti, A. Rindi // 12th IFToMM World Congress. Besancon, France. — 2007. Pp. 1-6.
18. Cole C. Longitudinal train dynamics / C. Cole, S. Iwnicki (ed.) // Handbook of railway vehicle dynamics. Taylor & Francis. — London, 2006. — Pp. 239-278.
19. Wu Q. A review of dynamics modelling of friction draft gear / Q. Wu, C. Cole, S. Luo et al. // Veh Syst Dyn. — 2014. № 52(6). Pp. 733-758.
Благодарность РФФИ
Исследование выполнено при финансовой поддержке РФФИ, НТУ «Сириус», ОАО «РЖД» и образовательного фонда «Талант и успех» в рамках научного проекта № 20-37-51001.
Acknowledgments. The reported study was funded by RFBR, Sirius University of Science and Technology, JSC Russian Railways and Educational Fund «Talent and success», project number 20-37-51001.
TRANSPORT AUTOMATION RESEARCH, 2022, Vol. 8, No. 4, рр. 377-388 DOI: 1020295/2412-9186-2022-8-04-377-388
Design Principles and Models of Automatic Control System for Second Locomotive at Virtual Assessment
Information about authors
Baranov L. A., Doctor in Engineering, Professor, Professor of the Department of "Information Management and Protection". E-mail: baranov.miit@gmail.com
Bestemyanov P. F., Doctor in Engineering, Professor, Professor of the Department of "Automation, Remote Control and Communication in Railway Transport". E-mail: bestpetr4@gmail.com
Balakina E. P., PhD in Engineering, Associate Professor of the Department of "Information Management and Protection". E-mail balakina_e@list.ru Pudovikov O. E., Doctor in Engineering, Associate Professor, Professor of the Department of "Electric Trains and Locomotives". E-mail: olegep@mail.ru
Russian University of Transport (MIIT), Moscow
Abstract: Automatic control system model of second train locomotive at virtual assessment as well as control quality analysis approaches have been considered. The controlled value represents being calculated in train traffic process the distance between"tail"of a first train and "head"of a second one, providing traffic safety at emergency braking of a first train. This distance is called "virtual replacer length". Formalization has been introduced allowing to account for synchronization absence at transfer of information on the coordinate and speed of a first locomotive and the control digital system work cycle of a second one. Speed and coordinate measurement errors have been considered which're connected with lagging at information transfer and possibility of failing in decoding of received by radio channel information by second locomotive receiver. At virtual length assessment, various ways of its upper assessment have been consid-
ered, the way has been chosen, taking into account motion speeds and the coordinates of both trains, emergency braking path of a first train, the possibility of emergency braking of a first train at the beginning of the cycle of information transfer, measurement methodological and instrumental errors. Digital proportional integral differentiating control law has been considered as an operator, transforming mismatches at the choice of a control of a second locomotive. Automatic control system model has been developed which relates to impulse system class and which feature is the calculation of control action in train traffic process. The parametric synthesis of automatic control system is recommended to be held by simulation modelling methods. Control quality criteria have been chosen.
Keywords: virtual replacer; emergency braking; service braking; automatic control system; radio channel; virtual replacer length; system model; simulation modelling; control law; measurement errors; virtual replacer length calculation; control quality criteria.
References
1. Baranov L. A., Erofeev E. V., Astrakhan V. I., Golovicher Ya. M., Maksimov V. M., Solyus P. G. Sistema avtomaticheskogo i telemekhanicheskogo upravleniya elektro-podvizhnom sostavom [System of automatic and telemechanical control of electric rolling stock]. Moscow: Transport Publ., 1984. 311 p. (In Russian)
2. Baranov L. A. 35 let kafedre «Avtomatika i telemekhanika» [35 years of the department "Automation and telemechanics". Automatic control of technological processes in transport]. Avtomaticheskoe upravlenie tekhnologicheskimi protsessami na transporte/Yubileynyy sbornik nauchnykh trudov [Anniversary collection of scientific papers]. I. 892. Moscow: MIIT Publ., 1996, pp. 3-8. (In Russian)
3. Bushuev S. V., Gundyrev K. V., Golochalov N. S. Povyshenie propusknoy sposob-nosti uchastka zheleznoy dorogi s primeneniem tekhnologii virtual'noy stsepki [Increasing the capacity of a railway section using virtual coupler technology]. Avtomatika na transporte [Transport automation research]. I. 1, 2021, pp. 7-20. DOI: 10.20295/2412-9186-2021-7-1-7-20. (In Russian)
4. Klimova E. E., Pilinushka L. E., Ryabov V. S. Tekhnologiya virtual'noy stsepki poezdov kak instrument, povysheniya propusknoy i provoznoy sposobnosti linii [The technology of virtual train coupling as a tool to increase the throughput and carrying capacity of the line]. Transportnaya infrastruktura sibirskogo regiona: Materialy desyatoy Mezhdunarodnoy nauchno-prakticheskoy konferentsii [Transport infrastructure of the Siberian region: Proceedings of the tenth International scientific and practical conference]. 2019, pp. 60-64. (In Russian)
5. Flammini F., Marrone S., Nardone R., Petrillo A., Santini S., Vittorini V. Towards Railway Virtual Coupling. - International Conference of Electrical Systems for Aircraft, Raiway, Ship Propulsion and Road Vehicles and International Transportation Electrification Conference. 2019.
6. Mitchell I., Goddard E., Montes F., Stanley P., Muttram R., Coenraad W., Pore J., Andrews S., Lochman L. ERTMS Level 4, Train Convoys or Virtual Coupling. Institution of railway signal engineers, IRSE news issue 219. 2016.
7. Goikoetxea J. Roadmap Towards the Wireless Virtual Coupling of Trains. Springer International Publishing Switzerland, 2016, pp. 3-9.
8. Nikitin A. B., Kushpil'I. V. Rezul'taty issledovaniya tekhnicheskikh sredstv kontrolya tselostnosti poezdov [Results of the study of technical means for monitoring the integrity of trains]. Avtomatika na transporte [Transport automation research]. I. 4, 2020, pp. 411-434. DOI 10.20295/2412-9186-2020-6-4-411-434. (In Russian)
9. Quante F. Kontrol'polnosostavnosti gruzovykh poezdov [Control of the completeness of freight trains]. Zheleznye dorogi mira [Railways of the world]. 2005, I. 2, pp. 43-48. (In Russian)
10. Bestem'yanov P. F. Energy efficient algorithms assessment of the Rail circuit operation. Russian Electrical Engineering. 2017, vol. 88, I. 9, pp. 557-562. (In Russian)
11. Baranov L. A. Kvantovanie po urovnyu i vremennaya diskretizatsiya v tsifrovykh sistemakh upravleniya [Level quantization and temporal discretization in digital control systems]. Moscow, Energoatomizdat Publ., 1990. P. 306. (In Russian)
12. Baranov L. A., Golovicher Ya. M., Erofeev E. V., Maksimov V. M. Mikroprotsessornye sistemy avtovedeniya poezdov [Microprocessor-based systems of automatic train guidance]. Moscow: Transport Publ., 1990, p. 272. (In Russian)
13. Rozenfel'd V. E., Isaev I. P., Sidorov N. N. Teoriyaelektricheskoytyagi[Theory of Electric Traction]. Moscow: Transport Publ., 1995, p. 294. (In Russian)
14. Baranov L. A., Erofeev E. V., Mileshin I. S., Chin' L. M. Optimizatsiya upravleniya dvizheniem poezdov [Optimization of train traffic control]. Moscow: MIIT Publ., 2011, p. 164. (In Russian)
15. Lisitsyn A. L. Nestatsionarnye rezhimy tyagi (Tyagovoe obespechenie perevozoch-nogoprotsessa) [Non-stationary traction modes (Traction support for the transportation process)]. Moscow: Intekst Publ., 1996. 159 p. (In Russian)
16. Gao, Gj., Chen, W., Zhang, J. et al. Analysis of longitudinal forces of coupler devices in emergency braking process for heavy haul trains. J. Cent. South Univ. 2017, pp. 2449-2457. DOI: 10.1007/s11771-017-3656-9
17. Pugi L., Fioravanti D., & Rindi A. Modelling the longitudinal dynamics of long freight trains during the braking phase. 12th IFToMM World Congress. Besancon, France. 2007, pp. 1-6.
18. Cole C. Longitudinal train dynamics. In: Iwnicki S (ed) Handbook of railway vehicle dynamics. Taylor & Francis, London. 2006, pp 239-278.
19. Wu Q., Cole C., Luo S., Spiryagin M. A review of dynamics modelling of friction draft gear. Veh Syst Dyn. 2014, I. 52(6), pp. 733-758.