Основные понятия анализа опасностей. Отказ, вероятность отказа

Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже

Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.

Подобные документы

Причины техногенных аварий. Аварии на гидротехнических сооружениях, на транспорте. Краткая характеристика крупных аварий и катастроф. Спасательные и неотложные аварийно-восстановительные работы при ликвидации крупных аварий и катастроф.

реферат , добавлен 05.10.2006


Количественная оценка полного риска эксплуатации опасных производственных объектов с помощью математического ожидания ущерба. Формулы расчёта риска аварии, вероятности события, связанного с причинением вреда человеку и окружающей природной среде.

статья , добавлен 01.09.2013


Признаки, позволяющие отнести событие к чрезвычайной ситуации техногенного характера. Причины производственных аварий. Пожары, взрывы, угрозы взрывов. Аварии на коммунальных системах жизнеобеспечения, на очистных сооружениях. Внезапное обрушение зданий.

презентация , добавлен 09.03.2015


Классификация чрезвычайных ситуаций. Краткая характеристика аварий и катастроф, характерных для Республики Беларусь. Аварии на химически опасных, пожаро- и взрывоопасных объектах. Обзор стихийных бедствий. Возможные чрезвычайные ситуации для г. Минска.

реферат , добавлен 05.04.2015


История и виды аварий на гидродинамически опасных объектах, их причины и последствия. Затопление прибрежных территорий в результате разрушения гидротехнических сооружений (плотин и дамб). Меры по уменьшению последствий аварий на опасных объектах.

реферат , добавлен 30.12.2010


Основное понятие об авариях, примерный их перечень. Человеческий фактор как одна из причин аварий. Анализ аварий на шахте "Западная-Капитальная" (Ростовская обл., г. Новошахтинск), шахтах "Ак Булак комур", "Комсомольская", "Юбилейная", "Ульяновская".

реферат , добавлен 06.04.2010


Сущность техногенных аварий. Анализ количества чрезвычайных происшествий и аварий на коммунально-энергетических системах жизнеобеспечения в Республике Хакасия. Динамика аварий на коммунально-энергетических системах в городских муниципальных образованиях.

курсовая работа , добавлен 09.07.2011


Организационные основы осуществления мероприятий по предупреждению и ликвидации последствий аварий и катастроф природного и технического характера. Функционально-организационные структуры поисково-спасательной службы по делам гражданской обороны.

отчет по практике , добавлен 03.02.2013

ОСНОВЫ РАСЧЕТА РАСЧЕТА НАДЕЖНОСТИ ТЕХНИЧЕСКИХ СИСТЕМ ПО НАДЕЖНОСТИ ИХ ЭЛЕМЕНТОВ

Целевое назначение и классификация методов расчета

Расчеты надежности - расчеты, предназначенные для определения количественных показателей надежности. Они проводятся на различных этапах разработки, создания и эксплуатации объектов.

На этапе проектирования расчет надежности производится с целью прогнозирования (предсказания) ожидаемой надежности проектируемой системы. Такое прогнозирование необходимо для обоснования предполагаемого проекта, а также для решения организационно-технических вопросов:
- выбора оптимального варианта структуры;
- способа резервирования;
- глубины и методов контроля;
- количества запасных элементов;
- периодичности профилактики.

На этапе испытаний и эксплуатации расчеты надежности проводятся для оценки количественных показателей надежности. Такие расчеты носят, как правило, характер констатации. Результаты расчетов в этом случае показывают, какой надежностью обладали объекты, прошедшие испытания или используемые в некоторых условиях эксплуатации. На основании этих расчетов разрабатываются меры по повышению надежности, определяются слабые места объекта, даются оценки его надежности и влияния на нее отдельных факторов.

Многочисленные цели расчетов привели к большому их разнообразию. На рис. 4.5.1 изображены основные виды расчетов.

Элементный расчет - определение показателей надежности объекта, обусловленных надежностью его комплектующих частей (элементов). В результате такого расчета оценивается техническое состояние объекта (вероятность того, что объект будет находиться в работоспособном состоянии, средняя наработка на отказ и т.п.).

Рис. 4.5.1. Классификация расчетов надежности

Расчет функциональной надежности - определение показателей надежности выполнения заданных функций (например, вероятность того, что система очистки газа будет работать заданное время, в заданных режимах эксплуатации с сохранением всех необходимых параметров по показателям очистки). Поскольку такие показатели зависят от ряда действующих факторов, то, как правило, расчет функциональной надежности более сложен, чем элементный расчет.

Выбирая на рис 4.5.1 варианты перемещений по пути, указанному стрелками, каждый раз получаем новый вид (случай) расчета.

Самый простой расчет - расчет, характеристики которого представлены на рис. 4.5.1 слева: элементный расчет аппаратурной надежности простых изделий, нерезервированных, без учета восстановлений работоспособности при условии, что время работы до отказа подчинено экспоненциальному распределению.

Самый сложный расчет - расчет, характеристики которого представлены на рис. 4.5.1 справа: функциональной надежности сложных резервированных систем с учетом восстановления их работоспособности и различных законов распределения времени работы и времени восстановления.
Выбор того или иного вида расчета надежности определяется заданием на расчет надежности. На основании задания и последующего изучения работы устройства (по его техническому описанию) составляется алгоритм расчета надежности, т.е. последовательность этапов расчета и расчетные формулы.

Последовательность расчета систем

Последовательность расчета системы представлена на рис. 4.5.2. Рассмотрим основные ее этапы.

Рис. 4.5.2. Алгоритм расчета надежности

Прежде всего четко следует сформулировать задание на расчет надежности. В нем должны быть указаны: 1) назначение системы ее состав и основные сведения о функционировании; 2) показатели надежности и признаки отказов, целевое назначение расчетов; 3) условия, в которых работает (или будет работать) система; 4) требования к точности и достоверности расчетов, к полноте учета действующих факторов.
На основании изучения задания делается вывод о характере предстоящих расчетов. В случае расчета функциональной надежности осуществляется переход к этапам 4-5-7, в случае расчета элементов (аппаратурной надежности) - к этапам 3-6-7.

Под структурной схемой надежности понимается наглядное представление (графическое или в виде логических выражений) условий, при которых работает или не работает исследуемый объект (система, устройство, технический комплекс и т.д.). Типовые структурные схемы представлены на рис. 4.5.3.

Рис. 4.5.3. Типовые структуры расчета надежности

Простейшей формой структурной схемы надежности является параллельно-последовательная структура. На ней параллельно соединяются элементы, совместный отказ которых приводит к отказу
В последовательную цепочку соединяются такие элементы, отказ любого из которых приводит к отказу объекта.

На рис. 4.5.3,а представлен вариант параллельно-последовательной структуры. По этой структуре можно сделать следующее заключение. Объект состоит из пяти частей. Отказ объекта наступает тогда, когда откажет или элемент 5, или узел, состоящий из элементов 1-4. Узел может отказать тогда, когда одновременно откажет цепочка, состоящая из элементов 3,4 и узел, состоящий из элементов 1,2. Цепь 3-4 отказывает, если откажет хотя бы один из составляющих ее элементов, а узел 1,2 - если откажут оба элемента, т.е. элементы 1,2. Расчет надежности при наличии таких структур отличается наибольшей простотой и наглядностью. Однако не всегда удается условие работоспособности представить в виде простой параллельно-последовательной структуры. В таких случаях используют или логические функции, или графы и ветвящиеся структуры, по которым оставляются системы уравнений работоспособности.

На основе структурной схемы надежности составляется набор расчетных формул. Для типовых случаев расчета используются формулы, приведенные в справочниках по расчетам надежности, стандартах и методических указаниях. Прежде чем применять эти формулы, необходимо предварительно внимательно изучить их существо и области использования.

Расчет надежности, основанный на использовании параллельно-последовательных структур

Пусть некоторая техническая система D составлена из n элементов (узлов). Допустим, надежности элементов нам известны. Возникает вопрос об определении надежности системы. Она зависит от того, каким образом элементы объединены в систему, какова функция каждого из них и в какой мере исправная работа каждого элемента необходима для работы системы в целом.

Параллельно-последовательная структура надежности сложного изделия дает представление о связи между надежностью изделия и надежностью его элементов. Расчет надежности ведется последовательно - начиная от расчета элементарных узлов структуры к ее все более сложным узлам. Например, в структуре рис. 5.3,а узел, состоящий из элементов 1-2 - элементарный узел, состоящий из элементов 1-2-3-4, сложный. Эта структура может быть сведена к эквивалентной, состоящей из элементов 1-2-3-4 и элемента 5, соединенных последовательно. Расчет надежности в данном случае сводится к расчету отдельных участков схемы, состоящих из параллельно и последовательно соединенных элементов.

Система с последовательным соединением элементов

Самым простым случаем в расчетном смысле является последовательное соединение элементов системы. В такой системе отказ любого элемента равносилен отказу системы в целом. По аналогии с цепочкой последовательно соединенных проводников, обрыв каждого из которых равносилен размыканию всей цепи, мы и называем такое соединение "последовательным" (рис. 4.5.4). Следует пояснить, что "последовательным" такое соединение элементов является только в смысле надежности, физически они могут быть соединены как угодно.

Рис. 4.5.4. Блок-схема системы с последовательным соединением элементов

С позиции надежности, такое соединение означает, что отказ устройства, состоящего из этих элементов, происходит при отказе элемента 1 или элемента 2, или элемента 3, или элемента n. Условие работоспособности можно сформулировать следующим образом: устройство работоспособно, если работоспособен элемент 1 и элемент 2, и элемент 3, и элемент n.

Выразим надежность данной системы через надежности ее элементов. Пусть имеется некоторый промежуток времени (0,t ), в течение которого требуется обеспечить безотказную работу системы. Тогда, если надежность системы характеризуется законом надежности Р(t), нам важно знать значение этой надежности при t=t , т.е. Р(t ). Это не функция, а определенное число; отбросим аргумент t и обозначим надежность системы просто Р. Аналогично обозначим надежности отдельных элементов P 1 , P 2 , P 3 , ..., P n .

Для безотказной работы простой системы в течение времени t нужно, чтобы безотказно работал каждый из ее элементов. Обозначим S - событие, состоящее в безотказной работе системы за время t ; s 1 , s 2 , s 3 , ..., s n - события, состоящие в безотказной работе соответствующих элементов. Событие S есть произведение (совмещение) событий s 1 , s 2 , s 3 , ..., s n:
S = s 1 × s 2 × s 3 × ... × s n .

