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

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

 

 

 

 

1

−μ

−μ

 

= (1 + μ2 )(1 2μ).

 

 

 

 

 

 

 

 

=

 

−μ 1

−μ

 

 

 

 

 

 

−μ

−μ

1

 

 

 

 

 

 

 

 

Кроме того

 

 

 

x = εx E(1 −μ2 )− EαT (1 + μ)2 ,

 

поэтому

 

 

 

 

 

 

 

 

 

 

 

 

E(1 −μ)

 

αET

 

 

 

σx =

x =

 

 

εx

 

 

 

 

 

 

 

.

(3.8)

 

 

(1 +μ)(1 2μ)

(1 2μ)

Обозначим через

k =

 

 

 

E

 

 

 

– объемный модуль упруго-

3(1 2μ)

сти, а c2

=

E(1 −μ)

квадрат

скорости звука.

Тогда

(1 + μ)(1 2μ)ρ

 

σx = ρc2εx 3αkT , откуда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

εx = u

=

σx + 3αkT .

(3.9)

 

 

 

 

 

 

x

 

 

 

 

ρc2

 

 

 

 

Из уравнения равновесия после его дифференцирования по x получаем

2σ

x

= ρ

2

u

= ρ

2ε

x

 

 

 

 

 

 

t2

 

x2

 

 

x

 

t2

или, после подстановки (3.9),

2σx =

1

 

2σx +

α

2T .

(3.10)

 

 

x2

c2

t2

3 k

t2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для определения температуры в стенке необходимо решить уравнение теплопроводности (3.7) для температурного поля с объемными источниками тепла. Начальные и граничные условия имеют вид

t = 0 T (x,0)= 0 ;

T (x,t)→ 0 при x → ∞ ;

(3.11)

Tx (0,t)= 0 ,

где последнее условие показывает, что за время td распростране-

ния теплового импульса на глубину l можно пренебречь теплообменом на поверхности стенки.

49

Решение (3.7) с учетом (3.11) имеет вид

 

 

 

р

0

a

 

 

 

x

 

 

 

 

 

 

 

 

 

 

T (x,t)=

 

 

 

exp

 

 

[tH (t)

(t td )H (t td )],

 

 

λ

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где p0 = q l , или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

р0atd

 

 

 

x

 

 

t

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

T (x,t)=

 

 

 

 

exp

 

 

 

 

 

 

H (t)

 

 

 

1 H (t td ) .

λ

 

 

 

 

t

 

t

 

 

 

 

 

 

 

l

 

d

 

 

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При td 0 функция Хевисайда

H (t td ) H (t) , поэтому

T (x,t)= Tm exp(− x l)H (t), где температура поверхности

 

Tm = рλ0a td .

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

Q = р0ltd = qtd сохраняет фиксированное значение при td 0 .

Для решения волнового уравнения (3.10) поставим к нему следующие условия:

а) до воздействия импульса напряжения в конструкции равны нулю, т.е.: при t < 0 и x > 0 , σx (x,t)= 0 ;

б) на значительном удалении от поверхности напряжения стремятся к нулю: при x → ∞ и t > 0 , σ(x,t)→ 0 ;

в) скачкообразное изменение напряжений в момент начала воздействия теплового импульса:

 

t = 0 ;

σ(0,t)= 0 , σx (x,0 +)= −3αkT (x,0 +),

(3.12)

но так

как

du dx = 0 и

поэтому εx (x,0 +)= 0 ,

то

из (3.8)

∂σx (x,0 +)= 0 при x > 0 .

 

 

 

 

 

 

 

 

t

 

видно, что начальные напряжения при t=0 возрас-

Из (3.12)

тают скачком.

 

 

 

 

 

 

ξ = x l

 

η= ct l

 

Введя новые независимые переменные

и

в

(3.10):

 

 

 

 

 

 

 

 

 

 

 

∂σx =

∂σx ∂ξ =

1 ∂σx ;

2σx

=

1

2σx ;

2σx

= c2 2σx

;

 

x

∂ξ ∂x

l ∂ξ

x2

 

l 2 ∂ξ2

t2

l 2 ∂η2

 

50

 

 

 

 

 

 

 

 

1 2σ

 

 

 

 

 

 

1

 

2σ

 

 

 

3αk 2T

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

=

 

 

 

 

 

 

 

 

x

+

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

∂ξ

2

 

 

 

 

l

 

∂η

2

 

 

 

l

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂η

