Numerical simulation of electron-beam welding with a longitudinal oscillation of the beam based on the experimental determination of the keyhole form


Citar

Texto integral

Resumo

The research examines the process of electron-beam welding in the keyhole mode with the use of beam oscillations with different focusing regimes on the example of stainless steel - grade 502. We study the impact of longitudinal beam oscillation and their parameters on the shape of the keyhole, the character of the processes of heat and mass transfer and weld’s parameters, to develop methodological recommendations for electron-beam welding with oscillation. A numerical three-dimensional mathematical model of electron beam welding is presented. The model was developed on the basis of the combined solution of heat conduction equation and Navier-Stokes equation in the moving coordinate system with taking into account phase transitions at the interface of solid and liquid phase. The boundary conditions at the free surface of the weld pool and on the keyhole walls are presented by thermocapillary convection (Marangoni effect). As input parameters the form of keyhole and distribution of the beam energy on the keyhole walls based on experimental data on the parameters of the secondary signal in the plasma over the welding zone by using the method of a synchronous accumulation are used. The keyhole was approximated by an oblique elliptical cone with a spherical apex. Such an approach eliminates the need to take into account all the complex factors that affect the formation of the keyhole. Calculations of thermal and hydrodynamic processes have been carried on a cluster using a simulation package COMSOL Multiphysics. Comparison of the calculated and experimental cross sections of the welds showed good agreement. With the help of the proposed method during the research characteristic features of the beam power distribution for different focus modes were identified. The analysis of the significance of various factors on the formation of the keyhole geometry has been conducted.

Texto integral

Введение. В настоящее время широко известны успехи, достигнутые в области численного моделирования электронно-лучевой и лазерной сварки [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], может быть использована для определения геометрии канала проплавления.
×

Sobre autores

G. Permyakov

Perm National Research Polytechnic University

Email: gleb.permyakov@yandex.ru
29, Komsomolskiy Av., Perm, 614990, Russian Federation

D. Trushnikov

Perm National Research Polytechnic University

29, Komsomolskiy Av., Perm, 614990, Russian Federation

V. Belenkiy

Perm National Research Polytechnic University

29, Komsomolskiy Av., Perm, 614990, Russian Federation

T. Olshanskaya

Perm National Research Polytechnic University

29, Komsomolskiy Av., Perm, 614990, Russian Federation

Bibliografia

  1. Rai R., Palmer T. A., Elmer J. W., Debroy T. Heat transfer and fluid flow during electron beam welding of 304L stainless steel alloy, Welding Journal, March 2009, Vol. 88, No. 3, P. 54-61.
  2. Sudniky W., Radajz D. and Erofeew W. Computerized 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.
  3. Cho W.-I., Na S.-J., Thomy C., Vollertsen F. Numerical simulation of molten pool dynamics in high power disk laser welding, Journal of Materials Processing Technology, 2012, Vol. 212, Iss. 1, p. 262-275. doi: 10.1016/j.jmatprotec.2011.09.011.
  4. Trushnikov D., Belenkiy V., Schavlev V., Piskunov A., Abdulin A., Mladenov G. Plasma charge current for control and monitoring at electron beam welding with the beam oscillation, Sensors, 2012, Vol. 12(12), P. 17433-17445. doi: 10.3390/s121217433.
  5. Trushnikov D. N., Belen’kii V. Ya. Investigation of the formation of the secondary current signal in plasma in electron beam welding with oscillations of the electron beam, Welding International, 2013, Vol. 27, Iss. 11,
  6. P. 877-880. doi: 10.1080/09507116.2013.796645.
  7. Max J. Méthodes et techniques de traitement du signal et applications aux mesures physiques, 3rd ed.; Jean-Louis Lacoume, Paris, France, 1981, Vol. I,
  8. P. 10-50.
  9. Orfanidis S. J. Optimum Signal Processing. An Introduction, 2nd ed.; Prentice-Hall : Englewood Cliffs, NY, USA, 1996, P. 5-30.
  10. Bachmann M., Avilov V., Gumenyuk A., Rethmeier M. Experimental and numerical investigation of an electromagnetic weld pool support system for high power laser beam welding of austenitic stainless steel, Journal of Materials Processing Technology, 2014,
  11. Vol. 214, Iss. 3, P. 578- 591. doi: 10.1016/j.jmatprotec. 2013.11.013.
  12. Voller V. R., Prakash C. A fixed grid numerical modelling methodology for convection-diffusion mushy region phase-change problems, International Journal of Heat and Mass Transfer, 1987, Vol. 30, Iss. 8, P. 1709-1720. doi: 10.1016/0017-9310(87)90317-6.
  13. Brent A. D., Voller V. R., Reid K. J. Enthalpy-porosity technique for modeling convection-diffusion phase change: Application to the melting of a pure metal, Numerical Heat Transfer, 1988, Vol. 13, Iss. 3, P. 297-318. doi: 10.1080/10407788808913615.
  14. Zhang W., Kim C. H., DebRoy T. Heat and fluid flow in complex joints during gas-metal arc welding, Part I: Numerical model of fillet welding, Journal of Applied Physics, 2004, Vol. 95, Iss. 9, P. 5210-5219. doi: 10.1063/1.1699485.
  15. Kou S., Sun D. K. Fluid flow and weld penetration in stationary arc welds, Metalurgical Transactions A, 1985, Vol. 16, Iss. 2, P. 203-213. Doi: 10.1007/ BF02815302.
  16. Kim C. H., Zhang W., DebRoy T. Modeling of temperature field and solidified surface profile during gas metal arc fillet welding, Journal of Physics, 2003, Vol. 94, Iss. 4, P. 2667-2679. doi: 10.1063/1.1592012.
  17. De A., DebRoy T. Probing unknown welding parameters from convective heat transfer calculation and multivariable optimization, Journal of Physics D: Applied Physics, 2004, Vol. 37, No. 1, P. 140-150. Doi: 10.1088/ 0022-3727/37/1/023.
  18. De A., DebRoy T. A smart model to estimate effective thermal conductivity and viscosity in weld pool, Journal of Applied Physics, 2004, Vol. 95, Iss. 9,
  19. P. 5230-5240. doi: 10.1063/1.1695593.

Arquivos suplementares

Arquivos suplementares
Ação
1. JATS XML
Baixar

Declaração de direitos autorais © Permyakov G.L., Trushnikov D.N., Belenkiy V.Y., Olshanskaya T.V., 2015

Creative Commons License
Este artigo é disponível sob a Licença Creative Commons Atribuição 4.0 Internacional.

Este site utiliza cookies

Ao continuar usando nosso site, você concorda com o procedimento de cookies que mantêm o site funcionando normalmente.

Informação sobre cookies