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

Tensor-Gotman

.pdf
Скачиваний:
55
Добавлен:
09.02.2015
Размер:
1.52 Mб
Скачать

Произведение тензоров и векторов

Тензор символически записывается в виде

 

 

 

 

 

 

t11

t12

t13

 

 

t11

t12

 

 

 

 

 

 

 

 

Tik

 

 

 

=

t21

t22

t23

 

или

T = t21

t22

 

 

 

 

 

 

 

 

 

 

t31

t32

t33

 

 

 

t32

 

 

 

 

 

 

 

 

t31

Вектор записывается в виде

a = e1a1 +e2a2 + e3a3.

t13

t23 (11.11)

t33

(12.11)

Каждый тензор второго ранга можно представить в подобном виде, то есть, в виде трёх векторов или записать его как сумму трёх диад (диадика)

T = e1t1 + e2t2 +e3t3

(13.11)

Это получается следующим образом: пусть матрица – столбец ортов и матрица – строка имеют вид

e

 

 

 

 

 

 

1

 

 

E = (e1

, e2

, e3 )

 

E = e2

 

или

(14.11)

 

 

 

 

 

 

 

e3

 

 

 

 

 

 

а векторы тензора записаны в форме

 

 

 

 

 

t1 = e1t11 + e2t12 +e3t13

 

 

 

t2 = e1t21 +e2t22 +e3t23

 

 

(15.11)

t3 = e1t31 + e2t32 +e3t33

 

 

 

потому что произведение тензора на столбец ортов справа получается в виде

t

 

t

 

t

 

 

e

 

 

e t

+ e

t

 

+e

 

t

 

 

t

 

 

 

11

12

13

 

1

 

 

1 11

 

2 12

 

3 13

 

1

 

 

T E = t21

t22

t23

 

e2

 

= e1t21 +e2t22 + e3t23

 

= t2

 

(15.11)

t

31

t

32

t

33

 

e

3

 

e t

+ e

2

t

32

+e

3

t

33

 

t

3

 

 

 

 

 

 

 

 

 

1 31

 

 

 

 

 

 

 

 

 

Если умножить матрицу-столбец (14.11) слева на E , получим исходный тензор

e1 T = E T E = e2e3

 

 

e1t11 +e2t12 +e3t13

 

t11

 

 

e t

+e

2

t

22

+ e

3

t

23

=

t

21

 

 

 

1 21

 

 

 

 

 

 

 

 

e t

+e

2

t

32

+e

3

t

33

 

t

31

 

 

 

1 31

 

 

 

 

 

 

t12 t22 t32

t13

t23 (16.11)

t33

Умножение тензора Т скалярно на вектор а даёт выражение нового вектора a' a' = T a = e1 (t1 a) + e2 (t2 a) + e3 (t3 a) =

= e1 (t11a1 + t12a2 + t13a3 ) +e2 (t21a1 + t22a2 + t23a3 ) +

(17.11)

+ e3 (t31a1 + t32a2 +t33a3 )

Произведение тензора Т на вектор а с последующим свёртыванием

Определение 4.11 Скалярным произведением тензора Т на вектор а называется вектор a' , составляющие которого линейным образом выражаются через составляющие вектора а, причём коэффициентами являются компоненты тензора Т.

Замечание 6.11 Вектор a' = T a называется линейной векторной функцией вектора а.

Рассмотрим произведение, в котором производится свёртывание по первому индексу тензора Τik Ai . Это получается при умножении сопряженного (транспонированного)

тензора

51

t11 t21 t31 Tc = Tik = t12 t22 t32 t13 t23 t33

на вектор справа

t

a +t

21

a

2

+ t

31

a

3

 

t

t

21

11 1

 

 

 

 

 

 

11

 

Tc a = t12a1 +t22a2 + t32a3

 

= t12

t22

t

a +t

23

a

2

+t

33

a

3

 

t

t

23

13 1

 

 

 

 

 

 

13

 

или исходного тензора на вектор слева

t

31

 

a

 

 

 

 

1

 

= Tik a i

t32

 

a2

 

t33

 

 

 

 

 

