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

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

Кондратьев Никита Сергеевич

УПРУГОВЯЗКОПЛАСТИЧЕСКАЯ МОДЕЛЬ ДЛЯ ОПИСАНИЯ ДЕФОРМИРОВАНИЯ МНОГОФАЗНЫХ ПОЛИКРИСТАЛЛОВ В НЕИЗОТЕРМИЧЕСКИХ УСЛОВИЯХ

05.13.18 — Математическое моделирование, численные методы и комплексы программ

АВТОРЕФЕРАТ

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

Пермь — 2014

005560421

005560421

Работа выполнена в ФГБОУ ВПО «Пермский национальный исследовательский политехнический университет»

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

Трусов Петр Валентинович, доктор физико-математических наук, профессор

Официальные оппоненты:

Романова Варвара Александровна, доктор физико-математических наук, ведущий научный сотрудник, Институт физики прочности и материаловедения СО РАН, г. Томск Гаришин Олег Константинович, доктор физико-математических наук, старший научный сотрудник, Институт механики сплошных сред УрО РАН, г. Пермь

Ведущая организация:

Институт машиноведения УрО РАН, г. Екатеринбург

Защита состоится «27» января 2015 г. в 16:00 на заседании диссертационного совета Д 212.188.08 при ФГБОУ ВПО «Пермский национальный исследовательский политехнический университет» по адресу: 614990, г. Пермь, Комсомольский пр., д. 29, ауд. 423.

С диссертацией можно ознакомиться в научной библиотеке и на сайте Пермского национального исследовательского политехнического университета

Автореферат разослан «10» декабря 2014 г.

Ученый секретарь диссертационного совета Д 212.188.08, кандидат физико-математических наук, доцент У'" А. И. Швейкин

Перечень основных обозначений, сокращений и символов

Обозначения: ЛГ (Ны) - число систем скольжения (систем двойникования) для рассматриваемого типа решетки, Ку - число фасеток

рассматриваемого кристаллита, у'*' - скорость сдвига в к-й системе скольжения, 0,9 - абсолютная температура (К) на макроуровне и мезоуровне, - сдвиговое напряжение, критическое

напряжение сдвига в к-й системе скольжения (системе двойникования), Я(-) - функция Хэвисайда, Ь'*'^*'), ¿^'(п^.') - единичные векторы

направления скольжения (двойникования) и нормали плоскости скольжения (двойникования) к-й системы скольжения (двойникования) в актуальной конфигурации, сС - внешняя нормаль границы с соседним

зерном п, V(•),V(•) - операторы Гамильтона (набла-операторы) в отсчетной и актуальной конфигурации, Б, Пе, Б'", Л'11 — тензор деформации скорости, его упругая, пластическая и температурная составляющие (макроуровень), с1, (Iе, с1'л, &л - тензор деформации скорости, его упругая, пластическая и температурная составляющие (мезоуровень), Ш<1) - ориентационный тензор к-й системы скольжения, г*' - ориентационный тензор к-й системы двойникования, \У (\у) -тензор вихря макроуровня (мезоуровня), Е (о) - тензор напряжений Коши макроуровня (мезоуровня), £2(ю) - (антисимметричный) тензор спина, характеризующий вращение подвижной системы координат на макроуровне (мезоуровне), П (п) - тензор четвертого ранга упругих свойств макроуровня (мезоуровня), у0 - скорость сдвига при сдвиговых напряжениях, равных критическим, тпг - показатель скоростной чувствительности скольжения и двойникования, (<3Л|,) - энергия активации скольжения дислокаций (двойникования), 0„, бдду

критические значения температуры и энергии дефекта упаковки материала, при превышении которых следует учитывать механизм разупрочнения за счет возврата, Ц„,ЦГС - безразмерный параметр, устанавливаемый в ходе процедуры идентификации модели, <21, -энергия активации диффузии, ®*гс - интенсивности накопленных

неупругих деформаций, критическая величина интенсивности неупругих деформаций и критическая температура, после которой необходимо учитывать процесс рекристаллизации, аи'!) - матрица параметров зернограничного упрочнения; т40) - напряжения зернограничного упрочнения, к - газовая постоянная.

Сокращения: ПО - представительный объем, СС - система скольжения, СД - система двойникования, ВДС - внутризеренное дислокационное скольжение, РД - решеточная дислокация, ДОН дислокация ориентационного несоответствия.

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

Актуальность работы. Многочисленные экспериментальные и теоретические исследования интенсивного упругопластического деформирования поликристаллов свидетельствуют о существенной перестройке мезо- и микроструктуры материала, которая в значительной степени определяет поведение материала на макроуровне. Вследствие этого в настоящее время актуальной является задача построения моделей неупругого деформирования моно- и поликристаллов, способных описывать изменения (эволюцию) дефектной структуры. Первые попытки построения таких моделей были предприняты еще в середине XX века в работах Дж. Тейлор, Дж. Бишоп, Р. Хилл, Т.Г. Линь и др. Ощутимый вклад в развитие указанного направления внесли советские и российские ученых: Я.Д.Вишняков, О.А.Кайбышев, Р.З.Валиев, В.А.Лихачев, В.В.Рыбин, Т.Д.Шермергор, С.Д.Волков, В.Е.Панин и др.

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

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

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

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

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

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

Научная новизна работы заключается в следующем:

- предложена структура двухуровневой, включающей макро- и мезоуровни, модели физической теории пластичности для описания поведения многофазных материалов, в рамках которой: 1) предложен способ описания взаимодействия дислокаций с границами зерен, учитывающий сдвиги в соседних кристаллитах; 2) предложен способ описания процессов разупрочнения при высокотемпературном деформировании многофазных материалов — динамической рекристаллизации и возврата;

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

