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

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

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

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

131

где

О =

ЕН*

12(1 ^ 2)*

Из уравнений (7) определяют поперечные силы (?* и <3У:

(Ю )

здесь у 2 — оператор Лапласа,

 

 

 

 

а3 .

 

а*

 

 

(И)

 

 

 

у-==

 

 

 

 

 

Подстановка выражений (10) в уравнение равновесия (б)

приводит

к следующей

зависимости:

 

 

 

 

 

 

 

 

 

 

а

а

о»

.

^

 

( 12)

 

 

 

- г * ш= т г + т - ’

 

 

 

 

 

 

где у 4 — двойной оператор

Лапласа,

 

 

 

а4

 

 

У4 = 7 2Уа =

- | г + 2

 

а4

(13)

 

 

дх*ду* +

ду*

Исходя из соотношений (8), выразим деформации через напряжен

 

 

 

= ~]7 (ст.г — \Оу)]

 

 

 

 

&у = - ^ { а у — \ох)\

 

(14)

 

 

 

 

 

 

 

 

 

 

 

 

 

У = 2т (1 +

 

у)

 

 

 

ций (4), получим

 

 

 

 

 

 

 

 

 

_1_ Г &а*

2 а ч

. +

 

V

(

дх*

 

 

Е 1 ду*

дхду

^

дх2

 

\

 

 

 

 

 

 

 

1

а2м>

 

 

(15)

 

 

 

 

 

/? ' дх2 *

 

 

 

 

 

 

 

 

Если ввести функцию напряжений в

срединной поверхности Ф

по формулам

а2Ф

 

 

а*Ф

 

 

а2Ф

 

 

 

 

 

 

 

° * = - 3 7 г ;

аи = ~ № ’ х = ~~я7д77>

(16>

ду*

дхду •

 

132

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

 

то зависимости (12) и (15) можно представить в виде

 

 

 

д*Ф

я

07)

Н

 

дх2 + Л ;

 

у4Ф =

1

д2ьи

(18)

 

 

 

При решении задач устойчивости в уравнение (17) нужно подставить

вместо ц фиктивную поперечную нагрузку <Д равную сумме дополни­ тельных проекций основных усилий рх»/?«/, 5 на направление нормали (усилие рх действует вдоль оси х, усилие ру — вдоль касательной к ли­ нии (/, усилия 5 — касательные):

Положительными считают усилия, способствующие увеличению параметров кривизны; в отношении рх и ру положительными считают усилия сжатия. Подставляя выражение (19) в формулу (17), получим

О .

 

1

д2Ф

 

даш

 

 

д2&

п д2и>

(

20)

А

^ 10

~

7?

дх2

 

^х дх2

 

 

ду1

8 дхду

 

 

 

 

 

Применим к

уравнению

(20) оператор

 

V 4, а

к уравнению (18) —

оператор

 

 

тогда выражения (20) и (18) приводятся к одному

разрешающему

уравнению

 

 

 

 

 

 

 

 

 

й

,

Е

д*м

,

_ _4 /

)

\

+

 

4 / а-ш \

 

 

А - у '»

+

 

" з р - +

Ч ^

)

РЛ *[: щ Г ) +

 

 

 

\ - ш

 

+

РЛ*

 

 

+МЙ-)-«: 0. (21)

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

гают

Мх = Н =

0; учитывают

только усилия

0 у и

Му. При

таких

условиях справедливы

следующие

соотношения:

 

_

ди

____/ д2ад

 

щ \

 

д о ___ _

ди_ ___ до

х ~

дх 1

ху

\ ду2

Я2 /

1

ду ~~ К '

ду

дх

Отсюда

 

д2ел

____ 1_

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ду2

~

 

У?

дх2

 

 

Уравнение

совместности

деформаций получает вид

 

 

 

 

о

д*ех

 

1

дЧх

_

д2ку

 

(22)

 

 

к

ду*

 

/? '

