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

книги / Оболочки и пластины

..pdf
Скачиваний:
21
Добавлен:
19.11.2023
Размер:
71.66 Mб
Скачать

D

K

T Л * 4 **г

2Я» У

Х х у R*

ьтг

 

 

 

 

 

Eh

 

 

 

 

 

+ Т

* (’—

“ ) <2 - ^> .*3,

~

+ 5, -£ - ) .

(5.8,18)

причем функция ф выражена через (5,8,6), а х нмеет значение

 

 

1 + / А +

 

ЗАсо

. (5,8,19)

 

 

 

 

 

 

 

4( 1- 0+ -5- а К *) .

Формулы (5,8,18) значительно упрощаются при условии (5,8,13), когда

S, =0, в, = —; 0,.= —

«М, = - Ю (-fiL + J L -).

6Т, = - £ / . ( 1 - «

+ Y

•>Vk) + . (5,8,20)

Для решения задачи об устойчивости остается написать дифферен­

циальное уравнение равновесия:

 

 

 

 

d26Mj ,

rp

d?w

, 6Т2

0.

(5,8,21)

т

ч

, «

-г ■

dx2

 

dx2

 

 

 

Она легко интегрируется как в общем случае,

когда

6Mi и бГ2 опреде­

ляются формулами (5,8,18), так и в случае 5У= 0. Внося сюда (5,8,20), получим

d4w

2аи/2

л .

да = 0,

(5,8,22)

"d*4" +

kER2 V 3

dx2

 

 

где через £ обозначена гибкость1

3ft

£

h '

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

да = С sin ах,

при этом наименьшее значение критической силы получается из условия

£

Исследование других случаев устойчивости основано «а применении либо уравнения (5,7,22), либо (5,8,21) при значениях 6МЬ бМ 2 согласно (5,8,18); оно вполне, аналогично исследованию соответствующих упругих задач, поскольку дифференциальное уравнение (5,8,21) линейно и со­ держит только четные производные от да.

1 При этом отброшены члены порядка hj2 сравнительно с единицей.

Весьма 'всесторонне задача устойчивости круговой цилиндрической оболочки исследована В. И. Королевым; им обследована устойчивость оболочки 'при осевом сжатии, при одновременном действии внутреннего давления и осевого сжатия, при осевом сжатии и наличии предваритель­

ного внутреннего давления и т. д.

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

§9. ПРИБЛИЖЕННЫЙ МЕТОД ИССЛЕДОВАНИЯ УСТОЙЧИВОСТИ ОБОЛОЧЕК

СУЧЕТОМ ФИЗИЧЕСКОЙ И ГЕОМЕТРИЧЕСКОЙ НЕЛИНЕЙНОСТИ

Изучение устойчивости оболочки «в большом» учитывает только геометрическую нелинейность, что означает, как отмечено выше, удер­ жание в рядах для выражения деформаций через перемещения квадра­ тичных Членов, например:

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

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

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

ки. Такой учет

приводит к весьма громоздким вычислениям.

 

Задача

несколько упрощается,

если строить теорию, основываясь

на гипотезе

о

нелинейно-упругом Материале, предполагая совпадение

законов нагрузки

и разгрузки.

 

 

 

 

 

При такой

постановке

задачи

упруго-пластические

свойства

мате­

риала не рассматриваются.

Тем не менее

результаты

решения

могут

быть применены

к широкому классу

материалов (например, сплавы,

пластмассы, сталь в зоне упрочнения

в случае активной дефоомаиии

и др.).

 

 

 

 

 

 

 

*

 

Приведем

результаты

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

И. А. Лукаша [27], который

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

гофа Лява, а материал предполагает нелинейно-упругим полагая совпадающими законы:

что имеет место для несжимаемого материала (v = 0,5).

Зависимость

{5,9,1) может’быть записана в достаточно общем виде:

 

 

c r - ' f Л,Л.

 

 

/=1

 

хде А {

и k i —(некоторые константы. При рассмотрении частных случаев

этой зависимости

 

 

а) а = А г к

(5 ,9,2)

А и k

определяют опытным'путем из диаграммы растяжения (сжатия),

 

б) а = Се (1 — ае)

