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

книги / Механика сплошной среды. Законы сохранения

.pdf
Скачиваний:
1
Добавлен:
12.11.2023
Размер:
3.31 Mб
Скачать

Исследователи по разным причинам не работают с какой-либо одной мерой деформации. От одной меры всегда можно перейти к другой, так что все меры одинаково ценны, и нет оснований какой-либо одной присваивать наименование «истинная». На рис. 1.4 для качественного сравнения приведены зависимости различных мер деформации от кратности при

одноосном удлинении среды: 1 — меры Грина 21 ( 2 1) , 2 — левой меры относительных удлинений 1, 3 — логарифмической меры ln , 4 — правой меры относительных удлинений 1 1 , 5 — меры Альманси 21 (1 12 ) .

Рис. 1.4. Зависимости мер деформаций от удлинения

Следует заметить, что все меры деформации неодинаково ведут себя при растяжении и сжатии. Например, при деформировании стержня со

скоростью

деформации

/ lˆ / lˆ 1

с-1 (знак плюс

соответствует

растяжению, минус — сжатию) зависимость удлинения

от параметра

нагружения

exp( t)

изображена на

рис. 1.5 (кривая 1 соответствует

растяжению, кривая 2 — сжатию).

21

Рис. 1.5. Зависимость удлинения при растяжении и сжатии

1.7. Изменение материальной площадки и объема

Рассмотрим, как изменяется малая материальная площадка при ее аффинной деформации. Возьмем в K0 пару произвольных малых материальных отрезков, dx и dx , исходящих из некоторой материальной точки, и рассмотрим построенный на этих отрезках малый материальный

параллелограмм.

 

Если

 

 

площадь

 

 

параллелограмма

 

 

обозначить

dS , а

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Kt

мы получим

единичную нормаль к его плоскости n, то ndS dx

dx

. В

материальный

 

 

параллелограмм

 

 

 

 

ˆ

 

 

ˆ

 

ˆ

 

 

ˆ

,

 

 

 

где

 

 

 

 

F dx ,

 

 

 

 

 

 

ndS dx

dx

 

 

 

 

 

 

dx

 

 

 

 

ˆ

 

и

 

 

— его площадь и единичная нормаль. Свяжем

dx F dx , а

dS

 

n

 

 

 

ˆ

 

ˆ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

векторы ndS и ndS :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

nˆd ˆ dxˆ dxˆ dx F T F dx dx

ei e

e

e j dx

 

 

 

 

 

 

 

ˆ ˆ

 

 

 

 

 

 

 

 

 

 

ˆ

 

 

 

 

 

 

 

 

i

 

 

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ˆ T

 

 

 

 

 

 

d

i

d

 

 

j

ei

e j

 

 

 

 

Jd

i

d

 

 

j

 

ijk e

k

 

 

gdˆ

i

d

j

 

 

e

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ijk F

 

 

 

 

 

 

gˆ d d F

 

( ijk e )

 

 

 

 

d d F

(ei

e j )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

j

ˆ T

 

 

 

 

k

 

 

 

 

 

ˆ

 

i

 

j

ˆ T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ˆ T

 

 

 

 

 

 

 

 

 

 

ˆ T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ˆ

 

F

 

(dx dx )

JF nd ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g

 

 

 

 

 

 

 

 

 

 

где g ||gij ||,

g ||gij || — квадраты объемов параллелепипедов, построенных

на локальных материальных базисах рассматриваемой материальной точки в

K

и

K

,

J

 

g

— якобиан преобразования

 

материальных

координат.

 

 

 

0

 

t

 

 

 

g

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Короче говоря,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

1

T

 

 

 

 

 

 

 

 

 

 

nd ,

nd J F

 

(1.55)

 

 

 

 

 

 

 

nd JF

 

nd .

 

 

 

 

 

 

 

 

 

 

22

 

 

 

 

Перемножая любое из соотношений (1.55) скалярно само на себя и используя

(1.29), (1.30), получим

d ˆ / d J (n U 2

n)1/2 ,

d / d ˆ J 1(nˆ V 2

nˆ)1/2 .

(1.56)

