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

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

.pdf
Скачиваний:
5
Добавлен:
13.11.2023
Размер:
24.79 Mб
Скачать

Устойчивость оболочек в пределах упругости

191

Исследование задачи в предположении, что ребра, окаймляющие панели, остаются прямолинейными и вместе с тем свободно сближаются между собой, приводит для квадратной панели в первом приближении к следующему уравнению (при V = 0,3):

я* = 7,5^э _ 2 , 0 6 + (0,154Л*2 + 22)

(262)

В случае несмещающихся ребер зависимость имеет вид

Я* = 28,9^3 — 6,1

+ (0,56*2 + 22)

(263)

Для удлиненной цилиндрической панели > Ь) может быть полу­ чено точное решение задачи [2]. Окончательные результаты для верх­ него и нижнего критических давлении для панелей различной кри­

визны н при разных условиях закрепления приведены в табл. 3. Под ^

н я*2 понимают соответственно параметры верхнего и нижнего значений критических нагрузок:

3. Критические нагрузки для удлиненных цилиндрических панелей

Параметр

Шарнирное закрепление

Защемление по длинным

ло длинным кромкам

 

кромкам

КрП Ш 131Ш

 

 

 

 

к

 

 

*

 

 

*1

<7о

ч

 

5

4,29

3.94

 

 

10

19,4

—2, 60

__

20

60,6

—26,5

95,0

47,4

30

95.5

—44,5

40

255

65,2

Коническая панель

Рассмотрим случай пологой конической панели, круговой в плане

при подвергающейся действию поперечной нагрузки интен­

сивностью д (рис. 43, а). Примем, что панель защемлена по контуру

а)

I)

6)

Рис. 43

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

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

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

Устойчивость оболочек в пределах упругости

197

получить лишь после исследования других участков кривой д (0 и определения нижнего критического давления. Окончательные резуль­ таты вычислений с помощью электронных цифровых машин даны на рис. 48. Здесь показаны зависимости между параметром нагрузки ф —

— дЯЧЕ№ н безразмерной стрелой прогиба в центре ^ == ■- ™-ах для

п

панелей различной начальной кривизны, жестко защемленных по кон­

туру. На графике отмечены значения верхнего и нижнего критических давлении.

Максимальное значение -^-на рис. 48 составляет 4, что соответствует

р #

26. Кривые по результатам вычислений, относящихся к панелям не­

 

сколько

большего

подъема (36 <<

 

< р<81),

показаны на рис. 49. Для

 

7

1

1

1

 

 

 

 

 

&

7ой"

 

 

V

дпя по/зной ч

 

 

 

 

 

 

 

 

г ~ ' " 41фС//0/

 

 

А

I

 

*

 

 

 

’1г

 

 

 

 

 

р

I 'II Г

 

 

 

 

X

 

Я

 

 

 

 

 

 

 

 

О

10

70 30 ЬО

50 60

70 80 р

 

Рис. 49

 

 

Рис.

50

 

р =

49 и р = 64 было получено несколько различных ветвей равновес­

ных состояний оболочки, не связанных между собой. Одни из этих вет­ вей соответствуют, по-видимому, устойчивым, а другие — неустойчи­ вым формам равновесия. На рис. 49 по оси ординат отложен параметр

•у = 6 (1 — V2) У З (Г — V2) д*, или

при

V = 0,3 у & 9д*.

(291)

Окончательные

данные для определения

дв и дн при 0 <

р <

80

приведены на рис.

50. Значения г)9 и щ

характеризуют отношения

дв

и дн к величине дов для полной сферы. Кривая для г)в почти совпадает с кривой 1 (рис. 46), полученной другими методами. На рис. 50 нане­ сены также данные экспериментов; ромбы соответствуют нагружению оболочки маслом, остальные — нагружению воздухом. Можно считать, что в первом случае выпучивание происходило при постоянном объеме, а во втором -—при постоянном давлении. Наиболее резкий хлопок на­ блюдался во втором случае; падение нагрузки в опытах показано стрелками. Экспериментальные значения т] лежат в вилке, образуемой г\в и Т1я. В практических расчетах можно пользоваться графиком на рис. 50, для тщательно изготовленных оболочек значения т)н определяют непосредственно по графику, а при относительно большой начальной погиби необходимо уменьшать Т|Л приблизительно на 40%. Для случая шарнирного закрепления оболочки можно получить ориентировочные данные, сравнивая графики рис. 46 и 50.

193

Устойчивость оболочек

ЦИЛИНДРИЧЕСКИЕ

ОБОЛОЧКИ ЗА ПРЕДЕЛАМИ УПРУГОСТИ

Формулой (43) для верхнего критического напряжения рв при осевом сжатии круговой цилиндрической оболочки можно пользоваться при условии, что величина рв не превышает предела пропорциональности

материала: рв

опч. Это условие можно записать в виде

 

 

 

