УЧЕНЫЕ ЗАПИСКИ КАЗАНСКОГО УНИВЕРСИТЕТА. СЕРИЯ ФИЗИКО-МАТЕМАТИЧЕСКИЕ НАУКИ
2019, Т. 161, кн. 4
С. 509-525
ISSN 2541-7746 (Print) ISSN 2500-2198 (Online)
УДК 539.3
doi: 10.26907/2541-7746.2019.4.509-525
КОНТИНУАЛЬНАЯ МОДЕЛЬ ПОВРЕЖДЁННОЙ СРЕДЫ ДЛЯ ОПИСАНИЯ ПРОЦЕССА РАЗРУШЕНИЯ
ПРИ ПОЛЗУЧЕСТИ
Научно-исследовательский институт механики Национального исследовательского Нижегородского государственного университета им. Н.И. Лобачевского, г. Нижний Новгород, 603950, Россия
Рассмотрены основные закономерности процессов деформирования и деградации начальных прочностных свойств конструкционных материалов (металлов и их сплавов) по механизму длительной прочности. С позиции механики повреждённой среды (МПС) развита математическая модель, описывающая процессы неупругого деформирования и накопления повреждений при ползучести. Модель МПС состоит из трёх взаимосвязанных частей: соотношений, определяющих неупругое поведение материала с учётом зависимости от процесса разрушения; уравнений, описывающих кинетику накопления повреждений; критерия прочности повреждённого материала.
Разработана экспериментально-теоретическая методика нахождения материальных параметров развитых определяющих соотношений МПС. Определение материальных параметров и скалярных функций модели МПС, оценка её достоверности и определение границ применимости проведены на основе анализа процессов деформирования и разрушения лабораторных образцов в условиях мягкого нагружения (контролируемых напряжений).
Приведены результаты экспериментальных исследований кратковременной ползучести жаропрочного сплава ВЖ-159 при постоянных значениях температуры и разных значениях задаваемых в образцах напряжений.
Проведены численные исследования процесса деформирования и накопления повреждений и выполнено сравнение полученных численных результатов с данными натурных экспериментов. Результаты сопоставления расчётных и экспериментальных данных позволили сделать вывод о достоверности предложенных определяющих соотношений МПС при деградации начальных прочностных свойств конструкционных материалов по механизму длительной прочности.
Ключевые слова: нестационарная ползучесть, длительная прочность, моделирование, определяющее соотношение, механика повреждённой среды, температура, повре-ждённость, материальный параметр, численный эксперимент, натурный эксперимент
Ресурс элементов конструкций, работающих в условиях повышенных температур Т (Т > 0.5 Тпл, Тпл - температура плавления) и циклических механических воздействий, определяется в основном физическими процессами деградации начальных прочностных свойств конструкционных материалов (процессами малоцикловой усталости и накопления повреждений в результате ползучести), которые приводят к одному из наиболее опасных типов разрушения - хрупкому разрушению конструкций, изготовленных из пластичных материалов.
И.А. Волков, Л.А. Игумнов, С.Ю. Литвинчук
Аннотация
Введение
Существующие на сегодняшний день нормативные методы оценки ресурса элементов конструкций не учитывают реальных процессов, протекающих в материале. Упругий расчёт, используемый при нормативном подходе, не позволяет учесть реальные характеристики вязкопластического деформирования материала, от которых в значительной степени зависит ресурс элементов конструкций. В общем случае прочность конструкций должна учитывать время и историю нагружения. Как следствие этого, критерий разрушения будет тесно связан с определяющими соотношениями, описывающими процесс разрушения [1].
В связи с этим становится необходимой разработка новых методов оценки ресурса элементов конструкции на базе соответствующих уравнений термовязкопла-стичности, уравнений накопления повреждений и критериев разрушения со всесторонним их обоснованием посредством проведения соответствующих натурных и численных экспериментов на лабораторных образцах и численного анализа процессов деформирования и разрушения элементов конструкций в эксплуатационных условиях [1, 2].
Зависимость долговечности материала от длительности (частоты) цикла и от наличия и длительности выдержек в цикле, протекающих на фоне повышенных постоянных или меняющихся температур, обусловлена развивающимися деформациями ползучести и вызванными этими деформациями повреждениями, которые в противоположность усталостным повреждениям, развивающимся по телу зерна, накапливаются по границам зерен и приводят к межкристаллитному разрушению.
Наиболее распространенным методом исследования длительной прочности является построение экспериментальной кривой: «растягивающее напряжение / время до разрушения» в экспериментах по одноосному растяжению образцов при постоянной температуре [3]. На кривых ползучести имеется три характерных участка. Первый из них соответствует неустановившейся ползучести, причём скорость деформации ползучести убывает; вторым является участок установившейся ползучести, когда скорость деформации ползучести практически постоянна. С точки зрения исследования характеристик процесса накопления повреждений наибольший интерес представляют время начала третьего участка, скорость ползучести на установившемся участке, закон изменения скорости ползучести на третьем участке, время до разрушения, и деформация ползучести к моменту разрушения. Наличие третьего неустановившегося участка на кривой ползучести вызвано влиянием накопленной повреждённости на деформационные характеристики материала (на скорость ползучести), а закон изменения скорости ползучести па этом участке определяется скоростью накопления повреждения в материале на этом участке. Начало третьего участка есть граница начала влияния накопленной по-вреждённости на деформационные характеристики материала.
Три стадии деформации ползучести являются ни чем иным, как проявлением на макроуровне последовательных изменений в материале. После продолжительного начального периода с преобладанием вязкого деформирования характерное для третьей стадии быстрое увеличение скорости деформации является прямым следствием ускоренного роста микродефектов с последующим их слиянием, приводящим к разрушению при ползучести. Применительно к конструкционным материалам очень важен вопрос о взаимосвязи процессов деформирования и разрушения, поскольку наличие такой взаимосвязи даёт возможность искать пути оценки повреждённости материала при ползучести по деформационным характеристикам. Экспериментальные исследования показали наличие чётко выраженной связи кинетических процессов деформирования и накоплениями повреждений при ползучести: существует линейная связь между деформацией ползучести на всех стадиях и объёмной доли микродефектов, на третьей стадии наблюдается прямо пропорциональная зависимость между приростом объёмной доли микропор и увеличением
скорости ползучести [4, 5]. Разделить процессы зарождения и роста несплошно-стей и деформирования на ведущий и ведомый не представляется возможным. Эти процессы взаимосвязаны и влияют друг на друга.
Следовательно, экспериментальное исследование процесса ползучести на третьем участке кривой имеет существенное значение для исследования закономерностей процесса накопления повреждений, хотя в реальных циклических процессах деформирования материала с выдержками различной длительности разрушение (образование макротрещины) вследствие накопления повреждений происходит при малых деформациях ползучести (но при больших её траекториях). Поскольку для реальных условий эксплуатации конструкционных материалов начало установившегося участка достигается за несколько десятков часов, а общее время до разрушения составляет от нескольких сотен до тысяч часов, возникает проблема проведения испытаний на достаточно большой временной базе. Для её решения необходимо установить соответствие между испытаниями на кратковременную ползучесть при повышенных температурах и характеристиками ползучести и длительной прочности для большой временной базы при более умеренных температурах. Существует несколько различных зависимостей, устанавливающих это соответствие [6-10].
Большинство исследований конструкционных материалов на ползучесть и длительную прочность проведено при испытаниях на растяжение цилиндрических образцов при постоянных напряжениях и температурах. Однако реальные задачи ползучести и длительной прочности элементов конструкций связаны с гораздо более сложными условиями, например с переменностью напряжений, с изменением вида напряжённого состояния и циклически меняющимися температурами. Поскольку ползучесть металлов и длительная прочность весьма чувствительны к изменению напряжённого состояния и температуры, которая влияет на все их параметры, возникает проблема исследования закономерностей неизотермической ползучести и длительной прочности при различных видах напряжённых состояний и разработки методов оценки ползучести и длительной прочности материала конструктивных элементов в условиях многоосного напряжённого состояния на основании экспериментальных данных, полученных в результате испытаний при одноосном напряжённом состоянии.
Кривые долговечности, как правило, построены для постоянных напряжений (нагрузок), в то время как на практике напряжения и температуры меняются по довольно сложным законам. Вследствие этого возникает проблема суммирования повреждений для весьма нерегулярных историй нагружения. При оценке времени разрушения в условиях переменных напряжений и температур основным вопросом является возможность проведения расчётов по результатам испытаний при постоянных напряжениях и температурах. В технических приложениях в основном для этих целей используют правило линейного суммирования повреждений, согласно которому принимается, что все части использованной долговечности независимы друг от друга. Разрушение происходит тогда, когда сумма частей долговечности равна единице. Этот закон широко используется в инженерных расчётах для случаев работы конструкций в условиях переменных температур и переменных нагру-жений. Попытки проверить точность расчётов отдельно для изменений температур и напряжений очень малочисленны, хотя даже беглое рассмотрение кинетики процессов повреждения позволяет увидеть очевидное различие. Существует два основных недостатка имеющихся экспериментальных программ. Первый недостаток заключается в том, что к повреждению при ползучести относятся скорее как к теоретическому понятию, чем как к реальному физическому явлению, поддающемуся измерению. Второй недостаток связан с тем, что эксперименты планируются с целью проверки лишь какого-либо конкретного закона суммирования повреждений
и не ставят себе задачу выявления общих закономерностей этого процесса. В результате для каждой серии испытаний (изменений температур или напряжений) можно определить только, является ли закон суммирования пригодным, но нельзя установить причину непригодности. Анализ имеющихся экспериментальных результатов позволяет установить значительное влияние истории изменения напряжения и температуры на длительную прочность и сделать вывод о неприменимости правила линейного суммирования повреждений в условиях изменяющихся напряжений и температур. Была установлена хорошая корреляция между твердостью при комнатной температуре подверженного ползучести образца и длительной прочностью. Сравнение с линейным законом суммирования относительных долго-вечностей приводит к следующим отклонениям [11]:
- накопление повреждений при ползучести в испытуемой сплаве в общем случае не является линейной функцией времени предварительного воздействия;
- довольно сложно сравнить состояние повреждений для одинаковых по времени частей долговечности, так как отдельное влияние каждого из факторов (температуры или напряжения) на процесс зарождения и роста пор не изучено. Однако установлено, что для каждой заданной части долговечности суммарный объем пор возрастает с уменьшением напряжения (при постоянной температуре) и с уменьшением температуры (при постоянном напряжении). Для одного и того же воздействия в течение части долговечности степень размягчения (объем пор) возрастает с уменьшением напряжения, но почти не чувствительна к температуре;
- не существует какого-либо единого значения повреждения при ползучести, когда разрушение наступает всегда. Каждая совокупность испытаний или условий эксплуатации дает свое характерное значение повреждения при разрушении.
В настоящей работе с позиции механики повреждённой среды развита математическая модель МПС, описывающая процессы накопления повреждений в конструкционных материалах (металлах и их сплавах) при механизме деградации, обусловленном развивающимися деформациями ползучести. Проведены численные исследования процессов длительной прочности жаропрочного сплава ВЖ-159 и выполнено сравнение полученных численных результатов с данными натурных экспериментов.
1. Определяющие соотношения механики повреждённой среды
Основные положения рассматриваемого варианта МПС заключаются в следующем:
- материал среды изначально изотропен и повреждения в нём отсутствуют (учитывается только анизотропия, вызванная процессами деформирования; анизотропия упругих свойств, вызванная процессами повреждённости материала, не учитывается);
- компоненты тензоров деформаций в^ и скоростей деформаций вц представляют сумму «мгновенной» и «временной» составляющих. «Мгновенная» составляющая состоит из упругих компонент ве, в- , не зависящих от истории деформирования и определяющих конечное состояние процесса, и пластических компонент вРрз, в®Т, зависящих от истории процесса деформирования. Временная составляющая (деформаций ползучести вц , в ^) описывает временную зависимость процессов деформирования при низких скоростях нагружения;
- эволюция эквипотенциальных поверхностей ползучести описывается изменением её радиуса Сс и перемещением её центра рс;
- изменение объема элемента тела упруго, то есть ври = вС =0;
- рассматриваются процессы деформирования, характеризующиеся малыми деформациями;
- единственным структурным параметром, характеризующим на макроуровне степень повреждённости материала, является скалярный параметр ш - «повре-ждённость» (шо < ш < ш^ );
- влияние уровня накопленной повреждённости на процесс деформирования материала учитывается путём введения эффективных напряжений ас, .
Модель повреждённой среды состоит из трёх взаимосвязанных частей:
1) соотношений, определяющих вязкопластическое поведение материала с учётом зависимости от процесса разрушения;
2) уравнений, описывающих кинетику накопления повреждений;
3) критерия прочности повреждённого материала.
1.1. Соотношения термоползучести. Связь между тензором напряжений и тензором упругих деформаций на базе уравнений термоупругости
а = 3К [е - аТ] ,а = 3К
е — а Т — аТ
К
+ ка
с
а, = 20е£, а, = 2Се, + - а
е'е = еес-ег, ег, ег, ,
(1)
где а , е - шаровые, а а . , , е ., - девиаторные компоненты тензоров напряжений а, и деформаций е, соответственно; 0(Т) - модуль сдвига, К(Т) - объёмный модуль упругости, а(Т) - коэффициент температурного расширения - функция температуры Т .
Для описания процессов ползучести введем в пространстве напряжений семейство эквипотенциальных поверхностей ползучести Рс (поверхностей, на которых скорость деформации ползучести постоянна), имеющих общий центр р, и различные радиусы Сс, определяемые текущим напряжённым состоянием [12, 13]:
= Б,Б, - С2 =0, Б, = а, - рС,, г = 0,1, 2,...
Согласно ассоциированному закону
л«
дР(
ес = Л дРс = \ Бс
е, ~ Лс дБс. ~ Л,' (с)
(2)
(3)
где Лс соответствует текущей поверхности ), определяющей текущее напряжённое состояние Б,.
с,
Среди этих эквипотенциальных поверхностей можно выделить поверхность с радиусом Сс, соответствующую нулевой скорости ползучести
Р(0) = БС,Б, - Сс =0, Б, = а, - рс,,
(4)
где Б, и а, - совокупность напряжённых состояний, отвечающих (с определённым допуском) нулевой скорости ползучести.
Принимаем, что ,
С с = С с (Хс, Т), X с = \/ 3 ес, ес,, хс = I хс
(5)
Лс = Лс (фс,Т) , фс
Лс =
\/ЦЩ - С с
Сс
0, фс < 0;
Л фс > 0,
где Сс и Лс - экспериментально определяемые функции температуры Т.
Рис. 1. Зависимость \с от времени процесса t при SU = const при многоосном деформировании по лучевой траектории
Эволюционное уравнение для изменения координат центра поверхности ползучести есть [2]
(6)
Pi
С • С
9i eij
с с
92 Pij XX i,
где дС и д2 > 0 - экспериментально определяемые материальные параметры.
Конкретизируя соотношение (3), закон градиентальности представим в виде [14]
ej = Xi (фс,Т), S, откуда для XXс имеем [16]
Xi фс,
qc
Sij
Xc
VSij Sij— C i
Ci
i
Sij,
(7)
XX i
3
:rXi
Qi Qi Г*
Sij Sij — C i
(8)
Зависимость xc от t при SU = const при многоосном деформировании по лучевой траектории имеет вид, показанный на рис. 1.
На кривой Xc(t) (рис. 1) с определённой долей условности можно выделить три участка:
(I) участок неустановившейся ползучести (0 — xI"^] , на котором скорость деформации ползучести XXс убывает;
(II) участок установившейся ползучести (xi1 — xi2)), на котором скорость деформации ползучести Xс практически постоянна (Xс = const);
(III) участок неустановившейся ползучести [Xc2),X(3) , на котором деформации ползучести быстро растут (предшествует разрушению) и Xс резко возрастает.
Длины участков существенным образом зависят от величины SU = const.
Из (8) для трёх участков кривой ползучести выражение для Xc примет вид [13]:
Xc=
0,
Xi, X* XF,
Фс < ci или Xi = 0; 0 < Xi < x(1); x(1) < Xi < x(2);
X(2) < Xi < X(3),
(9)
где
Xi = x(0) 1 — + X(1) ei
i(1) "11 ,
ei(1): -11
XI I =3 Xi =2
¿Уст eii
aii— 2 Pii— a<
xII I = X"
Xi (1 — -У
3
получают из экспериментов при одноосном напряжённом состоянии лабораторного образца [13].
В формулах (9) введены обозначения: аС°') и - значение Ас в начальной
с(1) с(2)
и конечной точках первого участка кривой ползучести материала; е^ , е^ и
с(3)
е^ - границы участков кривой ползучести при одноосном напряжённом состо-
•с(нач) п
янии; ^11 - скорость деформации ползучести в начальный момент времени, <е1(1уст) - скорость деформации ползучести на участке установившейся ползучести (2-й участок кривой ползучести); ш - величина повреждённости материала; ас = у/2/3 Сс - предел ползучести при одноосном напряжённом состоянии; гс -параметр материала [12, 13].
Уравнения (1)-(9) описывают неустановившиеся и установившиеся участки кривой ползучести при разных уровнях напряжений и основные эффекты процесса ползучести при знакопеременном напряжении.
На стадии развития рассеянных по объему повреждений наблюдается влияние повреждённости на физико-механические характеристики материала. В первом приближении это влияние может быть учтено на базе концепции деградирующего континуума (введения эффективных напряжений) [2]
С < ~ г ( ) К
а'з= Р1(ш)аГ = ^ , 5 = Е2(ш)а = а, (10)
где С, К - эффективные модули упругости, определяемые по формулам Маккензи
£ = ^ - ш) (1 - "КМ*)- <11)
К = 4СК шгшш-- (12)
Аналогичным образом определяется эффективная внутренняя переменная рЗ :
С
рЗ = ¥1(ш)р% = с рЗ. (13)
С
1.2. Эволюционные уравнения накопления повреждений. Экспериментальный и теоретический анализ процессов повреждённости материала позволяет представить эволюционное уравнение накопления повреждений в элементарном объёме материалов в следующем общем виде [2, 4, 5]:
ш = /1 (0) /2 (в) /3 (ш) /4 (Й) (Й> (14)
где функция /1(0) описывает влияние кривизны траектории деформирования, /2 (в) - вид (объёмность) напряжённого состояния, /з(ш) - уровень накопленной повреждённости, /4 (й) - накопленную относительную энергию, затраченную на образование дефектов.
Учитывая тот факт, что к настоящему времени не имеется достаточно надежных систематизированных экспериментальных данных, характеризующих ползучесть материалов вплоть до разрушения в необходимом диапазоне рабочих нагрузок и температур и значительный разброс экспериментальных данных, эволюционное уравнение при ползучести необходимо формулировать в наиболее «простой» форме [5, 14-17]
ш = -с+г /с (в) (1 - «г* , (15)
' С + 1
где
Жс - Ж а
= , (16)
{¿с) = ^, Жс = ^ в3, (17)
Ж
Ж I
¡с (в) = ехр(ксв), (18)
где ш - значение величины повреждённости при ползучести; - значение энергии разрушения в конце стадии зарождения рассеянных повреждений при ползучести, а Ж/ - значение энергии, соответствующее образованию макроскопической
трещины; ¡с(в) - функция параметра объёмности напряжённого состояния в =
г
= а/аи; Шс = J Жс & - энергии, идущие на образование рассеянных поврежде-
0
ний при ползучести; кс, ас, гс - материальные параметры, зависящие от температуры Т .
1.3. Критерий прочности повреждённого материала. В качестве критерия окончания фазы развития рассеянных микроповреждений (стадии образования макротрещины) принимается условие достижения величины повреждённости своего критического значения
ш = ш} < 1. (19)
Интегрируя эволюционное уравнение накопления повреждений (15)-(18) с учётом определяющих соотношений термовязкопластичности (1)-(9) и критерия разрушения (19) по известной истории термомеханического нагружения в данном элементарном объёме материала, можно определить момент образования макроскопической трещины при ползучести.
Определение основных характеристик процесса вязкопластического деформирования повреждённых материалов (параметров состояния), которые в общем случае описываются тензорами а^ , в3, вР , рР, в3 и скалярами х, Ср, С с, Т и ш, может осуществляется при соответствующей формулировке определяющих соотношений МПС и линеаризации алгоритма нахождения Л, и сводится к формулировке определяющих соотношений МПС в приращениях, которые зависят от выбранного шага Дг. Шаг по времени Дг может корректироваться при прохождении сложных участков траектории деформирования в течении всего расчётного времени при условии устойчивости вычислений. Такой подход [18] наиболее удобен при решении краевых задач механики деформируемого твердого тела и используется в настоящей работе.
В общем случае напряжения, пластические деформации и деформации ползучести определяются интегрированием уравнений термоползучести (1)—(13) четырёхточечным методом Рунге-Кутты с коррекцией девиатора напряжений и последующим определением напряжений согласно уравнений термопластичности [5] с учётом средней скорости деформации ползучести в момент времени: гп+1 = гп + Дг.
2. Методика определения параметров определяющих соотношений МПС
Для практического применения уравнений термоползучести (1)—(18) необходимо иметь следующую информацию:
- зависимости 0(Т), К(Т), а(Т), где Т - температура;
- зависимость текущего радиуса поверхности ползучести, соответствующего нулевой скорости ползучести С c = С c (хc, T);
- зависимость A c = Ac (T) для различных участков кривой ползучести;
- зависимости модулей кинематического упрочнения gf(T), g2(T).
Материальные параметры эволюционных уравнений термоползучести определяются из базовых экспериментов [2].
В качестве основных базовых экспериментов принимаются эксперименты на одноосное растяжение-сжатие цилиндрических лабораторных образцов.
Основные типы базовых экспериментов - изотермические при постоянных базовых температурах Tj (j = 1, 2,...).
Типы образцов - цилиндрический сплошной и цилиндрический трубчатый. Выбранные типы образцов обеспечивают однородное распределение полей напряжений деформаций и температур в пределах рабочей части, исключают возможность потери устойчивости и формоизменении образца при знакопеременном на-гружении, максимально исключают влияние концентраторов на напряжённо-деформированное состояние при переходе от рабочей части образца к утолщённым местам [2, 12, 13].
Для определения модулей кинематического (анизотропного) упрочнения gf(T) и g2 (T) и зависимости для радиуса поверхности ползучести, соответствующего нулевой скорости ползучести С c = С c (х c,T) образец нагревается до значения температуры «базового» эксперимента T = Tj = const и проводятся испытания на ползучесть при одноосном напряжённом состоянии по схеме «мягкого» нагру-жения.
В начале образец нагружается до величины напряжения ст(1 в точке 1 (рис. 2, а). Этот уровень напряжений выбирается из анализа имеющегося веера кривых ползучести, полученных при «базовой» температуре T = Tj (кривая ползучести, соответствующая нулевой скорости ползучести). В результате релаксации процесс заканчивается в точке 2 (напряжение ст^), где скорость деформации ползучести близка к нулю.
Далее образец нагружается до напряжения обратного знака ст^ (точка 3 на рис. 2, а), и в результате процесса релаксации оказывается в точке 4. Таким образом, напряжения С(0)+ (точка 2) и С(0) (точка 4) характеризуют (с определённым допуском на остаточную деформацию) начальную верхнюю и нижнюю границы поверхности ползучести, соответствующие нулевой скорости ползучести.
Для определения трансформации поверхности ползучести на том же самом образце при заданном напряжении ст * 1 = const проводится ряд аналогичных действий после достижений назначенных уровней деформаций ползучести е^/ ),
ec(2) c(m)
е 1 1 , . . . , е 1 1 .
Полученный таким образом набор точек 2, 7, 12 17 и т. д. характеризует изменение верхней (при растяжении) границы поверхности ползучести в зависимости от накопленной деформации ползучести.
Точки 4, 8, 13, 19 и т. д. соответственно характеризуют изменение нижней (при сжатии) границы поверхности ползучести.
Таким образом, по результатам эксперимента при базовых постоянных температурах T = Tj найдены
- геометрическое место пределов ползучести при растяжении с заданным допуском на остаточную деформацию;
- геометрическое место обратных пределов ползучести при сжатии.
Рис. 2. Испытания на ползучесть при одноосном напряжённом состоянии по схеме «мягкого» (контролируемое напряжение) нагружения: а) экспериментальная информация для радиуса поверхности ползучести, соответствующей нулевой скорости деформации; Ь) экспериментальная информация для нахождения модулей кинематического (анизотропного) упрочнения
Зависимость радиуса поверхности ползучести соответствующей нулевой скорости деформации ползучести определяется по формуле
^ = J2_ a m)+ + * [?)~
m = 0,1, 2 .
Для определения модулей кинематического (анизотропного) упрочнения gC(T) и gC(T) необходимо проинтегрировать соотношение (6) при T = Tj = const. Имеем
Pi 1
^ (1 - e-g2e1i
gl
где e - основание натуральных логарифмов; g\ - тангенс угла наклона касательной к кривой p 11 ~ e 1 1 в начале координат (рис. 2.); р^ах = - предельное асимп-
тотическое значение ph при данной температуре T = Tj. Отсюда определяются модули gl и g2 анизотропного (кинематического) упрочнения.
Экспериментальное определение материальных параметров эволюционных уравнений накопления повреждений производится на второй стадии процесса накопления повреждений, с которой начинается влияние повреждённости на физико-механические характеристики материала, при одновременном расчёте экспериментальных процессов деформирования на этой стадии с использованием соотношений термовязкопластичности. Метод заключается в том, что все отклонения результатов численного моделирования процессов деформирования без учёта влияния повреждённости от экспериментальных на второй стадии процесса накопления повреждений приписываются влиянию повреждённости ш.
Приближенно границы W a, W f могут быть определены из испытаний на усталость при заданной амплитуде напряжений (или деформаций) на момент начала разупрочнения материала.
Для определения параметров эволюционного уравнения (15)-(18) при ползучести используется третий участок кривых ползучести efi(<7ii,Tj) при различных постоянных напряжениях и температурах. При этом значительное упрощение достигается, если имеет место подобие кривых ползучести [2]. В этом случае в качестве базовой кривой выбирают относительную кривую при T = Tj. Известное от-
c(3)
ношение скорости ползучести на третьем участке е^ к скорости установившейся ползучести еЦ^ позволяет определить параметр rc (ш, T) как функцию ш и T.
3. Результаты исследований
В [19] приведены результаты экспериментальных исследований процессов кратковременной высокотемпературной ползучести жаропрочного сплава ВЖ-159, которые были получены в лаборатории испытаний физико-механических свойств конструкционных материалов Научно-исследовательского института механики Нижегородского государственного университета им. Н.И. Лобачевского под общим руководством Д.А. Казакова. Эксперименты выполнены на универсальной испытательной машине ZWICKZ030 (Германия), с нагревательным устройством MAYTEC-HT080/1. Измерительная аппаратура включает в себя: датчик силы класса 1 по ISO 7500-1 (отклонение от задаваемого усилия в пределах одного опыта не превышает 15Н) и высокотемпературный датчик продольной деформации (тип -РМА-12/У7-1). Для контроля температуры внутри высокотемпературной печи по трём зонам и на образце используются три термопары (тип К). Поддерживание заданной температуры (отклонение не превышает 2 °C) осуществляется электронным блоком управления НТО-08/1.
Опыты были выполнены на цилиндрических образцах сплошного сечения с длиной рабочей части l = 50 мм и диаметром d = 8 мм в условиях одноосного растяжения для различных уровней нормальных напряжений стц и температур T.
Расчёт процесса кратковременной ползучести до образования макротрещины по механизму длительной прочности проводился с помощью программы EXPMODEL, предназначенной для численного моделирования неизотермического вязкопласти-ческого деформирования и накопления повреждений в конструкционных материалах (металлах и их сплавах) при нерегулярном нестационарном термомеханическом нагружении [20]. Физико-механические характеристики жаропрочного сплава ВЖ-159 и материальные параметры модели МПС следующие:
при T = 750 °C
K = 137000 МПа, G = 64000 МПа, Cc = 265 МПа, A^ = 0.0005 (МПа-ч)"1, A^ = 0.0005 (МПа-ч)"1, g<c = 1100 МПа, = 150 МПа, Wf = 57.5 МДж/м3, W = 14.3 МДж/м3, ac = 1, rc = 0.1, кс = 1, х^ =0, х^ = 0.025, = 0.83;
при T = 850 °C
K = 113000 МПа, G = 52000 МПа, Cc = 85 МПа, A^ = 0.00068 (МПа-ч)"1, A^ = 0.00068 (МПа-ч)"1, g1 = 1000 МПа, g2 = 150 МПа, Wf = 14 МДж/м3, W^ = 5 МДж/м3, ac = 1, rc =0.1, кс = 1, х(1) =0, х(2) = 0.06, = 0.83.
На рис. 3-6 представлены кривые ползучести. Здесь сплошными линиями отмечены результаты численного моделирования с использованием определяющих соотношений МПС (1)—(19), а пунктирной - соответствующие экспериментальные результаты. Видно качественное и количественное совпадение опытных и расчётных данных.
На рис. 7-11 приведены зависимости, полученные из эксперимента, проведенного при режиме нагрузки, характеристики которого представлены на рис. 6.
<*1
Т=750 С
0~11 =450 МПа
"11 =350 МПа
¡Г'
/
/
«11
1 1
=850 С
/ (Гц =150 МПа 1
и
1
°11=120 МПа
зо ¡.час
I,час
Рис. 3. Кривые ползучести при температуре Т = 750 °С и напряжениях ац = = 350 МПа и ап = 450 МПа
Рис. 4. Кривые ползучести при температуре Т = 850 °С и напряжениях ац = = 120 МПа и ац = 150 МПа
«1]
Т=850 С 1
СТц =150 МПа \ II
СТ, ] = 120 МПа -----
¡.час
Рис. 5. Кривые ползучести при температуре Т = 750 °С и переходе с уровня напряжений ац = 350 МПа на уровень ац = = 450 МПа
Рис. 6. Кривые ползучести при температуре Т = 850 °С и переходе с уровня напряжений ац = 120 МПа на уровень ац = = 150 МПа
В целом, анализируя полученные численные результаты, можно отметить, что предложенная модель повреждённой среды качественно и количественно описывает основные эффекты наблюдаемые при нестационарной ползучести конструкционных материалов (металлах и их сплавах) и деградацию начальных прочностных свойств материалов при деградации по механизму длительной прочности.
Заключение
Развита математическая модель механики повреждённой среды, описывающая процессы неупругого деформирования и накопления повреждений в конструкционных материалах (металлах и их сплавах) при деградации начальных прочностных свойств материалов по механизму длительной прочности.
Методом численного моделирования и сравнения полученных результатов с опытными данными проведена оценка достоверности определяющих соотношений МПС при ползучести, которая позволяет сделать вывод о достоверности развитых модельных представлений и точности разработанной методики определения материальных параметров, входящих в указанные соотношения.
Благодарности. Работа выполнена при финансовой поддержке Министерства науки и высшего образования Российской Федерации (Соглашение № 075-11-2019050).
Рис. 7. График зависимости эффективных напряжений 5ц от времени процесса Ь
Рис. 8. График зависимости pïi от времени процесса t
Рис. 9. График зависимости эффективного радиуса поверхности ползучести нулевого уровня Сс от времени процесса Ь
Рис. 10. График зависимости энергии разрушения при ползучести Шс от времени процесса Ь
Рис. 11. График зависимости величины повреждённости ш от времени процесса t
Литература
1. Митенков Ф.М., Камдалов В.Ф., Коротких Ю.Г. Методы обоснования ресурса ядерных энергетических установок. - М.: Машиностроение, 2007. - 448 с.
2. Волков И.А., Коротких Ю.Г. Уравнения состояния вязкоупругопластических сред с повреждениями. - М.: Физматлит, 2008. - 424 с.
3. Локощенко А.М. Ползучесть и длительная прочность металлов. - М.: Физматлит, 2016. - 504 с.
4. Lemaitre J. Damage modelling for prediction of plastic or creep fatigue failure in structures // Trans. 5th Int. Conf. SMiRT-5. - Amsterdam: North Holland, 1979. - V. L. -Paper L5/1 b, P. 1-8.
5. Murakami S. Mechanical description of creep damage and its experimental verification // J. Mec. Theor. Appl. - 1982. - V. 1, No 5. - P. 743-761.
6. Manson S.S., Ensign C.R. A quarter-century of progress in the development of correlation and extrapolation methods for creep rupture data // J. Eng. Mater. Technol. - 1979. -V. 101, No 4. - P. 317-325. - doi: 10.1115/1.3443696.
7. Le May I. Developments in parametric methods for handling creep and creep-rupture data // J. Eng. Mater. Technol. - 1979. - V. 101, No 4. - P. 326-330. - doi: 10.1115/1.3443697.
8. Larson F.R., Miller J.A. A time-temperature relationship for rupture and creep stresses // Trans. ASME. - 1952. - V. 74, No 5. - P. 539-605.
9. Никитенко А.Ф. Экспериментальное обоснования гипотезы существования поверхности ползучести в условиях сложного нагружения. Сообщение 1 // Проблемы прочности. - 1984. - № 8. - С. 3-8.
10. Никитенко А.Ф. Экспериментальное обоснования гипотезы существования поверхности ползучести в условиях сложного нагружения. Сообщение 2 // Проблемы прочности. - 1984. - № 8. - С. 8-10.
11. Вудфорд Д.А. Повреждение при ползучести и концепция остаточной долговечности // Труды Амер. о-ва инж.-мех. Сер. Д. Теорет. основы инж. расчетов. - 1979. -Т. 101, № 4. - С. 1-8.
12. Волков И.А., Игумнов Л.А., Коротких Ю.Г. Прикладная теория вязкопластич-ности. - Н. Новгород: Изд-во ННГУ, 2015. - 318 с.
13. Волков И.А., Игумнов Л.А., Казаков Д.А., Шишулин Д.Н., Сметанин И.В. Определяющие соотношения нестационарной ползучести при сложном напряженном состоянии // Проблемы прочности и пластичности. - 2016. - Т. 78, № 4. - С. 436-451. -doi: 10.32326/1814-9146-2016-78-4-436-451.
14. Боднер С.Р., Линдхолм У.С. Критерий приращения повреждения для зависящего от времени разрушения материалов // Труды Амер. о-ва инж.-мех. Сер. Д. Теорет. основы инж. расчетов. - 1976. - Т. 100, № 2. - С. 51-58.
15. Леметр Ж. Континуальная модель повреждения, используемая для расчёта разрушения пластичных материалов // Труды Амер. о-ва инж.-мех. Сер. Д. Теорет. основы инж. расчетов. - 1985. - Т. 107, № 1. - C. 90-98.
16. Пежина П. Моделирование закритического поведения и разрушения диссипативного твердого тела // Труды Амер. о-ва инж.-мех. Сер. Д. Теорет. основы инж. расчетов. -
1984. - Т. 106, № 4. - С. 107-117.
17. Локощенко А.М. К выбору критерия длительной прочности при сложном напряженном состоянии // Проблемы прочности. - 1989. - № 9. - С. 3-6.
18. Бантхья B., Мукерджи C. Об усовершенствованной схеме интегрирования по времени для системы определяющих соотношений неупругой деформации с нелинейностью жёсткого типа // Труды Амер. о-ва инж.-мех. Сер. Д. Теорет. основы инж. расчетов. -
1985. - № 4. - С. 54-60.
19. Капустин С.А., Казаков Д.А., Чурилов Ю.А., Галущенко А.И., Вахтеров А.М. Экспериментально-теоретическое изучение поведения изделий из жаропрочного сплава в условиях высокотемпературной ползучести // Проблемы прочности и пластичности. - 2008. - Вып. 70. - С. 100-111.
20. Волков И.А., Игумнов Л.А., Казаков Д.А., Емельянов А.А., Тарасов И.С., Гусева М.А. Программная реализация процессов вязкопластического деформирования и накопления повреждений в конструкционных сплавах при термомеханическом на-гружении // Проблемы прочности и пластичности. - 2016. - Т. 78, № 2. - С. 188-207. -doi: 10.32326/1814-9146-2016-78-2-188-207.
Поступила в редакцию 01.10.19
Волков Иван Андреевич, доктор физико-математических наук, главный научный сотрудник
Научно-исследовательский институт механики Национального исследовательского Нижегородского государственного университета им. Н.И. Лобачевского
пр. Гагарина, д. 23, к. 6, г. Нижний Новгород, 603950, Россия E-mail: [email protected]
Игумнов Леонид Александрович, доктор физико-математических наук, профессор, директор
Научно-исследовательский институт механики Национального исследовательского Нижегородского государственного университета им. Н.И. Лобачевского
пр. Гагарина, д. 23, к. 6, г. Нижний Новгород, 603950, Россия E-mail: [email protected]
Литвинчук Светлана Юрьевна, кандидат физико-математических наук, ведущий научный сотрудник
Научно-исследовательский институт механики Национального исследовательского Нижегородского государственного университета им. Н.И. Лобачевского
пр. Гагарина, д. 23, к. 6, г. Нижний Новгород, 603950, Россия E-mail: [email protected]
ISSN 2541-7746 (Print) ISSN 2500-2198 (Online) UCHENYE ZAPISKI KAZANSKOGO UNIVERSITETA. SERIYA FIZIKO-MATEMATICHESKIE NAUKI (Proceedings of Kazan University. Physics and Mathematics Series)
2019, vol. 161, no. 4, pp. 509-525
doi: 10.26907/2541-7746.2019.4.509-525
A Continual Model of Damaged Media for Describing the Creep Failure Process
I.V. Volkov* , L.A. Igumnov**, S.Yu. Litvinchuk***
Research Institute of Mechanics, Lobachevsky National Research State University, Nizhny Novgorod, 603950 Russia E-mail: *[email protected], **[email protected], ***[email protected]
Received October 1, 2019 Abstract
The main laws of deformation and degradation processes of the initial strength properties of structural materials (metals and their alloys) according to the long-term strength mechanism were considered. From the viewpoint of the mechanics of damaged media (MDM), a mathematical model was created to describe inelastic deformation and damage accumulation during the creep process. An experimental and theoretical methodology for determining material parameters of the derived defining relations of MDM was developed. The material parameters and scalar functions of the MDM model, its reliability and the scope of its applicability were determined based on the analysis of the processes of deformation and failure of laboratory specimens in the conditions of soft loading (controlled stresses). The results of the experimental study of short-term creep of the VZh-159 heat-resistant alloy were presented. The process of deformation and damage accumulation was numerically analyzed; the obtained numerical results were compared with the data of the full-scale experiments. Comparison of the numerical and experimental data shows that the introduced defining relations of MDM
adequately describe the response of materials in the conditions of degradation of the initial strength properties of structural materials according to the long-term strength mechanism.
Keywords: nonstationary creep, long-term strength, simulation, defining relations, mechanics of damaged media, temperature, damage, material parameter, numerical experiment, full-scale experiment
Acknowledgments. This work was supported by the Ministry of Science and Higher Education of the Russian Federation (project no. 075-11-2019-050).
Figure Captions
Fig. 1. Dependences of \c on process time t at SU = const under multiaxial deformation along the ray path.
Fig. 2. Uniaxial stress state creep tests under soft (controlled stress) loading: (a) experimental data for the creep surface radius of zero strain rate; (b) experimental data for finding the kinematic (anisotropic) hardening modules.
Fig. 3. Creep curves at temperature T = 750 °C and stresses an = 350 MPa and on = 450 MPa.
Fig. 4. Creep curves at temperature T = 850 °C and stresses an = 120 MPa and an = 150 MPa.
Fig. 5. Creep curves at temperature T = 750 °C and transition from stress level an = = 350 MPa to level an = 450 MPa.
Fig. 6. Creep curves at temperature T = 850 °C and transition from stress level an = = 120 MPa to level an = 150 MPa.
Fig. 7. Graph of effective stresses on versus process time t.
Fig. 8. Graph of pn versus process time t.
Fig. 9. Graph of effective radius of zero level creep surface Cc versus process time t.
Fig. 10. Graph of creep failure energy Wc versus process time t.
Fig. 11. Graph of damage degree w versus process time t.
References
1. Mitenkov F.M., Kaidalov V.F., Korotkikh Yu.G. Metody obosnovaniya resursa yadernykh energeticheskikh ustanovok [Methods of Substantiation of Nuclear Power Plants' Life]. Moscow, Mashinostroenie, 2007. 448 p. (In Russian)
2. Volkov I.A., Korotkikh Yu.G. Uravneniya sostoyaniya vyazkouprugoplasticheskikh sred s povrezhdeniyami [Equations of the State for the Viscoelasticoplastic Media with Defects]. Moscow, Fizmatlit, 2008. 424 p. (In Russian)
3. Lokoshchenko A.M. Polzuchest' i dlitel'naya prochnost' metallov [Creep and Long-Term Strength of Metals]. Moscow, Fizmatlit, 2016. 504 p. (In Russian)
4. Lemaitre J. Damage modelling for prediction of plastic or creep fatigue failure in structures. Trans. 5th Int. Conf. SMiRT-5. Amsterdam, North Holland, 1979, vol. L, art. L5/1 b, pp. 1-8.
5. Murakami S. Mechanical description of creep damage and its experimental verification. J. Mec. Theor. Appl., 1982, vol. 1, no. 5, pp. 743-761.
6. Manson S.S., Ensign C.R. A quarter-century of progress in the development of correlation and extrapolation methods for creep rupture data. J. Eng. Mater. Technol., 1979, vol. 101, no. 4, pp. 317-325. doi: 10.1115/1.3443696.
7. Le May I. Developments in parametric methods for handling creep and creep-rupture data. J. Eng. Mater. Technol., 1979, vol. 101, no. 4, pp. 326-330. doi: 10.1115/1.3443697.
8. Larson F.R., Miller J.A. A time-temperature relationship for rupture and creep stresses. Trans. ASME, 1952, vol. 74, no. 5, pp. 539-605.
9. Nikitenko A.F. Experimental verification of the hypothesis of the existence of surface creep under complex loading conditions. Report No. 1. Strength Mater., 1984, vol. 16, no. 8, pp. 1063-1068. doi: 10.1007/BF01530274.
10. Nikitenko A.F. Experimental verification of the hypothesis of the existence of surface creep under complex loading conditions. Report No. 2. Strength Mater., 1984, vol. 16, no. 8, pp. 1069-1071. doi: 10.1007/BF01530275.
11. Woodford D.A. Creep damage and the remaining life concept. J. Eng. Mater. Technol., 1979, vol. 101, no. 4, pp. 311-316. doi: 10.1115/1.3443695.
12. Volkov I.A., Igumnov L.A., Korotkikh Yu.G. Prikladnaya teoriya vyazkoplastichnosti [Applied Theory of Viscoplasticity]. Nizhny Novgorod, Izd. NNGU, 2015. 318 p. (In Russian)
13. Volkov I.A., Igumnov L.A., Kazakov D.A., Shishulin D.N., Smetanin I.V. Constitutive relations of unsteady creep in a complex stress state. Probl. Prochn. Plast., 2016, vol. 78, no. 4, pp. 436-451. doi: 10.32326/1814-9146-2016-78-4-436-451. (In Russian)
14. Bodner S.R., Lindholm U.S. An incremental criterion for time-dependent failure of materials. J. Eng. Mater. Technol., 1976, vol. 98, no. 2, pp. 140-145. doi: 10.1115/1.3443356.
15. Lemaitre J. A continuous damage mechanics model for ductile fracture. J. Eng. Mater. Technol., 1985, vol. 107, no. 1, pp. 83-89. doi: 10.1115/1.3225775.
16. Perzyna P. Constitutive modeling of dissipative solids for postcritical behavior and fracture. J. Eng. Mater. Technol., 1984, vol. 106, no. 4, pp. 410-419. doi: 10.1115/1.3225739.
17. Lokoshchenko A.M. Criteria for determining the long-term strength under conditions of complex loading. Strength Mater., 1989, vol. 21, no. 9, pp. 1121-1124. doi: 10.1007/BF01529282.
18. Banthia V., Mukherjee S. On an improved time integration scheme for stiff constitutive models of inelastic deformation. J. Eng. Mater. Technol., 1985, vol. 107, no. 4, pp. 282285. doi: 10.1115/1.3225820.
19. Kapustin S.A., Kazakov D.A., Churilov Yu.A., Galushchenko A.I., Vakhterov A.M. Experimental and theoretical study of the behavior of structural parts from heat-resistant alloy under high-temperature creep. Probl. Prochn. Plast., 2008, no. 70, pp. 100-111. (In Russian)
20. Volkov I.A., Igumnov L.A., Kazakov D.A., Emel'yanov A.A., Tarasov I.S., Guseva M.A. Software implementation of viscoplastic deformation and damage accumulation processes in structural alloys under thermal-mechanical loading. Probl. Prochn. Plast., 2016, vol. 78, no. 2, pp. 188-207. doi: 10.32326/1814-9146-2016-78-2-188-207. (In Russian)
Для цитирования: Волков И.А., Игумнов Л.А., Литвинчук С.Ю. Континуальная / модель повреждённой среды для описания процесса разрушения при ползучести // \ Учен. зап. Казан. ун-та. Сер. Физ.-матем. науки. - 2019. - Т. 161, кн. 4. - С. 509-525. -doi: 10.26907/2541-7746.2019.4.509-525.
For citation: Volkov I.V., Igumnov L.A., Litvinchuk S.Yu. A continual model of da-/ maged media for describing the creep failure process. Uchenye Zapiski Kazanskogo \ Universiteta. Seriya Fiziko-Matematicheskie Nauki, 2019, vol. 161, no. 4, pp. 509-525. doi: 10.26907/2541-7746.2019.4.509-525. (In Russian)