Как найти тензор поворота

Рассмотрим два правых декартовых базиса Соотношение между ними вполне определено направляющими косинусами Но можно написать и так:

называется тензором поворота.

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

Третий инвариант (детерминант, определитель)

Важное свойство тензора поворота выражается равенством

т. е. обратный тензор. Но не только тензоры поворота обладают таким свойством. Если в (8.1) один из базисов правый, а другой левый, получим комбинацию поворота с отражением; при этом будет

Любой тензор второго ранга имеет по крайней мере одно вещественное главное значение (корень Для тензора поворота оно Соответствующая главная ось называется осью поворота; известная теорема Эйлера в том и состоит, что такая ось существует. Если — орт этой оси, величина угла поворота, то

Доказывается эта формула так. На оси поворота в перпендикулярной плоскости

Пусть теперь тензор поворота есть функция времени: . В этом случае можно ввести вектор угловой скорости следующим образом. Продифференцируем тождество (8.3) по времени:

Тензор оказался антисимметричным, и в соответствии с (6.2) его можно представить в виде Значит,

Помимо этого общего представления, вектор допускает и другие. Подставив (8.4) в (8.5), получим

Здесь вектор со разложен по трем взаимно перпендикулярным направлениям При неподвижной оси поворота приходим к простейшему

Еще одно представление со связано с его компонентами (8.2). Поскольку и

Отметим и формулы

Далее повсеместно будет использоваться сходная с дифференцированием операция варьирования. Не отсылая читателя к курсам вариационного исчисления, ограничимся представлением о вариации величины х как о задаваемом нами бесконечно малом приращении, совместимом с ограничениями — связями. Если ограничений для х нет, то произвольна. Но когда функция независимого аргумента у, следует считать

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

Варьируя тождество (8.3), получим Этот тензор антисимметричен, и потому выражается через свой сопутствующий вектор Приходим к соотношениям

полностью аналогичным (8.5). Вектор малого поворота нельзя считать “вариацией о”, это единый символ (в отличие от

Малый поворот определяется вектором но конечный поворот также допускает векторное представление: согласно (8.4), вся информация о тензоре заключена в векторе Когда за одним поворотом следует другой, их векторы не складываются в результирующем повороте. Однако это не лишает их векторного характера.

Тензор поворота и кватернион

Вот мы и добрались до самого интересного
места этой статьи. А именно, до связи
тензора поворота с еще одним замечательным
объектом – кватернионом поворота. Для
начала перепишем формулу (54) для
повернутого вектора a′ используя
функции половинного аргумента и
вернувшись обратно к записи через вектор
m

a′ = a + sinθm × a + (1 – cosθ)m × (m × a) = = a + 2sin(θ/2)cos(θ/2)m × a + 2sin2(θ/2) m × (m × a).

(59)

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

λ0 = cos(θ/2), λ = sin(θ/2)m.

(60)

Полезное свойство данных обозначений

(λ0)2 + |λ|2 = (λ0)2 + λ · λ = 1.

(61)

Выражение (59) можно переписать в виде

a′ = a + 2λ0λ × a + 2λ × (λ × a) = = (λ0)2a + λ(λ · a) + 2λ0λ × a + λ × (λ × a).

(62)

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

Q(θ) = Q(φ)·Q(ψ).

(63)

Для этого воспользуемся формулой (52),
записанной также через функции половинного
аргумента и введенные обозначения (60)

P = 2sin2(θ/2)mm + 2[cos2(θ/2) – 1/2]1 + 2sin(θ/2)cos(θ/2)m × 1 = = 2{λλ + [(λ0)2 – 1/2]1 + λ0λ × 1}

(64)

Сопоставим векторам φ и ψ
пары

μ0 = cos(|φ|/2), μ = sin(|φ|/2)φ /|φ |,	(65)
ν0 = cos(|ψ|/2), ν = sin(|ψ|/2)ψ /|ψ|.	(66)

Подставив данные выражения в (63), получим

Q(θ) = Q(φ)·Q(ψ) = = 4{μμ + [(μ0)2 – 1/2]1 + μ0μ × 1}·{νν + [(ν0)2 – 1/2]1 + ν0ν × 1}= = 4{μν(μ · ν) + μμ[(ν0)2 – 1/2] + ν0μμ × ν + + νν[(μ0)2 – 1/2] + [(μ0)2 – 1/2] [(ν0)2 – 1/2]1 + [(μ0)2 – 1/2] ν0ν × 1 + + μ0μ × νν + [(ν0)2 – 1/2]μ0μ × 1 + μ0ν0μ × 1 × ν} = = 2{λλ + [(λ0)2 – 1/2]1 + λ0λ × 1}

(67)

Напомним, что для любых векторов b
и c справедливы равенства

b × 1 = 1 × b = – (b × 1)T, tr(b × 1) = 0, tr(b × cc) = (b × c) · c = 0,

tr(cc) = c · c, b × 1 × c = cb – (c · b)1, tr(b × 1 × c) = – 2(c · b),

(b × 1)× = –2b, (b × cc)× = (b × c) × c = c(b · c) – b(c · c), (cc)× = 0, (b × 1 × c)× = c × b.

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

λ0 = μ0ν0 – μ · ν,	(68)
λ = μ0ν + ν0μ + μ · ν.	(69)

Далее можно ввести множество объектов,
представляющих собой пару Λ = (λ₀, λ),
на котором задан закон композиции
элементов (операция умножения), в
соответствии с формулами (68, 69).

Можно убедиться, что данный закон
умножения является ассоциативным
(Λ·Μ)·Ν = Λ·(Μ·Ν), дистрибутивным
(Λ + Μ)·(Ν + Ρ) = Λ·Ν + Λ·Ρ + Μ·Ν + Μ·Ρ
и удовлетворяет условию (αΛ)·(βΜ) = αβΛ·Μ,
где α, β– произвольные вещественные
числа. Операция суммы вводится следующим
образом Λ + Μ = (λ₀, λ) + (μ₀, μ) = (λ₀ + μ₀, λ + μ).
В математике данное множество называется
алгеброй кватернионов, а сами
объекты кватернионами. Выше
нами были рассмотрены кватернионы
единичной длины (λ₀)² + λ · λ = 1,
которые наравне с тензором поворота
могут использоваться для задания
ориентации тела в трехмерном пространстве.
В общем случае кватернионы могут иметь
произвольную длину. Более подробно о
кватернионах и о том, как их использовать
для описания поворотов, будет написано
Zed’ом в следующем номере журнала.

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

• #
• #
• #
• #
• #
• #
• #
• #
• #
• #
• #

det(P – 1) = det[P·(1 – P T )] =
= det(P)det(1 – P T ) = det(1 – P T ) =
= det(1 – P) = – det(P – 1).

Получили, что det(P – 1) = –det(P – 1), а это возможно, если верно равенство (37). Итак, ненулевое решение уравнения (36) существует. Понятно, что если m есть решение (36), то и (–m) есть решение этого уравнения. Однако оба эти вектора определяют одну и ту же прямую, натянутую на вектор m.

Далее методом «от противного» можно показать, что для тензора поворота P, отличного от 1, существует только один неподвижный вектор ± m.

Выберем теперь правую ортонормированную тройку векторов d₁ = m, d₂, d₃. Согласно (32) мы можем записать

P = D1d1 + D2d2 + D3d3 = D1m + D2d2 + D3d3,

(39)

где тройка D₁, D₂, D₃ также является правой ортонормированной. Умножив равенство (39) скалярно на m справа и учтя (35), получим

Таким образом, равенство (39) принимает вид

P = mm + D2d2 + D3d3.

(41)

D2 = cosθd2 + sinθd3, D3 = – sinθd2 + cosθd3.

(42)

Подставляя эти разложения в (41) и проводя элементарные преобразования, получаем

P = mm + (cosθd2 + sinθd3)d2 + (–sinθd2 + cosθd3)d3 = = mm + cosθ(d2d2 + d3d3) + sinθ(d3d2 – d2d3) = = mm + cosθ(d1d1 + d2d2 + d3d3 – mm) + sinθm × (d1d1 + d2d2 + d3d3) = = mm + cosθ(1 – mm) + sinθm × 1

(43)

Последним равенством и заканчивается доказательство теоремы Эйлера. Обратим внимание, что замена m на – m влечет за собой замену θ на – θ. При этом сам тензор P не меняется.

Пусть вектор a параллелен оси поворота (лежит на оси поворота), т.е. a = αm, то P · a = a и вектор a не меняется при действии на него тензора P. Пусть a ортогонален оси поворота: a · m = 0. Тогда

a′ = P · a = cosθa + sinθm × a.

(44)

Этот результат изображен на рис. 3а.

Видно, что действие тензора поворота на вектор a, ортогональный оси поворота, сводится к повороту на угол θ вокруг m. Если θ > 0, то поворот производится против часовой стрелки при взгляде с конца m. Поворот произвольного вектора a показан на рис. 3б, при этом проекция вектора a на m сохраняется, а часть вектора a, ортогональная a, поворачивается на угол θ вокруг m.

Композиция поворотов. Правило квазикоммутативности

a ′′ ≡ P2 · a′ = P2·P1 · a = P · a , P = P2·P1.

(45)

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

Это означает, что если произвести повороты в другом порядке, сначала P₂ потом P₁, то результат будет другой.

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

Q(φm)·Q(θn) = Q(φm)·Q(θn)·1 = = Q(φm)·Q(θn)·Q T (φm)·Q(φm).

(47)

Q(φm)·Q(θn)·Q T (φm) = = n′n′ + cosθ(1 – n′n′) + sinθn′ × 1,

(48)

Q(φm)·Q(θn) = Q(θn′)·Q(φm).

(49)

То есть результирующий поворот можно осуществить двумя способами. Либо сначала повернуть на угол θ вокруг n, а затем на угол φ вокруг m. Либо сначала на угол φ вокруг m, а затем на угол θ, но уже вокруг другой повернутой оси n′.

Q(φm)·Q(θn)·Q T (φm) = Q(θn′), n′ = Q(φm) · n

(50)

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

Всем известен пример из детства с юлой. В механике этот объект называется волчком Лагранжа (см. рис. 4) и представляет собой осесимметричное тело вращения с одной точкой опоры. Юла одновременно вращается вокруг собственной оси n′, которая в свою очередь вращается вокруг неподвижной вертикали m, проходящей через точку опоры.

На первый взгляд кажется, что наиболее удобной формой задания тензора поворота юлы является P = Q(θn′)·Q(φm). Но ведь вектор n′ непостоянный, а это в ряде случаев может существенно усложнить моделирование. Но так как n′ = Q(φm) · n, то, пользуясь правилом квазикоммутативности, мы можем написать P = Q(φm)·Q(θn), где вектор n является константой и определяет направление оси юлы в начальный момент. Второй вариант уже существенно менее затратен для моделирования. Тем не менее, при постановке задач первый способ оказывается более выигрышным.

Вектор поворота

Тензор поворота, записанный в форме P = mm + cosθ(1 – mm) + sinθm × 1, не имеет каких-либо ограничений и может применяться всегда. Тем не менее, интуитивным образом поворота является не тензор второго ранга, а некий вектор θ, называемый вектором поворота

Его направление совпадает с осью поворота m, а его модуль равен модулю угла поворота θ. В правоориентированной системе отсчета θ > 0 при взгляде с конца m, если поворот происходит против часовой стрелки. Поскольку поворот выражает реальное явление, то наличие нескольких способов его описания возможно только при условии, что способы связаны взаимно однозначным образом. Проще всего показать, что тензор поворота однозначно может быть вычислен по вектору поворота, что следует непосредственно из представления тензора поворота доставляемого теоремой Эйлера

P = [(1– cosθ )/θ 2 ]θθ + cosθ1 + [sinθ/θ]θ × 1.

(52)

Видно, что последнее выражение нечувствительно к знаку θ. Поэтому величину θ в (52) можно рассматривать как модуль вектора θ, что позволяет без дополнительных соглашений однозначно вычислять тензор поворота по вектору поворота. В теории поворотов доказывается теорема, что представление (52) является ничем иным как тензорной экспонентой тензора θ × 1

Вам это ничем не напоминает комплексную экспоненту exp(θi), которая является оператором поворота векторов на комплексной плоскости вокруг оси перпендикулярной к ней? Общность явлений порождает общность описаний, правда со своими особенностями.

Не будем приводить здесь формулы для вычисления вектора поворота по тензору поворота, скажем лишь, что это можно также однозначно сделать. В заключение этого пункта приведем формулу для повернутого вектора a′

a′ ≡ P(θ) · a = a + [sinθ/θ]θ × a + [(1– cosθ )/θ 2 ]θ × (θ × a),

(54)

Теорема о представлении тензора поворота

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

Теорема. Любой поворот P(θ) может быть осуществлен в виде композиции поворотов вокруг произвольно выбираемых и фиксированных во времени осей m и n, а также оси e, выбираемой по правилу e = m × n / |m × n|,

P(θ) = Q(ψm)·Q( J e)Q(φn).

(55)

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

Если m = n, то углы ψ, J , φ называются углами Эйлера (углы прецессии, нутации и собственного вращения). Вектор e в этом случае произвольный, удовлетворяющий условию e · n = 0. Тензор поворота имеет вид

P(θ) = Q(ψm)·Q( J e)·Q(φm).

(56)

Однако некоторые читатели могут возразить: «Позвольте, но в учебниках по механике написано, что при использовании углов Эйлера оси как-то поворачиваются, а здесь они неподвижны…». А мы в ответ, в учебниках все правильно написано, только во всех трактовка однобокая. Вспомним правило квазикоммутативности поворотов. По нему мы можем записать

P(θ) = Q(φm′′)·Q( J e′)·Q(ψm),

(57)

Видим, что в классическом представлении повороты производятся в обратном порядке и вокруг вращающихся векторов e′ и m′′. Последнее обстоятельство наиболее часто вызывает трудности в понимании и последующем восприятии углов Эйлера у многих студентов.

На практике помимо углов Эйлера используется масса других систем углов. Наиболее известны «самолетные» и «корабельные» углы. О них мы расскажем в следующих статьях. А сейчас рассмотрим небольшой пример.

Многие, наверное, смотрели фильм «Косильщик лужаек». Там есть одна сцена, где главного героя раскручивают в непонятной машине, представляющей собой три скрепленных кольца, во внутреннее кольцо был помещен наш герой. В механике подобный механизм называется трехстепенным гироскопом и схематично изображен на рис. 5. Три кольца – внутреннее зеленое, среднее красное, наружное синее – соединены так, что внутреннее кольцо может свободно вращаться относительно среднего вокруг оси m′′, среднее кольцо вращается относительно внешнего вокруг оси e′, внешнее же кольцо вращается вокруг неподвижной оси n.

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

P = Q(φm′′)·Q( J e′)·Q(ψn), e′ = Q(ψn) · e, m′′ = Q( J e′)·Q(ψn) · m = Q(ψn)·Q( J e) · m,

(58)

где e, m – орты направления осей вращения в отсчетном положении (в начальный момент времени).

Обратите внимание, что мы не стали долго думать, как же подстроить эйлерову систему углов для нашей задачи, а, просто глядя на задачу, сходу написали ответ. Именно последнее обстоятельство является тем огромным преимуществом прямого тензорного языка при описании поворотов, которое заставляет приложить минимум усилий для его изучения с целью максимальной практической отдачи. Рекомендуем читателю получить остальные 5 вариантов представления тензора P.

Тензор поворота и кватернион

Вот мы и добрались до самого интересного места этой статьи. А именно, до связи тензора поворота с еще одним замечательным объектом – кватернионом поворота. Для начала перепишем формулу (54) для повернутого вектора a′ используя функции половинного аргумента и вернувшись обратно к записи через вектор m

a′ = a + sinθm × a + (1 – cosθ)m × (m × a) = = a + 2sin(θ/2)cos(θ/2)m × a + 2sin 2 (θ/2) m × (m × a).

(59)

λ0 = cos(θ/2), λ = sin(θ/2)m.

(60)

Полезное свойство данных обозначений

(λ0) 2 + |λ| 2 = (λ0) 2 + λ · λ = 1.

(61)

Выражение (59) можно переписать в виде

a′ = a + 2λ0λ × a + 2λ × (λ × a) = = (λ0) 2 a + λ(λ · a) + 2λ0λ × a + λ × (λ × a).

(62)

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

Для этого воспользуемся формулой (52), записанной также через функции половинного аргумента и введенные обозначения (60)

P = 2sin 2 (θ/2)mm + 2[cos 2 (θ/2) – 1/2]1 + 2sin(θ/2)cos(θ/2)m × 1 = = 2<λλ + [(λ0) 2 – 1/2]1 + λ0λ × 1>

(64)

Сопоставим векторам φ и ψ пары

μ0 = cos(|φ|/2), μ = sin(|φ|/2)φ /|φ |,	(65)
ν0 = cos(|ψ|/2), ν = sin(|ψ|/2)ψ /|ψ|.	(66)

Подставив данные выражения в (63), получим

Q(θ) = Q(φ)·Q(ψ) = = 2<λλ + [(λ0) 2 – 1/2]1 + λ0λ × 1>

(67)

Напомним, что для любых векторов b и c справедливы равенства

b × 1 = 1 × b = – (b × 1) T , tr(b × 1) = 0, tr(b × cc) = (b × c) · c = 0,

tr(cc) = c · c, b × 1 × c = cb – (c · b)1, tr(b × 1 × c) = – 2(c · b),

(b × 1)_× = –2b, (b × cc)_× = (b × c) × c = c(b · c) – b(c · c), (cc)_× = 0, (b × 1 × c)_× = c × b.

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

λ0 = μ0ν0 – μ · ν,	(68)
λ = μ0ν + ν0μ + μ · ν.	(69)

Далее можно ввести множество объектов, представляющих собой пару Λ = (λ₀, λ), на котором задан закон композиции элементов (операция умножения), в соответствии с формулами (68, 69).

Можно убедиться, что данный закон умножения является ассоциативным (Λ·Μ)·Ν = Λ·(Μ·Ν), дистрибутивным (Λ + Μ)·(Ν + Ρ) = Λ·Ν + Λ·Ρ + Μ·Ν + Μ·Ρ и удовлетворяет условию (αΛ)·(βΜ) = αβΛ·Μ, где α, β– произвольные вещественные числа. Операция суммы вводится следующим образом Λ + Μ = (λ₀, λ) + (μ₀, μ) = (λ₀ + μ₀, λ + μ). В математике данное множество называется алгеброй кватернионов, а сами объекты кватернионами. Выше нами были рассмотрены кватернионы единичной длины (λ₀) 2 + λ · λ = 1, которые наравне с тензором поворота могут использоваться для задания ориентации тела в трехмерном пространстве. В общем случае кватернионы могут иметь произвольную длину. Более подробно о кватернионах и о том, как их использовать для описания поворотов, будет написано Zed’ом в следующем номере журнала.

Вместо заключения

Выше был рассмотрен достаточно простой с формальной точки зрения и в то же время интуитивно понятный математический объект как тензор конечного поворота. Мы надеемся, что наша статья приоткрыла для читателя завесу таинственности и непонятности в вопросе описания произвольных поворотов в трехмерном пространстве. Как было сказано во введении, мы намеренно уклонились от вычислительных аспектов данной проблемы, желая разбудить у читателя интуицию. В следующей статье мы расскажем о координатной форме представления тензора поворота – матрице поворота. Что позволит читателю реализовать в программном коде примеры из данной статьи и множество собственных еще не родившихся идей. При этом он будет делать это, мы надеемся, с глубоким пониманием происходящего.

Тем не менее, читатель уже сейчас имеет такое мощное вычислительное средство как кватернион поворота. Этот объект интуитивно также прост для восприятия, как и тензор поворота. На сегодняшний день он является реальным конкурентом классическим матрицам поворота. Поэтому для всестороннего понимания поворотов читатель должен в итоге уметь свободно обращаться с двумя объектами: тензором поворота и кватернионом поворота.

ВЕКТОРЫ И ТЕНЗОРЫ. ИХ ПРЕОБРАЗОВАНИЯ ПРИ ПОВОРОТАХ СИСТЕМЫ КООРДИНАТ

Пусть и две декартовы системы координат, повернутые друг относительно друга, с базисными векторами (ортами) , , образующими правые ортонормированные тройки. Поскольку системы координат и декартовы, то и .

Здесь – символ Кронекера.

Произвольный вектор можно разложить подобно радиус-вектору по ортам обеих систем координат:

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

Величины и называются компонентами вектора и являются ортогональными проекциями данного вектора на орты и :

Установим связь между проекциями вектора на различные базисные орты:

где – матричные элементы матрицы поворота . Если объединить компоненты в одностолбцовую матрицу , а компоненты в одностолбцовую матрицу , то закон преобразования компонент вектора можно записать в матричных обозначениях:

Задание. Убедиться в справедливости последнего равенства, раскрыв в явном виде произведение матриц.

Докажем, что матрица ортогональна, т.е. :

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

Задание. Докажите, что

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

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

Такое определение вектора допускает обобщение на случай величин с числом компонент, большим трех. Так возможны девятикомпонентные величины, компоненты которых нумеруются двумя векторными индексами , каждый из которых пробегает независимо значения 1,2,3. Возможны 27-и компонентные величины, компоненты которых нумеруются тремя векторными индексами. Наконец, возможны – компонентные величины, компоненты которых нумеруются N векторными индексами (векторные индексы независимо пробегают множество значений 1,2,3). Если компоненты этих многокомпонентных величин преобразуются по законам:

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

Вопрос. Чему равно число независимых компонент тензора ранга N в случае двумерного пространства?

Компоненты тензора второго ранга естественно объединяются в квадратную матрицу со следующим законом преобразования матричных элементов:

где и квадратные матрицы с матричными элементами и .

Очевидно, что вектор является тензором первого ранга, а скаляр – нулевого. Ранг тензора также называют тензорной размерностью, или валентностью.

2.1 Найти матрицу преобразования системы декартовых координат на плоскости при повороте на угол .

Решение задачи 2.1 Матричные элементы искомой матрицы вычисляются как скалярные произведения , здесь индексы i,j принимают только два значения: 1 или 2. Так как все орты по определению имеют единичные модули, каждое скалярное произведение равно косинусу угла между соответствующими ортами. Нарисуйте на листе бумаги пояснительный чертеж и убедитесь, что углы между парами базисных орт и одинаковы и равны углу поворота . Поэтому . Угол между ортами равен , и соответственно . Угол между ортами равен , поэтому .

2.1.1. Убедиться, что определитель матрицы равен 1.

2.1.2 Убедиться, что матрица ортогональна, т.е. , где -транспонированная матрица, а -единичная матрица.

2.1.3 Убедиться, что – матрица поворота на угол совпадает с произведением матриц и , которые являются матрицами поворота на углы и соответственно.

2.1.4. Убедиться, что матрица поворота на угол совпадает с матрицей , где – матрица поворота на угол .

• 2.2 Найти матрицу поворота в трехмерном пространстве относительно заданной координатной оси на угол .
• 2.2.1 Вокруг оси Oz

Решение задачи 2.2.1 Очевидно, что базисные орты , повернутой вокруг оси Oz системы координат, лежат в Oxy плоскости исходной координатной системы. Выше (см. задачу 2.1) мы уже вычислили скалярные произведения для i,j=1 и 2. Фактически мы нашли соответствующие им матричные элементы искомой матрицы повороты в трехмерном пространстве: , . Для нахождения остальных матричных элементов заметим, что базисные орты ортогональны орту , поэтому . После выполнения поворота вокруг оси Oz направление аналогичной оси новой системы координат не изменится, т.е. орт . Оставшиеся матричные элементы вычисляются тривиально: (j=1,2,3). Выпишем явный вид матрицы поворота вокруг оси Oz:

2.2.2 Вокруг оси Ox

Решение задачи 2.2.2 во многом аналогично решению предыдущей задачи.

Приведем в качестве ответа явный вид искомой матрицы поворота:

2.2.3 Вокруг оси Oy

2.3 Найти матрицу поворота в трехмерном пространстве на углы Эйлера. Углы Эйлера определены следующим образом: вначале проводится поворот на угол вокруг оси , затем производится поворот на угол вокруг новой оси , а после этого производится поворот на угол вокруг новой оси .

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

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

Матрица поворота составлена из косинусов углов между ортами новой исходной, и новой повернутой координатных систем. Для ее вычисления мы фактически должны повторить решение задачи 2.2.2 и получить в результате ту же матрицу с заменой угла на . Давайте примем систему координат

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

Матрица составлена их косинусов углов между соответствующими ортами. Она совпадает с матрицей поворота вокруг оси Oz, найденной в ходе решения задачи 2.2.1, с заменой угла на . Итого: . Приведем для справки явный вид матрицы поворота на углы Эйлера , , .

• 2.3.2 Доказать, что .
• 2.3.3 Выразить матрицу обратного преобразования через произведение матриц поворотов вокруг осей Ox и Oz.
• 2.4 Найти матрицу для следующих углов Эйлера:
• 2.4.1

• 2.5 В случае двумерного пространства вычислить компоненты вектора в системе координат повернутой на угол по сравнению с исходной. Компоненты вектора и угол следующие:
• 2.5.1

• 2.6 В случае двумерного пространства вычислить компоненты тензора второго ранга в системе координат, повернутой на угол по сравнению с исходной. Компоненты тензора и угол следующие:
• 2.6.1

• 2.7 В трехмерном пространстве заданы компоненты вектора. Найти компоненты вектора в системе координат, повернутой на угол вокруг оси Ox по сравнению с исходной. Компоненты вектора и угол следующие:
• 2.7.1

Решение задачи 2.7 дается общей формулой:

Для конкретного варианта, указанного в пункте 2.7.1 получаем

• 2.8 В трехмерном пространстве заданы компоненты вектора. Найти компоненты вектора в системе координат, повернутой на угол вокруг оси Oy по сравнению с исходной. Компоненты вектора и угол следующие:
• 2.8.1

• 2.9 В трехмерном пространстве заданы компоненты вектора. Найти компоненты вектора в системе координат, повернутой на угол вокруг оси Oz по сравнению с исходной. Компоненты вектора и угол следующие:
• 2.9.1

2.10 В случае двумерного пространства найти компоненты тензора в системе координат, повернутой относительно исходной на угол .

Изменение компонент тензоров при повороте осей координат

Рассмотрим общие свойства тензоров второго ранга. По определению тензор второго ранга связывает два физических вектора, т. е. векторы, имеющие длину и направление. Пусть вектор b связан с вектором а с помощью тензора [Г]. Как изменятся компоненты векторов а и b при повороте осей координат? На рис. 20.9 показан вектор а в системе координат луг. Пусть

Повернем оси координат в положение xyz . Векторы а и b при этом не изменятся, но изменятся их компоненты. Связь между ними станет

Рис. 20.9. Компоненты вектора а в системах координатxyz и xyz

Проекции вектора а в исходной системе координат в новой системе

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

Здесь, например, 1_у– —косинус угла между исходной осью у и новой осью z •

Проекции вектора а на новые оси координат xyz получим как сумму проекций старых компонент вектора а на новые оси:

То же выражение в матричной форме

Аналогично, преобразовав вектор b , получаем b =[L b . Произведем обратное преобразование последних выражений:

Второе из выражений (20.21) подставим в (20.20): [L T ]

1 а*. Левую и правую части последнего выражения умножаем слева на [ll |:

В этом выражении — единичная матрица — эквивалент умножения матрицы на единицу.

Матрица направляющих косинусов [L] — ортогональная. Для нее

IL T ] = [L]-‘ . Тогда ll 1 ‘=[/,]. В итоге выражение (20.22) приобретет вид

Сравнивая последнее выражение с выражением (20.20), видим, что

— формула подобных преобразований, (20.23)

где [Т] и [Т ] — любой тензор второго ранга в исходной и повернутой системах координат; L — матрица направляющих косинусов (углов между исходными и повернутыми осями координат).

[spoiler title=”источники:”]

http://studwood.ru/1667881/matematika_himiya_fizika/vektory_tenzory_preobrazovaniya_povorotah_sistemy_koordinat

http://studme.org/275963/tehnika/izmenenie_komponent_tenzorov_povorote_osey_koordinat

[/spoiler]