Предположим, что элементы s 1 , s 2 , s 3 , ..., s n отказывают независимо друг от друга (или, как говорят применительно к надежности, "независимы по отказам", а совсем кратко "независимы"). Тогда по правилу умножения вероятностей для независимых событий Р(S)=P(s 1)× P(s 2)× P(s 3)× ...× P(s n) или в других обозначениях,
Р = Р 1 × Р 2 × Р 3 × ... × Р n .,(4.5.1)
а корочеP = ,(4.5.2)
т.е. надежность (вероятность работоспособного состояния) простой системы, составленной из независимых по отказам, последовательно соединенных элементов, равна произведению надежностей ее элементов.

В частном случае, когда все элементы обладают одинаковой надежностью P 1 =P 2 =P 3 = ... =P n , выражение (4.5.2) принимает вид
Р = P n .(4.5.3)

Пример 4.5.1. Система состоит из 10 независимых элементов, надежность каждого из которых равна Р=0,95. Определить надежность системы.

По формуле (4.5.3) Р = 0,95 10 » 0,6.

Из примера видно, как резко падает надежность системы при увеличении в ней числа элементов. Если число элементов n велико, то для обеспечения хотя бы приемлемой надежности Р системы каждый элемент должен обладать очень высокой надежностью.

Поставим вопрос: какой надежностью Р должен обладать отдельный элемент для того, чтобы система, составленная из n таких элементов, обладала заданной надежностью Р?

Из формулы (4.5.3) получим:
Р = .

Пример 4.5.2. Простая система состоит из 1000 одинаково надежных, независимых элементов. Какой надежностью должен обладать каждый из них для того, чтобы надежность системы была не меньше 0,9?
По формуле (4.5.4) Р = ; lgР = lg0,9 1/1000 ; Р » 0,9999.

Интенсивность отказов системы при экспоненциальном законе распределения времени до отказа легко определить из выражения
l с = l 1 + l 2 + l 3 + ... + l n ,(4.5.4)
т.е. как сумму интенсивностей отказов независимых элементов. Это и естественно, так как для системы, в которой элементы соединены последовательно, отказ элемента равносилен отказу системы, значит все потоки отказов отдельных элементов складываются в один поток отказов системы с интенсивностью, равной сумме интенсивностей отдельных потоков.

Формула (4.5.4) получается из выражения
Р = P 1 P 2 P 3 ... P n = ехр{-(l 1 + l 2 + l 3 + ... + l n )}.(4.5.5)
Среднее время работы до отказа
Т 0 = 1/ l с .(4.5.6)

Пример 4.5.3. Простая система S состоит из трех независимых элементов, плотности распределения времени безотказной работы которых заданы формулами:

при 0 < t < 1 (рис. 4.5.5).

Рис. 4.5.5. Плотности распределения времени безотказной работы

Найти интенсивность отказов системы.
Решение. Определяем ненадежность каждого элемента:
при 0 < t < 1.

Отсюда надежности элементов:
при 0 < t < 1.

Интенсивности отказов элементов (условная плотность вероятности отказов) - отношение f(t) к р(t):
при 0 < t < 1.
Складывая, имеем: l с = l 1 (t) + l 2 (t) + l 3 (t).

Пример 4.5.4. Предположим, что для работы системы с последовательным соединением элементов при полной нагрузке необходимы два разнотипных насоса, причем насосы имеют постоянные интенсивности отказов, равные соответственно l 1 =0,0001ч -1 и l 2 =0,0002ч -1 . Требуется вычислить среднее время безотказной работы данной системы и вероятность ее безотказной работы в течение 100ч. Предполагается, что оба насоса начинают работать в момент времени t =0.

С помощью формулы (4.5.5) находим вероятность безотказной работы P s заданной системы в течение 100ч:
P s (t)= .
P s (100)=е -(0,0001+0,0002) × 100 =0,97045.

Используя формулу (4.5.6), получаем

ч.

На рис. 4.5.6 представлено параллельное соединение элементов 1, 2, 3. Это означает, что устройство, состоящее из этих элементов, переходит в состояние отказа после отказа всех элементов при условии, что все элементы системы находятся под нагрузкой, а отказы элементов статистически независимы.

Рис. 4. 5.6. Блок-схема системы с параллельным соединением элементов

Условие работоспособности устройства можно сформулировать следующим образом: устройство работоспособно, если работоспособен элемент 1 или элемент 2, или элемент 3, или элементы 1 и 2, 1; и 3, 2; и 3, 1; и 2; и 3.

Вероятность безотказного состояния устройства, состоящего из n параллельно соединенных элементов определяется по теореме сложения вероятностей совместных случайных событий как
Р=(р 1 +р 2 +...р n)-(р 1 р 2 +р 1 р 3 +...)-(р 1 р 2 р 3 +р 1 р 2 р n +...)-... ± (р 1 р 2 р 3 ...р n).(4.5.7)
Для приведенной блок-схемы (рис. 4.5.6), состоящей из трех элементов, выражение (4.5.7) можно записать:
Р=р 1 +р 2 +р 3 -(р 1 р 2 +р 1 р 3 +р 2 р 3)+р 1 р 2 р 3 .

Применительно к проблемам надежности, по правилу умножения вероятностей независимых (в совокупности) событий, надежность устройства из n элементов вычисляется по формуле
Р = 1- ,(4.5.8)
т.е. при параллельном соединении независимых (в смысле надежности) элементов их ненадежности (1-p i =q i) перемножаются.

В частном случае, когда надежности всех элементов одинаковы, формула (4.5.8) принимает вид
Р = 1 - (1-р) n .(4.5.9)

Пример 4.5.5. Предохранительное устройство, обеспечивающее безопасность работы системы под давлением, состоит из трех дублирующих друг друга клапанов. Надежность каждого из них р=0,9. Клапаны независимы в смысле надежности. Найти надежность устройства.

Решение. По формуле (4.5.9)Р=1-(1-0,9) 3 =0,999.

Интенсивность отказов устройства состоящего из n параллельно соединенных элементов, обладающих постоянной интенсивностью отказов l 0 , определяется как

.(4.5.10)

Из (4.5.10) видно, что интенсивность отказов устройства при n>1 зависит от t: при t=0 она равна нулю, при увеличении t, монотонно возрастает до l 0 .

Если интенсивности отказов элементов постоянны и подчинены показательному закону распределения, то выражение (4.5.8) можно записать

Р(t) = .(4.5.11)

Среднее время безотказной работы системы Т 0 находим, интегрируя уравнение (4.5.11) в интервале :

Т 0 =
=(1/ l 1 +1/ l 2 +…+1/ l n )-(1/(l 1 + l 2 )+ 1/(l 1 + l 3 )+…)+(4.5.12)
+(1/(l 1 + l 2 + l 3 )+1/(l 1 + l 2 + l 4 )+…)+(-1) n+1 ´ .

В случае, когда интенсивности отказов всех элементов одинаковы, выражение (4.5.12) принимает вид

Т 0 = .(4.5.13)

Среднее время работы до отказа также можно получить, интегрируя уравнение (4.5.7) в интервале

Пример 4.5.6. Предположим, что два одинаковых вентилятора в системе очистки отходящих газов работают параллельно, причем если один из них выходит из строя, то другой способен работать при полной системной нагрузке без изменения своих надежностных характеристик.

Требуется найти безотказность системы в течение 400ч (продолжительность выполнения задания) при условии, что интенсивности отказов двигателей вентиляторов постоянны и равны l =0,0005ч -1 , отказы двигателей статистически независимы и оба вентилятора начинают работать в момент времени t=0.

Решение. В случае идентичных элементов формула (4.5.11) принимает вид
Р(t) = 2еxp(- l t) - еxp(-2 l t).
Поскольку l = 0,0005 ч -1 и t = 400 ч, то
Р (400) = 2еxp(-0,0005 ´ 400) - еxp(-2 ´ 0,0005 ´ 400)=0,9671.
Среднюю наработку на отказ находим, используя (4.5.13):
Т 0 = 1/l (1/1 + 1/2) = 1/l ´ 3/2 = 1,5/0,0005 = 3000 ч.

Рассмотрим самый простой пример резервированной системы - параллельное соединение резервного оборудования системы. В этой схеме все n одинаковых образцов оборудования работают одновременно, и каждый образец оборудования имеет одинаковую интенсивность отказов. Такая картина наблюдается, например, если все образцы оборудования держатся под рабочим напряжением (так называемый "горячий резерв"), а для исправной работы системы должен быть исправен хотя бы один из n образцов оборудования.

В этом варианте резервирования применимо правило определения надежности параллельно соединенных независимых элементов. В нашем случае, когда надежности всех элементов одинаковы, надежность блока определяется по формуле (4.5.9)

Р = 1 - (1-р) n .
Если система состоит из n образцов резервного оборудования с различными интенсивностями отказов, то
P(t) = 1-(1-p 1) (1-p 2)... (1-p n).(4.5.21)

Выражение (4.5.21) представляется как биноминальное распределение. Поэтому ясно, что когда для работы системы требуется по меньшей мере k исправных из n образцов оборудования, то
P(t) = p i (1-p) n-i ,где .(4.5.22)

При постоянной интенсивности отказов l элементов это выражение принимает вид

P(t) = ,(4.5.22.1)

где р = еxp(-l t).

Включение резервного оборудования системы замещением

В данной схеме включения n одинаковых образцов оборудования только один находится все время в работе (рис. 4.5.11). Когда работающий образец выходит из строя, его непременно отключают, и в работу вступает один из (n -1) резервных (запасных) элементов. Этот процесс продолжается до тех пор, пока все (n -1) резервных образцов не будут исчерпаны.

Рис. 4.5.11. Блок-схема системы включения резервного оборудования системы замещением
Примем для этой системы следующие допущения:
1. Отказ системы происходит, если откажут все n элементов.
2. Вероятность отказа каждого образца оборудования не зависит от состояния остальных (n -1) образцов (отказы статистически независимы).
3. Отказывать может только оборудование, находящееся в работе, и условная вероятность отказа в интервале t, t+dt равна l dt; запасное оборудование не может выходить из строя до того, как оно будет включено в работу.
4. Переключающие устройства считаются абсолютно надежными.
5. Все элементы идентичны. Резервные элементы имеют характеристики как новые.

Система способна выполнять требуемые от нее функции, если исправен по крайней мере один из n образцов оборудования. Таким образом, в этом случае надежность равна просто сумме вероятностей состояний системы, исключая состояние отказа, т.е.
Р(t) = еxp(- l t) .(4.5.23)

