1. Введение. Современный этап развития биоинженерии, трансплантологии, аддитивного производства биоматериалов ставит перед наукой, в частности перед механикой континуума, проблему построения адекватных математических моделей деформирования материалов со сложной внутренней микроструктурой зернистых, волокнистых, гранулированных, дипольных, полимерных, сотовых и т.д.
Одной из простейших и исторически первой математической моделью, в которой оказывается возможным введение нано/микромасштабной длины, является модель микрополярного упругого континуума, предложенная в работах Коссера [1]. В микрополярной теории упругости трансляционные перемещения каждого элемента континуума сопровождаются микроповоротами (спинорными перемещениями). Спинорные и трансляционные перемещения считаются независимыми [24].
Современные исследования [5] ясно указывают на то, что важной особенностью биоматериалов является чувствительность их определяющих постоянных и тензоров к преобразованиям, изменяющим ориентацию пространства, инверсиям и зеркальным отражениям. Как следствие, возникает необходимость классификации микрополярных упругих моделей с точки зрения изменения ориентации пространства.
Удобной, легко воспринимаемой формой представления тензоров четвертого и третьего рангов посредством блоков двумерных квадратных и прямо угольных матриц и связей между элементами указанных матриц являются обобщенные фигуры Ная [6, с. 113115]. Метод построения фигур Ная описан в его классических работах по кристаллографии [610]. В цикле статей авторов [1113] в матричной форме были построены определяющие уравнения для общего анизотропного, гемитропного и изотропного микрополярного упругого континуума. Показано, что определяющие тензоры для анизотропного случая содержат 171 независимую определяющую константу, гемитропного 9, изотропного 6. В представляемой работе подход, предложенный Наем, используется для получения определяющих уравнений ультраизотропного микрополярного упругого тела, т.е. такого тела, компоненты определяющих тензоров которого инвариантны относительно любых мыслимых преобразований пространства. Показано, что определяющие тензоры данного материала содержат лишь 4 независимые определяющие постоянные: модуль сдвига, коэффициент Пуассона, характерная нано/микродлина и еще одна, не имеющая физической размерности, постоянная.
2. Тензоры/псевдотензоры с постоянными компонентами. Ковариантно постоянные тензорные/псевдотензорные поля играют важную роль как при выводе определяющих уравнений, так и при преобразовании дифференциальных уравнений механики сплошных сред [1418]. Важным типом ковариантно постоянных тензоров/псевдотензоров являются тензоры/псевдотензоры с постоянными компонентами [19, с. 164].
Тензором/псевдотензором с постоянными компонентами [19, с. 164] называется тензор/псевдотензор, сохраняющий неизменными все свои компоненты при любых линейных преобразованиях координатного репера: самые важные из них повороты, преобразования гомотетии, центральная инверсия, зеркальные отражения.
В монографии [19, с. 164176] изложен общий алгоритм построения тензоров/псевдотензоров с постоянными компонентами для целых положительных/отрицательных алгебраических весов. Например, общий вид псевдотензора null с постоянными компонентами целого отрицательного веса представляется формулой:
null (2.1)
где r число ковариантных индексов, s число контравариантных индексов, N размерность пространства, g вес (целое отрицательное число), λP (P = 1, 2, ..., r !) произвольные постоянные (абсолютные инварианты), P перестановка в ряде индексов
В формуле (2.1) по ковариантным индексам, заключенным в фигурные скобки, производятся всевозможные перестановки. Число ковариантных, контравариантных индексов и вес псевдотензора должны удовлетворять ограничению
(2.2)
откуда следует взаимное ограничение
(2.3)
Если условие (2.2) не выполняется, то псевдотензор null с постоянными компонентами сводится к нулевому.
Отметим, что псевдотензорное поле null с постоянными компонентами является ковариантно постоянным и удовлетворяет псевдотензорному уравнению
null (2.4)
при условии, что
(2.5)
Обратим внимание, что псевдотензоры вида (2.1) не составляют полного набора ковариантно постоянных абсолютных тензоров. Примеры ковариантно постоянных тензоров и псевдотензоров подробно обсуждались в работах [1721]. Среди них фундаментальный ориентирующий псевдоскаляр e и его алгебраические степени, псевдотензорные единицы обобщенные d-символы, ò-символы, e-тензоры, метрические тензоры g kh, g hk, которые часто используются в микрополярных теориях механики сплошных сред [2226].
Рассмотрим важный для дальнейшего изложения пример. Абсолютный тензор четвертого ранга Csmil с постоянными компонентами согласно (2.1) можно представить в виде:
(2.6)
где a = C1212 и c = C2112 абсолютные инварианты (абсолютные скаляры).
Уравнение (2.6), справедливое в любой системе координат, в декартовых координатах можно представить следующим образом:
(2.7)
Рассмотрим ограничение (2.2). Тогда для тензора четвертого ранга, заданного в N-мерном пространстве, можно получить систему уравнений для определения числa ковариантных r и контравариантных s индексов в виде:
(2.8)
Выражая r и s, получим:
(2.9)
Решениями системы (2.9) должны быть целые неотрицательные числа, откуда немедленно заключаем, что произведение N | g | должно быть четным неотрицательным целым числом, удовлетворяющим неравенству
(2.10)
Анализируя неравенство (2.10), заинтересованный читатель может заключить, что в трехмерном пространстве и в пространствах с размерностью выше пятой невозможно построить псевдотензор четвертого ранга с постоянными компонентами ненулевого алгебраического веса. Однако в пространствах размерности 2 и 4 такая возможность существует. Все сказанное имеет место только когда [1]
3. Анизотропное микрополярное упругое тело. В дальнейшем все рассуж дения будем проводить в декартовой прямоугольной системе координат. Рассмотрим микрополярный упругий потенциал напряжений в расчете на единицу инвариантного элемента объема [2732], с естественными асимметричными тензорными аргументами (не разделяя пока на симметричную и антисимметричную конституэнты:
(3.1)
где асимметричный тензор деформации; - тензор изгибакручения.
В случае линейного анизотропного микрополярного упругого тела квадратичная энергетическая форма упругого потенциала напряжений записывается в виде [3335]:
null (3.2)
где null (c = 1, 2, 3) определяющие тензоры анизотропного микрополярного упругого тела. Для определяющих тензоров выполняются следующие условия:
1. null и null симметричны при перестановке пар индексов ( ).
2. null не симметричен при перестановке пар индексов ( ).
Выделим симметричные и антисимметричные составляющие асимметричных тензоров деформаций и изгибакручения в виде суммы, имеем:
(3.3)
Определим векторы, ассоциированные с антисимметричными составляющими асимметричных тензоров деформаций и изгибакручения согласно
(3.4)
Формулы, обратные к (1.4), записываются в виде:
(3.5)
Определяющие тензоры четвертого ранга в представлении (3.2) запишем в виде:
null (3.6)
Подстановка соотношений (3.3) с учетом (3.6) в энергетическую форму упругого потенциала (3.2) приводит к формуле:
null (3.7)
Раскрывая скобки и приводя подобные слагаемые в (3.7), имеем десятичленную форму:
(3.8)
Откуда, учитывая (1.5) и принимая во внимание
null (3.9)
после ряда преобразований получим:
(3.10)
Для упрощения дальнейшего изложения материала настоящей статьи введем в рассмотрение следующие обозначения:
null (3.11)
В новых обозначениях (3.11) энергетическая форма упругого потенциала (3.10) примет вид:
(3.12)
4. Гемитропное микрополярное упругое тело. Микрополярное тело называется гемитропным, если компоненты трех его определяющих тензоров четвертого ранга не изменяются при поворотах координатного репера, но, вообще говоря, изменяются при зеркальных отражениях и инверсиях трехмерного Евклидова пространства, т.е. полуизотропны.
В этом случае определяющие тензоры гемитропного микрополярного упругого тела могут быть представлены в виде [36, с. 70]:
null (4.1)
или
null (4.2)
где (c = 1, 2, 3) определяющие постоянные.
С учетом представлений (4.2) соотношения (3.11) запишутся в форме:
null (4.3)
5. Изотропное микрополярное упругое тело. Микрополярное тело называется изотропным, если компоненты его определяющих тензоров не изменяются не только при поворотах координатного репера, но и при зеркальных отражениях и инверсиях трехмерного Евклидова пространства. Последнее означает, что
null (5.1)
Тогда для определяющих тензоров (3.11) получим:
null (5.2)
6. Ультраизотропное микрополярное упругое тело. Изотропное микрополярное упругое тело назовем ультраизотропным, если компоненты двух его определяющих тензоров вообще не изменяются ни при каких преобразованиях трехмерного Евклидова пространства, т.е. являются тензорами с постоянными компонентами.
Учитывая представление (2.7) для тензоров четвертого ранга с постоянными компонентами, определяющие тензоры (5.1) можно представить в виде:
null (6.1)
или
null (6.2)
В обозначениях (3.11) имеем:
null (6.3)
7. Матричное представление определяющих тензоров микрополярного упругого тела. Удобной, легко воспринимаемой формой представления тензоров четвертого, третьего и второго рангов посредством блоков двумерных квадратных и прямоугольных матриц и связей между элементами указанных матриц являются обобщенные фигуры Ная [6, с. 113115]. Метод построения фигур Ная подробно описан в классических монографиях и физических стандартах по кристаллографии [610]. Указанный метод позволяет без труда установить количество независимых определяющих постоянных, характеризующих мик рополярный континуум, и возможных наличия/отсутствия алгебраических связей между ними.
Преобразование определяющих тензоров четвертого, третьего и второго рангов, участвующих в записи микрополярного упругого потенциала (3.12), к виду двумерных матриц можно произвести при помощи соответствующей замены индексов согласно правилу (см. табл. 1), предложенному в монографии [6, с. 113115], т.е. каждой паре тензорных индексов соответствует матричный индекс, обозначенный прописной латинской буквой:
(7.1)
Следует помнить, что метаиндексы a и b не являются тензорными индексами, а нумеруют определяющие тензоры.
Таблица 1. Соответствие пар тензорных и матричных индексов.
пары тензорных индексов (is,lm) | 11 | 22 | 33 | 23, 32 | 31, 13 | 12, 21 |
матричные индексы (K,N) | 1 | 2 | 3 | 4 | 5 | 6 |
Ансамбль двумерных матричных блоков Ная для определяющих тензоров линейного анизотропного микрополярного упругого тела представлен на рис. 1.
Рис. 1. Ансамбль блоков Ная для общих анизотропных микрополярных упругих определяющих тензоров
Фигуры Ная для определяющих тензоров анизотропного, гемитропного и изотропного микрополярных упругих тел приведены в работе [11].
Рассмотрим преобразование определяющих тензоров ультраизотропного микрополярного тела к виду двумерных матриц. В этом случае ненулевыми оказываются два тензора четвертого ранга EI(is)(lm), EII(is)(lm) и два тензора второго ранга EIV(is), EV(is). Заменяя пары индексов у тензоров Ea(is)(lm) (a = I,II) согласно табл. 1, получим:
null (7.2)
Матрицы ненулевых определяющих тензоров второго ранга Eb(is) (b = IV,V) для ультраизотропного микрополярного упругого тела примут вид:
null (7.3)
Метаиндексы a и b в (7.3) связаны соотношением a = b III.
Двумерные клетки Ная определяющих тензоров ультраизотропного мик рополярного упругого тела:
1. null представлены на рис. 2.
2. null представлены на рис. 3.
Рис. 2. Клетки Ная определяющих тензоров ультраизотропного микрополярного упругого континуума: a) null, b) null
Фигура Ная для ультраизотропного микрополярного упругого тела составляется из клеток (см. рис. 2 и 3). В результате несложных геометрических преобразований получим фигуру Ная (см. рис. 4).
Обозначения на фигурах Ная: жирными отрезками соединены равные компоненты, Z означает симметрию числовых значений элементов фигуры относительно главной диагонали, ● отличные от нуля компоненты, компоненты, вычисляемые по следующему правилу null " компоненты, вычисляемые по следующему правилу null нулевые компоненты.
Заключение. В работе выполнено обобщение метода построения двумерных фигур Ная для микрополярных упругих тел. Указанный подход позволяет наиболее простым образом выяснить количество независимых определяющих констант для микрополярных упругих тел и установить возможные связи между ними. В рамках настоящего исследования построены двумерные фигуры Ная для ультраизотропного микрополярного упругого тела. Показано, что определяющие тензоры данного материала содержат лишь 4 независимые определяющие постоянные: модуль сдвига, коэффициент Пуассона, характерная нано/микродлина и еще одна, не имеющая физической размерности, постоянная.
Финансирование. Работа выполнена по теме государственного задания ( госрегистрации 124012500437-9).