Материал: Погорелов В.И.-Прочность и устойчивость тонкостенных конструкций

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

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

Pmax = 2,67 Mnd при y0 = d ; Pmax = 4 Mnd при y0 = 3d / 4 .

Рис. 34. Эпюры напряжений в болтах

Расчет внутренних усилий в шпангоуте проводится для каждой из сил, создаваемых единичным болтом, а затем результаты суммируются в каждом сечении. На рис. 35 показаны направления внутреннего изгибающего момента М, крутящего Н и перерезывающей силы Q, возникающих в сечении шпангоута под действием сосредоточенной силы Pi . Там же приводятся их эпюры в зависимости от угла ϕ. Расчет внутренних усилий ведётся по

следующим зависимостям [26]:

 

M = kM Pi r ,

где r = d / 2 – радиус окружности болтов;

kM

=

1

[ϕsin ϕ+

1

cosϕ −1]

;

2π

2

 

 

 

 

 

Рис. 35. Внутренние усилия и моменты в шпангоуте

69

Q = kQ Pi , где kQ = 21π[ϕ− 2sin ϕ] ;

H = kH Pi r , где kH = 21π[ 23 sin ϕ−ϕ−ϕcosϕ] .

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

6. НАГРЕВ ГОЛОВНОГО ОТСЕКА

6.1. Тепловой апогей

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

q =

1

с

 

ρvс

 

(T T

 

) ,

(6.1)

 

2

 

f

 

p

r

w

 

 

где cf – коэффициент трения; cp – теплоемкость газа при постоянном давлении; Tr,, Tw температура восстановления и температура поверхности cтенки соответственно, видно, что функция q должна иметь экстремум, так как плотность увеличивается, а скорость уменьшается. Найдем высоту, на которой тепловой поток принимает максимальное значение. Выразим входящие в формулу (6.1) переменные величины через высоту H. Аппроксимируем плотность от высоты зависимостью, приведенной ранее в подразд. 1.1, полагая

 

 

 

 

 

 

 

 

ρ = ρ e−βH ,

где β = 0,15 103 (1/м),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

а

изменение

скорости по

формуле

(1.5). Далее,

так

как

T

= T

 

 

+ r

γ −1

M2

 

 

 

T T = T T + rT

 

γ −1

M2

 

1

 

 

 

 

, то

 

 

 

 

 

r

 

 

 

 

2

 

 

 

 

 

r

w

w

2

 

 

rT

γ −1

M2

,

 

 

где

разность

(T

T

)

обычно

значительно

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

w

 

 

 

 

 

 

 

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

70

M2 =

v 2

=

v 2

=

 

 

v 2 γ

=

v 2

 

 

2

γRT

γс

p

(γ −1)T

(γ −1)с

p

T

 

a

 

 

 

 

 

 

и тогда

 

 

rT (γ −1)v2

 

 

 

T

T

=

=

rv2

 

 

.

2(γ −1)сpT

 

r

W

 

 

2сp

 

 

 

 

Подставим теперь полученные соотношения в (6.1) и полу-

чим

q =

1

сf ρvc p

rv 2

=

1

сf rρv 3

2

2сp

4

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

q =

1

f

ρ v3 exp(−βH 3

Бe−βH

) .

 

 

4

 

 

3 a

 

 

 

 

β

 

 

 

 

 

 

 

 

Находим производную по H и приравниваем ее нулю, считая что r,сf , и Б от высоты не зависят:

dHdq = 14 f ρ3va3 exp(−βH 3βБe−βH ) (−β + 3βБe−βH β) = 0 ,

откуда 1+3

Б

e

−βH

= 0 и высота, на которой возникает экстре-

β

 

 

 

 

 

 

 

 

мум теплового потока, равна:

 

 

 

 

 

 

 

Hm =

1 ln

.

(6.2)

 

 

 

 

 

β

β

 

На этой высоте возникает максимальный тепловой поток к поверхности головного отсека, входящего в атмосферу. Его величина зависит от массы и аэродинамического сопротивления отсека.

Конечно, высота, на которой достигается экстремум q, определена приближенно, так как использовалась простейшая формула для теплового потока к пластинке в несжимаемом газе при Pr=1, однако если формулой (6.2) воспользоваться только для определения высоты, а при расчете тепловых потоков учесть все необходимые факторы, то можно получить более точные результаты.

