CC BY
90
УДК 681.51
СТРАТИФИЦИРОВАННАЯ МОДЕЛЬ ДЛЯ ОЦЕНКИ И АНАЛИЗА ЭФФЕКТИВНОСТИ ФУНКЦИОНИРОВАНИЯ СЛОЖНЫХ ТЕХНОЛОГИЧЕСКИХ СИСТЕМ СО МНОГИМИ СОСТОЯНИЯМИ
С.Ш. Баласанян
Государственный инженерный университет Армении, г. Капан, Армения E-mail: [email protected]
Предлагается обобщенная стратифицированная модель для оценки и анализа эффективности функционирования сложных технологических систем со многими работоспособными состояниями, учитывающая влияние надежности элементов. Показано, что применение предложенной модели в качестве схемы формализации при построении имитационной модели функционирования сложных технологических систем позволяет существенно упростить и ускорить построение модели системы, обеспечивая удобство проведения имитационных экспериментов. Применение предложенной модели иллюстрировано на примере моделирования технологической системы измельчения руды.
Ключевые слова:
Состояние, стратифицированная, страта, надежность, измельчение. Key words:
State, stratified, strata, reliability, grinding.
Благодаря структурной и функциональной избыточности сложные технологические системы обладают многими работоспособными состояниями. В отличии от бинарных систем, обладающих лишь двумя возможными состояниями (работоспособным и неработоспособным), для таких систем (multi-state systems) [1] практически невозможно определить общепринятое понятие отказа. Поэтому для систем со многими состояниями вместо надежности вводится понятие технической эффективности, оценка которого производится с помощью специально выбранных показателей эффективности, учитывающих последствие влияния отказов элементов системы на качество ее функционирования.
Методологической основой существующих методов оценки и исследования эффективности функционирования сложных технологических систем (СТС) со многими состояниями служит концепция системного подхода. Эта концепция в данном случае проявляется в том, что показатель эффективности рассматривается как функционал от процесса функционирования системы.
В рамках указанного подхода оценка эффективности функционирования сложных систем основывается на использовании модели процесса изменения работоспособности элементов системы. Ее сущность заключается в следующем. Формально каждый элемент E0 (г=1,и) системы в любой момент времени может находиться в одном из возможных состояний y0(t)eY0, каждое из которых характеризуется определенным уровнем работоспособности. Совокупность состояний элементов y°(t)=(y10(t),y20(t),_,ym°(t)) в произвольный момент времени однозначно определяет состояние системы.
С течением времени, под влиянием внешних и внутренних случайных факторов, элементы системы переходят из одного состояния в другое. В результате происходит последовательная смена состояний системы в целом.
Случайный «-мерный процесс Р(0={У(0} рассматривается как формализованный процесс изменения работоспособности элементов системы и описывает ее поведение во времени. Каждой реализации у°(0 процесса соответствует определенная траектория в пространстве состояний системы У0 = ^УЕсли обозначить через Р^^+в) 1=1
вероятность того, что формализованный процесс Р(0 в интервале времени У, t+6] имел реализацию а через еДв) - условный показатель эффективности функционирования системы для этой реализации, то показатель эффективности функционирования системы может быть определен как математическое ожидание условного показателя
Е (1 ,в) = | еу„(1,в)с1Руа(1, 1 + в). (1)
У" (1,1+в)
Несмотря на кажущуюся простоту записи, соотношение малопригодно для расчетов из-за чрезмерной трудности определения е°(^в), Ру^^+в) и может быть использовано лишь для оценки эффективности функционирования систем с небольшим числом состояний. Поэтому для оценки и исследования СТС с большим числом состояний применяется метод имитационного моделирования. Используемая при этом имитационная модель представляет собой композицию трех моделей. Одна из них воспроизводит формализованный процесс изменения работоспособности элементов СТС ¥*^), вторая имитирует процесс функционирования системы под воздействием процесса Р(0, а третья служит для формирования значения условного показателя эффективности еД^в) для реализованных в интервале времени [^+в] траекторий системы У^^+в) еР^+в).
Одной из основных проблем, препятствующих построению имитационной модели функциониро-
вания СТС, является преодоление большой размерности (числа состояний) системы. Так, при построении модели системы, состоящей из 20 бинарных элементов, число рассматриваемых состояний превышает 106.
Применение известных методов и способов упрощения моделей в данном случае либо недопустимо снижает точность моделирования, либо резко ограничивает ее возможности.
Одним из возможных путей решения описанной проблемы является применение концепции иерархического многоуровневого (стратифицированного) описания систем [2]. Ее сущность заключается в том, что процесс функционирования системы Р(/) описывается не одной моделью большой размерности, а иерархией сравнительно простых моделей 8м, ^ие/={1,2,...,£}, каждая из которых описывает процесс функционирования исследуемой системы на определенном уровне абстрагирования (детализации).
В настоящей работе предлагается обобщенная стратифицированная модель, позволяющая для определенного класса СТС успешно преодолеть сложность их формализации.
Придерживаясь терминологии, примененной в [2] предлагаемая модель названа строго-иерархической стратифицированной системой (СИСС),
а уровни детализации, включающие стратифицированное описание - стратами. Каждая страта 8м СИСС представляет собой упорядоченную совокупность несвязанных между собой конечного числа элементов, ге/м={1,2,...,рм}, являющихся временными статическими системами без памяти [3].
Структуру связей между элементами СИСС можно представить в виде ориентированного графа (рис. 1).
На каждом м-м уровне описания (страте) моделируемая система разбивается на подсистемы Ем, ¡&1м м-го уровня.
Каждый элемент 8м описывает поведение подсистемы Е м в зависимости от поведения входящих в ее состав и влияющих на ее функционирование подсистем (м—1)-го уровня. В целом каждая страта 8м описывает поведение моделируемой системы, состоящей из подсистем м-го уровня, в зависимости от поведения подсистем (м— 1)-го уровня. Наивысшая страта 8к описывает поведение системы, рассмотренной как одно целое, в зависимости от поведения подсистем предпоследнего уровня. Таким образом, на каждой м-й страте описывается поведение подсистем м-го уровня с точки зрения их внутреннего механизма, откладывая описание их взаимодействия для (м+1)-й страты.
Переменные, с помощью которых описывается функционирование системы на каждой страте, представляют собой выходные характеристики подсистем соответствующего уровня. Чем выше расположена страта, тем более обобщенными переменными описывается функционирование системы. Чем ниже спускаемся по иерархии, тем более детальным становится описание системы. Компоненты системы, являющиеся ее подсистемами на данной страте, более подробно раскрываются на нижних стратах. Обращаясь к нижним стратам, можно детально объяснить поведение отдельных подсистем при изменении работоспособности элементов системы. С другой стороны, чем выше поднимаемся, тем яснее становится последствие влияния изменения работоспособности элементов системы на качество ее функционирования. На каждой страте выявляются новые свойства системы и связи между ее подсистемами, которые не существуют на других стратах.
Математическое описание СИСС включает формальные описания ее элементов и взаимодействия между ними. Как математический объект каждый элемент 8?, являющийся временной статической системой без памяти [3], характеризуется множествами моментов времени Т? (Т? - подмножество множества действительных чисел), входных сигналов X?, выходных сигналов У? (X? и У? - произвольные множества) и оператором Н?. Входные и выходные сигналы являются функциями времени, их значения в момент времени /е Т? обозначаются х?(/) и у?(/) соответственно. Функционирование каждого элемента 8? при заданных множествах Т?, X? и У? детерминировано и описывается однозначным оператором ух?(0=Н?(х?(0), реализующим отображение Х?^У?.
Взаимодействие между элементами и стратами СИСС, а также между системой и внешней средой рассматривается в рамках механизма обмена сигналами [4]. Для построения математической модели сопряжения элементов и страт в СИСС, помимо сформулированных в [4] предположений, введены еще следующие:
• входной сигнал х?(/) еХ?, поступающий к страте 8? в момент времени /е Т, можно рассматривать как упорядоченную совокупность сигналов x:(t)еX:, 1е1?, одновременно поступающих на входы элементов 8?, 1е1? этой страты: х?(/) = (х 1?(/),х2?(/),..,х/(/)). Аналогично
г(0=(М0,У2?(0,---у;:(0). "
• элементы любой страты 8? СИСС могут принимать только сигналы, выдаваемые элементами страты 8?-1.
В силу принятых в [4] предположений каждому входному контакту к? е К? = ^ К? страты 8?
1е1 ?
(К? - множество входных контактов элемента 8? страты 8?) можно поставить в соответствие единственный выходной контакт с? е С?-1 = ^ С?-1
¡е1
страты 8?-1, соединенный с ним элементарным ка-
налом связи. Следовательно можно ввести однозначный оператор С?=Я?(к?), 1е1?,)е1?-\ пеИ?, 1е-Ь?-\ реализующий отображение К?^С?-1. Область значений оператора Я? совпадает со множеством выходных контактов страты 8?-1. Я? будем называть оператором сопряжения страт 8? и 8?-1вСИСС. Множества К? и С?-1вместе с оператором Я? составляют математическую модель сопряжения страт 8? и 8?-1 в СИСС. Оператор Я? задается бинарной ма-
(р- Л (р? Л
трицей А? размером X Т X X п?
1м ; К1=1 /
которой соответствуют входным контактам страты 8?, а строки - выходным контактам страты 8?-1.
Матрица А? по сравнению с таблицей сопряжения [4] обладает тем преимуществом, что позволяет применять аппарат теории матриц как при ее построении, так и при компьютерной имитации взаимодействия между стратами СИСС [5]. Таким образом, СИСС задается упорядоченной совокупностью конечного числа элементов 8?, ?е1, 1е1?, множеств К?, С?-1, ?е1 и семейством операторов Я?, ?е1.
С применением СИСС построена имитационная модель технологической системы измельчения молибденовой руды (ТСИР). При этом, учитывая характер прикладных задач, удобство описания, а также конструктивные и функциональные особенности ТСИР были выбраны следующие уровни описания (страты): страта измельчительных агрегатов, страта измельчительных комплексов и системная страта.
В результате применения предложенной формализации построен трехуровневый граф 0(Е,У), на основе которого построены модели Е? подсистем различных уровней Е?, 1е1?, ?е/={1,2,3} и модель их сопряжения следующим образом.
Страта измельчительных агрегатов.
На данной страте ТСИР представлена совокупностью р1=15 подсистем первого уровня - измельчительных агрегатов Е/, /еТ1={1,2,_,р1}, полученных в результате декомпозиции системы, состоящей из элементов ЕД /е7°={1,2,...,37}. Каждая подсистема Е 1 первого уровня рассмотрена как объект с двумя возможными состояниями - работоспособным (1) и отказа (0). Случайный процесс У1 (/) изменения состояний каждой подсистемы Е 1 вследствие изменения работоспособности элементов рассматривается как формализованный процесс ее функционирования.
Каждая подсистема Е 1 формализована в виде временной статической системы без памяти 8 1, описываемой множествами Т1, X1, У'иоператором.
Множество У 1 выходных сигналов включает два элемента (0 - состояние отказа, 1 - работоспособное состояние). Оператор Щпредставляет собой логическую функцию работоспособности подсистемы Е 1, описывающую условие работоспособности данной подсистемы в зависимости от работоспособности входящих в ее состав элементов.
Страта измельчительньх комплексов.
На этой страте ТСИР представлена в виде совокупности р2=5 измельчительных комплексов -подсистем второго уровня Е2, ¿еР={1,2,...,р2}.
Каждый измельчительный комплекс, в зависимости от работоспособности входящих в его состав или влияющих на его работу измельчительных агрегатов, может работать в одном из технологически возможных режимов, по определенной схеме измельчения. Каждый режим работы измельчи-тельного комплекса ^характеризуется конкретным набором значений его выходных характеристик у2=(у2и, у\2, угде у2и - производительность комплекса в т/ч; у2а - процентное содержание класса <80 мкм в измельченном продукте (пульпе); у\3 - плотность выходной пульпы в г/л.
Обозначим через х,2(0=(х2;2(0,...,х2;л(0) совокупность выходных характеристик состояний в момент / тех измельчительных агрегатов, от работоспособности которых зависит выбор структуры и режима работы измельчительного комплекса Е2. В каждый момент /е 71, когда происходит изменение работоспособности измельчительных агрегатов, в соответствии со значением вектора х2(/) происходит перестроение структуры комплекса Е2 ипереход в новый режим работы, характеризуемый новым значением вектора у,2(/).
Каждый измельчительный комплекс Е2, формализованный в виде временной статической системы без памяти БД характеризуется четверкой (77ХДУ2,Н). В соответствии с формальным описанием СиСС каждое множество 7Д /еРсовпадает с множеством 71=Т0. Входной сигнал х2(/) каждой Б2 представляет собой упорядоченную совокупность выходных сигналов соответствующих подсистем первого уровня.
В качестве выходного сигнала любого элемента ¿^рассматривается вектор выходных характеристик у2(0=(у2и(0,у2;,2(0,у2;,з(0) соответствующего из-мельчительного комплекса Е 2.
Множество Уэлемента Б2представляет собой совокупность значений вектора выходных характеристик подсистемы Е 2, соответствующих всевозможным режимам ее работы.
Каждый оператор ^реализующий отображение Х2^У2представляет собой бинарную поисковую признаковую структуру - алгоритм, который по заданному входному сигналу х2(0 еХД описывающую работоспособность входящих в состав подсистемы Е2и влияющих на ее функционирование подсистем первого уровня, осуществляет поиск и выдачу соответствующего выходного сигнала у2(0еУ2.
Системная страта.
На данной страте ТСИР рассматривается как одна крупная подсистема Е3 третьего уровня, представляющая собой смеситель (пульподелитель), входом которого является совокупность выходных продуктов (пульпы) всех измельчительных комплексов, а выходом - объединенная выходная пульпа, поступающая на флотацию.
Обозначим через у13(/), у23(0, у33(0 значения выходных характеристик подсистемы Е в момент времени соответственно ее производительность (расход твердого), процентное содержание класса <80 мкм и плотность объединенной пульпы. Соответствующие характеристики входных потоков, представляющие собой выходные характеристики измельчительных комплексов обозначим через
х3и(0=у2и(0, х3а(0=у2а(0, х\з(0=у2у(0, ¿еР:
Очевидно, что общая производительность ТСИР у13(/) и процент содержания класса <80 мкм у23(0 в объединенной пульпе связаны с характеристиками отдельных комплексов следующими соотношениями:
У?С) = £ *3дС),
(2)
уЗС) = £*3д(0 ■ хЩ) £х/д(0. (3)
1=1 / 1=1
Для определения зависимости плотности объединенной пульпы у33(/) от выходных характеристик измельчительных комплексов введем следующие обозначения: Б, V, 01 - соответственно плотность, объем и масса твердой фазы пульпы, выработанной ¿-м измельчительным комплексом за единицу времени; Б, V, 0 - соответственно плотность, объем и масса твердой фазы объединенной пульпы, выработанной системой за единицу времени.
Очевидно, что для каждого ¿-го потока пульпы имеет место следующее соотношение:
оу,=а+у - ан„) йе, (4)
где йг - плотность твердой фазы; йе - плотность воды.
Из соотношения (4) находим:
у -(5)
1 dp Щ - <) и
Подставляя выражение (5) в следующее соотношение, связывающее плотность объединенной пульпы с характеристиками отдельных потоков, получим:
в = ау = £ Бу/ £ V =
1=1 / 1=1
=£в ■ а/(в -ае)/£а/ов -4 •
¡=1 / ¡=1
Принимая йе= 1000 г/л и учитывая непрерывный характер процесса смешивания, окончательно получим:
^ в ■ а а
о=£
£-
■=1 в1 -1000/ ^ в -1°°° = £ ХЛ*) ■ х'з (*) /£_ ХдС)
I „3 1 ппп ¿-I
уУ =
^ х3з(*) -1000/ % х3 з(*) -1000
(6)
Соотношения (2, 3, и 6) соответственно описывают операторы Н13, Н23, Н33, составляющие оператор Н3, который реализует отображение Х^Р.
1=1
J_I
1,1 1,2
1
1,3
1,1 1,2 1,3 1,4 1,5 1,6 1,7 1,8 1,9 1,10 1,11 1,12 1,13 1,14 1,15
Рис. 2. Структура имитационной модели
На основании блок-схемы СИСС (рис. 2) построена математическая модель сопряжения элементов (страт) (X", С"-1, Л"), "=1,3.
Построение имитационной модели ТСИР осуществлено с использованием языка С++.
Разработанная имитационная модель использована для решения следующих прикладных задач:
• оценка степени снижения эффективности функционирования ТСИР из-за ненадежности оборудования;
• оценка влияния ненадежности отдельных элементов ТСИР на эффективность ее функционирования;
• выбор наилучшего варианта резервирования грунтовых насосов;
• определение наилучшей структуры резервных связей между измельчительными агрегатами и рациональной стратегии перестроения структуры ТСИР;
• разработка алгоритма ситуационного управления структурой ТСИР;
• определение оптимального числа ремонтных бригад и рациональной дисциплины обслуживания измельчительного оборудования.
В результате имитационного моделирования ТСИР установлено, что из-за ненадежности из-мельчительного оборудования эффективность функционирования системы снижается на 14,5 %.
Применение СИСС в качестве математической схемы формализации при построении имитационной модели ТСИР позволило: а) свести описание системы к описанию ее частей и их сопряжения, предоставляя возможность аналитического описания функционирования системы на первой и третьей стратах; б) уменьшить объем натурных экспериментов и трудоемкость построения модели; в) решать с помощью модели широкий круг прикладных задач на различных уровнях детализации; г) легко переходить к различным модификациям модели системы в процессе имитационных экспериментов; д) распараллеливать процесс имитации, что при реализации модели на многопроцессорных компьютерах позволяет существенно сократить время моделирования.
СПИСОК ЛИТЕРАТУРЫ
1. Levitin G., Amerib V. Multi-state systems with static performance dependent fault coverage // Reliability Engineering and System Safety. - 2008. - V. 93. - № 11. - P. 1730-1739.
2. Месарович М., Мако Д., Такахара И. Теория иерархических многоуровневых систем. - М.: Мир, 1973. - 344 с.
3. Месарович М., Такахара И. Общая теория систем: математические основы. - М.: 1978. - 311 с.
4. Бусленко Н.П. Моделирование сложных систем. - 2-е изд. -М.: Наука, 1978.- 400 с.
5. Баласанян С.Ш. Об одном методе стратифицированной формализации сложных технологических систем с учетом надежности элементов // Моделирование - 2010: Матер. Междунар. научной конф. - Киев, 2010. - Т. 1. - С. 143-150.
Поступила 07.02.2011 г.
