Файл: Волков Е.Б. Основы теории надежности ракетных двигателей.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 11.04.2024

Просмотров: 214

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

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

Р (Д )т

Рт

Р(Д|Л>)

Из такого распределения требований вытекает практическое следствие: данный (конкретный) элемент системы может счи­ таться приемлемым, если все его выходные характеристики ле­ жат в допусках. Если же одна или несколько характеристик вышли из допуска, то для выработки суждения необходимо вы­ числить значение 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