Механизм возникновения колебаний глубины канала при лучевых способах сварки

Размещено на http://www.allbest.ru/

Тульский государственный университет

Механизм возникновения колебаний глубины канала при лучевых способах сварки

Васильев А. А, Ерофеев В. А, Судник В. А.

Аннотация

канал лучевой сварка температура

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

Введение

При лучевых способах сварки характерным дефектом являются корневые пики глубины проплавления и пустоты, рис. 1. Возникновение этих пустот объясняют повторяющимися процессами, приводящими к периодическому схлопыванию горловинной части парогазового канала (ПГК) [1]. На рис. 1 показан шлиф участка окончания шва с линейным уменьшением тока луча при электронно-лучевой сварке (ЭЛС) Контракт между MTU Aero Engines, Мюнхен, Германия и Тульским государственным университетом, 1999-2000 титанового б-сплава Ti 6-2-4-2, что показывает независимость появления корневых дефектов от вводимой мощности.

Рис. 1 Продольное сечение шва титанового сплава Ti 6-2-4-2

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

Каплан, Мацунава и др. [4] предложили аналитическую модель схлопывания или коллапса ПГК при импульсной лазерной сварке, расчеты по которой сравнили с результатами просвечивающей рентгеновской микроскопии при облучении жидкого цинка. Хорошее соответствие теоретического получения ПГК и экспериментом позволила авторам предложить модель процесса, в которой учитывалась повторная конденсация при быстром охлаждении потоком защитного газа в ПГК в период паузы.

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

Туричин и др. [6] разработали аналитическую динамическую модель, основанную на вариационных принципах и формализме механики Лагранжа, лазерной и гибридной сварки с учетом схлопывания или спайкинга ПГК. Авторы свели модель к системе обыкновенных дифференциальных уравнений. Они показали, что нестабильное поведение ПГК при сварке волоконным лазером связано с величиной радиуса луча и распределением интенсивности.

Ерофеев [7] разработал упрощенную теоретическую модель схлопывания ПГК при ЭЛС, в которой глубина проплавления и формирование ПГК определялась по равенству значения температуры ПГК температуре кипения легкоиспаряющегося компонента сплава. Решение уравнений модели имело колебательный характер относительно глубины ПГК, но вследствие высокой частоты колебаний и инерционности тепловых процессов температурное поле за пределами сварочной ванны за период колебаний не изменялось. Соответствие расчетных профилей поперечного сечения шва макрошлифам достигалось введением калибровочных коэффициентов для учета влияния перемешивания расплава и погрешностей данных по температуре кипения (интенсивного парообразования). Значения коэффициентов определены для конкретного сплава при сравнении расчетных и опытных значений ширины шва и глубины проплавления.

Экспериментальные наблюдения частот образования корневых пиков показывают, что эти они лежат в диапазоне частот от единиц Гц [8] до нескольких килогерц [9].

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

Физические явления

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

Луч падает на металл и заглубляется, рис. 2, а. Энергия луча вызывает интенсивное испарение на дне ПГК, и образуется поток пара V рис. 2, б. Реактивная отдача пара создает давление на жидкий металл. Луч постепенно заглубляется в образовавшийся ПГК, при этом поверхность, на которую падает луч, заглубляется в толщину металла, рис. 2, б. Часть тепла уходит в металл, часть расходуется на испарение. При этом боковые стенки ПГК удерживаются давлением потока пара со его дна P_v, рис. 2, в. Стенки ПГК охлаждаются теплоотводом вглубь металла. Поток пара препятствует охлаждению, так как пар конденсируется на поверхности стенок ПГК. Пока тепловой энергии потока пара достаточно для поддержания температуры стенок на уровне температуры, при которой парциальное давление пара вместе с реактивным давлением отдачи паров достаточно для противодействия силам, стягивающим ПГК, последний устойчиво существует и заглубляется внутрь металла. При этом количество пара, выходящего из ПГК, уменьшается по мере заглубления, так как все большая часть этого пара конденсируется на боковой поверхности ПГК.

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

Рис. 2 Стадии формирования ПГК при лучевых способах сварки (a-e) согласно Ерофееву [7]

Процессы испарения характеризуются распределением мощности луча в ПГК P_l, рис.3, а, которая выделяется на его дне. Часть P_vap(z) этой мощности переносится паром и поглощается стенками ПГК Непоглощенная мощность пара P_ext выносится из ПГК.

Если мощность всего потока пара будет поглощена стенками ПГК, рис. 3, б и исчезнет давление пара в горловине ПГК, силы поверхностного натяжения начнут стягивать горловину. Мощность луча при схлопывании быстро перераспределяется со дна ПГК (P_l1) на его горловину (P_l2). Несмотря на это перераспределение ПГК захлопнется вследствие инерционности движения расплава.

Рис. 3 Распределение мощности в ПГК: а) устойчивый и б) неустойчивый ПГК

Математическая модель

Ввиду малого радиуса ПГК при его значительной длине допустимо считать цилиндром радиуса r_c, конец которого является полусферой, а горловина имеет тороидальную форму с радиусами, равными ширине сварочной ванны Z_R и разности Z_R-r_c,

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