- разработаны алгоритмы и комплекс программ, реализующих предлагаемую математическую модель.

На защиту выносятся:

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

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

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

-математическое описание разупрочнения двухфазных материалов при высокотемпературном деформировании в двухуровневой физической модели;

-подход к идентификации параметров модели, учитывающей эволюцию мезоструктуры;

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

-Комплекс программ, реализующий предлагаемый алгоритм, и результаты моделирования.

Практическая ценность работы заключается в возможности использования предлагаемой модели и разработанного комплекса программ для анализа процессов интенсивной пластической деформации, учитывающих изменение внутренней микро- и мезоструктуры, а также - прогнозирования физических и механических характеристик материалов на макроуровне. Получено государственное свидетельство о регистрации программы для ЭВМ №2013619775 от 14.10.2013 [9].

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

Апробация работы. Результаты диссертационной работы докладывались на конференциях: ХУШ-ХХШ Всероссийских школах конференциях молодых ученых «Математическое моделирование в естественных науках» (Пермь, ПНИПУ, 2009-2013), XVII и XVIII Зимних школах-конференциях по механике сплошных сред (Пермь, 2011, 2013); VII Российской научно-технической конференции «Механика микронеоднородных материалов и разрушение» (Екатеринбург, 2012); VII Международной конференции «Микромеханизмы пластичности, разрушения и сопутствующих

явлений» (Тамбов, 2013); первой Всероссийской научной конференции молодых ученых с международным участием «Перспективные материалы в технике и строительстве» (Томск, 2013). Работа полностью докладывалась и обсуждалась на семинарах кафедры математического моделирования систем и процессов ПНИПУ (рук. проф. П.В.Трусов), Института механики сплошных сред УрО РАН (рук. академик РАН В.П. Матвеенко), кафедры механики композиционных материалов и конструкций ПНИПУ (рук. проф. Ю. В. Соколкин).

Публикации. Основные результаты исследований диссертационной работы представлены в 20 публикациях, наиболее значимые из которых перечислены в списке [1-8], из них - 5 статей ([1,2,6-8]) опубликованы в журналах, рекомендованных ВАК.

Личный вклад автора - постановка задачи (совместно с научным руководителем), реализация программ на ЭВМ, проведение вычислений, анализ результатов.

Структура и объем работы. В диссертационную работу включено введение, список сокращений и основных обозначений, четыре главы, заключение, список цитированной литературы. Диссертация изложена на 169 страницах, содержит 21 рисунок, 3 таблицы и 171 наименование библиографического списка.

Основное содержание работы

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

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

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

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

определяющие соотношения (уравнения состояния), эволюционные уравнения и замыкающие уравнения.

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

скольжения (СС) в кристаллите и являются взаимно ортогональными. Критерием активации к-й системы скольжения (СС) является закон Шмида:

т^п?>.0.ь:*>=<£>. (1)

Для описания скоростей сдвигов по СС с учетом влияния температуры на процесс скольжения используются вязкие соотношения вида:

Г т Ч""'

х(4) Ч " У

ехр^^ЧО- (2)

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

