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

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

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

Подставляя теперь в выражения (5,1,9) для

Т и ...,

N ь N 2 значения

Ох, Оу,..., тХу из (5,1,7), учитывая (5,1,8) и (5,1,1),

имеем

 

Л_

 

 

2

2

 

 

7*1 =

£г

 

d z =

1 — V2 V дхг

дф

 

)

 

 

 

Eh

(ef0

+

ve° ).

 

 

 

1— v2

 

 

 

 

x x

 

y y }

 

Внося сюда значения

e \x и &

из (5,1,1), получим

 

Т г = - £ h

(е°

+ ve° ) =

Eh

 

Г

 

д*

+ — ( —

Y — kjW +

1 — v2

хх

уу

1— v2

L

 

2 \ дх

)

Аналогично найдем Т2

и S = S \ = S 2.

 

 

 

Таким образом, получаем важные формулы

 

• Гди, ,

v dva

-L |\

 

У

dw у

v2 [ дх

 

ду

 

2 V

 

 

 

f dw

V _L

v

/ ^

 

'iS . 4- v * ft+ _L 1

 

, д у

 

дх

2 '^ ду J + Т [ дх

 

s =

 

£ Л

Г-ди0

dv0

,

dw

4 = л

 

2 ( l « * v ду)

'V дх

 

дх

ду )

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

определены ниже.

Исключая из равенств (5,1,12) неизвестные перемещения срединно­ го слоя и0 и V 0 , можно выразить усилия Т ь Г2, S через прогибы w .

Предварительно заметим, что из (5,1,12) можно получить такие со­ отношения:

Пользуясь этими |равенствами, составим уравнение

(Тг - vT2) +

(Г2 - v7\) - 2 (1 + v)

д хд у

d2w

i

d2w

^

~дф

2

дх2

Г

Но согласно правилам дифференцирования

д2

\

±

(

dw

- V

I

J

д

Г

dw

'

d2w

1\

 

-

(

'

d2w

'\ 2

,

dw

d3w

ду2

L

2

ч

)

 

ду

1

^

дхду

Г

 

\

. дхду

.)

+

~дх

'' дхду2 '

 

 

 

 

~дх

 

 

 

 

 

 

J

 

 

 

д2

г ± { . dw • V

I

J

д

Г

dw

 

& w

'}

 

-

(

d*w

'\ 2

,

dw

d3w

дх2

L 2

 

\

ду

)

 

дх

L

I T

 

дхду

И

 

 

. дхду

,)

ду

дх2ду

 

 

д2

( . dw

 

dw Л _

 

 

f

dbw

 

 

dw

+

dw

 

dPw

\

_

 

 

дхду

\

дх

 

 

ду Г

 

дх

'\

дхду

 

 

ду

"& Г

ду2

)

 

 

 

 

 

d3w

 

 

 

 

 

 

 

 

 

 

 

dPw

dw

 

d3w

 

 

 

 

дх2ду

 

 

 

 

 

 

 

 

 

 

 

 

ду2

дх

 

дхду2 '

Внося это

в предыдущее уравнение,

получим

 

 

 

 

 

 

 

 

д2

( Т

г - у Т

г)

 

д2

