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

книги из ГПНТБ / Болошин, Н. Н. Надежность работы технологических узлов и оборудования обогатительных фабрик

.pdf
Скачиваний:
29
Добавлен:
19.10.2023
Размер:
7.14 Mб
Скачать

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

Отказы

также разделяются на п ол ные , которые приводят к

полной потере работоспособности объекта, системы, и

ч а с т и ч ­

ные, при

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

системы.

Полная классификация отказов приведена в работах

[5, 22, 23.

39].

 

 

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

Главным показателем надежности является б е з о т к а з н о с т ь , которая определяет свойство системы непрерывно сохранять ра­ ботоспособность в течение некоторой наработки без вынужден­ ных перерывов; безотказность количественно оценивается вероятно­ стью безотказной работы либо косвенными вероятностными показа­ телями — интенсивностью отказов, наработкой на отказ и другими показателями. Понятие безотказности часто используется для оп­ ределения «надежности» в узком смысле.

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

ского обслуживания и

ремонтов.

Д о л г о в е ч н о с т ь ю

называется свойство системы сохранять

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

Для количественной оценки надежности работы систем в тео­ рии надежности применяют следующие характеристики *.

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

Н а р а

б о т к о й

м е ж д у о т к а з а м и называется продолжи­

тельность

времени

между двумя последовательно возникшими от­

казами. Среднее значение (математическое ожидание) наработки

между отказами называется н а р а б о т к о й

на о т к а з .

К о э ф ф и ц и е н т о м т е х н и ч е с к о г о

и с п о л ь з о в а н и я

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

