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

книги из ГПНТБ / Дружинин Г.В. Надежность электрических схем авиационных систем

.pdf
Скачиваний:
26
Добавлен:
29.10.2023
Размер:
6.06 Mб
Скачать

Введением начального условия р(0) •

I накладываем ограничение,

чтобы к началу эксплуатации элемент был исправным.

По формуле (1.6) можно, имея

статистическую кривую >.(£),

найти вероятность исправной работы элемента в течение заданного времени I или найти с заданной вероятностью р время / исправной работы элемента.

p(ti) p(hM

В практических задачах иногда бывает нужно знать вероятность исправной работы в течение времени (П, t2), а нс времени (0, t2), т. е. ставится вопрос о переносе начала отсчета времени (рис. 1.2). Для этого случая запишем формулу (1.6) в виде:

 

 

-

J). (t)at -

.1 i (0 at

 

 

 

 

р [U) — t

0

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

и

>(t)at

 

 

 

 

 

 

•- ,f

 

 

 

 

P(h) = P V\ ) z

Л

 

 

(1-7)

Вместе с тем, согласно теореме умножения вероятностей,

 

 

 

0

P{t\)

 

 

 

 

(1 -8)

где p{t2)

— вероятность

исправной

работы

элемента

в

течение

 

времени (0, 12);

 

 

 

 

 

 

p(t\)

— вероятность

исправной

работы

элемента

в

течение

 

времени (0,

Л);

 

 

 

 

 

 

p(t2H\)

— условная вероятность

исправной работы

 

элемента

 

в течение времени (7Ь /2), найденная в предположении,

 

что в момент времени Л элемент был исправен.

Таким образом, вероятность исправной работы в течение вре­ мени (/ь t2) элемента, который был исправен к началу этого периода-, будет

t,

- I ). (<) at

(1.9)

10

Рассмотрим систему, состоящую из т разнородных элементов.

В соответствии с (1.2) и (1.6) функция

надежности системы будет

определяться выражением

т

t

 

 

 

 

 

 

Л

 

(/) dt

P( t ) :

У-10

 

;1 -1 о)

или

 

 

 

 

Р (t) =

-

f Л (/)