a3

 

 

 

t

t

 

11

12

a T = (a1 a2

a3 ) t21

t22

 

 

t32

 

t31

t

 

 

a t

+ a

2

t

21

+ a

3

t

31

 

13

 

 

1 11

 

 

 

 

 

t23

 

= a1t12 +a2t22 + a3 t32

 

t

33

 

a t

+ a

2

t

23

+ a

3

t

33

 

 

 

 

1 13

 

 

 

 

 

Тогда произведение вектора а на тензора Т получается в виде

 

 

 

 

 

 

 

 

 

a' = a T =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

) (18.11)

= e

(t

a +t

21

a

2

+t

31

a

3

) +e

2

(t

a +t

22

a

2

+t

32

a

3

) +e

3

(t

a + t

23

a

2

+t

33

a

3

1

11

1

 

 

 

 

12

1

 

 

 

 

13

1

 

 

 

 

Замечание 7.11 Из получения произведений тензора Т на вектор а видно, что можно использовать исходный и сопряжённый тензор для получения одного и того же выражения, поэтому справедливо равенство

a T = Tc a

(19.11)

Свёртывание по второму индексу произведения тензора Т на вектор а

Для этого тензор Т умножается на вектор а справа

t a

+t a

2

+ t a

3

 

t

 

11 1

12

 

13

 

 

11

t21a1 +t22a2 + t23a3

 

= t21

t

a

+t

32

a

2

+ t

33

a

3

 

t

31

 

31 1

 

 

 

 

 

 

 

 

t

t

 

a

 

 

12

13

 

1

 

 

t22

t23

 

a2

 

= Tik ak

t32

 

 

 

 

 

t33

a3

 

 

Тогда произведение тензора Т на вектор а справа получается в виде

 

 

 

 

 

 

 

 

a' = T a =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

) (20.11)

= e

(t

 

a

+t

a

2

+t

 

a

3

) + e

2

(t

21

a +t

22

a

2

+t

23

a

3

) +e

3

