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

книги из ГПНТБ / Волков Е.Б. Основы теории надежности ракетных двигателей

.pdf
Скачиваний:
55
Добавлен:
25.10.2023
Размер:
12.65 Mб
Скачать

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

ний к нахождению выходных характеристик системы

[у{, у2, . . .,

.. ., (/;) в заданных допусках:

 

v= 1, /.

(2.34)

Изложенное позволяет сформулировать рассматриваемую зада­ чу следующим образом. Общее выражение (2. 33) для показате­ ля надежности системы требуется представить в виде некото­

рой композиции вероятностей вида P, = tp(Py), зависящих от

у, каждая из которых относится к /-му элементу (например,

.V

Р = * П Ру, где х < 1 ). Пусть такая композиция найдена. Тогда

; - 1

функция вида

Ру= т(Ру)

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

Принципиальное решение задачи о декомпозиции функций многих переменных дано в работе [42] в виде теоремы: для каж­ дого f e M существует N(2N+ 1 ) таких функций

?,-у(*у), (у=1, 2 , . . . , yV; г = 1, 2........

2N, 2yV + 1),

в которых:

 

 

 

а) все функции ср,-;, непрерывны на интервале [0, 1];

 

б)

для любой функции /(.V1, х2, . . ., .V.V), непрерывной на

0 ^'л'], х2. . ., л'л-<1,

существуют 2/V+1 функций г|-; (/=1,

2/V+1),

каждая из которых непрерывна на /?<‘>, причем f( x ,, х п, ..

., xN)=

2/V+1

/ /V

\

 

 

= ^

^

9 1 ] (x i)

» т- е. каждую функцию N действительных

i=i

\J=1

/

 

 

переменных можно представить в виде суммы 2N+ 1 суперпозизий непрерывных функций одного переменного и суммы N не­ прерывных функций одного переменного. При этом функции ф* не зависят от /, а функции фг, напротив, однозначно определя­ ются функцией /. Частный случай задачи рассмотрен также в- работе [98], где показано, что функция вида (2. 33) представ­ ляется как

N

(2. 35)

г * )= П 1

У=1

 

70

где

 

 

 

 

 

 

 

 

 

 

 

P;= p (2;);

j =

1, W,

 

 

 

если удовлетворяется соотношение

 

 

 

 

 

 

 

dP

dP

_ p

d-P

 

 

 

 

 

d*,-

dzy

 

dz;dzj

 

 

 

являющееся для

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

(2.35)

необходимым

и

доста­

точным.

 

приближенный

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

функции

Ниже дается

Р = Р(*)

в виде,

удобном для решения рассматриваемых

здесь

прикладных задач.

 

 

 

 

 

 

 

С целью получения выражения для показателей Pj надеж­

ности элемента рассмотрим вначале вероятность Р(В).

 

 

 

 

No

 

 

 

 

 

 

Обозначая через Л0= р

Л„

 

 

 

 

 

 

 

 

v= 1

 

 

выходных

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

где No= k\ + k2+ . .. +kN— число всех

N элементов системы, и учитывая,

что

согласно

соотношению

(1. 6 ) /1 0=

Л'„

получаем

 

 

 

 

 

 

 

U А ,

 

 

 

 

 

 

 

 

v - l

 

 

 

 

 

 

 

 

 

Р(В) = Р (А,) Р |Л0)+ [ 1- Р (Л)] Р (В\А0)=

 

 

=

Р (5|Л0) {Р (А . Ж 1- Р (Л0)] Р (£| л 0) р- 1 (£|А,)}

 

пли

 

 

 

 

 

 

 

 

 

Р , в ) - Р т ) { р ( п ; л ) + Р № п в ) р ^ } -

=р(в|л"'К .г-)+ р( 5 х)тйшу!-

где

 

 

 

 

 

 

Л,

П В,

 

откуда

 

 

 

 

 

 

 

г

No

 

 

 

 

 

No

 

 

P(B) =

P(fi|A>)

Р

 

,

 

1

 

 

П

Р(ЯИо) L

 

 

V 1

;

 

\ Х\

) '