Я ^

1

Е

(292)

 

 

Ь

К 3 ( 1 — V2)

°пц

 

 

 

Для дуралюмина,

при

V = 0,3, алц= 2 * 1 0 3 дан1см2, Е = 7Х

X 10б дан1см2 должно быть

>210. Формула для нижнего критиче­

ского напряжения

 

 

 

 

 

д

 

Рн= о ,т в -^ -

 

 

 

> 70.

 

 

применима, если

 

 

Между тем,

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

в пластической

области.

 

 

 

Ниже приведены результаты решений задач о выпучивании оболочек за пределами упругости, полученные при рассмотрении устойчивости в малом. В практических расчетах следует пользоваться этими реше­ ниями с учетом экспериментальных данных. Эксперименты показывают, что при слабо развитых пластических деформациях необходимо так же, как и в пределах упругости, отличать устойчивость оболочек в малом н в большом. Поэтому рекомендуется при проведении расчетов исполь­ зовать верхние критические значения нагрузок, умноженные на попра­ вочные коэффициенты, учитывающие возможность выпучивания в боль­ шом. Когда пластические деформации значительны, можно вести рас­ четы лишь на устойчивость в малом.

Рассмотрим конкретные случаи. Используемые ниже величины Ес, Ек, фс, Фи определяют по формулам (54), (55) гл. 2; значения этих ве­

личин приведены

в табл. 13 и 14 гл. 2.

ш а р н и р н о

К р у г о в а я

з а м к н у т а я о б о л о ч к а ,

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

ного параметра сжимающей нагрузки Р —

[11:

(293)

где т], О, Я — безразмерные параметры,

 

п2к 9 А_ тпЯ

_ Фк

Цилиндрические оболочки за пределами упругости

199

— число полуволн по образующей; п — число полных волн вдоль

окружности). При осесимметричной форме потери устойчивости будет

Р = Еус

 

 

+

(т+т^Нт*-)’ ]'

(2941

Из условия минимизации полученного выражения по —

опреде­

ляем параметр верхнего критического

напряжения

 

 

 

Р » = -д - КфжФс.

(295)

отсюда

находим

 

 

 

 

 

 

 

Рв-- 4

^

4

*

(295а)

При

фл = Фс =

1 будет

2

_

Н

 

 

 

_

(296)

 

 

Рв~~

Ъ

Е

К '

 

 

 

что вытекает из формулы (43), полученной для упругой области, если принять V = 0,5.

По теории секущего модуля

_ 2 _ Н

(297)

р' ~ 3 Ес К

Сравнение выражений (295а) и (297) показывает, что по теории де­ формаций без учета эффекта разгрузки получается более низкое значе­ ние для рв.

Применение теории течения приводит к формуле

(298) Обратимся к более общему случаю несимметричной формы потери

устойчивости оболочки. Минимизируя

функцию (293) по т),

получим

з ^ х

 

(299)

т| =

 

Подстановка формулы (299) в выражение (293) дает

 

Рв =-3 - ^ 1 Кфкфс

(300)

где

 

 

12(1 — Х)03

(301)

<*» = ] / ч - (1 + ЗХ ) « 4 +

4(ЗХ — 1)03 + 4 •

200

Устойчивость оболочек

Выражение

(293) получает минимум

при

 

2

(302)

 

к г + ж

 

*

Критическое

напряжение при несимметричном выпучивании

Р. = - ~ - Е ^ с У ^ '

 

 

 

 

(303)

Из сравнения выражений (295) и (300) видно,

что при \

ф

1 будет

1, и выпучивание оболочек в малом при наличии пластических

ог

По т е о р и и д е ф о р м а ц и й

 

 

 

 

 

0,8

 

 

 

 

 

0,6

 

 

 

 

 

<№

|

о 5 л а с т $

 

 

1 э к с п е р и ? г е н -™

 

0,2

!1

т а л ь н ы .х д а н /1ЫХ

0

с5 0

1 0 0

1 5 0

2 0 0

Я

 

 

Рис.

62

 

 

деформаций должно происходить по осесимметричной форме, следова­ тельно, при проведении практических расчетов следует пользоваться формулой (295а).

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

к радиусу = 2. Значение изменялось в пределах 20—130. По-

теря устойчивости при слабо развитых пластических деформациях со­ провождалась хлопком; при этом вмятины имели ромбовидное очерта­

ние, как в упругой области. С уменьшением число волн по окруж­

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

На рис. 51 сплошная линия отвечает теории деформаций, штрихпунктирная — теории течения. Теоретические кривые мало отличаются одна от другой. Штриховая линия соответствует решению в упругой области (V = 0,5). Область экспериментальных точек ограничена тон­ кими линиями.

Результаты экспериментов в виде зависимости коэффициента а =

= -Е*Р. от величины

—- показаны на рис. 52; рв вычислено по теории

Рв

Л