(t

31

a

+ t

32

a

2

+t

33

a

3

1

 

 

11

1

 

12

 

13

 

 

 

 

1

 

 

 

 

 

1

 

 

 

 

 

Это получается также, если умножить сопряжённый тензор Т с на вектор а слева

(a a

 

 

 

)

t11

 

t21

t31

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

a

3

t

12

 

t

22

t

32

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

13

 

t

23

t

33

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= (a1t11 + a2t12 + a3t13

a1t21 + a2t22 + a3t23

a1t31 + a2t32 + a3t33 )

 

 

 

 

 

 

 

a' = a Tc

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= e1 (t11a1 + t12a2 + t13a3 ) + e2 (t21a1 + t22a2 + t23a3 ) +e3 (t31a1 + t32a2 + t33a3 )

Пример 2.11 Найти произведение aij x j тензора второго ранга aij , матрица которого в некотором базисе равна

 

2

0

3

 

 

 

5 1 2

 

,

(aij ) =

 

 

4

5

7

 

 

 

 

 

и тензора первого ранга xi (вектора), который в том же базисе имеет компоненты

(xi ) = ( 2 1 4 ),.

Решение.

52

a1 j x j = a11x1 + a12 x2 + a13 x3 = 2 2 + 0 1 + 3 4 =16 a2 j x j = a21x1 + a22 x2 + a23 x3 = 5 2 +1 1 + 2 4 =19 a3 j x j = a31x1 + a32 x2 + a33 x3 = 4 2 +5 1 + 7 4 = 41

Ответ: Вектор(16, 19, 41)

 

Умножение вектора а на тензор Т

 

a" = a T = (a e1) t1 + (a e2 ) t2 +a e3 t3

(21.11)

Координаты вектора a" получаются в виде

 

a1" = a1t11 + a2t21 + a3t31

 

a"2 = a1t12 + a2t22 + a3t32

(22.11)

a"3 = a1t13 + a2t23 + a3t33

 

Геометрическая интерпретация произведения вектора а на тензор Т

 

Произведение a T так составлено по векторам t1 , t2 t3 , как вектор a

составлен из

основных ортов i1, i2 i3 .

Для получения геометрической интерпретации ограничимся случаем двумерного пространства.

a = i1a1 + i2a2

a' = a1 t1 + a2 t2

Построим на взаимно перпендикулярных ортах i1 и i2 квадратную решётку из

растяжимых прутьев, соединённых шарнирами, как показано на рисунке 1.11. Теперь растянем стержни и повернём так, чтобы квадраты перешли в параллелограммы, как показано на рис.2.11. Тогда вектор а перейдет в новый вектор a’.

 

а

i2

 

O

i1

Рис.1.11 Вектор а в

 

системе i1, i2

Произведение тензоров

a

t2

O’ t1

Рис.2.11 Вектор а’ в

системе t1, t2

Пусть Aik и Bik - компоненты двух тензоров второго ранга. Их произведение будет иметь компоненты

Ciklm = Aik Blm

(23.11)

Числа Ciklm образуют тензор 4 –го ранга. Д о к а з а т е л ь с т в о:

53

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A'

 

=α

i'l

α

k 'm

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ik

 

 

 

 

lm

 

 

 

 

 

 

 

 

 

 

 

(24.11)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

'

 

=α

 

α

 

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i'l

k 'm

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ik

 

 

 

 

lm

 

 

 

 

 

 

 

 

 

 

 

 

 

C'

iklm

= A'

 

B'

= α

i'n

α

k'

p

α

k'r

α

 

A

 

B

r s

= α

i'n

α

k'

p

α

k'r

α

m's

C

n prs

(25.11)

 

 

 

ik

 

ik

 

 

 

 

 

m's

n p

 

 

 

 

 

 

 

 

Определение

5.11

Операция

образования

компонент

 

Ciklm

 

называется

внешним

умножением тензоров Aik и Bik .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Замечание 8. 11

Тензорное произведение некоммутативно.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ciklm = Aik Blm Clm ik

= Alm Bik

 

 

 

 

 

 

 

 

 

 

(26.11)

Определение 6.11 Произведением нескольких тензоров называется тензор, компоненты

которого равны

произведению компонент сомножителей.

При этом ранг произведения

равен сумме рангов сомножителей.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Определение 7.11 Внешним произведением тензоров произвольного ранга называется

новый тензор, у которого компоненты образованы умножением каждой компоненты

одного тензора на каждую компоненту второго.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a) aib j = Tij

 

b) ai Flk

= αilk

 

 

c) DijTlm

= Φijlm

 

e) εijkυm =θijkm

 

 

(27.11)

Замечание 9.11 Как видно из этих примеров, внешние произведения получаются путём

написания умножения тензоров друг за другом.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Замечание 10.11

 

Внешнее произведение двух векторов образует одну диаду.

 

 

§ 12.

Главные значения и главные направления тензора второго ранга.

Основные понятия

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рассмотрим произвольный тензор 2-го ранга Τik

. Если этот тензор умножить на вектор

Ak

и произвести

 

свёртывание

по

 

индексу вектора и одному индексу тензора, то в

результате получим некоторый вектор Bi

с компонентами

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Bi = Τik Ak

 

 

 

 

 

 

 

 

 

 

 

 

Ak ,

 

(1.12)

 

Тензор Τik , будучи умножен скалярно на некоторый вектор

преобразует его в

новый вектор в том смысле, что из компонент вектора

Ak

определённым действием

получаются компоненты другого вектора – вектора

 

 

 

 

 

 

 

 

 

 

 

 

 

Bi . Вектор Bi

вообще отличен от Ak

по величине и

 

 

 

 

 

 

 

 

 

'

 

x3

p1

направлению. Таким образом,

 

тензор

 

при

 

 

 

 

 

 

 

 

x3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n1

умножении на вектор изменяет длину этого

 

p3

n3

 

 

 

 

 

 

 

вектора и поворачивает его оси.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Задача заключается в том, чтобы найти для

 

 

 

 

 

 

 

 

 

 

 

 

x2'

данного тензора

 

Τik такие векторы

Ak ,

которые

 

 

 

 

 

 

 

 

 

 

 

 

