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

книги / Механика сплошных сред (теоретические основы обработки давлением композитных материалов с задачами и решениями, примерами и упражнениями)

..pdf
Скачиваний:
5
Добавлен:
12.11.2023
Размер:
9.02 Mб
Скачать
1.2. КИНЕМАТИКА
Рис. 7. Схема перемещения материаль9 ной частицы m

Вектор перемещения этой частицы имеет вид

U 4e1 18e2 3e3 .

Из (1.2.4) имеем E U L . Отсюда определяем радиус-вектор пространственной точки, в которую переместилась заданная свои-

ми исходными координатами материальная частица

E5e1 20e2 6e3 .

Задача 1.2.1.3. В координатах Л. Эйлера задан вектор перемещения

UE3E12e1 E2E32e2 E1E22e3 .

Определить исходные координаты материальной частицы m, которая в результате движения попадает в пространственную точку nχ (14, 17, 1).

Решение. Радиус-вектор заданной пространственной точки

E 14e1 17e2 e3 .

Вектор перемещения материальной частицы, которая попадает в эту точку (рис. 7), определяется формулой (1.2.4), из которой имеем

LE U .

Отсюда определяем радиус-вектор, характеризующий положение материальной частицы до движения (рис. 7):

L2e1 8e2 e3 .

1.2.2.Описание движения в лагранжевых координатах

Назначим базис ei , связанный с материальным множеством координат Li. Наглядным примером координирования материальных частиц m в таком базисе является координатная сетка, наносимая на образец до исследования его движения в начальный момент времени t = t0 (рис. 8, а). Обычно исследователь осуществляет наблюдение за движением изучаемого тела, фиксируя положения его отдельных материальных частиц m в определенные моменты времени t в

31

1. МЕХАНИКА СПЛОШНЫХ КОМПОЗИТНЫХ СРЕД
отдельных пространственных точках n. В этом случае, находясь в пространственных координатах Ek, он видит искажение в процессе движения лагранжевых изоповерхностей или изо-
линий Li = const (рис. 8, б). Поэтому
Рис. 8. Координатная сетка до (а) и после (б) де9 эйлеровы координаты иногда называ-
формации
ют отсчетом наблюдателя. В координатах же Li конфигурации объектов Li = const, в том числе и форма образца, остаются неизменными и происходит
изменение конфигурации изоповерхностей или изолиний Ek = const. Рассмотрим движение двух материальных частиц m1 и m2 (рис. 9, а) в лагран-
жевом базисе ei . В качестве меры изменения длины элементарного волокна dL (расстояния между m1 и m2) будем использовать коэффициент изменения длины
волокна

Φ || dE ||

 

|| dL || .

(1.2.25)

Попадая из пространственных положений n1 и n2 в пространственные положения n1χ и n2χ соответственно, материальные частицы m1 и m2 могут сближаться (Φ < 1) или расходиться (Φ > 1). При неизменном расстоянии между этими частицами Φ= 1. С помощью (1.2.12) вычислим квадрат длины волокна dE в момент времени t:

Рис. 9. Движение в лагранжевых координатах:

а – произвольная ориентация вектора dL ; б – вектор dL параллелен оси L1

32

 

 

 

1.2. КИНЕМАТИКА

|| dE ||2 = dE dE = [(E ) dL] [(E ) dL] .

(1.2.26)

В этом пункте оператор У.Р. Гамильтона (П1.74) =

 

 

.

 

Lj

 

Упражнение 1.2.1. Доказать, что, жонглируя индексами*, соотношение (1.2.26)

можно привести к виду

 

 

 

|| dE ||2 = T (dL dL) ,

(1.2.27)

 

 

G

 

где тензор

 

 

 

 

TG = [[Gik ]] = ( E) c(E )

(1.2.28)

называется тензором конечных деформаций А. Е. Грина

 

Учитывая, что единичный вектор, совпадающий с dL

по направле-

нию, =

dL

, запишем значение квадрата коэффициента (1.2.25):

|| dL ||

 

 

 

 

 

χ = TG ( ).

(1.2.29)

Формула (1.2.29) позволяет дать физическую интерпретацию диагональным компонентам тензора А.Е. Грина (1.2.28). Каждая такая компонента