(5 ,9 ,3 )

постоянные С, пг и а можно найти по условной диаграмме растяжения. Геометрическая нелинейность учитывается, как и -в § 8 этой главы, вводом квадратичных членов в выражениях деформаций через переме­

щения срединной поверхности:

 

 

ди

 

1

f

dw \ 2

— kjW,

 

 

 

 

~дх

' "2V"&rJ

 

 

 

 

 

до

,

1

/

dw \ 2

k 2W,

 

(5,9,4)

 

 

ду

 

2

\

ду

)

 

 

^ху

 

ди

 

 

dv

^

 

dw

dw

 

 

 

 

¥

 

 

дх

 

 

 

дх

ду

 

 

 

 

 

 

 

 

 

 

 

 

 

На основании гипотезы Кирхгофа—Лява

деформации элемента

«оболочки на расстоянии г от срединной поверхности будут:

 

е* = ех +

zx,,

 

 

 

 

 

 

)•

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ег/ =

е у +

 

 

 

 

 

 

 

 

) '

 

 

(5,9,5)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

е х у ех у + %2 Х ху>

(

.

х

у

-

d*w

\

 

\ у

 

 

дхду /

 

 

 

 

 

 

 

 

 

 

 

 

 

Внося (5,9,4) и (5,9,5) в выражение для интенсивности деформаций

(при v = 0,5)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

f

| /

е* +

еу +

 

 

+

Т

е*»’

(5,9,6)

 

 

 

получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

еп ~

у

з“

Г

 

4- Ь22 -(- 6322,

 

(5,9,7)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= ех ^

 

ехеУ +

е1у'

 

 

 

h

=

x l

+ x l

+ х хк у + *1у-

 

 

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

 

 

 

 

 

 

t

=

—— (— 1'С t

+ О»

 

(5,9,9)

имеем

 

 

h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ех =

“f" 7t±t,

Еу

Су -j-

 

®ху

^12^»

(5,9,10)

где

 

 

 

 

 

 

 

 

 

«1 =

Y ;

*2 = Т

*12 =

 

(5,9,11)

 

 

Далее

 

 

 

 

 

 

 

 

 

 

 

= -^7=“

+ V

+

V 2.

 

(5,9,12)

где

 

 

 

 

 

 

 

 

 

 

^ = ^ + в2 + вл + 4 - ^ *

 

 

Ь2 =

2 ^

+

2еух2 +

ехк 2 +

еухх + у

etl/x12,

(5,9,13)

 

63 =

«2 + И2 + XlXa +

X2

 

 

 

_Д_.

 

 

Энергия деформации нелинейно-упругого

и

упруго-пластического тела

согласно {25] дается формулой

 

 

 

 

 

 

v = n

i [0f <vfe> + ' T

] ‘tai!^

(5,9,14)

 

 

где

t = . , + e _ e , + . , + t,.

Для несжимаемого тела объемную деформацию 0= 0 шем в виде

 

v = Ш [ К

* - ] ,h“, y d z

 

 

0

 

или

с учетом au= a u(eu)

 

 

 

V = Ш

[ I f ^

^ и ]d* ^ 2-

 

 

6

 

Работа внешних сил № определяется по формуле

 

^ = | |

(?*“ + <7yu + q2w) d x d y

(Чх,

Чу, 4z — составляющие нагрузки по направлениям х,

и (5,9,14) запи-

(5,9,15)

(5,9,16)

(5,9,17)

у, z ) .

Приравнивая вариацию внешней работы вариации работы внутрен­ них сил, получим соотношение

б W = 6V,

(5,9,18)

которое связывает между собой четыре функции.

1. Функцию F 0 срединной поверхности оболочки при начальном на­ гружении (эта функция входит в (5,9,18) через кривизны пологой обо­ лочки) :

d>F0

PFo

.

= CPFQ

дх2

ду* ’

ху

дхду

входящих в (5,9,4) и (5,9,6).

 

 

 

2. Функцию F q срединной поверхности

оболочки при нагрузке q.

В уравнение (5,9,18) эта функция входит через перемещения и, v, w.

3.Функцию нагрузки q (х, у, z ) .

4.Функцию о(еи), которая описывает физические свойства мате­ риала. •

Пусть заданы функция ou=o(eu), форма оболочки до нагружения

F o = F 0 (x, у ) и функция

нагрузки q ( x , y ) ; требуется определить форму

поверхности оболочки

после нагружения. Рассмотрим приближенное

решение этой задачи, предполагая степенную зависимость напряжений от деформаций. Подставив (5,9,2) в (5,9,15), получим

V =

^

^ ' d x d y d z .

(5,9,19)

Здесь тройной интеграл

распространен на

весь объем оболочки.

При k = \ и А — Е получим энергию деформации линейно-упругого тела:

 

у

= ±

^

е \

d x dy d z .

 

(5,9,20)

При fe=.0, A — a s будем

иметь энергию деформации жестко-пластическо­

го тела:

 

 

 

 

 

 

 

 

 

v

=

 

Ш e«d x d y d z -

 

(5,9,21)

Внесем © (5,9,19) значение интенсивности деформации из (5,9,12)

и, учи­

тывая (5,9,9), получим

выражение для энергии деформации:

 

2kAh

 

+ а

+ Ь

+ 1

 

_*±i

 

 

 

 

 

(bx +

2

dxd yd t .

(5,9,22);

 

k + \

 

 

Ъ4 + b 3t2)

 

 

 

 

 

 

 

 

(/г + 1) 3

 

z_a

—b

—1

 

 

 

 

Введя обозначения

 

 

 

 

 

 

 

 

 

в2 =

b± -f~ bot -{- b3t2,

 

(5,9,23)

энергию деформации запишем в виде

 

 

 

 

у = ----- 2M/t к + г Г

^

Г eb+' dxdydt.

 

(5,9,24)

 

(*+1)3

2

Д

^

Д

 

 

Интегрирование по t можно произвести по формуле Симпсона:

,*+. d t = А

[

 

 

+

2eg+> j = -i- (e*+* +

 

+

4<frH), (5,9,25)

тде значения

 

 

 

 

 

 

 

 

 

 

 

е а =

V b x +

b2 +

b 3>

е0 =

1f b ,

еа = V ^

b^ +

b^

(5,9,26)

получаются из (5,9,23)

после подстанов-ки в нее соответственно

 

 

 

t = + \ ,

t = 0, f = — 1.

 

 

 

 

Окончательно для энергии деформации получим

 

 

 

 

V “

i k - T T T

?

?

( е £

+ | + с 5 +<1

+') ы

У’

( 5 ’9 - 2 7 >

 

3—

 

 

-

 

 

 

 

 

 

здесь двойной интеграл

берется по всей

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

оболочки. Вы­

полняя интегрирование,

получим

 

 

 

 

 

 

I =

^ F d x d y =

±

F O +

± ( F 1 + F2 +

F S + F J ,

(5,9,27')

 

s d

 

 

 

 

 

 

 

 

 

 

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

1 =

(8F 0 + F , + F 2 + F 3 +

F a).

(5,9,28)

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

 

 

 

 

 

 

 

 

г

 

ш/

пх

 

 

 

 

F

=-■=cos ——cos---- ,

 

 

 

 

 

 

2а

2а

 

 

 

ТО

 

 

 

 

 

 

 

/ = Г

Г

cos

cos

d x d y

=

л2

(5,9,29)

.J

J

 

v

■будет точным значением интеграла. По формуле

(5,9,28)

получим

Fo =

4 г >

Fi =

F2 =

F4 = 0; F S =

i;

 

' - f ( * T

+ 0 - i - -

 

 

что по сравнению с (5,9,29) дает ошибку в 3%. Учитывая (5,9,27), вы­ ражению для энергии можно придать вид

2k

( * + 1) 3 2

Здесь

 

F. = е'1+• +

+ 4^+>, i = О,

 

1,2, 3, 4,

(5,9,31)

 

 

 

а величины <?в, ен, ео определяют по формулам

(5,9,26). Если действует

только

одна поперечная

нагрузка, то работа

внешних сил

будет

равна

 

 

 

 

+ а

-\-Ь

 

 

 

 

 

 

 

 

U7 = j'

1* q w d x d y .

 

 

(5,9,32)

 

 

 

 

—а —b

 

 

 

 

 

Полная энергия системы

 

 

 

 

 

 

 

 

 

 

U = U ( u , v , w ) =

V — W

 

(5,9,33)

Представим перемещения в виде рядов:

 

 

 

 

 

 

 

« = £

c iuh

v = Y

i C'ioi ,

W =

£

W,W„

(5,9,34)

где Ui,

и W{ — задаваемые

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

граничным

условиям, а С*, С / и ш, — искомые коэффициенты.

Внося

(5,9,34)

©

(5,9,33) и сбставляя условия экстремума энергии

U

 

 

 

^

=

0 ; ^

= 0 , ^ =

 

0,

(5,9,35)

 

 

dQ

ас!-

 

дан

 

 

 

 

получим системы

алгебраических уравнений для определения коэффи­

циентов C i , C i и

Wi. Эта система нелинейных уравнений о дробными

показателями может быть решена приближенными

или графическими

методами. После определения

коэффициентов С*, Ci

и w { «нетрудно из

соотношений (5,9,5) и (5,9,4)

найти

функции перемещений и дефор­

маций.

очень

важно надлежащим образом вы­

При решении этой задачи

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

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

Щ = щ (*) щ (У); Vi = Vj (x)Vi (у); щ = щ (х) w (у).

(5,9,36)

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

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

 

e- = 7

T (8* + 8*)-

(5,9,37)

Внося сюда значения деформаций из (5,9,5), получим

 

е„ ^

(ех + хд.г +

еу + к уг) = у = - (е + кг);

(5,9,38)

здесь

 

 

 

Внося (5,9,38) в (5,9,19), получим следующее выражение для рабо­ ты внутренних сил:

+ J L ^ 2

V =

Л2*+1

 

 

+

Х2)*+' d xdy d z ,

(5,9,40)

 

*+1

 

 

 

 

 

 

 

 

 

 

 

 

( Л +

2

_л_

 

 

 

 

 

 

1 )3

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

или после интегрирования по z

 

 

 

 

 

 

 

 

 

/

h

Y + 2

f

h

\A + 2

 

2*+,

г

г { е +

* Т

)

-

( е-

ит

)

(5,9,41)

v =

*+l_

 

 

 

 

 

 

d x d y,

(fc+l)(A+2)3

2

 

 

 

 

 

 

 

 

где интегрирование распространено по всей площадке опорного плана оболочки. Найдем зависимосгь между нагрузкой и прогибом. Ограничи­ ваясь в (5,9,27) одной нулевой точкой, получим

 

 

 

 

(

,

h \ Л+ 2

(

h \ k+ 2

 

 

V =

 

2*+4

^ • + >т )

—у 1 - ”* Т )

(5,9,42)

 

 

 

k j ±

 

 

 

*0

 

 

 

 

 

 

 

 

 

 

 

 

(Л+1)(Л + ,2)3

2

 

 

 

 

 

 

Здесь во и к — значения -величин

(5,9,39)

в точке

( а / 2, 6/2). Далее, учи­

тывая

(5,9,4),

будем

иметь

 

 

 

 

 

 

 

 

 

е ^ =

 

 

 

 

 

Шп

(5,9,43)

 

 

С,) а 0 (kx -f- k y) w 0w 0

-|— — е0,

где а0, w 0, b0

н е — неизвестные коэффициенты.

 

 

Здесь

принято

 

 

 

 

 

 

 

 

 

 

 

и =

Си, и = C xv,

w

= w 0w .

 

(5,9,44)

Далее допустим, что связь между коэффициентами С, С, и w 0 такая же, как и в упругой теории. Положим

Щ

с = u>oMi Y а2>

 

(5,9,45)

WQ

 

C I = w 0k b i -----— b2,

 

где a u a 2, b i, b 2 — подлежащие определению коэффициенты. Из

(5,9,45)

имеем

 

С + C i — w 0 (кха г + k yb t) -----— (а2 -f b2).

(5,9,46)

Внося (5,9,46) в (5,9,44), получим

 

wl

Здесь

S = а0 (kxCii kyb\) —•.

b0— («г + ьг) а 0 = т.

(5,9,48)

kx + k y

 

 

Принимая во 'внимание (5,9,17) и (5,9,29) для работы внешних сил и учитывая приведенные зависимости, получим следующее выражение для полной энергии системы:

 

 

(

, Л NH-2

/

ft y+2

 

и =

2k+iAab

( e

+ х т )

- ( ев^ »

т )

ш ъ

 

k+3

 

*0

 

qw o ■

 

(ft 4- 1) (ft -b 2)

2

 

 

 

(5,9,49)

 

3

 

 

 

взяв производную no Wo, получим связь между нагрузкой и прогибом, которая после подстановки значений (5,9,44) и (5,9,47) примет следую­ щий вид (формула написана для цилиндрической квадратной в плане оболочки при /г„= 0, Ь = а) :

2*яМ

. В*+! (Сю, (ft -f 2) -

В] -

В*+1 [Сш0 (ft + 2) - В) /с п спч

q -

Ш

 

Ж

. (t>,y,W)

(ft + 1)(ft+ 2) 3

2

 

 

 

где

 

 

 

 

 

В — w 0kxs Н

~~ т

— — С0до0,

 

 

2

 

 

 

В = w 0kxs +

т ~ ~

C 0w 0,

(5,9,51)

С = w 0b0 w 0kx + Св,

С = Ш060 — й»о^л- — \ С0,

величины s, т , b, С и ш0 по-прежнему неизвестны и подлежат опреде­ лению.

Введем безразмерные параметры:

_

д (2а)4 .

w0

_ t .

 

£ft*

ft

(5,9,52)

_

1

-

2/o

_

2g0ft

 

R

^

a*

 

a2

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

св

а2

Подставив (5,9,53) и (5,9,52) в (5,9,50), получим

 

 

 

s =

5 [ 2s5»

+ f £ + l

]

E= - ? ' ,'5;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(5,9,54)

 

 

 

С = А ( 6. | _ . 2^ „ _ ! ) = А с-;

 

 

 

 

 

С = А ( а д - 2^ 1л_ А ) = А ?

 

 

 

 

После подстановки (5,9,52) и (5,9,54)

в (5,9,50)

будем иметь

где

Р =

D k { { Ь У + 1[С* (k + 2 ) - Ь'\ -

( Ь У + '

[С* ф

+

2) -

&*]}£*,

(5,9,55)

 

 

 

 

6* = 2s£0+ -£-£ + - |- ;

 

 

 

 

 

 

 

 

6* =

2S| 0-f y l

- -

1 ;

 

 

 

((5,9,56)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

с* =

60c,— 2 w 0lo ■

■'О .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c* =

b0l - 2 w

0l 0 - ^

f - ,

 

 

 

 

 

 

 

 

A , = 4

JI2.*2*+4

 

а

\2(1 -Л )

(5,9,57)

 

 

 

 

 

 

4+3 -

Ш

 

 

 

 

 

 

(£+!)(£+2)3 ^ C0

 

 

Формула

(5,9,55) выражает общую зависимость между безразмер­

ной 1нагрузкой_р и_пр_огибом £ с точностью до неизвестных пока коэф­

фициентов m ,

Ь0, S ,

w 0, С0Эти коэффициенты

найдем из следующих

двух

условий.

 

o s и £о = 0 (формула (5,9,55) должна давать решение

для

1. При k =

0, A =

жестко-пластической пластины:

 

 

 

 

 

 

 

 

 

 

 

 

я 4 -4

osw0h

 

 

 

 

(5,9,58)

 

 

 

 

 

1 6 / 3 a2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2.

При &= 0 и Л = £

формула

(5,9,55)

должна дать

решение для

упругой

оболочки:

 

 

 

 

 

 

 

 

 

 

 

 

 

р = с^з + аа| а10 +

Озй§ +

а4Е.

 

(5,9,59)

Составляя из формулы (5,9,55) эти условия 1) и 2) и приравнивая соответствующие коэффициенты, получим систему пяти уравнений для определения пяти неизвестных коэффициентов. Решая эту систему, най­

дем следующие значения коэффициентов m , b 0, S , ш0, С0: