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

Механика сплошной среды

..pdf
Скачиваний:
33
Добавлен:
15.11.2022
Размер:
13.71 Mб
Скачать

з

р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,