Далее, подставляя (1.56) в (1.55), можно связать нормали:

 

 

nˆ (n U

2

n)

1/2

ˆ T

n ,

n (nˆ V

2

nˆ)

1/2

F

T

nˆ .

(1.57)

 

 

F

 

 

 

Наконец, следует установить связь локальных материальных объемов в K0 и Kt . Проще всего для этого рассмотреть малый материальный шар

радиуса dr в K0 , который, как нам известно, превращается в Kt

в малый

материальный эллипсоид с полуосями 1dr, 2dr и 3dr . Отсюда

 

 

dV

 

 

 

 

 

.

(1.58)

 

 

2

3

 

dV

1

 

 

 

 

 

 

 

 

 

 

 

 

Эта величина называется кратностью изменения объема. Из (1.58) и

определений (1.36), (1.37), (1.39) следует

dV

||F|| ||U|| ||V|| J .

(1.59)

dV

 

 

Можно выразить данную величину и через меры деформаций (1.46), (1.47), (1.50), (1.51):

 

ˆ

 

 

 

 

 

 

 

 

 

 

 

dV

(1

2I (G) 4I

 

 

(G) 8I

 

(G))1/2

,

(1.60)

 

 

 

2

 

 

dV

 

 

1

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ˆ

 

 

 

 

 

 

 

 

 

 

 

 

dV

(1 2I ( A)

4I

 

 

( A) 8I

( A)) 1/2

,

(1.61)

 

2

dV

 

 

1

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ˆ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dV

 

 

 

 

 

ˆ

 

 

 

 

 

 

dV

exp I1

(H )

exp I1(H ) .

 

(1.62)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1.8. Перемещения и дисторсии. Малые деформации

Вектор перемещения точки тела из отсчетной конфигурации в

текущую определяется следующим образом:

 

u x x .

(1.63)

С помощью этого векторного поля на теле аффиноры могут быть записаны как

F x x u I u ,

(1.64)

F x x u I u ,

где диады u и u называют тензорами дисторсий. Эти тензоры несут в себе столько же информации, сколько и аффиноры, поэтому с их помощью можно записать все меры деформаций. Например, меры Грина и Альманси записываются как

23

C

1

( u u u u ),

 

2

(1.65)

 

 

 

A

 

ˆ

ˆ ˆ ˆ

 

 

2

( u u u u ).

 

 

1

 

 

 

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

 

y(

 

 

i )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

базисом e

( y( i ) x( i ) — радиус-вектор точки пространства,

i

 

 

 

 

 

i

i

 

 

 

 

 

 

координаты точки пространства, в которой первоначально находилась материальная точка с координатами i ) и набла-вектором e i i ei i ,

без привлечения текущей системы материальных координат. Аргументами всех определенных на континууме полей берут пространственные координаты i . Близость конфигураций подразумевает |u| 1 и, что более важно, |u | 1. Все меры деформации в этом случае до малых более высокого порядка по u сведутся к тензору

 

1

( u u ) ,

(1.66)

 

2

 

 

называемому тензором малых деформаций. Относительно тензоров Грина и Альманси в этом можно убедиться с помощью (1.65). Оператор

def ( )

1

( ( ) ( ) )

(1.67)

 

2

 

 

называют деформатором. С его помощью можно записать

 

def u.

(1.68)

Антисимметричная часть тензора малой дисторсии

 

 

1

(u u)

(1.69)

 

2

 

 

 

 

называется тензором малых поворотов. Из (1.66) и (1.69)

 

u .

(1.70)

Вспоминая (1.7) и (1.63), нетрудно увидеть

 

du ( ) dx ,

(1.71)

то есть смысл тензоров малых деформаций и поворотов как линейных операторов. Относительное перемещение du концов малого материального

отрезка dx вызывается его деформацией dx

и поворотом dx .

Антисимметричному тензору малых поворотов

ставят в соответствие

аксиальный вектор

 

 

 

 

 

 

 

 

1

Є: ,

(1.72)

 

 

2

 

 

 

 

24

 

обратно
Є .
С использованием (1.69) этот вектор малых поворотов
поле перемещений
21 u.

(1.73)

выражается через

