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

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

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

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

Добавлен: 23.06.2024

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

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

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

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

 

 

 

 

 

 

}

(ЗЛО)

где

 

 

 

 

£ ■ = /-* ;

A'ir( l , p ) = A jr( l - t p ) , /г =

11,

12,

21;

^ 2 2 (S, Р) =

lp/vіо (/ — 5)] + [Р/У20 (/ —È)1 - М 8а (/— S. Р).

Действительно, вектор-функция (Ф;1 (х,

р),

j =

1, 2),

дающая решение задачи [(3.7), (3.8)], пропорциональна вектор-функции (Ф'д (I — х, р), / = 1, 2) — решению за­ дачи [(3.10), (3.9)1 с коэффициентом пропорциональности

Фи (I,

р), поэтому должно выполняться

соотношение

 

Фц (I,

р) Ф 'п {I, р) = 1

 

(3.11)

и, следовательно, при

данном значении

р должно

быть

Ф 'и (/,

р) Ф 0. При выполнении этого

условия (ФЛ (х,

р), / =

1, 2) должно иметь вид

 

 

 

Фі\{х,р) = Ф'ц{1— х.р)ІФ[х(1,р),

/ = 1, 2 .

(3.12)

Обратно, вектор-функция (3.12) удовлетворяет уравнению (3.7) и граничным условиям (3.8); решение единственно в силу единственности решения задачи [(3.7), (3.8)]. Ус­ ловие Ф 'п (/, р) Ф 0 далее будем считать выполненным.

Таким

образом, искомые

передаточные

функции

от

к

üj (х,

t) (j — 1, 2)

выражаются через

решение

задачи

[(ЗЛО),

(3.9)] следующим

образом:

 

 

 

 

 

 

 

 

И77-і (х, р) =

 

 

 

 

=

[Ф/'і (I — X, p)l®'11(l>Р)]ехр( —ртДх)),

/ = 1 , 2 .

(3.13)

(/

Воздействие Ф2г.. Передаточные функции о т ^ г к i (х, t)

=

1, 2) определяются аналогично (3.13)

по формуле

 

 

 

 

й ^ (* , Р) =

 

 

 

 

 

 

=

[Ф/а {ж, р)/Ф ; 2У, Р)\ ехр [—рт2(х)],

/ = 1 , 2 ,

(3.14)

где

' , 2 (х, р), / =

1,

2] — решение

системы дифферен­

циальных уравнений

 

 

 

 

 

 

 

 

 

(^Ф1а/ 5 х ) - | - ( х ,

р) Ф1г —Л12 (х,

р) Ф24,

1

/g

 

 

( — d<t>;2/dx) + Л22(х, р )ф ; 2= Ап (х,

р)Ф1’Ъ

]

 

77


G граничными условиями

 

 

 

 

 

 

Ф12 (0. Р) = 0; Ф ;2 ( 0 , Р ) - 1 .

 

(3.16)

В

выражениях

(3.14) и

(3.15)

А"п (х, р)

=

(р/и10 (х)) +

+

(р/ѵso (х)) +

Л і

(х, р);

т2 (х) =

\d\lvw {I).

Здесь долж-

но

быть Ф' 22 (/,

р)ф 0.

 

X

р)/Ф' 22 (/, р),

Функция (Ф'^2 (х,

/ =

1, 2) является обобщенной

производной

регулярного

оригинала и теряет это свойство при умножении на ехр ра, где а > 0 .

Воздействие ѵ. Система уравнений

относительно пере­

даточных функций от V к

(х, t) (/ =

1, 2) имеет вид

dWlnldx+W in [(р/ц10 ( х ) Ь М п (X, р)] =

~(Х, р) ф Ais (х, р) 1Е2п,

dW2n/ d x + W 2n [(р/Оао (х)) + і422 (х, Р)} =

=2 (х, р) + Л21(х, р) ІКІГ1

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

 

Г 1п (0 ,р ) = 0 ; W2n(l,p) = 0.

(3.18)

[Можно показать, что вектор-функция (Ц7;-п (х, р), / =

1,

2)

является

обобщенной производной регулярного оригинала

и теряет

это свойство при умножении на ехр ра, где а >

0 .

Это означает, в соответствии с физическими представле­

ниями,

отсутствие

чистого

запаздывания реакции (öy- (х,

t),

/ =

1,

2) на воздействие

по ѵ.] Решение краевой задачи

[(3.17), (3.18)] удобно свести к решению задачи

Коши с ну­

левыми

условиями при X =

0 для той же системы уравне­

ний.

Обозначая решение

последней

задачи

[Ф7П (х, р),

/ =

1,2], получаем для решения задачи (3.18) соотношение

 

 

w Jn (X, р) = ф]п(х, р)— [Ф/2 (х, р)/ф ; 2 (/, р)] X

 

 

 

 

