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

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

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

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

151

давления <7, равномерно распределенного по боковой поверхности, и кручения парами Мк, приложенными по торцам, расчетные формулы имеют вид

где <7о, о* Яо. н — верхнее и нижнее критические давления в случае дей­ ствия только внешнего давления; 50, в и 50> н — верхнее и нижнее крити­ ческие напряжения при простом кручении.

С о в м е с т н о е д е й с т в и е о с е в о г о с ж а т и я , в н е ш ­ н е г о д а в л е н и я и к р у ч е н и я . Комбинируя выражения (97) и (99), получим формулы для определения верхних и нижних кри­ тических напряжений

(100)

З а м к н у т а я о б о л о ч к а , п о д в е р г а ю щ а я с я с о в ­ м е с т н о м у д е й с т в и ю в н е ш н е г о д а в л е н и я , к р у ­ ч е н и я и и з г и б а. При решении задачи об устойчивости оболочки при поперечном изгибе отмечалось, что при сравнительно большой длине ( / . > 4 #) основное значение имеет потеря устойчивости типа сжатия с образованием мелких вмятин в сжатой зоне. Поэтому при рас­ чете на комбинированную нагрузку можно исходить из формул типа (100), подставляя вместо р величину максимального напряжения рг при изгибе. В случае оболочки малой длины < 4/?) должно произойти выпучивание типа кручения с концентрацией вмятин в центральной зоне, при этом необходимо использовать формулу (99), подставляя вместо 5 величину, равную сумме касательных усилий, вызванных кру­ чением, и максимальных касательных усилий от изгиба.

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

152

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

отношение параметра рн =* рн — - нижнего критического напряжения

сжатия при наличии внутреннего давления к параметру р 0>в верхнего критического напряжения для простого сжатия, определяемого по

формуле (52), по оси абсцисс — параметр внутреннего давления д =

=

^

К 3 (1 — V2). Значению д = 0 соответствует рн = 0,18,

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

вдоль образующей.

З а м к н у т а я о б о л о ч к а п р и с о в м е с т н о м д е й ­ с т в и и в н у т р е н н е г о д а в л е н и я и к р у ч е н и я . Верх­ ние критические нагрузки для оболочек средней длины определяют по графику на рис. 17. Указанные на графике величины вычисляют по

формулам: — (83); </ — (65); 50] в — (82); в — (72).

Расчетной является сплошная линия, штриховая линия соответ­ ствует приближенной параболической зависимости

относящейся лишь к большим значениям ц. График показывает, что при

увеличении внутреннего давления щпроисходит некоторое возрастание

верхнего критического напряжения кручения $<,. На графике нанесены экспериментальные точки (черный кружок относится к стальному об­ разцу, остальные — к дуралюминовым). По графику можно вести прак­ тические расчеты, вводя поправочный коэффициент а, несколько более высокий, чем при простом кручении (см. стр. 148), например а = 0,85

при -5 - ^ 2 5 0 .

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

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

153

максимальным нормальным

напряжениям относился

оболочкам

сравнительно большой длины

4^; расчет по максимальным ка­

сательным напряжениям применяют при Ь < 4Я. При комбинирован­ ном нагружении, относящемся к первому случаю, нижние критические значения максимального напряжения изгиба можно определять по

графику на рис. 16, повышая рн (при ? < 0,5) примерно на 25%. Во вто­ ром случае при расчете следует пользоваться графиком на рис. 17, исходя из рекомендаций для случая совместного действия кручения и внутреннего давления.

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

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

Рассмотрим устойчивость конструктивно анизотропных оболочек. Примем, что главные направления жесткости ортотропной цилиндри­ ческой оболочки х, у совпадают с образующей и дугой поперечного сечения. Упругие свойства ортотропиых оболочек характеризуются четырьмя независимыми величинами: модулями упругости Ех и Ег по направлениям х, у, модулем сдвига О и коэффициентом Пуассона ух, отвечающим поперечной деформации вдоль дуги. Второй коэффициент V.,, соответствующий поперечной деформации по направлению х, связан

ссоотношением

154

 

 

 

 

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

 

 

 

 

Уравнения равновесия (5) и (6) в усилиях Мх>Nу, Т и моментах Мх