71

6.2. Конструктивные пути снижения нагрева

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

счет изменения его характеристик. Если обозначить через Q

искомое количество тепла, то производная dQdt представляет

собой тепловой поток, и поэтому можно записать следующее равенство:

q = dQ dt =1 4(c f rρv 3 ) .

(6.3)

Перейдем в полученном обыкновенном дифференциальном уравнении от независимой переменной t к H. Если считать, что при входе в атмосферу головной отсек движется по прямой, на-

клоненной к плоскости горизонта под углом наклона ϑа вектора скорости в конце активного участка траектории, то тогда

dHdt = −v sin ϑa ,

а уравнение (6.3) можно записать так:

dQ = 1 c f rρv 2 dH 4 sin ϑa

или после подстановки выражений для плотности и скорости в зависимости от высоты (см. поразд. 1.1):

 

 

 

 

 

 

 

 

 

 

 

 

1 c f rρЗva2 −βH 2

Б

 

 

 

 

 

 

 

 

 

 

dQ

 

 

 

 

 

 

 

 

e−βH

 

 

 

 

 

 

 

 

 

 

 

 

 

β

.

 

 

 

 

 

 

 

 

= −

 

 

 

 

 

 

 

 

 

 

e

 

 

 

 

 

 

 

 

 

 

 

 

 

dH

 

4

sin ϑa

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

f

rρ v

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Обозначая

А=

 

 

3

а

 

 

 

 

и

вводя независимую переменную

4sin ϑa

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z = e−βH ;

dz = −βe−βH dH = −βzdH

 

получим

 

 

 

 

 

dQ

 

 

A

 

2

Бz

 

 

 

 

QΣ

 

 

A

1 2

Бz

 

 

 

 

 

 

β

 

 

 

 

 

 

 

 

β

 

 

 

 

=

 

e

 

 

 

 

 

 

 

или

 

dQ

=

 

 

 

 

e

 

 

dz ,

 

 

dH

β

 

 

 

 

 

 

 

 

 

β

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

0

 

 

 

 

 

откуда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Бz

 

 

 

 

 

 

 

 

 

 

 

2

Б

 

 

 

 

 

 

 

 

 

 

 

 

A

β

 

 

 

 

A

 

 

 

 

 

β

 

 

 

 

 

 

QΣ = −

 

 

 

 

 

 

 

e

 

 

 

=

 

 

 

 

1

e

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

72

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

атмосферу. При 2(Б β) < 0,4

экспоненту можно разложить в ряд,

ограничившись

только

линейными

слагаемыми,

т.е.

e2(Б β) 12(Б β)

и QΣ ( А β) , т.е. количество тепла не зави-

сит от конструктивных характеристик головного отсека. Если

2Б/β ≥ 2,4, то e2(Бβ) 0 и можно уменьшить количество поступающего тепла за счет увеличения конструктивного параметра

 

сx ρ3Sm

 

 

Б=

1

. Этого можно достичь за счет увеличения cx

и

 

 

2m0 sinϑa

1

 

 

 

 

площади миделя Sm или снижения массы головного отсека m0. Большой коэффициент силы лобового сопротивления cx1 возни-

кает у плохообтекаемых тел, перед которыми образуется мощный скачок уплотнения и резко возрастает волновая составляющая сопротивления. При этом скорость ГО у поверхности Земли мала. Можно получить умеренные значения cx1 и сверхзвуковые ско-

рости движения ГО у Земли, если выполнить его в виде комбинации конусов с притуплением в носовой части. Другой возможный вариант – создание ГО с регулируемым cx1 , что достигается за

счет запрограммированного сбрасывания аэродинамических стабилизаторов или носовых обтекателей.

6.3. Прогрев теплозащиты до температуры уноса

Процесс нагрева ТЗП можно разделить на три этапа [8]. На первом этапе происходит движение тепловой волны до того момента, когда температура ТЗП на поверхности, обращенной к воздуху (рис. 36), становится

равной

температуре

уноса

 

покрытия. Далее

начинается

 

второй этап нагрева, на

 

котором

 

происходит

 

интенсивный унос

покрытия.

 

Для ГО унос покрытия начина-

 

ется на

высотах

40-50 км и

Рис. 36. Распределение температуры

заканчивается на 10-15 км, где

в ТЗП на первом этапе

тепловой

поток

к

ТЗП73