(T t

- v

T

J

-

2 (1 +v )

 

d2S

 

 

 

 

ду2

 

 

 

 

 

дх2

 

 

 

 

 

 

 

 

 

 

дхду

 

 

 

=

 

E h

 

 

 

 

 

d*w

 

d2w

 

 

^

 

d2w __^

 

d2w *|

(5,1,14)

 

 

 

 

 

 

 

 

~дх?

 

~ду2

 

 

 

'~ д у 2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дх2 J '

Уравнение (5,1,14) дает искомую связь между усилиями, действующи­

ми в срединной плоскости оболочки, и ее прогибом

до.

 

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

 

 

T x =

дг<р

 

 

s

= - h

- ^ Ц

(5,1,15)

h

 

 

 

ду2 ’

 

 

 

 

дхду

 

уравнение (5,1,14) после подстановки значений (5,1,15) примет вид

- j V2V2<P+ К

d2w

d2w

 

d2w

d2w

= 0,

(5,1,16)

~дф +

^2 ~д&

 

"дх2

~ду*

 

 

 

где

 

 

 

 

 

 

 

 

V2V2(P =

-тт- +

2

о

+

 

 

 

 

дх4

 

дх2ду2

 

 

 

V2V2 — бигармонический оператор

Лапласа.

 

 

Теперь перейдем к построению уравнения равновесия.

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

гивающие усилия h p ( x )

и h r ( у ) ,

приложенные

нормально к кромкам

оболочки, и сдвигающие усилия Их, действующие вдоль кромок.

Выделим

 

из

оболочки

двумя

 

 

 

 

 

 

 

парами взаимно

перпендикулярных

 

 

/

 

 

 

 

плоскостей, параллельных до де-

 

 

dN . )

.

Ni dy r d

формации оболочки плоскостям х о г

 

 

+

 

 

 

У

и y o z , элемент

h d x d y .

 

грани

 

 

 

 

 

 

 

Усилия,

действующие на

(T^ ^ I i d;] dx

 

 

 

 

этого элемента, показаны на рис. 5.4

' г

ду

 

 

 

 

 

При изгибе

срединного

слоя

 

d^ d y \ d i ^ ^ ^ r

 

 

эти усилия

поворачиваются

в про-'

 

у

*

 

|

 

 

странстве. Проектируя их на под­

 

 

 

 

 

 

 

вижные координатные оси и отбра­

 

(T<+ f r dxh

f a

l x

dxh

сывая малые величины более высо­

 

 

 

 

 

 

 

кого порядка, получим

 

 

 

 

 

 

 

Рис. 5.4

 

 

+

dx)dy~Tldy+(Sa+^ d

i d y )

 

d x ~

S id x

N j d y — .d x — ( N 1 + - ^ d x ) d y —

 

d*w

d x = 0,

 

 

 

1

2

dx2

 

\

1

dx

J

2

 

dx2

 

T 2d x + ^ T 2 + - ^ - d y ^ d x + ( s t + ^ - d x ^ d y - S xd y -

(^2 +

~ r d y ^ d x 4 -

• -45- d y — №

4 r

• - Й - dy = °>

 

 

 

dy

J

 

dy*

 

 

 

 

dy*

 

+ ( г ’ + ^

+ г 0 Ч т § Г + Т ^ ) +

+

H—~

d x -f-

d y —

d2w

-I-

 

 

 

dxdy

 

+ ( s * + ^ d y +

S2)

~ S k d y + ( Ni +

ddx

+

*** N l ) dy

+ ( N 2 + ^ - d y — N ^ d x + q ( x , y ) d x d y = 0.

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

dTi

 

dS2

= А \

дх

ду

дх"

дТ2

^

 

хг

ду

 

дх

— -А2 ду*

Третье уравнение после раскрытия скобок будет

 

 

 

 

 

 

 

 

 

d j\

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

+

Т

1

d2w

< ^ У +

-

du .

 

,

1

d2w

 

- f j r

W J x

+

27 ytr ■f

^

dy +

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

dT2

d y-

1

d*w

d x d y

-f 2 S j d y

— -M L d x +

 

 

dy

 

2

dy*

 

 

 

2

dxoy

 

 

 

*Si

d x d y ±

^

- d

x

+

- ^ - d x d y - ^ -

+

 

1

d2w

+

2 S id x

d y

dx

57

2

dxdy

 

 

dx

 

2Ri

 

 

2

dxdy

 

+ J h - dy d x dy

d yd x

dy

2R 2

dy

 

,

dN* d y d x +

 

'

di.

1

d*w

d y +

Щ г~ d xd y +

2

dxdy

" '

dx

q d x d y =

0,

(5,1,19)

Учитывая соотношение. (5,1,10) и отбрасывая величины третьего порядка малости, после сокращений получим

 

1

,

r

d*w

. т

1

т

d*w .

о

 

+

 

 

 

Т 1 — - - г Т 1 — — + т2 - — + 1 2 - Г г - + ^1

dxdy

+

 

 

 

Ri

'

*

dx*

 

 

R i

 

d f

 

 

 

 

 

 

 

4* S z

&w

,

dNx

+ -

dNt_ + ( } = 0 г

 

(5,1,20)

 

 

 

dxdy

+

dx

dy

 

 

 

 

 

ИЛИ

 

 

 

Tt

 

r 2

 

d2w

d2w

2S d2w

 

 

JW l_ +

=

 

 

 

T ,

 

 

dx "t~

dy

 

4

Ri

 

R i

 

dx2

dy

 

dxdy

 

 

Составляя уравнения

моментов

этих усилий

относительно

оси <xv

(рис. 5.5),

имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

(Мх+ - M i -

dx)

dy -

M xdy -

Нфх + ( Я , +

 

 

d y ) d* -

 

 

 

 

-

 

+

- f ; -

d* +

iVa) d y - ^ - d x

= 0 ,

 

(5,1,21)

 

 

 

 

 

 

 

что после упрощения

и учета

(5,1,11)

 

 

 

 

 

 

 

дает

 

dH

 

t-\ d

,

 

 

 

 

 

 

 

 

 

»,

dMi .

 

 

 

 

 

 

 

 

 

 

 

dx

dy

 

dx

v

 

 

 

 

+ i^ r

 

 

 

Аналогично получим

 

(5,1,22)

 

 

L

d x ) d u

 

 

 

 

 

 

 

 

 

 

 

A 2 =

dx

 

dy

 

 

 

 

 

 

 

 

 

 

 

d#

 

 

 

 

 

 

 

 

 

 

Соотношения (5,1,22)

выражают

ра­

 

 

 

 

 

 

 

венство нулю (главного момента

всех

сил, действующих на рассматриваемый элемент оболочки. С помощью

равенств (5,1,22) можно исключить из

(5,1,21)

перерезывающие усилия

N \ и N 2. Д ля этого внесем в (5,1,21)

и ЛГ2 из

(5,1,22). Имеем

 

д2М г

+ 2 - ^ L + ^ . + T J k 1 + J * ! L ) + T 2^ 2 + — )

 

 

 

дх2

 

 

 

 

 

 

 

 

ду2

 

1 \

 

дах2 )

