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

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

жении его импульсом давления. Параметры ударной волны можно определить с помощью соотношений, приведенных в [2, 10].

3.2. Напряжения в многослойном корпусе при ударной нагрузке

Рассмотрим многослойный корпус ракеты в виде конструкционной стенки, к которой крепятся несколько слоев теплозащитного материала. Слои изготовлены из изотропных и однородных упругих материалов с различными физико-механическими свойствами. При этом затвердевший клей между слоями может также рассматриваться в качестве одного из слоев композиции [16].

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

∂σx

= ρ

2u

;

εx =

u

;

σx = (λ + 2G)

u

, (3.1)

 

x

x

x

 

t2

 

 

 

 

Рис. 19. Многослойная стенка корпуса ракеты

44

где u – перемещение в направлении, перпендикулярном поверхности слоев; ρ – плотность материала; λ = μG(1)(12μ) –

константа Ляме; G = E2(1) – модуль сдвига; σx ,εx – напря-

жение и относительная деформация в направлении, перпендикулярном поверхности пакета.

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

2u

=

1 2u

,

(3.2)

x2

c2

t2

 

 

 

где c2 = (λ + 2G)ρ – скорость распространения продольных волн

в конструкционном материале.

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

Запишем начальные и граничные условия, для решения системы уравнений вида (3.2).

1.При t = 0 среда находится в покое и поэтому для всех сло-

ев u = 0; ut = 0 .

2.При x=0 (внутренняя поверхность) σx = 0 .

3.При х = δ (наружная поверхность) σx = −p(t) .

4.Напряжения в плоскостях раздела слоев и скорости движе-

ния этих плоскостей одинаковы, т.е.

(σx )i = (σx )i+1 ;

(∂u t = 0)i = (∂u t = 0)i+1 .

 

Волновое уравнение (3.2) гиперболического типа, и его можно решать с помощью метода характеристик. Приведем его сначала к каноническому виду, воспользовавшись следующим пре-

образованием независимых переменных: ξ = ξ(x,t);η = η(x,t) . В новых переменных уравнение (3.2) запишется так:

 

∂ξ 2

1 ∂ξ

2

2u

+ 2

2u

 

∂η ∂ξ

1 ∂η ∂ξ

+

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

c

2

∂ξ

 

 

x

c

2

 

x

 

 

t

 

 

 

 

 

∂ξ∂η ∂x

 

 

t t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂η 2

 

1

 

 

∂η 2

 

2u

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

= 0 ,

 

 

 

(3.3)

 

 

 

 

 

c2

 

 

∂η2

 

 

 

 

 

 

 

 

 

x

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

45

откуда, ввиду гиперболичности исходного уравнения,

 

∂ξ 2

1

∂ξ

2

 

 

 

 

 

= 0 ,

c2

 

x

 

 

t

 

или

∂ξ

1 ∂ξ

= 0 ,

∂ξ

+

1 ∂ξ

= 0 .

(3.4)

x

 

c t

 

x

 

c t

 

 

Записывая систему обыкновенных дифференциальных уравне-

ний, соответствующую (3.4), получаем

dx = −cdt ,

dx = cdt . Ее

интеграл равен: x + ct = const ,

x ct = const. Полагая

 

 

ξ = x ct ;

 

 

η= x + ct

 

(3.5)

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

 

 

 

 

 

 

 

 

 

2u

 

= 0 ,

 

 

 

 

 

 

 

 

∂ξ∂η

 

 

 

 

 

 

 

 

 

 

 

 

 

 

откуда

 

σx

 

 

 

 

 

u = u x + u

 

t

=

 

1 v = α (η) ,

 

 

 

∂ξ ∂x ∂ξ ∂t ∂ξ λ + 2G c

1

 

 

 

и

 

 

σx

 

 

 

 

 

u = u

x

+ u

t

=

 

+

1 v = β (ξ) ,

 

 

 

 

∂η ∂x ∂η ∂t ∂η λ + 2G c

1

 

 

 

или, с учетом выражения для скорости звука,

 

 

 

 

σx −ρcv = α(ξ) ,

 

 

(3.6)

 

 

σx cv (η) ,

 

 

 

гдеv = ∂ut – скорость движения частиц среды. Таким образом, уравнения характеристик имеют вид

dxdt = c ; dxdt = −c ,

а вдоль них выполняются следующие дифференциальные условия:

dσx −ρcdv = 0 , (ξ = const) ; dσx cdv = 0 , (η = const) .

Используя выражение (3.6), получаем

σx = (α+β)2 , v = −(α−β)2ρc ,

где α и β – константы интегрирования.

Расчет начинается с построения сетки характеристик для каждого из слоев композиции в плоскости x, t с помощью (3.5). Точки пересечения линий ξ=const и η=const служат узловыми,

46

в которых определяются напряжения и скорости движения частиц среды.

Если в узле известны σx и v , то формулы позволяют определить константы α(η) и β(ξ) , постоянные вдоль линий ξ=const и η=const соответственно, сходящихся в этом узле. Далее по известным значениям α и β в двух узлах можно определить σx и v в соседнем узле и т.д. При этом, конечно, используются соответствующие начальные и граничные условия.

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

Возникновение напряжений, превышающих допустимые, нежелательно из-за cнижения теплозащитной и несущей способности конструкции.

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

Тепловым ударом называется быстро изменяющееся во времени внешнее воздействие, создающее значительные градиенты температуры по толщине конструкции и напряжения. Будем предполагать, что материал сохраняет свои упругие свойства [20].

При тепловом ударе в материале возникают нежелательные эффекты: растрескивание, раскалывание и крошение, отделение слоев многослойной конструкции.

Для составления математической модели воспользуемся следующими допущениями:

тепловой поток q постоянный и действует в течение малого промежутка времени 0 < t td;

все поступающее тепло аккумулируется в слое материала толщиной l;

объемные источники тепла распределены в слое l по экспоненциальному закону:

p(x,t) = (q l)exp(x l)[H (t) H (t td )],

0,t < 0 ;

где функция Хевисайда H (t) =

1,t 0

47

– до начала теплового воздействия температура конструкции равна нулю.

Уравнения термоупругости для одномерного деформированного состояния записывают в следующем виде:

1) уравнение равновесия

∂σx

= ρ

2u

, где σx ,u – напряже-

 

 

x

 

t2

ния и перемещение в направлении нормали к поверхности; ρ – плотность материала;

2) геометрические уравнения εx = ∂u x ;

εy = 0 ; εz = 0 , где

εx ,εy ,εz – относительные деформации в

направлении осей

x, y, z ;

 

3)закон Гука:

εx = E1 x −μ(σy + σz )]+ αT (x),

0 = E1 y −μ(σx + σz )]+ αT (x), 0 = E1 z −μ(σy + σx )]+ αT (x),

где α,T – коэффициент линейного расширения и температура; E,

μ– модуль упругости и коэффициент Пуассона;

4)уравнение теплопроводности в случае объемных источников тепла:

ρc

T

=

 

λ

T

+ p(x,t) ,

t

 

 

 

 

 

 

x

 

x

 

или, если принять постоянным коэффициент теплопроводности,

T

= a

2T

+

p(x,t)

,

(3.7)

t

x2

ρc

 

 

 

 

где коэффициент температуропроводности

a = λ ρс.

Преобразуем сначала уравнение равновесия, выразив σx , через εx из физических уравнений:

σx −μσy −μσz = Ex − αT );

μσx + σy −μσz = −EαT ;

μσx −μσy + σz = −EαT .

Вычисляем определитель системы уравнений

48