X ехр [т2 (/)—т2 (х)]Ф2п (/, р),

/ — 1,

2,

(3.19)

где [Ф'у2 (х, р), / =

1,2] — решение задачи [(3.15),

(3.16)],

о котором выше предполагалось Ф' 22 (/,

р) Ф 0. Единствен­

ность решения задачи [(3.17), (3.18)] вытекает из единствен­ ности решений задач Коши [(3.15), (3.16)] и (3.17) с нулевы­ ми граничными условиями.


Воздействия уг (г — 1, 2). Система уравнений относи­ тельно передаточных функций от утк -öy (х, р) (/ = 1, 2) имеет вид

öWiGrldx + WlGr[(plvw (х)) А Ап (х, р)] =

=В, АХ, р) + А,, (х, р) Wogг\

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

U^ior(0, р) = 0 , W 2or(l, р)= = 0 г = 1, 2 .

(3.21)

Аналогично [(3.17), (3.18)] решение задачи [(3.20), (3.21)]

представляется в виде

WjQr (х, р) = Фусг (х, р) — [Фу-2 (х, р)/Ф'2(/, р)] X

X ехр [т2(/) —т2 (х)] Ф2ог(/, р), / , / - = 1 , 2 ,

(3.22)

где (Ф/Gr (X, р), /

= 1, 2) — решение системы (3.20) с ну­

левыми при х =

0

граничными условиями.

 

Воздействие # T

является функцией х; при каждом фик­

сированном

виде функции От (х, t) изображение

реакции

( (х, t), /

= 1,

2)

находят интегрированием системы урав­

нений

 

 

 

 

d-öj/dx +

fti [(plvw (х))ф Ап (х, р)]

 

 

= Zx(х, р) # т + А12 (х , р) ^ 2;

(3.23)

— ö ö 2/ ö x - f # 2

[(р/^20 {х )) + А 22{х , р)]

 

=Л22(х, р)Ъх

сграничными условиями

■öi (0 , р) = 0 ,

(/, р) = 0 .

(3.24)

Решение краевой задачи [(3.23), (3.24)] производится так же, как решение задачи [(3.17), (3.18)].

Численное решение систем уравнений для передаточных функций. Таким образом, отыскание изображений реакций (Oj- (х, t) } — 1, 2) на входные воздействия сводится к ре­ шению нескольких задач Коши для однородной и неодно­ родной систем уравнений с комплексными и переменными

79


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

dw1l d x ~ a 11w1 + a12iW'i + b1\ 1

^

dWildx^a^Wi + a^Wz + bz \

 

эти формулы имеют вид

 

 

wie = f+ (А) wu +

/- (А)

Wu + an Wn+Ьг

 

+

/о(А) (а12 bl-- a22 h)\

(3.26)

W2e= /+ (A) W2i +

/_ (A) ^ 22~ an да2г+ an wu +b2

 

+ /o (A) (a2i[bi — an b2),

где wji (/ = 1, 2) — значения координат wj (x) в начале участка; wje (/' = 1, 2) — значения их в конце участка;

А—длина участка; /+ (А) = (exp q+А + exp q_А)/2;

 

....

(е^-ьА—еч- Л)

при

b=j=О,

(А ) =

2 /

Ь

 

 

 

 

 

 

 

Ае?+ Л

при

b — О,

1

 

 

а 11 + а 22

/- (А)

/+ (А) при

а 11 а 2 2 ~ а 12 а 21

т

 

 

^11 а22

®І2 а2І>

 

 

 

 

 

 

 

j

— А

Р ? +

А _

ПрИ СІц #22 :—^12 ^21

 

а11+ а22

а11+ а22

V

 

 

 

 

 

И

1^22 Ч=

 

 

80