CC BY
11
Math-Net.Ru
В. П. Радченко, А. Е. Либерман, О. Л. Блохин, Релаксация остаточных напряжений в поверхностно упрочненном вращающемся цилиндре в условиях ползучести, Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2022, том 26, номер 1, 119-139
001: https://doi.org/10.14498/vsgtu1884
Использование Общероссийского математического портала Math-Net.Ru подразумевает, что вы прочитали и согласны с пользовательским соглашением
http://www.mathnet.ru/rus/agreement
Параметры загрузки: IP: 94.29.25.84 7 июля 2022 г., 11:19:18
Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки. 2022. Т. 26, № 1. С. 119-139 ISSN: 2310-7081 (online), 1991-8615 (print) d https://doi
УДК 539.376:621.787
Релаксация остаточных напряжений в поверхностно упрочненном вращающемся цилиндре в условиях ползучести
В. П. Радченко1, А. Е. Либерман1, О. Л. Блохин2
1 Самарский государственный технический университет, Россия, 443100, Самара, ул. Молодогвардейская, 244.
2 Зения АИ,
Республика Беларусь, 220100, Минск, ул. Кульман, 9.
Аннотация
Разработана методика расчета релаксации остаточных напряжений в консольно закрепленном вращающемся цилиндре после процедуры поверхностно-пластического деформирования в условиях ползучести, учитывающая влияние ступенчатого изменения параметров температурно-силового нагружения (разгрузки). Задача моделирует напряженно-деформированное состояние поверхностно упрочненного цилиндра (стержня), торцевое сечение которого жестко закреплено на диске, вращающегося с постоянной угловой скоростью.
Методика включает метод реконструкции полей остаточных напряжений и пластических деформаций и метод расчета релаксации остаточных напряжений в процессе ползучести вращающегося цилиндрического стержня. Поскольку вызванные вращением растягивающие напряжения по длине стержня не изменяются во времени, в каждом поперечном сечении решается задача о релаксации остаточных напряжений для растягиваемого стержня при постоянном напряжении.
Выполнено детальное численное исследование влияния количества оборотов на скорость релаксации остаточных напряжений для упрочненного дробью цилиндрического образца радиусом 3.76 мм из сплава ЭИ698 при температуре 700 °C.
Анализ результатов расчетов позволил установить нетривиальный эффект, заключающийся в том, что релаксация остаточных напряжений
Научная статья
© Коллектив авторов, 2022 © СамГТУ, 2022 (составление, дизайн, макет)
3 ©® Контент публикуется на условиях лицензии Creative Commons Attribution 4.0 International (https://creativecommons.org/licenses/by/4.0/deed.ru) Образец для цитирования
Радченко В. П., Либерман А. Е., Блохин О. Л. Релаксация остаточных напряжений в поверхностно упрочненном вращающемся цилиндре в условиях ползучести // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2022. Т. 26, № 1. С. 119-139. EDN: GFBZBC. https://doi.org/10.14498/vsgtu1884. Сведения об авторах
Владимир Павлович Радченко ft https://orcid.org/0000-0003-4168-9660
доктор физико-математических наук, профессор; заведующий кафедрой; каф. прикладной
математики и информатики; e-mail: [email protected]
Александр Евгеньевич Либерман © https://orcid.org/0000-0002-9185-2131
аспирант; каф. прикладной математики и информатики; e-mail: [email protected]
Олег Леонидович Блохин© https://orcid.org/0000-0002-0167-4758
инженер по машинному обучению; e-mail: [email protected]
.org/10.14498/vsgtu1884
в сечениях, находящихся под действием осевых растягивающих напряжений вследствие вращения, происходит менее интенсивно, чем в «хвостовом» сечении, где осевая нагрузка от вращения равна нулю. Полученные в работе результаты могут быть полезными при оценке эффективности поверхностно-пластического упрочнения деталей в условиях высокотемпературной ползучести.
Ключевые слова: остаточные напряжения, поверхностное пластическое упрочнение, вращающийся цилиндр, ползучесть, релаксация.
Получение: 14 сентября 2021 г. / Исправление: 28 декабря 2021 г. / Принятие: 24 января 2022 г. / Публикация онлайн: 21 марта 2022 г.
Введение. Технологии поверхностного пластического деформирования (ППД) деталей ответственных конструкций, эксплуатирующихся в условиях температурно-силовых воздействий, являются одним из резервов повышения их ресурса вследствие создания в приповерхностном слое поля сжимающих остаточных напряжений (ОН), препятствующих развитию различного рода микродефектов в материале. Положительное влияние ОН отмечается для ряда характеристик надежности: повышается предел выносливости при много-и малоцикловом нагружении и микротвердость, улучшаются трибологиче-ские характеристики и другие показатели [1—10]. Поскольку теоретическая (расчетная) оценка этих показателей связана с наличием информации о законах распределения ОН по глубине упрочненного слоя, одной из первых задач является решение проблемы реконструкции (формирования) напряженно-деформированного состояния ОН после упрочнения. Подавляющее большинство публикаций в этом направлении посвящено экспериментальным методам определения ОН и их использованию в дальнейших расчетах в поле внешних термомеханических нагрузок. Однако эти методы позволяют определить одну или две компоненты тензора ОН и не позволяют установить компоненты тензора пластических деформаций [3,4,10-12]. Существенное продвижение в области теоретических методов реконструкции ОН связано с возросшими возможностями математического моделирования технологических процессов на основе современного программного обеспечения в многочисленных коммерческих пакетах, базирующихся на методе конечных элементов (МКЭ). В этом направлении получен ряд впечатляющих результатов. В частности, в [13] за счет варьирования технологических параметров обработки поверхности дробью (длительность обработки, скорость и направление микрошариков и др.) на основе МКЭ с помощью 3Б-моделирования удалось получить поля ОН и исследовать влияние технологических параметров на их формирование. Аналогичные подходы реализованы и в работах [14-20].
Кроме самостоятельной ценности решение проблемы реконструкции ОН имеет важное значение для постановки краевых задач ползучести поверхностно упрочненных элементов конструкций с начальным напряженно-деформированным состоянием и развития методов их решения. Величина и характер распределения ОН в процессе их релаксации являются ключевой информацией для расчета остаточного ресурса по параметрическим критериям отказа при эксплуатации упрочненных элементов конструкций по техническому состоянию.
Анализ публикаций в этом направлении свидетельствует, что тематика, связанная с релаксацией ОН в условиях ползучести, находится в стадии становления. Так, в [21] применительно к сплошным поверхностно упрочненным цилиндрическим образцам разработан прямой метод решения краевой задачи релаксации остаточных напряжений, на основании обобщения которого решен ряд задач для полых и сплошных упрочненных цилиндрических деталей в условиях сложного нагружения внешними квазистатическими нагрузками [22-25].
В настоящей работе ставится задача разработки метода решения релаксации ОН в поверхностно упрочненном сплошном цилиндрическом образце при вращении в условиях ползучести. Первая попытка решения такой задачи предпринята в [26], однако для построения методики расчета использовался приближенный метод [27], согласно которому полагалось, что поскольку толщина поверхностного упрочненного слоя мала по сравнению с радиусом цилиндрического образца, этот слой не оказывает существенного влияния на жесткость и деформируемость самого цилиндра. Поэтому тонкий упрочненный слой можно представить «наклеенным» на цилиндр и деформирующимся с ним в режиме «жесткого» нагружения под действием осевой растягивающей силы. Однако оценить погрешность этого метода теоретически не представляется возможным. В настоящей работе решение построено на основании точного (но численного) метода решения соответствующей краевой задачи.
1. Постановка задачи. Рассматривается поверхностно упрочненный сплошной вращающийся цилиндр в условиях ползучести. Целью работы является разработка метода расчета релаксации наведенных остаточных напряжений после процедуры поверхностного пластического деформирования вследствие ползучести материала цилиндра в поле массовых сил. Рассматриваемый цилиндрический стержень длиной I и радиуса а консольно закреплен в сечении г = и вращается вокруг оси АА\ с угловой скоростью ш, при этом используемая цилиндрическая система координат (г, в, х) с центром в точке О жестко связана с вращающимся упрочненным цилиндром (рис. 1). В данной постановке задача моделирует напряженно-деформированное состояние поверхностно упрочненного цилиндра (стержня), торцевое сечение которого жестко закреплено на диске радиуса Я1, вращающегося относи-
С
о
\
7ТГ\ Л
ww W
dz
N
Ri
r2
Рис. 1. Растяжение сплошного цилиндра от действия центробежной силы N при вращении
с угловой скоростью ш [Figure 1. Stretching of a solid cylinder from the action of centrifugal force N during rotation
with angular velocity w]
тельно оси АА\ с постоянной угловой скоростью ш. При вращении возникает неоднородное осевое напряженное состояние за счет переменной продольной нагрузки N (г) (рис. 1), под действием которой развиваются деформации ползучести в каждом сечении образца и происходит релаксация остаточных напряжений в тонком приповерхностном упрочненном слое. При построении расчетной модели вводится гипотеза плоских сечений.
Схема решения поставленной задачи состоит из реализации следующих этапов:
1) реконструкция полей остаточных напряжений и пластических деформаций после процедуры упрочнения при нормальной («комнатной») температуре Т0;
2) перерасчет полей остаточных напряжений при приложении темпера-турно-силовых стационарных нагрузок от вращения и изменении температуры со значения Т0 (модуль Юнга материала Е0) до температуры «эксплуатации» Т\ (Т\ > То, модуль Юнга при этой температуре Е\);
3) расчет релаксации остаточных напряжений в упрочненном цилиндрическом образце вследствие ползучести при температуре Т\.
2. Напряженно-деформированное состояние вращающегося цилиндра. Предварительно найдем закон изменения осевого напряжения во вращающемся неупрочненном цилиндре. Согласно [28] центробежная сила, действующая на элемент объема йх (рис. 1), имеет вид
сШ = ш2(Кг + х)рР (1)
где ш — угловая скорость; р = ^/д — плотность материала; 7 — удельный вес; д — ускорение свободного падения; Р(г) —площадь поперечного сечения.
При определении центробежной силы элемент цилиндра рассматривается как точечная масса, сосредоточенная в центре тяжести сечения. Тогда из (1) имеем ^
N(г) = рш2 ! № + £)Р(£№,
J X
откуда для напряжения и = а (г) получаем
^ = Ш = РШ2Ш/г (*! + 0Р(2)
Так как площадь сечения рассматриваемого кругового цилиндра постоянна по длине, из (2) находим
I2- г2"
a(z) = рШ2(R\(l — z) + )
Далее, приняв во внимание, что I = R2 — R\ — длина цилиндра, и сделав необходимые преобразования, получим окончательную формулу распределения растягивающих напряжений во вращающемся круглом цилиндре постоянного сечения при ш = const:
1
2'
( Ri + z \ 2
V R2
ф) = -pu2R2 1 — V , 0 ^ 2 ^ R2 — Ri. (3)
3. Методика реконструкции напряженно-деформированного состояния после поверхностного упрочнения. Основы этой методики заложены в работе [27] для изотропного поверхностно пластического упрочнения поверхности образца (обработка дробью, пневмо- и гидродробеструйное упрочнение и др.), а в дальнейших работах [22-25] выполнено ее обобщение на другие виды упрочнения (обкатка роликом, алмазное выглаживание и др.), приводящие к анизотропному характеру распределения остаточных пластических деформаций в направлениях осей виг цилиндрической системы координат. Обозначим через а?68, агв68, а^68 радиальное, окружное и осевое остаточные напряжения, а через дг, дв, дг — соответствующие компоненты тензора остаточных пластических деформаций после процедуры упрочнения. Недиагональными компонентами тензоров остаточных напряжений и пластических деформаций будем пренебрегать в силу их малости по сравнению с диагональными компонентами.
В предположении, что в области сжатия поверхностного слоя вторичные пластические деформации отсутствуют, с использованием экспериментально определенной компоненты ав68 = ав68(г) в работах [22-25,27] для остальных компонент тензоров остаточных напряжений и пластических деформаций получены следующие зависимости:
1 Г
а?68 = - - аГ(СШ; (4)
•(г) = -¿у'1е-1 [а"'К) + (1+а)аГ(т-
-ЩГ-) К1 -»)аГ(г)* = ; (5)
дх(г) = аде (г), дг (г) = -(1 + а) дв (г); (6)
^ = ^ [ ^ (г) - Е + °Г(г)])<1г; (7)
а?8(г) = Ео(е0г - д(г)) + ц,(аГ(г) + а^8(г)), (8)
где Ег — модуль Юнга при температуре упрочнения; у — коэффициент Пуассона; а — феноменологический параметр анизотропии упрочнения. В случае пневмодробеструйной обработки микрошариками а = 1 и эпюры напряжений ав68, а^ практически совпадают [27].
Схема расчета полей остаточных напряжений и пластических деформаций в сплошном цилиндре после процедуры упрочнения его поверхности (в момент времени Ь = 0 - 0) имеет вид [22-25,27]
аГ(г) ( а?68(г) ( дв(г) ( дг(г), дг(г) ( ег2 ( а?8(г) (9)
(числа над стрелками означают номера формул, по которым вычисляется соответствующая величина). Из схемы (9) следует, что в конечном итоге компоненты а?68, а^8, дг, дв, дг определяются через агв68 при а = 1.
Таким образом, феноменологическая модель (9) основана на имеющихся экспериментальных данных о распределении компоненты ав68 = авв8(г). Однако экспериментально ее можно определить только в тонком упрочненном
слое (области сжатия), поэтому эти данные необходимо экстраполировать на всю область 0 ^ г ^ а. Для этого можно использовать аппроксимацию вида
-а — г)
аГ(г) = ао — а1 —, (10)
где сто, ст\, Ь — параметры аппроксимации, методика идентификации которых хорошо известна [22-25,27].
4. Расчет напряженно-деформированного состояния при мгновенном приложении температурно-силового нагружения. Рассматривается следующая температурно-силовая схема нагружения образца. Сначала образец мгновенно прогревается с температуры упрочнения То (модуль Юнга Ео) до температуры «эксплуатации» Т\ (модуль Юнга Е\). Затем цилиндр начинает вращаться вокруг оси АА\ с угловой скоростью ш. Предполагается, что при температурно-силовом нагружении не возникают вторичные пластические деформации, т. е. происходит упругая догрузка упрочненного образца. Поскольку ^ = ^-г) не зависит от температуры, соотношение (5) для температуры Т\ можно записать в виде
* -" = г-" /'I + -1+
1+р Ег
Ег(1 + ар) Ео
[-1 — р)аГ(г) — раГ(г)]. (11
Соотношение (11) по форме будет аналогично (5), если все эпюры остаточных напряжений после процедуры упрочнения умножить на коэффициент Е\/Ео. Таким образом, получаем распределения остаточных напряжений при температуре Т\. Учет напряженного состояния (3), вызванного вращением, позволяет найти начальные напряжения и деформации в упрочненном цилиндре в начальный момент времени £ = 0:
^-г, г, 0) = ^аТ(г) + \рш2К1
1 — ( )2
Я2
Ео 2 " '2'
Е1,
Ео
ег(г, г, 0) = — [-1 + р)аг(г, г, 0) — ра*-г, г, 0)] + дг-г).
Ог-г, г, 0) = — аХ*-г), г = г, в;
Е\
Здесь и далее через а^ = Стг-г,г,1), ^ = £г-г,г,1), г = г, в, г, обозначаются зависимости для компонент тензоров напряжений и полных деформаций; ст*-г, г, {) = аг-г, г, {) + -г, г, {) + ах-г, г, {).
5. Методика расчета кинетики напряжений во вращающемся поверхностно упрочненном цилиндре в условиях ползучести. Решить поставленную задачу аналитически невозможно в силу нелинейности любого закона ползучести для металлических материалов, а получение компонент тензоров напряжений (ст^, атгев) и пластических деформаций (дг, дг, д$) из (4)-(8) возможно только численными методами. С учетом неоднородности напряженного состояния по координате г € [0,1] при вращении выполняется
дискретизация цилиндра по этой координате: 0 = z0 < Z\ < ■ ■ ■ < zn = l с постоянным шагом А z = I/N. В дальнейшем упрочненный образец рассматривается как составной стержень из N элементарных цилиндрических стержней с образующей А z, при этом в пределах каждого из них напряженное состояние по координате z (Zi-\ ^ z ^ Zi, г = 1,N) можно считать однородным, а осевое растягивающее напряжение az(z), определяемое по (3),— постоянным.
В связи с этим расчет релаксации остаточных напряжений в каждом сечении вращающегося упрочненного цилиндра можно выполнить автономно. Для этой цели можно использовать разработанную ранее методику [21,22] для расчета кинетики напряженно-деформированного состояния в поверхностно упрочненном цилиндре в условиях ползучести при постоянной растягиваемой нагрузке (постоянном напряжении az (z), определяемого из соотношения (3)).
Постановка краевой задачи включает следующие соотношения для к-го элементарного стержня (к = 1,N): - уравнения равновесия
dar (г, zk, t) , . , . . .
г-d;,--+ar(r,Zk, ч = ae(r, zk, t); (12)
Г az (r, Zk, t)rdr = , (13)
J 0
где ar(r, Zk, t), ae(r, Zk, t), az(r, Zk, t) —радиальная, окружная и осевая компоненты тензора напряжений в цилиндрическом образце в сечении ^ = Zk;
- уравнение совместности деформаций
r^M + ^(r, Zk, t) = er(r, Zk, t), (14)
где r( , k, ), e( , k, ) — радиальная и окружная компоненты тензора полных деформаций;
- гипотеза плоских сечений:
£z(Г, Zk, t) = e*z(Zk, t); (15)
- краевые условия:
dar (r, Zk, t)
lim---= 0; ar(a, zk, t) = 0. (16)
r^o dr
Используя закон Гука, соотношение (15) и зависимость для полной деформации
£i(г, Zk, t) = ei(г, Zk, t) + qi(r)+ pi(r, Zk, t), i = r,d,z, (17)
где i и i — компоненты упругой деформации и деформации ползучести, находим зависимость для осевой компоненты тензора деформаций:
az (r, Zk, t) = Ei[e** (zk, t) - Qz (r) - Pz (r, Zk, +ß[ae (r, Zk, t) + ar (r, Zk, , (18)
при этом для компонент деформаций ползучести имеем начальные условия Рг(г, хк, 0) = 0, г = г, в, х.
В работе [21] (и последующих работах [22-25]) на основе уравнений (6), (12)-(18) получено дифференциальное уравнение для компоненты «г(г, Хк, ^ в виде
г2 (12аг (г,хк, I) + зг = д(г ^ ^
с граничными условиями (16), где
Ех 1 -/л2
2 + а
д(г, Хк, ^ = --2 —— (г) + рг ( г, Хк, £ ) - рв(г, Хк, £ ) -
11 + а
'йрв(г, Хк, ^ <рх (г, Хк, 1) \ г <!дг (г)
- П-ХТ^" + Р-+ + ар)
<1г <1г ) 1 + а <г
решение которого записывается следующим образом:
га 1 (
/0
Далее из (12) с учетом (19) находим окружную компоненту напряжений:
<«г (г, хк, 1)
<гг (г, Хк, Ь) = - £ ^з^ 9(V, , ^ф(19)
<гв(г, Хк, ^ = аг(г, Хк, ^ + г-
<
Для определения (г, Хк, ^ из (18) необходимо знать е*(Хк, £). Подставляя (18) в уравнение равновесия (13), получаем уравнение относительно величины е*(Хк, £), выражая которую, находим
£**(Хк, £) = ^ + ^ У ^(0+ Рг(Г, Хк, ^ - Е(«г(Г, Хк, г) + ав(г, Хк, г))^ гйг,
™ 1 ° 0 (20) а далее с использованием (20) определяем (г, Хк, ^ по формуле (18).
Расчет осуществляется на заданном временном интервале Ь € [0, ¿*], где ¿* — финишное значение. Далее при £ = ¿* осуществляется температурная разгрузка образца с температуры эксплуатации Тх до температуры Т0 (предполагается, что этот процесс осуществляется мгновенно). В этом случае полученные к моменту времени зависимости для остаточных напряжений умножаются на коэффициент Е0/Ех. Затем производится силовая разгрузка образца, т. е. полагается ш = 0. Здесь в каждом сечении образца г = Хк происходит ступенчатое изменение осевой компоненты тензора напряжений на величину, определяемую соотношением (3). Компоненты «г(г, Хк, ) и «в(г, Хк, ) реагируют лишь на ступенчатое изменение температурного режима с Тх до То. В итоге получаем поля остаточных напряжений при температуре Т0 после ползучести за время ¿*.
6. Выбор феноменологической модели ползучести. Для реализации методики расчета релаксации остаточных напряжений в условиях ползучести требуется использование какой-либо теории ползучести. В данной работе используется реологическая модель, предложенная для описания первой и второй стадий ползучести в работе [29] и обобщенная на третью стадию
ползучести [30], основные соотношения которой при сложном напряженном состоянии имеют следующий вид:
^гз — + Р^З •
р.. — 1±Иа.. _ а • с-г] — Е °г3 Еи%Зи •
Р^з — йгз + ^г^ + мг^ •
/ -\ т1—1 1 /з 1 Г \
йк(I) — \к ак(-)п2-11 [(1 + ^)аг3 - - (I)
'¡к-
(щи № —£ ^ (^, к
»к* ад — (1+ж) & ад - +рШ+Ркз ад)
ёки ад —
хк
111ь; т И22\Ч "Т" Н33\ Ък{-)П2-1 -Рки , [••• \ аии > 0, 0, [■■■ ]аии < 0;
^ — а% (1 + ш); ш — а(Б0 .
(21) (22)
(23)
(24)
(25)
(26) (27)
Здесь е^, е^, р^ — полная, упругая деформация и деформация ползучести соответственно; йг3, , — вязкоупругая, вязкопластическая и вязкая составляющие деформации ползучести; аг3, а^ —соответственно компоненты истинного и номинального тензоров напряжений; Е, у — упругие константы материала; 5, Б — соответственно интенсивности тензоров истинных и номинальных напряжений; Хк, ак, Ьк, с, П2, Ш1, а* — константы модели, при помощи которых описываются первая и вторая стадии ползучести материала и ее обратимая после разгрузки часть; , — коэффициенты Пуассона для обратимой и необратимой компонент деформаций ползучести; рк — активные вязкопластические деформации, наблюдаемые при отсутствии пуас-соновского сужения материала; Б) аппроксимируется соотношением вида 7(5о) — , где а1, т3 — константы модели, описывающие разупрочнение
материала вследствие деформации ползучести; по повторяющемуся индексу выполняется суммирование.
Поскольку вязкопластическая деформация рассчитывается в главных осях, отсутствует суммирование по индексу V в формулах (26).
В работе [30] приведен критерий рассеянного (объемного) разрушения энергетического вида для определения времени разрушения материала ¿*:
Г** аг^Аргу 1о А? (Б0)
— 1,
(28)
где величина А? в общем случае имеет аппроксимацию вида А?(5о) — ,
где а а и ша — постоянные параметры модели. Отметим, что все параметры реологической модели могут быть найдены при наличии одоосных экспериментальных кривых ползучести вплоть до разрушения по методике работы [30]. В этой же работе приведены реологические параметры для некоторых материалов.
7. Результаты расчетов и их анализ. В модельных расчетах использовался цилиндрический поверхностно упрочненный образец радиуса a = = 3.76 мм из сплава ЭИ698 (ХН73МБТЮ) после процедуры дробеструйной обработки (в соотношениях (6) величина a = 1) при нормальной («комнатной») температуре То = 26 С (модуль Юнга материала при этой температуре Ео = 2 ■ 105 МПа), помещенный затем в поле массовых сил от вращения при Ti = 700С (модуль Юнга Ei = 1.52 ■ 105 МПа). Численные значения геометрических параметров (см. рис. 1) следующие: Ei = 517 мм, = 667 мм.1 Значение коэффициента Пуассона р = 0.33 полагалось независимым от температуры. Расчеты для оценки кинетики остаточных напряжений во вращающемся образце в условиях ползучести выполнялись в пяти сечениях по длине I = 150 мм при значениях z = 0 (сечение 1), z = 37.5 мм (сечение 2), z = 75 мм (сечение 3), z = 112.5 мм (сечение 4) и z = 150 мм (сечение 5) при пяти величинах угловой скорости ш = {1500,1750, 2000, 2250, 2500} (об/мин). В каждом сечении рассчитывалась величина осевого напряжения a = o(z) от вращения по формуле (3) (в расчетах для сплава ЭИ698 использовалось значение 7 = 8.16 г/см3), которое остается постоянным во времени при фиксированной величине координаты z, и при этих значениях осуществлялся расчет релаксации остаточных напряжений по приведенной выше методике. В качестве примера на рис. 2 приведены графики осевых напряжений по длине образца в зависимости от угловой скорости в соответствии с формулой (3).
z, cm
Рис. 2. Распределения осевой нагрузки a(z) по длине цилиндрического образца от угловой скорости ш\ 1 — ш = 1500 об/мин; 2 — ш = 1750 об/мин; 3 — ш = 2000 об/мин; 4 — ш =
= 2250 об/мин; 5 — ш = 2500 об/мин
[Figure 2. Distributions of the axial load <t(z) along the length of the cylindrical sample on the angular velocity ш\ 1 — ш = 1500 rev/min; 2 — ш = 1750 rev/min; 3 — ш = 2000 rev/min; 4 — ш = 2250 rev/min; 5 — ш = 2500 rev/min]
В силу неоднородности полей остаточных напряжений не только по координате г, но и по координате г, осуществляется дискретизация по переменной г (0 = г0 < г1 < г2 < ■ ■ ■ < гм = а) с шагом Дп = п — п-1, i = 1, М. Реализация рассмотренной схемы решения в условиях ползучести осуществ-
Значения R\ и R2 соответствуют длине лопатки I = R2 — Ri = 150 мм и радиусу диска Ri одной из ступеней газотурбинного двигателя.
ляется численно шагами по времени. Для этого проводится дискретизация по переменной t (ti = t0 + ¿At, i = 0,1, 2,...; to = 0) и расчет компонент деформации ползучести на основании реологической модели (21)—(28) выполняется методом Эйлера численного интегрирования дифференциальных уравнений.
Экспериментальные данные для диаграммы одноосной ползучести сплава ЭИ698 при температуре 700 °С приведены в работе [30]. На основании этих данных в этой же работе вычислены все параметры реологической модели (21)—(28), численные значения которых приведены в таблице.
Значения параметров модели (21)—(28) для описания деформации ползучести сплава ЭИ698 (ХН73МБТЮ) при температуре 700 °С [30] [Values of model parameters (21)-(28) for describing the creep deformation of the EI698 (KhN73MBTYu) alloy at a
temperature of 700 С [30]]
a*, MPa k Afc, h-1 Ik , x10-4 bk, x10-4 c, x10-5 П2 m1 tt1, MPa-1-m1 m3 a a, MPa-1-mA mA
490.5 1 0.2 2.96 4.44 2.51 2.9 10.96 9.56 • 103 2.03 12.2 0
На первом этапе выполняется реконструкция полей остаточных напряжений после упрочнения образца из сплава ЭИ698 дробеструйной обработкой поверхности. Как указано выше (см. п. 3), исходной информацией является экспериментальная информация о распределении величины ст^^Н) (Н = а — г — глубина упрочненного слоя), которая приведена в монографии [27] для образца такого же радиуса а = 3.76 мм и представлена на рис. 3 точками. На основе этих экспериментальных данных найдены параметры аппроксимации (10): а0 = 11.58 МПа; а\ = —611.58 МПа; Ь = 0.8 ■ 10-3 м. На основе аналитической зависимости для &ЦеБ(г) по схеме (9) с использованием зависимостей (4)-(8) получены законы распределения для а^ез(Н), ст£е£|(Н) и Цг(К), г = г,0,х. На рис. 3 сплошными линиями представлены расчетные значения для всех компонент тензора остаточных напряжений. Как следует из представленных экспериментальных данных, зависимости для ах (Н) и (Н) близки, а величина аг(Н) в области сжатия материала на 3-4 порядка (по модулю) меньше значений для аг (Н) и (Н).
о
£ -200 sb -400 -600
0 0.05 0.1 0.15 0.2
h, mm
Рис. 3. Расчетные компоненты тензора остаточных напряжений (сплошные линии) и экспериментальные данные (точки) для упрочненного образца из сплава ЭИ698 (ХН73МБТЮ) по глубине упрочненного слоя h = а — г [Figure 3. Calculated components of the residual stress tensor (solid lines) and experimental data (dots) for a hardened sample from EI698 (KhN73MBTYu) alloy over the depth of the
hardened layer (h = a — r)[
В качестве примера детально рассмотрим результаты расчетов релаксации остаточных напряжений в цилиндрическом образце при ш — 2000 об/мин в сечениях, соответствующих г — 0 («корневое» сечение), г — 75 мм («центральное» сечение — сечение 3), г — 150 мм («хвостовое» сечение). Осевые напряжения от вращения, рассчитанные по формуле (3), имеют следующие значения: аг(0) — 317.9 МПа; аг(75) — 169.03 МПа; аг(150) — 0. На рис. 4 приведены зависимости для остаточных напряжений а в — а в (г, Ь) в процессе ползучести упрочненного образца при температуре 700 С в этих сечениях по глубине упрочнения к — а — г в области сжатия материала на временной базе £ € [0, 300] час, а на рис. 5 представлена аналогичная информация для осевой компоненты тензора остаточных напряжений ах — ах (г, ¿).
Отметим, что во всех расчетах при использовании критерия разрушения (28) его левая часть всегда принимала значения меньшие единицы, т.е. разрушения материала не наблюдалось.
Анализ зависимостей на рис. 4 и 5 свидетельствует, что релаксация остаточных напряжений в сечениях, находящихся под действием осевых растягивающих напряжений вследствие вращения, происходит менее интенсивно, чем в «хвостовом» сечении при г — 150 мм, где осевая нагрузка от вращения равна нулю. Это следует из сравнения финишных кривых с маркером 4 для а& — ае(г, Ь) с исходной эпюрой с маркером 1 (см. рис. 4, а и 4, с) и кривых для ах — ах(г, Ь) на рис. 5, а и 5, с (маркеры 1 и 6 соответственно). Другими словами, приложение растягивающих осевых напряжений приводит к замедлению процесса релаксации остаточных напряжений по сравнению с условиями чистой термоэкспозиции в сечении г — 150 мм, где ползучесть происходит только под действием остаточных (собственных) напряжений.
Отметим, что аналогичный нетривиальный вывод получен и в работе [31], в которой показано, что при всестороннем растяжении поверхностно упрочненной пластины в условиях плоского напряженного состояния наблюдается замедление процесса релаксации остаточных напряжений при ползучести, в отличие от условий термоэкспозиции (чисто температурное нагружение) без внешних силовых нагрузок.
Второй важный вывод состоит в том, что за время ползучести в течение 300 часов остаточные напряжения (по модулю) уменьшились на 16-20 % в сечениях по длине вращающегося образца, что следует из сравнения кривых с маркерами 1 и 4 на рис. 4 и с маркерами 1 и 6 на рис. 5. Аналогичные результаты получены и при других рассмотренных значениях числа оборотов. С точки зрения эффективности поверхностно пластического упрочнения деталей, эксплуатирующихся в условиях высокотемпературной ползучести (по крайней мере при температуре 700 С и угловой скорости вращения 2000 об/мин), эти результаты носят позитивный характер.
Полученные результаты могут быть использованы в авиадвигателестрое-нии для оценки эффективности поверхностно пластического упрочнения деталей авиационной техники, поскольку никелевый сплав ЭИ698 (ХН73МБТЮ) широко применяется для изготовления лопаток газовых турбин, термообра-ботанных и обточенных штамповок дисков, дефлекторов, лабиринтов и других деталей с рабочей температурой до 750 С, а также крепежных и других деталей, работающих при уровне температур порядка 800 С.
h, mm а
h, mm b
h, mm с
Рис. 4. Кинетика компоненты остаточных напряжений a в (r,t) в условиях ползучести при ш = 2000 об/мин в сечениях г = 0 (a), г = 75 мм (b) и г = 150 мм (с). Маркеры: 1 — после процедуры упрочнения при Т = 26 С в момент времени t = 0 — 0; 2 — после температурно-силовой нагрузки при Т = 700 С в момент времени t = 0 + 0; 3 — после ползучести при температурно-силовой нагрузке при Т = 700 С в момент времени t = 300 — 0 ч; 4 — после температурно-силовой разгрузки при Т = 26 С в момент времени t = 300 + 0 ч
[Figure 4. Kinetics of the residual stress component a g (r,t) (calculated values) under creep conditions at ш = 2000 rev/min in sections г = 0 (a), г = 75 mm (b), and г = 150 mm (c). Labels: 1 — after the hardening procedure at T = 26 С at time t = 0 — 0; 2 — after temperature-force load at T = 700 С at time t = 0 + 0; 3 — after temperature-force load at T = 700 С under creep conditions at time t = 300 — 0 h; 4 — after temperature-force unloading at T = 26 С at
time t = 300 + 0 h]
h, mm а
h, mm b
« -200 cu
ь -400
-600
2,3
\1
0.05
0.1 h, mm
0.15
0.2
Рис. 5. Кинетика компоненты остаточных напряжений аг (г, £) в условиях ползучести при ш = 2000 об/мин в сечениях г = 0 (а), г = 75 мм (Ь) и г = 150 мм (с). Маркеры: 1 — расчетные значения после процедуры упрочнения при Т = 26 С в момент времени Ь = 0 — 0; 2 — после температурной нагрузки (Т = 700 С) в момент времени t = 0 + 0; 3 — после силовой нагрузки при Т = 700 С в момент времени £ = 0 + 0; 4 — после ползучести при темпера-турно-силовой нагрузке при Т = 700 С в момент времени £ = 300 — 0 ч; 5 — после силовой разгрузки при Т = 700 С в момент времени £ = 300 + 0 ч; 6 — после температурной разгрузки (Т = 26 С) в момент времени £ = 300 + 0 ч
[Figure 5. Kinetics of the residual stress component az (r,t) (calculated values) under creep conditions at ш = 2000 rev/min in sections г = 0 (a), г = 75 mm (b), and г = 150 mm (c). Labels: 1 — after the hardening procedure at T = 26 °C at time t = 0 — 0; 2 — after thermal load (T = 700 C) at time t = 0 + 0; 3 — after force load at T = 700 C at time t = 0 + 0; 4 — after temperature-force load at T = 700 C under creep conditions at time t = 300 — 0 h; 5 — after force unloading at T = 700 C at time t = 300 + 0 h; 6 — after thermal unloading (T = 26 C) at
time t = 300 + 0 h]
Заключение. В работе разработана методика расчета релаксации остаточных напряжений в консольно закрепленном вращающемся поверхностно упрочненном цилиндрическом образце в условиях ползучести с начальным напряженно-деформированным состоянием после процедуры упрочнения, учитывающая влияние ступенчатого изменения параметров температур-но-силового нагружения (разгрузки) и на распределение остаточных напряжений и их кинетику в процессе ползучести.
Выполнены модельные расчеты для оценки скорости релаксации остаточных напряжений в зависимости от количества оборотов для упрочненного дробью цилиндрического образца радиусом 3.76 мм из сплава ЭИ698 при температуре 700 С. Анализ результатов расчетов позволил установить нетривиальный эффект, заключающийся в том, что релаксация остаточных напряжений в сечениях, находящихся под действием осевых растягивающих напряжений вследствие вращения, происходит менее интенсивно, чем в «хвостовом» сечении, где осевая нагрузка от вращения равна нулю. Полученные в работе результаты могут быть полезны при оценке эффективности поверхностно-пластического упрочнения деталей в условиях высокотемпературной ползучести.
Конкурирующие интересы. Заявляем, что в отношении авторства и публикации этой статьи конфликта интересов не имеем.
Авторский вклад и ответственность. Все авторы принимали участие в разработке концепции статьи и в написании рукописи. Авторы несут полную ответственность за предоставление окончательной рукописи в печать. Окончательная версия рукописи была одобрена всеми авторами.
Финансирование. Работа выполнена при поддержке Российского научного фонда (РНФ 19-19-00062, Самарский государственный технический университет).
Библиографический список
1. Биргер И. А. Остаточные напряжения. М.: Машгиз, 1963. 232 с.
2. Гринченко И. Г. Упрочнение деталей из жаропрочных и титановых сплавов. М.: Машиностроение, 1971. 120 с.
3. Сулима А. М., Шувалов В. А., Ягодкин Ю. Д. Поверхностный слой и эксплуатационные свойства деталей машин. М.: Машиностроение, 1988. 240 с.
4. Ножницкий Ю. А., Фишгойт А. В., Ткаченко Р. И., Теплова С. В. Разработка и применение новых методов упрочнение деталей ГТД, основанных на пластическом деформировании поверхностных слоев // Вестник двигателестроения, 2006. №2. С. 8-16.
5. Dai K., Shaw L. Analysis of fatigue resistance improvements via surface severe plastic deformation// Intern. J. Fatigue, 2008. vol.30, no. 8. pp. 1398-1408. https://doi.org/ 10.1016/j.ijfatigue.2007.10.010.
6. James M. N., Hughes D. J., Chen Z., et al. Residual stresses and fatigue performance // Eng. Failure Anal., 2007. vol. 14, no. 2. pp. 384-395. https://doi.org/10.1016/j.engfailanal. 2006.02.011.
7. Majzoobi G. H., Azadikhah K., Nemati J. The effect of deep rolling and shot peening on fretting fatigue resistance of Aluminum-7075-T6 // Mater. Sci. Eng. A, 2009. vol. 516, no. 1-2. pp. 235-247. https://doi.org/10.1016/j.msea.2009.03.020.
8. Soady K. A. Life assessment methodologies incorporating shot peening process effects: mechanistic consideration of residual stresses and strain hardening. 1. Effeact of shot peening on fatigue resistance// Mater. Sci. Technol., 2013. vol.29, no. 6. pp. 673-651. https://doi. org/10.1179/1743284713Y.0000000222.
9. 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// Mater. Design., 2012. vol.35. pp. 741-748. https://doi.org/10.1016/j.matdes.2011.09.055.
10. Павлов В. Ф., Кирпичев В. А., Вакулюк В. С. Прогнозирование сопротивления усталости поверхностно упрочненных деталей по остаточным напряжениям. Самара: Самар. науч. центр РАН, 2012. 125 с.
11. Радченко В. П., Павлов В. Ф., Саушкин М. Н. Исследование влияния анизотропии поверхностного пластического упрочнения на распределение остаточных напряжений в полых и сплошных цилиндрических образцах // Вестн. Перм. нац. иссл. политехн. унта. Механика, 2015. №1. С. 130-147. https://doi.org/10.15593/perm.mech/2015.1.09.
12. Radchenko V. P., Saushkin M. N., Bochkova T. I. Mathematical modeling and experimental study of forming and relaxation of residual stresses in plane samples made of EP742 alloy after ultrasonic hardening under high-temperature creep conditions // PNRPU Mechanics Bulletin, 2018. no. 3-4. pp. 88-98. https://doi.org/10.15593/perm.mech/eng.2018.3.09.
13. Chen H., Wang S., Lu S., et al. Simulation and experimental validation of residual stress and surface roughness of high manganese steel after shot peening // Procedia CIRP, 2018. vol.71. pp. 227-231. https://doi.org/10.1016/j.procir.2018.05.066.
14. Isa M. R., Sulaiman S. N., Zaroog O. S. Experimental and simulation method of introducing compressive residual stress in ASTM A516 grade 70 steel// Key Eng. Mater., 2019. vol. 803. pp. 27-31. https://doi.org/10.4028/www.scientific.net/KEM.803.27.
15. Киселев И. А., Жуков Н. А., Васильев Б. Е., Селиванов А. Н. Учет остаточных напряжений при расчетах прочности элементов замковых соединений. Часть 1. Моделирование дробеструйной обработки // Известия высших учебных заведений. Машиностроение, 2018. №11. С. 49-59. https://doi.org/10.18698/0536-1044-2018-11-49-59.
16. Meguid S. A., Maicic L. A. Finite element odeling of shot peening residual stress relaxation in turbine disk assemblies // J. Eng. Mater. Technol., 2015. vol. 137, no. 3, 031003. https:// doi.org/10.1115/1.4030066.
17. Gallitelli D., Boyer V., Gelineau M., et al. Simulation of shot peening: From process parameters to residual tress fields in a structure // Comptes Rendus Mécanique, 2016. vol. 344, no. 4-5. pp. 355-374. https://doi.org/10.1016/j.crme.2016.02.006.
18. Zimmerman M., Klemenz M., Schulze V. Literature review on shot peening simulation// Int. J. Comput. Mater. Sci. Surf. Eng., 2010. vol.3, no. 4. pp. 289-310. https://doi.org/ 10.1504/ijcmsse.2010.036218.
19. Purohil R., Verma C. S., Rana R. S., et al. Simulation of shot peening process// Materials Today: Proceedings, 2017. vol.4, no. 2 A. pp. 1244-1251. https://doi.org/10.1016/j. matpr.2017.01.144.
20. Jebahi M., Gakwaya A., Lévesque J., et al. Robust methodology to simulate real shot peening process using discrete-continuum coupling method// Int. J. Mech. Sci., 2016. vol. 107. pp. 21-33. https://doi.org/10.1016/j.ijmecsci.2016.01.005.
21. Радченко В. П., Саушкин М. Н. Прямой метод решения краевой задачи релаксации остаточных напряжений в упрочненном изделии цилиндрической формы при ползучести // ПМТФ, 2009. Т. 50, №6. С. 90-99.
22. Радченко В. П., Кочеров Е. П., Саушкин М. Н., Смыслов В. А. Экспериментальное и теоретическое исследование влияния растягивающей нагрузки на релаксацию остаточных напряжений в упрочненном цилиндрическом образце в условиях ползучести // ПМТФ, 2015. Т. 56, №2. С. 169-177. https://doi.org/10.15372/PMTF20150217.
23. Деревянка Е. Е., Радченко В. П., Цветков В. В. Релаксация остаточных напряжений в поверхностно упрочненном цилиндре в условиях ползучести при жестких ограничениях на линейные и угловые деформации// Изв. РАН. МТТ, 2021. №3. С. 118-127. https://doi.org/10.31857/S057232992103003X.
24. Радченко В. П., Цветков В. В., Саушкин М. Н. Релаксация остаточных напряжений в упрочненном цилиндре в условиях ползучести при нагружении осевой силой, крутящим моментом и внутренним давлением// ПМТФ, 2020. Т. 61, №4. С. 96-107. https:// doi.org/10.15372/PMTF20200412.
25. Радченко В. П., Саушкин М. Н., Цветков В. В. Влияние термоэкспозиции на релаксацию остаточных напряжений в упрочненном цилиндрическом образце в условиях ползучести// ПМТФ, 2016. Т. 57, №3. С. 196-207. https://doi.org/10.15372/ PMTF20160320.
26. Просвиркина Е. А., Саушкин М. Н. Релаксация остаточных напряжений в поверхностно упрочненном слое сплошного вращающегося цилиндра в условиях ползучести / Труды Третьей Всероссийской научной конференции (29-31 мая 2006 г.). Часть 1: Математические модели механики, прочности и надежности элементов конструкций / Матем. моделирование и краев. задачи. Самара: СамГТУ, 2006. С. 192-199.
27. Радченко В. П., Саушкин М. Н. Ползучесть и релаксация остаточных напряжений в упрочненных конструкциях. М.: Машиностроение-1, 2005. 226 с.
28. Биргер И. А., Шорр Б. Ф., Иосилевич Г. Б. Расчет на прочность деталей машин. М.: Машиностроение, 1979. 702 с.
29. Самарин Ю. П. Уравнения состояния материалов со сложными реологическими свойствами. Куйбышев: Куйб. гос. унив., 1979. 84 с.
30. Радченко В. П., Еремин Ю. А. Реологическое деформирование и разрушение материалов и элементов конструкций. М.: Машиностроение-1, 2004. 265 с.
31. Радченко В. П., Бербасова Т. И., Шишкин Д. М. Релаксация остаточных напряжений в поверхностно упрочненном призматическом образце в условиях ползучести при двухосном нагружении// ПМТФ, 2021. Т. 62, №5. С. 184-194. https://doi.org/ 10.15372/PMTF20210518.
Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki
[J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2022, vol. 26, no. 1, pp. 119-139
d https://doi.org/10.14498/vsgtu1884
ISSN: 2310-7081 (online), 1991-8615 (print)
MSC: 74A10, 74D10
Relaxation of residual stresses in a surface-hardened rotating cylinder under creep conditions
V. P. Radchenko1, A. E. Liberman1, O. L. Blokhin2
1 Samara State Technical University,
244, Molodogvardeyskaya st., Samara, 443100, Russian Federation.
2 Zenia-AI,
9, Kulman st., Minsk, 220100, Belarus.
Abstract
A technique for calculating the relaxation of residual stresses in a can-tilevered rotating cylinder after the procedure of surface plastic deformation under creep conditions has been developed, taking into account the effect of a stepwise change in the parameters of temperature-force loading (unloading). The problem simulates the stress-strain state of a surface-hardened cylinder (rod), the end section of which is rigidly fixed on a disk rotating at a constant angular velocity.
The technique includes a method for reconstructing the fields of residual stresses and plastic deformations and a method for calculating the relaxation of residual stresses during creep of a rotating cylindrical rod. Since the tensile stresses caused by rotation along the length of the rod do not change in time, the problem of relaxation of residual stresses for a stretched rod at constant stress is solved in each cross section.
A detailed numerical study of the effect of the number of revolutions on the rate of relaxation of residual stresses was performed for a shot-hardened cylindrical sample with a radius of 3.76 mm made of EI698 alloy at a temperature of 700°C.
Research Article
© Authors, 2022
© Samara State Technical University, 2022 (Compilation, Design, and Layout)
3 ©® The content is published under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/)
Please cite this paper in press as:
Radchenko V. P., Liberman A. E., Blokhin O. L. Relaxation of residual stresses in a surface-hardened rotating cylinder under creep conditions, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2022, vol. 26, no. 1, pp. 119-139. https://doi.org/10.14498/vsgtu1884 (In Russian). Authors' Details:
Vladimir P. Radchenko A https://orcid.org/0000-0003-4168-9660
Dr. Phys. & Math. Sci., Professor; Head of Dept.; Dept. of Applied Mathematics & Computer Science; e-mail: [email protected]
Aleksandr E. Liberman © https://orcid.org/0000-0002-9185-2131 Postgraduate Student; Dept. of Applied Mathematics & Computer Science; e-mail: [email protected]
Oleg L. Blokhin © https://orcid.org/0000-0002-0167-4758 Machine Learning Engineer; e-mail: [email protected]
Analysis of the calculation results allowed to establish a non-trivial effect, which consists in the fact that the relaxation of residual stresses in sections subjected to axial tensile stresses due to rotation occurs less intensively than in the "tail" section, where the axial load from rotation is zero. The results obtained in this work can be useful in evaluating the effectiveness of surface-plastic hardening of parts under high-temperature creep conditions.
Keywords: residual stresses, surface plastic hardening, rotating cylinder, creep, relaxation.
Received: 14th September, 2021 / Revised: 28th December, 2021 / Accepted: 24th January, 2022 / First online: 21st March, 2022
Competing interests. We declare that we have no competing interests.
Author's Responsibilities. We take full responsibility for submit the final manuscript
to print. We approved the final version of the manuscript.
Funding. This study was supported by the Russian Science Foundation (RSF 19-19-00062,
Samara State Technical University).
References
1. Birger I. A. Ostatochnye napriazheniia [Residual Stresses]. Moscow, Mashgiz, 1963, 232 pp. (In Russian)
2. Grinchenko I. G. Uprochnenie detalei iz zharoprochnykh i titanovykh splavov [Hardening Parts Made of High-Resistant and Titanium Alloys]. Moscow, Mashinostroenie, 1971, 120 pp. (In Russian)
3. Sulima G. N., Shuvalov V. A., Yagodkin Yu. D. Poverkhnostnyi sloi i ekspluatatsion-nye svoistva detalei mashin [Surface Layer and Performance of Machine Parts]. Moscow, Mashinostroenie, 1988, 240 pp. (In Russian)
4. Nozhnitskii Yu. A., Fishgoit A. V., Tkachenko R. I., Teplova S. V. Development and application of new GTE parts hardening methods based on the plastic deformation of the surface layers, Vestn. Dvigatel., 2006, no. 2, pp. 8-16 (In Russian).
5. Dai K., Shaw L. Analysis of fatigue resistance improvements via surface severe plastic deformation, Intern. J. Fatigue, 2008, vol.30, no. 8, pp. 1398-1408. https://doi.org/ 10.1016/j.ijfatigue.2007.10.010.
6. James M. N., Hughes D. J., Chen Z., et al. Residual stresses and fatigue performance, Eng. Failure Anal., 2007, vol. 14, no. 2, pp. 384-395. https://doi.org/10.1016/j.engfailanal. 2006.02.011.
7. Majzoobi G. H., Azadikhah K., Nemati J. The effect of deep rolling and shot peening on fretting fatigue resistance of Aluminum-7075-T6, Mater. Sci. Eng. A, 2009, vol. 516, no. 1-2, pp. 235-247. https://doi.org/10.1016/j.msea.2009.03.020.
8. Soady K. A. Life assessment methodologies incorporating shot peening process effects: mechanistic consideration of residual stresses and strain hardening. 1. Effeact of shot peening on fatigue resistance, Mater. Sci. Technol., 2013, vol. 29, no. 6, pp. 673-651. https://doi.org/ 10.1179/1743284713Y.0000000222.
9. 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, Mater. Design., 2012, vol.35, pp. 741-748. https://doi.org/10.1016/j.matdes.2011.09.055.
10. Pavlov V. F., Kirpichev V. A., Vakuluk V. S. Prognozirovanie soprotivleniia ustalosti poverkhnostno uprochnennykh detalei po ostatochnym napriazheniiam [Prediction of Fatigue Resistance of Surface Reinforced Parts by Residual Stresses]. Samara, Samara Sci. Center of RAS, 2012, 125 pp. (In Russian)
11. Radchenko V. P., Pavlov V. Ph., Saushkin M. N. Investigation of surface plastic hardening anisotropy influence on residual stresses distribution in hollow and solid cylindrical specimens, PNRPU Mechanics Bulletin, 2015, no. 1, pp. 130-147 (In Russian). https://doi. org/10.15593/perm.mech/2015.1.09.
12. Radchenko V. P., Saushkin M. N., Bochkova T. I. Mathematical modeling and experimental study of forming and relaxation of residual stresses in plane samples made of EP742 alloy after ultrasonic hardening under high-temperature creep conditions, PNRPU Mechanics Bulletin, 2018, no. 3-4, pp. 88-98. https://doi.Org/10.15593/perm.mech/eng.2018.3.09.
13. Chen H., Wang S., Lu S., et al. Simulation and experimental validation of residual stress and surface roughness of high manganese steel after shot peening, Procedia CIRP, 2018, vol.71, pp. 227-231. https://doi.Org/10.1016/j.procir.2018.05.066.
14. Isa M. R., Sulaiman S. N., Zaroog O. S. Experimental and simulation method of introducing compressive residual stress in ASTM A516 grade 70 steel, Key Eng. Mater., 2019, vol. 803, pp. 27-31. https://doi.org/10.4028/www.scientific.net/KEM.803.27.
15. Kiselev I. A., Zhukov N. A., Vasilyev B. E., Selivanov A. N. Modeling of residual stresses when calculating strength of lock joint elements. Part 1. Modeling of the shot peening process, Proceedings of Higher Educational Institutions. Machine Building, 2018, no. 11, pp. 49-59 (In Russian). https://doi.org/10.18698/0536-1044-2018-11-49-59.
16. Meguid S. A., Maicic L. A. Finite element odeling of shot peening residual stress relaxation in turbine disk assemblies, J. Eng. Mater. Technol., 2015, vol. 137, no. 3, 031003. https:// doi.org/10.1115/1.4030066.
17. Gallitelli D., Boyer V., Gelineau M., et al. Simulation of shot peening: From process parameters to residual tress fields in a structure, Comptes Rendus Mécanique, 2016, vol. 344, no. 4-5, pp. 355-374. https://doi.org/10.1016Zj.crme.2016.02.006.
18. Zimmerman M., Klemenz M., Schulze V. Literature review on shot peening simulation, Int. J. Comput. Mater. Sci. Surf. Eng., 2010, vol.3, no. 4, pp. 289-310. https://doi.org/ 10.1504/ijcmsse.2010.036218.
19. Purohil R., Verma C. S., Rana R. S., et al. Simulation of shot peening process, Materials Today: Proceedings, 2017, vol.4, no. 2 A, pp. 1244-1251. https://doi.org/10.1016/j. matpr.2017.01.144.
20. Jebahi M., Gakwaya A., Lévesque J., et al. Robust methodology to simulate real shot peening process using discrete-continuum coupling method, Int. J. Mech. Sci., 2016, vol. 107, pp. 21-33. https://doi.org/10.1016/j.ijmecsci.2016.01.005.
21. 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, J. Appl. Mech. Tech. Phys., 2009, vol.50, no. 6, pp. 989-997. https://doi.org/10.1007/ s10808-009-0133-8.
22. 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, J. Appl. Mech. Tech. Phys., 2015, vol.56, no. 2, pp. 313-320. https://doi.org/10.1134/S0021894415020170.
23. Radchenko V. P., Tsvetkov V. V., Derevyanka E. E. Relaxation of residual stresses in a surface-hardened cylinder under creep conditions and rigid restrictions on linear and angular deformations, Mech. Solids., 2020, vol.55, no. 6, pp. 898-906. https://doi.org/10.3103/ S0025654420660024.
24. Radchenko V. P., Tsvetkov V. V., Saushkin M. N. Residual stress relaxation in a hardened cylinder under creep, loaded by an axial force, torque and internal pressure, J. Appl. Mech. Tech. Phys., 2020, vol.61, no. 4, pp. 583-592. https://doi.org/10.1134/ S0021894420040124.
25. Radchenko V. P., Saushkin M. N., Tsvetkov V. V. Effect of thermal exposure on the residual stress relaxation in a hardened cylindrical sample under creep conditions, J. Appl. Mech. Tech. Phys., 2016, vol.57, no. 3, pp. 559-568. https://doi.org/10.1134/ S00218944160302022.
26. Saushkin M. N., Prosvirkina E. A. Relaxation of residual stresses in a surface-hardened layer of a solid rotating cylinder under creep conditions, In: Proceedings of the Third All-Russian Scientific Conference (29-31 May 2006). Part 1, Matem. Mod. Kraev. Zadachi. Samara, Samara State Technical Univ., 2006, pp. 192-199 (In Russian).
27. Radchenko V. P., Saushkin M. N. Polzuchest' i relaksatsiia ostatochnykh napriazhenii v up-rochnennykh konstruktsiiakh [Creep and Relaxation of Residual Stresses in Hardened Structures]. Moscow, Mashinostroenie-1, 2005, 226 pp. (In Russian)
28. Birger I. A., Shorr B. F., Iosilevich G. B. Raschet na prochnost' detalei mashin [Calculation of the Strength of Machine Parts]. Moscow, Mashinostroenie, 1979, 702 pp. (In Russian)
29. Samarin Yu. P. Uravneniia sostoianiia materialov so slozhnymi reologicheskimi svoistvami [Equation of State of Materials with Complex Rheological Properties]. Kuibyshev, Kuibyshev State Univ., 1979, 84 pp. (In Russian)
30. Radchenko V. P., Eremin Yu. A. Reologicheskoe deformirovanie i razrushenie materialov i elementov konstruktsii [Rheological Deformation and Fracture of Materials and Structural Elements]. Moscow, Mashinstroenie-1, 2004, 265 pp. (In Russian)
31. Radchenko V. P., Berbasova T. I., Shishkin D. M. Relaxation of residual stresses in a surface-hardened prismatic sample subjected to biaxial loading under creep conditions, J. Appl. Mech. Tech. Phys., 2021, vol.62, no. 5, pp. 861-869. https://doi.org/ 10.1134/S0021894421050187.