Для повышения надежности работы системы (например, кон­

* Приводятся характеристики, используемые для определения надежности восстанавливаемых систем.

10

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

иненагруженный резерв, когда резервные элементы подготовлены

квключению (в таком режиме работают дробилки и насосы).

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

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

Для подробного ознакомления с теорией вероятностей может быть рекомендован курс теории вероятностей Е. С. Вентцель [43].§

§ 2. Некоторые сведения из теории вероятностей

Понятия и определения надежности тесным образом связаны с

понятиями с о б ы т и я и с

л у ч а й и о й вел и ч и н ы [43].

Под с о б ы т и е м

в

теории вероятностей понимают всякий

факт, который в результате опыта может произойти или не про­ изойти. Возможность появления некоторого события (например, отказа) изменяется числом Р, называемым вероятностью этого со­ бытия ( 0 ^ Я ^ 1 ) .

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

заполняют некоторый промежуток,

называются н е п р е р ы в -

н ы м и.

 

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

отказавших технологических узлов;

примерами непрерывных слу­

чайных величин являются время безотказной работы или время восстановления.

II

Пусть дискретная случайная величина Т может принять в ре­ зультате опыта значения t\, t2, t3, .. tn. Каждое из них возможно, поэтому существует некоторая вероятность того, что в результате опыта величина Т примет данное значение.

Отношение числа опытов т, в результате которых случайная

величина

Т приняла

значение

к общему числу произведенных

опытов п

называется

ч а с т о т о й

появления события T — tj. Ча­

стота — является случайной величиной и меняется в зависимости

П

от количества произведенных опытов. Но при большом числе опы­ тов она имеет тенденцию к стабилизации около некоторого значе­ ния pi, называемого вероятностью события 7 = /,-. Поэтому послед­ нюю на' практике вычисляют по формуле

р,. =

Р ( Т = ^ ) = ^ .

(1)

_

а

 

(

 

Сумма вероятностей всех возможных значении случайной вели­

чины

 

 

 

£ р ; = П

(2)

 

i=i

 

Эта суммарная вероятность распределена определенным обра­ зом между отдельными значениями. С вероятностной точки зре­ ния случайная величина полностью описывается заданием этого распределения, т. е. точным указанием, какой вероятностью обла­ дает каждое из событий T = ti при любом i от 1 до п.

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

Простейшей формой задания закона распределения дискретной случайной величины является ряд распределения, в котором пере­ числены возможные ее значения и соответствующие им вероятно­

сти. Графическое изображение ряда

распределения называется

м н о г о у г о л ь н и ком р а с п р е д е л е н и я.

Для примера ниже приведен ряд

распределения случайной ве­

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

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

 

 

Количество от-

0

1

2

3

4

5 6

7

8

9

10

казавших кону­

 

 

 

 

 

 

 

 

 

 

сов, шт.

2,9

12,0

24,0

27,0

20,0

8,0

5,0

1,0 0,08

0,02

0,01

Вероятность от-

каза, 96

 

 

 

 

 

 

 

 

 

 

12

Закон распределения количественно выражается в двух фор­ мах:

1. Как для непрерывной, так и для дискретной случайной ве­ личины удобно пользоваться вероятностью события T<t, где I — текущая переменная. Вероятность такого события, зависящая от значения /, называется ф у н к ц и е й р а с п р е д е л е н и я слу­ чайной величины Т:

Р {t) — Р(Т <

t).

(3)

Такая функция распределения

иногда

называется

и и т е г р а л ь-

н о й. *

есть Неубывающая функция, а зна­

Функция распределения Р (t)

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

Р(— оо) = 0 и Р (+ оо ) = 1.

2. Для непрерывной случайной величины наиболее часто при­ меняется производная функция распределения — п л о т н о с т ь р а с п р е д е л е н и я случайной величины Т :

p(t) = P'(t).

(4)

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

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

Вероятность нахождения

величины Т в интервале от

а до р

выражается через плотность распределения:

 

Р (а <

t < Р) = J’ р (/) dt.

(5)

 

а

 

Плотность распределения является неотрицательной функцией; площадь под очерчивающей ее кривой (кривой распределения) и осью абсцисс равна единице:

- р С О

 

 

J p ( 0 * = l -

(6)

---00

 

 

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

определению (4),

выражается

через плотность распределения следующим образом:

 

P{t)=-- \

p{f)dt.

(7)

* В практике обогащения такую функцию называют кумулятивной.

13

Геометрически P(t) есть площадь под кривой распределения и осью абсцисс, лежащая левее ординаты, проходящей через точку t. Та­ ким образом, для описания случайной величины используются:

для дискретной — функция распределения и ряд распределения (графически — многоугольник распределения);

для непрерывной — функция распределения и плотность распре­ деления (графически— кривая распределения) [5, 43].

Рис. I. Кумулятивные кривые и кривые плотности распределения времени без­ отказной работы и времени восстановления подбункерного узла на Тырныаузскоп фабрике:

а — кумулятивная кривая времени безотказной работы; б — кривая плотности распределения времени безотказном работы; а — кумулятивная кривая времени восстановления; г — кривая плотности распределения времени восстановления

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

.Наиболее распространенными числовыми характеристиками слу­ чайной величины являются математическое ожидание и дисперсия ((или среднее квадратическое отклонение). .

Если случайная величина Т принимает значения tь U, t3,.. ., tn соответственно с вероятностями р\, рг, Рз,- ■ рп, то сумма произ­

14

ведений возможных значении случайной величины на вероятности

этих значений называется

м а т е м а т и ч е с к и м

о ж и д а н и е м :

П

+ О Э

tp {t) dt.

(8)

М (Г) = т , = 2

tipГ, M ( D = mt = J

f — I

— CO

 

 

Первая формула справедлива для дискретных случайных величин, вторая — для непрерывных.

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

Для оценки рассеивания случайной величины вокруг матема­ тического ожидания используются дисперсия и среднее квадрати­

ческое отклонение.

случайной величины называется математиче­

Д и с п е р с и е й

ское ожидание квадрата отклонения случайной величины от своего математического ожидания.

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

формуле

 

 

 

 

 

D (Т) =

М[Т — М (Г)]2 =

£ & -

mtf pt,

(9)

 

 

 

1=1

 

 

а для непрерывной случайной величины по формуле

 

D{T)=--M[T — М(Т)]2=

-{-С О

(ti — mt)2p(f)dt.

(10)

J

 

 

— о о

 

 

 

Дисперсия имеет

размерность

квадрата

случайной

величины,

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

Для этого определяют с р е д н е е к в а д р а т и ч е с к о е

от-

к л о и е н и е случайной величины:

 

о(Т) = а( = у Щ Т ) .

(11)

При исследовании надежности работы технических устройств важнейшим является вопрос об установлении закона распределе­ ния случайных величин: времени безотказной работы и времени восстановления — и соответствия их теоретическим законам.

Из теоретических законов при исследовании надежности рабо­ ты технологических узлов и оборудования обогатительных фабрик представляют интерес следующие.

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

_

t_

 

Р (0 = е

г ,

(12)

15

где Т — математическое ожидание случайной величины.

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

ствия большого

количества

случайных

действующих факторов.

Н о р м а л ь н ы й з а ко н .

Случайная

величина

нормально рас­

пределена, если

ее плотность

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

вид

 

 

(t-ту

 

 

p(t) = -----1= - е 2а-

(13)

 

а

V

 

 

Постоянные Т и о больше 0 и могут быть любыми: Т — математи­ ческое ожидание, а —- среднее квадратическое отклонение случай­ ной величины.

3 а ко н В е й б у л л а. Случайная величина имеет распределе­ ние Вейбулла, если функция распределения описывается уравне­ нием

Р(() = е

(14)

где Т — математическое ожидание случайной

величины; п — пока­

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

Наряду с перечисленными законами при исследовании надеж­ ности используются и другие распределения: равномерное, гамма, Релея и др. [5, 43, 49]. На рис. 2 приведены кривые теоретических функций распределения, используемых при изучении надежности.

Исследование надежности работы обслуживаемых восстанав­ ливаемых объектов производится на основе изучения отказов и последующих восстановлений, которые представляют собой потоки случайных событий **.

На длительном промежутке времени работа технологических узлов и оборудования может рассматриваться как процесс смены состояний. Узел, проработав случайное время //, отказывает и пребывает в состоянии отказа время t". После восстановления си­ стема работает время t2' и опять отказывает на время to" и т. д. (рис. 3).

Длительности промежутков времени работы t'„ и отказа 1п" представляют собой случайные величины:

* В теории надежности уравнение этого закона было предложено Вейбуллом в 1951 г. [5]. В 1934 г. аналогичная формула была предложена Розиным и Раммлером для характеристики распределения отдельных классов крупности сы­

пучих материалов [44—46]. Исследование надежности работы

технологических

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

лежит в преде­

лах 0,5—1,2.

 

** Изучение потоков случайных событий производится на основе теории мас­ сового обслуживания.

(6

периоды работы (,/ независимы друг от друга и от tn", распре­ делены по закону с математическим ожиданием Г, и дисперсией

0i2:

 

Рх(0 = р (tn > t ) \ ТХ= М ( Q ; а? = D {Q ;

(15).

периоды отказов tn" также независимы друг от друга и от

рас­

пределены по закону с математическим ожиданием Т2 и диспер­ сией ог2:

Р2(/) = р {irl > t) ; Т.г = М {t„ ); а 2= D{t„).

( 16>

Рис. 2. Кумулятивные кривые теоретических функций распределения:

а _ равномерное;

б — экспоненциальное; в — экспоненциально-степенное (Розина—

 

Раммлера, Вейбулла); г •— нормальное

Поток отказов,

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

пренебрежимо мало, превращается в поток отказов без восстанов­ ления, такой поток определяется распределением времени безот­

казной работы.

 

 

а

Количественной характери­

 

стикой потока отказов и вос­

 

становлений является интен­

 

сивность потока, которая рав­

 

на среднему числу отказов вос­

 

станавливаемой

системы

за

 

единицу времени, взятому в

 

данный момент

времени.

.Эта

Рис. 3. Поток отказов:

характеристика

является

диф­

а — с конечным време«ва1...васстаноцлеинп;

ференциальной

характеристи­

б — с пренебрежимо м!лым времен'Ш’ТГОг-^-— -

становления1 ОС. пу-5л«-:ц|. р.

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

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

Поток отказов характеризуется случайной величиной V(t), рав­ ной количеству отказов, происшедших за определенное время. Ве­ роятность того, что количество отказов V ( t ) ^ n , определяют сле­ дующим образом:

■Р\v (t) >п] = р [tn<t]=--.p \t[ + h + . ■ . + t'n<t\= Qn(t), (17)

где Qn(0 представляет собой закон распределения времени без­ отказной работы.

Тогда

Рп(0 =•■Р [V(/) -- п) = Qn{t) - Qn+, (0

(18)

и, в частности, вероятность, что не наступит ни одного отказа (ве­ роятность безотказной работы на интервале t ) :

^0 (0 = Я [1/ (/) =

0] == 1 — Q (0-

 

(19)

Функция восстановления Я (t) определяется следующим обра­

зом:

 

 

 

 

 

со

со

 

Q,.,.!_1(/)] =

со

Qn (t). (20)

П {() = М [V (0] = 2

пРп(0 = V „ [Q„ (() -

2

п = \

л= 1

 

 

/1=1

 

Среднее число отказов на участке

to)

определяется

разностью

Я (to)H(ty). Если интенсивность потока отказов представляет со­ бой дифференциальную характеристику, а функция восстановле­

ния интегральную,

то они

связаны соотношениями:

 

 

 

со

 

(21)

h (t) = Я' (t); Я(/ )=\ ' h{t)dt,

 

но так как

 

b

 

 

 

 

 

 

СО

Qn(0, a

Qn (t) = f n0f), ТО h (t) =

со

(22)

Я (0 = V

v f n (t).

/1=1

 

 

n=1

 

Величины Pn (t), H(t) и h(t) используются при решении практиче­ ских задач надежности, связанных с определением вероятности ■.обеспечения 100%-ной производительности, планируемого количе­ ства ремонтов, определения количества запасного оборудования и т. п.

Для экспоненциального закона вышеприведенные характери­

стики определяются по выражениям

 

 

= Р [Е(0 = n] =

( - i - )

е ~

(23)

Я(*) = ^11- ;

т =

11

(24)

 

18

Для других законов эти характеристики определяются слож­ ными формулами и не выражаются в конечом виде. На практикеИнтерес представляет определение характеристик потока иа боль­ шом отрезке времени, т. е. предельное значение этих характери­ стик. Для этих условий в теории надежности предложены простые приближенные формулы, справедливые для любого закона распре­ деления [5, 26, 51]. Для большого промежутка времени среднее число отказов, приходящееся на единицу времени, близко к вели­ чине, обратной среднему времени жизни элемента:

lim ЛЮ. — Пт h(t) — —— .

(25)

t - ю э L i ->-со

Т х

 

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

Для интервала времени t\ + t2 , если время безотказной работы распределено непрерывно, справедливо

lim [ft +

4) — Я & )] =

 

(26)

/-►СО

 

 

11

 

 

Если время безотказной работы Г,

имеет дисперсию а2, то

Пт

Я (0 — -ф-

а2

1

 

 

a r f _ Т ’

 

 

t -*-СО l_

I I

 

 

тогда для большого времени верны формулы

 

 

 

а2

1

— -

1 < Я ( 0 <

(27)

 

2т\

2

 

Тг

 

П

 

На большом интервале времени количество отказов

V(/)

сходится

с величиной М[Е(/)]

и определяется:

 

 

 

 

 

Ti

D { V { t ) ] ^ ^ - .

 

(28)

 

 

 

Т ъ,

 

 

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

 

 

а V 1

< v ® <

ф

и

a-/

t

(29)

 

 

7’3/2

т 3/=

 

 

Т

 

где

И а определяется по таблицам

квантилей

нормального зако-

 

~

 

 

 

 

 

 

 

на *.

 

 

 

 

 

 

 

М.,

* Б. В.

Г н е д е н к о и др.

Математические

методы

в

теории надежности.

«Наука»,

1965.

 

 

 

 

 

 

19

Соседние файлы в папке книги из ГПНТБ