Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Бескровный Н.Т. Экономика и оптимизация надежности и ремонта горношахтного оборудования

.pdf
Скачиваний:
8
Добавлен:
24.10.2023
Размер:
11.13 Mб
Скачать

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

При наличии значительного числа наблюдений обработку исход­ ных данных удобно проводить с помощью ЭВМ по специальным про­ граммам [6]. Это позволяет сократить время на проведение исследо­ ваний.

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

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

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

По данным литературы [21], [47], наработка на отказ многих невосстанавливаемых изделий — элементов средств автоматики и радиоэлектронной аппаратуры, а также тканых сит для углеобога­ тительных фабрик, радиально-поршневых насосов НП-120 в гидра­ влических механизмах подачи [1], магнитных логических элементов типа ЭПМ и основной части транзисторных элементов серии ЭТ аппаратуры горной автоматики [51], подчиняется экспоненциаль­ ному распределению. Это распределение используется при рассмо­ трении внезапных отказов в тех случаях, когда явления износа п старения настолько слабо выражены, что ими мояшо пренебречь.

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

Наработка на отказ у многих невосстанавливаемых изделий имеет распределение Вейбулла и логарифмически нормальное. К этим изделиям относятся, например, подшипники качения, рабочие колеса шламовых насосов, некоторые типы электронных ламп, некоторые полупроводниковые приборы, некоторые изделия, у которых отказ наступает вследствие усталостного разрушения. В работе [2] отме­ чено, что наработка на отказ основных узлов шахтных вагонедок хорошо согласуется с распределением Вейбулла.

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

10

Рис. 1. Характеристики кривой нормаль­ ного распределения:
- среднеквадратическое отклонение

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

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

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

Такое распределение более всего удовлетворяет слу­ чаю, когда узлы и детали подвержены как постепен­ ным, так и внезапным от­ казам [36]. В реальных условиях весьма вероятны и более сложные компози­ ции и суперпозиции раз­ личных законов.

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

дальнейшем будем рекомендовать таблицы Я. Б. Шора и Ф. И. Кузь­ мина [72].

Нормальное распределение. Плотность вероятности нормального распределения находится по уравнению

л

(х-а)*

( 1)

=

,

а У 2 л

 

 

где а — математическое ожидание;

о2 — дисперсия.

 

Изучение нормальной кривой показывает, что она определяется некоторыми основными характеристиками, приведенными на рис. 1. На этом рисунке видно, что 99,73% площади под кривой лежит в пределах шести квадратических отклонений, т. е. по три средних квадратических отклонения по каждую сторону от среднего значе­ ния. Аналогично 95,49% площади лежит в пределах четырех сред­ них квадратических отклонений и 68,26% в пределах двух подобных

11

отклонений. Кроме того, можно заметить, что среднеквадратическое отклонение равно х/6 размаха. Это очень хорошая оценка, так как на рисунке видно, что за трехсигмовыми пределами остается по 0,135% площади на каждую сторону.

В практических приложениях удобно пользоваться нормирован­

ным (а = 1) и центрированным

= 0) распределением плотности

вероятности

 

 

 

 

( 2)

Из уравнения (2) следует, что при отрицательных значениях

надо пользоваться уравнением

 

 

Ф о ( я ) =

Ф о ( - я ) .

( 3 )

Для центрированного и нормированного распределения плот­ ность вероятности табулирована и приведена в [72, табл. 1].

Переход к / (х) на основании уравнений (1), (2) можно осущест­

вить по формуле

 

/ ( * ) = 4 ф о ( - ^ ) .

(4)

Пример. Пусть заданы а — 50; т = 10; х = 40. Требуется опре­ делить плотность вероятности. В соответствии с уравнениями (4), (3), и [72, табл. 1] находим

/<4°>= і <Р.

) -

^

/< -1 ) = JS t (i)-

i s ■0,2420 = 0,0242.

Функция нормального распределения находится но

уравнению

 

F(x)

 

1

Г _ Іѵ-а>*

dx.

(5)

 

 

ст V 2п

J е га’

 

 

 

 

 

 

Для нормированного и центрированного распределения имеем функцию

 

л

2 dx.

 

F0(x)

 

( )

J*

 

6

 

 

 

Эта функция табулирована и приведена в [72, табл. 1. 21. Из уравнения (6) следует, что

F0( - x ) = i - F 0(x).

(7)

Переход к F (х) на основании уравнений [см. (5), (6)1 можно осуществить по формуле

F(x) = F0( ^ ) .

(8)

12

Если наработка я до отказа распределена по нормальному закону, то вероятность отсутствия отказа на промежутке от 0 до ж нахо­ дится по уравнению

у

 

OD

 

*(*)“ j f ( x ) ä x ~ l - F ( x ) = F0( ^ L ) .

(9)

Пример. Пусть наработка до отказа имеет нормальное распреде­ ление с параметрами а — 500 ч и о — 100 ч.

Определить вероятность безотказной работы для наработки х = = 200 ч и X = 700 ч.

По уравнению (9) находим

Р (20О)= ^ ( 50° - о200 ) = ^ о(3).

На основании [72, табл. 1.2] находим

Р(3) = 0,998650.

Аналогично находим

P (700) = ^(500^700 ) = JP ( _ 2)I

На основании уравнения (7) и [72, табл. 1.2] находим

Р (700) = 1- F0(2) = 1 - 0,97725 = 0,02275.

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

Интенсивность отказов находится по уравнению

табулированная функция [72, табл.. 1.4].

Пример. Пусть наработка

до отказа имеет нормальное распределе­

ние с параметрами а = 800

ч,

а = 150 ч.

Определить интенсивность

отказа для

х — 100 ч, х = 1000 ч.

По уравнению (10) и [72,

табл. 1.4]

находим

Ч 100>==®г і^ш1)

0,00000637 =0,4246- ІО-7І .

150

Аналогично находим

800—1000

) 150 ^

