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

Лабораторная работа № 3

.pdf
Скачиваний:
24
Добавлен:
10.06.2015
Размер:
593.58 Кб
Скачать

МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РОССИЙСКОЙ ФЕДЕРАЦИИ

ФЕДЕРАЛЬНОЕ АГЕНСТВО ПО ОБРАЗОВАНИЮ

К Г Э У

Государственное образовательное учреждение высшего профессионального образования

 

«КАЗАНСКИЙ ГОСУДАРСТВЕННЫЙ ЭНЕРГЕТИЧЕСКИЙ УНИВЕРСИТЕТ»

Кафедра «Промышленная теплоэнергетика установки и системы теплоснабжения»

Экз. № ______

ЛАБОРАТОРНАЯ РАБОТА

по учебной дисциплине " Надежность систем теплоэнергоснабжения

промышленных предприятий "

Лабораторная работа № 3

ОПРЕДЕЛЕНИЕ ЗАКОНА НАДЁЖНОСТИ НЕВОССТАНАВЛИ-

ВАЕМЫХ ТЕХНИЧЕСКИХ ОБЪЕКТОВ ПО ПОЛНОСТЬЮ ОПРЕДЕ-

ЛЁННОЙ ВЫБОРКЕ

Разработал: ассистент М.В. Акутин, д.т.н. Ю.В. Ваньков

Обсуждена на заседании кафедры протокол № _____

от "_____" ___________2012 г.

Казань - 2012 г.

1

Лабораторная работа №3. Определение закона надѐжности невосста-

навливаемых технических объектов по полностью определѐнной выборке

Учебные цели:

1.Исследовать эмпирический закон распределения случайной величины.

2.Освоение порядка построения гистограммы.

3.Изучить процедуру использования критерия согласия χ².

2

Наблюдения за эксплуатацией большого количества однотипных элементов в технических объектах различного назначения, работающих в примерно одинаковых условиях и на схожих по интенсивности нагрузках, показали, что количество отказов этих элементов различно в разные периоды наработки. Можно выделить три периода наработки, которые заметно отличаются количеством отказов и интенсивностью их проявления. На начальном этапе эксплуатации количество отказов постепенно снижается (приработка), затем следует период эксплуатации с примерно постоянным количеством отказов в одинаковые интервалы времени и при дальнейшем увеличении наработки происходит нарастание числа отказов вследствие быстрого накопления повреждений (например, вследствие износа). Эти периоды эксплуатации характеризуются определѐнным значением интенсивности отказов (t), показанной на рис.1.

Рис.1 Типичная - характеристика

Кривую на рис.1 принято называть - характеристикой. Как показано на рис.1 первый участок - характеристики отличается повышенным уровнем интенсивности отказов. В этот период (рис.1,I) происходит приработка составных элементов объекта, устранение мелких дефектов изготовления и сборки, выявление и устранение отступлений от технической документации. Для данного этапа эксплуатации объекта иногда используют термин «этап выжигания дефектов». При отработке конструкции и технологии изготовления, а также при совершенствовании выходного контроля интенсивность отказов на начальном этапе эксплуатации может быть снижена.

При нормальной эксплуатации (рис.1,II) интенсивность отказов практически не меняется с наработкой. Причинами отказов в этот период бывают, как правило, воздействие на объект неучтѐнных при проектировании факторов или эксплуатация объектов в условиях, действие которых трудно предусмотреть заранее.

3

Третий период эксплуатации (рис.1,III) характеризуется возрастанием интенсивности отказов. Это обычно связано с предельным накоплением повреждений в материалах основных деталей, а также с износом или старением конструкционных материалов. Этот период обычно связан с повышенными материальными затратами а поддержание работоспособности объектов. Основные из них идут на запасные части и на восстановление изношенных деталей.

В практических расчѐтах надѐжности - характеристики отдельных элементов имеют большое значение. Часто модели надѐжности сложных технических систем полностью создаются на использовании подобных характеристик. Во многих случаях решение вопросов обоснования сроков технического обслуживания, объѐмов регламентных работ невозможно без знания - характеристик комплектующих элементов. Эти характеристики в большинстве своѐм получаются на основе статистической обработки экспериментальных данных или результатов эксплуатации объектов – аналогов. Одним из важных моментов такой обработки является выявление теоретического закона распределения, соответствующего реальной статистике наработок между отказами объекта или его элементов.

Анализ типовых - характеристик показывает, что каждый из этапов эксплуатации (рис.1) может быть поставлен в соответствие определѐнному теоретическому закону распределения. Обычно период нормальной эксплуатации соответствует экспоненциальному закону, а этап старения и повышенного износа может быть описан нормальным законом, законом Вейбулла или некоторыми другими. Приработка элементов может быть представлена многими теоретическими законами распределения и в том числе законом Вейбулла. Таким образом, изучение особенностей основных теоретических законов распределения может способствовать выявлению закономерностей изменения показателей надѐжности реальных объектов.

