
книги из ГПНТБ / Донов А.Е. Динамика полета летательных аппаратов
.pdfп |
Р cos (а + ?дв) |
_ |
Q |
(252) |
G |
|
С |
||
|
|
|
||
или для малых значений а и удв |
|
|
|
|
|
(а + <P<J + |
О |
ао) |
(253) |
|
|
P - Q |
(254) |
|
G |
||
|
§ 2. ПРОГРАММИРОВАНИЕ АКТИВНОГО УЧАСТКА ИЛИ УЧАСТКА ВЫХОДА НА МАРШ
Для того чтобы заставить направление движения центра тя жести аппарата на активном участке изменяться должным обра зом, нужно при помощи управляющей аппаратуры реализовать какую-то определенную программу движения. Это можно сде
лать только при помощи центростремительного ускорения j „,
создаваемого силой Рп, так как только при наличии j n напра вление движения будет поворачиваться. Угловая скорость этого поворота (по отношению к земле) в данном случае определится
следующим |
образом: |
|
|
|
|
|
|
|
||
J n _- |
п |
~ |
cos I |
ё |
Р sin (а + yd„) |
Y |
(255) |
|||
v |
~~ |
Gv ® |
v |
Gv |
Gv |
|||||
Таким образом, |
мы видим, |
что |
состоит из трех составляющих, |
|||||||
обусловленных |
|
наличием: |
силы |
веса, |
тяги |
и аэродинамической |
||||
подъемной |
силы |
Y. Следовательно, какова |
бы ни была система |
|||||||
управления, |
получить требуемое |
изменение направления полета |
||||||||
можно только |
путем |
воздействия на |
эти |
составляющие. |
Осу |
|||||
ществить это воздействие можно многими |
способами. Рассмот |
|||||||||
рим наиболее характерные из них. |
|
|
|
|
||||||
|
Вращение аппарата вокруг центра тяжести |
|
|
|||||||
Изменение |
угловой |
скорости |
вращения |
аппарата |
вокруг его |
центра тяжести О приводит к изменению угла атаки, благодаря
чему происходит изменение |
нормальной составляющей тяги и |
||
подъемной силы Y, а |
следовательно, |
и изменение последних |
|
двух составляющих <v |
Так, |
например, |
если мы на фиг. 12, <5 |
будем увеличивать угловую скорость вращения аппарата, то это
вызовет увеличение угла атаки а, |
которое повлечет за |
собой |
увеличение нормальной составляющей тяги и подъемной силы Y. |
||
Следовательно, здесь увеличение |
вызовет увеличение |
и^. |
120
Данный способ воздействия на wv будет эффективен как в атмосфере, так и в пустоте. Только в атмосфере рассматриваемое воздействие осуществляется и через тягу и через подъемную силу, а в пустоте — только через тягу.
Для реализации исследуемого способа воздействия на необходимо иметь средства, позволяющие изменять угловую скорость вращения аппарата . В качестве таких средств обычно используются газовые или аэродинамические рули, или же те и другие вместе. Но можно для той же цели использовать и си стему быстровращающихся масс, помещенных внутри аппарата. Влияя тем или другим способом на угловую скорость вращения этих масс, можно изменять нужным образом и угловую ско рость шг вращения аппарата.
Изменение направления тяги двигателя
Поворачивая один или несколько двигателей, то есть изме няя угол <pde, можно непосредственно влиять на составляющую «>„, обусловленную тягой двигателя. Так как в общем случае сум марная тяга не проходит через центр тяжести аппарата, то из менение ее направления в общем случае вызывает также изме нение момента, создаваемого тягой относительно центра тяжести аппарата, а, значит, и изменение ш2. Таким образом, в данном случае на основное влияние на <ov через угол удв накладывается вторичное влияние через угол атаки.
Изменение величины тяги двигателей
Изменение величины тяги двигателя оказывает непосредст
венное воздействие на вторую составляющую для |
|
в формуле |
|||||
(255). Особенно сильно это воздействие будет |
в том случае, |
||||||
когда тяга составляет большой угол с направлением полета. |
|||||||
Однако имеется вторичное |
влияние изменения |
тяги на %. |
|||||
Это |
вторичное |
влияние |
возникает из-за того, |
что |
изменение |
||
тяги |
непосредственно связано с изменением скорости |
v и веса G |
|||||
(через |
расход |
горючего). |
Между тем названные |
сейчас пара |
|||
метры |
входят во все три |
члена |
выражения для |
|
Рассматри |
ваемое вторичное влияние сохраняется даже в том случае, когда отсутствует атмосфера, а тяга направлена по скорости движения, но имеется сила веса, направленная под углом к этой скорости. Влиять на % через величину тяги Р проще всего посредством регулирования расхода горючего.
Работа агрегатов, реализующих воздействие на % (рули, по
воротные двигатели, |
вращающиеся массы и т. |
п.), всегда связана |
|
с работой сложного |
комплекса |
специальных |
приборов и меха |
низмов. В автономной системе |
управления Еесь этот комплекс |
||
размещается в самом |
летательном аппарате. |
В неавтономной |
системе часть приборов и механизмов находится на земле. Если
121
система управления полностью автоматизирована, агрегаты, реа лизующие воздействие на u>„, работают таким образом, что в полете воспроизводится какой-нибудь заранее рассчитанный и „заложенный в аппаратуру системы управления параметр дви жения, который называется параметром программирования или параметром регулирования.
Параметрами программирования могут быть: угол, составлен ный осью ракеты с неподвижным направлением в простран стве; угол, составленный скоростью полета с этим неподвижным
направлением; |
проекция |
перегрузки на какое-нибудь направле |
||
ние, не |
совпадающее |
с |
направлением движения (например, |
|
проекция |
на нормаль к |
оси ракеты); угол атаки; отклонение |
||
руля; угол, составленный |
осью поворотного двигателя с осью |
|||
ракеты; высота |
полета; |
секундный расход горючего, если воз |
||
действие |
на сог, достигается |
изменением величины тяги, направ |
ленной под углом к скорости; угол между линией, соединяющей центр тяжести аппарата с неподвижной точкой земной поверх ности, и горизонталью, если аппарат движется вдоль радио луча, и т. д. Параметр программирования можно задать в виде зависимости от времени, скорости, высоты полета, количества израсходованного горючего и т. д.
Не всякий закон изменения параметра программирования, заложенный“ в аппаратуру системы управления, при попытке его реализации в полете приведет к удовлетворительным ре зультатам. Произойдет эго по следующим причинам:
1) программу нельзя реализовать из-за того, что значения аэродинамических параметров (например, коэффициента су и угла атаки) не могут быть обеспечены аэродинамическими фор мами летательного аппарата;
2) программу нельзя реализовать из-за того, что система управления не сможет обеспечить реализацию. Например, „за ложенной" программе могут соответствовать недопустимо боль шие потребные отклонения руля или же недопустимо большие скорости этого отклонения;
3) программу нельзя реализовать из-за того, что не будет обеспечена прочность летательного аппарата. Например, могут получиться чрезмерно большие потребные значения перегрузки или скоростного напора, которые приведут к разрушению аппа рата или порче его аппаратуры;
4) программу нельзя реализовать из-за того, что аэродинами ческий нагрев летательного аппарата будет чрезмерно велик;
5) программа не обеспечит вывода летательного аппарата в нужный режим полета;
6) летные характеристики, соответствующие „заложенной" программе, не являются оптимальными или же вообще не обес печивают выполнение задачи, возлагаемой на летательный аппа рат.
122
Таким образом, программа должна удовлетворять ряду усло вий, содержащихся в негативной форме в перечисленных сейчас пунктах.
В наиболее строгой постановке расчет программы обычно сводится к решению весьма сложной вариационной задачи. Ее сложность обусловливается наличием большого количества рассмотренных нами выше разнородных физических условий, которым должна удовлетворять функция, характеризующая программу. Вследствие этой сложности, не говоря уже о том, чго в строгой постановке (то есть с учетом всех возможных фи зических условий) само составление дифференциальных уравне ний экстремалей весьма затруднительно, численное интегриро вание получившейся системы уравнений, содержащей большое число неопределенных параметров, оказывается настолько гро моздким, что даже использование новейших достижений вычис лительной техники не приводит в данном случае к сколько-ни будь эффективным результатам. Поэтому основным эффектив ным способом расчета программы следует считать метод под бора более или менее подходящей программы. Этот метод состоит в том, что задаются различные варианты достаточно правдоподобных и совершенно определенных программ. Затем для каждого варианта рассчитывается движение, причем в про цессе расчета бракуются все варианты, не удовлетворяющие сформулированным выше условиям. После этого из всех рас смотренных вариантов выбирается наилучший.
Метод подбора можно существенно усовершенствовать, если придать ему целесообразную математическую форму. Так, на пример, искомую функцию, определяющую программу, можно разыскивать в виде заданной аналитической зависимости, содер жащей ряд неопределенных параметров. В данном случае реше ние задачи сводится к расчету движения для различных значе ний этих параметров с выбрасыванием в процессе расчета не пригодных вариантов и к выбору самого оптимального варианта.
При использовании рассматриваемого сейчас метода расчета программы аналитическую зависимость, содержащую неопре деленные параметры, целесообразно задавать не для самого параметра программирования, а для параметра О, непосредст венно характеризующего наклон скорости к горизонту, нужный
же для |
.закладки" в аппаратуру системы управления |
закон |
изменения параметра программирования рассчигывать»по |
пара |
|
метру 0. |
Эту аналитическую зависимость для 0 мы в дальней |
|
шем будем задавать как функцию времени t: |
|
|
|
0 = 6(t). |
(256) |
Чем больше будет параметров в рассматриваемой аналити ческой зависимости, тем лучше можно подобрать программу, но зато объем вычислительной работы будет больше. При бес
123
конечном числе параметров найденное таким образом решение совпадет с точным решением вариационной задачи.
Одним из самых основных параметров программы является время ее отработки аппаратурой, поэтому зависимость угла О
В В
от |
времени целесообразно разыскивать |
в виде |
|
0 = 6(т), |
(257) |
где |
— относительное время: |
|
X= Y • |
(258) |
Варьируя время отработки программы Т для одной и той же функции б (т), мы будем получать различные программы, отли чающиеся друг от друга только степенью растянутости прохож дения одних и тех же значений 0 во времени. На фиг. 88 изо бражен вид программы 0 (т) для участка АВ перехода от на клонного или вертикального участка движения с постоянным
124
значением угла б к участку горизонтального полета (фиг. 89, а, 6) и соответствующие графики зависимости 6от /
для различных значений Т (фиг. 90).
Зависимость б(х) удобно разыскивать в виде полинома
О(х) = о0+ 4-а2-2 + а„хл. |
(260) |
Удобство в данном случае состоит в том, что ряд коэффициен
тов этого полинома |
может |
быть |
определен |
сразу |
из условий |
|||||
„стыковки" |
движения |
на |
участке |
отработки |
данной программы |
|||||
с соседними |
участками |
движения. |
Рассмотрим |
на |
конкретном |
|||||
примере программы, |
изображенной |
на |
фиг. |
88, |
в чем состоят |
|||||
эти условия „стыковки". |
|
|
|
А |
траектории, |
показанной |
||||
Если при переходе через точку |
||||||||||
на фиг. 89, |
на летательный |
аппарат |
не |
действует |
ударная на |
грузка, мгновенно изменяющая скорость по величине и направ
лению, то угол |
б при переходе через эту |
точку не должен ис |
|
пытывать скачкообразного изменения. Поэтому для |
точки А |
||
разыскиваемой |
программы (257), то есть |
при т = 0, |
будет вы |
полняться условие |
|
|
|
|
0= 1 |
|
(261) |
Рассмотрим теперь производную т от функции б:
|
dG _db |
_d t _ |
■ |
|
(262) |
||
|
с/т |
|
dt |
с/х |
|
|
|
|
|
|
|
|
|||
Пользуясь формулой (189), |
имеем |
|
|
|
|||
f |
= n = - e |
l n |
r |
+ s l |
11* |
cos 6. |
(263) |
с/х |
|
V |
у |
V |
Vкост2 |
|
|
Если при прохождении центра тяжести аппарата через точку А нет перехода от одной ступени работы силовой установки к другой, который обычно сопровождается скачкообразным изме нением массы, тяги и ' аэродинамических сил, то перегрузка пу в точке А также не будет испытывать скачкообразного измене ния. Следовательно, на основании формулы (263) можно утвер
ждать, |
что не произойдет |
также |
и скачкообразного изменения |
|||
|
|
„ |
db „ |
на участке, предшествовавшем участку |
||
производной |
— . Но |
|||||
АВ, |
эта |
производная |
равна нулю, так как мы предположили, |
|||
что |
на |
этом |
участке |
б = |
const. |
Следовательно, при х = 0 для |
125
разыскиваемой программы (257), кроме условия (261), будет выполняться еще следующее условие:
# = |
(264) |
ах
Проводя аналогичные рассуждения и сделав те же самые предположения об отсутствии приложения ударной нагрузки и скачкообразного изменения перегрузки пу для перехода центра тяжести аппарата через точку В, мы получим для искомой про граммы (257) еще два условия: при т = 1
0 = |
0, |
(265) |
% = |
0- |
(266) |
Если взять полином (260) третьей степени, то условий (261), (264), (265) и (266) окажется достаточно для определения всех четырех коэффициентов а0, а,, а2, ал, так как эти условия при ведут к следующей линейной системе уравнений для их опре деления:
й0= 6о |
|
||
а, |
= 0 |
|
|
а 0+ <7, f а2+ а 3— О |
(267) |
||
Д} + 2и<>+ За3= 0 |
|
||
Решая эту систему, находим |
|
||
«о = °о |
|
|
|
ах= 0 |
(268) |
||
Д-2z:z:::r |
60Q |
||
|
|||
« з = 2 6 |
, , |
|
|
Следовательно, |
|
|
|
О(т) = 60— 30(Д2+ 20от3. |
(269) |
Если бы мы не воспользовались предположением об отсут ствии скачкообразного изменения перегрузки в точках А и В, то вместо условий (261), (264), (265) и (266) получились бы не сколько иные, более сложные. Но четыре коэффициента про граммного полинома по-прежнему определились бы из системы четырех линейных уравнений.
126
Выражение (269) не является единственным полиномом, удов летворяющим условиям (261), (264), (265) и (266). Чтобы в этом убедиться, рассмотрим такой полином Ф ( т ) :
Ф (Т) = 1’ (1 - т2) (Ь0 + blz + ... + b„ х » ) , |
(270) |
где Ьл, bu ...b„ — произвольные постоянные коэффициенты. Не |
||||
посредственная проверка показывает, что х - |
0, и при |
'=1 |
||
Ф = 0 | |
|
|
|
|
d± |
=0 • |
|
|
(271) |
d x |
I |
|
|
|
Поэтому, если мы составим сумму |
|
|
|
|
6= 60- 30/-’ + 20от3+ |
(1_ x f (b0 + Ь , х + ... + Ьат»), |
(272) |
то 0, определяемое из формулы (272), будет удовлетворять
условиям (261), |
(264), |
(265) и (266) |
при любых значениях коэф |
||
фициентов bQ, |
bl . . . b n. Подбирая |
соответствующим образом |
|||
эти коэффициенты, можно получить |
программу, |
которая будет |
|||
лучше удовлетворять |
сформулированным |
выше |
требованиям, |
||
чем программа (269). |
|
|
участке АВ фиг. 89 |
||
Исследуя для программы движения на |
условия „стыковки", мы рассмотрели только значения самой
функции 0 и ее |
первой |
производной. |
Однако |
анализ |
можно |
|||||||||||||
было бы продолжить |
дальше и определить производные более |
|||||||||||||||||
высокого порядка. При этом |
|
нужно |
было |
бы |
воспользоваться |
|||||||||||||
такими очевидными соображениями, что |
аппарат |
и управляю |
||||||||||||||||
щий орган |
(газовый |
руль, |
поворотный |
двигатель |
и т. |
п.) не |
||||||||||||
могут мгновенно изменить свое угловое |
положение |
или угло |
||||||||||||||||
вую скорость вращения. Точно |
так |
же |
не могут мгновенно из |
|||||||||||||||
мениться усилия, |
приводящие |
в движение управляющие органы. |
||||||||||||||||
В результате |
получились бы |
условия |
для значений |
производ |
||||||||||||||
ных более высокого порядка, |
чем первый, и этим |
условиям не |
||||||||||||||||
удовлетворяли бы ни программа (269), ни программа (272). |
||||||||||||||||||
Поэтому |
ни программа (269), |
ни программа (272) |
не являются |
|||||||||||||||
программами, |
удовлетворяющими |
всем |
условиям |
„стыковки", |
||||||||||||||
и если бы мы их |
„заложили" |
в |
управляющую аппаратуру, то |
|||||||||||||||
после |
прохождения |
точки |
|
А |
аппарат |
не выполнил |
бы точно |
|||||||||||
эти программы, |
а |
после |
прохождения |
точки |
В не |
пошел бы |
||||||||||||
строго |
по горизонтали. |
Получились |
бы отклонения от ожидае |
|||||||||||||||
мого движения, особенно сильные |
в районах |
точек |
Л и В. |
|||||||||||||||
Если у аппарата |
|
обеспечена |
достаточная |
устойчивость, |
эти от |
|||||||||||||
клонения должны быстро затухнуть. |
|
|
|
|
|
|
|
|||||||||||
Но если программы (269) и (272) и не являются наилучшими |
||||||||||||||||||
с точки |
зрения |
обеспечения |
|
точности их |
выполнения, |
то для |
||||||||||||
оценочных |
расчетов и |
анализа |
летных |
характеристик |
аппарата, |
127
а также выбора его конструктивных параметров эти программы вполне пригодны.
Анализируя вопросы программирования активного участка, мы до сих пор считали тягу двигателей и закон расходования го рючего заданными и рассматривали только вопросы програм мирования искривления траектории. Между тем при разыскании оптимальных условий выполнения активного участка в целом вопросы программирования тяги и расхода горючего являются также чрезвычайно важными.
В простейшем случае при анализе этих вопросов названные величины обычно считаются постоянными по времени для каж дой ступени работы ракетной системы. Варьируются только значения этих постоянных величин. Однако часто бывает по лезно считать их переменными в пределах работы каждой сту пени. При этом опять-таки в первом приближении каждая сту пень может быть разбита на участки, на которых тягу и расход горючего можно считать постоянными по времени.
§ 3. АНАЛИЗ ИЗМЕНЕНИЯ СКОРОСТИ ПОЛЕТА НА АКТИВНОМ УЧАСТКЕ ИЛИ УЧАСТКЕ ВЫХОДА НА МАРШ
Возьмем одну из ступеней работы силовой установки. Сум марную тягу на этой ступени будем считать постоянной:
Р = const. |
(273) |
Суммарный секундный расход горючего ch также примем по стоянным:
ch — const. |
(274) |
При сделанных предположениях относительно Р и ch будет постоянна также удельная тяга Руд:
Считая угол, составленный тягой с направлением движения, малым, мы сможем воспользоваться для анализа изменения скорости формулами (250) и (254). На основании этих формул находим
dv = |
|
dt ~ S'sin 6dt. |
(276) |
||
Принимая за начало отсчета |
времени |
начало работы |
ступени, |
||
получим, выполняя интегрирование, |
|
|
|||
v — v 0 = g |
Р |
dt —g |
l sin S dt. |
(277) |
|
D |
|||||
о |
|
о |
|
||
|
|
|
128
Применим |
к первому |
интегралу найденного выражения для |
||||||
v — v 0 известную теорему |
о среднем значении интеграла: |
|
|
|||||
|
v — v0= |
g 1— |
ср |
G dt — g |
l sin 6dt, |
|
|
(278) |
|
|
|
|
|
|
|
|
|
( |
Q \ |
|
|
|
отношения |
Q |
. |
Вос |
где I |
-p ) — некоторое среднее значение |
|
пользовавшись условиями (273) и (274), мы сможем взять пер вый интеграл правой части соотношения (278):
dt
|
|
|
G0 |
|
|
|
P . |
G0 — c,,t |
|
D i G |
n . |
|
(279) |
I n |
r |
“ |
Руд |
Руд |
• |
|
t'A |
'- 'П |
|
' - 'f l |
|
|
|
где
G
V4~ 7T (280)
G 0 ‘
Обозначим через T продолжительность работы ступени. Эта продолжительность всегда может быть вычислена, если даны удельная тяга Руд, значение р. в конце работы ступени ^ и
тяговооруженность в начале работы ступени ср0= v-/« :
|
7’= —— ^ (£о_= ^о |
|
= |
^ k ( i -jx,). |
(281) |
||
|
|
Р |
Ch |
|
То |
|
|
Вводя |
в рассмотрение |
безразмерное |
время |
т, второй интеграл |
|||
формулы (278) представим в следующем виде: |
|
||||||
|
Г |
т |
|
|
z |
|
|
|
J sin 0 dt = Т J |
sin б rf-c = — |
|
J |
sin б d~. |
(282) |
|
Теперь мы можем |
выражению |
для v — v0 придать следую |
|||||
щий вид: |
|
|
|
|
|
|
|
|
v — v0 — — gP.Уд |
ср |
In р — gT I sin б rfx = |
|
|||
|
|
|
J |
|
|
|
|
= |
- g p » \ i - ( |
# |
1п(Л— ^ ~ ( 1 - |
|1,) JsinerfT. |
(283) |
||
|
|
ср |
- |
|
|
|
|
9 А. Е. Донов |
129 |