2 \ 2

 

d*/2

)

 

 

 

 

 

 

 

 

 

 

 

2S

д2ш

+

9 =

0.

 

 

 

 

 

 

(5,1,23)

 

 

 

 

 

 

 

 

 

 

dxfy

 

 

 

 

 

 

 

 

 

 

 

Уравнения (5,1,18) и (5,1,23) дают искомую совокупность основных

уравнений

равновесия. Правые

части

уравнений

(5,1,18)

можно

поло-

жить

равными 'нулю,

так

 

 

 

 

d2w

 

d2w

имеют

порядок

 

как величины —— ,

ду2

 

 

 

 

h?w

 

 

h3w

 

 

 

 

 

дх2

9

 

 

 

 

 

 

W

\ т

 

 

 

 

то

правые

части

уравнений

(5,1,18)

----, а

N

 

 

, N 2

 

 

 

 

ъ2

 

~b*~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

имеют порядок

h3w2

h3w*

, в то время

как

дТг

 

dS

 

 

и другие

 

 

 

 

Ьъ

 

а5

 

hw2

 

 

 

 

дх

 

ду

 

 

 

 

 

будут

иметь порядок

hw2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

63

 

а3

 

 

 

 

 

 

 

 

 

 

 

 

Вводя по формулам (5,1,15) функцию напряжений, легко видеть,

что уравнения

(5,1,18)

 

будут

тождественно

удовлетворены. Учитывая

выражения

(5,1,11) для

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

уравнение (5,1,23) перепишем так:

 

 

 

 

 

 

 

 

 

 

 

 

