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

Курс механики сплошных сред

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

новка из (1, 2, 3), и равно —1, если перестановка нечетная; во всех других случаях ziJH равно нулю.

Введем обозначения:

I 6"

 

«1г 1

 

 

bq }■

. ^кр

b q

e

j

l

6/7p

®ijkpqr

 

V

 

UPQJ h

ip

IQ

в которых правые части являются определителями матриц, все эле­ менты которых являются символами Кронекера. Приведем разложения определителей по элементам:

®Чрч ~

 

®IJkpqr ~ ^ifijkp q Н" ^jfiklpq Н" ^kr^i/pq-

 

Теперь легко проверить соотношения

 

6,7 = 3, 8l/py = 2blp, ^i/kpqk — ^tjpqf

(41)

в которых, конечно, по повторяющимся индексам в левых частях следует производить суммирование.

Докажем тождество

^UbPpqr =

®i/kpqr•

(42)

Обе части равенства равны 1

при i = p= 1, j = q = 2,

k = r = 3.

Они продолжают оставаться равными, если осуществить перестановку любых двух из индексов t, /, k или из р, q, г. Таким образом, обе части равенства равны, когда i, /, k и р, q, г являются произ­ вольными перестановками из 1, 2, 3. В противном случае они, оче­ видно, оба равны нулю, что и требовалось доказать.

Используя (41), получаем .тождества

 

e,/fteMft = 6(/,6/9- 6 /(J6^ )

(43)

EiJkRp/k — 26,-р, ei/k&ijk — 6 . I

Видно, что обозначение г1/к удобно для сокращенной записи некоторых формул. Вот один из примеров. Если А —матрица, состав­ ленная из элементов А1)у то легко проверить соотношения:

Rl/kAlpA/qAkr = Bpqr det(A),

det А = -g- tijffipqrAipAfqAkr =

^i/kpqr^ip^/q^krt

(44)

где det (А) — определитель матрицы А.

 

единичный

П 1.3.2. Тензорное представление. Выше был введен

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

Легко видеть, что b,/pq и 6,7ftp(jr — соответственно компоненты тензоров четвертого и шестого рангов в любом ортонормированном базисе, называемые иногда тензорами антисимметризации прост-

ю*

ранства Е , которые после выбора базиса можно записать в виде: bijpqei ® e J® e p® e q,

ЬцьрчгЪ® eJ ® e lt® e p ® e q® ег.

Действительно, эти тензоры строятся из единичного тензора с помощью операций типа: тензорное умножение на единичный тен­ зор, перестановка, линейная комбинация при постоянных коэффи­ циентах (равных +1 или —1).

Название тензор антисимметризации объясняется на следующем примере. Компоненты

(45)

с точностью до множителя 1/2, равны компонентам тензора, полу­ ченного умножением с двойной сверткой первого тензора антисим­ метризации и, тензора, полученного перестановкой индексов тензора второго ранга Т. Но этот тензор как раз и совпадает с антисиммет­ ричной составляющей тензора Т (П1.4.2), что можно проверить, определив по формуле (43) его компоненты.

Представляет интерес выяснить, как ведет себя тензор

0 = ziJke l ® e J ® e k

(46)

при замене ортонормированного базиса. В базисе е\ его компоненты

по

(36)

и (44)

 

 

 

 

 

 

 

 

 

ppip 4 j prk*iih = bpv det(P).

 

 

±1

Так

как

матрица Р —ортогональная,

то

ее определитель

равен

и снова

получаем грдп

если

det(P) = + l (в этом случае

гово­

рят, что базис ё[ имеет ту

же

ориентацию,

что

и базис et);

знак

минус

перед

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

что базис

е\

имеет

противоположную

ориентацию. Таким образом, величина г^к представляет компоненты тензора во всех ортонормированных базисах, имеющих одинаковую ориентацию —это тензор ориентации этих базисов. Но в базисе обратной ориентации компоненты этого же тензора будут противо­ положными.

Этот же результат можно выразить по-другому через понятие псевдотензора или аксиального тензора, понимая под ним математический объект, который во всех ортонормированных базисах одной определенной ориентации представлен одинаково, а во всех базисах, ориентация которых не совпадает с зафиксиро­ ванной, имеет противоположные компоненты. Таким образом, для любого орто­ нормированного базиса функция г рдг определяет компоненты некоторого псевдо­

тензора. Оперируя с тензорами или псевдотензорами, как в П1.2.4, получим либо тензоры либо псевдотензоры.

Если зафиксировать ориентацию пространства, иными словами, если условиться иметь в виду только такую замену базисов, при которой ориентация остается неизменной, то тогда можно не делать различия между тензорами и псевдотензорами. Именно так и будет в большинстве случаев.

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

ГМ.3.3. Приложения, а) Пусть заданы три вектора А, В, С.

Образуем тензорное произведение С ® В ф А и умножим с тройной

сверткой тензор О на полученное произведение.

При использовании

компонентов А, В и С в базисе е, результат запишем в виде

еиЛ В ;С к.

(47)

Получили, таким образом, псевдоскаляр, который совпадает со смешанным произведением трех векторов А, В и С, равным опреде­ лителю матрицы компонентов этих векторов. Выражение (47) имеет вполне определенное значение во всех базисах с данной ориентацией и противоположнбе значение в базисах обратной ориентации.

б) Пусть

дан антисимметричный тензор второго

ранга * Я

ij = — Qji).

Формула

 

 

 

со,- = -у eiJkQkj

 

(48)

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

аксиального

вектора ю.

 

 

Умно­

Обратим внимание на то, что ю1= й 32, со2 = Я18, (о3 = Я22.

жив обе части равенства (48) на epqi и воспользовавшись

тож­

дествами (43),

получим

 

 

 

= ~2 ^ P q A ijtQ h j ~ ~2 ^ p q fijk fo k j = ®qp

 

 

ИЛИ

 

 

 

 

О//*»//*0»*-

 

(49)

Таким образом, видно, что антисимметричному тензору можно поставить в соответствие аксиальный вектор и, наоборот, аксиаль­ ному вектору можно поставить в соответствие антисимметричный тензор.

в) Формула У &ijk^kj определяет аксиальный вектор, соответст­

вующий

антисимметричной составляющей тензора второго ранга Т

В частном случае, когда

тензор — это тензорное

произведение

2 В 0 Л ,

получаем векторное

произведение А А В, компоненты ко­

торого

 

 

 

 

 

*-ijhAiBk-

(50)

П1.4. ТЕНЗОРЫ ВТОРОГО РАНГА

П1.4.1. Предварительные замечания. Если не будет специ­ ально оговорено противное, речь пойдет только о тензорах второго ранга, которые по-прежнему будем обозначать жирными прописны­ ми буквами А, В . .. В ортонормированном базисе е,- тензоры пред­ ставляются матрицами А, В . .. или компонентами тензоров Ay, By.

* В этом примере встречается элементарное определение, данное в П1.4.2.

В другом ортонормированном базисе е\ эти же тензоры будут представлены другими матрицами А*, В* или компонентами тент зора A*ij, B*Cj . . . Если Р —ортогональная матрица перехода от одного базиса к другому, то заданные тензоры связаны соотношениями, полученными в П1.2.3. В частности,

А* = РАРТ, А = РТА*Р.

(51)

Любое соотношение между тензорами А, В ... в любом ортонор­ мированном базисе выражается соотношением между квадратными

матрицами, представляющими эти

тензоры, и

соотношение между

матрицами имеет один и тот же

вид

во

всех

базисах,

т.

е. если

его требуется получить в базисе

е \ ,

то

в соотношении

в

базисе et

следует заменить А на А*, а В на В*.

Обратно, матрицы А и А* соответствуют одному и тому же

тензору

в базисах е\ и е ь

если эти

матрицы связаны формула­

ми (51).

Если некоторое

соотношение,

связывающее матрицы А,

В . . . , остается инвариантным относительно замены базиса (51), а ортогональная матрица Р —произвольна, то говорят, что данное соотношение является внутренним или объективным; в этом слу­ чае оно может быть интерпретировано непосредственно как соот­ ношение между тензорами Л, В . . . , представленными в базисе е{ матрицами А, В . . .

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

П1.4.2. Простейшие определения. Некоторые определения, отно­ сящиеся к матрицам, элементарные свойства которых полагают

известными,

могут быть

сразу же

отнесены

к тензорам второго

ранга.

 

 

 

 

 

 

Выше было дано определение тензора, полученного перестанов­

кой

индексов

из тензора,

например,

Л —транспонированного

тен­

зора,

который обозначим

Лт В любом ортонормированном

базисе

этот

тензор

представлен

транспонированной

матрицей Ат ,

т.

е.

матрицей, полученной транспонированием матрицы А, относящейся

к тензору Л в рассматриваемом базисе (Ajj = A/l).

Тензор называется симметричным, если он равен своему транс­ понированному тензору (т. е. не меняется при перестановке ин­ дексов).

Симметричные тензоры образуют шестимерное векторное под­ пространство в векторном пространстве девяти измерений тензоров второго ранга.

Тензор называется антисимметричным, если он равен транс­ понированному с обратным знаком.

Четная, или симметричная, составляющая и нечетная, или ан­ тисимметричная, составляющая тензора Л определяются соответст­ венно формулами

AS = ± (A + AT), Ла = ± 0 4 - Л Т ),

(52)

так что

A

АзЛ~Аа, А Т As А а%

(53)

Эти же соотношения

верны и для матриц в

любом ортонорми­

рованием базисе и дают каноническое разложение тензора на анти­ симметричную и симметричную составляющие.

Произведением двух тензоров А

и В

называют

произведение о

однократной

сверткой обоих тензоров —тензор

С.

В любом базисе

et тензор С

представлен матрицей

С,

равной

произведению мат­

рицы А на матрицу В*:

 

 

 

 

 

С = АВ, С/у = AikBkJ.

 

(54)

Тензор, получаемый транспонированием произведения, пред­ ставляет собой произведение в обратном порядке транспонирован­ ных тензоров; это верно как в отношении тензоров, так и в отно­ шении матриц.

Заметим, что произведение двух симметричных тензоров не будет симметрич­ ным тензором.

Определителем тензора А называется определитель

det (А) = det (Л/у).

(55)

Для обоснования такого определения (которое в отличие от предыдущих дается формулой, относящейся к матрице) проверим, что матрица А*, представляющая тензор А в другом базисе, имеет тот же определитель. Это утверждение следует из формулы (51), так как матрица и транспонированная матрица имеют один и тот же определитель, и в силу правила перемножения определителей имеем

det (А*) = det (A) (det (Р))2 = det (А),

поскольку Р — ортогональная матрица и, следовательно, согласно

(23) (det (Р))а= 1.

Говорят, что тензор невырожденный, если его определитель отличен от нуля; тензор называют вырожденным, если его опре­

делитель равен нулю.

Тензоры А и В называются взаимно обратными или взаимными,

если тензоры-произведения А В и В Л равны каждый

единичному

тензору.

 

В базисе е имеем тогда

 

А В= 1, ВА = 1,

(56)

эти равенства являются объективными, что легко проверяется: А*В* = РАРТРВРТ = РАВРТ = РРТ = 1.

В этом определении тензоры А и В предполагаются невырож­ денными, так как

det (A) det (В) = 1.

♦Напомним, что произведения двух матриц записываются без знака умноже­ ния между буквами.

Если А — заданный невырожденный тензор, то первого из ра­ венств (56) достаточно для однозначного определения тензора В.

Действительно,

девять

уравнений

AikBkJ=

при последователь­

ной подстановке

/ = 1,

/ = 2, /' = 3

представляют собой три системы

из трех линейных уравнений, в которых неизвестными являются компоненты вектора-столбца матрицы В. Определитель det (А) в каждой системе не равен нулю, так как тензор А — невырожден­ ный. Определенная таким образом матрица удовлетворяет тогда второму из соотношений (56), так как после умножения первого равенства справа на А, имеем

АВА = А,

A (BA— 1) = 0.

 

В связи с тем что det (А) =^0,

то, рассуждая как и

раньше, на­

ходим, что множитель справа

(BA— 1) тождественно

равен нулю,

что и требовалось доказать.

 

обозначается

Тензор, обратный невырожденному тензору А,

через А-1 и представляется матрицей А-1, обратной

матрице А.

Объект, обратный произведению А В, обозначается через В-1А-1. Такие определения дают возможность построить последователь­ ные степени тензора А; они могут быть представлены как степени матриц А" (п — положительное или отрицательное, А° = 1). Если А—симметричный тензор, то степени от А также симметричные

тензоры.

П1.4.3. Инварианты тензора второго ранга. Определение. Вся­ кая скалярная функция тензора второго ранга называется инва­ риантом этого тензора. Иными словами, инвариант определяется отображением векторного пространства тензоров второго ранга на множество действительных чисел.

Пусть дана матрица Т, представляющая тензор Т в ортонорми­ рованием базисе б/. Инвариант f является скалярной функцией Т,

которая должна принимать одно и то же значение

при замене Т

на Т* — матрицу,

представляющую Т в другом ортонормированном

базисе е]. В силу

(35) функция [ удовлетворяет тождеству

 

/ (РТРТ) = / (Т)

(57)

при любой ортонормированной матрице Р.

удовлетворяю­

Матричная функция со скалярными значениями,

щая тождеству (57), называется изотропной. Можно сказать также, что, будучи представленной в некотором базисе как функция ком­

понентов Т lJy эта функция не должна менять своего

значения

при

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

 

Выше отмечалось, что след тензора Т

является

инвариантом

(40). Точно так же след t r { r ,z} (п — целое

положительное или

от­

рицательное число) является инвариантом тензора Т.

След тензора второго ранга также является следом матрицы, представляющей тензор в ортонормированном базисе.

Отметим, в частности, три инварианта:

7 ,= tr{ T } = :r„, )

Т\\ — tr {Т2} =» Т ijTJh

Tm = tr {Т8} = Г tjT jhTkl. ,

Выше [см. определение (55)] отмечено также, что определитель detT является инвариантом тензора Т. При любом фиксированном значении X det(T — АЛ) также является инвариантом тензора Т. Таким образом, у многочлена

Р(X) = det (Т- XI) = det (Ти-Х Ь и),

(69)

называемого характеристическим полиномом тензора Т, все коэф­ фициенты — инварианты тензора Т (их называют элементарными инвариантами). Запишем этот полином в виде

 

 

Р(Х) = Т ш - Х Т и + Х^Т1- Х \

(60)

Раскрывая определитель (59) [можно, например,

использовать

(44)

и соотношения (41)], найдем:

 

 

 

Т I = &ijTи = Тih

 

^11

jq

2 filp ^fq ^lq^/p)T<p'T/q 2 ^ цТ Ц

T y T j t),

 

III

: -Q^UkpqrTipTjqTkr — det (T).

 

Инварианты (61) могут быть выражены как функции инвари­ антов (58). Для этой цели воспользуемся для Тш разложением определителя 6{jkpqr, который находится по известным правилам вычисления определителей третьего ранга:

^ijkpqr =* p^jq^kr + bJpbhqb lr - f bkpblqbjr b irbJqbkp §Jrbkqbip— §krbiqbjp .

И теперь нужные соотношения получаются автоматически:

Тг - Г„ Т и = \ - 7 „ ) , Гш = | (Т \~ 3Т.Г,, + 2ГШ). (62)

П1.4.4. Собственные векторы и значения симметричного тен­ зора второго ранга, а) Общее определение. Ненулевой вектор а на­ зывают правым собственным вектором тензора Т, если вектор Т а

коллинеарен а. Ненулевой вектор b называют левым собственным вектором тензора Т, если вектор Ь Т коллинеарен Ь.

Всякий вектор, пропорциональный собственному вектору, также является собственным, а направление собственного вектора—глав­ ным направлением тензора Т (правым или левым).

По определению, векторы а и b должны удовлетворять урав­ нениям

(7 —ц/) а = 0, Ь (Т— v/) = 0.

Каждое из уравнений, записанных в базисе eh дает линейную однородную систему, которая должна выполняться при подстановке

компонентов векторов а и Ь. Так как векторы а и Ь отличны от нуля, то ц н v будут корнями характеристического полинома Р (Я).

Эти корни называют собственными значениями

тензора

Т и обо­

значают

ts, /3.

 

 

 

 

рассмот­

б) Случай симметричного тензора. Далее ограничимся

рением только симметричных тензоров. В этом

случае

любой пра­

вый собственный вектор будет одновременно и

левым

собственным

вектором

и наоборот. Кроме того,

можно доказать

следующий ре­

зультат.

 

симметричного

тензора —всегда

Теорема. Собственные значения

действительны. Кроме того, существует по меньшей мере один орто-

нормированный базис е{, состоящий из трех взаимно ортогональных собственных единичных векторов, для которого

T = t1el ® e l + tiei ® e i + taea® e e.

(63)

Иными словами, матрица Т, представляющая тензор Т в этом базисе—диагональная, а ее диагональные элементы как раз явля­ ются собственными значениями тензора.

Напомним в общих чертах доказательство этого классического утверждения. Рассмотрим общий случай векторного евклидова пространства п измерений; дока­

зательство проведем методом индукции по размерности пространства.

В

одномер­

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

результат

очевиден,

так

как здесь любой тензор

имеет вид

t e ® e , где е — единичный вектор этого пространства. Допустим,

что формула (63)

справедлива для пространства п1 измерений, и докажем, что она также

верна

для п измерений.

 

 

 

 

 

 

 

 

 

 

 

Пусть X — вектор пространства Ел, х

единичный вектор, направленный

вдоль X , Т (Ху

Y)— билинейная

симметричная

форма,

соответствующая

рассмат­

риваемому тензору

Т. Положим

 

 

 

 

 

 

 

 

 

 

 

Т (Х %X)

T u X iX j

’ l< / = 1,2....... ”•

 

 

 

^

 

 

х

х

~~ b ifX iX j

 

 

 

Множество

единичных

векторов

х

из

Ел

компактно,

следовательно,

функция Q (х)

имеет на этом множестве по меньшей мере один

максимум,

кото­

рый обозначим t±t и пусть в \— одно из значений x t при котором этот максимум достигается.

Пусть Е? " 1 векторное евклидово пространство с размерностью л — 1, орто­

гональное е\ Н вектор из Е ?"1, ортогональный е±. Положим Х = е± -\ Н и вы­ числим значение выражения (64) на этом векторе. Получаем

T f a + M f , е + к Н )

к + 2 Х Т (в й Н ) + к 2Т (Н , И)

 

 

 

1+ к2Н Н

 

 

1+ к 2Н н

 

 

 

 

Это выражение,

если

его

рассматривать

как

функцию к, достигает макси­

мума, очевидно, при Ы О

и любом

Н из Е?” 1.

 

 

 

 

 

Следовательно, Т (ёъ

Н) =

0 при любом

Н из

Е ?"1.

 

 

 

 

В таком случае тензор

 

 

 

 

 

 

 

 

 

 

 

Т — txe i ® е±

 

 

 

 

 

в базисе Ел, образованном из

е\ и

ортонормирэванного,

но

произвольного

ба­

зиса ер = 1, 2, .

п— 1) пространства Е?” 1,

выражается

однозначно в

виде

линейной комбинации тензорных произведений двух векторов

ер и, следовательно,

его можно рассматривать

как симметричный тензор из Е?” 1

По

предположению,

теорема справедлива для

п— 1 измерения, следовательно,

можно

найти ортонор­

мированный репер et.-.en в Еп, в котором— тензор Т запишется в таком виде:

h*i ® *i+*2*а ® е2+ ---+tnen® еп.

Заметим теперь, что в системе осей £?/ тензор представлен диагональной мат­ рицей и, следовательно, величины — собственные значения тензора Т, в связи

с чем можно отождествить t{ и //. Теорема, таким образом, полностью доказана. Заметим, что, если

Н^ ^2 ^ ts

^П»

 

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

 

tp = sup Q (x)t

 

 

где x пробегает компактное множество:

 

 

х х = 1, х-е1 = х-е2= .. .x-ep_i=zQ.

t2, t3 различа­

Если в выражении (63) собственные значения

ются, то всякий собственный вектор коллинеарен

одному из векто­

ров ё{.

 

 

В самом деле, уравнения, из которых определяются компоненты

а{ некоторого

собственного вектора а,

в собственном ортонормиро-

ванном базисе

е {

запишутся в таком

виде:

 

(^1

И')^1==0» (^2

(ta р,)дз1=а0,

(65)

и если

то

а2 = а3==0.

 

 

Если два (и только два) из собственных значений равны, на­ пример, t1^=t2^ t 3y то собственному значению р = / 8 соответствуют

собственные векторы, коллинеарные е8. Собственному (дважды пов­ торяющемуся) значению t1= t2 соответствует любой вектор, ортого­

нальный ё8, так как из формулы (65) следует, что а8 = 0. В этом случае данный тензор— это тензор вращения с осью е3. Любая пара единичных векторов, ортогональных ё3, составляет вместе с векто­

ром е3 собственный

базис тензора

Т

 

При t1 = t2 = t9

имеется одно

лишь (трехкратное) собственное

значение. В этом случае тензор шаровой.

 

в) Следствия. В

данном главном базисе

ё 2У е3 тензор Т пол­

ностью определяется собственными значениями tlt t2, /8, следова­

тельно,

любой инвариант тензора

Г —функция

/1Э t2, t3. Кроме того,

как это

видно из преобразований

базиса типа

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

ров ег, е2, е8, любой инвариант также симметричная функция этих переменных. Элементарные инварианты, в частности, очень просто выражаются через собственные значения:

Тх= tx+ 12+ t3, T u =*t1t2-\-t2t3-\-t3ti, Tlu=*txt2t3.

(66)

Так как эти инварианты полностью определяют (с точностью до порядка) собственные значения tlt /2, t3, то любой инвариант — функция Ти Тп, Т ш. Таким образом, доказано, что изотропная функция матрицы Т, удовлетворяющая тождеству (57), необходимо имеет вид

/ =*=ф (7"I, Т п, Гш).

(67)