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

Применение дифракционных оптических элементов – фокусаторов излучения [1, 2] при лазерной обработке материалов предоставляет возможность создать требуемый комплекс свойств конструкционных материалов в зоне термического влияния. В этом случае форма лазерного пятна и распределение плотности мощности излучения рассматриваются в качестве основных параметров режима обработки, выбор параметров лазерного источника и оптической системы проводится в соответствии с результатами решения обратной задачи теплопроводности. В работах [3, 4] разработана методика расчета пространственного распределения мощности лазерного излучения для формирования требуемого энергетического воздействия на технологические объекты, применение которой позволяет создать более равномерное температурное поле по длине движущегося полосового источника.

Целью данной работы являются экспериментальные исследования температурных полей в конструкционном материале – хромоникельмолибденовой стали 40ХНМА при воздействии лазерных потоков, сформированных фокусаторами излучения.

Применение инфракрасного термометра «Кельвин» в составе системы бесконтактной диагностики температур

При проведении исследований для формирования излучения СО 2 -технологической установки ЛК-1300 «Хебр-1А» мощностью Q =10³ Вт применялся отражающий дифракционный оптический элемент [5]:

где Z ( U , V ) - высота рельефа в точке ( U , V ) оптического элемента; U , V – координаты точки отражающего оптического элемента в системе координат с осью OU , направленной противоположно проекции падающего луча и имеющей начало в центре расположенной в области G : U ² cos 0 + V ² = R ² отражающей поверхности; R – максимальный радиус фокусируемого пучка; L – длина фокального отрезка; f – фокусное расстояние оптического элемента; 0 - угол между оптической осью излучения и нормалью к плоскости оптического элемента; mod( h ) _{m X /(2 cos 0 )} - функция, равная наименьшему положительному остатку от деления h на m X/(2 cos 0 );

X - длина волны излучения; m - целое число; r -параметр фокусируемого пучка Гауссовского распределения плотности излучения, в котором интенсивность на расстоянии r от центра уменьшается в e раз по сравнению с интенсивностью в центре пучка; ( 9 ) - значение двойного интеграла вероятности в области G для фокусируемого пучка гаус-

совского распределения плотности излучения:

(9)* = 1JJexp

^П G

(U2 + V2)

d U d V .

Z ( U , V ) = mod-^ _m

V , /           m /(2 cos 0 )

