Full Text

Введение Важной и актуальной проблемой в механикенеоднородных материаловявляется разработка моделей механики неоднородных средсо случайной структурой, учитывающих особенности реальной структуры композиционных материалов, процессы нелинейного деформирования, появление и развитие областей разрушения в структуре микронеоднородных материалов задолго до полного разрушения конструкций. Доказано, что влияние распределения включений играет решающую роль в макроскопическом поведении неоднородных сред [1-4]. Таким образом, создание точных математических моделей позволяет предсказывать эффективное механическое поведение с учетом особенностей микроструктуры, таких как размер, объемная доля, ориентация, дисперсия и кластеризация составляющих компонентов. Распространенным направлением в микромеханике материалов со случайной структурой являются методы статистической механики для систем многих частиц, в которых используются многоточечные статистики для напряжений, деформаций и полей смещений, для описания взаимодействия неоднородностей. Основы статистической механики были заложены в работах М. Берана, В.А. Ломакина, Т.Д. Шермегора [6, 8, 9, 11, 12] и в дальнейшем широкое распространение получили во многих научных школах [1, 5, 7, 10, 13, 14]. Представленные характеристики полей учитывают хаотичность расположения элементов в структуре, а также статистическую дисперсию свойств компонентов. Например, эффективные упругие свойства могут быть оценены на основе статистических моментов первого порядка полей напряжений и структурных моментных функций [7, 8, 13, 15]. Моменты высшего порядка используются для исследования распределения напряжений и деформаций в компонентах. Методы, связанные с описанием гетерогенных структур с помощью корреляционных функций, также используются для отождествления моделей микроструктуры с экспериментально анализируемыми образцами материалов [16-18]. Анализ деформирования и разрушения в рамках статистических подходов основан на разработке многоуровневых математических моделей и численных алгоритмов для решения задачи микромеханического описания процессов деформации представительного объема материала [5-8, 13, 14, 19]. Параметры таких процессов зависят от конкретных характеристик полей структурных напряжений и деформаций, которые могут быть определены из решения стохастических краевых задач, где уравнения и граничные условия содержат случайные величины. Целью данной работы является разработка аналитического инструментария для анализа полей микроструктурных напряжений и деформаций в многокомпонентных неоднородных средах, позволяющего учитывать геометрические и физико-механические свойства компонентов, на основе построения корреляционных функций полей напряжений и деформаций. В разработанных ранее подходах в качестве статистических характеристик полей деформирования использовались статистические моменты, которыеявляются константами для конкретного представительного объема при заданном нагружении [6-8, 20-28]. При этом для оценки пространственного взаимодействия между микроструктурными компонентами используется набор структурных моментных функций различных порядков.Разработанная в данной работе аналитическая модель позволяет поставить в соответствие структурным моментным функциям корреляционные функции полей напряжений и деформаций, представляющие собой моментные функции второго порядка. 1. Постановка и решение стохастических краевых задач неоднородных сред Аналитические выражения для корреляционных функций полей напряжений и деформаций находятся на основе решения стохастических краевых задач в упругой и упругопластической постановке. Стохастическая краевая задача теории упругости для композитов со случайной структурой в некоторой области (представительном объеме) записывается следующим образом: , , , где - тензор структурных модулей упругости; - компоненты симметричного тензора малых макродеформаций, с помощью которого определяется нагружение представительного объема; - радиус-вектор с компонентами - поверхность представительного объема. Между компонентами выполняются условия идеального контакта , . Стохастическая краевая задача упругости композитов с помощью функции Грина сводится к интегродифференциальному уравнению относительно пульсаций перемещений [6-8, 22, 24], которое в первом и втором приближении для двухфазных материалов имеют следующий вид: , где функция Грина; обозначение производной ; флуктуация случайной индикаторной функции в точке; - разность тензоров структурных модулей упругости включений и матрицы; - поле флуктуаций перемещений в представительном объеме, которое является решением краевой задачи теории упругости или упругопластичности в стохастической постановке. С учетом разложения полей деформаций и напряжений на среднюю составляющую и пульсацию для пульсаций структурных напряжений и деформаций справедливы следующие соотношения: ё , . Тогда получаем , . Здесь комбинации типа при последующем осреднении образуют структурные моментные функции, значения которых могут быть получены для геометрии исследуемых представительных объемов. Выражения - используются для получения аналитического представления корреляционных функций полей напряжений и деформаций. Использование второго приближения предполагает использование структурных моментных функций более высокого порядка, при этом значительно увеличивается время расчета. Установлено, что для исследуемого класса пористых композитов вклад, вносимый в результаты вторым приближением, является незначительным (разница с результатами в первом приближении составляет 2-4%), поэтому с целью экономии вычислительных ресурсов представляется целесообразным использование первого приближения . 2. Аналитические выражения для многоточечных корреляционных функций полей напряжений и деформаций Напряженно-деформированное состояние представительного объема композита характеризуют моменты и моментные функции полей микроструктурных напряжений и деформаций. Существуют безусловные и условные статистические характеристики, предназначенные для описания механического поведения композита как макрооднородного материала и его компонентов в отдельности [7, 8]. Выражения для многоточечных условных и безусловных моментных функций стохастических полей напряжений и деформаций в представительных объемах композитов позволяют описывать взаимосвязь полей в удаленных на различное расстояние точках и при необходимости вычислять значения моментов случайных полей в произвольных точках представительного объема. В общем виде моментные функции n-го порядка полей напряжений могут быть получены из выражения (11) , , , (12) , , где оператор осреднения есть оператор математического ожидания. Для случайных однородных эргодических полей осреднение по объему совпадает со статистическим осреднением. - индикаторная функция для представительного объема, равная 1, если радиус-вектор находится в компоненте С, и 0 - в остальных случаях, - флуктуация индикаторной функции, величина определяет объемную долю компоненты в представительном объеме. В данных формулах случайные поля флуктуаций напряжений выражаются через поле структурных модулей упругости и флуктуации полей деформаций формулой (10). Для многокомпонентного композита в общем виде флуктуации структурных модулей упругости в i-м компоненте выражаются в следующем виде в упругом (13) и упругопластическом случае (14): (13) (14) где , - девиаторная часть тензора деформаций. Непосредственно получить решение упругопластической краевой задачи обычно не удается, так как связь тензоров напряжений и деформаций является физически нелинейной, и применить к ней традиционные методы механики затруднительно. Поэтому системе уравнений структурно-феноменологической модели, как и в упругой задаче, ставится в соответствие система уравнений для среды с эффективными свойствами. Чтобы воспользоваться статистическими методами, физические уравнения (которые записываются через определяющие соотношения в инвариантной форме: ) необходимо линеаризовать при определенных допущениях. Для решения поставленной задачи применяется вариант метода статистического осреднения, в пределах объема каждого компонента композита пренебрегаются пульсации деформаций, стоящие под знаком материальных функций С учетом формулы (10) и с помощью процедуры осреднения получены выражения для моментных функций полей напряжений, описывающих представительный объем как макрооднородный. Так, соотношение для безусловной корреляционной функции полей напряжений имеет вид (15) С учетом разложения полей деформаций и напряжений на средние составляющие и пульсации для многокомпонентной неоднородной среды корреляционные функции (моментные функции второго порядка) полей деформаций и напряжений в компоненте С записываются в следующем виде: (16) (17) Средние значения тензора деформаций и напряжений в компоненте С принимают вид [7, 8] (18) (19) Искомые моментные функции полей напряжений и полей деформаций как в компонентах неоднородной среды (16), (17), так и гомогенизированном представительном объеме (15) зависят от ряда безусловных и смешанных функций различных порядков, в которые входят флуктуации полей деформации , а также многоточечные структурные моментные функции которые содержат статистическую информацию о геометрии неоднородной структуры композитов. Для многокомпонентной среды структурные моментные функции строятся для каждой компоненты в отдельности: (20) Значения структурных моментных функций могут быть получены на основе анализа геометрической модели представительного объема [16, 18, 24]. В качестве примера далее приведены некоторые безусловные корреляционные функции полей деформаций, входящие в выражения для корреляционных функций полей напряжений и полученные с использованием первого и второго приближения решения стохастической краевой задачи в случае двухфазного композита: (21) (22) Таким образом, получены выражения для вычисления моментных функций первого и второго порядков полей напряжений и деформаций как в целом так и в компонентах композита с использованием структурных моментных функций высших порядков. 3. Трехмерные модели внутренней полидисперсной структуры текстурированных матричных композитов с эллипсоидальными включениями Для реализации методики были построены трехмерные модели внутренней полидисперсной структуры исследуемых материалов - двухкомпонентных текстурированных матричных композитов с эллипсоидальными включениями. В отличие от хаотически армированных композитов рассматривается текстурированный композит с эллипсоидальными включениями, имеющими определенную пространственную ориентацию. Предполагается, что включения расположены в матрице равномерно. Также было введено условие, согласно которому включения не могут пересекаться. Для создания геометрической модели был использован алгоритм, позволяющий синтезировать полидисперсные эллипсоидальные включения в представительном объеме заданного размера. Для проверки выполнения условия отсутствия пересечений эллипсоидов каждое включение описывается параллелепипедом, для которого проверяется пересечение с другими, уже существующими. При условии, что все эллипсоиды находятся в одной системе координат, необходимо отдельно сравнить минимумы и максимумы координат вершин параллелепипедов отдельно по трем осям. Так, если первый параллелепипед определен значениями координат , , , , , , а второй , , , , , , то эллипсоиды пересекаются, если выполняются все неравенства: , , , , , . На рис. 1приведены изображения некоторых полученных моделей представительных объемов с различной объемной долей эллипсоидальных включений. а б в Рис. 1. Представительные объемы структур с эллипсоидальными включениями с различной объемной долей: а - 5 %, количество включений 77; б - 10 %, количество включений 158; в - 15 %, количество включений 245 Fig. 1. Representative volumes of structures with ellipsoidal inclusionsof different volume fractions: а - 5 %, the number of inclusions is 77; b - 10 %, the number of inclusions is 158; c - 15 %, the number of inclusions is 245 Для исследуемых структур с эллипсоидальными включениями получены значения структурных моментных функций высших порядков, а также рассчитаны коэффициенты аналитических аппроксимирующих выражений для структурных моментных функций, необходимых для вычисления моментных функций полей напряжений и деформаций. В качестве выражений для аппроксимации использованы следующие зависимости [23, 27]: (23) (24) где использовано для аппроксимации функций второго порядка; - для функций порядка выше второго. Для построения моментных функций высшего порядка шаг между радиусами-векторами принимался равным. На рис. 2 отображены графики аппроксимированных моментных функций различных порядков для структуры с объемной долей включений 10 %. Рис. 2. Моментные функции различных порядков для структуры с объемной долей эллипсоидальных включений 10%: - моментная функций второго порядка; - моментная функция третьего порядка; - моментная функция четвертого порядка; - моментная функция пятого порядка Fig. 2. Correlation functions of different orders for a structure with a volume fraction of ellipsoidal inclusions amounting to 10 %: - is the second order correlation function; is the third order correlation function; - is the fourth order correlation function; - is the fifth order correlation function 4. Значения моментных функций полей деформирования в упругом случае при различных условиях нагружения и физико-механических характеристиках компонентов композитов Для получения значений условных (в компонентах композита) и безусловных (для представительного объема как целого) моментных функций полей напряжений и деформаций в представительных объемах и их компонентах использовано решение стохастической упругой и упругопластической краевой задачи в первом и втором приближении. В качестве входных параметров для вычислений выступают константы, характеризующие свойства компонентов представительного объема и входящие в функции Грина в подынтегральных выражениях, а также геометрические свойства компонентов, определяемые структурными моментными функциями. Нагружение представительного объема задается в виде граничных условий стохастической краевой задачи. Алгоритм численной реализации методики заключается в вычислении значений интегралов в выражениях для корреляционных функций при необходимых комбинациях индексов в зависимости от расстояния между радиусами-векторами и . При этом значения моментной функции находятся путем осреднения моментов для заданного расстоянии между радиусами-векторами , вычисленных в дискретных точках представительного объема V. Точность вычисления значений моментной функции определяется частотой (шагом) дискретизации представительного объема. Для реализации методик численного вычисления значений моментных функций полей деформирования в упругом и упругопластическом случае при различных условиях нагружения представительных объемов и физико-механических характеристиках компонентов созданы программные инструменты, реализованные в виде программного кода в среде Wolfram Mathematica, с использованием возможностей встроенных функций и параллельных вычислений, а также в среде С++. Используемые численные методики аналогичны методикам при вычислении моментов полей напряжений и деформаций, которые подробно рассмотрены в работах [24, 28]. На приведенных ниже рисунках представлены некоторые результаты для моментных функций полей деформаций в представительном объеме пористого композита. Представлены дискретные значения и их интерполяция. Для расчетов был выбран частный случай деформирования как чистый сдвиг. Для упругой краевой задачи был задан тензор макродеформаций на границе представительного объема в виде граничных условий (4), равный . Упругие свойства матрицы заданы следующими величинами: модуль упругости матрицы пористого композита = 2·105 МПа; коэффициент Пуассона = 0,3. На рис. 3-4 представлена нормированная корреляционная функция полей деформаций для пористого материала с различным содержанием сферических пор (р = 0,20 и р = 0,24). Нормирование проводилось делением значений построенной корреляционной функции на дисперсию полей деформаций, т.е. на корреляционную функцию при = 0. В качестве размеров используем условные единицы. Вид полученных моментных функций полей деформаций в целом схож с видом структурных моментных функций (см. рис. 2), поэтому для аппроксимации были протестированы те же исходные выражения (25), (26), что были использованы для аппроксимации структурных функций. Проведено исследование моментной функции полей деформаций для структур с эллипсоидальными порами от объемного содержания (рис. 5) и от вариации разброса максимального размера полуоси эллипсоидальных включений (рис. 6). Так как наибольшие значения корреляционной функции полей деформаций возникают при =0, а это дисперсии полей деформаций, то при решении краевой задачи в каких-то частных случаях достаточно учитывать только дисперсии для характеристики полей деформаций. Рис. 3. Нормированная интерполированная моментная функция деформации для структуры с объемной долей 20 % при чистом сдвиге; значения нормированной моментной функции деформаций; интерполированнаямоментнаяфункциядеформаций Fig. 3. Normalized interpolated correlation function of strain for a structure with a volume fraction of 20 % at pure share; are values of normalized correlation function of strains; is interpolated correlation function of strains Рис. 4. Нормированная моментная функция деформации для структуры с объемной долей 24 % при чистом сдвиге; значения нормированной моментной функции деформаций; интерполированнаямоментнаяфункциядеформаций Fig. 4. Normalized correlation function of strain for a structure with a volume fraction of 24 % at pure share; are values of normalized correlation function of strains; is interpolated correlation function of strains Рис. 5. Интерполированная ( ) и аппроксимированная ( ) моментные функции деформаций для структур с различным разбросом максимального размера полуоси эллипсоидальных включений: а - от 6 до 22; б - от 9 до 13; в - от 4 до 24 Fig. 5. Interpolated ( ) and approximated ( ) correlation function of strains for structures with a different range of maximum sizes of semi-axis ellipsoidal inclusions: а - is from 6 to 22; b - from 9 to 13; c - from 4 to 24 Рис. 6. Интерполированная ( ) и аппроксимированная ( ) моментные функции деформаций для структур с различной объемной долей включений: а - 4 %; б - 10 %; в - 13 % Fig. 6. Interpolated ( ) andapproximated ( ) correlation function of strains for structures with a different range of volume fraction of inclusions: а - 4 %; b - 10 %; c - 13 % В качестве примера решения упругопластической задачи был произведен расчет также для пористого материала с матрицей, обладающей следующими свойствами: модуль упрочнения ГПа, модуль сдвига ГПа, коэффициент объемного сжатия ГПа, . В качестве примера представлены моментные Ряд 3 Ряд 2 Ряд 1 Рис. 7. Нормированные моментные функции полей деформаций в упругопластической задаче в разных точках деформирования для структур со сферическими включениями с разным объемным содержанием р = 0,20; 0,24; 0,28 при Fig. 7. Normalized correlation functions of strain fields in an elastoplastic problem at different straining points for structures with spherical inclusions of different volume fractions р = 0,20; 0,24; 0,28 at функции полей напряжений и деформаций. Проведено исследование в зависимости от объемной доли включений. Корреляционные функции полей деформаций для структур с содержанием пор от 15 до 28 % можно аппроксимировать одним выражением с одинаковыми константами (рис. 7) Выводы Были рассмотрены частные случаи представительных объемов многокомпонентных композитов со случайным расположением эллипсоидальных включений. Разработана и реализована методика вычисления значений корреляционных функций полей деформирования в упругом и упругопластическом случае при различных условиях нагружения и физико-механических характеристиках компонентов композитов. Для получения значений условных и безусловных корреляционных функций полей напряжений и деформаций в представительных объемах и их компонентах использовано решение стохастической краевой задачи в полном корреляционном приближении. В качестве входных параметров для вычислений выступают константы, характеризующие свойства компонентов представительного объема, и входящие в функции Грина в подынтегральных выражениях, а также геометрические свойства компонентов и представительного объема, определяемые структурными моментными функциями. Построены корреляционные функции напряжений и деформаций для разреженных структур со сферическими и эллипсоидальными полыми включениями. Проведено исследование и выбор аппроксимирующих зависимостей для полученных моментных функций. Полученные численные результаты могут использоватьсяв стохастических моделях для предсказания момента начала разрушения микроструктурных компонентов материала.

About the authors

M A Tashkinov

Perm National Research Polytechnic University

N V Mikhailova

Perm National Research Polytechnic University

References

Statistics

Views

Abstract - 186

PDF (Russian) - 67