Научный журнал
Современные наукоемкие технологии
ISSN 1812-7320
"Перечень" ВАК
ИФ РИНЦ = 1,021

РАСЧЕТ ТЕМПЕРАТУРНЫХ ПОЛЕЙ В ПРОЦЕССЕ ЛАЗЕРНОЙ СВАРКИ СТАЛЬНЫХ И АЛЮМИНИЕВЫХ ЛИСТОВ ВНАХЛЕСТ

Гоц А.Н. 1 Гусев Д.С. 2 Гуськов В.Ф. 1 Люхтер А.Б. 1
1 ФГБОУ ВО «Владимирский государственный университет имени Александра Григорьевича и Николая Григорьевича Столетовых»
2 ООО «ИЦ при ВлГУ»
В работе предложена методика выбора источника теплоты и уточнены начальные и граничные условия теплообмена в процессе лазерной сварки стального и алюминиевого листов внахлест. Для моделирования температурных полей в процессе лазерной сварки размеры и форму сварочной ванны выбирают по микрошлифу поперечного сечения сварного шва. Расчет температурных полей свариваемых деталей основывается на решении уравнения нестационарной теплопроводности. Введенная формализация источника теплоты в процессе лазерной сварки стальных и алюминиевых листов внахлест позволила обеспечить корректное задание начальных и граничных условий теплообмена при численном решении методом конечных элементов уравнения нестационарной теплопроводности. Теплосодержание сварочной ванны определяется удельным теплосодержанием металлов свариваемых деталей, температурой плавления наиболее тугоплавкого металла и геометрическими размерами сварочной ванны. Выявлены условия теплообмена, оказывающие превалирующее влияние на отдельные стадии процесса теплопроводности при лазерной сварке листов, что позволяет осуществить целенаправленный выбор начальных и граничных условий теплообмена при нагреве и охлаждении деталей. При экспериментальных исследованиях изменение теплового состояния образца при лазерной сварке регистрировалось хромель-алюмелевыми термопарами. Получена удовлетворительная сходимость результатов расчетных и экспериментальных исследований. Результаты исследований теплового состояния свидетельствуют о значительной температурной неравномерности, возникающей в образце в процессе лазерной сварки, что приведет к появлению в нем остаточных напряжений.
источник теплоты
сварочный шов
температурное поле
сварочная ванна
теплосодержание
граничные условия
1. Sudnik W., Radaj D., Erofeew W. Computerised simulation of laser beam welding, modelling and verification. Journal of Physics D: Applied Physics. 1996. vol. 29. no. 11. P. 2811–2817. DOI: 10.1088/0022-3727/29/11/013.
2. Schwenk C. Modeling of thermomechanical phenomena in fusion welding. Welding Fundamentals and Processes. 2011. vol. 6. P. 830–841. DOI: 10.13140/2.1.2440.9602.
3. Кархин В.А., Хомич П.Н., Иванов С.Ю. Модели источников теплоты для прогнозирования тепловых полей при сварке плавлением // Известия Тульского государственного университета. Технические науки. 2010. № 4–1. C. 241–254.
4. Слезкин Д.В., Цвелев Р.В., Ерофеев В.А., Масленников А.В., Судник В.А. Методика расчетно-экспериментального определения эффективных мощности и радиуса теплового потока дуги при сварке стали плавящимся электродом в смеси защитных газов // Известия Тульского государственного университета. Технические науки. 2012. № 9. С. 172–181.
5. Судник В.А., Ерофеев В.А., Масленников А.В. Методика определения характеристик эквивалентного источника теплоты для выполнения расчётов деформаций при сварке // Известия Тульского государственного университета. Технические науки. 2015. № 6–2. C. 32–43.
6. Люхтер А.Б., Шлегель А.Н., Гусев Д.С., Самарин С.С. Исследование процесса формирования сварного соединения, полученного лазерной сваркой внахлест алюминия АМг2М и стали Ст3 // Материаловедение. 2018. № 1. С. 26–30.
7. Колесников К.С., Румянцев В.В., Леонтьев А.И., Полежаев Ю.В. и др.; под общ. ред. Колесникова К.С., Леонтьева А.И. Машиностроение. Энциклопедия. Теоретическая механика. Термодинамика. Теплообмен. М.: Машиностроение, 1999. Т. 1–2. 600 с.

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