В качестве примера рассмотрим систему, состоящую из двух резервных образцов оборудования, включаемых замещением. Для того чтобы эта система работала, в момент времени t, нужно, чтобы к моменту t были исправны либо оба образца, либо один из двух. Поэтому
Р(t) = еxp(- l t) =(exp(- l t))(1+ l t).(4.5.24)

На рис. 4.5.12 показан график функции Р(t) и для сравнения приведен аналогичный график для нерезервированной системы.

Рис. 4.5. 12. Функции надежности для дублированной системы свключением резерва замещением (1) и нерезервированнойсистемы (2)

Пример 4.5.11. Система состоит из двух идентичных устройств, одно из которых функционирует, а другое находится в режиме ненагруженного резерва. Интенсивности отказов обоих устройств постоянны. Кроме того, предполагается, что в начале работы резервное устройство имеет такие же характеристики, как и новое. Требуется вычислить вероятность безотказной работы системы в течение 100 ч при условии, что интенсивности отказов устройств l =0,001 ч -1 .

Решение. С помощью формулы (4.5.23) получаем Р(t) = (exp(- l t))(1+ l t).

При заданных значениях t и l вероятность безотказной работы системы составляет

Р(t) = е -0,1 (1+0,1) = 0,9953.

Во многих случаях нельзя предполагать, что запасное оборудование не выходит из строя, пока его не включат в работу. Пусть l 1 - интенсивность отказов работающих образцов, а l 2 - резервных или запасных (l 2 > 0). В случае дублированной системы функция надежности имеет вид:
Р(t) = ехр(-(l 1 + l 2 )t) + ехр(- l 1 t) - ехр(-(l 1 + l 2 )t).

Данный результат для k=2 можно распространить на случай k=n. Действительно

Р(t) = ехр(- l 1 (1+ a (n-1))t) (4.5.25)
, где a = l 2 / l 1 > 0.

Надежность резервированной системы в случае комбинаций отказов и внешних воздействий

В некоторых случаях отказ системы возникает вследствие определенных комбинаций отказов образцов входящих в систему оборудования и (или) из-за внешних воздействий на эту систему. Рассмотрим, например, метеоспутник с двумя передатчиками информации, один из которых является резервным или запасным. Отказ системы (потеря связи со спутником) возникает при выходе из строя двух передатчиков или в тех случаях, когда солнечная активность создает непрерывные помехи радиосвязи. Если интенсивность отказов работающего передатчика равна l , а j - ожидаемая интенсивность появления радиопомех, то функция надежности системы
Р(t) = еxp(-(l + j )t) + l t еxp(-(l + j )t).(4.5.26)

Данный тип модели также применим в случаях, когда резерв по схеме замещения отсутствует. Например, предположим, что нефтепровод подвергается гидравлическим ударам, причем воздействие незначительными гидроударами происходит с интенсивностью l , а значительными - с интенсивнностью j . Для разрыва сварных швов (из-за накопления повреждений) трубопроводу следует получить n малых гидроударов или один значительный.

Здесь состояние процесса разрушения представляется числом ударов (или повреждений), причем один мощный гидроудар равносилен n малых. Надежность или вероятность того, что трубопровод не будет разрушен действием микроударов к моменту времени t равна:

Р(t) = еxp(-(l + j )t) .(4.5.27)

Анализ надежности систем при множественных отказах

Рассмотрим метод анализа надежности нагруженных элементов в случае статистически независимых и зависимых (множественных) отказов. Следует заметить, что этот метод может быть применен и в случае других моделей и распределений вероятностей. При разработке этого метода предполагается, что для каждого элемента системы существует некоторая вероятность появления множественных отказов.

Как известно, множественные отказы действительно существуют, и для их учета в соответствующие формулы вводится параметр a . Этот параметр может быть определен на основе опыта эксплуатации резервированных систем или оборудования и представляет собой долю отка ов, вызываемых общей причиной . Другими словами, параметр а можно рассматривать как точечную оценку вероятности того, что отказ некоторого элемента относится к числу множественных отказов. При этом можно считать, что интенсивность отказов элемента имеет две взаимоисключающие составляющие, т. е. l = l 1 + l 2 , где l 1 - постоянная интенсивность статистически независимых отказов элемента, l 2 - интенсивность множественных отказов резервированной системы или элемента. Поскольку a = l 2 / l , то l 2 = a/ l , и следовательно, l 1 =(1- a ) l .

Приведем формулы и зависимости для вероятности безотказной работы, интенсивности отказов и средней наработки на отказ в случае систем с параллельным и последовательным соединением элементов, а также систем с k исправными элементами из п и систем, элементы которых соединены по мостиковой схеме.

Система с параллельным соединением элементов (рис. 4.5.13) - обычная параллельная схема, к которой последовательно подсоединен один элемент. Параллельная часть (I) схемы отображает независимые отказы в любой системе из n элементов, а последовательно соединенный элемент (II) - все множественные отказы системы.

Рис. 4.5.13. Модифицированная система с параллельным соединением одинаковых элементов

Гипотетический элемент, характеризуемый определенной вероятностью появления множественного отказа, последовательно соединен с элементами, которые характеризуются независимыми отказами. Отказ гипотетического последовательно соединенного элемента (т.е. множественный отказ) приводит к отказу всей системы. Предполагается, что все множественные отказы полностью взаимосвязаны. Вероятность безотказной работы такой системы определяется как R р ={1-(1-R 1) n } R 2 , где n - число одинаковых элементов; R 1 - вероятность безотказной работы элементов, обусловленная независимыми отказами; R 2 - вероятность безотказной работы системы, обусловленная множественными отказами.

l 1 и l 2 выражение для вероятности безотказной работы принимает вид

R р (t)={1-(1-e -(1- a ) l t ) n }e - al t ,(4.5.28)
где t - время.

Влияние множественных отказов на надежность системы с параллельным соединением элементов наглядно демонстрируется с помощью рис. 4.5.14 – 4.5.16; при увеличении значения параметра a вероятность безотказной работы такой системы уменьшается.

Параметр a принимает значения от 0 до 1. При a = 0 модифицированная параллельная схема ведет себя как обычная параллельная схема, а при a =1 она действует как один элемент, т. е. все отказы системы являются множественными.

Поскольку интенсивность отказов и среднее время наработки на отказ любой системы можно определить с помощью (4.3 .7 ) и формул
,
,
с учетом выражения для R р (t ) получаем, что интенсивность отказов (рис. 4.5.17) и средняя наработка на отказ модифицированной системы соответственно равны
,(4.5.29)
,где .(4.5.30)

Рис. 4.5.14. Зависимость вероятности безотказной работы системы с параллельным соединением двух элементов от параметра a

Рис. 4.5.15. Зависимость вероятности безотказной работы системы с параллельным соединением трех элементов от параметра a

Рис. 4.5.16. Зависимость вероятности безотказной работы системы с параллельным соединением четырех элементов от параметра a

Рис. 4.5.17. Зависимость интенсивности отказов системы с параллельным соединением четырех элементов от параметра a

Пример 4.5.12. Требуется определить вероятность безотказной работы системы, состоящей из двух одинаковых параллельно соединенных элементов, если l =0,001 ч -1 ; a =0,071; t=200 ч.

Вероятность безотказной работы системы, состоящей из двух одинаковых параллельно соединенных элементов, для которой характерны множественные отказы, равна 0,95769. Вероятность безотказной работы системы, состоящей из двух параллельно соединенных элементов и характеризуемой только независимыми отказами, равна 0,96714.

Система с k исправными элементами из п одинаковых элементов включает в себя гипотетический элемент, соответствующий множественным отказам и соединенный последовательно с обычной системой типа k из n, для которой характерны независимые отказы. Отказ, отображаемый этим гипотетическим элементом, вызывает отказ всей системы. Вероятность безотказной работы модифицированной системы с k исправными элементами из n можно вычислить по формуле

,(4.5.31)

где R 1 - вероятность безотказной работы элемента, для которого характерны независимые отказы; R 2 - вероятность безотказной работы системы с k исправными элементами из n , для которой характерны множественные отказы.

При постоянных интенсивностях l 1 и l 2 полученное выражение принимает вид

.(4.5.32)

Зависимость вероятности безотказной работы от параметра a для систем с двумя исправными элементами из трех и двумя и тремя исправными элементами из четырех показаны на рис. 4.5.18 - 4.5.20. При увеличении параметра a вероятность безотказной работы системы уменьшается на небольшую величину (l t).

Рис. 4.5.18. Вероятность безотказной работы системы, сохраняющей работоспособность при отказе двух из n элементов

Рис. 4.5.19. Вероятность безотказной работы системы, сохраняющей работоспособность при отказе двух из четырех элементов

Рис. 4.5.20. Вероятность безотказной работы системы, сохраняющей работоспособность при отказе трех из четырех элементов

Интенсивность отказов системы с k исправными элементами из n и средняя наработка на отказ могут быть определены следующим образом:

,(4.5.33)

где h = {1-e -(1-b )l t },

q = e (r a -r- a ) l t

.(4.5.34)

Пример 4.5.13. Требуется определить вероятность безотказной работы системы с двумя исправными элементами из трех, если l =0,0005 ч - 1 ; a =0,3; t =200 ч.

С помощью выражения для R kn находим, что вероятность безотказной работы системы, в которой происходили множественные отказы, составляет 0,95772. Отметим, что для системы с независимыми отказами эта вероятность равна 0,97455.

Система с параллельно-последовательным соединением элементов соответствует системе, состоящей из одинаковых элементов, для которых характерны независимые отказы, и ряда ветвей, содержащих воображаемые элементы, для которых характерны множественные отказы. Вероятность безотказной работы модифицированной системы с параллельно-последовательным (смешанным) соединением элементов можно определить с помощью формулы R ps ={1 - (1-) n } R 2 , где m - число одинаковых элементов в ответвлении, n - число одинаковых ответвлений.

При постоянных интенсивностях отказов l 1 и l 2 это выражение принимает вид

R рs (t) = e - bl t . (4.5.39)

(здесь А=(1- a ) l ). Зависимость безотказной работы системы R b (t) для различных параметров a показана на рис. 4.5.21. При малых значениях l t вероятность безотказной работы системы с элементами, соединенными по мостиковой схеме, убывает с увеличением параметра a .

Рис. 4.5.21. Зависимость вероятности безотказной работы системы, элементы которой соединены по мостиковой схеме, от параметра a

Интенсивность отказов рассматриваемой системы и средняя наработка на отказ могут быть определены следующим образом:
l + .(4.5.41)

