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

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

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

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

Доказательство эквивалентности задач I и II можно найти в книге Куранта и Гильберта [80].

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

малом изменении усилий и моментов 6ГЬ бГ2..., при

котором

функции Г1 + 6 Г1,

Гг+бГг,

удовлетворяют уравнениям равновесия типа

(2,16,9)

и

статическим усло­

виям типа (2,16,10), вариация функционала (2,16,12) равна нулю.

 

 

(2,3,16).

Уравнениями Эйлера задачи II являются уравнения

неразрывности

Докажем это утверждение для того же частного случая

линейной

деформации

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

Итак, требуется составить уравнения Эйлера и естественные граничные условия

для функционала

 

 

Эц =

§ W d V -2 ,

(2,16,14)

функции сравнения которого Ти Т2,

подчинены условиям

(2,16,9) и (2,16,10), и на

контуре заданы перемещения

 

 

и = v = w = 0 .

(2,16,15)

Из равенства (2,16,15) следует равенство нулю работы сил реакций на переме­ щениях контура, т. е.

2 = 0 .

Таким образом, требуется найти минимум функционала

Эц =

 

 

(T&i ф

 

Ф 5б12 •ф'

 

М 2х.2 ■ф’ 2f/Xjji) fdxdO

при дополнительных условиях

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dTi

dS

4-'■‘71 =

0 ,

 

 

 

 

 

 

 

 

Г - дх * д0

 

 

 

 

 

 

 

3S

дТ2

дН

1

а м1,

 

 

 

 

 

 

 

дх

^ 60

" * 2

 

ь

 

■ ф Щ - О ,

 

 

 

 

 

 

 

 

 

CD

 

 

 

 

 

 

 

т

дх*

 

1

д*М2,

2

д»Н

^

0.

 

 

 

 

а

г2

дО*

 

г

дхдО

 

Решать задачу будем методом множителей Лагранжа [18].

 

на

М * . 0)

Умножим первое из

уравнений

(2,16,9)

на

Х\(х,

0), второе

на Я3 (х, 0)

(Я< — множители Лагранжа)

и образуем

функционал

 

 

 

Э

 

2

я(Tie,. ф Г,е4 ф Se,2 ф

M ix, ф* у V I ф 2ЯХ12) rdxdQ ф“

г

 

дтг

as

Ф rqi>| "Ф^2 (f

dS

 

дТ2

дН

 

1

дМ%

+И1М

/

1

ае

diML

 

1

д*М2

2

д*Н

*

у

60

\

 

дх

 

 

 

/

ах

*

ае

* 2 ах

г

л

1 1

_

 

 

 

 

 

 

 

 

 

| rdxdQ.

"Г ^3 1

 

*2

дх2

 

Г2

аеа

 

г

ахае

 

Яп) _

 

 

 

 

 

 

 

 

Преобразуем члены, содержащие производные Ти Т2, ..., например,

 

ИXi

 

 

 

 

 

 

 

01

 

 

 

 

 

 

х0

 

 

 

 

 

 

 

(2,16,17)

(2,16,9)

и третье

■ф- r q ^ ф

(2,16,18)

после чего проварьируем (2,16,18) с учетом связи между деформациями и усилиямимоментами (2,7,10), считая функции сравнения Ти Т2, и их производные постоян­ ными на контуре. В результате имеем

 

б^ =

И

К 81"

■r I t e ) 6Ti *

 

{ч ~

1 »

+ ~ 7 h ) &Tt

 

(£ i a _ r " ^ ' _ l e

L) 6 S 4 '

 

 

 

 

 

 

 

 

 

 

 

 

) вм* +

 

Ф (V2X12 — 2

dx

 

r

dxdQ

) б я ]J rdxdQ =

0.

 

Считая вариации 6T\y 6T2,

независимыми,

получим

уравнения

Эйлера для

функцио-

нала

(2,16,18):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дк±

 

 

dk2

 

1

 

 

 

 

 

 

 

б! = Г дх

е2 — дд

 

^3»

 

 

(2,16,19)

 

 

 

 

 

дк2

 