Экспоненциальное распределение

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

t

d

P t e 0

при = const превращается в зависимость, соответствующую экспоненциальному закону распределения

P t

e

t

 

,

(1)

 

 

плотность которого равна

 

 

 

f t

e

t .

(2)

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

Постоянство интенсивности отказа позволяет получить выражение

4

dP t dt

const

 

P t

 

.

(3)

 

Отрицательное значение производной в (3) свидетельствует о том, что с увеличением наработки происходит снижение темпа уменьшения вероятности

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

 

 

 

 

Вероятность отказа при

экспоненциальном

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

Q t

1

e

t

(4)

 

Графики основных зависимостей экспоненциального распределения показаны на рис.2.

Рис.2 Графики экспоненциального распределения

Математическое ожидание, равное наработке на отказ, определится для экспоненциального распределения из выражения

 

T

e

t dt 1/

 

 

 

 

0

.

 

(5)

Как следует из (5), средняя наработка на отказ экспоненциального рас-

пределения обратно пропорциональна

интенсивности отказов. Это позволяет

записать,

 

 

 

 

что

P t

e t T

и

Q t 1 e

t T

(6)

 

 

 

.

Из (6) следует, что при экспоненциальном законе надѐжности достаточно знать среднюю наработку на отказ, чтобы определить вероятность безотказной работы объекта в любой момент времени.

Важной особенностью экспоненциального закона распределения является то обстоятельство, что вероятность безотказной работы в каком-то интервале времени не зависит от предшествующей наработки объекта, а зависит только от величины самого интервала. Допустим, что объект благополучно проработал время t и нас интересует вероятность безотказной работы на интервале t . Эту условную вероятность обозначим P( t/t). В то же время вероятность работоспособного состояния объекта к моменту (t+ t) можно рассматривать как вероятность сложного события, заключающегося в том, что объект безотказно отработал время t (одно событие) и затем также безотказно время t (второе собы-

5

тие). Тогда по правилу умножения вероятностей P(t+

t) = P(t) P( t/t), откуда

P( t/t)= P(t+ t)/ P(t).

 

 

 

 

 

 

 

 

 

Для экспоненциального закона получаем

 

 

 

t

 

e

t

 

t

t

 

P

 

 

 

 

e

 

t

 

 

 

 

 

 

 

e

t

 

.

(7)

 

 

 

 

 

 

 

 

 

 

 

 

Выражение (7) даѐт основание считать, что независимость вероятности безотказной работы объекта от предыстории нагружений может быть объяснена возникновением отказов только определѐнного класса, а именно внезапных отказов.

Распределение Вейбулла

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

 

 

m

m 1

 

 

t m t 0

 

f t

 

 

t

 

 

 

e

 

 

 

 

 

t 0

 

 

 

 

 

.

(8)

При 1/t0 = и m = 1 уравнение (8)

 

превращается в плотность экспонен-

циального распределения. Величина 1/t0

 

определяет масштаб, а m – асиммет-

рию (форму) распределения.

 

 

 

 

 

 

 

 

 

 

 

После интегрирования (8) от 0 до t получаем

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

F(t), равную Q(t) :

 

 

 

 

 

 

 

 

 

tm

 

 

 

 

 

 

 

 

 

 

 

 

 

Q t

 

1 e

 

t0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

(9)

Следовательно,

 

 

 

 

tm

 

 

 

 

 

 

 

 

 

t0

 

 

 

P t

e

 

.

(10)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отношение плотности

(8) и вероятности (10) даѐт интенсивность отка-

зов

 

 

 

 

 

 

 

 

 

 

 

 

t

 

m

t

m 1

 

 

 

 

 

 

 

 

 

 

 

t 0

 

 

.

(11)

 

 

 

 

 

 

 

 

Основные графики распределения Вейбулла показаны на рис.3.

а)

 

 

 

 

 

 

 

 

 

 

б)

Рис.3 Графики распределения Вейбулла

6

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

Средняя наработка на отказ при распределении Вейбулла определяется из условия

t m

T e t 0 dt 0

и равна

где

 

T t1 m

1

1

 

 

 

 

 

 

 

 

 

 

0

m

 

 

 

 

 

 

,

(12)

 

 

 

 

 

1

1

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

- гамма – функция;

 

 

 

 

1

1

 

u1 m e u du

 

 

 

 

 

 

 

 

m

0

.

 

 

 

 

 

 

 

 

 

 

Нормальное распределение

Двухпараметрическое нормальное (гауссово) распределение исключительно широко применяется в практических задачах теории надѐжности. Параметрами этого распределения является – математическое ожидание случайной величины и - среднеквадратическое отклонение. Плотность нормального распределения определяется зависимостью

 

1

 

 

x

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f (x)

 

e 2 2

 

 

 

 

 

 

 

 

 

 

2

.

