Reference shape of bodies with enhanced kinematics. Part II. Second gradient and microstructure
- Authors: Koifman K.G.1
-
Affiliations:
- Bauman Moscow State Technical University
- Issue: Vol 29, No 4 (2023)
- Pages: 54-76
- Section: Mechanics
- URL: https://journals.rcsi.science/2541-7525/article/view/310428
- DOI: https://doi.org/10.18287/2541-7525-2023-29-4-54-76
- ID: 310428
Cite item
Full Text
Abstract
The work develops differential-geometric methods for modeling finite incompatible deformations of hyperelastic solids with enhanced kinematics. The response of such bodies, along with the standard kinematic field represented by the deformation gradient, is characterized by additional tensor fields. As such, the paper considers: 1) the second deformation gradient and 2) the tensor field of the second rank, modeling the microstructure of the body. For each of these two cases, compatibility conditions are obtained and their geometric interpretation is proposed. Geometry is synthesized on the material manifold representing a body with enhanced kinematics. The corresponding affine connection has non-zero torsion and curvature, which can be useful for modeling a body with dislocations and disclinations.
Full Text
1. Предварительные сведения
. Настоящая статья продолжает работу [1], в которой представлен геометрический метод моделирования несовместных деформаций для гиперупругих тел и рассмотрены его особенности на примере простого материала. Несмотря на то что неевклидова отсчетная форма для тела из простого материала является классической и различные способы ее построения рассмотрены в статьях [2–10], подход к синтезированию неевклидовой формы, предложенный в исследовании [1], несколько отличается от этих способов. Он является комбинацией рассуждений Кренера [2], в рамках которых геометрия определяется на основе условий совместности, с идеей локальной разгрузки, предложенной в работах[1] [9; 10] для формализации локальных деформаций. Цель настоящей статьи — синтезирование неевклидовой формы для тел с расширенной кинематикой.
В работе используются основные структуры и общие формулы, определенные в [1]. В частности, есть евклидово физическое пространство [1, формула (2.1)] с ассоциированным векторным пространством . Ортонормированный базис в последнем обозначается через . Символ обозначает промежуточную форму [1, формула (2.8)], — ее подлежащее многообразие[2] и т. д. В случае необходимости приводятся ссылки на соответствующие формулы из первой части работы.
. Метод синтеза неевклидовой формы, используемый в работе, аналогичен методу подвижного репера, предложенному Картаном для построения неевклидовых пространств [11]. Остановимся более подробно на этой аналогии.
В случае простого материала переход к натуральному состоянию определяется полем локальных деформаций [1, формулы (3.8) и~(3.11)], условие совместности которых имеет вид . Это означает, что в случае выполнения последнего условия (и односвязности формы ) поле является градиентом некоторой глобальной деформации из промежуточной формы в глобально натуральную форму. Если же локальные деформации несовместны, то равенство нулю ротора не выполняется. Скажем тогда, что имеется источник несовместности, представленный тензорным полем второго ранга . В этом случае уравнение, характеризующее структурную неоднородность тела, имеет вид
(1.1)
Подход, предложенный Кренером, заключается в преобразовании условия совместности с тем, чтобы результат этого преобразования можно было интерпретировать геометрически. Действительно, в соответствии с [1, и ] условие, выражающее равенство нулю ротора , эквивалентно равенству нулю кручения специфической связности на многообразии , коэффициенты которой определяются формулами [1, (3.18)]. Следовательно, совместность деформаций эквивалентна утверждению, что связность евклидова. Но тогда появление источника несовместности равносильно изменению геометрии на многообразии : из евклидовой она переходит в неевклидову, представленную тензором кручения . В таком случае уравнение (1.1) преобразуется в равенство , характеризующее геометрическую структуру неевклидова пространства над .
С другой стороны, Картаном были получены структурные уравнения, характеризующие в общих чертах геометрию неевклидовых пространств [11–13]. Новизна идеи Картана заключалась в отказе от использования криволинейных координат и переходе к полям базисов более общего вида. Действительно, криволинейные координаты накладывают жесткие ограничения на виды локальных базисов за счет замены переменных (три функции определяют все поле базисов в трехмерном пространстве). Если же, отказываясь от использования замены переменных, перейти к неголономным базисам, то появляются дополнительные функциональные степени свободы, распоряжаясь которыми можно прийти к разнообразным геометриям. В явном виде поле базисов — подвижный репер — может быть определено по заданному полю обратимых линейных преобразований[3] трансляционного пространства в соответствии с равенствами , . Если через обозначить репер пространства , сопряженный к , а через — метрические коэффициенты относительно подвижного репера, то структурные уравнения Картана для евклидова пространства имеют вид
(1.2)
где — операция внешнего произведения [14]. Поля 1-форм , , относительно которых записаны соотношения (2), определяют взаимные искажения элементов репера при переходе от точки к точке.
Уравнения (1.2) можно рассматривать как «условия совместности» для заданных искажений подвижного репера. Действительно, в евклидовом пространстве реализуется не произвольный набор полей , а лишь тот, который удовлетворяет равенствам (1.2). Переход к пространствам более общего вида можно осуществить, определив семейства 1-форм неметричности , а также 2-форм кручения и кривизны , и подставив их в правые части уравнений (1.2):
(1.3)
Уравнения (1.3) определяют структуру пространства произвольной аффинной связности. В этой связи можно ассоциировать поля , и с величинами, выражающими несовместность полей с евклидовой геометрией.
Рассмотренная аналогия между источниками несовместности и тензорными полями кручения, кривизны и неметричности демонстрирует единство подходов, используемых в настоящей работе и в общей теории пространств аффинной связности.
2. Синтезирование неевклидовой отсчетной формы для среды второго градиента
2.1 Второй градиент деформации
. Квадратичное приближение деформации. Синтезированию поля локальных деформаций предпошлем рассмотрение кинематики среды второго градиента. Предположим, что упругий потенциал относительно промежуточной формы является отображением
(2.1)
Поэтому, если — произвольная деформация, то отклик тела в точке характеризуется равенством [15; 16]
В списке аргументов, наряду с первым градиентом деформации , представлен второй градиент , определяемый равенством . Следовательно, в окрестности точки выполнено разложение по формуле Тейлора второго порядка [17]:
(2.2)
в котором — достаточно малый вектор, т. е. . Здесь и в дальнейшем подразумевается отождествление второго градиента с билинейным отображением в силу естественного изоморфизма[4] [17] .
Обобщенная теорема Шварца [17] влечет, что второй градиент симметричен, т. е. . Поэтому его можно восстановить по значениям , в соответствии с формулой, известной из теории билинейных отображений [18]:
В свою очередь, из (2.2) вытекает следующее равенство:
доказывающее единственность второго градиента как симметричного билинейного отображения, для которого выполнено соотношение (2).
. Представления второго градиента. В прямоугольных координатах деформации соответствует представление [1, формула (2.10)], а первому градиенту деформации — разложение [1, формула (3.3)]. Поскольку второй градиент деформации можно записать как
то его разложение имеет вид[5]
Вместе с тем несмотря на то что в криволинейных координатах и первый градиент имеет простое разложение [1, формула (3.4)], представление, соответствующее второму градиенту, является более сложным. Для его получения рассмотрим следующее утверждение. Пусть
— тензорное поле, представленное в паре координат и . Тогда для его градиента
в точке справедливо разложение[6]
(2.3)
Здесь — компоненты градиента деформации [1, формула (3.4)], а и — символы Кристоффеля, отвечающие криволинейным координатам и , т. е.
Доказательство. Используя обобщенное правило дифференцирования произведения [17], получаем равенство[7]
(2.4)
В нем, согласно определению символов Кристоффеля,
.
Подставляя полученные выражения в (2.4) и заменяя соответствующим образом индексы суммирования, приходим к формуле (2.3).
В случае, когда , формула (2.3) приводит к равенству
(2.5)
являющемуся искомым разложением второго градиента. В нем
Следовательно, в силу перестановочности повторных производных, .
В работе используется следующий частный вид формулы (2.5). Выберем в качестве прямоугольные координаты . Тогда , что дает
(2.6)
В этом случае компоненты второго градиента выражаются через компоненты первого градиента и их производные, а также через символы Кристоффеля системы координат .
Замечание 1. Евклидова структура физического пространства и отождествления по изоморфизму, индуцируемые ею, позволяют скрыть истинную природу полей, используемых в настоящей работе. Действительно, предположим заданными криволинейные координаты на форме , а координаты на образах выберем прямоугольными. Тогда
- Деформация есть тройка скалярных полей , определенных равенством
- Градиент деформации есть тройка 1-форм , заданных как
То, что градиент деформации не является тензором в классическом понимании этого термина, отмечалось в работе [19, с. 245] и позднее в монографиях [20, 21].
- Второй градиент есть совокупность трех тензоров второго ранга .
В отличие от представления полей , , в исходной (евклидовой форме), их представление в соответствии с п. 1)–3) является общим и может быть непосредственно перенесено на произвольные гладкие многообразия.
Несмотря на то что в рамках настоящей работы такое описание полей избыточно[8], к нему придется прибегнуть уже в случае рассмотрения деформирования материальных двумерных поверхностей в евклидовом пространстве.
2.2 Семейство форм и условие совместности
. Гипотеза локальной разгрузки. Подобно случаю среды первого градиента [1, ], предположим справедливой гипотезу локальной разгрузки. Пусть фиксированы тензоры и , характеризующие натуральное состояние, и пусть задана промежуточная форма вместе с упругим потенциалом (2.1). Предположим далее, что определено семейство деформаций , для которого в любой точке выполнено условие[9]
(2.7)
являющееся расширением [1, формула (3.6)]. Здесь — первый градиент деформации , а , соответственно, второй градиент.
. Синтезирование локальных деформаций и гипердеформаций. Согласно второй из формул [1, (3.5)], в которой положим , первый градиент деформации в точке имеет представление
(2.8)
Из него, согласно равенству [1, (3.8)], синтезируется поле локальных деформаций , с разложением
(2.9)
Вместе с тем, поскольку рассматривается среда второго градиента, одного лишь поля недостаточно для описания локальной разгрузки. Нужно еще поле, представляющее второй градиент. С этой целью определим вторые градиенты деформаций :
которые, в соответствии с представлением (2.8) и формулой (2.6), имеют разложение
(2.10)
Синтезируем теперь новый тензор по формуле
(2.11)
и назовем его локальной гипердеформацией. Совокупность таких тензоров образует тензорное поле , которое предположим гладким.
Полагая теперь в формуле (2.10) , приходим к следующему представлению тензора гипердеформаций (2.11):
(2.12)
где
(2.13)
а — компоненты разложения (2.9). Заметим, что в силу перестановочности повторных производных, числа обладают свойством симметрии: .
Таким образом, локальная разгрузка среды второго градиента характеризуется парой тензорных полей , значения которых определяются согласно разложениям (2.9) и (2.12). При фиксированных криволинейных координатах им соответствует поле матриц ,
где
.
Кроме того, в соответствии с (2.7), выполнено свойство:
(2.14)
. Восстановление семейства деформаций. Подобно случаю первого градиента, рассмотрим задачу восстановления семейства деформаций по заданной паре гладких тензорных полей и , первое из которых имеет в качестве значений обратимые линейные преобразования, а второе — симметричные тензоры. Кроме того, предполагается, что эти поля удовлетворяют свойству (2.14).
Зафиксировав точку , изменим, если нужно, криволинейные координаты так, чтобы точке отвечали их нулевые значения[10]. Пусть и — соответствующие функции из разложений (2.9) и (2.12). Тогда относительно пары координат и определим отображение согласно правилу:
(2.15)
где — фиксированная тройка чисел.
Будучи заданным как квадратичная форма, построенное отображение является гладким. Кроме того, в силу обратимости матрицы , теорема об обратной функции [17] гарантирует существование окрестности нуля, в которой (2.15) является диффеоморфизмом. Поэтому, если через обозначить соответствующую окрестность точки в , то придем к деформации формы , являющейся частью формы , в некоторую другую форму . Отображение (2.15) является координатным представлением ; кроме того, по построению
что влечет равенства и .
Повторяя проделанную процедуру для всех точек формы , приходим к семейству деформаций , по которому синтезируется пара . Кроме того, для построенного семейства выполняется свойство (2.7), т. е. деформации являются разгрузочными. Вместе с тем в отличие от изначально определенных разгрузочных деформаций, заданных на всей промежуточной форме , полученные деформации определены лишь на ее частях. Это обстоятельство не противоречит общей методологии настоящей работы, поскольку рассматриваемые части конечны и потому также являются формами.
. Совместность локальных деформаций. В случае среды второго градиента будем называть пару тензорных полей, обладающую свойством (2.14), совместной, если 1) существует деформация из промежуточной формы в некоторую форму , для которой выполнено равенство , и 2) справедливо соотношение . Таким образом, для выполняется свойство
вытекающее из (2.7). В этом смысле форма является глобальной натуральной.
Заметим, что если пара совместна, то, в частности, совместно поле локальных деформаций . При этом потенциал [1, формула (3.1)], относительно которого определяется натуральное состояние, индуцирован из :
Вместе с тем из совместности в общем случае не вытекает совместность пары . По этой причине уместно расширить терминологию:
- пара имеет первый порядок совместности, если для некоторой деформации выполнено свойство 1), однако свойство 2) не выполнено,
- пара имеет второй порядок совместности (или просто совместна), если выполнены оба свойства 1) и 2),
- пара несовместна, если ни одно из свойств 1), 2) не выполнено.
Таким образом, условие совместности [1, формула (3.15)] для случая простого материала дополняется вторым условием, представленным свойством 2), что приводит к совокупности равенств:
(2.16)
Их выполнение является необходимым (а в случае односвязности и достаточным) для первого и второго порядка совместности пары . В координатной форме, согласно[11] (2.9) и (2.12),
Вместе с тем первое из условий (2.16) избыточно, если рассматривать второй порядок совместности. Действительно, симметрия функций по нижним индексам дает
Но это и означает, что . Поэтому второму порядку совместности отвечает следующее условие:
(2.17)
2.3. Геометрическая интерпретация условия совместности
9. Поле . Условие совместности (2.17) представим в следующем виде. Принимая во внимание, что значения поля являются обратимыми линейными преобразованиями, домножим обе части (2.17) на и просуммируем по :
(2.18)
Определим скалярные поля в соответствии с [1, формула (3.18)] и новые поля по формуле
(2.19)
Тогда равенство (18) принимает вид
(2.20)
В первой части работы [1, предложение 1] уже было показано, что поля являются коэффициентами некоторой связности. Установим, что аналогичное свойство выполняется и для полей .
. Свойства поля . Доказательство того, что поля являются коэффициентами некоторой связности, предварим определением закона преобразования полей . Для этого продифференцируем обе части соотношения [1, формула (3.20)] по новым координатам . Если через обозначить это дифференцирование, то полученная формула примет вид
(2.21)
Полагая в выражении (2.21) и принимая во внимание определение (2.13), получаем искомый закон преобразования полей :
(2.22)
Обозначая через многообразие, над которым определена промежуточная форма , приходим к следующему утверждению:
Предложение 1. Скалярные функции являются коэффициентами некоторой аффинной связности на многообразии .
Доказательство. Докажем утверждение методом, аналогичным тому, как было доказано предложение 1 в [1]. Если через обозначить функции (2.19), определенные относительно координат , то в соответствии с формулами [1, формула (3.21)] и (2.22) получаем
Последнее выражение приводит к закону преобразования коэффициентов связности в координатном репере, что и доказывает предложение.
Связность, соответствующая полям , обозначается через . Из симметрии полей по нижним индексам следует, что и полученная связность симметрична, т. е. ее кручение равно нулю: .
Получим другие выражения для функций . С этой целью восстановим по паре синтезирующее ее семейство деформаций . В свою очередь, семейству отвечает семейство первых градиентов и обратных градиентов с разложениями
соответственно. Их компоненты связаны соотношениями
(2.23)
Дифференцируя обратный градиент в точке , принадлежащей области определения деформации , получаем равенство
(2.24)
Заметим, что рассуждения, приводящие к последнему равенству, аналогичны использованным при выводе формулы (2.6). Полагая теперь , приходим к тензору третьего ранга
который по аналогии с имплантом [1, формула (3.12)] назовем гиперимплантом. Из (2.24) следует разложение для гиперимпланта:
Здесь
(2.25)
а — компоненты поля имплантов .
Для того чтобы связать поля и , рассмотрим второе из равенств (2.23). Дифференцируя обе его части по , получаем
что при дает
Из последнего равенства и вытекает, что
(2.26)
Следовательно, наряду с парой , состоящей из полей локальных деформаций и гипердеформаций, можно рассматривать пару , состоящую из импланта и гиперимпланта. В терминах элементов этой пары связность определяется формулой [1, (3.26)], а связность , в силу равенства (2.26), — формулой
(2.27)
Формулу (2.27) можно записать иначе. Для этого, используя правило дифференцирования сложной функции, получим
в произвольной точке . Переходя теперь к и учитывая (2.25), придем к равенству
(2.28)
где
Наконец, учитывая равенство (2.28), формулу (2.27) можно преобразовать к виду:
В таком виде связность определяется в монографии [6, формула (2.59)].
. Тензор неоднородности. Связности и определяют новое поле
(2.29)
с компонентами . Поле (2.29) является тензорным полем третьего ранга, в чем можно убедиться, сопоставив законы преобразования коэффициентов связностей и : слагаемые, соответствующие вторым частным производным, при вычитании взаимно уничтожаются. Следуя терминологии, предложенной в работе [6], назовем поле тензором неоднородности.
В общем случае тензор неоднородности не является симметричным, а его антисимметричная часть
(2.30)
пропорциональна кручению связности . Действительно,
что влечет
Таким образом, тензор (2.29) можно представить в виде
(2.31)
где — соответствующая симметричная часть.
. Условие совместности в терминах тензора неоднородности. Согласно определению (2.29) тензора неоднородности, условие совместности (2.20) может быть записано в окончательном виде:
(2.32)
Это необходимое (а в случае односвязности и достаточное) условие совместности пары .
Заметим, что в случае выполнения условия (2.32) из равенства (2.31) следует, что , т. е. поле локальных деформаций совместно. Этот результат находится в полном соответствии с установленном выше свойством, согласно которому совместность второго порядка влечет совместность первого порядка. Таким образом, в терминах полей кручения и неоднородности предложенную классификацию совместности деформаций можно записать следующим образом:
- пара имеет первый порядок совместности, если , но ,
- пара имеет второй порядок совместности, если ,
- пара несовместна, если и .
2.4. Неевклидова отсчетная форма
. Материальная связность. Связности Вайценбока [1, формула (3.18)], Леви-Чивита [1, формула (3.32)] и (2.19) могут служить примерами материальных связностей на . Каждой из них отвечает своя неевклидова отсчетная форма. Вместе с тем можно рассматривать и комбинации этих трех связностей, образуя разнообразные поля. Действительно, например, подберем такую связность , чтобы были выполнены следующие условия:
- кручение совпадает с кручением ,
- кривизна отлична от нуля,
- в случае совместных деформаций связность является евклидовой связностью.
Простейшим вариантом является связность, определяемая формулой
(2.33)
Здесь — тензор неоднородности (2.29). В явном виде, принимая во внимание, что компоненты материальной метрики определены равенством , а также используя формулы [1, (3.18)] и (2.19), формулу (2.33) можно записать как
(2.34)
Вклад второго градиента в геометрическую структуру представлен полями , которые наряду с определяют локальную гипердеформацию .
Поскольку первое слагаемое в (2.33) является связностью, а второе — тензором третьего ранга, то, составляя закон преобразования для полей , можно убедиться в том, что они действительно являются коэффициентами связности на . Полученная связность удовлетворяет требованиям а)– в), указанным выше. В самом деле, поскольку связности и симметричны, то
то есть кручение связности совпадает с кручением связности . Кроме того, кривизна связности отлична от нуля, будучи представленной в компонентах выражением
Здесь — кривизна, порожденная тензором Леви-Чивита, а — формальное выражение, которое представляет «кривизну», порожденную тензором . Наконец, если деформации совместны, то и потому связность совпадает со связностью Леви-Чивита. Но связность Леви-Чивита в этом случае сама порождается совместными деформациями, и потому ее кривизна равна нулю. Следовательно, в случае совместных деформаций является евклидовой связностью.
Связность можно в определенном смысле рассматривать как модификацию связности Леви-Чивита для случая второго градиента. Аналогичная модификация связности Вайценбока приводит к полю
(2.35)
В явном виде поле (2.35) представлено выражением
(2.36)
Рассуждениями, аналогичными тем, что касались поля , можно установить, что поле также является связностью на . Эта связность удовлетворяет свойствам б) и в), и с точностью до множителя удовлетворяет свойству а):
В рамках физических приложений геометрия, определяемая связностью (2.33) (или (2.35)), могла бы быть полезна в случае одновременного наличия дислокаций, дисклинаций и точечных дефектов[12] в кристалле. Отметим, что, по-видимому, первое указание на то, что в рамках второго градиента можно эффективно смоделировать наличие дислокаций и дисклинаций, содержится в работе [22].
. Структура неевклидовой формы. Связность , определенная согласно формуле (2.33), материальная метрика , заданная формулой [1, формула (3.28)], и соответствующая метрике форма объема задают на многообразии структуру пространства с неевклидовой связностью:
(2.37)
Пространство является искомой неевклидовой отсчетной формой над материальным многообразием, соответствующей теории второго градиента. Его геометрия характеризуется полями кручения , кривизны и неметричности . Подобно тому как это было сделано в работе [10], можно показать, что эти поля не зависят от выбора промежуточной формы .
Заменяя материальную связность (2.33) на связность (2.35), от неевклидовой формы (2.37) можно перейти к неевклидовой форме
(2.38)
Форму также следует считать искомой, поскольку в рамках построений настоящей статьи нет никакого способа предпочесть форму (2.37) форме (2.38). Следовательно, здесь возникает та же ситуация, что и для простого материала, где в качестве материальной можно выбирать как связность Вайценбока, так и связность Леви-Чивита.
2.5. Пример синтезирования неевклидовой формы
. Семейство деформаций. Проиллюстрируем рассуждения на примере центрально-симметричного деформирования промежуточной формы
, (2.39)
развивая соответствующий пример из [1]. Именно предположим, что определено семейство деформаций , для которого выполнены условия:
а) при каждом значении параметра деформация центрально-симметрична. Кроме того, для всех точек сферы справедливо свойство:
б) в сферических координатах каждая деформация имеет представление [1, (3.35)].
. Синтезирование поля локальных гипердеформаций. Используя матрицу компонент первого градиента [1, формула (3.37)], получим гиперматрицу частных производных от этих компонент. Ей соответствуют следующие три -матрицы:
Полагая теперь, в соответствии с (2.13), и заменяя затем все вхождения на , приходим к матрицам:
(2.40)
которыми, наряду с матрицей [1, (3.38)], определяется поле локальных гипердеформаций .
. Синтезирование геометрии. Метрический тензор , связность Леви-Чивита и связность Вайценбока синтезированы в рамках примера из статьи [1] и представлены формулами [1, (3.39)], [1, (3.40)] и [1, (3.42)] соответственно. В свою очередь, связность определяется по (2.19) и в силу формул [1, (3.38)] и (2.40) ее отличные от нуля коэффициенты имеют вид:
Полученным выражениям для коэффициентов связностей и отвечает тензор неоднородности с компонентами
(2.41)
Следовательно, подобно случаю первого градиента, условие совместности поля локальных деформаций представлено равенством .
Подстановка соотношений [1, (3.40)] и (2.41) в общую формулу (2.33) приводит к следующим выражениям для коэффициентов связности :
(2.42)
Кручение связности совпадает с кручением связности и представлено выражениями [1, формула (3.43)]. Помимо этого, в рамках рассматриваемой модельной задачи кривизна связности совпадает с кривизной связности Леви-Чивита . В частности, скалярная кривизна представлена выражением [1, формула (3.41)]. Наконец, тензор неметричности связности [1, формула (3.34)] характеризуется следующими компонентами:
Соотношения [1, (3.39)] и (2.42) определяют неевклидову форму (2.37) частного вида. Мера несовместности деформаций, представленная тензором , одновременно характеризует отклонение геометрии от евклидовой.
3. Синтезирование неевклидовой отсчетной формы для среды с микроструктурой
3.1. Микроморфный континуум
. Деформации. В заключение рассмотрим микроморфный континуум — среду, точки которой, наряду с трансляционными степенями свободы, имеют дополнительные степени свободы, т. е. обладают внутренней структурой. Деформация континуума представлена парой , где — макродеформация, соответствующая деформации в рамках теории простого материала, а — микродеформация, являющаяся полем обратимых линейных отображений, характеризующим изменение внутренней структуры частиц среды. Предположим, что упругий потенциал относительно промежуточной формы представлен отображением [23]
(3.1)
Тогда отклик тела в точке характеризуется равенством
где — градиент макродеформации в точке . Отметим, что в частном случае, если положить , то придем к среде второго градиента.
. Координатное представление деформаций. В паре криволинейных координат и , ассоциированных с формами и соответственно, градиент макродеформации представлен разложением [1, (3.4)]. Аналогичным образом в виде двухточечного разложения представим микродеформацию :
Тогда, в соответствии с формулой (2.3), для градиента справедливо разложение
В дальнейшем в качестве координат выбираются прямоугольные координаты .
Замечание 2. Тело с микроструктурой деформируется более сложным образом, чем простой материал. Как и в случае последнего, деформация тела с микроструктурой сопровождается изменением формы, т. е. области физического пространства, занимаемой телом. Это изменение определяется отображением . Однако (и в этом состоит отличие от простого материала) наряду с изменением формы меняются и положения элементов внутри микроструктуры [24], моделируемой касательным пространством к в рамках первого приближения. Именно поэтому, независимо от , задано двухточечное тензорное поле , определяющее изменение положений элементов микроструктуры.
3.2. Семейство форм и условие совместности
. Гипотеза локальной разгрузки. В случае микроморфного континуума гипотеза локальной разгрузки принимается в следующем виде. Пусть фиксированы тензоры и , и пусть выбрана промежуточная форма вместе с упругим потенциалом (3.1). Далее предположим, что определено семейство пар , состоящих из макродеформаций и микродеформаций , такое, что в любой точке выполнено условие
(3.2)
Здесь — градиент макродеформации .
. Синтезирование локальных макродеформаций и микродеформаций. Опираясь на определение [1, (3.8)], синтезируем поле локальных макродеформаций . В каждой точке локальная макродеформация имеет разложение [1, (3.11)], совпадающее с разложением локальной деформации в случае среды первого градиента.
По семейству синтезируем новые тензорные поля. Прежде всего, положим
(3.3)
Совокупности обратимых тензоров отвечает поле , которое будем полагать гладким. В координатном представлении,
(3.4)
Далее, синтезируем тензор третьего ранга
и назовем его локальной микродеформацией. Приходим к полю , которое также предположим гладким. Отметим, что в отличие от второго градиента макродеформации , это поле в общем случае не является симметричным по нижним индексам. Формула (2.3) влечет следующее представление для , подобное (2.12):
(3.5)
в котором, однако, величины
(3.6)
уже несимметричны по нижним индексам.
Таким образом, локальная разгрузка микроморфной среды характеризуется тройкой тензорных полей , определяемых разложениями [1, (3.11)], (3.4) и (3.5). При фиксированных криволинейных координатах им соответствует поле матриц ,
где
Кроме того, в соответствии с (3.2) выполнено свойство:
(3.7)
во всех точках .
. Совместность локальных деформаций. Назовем тройку тензорных полей, обладающих свойством (3.7), совместной, если 1) существует деформация из промежуточной формы в некоторую форму , для которой выполнено равенство , и 2) существует поле обратимых линейных преобразований такое, что и . В этом случае для пары выполняется свойство
вытекающее из (3.7), и потому форма является глобальной натуральной.
Как и в случае среды второго градиента, если тройка совместна, то, в частности, совместно поле локальных деформаций . Потенциал [1, (3.1)], относительно которого определяется натуральное состояние, индуцирован из :
Однако из совместности в общем случае не вытекает совместность всей тройки . Подобно среде второго градиента, можно соответствующим образом расширить терминологию.
На основании условий 1) и 2), условие совместности для среды с микроморфной кинематикой можно представить в виде совокупности равенств:
В координатной форме, согласно [1, (3.11)], (3.4) и (3.5),
(3.8)
(3.9)
Заметим, что в отличие от случая второго градиента полученные соотношения независимы.
3.3. Геометрическая интерпретация условия совместности
. Поля и . Выполним следующее преобразование условия совместности (3.8)–(3.9). Поскольку матричные поля и обратимы, домножим обе части равенства (3.8) на , а обе части (3.9) — на , произведя в обоих случаях суммирование по . Тем самым придем к следующей совокупности соотношений, эквивалентных исходным:
(3.10)
Задавая поля по формуле [1, формула (3.18)] и вводя новые поля по формуле, аналогичной [1, формула (3.18)],
(3.11)
и поля по аналогии с (2.19) как
(3.12)
можно переписать равенства (3.10) в лаконичном виде:
(3.13)
. Кручение и тензор неоднородности. Каждое из полей , и задает связность на . Действительно, для первого поля это уже было установлено при рассмотрении среды первого градиента; второе поле рассматривается аналогично. В случае третьего поля можно дословно повторить доказательство предложения 1. Отметим, что, в отличие от своего аналога из теории второго градиента, поля (3.12) несимметричны.
Таким образом, равенства (3.13) можно записать как
(3.14)
где — кручение связности , а
(3.15)
— аналог тензора неоднородности. Тем самым показано, что в случае микроморфного континуума меры несовместности представлены независимыми полями кручения и неоднородности .
Тензор неоднородности не симметричен по нижним индексам, а его антисимметричная часть с компонентами (2.30) пропорциональна разности кручений связностей и . Действительно,
что и приводит к равенству
3.4. Неевклидова отсчетная форма
. Материальная связность. Связности Вайценбока [1, формула (3.18)], Леви-Чивита [1, формула (3.32)], и новые связности (3.11), (3.12) служат примерами материальных связностей на , каждой из которых отвечает своя неевклидова отсчетная форма. Наряду с этим, можно рассматривать и комбинации этих связностей. Примером может служить связность , заданная формулой
(3.16)
аналогичной (2.33). В явном виде, принимая во внимание формулу (3.15), а также определения (3.11) и (3.12), получаем выражение
(3.17)
Кручение связности определяется равенством
а кривизна в общем случае отлична от нуля.
Альтернативную связность на можно получить, если вместо связности Леви-Чивита использовать связность Вайценбока . Действительно, положим
(3.18)
По форме построенная связность напоминает связность (2.35), однако, если раскрыть определения входящих в нее полей, то выражение для ее коэффициентов имеет иной вид:
(3.19)
Кручение связности (3.18) определяется формулой
Заметим, что если локальные деформации совместны (т. е. выполнены равенства (3.14)), то связность, определяемая формулой (3.16) (или (3.18)), является евклидовой. Поэтому она также может быть полезна для моделирования несовместных деформаций, возникающих в средах с дефектами.
Замечание 3. В частном случае, когда микродеформация полагается равной градиенту деформации, теория среды с микроструктурой сводится к теории среды второго градиента. Тогда формула (3.17) переходит в формулу (2.34), а формула (3.19) — в формулу (2.36).
. Структура неевклидовой формы. Связность (3.16), материальная метрика [1, формула (3.28)], и соответствующая метрике форма объема вводят на многообразии структуру пространства с неевклидовой связностью:
(3.20)
Пространство является примером неевклидовой отсчетной формы над материальным многообразием, соответствующей континууму с микроструктурой.
Если вместо материальной связности (3.16) использовать связность (3.18), то придем к альтернативному представлению неевклидовой формы:
(3.21)
Как и форма (3.20), синтезированная форма (3.21) в равной степени пригодна для моделирования глобального натурального состояния тела с микроструктурой. Выбор определенной формы требует иных соображений, выходящих за рамки статьи.
3.5. Расширенное описание для тела с микроструктурой
. Материальное и физическое расслоения. В силу наличия дополнительного тензорного поля , не связанного с преобразованием форм , кинематика сред с микроструктурой отличается от кинематики простого материала и сред второго градиента. В этой связи целесообразно поставить вопрос о том, можно ли, поменяв соответствующим образом описание форм, представить деформацию среды с микроструктурой как изменение форм, понимаемых в расширенном смысле. Ответ на этот вопрос положительный.
Действительно, естественным формализмом, пригодным для построения расширения форм, является теория расслоений, восходящая к работам Фельдбау, Эресманна и Стинрода[13] [27–29]. В рамках этого формализма над отсчетной формой и физическим пространством выстраиваются тотальные пространства локально тривиальных расслоений с одним и тем же типовым слоем . В явном виде эти расслоения определены структурами [30]
(3.22)
— материальное расслоение и
(3.23)
— физическое расслоение, где и — гладкие многообразия, представляющие тотальные пространства расслоений, а и — гладкие сюръективные отображения. Многообразие — общий элемент структур и — представляет типовой слой расслоения.
Предполагается, что расслоения (3.22) и (3.23) локально тривиальны. Это означает (случай расслоения идентичен рассмотренному далее), что для любой точки существуют окрестность в и диффеоморфизм , называемый локальной тривиализацией, такие, что , где — проекция декартова произведения на первый сомножитель. Таким образом, локально диффеоморфизм должен иметь представление для некоторого гладкого отображения , свойства которого зависят от частного вида расслоения. Примеры материальных расслоений и их физическую интерпретацию можно найти в работе [31].
. Дополненные деформации. Тотальные пространства и расслоений (3.22) и (3.23) над отсчетной формой и физическим пространством определяют дополненную отсчетную форму и дополненное физическое пространство соответственно. Каждой форме соответствует дополненная форма , определяемая как прообраз
Дополненная деформация определяется как отображение , удовлетворяющее свойству , где — классическая деформация. Таким образом, в специально подобранных локальных координатах на тотальных пространствах можно записать
Здесь — локальные координаты на форме , а — локальные координаты на типовом слое . Тильдами сверху обозначены координатные представления соответствующих отображений. Отображение определяет тип микроструктуры.
В рамках описания микроструктуры, используемого в настоящей работе, в качестве расслоений (3.22) и (3.23) достаточно выбрать касательные расслоения. Несмотря на их тривиальность в случае трехмерных тел, они перестают быть тривиальными при рассмотрении материальных поверхностей. При таком выборе представление дополненной деформации есть не что иное, как пара [31] . Этим получен желаемый результат: формализм касательных расслоений позволяет интерпретировать пару как одно точечное отображение между касательными пространствами.
Использование формализма расслоений приводит к необходимости модифицировать геометрический язык для описания неевклидовой отсчетной формы. Вместо теории аффинных связностей применяется более общая теория связностей на расслоениях. Подробно этот вопрос рассмотрен в работе [32].
3.6. Пример синтезирования неевклидовой формы
. Семейство деформаций. Завершим основные построения настоящего раздела модельным примером. Предположим, что промежуточная форма является полым шаром (39) и, что на ней задано семейство , состоящее из макродеформаций и микродеформаций . Это семейство подчиняется следующим условиям:
а) для любого макродеформация является центрально-симметричной и имеет представление [1, (3.35)] в сферических координатах ;
б) для каждого микродеформация в точках сферы является поворотом материальных элементов вокруг базисного вектора на угол ;
в) для всех точек сферы справедливы соотношения:
Явный вид тензора , задающего поворот вокруг оси, определяемой единичным вектором , на угол , может быть определен в соответствии с формулой Родрига [33; 34]. Относительно диадного базиса тензор представлен матрицей
где — единичная матрица, а
— матрица, составленная из компонент вектора . Тогда матрица относительно смешанного диадного базиса определяется по формуле
где — матрица перехода от базиса к базису в точке .
Окончательно
(3.24)
. Синтезирование поля локальных деформаций. Поле локальных макродеформаций определено в [1] и представлено формулой [1, (3.38)]. Осталось определить поля и , соответствующие локальным микродеформациям. Первое из них синтезируется по формуле (3.3). Полагая в (3.24) и заменяя затем все переменные на , получаем
(3.25)
Для синтезирования поля требуется еще определить скалярные поля в соответствии с (3.6). С этой целью получим гиперматрицу , дифференцируя элементы матрицы (3.24). В результате придем к следующим матрицам:
Из них при вытекают представления для искомых функций:
(3.26)
. Синтезирование геометрии. Поскольку метрический тензор , связность Леви-Чивита и связность Вайценбока определяются по полю локальных макродеформаций , то их выражения известны и представлены формулами [1, (3.39), (3.40), (3.42)] соответственно. Новыми являются выражения для связностей , и результирующей связности .
Коэффициенты связности определяются в соответствии с общей формулой (3.11). Принимая во внимание выражение (3.25) для поля , приходим к равенствам (приведены только те функции, которые отличны от нуля)
(3.27)
Из представленных выражений следует, что связность несимметрична, т. е. обладает ненулевым кручением. Отметим, что последнее отлично от кручения связности .
Связность определяется согласно формуле (3.12). Учитывая соотношение (3.25) для поля и формулу (3.26), определяющую частный вид полей , получаем следующие соотношения:
(3.28)
Несмотря на то что большинство соответствующих коэффициентов связностей и совпадает, между ними все же есть отличные друг от друга. Например, , в то время как . Мера отличия связностей друг от друга представлена тензором неоднородности (3.15), отличные от нуля компоненты которого, в соответствии с (3.27) и (3.28), имеют вид:
(3.29)
Следовательно, с учетом выражений [1, (3.43)] для кручения связности условия совместности (3.14) локальных деформаций в рассматриваемом модельном случае имеют вид
Подстановка соотношений [1, (3.40)] для коэффициентов связности Леви-Чивита и равенств (3.29) в общую формулу (3.16) приводит к следующим выражениям для коэффициентов связности :
(3.30)
Кручение связности отлично от нуля и представлено следующими компонентами:
Кривизна связности также отлична от нуля. В частности, скалярная кривизна имеет вид [1, (3.41)]. Неметричность равна нулю.
Формулы [1, (3.39)] и (3.30) определяют неевклидову форму (3.20). Меры несовместности деформаций, представленные полем кручения связности и тензором неоднородности , одновременно характеризуют отличие геометрии от евклидовой.
Заключение
Геометрический метод, изложенный в статье [1], обобщен в настоящей работе на тела с расширенной кинематикой. Показано, что благодаря дополнительным полям и , определяемым спецификой расширенной кинематики, в таких телах имеется больше вариантов для синтезирования геометрии. Сводка результатов представлена в таблице.
Таблица. Различные материальные связности
Table. Various material connections
| ||||||
Простой материал | [1, (3.18)] | [1, (3.32)] | - | - | - | - |
Второй градиент | [1, (3.18)] | [1, (3.32)] | - | (2.19) | (2.34) | (2.36) |
Микроструктура | [1, (3.18)] | [1, (3.32)] | (3.11) | (3.12) | (3.17) | (3.19) |
Ее строки содержат ссылки на явные представления материальных связностей. Во всех трех случаях кинематики, рассмотренных в работе, связности Вайценбока и Леви-Чивита вычисляются по одним и тем же формулам. Связности и могут быть определены лишь для тел с расширенной кинематикой, поскольку вычисление их коэффициентов требует знания дополнительных полей и . Наконец, связности и служат расширениями связностей Леви-Чивита и Вайценбока. В рамках физических приложений они могут быть полезны для геометрического моделирования тел с дислокациями и дисклинациями, поскольку их кручение и кривизна — плотности дислокаций и дисклинаций, — отличны от нуля.
Благодарности
Автор работы благодарит С.А. Лычева и А.Л. Левитина за обсуждения и критические замечания по содержанию статьи.
1 См. также [1, и ].
2 В рамках классической механики континуума нет необходимости выделять третью форму, наряду с отсчетной и актуальной формами. Вместе с тем в настоящей работе третья форма, называемая промежуточной, используется для синтезирования геометрии на материальном многообразии. Формально, подобно любой другой форме тела (см. [1, 7]), промежуточная форма представлена в виде упорядоченной совокупности , где — подлежащее многообразие (носитель геометрической структуры), а , , обозначают метрику, форму объема и связность над многообразием , индуцированные из евклидова пространства . При синтезировании неевклидовой отсчетной формы последние три поля заменяются на более общие геометрические поля.
3 Хотя в оригинальных работах Картана рассматривались только вращения, можно предположить, что эти преобразования являются достаточно общими.
4 В явном виде этот изоморфизм может быть определен следующим образом. Если — элемент , то ему отвечает элемент пространства , заданный равенством . Искомый изоморфизм является соответствием .
5 Здесь обозначает декартову арифметизацию [1, формула (2.2)].
6 Отображение есть координатное отображение, соответствующее криволинейным координатам на форме [1, формула (2.13)].
7 Формула (2.4) — производная тензорного поля — содержит производные скалярных и векторных базисных полей. Производные базисных полей, в свою очередь, выражаются через исходные базисные поля; коэффициенты разложения есть символы Кристоффеля соответствующей системы координат.
8 Поскольку физическое пространство евклидово, а размерности форм и пространства совпадают.
9 Несмотря на то что в рамках настоящей работы тензоры и предполагаются постоянными, ход рассуждений не изменится, если эти тензоры заменить на переменные тензоры, зависящие от точки . В таком случае правые части соотношений (2.7) и им подобным следует заменить на и соответственно. Такая замена целесообразна, если натуральное состояние определяется в рамках действия некоторого внешнего поля, которым нельзя пренебречь.
10 Это всегда можно сделать, используя операцию сдвига в .
11 Для вычисления производной можно воспользоваться равенством (2.3), где , а .
12 Поскольку неметричность [1, формула (3.34)] связности в общем случае отлична от нуля.
13 По-видимому, первое определение расслоения было дано в работе Уитни [25], хотя частные случаи использовались и ранее, например, Зейфертом [26].
About the authors
Konstantin G. Koifman
Bauman Moscow State Technical University
Author for correspondence.
Email: koifman.konstantin@gmail.com
ORCID iD: 0000-0002-7891-9995
tutor in mathematics
Russian Federation, 5, 2nd Baumanskaya Street, Moscow, 105005References
- Koifman K.G. Reference shape of bodies with enhanced kinematics. Part I. Geometric methods. Vestnik Samarskogo Universiteta. Estestvennonauchnaya Seriya Vestnik of Samara University. Natural Science Series, 2023. vol. 29, no. 4. pp. 26–53.
- Kröner E. Allgemeine Kontinuumstheorie der Versetzungen und Eigenspannungen. Archive for Rational Mechanics and Analysis, 1959, vol. 4, no. 4, pp. 18–334. Available at: https://libcats.org/book/789336. (In Russ.)
- Noll W. Materially uniform simple bodies with inhomogeneities. Archive for Rational Mechanics and Analysis, 1967, Vol. 27, No. 1, Pp. 1–32. DOI: https://doi.org/10.1007/BF00276433.
- Wang C.-C. On the geometric structures of simple bodies, a mathematical foundation for the theory of continuous distributions of dislocations. Archive for Rational Mechanics and Analysis, 1967, Vol. 27, No. 1, Pp. 33–94. DOI: https://doi.org/10.1007/BF00276434.
- Miri M., Rivier N. Continuum elasticity with topological defects, including dislocations and extra-matter. Journal of Physics A: Mathematical and General, 2002, Vol. 35, Number 7, Pp. 1727–1739. DOI: https://doi.org/10.1088/0305-4470/35/7/317.
- Epstein M., El anowski M. Material inhomogeneities and their evolution: A geometric approach. Berlin Heidelberg: Springer Science & Business Media, 2007. 261 p. DOI: https://doi.org/10.1007/978-3-540-72373-8.
- Yavari A., Goriely A. Riemann–Cartan Geometry of Nonlinear Dislocation Mechanics. Archive for Rational Mechanics and Analysis, 2012, Vol. 205, no. 1, Pp. 59–118. DOI: https://doi.org/10.1007/s00205-012-0500-0.
- Lychev S.A., Manzhirov A.V. The mathematical theory of growing bodies. Finite deformations. Journal of Applied Mathematics and Mechanics, 2013, vol. 77, issue 4, pp. 421–432. DOI: https://doi.org/10.1016/j.jappmathmech.2013.11.011. EDN: https://elibrary.ru/wqyump. (In English; original in Russian)
- Lychev S.A., Koifman K.G. Contorsion of Material Connection in Growing Solids. Lobachevskii Journal of Mathematics, 2021, Vol. 42, No. 8, Pp. 1852–1875. DOI: https://doi.org/10.1134/S1995080221080187.
- Lychev S.A., Koifman K.G. Reference shape of bodies with finite incompatible deformations. Vestnik of Samara University. Natural Science Series, 2022, vol. 28, no. 3–4, pp. 53–87. DOI: https://doi.org/10.18287/2541-7525-2022-28-3-4-53-87. (In Russ.)
- Kartan E. Spaces of affine, projective and conformal connection. Kazan: Izd-vo Kazanskogo universiteta, 1962, 210 p. Available at: https://libcats.org/book/444677. (In Russ.)
- Kartan E. Riemannian geometry in an orthogonal frame. Moscow: Nauka, 1960, 207 p.. URL: https://knigogid.ru/books/1911053-rimanova-geometriya-v-ortogonalnom-repere/toread?ysclid=lp893xl1f8549212671. (In Russ.)
- Kartan E. Geometry of Riemannian spaces. Moscow: Knizhnyi dom «Librokom», 2010, 248 p. Available at: https://reallib.org/reader?file=444675. (In Russ.)
- Lee J.M. Introduction to Smooth Manifolds. New York: Springer, 2012, 708 p. DOI: https://doi.org/10.1007/978-1-4419-9982-5.
- Toupin R.A. Theories of elasticity with couple-stress. Archive for Rational Mechanics and Analysis, 1964, vol. 17, pp. 85–112. DOI: https://doi.org/10.1007/BF00253050.
- Kalpakides V., Agiasofitou E. On Material Equations in Second Gradient Electroelasticity. Journal of Elasticity, 2002, vol. 67, pp. 205–227. DOI: https://doi.org/10.1023/A:1024926609083.
- Schwartz L. Analysis. Vol. 1. Moscow: Mir, 1972, 824 p. Available at: https://klex.ru/14z8. (In Russ.)
- Postnikov M.M. Lections in geometry. Semester II. Linear algebra. Moscow: URSS, 2017, 400 p. Available at: http://alexandr4784.narod.ru/pmmgeo2.html. (In Russ.)
- [Truesdell C., Toupin R. The Classical Field Theories. In: In: Flugge, S. (eds) Principles of Classical Mechanics and Field Theory / Prinzipien der Klassischen Mechanik und Feldtheorie. Encyclopedia of Physics / Handbuch der Physik, vol. 2 / 3 / 1. Springer, Berlin, Heidelberg. DOI: https://doi.org/10.1007/978-3-642-45943-6_2.
- Maugin G.A. Material inhomogeneities in elasticity. New York: CRC Press, 1993, 292 p. DOI: https://doi.org/10.1201/9781003059882.
- Marsden J.E., Hughes T.J. Mathematical foundations of elasticity. New York: Courier Corporation, 1994. 576 p. URL: https://archive.org/details/mathematicalfoun00mars.
- El anowski M., Epstein M. The symmetry group of second-grade materials. International Journal of Non-Linear Mechanics, 1992, vol. 27, issue 4, pp. 635–638. DOI: https://doi.org/10.1016/0020-7462(92)90068-I.
- Lychev S.A. On conservation laws of micromorphic nondissipative thermoelasticity. Vestnik of Samara State University, 2007, no. 4(54), pp. 225–262. Available at: https://www.elibrary.ru/item.asp?id=9505071. EDN: https://elibrary.ru/hzzzon. (In Russ.)
- Mindlin R.D. Micro-structure in linear elasticity. Archive for Rational Mechanics and Analysis, 1964, Vol. 16, Pp. 51–78. DOI: https://doi.org/10.1007/BF00248490.
- Whitney H. Sphere spaces. Proceedings of the National Academy of Sciences, 1935, vol. 21, no. 7, pp. 462–468. DOI: http://doi.org/10.1073/pnas.21.7.464.
- Seifert H. Topologie Dreidimensionaler Gefaserter Räume. Acta Mathematica, 1933, vol. 60, pp. 147–238. DOI: https://doi.org/10.1007/BF02398271.
- Feldbau J. Sur la classification des espaces fibrés. Comptes rendus de l’Académie des Sciences, 1939, vol. 208, pp. 1621–1623. Available at: https://zbmath.org/0021.16304.
- Ehresmann C. Sur les espaces fibrés différentiables. Comptes rendus de l’Académie des Sciences, 1947, vol. 224, pp. 1611–1612. Available at: https://zbmath.org/0029.42001.
- Steenrod N.E. The topology of fibre bundles. Moscow: Izdatel’stvo inostrannoi literatury, 2010, 272 p. Available at: https://reallib.org/reader?file=444395&pg=1. (In Russ.)
- Postnikov M.M. Lections in geometry. Semester IV. Differential geometry. Moscow: URSS, 2017, 504 p. Available at: http://alexandr4784.narod.ru/pmm4.html. (In Russ.)
- Lychev S.A., Koifman K.G., Digilov A.V. Nonlinear dynamic equations for elastic micromorphic solids and shells. Part I. Vestnik of Samara University. Natural Science Series, 2021, vol. 27, issue 1, pp. 81–103. DOI: https://doi.org/10.18287/2541-7525-2021-27-1-81-103.
- Nguyen V.H., Casale G., Le Marrec L. On tangent geometry and generalised continuum with defects. Mathematics and Mechanics of Solids, 2022, Vol. 27, no. 7, Pp. 1255–1283. https://doi.org/10.1177/10812865211059222.
- Lurie A.I. Analytical Mechanics. Moscow: GIFML, 1961, 824 p. Available at: https://eqworld.ipmnet.ru/ru/library/books/Lure1961ru.djvu. (In Russ.)
- Liang K.K. Efficient conversion from rotating matrix to rotation axis and angle by extending Rodrigues’ formula. arXiv, 2018. DOI: https://doi.org/10.48550/arXiv.1810.02999.
Supplementary files
