CC BY
90
УДК 539.5
ВЛИЯНИЕ ПОРОВОЙ СТРУКТУРЫ ХРУПКОЙ КЕРАМИКИ НА РАЗРУШЕНИЕ ПРИ ДИНАМИЧЕСКОМ НАГРУЖЕНИИ
В.А. Скрипняк, Е.Г. Скрипняк, А.А. Козулин, Е.Г. Пасько, В.В. Скрипняк, М.В. Коробенков
Томский государственный университет E-mail: [email protected]
Исследовано влияние поровой структуры на разрушение хрупких керамических материалов при динамическом сжатии. Показано, что поровые кластеры оказывают существенное влияние на размеры фрагментов при разрушении пористых керамических материалов. Сделан вывод о том, что при ударно-волновом нагружении наличие кластеров в поровой структуре приводит к образованию продольных трещин.
Ключевые слова:
Конструкционная керамика, пористость, динамическое нагружение, локализация деформации, микроповреждения, компьютерное моделирование. Key words:
Structural ceramics, porosity, dynamic loading, localization of deformation, microdamages, computer simulation.
Пористые керамические материалы на основе А1203, ZrO2, SiC, В4С широко применяют в инженерной практике как теплозащитные и конструкционные материалы. Предсказание их механического поведения при динамическом нагружении представляет существенную сложность. Это связано с влиянием параметров формы и распределения пор на мезоскопическом уровне на закономерности деформации и разрушения материалов [1-3]. Исследования влияния поровой структуры на механическое поведение пористой керамики интенсивно ведутся как экспериментальными [1], так и теоретическими методами [2-6]. В данной работе для изучения локализации деформации при динамическом нагружении с амплитудами до 10 ГПа предложена модификация метода [4].
Распределение пор в керамике определяется технологическими процессами ее получения. Современные методы рентгеновской томографии позволяют определять пространственное распределение макро- и мезопор. Для получения данных о кластерах пор используются результаты статистического анализа распределения пор на поверхности образцов. Данные о поровой структуре учитываются при создании моделей для компьютерного моделирования процессов деформации и разрушения керамических материалов под действием нагрузок. Под поровой структурой принято понимать совокупность параметров, характеризующих распределение, форму и размеры пор в объеме тела (коэффициент формы пор - отношение объема сферы, вписанной в пору, к объему сферы, описанной около поры, эффективный размер пор йр среднее расстояние между близлежащими порами Ьр, эффективный радиус кластера пор Хс/, среднее расстояние между кластерами пор Ьс).
Модель элементарного объема керамики с размерами пор, существенно превышающими размеры зерна, можно рассматривать в виде двухфазной среды, состоящей из конденсированной фазы и полостей пор, заполненных газом.
Для заданного относительного объема пор а и размера пор йг определяется число пор, которые должны содержаться в элементарном объеме. Удобно создавать модель, располагая поры в вершинах, ребрах или гранях многогранников Вороного, с помощью которых может быть создан элементарный объем. При равномерном распределении пор они располагаются в вершинах правильных тетраэдров. В этом случае обеспечивается одинаковое расстояние между близлежащими порами. Количество пор N в элементарном объеме модельной среды с характерным размером Ь соответствует их относительному объему
Vp
а= — =
4 N
3 nL (rp )
где г{=й{/2; ¥р - суммарный объем пор; У=Ьг -объем элементарного объема пористой среды.
В монодисперсном приближении является постоянной величиной, а относительный объем пор определяется соотношением:
а = (ф)пЫ (гр/Ь)3.
Расстояние между порами / и у определяется по формуле [7]
ц. = {[X(° - X(15 ][X('5 - Xи5 ]}1/2,
где Xй - радиус-вектор центра /-ой поры.
Минимальное расстояние между ближайшими порами Ь{ в элементарном объеме при отсутствии перекрытия их объемов определяется как 2гр. При фиксированном количестве пор N в элементарном объеме, в рамках предположения о равномерном пространственном распределении пор, минимальное расстояние между порами определяется как
LP *
бУ2
N - 3
1/3
b.
При наличии поровых кластеров, минимальное расстояние между порами в кластере обозначено а расстояние между кластерами (рис. 1, б). Кластеры пор удобно вводить внутри многогранника. Для формирования распределения кластеров пор с расстоянием между ними, кратном характерному размеру многогранников, может быть использована дискретная модель расчетного объема. Метод может быть использован для построения модельной структуры пористости для задачи в двумерной постановке.
На рис. 1 показаны два типа модельных элементарных объемов (с одинаковыми пористостью а=0,11 и коэффициентом формы 0,74, размерами поперечного сечения 300x172 мкм). На рис. 1, а, структура характеризуется равномерным распределением пор со средними размерами 25,2 мкм, минимальным расстоянием между центрами близлежащих пор ¿~100 мкм. На рис. 1, б, структура содержит кластеры из 6 пор с размерами 10,3 мкм. Расстояние между порами в кластере - 12,6 мкм, расстояние между кластерами - 100 мкм.
Элементарный объем керамического материала представляет многосвязную пространственную область, в которой расположены полости пор, окруженные конденсированной фазой (рис. 1). В случае, когда размеры полостей многократно превышают характерные размеры зерен или кристаллитов, для описания механического поведения конденсированной фазы может быть использован континуальный подход механики повреждаемых сред [6]. Введем систему отсчета, расположив оси декартовой системы координат ОХ1, ОХ2 в плоскости модели. Кинематика конденсированной фазы керамики описывается. компонентами тенз.оров скоростей деформации ¿¡¡и изгиба-кручения
. 1
¿У= 2
ды. ды:
дХ: дХ:
. 1
а = — 2
ды: дыу
дх, дх.
V : '
где щ - компоненты вектора эффективной массовой скорости в конденсированной фазе.
В лагранжевой системе отсчета уравнения сохранения массы, импульса и энергии имеют вид
дО
дхУ
- = р-
Сы: Л
С р ды -С1 дх..
СЕ .
где Оц - компоненты тензора напряжений в конденсированной фазе среды; р - массовая плотность; Е - удельная внутренняя энергия на единицу массы.
Определяющее уравнение в повреждаемой среде представим в виде
О: = (1 - Б)[-Р8у + ],
где Р - давление; ^ - девиатор тензора напряжения; ОО - символ Кронекера; D=Д¿1'/¿f - параметр поврежденности среды; Дер=\'0ерЖ, е1=Б1(Р'+Т*)Д; ¿р - интенсивность скорости пластической деформации; Бь Бг - постоянные материала, Т*=(Т-Т)/(ТВ-Т), Тт - температура плавления конденсированной фазы; Т=293 К; Р *=Р/РНЕ1, РНЕ1 - давление, соответствующее пределу упругости Гюгонио.
Механическое поведение конденсированной фазы описано моделью Джонсона-Холмквиста [8]. Давление Р в пределах интенсивностей ударных волн до 10 ГПа может быть рассчитано с помощью уравнения состояния
• Р=Х1в+Х2в2+Х3в'+Гр0Е, при сжатии 6>0;
• Р=К19, при растяжении 0<0,
где К1, К2, К3 - постоянные материала, 0=(р/р0)-1, Г - коэффициент Грюнайзена.
Девиатор напряжения вычисляется в рамках модели Друккера-Прагера с пластическим потенциалом Джонсона-Холмквиста [8]
а
г
= 2/
■с 1■ О
¿У--¿кк О у
3 -
где / - модуль сдвига конденсированной фазы.
Производная по времени Яумана тензора
а
вычисляется в соответствии с определением
dSч . . .
= S ц - Sл - ^ Щ ■
ч ч ч ' ч
\ д8
да ч
2 SчSч
1/2
а, = А(Р + Т )т (1 + С 1п ¿),
где А, С, т - постоянные материала, е=ёи/1 [с-1] -нормированная интенсивность тензора скорости деформации; X - коэффициент.
Использован критерий локального разрушения конденсированной фазы, принятый в рамках подхода механики сред с внутренними повреждениями
(Л=1).
Начальные условия соответствуют состоянию с нулевыми полями механических параметров. На-гружение элементарного объема ударным импульсом определяется граничными условиями. На верхней поверхности модельного объема задается условие нагружения импульсом давления. На нижней поверхности модельного объема задаются трансляционные граничные условия. На боковых границах задано отсутствия смещений в направлении, перпендикулярном движению волны. Такие граничные условия соответствуют макроскопической деформации среды в плоской ударной волне. Полости пор, как правило, заполнены газом. Вследствие значительного различия сжимаемости газов и конденсированных фаз керамических соединений
влиянием сжимаемого газа на сопротивление коллапсу пор при динамическом нагружении можно пренебречь. Тогда на поверхности пор граничные условия соответствуют таковым для свободных поверхностей.
Для конденсированной фазы на примере А1203-керамики приняты следующие значения коэффициентов [3]: р0=3,96-103 кг/м3; К1=251 ГПа; К2=0; К3=0; Г=1; ц=161 ГПа; А==0,001; Я2=1; Рт=1,46 ГПа; А=0,93 ГПа; С=0; т=0,6.
Задача решается в квазидвумерной постановке. Задача для пористого модельного объема среды сформулирована в 3^-постановке, в которой предполагается отсутствие смещений в направлении оси х3.
Решение осуществлялось конечно-разностным методом с использованием лагранжевого решателя программного комплекса ANSYS АитоиШ [10]. Устойчивость и сходимость численных решений достигались выбором пространственного и временного шага дискретизации.
Разрушению керамических материалов предшествует локализация деформации и развитие микроповреждений в конденсированной фазе [3-5]. Формирование зон локализации при интенсивном динамическом нагружении обусловлено возникающей неоднородностью поля напряжений. На рис. 2 показаны расчетные поля сдвиговых напряжений в модельных объемах керамики в области ударного перехода, а на рис. 3 - соответствующие поля ги-
Рис. 2. Распределение интенсивности напряжений в элементарных объемах структурированных оксид алюминиевой керамики при нагружении плоской ударной волной с амплитудой 5 ГПа
а б
Рис. 3. Поле гидродинамического давления за фронтом ударной волны в структурированной пористой оксид алюминиевой керамике
дродинамических давлении в момент времени 14 нс. Результаты расчетов показывают, что поры являются мезоскопическими концентраторами сдвиговых напряжении. Форма изолированных пор определяет характер распределения напряжении в объёме материала и коэффициент концентрации напряжении. С увеличением относительного расстояния от поры Ь{/й{>5 коэффициент концентрации сдвиговых напряжении уменьшается.
При сжатии элементарного объема в ударнои волне неупругая деформация локализуется вокруг пор и в полосах локализации. Полосы локализации имеют три характерных ориентации - вдоль, поперек и под углами, близкими к 45° к направлению распространения фронта волны. Формирование полос, имеющих поперечную ориентацию, обусловлено изменением формы пор и сопровождается резким ростом поврежденности в конденсирован-нои фазе. Под деиствием градиента гидростатического давления во фронте протекает трансформация формы пор из сферическои в эллипсоидную и далее в дискообразную. Изменение формы пор вызывает резкое увеличение концентрации напряжении в плоскости, ортогональнои направлению распространения волны. Полученные данные согласуются с выводами, представленными в [9].
Изменение формы сопровождается уменьшением объема пор за счет перемещения конденсиро-ваннои фазы. Из-за ограниченнои пластичности керамических материалов коллапс пор в керамических материалах происходит за счет сдвигов и разворотов микрофрагментов. Фрагменты разделены полосами локализованнои деформации. Сдвиги фрагментов обуславливают кратковременное локальное уменьшение гидродинамических давлении вокруг пор, рис. 3. Локальная релаксация давлении вызывает появление градиентов массовои скорости в направлении, ортогональном направлению распространению волны. Появление распределения массовых скоростеи в зоне ударного перехода приводит к формированию полос, ориентированных параллельно направлению распространения волны.
Локальное разрушение конденсированнои фазы керамики за фронтом ударнои волны показано на рис. 4.
Формирование зон разрушения (макроскопических трещин) сопровождается релаксации сдвиговых напряжении. Зоны разрушения разделяют элементарныи объем на фрагменты. Размеры фрагментов, сформировавшихся при сжатии в ударнои волне, зависят от скорости деформации и параметров поровои структуры. Средние размеры фрагментов, образовавшихся при сжатии пористои керамики, уменьшаются при увеличении амплитуды нагрузки и скорости деформации в объеме материла. Увеличение размеров пор при сохранении инте-гральнои пористости приводит к образованию распределения фрагментов по размерам. Из приведенных на рис. 4 результатов следует, что размеры фрагментов, образующихся при интенсивном динамическом нагружении, могут различаться на несколько порядков. Для оценки среднего размера фрагментов I в [3] предложена формула
I = К1С / рае0
,2/3
где К1С - трещиностоикость для конденсированнои фазы, а - скорость изменения относительного объема пор, скорость деформации е0.
Трещиностоикость высокоплотнои керамики К1С [МПа-м1/2] при гомологических температурах Т/Тт<0,1 может варьироваться в пределах 4...5 для А1203, 6...10 для ZrO2, 3...4,5 для SiC, 3...3,6 для В4С в зависимости от размеров зерна и содержания химических примесеи. Изменение средних размеров фрагментов для керамических материалов одного класса после нагружения в сходных условиях не превышает 30 %. Более широкое распределение размеров фрагментов, которое наблюдают в экспериментах, может быть обусловлено наличием кластеров пор в поровои структуре материала. В процессе ударного сжатия керамики с уменьшением размеров пор снижается а и увеличиваются размеры фрагментов и уменьшается суммарная площадь разрушения. В результате поверхностная энергия разрушения уменьшается.
Выводы
1. В условиях ударно-волнового нагружения деформации конденсированнои фазы пористои керамики локализуется на мезоскопическом уровне в полосах, имеющих три характерные
ориентации - вдоль и поперек и под углами, близкими к 45° к направлению распространения фронта волны.
2. Наличие кластеров пор в материале приводит к образованию узких и прямолинейных полос. Формирование полос сопровождается релаксацией сдвиговых напряжений и ростом локальной поврежденности конденсированной фазы в полосах, что приводит к неоднородности и не-
СПИСОК ЛИТЕРАТУРЫ
1. Louro H.L., Meyers M.A. Effect of stress state and microstructural parameters on impact damage of alumina-based ceramics // J. Mat. Sci. - 1989. - V. 24. - № 3. - P. 2516-2532.
2. Maiti S., Rangaswamy K., Geubelle P.H. Mesoscale analys of dynamic fragmentation of ceramics under tension// Acta Materialia. -2005. - V. 53. - № 3. - P. 823-834.
3. Скрипняк В.А., Скрипняк Е.Г., Жукова Т.В. Повреждаемость керамических покрытий и конструкционной керамики при интенсивном импульсном нагружении // Химическая физика - 2002. - Т. 21. - № 9. - С. 76-82.
4. Скрипняк В.А., Скрипняк Е.Г., Каракулов В.В. Моделирование ударно-волнового нагружения наноструктурной керамики и керамических композитов // Экстремальные состояния вещества. Детонация. Ударные волны. Матер. Межд. конф. IX Харитоновские тематические научные чтения-2007. - Са-ров: РФЯЦ ВНИИЭФ, 2007. - С. 374-398.
5. Псахье С.Г., Смолин А.Ю., Стефанов Ю.П. и др. Моделирование поведения сложных сред на основе комбинированного дискретно-континуальногог подхода // Физическая мезомеха-ника. - 2003. - Т. 6. - № 6. - С. 11-21.
стационарности поля напряжений на мезоско-пическом уровне. 3. Макроскопические прочностные характеристики исследованных материалов, в которых микропоры образуют кластеры, ниже. Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (проект 08-08-12055) и Минобрнауки РФ (программа АВЦП «Развитие научного потенциала высшей школы», проекты 2.1.1/6909, 2.1.1/5993).
6. Curran D.R., Seaman L. Simplified models offracture and fragmentation // Ed. L. Davison, D.E. Grady, M. Shahinpoor. High-pressure shock compression of solids II - dynamic fracture and fragmentation. - Berlin: Springer, 1996. - P. 340-365.
7. Макмиллан Н. Идеальная прочность твердых тел // Атомистика разрушения. Сб. статей. Пер. с англ. / Сост. А.Ю. Ишлин-ский. - М: Мир. 1987. - 248 с.
8. Johnson G.R., Holmquist T.J. A computational constitutive model for brittle materials subjected to large strains, high strain rates and high pressure // Shock-Wave and High Strain Rate Phenomena in Materials / Ed. M.A. Meyers, L.E. Murr, K.P. Staudhammer, M. Dekker. - N.Y.: AIP Press, 1992. - 1075 p.
9. Paskaramoorthy R., Meguid S.A. On the dynamic behaviour of porous materials // Int. J. of Solids and Structures. - 2000. - V. 37. -№ 16. - P. 2341-2358.
10. ANSYS AUTODYN Explicit software for nonlinear dynamics. Theory manual [Электронный ресурс]. - режим доступа: http://www.autodyn.org. - 29.06.2009.
Поступила 29.06.2009 г.