Gii =

Em Em

(1.2.30)

Li

Li

 

 

равна квадрату коэффициента χi изменения длины элементарного волокна до начала движения, параллельного оси Li. Действительно, положив, напри-

мер, dL = dL1e1 (рис. 9, б), получим 1 =1 и

2 = 3 = 0 . Это приведет к тому, что

χ2 = χ12 в (1.2.29) будет равно G11. Аналогичным образом это можно показать

для волокон dL , параллельных осям

L2 ( 1 = 0; 2 =1; 3 = 0) или L3

( 1 = 0; 2 = 0; 3 =1) . В общем случае

 

*В одночленах с повторяющимися индексами вид буквы, используемой для обозначения одинаковых индексов, не имеет значения. Поэтому такие индексы называют немыми. Замена одной буквы немых индексов на другую называется жонглированием индексами.

33

1. МЕХАНИКА СПЛОШНЫХ КОМПОЗИТНЫХ СРЕД

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Φi

 

 

Gii .

 

 

 

 

 

 

 

(1.2.31)

Таким образом, диагональные (m = k = i) компоненты Gmk тензора (1.2.28)

характеризуют изменение линейных размеров волокна окрестности материаль-

ной частицы, параллельного при t = t0 координатной оси Li, и называются ли

нейными конечными деформациями тензора А. Е. Грина.

 

 

 

 

 

Наиболее часто встречающуюся в практике обработки металлов давлением

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

G

dE d L

F 1

 

 

 

 

 

 

 

 

 

d L

 

 

 

 

 

 

 

(1.2.32)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

также можно записать с помощью компонент тензора (1.2.28).

 

 

 

Для волокна, параллельного при t = t0 оси Li,

 

 

 

 

 

 

 

 

 

 

 

 

Γi

 

Gii 1.

 

 

 

 

 

 

 

(1.2.33)

Теперь рассмотрим движение двух элементарных волокон, например, dLχ и

dLχχ , параллельных при t = t

0

осям L

1

и L

2

(рис. 10). В результате движения к

моменту времени t элемент dLc

dL e

займет положение dEχ , а второй элемент

 

 

 

 

 

1 1

 

 

 

 

 

 

 

 

 

 

 

 

 

dLcc

dL e

– положение dEχχ

. Определим угол

 

между векторами dEχ

и dEχχ .

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Косинус искомого угла найдем как скалярное произведение единичных векто-

ров, совпадающих по направлению с векторами dEχ

и dEχχ :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dEc

dEcc

 

 

 

 

 

 

 

 

 

 

 

 

 

cos

= || dEc||

|| dEcc|| .

(1.2.34)

 

 

 

 

 

 

 

 

 

 

В

 

соответствии

с

(П1.83)

 

 

 

 

 

 

 

 

dEχ

ωEi dL1ei

и dEχχ

 

ωEk dL2ek .

 

 

 

 

 

 

 

 

 

 

 

ωL1

 

 

 

 

ωL2

 

 

 

 

 

 

 

 

 

После подстановки этих результатов

 

 

 

 

 

 

 

 

в (1.2.34) и несложных преобразова-

 

 

 

 

 

 

 

 

ний, учитывая (1.2.28) и (П 1.12), по-

 

 

 

 

 

 

 

 

лучим cos

=

G12

.

 

 

Рис. 10. Схема к вычислению сдвиговой деформа9

G11

G22

 

 

 

 

 

 

 

 

 

 

 

ции в лагранжевых координатах

 

 

 

 

 

 

 

 

 

 

 

 

 

34

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1.2. КИНЕМАТИКА

Введем угол

ϑ

 

Σ

 

, дополняющий угол до прямого угла. Тогда

 

 

 

 

 

12

2

 

 

 

 

G12

 

 

 

 

sinϑ12

=

 

. Обобщая полученный результат, синус угла ϑik, дополняю-

G11 G22

 

 

 

 

 

 

щего до прямого угол между волокнами, которые при t = t0 были параллельны

осям Li и Lk при i ζ k, определяем по формуле

sin Jik

Gik

 

.

(1.2.35)

