CC BY
106
М оделирование систем
УДК 681.513.8
УПРАВЛЕНИЕ БЕЗОПАСНОСТЬЮ И СТОЙКОСТЬЮ СЛОЖНЫХ СИСТЕМ В УСЛОВИЯХ ВНЕШНИХ ВОЗДЕЙСТВИЙ1
А. А. Кочкаров, Г. Г. Малинецкий
Институт прикладной математики им. М. В. Келдыша, г. Москва
Построена математическая модель распространения внешних негативных воздействий по структуре сложной системы. Введены качественные и количественные характеристики, оценивающие подверженность элементов системы негативным воздействиям в зависимости от положения элементов в структуре системы. В ходе исследования модели выявлено несколько синергетических эффектов.
ВВЕДЕНИЕ
Моделирование сложных систем позволяет исследовать особенности их функционирования в различных условиях, наделять их требуемыми характеристиками и снижать риск возникновения чрезвычайных ситуаций (ЧС). Закономерен вопрос — возможно ли в построенной математической модели сложной системы учесть каждый из ее многочисленных элементов?
Рассмотрим проблему с позиций теории самоорганизации — синергетики и теории управления рисками. В математической модели исследуемой системы должны быть представлены основные элементы, по поведению, качеству и эффективности функционирования которых можно достоверно судить о всей системе. Такой подход в исследованиях, когда без детального представления сложных систем, процессов и явлений в них протекающих, принято называть системным синтезом [1]. О результативности этого подхода можно судить по многим работам [1—4]. Его плодотворность наглядно подтверждает цикл работ научной школы В. В. Кульбы [3], посвященный исследованиям по управлению рисками. В этих работах для моделировании поведения систем со сложной структурой были применены методы теории взвешенных ориентированных графов. Такой подход уже позволил обнаружить ряд синергетических эффектов [4] в поведении систем со сложной структурой.
1 Работа выполнена при поддержке РФФИ (проект № 0401-00510).
В настоящей работе этот подход реализован в виде вероятностно-детерминистической модели, описывающей распространение внешних воздействий (различного характера) по элементам исследуемой системы. В основе модели лежат формализации структуры системы в виде ориентированного графа и внешнего негативного влияния на систему в виде импульсного воздействия.
1. НАДЕЖНОСТЬ, ЖИВУЧЕСТЬ И СТОЙКОСТЬ
Под системой в кибернетике принято понимать объединение любых элементов, рассматриваемых как связное целое. Факт непосредственного (без посредников) взаимодействия между двумя элементами системы и определяет наличие связи между ними. Общую картину связей между всеми элементами системы отражает структура системы.
С точки зрения концепции безопасности [5], всякую сложную техническую систему следует изучать с трех основных позиций: надежности системы, живучести системы и ее безопасности.
Надежность [6] — свойство системы сохранять в течение определенного промежутка времени значения параметров, характеризующих функционирование системы.
Теория надежности использует аппарат теории вероятностей и математической статистики. Как правило, для оценки возможности возникновения опасного для окружающей среды состояния системы используется дерево событий (отказов) (рис. 1). Дерево событий (отказов) — это диаграммное представление всех событий (отказов), последовательное и (или) совместное появление которых в системе приводит к некоторому главному
Рис. 1. Двухполюсный граф " и дерево отказов $.
Вершины v1 и v3 — полюсные вершины. На вершину v1 подается сигнал, который должен достичь вершины г3. Главное событие / — недостижение сигнала из вершины v1 вершины v3. Промежуточные события I = {1, 2, 3, 4}, — недостиже-
ние сигналом вершины V. Пунктиром изображены резервное соединение (на графе #) и соответствующее ему событие (на дереве П). Считая надежность (вероятность безотказной работы в течение некоторого промежутка времени) для всех вершин графа равной 0,9, получаем, что надежность функционирования коммуникационной сети в виде графа # без резервного соединения равна \( /) = 1 - (1 - \(/3/4/1))(1 - \(/3/2/1)) х х (1 - \(/3/2/4/1)) = 1 - (1 - 0,93)(1 - 0,93)(1 - 0,94) « 0,975,
с резервным — \( / = 1 - (1 - \(/3/4/1))(1 - \(/3/2/1))(1 -
- \(/3/2/4/1))(1 - ДЗД = 1 - (1 - 0,93)(1 - 0,93)(1 - 0,94)(1 -
- 0,92) I 0,995. Событие / — событие противоположное событию /, т. е. / — функционирование коммуникационной сети #
событию (возможно, потенциально опасному происшествию).
В сложных многоэлементных системах к потенциально опасному происшествию могут привести последовательные и (или) совместные отказы различных элементов системы. Поэтому для повышения надежности элементов (вероятности безотказной работы) системы, и как следствие, надежности самой системы, применяются различные методы резервирования [6] (см. рис. 1).
Живучесть—свойство системы, характеризующее ее способность функционировать под влиянием внешних воздействий (нагрузок), возбуждаемых в окружающей систему среде.
Изучение живучести систем возможно на основе вероятностных моделей, в рамках математической теории надежности [6], и на основе детерминистических, в рамках механики катастроф [5].
Вероятностную модель, описывающую живучесть системы, принято называть прочностной моделью. Под действием внешней нагрузки прочность системы постепенно уменьшается до тех пор, пока система не откажет. Внешние нагрузки описываются случайной величиной (функцией), и как правило, не приводят к скачкообразному изменению прочности системы.
Детерминистическая модель живучести системы лежит в основе механики катастроф [5]. Объектом исследования механики катастроф являются системы, испытывающие постоянные внешние воздействия (нагрузки).
Простым примером таких систем служат инженерные конструкции. В рамках механики катастроф исследуются процессы накопления повреждений, достижения предельного (критического) состояния, реакции элементов конструкций на внешние воздействия.
Особое место в механике катастроф занимает изучение процесса закритического поведения элементов конструкций (систем), который и приводит к тем или иным нежелательным событиям (авариям, катастрофам и т. д.).
Внешние воздействия делят на воздействия природного [1] и техногенного [5] характера. Ко вторым относятся и воздействия, вызываемые умышленными действиями человека. В зависимости от интенсивности и мощности оказываемых на систему воздействий рассматриваются нормативные (проектные) и экстремальные (сверхнормативные) нагрузки. В первом случае изучается живучесть системы в штатных (нормальных) условиях функционирования, когда переход в аварийное состояние возможен при длительном накоплении системой повреждений и достижении предельного (критического) состояния. Во втором случае изучается живучесть системы, когда возможен относительно быстрый переход в аварийное состояние — форс-мажорные обстоятельства.
Механика катастроф занимается не столько изучением различного рода воздействий, сколько созданием аппарата перехода от воздействий к расчетным действующим нагрузкам.
Живучесть и надежность систем являются теми характеристиками, которые позволяют оценить риск возникновения ЧС при эксплуатации сложных технических систем. Используя эти критерии, можно обеспечить безопасность систем при ЧС или наделить системы необходимыми качественными характеристиками, не допускающими возникновения ЧС. В схеме, представленной на рис. 2, надежность и живучесть описывают переход от первого этапа ко второму. Живучесть системы предполагает тщательное описание поведения систем (в отличие от надежности) при имеющихся внешних воздействиях на систему как в докритической области (до ЧС), так и в закритической (при развитии ЧС). Третий этап предполагает изучение возможных последствий ЧС на окружающую систему среду и лежит в области обеспечения безопасности систем.
Согласно рассмотренной схеме безопасность системы можно обеспечить несколькими способами: не допустить развитие ЧС в системе, не допустить выхода ЧС за пределы системы и свести к возможному минимуму влияние аварий на окружающую систему среду.
Напомним, структура системы [7] отражает общую картину причинно-следственных взаимосвязей элементов системы. Во многих случаях сложность системы оп-
Иннциир ое ание ЧС (I этап)
\ Развитие ЧС (П этап) К
Предельное \ состояние / )
V /
Выход ЧС з; пределы системы (Ш этап)
Рис. 2. Схема развития чрезвычайных ситуаций
ределяется сложностью ее структуры. Для подобных технических систем исследование инициирования ЧС (I этап, см. рис. 2) внешними воздействиями имеет особое значение. Важно оценивать, как от структуры системы зависит достижение системой предельного состояния (критического уровня), за чертой которого риск возникновения ЧС резко возрастает.
В рамках модели, предлагаемой в настоящей работе, сложная техническая система считается подвергнутой влиянию внешних воздействий. Это соответствует попаданию системы в зону “форс-мажорных обстоятельств” (под влияние ненормативных нагрузок).
Стойкостью системы назовем ее способность противостоять внешним воздействиям и функционировать в штатном режиме на этапе инициирования ЧС, т. е. в до-критической области функционирования системы. Другими словами, стойкость — это живучесть системы в до-критической области функционирования, под влиянием внешних ненормативных воздействий (нагрузок). Поэтому основной характеристикой стойкости системы будет служить время достижения системой предельного состояния. Увеличение этого промежутка времени будет способствовать уменьшению риска развития ЧС в системе.
2. НАДЕЖНОСТЬ КАК ПОКАЗАТЕЛЬ КАЧЕСТВЕННОГО СОСТОЯНИЯ ЭЛЕМЕНТА СИСТЕМЫ
Существуют различные виды внешних воздействий, оказывающих влияние на систему. Воздействия могут быть механические, термические, электромагнитные, биохимические, радиационные, гидродинамические и др. Как правило, основными считают первые три из перечисленных. Всякое воздействие вызывает повышение нагрузки на те или иные элементы системы. Повышение нагрузки на отдельно взятый элемент системы отражается на состоянии этого элемента выполнять свои функции.
Установление связи между нагрузками а, , а{ , ..., а,- , вызываемыми воздействиями различного рода, как
П
уже отмечалось, является одной из важных задач механики катастроф. Установление такой связи позволяет ввести функционал
] = Ф(а,1 (t), аІ2(t), ..., аПп (t)), (1)
задающий зависимость показателя качественного состояния ] некоторого элемента системы от нагрузок ,
а, , ..., а, , вызванных внешними воздействиями в мо-
12 п
мент времени Р. Для механических систем, все элементы которых являются механическими, функционал (1) имеет вид ]м = Ф(аП(Р), а*(Р)), где показателем качественного состояния ]м элементов системы выступают повреждения, накопленные элементом к моменту времени Р. Деформация аП и напряжение а*, возникающие у элемента системы, служат количественными представлениями полученных при внешних воздействиях нагрузок и поражений.
Как показывает практика, сложные системы состоят из взаимодействующих элементов различной природы (механической, радиоэлектронной и др.). Поэтому не для всех сложных систем удается точно воспроизвести
функционал, определяемый выражением (1), в левой части которого как показатель качественного состояния элемента системы выступают накопленные за время внешнего поражающего воздействия повреждения. В такой ситуации как показатель качественного состояния элемента разумно использовать его надежность. Надежность элемента уменьшается при получении им повреждений, вызванных внешними воздействиями.
Функционирование системы предполагает постоянное взаимодействие образующих его элементов. Поэтому внешние воздействия, причинив повреждения одним элементам системы, окажут влияние на показатели качественного состояния (надежности) и элементов, не подвергнувшихся влиянию этого воздействия непосредственно; т. е. повышение нагрузки будет наблюдаться не только у тех элементов системы, которые оказались под влиянием внешних воздействий, но и у элементов, взаимодействующих (связанных) с ними. Это приведет к понижению надежности последних. Надо отметить, что речь идет только о нагрузках, которые не являются нормативными, т. е. учитываемыми при проектировании и эксплуатации системы. Нормативные нагрузки не должны приводить к повреждениям, резко ухудшающим надежность элементов системы.
3. МАТЕМАТИЧЕСКАЯ МОДЕЛЬ РАСПРОСТРАНЕНИЯ ПОРАЖАЮЩИХ ВОЗДЕЙСТВИЙ ПО СИСТЕМЕ
Следуя известным подходам к повышению надежности систем и их элементов, не представляется возможным продублировать все элементы системы, попадающей под внешнее влияние. Поэтому требуется изучить реакцию системы на “стороннее” влияние, найти наиболее уязвимые “места” системы и рекомендовать их к резервированию.
Опыт исследования многих сложных систем показывает, что на начальном этапе анализа их элементы целесообразно представлять в виде вершин графа, наделенных определенными свойствами, а взаимодействие описывать с помощью ребер. Будем считать тождественными следующие понятия: граф системы и структура системы, вершина графа и элемент системы, ребро графа и связь между элементами системы.
Длявсякого конечного графа будем пользоваться
обозначением — # = (V, Е), где V = {у}, і = 1, п — множество вершин, а Е = {е = (у, и)} — множество его ребер [8].
Распространение воздействия от одного элемента системы к другому будем задавать на графе системы ориентированным ребром — ребром с заданными началом и концом. Ориентированное ребро часто называют дугой, а граф с дугами — орграфом [8]. Орграф структуры моделируемой системы не будет иметь петель (т. е. дуг, конец и начало которых совпадает).
Надежностью элемента системы будем считать вероятность ,(Р < Т) того, что элемент будет работоспособен в течение времени Т с момента начала эксплуатации. Но надежностей элементов, приписываемых соответствующим вершинам графа системы, недостаточно для полного формального представления системы. Воздействие при прохождении от одного элемента к другому теряет свою “силу”. Надежность элемента системы и
долю уменьшения воздействия при переходе от одного элемента к другому можно получить экспериментально или путём экспертного анализа.
Таким образом,
на орграфе О = (2, Е) системы для вершины vi е 2, Е е {1, 2, ..., п} весом = Pv. (Р < Т)
является значение надежности элемента системы, соответствующего вершине V,. А весом w(vi, V) = Ъу, У е {1, 2, ..., п}, , * У, дуги (V,, V) е Е, причем со знаком “+”, является число 0 < Р 1, равное сохранившейся доле передаваемого воздействия, при переходе от вершины V,- к вершине Vj■
Процесс изменения весов вершин графа системы можно отразить следующим правилом, называемым импульсным воздействием. Импульсное воздействие определяется импульсом /тр.(Р), у е {1, 2, ..., п}, в дискретном времени Р = 0, 1, 2, ..., который задается отношением
гтрДР) = wj(f)/wj(t - 1) при Р > 0. (3)
Тогда для Р О 0 для ,-й вершины графа О определим импульсное воздействие
V,-
Wi(t + 1) = Wi.(t) — еу'трДР), (4)
к = 1
или
ёе§ V,-
трр + 1) = — еу'тр/Р), (5)
к = 1
полагая, что dеgvi — число входящих в вершину V,- дуг.
Формулы (3)—(5) задают изменения весов вершин графа О = (2, Е), тем самым определяя динамику распространения внешних воздействий по системе.
Автономное импульсное воздействие на взвешенном орграфе # определим по правилу (3) с вектором начальных значений 3(0) = (w1(0), w2(0), ..., wп(0)) и вектором импульсов
1тр(0) = (,'дар1(0), /тр2(0), ..., трп(0)), (6)
задающим импульс imрj (0) в каждой вершине V- в момент времени Р = 0. Автономное импульсное воздействие в паре с вектором начальных значений описывает состояние системы в начальный момент времени, когда под влияние внешних поражающих воздействий попадают все элементы системы.
Автономное импульсное воздействие, в котором вектор ,тр(0) = (1,1, ..., гтр,.(0), ..., 1), гтр,.(0) > 0, имеет только ,-ю отличную от единицы компоненту, назовем простым воздействием с начальной вершиной V- е 2. Простое импульсное воздействие описывает состояние системы в начальный момент времени, когда внешнее воздействие поражает один из элементов системы. А именно тот, который соответствует ,-й вершине графа системы.
В соответствии с описанным импульсным воздействием на орграфе, можно ввести различные критерии (признаки) достижения системой предельного состоя-
ния. К примеру, можно считать, что система находится в предельном состоянии, если надежность хотя бы одного из наиболее значимых элементов системы, например, элемента V,- е 2 ниже некоторого допустимого уровня. Этот уровень будем называть критическим уровнем надежности элемента и обозначим его как сг(^). Введенный критерий четко разделяет докритическое и за-критическое состояния элемента системы. Если надежность элемента ниже критического уровня, то элемент не в состоянии выполнять возложенных на него функций или функционировать требуемое время.
Представление исследуемой системы в виде взвешенного по правилу (2) графа О = (2, Е) и формализация внешнего влияния на систему как автономного импульсного воздействия (3)—(6) определяет модель распространения поражающих воздействий по системе.
Исследование построенной модели необходимо для решения важной задачи — выяснить, как внешнее воздействие распространяется по структуре системы и влияет на качественное состояние ее элементов.
Поставленная задача требует введения ряда параметров, отражающих подверженность элемента системы внешним воздействиям в зависимости от “положения” элементов в структуре системы.
Структурной уязвимостью у|(и) вершины и е 2назовем число путей, концом которых является вершина и.
Структурная уязвимость элемента дает качественную оценку его расположения в структуре системы. Она позволяет судить о том, насколько безопасно расположение элемента в структуре системы относительно других элементов во время поражающих воздействий.
Важно знать, как к окончанию времени распространения импульсного воздействия по графу системы изменились надежности элементов системы. Предельной надежностью вершины и назовем значение надежности соответствующего ей элемента системы на момент окончания времени воздействия и обозначим ее через Ьг(и).
Можно подсчитать сумму длин всех путей, концом которых является вершина и. Обозначим эту сумму через рэ(и) и назовем её мерой структурной уязвимости вершины и.
Конденсацией [9] графа О = (2, Е) будем называть граф О * = (2*, Е *) полученный из графа О стягиванием в некоторую вершину каждого контура.
Теорема 1. Предельная надежность Ьг(и) вершины и е Е графа О * = (2*, Е *) с равными весами е для всех ребер из Е* при автономном импульсном воздействии с начальным импульсом тр0, одинаковым для всех вершин из 2*, определяется формулой
Ьг(и) = wuimр0<"u'> +1 ерэ(и), (7)
где wu — надежность вершины и в начальный момент автономного импульсного воздействия.
Доказательство. Предельная надежность вершины и е 2* графа О * = (2*, Е *) определяется импульсными воздействиями, отходящими от вершин его множества уязвимости (множества вершин, от которых импульсное воздействие доходит до вершины и). Каждое импульсное воздействие, проходя путь от вершины V е 2* (важно, что таких путей может быть несколько) до вершины и, уменьшается в и) раз, /(V, и) — длина пути
(V, и), и становится равным imр0еl(v’ и). По этому же принципу высчитываются значения всех импульсных воздействий, подходящих в вершине и. А их число, как известно, равно VI (и) — числу путей, концом которых является вершина и. Перемножив, в соответствии с выражением (4), значения всех импульсных воздействий, подходящих к вершине и, получим искомую предельную
и / \ • • VI (и) рэ(и) • VI(и) +1 рэ(и)
надежность Ьг(и) = Wuimрo ,тр0 е = Wuimрo & .
Показатель степени рэ(и) множителя ерэ(и) появляется в произведении (7) как результат сложения длин всех путей, концом которых является вершина и. Множитель тр0 появляется в произведении (7) без уменьшающего его сомножителя е, поскольку одно импульсное воздействие подходит к вершине и непосредственно. ♦
4. ОБ ОСОБЕННОСТЯХ ПОСТРОЕННОЙ МОДЕЛИ
Контуры обратной связи графа # = (V, Е) из модели (2)—(6) потребовали отдельного рассмотрения по следующей причине. Использование конденсации #* модели (2)—(6) позволило заметить, что веса вершины конденсации, соответствующие контурам самого графа # системы, изменяются при отсутствии внешних импульсных воздействий. Это позволило считать контуры графа системы внутренними источниками импульсных воздействий. “Зацикливание” импульсного воздействия в контуре графа системы приводит к повторному изменению весов вершин самого контура и может оказать влияние на веса других вершин графа системы (рис. 3). Внутренние источники импульсных воздействий на графах систем приводят к появлению остаточного эффекта. Для произвольного графа # = (V, Е) эффект выхода из строя элементов системы по происшествию Р = /(/) шагов распространения импульсного воздействия по структуре системы в виде графа # назовем остаточным, где / с # — путь максимальной длины.
“Затухание” внутреннего источника возмущений возможно в двух случаях. Первый — выход из строя контурообразующего элемента остановит циркуляцию импульсного воздействия по контуру. Второй происходит, когда нагрузки, вызываемые импульсными воздействиями, уменьшаться до нормативных, т. е. не будут приводить к резкому снижению надежностей элементов системы. Пример (см. рис. 3) показывает, что если нейтрализация контура не производится искусственно — путем изменения структуры системы, то контур перестанет быть источником импульсных воздействий с момента выхода из строя одного из контурообразующих элементов системы. При этом важно знать, как долго контур будет внутренним источником возмущений.
Началом и окончанием действия контура Scon как внутреннего источника возмущений на графе G будем считать момент первого и, соответственно, последнего изменения надежности вершин (или вершины) образующих контур Scon. Время действия контура Scon в качестве внутреннего источника возмущений будем обозначать как T(S ).
con
Простой контур [9] Scon, вершины которого не являются началами дуг, не входящих в его образование, назовем изолированным.
Теорема 2. Время действия изолированного контура Scon в качестве внутреннего источника возмущений на графе G = (2, Е) с равными весами s для всех ребер из Е и равными весами w0 для всех вершин контура Scon при простом импульсном воздействии imp0, приложенному к одной из вершин контура Scon, определяется формулой
ln Cr(Scon)
S0p0
T(/con) = —"i u u , где cr(Sc0n) — критический уровень
надежности вершин контура, (К — 1) — число полных обходов импульсного воздействия по всем вершинам контура.
U2
Рис. 3. Граф с единственным контуром, образуемым вершинами {и4, и5, иб} (а); изменение (б) весов нескольких вершин графа
при простом импульсном воздействии 1шр(0) = (1; 1; 1; 1; 1; 0,9).
При установленном критическом уровне надежности сг(и) = 0,6, одинаковым для всех / = 1,6 элементов системы со структурой, соответствующей графу Оех, из строя сначала выйдут элементы и2 и и3, только затем контурообразующий элемент и6, при том, что в начальный момент времени надежность первых элементов была выше
74
CONTROL SCIENCES № 5 • 2ПП5
Доказательство. Рассмотрим орграф # = (V, Е) с изолированным контуром Бсоп с #, удовлетворяющий условиям настоящей теоремы. Поскольку простое импульсное воздействие будет приложено к одной из вершин изолированного контура Бсоп, то оно будет распространяться только по вершинам контура. А это позволяет исследовать процесс распространения импульсного воздействии по контуру Бспп, пренебрегая остальными вершинами и дугами графа #
Благодаря структурным особенностям контура Бсоп, простое импульсное воздействие, приложенное к одной из вершин этого контура, будет распространяться только в одном направлении: к вершине, смежной с вершиной, уже достигнутой воздействием. Это, напомним, позволяет легко проследить изменение уровня распространяющегося по контуру импульсного воздействия. В момент времени Р величина импульсного воздействия при передачи от одной вершине к другой, согласно условию теоремы, будет уменьшена в є раз, а по отношению к начальному моменту времени — в ъ раз, и будет равна еітр0. При этом вес каждой вершины контура зависит от того, сколько раз эта вершина испытала на себе влияние импульсного воздействия. Таким образом, можно вывести зависимость
w(P) = wОs ітр0
Щ
(8)
где к(Р) = [Р/Н(/соп)] + 1 — число импульсных воздействий, испытанных вершиной контура к моменту времени Р (квадратными скобками обозначена целая часть числа), а 1(/соп) — длина контура Бсоп. Вообще говоря, слагаемое [Р/Н(/соп)] показателя степени в правой части зависимости (8) есть число полных обходов зависимости [9] импульсного воздействия по всем вершинам контура /соп, а [Р/Н(/соп)@ + 1 — соответственно, номер текущего обхода.
Подставив в зависимость (8) значение критического уровня вершин контура, получим формулу
(9)
для определения времени окончания распространения простого импульсного воздействия по контуру Бсоп Выразим из формулы (9) время Т(/соп) действия контура как внутреннего источника возмущений:
1п СГ(/соп)
0(/ ) = w О ітР О
( соп) К 1 п є
(1О)
Формула (1О) определяет время действия изолированного контура в качестве внутреннего источника возмущений для исследуемой системы. ♦
5. АЛГОРИТМ ПОВЫШЕНИЯ СТОЙКОСТИ СИСТЕМЫ
Рассмотрим систему со структурой в виде ориентированного графа # = (V, Е). Предельное состояние системы определяется предельным состоянием некоторого
множества его элементов V- е 2, , е {1, 2, ■■■, пг}. Пусть система попадает под влияние импульсного воздействия 1тр(0). Для предотвращения перехода системы в предельное состояние необходимо над каждым из элементов V,- проделать следующие операции.
1. Подсчитать предельную надежность Ьг( V ) (см. теорему 1). В случае если среди элементов V,, , е {1, 2, ■■■, пп}, есть элементы контуров обратной связи или элементы, смежные с контурами, подсчет величины Ьг( V,) необходимо вести с учетом времени действия этих контуров.
2. Повысить надежность элемента V- на разницу wi.(0) - Ьг( V,). Повысить надежность элемента V- можно любым из известных методов резервирования [6].
К примеру, для графа системы, изображенного на рис. 1, предельным будем считать состояние, когда надежность элемента Vз опустится ниже отметки 0,4. Тогда предельное состояния этой системой достигается при импульсном воздействии 1тр1(0) = (1; 0,75; 1; 0,75) за время Р1 = 3, а при 1тр2(0) = (0,75; 1; 1; 1) за время Р1 = 3 (учитывая и резервное соединение). И в первом, и во втором случаях время, определяющее стойкость системы, равно трем единицам, но во втором случае система переходит в предельное состояние при меньшем воздействии по сравнению с первым. Для предотвращения перехода системы в предельное состояние под влиянием импульсного воздействия 1тр1(0) достаточно повысить надежность элемента Vз до 0,95.
ЗАКЛЮЧЕНИЕ
Предложенная математическая модель распространения внешних воздействий по системе позволяет объяснить ряд явлений, наблюдаемых в сложных технических системах при попадании их в условия внешних воздействий (форс-мажорные обстоятельства). В частности, становится ясно, как при внешних воздействиях выходят из строя (достигают предельного состояния) элементы, не получавшие импульсного воздействия непосредственно. Показано, что причиной выхода из строя элементов системы вне зоны форс-мажорных обстоятельств (после окончания влияния внешних воздействий) являются внутренние источники возмущений. Существенная особенность построенной модели состоит в демонстрации возможности выхода из строя (перехода в предельное состояние) при распространении импульсных воздействий по системе наиболее надежных элементов. Этот факт красноречиво подчеркивает прямую зависимость надежности элемента от его положения в структуре, а также зависимость стойкости всей системы от выбранной при проектировании структуры. Обнаруженные в ходе исследования модели явления относятся к классу синергетических эффектов [9], прогнозирование которых на начальных этапах исследования нетривиально.
Использование предложенных в работе характеристик (структурная уязвимость, предельная надежность) для оценки состояния элементов при внешних воздействиях, позволяет, с одной стороны, решить оптимизационную задачу — резервирование наиболее подверженных
Рис. 4. Множество уязвимости и граф уязвимости вершины р3 графа в
Граф уязвимости & = (1, 3) вершины у3 на графе О нарисован жирными стрелками. Множество уязвимости 1 = {и1, и2, и3, У2, у3} вершины у3 очерчено пунктирной линией. Внешние воздействия, приложенные ко всем или некоторым элементам множества уязвимости вершины у3 в состоянии привести к выходу из строя элемента системы, соответствующего этой вершине
внешнему влиянию элементов системы. С другой стороны, искусственно направленными воздействиями на множество уязвимости (рис. 4) некоторого элемента системы можно вывести его из строя.
Руководствуясь предложенной моделью, можно говорить о стойкости системе к конкретным воздействия. В случае, когда класс воздействий за время их распространения по структуре системы не переводит ее в предельное состояние, система можно считать абсолютно стойкой к этому классу воздействий.
Модель позволяет предложить рекомендации по наделению системы требуемым уровнем стойкости к внешним поражающим воздействиям с целью уменьшения риска [10] возникновения ЧС в системах, предназначенных для работы в условиях внезапных внешних воздей-
ствий, а при заданных внешних воздействиях — определить значение стойкости с целью предупреждения возможности возникновения ЧС.
ЛИТЕРАТУРА
1. Курдюмов С. П., Малинецкий Г. Г. Синергетика и системный синтез // Новое в синергетике: взгляд в третье тысячелетие. — М.: Наука, 2002.
2. Новое в синергетике: взгляд в третье тысячелетие / Под ред. Г. Г. Малинецкого, С. П. Курдюмова. — М.: Наука, 2002.
3. Методы формирования сценариев развития социально-экономических систем // В. В. Кульба, Д. А. Кононов С. А. Косяченко, А. Н. Шубин. — М.: СИНТЕГ, 2004.
4. Нестационарным структуры и диффузионный хаос // Т. С. Ахромеева, С. П. Курдюмов, Г. Г. Малинецкий, А. А. Самарский. — М.: Наука, 1992.
5. Безопасность России. Правовые, социально-экономические и научно-технические аспекты. Функционирование и развитие сложных народнохозяйственных, технических, транспортных систем, систем связи и коммуникаций / Под ред. К. В. Фролова. — М.: МГФ “Знание”, 1998.
6. Острейковский В. А. Теория надежности. — М.: Высшая школа, 2003.
7. Архипова Н. И., Кульба В. В. Управление в чрезвычайных ситуациях. — М.: РГГУ, 1998.
8. Малинецкий Г. Г. Хаос. Структуры. Вычислительный эксперимент: Введение в синергетику. — М.: Наука, 1997.
9. Лекции по теории графов // В. А. Емеличев, О. Н. Мельников, В. И. Сарванов, Р. И. Тышкевич. — М.: Наука, 1990.
10. Управление риском / В. А. Владимиров, В. В. Кульба, Г. Г. Малинецкий и др. — М.: Наука, 2000.
в (095) 250-79-71; 250-79-16
Е-таіІ: Azret_Kochkarov@mail.ru
GMalin@keldysh.ru □
ттшш Ф'ттмжш <mmk
Статья представляется в редакцию по электронной почте или в одном экземпляре на бумаге, к которому необходимо приложить дискету с текстом, идентичным напечатанному. Аннотация, название статьи и фамилии авторов желательно перевести на английский язык. Объем оригинальной статьи не должен превышать 12, обзорной — 18 страниц текста. Текст печатается через полтора интервала с одной стороны бумаги формата А4, страницы нумеруются. В электронной форме текст должен быть набран в редакторе не ниже Word 97 шрифтом Times New Roman, 12 пунктов.
Рисунки должны иметь расширение, совместимое с Word 97 (Рисунок MS Word (толщина линий — не менее З пкс), редакторы CorelDraw, Photoshop, Illustrator и т.п.).
Фотографии должны быть предельно четкими, на глянцевой бумаге или в формате TIFF с разрешением 3OO dpi (цветовая кодировка CMYK). Все буквенные обозначения на рисунках необходимо пояснить в основном или подрисуночном текстах.
Математические формулы следует набирать в формульном редакторе MathTypes 5.O Equation или MS Equation 3.O, греческие и русские буквы набирать прямым шрифтом, латинские — курсивом. Нумеровать следует только те формулы, на которые есть ссылки в последующем изложении.
Список литературы составляется в порядке цитирования и оформляется в соответствии с действующим ГОСТом. Ссылки на литературу в тексте отмечаются порядковыми номерами в квадратных скобках.
Перед заголовком статьи желательно указать индекс УДК, а в ее конце — номер контактного телефона и электронный адрес.
