CC BY
18
2021 Математика и механика № 70
УДК 539.3
Б01 10.17223/19988621/70/9
М.Ю. Соколова, Д.В. Христич
КОНЕЧНЫЕ ДЕФОРМАЦИИ НЕЛИНЕЙНО УПРУГИХ АНИЗОТРОПНЫХ МАТЕРИАЛОВ1
Рассмотрен вариант соотношений нелинейной упругости для анизотропных материалов, по типу симметрии относящихся к кристаллам кубической син-гонии. В предложенной модели учитывается физическая нелинейность в поведении таких материалов при конечных деформациях. Исходя из представления упругого потенциала как тензорного многочлена по деформациям, получены соотношения для напряжений, содержащие константы упругости второго и третьего порядков.
Ключевые слова: анизотропия, гиперупругость, конечные деформации, кубические материалы, тензорные базисы, инварианты
Будем рассматривать анизотропные материалы, обладающие симметрией упругих свойств, присущей кристаллам кубической сингонии [1 - 4]. Это означает, что упругие свойства таких материалов удовлетворяют условиям симметрии, присущей точечной группе объемно- или гранецентрированного куба. Группа симметрии кубических материалов характеризуется наличием трех поворотных осей четвертого порядка, четырех поворотных осей третьего порядка и шести осей симметрии второго порядка. Порождающими элементами группы симметрии кубических материалов являются три поворота на угол 90° вокруг поворотных осей четвертого порядка [3, 4].
По своим свойствам кубические материалы близки к изотропным материалам. Известно [2, 3], что под действием гидростатического давления сфера из анизотропного материала в общем случае становится эллипсоидом. В случаях изотропного и кубического материалов при воздействии гидростатического давления сферы остаются сферами, что не позволяет в таком опыте различить эти материалы.
Линейные упругие кубические материалы в рамках обобщенного закона Гука описаны в работах [2, 6-10]. В этих работах получены структурные представления тензоров упругости четвертого ранга, инвариантных относительно описанной группы симметрии. Тензоры упругости, записанные в произвольной (лабораторной) системе координат, имеют в общем случае 21 ненулевую компоненту, которые не являются независимыми. В работах [6-9] проведен анализ зависимости модулей упругости кубического материала и коэффициента Пуассона от направления растяжения образца. Для кубического материала может быть определена такая система координат, названная в статье [11] канонической, в которой тензор упругих свойств имеет три ненулевые независимые константы.
Нелинейные модели поведения кубических материалов могут учитывать либо геометрическую, либо физическую нелинейность. В наиболее сложных моделях
1 Работа выполнена при частичной поддержке гранта Президента Российской Федерации (проект МД-1803.2019.1) и РФФИ (проект № 18-31-20053).
необходимо учитывать геометрическую и физическую нелинейности одновременно. Построению моделей кубических материалов, учитывающих физическую нелинейность при конечных деформациях, посвящена работа [12]. В этой статье автор записывает нелинейные определяющие соотношения для кубического материала на основе девяти тензорных генераторов, построенных по тензору конечных деформаций Коши - Грина и полученных в работе [5].
Если при упругом деформировании в рассматриваемых материалах наблюдаются нелинейные эффекты даже в области малых деформаций, то требуется построить физически нелинейные определяющие соотношения. Решению этой задачи посвящены статьи [13, 14]. В этих работах для записи нелинейных определяющих соотношений используется разложение в ряд упругого потенциала с сохранением членов второй и третьей степеней относительно тензора малых деформаций. В статьях [13, 14] определена структура тензоров упругости шестого ранга, которые в канонической системе координат содержат шесть ненулевых независимых констант. Геометрически и физически нелинейная модель кубического материала предложена в [15].
В отличие от результатов, полученных в работах [13-15], авторы данной статьи используют разложение тензоров упругости четвертого и шестого рангов по собственным упругим состояниям кубического материала [13, 14]. Это позволяет записать упругий потенциал для кубического материала в виде функции инвариантов тензоров, которые являются проекциями тензора деформаций Коши - Грина в собственные подпространства кубического материала. Получающиеся при этом выражения для тензора напряжений отражают взаимное влияние процессов, происходящих в различных собственных подпространствах рассматриваемого материала, и позволяют для кубического материала описать эффекты второго порядка.
1. Общий подход к записи упругого потенциала деформаций в случае конечных деформаций кубического материала.
Строение тензоров упругости
Рассмотрим гиперупругий анизотропный материал, для которого можно записать упругий потенциал. В качестве такового используем удельную (отнесенную к объему) потенциальную энергию деформации. Дифференциал удельной потенциальной энергии деформаций может быть представлен в виде
йШ = Т: йг,
где Т - энергетический тензор напряжений, связанный с тензором напряжений Коши 8 соотношением
т 1 ё¥ Т = Ф-т • Е • Ф-1, £ = — 8 ,
Ф - аффинор деформаций, г = 1 (ф • Фт - е) - тензор деформаций Коши - Грина,
Е - единичный тензор, двоеточие означает свертывание тензоров [3, 16]. Из этого следует возможность определения напряжений при конкретизации вида удельной потенциальной энергии деформаций по формулам
дШ
т=~г. (1
дг
Такой подход к построению соотношений гиперупругости для изотропных и анизотропных материалов использовался в работах [16 - 18].
Представим тензорную функцию Ж (е) в виде ряда по степеням тензора деформаций Коши - Грина
Ж = Ж0 + Л0 : е+—N :: ее+—Ь::: еее +..., (2)
0 0 2! 3!
причем Ж0 = 0, Л0 = 0, если начальное состояние является ненапряженным.
В выражении (2) использовано произведение ее = е ] ек1ёё]ёкё1, где
ёё]ёкё1 = е ® ё] ® ёк ® ё1 - полиада, образованная векторами ортонормированно-
го базиса ё^, /' = 1,2,3, и произведение еее = егуек1 етпёё]ёкё1ётёп, где
ёгё]ёкё1ётёп = ё ® ё] ® ёк ® ё1 ® ёт ® ёп .
Сохраним в представлении (2) только первые два ненулевых члена:
Ж = — N :: ее+—Ь ::: еее , (3)
2! 3!
тогда из соотношений (1) и (3) следует выражение для напряжений
Т = N: е +1Ь :: ее . (4)
2
В выражениях (3) и (4) N и Ь - тензоры упругих констант четвертого и шестого рангов соответственно, которые удовлетворяют условиям внутренней симметрии
N]к1 = N]Ш = N]1к = ^кИ] , Ц}к1тп = ^]Мтп = ^у1ктп = ^уШт = ^утМ = ^Ыутп . (5)
Структура тензоров N и Ь для кубического материала известна [1 - 4, 16]. Тензор N содержит три независимые константы упругости второго порядка, а тензор Ь - шесть независимых констант упругости третьего порядка. Наименьшее число ненулевых компонент тензоры упругости имеют в системе канонических осей анизотропии материала [11, 17]. В произвольной (лабораторной) системе координат тензоры упругости кубического материала имеют произвольный вид. Главные оси анизотропии по В.В. Новожилову [11, 16] определяются как главные оси тензора напряжений при всестороннем сжатии, однако для кубического материала все главные значения тензора напряжений в этом случае равны, а главные векторы могут быть выбраны произвольно. Канонические оси анизотропии материала всегда совпадают с главными осями анизотропии, определенными по В.В. Новожилову. Однако в кубическом материале главные оси анизотропии, найденные из эксперимента, могут и не совпадать с каноническими осями. Тогда возникает необходимость определить взаимную ориентацию лабораторной системы координат с ортонормированным базисом к и системы канонических осей координат с базисом а .
Поскольку реализация эксперимента по всестороннему сжатию образца затруднена, в работе [11] предложено заменить этот эксперимент на три опыта по сжатию кубического образца в трех взаимно перпендикулярных направлениях. В этих опытах обязательным является измерение всех компонент тензора деформаций в лабораторной системе координат. Эквивалентность таких экспериментов
является следствием линейности связи между деформациями и напряжениями в области малых деформаций. В этом случае речь идет об определении начального положения осей анизотропии материала.
Пусть опыты по сжатию кубических образцов проводятся в лабораторной системе координат с ортонормированным базисом к1. Взаимная ориентация векторных базисов а1 и к' определяется ортогональным тензором поворота О. Для компонент тензора О = qijk1k1 выполняются тождества
+ + ^31 = 1 + ?22 + ^32 = 1 ?1з + ?23 + ?3з = 1,
Я\\Я\2 + ?21?22 + ?31?32 = 0 Ч\\Ч\3 + ?21?23 + <?31?33 = 0, (6)
?13?12 + ?23?22 + ?33?32 =
В работе [11] показано, что для определения положения канонических осей анизотропии в кубическом материале достаточно двух экспериментов на сжатие кубических образцов. В первом эксперименте тензор напряжений определяется как Т1 = —к 1к1. Пусть С - тензор упругих податливостей, обратный к тензору упругости N . Измеряемые деформации выражаются через константы податливости и компоненты тензора О следующим образом:
811 = [01 (С1111 -(С1122 + 2С1212 )) + С1122 + 2С1212 ] ,
812 = [07 (С1111 -(С1122 + 2С1212 ))] , 822 = [04 (С1111 -(С1122 + 2С1212 )) + С1122 ] ,
% = -/ [б8 (С1111 -(С1122 + 2С1212 ))] , (7) е33 = [05 (С1111 -(С1122 + 2С1212 )) + С1122 ] , 823 = [09 (С1111 -(С1122 + 2С1212 ))] ,
где обозначено
01 = + ?142 + Й , 04 = ?П?21 + ?122?22 + ?1з?23 ,
05 = ?П?321 + ?12?322 + ?13?323 , 07 = ?21?П + ?22Й + ?23?1з ,
3 3 3 2 2 2
08 = ?31?11 + ?32?12 + ?33?13, 09 = ?21?31?11 + ?22?32?12 + ?23?33?13.
Во втором эксперименте тензор напряжений Т2 = -/к2к2 и измеряемые компоненты тензора деформаций также выражаются через константы податливости и компоненты тензора О :
811 = [04 (С1111 -(С1122 + 2С1212 )) + С1122 + 2С1212 ] , 812 = [010 (С1111 -(С1122 + 2С1212 ))], 822 = [02 (С1111 -(С1122 + 2С1212 )) + С1122 ] ,
823 =-/ [011 (С1111 -(С1122 + 2С1212 ))] (8)
833 = [06 (С1111 -(С1122 + 2С1212 С1122 ] ,
813 = [012 (С1111 -(С1122 + 2С1212 ))] ,
где обозначено
б2 = Ч21 + Ч22 + Ч23, бб = Ч31Ч21 + ?32?22 + ?33?23' 010 = ?11?21 + Ч12Ч22 + ?13?23'
3 3 3 2 2 2
011 = Ч31Ч21 + Ч32Ч22 + Ч33Ч23 , 09 = Ч11Ч31Ч21 + Ч12Ч32Ч22 + Ч13Ч33Ч23 •
Для отыскания девяти компонент тензора О используем шесть соотношений (б) и четыре независимых соотношения из (7), (8):
[07 (С1111 -(С1122 + 2С1212 ))] = &12 , —[08 (С1111 -(С1122 + 2С1212 ))] = &13 ,
— [09 (С1111 -(С1122 + 2С1212 ))] = &23 ,
— [01О (С1111 -(С1122 + 2С1212 ))] = &12 •
Исключая из них множитель (С1111 -(С1122 + 2С1212)), получим три уравнения относительно компонент тензора О :
ЧиЧп + Ч32Ч132 + Ч33Ч133 = зг3(?21?П + Ч22Ч132 + Ч23Ч133 ), &12
Ч21Ч31ЧП + Ч22Ч32 Ч122 + Ч23Ч33Ч123 = I21 (Ч21Ч31 + Ч22 Ч132 + Ч23Ч133 ), (9)
& &12
3 3 3 _ ь12 / 3 3 3 \
Ч11Ч21 + Ч12 Ч22 + Ч13 Ч23 уЧ21Ч11 + Ч22 Ч12 + Ч23Ч13 ^
&12
Численное решение системы уравнений (б), (9) позволяет найти компоненты тензора О, то есть определить ориентацию канонической системы координат в кубическом материале относительно лабораторной системы координат по измеряемым в опытах деформациям.
Вычисления показали, что в случае кубических кристаллов канонические оси анизотропии совпадают с их кристаллографическими осями [100], [010] и [001] (обозначения из [1, 4]), а для композитных материалов или древесины их положение совпадает с преимущественными структурными направлениями.
Введем в рассмотрение тензорный базис, образованный диадами базисных векторов канонических осей симметрии кубического материала а{ :
А1 = а1а1, А2 = а2а2 , А3 = а3а3, А4 = (а3а2 + а2а3),
л/2
А5 = —^(а2а3 + а3а2), Аб = (З3а1 + а1а3). (10)
л/2 л/2
Базис (10) нормируется соотношением: Аг : А1 =5У , где 5У - дельта Кронекера.
Наряду с базисом (10) рассмотрим тензорный базис А.А. Ильюшина 1а (а = 0,1,...,5) с базисными тензорами [11, 16, 17]:
I0 = —^г(а + а2а2 + а3а3), I1 = —^(2а3а3 - а1а1 - а2а2), I2 = —^(31а1 - а2а2),
V б л/ 2
I3 =—^ (а3а2 + а2 а3), 14 =—^ (а2 а3 + а3а2), I5 =—^ (¿23 а3 + а1а3). (11) 2 2 2
Тензорный базис (11) также нормирован: Iа : Iе = 5ав.
Тензор деформаций £ = е 1]-а1а]- как симметричный тензор второго ранга можно разложить по базисам (10) и (11). Эти разложения имеют вид
£ = е11А1 + е22 А2 + е33А3 +л/2е12 А4 + л/2е23А5 + \/2е31А6 (12)
и £ = е 010 + е111 +е212 + е 313 + е 414 +е515, (13)
где еа = £: 1а .
Между коэффициентами разложений (12) и (13) имеется связь:
е0 = (11 + е22 + е33 ) , е1 = (2833 - е11 - е22 ) , е2 = (11 - е22 ) ,
е3 = л/2е12 , е4 = >/2е23, е5 = л/2е31. Обратная связь имеет вид
= __1_ = __1___1_ = ¡2
е11 = Т3е0 -Тбе1 , е22 =^не° -Тбе1 ^^ е33 = Т3е0 "V3^ == =
е12 = л/2 е3 , = •ч/^ е31 = По тензорам второго ранга (10) образуем тензоры четвертого и шестого рангов, обладающие внутренней симметрией (5):
А«Р = 2 (а Ав + Ав Аа),
АаРу = 1 (АаАвАу + АвАаАу + Ау Аа Ав + АаАуАв + АвАуАа + Ау Ав Аа), (14) где а, р,у = 1,2,...,6 .
Из работ [1 - 4, 15, 16] известно разложение тензоров упругости N и Ь по базисам (14). Для кубического материала эти разложения имеют вид
N = п11 (А11 + А22 + А33) + п12 (А12 + А23 + А13) + п44 (А44 + А55 + А66); (15)
Ь = с, (А111 + А222 + А333 ) + с2 (А155 + А266 + А344) + + с3 (А112 + А113 + А122 + А133 + А223 + А233 ) + + с4 (А144 + А166 + А255 + А244 + А355 + А366 )5А123 + с6А456, (16)
где п11, п12, п44 - ненулевые компоненты тензора упругости N - константы упругости второго порядка; с1, с2, с3, с4, с5, с6 - константы упругости третьего порядка.
По тензорам (11) построим базисы, состоящие из тензоров четвертого (1ав) и :того (!а|3у) рангов:
1ав= 2 ( + 1в1а), (17)
1аРу = -6 ((Г + 1в1а V + V 1а1в + 1а1у Iе + Йу 1а + V 1в1а), (18)
6
где а, р, у = 0,1,...,5 .
Базисы (17) и (18) нормируются соотношениями:
1ав :: 1уе = 2(5ау5ре+5ае5ру),
/аРу:::15е? =1(5а?5ре5у5 + 5ае5к5у5 +5а55к5уе +5а?5р55уе +5а55ре5у? +5ае5р55у?).
В работах [11, 17] получено разложение тензора упругости N по базису (5) в виде
N = «(1)1°° + п(2) (I11 +122 ) + п(3) (I33 +144 +155) . (19)
Коэффициенты в представлении (19) связаны с константами упругости второго порядка (15) соотношениями
п(1) = п11 + 2п12, п(2) = п11 -п12, п(3) = п44.
Тензор констант упругости третьего порядка Ь может быть разложен по базису (18), однако в своем разложении он должен содержать только такие комбинации тензоров 1аРу, которые являются инвариантными относительно группы ортогональных преобразований кубической сингонии. Методом прямой проверки установлено, что для кубического материала имеется шесть инвариантных комбинаций тензоров (18):
в(1) = 1000 в(2) = 1011 +1022 в(3) = 1033 + ^44 + 1055
В(4) =-^(1ш -31122), в(5) =4= (I144 +1155 -21133255 -1244), в(6) =-Ъ345.
' 46х ' ' л/2
Разложение тензора Ь по базису (20) имеет вид
Ь = § Ъ в( я). (21)
я=1
В выражении (21) коэффициенты Ъя связаны с константами упругости третьего порядка кубического материала (16) соотношениями
Ъ1 = —^ (с1 + 6с3 + 2с5), Ъ2 = л/3 (с - с5), Ъ3 = 2\/3 (с2 + 2с4), л/3
Ъ4 = с - 3с3 + 2с5, Ъ5 = 6 (с4 - с2), Ъ6 = 24с6 .
Разложения тензоров упругости (19) и (21) далее использованы для записи упругого потенциала для кубического материала и получения нелинейных определяющих соотношений.
2. Определяющие соотношения в собственных подпространствах кубического материала
Если в представлении для удельной потенциальной энергии деформаций (2) ограничиться только членом второго порядка, то соотношения (4) принимают вид обобщенного закона Гука:
Т = N : £ . (22)
Учтем в соотношениях (22) представление (19) и получим
Т = ««I00 : £ + «(2) (I11 +122): £ + «(3) (I33 +144 +155): £
или Т = «(1)£(1) + «(2)£(2) + «(3)£(3) . (23)
Из соотношений (19) и (23) следует, что тензорный базис (17) является собственным для кубического материала. Понятие о собственных тензорах и собственных состояниях введено в работах Рыхлевского, Ковина и Махрабади [19, 20]. Собственные тензоры для кубического материала получены в работах [11, 12, 1б, 21].
Тензоры деформаций £(2), £(3) являются собственными упругими состояниями кубического материала и принадлежат трем собственным подпространствам:
- одномерному (1Б) - = е010;
- двумерному (2Б) - £(2) = е^1 + е^2 ; (24)
- трехмерному (3Б) - £(3) = е^3 +&Д4 +е5!5. Спроектируем тензор Т в те же собственные подпространства:
Т(1) = ТI0, Т(2) = Щ1 + T212, Т(3) = 72I3 + T414 + Т515, (25)
где Та= Т: Г* .
В соответствии с (23) закон Гука можно записать в виде
Т = Х Т(,) = £ п( к) £(к). (2б)
к=1 к=1
Из представления (2б) следует, что в рамках линейной упругости тензоры напряжений Т(а) и деформаций £(а) в каждом собственном подпространстве соос-
ны и пропорциональны.
Конкретизируем соотношения (4) для кубического материала. Для этого выясним свойства тензоров В(ж), входящих в разложение (21) тензора Ь . Эти свойства проявляются при вычислении сверток тензоров В( ж) с собственными тензорами деформаций (24):
£ : В( ^) : £ = £ : В^ ) : £
£(а) : В : £(Р) " £(Р) : В : £(а) .
Из этих произведений ненулевыми являются следующие девять:
£(1) : В() : £(1) = , £(1) : В( ) : £(2) = Ш(£(1))£(2) ,
£(2) : В( ) : £(2) = (£22) , £(1) : В( ) : £(3) = Щ — (£(1))£(3) ,
£(3) : В( ) : £(3) = (£(3))I , £(2) : В( ) : £(2) = 0(2)2 , (27)
(5) = 2 (5) = 2 (б) = 1
£(2) : В : £(3) = 1 Р(3), £(3) : В : £(3) = 3 0(3)2, £(3) : В : £(3) = 3 0(3)3.
В соотношениях (27) обозначено — (А) = А: Е - первый инвариант тензора А ; 0(2)2 - проекция тензора £22) во второе собственное подпространство:
£22) =^3(&2 +е2)I0 + 0(2)2, 0(2)2 (-е2)I1 -^е^^^ ^ (28)
0(3)2 и 0(3)3 - проекции тензора £23) во второе и третье собственные подпространства:
£23) = ( + е 4 + ) I0 + 0(3)2 + 0(:
(3)2
{3) - ^3^3 ^4 I "г°(3)^ °(3)3,
^ (2е32-е?-е2 )I1 - ^ (е4-е52 )I2, (29)
>(3)3 = 3 е4е5I + е3е5I + ^^ ;
тензор
р(3) = ^ + 2 (зб е1 -72е 2 К+2 (зб е1 +72 е2 Jе5I5. (30)
В данной модели выражение для удельной потенциальной энергии деформации имеет вид
Ж = 1 ((1)./2 (£(1)) + «(2)./2 (£(2) ) + «(3)./2 (£(3))) +
+ б^Э /1(£(1)) (9Й1-2 (£(1) ) + Ь2-2(£(2) ) + Ь3-2(£(3) )) +
+ б -3 (£(2) ) (( - 2Ь5 ) + 31б -3 (£(3) ) (Ьб - 12Ь5 ) + 3 Ь5-3 (£(2) + £(3) X С3)
где введены обозначения для второго и третьего инвариантов тензора второго ранга: /2 (А) = А: А , /3 (А) = det А.
Запишем выражения для инвариантов тензоров, входящих в соотношения (27) и (31), через коэффициенты еа разложения (13). Инварианты тензоров £(3), £(2),
£(3) запишутся как
-/1(£(1)) = 1 е0, -/2(£(1)) = 1 е2 , -/3(£(1)) = Ще0 ,
-/1(£(2)) = 0, -2(£(2)) =е12 + е2 , /3(£(2)) е1 (е12 - 3е2 ), С32»
,2,2,2-- - 1
-/1(£(3)) = ^ -/2(£(3)) =е^ + е2 +е2 , -/3(£(3)) = /Ге3е4е5 .
Ж
Инварианты тензоров е/^ , е/2) , е/3) :
=•/1(£(22))=3( + е2)' •/1(£2З)>=3(( + е4 +е2)• (33)
Смешанный инвариант тензоров Е(2), Е(3):
• (е(2) + Е(3) ) = Е(2) : 0(3)2 + •з (е(2) ) + •З (Е(3) ) = Е(3) : Р(3) + •3 (е(2) ) + ^3 (е(3) ) =
=386 ( - 3е2 )+;/28З8485 + 2/6 (+82 - 282)+2/2 (- е2 )• (34)
Присутствие смешанного инварианта (34) в выражении для удельной потенциальной энергии деформаций позволяет учитывать взаимное влияние процессов, происходящих в собственных подпространствах 2Б и 3Б. Отметим, что выражения для смешанного инварианта (34) упрощаются в случаях, когда деформации во втором или третьем собственном подпространствах отсутствуют, поскольку при этом Е(2) : 0(3)2 = Е(3) : Р(3) = 0. Отметим, что
если Е(3) = 0, то
если Е(2) = 0 , то
• (е(2) + Е(3) ) = •з (е(2) ) = 3^6 ( 382 ) '
• (е(2) + Е(3)) = •^3(е(3)) = г: 838485 •
л//1
Подставляя (19) и (21) в определяющие соотношения (4), получим следующую форму связи между напряжениями и деформациями:
Т(1) =
(П +Ь180 )о + 3 Ь2 (( + 8 2 )+3 Ь3 (( +82 + 8 2)
Т(2) = (п2 + 1Ь280 ] Е(2) + Ь40(2)2 + 2Ь50(3)2 , (35)
3 2 1
Т(3) = п Е(3) + 3Ь5Р(3) + 3Ь6 0(3)3 .
Соотношения (35) содержат деформации во второй степени и являются физически нелинейными соотношениями. В этом случае в неодномерных собственных подпространствах 2Б и 3Б тензоры напряжений и деформаций Т(2) и Е(2), Т(3) и
Е(3) перестают быть соосными. Во втором подпространстве отклонение от соосности тензоров Т(2) и Е(2) связано как с появлением составляющей напряжений вдоль тензора 0(2)2, так и составляющей вдоль 0(3)2. В третьем подпространстве отклонение от соосности тензоров Т(3) и Е(3) связано с появлением составляющих напряжений вдоль тензоров 0(3)3 и Р(3), причем последний обращается в ноль, если Е(2) = 0 .
Анализ показывает, что полученные соотношения (35) не удовлетворяют обобщению частного постулата А.А. Ильюшина на анизотропные материалы, сформулированному в работах [11, 1б, 17], и учитывают взаимное влияние процессов, происходящих в различных собственных подпространствах.
3. Анализ нелинейных эффектов, описываемых моделью
Соотношения (35) описывают связь между конечными деформациями кубического материала и напряжениями. Для анализа нелинейных эффектов, описываемых соотношениями (35), будем считать деформации малыми. В этом случае собственные состояния кубического материала £(3), £(2), £(3) имеют простой физический смысл. Первое собственное состояние £(3) является чисто объемным деформированием. Второе состояние £(2) соответствует формоизменению, происходящему в главных осях анизотропии материала и связанному только с изменением длин материальных волокон. Третье собственное состояние £(3) соответствует
чистым сдвигам в координатных плоскостях.
Рассмотрим процесс деформирования, целиком расположенный в первом собственном подпространстве: £(3) =е010, £(2) =£(3) =0. В соответствии с (35) в ответ на такие деформации в кубических материалах появляются напряжения Т(3) =(« + Ь1е0 )е010, которые являются гидростатическими, нелинейно зависящими от объемных деформаций е0. Касательные напряжения в главных осях анизотропии не появляются.
Пусть процесс деформирования целиком расположен во втором собственном подпространстве £ = £(2) =е1!1 +е2!2, £(3) = £(3) = 0. Такому процессу соответствует изменение длин волокон, расположенных вдоль главных осей анизотропии, без сдвигов. В это случае в соответствии с (35) возникающие напряжения имеют вид Т = Т(3) + Т(2), причем
Т(1) = 1Ь2 (( +е2 ) Т(2) = («2 + 3Ь2е0 ) £(2) + Ь40(2)2 ,
то есть соотношения описывают нелинейную зависимость напряжений от деформаций. В процессе формоизменения £(2) появляются гидростатические напряжения. Как и в первом случае, касательные напряжения Т(3) в таком процессе не появляются.
Если процесс деформирования £ = £(3) =е3I3 +е414 +е5I5, £(3) = £(2) = 0 заключается в чистых сдвигах хотя бы в одной из плоскостей, содержащих канонические оси анизотропии, то в соответствии с (35) возникают напряжения
Т = Т(1) + Т(2) + Т(3) , причем
Т(1) = 3 Ь3 (е5 +е2 +е5 )I0, Т(2) = 2Ь5 ^^(3)2 , Т(3) = «3£(3) +3 Ьб0(3)3 ,
то есть соотношения (35) описывают нелинейную зависимость касательных на-
пряжений T(3) от сдвиговых деформаций Е(3) и прогнозируют возникновение нормальных напряжений, в том числе и гидростатических.
Заключение
Для кубического материала получено разложение тензоров упругости четвертого и шестого рангов по тензорным базисам в собственных подпространствах. Из условия существования упругого потенциала (удельной потенциальной энергии деформации) получены соотношения между напряжениями и конечными деформациями, содержащие деформации во второй степени. Выписаны выражения для напряжений в каждом из собственных подпространств кубического материала.
Предложенный вариант соотношений позволяет учесть взаимное влияние процессов в различных собственных подпространствах кубического материала. Если в предложенных соотношениях обнулить константы упругости третьего порядка, то они будут сведены к линейным соотношениям закона Гука для кубического материала. При сохранении в выражениях упругих констант третьего порядка показана непропорциональность тензоров напряжений и деформаций в каждом из собственных подпространств.
ЛИТЕРАТУРА
1. Сиротин Ю.И., Шаскольская М.П. Основы кристаллофизики. М.: Наука, 1979. 639 с.
2. Лехницкий С.Г. Теория упругости анизотропного тела. М.: Наука, 1977. 416 с.
3. ЧерныхК.Ф. Введение в анизотропную упругость. М.: Наука, 1988. 192 с.
4. Newnham R.E. Properties of Materials: Anisotropy, Symmetry, Structure. New York: Oxford University Press, 2005. 391 p.
5. Xiao H. А new representation theorem for elastic constitutive equations of cubic crystals // Journal of Elasticity. 1999. V. 53. P. 37-45. DOI: 10.1023/A:1007591025837.
6. Paszkiewicz T., Wolski S. Elastic properties of cubic crystals: Every's versus Blackman's diagram // Journal of Physics: Conference Series. 2008. V. 104. 012038. DOI: 10.1088/17426596/104/1/012038.
7. Knowles K.M. The biaxial moduli of cubic materials subjected to an equi-biaxial elastic strain // Journal of Elasticity. 2016. V. 124. P. 1-25. DOI: 10.1007/s10659-015-9558-x.
8. Knowles K.M., Howie P.R. The directional dependence of elastic stiffness and compliance shear coefficients and shear moduli in cubic materials // Journal of Elasticity. 2015. V. 120. P. 87-108. DOI: 10.1007/s10659-014-9506-1.
9. Norris A. Poisson's Ratio in Cubic Materials // Proceedings: Mathematical, Physical and Engineering Sciences. 2006. V. 462. No. 2075. Р. 3385-3405. DOI: 10.1098/rspa.2006.1726.
10. Duffy T. Single-crystal elastic properties of minerals and related materials with cubic symmetry // American Mineralogist. 2018. V. 103. Iss. 6. P. 977-988. DOI: 10.2138/am-2018-6285.
11. Соколова М.Ю., Христич Д.В. Программа экспериментов по определению типа начальной упругой анизотропии материала // Прикладная механика и техническая физика. 2015. Т. 56. № 5. С. 205-213. DOI: 10.15372/PMTF20150519.
12. Wright T.W. Bootstrap elasticity III: minimal nonlinear constitutive representation for cubic materials // Journal of Elasticity. 2015. V. 120. No. 1. P. 109-119. DOI: 10.1007/s10659-014-9507-0.
13. Claiton J.D. Nonlinear Elastic and Inelastic Models for Shock Compression of Crystalline Solids. Springer, 2019. 452 p.
14. Kambouchev N., Fernandez J., Radovitzky R. A poly convex model for materials with cubic symmetry // Modelling and Simulation in Materials Science and Engineering. 2007. V. 15. No. 5. P. 451-468. DOI: 10.1088/0965-0393/15/5/006.
15. Kube C.N., Turner J.A. Estimates of Nonlinear Elastic Constants and Acoustic Nonlinearity Parameters for Textured Polycrystals // Journal of Elasticity. 2016. V. 122. No. 2. P. 157— 177. DOI: 10.1007/s10659-015-9538-1.
16. Маркин А.А., Соколова М.Ю. Термомеханика упругопластического деформирования. М.: ФИЗМАТЛИТ, 2013. 320 с.
17. Соколова М.Ю., Христич Д.В. О симметрии термоупругих свойств квазикристаллов // Прикладная математика и механика. 2014. Т. 78. Вып. 5. С. 728-734.
18. Козлов В.В., Маркин А.А. Апробация определяющих соотношений нелинейной теории упругости при осевом сдвиге полого цилиндра // Вестник Томского государственного университета. Математика и механика. 2020. № 63. С. 102-114. DOI: 10.17223/ 19988621/63/9
19. Рыхлевский Я. О законе Гука // Прикладная математика и механика. 1984. Т. 48. Вып. 3. С. 420-435.
20. Mehrabadi M.M., Cowin S.C. Eigentensors of linear anisotropic elastic materials // The Quarterly Journal of Mechanics and Applied Mathematics. 1991. V. 44. Iss. 2. P. 331. DOI: 10.1093/qjmam/44.2.331.
21. Остросаблин Н.И. Об уравнениях линейной теории упругости // Прикладная механика и техническая физика. 1992. Вып. 3. С. 131-140.
Статья поступила 29.02.2020
Sokolova M.Yu., Khristich D.V. (2021) FINITE STRAINS OF NONLINEAR ELASTIC ANISOTROPIC MATERIALS. Vestnik Tomskogo gosudarstvennogo universiteta. Matematika i mekhanika [Tomsk State University Journal of Mathematics and Mechanics]. 70. pp. 103-116
DOI 10.17223/19988621/70/9
Keywords: anisotropy, hyperelasticity, finite strains, cubic materials, tensor bases, invariants.
Anisotropic materials with the symmetry of elastic properties inherent in crystals of cubic syngony are considered. Cubic materials are close to isotropic ones by their mechanical properties. For a cubic material, the elasticity tensor written in an arbitrary (laboratory) coordinate system, in the general case, has 21 non-zero components that are not independent. An experimental method is proposed for determining such a coordinate system, called canonical, in which a tensor of elastic properties includes only three nonzero independent constants.
The nonlinear model of the mechanical behavior of cubic materials is developed, taking into account geometric and physical nonlinearities. The specific potential strain energy for a hyperelastic cubic material is written as a function of the tensor invariants, which are projections of the Cauchy-Green strain tensor into eigensubspaces of the cubic material.
Expansions of elasticity tensors of the fourth and sixth ranks in tensor bases in eigensubspaces are determined for the cubic material. Relations between stresses and finite strains containing the second degree of deformations are obtained. The expressions for the stress tensor reflect the mutual influence of the processes occurring in various eigensubspaces of the material under consideration.
Financial support. The reported study was partially funded by the grant from the President of the Russian Federation according to the research project MD-1803.2019.1 and by the grant from RFBR according to the research project No. 18-31-20053.
Marina Yu. SOKOLOVA (Doctor of Physics and Mathematics, Associate Professor, Tula State University, Tula, Russian Federation). E-mail: m.u.sokolova@gmail.com
Dmitriy V. KHRISTICH (Doctor of Physics and Mathematics, Associate Professor, Tula State University, Tula, Russian Federation). E-mail: dmitrykhristich@rambler.ru
REFERENCES
1. Sirotin Yu.I., Shaskol'skaya M.P. (1982) Osnovy kristallofisiki [Fundamentals of crystal physics]. Moscow: Mir.
2. Lekhnitskiy S.G. (1963) Theory of Elasticity of an Anisotropic Elastic Body. San Francisco: Holden Day.
3. Chernykh K.F. (1999) An Introduction to Modern Anisotropic Elasticity. Begell House.
4. Newnham R.E. (2005) Properties of Materials: Anisotropy, Symmetry, Structure. New York: Oxford University Press.
5. Xiao H. (1999) А new representation theorem for elastic constitutive equations of cubic crystals. Journal of Elasticity. 53. pp. 37-45. DOI: 10.1023/A:1007591025837.
6. Paszkiewicz T., Wolski S. (2008) Elastic properties of cubic crystals: Every's versus Blackman's diagram. Journal of Physics: Conference Series. 104. Article 012038. DOI: 10.1088/1742-6596/104/1/012038.
7. Knowles K.M. (2016) The biaxial moduli of cubic materials subjected to an equi-biaxial elastic strain. Journal of Elasticity. 124. pp. 1-25. DOI: 10.1007/s10659-015-9558-x.
8. Knowles K.M., Howie P.R. (2015) The directional dependence of elastic stiffness and compliance shear coefficients and shear moduli in cubic materials. Journal of Elasticity. 120. pp. 87-108. DOI: 10.1007/s10659-014-9506-1.
9. Norris A. (2006) Poisson's Ratio in Cubic Materials. Proceedings: Mathematical, Physical and Engineering Sciences. 462(2075). pp. 3385-3405. DOI: 10.1098/rspa.2006.1726.
10. Duffy T. (2018) Single-crystal elastic properties of minerals and related materials with cubic symmetry. American Mineralogist. 103(6). pp. 977-988. DOI: 10.2138/am-2018-6285.
11. Sokolova M.Yu., Khristich D.V. (2015) Program of experiments to determine the type of initial elastic anisotropy of material. Journal of Applied Mechanics and Technical Physics. 56(5). pp. 913-919. DOI: 10.1134/S0021894415050193.
12. Wright T.W. (2015) Bootstrap elasticity III: minimal nonlinear constitutive representation for cubic materials. Journal of Elasticity. 120(1). pp. 109-119. DOI: 10.1007/s10659-014-9507-0.
13. Claiton J.D. (2019) Nonlinear Elastic and Inelastic Models for Shock Compression of Crystalline Solids. Springer.
14. Kambouchev N., Fernandez J., Radovitzky R. (2007) A polyconvex model for materials with cubic symmetry. Modelling and Simulation in Materials Science and Engineering. 15(5). pp. 451-468. DOI: 10.1088/0965-0393/15/5/006.
15. Kube C.N., Turner J.A. (2016) Estimates of nonlinear elastic constants and acoustic nonlinearity parameters for textured polycrystals. Journal of Elasticity. 122(2). pp. 157-177. DOI: 10.1007/s10659-015-9538-1.
16. Markin A.A., Sokolova M.Yu. (2015) Thermomechanics of elastoplastic deformation. Cambridge: Cambridge International Science Publishing.
17. Sokolova M.Yu., Khristich D.V. (2014) The symmetry of the thermoelastic properties of quasicrystals. Journal of Applied Mathematics and Mechanics. 78(5). pp. 524-528. DOI: 10.1016/j.jappmathmech.2015.03.012.
18. Kozlov V.V., Markin A.A. (2020) Aprobatsiya opredelyayushchikh sootnosheniy nelineynoy teorii uprugosti pri osevom sdvige pologo tsilindra [Testing of defining relations of nonlinear theory of elasticity in an axial strain of a hollow cylinder]. Vestnik Tomskogo gosudarstvennogo univrsiteta. Matematika i mekhanika - Tomsk State university Journal of Mathematics and Mechanics. 63. pp. 102-114. DOI: 10.17223/19988621/63/9.
19. Rychlewski J. (1984) On Hooke's law. Journal of Applied Mathematics and Mechanics. 48(3). pp. 303-314. DOI: 10.1016/0021-8928(84)90137-0.
20. Mehrabadi M.M., Cowin S.C. (1991) Eigentensors of linear anisotropic elastic materials. The Quarterly Journal of Mechanics and Applied Mathematics. 44(2). p. 331. DOI: 10.1093/ qjmam/44.2.331.
21. Ostrosablin N.I. (1992) Equations of the linear theory of elasticity. Journal of Applied Mechanics and Technical Physics. 33(3). pp. 438-446. DOI: 10.1007/BF00851743.
Received: February 29, 2020