Пример 4.5.14. Требуется вычислить вероятность безотказной работы в течение 200 ч для системы с одинаковыми элементами, соединенными по мостиковой схеме, если l =0,0005 ч - 1 и a =0,3.

Используя выражение для R b (t), находим, что вероятность безотказной работы системы с соединением элементов по мостиковой схеме составляет примерно 0,96; для системы с независимыми отказами (т.е. при a =0) эта вероятность равна 0,984.

Модель надежности системы с множественными отказами

Для анализа надежности системы, состоящей из двух неодинаковых элементов, для которых характерны множественные отказы, рассмотрим такую модель, при построении которой были сделаны следующие допущения и приняты следующие обозначения:

Допущения (1) множественные отказы и отказы других типов статистически независимы; (2) множественные отказы связаны с выходом из строя не менее двух элементов; (3) при отказе одного из нагруженных резервированных элементов отказавший элемент восстанавливается, при отказе обоих элементов восстанавливается вся система; (4) интенсивность множественных отказов и интенсивность восстановлений постоянны.

Обозначения
P 0 (t) - вероятность того, что в момент времени t оба элемента функционируют;
P 1 (t) - вероятность того, что в момент времени t элемент 1 вышел из строя, а элемент 2 функционирует;
P 2 (t) - вероятность того, что в момент времени t эл мент 2 вышел из строя, а элемент 1 функционирует;
P 3 (t) - вероятность того, что в момент времени t элементы 1 и 2 вышли из строя;
P 4 (t) - вероятность того, что в момент времени t имеются специалисты и запасные элементы для восстановления обоих элементов;
a - постоянный коэффициент, характеризующий наличие специалистов и запасных элементов;
b - постоянная интенсивность множественных отказов;
t - время.

Рассмотрим три возможных случая восстановления элементов при их одновременном отказе:

Случай 1. Запасные элементы, ремонтный инструмент и квалифицированные специалисты имеются для восстановления обоих элементов, т. е. элементы могут быть восстановлены одновременно .

Случай 2. Запасные элементы, ремонтный инструмент и квалифицированные специалисты имеются только для восстановления одного элемента, т. е. может быть восстановлен только один элемент.

Случай 3 . Запасные элементы, ремонтный инструмент и квалифицированные специалисты отсутствуют, и, кроме того, может существовать очередь на ремонтное обслуживание.

Математическая модель системы, изображенной на рис. 4.5.22, представляет собой следующую систему дифференциальных уравнений первого порядка:

P" 0 (t) = - ,
P" 1 (t) = -(l 2 + m 1 )P 1 (t)+P 3 (t)

Рис. 4.5.22. Модель готовности системы в случае множественных отказов

Приравнивая в полученных уравнениях производные по времени нулю, для установившегося режима получаем

- ,
-(l 2 + m 1 )P 1 +P 3 m 2 +P 0 l 1 = 0,

-(l 1 + m 2 )P 2 +P 0 l 2 +P 3 m 1 = 0,

P 2 = ,

P 3 = ,

P 4 = .

Стационарный коэффициент готовности может быть вычислен по формуле

Опасность технических систем. Отказ, вероятность отказа.

Определение опасности

Опасность - центральное понятие как сферы безопасности жизнедеяте­льности в техносфере, так и промышленной безопасности. Под опасностью понимаются явления, процессы, объекты, способные в определенных усло­виях наносить вред здоровью человека, ущерб окружающей природной сре­де и социально-экономической инфраструктуре, т. е. вызывать нежелатель­ные последствия непосредственно или косвенно. Другими словами, опасность - следствие действия некоторых негативных (вредных и опас­ных) факторов на определенный объект (предмет) воздействия. При несоот­ветствии характеристик воздействующих факторов характеристикам объ­екта (предмета) воздействия и появляется феномен опасности (например, ударная волна, аномальная температура, недостаток кислорода в воздухе, токсичные примеси в воздухе и т. п.).

Опасность - свойство, внутренне присущее сложной технической сис­теме. Она может реализоваться в виде прямого или косвенного ущерба для объекта (предмета) воздействия постепенно или внезапно и резко - в ре­зультате отказа системы. Скрытая (потенциальная) опасность для человека реализуется в форме травм, которые происходят при несчастных случаях, авариях, пожарах и пр., для технических систем - в форме разрушений, по­тери управляемости и т. д., для экологических систем - в виде загрязнений, утрате видового разнообразия и др.

Определяющие признаки - возможность непосредственного отрицате­льного воздействия на объект (предмет) воздействия; возможность наруше­ния нормального состояния элементов производственного процесса, в резу­льтате которого могут возникнуть аварии, взрывы, пожары, травмы. Наличие хотя бы одного из указанных признаков является достаточным для отнесения факторов к разделу опасных или вредных.

Количество признаков, характеризующих опасность, может быть увели­чено или уменьшено в зависимости от целей анализа.

Анализ реальных аварийных ситуаций, событий и факторов и человече­ская практика уже сегодня позволяют сформулировать ряд аксиом об опас­ности технических систем:

Аксиома 1. Любая техническая система потенциально опасна. Потенциа­льность опасности заключается вскрытом, неявном характере и проявляет­ся при определенных условиях. Ни один вид технической системы при ее функционировании не обеспечивает абсолютной безопасности.

Аксиома 2. Техногенные опасности существуют, если повседневные пото­ки вещества, энергии и информации в техносфере превышают пороговые значе­ния. Пороговые, или предельно допустимые, значения опасностей устанав­ливаются из условия сохранения функциональной и структурной целостности человека и природной среды. Соблюдение предельно допусти­мых значений потоков создает безопасные условия жизнедеятельности че­ловека в жизненном пространстве и исключает негативное влияние техно­сферы на природную среду.

Аксиома 3. Источниками техногенных опасностей являются элементы техносферы. Опасности возникают при наличии дефектов и иных неисправ­ностей в технических системах, при неправильном использовании техниче­ских систем. Технические неисправности и нарушения режимов использо­вания технических систем приводят, как правило, к возникновению травмоопасных ситуаций, а выделение отходов (выбросы в атмосферу, сто­ки в гидросферу, поступление твердых веществ на земную поверхность, энергетические излучения и поля) сопровождается формированием вред­ных воздействий на человека, природную среду и элементы техносферы.

Аксиома 4. Техногенные опасности действуют в пространстве и во време­ни. Травмоопасные воздействия действуют, как правило, кратковременно и спонтанно в ограниченном пространстве. Они возникают при авариях и катастрофах, при взрывах и внезапных разрушениях зданий и сооруже­ний. Зоны влияния таких негативных воздействий, как правило, ограниче­ны, хотя возможно распространение их влияния и на значительные терри­тории, например, при аварии на ЧАЭС.

Для вредных воздействий характерно длительное или периодическое не­гативное влияние на человека, природную среду и элементы техносферы. Пространственные зоны вредных воздействий изменяются в широких пре­делах от рабочих и бытовых зон до размеров всего земного пространства. К последним относятся воздействия выбросов парниковых и озоноразруша- ющих газов, поступление радиоактивных веществ в атмосферу и т. п.

Аксиома 5. Техногенные опасности оказывают негативное воздействие на человека, природную среду и элементы техносферы одновременно. Человек и окружающая его техносфера, находясь в непрерывном материальном, энергетическом и информационном обмене, образуют постоянно действую­щую пространственную систему «человек - техносфера». Одновременно существует и система «техносфера - природная среда». Техногенные опас­ности не действуют избирательно, они негативно воздействуют на все со­ставляющие вышеупомянутых систем одновременно, если последние ока­зываются в зоне влияния опасностей.

Аксиома 6. Техногенные опасности ухудшают здоровье людей, приводят к травмам, материальным потерям и к деградации природной среды.

1.2 Определение надёжности. Отказ, вероятность отказа.

Работа любой технической системы может характеризоваться ее эффек­тивностью, под которой понимается совокупность свойств, определяющих способность системы успешно выполнять определенные за­дачи.

В соответствии с ГОСТ 27.002-89 под надежностью понимают свойство объекта сохранять во времени в установленных пределах значения всех па­раметров, характеризующих способность выполнять требуемые функции в заданных режимах и условиях применения, технического обслуживания, ремонтов, хранения и транспортировки.

Надежность в общем случае - комплексное свойство, включающее та­кие понятия, как безотказность, долговечность, ремонтопригодность, сохра­няемость. Для конкретных объектов и условий их эксплуатации эти свойст­ва могут иметь различную относительную значимость.

Безотказность - свойство объекта непрерывно сохранять работоспособ­ность в течение некоторой наработки или в течение некоторого времени.

Отказ объекта - событие, заключающееся в том, что объект полностью или частично перестает выполнять заданные функции. При полной потере работоспособности возникает полный отказ, при частичной - частичный. Понятия полного и частичного отказов каждый раз должны быть четко сформулированы перед анализом надежности, поскольку от этого зависит количественная оценка надежности.

Причины возникновения отказов происходят из-за:

Конструктивных дефектов;

Технологических дефектов;

Эксплуатационных дефектов;

Постепенного старения (износа).

Наработка до отказа - вероятность того, что в пределах заданной нара­ботки отказ объекта не возникнет (при условии работоспособности в нача­льный момент времени).

Для режимов хранения и транспортировки может применяться анало­гично определяемый термин «вероятность возникновения отказа».

Средняя наработка до отказа - математическое ожидание случайной наработки объекта до первого отказа.

Средняя наработка между отказами - математическое ожидание слу­чайной наработки объекта между отказами.

Обычно этот показатель относится к установившемуся процессу эксплу­атации. В принципе средняя наработка между отказами объектов, состоя­щих из стареющих во времени элементов, зависит от номера предыдущего отказа. Однако с ростом номера отказа (т.е. с увеличением длительности эксплуатации) эта величина стремится к некоторой постоянной, или, как говорят, к своему стационарному значению.

Средняя наработка на отказ - отношение наработки восстанавливае­мого объекта за некоторый период времени к математическому ожиданию числа отказов в течение этой наработки.

Этим термином можно назвать кратко среднюю наработку до отказа и среднюю наработку между отказами, когда оба показателя совпадают.

Интенсивность отказов - условная плотность вероятности отказа невосстанавливаемого объекта, определяемая для рассматриваемого момен­та времени при условии, что до этого момента отказ не возник.

Параметр потока отказов - плотность вероятности возникновения отказа восстанавливаемого объекта, определяемая для рассматриваемого момента времени.

Параметр потока отказа может быть определен как отношение числа отказов объекта за определенный интервал времени к длительности этого интервала при ординарном потоке отказов.

