книги / Общая теория анизотропных оболочек
..pdf422 |
КОЛЕБАНИЯ И УСТОЙЧИВОСТЬ ОБОЛОЧЕК |
[ГЛ. 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,
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)
§ 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) EQ— E1( 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).
Считается, что магнитные и диэлектрические проницаемости газа и материала оболочки изотропны и равны единице. Внутри )болочки предполагается справедливость уравнений Максвелла цчя вакуума. Электромагнитные эффекты в теле оболочки вовсе зе рассматриваются.