ду2

~

дх2 "

 

 

 

 

 

 

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

133

Приходящуюся на единицу площади внешнюю нагрузку вдо#ь оси х, касательной к линии у и оси г, обозначим через Яу> Яг. Тогда урав­ нения равновесия запишутся

д а ,

д т

 

 

1

 

д т

, дОу_____1

дМу .

Яу

0 .

-д Г + ~д? +

^

= 0-

 

дх

 

ду

 

 

ду

К

'

 

 

 

 

дЧ1и

 

 

 

 

 

 

 

 

 

 

 

 

 

ду2

+ х

0*+?*=а

 

 

 

Объединяя их,

получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

Ш и

,

*<*,

 

(23)

 

*

ду4

+

Я

* д</2

+

 

дх*

~ '

 

где

 

 

 

 

 

 

_

 

 

,

 

 

 

, А *

 

 

 

 

 

 

р

 

дЯх

 

 

 

 

 

(24)

 

 

 

^

"

" д Г +

д*,

*

ду*

 

 

 

 

 

 

Соотношения

закона

Гука

имеют

вид

 

 

 

 

 

 

 

Е

 

 

 

 

 

 

 

Я/13

 

 

 

Ох = 1 — Vт с, & Еех

М ;,= -

12(1— V2)

 

(25)

тч

 

 

 

 

 

Х

О

 

У

 

 

 

 

 

Вводя

переменные а = —

, р =

- - и пользуясь оператором

 

 

 

 

 

 

 

д4

,

д2

 

 

 

 

(26)

 

 

 

 

 

й

= - ар4

+

ар3

 

 

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

и (23) приводим к виду

 

 

 

 

 

й<т,

12/?

 

 

_

 

 

 

 

 

 

 

(27)

Л3

*

да2

"

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дЦх .

д<7у

дгдг

 

 

(28)

 

 

 

 

 

 

да +

ар

ар2 *

 

 

 

 

 

 

 

 

 

 

 

Исключая ах и используя выражение (25) для Му, получаем, после

исключения оператора + 1 )• следующее дифференциальное урав­

нение для решения задач устойчивости:

(29)

Уравнением (29) следует пользоваться при исследовании устойчи­ вости оболочек средней и особенно большой длины в случае слабо вы­ раженного волнообразования по длине оболочки.

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

134

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

Дополняя соотношения (2) нелинейными членами, получим сле­ дующие выражения для деформаций в срединной поверхности:

 

 

_

ди

до

дао

дао

 

 

 

 

^~~

ду

' дх

дх

 

ду

 

 

Уравнение совместности деформаций принимает вид

дЧх

, дЧу

 

д2у

 

1

 

 

1

д2ао

ду2

дх2

дх ду

 

~ г

М в *

 

(31)

 

 

 

где 1 (и;, ао) — оператор,

 

 

 

 

 

 

 

 

С {ао, но) — 2 Г д2ао

д2ао

/

д2ю

\ 2~|

(32)

 

 

 

[

дх2 *

ду1

\

дхду

) ] ‘

 

Изменения кривизн и кручение срединной поверхности определяют по формулам (3).

Соотношения Гука (8) остаются прежними. Первые два уравнения равновесия отвечают соотношениям (5), третье уравнение равновесия

принимает форму

 

 

 

 

 

дОх ,

д()у .

, д2ао

{ 1

д2ао \

 

~дГ +

+

Ь ё +

а* \~ К

+ ~ д ^ ) +

 

 

+

2 т й - Ц

г +

‘' = ° -

(33)

 

 

Для поперечных сил справедливы соотношения (10). Подставим эти выражения в уравнение (33), а соотношения (14) — в уравнение (31),

тогда, вводя функцию напряжений по формулам

(16), уравнения (33)

и (31) приводим к следующему виду:

 