1.770

=0,118-IO'7- .

150

' 150

H

13

Экспоненциальное распределение. Распределение случайной поло­ жительной величины называется экспоненциальным, если его плот­ ность вероятности имеет вид

f ( t ) = h r u ,

(12)

где X — постоянная (параметр распределения).

Функция экспоненциального распределения находится по уравне-

ншо F{t) = 1—е-м. (13)

Если t — наработка, то вероятность безотказной работы на интер­ вале от нуля до t будет

P(*) = e-w = exp(-to).

(14)

Знак ехр читается экспонента и обозначает число е в степени, приведенное за ним в скобках.

Вычисление по уравнению (14) удобно производить по функции

ехр (—я) [72, табл. 2.1].

 

Пример. Допустим, интенсивность отказов реле X = 0,0006

.

Найти вероятность безотказной работы за 300 ч.

 

По уравнению (14) и [72, табл. 2.1] находим

 

Р (300) = ехр (-0,0006 • 300) = ехр (-0,18) = 0,8353.

 

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

M(t) = Tcp = - f ,

(15)

» ч о - і -

(16)

 

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

По формуле (15) находим

Гср“ - 7^

= 1667 ч.

 

0,0006

 

По формуле (16) находим

 

б2:

1

=2777778 ч2.

0,00062

Коэффициент вариации для экспоненциального распределения

равен единице

а (О

 

 

V(t):

1.

(17)

M(t)

Распределение Вейбулла. Случайная положительная величина t имеет распределение Вейбулла, если плотность вероятности ее запи­ сывается в виде

/ ^ =т ( т ) й_1ехр [ - ( т ) Т ’

(18)

где а и Ъ — положительные постоянные.

14

Плотность вероятности, умноженная на постоянную а, табули­ рована для различных значений Ъ и t/a [72, табл. 3.3].

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

Пример. Найти плотность распределения в точке t = 200 ч, если наработка изделия до отказа распределена по закону Вейбулла с параметрами а = 900 ч, Ъ = 2,5.

Определяем \

По [72, табл. 3.3] для этого значения и Ъ = 2,5 находим а/ (t) = 0,8380.

Отсюда определяем

/(« ) -- 2^ - = 0,0°0931.

Функция распределения Вейбулла записывается в виде

F(t) = l - e x p

(19)

На основании этого уравнения можно записать выражение для определения вероятности отсутствия отказа на промежутке от 0 до t.

P(t) = l ^(0=*exp [ — (-7 )* ].

(20)

Значения вероятности Р (t) в функции от t/a и Ъ табулированы

и приведены в

[72, табл. 3.1].

 

примера

найти наработку

Пример. В

условиях

предыдущего

до отказа, отвечающую вероятности 0,3679.

 

1.

По [72, табл. 3.3] для t

= 2,5 и Р

(t)

= 0,3679, находим t/a =

Отсюда t — \а =

1 -900 =

900 ч.

 

[(18), (20)] находится

по

Интенсивность

отказов

на основании

уравнению

 

 

 

 

 

 

 

 

 

 

- К

г

) “ 1-

<21>

Для математического ожидания и среднего квадратического откло­ нения величины t, распределенной по закону Вейбулла, справед­ ливы формулы

 

tcp = аКъ,

 

(22)

где

o(t) = aCB,

(23)

 

 

 

 

-

х , = г ( і +

| ) ,

 

(24)

С ,= ] / г ( і +

| ) - К

| .

(25)

Знак Г означает

гамма-функцию.

 

 

15

Коэффициент вариации определяется по уравнению

(

S7-

<26>

Значения величины К в, Св и ѵ в функции от Ъ приведены в [72, табл. 3.5].

Гамма-распределение. Случайная положительная величина под­ чиняется гамма-распределению, если ее плотность вероятности имеет вид

/(*) = X ^ g J e x p ( - ;U ),

■ (27)

где к и т — положительные постоянные, причем пъ — целое число. Это уравнение при х = 2kt можно записать в следующем виде

f(*>=W w *m-1ехР ( - f ) •

(28)

Функция / (х) зависит только от одного параметра т, табулиро­ вана и приведена в [72, табл. 4.1].

Переход к f (t) осуществляется по формуле

f{t) = 2kt(x).

(29)

Вероятность безотказной работы

 

СО

 

P(t)= J fit)dt

(30)

t

 

может быть вычислена по уравнению

 

P(t) = Po(x).

(31)

Значения Р 0 (х) в зависимости от я и т приведены в [72, табл. 4.2]. При помощи функций (29) и (31) можно вычислить интенсив­

ность отказов

X(t) =

n t )

2М(х)

(32)

Pit)

