УДК 629.46.004.67:51.001.57
В. В. МЯМЛИН (ДИИТ)
ИСПОЛЬЗОВАНИЕ ТЕОРИИ КУСОЧНО-ЛИНЕЙНЫХ АГРЕГАТОВ ДЛЯ ФОРМАЛИЗАЦИИ РАБОТЫ РЕМОНТНЫХ МОДУЛЕЙ ПОТОЧНОЙ ВАГОНОРЕМОНТНОЙ ЛИНИИ С ГИБКОЙ ТРАНСПОРТНОЙ СИСТЕМОЙ
Розглянуто метод формалiзацil процесу роботи ремонтного модулю потоково! лшп для ремонту вагошв з гнучкою транспортною системою. В якосп математично! схеми процесу функцюнування модулю запропо-нована окрема схема кусково-лшшного агрегату з теорп агрегатiв. Ця схема може бути покладена в основу розробки iмiтацiйних моделей для проектування та дослвдження потокових вагоноремонтних лiнiй ново! генераци.
Рассмотрен метод формализации процесса работы ремонтного модуля поточной линии для ремонта вагонов с гибкой транспортной системой. В качестве математической схемы процесса функционирования модуля предложена частная схема кусочно-линейного агрегата из теории агрегатов. Данная схема может быть положена в основу разработки имитационных моделей для проектирования и исследования поточных вагоноремонтных линий нового поколения.
The method of formalization of the repair module work process of the car-repair flow line with flexible transport system is considered. As the mathematical scheme of process of the module functioning, the specific scheme of a piecewise-linear unit from the theory of units is offered. The given scheme can be put in a basis of development of imitating models for designing and research of the car-repair flow lines of new generation.
При проектировании и строительстве существующих депо, использующих «типовые» поточные линии для ремонта вагонов, недостаточно уделялось внимания вопросам предпро-ектного анализа и моделирования технологических процессов. Поэтому неэффективность принятых проектных решений обнаруживалась только в ходе эксплуатации таких линий. А это уже затраченные огромные средства. Невооружённым глазом видна нерациональность такого подхода к разработке сложных производственных систем (СПС) как по времени, так и по затратам.
Особенности вагоноремонтного производства требуют учёта большого числа случайных факторов, оказывающих существенное влияние на работу поточных линий. Задача состоит в том, чтобы своевременно выявить эти факторы и уже с их учётом создавать такие производственные системы, которые могли бы легко адаптироваться к постоянно изменяющейся ремонтной среде. Существующие в настоящее время поточные линии для ремонта вагонов очень плохо приспособлены к условиям стохастического производства. Они очень чутко реагируют на любые колебания трудоёмкостей ремонтных работ. Вместе с тем, сегодня уже возможны и другие структурные варианты поточных линий, например, гибкие. В случае гибкой поточной линии появляется дополнительная свобода действий, которая проявляется в зна-
чительном увеличении числа альтернативных вариантов принятия управленческих (проектных) решений и согласования взаимодействия между отдельными производственными элементами. Это позволяет сглаживать многие производственные ситуации, вызванные нестабильностью ремонтной среды.
В связи с этим остро встаёт вопрос, связанный с более детальными проработками технологических решений ещё на стадии проектирования объектов. Одним из самых мощных инструментов на этой стадии должно явиться моделирование. По сути дела, моделирование -это имитация практической апробации будущего производства. Ранее автором уже рассматривались некоторые подходы к имитационному моделированию работы поточных вагоноремонтных линий с разными структурами [1-5].
При разработке модели имеют место различные степени абстрагирования. Остановимся только на крайних вариантах. Нижним уровнем выступает физическое моделирование, при котором функционирование объекта проверяется на модели, отличающейся от оригинала, в лучшем случае, размерами. Физические модели с учётом требований, вытекающих из теории подобия, соответствуют неплохой адекватности исследуемому объекту, а, следовательно, и высокой достоверности результатов моделирования. Общими недостатками таких моделей являются большая стоимость, отсутствие универ-
сальности и длительная продолжительность эксперимента.
Верхним уровнем является математическое моделирование. Оно позволяет избежать существенных затрат, значительно сократить время проектирования, исключить метод натурных проб и ошибок. Поэтому для описания процессов, происходящих на поточной линии, мы будем пользоваться исключительно математическими моделями.
Согласно технологическому процессу, каждый ремонтируемый вагон должен последовательно пройти через ряд специальных позиций (участков): мойки, диагностики, разборки, ремонта, сборки, окраски и сушки. Кроме того, вагоны должны перемещаться при помощи трансбордерных тележек и где-то находиться в период ожидания выполнения следующего комплекса работ. Таким образом, гибкая поточная линия для ремонта вагонов представляет собой сложную организационно-технологическую систему, состоящую из отдельных элементов (подсистем), которые находятся в постоянном взаимодействии друг с другом.
При расчленении поточной линии как сложной системы на отдельные элементы будем считать, что для достаточной глубины проработки в качестве отдельных элементов целесообразно принимать всевозможные технологические модули. Такими модулями могут быть «ремонтный модуль» (РМ), «транспортный модуль» (ТМ), «модуль для ожидания» (ОМ).
В данной работе мы остановимся пока только на математической схеме функционирования ремонтного модуля. Под «ремонтным модулем» (РМ) будем понимать часть территории специализированной позиции, необходимой для размещения одного вагона, оснащённой специальным оборудованием, укомплектованной производственным персоналом и предназначенной для выполнения определённого комплекса технологических операций. На каждой позиции может находиться несколько однотипных, взаимозаменяемых РМ. Ремонтные модули одной позиции полностью идентичны. В настоящее время возможны различные структурные схемы организации поточных линий. В зависимости от вариантов компоновки поточной линии РМ могут быть связаны между собой напрямую (используется традиционный грузо-ведущий конвейер) либо через трансбордерные тележки. В варианте жёсткой поточной линии количество РМ соответствует фронту работы (число одновременно ремонтируемых вагонов). Для гибкой поточной линии количество РМ всегда превышает фронт работы поточной линии.
Любые процессы физической природы принято рассматривать в пространстве и во времени. Даже в самом примитивном случае механического движения возникает вопрос о перемещении, которое определяется изменением положения тела в пространстве в течение определённого интервала времени. Рассматривая эту концепцию в более широком смысле, и распространяя её на другие явления, в том числе производственные, под «положением» будем понимать «состояние»: г(^) - состояние системы в момент времени t, а под «перемещением» - «процесс». Процесс представляет собой постоянный «переход» системы из одного состояния в другое.
Рассматриваемая нами система - поточная линия для ремонта вагонов - носит дискретный стохастический и динамический характер.
В дальнейшем, наравне с уже приведенной терминологией, будем также пользоваться известной терминологией теории массового обслуживания.
Поведение системы будем рассматривать на некотором интервале времени [0,7]. Будем считать, что характеристики системы изменяются во времени дискретно. Эти условия отражают динамический характер функционирования системы. Авторы [6] определяют динамическую систему как структурированный объект, в который периодически можно вводить и из которого можно выводить вещество, энергию и информацию. Это свидетельствует о том, что такая система обязательно должна иметь вход и выход.
Описание такой системы в виде единого вероятностного процесса было бы слишком громоздким. Поэтому напрашивается решение, связанное с моделированием каждой отдельной подсистемы и дальнейшей увязкой их между собой. Анализ такой системы, определение её структуры и основных параметров возможно за счёт использования математических имитационных моделей.
При выборе схемы формализации процесса функционирования СПС очень важно учесть два принципиальных момента: 1. Получить наиболее простую модель процесса функционирования системы; 2. Добиться как можно более точных результатов расчёта. В нахождении компромисса между этими двумя противоречивыми условиями и состоит задача оптимального моделирования.
Унифицированная имитационная модель функционирования различных систем может быть представлена в виде стандартной математической схемы - агрегата [7-8]. Такая схема предназначена для изучения сложных структур, состоящих из элементов, представленных ди-
намическими системами в широком понятии: непрерывно-детерминированными, дискретно-детерминированными, дискретно-стохастическими и непрерывно-стохастическими.
В качестве математического описания модулей для ремонта вагонов попытаемся воспользоваться моделью кусочно-линейного агрегата (КЛА) [7-9]. КЛА представляет собой частный случай динамической стохастической системы с дискретным казуальным вмешательством.
КЛА принадлежит к категории объектов, которые можно представлять в виде преобразователя информации (рис. 1). Он функционирует во времени t е Т, может получать входные сигналы х со значениями из некоторого множества X, посылать выходные сигналы у со значениями из некоторого множества У и находиться в определённые моменты времени в некотором состоянии г из некоторого множества 2.
Х(1)
У(0
Рис. 1. Общий вид кусочно-линейного агрегата
Динамике КЛА присущ «событийный» характер. Все события будем делить на две группы: внутренние и внешние. Внутренние события состоят в достижении траекториями КЛА некоторого подмножества состояний г е 2, внешние - в получении входного сигнала.
Как уже было отмечено ранее, КЛА имеет вход и выход. На вход КЛА в дискретные моменты времени ^ поступают входные сигналы. Будем считать, что длительность сигнала очень незначительная - практически мгновенная. В качестве поступающих сигналов может быть информация о том, что перед обслуживающим каналом появилось требование (ремонтируемый вагон), которое нуждается в обслуживании или, например, что требование уже поступило в обслуживающий канал, и он может приступать к обслуживанию.
На выход КЛА поступают выходные сигналы. Выходной сигнал у(?') принадлежит множеству У, у(^ е У. В качестве выходного сигнала может быть, например, информация о том, что канал закончил обслуживание требования или, например, требование покинуло канал (РМ свободен). В нашем случае требование не всегда сразу же покидает канал после обслуживания - должно быть свободное место в очереди к следующему каналу.
КЛА в промежутке между дискретными значениями времени поступления сигналов может находиться в одном из возможных со-
стояний: 1. Агрегат находится в стадии ожидания поступления требования (вагона); 2. Агрегат находится в стадии обслуживания требования; 3. Агрегат уже обслужил требование и находится в стадии ожидания, когда оно покинет агрегат; 4. Требование находится в агрегате, но агрегат его не обслуживает по причине отказа (поломки); 5. Агрегат находится в стадии ремонта, требование в агрегате отсутствует.
Ремонтный модуль (РМ) представляет собой одноканальную систему массового обслуживания (СМО) с неограниченным ожиданием. Правда, в отличие от обычных «классических» СМО, когда требование может без особых проблем стать в очередь на обслуживание или легко покинуть систему после обслуживания, данная СМО имеет ряд специфических особенностей. Например, учитывая большие габариты требований (ремонтируемые вагоны), не всегда может быть в наличии необходимое количество мест для ожидания или, например, не всегда может быть готовым транспортное средство, предназначенное для перемещения требований. Это всё накладывает дополнительные условия на её функционирование. Есть и другие особенности.
Представим СМО в виде кусочно-линейного агрегата. КЛА являет собой частный случай стохастической динамической системы с дискретным казуальным вмешательством.
КЛА функционирует следующим образом. В дискретные моменты времени ^ перед КЛА появляются требования на обслуживание. Каждое требование характеризуется параметром а. Если агрегат свободен, то требование начинает обслуживаться. Если агрегат занят, то требование становится в очередь. В общем случае, когда имеется несколько агрегатов, расположенных параллельно, то требование ожидает, какой из агрегатов освободиться раньше.
Состояние КЛА в момент времени t е Т будем обозначать вектором г(?) е 2, где 2 - множество состояний КЛА. Переход КЛА из состояния г1{() в состояние г2(^) и так далее происходит за очень непродолжительный момент времени, практически равный нулю, т.е. - скачкообразно. Момент скачка определяется входными сигналами х(0 е X и внутренними параметрами самого КЛА, g(t) е G.
Будем считать, что в начальный момент времени ^ е Т система пребывает в начальном состоянии гн е 2, где гн - точка, находящаяся внутри замкнутого «-мерного эвклидова пространства 2. Проследим за изменением состояния системы от действия внутренних причин. Под действием этих причин система может переходить из состояния гн в другие состояния е 2, совершая движения г(0, где t > Точка
^ будет перемещаться в пространстве до тех пор, пока не выйдет на границу области ^ Момент времени выхода на границу обозначим через V.
Для всех V совокупность (V, zt) таких, что ^ < V < V, а ^ = ,г(0 представляет собой фрагмент движения на отрезке (н, V). Этот фрагмент движения назовём перемещением точки ^ внутри пространства ^ Для задания этого перемещения следует указать соотношения, которые определяют ^ для V е (н, V) по заданным значениям tн и гн.
тт *
Для определения момента V и самого состояния системы в этот момент z(t) следует решить совместно уравнения движения точки и уравнения, которые описывают границу пространства состояний Z.
Моменты времени получения входных или
выдачи выходных сигналов назовём «кон* А
трольными» моментами времени V . А состояния КЛА в контрольные моменты времени -«критическими» состояниями z(t). После «критического» состояния КЛА скачком может перейти в новое состояние.
В основу формализации общей схемы процесса функционирования системы положим следующие принципы: 1. Система функционирует во времени и в каждый момент времени может находиться в одном из возможных состояний; 2. На вход системы поступают сигналы; 3. Система может подавать выходные сигналы; 4. Состояние системы в некоторый момент времени определяется предыдущим состоянием системы, а также входными сигналами, которые поступили в данный момент времени и ранее; выходной сигнал в некоторый момент времени определяется состояниями системы и входными сигналами, относящимися к настоящему и предыдущим состояниям. Попытаемся дать формальную интерпретацию каждому принципу.
КЛА в дискретные моменты времени V вых подаёт выходные сигналы у. Выходной сигнал у принадлежит множеству У, у е У и определяется по состоянию агрегата z(t) при помощи оператора выхода Е.
Кроме состояния z(t), будем также рассматривать и состояние z(t+0). Условимся, что для всякого (V1 > V) момент ^+0) находится в полуинтервале (V, Для любого момента времени V состояние агрегата z(t) может быть получено по предыдущим состояниям с помощью случайного оператора перехода Р.
В начальный момент времени 4 система находится в состоянии zн. Будем считать, что процесс функционирования агрегата в момент поступления входного сигнала х описывается подоператором перехода Р1. Поэтому в момент
прихода в агрегат V вх е T входного сигнала х состояние системы можно определить следующим образом: z(t*вх+0) = Р^вх, z(t*вх),х]. *
Если полуинтервал времени (V вх1, V вх2] не
содержит моментов поступления входных сиг* * *
налов, кроме V вх2, то для V е (V вх1, V вх2] состояние КЛА определяет подоператор Р2, z(t) = Р2[^
^'вх1, z(t*вх1+0)].
Оба подоператора Р1 и Р2 рассматриваются как совокупный оператор перехода КЛА в новое состояние.
Рассмотрим работу оператора выхода Е. Представим его в виде двух подоператоров Е1 и Е2. Подоператор Е1 определяет моменты выдачи выходных сигналов, а подоператор Е2 - их содержание.
Момент подачи выходного сигнала определяется следующим порядком. Пусть задано подмножество состояний агрегата Zy, при достижении которых он должен посылать выходной сигнал. Подоператор Е1 представляет характеристическую функцию подмножества 2у: если z(t) е Ху, то Е1 = 1, а если z(t) ё 2у, то Е1 = 0. Поэтому математическая модель формирования выходного сигнала будет иметь вид
y(t') = Е1^ вых, ^ вых), £у] Е2^ вых, ^ вых)].
Учитывая всё выше изложенное, под кусочно-линейным агрегатом (КЛА) будем иметь ввиду некий объект, представляемый совокупностью множеств Т, X, У, Z, Ху, G и случайными операторами Р и Е.
Известно, что каждая ремонтная позиция представляет собой систему массового обслуживания [4]. В общем случае СМО состоит из двух подсистем: подсистемы ожидания обслуживания и подсистемы самого обслуживания. Для нашего случая, когда требованиями выступают крупногабаритные физические объекты -вагоны, должны быть предусмотрены соответствующие площади, на которых могли бы размещаться вагоны в период ожидания обслуживания и площади, на которых осуществляется само обслуживание.
Представим схему КЛА также в виде системы, состоящей из двух подсистем: подсистемы ожидания и подсистемы обслуживания. Подсистема обслуживания является основной подсистемой, а подсистема ожидания - вспомогательной. Каждое требование, чтобы получить обслуживание, обязательно должно пройти через подсистему обслуживания. В подсистему ожидания, требование попадает только в том случае, если в момент его поступления в систему, подсистема обслуживания оказывается занятой.
Ремонтный модуль представим в виде одно-канальной системы массового обслуживания с неограниченным временем ожидания. Требова-
ние будет находиться в очереди до тех пор, пока не попадёт на обслуживание. Покидать систему без обслуживания требование не имеет права. Продолжительность обслуживания требования является случайной величиной т с заданной плотностью распределения / (т).
Состояние КЛА г^) будем образно представлять точкой в многомерном пространстве. Изменение хотя бы одной координаты мгновенно изменяет положение точки и, таким образом, изменяет состояние всей системы.
При описании ремонтных модулей поточной линии в виде кусочно-линейных агрегатов будем считать, что основное состояние соответствует количеству требований (вагонов), находящихся в системе. А вектор дополнительных координат содержит информацию, которая необходима для вычисления дальнейшего протекания процесса г(().
Возьмём в качестве состояния КЛА пару V, 2Ч, т.е. 2 = (V, г„). Параметр V назовём дискретной составляющей состояния или основным состоянием, а zV - вектором вспомогательных координат. Дискретная составляющая
V показывает общее количество требований (вагонов), находящихся в системе (на обслуживании и в ожидании обслуживания), V = (0, 1, 2,..., к). Максимальная величина этого параметра в общем случае определяется количеством обслуживающих каналов п на позиции и вместимостью зоны ожидания перед позицией I; к = (п + I). Параметр V может изменяться с единичной скоростью в сторону уменьшения или увеличения. Вектор г,, может иметь следующие координаты: г, = (, ^), где
^1, ) и С = (С1, С 2,...С I), а
^ - время, оставшееся до окончания обслуживания требования первым каналом;
- время, оставшееся до окончания обслуживания требования вторым каналом;
Ъ,п - время, оставшееся до окончания обслуживания требования п-м каналом;
время, оставшееся до начала обслуживания первого требования в очереди;
^ 2- время, оставшееся до начала обслуживания второго требования в очереди;
- время, оставшееся до начала обслуживания 1-го требования в очереди;
В момент времени, когда Е,к станет равным нулю (к = 1, 2,., п), т.е. обслуживание требования к каналом окончено, происходит изменение состояния системы (скачок). В этот момент дискретный параметр уменьшается на единицу
V = V - 1(обслуженное требование покидает си-
стему). КЛА посылает выходной сигнал о том, что обслуживающий канал свободен. Аналогично, когда станет равным нулю (первое требование из очереди начало обслуживаться), также происходит изменение состояния системы (скачок).
Таким образом, математическая схема кусочно-линейного агрегата может быть положена в основу при разработке имитационных моделей работы гибких поточных линий для ремонта вагонов.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Дударев, А. Е. Применение имитационного моделирования для анализа функционирования поточных вагоноремонтных линий с гибкими связями между производственными участками [Текст] / А. Е. Дударев, В. В. Мямлин. - Д., 1986. - 16 с. - Деп. в ЦНИИ ТЭИ МПС 30.06.86, № 3582.
2. Мямлин, В. В. Использование ЭВМ для анализа функционирования различных поточных линий для ремонта вагонов [Текст] / В. В. Мямлин // Железнодорожный транспорт. Серия «Вагоны и вагонное хозяйство». Ремонт вагонов. - Вып. 1. - М.: ОИ/ЦНИИ ТЭИ МПС, 1989. - С. 1-11.
3. Мямлин, В. В. Применение теории сложных систем к исследованию работы поточных линий для ремонта вагонов [Текст] / В. В. Мямлин // Транспортные связи. Проблемы и перспективы: Докл. Межд. науч.-практ. конф. (Днепропетровск, 29.05 - 30.05.2008). - Д., 2008. - С. 8.
4. Мямлин, В. В. Моделирование работы гибких поточных линий для ремонта вагонов как многофазных многоканальных систем массового обслуживания [Текст] / В. В. Мямлин // Проблемы и перспективы развития железнодорожного транспорта: Тезисы 68 Межд. науч.-практ. конф. (Днепропетровск, 22.05-23.05.2008). - Д., 2008. - С. 51-52.
5. Мямлин, В. В. Совершенствование поточного метода ремонта вагонов за счёт гибкости транспортной системы между технологическими модулями [Текст] / В. В. Мямлин // Залiзн. трансп. Украни. - 2008. - № 4. - С. 15-17.
6. Калман, Р. Очерки по математической теории систем [Текст] / Р. Калман, П. Фалб, М. Ар-биб. - М.: Мир, 1971. - 214 с.
7. Бусленко, Н. П. Моделирование сложных систем [Текст] / Н. П. Бусленко. - М.: Наука, 1978. - 400 с.
8. Бусленко, Н. П. Лекции по теории сложных систем [Текст] / Н. П. Бусленко, В. В. Калашник, И. Н. Коваленко. - М.: Советское радио, 1973. -440 с.
9. Бусленко, В. Н. Автоматизация имитационного моделирования сложных систем [Текст] (Серия «Библиотечка программиста») / В. Н. Бусленко. - М.: Наука, 1977. - 240 с.
Поступила в редакцию 29.08.2008.