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

Техническая диагностика ПС / КОЧЕРГА_ВГ_УП_НАДЕЖНОСТЬ_ТЕПЛОВОЗОВ

.PDF
Скачиваний:
179
Добавлен:
09.06.2015
Размер:
605.28 Кб
Скачать

 

 

Таблица 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

= e1 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