Reference shape of bodies with enhanced kinematics. Part II. Second gradient and microstructure

Cover Page

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. Настоящая статья продолжает работу [1], в которой представлен геометрический метод моделирования несовместных деформаций для гиперупругих тел и рассмотрены его особенности на примере простого материала. Несмотря на то что неевклидова отсчетная форма для тела из простого материала является классической и различные способы ее построения рассмотрены в статьях [2–10], подход к синтезированию неевклидовой формы, предложенный в исследовании [1], несколько отличается от этих способов. Он является комбинацией рассуждений Кренера [2], в рамках которых геометрия определяется на основе условий совместности, с идеей локальной разгрузки, предложенной в работах[1] [9; 10] для формализации локальных деформаций. Цель настоящей статьи — синтезирование неевклидовой формы для тел с расширенной кинематикой.

В работе используются основные структуры и общие формулы, определенные в [1]. В частности, E есть евклидово физическое пространство [1, формула (2.1)] с ассоциированным векторным пространством V. Ортонормированный базис в последнем обозначается через cii=13. Символ SR обозначает промежуточную форму [1, формула (2.8)], SR — ее подлежащее многообразие[2] и т. д. В случае необходимости приводятся ссылки на соответствующие формулы из первой части работы.

2. Метод синтеза неевклидовой формы, используемый в работе, аналогичен методу подвижного репера, предложенному Картаном для построения неевклидовых пространств [11]. Остановимся более подробно на этой аналогии.

В случае простого материала переход к натуральному состоянию определяется полем локальных деформаций H:  SREndV [1, формулы (3.8) и~(3.11)], условие совместности которых имеет вид curl H=0. Это означает, что в случае выполнения последнего условия (и односвязности формы SR) поле H является градиентом некоторой глобальной деформации из промежуточной формы в глобально натуральную форму. Если же локальные деформации несовместны, то равенство нулю ротора не выполняется. Скажем тогда, что имеется источник несовместности, представленный тензорным полем второго ранга η. В этом случае уравнение, характеризующее структурную неоднородность тела, имеет вид

curl H=η. (1.1)

Подход, предложенный Кренером, заключается в преобразовании условия совместности с тем, чтобы результат этого преобразования можно было интерпретировать геометрически. Действительно, в соответствии с [1, 20 и 22] условие, выражающее равенство нулю ротора H, эквивалентно равенству нулю кручения T специфической связности Γ на многообразии SR, коэффициенты которой определяются формулами [1, (3.18)]. Следовательно, совместность деформаций эквивалентна утверждению, что связность Γ евклидова. Но тогда появление источника несовместности η равносильно изменению геометрии на многообразии SR: из евклидовой она переходит в неевклидову, представленную тензором кручения T0. В таком случае уравнение (1.1) преобразуется в равенство T=T0, характеризующее геометрическую структуру неевклидова пространства над SR.

С другой стороны, Картаном были получены структурные уравнения, характеризующие в общих чертах геометрию неевклидовых пространств [11–13]. Новизна идеи Картана заключалась в отказе от использования криволинейных координат и переходе к полям базисов более общего вида. Действительно, криволинейные координаты накладывают жесткие ограничения на виды локальных базисов за счет замены переменных (три функции определяют все поле базисов в трехмерном пространстве). Если же, отказываясь от использования замены переменных, перейти к неголономным базисам, то появляются дополнительные функциональные степени свободы, распоряжаясь которыми можно прийти к разнообразным геометриям. В явном виде поле базисов zii=13 — подвижный репер — может быть определено по заданному полю Ω: EAutV обратимых линейных преобразований[3] трансляционного пространства V в соответствии с равенствами zi=Ωci, i=1,2,3. Если через ϑii=13 обозначить репер пространства V, сопряженный к zii=13, а через gij=gzi, zj — метрические коэффициенты относительно подвижного репера, то структурные уравнения Картана для евклидова пространства имеют вид

gmkωjm+gmjωkmdgjk=0,dϑi+ωjiϑj=0,dωji+ωkiωjk=0 (1.2)

где — операция внешнего произведения [14]. Поля 1-форм ωji, i,j=1,2,3, относительно которых записаны соотношения (2), определяют взаимные искажения элементов репера zii3 при переходе от точки к точке.

Уравнения (1.2) можно рассматривать как «условия совместности» для заданных искажений подвижного репера. Действительно, в евклидовом пространстве реализуется не произвольный набор полей ωji, а лишь тот, который удовлетворяет равенствам (1.2). Переход к пространствам более общего вида можно осуществить, определив семейства 1-форм неметричности Qij, а также 2-форм кручения Ti и кривизны Rji, и подставив их в правые части уравнений (1.2):

gmkωjm+gmjωkmdgjk=Qjk,dϑi+ωjiϑj=Ti,dωji+ωkiωjk=Rji. (1.3)

Уравнения (1.3) определяют структуру пространства произвольной аффинной связности. В этой связи можно ассоциировать поля Qij, Ti и Rji с величинами, выражающими несовместность полей ωji с евклидовой геометрией.

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

2. Синтезирование неевклидовой отсчетной формы для среды второго градиента

2.1 Второй градиент деформации

3°. Квадратичное приближение деформации. Синтезированию поля локальных деформаций предпошлем рассмотрение кинематики среды второго градиента. Предположим, что упругий потенциал относительно промежуточной формы SR является отображением

SR×EndV×HomV; EndVX, F1, F2W^2X, F1, F2. (2.1)

Поэтому, если γDeformSR;S — произвольная деформация, то отклик тела в точке XSR характеризуется равенством [15; 16]

W=W^2X, F1X, F2X

В списке аргументов, наряду с первым градиентом деформации F1XEndV, представлен второй градиент F2XHomV; EndV, определяемый равенством F2X:=DXF1. Следовательно, в окрестности точки X выполнено разложение по формуле Тейлора второго порядка [17]:

γX+h=γX+F1Xh+12F2Xh, h+oh2, (2.2)

в котором hV — достаточно малый вектор, т. е. X+hSR. Здесь и в дальнейшем подразумевается отождествление второго градиента с билинейным отображением V×VV в силу естественного изоморфизма[4] [17] HomV; EndVvecL2V,V;V.

Обобщенная теорема Шварца [17] влечет, что второй градиент симметричен, т. е. F2Xv, u=F2Xu, v. Поэтому его можно восстановить по значениям F2Xv, v, vV в соответствии с формулой, известной из теории билинейных отображений [18]:

F2Xu, v=F2Xu+v, u+vF2Xu, uF2Xv, v2.

В свою очередь, из (2.2) вытекает следующее равенство:

F2Xv, v=2lim       s0vecγX, γX+svsF1Xvs2,

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

4°. Представления второго градиента. В прямоугольных координатах xii=13 деформации γ соответствует представление [1, формула (2.10)], а первому градиенту деформации F1X — разложение [1, формула (3.3)]. Поскольку второй градиент деформации можно записать как

F2X=XJF1XcJ,

то его разложение имеет вид[5]

F2X=2xiXIXJCoorXcicIcJ.

Вместе с тем несмотря на то что в криволинейных координатах QII=13 и qii=13 первый градиент F1 имеет простое разложение [1, формула (3.4)], представление, соответствующее второму градиенту, является более сложным. Для его получения рассмотрим следующее утверждение. Пусть

SRXAX=AJiXeiγXEJXEndV

— тензорное поле, представленное в паре координат QII=13 и qii=13. Тогда для его градиента

DXA=QKAXEKX

в точке XSR справедливо разложение[6]

DXA={KAJiσR+AJlFKjΓqjliγALiΓQKJL}XeiγXEJXEKX (2.3)

Здесь FJi — компоненты градиента деформации [1, формула (3.4)], а ΓQJKI и Γqjki — символы Кристоффеля, отвечающие криволинейным координатам QII=13 и qii=13, т. е.

ΓQJKI=EIJEK и Γqjki=eijek.

Доказательство. Используя обобщенное правило дифференцирования произведения [17], получаем равенство[7]

DXA=KAJIσRXeiγXEJX+AJiKeiγσRXEJX+AJieiγXKEJσRXEKX(2.4)

В нем, согласно определению символов Кристоффеля,