Вероятностью безотказной работы Р(t) называется вероятность того, что при определённых условиях эксплуатации в заданном интервале времени или в пределах заданной наработки не произойдёт ни одного отказа:

Так как безотказная работа и отказ являются событиями несовместными и противоположными, то между ними справедливо следующее соотношение:

Так как Q (t) есть закон распределения случайной величины (отказов), то зависимость между возможными значениями непрерывной случайной вели­чины Т и вероятностями попадания в их окрестность называется её плотностью вероятности.

Частота отказов a(t) есть плотность вероятности времени работы изделия до первого отказа:

Интенсивностью отказов называется отношение числа отказавших изде­лий в единицу времени к среднему числу изделий, исправно работающих в данный отрезок времени. Вероятностная оценка этой характеристики находится из выражения:

Средней наработкой до первого отказа называется математическое ожидание М[t] времени работы изделия до отказа. Как математическое ожидание, Т ср вычисляется через частоту отказов (плотность распределения времени безотказной работы):

так как t > 0 и P(0) = 1, а P(∞) = 0, то

Зная один из показателей надёжности и закон распределения отказов, можно вычислить остальные характеристики надёжности с учётом следую­щих формул:

Опыт взаимодействия человека с техническими системами позволяет идентифицировать травмиру­ющие и вредные факторы, а также выработать ме­тоды оценки вероятности появления опасных ситуа­ций. Прежде всего, это накопление статистических данных об аварийности и травматизме (табл. 1), различные способы преобразования и обработки ста­тистических данных, повышающие их информатив­ность. Недостатком этого метода является его ог­раниченность, невозможность экспериментирования и неприменимость к оценке опасности новых тех­нических средств и технологий.

Значительное развитие и практическое примене­ние получила теория надежности. Надежность это свойство объекта сохранять во времени в уста­новленных пределах значения всех параметров, по­зволяющих выполнять требуемые функции. Для ко­личественной оценки надежности применяют веро­ятностные величины.

Таблица 1

Ши­рокое распространение получила диаграмма ветвя­щейся структуры, называемая «дерево событий». Рассмотрим процедуру построения дерева, его качественный и количественный анализ на приме­ре.

Будем считать, что для гибели человека от элек­трического тока необходимо и достаточно включе­ние его тела в цепь, обеспечивающую прохождение смертельного тока. Следовательно, чтобы произо­шел несчастный случай (событие А), необходимо од­новременное выполнение по крайней мере трех ус­ловий: наличие потенциала высокого напряжения на металлическом корпусе электроустановки (собы­тие Б), появление человека на заземленном прово­дящем основании (событие В), касание человека корпуса электроустановки (событие Г).

В свою очередь событие Б может быть следстви­ем любого из событий - предпосылок Д и Е, на­пример, нарушение изоляции или смещение неизо­лированного контакта и касание им корпуса. Собы­тие В может появиться как результат предпосылок Ж и З, когда человек становится на заземленное проводящее основание или касается телом зазем­ленных элементов помещения. Событие Г может явиться одной из трех предпосылок И, К и Л - ремонт, техобслуживание или работа установки.

Анализ дерева событий состоит в выявлении ус­ловий, минимально необходимых и достаточных для возникновения или невозникновения головного со­бытия. Модель может давать несколько минималь­ных сочетаний исходных событий, приводящих в совокупности к данному происшествию. В данном примере имеются двенадцать минимальных аварий­ных сочетаний: ДЖИ, ДЖК, ДЖЛ, ДЗИ, ДЗК, ДЗЛ, ЕЖИ, ЕЖК, ЕЖЛ,ЕЗИ, ЕЗК, ЕЗЛ и три минималь­ных секущих сочетания, исключающих возможность появления происшествия при одновременном отсут­ствии образующих их событий: ДЕ, ЖЗ, ИКЛ.

Аналитическое выражение условий появления исследуемого происшествия имеет вид А = (Д + Е) (Ж + 3)(И + К + Л). Подставив вместо буквенных символов вероятности соответствующих предпосы­лок, можно получить оценку риска гибели челове­ка от электрического тока в конкретных условиях.

Например, при равных вероятностях Р(Д) = Р(Е) = = ...Р(Л) = 0,1 вероятность гибели человека от элек­трического тока в рассматриваемом случае

Р(А)=(ОД+0,1)(0,1+ОД)(0,1+о,1+ОД)=0,012.

Таким образом, может быть рассчитана вероятность несчастного случая или аварии на производстве.

Анализ причин появления опасности для человека при его взаимодействии с техническими системами позволяет выделить причины - организационные и технические. Для устранения организационных причин совершенствуется технологический процесс, уточняются процедуры подготовки и контроля операторов. При этом техническая система рассматривается как замкнутая система, взаимодействующая с окружающей средой. В этом случае под окружающей средой понимается комплекс условий на каждом этапе жизненного цикла системы. В комплекс условий включаются все возможные факторы, воздействующие на систему, в том числе профессионализм конструкторов, технологические факторы производственного процесса изготовления, режимы эксплуатации (электрические, тепловые и др.). Объективной закономерностью является то, что при переходе от этапа к этапу в жизненном цикле технической системы количество воздействующих на систему факторов возрастает, увеличивается в связи с этим и степень жесткости влияния. Это ведет к уменьшению надежности и увеличению опасности в цепочке «человек - техническая система - окружающая среда», что делает задачу обеспечения безопасности технических систем чрезвычайно сложной.

На практике необходимый уровень безопасности технических средств и технологических процессов устанавливается системой государственных стандартов безопасности труда (ССБТ) с помощью соответствующих показателей. Стандарты формируют общие требования безопасности, а также требования безопасности к различным группам оборудования, производственных процессов, требования к средствам обеспечения безопасности труда.

Нормативные показатели безопасности во всех сферах труда разрабатываются в соответствии с санитарными нормами и вводятся посредством соответствующих государственных стандартов (ГОСТ). Так, например, внедрение новой техники увеличило интенсивность шума и вибрации и расширило диапазон частот в ультра и инфразвуковых частях спектра колебаний. Это вызвало необходимость разработки и включения в ГОСТ нормативов допустимых уровней ультра и инфразвука на производстве.

Соответствующие нормативы, гарантирующие безопасное взаимодействие человека с техническими системами и технологическими процессами, установлены для электромагнитных полей, электрического напряжения и тока, излучений оптического диапазона, ионизирующих излучений, химических, биологических и психофизических опасных и вредных факторов. При разработке технических средств и технологий применяются все возможные меры для снижения опасных и вредных факторов ниже предельно допустимого уровня. Для каждого технического средства разрабатываются правила эксплуатации, гарантирующие безопасность при их выполнении. Для каждой технологической операции также разрабатываются правила техники безопасности.
2 Качественный и количественный анализ опасностей

Качественный анализ опасностей

Качественные методы анализа опасностей включают:

Предварительный анализ опасностей;

Анализ последствий отказов;

Анализ опасностей с помощью "дерева причин";

Анализ опасностей методом потенциальных отклонений;

Анализ ошибок персонала;

Причинно-следственный анализ.

В результате анализа аварийной (потенциальной) опасности могут быть определены следующие показатели:

Индивидуальный риск;

Социальный риск;

Структура поражённых по степени тяжести;

Вид поражений;

Материальный ущерб и др.

Наиболее распространённым методом анализа безопасности является метод построения "деревьев отказов (ошибок)". В терминологии теории построения и анализа "деревьев отказов " выход из строя определённых элементов, например, нарушение герметичности резервуара со сжиженным углеводородным газом с последующим образованием облака топливовоздушной смеси и его взрывом, классифицируется как внешнее нежелательное событие (ВНС).

В строящихся деревьях, как правило, имеются ветви опасностей. Многоэтажный процесс ветвления "дерева" требует введения ограничений с целью определения его пределов. Логические операции принято обозначать соответствующими символами (см.табл.2).

Таблица 2 - Символы событий

Построение "дерева причин", "дерева отказов" является эффективной процедурой выявлении причин различных нежелательных событий (аварий, травм, пожаров, дорожно-транспортных происшествий) и экспертизой безопасности оборудования и процессов.

Рисунок 2

А - отказ средств борьбы со взрывом; Б - образование облака ТВС; В - разгерметизация каждой ёмкости; Г - инициация взрыва; Д- факел, печь; Е - автотранспорт; З - электродвигатель; Ж - огневые работы; И - удар предмета; К - разрушение резервуара; Л - разрушение трубопровода; М -разгерметизация арматуры; Н - температура; О - скорость ветра; П - состояние атмосферы.

Любой объект техносферы потенциально опасен. Всегда существует возможность происшествия: инцидента, аварии, катастрофы.

Инцидент – событие в результате которого возникает или может возникнуть несчастный случай.

Аварией обычно считается происшествие, в результате которого повреждена техника, без гибели людей.

Крупная авария , повлекшая за собой человеческие жертвы, значительный материальный ущерб, загрязнение среды, считается катастрофой.

Причины происшествия могут быть внутренними (отказы техники, ошибочные действия персонала) и внешними (транспортные аварии при перевозке опасных грузов, противоправные действия, природная среда и др.).

Опасность объекта техносферы – это его свойство, состоящее в возможности в процессе эксплуатации при определенных обстоятельствах причинять ущерб человеку, организации, окружающей природной среде .

Экономический ущерб , который может быть причинен объектом, называется потенциалом угрозы . Верхний предел потенциала угрозы обозначается как потенциал опасности технического объекта.

По потенциалу опасности промышленные объекты делятся на неопасные и опасные. Ущерб от потенциально опасных объектов наступает в случае аварий. В РФ опасные объекты подлежат регистрации в государственном реестре, обязаны заниматься декларированием безопасности и страхованием ответственности за ущерб третьим лицам.

Согласно федеральному закону «О промышленной безопасности опасных производственных объектов» от 21.07.97 № 116-ФЗ различают 5 групп объектов по виду опасности :

1) опасные вещества (воспламеняющиеся, окисляющиеся, горючие, взрывчатые, токсичные);

2) давление (более 0,07 МПа), температура нагрева воды (более115°);

3) высота (грузоподъемные машины, эскалаторы, фуникулеры, канатные дороги);

4) расплавы черных и цветных металлов;

5) подземные условия (горные работы).

По природе образующихся в результате аварии опасных факторов выделяют 6 групп потенциально-опасных объектов:

1) ядерно и радиационно опасные;

2) химически опасные;

3) пожаровзрывоопасные;

4) биологически опасные;

5) гидродинамически опасные;

6) объекты жизнеобеспечения.

