CC BY
18
Механика
DOI: 10.18721/JPM.14114 УДК 539.3
СООТНОШЕНИЯ МЕЖДУ МОДУЛЕМ ЮНГА И КОЭФФИЦИЕНТОМ ДИФФУЗИИ ДВУХФАЗНОГО МАТЕРИАЛА
К.П. Фролова
Институт проблем машиноведения РАН, Санкт-Петербург, Российская Федерация;
Санкт-Петербургский политехнический университет Петра Великого, Санкт-Петербург, Российская Федерация
В работе устанавливается взаимосвязь между изменениями эффективных упругих и диффузионных свойств двухфазного композита через микроструктурные параметры. Предполагается, что в материале присутствуют одинаковые по форме неоднородности. Представлен вывод соотношений в явном тензорном виде. При установлении взаимосвязи между эффективными свойствами учитывается эффект сегрегации, заключающийся в скачке концентрации растворенного вещества на границе раздела матрица/неоднородность. Полученные соотношения целесообразно использовать для определения одних эффективных свойств через другие, когда микроструктура материала неизвестна. Установленная взаимосвязь проверена для изотропного материала с порами. Найденные приближенные соотношения сравниваются с точными, полученными для конкретной микроструктуры.
Ключевые слова: эффективный модуль Юнга, эффективный коэффициент диффузии, взаимосвязь между свойствами композита
Ссылка при цитировании: Фролова К.П. Соотношения между модулем Юнга и коэффициентом диффузии двухфазного материала // Научно-технические ведомости СПбГПУ. Физико-математические науки. 2021. Т. 14. № 1. С. 177-189. DOI: 10.18721/JPM.14114
Статья открытого доступа, распространяемая по лицензии CC BY-NC 4.0 (https://creative-commons.org/licenses/by-nc/4.0/)
CROSS-PROPERTY CONNECTIONS BETWEEN YOUNG'S MODULUS AND DIFFUSION COEFFICIENT OF TWO-PHASE COMPOSITE
K.P. Frolova
Institute for Problems in Mechanical Engineering RAS, St. Petersburg, Russian Federation;
Peter the Great St. Petersburg Polytechnic University, St. Petersburg, Russian Federation
The paper interrelates changes in the effective elastic and diffusion properties of a two-phase composite using microstructural parameters. It is suggested that there are some inhomogeneities identical in shape in the material. The development of the cross-property connections in the explicit tensor form has been presented. The segregation effect, being a constant jump in concentration of particles of the solute flux at the matrix/inhomogeneity interface, was taken into account. It is a good practice to apply the derived cross-property relations to finding some effective properties of material using others when the material's microstructure is unknown. The obtained expressions were put to the test for isotropic material with pores; the approximate correlations were compared with exact ones found for the specific microstructure.
Keywords: effective Young's modulus, effective diffusion coefficient, cross-property connection, segregation effect
Citation: Frolova K.P., Cross-property connections between Young's modulus and diffusion coefficient of two-phase composite, St. Petersburg Polytechnical State University Journal. Physics and Mathematics. 14 (1) (2021) 177-189. DOI: 10.18721/JPM.14114
This is an open access article under the CC BY-NC 4.0 license (https://creativecommons.org/ licenses/by-nc/4.0/)
Введение
Нахождение соотношений между различными эффективными свойствами гетерогенных материалов является фундаментальной проблемой механики [1, 2]. Данные математические выражения устанавливают зависимости между изменениями физических свойств, различных по природе, но связанных с одними и теми же особенностями микроструктуры материала. Значимость данной проблемы обусловлена тем, что дает возможность определять одни эффективные свойства через другие при отсутствии полной информации о микроструктуре композита.
Поиск взаимосвязей между изменениями различных свойств обсуждается в литературе, начиная с 1960-х годов. Детальный обзор данной проблемы представлен в монографии [3]; согласно представленной информации, опубликованные исследования можно условно разделить на четыре направления: качественные изыскания, установление эмпирических зависимостей, определение диапазонов изменения характеристик материалов,
нахождение соотношений для материалов с изолированными неоднородностями в явном виде.
В настоящее время имеется лишь ограниченное количество работ, опубликованных по четвертому направлению.
Взаимосвязь эффективных упругих и проводящих свойств материалов в явном виде была установлена в работах [4, 5]. Полученные авторами выражения можно использовать для описания в общем случае анизотропных материалов с изотропной матрицей; при этом показано, что точность предлагаемых соотношений зависит от формы неоднородностей и разницы в упругих свойствах компонентов среды. В статье [6] получены связующие уравнения между
прочностными характеристиками материала и теплопроводностью металлического композита с чешуйками графита. В статьях [7, 8] получены соотношения между эффективным тепловым расширением и эффективной теплопроводностью, точность которых зависит от формы неоднородностей и разницы в теплопроводности фаз композита. В статьях [9, 10] получены взаимные зависимости между эффективным коэффициентом диффузии и теплопроводностью металл-алмазного композита. В статье [11] в явном виде получены соотношения между тепловыми и электропроводящими свойствами композита.
Процесс диффузии существенно отличается от процесса теплопроводности, рассмотренного в статьях [4, 5], тем, что концентрация, как правило, испытывает скачок на границе раздела двух фаз (матрица — неоднородность), тогда как температура является непрерывной функцией [12]. Явление накопления частиц диффундирующего вещества на границе раздела фаз или внутри включений известно как эффект сегрегации [13].
Данная работа посвящена установлению взаимосвязи между эффективными упругими и диффузионными свойствами двухфазного композита с неоднородностями одинаковой формы; при этом учитывается эффект сегрегации.
Полученные соотношения проверяются для случая изотропного материала с порами; при этом определяется связь между изменениями модуля Юнга и коэффициента диффузии.
Соотношения между эффективными упругими и диффузионными свойствами двухфазного материала
Соотношения между эффективными упругими и проводящими свойствами, полученные в статьях [4, 5], основаны на том,
что изменение этих различных по природе свойств контролируется одними и теми же микроструктурными параметрами. При выводе соотношений использовалось два основных предположения:
неоднородности имеют форму сфероидов;
эффективные свойства определяются без учета взаимодействия неоднородностей.
При этом в статьях [4, 5] показано, что полученные соотношения выполняются (в том числе) для материалов с неоднородностями, форма которых отлична от сфероидальной, и при больших концентрациях, чем допускает метод гомогенизации без учета взаимодействия неоднородностей. Это объясняется тем, что форма и концентрация неоднород-ностей влияют на упругие и проводящие свойства в одинаковой степени.
Основываясь на аналогии между уравнениями диффузии и теплопроводности, воспользуемся для определения связи между эффективными упругими и диффузионными свойствами методикой, изложенной в работах [4, 5]. Согласно данной методике, в рассмотрение вводятся тензоры вклада в искомые свойства, которые выступают как основные микроструктурные параметры [3].
При определении эффективных упругих свойств вводится тензор вклада неоднородности в податливость — тензор 4-го ранга Н, описывающий дополнительную деформацию, возникающую в представительном объеме ввиду присутствия изолированной неоднородности. Указанный тензор зависит от формы неоднородности и разности упругих свойств матрицы и неоднородности. В случае сфероидального включения тензор 4-го ранга Н является трансверсально-изотропным (ось симметрии совпадает с осью симметрии сфероида) и может быть представлен как линейная комбинация элементов тензорного базиса Тр Т2,..., Т6:
н = ^КТ,
к=1
Т = ее,
Т2 = 2 ((ееУт+(ее)Гт-ее) ,
т3 = епп, т4 = ппе,
Т5 = 4 (пеп + ( пеп )(,
+(епп )(1,4 )+(епп);
2)(3,4)
2,3)
Т6 = пппп,
(0 — проектор на плоскость, ортогональную орту оси симметрии сфероидальной неоднородности п; 0 = I — пп (I — единичный тензор 2-го ранга)).
Эффективный тензор податливости можно определить в рамках метода гомогенизации без учета взаимодействия неоднородно-стей следующим образом:
8е# = +1 ^ Н к,
V'
(2)
где 80 — тензор податливости матрицы, V — представительный объем, Vk — объем к-ой неоднородности.
При определении эффективных упругих свойств материала с одинаковыми по форме неоднородностями, выражение (2) можно переписать в виде
8е# = 8° +р(^П + ЩI) + + Щ (I ю+ю I) + Щ (I • ю + ю1)+ (3)
+ щ п
где
I = 2 ((II )(1,4) + (II )Г
(2,4)
(1)
где
— единичный тензор 4-го ранга; параметры Щ ( = 1, 2, ..., 5) выражаются через входящие в выражение (1) коэффициенты к. как
Щ = к,-к2/2, Щ = к2, Щ = - 2К + к2 + 2к3,
Ж4 = - 2И2 + И5, И
ж5 = А + у - 2И3 - И5 + И6,
П = VТК (пппп)(к>, V к
ю = VТК (пп)(к',
р==1 Тк
(4)
V
к
р — объемная доля неоднородностей.
Аналогичным образом определяются эффективные диффузионные свойства материала, состоящего из матрицы с шаровым тензором диффузии Б0 = А01 и неоднородностей с Б1 = А^ [14]. В рассмотрение вводится либо тензор вклада неоднородности в диффузию — НА (тензор 2-го ранга), определяющий дополнительный массовый поток, вызванный присутствием в материале неоднородности, либо тензор вклада в сопротивляемость НАК (тензор 2-го ранга), Н БК = - Н А/А02.
Предполагается, что и матрица, и неоднородность подчиняются линейному закону Фика. При этом считается, что нормальная компонента потока диффундирующего вещества непрерывна при переходе через границу матрицы (со знаком плюс) и неоднородности (со знаком минус), а концентрация испытывает скачок
х )| = 5С ( х )|
>дК- ■
тором частицы диффундирующего вещества скапливаются на границе раздела двух фаз, тогда как случай 5 < 1 соответствует оседанию частиц внутри неоднородностей [14]. При исследовании матрицы с порами, только случай 5 < 1 представляет интерес с физической точки зрения.
В случае сфероидальной неоднородности тензоры вклада определяются следующим образом:
НА = А0 [В1 (I - пп ) + В2пп
Н
БК
~[В (I - пп) + В2пп],
(6)
при этом коэффициенты В1 и В2 зависят от формы неоднородности, разницы в коэффициентах диффузии матрицы и неоднородности, а также от параметра сегрегации.
Эффективный тензор диффузии вводится в рамках метода гомогенизации без учета взаимодействия неоднородностей как
^ = АД+V ТКк Н А, К к
(7)
(5)
где 5 — параметр сегрегации.
Наличие эффекта сегрегации приводит к принципиальному отличию процесса диффузии от процесса теплопроводности, рассмотренного в статьях [4, 5]. Параметр сегрегации показывает соотношение концентраций на границе неоднородности и внутри нее и равен единице в случае непрерывной функции концентрации при переходе через границу раздела матрица/неоднородность. Случай 5 > 1 соответствует материалу, в ко-
а эффективный тензор сопротивляемости имеет вид
(^ )-1=а I+V тКнк . (8) А0 К к
Эффективные диффузионные свойства материала с одинаковыми по форме неодно-родностями можно выразить через объемную долю неоднородностей р и тензор второго ранга ю, определяемые выражениями (4), следующим образом:
—-1 = I - А0 ()-1 А0
= В1р1 +(В2 -В1)ю.
(9)
Установление явной взаимосвязи между эффективными упругими и диффузионными свойствами возможно, если они выражаются через одни и те же микроструктурные параметры, т. е. скалярный параметр р и тензор-
а
ный параметр ю. Таким образом, для получения соотношений между эффективными свойствами, в выражении для эффективного тензора податливости (3) необходимо избавиться от слагаемого, содержащего тензор 4-го ранга П. Согласно утверждениям в статьях [4, 5], это возможно за счет корректировки коэффициентов Щ, Щ2, Щ, Щ4:
8е// = 8° + -1 [р (^ • I + ,21 )
+^ ( Гю + юТ ) + ^ (1ю+ю1)
где Е0 — модуль Юнга матрицы, коэффициенты я. (. = 1, 2, 3, 4) определяются как
¿1 = Е° (К - к2/2), ¿2 = Е°К, ¿3 = Е° (-2/ + к + 2/),
¿4 = Е° ( 2к2 + К5 ),
а коэффициенты /г. (. = 1, 2, ..., 5) выражаются через к. следующим образом:
/г. = к (1 - 5 ¿щп/ )
при I = 1, 2, 6; /г. = / (1 + 5 в1§пкг.) при I = 3, 5; 5 =
К + к2 / 2 - 2к3 - к5 + к6
/ + к!/2 + 2 |к3| + |к5| + |к6|
Для вывода соотношений между эффективными упругими и диффузионными свойствами материала для общего случая распределения неоднородностей по ориентациям, выразим тензор ю и параметр р через тензоры и (Бе/ ) , соответственно:
ю =
' 1
(^2 - В)
V Ц
Б
е//
В1
(В2 - В)
\е//
р!,
(11)
Р =
1гБе/ - 3Ц° Ц (2ВХ + В°) :
ю =
(1°)
(В2 - В!)
I - А (Бе//)-
В
(12)
(В2 - В!)
Р^
Р =
3 - Б°гг (Бе/
2В] + В2
В результате подстановки этих выражений в равенство (10) можно получить два соотношения, устанавливающих явную взаимосвязь между эффективной податливостью и эффективными диффузионными свойствами:
Е° (8е// - 8° ) = = (а^ • I + а21) х
х( 1гБе# / Ц° - 3) + +а3 [(Бе// / Ц° -1)I + (13)
+Х (Бе/ / Ц° -1)] + +а4 [(Бе/ / Ц° -1)-1 + +1-(Бе// / Ц° -1)] при использовании соотношений (11) и
Е° (8е//- 8° ) = = (ахХ •Х + а21 )х
х(3 - Ц°1х ( Бе//)-1) +
1
+а„
+а,
(I - А0 ( В* )-1 +1 (I - А ( В* )-1) (I - А ( В* )-1 )•. (I - А0 ( В* )-1)
(14)
+
I +
+ЛI -
при использовании соотношений (12); здесь
_ ( В2 - В!)-5з В!
1 (В2 -Вх)(2ВХ + В2У
а =
а 2 =
*2 (В2 - В,)-54 В (В2 - В,)(2В + В2 у
аз =
а 4 =
2 (В2 - В,)' 2 ( В2 - В, )■
Соотношения между эффективным модулем Юнга и эффективным коэффициентом диффузии материала с порами
Для определения взаимосвязи между эффективными упругими модулями и эффективными коэффициентами диффузии необходимо использовать компонентное представление выражений (13) и (14). В случае изотропного материала будут иметь место следующие соотношения:
Е
Еп
■ [1 + (3а, + 3а 2
+2а3 + 2а 4)
(Ае1 ^
V А0
(13) и
Еп
■ [1 + (3а, + 3а 2
(16)
+2а3 + 2а 4)
(15)
при использовании тензорного выражения
при использовании тензорного выражения
(14).
Соотношения (15) и (16) выражают зависимость одного эффективного модуля от другого на макроуровне. Необходимая информация о микроструктуре сводится к определению формы отдельной неоднородности, поскольку от нее зависят коэффициенты а.. В общем случае соотношения (15) и (16) являются приближенными, поскольку они получены с использованием приближенного выражения для эффективного тензора податливости (10).
Изотропному материалу, содержащему неоднородности, соответствуют два случая: они имеют сферическую форму, они распределены по ориентациям произвольным образом.
Соотношения (15) и (16) являются точными в первом случае, поскольку 5 = 0 и, следовательно, приближенное выражение для эффективного тензора податливости (10) совпадает с точным. Во втором случае 5 Ф 0 и выражения (15), (16) являются приближенными.
Исследуем два этих случая и оценим точность полученных приближенных соотношений.
Определим по формулам (15), (16) связь между эффективным модулем Юнга и эффективным коэффициентом диффузии материала со сферическими порами и сравним полученное выражение с прямо вычисляемым.
В случае матрицы с порами коэффициенты тензора вклада неоднородности в податливость И имеют вид
5
3
Цб 2Д' К = 1 Ц2' к 4з Д'
- д± Д , К 4 45 , Кб II Д|2
В =
1 -к
(17)
где А = 2Ц д6 — ц) и
д = |[4к-1 -2(3к-1)/ -2к/ ], д2 = 2|[1 -(2-к) /-к/ ], 4з = д, = 2|[( 2к-1) / + 2к/ ],
45 = 4|[ / + 4к/ ],
46 = 81к[ ./с - >1 ]> к = (1-у)/2;
ц — модуль сдвига матрицы; V — коэффициент Пуассона матрицы; /, / — функции, зависящие от формы сфероидальной неоднородности, т. е. от отношения длин полуосей сфероида у = а3/а (а — полуось вращения) как
/о =
1 - 8
2 (1 -Г-2)'
/ = [(2+у-2) 8-3у
1 - 2 (1 - ¿к) /с
(18)
где А = Б0/Бг
Тогда выражения (15) и (16) сведутся соответственно к равенствам
Е
Еп
1, (1 -у)(9 + 5У .
2(7 - 5у)
1 + 2sк
1 -к
Б -1
V Бс У
Е
Е
1 , (1 -у)(9 + 5У) .
2(7 - 5У)
(19)
1 + 2sк
1-к
1 - а
Б
где
Чтобы получить прямую взаимосвязь между эффективными модулями, выразим объемную долю неоднородностей непосредственно через эффективный коэффициент диффузии, определяемый с помощью метода гомогенизации без учета взаимодействия пор, и подставим в точное выражение для эффективного тензора податливости, определяемого также без учета взаимодействия.
Так, тензоры вклада сферической неоднородности в диффузию и в сопротивляемость определяются выражениями
8 =
,2
1
1 л/1-Г , агСап --, у < 1
-1п
У
у+УУ2-1
, у> 1.
Коэффициенты тензора вклада в диффузию имеют вид
В =
1 -к
¿к + (1 -¿к) /0'
НБ = « ОД
Н
БК
2 ¿к + 1
3 (1 -к)
( 2^к + 1) Бс
(20)
I.
(18)
Подстановка первого и второго выражений (20) соответственно в уравнение для эффективного тензора диффузии (7) и эффективного тензора сопротивляемости (8) позволит выразить объемную долю неодно-родностей через эффективные диффузионные свойства.
Последующая подстановка полученного выражения для р в уравнение для эффективного тензора податливости (2), с учетом формулы для изотропного тензора вклада сферической поры в податливость, т. е.
Н
15 (1 -V)
2| 1
1
1
7 - 5^
10 (1 + V) 3
.-11 • 1'
I • I +
зывающее эффективный модуль Юнга с эффективным коэффициентом диффузии:
Е
Л
Еп
1 +
Ео^ (
15п
Л"
V По
-1
У.
(22)
(21)
а при использовании этого же метода в терминах тензоров вклада в диффузию выражение имеет вид:
Е
Л
приведет к результату, в точности совпадающему с выражениями (19).
Таким образом, выполнена проверка полученных приближенных соотношений для случая микроструктуры, при которой они являются точными.
Во втором случае, соответствующем изотропному материалу (при произвольном распределении неоднородностей по ориен-тациям), соотношения (15) и (16) являются приближенными. При этом точную взаимосвязь между эффективными упругими и диффузионными свойствами можно установить для случая произвольного распределения сфероидов по ориентациям независимо, по аналогии с тем, как это делалось выше для случая сферических неоднородностей.
Для количественной и качественной оценок полученных нами приближенных соотношений, сравним их с точными, определенными для конкретной микроструктуры.
Для установления точных соотношений определим объемную долю неоднородностей через эффективный коэффициент диффузии изотропного материала и затем подставим полученное выражение в уравнение для определения эффективного модуля Юнга. Адекватность такого подхода обусловлена изотропией эффективных тензоров, в результате которой единственным общим микроструктурным параметром, определяющим эффективные свойства материала, является скалярный параметр р. При использовании метода гомогенизации без учета взаимодействия в терминах тензоров вклада в диффузию, получается следующее выражение, свя-
Е
(23)
где = 8^ + 4И2 + 4И3 + 2к5 + 3к6, П = 2Я/3 + В2/3.
Сравнение результатов, полученных с помощью приближенных соотношений (15), (16), с результатами, полученными на основании точных соотношений (22) и (23) соответственно, представлены на рис. 1 для случая X = 0, что соответствует П1 ^ ю (в этом случае параметр сегрегации, как следует из соотношения (18), не играет роли), Е0 = 208 ГПа.
Видно, что соотношения (15) и (16) достаточно точно описывают связь между эффективными свойствами на макроуровне. При этом эффективные значения, полученные через тензоры вклада в сопротивляемость диффузии (рис. 1, а) и через тензоры вклада в диффузию (рис. 1, Ь), различаются. Это связано с тем, что соотношения получены без учета взаимовлияния неоднородностей.
Для получения точных соотношений с учетом взаимовлияния неоднородностей обратимся к методу гомогенизации Максвелла [3]. Результат, полученный с помощью этого метода через тензоры вклада в диффузию, совпадает с результатом, полученным через тензоры вклада в сопротивляемость и, таким образом, результат однозначен, в отличие от метода гомогенизации без учета взаимодействия.
Модуль Юнга материала со сфероидальными порами определялся нами в статье [15], а эффективный коэффициент диффузии такого материала — в статье [16]. Для получения точных соотношений необходимо
Е&!Е<> 1.0
0.9
0.8
0.7
0.6
а)
У = 5
У — I
1.0 1.2 1.4 1.6 1.8 2.0
2.2 2.4
Е^/Ео 1.0
0.9 0.8 0.7 0.6 0.5 0.4
Ъ)
■1 ■ 1
Г =
- у = хТ
у= 1
~~! ! ! 1 Г~
1.0 1.2 1.4 1.6 1.8 2.0
2.2 2.4
Рис. 1. Зависимости относительного эффективного модуля Юнга от относительного эффективного коэффициента диффузии при использовании соотношений податливость-сопротивляемость диффузии (а) и податливость-диффузия (Ь) для различных значений отношения длин полуосей сфероида у.
Сравниваются приближенные соотношения (сплошные линии) с точными (пунктиры)
подставить выражение для объемной доли включений
3 Б/ - а
(Бе// + 2Бс )п
в выражение для эффективного модуля Юнга.
Сравнение установленных в данной работе приближенных соотношений с точными соотношениями, полученными в рамках схемы Максвелла, представлено на рис. 2 для случая А = 0, Е0 = 208 ГПа.
Анализ графиков на рис. 2 позволяет заключить, что предпочтительнее использовать соотношения, построенные с помощью тензоров вклада в диффузию, поскольку они лучше совпадают с точным соотношением, построенным с использованием схемы Максвелла.
Исследуя влияние сегрегации на взаимосвязь между эффективными свойствами материала, учтем, что коэффициент диффузии пор намного больше коэффициента диффузии матрицы (Б1 >> Б0), но при этом не бесконечен; такое предположение соответствует реальному материалу. Примем для определенности, что А = 0,01, Е0 = 208 ГПа. Зависимость от параметра сегрегации неявно присутствует в выражениях (15), (16) только в коэффициенте п, поскольку
3ах + 3а 2 + 2а3 + 2а4 = 3 (¿1 + ¿2) + ¿3 + ¿4 3п
этот же коэффициент можно выделить и в самом выражении для эффективного тензора диффузии.
Зависимость коэффициента п от параметра сегрегации при разных значениях отношения длин полуосей сфероида показана на рис. 3.
Видно, что сегрегация не оказывает заметного эффекта при наличии сферических пор. В случае сплюснутых и вытянутых сфероидальных пор влияние сегрегации более выражено, что сказывается и на эффективных коэффициентах диффузии [14, 16]. При этом, согласно зависимостям, изображенным на рис. 4, сегрегация не оказывает заметного влияния непосредственно на связь между эффективными свойствами (графики построены для соотношений, основанных на связи податливость-диффузия). Таким образом, еще одним преимуществом использования установленной взаимосвязи между эффективными свойствами в явном виде является ее независимость от того, оседают частицы диффундирующего вещества в порах или нет — достаточно информации об
Е&/Еа 1.0 0.9 0.8 0.7 0.6 0.5 0.4
1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4
Е^/Па
Рис. 2. Сравнение приближенных соотношений податливость-диффузия (сплошные линии) и податливость-сопротивляемость (точечные линии) с точными, полученными в рамках схемы Максвелла (пунктирные линии)
у = 0.1
7 = 5
7=1
0.0 0.2 0.4 0.6 0.8 1.0
£
Рис. 3. Зависимость коэффициента п от параметра сегрегации при разных значениях отношения длин полуосей сфероида у
Е&/Ео
■1 1 1 1 ■
г= 5
; У = = 0.1
1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4
Е&ИУъ
Рис. 4. Зависимости относительного эффективного модуля Юнга от относительного эффективного коэффициента диффузии при двух значениях отношения длин полуосей сфероида у и параметра сегрегации з: 1,00 (сплошная линия) и 0,01 (штрихпунктир)
эффективных свойствах материала (упругих или диффузионных).
Заключение
В работе в явном виде получены соотношения, устанавливающие связь между эффективными упругими и диффузионными свойствами материала. Полученные зависимости целесообразно использовать для определения изменений одних свойств через другие, когда неизвестна микроструктура материала.
При выводе соотношений использовалась известная (из литературы) методика, изложенная в терминах тензоров вклада для определения взаимосвязи между эффективными упругими и теплопроводящими свойствами. В полученных в настоящей статье зависимостях учтен эффект сегрегации, принципиально отличающий процесс диффузии от процесса теплопроводности.
Подробно исследован случай изотропного материала с порами, в которых может оседать диффундирующее вещество. Показано, что для этого случая установленные соотношения достаточно точно описывают
зависимость эффективного упругого модуля Юнга от эффективного коэффициента диффузии.
На основании сравнения указанных результатов с полученными в рамках метода Максвелла, который учитывает взаимовлияние множественных неоднородностей, сделан вывод о предпочтительном применении выражений, связывающих податливость материала с диффузионной проводимостью, а не с сопротивляемостью.
Проанализировано влияние сегрегации на взаимосвязь между эффективными свойствами материала. Показано, что учет сегрегации через задание соответствующего параметра не оказывает в данном случае заметного влияния, что доказывает преимущество использования полученных соотношений при определении одних свойств через другие на макроуровне, поскольку отпадает необходимость эмпирического определения величины параметра сегрегации.
Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (грант № 20-08-01100).
СПИСОК ЛИТЕРАТУРЫ
1. Berryman J.G., Milton G.W. Microgeometry of random composites and porous media // Journal of Physics D: Applied Physics. 1988. Vol. 21. No. 1. Pp. 87-94.
2. Gibiansky L.V., Torquato S. Connection between the conductivity and bulk modulus of iso-tropic composite materials // Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences. 1996. Vol. 452. No. 1945. Pp. 253-283.
3. Kachanov M., Sevostianov I. Micromechanics of materials, with applications. Berlin, Germany: Springer, 2018. Vol. 249. 712 p.
4. Sevostianov I., Kachanov M. Explicit cross-property correlations for anisotropic two-phase composite materials // Journal of the Mechanics and Physics of Solids. 2002. Vol. 50. No. 2. Pp. 253-282.
5. Sevostianov I., Kachanov M. Connections between elastic and conductive properties of het-
erogeneous materials // Advances in Applied Mechanics. Vol. 42. Edited by H. Aref, E. van der Giessen. USA: Elsevir, 2009. Pp. 69-252.
6. Mazloum A., Oddone V., Reich S., Sevostianov I. Connection between strength and thermal conductivity of metal matrix composites with uniform distribution of graphite flakes // International Journal of Engineering Science. 2019. Vol. 139. June. Pp. 70-82.
7. Sevostianov I. On the thermal expansion of composite materials and cross-property connection between thermal expansion and thermal conductivity // Mechanics of Materials. 2012. Vol. 45. February. Pp. 20-33.
8. Mazloum A., Sevostianov I. Connections between anisotropic tensors of thermal conductivity and thermal expansion coefficients // International Journal of Engineering Science. 2018. Vol. 122. January. Pp. 1-13.
9. Anisimova M.A., Knyazeva A.G., Sevostia-
nov I. Connection between diffusion coefficient and thermal conductivity of a metal matrix composite // IOP Conference Series: Materials Science and Engineering. 2017. Vol. 175. No. 1. P. 012051.
10. Anisimova M., Knyazeva A., Sevostianov I.
Effective thermal properties of an aluminum matrix composite with coated diamond inhomogeneities // International Journal of Engineering Science. 2016. Vol. 106. September. Pp. 142-154.
11. Mazloum A., Kovacik J., Emmer S., Sevostianov I. Copper-graphite composites: Thermal expansion, thermal and electrical conductivities, and cross-property connections // Journal of Materials Science. 2016. Vol. 51. No. 17. Pp. 7977-7990.
12. Zhang Y., Liping L. On diffusion in heterogeneous media // American Journal of Science. 2012. Vol. 312. No. 9. Pp. 1028-1047.
13. Kaur I., Mishin Y., Gust W. Fundamentals of grain and interphase boundary diffusion. 3rd re-
vised and enlarged edition. John Wiley & Sons. Inc., 1995. 528 p.
14. Knyazeva A.G., Grabovetskaya G.P., Mishin I.P., Sevostianov I. On the microme-chanical modelling of the effective diffusion coefficient of a polycrystalline material // Philosophical Magazin. 2015. Vol. 95. No. 19. Pp. 2046-2066.
15. Фролова К.П. Определение эффективного модуля Юнга среды с микроструктурой, характерной для водородной деградации // Научно-технические ведомости СПбГПУ. Физико-математические науки. 2020. Т. 13. № 2. C. 160-174.
16. Frolova K.P., Vilchevskaya E.N. Effective diffusion coefficient of a porous material applied to the problem of hydrogen damage (Chapter 7) // Polyanskiy V.A., Belyaev A.K. (Eds.). Advances in Hydrogen Embrittlement Study. Springer, Cham., 2021. Advanced Structured Materials. Vol. 143. Pp. 113-130.
Статья поступила в редакцию 19.01.2021, принята к публикации 01.02.2021.
СВЕДЕНИЯ ОБ АВТОРЕ
ФРОЛОВА Ксения Петровна — младший научный сотрудник Института проблем машиноведения РАН, Санкт-Петербург, Российская Федерация; ассистент Высшей школы теоретической механики Санкт-Петербургского политехнического университета Петра Великого, Санкт-Петербург, Российская Федерация.
199178, Российская Федерация, г. Санкт-Петербург, Большой проспект В.О., 61 [email protected]
REFERENCES
1. Berryman J.G., Milton G.W., Microgeometry of random composites and porous media, Journal of Physics D: Applied Physics. 21(1) (1988) 87-94.
2. Gibiansky L.V., Torquato S., Connection between the conductivity and bulk modulus of isotropic composite materials, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences. 452 (1945) (1996) 253-283.
3. Kachanov M., Sevostianov I., Microme-chanics of materials, with applications, Vol. 249, Springer, Berlin, Germany, 2018.
4. Sevostianov I., Kachanov M., Explicit cross-property correlations for anisotropic two-phase composite materials, Journal of the Mechanics and Physics of Solids. 50 (2) (2002) 253-282.
5. Sevostianov I., Kachanov M., Connections between elastic and conductive properties of heterogeneous materials, In the book: Advances in Applied Mechanics, Edited by H. Aref, E. van der Giessen, Elsevir, USA. 42 (2009) 69-252.
6. Mazloum A., Oddone V., Reich S., Sevostianov I., Connection between strength and thermal conductivity of metal matrix composites
with uniform distribution of graphite flakes, International Journal of Engineering Science. 139 (June) (2019) 70-82.
7. Sevostianov I., On the thermal expansion of composite materials and cross-property connection between thermal expansion and thermal conductivity, Mechanics of Materials. 45 (February) (2012) 20-33.
8. Mazloum A., Sevostianov I., Connections between anisotropic tensors of thermal conductivity and thermal expansion coefficients, International Journal of Engineering Science. 122 (January) (2018) 1-13.
9. Anisimova M.A., Knyazeva A.G., Sevostia-nov I., Connection between diffusion coefficient and thermal conductivity of a metal matrix composite, IOP Conference Series: Materials Science and Engineering. 175 (1) (2017) 012051.
10. Anisimova M., Knyazeva A., Sevostianov I., Effective thermal properties of an aluminum matrix composite with coated diamond inhomogeneities, International Journal of Engineering Science. 106 (September) (2016) 142-154.
11. Mazloum A., Kovacik J., Emmer S., Sevostianov I., Copper-graphite composites: Thermal expansion, thermal and electrical conductiv-
Received 19.01.2021, accepted 01.02.2021.
ities, and cross-property connections, Journal of Materials Science. 51 (17) (2016) 7977-7990.
12. Zhang Y., Liping L., On diffusion in heterogeneous media, American Journal of Science. 312 (9) (2012) 1028-1047.
13. Kaur I., Mishin Y., Gust W., Fundamentals of grain and interphase boundary diffusion, 3rd revised and enlarged edition, John Wiley & Sons. Inc., 1995.
14. Knyazeva A.G., Grabovetskaya G.P., Mishin I.P., Sevostianov I., On the micromechanical modelling of the effective diffusion coefficient of a poly-crystalline material, Philosophical Magazin. 95 (19) (2015) 2046-2066.
15. Frolova K.P., Determination of the effective Young's modulus of medium with microstructure typical for hydrogen degradation, St. Petersburg Polytechnical State University Journal. Physics and Mathematics. 13 (2) (2020) 160-174.
16. Frolova K.P., Vilchevskaya E.N., Effective diffusion coefficient of a porous material applied to the problem of hydrogen damage, Ch. 7, In the book: Polyansky V.A., Belyaev A.K. (Eds.). Advances in Hydrogen Embrittlement Study, Springer, Cham., 2021, Advanced Structured Materials, Vol. 143, Pp. 113-130.
THE AUTHOR
FROLOVA Ksenia P.
Institute for Problems in Mechanical Engineering RAS
Peter the Great St. Petersburg Polytechnic University
61 Bolshoi Ave. of V. Isl., St. Petersburg, 199178, Russian Federation
© Санкт-Петербургский политехнический университет Петра Великого, 2021
