К ТЕОРИИ УПРУГОСТИ МИКРОНЕОДНОРОДНЫХ СРЕД, УЧИТЫВАЮЩЕЙ СТОХАСТИЧЕСКИЕ ИЗМЕНЕНИЯ СВЯЗНОСТИ СОСТАВЛЯЮЩИХ КОМПОНЕНТОВ

Сараев Л.А. К теории упругости микронеоднородных сред, учитывающей стохастические изменения связности составляющих компонентов // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2021. - № 2. С. 132-143. БО!: 10.15593/регш.шесЬ/2021.2.12

Saraev L.A. On the theory of elasticity of micro-inhomogeneous media with account for stochastic changes in the connectivity of constituent components. PNRPUMechanics Bulletin, 2021, no. 2, pp. 132-143. DOI: 10.15593/perm.mech/2021.2.12

ВЕСТНИК ПНИПУ. МЕХАНИКА № 2, 2021 PNRPU MECHANICS BULLETIN

https://ered.pstu.ru/index.php/mechanics/index

БО1: 10.15593/регш.шесЬ/2021.2.12 УДК 539.3

К ТЕОРИИ УПРУГОСТИ МИКРОНЕОДНОРОДНЫХ СРЕД, УЧИТЫВАЮЩЕЙ СТОХАСТИЧЕСКИЕ ИЗМЕНЕНИЯ СВЯЗНОСТИ СОСТАВЛЯЮЩИХ КОМПОНЕНТОВ

Л.А. Сараев

Самарский национальный исследовательский университет им. академика С.П. Королева, Самара, Россия

О СТАТЬЕ АННОТАЦИЯ

Предложена математическая модель расчета эффективных модулей упругости микронеоднородного двухкомпонентного изотропного композиционного материала, связность компонентов которого изменяется случайным образом в зависимости от уровня их относительных объемных содержаний. Сформулировано стохастическое уравнение для параметра связности составляющих компонентов, согласно которому при увеличении объемного содержания наполнителя отдельные включения образуют структуры матричной смеси в виде взаимопроникающих каркасов, а затем превращаются в новую связующую матрицу с отдельными включениями из материала остатков старой матрицы. Алгоритм численного решения этого стохастического дифференциального уравнения построен в соответствии с методом Эйлера - Маруямы. Для каждой реализации этого алгоритма строятся соответствующие стохастические траектории для случайной функции связности составляющих компонентов композиционного материала. Разработан вариант метода расчета математического ожидания случайной функции связности составляющих компонентов и получено для него соответствующие дифференциальное уравнение. Показано, что численное решение этого уравнения и среднее значение функции фактора производства, вычисленное по всем реализациям стохастических траекторий, дают близкие численные значения. Установлены новые макроскопические определяющие соотношения для микронеоднородных материалов с изменяемой микроструктурой и вычислены их эффективные модули упругости. Отмечено, что формулы для этих эффективных модулей упругости обобщают известные результаты для изотропных композиционных материалов. Значения эффективных модулей упругости, построенных по полученным в статье выражениям, лежат внутри вилки Хашина - Штрикмана для нижней и верхней границы эффективных модулей упругости композиционных материалов. Численный анализ разработанных моделей показал хорошее соответствие известным экспериментальным данным.

©ПНИПУ

Получена: 8 февраля 2021 г. Принята: 31 мая 2021 г. Опубликована: 12 июля 2021 г.

Ключевые слова:

включения, композит, компоненты, макроскопические свойства, матрица, матричная смесь, модули упругости, связность, статистическая однородность, микроструктура, эргодичность, эффективные соотношения.

© Сараев Леонид Александрович - д.ф.-м.н., проф., зав. каф., e-mail: saraev_leo@mail.ru, : 0000-0003-3625-5921.

Leonid A. Saraev - Doctor of Physical and Mathematical Sciences, Professor, e-mail: saraev_leo@mail.ru, : 0000-0003-3625-5921.

Эта статья доступна в соответствии с условиями лицензии Creative Commons Attribution-NonCommercial 4.0 International License

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

ON THE THEORY OF ELASTICITY OF MICROINHOMOGENEOUS MEDIA WITH ACCOUNT FOR STOCHASTIC CHANGES IN THE CONNECTIVITY OF CONSTITUENT COMPONENTS

L.A. Saraev

Samara National Research University, Samara, Russian Federation

ABSTRACT

The paper proposes a mathematical model aimed at calculating the effective elastic moduli of a micro-inhomogeneous two-component isotropic composite material, which components are connected randomly depending on the level of their relative volumetric contents. A stochastic equation is formulated for the connectivity parameter of the constituent components, according to which, with an increase in the volumetric content of the filler, individual inclusions build the structures of the matrix mixture in the form of interpenetrating frameworks, and then turn into a new binding matrix with individual inclusions from the material of the rest of the old matrix. The algorithm for the numerical solution of this stochastic differential equation is constructed in accordance with the Euler-Maruyama method. For each implementation of this algorithm, the corresponding stochastic trajectories are constructed for the random connectivity function of the constituent components of the composite material. A variant of the method aimed at calculating the mathematical expectation of a random connectivity function of the constituent components has been developed and the corresponding differential equation has been obtained for it. It is shown that the numerical solution of this equation and the average value of the production factor function calculated for all realizations of stochastic trajectories give close numerical values. New macroscopic constitutive relations are found for microinhomogeneous materials with a variable microstructure and their effective elastic moduli are calculated. It is noted that the formulas for these effective elastic moduli generalize the known results for isotropic composite materials. The values of the effective elastic moduli, constructed according to the expressions obtained in the paper, lie within the Khashin-Shtrikman range for the lower and upper bounds of the effective elastic moduli of the composite materials. The numerical analysis of the developed models showed a good agreement with the known experimental data.

©PNRPU

ARTICLE INFO

Received: 8 February 2021 Accepted: 31 May 2021 Published: 12 July 2021

Keywords:

inclusions, composite, components, macroscopic properties, matrix, matrix mixture, elastic moduli, connectivity, statistical homogeneity, microstructure, ergodicity, effective relations.

Введение

Развитие методов расчета эффективных характеристик перспективных дисперсно-упрочненных композиционных материалов со случайной микроструктурой представляет собой актуальное направление современной механики деформируемого твердого тела. Подобные конструктивные материалы, наполненные высокопрочными и высокомодульными дисперсными частицами, структурно неоднородные порошковые стали и сплавы проявляют высокие механические и прочностные свойства, благодаря которым они широко применяются в машиностроении, в теплоэнергетическом комплексе, медицине и других отраслях.

