Курс механики сплошных сред
..pdfновка из (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) |