Gii

 

 

Gkk

 

Ранее (1.2.31) было установлено, что диагональные компоненты тензора (1.2.28) ответственны за изменение линейных размеров окрестности материальной частицы. Из (1.2.35) следует, что при Gik = 0 (i ζ k) не происходит изменения исходного прямого угла. Следовательно, изменение угловых размеров между волокнами этой окрестности, параллельными при t = t0 координатным осям Li и Lk, в основном характеризуется боковыми компонентами тензора (1.2.28), которые называются угловыми (сдвиговыми) конечными деформациями тензора А. Е. Грина.

Рассмотрим объем элементарного прямоугольного параллелепипеда с реб-

χ

dL

χχ

и dL

χχχ

:

 

 

 

 

 

 

рами dL ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d:

c

cc

 

ccc

(1.2.36)

 

 

 

 

 

 

(dL

udL )

 

dL .

 

 

 

 

 

 

 

L

 

 

 

 

 

Пусть в момент времени t = t0 его ребра будут параллельны осям L1, L2 и L3

соответственно (рис. 11). Тогда

 

 

 

 

d:L = dL1dL2dL3. К моменту време-

 

 

 

 

ни t в общем случае это тело станет

 

 

 

 

непрямоугольным параллелепипе-

 

 

 

 

дом с ребрами dEχ,

dEχχ и dEχχχ соот-

 

 

 

 

ветственно, а его объем будет

 

 

 

 

 

d:

(dEcudEcc) dEccc .

(1.2.37)

 

 

 

 

E

 

 

 

 

 

 

 

 

 

 

 

Используя (П1.14), (П1.83), (1.2.13)

 

 

 

 

и (1.2.14), из (1.2.36) и (1.2.37) находим

 

 

 

 

соотношение между объемами

 

 

 

 

 

Рис. 11. Схема к вычислению изменения объема в лагранжевых координатах

35

1. МЕХАНИКА СПЛОШНЫХ КОМПОЗИТНЫХ СРЕД

dΩE = JLdΩL.

(1.2.38)

Отсюда при dΩE = dΩL получаем условие постоянства объема в материальных координатах:

JL = 1.

(1.2.39)

Таким образом, выполненным анализом установлен физический смысл компонент тензора (1.2.28) и якобиана (1.2.13): диагональные компоненты этого тензора характеризуют изменение линейных размеров окрестности движущейся материальной частицы, боковые компоненты определяют изменение угловых размеров этой окрестности, а якобиан преобразования координат равен отношению текущего значения ее объема к его начальному значению.

При анализе движения окрестности частицы m удобно вместо тензора (1.2.28) использовать лагранжев тензор конечных деформаций

TL =[[Lij ]] =

1

(TG Tδ).

(1.2.40)

2

 

 

 

В связи с тем, что в (1.2.40) постоянный сомножитель перед скобками и единичный тензор не вносят дополнительной физической информации, физический смысл компонент Lij и Gij полностью совпадает. Диагональные компоненты тензора (1.2.40) называются линейными конечными деформациями тензора Ж. Лаг ранжа, а боковые – угловыми (сдвиговыми) конечными деформациями тензора

Ж.Лагранжа.

Тензор (1.2.40) можно представить в виде функции вектора перемещения.

Для этого из (1.2.11) находим

E(Li , t) = U (Li , t) + L .

(1.2.41)

Подстановкой (1.2.41) в (1.2.28), учитывая, что в лагранжевых координатах

L = L = Tδ , из (1.2.40) получим

TL =

1

( U +U + U cU ).

(1.2.42)

2

 

 

 

Задачи к пп. 1.2.2

Задача 1.2.2.1. Положение материальной частицы в некотором объеме сплошной среды к определенному моменту времени характеризуется радиус-вектором

E = L1e1 + (L2 + 3L3 )e2 + (L3 + 3L2 )e3 .

36

1.2. КИНЕМАТИКА

Определить компоненты тензоров конечной деформации А. Е. Грина и Ж. Лагранжа и проанализировать их физический смысл.

Решение. Компоненты Gik тензора конечных деформаций А. Е. Грина TG найдем подстановкой в (1.2.29) заданного закона движения:

 

 