бы не поворачивались этим тензором, а только

 

x1

 

 

 

 

 

 

 

 

 

x2

изменяли длину (рис. 4.12).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n2

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

'

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p2

 

 

 

 

 

 

 

 

Τik Ak

= λAi ,

 

 

 

 

 

(2.12)

 

 

 

 

 

x1

 

 

 

 

 

 

 

 

где λ - скаляр.

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 4.12

Главные оси тензора

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

напряжений

pik

в точке М. На

 

Физический

смысл

этой

 

задачи

виден

на

 

 

 

 

 

 

 

 

площадках, перпендикулярным к

некоторых примерах.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

осям

x' , x'

, x'

, касательные

Задача 1.12 Напряжение на площадке с нормалью

 

 

 

 

 

 

 

 

 

1

 

2

 

3

 

 

n равно

 

 

 

 

 

pnk

= pik ni

 

 

 

 

 

(3.12)

 

 

 

 

 

 

напряжения равны нулю.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

54

причём вообще вектор pn не параллелен орту n , то есть, на каждой площадке есть как

нормальные, так и касательные напряжения.

Интерес представляют такие площадки, на которых есть только нормальные напряжения, а касательные равны нулю. Для этих площадок

pn || n или λn = pn = pi ni

Тогда ориентация этих площадок, которую дают орты n , определится из системы уравнений

λnk = pik ni .

(4.12)

Задача 2.12 Если диэлектрические свойства среды определяются тензором Tik , то

возникает вопрос, как следует направить электрическое поле Е, чтобы электрическая индукция D была направлена по вектору напряжённости Е?

Общий вид зависимости D и Е имеет линейный характер

Di = Tik Ek

(5.12)

Поставленная задача требует отыскания векторов Ei , удовлетворяющих уравнениям

λ Ei = Tik Ek

(6.12)

Определение 1.12 Если существуют для тензора Tik

векторы Ak , удовлетворяющие

уравнениям

 

λ Ai = Tik Ak ,

(7.12)

то направления, определяемые этими векторами

Ak , называются главными

(собственными) направлениями тензора Tik .

Определение 2.12 Оси главных направлений называются главными осями тензора. Определение 3.12 Значения компонент тензора в координатной системе главных осей называются главными значениями.

Определение главных направлений и главных значений тензора Tik

Согласно (7.12) компоненты вектора А, определяющие оси тензора Tik , удовлетворяют системе трёх уравнений:

Tik Ak

λ Ai = (Tik λδik ) Ak = 0

(8.12)

или

 

 

 

(T11 λ) A1 +

T12 A2 +

T13 A3 = 0

 

T21 A1 + (T22 λ) A2 +

T23 A3 = 0

(9.12)

T31 A1 +

T32 A2 + (T33 λ) A3 = 0

 

Эта однородная система служит

для определения A1 , A2 , A3 .

При этом, ищется

отличное от нуля, или, нетривиальное решение этой системы. Однородная система уравнений имеет нетривиальное решение только в том случае, когда её определитель равен нулю. Таким образом, для определения главных значений имеется уравнение

(T11 λ)

T12

T13

 

 

 

 

T21

(T22 λ)

T23

= 0

(10.12)

T31

T32

(T33 λ)

 

 

Уравнение (10.12) представляет собой кубическое уравнение относительно λ .

Определение 4.12 Уравнение (10.12) называется характеристическим уравнением

тензора Tik .

55

Замечание 1.12 Корни кубического уравнения (10.12) в общем случае могут быть не все действительными, и тогда из этого уравнения нельзя найти главные направления тензора Tik .

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

Например, если известны ковариантные компоненты тензора Tik , то уравнение, определяющее собственные векторы А , имеет вид:

 

 

 

 

λ A = T

Ak .

 

(11.12)

 

 

 

 

i

ik

 

 

 

 

Отсюда в силу

A = g

ik

Ak

получим

систему линейных

однородных уравнений

 

i

 

 

 

 

 

 

 

относительно Ak

вида

 

 

 

 

 

) Ak

 

 

 

 

 

 

