Pull to refresh

Магия тензорной алгебры: Часть 3 — Криволинейные координаты

Reading time 8 min
Views 64K

Содержание


  1. Что такое тензор и для чего он нужен?
  2. Векторные и тензорные операции. Ранги тензоров
  3. Криволинейные координаты
  4. Динамика точки в тензорном изложении
  5. Действия над тензорами и некоторые другие теоретические вопросы
  6. Кинематика свободного твердого тела. Природа угловой скорости
  7. Конечный поворот твердого тела. Свойства тензора поворота и способ его вычисления
  8. О свертках тензора Леви-Чивиты
  9. Вывод тензора угловой скорости через параметры конечного поворота. Применяем голову и Maxima
  10. Получаем вектор угловой скорости. Работаем над недочетами
  11. Ускорение точки тела при свободном движении. Угловое ускорение твердого тела
  12. Параметры Родрига-Гамильтона в кинематике твердого тела
  13. СКА Maxima в задачах преобразования тензорных выражений. Угловые скорость и ускорения в параметрах Родрига-Гамильтона
  14. Нестандартное введение в динамику твердого тела
  15. Движение несвободного твердого тела
  16. Свойства тензора инерции твердого тела
  17. Зарисовка о гайке Джанибекова
  18. Математическое моделирование эффекта Джанибекова


Введение


Читая отзывы к своим статьям, понял, что я излишне перегрузил читателя теоретическими вводными. Прошу за это прощения, признаться честно, я сам далек от формальной математики.

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

А пока что рассмотрим некоторые теоретические основы. Добро пожаловать под кат.

1. Матрица Якоби и локальная метрика. «Жонглирование» индексами


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

Простейший жизненный пример криволинейной системы координат — географические координаты (\varphi, \lambda, h) — широта, долгота и высота над уровнем моря, по которой определяется положение объектов вблизи поверхности Земли. Криволинейные координаты широко применяются в астрономии. В механике примером таких координат могут служить обобщенные координаты механической системы, однозначно определяющие её положение в пространстве с учетом геометрии наложенных на систему связей. На этом и строится аналитическая механика.



Рис. 1. Криволинейные координаты в трехмерном пространстве

Рассмотрим криволинейные координаты, заданные в трехмерном евклидовом пространстве (рисунок 1). Пусть положение точки задается в этих координатах вектором

\mathbf{q} = \begin{bmatrix} q^1 \\ q^2 \\ q^3 \end{bmatrix} \quad (1)

и декартовы координаты точки связаны с (1) соотношением

\vec{r} = \vec{r}\left(q^1, q^2, q^3\right) \quad (2)

или, в компонентной форме

x^i\left(q^1, q^2, q^3\right) \quad (3)

Рассмотрим частную производную \cfrac{\partial \vec{r} }{\partial q^1}. Результат такого дифференцирования — это вектор, направленный по касательной к координатной линии q^1. Дифференцируя (2) по всем криволинейным координатам получим тройку векторов

\vec{e}_1 = \cfrac{\partial \vec{r} }{\partial q^1}, \quad \vec{e}_2 = \cfrac{\partial \vec{r} }{\partial q^2}, \quad \vec{e}_3 = \cfrac{\partial \vec{r} }{\partial q^3} \quad (4)

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

Векторы (4) собирают в матрицу

\mathbf{J} = \begin{bmatrix} \cfrac{\partial \vec{r} }{\partial q^1} && \cfrac{\partial \vec{r} }{\partial q^2} && \cfrac{\partial \vec{r} }{\partial q^3} \end{bmatrix} = \begin{bmatrix} \cfrac{\partial x^1 }{\partial q^1} && \cfrac{\partial x^1 }{\partial q^2} && \cfrac{\partial x^1 }{\partial q^3} \\ \cfrac{\partial x^2 }{\partial q^1} && \cfrac{\partial x^2 }{\partial q^2} && \cfrac{\partial x^2 }{\partial q^3} \\ \cfrac{\partial x^3 }{\partial q^1} && \cfrac{\partial x^3 }{\partial q^2} && \cfrac{\partial x^3 }{\partial q^3} \end{bmatrix} \quad (5)

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