Различают следующие виды опасных техногенных явлений: транспортные аварии, пожары, взрывы, химические аварии, радиационные аварии, гидродинамические аварии, разрушения зданий.

Оценка риска аварии – процесс, используемый для определения вероятности (или частоты) и степени тяжести последствий реализации опасностей аварий для здоровья человека, имущества и (или) окружающей природной среды.

Оценка риска включает анализ вероятности (или частоты), анализ последствий и их сочетания.

Риск аварии – мера опасности, характеризующая возможность возникновения аварии на опасном производственном объекте и тяжесть ее последствий.

Основными количественными показателями риска аварии являются:

· технический риск – вероятность отказа технических устройств с последствиями определенного уровня (класса) за определенный период функционирования опасного производственного объекта (определяют методами теории надежности);

· индивидуальный риск – частота поражения отдельного человека в результате воздействия исследуемых факторов опасности аварий. Рекомендуется оценивать индивидуальный риск отдельно для персонала объекта и для населения прилегающей территории или, при необходимости, для более узких групп, например, для рабочих различных специальностей;

· потенциальный территориальный риск (или потенциальный риск) – частота реализации поражающих факторов аварии в рассматриваемой точке территории;

· коллективный риск – ожидаемое количество пораженных в результате возможных аварий за определенное время;

· социальный риск , или F/N-кривая (в зарубежных работах – кривая Фармера), – зависимость частоты возникновения событий (F), в которых пострадало на определенном уровне не менее N человек, от этого числа N. Характеризует тяжесть последствий (катастрофичность) реализации опасностей. Под N можно понимать и общее число пострадавших, и число смертельно травмированных или другой показатель тяжести последствий. Критерий приемлемого риска будет определяться не числом для отдельного события, а кривой, построенной для различных сценариев аварии с учетом их вероятности. В настоящее время общераспространенным подходом для определения приемлемости риска является использование двух кривых, когда, например, в логарифмических координатах определены F/N-кривые приемлемого и неприемлемого риска смертельного травмирования. Область между этими кривыми определяет промежуточную степень риска, вопрос о снижении которой следует решать, исходя из специфики производства и региональных условий;

· ущерб от аварии – потери (убытки) в производственной и непроизводственной сфере жизнедеятельности человека, вред окружающей природной среде, нанесенные в результате аварии на опасном производственном объекте и исчисляемые в денежном эквиваленте».

Надежность технической системы – это вероятность сохранения системой работоспособности в течение определенного времени.

Техническую систему будем представлять в виде сложной системы следующей иерархии:

Техническая система;

Устройства;

Элементы.

Техническая система –– совокупность взаимосвязанных элементов (объектов, устройств), обеспечивающих выполнение конкретных практических задач.

Устройством называется законченная конструкция, которая, являясь частью системы, имеет самостоятельное целевое назначение.

Элементы –– это части системы или устройства, которые выполняют в нем определенные функции и не могут иметь самостоятельного (вне связи с другими элементами или устройствами) применение.

Причины недостаточной надежности технических систем.

Основными причинами отказов технических систем являются внезапные (случайные) отказы, отказы вследствие ухудшения характеристик элементов (старение, износ), а также по вине скрытых производственных дефектов, характерных для начального периода эксплуатации, или нарушения условий эксплуатации .

Возрастание интенсивности отказов технических систем связано, как правило, с ужесточением условий их функционирования (эксплуатации) и с недостаточной квалификацией обслуживающего персонала.

В целом, все причины, приводящие к снижению надежности технических систем, можно разделить на следующие: конструктивные, производственные, эксплуатационные, организационные.

Конструктивные причины : низкая надежность элементной базы, неправильный выбор элементов, неудачное схемно-компоновочное решение, недостаточная унификация элементов, недостаточная отработка технологий на этапах испытаний.

Производственные причины : нарушение качества материалов, недостаточный контроль входных параметров, недостаточная отработка технологии производства и сборки устройств, общая низкая культура производства.

Эксплуатационные причины : низкая квалификация технического персонала, низкая эффективность контрольно––проверочной аппаратуры, нарушение условий эксплуатации.

Организационные причины : отсутствие требований по поддержанию заданных показателей надежности, несоответствие заводских испытаний реальным условиям эксплуатации, неритмичность эксплуатации.

1.3. Цена надежности.

Стоимость технической системы, как правило, определяется стоимостью ее создания (строительства) и стоимостью эксплуатации системы и зависит от надежности системы.

C å (Р)= C 0 (Р) + C э (Р)

где C å (Р) –– общая стоимость технической системы;

C 0 (Р) –– стоимость создания технической системы;

C э (Р) –– стоимость эксплуатации технической системы,

Р –– надежность системы.

Расходы, связанные с созданием технической системы, являются функцией требований к ее надежности. Чем выше требования к надежности системы, тем выше ее стоимость, т.е. функция C 0 (P) –– неубывающая функция надежности системы (рис. 1.2).

Стоимость эксплуатации технической системы также зависит от ее надежности, но в обратной зависимости. Чем выше надежность системы, тем ниже стоимость ее эксплуатации, чем надежнее создана система, тем меньше средств она требует на свое поддержание в исправном состоянии.

Рис. 1.2. Изменение суммарной стоимости С S системы в зависимости от ее надежности Р, затрат на эксплуатацию С Э и создание системы С 0 .

Рациональное распределение средств на повышение надежности технических систем на этапе проектирования, изготовления, испытания и эксплуатации может привести к существенной экономии суммарных расходов обеспечения функционирования системы. Зачастую распределение средств, выделенных на снижение негативных последствий аварий принимается за управление риском.

Анализ опасностей и оценки риска аварий на опасных производственных объектах (далее - анализ риска аварий) представляют собой совокупность научно-технических методов исследования опасностей возникновения, развития и последствий возможных аварий , включающую планирование работ, идентификацию опасностей аварий, оценку риска аварий, установление степени опасности возможных аварий, а также разработку и своевременную корректировку мероприятий по снижению риска аварий .

Риск аварии - мера опасности, характеризующая возможность возникновения аварии на опасном производственном объекте и соответствующую ей тяжесть последствий. В анализе риска аварий в качестве основных количественных показателей опасности (показателей риска) рекомендуется использовать:

технический риск - вероятность отказа технических устройств с последствиями определенного уровня (класса) за определенный период функционирования опасного производственного объекта;

индивидуальный риск - ожидаемая частота (частота) поражения отдельного человека в результате воздействия исследуемых поражающих факторов аварии;

потенциальный территориальный риск (или потенциальный риск)– частота реализации поражающих факторов аварии в рассматриваемой точке на площадке опасного производственного объекта и прилегающей территории;

коллективный риск (или ожидаемые людские потери) - ожидаемое

количество пораженных в результате возможных аварий за определенный период времени;

социальный риск (или риск поражения группы людей) - зависимость частоты возникновения сценариев аварий F, в которых пострадало на определенном уровне не менее N человек, от этого числа N. Характеризует социальную тяжесть последствий (катастрофичность) реализации совокупности сценариев аварии и представляется в виде соответствующей F/N-кривой;

ожидаемый ущерб - математическое ожидание величины ущерба от возможной аварии за определенный период времени;

материальный риск (или риск материальных потерь) - зависимость частоты возникновения сценариев аварий F, в которых причинен ущерб на определенном уровне потерь не менее G, от количества этих потерь G. Характеризует экономическую тяжесть последствий реализации опасностей аварий и представляется в виде соответствующей F/G-кривой.

Для некоторых ситуаций определены допустимые значения риска, например, нормативные значения пожарного риска для производственных объектов.

Величина индивидуального пожарного риска в зданиях, сооружениях, строениях и на территориях производственных объектов не должна превышать одну миллионную в год.

Риск гибели людей в результате воздействия опасных факторов пожара должен определяться с учетом функционирования систем обеспечения пожарной безопасности зданий, сооружений и строений.

Для производственных объектов, на которых обеспечение величины индивидуального пожарного риска одной миллионной в год невозможно в связи со спецификой функционирования технологических процессов, допускается увеличение индивидуального пожарного риска до одной десятитысячной в год. При этом должны быть предусмотрены меры по обучению персонала действиям при пожаре и по социальной защите работников, компенсирующие их работу в условиях повышенного риска.

Величина индивидуального пожарного риска в результате воздействия опасных факторов пожара на производственном объекте для людей, находящихся в селитебной зоне вблизи объекта, не должна превышать одну стомиллионную в год.

Величина социального пожарного риска воздействия опасных факторов пожара на производственном объекте для людей, находящихся в селитебной зоне вблизи объекта, не должна превышать одну десятимиллионную в год.

Наиболее распространенной методикой количественной оценки риска является мультипликативная форма его представления :

где: - значение риска;

– вероятность возникновения хотя бы одной аварии за рассматриваемый период работы объекта или технической системы;

Ожидаемый ущерб при возникновении аварии;

– средства, выделяемые на снижение риска;

– средства, выделяемые на снижение вероятности реализации аварии;

– средства, выделяемые на снижение ожидаемого ущерба в случае возникновения аварии.

Основной принцип управления риском состоит в приоритетном максимальном снижении вероятности возникновения аварии, и во вторую очередь – забота о сокращении ожидаемого ущерба.

Возникновение аварии очевидным образом непосредственно связано с надежностью технических систем, ее устройств или элементов.

Следует заметить, что вероятность возникновения хотя бы одной аварии имеет существенное отличие от вероятности возникновения ровно одной аварии.

2. Основные понятия теории надежности.

Теория надежности в большинстве случаев оперирует случайными величинами, поэтому большая часть понятий и определений связана с понятийным аппаратом теории вероятностей.

Отказ –– полная или частичная утрата работоспособности элементом, устройством или рассматриваемой технической системой.

Исправность –– состояние системы, при котором она в данный момент времени соответствует всем требованиям, установленным в отношении как основных параметров, так и «второстепенных».

Работоспособность –– состояние системы, при котором она в данный момент времени соответствует всем требованиям, установленным в отношении ее основных параметров.

Безотказность –– свойство системы сохранять работоспособность в течение заданного интервала времени в определенных условиях эксплуатации.

Неисправность –– состояние системы, при котором она в данный момент времени не соответствует хотя бы одному из требований, установленных в отношении как основных параметров, так и «второстепенных».

2.1. Основные количественные характеристики надежности и связь между ними.

Основной количественной характеристикой надежности является вероятность безотказной работы, определяемая как вероятность P(t) нахождения системы в исправном состоянии в течение времени T ³ t, где Т –– случайная величина продолжительности работы системы до отказа, t –– детерминированная величина текущего времени или его конкретное значение :