Для расчетов эффективных характеристик упругих свойств композиционных материалов существует ряд классических, структурно-феноменологических методов и подходов, использующих модули упругости компонентов и их объемные содержания [1-9].

На основе этих методов были разработаны их модификации и построены модели композиционных материалов с различной степенью связности составляющих компонентов [10-19].

Эти модели достаточно адекватно описывают упругие свойства дисперсно-упрочненных композиционных материалов, содержащих хаотически распределенные

высокопрочные и высокомодульные частицы и нано-структурные элементы, пластиков, наполненных стеклянными микросферами и т.д. [20-26].

Кроме того, подобного рода модели способны давать достаточно точные оценки модулей упругости пористых металлов, металлокерамики, полимеров и строительных материалов [27-36].

Во всех отмеченных выше моделях дисперсные частицы или поры представляют собой отдельные включения, распределенные в связующей матрице. Значения эффективных модулей упругости зависят только от значений модулей упругости составляющих компонентов и их объемных содержаний. Экспериментальные исследования упругих свойств конструкционных композиционных материалов показывают, что при технологическом изготовлении образцов композитов путем прессования, брикетирования, спекания и пропитки расплавом связующего компонента, при одних и тех же значениях модулей упругости составляющих компонентов и их объемных содержаний, макроскопические эффективные модули упругости изготовленных образцов композитов практически всегда проявляют определенный разброс их свойств [37, 38].

Это связано с особенностями изготовления образцов композитов. Многочисленные исследования внутренней структуры дисперсно-упрочненных композитов

показывают, что при достаточно большом объемном содержании наполнителя отдельные включения, спекаясь, образуют структуры матричной смеси в виде взаимопроникающих каркасов. В предельных случаях дисперсные частицы могут образовать новую связующую матрицу с отдельными включениями из материала остатков исходной матрицы [39].

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

1. Постановка задачи

Пусть упругий изотропный двухкомпонентный композиционный материал занимает объем V, ограниченный поверхностью 5 . Законы Гука материалов его компонентов имеют вид

,

Здесь ,,, = с,.,. -1 • 5,, с рр; еЧ = еЧ. -1 • 5,Ч • е рр; ст.,., еЧ

Ч Ч ^ Ч рр ' Ч . 3 . Р .

- тензоры напряжений и деформаций, V - объемы составляющих компонентов

- радиус-вектор координат, ц - модули сдвига, К -объемные модули компонентов. Для определенности будем считать, что ц > ц и К2 > К .

Геометрические особенности структуры такого композита могут быть описаны случайными индикаторными функциями координат к (г), каждая из которых равна единице в точках объема V и равна нулю в точках вне этого объема. Для описания связности составляющих компонентов композиционного материала введем неопределенные модули упругости тела сравнения ц и К, удовлетворяющие неравенствам

ц <ц<ц; К <к<К [1].

С помощью функций к (г) и модулей ц и К закон Гука (1) запишется в виде

5ч

г))• е. (г) г ))•£ (г I

где [ц, ] = ц,-ц; [К, ] = К - К.

Все физические и геометрические поля рассматриваемой задачи предполагаются случайными, статистически однородными и эргодическими полями, и их математические ожидания заменяются средними значениями по соответствующим объемам [1-4]

(/> = \ I г (г) • &; (Д = У •//(г) • <ъ.

' V К

Здесь угловыми скобками обозначена операция статистического осреднения.

Для установления макроскопических определяющих уравнений и вычисления эффективных характеристик рассматриваемого композиционного материала необходимо найти связь между макронапряжениями

^с) и макродеформациями ^е^ , которая в самом

общем случае имеет вид

где Е*и - тензор эффективных модулей упругости композиционного материала.

Для установления макроскопических определяющих уравнений

(Ч = 2 ц • {еч) + 2 ([ц {еч) 1 +[ц2 ]• С2 • {еч) 2 ) , (с рр) = 3 ^ (е рр) + 3 ([ К1 ]• сг (е рр) 1 +[ ^ ]• V (е №) 2),

где с = ^ - объемные содержания составляющих

компонентов. Таким образом, для решения задачи необходимо выразить деформации ^ е;.) через макроскопические деформации ^ ) .

Получение соотношений (4) достигается статистическим осреднением системы уравнений, состоящей из закона Гука (2), уравнений равновесия

Ср, р

и формул Коши

2 е.

связывающих компоненты тензора деформаций (г) с компонентами вектора перемещений и (г). Граничные условия отсутствия флуктуаций величин на поверхности 5 объема V

/(г)| г^ =(/) .

2. Модули упругости композиционного материала со стохастической связностью составляющих компонентов

Исключение из системы уравнений (2), (5), (6) компонентов тензоров напряжений и деформаций дает систему равновесия микронеоднородной среды в перемещениях [5-7]

цКрр (г) + |К + £ • и'р,р,(г) -х,р,р (г) = 0

где

т.. = I. +---5.. -X ,

ч 9 3 ч рр

^ =-2-Ц-((от, - 1).К1 +(т2 - 1).к2)• еч,

Трр =-3-К •(( К - 1)-к, +(к2 - 1)-К2 )-е р

т = *, т = К = К = —,

1  2 -- 1 к 2 к

Ц

Ц

штрихами обозначены флуктуации величин в полном объеме V.

С помощью тензора Грина

С. (г ) =

1

8 -л-ц

3 • К + 5 • ц о • г--—• г

К ,рр 3• К + 8-ц д

содержащего неопределенные модули упругости тела сравнения ц и К , система уравнений (8), (9) сводится к системе интегральных уравнений, ядрами которых служат вторые производные тензора Грина [8]

еЧ (г) = | Олч (г - г,) Ч (г,) • й г,.

Величины находятся из известных соотноше-

ний [8, 9]

Ч =Ы + с-1

Для вычисления случайных моментов • е' ^ необходимо умножить уравнения (10) на к'(г), усреднить их по полному объему V и воспользоваться свойством изотропности функции к' (г) [8]

1 + а •

еч) 2 =

рр / 2

1+а •

. 1+у

"1 + У-((. -1) + с •(.2 -. ))•

где

2 4 - 5-V 1 1 +у 3 • К - 2 • ц

а =---; у =---; v = -

15 1 -V 3 1 -V' 2-(3 - К + 2-ц)'

Подставляя формулы (12) в выражения (4), получаем макроскопический закон для рассматриваемого композита

[Ы = 2 -Ц*е-ч),

К-рр> = 3-К* рр),

где

Ц = Ц-

с -

т +—-^--

1+а -

