Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Волков Е.Б. Основы теории надежности ракетных двигателей

.pdf
Скачиваний:
55
Добавлен:
25.10.2023
Размер:
12.65 Mб
Скачать

 

Р(0) =Pi ■Р2 ■Рз • Р4• Р5= 0,2 0,2 ■0,9 • 0,9 -.0,125=0,004,

 

откуда

 

 

0,5249

-0,004 =

0,011;

 

 

 

Р(1) = -

 

 

 

вычисляется значение

0,2000

 

 

 

 

 

С (1)=С (0)+С ,

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С (0) = V

С}°) = 3 + 1 +

2 +

1 + 15 =

23; С/ = С, =

3,

 

откуда С(1) =23 + 3=26.

шагу

№ 2 [максимальная среди

оставшихся улО)

4.

Выбирается у по

для кф1 или у»(2)];

 

 

 

 

величина Р;,(1) или

Р,(2),

по табл. 8.1 отыскивается соответствующая

если шаг j\s 2 имеет у ,(2);

 

 

 

 

 

 

— вычисляется значение

 

 

 

 

 

 

Р (2) =

РП 1) Р(1)

ИЛИ Р( 2) =

Р / ( 2 )

Р ( 1 ) , если шаг № 2

имеет

 

 

РИО)

 

 

Р / ( 1 )

 

 

 

 

вычисляется значение

 

 

 

 

 

 

 

С (2)= С (1)+ С Л[или С(2)=С(1)+С,-.

 

 

если шаг Jx’s 2 имеет

i (2)]-

производятся на последующих шагах.

По ре­

Аналогичные вычисления

зультатам расчета для рассматриваемого примера составлена табл. 8.3, ко­

торая позволяет решить обе поставленные задачи оптимального резервиро­

вания.

1-я задача

(минимум стоимости), то процесс прекращает­

Если решается

ся на таком шаге

М, когда

выполняется условие Р(М—1)<Ртр<Р(А1,)-

Пусть требуемая надежность системы рассматриваемой в примере Ртр>0,47.

Из табл. 8.3 определяем, что оптимальная система

следующая: «^=4; пн —

=3; ш3=1;

т.|=1;

т5= 3 и стоимость ее — 69 единиц.

Таблица 8.3

 

 

 

 

 

 

м

mi

то

т3

 

т5

Р(Л4)

С (/И)

0

0

0

0

0

0

0,004

23

1

0

1

0

0

0

0,011

26

2

0

1

0

0

1

0,026

31

3

1

1

0

0

1

0,048

34

4

1

2

0

0

1

0,072

37

5

1

2

0

0

2

0,115

43

6

2

2

0

1

2

0,156

45

7

2

2

0

1

2

0,171

46

8

2

2

1

1

2

0,187

47

9

3

2

 

1

2

0,226

50

10

3

2

1

1

2

0,300

55

11

3

3

1

1

3

0,356

58

12

4

3

1

1

3

0,400

61

13

4

3

1

1

3

0,472

69

14

5

3

1

1

4

0,518

72

350

Если решается 2-я задача (максимум надежности при ограничении по стоимости), то процесс необходимо закончить на таком шаге М, когда вы­ полняется условие

Стр<С (М +1).

Если, например, стоимость системы не должна превышать СТр = 60 еди­ ниц, то процесс надо! закончить на шаге № II. Характеристики оптимальной системы будут: Р = 0,356, inl =3, m2= 3, т 3= 1 , /7Ц = 1, тъ = 3.

8. 2. РЕЗЕРВИРОВАНИЕ МОЩНЫХ ДВИГАТЕЛ ЬНЫХ УСТАНОВОК

8. 2. 1. Методы резервирования

Существуют различные сложные системы, которые состоят из большого количества идентичных как по устройству, так и по вы­ полняемым функциям подсистем. К подобным системам можно отнести двигательные установки мощных ракет-носителей, напри­ мер, «Сатурн-V» [52].

Такая двигательная установка состоит из п одинаковых бло­ ков, причем блоком является автономный двигатель, имеющий свою систему подачи и регулирования и работающий независи­ мо от других двигателей.

Тяга двигательной установки ракеты-носителя, особенно пер­ вой ступени, в зависимости от решаемых задач может находить­ ся в пределах 10 000—80 000 кН и даже более. Так, тяга первой ступени ракеты-носителя «Сатурн-V» составляет примерно 34000 кН. Поскольку, двигательную установку, состоящую из од­ ного блока, с такой тягой по причинам экономического и техни­ ческого характера создать крайне трудно, то тяга одного блока в настоящее время ограничена по величине. Например, двигатель F-1 имеет тягу 6800 кН.

