Механика сплошной среды
..pdfз
р2= ^ ' = 1 > г& ,
(4.57).
3
р2- | 2= 2 е г/^ = 2 £ е,&. i=1
Величины gu егсвязанные соотношениями (4.54)
ег = -^-(gi— !) ( / = 1 . 2 . 3 ) . |
(4.58) |
называются главными компонентами тензоров ^ и <§ и они для: точки х = const в момент t =const не зависят от выбора системы координат, т. е. инвариантны относительно преобразований бази са ег. Из (4.54), (4.56) следует, что
Л1=е*/6'7= е 11+ 622+ 633= 6! + е2+ е3,
1е2=■— (ill~ tifiii) = + е2е3 + e3ex, (4.59)
/ез== |^ij I ==
также являются инвариантами ортогональных преобразований ба зиса ег*. В любой фиксированный момент t в точке х = const во всех системах координат инварианты имеют одинаковые числовые значения *.
Итак, в любой точке х тела в любой заданный момент времени t= tx существует главный физический репер-тройка ортогональных
физических волокон, который |
в |
другой |
заданный («начальный») |
момент времени t= t0 состоял |
из |
тех же |
ортогональных волокон. |
Этому факту соответствует полярное представление аффинора де формации А через тензор чистой деформации Е и ортогональный
тензор О, |
__ |
|
|
|
в виде |
E = V W , |
ОтО= 1 |
(4.60') |
|
А=ОЕ, |
Ат = ЕОт. |
(4.60) |
||
|
Существованиепреобразования (4.60), ортогональностьодного и того же физического репера в два фиксированных момента време ни t\, t0 нередко ошибочно трактуются как существование непре рывного процесса вращения физического ортогонального репера
во времени. Физические волокна, образующие ортогональный ре пер в моменты t0 и fi, не ортогональны на интервале t0< t< tu Скорости вращения главных осей деформаций относительно этих
* Явные выражения главных компонент симметричного тензора через его инварианты даны в § 6.
осей не совпадают со -скоростями вращения физических волокон, направленных по этим осям.
Исключение из отмеченного правила непосредственно вытекает из (4.52), (4.55). Главные векторы X и главные оси деформаций не будут вращаться относительно вещества, т. е. физический глав
ный репер будет совпадать с репером главных |
осей деформации |
в любой момент времени при условии |
|
■f(x, t)= k*(t)3'(x), |
(4.61) |
причем gx(x, t) =k2(t)gx'(x) и аффинор Л(х, |
t) = k (l)0 (х, t) X |
ХУ^'(х). Условие (4.61) обеспечивает совпадение с направления ми главных физических волокон не только главных осей дефор маций, но и главных осей любого линейного (см. § 9) оператора по времени над тензором деформации.
§ 5. КОНЕЧНЫЕ, МАЛЫЕ И БЕСКОНЕЧНО МАЛЫЕ ДЕФОРМАЦИИ
Уравнения совместности деформаций. Шесть компонент тензо ра деформаций ец или метрического тензора gij = 6,j + 2cij в ок рестности любой фиксированной физической точки х среды могут как угодно независимо изменяться с течением времени, т. е. зада ние шести произвольных функций времени ец(х, t) возможно, и деформация окрестности точки при этом будет аффинной. Но если бы мы задали для всех точек среды хотя бы в какой-нибудь мо
мент времени t компоненты ец или ga как произвольные |
непре |
|||
рывно дифференцируемые функции |
координат, т. е. произвольно |
|||
задали бы поле тензора деформации, |
то деформации |
оказались |
||
бы несовместными, |
перемещение — неоднозначным, |
т. е. между |
||
соседними частями |
образовались бы щели или различные |
физи |
ческие объемы заняли бы одну и ту же область пространства. Та кая возможность исключена благодаря свойству закона движения х= х (х, t) = х + u (х, t), а именно непрерывной взаимно однознач ной зависимости между х и х для любого t и существованию про изводных. Компоненты тензора ец (или ga) получаются путем дифференцирования вектора х(х, t), т. е. шесть скалярных функ ций ец выражены через три и -. Значит, между ег;- должны суще ствовать соотношения, полная система которых представляет уравнения совместности деформаций. По существу они должны быть следствием независимости порядка дифференцирования век
тора х типа Эщ= э^* , так как ga=9i3j, а векторы |
выражаются |
через один вектор: э*=де,*. |
|
* Для краткости дифференцирование любой функции ф(х, t) по х' обозна чаем ф,( = <?ф/<?Х‘.
Производные векторов репера э* по координатам могут быть также разложены по векторам базиса э* или эг:
|
^ - |
a , l/s=r&.afc= r IMs*f |
(5.1) |
где вследствие симметрии 3i:i= 3id = x tii |
|
||
|
|
rjy= r f f» r w . ^ , |
( 5 . 2 ) |
|
|
Гг/л—Г/J.I= Tijghi\ |
|
|
|
|
|
Г — символы |
Кристоффеля 1-го и 2-го рода; они симметричны по |
||
первым двум |
нижним |
индексам и определяются |
по заданному |
метрическому тензору ga(x, t). В самом деле, дифференцируя по
xh равенство ёц = з^э ^ имеем |
на основании (5.1), |
(4.20') |
|
||||||
|
|
|
|
+ |
|
|
|
|
(5.3) |
Аналогично |
|
|
|
|
|
|
|
|
|
|
|
|
ё/*,г= |
Г/г.А + |
Г*г./, |
|
|
|
|
откуда |
|
|
|
|
|
|
|
|
|
Гг/.А=-^- |
|
£г/,й)==е,*>/"Ь |
|
гИ,к, |
(5.4) |
||||
Гг/.й и Г?/ |
не являются тензорами, что следует |
из |
(5.1). |
получа |
|||||
Математически |
уравнения |
совместности деформаций |
|||||||
ются как условия |
интегрируемости |
системы |
дифференциальных |
||||||
уравнений |
(5.1) |
относительно |
з{, если |
заданы |
функции |
Г*/ (х). |
|||
Дифференцируя |
(5.1) по х1, имеем |
|
|
|
|
|
|||
|
|
|
дх* дх1 |
|
|
|
|
|
(5.5) |
|
|
|
|
|
|
|
|
|
|
Обозначим |
равную нулю при 3{ = Xi, |
разность 3itjl—э,,г/ следую |
|||||||
щим образом: |
|
|
|
|
|
|
|
|
|
|
|
|
д2Э[ |
7 7 7 - ^ |
|
|
(5.6) |
||
|
|
|
дх^дх1 |
|
|
||||
|
|
|
дх1 dxi |
|
|
|
|
|
Из (5.5) видно, что
где значок [/7] означает альтернирование выражения, взятого в скобках, по индексам /, /, т. е.
п» _ дТ'н |
Гп -pm |
дГ™ |
рМ pm |
|
|
il |
(5.7) |
||||
А/7г= |
|
i j l n l ------— ------ 1 U*-nj. |
|||
дх1 |
|
|
dxJ |
|
|
С помощью формул (5.3), |
(5.7) введем |
|
|
||
R m = Rmgmk= |
( ^ |
— |
|
(5.8) |
|
Отсюда очевидно, что справедливы формулы |
|
||||
Rjlih=: |
Rliik== |
Rjlhi== Rihjh |
(5.8') |
||
|
|
|
|
|
|
Rjlih = |
0 при /= / . |
|
|
||
Поэтому можно доказать, |
что среди всех величин |
Rjuh имеется |
только шесть независимых. Для этого следует обратить внимание на порядок расположения четырех индексов выражения Rjuk и учесть, что оно согласно (5.8) равно альтернации по первым двум
индексам выражения, стоящего в скобках (5.8). Поэтому |
тожде |
||
ственно не равны |
нулю согласно |
(5.8') только компоненты Rmnpq, |
|
в которых ни два |
первых (т, п), |
ни два последних (р, q) |
не рав |
ны между собой. Но среди четырех индексов, принимающих каж
дый значения |
(1, 2, 3), |
два |
равны между |
собой. |
Следовательно, |
||||
только следующие две |
группы последовательно |
расположенных |
|||||||
индексов сфра и аРРу, где аф ^ф уф а, принимают значения |
(1, 2, |
||||||||
3) и могут дать независимые и нетождественно |
равные нулю вы |
||||||||
ражения |
(5.8). Первая |
группа дает комбинации |
1221: 2332, |
3113, |
|||||
1331, 3223, 2112, причем либо три первых, либо |
три |
последних |
|||||||
приводят |
к |
независимым |
компонентам, |
например |
# 1221, |
^ 2332, |
|||
Rzm; вторая группа сфРу из шести возможных |
числовых комби |
||||||||
наций 1223, 2331, 3112, 1332, 3221, 2113 дает либо |
первые, |
либо |
|||||||
последние три |
независимые |
компоненты, |
например |
R m 3>-^гззь |
^3112-
УСЛОВИЯ совместности деформаций, являющиеся тождествами,
если метрический тензор 9 или тензор деформации |
|
<§ выражены |
||||
через вектор *(х, t) или вектор u(x, t), имеют вид |
|
|
||||
D __ 5 r aPa |
Пт г |
5 г рра |
, |
г> |
__л |
|
Аарра— |
А арАa pm |
|
г А ррАаа т — |
|
||
|
|
|
|
|
|
(5.9) |
R*№ = |
---- r ^ r Pvm----^ - |
+ r ^ r avm = |
0, |
|||
аф $фуфа, и принимают значения 1 |
2, 3. |
|
|
|
|
Деформацию сплошной среды в эйлеровом пространстве х за бесконечно малое время dt в любой фиксированный момент t= t0 можно рассматривать с точки зрения Лагранжа, если поле векто ра скорости v(x, t) задано и если в момент времени t'=to + dt определить перемещение
и' (х, |
t')=v(x> |
V)dt |
(5.10) |
и координаты частиц |
|
|
|
х ’^ х |
+ и' (л:, |
V). |
(5.11) |
В пределах интервала dt при данном t{xu *2, *3) будет лагранжевой ортогональной системой координат в репере ег. Тензор беско нечно малых деформаций среды за время dt обозначим v^dt=z'i^ причем Vij называется тензором скоростей деформаций среды в
эйлеровом |
пространстве. |
находится по формулам (4.49) путем |
Тензор |
деформаций е'г?- |
|
замены х-+х\ х->х', |
U- H I', |
причем нелинейные слагаемые суть бесконечно малые высшего порядка. Метрический тензор ga находится из (4.48) и равен
(5.12')
т. е. при dt-+0 он равен б*j. Внося в (5.12) выражения и' из (5.10), деля на dt и отбрасывая бесконечно малые величины, получим выражения компонент тензора скоростей деформаций в эйлеро вом пространстве
(5.13)
С другой стороны, можно рассматривать малые относительные перемещения точек среды, т. е. такие, что в любой момент време ни t>to {to — начальный момент) перемещение точек u(x, t) = =х—х обладает свойством
- ^ - < 6 , 6 < 1 |
(5.14) |
dxj |
|
для всех i>/ и в любой точке х. В этом случае можно |
исключить |
переносное движение среды (связывая систему координат х с ка
ким-нибудь физическим |
волокном и плоскостью |
в точке х=0); |
из (5.14) получим, что |
в подвижной системе сам |
вектор переме |
щения и будет малым порядка б сравнительно с размерами обла сти, занятой средой. Разница между этим и рассмотренным выше случаем в том, что в случае (5.10) dui//dxj = dvildxjdt является бесконечно малой величиной.
От (5.14) отличается случай малой деформации, когда eij со храняет полное выражение (4.49), но удовлетворяет условию
Ы < б , б<1, |
(5.15) |
отличие в том, что при малых деформациях система координат может быть не привязана к физическим волокну и плоскости.
Мы рассмотрим случай малых порядков б перемещений, затем автоматически распространим все результаты на бесконечно ма лые деформации. Неточность всех формул теории малых дефор маций будет порядка б сравнительно с единицей. Соответствую щие формулы в случае бесконечно малых деформаций будут точ ными. Вектор относительного перемещения 8 точки N относитель но М в момент t, очевидно, равен
6 = P - S = * ( x + S, t ) - x ( x , |
дх |
1. |
|
|
|
|
|
Обозначая D=A—I тензор дисторсии, получим |
согласно |
(4.4), |
|
(4.7) |
|
|
|
6 = D |= ~Ui— |
|
|
|
OX |
OXj |
|
(5.16) |
|
|
|
|
Пренебрегая малыми порядка |
б2 сравнительно |
с б, из |
(4.48), |
(4.49) получим для компонент метрического тензора ^ и тензора деформаций 6 выражения компонент
%= (gii), |
gu= Sij + 2et/, |
|
(5.17) |
1= (е„), е17= е /г= ^ |
( ^ - + -^ -) ( / , / = 1 , 2 , 3 ) . |
Соотношения (5.17) представляют формулы Коши, выражаю щие тензор малой деформации через вектор перемещения и. Со ответствующие выражения компонент тензора скоростей дефор мации Vij в эйлеровом пространстве, как видим, получаются из (5.17) простой заменой и на v и х на х и приводят к формулам Стокса:
Vil= Vii |
- L |
( дщ |
dv/ ) |
(5.18) |
|
2 |
U x , + |
дх{ ) |
|
Малые и бесконечно малые деформации аддитивны в том смысле, что если даны поля перемещения u1(х, t) и u2(x, t) с со ответствующими деформациями е{., г\}, вычисляемыми по фор
мулам Коши (5.17), то полю перемещения u = u1-fu2 соответ ствуют деформации Sij, равные сумме соответствующих деформа ций:
и(х, /) = их(х, 0 + “2(х, 0. 8„=в‘/ + в*. |
(5.19) |
Аналогичное верно для двух полей эйлеровом пространстве при любых формул Стокса
скоростей vl(x, t), v2(x, t) в деформациях на основании
\ { х , /)=У х(л;, f) + va(*. /), Vil= v\. + v2ir |
(5.20) |
Как уже отмечалось в § 4, можно установить кинематический смысл компонент тензора малых деформаций. Из формулы (4.38) находим с рассматриваемой точностью
(!+ е а)2= 1 +2еаа; еа= еаа.
Косинус угла 0ар между координатными волокнами а и р , равный нулю до деформации, после деформации согласно формуле (4.44) равен
cos 0ар=2еор (1—еаа—ерр)=2бар,
т. е.
— 0<xp=sin ^ — 0ар) = 2 еар-
Следовательно, компоненты тензора малой деформации с одина ковыми индексами суть относительные удлинения координатных
волокон, а удвоенные компоненты малой деформации |
со смешан |
|||||
ными индексами суть уменьшения прямых углов |
между |
парами |
||||
координатных волокон, называемые сдвигами. |
|
|
малой |
|||
Координатный вектор-волокно |
(£)а = £аеа в результате |
|||||
деформации станет вектор-волокном (р)а, направляющие |
косину |
|||||
сы которого согласно (4.46) определяются формулами |
|
|
||||
|
|
|
|
|
|
(5.21) |
(P)a=l(P)al4eI= |
| a (6to + - g - ) e ,, |
a = l , |
2, |
3. |
|
|
Разложение тензора дисторсии на тензор чистой деформации и |
||||||
тензор поворота. Рассмотрим тождество |
|
|
|
|
||
D = g + a>, |
co=(©i/); |
|
|
|
(5.22) |
|
дщ-= ег/+ со17) |
|
1 ( |
ди |
дио |
\ |
|
|
|
|||||
coap= -copa= T ^ - 3 - — ^й - j. |
|
|||||
дх; |
|
|
|
|
|
|
в котором со называется тензором поворота. Тензор со антисимме
тричен, и потому он может быть выражен через вектор |
о, назы |
||
ваемый ротором вектора и, или вектором |
поворота окрестности |
||
точки х, |
|
|
|
о)= -j- rot 11= © ^ . |
|
(5.23) |
|
Компоненты вектора со связаны с компонентами тензора со со |
|||
отношениями |
|
|
|
1 ( диу |
ди^ |
\ |
(5.24) |
|
|
|
|
На основании разложения (5.22) вектор |
относительного пере |
||
мещения 6 (5.16) представим в виде суммы |
|
|
|
6 = 6 е + 6“, 68=6fel= e //^/el., |
60)=col7g/ei= coxS, |
(5.25) |
причем 6е и 6“ называются соответственно перемещением, связан ным с чистой деформацией, и перемещением, связанным с поворо том со окрестности точки М. Равенство 6Швекторному произведе
нию ©XI проверяется |
непосредственно |
на |
основании |
(5.24) и |
||
определения |
|
|
|
|
|
|
|
|
|
е3 |
|
|
|
|
© х |= ©1 ©2 ©з |
|
|
(5.26) |
||
|
|
£l |
|
|
|
|
Метрический тензор с точностью до б2 равен |
|
|
||||
£ = / + 2?, |
g = -L (D + DT) |
|
(5.27) |
|||
и, следовательно, в соответствии (4.59) |
|
|
|
|||
1+ 2eu |
2e12 |
2e13 |
|
|
|
|
g = \g i/\ = |
2e2i |
1 -f- 2e22 |
2e23 |
= |
l + 2 / ei. |
(5.28') |
|
2e3i |
2^32 |
1 + 2e33 |
|
|
Поэтому относительное объемное расширение частицы будет рав но
~ — V g —1 = /e i= e 11 + 822 + e83= e <i= divu . (5.28)
Закон сохранения массы р(1 + 6 )= р 0 дает выражение плотности р через 0 и начальную плотность р0:
В частности, если можно пренебречь изменением плотности, то условие объемной несжимаемости вещества имеет вид
0 = |
+ ^ |
+ Jfc3_= divu==0. |
(5.29) |
дхх |
дх2 |
дх3 |
|
Разложение тензора малой деформации на девиатор и шаро вой тензор. Девиатором деформаций называется тензор с компо нентами Eij,
S i j = E i j — ебг-j, е = 1/30, |
(5.30) |
причем ебг-j называется шаровым тензором, так как соответствую щая ему поверхность есть сфера
ебt/X tXf = е (Х\ + Х\ + Х\) = const.
Первый |
инвариант девиатора |
деформаций |
ЛГ = е^--6^ = 0, вто |
||
рой — |
|
|
|
|
|
|
—2/2i = эг= е г/ег/==е2 +1* + Ц, |
(5.31) |
|||
где еь е2 |
ез — главные компоненты девиатора |
деформаций. |
|
||
Относительное перемещение 6е |
(5.25) |
в главных, осях деформа |
|||
ций имеет |
выражение (в главном |
репере |
ег0) |
6е = ег-£гего. Рассмат |
ривая октаэдрическую площадку (отсекающую равные отрезки по
главным |
осям) и волокно |
равно наклоненное к главным осям, |
так что |
li= \\\lY 3, найдем |
его относительное удлинение |
еЪ= &!/£а= 1/3 (е1"Ьб2+ е3)= е=0/3.
Сдвиг конца волокна % (линейный сдвиг октаэдрической площад
ки) равен V(6e)2—| 2е|, относительный сдвиг получается деле нием на \%\. Удвоенный относительный сдвиг октаэдрической пло щадки называется октаэдрическим сдвигом yv:
yv= 2 У(бЕ)2- 1 2е|/1 || = 2 Kl/3(e? + е^ + е з )-е а7
На основании |
( 5 . 3 1 |
)получаем |
|
|
|
|
|
|
= Y |
v = |
2 /V e3i + Ц +~sl= |
2 / 3 V |
( e x - |
e 2) 2 + |
( e a- |
e 3) a + (e |
|
|
|
|
|
|
|
|
( 5 . 3 2 ) |
|
T . e. пропо р ц и о н а л е н |
к о р н ю из |
вто р о го |
и н в а р и а н та д евиатора |
де- |
||||
формации |
с обратным знаком. |
Величина э-Л /г^ри |
называется |
|||||
модулем |
девиатора, |
или интенсивностью деформации. |
Аналогии- |
ные формулы верны для девиатора скоростей деформаций (Vk— главные компоненты)
Vij= va— 1/Зи/а6</, |
1/3 [(«!—t£)a + {щ— vzf + |
|
(5.32') |
Условия совместности Сен-Венана для малых деформаций по
лучаются из (5.4), (5.9). |
Символы Кристоффеля |
(учитывая ga= |
|||
= 6,j + 2e,j) будут малыми |
порядка б, а их произведения — поряд |
||||
ка б2. В результате получаем шесть уравнений |
|
||||
дГдр» _ ^грр« |
_££ору__ _f£ggv_ |
(5.33') |
|||
дхе |
д\ ' |
дхе |
дха |
||
’ |
из которых следуют шесть условий Сен-Венана:
2 |
|
__ _32еи_ |
дае22 |
|
|
|
|
Зх3 дх2 |
дх\ |
Зх2 |
|
|
|
32elt |
д / |
dea3 | |
Зе81 ^ |
де1г \ |
(5.33) |
|
Зх3 Зх2 |
дх1 \ |
dxi |
Зх2 |
Зх3 ) |
||
|
(в каждой из выписанных формул следует сделать круговую пе рестановку индексов 1, 2, 3). Естественно, что если в (5.33) вне сти выражения ец по формулам Коши (5.17), получатся тождест ва. Иногда их называют тождествами Сен-Венана.
Все полученные выше для малых деформаций среды формулы при указанной ранее замене х-*-х, u->-v, ei3-*-Vij справедливы для мгновенных деформаций ее в эйлеровом пространстве, но имеют соответствующую трактовку. Заменяя и на v(x, t), е - - на viit х на х, из (5.23) и (5.24) находим
й(л:, /)= -^ -го 1 у = й гег, |
(5.34) |
причем S2 называется вектором вихря — вектором |
мгновенной уг |
ловой скорости вращения частицы среды. Тензор |
(5.18) пред |
ставляет тензор скоростей деформаций частицы, т. е. компоненты
vn, v22 V33 суть |
мгновенные скорости |
относительных |
удлинений |
|
координатных волокон, взятых по осям х,, |
а удвоенные смешан |
|||
ные компоненты |
2u,2l 2v23 2v3, — скорости |
сдвигов (скольжения) |
||
координатных площадок. (Их механический смысл |
понятен из |
|||
рис. 5.1.) |
малых деформаций |
соответствуют |
инварианты |
|
Инвариантам |
тензора скоростей деформаций
Л = 1>х1+ v22 + t'33= y<A/=div v,