Численное моделирование процесса электронно-лучевой сварки с продольной осцилляцией луча на основе экспериментально определенной формы канала проплавления

Полный текст

Введение. В настоящее время широко известны успехи, достигнутые в области численного моделирования электронно-лучевой и лазерной сварки [1-3]. Однако все результаты связаны с моделированием сварки статическим пучком. Законченные динамические модели сварки с осцилляциями пучка отсутствуют до сих пор. Комплексный характер, высокая скорость протекания процессов, высокие значения градиентов температур и многофакторность процессов делают непосредственное численное моделирование в условии периодических воздействий крайне затруднительным даже с использованием современных вычислительных ресурсов [1]. Основная трудность при моделировании - определение формы канала проплавления. Для решения этой проблемы ранее была предложена методика экспериментального определения формы канала проплавления при электронно-лучевой сварке с осцилляцией пучка по параметрам вторичного сигнала [4]. Кроме формы канала проплавления указанная методика дает дополнительную информацию о процессах в канале проплавления. В частности, данная методика позволяет определить распределение энергии пучка по стенкам канала проплавления [5]. Такой подход избавляет от необходимости учета всех сложных факторов, оказывающих воздействие на формирование канала проплавления. Целью данного исследования является использование данных параметров в качестве входных для проведения математических расчётов, определение влияния различных осцилляций пучка и их параметров на форму канала проплавления, протекание процессов тепломассопереноса и параметры швов для разработки методических рекомендаций. Расчёты тепловых и гидродинамических процессов производились на кластере с использованием пакета моделирования COMSOL Multiphysics. В результате расчётов была получена геометрия сварного шва, поля температур и скоростей течения расплава. Математическая модель. На рис. 1 показана схема модели. Параметры процесса и физические параметры стали 15Х5М представлены в табл. 1 и 2 соответственно. Следующие упрощения были сделаны при моделировании процесса сварки: - ламинарное течение расплава; - нормальное распределение энергии в пучке; - линейный, вдоль оси X, источник тепла; - моделирование фазового перехода с помощью скачка вязкости и введения скрытой теплоты фазового перехода; - форма и размер канала проплавления были определены из экспериментальных данных; канал аппроксимировали усечённым конусом со сферической вершиной. Реконструкция формы канала. Реконструкция формы канала проплавления и получение данных о распределении энергии в нем осуществлялись посредством регистрации вторичных излучений из зоны сварки методом синхронного накопления [6; 7]. На рис. 2 представлены результаты реконструкции формы канала при электронно-лучевой сварке с осцилляцией вдоль стыка на трех режимах фокусировки [3; 4]. Для наглядности масштаб вдоль вертикальной и горизонтальной осей выбран неодинаковым. Видно, что при изменении режима фокусировки изменяется характер взаимодействия пучка со стенками канала проплавления. При сварке на недофокусированном режиме электронный пучок взаимодействует в основном с передней стенкой канала проплавления, на перефокусированном - с передней и задней стенками практически одинаково. Полученная форма канала проплавления и распределение энергии пучка в нем использовались при моделировании в качестве входных параметров. Рис. 1. Геометрия модели, начальные и граничные условия Таблица 1 Параметры процесса Тип осцилляции Параметры пучка Параметры осцилляции Скорость сварки Vсв, мм/с Мощность Q, кВт Min радиус, мм Амплитуда, мм Частота, Гц Продольная 3 0,25 0,75 645 5 Таблица 2 Физические характеристики стали 15Х5М Характеристика Обозначение Величина Начальная температура T0 [K] 293,15 Ликвидус TL [K] 1730 Солидус TS [K] 1700 Скрытая теплота плавления Lf [Дж·кг-1] 700 Динамическая вязкость μ [Па·с] 0,007 Плотность ρ [кг·м-3] 7750 Удельная теплоёмкость c [Дж·кг-1·K-1] 483 Теплопроводность λ [Вт·м-1·K-1] 37 а б в Рис. 2. Результаты реконструкции канала проплавления (продольная проекция) по параметрам вторичного сигнала при электронно-лучевой сварке с осцилляцией вдоль стыка для различных режимов фокусировки: а - перефокусированный; б - острый фокус; в - недофокусированный Основные уравнения. Модель основана на совместном решении тепловой и гидродинамической задачи обтекания жидким металлом канала проплавления в трехмерной постановке. Рассмотрим основные уравнения для расчетов [1; 8]. Течение жидкого металла описывается уравнениями непрерывности и движения для несжимаемой жидкости в подвижной системе координат: (1) (2) где r - плотность; - вектор скорости течения расплава; m - динамическая вязкость; P - давление; F - слагаемое, представленное в виде (3) где g - ускорение свободного падения; a - коэффициент термического расширения; T - абсолютная температура; Tref - исходная температура, принятая как температура солидуса (TS); Vсв - скорость сварки. Первый член уравнения представляет термогравитационную конвекцию, второй член уравнения представляет диссипацию энергии в двухфазной зоне в соответствии с уравнением Кармана-Козени [9; 10], где B - небольшая вычислительная постоянная - введена для избегания деления на ноль; C - константа учёта морфологии двухфазной зоны (в данных исследованиях использовалась величина, равная 1,6·104 [10]); fL - жидкая составляющая [11]: (4) где TL и TS - температура ликвидуса и солидуса соответственно. Последний член уравнения (3) описывает движение заготовки относительно электронного пучка. Распределение тепловой энергии источника в заготовке может быть описано с помощью дифференциального уравнения переноса энергии: (5) где T - абсолютная температура; - коэффициент температуропроводности; u - скорость течения расплава; c - теплоёмкость; r - плотность; Q - подводимая тепловая мощность. Скрытая теплота плавления и кристаллизации учитывалась введением эффективной теплоемкости (6) где С0 - теплоемкость в зависимости от температуры; Hf - скрытая теплота плавления; Tmelt - температура плавления, которая принята средней в интервале от температуры солидуса до ликвидуса. На верхней поверхности сварочной ванны и на стенках канала проплавления действует эффект Марангони. Для этих поверхностей в проекциях на оси X, Y и Z можно записать [12-15]: (7) где u, v, w - компоненты скорости вдоль осей X, Y и Z соответственно; m - динамическая вязкость; fL - жидкая составляющая; - температурный коэффициент поверхностного натяжения; T - абсолютная температура. Результаты. На рис. 3 представлены экспериментальные и расчётные сечения сварных швов при сварке с продольной осцилляцией для различных режимов фокусировки. Экспериментальные и расчётные данные для острого фокуса и перефокусированного режима хорошо согласуются. Различия между расчётами и экспериментом для перефокусированного режима вызваны ошибками при определении глубины канала проплавления. Проведен анализ значимости влияния различных факторов на формирование геометрии проплавления и на вероятность формирования дефектов при электронно-лучевой сварке с осцилляцией электронного пучка. Одной из причин специфической формы проплавления с уширением в верхней части являются конвективные явления в жидкой фазе (характерные значения чисел Марангони и Рэлея составили 5∙103 и 2∙103). Интенсивность этих явлений зависит от распределения вводимой энергии в канале проплавления и, как следствие, от режима фокусировки. Установлено, что интенсивность конвективных течений, обусловленных эффектом Марангони, повышается при переходе от недофокусированного режима к перефокусированному. Интенсивность термогравитационной конвекции при этом несколько снижается. Выводы. Математическая модель позволяет проанализировать влияние параметров осцилляций на геометрию швов через изучение влияния геометрии канала на тепломассоперенос в жидкой ванне. Методика, описанная в статьях [4; 5], может быть использована для определения геометрии канала проплавления.

Об авторах

Г. Л. Пермяков

Пермский национальный исследовательский политехнический университет

Email: gleb.permyakov@yandex.ru
Российская Федерация, 614990, г. Пермь, Комсомольский проспект, д. 29

Д. Н. Трушников

Пермский национальный исследовательский политехнический университет

Российская Федерация, 614990, г. Пермь, Комсомольский проспект, д. 29

В. Я. Беленький

Пермский национальный исследовательский политехнический университет

Российская Федерация, 614990, г. Пермь, Комсомольский проспект, д. 29

Т. В. Ольшанская

Пермский национальный исследовательский политехнический университет

Российская Федерация, 614990, г. Пермь, Комсомольский проспект, д. 29

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

Дополнительные файлы