дк\

 

 

 

 

д*к3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8ia= Г ——4

----

 

 

X1

 

 

 

 

 

 

 

 

 

 

 

дх

"

а е

 

 

 

 

ах2 ’

 

 

 

 

 

акг к

 

 

 

 

 

 

 

 

 

 

1

д*к3

 

 

 

 

дк2

1

д*к3

 

 

Г*2 = -ГГ- 'f

 

 

 

*1 2 =

дх

 

а е 2дхдв*

 

 

 

а

е

 

 

 

 

 

 

 

 

 

Чтобы получить уравнения Эйлера для

функционала

 

(2,16,17),

исключим из

(2,16,18)

к\у к2у

к3. Получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дх*

дх12

 

 

 

 

 

 

 

 

 

 

 

 

дх

а

е

 

= 0 ,

 

 

 

 

 

 

 

де12

 

1

дех

 

дх12

 

дхг

 

 

 

 

 

дх

 

г

а е

 

 

 

дх

' а ге = о.

 

 

 

 

 

 

 

дЧ2

3*812

 

 

1

 

а2е1

 

 

 

 

 

Xi4~ г

дх2

дхдв

 

Г

'

а

е 2= 0 ,

 

 

(2,16,20)

которые полностью совпадают с уравнениями неразрывности срединной поверхности (2,3,16), если в последние подставить условия (2,16,2).

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

Рассмотрим панель пологой

оболочки x i ^ x ^ . x 2i y i ^ y ^ y 2, нагруженную нор­

мальным давлением qy с граничными условиями

 

dw

и =

v = 0 на кромке х =

хп,

w = —— =

дх

 

 

(2,16,21)

dw

 

 

и =

v = 0 на кромке у = уп , п = 0, 1 .

w = — =

<>У

 

 

 

Зависимость между перемещениями и деформациями

описывается формулами (2,16,22).

Вариационная задача I формулируется следующим образом: составить уравнения

Эйлера для функционала

 

 

 

 

 

3 , = j Wdv Л ,

(2,16,11)

на функции сравнения которого

 

 

с

о

о

 

ехх* еху* ®уу» **» Kyt KXyt u0 t v0 t w

наложены условия

дх

kiW +

- 2

(

. 1 г)

БХ Х =

0

ди0

 

 

 

1

/

' dw V

0

dv0

k2W

 

1 i/ dw

8УУ =

ду

 

2

(

 

 

 

 

 

 

 

о

dv0

.

ди0

■ф*

dw

dw

Еху =

дх

Ф

л

дх

(2,16,22)

 

 

ду

 

~ду~’

d2w

дх2

d2w ду* *

d2w

*ху

дхду

и граничные условия (2,16,21).

Варьируем функционал (2,16,11), используя при этом зависимость между уси­ лиями-моментами и деформациями (2,7,9). Подставив затем в выражение для вариа­

ции значения деформаций из (2,16,22), будем иметь

 

 

 

бэ‘ " И

!г‘* [ 1ST ~

+ Т

С а Г )

]

 

d*w

d*w

— 2Я6 *

d2w

qdw'j dxdy = 0.

(2,16,23)

— Мг6

■М2б „

 

 

дх2

ду2

 

дхду

 

 

 

Преобразуем, как это делали выше,

слагаемые, содержащие

производные

перемеще­

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

граничных условий (2,16,21). Имеем

 

 

6S

дТ2

, в

 

 

 

 

 

 

 

 

 

 

 

 

 

ду

bv0

 

 

 

 

 

 

 

 

 

дх2 J

 

 

 

 

дТг

dS

\

dw

 

Л

ду2

 

дх ^

ду

)

дх

_ f d S

 

дТ2 \

dw

d2w

 

 

 

\ дх

^

ду )

ду

дхду

 

 

 

д2М 1

д2М2

о

д2Н

I .

) . .

 

 

(2,16,24)

дх2

ду2

 

дхду

J

J

 

 

 

 

 

 

Отсюда, приравнивая нулю множители при вариациях био, био и бw, получим уравне­

ния Эйлера для функционала

