Построение тепловых моделей при электронно-лучевой сварке методом функций Грина Текст научной статьи по специальности «Физика»

ВЕСТНИК ПНИПУ

2017 Машиностроение, материаловедение Т. 19, № 3

Б01: 10.15593/2224-9877/2017.3.04 УДК 621.7

Т.В. Ольшанская1, Е.М. Федосеева1, Е.Г. Колева2

1 Пермский национальный исследовательский политехнический университет, Пермь, Россия 2 Химико-технологический и металлургический университет, София, Болгария

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

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

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

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

Установлено, что при решении тепловой задачи форма источника тепла должна имитировать «экранированное» состояние парогазового канала, что позволит при расчетах более точно передать «кинжальную» форму проплавления

По предложенной модели произведен расчет геометрии сварного шва в поперечном сечении для стали 38Х2НМЮА. Расчет показал, что предложенная модель источника тепла для электронно-лучевой сварки статическим лучом и полученное на ее основе решение тепловой задачи позволяют наиболее точно описывать геометрию сварного шва.

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

T.V. Olshanskaya1, E.M. Fedoseeva1, E.G. Koleva2

1 Perm National Research Polytechnic University, Perm, Russian Federation 2 University of Chemical Technology and Metallurgy, Sofia, Bulgaria

CREATION OF THERMAL MODELS AT ELECTRON BEAM WELDING BY METHOD OF FUNCTIONS OF GREEN

Development of new control systems of an electronic beam and programming of the modes of welding has allowed to expand considerably possibilities of improvement of the processing methods of electron beam welding providing almost unlimited variety of the modes of impact of a beam on a surface of the processed detail. On the basis of it, there is a need for development of the generalized system approach for the solution of thermal tasks at electron beam welding. The differential equation of heat conductivity presented in work, is mathematical model of the whole class of the phenomena of heat conductivity and has an infinite set of decisions. For receiving from a set of decisions of one private, characterizing concrete process, the additional data of of a condition of unambiguity containing body geometry, boundary and entry conditions are necessary. Thermal processes at electron beam welding are simulated and expression which is initial for creation of the thermal models considering different types of dynamic positioning of an electronic beam is received.

The possibility of application of a method of functions of Green as one of universal ways of creation of thermal models at electron beam welding with various dynamic positioning of a beam is presented.

On the basis of the provided analysis of the physical processes happening at electron beam welding the conclusion is made that formation of deep pro-melting represents a number of elementary cycles, each of which includes evaporation of a layer of metal and the subsequent shielding of an electronic beam.

It is established that at the solution of a thermal task the form of a source of heat has to imitate a "screened" condition of the steam-gas channel that will allow to transfer more precisely when calculating a "dagger" form of pro-melting

On the offered model calculation of geometry of a welded seam in cross section for steel 38H2NMYuA is made. Calculation has shown that the offered heat source model for electron beam welding by a static beam and the solution of a thermal task received on its basis, allow to describe geometry of a welded seam most precisely.

Keywords: electron beam welding, thermal tasks, form of a welded seam, thermal processes, beam oscillation, physical processes, channel of pro-melting, welding bathtub, welded seam, static beam.

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

Тепловые процессы при сварке чаше всего описываются с помощью дифференциального уравнения теплопроводности в подвижной системе координат с неподвижным источником теплоты:

Эх

Г д 2Т д 2Т д 2Т ^

лгдТ 1 жг,

+ + — ¥ (л, у, х, т), дх ср

(1)

= а

V

Эх2 + Эу2 + Эг2

+

где ¥(х, у, х, т) - функция источника теплоты.

Одним из методов решения задач теплопроводности является метод функций Грина (метод источников). Данный метод позволяет получать аналитические решения в зависимости от различных краевых задач и характера сварочного источника энергии, рассматривать закономерности изменения температурного поля в материале практически при любой продолжительности воздействия источника теплоты, свести многообразие способов нагрева материалов к некоторым схемам, охватывающим основные особенности процессов нагрева [1, 2].

Учитывая, что в последние годы, с развитием новых систем управления электронным лучом и программирования режимов сварки, значительно расширились возможности совершенствования технологических приемов электронно-лучевой сварки, обеспечивающих почти неограниченное разнообразие режимов воздействия луча на поверхность обрабатываемой детали (осцилляция луча, многофокусная и многолучевая сварка, комбинированные способы сварки и др.) [3], возникает необходимость в разработке обобщенного системного подхода для решения тепловых задач при электронно-лучевой сварке. В этом случае аналитический метод с помощью функций Грина является наиболее подходящим и универсальным.

Дифференциальное уравнение теплопроводности (1) является математической моделью целого класса явлений теплопроводности и имеет бесконечное множество решений. Чтобы получить из этого множества одно частное решение, характеризующее конкретный процесс, необходимо иметь дополнительные данные - условия однозначности, содержащие геометрию тела, граничные и начальные условия. Для моделирования тепловых процессов при электронно-лучевой сварке примем следующие условия однозначности:

- геометрия свариваемого тела - бесконечная пластина толщиной 8:

- граничные условия смешанного типа: по осям х и у граничные условия первого рода равны 0, по оси г граничные условия второго рода равны 0:

< х < < у < 0 < г <8;

дТ дТ

—(л, у,0, Г) = 0, — (л, у, 8, Г) = 0; дг дг

начальные условия Т (л, у, г,0) = 0.

Интегральное решение задачи теплопроводности (1) методом источников представляет собой совокупность стандартизованной функции Грина для конкретных геометрических и граничных условий и функции источника тепла:

Т (л, у, г, т) = {{{{ О (л, у, г, л', у', г', т) ехр

С V (л- х ) V2т^

4а

2а

х¥ (л, у, г, т)дл'ду'дг'дт, (2)

где О (л, у, г, л , у , г, т) - стандартизованная функция Грина; ¥ (л, у, г, т) - математическое описание функции источника тепла с помощью дельта-функции Дирака; л', у', г - координаты источника тепла.

Решение конкретной тепловой задачи сводится к нахождению стандартизованной функции Грина с учетом граничных условий и определению математического описания функции для заданного источника тепла. Известно, что функция Грина допускает неполное разделение переменных (она разделяется по пространственным переменным л, у, г, но не разделяется по времени т), т.е. может быть представлена в виде произведения одномерных функций Грина, выбираемых исходя из краевых условий:

О (л, у, г, л', у', г', т) = Ох (л, л, т) Оу (у, у', т) О2 (г, г ,т).

(3)

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

Ол (л л', т):

1

24 пат Оу (у,у',т) = 1

ехр

С ( ' , Т7 \2 ^

(л - л + Vт) 4ат

2

пат

г ехр

С I '\2 ^

(у - у )

4ат

(4)

(5)

1

^ (г,т) = :тт=тх

2\] па т

Г

■Т

Г (

ехр

(г - г' + 2п8) 4ат

2Л

+ ехр

(г + г' + 2 п8) 4ат

Л ^

(6)

у У

Таким образом, в общем виде интегральное решение задачи теплопроводности в подвижной системе координат (2) для бесконечной пластины будет следующим:

Т (x, ^ г, т) = ЦЦ

1

8 (пат |

-ехр

Г / ' ,„ \2 ^ (х - х + Ут)

4ат

ехр

( /\2 ^ (у- у )

4ат

Г

ехр

(г - г' + 2п8) 4ат

2

+ ехр

(г + г' + 2п8)

4ат

У У

х ¥ (х, у, г, т)Эх'Эу'Эг'Эт.

(7)

Полученное выражение (7) является исходным для построения тепловых моделей, учитывающих различные виды динамического позиционирования электронного луча. Задачей построения решения тепловых задач становится определение вида источника тепла для каждого конкретного случая и его математического выражения ¥ (х, у, г, т).

В работах Н.Н. Рыкалина предложено распределение температур при электронно-лучевой сварке с глубоким проплавлением рассматривать как результат действия двух источников тепла: точечного на поверхности и линейного, ограниченного по глубине величиной пластины [1, 4]. При тепловой мощности луча Q мощность точечного источника принимается q1 = kQ, линейного д2 = (1 - k)Q, где к - коэффициент распределение энергии. Решение тепловой задачи аналитическим методом для электронно-лучевой сварки с применением комбинированного источника тепла нашло продолжение в работе [2]. Предложенная тепловая модель для бесконечной пластины строится методом функций Грина, источник тепла является нормально распределенным по поверхности и линейным по глубине (рис. 1).

п = — со

п=—со

Рис. 1. Схема источника тепла при электронно-лучевой сварке [2]

Математическое описание функции источника тепла в данном случае имеет следующий вид:

р (*г) = ф {^)8(/)8(г')Е (т)+к 8(/Ж/ )8(г' )е м),(8)

Е (')= I1 пРи 0 * г' * к Е (т) = I1 пРи *0

У ' |0 при к < г < 0, 1 ; |о при т>',

0 4аК ' К сС2'

где д - мощность электронного луча; п - коэффициент полезного действия; с - удельная теплоемкость; р - плотность металла; к1 и к2 - коэффициенты, учитывающие распределение мощности луча между поверхностным источником и линейным; к - глубина линейного источника; Е(г') и Е(т) - единичные функции. Для имитации воздействия нормально-кругового источника рассчитывается время действия фиктивного источника '0 через коэффициент сосредоточения К для заданного диаметра электронного луча с.

Интегральное решение задачи теплопроводности для электроннолучевой сварки статическим лучом бесконечной пластины толщиной 8 имеет вид [2]

Т (х, у, г, т) =

кт

4ср((а) ^ ((Г)

-ехр

2 >

ехр

(г + 2п8) 4ат

2Л

С т +

к2 дп г1

8кпХ^ т

'0

I _ехр

•> т

(х + V т)2 + у 4ат

Г (х + Vт)2 + у2 Л 4ат

п=—со

Автором работы [2] установлено, что при мощности электронного луча от 5 до 10 кВт для легированных сталей среднее значение коэффициентов к1 и к2 составляет 0,3 и 0,7 соответственно, при коэффициенте полезного действия п = 0,9. Однако последующие многочисленные расчеты для легированных сталей и цветных металлов показали, что тепловая модель для электронно-лучевой сварки статическим лучом (9) не всегда дает удовлетворительные результаты при сравнении геометрии сварного шва, полученной рассчетно, с экспериментальной.

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

Глубокое проплавление при электронно-лучевой сварке сопровождается достаточно сложными тепловыми, гидро- и газодинамическими, плазменными и другими процессами [1, 4, 5]. При воздействии концентрированного потока энергии на поверхность тела с удельной мощностью выше критической 109-1010 Вт/м2 часть потока энергии частично отражается от поверхности, а остальная часть поглощается в тонком поверхностном слое металла, вызывая его нагрев, последующие плавление и испарение. Все существующие гипотезы, объясняющие феномен «кинжального» проплавления, исходят из единого мнения, что глубокое проникновение электронов в металл вызвано образованием в жидком металле сварочной ванны узкого и глубокого канала проплавления (парогазового канала) [6].

Формирование парогазового канала сопровождается достаточно сложными процессами, такими как локальный перегрев металла, стохастическое перемещение зоны взаимодействия луча с металлом по глубине канала проплавления, периодическое смыкание канала в верхней части, образование волн на поверхности расплава и др. [7]. Теоретически и экспериментально установлено, что канал проплавления и окружающий его жидкий металл имеют переменную во времени форму, а процесс формирования глубокого проплавления при электронно-

(9)

лучевой сварке имеет автоколебательный характер [1, 8]. Было выдвинуто множество гипотез об автоколебательных процессах в канале проплавления. Многие из этих гипотез не получили однозначного подтверждения или опровержения до сих пор. Приведем лишь некоторые значимые из них.

В работах коллектива академика Н.Н. Рыкалина выдвигается гипотеза о взрывном характере взаимодействия электронного луча с металлом, обусловленном большими скоростями ввода энергии при электронно-лучевой сварке, превышающими скорость отвода тепла за счет теплопроводности [9-11]. Согласно этой теории процесс углубления канала проплавления в металле является прерывистым. В результате периодического взрывного вскипания металла давление пара над поверхностью и в канале возрастает. При некоторой критической плотно*

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

Периодическое смыкание канала проплавления может происходить и в результате автоколебательных процессов в жидком металле сварочной ванны, связанных с кипением и конденсацией паров металла в канале [12-14], частота колебаний составляет порядка 1-100 Гц. Движение расплавленного металла по стенкам канала имеет периодический характер и осуществляется объемными волнами под действием силы гидростатического давления и поверхностного натяжения, противодействующих реакции отдачи паров. При этом скатывающаяся на

дно канала волна оголяет твердые слои металла на передней стенке, после чего процессы плавления и переноса повторяются [15, 16]. Зона интенсивного испарения сосредоточивается на дне канала, боковая поверхность которого выводится из зоны действия луча при заглублении в металл. Боковые стенки парогазового канала удерживаются за счет реакции отдачи паров. При этом стенки канала охлаждаются за счет теплоотвода в твердый металл, и поток пара конденсируется на поверхности стенок канала. В момент, когда в верхней части канала давление отдачи паров будет минимальным за счет конденсации и меньше суммы капиллярного и гравитационного давления жидкости, может происходить стягивание входного отверстия силами поверхностного натяжения (смыкание канала проплавления жидкостью) [12-14]. В момент смыкания канала жидкостью также наблюдается экранирование электронного луча.

а б в

Рис. 2. Распределение жидкой фазы в канале по мере увеличения его глубины: а - р > (ро + ра); б - р > (ро + ра); в - р < (ро + ра); 1 - жидкая фаза; 2 - пар; р - давление пара в канале; ро - гидростатическое давление жидкого металла; ра - давление поверхностного натяжения [4]

Какими бы причинами ни вызывалось периодическое, частичное или полное экранирование луча, оно приводит к рассеиванию электронного луча на стенки, которое может происходить близко у поверхности на какой-то глубине. Удельная мощность луча в момент его рассеяния на стенки падает на 1-2 порядка (пропорционально площади внутренней поверхности канала), поэтому обеспечивается практически только плавление кромок [4, 7], приводящее к формированию расширенной верхней части. Также определенное влияние на формирование

расширенной верхней части геометрии проплавления может оказывать эффект Марангони, как показано в работах [17, 18], где численным экспериментом продемонстрированы решения для обтекания стационарного канала проплавления при электронно-лучевой сварке потоком жидкости при различных режимах течения.

На основе приведенного анализа физических процессов, происходящих при электронно-лучевой сварке, можно сделать заключение, что глубокое проплавление представляет собой ряд элементарных циклов, каждый из которых включает испарение слоя металла и последующее экранирование электронного луча. Таким образом, при автоколебательных процессах проплавления парогазовый канал может находиться в двух квазипериодических состояниях: первое состояние соответствует открытому каналу, когда давление паров металла значительно превышает сумму гидростатического давления и поверхностного натяжения (рис. 2, а); второе соответствует экранированию электронного луча, сопровождающемуся уменьшением удельной мощности луча и рассеиванием его на стенки канала (рис. 2, б).

При решении тепловых задач электронно-лучевой сварки комбинированные источники тепла, предложенные в работах Н.Н. Рыкалина [1, 4] и работе В.М. Язовских [2] (см. рис. 1, формула (9)), имитируют первое - «открытое» состояние парогазового канала с учетом расширения его верхней части. Следует отметить, что в работах Н.Н. Рыка-лина представлен расчет времени образования канала проплавления как суммарное время циклов «чистого испарения» и экранирования луча. Результаты расчетов показали, что в общем балансе времени образования канала на процесс «чистого испарения» уходит около 5-10 % времени, тогда как рассеивание занимает остальные 90-95 % времени. Исходя из этого при решении тепловой задачи форма источника тепла должна имитировать второе - «экранированное» состояние парогазового канала, что позволит при расчетах более точно передать «кинжальную» форму проплавления. Можно выделить два основных момента, позволяющих описать подобный источник тепла:

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

2. Экранирование происходит не на самой поверхности, а на некоторой глубине от нее, и одновременно в канале проплавления продолжают происходить процессы, вызывающие нагрев основного металла (конденсация пара и теплоотвод в твердый металл). В связи с этим линейный по глубине источник тепла может быть представлен как действующий под поверхностью на каком-то расстоянии.

Таким образом, глубокое проплавление при электронно-лучевой сварке можно представить как результат действия комбинированного источника тепла, непрерывно действующего в течение определенного отрезка времени t (рис. 3) и состоящего:

- из нормально распределенного на поверхности кругового источника радиусом г, вводимого в начале координат;

- линейного по глубине источника (вдоль оси Z, длиной 2Н), действующего под поверхностью на расстоянии Н1, вводимого в координате (0, 0, 5).

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

Рис. 3. Схема комбинированного источника тепла при электронно-лучевой сварке статическим лучом: г - радиус поверхностного источника тепла; Н1 и Н2 - глубина залегания и действия линейного источника тепла соответственно; 5 - точка ввода линейного источника тепла длиной 2Н, где Н = (Н2- НО/2,

5 = Н + Й1

Математическое описание функции источника тепла (см. рис. 3) имеет следующий вид:

^(х, у, г, т)

Е(х')В(у)8(г)Е(т) + |8(х')8(/)£(г)Е(т)], (10)

Г1 ПРИ -г < X < г, ' Г1 при -г < у' < г, Е (х ) = 1 , Е (у ) = 1 ,

10 при г < х < - г, 10 при г < у < -г,

Г1 при -Н < г < Н, Г1 при^ <т<',

Е(г ) = 1 , Е(т) = 1 о

у ' [0 при Н < г <-Н, у / [0 при т>', 0 4аК

'0 ='

К =

й2

Для получения интегрального решения тепловой задачи подставляем функцию источника (10) в уравнение (7), в котором уже учтены краевые условия для бесконечной пластины в подвижной системе координат. Интегральное решение тепловой задачи при электроннолучевой сварке статическим лучом в аналитической форме имеет следующий вид:

Т (х, у, г, т) =

\±

16г2фл/Ла * л/т

кт

ей"

у + г 2л/0т,

■erf

у - г

г(х+г+Ут] х-г+Ут' егг -— I - егг -—

ч 2у1 ат У ^ 24 ат _

2 ехр - ' дт +

V

4ат

у

+

к2?П Г 1

/ -

т

16НяА т

'0

ехр

(х + У т)2 + у

2 ]

у

2

- ей"

ей"

г - 5 + Н + 2п8^ 2л/ат ,

+ erf

4ат

( г + 5 + Н + 2п8

V

2л/ат

( г - 5 - Н + 2п8 ] г + 5 - Н + 2п8 ]

2>/ат

- erf

2л/ат

Эт.

(11)

Полученное выражение (11) может быть хорошо реализовано в математическом пакете Ма&САБ. Для этого необходимы следующие исходные данные: теплофизические свойства свариваемого материала, параметры режима сварки и параметры источника тепла (см. рис. 3): к1 и к2 - коэффициенты, учитывающие распределение мощности электронного луча между поверхностным источником и линейным; г - ра-

диус поверхностного источника тепла; 2к — длина линейного источника; 5 - координата ввода линейного источника. Из схематического изображения источника тепла видно, что для определения к и 5 необходимо знать расстояние от поверхности кь на которое заглублен источник, действующий по оси Z, и глубину действия этого источника к2, представляющую собой глубину проплавления.

Судить о глубине линейного источника к2, действующего по оси Z, можно по расчетной глубине проплавления, которая связана с параметрами электронно-лучевой сварки критериальным уравнением. В настоящие время существует большое количество работ, посвященных определению глубины проплавления. Ниже приведены некоторые наиболее значимые математические модели для определения глубины проплавления при сварке статическим лучом:

1-0,5*

Н =(пи )Ч°В5* -1'

(2гл )-°,5к, к = 0,68(ХГпл )°'15 [19], (12)

0,15

Н =

Ш

г V

л св

Г 5ХТср

срТср+р^пЛ +

5ХТ Л

-1

+

ср

2ГлVсв у

[6],

Н = Ш ^

' 2,24Х ^

V ф^в Г у

(( - Т° )) [6]

(13)

(14)

Н

пэи1 - 4г 2еаТс

ср

г V

л св

фТпл + Р^пл + 5АТШ

2а

- + -

2ГлVсв у

[6],

(15)

Н

Ш

2Д^св р(Т + ^ )1 + 2,5^

1 1

— + —

[4],

а

г V

л св у

(16)

где Н - глубина проплавления, рассчитанная по соответствующим формулам, см; и - ускоряющее напряжение, В; I - ток луча, А; Пэ -эффективное КПД; Vcв - скорость сварки, см/с; гл - радиус электронного луча, см; с - теплоемкость, Дж/(г-°С; р - плотность, г/см3; X - коэффициент теплопроводности, Дж/(см-с-К); а - коэффициент температуропроводности, см2/с; Тпл - температура плавления, К; Ткип - темпера-

1

пл

тура кипения, К; Тср - средняя температура сварочной ванны, К, Тср = 1,5Тпл; Т0 - начальная температура, К; Ьпл - скрытая теплота плав-

3 3

ления, Дж/см ; Ькип - скрытая теплота парообразования, Дж/см ; 8 - коэффициент черноты, 8 = 0,93; £ - доля металла, испаренного из сварочной ванны, £ = 0,01; а - постоянная Стефана-Больцмана, а = 5,6810-8 Вт-м-2-К-4 .

Формулы (12)-(16) могут быть использованы и для предварительного определения основных параметров режима сварки (тока луча, ускоряющего напряжения и скорости сварки) для заданной глубины проплавления.

Параметр Н1 - расстояние от поверхности, на которое заглублен источник, действующий по оси в формуле (11) определяется следующим образом. Известно, что при воздействии электронного луча на материал распределение плотности мощности луча близко к нормальному закону Гаусса по поверхности и по объему, с максимумом под поверхностью нагреваемого тела (рис. 4).

Рис. 4. Схема теплового источника при электронно-лучевом воздействии: 1 - распределение плотности мощности луча по поверхности тела;

2 - по объему [4]

Как предложено в работе [4], глубина максимального энерговыделения гтах принимается

гтах = 0,758, (17)

где 8 - максимальная глубина проникновения электронов в материал, определяемая по формуле Шонланда:

ю и2

8 = 2,35-10-12—, (18)

Р

где и - ускоряющие напряжение, В; р - плотность материала, кг/м .

В первом приближении параметр к1 в формуле (11) может соответствовать глубине максимального энерговыделения гшах и определяться из формул (17) и (18).

Радиус поверхностного источника тепла г в формуле (11) (см. рис. 3) определяется из условия экранирования электронного луча. Поскольку в момент рассеивания электронного луча на стенки канала удельная мощность электронного луча (д2 = Ш/лг2) уменьшается в к раз (на 1-2 порядка), а вводимая мощность луча не изменяется (д = и!), то увеличение размера зоны в канале проплавления, на которую приходится рассеивание электронного луча, будет определяться как

г = лкл, (19)

где г - радиус «рассеянного» луча, соответствующий радиусу поверхностного источника тепла в формуле (11); к изменяется в пределах 10-100 раз, конструктивно для сталей может быть взят в пределах 50-60; гл - радиус электронного луча.

На рис. 5 представлен результат расчета распределения температур по формуле (11) для стали 38Х2НМЮА, режим сварки: и = 28 кВ; 1луча = 0,2 А; vcв - 16 м/ч; свариваемая толщина 20 мм. Распределение температур приведено в трех координатных плоскостях, при этом в продольном сечении показаны отдельно графики для поверхностного (рис. 5, а) и линейного по глубине источников тепла (рис. 5, б) и итоговое распределение температур (рис. 5, в).

Для сравнения был произведен расчет геометрии сварного шва в поперечном сечении для стали 38Х2НМЮА по формуле (9), где источник тепла поверхностный нормально распределенный и линейный по глубине, и по формуле (11) с предложенным разделенным источником тепла (рис. 6). Геометрия сварного шва в поперечном сечении определялась по изотерме кристаллизации, точками показаны экспериментально измеренные значения ширины шва. Из рисунка видно, что предложенная модель источника тепла для электронно-лучевой сварки статическим лучом (3) и полученное на ее основе решение тепловой задачи (11) позволяют наиболее точно описывать геометрию сварного шва.

О -2-4-6 в

-10 -8-6^-2 0 2 4 6 д

У, мм

4 2 0 -2 -4 -6

г

Рис. 5. Распределение температур при электронно-лучевой сварке стали 38Х2НМЮА: плоскости XOZ (а, б, в); плоскость ХОУ на поверхности (г); плоскость УOZ (д); а, г, д - комбинированный источник тепла; б - действие поверхностного источника тепла; в - линейный по глубине источник

а б в

Рис. 6. Геометрия (а, б) и макроструктура (в) сварного шва стали 38Х2НМЮА, режим сварки: и = 28 кВ; /луча = 0,2 А; усв = 16 м/ч; а - расчет по формуле (9): к1 = 0,3, к2 = 0,7, Н = 20 мм; б - расчет по формуле (11): к1 = 0,3, к2 = 0,7, Н1 = 0,4 мм, Н2 = 20 мм, г - 3,7 мм

При электронно-лучевой сварке осцилляция луча по траекториям различного вида в большинстве случаев применяется для устранения характерных дефектов, имеющих место при сварке статическим лучом. Наиболее широко используются следующие виды развертки луча: продольная и поперечная, х-образная, а также перемещения луча по эллипсу, дуге и окружности. Амплитуда осцилляции луча чаще всего лежит в пределах 1-3 мм, частота - от 50 Гц до 1 кГц и зависит от вида свариваемого материала. При достаточно малых значения частоты и амплитуды осцилляции луча размеры парогазового канала практически не изменяются, при более высоких значениях происходит его расширение. Осцилляция электронного луча оказывает воздействие на паро-и гидродинамические процессы в канале проплавления, повышая его устойчивость, и, как следствие, приводит к изменению конфигурации сварочной ванны. Однако в форме поперечного сечения сварного шва сохраняется расширение в верхней части, т.е. можно предположить, что при данных режимах осцилляции луча процессы, связанные с экранированием электронного луча, сохраняются, хотя и происходят в меньшей степени.

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

Обычно применяемая на практике частота осцилляции луча для различных видов развертки выбирается достаточно большой (от 400 Гц), поэтому источник тепла можно рассматривать как постоянно действующий какой-то определенной формы.

Для построения тепловых моделей выбраны продольная, поперечная и х-образная траектории осцилляции луча. При сварке с продольными и поперечными колебаниями электронного луча источник тепла состоит из поверхностного линейного, действующего вдоль соответствующих осей X и У, и плоского, действующего по глубине вдоль оси Z по ширине вдоль соответствующих осей X и У (рис. 7, а, б). Осцилляция луча по х-образной траектории представляет собой суперпозицию продольных и поперечных колебаний, поэтому источник имеет прямоугольную форму (рис. 7, в).

а

б

Математическое описание источника тепла: - поперечные колебания луча (см. рис. 7, а)

^ (х, у, г, т) = ^ Ф

С

к

4гАх

Е (х')Е1 (у'Щг')В (т) +

+ •

4 АН

8(х')Е(у')Е(г )Е(т)1;

продольные колебания луча (см. рис. 7, б)

^ (х, у, г, т) = ^ ср

(

к

V 4гВх

Е (X)Е(У)8(г')Е(т) +

+ •

4ВН

Е (X )8(у )Е (г' )Е (т) 1;

■ х-образные колебания луча (см. рис. 7, в)

^ (X, у, г, т) = М Ф

к

4 А! В1

Е (X)Е (у")8(г')Е(т) +

+

8 АВН

Е (X )Е (у' )Е (г' )Е (т)1,

где Е1 (х ) =

(1 при -В1 < х < В1, , (1 при -А1 < у' < А1

0 при В1 < х < -В1,

Е1 (у' ) =

1

0 при А1 < у <- А1;

Е ) = (1 при -В < х' < В, (х ) V0 при В < х' <-В,

(1 при -Н < г < Н,

е (г )=

0 при Н < г < -Н,

Е (у' ) Е (т) =

(1 при -А < у' < А,

0 при А < у < - А;

'1 приг0 <т< t, 0 при т> t,

(20)

(21)

(22)

где А, В - амплитуды колебаний источника вдоль осей Г и Х, действующего на глубине; А1, В1 - амплитуды колебаний поверхностного источника (В1 = В + г; А1 = А + г); г - радиус поверхностного источника.

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

предлагается использовать прямоугольную форму источника и вводить действие фиктивного источника со временем ¿0. Благородя этому форма пятна нагрева становится близкой к круговой, а распределение мощности по пятну нагрева отличается значительно меньшей кривизной по сравнению с нормально-круговым источником тепла.

Подставив функции источников тепла (20)-(22) в уравнение общего вида интегрального решения задачи теплопроводности в подвижной системе координат для бесконечной пластины (7), получаем следующие выражения интегрального решения задач теплопроводности для электронно-лучевой сварки с осцилляцией луча:

- с продольными колебаниями электронного луча

Т (х, у, г, т) =

ей"

кт

г ±

16гВ1ср4ла; 4т

(

ей"

х + В + V т

(

+

у + г 240т,

-Г -

32Вкср4ш 14т * Ё

- ей"

240т ^

240т

Ё ехР

-erf

х - В + V т

(

ей"

х + В + V т 240т ,

- erf

erf

-erf

( г - 5 + к + 2п84 . 240т ,

5 - к + 2п8 ^ 240т

+ erf

( (г + 2п8)2 ^

4ат

ч

( х - В + V т ^

V 2л/0т у

( г + 5 + к + 2п8

240т Эт +

ехр

у

2

4ат

-erf

240т

( г + 5 - к + 2п8 ^ 240т

с поперечными колебаниями электронного луча

Т (х, у, г, т):

erf

\6rAfp4na; 4т

Г_1

К т

(

erf

х + г + V т

(

+

к2ЯП

32Акср4т ' 4т

У + А! 240т

- ей"

(у-Л'

240т,

ч 240т Ё ехр

- ей"

( (г + 2п8)2 ^

4ат

Эт;

х - г + V т 240т ,

Эт +

(23)

ехр

( (х + Vт)2 ^ 4ат

erf

у + А 240т у

-ей"

( У - А'

240т

Т

ей"

2у[ах

( г - 5 - Н + 2п8 2>/От

5 + Н + 2п8 ^ _( г + 5 + Н + 2п8^ - + erf -

■erf

2Тот

( г + 5 - Н + 2п8 2у[ах ,

Эх;

(24)

с х-образными колебаниями электронного луча

Т (х, у, г, т) =

erf

кт

16гА1Б1ср4ла; >/т

■к /г

У + А

г

ей"

х + Б + У т

2>/От

- ей"

х-

Б1 + У т

2>/ат

х

- ей"

У - А 2>/От.

Т ехР

(г + 2п8)

2Л

4ат

Эт +

+

к2?П

64 АБНср

Г

erf

х + Б + У т

С

2л/ах

-erf

х - Б + У т 2Тат ,

erf

у + А 24ог,

- erf

У - А 2л/ах,

- erf

Т

erf

■5 + Н + 2п8 ^ _( г + 5 + Н + 2п8^ ■ + ей7 ■

24аГ

24аГ

( г - 5 - Н + 2п8 2л/ат ,

- erf

( г + 5 - Н + 2п84 24аГ ,

Эт.

(25)

Полученное решение для сварки с х-образной осцилляцией электронного луча может быть использовано в случае электронно-лучевой сварки с разверткой по круговой и эллиптической траекториям. На рис. 8 продемонстрировано применение формулы (25) для определения формы проплавления (сварного шва) и построения термических циклов при электронно-лучевой сварке с х-образной осцилляцией луча. Сварка проводилась на стали 35Х2НМ толщиной 20 мм при мощности электронного луча 5600 Вт, скорость сварки составляла 10 м/ч, ток фокусировки обеспечивал максимальную глубину проплавления, частота х-образных колебаний 400 Гц, амплитуда колебаний 1,8 мм. Термические циклы сварки представлены для зоны перегрева (Ттах = 1300 °С) для различных участков по глубине сварного шва. Расчет производился в математическом пакете МаЛСАБ.

в

Рис. 8. Макроструктура сварного шва (а), расчетная форма сварного шва (б) и термические циклы сварки (в) для зоны перегрева различных участков по глубине сварного шва, сталь 35Х2НМ, электронно-лучевая сварка с х-образной

осцилляцией луча

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

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

1. Лазерная и электронно-лучевая обработка материалов: справочник / Н.Н. Рыкалин, А.А. Углов, И.В. Зуев, А.Н. Кокора. - М.: Машиностроение, 1985. -496 с.

2. Язовских В.М. Математическое моделирование и инженерные методы расчета в сварке: в 2 ч. Ч. 2. Тепловые процессы при сварке и моделирование в пакете MathCAD. - Пермь: Изд-во Перм. гос. техн. ун-та, 2008. - 119 с.

3. Ольшанская Т.В., Саломатова Е.С. Обзор современных способов управления электронным лучом при электронно-лучевой сварке // Вестник Пермского национального исследовательского политехнического университета. Машиностроение, материаловедение. - 2016. - Т. 18, № 4. - С. 169-187.

4. Рыкалин Н.Н., Углов А.А., Зуев И.В. Основы электронно-лучевой обработки материалов. - М.: Машиностроение, 1978. - 239 с.

5. Назаренко O.K., Кайдалов А.А. Электронно-лучевая сварка. - Киев: Наук. думка, 1987. - 256 с.

6. Кайдалов А.А. Электронно-лучевая сварка и смежные технологии. -2-е изд., перераб. и доп. - Киев: Экотехнология, 2004. - 260 с.

7. Kaplan A.F.H., Normanand P., Eriksson I. Analysis of the keyhole and weld pool dynamics by imaginge valuation and photodiode monitoring // Proceedings of LAMP 2009 -the 5th International Congress on Laser Advanced Materials Processing. - 2009. - Р. 1-6.

8. Автоколебательный режим нагрева металла при ЭЛС / Н.Н. Рыкалин, А.А. Углов, И.В. Зуев, В.И. Скобельник, С.В. Селищев // Электронно-лучевая сварка. -М.: МДНТП, 1993. - С. 77-82.

9. Рыкалин Н.Н., Зуев И.В., Углов А.А. Оценка глубины проплавления при электронно-лучевой сварке // Физика и химия обработки материалов. - 1972. - № 1. -С. 9-14.

10. Зуев И.В., Рыкалин Н.Н., Углов А.А. О колебаниях глубины проплавления при электронно-лучевой сварке // Физика и химия обработки материалов. - 1975. -№ 1. - С. 136-141.

11. Рыкалин Н.Н., Углов А.А., Зуев И.В. Механизм сварки и обработки электронным лучом // Сварка электронным лучом. - М.: МДНТП, 1974. - С. 13-19.

12. Формирование эквивалентного источника теплоты для расчетов деформаций конструкций на основе теоретического моделирования воздействия электронного луча на металл / В.А. Ерофеев, Р.В. Логвинов, В.М. Нестеренков, В.В. Плошихин // Известия Тул. гос. ун-та. Технические науки. - 2008. - Вып. 4. - С. 155-166.

13. Васильев А.А., Ерофеев В.А., Судник В.А. Возникновение колебаний глубины парогазового канала при лучевых способах сварки // Известия Тул. гос. ун-та. Технические науки. - 2014. - Вып. 7. - С. 130-140.

14. Ерофеев В.А. Моделирование особенностей формирования парогазового канала при электронно-лучевой сварке // Сварка и диагностика. - 2009. - № 4. -С. 12-18.

15. Назаренко O.K., Кайдалов А.А., Акопьянц К.С. Периодический перенос металла при электронно-лучевой сварке // Автоматическая сварка. - 1973. - № 8. -С. 75-76.

16. Назаренко О.К., Кайдалов А.А., Акопьянц К.С. О нестабильности глубины проплавления при электронно-лучевой сварке // Сварка электронным лучом. - М.: МДНТП, 1974. - С. 26-30.

17. Cho J.H., Na S.J. Three-dimensional analysis of molten pool in GMA-laser hybrid welding // Welding Journal. - 2009. - № 2. - P. 35-43.

18. Heat transfer and fluid flow during electron beam welding of 304L strainless steel alloy / R. Rai, T.A. Palmer, J.W. Elmer, T. Debroy // Welding Journal. - 2009. - № 3. -P. 54-61.

19. Ерохин А.А. К расчету глубины проплавления при электронно-лучевой сварке // Сварочное производство. - 1976. - № 4. - С. 4-6.

References

1. Rykalin N.N., Uglov A.A., Zuev I.V., Kokora A.N. Lazernaia i elektronno-luchevaia obrabotka materialov [Laser and electron beam processing of materials]: spra-vochnik. Moscow: Mashinostroenie, 1985, 496 p.

2. Iazovskikh V.M. Matematicheskoe modelirovanie i inzhenernye metody rascheta v svarke [Mathematical modeling and engineering methods of calculation in welding]. V 2 chastiakh. Chast' 2. Teplovye protsessy pri svarke i modelirovanie v pakete MathCAD. Perm': Izdatel'stvo Permskogo gosudarstvennogo tekhnicheskogo universiteta, 2008, 119 p.

3. Ol'shanskaia T.V., Salomatova E.S. Obzor sovremennykh sposobov upravleniia elektronnym luchom pri elektronno-luchevoi svarke [The review of modern ways of management of an electronic beam at electron beam welding]. Vestnik Permskogo natsion-al'nogo issledovatel'skogo politekhnicheskogo universiteta. Mashinostroenie, materialove-denie, 2016, vol. 18, no. 4, pp. 169-187.

4. Rykalin N.N., Uglov A.A., Zuev I.V. Osnovy elektronno-luchevoi obrabotki materialov [Bases of electron beam processing of materials]. Moscow: Mashinostroenie, 1978, 239 p.

5. Nazarenko O.K., Kaidalov A.A. Elektronno-luchevaia svarka [Electron beam welding]. Kiev: Naukova dumka, 1987, 256 p.

6. Kaidalov A.A. Elektronno-luchevaia svarka i smezhnye tekhnologii [Electron beam welding and adjacent technologies]. 2-e izdanie, pererabotannoe i dopolnennoe. Kiev: Ekotekhnologiia, 2004, 260 p.

7. Kaplan A.F.H., Normanand P., Eriksson I. Analysis of the keyhole and weld pool dynamics by imaginge valuation and photodiode monitoring. Proceedings of LAMP 2009 -the 5th International Congress on Laser Advanced Materials Processing, 2009, pp. 1-6.

8. Rykalin N.N., Uglov A.A., Zuev I.V., Skobel'nik V.I., Selishchev S.V. Avto-kolebatel'nyi rezhim nagreva metalla pri ELS [The self-oscillatory mode of heating of metal at electron-beam weldings]. Elektronno-luchevaia svarka. Moscow: MDNTP, 1993, pp. 77-82.

9. Rykalin N.N., Zuev I.V., Uglov A.A. Otsenka glubiny proplavleniia pri elek-tronno-luchevoi svarke [Pro-melting depth assessment at electron beam welding]. Fizika i khimiia obrabotki materialov, 1972, no. 1, pp. 9-14.

10. Zuev I.V., Rykalin N.N., Uglov A.A. O kolebaniiakh glubiny proplavleniia pri elektronno-luchevoi svarke [About fluctuations of depth of pro-melting at electron beam welding]. Fizika i khimiia obrabotki materialov, 1975, no. 1, pp. 136-141.

11. Rykalin N.N., Uglov A.A., Zuev I.V. Mekhanizm svarki i ob-rabotki elektron-nym luchom [Mechanism of welding and processing by an electronic beam]. Svarka elektronnym luchom, Moscow: MDNTP, 1974, pp. 13-19.

12. Erofeev V.A., Logvinov R.V., Nesterenkov V.M., Ploshikhin V.V. Formirovanie ekvivalentnogo istochnika teploty dlia raschetov deformatsii konstruktsii na osnove teo-reticheskogo modelirovaniia vozdeistviia elektronnogo lucha na metall [Formation of an equivalent source of warmth for calculations of deformations of designs on the basis of theoretical modeling of impact of an electronic beam on metal]. Izvestiia Tull'skogo gosu-darstvennogo. universiteta. Tekhnicheskie nauki, 2008, iss. 4, pp. 155-166.

13. Vasil'ev A.A., Erofeev V.A., Sudnik V.A. Vozniknovenie kolebanii glubiny parogazovogo kanala pri luchevykh sposobakh svarki [Emergence of fluctuations of depth of the steam-gas channel at beam ways of welding]. Izvestiia Tul'skogo gosudarstvennogo universiteta. Tekhnicheskie nauki, 2014, iss. 7, pp. 130-140.

14. Erofeev V.A. Modelirovanie osobennostei formirovaniia parogazovogo kanala pri elektronno-luchevoi svarke [Modeling of features of formation of the steam-gas channel at electron beam welding]. Svarka i diagnostika, 2009, no. 4, pp. 12-18.

15. Nazarenko O.K., Kaidalov A.A., Akop'iants K.S. Periodiche-skii perenos met-alla pri elektronno-luchevoi svarke [Periodic transfer of metal at electron beam welding]. Avtomaticheskaia svarka, 1973, no. 8, pp. 75-76.

16. Nazarenko O.K., Kaidalov A.A., Akop'iants K.S. O nestabil'nosti glubiny proplavleniia pri elektronno-luchevoi svarke [About instability of depth of pro-melting at electron beam welding]. Svarka elektronnym luchom. Moscow: MDNTP, 1974, pp. 26-30.

17. Cho J.H., Na S.J. Three-dimensional analysis of molten pool in GMA-laser hybrid welding. Welding Journal, 2009, no. 2, pp. 35-43.

18. Rai R., Palmer T.A., Elmer W.J., Debroy T. Heat transfer and fluid flow during electron beam welding of 304L strainless steel alloy. Welding Journal, 2009, no. 3, pp. 54-61.

19. Erokhin A.A. K raschetu glubiny proplavleniia pri elektronno-luchevoi svarke [To calculation of depth of pro-melting at electron beam welding]. Svarochnoe proizvod-stvo, 1976, no. 4, pp. 4-6.

Получено 05.06.2017

Об авторах

Ольшанская Татьяна Васильевна (Пермь, Россия) - кандидат технических наук, доцент кафедры сварочного производства, метрологии и технологии материалов Пермского национального исследовательского политехнического университета; е-mail: [email protected].

Федосеева Елена Михайловна (Пермь, Россия) - кандидат технических наук, доцент кафедры сварочного производства, метрологии и технологии материалов

Пермского национального исследовательского политехнического университета; е-mail: [email protected].

Колева Елена Георгиевна (София, Болгария) - доктор наук, Химико-технологический и металлургический университет; e-mail: [email protected].

About the authors

Tatyana V. Olshanskaya (Perm, Russian Federation) - Ph. D. in Technical Sciences, Associate Professor, Department of Welding Production, Metrology and Technology of Materials, Perm National Research Polytechnic University; e-mail: [email protected].

Elena M. Fedoseeva (Perm, Russian Federation) - Ph. D. in Technical Sciences, Associate Professor, Department of Welding Production, Metrology and Technology of Materials, Perm National Research Polytechnic University; e-mail: fedoseeva.pstu@ mail.ru.

Elena G. Koleva (Sofia, Bulgaria) - Doctor of Sciences, University of Chemical Technology and Metallurgy; e-mail: [email protected].