
книги из ГПНТБ / Болошин, Н. Н. Надежность работы технологических узлов и оборудования обогатительных фабрик
.pdfционного свойства, носят случайный характер появления во вре мени, не зависят от длительности эксплуатации. Постепенные отка зы зависят от длительности эксплуатации, проявляются в постепен ном количественном изменении состояния объекта и связаны с яв лением износа или старения. Характер отказов предопределяет за кон распределения отказов во времени.
Отказы |
также разделяются на п ол ные , которые приводят к |
|
полной потере работоспособности объекта, системы, и |
ч а с т и ч |
|
ные, при |
которых возможно частичное использование |
системы. |
Полная классификация отказов приведена в работах |
[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 2л |
|
|
Постоянные Т и о больше 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