Проблемы машиностроения и надежности машин, 2019, № 5, стр. 89-97

ОСОБЕННОСТИ ТЕРМОРЕЛАКСАЦИИ ОСТАТОЧНЫХ НАПРЯЖЕНИЙ, ВОЗНИКАЮЩИХ ПРИ ЛАЗЕРНО-УДАРНО-ВОЛНОВОЙ ОБРАБОТКЕ ЖАРОПРОЧНЫХ МАТЕРИАЛОВ

Г. Ж. Сахвадзе 1*, Р. З. Кавтарадзе 2, Т. М. Натриашвили 3, Г. Г. Сахвадзе 14

1 Институт машиноведения им. А.А. Благонравова РАН
г. Москва, Россия

2 Московский государственный технический университет им. Н.Э. Баумана
г. Москва, Россия

3 Институт механики машин им. Р. Двали
г. Тбилиси, Грузия

4 Московский государственный университет им. М.В. Ломоносова
г. Москва, Россия

* E-mail: sakhvadze@mail.ru

Поступила в редакцию 25.03.2019
Принята к публикации 25.06.2019

Полный текст (PDF)

Аннотация

В настоящей статье исследуются особенности терморелаксации сжимающих остаточных напряжений, сгенерированных при лазерно-ударно-волновой обработке, высоколегированной жаропрочной стали 37Х12Н8Г8МФБ. Проведено конечноэлементное моделирование процесса терморелаксации сгенерированных сжимающих остаточных напряжений. Получены закономерности влияния термических воздействий на характер перераспределения сжимающих остаточных напряжений в температурном диапазоне от 200 до 650°С. Сравнительный анализ показал, что полученные в результате конечноэлементного моделирования данные хорошо коррелируются с известными в литературе экспериментальными данными.

Ключевые слова: лазерно-ударно-волновая обработка, терморелаксация, остаточные напряжения, конечноэлементное моделирование

Введение. Все известные способы механической обработки поверхностей, такие как глубокая прокатка, лазерно-ударно-волновая обработка (ЛУВО), дробеструйная обработка и др. оказывают положительное влияние в повышении качества обрабатываемых поверхностей [1]. Это достигается благодаря сгенерированным при этих процессах полям сжимающих остаточных напряжений (СОН), которые блокируют нежелательные процессы появления и распространения усталостных трещин, тем самым повышая их усталостные свойства. Для конструкций, которые во время эксплуатации испытывают воздействия повышенных механических и/или температурных нагрузок, эти полезные эффекты имеют настоящую ценность только в том случае, если сгенерированные приповерхностные СОН остаются стабильными и, если они не полностью релаксируют со временем. Известно, что по сравнению с традиционными способами упрочнения поверхностей, именно технология ЛУВО создает глубокие и термостабильные слои СОН для целого ряда материалов [27].

Высоколегированные жаропрочные стали 37Х12Н8Г8МФБ, исследуемые в настоящей статье, в основном используется в ответственных компонентах авиационных двигателей, которые работают в условиях высокой температуры и всегда подвержены циклическим нагрузкам. Обычно повреждения компонентов авиационных двигателей происходят из-за возникновения микротрещин на поверхности. Сжимающие остаточные напряжения, сгенерированные при ЛУВО, могут существенно улучшить усталостные характеристики металлических материалов, и они будут сопротивляться возникновению и распространению трещин. Однако эти СОН могут быстро релаксироваться при повышенных температурах. В литературе очень мало данных о влиянии ЛУВО на механические свойства стали 37Х12Н8Г8МФБ при повышенных температурах. Поэтому очень важно изучение термической стабильности СОН, возникающих в исследуемых материалах при обработке по технологии ЛУВО. В статье исследуется процесс терморелаксации СОН при повышенных температурах методом конечных элементов.

Исследуемый материал. С целью проведения сравнительного анализа и проверки адекватности разработанной численной модели расчета, конечноэлементным моделированием (КЭМ) было исследовано для материала, который использовался в экспериментах [8] – высоколегированной жаропрочной стали 37Х12Н8Г8МФБ (ее зарубежный аналог – Iron GH2036 alloy). Физико-механические характеристики стали приведены в табл. 1. Для моделирования технологии ЛУВО и процесса терморелаксации исследовался образец размером (30 × 30 × 5) мм.

Таблица 1.

