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

книги из ГПНТБ / Крауз С.В. Основы технической эксплуатации авиационного оборудования I. Элементы теории надежности оборудования летательных аппаратов. II. Средства и методы контроля и подготовки авиационного оборудования

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

в (1.11) значение f(t), получим соотношение, связывающее в еди­ ную формулу надежность, частоту и интенсивность отказов:

( 1. 12)

Частота отказов f(t) является дифференциальным законом рас­ пределения, т. е.

=

(1.13)

t

Пользуясь формулами (1.12) и (1.13), можно составить выражение, связывающее интенсивность отказов и частоту:

*(*) = _ / ! * ] — .

(1.14)

? / ( 0 dt t

Интегрирование уравнения (1.11) позволяет получить весьма важное соотношение между надежностью и интенсивностью отка­ зов в работе

-

f 40 at

.

(1.15)

p ( t ) = е

«

Подставляя выражение (1.15) в соотношение (1.12), получим формулу для частоты отказов:

t

- J X (г)dt

f[t) = l{t)e °

.

(1.16)

Математическое ожидание времени безотказной работы или, что то же самое, среднее время наработки исследуемых объектов до появления отказа может быть выражено следующим образом:

ОС

Tcp^ $ tf(t)dt .

(1.17)

о

 

Если заменить f(t) согласно (1.16), то среднее время безот­ казной работы выразится следующей формулой:

t

ос

— ( X (рdt

 

 

г ср = j

tl(t) е *

dt.

(1.18)

о

 

 

 

Таким образом, между двумя наиболее распространенными характеристиками надежности ^ (t) и Гср существует однознач­

ная зависимость.

Важной характеристикой надежности также является среднее квадратичное отклонение з величины времени безотказной ра­

10

боты рассматриваемого объекта, показывающее степень его откло­ нения относительно среднего значения времени наработки

a==y r J { t ~ T^ f [ t ) d t

(1Л9)

Если интенсивность отказов X является постоянной величи­ ной, тогда можно считать, что среднее время наработки на один отказ Та может быть определено следующим путем:

Тп _ nicptp_^

(|.20)

п

 

где тср— среднее число систем безотказно работающих в интерва­ ле времени /р, а так как интенсивность отказов в данном случае мо­ жет быть выражена посредством формулы

п

(1.21)

к =

^ср

 

тогда, следовательно,

 

x = _ L .

(1.22)

Т

 

1 п

 

Учитывая, что X =< const, и пользуясь равенсвтом

( 1.15), можно

написать

(1.23)

Р [t) = e - u

или, подставляя (1.22), получить следующую формулу: t

p ( t ) = e 7",

(1.24)

выражающую .'собой экспоненциальную зависимость надежности от времени наработки.

ЗАМЕНА НЕИСПРАВНЫХ ЭЛЕМЕНТОВ ОБОРУДОВАНИЯ ВНОВЬ УСТАНАВЛИВАЕМЫМИ.

СРЕДНЯЯ ЧАСТОТА ОТКАЗОВ С УЧЕТОМ ЗАМЕНЫ ЭЛЕМЕНТОВ И СРЕДНЕЕ ВРЕМЯ НАРАБОТКИ МЕЖДУ ДВУМЯ ОТКАЗАМИ

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

11

тов служит величина математического ожидания времени наработ­ ки рассматриваемого элемента системы или объекта между от­ казами.

Под средней частотой отказов элементов' или объектов обору­ дования /сР (t) понимается количество отказов в единицу времени, отнесенное к первоначально взятому количеству объектов N:

/с р ( 0 =

n(t)

(1.25)

т

 

где n(t) — число отказов за время С

на самолете элементы или

Для случая, когда все отказавшие

объекты испытываемой партии оборудования заменены новыми, n(t) = n0(t) 4- л, (t),

где n0(t) — число отказов первоначально взятых объектов;

ni(t) — число отказов вновь установленных объектов взамен поврежденных,

и частота отказов этого оборудования может быть выражена сле­ дующей формулой:

/ , ( , „

” 4 ' ) + Е Л ).

(1.26)

Р

Л7

Л t

 

Средняя частота отказов первоначально взятых объектов обо­

рудования равна:

 

 

 

JU\/ о

(J0 =

^7 .

(1-27)

Для отказавших n\(t) объектов первоначально взятыми объек­ тами следует считать nQ(t) вновь установленных объектов взамен поврежденных ранее из партии в N объектов.

Следовательно, является справедливым следующее выражение:

/ „ « > = "Д.

(1.28)

n0(t) t

 

Согласно (1.27) и (1.28)

 

(1.29)

nl{ t ) = f 2[t) N f ;

тогда

 

 

<hit)

 

(1.30)

Nt

 

 

Подставляя значение (1.27) и (1.30)

в (1-26),

получим

/с р (^ )= /о (0 +

^/,2(0

(1.31)

или, подставляя в (1.31) выражение (1.16), окончательно получим

следующую зависимость:

t

 

 

t

,dt

 

- j ).(/)<//

+ tt-2(t)e 0

(1.32)

f ef{t) = k(t)e 0

.

12

Среднее время наработки между двумя отказами Т тем боль­ ше, чем меньше частота отказов. Если величина частоты отказов не остается постоянной, а изменяется в зависимости от времени, истекшего с начала эксплуатации самолета или системы, то, оче­ видно, и величина среднего времени наработки между отказами Т также зависит от времени, истекшего с начала эксплуатации.

Строго теоретическое вычисление величины среднего времени наработки между двумя отказами самолетной системы связано с предварительным определением вероятности замены деталей после первого, второго и т. д. отказа. Для этой цели необходимо распо­ лагать данными распределения частоты отказов всех элементов системы и данными условных вероятностей отказов системы. На­ пример, если необходимо найти среднее время наработки системы между первым и вторым отказом, то для вновь замененного эле­ мента, не подвергавшегося ранее рабочему режиму, начало отсче­ та функции f(t) следует брать от нуля, а для всех остальных эле­ ментов начало отсчета следует брать от Т\ — среднего времени на­ работки до первого отказа. По условным значениям 7\ можно най­ ти математическое ожидание времени наработки системы между первым и вторым отказом 7Y Учитывая, что современные системы авиационного оборудования содержат очень большое количество составных элементов и деталей, допустимо без большой погреш­ ности пренебречь влиянием замены одного из элементов системы на общий закон распределения вероятностей времени повреждения системы. Тогда сразу отпадет необходимость производить громозд­ кие вычисления, связанные с условными вероятностями, а точ­ ность вычисления при этом упрощении практически не изменяется.

Математическое ожидание появления второго отказа или, что то же самое, среднее время наработки системы между первым и вто­ рым отказом определяется по формуле.

со

 

г 2= \ { t - r x) f [ t - T x)dt,

(1.33)

г,

 

где f(t Ti) — плотность вероятности времени появления второго отказа после момента появления первого отказа или, что то же самое, частота появления отказа

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

ивторого отказов.

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

время наработки системы между

первым

и вторым •отказами,

находят

математическое ожидание

времени

появления третьего

отказа,

аналогично четвертого и затем k-ro

отказа:

 

 

 

 

 

Тъ= \ Ц - Т 2)/а-Т^<М-,

(1.34)

 

 

г,

 

 

Tk=

]

{ t - T k^ ) f ( t - T k^)dt.

(1.35)

 

1_

,

 

 

13

Производя замену переменной t Tk_1— ~, выражение (1.35) можем представить в следующем виде:

I

сг

(1.36)

 

Tk= ^ / ( - z) dz .

Величина Тк — среднее время наработки аппаратуры сложных систем на один отказ в настоящее время является основным критерием практической оценки эксплуатационной надежности сложных систем.

УРАВНЕНИЯ КЛАССИЧЕСКОЙ КРИВОЙ ЧАСТОТЫ ОТКАЗОВ И НАДЕЖНОСТИ

Классическую кривую плотности распределения отказов (фиг. 1.4) можно рассматривать как совокупность трех кривых рас­ пределения: А — кривой, огибающей точки распределения, соответ­ ствующие целым значениям аргумента, подчиняющимся закону Пуассона:

ГЛ1) = Ц ~ е~Та,

(1-37)

где Та — математическое ожидание для данного .распределения;

Б— прямой параллельной оси абсцисс, характеризующей равномерную плотность распределения

fb[t) = const;

(1.38)

В — кривой, подчиняющейся .нормальному закону распреде­ ления

« - V

 

f A t) = — ^ = e

.

(1.39)

0ВV

 

 

Поскольку площадь, ограниченная

огибающей

кривой f(t)

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

п» .

_ ГЧ.

_

Яв.

•>

(1.40)

"ТТ5

45 —

N

Чв—

N

N

 

 

 

 

+

Яб "f йв = N

 

 

 

(1.41)

14

15

и, следовательно,

 

 

 

 

 

 

 

 

 

 

<7» +

<?б + 9в — 1-

 

(1-42)

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

выражена следующей формулой:

 

 

 

 

 

 

f{t)

= / а (7)

/ б {t)

+ / в №

,

(1.43)

где / а (П; fs(ty,

/в (t)

— частоты

появления

отказов,

соотвег-

ствующие кривым распределения А,

Б и В и равные:

 

 

/

а ( 0

= <?а ур е Га;

 

(1.44)

 

fb{t)

= q 6f b{t)\

 

 

(1.45)

 

 

 

 

 

 

« - v

 

 

/ . « > -

, V

c

 

(1.46)

 

 

 

 

Зв У

2-

 

 

 

 

Подставляя

эти выражения в (1.38), получим уравнение ре-

зультирующей кривой Г:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

« - У

 

f(t) = ^ ~ - e ~ Ta +

?б /б (0

+

2V >

(!-47)

 

г!

 

 

 

 

5в к *,т-

 

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

плотности

вероят­

ности времени появления отказов аппаратуры авиационного обо­

рудования.

>

 

 

 

 

 

У р а в н е н и е н а д е ж н о с т и .

Путем применения интеграль­

ного закона распределения формула (1.47)

может быть преобразо­

вана для получения

зависимости

функции

распределения,

пред­

ставляющей в нашем случае величину ненадежности.

 

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

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

величин. Поэтому:

— для

первого

члена, представляющего

выражение

закона

Пуассона, интегральный закон может быть выражен в виде следую­ щей формулы:

Ч^ Т\

(1.48)

о

— для второго члена, представляющего выражение закона равномерной плотности распределения, если учесть, что для слу­ чая, когда </<С1> ^ то согласно (1.4) и (1.24) можно написать

qu(*) = 1 - е Гб;

(1.49)

16

— для третьего члена, представляющего выражение нормаль­ ного закона распределения:

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.50)

а так как

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f i t )

= ----- —

е

2о *

 

 

 

(1.51)

 

 

 

 

 

 

 

 

 

 

 

 

 

V 2"

 

 

 

 

 

 

 

то,

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

 

 

 

 

( / - V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

<7ш(-) =

-----7 ~

 

 

 

 

 

 

 

(1.52)

 

 

 

 

 

V 2"

 

 

 

 

 

 

 

 

где

т

— время появления отказа;

 

 

 

появления

отказа;

 

Тв — математическое

ожидание времени

 

ав — среднее

квадратичное отклонение срока

службы

объек­

 

 

та от его среднего срока службы.

 

 

 

 

 

Из теории вероятностей известно, что

 

 

 

 

 

 

 

 

 

 

 

«

« -

ГЪУ

 

 

 

 

 

 

 

ч 1 °°)

з 1/2-

j e

2V

d t =

1;

 

 

(1.53)

 

 

 

 

 

 

 

 

 

 

 

 

тогда можно записать

 

 

ОС

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

°°

 

i t - т у

 

 

1

 

t- T „

 

 

1 -

ЧИ! W =

 

 

 

dt-

 

I

*

^ dt

(1.54)

 

 

 

 

 

 

 

V *

 

 

 

 

 

 

 

 

 

 

 

или, что то же,

 

 

 

 

 

i t - ту

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 -

Ч\\1(т) =

 

е

2V

dt .

 

 

(1.55)

 

 

 

 

 

V'2*

 

 

 

 

 

 

 

 

Если приравнять

t -

Т.

х и dt

 

dx и таким

путем пе-

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

7"в ,

то

выражение

(1.53)

примет вид:

00

(/—Т )’

 

 

со

х i

 

 

 

 

 

 

I

 

 

dt ■= —■_ !

 

 

dx = 1,

 

(1.56)

 

 

 

 

 

 

 

 

 

 

V 2'

 

 

 

 

)/'2^

 

 

 

 

 

 

а так как среднее выражение равенства

(1.53)

симметрично отно­

сительно оси ординат,

то

можно написать

 

 

 

 

 

 

 

 

 

 

 

£1

 

1

 

 

 

 

 

 

 

 

 

 

е

2

dx

 

 

 

 

(1.57)

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2. С. В.

Крауз и др.

 

 

 

 

 

/

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ОС. ‘ ПУБЛИ

'- : . т р ■■ ■ .

- . о

 

или, что то же самое,

X

1

(1.58)

\'~1к

X

Первый член правой части выражения (1.58) представляет со­ бой функцию Лапласа:

J*<? 2 dx = Ф ! дс).

(1.59)

Второй член правой части выражения (1.58) является правой частью выражения (1.55) и, следовательно,

~4 гг U

2 d x = — I — - Гв‘

(1-60)

У

J

J

 

 

 

Л'

"

 

 

тогда, согласно

(1.58),

(1.59) и (1.60),

 

 

(1.61)

(1.62)

Умножая полученные выражения ненадежности [см. (1.48), (1.49), (1.50)] на соответствующие коэффициенты и суммируя их, получим интегральный закон распределения отказов в виде следую­ щей формулы:

=

е T^ + q \ \ - e Т\ - Q,

(1.63)

о

иуравнения надежности

р( х )

о

(1.64)

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

службы

в месяцах и годах, время налета или наработки в ча­

сах и т.

д.

18

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

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

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

Если для какого-либо элемента или объекта авиационной тех­ ники g размеры интервалов определяющих параметров равны ве­ личине ад, тогда для этого элемента все основные характеристики

надежности

могут

быть записаны

с соответствующим

индексом:

^ 3' а!г’ Л

ит. д.,

а формула (1.47)

 

для этого случая перепишется

в следующем виде:

 

 

 

 

 

 

 

 

'

п f

'

I

Ч'

(1.65)

 

У

1"

Я6 g16

 

, 2ir

 

 

 

 

 

 

 

 

Если размеры интервалов определяющих параметров кривой распределения системы в-целом равны значению А, тогда коэф­ фициент приведения размеров интервалов распределения отдельных элементов k■размерности интервалов распределения системы мо­ жет быть получен из следующего равенства:

( 1. 66)

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

О*

19

Соседние файлы в папке книги из ГПНТБ