т<*> = п<*> . о . ¿(*> = т<*> ич

По аналогии с соотношением для скоростей сдвига по механизму ВДС (1) вводится соотношение скорости изменения объемной доли двойников:

■(*).

т( ) I Г С

Чк в/"-' (4)

о, т«<0.

В п. 2.3 рассматриваются вопросы связи параметров мезоуровня с «родственными» характеристиками макроуровня; исследуется согласование определяющих соотношений различных масштабных уровней. Связь соседних масштабных уровней осуществляется посредством включения в определяющие соотношения

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

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

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

В п. 3.2 рассматривается описание упрочнения СС поликристаллических материалов за счет границ кристаллитов в физических многоуровневых моделях. Полагается, что граница представляет собой двумерную специфическую область, отделяющую различные однородные части кристалла (зерна, фазы, двойники). Данный дефект поликристалла является труднопреодолимым препятствием для мобильных решеточных дислокаций (РД), что обусловлено в первую очередь резкой сменой ориентаций СС при переходе через границу (А.Н. Орлов, В.Н. Перезвенцев и др.). Принимается следующий механизм взаимодействия подвижных дислокаций с границей кристаллита: РД рассматриваемого кристаллита переходит в энергетически более выгодную СС соседнего (А.Н. Орлов, В.В.Рыбин), оставляя в границе дислокацию ориентационного несоответствия (ДОН). Следующая решеточная дислокация, скользящая по той же СС кристаллита, будет испытывать дополнительное сопротивление за счет поля упругих напряжений ранее образовавшейся ДОН. Предлагается математическое описание увеличения критических напряжений за счет взаимодействия дислокаций с границей кристаллита.

Граница кристаллита аппроксимируется плоскими участками (фасетками), к-я фасетка границы обозначается Sí и определяется нормалью q*. Рассмотрению подлежат СС, решеточные дислокации которых выходят на данную фасетку границы; определяется СС в соседнем кристаллите, в которую перейдет дислокация из рассматриваемой. Для этого предлагается энергетический критерий, в основу которого закладывается минимум «несовместности поверхностных деформаций», обусловленных пластическими сдвигами в соседних кристаллитах [5]. Вектор Бюргерса ДОН определяется разностью векторов Бюргерса решеточных дислокаций, вклад в

сопротивление вносит составляющая, лежащая в плоскости границы, поэтому рассматривается краевая составляющая этой дислокации. В предположении, что ДОН можно отнести к рассматриваемому кристаллиту, получены качественные зависимости скорости изменения критических напряжений на СС [4].

При выводе эволюционного уравнения для критических напряжений сдвига, обусловленных влиянием границы на подвижные дислокации, использовались гипотезы, приведенные ниже. Согласно соотношению Орована, скорость поступления мобильных дислокаций к границе пропорциональна скорости сдвига у^ для рассматриваемой СС (с учетом геометрии фасеток кристаллита), в силу чего скорость образования ДОН можно считать пропорциональной у^. Принимается, что зависимость от геометрических характеристик фасеток характеризуется отношением площади части границы (фасетки границы), пересекаемой рассматриваемой СС, ко всей площади границы 5 текущего кристаллита и мерой разориентации СС текущего и соседнего кристаллита [5]. Скорость изменения добавочного

критического напряжения сдвига т^, обусловленная ДОН, для данной

_/-й СС рассматриваемого /-го кристаллита, дислокации которой взаимодействуют с тя-ым кристаллитом через А>ую фасетку границы, запишется в виде [4,7]:

Касательные напряжения от текущей ДОН учитываются для всех СС (по отношению к которым напряжения ДОН являются барьерными); скорость изменения барьерных напряжений ДОН по всем фасеткам границы к-\,К} (одновременно по всем соседним кристаллитам т),

опуская в левой части (5) индекс текущего кристаллита /', тогда записывается в виде:

В п. 3.3 рассматриваются механизмы разупрочнения в материалах с несколькими фазами. Приводятся соотношения, описывающие уменьшение критических напряжений на СС за счет динамического возврата и рекристаллизации [1]:

(5)

(6)

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

В п. 4.1 приводится система уравнений двухуровневой математической модели, описывающей неупругое деформирование двухфазного поликристалла в процессах глубоких пластических деформаций: задается схема деформирования на макроуровне, которая передается на мезоуровень (кристаллит, фаза, зерно) посредством расширенной гипотезы Фойгта, предполагающей равенство градиентов перемещения в рамках ПО. Определяющим соотношением обоих масштабных уровней (макро и мезо) является закон Гука в скоростной релаксационной форме; на макроуровне тензор упругих свойств и кососимметричный тензор, описывающий скорость ротации ПО поликристалла определяются осреднением аналогичных величин с мезоуровня; скорость пластических деформаций макроуровня определяется из нижележащего мезоуровня осреднением скоростей сдвигов по активным СС и изменения объемной доли двойников по системам двойникования (СД), для нахождения которых, в свой черед, применяется кинематическое соотношение вязкого типа. Делается предположение о квазистатическом характере процессов деформирования, но которые протекают при средних и высоких гомологических температурах. Проверяется математическая корректность и замкнутость полученной системы уравнений.

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

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

Параметры, требующие определения, условно делились на два типа; часть параметров предполагалась известными го работ других исследователей, либо находилась с помощью простых процедур, параметры второго типа определяются из решения следующей задачи оптимизации [2]:

найти оптимальное управление и = при заданных и

фиксированных параметрах состояния

х = {п1Ш,п1122,п1212, х[% у0,/0, ти5, ты, к{т}, А{Ш)], доставляющее минимум функционалу:

2(а?(х,п)-о?)2-»тш,

(8)

при ограничениях типа равенства [2]:

о = с: (п - п'") + ю а - а ю, с! " = £ + 2 /''Ч,*'

Уу = УУ,

1 V, 7 '

к

т = 1( Уут - У у) +у ™ (п<*>Ь<*> - Ь 1*>п<*>),

^ = (1 -+ /£ «т^, =у/Е

+(1-9,

А»=йг

1=1 х„-1

V 1о "

/(Ц=

г,

О,

Тт

>0,

(9)

Здесь а^ - значение напряжений, определяемое по аппроксимированной экспериментальной диаграмме нагружения при интенсивности деформации е; а^ - интенсивность напряжений, рассчитанная с помощью модели; индекс / обозначает номер точки на аппроксимированной экспериментальной диаграмме нагружения.

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

используемые на стадиях идентификации и верификации модели. С целью определения материальных параметров и проверки модели рассмотрен эксперимент по одноосной осадке материала в различных кристаллографических направлениях: на этапе идентификации - в направлении [100], верификации- [011]. Можно отметить удовлетворительное соответствие между модельными и натурными данными.

« 5 о и

Ш

Г .......

1

т л ...............

Ж 1 * /

1

£ вм 4М Ш §£* ф. 5 11ктем«|«вшк1Ь деф»|ВДЖ№£1 с

Цтскашноси* дсфорЛыддИ

Рис. 1. Диаграмма сжатия монокристалла тантала в направлении [100] (а, стадия идентификации), диаграммы сжатия монокристалла тантала в направлении [011] (б,

стадия верификации);-----численные результаты ,--эксперимент

В п. 4.4 приводятся основные результаты численного моделирования деформирования осадкой, стесненной осадкой и простым сдвигом двухфазных материалов [3, 6]. На рис. 2 показана «а-е» диаграмма для указанных нагружений дуплекс стали с учетом упрочнения за счет границ кристаллита и без него.

ГЛа

0.20- •

0ЛС-о.оз-

0.04-

<5.84- -8.»

0.02 0.« 0.06 ОМ 0.10 Ей ГПа

о.ог е.64 о.ов 0.08 о.ю

0.02 0.С4 с.Оо е.оа о.ю

Рис. 2. Диаграмма нагружения при простом сдвиге (--с учетом

упрочения за счет границ кристаллитов,-----без учета): а) простой сдвиг, б)

осадка, в) стесненная осадка

На рис.3 показана эволюция критических напряжений на одной из активных СС выделенного из поликристалла зерна и накопленного сдвига в случаях с учетом упрочнения за счет границ и без него.

• ГПа

б)

0.05 е. В С. 5 0.2Й " 0 05 С. й

Рис. 3. Эволюция критических напряжений на одной из активных СС кристаллита аустенитной фазы и накопленного сдвига в случае с учетом упрочнения за счет границ (пунктирная линия) и без него (сплошная линия) В п. 4.5 представлен анализ результатов неизотермического деформирования многофазных материалов. На рис.4 изображены диаграммы одноосного сжатия вдоль оси О.т1 фиксированной лабораторной системы координат ПО дуплекс стали марки БАБ 2507 при температурах 800°С, 900°С и 700°С; маркерами «•» отмечены экспериментальные данные [2]

о.,МПа

350 300 250 200 150 100 50

а)

