Фазовые превращения в материалах с включениями

Филиппов, Роман Александрович

Фазовые превращения в материалах с включениями тема автореферата и диссертации по механике, 01.02.04 ВАК РФ

Филиппов, Роман Александрович АВТОР

кандидата физико-математических наук УЧЕНАЯ СТЕПЕНЬ

Санкт-Петербург МЕСТО ЗАЩИТЫ

2013 ГОД ЗАЩИТЫ

01.02.04 КОД ВАК РФ

Диссертация по механике на тему «Фазовые превращения в материалах с включениями»

Автореферат диссертации на тему "Фазовые превращения в материалах с включениями"

На правах рукописи

005538830

Филиппов Роман Александрович

ФАЗОВЫЕ ПРЕВРАЩЕНИЯ В МАТЕРИАЛАХ С ВКЛЮЧЕНИЯМИ

Специальность 01.02.04 — Механика деформируемого твердого тела

Автореферат 21 НОЯ 2013

диссертации на соискание учёной степени кандидата физико-математических наук

Санкт-Петербург — 2013

005538830

Работа выполнена в Федеральном государственном бюджетном учреждении науки Институт проблем машиноведения Российской Академии паук (ИПМаш РАН)

Научный руководитель: доктор физико-математических наук,

Фрейдин Александр Борисович

Официальные оппоненты: доктор физико-математических наук,

профессор Волков А. Б.

кандидат физико-математических наук, доцент Устинов К.Б.

Ведущая организация: Федеральное государственное бюджетное

учреждение науки

Институт прикладной механики РАН, г.Москва

Защита состоится 12 декабря 2013 года в 14ю часов на заседании диссертационного совета Д 002.075.01 при Институте проблем машиноведения РАН по адресу: 199178, Санкт-Петербург, Большой пр. В.О., 61.

С диссертацией можно ознакомиться в ОНТИ Института проблем машиноведения РАН.

Автореферат разослан 12 ноября 2013 года.

Ученый секретарь

диссертационного совета Д 002.075.01 доктор технических наук, профессор

В.В. Дубаренко

Общая характеристика работы

Актуальность темы.

Работа посвящена развитию моделей механики, описывающих фазовые превращения в материалах с включениями при термомеханических воздействиях.

Исследования фазовых превращений в деформируемых телах представляют актуальное направление современной механики материалов. Применительно к дисперсным композитным материалам следствием фазовых (структурных) превращений вокруг частиц является формирование переходных слоев измененного материала матрицы и, как следствие, изменение Примерами являются сплавы на основе титана, которые испытывают мартенсит-ные превращения и в которых при определенных термомеханических условиях обработки выделяются частицы, не претерпевающие фазовое превращение и способные инициировать вокруг себя (в матрице) фазовые превращения. Фазовые превращения могут также происходить внутри частиц. Примером являются фазовые превращения частиц диоксида циркония в керамических композитах, приводящие к эффектам трансформационного упрочнения.

В диссертационной работе рассмотрен двухкомпонентный композитный материал одна из компонент которого (материал включений или матрицы) может находится в двух фазовых состояниях, различающихся модулями упругости и собственной деформацией превращения. Исследованы две задачи, в которых устанавливается связь между внешними термомеханическими воздействиями и фазовым превращениями.

Первая задача - задача описания фазовых превращений, идущих в матрице вокруг включений. В результате формирования областей новой фазы возникают переходные слои, увеличивающие эффективный объем включений и, как следствие, возрастает влияние объемной доли включений на упругие свойства композитного материала.

Влияние переходных слоев, окружающих включения, на эффективные упругие свойства композитов обсуждалось с использованием различных приближений в работах Р.В. Гольдштейна, М. Качанова, И. Севостьянова, К.Б. Устинова, Ф. Бонднолн, С. Бутлеба и др. В работах С.А. Лурье развита градиентная модель межфазпого слоя для описания свойств материала в окрестности включения.

В отмеченных выше работах влияние переходных слоев рассматривалось безотносительно механизмов их формирования. В данной работе исследуется возникновение переходных слоев как равновесных областей новой фазы. Задача о двухфазной конфигурации исследуется в постановке, развитой в

работах М.А. Гринфельда, В.А. Еремеева, Л.М. Зубова, Н.Ф. Морозова, В.Г. Осмоловского, А.Л. Ройтбурда, А.Б. Фрейдина, Р. Абейаратне, М.Е. Гёртина, Р.Д. Джеймса, Дж. Ноулса и др., согласно которым равновесная межфазная граница должна удовлетворять дополнительному термодинамическому условию, что, в свою очередь приводит к зависимости размеров области новой фазы от напряжений. В результате объемная доля получающихся составных включений становится управляемым параметром, зависящим от термомеханических воздействий на технологической стадии изготовления композита.

Вначале исследовано развитие сферического слоя новой фазы вокруг сферического изолированого включения. Затем рассмотрены фазовые превращения вокруг взаимодействующих сферических включений в композитном материале. Построены зависимости равновесного радиуса составных включений от внешней деформации. Исследована устойчивость межфазных границ, а также энергетические изменения и перераспределение напряжений вследствие развития областей новой фазы. Исследована возможность увеличения эффективных модулей упругости композита в результате формирования устойчивых равновесных областей новой фазы. Отметим, что исследование устойчивости межфазных границ остается одной из центральных проблем механики материалов, претерпевающих фазовое превращение.

Во второй задаче рассмотрены фазовые превращения взаимодействующих включений в дисперсном композитном материале. Разработанная модель, позволяет определять диапазон критических размеров включений диоксида циркония (Zr02) в эффекте трансформационного упрочнения керамик. Эффект трансформационного упрочнения наблюдается в керамиках с диспергированными в них метастабильными включениями ZrOi- Включения ZrO2 могут претерпевать фазовое превращение мартеиситного типа, сопровождающееся положительной собственной деформацией превращения. Если фазовое превращение включений инициируется полем напряжений, создаваемым трещиной, то перераспределение напряжений, вызванное увеличением размеров частиц, приводит к замедлению или блокированию роста трещины, в чем и заключается эффект трансформационого упрочнения.

с включениями, находящимися в разных фазовых состояниях, с учетом изменений химической энергии, энергии деформаций и поверхностной энергии. Однако без учета дополнительных факторов в критерии превращения превалирующая роль изменения химической энергии приводит к незначительному влиянию напряжений на критический размер включений и, как следствие, очень узкому диапазону размеров включений, обеспечивающих трансформационное упрочнение, что ставит иод сомнение адекватность такой постановки задачи. В настоящей работе в модель вводится энергетический барьер, преодоление которого служит условием фазового превращения, и предлагаются методики определения величины энергетического барьера на основе экспериментальных данных, а также учитывается взаимодействие включений. В результате определяется диапазон размеров включений Zr02, обеспечивающих эффективное трансформационное упрочнение.

Основной целью диссертационной работы является разработка и реализация моделей, описывающих влияние внешних тсрмомеханических воздействий на фазовые превращения в материалах с включениями.

Задачами работы являются:

1. Разработка и исследование модели, описывающей развитие переходных слоев в дисперсных композитных материалах как областей новой фазы. Исследование условий зарождения и устойчивого роста равновесных сферических границ областей новой фазы и их влияния на локальные поля напряжений и эффективные модули упругости.

2. Разработка и исследование модели, описывающей фазовые превращения включений в дисперсных композитах, позволяющей оценивать напряжения, индуцированные фазовым превращением частиц и определить диапазон размеров включений диоксида циркония в эффекте трансформационного упрочнения керамик.

Научную новизну диссертации представляют следующие положения, выносимые на защиту:.

1. Впервые поставлена и исследована задача о формировании переходных слоев в материале с включениями как термодинамически равновесных областей новой фазы.