Му, Н, приходящихся на единицу длины контура, будут:

 

 

 

дМх

 

 

 

В Т

|

 

 

0 ,

 

(101)

 

 

дх

- +

ду

= 0;

дх

+

ду

 

*

 

 

 

дШх

 

2

дх ду +

д,?

+

 

 

4

 

(102)

 

 

дх3

+

«

+

°'

 

 

 

Зависимости закона Гука для деформаций в срединной поверхности

представим в

виде

 

 

 

 

 

 

 

 

 

 

 

 

 

Е*

 

в,

а.*,

о

 

мх

 

 

 

 

е|/

ву

62

0

 

Му

;

 

(103)

 

 

 

V

 

0

0

 

Ьа

 

т

 

 

 

здесь

приняты обозначения

 

 

 

 

 

 

 

 

Л

1

А

 

1

6„ =

-

Уа

 

У1

;

= он

 

= _ЁуГ ’

 

 

 

Я,Л

 

ЕгН

 

 

 

 

 

 

 

 

 

 

 

 

 

(104)

Например, для деформации вх матрица (103) разворачивается в виде

ех = М * + ■М , + о . г =

.

Используя форму матрицы (103), запишем моменты

м х

А

А

0

 

Му =

А

о<%

0

1

н

0

0

2Оа

X

где

ДЛ 3

 

 

А Л3

п _

 

г» _

А — 12(1 — у ,и ,)'

2

12 (1 - г ^ ) ’

Ос =

^

3; А р— Г?1/У2 = Д л^.

(105)

(106)

(107)

Изменения кривизн и*, %ц, я определяют по формулам (3). Введем функцию усилий в срединной поверхности <р по формулам

д2Ф

^

_ д2ср

д2(р

Мх =

 

 

Г = — ■дхду

Функция ф = ФН,

где Ф — функция

напряжений

типа (16). Уравнение (102) можно записать в виде

3%

6%

+

 

д2ср

° 1~д*

 

ду4

I * 2"

 

 

здесь

А

= Оу +

2О0.

 

 

 

(108)

в формулах

(109)

(ПО)

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

155

Уравнение (109) переходит в уравнение (23) гл. 2 для ортотропных пластинок. Уравнение совместности деформаций (4) получает вид

д4Ф , х

| 1

д2и)

= 0 . (111)

64 а*« + 263 дх*ду*

1 Г

дх-

Вводя операторы

 

д4

 

а4

 

Уй =

*!2 дх*

“Ь 2^3

ах2Я , А “Ь ^

д у 4

У4о

д4

+ 2Д

а4

а4

 

дх4

 

дх* ду* + ° 2 ду* •

представим (109) и (111) в форме 4 1 а2ф _

( 112)

а2ш

^< Р + 4 " дх2 = 0 .

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

При решении задач об устойчивости в малом в первое из уравне­ ний (112) следует подставлять интенсивность фиктивной поперечной нагрузки по формуле (19), тогда оно примет вид

4

1

д2Ф

, (

д2т

дгт , л а2и» \

л

 

Я '

ал4

+ Л ( Р*

дхг

+ р и~ду*~ + 2$~ Ш ^ )

= 0 -

Способом, описанным при получении уравнений (34) и (35), прихо­ дим к следующим уравнениям, относящимся к ортотропным оболочкам большого прогиба (при ^ = 0):

ч ),ш = Ц а ,,

ф ) + - 1 - , - 0

- ;

 

4_

1

.

,

.

1

д2и)

(113)

у вф —

х

Мв> ш)

/ Г '“ а*4" ’

 

где I — оператор по формулам

(32)

и (36).

 

 

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

V I)(И ’ - И 'о ) « М » .

<Р ) + Х " 0 " >

 

УбФ --------3- [М ® .

«0 — 1. (щ, О10)) —

(114)

1 а4 (ш — щ)

Ядх*

здесь хв) и ш0 — функции полного и и чального прогибов, как и в урав­ нениях (38), (39).

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

Энергия деформации в срединной поверхности и энергия изгиба

для ортотропной

оболочки

равны

 

^ с = 4 "

П

+

(115)

 

р

 

 

+

