Расчет температурных полей в конструкционных материалах в программном комплексе Star-CD Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

Физика

УДК 621.791.72

РАСЧЕТ ТЕМПЕРАТУРНЫХ ПОЛЕЙ В КОНСТРУКЦИОННЫХ МАТЕРИАЛАХ В ПРОГРАММНОМ КОМПЛЕКСЕ STAR-CD

© 2008 С.П. Мурзин, А.В. Меженин, Е.Л. Осетров

Самарский государственный аэрокосмический университет

Проведены расчеты температурных полей в программном комплексе STAR-CD для технологического объекта из конструкционного материала - листового низколегированного титанового сплава ОТ4-1. Проведенные исследования показали, что результаты расчета температурных полей в программном комплексе STAR-CD в листовом низколегированном титановом сплаве ОТ4-1 при воздействии движущихся распределенных поверхностных источников энергии удовлетворительно (погрешность не превышает 15...20%) коррелируют с результатами экспериментальных исследований. Это позволяет сделать вывод о применимости программного комплекса STAR-CD для расчета температурных полей в конструкционных материалах.

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

При обработке изделий лазерным излучением характерны большие градиенты температур, смещение точек структурных и фа-

зовых переходов, изменение в процессе обработки оптических и теплофизических свойств материалов. Принятая для расчетов температурных полей схема нагрева представлена на рис. 1. На рис. 1 обозначено: Н1, Н2 - ширина и толщина исследуемого объекта соответственно.

При расчете температурных полей целесообразно использовать программный комплекс STAR-CD, который рассчитан на применение в рамках интегрированной САЕ среды и предназначается для исследований, проектирования, разработки и оптимизации задач, связанных с моделированием течения жидкостей и газов, термодинамических явлений, а также для расчета температурных процессов. Для численного решения задач в STAR-CD реализован метод контрольного

Рис. 1. Схема расчета температурных полей при лазерной обработке в декартовой системе координат

(х, y, z)

767

Известия Самарского научного центра Российской академии наук, т. 10, №3, 2008

объема. Одним из важных свойств этого метода, подробно рассмотренного в работах [3, 4], является точное интегральное сохранение величины энергии на всей расчетной области при любом минимальном числе узловых точек. Это позволяет значительно ускорить расчеты при экономии ресурсов памяти персонального компьютера. Для расчета температурных полей в программном комплексе STAR-CD используется уравнение сохранения энтальпии:

f+| (pUh+^ )+| p^)+

d / \ dP dP dP

+ — (pWh + Fhz ) = —+ U — + V— +

dz

dt

dx dy

dP

+ W— + t dz

dV

+ T21 + T22

dx

dU

■ + t

dW

+ T32 - + T33

dx dV dy dW

12

+ t

dU

dy

dV

dz

■ + t

13

d_U_

dz

+

■ + t

dW (1)

— + dx

dy

dz

+ S„

где p - плотность; h - термическая энтальпия; t - время; U, V, W - значения скоростей в направлении x,y, z соответственно; Fh,x,Fhy,Fhz - поток диффузной энергии в направлении x,y,z соответственно; P - давление; t - тензор напряжений; Sh - источник энергии (величина, учитывающая подвод и отвод энергии).

Для схемы, представленной на рис. 1, уравнение сохранения энтальпии (1) примет вид:

dph d ( ) , dFh,x , dFh,z п

ИГ+н(рп+Fhy )+“аГ+ИТ =0- (2)

zr , dT dT dT

hx = k lx;FKy = k ~dy;F z = k^, (3)

dz

где k - коэффициент теплопроводности; t - температура.

h = cpT - c0To, (4)

где t - температура; сp - удельная теплоемкость при постоянном давлении для температуры t ; c°P - удельная теплоемкость при постоянном давлении для начальной температуры То;

Подставив значения выражений (3) и (4) в уравнение (2) получим уравнение математической модели тепловой задачи при лазерной обработке в неподвижной системе координат:

_3_

dx

k(T) — dx

+i lk(T)

dT_

dy