- Для случая сферической симметрии определен равновесный радиус межфазной границы в зависимости от внешней деформации и модулей упругости фаз композита.

- Развит новый метод исследования устойчивости межфазной границы,

основанный на соотнесении деформаций па межфазпой границе с гра-

ницами зоны фазовых переходов. Определены диапазоны радиусов межфазной границы и внешних деформаций, в которых сохраняется устойчивость.

- Впервые построены зависимости эффективных упругих модулей композита с матрицей, претерпевающей фазовое превращение, от управляющей внешней деформации. Для изотропного композита показана возможность увеличения эффективного объемного модуля упругости при условии уменьшения эффективного модуля сдвига.

- Показано, что формирование переходного слоя новой фазы приводит к сбросу энергии деформации в матрице и, как следствие, к уменьшению концентрации напряжений. Определены ограничения на параметры материала, при выполнении которых возникновение переходного слоя приводит к сбросу энергии деформаций во всех компонентах композита.

- На примере двухфазных центрально-симметричных полей деформаций показано принципиальное различие формирования равновесных устойчивых областей новой фазы в однородном материале и вокруг включений: возникновение областей новой фазы в однородном материале возможно, только если модуль сдвига в результате фазового превращения возрастает, а возникновение областей новой фазы вокруг включений возможно, только если модуль сдвига уменьшается.

2. На основе новой модели трансформационного упрочнения, учитывающей энергетический барьер, который должен быть преодолен на пути превращения, найдены диапазоны размеров включений Zr02, обеспечивающих эффективное трансформационное упрочнение в керамиках с учетом упругого взаимодействия включений и термоусадочных напряжений. Предложены методики определения величины барьера на основе экспериментальных данных. Теоретически показано, что величины критических радиусов уменьшаются с ростом объемной концентрации включений ¿?гОг-

Научно-практическая значимость работы заключается в создании и исследовании моделей механики деформируемого тела, учитывающих влияние термомеханических воздействий на формирование областей новой фазы в материалах с включениями, в том числе в дисперсных композитах. Результаты исследования могут быть использованы при разработке новых материалов, которые за счет управляемого фазового перехода обладают повышенными жесткостью и трещиностойкостью.

Достоверность полученных результатов обеспечивается строгой постановкой математических задач, применением математически обоснован-

пых методов решения, использованием в численных процедурах надежных алгоритмов и программ, совпадением численных результатов с полученными для частных случаев аналитическими результатами.

Апробация работы. Основные результаты исследований, представленные в диссертационной работе, докладывались и обсуждались на International Workshop on NanoBiotcchnologies (Санкт-Петербург, 2006), международной школе-конференцни "Advanced Problems in Mechanics" (Санкт-Петербург, 2005, 2007, 2009, 2011, 2013), Всероссийском съезде по теоретической и прикладной механике (Нижний Новгород, 2011), Europian Solid Mechanics Conference (Грац, Австрия, 2012), BaltMatTrib 2012 (Tallinn, Estonia, 2012), XVIII Зимней школе по механике сплошных сред (Пермь, 2013).

Полностью результаты диссертации обсуждались на семинарах ИПМаш РАН, лаборатории неклассических моделей композиционных материалов и конструкций ИПРИМ РАН и кафедры теоретической механики СПбГПУ.

Публикации. По теме диссертации автором опубликовано 5 статей, из них 4 - в изданиях, рекомендованных ВАК России, и тезисы 10 конференций.

Личное участие автора в работах, написанных в соавторстве, состоит в получении аналитических и численных решений поставленных задач и их исследовании.

Структура и объём работы. Диссертация состоит из введения, четырех глав, заключения и списка литературы. Работа изложена на 90 страницах машинописного текста, содержит 22 рисунка, список использованных источников из 104 наименований.

Исследования поддерживались грантами РФФИ (07-01-00525-а, 10-01-00670-а, 11-01-16070-моб_з_рос, 12-01-09285-моб_з, 13-08-00553-а, 13-01-00687-а), программой ОЭММПУ РАН №13 (рук. акад. РАН И.Г.Горячева) и программой фундаментальных исследования госакадемий РФ №23 (рук. акад. РАН И.Г.Горячева и акад. РАН Н.Ф.Морозов), гранта Президента Российской Федерации для государственной поддержки ведущих научных школ Российской Федерации НШ-3776.2010.1 (рук. акад. РАН Н.Ф.Морозов), Министерством образования и науки РФ (договор 14.В25.31.0017), программой безвалютного обмена между РАН и АН Эстонии.

Содержание работы

Во введении обоснована актуальность темы диссертации, дана общая характеристика работы, приводится обзор публикаций по теме диссертации, указаны основные цели работы, кратко изложена структура диссертации, охарактеризована ее научная новизна, научная и практическая значимость, сформулированы основные положения, выносимые на защиту.

В первой главе исследовано развитие переходного слоя повои фазы вокруг изолированного включения. Для случая сферической симметрии определены локальные поля напряжений и перераспределение плотности энергии деформаций в результате фазового превращения.

В параграфе 1.1 приведена математическая постановка задачи определения равновесной межфазной границы при бездиффузионном фазовом превращении мартенситного типа, сопровождающегося собственной деформацией превращения и изменением модулей упругости1.

Задача описания равновесной двухфазной деформации упругого тела состоит в нахождении межфазной границы Г и соответствующего поля перемещений и(х), достаточно гладкого при х $ Г, непрерывного на Г и удовлетворяющего условиям равновесия и граничным условиям:

х^Г: V • сг = 0, в = const (1)

х 6 Г : [и] = 0, Н-п = 0 (2)

[/] - И : [в] = 0, (<т)=Л(о-+ + <х_) (3)

где е и <т - тензоры деформации и напряжения, в - температура, / = /(е, в) - плотность свободной энергии Гельмгольца, п - единичный вектор нормали к Г, индексами "—" и "+" обозначены величины, относящиеся к исходной и новой фазе матрицы соответственно, квадратными скобками обозначен скачок величины на Г, [ • ] = (■)+ — (•)-• Фазы матрицы "+" и "—" и различаются тензорами модулей упругости С± и тензорами деформаций в ненапряженном состоянии. Если ер_ = 0, то ер = е+ - собственная деформация превращения. Массовые силы, термоупругие напряжения и поверхностная энергия межфазной границы не учитываются.

Условия (1), (2) - обычные условия равновесия составного тела. Условие термодинамического равновесия (3) связано с дополнительной степенью свободы, порождаемой неизвестной межфазной границей2. Подчеркнем, что

1 Кубланов Л.Б., Фрейдин А.Б. Зародыши твердой фазы в деформируемом материале //

ПММ. 1988. Т. 52. С. 493-501.

2 Гринфельд М.А. Об условиях термодинамического равновесия фаз нелинейно-упругого материала // Докл. АН СССР. 1980. Т. 251. № 4. С. 824-827.

условия (1)-(3) являются только необходимыми условиями минимума энергии двухфазного тела. Не всякая двухфазная деформация, удовлетворяющая этим условиям, устойчива. Поэтому найденное равновесное решение нужно дополнительно исследовать на устойчивость.

Плотность свободной энергии / материала матрицы (рис. 1) и включений задана квадратичными зависимостями:

/(М) = тт{/-(М),/+(М)}

(4)

(е -4)

С,

(е-4),

Ме) = -е:С(:е (5)

где - плотности свободной энергии фаз в ненапряженном состоянии, верхние и нижние знаки и индексы ± соответствуют друг другу, ер = 0, ер = £р -деформация превращения, С; - тензор модулей упругости материала включений. Определяющие соотношения фаз матрицы и материала включения принимают вид:

<г±(е) = С±:(е-4), <г = С4 : £ (6)

Отметим, что одно из ненапряженных состояний может быть только гипотетическим: если /+(ер) > /~{ер), то деформация е — ер принадлежит области определения фазы "—", а не фазы "+". Далее будем говорить, что имеем дело со случаем а, если /+(ер) < /~(ер) (рис. 1а), и со случаем 6, если /+(ер) > /-(ер) (рис. 16). Реализация того или иного случая зависит от параметра 7, который в свою очередь зависит от температуры.

Условие непрерывности усилия (2) с учетом условия непрерывности перемещений позволяет записать термодинамическое условие (3) в виде зависимости между нормалью к межфазной границе и деформацией на одной из сторон границы:

/ :

(Ъ)

Рис. 1: Зависимости /(е, в):

(а) /+(е") < /-(£"),

(б) Г(е') > Г{е")

Х(Ч±, 7*>п) = 7* + ^ Ч±: (Сх 1 ± Кт(п)) : = О

(7)

вг

-1

где Ч± = : е± - С+ : ер,С! = С+ - С_, Ъ = 7 ' Кт(п) = {п ® Ст(п) ® п}8, в^п) = (п • Ст • п)-1.

В параграфе 1.2 рассмотрен сферический слой новой фазы радиуса Яг, формирующийся вокруг включения радиуса Л, в неограниченной матрице в однородном ноле деформации е0 = ^оЕ (рис. 2). Материалы включений

и фаз матрицы - изотропные, кг,±, Ц{,± - модули объемного сжатия и сдвига включения и фаз матрицы. Принимается, что собственная деформация превращения является шаровой £+ = (?9Р/3)Е, где Е - единичный тензор.

В сферических координатах перемещения и = и(г)ег и деформации еГ, £д = определяются формулами

О Л , 21)

и = Аг + —, ег = — = А--

г аг гл

и . В

= = - = А+- (8)

где зависимости констант А — {А{, А+, А-} и Б = {Д, 1) \ , £>_} от внешней деформации Со и относительного обратного радиуса межфазной границы р = — € [0,1] определены включение в матрице претерпеваю- /1р

щей фазовое превращение 113 Усл0В11И

£г |г-юс= ^о/з = го, и |г=0< оо; г = Я;, г = Дг : м = о, [аг] = О

Рис. 2: Изолированное сферическое

Нетрудно убедиться в том, что й+<Г_

Л = + = = А = О,

Б. _ п _ „ _ 3(Ы+г-ЯЬ) , _

= I ¿у = К+-4/Л г)1

2± = 5+ = Г>%е0 + 1У°, Щ-

г = (1+ - 12//х г, г = р3, = 3к± + 4/х+, = 3к± + Ац. (1± = 3к{+ 4(1±, к\ = к+ — кА^ = - ц+- £ =

После подстановки (8) при г = Яг и п = еГ в термодинамическое условие равновесия (7) получено уравнение для определения зависимости г* = рЦ^о) ~ куба обратного равновесного радиуса от внешней деформации

<1_ к\ Ъ = 7 - к-к+(,вр)2/2к\ 10

На рис. 3 приведена зависимость о) = Т^/Щ- при к+ > к-, ц+ < /1_ и х)р < 0, где - равновесный радиус межфазной границы. Величина р, монотонно меняется от 1 до 0 при увеличении сжимающих деформаций по модулю е0 от когда происходит зарождение фазы "+" на поверхности включения, н до когда матрица полностью переходит в фазовое состояние "+". Случай ¡х+ > Ц- не обсуждается, так как будет показано что при таком выборе параметров сферические межфазные границы неустойчивы.

В параграфе 1.3 исследованы про-р» странственпые перераспределения плотности ^—г- 1 энергии деформаций, окружных напряжений | ад и интенсивности касательных напряжений

/ ; т = \ад — аг\/2 в теле с включением вследствие

/| | формирования равновесного переходного слоя

Г, \ ; новой фазы.

I \ ; | Параметры фаз композита (Табл. 1) вы-

—'т-1—1-■—'—- бирались, исходя из следующих соображений:

£' К X 0 £0 а„ ^ п г- с -

„ „ , . ур < 0 обеспечивает уменьшение свободной

гнс. о: зависимость р,(£о) при

> к и < 0 энергии тела в результате фазового перехода

при сжатии £о < 0; /л+ < ц- - обеспечивает устойчивость сферических границ но отношению к возмущению радиуса (см.Параграф 2.1); к+ > к- - выполнение этого условия позволяет обеспечить рост эффективного объемного модуля сжатия композита при условии падения эффективного модуля сдвига < /х_; кроме того при проведении расчетов материал влючений был выбран более жестким, чем материал матрицы. Величины модулей упругости и параметра у в Табл. 1 даны в одних тех же единицах измерения. Случаи а и Ь отличаются только величиной 7.

Таблица 1: Параметры материалов фаз матрицы и включения для случая а и Ъ

т к+ /'+ к. 7

а 500 170 140 6 20 9 -0.08 0.02

Ь 500 170 140 6 20 9 -0.08 0.07

В работе показано, что в обоих случаях а и Ъ в результате возникновения равновесного переходного слоя новой фазы происходит локальный сброс плотности энергии деформаций как в матрице впе составного включения, так и в исходном включении. В случае Ъ плотность энергии деформаций уменьшается также и в области переходного слоя.

Случай Ь также является более предпочтительным, чем случай о, так как в случае Ъ окружные напряжения ад во всех областях сжимающее, а

т падает и в переходном слое, и в матрице. В случае а в переходном слое напряжение сте растягивающее, а т после фазового превращения возрастает. Величина интервала превращения — £s| больше в случае Ь.

Во второй главе исследована устойчивость межфазной сферической границы. Условия (1)-(3) являются условиями экстремума функционала энергии по отношению к возмущениям поля перемещений и межфазной границы. Поэтому найденные равновесные решения могут оказаться неустойчивыми.

В параграфе 2.1 исследована устойчивость межфазной границы но отношению к возмущению радиуса границы.

Формула Эшелби для энергии взаимодействия включения с внешним полем в случае сферической границы области новой фазы с радиусом Яг и с учетом (9) принимает вид:

^ = (7 + - Зв0) + f (З^о - t)A+) + - г)

Зависимость E{z) имеет экстремум при г £ [0,1]. Если < ß-, то это минимум. Минимальность E(z) означает, что найденные равновесные межфазные границы являются устойчивыми по отношению к возмущениям радиуса границы. Значения Е в точках минимума - отрицательные, что свидетельствует о том, что возникновение равновесной области новой фазы уменьшает свободную энергию тела. Если > то экстремум соответствует максимуму E(z), что означает неустойчивость равновесной межфазной границы. Поэтому случай > далее не рассматривается.

Устойчивость по отношению к возмущениям радиуса границы является только необходимым условием устойчивости границы, поэтому в параграфе 2.2 и параграфе 2.3 исследована устойчивость найденных равновесных центрально-симметричных решений относительно осесимметричных возмущений и произвольных возмущений соответственно.

Следуя кинетическому подходу к анализу устойчивости1 перемещения точек тела и положение межфазной границы в возмущенном состоянии в сферической системе координат задаются формулами:

u= u*(r) + w(r,0), гг= 7Гг + г?(0), w (r,ff)=wr{r,e)eT + we{r,ff)ee

В приближении линейной термодинамики необратимых процессов вместо условия равновесия (7) для возмущенного состояния имеем:

vTn = -*x{e±, n) =-xx(q±,n), х > О (11)

■Еремеев В.А., Фрейдин А.Б., Шарипова Л.Л. Об устойчивости равновесия двухфазных

упругих тел // Прикл. математика и механика. 2007. Т. 71. № 1. C.G6-92.

где - скорость движения границы Г в направлении нормали, х ~~ термодинамическая сила , х - кинетический коэффициент.

Решение линеаризованной возмущенной краевой задачи относительно возмущения границы искалось в виде полиномов Лежаидра т](в) = (пРп(со$ в). Тогда кинетическое уравнение преобразуется к виду уравнения эволюции п- й моды возмущения амплитуды

L О

-0.1

-0.2

Кх

= п = 2,3,... (12)

е* -0.023 -0.021

Рис. 4: Зависимость L

£ £о = Ьа(ео)

где Ьп{ео) - кинетический коэффициент, зависящий от модулей упругости фаз, собственной деформации превращения и внешней деформации. Корни уравнения Ьп{е0) = 0 - точки бифуркации возмущенной задачи. Отрицательные значения Ьп соответствуют экспоненциальному затуханию со временем амплитуды возмущений межфазиой границы. Положительные значения Ьп соответствуют росту возмущений н означают неустойчивость.

На рис. 4 приведена зависимость Ь — Ь^ от внешней деформации £0 для случая Ь. В работе показано, что в обоих случаях а н Ъ в интервале превращения при некоторой внешней деформации 6 (г/, г") выполняется равенство Ь(е^) = 0 (точка К на рис. 4), то есть происходит потеря устойчивости межфазиой границы по отношению к осссимметричным возмущениям.

Осесимметричные возмущения не исчерпывают все возможные возмущения. В связи с этим в работе предложен метод оценки устойчивости по отношению к произвольным, в том числе негладким возмущениям межфазной границы, основанный на соотнесении деформаций иа межфазной границе с границами зоны фазовых переходов.

Термодинамическое условие (7) может быть удовлетворено не при всех деформациях. Это привело к понятию зон фазовых переходов (ЗФП)1 - областей, образованных в пространстве деформаций всеми деформациями, которые могут существовать в данном теле на равновесных межфазных границах. ЗФП определяется исключительно параметрами материала и является своего рода паспортом материала, претерпевающего фазовые превращения.

'Фрейдин А.Б., Чискис A.M. Зоны фазовых переходов в нелинейно-упругих изотропных материалах // Изв. АН. МТТ. 1994. №4. С. 91-109 (Ч. 1)

В недавних работах Ю. Грабовского и Л. Трускиновского1 и М.А. Ап-тимонова, A.B. Черкаева и А.Б. Фрейдина2 показано, что необходимым условием устойчивости межфазной границы но отношению к произвольным, в том числе негладким возмущениям границы, как и условием соответствия двухфазной деформации глобальному минимуму энергии, является принадлежность деформаций на межфазной границе внешним границам ЗФП. Это позволяет использовать ЗФП для отбраковки неустойчивых решений, не конкретизируя класс возмущений.

На рис. 5 показаны сечення ЗФП плоскостью осесимметричных деформаций £\ = £т, £ — £g = в области пересечения внешних и внутренних границ ЗФП для случая Ь. ЗФП состоит из двух подзон "+" (слева от пунктирной линии) и "—" (справа от пунктирной линии), образованных деформациями, которые могут существовать на межфазных границах со стороны фаз "+" и "—" соответственно. Внешние и внутренние границы подзон пересекаются в точках X и X', соответствующих состоянию всестороннего

£i

-0.03

-0.04

N А

X' V ^

\ ^

Л. \

V \

-0.04 -0.03 е Рис. 5: Сечения ЗФП

растяжения/сжатия фаз

и "+" соответственно.

В рассматриваемом случае деформации на межфазных границах принадлежат внешним границам ЗФП при выполнении условия (£g — £~) |г > 0. Точка X на рис. 3 и 4 соответствует внешней деформации, при которой па межфазной границе = £~, то есть при которой происходит потеря устойчивости равновесной сферической межфазной границе по отношению к произвольным возмущениям.

Пары сосуществующих деформаций на равновесных сферических межфазных границах при различных внешних полях £oi и £02 обозначены на рис. 5 квадратиками и крестиками соответственно. Величины £qi и £02 выбраны так, что деформации на межфазной границе, соответствующие £оь находятся на внешних границах ЗФП, то есть устойчивы по отношению к произвольным возмущениям, а в случае £q = £02 одна из деформаций попадает на внутреннюю границу, что означает неустойчивость. Но устойчивость по отношению к осесимметричпым возмущениям при £q = £q2 все еще сохраня-

'Grabovsky Y., Truskinovsky L. Roughening instability of broken extremals // Arch. Rat.

Mech. Anal. 2011. Vol. 200. № 1. P. 183-202.

2М.А.Антимонов, А.В.Черкаев, А.Б.Фрейдин. Оптимальные микроструктуры и точная нижняя граница энергии упругих композитов из двух изотропных фаз // Научно-технические ведомости СПбГПУ. Физико-математические науки. 2010. №3. С. 112-122

ется: на рис. 3 и 4 деформация е{)2 занимает положение между деформациями, соответствующими точкам К н X.

Анализ устойчивости показал, что при к+ > деформации на сферической межфазной границе, начиная с некоторого радиуса границы, перестают соответствовать внешним границам ЗФП. Это означает ограниченность устойчивого развития области новой фазы: начиная с некоторого радиуса межфазная граница теряет устойчивость. Отметим, что для случая b интервал внешних деформаций, при которых равновесная сферическая межфазная граница остается устойчивой, шире по сравнению со случаем а. Поскольку случаи а и b отличаются только величиной параметра 7, зависящим от температуры, то интервал устойчивых радиусов межфазной границы также зависит от температуры.

В заключении этой главы па классе центрально-симметричных двухфазных деформаций, показано принципиальное различие формирования устойчивых равновесных межфазных границ в однородном материале и вокруг включения: возникновение устойчивых равновесных областей новой фазы в однородном материале возможно, только если модуль сдвига в результате фазового превращения возрастает, а для фазового превращения вокруг включения - когда модуль сдвига убывает.

В третьей главе исследована возможность формирования устойчивых мезфазных границ в композитном материале, состоящем из матрицы, способной претерпевать фазовый переход, и изотропно распределенных взаимодействующих сферических включений. Исследована возможность увеличения эффективных модулей упругости композита в результате формирования устойчивых переходных слоев новой фазы.

Композит подвергнут всестороннему сжатию, которому соответствует средняя но объему деформация £q = £оЕ. Начиная с некоторой деформации вокруг включений формируются области новой фазы.

Следуя основной гипотезе метода эффективного поля1 считаем, что каждое включение в композите ведет себя как изолированное включение в однородной среде с модулями упругости С_ под действием эффективного внешнего поля е», которое складывается из поля Ео и полей, наведенных окружающими включениями. Для случая изотропного распределения включений в пространстве эффективное поле определятся формулой:

et = (l-pA:P0rl:(e0-pA:Pl) (13)

'Канаун С.К., Левин В.М. Метод эффективного поля в механике композитных материалов. Петрозаводск. Издательство Петрозаводского ун-та. 1993. 538 с.

где

1 — аз_„„ 5 — 2.., 1__ч 3к- +

А = —--ЕЕ Н——-(I — -ЕЕ), ае_ - р

3 ' Зк- + 4ц-

Ро = Р1(к±,и Д±,;,Ро/р+)ЕЕ + Рг(fci.ii ц±*,Ро/р+) (^1 - ^ЕЕ^ (14)

Р£ = ^Р0рЕ, Р0р = ^ + А:П = А- - ** (15)

I - единичный тензор четвертого ранга, /ун = 1/2+ ро -

объемная концентрация включений, - объемная концентрация областей,

занятых новой фазой, р - объемная концентрация составных включений, р0+р+=р.

Из (13) следует, что при гидростатическом нагруженпи эффективное поле является шаровым:

е, = е.Е, £. - , Ро = (к1А%р+ - кйА!ро) (16)

где Л|, АЕ+, Лf и Л+ определяются соотношениями (9) при замене в них г на,ро/р+. Поле перемещений и деформаций на межфазной границе типичного включения определяется формулами для изолированного включения (8), находящегося в эффективном поле (16).

Отмечается, что в рамках приближения эффективного поля возможно рассмотрение двух типов потери устойчивости: относительно возмущений концентрации и относительно возмущений "эффективной" межфазной границы при фиксированном эффективном поле. Во втором случае проверка устойчивости композита с переходными слоями новой фазы сводится к проверке устойчивости межфазной границы в случае изолированного включения, находящегося в фиксированном эффективном внешнем поле, что позволяет ограничиться соотнесением деформаций на эффективной межфазной границе с ЗФП.

Средняя плотность энергии взаимодействия ё областей новой фазы с внешним полем имеет вид

ё = (27 + 4(1?" - 34) + 9/С1Л+4 - 3*Л+ + 12щВ+О0.) л _ по = ^д*" р* -

+ яу ' <г-Щ" 0 ,1- + Зр(1ки

где - эффективное поле в случае композита с однофазной матрицей (когда Р+ = 0).

0.0001

0.0002

Зависимость средней плотности энергии взаимодействия от объемной концентрации составных включений при заданном значении средней деформации для случая < Ц- приведена на рис. 6. Минимум на кривой соответствует равновесной концентрации р = р* и означает устойчивость структуры по отношению к возмущениям концентрации р. Отрицательность энергии взаимодействия означает, что возникновение равновесных переходных слоев энергетически выгодно.

В работе получено выражение для тензора эффективных модулей С», который связывает средние напряжения ¿г и средине деформации ё = £о

Ро 0.1 0.2 Р*

Рис. 6: Зависимость ё = ё(р,)

& = С* : ео + <тр

где С*

К ЕЕ + 2ц* (I - -ЕЕ

к* = к- +

pPi

- тензор эффективных модулей упругости РР'2

1-pPl

1 — зе_

2(1

5 — 2ее_^

рР2—7Z-)

1.12 1.08 1.04 1

- ✓ * h*

/ Е„

Ео ц*

р,- ~ 15/i-

Индуцированные собственной деформацией напряжения

«г* = р (I - Ро : (I - РА : Ро)"1 : А) : Pg

0.02 0.06 0.1 Р*ц* Рис. 7: к'/кд, Д*//4о>

На рис. 7 приведены графики зависимостей относительного эффективного модуля сжатия к,/kg, относительного модуля сдвига и относительного модуля Юнга Е*/Eq от концентрации слоистых включений для случая < fi-. Здесь fcg, Pq, Eq - эффективные модули исходного композита с включениями в однофазной матрице в фазе "—".

Пунктирные части линий на рис. 7 для эффективных модулей и на рис. 6 для энергии взаимодействия соответствуют концентрациям, которые являются равновесными и существуют при таких внешних деформациях, что на межфазной границе эффективного составного включения локальные деформации со стороны фазы "—" уходят с внешних границ ЗФП, что может рассматриваться как потеря устойчивости. Отмстим, что это происходит при

большем радиусе эффективного слоистого включения в композите, чем в случае изолированного включения, что позволяет говорить о стабилизирующей роли взаимодействия включений. Также по сравнению с изолированным включением происходит возрастание интервала деформаций, в котором происходит превращение.

Важно отметить, что фазовое превращение не может привести к росту всех эффективных модулей упругости. В рассматриваемом случае сферической симметрии модуль сдвига новой фазы должен быть меньше модуля сдвига исходной фазы. Поэтому возможно увеличение только эффективного объемного модуля упругости.

В четвертой главе исследована модель, позволяющая определять диапазон критических размеров метастабильных включений Zr02 в эффекте трансформационного упрочнения керамик с учетом впервые вводимого энергетического барьера. Разработаны методики определения величины энергетического барьера.

Рассмотрен композитный материал со взаимодействующими сферическими включениями ¿ГгОг радиуса И, которые могут находиться в аустенит-ном А или мартенентном М состояниях в неограниченной матрице. Переход Л —М сопровождается положительной собственной деформацией превращения ер = (г9р/3)Е. Задача решена при условии идеального контакта на границе матрицы и включений. Исследован случай, когда все включения находятся в одном и том же фазовом состоянии в однородном поле внешних напряжений и температур.

Объемная плотность свободной энергии Гельмгольца и определяющие соотношения термоунругих фаз включения и матрицы заданы формулами

Де,Т) = тт{/л(е,Т),Ме,Т)}

л,м

/А(е, Т) = ¡£(Т) + - аААТЕ) : Сл : (е - аААТЕ)

/м(е, Т) = (Т) + \{е - е? - амАТЕ) : Сл/ : (е - ер - ал/ДТЕ) (17) <тм = Сл/:(е - £р - ал/ДТЕ),аА,о = Сло:(е - «лоЛТЕ)

где Сд/;д,о и ад/,д,о - тензоры модулей упругости и коэффициенты линейного расширения материалов включений, находящихся в разных фазовых состояниях и матрицы, соответственно, индексом 0 обозначены величины, относящиеся к матрице, Т) - объемные плотности химических энергий включений в разных фазовых состояниях, АТ - перепад температур.

Принимается, что при остывании переход включений из аустенитного

в мартеиситное состояние происходит, когда разность объемных плотностей

энергий Гиббса композита с включениями в аустепитном О а и мартснситном состояниях превышает некоторый энергетический барьер Z. Для случая сферической симметрии условие перехода А М принимает вид

4 4

<?л - <3л/ = -ттЛ,3(у(Т) + ф(ер, АТ, аех,р)) + 4тгпЯ2£ >г = -тгпВ?г (18)

о о

где 7(Т) = /\}{Т) — /д®(Г), £ = Си — Сл ~ разность плотностей поверхностной энергии включений в разных фазовых состояниях, ф(ер, АТ, стех,р) - объемная плотность изменения потенциальной энергии деформаций, порождаемая собственной деформацией превращения ер, внутренними термоусадочными напряжениями, возникающими при изменении температуры композита на АТ, и внешним полем напряжений асх, р = |7гпД3 - объемная концентрация включений, тг - количество включений в единице объема композита, 2 - объемная плотность энергии барьера. Природа энергетического барьера может быть связана с перестройкой кристаллической решетки при фазовом превращении и (или) с энергозатратами на образование зародышей новой фазы, величина которых зависит от плотности центров зародышеобразования и стабилизирующих добавок.

Из (18) следует, что выражения для оценки радиусов включений, обеспечивающие пребывание включений в метастабпльном состоянии при остывании керамики до эксплутацнонных температур (11т) и переход в стабильное состояние во внешнем поле напряжений (Я/), а также относительный диапазон размеров включений определяются формулами

Ит ^

d =

~/(Т) + ф(£р,АТ,аех = 0,р) + г

Rf = (Г) + ф{ер, AT, (Тех = + (19)

RT-Rf _ Ф(ер, AT, аех = af,p) - ф(ер,АТ,аех = 0,р)

Rt 1{Т) + ip(eP, AT, <тех = af,p) + z

где ip(ep, AT,aex = 0,p) и ip(ep, aex = <т/,р) - объемные плотности изменения потенциальной энергии деформаций при фазовом превращении включений в отсутствие внешних напряжений и во внешнем поле напряжений <7/, эквивалентном полю, наводимому трещиной.

В приближении эффективного поля с использованием формулы Эшел-би для энергии взаимодействия включения с внешним полем получено выражения для гр

ф(ер, АТ,асх,р) = -

-{ам - <т\,): (ер + (ам - а0) ATE) + (стА - <r\): (аА - aQ) ATE) (20)

где Вд = С^1 — Сд \ В\7 = Сд/ — Сд1 - разности тензоров податливости, а* и сг\[ - эффективные поля внешних напряжений для матрицы с включениями в разных фазовых состояниях. Для случая изотропных фаз и шарового внешнего поля, сгех = <тС1Е, эффективные внешние напряжения и напряжения внутри включений определены соотношениями

, _ к0(аех(ЗкА + 4/х0) + 12рр0кААТ(аА - а0))

Л ~ кй(ЗкА + 4/л0) + 4рц0кА{кА - к0)

к0(аех(3кы + 4/х0) + 4РЦаки^ + ЗЛТ(аЛ/ - а0))) сгы =-^^-—-:-г-^т-—-^

f м =

ка(3км + 4/х0) + 4рц0км(км - к0) = кА{осх{Зко + 4ц0) - 12ц0к0АТ(аА - а0)(1 - p))g Л ко(ЗкА + 4цо) + 4рр0кА(кА - к0)

Ыаех(3к0 + 4/i0) - 4/x0fco(l - Р)(#р + ЗЛТ(ал/ - ар)))

ко(Зкм + 4¿¿0) + 4рц0км{км - к0)

Таблица 2: Термоупругие параметры материалов фаз включений и матрицы

a,K~l E, ГПа v к, ГПа ti, ГПа

A, Zr02 12 * 1(T6 207 0,25 138 83

M, Zr02 12 *.10"6 207 0,25 138 83

0, A12OÍ 8 *10"6 398 0,23 245 161

0, WC 5,3* 10"6 680 0,2 377 283

Собственная деформация превращения равна др = 0,05, критическая внешняя нагрузка выбрана равной 400 МПа - пределу прочности Л120з на растяжение, АТ = lOOO/'t'. Параметры при которых производился расчет, приведены в Таблице 2. Зависимость 7 (Т) аппроксимирована линейной функцией температуры 7 (Т) = -q(l-T/Tt), где q = 251,74 МПа-тепловой параметр превращения, Т* = \322K - температура, при которой jf¡{T„) = fAl{T»).

Объемная плотность энергетического барьера z оценивается двумя способами. Первый способ основан на использовании результатов экспериментального определения относительной ширины диапазона критических радиусов включений ZrO2 в Л/гОз-матрнце. Для объемной концентрации включений р = 0,1 экспериментально определенные величины критических радиусов составляют R¡ = 50нм, Rt = ЗООнм1, откуда следует, что d = 5/6. Из третьей формулы в (19), с использованием (20), (21) и данных из Табл.2, получена оценка величины плотности энергетического барьера 2 = 147 МПа.

^.Н. De Aza, J. Chevalier, G. Fantozzi, M. Schehl, R. Torrecillas. Crack growth resistance of alumina, zirconia and zirconia toughened alumina ceramics for joint prostheses // Biomaterials.-2002,- Vol. 23,- P. 937-945.

Подчеркнем, что при таком определении параметра г не требуется знать величину трудноопределимого в прямых экспериментах изменения плотности поверхностной энергии. Вычисленная при г = 147 МПа по первой формуле (19) величина изменения плотности поверхностной энергии £ = 0.31 Дж/м2, что хорошо согласуется с приводимой в литературе величиной изменения поверхностной энергии при мартенситном превращении Zr02■

Величина г также оценивалась с использованием экспериментальных данных1, устанавливающих связь между величинами критического радиуса Лу кристаллов свободного от напряжений (■ф — 0) порошка Zr02 в зависимости от содержания стабилизирующего оксида У2О3. Принималось, что добавка У2О3 не влияет на механические свойства Zr02 и изменение поверхностной энергии £ при мартенситном превращении. Из (19) следует формула для вычисления барьера в зависимости от содержания 1*>Оз

z> 120

МПа

100

Y2O2, моль%

Рис. 8: Зависимость ¿(УгОз) для свободного от напряжений порошка ZrOi.

г(У2Оз) = _7(Т.(У20з))

Rt(Y203)

(22)

На рис. 8 представлена зависимость 2(У20з) для свободного от напряжений порошка Zr02- Величина г растет с увеличением содержания оксида иттрия, причем зависимость г(У20з) быстро возрастает при увеличении содержания У20з с 0 до 1 моль-% и затем слабо меняется: расхождение между величинами z(l — то1У2Оз) и z{ 1,5 — то1У2Оз) составляет 4%. При содержании Y2O3, равном 1 моль%, наблюдается наилучшее соответствие величин барьера, полученных первым (z = 147 МПа) и вторым (z = 151 МПа) способами.

На рис. 9 представлены зависимости относительной ширины диапазона критических радиусов d мстастабильпых включений ZrO2 от z в AI2O3- и WC-матрицах соответственно при объемной концентрации включений ZrO2 10% при Т=293 К. Вертикальной штрихпунктирной линией на Рис. 9 отмечена абсолютная величина изменения суммы плотности химический энергии и энергии деформаций |7(Г) + ф\ , которая является асимптотой для величин критических радиусов. При г —>■ |7(Т) + ф\ Rtj —»■ 00.

'Aran Suresli, Merrilea J. Mayo, Wallace D. Porter, Claudia J. Rawn. Crystallite and Grain-Size-Dependent Phase Transformations in Yttria-Doped Zirconia // Journal of the American Ceramic Society.-2003.- Vol. 86.- No. 2. - P. 3C0-362

При 2 = 0 величина ¿{г) не превышает 0,1 для обеих матриц. Получить экспериментально столь узкий диапазон размеров включений гю2, обеспечивающих эффект трансформационного упрочнения, практически невозможно из-за сложности предсказания окончательного размера включений Zr02 после производства.

При 2 = 147 МПа величина ¿(г) составляет 0,0 и 0,3 для у4/20:)- и IVС-матриц соответственно. Более широкий диапазон критических радиусов включений Zr02 в Л/гО.-гматрице по сравнению с И'С-матрицей обусловлен различием коэффициентов липейпого расширения матриц. В случае Л/гОз-матрЩЫ формирующиеся во включениях Zr02 растягивающие термоусадочные напряжения в меньшей степени нивелируют положительную собственную деформацию превращения, что приводит к большей величине изменения плотности энергии деформации -ф по сравнению с И^С-матрицей и, как следствие, к более широкому диапазону критических радиусов. Показано, что более широким диапазоном и большими величинами критических радиусов включений 2т0г будут обладать включения, диспергированные в жесткой матрице с большим коэффициентом линейного расширения (близким к коэффициенту линейного расширения Zr02)■ Показано также, что следствием учета упругого взаимодействия включений является уменьшение обоих критических радиусов включении Zr02 при увеличении объемной концентрации включений, что для критического радиуса Ят согласуется с экспериментальными данными.

В заключении приведены основные результаты и выводы по диссертационной работе.

0,8 0,6

0,4 0,2

0,8 0,6

0,4 0,2

А120з , 1=29 ЗК II

"50 ТОО 150 2, МПа

\УС, Т=2 II

0 40 80 120 160 180 2, МПа

Рис. 9: Зависимости <?(г) мета-стабильных включений 2гС?2 в Л/2О3- и И'С-ыатрицах.

- Развит новый метод исследования устойчивости межфазной границы, основанный на соотнесении деформаций на межфазной границе с границами зоны фазовых переходов. Определены диапазоны радиусов межфазной границы и внешних деформаций, в которых сохраняется устойчивость.

2. В результате решения задачи о фазовых превращениях взаимодействующих включений в дисперсном композитном материале развита новая модель трансформационного упрочнения керамик, в которой помимо обычно учитываемых химической и поверхностной энергий и термоусадочных напряжений вводится дополнительный энергетический барьер, который должен быть преодолен на пути превращения, и учитывается упругое взаимодействие включений. Введение барьера позволило получить адекватную величину диапазона размеров включений ZrO■l, обеспечивающих трансформационное упрочнение в керамиках. Предложены методики определения величины барьера на основе экспериментальных данных. Теоретически показано, что величины критических радиусов уменьшаются с ростом объемной концентрации включений Zr02■ Проведены расчеты диапазонов радиусов частиц для .4/2(9;г и И'С-матриц.

Публикации автора по теме диссертации

Публикации в изданиях из перечня ВАК:

1. R.A. Filippov, А.В. Freidin, E.N. Vilchevskaya. On new phase intermediate layers in nanocomposites as a source of increasing the elastic moduli // Nanosystems: Physics, Chemistry, Mathematics. — 2011— Vol. 2.— №2 — C. 107-118.

2. P.A. Филиппов. Переходные слои в композитах как области новой фазы // Вестник Нижегородского государственного университета пм. Н.И.Лобачевского.- 2011.- Вып. 4(5).- С. 545-546.

3. Е.Н. Вильчевская, Р.А. Филиппов, А.Б. Фрсйдин. О переходных слоях в композитных материалах как областях новой фазы // Известия РАН. Механика твердого тела,— 2013.— №1.— С. 128-159.

4. А.В. Freidin, R.A. Filippov, I. Hussainova, E.N . Vilchevskaya. Critical radius in the effect of transformation toughening of zirconia doped ceramics and cermets // Key Engineering Materials.- 2013,- Vol. 527,- C. 68-73.

Прочие публикации:

5. R.A. Filippov, A.B. Freidin, E.N. Vilchevskaya. On transition layers in nanocomposites as new phase areas // Proc. of the XXXIX Summer School APM-2011 (Advanced Problems in Mechanics). 2011. IPME RAS, St. Petersburg.- P. 158-163.

Подписано в печать 08.11.13 Формат 60x84 1/16 Цифровая Печ. л. 1.0 Тираж 100 Заказ 16/11 печать

Отпечатано в типографии "Фалкон Принт" (197101, г. Санкт-Петербург, ул. Большая Пушкарская, д. 54, офис 2)

Текст научной работы диссертации и автореферата по механике, кандидата физико-математических наук, Филиппов, Роман Александрович, Санкт-Петербург

ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ УЧРЕЖДЕНИЕ НАУКИ ИНСТИТУТ ПРОБЛЕМ МАШИНОВЕДЕНИЯ РОССИЙСКОЙ

АКАДЕМИИ НАУК

На правах рукописи

04201450716 Филиппов Роман Александрович

ФАЗОВЫЕ ПРЕВРАЩЕНИЯ В МАТЕРИАЛАХ С ВКЛЮЧЕНИЯМИ

01.02.04 - Механика деформируемого твердого тела

Диссертация

на соискание ученой степени кандидата физико-математических наук

Научный руководитель:

доктор физико-математических наук

А.Б. Фрейдин

Санкт-Петербург - 2013

Содержание

Стр.

ВВЕДЕНИЕ 4

Глава 1. РАЗВИТИЕ ПЕРЕХОДНОГО СЛОЯ НОВОЙ ФАЗЫ ВОКРУГ ИЗОЛИРОВАННОГО ВКЛЮЧЕНИЯ 15

1.1. Постановка задачи определения равновесной межфазной границы

при бездиффузионном фазовом превращении мартенситного типа . 15

1.2. Равновесный радиус межфазной границы...............19

1.3. Перераспределения плотности энергии деформаций, окружных напряжений и интенсивности касательных напряжений вследствие формирования равновесного слоя новой фазы............23

1.4. Выводы по первой главе.........................27

Глава 2. УСТОЙЧИВОСТЬ СФЕРИЧЕСКОЙ МЕЖФАЗНОЙ ГРАНИЦЫ 28

2.1. Устойчивость межфазной границы по отношению к возмущению радиуса..................................28

2.2. Устойчивость межфазной границы по отношению к осесимметрич-ным возмущениям............................31

2.3. Устойчивость межфазной границы по отношению к произвольным возмущениям...............................35

2.4. Устойчивость межфазных границ, возникающих в однородном теле, и индуцированных неоднородностью................40

2.5. Выводы по второй главе.........................44

Глава 3. ДИСПЕРСНЫЙ КОМПОЗИТ С хМАТРИЦЕЙ ПРЕТЕРПЕВАЮЩЕЙ ФАЗОВЫЕ ПРЕВРАЩЕНИЯ 46

3.1. Метод эффективного поля .......................47

3.2. Устойчивость равновесной микроструктуры .............52

3.3. Эффективные модули..........................56

3.4. Выводы по третьей главе........................58

Глава 4. ТРАНСФОРМАЦИОННОЕ УПРОЧНЕНИЕ КЕРАМИК. ДИАПАЗОН КРИТИЧЕСКИХ РАЗМЕРОВ ВКЛЮЧЕНИЙ гЮ2 60

4.1. Постановка задачи определения диапазона критических радиусов

с учетом энергетического барьера...................65

4.2. Вычисление изменения энергии деформаций.............69

4.3. Выбор параметров модели и вычисление напряжений........71

4.4. Оценка величины энергетического барьера..............74

4.5. Результаты расчетов критических радиусов..............79

4.6. Выводы по четвертой главе.......................84

ЗАКЛЮЧЕНИЕ 85

СПИСОК ЛИТЕРАТУРЫ 88

ВВЕДЕНИЕ

Актуальность

Работа посвящена развитию моделей механики, описывающих фазовые превращения в материалах с включениями при термомеханичеекпх воздействиях.

Исследования фазовых превращений в деформируемых телах представляют актуальное направление современной механики материалов.

Применительно к дисперсным композитным материалам следствием фазовых (структурных) превращений вокруг частиц является формирование переходных слоев измененного материала матрицы и, как следствие, изменение деформационно-прочностных свойств композита. Примерами таких композитов являются сплавы на основе титана, которые испытывают мартенситные превращения и в которых, при определенных термомеханических условиях обработки, выделяются частицы, не претерпевающие фазовое превращение и способные инициировать его в матрице вокруг себя [39]. Также одним из требующих объяснения фактов является отмеченная в литературе возможность изменения деформационно-прочностных свойств нанокомпозитов (увеличения модулей упругости, упрочнения) при малой объемной доле включений [50,94,98,101].

Фазовые превращения могут также происходить внутри частиц. Примером являются фазовые превращения частиц диоксида циркония в керамических композитах, приводящие к эффектам трансформационного упрочнения [25,47,56,65,81].

воздействиями и фазовыми превращениями.

Первая задача - задача описания фазовых превращений, идущих в матрице вокруг включений. В результате формирования областей новой фазы в матрице возникают переходные слои вокруг включений, увеличивающие эффективный объем включений и, как следствие, возрастает влияние объемной доли включений на упругие свойства композитного материала.

В ряде работ [54,68,93,98,103] было обращено внимание на влияние на механические свойства композита состояния границ, разделяющих включения и матрицу. Но классические теории механики композитных материалов исходят из непрерывности перемещений на границах, т. е. идеальной адгезии, и даже в этом случае неспособны объяснить изменение модулей упругости. Остается признать, что влияние на деформационно-прочностные свойства композита может оказывать не только состояние границы раздела, но и состояние окружающего частицу материала [57,68,93,103].

Влияние переходных слоев вокруг частиц на эффективные упругие свойства нанокомпозитов обсуждалось с использованием различных приближений в работах [3,10,48,49]. В работе [97] проведено сравнение различных методов оценки эффективных свойств паиокомпозита с переходными слоями: аппроксимация без учета взаимодействия включений, с использованием дифференциальной схемы Мори - Танака, с помощью оценок Хашина - Штрикмана для композита и метода эффективного ноля. Показано, что при большой разнице в упругих свойствах включения и матрицы аппроксимация без учета взаимодействия включений дает существенную ошибку. В работах [41,95] исследовалось влияние микроструктуры на эффективные физико-механические свойства, в частности на поведение материала под действием термических напряжений. В работах [7,23,87-89] развивается градиентная модель межфазного слоя для описания свойств материала в окрестности жесткого включения. Показывается, что градиентные модели могут достаточно адекватно моделировать эффекты аномального увеличения эф-

фективных свойств материала и уменьшения концентрации напряжений в окрестности жесткой фазы, т. е. эффект разгрузки матрицы в окрестности жесткого включения.

В приведенных выше примерах работ, посвященных исследованию влияния переходных слоев на механические свойства композитов, это влияние обсуждалось безотносительно механизмов формирования переходных слоев.

В диссертационной работе исследуется возможность связать влияние малых объемных долей частиц с формированием переходных слоев, возникающих вокруг частиц вследствие фазового превращения материала матрицы. Отметим, что фазовое превращение является очевидным примером самопроизвольного термодинамически выгодного процесса, идущего в материале постольку, поскольку этот процесс уменьшает энергию тела.

При описании фазовых превращений с позиции механики деформируемого твердого тела можно выделить два подхода. В первом подходе в феноменологическую модель, описывающую фазовое превращение, вводятся дополнительные параметры, характеризующие особенности микроструктуры и протекающих процессов в среднем (см. работы А.Е. Волкова [4], В.А. Лихачева [20], В.Г. Малинина [21], Г.А. Малыгина [24], A.A. Мовчана [26-28] , А.И. Разова, К. Баттачарьи [96], К. Лекслен [83-85]). Формулируются определяющие соотношения для этих параметров. Данный подход позволяет выявить важные особенности фазового превращения, исключая при этом явное рассмотрение межфазных границ, локальных нолей напряжений и деформаций.

Второе направление основано на рассмотрении фазовых превращений с учетом условий равновесия на границе фаз деформированного материала, а также включает детальное описание возникающих под напряжением двухфазных структур [9,11,16,29,30,32,34,37,40,58,72,75].

Задача о двухфазной конфигурации исследуется в русле второго подхода, согласно которому равновесная межфазная граница должна удовлетворять дополнительному термодинамическому условию, что, в свою очередь

приводит к зависимости размеров области новой фазы от напряжений. В результате объемная доля получающихся составных включений становится управляемым параметром, зависящим от термомеханических воздействий на технологической стадии изготовления композита.

Керамики обладают многими полезными свойствами - высокой твердостью и температурой плавления, низкой тепло- и электропроводностью, слабой химической активностью в агрессивных средах Серьезным недостатком керамик, ограничивающим области их применения, является хрупкость. Например, керамика на основе оксида алюминия А^О^, имея высокую твердость, является хрупким материалом с низкой трещипостойкостью (см., напр., [47,77]). Прочность и ударная вязкость керамик могут быть увеличены в результате диспергирования в них частиц диоксида циркония Zr02 [25,33,47,77,81]. Диоксид циркония может претерпевать фазовое превраще-

ние мартенситного типа. При переходе из высокотемпературной аустепитной в низкотемпературную мартенситную фазу происходит превращение тетрагональной кристаллической решетки (£ — Zr02) в моноклинную (т — Zr02), которое сопровождается положительной объемной (3-5%) и сдвиговой (1315%) деформациями превращения. Если фазовое превращение частиц инициируется полем напряжений, создаваемым трещиной, то перераспределение напряжений, вызванное увеличением размеров частиц, и энергозатраты на микрорастрескивание материала матрицы вокруг перешедших в новое фазовое состояние частиц приводят к замедлению или блокированию роста трещины. В результате прочность керамики увеличивается, в чем и заключается эффект трансформационного упрочнения.

Для эффективного трансформационного упрочнения размеры включений должны принадлежать определенному диапазону, а именно быть такими, чтобы включения оставались в метастабилыюм высокотемпературном аусте-нитном фазовом состоянии, несмотря на остывание керамики до эксплуатационных температур, но переходили в низкотемпературное мартенситное состояние в поле напряжений распространяющейся трещины [47,53,56,65,77,81]. Обычно, начиная с пионерских работ Р. Гарвье [65,66] и А. Эванса [56] критические размеры включений 2г02 определяются в результате сравнения энергий Гиббса керамики с включениями, находящимися в разных фазовых состояниях, с учетом изменений химической энергии, энергии деформаций и поверхностной энергии. Однако без учета дополнительных факторов в критерии превращения превалирующая роль изменения химической энергии приводит к незначительному влиянию напряжений на критический размер включений и, как следствие, очень узкому диапазону размеров включений, обеспечивающих трансформационное упрочнение, что ставит под сомнение адекватность такой постановки задачи.

В диссертационной работе при сравнении плотностей энергий Гиббса композитного материала с частицами в разных фазовых состояниях предлага-

ется принять во внимание, что уменьшение энергии Гиббса является только необходимым условием возможности фазового превращения. Реализация этой возможности связана с кинетикой фазового превращения, зависящей от энергетических барьеров, которые должны быть преодолены на пути превращения частицы (см., напр., [22]). Одним из таких барьеров может быть активационный барьер, связанный с перестройкой кристаллической решетки, которая затрудняется при введении стабилизирующих добавок. В литературе (см. [15, 73, 79,106]) также отмечалось, что для осуществления фазового превращения необходимы дополнительные энергозатраты на образование зародышей новой фазы, причем нуклеационный барьер, связанный с этими энергозатратами, зависит от плотности центров зародышеобразования и стабилизирующих добавок. Но несмотря на то что на существование энергетических барьеров указывалось и ранее, модели, явным образом учитывающие эти барьеры, развиты не были.

В настоящей работе исследуется, как учет энергетического барьера в критерии фазового превращения увеличивает влияние изменения потенциальной энергии деформаций па критические размеры частиц, что, в свою очередь, приводит к более широкому диапазону их размеров. Величина барьера рассматривается как дополнительный параметр модели, для оценки которого предлагаются два способа. В первом способе используются результаты экспериментального определения относительной ширины диапазона критических радиусов, а во втором - экспериментально полученная зависимость критического радиуса от содержания стабилизирующей добавки. Кроме того, в отличие от проведенных ранее исследований, при вычислении изменения объемной плотности потенциальной энергии деформаций в диссертационной работе учитывается вклад упругого взаимодействия частиц друг с другом. Для демонстрации влияния термоупругих параметров матрицы приводятся расчеты критических размеров частиц в зависимости от энергетического барьера для А/2О3- и И^С-матриц, различающихся коэффициентами линейного

расширения и модулями упругости. Определяется диапазон размеров включений ^гОг, обеспечивающих эффективное трансформационное упрочнение.

Целью диссертационной работы является разработка и реализация моделей, описывающих влияние внешних термомеханических воздействий на фазовые превращения в материалах с включениями.

Задачами работы являются:

Научную новизну диссертации представляют следующие положения, выносимые на защиту:.

- Показано, что формирование переходного слоя повой фазы приводит к сбросу энергии деформации в матрице и, как следствие, к уменьшению концентрации напряжений. Определены ограничения на параметры материала, при выполнении которых возникновение переходного слоя приводит к сбросу энергии деформаций во всех компонентах композита.

2. На основе новой модели трансформационного упрочнения, учитывающей энергетический