( 4 ? ) ’ + , 0 ' ' ( а ? 3 » Я л

< " 6>

где Р = 2лЯЬ — площадь поверхности оболочки. Рассмотрим решения конкретных задач.

У с т о й ч и в о с т ь о р т о ­ т р о п но и о б о л о ч к и пр и

о с е в о м

с ж а т и и

у с и ­

л и я м и

IV, р а в н о

м е р н о

р а с п р е д е л е н н ы м и по т о р ц а м о б о л о ч к и . При­ нимаем следующее аппроксими­ рующее выражение для прогиба:

XV =

с /

пх

я у

 

,

 

/ ( СОЗ

---- С05 -г2- +

 

 

\

 

 

 

.

2пх ,

,

у .

,

\

-+- а соз

— \-Ь соз у-

+

А

),

 

 

*и

 

 

/

 

 

 

 

 

(117)

где /Л-, 1у — длины полуволн вдоль образующей и по дуге; а, Ь и й — постоянные величины. Рассмотрим нелинейную задачу (решение линей­ ной задачи получается как предельный случай). Выражение (117) подставляют во второе из уравнений (ИЗ); в результате интегрирова­ ния его определяют функцию ср. Далее вычисляют полную энергию системы; минимизация энергии по пяти параметрам приводит к системе пяти нелинейных уравнений. Результаты вычислений с помощью элек­ тронных цифровых машин показаны на рис. 18 П ]. На графике (рис. 18) обозначено:

^ +

~2~

гх

ас

Р Г ;

(118)

” 7

аи==7 Ш ;

е0, в — критическая деформация кольцевого выпучивания; у — харак­ теристика ортотропности оболочки»

У =

А Л

(119)

Р262

 

при относительно сильных стрингерах величина у значительна, при сильных шпангоутах — мала; N а — верхнее критическое усилие сжа­ тия, отвечающее случаю осесимметричного кольцевого выпучивания; N — верхнее критическое усилие, соответствующее тому случа!

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

157

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

угольные клетки.

Величина

М0, а будет:

 

 

 

 

 

 

 

^ • в = ~ ж

У ^ ; -

 

 

(120)

 

 

 

 

 

 

 

Введем

обозначение

 

 

 

 

 

 

 

 

 

 

ту

N11 1 /~ Ь Г

 

 

(121)

 

 

 

 

* « • « - —

V “о т*

 

 

 

 

 

 

 

 

 

тогда величине Л/0)в будет соответствовать параметр

М0<в =

Для

ортотропных оболо­

Ун 1

ас = + ;

«и = 0,5

 

чек точки

бифуркации

«кле­

 

точного» типа

и «кольце­

«6

 

 

 

 

вого» типа Л^о, оне совпадают.

Оболочки,усиленньи

 

Выпучивание’в малом должно

 

0.8

стрингеримиу г

 

сопровождаться образованием

0,6 :-е

\

 

 

 

клеток,

причем с

возраста­

 

 

 

ииему отношение

ло

па-

ДО

^Оболочки,усиленные

——

 

 

 

Л'о.о

 

0.2 Изотропныё\

шпангоутами

дает. Закритическне диаграм­

 

 

1—

мы относятся к выпучиванию

оболочки —1___ 1_1 1

«по ромбам»; характер их рез­

0,64-

4 1

а

!

± х

ко меняется в зависимости от

 

 

4 9 16

64

у. В случае усиленных стрин­

 

Рис.

19

 

 

геров (при

значительных у)

 

 

 

 

 

диаграмма имеет за точкой бифуркации обычный падающий участок.

С уменьшением у наклон этого

участка изменяется; начиная с

< — , получается не падающая,

а восходящая характеристика, сме­

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

На рис. 19 показано изменение отношения нижнего критического усилия к верхнему в зависимости от у при ас = 4 и аи = 0,5. Эффект нелинейности оказывается наибольшим для изотропных оболочек. Этот эффект несколько ослабевает в случае оболочки, усиленной в про­ дольном направлении, и резко падает для оболочек, имеющих попереч­ ные подкрепления.

П о д к р е п л е н н а я о б о л о ч к а п о д в е р г а е т с я о д ­

н о в р е м е н н о м у

д е й с т в и ю

о с е в о г о

с ж а т и я

и

в н у т р е н н е г о

д

а в л е н и я .

Результаты решения этой задачи

показаны

на рис. 20.