а., МПа

б)

0,05 0,10 0,15 0,20 0,25

Рис. 4. Диаграмма одноосного сжатия представительного

макрообъема дуплекс стали марки 8АР2507 при 9, °С: 800 (а), 900 (б) и 700 (в); маркерами «•» отмечены данные натурного эксперимента

0,05 0,10 0,15 0,20 0,25 е„

На рисунке 5 приведены гистограммы распределения интенсивности напряжений для зерен аустенитной (средние значение -390 МПа, дисперсия - 22 МПа) и ферритной (средние значение -203 МПа, дисперсия - 19 МПа) фаз в момент окончания процесса одноосного сжатия ПО дуплекс стали марки 8АР2507 при температуре 800°С.

ст., МПа 500 ;

400 ■

а)

ЛЕ

б)

250 300 350 400 о„ МПа 140 160 180 200 220 о,,МПа

Рис. 5. Гистограмма распределения интенсивности напряжений в аустенитной (а) и ферригной (б) фазах в момент окончания одноосного сжатия дуплекс стали SAF2507

при температуре 800°С В заключении сформулированы основные результаты диссертационной работы:

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

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

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

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

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

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

Список работ, опубликованных автором по теме диссертации

1. Трусов П.В., Кондратьев Н.С. Двухуровневая модель для описания неизотермического деформирования двухфазных поликристаллов // Вычислительная механика сплошных сред = Computational continuum mechanics. - 2014. - Т. 7. № 2. - С. 181-199. (перечень ВАК)

2. Кондратьев Н.С., Трусов П.В. Математическая модель для описания деформирования ОЦК-монокристаллов, учитывающая двойникование // Вычислительная механика сплошных сред =

Computational continuum mechanics. - 2011. - T.4, № 4. - C. 20-33.(перечень ВАК)

3. Кондратьев H.C., Трусов П.В. Моделирование неупругого деформирования поликристаллических материалов с учетом упрочнения за счет границ кристаллитов// Вестник Пермского университета. Серия: Физика, 2012. -№ 4 (22). - С. 92-100.

4. Кондратьев Н.С., Трусов П.В. Описание упрочнения систем дислокационного скольжения за счет границ кристаллитов в поликристаллическом агрегате // Вестник Пермского национального исследовательского политехнического университета. Механика. - Пермь: Изд-во ПНИПУ, 2012. - № 3. -С. 79-98.

5. Кондратьев Н.С., Трусов П.В. О мере разориентации систем скольжения соседних кристаллитов в поликристаллическом агрегате // Вестник Пермского национального исследовательского политехнического университета. Механика. - Пермь: Изд-во ПНИПУ, 2012. -№2. -С. 112- 127.

6. Кондратьев Н.С., Трусов П.В. Физическая модель неупругого деформирования двухфазных поликристаллов // Вестник Тамбовского университета. Серия: Естественные и технические науки. - 2013. - Т. 18. № 4-2. - С. 1873-1874. (перечень ВАК)

7. Трусов П.В., Кондратьев Н.С. Описание неупругого деформирования двухфазных поликристаллических материалов // Деформация и разрушение материалов. - 2013. - № 6. - С. 8-15. (перечень ВАК)

8. Швейкин А.И., Бразгина О.В., Кондратьев Н.С. Моделирование эволюции структуры ГЦК-, О ЦК- и ГПУ-поликристаллов при неупругом деформировании // Вестник Нижегородского университета им. Н.И. Лобачевского. - Н. Новгород: Изд-во ННГУ им. Н.И. Лобачевского, 2011. - №4. Часть 4. - С. 1859-1861. (перечень ВАК)

9. Кондратьев Н.С., Трусов П.В., Швейкин А.И. «Реализация двухуровневой модели неупругого деформирования ГЦК-поликристаллов для применения в пакете Abaqus» («Модель ГЦК-поликристалла для Abaqus»). - Свидетельство о государственной регистрации программы для ЭВМ№ 2013619775 от 14.10.2013.

Подписано в печать 27.11.2014.

Формат 60 х 90/16. Набор компьютерный.

_Усл. печ. л.1. Тираж 100 экз. Заказ № 1609/2014. _

Отпечатано с готового оригинал-макета в типографии издательства Пермского национального исследовательского политехнического университета.

Адрес: 614990, г. Пермь, Комсомольский пр-т, 29, к. 113. Тел. (342) 219-80-33.