CC BY
56
В.С. СМОРОДИН
ВЕРИФИКАЦИЯ ИМИТАЦИОННЫХ МОДЕЛЕЙ ТЕХНОЛОГИЧЕСКИХ ПРОЦЕССОВ ПРОИЗВОДСТВА С ПЕРЕМЕННОЙ СТРУКТУРОЙ________________________________________________
Abstract: Technological processes of manufacture, described by means of the probabilistic network graphs, are considered. Particularities of imitation and way of formalization of technological process are given. Strategy of verification of simulation models of technological processes of manufacture with the variable structure is offered.
Key words: technological processes of manufacture, strategy of verification, simulation models with the variable structure.
Анотація: Розглядаються технологічні процеси виробництва, які описуються імовірносними сітковими графіками, особливості імітації і засоби формалізації технологічного процесу. Запропонована методика верифікаціі імітаційних моделей технологічних процесів виробництва із змінною структурою.
Ключові слова: технологічні процеси виробництва, методика верифікаціі, імітаційні моделі змінної структури.
Аннотация: Рассматриваются технологические процессы производства, описываемые с помощью вероятностных сетевых графиков, особенности имитации и способы формализации технологического процесса. Предложена методика верификации имитационных моделей технологических процессов производства с переменной структурой.
Ключевые слова: технологические процессы производства, методика верификации, имитационные модели переменной структуры.
1. Введение
Проблема верификации имитационной модели (ИМ) объекта исследования особенно актуальна при разработке ИМ вероятностных технологических процессов производства (ВТПП) из-за непредсказуемого характера взаимодействий агрегатов AMTXOiJ и ASOBi [1]. В тех случаях, когда
ВТПП реализуется на оборудовании, при работе которого могут происходить отказы и аварии, актуальна задача определения с помощью имитационной модели степени влияния надежностных характеристик оборудования и состава ресурсов предприятия на реализацию ВТПП.
Для решения данной задачи требуется верификация имитационной модели производственной системы с переменной структурой. Вероятностный технологический процесс обладает следующими особенностями: наличием возможности отказов оборудования, приводящих к авариям и конкуренции MTXOij за ресурсы, что в результате ухудшает временные
характеристики реализации ВТПП; наличием графовой структуры связей микротехнологических
операций (MTXOj , где i, j = І, n , n - количество узлов графа), реализующих ВТПП, которая
изменяется вследствие возникновения отказов и аварий оборудования в процессе его функционирования; необходимостью исследования во времени реализации технологического процесса динамики использования ресурсов предприятия; необходимостью оценки возникновения предаварийной ситуации для досрочного перевода оборудования ВТПП на резервирование или профилактику.
Эти особенности определяют актуальность решения следующих задач:
- верификация имитационной модели вероятностного сетевого графика (ВСГР), отображающего параллельно-последовательный процесс выполнения технологических операций, для получения доверия к модели;
- определение значения критического времени Tkp реализации ИМ ВСГР;
- оценка значений таких откликов имитационной модели, как суммарные затраты финансов
(У mtij). Важно при этом анализировать также результат влияния состава и размеров ресурсов
ВТПП при И -м варианте его организации (И = 1, Ы) на значения перечисленных откликов модели. Решение поставленных задач осуществляется в рамках технологических процессов производства, описываемых с помощью вероятностного сетевого графика. В работе приводятся особенности имитации, способы формализации подобных технологических процессов и предложена методика верификации их имитационных моделей.
2. Формализация вероятностного технологического процесса производства
Во время функционирования вероятностного технологического процесса производства каждая
микротехнологическая операция из множества {МТХОц, i,] = 1, п} требует для выполнения предоставления на время её реализации десяти типов ресурсов: ресурса времени реализации (К = 10) ; стоимости Сц (затрат) на выполнение операции (К = 9) ; комплектующих изделий ко^ (К = 8); материалов mtij (К = 7); количества бригад п6у (К = 6) и исполнителей п5ц (К = 5) ; количества общих ресурсов пАЦ (К = 4) и места на общих ресурсах ¥4у предприятия; количества ресурсов индивидуального использования п3у (К =3); количества устройств оборудования общего пользования п2у (К = 2) и места на этом оборудовании У2у; количества устройств оборудования индивидуального пользования % (К =1) .
Запросы ресурсов (т^ , Сц , коц , mtij, У4^, У2у ) являются случайными величинами. Запросы количества ресурсов (ZPij (п^ + п6гу)) - постоянными и индивидуальными величинами для каждой МТХОу . Случайные величины задаются с помощью соответствующих функций распределения (Ещ (г)). Таким образом, структура запросов И -го варианта реализации ВСГР множеством {МТХОу} описывается множеством ОИ. В качестве параметров ИМ ВСГР используются
множество {пОК} общего количества ресурсов К -го типа и множество {¥ОК} начального размера
ресурсов общего пользования.
В вектор откликов УИ для И -го варианта имитационной модели ВСГР входят следующие
компоненты: критическое время его реализации (ТКрь), суммарная стоимость уСуИ выполнения
множества {МТХОц}, суммарные запросы комплектующих изделий у коуИ , суммарный расход
материалов у mtijh. Статистиками имитации в I -й реализации ИМ ВСГР являются множество
общее количество используемых комплектующих изделий (Уко^) и материалов
i, 1
ранних и поздних сроков свершения событий, а также резервов выполнения SOBi ({tpil, tm , R}); множество характеристик времени выполнения MTXOiJ ({t ijl, tnijl, t ijl, tno ijl}) [2]. Таким
образом, на входе имитационной модели имеются множества {Gh}, {nOR}, {VOR}, а на выходе модели определяются значения компонентов вектора Yn.
Критерием качества при анализе результатов моделирования ВСГР является функция
min Yh =j({n0R}, {Gh}) . (1)
Оценка h -го варианта организации ВТПП получается путем усреднения векторов откликов, полученных из всех l -х реализаций структуры ВСГР согласно процедуре Монте-Карло (l = 1, N).
3. Особенности реализации имитационной модели вероятностного сетевого графика
Анализ функционирования ВТПП осуществляется с помощью имитационной модели агрегатного типа, которая реализуется благодаря системе автоматизации имитационного моделирования (САИМ) [2].
САИМ реализует агрегатный способ имитации, когда MTXOy и SOBi, используемые при формализации ВТПП, заменяются соответствующими агрегатами-имитаторами AMTXOy и ASOBi в
ИМ ВСГР. Первой особенностью имитации с помощью САИМ является использование процедуры Монте-Карло [3]. Для l -ой реализации в режиме прямой имитации (модельное время растет от 0 до Ткрl) каждому агрегату AMTXOy выделяются ресурсы предприятия на время выполнения
MTXOy, что имитируется соответствующими операторами ожидания WAIT (тУ]-). Перед
выполнением операции в моменты свершения события SOBi фиксируется статистика расхода
ресурсов (у CiJh, у koiJh , у mtiJh) множеством {MTXOi}-} и статистика свершения событий
i, 1 i, 1 i, 1
завершения выполнения операций {AMTXOi}-} . Второй особенностью имитации является то, что
множества {AMTXOy} и {ASOBi} реализуются только двумя универсальными реентерабельными программами-имитаторами. Поэтому алгоритмы программ-имитаторов агрегатов AMTXOy с помощью множества параметров {Gh} формируют запросы на ресурсы ВТПП
и взаимодействуют посредством сигналов с агрегатами ASOBi. Для каждого агрегата AMTXOy и
ASOBi в базе данных САИМ отводится место для хранения запросов ресурсов предприятия и статистик имитации. Имеет место подобие ВТПП с имитационной моделью ВСГР за счет замены MTXOjj и SOBi на соответствующие агрегаты-имитаторы [4].
4. Методика верификации имитационной модели
Верификация ИМ, построенной с помощью системы автоматизации имитационного моделирования (САИМ) агрегатного типа, реализуется следующей последовательностью этапов.
На первом этапе исследования формируется структура имитационной модели ТПОП путем замены МТХОу на ЛМТХОу, БОБу на ASOBi, подготовки таблиц коммутации агрегатов
согласно ВСГР и последующей актуализации таблиц запросов ресурсов агрегатами в базе данных САИМ.
На втором этапе осуществляется подготовка параметров имитации на основе натурных экспериментов с прототипами ТПОП либо экспертным заданием их значений в тех случаях, когда ТПОП проектируется. Результаты измерений значений запросов ресурсов используются для
построения функций распределения. Количество пК ресурсов К -го типа (К = 1,10) , используемых
каждой МТХОу, задается технологом, исходя из априорной информации (чаще всего нормативных
расходов ресурсов каждой МТХОу) с помощью соответствующей функции распределения
(ЕуК (г)). Функции распределения формируются с помощью гистограмм с последующей их
аппроксимацией стандартными распределениями.
На третьем этапе осуществляется «запитка» имитационной модели наборами параметров средствами САИМ, в результате которой вводятся исходные значения параметров в базу данных для И -го варианта реализации имитационной модели ТПОП. Компоновка элементов имитационной модели при этом осуществляется из универсальных «заготовок» программ-агрегатов, находящихся в библиотеке моделей ЫБСОМР системы автоматизации имитационного моделирования. В итоге формируется программа имитационной модели ВСГР с параметрами пространства начальных значений, которые отображаются в виде трех множеств: { пок }, { Уок }, { С }.
На четвертом этапе осуществляется верификация ИМ ВСГР реальному ТПОП. Проводится испытание очередной версии ИМ ВСГР, которое состоит из двух фаз. На первой фазе испытания проводится верификация имитационной модели, представляющая собой завершение этапа создания программы ИМ ВСГР.
Верификация программы имитационной модели ВСГР обычно осуществляется в серединной точке пространства параметров. Основная трудность верификации состоит в том, что запросы ресурсов МТХОу формируются с помощью функций распределения ЕуК (г) с
использованием жребиев 2-го типа [3]. Верификация реализуется следующей последовательностью шагов этапа 4. На шаге 1 этапа 4 проверяется правильность отображения динамики взаимодействия управляющей программы моделирования с программами агрегатов ЛМТХОу и ЛSOBi согласно технологии «проектирования единичной нити» [5]. С этой целью
формируется пять вариантов одной реализации ИМ ВСГР для случая, когда функции распределения представляют собой нормальное распределение с параметрами N(ау , Бу), где
ау и Бу - соответственно среднее значение и средне-квадратичное отклонение запроса ресурсов
ТПОП (Ту , Су , коу, mtу). В варианте 1 для всех запросов ЛМТХОу по функциям распределения
формируются конкретные значения запросов, равные а у - Бу /2 . В варианте 2 значения запросов
ресурсов от ЛМТХОу равны ау - Бу / 4 . Для вариантов 3, 4 и 5 устанавливаются значения запросов
ресурсов ЛМТХОу , соответственно равные ау , ау + Б у /4 и ау + Б у /2 . Таким образом, в текущей
реализации ИМ ВСГР значения вероятностных запросов ресурсов предприятия для реализации {ЛМТХОу} заменяются детерминированными величинами на всем диапазоне нормального
распределения N (ау, Бу) их значений. На шаге 2 этапа 4 каждый из пяти вариантов ВСГР
представляет собой сетевой график с постоянными значениями запросов ресурсов (у , Сук, квук ,
тґук ), где к - номер варианта. Поэтому на основании методики анализа параметров сетевого
графика рассчитываются отклики Уі = (Ткр1, IС у • I к°і]і, I ті уі) отладочного варианта
имитационной модели. На шаге 3 этапа 4 проводится имитационный эксперимент с помощью программы ИМ ВСГР для тестового примера. Проверяется близость (с точностью до заданной
значениями вектора откликов У{, полученными с помощью программы имитационной модели ВСГР. При выполнении неравенств
считаем, что ИМ ВСГР для I -го варианта (I = 1,5 ) верифицирована. Проверка близости расчетных значений вектора откликов с откликами ИМ ВСГР для всех пяти вариантов позволяет убедиться, что «единичная нить» ИМ ВСГР верифицирована.
На шаге 4 этапа 4 проверяется поведение ИМ ВСГР при «большой нагрузке», которая в данном случае представляет собой имитацию на тестовом варианте с использованием жребиев 3го типа [3], обеспечивающих розыгрыш вероятностных характеристик I-й реализации ВСГР (Ту, ,
Су{, коу{, mtyl ). В итоге имитации тестового варианта программы ИМ ВСГР определяется вектор
откликов У1 = (Ткр1, у Су,, У коу1, у mty■l) при I = 1, N имитационной модели согласно процедуре
и у ' , у ' , у
Монте-Карло. Ожидаемое значение компонентов вектора отклика у; = (Т*р,, С*, ко*, т^) ИМ ВСГР вычисляется по формулам
ошибки имитации є) расчетных значений вектора откликов у; = (Т*кр1, С*, к°*, ті*) со
(2)
и проверяется близость расчетных значений вектора откликов Y* = (T*pi. C*. ko*. mt*) со
значениями вектора Yl = (Tkpl. ICiJl. IkoiJl. Imj), полученными в ходе имитации.
г. ] г. ] U ]
Таким образом, выполнение неравенств (2) для всех откликов имитационной модели позволяет убедиться в том, что тестовый вариант ИМ ВСГР верифицирован в случае розыгрыша запросов ресурсов «большая загрузка» по функциям распределения FiJR (z).
Вторая фаза этапа 4 представляет собой испытание верифицированной версии программы ИМ ВСГР с параметрами, значения которых соответствуют серединной точке пространства параметров, согласно методикам, изложенным в работе [б]. Испытание включает в себя оценку точности имитации, проверку длины переходного периода имитации, анализ устойчивости имитационной модели, оценку чувствительности откликов модели к изменениям параметров и оптимизацию состава параметров откликов модели.
Проверка близости средних значений откликов имитационной модели соответствующим характеристикам реального ВТПП осуществляется с помощью САИМ автоматически. По завершении данного этапа делается вывод, что тестовый вариант ИМ верифицирован и прошел все процедуры испытания модели, поэтому имитационной модели ВСГР можно доверять.
5. Выводы
Предложенная методика верификации имитационной модели технологического процесса производства с переменной структурой в процессе функционирования технологического цикла позволяет решать задачу анализа степени влияния надежностных характеристик оборудования и состава ресурсов предприятия для тех случаев, когда технологический процесс реализуется на оборудовании, при работе которого могут происходить отказы и аварии и тем самым может являться потенциальным источником техногенной опасности. Подобная методика имеет перспективу использования, поскольку обеспечивает возможность автоматизации проведения расчетов при оценке надежности и анализе безопасности функционирования производства.
СПИСОК ЛИТЕРАТУРЫ
1. Смородин В.С. Методика контроля и принятия решений при имитационном моделировании технологических процессов опасного производства // Проблемы управления и информатики. - 2006. - № б. - С. 79-91.
2. Смородин В.С. Агрегатно-процессный стенд имитации для контроля реализаций технологии опасного производства // Проблеми програмування. - 2006. - № 4. - С. 7З-8З.
3. Максимей И.В., Серегина В.С. Задачи и модели исследования операций: Учебное пособие. - Гомель: БелГУТ, 1999. - Ч. 2: Методы нелинейного и стохастического программирования. - 10З с.
4. Максимей И.В. Технология имитации и обработки результатов в системе автоматизации моделирования вероятностных технологических процессов производства / И.В. Максимей, А.Н. Гончаров, В.С. Смородин // Реєстрація, зберігання і обробка даних (Data Recording, Storage & Processing). - 2005. - Т. 7, № З. - С. 71-87.
6. Максимей И.В. Имитационное моделирование на ЭВМ. - М.: Радио и связь, 1988. - 222 с.
Стаття надійшла до редакції 20.09.2007
