книги / Механика сплошной среды. Т.4 Основы механики твёрдых сред
.pdf§ 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 -  | 
7га  | 
	7гa  | 
	7rGа  | 
	т г и а  | 
  | 
	
  | 
	
  | 
	( 2.8. 111)  | 
Для приложений важное значение имеет интенсивность тензора напряже ний (2.5.65), которую в данном случае находим по формуле
— \ / 2(сг?3 +  | 
	з)  | 
	
  | 
	Ш 1 г  | 
|
—  | 
	4  | 
|||
13 1  | 
	23'  | 
	—  | 
	7га  | 
|
Максимальное значение интенсивности достигается на внешней поверх ности цилиндра при г = а и составляет 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) в базисе га имеет еле-