(13)

 

 

 

 

Функция распределения F(x) (рис.4) при нормальном законе определяет-

ся интегралом от плотности f(x)

с пределами интегрирования от -

до + .

Случайная величина t как и во всех задачах надѐжности имеет смысл наработки объекта и поэтому определена на положительной полуоси чисел, а нормальный закон, как уже отмечалось, определѐн на всей числовой оси от - до + . В связи с этим в теории надѐжности рассматривают усечѐнный нор-

мальный закон, плотность которого f (t) определяется путѐм умножения (13) на постоянный множитель

f (t) c f (t) ,

7

 

c

1

 

 

 

 

 

 

где

F(b)

F(a)

a,b – левая и правая границы усечѐнного рас-

 

 

 

 

пределения. F(a),F(b) –

значения функций распределения нормального закона

на левой и правой границах усечения.

Смысл постоянного множителя с становится ясным при рассмотрении графика плотности нормального распределения, представленного на рис.5.

Рис.4 Плотность и функция распределения нормального закона при математическом ожидании =0 и среднеквадратическом отклонении =1.

x x /

Рис.5 Плотности нормального и усечѐнно – нормального распределений при = 2, = 3. Параметры усечения a=0, b=10.

Известно, что площадь под кривой плотности распределения всегда

 

 

 

 

 

 

f t dt 1

 

должна быть равна единице, то есть в данном случае 0

 

 

. Как показа-

но на рис.5 для обеспечения этого условия кривую плотности f

t усечѐнно-

го нормального закона приходится сдвигать вверх и вправо путѐм умножения исходной плотности нормального закона на постоянный множитель. Соответ-

8

ственно будут меняться основные параметры: математическое ожидание и среднеквадратическое отклонение. Расчѐты показывают, что при отношении /< 0.5 (коэффициент вариации) постоянный множитель c для усечѐннонормального закона близок к единице. Поэтому во многих практических задачах теории надѐжности пользуются параметрами нормального закона распределения случайной наработки объекта до отказа. При этом математическое ожи-

дание отождествляют со средней наработкой до отказа Т0.

Вероятность безотказной работы при нормальном распределении равна

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

2

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

T0

 

 

 

 

 

 

 

P t

 

 

 

 

 

 

 

 

e

2 2

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

.

 

 

(14)

Вероятность отказа рассчитывается по формуле (при с

1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

t

 

 

2

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

T0

 

 

 

 

Q t

 

 

 

 

 

 

 

 

 

 

 

 

e

2

2

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

(15)

Интенсивность отказов определяется отношением плотности к вероятно-

сти безотказной работы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T0

 

 

 

 

 

 

 

 

 

 

 

t

 

e

2

2

 

 

 

 

 

 

t

T0

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

 

 

2 2

 

dt

(16)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Интегралы

в выражениях (14) (16) не выражаются через элементарные

функции. Обычно

они представляются через интеграл вероятности от парамет-

ра

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z = (t-T0)/

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

e z

 

 

2dz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

(17)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

для которого составлены таблицы.

С учѐтом (17) вероятность безотказной работы при нормальном законе

определяется по формуле

 

 

 

 

P t

1

t T0 /

.

(18)

Исходная информация: - выборка чисел, представляющих собой наработки объектов до отказов;

Алгоритм выполнения работы:

1. В основе работы лежит статистическая обработка заданного вариационного ряда наработок до отказа, которая включает:

вычисление статистических оценок параметров распределения

T

1 N

ti

 

 

 

 

N i 1

 

 

,

(19)

9

1

N

 

 

 

 

ti

T .

 

N 1 i 1

(20)

построение гистограммы распределения числа отказов (рис.6)

Рис.6 Гистограмма частоты отказов

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

k = 1 + 3,3 lg N ,

где N - объѐм выборки.

вычисляют критерий χ2 (критерий «Хи-квадрат» Пирсона)

 

k

l

 

M

2

 

 

2

i

i

 

 

 

 

 

 

 

i 1

 

 

M i ,

(21)

где li - количество отказов в i-й группе наблюдений; Mi – математическое ожидание числа отказов в i-й группе наблюдений при принятой гипотезе о виде теоретического закона распределения; k - количество групп наблюдения.

Процедура использования критерия согласия χ 2 обычно такова. По формуле (21) оценивается мера расхождения между статистическими данными и оценкой их теоретического значения. Затем определяется число степеней свободы r и с помощью таблиц квантилей распределения χ 2 находится вероятность того, что величина, имеющая χ 2 распределение со степенями свободы r, превзойдѐт подсчитанное по (21) значение. Когда эта вероятность χ мала, проверяемая гипотеза отбрасывается как неправдоподобная.

Если обозначить через

2 ,r табличное значение квантили распределе-

ния χ 2 со степенями свободы r, то в случае

2

2,r можно с вероятностью

не менее χ утверждать, что статистическая гипотеза может быть принята.

10