Горелик А.В., Ермакова О.П.
ПРАКТИКУМ ПО ОСНОВАМ ТЕОРИИ НАДЕЖНОСТИ
Учебное пособие для студентов специальности 190901 «Системы
обеспечения движения поездов»
Москва - 2010
Рецензенты: Митрохин В.Е. , доктор технических наук, профессор
Зав. кафедры «Системы передачи информации» ОмГУПС.
Шалягин Д.В., доктор технических наук, профессор,
начальник отдела ПКТБ ЦШ филиала ОАО «РЖД».
Горелик А.В., Ермакова О.П. Практикум по основам теории
надежности.- М.:
В учебном пособии содержатся задачи по всем основным разделам
дисциплины «Основы теории надежности», которая входит в перечень
базовых дисциплин федерального государственного образовательного
стандарта по специальности 190901 «Системы обеспечения движения
поездов».
В него включены задачи по расчету показателей надежности
невосстанавливаемых и восстанавливаемых изделий при основном и
резервном соединении элементов, а также задачи по оценке показателей
надежности по данным об их отказах, полученным по результатам
испытаний. В каждой главе приведены краткие сведения из теории
надежности, типовые примеры и задачи с ответами. Пособие включает как
простые задачи, полезные при первоначальном изучении теории надежности,
так и более сложные, решение которых будет способствовать выработке
практических навыков.
Книга предназначена для студентов специальности 190901 «Системы
обеспечения движения поездов»
ВВЕДЕНИЕ
Проблема надежности продолжает оставаться одной из основных в
современной технике. Это объясняется тем, что постоянно
совершенствуются и усложняются технические устройства и системы,
применяемые в практической сфере, и одновременно повышаются
требования к надежности их функционирования. Инженеры, физики и
математики приложили немало совместных усилий для разработки
современной теории надежности.
Наилучшим методом, изучения основ теории надежности и
закрепления теоретических знаний является решение практических задач.
Между тем в отечественной литературе задачники по теории надежности уже
давно не издавались.
Авторы решили восполнить этот пробел и подобрали, составили и
систематизировали задачи по основным разделам дисциплины «Основы
теории надежности».
Цель учебного пособия «Практикум по основы теории надежности» -
помочь читателю изучить теорию надежности и приобрести навыки
применения ее результатов к решению различных прикладных вопросов.
Учебное пособие содержит четыре главы. В каждой главе по три
параграфа. В первом из них даются сведения из теории и расчетные
формулы, во втором – типовые примеры с решениями, в третьем – задачи с
ответами. Числовые ответы являются приближенными, хотя, как правило,
получены при выполнении расчетов с использованием программы Mathcad.
В книге помещены задачи различной трудности – от простых до
достаточно сложных. Причем наряду с оригинальными задачами,
предложенными авторами, в сборник включены типовые примеры и задачи,
взятые из источников, указанных в списке литературы, промежуточные
результаты вычислений и ответы к типовым задачам и примерам уточнены с
помощью современных математических программ. Простые задачи могут
решаться студентами при первоначальном изучении теории надежности. К
простым задачам относятся задачи первой и второй глав. Более сложные
задачи приведены в третьей и четвертой главах.
Учебное пособие, конечно, не охватывает всех разделов теории
надежности. Так в нем не содержатся задачи оценки надежности систем по
результатам испытаний ее компонентов и контроля. Однако, по мнению
авторов, решение предложенных задач будет способствовать углублению
знаний теории и выработки практических навыков.
Книга предназначена для студентов специальностей 190402
«Автоматика телемеханика и связь на железнодорожном транспорте» и
190401 «Электроснабжение железных дорог».
Показатель надежности – характеристика одного или нескольких
свойств, составляющих надежность объекта.
К числу наиболее широко применяемых показателей надежности,
которые регламентирует ГОСТ 27.002-89 (приложение 1), относятся:
Характеристикой надежности будем называть количественное
значение показателя надежности конкретного изделия.
Выбор количественных характеристик надежности зависит от вида
изделия.
Основные показатели надежности можно разбить на две группы:
Невосстанавливаемыми называются такие изделия, для которых в
рассматриваемой ситуации проведение восстановления работоспособного
состояния не предусмотрено в нормативно-технической и (или)
конструкторской документации. Если происходит отказ такого изделия, то
выполняемая операция будет сорвана, и ее необходимо начинать вновь в том
случае, если отказ можно устранить. К таким изделиям относятся изделия
однократного действия, такие как ракеты, управляемые снаряды,
искусственные спутники Земли, а также системы многократного действия,
такие как системы управления воздушным и железнодорожным движением,
системы управления химическими, металлургическими и другими
ответственными производственными процессами.
Восстанавливаемыми называются такие изделия, для которых в
рассматриваемой ситуации проведение восстановления работоспособного
состояния предусмотрено в нормативно-технической и (или)
конструкторской документации. Если произойдет отказ такого изделия, то он
вызовет прекращение функционирования изделия только на период
устранения отказа. К таким изделиям относятся: телевизор, агрегат питания,
локомотив, автомобиль и т. п.
На рис. 1.1 показан временной график работы невосстанавливаемых и
восстанавливаемых изделий.
Н.о К.о Н.о К.о
tHP tH. p
t П1 tП 2
Рисунок 1.1. Временной график работы невосстанавливаемых и
восстанавливаемых изделий:
Рассмотрим следующую модель испытаний. На испытании находится
No изделий, и испытания считаются законченными, если все они отказали.
Причем отказавшие изделия отремонтированными или новыми не
заменяются. Тогда показателями надежности данных изделий являются:
- вероятность безотказной работы P(t);
Вероятность безотказной работы выражает вероятность того, что в
пределах заданной наработки отказ изделия не возникнет.
Вероятность безотказной работы по статистическим данным об отказах
оценивается согласно выражению
P( t)=NO ■- n(t), (1.1)
NO
где n(t) - количество изделий, отказавших к моменту времени t, при их
исходном количестве No ; P(t) - статистическая оценка вероятности
безотказной работы. При большом числе изделий NO статистическая оценка
P(t) практически совпадает с вероятностью безотказной работы P(t). На
практике иногда более удобной характеристикой является вероятность отказа
Q (t).
Вероятность отказа - вероятность того, что при определенных
условиях эксплуатации в заданном интервале времени произойдет хотя бы
один отказ. Отказ и безотказная работа являются событиями
несовместимыми и противоположными, поэтому
Q( t) = nt), Q (t) = 1 - P (t).
NO
(1.2)
Частотой отказов называется отношение числа отказавших изделий в
единицу времени к первоначальному числу изделий при условии, что все
вышедшие из строя изделия не восстанавливаются.
Частоту отказов по статистическим данным об отказах оценивается
согласно выражению
a( t) =
n (At)
NO A t ’
(1.3)
где n (At) - число отказавших изделий в интервале времени
от t -
A t
до
A t
t +--
2
Частота отказов есть плотность вероятности (или закон
распределения) времени работы изделия до первого отказа. Поэтому
dP (t)_ dQ (t)
a( t) = = ,
dt dt
t t ^
0 0 t
Интенсивностью отказов называется отношение числа отказавших
изделий в единицу времени к среднему числу изделий, исправно работающих
в данный отрезок времени.
Согласно определению интенсивность отказов по статистическим
данным об отказах определяется
N + Nm
где Nср = i 2 i+1
Д t) =
n (At)
IT 4 ,
Nср At
(1.5)
- среднее число исправно работающих изделий в
интервале At; Ni- число изделий, исправно работающих в начале интервала
At; Ni+1 - число изделий, исправно работающих в конце интервала At.
Интенсивность отказов есть условная плотность вероятности
возникновения отказа невосстанавливаемого изделия, определяемая для
рассматриваемого момента времени при условии, что до этого момента
отказ не возник.
Вероятностная оценка этой характеристики находится из выражения
X (t) = ^(tl. (1.6)
P (t)
Интенсивность отказов и вероятность безотказной работы связаны
между собой зависимостью
t
-f Л( t) dt
P(t) = e 0 (1.7)
Средней наработкой до отказа ТO называется математическое
ожидание наработки изделия до первого отказа.
Как математическое ожидание ТO вычисляется через частоту отказов
(плотность распределения времени безотказной работы):
TO = М [t ]= J ta(t)dt. (1.8)
-да
Так как t положительно и P(0) = 1, а P(да) = 0, то
TO = jP(t)dt. (1.9)
0
По статистическим данным об отказах средняя наработка до первого
отказа вычисляется по формуле
NO
X ti
~
To = (1.10)
NO
где ti – время безотказной работы i- го изделия, NO - число испытуемых
изделий.
Как видно из формулы (1.10), для определения средней наработки до
первого отказа необходимо знать моменты выхода из строя всех испытуемых
~
изделий. Поэтому для вычисления TO пользоваться указанной формулой
неудобно. Имея данные о количестве вышедших из строя изделий ni в
каждом i-ом интервале времени, среднюю наработку до первого отказа
лучше определить из выражения
m
∑ nit срi
i=i
O NO
~
.
(1.11)
В выражении (1.11) tсрi и m находятся по следующим формулам
t =
срi
ti-i + ti
2
t
, m= n,
∆t
где ti-1 - время начала i-ого интервала; ti - время конца i-ого интервала; tn -
время, в течение которого вышли из строя все изделия; Лt = ti-1 - ti -
интервал времени.
При изучении надежности технических устройств часто применяются
следующие законы распределения времени безотказной работы:
экспоненциальный, нормальный, Релея, Гамма, Вейбулла.
В табл. 1.1 приведены выражения для оценки количественных
характеристик надежности изделий при указанных законах распределения
времени безотказной работы.
Таблица 1.1. Основные соотношения количественных характеристик
надежности при различных законах распределения времени безотказной
работы
Закон |
Частота |
Вероятность |
Интенсивность |
Экспоненциальный |
λe-λt |
e-λt |
λ = const |
Релея |
t2 — е 2°2 |
_ t2 e 2°2 |
t О2 |
Вейбулла |
20 ktk-1 e ■ Л° tk |
e Л'tk |
Л0 ktk-1, |
Нормальный |
(t - T )2 1 2^ |
t -(t-T)2 1 - ∫ e 2σ dt 2πσ 0 |
t (t-T)2 1 - 1 ∫ e 2σ2 2πσ 0 |
J 2na | |||
(t-T)2 2πσ | |||
Гамма |
Л e - *0 ■ (k — 1)! |
e ~ Л° t^*0 t) |
*0(*0t)k-1 (k -1)! kf*< |
Из выражений для оценки количественных характеристик надежности
видно, что все характеристики являются функциями времени. На рис. 1.2
приведены зависимости количественных характеристик надежности изделий
г
д
Рисунок 1.2. Зависимость количественных характеристик от времени
Рассмотренные показатели надежности позволяют достаточно полно
оценить надежность невосстанавливаемых изделий. Они также позволяют
оценить надежность восстанавливаемых изделий до первого отказа.
Наличие нескольких показателей вовсе не означает, что всегда нужно
оценивать надежность изделия по всем показателям.
Наиболее полно надежность изделия характеризуется частотой отказов
a(t). Это объясняется тем, что частота отказов является плотностью
распределения, а поэтому несет в себе всю информацию о случайной
величине – времени безотказной работы.
Средняя наработка до первого отказа является достаточно наглядной
характеристикой надежности. Однако применение этого показателя для
оценки надежности сложной системы ограничено в тех случаях, когда:
- время работы системы гораздо меньше среднего времени безотказной
работы;
однопараметрический;
Интенсивность отказов - наиболее удобная характеристика надежности
простейших элементов, так как она позволяет более просто вычислить
количественные характеристики сложной системы.
Наиболее целесообразным показателем надежности сложной системы
является вероятность безотказной работы. Это объясняется следующими
особенностями вероятности безотказной работы:
В табл.1.2 показана взаимосвязь между показателями надежности.
Таблица 1.2 Взаимосвязь между показателями надежности
Известно | Требуется определить | |||
P (t) | Q (t) | f (t) | ^( t) | |
P (t) | - | 1 - p (t) | dP (t) | _ 1 dP (t) P (t) dt |
Q (t) | 1 - Q (t) | - | dQ (t) | 1 dQ (t) 1 - Q (t) dt |
f (t) | ∞ J f (t )dt t | t J f (t) dt 0 | - | f (t) |
Л( t) | t -J A( t) dt | t -J Л( t) dt | t -J Л( t) dt | - |
ПОКАЗАТЕЛИ НАДЕЖНОСИ ВОССТАНАВЛИВАЕМЫХ ИЗДЕЛИЙ
Рассмотрим следующую модель испытаний. Пусть на испытание
поставлено N изделий. Отказавшие изделия немедленно заменяются
исправными (новыми или отремонтированными). Испытания считаются
законченными, если число отказов достигает величины достаточной для
оценки надежности с определенной доверительной вероятностью. Если не
учитывать время необходимое на восстановление системы, то
количественными характеристиками надежности могут быть параметр
потока отказов ^( t) и средняя наработка на отказ T.
Параметр потока отказов есть отношение среднего числа отказов
восстанавливаемого объекта за произвольно малую его наработку к
значению этой наработки.
Согласно определению
M |n(t + Лt)|-M|n(t)|
^(t) = lim ——-----—---- v ' = M' | n (t) |,
где Mn(t + Лt)| и Mn(t)| - математическое ожидание числа отказов за время
t + Л t, t.
По статистическим данным параметром потока отказов называется
отношение числа отказавших изделий в единицу времени к числу
испытываемых изделий при условии, что все вышедшие из строя изделия
заменяются исправными (новыми или отремонтированными).
Согласно определению
_. . n(Лt)
~(t) = —,
N⋅∆t
где n(Лt) - число отказавших изделий в интервале времени от
(1.12)
∆t
t-— до
2
t + ~; N — число испытываемых изделий; Лt - интервал времени.
Параметр потока отказов обладает следующими важными свойствами:
1
t ^ ® стремится к —. Это означает, что при длительной эксплуатации
TO
ремонтируемого изделия поток его отказов независимо от закона
распределения времени безотказной работы становится стационарным.
Однако это вовсе не означает, что интенсивность отказов есть величина
постоянная;
N
параметров потоков отказов элементов, т.е. ^с(t) ^ ^^i(t).
i=1
Средняя наработка на отказ есть отношение наработки
восстанавливаемого изделия к математическому ожиданию числа его
отказов в течение этой наработки.
Эта характеристика определяется по статистическим данным об
отказах по формуле
n
_ Е ti
T = i=L, (1.13)
n
где ti - время исправной работы изделия между (i - 1)-м и i-м отказами; n -
число отказов за некоторое время t.
Из формулы (1.13) видно, что в данном случае наработка на отказ
определяется по данным испытания одного образца изделия. Если на
испытании находится N образцов в течение времени t, то наработка на отказ
вычисляется по формуле
N nj
T = У , (1.14)
j=1
где tij - время исправной работы j-го образца изделия между (i - 1)-м и i-м
отказом; nj- - число отказов за время tj-го образца.
Наработка на отказ является достаточно наглядной характеристикой
надежности, поэтому она получила широкое распространение на практике.
Параметр потока отказов и наработка на отказ характеризуют
надежность ремонтируемого изделия, но не учитывают времени,
необходимого на его восстановление. Поэтому они не характеризуют
готовности изделия к выполнению своих функций в нужное время. Для этих
целей вводятся такие показатели, как коэффициент готовности и
коэффициент вынужденного простоя.
Коэффициент готовности есть вероятность того, что изделие
окажется в работоспособном состоянии в произвольный момент времени,
кроме планируемых периодов, в течение которых применение объекта по
назначению не предусматривается.
Эта характеристика определяется по статистическим данным как
отношение времени исправной работы к сумме времен исправной работы и
вынужденных простоев изделия, взятых за один и тот же календарный срок.
Эта характеристика обозначается KГ .
Согласно данному определению
К Г = -P-, (1.15)
t P + t П
где tP – суммарное время исправной работы объекта; tП – суммарное время
вынужденного простоя.
Времена tP и tП вычисляются по формулам
tP=itpi, t п=Ь ni, (1.16)
i=1 i=1
где tPi - время работы изделия между (i – 1)-м и i-м отказом; tПi - время
вынужденного простоя после i-го отказа; n – число отказов (ремонтов)
изделия.
Для перехода к вероятностному показателю величины tP и tП
заменяются математическими ожиданиями времени между соседними
отказами и времени восстановления соответственно.
Тогда
К
T
т+Гв
(1.17)
где T - среднее время наработки на отказ; TВ - среднее время
восстановления.
Среднее время восстановления есть математическое ожидание
времени восстановления работоспособного состояния объекта.
Как математическое ожидание ТB вычисляется через частоту
восстановления (плотность распределения времени восстановления):
TB = J ta B (t) dt,
0
(1.18)
где a - частота восстановления, равная a (t) = ^(t), (S(t) - вероятность
B в dt
восстановления).
По статистическим данным среднее время восстановления вычисляется
по формуле:
ni
~
T — I = 1
(1.19)
ni
Коэффициентом вынужденного простоя называется отношение
времени вынужденного простоя к сумме времен исправной работы и
вынужденных простоев изделия, взятых за один и тот же календарный срок.
Согласно определению
~ _ t П
(1.20)
или переходя к средним величинам,
Kn = (1.21)
П T + T В
Коэффициент готовности и коэффициент вынужденного простоя
связаны между собой зависимостью
Kn = 1 - K г •
(1.22)
При анализе надежности восстанавливаемых систем обычно
коэффициент готовности вычисляют по формуле
К
TO
To + Тв
(1.23)
Формула (1.23) верна только в том случае, если поток отказов
простейший, и тогда T = To •
Часто коэффициент готовности, вычисленный по формуле (1.23)
отождествляют с вероятностью того, что в любой момент времени
восстанавливаемая система исправна. На самом деле указанные
характеристики неравноценны и могут быть отождествлены при
определенных условиях.
Поэтому вводится понятие функции готовности, которая определяется
из выражения:
Pr (t ) = K Г (t) = ^— + -A- e ~(A+»t
A + p A + p
t
Pr (t ) = K Г (t) = K Г + (1 - K Г) e K Г ■ TB
(1.24)
; 1 1
где A = ; p =
TO TВ
Кроме установившегося коэффициента готовности часто используется
среднее значение этого коэффициента за интервал времени 0...ti
ti
K Г = - f K Г (t) dt
ti 0
(1.25)
Коэффициент оперативной готовности KоГ - вероятность того, что
объект окажется в работоспособном состоянии в произвольный момент
времени, кроме планируемых периодов, в течение которых применение
объекта по назначению не предусматривается, и, начиная с этого момента,
будет работать безотказно в течение заданного интервала времени
Согласно определению
Kог = Kг • e-At (1.26)
Коэффициент оперативной готовности оценивает не только готовность
системы к выполнению заданных функций, но и способность выполнять
функции определенной временной протяженностью.
Задачи, которые встречаются при определении количественных
характеристик надежности, могут быть разделены на следующие группы:
В настоящей главе при определении количественных характеристик
надежности технических устройств по статистическим данным об их отказах
не учитывается достоверность полученных результатов.
Кроме того, следует иметь в виду, что частота, интенсивность и
параметр потока отказов, вычисленные по формулам (1.3), (1.5), (1.12)
являются постоянными в диапазоне интервала Лt, а функции а(t), Л(t),
S)( t) - ступенчатыми кривыми или гистограммами. Для удобства изложения
в дальнейшем при решении задач на определение частоты, интенсивности и
параметра потока отказов по статистическим данным об отказах изделий
ответы относятся к середине интервала Лt. При этом результаты вычислений
графически представляются не в виде гистограмм, а в виде точек,
отнесенных к середине интервалов Лt и соединенных плавной кривой.
Рассмотрим типовые примеры.
Пример 1.1. На испытание поставлено 1000 однотипных электронных
ламп. За 3000 ч отказало 80 ламп. Требуется определить вероятность
безотказной работы и вероятность отказа электронных ламп в течение 3000 ч.
Решение. По формулам (1.1) и (1.2) определяем
P(3000) =
NO - n(t) _ 1000 - 80
NO ~ 1000
2(3000) =
n^ = -80- = 0,08
NO 1000
или
2(3000) = 1 - P( t) = 1 - 0,92 = 0,08
Пример 1.2. На испытание поставлено 1000 однотипных ламп. За
первые 3000 ч работы отказало 80 ламп, а за интервал 3000 ч – 4000 ч
отказало еще 50 ламп. Определить частоту и интенсивность отказов
электронных ламп в промежутке 3000 ч – 4000 ч работы.
Решение. По формулам (1.3) определим частоту отказов
a(3500) =
n(Лt) _ 50
NOЛt " 1000 -1000
=5-10-5 (1/ч)
Определим среднее число исправно работающих изделий в интервале
Л t.
N
ср
Ni + Ni+1 _ 920 + 870
2 = 2
= 895 (шт)
По формуле (1.5) находим интенсивность отказов
λ~(3500) =
n(∆t)
Nср∆t
50
895 ⋅1000
= 5,6 ⋅10-5(1/ч)
Пример 1.3. На испытании находилось 1000 образцов
неремонтируемой аппаратуры. Число отказов n(∆t) фиксировалось через
каждые 100 ч работы (∆t = 100 ч). Данные об отказах приведены в табл. 1.3.
Требуется вычислить количественные характеристики и построить
зависимость характеристик от времени.
Таблица 1.3. Данные об отказах к примеру 1.3
∆ti, ч |
n(∆ti) |
∆ti, ч |
n(∆ti) |
∆ti, ч |
n(∆ti) |
0 – 100 |
50 |
1000 - 1100 |
15 |
2000 – 2100 |
12 |
100 – 200 |
40 |
1100 - 1200 |
14 |
2100 – 2200 |
13 |
200 – 300 |
32 |
1200 – 1300 |
14 |
2200 – 2300 |
12 |
300 – 400 |
25 |
1300 – 1400 |
13 |
2300 – 2400 |
13 |
400 – 500 |
20 |
1400 – 1500 |
14 |
2400 – 2500 |
14 |
500 -600 |
17 |
1500 – 1600 |
13 |
2500 – 2600 |
16 |
600 – 700 |
16 |
1600 – 1700 |
13 |
2600 – 2700 |
20 |
700 – 800 |
16 |
1700 – 1800 |
13 |
2700 – 2800 |
25 |
800 – 900 |
15 |
1800 – 1900 |
14 |
2800 – 2900 |
30 |
900 - 1000 |
14 |
1900 - 2000 |
12 |
2900 - 3000 |
40 |
Решение. Аппаратура относится к классу невосстанавливаемых
~ ~~
изделий. Поэтому показателями надежности будут P (t) , α~(t) , λ(t) , TO .
Вычислим P(t) . На основании формулы (1.1) имеем
P~(100) =
NO - n(100)
1000 - 50
1000
P~(200) =
NO - n(200)
NO
1000 - 90
1000
P~(3000) = NO
- n(3000) 1000 - 575
NO
1000
= 0,425.
Для расчета характеристик α(t) и λ(t) применим формулы (1.3) и (1.5),
тогда
α~(50) =
n(∆t)
NO∆t
50
1000⋅100
=5⋅10-4
(1/ч),
α~(150) = n(∆t) = 40
NO∆t 1000 ⋅100
= 4 ⋅ 10-4
(1/ч),
£(2950) = n^ =---40---= 4 -10-4 (1/ч),
NO A t 1000 -100
1(50) = n(At) =---
Nср At 100 -
50-------v = 5,13 -10-4 (1 /ч),
^ 1000 + 950 ^ , v ’
I 2 )
/1(150) =
n (At) _ 40
NсрAt = 100-f 950+910 1
I 2 )
= 4,3 -10-4 (1/ч),
/(2950) = n^)- =---
Nср At100 -
40
—°------ = 9 -10-4 (1/ч)
f 465+4251
I 2)
~ ^Z
Значения P(t), a(t), /(t), вычисленные для всех Ati, приведены в
табл. 1.4, а их зависимости от времени на рис. 1.3 и 1.4.
Вычислим среднее время безотказной работы, предположив, что на
испытании находились только те образцы, которые отказали.
Расчет проведем, используя формулу (1.11), учитывая, что
t
m = —
A t
3000
100
= 30 и NO = 575, имеем
m
E nitcpi
т i=1
TO =
~
NO
50 - 50 + 40-150 + 32 - 250 +... + 30 - 2850 + 40 - 2950
575
= 1400 (ч).
Полученное значение средней наработки до первого отказа является
заниженным, т.к. опыт был прекращен после отказа 575 образцов из 1000,
поставленных на испытание.
~ ~
Таблица 1.4. Вычисленные значения P (t), a( t), /(t) к примеру 1.3
A ti, ч | P( t) | £(t) -10-4, (1/ч) | /(t) -10-4 , (1/ч) |
0 – 100 | 0,950 | 5 | 5,14 |
100 – 200 | 0,910 | 4 | 4,3 |
200 – 300 | 0,878 | 3,2 | 3,58 |
300 – 400 | 0,853 | 2,5 | 2,89 |
400 – 500 | 0,833 | 2 | 2,38 |
500 -600 | 0,816 | 1,7 | 2,06 |
600 – 700 | 0,800 | 1,6 | 1,98 |
700 – 800 | 0,784 | 1,6 | 2,02 |
800 – 900 | 0,769 | 1,5 | 1,93 |
900 – 1000 | 0,755 | 1,4 | 1,84 |
1000 - 1100 | 0,740 | 1,5 | 2 |
1100 - 1200 | 0,726 | 1,4 | 1,91 |
1200 – 1300 | 0,712 | 1,4 | 1,95 |
1300 – 1400 | 0,699 | !.3 | 1,84 |
1400 – 1500 | 0,685 | 1,4 | 2,02 |
1500 – 1600 | 0,672 | 1,3 | 1,92 |
1600 – 1700 | 0,659 | 1,3 | 1,95 |
1700 – 1800 | 0,646 | 1,3 | 2 |
1800 – 1900 | 0,632 | 1,4 | 2,2 |
1900 – 2000 | 0,620 | 1,2 | 1,92 |
2000 – 2100 | 0,608 | 1,2 | 1,95 |
2100 – 2200 | 0,595 | 1,3 | 2,17 |
2200 – 2300 | 0,583 | 1,2 | 2,04 |
2300 – 2400 | 0,570 | 1,3 | 2,25 |
2400 – 2500 | 0,556 | 1,4 | 2,48 |
2500 – 2600 | 0,540 | 1,6 | 2,9 |
2600 – 2700 | 0,520 | 2 | 3,76 |
2700 – 2800 | 0,495 | 2,5 | 4,9 |
2800 – 2900 | 0,465 | 3 | 6,24 |
2900 - 3000 | 0,425 | 4 | 9 |
~
Рисунок 1.3. Зависимость P от t к примеру 1.3.
Рисунок 1.4. Зависимость α~ и λ~ от t к примеру 1.3.
Пример 1.4. В течение некоторого периода времени производилось
наблюдение за работой одного экземпляра радиолокационной станции. За
весь период наблюдения было зафиксировано 15 отказов. До начала
наблюдения станция проработала 258 часов, а к концу наблюдения наработка
станции составила 1233 часа. Требуется определить среднюю наработку на
отказ T .
Решение. Наработка радиолокационной станции за наблюдаемый
период равна
t = 12 -11 = 1233 - 258 = 975 (ч)
n
Принимая ^ ti = 975 ч, по формуле (1.13) находим среднюю наработку
i=1
на отказ:
n
t
~=а:=975 = 65 (ч).
n 15
Пример 1.5. Производилось наблюдение за работой трех экземпляров
однотипной аппаратуры. За период наблюдения было зафиксировано по
первому экземпляру 6 отказов, по второму и третьему – 11 и 8 отказов
соответственно. Наработка первого экземпляра составила 181 ч, второго –
329 ч и третьего 245 ч. Требуется определить наработку аппаратуры на отказ
~
Решение. Определим суммарную наработку трех образцов аппаратуры:
N nj
t У = ЕЁ tj = 181 + 329 + 245 = 755 (ч)
j=1i=1
Определим суммарное количество отказов:
N
n Е = Е nj = 6 +11 + 8 = 25 (отказов).
j=1
Находим среднюю наработку на отказ по формуле (1.14)
~
N nj
j=1i=1
N
j=1
t^=
n Е
755
25
= 30,2 (ч).
Пример 1.6. Система состоит из 5 приборов, причем отказ любого
одного из них ведет к отказу системы. Известно, что первый прибор отказал
34 раза в течение 952 ч работы, второй -24 раза в течение 960 ч работы, а
остальные приборы в течение 210 часов работы отказали 4, 6 и 5 раз
соответственно. Требуется определить наработку до отказа системы в целом,
если справедлив экспоненциальный закон надежности для каждого из пяти
приборов.
Решение. Для решения данной задачи воспользуемся следующими
соотношениями:
n
i=1
To = —.
O Л
Определим интенсивность отказов для каждого прибора:
~ 34 ~24
Л~ =---= 0,0357 (1/ч), Л =---= 0,025 (1/ч),
1 952 2960
~4 6 5
^,4,5 = “W = 0,0714 (1/ч).
Интенсивность отказов системы равна
Л =2 Л = Л + Л + Л3,4,5 = 0,0357 + 0,025 + 0,0714 = 0,1321 (1/ч)
i=1, ,
Средняя наработка до отказа системы
To = — = 1 = 7,57 (ч)-
O Лс 0,1321
Пример 1.7. За наблюдаемый период эксплуатации в аппаратуре было
зафиксировано 8 отказов. Время восстановления составило t1 = 12 мин, t2 = 23
мин, t3 = 15 мин, t4 = 9 мин, t5 = 17 мин, t6 = 28 мин, t7 = 25 мин, t8 = 31 мин.
~
Определить среднее время восстановления аппаратуры TВ .
Решение.
ni
8
~
Т — 1 =1
TB =
ni
- i=1
600
88
= 75 (мин).
Пример 1.8. Аппаратура имела среднюю наработку на отказ 65 ч и
среднее время восстановления 1,25 ч. Требуется определить коэффициент
готовности.
Решение. По формуле (1.17) имеем
К
T
65
Г Т + TB 65 +1,25
= 0,98
Пример 1.9. Пусть время работы элемента до отказа подчинено
экспоненциальному закону распределения с параметром λ= 2,5 ⋅ 10-5 1/ч.
Требуется вычислить количественные характеристики надежности
элемента P(t), α(t), TO , если t = 500, 1000, 2000 ч.
Решение. Используем формулы для P(t), α(t) и TO , приведенные в
табл.1.1.
Вычислим вероятность безотказной работы:
P(t) = e-λ⋅t = e-2,5⋅10-5⋅t .
P(500) = e-2,5⋅10-5⋅500 =0,9875;
P(1000)=e-2,5⋅10-5⋅1000 =0,9753;
P(2000) = e-2,5⋅10-5⋅2000 = 0,9512.
Вычислим частоту отказа:
α(t)=λ(t)⋅P(t)=2,5⋅10-5 ⋅ e-2,5⋅10-5⋅t
α(500) =2,5⋅10-5 ⋅ e-2,5⋅10-5⋅500 =2,5⋅10-5⋅0,9875=2,469⋅10-5(1/ч);
α(500)=2,5⋅10-5⋅e-2,5⋅10-5⋅1000 =2,5⋅10-5 ⋅0,9753=2,439⋅10-5(1/ч);
α(500) =2,5⋅10-5 ⋅ e-2,5⋅10-5⋅2000 =2,5⋅10-5 ⋅0,9512=2,378⋅10-5(1/ч).
Вычислим среднюю наработку до первого отказа:
1
TO =λ
1
2,5⋅10-5
= 40000 (ч).
Пример 1.10. При эксплуатации системы было зарегистрировано n = 40
отказов. Распределение отказов по группам элементов и время, затраченное
на восстановление, приведены в табл. 1.5 Необходимо найти величину
среднего времени восстановления системы.
Решение. Определяем среднее время восстановления аппаратуры по
группам элементов.
Для полупроводниковых приборов:
ni
∑ti
i =1
TB =
~
ni
8
∑ti
= i=1
8
600
8
= 75 (мин)
Аналогично:
Рассчитаем среднее время восстановления системы по формуле:
~m
TBC = E tei ’ mi ,
i=1
где tвi - среднее время восстановления элементов i-й группы; mi - вес отказов
по группам элементов.
TBC = 0,2 • 75 + 0,25 • 76 + 0,1 -113 + 0,35 • 50 + 0,1-120 = 74,8 (мин).
Таблица 1.5. Количество зарегистрированных отказов по группам к
примеру 1.10
Группа элементов | Кол-во отказов по | Вес отказов по n группе mi = — | Время | |
ППП | 8 | 0,2 | 80 59 43 99 73 | |
Резисторы и конденсаторы | 10 | 0,25 | 61 58 54 94 | |
Реле, | 4 | 0,1 | 102 98 124 128 | |
ЭВП | 14 | 0,35 | 60 64 56 36 65 44 42 |
33 32 23 86 75 61 23 | |||
Прочие элементы | 4 | 0,1 | 125 |
133 | |||
115 | |||
107 |
1.3. Задачи
В настоящем разделе приведены задачи по расчету количественных
показателей надежности, предлагаемые для самостоятельного решения. Эти
задачи легко решить, используя типовые примеры, приведенные в
предыдущем разделе.
Ответ: P(3000) = 0,5; P(3100) = 0,25; P(3050) = 0,375; a(3050) = 2,5 -10-3
1/ч; λ~(3050) = 6,67⋅10-3 1/ч.
Ответ: PA (100) = 0,9; PB (100) = 0,8; PC (100) = 0,5; aA (50) = 10-31/ч;
α~B(50)=2⋅10-3 1/ч; α~C(50)=5⋅10-31/ч; λ~A(50)=1,05⋅10-3 1/ч;
λ~B(50)=2,2⋅10-3 1/ч; λ~C(50) =6,67⋅10-3 1/ч.
отказы через каждые 100 часов работы. Данные об отказах приведены в табл.
1.6.
Таблица 1.6. Исходные данные к задаче 1.4
∆ti, ч | n (Л ti) | ∆ti, ч | n (Л ti) |
0 - 100 | 45 | 800 - 900 | 16 |
100 - 200 | 40 | 900 - 1000 | 16 |
200 - 300 | 35 | 1000 - 1100 | 15 |
300 - 400 | 32 | 1100 - 1200 | 14 |
400 - 500 | 28 | 1200 - 1300 | 15 |
500 - 600 | 25 | 1300 - 1400 | 13 |
600 - 700 | 20 | 1400 - 1500 | 14 |
700 - 800 | 17 | 1500 - 1600 | 13 |
Необходимо определить показатели надежности: вероятность
безотказной работы, вероятность отказа, частоту отказов, интенсивность
отказов, среднее время наработки до отказа. Построить график зависимости
вероятности безотказной работы от времени
~z ~Z ~
Ответ: P (t), a( t), Л( t) см. табл. 1.7 и рис.1.5, To = 575,42 ч
Таблица 1.7. Вычисленные значения P( t), а( t), Л( t) к задаче 1.4
∆ti, ч | P( t) | α~(t) ⋅10-4 1/ч | λ~(t)⋅10-4 1/ч |
0 - 100 | 0,97 | 2,81 | 2,8 |
100 - 200 | 0,947 | 2,5 | 2,61 |
200 - 300 | 0,925 | 2,19 | 2,34 |
300 - 400 | 0,905 | 2 | 2,186 |
400 - 500 | 0,8875 | 1,75 | 1,95 |
500 - 600 | 0,872 | 1,56 | 1,776 |
600 - 700 | 0,859 | 1,25 | 1,444 |
700 - 800 | 0,849 | 1,06 | 1,244 |
800 - 900 | 0,839 | 1 | 1,185 |
900 - 1000 | 0,829 | 1 | 1,199 |
1000 - 1100 | 0,819 | 0,94 | 1,138 |
1100 - 1200 | 0,811 | 0,875 | 1,074 |
1200 - 1300 | 0,801 | 0,94 | 1,163 |
1300 - 1400 | 0,793 | 0,813 | 1,019 |
1400 - 1500 | 0,784 | 0,875 | 1,109 |
1500 - 1600 | 0,776 | 0,813 | 1,041 |
P(t)
о .в
0.4
0.2
О 500 ЮОО 1500 2000
Рисунок 1.5 Зависимость P от t (к задаче 1.4)
построить графики этих функций.
Таблица 1.8. Исходные данные к задаче 1.5
∆ti , ч | 1 группа | 2 группа | 3 группа | ∑ n(ti ) |
0 - 25 | 4 | 6 | 5 | 15 |
25 - 50 | 8 | 9 | 8 | 25 |
50 - 75 | 6 | 5 | 7 | 18 |
75 - 100 | 3 | 4 | 5 | 12 |
100 - 150 | 5 | 5 | 6 | 16 |
150 - 200 | 4 | 3 | 3 | 10 |
200 - 250 | 1 | 3 | - | 4 |
250 - 300 | 2 | 2 | - | 4 |
300 - 400 | 3 | 4 | - | 7 |
400 - 550 | 5 | - | - | 5 |
P P
Ответ: P(t), а(t), Л(t) см. табл.1.9. и рис. 1.6 и 1.7.
P P
Таблица 1.9. Вычисленные значения P (t), а( t), A( t) к задаче 1.5
∆ti, ч | P( t) | α~(t)⋅10-3 1/ч | λ~(t) ⋅10-3 1/ч |
0 - 25 | 0,95 | 2 | 2,05 |
25 - 50 | 0,867 | 3,3 | 3,67 |
50 - 75 | 0,807 | 2,4 | 2,87 |
75 - 100 | 0,767 | 1,6 | 2,03 |
100 - 150 | 0,713 | 1,1 | 1,44 |
150 - 200 | 0,68 | 0,67 | 0,957 |
200 - 250 | 0,67 | 0,4 | 0,588 |
250 - 300 | 0,65 | 0,4 | 0,606 |
300 - 400 | 0,615 | 0,35 | 0,558 |
400 - 550 | 0,59 | 0,33 | 0,542 |
~
Рисунок 1.7 Зависимость α~ и λ от t (к задаче 1.5)
~
концу наблюдения составило t2 ч. Требуется найти наработку на отказ T .
Исходные данные для решения задач и ответы приведены в табл. 1.10.
Таблица 1.10. Исходные данные и ответы к задачам 1.6 – 1.9
Номер задачи | Исходные данные | Ответ T , ч | ||
t1 , ч | t2 , ч | n(∆t) | ||
1.6 | 350 | 1280 | 15 | 62 |
1.7 | 400 | 1600 | 3 | 400 |
1.8 | 1000 | 6400 | 9 | 600 |
1.9 | 770 | 4800 | 7 | 575 |
1.10 – 1.13. В течение некоторого времени проводилось наблюдение за
работой NO экземпляров восстанавливаемых изделий. Каждый из образцов
проработал ti ч и имел ni отказов. Требуется определить наработку на отказ
T по данным наблюдения за работой всех изделий.
Исходные данные для решения задач и ответы приведены в табл. 1.11.
Таблица 1.11. Исходные данные и ответы к задачам 1.10 – 1.13
Номер | Исходные данные | Ответ | |||||||||
n1 | t1 ,ч | n2 | t2 ,ч | n3 | t3 , ч | n4 | t4 , ч | n5 | t5 , ч | ||
1.10 | 1 | 300 | 3 | 600 | 2 | 400 | - | - | - | - | 21,6 |
1.11 | 3 | 90 | 6 | 270 | 4 | 140 | 5 | 230 | 3 | 180 | 43 |
1.12 | 12 | 960 | 15 | 1112 | 8 | 808 | 7 | 1490 | - | - | 104 |
1.13 | 8 | 144 | 5 | 150 | 4 | 112 | 8 | 216 | - | - | 26 |
1.14 – 1.17. Электронная аппаратура состоит из k групп элементов. В
процессе эксплуатации зафиксировано n отказов. Количество отказов в j-й
группе равно nj ; среднее время восстановления элементов j-й группы равно
~
tj мин. Требуется вычислить среднее время восстановления аппаратуры TВ .
Исходные данные для решения задач и ответы приведены в табл. 1.12.
Таблица 1.12. Исходные данные и ответы к задачам 1.14 – 1.17
Номер | Исходные данные | Ответ | |||||||||||
k | n | n1 | t1 | n2 | t2 | n3 | t3 | n4 | t4 | n5 | t5 | ||
1.14 | 5 | 12 | 1 | 20 | 4 | 20 | 3 | 16 | 2 | 36 | 2 | 40 | 28,3 |
1.15 | 5 | 40 | 5 | 15 | 8 | 25 | 12 | 60 | 6 | 40 | 9 | 20 | 35,4 |
1.16 | 4 | 9 | 2 | 37 | 1 | 480 | 2 | 60 | 4 | 25 | - | - | 86 |
1.17 | 5 | 18 | 3 | 72 | 5 | 40 | 4 | 36 | 2 | 120 | 4 | 60 | 57,8 |
1.18 – 1.21. Изделие имеет среднюю наработку на отказ T и среднее
время восстановления TВ . Необходимо определить коэффициент готовности
изделия.
Исходные данные для решения задач и ответы приведены в табл. 1.13.
Таблица 1.13. Исходные данные и ответы к задачам 1.18 – 1.21
Номер | Исходные данные | Ответ KГ | |
T, ч | TВ, ч | ||
1.18 | 230 | 12 | 0,95 |
1.19 | 556 | 23 | 0,96 |
1.20 | 556 | 2,5 | 0,995 |
1.21 | 430 | 8 | 0,98 |
Ответ: P(6) = 0,995, а(100) = 0,75-10—3 1/ч, TO = 1220 ч.
закон надежности. Требуется рассчитать интенсивность отказов и частоту
отказов сети для момента времени 120 ч.
Ответ: λ=0,83⋅10-3 1/ч, α(120)=0,747⋅10-3 1/ч.
Элементы |
Количество отказов, ni |
Вес отказов по группе, n m= i n |
Время восстановления, t , мин |
Суммарное время восстановления по группе, ti , мин |
Полупроводник овые приборы | 5 | 0,25 | 80 20 22 18 23 | 163 |
Реле | 2 | 0,1 | 10 12 | 22 |
Резисторы (конденсаторы) | 10 | 0,5 | - | 230 |
Пайка | 3 | 0,15 | - | 42 |
~
Ответ: TB = 0,38 ч, KГ = 0,9989.
а(t), A(t) и To интегральных микросхем для времени их работы
t = 500, 1000, 2000 ч при условии, что параметр распределения о = 1000 ч,
время работы ИМС до отказа подчиняется закону распределения Релея.
Ответ: Результаты расчетов приведены в табл. 1.15.
Таблица 1.15. Результаты расчета к задаче 1.25
t, ч |
p(t) |
α(t)⋅10-4, 1/ч |
λ(t) ⋅10-3, 1/ч |
TO, ч |
500 |
0.88 |
4.14 |
5 |
0.0396 |
1000 |
0.606 |
6.06 |
1 |
0.0396 |
2000 |
0.135 |
2.7 |
2 |
0.0396 |
1.26 - 1.45. На испытание поставлено N0 изделий. За время t ч вышло
из строя n(t) штук изделий. За последующий интервал времени Лt вышло из
строя n(Лt) изделий. Необходимо вычислить вероятность безотказной
работы за время t и t + Лt, частоту и интенсивность отказов на интервале
Л t.
Исходные данные для решения задачи и ответы приведены в табл.1.16.
Таблица 1.16. Исходные данные и результаты расчета к задачам 1.26 – 1.45
Номер | Исходные данные | Ответы | |||||||
N0 | t | Л t | n(t) | n (Л t) | P( t) | P( t + Л t) | ~(t) | Л( t) | |
1.26 | 1000 | 8000 | 1000 | 290 | 50 | 0,71 | 0,66 | 0,00005 | 0,000073 |
1.27 | 1000 | 14000 | 1000 | 540 | 50 | 0,46 | 0,41 | 0,00005 | 0,000115 |
1.28 | 100 | 5000 | 100 | 10 | 10 | 0,9 | 0,8 | 0,001 | 0,001176 |
1.29 | 100 | 4000 | 200 | 20 | 20 | 0,8 | 0,6 | 0,001 | 0,001429 |
1.30 | 100 | 3000 | 1000 | 20 | 20 | 0,8 | 0,6 | 0,0002 | 0,000286 |
1.31 | 1000 | 800 | 100 | 216 | 15 | 0,784 | 0,769 | 0,00015 | 0,000193 |
1.32 | 1000 | 2300 | 100 | 417 | 13 | 0,583 | 0,57 | 0,00013 | 0,000225 |
1.33 | 1000 | 1200 | 100 | 274 | 14 | 0,726 | 0,712 | 0,00014 | 0,000195 |
1.34 | 1000 | 900 | 100 | 231 | 14 | 0,769 | 0,755 | 0,00014 | 0,000184 |
1.35 | 1000 | 10000 | 1000 | 370 | 40 | 0,63 | 0,59 | 0,00004 | 0,000066 |
1.36 | 1000 | 15000 | 1000 | 590 | 40 | 0,41 | 0,37 | 0,00004 | 0,000103 |
1.37 | 1000 | 21000 | 1000 | 840 | 50 | 0,16 | 0,11 | 0,00005 | 0,00037 |
1.38 | 1000 | 11000 | 1000 | 410 | 40 | 0,59 | 0,55 | 0,00004 | 0,00007 |
1.39 | 1000 | 1300 | 100 | 288 | 13 | 0,712 | 0,699 | 0,00013 | 0,000184 |
1.40 | 1000 | 1900 | 100 | 368 | 12 | 0,632 | 0,62 | 0,00012 | 0,000192 |
1.41 | 1000 | 2700 | 100 | 480 | 25 | 0,52 | 0,495 | 0,00025 | 0,000493 |
1.42 | 100 | 6000 | 500 | 50 | 20 | 0,5 | 0,3 | 0,0004 | 0.001 |
1.43 | 100 | 4000 | 1000 | 10 | 6 | 0,9 | 0,84 | 0,00006 | 0,000068 |
1.44 | 100 | 10000 | 1000 | 25 | 5 | 0,75 | 0,7 | 0,00005 | 0,000068 |
1.45 | 1000 | 1400 | 100 | 301 | 14 | 0,699 | 0,685 | 0,00014 | 0,000202 |
отказать за один год работы, если всего работало 250 реле?
Ответ: n (t) = 22.
Ответ: P (t) = 0,82, To = 5000 ч.
Ответ: P(1000) = 0,87: а(1000) = 0,25 • 10-3 1/ч; Л(1000) = 0,29 •Ю-3 1\ч
Ответ: P(150) = 0,92: α(150)=11,9⋅10-4 1/ч; λ(150)=12,9⋅10-4 1\ч
НЕВОССТАНАВЛИВАЕМЫХ НЕРЕЗЕРВИРОВАННЫХ СИСТЕМ
Задача расчета показателей надежности сложной системы состоит в
том, чтобы определить ее показатели надежности, если известны показатели
надежности отдельных элементов и структура системы, т.е. характер связей
между элементами с точки зрения надежности.
Наиболее простую структуру имеет нерезервированная система,
состоящая из n элементов, у которой отказ одного из элементов приводит к
отказу всей системы (рис.2.1). В этом случае система имеет логически
последовательное соединение элементов.
1
2
n
Рисунок 2.1. Логически последовательное соединение элементов
При расчете показателей надежности таких устройств предполагается,
что отказ элемента является событием случайным и независимым.
Тогда вероятность безотказной работы изделия в течение времени t
равна произведению вероятностей безотказной работы ее элементов в
течение того же времени. Так как вероятность безотказной работы элементов
в течение времени t можно выразить через интенсивность отказов в виде
(1.7), то расчетные формулы для вероятности безотказной работы
технического устройства при последовательном соединении элементов
можно записать следующим образом:
n
Pc(t) = Pi( t) • Р 2( t) • .•. • Pn (t) = П Pi(t),
i=1
(2.1)
PC(t) = e 0 • e о •... • e о = e i=ю
Выражения (2.1) наиболее общие. Они позволяют определить
вероятность безотказной работы изделия до первого отказа при любом законе
изменения интенсивности отказов во времени.
На практике наиболее часто интенсивность отказов изделий является
величиной постоянной. При этом время возникновения отказов подчинено
экспоненциальному закону распределения, так как для нормального периода
работы аппаратуры справедливо условие λ = const.
В этом случае выражения для количественных характеристик примут
вид:
t
-^
Pc (t ) = e-2 = eT ; 2 = X 2,.; ac (t) = ле-'; /ол = -(2.2)
i=1
Если все элементы данного типа равнонадежны, интенсивность отказов
будет
r
2c = X N-,(2.3)
i=1
где: Ni - число элементов i-го типа; r – число типов элементов.
На практике очень часто приходится вычислять вероятность
безотказной работы высоконадежных систем. При этом произведение 2ct
значительно меньше единицы, а вероятность безотказной работы P(t) близка
к единице.
В этом случае количественные характеристики надежности можно с
достаточной для практики точностью
приближенным формулам:
r
i=1
вычислить
по следующим
r
- = X Ni-i
i=1
;
(2.4)
1
T о. c . = —
2c
1
r
X N-
i=1
; ac (t) = - (1 - 2ct).
Вычисление количественных
характеристик
надежности по
приближенным формулам не дает больших ошибок для систем, вероятность
безотказной работы которых превышает 0,9, т.е. для -t < 0,1.
При расчете показателей надежности систем часто приходится
перемножать вероятности безотказной работы отдельных элементов расчета
и возводить их в степень и извлекать корни. При значениях вероятность P(t),
близких к единице, эти вычисления можно с достаточной для практики
точностью выполнить по следующим приближенным формулам:
nn
PC(t)=∏pi(t)≈1-∑qi(t),
i=1 i=1
PC(t)=pin(t)=1-nqi(t), (2.5)
n pi(t) =1-qi(t),
n
где qi (t) - вероятность отказа i-го блока.
В зависимости от полноты учета факторов, влияющих на работу
изделия, различают прикидочный, ориентировочный и окончательный расчет
показателей надежности.
ПРИКИДОЧНЫЙ РАСЧЕТ ПОКАЗАТЕЛЕЙ НАДЕЖНОСТИ
Прикидочный расчет основывается на следующих допущениях:
Прикидочный расчет показателей надежности выполняется в
следующий случаях:
Прикидочный расчет показателей надежности позволяет судить о
принципиальной возможности обеспечения требуемой надежности изделия.
ОРИЕНТИРОВОЧНЫЙ РАСЧЕТ ПОКАЗАТЕЛЕЙ НАДЕЖНОСТИ
Ориентировочный расчет показателей надежности учитывает влияние
на надежность только количества и типов, входящих в систему элементов, и
основывается на следующих допущениях:
Для определения показателей надежности изделия необходимо знать:
Выбор λi для каждого типа элементов производится по
соответствующим таблицам, приведенных в справочниках по надежности.
Таким образом, при ориентировочном расчете показателей надежности
необходимо знать структуру системы, номенклатуру применяемых элементов
и их количество.
Ориентировочный метод расчета показателей надежности используется
на этапе эскизного проектирования после разработки принципиальных
электрических схем изделий. Он позволяет наметить пути повышения
надежности изделия и производится по формулам (2.2)…(2.4).
ОКОНЧАТЕЛЬНЫЙ (ПОЛНЫЙ) РАСЧЕТ ПОКАЗАТЕЛЕЙ
НАДЕЖНОСТИ
Полный расчет показателей надежности изделия выполняется тогда,
когда известны реальные режимы работы элементов после испытания в
лабораторных условиях основных узлов и макетов изделия.
Элементы изделия находятся обычно в различных режимах работы,
сильно отличающихся от номинальной величины. Это влияет на надежность
как изделия в целом, так и отдельных его составляющих частей. Выполнение
окончательного расчета показателей надежности возможно только при
наличии данных о коэффициентах нагрузки отдельных элементов и при
38
наличии графиков зависимости интенсивности отказов элементов от их
электрической нагрузки, температуры окружающей среды и других
факторов, т.е. для окончательного расчета необходимо знать зависимости
4 = f (кнТ ,...).
Эти зависимости приводятся в виде графиков либо их можно
рассчитать с помощью, так называемых поправочных коэффициентов
интенсивности отказов ki , позволяющих учесть влияние различных факторов
на надежность изделия.
Для определения показателей надежности изделия необходимо знать:
1) число элементов с разбивкой их по типам и режимам работы;
В общем случае Ai зависит от следующих воздействующих факторов:
Знание зависимости интенсивности отказов Ai от воздействующих
факторов является необходимым для правильного использования элементов с
целью получения заданной вероятности безотказной работы за время t.
При разработке и изготовлении элементов обычно предусматриваются
определенные, так называемые «нормальные» условия работы: температура
+25±10° С, номинальный электрический режим, отсутствие механических
перегрузок и т. Д. Интенсивность отказов элементов в «нормальном» режиме
эксплуатации называется номинальной интенсивностью отказов AHi.
Интенсивность отказов элементов при эксплуатации в реальных
условиях λi равна номинальной интенсивности отказов λнi , умноженной на
поправочные коэффициенты ki , т.е.
n
λi=λнi⋅k1⋅k2⋅...⋅kn =λнi∏ki,
i=1
где: λнi - интенсивность отказов элемента, работающего в нормальных
условиях при номинальной электрической нагрузке; k1,k2,...,kn-
поправочные коэффициенты, зависящие от различных воздействующих
факторов.
Полный расчет показателей надежности применяется на этапе
технического проектирования изделия.
Пример 2.1. Система состоит из 12600 элементов. Интенсивность
отказа элемента λ= 0,32⋅10
-
6 1/ч.
Необходимо определить вероятность безотказной работы систем в
течение времени t = 50 ч.
Решение. Интенсивность отказов системы по формуле (2.3) равна:
λc =λ⋅N =0,32⋅10-6 ⋅12600=4,032⋅10-3 (1/ч).
Тогда на основании (2.2)
PC(50)=e-λct = e-4,032⋅10
3⋅50 =0,82.
Пример 2.2. Используя данные примера 2.1, вычислить среднюю
наработку до первого отказа системы.
Решение. Средняя ]
вычисляется по формуле
примера 2.1, получим
11
Τ==
O.C λc 4,032⋅10-
наработка до первого отказа системы ΤO.C
= 250 (ч).
Пример 2.3. Система состоит из трех блоков, средняя наработка до
первого отказа которых равна T1 = 160 ч, T2 = 320 ч, T3 = 600 ч. Для блоков
справедлив экспоненциальный закон надежности. Требуется определить
среднюю наработку до первого отказа системы.
Решение. Воспользуемся формулой (2.2) для средней наработки до
первого отказа системы. В нашем случае интенсивность отказов системы
равна
c 1 2 3 T1 T2 T3 160 320 600
Τ
1
O.C = А
λc
1
0,011
= 91 (ч).
Пример 2.4. Система состоит из двух устройств. Вероятности
безотказной работы каждого из них в течение времени t = 100 ч. равны:
/>1(100) = 0,95; p2(100) = 0,97. Справедлив экспоненциальный закон
распределения надежности. Необходимо найти среднюю наработку до
первого отказа системы.
Решение. Найдем вероятность безотказной работы системы по формуле
(2.1)
PC(t) = Р1( t) • Р 2( t).
Отсюда PC (100) = p1 (100) • p2 (100) = 0,95 • 0,97 = 0,92.
Найдем интенсивность отказов системы, воспользовавшись формулой
PC (100) = 0,92 = e -1 = e -1 '10°.
Из этого выражения найдем Лс -100.
1c-100 = ln0,92«0,083 или 1c = 0,83-10-3 (1/ч).
Среднее время наработки до первого отказа
1
T O. C. = —
Lc
1
0,83 -10"3
= 1200 (ч).
Пример 2.5. Вероятность безотказной работы одного элемента в
течение времени t равна p(t) = 0,9997. Требуется определить вероятность
безотказной работы системы, состоящей из N = 100 таких элементов.
Решение. Вероятность безотказной работы системы равна
PC (t) = PN (t) = (0,9997)100. Вероятность p (t) близка к единице, поэтому для
вычисления вероятности безотказной работы системы воспользуемся
формулой (2.5). В нашем случае q(t) = 1 - p(t) = 1 - 0,9997 = 0,0003. Тогда
вероятность
безотказной работы
системы
PC (t) = 1 - Nq(t) = 1 -100 • 0,0003 = 0,97 .
Пример 2.6. Вероятность безотказной работы системы в течение
времени t равна Pc (t) = 0,95. Система состоит из N = 120 равнонадежных
элементов. Необходимо найти вероятность безотказной работы элемента.
Решение. Вероятность безотказной работы элемента будет равна
pi (t) = NPc (t). Так как вероятность безотказной работы системы близка к
единице, то вычисления p(t) удобно выполнить по формуле (2.5). Для
нашего случая Qc (t) = 1 - Pc (t) = 1 - 0,95 = 0,05.
Тогда
p (t) = NPC (t) = 1 - QCt) = — = 0,9996.
C N 120
Пример 2.7. Система состоит из пяти приборов, вероятность исправной
работы которых в течение времени t = 100 ч равны: p1(100) = 0,9996;
p2(100) = 0,9998; p3(100) = 0,9996; p4(100) = 0,999; p5(100) = 0,9998.
Требуется определить частоту отказов системы в момент времени t = 100 ч.
Предполагается, что отказы приборов независимы и для них
справедлив экспоненциальный закон распределения надежности.
Решение. По условию задачи отказы приборов независимы, поэтому
вероятность безотказной работы системы равна произведению вероятностей
безотказной работы приборов. Тогда по формуле (2.5) для случая
высоконадежных систем имеем
5
Pc (t) = 1 -Z qi( t).
i=1
Определим вероятность отказа каждого блока:
q1 (100) = 1 - p1 (100) = 1 - 0,9996 = 0,0004,
q 2 (100) = 1 - p2(100) = 1 - 0,9998 = 0,0002,
q 3 (100) = 1 - p3 (100) = 1 - 0,9996 = 0,0004
q 4 (100) = 1 - p4(100) = 1 - 0,999 = 0,001,
q 5 (100) = 1 - p5 (100) = 1 - 0,9998 = 0,0002.
Тогда
Pc (100) = 1 - (0,0004 + 0,0002 + 0,0004 + 0,001 + 0,0002) = 0,9978.
Так как вероятность безотказной работы системы близка к единице, то
в соответствии с формулой (2.4) для PC (t) интенсивность отказов системы
можно вычислить из выражения
отсюда
Pc (t) = 1 - Act,
1 — Pc (t )
t
Подставляя значения Pc (100) и время t = 100 ч, получим
1 - 0,9978
100
= 2,2 -10—5 (1/ч).
Тогда частота отказов в соответствии с формулой (2.4) будет
ac (t) « Ac (1 - Act) = 2,2 -10-5(1 - 2,2 -10-5 -100) = 2,195 -10-5 (1/ч).
Пример 2.8. Все элементы электронного усилителя работают в
нормальный период эксплуатации, т.е. A = const. Усилитель должен
непрерывно работать в течение 10 ч. Из схемы известно, что усилитель
состоит из 2 ламп, 8 резисторов и 6 конденсаторов. Интенсивность отказов
всех элементов указана в табл. 2.1. Требуется рассчитать вероятность
безотказной работы системы PC (t) и среднюю наработку до первого отказа
TO.C.
Решение. Для выполнения ориентировочного расчета надежности
усилителя рассчитаем интенсивность отказов компонентов по группам, и
рассчитанные значения занесем в табл. 2.1.
Таблица 2.1. Таблица расчета интенсивности отказов к примеру 2.8
Наименование и | Обозначение по | Количество | Интенсивность Ai -10-5 1/ч | Интенсивность Ni • Ai -10-5 1/ч |
Резистор | R1, R2 | 2 | 0,09 | 0,18 |
Резистор | R3, R4 | 2 | 0,12 | 0,24 |
Резистор | R5, R6 | 2 | 0,10 | 0,20 |
Резистор | R7, R8 | 2 | 0,10 | 0,20 |
Конденсатор | C1, C2 | 2 | 0,03 | 0,06 |
Конденсатор | C3, C4 | 2 | 0,13 | 0,26 |
Конденсатор | C5, C6 | 2 | 0,09 | 0,18 |
Лампы | Л1, Л2 | 2 | 9 | 18 |
Интенсивность отказов усилителя рассчитаем по формуле (2.3):
Ac = (0,18 + 0,24 + 0,2 + 0,2 + 0,06 + 0,26 + 0,18 +18) -10-5 = 19,32 -10-5 (1/ч).
Pc (10) = e- Act = e-19,32-10 -5-10 = 0,9981
Τ
1
O.C - г
1
19,32 -10=5
- 5176 (ч).
Пример 2.9. Изделие состоит из N - 3 групп приборов. Отказы первой
группы подчинены экспоненциальному закону надежности с интенсивностью
отказов Л - 0,2 -10—3 1/ч. Отказы приборов второй группы - закону Релея с
параметром & —1000 ч и отказы приборов третьей группы - закону Вейбулла
с параметрами Л0 - 0,1 -10—4 1/ч и k —1,5. Требуется определить вероятность
безотказной работы изделия в течение времени t - 500 часов.
Решение. Найдем вероятность безотказной работы каждой группы
приборов за время t - 500 часов.
Д(500) - e-Лt - e-0,2'10-3 '500 - 0,905;
t2 (500)2
—
P2(500) - e 2& - e 2'(1000) - 0,883;
P3(500) - e-Л0tk - e-0,1'10-4-5001,5 - 0,895.
PC(500) - P1 (500) - P2(500) - P3 (500) - 0,716.
В настоящем разделе приведены задачи на расчет показателей
надежности невосстанавливаемых изделий при основном соединении
элементов. Эти задачи легко решать, используя типовые примеры,
рассмотренные в разделе 2.2.
Ответ: PA(200) - 0,27; TO.A -151,5 ч.
Ответ: PC(50) = 0,953; TO.C = 1040 ч.
Ответ: Рпр (t) = 0,93.
Ответ: Ризд (t) = 0,819.
Ответ: PK (t) = 0,595.
Исходные данные для решения задач и ответы приведены в табл. 2.2.
Таблица 2.2. Исходные данные и ответы к задачам 2.7 – 2.9
Номер | Исходные данные | Ответ | |||||
N групп | Л-10-3, | а, ч | Л0 -10-4, | k | t, ч | ||
2.7. | 3 | 0,1 | 1200 | 0,03 | 1,5 | 1000 | 0,54 |
2.8. | 2 | - | 1000 | 1,6 | 1,3 | 500 | 0,53 |
2.9. | 2 | 0,09 | - | 1,3 | 1,3 | 120 | 0,93 |
2.10. – 2.12. Система состоит из N блоков, средняя наработка до
первого отказа которых равна To.i . Для блоков справедлив экспоненциальный
закон надежности. Требуется определить среднюю наработку до первого
отказа системы.
Исходные данные для решения задач и ответы приведены в табл. 2.3.
Таблица 2.3. Исходные данные и ответы к задачам 2.10 – 2.12
Номер | Исходные данные | Ответ TO.C , ч | |||||
N блоков | TO1 , ч | TO2 , ч | TO3 , ч | TO4 , ч | TO5 , ч | ||
2.10. | 3 | 150 | 750 | 500 | - | - | 100 |
2.11. | 4 | 1600 | 1800 | 2000 | 2200 | - | 468 |
2.12. | 5 | 20 | 30 | 40 | 50 | 60 | 6,9 |
2.13 – 2.15. Система состоит из N блоков. Вероятность безотказной
работы каждого блока в течение времени t равна pi (t) . Для блоков
справедлив экспоненциальный закон надежности. Требуется определить
среднюю наработку до первого отказа системы.
Исходные данные для решения задач и ответы приведены в табл. 2.4.
Таблица 2.4. Исходные данные и ответы к задачам 2.13 – 2.15
Номер | Исходные данные | Ответ TO.C , ч | ||||||
N | t, ч | p1(t) | p2(t) | p3(t) | p4(t) | p5(t) | ||
2.13. | 3 | 1000 | 0,97 | 0,98 | 0,96 | - | - | 10640 |
2.14. | 5 | 240 | 0,9 | 0,8 | 0,85 | 0,7 | 0,75 | 210 |
2.15. | 4 | 10 | 0,94 | 0,95 | 0,97 | 0,98 | - | 62,5 |
2.16 – 2.18. Система состоит из N элементов. Вероятность безотказной
работы одного элемента в течение времени t равнаp(t) . Требуется
определить вероятность безотказной работы системы.
Исходные данные для решения задач и ответы приведены в табл. 2.5.
Таблица 2.5. Исходные данные и ответы к задачам 2.16 – 2.18
Номер задачи | Исходные данные | Ответ PC(t) | |
p(t) | N элементов | ||
2.16 | 0,9999 | 1000 | 0,9 |
2.17 | 0,9998 | 50 | 0,99 |
2.18 | 0,9996 | 100 | 0,96 |
2.19 – 2.21. Вероятность безотказной работы системы в течение
времени t равна PC (t). Система состоит из N равнонадежных элементов.
Необходимо найти вероятность безотказной работы элемента.
Исходные данные для решения задач и ответы приведены в табл. 2.6.
Таблица 2.6. Исходные данные и ответы к задачам 2.16 – 2.18
Номер задачи | Исходные данные | Ответ p (t) | |
PC(t) | N элементов | ||
2.19. | 0,97 | 200 | 0,99985 |
2.20. | 0,95 | 300 | 0,99983 |
2.21. | 0,98 | 1000 | 0,99998 |
2.22 – 2.24. В изделии могут быть использованы только те элементы,
интенсивность отказов которых равна А. Изделие имеет число элементов N.
Требуется определить среднюю наработку до первого отказа и вероятность
безотказной работы в конце первого часа.
Исходные данные для решения задач и ответы приведены в табл.2.7.
Таблица 2.7. Исходные данные и ответы к задачам 2.16 – 2.18
Номер задачи | Исходные данные | Ответ | ||
А, 1/ч | N элементов | PC(t) | TO.C , ч | |
2.22. | 1-10-5 | 500 | 0,995 | 200 |
2.23. | 1-10-5 | 2500 | 0,975 | 40 |
2.24. | 1-10-3 | 100 | 0,905 | 10 |
надежности. Определить среднюю наработку до первого отказа системы.
Ответ: Toc = 1223,2 ч.
Ответ: Рс(200) = 0,52.
Ответ: Рс(100) = 0,247, TO.C = 72,97 ч.
отказа системы в целом, если для приборов каждого типа справедлив
экспоненциальный закон надежности.
Ответ: Toc = 101 ч.
Ответ: То = 22410 ч
Ответ: Q(10) = 0,259
Ответ: Pc (100) = 0,333
Ответ: q(10) = 0,3-10-5
Ответ: n=5
1/ч, λ2 = 0,3 ⋅ 10-3 1/ч, λ3 = 0,8 ⋅ 10-3 1/ч. Закон надежности считаем
экспоненциальным. Найти среднюю наработку до первого отказа устройства.
λ3
Рисунок 2.2. Схема расчета надежности (к задаче 2.34).
Ответ: TO.C = 625 ч
Расчет следует выполнить по данным о надежности элементов,
приведенным в приложении 2.
Исходные данные для решения задачи приведены в табл.2.8.
Таблица 2.8. Исходные данные к задачам 2.35 – 2.40
Номер | Резисторы | Конденсаторы | Транзисторы | Время | |||
тип | шт | тип | шт | тип | шт | t ,ч | |
2.35 | ВС-0,25 | 20 14 | Слюдяные | 10 12 | Германиевые | 5 4 | 650 |
2.36 | МЛТ-0,5 МЛТ-1 МЛТ-2 | 32 19 | Керамические | 8 15 | Германиевые | 4 11 | 260 |
2.37 | МЛТ-0,25 | 40 | Керамические | 30 49 | Кремниевые | 10 7 | 5000 |
2.38 | МЛТ-0,25 МЛТ-0,5 | 21 5 3 | Слюдяные | 8 | Германиевые | 6 4 | 2000 |
2.39 | МЛТ-0,25 МЛТ-0,5 | 13 10 | Керамические | 18 | Кремниевые | 5 9 | 250 |
2.40 | МЛТ-0,25 | 25 2 | Керамические | 27 | Германиевые | 8 | 500 |
высокочастотные низкочастотные | 12 |
Ответ: Результаты расчетов приведены в табл. 2.9.
Таблица 2.9. Результаты расчетов к задачам 2.35 -2.40
Номер | 2.35 | 2.36 | 2.37 | 2.38 | 2.39 | 2.40 |
P(t) | 0,9516 | 0,8688 | 0,3719 | 0,8772 | 0,9801 | 0,9234 |
TO, ч | 13089 | 1848 | 5056 | 15267 | 12422 | 6274 |
Ответ: t = 1050 ч
Ответ: «(1000) = 0,37 -10-3 1/ч
Ответ: Q (100) = 0,01, Q (300) = 0,03
Ответ: P (200) = 0,96,
Ответ: t = 22,95 ч
Исходные данные для решения задачи и ответы приведены в табл.2.10
Таблица 2.10. Исходные данные и ответы к задачам 2.50 -2.55.
Номер | Исходные данные | Ответы | ||
n , шт | t, ч | PC (t) | λ, 1/ч | |
2.50 | 3 | 2000 | 0.98 | 3,3⋅10-6 |
2.51 | 5 | 500 | 0,997 | 1,2 ⋅10-6 |
2.52 | 2 | 100 | 0,996 | 0,2 ⋅10-4 |
2.53 | 10 | 10000 | 0,999 | 0,1⋅10-7 |
2.54 | 4 | 1000 | 0,998 | 5⋅10-8 |
Резервирование применение дополнительных средств и (или)
возможностей с целью сохранения работоспособного состояния изделия при
отказе одного или нескольких его элементов.
В этом случае отказ наступает только после отказа основных и всех
резервных элементов. При этом возможно резервирование на уровне всей
системы в целом (общее резервирование) или на уровне отдельных ее
элементов (раздельное резервирование).
На практике применяются способы резервирования, приведенные на
рис. 3.1. Схемные реализации различных способов резервирования показаны
на рис. 3.2.
Рисунок 3.1. Способы резервирования
Общее резервирование - резервирование, при котором резервируемым
элементом является все изделие в целом (рис.3.2,а). Раздельное
резервирование — резервирование, при котором резервируемыми являются
отдельные элементы изделия или их группы (рис.3.2,б).
Основным параметром резервирования является его кратность.
Кратность резерва - отношение числа резервных элементов изделия к числу
резервируемых ими основных элементов изделия, выраженное
несокращенной дробью.
В зависимости от кратности резервирование подразделяется на
резервирование с целой и дробной кратностью. Схемные обозначения обоих
видов резервирования при постоянном включении резерва одинаковы. Для
их различия на схеме указывается кратность резервирования m.
При резервировании с целой кратностью величина m есть целое число,
при резервировании с дробной кратностью – дробное несокращаемое число.
2
1
1
2
n
n
a)
б)
2
1
2
n
n
2
n
е)
k
m = —
n
1
г)
в)
к
n - k
m=---
k
д)
Рисунок 3.2. Схемные реализации различных способов резервирования
а) общее постоянное с целой кратностью; б) раздельное постоянное с
целой кратностью; в) общее замещением с целой кратностью; г)
раздельное замещением с целой кратностью; д) общее постоянное с
дробной кратностью; е) раздельное замещением с дробной кратностью.
4
Например, m = — означает наличие резервирования с дробной кратностью,
при котором число резервных элементов равно четырем, число основных –
двум, а общее число элементов – шести. Сокращать дробь нельзя, так как
4
если m = — = 2, то это означает, что имеет место резервирование с целой
кратностью, при котором число резервных элементов равно двум, а общее
число элементов равно трем.
По способу включения резервирование разделяется на постоянное и
резервирование замещением. Постоянное резервирование – резервирование
без перестройки структуры изделия при возникновении отказа его
элемента. Резервированием замещением – резервирование, при котором
функции основного элемента передаются резервному только после отказа
основного элемента.
При включении резерва по способу замещения резервные элементы до
момента включения в работу могут находиться в трех состояниях:
Приведем основные расчетные формулы для указанных выше видов
резервирования.
Pc( t) = 1 -
-| m+1
n
1-П Pi( t)
i=1 _
,
(3.1)
где pi (t) - вероятность безотказной работы i-го элемента в течение времени
t; n – число элементов основной системы или любой из резервных систем; m
– кратность резервирования (число резервных цепей).
При экспоненциальном законе распределения надежности, когда
Pi(t ) = e'"t,
PC (t) = 1 - (1 - e- ^)m+1, (3.2)
mm
(3.3)
(3.4)
1 c (t) =
1 o (m + 1)e - 1 ot (1 - e ~ 1 ot m
1 -(1 - e-1 ot m+1
n
где Ao = ^ Ai - интенсивность отказов основной системы или любой из
i=1
резервных систем; TO - среднее время наработки до первого отказа основной
системы или любой из резервных систем.
При резервировании неравнонадежных изделий
mm
Pc(t) = 1 - п qi (t) = 1 - П [1 - Pi (t)], (3.5)
i=0 i=0
где qi(t), pi (t) - вероятность отказов и вероятность безотказной работы i-го
изделия в течение времени t соответственно.
n
Pc (t) = П {1 - [1 - Pi( t)] mi }, (3.6)
i=1
где pi (t) - вероятность безотказной работы i-го элемента в течение времени
t; mi - кратность резервирования i-го элемента; n – число элементов основной
системы.
При экспоненциальном законе распределения надежности, когда
Pi (t) = e ,
n
Pc (t) = П{1 - [1 - e-Ai" ] mi+1}. (3.7)
i=1
При равнонадежных элементах и одинаковой кратности
резервирования
Pc(t ) = 1 -(1 - e 1 tm
т - fp^-^ m 1
0 A(m + 1) i=o vi(vi + 1)...(vi + n - 1)
AC (t ) =
n (m + 1)Ae- At (1 - e ~ At m
1 -(1 - e ’At m+1
(3.8)
(3.9)
(3.10)
i +1
где Л - интенсивность отказа одного элемента системы; v =----
m +1
.
При экспоненциальном законе надежности и ненагруженном состоянии
резерва
Pc (t ) = e - Л f &£,
i=0 i!
I Ло W>t)m
Л (t )= mU^i
m!Y WL
i'
i=0
(3.11)
(3.12)
(3.13)
где Ло, To - интенсивность отказа и средняя наработка до первого отказа
основной (нерезервированной) системы.
При экспоненциальном законе надежности и облегченном состоянии
резерва
Pc (t) = e - wt
ma
1+ f 7(1 - e
i=1 i!
. - wit
,
(3.14)
m
TOC = T f —7,
Ло i=ol + ik
(3.15)
i-1 ( 1 \
где ai = П j+у ;
j=0k Л1 J
k = i; Л
Ло 1
интенсивность отказов резервного
устройства до замещения.
При нагруженном состоянии резерва формулы для PC (t) и TO.C
совпадают с (3.2) и (3.3) соответственно.
(рис.3.2,г):
n
Pc (t )=n Pi(t) (3.16)
i=1
где Pi (t) - вероятность безотказной работы системы из-за отказов элементов
i-го типа, резервированных по способу замещения. Вычисляется Pi (t) по
формулам общего резервирования замещением [формулы (3.11), (3.14)].
При равной надежности всех элементов:
к C (t ) =
E i
к o (кt) m
m
m! E
i=0
(Mm
i!
,
(3.17)
(3.18)
n - k
Pc (t) = E CP—i (t)' [1 — P (t )]i, (3.19)
i=0
где P (t) - вероятность безотказной работы основной или любой из резервных
систем; n – общее число элементов расчета резервированного соединения, k –
число элементов необходимое для нормальной работы, n - к - число
резервных элементов. Кратность резервирования m =
n - к n!
----, Ci =--------
k n i!(n - i)!
число комбинаций из n элементов по i.
При экспоненциальном законе распределения функции надежности
PC(t)= E'Cin--kn-i)t (1 - e-^t) ,
i=0
__ i n - к i
к i=0 к + I
(3.20)
(3.21)
.
Pc(t) = e-"ot E^ (3.22)
i=0 i!
Этот вид резервирования применяется, если все элементы системы
выполняют одинаковые функции. Основная система имеет n элементов, а k
k
элементов находятся в холодном резерве. Кратность резервирования m = —.
n
Надежность данной системы равна надежности системы с общим
резервированием с замещением, но в то же время имеет в n раз меньше
резервных элементов. Однако переключающие устройства при этом
усложняются.
Вероятность безотказной работы резервированной системы можно
вычислить, суммируя вероятности всех благоприятных гипотез, т.е.
k
Pc (t) = Z Pj( t), (3.23)
J=1
где p j (t) - вероятность j-й благоприятной гипотезы, k – число
благоприятных гипотез.
При расчете показателей надежности дискретных устройств
управления требуется учитывать характер возникающих в них
неисправностей (например, в безопасных системах). Поэтому в таких
устройствах отдельно подсчитывается вероятность отсутствия ложного
сигнала «0» - Po(t) и вероятность отсутствия ложного сигнала «1» - Д(t). На
рис.3.3 показана избыточная структура, содержащая n идентичных каналов
передачи информации.
Рис.3.3 Избыточная структура
Число n называется избыточностью структуры. На входы каналов
поступают двоичные переменные xi...xn. Выходы каналов подключены к
входам одного восстанавливающего органа (ВО). Такая структура называется
(n,1) - структурой. Если все элементы структуры исправны, то сигналы yi на
выходах каналов равны и на выходе восстанавливающего органа
формируется сигнал
Если в одном или нескольких каналах возникает неисправность, то на
выходах появляются ошибки типа 0 ^ 1 или 1 ^ 0. Восстанавливающий
орган корректирует (исправляет) определенное число ошибочных сигналов
yi и реализует в этом случае некоторую пороговую функцию:
y = Mnd = f (У1, y 2- yn) =
1, апёе fy
i=1
0, апёе Xy/
i =1
3.24
где n - избыточность (число входов восстанавливающего органа), yi -
значение логической переменной, d - порог восстанавливающего органа.
Q , n +1 _ _ _ ~ _
При нечетном n, если порог d = — восстанавливающий орган
реализует мажоритарную логическую функцию.
При расчете показателей надежности (n ,1) - структур (рис.3.3) будем
предполагать, что ВО абсолютно надежен и известны pо(t) - вероятность
отсутствия в каналах ошибки типа 1 ^ 0 и Р1(t) - вероятность отсутствия в
каналах ошибки типа 0 ^ 1. Тогда вероятности отсутствия ложных сигналов
«0» и «1» можно определить по формулам:
/Дt) = Cd ■ (1 — p0(t))n—d + Cn+1 ■ (p0(t))d+1 • (1 — p0(t))n—(d+1) + -+Cn ■ (p0(t))n =
n
= Xcn ■ (p0(t))' ■ (1 — p0(t))n—‘ 3.25
i=d
P1(t) = Cn—(d—11 ■ (p1(t))n—(d —" ■ (1 — p1(t))d—1 + Cn—(d—2) ■ (p1(t))n—(d—2) ■ (1 - p1(t))d—2 +
n
+... + Cn ■ (p1(t))n = X Cn ■ (p1(t))i ■ (1 — p1(t))n-i 3.26
i=n—(d—1)
Часто для расчетов удобнее пользоваться выражениями для
вероятности отказов ^0, Q1. Эти выражения получают по аналогии с
формулами (3.25) и (3.26):
n
Q(t) = XCn ■ (q0(t))i ■ (1—q0(t))n—‘ 3.27
i=n—(d—1)
n
Q1(t) = X Cn ■ (q1(t))' ■ (1—q1(t))n — ‘ 3.28
i=d
Реально надежность (n ,1) - структур зависит от надежности ВО. Так
как ВО является последним сравнивающим элементом структуры, то
надежность (n ,1) - структуры в целом не может быть выше надежности ВО.
В приложении 3 приведены выражения для расчета показателей
надежности некоторых ( n,1 ) - структур.
Пример 3.1. Дана система, схема расчета показателей надежности
которой изображена на рис.3.4. Необходимо найти вероятность безотказной
работы системы при известных вероятностях безотказной работы ее
элементов (значения вероятностей указаны на рисунке).
а б вг
Рисунок 3.4. Схема расчета показателей надежности (к примеру 3.1).
Решение. Из рисунка 3.4 видно, что система состоит из двух (I и II)
неравнонадежных устройств.
Устройство I состоит из четырех узлов:
а – дублированного узла с постоянно включенным резервом, причем
каждая часть узла состоит из трех последовательно соединенных (в смысле
надежности) элементов;
б – дублированного узла по способу замещения;
в – узла с одним нерезервированным элементом;
1
Устройство II представляет собой нерезервированное устройство,
надежность которого известна.
Так как оба устройства неравнонадежны, то на основании формулы
(3.5) имеем
m
PC(t)=1-∏(1-pi(t))=1-[1-pI(t)]⋅[1-pII(t)].
i=0
Найдем вероятность pI (t). Вероятность безотказной работы
устройства I равна произведению безотказной работы всех узлов, т.е.
Pl(t) = Ра (t) • Рб (t) • Рв (t) • Рг (t) .
В узле а число элементов основной и резервной цепи n = 3, а кратность
резервирования m = 1, тогда на основании формулы (3.1), имеем
Ра (t) = 1 -
3
1 -П Pi(t)
i=1
2
= 1 - (1 - 0.93]2 = 0.93.
В узле б кратность общего резервирования замещением m = 1, тогда на
основании формулы (3.11) имеем
рб (t) = e -A-t ]m ("t)- = e -A-t (1 + A-1) = 0.9 • (1 + 0,1) = 0.99.
i=0 i!
В узле г применено резервирование с дробной кратностью. В этом
случае число основных и резервных систем n = 3, число систем,
необходимых для нормальной работы k = 2. Тогда на основании формулы
(3.19)
п _ к 1-5
Рг(t) = XC'nPn-1 (t)-[1 -р(t)]' = Xс3р '(t)■ [1 -р(t)]i = 3p2(t)-2p3(t) =0,972
i=0 i=0
Вероятность безотказной работы устройства I будет
Pi (t) = Ра (t) - Рб (t) - Рв (t) ■ рг (t) = 0,93 ■ 0,99 - 0,97 ■ 0,972 = 0,868.
Тогда вероятность безотказной работы резервированной системы будет
m
Pc (t) = 1 -П (1 - Pi (t)) = 1 - [1 - Pi (t)] ■ [1 - Pii (t)] = 1 - (1 - 0,868)(1 - 0,9) = 0,987.
i=0
Пример 3.2. Вероятность безотказной работы преобразователя
постоянного тока в переменный в течение t = 1000 ч равна 0,95. Для
повышения надежности системы электроснабжения на объекте имеется такой
же преобразователь, который включается в работу при отказе первого.
Требуется рассчитать вероятность безотказной работы и среднюю наработку
до первого отказа системы, состоящей из двух преобразователей.
Решение. В данном случае имеет место общее резервирование
замещением кратности m=1. Для расчета вероятности безотказной работы
воспользуемся формулой (3.11):
Pc (t) = e- Xm ^ = e" "0t (1 + Aot).
i=0 i!
По условию задачи вероятность безотказной работы основной системы
PO (1000) = e-Aot = 0,95, тогда Aot = ln0,95 = 0,05. Подставив эти значения в
формулу, получим:
PC (1000) = 0,95(1 + 0,05) = 0,9975.
Среднюю наработку до первого отказа системы рассчитаем по формуле
(3.12):
Так как в течение времени t = 1000 ч Aot = 0,05, то
Ао = ^“5 = 0’^5 = 0,5 -10-4 (1/ч), а средняя наработка до первого отказа
11
нерезервированного преобразователя То = — =-----
Ао 0,5 -10
—
т = 20000 (ч).
Тогда средняя наработка до первого отказа резервированной системы
То. C = 2То = 40000 (ч).
Пример 3.3. Система состоит из 10 равнонадежных элементов, средняя
наработка до первого отказа элемента То = 1000 ч. Предполагается, что
справедлив экспоненциальный закон надежности для элементов системы и
основная и резервная системы равнонадежны. Необходимо найти среднюю
наработку до первого отказа То с системы, а также частоту отказов ас (t) и
интенсивность отказов Ас (t) в момент времени t = 50 ч в следующих
случаях:
Решение. По условию задачи справедлив экспоненциальный закон
надежности для элементов, поэтому средняя наработка до отказа основной
системы
1 1 1 TO 1000
То о — То 1 — — — — — —100 (ч).
I Ai
i—1
Тогда на основании формулы (3.3) для средней наработки до первого
отказа при постоянно включенной одной резервной системе имеем
m1
TO.C = TO.2 = TO.O I . = TO.O
i=0 i + 1
\
1 + ' ] = 150 (ч).
2 7
При дублировании системы по методу замещения
TO. C — TO .3 — To . o (m +1) — 2То O — 200 (ч).
В случае нерезервированной системы интенсивность отказов не
зависит от времени и равна сумме интенсивностей отказов элементов.
Найдем интенсивность и частоту отказов системы в момент времени
t = 50 ч для случая 1:
10 10 1 1
А(50) = 2 2-=2 = ^ = 0,01 (1/ч).
i=1 i=1 TO. i TO. O
«1(50) = 21(50) • P (50) = 21(50) • e -21(5°)•5° = 0,01 • e ■°,°1•5° = 6-10-3 (1/ч).
В случае дублированной системы интенсивность и частота отказов
могут быть найдены по известной вероятности безотказной работы системы.
В рассматриваемом случае число элементов нерезервированной системы
n = 10, кратность резервирования m = 1. Тогда на основании формул (3.2) и
(3.11) имеем
P2 (t) = 1 - (1 - e-2ot)m+1 = 2e~2ot - e~22ot,
P3 (t) = e2 2(2)- = e-2 (1 + 20t)
i=0 i!
Найдем интенсивность и частоту отказов для случаев 2 и 3:
«2 (t) = — Pt) = 220e -20 (1 - e - 201),
dt
1 M = O(t) 22„e -2ot (1 - e ) 22O (1 - e2)
a3 (t) = -Pt) = 22 • t • e2,
Подставляя в полученные выражения исходные данные, будем иметь:
a 2(50) = 4,8 •Ю-3 (1/ч); 22(50) = 5,7 •Ю-3 (1/ч);
a3(50) = 3 • 10-3 (1/ч); 23(50) = 3,3 3 40-3 (1/ч);
Пример 3.4. Для повышения надежности усилителя все его элементы
дублированы. Предполагается, что все элементы подвержены лишь одному
виду отказов и последствия отказов отсутствуют. Необходимо найти
вероятность безотказной работы усилителя в течение t = 5000 ч. Состав
элементов нерезервированного усилителя и данные о интенсивности отказов
его элементов приведены в табл.3.1.
Таблица 3.1. Данные о интенсивности отказов элементов ( к примеру 3.4)
Элементы | Количество | Интенсивность отказов элементов A , |
Транзисторы | 1 | 2,16 |
Резисторы | 5 | 0,23 |
Конденсаторы | 3 | 0,32 |
Диоды | 1 | 0,78 |
Катушки индуктивности | 1 | 0,09 |
Решение. В данном случае имеет место раздельное резервирование с
кратностью m = 1, число элементов нерезервированной системы n = 11.
Тогда используя данные табл.3.1, на основании формулы (3.7) получим:
PC (5000) = П {1 - [1 — e-Ai '500° ]2}.
i=1
Так как Ai << 1, то для приближенного вычисления показательную
функцию можно разложить в ряд и ограничиться первыми двумя членами
разложения 1 - e - Ai"5°°° = 5000Аг-. Тогда
11 11 11
PC (5000) = П{1 — [1 — 5000Ai ]2} = 1 -^ (5000A)2 = 1 - 50002 ^ A2 = 0,985
i=1 i=1 i=1
Пример 3.5. Схема расчета показателей надежности устройства
приведена на рис.3.5:
Рисунок 3.5. Схема расчета показателей надежности (к примеру 3.5).
Предполагается, что последствие отказов отсутствует и все элементы
расчета равнонадежны. Интенсивность отказа элементов А = 1,35 • 10-31/ч.
Требуется определить наработку до первого отказа резервированной
системы.
Решение. В данном случае имеет место раздельное резервирование
равнонадежных устройств с постоянно включенным резервом. Число
элементов нерезервированной системы n=2, кратность резервирования m=1.
Для вычисления средней наработки до первого отказа воспользуемся
формулой (3.9):
т (n -1)! Л 1 = 1 f 1 .
O.C Л(m + 1)“0 v, (v, + 1)...( v, + n - 1) 2Л=о v (v^ 1)
Так как
i +1 i +1
m + 1. = T,
то v0 =
1
2, V1 = 1 ’
Тогда
T 14 1
TO C = 22 3 + 2
1
11
122
11
12 -1,35 -10-3
= 680 (ч)
Пример 3.6. Схема расчета показателей надежности приведена на
рис.3.6. Интенсивности отказов элементов имеют следующие значения:
21 = 0,23-10-3 1/ч: 22 = 0,5-10-4 1/ч: 23 = 0,4-10-3 1/ч. Предполагаем, что
последствие отказов элементов отсутствует. Необходимо найти среднюю
наработку до первого отказа устройства.
Рисунок 3.6. Схема расчета показателей надежности (к примеру 3.6).
Решение. Готовой формулы для средней наработки до первого отказа в
рассматриваемом случае нет. Поэтому необходимо воспользоваться
соотношением:
Tc Ч Pc (')dt ■
0
Найдем выражение для вероятности безотказной работы устройства.
Очевидно PC (t) = pI (t) - pI (t) - рш (t), где
Pi (t) = 1 - [1 - p,( t )]2 = 2 A( t) - p2( t);
Piii (t) = 1 - [1 - p з( t )]2 = 2 p з (t) - p 2 (t) ■
Тогда, подставляя значения pI (t) и pIII (t) в выражение для Pc (t),
получим
Pc (t) = [2 P1 (t) - p2 (t)] - p 2 (t) - [2 p з (t) - p з2 (t)] =
Так как p1 (t) = e 21 t; p2 (t) = e 22 t; p3 (t) = e 23 t, то
p (t) = 4e-(2i+ 22+ 2з)t - 2e-(22з+ 22+ 2i)t - 2e_(22i+ 22+ 2з)t + e-(22i+ 22+ 223)t
•
To. C
^
= J PC (t) dt =
0
4
2 2
—
1
+--.
221 + 22 + 223
Подставляя в выражение для To с значение интенсивностей отказов из
условия задачи, получаем
To. с =
4 2 2
10-3(0,23 + 0,05 + 0,4) 10-3(0,8 + 0,05 + 0,23) 10-3(0,46 + 0,05 + 0,4)
---з------1-----------= 2590 (ч)
10-3(0,46 + 0,05 + 0,8)
Пример 3.7. Система имеет кратность общего резервирования m=5.
Основная нерезервированная система содержит четыре равнонадежных
элемента с логически последовательным соединением. Интенсивность отказа
одного элемента 2 = 0,2 -10-3 (1/ч). Определить характеристики надежности
системы за 1000 ч работы.
Решение. Определим интенсивность отказов основной системы по
формуле:
2o = n • 2 = 42 = 4 • 0,2-10-3 = 0,8-10-3 (1/ч).
Вероятность безотказной работы системы определим по формуле (3.2):
PC (t) = 1 - (1 - e- 2ot) m+1 = 1 - (1 - e -0,8 )6 = 0,972.
Среднее время наработки до первого отказа и интенсивность отказов
системы соответственно равны
m
2o i=0 i + 1
1
0,8 •Ю"3
(
1 +
11111Л
2 3 4 5 6 J
= 3062
(ч).
2 C (t) =
2 o (m + 1)e "2 ot (1 - e ~ 2 ot m
1 -(1 - e'2 ot m+1
0,8 •10-3 •6 e ~ 71 -e"0,8)5 = 0.11 -10-3
1 - (1 - e "°,8)6
Пример 3.8. Вычислительная система построена из 500 однотипных
блоков с интенсивностью отказа 2 = 0,3 -10 6 1/ч. В скользящем холодном
резерве находятся пять таких же блоков, которые могут заменить любой из
отказавших блоков. Определить показатели надежности системы за 10000
часов работы.
Решение. Определим показатели надежности системы, используя
формулы
Pc (t ) = e~Aot f
i=0 i!
-zA"
i=0 i!
O.C O
Ao
Определим интенсивность отказов основной системы за время
t = 10000 ч.
Aot = nAt = 0,3 -10
Тогда
—I
6 • 500 • 10000 = 1,5.
PC (10000) = e—1,5
2! 3! 4! 5!
= 0,2231 • 4,4617 = 0,9954
A = 1,5 -10 (1,5) = 0,21^10-6 (i/ч)
C 5!-4,4617
Toc = 5 +1 4 = 40000 (ч).
1,5 -10-4
—
Пример 3.9. Для повышения точности измерения некоторой величины
применена схема группирования приборов из пяти по три, т.е. результат
измерения считается верным по показанию среднего (третьего) прибора.
Требуется найти вероятность и среднее время наработки до первого отказа
такой системы, если интенсивность отказов каждого прибора A = 0,4 -10-3
1/ч, последствие отказов отсутствует, а время, в течение которого система
измерения должна быть исправна, t = 500 ч.
Решение. Решим эту задачу двумя способами.
Способ 1. В данном случае измерительная система отказывает в том
случае, если откажут из приборов три и более, т.е. имеет место общее
резервирование дробной кратности, когда общее число приборов n = 5, число
приборов, необходимых для нормальной работы k = 3, а кратность
2
резервирования m = —
.
Используя формулу (3.19), получим
n - k 2 с
Pc (t) = Z C'nPn - i (t) • [1 — P (t )]i = Z C5p5—' (t)' [1 — P (t )]i = 6 P 5( t) —15 p 4( t) +10 p 3( t)
i=0 i=0
Так как вероятность безотказной работы одного прибора в течение времени
t = 500 ч будет
p (500) = e -2 = e "0’4d0 -3 '500 = 0,82, то
PC (500) = 6p5(t) -15p4(t) +10p3(t) = 15 • 0,825 -15 • 0,824 +10 • 0,823 = 0,95.
Средняя наработка до первого отказа системы вычисляется по формуле
(3.21):
T = 1n-k 1 = 1|J_ = 1 (1 +1 +1V 47
O.C Л Z0k + i 2^03 + i A3 4 5J 602
47
60 • 0,4 • 10-3
= 1958 (ч)
Способ 2. Благоприятными ситуациями в рассматриваемом случае
являются следующие:
Вероятность безотказной работы измерительной системы на основании
(3.23) будет
Pc (t) = p1( t) + p 2( t) + p 3( t) = 6 p 5( t) -15 p 4( t) +10 p3, что совпадает с
ответом, полученным при первом способе решения.
Средняя наработка до первого отказа вычисляется в данном случае
путем интегрирования Pc (t) по всей временной оси:
^ ^
To.с = J Pc (t) = J [6p 5(t) -15 p 4(t) +10p3 (t)]dt
0 0
TO .C
z z z
= 6 J e~52tdt -15 J e-42tdt +10 J e~2tdt =
0 0 0
f 6
15
15 10 Л
4 + 3 /
1
47
602
что
совпадает с ответом, полученным при первом способе решения.
Пример 3.10. Система, представляющая мажоритарную (3,1) –
структуру, состоит из трех логически последовательно соединенных блоков,
которые образуют три вычислительных канала. Известны вероятности
отсутствия ложного сигнала «0» каждого из блоков, которые соответственно
равны: p 01 = 0,93, p02 = 0,95, p 03 = 0,98. Система считается
работоспособной, если работоспособны два ее канала. Найти вероятность
отсутствия ложного сигнала «0» системы в целом, если ВО абсолютно
надежный.
Решение. Вероятность отсутствия ложного нуля для (3,1) – структуры,
в которой логическая формула ВО (система голосования) M 32 , равна:
p,( t) = 3 p 02( t) - 2 p 0( t).
Учитывая, что канал образуют три последовательно соединенных
элемента, получим окончательную формулу расчета вероятности отсутствия
ложного сигнала «0» системы:
P0.C (t) = 3 • [p01(t) • p02 (t) • p03 (t)]2 - 2 • [p01(t) • p02 (t) • p03 (t)]3 •
Подставив значения вероятностей отсутствия ложного сигнала «0»
каждого блока в выражение для P0.C (t) получим:
P0.C(t) = 3(0,93 • 0,95 • 0,98)2 - 2(0,93 • 0,95 • 0,98)3 = 0,9506
Пример 3.11. Система, представляющая мажоритарную (3,1) –
структуру, состоит из трех логически последовательно соединенных блоков,
которые образуют три канала передачи информации. Известны вероятности
отсутствия ложного сигнала «1» каждого из блоков, которые соответственно
равны: p11 = 0,95, p12 = 0,96, pn = 0,98. Система считается
работоспособной, если работоспособны два ее канала. Найти вероятность
отсутствия ложного сигнала «1» системы в целом, если вероятность
безотказной работы ВО равна p^ = 0,999.
Решение. Вероятность отсутствия ложной единицы для (3,1) –
структуры, в которой логическая формула ВО M 32 , равна:
P( t) = 3 p2( t) - 2 p3( t).
Учитывая, что канал образуют три последовательно соединенных
элемента, а также последовательное соединение ВО, получим окончательную
формулу расчета вероятности отсутствия ложного сигнала «1» системы
Pi.с(t) = {3 • [P11 (t) • P12(t) • P13(t)]2 - 2 •[P11 (t) • P12(t) • P13(t)]3} • Pai .
Подставив значения вероятностей отсутствия ложного сигнала «1»
каждого блока и значение вероятности безотказной работы ВО в выражение
для Pq. c (t) получим:
P1.C(t) = {3(0,95 • 0,96 • 0,98)2 - 2(0,95 • 0,96 • 0,98)3} • 0,999 = 0,9675
Рисунок 3.7. Схема расчета показателей надежности (к задаче 3.1).
Ответ PC (t) = 0,997
Рисунок 3.8. Схема расчета показателей надежности (к задаче 3.2).
Ответ PC (t) = 0,991, QC (t) = 0,009
Рисунок 3.9. Схема расчета показателей надежности (к задаче 3.3).
Ответ Pc (t) = 0,95,
p1 p2
p3 p4
Рисунок 3.10. Схема расчета показателей надежности (к задаче 3.4).
Ответ Pc (t) = 0,944,
Ответ Toc = 1200 ч.
Рисунок 3.11. Схема расчета показателей надежности (к задаче 3.4).
Ответ Pc (t) = 0,998
работы изделия в течение времени t = 100 ч и среднюю наработку до первого
отказа, если для элементов справедлив экспоненциальный закон надежности.
Рисунок 3.12. Схема расчета показателей надежности (к задаче 3.7).
Ответ PC(100) = 0,99, TO.C = 1500 (ч)
Рисунок 3.13. Схема расчета показателей надежности (к задаче 3.8).
Ответ PC (100) = 0,994, TO.C = 1760 (ч)
вероятность отказа устройства за 100 часов работы, если справедлив
экспоненциальный закон надежности.
Ч 21 —Т
! ° ______ ” 22 f
Ч 21 —I
------------ 2з Ч
Рисунок 3.14. Схема расчета показателей надежности (к задаче 3.9).
Ответ PC (100) = 0,999, Q(100) = 0,001
Ответ PC (1) = 0,999
Рисунок 3.15. Схема расчета показателей надежности (к задаче 3.11).
Ответ Toc = 1382 ч.
равны TO.1 и TO.2 . Справедлив экспоненциальный закон надежности для
элементов. Необходимо найти среднюю наработку до первого отказа
системы.
Ответ TO.C
3TO .1 • TO .2 ч
2(T0.1 + TO .2)
Ответ PC(100) = 0,98
надежности. Определить вероятность безотказной работы системы за 500 ч
работы.
P=0,9
P=0,9
P=0,9
Рисунок 3.16. Схема расчета показателей надежности (к задаче 3.14).
Ответ PC (500) = 0,983
Ответ ToC = 100 ч.
Рисунок 3.17. Схема расчета показателей надежности (к задаче 3.16).
a) изделие нерезервированно и средние наработки до первого отказа
элементов равны To1 = To 2 = 300 ч;
Какой из вариантов более предпочтителен с точки зрения надежности,
если надежность изделия оценивать средней наработкой до первого отказа.
Ответ. Более предпочтителен вариант b, так как для варианта а: Toc = 150
ч, для варианта b: To с = 200 ч, для варианта с: То с = 180
надежности имеется возможность либо облегчить режимы работы элементов
и тем самым снизить интенсивность отказов изделия вдвое, либо
дублировать изделие при постоянно включенном резерве без облегчения
режимов работы элементов.
Какой способ более целесообразен, если надежность изделия оценивать
средней наработкой до первого отказа?
Ответ. Более целесообразно облегчить режимы работы элементов, так как в
этом случае среднее время безотказной работы изделия возрастет вдвое, а
при дублировании – только в 1,5 раза.
Ответ. Более целесообразно дублировать изделие, так как при дублировании
вероятность безотказной работы Рс (20) = 0,93, а при облегченном режиме
работы элементов Рс (20) = 0,85.
Ответ То.. = 20480 ч.
b) скользящее резервирование при одном резервном элементе, находящемся
в ненагруженном состоянии.
Какой из вариантов более целесообразен с точки зрения надежности, если
интенсивность отказов элемента Л ?
а)
Рисунок 3.18. Схема расчета показателей надежности (к задаче 3.20).
Ответ. Pa (t) = Рь (t) = e~2Л't (1 + 2Л • t), т.е. варианты равноценны.
имеет среднюю наработку до первого отказа, равную
То = 1 • Для
Л
повышения надежности применено скользящее резервирование, при котором
m резервных элементов находятся в ненагруженном режиме. Необходимо
найти среднюю наработку до первого отказа системы.
Ответ Tn r = —------.
O.C. N
Ответ n = 4.
1
резервирование с кратностью m = —.
Ответ Pc (t) = 0,919
«холодный» резерв. Определить вероятность безотказной работы системы за
один час времени работы с учетом резервирования.
Ответ PC (1) = 0,999
Ответ n = 5.
Ответ PC (2000) = 0,995
Рисунок 3.19. Схема расчета показателей надежности (к задаче 3.27).
Ответ PC(100) = 0,998
3.28. Найти вероятность безотказной работы системы, если вероятность
безотказной работы элемента p (t) = 0,3. Для элемента применено
1
резервирование с кратностью m =
Ответ PC (t) = 0,084
Рисунок 3.20 Схема расчета показателей надежности (к задаче 3.29).
Ответ PC (t) = (P + P, - P • P,) • P,
Ответ: а) P(100) = 0,91; б) P(100) = 0,992 ; в) P(100) = 0,997
вероятность безотказной работы в течение времени t = 10 ч и среднюю
наработку до отказа системы.
Ответ PC (10) = 0,99 TOC = 150 ч
отказав Л = 5 -10-3 1/ч. Кратность резервирования m = 3. Рассчитать
вероятность безотказной работы в течение времени t = 50 ч и среднюю
наработку до первого отказа системы.
Ответ PC (50) = 0,99 TO.C = 330 ч
Рисунок 3.21 Схема расчета показателей надежности (к задаче 3.33).
Ответ Лс (t) =
8Л - e ~4Лt (1 - e ~4Лt)
1 - (1 - e" 4 Лt )2 +
Л(10 - e ~4Лt)
о - 4 Лt
2 - e
с постоянно включенным резервом.
Ответ: а) Л <----C—; б) Л <----——--------; в)
nt nt
- ln[ - J-nPM J
t
Ответ: РП(100) = 0,992; То.П = 1400 ч
работы в течение времени t = 100 ч и среднюю наработку до отказа.
Ответ: РП(100) = 0,98; ТО.П = 1000
Ответ: Рб (t) = 33^ = 0,96; б) Рб = 1 - 1 -3/Ру (t) = 0,81;
в) Рб (t) = 3Л -1- - Ру (t) = 0,92
Ответ: РУ(40) = 0,98; ТОу = 208 ч
с нагруженным резервом и дробной кратностью «два из трех» обеспечило
повышение надежность вычислительного устройства.
Ответ РБ (t) > 0,5
Рисунок 3.22 Схема расчета показателей надежности (к задаче 3.42).
Ответ Pc (t) = 0,821
3.43. Определить вероятность безотказной работы системы, показанной
на рис. 3.23, если известны вероятности безотказной работы ее элементов
pi( t) = 0,9 p 2( t) = 0,95 p з( t) = 0,8 p 4( t) = 0,97 p 5( t) = 0,97.
Рисунок 3.23 Схема расчета показателей надежности (к задаче 3.43).
Ответ Pc (t) = 0,849
3.44. Устройство состоит из двух равнонадежных блоков соединенных
последовательно. Интенсивность отказов одного блока Л = 0,01 1/ч. Для
повышения надежности предлагается ввести дублирование с нагруженным
резервом. Определить какой вид резервирования оптимальный, если время
работы устройства t = 10 ч.
Ответ: раздельное резервирование лучше, т.к. для общего резервирования
Рс (10) = 0,967 To с = 75 ч; для раздельного резервирования Рс (10) = 0,982
To , с = 91,66 ч
3.45. Устройство состоит из четырех элементов. Интенсивность отказов
одного элемента Л = const = 0,01 1/ч. Для повышения надежности
предлагается ввести резервирование устройства с постоянно включенным
резервом и кратностью m = 4. Определить вероятность безотказной работы
системы в течение времени t = 50 ч и среднее время наработки до первого
отказа.
Ответ: PC(50) = 0,517; TO.C = 57 ч
Ответ: а) То У = 1500 ч; б) Т0У = 2000 ч
Ответ: а) То.. У = 150 ч; б) То.У = 200 ч
Ответ РБ (100) = 0,96
Требуется определить вероятность отказа блока в течение времени t = 250 ч
для следующих вариантов: а) резерв отсутствует; б) имеется один приемник в
ненагруженном скользящем резерве.
Ответ: а) QE (250) = 0,18; б) QБ (250) = 0,015
равна Ло = const = 10 3 1/ч. Определить вероятность безотказной работы
РР (t) в течение времени t = 100 ч и среднюю наработку до отказа
радиопередающего канала.
Ответ: РР(100) = 0,99; TO,Р = 3000 ч
Ответ: Ро( t) = 0,896
Ответ: Ро( t) = 0,999
Ответ: Р0(t) = 0,861
Ответ: Р0( t) = 0,945
голосования «1 из 2». Определить вероятность появления сигнала ложного
«0» системы.
Ответ: Po(t) = 0,826
Восстанавливаемое изделие это изделие, для которого в
рассматриваемой ситуации проведение восстановления работоспособного
состояния предусмотрено в нормативно-технической и (или)
конструкторской документации. Восстановление возможно с прекращением
выполнения изделием своих функций и без нарушения выполнения своих
функций.
При разработке сложной электронной аппаратуры, цифровой техники и
систем автоматического управления практически встречаются следующие
случаи восстанавливаемости изделий.
Основной круг задач, рассматриваемых при расчете надежности
восстанавливаемых систем, относится к следующей ситуации. Исправное
изделие начинает эксплуатироваться в момент времени t = 0 и, проработав
случайное время Xi , выходят из строя. На ремонт требуется случайное время
Yi . Этот процесс продолжается в течение всего срока службы изделия,
причем величины Xi и Yi (i = 1,2...) независимы.
Рассмотрим наиболее распространенный метод расчета надежности
восстанавливаемых систем.
МЕТОД, ОСНОВАННЫЙ НА ИСПОЛЬЗОВАНИИ КЛАССИЧЕСКОЙ
ТЕОРИИ ВЕРОЯТНОСТЕЙ
При решении большого класса задач удобно исходить из вероятностей
нахождения системы в том или ином состоянии. В общем случае число таких
состояний будет больше двух, но при решении задач теории надежности
обычно приходится иметь дело с конечным числом состояний.
Пусть в момент времени t система находится в состоянии i. Если
вероятность Pij (t, t + Аt) перехода системы за время Лt из состояния i в
состояние j не зависит от поведения системы до момента t, то такой
случайный процесс называется марковским процессом. Если эта вероятность
также не зависит от момента t , то имеет место однородный марковский
процесс.
Для этого случая можно найти характеристики надежности путем
решения дифференциальных уравнений.
Рассмотрим очень важный для теории надежности случай, когда
потоки отказов и восстановлений являются простейшими: А = const,
и = const. Это значит, что производительность труда ремонтника постоянна
и не зависит от времени. Поэтому время восстановления имеет
экспоненциальный закон распределения F (t) = 1 - e ~ ^; TB = —.
Например, пусть требуется найти коэффициент готовности KГ (t).
Если система исправна будем говорить, что она находится в состоянии (Gq),
если неисправна и восстанавливается - в состоянии (Q). Обозначим
вероятности нахождения системы в момент (t, t + А t) в этих состояниях через
Po( t), Pq (t + А t) и P1( t), P1(t + А t) соответственно. Естественно, что
Po(t) + P1(t) = 1 и KГ = Pq(t). Обозначим также через PQ1 (Лt) и P1Q(Лt)-
условную вероятность того, что в момент времени t система находится или в
состоянии Gq или в состоянии G1, а в момент времени t + Аt или в состоянии
G1 или в состоянии Gq , т.е. за интервал времени Аt произошел отказ
(восстановление) системы. Тогда
Будем считать, что за время At может произойти только один отказ или
только одно восстановление. Тогда на интервале At могут произойти четыре
несовместимых события: А1 (Go, G0) - в момент времени t система находилась
в состоянии G0 , в момент времени она осталась в том же состоянии,
т.е. отказа не произошло; A2 (Go, G1) - отказ произошел; А3 (G1, G0) -
восстановление произошло; А4 (G1, G1) - восстановление не произошло. Тогда
или
Po (t + A t)-Po (t'
A t
A t
A’, (t)+vP (t)
= APo (t)+vP (t).
Положим At ^ 0. Тогда получим систему дифференциальных
уравнений:
dPA = -APo (t)+vP (t)
dt
dP^ = APo(t)-vP(tt) ’
I dt
(4.1)
Система дополняется условием P (t)+ P (t) = 1.
Решение системы (4.1) при начальных условиях Po(t) = 1 и P1 (t )= 0 (в
начальный момент времени система работоспособна) имеет вид:
Po (t ) = K г (t) =
P (t )= K П (t) =
VA
—+e
A + v A + v
AA
---e
.-(a+v) t
,-(a+v )t
(4.2)
Если в начальный момент времени система неработоспособна, то
Po (o) = o, Pi (o) = 1 и решение системы имеет вид
v v e-(A+v )t
Pl (t ) = K П (t) =
A + v e-(A+v )t
(4.3).
При t ^ от независимо от начального состояния системы (Go или G1)
вероятности Ро( t) = K Г (t), P^t ) = K П (t) стремятся к постоянным значениям
K Г = Д ; K П = Д .
(4.4)
Это означает, что при экспоненциальных законах распределения
времени наработки на отказ и времени восстановления, случайный процесс
работы восстанавливаемой системы стабилизируется, и вероятность застать
систему работоспособной в произвольный момент времени остается
постоянной. Система с указанными свойствами называется эргодической, а
сам процесс – марковским случайным процессом.
При малых значениях величины (А + a)t, когда e-(А+A)t = 1 - (А + a)t,
из уравнения (4.2) получим:
к Г (t) = ~A+А [1 - (А + A) t ] = 1 - At;
K п (t) =
(4.5)
Отсюда следует, что в начальный период эксплуатации коэффициент
готовности примерно равен вероятности безотказной работы P(t), а
вероятность P1(t) - вероятности отказа Q(t) .
При нескольких работоспособных состояниях
n
к Г = Е jt). (4.6)
j'=1
где n - число работоспособных состояний; Pj (t) - вероятность j -го
работоспособного состояния.
Часто число неработоспособных состояний значительно меньше числа
работоспособных. При этом удобнее вычислить коэффициент простоя:
m+1-n
кП (t) = Е P( t), (4.7)
i=1
где Pi(t) - вероятность i-го неработоспособного состояния; m +1 - общее
число состояний.
Величина ошибки AK Г (t) при использовании коэффициента
готовности для оценки вероятности нахождения в работоспособном
состоянии вычислении при KГ > 0,9 составит:
KГ
AKГ(t) - (1 - KГ)KГ-KГ -100% (4.8)
МЕТОД, ОСНОВАННЫЙ НА ИСПОЛЬЗОВАНИИ ТЕОРИИ ГРАФОВ
Метод состоит в следующем:
Рассмотрим эту методику на частном примере.
Пусть дано некоторое устройство, для повышения надежности
которого применено резервирование. Известны:
Предполагается, что при замене работающего элемента на резервный
перерыва в работе устройства не происходит, поэтому отказ устройства
наступает при одновременной неработоспособности основного и всех
резервных элементов (состояние m +1).
Необходимо вычислить вероятность безотказной работы P(t) в течение
времени t и вероятность того, что резервированное устройство будет
исправно в любой момент времени t ( KГ (t)).
Решение. Сделаем следующие допущения:
Рассмотрим случай ненагруженного резерва с абсолютно надежным
переключателем и одной ремонтной бригадой, обслуживающей устройство
(ограниченное восстановление)
При указанных выше допущениях функционирование
резервированного устройства можно представить графом, изображенным на
рис. 4.1.
λ λλλ
G о.. . . . ч Gj L k . . .4Gm+1
µ µµ
λλλλ
Рисунок 4.1. Граф состояний резервированной восстанавливаемой
системы в ненагруженном резерве при ограниченном а) и неограниченном б)
восстановлении
Узлам графа соответствуют состояния системы (0, 1, 2, …, m+1), а
ветвям – возможные переходы из одного состояния в другое.
Искомая система дифференциальных уравнений может быть
составлена с помощью графа по следующим правилам. Производная
вероятности состояния равна сумме стольких слагаемых, сколько стрелок
связано с этим состоянием. Каждое слагаемое равно произведению
интенсивности потока событий, переводящего систему по данной стрелке, на
вероятность того состояния, из которого исходит стрелка. Слагаемое имеет
знак минус, если стрелка исходит из данного состояния, а знак плюс – если
стрелка направлена в данное состояние. Полученная система уравнений
называется системой уравнений Колмогорова.
Для графа состояний (рис.4.1,а) получим следующую систему
дифференциальных уравнений.
Pt)=-AP0( t)+ri1( t)
dt
....
dPi (t)
~ri = APj—1 (t) - (A + a)Pj (t) + ^Pj+i (t)
t
(4.9),
. . . .
driri = APm (‘) - riPm+1( t)
dt
где J = 1,...,m.
Надежность восстанавливаемого изделия, как правило, определяется
при условии, что в момент включения все элементы исправны. Тогда
Pj (0) =
1 j = 0,
0, j * 0.
Начальные условия для рассматриваемого графа состояний: /0(0) = 1;
P1(0) = ... = Pj (0) = Pm+1(0) = 0 .
Система решается с помощью преобразований Лапласа или
dPj (t)
численными методами. При t ^ да производные —---> 0 и система (4.9)
dt
превращается в однородную систему линейных уравнений:
'0 = -APg( t) + riP1( t)
....
0 = APj-1 (t) - (A + riPj (t) + ^Pj+1 (t)
....
0 = APm (t) - ri^m+1( t)
Для решения системы (4.10) необходимо добавить уравнение
m+1
Z PJ =1
(4.10)
(4.11)
J=9
В результате решения системы (4.10) совместно с уравнением (4.11)
получим установившиеся значения коэффициентов простоя и готовности:
1
K П Pm+1
А
-
K Г = 1 - Pm+1 = 1
1
.
(4.12)
А
—
Если та же система обслуживается (m +1) ремонтными бригадами
(неограниченное восстановление), то очередь на ремонт отсутствует. Граф
состояний для ненагруженного резерва и неограниченного восстановления
представлен на рис.4.1,б. В результате решения системы уравнений
совместно с нормировочным уравнением (4.11) получим
K П = Pm+1 =
1
m+'(m +1)! Г а
-^
j=Q j!
\ m+1-j
K Г = 1 — Pm+1 = 1
1
(4.13).
m+ 1( m +1)! ( а
-^
j=Q j!
\ m+1-j
Графы состояний для системы, состоящей из одного основного и m
элементов в нагруженном резерве, представлены на рис.4.2.,а для
ограниченного восстановления и б для неограниченного.
m+1)2
G0
А
Gj
А ^^^А
a)
G0
.^‘
(m+2-/)^rA(m+1—/_Л
Gj
Л
Gm+.
(m +1)а
б)
Рисунок 4.2. Граф состояний резервированной восстанавливаемой
системы в нагруженном резерве при ограниченном а) и неограниченном б)
восстановлении
Составляя систему уравнений и решая ее совместно с нормировочным
уравнением (4.11), получим
для ограниченного восстановления
K П Pm+ТГ
K Г = 1 - Pm+1 = 1
1
m+11 Л V
j=0 j! к л у
(4.14)
для неограниченного восстановления
K г = 1 — Pm+1 = 1
Л
K П = Pm+1 =
\ m+1
\ m+1
m
j=0
m +1- j
m +1
\m+1-j (
(4.15)
Выражение (4.15) представляет собой вероятность случайного
исхода,
имеющего биноминальное распределение. Это объясняется независимостью
отказов и восстановлений элементов.
Если изделие состоит из одного элемента и m резервных элементов в
нагруженном резерве с разными интенсивностями отказов, то число
возможных состояний равно 2m+1 и граф состояний имеет разветвленный
вид. На рис.4.3 представлен граф состояний для двух элементов с
различными интенсивностями отказов и неограниченным восстановлением.
G
G
G
G
1
Рисунок 4.3. Граф состояний резервированной восстанавливаемой
системы, состоящей из двух элементов с различными интенсивностями
отказов
В этом случае
Kг(t) = 1 - Д(t); Kп(t) = Рз(t).
При t >/ KГ = 1 — P3; KП = Р3.
Рассмотрим резервированные системы, для которых отказы
недопустимы, но ремонт отказавшего элемента производится во время
выполнения задачи. Если система состоит из основного элемента и m
элементов в нагруженном резерве, то для случая ограниченного
восстановления граф состояний системы представлен на рис.4.4.
(m+1)2 (m+2-Mm+j 2z^. Z^
G-G^T LI
Рисунок 4.4. Граф состояний резервированной восстанавливаемой
системы, состоящей из одного и m элементов в нагруженном резерве.
Состояние m +1 является поглощающим
При попадании системы в состояние m +1 происходит отказ системы,
который недопустим и приводит к невыполнению поставленной задачи.
Вероятность безотказной работы системы при условии, что в
начальный момент времени t = 0 в системе нет неисправных элементов, т.е.
P0(0) = 1; Д(0) =... = Pm+1(0) = 0 равна
m
P'(ti) = Z Pj( ti). (4.16)
i=0
Вероятность отказа системы в течение времени выполнения задачи
также является условной вероятностью и соответственно равна
Q'(ti) = Pm+1( ti). (4.17)
Важным показателем является средняя наработка до отказа
^ ^ m
TO = J P‘(t)dt = J £ pj (t)dt. (4.18)
0 0 i=0
Примеры резервированных систем, отказы которых недопустимы во
время функционирования: система управления выведением на орбиту
космического аппарата, радиолокационная станция, обеспечивающая
посадку самолета, систему автоматического управления тормозами тепловоза
и т.п. Конструктивное исполнение таких систем позволяет производить
ремонт непосредственно во время применения.
При решении системы уравнений, составленных по графу рис.4.4 с
помощью преобразований Лапласа, целесообразно использовать правило
облегчающее расчет.
Для определения средней наработки до отказа достаточно найти
преобразование Лапласа вероятности безотказной работы P (s) и подставить
в него s = 0.
МЕТОД ПРЕОБРАЗОВАННЫХ СЕТЕЙ
Основная идея этого метода заключается в том, что сети придается
такая конфигурация, для которой можно применять известные расчетные
формулы. При этом исходная сеть преобразуется дважды: в сеть с заведомо
более высокой и в сеть с заведомо более низкой надежностью. В результате
расчетов показателей надежности двух полученных сетей получают
верхнюю и нижнюю границы интервала показателей надежности, внутри
которого находится значение показателя надежности рассматриваемой сети.
Преобразование исходной сети заключается в поэтапном исключении
из рассмотрения отдельных линий путем их закорачивания или разрыва. При
этом разрыв любой линии связи уменьшает надежность сети или переводит
ее в неработоспособное состояние, а закорачивание линии связи
(равносильно объединению двух узлов) повышает надежность сети.
В ходе преобразования сети намечают начальный и конечный узлы,
поочередно выбирают линии связи и преобразуют сеть до тех пор, пока не
будут получены последовательно-параллельные цепи. Правил выбора тех
или иных линий связи не существует.
Рассмотрим применение этого метода на примере сети, показанной на
рис.4.5. Для простоты будем считать узлы сети абсолютно надежными, а
линии связи равнонадежными.
б)
в)
18 18
Рис.4.5 Преобразование сети:
г)
ж)
а) - исходная схема; б) - первый этап преобразования путем удаления ветви 5;
в) - второй этап преобразования путем удаления ветви 4; г) - эквивалентная
логическая схема для расчета нижней границы показателей надежности; д) -
первый этап преобразования путем закорачивания ветви 5: е) – второй этап
преобразования путем закорачивания ветви 4; ж) – эквивалентная схема для
расчета верхней границы показателей надежности.
Для получения сети с заведомо более низкой надежностью необходимо
последовательно разорвать линии 5 и 4 (рис.4.5,б, в). В результате этих
преобразований получим эквивалентную логическую схему для расчета
показателей надежности (рис.4.5, г)
Для получения сети с заведомо более высокой надежностью
необходимо последовательно закоротить линии 5 и 4 (рис.4.5,д, е). В
результате этих преобразований получим вторую эквивалентную логическую
схему для расчета показателей надежности.
В рассматриваемом примере использование метода марковских цепей
для расчета показателей надежности приведет к сложным решениям из-за
большого числа состояний системы, поэтому используют простой
приближенный метод расчета, согласно которому для последовательного
включения подсистем:
n
Ka = 1 - n + E К Ai; (4.19)
i=1
для параллельного включения:
m
Ka = 1-П (1 — KAi), (4.20)
i=1
где KÃ - коэффициент готовности последовательной (параллельной) группы
из n(m) подсистем.
Если в системе применяется скользящее резервирование подсистем, то
для расчета коэффициента готовности вместо формулы (4.20) применяется
формула:
m
Ka = Z cmKAC (1 — KAi) m-i (4.21)
i=k
В табл.4.1 приведены нормированные значения показателя
коэффициент готовности KГ телекоммуникационных сетей.
Таблица 4.1. Нормированные значения показателя коэффициент готовности
Тип сети электросвязи |
Норма |
Сеть междугородной и международной телефонной |
Не менее 0,999 |
Сеть зоновой телефонной связи |
Не менее 0,9995 |
Сеть местной телефонной связи |
Не менее 0,9999 |
Телеграфная сеть связи и сеть Телекс |
Не менее 0,9999 |
Сеть передачи данных |
Не менее 0,99 |
4.2. Типовые примеры и их решения
Пример 4.1. Коэффициент простоя K П = 0,1. Интенсивность
восстановления ^ = 0,1 -10—3. Определить вероятность отказа системы за 100
ч работы, если справедлив экспоненциальный закон надежности.
Решение. Для решения данной задачи воспользуемся следующими
соотношениями:
кп = 1 - к г ;
K
Г =
;
Q (t) = 1 - e- Л
Определим интенсивность отказов системы
л= д(1 - К г ) = дК п
К г 1 - КП
0,1 -10-3 - 0,1
1 - 0,1
= 0,1-10-4
(1/ч)
Определим вероятность отказа системы за время t = 100 ч работы
Q (100) = 1 - e - 0,1-10 —4-100 = 0,001
Пример 4.2. Аппаратура имела среднюю интенсивность отказа
Л = const = 0,001 1/ч и среднее время восстановления ТВ = 1 ч. Определить
коэффициент готовности KГ и вероятность отказа аппаратуры за 1 ч работы.
Справедлив основной закон надежности.
Решение. Для решения данной задачи воспользуемся следующими
соотношениями:
к - Д ■
1
Q (t) = 1 - e-Л
.
Определим коэффициент готовности аппаратуры
Кг = -^- = —=--=--1--= 0,999
Г Л + д Л-TB +1 0,001 +1
Определим вероятность отказа аппаратуры за время t = 1 ч
Q (t) = 1 - e-Л = 1 - e-0,001 = 0,001
Пример 4.3. Коэффициент готовности изделия KГ = 0,9. Среднее
время восстановления ТВ = 100 ч. Найти вероятность безотказной работы
устройства за 10 ч, если справедлив экспоненциальный закон надежности.
Решение. Для решения данной задачи воспользуемся следующими
соотношениями:
K
Д .
1
P (t) = e ~Л
.
Определим интенсивность отказов системы
, М(1 — K Г) 1 — K Г 1 — 0,9
KГ
TB • KГ 100 ■ 0,9
= 0,001 (1/ч)
P (10) = eЛ0 = e "0,01 = 0.99
Пример 4.4. Известно, что интенсивность отказов устройства Л = 0,02
1,ч, а среднее время восстановления ТВ = 10 ч. Требуется вычислить
коэффициент готовности KГ и функцию готовности PГ (t) устройства.
Решение. Для
решения задачи воспользуемся следующими
соотношениями:
1
TO = Л ’
кг = - TO- ;
Г То + Тв
t
—
K Г (t) = K Г + (1 — K Г) e K Г' ТВ .
Определим среднюю наработку до первого отказа устройства
1
1
0,02
= 50 ч
Тогда коэффициент готовности устройства К
TO
То + Тв
50
50 +10
= 0,83.
Функция готовности
tt
к Г (t) = к Г + (I — K Г) e K Г ■ ТВ = 0,83 + (1 — 0,63) e °,8340 = 0,83 + 0,17e "0,12 t
Пример 4.5. В случае отказа непрерывно работающего устройства оно
некоторое время находится в нерабочем состоянии – ремонтируется. После
ремонта устройство немедленно начинает работать. Интенсивность отказа
устройства Л = 0,001 1/ч, а интенсивность восстановления м = 0,1 1/ч.
Определить вероятности появления возможных состояний станции в
моменты времени t1 = 5 и 12 = 30 ч. Рассмотреть два случая: а) устройство в
момент начала функционирования (t = 0) работоспособно; в) устройство в
момент начала функционирования (t = 0) неработоспособно.
Решение. Устройство может находиться в двух состояниях:
работоспособном и неработоспособном. Для определения вероятностей этих
состояний воспользуемся выражениями (4.2) в случае, если устройство в
момент времени t = 0 работоспособно и (4.3) - если неработоспособно. Тогда
для случая а)
Р0 (t ) = K Г (t) = -Д + |
Л _ e Л + д |
- (Л+д ) t = 0,1 + 0,001 + 0,1 |
0,001 - e 0,001 + 0,1 |
- 0,505 = 0,996 |
Л Р1 (t ) = K П (t) =-- |
Л e |
-лД)t _ 0,001 |
0,001 |
- 0,505 = 0,04 |
e | ||||
Л + д |
Л + д |
0,001 + 0,1 |
0,001 + 0,1 | |
Р0 (t ) = K Г (t) = -^- + |
Л _ e |
- (Л+Д)t = 0,1 |
0,001 - e |
■3,03 = 0,90 5 |
Л + д |
Л + д |
0,001 + 0,1 |
0,001 + 0,1 | |
Р1 (t ) = K П (t) = --- |
Л e |
-(л+Д)t _ 0,001 |
0,001 |
-3,03 = 0,95 |
e | ||||
Л + д |
Л + д |
0,001 + 0,1 |
0,001 + 0,1 | |
для случая в) | ||||
Ро (t ) = Kr (t ) = -^- - 0Г |
Д- e - |
(Л+Д)t _ 0,1 |
0,1 „- |
0,505 = 0,9 8 6 |
e | ||||
Л + д |
Л + д |
0,001 + 0,1 |
0,001 + 0,1 | |
Р1 (t ) = K П (t) =-----+ |
Д e" |
(Л+Д )t= 0,001 |
0,1 - e |
'0,505 = 0,014 |
Л + д |
Л + д |
0,001 + 0,1 |
0,001 + 0,1 |
/>0 (t ) = K Г (t ) =
Д Д е~Л+^ )t
Л + д Л + д
Р1 (t) = K П (t) = ^— + -Д e -(Л+Д )t
Л + д Л + д
----0,1---0,1----e ~3,03 = 0,99
0,001 + 0,1 0,001 + 0,1
0,001 0,1
-----------+-----------e 3,03 = 0,01
0,001 + 0,1 0,001 + 0,1
Пример 4.6. В случае отказа непрерывно работающего устройства оно
некоторое время находится в нерабочем состоянии – ремонтируется. После
ремонта устройство немедленно начинает работать. Интенсивность отказа
устройства Л = 0,001 1/ч, а интенсивность восстановления д = 0,1 1/ч.
Определить среднее значение коэффициента готовности в течение наработок
t = 5 ч и t = 30 ч для случая, когда в момент начала функционирования
устройства работоспособно.
Решение. Среднее значение функции готовности вычисляется по
формуле:
Подставив выражение
ti
K Г (t) — - f K Г (t) dt.
ti 0
для K Г (t) согласно (4.2) и выполнив
интегрирование получим:
ti
K г (t) — - JK г (t) dt — -
ti
1
ti
0
t
; i 01Л + д
+ Л e ’(Л+д) ti 1 dt —
(1 — e—( Л + д ) ti )
Для заданных числовых значений t KГ(5) - 0,991 KГ(30) - 0,99.
Пример 4.7. Определить величину ошибки AKГ при вычислении
вероятности нахождения восстанавливаемой системы в работоспособном
состоянии в момент времени, равный наработке до отказа То, если K Г — 0,6,
а распределении времени до отказа и времени восстановления подчиняется
экспоненциальному закону. В начальный момент времени t — 0 система
работоспособна.
Решение. При экспоненциальном законе распределения времени до
отказа и времени восстановления коэффициент готовности можно найти по
формуле
K Г -
д
TO
Подставив в выражение для функции готовности – вероятности
нахождения системы в работоспособном состоянии t — То найдем
соотношение между функцией готовности и коэффициентом готовности.
1
— I
—
(Л+д > — K Г + (1 — K Г) e 1—K Г
Ошибка при использовании коэффициента готовности для оценки
вероятности нахождения в работоспособном состоянии в момент времени
t — То
1
AK Г — K Г (t) — K Г — (1 — K Г) e 1—K Г — 0,4 • 0,082 — 0,033
Пример 4.8. Поездная радиостанция включает в себя приемный и
передающий блоки, интенсивность отказов которых одинаковы и равны
к = 10 2 1/ч. Интенсивность восстановления ц = 2 1/ч. Станцию обслуживает
одна ремонтная бригада. При неработоспособности любого из блоков
радиостанция неработоспособна. При этом работоспособный блок
выключается, и в нем не могут происходить отказы. Определить значения
коэффициентов готовности и простоя радиостанции.
Решение. Число возможных состояний поездной радиостанции равно
двум: Gо - работоспособное состояние (оба блока работоспособны); G1 -
неработоспособное (один из блоков неработоспособен). Граф состояний
представлен на рис.4.6.
22
Ц
Рисунок 4.6. Граф состояний (к примеру 4.8)
По графу составляем систему дифференциальных уравнений:
<
dPo (t)
dt
dW
dt
= -2кРо (t) + p₽i (t)
= 2кРо (t) -pP (t)
При t > v перейдем к системе алгебраических уравнений:
- 2кР0 (t) + pPi (t) = 0
1 2кРо (t) -рР (t) = 0
Заменив одно уравнение нормировочным условием: Р0 + Pi = 1 имеем:
- 2кР0 (t) + pPi (t) = 0
1 Р0 + Pi = 1
В результате решения этой системы уравнений получим:
2
2 + 0,02
= 0,99; K п = Pi =
0,02
2 + 0,02
= 0,0i.
Пример 4.9. Восстанавливаемое устройство, время безотказной работы
и время восстановления которого подчиняются экспоненциальному закону
распределения, имеет коэффициент готовности KГ = 0,9. Определить
вероятность нахождения устройства в работоспособном состоянии KГ (t) в
момент времени ti = 50 ч, если наработка до первого отказа То = 500 ч.
Решение. Зная коэффициент готовности и среднее время наработки до
первого отказа, найдем среднее время восстановления, используя формулу
(1.23)
Kr = _ TO- = 500_ = 0,9
Г To + Tb 500 + Tb
500 - 450
---------= 56 ч.
0,9
Используя выражение (4.2) найдем соотношение между функцией
готовности и коэффициентом готовности устройства.
_ _ (TO + TB \t
-t
Kr(t) = _ 0_ + _ B_ e ^TO'TB ^ = Kr + (1 -Kr)e KГ'TB
Г To + TB To + TB ГГ
Для t = 50 ч
50
KГ(t) = 0,9 + (1 - 0,9)e °,9'56 = 0,9 + 0,1 • 0,37 = 0,94.
Пример.4.10. Для питания электростанции используется
электроагрегат с двумя генераторами, каждый из которых обладает
производительностью, достаточной для нормальной работы. Эти генераторы
работают поочередно. При отказе работающего генератора в работу
включается резервный генератор, а отказавший отключается и
ремонтируется. Отказ электроагрегата состоит в прекращении питания
радиостанции.
Конструкции электроагрегата допускает одновременный ремонт обоих
генераторов, имеется нужное число ремонтников.
Определить коэффициент готовности электроагрегата, если наработка
на отказ генератора TO в 5 раз больше среднего времени устранения отказа
TB . Предполагается, что время безотказной работы и время восстановления
подчиняются экспоненциальному закону распределения.
Решение. Электроагрегат может находиться в одном из трех состояний:
Обозначим вероятности указанных состояний в момент времени t
через Po (t), P1 (t), P2 (t). Эти вероятности при t > ® имеют пределы Po, P,
P2. Поскольку для рассматриваемого агрегата переход из состояния Gо в
состояние G1 не нарушает его работоспособности K Г = Po + Pi.
Составим граф состояний (рис.4.7) и соответствующую этому графу
систему дифференциальных уравнений.
AA
G0 G1G2
А2a
Рисунок 4.7. Граф состояний (к примеру 4.10)
' dPo( t)
dt
dPi( t)
dt
dP2( t)
. dt
= -APo( t) + APi( t)
= APo(t) - (A + a)P^t) + 2^P2(t)
= AP1( t) - 2^( t)
Для определения установившихся значений Po и Pi положим все
производные равными нулю. Получим систему линейных уравнений, которая
добавляется характеристическим уравнением Po + Pi + P2 = 1:
0 = -APo( t) + APi( t)
0 = APi( t) - 2aP2( t)
В результате решения системы уравнений получим:
K = 2А2 , 2Aa
Г 2a 2 + 2Aa + A2 2a 2 + 2Aa + A2
При A = 0,2a получим KГ = 0,98.
Пример 4.11. Преобразователь состоит из рабочего блока и блока в
ненагруженном резерве. Распределения времен между отказами и
восстановлениями подчиняются экспоненциальному закону с параметрами
A = 8-10
3
1/ч, a = 0,8 1/ч. Определить значения коэффициента простоя и во
сколько раз уменьшается величина коэффициента простоя преобразователя
при применении неограниченного восстановления по сравнению с
ограниченным.
Решение. Для определения коэффициентов простоя для случаев
ограниченного и неограниченного восстановления воспользуемся
соответственно выражениями (4.12) и (4.13). Число возможных состояний
равно трем.
Для ограниченного восстановления
K
П
= р , =
1 m+1
1
1
2
Для неограниченного восстановления
K
П
1 m+1
1
m+1( m +1)! (1
Е^
j=о j!
m+1- j
1
2
Для рассматриваемого примера справедливо соотношение 1 >> Л и,
полученные выражения, могут с достаточной для практики точностью,
определены приближенно:
22
KПО » — *10-4; Kпн * * 0,5 -10-4.
12 212
При применении неограниченного восстановления по сравнению с
ограниченным величина коэффициента простоя уменьшилась в два раза.
Пример 4.12. Преобразователь состоит из рабочего блока и блока в
нагруженном резерве. Распределения времен между отказами и
восстановлениями подчиняются экспоненциальному закону с параметрами
Л = 8 -10-3 1/ч, 1 = 0,8 1/ч. Определить значения коэффициента простоя и во
сколько раз уменьшается величина коэффициента простоя преобразователя
при применении неограниченного восстановления по сравнению с
ограниченным.
Решение. Для определения коэффициентов простоя для случаев
ограниченного и неограниченного восстановления воспользуемся
соответственно выражениями (4.14) и (4.15). Число возможных состояний
равно трем.
Для ограниченного восстановления
K П = Pm+1 =
1
m+11 A jjV
^ 7 А
j=о j! к Л У
1
2
2 Л2
А 2 + 2 Ла + 2 Л2 .
Для неограниченного восстановления
K П = Pm+1 =
\ m+1
2
.
При а >> Л можно записать
2 Л2
K п . о * г * 2'10 ; K п . н
^2 * 10-4
А
Выигрыш при применении неограниченного восстановления также
получился равным двум.
Пример 4.13. радиоприемное устройство, состоящее из рабочего блока
и блока в нагруженном резерве, рассчитано на непрерывную круглосуточную
работу. Через три часа после включения это устройство может получить
команду на перестройку режима работы. Интенсивность отказов и
восстановлений каждого блока равны Л = 8 • 10=3 1/ч, а = 0,2 1/ч. Имеются
две дежурные ремонтные бригады. Определить вероятность застать
радиоприемное устройство в неработоспособном состоянии через три часа
работы (функцию простоя) и значение коэффициента простоя.
Решение. Радиоприемное устройство в любой момент времени может
находиться в одном из трех состояний:
G0 - оба блока работоспособны;
G1 - один блок неработоспособен;
G2 - оба блока неработоспособны.
При нахождении в состояниях G0 и G1 устройство работоспособно, в
состоянии G2
-
неработоспособно. Граф состояний устройства с
соответствующими интенсивностями переходов показан на рис.4.8.
Л
G0 G1 G2
А 2 а
Рисунок 4.8. Граф состояний (к примеру 4.13)
Система дифференциальных уравнений, составленная по этому графу,
имеет вид:
dPi^ = -2ЛPo( t) + цр( t)
dt
■ d^ = 2P (t) - (Л + ц) P (t) + 2рРг (t)
dt
dP^ = Лр( t) - 2 Ptt)
dt
Для определения функции простоя устройства решим эту систему при
начальных условиях Po (0) = 1; P (0) = P2 (0) = 0. Переходя к изображениям,
получим систему алгебраических уравнений:
'(s + 2 Л) P0° (s) - X (s) = 1
Для получения величин Pi (s) используем правило Крамера:
А-
P( s) = ^ ’
где А - определитель, элементами которого являются коэффициенты при
О
О
О
P0 (s), P (s), P2 (s); Ai - определитель, который образуется из А путем
замены i -го столбца коэффициентами правой части системы.
В примере требуется определить функцию простоя P2(t), поэтому
найдем определители А и А 2:
s + 2Л | - ц | 0 | s + 2Л | - ц | 1 | ||
А = | - 2Л | s + Л + ц | -2Л | А 2 = | - 2Л | s + Л + ц | 0 |
0 | -Л | s + 2ц | 0 | -Л | 0 |
Следовательно,
P," (s ) = А =---------------2 Л------------у .
Переходя от изображений к оригиналам, получим:
; 2
K (t) = P2 (t) =--------- [1 - 2e"(л+ц)t + e"2(л+ц)t ].
Используя это выражение, определим коэффициент простоя при t ^ ®
K
П =
(Л + a)2.
Подставив числовые значения, получим:
K П (3) = 2 -10-4, K П = 1,5 -10-3.
Пример.4.14. Поездная радиостанция включает в себя приемный и
передающий блоки, интенсивность отказов которых одинаковы и равны
Х = 10-2 1/ч. Интенсивность восстановления р = 2 1/ч. Станцию обслуживает
одна ремонтная бригада. При неработоспособности любого из блоков
радиостанция неработоспособна. При этом работоспособный блок не
выключается, и в нем могут происходить отказы. Определить значения
коэффициентов готовности и простоя радиостанции.
Решение. Число возможных состояний поездной радиостанции равно
трем: G0 - работоспособное состояние (оба блока работоспособны); G1 -
один из блоков неработоспособен; G2 - оба блока неработоспособны. Граф
состояний представлен на рис.4.9.
2ЛЛ
G0 G1
А^
Рисунок 4.9. Граф состояний (к примеру 4.14)
Этому графу соответствует система дифференциальных уравнений:
dP^ = -2ЛР0( t) + цР1( t)
dt
• ^ = 2 ЛР0 (t) - (Л + a)Pi (t) + ^ (t)
dt
^ = Лр( t) - р t)
dt
При t ^ да переходим к системе алгебраических уравнений:
- 2ЛР0( t) + Ai( t) = 0
<2ЛР0(t) - (Л + a)Pi(t) + Рt) = 0
ЛР1( t) - аР2( t) = 0
При решении этой системы используем нормировочное уравнение:
Р0 + Р1 + Р, = 1,
которое может заменить любое из уравнений системы. В результате решения
системы уравнений либо подстановкой, либо по правилу Крамера получим:
Po = —; P = —; P = .
Коэффициент готовности поездной радиостанции равен:
2
------------------= 0,99.
4 + 0,04 + 0,0001
Коэффициент простоя
Kn = P + P = ,2—— + 2—2 ,
0,04 + 0,0002
4 + 0,04 + 0,0001
= 0,01.
Пример.4.15. Определить коэффициент готовности KÃ и вероятность
безотказной работы за время t = 100 ч сети (рис.4.5,а), если интенсивность
отказа линий связи — = 0,2 •Ю"4 1/ч, а среднее время восстановления
T- = 1000 ч. Узлы сети считать абсолютно надежными, а линии связи
равнонадежными. Для устройств справедлив экспоненциальный закон
надежности.
Решение. Для решения задачи воспользуемся эквивалентными
логическими схемами для расчета показателей надежности, приведенными на
рис.4.5,г, ж. Схема (рис.4.5,г) соответствует общему резервированию с
постоянно включенным резервом. Минимальную вероятность безотказной
работы сети за время t можно определить по формуле:
Г -2\3
P^ (t) = 1 - (1 - e-2—• t)3 = 1 -11 - e-0,4'10 I = 0,999999936
Используя формулы (4.19) и (4.20), определим минимальный
коэффициент готовности сети
n
KAJ = 1 - [1 - (1 - n + £ K,r,
i=1
где KÃi - коэффициент готовности отдельной линии связи.
По данным об интенсивности отказов и времени восстановления линии
связи определим коэффициент готовности отдельной линии связи KÃi
1
Подставим полученное выражение в формулу для KÃ.Í получим
KaJ = 1 - [1 — (1 — n + £кЛ)]3 = 1 —
i=1
( n
1
3
= 0,99994
Схема (рис.4.5,ж) соответствует раздельному нагруженному
резервированию. Максимальную вероятность безотказной работы сети
определим по формуле
Рл (t) = |1 - (1 - e- Л )3 ]2
2
1 -[1 - e -°'2'10-2) = 0,999999984
Максимальный коэффициент готовности сети определим по формуле:
1 - 1
1
3
= 0,999985
Таким образом, для исходной сети вероятность безотказной работы
находится внутри интервала Pj (t)...Рл(t), а коэффициент готовности в
интервале Кл./ ...Кл.л.
Ответ Q(500) = 0,012
Ответ KГ = 0,994; Q(10) = 0,03
Ответ KП = 0,048; Q(1) = 0,049
Ответ P(100) = 0,591
Ответ P(10) = 0,9
Ответ т = 0,105 1/ч
Ответ: KП = 0,058; TO = 37,7 ч
Ответ: KГ = 0,83; KОГ = 0,83е=0,02t
Ответ: KГ(10) = 0,867
Необходимо определить вероятность застать систему в исправном состоянии
в момент времени t = 12 ч.
Ответ; KГ (12) = 0,987
Ответ: KГ(20) = 0,85
Ответ: KГ (t) = 0,91 + 0,09e“0,11t
Ответ: Kп = 2,5 • 10-3
Ответ: Kп = 8 • 10-4
ненагруженном резерве. Ремонт производит одна бригада. Среднее время
ремонта TB = 0,5 ч. Определить коэффициент простоя усилителя с резервом.
Ответ: Kп = 9 40 6
независимо друг от друга. Определить функцию готовности.
2
Ответ: Kr(t) = 1--у [1 - 2e"(Л+“)t + e"2(л+“)t
Г Л + А)2
интенсивность отказов устройства, чтобы удовлетворялось условие
KП < 2-10-4.
Ответ: Л << uJ —П = 2 -10-2 1/ч
2
Ответ: Л << аКП = 3 -10-4 1/ч
Ответ: Kп = 2 -10-4
Ответ: Kп = 10 4
резервом. Ремонт производится одной ремонтной бригадой с
интенсивностью р = 2 1/ч. Определить значение коэффициента простоя.
Ответ: Kп = 10-4
Ответ: Kп = 0,5 -10-4
Ответ: То = ^^^ = 2,6 -104 ч
O 2X2
до отказа системы.
1 (
Ответ: TO = —
= 51000 ч
Ответ: K Г = P0 + P1 + P2 = 0,99 + 0,99 -10-2 + 0,495 -10-4 = 0,99995
Ответ KГ = 0,99995
Ответ KГ = 0,9999
Ответ KГ = 0,899
Ответ KГ = 0,999
To .1, ц
TO.2, ц
TO .1, ц
TO.2, ц
Рисунок 4.10. Схема расчета надежности (к задаче 4.30)
Ответ Кг = 0,438
Для элементов справедлив экспоненциальный закон надежности.
Рисунок 4.11. Схема расчета показателей надежности (к задаче 4.31)
Ответ Кг = 0,9875
Рисунок 4.12. Схема расчета показателей надежности (к задаче 4.32)
Ответ Кг = 0,7426
4.33. Схема расчета показателей надежности восстанавливаемой
системы приведена на рис.4.13. Необходимо найти коэффициент готовности
системы, если известны интенсивности отказов и среднее время
восстановления элементов. Для элементов справедлив экспоненциальный
закон надежности.
Рисунок 4.13. Схема расчета показателей надежности (к задаче 4.33)
Ответ: Кг =
1 , 3А вв 2 -1
АТ.1 + 1 (А Тв 2 + 1)3
-1
системы, если известны среднее время наработки до отказа и среднее время
восстановления элементов. Для элементов справедлив экспоненциальный
закон надежности.
Рисунок 4.14. Схема расчета показателей надежности (к задаче 4.34)
т ■ т
1O2 Вв2
Т„+ Т„3 Т„.
Ответ: Kr = -OL———/1 +
Г (То, + Тв 1)4
(To 2 + Тв 2)2
+ Т--2
То з + Тв з
закон надежности.
Рисунок 4.15. Схема расчета показателей надежности (к задаче 4.35)
Ответ: 0,999913 < KГ < 0,999985
4.36. Определить коэффициент готовности K сети (рис.4.15), если
среднее время наработки до первого отказа линии связи То = 500 ч, а
интенсивность восстановления ^ = 0,1 -10-2 1/ч. Узды сети считать
абсолютно надежными, а линии связи равнонадежными. Для устройств
справедлив экспоненциальный закон надежности.
Ответ: 0,134827 < KГ < 0,279923
4.37. Определить коэффициент готовности K сети (рис.4.16), если
среднее время наработки до первого отказа линии связи То = 1000 ч, а
среднее время восстановления Тв = 1 ч. Узды сети считать абсолютно
надежными, а линии связи равнонадежными. Для устройств справедлив
экспоненциальный закон надежности.
Рисунок 4.16. Схема расчета показателей надежности (к задаче 4.37)
Ответ: 0,99999999979 < Кг < 0,99999999984
Ответ: 0,9999996385 < KГ < 0,999999971
Рисунок 4.17. Схема расчета показателей надежности (к задаче 4.39)
Ответ: 0,9999999991235 < Кг < 0,9999999993864
4.40. Определить коэффициент готовности KГ сети (рис.4.17), если
интенсивность отказа линии связи А = 0,7 -10-3 1/ч, а среднее время
восстановления ^ = 0,5 1|ч. Узды сети считать абсолютно надежными, а
линии связи равнонадежными. Для устройств справедлив экспоненциальный
закон надежности.
Ответ: 0,9999999509 < KГ < 0,9999999972
ЛИТЕРАТУРА
1.
Вентцель Е.С. Теория вероятностей. М.: Наука, 1999
Основы надежности устройств ЖАТС. Учебное пособие – М.: РГОТУПС,
2003. 89 с.
задачах. Учебное пособие – М.: МИИТ,2009.-98 с.
и определения» - М.: Издательство стандартов, 1990.
железнодорожной автоматики и телемеханики. 2-е изд., перераб. и доп.- М.:
Транспорт, 1986. 144 с.
1977. 536 с.
машин и систем: Учеб.пособие для вузов по спец. «Вычислительные
машины, комплексы, системы и сети».- М.:Высш.шк., 1989.-216 с.
назначения. Справочник.- М.:ВНИИ Электростандарт, 1989. 196 с.
систем железнодорожной автоматики, телемеханики и связи. Учебное
пособие для вузов ж.д. трансп./ Под ред. Вл.В. Сапожникова.-
М.:Маршрут,2003. 263 с.
обслуживания информационных систем железнодорожного транспорта/ Под
ред. Н.Г. Сачкова: Учебное пособие для вузов ж. – д. транспорта.-
М.:Маршрут,2005.-416 с.
Маликова.- М.: Советское радио, 1972. 408 с.
Учебное пособие/ Под. ред. Г.В. Дружинина.- М.: Энергия, 1976.
вычислительных систем.- М.: Радио и связь, 1991. 132 с.
технических мероприятий по повышению надежности действующих
устройств сигнализации, централизации и блокировки.- М.:МПС, 1989. 80 с.
1984. 79 с.
ПРИЛОЖЕНИЕ 1
ТЕРМИНЫ И ОПРЕДЕЛЕНИЯ, ИСПОЛЬЗУЕМЫЕ В ТЕОРИИ
НАДЕЖНОСИМ
Таблица П1.1
Термины |
Определения |
ОБЩИЕ
Е Dependability
Исправность
Неисправность
|
ПОНЯТИЯ Свойство объекта сохранять во времени в Пр и м еч ани е . Надежность является Свойство объекта непрерывно сохранять Свойство объекта сохранять работоспособное Свойство объекта, заключающееся в Свойство объекта сохранять значения Состояние объекта, при котором он Состояние объекта, при котором он не |
Термины |
Определения |
8. Работоспособное состояние Работоспособность E. State of serviceability |
Состояние объекта, при котором значения |
технической и (или) конструкторской документации | |
9. Неработоспособное состояние Неработоспособность E. State of no serviceability | Состояние объекта, при котором значение |
10. Предельное состояние E. Marginal state | Состояние объекта, при котором его |
11. Критерий предельного E. Marginal state criterion | Признак или совокупность признаков (или) конструкторской документации |
12. Дефект E. Defect | По ГОСТ 15467—79 |
13. Повреждение E. Damage | Событие, заключающееся в нарушении |
14. Отказ E. Failure | Событие, заключающееся в нарушении |
15. Критерий отказа E. Failure criterion | Признак или совокупность признаков (или) конструкторской документации |
16. Техническое обслуживание E. Maintenance | По ГОСТ 18322—78 |
работоспособного состояния Восстановление E. Restoration
| По ГОСТ 18322—78 По ГОСТ 18322—78 |
Термины | Определения |
19. Обслуживаемый объект | Объект, для которого проведение |
20. Необслуживаемый объект | Объект, для которого проведение, |
21. Восстанавливаемый объект E. Restorable item | Объект, для которого в рассматриваемой |
22. Невосстанавливаемый объект E. No restorable item | Объект, для которого в рассматриваемой |
23. Ремонтируемый объект E. Repairable item | Объект, для которого проведение ремонтов (или) конструкторской документации |
24. Неремонтируемый объект E. No repairable item | Объект, для которого проведение ремонтов |
25. Показатель надежности E. Reliability index | Количественная характеристика одного или |
26. Единичный показатель E. Simple reliability index | Показатель надежности, характеризующий |
27. Комплексный показатель E. Integrated reliability index | Показатель надежности, характеризующий |
28. Наработка E. Operating time | Продолжительность или объем работы |
29. Наработка до отказа E. Time to failure | Наработка объекта от начала его |
Термины | Определения |
30. Наработка между отказами E. Time between failures | Наработка объекта от окончания |
31. Технический ресурс Ресурс | Наработка объекта от начала его |
32. Срок службы | Календарная продолжительность от начала |
33. Срок сохраняемости E. Storageability time | Календарная продолжительность хранения Примечание. Различают сроки |
34. Время восстановления Время восстановления E. Restoration time | Продолжительность восстановления работоспособного состояния объекта |
ПОКАЗАТЕЛИ Показатели 35. Вероятность безотказной E. Survival probability | НАДЕЖНОСТИ Вероятность того, что в пределах заданной Примечание к терминам 35 - 40. Для вероятность безотказного хранения (транспортирования), среднее время хранения |
36. Средняя наработка до отказа E. Mean time to failure | Математическое ожидание наработки объекта до первого отказа |
37. Гамма - процентная E. γ-percentile time to failure | Наработка, в течение которой отказ объекта |
Термины | Определения |
38. Средняя наработка на отказ Наработка на отказ E. Mean time between failures |
Отношение наработки восстанавливаемого |
39. Интенсивность отказов Ндп. Опасность отказов E. Failure rate | Условная плотность вероятности |
40. Параметр потока отказов Ндп. Средняя частота отказов Интенсивность потока отказов Е. Failure rate | Отношение среднего числа отказов |
Показатели 41. Средний ресурс | долговечности Математическое ожидание ресурса. Примечание к терминам 41 - 46. В объекта |
42. Гамма - процентный ресурс | Наработка, в течение которой объект не |
43. Назначенный ресурс | Суммарная наработка объекта, при достижении которой применение по назначению должно быть прекращено |
44. Средний срок службы E. Mean lifetime | Математическое ожидание срока службы |
45. Гамма - процентный срок E. γ - percentile lifetime | Календарная продолжительность от начала |
Термины | Определения |
46. Назначенный срок службы E. Specified lifetime | Календарная продолжительность |
Показатели рем
Вероятность восстановления
восстановления работоспособного состояния Среднее время восстановления
| быть прекращено онтопригодности Вероятность того, что время восстановления работоспособного состояния Математическое ожидание времени восстановления работоспособного состояния. Примечание . Аналогично определяется восстановления работоспособного состояния» |
Показатели сохраняемости | |
Комплексные
Ндп. Коэффициент
| Математическое ожидание срока Срок сохраняемости, достигаемый показатели надежности Вероятность того, что объект окажется в Вероятность того, что объект окажется в Отношение математического ожидания |
Термины | Определения |
54. Коэффициент планируемого | Доля периода эксплуатации, в течение |
ремонте | |
55. Коэффициент сохранения | Отношение значения показателя |
ОТК 56. Независимый отказ E. Primary failure | АЗЫ Отказ объекта, не обусловленный отказом |
| Отказ объекта, обусловленный отказом Отказ, характеризующийся скачкообразным |
59. Постепенный отказ E. Gradual failure | Отказ, характеризующийся постепенным |
60. Перемежающийся отказ E. Intermittent failure | Многократно возникающий |
61. Конструкционный отказ E. Design-error failure | Отказ, возникший в результате |
62. Производственный отказ E. Manufacture error failure | Отказ, возникший в результате |
| Отказ, возникший в результате нарушения Явления, процессы, события и состояния, Явления, процессы, события и состояния, |
Термины | Определения |
РЕЗЕРВИРОВАНИЕ
66. Резервирование I Применение дополнительных средств и
E. Redundancy | (или) возможностей с целью сохранения отказе одного или нескольких его элементов |
67. Основной элемент | Элемент структуры объекта, необходимой |
68. Резервный элемент | Элемент объекта, предназначенный для |
69. Резервируемый элемент E. Element under redundancy | Основной элемент, на случай отказа |
70. Резерв Ндп. Избыточность E. Reserve | Совокупность дополнительных средств и |
71. Нагруженный резерв Ндп. Горячий резерв E. Loaded reserve | Резерв, который содержит один или |
72. Облегченный резерв Ндп. Теплый резерв E. Reduced reserve | Резерв, который содержит один или |
73. Ненагруженный резерв Ндп. Холодный резерв E. Unloaded reserve | Резерв, который содержит один или |
74. Восстанавливаемый резерв E. Restorable reserve | Резерв, который содержит один или |
75. Невосстанавливаемый резерв E. No restorable reserve | Резерв, который содержит один или несколько |
76. Кратность резерва Кратность Е. Redundancy multiplicity | Отношение числа резервных элементов |
Термины | Определения |
77. Структурное резервирование E. Structural redundancy |
Резервирование с применением резервных |
78. Временное резервирование E. Temporary redundancy |
Резервирование с применением резервов |
79. Информационное E. Information redundancy |
Резервирование с применением резервов |
80. Функциональное E. Functional redundancy |
Резервирование с применением функциональных резервов |
81. Нагрузочное резервирование |
Резервирование с применением нагрузочных резервов |
82. Общее резервирование E. Whole-system redundancy |
Резервирование, при котором |
83. Раздельное резервирование E. Inherent redundancy |
Резервирование, при котором |
84. Постоянное резервирование E. Constant redundancy |
Резервирование без перестройки структуры |
85. Динамическое E. Dynamic redundancy |
Резервирование с перестройкой структуры |
86. Резервирование замещением E. Stand-by redundancy |
Динамическое резервирование, при котором |
87. Скользящее резервирование |
Резервирование замещением, при котором |
88. Смешанное резервирование E. Combined redundancy |
Сочетание различных видов резервирования в одном и том же объекте |
89. Дублирование Е. Duplication |
Резервирование с кратностью резервирования один к одному |
ПРИЛОЖЕНИЕ 2
ИНТЕНСИВНОСТИ ОТКАЗОВ ЭЛЕМЕНТОВ
Таблица П2.1
Интенсивности отказов резисторов в номинальном режиме
(T = +20 C; KH = 1)
Типы | Номинальная мощность рассеивания PHOM , Вт | |||||||||||||
резисторов | 0,25 | 0,5 | 1 | 2 | 5 | 10 | 15 | 20 | 25 | 30 | 50 | 60 | 75 | 100 |
Интенсивность отказов 2-10 | 6 1/ч | |||||||||||||
МЛТ | 0,4 | 0,5 | 1,0 | 1,6 | - | - | - | - | - | - | - | - | - | - |
ТВО | 0,4 | 0,45 | 0,8 | 1,4 | 2,2 | 3,0 | - | 4,0 | - | - | - | 6,0 | - | - |
МОУ | 0,5 | 0,55 | 1,1 | 1,5 | 2,3 | 3, | - | - | 4,2 | - | 5,5 | - | - | 10 |
МУН | 0,8 | 0,6 | 1,2 | 2,0 | - | -1 | - | - | - | - | - | - | - | - |
УНУ | 0,6 | 0,7 | 1,2 | 1,7 | 2,3 | 3,0 | - | - | 4,8 | - | 8,0 | - | - | 12 |
КЭВ | 0,6 | 0,75 | 1,3 | 1,75 | 2,4 | 3,1 | - | - | 5,0 | - | - | - | - | - |
ВС | 0,7 | 0,8 | 1,35 | 1,8 | 2,5 | 3, | - | - | - | - | - | - | - | - |
УЛИ | 0,6 | 0,65 | 1,3 | - | - | -3 | - | - | - | - | - | - | - | - |
БЛП | 0,6 | 0.75 | 1,4 | - | - | - | - | - | - | - | - | - | - | - |
СПО | 0,6 | 0,7 | 1,15 | 1,8 | - | - | - | - | - | - | - | - | - | - |
СП | 0,7 | 0,8 | 1,3 | 2,0 | - | - | - | - | - | - | - | - | - | |
ПТН | - | 1,1 | 1,4 | 1,8 | - | - | - | - | - | - | - | -- | - | - |
ПКВ | - | 1,2 | 1,5 | 2,0 | 2,5 | - | - | - | - | - | - | - | - | - |
ПЭВ | - | 1,6 | 1,5 | 2,0 | 2,5 | - | - | - | - | - | - | - | - | - |
ПТП | - | - | 2,2 | 2,6 | 3,0 | - | -- | - | - | - | - | - | - | - |
РП | - | - | - | 3,0 | - | - | - | - | 4,7 | - | - | - | 8,5 | - |
Таблица П2.2
Интенсивности отказов конденсаторов в номинальном режиме
T = +20 C; KH = 1)
Тип конденсатора |
Интенсивность отказов Л-10 6 1/ч |
Бумажные |
1,8 |
Металлобумажные |
2,0 |
Слюдяные |
1,2 |
Стеклянные |
1,6 |
Керамические |
1,4 |
Пленочные |
2,0 |
Электролитические алюминиевые |
2,4 |
Электролитические танталовые |
2,2 |
Оксидно-электролитические |
1,3 |
Подстроечные |
2,0 |
Таблица П2.3
Интенсивности отказов полупроводниковых приборов и компонентов
электронных схем в номинальном режиме
Полупроводниковые приборы | Интенсивность отказов |
Диоды | |
Выпрямительные точечные | |
германиевые | 0,7 |
кремниевые | 2,0 |
Выпрямительные микроплоскостные германиевые | 0,7 |
Выпрямители плоскостные германиевые | 5,0 |
Импульсные точечные германиевые | 3,0 |
Импульсные плоскостные мезадиоды | |
германиевые | 2,0 |
кремниевые | 2,5 |
Импульсные сплавные кремниевые | 0,6 |
Управляемые кремниевые | 5,0 |
Стабилитроны кремниевые | 5,0 |
Варикапы кремниевые | 5,0 |
Выпрямительные столбы кремниевые | 5,0 |
Микромодульные | |
германиевые | 4,2 |
кремниевые | 4,5 |
Диоды излучающие | 0,58 |
Диоды и диодные сборки | 0,7 |
Тиристоры | 0,45 |
Оптопары | 1,0 |
Транзисторы | |
Маломощные низкочастотные | |
германиевые | 3,0 |
кремниевые | 4,0 |
Мощные низкочастотные германиевые | 4,6 |
Маломощные высокочастотные германиевые | 2,6 |
Мощные высокочастотные | |
германиевые | 5,0 |
кремниевые | 1,7 |
Микромодульные германиевые | 1,0 |
Микросхемы | |
гибридные | 1,0 |
полупроводниковые | 0,8 |
однокристальные ЭВМ и микропроцессоры | 1,0 |
оптоэлектронные | 0,22 |
Коммутационные устройства | 2,8 |
Разъемы | 9,1 |
Соединения пайкой | 0,001 |
Провода соединительные | 0,012 |
Кабели | 0,22 |
Предохранители плавкие | 0,16 |
Контакты магнитоуправляемые | 0,03 |
Таблица П2.4
Интенсивности отказов трансформаторов и моточных изделий (дроссели,
катушки, обмотки) в номинальном режиме
Трансформаторы и моточные изделия |
Интенсивность отказов λ⋅ 10-6 1/ч |
Автотрансформаторы Трансформаторы |
5,0 |
силовые |
3,0 |
высоковольтные |
4,0 |
импульсные |
0,5 |
межкаскадные |
0,11 |
Дроссели |
1,0 |
Катушки индуктивности |
0,5 |
Электродвигатели |
80 |
Таблица П2.5
Интенсивности отказов устройств автоматики и телемеханики
Тип устройства | Интенсивность отказов |
Преобразователь ПЧ 25/50 | 0,9 |
Блок БПШ | 0,9 |
Известительное реле | 1,2 |
Схемы кодирования в хвост удаляющегося поезда | 1,2 |
Электродвигатель | 1,5 |
Пневматическая обдувка | 1,5 |
Кодовый путевой трансмиттер | 1,8 |
Сопротивление: | |
цепи ПХ-ОХ относительно земли | 1,8 |
изоляции линии между релейным шкафом и | 2,7 |
дроссельным трансформатором | |
изоляции кабеля между релейным шкафом и КЯ в | 2,7 |
высоковольтной линии | |
Изоляция монтажа кабеля | 2,9 |
Основной источник питания (ПХ-ОХ) на посту ЭЦ | 2,9 |
Резервный источник переменного тока и предохранитель | 2,9 |
Схемы изменения направления | 2,9 |
Релейные схемы ЭЦ | 2,9 |
Перекрытие входных, поездных и маневровых сигналов | 3,3 |
Емкость конденсаторов | 4,3 |
Аккумуляторы батарей: | |
рабочих | 4,5 |
станционных при максимальной нагрузке | 4,5 |
входного светофора | 4,5 |
Электрические стрелки | 4,6 |
Лампы сигналов автошлагбаума | 5,5 |
Цепь двигателя автошлагбаума | 5,5 |
Схемы открытия (закрытия) автошлагбаума | 5,5 |
Автошлагбаумы | 5,5 |
Стабилитронный блок | 5,7 |
Реле: | |
путевые |
6,4 |
ДСШ |
6,4 |
Предохранители |
6,4 |
Цепи двойного снижения напряжения |
6,4 |
Трансмиттерное реле |
8,9 |
Рельсовая цепь |
9,0 |
Дешифраторные ячейки |
9,2 |
Лампы светофоров: | |
входных |
14,9 |
станционных |
14,9 |
Устройства звуковой и световой сигнализации |
14,9 |
Изолирующие стыки |
25,5 |
Реле контроля и защиты |
0,6142 |
Контакторы |
0,0287 |
Автоматы защиты и выключатели автоматические |
0,303 |
ПРИЛОЖЕНИЕ 3
ВЫРАЖЕНИЯ ДЛЯ РАСЧЕТА НАДЕЖНОСТИ МАЖОРИТАРНЫХ
СТРУКТУР
n | d | Логическая | Надежность (n,1) - структур | |
P0(t) | P1(t) | |||
2 | 1 | M21 | 2 Ро - P0 | p12 |
2 | M22 | p02 | 2 Р1 - P12 | |
3 | 1 | M31 | 3 Р 0 - 3 Р 2 + Р 0 | p13 |
2 | M32 | 3 p 2- 2 p 0 | 3 Р2 - 2 Р3 | |
3 | M33 | p03 | 3 Р1 - 3 Р1 + Р3 | |
4 | 1 | M14 | 4 Р 0 - 6 Р 2 + 4 Р 0 - Р 4 | p14 |
2 | M42 | 6 Р2 - 8 Р 3 +3 Р4 | 4 Р13 - 3 Р14 | |
3 | M43 | 4 Р 0 - 3 Р 0 | 6 Р12 - 8 Р13 + 3 Р14 | |
4 | M44 | p04 | 4 Р1 - 6 Р12 + 4 Р13 - Р14 | |
5 | 1 | M51 | 5 Р 0 - 10 Р2 + 10 Р 0 - 5 Р 0 + + Р 0 | p15 |
2 | M52 | 10Р2 - 20Р0 +15 Р0 - 4Р0 | 5 Р14 - 4 Р15 | |
3 | M53 | 10 Р 0 - 15 Р 4 + 6 Р 5 | 10Р3 -15 Р4 + 6Р5 | |
4 | M54 | 5 Р 4 - 4 Р 0 | 10Р12 - 20Р13 +15Р14 - 4Р15 | |
5 | M55 | p05 | 5 Р1 -10Р2 +10Р3 - 5 Р4 + + Р15 |
ПРИЛОЖЕНИЕ 4
ФОРМУЛЫ ПРЕОБРАЗОВАНИЙ ЛАПЛАСА
Таблица П4.1
Изображение | Оригинал | |
as1 (x) + Ps 2 (x) | aP1( t ) + e 2( t ) | |
- p (0) + sa (x) | p'(t) | |
s (x) C xx | t j p (x) dx + C 0 | |
s [ X} b > 0 | p (bx) | |
e’bxs(x), b > 0 | P(t- b) | |
s1( x ) s 2( x ) | t j p1 (x) p 2 (t - x) dx 0 | |
sk (x) | (-1) ktkp (t) | |
s (x - a) | eap (t) | |
—, n = 1,2,... | tn-1 (n -1)! | |
----1----, n = 1,2,... | . n-1 t ebt (n -1)! | |
1 2 , z2 | 1 • ь — sin bt | |
1 x 2 - b 2 | cos bt | |
x 2 l2 x - b | 1 sh bt | |
1 x (x 2 + b 2) | chbt | |
1 x (x - b) | Aebt + K | A = 1, K = -1 |
x + d x (x - b) | , b + d _ d A= b ’K = -b | |
1 (x - b)(x - c) | Aebt + Bect | A = —, B = — |
x + d (x - b)(x - c) | , b + d „ c + d b — c c — d |
Продолжение табл.П4.1
Изображение | Оригинал | |
1 x (x - b)(x - c) | Aebt + Bec + K | A =--1--, B = —1—, b (b - c) c (c - b) K = — |
x + d x (x - b)(x - c) | Aebt + Bect + K | . b + d „ c + d A =-------, B =-------, b (b - c) c (c - d) K = — |
x2 + gx + d | A = b2 + gb + d B = c2 + gc + d k = — | |
1 | Aeb + Bec + Cedt | 1 |
(x - b)(x - c)(x - d) | (b - c)(b - d), B =------1------, (c - b)(c - d) C =-----1----- (d - b)(d - c) | |
x + e | . b + e A =------------. | |
(x - b)(x - c)(x - d) | (b - c)(b - d) B =---c+e---, (c - b)(c - d) c = d + e (d - b)(d - c) | |
2 . . x + gx + e (x - b)(x - c)(x - d) | A= b2 + gb + e C d2 + gd + e = (d - b)(d - c) |
СОДЕРЖАНИЕ
ВВЕДЕНИЕ 4
НЕРЕЗЕРВИРОВАННЫХ СИСТЕМ 35
РЕЗЕРВИРОВАННЫХ СИСТЕМ 52
СИСТЕМ 83
ЛИТЕРАТУРА 117
ПРИЛОЖЕНИЕ 1 119
ПРИЛОЖЕНИЕ 2 129
ПРИЛОЖЕНИЕ 3 133
ПРИЛОЖЕНИЕ 4 134
СОДЕРЖАНИЕ 136
Комментарии (0)