Физико-механические характеристики высоколегированной жаропрочной стали 37Х12Н8Г8МФБ

Характеристика Значение Размерность
Предел текучести 600 МПа
Предел прочности 1000 МПа
Модуль упругости 203 ГПа
Твердость по Бринеллю 269 МПа
Коэффициент температурного расширения 15.9 × 10–6 1/град
Коэффициент теплопроводности 17.14 Вт/(м град)
Теплоемкость 440 Дж/(кг град)

Схематизация технологии ЛУВО и возникновения СОН. Схематично общая методология технологии лазерно-ударно-волновой обработки и возникновения СОН показана на рис. 1.

Рис. 1.

Принципиальная схема технологии ЛУВО (а), схемы возникновения пластических деформаций в приповерхностной зоне во время ЛУВО (б) и ОН после ЛУВО (в): 1 – лазерный импульс; 2 и 3 – прозрачный и непрозрачный слои; 4 – ударная волна; 5 – исследуемый материал; 6 – плазма; 7 – зона воздействия лазерного импульса; 8 и 9 – сжимающие и растягивающие пластические деформации; 10 – релаксация деформируемого материала и возникновение СОН.

Обрабатываемая поверхность (рис. 1а, позиция 5) покрывается непрозрачным для лазерного излучения слоем с низкой температурой испарения: черной краской, металлической фольгой или лентой 3. Сверху находится прозрачный слой 2, роль которого наиболее часто выполняет вода или стекло. Энергия лазерного импульса 1 поглощается непрозрачным слоем, что приводит к его нагреву, испарению и формированию высокотемпературной плазмы, ограниченной с одной стороны поверхностью исследуемого материала, а с другой – прозрачным слоем, сдерживающим распространение плазмы. Вследствие ограниченности объема, давление резко возрастает до высоких значений и быстро распространяется в материал в виде ударной волны 4, которая способствует появлению в материале сжимающих остаточных напряжений. Если эти напряжения превосходят упругий предел Гюгонио, то материал деформируется пластически [2]. Более подробно протекающие при ЛУВО физико-механические процессы описаны в [4].

Схематизация численного моделирования. Предложенная в настоящей статье стратегия конечноэлементного моделирования определения терморелаксации СОН, сгенерированных при ЛУВО, состоит из трех шагов (рис. 2).

Рис. 2.

Схема моделирования технологии ЛУВО и определения терморелаксации остаточных напряжений.

Область КЭМ обведена пунктирными линиями: 1) вычисление собственных деформаций, вызванных с применением технологии ЛУВО путем решения явной задачи динамического анализа посредством пакета LS-DYNA; 2) далее решается неявная задача ANSYS, где полученные на шаге 1 собственные деформации вводятся в качестве входных данных, и выполняется решение задачи для определения полных СОН, вызванных технологией ЛУВО; 3) полученные на шаге 2 СОН являются входными данными для решения неявной сопряженной термомеханической задачи в ANSYS, где добавляется термонагружение по поверхности исследуемого материала и вычисляется терморелаксация СОН.

Заметим, что при моделировании процессов ЛУВО возникают большие вычислительные затраты. Чтобы существенно уменьшить их, используют известный метод собственных деформаций [9]. Более подробно особенности численного моделирования, дискредитации по сетке и по времени приведены в [10].

Определяющее соотношение. При ЛУВО из-за воздействия большого по величине и кратковременного по времени давления скорости деформации в исследуемом материале достигают величин 106 c–1 и более, поэтому при расчетах напряженно-деформированного состояния в качестве определяющего соотношения используется модель пластического деформирования материала Джонсона–Кука [11]

(1)
$\sigma = (A + B{{\bar {\varepsilon }}^{n}})(1 + C\dot {\bar {\varepsilon }}*)[1 - {{(T*)}^{m}}],$
где ${\sigma }$ – эквивалентное напряжение по Мизесу, ${\bar {\varepsilon }}$ – эквивалентная пластическая деформация, $\dot {\bar {\varepsilon }}{\kern 1pt} {\text{*}}$ = $\dot {\bar {\varepsilon }}{\text{/}}{{\dot {\bar {\varepsilon }}}_{0}}$ – безразмерная эквивалентная скорость деформации (${{{\dot {\bar {\varepsilon }}}}_{0}}$ = 1.0 с–1), T* – гомологическая температура, связь которой с абсолютной температурой Т определяется следующим образом
$T* = {\text{(}}T - {{T}_{0}}{\text{)/(}}{{T}_{m}} - {{T}_{0}}{\text{)}},$
где T0 – комнатная температура, Tm – температура плавления материала образца.

