Файл: Волков Е.Б. Основы теории надежности ракетных двигателей.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 11.04.2024
Просмотров: 238
Скачиваний: 0
Р Г = Р^П Л^-jc помощью формулы (2.87). При этом по оси z
следует взять две точки z = 0 |
и |
z = 6, |
так |
как |
случайная |
|
функция и (г)=7д(г)—Т (z, х) |
монотонна на |
интервале |
[0, б], |
|||
а при z = 6 имеет разрыв первого |
рода. |
Следовательно, |
7Vi=2. |
|||
По осп л- можно взять одну точку х = л'цр, |
соответствующую кри |
|||||
тическому сечению сопла, если |
считать, что величина Тд по х |
|||||
является одинаковой во всех сечениях и что Т(х, г) |
монотонно |
возрастает в интервале [0, .vKp], а затем монотонно убывает. Сле
довательно, в первом приближении |
N\ = 2, N2 = 1 и N = N\N2 — 2, |
|||
•а задача определения Рг состоит в вычислении вероятностей |
||||
Рх= Р 1Рд(0)-Р[л-кр, |
0) > |
0) = Р (иТ1> 0); |
|
|
. Ра = Р{Рл(8 )-7 ’К-р, 8 ) > 0 ) = Р ( « „ > 0 ) , |
|
|||
■где |
|
|
ТЛ Ь ) - Т ( х кр, 8), |
|
ип = ТЛ 0 ) - Т ( х кр, 0); |
ит, = |
|
||
.и коэффициента корреляции qu u |
. |
Таким образом, |
искомое |
|
значение Рт определится по формуле |
|
|
|
|
P r = p i P 2+ ( P m - P i P 3) |
a resin QuT 1ur 2 - |
(5 - 13) |
Пусть случайная величина и-п имеет нормальное распределение. Тогда вероятность надежного охлаждения в /-м сечении
|
|
|
|
|
|
|
Pi = P{iiTi> °)= f 'Uli)- |
|
|
(5.14) |
|||||||
Здесь F(hi) — функция Лапласа; |
|
|
|
|
|
|
|
||||||||||
|
|
|
//, |
|
1 |
г . |
i |
т. |
|
‘Пт. ~ 1 |
при f= 1,2, |
(5.15) |
|||||
|
|
|
|
~2 |
I |
|
|
|
9 |
||||||||
|
|
|
|
|
|
1 Д |
I |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
G j . |
|
, |
1I |
|
|
|
+ LC |
|
|
|
|
|
|
|
|
|
|
1ЛL |
|
|
|
|
|
|
|
|
|
||
где |
|
|
|
|
|
|
7 \ |
и а*2 — среднее значение и дисперсия ве- |
|||||||||
|
|
|
|
|
|
|
_ |
_ |
личины Г,-; |
|
|
запаса |
по тем |
||||
|
|
|
|
|
v\Ti = TRi/Ti — среднее |
значение |
|||||||||||
|
|
|
|
|
|
|
|
|
|
пературе; |
|
|
|
|
|||
-vT — cr |
|
T |
|
и v T. = 5Т. Тi ~ коэффициенты |
вариации |
случай- |
|||||||||||
Для |
|
определения |
_ |
|
ных величин TAi и 7Y |
|
|||||||||||
|
T A i и |
7 \ |
|
воспользуемся |
выражениями |
||||||||||||
(5.10), (5.11), (1.72) и (1.73). Тогда |
|
|
|
|
|
||||||||||||
Т |
- |
— |
Г |
|
} |
|
|
|
|
|
|
Т |
- = |
Т |
I |
|
|
1 Д |
1 ------ 1 г.д |
|
при z е= [0, |
6], |
|
1 А 1 |
1 Х .Д |
|
Ъ, |
||||||||
о |
■ |
2 |
|
|
|
2 |
|
2 |
1 при г = |
||||||||
|
|
|
|
|
|
|
|
|
|
от |
= 3 Т |
|
|
|
|||
д |
|
= атп„ |
|
|
|
|
|
|
|
|
1 |
|
I |
|
|
||
|
|
|
г . д |
|
|
|
|
|
|
|
Д |
Х . Д |
|
|
|||
т |
|
Д |
х . а |
и 4 |
, |
°гх — средние |
значения и дисперсии вели- |
||||||||||
1 Г . Д ’ |
|
|
|
|
|
|
|
используемых |
|||||||||
тГ . д |
|
и |
Тх.л |
■— |
определяются |
из соотношений, |
|||||||||||
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
182
в прочностных расчетах, а также из справочных или экспери ментальных данных. Найдем величины Ti и ^..соответствующие
фиксированным значениям х и 2 (обозначим их через х,-, г,-). Точ нее, I — здесь номер пары значений х и z; так, если N2=l, то такими парами будут (Х|, z\) и (хь z2). Из приведенных соотно шений получаем
|
T l -----T ( X h |
Z ; ) |
— 7’Х-П;-Ь Qs;(B— zn |
|
|
|
||||||
От/ ~ От |
L T |
!K |
+ |
( t )'-! + т |
< * - ’ ■> |
"X |
||||||
1 x.n i |
V- |
|||||||||||
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Qsi |
|
1 |
|
|
5 |
1 |
|
|
|
|
|
|
|
|
|
+ |
|
|
|
|
|||
|
|
|
|
(2г 1* |
Г |
+ - |
|
|
|
|
||
|
|
|
|
|
X |
Иж |
|
|
|
|
||
|
|
|
|
1 |
|
|
|
( ° 2г Г |. + °гж/)+ |
|
|
||
|
|
|
a |
|
1 |
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
I |
|
|
а ж1 |
|
|
|
|
|
|
|
|
|
oil |
|
*а.г / |
| |
| |
о 2 |
| |
1 |
_ 2 |
|
|
|
|
|
Г4 |
I |
T l T |
г4 |
1 -4 |
Ja>K/ |
|||
|
|
|
|
|
|
/ |
|
Л |
Л |
C |
t . ; |
|
^ r i — ^ 0 |
+ s rSrQ o ^ r |
a r г; |
|
Г ж ,■— |
Т’вх + |
V,-; |
|
|
o , , = S y ; 4 |
|
Ur |
\ |
r |
|
( » ; . + < + < + 1 |
||||
1 1 ж / |
1 |
|
|
|
|
|
' |
|||
2 |
вх + |
< ; |
a ^ . = |
ar i (^?p + |
0,67T;2- + 3 > 2D); |
|||||
Ж X |
вх |
|
|
|
|
|
|
|
|
|
3a.|; |
|
аж i |
(^СрЖТ 0-.04z),)|_-{- 0 ,0 4 1>оэ . -|- |
|
||||||
|
+ 0 ,6 4 ^ ж , + |
0,64г& ; + |
- 1 |
4 r ) |
; |
|
|
|||
|
|
у,. = я |
|
Qsi D d x |
|
|
|
|
||
|
|
|
|
|
СжП1Ж |
|
|
|
|
|
v = x/ax — коэффициент |
вариации |
х; |
(общее обозначение). |
|||||||
Величины, входящие в приведенные |
соотношения, |
содержат |
||||||||
|
__ |
— |
|
и дисперсии |
0 |
9 |
■ |
соответст |
||
средние значения То, ar ;,... |
oj-o, |
\ |
||||||||
вующих аргументов. Некоторые из них (б, Пж |
5 Ж,-) |
являются |
185
«первичными», т. е. аналитически не выражаются с помощью других аргументов, а находятся непосредственно по справочным, экспериментальным данным п данным по технологическим допу скам. Другие же, в свою очередь, представляют функции слу
чайных аргументов (Го, т, icy).
Приведенные соотношения (5. 14) и (5. 15) позволяют решить задачу определения вероятности Р,-.
Для вычисления коэффициента корреляции q£j- между uTi и мт} используем соотношение (4.9), из которого следует, что
|
*т.ат- |
|
|
J , . |
- а - |
|
|
|
|
°Т-°Т ■ |
|
|
Г' Т] п |
I |
7 д I Г л / |
|
|
|
|
|
|||
е,г |
|
|
оГ д |
|
|
|
' i 1л I |
б7-;Г |
|||
t .-g ; |
QTiTj “Г |
|
|
i |
T !I |
j |
а,оу |
||||
|
‘ |
J |
|
aiaj |
|
|
л J’ |
||||
|
Ji u J |
|
|
|
|
|
|
|
|
|
а7. равно сгтj при замене индекса / на у:
г |
~ |
3гЛ(. |
3j — |/ s j - . - f о]Д .—2зг^ ДА ; г д г |
з, = 4 |
37';- |
Упрощая задачу, в первом приближении можно считать, что случайные величины Г,- и Гд,-; Г,- п Гд ц Г,- и Гд Tj и Гд ,■незави симы. Тогда
|
б ы 7- . Н у . ' |
О;: |
|
ar°TrTiTj |
: |
1 |
Д 7 |
1f л; |
|
(5. 16) |
|
|
' i ' 1 |
|
|
|
1/4,. + |
4 |
Д(.К |
|
|
||
|
|
|
|
|
|
|
|
+ °Г |
|
|
|
Для |
определения величин |
|
дг.Г;. н дГ(.г,(;. |
используем |
соотноше |
||||||
ние |
(1.80) п выражения |
для Г и Гд. |
Из |
соотношения (1.80) |
|||||||
видно, что коэффициенты корреляции |
п-.^. и gr |
j.^. |
находятся |
||||||||
в функциональной |
зависимости от производных по аргументам, |
||||||||||
общим для Г,- и Ту, |
Г,, |
и ГЛ;.. |
|
|
|
|
|
|
|||
Кроме того, |
|
|
|
|
|
|
|
|
|
|
|
|
Qr,Tj '■ |
о_ g |
(а13Г „ 'г / 35+ - ' •)’ |
|
|
||||||
где |
|
|
Г/ |
7-/ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
а, |
|
|
|
|
|
|
|
_L |
JL |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
аг;. |
X |
ж; |
£Zi>—
т. е. каждый коэффициент есть произведение производных от Г£ и Tj по соответствующему общему аргументу. Производные берутся в точке, где все аргументы равны своим средним значе
ниям (это подчеркивается введением индекса «0»), Таким обра зом,
184
0., |
„ |
|
( а I <4 |
fl 5; + . . . ) t 5 г а7- |
От- |
т + |
01 Оо |
(5. 17)! |
|||
|
= > |
7„°___11= — |
д| |
'гД= |
- -r -'-f4— |
|
|||||
“г, |
|
г2 |
|
|
|
|
I 1+ °1 1 |
1 Оо |
|
||
где |
|
|
1/'аГ, + аГ I ]/"°Г,+°Г., |
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
у-,./3?',: ^2 — Зг о/3?'г! |
Зг .] > 0 ; |
= Г т > 0 - |
|
||||
Подставляя выражения |
(5.15) и (5.17) в зависимости |
(5.14) |
|||||||||
и (5. 13), |
находим искомое значение Рт = Р(т4|). |
|
в вы |
||||||||
Остановимся на рассмотрении события А2, входящего |
|||||||||||
ражение |
(5. 8) |
и состоящего в выполнении условий по прочности: |
|||||||||
и устойчивости для случая, когда событие А\ |
выполняется. В ра |
||||||||||
ботах [2, |
83, 27] |
детально исследованы |
методы расчета |
камеры- |
|||||||
сгорания ЖРД |
на прочность, устойчивость и колебания |
и уста |
новлена необходимость при проведении таких расчетов рассмот рения нескольких расчетных сечений камеры, а в общем слу чае— всей конструкции камеры, времени работы и эксплуатациидвигателя. Следовательно, в выражении (5.8) щ и и2 представ ляют собой случайные поля четырех переменных u\ = ii\(x,y,z,x)\. и2 = и2(х, у, z, т), где х, у, z — координаты, т — время. Введение трех координат обусловлено тем, что оболочка камеры как пра вило является двуслойной.
В труде [83] показано, что при расчете на прочность в качестве-
условия успешного |
функционирования камеры целесообразна |
||
принять следующее: |
|
|
|
|
•п = Е Ш > |
1 |
|
или |
. Рл |
|
|
|
х ) > р л(х, т), |
(5. 18) |
|
Pp= P{aR) = fp{x , У, z, |
|||
где рл= p{AR) = рр(х, |
у, z, т) — допустимое значение |
давления, |
вычисляемое в функции от допустимого увеличения радиуса ка меры (от допускаемой деформации) или в функции от коорди нат и времени; рл {х, т) — давление в камере двигателя в функ ции от длины камеры и времени работы двигателя.
Используя полученные в работах [2, 7, 83 и др.] соотношення-
для |
определения рр, можно найти значения рр для любых х, у, г |
и х. |
Тогда аналогично изложенному выше находим вероятность- |
события А 2\, |
состоящего в выполнении условия |
«,>0 |
по проч |
||
ности: |
|
|
|
|
|
Рп = Р ( Ж 31) = |
Р ;/7р (.х-, у , Z, |
х) — рл(х,х) = и1(х, |
У , Z , |
г ) > 0 } |
= |
|
= |
Р ^ П |
|
(5. |
19) |
где А2н= {«1 *>0}— событие, состоящее в выполнении условия по прочности для фиксированной «чет верки» дискретных значений (,v,-, у,-, zit т,-);
185
v — номер комбинации (л*,-, yit Zj, п ) ;
Л/'1 = N i N 2JVaN .i— общее число комбинаций |
(общее число то |
|
чек) ; |
|
|
JVi, N 2 , N 3, ЛД — числа |
точек на осях х, у, z и т соответст |
|
венно. |
|
|
Далее для определения Рп используется изложенный уже ме |
||
тод [см. соотношения (5.6), |
(2.86)—(2.90)]. |
|
Расчет на устойчивость |
камеры сгорания |
обычно сводится |
к рассмотрению устойчивости неохлаждаемого |
насадка, нагру |
женного на нерасчетном режиме работы внешним перепадом
давления. Порядок приближенного расчета вероятности |
Ру = |
||
= Р(»2> 0 ) |
аналогичен приведенному выше и приводит к |
соот |
|
ношению |
|
|
|
|
Ру = Р(„а> 0 ) = Р (Л аа)= р ( п |
Лаа,) , |
(5.20) |
где событие |
A2ii—{u2i'^>0}i=], N 2 и N 2 |
— число точек, рас |
сматриваемых при вычислении вероятности P (h2>0). Следова тельно, соотношение (5. 8 ) может быть записано в виде
Рпр= Р(..43)= Р(Л 21 П Л22)= |
Л', |
/Vа |
(5.21) |
||
Р П |
Ай11 п л 2 21 |
||||
|
|
|
;=1 |
|
|
Таким образом, из соотношений |
(5.19), (5.20) и равенства |
||||
Л = Л] Р) А2 получаем следующее |
выражение для показателя на |
||||
дежности камеры сгорания: |
|
|
|
|
|
PM= p H i П Л3) = Р |
|
|
Аг, |
N2 |
|
1 |
Ли П |
Л,- ГП Л |
|
||
|
|
/=i |
1=1 |
|
или
(5. 22)
Поскольку события Л2н и Л22; уже определены выше как соот ветствующие условию успешного охлаждения (пример 2.5), то выражение (5. 22) может быть переписано в виде
Рк.с = Р ( П |
Р ( n А2 |
п А2з,- j = РГР, |
(5. 23) |
\ / =1 |
/ 'i=i |
/=1 |
|
где Рг — величина, определяемая из формулы (5. 13). Согласно выражению (2.86)
N
Рк.= = Рт-РпР = |П Р ( И /> 0 ) + Рт- П p f c > 0 ) K N X /=1
186