
книги / Оболочки и пластины
..pdf
in R |
• |
/ |
где п = ----- , |
причем i |
и k принимают целые положительные значе |
ния. Здесь i показывает число полуволн в меридиональном направле нии (/— длина оболочки), k дает число полуволн в окружном направ лении. Решение в форме (4,12,6) соответствует тому случаю, когда на
концах оболочки |
(z = 0 и z = l) исчезают как радиальные, так |
и окруж |
|
ное перемещения |
(иФ 0). |
что |
уравнения |
Легко убедиться непосредственной подстановкой, |
|||
(4,12,5) тождественно удовлетворяются, если функции |
U(t), |
V(t), W(t) |
определяются из системы обыкновенных дифференциальных уравнений: |
|
|||||||||||
_ т ( 1 - ^ |
J |
^ |
+ ( n t + ± J Lk, \ u _ |
l ± |
v nkV - v n W - |
|
||||||
Eh |
|
dt* |
|
\ |
|
2 |
) |
2 |
|
|
|
|
|
|
— |
1~ * |
--1~ Тг n(kV + W) = |
0, |
|
|
|||||
|
|
|
|
Eh |
|
R * |
K |
/ |
|
|
|
|
- m ( 1 |
~ |
v2) |
|
E t—z |
пш |
+ —fk2 + |
2 |
|
) |
v + m л +Л |
|
|
Eh |
|
dt2 |
|
|
2 |
|
\ |
|
|
|
||
|
|
|
|
1- V |
2 |
JTi-Tg nkU = Q |
|
|
(4,12,7) |
|||
|
|
|
|
Eh |
R |
|
|
|
|
v |
’ |
|
m |
(1 — v2) |
|
dPW |
— vnU + |
kV + c2 {n2 + k2)2W — |
|
||||||
|
|
|
|
dt* |
|
|
|
|
|
|
|
|
|
|
- |
|
|
1 7 |
n « + |
7 (jfe*_ |
1 )] |
= |
0 . |
|
|
|
|
|
Eh |
R 2 1 |
1 |
v |
/J |
|
|
|
|
|
Система уравнений (4,12,7) может быть представлена в матричной |
||||||||||||
форме |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
tn~7T + ^ - T i S 1- T 2S2)f = 0, |
|
|
||||||||
|
|
|
dt2 |
|
|
|
|
|
|
|
|
|
где / — вектор с компонентами U, V, tt7, |
|
|
|
|
|
|||||||
|
|
|
п2 + |
1~ v k2 |
— L+-X. nk |
|
— vn |
|
||||
|
|
|
|
|
2 |
|
2 |
|
|
|
|
|
|
|
|
|
2 |
|
|
2 |
|
|
|
|
|
|
|
|
|
— v/г |
|
k |
с2 (n2+ |
/г2)2 |
|
|||
|
|
0 |
|
n k |
n |
|
|
0 |
|
— n k |
n |
|
S i= — |
— n k |
0 0 |
S ,= A n k |
|
0 |
0 |
|
|||||
R 2 |
0 |
|
0 |
|
|
R 2 |
0 |
|
0 |
k 2— 1 |
|
|
|
|
|
n 2 |
|
|
|
|
Частоты собственных колебаний незагруженной оболочки опреде ляются из уравнения
|£ — р2£ | = 0,
а критические параметры продольной сжимающей и радиальной нагру зок — из уравнений
= 0 и \R — qRS2\ = 0.
где
(я2 + & + |
I)2 (/X2 - f А?)2 + (1 — V) л 2 ( л 4 — б4) + |
- — |
с2 |
л4 |
|||
g (п, k) = |
|
|
|
|
|
(4,12,13) |
|
|
( л 2 + |
А2)2 |
|
|
|
||
|
|
|
|
|
|
||
.Введем обозначения: |
|
|
|
|
|
|
|
D g (n ,k ) = |
„ |
P g (n , Л) |
= ^ ^ |
D g (я, |
k) |
|
(4,12,14) |
mR4 |
|
я2Я2 |
|
(А2 — I)/?2 |
|
||
|
|
|
|
||||
Уравнение (4,12,12) может быть теперь записано в виде |
|
|
|||||
*L |
|
|
|
|
|
(4,12,15) |
|
dt2 |
|
|
|
|
|
||
|
|
|
|
|
|
||
где pt — 2nRTi», |
R |
Следовательно, |
задача сведена к извест- |
||||
|
|
|
|
|
|
|
ному уравнению.
В заключение отметим, что эти результаты легко обобщаются на случай ортотропной цилиндрической оболочки. В самом деле, пусть Е1, Е2 и V I , V2 — модули упругости и коэффициенты Пуассона в направ лениях p= const и a = const соответственно; G — модуль сдвига. Введем следующие дифференциальные операторы:
|
|
|
|
2 |
|
г |
а 2 |
|
, |
„ |
|
а 2 |
|
|
|
|
|
|
V° ~ |
0 да2 |
^ |
2 |
ар2 |
|
|
|
|
||
|
|
|
|
а4 |
|
+ £4 - |
* |
|
I |
£ — |
|
|
|
|
|
|
|
|
да4 |
|
|
да2Эр2 |
|
|
2 Эр4 ’ |
|
(4,12,16) |
||
|
|
|
|
д* |
|
д4 |
|
+ £ |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|||||
где |
|
|
|
да4■+ £3 |
da2 ар2 |
|
|
2 Эр4 |
’ |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
£ 0 — 2G (1 — |
|
“Ь ^ 2 v i> |
|
|
|
||||||
|
|
|
£ 3 = |
4G (1 — v xv 2) - f £ XV 2 |
+ |
£ 2V „ |
|
(4,12,17) |
||||||
|
|
|
£4 |
|
G |
£ I V2 — |
£ 2VI . |
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Уравнение, |
аналогичное |
(4,12,11), |
для |
ортотропной |
оболочки прини |
|||||||||
мает |
вид |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4у4ф |
, |
£A (l-viv«) |
. *Ф_ |
^ |
|
|
|||||
|
|
|
V‘V* |
+ |
С2 |
|
|
|
да4 |
|
|
|||
|
|
j ^ - ^ z |
^ |
|
^ |
|
|
|
|
|
|
0 |
(4 12,18) |
|
|
^ |
с2л |
\ 1 |
За2 |
|
ар2 |
|
|
|
а/2 у |
Vl |
|
v |
|
где |
Ф — функция, |
связанная |
с радиальным |
перемещением |
соотноше |
|||||||||
нием (4,12,9). |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
Дальнейшие преобразования вновь приводят к уравнению (4,12,15), |
коэффициенты которого для случая опертой оболочки определяются без труда [24].
Составляющие |
АХ и АУ равны нулю. Кроме внешнего давления (4,13,1) |
|||
и приведенной |
нагрузки (4,13,5) на |
оболочку |
действуют |
силы инерции |
|
— т ■-----д2и , — т ----- |
d2v , — т |
------dPw . |
(4,13,6) |
|
di2 |
dt2 |
dt2 |
|
Отбрасывая в согласии с принятым приближением тангенциальные со ставляющие сил инерции, найдем, что силы, действующие на оболочку, сводятся к радиальной нагрузке
Z = — (q0 + qtcos 00 — (<7о + qtcos 00 (V2 + 2)w — m |
. (4,13,7) |
Первое слагаемое дает равномерное сжатие оболочки и может быть отброшено, если под w (г|э, р, t) понимать отклонение от невозмущенного безмоментного Состояния. Тогда уравнение (4,13,2) принимает вид
Jc2(V 2 + I)2 + 1] ( V 3 + |
2) w + |
J?i±i^£L0Otf_ (V2 + 1 _ v) (v* + |
2) а, + |
||||
|
mR2 |
/ о |
, |
, |
ч d2w |
п |
(4,13,8) |
|
Eh |
(V |
+ |
1 — V) —— = |
о. |
||
|
|
|
|
dt2 |
|
|
|
Будем искать решения уравнения |
(4,13,8) в классе функций |
|
|||||
|
ш (ф ,М ) = /(0^(Ф ,Р), |
|
(4,13,9) |
||||
где f(t) — неизвестная |
функция |
времени; Р(ф, р)— решения |
дифферен |
||||
циального уравнения |
|
|
|
|
|
|
|
|
|
V*F + XF = 0, |
|
(4,13,10) |
удовлетворяющие граничным условиям для до (т. е. условиям непрерыв ности и однозначности на сфере). Подстановка в (4,13,8) дает после
•сокращения на /■'(ф, Р) |
|
|
[с2 (X — I)2 + 1] (К— 2) / — - (t?a + q‘cos 90 R |
(k — 1 + |
v) (A, — 2) / -f- |
2Eh |
|
|
+ ^ . (X _ 1 + v ) - g - |
= °. |
(4,13,11) |
Введем обозначения: |
|
|
* • = R ( t l ? + - i c l <x - 1>a + 4 |
<4 ’ 1 3 ' 12) |
и, внеся их в (4,13,11), получим
----- ?0 + ftyos,6'-. у = о, |
(4,13,13) |
Формулы (4,13,12) дают собственные частоты и критические силы, за висящие от неизвестного еще параметра К. Впрочем, один практически важный вопрос может быть решен до конца и без определения К. Гра-
ницы главных областей неустойчивости могут быть найдены по извест ным приближенным формулам. В частности, низшая граница
02 = 4р211 — ь + Т * |
(4,13,14) |
Яо |
|
Для практических приложений интересно знать огибающую ниж них границ областей неустойчивости. Предположим, что параметр X мо жет принимать любые вещественные положительные значения, т. е. до пустим, что уравнение (4,13,10) имеет сплошной спектр собственных значений. Для ограниченной области, какой является сфера, спектр соб ственных значений дискретный. Но в окрестности интересующих нас значений А, спектр уравнения (4,13,10) все же достаточно «густ», чтобы погрешность, следующая из сделанного допущения, была невелика.
В дальнейшем обозначим
|
|
|
02 = |
q(X). |
|
|
(4,13,15) |
||
где |
|
|
mR* 4 V ’ |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
q(k) = ( X - 2) с2 (A, — l)2 + 1 |
|
(<7o+ |
2 |
(4,13,16) |
||||
|
|
|
|
2Eh |
|||||
|
|
|
X — 1 4- v |
|
|
|
|
||
Для |
определения |
огибающей |
положим |
dq |
= 0, отсюда |
получим |
|||
уравнение для X. Рассмотрим случай |
|
|
dX |
|
значений |
||||
достаточно больших |
|||||||||
Л>1. Тогда |
|
|
|
|
|
|
|
|
|
|
q (X )^ c W |
+ t |
( ‘7в+ |
2 |
q ) R |
(4,13,17) |
|||
|
|
|
|
|
|
2Eh |
|
|
|
|
|
dq |
|
|
|
|
|
|
|
откуда корень уравнения -2 - = 0: |
|
|
|
|
|
|
|||
|
|
dX |
|
|
|
|
|
|
|
|
|
|
(яо+ 2 |
* ) |
' |
|
|
(4.13.18) |
|
|
|
|
|
4Ehe* |
|
|
|
||
и, следовательно, |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
||
|
02: |
4£/t |
1 — |
( . + - Н * R * |
(4.13.19) |
||||
|
mR2 |
|
16£^/i2c2 |
|
|||||
Введем обозначения: |
|
|
|
|
|
|
|
|
|
|
|
4Ehc |
= Я - |
Eh |
|
|
|
|
(4.13.20) |
|
|
R |
mR* |
|
|
|
|
||
|
|
|
|
|
|
|
|
||
Формула |
(4,13,19) принимает вид |
|
|
|
|
|
|
||
|
02=4р2 |
|
|
|
|
|
|
(4.13.21) |
Конечно, следует учитывать, что q** представляет собой приближенное (в смысле сделанных допущений) значение минимального критического давления. Действительно, учтя, что
с2 = — —Л*——
12/?2 (1 — v2) ’
получим хорошо известную формулу
2Eh? |
1 |
(4,13,22) |
|
R* |
' у/ 3 (1 — v2) |
||
|
Наконец, займемся определением параметра К в общем случае. Уравнение (4,13,10) на сфере
_ I |
_ |
p |
L |
дф , |
kF = 0 |
(4,13,23) |
sin “ф |
[_ |
дф |
\ |
sin2 г|) d(i2 |
|
приводит, как известно, к сферическим функциям. Решение уравнения (4,13,23) находим в виде
К ( ф , Р) = Р ( ф ) |
s in ftp |
(4,13,24) |
|
cos ftp.
Условие однозначности Р(ф, р) на сфере требует, чтобы ft было целым числом или нулем (ft = 0, 1,2...). Подстановка в (4,13,23) дает
— — |
Г — — ( s i n ф ) |
] + |
( ъ |
& ) р = 0. |
(4,13,25) |
sin гр |
iL avdV |
dtyагр Jj |
\ |
sin2 тр / |
|
Полагая х = созф, приводим уравнение (4,13,25) к виду
г[<1-ч-£Кх-гг5)я“0- <4’13’26)
Это известное уравнение для присоединенных полиномов Лежандра. Уравнение (4,13,26) имеет собственные значения
Х„ = п(п+ 1) |
(п = 0 ,1 ,2 ,...) . |
(4,13,27) |
Каждому собственному значению Я„ соответствует (га+1) |
собственных |
|
функций: |
|
|
р*(*) = (1 _ * 2)* |
~~ |
Рп{х) (ft = 0, 1,2, ...,п), |
(4,13,28) |
|
где |
dxk |
|
|
|
|
|
dn |
|
|
Рп(*) |
|
|
|
|
2Пп\ |
dxn |
|
||
Теперь возможно записать |
для |
уравнения (4,13,23) |
систему его |
|
решений: |
|
|
|
|
ft = 0 |
Р0(ф, Р) = Pn(cosij3), |
|
||
ft = 1 |
Р_!(ф, Р) = PJf> (cosф) sin Р, |
(4,13,29) |
||
|
р1(ф, Р) = Р„ *(cos ф) cos Р, |
|
»
»»
k = п F-n(ф, Р) = Я л) (cos ф) sin /гР,
Fn(ф, P) = Р%] (cosy) cosnp.
Известно, что полиномы Лежандра Рп(х) имеют в интервале измене ния ф(0, я) ровно п нулей. Присоединенные функции P {nk) (х) имеют со
ответственно (n—k) нулей. |
в нуль на |
2k |
меридианах, |
|
Так как sin&p и cos&p обращаются |
||||
.а Рп] (х) ввиду только что сказанного — на (п—k) |
широтах, то |
вся |
||
сфера разбивается на «клетки», внутри |
которых Е(ф, р) |
сохраняет |
по |
стоянный знак. Это значит, что число Я определяет вид формулы коле баний и, в частности, размеры «полуволн» в меридиональном и широт ном направлениях. Чем меньше размеры полуволн, тем больше пара метр Я. В этом случае разница между двумя соседними собственными значениями становится малой по сравнению с их величиной, чем оправ дывается допущение о непрерывности йзменения Я.
§ 14. ФЛАТТЕР ОБОЛОЧЕК И ПЛАСТИН. ПОСТАНОВКА ЗАДАЧИ
Ранее были исследованы автоколебания параметрического харак тера. Автоколебания пластин и оболочек другого типа возникают при взаимодействии их с потоками жидкостей и газов.
Здесь рассматривается весьма интересное и существенное для тех ники больших скоростей явление, называемое панельным флаттером и состоящее в том, что обшивка или другие тонкостенные элементы кон струкции типа пластин и оболочек, обтекаемые сверхзвуковым потоком жидкости или газа, при определенных критических скоростях приходят в колебательное движение с интенсивно нарастающими амплитудами, которое может привести конструкцию к разрушению.
Теоретическое исследование панельного флаттера в корректной, в физико-математическом смысле, постановке стало возможным после того, как в 1947 г. был установлен закон плоских сечений в аэродина мике больших сверхзвуковых скоростей [25].
Анализируя движения тонких твердых тел с большими сверхзвуко выми скоростями в различных средах, А. А. Ильюшин обнаружил сле дующее общее свойство, названное им законом плоских сечений: «Если вектор скорости какой-нибудь точки тела правильной аэродинамической формы есть V и если поперечные скорости других его точек порядка не более &V, то при установившемся и неустановившемся движениях тело вызывает в окружающей среде только поперечные возмущения, причем давление в любой точке поверхности тела, рассчитанное согласно
этому закону, может отличаться |
от |
истинного на |
величину порядка |
|
не больше |
|
|
|
|
1 + /2 = |
_L(V + _L') |
|
||
2М2 |
2 |
V |
М 2 J |
|
|
|
|
Vй |
|
по сравнению с единицей». Здесь М2 = -5—>1 есть число Маха, / ---------- |
||||
|
|
|
Щ |
Vo |
параметр Ильюшина, выражающийся через скорость движения тела V, наклон нормали к площадке е и скорость звука в невозмущенной сре де и0, т. е. скорость звука в газе на бесконечности; этот параметр имеет фундаментальное значение, поскольку в линеаризированной и нели-