
Техническая диагностика ПС / КОЧЕРГА_ВГ_УП_НАДЕЖНОСТЬ_ТЕПЛОВОЗОВ
.PDF
|
|
Таблица 4.1 |
Рекомендуемые законы распределения |
||
|
|
|
Объект контроля |
Наиболее характерный вид |
Закон распределения |
|
и (или) причина отказа |
|
Тепловоз, агрегат |
Смешанный |
Экспоненциальный |
Деталь |
Износ |
Нормальный |
|
Старение |
Нормальный |
|
Усталостное разрушение |
Нормальный |
|
Контактная усталость (подшипники) |
Вейбулла–Гнеденко |
|
Разрушение от перегрузки |
Экспоненциальный |
|
Задиры |
Экспоненциальный |
|
Прогары |
Вейбулла–Гнеденко |
|
Коррозия |
Нормальный |
|
Смешанный |
Экспоненциальный |
4.2. Экспоненциальный закон распределения
Экспоненциальный закон достаточно точно описывает надежность узлов при внезапных отказах, имеющих случайный характер. Попытки применить его для других типов и случаев отказов, особенно постепенных, вызванных износом и изменением физико-химических свойств элементов показали его недостаточную приемлемость. Вероятность безотказной работы определя-
ется по уравнению
P(l) = e−λl . |
(4.1) |
Зависимости изменения вероятности безотказной работы от пробега для экспоненциального закона распределения показана на рис. 4.1.
В свою очередь средняя наработка до отказа определяется как
L = |
1 |
. |
(4.2) |
|
|||
1 |
λ |
|
|
|
|
Дисперсия наработки до отказа (на отказ) определяется по уравнению
Д(l) = |
1 |
. |
(4.3) |
|
|||
|
λ2 |
|
Среднеквадратичное отклонение наработки до отказа (на отказ)
σ(l) = |
|
= |
1 |
= L . |
(4.4) |
|
Д(l) |
||||||
|
||||||
|
|
|
λ |
1 |
|
|
|
|
|
|
|
Таким образом, равенство σ(l) = L1 – характерный признак экспоненци- ального распределения.
41

Необходимо отметить, что величина интенсивности отказа при экспо- ненциальном распределении является величиной постоянной и не зави-
сящей от времени или наработки объекта l = const.
Тогда вероятность безотказной работы объекта на интервале пробега l=L1 при экспоненциальном законе распределения составит
− L1 |
= e−1 0,368. |
|
P(l) = e−λL1 = e L1 |
(4.5) |
Таким образом, средняя наработка до отказа это наработка, в течение которой вероятность безотказной работы достигнет величины 0,368.
Частота отказов для экспоненциального закона распределения будет
a(l) = l × e−λl . |
(4.6) |
Зависимости изменения частоты отказов и интенсивности отказов от про- бега для экспоненциального закона распределения показаны на рис. 4.2.
P(l) 1
0,8
0,6
0,4
0,2
0 |
0,5 |
1 |
l/L1 |
|
Рис. 4.1. Функция P(l)
для экспоненциального закона распределения
Рис. 4.2. Функции a(l) и λ(l)
для экспоненциального закона распределения
Пример. В результате испытания десяти топливных насосов высокого давления получены наработки их до отказа: 400, 440, 500, 600, 670, 700, 800, 1200, 1600, 1800 ч. Предполагая, что наработка до отказа топливных насосов подчиняется экспоненциальному закону распределения, необхо-
димо оценить величину интенсивности отказов λ, а также рассчитать ве- роятность безотказной работы за первые 500 ч и вероятность отказа в промежутке времени между 800 и 900 ч работы дизеля.
Во-первых, определим величину средней наработки топливных насосов до отказа по уравнению (2.6):
L = |
1 ål |
= 400 + 440 + 500 + 600 + 670 + 700 + 800 +1200+1600+1800 = 871 ч. |
|||
|
10 |
|
|
|
|
1 |
|
1 i |
|
|
|
N |
10 |
|
42

