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

книги / Механика сплошной среды. Т.4 Основы механики твёрдых сред

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

§ 2.9. Осесимметричные задачи теории упругости

161

тогда из (2.8.106) и (2.8.107) находим напряжения и перемещения в круговом цилиндре при кручении:

2M ix2

2Mix'

2м Д 2ж3

2М Д Ч 3

сг13= ---------- 0"2 3 =

--------------------- « 1 =

--------------------ад2 =

----------------------« 3 = 0 -

а

a

7rGа

т г и а

 

 

 

( 2.8. 111)

Для приложений важное значение имеет интенсивность тензора напряже­ ний (2.5.65), которую в данном случае находим по формуле

\ / 2(сг?3 +

з)

 

Ш 1 г

4

13 1

23'

а

Максимальное значение интенсивности достигается на внешней поверх­ ности цилиндра при г = а и составляет crwmax = 4М^/(тга3) — это значение пропорционально крутящему моменту МI и обратно пропорционально третьей степени радиуса цилиндра. □

§ 2 .9 . О с е с и м м е т р и ч н ы е за д а ч и т е о р и и у п р у г о с т и

2.9.1. Общие уравнения

Некоторые твердые тела обладают цилиндрической анизотропией — для них главный базис анизотропии са (х) совпадает не с ёа, а с физическим базисом цилиндрической системы координат er , е^, ez, и поэтому различен при переходе от одной точки тела к другой. Цилиндрической анизотропи­ ей обладают многие композиционные материалы, образованные намоткой, а также конструктивно-анизотропные материалы (подкрепленные металлокон­ струкции, вафельные оболочки и др.).

Пусть для линейно-упругого тела выполнены следующие условия:

1) соответствующая ему область V является телом вращения

(см.

т. 1,

п. 4.2.6), т. е. для него существует ось вращения, которую

будем

по­

лагать совпадающей с осью Ох3 и осью Oz цилиндрической системы координат Огфг;

2)тело является цилиндрически ортотропным во введенной цилиндриче­ ской системе координат;

3)векторы плотности массовых сил f' в V, поверхностных усилий \!пе на

части поверхности

перемещений ие

на части поверхности

тела

ортогональны к е^:

 

 

 

 

 

? - е ф = 0, t'ne

= 0,

ие

= 0,

(2.9.1)

иих компоненты в базисе ег,еф,ег не зависят от ф, а зависят только от

гя z :

f'

f'

t'

t'

UP

Г, 2.

(2.9.2)

Jr

Jz>

vner

vru

 

 

 

Тогда решение квазистатической задачи линейной теории упругости (2.6.63) также будет осесимметричным, т. е. компонента иф вектора перемеще­

162

Глава 2. Упругие среды с малыми деформациями

 

ний и будет нулевой, а две другие будут зависеть только от г и г:

 

ur, uz || r,z\

иф = О в V U £ .

(2.9.3)

 

Тензор деформации е, вычисляемый по соотношениям Коши (2.1.20),

записанным в цилиндрической системе координат (см. т. 1, упр.

11 к § 2.6),

имеет следующие ненулевые компоненты, также зависящие только от г и г:

 

_

диг

_ u r

_

duz

 

_

1 ( диг

duz\

£r

dr

г

8z

dz

8rz

2 Vdz

dr )

В силу ортотропии среды, тензор модулей упругости 4С в базисе ci = ег,

С2 = еф,

сз = ez имеет вид (2.6.13), а ненулевые определяющие соотношения

(2.6.2)