Рассмотрим двигательную установку с тягой /?д.у, скомпоно­ ванную из п отдельных блоков с меньшими тягами Ra, где

И =zRr.y/Ro.

Если все блоки (двигатели) работают независимо и облада­ ют одинаковой надежностью, то вероятность исправной работы такой двигательной установки определяется зависимостью

Рд.у= П рб,-

(8.20)

;=1

 

где Pq — вероятность исправной работы отдельного двигателя.

Нетрудно убедиться в том, что даже при высокой надежности отдельных блоков и при большом их количестве надежность дви­ гательной установки будет чрезвычайно низкая.

Пример. Пусть необходимо создать двигательную установку носителя с

тягой /?д.у= 100000 кН из двигателей с тягой У?а=5000 кН, обладающих очень высокой надежностью (Рб=0,990).

351

В этом случае количество блоков определится из соотношения

п = — — = 'У*

== 20 и с о г л а с н о ф о р м у л е ( 8 . 2 0 ) Р (.у = 0 , 9 9 0 20 = 0 , 8 1 5 .

Рй оООО

Обеспечить требуемую надежность многоблочной двигатель­ ной установки позволяет метод общего резервирования с дробной кратностью, с постоянно включенным или с «холодным» резер­ вом.

При постоянно включенном резерве все блоки, как основные так и резервные, работают с момента старта ракеты. При хо­ лодном резервировании резервные блоки находятся в нерабочем состоянии до тех пор, пока не откажут основные. Когда отказы­ вает один из основных блоков, он выключается и в работу вклю­ чается один из резервных. Холодное резервирование обладает недостатками:

увеличивается пассивная масса носителя за счет массы ре­ зервных блоков, которые в полете могут не использоваться;

ухудшаются динамические характеристики ракеты из-за резкого изменения перегрузок в моменты выключения отказав­ ших блоков п включения резервных;

усложняется система выключения и запуска двигателей в условиях полета.

Так как из-за указанных недостатков применение холодного резервирования для повышения надежности мощных ракет-носи­ телей следует считать нецелесообразным, то в дальнейшем рас­ сматривается только лишь общее резервирование с постоянно включенным резервом и дробной кратностью.

Общее резервирование с постоянно включенным резервом мо­ жет быть двух видов: с постоянным и с переменным коэффици­ ентами форсирования тяги

Кф= RoI-Rg,

где Яг, — тяга на номинальном режиме работы. Рассмотрим два случая.

1. Резервирование с К ф = 1=const

Все как основные, так и резервные двигатели работают на номинальном режиме, и тяга их в процессе выключения отказав­ ших блоков не изменяется. При отказе от одного до т включи­ тельно блоков в остальных блоках тяга не изменяется, а суммар­ ная тяга двигательной установки будет уменьшаться по ступен чатому закону:

/?д.у(/) = ( л - / ) Я б,

( 8 . 2 1 )

где /= 0 , 1, 2,..., т — количество отказавших блоков. Недостатком данного метода резервирования является зна­

чительное уменьшение тяги двигательной установки при отказе

352

блоков ii, как следствие, снижение коэффициента тяговооруженпости ракеты.

Очевидно, что в этом случае для обеспечения выполнения за­ дачи, если даже уменьшение коэффициента тяговооруженности находится в допустимых пределах, необходимо увеличивать вре­ мя работы двигательной установки.

Разновидностью метода резервирования при /Сф= 1 является использование естественного резерва блоков в установленные моменты времени [68].

В процессе полета из-за израсходования компонентов топли­ ва и уменьшения массы носителя увеличивается коэффициент

тяговооруженности, в результате че­

 

 

го в некоторый момент времени об­

 

 

разуется естественный резерв тяги,

 

 

который равен тяге

одного

блока.

 

 

■Следовательно, с этого момента вре­

 

 

мени двигательная установка будет

 

 

располагать одним резервным

бло­

 

 

ком.

коэффициента

 

 

Текущее значение

 

 

тяговооруженности определяется за­

 

 

висимостью

 

 

 

 

n R 6

 

( 8. 22)

 

 

!ЧТ):

 

 

 

Оо — аь'п

 

Рис. 8.3. Схема

образования

где G6— расход топлива через один

резерва тяговооруженности

блок;

 

 

масса ракеты.

 

 

Qo — стартовая

 

 

Соответственно начальное значение коэффициента тяговоору

женности

 

 

 

 

 

Р(0) = - ^ С .

 

(8.23)

 

 