Уравнение энергии

Во всех точках пространства моделирования состояние вещества определяется энтальпией, вычисляемой решением нестационарного уравнения сохранения энергии:

где Н - объемная энтальпия, t - время, T - температура, v_w - скорость сварки, - коэффициент теплопроводности, зависящий от температуры, рис. 4; - функция, описывающая объемное распределение мощности луча и пара по объёму ПГК вдоль его длины z.

Связь между энтальпией и температурой определяется соотношением:

где - доля жидкой фазы, - доля паровой фазы С- объёмная теплоёмкость, H_L,H_V -теплоты плавления и парообразования

Начальные условия:

Граничные условия:

- на верхней и нижней поверхностях учтена радиационная теплоотдача:

где - степень черноты поверхности, - постоянная Стефана-Больцмана и S - тощина листа.

- в бесконечности

(5)

Источник Q(z) описывает поглощение мощности луча на дне ПГК

расход энергии луча на образование пара

и распределение мощности потока пара P_V(z) вдоль ПГК для поддержания равновесия поверхности ПГК

где H_k - энтальпия металла, при которой давление пара обеспечивает равновесие ПГК, H_c - фактическое значение энтальпии на глубине z

Мощность теплоты пара передаётся стенкам ПГК при конденсации пара металла, истекающего со дна ПГК. По мере движения от дна к горловине ПГК, мощность потока пара убывает.

Поэтому из ПГК вытекает поток пара, тепловая мощность которого существенно меньше мощности потока пара на его дне.

Таким образом, источник теплоты описывается зависимостью:

Уравнение движения

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

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

где - поверхностное натяжение для пары расплав - пар при температуре поверхности ПГК T_c, - плотность расплава, z_R- глубина расположения горловины ПГК

В момент исчезновения потока пара ,

Длительность t_c схлопывания ПГК определяется по условию

Численное решение

При компьютерной имитации пространство искусственно ограничено и имеет форму параллелепипеда. Размеры области имитации и положение центра координат в ней находятся на основе предварительного расчёта размеров высокотемпературной зоны, в которой температура превышает значение, равное 0,2 температуры плавления. Уравнение теплопроводности решали методом конечных разностей на трехмерной равномерной сетке с шагом, равным 0,25 диаметра луча. Шаг времени был выбран по физическому условию недопустимости агрегатного изменения в объеме металла за один шаг сетки.

На начальном этапе принимается, что металл холодный, ПГК отсутствует, и мощность луча выделяется в поверхностном контрольном объеме. В цикле времени рассчитываются энтальпии температуры в узлах сетки. Затем определяются узлы ПГК, в которых температура больше температуры кипения. Определяется мощность потока пара со дна ПГК, которая затем перераспределяется в контрольные объемы по длине ПГК, начиная со дна. Если этой мощности хватает для всех узлов, лежащих на пути луча, то цикл моделирования продолжается, в ходе которого ПГК углубляется. Если нет, то моделируется схлопывание ПГК, т.е. глубина ПГК принимается равной нулю и начинается моделирование следующего цикла заглубления луча.

Свойства сплава

Имитацию проводили для титанового б-сплава Ti 6-2-4-2, T_s = 1230 °С и T_L = 1575 °С, теплофизические свойства которого принимали по данным работы [10]. Температура и энтальпия связаны между собой кусочно-линейной функцией T(H), рис. 4.

Рис. 4 Зависимость энтальпии и теплопроводности от температуры титанового сплава Ti 6-2-4-2

Результаты решения

Воспроизводили процесс формирования ПГК и сварочной ванны при ЭЛС листов толщиной 1,2 см, ускоряющем напряжении 140 кВ и токе 14,5 мА, диаметре луча 0, 75 мм и скорости сварки 1 см/с, рис. 6. Максимальная глубина ПГК достигала 10,3 мм. Температура стенок ПГК составила 2815 °С, а на дне ПГК - 2904 °С. Избыточное давление пара на дне ПГК равно 0,3 Н/см², или на ~ 3 % выше окружающего давления.

Рис. 6 Распределение температуры в продольном сечении (a), на поверхности листа (b) и в плоскости воздействия луча (c)

Численное решение показало, что процесс имеет колебательный характер, при котором глубина ПГК непрерывно изменяется, рис. 7(a). Процесс через некоторое время переходит в квазистационарное состояние, при котором амплитуда колебаний глубины ПГК стабилизируется. В периоде колебания мощность парообразования на дне ПГК изменяется незначительно. При этом мощность потока пара, истекающего из ПГК, изменяется от максимального, равного количеству образующегося пара, до нуля в момент схлопывания, рис.7 (b). Глубина ПГК периодически изменяется, рис. 8.

Рис. 7 Изменение глубины ПГК z_c и сварочной ванны z_l, мощности потока пара на дне ПГК P_vap и на выходе из ПГК P_ext

