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

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

186

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

Соотношения Гука получаются в виде приведенных ранее уравне­ ний (8) и (9). Первые два уравнения равновесия (5) также остаются прежними. Третье уравнение равновесия имеет следующую форму:

 

С V5У2® = Чхк ( Ь х +

)

+

 

 

+ аУн ( кУ + -^р г) +

2тА

д2т

+

(243)

 

дх ду

где д — интенсивность поперечной нагрузки.

 

напряжения

Вводя в

уравнение

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

деформаций (241)

по соотношениям упругости (14), получим

 

 

 

_1_ Гд2(?х

2 дН

д*ои

V

/

д2ах

+ 2

д2т

 

Е |_ ду2

дх ду

ду2

\

ду2

дхду +

 

 

 

_1_

Ь (и>,

“») — У ^ 1.

 

(244)

 

 

2

 

 

 

 

 

Выразим напряжения о*, Оу, т в равенствах (243), (244) через функ­ цию напряжений Ф по формулам (16); тогда придем к уравнениям

- у V4» =

1 (а>, ф) + Ук® +

;

(245)

- 1 - у 4ф =

-----(. (со, О))— у*».

 

(246)

Оператор Ь в применении к функциям ш, Ф определяется выраже­ нием (36).

Принимая в уравнениях (240)—(241), (245)—(246) кх = 0, ку =

приходим (при = 0) к приведенным ранее выражениям (30)—(31), (34)—(35)', отвечающим случаю круговой цилиндрической оболочки.

Допустим теперь, что оболочка имеет начальные прогибы о>0 (х, у). Обозначим через ау полный прогиб; выражения для деформаций (240) будут иметь вид

ди .

1

/д и >\2

1 (

дщ \ 2

 

 

 

 

 

- т

{ - д Г )

>

до

 

1

/ а ш \ 2

1 ( дщ \2

(247)

~дУ

 

 

 

2 - ( т г )

 

 

 

:

 

ди ,до_

дш_

дш_

дщ

дщ

 

 

V

~ду ^ дх

дх

* ду

дх

ду

'

 

Основные соотношения (245)—(246) принимают форму

 

- у

V4 (“>— к>0)

= Ь (а . ф) + у*Ф +

;

(248)

V4® = —

(“>. а>) — Ь (“V “'о)] — Ч к (®—®о)- (249>

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

187

Вводя в уравнения (248)—(249) вместо полного прогиба ш дополни­

тельный прогиб

= но — а>0, перепишем их в развернутом виде:

В 4

д2Ф / д*щ

д*щ '\

,

а=ф /

д2щ

д2щ \

-

дх- V ду*

) +

дф ( ' дх2 +

дх2

)

 

( д*щ

 

 

 

 

д-ф

 

(248а)

дхду

V дхду

1 дхду ) + кх

ду2 + ■ку дх3

+ “Г

 

/ д'2щ

\2

 

 

 

д2щ

 

 

Ф “

 

)

дх2

ду2 + 2 дх ду

дх ду

 

д*щ

д*ш0

д-ш.

д'2ы>п

 

а2ш,

г

д2щ

 

дх*

ду2

ду2

* дх'2 ~

 

к х ~дф

Ку

дх2

(249а)

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

Пусть контур оболочки прямоугольный, в этом случае удобно ввести декартовы координаты х, у, откладываемые в основной плоскости 1—4 вдоль сторон контура (рис. 40). Началь­ ное положение точек срединной поверхности определяется координатой г. Прогиб [^отсчи­ тывают от исходной срединной поверхности параллельно оси г. Уравнения, отвечающие данной трактовке, получим, положив в зави­ симостях (248а), (249а) кх = ку= 0 и вводя г вместо ш0; тогда придем к следующим соот­ ношениям:

2_

 

 

 

 

,

\

 

д2Ф

{ д-щ

&г \

Л

У

 

дх2 \ ду2

+

ду2

)

+

ду2

 

+

дх2') ~

 

 

 

о Д'Ф

(

д~У\.

.

 

\

,

Ч .

(250)

 

 

 

дх ду

\

дхду

дх ду )

'

Л 1

 

 

 

 

1

^ (Т ,

_

( &<г>1

 

3-Ш|

У в ч

, 0

дГ-щ

&г_

Е

7 ф

_

дхду )

 

дх2 '

ду2

' г

дх ду

дх ду

 

 

 

д2Ы)1

 

дЧ

 

 

дЧ

 

(251)

 

 

 

~д^~

 

ду1

 

ду°-

дх1

 

 

 

 

 

 

 

 

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

Если центральный угол 2<р, охватываемый наибольшим размером оболочки, настолько мал, что можно принять соз <р & 1, з1п ср = ф, то системы уравнений (248а)—(249а) и (250)—(251) являются эквива­ лентными.

Панель, прямоугольная в плане

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

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

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

следующие граничные условия для краев х =

0, х =

а:

а» = 0;

= 0 ; Ох = 0;

т = 0.

(252)

Аналогичные условия будут для краев у =

0, у =

Ь. Оболочка испы­

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

При решении задачи используем уравнения (248а)—(249а); при­

мем ю0 = 0 и запишем их в виде

 

 

 

Х = - % - ч * щ - Ц Ф ,

щ ) - ь х ^

- к у ^

- -----2 - = 0 ;

(253)

Г =~Ё~ V4® + -§ ~ 1

+

кх1 % г +

к»1& ~ = °-

(254)

Функции прогиба и напряжений аппроксимируем с помощью сле­ дующих выражений:

Щ

=

[ г 51П 31П

о

Ф = А х 51П — 51П

о

.

(255)

 

 

а

а

 

 

Первые

три

граничных

условия

типа (252) удовлетворяются на

всех кромках; четвертое условие — лишь «в среднем»; так, для кромок

ь

х = 0, х = а это условие вы няется в виде — I т йу = 0. Урав-

(X V

а

нения Бубнова-Галеркина в применении к двум уравнениям системы (253)—(254) можно записать

аЬ

|

|

X 51П ■— 51П

йх Ау = 0;

 

о

о

 

 

 

а Ь

 

 

 

|

|

У 5Й1 ■—-^ - й х й у = 0.

(256)

о

о

 

 

 

Подставив в уравнение (256) выражения (253)—(255), после инте­ грирования приходим к системе двух уравнений, из которых получаем следующую зависимость между нагрузкой и прогибом:

+

(257)

 

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

189

здесь введены безразмерные

 

параметры

 

 

 

кха2

 

ЬУЬ2

к

= кх -{- ку\

 

 

 

 

 

 

ку ~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5 =

-

 

 

' аЬ \2

 

 

(258)

 

 

 

 

 

 

 

 

Для квадратной

панели

(а = Ь)

при V = 0,3

будет

 

 

д* = 8.77С3 — 2,46/г*^2 + (0,154/г*2 +

22) ^

(259)

В случае круговой цилиндрической оболочки, если координату г/

отсчитывают вдоль дуги, надо принять кх = 0, ку =

И

параметр А;*

будет к* =

Д/Г Для сферической оболочки кх =

 

 

при с = 5

величина к* =

252

 

 

 

 

 

 

ДЛ

 

 

 

 

 

 

Графики

 

 

 

 

 

Л* приведены на

<7* (5) по формуле (259) для различных

рис. 41 штриховыми линиями. При малых к* величина д* непрерывно

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

производную-^—-. Приравнивая ее нулю, находят значения С» отве­

чающие верхней и нижней критическим нагрузкам:

1 \4

190

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

 

В точке,

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

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

я* ( х + х ) а

при Г = — к*

1^ 6(1 • - V 3)

ь 32

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

и Ян Для оболочек относительно

большой кривизны.

В уточненных решениях функ­ ции и»иФ выражаются в виде рядов

VI г

1 ях

ту

 

1.1

. тх .

ту

 

Ф =

(261)

' з1п----- зт -Ц г“

 

а

Ъ

 

й /

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

зывается существенным, начиная с к* «а 36. Уточненные решения дают несколько меньшие верхние критические нагрузки; при этом заметно возрастает нижняя критическая нагрузка.

На рис. 42 показана зависимость верхней и нижней критических нагрузок от параметра к* для панелей с различным отношением сторон по данным четвертого приближения при V = 0,3. Прощелкивание

будет при к* >

18 для X = 1, при к* >►20,4 для X. =

и ПРН >

> 30 для X =

0,5.

 

Если в решении исходить из аппроксимирующего выражения для Ф в виде Ф = С/ (х) V (у), где V и V — «балочные» функции, удовлетво­ ряющие условию защемления балки по краям, то четвертое граничное условие (252) может быть выполнено во всех точках края. При этом для квадратной панели в уравнении (259) заметно изменяется только коэф­ фициент при р ; он оказывается равным 7,48 вместо 8,77.

Источник: https://studfile.net/preview/19992756/