Затем рассчитываем величину интенсивности отказов l по уравнению
(4.2):
l = |
1 |
= |
1 |
= 1,15 ×10 |
-3 |
1 |
. |
L |
871 |
|
ч |
||||
|
|
|
|
|
|||
|
1 |
|
|
|
|
|
|
Величина вероятности безотказной работы топливных насосов при на- работке 500 ч по формуле (4.1) составит
P(500) = e-ll = e-1,15×10−3×500 = 0,563.
Вероятность отказа в промежутке между 800 и 900 ч работы насосов имеет вид:
Q(900) - Q(800) = [1 - e-ll2 ]- [1 - e-ll1 ]=
=[1 - e-1,15×10−3×900 ]- [1 - e-1,15×10−3×800 ]= 0,645 - 0,601 = 0,044.
4.3.Нормальный закон распределения
Нормальный закон распределения хорошо описывает отказы, возни- кающие в элементах конструкции тепловозов в процессе изнашивания, ус- талости и коррозии деталей, работающих на номинальных (расчетных) режимах нагружения.
Вероятность безотказной работы при нормальном законе распределе-
ния определяется по уравнению
æ |
l - L1 |
ö |
|
|
ç |
÷ |
(4.7) |
||
|
||||
P(l) = 1- F ç |
sТ |
÷, |
||
è |
ø |
|
где sТ – дисперсия; F – интеграл Лапласа.
Зависимость изменения вероятности безотказной работы от наработки показана на рис. 4.3.
Частота отказов и интенсивность отказов для нормального распреде- ления рассчитываются по уравнениям:
|
|
|
1 |
|
|
|
æ |
(l - L |
)2 ö |
|
|
|
|||
|
a(l) = |
|
|
|
|
|
|
ç |
|
1 |
÷ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
2 |
|
|
exp ç- |
|
2 |
÷; |
|
|
(4.8) |
||
|
|
|
|
|
2psТ |
è |
|
2 sТ |
ø |
|
|
|
|||
|
1 |
æ |
|
|
|
(l - L )2 ö |
é |
æ l - L |
öù-1 |
|
|||||
l(l) = |
|
|
|
ç |
- |
1 |
÷ |
|
ç |
1 |
÷ |
(4.9) |
|||
|
|
|
|
||||||||||||
2 |
|
exp ç |
2 |
÷ |
ê1- Fç |
sТ |
÷ú . |
||||||||
|
|
2psТ |
è |
|
|
|
2 sТ |
ø |
ë |
è |
øû |
|
Зависимости частоты отказов и интенсивности отказов для нормального закона распределения показаны на рис. 4.4.
43

P(l)
l
Рис. 4.3. Функции P(l) для нормального
закона распределения
Рис. 4.4. Функции a(l) и λ(l)
для нормального закона распределения
4.4. Закон распределения Вейбулла–Гнеденко
Закон распределения Вейбулла–Гнеденко получил широкое распростра- нение и используется применительно к системам, состоящим из рядов эле- ментов, соединенных последовательно с точки зрения обеспечения безотказ- ности системы. Например, системы, обслуживающие дизель-генераторную установку: смазки, охлаждения, питания топливом, воздухом и т. д.
Распределение Вейбулла–Гнеденко обычно задается двухпараметри-
ческой функцией частоты отказов
= b æ l öb−1 a(l) ç ÷
a è a ø
é |
æ l öb ù |
|
|||
expê- ç |
|
÷ |
ú, |
(4.10) |
|
|
|||||
ê |
è a ø |
ú |
|
||
ë |
|
|
|
û |
|
где b – параметр формы; a – параметр масштаба распределения. Иногда
параметр масштаба задается в виде a = a−b . В этом случае выражение (4.10) преобразуется к виду:
a(l) = |
b ×lb−1 |
é |
æ l öb ù |
|
|||
ab |
expê- ç |
|
÷ |
ú. |
(4.11) |
||
|
|||||||
|
ê |
è a ø |
ú |
|
|||
|
|
ë |
|
|
|
û |
|
Зависимость частоты отказов представлена на рис. 4.5.
Функция распределения вероятности безотказной работы определяет- ся следующим образом:
é |
æ l öb ù |
|
|||
P(l) = exp ê- ç |
|
÷ |
ú. |
(4.12) |
|
|
|||||
ê |
è a ø |
ú |
|
||
ë |
|
|
|
û |
|
Вид функции вероятности безотказной работы для распределения Вейбулла–Гнеденко показан на рис. 4.6.
44

