Файл: Волков Е.Б. Основы теории надежности ракетных двигателей.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 11.04.2024
Просмотров: 221
Скачиваний: 0
л'2 u 1/2 — действующее значение изгибающего момента п его допустимое
значение и т. д.
При этом случайные величины х,- и (/,■ не зависят от времени работы рас сматриваемой системы н распределены нормально со средними значениями .Vf
и tji и средними квадратическими |
отклонениями |
а,. , |
ац |
Коэффициенты |
||||||||
корреляции величин х,- и |
|
(/,■ равны |
гх .у.. |
По данным прочностных расчетов |
||||||||
с использованием формул (1.72) — (1.80) найдены: |
|
|
|
|
|
|||||||
n = In — xi’> |
|
|
+ зу. — 2ад-/°г//о-,•</,•; |
|
>ч — |
|
on- |
|||||
Здесь Qa — коэффициент |
корреляции величин г/,- |
и и, |
согласно выражению |
|||||||||
(I. 109), определяемый как |
|
|
|
|
|
|
|
|
|
|||
Qij ~ |
°y;aUj |
гУ;У 1+ |
°Л-г°.Г;- |
Cv.-.v- — |
аУ-°Х] |
ГУ[X ; — |
|
|
||||
о;ау |
а,-а? |
аj |
ajdj |
ГУ jX!‘ |
||||||||
|
' |
] |
|
' 1 |
‘ |
J |
|
J ' |
Требуется вычислить показатель надежности конструкции, представляемый в
виде |
вероятности |
Pi, 2 ■• • s=P(x'i<i/,b |
х2<</2, • • •; |
|
Xs<tjs), |
если |
/ti=/i2='2; |
|||||||
Л3=/г4=3;' /г5=1,5; Л6=/г7=1,8; Л3=1,9, а корреляционная |
матрица |
II, со |
||||||||||||
ставленная из расчетных значений д,-j, имеет вид |
|
|
|
|
|
|||||||||
|
|
|
|
0,9 |
0,8 |
0,5 |
0,0 |
0,6 |
0,7 |
0,8 |
|
|
||
|
|
|
|
1 |
0,7 |
0,6 |
0,5 |
0,0 |
0,8 |
0,9 |
|
|
||
|
|
|
|
|
1 |
0,5 |
0,6 |
|
0,7 |
0,8 |
0,9 |
|
|
|
|
|
lle/yll |
|
|
|
1 |
0,7 |
|
0,8 |
0,8 |
0,8 |
|
|
|
|
|
|
|
|
|
1 |
|
0,7 |
0,7 |
0,6 |
|
|
||
|
|
|
|
|
|
|
|
|
1 |
0,5 |
0,6 |
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
0, 8 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
Решение. |
Находим величины |
|
|
|
|
|
|
|
|
|
||||
|
Р,- = F (Л,-); |
Pi = |
Р2 = |
(2) = 0,977; |
Р3 = Р4 = F (3) = |
0,999; |
||||||||
Р5 = |
F( 1,5) = 0,933; |
Р6 = |
Р7 = F (1,8) = |
0,964; |
Р8 = |
F (1,9) = |
0,974; |
|||||||
Рт = 0,933; |
Д |
|
|
|
2-2 |
|
|
|
|
-j- arcsin 0,8 + |
...) = |
|||
П Р / = |
0,8о1; К — ------------ (arcsin 0,9 |
|||||||||||||
|
|
1=1 |
|
|
3,14-7-8 |
|
|
|
|
|
|
|
||
|
|
|
|
|
= |
0,472. |
|
|
|
|
|
|
|
|
Следовательно, |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
Pi, 2 , |
s=0,851 + (0,933 — 0,851) • 0,472= |
0,890. |
|
|
||||||||
Пусть в условиях iii> 0 все х, |
равны |
(хт= |
х2= .. . — xN=x)> |
т. е. нагрузка одна и та же, а величины г/,-, несущие способности: по отношению к нагрузке х, различны и независимы. Тогда ко
эффициент корреляции между |
и |
как видно из соотноше |
ния (1. 109), равен |
|
|
2 |
|
1 |
сх |
|
|
е,у |
|
|
102
где |
и Зд. — дисперсии несущей способности г/,- и нагрузки х. |
Если дисперсия о2, одинакона V /^[1, А^], то
_ 1
^1+ (o;//a.v)2
Следовательно в соотношении (2.91) для вычисления искомой вероятности Р(ы,->0, v i — \,N) величина
K N = nN {N — 1) У |
|
|
4 |
|
|
|
aresin |
+ |
4 |
' |
3!// Ф aUj' |
||
|
^<J |
V 4 |
, 1 / 4 |
|
||
|
|
|
|
|
+ °, |
|
или |
К N ■ ■-— arcsin [1 +(зу/зЛ.)2]-1, |
если |
V ( ,y \ . = |
5 , , = v |
||
|
тt |
|
|
|
* |
J |
|
Интересно отметить, что при a2v^ > 4 ; (дисперсия нагрузки су |
|||||
щественно больше |
дисперсии |
несущей |
способности) имеем |
Qi}— >-1 и Р (п;>О, V / = 1, N)— >-Рт , т. е. вероятность одновре менного выполнения условий w*>0 оказывается близкой к ми нимальной из вероятностей Р,-.
В рассмотренных выше примерах для успешного функцио нирования системы было необходимо, чтобы все компоненты век
тора uN= ( u \ . . . u N) |
были больше нуля. Но во всех этих |
слу |
|
чаях предполагается, |
что стоит одному из |
условий w,->0 (или |
|
л,-<ц{< 6,-) нарушиться (т. е. реализуется |
событие Л,-), как |
си |
стема выйдет из строя. Однако это не всегда выполняется и осо бенно, если система входит в более сложную систему. Наиболее общее из рассматриваемых нами выражений для показателя на дежности является выражение,(2.44):
N |
|
|
|
р = 1 - 2 т + 2 w u - 2 |
+ |
||
<-i |
i<j |
i<j<k |
|
“Ь ■• • “Ь( — 1)
Используя приведенные соотношения, выпишем для общего слу
чая выражения для величин qu q^, . . ., q |
2..... к'- |
||||
|
<h = P(A,)\ |
|
|
||
<7,у = q f l j + |
(qm - |
qflj) — a rcsin |
б л ; a j + s2; |
||
Qijk q i q f l k + { q m |
Я'Д |
|
з“Ьгз’ |
||
|
N |
|
N |
^ |
- |
|
П <h |
Qm |
|||
<71,21...1лг = |
П qi |
к n ~\~ |
|||
|
i«i |
|
/=1 |
|
|
103
Здесь
K n |
4 |
" V |
arcsin р , |
, |
|
||||
|
n N ( N — 1) jtU |
1 3 |
||
или |
|
<</ |
|
|
_nW(W- 1) _^ |
|
_ |
||
Л лг |
arcsin q,;-, |
|||
L |
l<]
если
Pl4,) = P (a, < « ,< & ,) ,
где q,-j — коэффициенты корреляции между «» и iif,
qm— минимальное из значений qit входящих в выражения
ДЛЯ qiit qijhi • - •> Я\, 2......... |
дг- |
Вычисление величин <7,-, </{,■,. . входящих в выражение (2.44), в ряде случаев сводится к определению вероятности непревышения случайной функцией или случайным полем некото рого уровня [см. соотношения (2. 12) — (2. 18)].
В настоящее время выполнен ряд фундаментальных исследо ваний [5, 46] задачи по определению вероятности
Р 1Тр)= Р [к (т)> 0], где t e [ 0 , тр]. |
(2.'93) |
Исходными допущениями для задачи в работе [5, 46] являются непрерывность реализации и;(т) случайной функции и (г) и дважды дифференцируемость корреляционной функции р(Дт) при Дт— »-0 (в стационарном случае). Основное внимание уде ляется при этом случаю, когда н(т) — нормальная стационарная случайная функция. Типовые допущения и метод решения зада чи в принятой постановке рассмотрим, распространяя на неста ционарный случай решение задачи, данное в работе [94].
Пусть Д — весьма малый интервал времени. Вероятность Р(т, Д) выброса за время Д функции и(т) за линию и(т ) = 0 представим в виде
Р(Т, А)= Р (А п £ ) = Р ( Д ) - Р ( Л П В),
где А — событие, состоящее в невозникновении выброса в мо
мент времени т; В — событие, состоящее в возникновении вы броса -в весьма малом интервале Д, следующем за моментом т. Последнее соотношение перепишем в виде
Р(т, д ) = | <a\u{T)\d [и (т)] —
о
00 со
— f ('?[«(*)> к(г + д)]г/[«(т)]£/[д(т-|-д)].
6 6
104
Здесь гр[д(т)] — плотность распределения и(т) в данный момент т; ф[ы(т), г/(т + А)]— плотность совместного распределения и(т) в моменты т п т+Д.
Допустим, что выбросы образуют нестационарный пуассонов ский поток, а и(т ) — непрерывна. Тогда согласно теореме рабо ты [87] параметр потока выбросов определяется как
Х(т)= lim Р(+ А)
д-о л
ИЛИ
l(j
|
— ? |
f «р[и(т), и ("£-}- л)] d \и (т)] d [w(t-f д)]| |
(2.94) |
||||
|
6 |
о |
|
|
|
Ja- o |
|
и, следовательно, |
|
|
|
|
|
||
|
|
Р(тр) = ехр | — |
f P А(т)Дг|, |
(2.95) |
|||
|
|
|
|
I |
6 |
J |
|
где А(т) |
находится из соотношения (2. 94). |
|
|
||||
Для |
случая, когда |
м(т) — нормальная |
случайная |
функция, |
|||
из соотношения (2. 94) |
получаем |
|
|
|
|||
|
|
Х ( т ) = — ( Д [Л (г)] — Р ь2)д-*0- |
|
||||
Здесь |
|
|
А(т) = «(т)/з(т); |
|
|
||
гг(т) и |
а ( г ) — среднее |
значение |
и среднее квадратическое от |
||||
клонение и (т) |
в данный момент т; |
|
|
|
|||
Pi.2= ^ |
(*i) F (Ла)+ | |
ф Л - 2, z)dz |
[см. 1.104)], |
||||
|
|
|
6 |
|
|
|
|
где |
|
//-!= |
/г (г); |
Л2 = |
Л (т-|-д); |
|
|
р = д(т, |
Д )—-коэффициент |
корреляции |
между |
значениями |
и(т) в данный момент времени т и в момент времени т+А.
Сучетом известного правила дифференцирования интеграла из соотношения (2. 95) можно получить выражение для расчета
Л,(т), которое, к сожалению, оказывается достаточно сложным. Частный случай (2. 94) рассмотрен в труде [94] и относит ся к задаче, исследованной ранее [46]. В этой задаче предпола
гается, что гДт)— непрерывная стационарная нормальная слу чайная функция, а корреляционная функция дважды дифферен
105
цируема в нуле. Кроме того, считается, что поток выбросов яв ляется пуассоновским и стационарным. При этих предположе ниях из соотношения (2. 95) получено
|
X(т) —Х = const = |
j ^ / " С„ ©(//.) |
(2.96) |
и |
Р(тр)= |
е_ ь р, |
(2.97) |
где |
Со — среднее число пересечений реализациями |
«,-(т) сред |
|
него значения й функции н(т); |
/г— й/а, а й н а — среднее зна |
||
чение и среднее квадратическое отклонение и(т); |
|
?1Л)=12л) 2 expf - - ^- j .
В практических приложениях соотношения типа (2.97) в (2.95) (помимо сложности последнего) вследствие большого числа принятых допущений оказываются не всегда применимы
ми. Так, выборочные реализации «,•(т) могут иметь |
разрывы, |
различную длительность, а корреляционная функция |
q ( A t ) мо |
жет оказаться недпфферепцируемой при Дт = 0. |
постановку |
В связи с этим рассмотрим другую возможную |
задачи, позволяющую в ряде случаев получать приближенное аналитическое решение с меньшим числом ограничивающих до пущений.
Пусть и(т) нестационарная случайная функция с возможны ми разрывами (первого рода) в отдельные моменты времени.
Выберем сечения T,- = const (/=1, N) таким образом, чтобы к интервалах [т,-, т,-+|] случайная функция и(т) могла быть аппрок симирована функцией, заданной на монотонных реализациях [это и есть правило С^из соотношения (1. 124) в данном случае].
Тогда, принимая основное допущение о том, что вид функции и(т) позволяет осуществить такой выбор сечений, обозначим че
рез Ai |
событие, состоящее в непересечении и{х) границы и(т ) — |
||||||||
= 0 в /-м сечении. Представим Р(тр) в виде |
|
|
|
||||||
|
|
P(tp)= P |
N |
\ |
N |
|
\ |
|
|
|
|
П Л,\Р |
П |
АЛ, |
|
||||
|
N |
1=1 |
')р(в /=1 |
|
|
|
|||
где |
А;— событие, |
состоящее |
в |
выполнении |
условия |
||||
А = П |
|||||||||
|
/=1 |
|
|
|
|
|
|
|
|
и (т;) = И;> 0 |
во всех N выбранных сечениях |
(/=1, А); |
В — со |
||||||
бытие, состоящее в том, что условие |
ц(т)>0 |
будет выполнено |
|||||||
внутри выбранных интервалов, т. е. |
|
|
|
|
|||||
|
|
v t e ( t „ |
т/+1), i = |
1, N. |
|
|
ЮТ