(T

λ g

ik

= 0 .

(12.12)

 

 

 

 

ik

 

 

 

 

Однако, для того, чтобы привести эту систему

к виду (9.12), необходимо пользоваться

смешанными компонентами

тензора

 

 

g i

= δ i

. Умножив

предыдущие уравнения

 

 

 

 

 

 

 

 

.k

.k

 

 

 

 

 

 

 

(i =1,2,3) на g il , пронумеровав по i , получим в силу уже известных формул

Aik = g il g km A ,

A

= g

il

g

km

Alm ,

A

= g

A l

= g

il

Al

 

 

 

lm

ik

 

 

 

 

ik

 

kl i

 

k '

 

 

Al k = g kl Ail ,

 

Al k = gil Alk ,

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

Ai

= g il A ,

Ai

= g

 

Ail ,

 

Aik = g il A k

= g kl Ai

.

 

 

 

kl

 

 

 

 

k

lk

k

 

 

 

 

 

 

l

 

l

 

 

 

 

следующие уравнения

(T i k

λ g i k ) Ak = 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(13.12)

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

(T 1 1 λ)

T 1 2

T 1 3

 

 

T 2 1

(T 2 2 λ)

T 2 3

= 0

(14.12)

T 3 1

T 3 2

(T 3 3 λ)

 

 

Замечание 3.12 Обычно рассматриваются только симметричные тензоры, потому что у них корни характеристических уравнений всегда действительные.

Остановимся на рассмотрении симметричных тензоров второго ранга, отнесённых к прямоугольным декартовым системам координат, так что Ti k = Tk i . В этом случае все

корни λ1 , λ2 , λ3 характеристического уравнения (13.12) вещественные

(действительные).

Действительно, пусть λ - какой-нибудь из корней уравнения (13.12) и пусть ему отвечают в силу системы уравнений (9.12) какие-то величины Ai , вообще комплексные.

Тогда умножив каждое (i =1,2,3) из тождеств

λ Ai = Tik Ak

(15.12)

на величины Ai , комплексно сопряжённые с Ai , и просуммировав по i , получим

 

 

 

 

 

 

 

λ Ai

 

 

= Tik Ak

 

 

 

 

 

 

(16.12)

 

 

 

 

Ai

Ai

 

 

Так как Ti k = Tk i , то

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T A

 

=

1

(T A

 

+T A

 

)=

1

(T A

 

+T

 

A )=

1

T

(A

 

+

 

A ).

A

A

A

A

A

A

A

2

2

 

ik k i

ik k i ik k i

ik k i k i k i

2 ik

k i k i

56

Отсюда видно, что сумма Tik Ak Ai вещественна, так как все Ti k вещественны и выражение в скобках вещественно. Вспомним, что произведение комплексного числа на сопряжённое

z z = (x + iy) (x iy) = x2 + y2

равно сумме квадратов вещественных чисел. Поскольку

A

 

=

 

A

 

2 +

 

A

 

2 +

 

A

 

2 - тоже вещественная величина, то из (15.12) следует, что корень

A

 

 

 

 

 

 

i

i

 

 

1

 

 

 

 

2

 

 

 

3

 

 

 

 

 

 

 

 

 

λ вещественная величина. При этом, конечно, все компоненты Ak

тоже вещественны (это

следует из λ Ai = Tik Ak ).

 

 

 

 

 

 

 

Итак, λ1 , λ2 , λ3 - вещественные числа.

 

 

 

 

 

Заметим,

что если

 

Ti k - ковариантные компоненты тензора, то,

исходя из равенства

λ A = T

 

Ak ,

 

 

умножая

 

его на

 

 

, суммируя по i и используя

T

 

= T

, аналогично

 

 

 

 

 

A

 

 

i

 

 

ik

 

 

 

 

 

 

 

 

 

 

 

 

 

i

i k

k i

 

предыдущему получим вещественность характеристического уравнения (13.12).

Главные значения и главные направления симметричных тензоров второго ранга

В дальнейшем будут рассмотрены только симметричные тензоры с действительными компонентами. Это несколько проще в математическом отношении, так как тензоры, важные для механики сплошной среды, обычно симметричны, то жертвуя немногим, целесообразно принять такое ограничение.

Для каждого симметричного тензора Ti j , заданного в некоторой точке пространства, и для каждого направления в этой точке (характеризуемого единичным вектором ni ) существует вектор, определяемый внутренним произведением (см. определение 6.11)

υi = Ti j n j

(17.12)

Здесь Ti j можно рассматривать как линейный векторный оператор, который ставит в

соответствие направлению ni вектор υi . Если направление таково, что

вектор υi

параллелен ni ,

то указанное внутреннее произведение выражается

скаляром,

умноженным на ni . В этом случае (так как υi = λ ni ) получается

 

 

Ti j n j = λ ni

(18.12)

и направление ni

называется главным направлением или главной осью тензора Ti j . С

помощью тождества ni = δi j n j соотношению (18.12) можно придать форму

 

 

(Ti j λδi j )n j = 0,

(19.12)

которое представляет систему трёх уравнений для четырёх неизвестных ni и λ ,

соответствующих каждому главному направлению. В развёрнутой записи система, которую следует решить, имеет вид

(T11 λ) n1 + T12 n2 + T13 n3 = 0

 

T21 n1 + (T22 λ) n2 + T23 n3 = 0

(20.12)

T31 n1 + T32 n2 + (T33 λ) n3 = 0

Это однородная система уравнений, поэтому при любом λ существует тривиальное решение ni = 0 , но наша цель состоит в том, чтобы получить нетривиальное решение, то есть отличное от нуля. Кроме того, не теряя общности, можно ограничится только

решениями, для которых

n

n

i

=1

(это значит, что n2

+ n2

+ n2

=1, что справедливо, так

 

i

 

 

1

2

3

 

как с самого начала вектор

 

ni предполагался единичным).

Для того, чтобы система

(19.12) или, что то же самое (20.12), имела нетривиальное решение, определитель из коэффициентов должен быть равен нулю, что можно записать так:

57

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ti j λδi

j

= 0 .

 

 

 

(21.12)

В развёрнутом виде это кубическое уравнение относительно λ :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λ3 I

T

λ2

+ II

T

λ III

= 0 ,

(22.12)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

которое

является

 

характеристическим

 

уравнением

тензора

Ti j , а его скалярные

коэффициенты соответственно равны

 

 

 

 

 

 

 

 

 

 

 

 

 

IT =Tii

=T11 +T22 +T33 (= λ1 +λ2 +λ3 ) = inv

 

 

 

 

 

 

(23.12)

II =

1

(T T

jj

 

 

T

T

j i

)

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

2

 

ii

 

 

 

i j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(24.12)

=

 

T22

 

 

T23

 

+

 

T11

T21

 

+

 

T11

T31

 

(= λ λ

2

+ λ λ

3

+ λ

2

λ ) = inv

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T32

 

 

T33

 

 

 

 

 

T12

T22

 

 

 

T13

T33

 

1

 

 

1

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T11 T12

 

T13

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

IIIT =detTi j

=

T21

 

T22

 

T23

(=λ1λ2λ3 ) =inv

 

 

 

 

 

 

(25.12)

 

 

 

 

 

 

 

 

 

 

 

 

T31

 

T32

 

T33

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

называются соответственно первым, вторым и третьим инвариантами тензора Ti j . Три корня кубического уравнения (20.12), обозначенные λ1 , λ2 , λ3 , называются главными значениями тензора Ti j . У симметричного тензора с действительными компонентами

главные значения действительны; если все они различны, то три главных направления взаимно ортогональны. В главных осях таблица из компонент тензора приводится к диагональной форме

 

λ

0

0

 

 

1

λ2

 

 

 

T =

0

0

 

(26.12)

 

0

0

 

 

 

 

λ3

 

Замечание 4.12 Если λ1 = λ2 , то диагональный вид тензора не зависит от выбора осей, соответствующих λ1 и λ2 , и нужно установить только главную ось, соответствующую

λ3 .

Замечание 5.12 Если все главные значения равны, то любое направление является

главным.

Если главные значения упорядочены, то их

принято обозначать λI , λII , λIII и

располагать в порядке убывания: λI

> λII

> λIII .

 

 

 

 

 

 

 

Преобразование системы Ox x

2

x

3

к системе главных осей Ox* x* x*

даётся элементами

таблицы

 

1

 

 

 

 

 

 

 

 

 

1

2

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x1

 

 

 

 

 

 

x2

 

x3

 

 

 

 

 

 

x*

a

 

= n

(1)

 

a

 

= n(1)

a

 

= n(1)

 

 

 

 

 

1

11

 

 

1

 

 

12

2

13

3

 

 

 

 

 

x*

a

21

= n(2)

 

a

22

= n(2)

a

23

= n(2)

 

 

 

 

 

2

 

 

 

 

1

 

 

 

2

 

3

 

 

 

 

 

x*

a

31

= n(3)

 

a

32

= n(3)

a

33

= n(3)

 

 

 

 

 

3

 

 

 

 

1

 

 

 

2

 

3

 

 

 

 

где ni( j) - направляющие косинусы j - того главного направления.

Пример 1.12 Найти главные направления и главные значения декартова тензора Т второго ранга, который представлен матрицей

58

 

 

3

1

0

[T

]= −1

3

0

ik

 

 

 

 

 

 

 

 

0

 

0

1

 

 

 

 

 

Решение. Для определения главных значений необходимо решить уравнение

3 λ

1

0

 

 

 

 

 

 

1

3 λ

0

 

= (1 λ)[(3 λ)

2

1]

= 0

 

 

 

 

0

0

1

 

 

 

 

 

 

λ

 

 

 

 

Это кубическое уравнение

λ3 7λ2 +14λ 8 = 0

Но совершенно очевидно, что один корень известен λ1 =1. Поделив кубическое уравнение на λ 1, получим

λ2 6λ +8 = 0

Это квадратное уравнение имеет два корня

λ2,3 = −

(6)

±

(6) 2

8 = 3 ± 9 8 = 3 ±1

 

 

 

 

2

 

2

 

λ2

= 2,

 

λ3 = 4

 

 

Замечание 6.12. Здесь использовано квадратное уравнение вида x2 + px + q = 0 , решение

которого равно x1,2 = −

p

±

 

p

2

 

 

2

 

 

q

 

 

 

 

2

 

 

 

Таким образом, найдены все три главные значения λ1 =1 λ2 = 2,

λ3 = 4

 

 

 

 

Определение главных направлений.

 

 

 

 

 

 

 

 

 

 

 

 

 

а) Пусть ni1

-

компоненты единичного вектора

главного

направления ,

соответствующего λ1 =1. Тогда два первых уравнения системы

 

 

 

 

 

 

 

 

(3 λ) n1

n1

+ 0 n31 =

0

 

2 n1

n1

+ 0 n31

= 0

 

 

 

 

 

 

 

 

 

1

2

 

 

 

 

1

 

2

 

 

 

 

 

 

 

 

 

 

n11 + (3 λ) n12 + 0 n31

= 0

n11 + 2 n12 + 0 n31 = 0

 

 

 

 

 

 

0 n1 + 0 n1 + (1 λ) n1

= 0

0 n1 + 0 n1 + 0 n1 = 0

 

 

 

 

 

 

 

1

 

2

3

 

 

 

 

1

 

2

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дают 2n1

n1

= 0,

n1 + 2n1 = 0,

откуда n1

= n1

= 0,

а

из

условия

n

n

i

=1 получим

 

1

2

 

 

1

2

 

1

 

2

 

 

 

 

 

 

i

 

 

n31 = ±1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

б) Для λ2

= 2,

система уравнений (4) имеет вид

 

 

 

 

 

 

 

 

 

 

 

 

 

(3 λ) n2

n2

+ 0 n32 = 0

 

n2

n2

+ 0 n32

= 0

 

 

 

 

 

 

 

 

 

1

2

 

 

 

 

1

 

2

 

 

 

 

 

 

 

 

 

 

n12 + (3 λ) n22 + 0 n32

= 0 n12 +n22 + 0 n32 = 0

 

 

 

 

 

 

0 n2 + 0 n2 + (1 λ) n2

= 0

0 n2

+ 0 n2 n2 = 0

 

 

 

 

 

 

 

1

 

2

3

 

 

 

 

1

 

2

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

даёт n2

n2 =

0, n2

+ n2

= 0, и n2

= 0. Таким образом,

n2

= 0,

а n2

= n2

= ±1 2 , так

1

 

2

 

1

2

 

3

 

 

 

 

 

 

3

 

1

 

2

 

как ni ni

=1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

в) Для λ3

= 4 из системы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3 λ) n3

n3

+ 0 n33 =

0

 

n3

n3 + 0 n33 = 0

 

 

 

 

 

 

 

 

1

2

 

 

 

 

 

1

 

2

 

 

 

 

 

 

 

 

 

n13 + (3 λ) n32 + 0 n33

= 0 n13 n32 + 0 n33 = 0

 

 

 

 

 

 

0 n3 + 0 n3 + (1 λ) n3

= 0

0 n3

+ 0 n3 3n3

= 0

 

 

 

 

 

 

 

1

 

2

3

 

 

 

 

1

 

 

2

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

59

получаем n3

n3

= 0 ,

n3

n3

= 0

3 n3

= 0 . Таким образом, n3

= 0, а n3

= n3

= m1 2 .

1

2

 

1

2

 

3

3

1

2

 

Ориентация главных осей xi' относительно исходной системы xi определяются направляющими косинусами, которые даны в следующей таблице

 

x1

 

x2

 

x3

x1'

0

 

0

 

±1

x2'

±1

2

±1

2

0

x3'

m1

2

±1

2

0

Ответ Из таблицы видно, что матрица преобразования такова:

0

 

0

 

±1

A = ±1

2

±1

2

0

±1

2

±1

2

0

Степени тензора второго ранга. Соотношения Гамильтона –Кэли

Непосредственным матричным умножением квадрат тензора Ti j получается как внутреннее произведение Ti k Tk j , куб – как произведение Ti k Tk mTm j и т.д. Таким образом,

если Ti j представлен в диагональной форме (26.12), то

n -я степень этого тензора (и

соответствующей матрицы) даётся формулой

 

 

 

 

 

 

 

 

λn

0

0

 

 

λn

0

0

 

 

(T)n =

 

1

λn

 

 

 

1

 

 

 

 

 

0

0

 

или T n =

0

λn

0

 

(27.12)

 

 

 

2

n

 

 

 

2

 

 

 

 

 

0 0 λ

 

 

 

 

n

 

 

 

 

3

 

 

0 0 λ3

 

 

 

 

 

 

 

 

 

 

 

 

 

Замечание 7.12 Сравнение (26.12) и (27.12)

 

показывает, что тензор Ti j

и все его целые

степени имеют одни и те же главные оси.

Все главные значения удовлетворяют уравнению (22.12), а матрица T n имеет диагональный вид (27.12), поэтому сам тензор Т будет удовлетворять уравнению (22.12). Таким образом,

T 3 I

T

T 2 + II

T

T III

T

E = 0

(28.12)

 

 

 

 

 

где Е – единичная матрица. и

Определение 5.12 Соотношение (28.12) называется соотношением Гамильтона – Кэли.

Если умножить каждый член соотношения (28.12) на Т по правилу перемножения матриц, то получается равенство

T

4 = I

T

T 3

II

T

T

2 + III

 

T

 

 

 

(29.12)

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

Подставляя T 3 из (28.12), получают

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T 4 = (I2

II

T

)T 2

+ (III

I

T

II

T

)T + I

III

T

E

(30.12)

T

 

 

 

 

 

T

 

 

 

T

 

 

Замечание 8.12 Продолжая эту процедуру, можно получить все целые положительные степени Т в виде линейных комбинаций Т2, Т и Е.

§ 13. Ковариантное дифференцирование тензоров

Ковариантный дифференциал тензора

Рассмотрим выражение дифференциала вектора а через дифференциалы его

компонент. В декартовой системе координат имеем

 

da = d (al il ) = il dal .

(1.13)

60

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]