Method of calculation of elastic effective properties of two-phase polydisperse media using multipoint statistical descriptors and the integral equations technique

Abstract


The aim of this study is to develop a new analytical approach to calculation of effective properties of elastic heterogeneous media based on multipoint approximations of solutions of stochastic boundary value problems. Prediction of macroscopic properties of heterogeneous media is associated with the need for a reliable description of their microstructural behavior, including the interaction between individual components. A number of analytical and numerical approaches have been developed to evaluate the effective properties of structurally inhomogeneous media. However none of them makes it possible to calculate the effective properties of such media with an absolute accuracy. One of the main limitations is imposed by taking into account the features of the microstructure of the medium, such as orientation, size, shape and distribution of inclusions, as well as the features of influence of the matrix on inclusions. In this paper, multipoint approximations of solutions of boundary value problems are used to calculate the effective characteristics, in which it is necessary to employ higher-order correlation functions, which allows to take into account the multiparticle interaction of microstructural elements with a higher extent. Analytical expressions for the calculation of the effective properties of structurally inhomogeneous media using multipoint higher-order approximations of solutions of stochastic boundary value problems in elastic formulations are obtained. A numerical comparison of the calculation results of relative effective characteristics of porous inhomogeneous polydisperse media with spherical inclusions of different volume fractions is performed. For the numerical solution of the integral-differential equations, a global adaptive strategy is applied in conjunction with the multidimensional integration rule and the IMT variable transformation rule to handle the singularity of the Green’s function. Some conclusions are made on the effectiveness and limitations of the proposed approach.

Full Text

Введение Многие современные материалы, используемые в высокотехнологичных приложениях, обладают неоднородной микроструктурой. Неоднородность структуры материала может быть определена целью его применения либо являться результатом возникновения нежелательных дефектов, обусловленных технологическими особенностями процесса производства. Установлено, что влияние на эффективные характеристики материала оказывает не только объемная доля содержащихся в нем компонентов, но и морфологическое строение микроструктуры. В связи с этим актуальной проблемой механики неоднородных сред является установление взаимосвязи между конфигурацией микроструктуры и механическими и физическими свойствами, реализованными в материале. В частности, необходимо развитие методов математического моделирования, способных оценить поведение неоднородных материалов на различных масштабных уровнях. Предсказание макроскопических свойств композиционных материалов сопряжено с необходимостью достоверного описания их микроструктурного поведения, включая взаимодействия между отдельными компонентами. В целом для создания эффективной методики с целью определения макроскопических свойств структурно-неоднородных сред необходимы: - геометрическая модель микроструктуры, для которой можно в явном виде произвести численный расчет с учетом ее геометрии и свойств компонент; - механическая модель, с достаточной точностью описывающая взаимодействие распределенных компонент и окружающей их матрицы. К настоящему времени разработан ряд аналитических и численных подходов для оценки эффективных свойств структурно-неоднородных сред. Одним из основополагающих исследований в области определения эффективных свойств неоднородных сред является решение для единственного эллипсоидального включения, полученное в работах Эшелби [1]. Границы эффективных свойств упругих неоднородных сред на основе вариационных принципов были получены в работах Фойгта и Рейсса, Хашина и Штрикмана [2], Берана [3]. Многие подходы базируются на методах среднего поля, впервые появившихся в середине 20-го века. Эти методы предлагают ряд алгоритмов для определения упругих констант, например метод Мори-Танака [4], метод двойного включения [5], а также другие [6-8]. Методы гомогенизации для нелинейных композитов берут начало в работах Тейлора и Хилла [9], в которых рассматривалась проблема определения макроскопического поведения пластически деформируемых поликристаллов. Самосогласованные модели применительно к упругопластическим композитам были предложены Е. Кренером [10], Б. Будянски [11] и Р. Хиллом [9], в которых предполагалось постоянство пластических деформаций. Позднее Дж.Р. Уиллисом [12] было получено нелинейное обобщение границ Хашина-Штрикмана. П. Понте-Кастанедой [13, 14] были предложены новые вариационные подходы, позволяющие улучшить оценки границ и получить более общие решения для нелинейных сред. Некоторые авторы предлагали модификации схем Мори-Танака для учета произвольно ориентированных эллипсоидальных включений [15], а также для волокон с заданной функцией распределения ориентации [16], в том числе их реализацию для упругопластического случая [17]. С развитием вычислительных возможностей появились работы, позволяющие производить расчет эффективных модулей для представительного объема путем прямого решения методом конечных элементов [18-20]. Несмотря на достигнутые успехи, прямой конечно-элементный расчет представительных объемов для определения эффективных свойств сильно зависит от качества и густоты сетки и достаточно затратен в плане вычислительных ресурсов при исследовании структур со случайной ориентацией и сложной геометрией включений. В отдельную задачу входит определение минимально допустимого размера представительного объема [21-25]. В работах Севостьянова и Качанова [26-28], а также Бема [29-32] и других [33-35] установлено влияние формы неоднородностей на макроскопические свойства. Другая группа аналитических моделей для определения макроскопического отклика неоднородных материалов базируется на использовании многоточечных статистических дескрипторов. Трехточечные оценки термоупругих свойств двухфазных материалов на основе набора многоточечных корреляций были получены в работах С. Торквато [36,37], при этом удалось получить границы, находящиеся внутри границ Хашина-Штрикмана. В отечественных работах в данном направлении широко распространен подход, получивший развитие в работах Т.Д. Шермергора и С.Д. Волкова [38, 39] и позволяющий определять макроскопические модули, а также структурные поля деформаций и напряжений в композитах стохастической структуры. Так, для исследования структурных полей деформирования и расчета эффективных характеристик микронеоднородных сред решаются стохастические краевые задачи, уравнения и граничные условия которых содержат случайные величины. Исходное стохастическое уравнение с помощью функций Грина преобразуется в интегродифференциальное, которое решается приближенно методом итераций. Результаты приближенного решения зависят в первую очередь от вида координатной зависимости статистических дескрипторов структурных модулей упругости. При этом используются различные предположения: - сингулярное приближение, которое заключается в том, что в интегральных уравнениях равновесия, ядрами которых являются вторые производные тензора Грина для изотропной неограниченной среды, удерживаются только сингулярные составляющие этих производных [38]. С помощью обобщенного сингулярного приближения, изменяя параметры, как частный случай можно получать различные приближенные решения для эффективных упругих констант (границы Фойгта, Рейса, а также Хашина-Штрикмана); - гипотеза сильной изотропии в том смысле, что макроскопические характеристики композита не зависят от многоточечных моментных функций структурных модулей упругости [40-42]; - предельная локальность моментных функций структурных модулей упругости: ограничиваются учетом дисперсий моментных функций [43]; - корреляционное приближение основано на пренебрежении моментами выше второго порядка. В общем случае корреляционная теория может быть применена при малых флуктуациях структурных модулей упругости, т.е. когда среднеквадратичные отклонения структурных модулей упругости являются малыми по отношению к их математическим ожиданиям [39, 43]; - одноточечное приближение является уточнением корреляционного приближения за счет использования одноточечных моментов всех порядков. С помощью одноточечного приближения нельзя описывать анизотропию свойств композиционных материалов, связанную с ориентацией включений, например, матрицы с ориентированными эллипсоидальными включениями. Таким образом, в большинстве работ инструментарий описания стохастических структур ограничивается, как правило, статистическими дескрипторами второго порядка, которые во многих случаях принимаются равными постоянным значениям. В данной работе для расчета эффективных характеристик предлагается использование многоточечных приближений решений краевых задач, в которых необходимым является учет в явном виде моментных функций высших порядков, что позволит в более полной мере учитывать многочастичное взаимодействие микроструктурных элементов. Использование высших приближений решений даст возможность уточнить границы, полученные с использованием описанных выше допущений. 1. Аналитические выражения для эффективных констант Будем считать, что микроструктурные компоненты двухфазной среды являются изотропными, т.е. (1) где нижний индекс I обозначает включения; M - матрицу; - тензоры структурных модулей упругости; - дельта Дирака; - постоянные Ламе и модули сдвига включений и матрицы соответственно, выражаемые через модуль упругости и коэффициент Пуассона как (2) где - модуль упругости включения; - коэффициент Пуассона. Аналогично для и . Макроскопические деформации и напряжения определяются путем осреднения по элементарному макрообъему : (3) Поля микроструктурных напряжений и деформаций являются случайными, зависят от пространственного радиуса-вектора и представляют собой кусочно-постоянные функции координат. При моделировании неоднородных структур со случайным разбросом включений, как правило, принимается гипотеза эргодичности, согласно которой осреднение по объему совпадает со статистическим осреднением: (4) В общем случае, когда , а - объемная доля включений, эффективные константы упругой неоднородной среды могут быть получены из решения системы линейных уравнений (5) где - тензор эффективных модулей упругости, компоненты которого находятся по формуле ; - компоненты постоянного заданного симметричного тензора малых макродеформаций, определяющего нагрузку представительного объема. В случае пористых материалов последнее уравнение сводится к следующей зависимости: (6) При вычислении структурных деформаций и напряжений используют явные представления поля через постоянные тензоры модулей упругости элементов структуры (1) и случайную индикаторную функцию , которая принимает значения 0 или 1 в зависимости от положения радиуса-вектора внутри представительного объема: (7) где и - тензоры модулей упругости включений и матрицы соответственно. Выражения для средних значений напряжений и деформаций в матрице могут быть получены на основе формул для условного математического ожидания, представленных в [39]: (8) (9) где , , и являются флуктуациями полей случайной индикаторной функции, напряжений, деформаций и перемещений соответственно. В общем случае флуктуации случайной величины в точке определяются как , где - оператор статистического осреднения. Выражение для статистического момента второго порядка , содержащего флуктуации индикаторной функции и поля напряжений, выражается следующим образом [44-46]: (10) где - разность тензоров структурных модулей упругости включений и матрицы; - центральный момент второго порядка (дисперсия) случайного поля индикаторной функции, выражаемый через объемную долю включений. Тогда, с учетом соотношения (5), общее выражение для компонент тензора эффективных модулей упругости принимает следующий вид: (11) (12) Считая макрооднородный представительный объем неоднородной среды изотропным, достаточно определить только две его упругие константы. Так, учитывая, что для изотропного тензора эффективных модулей упругости , значение эффективной константы можно получить из решения задачи о чистом сдвиге. (13) Значение второй эффективной константы находится из решения задачи о всестороннем растяжении: (14) (15) Таким образом, эффективные модули согласно предложенному методу определяются как сумма их средних значений и уточняющей добавки, связанной с упругими многочастичными взаимодействиями между неоднородностями. Для вывода аналитических выражений для смешанного статистического момента второго порядка , определяющего величину данной добавки, в данной работе предлагается использовать приближенные решения стохастической краевой задачи относительно перемещений, представленные в виде интегродифференциального уравнения и полученные на основе метода функций Грина. 2. Приближенные решения для эффективных констант на основе многоточечных статистик Как показано в работах [44-48], стохастическая краевая задача теории упругости для структурно-неоднородных сред с помощью функции Грина сводится к интегродифференциальному уравнению относительно флуктуаций перемещений: (16) где - функция Грина; - флуктуация поля структурных модулей упругости; - обозначение k-й компоненты радиуса-вектора ; - обозначение производной ; - приближение, в котором решается задача. Форма данного уравнения предполагает его итерационное решение. В первой и второй итерации решения уравнения (16) принимают следующую форму: (17) (18) где - флуктуация случайной индикаторной функции в точке; - поле флуктуаций перемещений в представительном объеме, которое является решением краевой задачи теории упругости в стохастической постановке. Выражения для момента с использованием решений (17) и (18) выражаются в следующем виде: (19) (20) Второе слагаемое в скобках в выражении (20) является добавкой, вносимой вторым приближением (18) к решению в первом приближении (17). Перемножение флуктуаций индикаторной функции приводит к появлению в выражениях в явном виде корреляционных и моментных функций, которые используются для оценки пространственного взаимодействия между микромасштабными структурными компонентами [49]. Моментная функция n-го порядка может быть определена как среднее скалярного произведения n случайных величин в различных положениях радиусов-векторов . Для флуктуации случайной индикаторной функции выражение для моментной функции в целом можно записать в виде (21) Моментные функции чувствительны к таким параметрам, как распределение, ориентация и форма микроструктурных составляющих. Теоретически микроструктурная морфология композитов может быть однозначно определена бесконечным числом моментных функций. Эти функции в том числе используются в качестве геометрических характеристик в задачах получения статистик высоких порядков для микроструктурных полей напряжений и деформаций [45]. Возвращаясь к расчету эффективных констант, в явном виде значение смешанного момента в выражении для эффективной константы (13) с использованием первого приближения (19) и с учетом свойств симметрии входящих в него тензоров записываем как (22) После преобразования получаем (23) Добавка, вносимая вторым приближением в (20), выражается следующим образом: (24) Таким образом, эффективная константа может быть определена как сумма среднего по представительному объему значения компоненты тензора , а также двух слагаемых, которые являются уточнениями, вносимыми первым и вторым приближением решения уравнения (16) соответственно: (25) Для определения эффективной константы запишем в явном виде момент : (26) После упрощений данное выражение приводится к виду (27) Аналогичным образом находится добавка, соответствующая второму приближению решения: (28) В окончательно виде формула для эффективной константы принимает вид (29) Таким образом, особенностью метода последовательных приближений является возможность уточнения значения эффективных констант в зависимости от морфологических особенностей, выражаемых через структурные моментные функции. Из выражений (19) и (20) видно, что каждое последующее приближение зависит от структурной моментной функции более высокого порядка, что позволяет в более полной степени учитывать морфологию микроструктуры в необходимых случаях. 3. Численная реализация метода Для численного решения интегралов, входящих в полученные уравнения, необходим явный вид функции Грина. Когда размеры представительного объема много меньше характерных размеров твердого тела, в случае изотропной среды в качестве функции Грина для трехмерного пространства в работах [10, 38] было получено следующее выражение: (30) При этом Для вычисления интегралов, входящих в выражения (25) и (29), необходимо использовать методы численного интегрирования. Интегрирование проводится по всей области статистической зависимости в сферических координатах, где координаты вектора определяются по формулам , , , , , , где - сторона представительного объема. Переход к сферическим координатам обусловлен избавлением от сингулярности функции Грина (30) по всем трем переменным декартовой системы координат и переходом к сингулярности только по одной переменной . Для достижения допустимой погрешности интегрирования применяется глобальная адаптивная стратегия в совокупности с многомерным правилом интегрирования и правилом преобразования переменных IMT для обхода сингулярности. Глобальная адаптивная стратегия основана на том, что из всех подобластей исходной области интегрирования выбирается подобласть с наибольшей погрешностью и делится пополам. Затем вычисляются значение интеграла и погрешность для каждой половины, и процедура повторяется для всего множества подобластей. После каждого такого шага пересчитываются глобальное значение интеграла и глобальная погрешность, которые представляют собой сумму значения интеграла и погрешности в каждой подобласти. Ожидается, что глобальная погрешность должна монотонно снижаться с увеличением числа областей. Многомерное правило интегрирования [50] определяется следующим образом. Для куба , задается набор точек, удовлетворяющий следующим условиям: любая точка из множества может быть получена перемещением или заменой знака координаты любой другой точки из этого же множества; все точки множества имеют одинаковый весовой коэффициент. Множество точек, соответствующее вышеописанным критериям, называется орбитой. Если правило интегрирования имеет орбит а i-я орбита имеет весовой коэффициент , интеграл приближенно вычисляется по формуле (31) Преобразование переменных IMT предложено в работах Ири, Моригути и Токасава [51] и основано на идее преобразования независимой переменной таким образом, что все производные нового подынтегрального выражения обращаются в ноль в конечных точках интервала интегрирования. В приведенном ниже примере выполнено численное сравнение результатов расчета относительных эффективных характеристик пористых неоднородных полидисперсных сред со сферическими включениями. Пример геометрической модели такой среды представлен на рис. 1, где красным цветом выделены включения (поры). Создание модели осуществлялось методом последовательного синтеза, согласно которому каждое новое включение помещается в произвольную точку внутри куба, при этом контролируется его пересечение с уже существующими включениями. В случае наличия пересечения данное включение отбрасывается и генерируется новое. В качестве включений использовались сферы различных радиусов, проверка пересечения между которыми осуществлялась путем расчета расстояния между центрами. Рис. 1. Представительный объем полидисперсной неоднородной среды с объемной долей включений 40 % Fig. 1. Representative volume element of a polydisperse heterogeneous medium with 40 % volume fraction of inclusions Моментные функции, характеризующие неоднородную среду, получены на основе разбиения трехмерного представительного объема сеткой с фиксированным шагом и кубической формой конечного элемента (так называемая воксельная сетка). Преимущество использования такого способа получения значений моментных функций связано еще и с тем, что информацию о представительных объемах неоднородных материалов в таком виде легко получать из экспериментальных исследований, например, на основе микрокомпьютерной томографии. На рис. 2 изображены нормированные моментные функции второго и третьего порядка вида , где - дисперсия n-го порядка индикаторной функции , полученные для изображенной на рис. 1 среды и используемые в выражениях (25) и (29). Рис. 2. Моментные функции второго и третьего порядка для двухфазной неоднородной среды с объемной долей включений 40 % Fig. 2. Second and third order correlation functions for a two-phase heterogeneous medium with 40 % volume fraction of inclusions На рис. 3 представлены значения отношения эффективных модулей и среды к модулям матрицы. Данные величины являются безразмерными и изменяются от 0 до 1. Также приведено сравнение с границами Хашина-Штрикмана [2] и сингулярного приближения, предложенного Т.Д. Шермергором [38]. Границы Хашина-Штрикмана эффективных констант для двухфазной неоднородной среды вычисляются по формулам (32) (33) где константы с индексом 1 относятся к компоненту с более жесткими свойствами. Рис. 3. Значения эффективных характеристик среды, полученные в первом приближении, в сравнении с границами Хашина-Штрикмана и сингулярного приближения Fig. 3. Values of the medium effective characteristics obtained in the first approximation in comparison with the Hashin-Shtrikman and singular approximation bounds Сингулярное приближение основано на приведении решения к интегральному выражению, содержащему вторую производную функции Грина, и затем выделению сингулярной составляющей из нее, сводя таким образом интеграл по сингулярной области к постоянному тензору [38]. Согласно сингулярному приближению, точное значение эффективных констант находится в границах, образованных следующими неравенствами: (34) (35) где ; ; . Все изложенные результаты, включая численное интегрирование, получены в программной среде Wolfram Mathematica. Как видно из результатов, использование предложенной методики с применением первого приближения позволяет уточнить значения эффективных констант для пористых структур с объемной долей выше 20 % по сравнению с сингулярным приближением, где не используется в явном виде информация о морфологии. Установлено, что для пористых структур с равномерным распределением включений добавка, вносимая вторым приближением решения, является несущественной. В то же время значимость этой добавки увеличивается с ростом отклонения распределения включений в представительном объеме от равномерного, который может быть обусловлен формированием кластеров, сильным разбросом размеров включений и т.д. Заключение В данной работе предложен метод вычисления эффективных характеристик полидисперсных неоднородных сред на основе применения многоточечных моментных функций и метода интегральных уравнений в виде приближенных решений стохастических краевых задач теории упругости. Входными параметрами данной модели являются упругие свойства компонент среды, а также формализованная в виде многоточечных моментных функций информация о морфологии среды. Полученные в явном виде выражения для эффективных констант макрооднородной изотропной среды представлены в виде средней составляющей и добавки, точность которой определяется приближением используемого решения и порядком входящих в него моментных функций. Интегродифференциальные уравнения, которые входят в добавку, могут быть решены численно с использованием многомерных методов интегрирования и правил обращения с сингулярностью функции Грина. Важным фактором, влияющим на работоспособность представленного подхода при использовании приближений высоких порядков, является точность используемых методов численного интегрирования. Это связано с тем, что при росте размерности интегралов в связи с увеличением порядка приближения возрастает значение фактора сингулярности производных функций Грина. При перемножении таких производных требуется тонкая отладка методов и правил интегрирования для «отсечения» зоны сингулярности без влияния на погрешность интегрирования. Так, для исследованных примеров пористых сред установлено, что использование второго приближения с описанными выше способами численного интегрирования приводит к недопустимой погрешности при объемной доле включений выше 50 %. Это связано с тем, что при увеличении объемной доли точка пересечения с осью графика абсцис моментных функций (рис. 2) смещается вправо, тем самым увеличивая скорость возрастания значений подынтегральной функции в области, близкой к сингулярности, тем самым увеличивая погрешность интегрирования. Таким образом, для получения надежных результатов использования методики с многоточечными статистиками высших порядков в качестве инструментов описания морфологии среды необходимо дополнительное исследование подходящих способов численного интегрирования.