v= l

2

р (л ; п л ; - ) + . . . + ( - 1)лг°-1 р ( п л

K J

 

 

 

 

 

 

 

v « l

Р (5 )—Р ( 5 |л 0) |P l,2..... Afo +

^ ? v 0v — 2

 

 

 

 

1

 

 

 

v= l

i < J

+

2 Qijlfiijk —

+

( —

1 ) Л,° _ 1 ^ 1 ,2 ........Aro®l,2.........Aro

i<j<k

71

или

 

Р(В) = Р( В\ А0)

 

Лго

 

 

 

 

 

1-> ]< 7v^lv +

 

 

 

 

 

 

 

 

 

I

v= l

 

 

 

 

'!

j ~

‘ ■+ (—1 VV°

^ 1,2..... JV„T)l,2......A'„j .

(2.36)

Здесь Pi,2,...,w „= P(^ u l— вероятность нахождения

в

допусках

всех N0— ki + li2 +

. . . + k N выходных

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

элементов

системы;

состоящее в

выходе одной (v-й)

характери­

/1 ,— событие,

стики из допуска;

 

 

 

 

 

 

 

 

 

 

 

<7v = P (/ U )—-вероятность выхода одной

(v-й)

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

из установленного допуска (v = 1, N0) ;

 

 

 

 

/-й) ха-

<7,-;=Р (JTi Г) -dj)— вероятность

выхода

двух

(ей

и

рактернстнк из установленного допуска;

 

 

 

 

<hjk=P(Ai П А] П Л ) . •••. <7i ,2.....лг,= Р^ р

 

(Д К )

 

0

р ( д к П А,) ш

 

 

 

 

Р(В|Л0)

4

Р(ДИо)

 

 

 

 

 

 

 

 

 

р

Л'„

 

 

 

 

 

 

 

 

 

 

i=1 <)

 

 

 

 

 

 

 

 

 

 

П

А

 

 

 

 

(2.3Т)

 

 

 

 

Р (B\AV)

 

 

 

 

 

 

 

 

 

 

 

 

 

Лv= I 0v; т1/;-=

1 0/у-; •••;

ть,2, ...,л'0=

1 Oi, 2.....n0\ (2.38)

р(Б|Л,), P ( 5 1Л,• П A j)

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

события

В при

выходе

одной (v-й) характеристики за пределы допуска, двух (£-й, j-й)

характеристик и т. д.

 

в виде

 

 

Представим соотношение (2. 36)

 

 

 

Р (Я) =

Р(Я|Д,) I -

д-

д .

(2.39)

 

V Q . +

Здесь

 

 

i-

1

 

 

 

 

 

 

 

 

Qj = J

Qipn ~ 5 4iv^v ■■• +

( -

1f J

,2,..

,kj 1ll,2.....kj

i =l

 

 

 

 

 

(2.40)

 

 

 

 

 

 

— вероятностная доля в общей вероятности отказа системы1,

обусловленная влиянием только /-го элемента;

 

A = ^ 9 i f t 1'l/ft— 2

<7/»хтЬ‘* ; с + - • • + ( — l)'Vo-1<7i.2,. . . , д ^ 'П гг.........N a

/<ft

1 < * < Х

 

 

 

 

 

 

 

 

 

 

 

(2.41)

72

— вероятностная доля в общей вероятности отказа системы,

обусловленная возможным

«взаимодействием» элементов

(ин­

дексы i, к, х и т.

д. относятся к

характеристикам

различных

элементов).

 

 

 

 

 

 

 

 

Соотношения

(2. 36) и (2. 39)

справедливы для любой струк­

туры элементов системы. Так, в частном

случае при

Tlv =

Tl =

-•• =

1li,2..... лг„= I нарушение любого из условий (2.32)

приводит

к невыполнению задач, возложенных на

систему

(«последова­

тельное соединение»). Тогда только выполнение всех

условий

(2.32)

влечет за

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

Следова­

тельно,

 

 

 

 

 

 

 

 

 

 

 

Р(5|До)=1 и Р ( В ) = 1 - 2 ? / +

2 > у - - - - +

 

 

 

 

 

 

 

 

/= 1

1<;

 

 

 

 

 

 

 

 

 

 

'No

 

 

 

 

 

 

( —

1)ЛЛ°|—V i,2 .........Л'0 = Р

П А ,

 

 

 

 

 

 

 

 

 

 

,i=i

 

 

 

При

независимых

Д,- из

последнего

соотношения

следует,

 

 

Л'о

 

 

 

 

 

 

 

 

что

Р ( В )= П рм где Р; =

Р(А,).

 

 

 

 

;- 1

Вдругом случае, при t],=T],; -=...=rili2 .....Л'„ - |= 0и r]i,2,...iv„=l,

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

характеристик за пределы допусков

означает отказ системы

(«параллельное соединение»). Тогда

по-прежнему Р(£?|До) = 1

 

Л'о

 

и Р ( 5 ) = 1 —<71,2.... n 0 — 1— Р( р|

Д/I- При независимых Д, име-

No

1 =

1

 

 

 

 

■емр ( 5 ) = 1 - П ? / ; ? / = 1 - Р /. i=i

Дальнейшее решение задачи о декомпозиции функции Р (В) должно связываться со структурной схемой и особенностями функционирования системы. Здесь для определенности рассмот­ рения примем следующие допущения, отражающие специфику исследуемых далее задач.

1. Система допусков выбрана таким образом, что при выпол нении условия (2.32) вероятность Р(В|До) превышает требуе­ мое значение Р(В)Т величины Р(В). Пусть Р(В )Т= 0,90, тогда согласно допущению вероятность Р(Д|До) выполнения системой возложенных на нее задач при нахождении всех ее характери­ стик в установленных допусках, должна превышать значение 0,9, т. е. должно выполняться Р |До) >0,9 [например, Р(В|До) = = 0,98—0,99]. Практически это допущение является вполне ра­ зумным. Вместе с тем при решении задачи о распределении тре­

73

буемого значения по элементам системы оно позволяет «норми­ ровать» значение

Р (Д )т

Рт

Р(Д|Л>)

Из такого распределения требований вытекает практическое следствие: данный (конкретный) элемент системы может счи­ таться приемлемым, если все его выходные характеристики ле­ жат в допусках. Если же одна или несколько характеристик вышли из допуска, то для выработки суждения необходимо вы­ числить значение Q, из соотношения (2.38) , а возможно и Р(В), с последующим сравнением их с требуемыми величина­ ми. При этом в зависимости от результатов сравнения могут возникнуть две ситуации: 1 ) данный элемент следует признать непригодным для использования в системе; 2 ) данный элемент, несмотря на выход за пределы допуска отдельных его характе­ ристик, может быть использован в системе.

2.Структурные и функциональные особенности системы та­

ковы, что на основе выражения

(2.39) возможно

представле­

ние

 

 

P(fi) = Pl5|A,) n ( l - Q y ) + 0 [ I - P T]

1.2.42)

J- 1

 

 

P ( 5 ) = P 0 VQ,-;

Р0= Р ( £ |А Д

 

;'-i

 

 

т. e. система близка в некотором смысле к последовательной или к параллельной (структура элемента системы произвольного ви­ да). Конкретный вид зависимости

Plfi) = 'MQy),

(2.43)

вытекает из результатов расчета коэффициентов тщ T},j,.. . на основе статистического моделирования уравнений движения си­ стемы.

Второе допущение, как и первое, не является слишком стес­ няющим. Действительно, пусть, например, система состоит из /V последовательных блоков, каждый из которых включает 1г$. параллельных элементов, т. е. достаточно сложна. Тогда Р(В|Ло) = 1 и соотношения (2.42) оказываются точными.

Изложенное позволяет в качестве искомого выражения для показателя надежности j-го элемента принять значение Pj = = 1 — Qj и осуществить распределение величины Рт по элемен­ там системы, т. е. найти такие величины PTtj, что при Р ;> Р Т, ? справедливо соотношение Р ( В )> Р Т.

7+

При рассмотрении конкретного элемента его номер / в систе­ ме не имеет значения, а соотношение для показателя Pj= 1— Qj надежности элемента может быть записано в виде

Р =

+

..... л-

(2.44)

 

/=1

i<j

 

Здесь

qit qtj, . ..,

qit2..... к — вероятность выхода i-го,

/г-го и

/-го п наконец всех характеристик системы за установленные до­ пуски в течение времени тр работы элемента. Допуски могут ■быть односторонними и двусторонними. В первом случае наибо­ лее общим выражением для <7,- является

 

<7;=

Р Х О *,

t) —аДх, у, z, ?)< 0},

 

(2.45)

где

U,-(•) и а, ( - ) — некоторое

четырехмерное

случайное

поле

(см.

пример 2 .3) и ограничивающее (допустимое) поле;

 

 

х €Е [0,

х 0]; у ЕЕ[0,

*/„];

z ЕЕ[0, z u];

х £Е [0,

тр].

 

В частных случаях

 

 

 

 

 

 

Д =

Р \U,(xt у ) - а , ( х , г / ) < 0 ) ;

 

1

6

 

0/ =

Р

— a ;(t)< 0 );

q, = P {Ut —a, < 0 ) .

J

 

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

зовать следующее достаточно общее выражение:

 

 

?/ =

Р (М *)<£//(•*,

«Лг, T)<a,(t)),

(2.47)

где о, (т)

и bi(x)

— неслучайные функции времени.

когда

Из выражений

(2. 47) легко получить частные случаи,

 

 

<7,-= Р (М*)<£/,(■*, уХ ^ - н )};

 

 

<7/=

Р {Мг)<£7(1;)<адт)];

(2. 48)

<7/=

Р 1Ь-, (т) <

U, < at (т));

q, = Р (b, < Ut < а;)

 

Важно подчеркнуть, что если допуски на характеристики элемен­ та известны, то вычисление величин qi, qц, qit2..... k молено осу­

ществить по данным проектирования и испытаний этого

эле­

мента, не обращаясь к системе в целом

(методы расчета

вели­

чин <7;,

<7;j, .. .,

<71,2,.... л

см. 2.3).

Для

расчета

по формуле

(2.44)

показателя надежности элемента

системы

необходимо

.дополнительно

знать

значения

коэффициентов

гц, ti,7, ...,

• • •, т)1,2,.... ft. которые целесообразно «выдавать»

разработчику

элемента наряду с допусками на

характеристики.

 

 

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

75

Коэффициенты щ, гр,-.......как видно из выражении (2.38), изме­ няются в пределах от нуля до единицы и учитывают, что выход за допуск данной характеристики x t или совокупности характе­ ристик элемента может еще не привести к невыполнению задач системой в целом, вследствие возможной взаимной компенсации отклонений, взаимного «дополнения» характеристик и т. д.

Пусть задана система уравнений движения (функционирова­

ния) рассматриваемой системы н решения £/v=tpv(x,-, г = 1, No) уравнений (No— число выходных характеристик всех элементов системы) должны удовлетворять требованиям (2.34). Наруше­ ние этих требований расценивается как отказ системы. Рассмот­ рим вначале метод определения вероятности Р(В|Ло) по резуль­ татам статистического моделирования процесса функционирова­ ния системы в целом.

Для осуществления такого моделирования необходимо знать функции плотности вероятности величин л',-. Здесь имеется ряд особенностей, на которых целесообразно остановиться. Пусть X) случайная величина с функцией плотности вероятности /о(х,-), когда нет неисправности в системе, влияющей на параметры или вид/(л',). Это состояние встречается с вероятностью Р„= 1—

,ьр

— | f u(y)dy, где /н(у) — функция плотности вероятности вре­

о

мени до возникновения неисправности.

Из формулы полной вероятности следует, что с учетом неис­ правности f(Xj) будет выражено так:

/(-*/) = Рн/о(Л'Л + \ Pf H(y)f{Xi\y)dy,

(2.49)

о

 

где / (-V,-1//)—функция плотности вероятности случайной величи­ ны при возникновении неисправности в момент г/е[0, тр], а тр — время работы элемента в системе.

Пусть, например,

и известны функции цДт), аг(т), где т — момент возникновения неисправности. Тогда

f (•*/)=

1 -\f2л о

76

В наиболее простом виде выражение (2. 49) удается предста­ вить в том случае, если неисправность независимо от времени ее

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

/(х,-). Пусть уста­

новлено,

что npii_ возникновении

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

/ (х*) = f* (.v,-) =

=/(х*|Ли), где An — событие,

состоящее

в

возникновении не­

исправности. Тогда

 

 

 

 

 

 

7 (* ,) = р„/о (•*,-)+(1 - р„) г

№ ) = р ш

/

(* ,и „ )+

 

 

+ Р ( А . ) / ( Л';Ин)-

 

 

(2.50)

Если таких неисправностей несколько, то с помощью соотно­

шения (1.31) легко показать, что

 

 

 

 

?(■*/)= У Р ( А „ )/(* ,И „ ,) -

у

Р ( лнч П Ли]) f

{х;\Ан,[\Лн])+

 

 

 

•><]

 

 

 

 

+ ... +

( - l K P ( n ^ „ vj / ^

 

 

 

 

+

P^ П

 

П

) 1

 

 

 

 

 

 

v -1

 

 

 

где

_ с — число неисправностей;

 

 

 

/T„v — событие, состоящее в возникновении v-й неисправ­

 

ности;

 

возникновении А.

f ( Xj \ A) — функция /(х,-) при

Если неисправности влияют на ху в различные моменты вре­ мени различно, то

стр

/(■*/) = 2 \ /» (?« )/(JCill/’)flf?*_ v=l 0

“ 2 f P J P / н Ci/v, y j ) f { X i \ y „ y j ) d y vd yj +

•i<jb о