(2,16,11):

 

 

 

 

 

 

 

 

дТг

dS

 

 

 

 

 

 

 

дх

Ф ду

- ° ’

 

 

 

 

 

 

dS

дТ2

= 0,

 

 

 

 

 

 

дх

ду

 

 

 

 

 

d*w \

I

d2w \

d2w

д2М\

д2М 2

_

д2Н

4 - 9 = 0 .

) - f Tt U + — ) + 2S - ^ 7 4 - - z f - ♦ ~ ^ Г ♦ 2

дхду

дх2

 

ду*)

дхду

дх2

ду2

 

 

(2,16,25)

причем последнее уравнение из (2,16,25) получено с учетом двух первых. Уравнения (2,16,25) совпадают с системой уравнений равновесия для пологой оболочки (2,15,14).

§ 17. НЕКОТОРЫЕ СВЕДЕНИЯ ИЗ ТЕОРИИ УПРУГОСТИ АНИЗОТРОПНЫХ СРЕД

В технике имеют широкое применение материалы со свойствами, зависящими от направления. Такие материалы называются анизотроп­ ными, к ним относятся: дерево, фанера, .прорезиненные ткани, стекло­ пластики и т. д.

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

Приведем основные уравнения теории упругости анизотропного те­ ла. Подробное изложение можно найти в книгах Лява [5] и Лехницкого [16].

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

[ р

'1

(

#11

#12

~ х х

 

 

еу у

 

 

#21

#22

е ZZ

 

 

 

в у г е ZX

#16 #26 7УУ

(2,17,1)

J y z

е х у 1

I1 #61 #62

) Уаху)

Разрешая эти уравнения, получим выражение компонент напряжения через компоненты деформации:

( °хх'

<b11

^12

^16

,е хх ]

УУУ

 

Ь<у

Ь,

еуу

У21 ^22

26

 

 

 

 

N N

Jy z

 

 

 

вуг

 

 

 

 

ezx

К а х у )

<A l

Ь 62

J

Кеху J

(2,17,2)

Обобщенный закон Гука содержит 36 постоянных ац. Если сущест­

вует упругий потенциал, то

 

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

нии обобщенного закона Гука

(2,17,1)

и (2,17,2) снижается до 21. Еще

больше упрощаются соотношения

(2,17,1),

если в строении материала

обнаруживается какая-либо симметрия.

 

 

Рассмотрим основные случаи упругой симметрии.

1.

В каждой точке тела

имеется

плоскость упругой симметрии,

ось z, перпендикулярная плоскости симметрии. В этом случае

 

#14 =

#15 =

#24 = #26 =

#34 =

#35 = #46 = #66 == О

и вместо

(2,17,1)

имеем

 

 

 

 

 

 

в х х =

а 11а хх+

#12° у у +

#13° z z + а и ° х у >

 

 

€ у у ~

#12#**

#22^ у у

#23#zz “Г #26®ху>

 

 

e zz =

#13#^** +

#23О уу +

a 2Sa

zz+ #36ajcу>

^уг ^44^ у г ~Ь ^45^гд:>

(2,17,3)

^ г х ^45^ у г Q 'bfP zx'

е ху — а 1 в ° х+х а и а у у+ а 3 в ° гг+ а ы а х у

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

Закон Гука для ортотропного тела имеет вид

ехх

 

 

еу у

=

а 1 2 °х х

N N

II

согН 3

вуг #44& y z

* z x

^55 ^zx

^ху ^

а1ва ху

+

а 12а уу

+

а13а22»

+

а 22° у у

+

^23^Z2»

+

& 22ру у

+

а 33° Z Z ’

»

 

 

 

У

 

 

 

(2,17,4)

Если в этих формулах выразим постоянные а.ц через технические постоянные (модули упругости и коэффициенты Пуассона), то получим закон Гука для ортотропного тела в форме

ехх

е уу =

N N

II

1

<Угг

V21

 

 

V31

 

-----

Ег

 

a yy

p

a « ’

Е1

хх

 

 

£ 3

 

 

V12

a xx +

1 a

-