\mathbf{J} = \cfrac{\partial \vec{r}}{\partial \vec{q}}

Легко догадаться, что если функция (2) линейна относительно компонент вектора \vec{q}, то её можно выразить матричным соотношением

\vec{r} = \mathbf{A}_{01} \vec{q}

то мы рассматриваем косоугольную систему координат, и матрица Якоби будет равна матрице преобразования от косоугольных координат к декартовым

\mathbf{J} = \mathbf{A}_{01}

Теперь, любой вектор, заданный в пространстве (тензор ранга (1, 0)) можно представить через его контравариантные компоненты в криволинейной системе координат

\vec{a} = a^i \, \vec{e}_i = a^i \, \cfrac{\partial \vec{r} }{\partial q^i} \quad (6)

Однако, компоненты вектора, из-за переменного базиса, будут зависеть от положения в пространстве точки приложения вектора. Кроме того, для того чтобы представление (6) существовало, надо чтобы векторы, составляющие базис были не компланарны. Из курса векторной алгебры нам известно, что векторы некомпланарны, если их смешанное произведение отлично от нуля. Отсюда возникает условие, которому должен удовлетворять определитель матрицы Якоби

\det (\mathbf{J}) \ne 0 \quad (7)

Данный определитель как раз определяет смешанное произведение векторов базиса.

Теперь вычислим ковариантные компоненты вектора \vec{a}. Для этого, в самой первой статье цикла, мы умножали вектор скалярно на соответствующий вектор базиса

\vec{a} = a_i = \vec{a} \cdot \vec{e}_1 = a^i \, (\vec{e}_i \cdot \vec{e}_1) = a^i \, \cfrac{\partial \vec{r} }{\partial q^i} \cdot \cfrac{\partial \vec{r} }{\partial q^i}

В той же, первой статье, мы определили, что ковариантные компоненты вектора связаны с контравариантными через метрический тензор g_{ij}

a_j = a^i g_{ij}

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

g_{ij} = \cfrac{\partial \vec{r} }{\partial q^i} \cdot \cfrac{\partial \vec{r}}{\partial q^j}

которое можно представить в матричной форме

\mathbf{g} = \mathbf{J}^T \, \mathbf{J}

Эту связь можно представить и в тензорной форме, но для этого придется ввести явно метрику для декартовых координат \hat{g}_{kl}

\hat{g}_{kl} = \begin{bmatrix} 1 && 0 && 0 \\ 0 && 1 && 0 \\ 0 && 0 && 1 \end{bmatrix}

Тогда, преобразование декартовой метрики в криволинейную будет выглядеть так

g_{ij} = J_i^k \, J_j^p \, g_{kl} \quad (8)

Выражение (8) вводит метрический тензор для криволинейных координат. Этот тензор зависит от положения точки в пространстве, поэтому говорят что он задан локально или определяет локальную метрику

Определившись с метрикой, мы можем записать правила преобразования контравариантных координат в ковариантные

a_j = g_{ij} \, a^i \quad (9)

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

a^k = a_i g^{ik} \quad (10)

В тензорном исчислении операции опускания (9) и поднятия (10) индексов называют «жонглированием» индексами.

Выписав соотношения (9) и (10) мы подразумевали, что матрицы g_{ij} и g_{ij} взаимно обратимы. Это возможно лишь в том случае, если

\det |g_{ij}| \ne 0

Данное условие выполняется для криволинейных координат, если матрица Якоби не вырождена, и это непосредственно следует из (8), так как

\det |g_{ij}| = \det |J_i^k J_j^l g_{kl}| = \det|J_i^k| \, \det|J_j^l| = J^2 \ne 0

то есть условие (7) выполняемое для всех точек пространства — достаточное условие невырожденности локальной метрики.

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

g_{ik} \, g^{kj} = \delta_i^j

где

\delta_i^j = \begin{cases} 1, \quad i = j \\ 0, \quad i \ne j \end{cases}

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

2. Взаимный базис


Введем векторы \vec{e}^{\, 1}, \, \vec{e}^{\,2}, \, \vec{e}^{\,3}, получаемые из векторов исходного базиса путем поднятия индекса

\vec{e}^{\,j} = g^{ji} \, \vec{e}_i \quad (11)

