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

 

 

ЮТ