+ ••• + (

f(Ul,

1/2. ••■. yc)f {Xitel, 1/2.•••y,)U dy^ +

 

 

 

 

 

v -1

 

 

+ /o (*() j

• •

j / н {уI,

Уа, • • • . Ус) П

dlJ->

(2.51)

 

 

 

 

V= 1

 

где /и(г/,,

г/j),.. ., fn(yi,

уг, ■■■, Ус)

•— совместные функции

плот­

ности вероятности времени до возникновения

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

(v-й и /-й)

и т. д.; fo{Xi)

— функция f(Xi) при

отсутствии

неис­

правностей.

При вычислении коэффициентов тр, гр-р ... следует использо­ вать функции плотности вероятности /(х,-), найденные с учетом

77

_ —>го

неисправностей. Дополнительно потребуется f (Х() множить на коэффициенты усечения Д,-. Будем говорить о наличии усече- н и я п е р в о г о т ип а, если этот коэффициент имеет вид

—х

Если коэффициент определяется как

 

—1

\ f\.x,)dx,

ПЛИ Е ;= \ f ( x , ) d X i

_ а

то будем говорить об у с е ч е н и и в т о р о г о т и п а . Если Д ,= = 1, то функцию J(Xi)Ei будем называть неусеченной. Пусть функции f (xfiEj для всех аргументов а',-, влияющих на выходные характеристики системы, известны. Тогда в соответствии с мето­ дом статистических испытании можно, используя иеусечеиные функции /;, найти вероятность Р (В) выполнения системой в це­ лом поставленной задачи:

 

 

 

 

 

 

 

( 7 ( У , ) П <1У1= ,

 

1^- 52)

где

f(yc) — совместная

функция

плотности

вероятности (yh

 

 

 

 

У* ■

Уд;

 

реализаций при

расчете

функции

 

п и т — общее число

 

 

 

 

y4 — y { X tvj,

A-jVo = (A'i . . . А'д-J, /= 1 , Р И ЧИСЛО СЛу-

 

 

 

 

чаев, когда условия (2.34) при моделировании вы­

 

 

 

 

полнялись.

 

 

 

 

 

 

 

Число реализаций п выбирается из условия, чтобы довери­

тельный интервал [Р, Р]

для вероятности Р = Р(В)

был не шире

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

 

 

 

 

 

— - н 2 < Р ( Я ) <

— + ех,

 

 

(2.53)

 

 

 

 

 

п

 

 

п

 

 

 

т-по

-

=

т

п

К

 

т

 

 

 

 

где

 

------ н,

г, = Р---------- заданные значения;

 

 

 

 

 

/;

-

 

 

п

 

 

 

 

Р и Р — верхняя

и нижняя

границы доверительного

интервала

для Р (В), определяемые

из соотношения (1.

174),

где следует

положить А'Гф= я т.

 

 

 

 

 

 

 

Для вычислений искомой вероятности Р (В/А0), входящей в

выражения для коэффициентов rp, rp-j, . . ., гр, 2..... и,

очевидно,

достаточно,

повторяя изложенную

процедуру

при

моделирова­

нии, воспользоваться функциями плотности вероятности /'(*,■)£,• с усечением первого типа.

78

Для нахождения вероятности Р(В|Л v) по смыслу ее опреде­ ления необходимо для моделировании случайной величины использовать усечение второго типа. Коэффициент усечения при этом берется из условия «выведения» за пределы допуска, т. е. за ту его границу a v или 6 V, выход за которую наиболее неблаго­ приятен для выполнения задачи системой в целом. В случае, когда заранее неясно, выход за какую границу допуска является более опасным, можно приближенно принять ту, при которой ко­ эффициент г| v оказывается большим. При моделировании осталь­ ных величин Xi (1фх) в случае вычисления1!; следует восполь­ зоваться иеусеченными плотностями.

Аналогично изложенному находятся вероятности

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

Следует подчеркнуть также такие три «привлекательных» стороны показателя (2.44).

1. Известно, что допуски на характеристики в ряде случаев приходится задавать приближенно, то излишне расширяя, то не­ сколько сжимая поле допуска. Это происходит вследствие сложности разработки рациональной системы допусков. Пред­ положим, что такая система в принципе существует, а фактиче­ ские поля допусков заужены. Тогда значения вероятности qi, qij,. . . выхода характеристик из допусков оказываются завы­ шенными, но уменьшаются коэффициенты rp, rp-j,. . ., так как ве­ роятности P(fi|A;), P(5|.4,j) увеличиваются. Если допуски рас­ ширены по сравнению с рациональными, то величины q{, qu, ...

уменьшаются, но возрастают тр, трр,.. . Следовательно, в некото­ рой области отклонения совокупности границ допусков от их ра­ циональных значений произведения ^тр> q^'rp-j,. .., входящие в соотношение (2.44), могут сохранить неизменные значения по сравнению со случаем, когда имелась бы рациональная система допусков. Это позволяет построить систему контроля по совокуп­ ности параметров в сочетании с принятой системой контроля по допускам для отдельных характеристик с приемлемым значени­ ем риска поставщика.

2 . Использование смешанных плотностей распределения вида (2. 49) — (2. 51) при определении q q ^ и rp, rp-j позволяет учесть влияние неисправностей элементов на выходные характеристики системы. Заметим, что деление всех нежелательных явлений на неисправность (любое отступление от технической документа­ ции [60]) и отказ — «сильная форма неисправности», при кото­ рой элемент выходит из строя, является удобной, но упрощеи-

79

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