= ^ (® .

(34>

= - 4 - 1.

(35)

Оператор 1 в применении к функциям ао, Ф:

 

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

135

 

Если оболочка до нагружения имеет начальные прогибы щ (х, у),

 

то выражения для деформаций получают вид

 

 

 

 

 

 

 

 

 

ди . \ ( дш у

1_ / дщ \ 2

 

 

 

*■*

_

дх +

2 \

дх

)

2 \

дх

)

 

 

 

до

1

 

.

1

/ дш \2

 

1

/

дщ \2

(

 

= - —

 

 

ди)

дни

дщ

 

 

 

(37)

 

_

 

 

 

дщ

 

 

 

 

^

ду "*™дх

*’ "дх™ ’

ду

дх

 

 

ду

 

 

 

здесь

ш — полный

прогиб.

Повторяя

вывод

 

основных

уравнен

 

вместо

выражений

(34) и

(35)

получим

 

 

 

 

 

 

 

(ге — а>о)=Ц

 

.

1

д2Ф

(38)

 

 

 

ф) +!«?■

;

 

-^Г V 4® = — - у - (*• (“ >, ю) —

(а»«, “ >„)] —

 

 

 

(ш —

а»0). (39)

 

При интегрировании приведенных выше линейных или нелинейных уравнении необходимо удовлетворить граничным условиям. Для тор­ цовых сечений замкнутой оболочки эти условия формулируются так: при шарнирном опирании оболочки по краям х = 0, х = /, (направле­ ние координатных осей соответствует рис. 3) для точек краев должно выполняться условие

при защемлении оболочки по краям

ш = 0, ^дх= 0.

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

ох = 0; т = 0.

В случае несмещающнхся кромок следует положить и = 0; о = 0.

Замкнутые круговые оболочки

Сжатие замкнутой оболочки вдоль образующей. Рассмотрим замкну­ тую круговую цилиндрическую оболочку длиной I , шарнирно опертую * по торцам, подвергающуюся сжатию вдоль образующей усилиями р, равномерно распределенными вдоль дуговых кромок (рис. 4).

Об о з н а ч е н и я : К — радиус срединной поверхности оболочки;

к— толщина оболочки. Исследуя устойчивость оболочки в малом, определяем верхнее критическое напряжение; при этом исходным

• Обозначения на схемах закрепления краев см. рнс. 2, гл. 2.

136

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

является дифференциальное уравнение (21) Для рассматриваемого случая оно принимает вид

О а . Е

д*ха

. /

Я3

дхл

(40)

 

Приведем первый вариант решения, в котором предполагаем, что поверхность оболочки после выпучивания является осесимметричной, т. е. что поперечные сечения остаются круговыми. В этом случае про­ гиб ш будет зависеть только от х\ уравнение (40) переходит в следующее:

 

 

Л Лх» + Р

йх*

 

Е

~дхг

=

0.

(41)

 

 

+

Ж

 

 

В соответствии

с граничными условия­

 

 

ми принимаем выражение для

прогиба

 

 

 

,

.

тпх

,

 

 

... .

 

 

ш = / $ т —

 

 

(41а)

 

 

где т —число полуволн изогнутой поверх­

 

 

ности по образующей оболочки.

в урав­

 

 

Подставляя

выражение (41а)

 

 

нение (41), находим

[1]

 

 

 

 

 

 

Р

2 -

 

I

А

 

 

(42)

 

 

Л

+

X2

 

 

 

 

 

 

*

тпК

я/? „

 

 

 

 

 

 

.

здесь А =

—. Приравниваем нулю производную от р

по Л;

а1.Х

при этом считаем т > 1. Получаем следующее выражение для

К = У \ 2 (1 - V 3)

(42а)

Подстановка выражения (42а) в формулу (42) приводит к следующему значению верхнего критического напряжения рв:

 

1

с

Л

(43)

Ра

]Л }(1— V2)

Я

 

При V = 0,3 имеем

 

 

 

 

 

р, = 0,605Е А .

 

(44)

Из выражения (42а)