(

К * = К -

К +

^2 -(.2 - . )-(1 + У( . - 1)) '

1 + У-((- 1) + <1 -(.2 - )),

есть эффективные модули упругости рассматриваемого композита.

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

Если, например, в общих формулах (14) положить ц = ц , К = К , то получится модель композита, в котором первый компонент играет роль связующей матрицы, а второй компонент - роль отдельных включений. Выражения для эффективных модулей упругости (14) принимают вид

(

Цх = Ц -

К* = К -

1+-

1+

с2 -

1+а - с -

С -

Л

где

2 4 - 5-V- 1 1 + VI 1 3 - К- - 2-ц.

=---1; у, =---1; V. =---!-1

15 1 -V 3 1 -V 1 2 3 - К + 2-Ц

Ц2 К2

т = — ; я = .

Ц,

К

Если же в общих формулах (14) положить Ц = Ц2 , К = К , то получится модель композита, в котором, наоборот, второй компонент играет роль связующей матрицы, а первый компонент - роль отдельных включений. В этом случае выражения для эффективных модулей упругости (14) принимают вид

Цн =Ц2 1 +

с -

т + а2 - с2 -

Кн = К 2 '

1+-

с, -

д + У2- с2-

где

2 4 - 5-V,

1 1 + V

15 1 -V,

, У 2 =г

3 1 -V,

1 3-к2 -2-ц2

2 3 - К + 2-Ц2

Для допущенных ранее неравенств ц2 > ц, и К > К формулы (15) представляют собой нижние границы, а формулы (16) - верхние границы Хашина -Штрикмана для изотропных эффективных модулей упругости. Эти границы следуют из известных теорем теории упругости о минимуме потенциальной энергии и минимуме дополнительной энергии и являются существенно более узкими, чем границы Фойгта и Рейсса

\

/

е

j mr < M < vf , \KR < K* < Kf ,

где

Vl -Ц2

MF C1 ' V1 + C2 ' V2 ' VR

Cl -V2 + с2 -Ml kp = с - K1 + c2 - K2, KR =

K1 - K2

с - K2 + с - K

Другим типом связности составляющих компонентов обладает микронеоднородная среда, у которой матрица взаимопроникающими компонентами. Модули упругости такого композиционного материала получаются из общих формул (14), если в них положить

|(ц) = с1 -Ц, + с2 - Ц2 , КК) = с, - К, + С2 - К2.

В этом случае выражения для эффективных величин (14) принимают вид

vC = (v)-

-(V2 ~M-1 )2

KC = K -

(v) + a-(v2 -V1 )-(с1 -с2)'

У-( K2 - K1 )2 - C1 - с2

K + У-

где

а= 2_ 4 - 5 .. = 1 ^ . , 3 -

15 1 -

Из формул (19) для эффективных модулей упругости видно, что каждый компонент обладает самостоятельной несущей способностью и не представляет собой набор отдельных включений. Так, например, при ц = К = 0 макроскопические величины ц , К тождественно в нуль не обращаются. Тот факт, что рассматриваемый композит представляет собой матричную смесь, подтверждает инвариантность формул (19) относительно одновременной перестановки величин (с, , ц, , К,) и (с2 , ц2 , К2).

Если в соотношениях (14) принять ц = ц* ; К = К * , то относительно эффективных модулей упругости многокомпонентного композита получаются уравнения так называемой самосогласованной модели

C1-V1

C2 -V2

vS +aS-

= 1,

C1 - K1

C2 - K2

K* +YS-

= 1,

где

2 4 - 5-v*

15 1 -v*

Ys =

1 1 + v*

3 1 -v*

3-K* -2-vS

6- К* + 2-vS

На рис. 1 показаны верхние и нижние границы Фойгта и Рейсса для эффективного модуля упругости сдвига, рассчитанные по формулам (18), верхние и нижние границы Хашина - Штрикмана, рассчитанные по формулам (15) и (16), и кривая модуля упругости сдвига для матричной смеси, рассчитанная по формулам (19).

Mi

20

10

Mr// 1

/Îh V Ml/ Mc У J Mf

0

0,5

1,0

Рис. 1. Графики функций верхних и нижних границ Фойгта и Рейсса для эффективных модулей упругости сдвига, рассчитанных по формулам

Fig. 1. Graphs of the functions of the upper and lower Voigt and Reuss boundaries for the effective shear elastic moduli calculated by formulas

Расчетные значения: 1; 0,35 ; ц2 = 20; V = 0,25.

В самом общем случае вид неопределенных модулей упругости тела сравнения ц и К может быть представлен соотношениями

[m = V„ = U1 (C1 )-M1 + U2 ( C2 ) - V2 ,

l-K = Ku = U1 ( с )-K + U2 (C2 )-K2

где величины щ (с), щ (с2) - монотонно возрастающие

функции, описывающие связность составляющих компонентов композитов и удовлетворяющие условиям

0 < u (с )< 1, u (0) = 0, U (1)= 1, 0<u (с)< 1, U (0) = 0, U(1) = 1,

U1 (с1 )+ U2 (с2 )= 1

Подстановка величин (21) в соотношения (12) позволяет получить формулы для усредненных по объему V2 деформаций

W2 =(1 + au -(Ч - 1))-^u {j, (в„)2 =(1 + Yu -(kU - OH-(S•

c

2

a

- с1 - с2

as =

VS =

Формулы (14) для эффективных модулей упругости принимают вид

(

Ци = Ц»

К* = К»

c2 •(m» -m» )-(l + g„ •(m» -1)) 1 + gu •((m -1) + Ci •(mU -m)) C2 •(k» -kU )•( 1 + Yu •(kU - 1 ))'

1 '1 + Yu •

m» +

kU +

где

g = А•4-5•v,  Y = L1±^V =_ 3• К -2•ц»

1 5 1 -v, 3 1 -v„' " 2•

и" и" JC» К»

m» = m'2 =0-; k» = ku» =

2  1 К 2 К

ц»

ц»

На практике изготовление дисперсно-упрочненных композиционных материалов обычно сопровождается прессованием, брикетированием, спеканием дисперсных частиц и пропиткой расплавом связующего компонента. При малых объемных содержаниях второго компонента композиционного материала с2 распределенные в объеме первого компонента V дисперсные частицы практически не влияют друг на друга, их взаимодействие между собой и связность будут минимальными. В этом случае основной вклад в несущую способность композиционного материала обеспечивает первый компонент, образуя связующую матрицу, в которой хаотически распределены дисперсные включения второго компонента, а функция связности щ (с2) будет принимать малые значения, сравнимые со значениями величины с . С ростом объемного содержания второго компонента его брикетируемые и спекаемые дисперсные частицы образуют структуры, которые вместе с материалом первого компонента формируют взаимопроникающие друг в друга каркасы. При этом взаимодействие дисперсных включений в этих каркасах существенно увеличивается и сопровождается соответствующим ростом функции связности щ (с) и ее

производной.

При объемных содержаниях второго компонента композита с , близких к единице, компоненты композиционного материала могут меняться ролями. Основной вклад в несущую способность композита начинает вносить материал второго компонента, формируя связующую матрицу, в которой хаотически распределяются включения из материала первого компонента, а функция связности щ (с2) начинает принимать значения, близкие к единице.

Многочисленные экспериментальные наблюдения показывают, что при технологическом изготовлении образцов дисперсно-упрочненного композиционного материала для одних и тех же значений объемных содержаний дисперсных частиц с эффективные модули

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

Очевидно, что изменение функции связности щ (с) должно соответствовать изменениям объемов компонентов V и V . Так, например, уменьшению объема первого компонента будет соответствовать увеличение функции связности щ (с2), и наоборот. Таким образом, естественно предположить, что положительное приращение функции связности Дщ > 0 будет пропорционально относительному приращению первого

ДК п

компонента Дю =--1 > 0.

V,

Относительное приращение Дю связано с приращением Дс2 соотношением

Дю = -

AV_ ДК2_ Д(с2 •К) _ Дс2 _ Дс2

V V с •к с 1 - с

Таким образом, переменные Дю и Дс2 связаны уравнением

d ю = ^-

1 - с,

или

= 1 - exp (-ю).

Здесь использовано очевидное начальное условие

с2 0= 0.

2 1ю=0

Приращение функции связности Дщ (ю) = = щ (ю+Дю)-щ (ю) на некотором промежутке Дю можно представить в виде трех слагаемых:

.А , .

Ды2 = Дм2 + Дм2 + Дм2

Здесь Д»А - частичное приращение функции связности на отрезке Дю за счет слабого взаимодействия отдельных частиц; Дм^ - частичное приращение функции связности на отрезке Дю за счет существенного взаимодействия отдельных частиц, образующих каркасные структуры; AuW - случайные колебания приращения функции связности на отрезке Дю за счет воздействия случайных внешних факторов при технологическом изготовлении образцов дисперсно-упрочненного композиционного материала.

Величины Д»А , Д^, Ды2 можно представить в виде

Д»А (ю) = A (1 - м2 ((!>))• Дю,

Дн^ (ю) = B • u (ю)(1 - м2 (го^Дю,

Ды2 (ю) = ст • м2 (ю)(1 - ы\ (ю))•6(ю)• Aw.

v

с

2

Здесь А - коэффициент, описывающий взаимодействие включений при их малых объемных содержаниях; В - коэффициент, описывающий интенсивное взаимодействие между собой включений и формирование из них групп каркасных образований; w - стандартный

винеровский процесс; Дн> = е(ю)^Дю , с - среднеквадратичное отклонение функции связности щ (ю) от своего математического ожидания; функция е(со) ~ А^(ОД) - случайная величина с нормальным законом распределения, нулевым средним значением (е) = 0 и единичной дисперсией ^е2) = 1 [40].

Множитель (1 - щ (ю)) соотношений (27) соответствует замедлению роста функции связности щ (ю) ,

при котором достаточно большая совокупность включений и связующая матрица образуют взаимопроникающие каркасы. Параметр X описывает интенсивность

стремления функции (1 - щ£(ю)) к своему предельному

нулевому значению.

Подстановка формул (27) в соотношение (26) дает

Дщ (ю) = ((А + В • щ («))• Дю+сщ (ю)•е(ю)• Дw)х

х

Предельный переход в уравнении (28) при Дю ^ 0 и Д-^ ^ 0 приводит к стохастическому дифференциальному уравнению диффузии Ито [40]

ёщ = Р (щ) • ё ю+ Q (щ) • ём>.

Здесь

Р

- коэффициент сноса уравнения (29),

Q (щ ) = <

•

- коэффициент волатильности уравнения (29). Начальное условие для уравнения (29) с коэффициентами (30) и (31) имеет вид

А 0)=0

Исключение из уравнений (25) и (29) вспомогательной переменной ю приводит к стохастическому дифференциальному уравнению диффузии Ито относительно объемного содержания второго компонента с

ёщ =

(А + В • щ )•-

- с,

1 - с

•

'2 J

Стохастическое дифференциальное уравнение (27) с начальным условием (26) решается численно с помо-

щью алгоритма, построенного в соответствии с методом Эйлера -Маруямы [40]

(г+1) (г) , Щ ' = щ ; +

[

«•вО

Дс

1 - с

.

(1 -щ1

При реализации алгоритма (34) на малом шаге Дс(г ' > 0 начиная с нулевого значения генерируется

случайное число е(,) = е (с(г )) и вычисляется последующие значение +1. Результатом работы этого алгоритма образуются последовательности |с(1)| и | г) |, которые на координатной плоскости образуют систему точек |с(,), щ1'! и соответствующую ей случайную траекторию. Каждый раз при повторении алгоритма (34) образуется новая случайная траектория, поскольку случайная величина е = е( с2) генерирует новые случайные значения.

В численных расчетах по алгоритму (34) отрезок с е[0,1] разбивался на п = 100 одинаковых частей

с постоянным значением шага Дс(,) = Дс = 0,01. Число реализаций случайной функции связности щ (с2) было принято т = 50. Следует отметить, что вблизи начальной точки |с(0) = 0,щ(0) = 0| и вблизи конечной точки

|с(п+1) = 1, щ(п+1) = 1| функция связности становится почти детерминированной.

Для вычисления математического ожидания функции связности М [щ ] необходимо выполнить статистическое усреднение стохастического дифференциального уравнения

ёМ [щ ]= М [(А + В • щ )•(l - 4)}-

и получить уравнение, содержащее статистические моменты

ёМ [щ ] 1 --—- =-х

ёс^ 1 с

х

Рекуррентное вычисление с помощью уравнения (35) статистических моментов вида М ] приводит к

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

В качестве такого допущения здесь предполагается, что флуктуации величины щ = щ (с2) относительно ее

+

X

2

щ

математического ожидания М [щ ] пропорциональны случайной величине е = е

щ -М[щ] = £-е .

Здесь ^ = М[щ]-(1 -М[щ- коэффициент

пропорциональности.

Тогда выражение для величины и%, с учетом формулы (36), принимает вид

= M [щ ]h

1 +

M [щ ]

Ограничиваясь в формуле (37) малыми флуктуаци-

< 1, рассмотрим три слагаемых сходяще-

M [щ ]

гося биномиального ряда

uh = M [щ ]'

( £ h-(h-l) е , ^

1 + А--г—г • s н---->---——г- • s +

M [щ ]

M [щ ]2

Вычисляя по формуле (38) средние величины

M [щл]« M [щ f M [щл+1 ]» M [щ f

i.в

2 M [щ ]2

'i+ 4^+1) .В2

2 M [щ ]2

[nwi]

и подставляя их в уравнение

-М [щ ] 1 - М [щ ]л , г ,ч

-^Т1=(А+5 - М [щ2])-

(1 - М щ ]л)2

W2

• M [щ f

1 - с

х

Начальное условие для уравнения (39) имеет вид

М[щ

Сравнение результатов численного решения задач Коши (39), (40) с численными значениями статистического среднего, вычисленного по всем т = 50 реализациям алгоритма (34), показывает их почти полное совпадение.

На рис. 2 представлены численные реализации решений алгоритма (34) в виде семейства стохастических

траекторий функции связности щ (с2) и кривая ее математического ожидания М [щ ] , построенная в соответствии с численным решением задачи Коши (3 9), (40). и2 (о X м [и ]

1,0

0,5

0,5

1,0

Рис. 2. Семейство стохастических траекторий функции связности щ

Fig. 2. The family of stochastic trajectories of the connection function щ

Расчетные значения: A = 0,1; B = 0,5 ; n = 100; m = 50; ст = 0,5 ; ^ = 1.

Эффективные стохастические модули упругости композиционного материала рассчитываются по формулам (23), (34) и (39), (40).

На рис. 3 представлены численные реализации рас-

w *

четов стохастических модулей упругости ц^ в виде семейства стохастических траекторий, построенных с помощью алгоритма (28), и кривая математического

ожидания стохастического модуля упругости M [ц* ] ,

построенная в соответствии с численным решением задачи Коши (39), (40). Все стохастические траектории

модуля ц^ и его математическое

ожидание M [ц* ]

расположены внутри вилки верхней и нижней границ Хашина - Штрикмана.

Расчетные значения: А = 0,1; В = 0,5 ; п = 100;

т = 50 ; а = 0,5 ; Х = 1; 1; V = 0,35 ; ц2 = 20 ;

V, = 0,25.

Применим полученные результаты для расчета модулей упругости композиционного материала, образованного из кобальта (Со) , упрочненного дисперсными

частицами карбида вольфрама (ЖС).

В работе [37] представлены экспериментальные данные по измерению эффективного модуля упругости Юнга двухкомпонентного композиционного материала

h

h

0

c

2

0

(Со-ЖС). З начения[ПW2] модуля упругости и коэффициента Пуассона первого компонента - кобальта -составляют Е = 206,27 ГПа и у1 = 0,30 . Значения модуля упругости и коэффициента Пуассона второго компонента - карбида вольфрама составляют Е = 703,18 ГПа и у2 = 0,22.

ц* (^ ) , Ы [ц* (€2 )]

0 0,5 1,0

Рис. 3. Семейство стохастических траекторий функции

модуля упругости ц* (с2), построенных с помощью

алгоритма (34), кривая математического ожидания

стохастического модуля упругости M |ц* (с2)] , построенная

в соответствии с численным решением задачи Коши (39), (40), и верхняя и нижняя границы Хашина - Штрикмана

Fig. 3. The family of stochastic trajectories of the elastic modulus function ц* (c2) constructed using algorithm (34), the mathematical expectation curve of the stochastic elastic modulus M |ц* (c2)] , constructed in accordance with the numerical

solution of the Cauchy problem (39), (40), and the upper and lower Khashin - Shtrikman bounds

На рис. 4. представлено сравнение численных реализаций расчетов стохастических модулей упругости

E*

в виде семейства стохастических траекторий, по-

строенных с помощью формул (23) и алгоритма (34).

700

450

200

/о Jfo

0

0,5

1,0

Рис. 4. Семейство стохастических траекторий функции модуля упругости Е*

и алгоритма (34). Светлыми кружками обозначены экспериментальные значения модуля упругости Юнга композита

Fig. 4. A family of stochastic trajectories of the elastic modulus function constructed E* (с2) using formulas (23) and algorithm

(34). The open circles indicate the experimental values of the Young's modulus of elasticity of the composite Расчетные значения: A = 0,2; B = 1,0 ; n = 100; m = 50; ст = 1,0 ; X = 1.

В работе [38] представлены экспериментальные данные по измерению эффективного модуля упругости Юнга еще одного двухкомпонентного композиционного материала (Co -WC). В этом случае значения модуля упругости и коэффициента Пуассона первого компонента - кобальта - составляют E = 207,93 ГПа и V[ = 0,31. Значения модуля упругости и коэффициента Пуассона второго компонента - карбида вольфрама -составляют E2 = 710,08 ГПа и v2 = 0,19 .

На рис. 5. представлено сравнение численных реализаций расчетов стохастических модулей упругости E* такого композиционного материала в виде семейства стохастических траекторий, построенных с помощью формул (23) и алгоритма (34).

E*

700

450

200

Рис. 5. Семейство стохастических траекторий функции модуля упругости E* (c2), построенных с помощью формул (23) и

алгоритма (34). Точками обозначены экспериментальные значения модуля упругости Юнга композита

Fig. 5. A family of stochastic trajectories of the elastic modulus function constructed EU (c2) using formulas (23) and algorithm (34).

The dots indicate the experimental values of the Young's modulus of elasticity of the composite

Расчетные значения: A = 0,25 ; B = 0,25; n = 100; m = 50; CT = !,°0 ; X = 1.

Заключение

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

c

c

U \ 2

c

2

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

В соответствии с методом Эйлера - Маруямы построен алгоритм численного решения этого стохастического дифференциального уравнения. Каждой реализации этого алгоритма соответствует своя стохастическая траектория случайной функции связности составляющих компонентов композиционного материала.

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

Библиографический список

1. Эшелби Дж. Континуальная теория дислокаций. - М.: Иностранная литература, 1963. - 248 с.

2. Hashin Z. Theory of mechanical behavior of heterogeneous media // Appl. Mech. Rev. - 1964. - Vol. 17, no. 1. - P. 1-10.

3. Hill R. A self-consistent mechanics of composite materials // J. Mech. Phys. Solids. - 1965. - Vol. 13, no. 4. - P. 213-222.

4. Хорошун Л.П. О методе определения упругих модулей армированных тел // Механика полимеров. - 1968. - № 1. -С. 78-87.

5. Yeh R.H.T. Variational principles of the elastic moduli of composite materials // J. Appl. Phys. - 1970. - Vol. 41, no. 8. -P. 3353-3356. DOI: 10.1063/1.1659424

6. Шермергор Т.Д. Теория упругости микронеоднородных сред. - М.: Наука, 1977. - 400 с.

7. Композиционные материалы: 8 т. Т. 2. Механика композиционных материалов / под ред. Дж. Сендецки. - М.: Мир, 1978. - 568 с.

8. Кристенсен Р. Введение в механику композитов. - М.: Мир, 1982. - 336 с.

9. Ванин Г.А. Микромеханика композиционных материалов. - Киев: Наукова думка, 1985. - 304 с.

10. Сараев Л.А. Математическое моделирование упруго-пластических свойств многокомпонентных композиционных материалов. - Самара: АНО СНЦ, 2017. - 222 с.

11. Алешин В.И. Эффективные упругие свойства двухфазных композитов // Журнал технической физики. - 2007. -Т. 77, № 9. - С. 54-60.

12. Ustinov K.B., Goldstein R.V. On application of classical Eshelby approach to calculating effective elastic moduli of dispersed composites // Int. J. Fract. - 2007. - Vol. 147, no. 1-4. -P. 55-66.

13. Zheng Q.-S., Du D.-X. An explicit and universally applicable estimate for the effective properties of multiphase composites which accounts for inclusion distribution // J. Mech. Phys. Solids. - 2001. - Vol. 49, no. 11. - P. 2765-2788.

14. Большаков В.И., Андрианов И.В., Данишевский В.В. Асимптотические методы расчета композитных материалов с учетом внутренней структуры. - Днепропетровск: Пороги, 2008. - 247 с.

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

Установлены новые макроскопические определяющие соотношения для микронеоднородных материалов с изменяемой микроструктурой и вычислены их эффективные модули упругости.

Значения эффективных модулей упругости, построенных по полученным в статье выражениям, лежат внутри вилки Хашина - Штрикмана для нижней и верхней границы эффективных модулей упругости композиционных материалов.

Численный анализ разработанных моделей показал хорошее соответствие известным экспериментальным данным.

15. Димитриенко Ю.И., Соколов А.П. Об упругих свойствах композиционных материалов // Математическое моделирование. - 2009. - Т. 21, № 4. - С. 96-110.

16. Рогачев Е.А., Зверев М.А., Суриков В.И. Прогнозирование упругих свойств полимерных композиционных материалов // Омский научный вестник. Серия: Приборы, машины и технологии. - 2010. - № 3 (93). - С. 28-31.

17. Димитриенко Ю.И., Соколов А.П. Многомасштабное моделирование упругих композиционных материалов // Математическое моделирование. - 2012. - Т. 24, № 5. - С. 3-20.

18. Abedini A., Chen Z.T. A micromechanical model of particle-reinforced metal matrix composites considering particle size and damage // Computational Materials Science. - 2014. -Vol. 85. - P. 200-205. DOI: 10.1016/i.commatsci.2014.01.012

19. Dominguez-Rodriguez G., Tapia A., Aviles F. An assessment of finite element analysis to predict the elastic modulus and Poisson's ratio of singlewall carbon nanotubes // Computational Materials Science. - 2014. - Vol. 82. - P. 257-263. DOI: 10.1016/j.commatsci.2013.10.003

20. Головин Н.Н., Зарубин В.С., Кувыркин Г.Н. Оценка эффективных упругих характеристик материалов, модифицированных фуллереном // Композиты и наноструктуры. -2011. - № 4. - С. 21-31.

21. Моделирование зависимостей физико-механических характеристик от параметров микро- и наноструктуры полимерных композиционных материалов / С.А. Лурье, Ю.М. Миронов, В.А. Нелюб, А.С. Бородулин, И.В. Чуднов, И.А. Буянов, Ю.О. Соляев // Наука и образование: научное издание МГТУ им. Н.Э. Баумана. - 2012. - № 6. - С. 38-60. DOI: 10.7463/0612.0431339

22. Modelling of mechanical properties of nanoparticle-filled polyethylene / В. Lehmann, A.K. Schlarb, K. Friedrich, M.Q. Zhang, M.Z. Rong // Int. J. Polym. Mater. - 2008. - Vol. 57, no. 1. - P. 81-100.

23. Bulk Nanostructured Materials: Non-Mechanical Synthesis / Y. Ivanisenko, A. Darbandi, S. Dasgupta, R. Kruk, H. Hahn // Advanced Engineering Materials. - 2010. - Vol. 12, no. 8. - P. 666-676. DOI: 10.1002/adem.201000131

24. Лурье С.А., Соляев Ю.О. Модифицированный метод Эшелби в задаче определения эффективных свойств со сферическими микро-и нановключениями // Вестник

Пермского государственного технического университета. Механика. - 2010. - № 1. - С. 80-90.

25. Зайцев А.В., Фукалов А.А. Эффективные модули объемного сжатия дисперсноупрочненных композитов со сплошными и полыми анизотропными сферическими включениями // Вестник Пермского государственного технического университета. Механика. - 2010. - № 4. -С. 46-54.

26. Study of Reinforcing Mechanisms in TRIP-Matrix Composites under Compressive Loading by Means of Micromechanical Simulations / S. Priiger, L. Mehlhorn, U. Muhlich, M. Kuna // Advanced Engineering Materials. - 2013. - Vol. 15, no. 7. -P. 542-549. DOI: 10.1002/adem.201200323

27. Влияние физикомеханических характеристик полых стеклянных микросфер на свойства сферопластиков / И.И. Соколов, М.Г. Долматовский, И.С. Деев, В.Я. Стеценко // Пластические массы. - 2005. - № 7. - С. 16-18.

28. Influence of Powder Size on Production Parameters and Properties of Syntactic Invar Foams Produced by Means of Metal Powder Injection Moulding / J. Weisel, N. Salk, U. Jehring, J. Baumeister, D. Lehmhus, M.A. Bayoumi // Advanced Engineering Materials. - 2013. -Vol. 15, no. 3. - P. 118-122. DOI: 10.1002/adem.201200129

29. Jin H.-J., Weissmuller J. Bulk Nanoporous Metal for Actuation // Advanced Engineering Materials. - 2010. - Vol. 12, no. 8. - P. 714-723. DOI: 10.1002/adem.200900329

30. Schmidt K., Becker J. Generating Validated 3D Models of Microporous Ceramics // Advanced Engineering Materials. - 2013. -Vol. 15, no. 1-2. - P. 40-45. DOI: 10.1002/adem.201200097

31. Giraud A., Sevostianov I. Micromechanical modeling of the effective elastic properties of oolitic limestone // Int. J. Rock Mech. Min. - 2013. - Vol. 62. - P. 23-27.

32. Effect of Pore Morphology on Deformation Behaviors in Porous Al by FEM Simulations / Y.J. Cho, W.J. Wook Jin Lee,

References

1. Eshelbi Dzh. Kontinual'naia teoriia dislokatsii [Continual theory of dislocations]. Moscow, Inostrannaia literatura, 1963, 248 p.

2. Hashin Z. Theory of mechanical behavior of heterogeneous media. Appl. Mech. Rev., 1964, vol. 17, no. 1, pp. 1-10.

3. Hill R. A self-consistent mechanics of composite materials. J. Mech. Phys. Solids., 1965, vol. 13, no. 4, pp. 213-222.

4. Khoroshun L.P. O metode opredeleniia uprugikh modulei armirovannykh tel [On the method for determining the elastic moduli of reinforced bodies], Mekhanika polimerov, 1968, no. 1, pp. 78-87.

5. Yeh R.H.T. Variational principles of the elastic moduli of composite materials. J. Appl. Phys, 1970, vol. 41, no. 8, pp. 33533356. DOI: 10.1063/1.1659424

6. Shermergor T.D. Teoriia uprugosti mikroneodnorodnykh sred [Theory of elasticity of micro-inhomogeneous media]. Moscow, Nauka, 1977, 400 s.

7. Sendetski Dzh. Kompozitsionnye materialy 8 t. T 2. Mek-hanika kompozitsionnykh materialov [Mechanics of composite materials]. Moscow, Mir, 1978, 568 p.

8. Kristensen R. Vvedenie v mekhaniku kompozitov [Introduction to the mechanics of composites]. Moscow, Mir, 1982, 336 p.

9. Vanin G.A. Mikromekhanika kompozitsionnykh material-lov [Micromechanics of composite materials]. Kiev, Naukova dumka, 1985, 304 p.

10. Saraev L.A. Matematicheskoe modelirovanie uprugoplas-ticheskikh svoistv mnogokomponentnykh kompozitsionnykh materi-alov [Mathematical modeling of elastoplastic properties of multicom-

S.K. Park, Y.H. Park // Advanced Engineering Materials. - 2013. -Vol. 15, no. 3. - P. 166-169. DOI: 10.1002/adem.201200145

33. Deqing W. Relation of Cell Uniformity and Mechanical Property of a Close Cell Aluminum Foam // Advanced Engineering Materials. - 2013. - Vol. 15, no. 3. - P. 175- 179. DOI: 10.1002/adem.201200135

34. Montero-Chacon F., Marin-Montin J., Medina F. Meso-mechanical characterization of porosity in cementitious composites by means of a voxel-based finite element model // Computational Materials Science. - 2014. - Vol. 90. - P. 157-170. DOI: 10.1016/i.commatsci.2014.03.066

35. Functionalized Metallic Hollow Sphere Structures / H. Goehler, U. Jehring, J. Meinert, R. Hauser, P. Quadbeck, K. Kuemmel, G. Stephani, B. Kieback // Advanced Engineering Materials. - 2014. - Vol. 16, no. 3. - P. 335-339. DOI: 10.1002/adem.201300057

36. Rapid Sintering of Porous Monoliths Assembled from Microbeads with High Specific Surface Area and Multimodal Porosity / Th.C. Schumacher, T.Y. Klein, L. Treccani, K. Rezwan // Advanced Engineering Materials. - 2014. - Vol. 16, no. 2. -P. 151-155. DOI: 10.1002/adem.201300220

37. Nishimatsu C., Gurland J. Experimental servey of the deformation of the hardductile two-phase alloy system W-Co // Trans. Amer. Soc. Metals. - 1960. - Vol. 52, No. 2. -P. 469-484.

38. A systematic investigation of elastic moduli of W-Co alloys / H. Doi, Y. Fujiwara, K. Miyake, Y. Oosawa // Met. Trans. -1970. - Vol. 1, no. 5. - P. 1417-1425.

39. Современные композиционные материалы / под ред. Л.И. Браутмана и Р. Крока. - М.: Мир, 1970. - 672 с.

40. Кузнецов Д.Ф. Стохастические дифференциальные уравнения: теория и практика численного решения. - СПб.: Изд-во Политехн. ун-та, - 2007. - 800 с.

ponent composite materials]. Samara, ANO SNTs, 2017, 222 p.

11. Aleshin V.I. Effektivnye uprugie svoistva dvukhfaznykh kompozitov [Effective elastic properties of two-phase composites]. Zhurnal tekhnicheskoi fiziki, 2007, vol. 77, no 9, pp. 54-60.

12. Ustinov K.B., Goldstein R.V. On application of classical Eshelby approach to calculating effective elastic moduli of dispersed composites. Int. J. Fract, 2007, vol. 147, no. 1-4, pp. 55-66.

13. Zheng Q.-S., Du D.-X. An explicit and universally applicable estimate for the effective properties of multiphase composites which accounts for inclusion distribution. J. Mech. Phys. Solids, 2001, vol. 49, no. 11, pp. 2765-2788.

14. Bol'shakov V.I., Andrianov I.V., Danishevskii V.V. Asimptoticheskie metody rascheta kompozitnykh materialov s uchetom vnutrennei struktury [Asymptotic methods for calculating composite materials taking into account the internal structure]. Dnepropetrovsk, Porogi. 2008, 247 p.

15. Dimitrienko Iu.I., Sokolov A.P. Ob uprugikh svoistvakh kompozitsionnykh materialov [On elastic properties of composite materials // Mathematical modeling]. Matematicheskoe modeliro-vanie, 2009, vol. 21, no. 4, pp. 96-110.

16. Rogachev E.A., Zverev M.A., Surikov V.I. Prognoziro-vanie uprugikh svoistv polimernykh kompozitsionnykh materialov [Prediction of elastic properties of polymer composite materials]. Omskii nauchnyi vestnik. Ser. Pribory, mashiny i tekhnologii, 2010, no. 3 (93), pp. 28-31.

17. Dimitrienko Iu.I., Sokolov A.P. Mnogomasshtabnoe modelirovanie uprugikh kompozitsionnykh materialov [Multiscale

modeling of elastic composite materials], Matematicheskoe mod-elirovanie, 2012, vol, 24, no, 5, pp, 3-20,

18, Abedini A,, Chen Z,T, A micromechanical model of particle-reinforced metal matrix composites considering particle size and damage, Computational Materials Science, 2014, vol, 85, pp, 200-205, DOI: 10,1016/i,commatsci,2014,01,012

19, Dominguez-Rodriguez G,, Tapia A,, Aviles F, An assessment of finite element analysis to predict the elastic modulus and Poisson's ratio of singlewall carbon nanotubes. Computational Materials Science, 2014, vol, 82, pp, 257-263, DOI: 10,1016/j,commatsci,2013,10,003

20, Golovin N,N,, Zarubin V,S,, Kuvyrkin G,N, Otsenka effektivnykh uprugikh kharakteristik materialov, modifitsiro-vannykh fullerenom [Estimation of effective elastic characteristics of materials modified with fullerene], Kompozity i nanostruktury, 2011, no, 4, pp, 21-31,

21, Lur'e S,A,, Mironov Iu,M,, Neliub V,A,, Borodulin A,S,, Chudnov I,V,, Buianov I,A,, Soliaev Iu,O, Modelirovanie zavisi-mostei fiziko-mekhanicheskikh kharakteristik ot parametrov mikro- i nanostruktury polimernykh kompozitsionnykh materialov [Modeling the dependences of physical and mechanical characteristics on the parameters of the micro- and nanostructure of polymer composite materials], Nauka i obrazovanie: nauchnoe izdanie MGTU im, N,E, Baumana, 2012, no, 6, pp, 38-60, DOI: 10,7463/0612,0431339

22, Lehmann V,, Schlarb A,K,, Friedrich K,, Zhang M,Q,, Rong M,Z, Modelling of mechanical properties of nanoparticle-filled polyethylene, Int, J, Polym, Mater, 2008, vol, 57, no, 1, pp, 81-100,

23, Ivanisenko Y, Darbandi A,, Dasgupta S,, Kruk R,, Hahn H, Bulk Nanostructured Materials: Non-Mechanical Synthesis, Advanced Engineering Materials, 2010, vol, 12, no, 8, pp, 666676, DOI: 10,1002/adem,201000131

24, Lur'e S,A,, Soliaev Iu,O, Modified Eshelby's method in the problem of determining effective properties with spherical micro and nanoinclusions, PNRPU Mechanics Bulletin, 2010, no, 1, pp, 80-90,

25, Zaitsev A,V,, Fukalov A,A, Effective moduli of bulk compression of dispersion-strengthened composites with solid and hollow anisotropic spherical inclusions, PNRPU Mechanics Bulletin, 2010, № 4, S. 46-54,

26, Priiger S,, Mehlhorn L,, Muhlich U,, Kuna M, Study of Reinforcing Mechanisms in TRIP-Matrix Composites under Com-pressive Loading by Means of Micromechanical Simulations, Advanced Engineering Materials, 2013, vol, 15, no, 7, pp, 542-549, DOI: 10,1002/adem,201200323

27, Sokolov I,I,, Dolmatovskii M,G,, Deev I,S,, Stetsen-ko V,Ia, Vliianie fizikomekhanicheskikh kharakteristik polykh stekliannykh mikrosfer na svoistva sferoplastikov [Influence of physicomechanical characteristics of hollow glass microspheres on the properties of spheroplastics], Plasticheskie massy, 2005, no, 7,

pp. 16-18.

28. Weisel J., Salk N., Jehring U., Baumeister J., Lehmhus D., Bayoumi M.A. Influence of Powder Size on Production Parameters and Properties of Syntactic Invar Foams Produced by Means of Metal Powder Injection Moulding. Advanced Engineering Materials, 2013, vol. 15, no. 3, pp. 118-122. DOI: 10.1002/adem.201200129

29. Jin H.-J., Weissmuller J. Bulk Nanoporous Metal for Actuation. Advanced Engineering Materials, 2010, vol. 12, no. 8, pp. 714-723. DOI: 10.1002/adem.200900329

30. Schmidt K., Becker J. Generating Validated 3D Models of Microporous Ceramics. Advanced Engineering Materials, 2013, vol. 15, no. 1-2, pp. 40-45. DOI: 10.1002/adem.201200097

31. Giraud A., Sevostianov I. Micromechanical modeling of the effective elastic properties of oolitic limestone. Int. J. Rock Mech. Min, 2013, vol. 62, pp. 23-27.

32. Cho Y. J., Wook Jin Lee W.J., Park S.K., Park Y.H. Effect of Pore Morphology on Deformation Behaviors in Porous Al by FEM Simulations. Advanced Engineering Materials, 2013, vol. 15, no. 3, pp. 166-169. DOI: 10.1002/adem.201200145

33. Deqing W. Relation of Cell Uniformity and Mechanical Property of a Close Cell Aluminum Foam. Advanced Engineering Materials, 2013, vol. 15, no. 3, pp. 175-179. DOI: 10.1002/adem.201200135

34. Montero-Chacon F., Marin-Montin J., Medina F. Me-somechanical characterization of porosity in cementitious composites by means of a voxel-based finite element model. Computational Materials Science, 2014, vol. 90, pp. 157-170. DOI: 10.1016/i.commatsci.2014.03.066

35. GoehlerH., JehringU., Meinert J., HauserR., Quad-beckP., Kuemmel K., Stephani G., Kieback B. Functionalized Metallic Hollow Sphere Structures. Advanced Engineering Materials, 2014, vol. 16, no. 3, pp. 335-339. DOI: 10.1002/adem.201300057

36. Schumacher Th.C., Klein T.Y., Treccani L., Rezwan K. Rapid Sintering of Porous Monoliths Assembled from Microbeads with High Specific Surface Area and Multimodal Porosity. Advanced Engineering Materials, 2014, vol. 16, No. 2, pp. 151-155. DOI: 10.1002/adem.201300220

37. Nishimatsu C., Gurland J. Experimental servey of the deformation of the hardductile two-phase alloy system W-Co. Trans. Amer. Soc. Metals, 1960, vol. 52, no. 2, pp. 469-484.

38. Doi H., Fujiwara Y., Miyake K., Oosawa Y. A systematic investigation of elastic moduli of W-Co alloys. Met. Trans, 1970, vol. 1. no. 5, pp. 1417-1425.

39. Brautmana L.I., Kroka I.R. Sovremennye kompozi-tsionnye materialy. Mir, 1970, 672 p.

40. Kuznetsov D.F. Stokhasticheskie differentsial'nye urav-neniia: teoriia i praktika chislennogo resheniia [Stochastic differential equations: theory and practice of numerical solution], Saint Petersburg, Izd-vo Politekhnicheskogo un-ta, 2007, 800 p.