( 1 ( U² cos ²0+ V² cos el      4 f

X X

U cos 0     V R ² - U ²

L J ^dU J

Пf (9)

Фокусатор имел следующие параметры: f =1 м; L ₀=12-10 " ³ м; R =25^10 " ³ м; r = 0,7 R ; 0 =45°.

Для определения температуры в зоне лазерного воздействия в диапазоне значений T =(600....1600) K использовался инфракрасный термометр «Кельвин», имеющий предел допустимой абсолютной погрешности А в диапазоне температуры окружающей среды (293±5) K: А = ( 1 + 0,01 | т ; + 273| ) K, где Т_Т -показания прибора (°С). Для дополнительной калибровки прибора установлена зависимость между регистрируемым потоком теплового излучения и глубиной закалки центральной области зоны термического влияния.

На рис. 1 а представлены расчетные зависимости глубины закалки технологического объекта из стали 40ХНМА с поглощающим покрытием МЦС-510 на основе метилцеллюлозы и силиката натрия при различной скорости движения полосовых энергетических источников:

q ( x , y ) = q₀ ( - 2,4568 ■ 10 ^-7 X ¹⁰ + 2,5128 ■ 10 ^-5 X ⁸ -

• - 7,7806 ■ 10 ^-4 X ⁶ + 6,8237 ■Ю ^-3 X ⁴ -

  • - 1,7760 ■Ю ^-2 X ² + 1,0142) ■                         (3)

■ ( - 1,0086 y ⁶ + 2,7886 У ⁴ - 2,7350 У ² +

+ 0.9967) [Вт/м ² ]

при q ₀=8,8422^10⁷ Вт/м²;

X € [ - 5,5 ■Ю 3 м; 5,5 ■Ю - ³ м]; y е [ - 10 - ³ м;10 ^{- 3} м].

q ( X , y ) = q ₀ (1,1513 ■ 10 ^-7 X ¹⁰ - 1,0036 ■ 10 ^-5 X ⁸ +

+ 2,4962 ■Ю^-4 X ⁶ - 2,0921 ■ 10 ^-3 X ⁴ +

+ 8,2705 ■Ю^-3 X ² + 0,9969) x                       (4)

x ( - 1,3516 У ⁶ + 3,3895 У ⁴ - 3,0153 У²

+ + 0,9966) [Вт/м ² ]

при q ₀=7,9698-10⁷ Вт/м²;

X е [ - 6,8■Ю ^{- 3} м; 6,8■Ю ^{- 3} м]; у е [ - 10 ^{- 3} м; 10 ^{- 3} м], где q 0 - плотность мощности в центре энергетического источника; x , у - пространственные координаты; x = x / 10 - ³ м , y = у /10 ³ м - безразмерные координаты. Формулы (3) и (4) получены в результате аппроксимации расчетного распределения плотности мощности в фокальной плоскости полиномами десятой степени [4].

Зависимость (3) описывает распределение плотности мощности в фокальной плоскости оптического элемента (1) при значении параметров фокусируемого пучка R_Ф = г_Ф =18^10^-3 м и мощности излучения Q =950 Вт, а зависимость (4) - при R ф = Г ф =22^10^-3 м.

При определении глубины закалки hz рассматривалась нелинейная задача теплопроводности при нагреве технологического объекта полосовым поверхностным энергетическим источником. Расчет проводится для технологического объекта, который перемещался с постоянной линейной скоростью и в положительном направлении оси Oy в системе координат (x, у, z), помещенной в центр энергетического источника мощностью Q . Принималось допущение: энергетический источник с заданными геометрией и пространственным распределением мощности является поверхностным и занимает область ф на поверхности исследуемого объекта. На остальной части поверхности задаются нелинейные условия теплообмена с окружающей средой.

После воздействия на образец лазерного излучения изготавливались микрошлифы поперечного сечения зоны термического влияния и проводились измерения глубины закаленного слоя. Характерные диаграммы регистрируемого сигнала для образца из стали 40ХНМА с поглощающим покрытием МЦС-510, а также результаты измерений глубины закаленных слоев образцов приведены на рис. 1 б. На рис. 1 обозначено: п, - относительные единицы ед показаний прибора; hz - глубина закалки.

Рис. 1. Расчетные зависимости глубины закалки h _z технологического объекта из стали 40ХНМА с поглощающим покрытием от температуры поверхности T при воздействии поверхностных энергетических источников (3), (4) (а); результаты измерения глубины зоны закалки h _z образцов из стали 40ХНМА с поглощающим покрытием МЦС-510 после проведения лазерного воздействия и интенсивности сигналов инфракрасного термометра n _ед , полученные на этих образцах (б). Скорость движения энергетических источников υ , 10^-3 м/с: кривая 1 – 1,3; кривая 2 – 0,7

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

Результаты экспериментальных исследований

В результате расчета и эксперимента определены: изотермы температурного поля и распределение температуры T по ширине зоны термического влияния на поверхности деталей из хромоникельмолибденовой стали 40ХНМА, имеющие место при движении полосовых энергетических источников мощностью Q =0,95∙103 Вт со скоростью υ =1,1∙10-2 м/с, представленные на рис. 2, 3.

Рис. 2. Изотермы температурного поля на поверхности детали из стали 40ХНМА при проведении лазерного термоупрочнения. Скорость движения полосовых энергетических источников с различным распределением мощности υ =1,1∙10^-2м/с. Мощность излучения СО₂-лазерной установки Q =950 Вт. пунктирная линия – результаты расчета, сплошная – эксперимент

Рис. 3. Распределение температуры T по ширине зоны термического влияния на поверхности стали 40ХНМА при проведении лазерного термоупрочнения полосовыми энергетическими источниками с различным распределением мощности. Скорость движения энергетических источников υ =1,1∙10^-2м/с. Мощность излучения СО₂-лазерной установки Q =950 Вт. пунктирная линия – результаты расчета, точки – эксперимент

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

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

Скорости нагрева и охлаждения материала на поверхности деталей из хромоникельмолибденовой стали 40ХНМА, имеющие место при движении полосовых энергетических источников мощностью

Q =0,9540 3 Вт со скоростью и =1,140 "2 м/с, приведены на рис. 4.

Рис. 4. Термические циклы (а), а также скорости нагрева и н и охлаждения и _охл (б) точек на поверхности технологического объекта из стали 40ХНМА.

x , 10 - м : кривая 1 - 0; кривая 2 - 2; кривая 3 - 4;

кривая 4 – 6; кривая 5 – 8. t _пр – момент прохождения центральной части лазерного пятна над расчетным сечением

При воздействии излучения, имеющего параметры фокусируемого пучка R_Ф = r_Ф =18∙10^-3 м, и мощность, равную Q =950 Вт, на движущийся со скоростью и =1,1•10^-2м/с технологический объект из стали 40ХНМА, величина максимальной температуры на поверхности в центре зоны термического влияния равна max |^ T ( x = 0, y , z )| z = h ] =1450 K. Ширина изотермы T =1100 K составляет b =6,5 - 10^-3 м.

1 = 1100 K

Изменение параметров фокусируемого пучка излучения до значений R_Ф = r_Ф =22∙10^-3 м приводит к уменьшению в 1,12 раза величины максимальной температуры в центре зоны термического влияния, которая составляет max |^ T ( x = 0, y , z )| z = h ] =1300 K.

Ширина изотермы T =1100 K увеличивается в 1,4 раза и равна b_T = _{1100 K}=9,2 - 10^-3 м. Результаты экспериментальных исследований температурных полей в обрабатываемом материале при воздействии движущихся распределенных поверхностных источников энергии удовлетворительно (погрешность 8 _T <  10^15%) коррелирует с результатами расчетов. Скорость охлаждения в диапазоне значений температур от T =1100 K до T =940 K составляет и _охл « 10'^2 - 10 ' K/с.

Заключение

Приведенные результаты экспериментальных исследований температурных полей в обрабатываемом материале – хромоникельмолибденовой стали 40ХНМА при воздействии движущихся распределенных поверхностных источников энергии удовлетворительно (погрешность 8 _T <  10^15%) коррелирует с результатами расчетов. Воздействие на объекты полосовых движущихся энергетических источников предоставляет возможность обеспечить выравнивание значений максимальных температур в поперечных сечениях зоны термического влияния в ходе термических циклов нагрева и охлаждения, а увеличение плотности мощности лазерного излучения по краям полосы компенсирует повышенный теплоотвод с периферийных участков. Дополнительное перераспределение мощности излучения к краям лазерного пятна достигается путем повышения доли энергии, формируемой периферийными зонами фокусатора, например, увеличением радиуса фокусируемого луча или использованием кольцевого распределения плотности потока. Достигается выравнивание значений максимальных температур вдоль линии перемещения энергетического источника в центре и на периферии зоны термического влияния с увеличением ширины требуемых изотерм максимальных температур.