Уравнение (1) содержит пять материальных констант, которые определяются эмпирическим путем: A – статический предел текучести, B – модуль деформационного упрочнения, n – показатель степени в законе деформационного упрочнения, C – коэффициент скорости деформаций, m – показатель степени в законе температурного разупрочнения. Материальные константы для высоколегированной жаропрочной стали 37Х12Н8Г8МФБ приведены в табл. 2 [12]. Заметим, что модель Джонсона–Кука в аналогичных задачах апробирована и дает удовлетворительные результаты [1319].

Таблица 2.

Константы для модели Джонсона–Кука для высоколегированной жаропрочной стали 37Х12Н8Г8МФБ

Константа Значение Размерность
A 900 МПа
B 1200 МПа
C 0.0092
n 0.6
m 1.27

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

(2)
$\frac{{{{{\sigma }}^{{res}}}}}{{{\sigma }_{0}^{{res}}}} = \exp [ - {{(A{{t}_{a}})}^{m}}],$
где ${\sigma }_{{\text{0}}}^{{res}}$ являются начальными СОН до отжига, ${{{\sigma }}^{{res}}}$ – остаточные напряжения при заданном времени отжига ${{t}_{a}}$ и при заданной температуре отжига Ta; m – численный параметр, в основном зависящий от доминирующего механизма релаксации; А – это функция, зависящая от свойств материала и температуры, которую можно определить следующим образом [8]
(3)
$A = B\exp \left( { - \frac{{\Delta H}}{{k{{T}_{a}}}}} \right),$
где B – постоянная величина для данного материала (для исследуемого нами материала она равна 2.05 × 1012 мин–1 [8]); k – постоянная Больцмана (для исследуемого материала она равна 8.617343 × 10–5 эВ/К [8]); ΔH – энтальпия активации материала при релаксации; Ta – температура отжига.

Для анализа уравнение (2) удобнее представить в виде

(4)
$\ln \left[ {\ln \left( {\frac{{\sigma _{0}^{{res}}}}{{{{\sigma }^{{res}}}}}} \right)} \right] = m~\ln ~A + m\ln {{t}_{a}}.$

На рис. 3 в качестве примера показан график зависимости $\ln [\ln (\sigma _{0}^{{res}}{\text{/}}{{\sigma }^{{res}}})]$ от lnta при заданной температуре отжига Та.

Рис. 3.

Влияние времени отжига ta на терморелаксацию СОН при разных температурах отжига Та: 1Ta = = 200°С; 2Ta = 300°С; 3Ta = 400°С; 4Ta = 500°С; 5Ta = 600°С; 6Ta = 650°С.

Эта зависимость почти линейная, при этом угол наклона линии численно равен mi [8]. За истинным значением m берется среднее значение наклонов всех линий. После этого из формулы (3) легко можно определить энтальпию активации процесса релаксации ΔH. Согласно полученным данным, значения m равны 0.22, 0.31, 0.311, 0.412, 0.512 и 0.513 при температурах 200°C, 300°C, 400°C, 500°C, 600°С и 650°С, соответственно. Таким образом среднее значение m составляет 0.38, а значение ΔH составляет около 1.88 эВ. Зная m и ΔH, из уравнения (3) можно определить, например, время отжига ${{t}_{a}}$, при котором достигается заданный уровень СОН (${{{\sigma }}^{{res}}}$).

Полученные результаты и их анализ. На рис. 4 показаны распределения поверхностных и приповерхностных СОН при разных температурах отжига Ta. Максимальные СОН наблюдаются в поверхностном слое, которые составляют около 510 МПа. Глубина упрочненного слоя со значительными СОН составляет примерно 0.62 мм. Исследования показывают, что СОН уменьшаются по мере удаления от поверхности. Кроме того, на определенной глубине (это – примерно 0.58 мм) появляются нежелательные растягивавшие остаточные напряжения (РОН).

Рис. 4.

Распределения СОН при разных температурах отжига Ta: (а) по поверхности (по оси x); (б) в приповерхностной области (по глубине – по оси z): 1 – после ЛУВО (начальное распределение СОН – до приложения термической нагрузки); 2 – Ta = 200°С; 3Ta = 300°С; 4Ta = 400°С; 5Ta = 500°С; 6Ta = 600°С; 7Ta = 650°С.