(2.1)
где W(T) –– вероятность реализации события, заключающегося в том, что отказ системы не произойдет ранее t.

Функция P(t) обладает следующими свойствами:

P(0) = 1, P(¥) = 0, P(t 2) £ P(t 1) при t 2 > t 1 ,

т.е. функция P(t) –– невозрастающая функция времени.

Эту функцию часто называют функцией надежности или просто надежностью технической системы.

Характеристикой, противоположной надежности, является вероятность отказа Q(t), как вероятность того, что устройство или техническая система откажет в течение времени T < t:

(2.2)
Свойства функции Q(t) : Q(0) = 0, Q(¥) = 1, Q(t 2) ³ Q(t 1) при t 2 > t 1 .

Таким образом, функция ненадежности Q(t) представляет собой функцию распределения времени исправной работы системы F(t) (П 2.20) Приложения 2.

Очевидно: Q(t) + P(t) = 1; Q(t)=1-P(t);

t

Q(t)

P(t)

(2.3)
т.е. система может находиться либо в исправном, либо в неисправном состоянии.

Рис. 2.1. Характер изменения функций надежности и отказов во времени.

Плотность f(t) распределения времени работы системы до отказа согласно определению плотности из классической теории вероятностей имеет вид (П 2.21) приложения 2:

(2.4)
Наряду с аналитическими методами определения различных параметров надежности широко используются статистические методы, с помощью которых определяются так называемые статистические характеристики надежности систем. Эти характеристики представляют собой результаты обработки экспериментальных данных или данных прямых наблюдений.

Поскольку в эксперименте невозможно произвести наблюдения при t®¥ или на бесконечно малом временном интервале Dt®0, а также при бесконечно большом числе испытываемых систем, то статистические характеристики следует рассматривать как оценочные или приближенные к теоретическим. Статистическая плотность отказов в теории надежности определяют в виде отношения:

(2.5)
где Dn(t, Dt) –– число отказавших элементов на интервале Dt (от момента времени t до момента t+Dt); N 0 –– общее число элементов, выставленных на испытания; Dt ––интервал времени проведения испытаний или наблюдений.

n(t,∆t)-n(t)=∆n

Вышедшие из строя

n(t)

t

n(t,∆t)

t+∆t

T

Рис. 2.2. Порядок вычисления количества отказавших элементов Dn.

Статистически функция распределения времени исправной работы Q*(t) системы оценивается как отношение числа устройств n(t), вышедших из строя за время от начала испытаний до некоторого момента t, к общему числу элементов или устройств N 0 , поставленных на испытание:

(2.6)
Величину Q*(t) называют частостью отказов, которая является оценкой функции распределения отказов или вероятности отказа.

Очевидно, что чем больше проведено число независимых испытаний, тем ближе величина частости к соответствующему значению вероятности. В теории вероятностей такой характер приближения одних величин к другим чрезвычайно употребителен и для его описания введен специальный термин –– сходимость по вероятности.

Согласно первой предельной теореме (закон больших чисел), последовательность случайных величин x n сходится по вероятности к величине J, если при сколь угодно малом e>0 вероятность неравенства

½x n - J½< e с увеличением n неограниченно приближается к единице.

Таким образом, можно утверждать, что с увеличением числа опытов частость события сходится к его вероятности по вероятности.

Одной из характеристик надежности технических систем является частота отказов, в дальнейшем обозначаемая a(t).

Частотой отказов а * (t) называется отношение числа отказавших образцов в единицу времени к числу образцов, первоначально установленных на испытание при условии, что отказавшие образцы не восстанавливаются и не заменяются исправными.

(2.7)
Выражение (2.9) является статистическим определением частоты отказов. Вместе с тем этой характеристике можно придать вероятностное определение.

Число отказавших образцов на интервале Dt может быть определено по формуле:

(2.8)
где N(t) –– число образцов, исправно работающих к моменту t; N(t+Dt) –– число образцов, исправно работающих к моменту t+Dt.

При достаточно большом числе образцов N 0 справедливы следующие соотношения:

(2.9)

где P(t)- вероятность сохранения работоспособности исследуемым элементом до момента t, т.е. надежность рассматриваемого элемента.
Подставляя (2.8) в (2.9) и учитывая (2.7), получим:

(2.10)
Устремляя к нулевому пределу интервал Dt, получим с учетом определения Dn(t, Dt), принятого в (2.5):

(2.11)
или

(2.12)
Последнее равенство подтверждает идентичность (2.7), (2.10).

Из выражения (2.12) следует утверждение, что частота отказов представляет собой плотность распределения времени работы системы до ее отказа.

Наиболее употребительной в теории надежности является такая характеристика, как интенсивность отказов :

(2.13)
т.е. l(t) является условной плотностью распределения вероятности исправной работы системы, вычисленной при условии, что к моменту t система была исправна.

Статистической интерпретацией интенсивности отказов l(t)* является отношение числа однотипных устройств Dn(Dt) , вышедших из строя в интервале времени Dt , к числу устройств N(t) из общего числа N 0 , поставленных на испытания, продолжающих к моменту времени t оставаться исправными, умноженному на длину интервала Dt, при условии, что отказавшие образцы не восстанавливаются и не заменяются исправными.

(2.14)
Разделив числитель и знаменатель (2.14) на N 0 , получим:

,

(2.15)
Таким образом, интенсивность отказов определяется как отношение частоты отказов к статистической оценке вероятности P(t) исправной работы рассматриваемого элемента или устройства

(2.16)

P(t)= (2.17)

Типичная кривая изменения интенсивности отказов технических систем представлена на рисунке 2.3.

l(t)

t 2

t

Рис. 2.3. Типовая зависимость интенсивности отказов технических систем от времени.

Как видно из рис. 2.3 кривая l(t) имеет три характерных участка. Первый участок (от 0 до t 1) –– участок приработки, второй участок (от t 1 до t 2) –– участок нормальной эксплуатации системы, третий временной интервал (от t 2 и далее) –– участок старения системы. Здесь уместно отметить, что в период нормальной работы системы (от t 1 до t 2), как правило, интенсивность отказов не зависит от времени, l=const.

Выражение (2.16) с учетом (2.11) приобретает вид:

(2.18)
Интегрируя (2.18) при начальном условии P(0) = 1, получим:

Поскольку т.е. С=1.

Нижний предел интегрирования равен 0, т.к. отсчет времени производится от момента включения системы в работу.

Верхний предел определяется аргументом функции P(t) т.е. значением аргумента t.

Окончательно получаем:

(2.19)
Выражение (2.19) определяет вероятность безотказной работы технических систем и является одним из основных в теории надежности.

Среднее время до отказа технической системы Т С определяется как его математическое ожидание с нижним пределом интегрирования, равным нулю, поскольку время не имеет отрицательных значений:

(2.20)
Статистической интерпретацией среднего времени до первого отказа является среднее арифметическое значение времени работы устройства до ее первого отказа:

(2.21)
где R i –– называют частостью времени отказов t i ; t i –– время работы i-го элемента до первого отказа; N 0 –– число элементов, поставленных на испытание.

Выражение (2.20) можно представить в ином виде, подставив выражение плотности f(t) согласно (2.4):

(2.22)
Произведем интегрирование по частям:

(2.23)
Дисперсия D[T] случайного времени Т безотказной работы системы:

(2.24)
где Т –– случайное время безотказной работы системы; Т с –– математическое ожидание времени работы системы до отказа; f(x) –– дифференциальный закон распределения случайного времени безотказной работы системы.

Для случая распределения случайной величины по закону Пуассона при постоянном значении интенсивности отказов :

(2.25)
Тогда формула для вычисления дисперсии D[T] может быть выведена на основании следующих преобразований с учетом того, что время не может быть отрицательным, т.е. в (2.22) нижний предел интегрирования равен нулю:

(2.26)
Для вычисления этих интегралов вводятся обозначения:

тогда, пользуясь формулой интегрирования по частям, имеем:

Подставляя это выражение в (2.24), получим:

(2.27)

2.2. Характеристики технических систем, используемые в теории надежности.

Коэффициент стабильности надежности K ст –– отношение значений вероятностей исправной работы устройства для двух произвольных периодов времени;

K ст (2.28)
Если коэффициент стабильности равен единице, то надежность системы на участке t остается неизменной.

На практике часто используется показатель изменения надежности:

(2.29)
P(t) –– вероятность нахождения системы в исправном состоянии в течение времени T ³ t.

Коэффициент стоимости эксплуатации K сэ –– отношение стоимости одного года эксплуатации системы С э к стоимости изготовления системы С 0:

(2.30)
В корректной постановке С э =С э (t) и чем больше срок эксплуатации системы, тем выше износ ее элементов и тем выше значение стоимости эксплуатации. Однако зачастую в инженерной практике принимают С э =const

Особого внимания заслуживает коэффициент эффективности системы.

(2.31)
где С 0 –– стоимость разработки (создания) системы; P(t) –– надежность технической системы; Cp(x) –– мгновенное значение прибыли; m(x) –– мгновенное значение платы за аренду (за загрязнение окружающей среды); Cэ(х) –– мгновенное значение расходов на эксплуатацию (ремонт) системы.

Рассмотрим расчетный пример

Пример 2.1.

Для расчетного примера приняты следующие значения величин:

Интенсивность отказов l=0.05 (1/год); стоимость С 0 =150 (усл. ед.); польза Ср(х)=40 (усл.ед.); амортизация m(x)=3 (усл. ед.); эксплуатация , где С Э =1(усл. ед.), К Э =0.5 (усл. ед.), что после интегрирования приводит к выражению:

Результаты расчетов по программе в Matlab:

%Вычисление рационального срока эксплуатации системы

c0=150;L=0.05;cp=40.*t;m=3.*t;ce=t+0.25.*t.^2;

B=c0+m+(1-p).*ce;

plot(t,Kf,"k-",t,K1,"k+","LineWidth",3)

Рис. 2.4. Изменение коэффициента эффективности системы и ее надежности.

Как видно из графиков рис. 2.4 срок рентабельности технической системы находится в пределах от 6 лет до 20 лет, т.е. определяется тем временем, в пределах которого числитель (2.31) превышает знаменатель или «прибыль» системы выше расходов на ее создание и эксплуатацию.

2.3. Марковские процессы, потоки событий. Элементы теории массового обслуживания.

Говорят, что в физической системе происходит случайный процесс, если она с течением времени может под влиянием случайных факторов переходить из состояния в состояние.

