книги / Механика сплошной среды. Т.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) в базисе га имеет еле-