E

2

 

 

E

2

 

 

G

 

= G

 

=

E E

+

E E

2 +

E

E

G

=

 

1

+

 

 

3

 

=1;

 

 

 

1

1

 

2

 

3

3 = 0;

 

 

 

 

 

12

 

21

 

 

L L

 

 

L L

 

L

L

11

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L1

 

 

 

 

L1

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

 

 

1

2

1

2

 

 

 

 

 

G

 

= G

= E1 E1

+

E2 E2

+ E3 E3

=

0;

 

 

 

 

 

 

 

 

13

 

31

 

L

 

L

 

 

L

 

L

L

 

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

3

 

 

 

1

 

3

 

1

3

 

 

 

 

 

 

 

 

 

 

 

 

G

=

 

E1

2

+

 

E2

2

+

E3

2

=10;

 

 

 

 

 

 

 

 

 

 

 

22

 

 

L2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L2

 

 

L2

 

 

 

 

 

 

 

 

 

 

 

G

 

= G

= E1 E1

+

E2 E2

+ E3 E3

=

6;

 

 

 

 

 

 

 

 

23

 

32

 

L

 

L

 

 

L

 

L

L

 

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

3

 

 

 

2

 

3

 

2

3

 

 

 

 

 

 

 

 

 

 

 

 

G

=

 

E1

2

+

 

E2

2

+

E3

2

=10.

 

 

 

 

 

 

 

 

 

 

 

33

 

 

L3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L3

 

 

L2

 

 

 

 

 

 

Таким образом,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

0

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

TG =

0

10

 

6

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

6

 

10

 

 

 

 

 

 

 

 

Компоненты тензора конечных деформаций А. Е. Грина имеют конкретный физический смысл. Диагональные i-е компоненты характеризуют изменение длины волокна, параллельного в начальный момент времени t0 лагранжевой оси Li, и с их помощью можно вычислить коэффициент изменения длины такого волокна (1.2.31). Подставляя в эту формулу компоненты полученного тензора

А. Е. Грина, устанавливаем, что волокно dL′, параллельное до движения оси L1

(рис. 12), в процессе движения не изменяет своей длины (dL′ = dE; χ1 = 1), а

37

1. МЕХАНИКА СПЛОШНЫХ КОМПОЗИТНЫХ СРЕД

 

волокна dLχχ и dLχχχ , параллельные до движе-

 

ния осям L2 и L3 соответственно, в процессе

 

движения

удлиняются

в

10

раз

 

( Ȥ2 Ȥ3

10 ; dEχχ ! dLχχ и dEχχχ ! dEχχχ ). Боко-

 

вые ik-е (i ζ k) компоненты тензора конечных

 

деформаций А. Е. Грина характеризуют изме-

 

нение угловых размеров между волокнами, па-

 

раллельными в начальный момент времени t0

 

лагранжевым осям Li и Lk. С помощью компо-

Рис. 12. Изменение волокон в лагран9

нент Gik вычисляется косинус угла

ik

между

жевых координатах

этими волокнами (1.2.34). В связи с тем, что для

 

полученного тензора А. Е. Грина G12 = G21 = 0 и

G13 = G31 = 0, отмечаем низменность прямого угла между волокнами, параллель-

ными в начальный момент времени t0 осям L1, L2 ( 12 = 0) и осям L1, L3 (

13 = 0)

соответственно. Прямой до движения угол между волокнами, параллельными осям L2 и L3, изменился (рис. 12):

cos 23

6

 

 

=

 

 

0,6 .

 

10

 

10

 

Отсюда 23 | 53о8χ.

Компоненты Lik тензора конечных деформаций Ж. Лагранжа TL связаны с компонентами тензора А. Е. Грина формулой (1.2.40), которая в скалярной форме имеет вид:

1

Lik 2 (Gik Γik ) ,

где Γik – компоненты единичного тензора (символ Л. Кронекера). Тогда тензор конечных деформаций Ж. Лагранжа

 

1

1 1

0

0

0

0

0

 

T

0

10 1

6

0

4,5

3

 

 

.

L

 

 

 

 

 

 

 

2

0

6

10 1

