Н а д е ж н о с т ь / poz115
.pdf
Элемент 3. Распределение Рэлея:
P (t) = e −λt 2 |
= e−8 10−8 t 2 |
3 |
|
Элемент 4. Экспоненциальное распределение:
P4 (t) = e −λt = e −0,0002 t
Элемент 5. Усечённое нормальное распределение:
|
|
|
|
|
|
|
|
t |
− m0 |
|
|
|
t |
− 2000 |
|||||
|
|
|
∞ |
|
( x −m0 ) |
2 |
|
|
|
|
|
|
|
0,5 −Ф0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
c |
|
|
|
|
0,5 −Ф0 |
σ0 |
|
|
|
|
900 |
|
||||||
|
|
|
− |
|
|
|
|
|
|
|
|
|
|
|
|
||||
P5 (t) = |
|
∫e |
2σ02 |
|
dx = |
|
= |
|
|
||||||||||
σ0 |
2π |
|
|
|
|
|
m0 |
|
|
0,5 +Ф |
|
2000 |
|
||||||
|
t |
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
0,5 +Ф0 |
σ0 |
|
|
|
|
0 |
|
900 |
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Табулируя эти функции от 0 до 2000 часов с шагом 100 часов, получим таблицу 1.7.
Таблица 1.7 – Вероятность безотказной работы элементов
t, час |
P1(t) |
P2(t) |
P3(t) |
P4(t) |
P5(t) |
Pc(t) |
0 |
1 |
1 |
1 |
1 |
1 |
1 |
100 |
0,980199 |
1 |
0,995696 |
0,9992 |
0,996918 |
0,972194 |
200 |
0,960789 |
0,999994 |
0,990256 |
0,996805 |
0,987730 |
0,936745 |
300 |
0,941765 |
0,999917 |
0,983464 |
0,992826 |
0,972604 |
0,894281 |
400 |
0,923116 |
0,999532 |
0,975087 |
0,987282 |
0,951817 |
0,845456 |
500 |
0,904837 |
0,998321 |
0,964883 |
0,980199 |
0,925741 |
0,790895 |
600 |
0,886920 |
0,995466 |
0,952605 |
0,971611 |
0,894839 |
0,731242 |
700 |
0,869358 |
0,989932 |
0,938013 |
0,961558 |
0,859646 |
0,667280 |
800 |
0,852144 |
0,980612 |
0,920884 |
0,950089 |
0,820755 |
0,600058 |
900 |
0,835270 |
0,966491 |
0,901022 |
0,937255 |
0,778801 |
0,530939 |
1000 |
0,818731 |
0,946799 |
0,878275 |
0,923116 |
0,734444 |
0,461577 |
1100 |
0,802519 |
0,921097 |
0,852542 |
0,907738 |
0,688351 |
0,393774 |
1200 |
0,786628 |
0,889326 |
0,823788 |
0,891188 |
0,641180 |
0,329303 |
1300 |
0,771052 |
0,851793 |
0,792053 |
0,873541 |
0,593567 |
0,269727 |
1400 |
0,755784 |
0,809123 |
0,757456 |
0,854875 |
0,546108 |
0,216247 |
1500 |
0,740818 |
0,762184 |
0,720202 |
0,835270 |
0,499352 |
0,169613 |
1600 |
0,726149 |
0,712001 |
0,680578 |
0,814810 |
0,453789 |
0,130105 |
1700 |
0,711770 |
0,659674 |
0,638951 |
0,793581 |
0,409845 |
0,097577 |
1800 |
0,697676 |
0,606303 |
0,595754 |
0,771669 |
0,367879 |
0,071540 |
1900 |
0,683861 |
0,552922 |
0,551479 |
0,749162 |
0,328179 |
0,051268 |
2000 |
0,670320 |
0,500461 |
0,506654 |
0,726149 |
0,290960 |
0,035911 |
|
|
|
11 |
|
|
|
В последнюю колонку записаны значения вероятностей безотказной работы системы, которые определяются произведением вероятностей безотказной работы элементов:
Pс(t) = P1(t) P2(t) P3(t) P4(t) P5(t)
На рисунке 1.4 показаны графики функций Pi(t), i=1,2,3,4,5, соответствующих вероятностям безотказной работы элементов. Номера графиков соответствуют номерам элементов. На рисунке 1.5 изображён график вероятности безотказной работы системы Pс(t).
Рисунок 1.4 – Вероятность безотказной работы элементов
Рисунок 1.5 – Вероятность безотказной работы системы
Из графиков видно различное поведение вероятностей безотказной работы элементов. Скорость убывания вероятностей зависит от вида и параметров закона распределения. В нашем случае медленнее всего убывает P(t) для
12
экспоненциального распределения и распределения Рэлея, т.е. при большом времени работы наиболее надёжными оказываются третий и четвёртый элементы системы.
Вычислим среднее время безотказной работы системы:
T1 = ∫∞ Pc (t)dt = ∫∞ P1 (t)P2 (t)P3 (t)P4 (t)P5 (t)dt
0 0
по формуле Симпсона:
T |
= |
h |
|
1 |
+ n−1 |
((3 |
+ (−1)k )P (kh)P (kh)P (kh)P (kh)P (kh)) |
, |
|||||
|
|||||||||||||
1 |
3 |
|
|
∑ |
|
1 |
2 |
3 |
4 |
5 |
|
|
|
|
|
|
|
k =1 |
|
|
|
|
|
|
|
|
|
где n – число точек;
h – шаг интегрирования, выбираемый из условия обеспечения требуемой точности.
Расчёты показывают, что для данных таблицы 1.7 Т1=976,3 час.
На рисунке 1.6 изображены графики интенсивностей отказов элементов. Кривая 4, соответствующая экспоненциальному закону, параллельна оси времени, т.к. имеет постоянную интенсивность отказа. Все остальные кривые интенсивностей отказов являются возрастающими функциями времени.
На рисунке 1.7 показан график интенсивности отказа системы, равной сумме интенсивностей отказов её элементов:
λc (t) = λ1 (t) + λ2 (t) + λ3 (t) + λ4 (t) + λ5 (t)
Рисунок 1.6 – Интенсивность отказов элементов
13
Рисунок 1.7 – Интенсивность отказа системы
Интенсивность отказа системы также является возрастающей функцией времени, что говорит о том, что система является стареющей, а закон распределения времени до её отказа не экспоненциальный.
Вычислим плотности распределения вероятностей времени безотказной работы элементов.
Элемент 1. Распределение Вейбулла:
|
αt |
α−1 |
|
|
t α |
|
|
2t |
|
|
t |
2 |
|
f (t) = |
e |
|
|
|
= |
|
e |
|
|
|
|||
|
|
|
1800 |
|
|||||||||
|
|
β |
|
||||||||||
|
|
|
|
− |
|
|
|
|
|
|
− |
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
||
βα |
|
|
|
|
|
|
18002 |
|
|
|
|
||
Элемент 2. Гамма-распределение:
f |
|
(t) = |
t |
α−1 |
e− |
t |
t |
6 |
e − |
t |
||
|
|
β |
= |
|
|
|||||||
|
|
|
300 |
|||||||||
2 |
βα Г(α) |
3007 Г(7) |
||||||||||
|
|
|
|
|
|
|
||||||
Элемент 3. Распределение Рэлея:
f3 (t) = 2λte −λt 2 = 2 8 10−8 te −8 10−8 t 2
Элемент 4. Экспоненциальное распределение:
f4 (t) = λ e −λt = 0,0002e −0,0002t
Элемента 5. Усечённое нормальное распределение:
|
|
|
c |
e − |
( t −m0 )2 |
|
1 |
|
|
|
e − |
( t −2000)2 |
f5 |
(t) = |
|
2σ02 = |
|
|
|
|
|||||
σ0 |
|
|
2000 |
2 9002 |
||||||||
|
|
2π |
|
|
|
+Ф |
|
|
|
|||
|
|
|
|
|
|
900 2π 0,5 |
0 |
|
|
|
|
|
|
|
|
|
|
|
900 |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
||
14
Табулируя плотности распределения от 0 до 2000 часов с шагом 100 часов, получим таблицу 1.8.
Таблица 1.8 – Плотности распределения времени безотказной работы элементов
t, час |
f1(t) |
f2(t) |
f3(t) |
f4(t) |
f5(t) |
0 |
0,0002 |
0 |
0,000038 |
0 |
0 |
100 |
0,000196 |
0 |
0,000048 |
0,000016 |
0,000062 |
200 |
0,000192 |
0 |
0,000061 |
0,000032 |
0,000122 |
300 |
0,000188 |
0,000002 |
0,000075 |
0,000048 |
0,000180 |
400 |
0,000185 |
0,000007 |
0,000092 |
0,000063 |
0,000235 |
500 |
0,000181 |
0,000019 |
0,000112 |
0,000078 |
0,000286 |
600 |
0,000177 |
0,00004 |
0,000134 |
0,000093 |
0,000331 |
700 |
0,000174 |
0,000072 |
0,000158 |
0,000108 |
0,000371 |
800 |
0,000170 |
0,000116 |
0,000185 |
0,000122 |
0,000405 |
900 |
0,000167 |
0,000168 |
0,000213 |
0,000135 |
0,000433 |
1000 |
0,000164 |
0,000227 |
0,000242 |
0,000148 |
0,000453 |
1100 |
0,000161 |
0,000288 |
0,000272 |
0,000160 |
0,000467 |
1200 |
0,000157 |
0,000347 |
0,000303 |
0,000171 |
0,000475 |
1300 |
0,000154 |
0,000402 |
0,000332 |
0,000182 |
0,000476 |
1400 |
0,000151 |
0,000450 |
0,000360 |
0,000191 |
0,000472 |
1500 |
0,000148 |
0,000487 |
0,000385 |
0,000200 |
0,000462 |
1600 |
0,000145 |
0,000514 |
0,000407 |
0,000209 |
0,000448 |
1700 |
0,000142 |
0,000530 |
0,000425 |
0,000216 |
0,000430 |
1800 |
0,000140 |
0,000535 |
0,000438 |
0,000222 |
0,000409 |
1900 |
0,000137 |
0,000531 |
0,000446 |
0,000228 |
0,000385 |
2000 |
0,000134 |
0,000517 |
0,000449 |
0,000232 |
0,000359 |
Графики, построенные по данным таблицы 1.8, представлены на рисунке 1.8.
Рисунок 1.8 – Плотности распределения времени до отказа элементов
Плотность распределения времени до отказа системы fc(t) изображена на рисунке 1.9. Для её изображения вычисления выполнялись по формуле:
fc (t) = λc (t)Pc (t)
15
Рисунок 1.9 – Плотность распределения времени до отказа системы
Из графика отчётливо видна неэкспоненциальность распределения времени до отказа нерезервированной системы, если законы распределения времени до отказа её элементов не являются экспоненциальными.
1.3 Задачи для самостоятельного решения
Задача 1.1. Техническая система состоит из n = 3 подсистем, которые могут отказать независимо друг от друга. Отказ каждой подсистемы приводит к отказу всей системы. Вероятность того, что в течение времени t первая подсистема проработает безотказно, равна 0,7; вторая – 0,9; третья – 0,8. Найти вероятность того, что в течение времени t система проработает безотказно. Найти вероятность отказа системы за время t.
Задача 1.2. Известно, что серийно выпускаемая деталь имеет экспоненциальное распределение времени до отказа с параметром λ=10-5 час-1. Деталь используется конструктором при разработке нового прибора. Назначенный ресурс прибора Тн=104 час.
Определить следующие показатели надёжности детали:
-вероятность отказа детали до момента Тн;
-вероятность того, что деталь безотказно проработает в течение времени Тн;
-вероятность того, что деталь безотказно проработает в интервале времени от
103 до 104 час.
Задача 1.3. Проектируется нерезервированная система, состоящая из элементов четырёх групп. Количество элементов каждой группы, а также интенсивности их отказов приведены в таблице 1.9.
16
Таблица 1.9 – Данные о числе элементов системы и интенсивности их отказов
Номер группы |
Число элементов |
Интенсивность отказа элемента, |
час-1 |
||
1 |
10 |
2·10-6 |
2 |
15 |
4·10-6 |
3 |
32 |
2,5·10-6 |
4 |
8 |
5·10-6 |
Определить:
-интенсивность отказа системы;
-среднее время безотказной работы;
- |
вероятность безотказной работы системы в течение времени t1 = 100 час, |
|
t2 = 1000 час и в интервале указанных наработок; |
- |
плотность распределения времени безотказной работы системы при наработке |
|
t2 = 1000 час. |
Задача 1.4. Система состоит из пяти элементов. Данные о их надёжности приведены в таблице 1.10.
Таблица 1.10 – Законы распределения времени до отказа элементов и их параметры
Вариант |
|
|
Элементы |
|
|
|
1 |
2 |
3 |
4 |
5 |
||
|
||||||
1 |
TN (390; 100) |
Г (9; 65) |
Exp (8·10-5) |
R (2·10-5) |
W (5; 200) |
|
2 |
R (1·10-5) |
W (4,5; 180) |
Г (8; 77) |
TN (400; 92) |
Exp (1·10-4) |
|
3 |
Г (10; 70) |
Exp (5·10-5) |
TN (375; 86) |
R (3·10-5) |
W (4,8; 190) |
|
4 |
TN (380; 100) |
R (1,6·10-5) |
W (7; 210) |
Exp (2·10-4) |
Г (9; 85) |
|
5 |
W (6; 195) |
TN (410; 95) |
Exp (2·10-5) |
Г (8; 75) |
R (2,5·10-5) |
Определить:
-вероятность безотказной работы системы;
-среднее время безотказной работы системы;
-интенсивность отказов системы;
-плотность распределения времени до отказа системы.
Решение представить в аналитическом виде, в виде графиков и таблиц.
Задача 1.5. Система состоит из пяти элементов с экспоненциальными законами распределения времени до отказа. Показателями их надёжности являются: P1(100) = 0,99; λ2 = 0,00001 час-1; Т3 = 8100 час, Т4 = 7860 час, λ5 = 0,000025 час-1.
Определить время t, в течение которого система будет исправна с вероятностью 0,92.
17
Задача 1.6. Система состоит из пяти элементов с постоянными интенсивностями отказов. Вероятность безотказной работы элементов в течение t
часов имеют следующие значения: P1(100) = 0,99; P2(200) = 0,97; P3(157) = 0,98; P4(350) = 0,95; P5(120) = 0,98.
Определить вероятность безотказной работы системы в течение 625 часов её функционирования, а также среднее время безотказной работы.
Задача 1.7. Время работы до отказа серийно выпускаемой детали распределено по нормальному закону с параметрами: m = 1000 час, σ = 250 час. Определить:
-вероятность того, что деталь проработает безотказно более 1200 часов;
-вероятность того, что наработка до отказа будет находиться в интервале [m−3σ, m+3σ];
-вероятность того, что, безотказно проработав до момента времени 1200 часов, деталь безотказно проработает и до 1500 часов.
Задача 1.8. Комплектующая деталь, используемая при изготовлении устройства, по данным поставщика имеет нормальное распределение времени до отказа с параметрами m = 4000 час, σ = 1000 час. Определить следующие показатели надёжности детали:
-наработку до отказа, соответствующую 90 % надёжности детали;
-вероятность того, что деталь имеет наработку, лежащую в интервале [2000; 3000];
-вероятность того, что деталь имеет наработку, большую, чем 4000 часов.
18
2 РАСЧЁТ ПОКАЗАТЕЛЕЙ НАДЁЖНОСТИ РЕЗЕРВИРОВАННЫХ НЕВОССТАНАВЛИВАЕМЫХ СИСТЕМ
2.1 Методы расчёта показателей надёжности
Критерии надёжности резервированных невосстанавливаемых систем те же, что и нерезервированных невосстанавливаемых систем.
Основными видами резервирования являются: общее постоянное, общее замещением, раздельное постоянное, раздельное замещением. Структурные схемы резервированных систем приведены на рисунке 2.1.
Приведём основные соотношения для показателей надёжности резервированных систем.
2.1.1 Общее резервирование с постоянно включенным резервом
Пусть Pi(t) – вероятность безотказной работы i-го элемента за время t, Qi(t) – вероятность отказа i-го элемента за время t, fi(t) – плотность распределения времени до отказа i-го элемента в момент времени t. Тогда вероятность безотказной работы, плотность распределения времени безотказной работы и интенсивность отказов системы с кратностью резервирования m определяются соотношениями:
m
Pc (t) = 1 − П(1 − Pi (t)) , (2.1)
i =0
m |
|
|
|
|
|
|
fc (t) = ∑(1 |
− P0 (t))...fi (t)...(1 − Pm (t)) , |
(2.2) |
||||
i =0 |
|
|
|
|
|
|
|
|
m |
m |
(t ) |
|
|
|
|
∑ f j ( t ) |
ПQ i |
|
||
λc ( t ) = |
j =0 |
i ≠ j |
|
, |
(2.3) |
|
|
m |
|
||||
|
|
1 − ПQ i (t ) |
|
|||
|
|
|
i =0 |
|
|
|
19
а – общее резервирование с постоянно включённым резервом; б – раздельное резервирование с постоянно включённым резервом; в – общее резервирование замещением; г – раздельное резервирование замещением
Рисунок 2.1 – Структурные схемы резервированных систем
20