Цель исследования: уточнение выбора источника теплоты и начальных и граничных условий теплообмена в процессе лазерной сварки стального и алюминиевого листов внахлест путем сравнения расчетных и экспериментальных температурных полей при расчете температурных полей методом конечных элементов (МКЭ).

Материалы и методы исследования

В качестве объекта исследования принята лазерная сварка опытного образца (рис. 1), осуществляемая непрерывным излучением без сквозного проплавления нижней алюминиевой заготовки и использования присадочного материала [6].

goc1a.tif

Рис. 1. Сварной образец: 1 – сталь 20, толщина 1,5 мм; 2 – алюминиевый сплав АМг2М, толщина 3,0 мм; 3 – сварной шов

goc2a.tif goc2b.tif

а) б)

Рис. 2. Модель сварочной ванны: а – микрошлиф поперечного сечения сварного шва: 1 – сталь 20; 2 – алюминиевый сплав АМг2М; 3 – сварной шов; б – модель сварочной ванны: bsv – ширина сварочной ванны, dp – диаметр пятна лазерного излучения

Микрошлиф поперечного сечения сварного шва и модель сварочной ванны, соответствующая поперечному сечению шва, показаны на рис. 2.

Расчет температурных полей свариваемых образцов (рис. 1) основывается на решении уравнения нестационарной теплопроводности методом конечных элементов:

goch01.wmf, (1)

где goch02.wmf – оператор Лапласа; T – температура детали; goch03.wmf – температуропроводность соответствующего материала; ρ, с, λ – соответственно плотность, теплоемкость, теплопроводность свариваемых образцов; τ – текущее время; Q – производительность внутренних источников теплоты (в расчетах Q = 0), для заданных условий однозначности, наибольшая сложность определения которых заключается в формировании начальных и граничных условий теплообмена.

При их назначении в целях упрощения целесообразно исходить из следующих допущений:

1. Нагревание свариваемых деталей образца осуществляется за счет теплосодержания расплава металлов в объеме сварочной ванны.

2. Охлаждение свариваемых деталей и расплава металлов происходит в результате теплопроводности и свободной конвекции в окружающую среду.

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

4. Скорость движения лазерного луча (сварочной ванны) vsv, диаметр пятна луча лазера dp и продолжительность нагрева от излучения τsv связаны соотношением

goch04.wmf. (2)

5. Термическое сопротивление контакта между соединяемыми деталями образца отсутствует.

6. Расплав металла в объеме сварочной ванны имеет теплофизические свойства наиболее тугоплавкого металла.

Первое допущение позволяет определить теплосодержание единицы длины сварного шва по формуле

goch05.wmf, (3)

где cs, ca – удельное теплосодержание при температуре плавления стали соответственно для стали и алюминия; ms, ma – массы стали и алюминия в сварной ванне; t0 = 293 K – начальная температура; rs, ra – удельная теплота плавления стали и алюминия; l – длина сварочного шва.

Для определения масс ms и ma расплавленных металлов поверхность сварного шва представляется в виде квадратов со стороной bsv, описывающих пятно лазерного луча, а глубина сварочной ванны определяется с использованием микрошлифа поперечного сечения сварного шва для используемого технологического режима сварки (рис. 2).

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

На рис. 3 схематично представлено перемещение объема сварочной ванны и изменение температуры расплава металла в ней по направлению сварки.

goc3a.tif goc3b.tif

а) б) в)

