Доклады БГУИР
2012 № 2 (64)
УДК 681.3;007.003;007.008;65.0
СИНТЕЗ ДИНАМИЧЕСКОЙ СТРУКТУРЫ ВЕРОЯТНОСТНЫХ ПРОИЗВОДСТВЕННЫХ СИСТЕМ
В С. СМОРОДИН
Гомельский государственный университет имени Ф. Скорины Советская, 104, Гомель, 246028, Беларусь
Поступила в редакцию 7 февраля 2012
Предложен новый подход к унификации методов системных исследований в области проектного моделирования сложных систем с вероятностными параметрами их функционирования на основе конечного множества математических моделей. В качестве инструмента реализации предлагаемого подхода используется динамическое имитационное моделирование, базирующееся на разработке динамических имитационных моделей вероятностных производственных систем.
Ключевые слова: классы сложных систем, вероятностные производственные системы, динамические имитационные модели, синтез оптимальной структуры.
Введение
В настоящее время накоплен большой опыт построения математических моделей объектов исследования из различных отраслей науки, техники, промышленности и экономики. Однако проявляющаяся последнее время тенденция к использованию общесистемных принципов и методов исследований [1] в различных областях знаний наталкивается на определенные трудности. Это связано, как известно, с целым комплексом различных причин: наличием больших объемов разнородной информации в различных областях знаний, использованием различных понятийных аппаратов, профессиональной разобщенностью исследователей и других.
Многочисленные попытки унификации системного подхода [2] при решении конкретных задач науки и практики привели к понятию сложной системы как многопараметрического объекта, представимого конечным множеством математических моделей, каждая из которых отражает конкретную группу свойств исходной системы. При этом появляется возможность выделить классы сложных систем со специфическими свойствами, на основании которых разрабатываются методологические принципы построения математических моделей, которые характеризуются единой математической терминологий и могут быть доступны специалистам различных предметных областей.
Такой подход дает основания рассчитывать в дальнейшем на создание необходимой базы, позволяющей исследователю работать с системами любой степени сложности, вне зависимости от ее физической сущности или ограниченности рамками определенной формализации. В тексте данной статьи к классам сложных систем относятся производственные и экономические системы, вычислительные сети, системы управления, системы энергообеспечения, а также другие системы, призванные обеспечивать безопасность жизнедеятельности и производства.
Системный анализ сложных систем положен в основу исследований производственных систем с вероятностными характеристиками их функционирования (вероятностных производственных систем) как класса сложных систем. При этом под вероятностными характеристиками функционирования могут пониматься характеристики надежности функционирования оборудования, задействованного в ходе реализации технологического цикла производства, харак-
теристики выполнения технологических операций, качественные характеристики используемых в процессе производства материалов и комплектующих изделий и другие.
В качестве инструмента для реализации предлагаемого подхода используется динамическое имитационное моделирование, базирующееся на разработке динамических имитационных моделей.
Краткий анализ состояния вопроса
Как показал анализ современного состояния разработок в области исследования и проектирования технологических объектов с вероятностными характеристиками их функционирования [3], проблема моделирования вероятностных производственных систем состоит в недостаточной результативности методов их исследования при увеличении количества учитываемых параметров, в особенности для тех случаев, когда структура технологического цикла изменяется в процессе функционирования объекта исследования.
Это связано, в первую очередь, с многообразием сложных технологических систем, в ходе реализации которых могут изменяться параметры их функционирования и структура технологического цикла; сложностью практических задач при оценке уровня надежности и безопасности потенциально опасных промышленных объектов; необходимостью учета человеческого фактора при выполнении работ на потенциально опасных объектах.
В соответствии с вышесказанным является актуальной задача обеспечения надежности и безопасности функционирования производственных систем, а также повышение эффективности их использования. Одним из подходов на пути решения данной проблемы является проектное моделирование рациональной структуры вероятностных производственных систем как на стадии проектирования новых, так и при эксплуатации существующих технологических линий.
Поэтому разработка новых методов, программных средств анализа и синтеза оптимальной структуры наиболее сложных вероятностных технологических процессов, использующихся при организации работ на ответственных объектах при наличии оборудования, которое может отказывать и служить источником возникновения аварий техногенного характера, являются весьма актуальными с научной и практической точек зрения.
Методика построения компонентов динамической имитационной модели вероятностной
производственной системы
Для построения компонентов динамической имитационной модели вероятностных производственных систем использовалась агрегатная система автоматизации моделирования вероятностных технологических процессов производства, реализующая агрегатный способ имитации сложных систем, а также способ формализации вероятностных производственных систем на основе вероятностных сетевых графиков (ВСГР) [4]. При этом используется шесть типов агрегатов-имитаторов:
АТОРу - агрегат-имитатор выполнения технологической операции (МТХО^);
ASOBi - агрегат-имитатор свершения i -го события в процессе имитации функционирования вероятностной производственной системы (ВПС);
АОВШ г - агрегат-имитатор функционирования оборудования индивидуального пользования с номером г;
АКАМГ и АОВОРг - агрегат-имитатор совместного использования оборудования общего пользования с помощью выделенного канала;
APROCk - процедура-имитатор с номером к ликвидации аварийной ситуации при выполнении имитационной модели (ИМ) ВПС.
Каждый из этих агрегатов представляет собой реентерабельную программу из конечного множества математических моделей, составляющих математическую модель исходной вероятностной производственной системы. Упомянутые модели агрегатов-имитаторов в базе данных модели имеют свой набор переменных и статистик моделирования в соответствии с версией модели. Количество версий равно числу агрегатов данного типа в ИМ ВПС. Динамика реализации алгоритмов этих агрегатов одинакова для каждого типа агрегатов, но каждая из версий
может находиться в различных состояниях, согласно временной диаграмме реализации агрегатов в модельном времени. Рассмотрим кратко динамику реализации алгоритма одной из версий этих агрегатов.
ЛМТХОу представляет собой четырехполюсный агрегат, имитирующий выполнение технологической операции ЫТХО^ в ИМ ВПС. В режиме прямой имитации приходит сигнал от ASOBi, который по соответствующим функциям распределения формирует значения параметров агрегата (т^ ; с^ ; { Vrij }; { т}; { ко^ }) в I - ой реализации ИМ. Затем каждый АТОЦ определяет индивидуальные запросы на ресурсы и оборудование в виде списков запросов (SPINRSijl; SPOBRijl; SPISPijl). Далее происходит обращение к системе распределения ресурсов, которая выделяет требуемые ресурсы на время имитации ЫТХОу . Когда при имитации АТОр) на выделенном ему оборудовании возникают опасные отказы, то интервал времени выполнения АТОр) будет увеличиваться до тех пор, пока не будут завершены восстановительные
работы на соответствующем агрегате-имитаторе оборудования. При этом проводится анализ возможных последствий случившегося опасного отказа оборудования путем розыгрыша с вероятностью Раг жребия первого типа «произошла авария», который позволяет идентифицировать наличие аварии при работе агрегата на оборудовании с номером г . Если возникла сложная авария, то для ее ликвидации используется «штатная» последовательность процедур ликвидации аварий (агрегатов-имитаторов) { АРЯОСк }. Каждый из них обладает своим оборудованием, составом исполнителей и действует по утвержденному алгоритму ликвидации последствий аварийной ситуации. Длительность тргк1 реализации АРЯОСк также является случайной
величиной и формируется по известной до начала имитации функции распределения Ф4к (трг) .
По окончании реализации этой последовательности агрегатов выполнение операции на отказавшем оборудовании продолжается с прерванного места, формируется признак «была авария» (%ак = 1), а время выполнения операции на оборудовании увеличивается на величину, равную суммарному времени выполнения всех процедур АРЯОСк. По окончании имитации последовательности { АРЯОСк } инициируется возврат на ATOPj , где и формируется сигнал, поступающий на г -й вход агрегата ASOBj. При этом передается признак %ак = 1 («была авария» на оборудовании номера к) , а агрегат ATOPj переходит в режим ожидания инверсного сигнала от ASOBj.
Агрегаты ASOBi являются многополюсными с различным числом входов и выходов. Выходы у ASOBi могут быть одиночными и «кустовыми». Из «кустовых» выходов агрегата формируются сигналы двух типов: действительный Sgd , разыгрываемый по вектору вероятностей { Pijkl } и (к — 1) фиктивных сигналов Sgf . Выходы ASOBi нумеруются, поэтому при адресации сигнала указывается номер события i и номер входа г в агрегат ASOBj . Только действительные сигналы Sgd , поступающие в режиме прямой имитации на вход ATOPj, инициируют его работу по изложенному алгоритму. Фиктивные сигналы Sgf обходят алгоритм выполнения ATOPj . При этом у агрегатов ASOBj используется еще один тип выходных «кустовых» сигналов, называемых резервными выходами ASOBj . С их помощью реализуется «технологическое резервирование» в ИМ ВПС. При появлении отказов оборудования во время выполнения ATOPj, которые посылают сигналы на входы ASOBj, внутри сигнала устанавливается признак аварии (лав). Если лав = 1 (это означает, что ASOBj инициируется в ситуации, когда на предыдущем ATOр, была авария), то в этих случаях инициируются «резервные» вы-
ходы. С этих выходов может формироваться любая комбинация сигналов Sgd и Sgf . Если лав = 0 на всех входах ASOBj, то это означает, что аварий до активизации агрегата не было, поэтому активизируются «штатные» ATOPjS.
Если хотя бы на один из входов ASOBj приходят сигналы от ATOPj, во время выполнения которых была авария, то это означает активизацию резервного кустового выхода. Механизм формирования действительных сигналов у «кустовых» выходов третьего типа основан на использовании булевой матрицы ||упг||, которую составляет технолог-эксперт до начала имитации. Наличие единицы на пересечении п -й строки с г -м столбцом в этой матрице означает необходимость включения резервных ATOPjr в том случае, когда во время выполнения АТОРпУ произошла авария (лави = 1). С помощью этого механизма кустовые выходы третьего типа становятся динамическими регуляторами подключения резервных ATOPjr, если на входы ASOBj поступают сигналы от АТОр^, во время выполнения которых на оборудовании, используемом ими, происходили аварии.
Имитация функций оборудования инициируется запусками АТОр^, который может использовать любое число имитаторов оборудования. За каждым АТОр^ закрепляется два типа
оборудования: индивидуального пользования с номером г, имитируемого агрегатами АОВШг; общего пользования, имитируемого агрегатами АОВОРг. На агрегатах оборудования
общего пользования АТОр^Р использует некоторый объем общего ресурса Уг2у на время его
имитации. По окончании имитации занятый объем ресурса возвращается агрегату оборудования общего пользования АОВОРГ. Поскольку общим оборудованием могут пользоваться несколько ATOPj , то синхронизация их доступа реализуется следующим образом: работа оборудования общего пользования имитируется одним агрегатом АОВОРГ и множеством агрегатов-имитаторов каналов общего пользования АКАМГ, число которых равно числу агрегатов ATOPj , использующих совместно ресурсы АОВОРГ. Агрегат АОВОРГ функционирует непрерывно (последовательно переходит из состояния S1 в состояние S2, а затем снова в состояние S1) в течении I -й реализации ИМ ВПС. Длительность нахождения оборудования в этих состояниях определяется при их смене соответственно функциям распределения: k -й интервал безотказного функционирования тБОИ находится по функции распределения Ф1г (тБО), определяя длительность нахождения общего оборудования в состоянии S1; k -й интервал восстановления работоспособности оборудования тВсш также находится по функции распределения Ф2г (хВО), определяя таким образом длительность нахождения общего оборудования в состоянии S2. Это чередование состояний агрегата продолжается до окончания I -й реализации ИМ. Далее весь процесс имитации АОВОРг продолжается аналогичным образом, но уже в (I + 1)-й реализации имитационной модели.
Поскольку АТОР) и АОВОРГ функционируют независимо, каждый согласно своим алгоритмам, то для синхронизации совместного использования несколькими АТОЦ^Р одного и того же АОВОРГ (когда он находится именно в состоянии S1) введено использование агрегата-канала АКАЛ г, имитирующего использование общего ресурса Угц . Количество таких АКАЛ г определяется числом использований АОВОРГ всеми агрегатами АТОР^.
Таким образом, с помощью комбинаций различных типов кустовых выходов и задания числа разветвлений каждого кустового выхода, эксперту-технологу на стадии проектирования
предоставляется возможность динамического регулирования выполнения множества { АТОЦ } в зависимости от сложившейся операционной обстановки.
Синтез динамической структуры технологического цикла на имитационной модели вероятностной производственной системы
Как видно из методики построения компонентов динамической имитационной модели, исходная вероятностная производственная система представляется в виде конечного набора взаимосвязанных математических моделей. Связь между компонентами математической модели исходной ВПС осуществляется посредством синхронизации взаимодействия агрегатов-имитаторов, входящих в состав моделей-компонентов.
На основании значений множества { }, агрегат события ASOBj осуществляет технологическое резервирование по изложенной выше методике. Кроме технологического резервирования, в ИМ ВПС предусмотрено автоматическое индивидуальное резервирование оборудования, которое включается при достижении фактической наработки оборудования его критического значения. До начала имитации эксперт-технолог устанавливает для каждого устройства порог Q0г этой наработки. Фактическая наработка QФr накапливается на индикаторе 1пёг путем добавления к значению этого индикатора при каждом использовании устройства некоторого приращения ДQфГ. С интервалом Дтй значения этих индикаторов посылаются системе принятия решений БРКЕБН для обеспечения оперативного управления процессом имитации.
Окончательный вариант организации рациональной структуры технологического цикла, который соответствует текущему состоянию оборудования, принимается на основе анализа статистической информации, собранной в процессе функционирования имитационной модели.
Важной статистикой реализации ВСГР, лежащего в основе формализации вероятностной производственной системы, является граф критических путей (GR _ КЯР11) в N реализациях к -го варианта ВПС, который получен следующим образом. Каждый элемент I -й реализации ВСГР критический путь KRPlк представляет собой последовательность чередования статистик выполнения агрегатов (ASOBi, ATOPj , ASOBj). При этом статистиками выполнения ASOBi являются четверки значений (i, tрil, ?п11, Ril = 0), а статистиками ATOPj служат тройки значений (у, т^, с^). Используя исходную структуру ВСГР, эти последовательности накладываются, а в итоге формируется вероятностный граф критических путей GR _KRPк . На выходе ASOBi подсчитывается количество запусков ATOPj (у ^) во всех N
реализациях исходного ВСГР. Затем происходит формирование вероятностей Pу запусков ATOPj с выходов ASOBi графа GR _KRPк . Таким образом, кроме четверки статистик ASOBi
(i, tрi, tпi, Ri), на каждом выходе этого агрегата формируется вектор { Pь } вероятностей запусков ATOPj , где 5 - номера компонентов этого вектора вероятностей.
Другой группой статистик { STкl} являются графики расхода ресурсов в модельном времени при I -й реализации ИМ ВСГР: ресурсов (GR1rl) с номером г , материалов и комплектующих изделий (GR2г1) с номером г, суммарного расхода финансовых средств (GR3l). Фиксируются также временные диаграммы использования: оборудования (DIAGll), исполнителей (DIAG2l), бригад исполнителей (DIAGъl).
Таким образом, в конце l -й реализации имитационной модели ВПС множество статистик { STкl } включает в себя следующие графики и диаграммы: { GRЫ }, { GR2г1 }, { GR3l }, { DIAGll },{ DIAG2l }, { DIAGъl }. Эти статистики запоминаются в виде файлов графиков и диаграмм в базе данных имитационной модели и, по решению эксперта-технолога, могут направляться в систему принятия решений (SPRESH) для выполнения анализа в цикле каждой l -й реализации имитационной модели. По завершении реализаций ИМ графики и диаграммы (в
модельном времени t0) накладываются друг на друга системой SPRESH и в виде набора параметров (06) передаются эксперту для выбора окончательного варианта структуры ВПС.
Обсуждение результатов
Новизна класса формализуемых систем производственных систем состоит в том, что в процессе функционирования технологического цикла может произойти изменение его структуры в результате управляющих воздействий аппаратуры управления или оператора технологического цикла вследствие сбоев или отказов функционирования задействованного оборудования. При этом предложенный подход позволяет решать следующие задачи:
- проектного моделирования надежных и безопасных производственных систем с учетом операций резервирования и ликвидации последствий сбоев и аварий оборудования;
- синтеза оптимальной структуры технологического цикла с учетом надежностных характеристик используемого оборудования на стадии проектирования;
- обеспечения динамической имитации класса производственных систем с вероятностными параметрами их функционирования на стадии проектного моделирования с использованием конечного множества математических моделей.
Все вышеперечисленное в совокупности дает возможность ускорения разработки новых проектов технологических комплексов и снижения себестоимости проектирования их оптимальной структуры при наличии элементов потенциальной опасности.
Заключение
Разработанные подход и программное обеспечение проектного моделирования производственных систем с вероятностными характеристиками позволяют: повысить надежность и безопасность производства существующих технологических систем за счет резервирования цепочек технологических операций при возникновении аварийный ситуаций; снизить вероятность возникновения аварийных ситуаций на производстве, на основе построения рациональной структуры технологического цикла на стадии проектирования; обеспечить резервирование потенциально опасных технологических операций при проведении пуско-наладочных работ и монтаж технологических линий, имеющих определенный срок предварительной эксплуатации.
SYNTHESIS OF DYNAMIC STRUCTURE FOR PROBABILISTIC PRODUCTIVE SYSTEMS
VS. SMORODIN Abstract
The new direction is offered near the unitization of methods of system researches in area of project design of the difficult systems with the probabilistic parameters of their functioning on the basis of finite set of mathematical models. As an instrument of realization offered approach a dynamic imitation design being based on development of dynamic simulation models of the probabilistic productive systems is used.
Список литературы
1. Бусленко Н.П. Лекции по теории больших систем. М., 1973.
2. КлирДж. Системология. М., 1990.
3. Смородин В.С. // Математичш машини i системи (Mathematical Machines and Systems). 2007. №1. С. 105-110.
4. Смородин В.С. // Проблемы управления и информатики. 2006. №5. С. 79-91.