книги / Надежность судовой электронной аппаратуры и систем автоматического управления
..pdf0,99, то повысить надежность на 17% можно, уменьшив число элементов с 40 до 20.
Однако уменьшение числа элементов в аппаратуре, как и уве личение их надежности, является сложной технической задачей.
Так как аппаратура может находиться либо* в исправном со стоянии, либо в состоянии отказа, то сумма вероятности безотказной работы Р (/) и вероятности повреждения Q (t) всегда равна единице
р (0 + Q (0 = 1.
На основании этого вероятность повреждения, аппаратуры определяется соотношением ^
Q (0 = 1 - P (t) = 1 -
— Pi (t) Рг |
i f ) |
........ |
PN it)- |
|
|
|
|
(33) |
|
Выражая |
Q (t) |
системы |
|
|
через вероятность поврежде |
|
|||
ния элементов, |
получим |
|
||
Q (0 — l - [ l - < |
7 i |
(t)] [1 - |
|
|
q%(01 ,. . ., |
[1 ~ ам (Й ]. |
определению условной вероят- |
||
|
|
|
Рис. 18. К |
|
|
|
|
(34) |
ности. |
Если надежность элементов одинакова, то формула (34) приобре тает вид
Q (/) = 1 — [1 — ^ |
(35) |
В частном случае, когда надежность оценивается за достаточно короткий промежуток времени, вероятность повреждений ‘эле ментов практически оказывается весьма малой величиной, т. е*
<7i (О <С 1; |
( / R 1; |
gN (0 < 1. |
В этом случае, пренебрегая членами высшего порядка малости, формулу (34) можно представить таким приближенным равенством:
Q ( 0 « 1 -{1 - [<7i (0 + <72 ( t ) + . . . + q N \ t )]} = J <7*( t ) . (36) i=l
Следовательно, при весьма малой вероятности повреждения элементов вероятность отказа аппаратуры за 'Небольшой проме-„ жуток времени приблизительно равна сумме вероятностей отказа ее элементов за тот же промежуток времени. При одинаковой на дежности элементов вероятность повреждения аппаратуры будет
Q (t) ^ Nq (t). |
(37) |
Если требуется определить вероятность безотказной^ работы изделия (рис. 18), проработавшего время tlt за следующий проме-
л* |
51 |
жуток времени от до / 2, то для этого необходимо воспользоваться условной вероятностью, равной отношению вероятности того, что
изделие безотказно проработает время от 0 до |
t 2 к вероятности |
безотказной работы изделия за время от 0 до |
т. е. |
Р ( к , ^ = Щ - |
(38) |
Определение условной вероятности поясняется следующим примером. Пусть нам известна функция надежности P (t). Тре буется определитьвероятность безотказной работы элемента, проработавшего время tlt за промежуток времени от t x до t 2-
Определив по Графику значения вероятностей Р (/2) = 0,4 и P (*i) = 0,6, вычисляем условную вероятность -
Я ft, к) = |
= 0,67. |
Условная вероятность для изделия определяется так же, как
идля элементов.
Взаключение отметим что такая количественная характери стика надежности, как вероятность безотказной работы, обладает следующими достоинствами:
1)позволяет характеризовать изменение надежности во вре
мени; 2) дает возможность достаточно наглядно судить о надежности
изделий, ибо чем больше вероятность исправной работы тех или иных элементов, тем они надежнее;
3)может быть использована для расчета надежности сложных систем до их изготовления, что очень важно при проектировании новой аппаратуры;
4) охватывает большинство факторов, существенно влияющих на надежность изделия, т. е. достаточно обширна.
§ 18Среднее время безотказной работы
Под средним временем безотказной (исправной) работы в тео рии надежности понимается математическое ожидание времени работы изделия до отказа.
Среднее время безотказной работы обозначается символом Тср. Зная дифференциальный закон распределения, можно опреде лить и математическое ожидание случайной величины, т. е. сред нее время исправной работы, которое определяется зависимостью
00 оо
гср = J fdQ (t) = j tq (t) dt.
--- 00 ----00
Так как такая случайная величина как время не может иметь отрицательных значений, то пределы интегрирования для нашего
• оо
Гс = J t v (t) dt.
о
Подставив в выражение (39) вместо Q' (t) производную от ве: роятности безотказной работы с обратным знаком и выполнив интегрирование по частям, получим
оо 00 оо
J |
tP' (t) dt |
tP(t) I + J P (t) dt. |
ü |
oo |
|
|
|
|
Учитывая, что |
tP (t) | = 0, |
можем записать |
|
о |
|
00
о
Выражение (40) устанавливает зависимость между средним временем и вероятностью безотказной работы. Из этой формулы видно, что Тср есть плошадь под кривой вероятности безотказной работы.
Среднее время безотказной работы однотипных элементов при близительно можно вычислить по формуле
(41)
где tt — время исправной работы i-го элемента; N 0 — число испытываемых элементов.
Чем больше N 0, тем точнее определяется среднее время исправ ной работы элементов. Как видно из формулы (41), для определе ния среднего времени исправной работы необходимо знать моменты выхода из строя всех элементов, над которыми проводится опыт. Поэтому для вычисления Тср пользоваться указанной формулой неудобно. Имея данные о количестве вышедших из строя элемен тов в каждом интервале времени, среднее время исправной работы лучше определять по формуле
m
(42)
В выражении (42) /ср. и m определяются по следующим фор
мулам:
4 _ tt-i 4- U
fcpi — 2 *
гДе tt_t — время начала i-то интервала; ~ti — время конца i-то интервала
где tk — время, в течение которого вышли из строя все элементы; A t — интервал времени.
Для пояснения способа определения среднего времени исправ ной работы рассмотрим следующий пример.
Пример. Определим среднее время исправной работы ламп 6Ж4, установленных в электронной вычислитель
ной машине. Для определения |
Тср |
воспользуемся -дан |
||
ными табл. 2 и формулой (42): |
|
|
|
|
_ |
20*500 -(-25* 1500 -j- • • • -j-20*25500 |
« n QQC |
, |
|
J cp — |
----------------------------- 7ÔÜÔ---------------------------- |
= |
1 1 б У 0 |
Ч а с * |
Среднее время исправной работы элементов является одной из наиболее наглядных количественных характеристик надежности.
юоо |
2000 тср зооо |
то |
ь |
Рис. 19. Распределение числа отказов |
элементов |
||
во |
времени. |
|
|
По ее величине, можно судить о надежности элементов, быстро определить количество запасных элементов для работы аппара туры в течение календарного времени и т. д. Однако эта харак теристика надежности имеет свои недостатки. Как математическое ожидание случайной величины, она не может полностью характе ризовать время работы элемента. Поясним это на примере.
Пример. Пусть у двух групп однотипных элементов среднее время исправной работы одинаково. Все элементы первой группы выходят из строя в интервале времени 2000—3000 час., распределение отказавших элементов п показано на рис. 19, а. Элементы второй группы выходят из строя в интервале времени 500—4100 час.; их распре деление приведено на рис. 19, б.
Несмотря на то, что среднее время исправной работы элементов обеих групп равно между собой, надежность элементов первой и второй групп различна. Если для эле ментов первой группы с уверенностью можно сказать,
что до 2000 час. работы отказов не будет, то для элемен тов второй группы этого утверждать нельзя.
Среднее время исправной работы удобно для оценки надеж ности простейших элементов, которые после отказа не ремонти руются. Его можно использовать и для оценки надежности слож ных устройств, но в последнем случае эта характеристика будет отражать величину среднего времени исправной работы множества систем до первого их отказа. Среднее время исправной работы систем при их использовании после ремонта в общем случае будет отличаться от Тср и совпадет с ним лишь в том случае, если от казы составных элементов системы не зависят от времени их ра боты, т. е. если можно считать, что элементы не изнашиваются
§ 19. Среднее время между соседними отказами (наработка на отказ)
Наработкой на отказ называется среднее время работы ремон тируемого изделия между двумя соседними отказами. Наработка на отказ обозначается tcp и определяется из выражения
|
|
|
П |
|
|
h + |
^2 + •••~Ь tn |
2 |
|
|
|
/=1 |
|
(43) |
|||
рсР |
п |
|
п |
||
где п — число отказов аппаратуры за время |
испытания |
t\ |
|||
ti т— время исправной |
работы |
аппаратуры между |
(i — 1) |
||
и i-м отказами аппаратуры. |
|
tcр будет стремиться |
|||
При достаточно большом числе |
отказов |
к среднему времени между двумя соседними отказами. Если ис пытание проводится не с одним, а с несколькими образцами ап паратуры, то среднее время между двумя соседними отказами можно определить из выражения
|
м |
|
|
*сР= |
м |
* |
(44) |
где t cP[ — среднее время исправной |
работы между |
двумя сосед |
ними отказами i-ro образца аппаратуры, вычисленное по формуле (43);
М — число испытываемых образцов.
Кроме формулы (44), можно пользоваться выражением |
|
||
м |
nj |
|
|
2 *2 |
|
|
|
*ср = |
У |
» |
(45) |
где п* — общее число отказов при испытании М образцов аппа ратуры
Л1
и* = 2 Щ, /=i
tj. — время исправной работы между двумя соседними отка
зами /-го образца аппаратуры; п,- — число отказов /-го образца аппаратуры.
Способ определения среднего времени между соседними отка зами на практике можно показать на следующем примере.
Пример. В результате испытания 36 образцов аппара туры была зафиксирована 231 неисправность. При этом каждый из испытываемых образцов исправно проработал
19час.
По формуле (45) рассчитываем среднее время между
соседними отказами
' |
19-36 |
„ |
^ср — |
по] — |
час. |
Среднее время между двумя соседними отказами аппаратуры, как и среднее время ее исправной работы, является одной из на глядных характеристик аппаратуры и они идентичны при соблю дении условий эргодичности.
§ 20. Частота отказов
Под частотой отказов понимается плотность вероятности времени работы изделия до первого отказа. Статистически оно определяется как отношение числа отказавших изделий в единицу времени к первоначальному числу испытываемых изделий при усло вии, что все вышедшие из строя изделия не восполняются, т. е. число испытываемых изделий во время испытания уменьшается.
Частота отказов обозначается символом a (t).
Если в процессе испытаний N Qэлементов на надежность фик сировать число отказов п (t), происшедших в определенные ин тервалы времени At, то частота отказов определяется из уравне ния
|
= |
<«> |
где п (t) — число |
отказавших изделий |
в интервале времени от |
. A t |
. . At |
|
t - - 5 - Д° * + -«-;
N 0 — первоначальное число испытываемых изделий; At — интервал времени.
Установим зависимость между частотой отказов и вероят ностью безотказной работы P (t). Число элементов, которые будут
исправно работать к моменту времени t, в среднем равно |
(47) |
N.(t) = N 0P(t). |
Число элементов, которые будут исправно работать к моменту времени t + At, равно
N (t + At) = N 0P {t + А/)..
Таким образом, число отказавших элементов можно опреде
лить |
как разность |
N (t) — N (t + At) |
= N о [P (0 — |
|||
|
п (t) = |
|||||
|
|
— P (t -b A01. |
|
(48) |
||
Подставив выражение |
(48) в (46), получим |
|
||||
|
|
„ /,Л _ |
N0[ P ( t ) - P ( t + At)) |
|
||
|
|
а У>— |
|
N0At |
|
|
Устремив At к |
нулю |
и перейдя к пределу, получим |
||||
|
a (t) — — lim- |
Р (t + At) — P (О _ |
(49) |
|||
|
At |
= - Р ' (о |
||||
|
|
дt |
|
|
|
|
или |
|
a (t) |
= —P' {t) = |
Q' (0. |
|
|
|
|
|
||||
Таким образом, функция частоты отказов a (t) есть производ |
||||||
ная |
от функции |
вероятности безотказной работы |
P (t), взятая |
с обратным знаком. Она характеризует скорость снижения ве роятности безотказной работы во времени.
Так как Q (/) есть закон распределения вероятностей времени исправной работы элемента Г, производная от которого представ ляет собой плотность распределения вероятностей случайной величины Г, то функция частоты отказов a (t) есть дифференци альный закон распределения случайной величины Т.
Выполнив интегрирование выражения (49) от нуля до /, по
лучим
t
Q (t) = J а (t) dt
P {t) = 1 — j a (t) dt. |
(50) |
о |
|
Таким образом, между частотой отказов, вероятностью без" отказной работы и вероятностью отказов при любом законе рас" пределения времени возникновения отказов существуют однознач ные зависимости.
Пример. Требуется определить зависимость частоты отказов от времени для ламп 6Ж4, установленных в ци фровой вычислительной машине.
Распределение |
отказавших ламп во времени дано |
в табл. 2 (см. стр. |
48). |
На основании формулы (46) и табл. 2 получаем:
а (50°) = тдддтщ} = °,2- КГ4,
а (1500) = |
0,25-10"4, |
а (2500) = |
0,35* 10~4, |
а (25 500) = |
0,2 *10"4. |
Достоинством этой количественной характеристики надеж ности является, во-первых, возможность судить по ее величине о числе изделий, которые могут выйти из строя на каком-то про межутке времени. Во-вторых, частота отказов позволяет весьма просто определить количество отказавших образцов аппаратуры в промежутке времени At.
Действительно, из выражения (46) можно определить необхо димое количество запасных элементов по формуле
п (t) = a'(t)N0 At. |
(51) |
Следовательно, если известно а (t), то известно и п (t), что по зволяет обосновать необходимое число запасных элементов для обеспечения нормальной эксплуатации аппаратуры.
Например, если известно, что a (t) = 0,2*10-4, а в устройстве работает 100 элементов; то за 500 час. работы из строя может выйти п = 0,2* 10"4 *500* 100 = 1, или один элемент.
Однако, используя величину частоты отказов, можно придти к неправильному представлению о надежности тех или иных элементов.
На первый взгляд кажется, что чем больше частота отказов, тем ненадежнее элементы. Предположим, что частота отказов ламп известна. На основании табл. 2 частота отказов в первую тысячу часов работы равна 0,2*10"4; этой же величине она равна на по следней тысяче часов. Следовательно, может показаться, что на дежность ламп в первую и последнюю тысячу часов их работы одинакова. На самом деле надежность ламп с течением времени падает: если в течение первой тысячи часов вышло из строя 20 ламп из 1000, то на последнем участке времени вышло из строя 20 ламп из 20. Следовательно, надежность ламп для последнего уча стка во много раз меньше, чем для первого.
Типичная кривая изменения частоты отказов аппаратуры во времени показана на рис. 20. Из кривой видно, что в работе аппа ратуры можно выделить три характерных участка.
За время от 0 до tx частота отказов вначале растет, а затем резко снижается. Это объясняется тем, что в начальный период эксплу атации число отказов повышено за счет элементов, имевших вну тренние дефекты.
Этот |
участок обычно |
называют п е р и о д о м п р и р а |
б о т к и |
э л е м е н т о в . |
Если элементы проходят предварительг |
ную тренировку, то этот участок отсутствует. На участке времени от t 1 до 12 частота отказов уменьшается незначительно. Этот уча сток характеризует нормальную работу элементов. Рост кривой частоты отказов на участке вре- а
мени t 2 — t3 объясняется механи ческим или .электрическим изно сом элементов. Падение кривой частоты отказов после t3 объяс няется не увеличением надежно сти элементов, а незначительным количеством исправно работающих к этому времени элементов, вслед ствие чего число отказавших эле ментов тоже будет небольшим.
Частота отказов полно характеризует надежность изделий до их первого отказа. Все отказавшие изделия в дальнейших испы таниях не участвуют. Поэтому частоту отказов можно использо вать лишь для оценки надежности тех изделий, которые после воз никновения отказов не ремонтируются и в дальнейшем не эксплу атируются.
§ 21. Опасность (интенсивность) отказов
Опасностью отказов называют отношение частоты отказов к вероятности безотказной работы невосстанавливаемого при испытании изделия, взятых для одного и того же момента вре мени; статистически определяется как отношение числа отка завших изделий в единицу времени к среднему числу изделий, ис правно работающих в данный отрезок времени.
Опасность отказов обозначается символом X (t) и определяется
по формуле |
|
|
|
|
|
|
|
|
|
|
|
. |
X W = J T W K Ï’ |
|
|
<52> |
|||
где п |
(t) — число отказавших |
изделий за |
промежуток |
времени |
|||||
|
|
|
. |
д/ |
, . |
д t |
Ï |
|
|
|
|
от t -----2“ |
ДО * + |
~ |
|
|
|||
|
At — интервал |
времени; |
Nt-i + Ni |
|
|
|
|||
|
|
|
N(t) = |
|
|
|
|||
|
|
|
|
|
2 |
|
|
|
|
где |
Ntmmi — число |
исправно |
работающих изделий |
в |
начале |
||||
|
интервала |
А*; |
|
|
|
изделий в |
конце ин |
||
|
Nt — число |
исправно работающих |
|||||||
|
тервала At. |
|
|
|
|
|
|
Выражение (52) является статистическим определением опас ности отказов. Для вероятностного представления этой характе ристики установим зависимость между опасностью отказов, вероятностью безотказной работы и частотой отказов.
Для нахождения этой взаимосвязи, подставив в формулу (52) значения N (t) и п (t) из (47) и (48), получим
г т |
_ Mo[P ( t ) - P ( t + At)) |
Л W |
N0P (/) Дt |
Устремив Лt к нулю и перейдя к пределу,
%(t) ■— — lim тP ( t + à t ) - P ( t ) |
|
о |
P (t) Дt |
получим
P'(t)
P (t)
или
(53)
Таким образом, опасность отказов выражается через отноше ние частоты отказов к вероятности безотказной работы.
После интегрирования в пределах от 0 до /, получим
J X(t) dt = — lnP(t),
При X = const вероятность исправной работы аппаратуры можно определить в виде
P (t) = в -« . |
(55) |
Из выражений (54) и (55) следует, что, зная опасность отказов элементов как функцию времени, можно определить вероятность исправной работы в течение времени t.
Среднее время исправной работы элементов также можно выразить через опасность отказов. Подставив в выражение (40) значения из (54), получим65
00
(56)
о
Если опасность отказов элементов не зависит от времени, то среднее время исправной работы будет равно
00