CC BY
38
УДК 621.791.72
Е.С. Саломатова, Д.Н. Трушников, А.И. Цаплин, В.Я. Беленький
Пермский национальный исследовательский политехнический университет,
Пермь, Россия
ПРИМЕНЕНИЕ ПАРАЛЛЕЛЬНЫХ ВЫЧИСЛЕНИЙ ДЛЯ ИССЛЕДОВАНИЯ ИСПАРЕНИЯ ПРИ ЭЛЕКТРОННО-ЛУЧЕВОЙ СВАРКЕ
Формирование парогазового канала при электронно-лучевой сварке (ЭЛС) сопровождается достаточно сложным процессами, такими как локальные перегревы металла с последующими взры-вообразными вскипаниями (абляцией), стохастическое перемещение зоны взаимодействия пучка электронов с металлом по глубине парогазового канала, периодическое смыкание канала в верхней части, образование волн на поверхности расплава и др. Для описания процесса ЭЛС до последнего времени в основном применяются упрощенные модели, в которых канал проплавления заменяется линейным источником тепла в комбинации с точечным поверхностным источником, располагающимся на поверхности свариваемого образца. Такие модели дают удовлетворительное описание картины на некотором удалении, но абсолютно не пригодны для исследования самого парогазового канала. В предлагаемой работе описывается оригинальная методика косвенного измерения давления паров и температуры в парогазовом канале при электронно-лучевой сварке. Эта методика основана на экспериментальном определении концентрации химических элементов в парах над зоной сварки и сопоставлении этих результатов с результатами моделирования. 2D-модель основана на совместном решении уравнений диффузии с учетом тепло- и массопереноса в расплаве. Расчеты производились на вычислительном кластере в пакете прикладного программного обеспечения Comsol Multiphysics с применением технологии параллельных вычислений.
Ключевые слова: электронно-лучевая сварка, температура на стенках парогазового канала, термодинамические расчеты, давление паров, диффузия элементов в расплаве, применение параллельных вычислений, тепломассоперенос в расплаве, экспериментальное определение концентрации химических элементов, программное обеспечение Comsol Multiphysics, определение потерь легирующих компонентов.
E.S. Salomatova, D.N. Trushnikov, A.I. Tsaplin, V.Ya. Belenkiy
Perm National Research Polytechnic University, Perm, Russian Federation
THE HIGH PERFORMANCE COMPUTING FOR RESEARCH
OF EVAPORATION DURING ELECTRON BEAM WELDING
The formation of vapor-gas channel in electron beam welding is accompanied by rather complex processes, such as local overheating of the metal followed by explosively Effervescence (ablation), stochastic movement zone of interaction of the electron beam with the metal vapor-gas channel depth, periodic closing of the channel in the upper part, the formation of waves on the surface melt and others.
to describe the process of EBW until recently mainly used simplified models, in which the channel penetration replaced linear heat source in combination with a surface point source located on the surface of the welded sample. Such models provide a satisfactory description of the picture at a certain distance, but it is absolutely not suitable for the study of the vapor-gas channel. In the present paper describes an original technique for indirect measurement of the vapor pressure and temperature in the vapor-gas channel in electron beam welding. This technique is based on an experimental determination of the concentration of chemical elements in the vapor above the welding zone and to compare these results with the simulation results. 2D model is based on the simultaneous solution of diffusion equations with the heat and mass transfer in the melt. The calculations were performed on the computer cluster in the package application software Comsol Multiphysics using the technology of parallel computing.
Keywords: electron-beam welding, the temperature on the walls of the vapor-gas channel, thermodynamic calculations, vapor pressure, diffusion of elements in the melt, the use of parallel computing, the heat-mass transfer in the melt, the experimental determination of the concentration of chemical elements, the software Comsol Multiphysics, the definition of loss of alloying components.
Введение
Электронно-лучевая сварка (ЭЛС) занимает особое место благодаря высокой концентрации мощности в сварочном электронном пучке и его способности глубокого проникновения в металл, что обусловливает ее широкое применение при изготовлении деталей ответственных изделий из различных сплавов. Важным преимуществом является также ведение процесса ЭЛС в вакууме (Р = 10-2 Па), что обеспечивает эффективную защиту зоны сварки от воздействия газовой среды.
При воздействии мощного концентрированного электронного пучка на металл давление паров, покидающих зону плавления, деформирует поверхность сварочной ванны, и в расплаве формируется глубокий и узкий парогазовый канал, жидкие стенки которого удерживаются давлением паров металла [1-4]. В парогазовом канале проходит интенсивное испарение материала, что может приводить к снижению концентрации некоторых легирующих элементов в расплавленном металле. Такой процесс имеет место для металлических сплавов, содержащих легкоплавкие легирующие элементы. Поэтому процессы испарения при ЭЛС представляют интерес при разработке теоретических моделей процесса формирования сварного шва [5, 6].
В работе [7] рассмотрены процессы испарения, в которых диффузия примесей в расплаве, протекающая в тонком слое у поверхности парогазового канала, является основным лимитирующим процессом для определения потери легирующих компонентов при ЭЛС. При этом не учитывается зависимость коэффициента диффузии от температуры расплава.
Целью данной работы является построение нелинейной теоретической модели, которая описывает неравновесные процессы в жидкой фазе металла при ЭЛС с учетом процессов испарения химических эле-
ментов со стенок канала, конденсации элементов сплава на стенки парогазового канала и диффузии элементов в расплаве.
Постановка задачи. Определяющее уравнение и краевые условия
Для построения модели были приняты следующие предположения: форма парогазового канала аппроксимировалась цилиндром; поверхность парогазового канала считалась изотермической, равной температуре окружающей среды; перенос тепла и массы вещества вдоль оси канала считался несущественным; коэффициент диффузии в расплаве принимался зависимым от температуры; истощение расплава легирующими элементами происходило в тонком слое у поверхности парогазового канала; из-за высокого давления в парогазовом канале скорость паров металла на выходе достигала скорости звука.
Таким образом, задача сводилась к двухмерной, а расчетная область содержала только жидкую фазу. Для решения задачи был использован численный метод, позволяющий учесть процессы испарения, конденсации и диффузии элементов в расплаве. Расчетная схема представлена на рис. 1.
Рис. 1. Расчетная схема: 1 - входная граница; 2 - боковая граница; 3 - выходная граница; 4 - ось симметрии; 5 - парогазовый канал
Система уравнений, описывающая процессы испарения при ЭЛС, включает уравнение переноса импульса (уравнение Навье - Стокса)
! - \ - УР -(и ■У)и = — + чАи,
У ' Р уравнение переноса энергии
реи-Ут=у-т?г, (1)
Применение параллельных вычислений для исследования испарения при ЭЛС уравнение переноса примеси
О .уд,. = А -V2 5,
где О - поле скоростей расплава; Р - давление в расплаве; р - плотность расплава; V - коэффициент кинематической вязкости; ср - теплоемкость при постоянном давлении; X - коэффициент теплопроводности; - концентрация ,-го компонента в расплаве; Т - температура на
д2 д2 дТ дТ
--1--; У =--1--
дх2 ду2 дх ду
центрация ,-го компонента в расплаве; Б, - коэффициент диффузии
,-го элемента расплава.
В соответствии с работой [8] коэффициент диффузии принимался
зависящим от температуры:
' а 1
стенках парогазового канала, К; V2 = +--~ ; У = ^—+ ; - кон-
Б г = Б0,. ехр
V 8,13 •
где Б0, - амплитудное значение диффузии ,-го элемента расплава; б -
энергия активации ,-го компонента сплава; численные значения для конкретной марки стали представлены в табл. 1.
Таблица 1
Составляющие характеристики коэффициента диффузии для легирующих компонентов стали 12Х18Н10Т [8]
Составляющие характеристики коэффициента диффузии Основные легирующие элементы стали 12Х18Н10Т
Сг Мп
Б0,, м2/с 1,318 ■ 10-7 2,515 ■ 10-7
б,, кДж/кг 3,92 3,99
Плотность потока массы ,-го легирующего элемента на стенке парогазового канала определяется суммой плотностей потоков, обусловленных испарением (¿¡е\) и конденсацией (/с,):
Ji = Jevi + Зс i.
Плотность потока, обусловленного испарением Jevi, прямо пропорциональна концентрации примеси:
ЗвУ; = 3 0,. ^
V 0,
где 3 0, - плотность потока испарения над чистым ,-м элементом; $0,,
- соответственно начальная и текущая концентрации ,-го элемента в сплаве, последнее определяется из решения уравнения (1).
В свою очередь, плотность диффузионного потока испарения ,-го элемента определяется в соответствии с работой [9] по формуле
Р (Т)
30, Л-
где д, - масса в молях ,-го элемента сплава; ЗУ, - мольный поток ,-го
элемента сплава; А - размерный коэффициент; р (Т) - давление ,-го
элемента расплава, определяется по методике, описанной в работе [10].
Плотность потока, конденсирующегося из парогазовой фазы на стенки канала проплавления, Л, пропорциональна средней концентрации элементов в парах:
= к^,
, V
$ 0,
где к - коэффициент конденсации, 0 < к < 1 (для насыщенного пара к = 1); - среднее значение концентрации на поверхности парогазового канала ,-го элемента.
Коэффициент к определяется из условия равенства суммарного потока массы на стенках парогазового канала ХЗ потоку массы, уносимой через его выходное сечение:
Е ( - Зе, )'$канала =^^,
где $канала - это площадь поперечного сечения парогазового канала; дт
— - поток массы паров металла, уносимых через выходное сечение дt
парогазового канала. Этот поток массы определяется из условия достижения на выходе из парогазового канала местной скорости звука [11]:
1 + у
йт г-— ( 2 ^
— Ч ТРкРк —7 ш ^ у +1
где у - показатель адиабаты; Рк и рк - давление и плотность паров в критическом сечении на выходе из парогазового канала соответственно.
Система дифференциальных уравнений имела следующие граничные условия: на входной границе задачи - постоянное значение концентрации и температуры £ (х = 0, Ь; у) = £0, и Т (х = 0, Ь; у) = Тп. На выходной границе - условие продолжимости (постоянство плотности теплового потока):
(х=Ь, у ) =
Боковая граница 2 является адиабатной:
дТ (х=а н)=0.
дх
По нормали к оси симметрии (у = 0) задавались нулевые потоки массы и энергии:
д-Т=0 и &=0.
Эх Эх
На поверхности парогазового канала, которой соответствует граница 5, задавалась постоянная температура
Т = Тп,
где Тп - температура на поверхности парогазового канала, и плотность потока массы
А ^=А.
дх
Результаты численного исследования
Численное исследование процесса испарения химических элементов со стенок канала, конденсации элементов сплава на стенки парогазового канала и диффузии элементов в расплаве проводилось для
2 2(7-1)
. £
канала 5
расчетной области при Н = 1,5 • 10-3 м; L = 3 • 10-3 м; г = 2,5 • 10-4 м; l = 1,6 • 10-3 м. Исходные данные представлены в табл. 2.
Основой экспериментальной методики являлось определение температуры в канале проплавления при ЭЛС по данным о химическом составе паров, образующихся над зоной сварки. Для экспериментальной оценки состава паров на электронно-лучевой установке с энергетическим агрегатом ЭЛА-6ВЧ фирмы 8БЬМ1 (Украина) производились проходы осциллирующим электронным пучком по образцу из аустенитной стали 12Х18Н10Т. Вблизи зоны воздействия электронного пучка на образец (на расстоянии 15 мм) была установлена пластина из стекла. Далее проводился рентгенофлуоресцентный анализ полученного покрытия, данные представлены в табл. 3.
Таблица 2
Данные, используемые для расчета при ЭЛС стали 12Х18Н10Т
№ п/п Характеристики Значение
1 Плотность, кг/м3 7810
2 Теплопроводность, Вт/(мК) 30,6
3 Теплоемкость при постоянном давлении, кДж/(кгК) 0,830
4 Динамическая вязкость, Па-с 110-3
5 Рассчитанный коэффициент конденсации 0,97
6 Начальная температура на стенках канала проплавления ^ К 1800
Таблица 3
Первоначальные и экспериментальные значения концентраций элементов в стали 12Х18Н10Т
Элементы Концентрация (доли)
в основном металле полученного покрытия
Железо 0,676 0,336
Марганец 0,02 0,350
Хром 0,19 0,240
Для решения задачи тепло- и массопереноса, включающей уравнения переноса импульса, энергии, массы, использовался пакет Сошбо1 МиШрЬувюБ. Одним из входных параметров модели являлась температура стенок парогазового канала. Для каждого значения температуры путем решения испарения, конденсации, диффузии, тепло- и массопе-реноса определялась концентрация легирующих элементов в расплаве,
давление и концентрация химических элементов в парах. Применение параметрического решателя СошБо1 МиШрЬувюБ для кластера, когда решение для каждого значения параметра температуры проводится на отдельном узле, позволило значительно сократить время расчетов. Для определения температуры парогазовой фазы в парогазовом канале результаты расчетов сопоставлялись с данными, полученными при экспериментальных исследованиях (рис. 2).
На рис. 2 представлены результаты расчетов химического состава парогазовой фазы в зоне ЭЛС (представлены результаты только для железа, хрома и марганца, так как содержанием остальных элементов можно пренебречь) и концентрации химических элементов, полученные экспериментально при рентгенофлуоресцентном анализе (горизонтальные линии).
............
2100 2150 2200 2250 2300 23502400 24502500 2550 2600 2650 27002750 Температура, К
Рис. 2. Зависимость концентрации железа, марганца и хрома на поверхности канала от температуры для стали 12Х18Н10Т: горизонтальными линиями обозначены экспериментальные значения концентраций; кривыми - расчетные
Значения температуры парогазовой фазы в канале определялись по пересечению расчетных и экспериментальных кривых концентрации элементов и составили 2320 К (для Fe), 2340 К (для Cr) и 2410 К (для Mn). Средняя арифметическая температура канала по этим эле-
ментам составляет 2360 К. Давление паров в канале при этой температуре - около 2 кПа (рис. 3).
На рис. 3 представлены кривые зависимости давления вышеуказанных элементов от температуры, полученные по разработанной математической модели, и зависимости, полученные в работах [9, 12] по равновесной модели, в соответствие с которой температура в канале составляет 2610 К. При этой температуре общее давление паров в канале составляет 11 кПа. Из графиков видно, что неравновесная модель дает меньшие значения температуры и давления паров в канале по сравнению с равновесной моделью.
Рис. 3. Расчетные зависимости давления паров железа, марганца и хрома в канале от температуры: сплошные линии - с учетом диффузионных процессов; штриховые - без учета диффузионных процессов
На основании проведенных исследований можно сделать следующие выводы:
1. Диффузионные процессы в свариваемом материале при ЭЛС в значительной мере влияют на концентрацию химических элементов в поверхностном слое расплава парогазового канала и, как следствие, в парах над зоной сварки. Результаты численного моделирования с учетом обеднения стенок канала легирующими элементами приводят к меньшим, по сравнению с равновесной моделью, температурам
(~2360 К). Указанное различие температур обусловливает значительно меньшие расчетные величины давлений паров в парогазовом канале (2кПа - при учете диффузионных процессов и 11кПа - без учета диффузии).
2. Установлено, что состояние пара в канале близко к насыщенному. Поток массы, уносимой из выходного сечения парогазового канала, составляет около 3 % от общего потока, испаряемого со стенок парогазового канала.
3. Разработанная модель позволяет рассчитать давления, температуру и концентрации химических элементов в парогазовой смеси над зоной сварки при ЭЛС. Модель также позволяет прогнозировать изменение химического состава сварного шва при ЭЛС.
Библиографический список
1. Ерофеев В. А., Логвинов Р.В., Нестеренков В.М. Особенности использования эквивалентного источника теплоты при учете деформаций и напряжений в процессе электронно-лучевой сварки // Сварка и диагностика. - 2010. - № 4. - С. 22-26.
2. Язовских В.М. Модель глубокого проплавления при электронно-лучевой сварке // Сварочное производство. - 2005. - № 8. - С. 16-20.
3. Numerical Modelling of the EBW process / V.A. Sudnik, V.A. Ero-feev, K.-H. Richter, K.-V. Heins // Proceedings of International Conference Computer Technology in Welding and Manufacturing & Mathematical Modeling and Information Techhnologies in Welding and Related Processes, Kiev, 6-8 June. - Kiev, 2000. - Р. 295-300.
4. Особенности структуры и свойств сварных соединений из толстолистовых сталей, выполненных электронно-лучевой сваркой / Т.П. Муравьева, В.К. Драгунов, А.П. Слива, А.Л. Гончаров // Сварочное производство. - 2010. - № 6. - С. 38-42.
5. Kaplan A.F.H., Norman P., Eriksson I. Analysis of the keyhole and weld pool dynamics by imaging evaluation and photodiode monitoring // Proceedings of LAMP2009 - the 5th International Congress on Laser Advanced Materials Processing. - 2009. - P. 1-6.
6. Numerical simulation of molten pool dynamics in high power disk laser welding / W.-I. Cho, S.-J. Na, C. Thomy, F. Vollertsen // Journal of Materials Processing Technology. - 2012. - № 212. - P. 262-275.
7. Development of a theory for alloying element losses during laser beam welding / U. Dilthey, A. Goumeniouk, V. Lopota, G. Turichin, E. Valdaitseva // Journal of physics d: applied physics. - 2001. - Vol. 34, № 1. - P. 81-86.
8. Ono Y., Matsumoto Sh. Diffusion of Chromium, Manganese, and Nickel in Molten Iron // Materials Transactions. - April 1974. - № 11. -P. 38-39.
9. Collur M.M., Paul A., Debroy T. Mechanism of Alloying element vaporization during laser welding // Metallurgical Transactions. - December 1987. - Vol. 18B. - Р. 733-746.
10. Саломатова Е.С., Трушников Д.Н., Беленький В.Я. Моделирование температуры в парогазовом канале при электронно-лучевой сварке // Тепловые процессы в технике. - 2013. - № 11. - С. 514-519.
11. Ландау Л. Д., Лифшиц Е.М. Теоретическая физика: учеб. пособие: в 10 т. Т. VI. Гидродинамика. - 3-е изд., перераб. - М.: Наука, 1986. - 736 с.
12. Саломатова Е.С., Беленький В.Я. Термодинамическая оценка влияния теплофизических свойств легирующих элементов на давление паров в канале проплавления при электронно-лучевой сварке высоколегированных сталей // Сварка и диагностика. - 2011. - № 2. - С. 22-25.
References
1. Erofeev VA., Logvinov R.V., Nesterenkov VM. Osobennosti is-polzovaniya ekvivalentnogo istochnika teploty pri uchete deformatsiy i napryazheniy v protsesse elektronno-luchevoy svarki [Features of the use of the equivalent heat source, taking into account the strain and stress in the process of electron beam welding]. Welding and Diagnostics, 2010, no. 4, pp. 22-26.
2. Yazovskikh VM. Model glubokogo proplavleniya pri elektronno-luchevoy svarke [Model of deep penetration in electron beam welding]. Welding Engineering, 2005, no. 8, pp. 16-20.
3. Sudnik V.A., Erofeev V.A., Richter K.-H., Heins K.-V. Numerical Modelling of the EBW process. Proceedings of International Conference Computer Technology in Welding and Manufacturing & Mathematical Modeling and Information Techhnologies in Welding and Related Processes, Kiev, 2000, 6-8 June, pp. 295-300.
4. Muraveva T.P., Dpagunov V.K., Sliva A.P., Goncharov A.L. Osobennosti stpuktupy i svoystv svarnykh soedineniy iz tolstolistovykh sta-ley, vypolnennykh elektponno-luchevoy svapkoy [Features structs and properties of welded joints of steel plates made electron-beam welding]. Welding Engineering, 2010, no. 6, pp. 38-42.
5. Kaplan A.F.H., Norman P., Eriksson I. Analysis of the keyhole and weld pool dynamics by imaging evaluation and photodiode monitoring. Proceedings of LAMP2009. The 5th International Congress on Laser Advanced Materials Processing, 2009, pp. 1-6.
6. 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, no. 212, pp. 262-275.
7. Dilthey U., Goumeniouk A., Lopota V., Turichin G., Valdaitseva E. Development of a theory for alloying element losses during laser beam welding. Journal of physics d. Applied physics, 2001, vol. 34, no. 1, pp. 81-86.
8. Ono Y., Matsumoto Sh. Diffusion of Chromium, Manganese, and Nickel in Molten Iron. Materials Transactions, April 1974, no. 11, pp. 38-39.
9. Collur M.M., Paul A., Debroy T. Mechanism of Alloying element vaporization during laser welding. Metallurgical Transactions, December 1987, vol. 18B, pp. 733-746.
10. Salomatova E.S., Trushnikov D.N., Belenkiy V.Ya. Modelirovanie temperatury v parogazovom kanale pri elektronno-luchevoy svarke [Modeling of temperature of vapor-gas channel in electron beam welding]. Thermal processes in engineering, 2013, no. 11, pp. 514-519.
11. Landau L.D., Lifshits E.M. Teoreticheskaya fizika [Theoretical Physics]. Moscow: Nauka, 1986. 736 p.
12. Salomatova E.S., Belenkiy V.Ya. Termodinamicheskaya otsenka vliyaniya teplofizicheskikh svoystv legiruyushchikh elementov na davlenie parov v kanale proplavleniya pri elektronno-luchevoy svarke vysokolegiro-vannykh staley [Thermodynamic evaluation of the influence of thermal properties of alloying elements on the vapor pressure in the channel penetration in electron beam welding of high-alloy steels]. Welding and Diagnostics, 2011, no. 2, pp. 22-25.
Об авторах
Саломатова Екатерина Сергеевна (Пермь, Россия) - ассистент кафедры «Сварочное производство и технология конструкционных материалов» ФГБОУ ВПО ПНИПУ (614990, г. Пермь, Комсомольский пр., 29, e-mail: [email protected]).
Трушников Дмитрий Николаевич (Пермь, Россия) - кандидат технических наук, доцент кафедры «Прикладная физика» ФГБОУ ВПО ПНИПУ (614990, г. Пермь, Комсомольский пр., 29, e-mail: [email protected]).
Цаплин Алексей Иванович (Пермь, Россия) - доктор технических наук, профессор, заведующий кафедрой «Общая физика» ФГБОУ ВПО ПНИПУ (614990, г. Пермь, Комсомольский пр., 29, e-mail: [email protected]).
Беленький Владимир Яковлевич (Пермь, Россия) - доктор технических наук, профессор кафедры «Сварочное производство и технология конструкционных материалов» ФГБОУ ВПО ПНИПУ (614990, г. Пермь, Комсомольский пр., 29, e-mail: [email protected]).
About the authors
Ekaterina S. Salomatova (Perm, Russian Federation) - Assistant Lecturer, Department of Welding Production and Technology of Constructional Materials, Perm National Research Polytechnic University (29, Kom-somolsky av., Perm, 614990, Russian Federation, e-mail: [email protected]).
Dmitriy N. Trushnikov (Perm, Russian Federation) - Ph. D. in Technical Sciences, Associate Professor, Department of Applied Physics, Perm National Research Polytechnic University (29, Komsomolsky av., Perm, 614990, Russian Federation, e-mail: [email protected])
Aleksey I. Tsaplin (Perm, Russian Federation) - Doctor of Technical Sciences, Professor, Department of General Physics, Perm National Research Polytechnic University (29, Komsomolsky av., Perm, 614990, Russian Federation, e-mail: [email protected]).
Vladimir Ya. Belenkiy (Perm, Russian Federation) - Doctor of Technical Sciences, Professor, Department of Welding Production and Technology of Constructional Materials, Perm National Research Polytechnic University (29, Komsomolsky av., Perm, 614990, Russian Federation, e-mail: [email protected]).
Получено 1.10.2014
