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

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

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

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

171

Исключая из этих уравнений функцию %, получим следующее раз­

решающее уравнение:

(-188>

4И*3^ М 4 -^>=°-

Уравнение (188) используют в качестве исходного.

Для случая шарнирно опертой оболочки примем функцию ф в виде

где

ф = Ах3 (1 + ах) 51П пх,

(189)

 

б +

V

 

 

 

а

 

(190)

 

8-1- V ’

 

 

 

 

 

при V =

0,3 а = —0,759.

 

 

 

 

Уравнение Бубнова-Галеркина, соответствующее дифференциаль­

ному уравнению (188), имеет вид

 

 

 

1

(* з 4 4 4 ^Л х + р * \

*

1 у * а х = °-

(,э|)

О

 

0

0

 

Отсюда находим выражение для ц:

 

 

 

\1 =

 

1

 

 

(192)

 

 

^ ф2 с1х

 

 

 

 

 

6

 

 

 

Для оболочки, жестко защемленной по краю, будет

= 0,

тогда

 

 

 

 

 

 

 

 

 

йх

 

 

ц = -

 

 

 

(193)

I**с1х

Выражение (189) подставляем в формулу (192). Из условия мини­ мума р, по р находим рщщ = 20,4; используя формулы (181), (182), (185), получим значение верхнего критического давления для оболочки, замкнутой в вершине:

’ « = 3 , Ш : ( А ) \ е а | / ГА

а

 

V (I — V*)»

 

 

 

 

 

 

•5

_з_

1

 

 

 

(194)

= 3,1ЪЕ

 

а) 2

 

_з_ *

(1 - V *)4

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

Решение задачи

методом конечных разностей приводит к следу­

ющей формуле для

критического давления

[11:

 

5

з

 

 

(195)

(1 — V2) 4

Усеченная коническая оболочка (рис. 28). Докритические усилия определяют по формулам

(196)

где 10 — расстояние по образующей от вершины до меньшего основа-

ния. Вводя безразмерную

координату х =

, приводим эти выраже-

ния к виду

 

ч

 

 

 

 

 

 

ре3ЕЛ

а

; Д/2 =

— дгр,е3ЕЛ

а. (197)

Величины е и р. определяют по формулам (181) и (183). Как и в слу­ чае оболочки, замкнутой в вершине, разрешающее уравнение после упрощений получается в виде уравнения (188). Но если в случае замкну­

Рис. 28

той оболочки величина х изменялась в пределах 0 ^ х «$:• 1, то теперь

она заключена в пределах ^ х ^ 1.

В результате решения задачи по методу конечных разностей для слу­

чая,

когда усеченная оболочка по

меньшему основанию защемлена»

а по большему — шарнирно

оперта,

получается следующая формула

для

верхнего критического

давления

[1]:

 

 

5

3

(198)

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

173

Коэффициент Сх берут по графику на рис. 29.

Формулы (194), (195), (198) применимы, если угол а не близок ни

лг

кнулю, ни к - у [в этом предположении получены исходные уравнения

(186)—(188)]. Как показывают эксперименты, этими формулами можно

пользоваться при значениях угла а, лежащих в пределах 20°

^ 80°;

если а^> 80°, то при

вычислении

ув кони-

 

ческую оболочку можно заменить цилиндри­

 

ческой с радиусом, равным среднему радиусу

 

конической

оболочки.

 

 

 

 

 

 

Испытания тщательно изготовленных зам­

 

кнутых и усеченных конических оболочек

 

при внешнем давлении показывают, что потеря

 

устойчивости сопровождается

хлопком. Для

 

тщательно

изготовленных

оболочек экспери­

 

ментальное

значение

критического давления

 

получается

несколько

выше

теоретического,

 

определенного по формуле

(195)

для замк­

 

нутых в вершине оболочек и по формуле (198)

 

для усеченных оболочек,

и ниже

значения,

 

получаемого по формуле

(194).

В

среднем

 

экспериментальная величина

критического

 

давления приблизительно

на

5%

выше зна­

 

чения, вычисленного по формуле (195). и

 

примерно

настолько

же

ниже

величины,

 

определенной по формуле (194).

 

 

 

Впрактических расчетах можно использовать формулы (195) и

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

Кручение

Приведем результаты решения задачи об устойчивости в малом усе­ ченной, свободно опертой конической оболочки, нагруженной равно­ мерно распределенными по торцовым сечениям сдвигающими усилиями. Координату 5 по-прежнему отсчитываем по образующей от вершины конуса; соответствующие расстояния от вершины до меньшего и боль­ шего оснований обозначим через /0 и /х, длину оболочки по образую­ щей — через /. Сдвигающее усилие на единицу длины контура в пло­

скости меньшего основания (при 5 =

/0) обозначим через

а в пло­

скости большего основания — через

Т

эти усилия выражаются че­

рез крутящий момент Мк по формулам

 

 

Т

______ МК

г р ______ МК

,.дд

2л (/о соз а)2»

11

2л (/х со$ а )21

 

где а — угол наклона образующей

к основанию.

 

Приближенная формула для критического сдвигающего усилия имеет вид [1]

ДМРа- 0 Е -~

С0 *б4 о.

(200)

Т1о = 4 (я г + $ а)

 

 

174

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

 

 

 

лк

 

Г я а + €г

1 Т

р = 1 + т г ; е4 = х -

 

[ б ( р * - Ш

(201)

 

С 2

 

 

 

Е = 1п ( 1+ т ) :

 

 

 

=2,61 ( ^ ) т

[ н - ® д а ( - 5 - ) т

« н - п > ] ;

 

 

1

 

 

(202)

 

 

 

 

 

 

4

2(1 +

*) *

 

 

Условия

применимости

приведенных формул:

 

 

е г ^ 0,1; 0 ,2 55^*^0,33 .

 

Ниже приведены некоторые значения коэффициента Са при V = 0,3:

е

0,12

0,10

0,07

0,ОБ

0,03

Со

3,75

3,87

3,95

4,08

4,31

Подкрепленные оболочки при действии внешнего давления

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

следующие обозначения:

 

 

 

24 (1

V2) у,,.

 

 

,2 = 1 + 2 4 Я

- у а )уР;

(203)

 

 

«>1=Ч

2&,6?

(1 — V»);

 

 

- Л

 

 

2 Ь ^

-(1—V2),

(204)

 

 

<йа= 1-

ин

где

и

— моменты инерции половины сечения ребра относительно

осей, лежащих в срединной поверхности,

, Р_

Ч * ? )3

+Т^в? (*+«?)*!

^ -------- \ 2 ~

_

? ) 3

 

12

+ Т « ( И б ? ) г;

 

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

упругости

175

Л — толщина оболочки;

Ь2 — значения

толщины

продольных

(рис. 30, а) и круговых (рис. 30, б) ребер соответственно; б{\ 6$ — вы­

сота ребер; 1Ъ (2 — шаг ребер.

Если оболочка подкреплена односторонними ребрами (сечение кру­ говых ребер показано на рис. 31), то формулы (203) и (204) примут вид

(1 = 1

, 12 (1 — V2)

гРт

<2=1

+

12(^

уа)^ - .

(205)

 

 

 

6,6?

 

а г =

1+

6, 6?

(206)

- ^

- ( 1- ^ ) .

здесь У{* — момент инерции сечения продольного ребра относительно

центральной оси профиля кругового сечения оболочки; У? — момент

инерции сечения кругового ребра относительно центральной оси про­ филя осевого сечения оболочки.

Верхнее критическое значение внешнего давления для шарнирно

опертой замкнутой

круговой

конической оболочки будет

[1]

 

 

 

5

з

 

 

л

и

; : ) 2

(*8 а ) 2з >

(207)

 

 

 

 

 

 

( ( - V 2) 4

 

 

Р, =

 

— V2

-

(208)

 

 

 

 

 

_

 

(1 V2) (02

 

(209)

 

На

 

V 2

 

 

 

 

 

Формула (207)

получена

в предположении, что угол а

не близок

ни к нулю, ни к

В результате решения задачи с помощью метода

конечных разностей получается формула (207) с коэффициентом 2,8 вместо 3,15.

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

Для усеченной оболочки верхнее критическое давление

(1 — V2) 4

Значения коэффициентов Съможно определять по графику на рис. 29 при тех же граничных условиях.

176

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

ОБОЛОЧКИ ВРАЩЕНИЯ

Сферические оболочки

Во многих областях техники применяют сферические оболочки

в виде сегментов, закрепленных по краю. Следует различать сегменты

сбольшим углом охвата, стрела подъема Н которых сравнима с радиу­ сом кривизны срединной поверхности К (рис. 32), и пологие панели, для которых Н < Будем рассматривать устойчивость сегментов большого подъема. При потере устойчивости такого сегмента образуются сравнительно мелкие вмятины, и критические нагрузки для него будут темн же, что и для полной сферической оболочки. Поэтому будем опре­ делять критические нагрузки для полной сферической оболочки, на­

ходящейся под действием внешнего равномерного давления 7.

Примем, что в пределах первич­ ной вмятины оболочка является по­ логой. Исходные уравнения линей­ ной задачи имеют вид [1]

 

° ’ 4“" =

У2<1>+

Чг,

 

 

 

 

1

,

 

I

( 211)

 

 

 

 

 

где ф — функция усилий в срединной

поверхности;

К — радиус

 

сре-

динной поверхности; Н— толщина оболочки;

О =

 

Е№

 

 

у 3

—— —----- у—;

 

и у 4 — операторы по формулам (11)

и

(13).

 

12(IV3) •

 

Исключая функцию <р, получим

 

 

 

 

 

 

 

 

Е1г

(212)

Найдем величину дг. Если внешнее давление <7, то нормальные усилия во всех нормальных сечениях оболочки

Сжимающие усилия N и напряжения о считают положительными, тогда

2 »

2Л ’

(213)

 

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

ция укорочения вдоль дуги любого нормального сечения будет е =

с другой стороны, эта величина

 

ст(1— у)

? /? (! — у ) .

Е

2ЕН »

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

177

отсюда находим

(214)

В дальнейшем через до будем обозначать дополнительный прогиб, имеющий место при выпучивании оболочки. Величина в уравне­ ниях (211), (212): дг= —/кту2до. Уравнение (212) принимает вид

О2

V*® + У2'*' = 0. (215)

Примем, что решение уравнения (215) должно удовлетворять соот­

ношению

 

у2до =

— Я,2до,

(216)

 

 

где Я — неопределенный

параметр.

 

 

 

Тогда из уравнения

(215) получим

 

 

о =

^ 12

г

Е

(217)

 

 

/I

+

/?2Я2 "

 

Минимизируя о по Я2,

находим

 

 

V

 

1

 

 

(218)

= - V 12(1 -V * ).

Далее, по формуле (217) определяем величину верхнего критиче­

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

 

 

 

 

1

Е 4 -

« 0 .6 0 5 ^ 4 - .

(219)

сгв

1/3(1 - V

2)

К

 

Напряжению по формуле

(219)

соответствует давление

 

9

Р

 

 

( 220)

Я*

 

 

К 3 ( 1 — V2)

 

 

 

 

Формула (219) для верхнего критического напряжения для сфери­ ческой оболочки в точности совпадает с выражением (43) для верхнего критического напряжения, отвечающим случаю осевого сжатия цилин­ дрической оболочки.

Определяя верхнюю критическую нагрузку, мы не рассматривали форму выпучивания оболочки. Исследование формы изогнутой поверх­ ности, проведенное в линейной постановке, показывает, что в случае, если выпучивание является осесимметричным, потеря устойчивости должна сопровождаться появлением одной сравнительно глубокой вмя­ тины и ряда кольцевых складок, размельчающихся по мере удаления от центра основной вмятины [1].

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

178

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

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

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

(221)

— у 4Ф -------- — Ь (ш,

где Ф — функция напряжений в срединной поверхности; <7— попереч­ ная нагрузка; Ь — оператор по формулам (32) и (36).

Воспользуемся полярными координатами г, ср. Совместим начало радиуса-вектора с центром вмятины. Если функции ш и Ф не зависят

от ф, то получаем

следующие выражения для операторов, входящих

в уравнения (221)

[1]:

 

_ 1 _ д _ / й т _ Л Ф \ #

 

= ~ 7 ' йг \ йг * йг ) *

 

(222)

Уравнения

(221)

примут вид

 

 

_

1

Л ( 4Ф\ йЛ /( йФА

йюш \ ,

дг

 

Ж'Чг\г Лг)+ йЛ йг '

/

А ’

(223)

 

 

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

 

179

Интегрируя эти уравнения по г,

получим после деления на г

 

й

й . 0 ,

1

 

йФ_

йш

 

|__ у .

 

 

 

 

 

 

Я

 

йг

йг

 

йг

+

Л

 

(224)

 

 

й

. в

 

Г

1

йт

[

/

йш

\2."1

 

 

 

 

 

 

 

 

~ й г^ ^ ~~

[

/? ’

йг ~т"

(

йг

)

] 9

 

 

под 4е

понимают функцию нагрузки,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

г

 

 

 

 

 

 

 

 

 

 

 

V ^ - у ^ Ч г й г .

 

 

 

 

 

 

(225)

 

 

 

 

 

 

г

О

 

 

 

 

 

 

 

Напряжения

в срединной

 

поверхности

в общем

случае

[I]

-

г

дг

о- _!_

&2ф .

 

_

д*ф .

 

_

д

/

\

дФ \

°г ~

г3

дф2 *

 

С(р ~~

дг2 *

Т ~

 

дг

\

г

дф ) '

В осесимметричных задачах бу-

. ____

 

 

 

с

 

 

_____ ,

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

 

_

1

( й2Ф _

V

йФ \

 

 

 

 

 

Т ’Ч Г ) '

 

 

и

 

находим

 

С другой стороны, еф =

-р-, отсюда

 

 

 

г

( й2Ф ___ йФ \

 

(227)

 

и ~~

Е

\ йг2

г

йг

/

 

 

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

ш = 0;

= 0 при г = с,

(228)

где с — радиус вмятины. Предполагалось, далее, что по краю вмятины отсутствуют радиальные перемещения:

и = 0 при г = с.

(229)

180

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

 

 

С учетом выражения (227) эго условие можно записать в виде

^2Ф

у

0 при г =

с.

(230)

йг*

г

57"

 

 

 

В полюсе при г -> 0 радиальные напряжения должны быть огра­ ничены по величине; отсюда из формул (226) вытекает четвертое гранич­

ное условие

 

ЛФ

 

 

 

 

 

0 при г = 0 .

 

(231)

 

 

Т— =

 

Выражение для

прогиба аппроксимируется в виде

 

 

 

 

 

 

(232)

'

 

___ —• л

Эта функция удовлетворяет гранич-

0,2 -

 

ным условиям

(228). Выражение (232)

-----— *

 

 

подставляют в

правую

часть второго

10

15

20

из уравнений (224). Интегрирование

результата подстановки

ведут с учетом

Р и с.

34

 

условий (230)

и (231),

Далее исполь­

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

дят зависимость а (&),

показанную

на рис. 34 [1].

Под ^ и с пони­

маются безразмерные величины:

 

 

 

 

 

 

 

 

ЯН2

(233)

ь

А

°

ЕН

2ЕН2

 

Минимальное значение

напряжения

 

 

ан = 0,155,

или

о* = 0,155Е-^-

(234)

при стреле прогиба ^ =

9,16. Нижнему критическому напряжению по

(234) соответствует давление

 

 

 

 

 

? « = 0 ,3 1 *

( - ^

) 2

(235)

Полученное значение ан почти в 4 раза меньше верхнего критического напряжения ств.

В результате другого, уточненного решения задачи по методу Ритца нижнее критическое напряжение оказалось

0,193

0,203.

Граничные условия, принятые при получении формулы (234), не вполне соответствуют истинному характеру сопряжения зоны вмя­ тины с остальной частью оболочки: защемление по контуру вмятины надо считать не жестким, а упругим. Решение задачи в таком предполо­ жении по методу Бубнова-Галеркина в первом приближении приводит