KeiγσRX=FKjXjeiσγX=FKjXΓqjilγXelγX,

KEJσRX=ΓQKLJXELX.

Подставляя полученные выражения в (2.4) и заменяя соответствующим образом индексы суммирования, приходим к формуле (2.3).

В случае, когда A=F1, формула (2.3) приводит к равенству

F2X={KFJiσR·+FJlFKjΓqjliγFLiΓQKJL}XeiγXEJXEKX (2.5)

являющемуся искомым разложением второго градиента. В нем

KFJi=2qiQJQK.

Следовательно, в силу перестановочности повторных производных, KFJi=JFKi.

В работе используется следующий частный вид формулы (2.5). Выберем в качестве qii=13 прямоугольные координаты xii=13. Тогда Γxjki=0, что дает

F2X={KFJiσRFLiΓQKJL}XciEJXEKX (2.6)

В этом случае компоненты второго градиента выражаются через компоненты первого градиента и их производные, а также через символы Кристоффеля системы координат QII=13.

Замечание 1. Евклидова структура физического пространства и отождествления по изоморфизму, индуцируемые ею, позволяют скрыть истинную природу полей, используемых в настоящей работе. Действительно, предположим заданными криволинейные координаты QII=13 на форме SR, а координаты на образах выберем прямоугольными. Тогда

  1. Деформация γ: SRS есть тройка γ1,γ2,γ3 скалярных полей γi:SR, определенных равенством

    γQ1,Q2,Q3=o+ciγiQ1,Q2,Q3

  1. Градиент деформации F1=ciFi есть тройка 1-форм F1,F2,F3, заданных как

    Fi=xiQIdQI,i=1,2,3.

    То, что градиент деформации не является тензором в классическом понимании этого термина, отмечалось в работе [19, с. 245] и позднее в монографиях [20, 21].

  1. Второй градиент F2=ciFi есть совокупность трех тензоров второго ранга F1, F2, F3.

В отличие от представления полей γ, F1, F2 в исходной (евклидовой форме), их представление в соответствии с п. 1)–3) является общим и может быть непосредственно перенесено на произвольные гладкие многообразия.

Несмотря на то что в рамках настоящей работы такое описание полей избыточно[8], к нему придется прибегнуть уже в случае рассмотрения деформирования материальных двумерных поверхностей в евклидовом пространстве.

2.2 Семейство форм и условие совместности

5°. Гипотеза локальной разгрузки. Подобно случаю среды первого градиента [1, 16°], предположим справедливой гипотезу локальной разгрузки. Пусть фиксированы тензоры N1EndV и N2HomV; EndV, характеризующие натуральное состояние, и пусть задана промежуточная форма SR вместе с упругим потенциалом (2.1). Предположим далее, что определено семейство γXXSR деформаций γX:SRSX, для которого в любой точке XSR выполнено условие[9]

W^2X, F1, F2F1F1=F1XX, F2=F2XX=N1, (2.7)

W^2X, F1, F2F2F1=F1XX, F2=F2XX=N2,

являющееся расширением [1, формула (3.6)]. Здесь F1X=DγX — первый градиент деформации γX, а F2X=DF1X, соответственно, второй градиент.

6°. Синтезирование локальных деформаций и гипердеформаций. Согласно второй из формул [1, (3.5)], в которой положим qii=13=xii=13, первый градиент деформации γX в точке YSR имеет представление

F1XY=F1XIiYciEIY. (2.8)

Из него, согласно равенству [1, (3.8)], синтезируется поле локальных деформаций , с разложением

H1X=HIiσRXciEIX,гдеHIiσRX=F1XIiX. (2.9)

Вместе с тем, поскольку рассматривается среда второго градиента, одного лишь поля H1 недостаточно для описания локальной разгрузки. Нужно еще поле, представляющее второй градиент. С этой целью определим вторые градиенты деформаций γX:

F2XY=DYF1X,

которые, в соответствии с представлением (2.8) и формулой (2.6), имеют разложение

F2XY={K F1XJiσR·F1XLiΓQKJL}YciEJYEKY (2.10)

Синтезируем теперь новый тензор H2XHomV; EndV по формуле

H2X:=F2XYY=X, (2.11)

и назовем его локальной гипердеформацией. Совокупность таких тензоров образует тензорное поле H2: SRHomV; EndV, которое предположим гладким.

Полагая теперь в формуле (2.10) Y=X, приходим к следующему представлению тензора гипердеформаций (2.11):

H2X={SKJiσRHLiσRΓQKJL}XciEJXEKX, (2.12)

где

SKJiσRX:=KF1XJiσRX, (2.13)

а HIi — компоненты разложения (2.9). Заметим, что в силу перестановочности повторных производных, числа SKJiσRX обладают свойством симметрии: SKJiσRX=SJKiσRX.

Таким образом, локальная разгрузка среды второго градиента характеризуется парой тензорных полей H1, H2, значения которых определяются согласно разложениям (2.9) и (2.12). При фиксированных криволинейных координатах QII=13 им соответствует поле матриц XHIiσRX,SJKiσRX,

где

detHIiσRX0,

SKJiσRX=SJKiσRX, XSR.

Кроме того, в соответствии с (2.7), выполнено свойство:

XSR:W^2X, F1, F2F1F1=H1X, F2=H2XN1, (2.14)

XSR:W^2X, F1, F2F2F1=H1X, F2=H2XN2.

7°. Восстановление семейства деформаций. Подобно случаю первого градиента, рассмотрим задачу восстановления семейства деформаций γXXSR по заданной паре гладких тензорных полей H1: SREndV и H2: SRHomV; EndV, первое из которых имеет в качестве значений обратимые линейные преобразования, а второе — симметричные тензоры. Кроме того, предполагается, что эти поля удовлетворяют свойству (2.14).

Зафиксировав точку XSR, изменим, если нужно, криволинейные координаты QII=13 так, чтобы точке X отвечали их нулевые значения[10]. Пусть HIiσRX и SJKiσRX=SKJiσRX — соответствующие функции из разложений (2.9) и (2.12). Тогда относительно пары координат QII=13 и xii=13 определим отображение 33согласно правилу:

xiQ1,Q2,Q3:=bXi+HIiσRXQI+SJKiσRXQJQK,i=1,2,3, (2.15)

где bXii=13 — фиксированная тройка чисел.

Будучи заданным как квадратичная форма, построенное отображение является гладким. Кроме того, в силу обратимости матрицы HIiσRX, теорема об обратной функции [17] гарантирует существование окрестности нуля, в которой (2.15) является диффеоморфизмом. Поэтому, если через NXSR обозначить соответствующую окрестность точки X в E, то придем к деформации γX:NXN~X формы NX, являющейся частью формы SR, в некоторую другую форму N~X. Отображение (2.15) является координатным представлением γX; кроме того, по построению

xiQIσRX=HIiσRX,2xiQJQKσRX=SJKiσRX,

что влечет равенства F1XX=H1X и F2XX=H2X.

Повторяя проделанную процедуру для всех точек формы SR, приходим к семейству деформаций γXXSR, по которому синтезируется пара H1, H2. Кроме того, для построенного семейства выполняется свойство (2.7), т. е. деформации γX являются разгрузочными. Вместе с тем в отличие от изначально определенных разгрузочных деформаций, заданных на всей промежуточной форме SR, полученные деформации определены лишь на ее частях. Это обстоятельство не противоречит общей методологии настоящей работы, поскольку рассматриваемые части конечны и потому также являются формами.

8°. Совместность локальных деформаций. В случае среды второго градиента будем называть пару H1, H2 тензорных полей, обладающую свойством (2.14), совместной, если 1) существует деформация γ0DeformSR;S0 из промежуточной формы SR в некоторую форму S0, для которой выполнено равенство H1=Dγ0, и 2) справедливо соотношение H2=DH1. Таким образом, для γ0 выполняется свойство

W^2X, F1, F2F1F1=DXγ0, F2=DX2γ0N1,

W^2X, F1, F2F2F1=DXγ0, F2=DX2γ0N2,

вытекающее из (2.7). В этом смысле форма S0 является глобальной натуральной.

Заметим, что если пара H1, H2 совместна, то, в частности, совместно поле локальных деформаций H1. При этом потенциал [1, формула (3.1)], относительно которого определяется натуральное состояние, индуцирован из W^2:

W^1X, F:=W^2X, F, H2X

Вместе с тем из совместности H1 в общем случае не вытекает совместность пары H1, H2. По этой причине уместно расширить терминологию:

  1. пара H1, H2 имеет первый порядок совместности, если для некоторой деформации γ0 выполнено свойство 1), однако свойство 2) не выполнено,
  2. пара H1, H2 имеет второй порядок совместности (или просто совместна), если выполнены оба свойства 1) и 2),
  3. пара H1, H2 несовместна, если ни одно из свойств 1), 2) не выполнено.

Таким образом, условие совместности [1, формула (3.15)] для случая простого материала дополняется вторым условием, представленным свойством 2), что приводит к совокупности равенств:

curlH1=0,DH1=H2. (2.16)

Их выполнение является необходимым (а в случае односвязности SR и достаточным) для первого и второго порядка совместности пары H1, H2. В координатной форме, согласно[11] (2.9) и (2.12),

JHKiKHJi=0JHKi=SJKi,i,J,K=1,2,3.

Вместе с тем первое из условий (2.16) избыточно, если рассматривать второй порядок совместности. Действительно, симметрия функций SJKi по нижним индексам дает

KHJi=SKJi=SJKi=JHKi.

Но это и означает, что curl H1=0. Поэтому второму порядку совместности отвечает следующее условие:

DH1=H2 или JHKi=SJKi,i,J,K=1,2,3. (2.17)

2.3. Геометрическая интерпретация условия совместности

9°. Поле Λ. Условие совместности (2.17) представим в следующем виде. Принимая во внимание, что значения поля H1 являются обратимыми линейными преобразованиями, домножим обе части (2.17) на H11iI и просуммируем по i:

H11iIJHKiH11iISJKi=0,I,J,K=1,2,3. (2.18)

Определим скалярные поля ΓJKI в соответствии с [1, формула (3.18)] и новые поля ΛJKI по формуле

ΛJKI:=H11iISJKi. (2.19)

Тогда равенство (18) принимает вид

ΓJKIΛJKI=0,I,J,K=1,2,3. (2.20)

В первой части работы [1, предложение 1] уже было показано, что поля ΓJKI являются коэффициентами некоторой связности. Установим, что аналогичное свойство выполняется и для полей ΛJKI.

10°. Свойства поля Λ. Доказательство того, что поля ΛJKI являются коэффициентами некоторой связности, предварим определением закона преобразования полей SJKi. Для этого продифференцируем обе части соотношения [1, формула (3.20)] по новым координатам Q~L. Если через ~L обозначить это дифференцирование, то полученная формула примет вид

~LFX~​ JiY=2QKQ~JQ~Lσ~RYFXKiY+QKQ~Jσ~RYQMQ~Lσ~RYMFXKiY. (2.21)

Полагая в выражении (2.21) Y=X и принимая во внимание определение (2.13), получаем искомый закон преобразования полей SJKi:

S~LJiσ~RX=SMKiσ~RXQMQ~Lσ~RXQKQ~Jσ~RX+HKiσ~RX2QKQ~LQ~Jσ~RX. (2.22)

Обозначая через SR многообразие, над которым определена промежуточная форма SR, приходим к следующему утверждению:

Предложение 1. Скалярные функции ΛJKI являются коэффициентами некоторой аффинной связности на многообразии SR.

Доказательство. Докажем утверждение методом, аналогичным тому, как было доказано предложение 1 в [1]. Если через Λ~JKI обозначить функции (2.19), определенные относительно координат Q~II=13, то в соответствии с формулами [1, формула (3.21)] и (2.22) получаем

Λ~JKI=H1~​ iIS~JKi=

=Q~IQLH1iLSPRiQPQ~JQRQ~K+HMi2QMQ~JQ~K=

=ΛPRLQ~IQLQPQ~JQRQ~K+Q~IQL2QLQ~JQ~K.

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

Связность, соответствующая полям ΛJKI, обозначается через Λ. Из симметрии полей SJKI по нижним индексам следует, что и полученная связность симметрична, т. е. ее кручение равно нулю: TΛ=0.

Получим другие выражения для функций ΛJKI. С этой целью восстановим по паре H1, H2 синтезирующее ее семейство деформаций γXXSR. В свою очередь, семейству γXXSR отвечает семейство первых градиентов F1XXSR и обратных градиентов F11XXSR с разложениями

F1X=F1XIiciEIиF11X=F11XiIEIci

соответственно. Их компоненты связаны соотношениями

F1XJiF11XjJδjiиF11XiIF1XJiδJI.(2.23)

Дифференцируя обратный градиент F11X в точке Y, принадлежащей области определения деформации γX, получаем равенство

DYF11(X)={KF11XjIσR+F11XJLΓQKLI}YEIYcjEKY. (2.24)

Заметим, что рассуждения, приводящие к последнему равенству, аналогичны использованным при выводе формулы (2.6). Полагая теперь Y=X, приходим к тензору третьего ранга

P2X:=DYF11(X)Y=X,

который по аналогии с имплантом [1, формула (3.12)] назовем гиперимплантом. Из (2.24) следует разложение для гиперимпланта:

P2X={RKjIσR+PjLσRΓQKLI}XEIXcjEKX

Здесь

RKjIσRX:=KF11XjIσRX,(2.25)

а PjL — компоненты поля имплантов P1=H11.

Для того чтобы связать поля RKjI и SJKI, рассмотрим второе из равенств (2.23). Дифференцируя обе его части по QK, получаем

F1XJiKF11XiI+F11XiIKF1XJi=0

что при Y=X дает

HJiRKiI+H1iISKJi=0

Из последнего равенства и вытекает, что

SKJi=HIiHJjRKjI. (2.26)

Следовательно, наряду с парой H1, H2, состоящей из полей локальных деформаций и гипердеформаций, можно рассматривать пару P1, P2, состоящую из импланта и гиперимпланта. В терминах элементов этой пары связность Γ определяется формулой [1, (3.26)], а связность Λ, в силу равенства (2.26), — формулой

ΛJKI=P11KjRJjI. (2.27)

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

F11XjIQK=F1XKiF11XjIxi

в произвольной точке Y. Переходя теперь к Y=X и учитывая (2.25), придем к равенству

RKjI=HKiR~ijI, (2.28)

где

R~ijI=F11XjIxiY=X,R~ijI=R~jiI.

Наконец, учитывая равенство (2.28), формулу (2.27) можно преобразовать к виду:

ΛJKI=P11JiP11KjR~ijI.

В таком виде связность Λ определяется в монографии [6, формула (2.59)].

11°. Тензор неоднородности. Связности Γ и Λ определяют новое поле

D:=ΓΛ (2.29)

с компонентами DJKI=ΓJKIΛJKI. Поле (2.29) является тензорным полем третьего ранга, в чем можно убедиться, сопоставив законы преобразования коэффициентов связностей Γ и Λ: слагаемые, соответствующие вторым частным производным, при вычитании взаимно уничтожаются. Следуя терминологии, предложенной в работе [6], назовем поле D тензором неоднородности.

В общем случае тензор неоднородности D не является симметричным, а его антисимметричная часть

DJKI=12DJKIDKJI (2.30)

пропорциональна кручению связности Γ. Действительно,

DJKI=ΓJKIΓKJI+ΓKJIΛKJI=TΓJKi+DKJi,

что влечет

2DJKI=TΓJKI.

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

D=12TΓ+Sym D, (2.31)

где Sym D — соответствующая симметричная часть.

12°. Условие совместности в терминах тензора неоднородности. Согласно определению (2.29) тензора неоднородности, условие совместности (2.20) может быть записано в окончательном виде:

DΓ,Λ=0 или DJKI=0I,J,K=1,2,3. (2.32)

Это необходимое (а в случае односвязности SR и достаточное) условие совместности пары H1, H2.