+1 [k(T)

dT_

dz

/т i i dT dT

- c'(T)-'sT+v ay ' =0

(5)

Проведем расчет температурного поля в пакете STAR-CD для технологического объекта из конструкционного материала ? листового низколегированного титанового сплава ОТ4-1 (160Ч100Ч2 мм), который перемещается с постоянной скоростью v = 0,6»10-2 м/с в положительном направлении оси Oy в системе координат (x-, y, z), помещенной в центр энергетического источника, описываемого функцией:

.,—10

q( x , y) = q0 (-2,4568 • 10-7 x +

C 8 л 6

+ 2,5128 • 10"5 x - 7,7806 • 10"4x +

o—4 o—2

+ 6,8237 • 10'3 x -1,7760 •10"2 x +

— 6 —4 (6)

+1,0142) • (-1,0086y + 2,7886y -

— 2 о

- 2,7350y + 0.9967) [Вт/м2]

при q0 =4,1884»107 Вт/м2;

x e [-5,5 40-3м; 5,5 -10-3м]; y e [-10-3м; 10 -3м],

а также энергетического источника, описываемого функцией:

,—10

q(x, y) = q0 (1,1513 • 10- x -

C 8 л 6

-1,0036 • 10"5 x + 2,4962 • 10'4x -

0—4 o—2

- 2,0921 • 10"3 x + 8,2705 • 10'3 x +

— 6 —4 (7)

+ 0,9969) • (-1,3516y + 3,3895y -

— 2 о

-3,0153y + 0,9966) [Вт/м2]

при q0 =3,7752-107 Вт/м2; x e [-6,8 40 3м; 6,8 40 3м]; y e [-10-3м; 10 -3м].

Зависимости коэффициента теплопроводности k(T) и удельной объемной теплоемкости cV(T) от температуры описывают-

768

Физика

ся квадратичными полиномами [5]:

k(T) = 10,23 + 0,0063233 T + 1,6775 • 10~6T2,(8)

cV(T) = 414,22 + 0,30039T + 0,00011797 T 2.(9)

Коэффициент теплового сопротивления принимаем: равным R =0,3 м2 K/Вт. На рис. 2 представлено температурное поле квазистационарного состояния в заданном сечении технологического объекта из титанового сплава ОТ4-1 при расчете в программном комплексе STAR-CD при воздействии энергетического источника (6).

При воздействии на технологический объект из листового титанового сплава ОТ4-1 энергетического источника (6) максимальная температура т вдоль линии перемещения энергетического источника имеет место при у ~ 1,08-10'3 м, т.е. через t« 0,18 с после прохождения центра источника над расчетной точкой. Величина максимальной температуры в центре зоны термического влияния равна max[T(x = 0, у ,z)| _ ]у=1314 K, где h§ - заданная глубина обработки, а при x = ± 440"3 м снижается в 1,47 раза. Ширина изотермы T = 1050 K составляет

ЬТ =1050 K =6,76^10 3 м.

При воздействии энергетического источника (7) практически выполняется условие равенства величин максимальных температур в центре зоны термического влияния и при x = ± 2-10"3 м ; ± 440"3 м:

Рис. 2. Температурное поле квазистационарного состояния в заданном сечении технологического объекта из титанового сплава ОТ4-1 при расчете в программном комплексе STAR-CD

max[T ( Х = 0) z ] у =

= max[T(x = ±2 • 10 3 = max[T(x = ±4 • 10 3

м)

\z=hg

м)

z=hs

] у ] у

1250 ± 12K.

Причем, по сравнению с воздействием энергетического источника (6) max[T(x = 0, у ,z)\ ]

\z=hs у

уменьшается в 1,12 раза, а

max[T(x = ±4• 10~3 м , у ,z)| ]

lz=hs у

увеличивается 1,28 раза. Ширина изотермы T = 1100 K увеличивается в 1,28 раза и принимает значение bT=1100 K =9,33• 10-3 м. Максимальная скорость нагрева точек на глубине hs = 0,4^10-3 м технологического объекта составляет ин « 7• 103 K/с. Скорость охлаждения в диапазоне значений температур от t = 1100 K до T =940 K составляет иол « 103...2403 K/с.

Результаты расчетов показывают, что формирующиеся температурные поля достаточно однородны по всей длине излучения полосового энергетического источника, за исключением переходных зон вблизи краев (± 1.. .2)^ 10-3 м. Обработка поверхности полосовым энергетическим источником имеет ряд технологических преимуществ. В этом случае отдельные зоны лазерного воздействия, в которых достигается одна и та же степень развития структурно-фазовых превращений располагаются параллельно поверхности. Тем самым полосовой тепловой источник может обеспечить наиболее однородное распределение микромеханических характеристик по обработанной поверхности. Учет изменений формы светового пятна и распределения интенсивности приведет к качественным отличиям в распределении температурных полей. Максимальные значения температур могут достигать не в центре зоны термического влияния, а у краев лазерного пятна.

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

769

Известия Самарского научного центра Российской академии наук, т. 10, №3, 2008

на поверхности технологического объекта из титанового сплава ОТ4-1 при воздействии поверхностных энергетических источников (6) (а) и (7) (б):

______- результаты расчета: 1 - у = 1,5*10'3м; 2-у = 4*10'3 м; 3-у = 6,5*10'3 м; 4 - у =9*10'3м.

Экспериментальные данные: *- у = 1,5-10'3 м; д- у = 4*10'3 м; 0 - у = 6,5*10'3 м; п - У =9*10'3м.

экспериментальной оценки распределения температуры на поверхности технологического объекта из титанового сплава ОТ4-1 при воздействии поверхностных энергетических источников (6) и (7) представлены на рис. 3.

Проведенные исследования показали, что результаты расчета температурных полей в программном комплексе STAR-CD в листовом низколегированном титановом сплаве ОТ4-1 при воздействии движущихся распределенных поверхностных источников энергии удовлетворительно (погрешность не превышает 15.. .20%) коррелируют с результатами экспериментальных исследований. Это позволяет сделать вывод о применимости программного комплекса STAR-CD для расчета температурных полей в конструкционных материалах. Применение полосовых энергетических источников для создания требуемого энергетического воздействия на технологические объекты позволяет увеличить ширину зоны обработки без перегрева ее центральных участков.

Заключение

Проведен расчет температурного поля в пакете STAR-CD для технологического объекта из конструкционного материала - листового низколегированного титанового сплава ОТ4-1. Формирующееся температурное поле достаточно однородно по всей длине излучения полосового энергетического источника, за исключением переходных зон вблизи краев. Показано, что полосовой тепловой источник

может обеспечить наиболее однородное распределение микромеханических характеристик по обработанной поверхности.

При проведении экспериментальных исследований для контроля температурных полей применялся инфракрасный термометр “Кельвин”. Проведенные исследования показали, что результаты расчета температурных полей в программном комплексе STAR-CD в листовом низколегированном титановом сплаве ОТ4-1 при воздействии движущихся распределенных поверхностных источников энергии удовлетворительно (погрешность не превышает 15.20%) коррелируют с результатами экспериментальных исследований. Это позволяет сделать вывод о применимости программного комплекса STAR-CD для расчета температурных полей в конструкционных материалах.

Работа выполнена при поддержке аналитической ведомственной целевой программы "Развитие научного потенциала высшей школы (2006-2008 годы)" и гранта РФФИ 08-0899080

СПИСОК ЛИТЕРАТУРЫ

1. Технология лазерной обработки конструкционных и инструментальных материалов в авиадвигателестроении / Р.Р. Латыпов, Н.Г. Терегулов, А.М. Смыслов, А.В. Любанов; под общ. ред. В. Ф. Безъязычного. М.: Машиностроение, 2007.

770

Физика

2. Григорьянц А.Г., Штанов И.К., Мисюров А.И. Технологические процессы лазерной обработки / Подред. А.Г. Григорьянца. М.: МГТУ, 2006.

3. Бек Дж., Блакуэлл Б., Сент-Клэр Ч. Некорректные обратные задачи теплопроводности: Пер. с англ. М.: Мир, 1989.

4. Патанкар С. Численные методы решения задач теплообмена и динамики жидкости: Пер. с англ. М.: Энергоатомиздат, 1984.

5. Murzin S.P. Increasing the efficiency of laser treatment of materials using elements of computer optics // Journal of Advanced Materials, 2003. Vol. 10. №2.

CALCULATION OF TEMPERATURE FIELDS IN CONSTRUCTIONAL MATERIALS IN THE PROGRAM COMPLEX STAR-CD

© 2008 S.P. Murzin, A.V Mezhenin, E.L. Osetrov

Samara State Aerospace University

Calculations of a temperature fields in the program complex STAR-CD for technological object from a constructional material - sheet low-alloy titanic ОТ 4-1 are carried out. The conducted researches have shown that results of calculation of temperature fields in the program complex STAR-CD in sheet low-alloy titanic ОТ4-1 at influence of the moving distributed superficial energy sources well (the error does not exceed 15 ... 20 %) correlate with results of experimental researches. It allows to draw a conclusion on applicability of a program complex STAR-CD for calculation of temperature fields in constructional materials.

771