2019. 15(6). 407-414 Строительная механика инженерных конструкций и сооружений Structural Mechanics of Engineering Constructions and Buildings
HTTP://JOURNALS.RUDN.RU/STRUCTURAL-MECHANICS
Расчет и проектирование строительных конструкций Analysis and design of building structures
DOI 10.22363/1815-5235-2019-15-6-407-414 научная статья
УДК 691.3
Эффективные модули двухфазных строительных композитов с зернистым заполнителем
В.Т. Ерофеев, А.С. Тюряхин, Т.П. Тюряхина, А.В. Тиньгаев*
Национальный исследовательский Мордовский государственный университет имени Н.П. Огарева, Российская Федерация,
430005, Саранск, ул. Большевистская, 68
История статьи:
Поступила в редакцию: 19 октября 2019 г. Доработана: 02 декабря 2019 г. Принята к публикации: 06 декабря 2019 г.
Для цитирования
Ерофеев В.Т., Тюряхин А.С., Тюряхина Т.П., Тиньгаев А.В. Эффективные модули двухфазных строительных композитов с зернистым заполнителем // Строительная механика инженерных конструкций и сооружений. 2019. Т. 15. № 6. С. 407-414. http://dxdoi.org/10.22363/1815-5235-2019-15-6-407-414
Аннотация
В книге Р.М. Кристенсена «Введение в механику композитов» (1982) приведена расчетная формула для объемного модуля полидисперсных композитов со сферическими включениями. Эта формула известна русскоязыиному читателю почти 40 лет, но, к сожалению, не используется в практике строительного материаловедения. Для выявления прикладных возможностей формула Р.М. Кристенсена видоизменяется и сводится к безразмерной функции k = k (w, n, 6), зависящей от трех безразмерных параметров, то есть зависимой от трех величин: w - объемной доли включения, n - отношения модуля сдвига материала матрицы к величине объемного модуля той же матрицы, 6 - отношения объемных модулей материалов матрицы и включения. Численные исследования этой функции выывляют, что в двухфазных зернистык композитах существенно сужается область значений эффективных модулей по сравнению с областью, ограничиваемой оценками Фойгта и Рейсса (в смысле верхней и нижней границ реальных значений). При этом нижняя оценка по Кристенсену совпадает с оценкой по Рейссу. Приведены численные и графически оформленные результаты на примерах исследования двух характерных групп композиционных материалов. Кроме того, безразмерная форма эффективного модуля позволяет построить в плоском пространстве k - w систему наглядный графических зависимостей функций k(w). При разных значениях 6 функцией k = k(w, n) отображается пучок криволинейных отрезков, которым задается положение плоской фигуры в плоском пространстве. Приведены примеры построения фигур для характерных областей значений функции k(n, 6, w).
Ключевые слова: двухфазная модель зернистого композита; сферическая форма фаз матрицы и заполнителя; эффективный объемный модуль упругости композита; вилка Фойгта - Рейсса для эффективного модуля
Введение
Интенсивное развитие строительных композиционных материалов (СКМ) стимулирует их тео-
Ерофеев Владимир Трофимович, д. т. н., профессор, заведующий кафедрой строительных материалов и технологий; академик Российской академии архитектуры и строительных наук; Scopus Author ID: 57080123300, eLIBRARY SPIN-код: 4425-5045.
Тюряхин Алексей Сергеевич, к. т. н., доцент кафедры прикладной механики. Тюряхина Татьяна Павловна, аспирант кафедры строительных материалов и технологий.
Тиньгаев Александр Васильевич, магистрант кафедры строительных материалов и технологий.
© Ерофеев В.Т., Тюряхин А.С., Тюряхина Т.П., Тиньгаев А.В., 2019
г————. This work is licensed under a Creative Commons Attribution 4.0 (ccJ^Ci^l International License
https://creativec0mm0ns.0rg/licenses/by/4.Q/
ретическое исследование. При этом развивается и углубляется теория СКМ, направленная на рост числа прикладных задач механики композитов [2-4; 6-17,]. В данном исследовании решается одна из таких задач.
В механике СКМ модель представительного объема (модель ПО) для зернистого композита достаточно часто имеет вид двухфазной модели, в которой в матрицу (толстостенную сферу) включен шарообразный заполнитель (рис. 1). Материал каждой из фаз представлен упругим твердым телом (сплошным однородным и изотропным).
Рис. 1. Диаметральное сечение модели ПО зернистого композита [Figure 1. The diametrical section of the model software granular composite]
Под действием равномерного наружного давления q модель ПО деформируется, в результате чего в точках раздела фаз возникает контактное давление p. Равновесное состояние такой неоднородной гетерогенной модели ПО (рис. 1) можно описать совокупностью однородных моделей, составленной из трех расчетных схем (рис. 2).
Рис. 2. Расчетные схемы модели ПО, представленные: а - эффективной однородной моделью (сплошной шар); б - моделью однородной матрицы (пустотелый шар); в - моделью однородного заполнителя (сплошной шар) [Figure 2. The design schemes of the software model, presented: а - effective homogeneous model (solid ball); б - a model of a homogeneous matrix (hollow ball); в - a model of a homogeneous filler (solid ball)]
Контактное давление p можно найти из условия равенства радиальных перемещений в смежных точках матрицы и заполнителя (рис. 2, б и в). Такая задача решается в элементарных алгебраических выражениях, если воспользоваться применительно к матрице (рис. 2, б) решением задачи Ламе для толстостенной сферы [1; 18]. Не приводя выкладок, в результате получим следующее выражение:
p = q
1 + 2Пм
(1+2Пм w)+2^ev
(1)
где м - объемная доля материала заполнителя, м = а3/Я3 (рис. 1); V - объемная доля материала матрицы, V = 1-м.
Кроме того, в формуле (1) приняты обозначения, использованные авторами работ [5] и [18]:
Пм =-
1 - 2к. 2G
1 + к.
3K„
e K
e=—= м
e. Kз
(2)
e.. =-
Ем 1 - 2к
= 3K g =
-^м' Ом
Ем 1+ К,
= 2G
(3)
где Ем и Vм - модуль Юнга и коэффициент Пуассона материала матрицы соответственно; Км - объемный модуль; Ом - модуль сдвига матрицы; ем -шаровой модуль объемных деформаций; gм - девиа-торный модуль сдвига матрицы; пм - отношение модулей матрицы; К и е3 - объемный и шаровой модули материала заполнителя соответственно; 0 -отношение объемных модулей материалов матрицы и заполнителя.
Если теперь с учетом (1) приравняем радиальные перемещения в наружных точках эффективной модели (рис. 2, а) и модели матрицы (рис. 2, б), то в результате получим (снова без выкладок) формулу для вычисления величины безразмерного эффективного модуля двухфазного ком-
позита:
к =
К =(1+2Пм w )+2^ev
k м V+е(2Пм + w)
(4)
где к - объемный модуль композита (К), выраженный в долях от модуля матрицы (Км).
При этом обратная формула для вычисления эффективного объемного модуля К примет вид
К = К м к = К
(1+2Пм w )+2^ev V+е(2Пм + w) '
(5)
Ранее в работе Р.М. Кристенсена [17] была получена формула, аналогичная формуле (5), которая в наших обозначениях имеет вид
К - Км
_м
Кз - К
w
(1 - w)(Кз - Км)
(6)
1+
Км + т G
В такой форме Кристенсеном с учетом парадигмы Хашина [19] получена формула для вычисления величины объемного модуля упругости (К) в полидисперсных композитах со сферическим включением (заполнителем). Простые преобразования формулы (6) с учетом обозначений (1) - (3) приводят нас снова к формулам (4) и (5). Следовательно, формулы (4), (5) и (6) идентичны друг другу. Способ получения формулы (4) посредством использования решения задачи Ламе менее трудоемок, так как требует выполнения только элементарных алгебраических преобразований. Кри-стенсен же получил формулу (6) более сложным образом, непосредственно используя уравнения линейной теории упругости.
1. Расчетные формулы эффективного модуля (к) композитов
В формулах (4) - (6) зафиксирован факт того, что материал моделей на расчетных схемах, имеющих форму шара (рис. 2, а и в), представлен только одной упругой константой (модулем К или Кз), а материал матрицы, имеющей форму толстостенной сферы (рис. 2, б) - двумя константами (модулями Км и Ом). Этот факт с точки зрения упругих свойств служит признаком особенностей определения эффективного модуля К, присущих только двухфазным композитам, в отличие от эффективного модуля других (многофазных) композитов.
Исследуя формулу (4), заметим, что параметры Пм и ^ имеют одинаковый и удобный интервал (0, 1) их значений. Однако значения модуля к и параметра 0 изменяются в другом интервале (0, да) чисел. Поэтому интервал (0, да) расчленим на 2 интервала: интервал малых (0, 1] и интервал больших [1, да) чисел. В таком случае получим два рабочих варианта формулы (4):
- вариант 1 - для чисел модуля (к > 1) при значениях 0 = 01 < 1:
к _ К_(1 + 2Лмw) + 2Лм91К ; 1 Км V + 9i(2nM + w) '
(7)
- вариант 2 - для чисел модуля (fe < 1) при значениях 9 = 92 > 1:
k _
Кг_ _(1 + 2Пмw) + 2^V Км V + 92(2Пм + w)
(8)
Особо отметим, что в формулах (7) и (8) параметры 01 и 02 имеют инверсивную связь (с коэффициентом инверсии, равным единице), так как их произведение 01 • 02 = 1. Это свойство облегчает анализ расчетных формул. Но предварительно вспомним гипотезы Фойгта [21] и Рейса [20].
2. Плоское пространство значений модулей к
Применяя правило смеси, сначала определим, согласно [21], эффективный объемный модуль (К) композита (рис. 1, 2). В результате получим зависимость по Фойгту:
Кф _К V + Кз w .
(9)
Затем, применив то же правило, определим, согласно [20], эффективную объемную податливость (1/К) материала и получим зависимость по Рейссу:
1 V w
-=-+-. (10)
Кр К К
В отличие от формул (4) и (6) в зависимостях (9) и (10) влияние модуля сдвига матрицы (Ом) не учитывается. Именно в этом и заключается главным образом приближенность оценок по Фойгту и по Рейссу.
Величины модулей КФ и КР соответственно служат, как правило, верхней и нижней границей значений реальных модулей К композита. Однако, анализ решения Кристенсена (6), представленного формулами (4), (7) и (8), позволяет уточнить данное правило (в смысле верхней и нижней оценки вилки Фойгта - Рейсса) применительно к двухфазной структуре зернистых СКМ. Но сначала выразим модули Кф и Кр в долях от модуля Км (матрицы) и получим формулы для безразмерных величин тех же модулей:
Кф Кр 1
кф= — ^ + -, кр =—- =-. (11)
Ф Км 0 Р Км V + 0^
мм
Применительно к расчетным формулам (7) и (8) выражения (11) примут соответствующий вид: - вариант 1 (при 0 < 1):
w 1
кФ1 _(1 - w)+—, кР1 _-;
Ф1 V 7 9/ Р1 (1 - w) + 91w
(12)
- вариант 2 (при 9 > 1):
кф2 _ (1 - w) +
w
92
k Р2
1
(1 - w) + 92 w
. (13)
В случаях, когда величина отношения 0 неизменна (6 = const), выражения (12) и (13) становятся функциями одного аргумента: kф = k$>(w) и kp = kp(w). Следовательно, эти функции могут быть представлены графически в плоском пространстве осей k и w (рис. 3). О наглядности такого пространства, использованного при анализе формул (4) - (6), свидетельствует содержание работы [6].
2
1
! \
91 V
0,5
Рис. 3. Область плоского пространства k - w, ограниченная осями k, 6 и w. Пунктиром вывделены 2 отрезка функции kф = kф(0, w) при значениях 6 = 61 и 6 = 62 [Figure 3. The region of flat space k - w bounded by the axes k, 6 and w. The dashed line marks 2 segments of the function kф = kф (6, w) for the values 6 = 61 and 6 = 62]
k
9
9
w
0
1
В обозначенной области пространства к - м функции модулей по Фойгту, например, аналогично формулам (12) и (13), примут вид
кФ1 = (1 - +— , кФ2 = (1 - +— .
01 02
Эти функции отображаются двумя семействами прямолинейных отрезков, исходящих из их общей точки с координатами (0, 1). На рис. 3 пунктирными линиями изображены два таких отрезка, по одному из каждого семейства. Правые концы выбранных отрезков имеют координаты (1, 02) и (1, 01) (на рис. 3 выделены жирными точками).
Функции модулей по Рейсу и по Кристенсе-ну (в той же области пространства к - м) имеют криволинейное очертание; их графики расположены ниже отрезков по Фойгту и непосредственно примыкают к этим отрезкам своими концами в точках (0, 1), (1, 02) или (1, 01) (рис. 3). Для более четкого описания этих функций рассмотрим два численных примера.
3. Сравнительный анализ оценок по Фойгту, Рейссу и Кристенсену
Проведем анализ расчетных формул (7), (8) и формул (12) и (13) для сравнения верхней и ниж-
ней оценок величины эффективного модуля упругости в двухфазных композитах. С этой целью выделим две группы композитов. В первой из них отношение объемных модулей компонентов 0 примем равным 0,25, то есть меньше единицы; во второй - больше единицы (0 = 4).
Задавшись шагом Дм = 0,125, вычислим модули кф и кр по формулам (12), (13) как для первой, так и для второй группы композитов. Результаты вычислений представлены в табл. 1.
Затем для этих же групп вычислим модули к по формулам (7) и (8) при значениях отношения Пм, равных 1, 'Л и 0, которым соответствуют значения коэффициента Пуассона Vм, равные 0, 0,2 и 0,5. Результаты отображены в табл. 2.
Сопоставив таблицы, отметим, что значения строки 3 в табл. 2 совпали со значениями строки 2 в табл. 1. Значения же строки 6 в табл. 2 в свою очередь повторяют значения строки 4 в табл. 1. Таким образом, нижняя граница значений эффективных модулей к, вычисляемых по Кристенсену (табл. 2) точно совпала с оценкой модулей кР по Рейсу (табл. 1). Однако оценки по Фойгту и по Кристенсену различны. Это следует из сопоставления других строк тех же таблиц.
Таблица 1
Значения модулей по Фойгту (кф) и по Рейссу (kp) для двух групп композитов [Table 1. The values of the modules according to Voigt (кф) and Reuss (kp) for two groups of composites]
Группа [Group]
Модуль [Module]
Объемная доля заполнителя (w) [Volume fraction of aggregate (w)]
0 0,125 0,25 0,375 0,5 0,625 0,75 0,875 1
1 1,375 1,75 2,125 2,5 2,875 3,25 3,625 4
1 1,103 1,231 1,391 1,6 1,882 2,286 2,909 4
1 0,906 0,813 0,719 0,625 0,531 0,438 0,344 0,25
1 0,727 0,571 0,471 0,4 0,348 0,308 0,276 0,25
по Фойгту [according to Voigt] по Рейссу [according to Reuss] по Фойгту [according to Voigt] по Рейссу [according to Reuss]
Таблица 2
Значения модулей к (по Кристенсену) в двух группах композитов (0 = 1/4 и 0 = 4) при трех значениях параметра Пм [Table 2. The values of the modules к (according to Christensen) in two groups of composites (0 = 1/4 and 0 = 4)
for three values of the parameter Пм]
Параметры Объемная доля заполнителя (w) [Options]_[Volume fraction of aggregate (w)]
e Пм 0 0,125 0,25 0,375 0,5 0,625 0,75 0,875 w^1
0,25 1 1 1,2 1,429 1,692 2 2,364 2,8 3,333 4
0,25 0,5 1 1,162 1,353 1,581 1,857 2,2 2,636 3,211 4
0,25 0 1 1,103 1,231 1,391 1,6 1,882 2,286 2,909 4
4 1 1 0,88 0,769 0,667 0,571 0,483 0,4 0,323 0,25
4 0,5 1 0,86 0,739 0,633 0,538 0,455 0,379 0,311 0,25
4 0 1 0,727 0,571 0,471 0,4 0,348 0,308 0,276 0,25
Для большей наглядности рассмотрим два графических примера.
Пример 1. Совместим в одном рисунке графические зависимости, построенные по данным табл. 1 и 2 для первой группы выбранных материалов
(рис. 4). При этом все значения эффективных модулей лежат в интервале чисел 1 < к < 4.
Концевые точки всех криволинейных графиков (рис. 4) совмещены с концевыми точками прямолинейного отрезка по Фойгту. Нижний график зна-
чений модулей по Рейсу дает нижнюю оценку значений эффективных модулей реальных композитов. Область плоского пространства, ограниченная верхним и нижним графиками, наглядно демонстрирует так называемую вилку Фойгта - Рейсса.
-кФ ~к, пм=1
к, пм=0,5 -кР
группы материалов (рис. 5), но со значениями модулей к < 1.
0 0,125 0,25 0,375 0,5 0,625 0,75 0,875 1 Объемная доля заполнителя, w
Рис. 4. Зависимости модулей fcsi, kp1 и ki от величины объемной доли w. Значения ki соответствуют значениям г|м, равным 1, / и 0
[Figure 4. Dependences of the modules кф1, kp1 and k1 on the volume fraction w. The values of k1 correspond to the values of Пм equal to 1, /, and 0]
На рис. 4 также отображены зависимости значений модулей k1 по Кристенсену - три криволинейных графика, соответствующих значениям коэффициента Пуассона Vm = 0, 0,2 и 0,5. При этом нижний график (со значением Vm = 0,5) точно совпал с графиком модулей kp1 (по Рейссу). Следовательно, нижняя оценка по Рейссу, определяемая второй из формул (12), является и частным случаем формулы (7) при значении Пм = 0 (или Vm = 0,5, что то же самое).
Оценка по Фойгту (верхний график) не совпадает с верхней оценкой по Кристенсену (2-й график сверху) (рис. 4), которая, согласно выражению (7), при значении Пм = 1 определяется формулой
_ (1 + 2w) + 29tV x 1 _ V + 0j(2 + w)
(14)
Значениям тах к1 отвечают материалы матрицы,
коэффициент Пуассона которых Vм = 0. Такие материалы лучше других материалов (со значениями Vм > 0) сопротивляются деформациям сдвига. Другими словами, чем меньше величина коэффициента Vм, тем выше жесткость (и прочность) композита.
Пример 2. Аналогично поступим и со второй группой выбранных композитов. Совместив в одном рисунке зависимости, построенные по данным табл. 1 и 2, получим графики, подобные графикам первой
0,95
0,85 * 0,75
К 0,65 ё
I 0,55
■со
0,45 0,35 0,25
-0-кФ к, пм=1 " " к, Пм=0,5
0 0,125 0,25 0,375 0,5 0,625 0,75 0,875 1 Объемная доля заполнителя, w
Рис. 5. Зависимости модулей kф2, kp2 и k2 от величины объемной доли w. Значения k2 соответствуют значениям параметра Пм, равным 1, / и 0
[Figure 5. Dependences of the modules kф2, kp2 and k2 on the volume fraction w. The values of k2 correspond to the values of the parameter г|м equal to 1, / and 0]
Последовательность в расположении графиков на рис. 5 сохраняется такой же, как и на рис. 4. Основные различия рисунков заключаются в том, что: 1) прямолинейный отрезок на рис. 4 изображен восходящим, а на рис. 5 аналогичный отрезок - нисходящим; 2) рис. 4 характеризует композиты первой группы (при 0 < 1), рис. 5 - композиты второй группы (при 0 > 1).
Зависимости по Кристенсену (2-й и 4-й графики на рис. 5), определяемые формулой (8), при значениях пм = 1 и пм = 0 принимают вид функций
max к2
min к2
(1 + 2w)+292У V + 92(2 + w)
1
V + 92w
(15)
(16)
Таким образом, зависимость (16) как частный случай формулы Кристенсена (8) совпадает с формулой Рейса (13).
Заключение
Расчетная формула (6) Р.М. Кристенсена [17] для величины эффективного объемного модуля полидисперсных композитов существует уже почти 40 лет и имеет достаточно широкую известность. Тем не менее, к сожалению, она фактически не используется в строительном материаловедении. Вышеизложенные численные исследования этой формулы позволили выявить и графически проиллюстрировать ряд несомненных ее достоинств.
Во-первых, относительно верхних и нижних оценок величины объемных модулей реальных композитов формула выявляет суженную вилку оценок по Кристенсену, по сравнению с общеизвестной вилкой фойгта - Рейсса».
Во-вторых, нижние оценки по Рейссу (13) и по Кристенсену (16) совпадают. Это уточняет положение вилки по Кристенсену в рамках вилки фойгта -Рейсса.
В-третьих, наши интерпретации формулы Кристенсена имеют вид элементарных алгебраических функций трех безразмерных параметров (Пм = 2Gh/3Kh, 6 = KM/Ks и w = a3/R3). Это существенно облегчает их численный анализ и практическое пользование ими.
И, наконец, в-четвертых, безразмерная форма представления величины эффективного модуля зернистого композита (k = K/K) позволяет построить в плоском пространстве k - w систему бесконечного множества графических зависимостей модуля k от количественного содержания объемной доли заполнителя. Это дает возможность в каждом конкретном композите вычислить место его эффективного модуля в плоском пространстве модулей, определяемых формулой Кристенсена, и в случае необходимости усилить (или ослабить) те или иные механических свойства композита.
Список литературы
1. Безухов Н.И. Основы теории упругости, пластичности и ползучести. М.: Высшая школа 1968. 512 с.
2. Бобрышев А.Н., Ерофеев В.Т., Козомазов В.М. физика и синергетика дисперсно-неупорядоченных конденсированных композитных систем. СПб.: Наука, 2012. 176 с.
3. Васильев В.В., Протасов В.В., Болотин В.В. и др. Композитные материалы: справочник. М.: Машиностроение, 1990. 512 с.
4. Гусев Б.В., Кондрашенко В.И., Маслов Б.П., Фай-сович А. С. формирование структуры композиционных материалов и их свойства. М.: Научный мир, 2006. 566 с.
5. Ерофеев В.Т., Тюряхин А.С., Ерофеева И.В. О связях параметров несущей фазы с эффективными параметрами в моделях зернистых композитов // Строительная механика и расчет сооружений. 2018. № 3. С. 7-17.
6. Erofeev V., Tyuryakhin A., Tyuryakhina T. Flat space of values of volume module of grain composite with spherical fill-lem // International Journal of Civial Engineering and Technology (IJCIET). 2019. Issue 8. Рр. 333-342.
7. Берлин А.А., Вольфсон С.А., Ошмян В.Г., Ени-колопов Н.С. Принципы создания композитных полимерных материалов. М.: Химия, 1990. 240 с.
8. Аскадский А.А., Голенева Л.М., Бычко К.А., Казанцева В.В., Константинов К.В., Алмаева Е.С., Клин-ских А.Ф., Коврига О.В. Градиентные полимерные материалы // Российский химический журнал. 2001. Т. 45. № 3. С. 123-128.
9. Duan K., Hu Xiao, Wittmann F.H. Boundary effect on concrete fracture and non-constant fracture energy distribution // Engineering Fracture Mechanics. 2003. No. 70. Рр. 2257-2268.
10. Сarpinteri A., Chiaia B., Cornetti P. On the mechanics of quasi-brittle materials with a fractal microstructure // Engineering Fracture Mechanics. 2003. No. 70. Рр. 2321-2349.
11. Ayatollahi M.R., Akbardoost J. Size effects on fracture toughness of quasi-brittle materials-A new approach // Engineering Fracture Mechanics. 2012. No. 92. Рр. 89-100.
12. Hu X., Guan J., Wang Y., Keating A., Yang S. Comparison of boundary and size effect models based on new developments // Engineering Fracture Mechanics. 2017. No. 22. Pp. 146-167.
13. Muralidhara S., Raghu Prasad B.K., Eskandari H., Karihaloo B.L. Fracture process zone size and true fracture energy of concrete using acoustic emission // Construction and Building Materials. 2010. No. 24. Pp. 479-486.
14. Muralidhara S., Raghu Prasad B.K., Karihaloo B.L., Singh R.K. Size-independent fracture energy in plain concrete beams using tri-linear model // Construction and Building Materials. 2011. № 25. Pp. 3051-3058.
15. Shafigullin L.N., Bobrishev A.A., Erofeev V.T., Treshchev A.A., Shafigullina A.N. Development of the recommendations on selection of glass-fiber reinforced poly-urethanes for vehicle parts // International Journal of Applied Engineering Reserch. 2015. Vol. 10. No. 23. Pp. 43758-43762.
16. Shafigullin L.N., Treshchev A.A., Hodorovich P.Y., Erofeev V. T. The Stress-Strain State Of Layered Ortho-tropic Conditional Half-Space Taking Into Account Different Resistance // Revista Publicando. 2017. Vol. 4. No. 13(2). Pp. 109-127.
17. Кристенсен Р.М. Введение в механику композитов. М.: Мир, 1982. 336 с.
18. Черкасов В.Д., Тюрякин А.С. Теория двухсвязных моделей микромеханики композитов: монография. Саранск: Изд-во Мордов. ун-та, 2009. 108 с.
19. Hashin Z. The elastic moduli of heterogeneous materials // J. Appl. Mech. 1962. Vol. 29. No. 1. Pp. 143-150.
20. Reuss A. Berechung der Fliessgrenze von Mischkristallen auf Grund der Plastizitatsbedingund // Z. Angew. Math. und Mech. 1929. Vol. 9. No. 1. Pp. 49-58.
21. Voigt W. Lehrbuch der Kristallphysik. Berlin: Teubner, 1928. 962 p.
RESEARCH PAPER
Effective modules of two-phase construction composites with grain filler
Vladimir T. Erofeev, Aleksej S. Tyuryahin, Tatyana P. Tyuryahina, Aleksandr V. Tingaev*
National Research Ogarev Mordovia State University (National Research University), 68 Bolshevistskaya St., Saransk, 430005, Russian Federation *[email protected]
Article history: Abstract
Received: October 19, 2019 In the book of R.M. Christensen, "Introduction to the Mechanics of Compo-
Revised: December 02, 2019 sites" (1982), a calculation formula is given for the bulk module of polydisperse
Accepted: December 06, 2019 composites with spherical inclusions. This formula has been known to the Russian-
speaking reader for almost 40 years, but unfortunately, it is not used in the practice of building materials science. To identify applied possibilities, R.M. Christensen's formula is modified and reduced to a dimensionless function k = k (w, n, 9), which depends on three dimensionless parameters, i.e., it depends on three quantities: w is the volume fraction of the inclusion, n - the ratio of the shear modulus of the matrix material to the volume modulus of the same matrix, 9 is the ratio of the volume moduli of the matrix materials and inclusion. Numerical studies of this function reveal that in two-phase granular composites, the range of effective moduli is significantly narrowed compared to the region limited by Voigt and Reuss estimates (in the sense of the upper and lower bounds of real values). At the same time, the lower Christensen score is the same as the Reuss score. Numerical and graphically presented results are given on the examples of the study of two characteristic groups of composite materials. In addition, the dimensionless form of the effective module allows to construct a system of visual graphic dependencies of the functions k(w) in a flat space k - w. For different values of 9, the function k = k (w, n) displays a bunch of curved segments, which sets the position of the plane figure in flat space. Examples of constructing figures for characteristic regions of the values of the function k (n, 9, w) are given.
Keywords: two-phase model of granular composite; spherical shape of the phases of the matrix and aggregate; effective bulk modulus of elasticity of the composite; Voigt - Reuss plug for an effective module
For citation
Erofeev V.T., Tyuryahin A. S., Tyuryahina T.P., Tingaev A.V. (2019). Effective modules of two-phase construction composites with grain filler. Structural Mechanics of Engineering Constructions and Buildings, 15(6), 407-414. http://dx.doi.org/10.22363/1815-5235-2019-15-6-407-414. (In Russ.)
References
1. Bezukhov N.I. (1968). Osnovy teorii uprugosti, plas-tichnosti i polzuchesti [Fundamentals of the theory of elasticity, plasticity and creep]. Moscow, Vysshaya shkola Publ. (In Russ.)
2. Bobryshev A.N., Erofeev V.T., Kozomazov V.M. (2012). Fizika i sinergetika dispersno-neuporyadochennyh kondensirovannyh kompozitnyh sistem [Physics and synergetics of dispersively disordered condensed composite systems]. Saint Petersburg, Nauka Publ. (In Russ.)
3. Vasil'ev V.V., Protasov V.V., Bolotin V.V. (1990). Kompozitnye materialy [Composite materials]: reference book. Moscow, Mashinostroenie Publ. (In Russ.)
4. Gusev B.V., Kondrashenko V.I., Maslov B.P., Fay-sovich A.S. (2006). Formirovanie struktury kompozicion-nyh materialov i ih svojstva [Formation of the structure of
Vladimir T. Erofeev, Doctor of Technical Sciences, Professor, Dean of the Faculty of Architecture and Construction; Academician of the Russian Academy of Architecture and Construction Sciences; Scopus Author ID: 57080123300, eLIBRARY SPIN-code: 4425-5045. Aleksej S. Tyuryahin, Doctor of Engineering, Associate Professor of Department of Applied Mechanics.
Tatyana P. Tyuryahina, graduate student of Faculty of Architecture and Civil Engineering.
Aleksandr V. Tingaev, undergraduate of Faculty of Architecture and Civil Engineering.
composite materials and their properties]. Moscow, Nauchnyj mir Publ. (In Russ.)
5. Erofeev V.T., Tyuryakhin A.S., Erofeeva I.V. (2018). On the connections of the carrier parameters with effective parameters in models of grain composites. Structural Mechanics and Analysis of Constructions, (3), 7-17. (In Russ.)
6. Erofeev V., Tyuryakhin A., Tyuryakhina T. (2019). Flat space of values of volume module of grain composite with spherical fill-lem. International Journal of Civial Engineering and Technology (IJCIET), (8), 333-342.
7. Berlin A.A., Vol'fson S.A., Oshmyan V.G., Eni-kolopov N.S. (1990). Principy sozdaniya kompozitnyh poli-mernyh materialov [Principles of creating composite polymer materials]. Moscow, Himiya Publ. (In Russ.)
8. Askadskij A.A., Goleneva L.M., Bychko K.A., Ka-zanceva V.V., Konstantinov K.V., Almaeva E.S., Klins-kih A.F., Kovriga O.V. (2001). Gradientnye polimernye materialy [Gradient materials]. Russian Chemical Journal, 45(3), 123-128. (In Russ.)
9. Duan K., Hu Xiao, Wittmann F.H. (2003). Boundary effect on concrete fracture and non-constant fracture energy distribution. Fracture Mechanics, (70), 2257-2268.
10. Carpinteri A., Chiaia B., Cornetti P. (2003). On the mechanics of quasi-brittle materials with a fractal microstructure. Engineering Fracture Mechanics, (70), 2321-2349.
11. Ayatollahi M.R., Akbardoost J. (2012). Size effects on fracture toughness of quasi-brittle materials - a new approach. Engineering Fracture Mechanics, (92), 89-100.
12. Hu X., Guan J., Wang Y., Keating A., Yang S. (2017). Comparison of boundary and size effect models based on new developments. Engineering Fracture Mechanics, (22), 146-167.
13. Muralidhara S., Raghu Prasad B.K., Eskandari H., Karihaloo B.L. (2010). Fracture process zone size and true fracture energy of concrete using acoustic emission. Construction and Building Materials, (24), 479-486.
14. Muralidhara S., Raghu Prasad B.K., Karihaloo B.L., Singh R.K. (2011). Size-independent fracture energy in plain concrete beams using tri-linear model. Construction and Building Materials, (25), 3051-3058.
15. Shafigullin L.N., Bobrishev A.A., Erofeev V.T., Treshchev A.A., Shafigullina A.N. (2015). Development of the recommendations on selection of glass-fiber reinforced polyurethanes for vehicle parts. International Journal of Applied Engineering Research, 210(23), 43758-43762.
16. Shafigullin L.N., Treshchev A.A., Hodorovich P.Y., Erofeev V.T. (2017). The Stress-Strain State of Layered Orthotopic Conditional Half-Space Taking into Account Different Resistance. Revista Publicando, 4 (13-2), 109-127.
17. Christensen R.M. (1982). Vvedenie v mekhaniku kompozitov [Introduction to the mechanics of composites]. Moscow, Mir Publ. (In Russ.)
18. Cherkasov V.D., Tyuryakhin A.S. (2009). Teoriya dvuhsvyaznyh modelej mikromekhaniki kompozitov [The theory of biconnected models of micromechanics of composites]: monograph. Saransk, Publishing House of Mordovia University. (In Russ.)
19. Hashin Z. (1962). The elastic moduli of heterogeneous materials. J. Appl. Mech., 29(1), 143-150.
20. Reuss A., Angew Z. (1929) Berechung der Fliessgrenze von Mischkristallen auf Grund der Plastizitatsbedingund. Math. und Mech., 9(1), 49-58.
21. Voigt W. (1928). Lehrbuch der Kristallphysik. Berlin, Teubner.