a(l) |
b<1 |
P(l) |
|||
|
|
|
|
||
|
b>1 |
|
b<1 |
||
|
|
b=1 |
|||
|
b=1 |
|
|
|
|
|
|
|
|
b>1 |
|
|
|
|
|
|
|
|
l |
|
l |
Рис. 4.5. Функция a(l) для распределения Рис. 4.6. Функция P(l) для распределения
|
|
Вейбулла–Гнеденко |
|
|
|
|
|
|
|
|
Вейбулла–Гнеденко |
|
|
|
||
Интенсивность отказов |
|
|
|
|
|
|
öb−1 |
|
|
|
||||||
|
|
|
|
l(l) = |
b |
æ l |
|
|
|
|||||||
|
|
|
|
|
|
ç |
|
÷ . |
(4.13) |
|||||||
|
|
|
|
a |
|
|||||||||||
Наработка до отказа |
|
|
è a |
ø |
|
|
|
|
||||||||
|
|
æ |
|
|
1 ö |
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
L1 |
= aГ ç1 + |
|
|
÷, |
(4.14) |
|||||||
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
è |
|
|
|
b ø |
|
|
|
||
где Г – табулированная гам- |
|
|
|
|
|
|
|
|
|
|
|
|
||||
ма-функция. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Зависимость интенсивности |
|
|
|
|
|
|
|
|
|
|
|
|
||||
отказов |
для |
распределения |
|
|
|
|
|
|
|
|
|
|
|
|
||
Вейбулла–Гнеденко представ- |
|
|
|
|
|
|
|
|
|
|
|
|
||||
лена на рис. 4.7. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Пример. |
Время |
простоя |
|
|
|
|
|
|
|
|
|
|
|
|
||
тепловозов в неплановых ре- |
|
|
|
|
|
|
|
|
|
|
|
|
||||
монтах по вине вспомогатель- |
|
|
|
|
|
|
|
|
|
|
|
|
||||
ного |
оборудования |
подчиня- |
|
|
|
|
|
|
|
|
|
|
|
|
||
ется |
закону |
распределения |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Вейбулла–Гнеденко с пара- |
|
|
|
|
|
|
|
|
|
|
|
|
||||
метрами b = 2 и a = 46. Тре- |
Рис. 4.7. Функция λ(l) для распределения |
|||||||||||||||
буется |
определить |
вероят- |
|
|
|
|
|
|
|
|
Вейбулла–Гнеденко |
|
|
|
||
ность |
выхода тепловозов из |
|
|
|
|
|
|
|
|
|
|
|
|
неплановых ремонтов после 24 ч простоя и время простоя, в течение ко- торого работоспособность будет восстановлена с вероятностью 0,95.
Из уравнения (4.12) найдем вероятность восстановления работоспо-
собности локомотива после простоя его в депо в течение суток
é |
æ |
24 |
ö2 |
ù |
= 1- 0,76 = 0,24. |
Q(24) = 1- P(24) = 1- exp ê- ç |
46 |
÷ |
ú |
||
ê |
è |
ø |
ú |
|
|
ë |
|
|
|
û |
|
45