0

3

4,5

 

 

 

 

С помощью компонент тензора конечных деформаций Ж. Лагранжа, так же как и по тензору А. Е. Грина, можно установить физическую картину деформа-

38

1.2. КИНЕМАТИКА

ции окрестности материальной частицы. Так, из полученного тензора TL видно, что деформация волокна, параллельного до движения оси L1, отсутствует, так как L11 = 0; волокна, параллельные осям L2 и L3, удлиняются (L22 > 0 и L33 > 0). Сдвиговая деформация отлична от нуля только между волокнами, параллельными в начальный момент времени t0 осям L2 и L3, так как L12 = L21 = 0,

L13 = L31 = 0 и L23 = L32 0.

Задача 1.2.2.2. По условиям задачи 1.2.2.1 определить изменение квадрата длины направленного волокна dL .

Решение. В результате движения направленного волокна dL от момента времени t0 к произвольному моменту времени t оно займет положение dE . Для

вычисления разности между квадратами модулей || dE ||2 || dL ||2 этих векторов воспользуемся формулой (1.2.29) для вычисления квадрата коэффициента изменения длины волокна, которая в скалярной форме имеет вид

 

 

χ2 = Gik

i

 

k ,

 

 

 

 

где направляющие косинусы вектора dL

 

 

 

 

 

 

 

 

 

 

j

=

 

dLj

 

.

 

 

 

 

 

 

 

 

|| dL

||

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда с учетом (1.2.29) имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

|| dE ||2 || dL ||2 = 2T

dL dL .

 

 

 

 

 

 

 

 

 

 

 

L

 

 

 

 

 

Теперь, используя результаты предыдущей задачи, получим

 

 

 

 

0

 

0

0

 

 

dL1

 

 

 

 

|| dE ||2 || dL ||2 =[[dL dL dL ]] 0

 

9

6

 

 

dL

= 9dL2

+12dL dL

+ 9dL2

1

2

3

 

 

 

 

 

 

2

2

2

3

3 .

0 6 9 dL3

Отсюда видно, что общее изменение длины волокна, так же, как это было установлено в предыдущей задаче, связано лишь с изменениями проекций вектора dL на вторую и третью лагранжевы оси координат (рис. 12).

Задача 1.2.2.3. Используя тензоры TG и TL задачи 1.2.2.1, определить главные компоненты тензора конечных деформаций А. Е. Грина и тензора конечных деформаций Ж. Лагранжа.

39

1. МЕХАНИКА СПЛОШНЫХ КОМПОЗИТНЫХ СРЕД

Решение. Составим характеристическое уравнение матрицы тензора конечных деформаций А. Е. Грина, полученного в задаче 1.2.2.1:

|Gij – λδij| = 0

или

(1 – λ)(10 – λ)(10 – λ) – 36(1 – λ) = 0.

Отсюда λ1 = 1 и

λ2 – 20λ + 64 = 0.

Решая это уравнение, получим λ2 = 16; λ3 = 4. Соблюдая условие

G1 G2 G3,

составим диагонализированный вид тензора конечных деформаций А. Е. Грина:

16 0 0

TG = 0 4 0 .

0 0 1

Тогда в соответствии с (1.2.40) тензор конечных деформаций Ж. Лагранжа в главных координатах

 

 

1

16 1

0

0

7,5

0

0

TL

=

 

0

4 1

0

= 0

1,5

0 .

2

 

 

0

0

11

0

0

0

 

 

 

Эти же результаты можно получить, решая полные кубические характеристические уравнения –λ3 + GIλ2 GIIλ + GIII = 0 и –λ3 + LIλ2 LIIλ + LIII = 0, где инварианты тензора конечных деформаций А. Е. Грина

GI = Gii = 21; GII = 12 (GI2 GijGji ) = 84; GIII = |Gij | = 64,

а инварианты тензора конечных деформаций Ж. Лагранжа

LI = Lii = 9; LII = 12 (LI2 Lij Lji ) = 11,25; LIII = |Lij | = 0.

Самостоятельно проанализировать физический смысл компонент тензоров конечных деформаций А. Е. Грина и Ж. Лагранжа в главных координатах.

40