Qo

 

 

Из-за увеличения

второго

члена

знаменателя в

уравнении

{8.22) (тяга блока считается постоянной) с некоторого момента Т] наступает условие Ajlii= jli(тГ]) —• цо, позволяющее выключить один блок. При выключении блока коэффициент тяговооружен­ ности уменьшится, а затем начинает расти пока не наступит сле­ дующее условие Др,2=р.(т2) — ро, позволяющее выключить вто­ рой двигатель и т. д. Процесс образования резервных двигате­ лей в полете показан на рис. 8. 3.

Количество резервных блоков, которые могут образоваться в процессе полета, приблизительно можно определить следующим образом.

В предположении, что допустимое значение коэффициента тяговооруженности при отключении последовательно от 1 до т двигателей должно быть не меньше начального р0, можно запи­ сать условие

353

!А= !Ао + Н'ш.

 

 

(8.24)

где р. = ------------------------ значение коэффициента

тяговоору-

<?о —06(« — "О v

создаваемое

резервными

женности,

блоками;

,

,

 

(п т) Re.

 

р.= —---------—5------ значение коэффициента тяговоору-

Qo— G o ( n — m) г

при

отключении т бло­

женности

ков.

Ц в уравнение (8. 24) и

Подставив выражения для цо, Цш>

разрешая его относительно т, получим

 

 

 

т = -

 

 

(8.25)

*Д-У

 

 

 

Г-о1'

Для осуществления такого резервирования необходимо знать моменты времени т*, в которые образуется естественный резерв блоков.

Приближенно указанные моменты времени можно определить следующим образом.

После отказа к — 1 двигателя в момент п, масса израсходо­ ванного топлива определяется по зависимости

 

Q/—1= [« — (^Т- Л] ^йХк^Г^й ^

т/>

(8- 26)

где

 

 

 

 

й = 0

 

 

 

 

 

 

 

 

 

т,. = т(/ = 0, 1,

2 , . . . , k — \)\

 

т0 =

0;

т; > 0

— момент

отказа

/-го двигателя.

 

 

 

 

 

 

Полагая, что в момент отказа

к — 1 двигателя выполняется

условие ц(т) = цо, после подстановки

уравнения (8.26)

в зави­

симость (8. 25), получим

 

 

 

 

 

 

 

Л)—

(/г — к) R 6

л—1

 

 

 

 

 

 

 

Qo — [и — (* + 1)] ОбЪ — ° б ^

Ъ

 

 

 

 

 

 

*=о

 

Разрешив

последнее уравнение

относительно т/„ получим

 

 

 

 

/г-1

 

 

(8. 27)

 

xk= ak ~ b k V t;,

 

где

 

 

 

/=i

 

 

 

 

 

l

 

.

 

 

 

bk

 

 

 

 

 

п — (А+ 1)

 

 

 

 

 

 

 

ак= ьк

Оо

п — /г

 

 

 

 

 

Н-0

 

 

 

 

 

06

 

 

 

 

354

Уравнение (8. 27) позволяет определить моменты времени Т/;, начиная с которых образуются резервные блоки.

2.Резервирование при Кф>1-

Вмомент запуска все двигатели имеют номинальную тягу. При отказе одного двигателя в остальных тяга увеличивается на

величину, обеспечивающую ту же номинальную тягу двигатель­ ной установки. При отказе второго двигателя до т включитель­ но каждый раз увеличивается тяга работающих и тяга двига­ тельной установки не изменяется.

Величина коэффициента форсирования тяги блоков опреде­

ляется из соотношения

 

К ф= ^ = ^ -

(/-О , 1, 2 , . . . , т).

R 6 n —

i

Величина УСф ограничивается прочностью конструкции и за­ пасами устойчивости рабочего процесса. Для конкретного двига­ теля Лф величина заданная.

При заданном значении Кф можно определить допустимую кратность резервирования

тКф 1

а= — —---------.

п Кф

Всоставных двигательных установках отказавшие двигатели воздействуют на окружающие и могут вызвать отказ последних.

Взависимости от вида отказа воздействие отказавшего двигате­ ля на соседние может проявляться по-разному. Основными воз­ действиями являются: механическое, тепловое, вибрационное, ре­ жимное и комбинированное. Отказы взрывного характера сопро­ вождаются разрушением конструкции и разлетом осколков н элементов, которые, ударяя в окружающие двигатели, могут их разрушить или вызвать недопустимые изменения режима ра­ боты.

При нарушении герметичности полостей, заполненных горя­ чим газом (газогенераторы, камеры двигателя, газоводы и др.), возникает истечение струи газа, обладающей значительной тепло­ вой и кинетической энергией. Воздействие таких струй на окру­ жающие двигатели вызывает прогары элементов конструкции.

При нарушении герметичности жидкостных магистралей воз­ никает истечение компонентов топлива в окружающую среду, что может вызвать пожар в двигательном отсеке и другие по­ следствия.

Таким образом, резервирование двигателями само по себе не может обеспечить надежное выполнение задачи при их отка­ зах. Для обеспечения работоспособности и надежности резерви­ рованной двигательной установки необходимо применять систе­

355

мы аварийной защиты, которые определяют состояния двигате­ лей и при возникновении у них аварийных ситуаций вы­ ключают их, тем самым защищая соседние двигатели от воз­ действия отказавших.

Поэтому в последующем рассматривается теория резервиро­ вания двигательных установок с учетом взаимодействия систе­ мы аварийной защиты с двигателями.

8.3.ВЕРОЯТНОСТНАЯ МОДЕЛЬ РЕЗЕРВИРОВАННОЙ СИСТЕМЫ

8.3. 1. Вероятности состояний системы

Совокупность блока с системой аварийной защиты будем на­ зывать о б о б щ е н н о й с и с т е м о й (ОС), а совокупность обоб­ щенных систем — о б ъ е к т о м з а ш и т ы (03). В соответствии с приведенным определением вся резервированная двигательная установка является объектом защиты, а составляющие ее двп-

Рис. 8.4. Модель функционирования обобщенной системы

356

гателп совместно с системами аварийной защиты — обобщенны­ ми системами.

Пусть 03 состоит из п равноправных и одинаковых по своим характеристикам ОС, из которых т резервных. Объект защиты выполняет возложенные на него функции только тогда, когда ра­ ботоспособны п т обобщенных систем, а отказавшие до т си­ стем включительно выключены без разрушения и воздействия их на окружающие системы. Оставшиеся п т обобщенных си­ стем работают на постоянном режиме или переключаются на новый.

Состояние объекта защиты в момент времени т определяется состояниями обобщенных систем. Процесс эволюции состояний объекта защиты во времени представлен на рис. 8. 4 и описы­ вается функциями £;(т);=о,1,2 ш, которые характеризуют фазо­ вое состояние г-й обобщенной системы в момент т. Под фазовым состоянием понимается вероятность исправного состояния неко­

торого числа обобщенных систем, так что

все

!о(т) ■— вероятность

нулевого состояния, когда исправны

п обобщенных систем;

первого фазового состояния, когда

ис­

£i (t) — вероятность

правны п — 1 обобщенных систем, а неисправная обобщенная выключена системой аварийной защиты;

£,„(т)— верятность т-го фазового состояния, когда работа­ ют а т обобщенных систем, а т неисправных отключены си­ стемой аварийной защиты (САЗ). В момент т объект защиты' выполняет заданные функции, если он находится в одном из т фазовых состояний.

Так как перечисленные состояния являются несовместимы­ ми, то вероятность появления ровно т состояний за интервал времени от 0 до ть т. е. вероятность исправного функциониро­ вания объекта защиты, определится зависимостью

т

р ™ о о = 2 *'(*)•

(8- 28)

1=0

 

Рассмотрим подробнее процесс эволюции состояний объекта защиты (см. рис. 8. 4) в течение промежутка времени Дт, к на­ чалу которого он находился в одном из г'-х фазовых состоянии.

За промежуток времени Дт в зависимости от состояния обоб­ щенных систем объект защиты может находиться в г-ом состоя­ нии, перейти в г + 1 состояние или перейти в состояние отказа.

Вероятность того, что объект защиты к моменту времени т +Дт будет находиться в г-ом фазовом состоянии, в котором он. находился к моменту т, определяется зависимостью

£/,-(-г+ лт) = [1 — <7г(т, т +

дт)] [1—г,-(т, т +

дт)]

£,• (т), (8.29)

где <7i (т, т + Дт) — условная

вероятность того,

что

объект защи­

ты за время Дт перейдет из г'-го

состояния в-

357

аварийное при условии, что в момент т он на­ ходился в исправном /-ом состоянии;

г,;(т, т + Дт) — условная вероятность того, что за время Дт произойдет ложное срабатывание системы аварийной защиты и исправная обобщенная система будет выключена.

В /+1 фазовое состояние объект защиты переходит в тех слу­ чаях, когда произойдет ложное срабатывание одной из систем аварийной защиты, или когда система аварийной защиты обна­ ружит возникновение аварийной ситуации в одной из обобщен­ ных систем и выключит двигатель.

Вероятность перехода объекта защиты нз состояния / в со­ стояние /+1 определится зависимостью

S /(/+ u (T - f Дт) = [гДт, т + дт) + ?дт, т -f дт)(1 — <7 „)]Е£(т ), (8.30)

где qH— условная вероятность того, что система аварийной за­ щиты не обнаруживает наступающего аварийного состояния.

Аналогично можно описать эволюцию всех 0, 1,2, ... , т фа­ зовых состояний объекта защиты.

Из фазового состояния т объект защиты переходит в состоя­ ние аварии тогда, когда происходит ложное срабатывание САЗ или когда объект находится в аварийном состоянии.

В этом случае объект становится ненадежным и не выпол­ нит поставленную задачу.

Введем обозначения:

’Ы И

т - f дт) =

[1 —<7,(т, т -f дт)][1 —г,(г, t - f дт)];

j ^

'Wo'+iA'C

т - f дт) =

[г, (т, т -f дт) +<7 ,дт, т - f дт)(1 —?„)].

J

С учетом введенных обозначений вероятности перехода объекта

защиты в разные состояния

(8. 29) — (8. 30) перепишутся в виде

(т-f ДТ)=Ф,7 (т, т - f дт)£,(т);

(8.32)

b(i + l ) (T-f ДТ) —

(т, т - f дт) £;(Т).

 

Вероятность /-го состояния объекта защиты определяется за­ висимостью

S/(т д~ Дт)—■£(/—и /(т ~f AT)_1Sii (т-f дт),

(8.33)

где

 

. S(/—I)/(т-f Дт) = Ф(/_1)/(т, т -f дт)5,-_1(т);

 

Ф(/-1)/(т, т -f дт) = гг_1(т, т -f дт) <7 j- _ 1 ( t , т -f дт) (1

<7Н).

Таким образом, на основании зависимости (8.33) можно так записать систему уравнений, выражающих эволюцию обобщен­ ной системы за промежутки времени Дт:

358

?0 ( Г + Д Т ) = ф„о(Т| t + A T ) S j ( T ) ;

 

е1(Т + ДТ)=

фи1(г, т + дт)*э(т:)+ <|>ц (т,

t - f AT)6i(t);

£m( t + ДТ)=

Ф(„г_ 1)т(Т, T-f A ^ - l W

+ 'W 'K Т+ДТ.) \m(t).

В системе (8.34), произведя предельный переход при Дт=0, получим систему дифференциальных уравнений, подобных урав­ нениям А. Н. Колмогорова [24], которые определяют вероятность исправной работы объекта защиты в момент времени т:

(г)

Фоо (т)£о (^);

 

 

 

dx

 

 

 

 

 

 

 

 

(т)

—Фо1 i1") ^1 ("О “1Фи(Т) “Л

)l