Рис. 3. Схема задания начальных и граничных условий теплообмена в сварном шве: а – схема для шага расчета τsv1; б – схема для шага расчета τsv2: буквами обозначены поверхности, контактирующие со стальным (А, Б, В) и алюминиевым (И, К, Л) листами, окружающей средой (Г, Д, Е) и соседней сварочной ванной (Ж); в – расчетный термический цикл для сварочных ванн в момент времени τi и τi + τsv

В момент времени τi температура в объеме ванны 1 мгновенно повышается до температуры плавления стали. В течение времени τsv осуществляется перенос теплоты от сварочной ванны 1 через ее граничные поверхности А, И, Ж в металл заготовки. По завершению времени τi+τsv в качестве источника теплоты выступает последовательно следующий 2 и другие объемы сварочной ванны по направлению сварки с аналогичным интервалом времени τsv и граничными поверхностями. Охлаждение сварного шва и образца при перемещении источника теплоты и окончании сварки осуществляется путем теплопроводности в металл образца и свободной конвекции в окружающую среду.

В соответствии с описанной моделью нагрева и охлаждения образца при сварке начальные условия теплообмена τ = 0: Т = f (х, y, z), где Т – температура плавления стали; x, y, z – координаты изотермической поверхности (поверхности сварочной ванны); Т = Т0, где Т0 = 293 K – начальная температура образца.

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

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

Теплофизические параметры материалов сварных образцов задаются в зависимости от температуры (таблица) по справочным данным [7].

Теплофизические характеристики стали 20 и АМг2М в зависимости от температуры

Характеристики

Материал

Температура, °С

0

100

300

600

1000

1530

Теплопроводность, Вт/(м·К)

Сталь 20

55

52

48

34

28

25

АМг2М

155

159

168

173

174

175

Коэффициент теплового расширения х10-6, 1/К

Сталь 20

11,1

12,3

13,8

15,1

13,3

13,3

АМг2М

23

24,9

28,8

30

30,5

31

Плотность, кг/м3

Сталь 20

7859

7834

7770

7659

7600

7600

АМг2М

2680

Удельное теплосодержание, Дж/(кг·К)

Сталь 20

325

358

446

579

628

760

АМг2М

904

963

1037

1228

1177

1177

Удельная теплота плавления, кДж/кг

Сталь 20

271

АМг2М

390

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

Изменение теплового состояния образца при сварке целесообразно регистрировать хромель-алюмелевыми термопарами. Схема установки термопар в образце и функциональная схема измерения температуры приведены на рис. 4.

goc4.tif

Рис. 4. Функциональная схема измерения температуры

Термопары устанавливались в срединной плоскости стального листа, расположенной перпендикулярно к сварному шву, в пяти сверлениях диметром 1 мм и глубиной 0,3 мм на расстоянии 0,3; 1,0; 1,7; 3,5; 6,0 мм от кромки сварного шва. Контакт термопары с образцом обеспечивался их приваркой к металлу с помощью разряда конденсатора. Экспериментальное измерение температур проводилось при сварке пяти образцов.

Сварка осуществлялась на лазерном роботизированном комплексе ЛРК-С с максимальной выходной мощностью лазера 3 кВт в лаборатории лазерной сварки научно-образовательного центра ВлГУ.

Непрерывная регистрация сигналов от термопар осуществлялась с помощью самописца цифрового многоканального S-Recorder L. Началом регистрации являлся сигнал включения лазерного источника, который поступал с контроллера робота и дублировался на вольтметр. Последующий переход от величины сигнала термопары к значениям температур осуществлялся через номинальные статические характеристики преобразования термопар по ГОСТ Р 8.585-2001.

Результаты исследования и их обсуждение