Для определения времени восстановления работоспособности локо- мотива с заданной величиной доверительной вероятности также исполь- зуем выражение (4.12)
é |
æ |
|
x ö |
2 ù |
|
|
é |
æ |
x ö2 |
ù |
|
|||
Q(x) =1 - exp ê- ç |
|
|
|
÷ |
ú |
= 0,95; |
exp ê- ç |
|
÷ |
ú |
= 0,05; |
|||
|
|
|
|
|||||||||||
ê |
è 46 ø |
ú |
|
|
ê |
è |
46 ø |
ú |
|
|||||
ë |
|
|
|
|
|
û |
|
|
ë |
|
|
|
û |
|
|
æ |
|
x ö2 |
|
|
|
|
|
|
|
|
|||
- 2,99 = ç |
|
|
|
÷ ; |
x = 79,5 ч » 3,3 сут. |
|
||||||||
46 |
|
|||||||||||||
|
è |
ø |
|
|
|
|
|
|
|
|
4.5. Закон распределения Рэлея
Закон распределения Рэлея используется в основном для анализа ра- боты элементов, имеющих ярко выраженный эффект старения (элементы электрооборудования, различного рода уплотнения, шайбы, прокладки, из- готовленные из резиновых или синтетических материалов). Функция час-
тоты отказов в этом случае определяется выражением
a(l) = |
2l |
é |
æ |
l ö2 |
ù |
|
||
|
|
exp ê- ç |
|
÷ |
ú, |
(4.15) |
||
S |
2 |
|
||||||
|
ê |
è |
S ø |
ú |
|
|||
|
|
|
ë |
|
|
|
û |
|
где S – параметр распределения.
Количественные характеристики безотказности в случае закона рас- пределения Рэлея имеют вид:
– интенсивность отказов
|
|
|
|
l(l) = |
2l |
, |
|
|
|
|
|
|
|
|
(4.16) |
|
|
|
|
|
|
|
|
|
|
|
|
||||
– вероятность безотказной работы |
S 2 |
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|||||
P(l), |
|
|
|
|
é |
æ l |
ö |
2 ù |
|
||||||
a(l), |
|
a(l) |
|
|
P(l) = expê- ç |
|
|
|
÷ |
ú, |
(4.17) |
||||
|
|
|
|
||||||||||||
(l) |
|
|
|
|
ê |
è S |
ø |
ú |
|
||||||
|
|
|
|
ë |
|
|
|
|
|
|
û |
|
|||
|
|
P(l) |
|
наработка до первого отказа |
|||||||||||
|
|
|
|
|
|
|
S |
|
|
|
|
|
|||
|
|
|
|
|
|
L = |
|
p |
. |
|
|
(4.18) |
|||
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
1 |
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(l) |
|
Характер изменения во времени |
|||||||||||
|
|
|
|
функций безотказности |
для |
закона |
|||||||||
|
|
l |
|
распределения Рэлея |
|
показан на |
|||||||||
|
Рис. 4.8. Функции P(l), a(l), λ(l) |
рис. 4.8. |
|
|
|
|
|
|
|
|
для распределения Рэлея
46

Пример. Для величины наработки 120 тыс. км требуется определить ве- роятность безотказной работы, интенсивность отказов и среднюю наработку до первого отказа катушки электромагнитного контактора. Известно, что нара- ботки контакторов до отказа по параметрам старения изоляции катушек мож- но описать функцией распределения Рэлея с параметром S = 260 тыс. км.
По выражениям (4.16)–4.18) находим:
|
|
|
3 |
|
é |
|
æ |
120 ×103 |
ö2 |
ù |
|
|
||||
|
P(120 ×10 |
) = exp |
ê |
- |
ç |
|
|
÷ |
|
ú |
= 0,81; |
|
||||
|
|
|
3 |
|
||||||||||||
|
|
ê |
ç |
260 ×10 |
÷ |
|
ú |
|
||||||||
|
|
|
|
|
|
è |
|
ø |
|
|
|
|||||
|
|
|
|
|
ë |
|
|
|
|
|
|
|
û |
|
|
|
|
2×260×103 |
|
|
1 |
|
|
|
|
|
× 260 ×103 |
|
|||||
l(120×103 ) = |
= 3,5×103 |
|
; |
L = |
|
|
p |
= 230 ×103 км. |
||||||||
|
|
|
|
|
|
|
|
|||||||||
(260×103 )2 |
|
|
|
|
км |
|
1 |
|
|
|
|
|
2 |
|
||
|
|
|
|
|
|
|
|
|
|
|
|
4.6. Определение закона распределения
При исследовании надежности тепловоза, его узлов и деталей возни- кает необходимость определения вида и аналитического выражения зако- на распределения длительности работы (наработки) до отказа. Имея ста- тистические данные об отказах узла, можно определить теоретический за- кон распределения пробега до отказа. В теории надежности используется
несколько методов определения законов распределения и расчета оценок их параметров по статистическим данным: метод моментов; метод разде- ляющих разбиений; графический метод.
Наиболее приемлемым для условий локомотивных депо является метод моментов, так как он применяется в случае, если известны реализации l1, l2, l3 …ln наработки до отказа n узлов. Для получения вариационного ряда их располагают в порядке возрастания. Далее диапазон пробегов ln – l1 делят на k интервалов l (желательно одинаковой величины)
Dl @ |
lmax − lmin |
, |
(4.19) |
|
|||
1 + 3,3lg n |
|
|
где lmax, lmin – соответственно наибольшее и наименьшее значение в ва- риационном ряду.
Затем подсчитывают число случаев попадания пробега до отказа n в каждый интервал, введя условную единицу xi:
xi = |
li |
, |
(4.20) |
|
Dl |
||||
|
|
|
где li – значение пробега для середины i-го интервала.
47