Заметим, что в случае выполнения условия (2.32) из равенства (2.31) следует, что TΓ=0, т. е. поле локальных деформаций H1 совместно. Этот результат находится в полном соответствии с установленном выше свойством, согласно которому совместность второго порядка влечет совместность первого порядка. Таким образом, в терминах полей кручения и неоднородности предложенную классификацию совместности деформаций можно записать следующим образом:

  1. пара H1, H2 имеет первый порядок совместности, если TΓ=0, но D0,
  2. пара H1, H2 имеет второй порядок совместности, если D=0,
  3. пара H1, H2 несовместна, если TΓ0 и D0.

2.4. Неевклидова отсчетная форма

13°. Материальная связность. Связности Вайценбока Γ [1, формула (3.18)], Леви-Чивита L [1, формула (3.32)] и Λ (2.19) могут служить примерами материальных связностей на SR. Каждой из них отвечает своя неевклидова отсчетная форма. Вместе с тем можно рассматривать и комбинации этих трех связностей, образуя разнообразные поля. Действительно, например, подберем такую связность Δ, чтобы были выполнены следующие условия:

  1. кручение Δ совпадает с кручением Γ,
  2. кривизна Δ отлична от нуля,
  3. в случае совместных деформаций связность Δ является евклидовой связностью.

Простейшим вариантом является связность, определяемая формулой

Δ:=L+D, или в компонентах ΔJKI:=LJKI+DJKI. (2.33)

Здесь D — тензор неоднородности (2.29). В явном виде, принимая во внимание, что компоненты материальной метрики G определены равенством GIJ=δijHIiHJj, а также используя формулы [1, (3.18)] и (2.19), формулу (2.33) можно записать как

ΔJKI=GIS2JGSK+KGJSSGJK+H11iIJHKiSJKi. (2.34)

Вклад второго градиента в геометрическую структуру представлен полями SJKi, которые наряду с HKi определяют локальную гипердеформацию H2.

Поскольку первое слагаемое в (2.33) является связностью, а второе — тензором третьего ранга, то, составляя закон преобразования для полей ΔJKI, можно убедиться в том, что они действительно являются коэффициентами связности на SR. Полученная связность удовлетворяет требованиям а)– в), указанным выше. В самом деле, поскольку связности L и Λ симметричны, то

TΔ=TΓ

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

ΔABCI=LABCI+DABCI+LBCEDAEI+DBCELAEILACEDBEIDACELBEI.

Здесь L — кривизна, порожденная тензором Леви-Чивита, а D — формальное выражение, которое представляет «кривизну», порожденную тензором D. Наконец, если деформации совместны, то D=0 и потому связность Δ совпадает со связностью Леви-Чивита. Но связность Леви-Чивита в этом случае сама порождается совместными деформациями, и потому ее кривизна равна нулю. Следовательно, в случае совместных деформаций Δ является евклидовой связностью.

Связность Δ можно в определенном смысле рассматривать как модификацию связности Леви-Чивита L для случая второго градиента. Аналогичная модификация связности Вайценбока Γ приводит к полю

Δ:=~Γ+D, или в компонентах Δ~JKI:=ΓJKI+DJKI. (2.35)

В явном виде поле (2.35) представлено выражением

Δ~JKI=H11iI2JHKiSJKi. (2.36)

Рассуждениями, аналогичными тем, что касались поля Δ, можно установить, что поле Δ~также является связностью на SR. Эта связность удовлетворяет свойствам б) и в), и с точностью до множителя удовлетворяет свойству а):

TΔ~=2TΓ

В рамках физических приложений геометрия, определяемая связностью (2.33) (или (2.35)), могла бы быть полезна в случае одновременного наличия дислокаций, дисклинаций и точечных дефектов[12] в кристалле. Отметим, что, по-видимому, первое указание на то, что в рамках второго градиента можно эффективно смоделировать наличие дислокаций и дисклинаций, содержится в работе [22].

14°. Структура неевклидовой формы. Связность Δ, определенная согласно формуле (2.33), материальная метрика G, заданная формулой [1, формула (3.28)], и соответствующая метрике форма объема μ=dVG задают на многообразии SR структуру пространства с неевклидовой связностью:

SR=SR, G,dVG,Δ (2.37)

Пространство SR является искомой неевклидовой отсчетной формой над материальным многообразием, соответствующей теории второго градиента. Его геометрия характеризуется полями кручения TΓ, кривизны Δ и неметричности QΔ. Подобно тому как это было сделано в работе [10], можно показать, что эти поля не зависят от выбора промежуточной формы SR.

Заменяя материальную связность (2.33) на связность (2.35), от неевклидовой формы (2.37) можно перейти к неевклидовой форме

S~R=SR, G,dVG,Δ~ (2.38)

Форму S~R также следует считать искомой, поскольку в рамках построений настоящей статьи нет никакого способа предпочесть форму (2.37) форме (2.38). Следовательно, здесь возникает та же ситуация, что и для простого материала, где в качестве материальной можно выбирать как связность Вайценбока, так и связность Леви-Чивита.

2.5. Пример синтезирования неевклидовой формы

15°. Семейство деформаций. Проиллюстрируем рассуждения на примере центрально-симметричного деформирования промежуточной формы

SR=XE:  Ri<Xo<Re, (2.39)

развивая соответствующий пример из [1]. Именно предположим, что определено семейство γρρ]Ri,Re[ деформаций γρ:SRSρ, для которого выполнены условия:

а) при каждом значении параметра ρ]Ri,Re[ деформация центрально-симметрична. Кроме того, для всех точек сферы Lρ=XE:  Xo=ρ справедливо свойство:

XLρ:W^2X, F1, F2F1F1=DXγρ, F2=DX2γρ=N1,

XLρ:W^2X, F1, F2F2F1=DXγρ, F2=DX2γρ=N2,

б) в сферических координатах r,θ,φ каждая деформация γρ имеет представление [1, (3.35)].

16°. Синтезирование поля локальных гипердеформаций. Используя матрицу компонент первого градиента F1ρ [1, формула (3.37)], получим гиперматрицу KF1ρJi частных производных от этих компонент. Ей соответствуют следующие три 3×3-матрицы:

KF1ρJ1=ωρf0''RsinΘcosΦωρf0'RcosΘcosΦωρf0'RsinΘsinΦωρf0'RcosΘcosΦωρf0RsinθcosΦωρf0RcosΘsinΦωρf0'RsinΘsinΦωρf0RcosΘsinΦωρf0RsinΘcosΦ,

KF1ρJ2=ωρf0''RsinΘsinΦωρf0'RcosΘsinΦωρf0'RsinΘcosΦωρf0'RcosΘsinΦωρf0RsinθsinΦωρf0RcosΘcosΦωρf0'RsinΘcosΦωρf0RcosΘcosΦωρf0RsinΘsinΦ,

KF1ρJ3=ωρf0''RcosΘωρf0'RsinΘ0ωρf0'RsinΘωρf0RcosΘ0000.

Полагая теперь, в соответствии с (2.13), R=ρ и заменяя затем все вхождения ρ на R, приходим к матрицам:

SJK1=ωRf0''RsinΘcosΦωRf0'RcosΘcosΦωRf0'RsinΘsinΦωRf0'RcosΘcosΦωRf0RsinθcosΦωRf0RcosΘsinΦωRf0'RsinΘsinΦωRf0RcosΘsinΦωRf0RsinΘcosΦ,SJK2=ωRf0''RsinΘsinΦωRf0'RcosΘsinΦωRf0'RsinΘcosΦωRf0'RcosΘsinΦωRf0RsinθsinΦωRf0RcosΘcosΦωRf0'RsinΘcosΦωRf0RcosΘcosΦωRf0RsinΘsinΦ,SJK3=ωRf0''RcosΘωRf0'RsinΘ0ωRf0'RsinΘωRf0RcosΘ0000,(2.40)

которыми, наряду с матрицей [1, (3.38)], определяется поле локальных гипердеформаций H2.

17°. Синтезирование геометрии. Метрический тензор G, связность Леви-Чивита L и связность Вайценбока Γ синтезированы в рамках примера из статьи [1] и представлены формулами [1, (3.39)], [1, (3.40)] и [1, (3.42)] соответственно. В свою очередь, связность Λ определяется по (2.19) и в силу формул [1, (3.38)] и (2.40) ее отличные от нуля коэффициенты имеют вид:

Λ111=f0''f0',Λ221=f0f0',Λ331=f0sin2Θf0',Λ122=Λ212=Λ133=Λ313=f0'f0,Λ332=sinΘcosΘ,Λ233=Λ323=cotΘ.

Полученным выражениям для коэффициентов связностей Γ и Λ отвечает тензор неоднородности с компонентами

D111=D122=D133=ω'ω. (2.41)

Следовательно, подобно случаю первого градиента, условие совместности поля локальных деформаций представлено равенством ω'=0.

Подстановка соотношений [1, (3.40)] и (2.41) в общую формулу (2.33) приводит к следующим выражениям для коэффициентов связности Δ:

Δ111=f0''f0'+2ω'ω,Δ221=f0ωf0'+f0ω'ωf0'2,Δ331=f0ωf0'+f0ω'sin2Θωf0'2,Δ122=Δ133=f0'f0+2ω'ω,Δ212=Δ313=f0'f0+ω'ω,Δ332=sinΘcosΘ,Δ233=Δ323=cotΘ. (2.42)

Кручение связности Δ совпадает с кручением связности Γ и представлено выражениями [1, формула (3.43)]. Помимо этого, в рамках рассматриваемой модельной задачи кривизна связности Δ совпадает с кривизной связности Леви-Чивита L. В частности, скалярная кривизна представлена выражением [1, формула (3.41)]. Наконец, тензор неметричности связности Δ [1, формула (3.34)] характеризуется следующими компонентами:

Q111=ωω'f0'2,Q122=ωω'f02,Q133=ωω'f02sin2Θ.

Соотношения [1, (3.39)] и (2.42) определяют неевклидову форму (2.37) частного вида. Мера несовместности деформаций, представленная тензором D, одновременно характеризует отклонение геометрии от евклидовой.

3. Синтезирование неевклидовой отсчетной формы для среды с микроструктурой

3.1. Микроморфный континуум

18°. Деформации. В заключение рассмотрим микроморфный континуум — среду, точки которой, наряду с трансляционными степенями свободы, имеют дополнительные степени свободы, т. е. обладают внутренней структурой. Деформация континуума представлена парой γ, M, где γDeformSR;S — макродеформация, соответствующая деформации в рамках теории простого материала, а M: SREndV — микродеформация, являющаяся полем обратимых линейных отображений, характеризующим изменение внутренней структуры частиц среды. Предположим, что упругий потенциал относительно промежуточной формы SR представлен отображением [23]

SR×EndV×EndV×HomV; EndVX, F, M1, M2W^3X, F, M1, M2. (3.1)

Тогда отклик тела в точке XSR характеризуется равенством

W=W^3X, FX, MXDXM

где FX=DXγ — градиент макродеформации γ в точке X. Отметим, что в частном случае, если положить M=F, то придем к среде второго градиента.

19°. Координатное представление деформаций. В паре криволинейных координат QII=13 и qii=13, ассоциированных с формами SR и S соответственно, градиент макродеформации F представлен разложением [1, (3.4)]. Аналогичным образом в виде двухточечного разложения представим микродеформацию M:

MX= MIiXeiγXEIX.

Тогда, в соответствии с формулой (2.3), для градиента DXM справедливо разложение

DXM=KMJiσR+MJlFKjΓqjliγMLiΓQKJLXeiγXEJXEKX.

В дальнейшем в качестве координат qii=13 выбираются прямоугольные координаты xii=13.

Замечание 2. Тело с микроструктурой деформируется более сложным образом, чем простой материал. Как и в случае последнего, деформация тела с микроструктурой сопровождается изменением формы, т. е. области физического пространства, занимаемой телом. Это изменение определяется отображением γ. Однако (и в этом состоит отличие от простого материала) наряду с изменением формы меняются и положения элементов внутри микроструктуры [24], моделируемой касательным пространством к SR в рамках первого приближения. Именно поэтому, независимо от γ, задано двухточечное тензорное поле M, определяющее изменение положений элементов микроструктуры.

3.2. Семейство форм и условие совместности

20°. Гипотеза локальной разгрузки. В случае микроморфного континуума гипотеза локальной разгрузки принимается в следующем виде. Пусть фиксированы тензоры N1EndV и N2HomV; EndV, и пусть выбрана промежуточная форма SR вместе с упругим потенциалом (3.1). Далее предположим, что определено семейство пар γX, MXXSR, состоящих из макродеформаций γX:SRSX и микродеформаций MX: SREndV, такое, что в любой точке XSR выполнено условие

W^3X, F, M1, M2FF=FXX, M1=MXX, M2=DXMX=N1,W^3X, F, M1, M2M2F=FXX, M1=MXX, M2=DXMX=N2. (3.2)

Здесь FX=DγX — градиент макродеформации γX.

21°. Синтезирование локальных макродеформаций и микродеформаций. Опираясь на определение [1, (3.8)], синтезируем поле локальных макродеформаций H: SREndV. В каждой точке XSR локальная макродеформация HX имеет разложение [1, (3.11)], совпадающее с разложением локальной деформации в случае среды первого градиента.

По семейству MXXSR синтезируем новые тензорные поля. Прежде всего, положим

AX:=MXX,XSR. (3.3)

Совокупности обратимых тензоров AX отвечает поле A: SREndV, которое будем полагать гладким. В координатном представлении,

A=AIiciEI. (3.4)

Далее, синтезируем тензор третьего ранга

BX:=DYMXY=X,XSR

и назовем его локальной микродеформацией. Приходим к полю B: SRHomV; EndV, которое также предположим гладким. Отметим, что в отличие от второго градиента макродеформации γ, это поле в общем случае не является симметричным по нижним индексам. Формула (2.3) влечет следующее представление для BX, подобное (2.12):

BX=SKJiσRHLiσRΓQKJLXciEJXEKX, (3.5)

в котором, однако, величины

SKJiσRX:=KMXJiσRX (3.6)

уже несимметричны по нижним индексам.

Таким образом, локальная разгрузка микроморфной среды характеризуется тройкой тензорных полей H, A, B, определяемых разложениями [1, (3.11)], (3.4) и (3.5). При фиксированных криволинейных координатах QII=13 им соответствует поле матриц XHIiσRX,AIiσRX, SJKiσRX,

где

detHIiσRX0,

detAIiσRX0,

SKJiσRXSJKiσRX,XSR.

Кроме того, в соответствии с (3.2) выполнено свойство:

W^3X, F, M1, M2FF=HX, M1=AX, M2=BX=N1,W^3X, F, M1, M2M2F=HX, M1=AX, M2=BX=N2 (3.7)

во всех точках XSR.

22°. Совместность локальных деформаций. Назовем тройку H, A, B тензорных полей, обладающих свойством (3.7), совместной, если 1) существует деформация γ0DeformSR;S0 из промежуточной формы SR в некоторую форму S0, для которой выполнено равенство H=Dγ0, и 2) существует поле обратимых линейных преобразований M0: SREndV такое, что A=M0 и B=DM0. В этом случае для пары γ0, M0 выполняется свойство

W^3X, F, M1, M2FF=DXγ0, M1=M0X, M2=DXM0=N1,

W^3X, F, M1, M2M2F=DXγ0, M1=M0X, M2=DXM0=N2,

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

Как и в случае среды второго градиента, если тройка H, A, B совместна, то, в частности, совместно поле локальных деформаций H. Потенциал [1, (3.1)], относительно которого определяется натуральное состояние, индуцирован из W^3:

W^1X, F:=W^3X, F, AX, BX

Однако из совместности H в общем случае не вытекает совместность всей тройки H, A, B. Подобно среде второго градиента, можно соответствующим образом расширить терминологию.

На основании условий 1) и 2), условие совместности для среды с микроморфной кинематикой можно представить в виде совокупности равенств:

curl H=0​​​​​​​​,

DA=​​​​​​​​B.

В координатной форме, согласно [1, (3.11)], (3.4) и (3.5),

JHKiKHJi=0 (3.8)

JAKiSJKi=0,i,J,K=1,2,3. (3.9)

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

3.3. Геометрическая интерпретация условия совместности

23°. Поля Σ и Λ. Выполним следующее преобразование условия совместности (3.8)–(3.9). Поскольку матричные поля HIi и AIi обратимы, домножим обе части равенства (3.8) на H1iI, а обе части (3.9) — на A1iI, произведя в обоих случаях суммирование по i. Тем самым придем к следующей совокупности соотношений, эквивалентных исходным:

H1iIJHKiH1iIKHJi=0,A1iIJAKiA1iISJKi=0,i,J,K=1,2,3. (3.10)

Задавая поля ΓJKI по формуле [1, формула (3.18)] и вводя новые поля ΣJKI по формуле, аналогичной [1, формула (3.18)],

ΣJKI:=A1iIJAKi, (3.11)

и поля ΛJKI по аналогии с (2.19) как

ΛJKI:=A1iISJKi, (3.12)

можно переписать равенства (3.10) в лаконичном виде:

ΓJKIΓKJI=0ΣJKIΛJKI=0,i,J,K=1,2,3. (3.13)

24°. Кручение и тензор неоднородности. Каждое из полей ΓJKI, ΣJKI и ΛJKI задает связность на SR. Действительно, для первого поля это уже было установлено при рассмотрении среды первого градиента; второе поле рассматривается аналогично. В случае третьего поля можно дословно повторить доказательство предложения 1. Отметим, что, в отличие от своего аналога из теории второго градиента, поля (3.12) несимметричны.

Таким образом, равенства (3.13) можно записать как

TΓJKI=0,DJKI=0, i,J,K=1,2,3, (3.14)

где TΓ — кручение связности Γ, а

D:=ΣΛ (3.15)

— аналог тензора неоднородности. Тем самым показано, что в случае микроморфного континуума меры несовместности представлены независимыми полями кручения TΓ и неоднородности D.

Тензор неоднородности D не симметричен по нижним индексам, а его антисимметричная часть с компонентами (2.30) пропорциональна разности кручений связностей Σ и Λ. Действительно,

DJKI=ΣJKIΛJKI=

=ΣJKIΣKJI+ΣKJIΛJKI=

=TΣJKI+ΣKJIΛJKI=

=TΣJKI+ΣKJIΛKJI+ΛKJIΛJKI=

=TΣJKITΛJKI+DKJI,

что и приводит к равенству

2DJKI=TΣJKITΛJKI.

3.4. Неевклидова отсчетная форма

25°. Материальная связность. Связности Вайценбока Γ [1, формула (3.18)], Леви-Чивита L [1, формула (3.32)], и новые связности Σ (3.11), Λ (3.12) служат примерами материальных связностей на SR, каждой из которых отвечает своя неевклидова отсчетная форма. Наряду с этим, можно рассматривать и комбинации этих связностей. Примером может служить связность Δ, заданная формулой

Δ:=L+D, (3.16)

аналогичной (2.33). В явном виде, принимая во внимание формулу (3.15), а также определения (3.11) и (3.12), получаем выражение

ΔJKI=GIS2JGSK+KGJSSGJK+A1iIJAKiSJKi. (3.17)

Кручение связности Δ определяется равенством

TΔ=TΣTΛ

а кривизна в общем случае отлична от нуля.

Альтернативную связность на SR можно получить, если вместо связности Леви-Чивита L использовать связность Вайценбока Γ. Действительно, положим

Δ~:=Γ+D. (3.18)

По форме построенная связность напоминает связность (2.35), однако, если раскрыть определения входящих в нее полей, то выражение для ее коэффициентов имеет иной вид:

Δ~JKI=H11iIJHKi+A1iIJAKiSJKi. (3.19)

Кручение связности (3.18) определяется формулой

TΔ~=TΓ+TΣTΛ

Заметим, что если локальные деформации совместны (т. е. выполнены равенства (3.14)), то связность, определяемая формулой (3.16) (или (3.18)), является евклидовой. Поэтому она также может быть полезна для моделирования несовместных деформаций, возникающих в средах с дефектами.

Замечание 3. В частном случае, когда микродеформация полагается равной градиенту деформации, теория среды с микроструктурой сводится к теории среды второго градиента. Тогда формула (3.17) переходит в формулу (2.34), а формула (3.19) — в формулу (2.36).

26°. Структура неевклидовой формы. Связность Δ (3.16), материальная метрика G [1, формула (3.28)], и соответствующая метрике форма объема μ=dVG вводят на многообразии SR структуру пространства с неевклидовой связностью:

SR=SR, G,dVG,Δ (3.20)

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

Если вместо материальной связности (3.16) использовать связность (3.18), то придем к альтернативному представлению неевклидовой формы:

S~R=SR, G,dVG,Δ~ (3.21)

Как и форма (3.20), синтезированная форма (3.21) в равной степени пригодна для моделирования глобального натурального состояния тела с микроструктурой. Выбор определенной формы требует иных соображений, выходящих за рамки статьи.

3.5. Расширенное описание для тела с микроструктурой

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

Действительно, естественным формализмом, пригодным для построения расширения форм, является теория расслоений, восходящая к работам Фельдбау, Эресманна и Стинрода[13] [27–29]. В рамках этого формализма над отсчетной формой B и физическим пространством E выстраиваются тотальные пространства локально тривиальных расслоений с одним и тем же типовым слоем F. В явном виде эти расслоения определены структурами [30]

ξB=B,B,πB,F (3.22)

— материальное расслоение и

ξE=E,E,πE,F (3.23)

— физическое расслоение, где B и E — гладкие многообразия, представляющие тотальные пространства расслоений, а πB:BB и πE:EE — гладкие сюръективные отображения. Многообразие F — общий элемент структур ξB и ξE — представляет типовой слой расслоения.

Предполагается, что расслоения (3.22) и (3.23) локально тривиальны. Это означает (случай расслоения ξE идентичен рассмотренному далее), что для любой точки XB существуют окрестность U в B и диффеоморфизм Φ:πB1UU×F, называемый локальной тривиализацией, такие, что prUΦ=πB, где prU: U×FU — проекция декартова произведения на первый сомножитель. Таким образом, локально диффеоморфизм Φ должен иметь представление Φp=πBp,φp для некоторого гладкого отображения φ:πB1UF, свойства которого зависят от частного вида расслоения. Примеры материальных расслоений и их физическую интерпретацию можно найти в работе [31].

28°. Дополненные деформации. Тотальные пространства B и E расслоений (3.22) и (3.23) над отсчетной формой B и физическим пространством E определяют дополненную отсчетную форму и дополненное физическое пространство соответственно. Каждой форме SE соответствует дополненная форма SE, определяемая как прообраз S=πE1S.

Дополненная деформация определяется как отображение κ:BS, удовлетворяющее свойству χπB=πEκ, где χ:BS — классическая деформация. Таким образом, в специально подобранных локальных координатах на тотальных пространствах можно записать

κ~Q1,Q2,Q3f1,,fm=χ~Q1,Q2,Q3;τQ1,Q2,Q3;f1,,fm

Здесь QII=13 — локальные координаты на форме B, а faa=1m — локальные координаты на типовом слое F. Тильдами сверху обозначены координатные представления соответствующих отображений. Отображение τ:3+mm определяет тип микроструктуры.

В рамках описания микроструктуры, используемого в настоящей работе, в качестве расслоений (3.22) и (3.23) достаточно выбрать касательные расслоения. Несмотря на их тривиальность в случае трехмерных тел, они перестают быть тривиальными при рассмотрении материальных поверхностей. При таком выборе представление дополненной деформации есть не что иное, как пара [31] γ, M. Этим получен желаемый результат: формализм касательных расслоений позволяет интерпретировать пару γ, M как одно точечное отображение между касательными пространствами.

Использование формализма расслоений приводит к необходимости модифицировать геометрический язык для описания неевклидовой отсчетной формы. Вместо теории аффинных связностей применяется более общая теория связностей на расслоениях. Подробно этот вопрос рассмотрен в работе [32].

3.6. Пример синтезирования неевклидовой формы

29°. Семейство деформаций. Завершим основные построения настоящего раздела модельным примером. Предположим, что промежуточная форма SR является полым шаром (39) и, что на ней задано семейство γρ, Mρρ]Ri,Re[, состоящее из макродеформаций γρ:SRSρ и микродеформаций Mρ: SREndV. Это семейство подчиняется следующим условиям:

а) для любого ρ]Ri,Re[ макродеформация γρ является центрально-симметричной и имеет представление [1, (3.35)] в сферических координатах r,θ,φ;

б) для каждого ρ]Ri,Re[ микродеформация MρX в точках X сферы Lρ={XE:  Xo=ρ} является поворотом материальных элементов вокруг базисного вектора ERX на угол τρ;

в) для всех точек сферы Lρ справедливы соотношения:

XLρ:W^3X, F, M1, M2FF=DXγρ, M1=MρX, M2=DXMρ=N1,

XLρ:W^3X, F, M1, M2M2F=DXγρ, M1=MρX, M2=DXMρ=N2.

Явный вид тензора MρX, задающего поворот вокруг оси, определяемой единичным вектором k=ERX, на угол τρ, может быть определен в соответствии с формулой Родрига [33; 34]. Относительно диадного базиса cicj тензор MρX представлен матрицей

MρXdE+sinτρK+1cosτρK.K,

где E — единичная матрица, а

K0k3k2k30k1k2k10

— матрица, составленная из компонент вектора k=kici. Тогда матрица MρX относительно смешанного диадного базиса ciEI определяется по формуле

MρX=MρXd.T,

где T — матрица перехода от базиса cii=13 к базису ER, EΘ, EΦ в точке X.

Окончательно

Mρ]Ji=​​sinΘcosΦRcosΘcosΦcosτρsinΦsinτρRsinΘcosτρsinΦ+cosΘcosΦsinτρsinΘsinΦRcosΘcosτρsinΦ+cosΦsinτρRsinΘcosΦcosτρcosΘsinΦsinτρcosΘRcosτρsinΘRsin2Θsinτρ. (3.24)

30°. Синтезирование поля локальных деформаций. Поле H локальных макродеформаций определено в [1] и представлено формулой [1, (3.38)]. Осталось определить поля и B, соответствующие локальным микродеформациям. Первое из них синтезируется по формуле (3.3). Полагая в (3.24) R=ρ и заменяя затем все переменные ρ на R, получаем

AJi=sinΘcosΦRcosΘcosΦcosτRsinΦsinτRRsinΘcosτRsinΦ+cosΘcosΦsinτRsinΘsinΦRcosΘcosτRsinΦ+cosΦsinτRRsinΘcosΦcosτRcosΘsinΦsinτRcosΘRcosτRsinΘRsin2ΘsinτR. (3.25)

Для синтезирования поля требуется еще определить скалярные поля SJKi в соответствии с (3.6). С этой целью получим гиперматрицу KMρJi, дифференцируя элементы матрицы (3.24). В результате придем к следующим матрицам:

KMρJ1=0cosΘcosΦcosτρsinΦsinτρsinΘcosτρsinΦ+cosΘcosΦsinτρcosΘcosΦRcosΦcosτρsinΘRcosΘcosτρsinΦ+cos2ΘcosΦsinτρsinΘsinΦRcosΘcosτρsinΦ+cosΦsinτρRsinΘcosΘsinΦsinτρcosΦcosτρ,

KMρJ2=0cosΘcosτρsinΦ+cosΦsinτρsinΘcosΦcosτρcosΘsinΦsinτρcosΘsinΦRcosτρsinΘsinΦRcosΘcosΦcosτρcos2ΘsinΦsinτρcosΦsinΘRcosΘcosΦcosτρsinΦsinτρRsinΘcosτρsinΦ+cosΘcosΦsinτρ,

KMρJ3=0cosτρsinΘsin2ΘsinτρsinΘRcosΘcosτρRsin2Θsinτρ000.

Из них при ρ=R вытекают представления для искомых функций:

SKJ1=0cosΘcosΦcosτRsinΦsinτRsinΘcosτRsinΦ+cosΘcosΦsinτRcosΘcosΦRcosΦcosτRsinΘRcosΘcosτRsinΦ+cos2ΘcosΦsinτRsinΘsinΦRcosΘcosτRsinΦ+cosΦsinτRRsinΘcosΘsinΦsinτRcosΦcosτR,

SKJ2=0cosΘcosτRsinΦ+cosΦsinτRsinΘcosΦcosτRcosΘsinΦsinτRcosΘsinΦRcosτRsinΘsinΦRcosΘcosΦcosτRcos2ΘsinΦsinτRcosΦsinΘRcosΘcosΦcosτRsinΦsinτRRsinΘcosτRsinΦ+cosΘcosΦsinτR, (3.26)

SKJ3=0cosτRsinΘsin2ΘsinτRsinΘRcosΘcosτRRsin2ΘsinτR000.

31°. Синтезирование геометрии. Поскольку метрический тензор G, связность Леви-Чивита L и связность Вайценбока Γ определяются по полю локальных макродеформаций H, то их выражения известны и представлены формулами [1, (3.39), (3.40), (3.42)] соответственно. Новыми являются выражения для связностей Σ, Λ и результирующей связности Δ.

Коэффициенты связности Σ определяются в соответствии с общей формулой (3.11). Принимая во внимание выражение (3.25) для поля A, приходим к равенствам (приведены только те функции, которые отличны от нуля)

Σ221=RcosτR,Σ231=RsinΘsinτR,Σ321=RsinΘsinτR,Σ331=RcosτRsin2Θ,Σ122=Σ133=1R,Σ132=sinΘτ'R,Σ212=cosτRR,Σ312=sinΘsinτRR,Σ332=cosΘsinΘ,Σ123=cscΘτ'R,Σ213=cscΘsinτRR,Σ233=Σ323=cotΘ,Σ313=cosτRR. (3.27)

Из представленных выражений следует, что связность Σ несимметрична, т. е. обладает ненулевым кручением. Отметим, что последнее отлично от кручения связности Γ.

Связность Λ определяется согласно формуле (3.12). Учитывая соотношение (3.25) для поля и формулу (3.26), определяющую частный вид полей SJKi, получаем следующие соотношения:

Λ221=RcosτR,Λ231=RsinΘsinτR,Λ321=RsinΘsinτR,Λ331=RcosτRsin2Θ,Λ122=Λ133=1R,Λ212=cosτRR,Λ312=sinΘsinτRR,Λ332=cosΘsinΘ,Λ213=cscΘsinτRR,Λ233=Λ323=cotΘ,Λ313=cosτRR. (3.28)

Несмотря на то что большинство соответствующих коэффициентов связностей Σ и Λ совпадает, между ними все же есть отличные друг от друга. Например, Λ132=0, в то время как Σ132=sinΘτ'R. Мера отличия связностей друг от друга представлена тензором неоднородности (3.15), отличные от нуля компоненты которого, в соответствии с (3.27) и (3.28), имеют вид:

D132=sinΘτ'R,D123=cscΘτ'R (3.29)

Следовательно, с учетом выражений [1, (3.43)] для кручения связности Γ условия совместности (3.14) локальных деформаций в рассматриваемом модельном случае имеют вид

ω'=0 и τ'=0

Подстановка соотношений [1, (3.40)] для коэффициентов связности Леви-Чивита и равенств (3.29) в общую формулу (3.16) приводит к следующим выражениям для коэффициентов связности Δ:

Δ111=ω'ω+f0''f0',Δ221=f0ωf0'+f0ω'ωf0'2,Δ331=f0sin2Θωf0'+f0ω'ωf0'2,Δ122=Δ212=Δ133=Δ313=ω'ω+f0'f0,Δ132=sinΘτ'R,Δ332=cosΘsinΘ,Δ123=cscΘτ'R,Δ233=Δ323=cotΘ. (3.30)

Кручение связности Δ отлично от нуля и представлено следующими компонентами:

T132=T312=sinΘτ'R,T123=T213=cscΘτ'R

Кривизна связности Δ также отлична от нуля. В частности, скалярная кривизна имеет вид [1, (3.41)]. Неметричность равна нулю.

Формулы [1, (3.39)] и (3.30) определяют неевклидову форму (3.20). Меры несовместности деформаций, представленные полем кручения связности Γ и тензором неоднородности D, одновременно характеризуют отличие геометрии от евклидовой.

Заключение

Геометрический метод, изложенный в статье [1], обобщен в настоящей работе на тела с расширенной кинематикой. Показано, что благодаря дополнительным полям AIi и SJKi, определяемым спецификой расширенной кинематики, в таких телах имеется больше вариантов для синтезирования геометрии. Сводка результатов представлена в таблице.

 