Потоком событий называется последовательность событий, наступающих одно за другим в случайные моменты времени.

Плотностью (интенсивностью) потока называется среднее число событий в единицу времени.

Большое значение имеют процессы, для которых состояние системы изменяется в случайные моменты времени. Особую роль играют такого рода процессы, для которых выполнены три условия:

стационарность,

отсутствие последействия,

ординарность.

Процессы, удовлетворяющие всем этим условиям, называются простейшими или однородными процессами Пуассона.

В перечисленные условия вкладывается следующий смысл.

Стационарность означает, что для любой группы из конечного числа непересекающихся промежутков времени вероятность наступления определенного числа событий на протяжении каждого из них зависит от этих чисел и от длительности промежутков времени, но не зависит от сдвига всех временных отрезков на одну и ту же величину. В частности, вероятность появления m событий в течение промежутка от t до t +Dt не зависит от t и является функцией только аргументов m, Dt.

Отсутствие последействия означает, что вероятность наступления m событий в течение интервала времени (t, t+Dt) не зависит от того, сколько раз и как появились события ранее. Это предположение означает, что условная вероятность появления m событий на промежутке (t, t+Dt) при любом предположении о наступлении событий до момента t совпадает с безусловной вероятностью. Отсутствие последействия означает взаимную независимость появления того или иного числа событий в непересекающиеся моменты времени.

Ординарность выражает собой требование практической невозможности появления двух и более событий за малый промежуток времени Dt. Точнее, это означает следующее: обозначим через вероятность появления более, чем одного события за этот малый промежуток времени. Тогда условие ординарности состоит в следующем:

Если P k (t) –– вероятность появления ровно k событий за время t:

то P 0 (t) –– можно интерпретировать как вероятность того, что длительность промежутка времени между двумя последовательными появлениями событий окажется большей t.

Если события образуют пуассоновский поток, то число m событий, попадающих на любой интервал времени (t 0 , t 0 +t) распределено по закону Пуассона:

(2.32)
где a –– математическое ожидание числа событий, попадающих на этот участок:

(2.33)
l(t) –– плотность (интенсивность) потока.

Если l(t)=const , пуассоновский поток называется стационарным пуассоновским или простейшим потоком.

Расстояние (временной интервал) Т между двумя соседними событиями в простейшем потоке есть непрерывная величина, распределенная по показательному закону с плотностью:

(2.34)
Для случайной величины Т, распределенной по показательному закону, справедливы характеристики:

(2.35)
В физической системе S происходит случайный процесс, если она с течением времени может под влиянием случайных факторов изменять свое состояние.

Случайный процесс, протекающий в системе, называется процессом с дискретным временем, если переходы системы из состояния в состояние возможны только в определенные моменты времени t 1 , t 2 ,…t n ,… Если переходы возможны в любые произвольные моменты времени, процесс называется процессом с непрерывном временем.

Случайный процесс с дискретным состоянием называется марковским , если все вероятностные характеристики в будущем зависят только от того, в каком состоянии находится этот процесс в настоящее время и не зависят от того, каким образом этот процесс протекал в прошлом. Будущее зависит от прошлого только через настоящее. Если процесс марковский, то все потоки событий, переводящие систему из состояния в состояние, являются пуассоновскими.

При анализе случайных процессов с дискретным состоянием удобно пользоваться геометрической схемой, называемой графом состояний , который изображает возможные состояния системы и возможные переходы из состояния в состояние.

Каждое состояние системы обозначается квадратом или кружком, а возможные переходы системы из состояния в состояние –– стрелками, соединяющими квадраты или кружки. Заметим (рис. 2.5), что стрелками отмечаются только непосредственные переходы системы из состояния в состояние.

Например, если система из состояния S 0 может перейти в состояние S 3 только через состояние S 1 или S 2 , то стрелками отмечаются только переходы из S 0 в S 1 и из S 0 в S 2 , а не из состояния S 0 в S 3 .

S 0

S 2

S 3

S 1

Рис. 2.5. Граф состояний системы.

В теории надежности чаще встречаются ситуации, когда переходы системы из состояния в состояние происходят в случайные моменты времени, которые заранее предсказать невозможно. Для описания таких процессов в ряде случаев может быть применена схема марковского процесса с дискретным состоянием и непрерывным временем.

Система S называется системой с дискретным состоянием, если она имеет счетное множество возможных состояний (число состояний можно перенумеровать) S 1 , S 2 ,…,S n ,… и переход из одного состояния в другое осуществляется скачком. В дальнейшем рассматриваются только системы с дискретным состоянием.

Состояние системы называется «состоянием без выхода», если из него невозможен переход ни в какое другое состояние.

Для описания случайного процесса, протекающего в системе, зачастую пользуются вероятностями состояний:

p 1 (t), p 2 (t),…,p n (t),

где p k (t) –– вероятность того, что в момент t система находится в состоянии S к.

Вероятности p k (t) удовлетворяют условию:

Введем в рассмотрение плотность l ij вероятностей перехода системы из состояния S i в состояние S j .

t

t+∆t

∆t

T

Рис. 2.6. Представление режима работы системы во времени.

Пусть система (рис.2.6) в момент t находится в состоянии S i . Рассмотрим элементарный участок Dt, примыкающий к моменту t.

Назовем плотностью вероятностей (или интенсивностью) перехода из состояния S i в состояние S j величину l ij как предел отношения вероятности перехода от состояния S i в состояние S j за время Dt к продолжительности этого промежутка времени Dt:

(2.36)
где P ij (Dt) –– вероятность того, что система, находившаяся в момент t в состоянии S i , за время Dt перейдет в состояние S j (справедливо только для i¹j).

При малом значении временного интервала Dt вероятность P ij (Dt) с точностью до бесконечно малых высшего порядка малости равна:

(2.37)
Если все интенсивности перехода l ij не зависят от времени, марковский процесс называют однородным , в противном случае –– процесс называется неоднородным.

Пусть нам известны все l ij для всех пар (S i , S j) . Построим граф состояний системы и против каждой стрелки поставим соответствующую плотность вероятности перехода (рис. 2.7.3). Такой граф называется размеченным графом состояний.

S i

S j

S k

l ij

l jk

Рис.2.7 Пример построения размеченного графа.

При наличии размеченного графа состояний системы, можно определить вероятности состояний P 0 (t), P 1 (t), P 2 (t)…как функции времени, а именно, эти вероятности удовлетворяют дифференциальным уравнениям Колмогорова.

Продемонстрируем методику вывода системы дифференциальных уравнений Колмогорова на конкретном примере (рис.2.8).

Пусть система имеет пять состояний S 0 , S 1 , S 2 , S 3 , S 4 . Поставим задачу найти одну из вероятностей состояния, например, P 0 (t). Это есть вероятность того, что в момент t система будет находиться в состоянии S 0 .

Придадим моменту t малое приращение Dt и найдем вероятность того, что в момент t+Dt система будет находиться в состоянии S 0 .

Реализация такого события возможна двумя путями:

а) система не изменит своего состояния за промежуток времени Dt;

б) система, находясь в момент t в состоянии S 3 , перейдет за Dt в состояние S 0 .

Вариант а) реализуется, если в момент t система с вероятностью P 0 (t) находилась в состоянии S 0 и не перешла из состояния S 0 в состояние S 1 . Вероятность последнего события может быть вычислена (для малых значений Dt) по формуле:

,

где P 0 (t) –– вероятность нахождения системы в момент t в состоянии S 0 , l 01 ×Dt –– вероятность перехода системы за промежуток времени Dt из состояния S 0 в состояние S 1 , (1-l 01 ×Dt) –– вероятность неперехода системы за интервал времени Dt из состояния S 0 в состояние S 1 .

l 01

S 0

l 30

l 42

S 2

S 3

l 13

S 4

l 12

l 34

S 1

Рис. 2.8. Фрагмент размеченного графа технической системы.

Вариант б) реализуется в том случае, если система в момент t находилась с вероятностью Р 3 (t) в состоянии S 3 и за интервал времени Dt перешла в состояние S 0:

,

где l 30 ×Dt –– вероятность перехода за малый интервал времени Dt системы из состояния S 3 в состояние S 0 .

Поскольку система в момент t+Dt могла находиться в состоянии Р 0 только или первым или вторым способом, то получаем:

(2.38)
или:

(2.39)
Рассмотрим состояние S 1 и выведем уравнение для определения вероятности P 1 (t) того, что в момент t+Dt система будет находиться в состоянии S 1 .

Реализация такого состояния возможна, если:

Система находилась в момент t в состоянии S 0 и за время Dt перешла в состояние S 1 . Вероятность такого перехода определяется произведением соответствующих вероятностей:

Система в момент t находилась в состоянии S 1 и за интервал Dt своего состояния не изменила, т.е. не перешла ни в состояние S 2 , ни в состояние S 3 . Оценим вероятность осуществления этого варианта.

Вероятность того, что система, находясь в состоянии S 1 , перейдет за время Dt в состояние S 2 или S 3:

Вероятность неперехода системы из состояния S 1 ни в одно из этих состояний:

Окончательно получим:

Или, при стремлении Dt к нулю, имеем окончательно:

(2.40)
Аналогичным образом могут быть получены зависимости системы дифференциальных уравнений Колмогорова для всех остальных состояний рассматриваемой системы.

В итоге получим систему дифференциальных уравнений:

(2.41)

Интегрирование этой системы дифференциальных уравнений при начальных, условиях, например,

дает искомые функции вероятностей состояний:

Все уравнения (2.41) построены по определенному правилу, зная которое можно выписывать систему для размеченного графа почти автоматически:

· в левой части каждого уравнения стоит производная ,

· в правой части содержится столько членов, сколько стрелок связано непосредственно с данным k –– м состоянием,

· член правой части уравнения имеет знак плюс , если стрелка ведет в данное состояние и знак минус , если стрелка выходит из данного состояния,

· каждый член правой части уравнения равен плотности потока событий, переводящего системупо данной стрелке, умноженной на вероятность того состояния, из которого исходит стрелка.

Эти правила составления системы дифференциальных уравнений Колмогорова справедливы для любой непрерывной марковской цепи.

Например.

Рис. 2.9. Размеченный граф системы с дискретным состоянием

и непрерывным временем.

Система дифференциальных уравнений такой системы имеет вид:

(2.42)

Начальные условия для интегрирования такой системы отражают состояние системы в начальный момент времени. Так, если в момент t=0 система была в состоянии S k , то полагают:

Число уравнений в системе может быть уменьшено на единицу, если учесть условие, что для любого t (для рассматриваемой системы).