V82 0

 

Ег

Ег

 

£a

 

Via

° x x

-2p^ - 0

УУ + - £ ~ a z z ’

 

Ex

 

 

 

 

 

£3

вУг ~ гО з а уг>

(2,17,5)

Р = ----- O'

ху п ^ху>

о12

причем

EIV21= £ 2V12,

(2,17,6)

•^2VS2= ^3V2S> £sV13 — ElV31>

здесь Ей E2, Еъ— модули Юнга в направлениях осей х, у, г, 0 2з, G13, G12— модули сдвига в плоскостях yoz, zox и хоу, Vij — коэффициенты Пуассо­ на, характеризующие поперечное сжатие при растяжении в направлении осей координат. Например, коэффициент vi2 характеризует сокращение

внаправлении у при растяжении в направлении х.

3.В каждой точке тела имеется ось симметрии вращения, иными

словами, все направления, перпендикулярные оси симметрии, эквива-

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

лентны. Плоскость, перпендикулярная оси симметрии, называется плос­ костью изотропии.

Закон Гука в случае, когда ось z совпадает с осью симметрии, принимает вид

*

II

“ i P x x

+

^ i P y y

-f-

^13^22»

 

 

еУУ ~

■& 2 р х х

+

О'Ъ'Рyy

+

a 23C z z »

II N N

° 1 3 а хх 4 - a 2 p y y + Я 3 3 СT2 2 ,

eyz

a iP y z >

 

 

 

 

 

 

а 1 4 р г х *

 

 

 

 

 

II

2 ( f l u - ---

1^ 1 2 ) & xу

 

 

Вводя технические постоянные, перепишем закон Гука так:

 

1

/

. ч

v' .

^хх

g

( ^ JO:

^ У У '

£ / ^*z»

(2,17,7)

еу у —

Р у у

V 0**)

^ZZJ

егг ------------

( a x x Jr

a y y ) +

a zZ’

 

 

ezx

G' JZx>

 

 

 

 

 

 

 

e

= ±

a

- i i i + J i , ,

 

 

(2,17,8)

 

 

ХУ

a

x y

~

P

x y '

 

 

 

где E и Er — модули Юнга

по

направлениям, лежащим

в плоскости

изотропии и перпендикулярным

к

ней, v — коэффициент

Пуассона в

плоскости

изотропии, v' — коэффициент

Пуассона,

характеризующий

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

при

растяжении в направлении,

лежащем

в плоскости изотропии,

G и G' — модули сдвига

в плоскости

изотропии и в плоскости, ей перпендикулярной.

изотропное тело. Закон Гука

4.

Полная

симметрия характеризует

имеет вид

 

 

 

 

 

 

 

 

 

 

 

ехх ~

 

\Рхх

V (р у у

°zz)l>

eyz

=

Gу Z'

 

 

еУУ ~

 

V*yy

V (Рхх “Ь

^zz)]»

e zx =

в 2х>

 

 

^ZZ =

g

[&ZZ

^ Р х х

“Ь

 

^ху ==

&х у >

(2*17,9)

где £ — модуль Юнга, G — модуль сдвига, v — коэффициент Пуассона. В случаях 1—4 рассматривались тела, обладающие прямолинейной

анизотропией.

Теперь рассмотрим тело, обладающее криволинейной анизотропией. Введем систему криволинейных координат g, т], g, таким образом, чтобы

координатные линии совпадали с линиями анизотропии, т. е. с линия­ ми, вдоль которых упругие 'Свойства тела постоянны. В этом случае обобщенный закон Гука запишется в матричной форме так:

л

ег\ц

41

ек

1 е6я J1