(тепловые деформации включены в f') в цилиндрических координатах:

 

 

<Jr — Xi£r + А6£ф + A$£z,

СГф— Аб£г + А2£ф + A\£z,

 

 

 

&z — А$£г + А4£ф + A^Sz,

arz = 2 A%£rz.

(2.9.5)

Здесь а г =

а г г , Оф

= Офф, a z =

crz z ,

аналогичные обозначения введены и для

компонент тензора деформации.

 

 

 

 

 

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

ские константы (2.6.16) и (2.6.17):

 

 

 

 

 

 

 

 

_

1

Щ2

 

Щз

_

 

Щ2

. 1

^23

 

 

Е[ аг

Е[ аф

Е\ a z ’

ф~

 

 

Еуаг+ Е2° ф

Е2 a z ’

 

. _

^13

 

^23

.

1

£ rz

 

1

(Jnr

 

 

(2.9.6)

z

^

 

^

 

ЕЛ

 

 

VEa

Ep

 

 

 

е 2

 

 

2G13

 

Согласно (2.9.3)-(2.9.5), все выписанные компоненты тензоров напряже­

ний и деформации зависят только от г и г, и

 

&гф £z<fi —0, &гф оzcf) —0 в V.

(2.9.7)

Уравнения равновесия в задаче (2.6.63), записанные относительно напря­ жений (форма (2.2.9а)), при переходе к цилиндрической системе координат

(см. т. 1,

упр.

12 к

§ 2.6)

содержат

только

два

тождественно ненулевых

уравнения:

 

 

 

 

 

 

 

дог +

+

+ о

, = о

 

 

о ,

or

г

 

oz

 

or

Г

OZ

Граничные условия из (2.6.63) с учетом (2.9.2) в осесимметричном случае

принимают вид

 

 

(oynr + arznz) |2ст = t'ner, (arznr + aznz) |Е^ = t'nez,

(2.9.9а)

= ^ег*

=

(2.9.96)

§ 2.9. Осесимметричные задачи теории упругости

163

Если ось вращения O z принадлежит телу В , то на части границы двумер­ ной области V' определения задачи должны быть заданы условия симметрии

г = 0: arz = 0, иг = 0.

(2.9.10)

Возможен частный случай, когда на боковой поверхности £ заданы толь­ ко поверхностные усилия t fne (при этом, как известно, необходимо задать еще перемещение тела как жесткого целого в одной точке), тогда следует использовать уравнения совместности деформаций (2.3.7), которые в цилин­ дрических координатах запишем с помощью выражений из т. 1, упр. 13 к § 2.6; в осесимметричном случае эти уравнения представляют собой только

четыре тождественно ненулевых соотношения:

 

pc e rz = 0,

 

д 2е г

d 2s z

2 ^ ^rz

= 0,

д2£ф

дг^ _

 

 

 

 

 

 

 

 

 

~ д! +

drz

d r d z

 

Г dz2

dr

 

dz

 

I д

/

2 9 6ft.

dsr

= 0,

d_(

_ djrstt,)

0.

(2.9.11)

-----

—-

d r d r

dz^ r

dr

 

г d r

 

d r

 

 

 

 

 

Дифференцируя второе из этих уравнений по г и подставляя его в первое

уравнение, получаем

 

 

d

 

 

 

 

 

 

 

з ? (г'

~

 

 

= 0

Приведем третье уравнение в (2.9.11) к виду

1 (<?

_

д (Г£^ \ -

0.

г

\ЬГ

 

О

=

 

 

d r

 

 

Тогда совместно с четвертым уравнением в (2.9.11) после интегрирования

имеем

ч

ег ---- =

д, д = const,

d r

 

причем константу d всегда можно принять равной нулю: д = 0.

В результате система (2.9.11) оказывается эквивалентной только двум

уравнениям совместности:

 

 

 

Щ

+ г ^

= 2d e rz

de

(2.9.12)

d z

d z 1

d z

Г ~dr

£ Ф ~ Sr ~

2.9.2. Постановка осесимметричной задачи в перемещениях

Подставляя выражения (2.9.4) и (2.9.5) в уравнения равновесия (2.9.7), получаем систему двух уравнений второго порядка относительно иг и uz\

Х \ и г ,г г

+ А%иг ,zz

+ (А5 А8 ) и z , r z

+ Ai ( и г ' г / r ) -

 

<

 

~ а2(«ГД2) + (А5 - A4 ) ( u z ,z / r ) + p f r = 0.

(2.9.13)

As u z , r r

+ A3U z ,zz

+ (А5 + A§ ) u r , r z

+ (А4 + As ) ( u r , z / r )^~

 

k

 

 

+ ^ 8 ( u z , r / r ) + P f ' z = 0-

 

Здесь ur>r = dur/dr, ur>z = dur/dz.

164

Глава 2. Упругие среды с малыми деформациями

Подставляя (2.9.4) и (2.9.5) в граничные условия (2.9.9) и (2.9.10), выра­ зим их через перемещения:

(пг{\\иг,г + А§иг/ г + A§uZjZ) + n z \% {ur ,z + ^ z ,r))^ = tneri

(jlz(^\Ur,r H- A6^r/r + A§UZ^Z) + T i r \ % ( u r

z + иг^У)

^nez’

(2.9.14)

I ZJU

 

I ZJU

 

u ez'i

 

 

u r \ y

 

— ^ e r ?

u z \ y

 

 

 

r = 0:

urjZ + uZjr = 0,

 

 

ur = 0.

 

(2.9.15)

Задачу (2.9.13)—(2.9.15) называют осесимметричной задачей в переме­ щениях.

2.9.3. Функции напряжения в осесимметричной задаче

Если массовые силы потенциальны, т. е. существует такая скалярная функция x(r,z), что

fr = X,r,

f'z = x,z,

(2.9.16)

то можно ввести две функции напряжения Ф(r,z) и y(r,z), удовлетворяющие соотношениям

<Jr = p ~ X , о-ф = Ф,г - r^tZZ - х, crz =

- X, ®rz = 1 ,z (2 .9 .1 7 )

(Ф и 7 никак не связаны с функциями (2.8.15) для двумерных задач). Непосредственной проверкой несложно убедиться, что выражения (2.9.17)

с учетом (2.9.16) тождественно удовлетворяют уравнениям равновесия (2.9.7). Подставляя (2.9.6) в (2.9.12), выражаем уравнения совместности через

напряжения:

j —1

 

/ 6, / 6

" О ! ~ 7 ,/6

"

С/7,/6

1

771

 

V "

,

 

" Z I ~ 7

/ 6/ 6

 

"

/6 ,/6 /6 /

п

 

(J г

 

= 0,

тр

(^ 2 ,2

^31 a r,z

1/32СГф^))

+

тр

 

{р'ф^гу

 

V2\<Jr,zz

 

^23a z,zz)

 

^ TZ,TZ

 

Ц/3

 

 

 

 

 

 

 

Ц/2

 

 

 

 

 

 

 

 

Ст

 

 

 

Г

 

/

 

 

 

Ч

 

1+^12

,

1+^21

 

Ц32

Ц31

чз

 

 

 

 

^21

г

 

 

 

 

 

(2.9.18)

тр

 

(&у,Гф,г

^23" Z O&z,r)~ Z 7 /

 

 

тр

~<Jr+

тр

®ф

"

тр

-az =0.

Ц/2

 

" Z ^ 7 , 7

 

 

 

 

 

Ц/1

~7-

'

Ц/2

 

-с/3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь использованы соотношения взаимности из (2.9.6) для коэффициентов Пуассона.

Подставляя выражения (2.9.17) в (2.9.18), получаем систему двух уравне­ ний относительно двух функций напряжений Ф и у, а (2.9.17) в (2.9.9а) — граничные условия для этих уравнений на поверхности £ двумерной области. В результате получим формулировку осесимметричной задачи в напряже­ ниях.

2.9.4. Задача Ламе

Рассмотрим случай, когда осесимметричное тело представляет собой по­ лый цилиндр с внешним и внутренним радиусами гд и г\ (рис. 2.9.1), на боковых поверхностях которого заданы равномерные давления ре2 и ре\,

§2.10. Оболочки и пластины

 

167

2 r fe+1

1 4 - r 2k

(2.9.29)

f = f 2: Оф(2)=Ре\к

1_2fc —Ре2&~щ:-

1

~ г\

1 -Г\

 

Рассмотрим предельный случай очень тонкого цилиндра, для которого h = = h/r2 <С 1, где h = Г2 — г\, тогда пренебрегая малыми по величине членами,

имеем 1 —f \ k = 1 (1 h)2k « 2kh,

1 + f \ k « 2,

и из (2.9.29)

получаем

°>(1) ~ 0-0(2) ~ (Pel ~Pe2)R/h,

h = r2 -

г ь Д = г2.

(2.9.30)

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

Поскольку максимальное значение радиальных напряжений в такой обо­ лочке max |<тг| = max{|pei|, |_ре21}, то при конечных значениях перепада дав­ лений (ре1—ре2) и малой толщине h значения тангенциальных напряжений &ф(\) оказываются существенно больше максимальных значений радиальных напряжений max |оу|. Таким образом, эти напряжения являются наиболее опасными при оценке возможной потери прочности оболочки (при &ф(\) > 0) или потери ее устойчивости (при сгф(\) < 0).

§ 2.10. Оболочки и пластины

2.10.1. Понятия оболочки и пластины

Одними из наиболее широко применяемых на практике типов конструкций являются оболочки и пластины — тела, у которых характерный размер в одном из направлений h существенно меньше двух других характерных размеров lj: h/li <С 1. Обычно полагают, что модели теории оболочек доста­ точно адекватны при соотношениях h/li ^ 1/20. Частным случаем оболочек можно считать и стержни — тела, у которых два характерных размера h\ и h2 намного меньше третьего характерного размера 1: hi/l 1. Оболочки, пластины (плоские оболочки) и стержни из различных материалов являются основой современных строительных конструкций: зданий, мостов, корпусов автомобилей, самолетов и кораблей, ракет и спутников. Уравнения теории оболочек широко используют при проектировании таких конструкций.

Дадим следующее геометрическое определение оболочки.

Определение 2.10.1.

О б о л о ч к о й в теории твердых сред называют

твердое тело V с 8 %,

если для него можно ввести ортогональные (вообще

говоря, криволинейные) координаты Х \ в которых это тело представля­ ет собой окрестность Vh двумерной поверхности HQ (см. т. 1, п. 3.2.19).

Если для тела V в качестве указанных координат Х г можно выбрать прямоугольные декартовы координаты х \ то такое тело называют пластиной.

Понятие окрестности поверхности Vh предполагает, что в выбранных координатах Х г две координаты X 1 являются координатами некоторой спе­ циальной поверхности SQ, называемой срединной поверхностью оболочки,

170

Глава 2. Упругие среды с малыми деформациями

2.10.2. Основные уравнения линейной теории упругости в ортогональных координатах

При построении моделей оболочек в качестве исходной рассмотрим квазистатическую задачу линейной теории упругости в форме (2.2.9):

' V- <r + pf = 0

в V,

(2.10.10а)

сг =

С • • е B F U E ,

(2.10.106)

< е =

(l/2)(V<8>u + V ® u T) B F U E ,

(2.10.10B)

n • сг|Ест = t ne,

u |Eu = ue.

(2.10.Юг)

Запишем эти уравнения в ортогональных координатах Х г, введенных выше для оболочек.

Используя выражение (т. 1, (2.6.41)) для дивергенции тензора второго ранга V • сг в ортогональных координатах, уравнение равновесия (2.10.10а) запишем в следующем виде:

{ H ^ H j ( j a a ^ a + ( Н а Щ с г а р ) ' Р + ( Н а Н р ( г а гу')

+ (7арНгуНар + (ia j H p H a j

- &ррЩНрасг1 1 НрН1а+ Н\Н2 Щр$а= 0,

а, /3,7 = 1,2, 3,

а ф /3 ф у ф а.

Здесь и далее обозначены производные по Х а :

(2. 10. 11)

 

(тарп = даа(3/ д Х \ На 1

= ЭНа/ д Х \

(2.10.12)

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

Соотношения Коши (2.10.10в) в координатах X 1 записываем, используя формулы (т. 1, (2.6.37)) для линейного тензора деформации над вектором перемещений и:

^

_ ^а,а

! Яа(3

 

 

 

 

в ГУ ГУ

-г-г-

“г

 

Z \ f l f i ' - П - а '

^ а ' T i p '

<21013>

 

Н а

Н а Н р

П а П j

В формулах (2.10.11)—(2.10.13) и далее обозначены физические компо­ ненты тензоров сг и в, а также вектора и в ортонормированном физическом

базисе га :

га = га/Н а,

а =1, 2, 3;

 

 

 

 

 

з

з

з

 

 

СГ= ^ 2 ^а/ЗГа^У/З,

£ = ^ £а13Га®Г/3,

U = ^ Щ аГа .

(2.10.14)

а,/3=1

a,f3 = 1

 

а= \

 

Предположим, что рассматриваемая

оболочка

является ортотропной, а

ее главные оси ортотропии Оса совпадают с осями Ога построенной для оболочки системы координат (векторы г/ = с/ являются касательными к ли­ ниям главных кривизн срединной поверхности оболочки), тогда компонентное представление определяющих соотношений (2.10.106) в базисе га имеет еле-