- D (

0%

 

d*w

\

— 2 D ( l — v)

d4w

 

 

d*w

+

v

d4w

+

 

 

~дх*~

 

dx2dy2J

 

 

 

 

дх2ду2

 

 

~dy*

 

дхгду*

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(■

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d V

 

 

- ^ L - h - ^ + q --^ 0,

 

+ h ^ ( k 1 + ^ - ) + h - ^ ( k , + - ^ ) - 2

 

ИЛИ

ду2 \

 

 

dx2 J

 

 

 

дх2

\

 

dy2

)

 

дхду

dxdy

 

 

 

 

 

 

 

 

 

 

 

 

d*w

^ g

&00

d*w

)

 

 

 

 

 

 

 

 

 

 

 

 

- D

(

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

~dy*

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx*

 

dx2dy2 +

 

 

 

 

 

 

 

+ h

а2ф

 

,

а2ф

dPw

 

 

k 9 а2Ф _j_

 

d*w

 

ааф

d2w

+

Q— 0 ,

i

 

I

dy2

dx2

 

 

 

 

— 2

 

dxdy

 

dy'

 

 

 

 

 

dx2

dx2

dy2

 

dxdy

 

 

 

что сокращенно запишем так:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D v V t o - A ^

J

 

2

дх*

О

-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

<Э2Ф

d*w

 

.

д2ф

д2до

о

а*ч>

 

)-< 7 =

0. (5,1,24)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ду*

дх2

 

 

дх2

' ду*

 

дхду

дхду

 

 

 

 

 

 

Уравнение равновесия (5,1,24) является вторым основным уравне­

нием теории гибких пологих оболочек.

основных

уравнений

(5,1,16)

и

Следует отметить, что при выводе

(5,1,24) координатные

 

оси о х

и о у

предполагались

совпадающими

с

главными кривизнами. В более общем случае,

когда

такого

совпаде­

ния нет, уравнения

(5,1,16)

и

(5,1,24)

будут содержать слагаемое

 

d2F

k 12 —

дхду 9

где F = F ( x , у ) — уравнение срединной поверхности оболочки. Этот член учитывает влияние кривизны кручения на деформирование и напряжен-

ное состояние оболочки. Полагая в (5,1.16) и .(5,1,24) начальные кривиз­ ны k\ = 0 и k 2= 0 , получим уравнения для пластин большого прогиба:

v V < p = £ r r — Y - — -— 1.

 

V V f

[v дхду )

дх'

дф

J'

 

 

 

 

dx°~

dy2

 

 

Dy2y2tiy — q = h Ldf

d2w

, д2ф

(pW

 

^2ф

d2w

дх2

dx2

dy2

 

дхду

дхду

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

1

2

2

\

 

U

 

I

и

 

 

I

дх2

 

 

 

 

 

 

 

 

 

 

 

Е

V

V Y

1

 

г

2

дхг

^

 

дф

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D y V ^

 

h ( k г J

■ + ь *

 

) -

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 ду2

 

 

 

 

 

 

 

д2ф

d^w

 

.

d2y

d2w

 

2

-

^ .

д2до )-< 7 = 0. (5,1,25)

 

 

 

% 2

дх2

 

дх2

ду2

 

 

дхду

 

дхду

 

 

 

 

 

 

Связь между усилиями, действующими в срединном

слое, и переме­

щениями определяется формулами (5,4,12):

 

 

 

 

 

 

 

 

 

 

 

Eh

Г

ди0

 

до0

 

 

/

dw '

 

V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

■4- V

 

 

 

 

 

 

.)‘ +

 

 

 

 

 

 

 

 

 

 

1 — V2 1.

дх

 

ду

 

 

\

дх

т

 

 

 

 

 

 

 

 

 

 

Eh

Г dvQ

 

ди0

1

1

/

dw '

 

V

 

 

 

 

 

 

 

 

 

1 — v2 [

дх ■+

V

ду

+

т

V ду

.)’ +

~2

 

 

 

 

 

 

(5,

,26)

 

 

 

 