Теперь возьмем и умножим (11) скалярно на вектор \vec{e}_k

\vec{e}^{\,j} \cdot \vec{e}_k = g^{ji} \, \vec{e}_i \cdot \vec{e}_k

но, мы знаем, что \vec{e}_i \cdot \vec{e}_k = g_{ik} — метрический тензор, поэтому, приходим к уравнению

\vec{e}^{\,j} \cdot \vec{e}_k = g^{ji} \, g_{ik}

\vec{e}^{\,j} \cdot \vec{e}_k = \delta_k^j \quad (12)

Если мы возьмем, например, вектор \vec{e}^{\,1}, то в силу (12) он перпендикулярен векторам \vec{e}_2 и \vec{e}_3 (его скалярное произведение с ними равно нулю), а скалярное произведение этого вектора на \vec{e}_1 — равно единице

Дальше возьмем и умножим (11) скалярно на \vec{e}^{\,k}

\vec{e}^{\,j} \cdot \vec{e}^{\,k} = g^{ji} \, \vec{e}_i \cdot \vec{e}^{\,k}

и в силу (12) это дает контравариантный метрический тензор

\vec{e}^{\,j} \cdot \vec{e}^{\,k} = g^{jk} \quad (13)

Система векторов \vec{e}^{\, 1}, \, \vec{e}^{\,2}, \, \vec{e}^{\,3} тоже образует базис, который называют взаимным или сопряженным с базисом \vec{e}_{\, 1}, \, \vec{e}_{\,2}, \, \vec{e}_{\,3}.

Снова рассмотрим вектор \vec{a}. Из соотношений (10) и (11) следует цепочка преобразований

\vec{a} = a^i \, \vec{e}_{\,i} = a_j \, g^{ji} \, \vec{e}_{\,i} = a_j \, \vec{e}^{\,j} \quad (13)

Умножим (13) скалярно на \vec{e}^{\,k}

\vec{a} \cdot \vec{e}^{\,k} = a_j \, \vec{e}^{\,j} \cdot \vec{e}^{\,k} = a_j \, g^{jk} = a^k

приходим к заключению, что любой вектор может быть разложен как по базису \vec{e}_{\,i} — тогда его компоненты будут контравариантные, так и по базису \vec{e}^{\,i} — компоненты будут ковариантными

\vec{a} = a^i \, \vec{e}_{\,i} = a_i \, \vec{e}^{\,i}

При этом, ковариантные компоненты — это скалярные произведение вектора на векторы базиса \vec{e}_{\,i}, а контравариантные компоненты — скалярные произведения вектора на векторы базиса \vec{e}^{\,i}

a_i = \vec{a} \cdot \vec{e}_{\,i}, \quad a^i = \vec{a} \cdot \vec{e}^{\,i} \quad (15)

что ещё раз иллюстрирует взаимность этих базисов.

Тут надо отметить, что векторы базиса \vec{e}_{\,i} получаются естественным путем — они касательны соответствующим координатным линиям и им можно приписать геометрический смысл. Что касается базиса \vec{e}^{\,i}, то его векторы не направлены по касательной координатным линиям, а перпендикулярны парам векторов касательного базиса. Такой базис иногда принято называть неголономным

3. Преобразование криволинейных координат. Формальное определение ковариантных и контравариантных компонент


Допустим, что мы работаем в криволинейной система координат, определенной вектором \vec{q}. Перейдем к другой системе координат, положение точек которой определяется вектором \vec{p}, таким, что преобразование от старой системы координат к новой определяется уравнениями

p^i = p^i(q^1,\, q^2, \, q^3) \quad (16)

Будем считать преобразование (16) обратимым, то есть допустим существование функции

q^i = q^i(p^1,\, p^2, \, p^3) \quad (17)

Для этого требуется, чтобы определитель матрицы Якоби

A_j^i = \frac{\partial p^i}{\partial q^j} \quad (18)

был отличен от нуля

\det |A_j^i| \ne 0

Тогда существует матрица B_j^i, обратная матрице (18), такая, что

A_k^j \, B_i^k = \delta_i^j

Матрица B_j^i является матрицей Якоби для преобразования (17). Тогда можно вычислить векторы нового базиса

