
книги из ГПНТБ / Смирнов, О. Р. Надежность судовых энергетических установок
.pdfравна 0,6827. Если таким же образом отложить отрезки длиной 2а, то вероятность наступления отказа в интервале (Т — 2а, Т + 2а) равна 0,9545. И, наконец, в интервале (Т — За, Т + За) эта вероят
ность |
имеет |
значение |
0,9973. |
Если ремонтные |
|
мероприятия или замену |
||
элементов |
установки |
|
производить |
в момент |
|
Т — а, |
то вероятность |
|
отказа до этого момен |
||
та, как это следует из |
||
предыдущего, |
равна |
5=3 |
Т= var |
«*0,16. Если же анало |
|||
гичные |
мероприятия, |
||||
|
5 = const |
||||
|
восстанавливающие ра |
||||
|
|
||||
|
|
ботоспособность элемен |
|||
|
|
та, осуществлять в мо |
|||
|
|
мент |
Т — 2а или Т — |
||
|
|
— За, |
то |
соответствую |
|
|
|
щие значения вероятно |
|||
Р и с . 33 . Х ар ак тер изм ен ен и я кривы х |
плотности |
сти отказа будут равны |
|||
«*0,023 и «*0,0014. От |
|||||
н орм альн ого зак он а . |
|
сюда видно,что в отличие от экспоненциального закона ремонтные, профилактические ме роприятия или замена элемента на аналогичный новый являются действенным методом повышения надежности, причем чем чаще они проводятся, тем выше безотказ
ность элемента. |
|
|
|
|
||
Отметим также, что выраже |
|
|
||||
ние (3.22) |
представляет |
собой |
|
|
||
плотность |
распределения |
слу |
|
|
||
чайной величины, которая мо |
|
|
||||
жет принимать значения на всей |
|
|
||||
числовой оси. Однако время |
|
|
||||
безотказной |
работы не |
может |
|
|
||
быть отрицательным. Рисунок 34 |
|
|
||||
показывает, |
Т |
|
ве |
|
|
|
что при — > 3 |
|
|
||||
роятностью отрицательных зна |
Р и с. 3 4 . |
З н а ч ен и е в ер оятн ости отк азов |
||||
в зав и сим ости от и н тер вала врем ени д л я |
||||||
чений t с достаточной, для прак |
||||||
|
н ор м ал ь н ого зак он а . |
тики точностью можно пренеб речь. Указанное соотношение Т и а часто выполняется для элемен
тов СЭУ, поэтому рассмотрим сначала |
количественные характерис- |
||
тики безотказности в предположении, |
т |
|
|
что — > 3 . |
|
||
Вероятность отказа |
|
|
|
Q (f) = j а (т) ск : |
( т -7 Т |
|
|
2а2 |
dx. |
(3.23) |
|
V |
а |
|
|
120
В выражении |
(3.23) |
сделаем |
^_f |
, тогда |
||
подстановку х = — |
||||||
|
|
|
при |
x — t |
х = —-— ; |
|
|
|
|
при г = —оо х — —оо; |
|
||
|
|
|
|
dx = |
adx, |
|
Теперь |
(3.23) можно переписать в виде |
|
||||
|
|
|
|
|
t — T |
|
|
|
|
|
|
|
( 3 ' 2 4 ) |
Можно |
видеть, |
что величина |
х также распределена нормально |
|||
с параметрами Т = |
0 и о — 1. Нормальное распределение с такими |
|||||
параметрами |
называются |
простейшим. |
функции |
|||
Интеграл |
(3.24) |
не выражается через элементарные |
и его вычисляют через специальные функции, выражающие определен- t2
ный интеграл от выражения е~{2 или е 2 (интеграл вероятностей). Для этих функций составлены таблицы (см. например, [42], [50]).
В теории надежности наиболее широко используют две |
функции: |
||||
|
X |
__ t 2 |
|
|
|
F ( x ) = ~ = - |
j e |
2 |
d t — значения |
функции распределения про- |
|
|
X |
е- |
стеишего |
нормального закона; |
|
|
|
|
|
||
|
|
2 |
dt — нормированная функция Лапласа. |
||
Очевидно, |
что |
|
|
|
|
|
|
|
F (х) = -i- -г Ф(х), |
|
|
|
|
F ( - x) = 1 - F ( x) = -F .-< $ (x). |
(3.25) |
||
Теперь вероятность отказа можно записать в виде |
|
||||
■ |
' |
|
|
. |
|
Вероятность безотказной работы в силу (3.25) определится выра
жением. |
|
‘ |
Интенсивность отказов |
|
|
|
|
_ ( t - T ) 2 |
1 /А _ |
_ |
е 2(72 |
121
Рассмотрим теперь случай, когда нельзя пренебречь вероятно*
стями отрицательных значений величины т < 3 ^ . В этом случае
для вычисления характеристик надежности используют усеченное
нормальное распределение, |
плотность которого |
|
|||||
|
|
|
аг (0 |
= са (О, |
|
|
|
где с — нормирующий множитель. |
|
|
|
||||
со |
|
|
|
|
|
|
|
Так как J |
аг (t) dt = 1, то величина с определяется из соотно- |
||||||
о |
|
|
|
|
|
|
|
шения |
|
|
|
|
а-т)г |
|
|
|
|
|
|
СО |
|
|
|
|
|
|
с |
2а2 |
|
|
|
|
|
|
|
dt = |
1. |
||
|
|
|
Уг 2 л |
а |
|||
|
|
|
|
|
|
||
Отсюда |
|
|
|
|
|
|
|
_ |
|
оУг2я |
_ |
]/2я |
_ |
1 |
|
^ |
00 _ |
U - T ) 2 |
|
со |
х 2 |
|
/ X \ |
|
J е |
2° 2 |
d t |
| е |
2 d x |
F \ ~ ( Г ) |
|
|
о |
|
|
г |
|
|
|
|
|
|
|
~а |
|
|
|
Таким образом, частота отказов в случае усеченного нормального |
|||||||
распределения |
имеет вид |
|
|
|
|
|
{ t - T ) 2
2 а 2
<h (t) =
Аналогично предыдущему может быть определен вид количест венных характеристик безотказности для усеченного закона. Эти характеристики приведены в табл. 12, а их графическая интерпре тация на рис. 35, б. Из этих данных видно, что интенсивность отказов невелика при малых t, однако растет с экспоненциальной скоростью с течением времени.
Следует отметить также, что чем меньше величина дисперсии а2, тем длительнее интервал времени, на котором вероятность насту пления отказа мала. Величина дисперсии во многом определяется однородностью качества изготовления элементов данного типа и ус ловий их эксплуатации.
Для описания времени безотказной работы в случае отказов, причиной которых является износ [16], используют гамма-распреде ление. Оно возникает в тех случаях, когда выполнены следующие условия:
—качество изготовления элементов однородно;
—нагрузка в процессе эксплуатации меняется в достаточно ши роких пределах;
—период приработки занимает небольшую часть времени экс плуатации.
Обычно эти условия в практике СЭУ оказываются выполнен ными. Так, изготовление деталей узлов трения, как правило, яв-
122
Таблица 12
Вид количественных характеристик надежности для различных законов распределения времени безотказной работы
Распределение времени безотказ ной работы
П о к а за т е л ь н о е
У сеч ен н о е н о р м ал ь н ое
Гам м а
В е й б у л л а
Частота отказов а ( t )
Хе ~ и
,( t - т у
1 |
с |
2 а ‘ ' |
Х Ы к ~ х е ~ м к
Вероятность безотказной работы р ( t )
е - и
' ( ^ )
i ~ Q
ё ~ и к
Интенсивность отказов
X (t)
X = con st
( t - т у
е2(12
'Г2’" р {т« ' )
Я, ( X t ) k ~ x
*
/ = о
X k t k ~ x
Среднее время безотказной работы Т
1 Д
т + ---------
v ™ ( | )
k
“ Г
г ( х + 0
1
|
|
|
|
|
|
|
|
|
|
|
х к |
|
|
Р эл ея |
t |
— |
С - |
|
t* |
|
|
|
t |
|
|
|
|
_ L |
р |
2а 2 |
|
е |
|
|
|
а 2 |
у |
ъ |
|
||
|
а 2 |
|
|
|
|
|
|
|
|||||
|
, |
|
(In t - и У |
|
|
|
|
|
(In < - ц ) 2 |
со |
(In t —р.)2 |
|
|
Л огар и ф м и ч еск и |
|
{ |
p — l n t |
\ |
1 |
-в |
202 |
|
|||||
1 |
е |
2<т2 |
\ |
о |
) |
- - L - Г е |
2° 2 |
Л |
|||||
н ор м ал ь н ое |
l/"2 n a t |
р |
^М- — In ^ |
||||||||||
■ V 2 n a t |
|
|
|
|
|
J |
( |
|
|||||
|
|
|
|
|
|
|
|
|
|
0 |
|
ляется серийным, что обеспечивает однородность их начального ка чества; нагрузка на трущиеся детали меняется вследствие изменения условий эксплуатации установки (переменные режимы, крен, диф ферент и т. п.), приработка трущихся пар занимает, как было пока зано, незначительную часть времени эксплуатации.
Р и с . |
3 5 . Х ар ак |
тер |
зав и сим остей |
п ок азател ей |
н адеж н ости от врем ени |
|
п ри |
различны х |
зак о н а х р асп р едел ен и я |
врем ени б езо т к а зн о й работы : |
|||
а — |
показательном ; б — норм альном ; |
в — л огар и ф м и чески -норм аль - |
||||
|
ном; |
г — В е й б у л л а |
и гамма; |
д — Р э л е я . |
Частота отказов в случае гамма-распределения определяется вы
ражением
бkfk—l
= |
(3-26) |
где Г (k) — гамма-функция;
Я. и k — параметры распределения.
124
Для целых k справедливо |
равенство |
Г ( к ) = |
{ k — 1)!. |
В этом случае частота отказов (3.26)
а ^ ==Я- ^ Т ) Т е" 4 |
<3-27) |
Из (3.26) и (3.27) видно, что при к — 1 гамма-распределение тождественно показательному. При возрастании k плотность гаммараспределения стремится к плотности нормального закона.
В рассматриваемом случае основные характеристики безотказ ности имеют вид:
t |
|
k |
t |
|
P ( t ) = l - \ a ( x ) d x = l - |
|
J%k~x e ~ * dx = |
||
о |
k-\ |
|
0 |
|
|
ШУ . |
|||
|
— g—xt |
|||
|
2 |
t! |
’ |
|
|
;=0 |
|
|
|
%{t). |
a(t) |
X (Xt)k-i |
||
P(t) |
k-i |
|||
|
(X— 1)! |
s |
(Xt)1 |
|
i ! |
|||
|
|||
|
|
T = f P(t)dt = ~ .
оA
Характер изменения в зависимости от времени этих величин для различных k приведен на рис. 35, г.
Гамма-распределение непосредственным образом связано с пока зательным. Обозначим через tK время от начала эксплуатации до воз никновения k -то отказа, учитывая, что время между отказами рас пределено показательно с параметром X. Тогда случайная величина tK имеет распределение (3.27).
Таким образом, гамма-распределение является естественным обоб щением схемы мгновенных повреждений, когда для отказа необхо димо накопление k элементарных повреждений. Здесь следует отме тить, что рассматриваемое распределение относится к периоду уста новившегося износа, который характеризуется [16] постоянством скорости, рельефа шероховатости и постепенным ростом зазора без изменения физической картины изнашивания.
Гамма-распределение использовано, например, в работе [45] при исследовании надежности отдельных деталей судовых двигателей внутреннего сгорания.
Рассмотрим гамма-распределение еще с одной стороны. Пусть некоторая система состоит из k элементов, причем система работо способна, пока исправен хотя бы один элемент системы (примером
125
могут служить резервированные вспомогательные механизмы СЭУ). Тогда плотность распределения времени безотказной работы системы а (t) представляет собой свертку плотностей соответствующих рас пределений применительно к элементам системы [79]:
а (0 = [а, (*)]** = «1 (0 * [а, (01(* -1)* = |
К (О]2* * К (0](* -2)* = • • •, |
где |
о° |
|
|
[О/ (О!2* = [«1 (01 * 1а2(01 = |
J <h (т) {t — Т) dx. |
— СО |
Если время безотказной работы каждого из элементов распреде
лено показательно с параметром X, то |
|
|
t |
|
|
[at (£)]2* = X2 j е-*-те-(*-х>>■dx = |
ХЧе~~и . |
|
о |
|
|
Аналогично |
|
|
и по индукции |
|
|
a(t) = iai( t ) r = ^ - e ^ |
. |
(3.28) |
Сравнивая выражения (3.26) и (3.28), можно видеть, что время до отказа системы (время до отказа k элементов) имеет гамма-распре деление.
Таким образом, время до отказа k элементов, если для каждого из них справедлив экспоненциальный закон надежности, имеет
гамма-распределение |
с параметрами %и k. |
Из изложенного, |
в частности, можно сделать вывод о том, что |
если экспоненциальный закон надежности справедлив для каждого из элементов резервированной системы, то он не действует в мас штабе всей системы в целом.
Будем теперь рассматривать процесс механического износа как процесс возникновения элементарных повреждений. При достаточно большом числе таких повреждений наступает отказ элемента. Отсюда становится понятным, что время работы до отказа по причинам, свя занным с механическим износом, хорошо согласуется с гамма-рас пределением. С ростом числа элементарных повреждений (напри мер, процесс общей коррозии) гамма-распределение становится близким к нормальному.
Другим распределением, широко используемым в теории надеж ности, является распределение Вейбулла.
Пусть имеется система, состоящая из большого числа элементов, причем отказ любого из них является отказом всей системы. Пусть также каждый из элементов имеет гамма-распределение времени ра
боты до отказа |
с параметрами, мало изменяющимися от элемента |
|
к элементу. В |
этом случае время безотказной |
работы системы до |
рого согласуется с распределением Вейбулла |
[35, 79]. |
126
Частота отказов с случае распределения Вейбулла имеет вид
|
|
a(t) = lk tkle~Mk |
0). |
(3.29) |
|
Здесь X и k — параметры распределения. |
= 1 распределение Вей |
||||
Из выражения (3.29) следует, что при k |
|||||
булла |
тождественно |
показательному распределению. |
безотказ |
||
' В |
рассматриваемом случае |
основные характеристики |
|||
ности имеют вид: |
|
t |
|
л |
|
|
t |
|
|||
|
Р ( t ) = 1—аJ(т)d x = 1—X k | т*-1е ~ и к d x = e ~ u k ', |
|
|||
|
о |
|
о |
|
|
|
|
*’V) = T § - = xktk- 1; |
(3.30) |
||
|
СО |
СО |
р |
|
|
|
T = \ P { t ) d t = \ |
е ~ м к d t = - |
|
|
|
|
0 |
0 |
|
я Т |
|
Характер изменения во времени этих величин для различных k изображен на рис. 35, г.
Распределение Вейбулла связано с показательным. Пусть вели чина х имеет распределение Вейбулла (3.29), тогда случайная вели чина хк [35] распределена показательно с параметром АЛ, что также свидетельствует о том, что при k = 1 распределение Вейбулла тождественно показательному.
Рассмотрим теперь распределение Вейбулла с других позиций. Пусть имеется система, состоящая из п элементов, соединенных, в смысле надежности, последовательно. Это значит, что отказ лю бого из элементов приводит к отказу всей системы. Предполагая, что время безотказной работы каждого из элементов распределено по одному и тому же закону с функцией распределения Q (t) и плот ностью a (t), можно утверждать, что время до отказа системы опре деляется законом распределения наименьшей порядковой стати
стики выборки объема п, т. е. |
[81 ] |
Qi (t) = 1 — [ 1 - Q(OP; |
a, (t) = n [1 - Q (OP-1 a (t). |
Можно показать, что для большого п при достаточно общих условиях для многих распределений, включая нормальное (исклю чение составляет показательное распределение), наименьшая по рядковая статистика имеет распределение одного из двух типов (учи тываются лишь распределения, полезные в теории надежности) [79]:
= |
«.<<>= ! - « - “ *• |
(3.31) |
Сравнение выражений (3.29) и (3.31) показывает, что второе из распределений (3.13) является распределением Вейбулла. Отсюда становится понятным, что время работы до отказа системы, состоя щей из большого числа последовательно соединенных элементов, во многих случаях хорошо согласуется с законом Вейбулла.
127
Отметим, что если экспоненциальный закон надежности справед лив для каждого из элементов, то он справедлив и для всей после довательной системы.
Распределение (3.31) называется распределением Гумбеля и ис пользуется в ряде работ для анализа надежности трубопроводов забортной воды СЭУ, отказы которых вызываются язвенной корро зией. Такое использование становится правомерным, если учесть, что отказ трубопровода' вызывается язвами наибольшей глубины, т. е. справедлива изложенная выше модель слабейшего звена.
Время безотказной работы согласуется с логарифмически-нор- мальным законом, если логарифм этого времени имеет нормальное распределение. Рассматриваемый закон широко используется [16] при обработке эксплуатационных и опытных данных по усталостной долговечности металлов, их длительной прочности и т. п. В част
ности, |
логарифмически-нормальное |
распределение используется |
||||||||||
в работе |
[41 ] для анализа отказов крепежных соединений. |
|||||||||||
Частота |
отказов, |
т. е. |
плотность распределения случайной вели |
|||||||||
чины т, |
в рассматриваемом случае определяется выражением |
|||||||||||
|
|
|
|
|
|
|
|
(In 2-д ) 2 |
|
|
||
|
|
|
|
a(t) = |
|
1 |
2 d2 |
|
(3,32) |
|||
|
|
|
|
ot V2n |
|
|
|
|
||||
Здесь а и р , |
— параметры распределения. |
|
|
|
|
|||||||
Количественные характеристики безотказности в случае лога- |
||||||||||||
рифмически-нормального |
распределения |
имеют |
вид |
|||||||||
|
|
|
|
|
|
|
|
I |
|
|
1 ( 1пт-ц\ |
|
Р(0 = 1 |
|
|
|
|
1 |
Г-1 |
|
dx = |
||||
|
|
|
|
о]Х2 я J |
|
т |
|
|||||
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
In t—\l |
|
|
|
|
|
|
|
|
|
‘ - T |
W |
I |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 ( ln(—д \2 |
|||
|
|
|
|
a ( 0 |
|
_ |
1 |
|
2 |
' |
a |
' |
|
|
|
4 0 |
|
|
|
|
|
|
|||
|
|
|
P (t) |
|
OtV 2 it |
|
|
|
|
|
||
|
|
|
|
|
0 |
- |
4 |
4 |
' |
|||
|
|
|
|
|
|
|
|
|||||
|
|
|
” |
|
|
|
” |
|
(In(-H)a |
|||
|
|
|
|
|
|
|
|
|
|
|
2 o2 |
dt. |
|
|
|
T = } i a W * = x h i j e |
|
|
|
||||||
|
|
|
|
|
|
|
||||||
Характер изменения во времени приведенных выше характеристик |
||||||||||||
изображен |
на |
рис. |
35, д. |
|
|
|
|
|
|
|
Обратимся к физической картине возникновения отказов, согла сующихся с логарифмически нормальным распределением. Рассмот рим случай возникновения в работающем элементе усталостной тре
щины. Фиксируем некоторые моменты времени tx |
< t 2 < • |
• • < |
tn |
и обозначим через х х < 3 х г < . . . < « хп размеры |
трещины |
в |
эти |
128
моменты времени. Предполагаем, что элемент отказывает, когда тре щина достигнет размера хп, и что увеличение трещины на проме жутке (^_x, tt) (i — 1 , 2 , . . ., п) пропорционально уже достигну тому за предыдущее время размеру трещины х ^ ъ т. е. [79]
xi xl-1 = $ixi-l- |
(3.33) |
Здесь бг- — независимые положительные случайные величины с про извольными законами распределения;
г = 1 , 2 , . . ., п.
Обозначим через х 0 начальный размер трещины в элементе (нару шение структуры, пустоты, инородные включения и т. п.).
В соответствии с (3.33) можно записать
хп = (1 + 6 и) хп- 1 — О + 5Я) (1 -{- 6 „_j)• • • (1 + 8 i) х0. (3.34)
Таким образом, величина хп получена в виде произведения неза висимых, положительных случайных величин и, следовательно, ло гарифм хп равен сумме логарифмов сомножителей. Но тогда, со гласно центральной предельной теореме, In хп имеет асимптотически нормальное распределение, т. е. величина хп распределена по лога рифмически нормальному закону с плотностью
а ( 0 = |
|
(In t—fX)2 |
1 |
2а2 |
|
|
|
|
|
to V 2л |
|
Ранее было отмечено, что если In т имеет нормальное распределение, то случайная величина т согласуется с Логарифмически-нормальным законом. Отметим здесь также, что аналогичная связь существует между распределениями Гумбеля и Вейбулла [81]. Так, если In т имеет распределения Гумбеля, то случайная величина т распределена по закону Вейбулла.
Распределение Рэлея используется в основном для элементов, имеющих ярко выраженный эффект старения. Частота отказов в этом случае определяется выражением
где о — параметр |
распределения. |
|
||
Количественные характеристики безотказности в случае закона |
||||
Рэлея имеют |
вид: |
|
|
i l |
|
t |
|
t |
|
|
|
2о*_ п _ |
||
P(t) = 1 — J а (т) dt — 1 — |
J |
re 2° 2 dx = 1 — J e~x dx = e 2° 2; |
||
0 |
|
|
0 |
0 |
|
|
|
00 |
CO t2 |
= |
= |
T = f j P ( t ) d t = j e- ^ d t = y j a . |
||
|
|
|
0 |
0 |
Выражение для X (t) показывает, что интенсивность отказов про порциональна времени работы. Характер изменения во времени харак теристик безотказности для закона Рэлея изображен на рис. 35, д.
9 о. Р. Смирнов |
129 |