1г a x i а 12 . . . а 1в ^ { % \

#21 #23 • • • #26 ^пя

(2,17,Ю)

°я Е

l # e i #62 • • • #66 JI К<%1 1

где ciij — постоянные.

Если координатные линии не совпадают с линиями анизотропии, то величины dij будут функциями координат.

Как и в случае прямолинейной анизотропии, наличие упругого по­ тенциала приводит к равенствам

в результате чего количество упругих постоянных снижается с 36 до 21. Еще больше упрощаются уравнения (2,17,10), если тело обладает упру­ гой симметрией.

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

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

Введем цилиндрическую систему координат, направив ось вдоль оси анизотропии.

Если в каждой точке тела существуют три плоскости упругой сим­ метрии — нормальная к оси анизотропии, радиальная и тангенциальная (случай ортотропного тела с цилиндрической анизотропией), то закон Гука запишется в форме

СЪ ■ч II

1

 

V21

- <т0е

V31

 

а гг

 

E i

Е3

 

 

E i

 

 

 

«ее =

V12

re

1

. 1

V23

-

Еъ

° Г Г

<

г-

^ 0

* “ ■

 

 

 

 

в ,

 

V13

 

 

V23

oee + —- <у22,

 

Ег

 

 

е2

t 3

(2,17,11)

 

1

>

 

 

 

£0Z= ——

 

 

 

 

 

t/23

 

 

 

 

 

1

^31 azr>

1

erQ= —— Ore> Gia

где Eit E2, E3— модули Юнга для растяжения — сжатия в направле­ ниях г, 0, г; V12— коэффициент Пуассона, характеризующий сжатие в •направлении 0 при растяжении в направлении г и т . п., G 2 з, G 3 1» G12 —

5 *

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

Остановимся вкратце на задаче нахождения значений упругих по­ стоянных ац и bij в соотношениях (2,17,1) и (2,17,2) при переходе к но­ вой системе координат. Выше указывалось, что величины а^- и Ьц по­ стоянны только в том случае, когда координатные линии совпадают с линиями анизотропии, т. е. когда упругие свойства материала постоян­ ны вдоль координатных линий. Таким образом, в самом общем случае преобразования криволинейных координат по формулам

Б = Б ( Б . л , » , л ' = л ' ( £ > л , 0 , . л . 06 ' = Е ' ( 6

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

 

 

 

 

 

Т а б л и ц а 2.1

*

2

3

4

5

6

 

 

 

 

 

 

1

®i

«2

«з2

2a2a3

2а3аь

2ct^(Z2

2

p?

P2

Рз

2p2p3

2p3Pi

2piPa

3

Yi

y\

Y3

2Y*Y3

2YSYI

2YiYa

4

P1Y1

PaYa

PsYs

PaYs -Ф- PsYa

PiYa P3Y1

PiYa + P2Y1

5

Yiai

Y2°2

Y3a3

Ya«3 4 Yi°a

Yi°34 Yaai

Yi“a 4 Yaai

6

OlPl

a2p2

ОзР»

ааРз азРг

aiP3,_^ a3Pi

aip2 -£• a2Pi

 

 

Таблица

2.2

В случае прямолинейной

анизо­

 

 

 

 

тропии упругие постоянные a,ij и

 

 

У

 

переходят в постоянные же величины

 

 

 

 

а'ц и Ьц

при ортогональном преобра­

x '

Ol

pi

Yi

зовании

декартовых координат.

Фор­

 

 

 

 

мулы преобразования в общем случае

 

a2

 

 

довольно громоздки, однако им мож­

y'

Pa

Y2

но придать компактный вид, если вве­

 

 

 

 

сти символы cjij — элементы табл. 2.1

2'

a3

Рз

Ys

[16]. Первый индекс

(i) указывает но­

 

 

 

 

мер строки, второй

(/) — номер столб­

ца.

Буквы аа, Ра, уь, обозначают направляющие косинусы преобразова­ ния из табл. 2.2.

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

запишутся так:

6 6

Qlj = ^ '^Q’mnQmiQnj' m = \m= 1

6

6

 

W tW j s s

bmn

(2,17,12)

WmWn Q m iQ n jt

m = 1 п = 1

где

wk = 1 при k = 1, 2, 3, wk = 2 при A = 4, 5, 6.

Вывод соотношений (2,17,12) можно 'найти в курсе Лява [5].

§ 18. ПОСТАНОВКА ЗАДАЧ ТЕОРИИ АНИЗОТРОПНЫХ ОБОЛОЧЕК

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

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

В самом деле, при выводе уравнений неразрывности срединной поверхности (2,3,16), соотношений между усилиями — моментами и на­ пряжениями (2,4,2) —(2,4,5) и уравнений равновесия (2,5,8), (2,5,10) соотношения закона Гука не использовались. Поэтому перечисленные выше уравнения без изменений переходят в теорию анизотропных обо­ лочек. Изменению подвергнутся лишь соотношения между усилиями — моментами и деформациями (2,7,9) и (2,7,10).

Выведем эти соотношения для однослойной ортотропной оболочки. Упругий потенциал ортотропного тела записывается в виде [5]:

W = — (ЬцвЦ + &22е22 + 633633 + ^44е23 + 655631 + 655612) +

623622633 - f - 6з 1е 3д е 11 + ^ 12^ 11^22

( 2 , 1 8 , 1 )

Для оболочки на основании гипотез Кирхгофа—Лява имеем

£33 = *31 = ^23 =

поэтому упругий потенциал ортотропной оболочки в обозначениях § 2.2 имеет 'вид

W = — (ЬцВ^ +&22е2^ +

+ 2612ei ^2 *),

(2,18,2)

а потенциальная энергия

 

 

 

ь_

 

 

2

(6ueiz)2 +

-f 655612* +

 

ABdadfi ^

 

_ h _

 

 

 

2

 

 

+ 2612e<*>e<*>)(l +

- ^ - ) ( l +

 

(2,18,3)

Интеграл слева берется по объему, занимаемому оболочкой, а двойной интеграле правой части — по срединной поверхности.

Вычислим внутренний интеграл правой части

(2,18,3), для чего под­

ставим

в «его вместо е)*1, e(2z\ е$ их

значения

из уравнений (2,3,8),

(2,3,10),

(2,3,14).

 

 

 

Пренебрегая величинами порядка k\h и k2h по сравнению с едини­

цей, получим

 

 

 

 

h_

 

 

 

2

 

 

 

 

J

(&цб(1г)г + . . . ) ( ! +

hz) (1 + м

dz =

 

_

h_

 

 

2

h_

~ f [ 6“ T T 0 (e‘ + *,‘‘), + s , s T T v <' ' + 2’<*) + h_

~ 2

+ 2b12(81~Ь 2Xl) (82+ 2Х2) +

+ ^ee ^ 0 — k±k2z2) е12 + 2 £ 1 f- (fz-L+ k2) — j zx12j' J dz =

=C138i -f- C22e2 -f- 2C12e1e2 -f- Ceeei2

++ D22^2 + 2D12X1X2 + 4D66Xi2,

где k\ и k2— главные кривизны срединной .поверхности,

СИ= НЬф Вц =

ьц.

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

J WdV = — (Cn 8i -f- С22е2 -г ^С12ггг2-f- Ceeei2 +

+ ^ 11^1 + D22^2 + (^D12yi1%2-f- 4DeeXi2) ABdadfi.

(2,18,4)

Вариация потенциальной энергии равна работе внутренних сил и моментов на перемещениях би, 6а, бw. Поэтому из (2,18,4) и (2,7,6) имеем

Jj [(Q liei “Ь Ql282) ^е1 “Ь (^2262 + Ql28l) ^е2

^вв812^12 “Ь

+ (Dii*i + £>12*2) + (А22И2+ Di2Xi) 6х2+ 4/Deex126x12] ABdadfi =

= jj (Т16е1 + Т2Ьг2+ S6e12+

+ Л426х2 + 2Ябх12) ABdadfi. (2,18,5)

Вариации беь бе2 независимы. Поэтому равенство (2,18,5) выпол­ няется тогда и только тогда, когда равны между собой множители при одинаковых вариациях. Отсюда имеем искомую связь между усилия­ ми — моментами и деформациями:

Ti =

Сив! + С12е2,

Мг = Dn xx + D12X2,

Т2=

С2282-f- С12&1»

М2 =

Z)22x2-(- T)12Xi,

S =

Сввв12,

Н =

(2,18,6)

2£>ввх12.