Результаты расчетного и экспериментальных исследований теплового состояния образца (рис. 5) при сварке приведены на рис. 6. Разработанная для моделирования теплового состояния трехмерная конечно-элементная модель образца со средними размерами конечного элемента в зоне термического влияния 0,05–0,2 мм, а в сварочном шве и околошовной области – 0,5 мм состояла из около 670 тыс. элементов и 3,88 млн узлов. Расчетные величины теплосодержания единицы длины сварного шва Qsv с принятой формой поперечного сечения (рис. 2) и времени нагрева τsv составляют для образца (рис. 1) и режима лазерной сварки со скоростью движения пятна излучения лазера диаметром dp = 2,56 мм, равной vsv = 18 мм/с соответственно Qsv = 535 Дж/см и τsv = 0,142 с.

goc5.tif

Рис. 5. Зависимость температуры стальной заготовки в виде функции времени при лазерной сварке образца: непрерывная линия – данные эксперимента, пунктирная – данные компьютерного моделирования при Qsv = 535 Дж/см: 1 – 0,3 мм от шва, 2 – 1,0 мм от шва, 3 – 1,7 мм от шва, 4 – 3,5 мм от шва, 5 – 6,0 мм от шва

goc6.tif

Рис. 6. Распределение температур в сварном образце в ходе и после сварки: а–г – через 0,142; 0,426; 1,278 и 1,664 с после начала сварки соответственно; д – тепловой след хвостовой части сварочной ванны; е – по длине сварного шва (точка В) в плоскости соприкосновения стального и алюминиевого листов при Qsv = 535 Дж/см по завершению сварки через: 1 – 0 с; 2 – 0,2 с; 3 – 0,4 с; 4 – 0,6 с; 5 – 1,4 с

Анализ результатов показывает, что относительная погрешность между значениями температур образца, полученных в расчете и эксперименте, при нагреве и охлаждении в процессе сварки не превышает соответственно 10 % и 7 %.

Согласно экспериментальным данным воздействие лазерного луча приводит к быстрому нагреву за 0,2 с околошовной зоны образца до максимальной температуры, затем в течение короткого промежутка времени τexp = 0,15 с температура практически не изменяется и при дальнейшем перемещении лазерного луча температура уменьшается с 1109 °С до 442 °С за 1 с. Следует отметить, что временной интервал неизменной максимальной температуры в непосредственной близости (0,3 мм) от сварного шва, полученный в эксперименте τexp = 0,15 с, практически не отличается от расчетного τsv = 0,142 с. После завершения сварки через 2 с наступал регулярный тепловой режим.

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

На рис. 6 представлены результаты расчетов теплового состояния образца в процессе лазерной сварки и на стадии охлаждения.

Распределение температур в сварном образце в процессе сварки подтверждает принятое допущение о форме сварочной ванны как источника теплоты (рис. 6, д).

Распределение температур по длине сварного шва в плоскости соприкосновения стального и алюминиевого листов (точка В поперечного сечения шва, рис. 6, е) позволяет оценить неравномерность температур, возникающую в процессе сварки. В частности, температура в точке В на участке сварки длиной 30 мм уменьшается с температуры плавления стали 1530 °С до температуры 440 °С, а через 1,4 с после завершения сварки температура выравнивается до уровня 280–300 °С.

Заключение

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

Работа выполнена при финансовой поддержке Министерства образования и науки Российской Федерации. Соглашение о предоставлении субсидии № 14.577.21.0158 от 28 ноября 2014 г. Уникальный идентификатор ПНИЭР RFMEFI57714X0158.


Библиографическая ссылка

Гоц А.Н., Гусев Д.С., Гуськов В.Ф., Люхтер А.Б. РАСЧЕТ ТЕМПЕРАТУРНЫХ ПОЛЕЙ В ПРОЦЕССЕ ЛАЗЕРНОЙ СВАРКИ СТАЛЬНЫХ И АЛЮМИНИЕВЫХ ЛИСТОВ ВНАХЛЕСТ // Современные наукоемкие технологии. – 2019. – № 2. – С. 32-37;
URL: http://top-technologies.ru/ru/article/view?id=37404 (дата обращения: 30.11.2020).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.074