(1.74)

Очевидно, что главные оси тензора малых деформаций совпадают с предельными главными осями pi эллипсоидов деформации.

Малые упругие деформации без поворотов и перемещений испытывают, например, детали строительных конструкций. Главные относительные удлинения сводятся к i 1 i до малых второго порядка по

i . При этом все шкалы деформации f ( i ) с той же точностью стремятся к

i . Поскольку локальными поворотами можно пренебречь, меры

деформации сводятся к тензору малых деформаций.

Существуют среды, физически неспособные испытывать деформации ( i 0) , но при этом имеют конечные локальные перемещения и повороты. Таковы, например, нерастяжимые гибкие нити и ткани. Надуваемые резиновые сферические оболочки, напротив, обнаруживают большие деформации без поворотов. Другие же среды имеют сравнимые (малые или немалые) деформации и повороты. При необходимости можно построить локальные меры деформаций и ротаций, наиболее подходящие рассматриваемому процессу.

1.9. Условия совместности деформаций

Итак, все меры деформаций определяются полем перемещений точек материального континуума. Попробуем сейчас решить обратную задачу — по заданному полю меры деформации определить поле перемещений. Этот вопрос потребуется далее для постановки краевой задачи в напряжениях. Потребуем только непрерывности этого поля меры деформации вплоть до производных второго порядка, а в остальном будем считать его вполне произвольным. Сначала решим эту проблему для малых деформаций. Математически она эквивалентна интегрированию системы шести дифференциальных уравнений (1.66)

u u 2

с заданными правыми частями требуемой непрерывности. В ортобазисе ki декартовой пространственной системы координат эти уравнения запишутся в виде

25

u1x1

u2x1

 

11

,

u2

22

,

u3

33

,

 

 

 

 

 

 

 

 

 

 

x2

 

 

x3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u1

2

 

,

u3 u1

2

 

,

u3

u2

2

 

.

 

12

13

 

23

 

x2

 

 

x1

x3

 

 

 

 

x2

x3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Интегрированием первых трех уравнений можно получить компоненты u1 ,u2 ,u3 с точностью до трех произвольных функций, которые следует найти из последних трех уравнений. Выясним, всегда ли можно найти по заданному

полю поле

u. Иначе говоря, мы должны найти условия существования

непрерывного

поля u, удовлетворяющего (1.66). Очевидно, что

(1.66)

эквивалентно (1.71) при условии (1.69). Из (1.71)

 

 

du ( ) dx ,

 

следует

 

 

 

M

 

 

u u0 ( ) dx ,

(1.75)

 

M 0

 

где u0 — перемещение точки M0 . Если материальное тело занимает односвязную область пространства, то необходимое и достаточное условие однозначности u, т.е. его независимости от пути интегрирования (напомним, что для Q dx это условие имеет вид QT 0 [1]) в (1.75), записывается в виде

( ) 0 .

(1.76)

Однако нам ничего не известно о существовании поля , кроме связи (1.69) его с искомым полем u. В последнем равенстве

Є Є Є (CII CIII )

(kik j kik j ki Iki ) i j k j ki I ,

поскольку

 

0

в силу (1.69). Равенство (1.76) тогда принимает вид

 

 

 

 

 

 

 

 

d ( ) dx ,

 

(1.77)

 

 

 

 

 

 

откуда

 

 

 

 

 

 

 

 

M

 

 

 

 

 

0

 

(1.78)

 

 

 

dx .

 

 

 

 

M0

 

 

Условие существования однозначного поля

 

в односвязной области

 

 

 

 

 

пространства эквивалентно условию независимости его от пути интегрирования:

0 .

(1.79)

Это и есть искомые нами условия совместности Сен-Венана. Эти ограничения задают в функциональном пространстве произвольных

26

достаточно гладких полей симметричного тензора на теле подпространство совместных деформаций. Уравнения (1.79) иногда записывают с использованием символического оператора несовместности

Ink :

Ink 0 .

(1.80)

Читателю следует убедиться, что

 

Ink def 0.

 

