Файл: Волков Е.Б. Основы теории надежности ракетных двигателей.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