получаем

длину

полуволны

=

- 5 ------- --------- У М .

(45)

 

V 12(1 — V3)

 

 

В случае весьма короткой оболочки, если

С 1) , нужно в фор-

муле (42) принять т = 1 и пренебречь вторым членом. Тогда

лЮ

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

137

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

. . тях

.

пи

(47)

1В)= [ 51П —

:

где т — число полуволн по образующей оболочки; п — число полных волн вдоль окружности. Исходим из дифференциального уравнения (40); подставляя в это уравнение выражение (47), получим (1]

 

И (

т ая 2

п2 у

 

 

 

Е

гп*л4

 

 

~ П \

~

+ 1 ? )

 

+

 

 

 

 

 

 

 

 

/ т2л2

 

п2 \ -

 

т2л2

 

(48)

 

 

 

 

 

 

 

 

 

 

 

 

 

Введем безразмерные параметры по формулам

 

 

 

РК .

 

 

 

 

 

 

 

 

Л2Л

(49)

Р = ЕК

 

 

 

л/,

 

»

 

Я

 

 

 

 

 

 

Параметры *0*

и г|

выражаются

через

длины полуволн

изогнутой

 

дуги

^

 

лЯ \

 

и

 

по

образующей

|

по формулам

 

 

> -

п

)

 

 

 

 

 

 

 

 

1ц .

 

 

 

 

 

 

 

 

 

 

 

■&=

 

Л

 

 

 

 

 

(50)

 

 

 

1х '

 

 

~

 

Ч

'

 

 

 

 

 

 

 

 

 

 

 

 

Параметр б1 характеризует

очертания

вмятины;

т| — длину полу­

волны 1У. Тогда вместо формул (48) можно записать

 

 

Р ~

 

1

 

(1 +

Ф2)2

 

Л +

б 2

(51)

 

 

 

(1 +

■&*)* V

12(1— -V2)

 

 

 

 

 

Считая числа т и п достаточно

большими, находим

минимум р

из условия

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

др_

 

0,

 

 

 

 

 

 

 

 

 

 

др

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1 +

б 2)3

 

 

 

 

 

 

 

 

Р =

Л-

 

 

 

 

 

 

 

 

 

 

 

 

Из этого условия находим значение р:

 

 

 

 

 

 

 

 

р =

К 1 2 (1 — V2)

 

 

 

и верхнее критическое значение

параметра р

 

 

 

Ре =

 

 

1—

■■

 

0,605.

 

(52)

 

^3 (1 — V2)

 

 

 

 

 

 

138

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

Приведенным решением не устанавливается однозначно форма волнообразования оболочки; вытекает лишь следующее условие, которому должны удовлетворять величины т| и д :

(53)

В случае выпучивания оболочки с образованием квадратных волн (д = 1) получаем г\ = 0,825; тогда

(54)

Верхнее критическое напряжение рв, определяемое в соответствии с формулой (52), в точности совпадает с формулой (43). Следовательно, потеря устойчивости оболочки в малом с образованием вмятин, распо­ ложенных в шахматном порядке, происходит при том же напряжении, что и в случае осесимметричного выпучивания.

Все приведенные выше формулы относятся к случаю п ^ 4. Прак­ тически это соответствует выпучиванию оболочки средней длины. Границы применимости теории оболочек средней длины в случае сжатия определяют, исходя из следующего условия, уточненного по сравнению с выражением (1):

(55)

Выражением (55) определяется область значений I для которых

можно считать справедливой формулу (43).

Рассмотрим случай

слабо выраженного волнообразования (л = 2,

п = 3), когда оболочка

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

При получении расчетных формул для этого случая используют урав­ нение (29). В выражении (28) следует положить

Принимая

где

из уравнения (29) получим [1]

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

139

здесь обозначено е =

и введен прежний параметр к =

. в

Из условия минимизации р по к находим к2 — е1/2л2 (п2— 1)/(1 — у2)1/2 и величину

 