m

 

 

 

 

 

и обозначение для безразмерных напряжений

Θ =

 

σx

 

, полу-

 

3αkT

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

чим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

2

 

 

 

 

2

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

Θ

 

Θ

 

 

 

 

 

T

 

 

 

Θ

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

=

 

 

2

+

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

2 +

 

 

 

 

 

exp

 

H (t) .

∂ξ

 

 

 

 

2

 

 

 

 

 

 

∂η

 

∂η

2

l

 

 

∂η

 

 

∂η

 

 

Tm

 

 

 

 

 

 

 

 

 

 

 

 

Решение этого уравнения имеет вид

Θ(ξ,η)= −exp(− ξ)ch η при ξ > η; Θ(ξ,η)= exp(− η)sh ξ при ξ < η.

Вфиксированный момент времени, когда η постоянно, при

ξ> η напряжения сжимающие, равны нулю при ξ → ∞ и имеют

максимальное значение при ξ = η (рис. 20). Здесь напряжения имеют скачок, так как при ξ < η они растягивающие и при ξ = 0

стремятся к нулю. Найдем максимальное напряжение, полагая

ξ = η.

Рис. 20. Температурные напряжения при тепловом ударе

Так как sh ξ = (eξ e−ξ)2 , то в растянутой зоне

Θ(ξ,ξ)= exp(− ξ)eξ 2e−ξ = 12 [1 exp(− 2ξ)].

При движении фронта внутрь тела η= ξ→ ∞ и Θ(ξ,ξ)→12 , поэтому наибольшее растягивающее напряжение равно:

σ

 

=1,5αT k ;

T

=

р0a

t

 

;

р

 

=

q

.

 

 

 

 

 

 

x

m

m

 

λ

d

 

 

0

 

l

51

Так как ch ξ = (eξ + e−ξ ) 2

то для сжимающих напряжений

при η = ξ

 

 

 

 

 

 

 

 

 

 

 

Θ(ξ,ξ)= −e−ξ((eξ + e−ξ) 2)= −1 2(1+ exp(−2ξ)) ξ > η ,

т.е. их максимальное значение также стремится к 1/2.

 

Оценим порядок

максимальных

напряжений в

стальной

стенке и твердом топливе.

 

 

 

 

 

 

Стальная стенка:

α =1,2 105 1

C ; T

= 40 C ;

μ = 0,3 ;

 

 

 

 

 

 

 

 

 

m

 

 

E = 2 1011 H м2 ; предел прочности – σb =1,2 109 H м2 .

Имеем k =

 

E

 

=

2 1011

=

5

1011 и σx

=1,2 108

H/м2, т.е.

3(1

2μ)

3 0.4

 

 

 

3

 

 

 

 

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

Твердое

топливо

[4].

Для

смесевого

топлива

α =1,2 104 1

C; T =40°C; μ = 0,35 ;

E = 2 107 H/м2. Получим

 

m

 

 

 

 

k = 2,22 107 ,

σx =1,6 105 H/м2.

 

 

 

Для баллиститного

топлива

α =1,2 104 1 C ;

T = 40 C ;

 

 

 

 

 

m

μ = 0,35 ; E = 2 108 H/м2. Получим k = 2,22 108 , σx =1,6 106 H/м2.

4. ГОЛОВНОЙ ОТСЕК В КОНЦЕ АКТИВНОГО УЧАСТКА

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

52

4.1. Устойчивость оболочек корпуса

Наибольшие осевые сжимающие нагрузки в оболочках, которые растянуты на участке входа в плотные слои атмосферы, возникают в конце активного участка. Поэтому для стабилизирующей юбки, например (рис. 21),

(G +G )n0

σ1 = f π1 δ 2 ψx1 , 2 r cos 2

где G1,G2 – веса грузов; n0x1 TGк – осевая перегрузка в конце активного участка; Gк – вес всей ракеты в конце активного участка; ψ2 – угол полураствора конической оболочки; Т – тяга двигателя.

Рис. 21. Нагрузки на головной отсек в конце активного участка

Критические напряжения потери устойчивости слабоконической оболочки ( ψ2 < 20 ), сжатой в осевом направлении, определяются по формуле σкр = (0.15 ±0.35)E (δR2 ) . Тогда коэффи-

циент запаса устойчивости проверяемой оболочки η = σкр / σ1 ≥ 1. Аналогично проверяются на устойчивость и цилиндрические

оболочки корпуса головного отсека.

53