Рис. 8 Продольные сечения сварочной ванны в моменты времени t₁ (a),. t₂(b), t₃(c) и строение шва (d) (1 - исходный металл, 2 - расплав, 3 - металл шва однократного плавления, 4 - металл шва периодического плавления

Моделирование схлопывания ПГК (11) показало, рис. 9, что радиус его горловины в начале схлопывания меняется медленно, но к концу сильно ускоряется.

Рис. 9 Изменение радиуса горловины ПГК после исчезновения потока пара при сварке титанового сплава Ti 6-2-4-2

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

Колебания глубины ПГК периодически изменяет глубину сварочной ванны, и металл корня шва периодически плавится и кристаллизуется, рис. 8. Возникновение пустот наиболее вероятно в этой области. Это позволяет оценивать вероятность появления пустот по толщине этой зоны периодического плавления.

Адекватность моделирования

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

Расчетная частота по данным рис. 8 составляет 15 Гц.

Эти значения соответствуют нижним значениям, приводимых в литературе [8].

Заключение

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

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

2. Разработана физико-математическая модель формирования ПГК при лучевых способах сварки, учитывающая процессы испарения и конденсации пара металла в ПГК.

3. Выполнено численное моделирование процессов испарения и конденсации пара металла в ПГК, которое подтвердило возникновение периодических колебаний глубины ПГК согласно предложенной гипотезе.

4. Показано, что предложенная модель дает понимание возможной причины возникновения нестабильности глубины ПГК и сварочной ванны.

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

1. Рыкалин Н. Н., Зуев И.В., Углов А.А. Основы электронно-лучевой обработки материалов. М.: Машиностроение. 1978. 239 стр. / Rykalin N.N., Zuev I.V., Uglov A.A. Basics of electron-beam treatment of materials (in Russian). Moscow. Mashinostroenie.1978. 240 p.

2. Kroos J., Gratzke U., Simon G. Towards a self-consistent model of the keyhole in penetration laser beam welding // Journal Physics D: Applied Physics, 1993, Vol. 26. P. 474-480.

3. Sudnik W., Radaj D., Erofeev W. Computerized simulation of laser beam welding, modeling and verification // Journal Physics D: Applied Physics, 1996, Vol. 29. P. 2811 - 2817.

4. Kaplan A. F. H., Mizutani M. Katayama S., Matsunawa A. Computerized Unbounded keyhole collapse and bubble formation during pulsed laser interaction with liquid zinc// Journal Physics D: Applied Physics, 2002, №35, P. 1218-1228.

5. Кривцун И. В., Сухоруков С. Б., Сидорец В. Н., Ковалев О. Б. Моделирование процессов испарения металла и газодинамики металлического пара в парогазовом канале при лазерной сварке // Автомат. сварка. - 2008.- №10. - С. 19 - 26 / Krivtsun I.V., Sukhorukov S.B., Sidorets V.N., Kovalev O. B. Evaporation processes modeling of metal and metal vapor in the gas dynamics of vapor-gas channel in laser welding. Automatic Welding. 2008, № 10. Pp. 19 - 26 (in Russian).

6. Туричин Г.А., Цибульский И.А., Валдайцева Е.А., Карасев М.В. Теория и технология гибридной сварки металлов больших толщин // Сб. трудов 6-й междунар. конф. "Лучевые технологии и применение лазеров". Изд-во СПбГПУ. - 2009. - С. 11 - 18./ Turichin G., Valdaitseva E., Pozdeeva E. et al. Theoretical investigation of dynamic behavior of molten pool in laser and hybrid welding with deep penetration, Paton welding journal, 2008 (7), p. 11-15.

7. Ерофеев В. А. Моделирование особенностей формирования парогазового канала при электронно-лучевой сварке // Сварка и диагностика, 2009, № 4, стр. 2-8./ Erofeev V. A. Modelling of formation features of a keyhole during electron beam welding. (In Russian).Welding and Diagnostics. 2009. №4. P. 12 - 18.

8. Cho M.H., Farson D.F., Kim J.I.// Control of chaos in laser-induced vapor capillaries// J. Laser Appl., 2003, Vol. 15, No. 3, P. 161-167.

9. Фольпп Й., Гатцен М., Фоллертсен Ф. Аналитическая модель динамики парогазового канала при лазерной сварке с глубоким проплавлением. Автоматическая сварка. 2013. №3. Стр. 15 - 18 / Volpp J., Gatzen M., Vollertsen F. Simplified analytical modeling of dynamic behavior of the keyhole for different spatial laser intensity distributions during laser deep penetration welding. Intern. Conf. Mathematical Modelling and Information Technologies in Welding and Related Processes. 2012, Katsiveli, Crimea, Ukraine. E.O. Paton Electric Welding Institute, Kiev. P. 27 - 30.

10. Sudnik V. A., Erofeev V. A., Richter K.-H., Heins K.-U. 2006. Numerical modelling of the EBW process. 16th Conf. Computer Technology in Welding and Manufacturing and 3rd Int. Conf. n Mathematical Modelling and Information Technologies in Welding and Related Processes. Kiev, Ukraine. PWI. P. 295-300.

Размещено на Allbest.ru