S =

 

Eh

 

/

ди0

 

до,

 

dw

 

dw

\

 

 

 

 

 

 

 

 

 

\

ду

 

 

—2—

1-------

 

ду

) '

 

 

 

 

 

 

 

 

 

 

2 (1 ■Ф- v)

 

 

дх

 

дх

 

 

 

 

 

 

В случае

 

сферической

оболочки

(k i = k 2= k )

уравнения

(5,1,16)

и

(5,1,25)

примут

вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

— v 2V2(P + k\/*w +

^

 

d2w

/

 

 

**

о ,

 

 

 

 

 

 

 

Е

 

 

 

 

 

 

 

дх2

 

ду2

\

дхду I

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dPw

 

о

д2^

 

d2w

\

 

 

 

 

 

 

 

 

\

ду2

 

дх2

дх2

 

— 2

дхду

дхду

J - < 7 =

0.

 

 

 

 

 

 

ду2

 

 

Для

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

оболочки

(fti=0,

fc2# 0 )

уравнения

будут

такими:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

- j W

< P + f c 2

дх2

\

dh)

d2w

( '

 

у

 

 

 

 

 

 

 

 

~д&

 

 

 

\ .

дхду

)

 

 

 

 

 

 

 

 

 

 

 

 

 

V

 

dPw

 

ааФ

 

dhs)

 

 

 

d2w

\

 

 

 

 

 

 

 

 

 

 

 

2

 

дх2 ■+

дх*

 

ду2

 

дхду

 

- 7

=

0 .

Если знаки

кривизн k {

и k 2

разные,

например, &i>0, k 2< 0 ,

то эти

уравнения

могут

быть

использованы

в

качестве

основных

уравнений

теории гибких оболочек отрицательной гауссовой

кривизны.

 

 

 

 

Опуская в уравнениях (5,1,16) и (5,1,25) нелинейные члены, полу­ чим основные уравнения оболочек малого прогиба:

VV<P +

k 1 —

r + k 2 —— = о,

E

ду2

дх2

к которым сводится линейная техническая теория оболочек.

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

Ш

Ъ 7 ,

а

 

' / / / / / / / ,

У ////Л

s х

в

Рис. 5.6

сительно прогибов здесь необходимо задавать еще по два условия на каждой кромке относительно функции напряжений ф.

Вместо граничных условий относительно ф можно задавать переме­ щения и и v, которые связаны с ш и ф условиями (5,1,26). В нелинейной теории изгиба пластинок понятия шарнирного опира-ния, жесткого за­ щемления и другие требуют некоторого уточнения. Например, схемы закрепления кромок (рис. 5.6) в линейной теории соответствуют поня­ тию шарнирного закрепления кромок. Однако условия работы этих пла­ стин при больших прогибах разные. Граничные условия схемы (а) мо­ гут быть записаны таким образом. Прогиб w и момент на кромках x='Const обращаются в нуль:

w = О, М ± = — D

d2w

d2w

= о,

+ v

 

 

~дх?~

^ду2 ) -

причем в силу первого из этих граничных условий второе обратится в более простое

d2w = 0. ~дх2

Кроме того, так как на контуре не приложено никаких внешних усилий,^ нормальные и касательные усилия на кромках равны нулю:

-*9L = 0,

T= - - ^ L = 0 .

ду2

дхду

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

32 П. М. Огнбалов, М. А. Колтунов

Следовательно, граничные условия для закрепления кромок по схе­ ме (б) запишутся так:

d2w =

о, и =

v = 0.

дх2

 

 

Эти условия говорят о том, что геометрия кромок оболочки при ее

изгибе остается прежней и они не

смещаются.

В третьем случае (в) закрепления

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

будут

 

 

и- const, V = const,

что выражает факт смещения кромок панели при'деформации оболочки параллельно самим себе. Возможны также комбинации этих способов шарнирного закрепления.

Рис. 5.7

Аналогично можно рассуждать и в случае защемленных кромок (рис. 5.7).

Граничные условия,- соответствующие схеме (с), можно записать

так:

 

w =

dw

Л

 

 

---- =

0,

 

агд> _

дх

 

 

т =

-

0,

OV =

= 0,

ду2

 

 

д)еду

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

(d)

 

 

 

w

dw

u =

v =

0.

=

дх

В первом случае (с) при деформации пластин кромки искривляются, во втором (d) — кромки остаются прямолинейными и не смещаются.

Влинейной теории в таком различии необходимости нет, и схемы

(с)и (d) можно считать эквивалентными.

Вслучае, когда на кромках пластины заданы перемещения и и V,

после интегрирования основных уравнений (5,1,16) и (5,1,25) по найден­ ной функции сриш следует составить общие выражения перемещений и и о, которые следует затем подчинить заданным граничным условиям на кромках. Другие условия здесь не рассматриваем.

§ 2. ПАНЕЛЬ ГИБКОЙ ПОЛОГОЙ ОБОЛОЧКИ

Рассмотрим задачи об изгибе и устойчивости панелей гибких поло­ гих оболочек (геометрическая нелинейность).

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

члены второго порядка малости относительно прогибов w = i ( x , у ) . Зада­ ча может быть сведена к интегрированию системы двух нелинейных

уравнений

 

(5,1,16)

и (5,1,25)

 

d2w| .

 

 

 

 

 

 

 

1

V

о

« Л

*

д2Ш .

,

d2w

02W

 

/

dPw

 

Ф =

У

Ф +

« 1

-------------Ь

^ 2 --------

дх2

ду*

 

\ дхду

 

Е

v

v

т

1

ду*

^

2

дх*

 

 

 

 

 

 

 

 

 

 

( к

^

+

k 1 -

°2Ф

\

 

 

 

 

 

 

 

 

 

 

\

дх*

 

ду*

)

 

 

 

 

 

<Э2ф

 

д2w

,

о2ф

0*W

 

с,

О2Ф

 

d*w

\

(5,2,2)

 

 

 

ду2

~~ д х +

дх2

ду*

 

дхду

дхду

= 0

 

 

 

 

)

 

при соответствующих граничных условиях.

упругих

оболочек решают

Большинство задач

нелинейной теории

приближенными методами, причем из-за сложности обычно ограничи­ ваются решением в рядах в первом приближении.

Здесь приведем решения задач об изгибе и устойчивости гибкой пологой оболочки, построенных в первом [3] и в более высоких прибли­ жениях [4] методом Бубнова — Галеркина.

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

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

Приведем некоторые соображения общего характера о решении задач этого типа. Пусть к оболочке, каким-то образом закрепленной на жестком кусочно-гладком кон­ туре, приложены произвольная поперечная нагрузка q(x, y) и нормальные к кромкам сжимающие или растягивающие напряжения, компоненты которых вдоль осей оу и ох будут р(х) и г (у). Требуется определить зависимость между действующими на обо­ лочку внешними усилиями и ее прогибами, не считая последние малыми. Для решения этой задачи необходимо проинтегрировать уравнения (5,2,1) и (5,2,2). Ввиду того что методы точного интегрирования этих уравнений еще не найдены, будем искать их при­ ближенные решения в виде рядов

ф =

А,пп [ип (х) vn (у) -

е (*) - х т .

т

п

fmnXmЕ

(*)

 

 

' - Е

Уп ),

тп

где Л шп и fmn — неизвестные постоянные, а функции Um (x), Vn (y)y Xni(*), Yn {y) выбирают заранее таким образом, чтобы удовлетворялись все статические и геометри­

ческие контурные условия.

Применяя метод Бубнова — Галеркина, внесем эти выражения для ф и w в урав­ нения (5,2,1) и (5,2,2), затем, умножая в соответствии с их физическим смыслом пер­ вое на -вариацию функции ф, а второе на вариацию функции w и учитывая независи­ мость вариаций параметров 6Amn между собой и вариаций параметров 6fmn между собой, проинтегрируем полученные выражения по области, ограниченной контуром обо­ лочки. Из полученной системы нелинейных алгебраических уравнений найдем неизвест­

ные параметры А тип И /тп-

Возможность применения метода Бубнова — Галеркина к операторам Ф и w изу­ чалась в работах И. И. Во-ровича [5] и других авторов. Об этом же рассказано в [6], где приводится доказательство, предложенное А. Р. Ржаницыным.

По физическому смыслу уравнения неразрывности деформаций i(5,2,l) и равнове­ сия (5,2,2) должны допускать по одному (решению и являться условиями экстремума некоторого функционала Q(o>, ф). Предположим, что такой функционал существует. Необходимые условия экстремума функционала Q(w, ф) по w и ф

SwQ (W, ф) =

о,

6ф<2 (гг), ф) =

0,

должны совпадать с уравнениями (5,2,2) и (5,2,1), т. е.

 

(ш, ф) = Ф (ш, ф),

6WQ (ш, ф) =

W (ш, ф).

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

 

 

 

в®[V ? (“>. ф)] = бшФ (ой,

ф),

[dwQ (ш, ф)] = вфи7 (гг), ф).

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

нения Ф = 0

и w = 0

являлись условиями

экстремума

функционала

Q

 

 

 

 

 

 

ашФ(гг>, ф) = 6ф№ (го, ф).

 

 

 

 

 

 

Нетрудно убедиться,

что для

уравнений

(5,2,1)

и (5,2,2)

это

равенство

выполняется.

В самом деле, уравнение Эйлера дает для первого из них:

 

 

 

 

6ШФ

дФ

д

Г дФ

\

 

 

д

Г дФ \

 

й2

 

/

дФ \

 

dw

дх

Ч

dwx

)

 

 

ду

\

dwy )

дх2

\

dwxx )

 

 

 

 

 

\

 

д2

(

дф

\

а

д2

 

(

дф

\

к

_ Е Г

&

/

_

а2ш

^ 2

 

Ф

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дхду \

dwxy

)

'

ду2

\

dwyy

J

 

 

|_

 

\

 

дУ2

/

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

УУ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

д2

{ _ & w _ \

 

 

d*w

 

 

+ 2 _ L A J ! L ] A J L

 

2Е

дх2ду2

 

 

 

 

дхду

\

дхду

)

 

 

ду2

\

дх2 )

 

 

 

 

 

 

 

 

 

 

Аналогично

получим

 

 

а

2

 

 

 

 

 

Л

 

а/ 2

 

 

 

\

,

 

 

 

 

h

Г

/

d2w

\

k

 

 

d2w

 

 

 

 

D

[

дх2

\

ду2

)

+

2

дхду

\

2

дхду

)

^

 

 

 

 

*

а2

/

d2w

 

\

“I

 

h

d4w

 

 

 

 

 

 

 

 

ду2

 

Vдх2

 

)

J =

2~D

дх2ду2'

 

 

 

 

 

Следовательно, это равенство выполняется с точностью до постоянного множителя. Умножив (5,2,1) на h/ED, получим точное выполнение условия применимости метода Бубнова — Галеркина к решению нелинейных задач упругих оболочек.

Рассмотрим решение задачи в первом приближении. Будем 'считать,

что функции £/(*),..., Y ( y )

выбраны так,

что

все

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

удовлетворены. Решение запишем в виде

 

 

( 5 .

<t = A[ U( x ) V{ y ) - Q{ x ) - X( y ) ] .

 

w = f X ( x ) Y ( y ) ,

 

 

(5,2,4)

где функции 0 (х) и Х ( у ) выбираются так, чтобы

 

 

6" (*) =

JL р (X),

Г (у) =

.1

г ( у ) .

(5,2,5)

Составляя уравнения Бубнова — Галеркина

 

 

 

Сj Ф ((/, V, X , Y , A ,

f , р , г) U V d x d y =

0,

(О)

\ \ W ( U , V, X, Y , A , f , р, г, q ) X Y d x d y = 0,

(G)