|

dx

 

 

 

 

(8.35)

 

 

 

 

 

I

dim(Т)

 

 

 

 

J

dx

 

 

 

 

где

dijiy)

_ii (t + At)

_

 

d x

д-с-,0

Ат

 

 

 

| o(t) — условная вероятность пребывания объекта защиты в момент т в нулевом фазовом состоянии, т. е. когда исправны все п обобщенных систем, при условии, что при Дт=0 система находилась в этом же состоянии;

|,-(т) — в общем случае условная вероятность пребывания объекта защиты в момент т в состоянии i при усло­ вии, что при Дт—0 система находилась в состоя­ нии i 1;

Ф/i (т)= Пт

<7;(т,

т4-^т)+г/ (т, t-J-Дт)—q-t { т, т-ЬЛт)г/(г, тЧ-Д?)"

 

Д^-*-0

 

 

 

Дтз

 

со

 

 

 

 

 

 

 

со

Ф(/—1)/ =

lim

г-1- 1

(1?, т + ДТ) + <?/-[ (Т, т 4- Дт) (1 — д н)

об

 

Дт-^0

 

 

 

Дт

 

 

являются функциями

перехода объекта

защиты из

состояния

г—1 в состояние i

и

характеризуют

интенсивность

перехода,

которая

зависит от

надежности двигателя и точности работы

системы аварийной защиты.

защиты, показанная на

Таким образом,

эволюция объекта

рис. 8. 4,

состоящего из п обобщенных систем, из которых т ре­

зервных,

описывается системой дифференциальных уравнений с

переменными коэффициентами и является случайным процессом. Этот случайный процесс характеризуется следующими свой­

ствами:

359

Соседние файлы в папке книги из ГПНТБ