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

книги из ГПНТБ / Батяев, Б. Г. Исследование надежности и точности линейных систем автоматического управления

.pdf
Скачиваний:
14
Добавлен:
19.10.2023
Размер:
5.67 Mб
Скачать

20

формулой (1.54)

©о ^ Jti й

£>(Т)=2 Jtp(Ofl/f-(f) =J2Jt е" dt -(-£) =(f)2. (1.62)

 

о

 

 

 

с

 

 

|

Принимая во внимание формулу

( I . 6 I ) ,

получим

 

 

 

 

 

 

 

 

 

(1.63)

Для экспоненциального закона надежности справедливо следующее

характеристическое

 

свойство

 

 

 

 

 

 

 

 

p ( t + 6 )

е

- x ( i + 6)

(1.64)

 

 

 

 

 

 

 

 

Приведем еще

оцно

замечательное

свойство

этого закона. Из выражения

( I . 6 I ) следует", что для любого t

 

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

—L____ (_ - —— ----i—-

- .—t----------- (1.65)

ait)

fi(t)

ait,)

pit,)

''‘

a(tn)

pita)

Покажем, что

если разность

t_______.— -

с > o

есть постоянная

 

 

 

 

act)

fid)

 

 

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

Действительно, пусть разность

act) pet)

есть-некоторая постоянная величина. Тогда

a i t ) = с a i t ) p e t ) .

21

Учитывая ( 1 . 3 5 ) , будем иметь

a(t)

(, ) _

с а а(*)

 

иди

рСО = Сo.(t) .

 

 

q,(t)

 

 

 

 

 

 

 

Откуда, используя

( I . I 6 )

и равенство

i- =

С

, подучим

 

d i

 

~

Л ’

 

(1 . 66)

 

 

 

 

Интегрируя полученное уравнение

от С

до

t ,

найдем

 

 

t

 

 

 

 

in p ( t ) ^ - \ x d t * - X t .

 

 

 

с

 

 

 

 

Таким образом,

 

 

 

 

 

 

 

p i t ) =

- A t

 

 

 

 

£

 

 

 

Закон надежности на основе гамма-распределеная Найдем характеристики надежности для случая, когда время воз­

никновения отказа имеет гамма-распределение. Для этого закона на­

дежности частота отказов выражается'формулой [2 8]

 

a W

= ~

r

W

“ ~

e ~ * ‘ >

 

(1 -67)

где cL

и А - параметры гамма-распределения,

принимающие положи­

тельные

значения.

элемента при целом числе d

 

 

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

будет равна

 

 

 

 

 

л tI

о£ •

- A t

d i .

 

 

$ ( * ) - U

m

t

=J*

Г ( х )

е

 

После замены переменной

At - %

ж вычисления последнего

интеграла

получим

 

,

 

 

 

 

 

 

 

 

( t ) - T7T\ J 2

e

d z = i - e

Ц

l !

(1 . 68)

 

ЛЧ)

 

 

 

 

 

i ~o

 

 

функцию надежности найдем по формуле ( 1 .4 ) , подставляя вместо q , ( i ) ее значение из (1 .6 8 )

/w=/-?«) =<rAtz;

1 А 0 ‘

 

 

(1 .69)

 

 

'

 

i =0

" а

 

 

 

 

На основании формул (1 .6 7 ) и (1 .6 9 ) определим интенсивность

отказов

 

 

 

 

 

 

 

 

 

