Full Text

Введение Использование многоуровневой концепции описания структурно-неоднородного материала [1-7] позволяет моделировать его механическое поведение на разных масштабных уровнях и уже на начальных стадиях деформирования выявить очаги концентрации напряжений. Это дает возможность глубже понять закономерности разрушения, являющегося в соответствии с общепринятыми представлениями [8-9] сложным многостадийным процессом возникновения и развития внутренних микродефектов. Известно, что стадия образования микроскопических трещин (стадия скрытого или рассеянного разрушения), феноменологически описываемая теориями механики поврежденности, косвенно характеризует накопление повреждений в микрообъемах сплошного материала [9-11]. Величину повреждений связывают с пластическими деформациями и оценивают по величине предельной накопленной деформации до разрушения, зависящей, в свою очередь, от истории изменения коэффициента жесткости напряженного состояния и показателя вида напряженного состояния Лоде-Надаи в процессе деформации. Таким образом, описание разрушения невозможно без получения данных об особенностях напряженно-деформированного состояния (НДС) материала и истории его изменения при развитой пластической деформации. Особенно актуально это для структурно-неоднородных сред, обладающих сложной нерегулярной внутренней структурой, оказывающей определяющее влияние на неоднородность деформационных и прочностных характеристик материала. Применительно к исследованию данных материалов существует два основных подхода решения задач, базирующихся на использовании стохастических методов [12-14]. Первый подход заключается в определении эффективных или вероятностных оценок параметров НДС с использованием моментных функций структуры материала. Информацию о структуре получают, используя образцы материала или модель случайной структуры и имитационное моделирование. Тогда со статистической точки зрения задача состоит в определении характеристик стохастических полей напряжений и деформаций в элементах структуры по ее известным статистическим свойствам [14]. В основу второго подхода положены численные методы (например, метод Монте-Карло), базирующиеся на формировании генеральной выборки случайных реализаций микроструктуры исследуемого материала, с тем чтобы вероятностные характеристики НДС для этой выборки совпадали с аналогичными характеристиками для некоторого представительного объема материала (макрообъема) [13]. При этом, как показала практика, для усреднения результатов достаточно иметь не менее 10 реализаций микроструктуры [5]. В настоящей работе на примере одного случайно выбранного фрагмента микроструктуры методологически показан способ численного исследования характера эволюции НДС структурно-неоднородного материала на макро- и микромасштабном уровнях одноосного растяжения и сжатия с учетом особенностей строения и реологии его компонентов. 1. Материал и методика исследований В качестве модельного материала использовали дисперсно-упрочненный металломатричный композит (ММК), матрицей которого является алюминиевый сплав А8, а наполнителем - частицы карбида кремния SiC двух типоразмеров 1-5 и 15-20 мкм, имеющие преимущественно форму неправильных призм. Микроструктура исследованного композиционного материала представлена на рис. 1 (снимок на растровом электронном Рис. 1. Микроструктура исследованного ММК Fig. 1. Microstructure of the studied metal matrix composite микроскопе). Экспериментально обнаружено прочное адгезионное взаимодействие между матрицей и частицами наполнителя [15]. Детально особенности микроструктуры и свойств заготовок из данного ММК описаны в работах [15-16]. Для проведения численных расчетов по исследованию механического поведения при нагружении ММК потребовалось создание его вычислительной модели, связанное непосредственно с разбиением моделируемого материала на отдельные области простой формы. При создании данной модели использовали опыт авторов, полученный ранее при разработке модели деформации сложнолегированной латуни [5]. Вычислительная модель ММК была создана на основе двухуровневого структурно-феноменологического подхода, связывающего решение задач на макро- и микромасштабном уровнях [1, 12]. Согласно этому подходу на первом (макромасштабном) уровне материал рассматривали как однородную, изотропную, изотропно упрочняющуюся среду с усредненными по объему свойствами, определяемыми экспериментально. На втором масштабном уровне - микроуровне, материал представили как микронеоднородную среду, состоящую из связных областей, моделирующих структурные составляющие материала. Считали, что в пределах структурных элементов модели остаются справедливыми все феноменологические уравнения и соотношения механики континуума. Основываясь на результатах исследований [15] при проектировании объема ММК на микроскопическом уровне полагали, что между частицами карбида кремния SiC и алюминиевой матрицей существует прочная адгезионная связь. Благодаря этому допущению геометрическая модель микрообъема ММК представляет собой трехмерный кусочно-однородный континуум, имитирующий алюминиевую матрицу, в которой располагаются частицы SiC, конфигурация которых близка к реальной форме, полученной по результатам стереологического анализа [15-16]. Рассматриваемый в данной статье микрообъем представлял собой куб с длиной ребра 30 мкм. Для учета влияния окружающих слоев материала вокруг микрообъема размещали буферный слой с усредненными механическими свойствами (макросвойствами), полученными по результатам испытаний стандартных образцов. Толщина буферного слоя принята равной линейному размеру микрообъема. Составленная таким образом вычислительная модель ММК на микроуровне представляет собой 3D-композицию из структурно-неоднородного микрообъема ММК в окружении буферного слоя и может рассматриваться на макроскопическом уровне как некий макрообъем композита с усредненными механическими свойствами. Моделирование нагружения вычислительной модели позволяет детально исследовать и описывать эволюцию НДС случайно выбранной микроструктуры композита. В этом случае граничные условия задаются микрообъему ММК в результате решения задачи на макроуровне, а выполняющие роль буфера слои материала позволяют более точно передать напряженно-деформированное состояние на микроуровень. Поскольку качество построения конечно-элементной дискретизации во многом определяет корректность результатов, получаемых при численном моделировании, то были разработаны приёмы построения трёхмерных сеток по геометрически нерегулярным структурам и создан программный комплекс [17], позволяющий получать 3D-модели объемов неоднородных материалов с учётом их сложной внутренней структуры в формате, подходящем для использования с конечно-элементным комплексом ANSYS. С помощью разработанного программного комплекса спроектировали 3D-модель микрообъема ММК, имеющую случайную нерегулярную структуру, представляющую собой включения частиц SiC в виде призм с фасками, окруженных алюминиевой матрицей (рис. 2). Рис. 2. Трехмерная вычислительная модель ММК Fig. 2. Three dimensional numerical model of metal matrix composite Реологические свойства алюминиевого сплава А8 и материала ММК задавали полученными экспериментально[2] кривыми деформационного упрочнения, представляющими зависимость напряжения течения от степени деформации (рис. 3). Рис. 3. Кривые деформационного упрочнения исследуемого ММК (1) и алюминиевого сплава А8 (2) Fig. 3. Curves of strain strengthening of the studied metal matrix composite (1) and aluminuim alloy А8 (2) Материал матрицы рассматривали как изотропную упругопластическую и пластически несжимаемую среду с изотропным деформационным упрочнением, особенности которого были установлены экспериментально при осадке образцов со скоростью деформации 1 с-1 при температуре 300 оС. Свойства буферного слоя, определенные при осадке образцов ММК на макроуровне в аналогичных условиях нагружения, соответствовали изотропной упруговязкопластической среде. Материал частиц карбида кремния полагали изотропным, подчиняющимся закону линейной упругости. Задавали модуль Юнга Е = 70 ГПа и коэффициент Пуассона = 0,34 для алюминиевого сплава А8 [18] и Е = 380 ГПа, = 0,19 для карбида кремния [19]. Для буферного слоя Е = 194 ГПа, = 0,26. Эти значения определили по правилу смесей [20] в зависимости от объемного содержания алюминиевого сплава А8 и карбида кремния в материале ММК. Численные расчеты механического нагружения ММК производили на его трехмерной геометрической вычислительной модели в квазистатичной постановке в программном комплексе ANSYS на вычислителе кластерного типа URAN ИММ УрО РАН. При разбиении Рис. 4. К заданию граничных условий при моделировании растяжения ( = 20 мкм, ), сжатия ( -17 мкм, ) Fig. 4. Setting boundary conditions under strain simulation ( = 20 microns, ), compression ( -17 microns, ) сетки использовали 268242 тетраэдральных конечных элемента типа SOLID 187, имеющих 10 степеней свободы. Из них 186223 элемента моделировали матрицу композита. Граничные условия задавали в перемещениях (где i = x, y, z) по соответствующим j-м граням буферного слоя (рис. 4). Такие граничные условия дают возможность рассматривать геометрическую модель ММК на макроуровне как 1/8 часть целого, находящегося в состоянии одноосного растяжения или сжатия вдоль вертикальной оси у за счет закрепления отдельных граней модели ММК, как плоскостей симметрии от перемещений в соответствующих направлениях. Предельную величину перемещений задавали из условия обеспечения величины степени макродеформации геометрической вычислительной модели, равной 0,2. Величину определяли как , где , - начальная и конечная высота геометрической вычислительной модели ММК. В результате моделирования были получены сведения об изменении значений компонент тензоров напряжений и приращений деформаций в узлах конечно-элементной сетки вычислительной модели ММК. По расчетным данным определяли коэффициент жесткости напряженного состояния и показатель вида напряженного состояния Лоде-Надаи , которые являются основными безразмерными параметрами, характеризующими напряженное состояние: , где - среднее нормальное (гидростатическое) напряжение; Т - интенсивность касательных напряжений, равная пределу текучести на сдвиг в пластической области. , где , , - главные напряжения. 2. Результаты исследований и их обсуждение Результаты численного моделирования растяжения и сжатия ММК на макроуровне как квазиоднородного материала показали формирование однородных полей характеристик НДС на каждом расчетном шаге. Установлено, что эффективные значения показателей напряженного состояния соответственно: = 0,577 и = -1 при растяжении; = -0,577 и = 1 при сжатии и остаются постоянными по мере нагружения. Расчетные значения показателей соответствуют известным литературным данным [10], что подтверждает адекватность результатов моделирования. Моделирование на макромасштабном уровне позволяет в среднем оценить уровень значений основных параметров НДС в том или ином виде нагружения, но не является информативным с точки зрения получения данных по выявлению локальных очагов концентрации напряжений и деформаций, наличие которых обусловлено в первую очередь микроструктурными особенностями материала. Соответственно отсутствует возможность отследить динамику изменения безразмерных параметров напряженного состояния с учетом неоднородного характера НДС. В этом смысле определяющее значение имеют результаты решения задачи на микроуровне. В отличие от однородных полей макронапряжений и макродеформаций, распределения параметров НДС на микроуровне характеризуются существенной неоднородностью на каждом расчетном шаге в рассмотренных видах нагружения. При этом результаты численных расчетов подтверждаются соответствующими экспериментальными данными и дают возможность их обоснования. Как показали испытания на растяжение образцов, процесс разрушения ММК инициируется и в дальнейшем определяется появлением и развитием трещин в алюминиевой матрице, в то время как наполнитель SiC демонстрирует высокую прочность, практически не деформируясь [15-16]. Численное моделирование позволило установить, что жесткие недеформируемые поверхности карбидов кремния являются причиной возникновения участков локальной пластической деформации и создают условия, схожие с теми, что возникают при растяжении тонкого пластичного слоя, расположенного между двумя жесткими поверхностями. По мере увеличения степени деформации объемная доля этих участков возрастает, что приводит к усилению неоднородности деформированного состояния на микроуровне (рис. 5). C другой стороны, твердые частицы наполнителя SiC, присутствующие в мягкой матрице, формируют в непосредственной близости от себя зоны концентрации растягивающих напряжений. В этом удалось убедиться путем определения коэффициента жесткости напряженного состояния k в узлах конечно-элементной сетки матрицы ММК. По расчетным данным были построены поля распределений k внутри матрицы в зависимости от степени деформации в каждом виде нагружения. Для удобства анализа данные поля были визуализированы средствами ANSYS с помощью специально разработанной программы. Рис. 5. Распределение накопленной деформации в матрице ММК по центральному сечению xy микрообъема ММК при сжатии со степенью деформации: (а) e = 0,04; (б) e = 0,2 Fig. 5. Distribution of accumulated strain in the metal matrix composite along the central section xy of the microvolume in the metal matrix composite under compression with strain rate: (а) e = 0,04; (b) e = 0,2 Известно, что при > 0 напряженное состояние характеризуется преобладанием растягивающих напряжений; при < 0 - сжимающих напряжений. При этом высокий уровень растягивающих напряжений способствует интенсивному пластическому разрыхлению и ускоряет процесс разрушения [9-10, 21]. На рис. 6, а показано распределение коэффициента жесткости напряженного состояния k в матрице по центральному сечению xy микрообъема ММК при растяжении со степенью деформации e = 0,2. Максимальные значения k, обусловленные близостью включений, имеют место непосредственно в характерных участках образования микротрещин, наблюдаемых в экспериментах. Подобный эффект концентрации растягивающих напряжений выявлен также при моделировании Рис. 6. Распределение коэффициента жесткости напряженного состояния в матрице ММК по центральному сечению xy микрообъема ММК при растяжении (а), сжатии (б) со степенью деформации e = 0,2 Fig. 6. Distribution of rigidity coefficient of strain state in the metal matrix composite along the central section xy of the microvolume in the metal matrix composite under stretching (а), compression (b) with strain rate e = 0,2 деформации сжатия (рис. 6, б). Положительные значения , в случае сжатия, также формируются непосредственно вблизи частиц карбида кремния и указывают на доминирующий характер растягивающих напряжений в этих областях. Причем по мере увеличения деформации объемная доля участков растягивающих напряжений увеличивается. Так при степени макродеформации e = 0,04 положительные значения определяются у 9 % узлов конечно-элементной сетки матрицы, при деформации e = 0,08 таких узлов уже 12 %, а при e = 0,2 положительные значения имеют более 25 % узлов матрицы ММК. Статистический анализ расчетных данных позволил установить, что для рассмотренных схем нагружения на каждом шаге деформации распределение значений в матрице подчинено нормальному закону. В качестве примера на рис. 7 показана плотность распределения значений при сжатии. Рис. 7. Плотность распределения коэффициента жесткости напряженного состояния матрицы ММК при сжатии в зависимости от степени деформации e: (1) e = 0,04; (2) e = 0,08; (3) e = 0,2 Fig. 7. Frequency distribution of rigidity coefficient of strain state of metal matrix composite under compression depending on strain rate e: (1) e = 0,04; (2) e = 0,08; (3) e = 0,2 При сжатии по мере увеличения деформации среднее значение по узлам конечно-элементной сетки матрицы возрастает. При растяжении наблюдается обратная закономерность: по мере нагружения среднее значение в матрице убывает. На величину среднего значения существенное влияние оказывает неоднородность напряженного состояния на микроуровне с наличием участков концентрации напряжений, а также изменение текущей картины напряженного состояния за счет повышения абсолютного уровня напряжений и их перераспределения в процессе нагружения. Статистическая обработка результатов расчета показателя вида напряженного состояния Лоде-Надаи показала, что по аналогии со значениями коэффициента жесткости напряженного состояния распределение значений в матрице ММК подчинено нормальному закону как при растяжении, так и при сжатии. При этом среднее значение показателя Лоде-Надаи практически не зависит от степени деформации и составляет при растяжении -0,19, при сжатии 0,17, в то время как амплитудные значения этого параметра на отдельных участках матрицы достигают величины ±0,99 во всех видах нагружения. Полученные распределения и , а также закономерности их изменения, в первую очередь объясняются тем, что микрообъем ММК не является представительным объемом исследуемого материала, а являет собой случайную реализацию его микроструктуры. Поэтому на микроуровне реализуется специфическое НДС, характерное лишь для этого конкретного фрагмента ММК. Статистическая выборка подобных фрагментов, и численное исследование нагружения каждого из них с использованием разработанной вычислительной модели, позволяет обобщить результаты моделирования и вывести общие закономерности эволюции НДС ММК на микроуровне. Выводы 1. На основе двухуровневого структурно-феноменологического подхода разработана 3D вычислительная модель деформации неоднородного материала, учитывающая особенности строения и реологии его составляющих компонентов. Выполнена численная реализация разработанной модели на примере моделирования механического поведения случайно выбранной микроструктуры металломатричного композита Al/SiC для одноосного растяжения и сжатия. 2. Установлено, что наиболее неблагоприятное с точки зрения разрушения напряженное состояние в мягкой матрице возникает в непосредственной близости от частиц карбида кремния, выступающих, с одной стороны, в роли концентраторов растягивающих напряжений. С другой стороны, жесткие поверхности карбидов кремния являются причиной возникновения участков локальной пластической деформации и создают условия, схожие с теми, что возникают при растяжении тонкого пластического слоя, расположенного между двумя жесткими поверхностями. Области концентрации растягивающих напряжений и локальной пластической деформации материала матрицы наиболее неблагоприятны с точки зрения разрушения независимо от вида нагружения. Именно там вероятнее всего зарождение первых трещин, которые по мере роста общей деформации будут распространяться по объему матрицы ММК. 3. Построены поля распределения основных безразмерных параметров напряженного состояния: коэффициента жесткости напряженного состояния и показателя вида напряженного состояния Лоде-Надаи в узлах конечно-элементной сетки матрицы ММК, позволяющие дать качественную и количественную оценку напряженного состояния выбранного фрагмента микроструктуры ММК для рассмотренных схем нагружения на каждом шаге деформации. Статистическая выборка и обработка результатов численного моделирования нагружения подобных фрагментов с использованием разработанной вычислительной модели позволяет установить закономерности эволюции НДС исследуемого композита на микроуровне. Авторы выражают благодарность главному программисту Института математики и механики УрО РАН А.С. Игумнову за помощь в проведении расчетов в программном комплексе ANSYS.

About the authors

S V Smirnov

Institute of Engineering Science

Email: svs@imach.uran.ru

A V Konovalov

Institute of Engineering Science

Email: avk@imach.uran.ru

M V Myasnikova

Institute of Engineering Science

Email: marina@imach.uran.ru

Yu V Khalevitsky

Institute of Engineering Science

Email: me@dijkstra.ru

A S Smirnov

Institute of Engineering Science

Email: smirnov@imach.uran.ru

References

Statistics

Views

Abstract - 216

PDF (Russian) - 100