Материал: Прочность, устойчивость, колебания. Т. 3

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

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

*/до

/ а®

'

\

,

(/Г2

 

Н

йг

г

»Й

1 V. |ТЗ

\ “ 5 Г +

Я

)

"й Г +

 

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).

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

Источник: https://tut-files.ru/previewfile/74208