Л _______ I_____

 

п8 — 1

(57)

 

Р" ~ УЗ (1 — у‘) '

"2 + 1

 

 

При V = 0,3

будет

 

 

 

 

Р .* 0,605

~

| .

(58)

В случае п =

2 находим рв = 0,363, что составляет 0,6

от «класси­

ческого» значения, определяемого по формуле (43). При /1 = 3 верхнее критическое напряжение составляет 0,8 от «классического» значения.

При

п2 > 1 получаем

значение ра 0,605, что совпадает с результа­

том

по формуле (43).

 

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

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

сравнима с

толщиной оболочки. Эти вмятины появляются обычно

в процессе

резко выраженного хлопка оболочки. Отсюда вытекает

необходимость решения задачи с позиции нелинейной теории.

В книге

[1] изложено несколько вариантов решения нелинейной

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

интегрирования этого уравнения определяют функцию

напряжений

в срединной поверхности Ф. Находят полную энергию системы

Э = 11е + 1)и - Ш ,

(59)

где 1!с — потенциальная энергия деформации срединной поверхности;

Уи — потенциальная

энергия изгиба; 117— работа внешних

сил. Ве­

личины Ци и Ус определяют по формулам

 

11с =

|

[(V2®)2 - (1 + V) Ц Ф , Ф)1 йх Лу,

(60)

 

|

|(у 2о03 — (1 — V) Ц ш . ш)| (1х <1у,

(61)

 

Р

 

 

где р = 2я/?/, — площадь поверхности оболочки. Далее полную энер­ гию системы варьируют по параметрам прогиба и находят диаграмму равновесных форм оболочки. В процессе выпучивания оболочки число н размеры вмятин являются переменными, поэтому диаграмма равно­ весных форм представляет собой огибающую серии кривых, отвечающих

140

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

тем или иным числам воли. Из сопоставления различных вариантов

диаграммы находят наименьшую величину р — нижнее критическое

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

по методу Ритца приводят к значениям параметра рн, лежащим в пре­ делах 0,182—0,334; позднее получены более низкие значения, доходя­

 

 

щие до 0,0427.

 

 

 

 

 

 

Во всех решениях по методу

 

 

Ритца нижнее

критическое

напря­

 

 

жение

не

зависит

от

отношения

 

 

4 с .

А. Алексеев,

пользуясь

 

 

методом

последовательных

прибли­

 

 

жений,

пришел к выводу, что вели-

 

 

чина рн падает с возрастанием

/? .

 

 

Обратимся к экспериментальным

 

 

данным

[1 ].

Область

эксперимен­

 

 

тальных значений ркр показана на

рис. 5. Величина рв = 0,605

действительно

является верхней

грани­

цей для реальных критических напряжений. Большая

часть

опытов

приводит к значениям

ркр, лежащим выше 0,18.

Некоторые экспери­

ментальные значения

лежат

ниже этой величины; иногда снижение

Л

 

 

 

также явную тенденцию

Ркр достигает 0,06—0,15. Рис. 5 показывает

к падению ркр с увеличением

, что вытекает

из решения С. А.

Алексеева. С увеличением

вероятность появления начальных про­

гибов должна повышаться, а это приводит

к снижению

средней

вели­

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

получаются значения р,

приведенные в табл. 1.

Из таблицы

следует,

что с увеличением

значение р

резко падает.

 

 

 

I. Значения р при осевом сжатии замкнутой

цилиндрической оболочки

Вероятность попа­

 

 

 

 

Я

 

 

 

дания эксперимен­

 

 

 

 

 

 

 

тальной точки

 

 

 

 

А

 

 

 

в вышележащую

250

|

500

750

1000

1500

2000

2500

область в %

90

0,18

1

0,16

0,14

0,13

0,11

1 0,09

0,08

99

0,14

|

0,12

0,10

0,08

0,07

| 0,065

0,06