Дадим геометрическую интерпретацию полученного условия. Сначала разделим еще недеформированную сплошную среду на бесконечно малые материальные кубики и разнесем их в пространстве, чтобы они не мешали друг другу. А затем сообщим каждому из них деформацию, соответствующую достаточно непрерывному, а в остальном произвольному полю деформаций, заданному на исходной конфигурации K0 материальной среды. Можно ли обратно склеить из деформированных параллелепипедов какое-либо непрерывное тело?

Чтобы ответить на этот вопрос, попробуем описать то, что мы сделали, в терминах дифгеометрии, при этом не ограничиваясь малыми деформациями. Пусть в K0 задана декартова система лагранжевых координат

i с однородной по этой

конфигурации

метрикой g .

На K

задается

 

 

 

 

 

 

 

 

 

 

 

 

ij

0

 

произвольное поле деформаций aij , в силу (1.23) изменяющее метрику:

 

gij

gij 2aij .

 

 

 

 

 

(1.81)

Эта метрика индуцирует ненулевые объект связности

 

 

 

 

 

 

ˆ

 

 

ˆ

 

ˆ

 

 

 

 

ˆ

1

 

g

ij

 

 

g

jk

 

 

 

 

 

gik

 

 

 

 

i jk

 

 

(

 

 

 

 

 

)

 

(1.82)

2

k

j

i

 

и тензор Римана ‒ Кристоффеля (тензор кривизны пространства)

 

 

 

 

 

2

ˆ

 

2 ˆ

 

 

 

2

ˆ

 

 

2 ˆ

 

 

 

 

 

ˆ

 

1

 

 

 

gil

 

 

jl

 

 

 

gik

 

 

 

jk

 

)

 

 

Rijkl

 