По формуле |
ni |
|
|
|
|
|
W (x ) = |
, |
|
|
(4.21) |
||
i |
n |
|
|
|
|
|
|
|
|
|
|
|
|
определяется частность отказов в каждом интервале пробега. |
|
|
||||
|
|
|
По |
форме |
построен- |
|
|
|
|
ной гистограммы частнос- |
|||
|
|
|
ти (рис. 4.9) можно выд- |
|||
|
|
|
винуть гипотезу о том, ка- |
|||
|
|
|
кой из известных теорети- |
|||
|
|
|
ческих законов распреде- |
|||
|
|
|
ления |
ближе |
всего |
по |
|
|
|
форме подходит к стати- |
|||
|
|
|
стическим данным. |
|
||
|
|
|
Параметры, |
входящие |
||
|
|
|
в аналитическое выраже- |
|||
|
|
|
ние закона распределения |
|||
Рис. 4.9. Пример гистограммы частности |
|
|
определяются |
приравни- |
||
|
|
ванием |
моментов теоре- |
|||
|
|
|
||||
тического распределения к соответствующим моментам статистического |
||||||
распределения. Для однопараметрических законов (экспоненциальный, Рэ- |
||||||
лея) приравнивается математическое ожидание; для двухпараметрических |
||||||
законов (нормальный, Вейбулла–Гнеденко), кроме этого, приравниваются |
||||||
дисперсии или среднеквадратические отклонения. |
|
|
|
|
||
Из полученных таким образом уравнений находят неизвестные пара- |
||||||
метры предполагаемых законов распределения. Точность соответствия |
||||||
теоретического закона распределения статистическому распределению |
||||||
проверяют критериями согласия (Колмогорова, Пирсона и др.). |
|
|
||||
Наиболее распространенным является критерий χ2 Пирсона, опреде- |
||||||
ляемый по зависимости |
ni − npi ) |
2 |
|
|
|
|
k |
|
|
|
|||
χ 2 = å ( |
, |
|
(4.22) |
|||
i=1 |
npi |
|
|
|
|
|
где i = 1, 2, …, k – интервалы группирования случайной величины χ в выбо- |
||||||
рочной совокупности объемом n значений; |
ni – частота попадания слу- |
|||||
чайной величины в i-й интервал; pi – теоретическая вероятность попада- |
||||||
ния в i-й интервал. |
|
|
|
|
|
|
При расчете критерия χ2 используют эмпирические (опытные) выбороч- |
||||||
ные данные ni и теоретические частоты npi, найденные по теоретической |
||||||
формуле закона. Количество интервалов группирования определяет число |
||||||
степеней свободы r = k – c – 1, где c – число неизвестных параметров за- |
48
кона. По величине и уровню значимости q находят величины χ2q . И если гипотеза верна, то при достаточно большом n
P{χ2 > χq2 }= q. |
(4.23) |
Следовательно, если найденная по опытным данным |
величина |
{χ2 > χ2q } при выбранном уровне значимости, то гипотеза о законе распре-
деления не противоречит опытным данным. Для критерия χ2 установлены весьма жесткие условия применения: n ≥ 200 и ni ≥ 18.
Метод разделяющих разбиений удобно применять в тех случаях, когда n элементов, включенных в работу, подвергаются проверке через некото- рый пробег Ln (например, до планового ремонта). Если за этот пробег от- казало m1 элементов, то отношение m1/n есть значение функции распреде- ления F(l=Ln).
При двухпараметрическом законе распределения, кроме этого, необхо- димо еще знать число отказавших элементов m2 за пробег 2Ln. Значения функции распределения вычисляются по формулам:
|
F(L |
n |
) = |
m1 |
, |
|
|
|
|
(4.24) |
|||
|
n |
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
||
F(2L |
) = |
m |
|
|
m |
2 |
æ |
|
m |
ö |
|
||
1 |
+ |
|
|
|
ç1 |
- |
1 |
÷. |
(4.25) |
||||
|
|
|
|
|
|
||||||||
n |
n |
|
|
n - m1 è |
|
n ø |
|
||||||
|
|
|
|
|
|
Приравнивая выражение функции теоретического распределения при l1=Ln и l2=2Ln к соответственным значениям, вычисленным по формулам (4.24) и (4.25), можно определить параметры закона распределения.
Графические методы определения законов распределения и их пара- метров, как правило, предполагают применение так называемых вероят- ностных бумаг. Вероятностные бумаги разрабатываются для каждого вида закона распределения. Они представляют собой прямоугольную систему координат с особыми масштабами по обеим осям для изображения графи- ка функции распределения случайной величины прямой линией.
КОНТРОЛЬНЫЕ ВОПРОСЫ
1.Назовите основные законы распределения случайных величин.
2.Назовите законы распределения случайных величин и соответст- вующие им наиболее характерные виды или причины отказов.
3.Расчет показателей надежности при экспоненциальном законе рас- пределения.
49
4.Чему равна вероятность безотказной работы объекта на интервале равном средней наработке до отказа при экспоненциальном законе рас- пределения?
5.Назовите последовательность действий при определении закона распределения случайной величины.
5. РАСЧЕТ НАДЕЖНОСТИ СЛОЖНЫХ СИСТЕМ
5.1. Общие сведения
Как отмечалось в предыдущих разделах, тепловоз является сложной системой, состоящей из элементов, соединенных друг с другом в опреде- ленной последовательности.
Элементом называется часть системы, не имеющая самостоятельного
эксплуатационного назначения и выполняющая в системе ограниченные функции. В свою очередь, системой называется совокупность совместно действующих устройств, обеспечивающих выполнение определенных практических задач.
В зависимости от постановки задач можно один и тот же объект рассмат- ривать как элемент или как систему. Например, при оценке надежности теп- ловоза его рассматривают как систему, состоящую из основных элементов: дизеля, тягового генератора, выпрямителя, охлаждающих систем, тяговых электродвигателей, вспомогательных устройств и систем регулирования и управления. В том случае, когда возникает необходимость в определении по- казателей надежности, например, тягового электродвигателя или генератора, которые по отношению к конструкции тепловоза являются его элементами, их рассматривают уже как систему, состоящую из основных элементов: якоря, остова, щеточного аппарата, полюсов и т. д. При оценке надежности якоря тя- гового электродвигателя его рассматривают тоже как систему, состоящую из вала, якорной втулки, сердечника, коллектора, секций обмотки и т. п.
По характеру влияния отказов на выполнение определенной задачи в процессе применения по назначению системы подразделяют на простые и сложные. Простая система при отказе элементов либо полностью прекра- щает выполнение своих функций, либо продолжает выполнять свои функ- ции в полном объеме, если отказавший элемент заменяется резервным. Сложные системы в результате наличия у них функциональной избыточ-
ности при отказе отдельных элементов и подсистем могут выполнять свои функции при некотором снижении характеристик эффективности. Отказ сложной системы определяется как событие, обусловленное выходом ха- рактеристик эффективности за нижний допустимый предел, и связанное с ним частичное или полное невыполнение задачи. Элементы в сложных системах могут иметь два типа соединения: резервное и основное.
50