About the authors

M A Tashkinov

Perm National Research Polytechnic University

References

  1. Eshelby J.D. The determination of the elastic field of an ellipsoidal inclusion, and related problems // Collected Works of J.D. Eshelby. - 2007. doi: 10.1007/1-4020-4499-2_18
  2. Hashin Z., Shtrikman S. A Variational approach to the theory of the effective magnetic permeability of multiphase materials // J. Appl. Phys. - 1962. - Vol. 33. - No. 10. - P. 3125-3131. doi: 10.1063/1.1728579
  3. Beran M. Statistical continuum theories // J. Rheol. (N. Y. N. Y). - 1965. - Vol. 9. - No. 1. - P. 339-355. doi: 10.1122/1.548991
  4. Mori T., Tanaka K. Average stress in matrix and average elastic energy of materials with misfitting inclusions // Acta Metall. - 1973. - Vol. 21. - No. 5. - P. 571-574. doi: 10.1016/0001-6160(73)90064-3
  5. Hori M., Nemat-Nasser S. Double-inclusion model and overall moduli of multi-phase composites // Mech. Mater. - 1993. - Vol. 14. - No. 3. - P. 189-206. doi: 10.1016/0167-6636(93)90066-Z
  6. Christensen R.M. Mechanics of composite materials. - New York: Wiley-Interscience, 1979. - 348 p.
  7. Berryman J.G. Measures of microstructure to improve estimates and bounds on elastic constants and transport coefficients in heterogeneous media // Mech. Mater. - 2006. - Vol. 38. - No. 8-10. - P. 732-747. doi: 10.1016/j.mechmat.2005.06.014
  8. Pierard O., Friebel C., Doghri I. Mean-field homogenization of multi-phase thermo-elastic composites: A general framework and its validation // Compos. Sci. Technol. - 2004. - Vol. 64. - № 10-11. - P. 1587-1603. doi: 10.1016/j.compscitech.2003.11.009
  9. Hill R. A self-consistent mechanics of composite materials // J. Mech. Phys. Solids. - 1965. - Vol. 13. - No. 4. - P. 213-222. doi: 10.1016/0022-5096(65)90010-4
  10. Kröner E. Statistical Modelling // Modelling Small Deformations of Polycrystals. - Dordrecht: Springer Netherlands, 1986. - P. 229-291. doi: 10.1007/978-94-009-4181-6_8
  11. Budiansky B. On the elastic moduli of some heterogeneous materials // J. Mech. Phys. Solids. - 1965. - Vol. 13. - P. 223-227. doi: 10.1016/0022-5096(65)90011-6
  12. Willis J.R. Variational and Related Methods for the Overall Properties of Composites // Advances in Applied Mechanics. - 1981. - Vol. 21. - No. C. - P. 1-78. doi: 10.1016/S0065-2156(08)70330-2
  13. Castañeda P.P. Exact second-order estimates for the effective mechanical properties of nonlinear composite materials // J. Mech. Phys. Solids. - 1996. - Vol. 44. - No. 6. - P. 827-862. doi: 10.1016/0022-5096(96)00015-4
  14. Castañeda P.P. New variational principles in plasticity and their application to composite materials // J. Mech. Phys. Solids. - 1992. - Vol. 40. - No. 8. - P. 1757-1788. doi: 10.1016/0022-5096(92)90050-C.
  15. Benveniste Y. A new approach to the application of Mori-Tanaka’s theory in composite materials // Mech. Mater. - 1987. - Vol. 6. - No. 2. - P. 147-157. doi: 10.1016/0167-6636(87)90005-6
  16. Allen D.H., Lee J.-W. The Effective Thermoelastic Properties of Whisker- Reinforced Composites as Functions of Material Forming Parameters // Micromechanics and Inhomogeneity. - New York: Springer New York, 1990. - P. 17-39. doi: 10.1007/978-1-4613-8919-4_2
  17. Dunn M.L., Ledbetter H. Elastic-plastic behavior of textured short-fiber composites // Acta Mater. - 1997. - Vol. 45. - No. 8. - P. 3327-3340. doi: 10.1016/S1359-6454(96)00401-6
  18. Computation of linear elastic properties from microtomographic images: Methodology and agreement between theory and experiment / C.H. Arns [et al.] // GEOPHYSICS. - 2002. - Vol. 67. - No. 5. - P. 1396-1405. doi: 10.1190/1.1512785
  19. Comparison of prediction on effective elastic property and shape optimization of truss material with periodic microstructure / J. Yan [et al.] // Int. J. Mech. Sci. - 2006. - Vol. 48. - No. 4. - P. 400-413. doi: 10.1016/j.ijmecsci.2005.11.003
  20. Ostoja-Starzewski M. Material spatial randomness: From statistical to representative volume element // Probabilistic Eng. Mech. - 2006. - Vol. 21. - No. 2. - P. 112-132. doi: 10.1016/j.probengmech.2005.07.007
  21. Drugan W.J., Willis J.R. A micromechanics-based nonlocal constitutive equation and estimates of representative volume element size for elastic composites // J. Mech. Phys. Solids. - 1996. - Vol. 44. - No. 4. - P. 497-524. doi: 10.1016/0022-5096(96)00007-5
  22. Monetto I., Drugan W.J. A micromechanics-based nonlocal constitutive equation and minimum RVE size estimates for random elastic composites containing aligned spheroidal heterogeneities // J. Mech. Phys. Solids. - 2009. - Vol. 57. - No. 9. - P. 1578-1595. doi: 10.1016/j.jmps.2009.05.005
  23. Determination of the size of the representative volume element for random composites: Statistical and numerical approach / T. Kanit [et al.] // Int. J. Solids Struct. - 2003. - Vol. 40. - No. 13-14. - P. 3647-3679. doi: 10.1016/S0020-7683(03)00143-4
  24. Khisaeva Z.F., Ostoja-Starzewski M. On the size of RVE in finite elasticity of random composites // J. Elast. - 2006. - Vol. 85. - No. 2. - P. 153-173. doi: 10.1007/s10659-006-9076-y
  25. Determination of the critical size of a statistical representative volume element (SRVE) for carbon reinforced polymers / D. Trias [et al.] // Acta Mater. - 2006. - Vol. 54. - No. 13. - P. 3471-3484. doi: 10.1016/j.actamat.2006.03.042
  26. Sevostianov I., Kachanov M., Zohdi T. On computation of the compliance and stiffness contribution tensors of non ellipsoidal inhomogeneities // Int. J. Solids Struct. - 2008. - Vol. 45. - No. 16. - P. 4375-4383. doi: 10.1016/j.ijsolstr.2008.03.020
  27. Kachanov M., Tsukrov I., Shafiro B. Effective Moduli of Solids With Cavities of Various Shapes // Appl. Mech. Rev. - 1994. - Vol. 47. - No. 1S. - P. S151. doi: 10.1115/1.3122810
  28. Evaluation of the effective elastic and conductive properties of a material containing concave pores / F. Chen [et al.] // Int. J. Eng. Sci. - 2015. - Vol. 97. - P. 60-68. doi: 10.1016/j.ijengsci.2015.08.012
  29. Influence of waviness and curliness of fibres on mechanical properties of composites / A.Y. Matveeva [et al.] // Comput. Mater. Sci. Elsevier B.V. - 2014. - Vol. 87. - P. 1-11. doi: 10.1016/j.commatsci.2014.01.061
  30. Böhm H.J., Rasool A. Effects of particle shape on the thermoelastoplastic behavior of particle reinforced composites // Int. J. Solids Struct. Elsevier Ltd. - 2016. - Vol. 87. - P. 90-101. doi: 10.1016/j.ijsolstr.2016.02.028
  31. Rasool A., Böhm H.J. Effects of particle shape on the macroscopic and microscopic linear behaviors of particle reinforced composites // Int. J. Eng. Sci. Elsevier Ltd. - 2012. - Vol. 58. - P. 21-34. doi: 10.1016/j.ijengsci.2012.03.022
  32. Klusemann B., Böhm H.J., Svendsen B. Homogenization methods for multi-phase elastic composites with non-elliptical reinforcements: Comparisons and benchmarks // Eur. J. Mech. A/Solids. Elsevier Masson SAS. - 2012. - Vol. 34. - P. 21-37. doi: 10.1016/j.euromechsol.2011.12.002
  33. Chen F.C., Young K. Inclusions of arbitrary shape in an elastic medium // J. Math. Phys. American Institute of Physics. - 1977. - Vol. 18. - No. 7. - P. 1412-1416. doi: 10.1063/1.523438
  34. Composite structure modeling and mechanical behavior of particle reinforced metal matrix composites / Y. Su [et al.] // Mater. Sci. Eng. A. - 2014. - Vol. 597. - P. 359-369. doi: 10.1016/j.msea.2014.01.024
  35. Automated modeling of random inclusion composites / M. Bailakanavar [et al.] // Eng. Comput. - 2012. - Vol. 30. - No. 4. - P. 609-625. doi: 10.1007/s00366-012-0310-x
  36. Torquato S. Random Heterogeneous Media: Microstructure and Improved Bounds on Effective Properties // Appl. Mech. Rev. - 1991. - Vol. 44. - No. 2. - P. 37. doi: 10.1115/1.3119494
  37. Torquato S. Morphology and effective properties of disordered heterogeneous media // Int. J. Solids Struct. Elsevier Science Ltd. - 1998. - Vol. 35. - No. 19. - P. 2385-2406. doi: 10.1016/S0020-7683(97)00142-X
  38. Шермергор Т.Д. Теория упругости микронеоднородных сред. - М.: Наука, 1977.
  39. Волков С.Д., Ставров В.П. Статистическая механика композитных материалов. - Минск: Изд-во Белорус. гос. ун-та, 1978. - 208 с.
  40. Ломакин В.А. Статистические задачи механики твердых деформируемых тел. - М.: Наука, 1970. - 139 с.
  41. Хорошун Л.П. Методы случайных функций в задачах о макроскопических свойствах микронеоднородных сред // Прикл. механика. - 1978. - Т. 14, № 2. - С. 3-17.
  42. Болотин В.В., Москаленко В. К расчету макроскопических постоянных сильно изотропных композиционных материалов // Механика твердого тела. - 1967. - Т. 3. - С. 106-111.
  43. Соколкин Ю.В., Ташкинов А.А. Механика деформирования и разрушения структурно неоднородных тел. - М.: Наука, 1984. - 116 с.
  44. Tashkinov M., Wildemann V., Mikhailova N. Method of successive approximations in stochastic elastic boundary value problem for structurally heterogenous materials // Comput. Mater. Sci. - 2012. - Vol. 52. - No. 1. - P. 101-106. doi: 10.1016/j.commatsci.2011.04.025
  45. Tashkinov M.A. Methods of Stochastic Mechanics for Characterization of Deformation in Randomly Reinforced Composite Materials // Mechanics of Advanced Materials. Ed. Silberschmidt V.V., Matveenko V.P. - Cham: Springer, 2015. - С. 43-78. doi: 10.1007/978-3-319-17118-0_3
  46. Tashkinov M. Micro-scale modeling of phase-level elastic fields of SiC reinforced metal matrix multiphase composites using statistical approach // Comput. Mater. Sci. Elsevier. - 2016. - Vol. 116. - P. 113-121. doi: 10.1016/j.commatsci.2015.10.047
  47. Tashkinov M. Modeling of elastic behavior of multicomponent composite materials based on approximate solution of stochastic boundary value problems // PNRPU Mech. Bull. - 2015. - No. 3. - P. 165-181. doi: 10.15593/perm.mech/2015.3.12
  48. Tashkinov M.A. Stochastic modelling of deformation process in elastoplastic composites with randomly located inclusions using high order correlation functions // PNRPU Mech. Bull. Perm National Research Polytechnic University. - 2014. - Vol. 2014. - No. 3. - P. 163-185. doi: 10.15593/perm.mech/2014.3.09
  49. Torquato S. Random Heterogeneous Materials // Applied Mechanics Reviews. - 2002. - Vol. 16. - XXI. - 703 p. doi: 10.1007/978-1-4757-6355-3
  50. Genz A.C., Malik A.A. An Imbedded Family of Fully Symmetric Numerical Integration Rules // SIAM J. Numer. Anal. - 2005. doi: 10.1137/0720038
  51. Davis P.J., Rabinowitz P. Approximate Integration over Infinite Intervals // Methods of Numerical Integration. Elsevier. - 1984. - P. 199-270. doi: 10.1016/B978-0-12-206360-2.50009-1

Statistics

Views

Abstract - 120

PDF (Russian) - 89

Cited-By


PlumX


Copyright (c) 2019 Tashkinov M.A.

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

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies