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

Радченко В.П., Саушкин М.Н., Бочкова Т.И. Математическое моделирование формирования и релаксации остаточных напряжений в плоских образцах из сплава ЭП742 после ультразвукового упрочнения в условиях высокотемпературной ползучести // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2016. - № 1. - С. 93-112. DOI: 10.15593/perm.mech/2016.1.07

Radchenko V.P., Saushkin M.N., Bochkova T.I. Mathematical modeling and experimental study of forming and relaxation of the residual stresses in plane samples made of EP742 alloy after the ultrasonic hardening under the high-temperature creep conditions. PNRPUMechanics Bulletin. 2016. No. 1. Pp. 93-112. DOI: 10.15593/perm.mech/2016.1.07

ВЕСТНИК ПНИПУ. МЕХАНИКА № 1,2016 PNRPU MECHANICS BULLETIN

http ://vestnik.pstu. ru/mechanics/about/inf/

Э01 10.15593/регш.шесЬ/2016.1.07 УДК 539.3 + 539.376

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ И ЭКСПЕРИМЕНТАЛЬНОЕ ИССЛЕДОВАНИЕ ФОРМИРОВАНИЯ И РЕЛАКСАЦИИ ОСТАТОЧНЫХ НАПРЯЖЕНИЙ В ПЛОСКИХ ОБРАЗЦАХ ИЗ СПЛАВА ЭП742 ПОСЛЕ УЛЬТРАЗВУКОВОГО УПРОЧНЕНИЯ В УСЛОВИЯХ ВЫСОКОТЕМПЕРАТУРНОЙ ПОЛЗУЧЕСТИ

В.П. Радченко, М.Н. Саушкин, Т.И. Бочкова

Самарский государственный технический университет, Самара, Россия

АННОТАЦИЯ

Выполнено комплексное расчетно-экспериментальное исследование остаточных напряжений в призматических образцах из сплава ЭП742 после ультразвукового упрочнения и температурной выдержки без нагрузки при температуре 650 °С в течение 100 часов. Выяснены закономерности распределения остаточных напряжений по толщине поверхностно-упрочненного слоя. Экспериментально установлено, что обработка ультразвуковым упрочнением образцов вызывает появление в поверхностном слое сжимающих остаточных напряжений, максимум которых наблюдается в подповерхностном слое, а спад - к поверхности. При температурной выдержке происходит релаксация наведенных сжимающих остаточных напряжений, уровень остаточных напряжений уменьшается в 1,4-1,6 раз и происходит смещение их максимума вглубь образца, однако толщина сжатого слоя сохраняется на уровне 200 мкм.

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

О СТАТЬЕ

Получена: 17 февраля 2016 г. Принята: 18 марта 2016 г. Опубликована: 30 марта 2016 г.

Ключевые слова:

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

© Радченко Владимир Павлович - доктор физико-математических наук, профессор, e-mail: radch@samgtu.ru. Саушкин Михаил Николаевич - кандидат физико-математических наук, доцент, e-mail: saushkin.mn@samgtu.ru. Бочкова Татьяна Игоревна - магистрант, e-mail: tanechka.bochkova@mail.ru.

Vladimir P. Radchenko - Doctor of Physical and Mathematical Sciences, Professor, e-mail: radch@samgtu.ru. Mikhail N. Saushkin - Ph.D. in Physical & Mathematical Sciences, Associate Professor, e-mail: saushkin.mn@samgtu.ru. Tat'iana I. Bochkova - Master Student, e-mail: tanechka.bochkova@mail.ru.

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

© ПНИПУ

A MATHEMATICAL MODELING AND EXPERIMENTAL STUDY OF FORMING AND RELAXATION OF THE RESIDUAL STRESSES IN PLANE SAMPLES MADE OF EP742 ALLOY AFTER THE ULTRASONIC HARDENING UNDER THE HIGH-TEMPERATURE CREEP CONDITIONS

V.P. Radchenko, M.N. Saushkin, T.I. Bochkova

Samara State Technical University, Samara, Russian Federation

ARTICLE INFO

ABSTRACT

Received: 17 February 2016 Accepted: 18 March 2016 Published: 30 March 2016

Keywords:

ultrasonic hardening, flat samples, residual stresses, mathematical model, high-temperature creep, residual stresses relaxation, experimental data

The authors carried out a complex (computational and experimental) research of residual stresses in prismatic samples made of EP742 alloy after ultrasonic hardening and he temperature of 650 °С for 100 hours unloaded. The laws of residual stresses distribution over the thickness of the surface-hardened layer were discovered. The experiments proved that the ultrasonic hardening of the sample caused compression residual stresses in the surface layer. The maximum stresses were observed in the subsurface layer and they decreased on reaching the surface. When the temperature was exposed, the induced compression residual stresses became relaxed, the level of residual stresses decreased by 1.4-1.6 times and its maximum was displaced deep into the sample, but the thickness of the compressed layer was kept to about 200 micron.

We developed the mathematical model of residual stresses forming in prismatic samples after surface plastic deformation and relaxation under high-temperature creep conditions. We used the hardened half-space as a model object for the prismatic sample, because the hardened layer was rather thin. We introduced Cartesian coordinate system, where x0y plane coincided with the half-space hardened surface, and 0z axis was directed into the depth of the hardened layer.

©PNRPU

Введение

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

и микротвердости поверхностных слоев при нормальных температурах отмечалась во многих работах как отечественных [1-6], так и зарубежных [7-12] исследователей.

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

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

Глубокая теоретическая проработка задачи для оценки напряженно-деформированного состояния в упрочненных цилиндрических деталях после процедуры упрочнения выполнена в работах [18, 19], а в последующих публикациях [20, 21] предложена математическая модель расчета остаточных напряжений с учетом эффектов деформационной анизотропии при поверхностном пластическом упрочнении. По всей видимости, отсутствие методов расчета начального напряженно-деформированного состояния в упрочненных конструкциях (после процедуры упрочнения) было основным сдерживающим фактором теоретического решения задачи релаксации наведенных остаточных напряжений в деталях в условиях высокотемпературной ползучести, хотя потребности, например, авиадвигателестроения и энергомашиностроения, стимулировали чисто экспериментальные исследования в этом направлении [22-28].

В публикации [29] впервые был представлен прямой численный метод решения краевой задачи кинетики остаточных напряжений в упрочненном цилиндрическом изделии в условиях ползучести, в работе [30] он был применен к задаче релаксации остаточных напряжений для упрочненного цилиндра при действии растягивающей нагрузки, а в [31] -при совместном действии растягивающей осевой нагрузки и крутящего момента. Однако теоретический подход решения задач такого рода для плоских образцов в настоящее время отсутствует. Это и является целью настоящей работы наряду с экспериментальной проверкой соответствующего метода решения краевой задачи для плоских образцов из сплава ЭП742, упрочненных УЗУ, в условиях высокотемпературной выдержки при температуре 650 °С и длительности 100 часов.

1. Методика расчета остаточных напряжений и пластических деформаций после процедуры упрочнения

Рассматривается призматический образец (прямоугольный параллелепипед), одна сторона которого подверглась процедуре поверхностно-пластического деформирования. Как правило, толщина упрочненного слоя составляет 100-200 мкм. Будем считать, что величины трех измерений призматического образца на 2-3 порядка больше толщины упрочненного слоя. Тогда в качестве модельной можно взять задачу об упрочненном полу-

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

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

При введенных ограничениях остаточные сжимающие напряжения и пластические деформации не зависят от координат х и у, а зависят лишь от переменной z. Тогда аг. =аг.qi = qi(^)(/ = х,у,z). В дальнейшем, где это возможно, для краткости будем

опускать в записи переменную z, при этом всегда будем считать, что 0 < z < +<».

Поскольку толщина упрочненного слоя составляет всего 100-200 мкм, его можно считать «наклеенным» на плоскость, ограничивающую полупространство. Тогда можно ввести гипотезу плоских сечений (x0z и y0z) для полных деформаций, а поскольку 0 < z < +<», математически это запишется следующим образом:

8х

где 8х(z) и 8 (z) - компоненты тензора полных деформаций.

Введем в рассмотрение гипотезу анизотропного упрочнения поверхности в виде

ях =щу,

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

Из условия пластической несжимаемости qx + qy + qz = 0 и (2) получаем

qz =-

Учитывая, что 8i (z) = е (z) + qi (z) (/ = х,у, z ), где е1 - компоненты тензора упругих деформаций, запишем соотношения (1) в развернутом виде:

Ех -Чау +аz )] + Ях = 0 Еу -У

где Е - модуль Юнга; V - коэффициент Пуассона. В этих соотношениях а2 = 0, что следует из условия равновесия

Эа _

&

= 0 (1хш аz (z) = 0).

В силу этого из (2) и (4) имеем

1+ av

у

a + v

Теперь из (2)-(5) нетрудно получить

q a

qX =--, qV =--, 4? =-.

X E(a + v) X y E(a + v) X z E(a + v) X

Таким образом, из формул (5), (6) следует, что если известна компонента тензора остаточных напряжений aX = aX (z), 0 < z <+оо, и коэффициент анизотропии a, то все

компоненты тензоров остаточных напряжений и пластических деформаций определяются через эти величины. Иными словами, чтобы определить напряженно-деформированное состояние в упрочненном полупространстве, достаточно знать экспериментальную диаграмму aX =aX (z) и величину a . В частных случаях, например при термопластическом упрочнении, упрочнении дробью, ультразвуковом упрочнении поверхности, величина a = 1 [20, 21] и формулы (5), (6) упрощаются:

1 -v 2

Поскольку в дальнейшем рассматривается процедура ультразвукового упрочнения, из (7) следует, что величина ay (z) и значения остаточных пластических деформаций

4 (z) (/ = х, у, z) после УЗУ определяется через зависимость aX = aX (z). Отсюда предлагается следующий феноменологический метод определения напряженно-деформированного состояния после упрочнения: величина ax =ax (z) определяется экспериментально (в настоящей работе на основе метода полосок [32]), а остальные компоненты находятся с использованием формул (7).

2. Методика идентификации параметров аппроксимации для компонент тензора остаточных напряжений aX = а х

В соответствии с модельным представлением упрочненного призматического образца упрочненным полупространством диаграмма aX = aX (z) с учетом естественного условия limz >|оaX(z) = 0 после процедуры поверхностного пластического упрочнения схематически должна иметь вид, представленный на рис. 1. При этом учитывается случай, когда максимум сжимающих остаточных напряжений находится как в подповерхностном слое ( z* ^ 0), так и на упрочненной поверхности ( z* = 0). Однако экспериментальный метод [32] позволяет определить значения aX =aX (z) лишь в тонком упрочненном слое 0 < z < z0, то есть в области сжатия материала (см. рис. 1). Но для решения задачи релаксации остаточных напряжений при ползучести необходимо иметь непрерывные поля остаточных напряжений и пластических деформаций во всей области интегрирования (во всем полупространстве). Поэтому первой задачей является аппроксимация экспериментальной эпюры aX = aX (z) для всех 0 < z < +о . Другими словами, необходимо выполнить

экстраполяцию с области сжатия 0 < z < z0, где определяются экспериментальные значения, на область z0 < z < , но таким образом, чтобы выполнялось условие самоуравновешенности остаточных напряжений

"Г*-"-'

|а х ( z № = 0,

поскольку к полупространству не приложены внешние силы.

Рис. 1. Схематическая экспериментальная

эпюра а х =ах

Выберем аппроксимацию а х = ах (z) в виде

а х (z) = а0ехр

' (z - z*)2 ^

а1 ехр

' (z - z*)2^

где а0, а1 и Ь - параметры, для определения которых (в первом приближении) можно использовать значения эпюры ах =ах (z) в характерных точках: а х (0) = а *, ах (z*) = аш;п, ах(z0) = 0 и условие (8); параметр I = 1 и имеет размерность длины, как и z . В результате с учетом (9) получаем следующую систему уравнений:

а0ехр(-z *2)-а1 ехр(-z *2 /Ь2) = а*,

а0 - а, = аш

а0 ехр(-(z0 - ^)2)- а1 ехр(-(z0 - ^)2 / Ъ2) = 0,

+<Х)

|а х (z )dz = 0.

Система уравнений (10) является переопределенной (3 параметра и 4 уравнения). Для устранения переопределенности можно отказаться от строгого выполнения одного из условий для эпюры ах = ах(г), например от ее прохождения через точку (z0,0). Тогда, от-

брасывая третье соотношение в (10) и подставляя (9) в (8), получим систему трех уравнений относительно трех неизвестных а0, а1 и Ь :

а0 ехр(-г *2) - а1 ехр(-г *2 /Ь2) = а*,

а0 -а1 = ашт ,

а0 | ехр(-(г - г*)2)й?г - а1| ехр(-(г - г*)2 / Ь2)й?г = 0,

о 0

2 х

которая решается численно с использованием функции егДх) = —ехр(-2)^.

у/п 0

Уточнение значений параметров а0, а1 и Ь (второе приближение) производится в результате минимизации отклонения расчетных данных ах (г1) от экспериментальных аХ (г-) в следующей норме:

Д =

1

Е"=х[а х

2Г.х[аХ

где гг. - точки дискретизации; п - количество точек дискретизации зависимости ах =а х (г). Для этой цели, начиная с первого приближения, варьируются величины а0, а1, а величина Ь определяется из третьего уравнения (11) (тем самым сохраняется самоуравновешенность эпюры ах =ах (г)) до достижения минимума функционалом (12). По

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

Таким образом, после идентификации параметров аппроксимации (9) величина ах = ах(г) будет иметь аналитическое представление для всех 0 < г < +оо, а значит, будут

иметь аналитические представления и остальные компоненты тензоров остаточных напряжений и пластических деформаций в соответствии с формулами (5), (6) при а^ 1 или (7) при а = 1.

3. Методика расчета кинетики напряжений в поверхностно-упрочненном полупространстве в условиях ползучести

Рассмотрим задачу расчета кинетики остаточных напряжений, наведенных в процессе поверхностного пластического деформирования полупространства, в условиях высокотемпературной ползучести материала при отсутствии внешних сил. Такой режим в дальнейшем будем называть термоэкспозицией (температурная выдержка без нагрузки). В качестве начальных условий краевой задачи используется напряженно-деформированное состояние, возникающее после процедуры упрочнения, которое сформировано в начальный момент времени / = 0 - 0 при температуре Тх (как правило, комнатная температура). Пусть при / = 0 происходит мгновенное прогревание полупространства с температуры Т1 до температуры Т2 (Т2 > 7Ц), при которой возникает процесс ползучести материала. Через

Е1 обозначим модуль Юнга при температуре Т1, а через Е2 - при температуре Т2. Предполагаем также, что новых пластических деформаций при температурном нагреве не возникает и они определяются равенствами (6) при Е = Е1. Температурные деформации в дальнейшем не учитываются, поскольку они дают просто равномерное увеличение объема и не приводят к изменению остаточных напряжений в силу равномерного прогревания тела. В дальнейшем будем рассматривать общий случай анизотропного упрочнения

В момент времени / = 0 + 0 при температуре Т2 в силу сохранения гипотезы плоских

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

( ) Е2(а + у) 1 + уа

ах (г) =--2-—ду (2X а у (2) =-ах (2).

1 -V2 а + у

Пусть теперь в течение времени / е [0,/*] образец выдерживается при температуре Т2,

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

в х (г,/) = Е- [а х (г) - Уа у (г,/)] + цх (г ) + рх (г) = 0,

2

(г,/) = Е- [а у (г,/) - уа х (г,/)] + цу (г ) + ру ( г,/) = 0.

Вуч " Е

2

Компонента деформаций вг (г, /) является «пассивной» и с учетом несжимаемости деформации ползучести рх + ру + р2 = 0 определяется из соотношения

в 2 (г,/) = - -у-[а х (г,/) + а у (г,/)] + д2 (г ) - ( (г) + ру (г,/)).

Е2

Решая систему уравнений (13) относительно ах (г, /) и а у (г,/), получаем соотношения, описывающие кинетику этих напряжений во времени вследствие ползучести:

Е г I М

а х (г,/) = [^ (г ) + рх ( г) + ( г ) + ру ( г))],

У Ё (14) а у (г) = [^у (г ) + ру (г,Х) + (г ) + рх ( г,Х ))].

Таким образом, если известны значения компонент тензора деформаций ползучести рх (г, /) и ру (г,/), то величины ах (г,/) и а у (г,/) определяются из (14).

Величины рх (г,/) и ру (г,/) вычисляются численно «шагами» по времени на основании выбранной теории ползучести, которая будет описана далее. Суть метода состоит в следующем. Пусть выполнена дискретизация по времени 0 = /0 <<••• </п = /* с шагом

= 1+х - (/ = 0,1,...,п-1) и нам известны значения ру (г,/) ( у = х,у, г ). Тогда на основании выбранной теории ползучести вычисляются приращения деформации ползучести

Др. (г,^.) за шаг времени Atj для всех г, и находятся значения р. (г,^.+1) = р. (г,ti) + Др,- (г,^). Далее с учетом а г (г) = 0 по формулам (14) определяются

Пусть в момент времени t = t * происходит температурная разгрузка от Т = Т2 до Т = Тх, при этом модуль Юнга мгновенно изменяется с Е2 на Ех. Тогда формулы (14) при Т = Тх принимают вид

Соотношения (15) и задают окончательные выражения для напряжений после ползучести (термоэкспозиции) в момент времени t = ^ + 0.

4. Экспериментальное исследование влияния термоэкспозиции на остаточные напряжения призматических образцов из сплава ЭП742 при температуре 650 °С после ультразвукового упрочнения

С целью экспериментального подтверждения изложенных выше методик расчета остаточных напряжений после процедуры упрочнения и их релаксации вследствие ползучести выполнен цикл экспериментальных исследований по влиянию термоэкспозиции на остаточные напряжения плоских образцов (прямоугольных параллелепипедов - брусков) с квадратным сечением 10 х 10 мм и длиной 100 мм из сплава ЭП742 после двустороннего ультразвукового упрочнения (УЗУ). Обработка противоположных граней образца осуществляется в специальном контейнере, наполненном микрошариками, которому через акустический трансформатор скорости (концентратор упругих колебаний) сообщаются колебания с малой амплитудой и высокой частотой (16-30 кГц). В результате хаотических соударений шариков с обрабатываемой поверхностью происходит ее упрочнение. Режимы УЗУ: расстояние от образца до концентратора (наконечника установки) - 7 мм; диаметр шариков - 1 мм; материал шариков - сталь ШХ 15; масса шариков - 4,6 г; амплитуда колебаний - 50 мкм; частота - 20 кГц; диаметр плоского концентратора - 40 мм; длительность обработки - 20 с (режим № 1), 40 с (режим № 2), 60 с (режим № 3), 80 с (режим № 4). После УЗУ часть упрочненных образцов (для всех режимов) использовалась для определения остаточных напряжений, другая часть подвергалась термоэкспозиции (температурной выдержке) при температуре 650 °С в течение 100 часов, и после температурной разгрузки определялись остаточные напряжения.

Методика определения остаточных напряжений в упрочненных образцах базируется на расчетно-экспериментальном методе полосок [1, 32]. Все образцы разрезали на электроэрозионном станке на две полоски толщиной к вдоль срединной плоскости, параллельной обработанным граням образца. Согласно методу полосок дополнительные остаточные напряжения за счет разрезки бруска учитываются измерением прогиба / (0) каждой полоски по схеме, изображенной на рис. 2, на базовой длине I. После измерения прогиба /(0) в пределах базовой длины I электрохимическим травлением удаляются слои материала толщиной а со стороны исследуемой поверхности полоски и замеряются

а,(гЛ+1)( . = Х у ).

ах (г^ * +0) = -2^ (г) + рх (г^ * +0) + (г) + ру (г^ * +0))], а у (г, t * +0) = (г) + ру (г ^ * +0) + -(( (г ) + рх (г ^ * +0))].

Еу

возникающие при этом прогибы /(а) (рис. 3). Измеренные прогибы /(0) и /(а) используются при вычислении остаточных напряжений [1, 32]:

а у (а) - уах (а) = -^/ (0)1 И - а | +

/2

+-

4Е 3/2"

(И - а)2 ^М йа

4(И - а)/(а) + 2|/£)й$

где Е - модуль Юнга; у - коэффициент Пуассона, И - толщина полоски до травления, 0 а. Состав раствора для травления и технология травления приведены в работе [20].

Рис. 2. Схема замера прогиба /

Рис. 3. Схема замера прогиба /

В случае УЗУ а у (а) = ах (а) и формула (16) принимает вид

а у (а) = ах (а) = ^|-8Е/(00 - а) +

+ •

4Е

3/2

(И - а)2 ^ - 4(И - а)/(а) + 21/(ВД

йа

Для увеличения разрешающей способности описанного метода при удалении слоев полоски травлением использовалась схема замера перемещений, основанная на консольном креплении образца (рис. 4). Связь между измеренными перемещениями полоски с удлинителем / (а)изм и прогибами / (а), используемыми в формулах (16) и (17), записывается с помощью интеграла Мора [2] в следующем виде:

/ (а) изм

4(1 + 2с / /)'

где с - длина удлинителя; I - базовая длина полоски.

Рис. 4. Схема замера перемещений / (а)и

На рис. 5 сплошными линиями приведены результаты экспериментального определения остаточных напряжений в образцах после процедуры УЗУ для всех четырех режимов обработки поверхности. Анализ этих данных свидетельствует, что сжимающие напряжения могут достигать величины минус 1100 МПа (режим № 1). Смена знака напряжений происходит на расстоянии а« 0,2 мм от поверхности. При всех режимах УЗУ наблюдается спад остаточных напряжений к поверхности образца. Увеличение времени обработки от 20 до 80 с приводит к уменьшению максимума сжимающих остаточных напряжений с повышением глубины его залегания и незначительным увеличением толщины упрочненного слоя.

На рис. 7 (сплошная линия 4) и на рис. 8 (сплошные линии) представлены экспериментальные зависимости остаточных напряжений ах = ах (г) после УЗУ при температуре 20 °С, последующей термоэкспозиции при температуре 650 °С в течение 100 часов и температурной разгрузке до 20 °С. Анализ этих данных свидетельствует, что влияние термоэкспозиции существенно сказывается на величинах остаточных напряжений, в частности, приводит к падению максимальных сжимающих напряжений в 1,4-1,6 раз. Характер распределения напряжений не претерпевает существенных изменений. Спад остаточных напряжений к поверхности сохраняется, максимальные напряжения тем выше, чем меньше время упрочнения. Примерно на 0,02 мм увеличивается глубина залегания сжимающих остаточных напряжений, причем для образцов, обработанных с большим временем, она выше.

5. Экспериментальная проверка математической модели формирования остаточных напряжений после УЗУ и их релаксации в условиях высокотемпературной ползучести

Для экспериментальной проверки методики расчета напряженно-деформированного состояния в упрочненном плоском образце (модельное представление - упрочненное полупространство) после процедуры УЗУ воспользуемся экспериментальными данными для четырех режимов упрочнения, представленных на рис. 5. Как уже отмечалось выше, обработка УЗУ относится к изотропному процессу поверхностного пластического упрочнения, поэтому в данном случае в соотношении (2) коэффициент анизотропии а =1 и в дальнейшем для оценки остаточных напряжений и пластических деформаций используются формулы (7). Поскольку упрочнение УЗУ происходит при нормальной температуре 20 °С, в (7) для сплава ЭП742 использовались справочные значения Е = 2,21-105 МПа и у = 0,3. Используя экспериментальные данные для ах =ах(г) на рис. 5, по методике, изложенной в п. 2, определялись все параметры аппроксимации (9). Для оценки первого приближения величин а0, а1 и Ь в (9) использовались значения характерных точек для

графиков г0, г *, а *, ат1п, приведенные в таблице. Расчетные значения зависимости ах =ах (г) для всех четырех режимов упрочнения приведены на рис. 5 штриховыми линиями. В последней строке таблицы приведены значения отклонения расчетных данных от экспериментальных, вычисленные на основании (12) в процентах. В целом согласование расчетных и экспериментальных данных удовлетворительное. В качестве примера на рис. 6 приведены расчетные зависимости для распределения компонент тензоров остаточных пластических деформаций для первого режима упрочнения.

-300

-600

-900

-1200

__

/ /У // // /

'// // // 1/ у

/г // /

/7 / / У / ' У

0 0,05 0,1 0,15 0,2 0,25

-300

-600

-900

-1200

// <

А // у

//

V У/

О 0,05 0,1 0,15 0,2 0,25

г, мм

г, мм

й -зоо с

-600

-900

-1200

/ / /

Л /

.....//•' 1 /I

......//

О 0,05 0,1 0,15 0,2 0,25 0,3 1, мм

в

«а -300

с

-600

-900

-1200

У

/ .....1/.. / / 1 / /

......./ /

// // // /./.......

чЧ /у

О 0,05 0,1 0,15 0,2 0,25 0,3

2, ММ

г

Рис. 5. Экспериментальные (сплошные линии) и расчетные штриховые линии эпюры остаточных напряжений ах = ах (2) в поверхностном слое (сплав ЭП742) при различных режимах УЗУ (см. таблицу): а - режим 1; б - режим 2; в - режим 3; г - режим 4

б

а

Рассмотрим теперь результаты численной реализации методики расчета релаксации остаточных напряжений в образцах после температурной выдержки при температуре 650 °С в течение времени 100 часов. Одним из основных этапов решения данной задачи является выбор теории ползучести. В работе [33] приведены экспериментальные данные по ползучести сплава ЭП742 при температуре 650 °С, из анализа которых следует, что при данной температуре существенной является деформация ползучести для тех уровней остаточных напряжений, которые возникают после процедуры упрочнения, вследствие чего и происходит релаксация остаточных напряжений.

В соответствии с [33] в качестве модели ползучести используется следующий вариант для одноосного напряженного состояния:

р(*)(*)+);

к=1

V (/) = {^к^к (О, V (0о(0 > 0,

У 1 0, Ук «)а

5 () = ак ()/ а **)п - у, (); (Г) = ()/ а **),

где ) = ^ ^=у (5) - вязкопластическая компонента деформации ползучести р(^) (описывает первую стадию ползучести); м>(Г) - деформация вязкого течения (описывает вторую стадию ползучести); 5, Xк, ак, п, с, т - параметры модели, методика идентификации которых приведена в [33]. Здесь а ** - обезразмеривающий коэффициент, который может выбираться произвольно исходя из соображений удобства.

Значения характерных точек эпюры ах = ах

и погрешностей аппроксимации для образцов из сплава ЭП742 в зависимости

от режимов УЗУ

Время обработки, с

20

Значения характерных точек эпюры а х = ах

г0, мм 0,189 0,193 0,204 0,210

г *, мм 0,037 0,040 0,058 0,073

а *, МПа -872,0 -848,5 -840,4 -840,0

°шт , МПа -1111,5 -1058,6 -1024,3 -1032,3

начения параметров аппроксимации (9

а0, МПа 119,2 103,1 140,9 155,2

а1, МПа 1230,7 1159,1 1167,3 1187,5

Ъ 0,097 0,101 0,121 0,131

Значения погрешностей аппроксимации, вычисленные на основании

Д,% 4,93 4,02 6,47 7,73

Модель (18) описывает деформацию ползучести в пределах первой и второй стадии, причем вся деформация является необратимой. В монографии [33] для модели (18) приведены следующие параметры для сплава ЭП742 при температуре 650 °С: 5 = 1,

а ** = 500 МПа, Х1 =Х = 0,022, а1 = а = 6,1 -10"3, п = 3,29, с = 0,722-10"6, т = 14,3.

Для сложного напряженного состояния с учетом 5 = 1 модель (18) обобщается следующим образом [33]:

Ру Ц) = и 1} () + wl] (г);

V»» () = (1+)р»» (г) - )+)+Р33С));

в (/) = |ХБ»» ('), Б»» ('К» (0 > 0,

»»

Б»» (г) = а(Б(г) / а * *) а»» ()/ а ** - Р»» (г);

^

а0(г) = а11(г) + а22(г) + а33(г),

где Ру (г) - тензор деформации ползучести; vij (г) и (г) - тензоры вязкопластической (необратимой) компоненты деформаций и деформации вязкого течения; Б (г) - интенсивность напряжений; | - коэффициент Пуассона для компоненты V»» (по рекомендации [33] можно использовать |" = 0,42); а**, 5, X, а, п, с, т - параметры, имеющие тот же смысл, что и в соотношениях (18) при 5 = 1. Обозначим в (19): р11 = рх, р22 = ру, р33 = рг,

а11 =ах , а22 =ау , а33 = а = ^ ^ = ^, ^2 = ^ , V33 = Vz , Р11 = Рх , Р22 = Ру , Р33 = Рг = 0;

по повторяющемуся индексу » суммирование в (19) не производится. При численной реализации приращения всех компонент деформации ползучести в соотношениях (19) вычислялись по методу Эйлера. В расчетах использовались следующие значения: при температуре Т1 = 20 °С Е = 2,21-105 МПа, при Т2 = 650 °С Е2 = 1,79-105 МПа, у = 0,3.

Подробно проанализируем результаты расчетов для первого режима упрочнения. На рис. 7 экспериментальная (сплошная линия) и расчетная (штриховая линия) зависимости, обозначенные цифрой 1, соответствуют ах = ах (2) после упрочнения при температуре 20 °С (г = 0 - 0); цифрой 2 отмечена эта же расчетная зависимость при ступенчатом изменении температуры с 20 °С до 650 °С (г = 0 + 0); цифрой 3 - расчетная зависимость ах =ах(г,Х * -0) после ползучести в течение г* = 100 часов при температуре 650 °С; цифрой 4 - распределение остаточных напряжений после температурной разгрузки с Т2 = 650 °С до Тх = 20 °С в момент времени г = г * + 0 (сплошная линия - эксперимент,

штриховая - расчет). Отклонение расчетных данных от экспериментальных в норме (12) для финишных зависимостей (г* = 100 часов) составляет А = 27 %.

На рис. 8 приведены финишные расчетные (штриховые линии) и экспериментальные (сплошные линии) зависимости ах = ах (г * + 0) после ползучести в течение 100 часов

при температуре 650 °С и последующей температурной разгрузки для режимов УЗУ 2-4. Отклонения расчетных данных от экспериментальных, рассчитанные по (12), составляют 27,4 % для режима № 2, 20,9 % - для режима № 3 и 37 % - для режима № 4.

Следует отметить, что «разница в 20 % по деформациям между кривыми ползучести, снятыми при одном и том же напряжении для данной температуры на образцах из материала одной плавки, не считается чрезмерно большой» [34, § 44, стр. 172-173], а в нашей работе специальные эксперименты на ползучесть сплава ЭП742 для исследуемых образцов не проводились, а использовались справочные данные для модели ползучести (18),

(19) из монографии [33]. Поэтому, несмотря на полученные отклонения расчетных данных от экспериментальных, результаты расчета релаксации остаточных напряжений следует признать удовлетворительными.

-Г О

2 -2

-4

-6

У 2

-300

-600

-900

0 0,05 0,1 0,15 0,2 0,25

-1200

: // 1 // , ...........'/л......'. / Л Ч /

3 4 ; п ь У/ ш у

Г \ ч _ / /

У/ // У

О 0,05 0,1 0,15 0,2 0,25

Рис. 6. Графики зависимостей qx = qx

о

^00 -600 -800

/

> / / / у

/ / /

¡У 1.........

/

0 0,05 0,1 0,15 0,2 0,25 г, мм

Рис. 7. Кинетика зависимости а х = а х

с -200 £

-400 -600 -800

/ /

/ / !/

у/ У....!

/

0 0,05 0,1 0,15 0,2 0,25 0,3 г, мм

б

а

0

й

с -200 £

^00 -600

"8000 0,05 0,1 0,15 0,2 0,25 0,3 г, мм в

Рис. 8. Экспериментальные (сплошные линии) и расчетные (штриховые линии) зависимости а х = а х (7) в момент времени г = 100 ч после ползучести и температурной разгрузки:

а - режим 2; б - режим 3; в - режим 4

у /

/* / / / // /

77 / / / /

......./ **/..... /

/

Выводы

Выполнены экспериментальные исследования по формированию остаточных напряжений после процедуры ультразвукового упрочнения призматических образцов из сплава ЭП742 и оценке степени их релаксации в условиях термоэкспозиции при температуре 650 °С в течение 100 часов. Выяснены закономерности распределения остаточных напряжений по толщине поверхностно упрочненного слоя. Показано, что термоэкспозиция в 1,4-1,6 раз уменьшает уровень сжимающих остаточных напряжений и вызывает смещение их максимума вглубь образца. Полученные экспериментальные данные играют важную самостоятельную роль в прикладных задачах, поскольку сплав ЭП742 используется для изготовления дисков турбины авиационных двигателей, эксплуатирующихся в том числе и при температуре 650 °С, а поверхностное упрочнение - штатная технология в авиадвигателестроении.

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

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

Работа выполнена при поддержке РФФИ (проект № 16-01-00249-a) и при поддержке Минобрнауки России в рамках базовой части государственного задания ФГБОУ ВО «СамГТУ» (код проекта: 1151).

Библиографический список

1. Биргер И.А. Остаточные напряжения. - М.: Машгиз, 1963. - 262 с.

2. Павлов В.Ф., Кирпичев В.А., Иванов В.Б. Остаточные напряжения и сопротивление усталости упрочненных деталей с концентраторами напряжений; Самар. науч. центр РАН. - Самара, 2008. - 64 с.

3. Остаточные напряжения и сопротивление усталости высоконагруженных резьбовых деталей / С.И. Иванов, В.Ф. Павлов, Б.В. Минин, В.А. Кирпичев, Е.П. Кочеров, В.В. Головкин; Самар. науч. центр РАН. - Самара, 2015. - 170 с.

4. Кравченко Б.А., Круцило В.Г., Гутман Г.Н. Термопластическое упрочнение - резерв повышения прочности и надежности деталей машин. - Самара: Изд-во Самар. гос. техн. ун-та, 2000. - 216 с.

5. Сулима Г.Н., Шувалов В.А., Ягодкин Ю.Д. Поверхностный слой и эксплуатационные свойства деталей машин. - М.: Машиностроение, 1988. - 240 с.

6. Марков А.И. Ультразвуковая обработка материалов. - М.: Машиностроение, 1989. - 237 с.

7. Prediction and characterization of residual stresses from laser shock peening / R.A. Brockman, W.R. Braisted, S.E. Olson, R.D. Tenaglia, A.H. Clauer, K. Langer, M.J. Shepard // International Journal of Fatigue. - 2012. - Vol. 36. - No. 1. - Р. 96-108. DOI: 10.1016/j.ijfatigue.2011.08.011

8. Dai K., Shaw L. Analysis of fatigue resistance improvements via surface severe plastic deformation // International Journal of Fatigue. - 2008. - Vol. 30. - No. 8. - Р. 1398-1408. DOI: 10.1016/j.ijfatigue.2007.10.010

9. Residual stresses and fatigue performance / M.N. James, D.J. Hughes, Z. Chen, H. Lombard, D.G. Hattingh, D. Asquith, J.R. Yates, P.J. Webster // Engineering Failure Analysis. - 2007. - Vol. 14. -No. 2. - Р. 384-395. DOI: 10.1016/j.engfailanal.2006.02.011

10. Majzoobi G.H., Azadikhah K., Nemati J. The effects of deep rolling and shot peening on fretting fatigue resistance of Aluminum-7075-T6 // Materials Science and Engineering A. - 2009. -Vol. 516. - No. 1-2. - Р. 235-247. DOI: 10.1016/j.msea.2009.03.020

11. Soady K.A. Life assessment methodologies incoroporating shot peening process effects: Mechanistic consideration of residual stresses and strain hardening: Part 1 - Effect of shot peening on fatigue resistance // Materials Science and Technology (United Kingdom). - 2013. - Vol. 29. - No. 6. - Р. 637651. DOI: 10.1179/1743284713Y.0000000222

12. Terres M.A., Laalai N., Sidhom H. Effect of nitriding and shot-peening on the fatigue behavior of 42CrMo4 steel: Experimental analysis and predictive approach // Materials and Design. - 2012. -Vol. 35. - Р. 741-748. DOI: 10.1016/j.matdes.2011.09.055

13. Применение ультразвука и взрыва при обработке и сборке / М.Ф. Вологин, В.В. Калашников, М.С. Нерубай, Б. Л. Штриков. - М.: Машиностроение, 2002. - 264 с.

14. Технологическое применение ультразвука в транспортном машиностроении / В.М. При-ходько [и др.]. - М.: Технополиграф центр, 2007. - 112 с.

15. Ультразвуковое поверхностное пластическое деформирование / В.Ф. Казанцев, Б. А. Куд-ряшов, Р.И. Нигметзянов, В.М. Приходько, Д.С. Фатюхин // Вестн. Харьк. нац. авт.-дорож. ун-та. - 2009. - № 46. - C. 7-9.

16. Гребенников М.А., Заличихие С.Д., Стебельков И.А. Физика и технология упрочнения деталей в поле ультразвука // Вестн. двигателестроения. - 2013. - № 1. - С. 72-74.

17. Александров М.К., Папшева Н.Д., Акушская О.Н. Ультразвуковое упрочнение деталей ГТД // Вестн. Самар. гос. аэрокосм. ун-та. - 2011. - № 3(27). - С. 271-276.

18. Радченко В.П., Саушкин М.Н. Математические модели восстановления и релаксации остаточных напряжений в поверхностно упрочненном слое цилиндрических элементов конструкций при ползучести // Изв. вузов. Машиностроение. - 2004. - № 11. - С. 3-17.

19. Радченко В.П., Саушкин М.Н. Ползучесть и релаксация остаточных напряжений в упрочненных конструкциях. - М.: Машиностроение-1, 2005. - 226 с.

20. Радченко В.П., Саушкин М.Н., Павлов В.Ф. Метод расчета остаточных напряжений и пластических деформаций в цилиндрических образцах с учетом анизотропии процесса поверхностного упрочнения // ПМТФ. - 2011. - Т. 52, № 2. - С. 173-182.

21. Радченко В.П., Павлов В.Ф., Саушкин М.Н. Исследование влияния анизотропии поверхностного пластического упрочнения на распределение остаточных напряжений в полых и сплошных цилиндрических образцах // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2015. - № 1. - С. 130-147. DOI: 10.15593/perm.mech/2015.1.09

22. On the effect of deeprolling and laser-peening on the stress-controlled low- and high-cycle fatigue behavior of Ti-6-Al-4V at elvated temperatures up to 550°C / I. Altenberber, R.K. Nalla, Y. Sano, L. Wagner, R.O. Ritchie // International Journal of Fatigue. - 2002. - Vol. 44. - Р. 292-302. DOI: 10.1016/j.ijfatigue.2012.03.008

23. Колотникова О.В. Эффективность упрочнения методами поверхностного пластического деформирования деталей, работающих при повышенных температурах // Проблемы прочности. -1983. - № 2. - С. 112-114.

24. Цейтлин В.И., Колотникова О.В. Релаксация остаточных напряжений в деталях турбин ГТД в процессе эксплуатации // Проблемы прочности. - 1980. - № 3. - С. 6-11.

25. Радченко В.П., Кирпичёв В.А., Лунин В.А. Влияние термоэкспозиции на остаточные напряжения образцов из сплава ЭП742 после ультразвукового упрочнения // Вестн. Самар. гос. техн. ун-та. Сер. Техн. науки. - 2012. - № 3(35). - С. 147-154.

26. Buchanan D.J., John R. Relaxation of shot-peened residual stresses under creep loading // Scripta Materialia. - 2008. - No. 3. - Р. 286-289. DOI: 10.1016/j.scriptamat.2008.03.021

27. Khadraoui M., Cao W., Castex L. Experimental investigations and modeling of relaxation behavior of shot peening residual stresses at high temperature for nickel base superalloys // Materials Science and Technology. - 1997. - Vol. 13. - No. 4. - Р. 360-367. DOI: 10.1179/026708397790302359

28. Xie L., Jiang C., Ji V. Thermal relaxation of residual stresses in shot peened surface layer of (TiB + TiC)/Ti-6Al-4V composite at elevated temperatures // Materials Science and Engineering: A. -2011. - Vol. 528. - No. 21. - Р. 6478-6489. DOI: 10.1016/j.msea.2011.04.075

29. Радченко В.П., Саушкин М.Н. Прямой метод решения краевой задачи релаксации остаточных напряжений в упрочненном изделии цилиндрической формы при ползучести // ПМТФ. -2009. - Т. 50, № 6. - С. 90-99.

30. Экспериментальное и теоретическое исследование влияния растягивающей нагрузки на релаксацию остаточных напряжений в упрочненном цилиндрическом образце в условиях ползучести / В.П. Радченко, Е.П. Кочеров, М.Н. Саушкин, В.А. Смыслов // ПМТФ - 2015. - Т. 56, № 2. - С. 169-177. DOI: 10.15372/PMTF20150217

31. Радченко В.П., Цветков В.В. Кинетика напряженно-деформированного состояния в поверхностно упрочненном цилиндрическом образце при сложном напряженном состоянии в условиях ползучести // Вестн. Самар. гос. техн. ун-та. Сер. Физ.-мат. науки. - 2014. - № 1(34). - С. 93108. DOI:10.14498/vsgtu1313

32. Иванов С.И. Определение остаточных напряжений в пластинках методом полосок // Вопросы прочности элементов авиационных конструкций; Куйб. авиац. ин-т. - Куйбышев, 1971. -С. 139-152.

33. Радченко В.П., Еремин Ю.А. Реологическое деформирование и разрушение материалов и элементов конструкций. - М.: Машиностроение-1, 2004. - 264 с.

34. Работнов Ю.Н. Ползучесть элементов конструкций. - М.: Наука, 1966. - 752 с.

References

1. Birger I.A. Ostatochnye napriazheniia [Residual Stresses]. Moscow: Mashgiz, 1963. 232 p.

2. Pavlov V.F., Kirpichev V.A., Ivanov V.B. Ostatochnye napriazheniia i soprotivlenie ustalosti uprochnennykh detalei s kontsentratorami napriazhenii [Residual Stresses and Fatigue Resistance of Hardened Parts with Stress Concentrators]. Samarskii nauchnyi tsentr Rossiiskoi akademii nauk, 2008. 64 p.

3. Ivanov S.I., Pavlov V.F., Minin B.V., Kirpichev V.A., Kocherov E.P., Golovkin V.V. Ostatochnye napriazheniia i soprotivlenie ustalosti vysokonagruzhennykh rez'bovykh detalei [Residual stresses and fatigue resistance of highly loaded threaded parts]. Samarskii nauchnyi tsentr Rossiiskoi akademii nauk, 2015. 170 p.

4. Kravchenko B.A., Krutsilo V.P., Gutman G.N. Termoplasticheskoe uprochnenie - rezerv povy-sheniia prochnosti i nadezhnosti mashin [Thermoplastic Hardening as a Reserve Increase of Strength and Reliability of Machine Parts]. Samarskii gosudarstvennyi tekhnicheskii universitet, 2000. 216 p.

5. Sulima G.N., Shuvalov V.A., Yagodkin Yu.D. Poverkhnostnyi sloi i ekspluatatsionnye svoistva detalei mashin [Surface Layer and Performance of Machine Parts]. Moscow: Mashinostroenie, 1988. 240 p.

6. Markov A.I. Ul'trazvukovaia obrabotka materialov [Ultrasonic Processing of Materials]. Moscow: Mashinostroenie, 1989. 237 p.

7. Brockman R.A., Braisted W.R., Olson S.E., Tenaglia R.D., Clauer A.H., Langer K., Shepard M.J. Prediction and characterization of residual stresses from laser shock peening. International Journal of Fatigue, 2012, vol. 36, no. 1, pp. 96-108. DOI: 10.1016/j.ijfatigue.2011.08.011

8. Dai K., Shaw L. Analysis of fatigue resistance improvements via surface severe plastic deformation. International Journal of Fatigue, 2008, vol. 30, no. 8. pp. 1398-1408. DOI: 10.1016/j.ijfatigue.2007.10.010

9. James M.N., Hughes D.J., Chen Z., Lombard H., Hattingh D.G., Asquith D.,Yates J.R., Webster P.J. Residual stresses and fatigue performance. Engineering Failure Analysis, 2007, vol. 14, no. 2, pp. 384395. DOI: 10.1016/j.engfailanal.2006.02.011

10. Majzoobi G.H., Azadikhah K., Nemati J. The effects of deep rolling and shot peening on fretting fatigue resistance of Aluminum-7075-T6. Materials Science and Engineering A, 2009, vol. 516, no. 1-2, pp. 235-247. DOI: 10.1016/j.msea.2009.03.020

11. Soady K.A. Life assessment methodologies incoroporating shot peening process effects: Mechanistic consideration of residual stresses and strain hardening: Part 1 - Effect of shot peening on fatigue resistance. Materials Science and Technology (United Kingdom), 2013, vol. 29, no. 6, pp. 637-651. DOI: 10.1179/1743284713Y.0000000222

12. Terres M.A., Laalai N., Sidhom H. Effect of nitriding and shot-peening on the fatigue behavior of 42CrMo4 steel: Experimental analysis and predictive approach. Materials and Design, 2012, vol. 35, pp. 741-748. DOI: 10.1016/j.matdes.2011.09.055

13. Vologin M.F., Kalashnikov V.V., Nerubai M.S., Shtrikov B.L. Primenenie ul'trazvuka i vzryva pri obrabotke i sborke [Application of ultrasound or explosion during processing and assembly]. Moscow: Mashinostroenie, 2002. 264 p.

14. Prikhod'ko V.M. [et al.] Tekhnologicheskoe primenenie ul'trazvuka v transportnom mashino-stroenii [The technological application of ultrasound for the transport engineering]. Moscow: Tekhnopoligraf tsentr, 2007. 112 p.

15. Kazantsev V.F., Kudriashov B.A., Nigmetzianov R.I., Prikhod'ko V.M., Fatiukhin D.S. Ul'tra-zvukovoe poverkhnostnoe plasticheskoe deformirovanie [The ultrasonic surface plastic deformation]. Vestnik Khar'kovskogo natsional'nogo avtomobil'no-dorozhnogo universiteta, 2009, no. 46, pp. 7-9.

16. Grebennikov M.A., Zalichikhie S.D., Stebel'kov I.A. Fizika i tekhnologiia uprochneniia detalei v pole ul'trazvuka [The physics and technology of details hardening in the ultrasound field]. Vestnik dvigatelestroeniia, 2013, no. 1, pp. 72-74.

17. Aleksandrov M.K., Papsheva N.D., Akushskaia O.N. Ul'trazvukovoe uprochnenie detalei GTD [The ultrasonic hardening of the turboshaft engine parts]. Vestnik Samarskogo gosudarstvennogo aerokosmicheskogo universiteta, 2011, no. 3(27), pp. 271-276.

18. Radchenko V.P., Saushkin M. N. Matematicheskie modeli vosstanovleniia i relaksatsii ostatochnykh napriazhenii v poverkhnostno uprochnennom sloe tsilindricheskikh elementov konstruktsii pri polzuchesti [Mathematical models of recovery and relaxation of residual stresses in a surface-hardened layer of cylindrical specimens under creep conditions]. Izvestiia vuzov. Mashinostroenie, 2004, no. 11, pp. 3-17.

19. Radchenko V.P., Saushkin M.N. Polzuchest' i relaksatsiia ostatochnykh napriazhenii v uproch-nennykh konstruktsiiakh [Creep and Relaxation of Residual Stresses in Hardened Structures]. Moscow: Mashinostroenie-1, 2005. 226 p.

20. Saushkin M.N., Radchenko V.P., Pavlov V.F. Method of calculating the fields of residual stresses and plastic strains in cylindrical specimens with allowance for surface hardening anisotropy. Journal of Applied Mechanics and Technical Physics, 2011, vol. 52, no. 2, pp. 303-310. DOI: 10.1134/S0021894411020180

21. Radchenko V.P., Pavlov V.Ph., Saushkin M.N. Issledovanie vliianiia anizotropii poverkhnost-nogo plasticheskogo uprochneniia na raspredelenie ostatochnykh napriazhenii v polykh i sploshnykh tsilindricheskikh obraztsakh [Investigation of surface plastic hardening anisotropy influence on residual stresses distribution in hollow and solid cylindrical specimens]. PNRPUMechanics Bulletin, 2015, no. 1, pp. 130-147. DOI: 10.15593/perm.mech/2015.1.09

22. Altenberber I., Nalla R.K., Sano Y., Wagner L., Ritchie R.O. On the effect of deeprolling and laser-peening on the stress-controlled low- and high-cycle fatigue behavior of Ti-6-Al-4V at elvated temperatures up to 550 °C. International Journal of Fatigue, 2002, vol. 44, pp. 292-302. DOI: 10.1016/j.ijfatigue.2012.03.008

23. Kolotnikova O.V. Effectiveness of hardening by methods of plastic surface deformation of components operating at high temperatures. Strength of Materials, 1983, vol. 15, no. 2, pp. 292-295. DOI: 10.1007/BF01523487

24. Tseitlin V.I., Kolotnikova O.V. Relaxation of residual stresses in gas-turbine engine parts. Strength of Materials, 1980, vol. 12, no. 8, pp. 982-984. DOI: 10.1007/BF00770521

25. Radchenko V.P., Kirpichev V.A., Lunin V.A. Vliianie termoekspozitsii na ostatochnye napriazheniia obraztsov iz splava EP742 posle ul'trazvukovogo uprochneniia [Influence of thermoexposition on residual stresses of specimens from EP742 alloy after the ultrasonic hardening]. Vestnik Samarskogo gosudarstvennogo tekhnicheskogo universiteta. Tekhnicheskie nauki, 2012, no. 3(35), pp. 147-154.

26. Buchanan D.J., John R. Relaxation of shot-peened residual stresses under creep loading. Scripta Materialia, 2008, no. 3, pp. 286-289. DOI:10.1016/j.scriptamat.2008.03.021

27. Khadraoui M., Cao W., Castex L. Experimental investigations and modeling of relaxation behavior of shot peening residual stresses at high temperature for nickel base superalloys. Materials Science and Technology, 1997, vol. 13, no. 4, pp. 360-367. DOI: 10.1179/026708397790302359

28. Xie L., Jiang C., Ji V. Thermal relaxation of residual stresses in shot peened surface layer of (TiB + TiC)/Ti-6Al-4V composite at elevated temperatures. Materials Science and Engineering: A, 2011, vol. 528, no. 21, pp. 6478-6489. DOI: 10.1016/j.msea.2011.04.075

29. Radchenko V.P., Saushkin M.N. Direct method of solving the boundary-value problem of relaxation of residual stresses in a hardened cylindrical specimen under creep conditions. Journal of Applied Mechanics and Technical Physics, 2009, vol. 50, no. 6, pp. 989-997. doi: 10.1007/s10808-009-0133-8

30. Radchenko V.P., Kocherov E.P., Saushkin M.N., Smyslov V.A. Experimental and theoretical studies of the influence of a tensile load on the relaxation of residual stresses in a hardened cylindrical specimen under creep conditions. Journal of Applied Mechanics and Technical Physics, 2015, vol. 56, no. 2, pp. 313-320. DOI: 10.1134/S0021894415020170

31. Radchenko V.P., Tsvetkov V.V. Kinetika napriazhenno-deformirovannogo sostoianiia v po-verkhnostno uprochnennom tsilindricheskom obraztse pri slozhnom napriazhennom sostoianii v uslo-viiakh polzuchesti [Kinetics of the stress-strain state of surface hardened cylindrical specimen under complex stress state of creep]. Vestnik Samarskogo gosudarstvennogo tekhnicheskogo universiteta. Seriia Fiziko-matematicheskie nauki, 2014, no. 1(34), pp. 93-108. DOI: 10.14498/vsgtu1313

32. Ivanov S.I. Opredelenie ostatochnykh napriazhenii v plastinkakh metodom polosok [Determination of residual stresses in the plates using the bent strip method]. Voprosy prochnosti elementov aviatsionnykh konstruktsii. Kuibyshevskii aviatsionnyi institut, 1971. pp. 139-152.

33. Radchenko V.P., Eremin Yu.A. Reologicheskoe deformirovanie i razrushenie materialov i ele-mentov konstruktsii [Rheological deformation and fracture of materials and structural elements]. Moscow: Mashinostroenie-1, 2004. 264 p.

34. Rabotnov Yu.N. Creep problems in structural members. Amsterdam, London, North-Holland Publ. Co., 1969, XIV+822 p.