Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Параметрическое проектирование расчёт и исследование траекторий дви..pdf
Скачиваний:
54
Добавлен:
15.11.2022
Размер:
15.12 Mб
Скачать

В) Расчет центровочных и моментных характеристику-й «сухой» субракеты,

у = 1, 2,3, •••>/*•

- расстояние от носика ракеты до плоскости крепления переходного отсека на (/- 1)-й

ступени

 

 

 

/'

- /

_ /вн

_ /

V - 1

*пгу

*сопл/-1

*заддну-1’

-определение коэффициентов ATi и для в общем случае усеченного конуса (переходный отсек) из табл. 2.4:

К, = /(£>у_, / Dj}, i = 1,2;

- положение ЦМ переходного отсека

*перех отс/ ~ ^/-1

^псрсх отс/^1»

- определение коэффициентов АГ3, К4, АГ5для эллиптических днищ из табл. 2.5:

K ,= ft{a lb \i =3,4,5;

- расстояние от носика ракеты до плоскости переходного отсека на (/-1)~й ступени

h^у-1 ^иерех отс/ ’

-положение ЦМ переднего днища ДУ

* д н перу - * f / “ 0 -

* 3 )^ДН перу ’

- положение ЦМ цилиндрической части ДУ

*цнлу = h

^цилу I

-расстояние от носика ракеты до плоскости крепления переходного отсека на (/+1)-й ступени

V , *цид/ ^ ^цил/ 1

- положение ЦМ заднего днища ДУ

^дн заду ” /у+1 + ( 1 " * 3Удн зад/ >

—определение коэффициентов ^ 1,^2 для усеченного конуса (сопловой блок), из табл. 2.4:

K = >;ldnj)i = h 2;

- положение ЦМ соплового блока ДУ

соплу — / у + 1 ^

^днзвду

(^сопл/

^сопл вн/ )>

- положение ЦМ «сухой» частиу'-й субракеты

■^cyV “ (^п е р о х о т с / ^ перех отс/

* д н перу ^ д н перу

^ ц и д /^ ц и д /

* д н з а д у ^ д н зад/

+ ^ с о г ц //Я с о п д / + ^ п г у ) ^ с

Ух /^

 

 

-моменты инерции соответствующих частейу-й субракеты относительно главных центральных плоскостей zoCfro и^оС^о:

C U = 0,25-mnepcK^[(D y /2)2+ (£,_, /2)г] ;

 

=Q>$'mmi(Djl2)1;

 

j£ b j= K sm ^(D j/2 ? ',

 

= 0,25- m ^ K ^ /2)* + (<*„ !2?\,

 

*C °

 

 

 

 

тУо°*о

~

Y m

12

 

 

и ЛН перу

 

i V 4 А7*дн пер;4 дм пер J

 

• С ° = 0 ,0 8 3 3 .т ^ С ;

7>оО*о

_ к m

/2

.

-

Y 'm 12

47 дн заду

* М ш дн задУ^дн заду ’

47 сопл/

* М т со1Ц /*с о п д /’

- моменты инерции тех же частейу-й субракеты относительно координатных осей хохо иуоУо (см. табл. 2.3)

] * 0*0

 

_

9 .

 

Т Ч 0хъ

 

47перех о т#

 

^

47перех о тс /»

 

Г*о*о

__ 9 .

Г ^ Д о .

 

 

и дилеру

 

^

 

47 дн перу»

 

 

7*о*о

9 .

 

7*о°*о .

 

 

,7цид/ ” ^

^цид/ >

 

 

7*о*б

_

9 ,

г

^

о .

 

 

47днзаду

^

^днзад/»

 

 

7*о*о

_

9 .

 

7го0хо .

 

 

47 с о т у

 

^

47 сопл/ >

 

ТУоУо

=

7 ^ 0

 

 

.

7*о<Э*о

47перех о т#

47перех отс/

~

47перех отс/ ’

7>о>'о

-

ТУо°Ч

.

 

/*о О *о .

47дн перу

 

47дн иер/ ~

47дн пер/ >

 

ТУоУо

_

7>оОхо

,

 

7*ол о •

47цн#

 

47шц/

 

' 47ЦИД/

5

7У0У0 _

тУоОь

I

 

7*оО*о .

47дн заду

 

47дн зад/

* 47днзад/ >

 

7>о>о

_

7 ^o0zo

1 Т Ч 0хй .

^СОПД/

 

17СОПЛ/

47сопл/

»

- моменты инерции «сухой» частиу - й субракеты относительно ее ЦМ

7*0*6

__

7*0*0

I

7*0*6

I

7*0*6

_L

7*0*0

7*6*0

_1_ 7*0*6 .

47 суjy

•”

47 перех о т#

^

47 дн пер у

^

и пял/

^

47 дн заду ^

47 сопл,'

"'л гу ’

=

t^ n q p e x отс/ + W nepex отс/ ( X cyxJ " Х перех о тс /)

1 + [ ^ д н

перу +

171m

пер/ ( * с у х /

* д н пер/ ) 1 +

“^[*^цшу

'«цид/ (*сух/ “

*цил/ )

1

[Л н заду

'«дн зад/ (*суд/ ~~ *дн зад/)

] "^*

 

" Н ^ с о п д /

+ '«соп лу (-^ с у х / “

* с о п д /)

1

[* A iry

^ п г у

( * с у х у ~

* п г у )

] 5

 

 

С) Расчет центровочных и моментных характеристикj - й субракеты,у = 1 , 2 , 3 , / ? ,

сучетом изменения на АУТ системы «топливо/ + ТЗП/», i = 1,2,

-значение относительной массы ракеты

М-теч// = 1 - ° Л С ^ - 1);

(если р ^ , < р^. = / юоу, то положить р ^ , = р*,);

- унесенная масса системы «топливо + ТЗП»

т х унесу/ = ™0у 0 — M'teiy/)»

- текущее время работы ДУ в точках / = 1, 2, 3, ...

'т с к/I = '« т у н е с // /'« с у м м у »

Г Д е '«сум м у = ' ” т у + < п / ;

- оставшаяся на борту ракеты масса системы «топливо + ТЗП»

'« о сту/

'«тО у

'« т у п е е // >

- значение массы ракеты в точках i = 1, 2, 3,...

'« у / = '« 0 /М 'те ц // >

—значение массы системы «топливо + ТЗП» в переднем днище, цилиндрической части и заднем днище ДУ соответственно

WTW = " гтплОУ( 1 - ^ / т у );

тттщ! ~ ^тиилО/О- t-mji llj)t

-положение ЦМ системы «топливо + ТЗП» в переднем днище, в цилиндрической части и в заднем днище ДУ соответственно

* т пд//

= /’#-0,375 /

V

д н п е р у J

 

л т цил//

^цид/ »

* т з д /,

= *)♦. +0,375./дн заду »

- положение общего ЦМ системы «топливо + ТЗП» ДУ

X y t ~~ (*т ДН перу/^T дн перу/ *т цилу/^т цид// *т дн э&д//'«т дн заду/ V '«т осту/ »

- положение ЦМj - й субракеты в точках / = 1, 2, 3,...

Х1, = (*«,»<„, +XVim

;

-моменты инерции системы «топливо + ТЗП» переднего днища, цилиндрической части и заднего днища ДУ относительно главных центральных плоскостей г^Охц и yoOzo соответственно

,= 0,2-т гютг^ / 2 ) 2;

=0,25'«тшц/([(Ду/2)2 +(^>нутру(/2)2] ;

(2.7)

/пнзм//= ®>2 -mTOT34f((Dy/ 2) ;

(в выражении (2.7) внутренний диаметр топливного заряда ^внутр// определяется из зависимости dIKmJ, = D] -2 - и /т ,

=0> 05938 • rnTmnbfj^mai!f]>

=0,0833

“УлиМ«у/ = 0,059387ЯТщ,,^/дн >

-моменты инерции системы «топливо + ТЗП» переднего днища, цилиндрической части и заднего днища ДУ относительно координатных осей *о*ои УоУо (см. табл. 2.3)

 

/Vo

_ 9 . Г*о°*о

/V»

_ 9 , Г*О0*О•

 

7*о*Ь

- 9 . 7*00*0

 

 

^ т д а п е р у /

 

А ^ т д н п е р у /5

 

 

 

** ^ т ц н д / / ’

 

17 т дн м д у/

^

47 т дн м д // ’

 

7 W 0

7^о0*о

.

7*о0*о •

ТУоУо

 

7Уо0*о .

7*о0*о

ТУоУо

_

7^0*0

I

7*о0*о

47тд н п е р у/

^ т д н п е р у /

~

*7 тдаперу/ 9

^ т ц ш у /

 

17тц ил // “

47т ц ш у / *

^тднзаду/

17т да зад// ~ ^тдаэад//9

- моменты инерции системы «топливо + ТЗП»у - й ДУ относительно ее ЦМ

 

 

 

7 * 0*0

7 * 0*0

 

I

7 * 0*0

I

7 * 0*0

 

 

 

 

 

 

 

 

47т/7

^тдаперу/ ~

^тцид/Г ~

^тдазад// >

 

 

 

 

*А//

[ ^ т да перу/

^тднперу/ (*т//

* т да перу/) ]

Urцпд//

цш у/ ( * т у / * т

цид/7)

1

+ 1Л дн заду/

да зад// (*ту/ ~

дн зад//)

] >

 

 

 

 

 

 

 

 

 

- моменты инерцииу-й субракеты относительно ее ЦМ

 

 

 

 

 

 

 

 

 

7 *0*0

_

7 *0*0 .

7 *0*0 .

>

 

 

 

 

 

 

 

 

 

47

у/

47 суд/ ^

^ т у

/

 

 

 

 

 

 

= [> С + » ^ ( * w

-* „ )2]+ К Г ° + * ,« * (* * -* „ )2] .

 

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

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

При проведении комплекса расчетов в обеспечение процесса определения

характеристик движения ракет на АУТ принято ступени и субракеты р-ступенчатой

составной ракеты нумеровать в восходящем порядке, начиная от основания ракеты к ее

вершине по схеме рис. 2.7 (ср. рис. 1.1), что исторически сложилось для расчетов летных характеристик ракет; при этом для 1-й субракеты (то есть ракеты на старте) рабочей ступенью будет 1-я ступень, а ее «полезной нагрузкой» - 2-я субракета вместе с переходным отсеком между 1-й и 2-й ступенями и т.д.

Расчет характеристик АУТ может производится, например, методом Рунге - Кутта 4-го порядка [49, с. 125-130] с постоянным заранее заданным шагом по времени.

Для проведения расчетов привлекаются:

аэродинамические характеристики, которые задаются в табличном виде;

моменты инерции и относительные координаты центра масс ракеты, задаваемые

втабличном виде для варианта, получаемого в процессе параметрического проектирования (подраздел 2.9).

Впроцессе проведения расчетов при необходимости определяются динамические коэффициенты, возмущающие силы и возмущающие моменты, действующие на ракету при ее движении на АУТ (подраздел 2.10.4).

Все модули блока расчета характеристик движения ракеты на АУТ и на ПУТ используются:

при определении характеристик рассеивания точек падения ПГ методом

статистических испытаний (подраздел 2.10.6);

при расчете проектно-конструкторских частных производных спроектированной

ракеты методом конечных приращений (подраздел 2.10.7);

для уточнения потерь конечной скорости ракеты в процессе параметрического проектирования (разделы 3-5).

Полная дальность полета ПГ определяется по характеристикам конца АУТ с

привлечением зависимостей эллиптической теории [4] или численным интегрированием

(раздел 7).

Рис. 2.7. Схема составной />-ступенчатой ракеты, принятая для расчета характеристик движения на АУТ

Численное моделирование движения ракеты на АУТ производится в предположении, что поле силы тяжести Земли - центральное (т.е. вектор силы тяжести направлен к центру Земли), а СУ - идеальна. Кроме того, предполагается плоская постановка задачи (т.е. движение ракеты осуществляется в плоскости, проходящей через центр Земли, точку старта и точку окончания АУТ). Вращение Земли не учитывается. Не учитывается и влияние составляющих от нормальных сил аэродинамических рулей.

Расчет траектории движения ракеты на АУТ производится с использованием системы дифференциальных уравнений, описывающих движение ЦМ ракеты в земной системе координат (рис. 2.8).

J

(2 .8)

У

где

bxj = (Xj cos0 + Yj sin0- Rxj coscp + RyJsinep)/m;

byJ= (Xj sin 0 - Yj cos 0 - RXJsin cp + RyJcos (p)/ m\

0’ =arCSin(yj/v;); W'=W^n(h)\

Wj = W'KWj\

a = p -9;

q = phVj /2; M = VJ/ah‘,

 

ph = p'h(\+Kphy,

c x =c;(\+KCly

 

 

rh=T^(\ + Kn y,

с у =с'у(\+к Суу

 

p*= pl(i+ ^J;

c d = cd(\ + Kcdy,

щ = щ 1(\+кту,

CT=C’T(\ + KCjy,

J t i n j ~ У в п у О + ^ J a , )>

R j ~

n ) 8 o ™ T j ~ P h ^ o J >

5биш= arcsin(-gS^C,

« a r c s i n ( - ^ - ^ ) ;

 

1 Cj

1 Cj A y

RXJ= Rj cos(56ara);

RyJ= R, sin(56мл);

0 = * /$Ш* ± ^ ~ gC°s6 ,

 

 

 

mjvj

wo, = « w 0 +

+ mTy(l + KmlJ) + тщ;

x(.t0J) =x0;

y(t0J) = y(l-

x(70j) = x0; y(t0J) = y0;

и при условиях в точках разделения ступеней на АУТ:

Ж ,) = Ж -1); Ж у )= Я О ;

Ж у)=*('Г); Ж у )= Ж 7"1);

7 —2, 3

, р ь

до момента, определяемого уравнением выключения последней (верхней) ДУ на АУТ. Примечание. Индекс «*» относит рассматриваемую характеристику к разряду номинальных. Коэффициенты Kwj, КриУ ..., Kmij являются случайными числами, назначаемыми датчиком случайных чисел в процессе статиспытаний. Для номинального расчета эти коэффициенты равны нулю за исключением Kwj, который равен единице.

Одновременно с интегрированием системы (2.8) осуществляется численное интегрирование системы дифференциальных уравнений:

А^у=-gsin0;

Д^сх/ ^-Xj/nij'i

У

(2.9)

А ^ ; = - p hS aj c o s a / m / i

при начальных условиях на каждом АУТ:

AKg(r0y) = AFCx(/0y) = A ^ ( r 0y) = 0;

j - 1,2,3, ...,р.

Система (2.9) определяет потери скорости на преодоление силы притяжения Земли, на аэродинамическое сопротивление атмосферы, на изменение тяги с высотой - раздельно по каждому АУТ р-ступенчатой ракеты. Далее рассчитывается суммарные потери конечной скорости:

A V * = i A K j - >

A V S * = t A W '

A V Cx = h A V C xj>

 

j~l

 

AK£ = AKg + AV& + AVSa.

Эти потери могут быть использованы, в частности, для анализа и корректировки некоторых результатов расчета (например, коэффициентов) в процессе параметрического проектирования.

Пояснения к рис. 2.8:

т. О - точка старта ракеты

При отсечке ДУ в текущий момент:

-т. В - верхняя точка траектории полета ракеты;

-т. С - точка падения ракеты на поверхность Земли;

-дуга Ок земного шара - дальность точки конца АУТ - /к;

-дуга кС земного шара - дальность ПУТ - /св.

Плоские системы координат:

-хО\у - скоростная (поточная) система координат

-х\0\у\ - связанная система координат

-хсОус- стартовая система координат

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

R sin 5б1Л(/-х т) + Уcos а(хд-х,)+ Х sinа(хд - хт),

откуда

 

 

8«.пя = arcsin г хд- х т

Y cos а Хд-хт X sina^ = arcsin

— —— (Уcos а + X sin а)

V

1-х,

(l-x T)R

или при условии малости углов атаки а (так называемые летные углы атаки)

5б1Л= arcsin

— —— (Ycosa + Xsma) » arcsin ' (Сд ~ С Ж

 

(l-x T)R

(1-ст)Л

Рис. 2.9. Схема сил, действующих на ракету

Расчет программы тангажа сучетом ограничений на углы атаки

Движение любой ракеты типа УБР или PH от точки старта до точки начала самостоятельного движения ПГ (боевой части УБР или ИКО для PH) осуществляется по программе угла тангажа закона изменения угла наклона оси ракеты к стартовому горизонту. Обычно в условиях подобной постановки задачи закон изменения угла тангажа задается в зависимости от времени. Таким образом, уравнения движения интегрируются при заданном законе изменения угла тангажа, что однозначно определяет совершенно конкретную реализацию траектории движения ракеты. В результате конечные значения фазовых характеристик траектории определяют шш дальность полета ПГ (для УБР) или характеристики орбиты ИКО (для PH).

Как правило, на программу тангажа <р(/) накладываются ограничения со стороны условий прочности, устойчивости, удобства приборной реализации и т.п. В первом приближении таким требованиям удовлетворяют классические ограничения [4, с. 279]:

-старт с определенной продолжительностью движения при постоянном (вертикальном или наклонном) угле тангажа;

- непрерывность ср(0, ф(0> ф(0 и ограниченность ф(Г);

-ограниченность нормальных (т.е.поперечных) нагрузок;

-нулевые углы атаки при скоростях, близких к звуковым (при прохождении звукового барьера);

-специальные условия, например околонулевые углы атаки при прохождении зоны максимальных скоростных напоров;

-обеспечение для УБР любой дальности из определенного диапазона с минимальным количеством программ ф(г) и т.п.

Перечисленным ограничениям в условиях проектных расчетов (в отличие от, например, расчетов попадающих траекторий для обеспечения контрольных или боевых пусков ракет) удовлетворяет программа ф(Г) вида:

-для одноступенчатых ракет форхмулы и графики (рис. 2.10, 2.11) изменения программы

для составных/т-ступенчатых ракет формулы и графики (рис. 2.12,2.13) изменения программы тангажа фпр(0 и углов атаки а(Г) имеют вид:

Рис. 2.12. График изменения программы тангажа српрСО Для многоступенчатых ракет

О(0 'i

 

 

 

 

 

 

 

 

 

 

 

ао(0

а(0=0

ai(0

_

а2(0 ^

 

 

--------- ►!*

,

 

.

.

*(0=о i

 

1

 

1

1

 

1

1

 

1

1

Л

 

t

 

'

1

 

1

1

 

 

1

 

—1

J

 

 

 

 

 

 

 

 

~

т

'

 

/

 

 

 

 

 

 

\

Г \

&тахО

 

\^шах1

\ ^шах2

 

 

\ j

 

 

 

 

 

1 (Участок наведении есть)'

L

а=фк-е

J

;

 

 

!

 

 

|1

^

!

1

£*i

A

a

!

! \ г \

и-тахл

,

| (Участка наведения нет)

J

L М О J Рис. 2.13. График изменения углов атаки а(Г) для одноступенчатых ракет

При этом для 2, 3,..., р-й ступеней составных ракет