Из количественного и качественного анализа данных, приведенных на рис. 5 следует, что релаксация СОН активно происходит при различных температурах в диапазоне от 200 до 650°С. Можно сделать важное заключение: они полностью не обнуляются даже при температуре 650°С. Максимальные СОН уменьшаются с 510 МПа (сразу после ЛУВО, до приложения термической нагрузки) до 401 МПа (при температуре отжига 200°С); до 370 МПа (при температуре 300°С); до 330 МПа (при температуре 400°С); до 270 МПа (при температуре 500°С); до 170 МПа (при температуре 600°С) и до 145 МПа (при температуре 950°С). Время отжига в нашем случае – всюду 60 мин. Эти уменьшения СОН в процентах составляют 21.4%, 27.5%, 35.3%, 47.2%, 66.6% и 71.5%, соответственно.

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

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

На рис. 5 показано влияние температуры отжига на характер перераспределения поверхностных и приповерхностных СОН, вызванных ЛУВО. Видно, что остаточные напряжения постепенно релаксируют с увеличением температуры отжига, и на определенной глубине (примерно 0.5 мм) становятся почти нулевыми.

Что касается учета фактора времени, можно сделать вывод, что СОН релаксируют в основном в течение первых 30 минут термообработки. Далее они имеют тенденцию стабилизироваться и терморелаксация практически прекращается.

Сравнительный анализ. Проверка адекватности разработанной численной модели

На рис. 6 показаны экспериментальные (кривая 2, взятая из [8]) и расчетные (кривая 1) зависимости распределения сжимающего остаточного напряжения от температуры отжига Та при заданной величине времени отжига ta (в нашем случае ta = 60 мин). Как видно (кривая 1) СОН плавно меняются с повышением температуры от 200 до 450°С (зона I), при этом наклон кривой (${{{\sigma }}^{{res}}}$Та) значительно увеличивается после 450°С, но после 600°С снова становятся гладким, что указывает на тот факт, что температура отжига в интервале от 450°C до 600°C (зона II) оказывает сильное влияние на процесс терморелаксации СОН исследуемого материала по сравнению с более низкой или более высокой температурами. Проведенные исследования показывают, что предел текучести уменьшается более заметно с повышением температуры именно в диапазоне от 450 до 600°C, т.е. процесс терморелаксации СОН в этом температурном интервале протекает более интенсивно по сравнению с другими интервалами. Сравнение с экспериментальными данными дает удовлетворительные результаты. Отсюда можно сделать вывод о том, что эффект влияния температуры и времени отжига на терморелаксацию СОН, сгенерированных технологией ЛУВО, может быть удовлетворительно описан с использованием модели Зинера–Верта–Аврами, использованного в нашем КЭМ.

Рис. 5.

Характер распределения СОН в зависимости от температуры отжига Ta на поверхности и в приповерхностной области: 1 – на поверхности; 2 – на глубине 0.1 мм; 3 – на глубине 0.2 мм; 4 – на глубине 0.35 мм; 5 – на глубине 0.5 мм.

Рис. 6.

Влияние температуры отжига Та (при времени отжига 60 мин) на характер распределения сжимающих остаточных напряжений ${{{\sigma }}^{{res}}}$, сгенерированных технологией ЛУВО: 1 – результаты конечноэлементного моделирования, 2 – экспериментальные результаты из [8], I – зона плавного изменения СОН, II – зона резкого изменения СОН.

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

Выводы. Технология лазерно-ударно-волновой обработки имеет уникальные преимущества в задаче улучшения механических свойств и усталостных характеристик материалов за счет создания сильных полей сжимающих остаточных напряжений в приповерхностной области. В статье описано исследование влияния термической обработки (в условиях изотермического отжига) на терморелаксацию сжимающих остаточных напряжений для высоколегированной жаропрочной стали 37Х12Н8Г8МФБ в температурном диапазоне от 200 до 650°С на основе имитационного конечноэлементного моделирования.

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