( I t

e

0

 

( 1. 11)

где обозначенная через Л (t) интенсивность выхода из строя системы равна сумме интенсивностей выхода из строя элементов:

т

 

•МО =

2

W

(Ы 2)

Для

системы, состоящей из

j= 1

 

 

одинаковых элементов,

 

Л (0

 

 

 

 

и соответственно

 

t

 

 

 

 

 

 

 

 

P{t) = e

- - m

i / .

(/) dt

(1.13)

 

 

ь

.

Для

периода времени (/ь 0)

формула (1.11)

будет иметь вид:

 

 

 

и

 

 

 

 

-

I

Л (t) dt

(1.14)

 

P[t'>it\)~ е

*'

 

t3

Вычисление интегралов j \ { t ) d t или | X(Х)йД можно произво-

 

 

х

л

под

дить численно или графически, планиметрированием площади

кривыми A (t)

или ХДО в пределах от 0 до 0-

кри­

Формулы

(1.6) —(1.14)

позволяют, имея статистические

вые Х;- (t) для

элементов,

найти

вероятность исправной работы

системы в течение заданного времени t или найти с заданной веро­

ятностью Р время

работы

системы

 

без повреждений.

 

 

 

Например, для вычисления веро­

 

ятности исправной

работы

системы

 

в течение времени (Д, t2) нужно

 

измерить площадь 5 (рис.

1.3) и

 

найти по таблице

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

 

функции:

 

Р и с .

1.3.

Р ( Ш = е - * -

При заданной вероятности Р сначала находится значение инте­

грала

 

 

 

 

|\\. { t ) d t = - In Р.

 

 

X

 

 

Затем, имея t\, по

графику

Л(П находим (Д — В),

при котором

площадь под соответствующим участком кривой равна —In Р.

И

Таким образом, вычисление функции надежности системы сво­ дится к суммированию площадей под графиками интенсивности выхода из строя элементов и нахождению по таблице экспоненци­ альной функции значений P(t).

Процесс вычисления P(t) значительно усложнится, если поль­ зоваться плотностью вероятности времени исправной работы эле­

ментов f(t) = —— — . функция надежности P(t) в этом случае, dt

в соответствии с формулами (1.1) — (1.2), будет определяться выра­

жением

 

т

(

 

тп

 

 

P(t) = U

[ l - q

[ t ) ] =U

[1 - j

( 1Л5>

j-i

(1.15)

J-i

о

вычисление

Уже по виду формулы

можно

заключить, что

функции надежности системы по плотностям вероятности времени исправной работы элементов в общем случае приводит к очень громоздким вычислениям. В общем случае пользоваться графи­ ками fj (0 при вычислении функции надежности системы P(t) практически возможно лишь при замене формулы (1.15) каким-либо приближенным выражением.

Вычисление функции надежности системы по плотностям веро­ ятности f j (t) элементов несколько упрощается лишь в случае, когда графики fj(t) достаточно точно совпадают с теоретическими зако­ нами распределения: нормальным законом или законом равной веро­ ятности и при этом систему можно разбить на небольшое число групп одинаковых элементов. Кроме того, при малых периодах работы

(при j A (t) dt 1), можно правую часть формулы (1.11) разло­

жить в степенной ряд и, пренебрегая членами в степени выше пер­

вой,

написать:

t

 

т

t

 

 

 

 

 

 

j \j(t)dt .

(U 6)

 

 

 

0

d t = i -

2

 

 

 

 

y-10

 

Когда, время

работы системы

мало

и

выполняется

условие

т

t

 

 

 

 

 

 

 

2

j

/,• (t) dt С

1, формулу

(1.15)

также можно переписать в виде

У-1 О

 

 

т

i

 

 

 

 

 

 

 

 

 

 

 

 

 

P(t) -

1 - 2

]'/f(t)dt.

(1.17)

J=io

Таким образом, при малых t можно не делать различия между

(t)

и

В этом

можно также

убедиться,

переписав формулу

(1.5)

в виде:

 

f ( t ) = \ ( t ) p ( t ) .

 

(1.18)

 

 

 

 

При малых / значения p(t) ^

1 и f(th ^

I (Л. При I =

0 p(t)

= 1

и f ( t ) = l ( t ) .

 

 

 

 

 

 

При

Согласно (1.18),

всегда справедливо

неравенство

1-( f ) > /(0 -

Х= const

f(t)

убывает с течением времени работы по экспо­

ненте.

 

 

 

 

 

 

12

Вычисление и использование интенсивности выхода из строя приводит к достаточно точному и простому решению задачи. Воз­ можность применения формул (1.6) —(1.14) ограничивается лишь хорошо согласующимся с практикой требованием практической невозможности совмещения двух и более первичных неисправностей в один и тот же момент времени.

Для применения изложенного метода оценки надежности эле­

мента необходимо

знать зависимость от времени интенсивности

выхода из строя

Эта характеристика должна определяться на

основании экспериментальных исследований. Графики должны ото­ бражать действительность, а не «подгоняться» под теоретические законы распределения, например, закон X= const. Только в этом случае можно использовать метод расчета надежности в качестве инструмента для повышения надежности элементов и систем. Неко­ торые примеры использования графиков Х(t) для решения ряда инженерных задач будут рассмотрены в § 1.4.

Вдальнейшем изложении графики>X(t) и АД) будут иногда называться «ламбда-характеристиками»:

X(t) — ламбда-характеристика элемента;

Л(t ) — ламбда-характеристика системы.

Внекоторых статьях надежность элементов и систем характери­ зуется средним сроком службы элемента или системы, под которым следует понимать математическое ожидание срока службы элемен­

та. Согласно определению математического ожидания,

ооос

^0= (

tdq(t) — — f

tdp (t).

(Ы 9)

о

о

 

 

Интегрируя по частям, имеем

ОО

оо

 

t0= - t p { t ) \

+ J p ( t ) d t .

(1.20)

оо

При конечных X функция надежности быстро уменьшается с тече­ нием времени. Поэтому существует предельное соотношение

 

lim t-p(t) = 0.

(1.21)

 

t ■' эо

 

Справедливость формулы (1.21) следует из того, что при t —то

и конечных X экспоненциальная функция p(t)

0 быстрее,

чем t -

то .

 

Из

(1.20) и (1.21) имеем при X —-- const:

 

 

K — l е htdt —

( 1. 22)

 

 

Поэтому при X -- const функция надежности может быть записана в виде

о (^)= е~и = е t„

(1.23)

13

Иногда средний срок Службы неудачно называют средним вре­ менем между отказами илисредней наработкой на один отказ. Для элемента неточность этих выражений состоит в том, что он может отказать только один раз.

Формула (1.23) часто называется «экспоненциальный закон надежности».

Учитывая (1.12), можно написать формулу для вычисления среднего срока службы системы ^с„ст по средним срокам службы элементов tj 0:

т

1

(1.24)

^СИСТ

Следует отметить, что средний срок службы элемента является понятием условным, так как обычно ни один элемент не работает столь продолжительное время. Элементы схем автоматики и радио электроники устаревают гораздо раньше, чем успеют проработать время, соизмеримое со своим средним сроком службы. Более того, значения tj0 обычно вычисляются по экспериментальным данным об отказах элементов в начальный период их эксплуатации. Поэтому tj о можно понимать как условный средний срок службы элемента, который имел бы место в действительности, если бы эле­ мент сохранял в течение всего срока службы ту интенсивность выхода из строя, которой он обладал в начальный период эксплуа­ тации.

Согласно формуле (1.23), средний срок службы системы (также

и элемента) можно понимать как время

работы t =

/0, в течение

которого

система

(элемент) остается исправной

с

вероятностью

P(ta) =

е"1 ^ 0,37.

среднего

срока службы

В дальнейшем

изложении понятиями

элемента и системы пользоваться не будем. Это объясняется как несколько путанным физическим смыслом термина «средний срок службы», так и практическим неудобством вычисления среднего срока службы системы по формуле (1.24). Кроме того, при опериро­ вании средними сроками службы элементов и систем эксперимен­ тальные данные подгоняются под теоретический закон распределе­ ния const, что также часто приносит вред (§ 1.4).

§ 1.3. ОБРАБОТКА СТАТИСТИЧЕСКИХ ДАННЫХ ОБ ОТКАЗАХ

Графики интенсивности выхода из строя A.(t) или плотности вероятности f(t) строятся по статистическим данным об отказах эле­ ментов (или систем). Статистические данные могут быть получены как в результате лабораторных испытаний, так и путем наблюдения за находящимися в реальных условиях эксплуатации элементами и системами. Рассмотрим способы получения и обработки эксперимен­ тальных данных в этих двух случаях.

При лабораторных испытаниях достаточно большая группа эле­ ментов работает в условиях, по возможности приближающихся

14

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

I

с п о с о б .

Разделим весь диапазон времени работы элеме

тов на интервалы

M i ~ t i t-, _ 1

и подсчитаем

количество

отка­

зов Д«г- элементов, приходящееся

на каждый г-тый интервал. Раз­

делим это число

на

число элементов,

оставшихся

исправными

к началу рассматриваемого интервала времени

Att,

и на

длину

этого

интервала.

В результате получим значение статистической

интенсивности выхода из строя, соответствующее

данному

интер­

валу

А/.

 

 

 

 

 

 

 

 

 

 

Л'С

 

 

 

(1.25)

 

 

 

[TV —n,-|

 

 

 

 

 

 

 

 

где

ni — обшее

(накопленное) число отказов элементов за

время

(О,

 

 

 

 

 

 

 

 

Здесь и в дальнейшем статистические величины будут отмечать­ ся звездочкой. Результаты вычислений сводятся в таблицу, в кото­ рой приведены интервалы Atf в соответствии с их порядковым номером и соответствующие им значения к*:

Д*1

at2

. . .

V

) * . . .

>4*

 

2

 

1>

Ч*

Эта таблица оформляется в виде ступенчатого графика рис. 1.4. Для этого по горизонтальной оси последовательно откладываются

интервалы Att и на каждом из них строится прямоугольник с высо­

той a.t*.

Соединив

полученные для каждого интервала времени

Att

значения

Хг* плавной кривой и переходя от статистических понятий

к вероятностным,

получим

график интенсивности

выхода

из

строя X(t).

 

 

 

 

II

с п о с о б . Ход вычислений совпадает с описанным выше, н

число отказавших

элементов

Ani , приходящееся на

интервал

вре­

15

мени А^., делится на число элементов, поставленных на испытание:

А«;

(1.26)

N- Д^-

Соединив полученные для каждого интервала времени At t значе­ ния / * плавной кривой и переходя от статистических понятий к веро­ ятностным, получим плотность вероятности времени исправной

работы f(t)

или,

другими словами,

плотность вероятности момента

повреждения.

и / / ,:

существует та

же

зависимость, чго и между

Между

ХД

вероятностными характеристиками

 

XД) и f(t). Действительно, раз­

делив в формуле (1.25)

числитель и знаменатель на N, получим

 

 

.

A^t

 

 

 

^_ ft*

 

 

 

N-Ltt

 

f*

(1.27)

 

 

 

i - Л .

i

-

ч*

р ?

 

 

 

 

 

 

 

N

 

 

 

 

 

и соответственно, перейдя от статистических характеристик к непре­ рывным вероятностным, получим формулу (1.18).

Так как для практических расчетов надежности обычно исполь­ зуется интенсивность выхода из строя, то при обработке результатов

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

форму­

лой (1.25).

 

У лиц, впервые приступающих к обработке данных об отказах,

обычно

возникают затруднения в определении числа

интерва­

лов А/,-,

на которые следует разбить период наблюдения

t. Число

интервалов не должно быть слишком большим, ибо при этом могут

обнаружиться случайные колебания значений

ХД

Вместе с тем оно

не должно быть слишком малым, ибо тогда мы

получим слишком

грубую кривую АД).

зависит от объема и однородности стати­

Длина интервала

стического материала.

При увеличении числа

наблюдений можно

выбирать все более короткие интервалы. Каких-либо правил на этот счет не существует, но необходимые навыки очень быстро приобре­ таются в процессе работы. Обычно интервал наблюдения разбивают на 7—20 интервалов.

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

интервал

АД (от t =

0 до t —- Л) па

несколько

подинтервалов

с тем, чтобы полнее выявить поведение кривой

АД)

на участке от

t — 0 до максимума,

который обычно

имеет

место

в начальный

период работы элементов.

 

 

 

При

исследовании

поведения элементов в

условиях реальной

эксплуатации сведения об их отказах обычно получаются в резуль­ тате наблюдения за группой работающих систем. В результате наблюдения необходимо иметь время исправной работы каждого из отказавших элементов. При наличии достаточного количества дан-

16

ных обработка их ведется по формуле (1.25). При этом под ^ п ( понимается число отказавших элементов, время исправной работы которых закончилось на интервале A*)- — ti При этом учи­ тываются отказы лишь установленных в системе перед началом периода наблюдения элементов.

В тех случаях, когда элементы или системы проработали раз­ личное время, приходится делить их па группы по числу наработан­ ных часов. При вычислении значений Х/|: для различных периодов работы элемента (системы) используется различное, все время уменьшающееся, число групп. При этом точность построения к(1) несколько снижается по мере роста /.

Рассмотрим примеры построения графиков X{t) по эксперимен­ тальным данным.

ПРИМЕР 1. Под наблюдением находились две группы рабо­ тающих однотипных устройств. Устройства первой группы нарабо­ тали 1100 часов; устройства второй группы - 500 часов. В устрой­ ствах первой группы имеются 40000 элементов определенного типа; в устройствах второй группы ----- 80000 элементов этого же типа. В процессе эксплуатации зафиксирован ряд отказов устройств из-за

выхода

из строя элементов данного типа. По

времени

исправной

работы

отказавшие элементы группируются

так, как

показано

в табл.

1.1.

 

 

Интервалы времени

[■работы (часы)

Число

элементов

1

группа

время

исправно!

работы которых

устройств

 

 

закончилось на

II

группа

данном

интервале

 

 

 

100-50

 

-200300

-300400

-400500

-500600

-600700

700-80»

Таблица

1.1

05-0

о

800-900

900-1 ООП

t000- 1 001

 

 

О

 

 

 

 

 

 

 

 

 

 

 

О

 

 

 

 

 

 

 

 

 

 

 

С4)

 

 

 

 

 

 

 

 

 

 

 

О

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

(3

21

44

23

14

12

9

11

8

11

9

10

14

43

88

49

27

22

 

 

 

 

_

 

 

 

 

 

 

 

Необходймо построить /.-характеристику, используя при этом

наиболее полно имеющиеся статистические данные.

 

 

При решении данной задачи значения /.**

для периода времени

работы 0500 часов вычисляются по статистическим

данным

для

обеих групп устройств, т. е. для 40000 -j—80000 =

120000 элементов.

Значения

для

периода времени 500—1100 часов работы вычис­

ляются, по статистическим данным устройств

первой

группы.

Ход

вычислений по формуле (1.25) иллюстрирует табл. 1.2.

 

По даннылццаф 1 и V табл. 1.2 строится пПЯ&срамма, д датем

график \(1)

(рчс.

1.5). Р " гос. публичная

*

 

-и: - '

 

 

?

*; •

. 2 НАУЧНО -ТЕХНИЧЕСКАЯ

I

^

 

 

2 Г. В. Дружинин Ц.

;- ‘Д Д. I Б ИБЛИОТЕКА СССР

 

 

4t&s.

17

 

И н т е р в а л ы в р е м е н и

1

 

 

р а б о т ы ( ч а с ы )

11

Ч и с л о о т к а з а в ш и х

 

н а д а н н о м и н т е р в а ­

 

л е э л е м е н т о в

111

О б щ е е ( н а к о п л е н -

 

н о е д о д а н н о г о и н ­

 

т е р в а л а ) ч и с л о о т к а ­

 

з о в

IV

Ч и с л о э л е м е н т о в ,

 

п р о д о л ж а ю щ и х

 

и с п р а в н о р а б о т а т ь

V

З н а ч е н и я с т а т и с т и -

 

ч е с к о й и н т е н с и в н о -

 

с т и в ы х о д а и з с т р о я

 

 

 

 

 

 

 

 

 

t

а б л и u а

1.2

 

 

 

 

 

 

 

 

 

 

 

 

о

 

 

 

 

 

 

 

 

 

 

 

о

о

 

О

О

о

о

о

 

О

О

о

о

о

 

о

о

гг

о

 

о

о

о

о

 

см

со

ю

 

оо

 

оз

 

1C

 

 

I

О

 

 

О

1

 

 

 

о

 

1

о

О

о

 

О

о

о

 

 

о

о

О

о

о

 

о

о

о

о

 

°

 

ю

»—«

см

С с

 

 

ю

ю

г--.

о ,

03

 

20

57

132

72

41

34

 

9

11

8

11

9

10

0

20

87

2И)

291

332

119

128

139

147

158

167

о •

О

СО

 

03

 

 

 

СМ

 

со

см

со

О

ОС

 

00

о

■ s:

 

00

г -

СО

to

со

со

О

оз

 

г—

г -

со

 

00

со

оо

00

00

6

Оз

оз

Оз

03

оз

СЗ

Оз

ОЗ

 

оз

СУЗ

03

см

 

 

 

 

 

 

со

СО

со

го

«О

со

СО

со

о со

 

со

со

 

«О

со

со

со со

 

О о

о

О о

о

 

О О

о

о

«г

о

 

 

 

 

 

*■“'

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

СО

с

О

 

00

 

СМ

г -

о

t—

ю

 

со

Оз

 

«о

со

СМ

 

CN

см

см

см

см"

см

ПРИМЕР И. Имеются приведенные в табл. 1.3 статистические данные об отказах трех групп одинаковых устройств.

Устройства разделены на группы по числу наработанных часов:

I

группа —

550 час.

II

группа —

400 час.

III

группа

—200 час.

В

каждой группепо 100 устройств. Построить Л-характеристику

рассматриваемого устройства.

Вычисления по формуле (1.25) сведены в табл. 1.4.

18

Интервалы времени

работы (часы)

Число

устройств, время

1 группа

устройств

исправной

работы которых

II

группа

закончилось на данном нн-

устройств

гервале

 

 

III

группа

 

 

устройств

 

 

 

 

 

о

О

 

О

-057 5

75—100

ю

 

CN

 

ю

 

ю

 

 

7

CN

1

 

 

1

 

ю

 

 

о

4

 

8

6

3

5

6

 

9

5

4

5

5

1 8

7

5

6

150-200

Т а б л и ц а

1.3

200-25)

 

-250300

300—400

тг

 

 

 

 

 

О

 

 

 

 

 

ю

 

 

 

 

 

ю

 

 

 

 

 

о

 

 

 

 

 

о

4

I

 

9

; 3

5

3

3

 

2

1

 

 

, 4

 

 

 

 

 

1

 

3

 

i

~

I

 

 

 

 

 

Таблица 1.4

 

 

1

 

I группа

 

 

 

 

II

группа

 

1 ш

 

 

 

 

 

 

 

 

 

 

 

 

 

[ группа

Mi

0 -2 5

2 5 -5 0 5 0 - 7 5

7 5 -

 

10 0 -

15 0 -

2 0 0 -

2 5 0 -

3 0 0 -

400

&Hi

15

25

18

100

150

200

250

300

400

550

12

 

16

10

4

 

3

7

5

« - S

м н

0

15

40

58

 

70

86

62

 

66

69

36

N

- n

3 0 0

285

2 6 0

242

 

230

214

138

 

134

131

64

 

 

СО

СО

СО

СО

 

СО

тГ

 

 

■4“

«■Г

к*

|

 

 

 

 

 

1

1

'

1

1

1

О

О

О

О

О

О

О

О

О

о

 

 

 

 

 

 

 

 

 

 

см

lO

о

О

 

1 ■4*

СО

сс

 

ю

со

СМ

 

 

СО

СЧ

СЧ

 

!

СП

to

 

V*

i d

i d

 

 

 

 

 

 

 

1

 

 

 

 

 

 

Построенная по данным табл. 1.4 Л-характеристика рассматри­ ваемого устройства изображена на рис. 1.6.

А , f/vee

f»**-

Оо

so

по

iso

гм iso >оо ш

*ео oso

то

t, *«*

 

 

Р и с .

1.6.

А-характеристика для примера

II.

 

2*

19