Интенсивность

 

внутреннего давления обозначим

Через

параметры

N и у будут:

 

 

 

 

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

На графике (рис. 20) изображена зависимость между сжимающим

усилием N и величиной приведенной деформации 7 +

Диаграммы относятся к оболочке, усиленной в осевом направлении. Сверху график ограничен линией, соответствующей случаю осесим­

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

от изотропной. В зависимости от параметра <7меняется и вид начального

участка закритической диаграммы. При малых значениях ц начальному участку диаграммы соответствует падающая характеристика. При зна­ чительном <7 получаются восходящие участки, при <7 = 4 начальный участок

получает форму петли. Можно сделать вывод,

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

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

Круговые панели

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

с

,

.

 

 

(122)

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

159

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

а» =

, .

гппх

ппу

(123)

г $1п

------- 51П — -г-

 

 

а

Ь

 

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

^

о а ,

 

Е

Ь*

(124)

Р ~~ НЬ1 П ^

/?а

я2Ф2 1

где

/ т2Ь

.

п? \

а

 

 

(125)

ф =

( —

+

Т

 

 

 

Минимизируя р по ф, находим

Подставляя формулу (126) в равенство (124), определяем верхнее

критическое напряжение

 

 

 

„Г'Е 4 ~ ;

(127)

К 3(1 — V2)

К ’

 

при V = 0,3 будет

Л

 

 

Ро = 0,605-Е•

 

(128)

Я

 

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

пологой панели следует в формуле (125) принять л = 1, т = - | - ; при

этом будет ф = 2. Формулы (127), (128) справедливы при условии

л/

12 (I — V2) — ? = - > 2 ,

'(129)

"

У

V ЯЛ

 

т. е. если охватываемый панелью центральный угол отвечает нера­ венству

 

(130)

0 > -

V * .

• / 12(1 — V2)

 

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

А л Ю

I /

Ь у-

(131)

Р‘ ~ 4 ЬЧ +

4л2 V

« ) Е'

 

160

 

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

 

 

 

Введем безразмерные

параметры нагрузки

и кривизны

 

 

р

Е

\ к

- к — —

___ —Ь

(132)

 

) '

 

Н

~

«А

 

Тогда из формулы (127) получим

 

 

 

 

 

 

*

 

 

1

 

к а 0,605*.

 

(133)

 

 

 

1 ^ 3

(1 —

V2)

 

 

 

 

 

 

 

 

 

Условие (130) применимости формулы (133) примет вид

 

 

 

 

^

 

а »1

'

 

 

(134)

 

 

 

"*

К 3 ( 1 — V3)

 

 

 

или приближенно к >

12.

к <

12)

из формулы

(131) получим

Для

пологой панели (при

 

 

 

 

я*

 

,

Л3

 

 

(135)

 

 

р ‘

3 (1

V3)

' 4я3

 

 

 

 

 

 

 

 

 

 

Ра = 3'6 +

395

 

 

 

(136)

 

 

 

 

 

 

 

При

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

больших прогибов

панели

используют урав­

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

лениях образуется лишь по одной полуволне ( т =

п = 1). Подставляя ю

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

напряжений Ф

в срединной

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

 

 

 

 

 

 

 

 

 

 

_

1

ф

_

=

1

( , /

2лх

2пи \

.

1

 

- ва

.

лх

^

 

 

_ _ ^ ^ С03_

+ С 0 5 _ ^ ) + - 5

 

У / - Г

ВШ —

X

 

 

 

 

 

 

 

 

пу

 

РУ2

 

 

 

 

 

(137)

 

 

 

 

 

 

 

X 51П —------

 

2Е

 

 

 

 

 

 

 

 

 

 

 

 

 

а

 

 

 

 

 

 

 

 

Статические

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

на

 

кромках: Оу =

0,

т =

0 при

у = 0, у =

Ь\ Ох — р,

т =

0 при х =

0 выполняются для

каждой

кромки лишь в среднем. Далее составляем уравнение Бубнова-Галер- кина

о о

(138) Сюда подставляем выражения для ш, Ф и выполняем интегрирова­

ние; затем находим параметр

[1]

 

 

*

* .

яа у*

10 .

(139)

Р

— Ра ^

з ^

Дтг.а

 

 

 

Зла