Список литературы

  1. Nikitin I., Besel M. Residual stress relaxation of deep-rolled austenitic steel // Scr. Mater. 2008. № 58. P. 239.

  2. Peyre P., Fabbro R., Merrien P. Laser shock processing of aluminium alloys // Application to high cycle fatigue behavior. Mater Sci. Eng. 1996. № 210. P. 102.

  3. Torres M., Voorwald H. An evaluation of shot peening, residual stress and stress relaxation on the fatigue life of AISI 4340 steel // Int. J. Fatigue. 2002. № 24. P. 877.

  4. Kwofie S. Description and simulation of cyclic stress–strain response during residual stress relaxation under cyclic load // Procedia Eng. 2011. № 10. P. 293.

  5. Xu S., Huang S., Meng X., Sheng J., Zhang H., Zhou J. Thermal evolution of residual stress in IN718 alloy subjected to laser peening // Optics and Lasers in Engineering. 2017. № 94. P. 70.

  6. Nikitin I., Scholtes B., Maier H.J. High temperature fatigue behavior and residual stress stability of laser-shock peened and deep rolled austenitic steel AISI 304 // Scr. Mater. 2004. № 50. P. 1345.

  7. Nalla P., Altenberger I., Noster U. On the influence of mechanical surface treatments-deep rolling and laser shock peening-on the fatigue behavior of Ti–6Al–4V at ambient and elevated temperatures // Mater. Sci. Eng. A. 2003. № 355. P. 216.

  8. Ren X.D., Zhou W.F., Ren Y.P., Xu S.D., Liu F.F., Yuan S.Q., Ren N.F., Huang J.J. Dislocation evolution and properties enhancement of GH2036 by laser shock processing: Dislocation dynamics simulation and experiment // Materials Science & Engineering A. 2016. V. 654. P. 184.

  9. Korsunsky A.M. Residual elastic strain due to laser shock peening: modelling by eigenstrain distribution // J. Strain Anal. Eng. Des. 2006. V. 41 (3). P 195.

  10. Zhong Z., Sagar B., Gokul R. Thermal relaxation of residual stress in laser shock peened Ti–6Al–4V alloy // Surf. Coat Technol. 2012. № 206. P. 4619.

  11. Zhong Z., Amrinder S.G., Dong Q. A finite element study of thermal relaxation of residual stress in laser shock peened IN718 superalloy // Int. J. Impact Eng. 2011. № 38. P. 590.

  12. Liao Y.L., Suslov S., Ye C. The mechanisms of thermal engineered laser shock peening for enhanced fatigue performance // Acta Mater. 2012. № 60. P. 4997.

  13. Сахвадзе Г.Ж. Особенности конечноэлементного моделирования остаточных напряжений, возникающих в материале при лазерно-ударно-волновой обработке, с использованием метода собственных деформаций // Проблемы машиностроения и надежности машин. 2018. № 4. С. 87.

  14. Сахвадзе Г.Ж. Конечноэлементное моделирование технологии многократной лазерно-ударно-волновой обработки материалов с использованием метода собственных деформаций // Проблемы машиностроения и надежности машин. 2018. № 5. С. 91.

  15. Сахвадзе Г.Ж., Гаврилина Л.В., Бармина О.В., Киквидзе О.Г. Повышение микротвердости материалов технологией двухсторонней лазерно-ударно-волновой обработки // Вестник научно-технического развития: Интернет-журнал № 7(107). 2016. С. 29. http://www.vntr.ru/nomera/2016-7107/

  16. Сахвадзе Г.Ж., Пугачев М.С., Киквидзе О.Г. Технология двухсторонней лазерно-ударно-волновой обработки материалов // Вестник машиностроения. 2016. № 10. С. 71.

  17. Сахвадзе Г.Ж., Кавтарадзе Р.З., Никабадзе М.У. Использование метода собственных деформаций при моделировании технологий однократной и многократной лазерно-ударно-волновой обработки материалов // Вестник машиностроения. 2018. № 7. С. 41.

  18. Sakhvadze G., Shokhin A., Kikvidze O. Residual Stresses distribution in Ti–6Al–4V titanium alloys during Laser Shock Processing // Vibroengineering PROCEDIA, October 2016. V. 8. P. 422.

  19. Vlastimir N., Miloš M., Dalibor P., Dejan J., Milan S. Parameters forecasting of laser welding by the artificial intelligence techniques // Facta Universitatis, Series: Mechanical Engineering. 2018. V. 16. № 2. P. 193. https://doi.org/10.22190/FUME180526025N

Дополнительные материалы отсутствуют.