Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Общая теория анизотропных оболочек

..pdf
Скачиваний:
2
Добавлен:
12.11.2023
Размер:
21.14 Mб
Скачать

§ 5]

ПОЛОГИЕ ОБОЛОЧКИ В ТЕМПЕРАТУРНОМ ПОЛЕ

421

Полагая

82^> 1, что имеет место для многих прикладных за­

дач, применим к уравнению

(5.20) метод

асимптотического

ин­

тегрирования.

 

 

 

 

 

Интегралы уравнения (5.20)

будем искать в виде (в дальней­

шем, где это допустимо, индексы т и п опускаются)

 

 

/ ( « )=

Ф

( *8)е8“,

<*>,

(5.24)

где

Ф (f, 8) = Ф0 (0 +

8-1Ф, (t) +

 

 

 

 

8_2Ф2 (*) + ...

 

— функция

интенсивности,

со (t)

— функция изменяемости.

 

Подставляя значение / (t)

из

(5.24)

в

(5.20) и производя

со'

кращение на е8ш, получим

 

 

 

 

 

8’4 * +

Ш

!] + s ( 2

 

 

 

 

 

 

 

+2B-W{ d/2

, d$j+i dta

d2u> ,

 

 

 

'

dt

dt"

 

 

 

 

J=0

 

 

 

 

 

 

 

 

+

®

л . [ + + ( у ) ' ] } = 0-

(5-25)

Потребовав, чтобы в уравнении (5.25)

коэффициенты при всех

степенях

8 обращались в нуль, при условии Фо=^=0 получим сле­

дующую бесконечную систему рекуррентных уравнений для опре­

деления со (f), Ф 0 (f),

Фл (t), . . .,

Ф .(f), . .

 

 

 

 

 

 

 

 

 

 

 

(5.26)

 

0

<гф0

dco

, л

 

<г2<в

п

 

(5.27)

 

^

dt

dt

1 Ф<>

— 0»

 

d4iJ

,

о <,Ф>+1

doi

I

л

dti ~

n

(5.28)

d«2

1 “

dt

dt

1 У+1

U

 

 

0 =

1,

2, 3,

...).

 

 

 

Решив (5.26), для функции изменяемости получим

*i = t \ \/tydt-\-C'v

(5.29)

со2 = — г j y/tj) сЙ-j-C'.

Внеся (5.29) в (5.27), найдем функцию интенсивности в нуле­ вом приближении:

ф 0 = с д н /*.

(5.30)

422

КОЛЕБАНИЯ И УСТОЙЧИВОСТЬ ОБОЛОЧЕК

[ГЛ. III

Имея (5.29) и (5.30), с помощью (5.28) можно найти следую­ щее, первое приближение. Каждое последующее приближение также найдется путем интегрирования уравнения (5.28).

Ограничимся рассмотрением нулевого приближения, ибо по­ дробные исследования показывают, что нулевое приближение в ряде случаев обеспечивает вполне удовлетворительную точность результатов.

Согласно (5.30), (5.29) из (5.24) в нулевом приближении получим

 

/ =

 

ехр^гВ j у/ф

+

ехр (— $ J

 

(5.31)

Постоянные интегрирования С

входящие в (5.31), с помощью

(5.18)

могут

быть

определены

из начальных

условий (5.17),

В

силу (5.31)

из (5.19) для

второй искомой функции имеем

F = — ( ~ ! г )

{ к ^ п 1-j- k2b2m 2)

exp (ia J \/<i dt^j -f-

 

 

-f- C2 exp ^— ib j

[А 22Ь*т* -j- А п а*п* -f-

 

 

 

 

+ {Ат +

2А1,)а ^ т Ч Г 1.

(5.32)

Имея значения

искомых функций

/ (t) и F

(t), с

помощью

(5.18)

легко

записать значения функции напряжений <р и функ­

ции перемещения w.

Таким образом, поставленная задача решена. С помощью фор­ мул (5.18), (5.31), (5.32) легко найти характер напряженно-

деформированного состояния оболочки,

когда температура ее из­

меняется во времени по произвольно заданному закону Т = Т ( t ) .

Рассмотрим случай, когда средняя температура оболочки из­

меняется по закону

 

 

T = Z fч± t ( 0 < * < * , ), Т = т ж ( f > g ,

(5.33)

где Т тлх — максимальная температура

оболочки,

tx — время,

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

температуры модули упругости материала оболочки зависят от температуры линейно или хорошо аппроксимируются линейной функцией, т. е.

Е{ = Е%(1 — Х'Т), £o = £ . ( g = const.

(5.34)

Приближенно будем считать, что и коэффициенты B(j. зависят от температуры линейно:

Вч == B<lj (1 - Х'Т), В?,. = Вч (д = const,

§ 5]

ПОЛОГИЕ ОБОЛОЧКИ В ТЕМПЕРАТУРНОМ ПОЛЕ

423

где

(5.35)

Тогда из (5.20) после некоторых преобразований, ограничиваясь основным тоном колебаний (т = п = 1), получим

^ + с’ (1 - Х *)/ = 0,

(5.36)

где

С2 _ в № 12ра*

 

(5.37)

X = ¥T,aJ t v |Х<|<1.

(5.38)

Укажем, что из физического содержания рассматриваемой за­ дачи следует, что 1—Xt > 0.

Уравнение (5.36) путем элементарной подстановки приводится к уравнению Бесселя и может быть интегрировано точно. Выпол­

нив интегрирование, получим

 

 

/ = ClV/ l- X * / Vs[ g ( l - X

f ) s/’] +

 

 

 

+

c 2s!1 -

и /-%[■§!■ (1 - X*)s/’] ,

(5.39)

где /v — функции Бесселя,

 

 

Постоянные интегрирования С,, подлежат определению

из на­

чальных

условий задачи

(5.17),

которые в данном случае

могут

быть записаны в виде:

 

 

 

при

t = 0

 

 

 

 

 

 

 

(5.40)

Удовлетворяя условиям (5. 40), получим

 

(5.41)

424

КОЛЕБАНИЯ И УСТОЙЧИВОСТЬ ОБОЛОЧЕК

[ГЛ. III

Подставляя значения С{ в (5.39), получим выражение для / (<) при т = 1, п=1 . Подставляя это выражение в (5.18), найдем нор­ мальное перемещение оболочки w (а, р, t):

Ш

к{ wo IH (ад Н ®

+ 1-'1. ( У i -ч, Ю 1-

 

- Т

f7-'/» ( У Н (С) -

Н (ад /-V, (Q]} sin т sin У *

(5-42)

где для сокращения записи введены следующие обозначения:

Таким же образом можно записать и выражение для искомой функции (р ( и, р, t).

Легко установить, что с увеличением t (т. е. с повышением температуры оболочки (5.33)) увеличиваются как амплитуда, так и условный период колебания нагреваемой оболочки. Обрат­ ную картину можно наблюдать при ее остывании, т. е. при повы­ шении жесткости оболочки.

Во многих прикладных задачах аргументы бесселевых функций, входящих в (5.42), имеют большие значения. Используя извест­

ную асимптотическую

формулу

 

 

 

 

 

 

 

J> ) =

/ ! ; C0S( e “ ” T

- T

) ’

 

из (5.42)

получим

 

 

 

 

 

 

 

w=

k[ш0 cos (С0 —

С) + ^ sin (С0 —

С)] sin ^ sin f .

(5.43)

Из (5.43)

найдем

перемещение

центра

оболочки

(а = а/2,

Р = 6/2):

 

 

 

 

 

 

 

 

 

 

1» ц =

/= 0 (1 -

>4)-Ч. sin { 1

11 -

(1 -

X ff’I - f » } ,

(5.44)

где

 

 

 

 

 

 

 

 

 

 

 

® = V

wo + ( 7 ) 2. . & =

aretg ^ r .

 

При

X=0 из (5.44) предельным переходом нетрудно получить

общеизвестное уравнение

гармонических

колебаний

 

 

 

 

и?щ= й> sin (ct -f- &).

 

 

(5.45)

Простой анализ показывает, что приближенное решение, полученное в нулевом приближении при асимптотическом интег­ рировании уравнения (5.36), совпадает с точным решением (5.42)

8 5]

ПОЛОГИЕ ОБОЛОЧКИ Б ТЕМПЕРАТУРНОМ ПОЛЕ

425

при асимптотическом представлении входящих в него бесселевых функций.

При наличии линейного затухания (с коэффициентом затуха­ ния е) взамен (5.36) будем иметь следующее уравнение:

 

 

& ' + 2 . « + е Ч 1 - Х 0 / = 0 ,

(5.46)

 

 

dPT

dt

 

 

решив

которое,

для

w ( и,

|3, t) получим

 

w = e -

i 2V/( Cl _

x 0

_2 { eC l

/ 1/a [ |

( l \t- - | ) 3/l] +

 

 

+

[Ii1 — xt -

Sf]}sinTsin?•

(5.47)

He вдаваясь в подробности, укажем, что в формуле для w имеется затухающий множитель, который, в зависимости от соот­ ношений величин в, с, X, может существенно изменить характер колебания оболочки. В частности, на­

пример (в

отличие

от

случая

е=0)

U

значительное повышение

температуры

 

может и не привести к увеличению

 

амплитуды

колебания.

 

 

 

 

2. Устойчивость при обтекании ор-

 

тотропной

цилиндрической

оболочки,

 

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

 

ной

■температуры. Рассмотрим

круго­

 

вую

цилиндрическую

оболочку

беско­

 

нечной длины, обтекаемую сверхзвуко­

 

вым потоком газа с невозмущенной

 

скоростью

U, направленной

по обра­

 

зующим оболочки (рис. 78). Будем полагать, что оболочка изго­ товлена из ортотропного материала, главные направления упру­

гости

которого совпадают с координатными линиями оболочки

а» Р>

Т-

Механические характеристики оболочки зависят оттемпературы, являющейся функцией лишь времени t, т. е. Т —Т (t).

В рассматриваемом случае исходные уравнения задачи, со­

гласно (5.12)

и

(4.2), в линейной постановке запишутся следую­

щим образом:

 

 

 

 

 

 

 

 

С 11

. / 1

о Ci2\

I

^22

1

(Y'w

^

Q0 da* ~|

\С 66

Q 0 ) д а Щ 2 ~ * ~

Q0

JR

d a*

O u S + 2 <c . . + 2 D » > ^ | i + ^ $ + s4 -? g + + i

(5.48)

 

+ ( 2^ + ' S ) i + %S u ^ = ° -

426

КОЛЕБАНИЯ И УСТОЙЧИВОСТЬ ОБОЛОЧЕК

[ГЛ. III

где,

как всегда,

 

 

 

 

Ca = Bik®h,

Dik= B{k(t)^2, Q0=

CuCB -C J|f

(5.49)

p—'fo/g — плотность

материала оболочки

(у0 — удельный вес,

g — ускорение силы тяжести), е — коэффициент затухания, рот —

давление невозмущенного

потока газа,

ат — скорость

звука

в невозмущенном потоке,

х — показатель политропы.

 

Исключая <р (а, р, t) из системы (5.48),

получим для

опреде­

ления искомого нормального перемещения w (а, р, t) следующее уравнение:

 

 

(5.50)

Решение уравнения (5.50) ищем в виде

 

 

 

(5.51)

Здесь / (t)

— некоторая комплексная

функция действительного

аргумента,

— волновое число,

X — длина полуволны по

образующим оболочки, п — целое число волн по окружности по­ перечного сечения оболочки.

Подставляя значение w из (5.51) в (5.50), для определения / (t) получим следующее дифференциальное уравнение второго порядка с переменными коэффициентами:

где

( 5 . 5 3 )

«Рсо __%Paag

( 5 . 5 4 )

fhaga 1о^аоо'

Подстановкой

( 5 . 5 5 )

§ 5] ПОЛОГИЕ ОБОЛОЧКИ В ТЕМПЕРАТУРНОМ ПОЛЕ

уравнение (5.52) приводится к виду

g 4 + w ( « ) * = o ,

где введены следующие новые обозначения:

 

 

 

»'

Ei (0)

 

 

 

12^0*(1 — VjVg) *

VИ\

El (t) hSgГА

,

12(1—VjVg) Д2 l»*"l

г w

El (0) Ri

i

**

д;J

12 (1 — Vj.v2 ) *To Г,-m .. тт I /. i (Л21 m= kR.

Здесь

427

(5.56)

(5.57)

(5.58)

A„ = ™ ‘ + 2[»< + 2| f § ( 1 -

E2(0 „4

Ei <»)"’

 

 

(5.59)

^ = ^ + [ 4 ( 0

2vi]m2« 2 +

« 4.

Ei (0) — модуль упругости главного направления 1—1 в началь­ ный момент £ = Со­

множитель второго члена уравнения (5.56) является боль­ шим параметром. Учитывая сказанное, уравнение (5.56) будем решать методом асимптотического интегрирования; решение урав­

нения (5.56)

будем искать

в виде

 

 

ф =

ф0(8, С)е8“ “ \

(5.60)

Фо(8- t) =

Ф0 (t) + П Ф Х(С) + 8-*Ф2 (< )+ • • • . Ф0(С )# 0 ,

(5.61)

где Фо — функция интенсивности, со (t) — функция изменяемости. Не вдаваясь в подробности, освещенные в предыдущем пункте настоящего параграфа, согласно (5.55), (5.56), (5.60) и (5.61)

в нулевом приближении для / (С) получим

/(f)=F-7.

 

(5.62)

 

t

 

РЮ = 8

+

(5-63)

 

о

 

где C1= a 14-i61, C2=a2-\-ib2— произвольные

комплексные по­

стоянные.

Имея значение / (С), по сути дела, можем считать, что постав­ ленная задача решена. Исследуя поведение функции / (С), мы мо­ жем установить значение критической скорости потока газа.

Характер функции / (С) существенно зависит как от скорости набегающего потока U, так и от коэффициентов упругости BiJe (Т)=В'.к (t), зависящих в свою очередь от температуры. Если

428

КОЛЕБАНИЯ И УСТОЙЧИВОСТЬ ОБОЛОЧЕК

1ГЛ. III

модуль амплитуды / (t)

в течение интересующего нас промежутка

времени

возрастает, то

будем считать, что оболочка

находится

в состоянии неустойчивости. В противном случае состояние обо­ лочки считается устойчивым. Скорость потока газа U*, при ко­ торой d\f (t)\/dt=0, является критической.

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

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

(Ян (t)= E 2 (t)= E (f), Gla (t)= E (t)/2 (1 + v), v1 = v2= v ) с

модулем

упругости E (T), линейно зависящим от температуры:

 

Е = Е0-Е * 0Т,

(5.64)

где Е0 и Ео — коэффициенты, определяемые из эксперимента. Будем полагать также, что средняя температура оболочки

изменяется по закону

 

Т =

Т = ТШ* (*>*!> .

(5.65)

гДе

— максимально достигаемая температура оболочки,

время, необходимое для достижения максимальной температуры. В силу (5.64) и (5.65) окончательно получим

Е = Е9- Е { ,

=

(5.66)

 

 

С 1

Тогда из (5.62) и (5.63), с учетом (5.57), (5.58) и (5.66), найдем цля / (£) и р (£) следующие выражения:

f(t) = |->/^-Я1«)-('+«^)<[(с08

 

 

 

 

 

9 +

2тс/'

 

 

 

 

 

 

3& +

2STC

, .

39 +

2STC'

 

 

 

>>]•

 

 

 

 

'«)= -5[г'Ч'COS------------ i Sin-----Г

)

-

 

 

 

 

 

390 +

2sn

. .

 

3&o +

2s7t

— r0/j( cos —— ---------i sin

 

2

 

 

2

 

 

 

 

) ] •

&(<) = - | [ r5/,cos

39 -f- 2*it

 

 

3&o + 2sn'

 

 

 

------2

]•

j, .

39 +

2sit

8/

. 390 4- 2STC"]

p.w=*£|jr h Sin--- 2 --------r0'>sin - 2-~

-----

 

(/ =

0,1, 2,3;

s =

0,1),

 

 

 

(5.67)

(5.68)

(5.69)

§ 5] ПОЛОГИЕ ОБОЛОЧКИ В ТЕМПЕРАТУРНОМ ПОЛЕ 429

где

ГЩ): ■Е ф [ВЩ) + А*и%

г0 =

г(0),

 

 

 

 

® (0 =

arctg

ли

«о

а (0),

А = 12уоЛ (1 — V2) тр,

 

,

(5.70)

д __EL

 

 

, 12(1 — у2)

т*

I

(тг 4 -га2)2 4 - 12Дг (1 ~

v2)

mi

 

 

а _ _ Я0 Я* I/

' П '

'

Л2

(j»2 +

 

n2)2j*

 

В (£) =

 

 

-

12ТоЛЯ (1 -

V2) (в +

1 ) 2.

 

Как было указано выше, критическая скорость определяется

из условия d\f (f)|/df=0,

которое в развернутой форме,

согласно

(5.67)—(5.70),

имеет

вид

 

 

 

 

 

 

 

{ г

~

[(е%(0 [д2 cos

_(_р2 sin p j

 

 

 

-f- Oj cos

— 6j sin PJ)2 -f- (e^(0 [fc2 cos PJ a2 sin p j -f-

 

 

 

 

&i cos Pj -f-« i sin P,)2]l/l} = 0.

(5.71)

Таким образом, мы получили уравнение, при помощи которого можно определить точное значение критической скорости U*. Однако, ввиду громоздкости уравнения (5.71), аналитическое определение U* затруднительно. Но если учесть, что модуль пер­ вого слагаемого функции / (£), т. е.

\ сА г~‘и е х Р [р2 (0 — (е + рУ2) £],

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

д

{г-'1ехр 1 -р 2 (£) -

(е +

р/2) £]} =

0.

(5.72)

Далее, имея в виду,

что при 0 ^ £ ^ £ г

 

 

 

B(t)> 0 ,

M

=

_ 8rv2Si n

| ,

0 < »

< | ,

 

р2 ( £ ) < 0 ,

из (5.72) получим

 

 

 

 

 

 

 

В ( 1 ) E QR CI

\ / V g a + A W * В

 

 

(5.73)

4(й2 +

Л2С/2)

 

(1 — v'2)

 

 

 

 

 

 

 

 

 

Отсюда легко определить значение критической скорости по­

тока газа U* в любой момент времени t (0 ^

t ^

£Д.

В большинстве случаев В2

A 2U2; тогда

из

(5.73) можно

получить

 

££= 1 _______1________ £,

 

 

 

 

 

 

 

(5.74)

 

 

Vо

4 (е+ JJ./2) EQE1( t ) ’

 

 

 

 

 

430

КОЛЕБАНИЯ И УСТОЙЧИВОСТЬ ОБОЛОЧЕК

[ГЛ. III

где

 

 

 

 

 

 

77* — 1 Г £ Ш

/а Г(т2 + » 2)2

, 12Д2 (1 - у2)

m2 17;(,

, 2е\

0

Д L Р* J

L m2

л”

А 2

( m 2 + n2)2J \

*” (j./

— критическая скорость флаттера, найденная по квазистатической с точки зрения температуры теории (т. е. когда в окончатель­ ной формуле критической скорости для несогреваемой оболочки вместо модуля упругости Е подставляется выражение ^-Е ^)).

Формула (5.74) показывает, что динамическая критическая скорость U* меньше критической скорости, найденной по квазистатической теории, и в случае достаточно малого затухания мо­ жет существенно отличаться от квазистатической скорости U*0.

§ 6. Флаттер цилиндрической оболочки в потоке сжимаемой проводящей жидкости в присутствии магнитного поля

Рассмотрим задачу устойчивости анизотропной круговой ци­ линдрической оболочки бесконечной длины, обтекаемой с внешней стороны сверхзвуковым потоком идеально проводящего невязкого газа с невозмущенной скоростью U, направленной по образующим ци­ линдра. Пусть во внутренней об­ ласти оболочки газ находится в по­

кое и имеет давление рт,

равное

давлению

невозмущенного

потока

обтекающего

газа.

Оболочка

нахо­

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

на­

пряженности

которого равен

Н 0

и направлен

параллельно

скорости

потока (рис. 79).

 

 

 

 

Пусть

оболочка

изготовлена

из

анизотропного материала,

который

имеет лишь одну плоскость упру­

гой симметрии, параллельную срединной поверхности

оболочки.

За координатную поверхность у= 0 (r= R ) принимается средин­

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

а,

р ( а —

вдоль образующей, р — по дуге окружности поперечного сечения). Система координат а, р выбрана так, что А =1, B = l, R1= co, RZ= R (рис. 79).

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