
книги из ГПНТБ / Батяев, Б. Г. Исследование надежности и точности линейных систем автоматического управления
.pdf20
формулой (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. |
с |
|
|
|
|
|
|
|
|