Файл: Динамика и управление ядерным ракетным двигателем [Текст] 1974. - 253 с.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 23.06.2024

Просмотров: 221

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

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

В случае несжимаемой жидкости в качестве аналога па­ дения напряжения естественно взять разность входного и выходного давлений рвХрвыХ. Аналог уравнений Кирх­ гофа в этом случае может быть записан (учитывая также наличие насосов в описываемом контуре) в следующем ви­ де:

 

 

N

п

 

 

 

(

/

 

где H j

— напор

і -го насоса; Ьі — константа, относящаяся

к і-й

ветви

контура

и определяемая

выражением

bi — hh/(2diP2S 2i), где

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

гидравлическо­

го сопротивления для /-го канала; dt — диаметр, /, — дли­ на; St — площадь сечения; р — плотность. Следует отме­ тить, что так как в наиболее сложных и разветвленных частях гидравлической схемы любого ракетного двигателя газ чаще всего можно считать практически совершенным и, кроме того, в этих трактах насосы, как правило, отсутству­ ют, то при описании газовых трактов ракетных двигателей чаще всего используется первый аналог уравнения Кирх­ гофа второго рода, а именно уравнение (3.90'.

В том случае, когда интерес представляет характер изменения давления в какой-либо точке тракта, возникает необходимость присоединения соответствующего допол­ нительного уравнения. Причем желательно, чтобы в этом уравнении интересующее давление рх было выражено по возможности только через одно из входных давлений тракта и не включало неинтересующих давлений. Такие уравнения могут быть написаны по аналогии с уравнениями Кирхгофа 2-го рода для электрической разомкнутой цепи («падение напряжения равно алгебраической сумме э. д. с. плюс алгебраическая сумма падений напряжения»). Дейст­ вительно, если просуммировать уравнения (3.89) для всех

ПО


элементов любой разомкнутой цепи газового тракта, то

получим уравнение

 

м

 

(р^)2- ( р Д 2 = 2 С ,0 ? П Р<,

(3.91)

І—1

 

являющееся указанным аналогом в случае

отсутствия

в описываемой цепи насосов. Относительно использования систем уравнений, состоящих из уравнений типа (3.88), (3.90) , (3.91), можно сказать следующее. Решения этих систем, дающие правильную информацию о величинах расходов, давлений и о направлении течения, в случае произвольного задания направлений течения не могут быть получены, если хотя бы одно из направлений задано неверно. Для получения правильных решений этих систем все направления течения должны быть заранее известны. Это связано с тем, что уравнение (3.89) и уравнения (3.90), (3.91) , получающиеся суммированием уравнений типа (3.89), верны только в случае, когда направления течения извест­ ны. В свою очередь, это обусловлено тем, что в уравнениях (3.90), (3.91) знак нелинейных членов сгС2;Г*ср t не зависит от знака G*. Именно наличие этих членов не позволяет аналитически определить истинное направление течения, если оно заранее неизвестно; в этом случае получаются решения, не дающие никакой информации (расход выража­ ется, например, комплексной величиной).

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

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

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

Этому требованию отвечает следующая форма

записи

уравнения импульса:

 

(P bx) 2- ( p L ,x) 2 = cG | G | T * p,

(3.92)

111


где |G| — модуль расхода, а индексы «вх», «вых» обозна­ чают вход в элемент и выход из него при принятом направле­ нии течения газа. Этому уравнению соответствуют уравне­ ния

1 сі G* I G/l 7’cp i = 0;

(3.93)

;= l

 

M

 

(pSx)a-(pJ?)s = ^ ct О, I Gt I r cpi,

(3.94)

i= 1

 

аналогичные уравнениям (3.90), (3.91), но отличающиеся тем, что в них знаки членов ctGi\G\lTtp г зависят от знака расхода G. При решении систем, состоящих из уравнений типа (3.88), (3.93), (3.94), направления течения газа могут быть заданы произвольно. Если после решения системы ка­ кие-либо значения Gt получились отрицательными, необ­ ходимо изменить направление движения газа в соответст­ вующем элементе на противоположное принятому. Вопрос о методах решения указанных систем уравнений требует специального рассмотрения.

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

цией уравнений (3.88), (3.90),

(3.91):

 

 

 

 

 

і= 1G o г І

~

0 ;

 

 

N

 

 

N

 

 

N

 

 

2£г ji +

i= 1

^cp i + zL kt рг = 0;

i= 1

 

 

 

 

i—\

 

M

 

 

 

 

 

 

 

 

•ГС* = (Pbxo!Px o)2 ^bx [ l/(p* o)2]

>J

 

 

 

 

 

 

 

 

( = 1

 

 

M

 

 

 

 

 

M

— [1/2 (p*o)2J

2

kt ^

p, - [ l/2 (p io )a] 2 kiVi-,

 

 

i= 1

 

 

 

 

i= 1

где kt = (pBXо/)2

(Рвых о/)2;

ii — (ßi — Gi0)/Gi0;

A

*

.

IT*

т*

ѵ

т

*

Пер г — (Л cp I

1

ср 0(7' 1

ср 0і>

 

 

n*x = (p*x- p Z 0)lp*x0-,

 

Л вх — (Рвх

 

рвх о)/Рвхо!

 

Р г = (с*

£ jo ) /c jO—

относительное изменение гидравлического сопротивления г'-го элемента тракта.

112


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

В заключение подчеркнем, что во всех приведенных выше

уравнениях

суммы

2 необходимо понимать

как алгебраи-

ческие, т. е.

знаки

І

нужно брать

членов, входящих под 2 ,

 

 

і

 

в соответствии с направлением течения газа.

СПИСОК ЛИТЕРАТУРЫ

1.Жирицкий Г. С. Авиационные газовые турбины. М., Оборон-

гиз, 1950.

2.Широков М. Ф. Физические основы газодинамики. М., Физ-

матгиз, 1958, с. 76.

3.Лойцянский Л. Г. Механика жидкости и газа. М., Физматгиз, 1959, с. 159, 479, 631, 634.

4.Абрамович Г. Н. Прикладная газовая динамика. М., «Наука»,

1969, с. 35, 172.

5.Жуков В. П., Креер Р. И. «Об одном эффективном методе иссле­ дования динамики сложных газовых трактов». — В сб. «Ав­ томатическое регулирование двигателей летательных аппара­ тов». Вып. 13, М., ЦИАМ, 1972, с. 209.

Г л а в а 4

ИССЛЕДОВАНИЕ ДИНАМИКИ

ИВОПРОСЫ УПРАВЛЕНИЯ ЯРД

СТВЕРДЫМ ДЕЛЯЩИМСЯ ВЕЩЕСТВОМ

§11. Исследование ЯРД как объекта управления вблизи номинального режима

Цифровой и аналоговый способы исследования дина мики. С точки зрения теории управления ядерный ракет­ ный двигатель представляет собой достаточно сложную систему с большим количеством динамических звеньев и с перекрестными связями между ними. Поэтому важным этапом разработки системы управления для ЯРД должно стать составление такой математической модели, которая была бы, с одной стороны, достаточно точной, чтобы отра­ зить принцип действия объекта и основные физические за­ кономерности в его звеньях, а с другой — достаточно простой, чтобы облегчить исследователю понимание проте-. кания процессов в объекте и выбор методов управления им. Особенно важно составление простой математической модели для применения современных аналитических методов проектирования системы управления [как для маршевого режима (см. § 12), так и для пускового (см. § 14)].

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

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

ризует динамику какого-либо звена и имеет вид

 

 

J

 

У

Wj(m)zj,

(4.1)

 

/= 1

 

1І4


где

у — выход

звена;

Zj

( / = 1-=-«/) — входы

звена;

Wj — передаточная функция

звена по отношению к г/,

о» >

0. Решение

системы

осуществляется для

каждого

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

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

вычислять

значения частотных характеристик

каждого

из звеньев,

либо заранее рассчитывать частотные

характе­

ристики звеньев с распределенными параметрами и вводить эти результаты в машину в виде таблицы.

Полученные частотные характеристики объекта следует аппроксимировать характеристиками возможно более простых типовых звеньев; результатом этой работы будет математическая модель объекта, к которой можно приме­ нять различные (аналитические или аналоговые) методы синтеза системы управления; один из аналитических мето­ дов синтеза изложен в § 12.

Аналоговый способ исследования динамики объекта опи­ рается на описание динамики каждого из звеньев обыкно­ венным дифференциальным уравнением (здесь речь идет лишь о применении стандартного аналогового оборудова­ ния). Поэтому для объекта, содержащего звенья с распреде­ ленными параметрами, необходимо предварительно заменить эти звенья какими-либо достаточно простыми звеньями с со­ средоточенными параметрами. Это удобно сделать, пред­ варительно рассчитав частотную характеристику каждого из этих звеньев; звенья, рассмотренные в § 3, 6, обычно уда­ ется с приемлемой точностью аппроксимировать типовыми звеньями не выше второго порядка.

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

115