Po ix)

Пример. Определить интенсивность отказов в точке t = 500 ч, если наработка изделия имеет гамма-распределение с параметрами

т = 3,

%=

0,004.

2•0,004*500 = 4.

По [72, табл. 4.1 и

Находим

X = 2kt =

табл. 4.2] для

х = 4, т = 3 определяем

 

 

 

 

/(ж) = 0,135, Р0 = 0,677.

По

формуле

(32) находим

 

 

 

 

Ч*) =

2 • 0,004 • 0,135 0,0016 ^ .

 

 

 

0,677

Ч

16

Математическое ожидание величины t, ее дисперсия и коэффи­ циент вариации находятся по уравнениям

«ср = ^ ( 0

= х *

(33)

 

 

 

(34)

a t)

_■ i_

(35)

v ( t )

V m

►cp

 

Логарифмически нормальное распределение. Положительная слу­ чайная величина у имеет логарифмически нормальное распределение, если ее логарифм распределен нормально.

Для плотности вероятности распределения имеем

f(y)-

фо (

^

Фо (

,

(36)

где М = 0,4343 — коэффициент перехода от натуральных к деся­ тичным логарифмам; ln yQ, lg yQ— математическое ожидание лога­ рифма случайной величины у\ аг, а 2 — среднеквадратическое отклонение логарифма случайной величины у, соответственно, для натуральных и десятичных логарифмов.

При этом имеют место соотношения

\gyQ= M \ n y 0,

 

 

 

(37)

а2=

Маи

 

 

 

(38)

Функция распределения величины у

записывается в виде

 

Іиу—Inyo \ _ я

(

lgy—Igyp ^

 

F ( y ) ^ o ( "Ln- ~ lnyo

) = F 0

(-

а 2

/ -

< 3 9 >

а ,

/

 

Значения функций ф0 (х) и FQ(х) находятся по таблицам для нор­

мального распределения [72, табл. 1.1

и

табл. 1.2].

 

 

Вероятность безотказной работы на протяжении наработки у

находится л о

уравнению

 

 

 

 

 

 

 

 

 

Р ( у ) - 1-

Л

ln у ~ In i/о

\ _

4

р

(

lg yУ-—lglgyp

\

(40).

( 1д‘,~ 1™

) °

1-

 

°

\

°2

) '

Интенсивность

отказов определяется

по

уравнению

 

1 /,л _

1

4

Г ln y -ln уо 1

М j

/

lgy-lgyo

\

(41)

где функция f x (у)

 

 

 

 

 

 

 

 

----J’

находится по

[72,

табл. 1.4]. '

 

 

Для математического ожидания ут дисперсии а\ и коэффициента

вариации ѵу справедливы равенства

 

 

 

 

 

 

 

 

Ут = Уо ехр

т

:уоехр

(-

 

)■

 

(42)

 

 

 

 

 

 

 

 

V2М2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(43)

1>„=1/'вхр(о;)-1='|/ехр

 

 

 

1.

 

(44)

2 Заказ 353

I

17

на отказ

Распределение Релея.

Случайная положительная величина t имеет распределение Релея, если плотность вероятности ее записывается в виде

где о2дисперсия. /(*)=-£•e*p[—ä f].

<45>

Функция этого распределения

имеет

вид

 

P(i) = l_ e x p

( - £

) .

(46)

Бели t — наработка, то вероятность безотказной работы на интер­ вале от 0 до t будет

Р(і) = ѳхр ( - ■ £ ? ) ■

(47)

Математическое ожидание и коэффициент вариации определя­

ются по формулам

 

M(t) = Tcp = o ] / - J ,

(48)

у = 1 т = Ѵ і г = О-79788-

(49)

Практически модель закона распределения обычно выбирается на основании некоторых общих соображений. Однако правильный выбор исходной теоретической модели закона распределения в зна­ чительной степени определяет необходимый объем статистических исследований, требуемых для оценки показателей надежности с заданной достоверностью. Например, когда было доказано, что закон распределения времени отказа полупроводниковых приборов лучше согласуется с законом Вейбулла, чем с экспоненциальным законом [49], стало возможным сокращать объемы выборочных испытаний в 4—6 раз.

Характер зависимостей показателей надежности во времени для рассмотренных законов распределения приведен на рис. 2. В случае нормального, Вейбулла (при 6 > 2), гамма (при т > 1) и логариф­ мически нормального распределений интенсивность отказов с тече­ нием времени ускоренно возрастает. В случае распределения Релея и Вейбулла (при Ъ — 2) интенсивность отказов возрастает пропор­ ционально. времени наработки. В случае распределения Вейбулла (при 1 <; Ъ<; 2) интенсивность отказов возрастает с замедлениемЭкспоненциальное, Вейбулла (при Ъ = 1) и гамма (при т = 1) распределения имеют постоянную интенсивность отказов.

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

2 *