л ( 0 =

a ( t )

x * i

(Ло 4

 

(1 .7 0 )

p ( t )

гп

aL~l^

 

 

 

 

 

ц

1*

/

 

 

 

 

 

 

 

г - о

 

 

 

 

Подставляя в выражение

(1 .3 5 )

значения а ( 7 )

и

<^{t)

,

получим

" a ( t )

 

 

 

 

 

 

 

м*\=

$ ( * )

 

 

 

 

 

 

( I .7 I )

 

 

 

 

 

 

 

 

 

Применяя формулу

( 1 .4 7 ),

вычислим математическое

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

мени безотказной

работы

 

 

 

 

 

 

СМ3

 

 

ос

 

 

 

Оо

 

 

f -

 

S j j g J A - t

- j f e )

 

-

0

 

 

о

 

 

 

 

 

 

откуда

 

 

 

ги+о

 

 

 

 

 

 

 

 

т=

cL

 

 

(1 .72)

 

 

 

Л Г(<*)

 

 

 

Для вычисления дисперсии воспользуемся формулой (1 .5 4 ) .

Получим

оо

j

а/ —4

 

 

 

 

 

 

Г(е<-+2) ( r ) Z -

d Z _ oi.

Л '

Л * Г ( с с )

'

~

A *

A* " a 1

23

\ИЛИ

В зависимости от изменения значений параметров d и Л мозшо

установить относительное влияние мгновенных и постепенных отказов

на продолжительность работы элемента.

 

 

 

 

 

При

d = /

из формул ( 1 . 6 7 ) ,

( 1 , 1 7 )

и ( 1 . 7 3 )

следует, что

а ( { ) - Л £

,

Т -

 

и £ )( Т )

= ( Т ) , т . е . гаммараспределение

сводится

к экспоненццел9ному закону.

 

 

 

 

 

 

Интенсивность отказов при d < I

уменьшается,

при d - /

постоян­

на,

при

d

I

увеличивается.

 

 

 

 

 

 

 

 

 

 

 

 

Закон надежности Вейбулла

 

 

 

 

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

формулой

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

—d iI

 

^

^

 

 

 

 

 

 

 

 

а (0 = оСЛв

t

 

7

 

 

 

( 1 . 7 4 )

где

d

и

Л- _

параметры распределения Вейбулла [2 8 ].

 

Подставляя

значение

« ( О

из выражения ( 1 . 7 4 ) в формулу ( I . I 7 ) ,

найдем вероятность отказа

 

 

 

 

 

 

 

 

 

 

 

 

 

г

 

 

 

!<*■

 

 

 

, <*-. t

-с-е

 

 

 

. v

(

. d -1 - Л * . ,

= ~ е

~ d t

 

 

^(0 = J^ca/

е

d t

 

( 1 . 7 5 )

 

функцию надежности выразим через вероятность

отк аза, а

затем

ш есто ^

( t)

подставим

ее

 

значение из ( 1 . 7 5 ) . Тогда

 

рсоч-уа) =е- x t .* (i.76)

Интенсивность отказов найдем на основании формул ( 1 . 7 4 ) и ( 1 . 7 6 )

M i )

a ( t )

,^ 4 - i- i

p t t )

( 1 . 7 7 )

 

 

 

24

Подставляя в выражение (1 .3 5 )

значения a ( t ) и <^(i) . получим

функцию ресурса надежности

 

а (0

еСЛ t dL-i

 

( 1 . 7 8 )

Применяя формулу ( 1 . 4 9 ) , вычислим математическое ожидание времени безотказной работы

оо

ос

л 4 d'

Т = \P(t)oLi

= \ е~

d t .

ОО

Произведя замену переменной = 2 , и вычисляя интеграл, получим

*

:

Г ( с ± )

• (1 .7 9 )

 

 

 

Дисперсию времени безотказной работы вычислим по формуле

( 1 . 5 4 )

ОО

оо

/

 

£ > ( T ) = z j { / > « ) < l t - ( f ) ‘ ~ z j t e

d t - ( f )

»

•(1.80)

oLA

Из последних выражений следует, что при оС - / распределение Вейбулла превращается в экспоненциальное, так как

л(0 = х ?

, й0 (т ) = - L ^ ( f ) 2-.

 

I

 

Закон надежности Рэлея

Рассмотрим основные характеристики надежности при распределе­

нии времени отказа по

закону Рэлея, В этом случае функция надежно ста

 

 

 

 

 

25

 

 

выражается

формулой

 

t z

 

 

 

 

 

е

 

 

 

 

P(t)

Z в *

 

( I . 81)

 

 

 

>

 

 

 

 

 

 

где <э -

параметр распределения

Рэлея.

Тогда вероятность

отказа

элемента равна

 

 

 

 

 

 

 

 

 

 

 

t*

 

 

 

 

$(t) = i - e

 

г в \

 

^-82)

Дифференцируя выражение

( 1 . 8 2 ) по t ,

получим частоту отказов

 

 

 

 

 

t *

 

 

a(t) = p'(t) = - ~ г

&

а.вз)

Применяя формулу

( 1 . 2 5 ) ,

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

 

 

 

л(0= а(0

t

 

( 1 . 8 4 )

 

 

 

p i t )

 

 

 

Принимая во

внимание формулы ( 1 . 8 2 )

и ( 1 . 8 3 ) , найдем функцию

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

 

 

 

 

 

 

 

\ a ( t )

 

 

t

 

 

 

J b it .

 

 

 

w

( 1 . 8 5 )

 

 

 

 

 

 

Математическое ожидание времени безотказной работы вычислим по формуле

оо

Y - ~ т г

7 ' ~ г ~

 

T = j p ( O d i = j e Z6dt=c>Je

(1.86)

На основании форцулы ( 1 . 5 4 )

вычислим дисперсию времени безот­

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

 

 

 

оо

 

оо

 

<£>(т) =2 jtpCOdt -(т)2~£J е

dt -(т)я=

(f)S

или

£ > (т ) = J ( 4 - ± ) = ( т / ( ^ р " / ) .

(1 .8 7 )

Нормальный закон надежности

 

 

Основные предпосылки применения нормального закона выражены в

центральной предельной теореме Л.М.Ляпупова [30]

. Нормальный закон

надежности характерен для постепенных отказов,

возникающих в

ре­

зультате износа элементов. Учитывая, что случайная величина Т при­ нимает неотрицательные значения, введем коэффициент пропорциональ­ ности в формулу частоты отказов.

 

Тогда

получим

 

 

 

 

( ± - т У

 

 

 

 

(1 .8 8 )

где

Те

и

б с -

параметры нормального закона распределения, а

Д

- коэффициент пропорциональности. Определим , используя

с в о й с т е о

(I .14)

частоты отказов. Для этого проинтегрируем (1 .3 8 ) в

пределах

от

0

до

Выполняя в интеграле справа замену переменной

- i , будем

иметь

 

бо

 

 

О-

в - Т с ) 2'

 

То

бо

где

 

 

 

к =

i

 

(1 .8 9 )

 

 

 

 

 

 

 

 

0,9

+ ф Г Х » 4,

 

 

 

 

 

?

I 0 \6 0 /

 

 

Пероятн.сгь отказа элемента 'определим по формуле ( I .IV ) ,

 

подставляя вместо ct(t)

ее

значение.

Получим

,

 

^ 0

i z l t i o jo

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

\9.(г£) +п {£

(1 .9 0 )

 

?<0=

[о7»"+ч(й)]

 

Подставляя в выражение

(1 .4 ) вместо q ( i )

ее значение,

получим

;футстд а

надалности

 

 

 

 

 

 

2( ' # ) + ^ ( й ]

ч * - « ( ч г )

( I . 91)

 

р и ) = 1

 

 

* >*<%($:)

 

 

 

 

 

pt+WJ.',]

 

данное

распределение удовлетворяет начальным условиям р ( о ) - !

и

q (о) - О , что подтверждает правильность выбора коэффициента про­

порциональности.

 

 

 

 

 

По формуле ( I .J 3 )

надаем интенсивность отказов

 

 

 

 

 

(_ ^ -т 0) й

 

 

 

0_

 

б 2 б /

 

(1 .9 2 )

 

M i ) =

 

 

 

 

 

/> ( 0

На основании формулы (1 .3 5 ) определим функцию ресурса надежности

2 Ш

^ ___________________ ^ ________________________________

( Т 93)

wт&Ым+ЧЦ*#)]

§2 . Надежность элемента многократного действия

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

Эти элементы предназначены для длительной работы, поэтому их назы­ вают элементами многократного действия [49] .

На обнаружение и восстановление неисправного элемента затрачи­ вается определенное время, которое называется временем восстановле­ ния. Величина времени восстановления Т зависит от большого числа факторов, которые не могут быть учтены полностью. Поэтому время вос­ становления элемента многократного действия является случайной ве­ личиной. Будем считать, что функция распределения времени восстанов­ ления V ( ' t ) задана и является непрерывной функцией времени.

Определение, 'функцияУ(т )~Р(Т<?), выражающая вероятность того, что элемент будет восстановлен до момента Т , называется вероят­ ностью восстановления.

Вероятность восстановления V ( 'r ) является неубывающей функ­ цией времени и заключена в пределах 0 & Y ( l ) & l . Функция V ( t ) яв­ ляется количественной характеристикой восстанавливаемости элемента. В связи с этим, она используется в качестве критерия надежности [48].

Наряду с вероятностью восстановления Y ( t ) часто употребляется и другая функция

(Нт) =j-V(r)., &-U

т .е . вероятность невосстановления элемента за время Т . Предпо­ ложим, что функция V ("£ ) дифференцируема и имеет непрерывную произ­ водную V ( Т ) , которая называется плотностью вероятности восстано­ вления

(2 .2 )

29

Плотность вероятности восстановления характеризует скорость' измене­

ния вероятности восстановления

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

.

По величине V ( r ) можно судить

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

которые могут быть

восстановлены за определенный промежуток времени. Очевидно, что ве­

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

% может быть определе­

на интегрированием функции 1>(г)

в пределах

от О

до

Z

 

Т

\

 

 

 

] [ { г ) = J v ( v ) d v .

,

 

( 2 . 3 )

 

о

 

 

 

 

Подставляя в выражение ( 2 . 1 ) вместо Mir) ее

значение

из (2-.3),

получим

 

 

 

 

 

в ( ъ ) = l - j v ( T ) c h : .

 

( 2 . 4 )

 

 

с

 

 

 

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

элемента G i ^ / r )

в промежутке времени ( Z , Z,)

при условии,

что до момента Z эле­

мент не был восстановлен. Для

этого необходимо воспользоваться фор­

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

именно

 

 

 

 

 

 

 

б- { % / z ) ~

 

 

(2*5)

Тогда вероятность

восстановления

элемента на интервале времени

( X , Zf ) при указанном условии будет равна

 

 

 

У ( ъ М = 1 - 6 - ( Ъ / г )

 

(2 .6 )

 

2 . Интенсивность

восстановления.

 

 

 

Определение, функция, выражающая плотность условной вероятно­

сти

восстановления в момент

Г

при условии, что до

этого момента

элемент не был восстановлен,

называется

интенсивностью восстановле­

ния.

 

 

 

 

 

Г, -Z+ д Z-

вычислим интенсивность

 

Полагая в формуле ( 2 . 6 )

восстановления

jJi ( г ) ,

как предел]

отношения

 

 

 

 

V ( * + * Щ

при

Д Т

о ,

 

 

 

 

 

Д X

 

 

 

 

 

 

а I X ) - и *

V ( ? + a ? /r )

1

 

бЮ -С С с+ л*)

- - т у ( 2 - 7)

Р ё —

=

 

 

J r ----------- =

J

6.

с

 

 

 

 

 

 

 

 

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