(

 

 

 

 

 

 

 

2

j k

i k

j l

i l

 

 

(1.83)

 

 

 

 

 

 

 

 

 

 

 

 

 

gˆ

ps

ˆ

 

ˆ

ˆ

ˆ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( p jk

s il

p jl

s ik ),

а это означает, что деформированный материальный континуум образует риманово пространство. Последнее сводится к аффинному пространству при условии

 

(1.84)

Rijkl 0.

Записывая это условие через aij с использованием (1.81) ‒

(1.83), мы и

получим условие совместности деформаций

 

27

 

 

ˆ

 

 

2a

 

2a jl

 

2a

 

 

2a jk

 

 

 

 

 

 

 

 

 

il

 

 

 

 

 

 

ik

 

 

 

 

 

 

 

 

 

 

 

 

Rijkl

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j k

i k

j l

i l

 

(1.85)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

gˆ ps ( A

 

A

A A

) 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p jk

s il

 

p jl s ik

 

 

где A

 

aij

aik

 

a jk

 

и [gˆ ps ] [

 

2a

 

] 1 , линейная часть которого (

 

 

 

 

ps

ps

i jk

 

k

j

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

aij ij

1) имеет вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 il

 

 

2 jl

 

 

2 ik

 

2 jk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 ,

 

(1.86)

 

 

 

 

 

j k

i k

 

j l

i l

 

полностью совпадающий

 

с

уравнением (1.79),

записанным в i

как в

декартовых пространственных координатах. Условие (1.84) требует, чтобы конфигурация Kt принадлежала пространству без кривизны, то есть аффинному пространству. В общем же случае задаваемая произвольным полем деформации конфигурация Kt оказывается принадлежащей риманову пространству или, иначе говоря, не существует поля перемещений, переводящих K0 в Kt в аффинном пространстве.

Меру деформации aij в (1.85), конечно, можно выразить через любой

из тензоров конечных деформаций. Условие (1.85) или (1.86) содержит 6 независимых уравнений в трехмерном случае, одно в двумерном и ни одного в одномерном (столько независимых компонент у тензора Римана ‒ Кристоффеля).

1.10. Определение перемещений по совместным деформациям

Ограничимся малыми деформациями. Используя (1.69), (1.75), (1.78), получим

 

 

 

M

 

M

 

M

 

 

 

 

 

 

 

 

( ) dx .

u u0 0

( x x0 )

 

dx

 

dx

 

 

 

 

M 0

 

M 0

 

M 0

 

Далее, сводя двойной

интеграл

к однократному: b ds s f (u)du

 

 

 

 

 

 

 

 

a a

ab f (v)(b v)dv , получим

M

u u0 0 ( x x0 ) ( x x) ( ) dx (1.87)

M0

формулу Чезаро. Очевидно, интеграл в (1.87) не зависит от пути. Перемещение, естественно, оказалось определенным с точностью до малого

28

перемещения среды как жесткого целого u0 0 (x x0 ), складывающегося

из трансляции точки M0 и поворота вокруг нее.

Формула Чезаро ввиду громоздкости подынтегральных выражений для вычисления перемещений не применяется, для чего имеется другой путь. Сначала находятся компоненты тензора дисторсии. Диагональные компоненты тензоров дисторсии и малых деформаций совпадают, а недиагональные находятся так. Из (1.66), (1.69), (1.70) в ортобазисе декартовой пространственной системы координат получаем

ij ,k

ik , j jk ,i .

Затем из (1.69)

 

 

ui , jk ij ,k ik , j jk ,i ,

откуда

 

 

M

 

 

ui, j ui0, j ( ij ,k

ik , j jk ,i )dxk .

M 0

 

 

Интеграл в последнем выражении берется по ломаной со звеньями, параллельными осям координат (примем M0 (0,0,0) , M (x1, x2 , x3 ) ):

 

 

 

x1

 

 

 

 

 

 

 

 

ui, j ui0, j

( ij,1 i1, j

j1,i )

dx1

 

 

 

 

 

0

 

 

 

 

 

x2 x3 0

 

 

 

 

 

 

 

 

 

 

 

 

x2

 

 

 

 

 

 

x3

 

 

 

 

( ij,2

i 2, j

j 2,i )

 

 

 

dx2 ( ij,3 i3, j j3,i )dx3 ,

i j.

 

 

 

 

x3

0

 

 

 

 

0

 

 

 

 

 

 

0

 

 

 

 

Далее

 

 

 

 

 

 

 

 

 

 

 

 

 

x1

 

 

 

 

x2

 

 

x3

 

ui

ui0 ui,1

 

 

 

 

dx1 ui,2

 

dx2 ui,3dx3 .

 

 

 

 

x2

x3

0

 

x3 0

 

 

 

 

0

 

 

 

 

0

 

 

0

 

29

2.Кинематика сплошной среды

2.1.Материальная производная

На деформируемом материальном континууме бывают определены некоторые физические поля (например, перемещений или сил), изменяющиеся со временем. Для определения скорости изменения значений этих полей в фиксированных материальных точках используется

материальная производная:

 

( )

 

 

 

 

 

 

( )

 

.

(2.1)

 

 

t

 

 

x

 

 

 

 

При лагранжевом способе описания движения все поля рассматриваются как функции лагранжевых координат (не зависящих от времени) и времени, и потому для любого такого поля c( x,t ) правомерны записи (2.1)

в виде

c( x,t) c( x,t)

 

dc( x,t)

.

 

 

 

t

 

dt

 

Скорость и ускорение движения точки континуума определяются как

 

v(x,t) xˆ ( x,t)

(2.2)

a( x,t) xˆ ( x,t) .

(2.3)

2.2. Эйлеров способ описания движения

Эйлеров способ описания движения континуума заключается в рассмотрении всех заданных на нем полей как функций пространственных координат и времени. Чтобы перейти от лагранжева описания к эйлерову, нужно знать, какая именно материальная точка находится сейчас в данной точке пространства. Связь (1.2)

x x(xˆ,t)

как раз и содержит в себе такую информацию. Подставляя в c( x,t) эту связь, получим c c(x(xˆ,t),t) cˆ( xˆ,t) — эйлерово описание поля c . Ранее мы уже встречались с разновидностью эйлерова описания, когда рассматривали приближение малых деформаций.

Пусть требуется записать материальную производную поля c cˆ( xˆ,t) , т.е. заданного эйлеровым способом как функции пространственных координат и времени. В соответствии с (2.1) для этого необходимо взять производную c cˆ( xˆ,t) по t при фиксированном x . Очевидно фиксировать

30