\vec{e}_{\,i}^{\,'} = \frac{\partial \vec{r}}{\partial p^i} = \frac{\partial \vec{r}}{\partial q^j} \, \frac{\partial \vec{q^j}}{\partial p^i} = \frac{\partial \vec{r}}{\partial q^j} \, B_i^j

Получаем связь между старым базисом и новым

\vec{e}_i^{\,'} = \vec{e}_k \, B_i^k \quad (19)

\vec{e}_i = \vec{e}_k^{\,'} \, A_i^k \quad (20)

Разложим вектор \vec{a} в новом базисе

\vec{a} = a^{'i} \, \vec{e}_{\,i}^{\,'}

и используя соотношение (19), напишем

\vec{a} = a^{'i} \, B_i^k \vec{e}_{\,k} = a^k \, \vec{e}_{\,k} \quad (21)

С учетом того, что векторы базиса линейно независимы, приравниваем коэффициенты при них в (21)

a^k = a^{'i} \, B_i^k

Теперь умножим обе части (21) на A_k^j

a^k \, A_k^j = a^{'i} \, B_i^k \, A_k^j

Учтем, что

B_i^k \, A_k^j = \delta_i^j

a^k \, A_k^j = a^{'i} \, \delta_i^j = a^j

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

a^{'j} = a^k \, A_k^j \quad (22)

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

Действительно, чтобы получить векторы нового базиса, мы использовали матрицу B_i^k по формуле (19). Чтобы получить контравариантные компоненты заданного в новом базисе вектора мы используем матрицу A_k^j

А теперь посмотрим, как преобразуется вектор, заданный своими ковариантными компонентами

a_j^{\,'} = \vec{a} \cdot \vec{e}_j^{\,'} = \vec{a} \cdot \vec{e}_j \, B_i^j = a_j \, B_i^j \quad (23)

Из (23) видно, что
Ковариантные компоненты вектора преобразуются тем же оператором, которым осуществляется преобразование базиса

Формулы (19), (22) и (23) и сформулированные определения, вынесенные в цитатный блок дают формальное определение контравариантных и и ковариантных координат и иллюстрируют разницу между ними. Можно сформулировать утверждение
Тензор ранга (1,0) преобразуется оператором обратным, используемому при преобразовании базиса, а тензор ранга (0,1) преобразуется тем же самым оператором, что используется при преобразовании базиса.

4. Ковариантная производная. Символы Кристоффеля 2-го рода


Предположим, что мы хотим продифференцировать вектор, заданный произвольными координатам по какой-то из координат. Что мы должны сделать? Давайте попробуем выполнить эту операцию

\frac{\partial \vec{a}}{\partial q^j} = \frac{\partial}{\partial q^j}(a^i \, \vec{e}_i) = \frac{\partial a^i}{\partial q^j} \, \vec{e}_i + a^i \frac{\partial \vec{e}_i}{\partial q^j} \quad (24)

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

\frac{\partial \vec{e}_i}{\partial q^j} = \Gamma_{ij}^k \, \vec{e}_k \quad (25)

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

\frac{\partial g_{ij}}{\partial q^k} = \frac{\partial \vec{e}_i}{\partial q^k} \cdot \vec{e}_j + \vec{e}_i \cdot \frac{\partial \vec{e}_j}{\partial q^k} \quad (26)

Подставим (25) в (26)

\frac{\partial g_{ij}}{\partial q^k} = \Gamma_{ik}^m \, (\vec{e}_m \cdot \vec{e}_j) + \Gamma_{jk}^m \, (\vec{e}_m \cdot \vec{e}_i)

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

\frac{\partial g_{ij}}{\partial q^k} = \Gamma_{ik}^m \, g_{mj} + \Gamma_{jk}^m \, g_{mi} \quad (27)

Прежде чем начать работать с (27), скажем, что искомые коэффициенты разложения симметричны относительно нижних индексов, так как проведя прямое дифференцирования базисного вектора приходим к выражению

\frac{\partial \vec{e}_i}{\partial q^j} = \frac{\partial^2 \vec{r}}{\partial q_i \partial q_j} = \Gamma_{ij}^k \, \vec{e}_k

откуда, в силу непрерывности рассматриваемых функций, заключаем, что

\Gamma_{ij}^k = \Gamma_{ji}^k \quad (28)

Теперь, в (27) переставим индексы i и k

\frac{\partial g_{kj}}{\partial q^i} = \Gamma_{ki}^m \, g_{mj} + \Gamma_{ji}^m \, g_{mk} \quad (28)

А теперь, переставим в (27) индексы j и k

\frac{\partial g_{ik}}{\partial q^j} = \Gamma_{ij}^m \, g_{mk} + \Gamma_{kj}^m \, g_{mi} \quad (29)

Теперь сложим (29) и (30) учтя при этом симметричность (28)

\frac{\partial g_{kj}}{\partial q^i} + \frac{\partial g_{ik}}{\partial q^j} = \Gamma_{kj}^m \, g_{mj} + \Gamma_{ij}^m \, g_{mk} + \Gamma_{ji}^m \, g_{mk} + \Gamma_{kj}^m \, g_{mi}

\frac{\partial g_{kj}}{\partial q^i} + \frac{\partial g_{ik}}{\partial q^j} = \Gamma_{kj}^m \, g_{mj} + 2\,\Gamma_{ij}^m \, g_{mk} + \Gamma_{kj}^m \, g_{mi} \quad (31)

Вычитаем (27) из (31), снова учитывая (28)

\frac{\partial g_{kj}}{\partial q^i} + \frac{\partial g_{ik}}{\partial q^j} - \frac{\partial g_{ij}}{\partial q^k} = \Gamma_{ki}^m \, g_{mj} + 2 \, \Gamma_{ij}^m \, g_{mk} + \Gamma_{kj}^m \, g_{mi} - \Gamma_{ik}^m \, g_{mj} - \Gamma_{jk}^m \, g_{mi}

\frac{\partial g_{kj}}{\partial q^i} + \frac{\partial g_{ik}}{\partial q^j} - \frac{\partial g_{ij}}{\partial q^k} = 2 \, \Gamma_{ij}^m \, g_{mk} \quad (32)

Умножаем (32) на \cfrac{1}{2} \, g^{mk}, и получаем окончательно

\Gamma_{ij}^m = \frac{1}{2} \, g^{mk}\left( \frac{\partial g_{kj}}{\partial q^i} + \frac{\partial g_{ik}}{\partial q^j} - \frac{\partial g_{ij}}{\partial q^k} \right) \quad (33)

Выражение (33) определяет так называемый символ Кристоффеля 2-го рода. Тогда

\frac{\partial \vec{a}}{\partial q^j} = \frac{\partial a^m}{\partial q^j} \cdot \vec{e}_m + a^i \, \Gamma_{ij}^m \, \vec{e}_m = \left(\frac{\partial a^m}{\partial q^j} + a^i \, \Gamma_{ij}^m \right)\cdot \vec{e}_m \quad (34)

Выражение, стоящее в скобках в (34) называется ковариантной производной контравариантных компонент вектора

\nabla_j \, a^m = \left(\frac{\partial a^m}{\partial q^j} + a^i \, \Gamma_{ij}^m \right) \quad (35)

Исходя из (35) мы должны понимать, что пытаясь дифференцировать по криволинейной координате, мы обязаны учитывать зависимость базиса от координат. Если метрика не зависит от положения точки приложения вектора в пространстве, то (35) превращается в частную частную производную, ибо все символы Кристоффеля равны будут нулю, из-за того что метрический тензор не зависит от координат. В любой косоугольной системе координат, и в их частном случае — декартовых координатах, символы Кристоффеля, согласно (33) равны нулю. А значит, согласно (35) ковариантрая производная от вектора по координате будет совпадать с его частной производной по этой координате, к чему мы приучены вобщем-то давно. Но если бы (33) был тензором, то он, будучи равен нулю, остался бы нулевым в любой другой системе координат. Но в криволинейных координатах (33) нулю не равны. А значит символы Кристоффеля не являются тензором. При преобразовании системы координат меняются компоненты, но не сущность тензора. Нулевой тензор должен быть таковым в любой системе координат.

Заключение


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

Продолжение следует…
Tags:
Hubs:
+29
Comments 8
Comments Comments 8

Articles