Таблица. Различные материальные связности

Table. Various material connections

 

Γ

L

Σ

Λ

Δ

Δ~

Простой материал

[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)

 

Ее строки содержат ссылки на явные представления материальных связностей. Во всех трех случаях кинематики, рассмотренных в работе, связности Вайценбока Γ и Леви-Чивита L вычисляются по одним и тем же формулам. Связности Σ и Λ могут быть определены лишь для тел с расширенной кинематикой, поскольку вычисление их коэффициентов требует знания дополнительных полей AIi и SJKi. Наконец, связности Δ и Δ~ служат расширениями связностей Леви-Чивита и Вайценбока. В рамках физических приложений они могут быть полезны для геометрического моделирования тел с дислокациями и дисклинациями, поскольку их кручение и кривизна — плотности дислокаций и дисклинаций, — отличны от нуля.

Благодарности

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

 

1 См. также [1, 16° и 17°].

2 В рамках классической механики континуума нет необходимости выделять третью форму, наряду с отсчетной и актуальной формами. Вместе с тем в настоящей работе третья форма, называемая промежуточной, используется для синтезирования геометрии на материальном многообразии. Формально, подобно любой другой форме тела (см. [1, 7°]), промежуточная форма представлена в виде упорядоченной совокупности SR=SR,gSR,ϵSR,SR, где SR — подлежащее многообразие (носитель геометрической структуры), а gSR, ϵSR, SR обозначают метрику, форму объема и связность над многообразием SR, индуцированные из евклидова пространства E. При синтезировании неевклидовой отсчетной формы последние три поля заменяются на более общие геометрические поля.

3 Хотя в оригинальных работах Картана рассматривались только вращения, можно предположить, что эти преобразования являются достаточно общими.

4 В явном виде этот изоморфизм может быть определен следующим образом. Если A — элемент HomV; EndV, то ему отвечает элемент LA пространства L2V,V;V, заданный равенством LAu, v:=Auv. Искомый изоморфизм является соответствием ALA.

5 Здесь Coor обозначает декартову арифметизацию [1, формула (2.2)].

6 Отображение σR есть координатное отображение, соответствующее криволинейным координатам QII=13 на форме SR [1, формула (2.13)].

7 Формула (2.4) — производная тензорного поля — содержит производные скалярных и векторных базисных полей. Производные базисных полей, в свою очередь, выражаются через исходные базисные поля; коэффициенты разложения есть символы Кристоффеля соответствующей системы координат.

8 Поскольку физическое пространство евклидово, а размерности форм и пространства совпадают.

9 Несмотря на то что в рамках настоящей работы тензоры N1 и N2 предполагаются постоянными, ход рассуждений не изменится, если эти тензоры заменить на переменные тензоры, зависящие от точки XSR. В таком случае правые части соотношений (2.7) и им подобным следует заменить на N1X и N2X соответственно. Такая замена целесообразна, если натуральное состояние определяется в рамках действия некоторого внешнего поля, которым нельзя пренебречь.

10 Это всегда можно сделать, используя операцию сдвига в 3.

11 Для вычисления производной DH1 можно воспользоваться равенством (2.3), где AIi=HIi, а q=x.

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, 105005

References

  1. 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.
  2. 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.)
  3. 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.
  4. 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.
  5. 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.
  6. 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.
  7. 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.
  8. 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)
  9. 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.
  10. 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.)
  11. 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.)
  12. 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.)
  13. Kartan E. Geometry of Riemannian spaces. Moscow: Knizhnyi dom «Librokom», 2010, 248 p. Available at: https://reallib.org/reader?file=444675. (In Russ.)
  14. Lee J.M. Introduction to Smooth Manifolds. New York: Springer, 2012, 708 p. DOI: https://doi.org/10.1007/978-1-4419-9982-5.
  15. 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.
  16. 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.
  17. Schwartz L. Analysis. Vol. 1. Moscow: Mir, 1972, 824 p. Available at: https://klex.ru/14z8. (In Russ.)
  18. 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.)
  19. [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.
  20. Maugin G.A. Material inhomogeneities in elasticity. New York: CRC Press, 1993, 292 p. DOI: https://doi.org/10.1201/9781003059882.
  21. Marsden J.E., Hughes T.J. Mathematical foundations of elasticity. New York: Courier Corporation, 1994. 576 p. URL: https://archive.org/details/mathematicalfoun00mars.
  22. 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.
  23. 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.)
  24. 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.
  25. 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.
  26. Seifert H. Topologie Dreidimensionaler Gefaserter Räume. Acta Mathematica, 1933, vol. 60, pp. 147–238. DOI: https://doi.org/10.1007/BF02398271.
  27. 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.
  28. 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.
  29. 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.)
  30. 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.)
  31. 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.
  32. 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.
  33. Lurie A.I. Analytical Mechanics. Moscow: GIFML, 1961, 824 p. Available at: https://eqworld.ipmnet.ru/ru/library/books/Lure1961ru.djvu. (In Russ.)
  34. 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

Supplementary Files
Action
1. JATS XML

Copyright (c) 2023 Koifman K.G.

Creative Commons License
This work is licensed under a Creative Commons Attribution 4.0 International License.

Согласие на обработку персональных данных с помощью сервиса «Яндекс.Метрика»

1. Я (далее – «Пользователь» или «Субъект персональных данных»), осуществляя использование сайта https://journals.rcsi.science/ (далее – «Сайт»), подтверждая свою полную дееспособность даю согласие на обработку персональных данных с использованием средств автоматизации Оператору - федеральному государственному бюджетному учреждению «Российский центр научной информации» (РЦНИ), далее – «Оператор», расположенному по адресу: 119991, г. Москва, Ленинский просп., д.32А, со следующими условиями.

2. Категории обрабатываемых данных: файлы «cookies» (куки-файлы). Файлы «cookie» – это небольшой текстовый файл, который веб-сервер может хранить в браузере Пользователя. Данные файлы веб-сервер загружает на устройство Пользователя при посещении им Сайта. При каждом следующем посещении Пользователем Сайта «cookie» файлы отправляются на Сайт Оператора. Данные файлы позволяют Сайту распознавать устройство Пользователя. Содержимое такого файла может как относиться, так и не относиться к персональным данным, в зависимости от того, содержит ли такой файл персональные данные или содержит обезличенные технические данные.

3. Цель обработки персональных данных: анализ пользовательской активности с помощью сервиса «Яндекс.Метрика».

4. Категории субъектов персональных данных: все Пользователи Сайта, которые дали согласие на обработку файлов «cookie».

5. Способы обработки: сбор, запись, систематизация, накопление, хранение, уточнение (обновление, изменение), извлечение, использование, передача (доступ, предоставление), блокирование, удаление, уничтожение персональных данных.

6. Срок обработки и хранения: до получения от Субъекта персональных данных требования о прекращении обработки/отзыва согласия.

7. Способ отзыва: заявление об отзыве в письменном виде путём его направления на адрес электронной почты Оператора: info@rcsi.science или путем письменного обращения по юридическому адресу: 119991, г. Москва, Ленинский просп., д.32А

8. Субъект персональных данных вправе запретить своему оборудованию прием этих данных или ограничить прием этих данных. При отказе от получения таких данных или при ограничении приема данных некоторые функции Сайта могут работать некорректно. Субъект персональных данных обязуется сам настроить свое оборудование таким способом, чтобы оно обеспечивало адекватный его желаниям режим работы и уровень защиты данных файлов «cookie», Оператор не предоставляет технологических и правовых консультаций на темы подобного характера.

9. Порядок уничтожения персональных данных при достижении цели их обработки или при наступлении иных законных оснований определяется Оператором в соответствии с законодательством Российской Федерации.

10. Я согласен/согласна квалифицировать в качестве своей простой электронной подписи под настоящим Согласием и под Политикой обработки персональных данных выполнение мною следующего действия на сайте: https://journals.rcsi.science/ нажатие мною на интерфейсе с текстом: «Сайт использует сервис «Яндекс.Метрика» (который использует файлы «cookie») на элемент с текстом «Принять и продолжить».