192 |
|
Устойчивость оболочек |
Представляя |
оболочку как искривленную пластинку, возьмем |
|
за основу уравнения (250)—(251). По рис. 43, б находим г = —— = аг, |
||
где а « |
а = |
и |
Выразим уравнения (250)—(251) в полярных коор |
||
динатах г, ф (рис. 43, в).
Производные от некоторой функции Ъ по я, у выражаются через
производные от 2 по г, ф в виде |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
дЪ__д1_, |
дЪ |
|
1 |
д! |
|
|
|
|
|
(264) |
|||||||
|
|
|
|
|
|
дх ~~ дг * |
ду |
|
г |
дф 1 |
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
дП |
_ |
|
д |
( \ |
|
|
дЪ \ > |
|
|
|
|||
|
|
|
|
|
дх1 ~~ дг2 1 |
|
дх ду |
|
дг |
\ |
г |
дц> ) |
7 |
|
|
|
||||||
|
|
|
|
|
|
д-1 |
|
_ 1_ а2 |
|
_1_ ^ 2 |
|
|
|
|
|
|
(265) |
|||||
|
|
|
|
|
|
ду2 |
г |
’ |
дг |
' |
г2 |
дфа * |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
Оператор |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
07 _ |
а°-2 |
|
ду- |
~ |
дг* |
+ |
г |
' |
дг |
|
|
г- |
|
аа2 |
|
|
|
||
|
|
|
|
- |
Зл:3 |
+ |
|
|
|
З ф - |
|
|
|
|||||||||
|
|
|
|
|
_ |
1 |
|
а |
/ |
а г \ |
, |
|
|
|
з^г |
|
|
|
|
(266) |
||
|
|
|
|
|
- |
г |
’ аг |
V |
* • ) + |
г'- ' а<ра • |
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||||
Используя |
зависимости |
(264) — (266) и вводя относительную |
|
коор - |
||||||||||||||||||
днпату |
р = |
|
приводим основные |
уравнения |
(250) — (251) |
к |
виду |
|||||||||||||||
в |
4 |
|
_ |
1 |
а |
/а® |
дм_\ . |
1 |
( |
а2Ф |
|
а^ |
|
а»Ф_ |
а2ш \ |
|||||||
л |
у |
ш - |
р |
' ар |
\ аР ‘ а р ) + |
р М |
зф2 |
‘ |
ар3 + |
зРз |
' зфз |
+ |
||||||||||
|
|
2 / |
а2Ф |
а^ |
|
аФ |
а2ву |
\ |
|
2 |
|
|
а2Ф |
_а^о_ |
|
|
||||||
|
" "рз" \ |
дрдф" ‘ аф |
“аф"4 ар аФ / |
р2 |
* ар аф ’ ар аф |
|
|
|||||||||||||||
|
|
|
|
|
2 |
аф |
дП) . |
с а |
|
Vа2ФVI/ |
|
|
д |
4. |
|
|
(267) |
|||||
|
|
|
|
|
р4 |
аф |
|
аф |
+ |
- |
|
ар2 |
+ |
л |
|
’ |
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
1 |
|
|
^ ( д2® \2 |
|
2 |
а^ |
|
а2ш |
|
|
|
||||||||
|
|
|
"2Г ^ |
= "р2" \ араф ) |
|
рз*'“аф |
араф |
|
|
|
||||||||||||
|
_1_ / |
д-ш |
) |
___ 1_ |
дни |
|
дгш ___ 1_ _ а2ад |
|
^ а2ср |
_ |
са_ д2ш _ |
|
||||||||||
■ |
"р4^\ |
аф |
р |
ар- * ар2 |
р2 |
ар2 |
|
|
аф2 |
|
р |
ар2 |
’ |
* |
||||||||
оператор |
Лапласа имеет вид |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
2 _ |
а2 |
|
|
1 |
а |
|
1 |
|
а2 |
|
|
|
|
|||
|
|
|
|
|
^ |
ар2 "• |
р |
’ ар |
*“ |
р2 |
|
аф2 |
* |
|
|
|
||||||
В соответствии с принятыми граничными условиями, имеем |
|
|
ш = 0; |
—0; ог —0; т = 0 при р = 1. |
(269) |
Устойчивость оболочек в пределах упругости |
193 |
Аппроксимируем прогиб выражением
ы) = Д (1 — р2)2 + /гр4 (1 — Р2)2 соз л<р. |
(270) |
Первый член в формуле (270) отвечает решению осесимметричной линейной задачи для защемленной пластинки, а второй — отображает изгиб по п волнам вдоль дуги с равными амплитудами к центру кри визны; эти два члена во взаимном сочетании характеризуют образова ние ряда вмятин, преимущественно направленных к центру кривизны.
В результате подстановки выражения (270) в правую часть уравне ния (268) находим
- I - у 4 ф = |
Я 0 + К п С О З Л ф - I - |
|
+ |
соз 2лф, |
(271) |
где # 0, /?Л — функции р, содер жащие параметры /ь / 2- Интег рирование равенства (271) про водим с учетом того, что в центре панели при р = 0 величина сгг должна быть ограниченной; при этом приходим к следующей за висимости:
А - ф = ф 0 -\- ф п соз Лф -|-
+ Ф-2п соз 2лф, |
(272) |
где Ф0, Фл, Ф2П — новые функции р, / ь / 2. Интегрируя, далее, уравне ние (267) по методу Бубнова-Галер кина, получим
7* = |
€&• + |
+ |
Сз^ + |
С ^ + |
СЙБ: |
(273) |
ВД* + |
В& + |
В Ц + |
ВД* + |
Бвр + |
5 в й а = 0; |
(274) |
здесь введены параметры
Коэффициенты С*, . . С6 и В1г . . Вв зависят от числа волн п. Исключая | из системы (273)—(274), определяем зависимость параметра поперечной нагрузки <7* от безразмерной величины Ё, характеризующей прогиб центра панели. Случай п = 0 соответствует осесимметричной форме прогиба.
Зависимость <7* (0 для -|г- = 100 при а = 0,1 и /1 = 0 (сплошная
линия) приведена на рис. 44. Кривые для п Ф 0 до 5 = 7,5 распола гаются выше кривой, отвечающей осесимметричной форме равновесия, а при I, > 7,5 — ниже этой кривой. На рис. 44 для примера приведена штриховая линия, соответствующая несимметричной форме выпучива ния при п = 10. Следовательно, несимметричная форма должна иметь место при прогибе в центре, равном 7,5 толщинам оболочки. В
приведенном примере Гпри-“ = 100^ огибающая семейства кривых,
194 |
Устойчивость |
оболочек |
соответствующих различным л, при К,> |
7,5 мало отличается от кривой, |
|
соответствующей осесимметричной форме. Нижнее критическое давле
ние оказывается близким к нулю. |
|
|
При |
жесткой заделке панели, когда точки кромки не смещаются, |
|
нижнее |
критическое давление при |
= 100 принимает отрицательное |
значение.
Исследование осесимметричного выпучивания конической панели со скользящей заделкой по контуру, испытывающей действие сосредо точенной в центре силы Р, направленной вдоль оси, приводит к сле дующей зависимости между нагрузкой Р и стрелой прогиба /:
Я* = 0,294^ — 0,332- у - г;2 + |
[о ,16 |
+ 1.47^ |
(276) |
|
где |
Р<? |
|
|
|
Я* = |
й = - |
|
|
|
|
ЯАа |
|
|
|
Сферическая панель
Рассмотрим устойчивость пологого сферического сегмента, нагружен ного равномерно распределенным внешним давлением ? (рис. 45).
Величины //, |
Р, |
с связаны соотношениями |
|
|
|
|
|
||
|
|
|
|
с* * 2РЯ; |
|
|
|
||
|
|
|
|
|
р р а |
Н |
|
(277) |
|
|
|
|
|
|
н |
* |
2 Т - |
|
|
|
|
|
|
|
|
|
|||
|
|
|
Решим |
осесимметричную |
задачу |
||||
|
|
|
о выпучивании |
пологой |
сферической |
||||
|
|
|
оболочки. В уравнениях (248а)—(249а) |
||||||
|
|
|
примем |
доо= 0 . Вводя полярные |
коор |
||||
Рис. |
45 |
динаты, |
считая |
до = |
ш (г), |
Ф = |
Ф (г), |
||
ний. При |
д = |
сопз! |
определяем первый интеграл этих урав- |
||||||
из уравнений (224) |
получим |
|
|
||||||
|
В |
а |
1 |
4 |
*/до |
/ а® |
' |
\ |
4Ф |
, |
(/Г2 |
|
Н |
йг |
г |
»Й |
1 V. |ТЗ |
\ “ 5 Г + |
Я |
) |
"й Г + |
2Л |
|
|
1 |
а |
1 |
а |
йФ |
|
йдо |
. г |
\ |
(/до |
|
|
Е |
йг |
г |
йг Г Аг ” |
( 2 ‘ |
|
|
+ Т Г ] йг |
|||
|
Вводя |
безразмерные |
параметры |
|
|
|
|
|
|
||
|
|
|
|
|
I |
. . |
|
|
1 |
М> |
|
|
|
|
|
Яр3 ' <Ш + в ; у * |
Е1ф }/"т' |
|
’ |
||||
Р |
1 |
|
|
Я = |
|
|
|
|
|
1 |
|
К т ‘ Л ’ |
|
|
|
т |
12 (1 — V3 ) |
||||||
|
|
|
2 Я |Л п |
|
|
||||||
(278)
(279)
(280)
(281)
Устойчивость оболочек в пределах упругости |
155 |
||||||||
приводим уравнения |
(278)—(279) |
к виду |
|
|
|
||||
|
и |
(а) = |
р ' (ау + |
Р 02); |
(282) |
||||
|
|
(V) = |
|
Р' (О3— а4), |
(283) |
||||
где С — линейный оператор, |
Л |
1 |
1 |
|
|
||||
|
1* _ |
0 — |
, _ . |
|
- |
0 |
(284) |
||
|
Ь ~ |
|
40 |
0 |
М |
|
|
||
Введем параметр к -- |
№ |
|
и величину |
а по формуле (233); |
вели |
||||
чины р' и Р связаны с |
этими |
параметрами |
соотношениями |
|
|||||
|
|
|
|
|
|
|
а |
|
(285) |
|
Р' |
Ут |
к' |
V т |
|||||
|
|
||||||||
Если панель защемлена по контуру и радиальные смещения точек |
|||||||||
контура отсутствуют, |
граничные |
условия |
получают вид |
|
|||||
|
а = |
О, |
у = |
0 при |
0 = |
0; |
(286) |
||
а = |
1, |
|
---- Vу = 0 |
при |
0 = 1. |
(287) |
|||
Ниже приведены окончательные результаты решений, полученных с помощью электронных цифровых машин. Данные, относящиеся копределению верхней критической нагрузки ?в». показаны на рис. 46. Кри
вая 1 соответствует жесткому защемлению по контуру» кривая 2 — шар нирному закреплению с несмещающдмися в плане точками контура, кривая 3 — случаю шарнирного, опирання со свободно смещающимся в плане контуром. По оси абсцисс отложена величина
р = |
Г2р' = |
V 12 (I— |
V2) |
= (^12(1 - V 2) |
* . |
(288) |
|
Используя |
условия |
(269), |
параметр р выразим в виде |
|
|
||
р = 2 ^ 1 2 (1 — V2) |
; |
при |
V = 0,3 будет р р*6,6 |
(289) |
|||
196 |
Устойчивость |
оболочек. |
|
||
По оси ординат отложена |
величина |
|
|||
|
К 3 ( 1 - Г г) Чв ( К V |
(290) |
|||
Ч » - |
2 |
' |
Е \ Л ) ' |
||
|
|||||
здесь <7* — значение верхнего критического давления для полной сферы, определяемое по формуле (220); обозначим его через д0,в. Тогда будет
|
а) |
6) |
в) |
|
|
Рис. 47 |
|
т|в = — |
Процесс |
вычислений значений г), |
соответствующих участ |
ие. в кам кривых, показанных на рис. 46 штриховой линией, сходился
настолько медленно, что относящиеся к указанным участкам значения нельзя считать оконча
тельными. |
защемленной |
|||
Для |
|
|||
панели |
(кривая 1): |
|
||
р < |
25 |
|
п * « 1 ; |
|
35 < р < 75 |
т|Л> |
1; |
||
р > 7 5 |
|
1]в < |
1. |
|
Такой |
вид |
кривой |
||
т)в (р) объясняется раз |
||||
личным |
|
характером |
||
волнообразования обо |
||||
лочки |
в |
зависимости |
||
от р. Формы |
волнооб |
|||
разования сферических панелей различной кри визны по данным экспе
риментов, проведенных над оболочками, изготовленными из магниевого сплава и жестко защемленными по контуру, показаны на рис. 47. Про гиб, максимальный в центре и монотонно уменьшающийся по направ лению к краю (рис. 47, а), имел место при малых р. При р > 20 макси мальный прогиб не в центре (рис. 47, б); в случае развитой вмятины максимум прогиба будет уже на расстоянии половины радиуса от центра. При р ;> 55 появлялся новый пик в центре (рис. 47, в).
При 15<Зр<*30 величина <7« при шарнирном закреплении (кри вая 2 рис. 46) лежит выше, чем при защемлении панели по контуру. Это объясняется, возможно, различным влиянием изгибных напряже ний, образующихся у края. В остальной области изменения р случаю шарнирно закрепленной оболочки соответствуют значения г|*. меньшие, чем для защемленной. Сильно падает значение Т1в при свободно смеща ющемся в плане контуре (кривая 3, рис. 46).
Выпучивание реальных оболочек происходит при нагрузке, лежащей, как правило, значительно ниже дв. Поэтому данные для расчета можно