Full Text

Введение В настоящее время процесс аддитивного формирования изделий методом послойной проволочной наплавки успешно реализован различными научными группами, однако остается нерешенным целый ряд фундаментальных проблем, связанных со структурным состоянием сплавов, формируемых в данном процессе: грубое дендритное строение и текстура изделий из титановых сплавов [1], значительная пористость изделий из сплавов на основе алюминия [2], а также грубое дендритное строение и ликвация в изделиях из сплавов на основе железа [3]. Эти структурные особенности материалов, получаемых методом многослойной проволочной наплавки, значительно снижают уровень и однородность по направлению механических свойств (в частности, предела текучести, предела прочности, ударной вязкости, ресурса по малоцикловой и многоцикловой усталости). Ряд исследователей предлагают решать эту проблему путем применения в процессе аддитивного формирования заготовок послойной наплавкой вибрационных воздействий на расплавленный металл с целью управления формированием структуры [4; 5]. На сегодняшний день известны различные способы введения вибраций в зону сварки, такие как, например, электродуговые вибрации [6; 7], тепловые эффекты когерентного лазерного луча [8], прямое введение вибраций в материал через волновод [9]. Данные методы позволяют улучшить свариваемость без необходимости термообработки, что приводит к сокращению технологического цикла, снижению себестоимости и времени на производство изделий. В работах [4; 5; 10-12] однозначно установлено, что использование механических колебаний в процессах дуговой сварки приводит к улучшению структуры и свойств сварных соединений, в то время как при сварке без механических колебаний сварные швы обладают более грубой структурой. Погружение волновода в расплавленный алюминий позволяет осуществлять целенаправленную кавитацию, что способствует получению «пеноматериалов» и приводит к перемешиванию при химическом легировании расплава [13]. Авторы работы [9] добились снижения размера зерна и повышения твердости в сварном шве и в зоне термического влияния на алюминиевых сплавах и стали, хотя применяемые подходы к управлению амплитудой колебаний и энергией неясны и сомнительны. В работе [14] источником колебаний выступал модулированный ток, переменное электромагнитное поле или импульсный режим генерации лазерного излучения. Ультразвуковое воздействие успешно применялось при наплавке образцов из никелевого сплава Inconel 625 с чередующейся столбчатой/равноосной структурой по высоте путем последовательного включения и выключения ультразвукового воздействия [15]. Также в работе [16] показано положительное влияние введения ультразвуковых воздействий на формирование микро- и макроструктуры наплавленного материала из жаропрочного сплава на никелевой основе. Можно заключить, что вышеупомянутые работы имеют в основном экспериментальный характер и в большей степени направлены на исследование структурных изменений. Вместе с тем вибрационные воздействие влияют на гидродинамические процессы [17; 18], в том числе на процессы формирования капли, формирование ванны расплава и кристаллизацию наплавляемого валика, что позволяет управлять тепломассопереносом в жидкости, процессом кристаллизации и формой валика [19]. Вибрации влияют на характер движения в жидкой капле, что изменяет ее форму, подобного эффекта можно добиться изменением поверхностных сил, таких как сила натяжения и сила Марангони. Численное моделирование может быть использовано как мощный инструмент исследования мезоуровневых процессов, происходящих в сварочной ванне, в частности, вибрационного воздействия на течение металла, формирование микро- и макроструктуры [19-21]. Так, в работе [19] использовался метод объема жидкости (VOF) для исследования влияния ультразвука на поле давления, температурное поле и поток жидкости в ванне расплава. Было показано, что акустический поток способствовал увеличению скорости потока и высокотемпературной диффузии, что привело к уменьшению градиента температуры. Кроме того, на поверхности расплавленного металла появилось множество капиллярных волн, указывающих на то, что состояние поверхности ванны расплава также явно изменилось под воздействием ультразвуковой вибрации. С улучшением смачиваемости и растекания расплава под воздействием ультразвука ширина слоя осаждения увеличилась, а высота и глубина уменьшились, но площадь поперечного сечения слоя осаждения осталась прежней. Ранее авторами работы были проведены численные исследования процесса проволочной наплавки с применением метода конечных объемов с учетом силы поверхностного натяжения, силы Марангони, давления паров, и др. [22]. Были сформулированы условия, при которых наблюдается непрерывный перенос металла. Однако численное моделирование указанного процесса сеточными методами занимает слишком много времени. Использование метода гидродинамики сглаженных частиц при численном моделировании процессов тепломассопереноса позволяет в ряде случаев значительно увеличить производительность численной реализации [23-26]. Так, метод SPH (метод сглаженных частиц) обладает большим потенциалом по производительности [27]. Авторами статьи ранее была проведена работа по моделированию процесса наплавки проволоки с применением метода SPH [28-30]. Модель позволяет рассчитывать нестационарные объемные распределения температур, скоростей течения расплава, давлений, форму и размеры расплавленной ванны, форму свободной поверхности расплавленного металла, форму и размеры наплавляемого валика. В процессе проволочной наплавки осуществляется оплавление проволочного материала тепловым источником, и дальнейшее стекание капли на подложку. Таким образом, простейшей системой, подлежащей моделированию при исследовании влияния вибраций на формирование капли при наплавке, является объем жидкой фазы, удерживаемый силами поверхностного натяжения на некотором цилиндрическом объекте. В дальнейшем будет исследован процесс с наклонной проволокой с добавлением уравнений, описывающих тепломассоперенос, влияние вибраций на каплю, лежащую на подложке, после чего станет возможным сделать полное описание наплавки непрерывно подаваемой проволоки при вибрационных воздействиях. В предлагаемой работе рассмотрено влияние вибрационных воздействий на характер движения жидкости в капле, которое отражается в изменении значения коэффициента поверхностного натяжения. Предложена математическая модель течения жидкости с учетом силы поверхностного натяжения в формализме метода гидродинамики сглаженных частиц. Показано, что метод позволяет напрямую учитывать влияние вибрационных воздействий через введение дополнительных граничных условий. Проведена верификация разработанной математической модели сравнением с известными результатами натурных экспериментов для капли воды. Проведена серия численных экспериментов по определению влияния вибрационных воздействий на значение коэффициента поверхностного натяжения для стали марки 12Х18Н10Т. 1. Математическая модель 1.1. Постановка задачи В процессе отрыва капли от цилиндрического подвеса (рис. 1), на который действует источник вибрационных воздействий, время отрыва капли зависит от диаметра цилиндрического подвеса (d), массы капли (m), амплитуды и частоты колебаний (A, ω) и геометрии капли, которая, в свою очередь, зависит от теплофизических характеристик материала, поверхностного натяжения и вязкости, определяющих характер течения жидкости в капле. Рис. 1. Принципиальная схема процесса отрыва капли от цилиндрического подвеса Fig. 1. Schematic diagram of the droplet detachment process from a cylindrical suspension Длина звуковых волн, соответствующая исследуемым колебаниям λ = с/ν (с - скорость звука, ν - частота колебаний), составляет λ ~ 0,1 м и более, что много больше размеров области жидкого металла (l0 ~ 0,01 м). Соотношение λ >> l0 выполняется до частот порядка 20 кГц и для всех размеров характерных областей жидкого металла при его наплавке (капля, свисающая с проволоки; ванна расплавленного металла на подложке). В связи с выполнением условия λ >> l0 можно пренебречь акустическими эффектами при описании влияния вибрационных воздействий на аддитивное формирование металла наплавкой проволоки [17]. При описании движения жидкости положим, что расплав металла является несжимаемой ньютоновской жидкостью, при этом не учитывается влияние температуры. Акустические эффекты при вибрационных воздействиях также не учитываются. Варьируемыми параметрами являются: диаметр цилиндрического подвеса, масса капли, частота вибрационных воздействий, амплитуда вибрационных воздействий. Введение вибрационных воздействий в жидкую фазу приводит к генерации или перестроению существующих течений. Высокочастотные вибрационные эффекты могут описываться осредненными уравнениями движения, в которых появляется дополнительная вибрационная сила. Влияние вибраций конечной частоты вводится в уравнения движения посредством модуляции силы тяжести либо через граничные условия, т.е. задается периодический закон для скорости движения какой-либо границы расчетной области. Посредством вибрационного воздействия можно эффективно управлять состоянием системы: распределением полей скоростей и температур, что приводит к локальным изменениям физических параметров системы, таких как поверхностное натяжение или вязкость. Влияние вибраций на физические параметры системы возможно, если имеет место зависимость этих параметров от скорости движения или температуры [17; 18; 31-33]. В работе рассматривается взаимодействие твердого и жидкого металла. Для этого выделяются две фазы: Ωl - жидкая и Ωs - твердая, «цилиндрический подвес», объединение которых представляет собой всю исследуемую область - Ω (см. рис. 1). Движение металлического расплава можно описать как движение вязкой несжимаемой жидкости. В общем случае система уравнений будет состоять из дифференциальных уравнений, описывающих эволюцию плотности ρ, скоростей u в форме балансовых законов (уравнения баланса массы и импульса): (1) где u - скорость, ρ - плотность, fν - вязкие силы, fs - сила поверхностного натяжения, g - ускорение свободного падения. Поскольку в центре внимания данной работы находится гидродинамика жидкой капли, делается предположение о жесткой и неподвижной твердой фазе. Плотность ρ и давление P связаны с помощью уравнения состояния [34] (2) где и - скорость звука и плотность при нулевом приложенном напряжении. Для несжимаемых жидкостей вязкие силы примут следующий вид: (3) где η - динамическая вязкость. Следуя континуальному подходу Брэкбилла и Коте [35], основанному на континуальной поверхностной силе (CSF), эффекты поверхностного натяжения рассматриваются как объемные силы в уравнении (1), распределенные по межфазному объему конечной ширины (по поверхностным частицам). Сила поверхностного натяжения складывается из нормальной и тангенциальной составляющих: (4) где α - коэффициент поверхностного натяжения, - кривизна поверхности, Систему дифференциальных уравнений в частных производных (1) необходимо дополнить начальными условиями: (5) 1.2. Метод гидродинамики сглаженных частиц Численное решение предложенной математической модели было реализовано с применением метода гидродинамики сглаженных частиц (SPH). Суть метода гидродинамики сглаженных частиц заключается в аппроксимации полевых переменных {ρ, u, T} в системе (1). Переход от континуального представления среды к дискретному предполагает, что вместо непрерывной функции f (r), характеризующей какую-либо полевую переменную, вводится ее дискретный аналог fi. Кусочно-постоянная величина fi определяется для каждой частицы i как сумма N величин fj из частиц окружения j в пределах заданной дистанции h: (6) где mj - масса частицы, ρj - плотность частицы, W - сглаживающая функция (ядро). Здесь и далее используются эквивалентные обозначения , где и . Формально суммирование в выражении (1) производится только по частицам внутри некоторой области, для которой , где h - параметр дискретизации метода SPH (радиус сглаживания). Преимуществом метода SPH в таком представлении является то, что производные по координатам не отражаются на массах частиц mj и значениях полевых величин fj, а дифференцирование происходит только по весовым функциям Wij. Функция ядра зависит от расстояния между частицами и радиусом сглаживания h. При этом градиент некоторой величины можно представить в следующем виде: (7) Функция-ядро W представляет собой аппроксимацию δ-функции Дирака. Далее рассмотрим свойства функции W. На ядро накладывается ограничение компактности (нормировки): (8) Стоит отметить, что при h→0 ядро превращается в δ-функцию: (9) Свойство радиальной симметрии функции W позволяет перейти от градиента к ее производной в следующем виде: (10) Из свойства (10) вытекает антисимметричность градиента Легко также заметить, что выражение для плотности частиц получается усреднением масс по весовым функциям (11) Следует отметить, что выражение (11) подразумевает некоторую среднюю массу в области радиуса Для остальных полевых величин осреднение будет взвешенным. Применяя дискретизацию метода сглаженных частиц к системе (1), получим выражения для i-й частицы: (12) где - сила давления, - вязкая сила, - сила поверхностного натяжения. Действие на частицу i является результатом суммирования всех вкладов межчастичного взаимодействия с соседними частицами j. Вязкие силы и силы давления в уравнении импульса (12) дискретизируются в соответствии с формулировкой, предложенной Адами и др. [36; 37]: (13) Для дискретизации поверхностных сил сначала необходимо определить поверхность. В формализме SPH для этого вводят функцию уровня с, которая тождественно равна единице для каждой частицы жидкой фазы. Тогда, вычисляя градиент этой функции, можем определить нормаль к поверхности: (14) Силу поверхностного натяжения можно разделить на два слагаемых: Первое слагаемое представляет собой пропорциональные кривизне нормальные силы поверхностного натяжения, а второе представляет собой тангенциальные силы Марангони, обусловленные градиентом поверхностного натяжения. Первое слагаемое имеет вид: (15) Второе слагаемое примет следующий вид: (16) где I - единичный тензор, - производная коэффициента поверхностного натяжения по температуре. Однако в работе не учитываются температурные эффекты, поэтому выражение (16) принимается равным 0. При интегрировании полученных дифференциальных уравнений для вязких течений временной шаг ограничен вязкой составляющей, в отличие от известного условия Куранта, как, например, в сеточных методах. Для стабильного численного моделирования шаг по времени должен удовлетворять следующим условиям [38]: 1) условие Куранта: 2) условие, возникающие при учете диффузии: 3) условие, возникающее при учете поверхностного натяжения: 4) условие, возникающее при учете термокапиллярной диффузии: Для численного решения системы уравнений (1) использование классических явных схем интегрирования (например, схемы Верле) в методе SPH представляется невозможным, поскольку силы вычисляются по скоростям, определенным в полушаге по времени. Эта ситуация может быть улучшена путем вычисления экстраполированных на полный шаг скоростей частиц: (17) Добавляя интегрирование уравнений для локальной плотности и внутренней энергии, запишем полную схему Верле с учетом экстраполяции следующим образом: 1. b. 2. Определение правых частей , в конце шага по времени: 3. Разделение эволюции во времени для плотности и энергии на два отдельных шага аналогично интегрированию скорости u, дает существенные преимущества по отношению к простой схеме Эйлера для ρ и T, которая, в свою очередь, как известно, приводит к нарушению балансовых уравнений. 2. Численные эксперименты и результаты В первой серии численных экспериментов была рассмотрена задача моделирования процесса отрыва капли воды от вертикального цилиндра (проволоки) без учета теплопроводности. Вибрации в систему вносились через проволоку. В работе [33] рассмотрено экспериментальное исследование поведения капли воды в условиях высокочастотных механических вибраций и влияние таких воздействий на параметры, характеризующие эффекты поверхностного натяжения. Результаты данной статьи были использованы для верификации разработанной модели (рис 2). Рис. 2. Зависимость коэффициента поверхностного натяжения от амплитуды скорости при различных частотах [33] Fig. 2. Surface tension coefficient dependence on the peak velocity at different frequencies [33] Для определения модельного параметра коэффициента поверхностного натяжения α использовались сталагмометрический (метод отрыва капли) и метод висячей капли. Сталагмометрический метод основан на предположении, что в момент отрыва капли сила тяжести, действующая на нее, равна силе поверхностного натяжения, приложенного к наиболее узкой части капли (рис. 3): (18) где σ - коэффициент поверхностного натяжения, m - масса капли, g - проекция ускорения свободного падения, d - диаметр наиболее узкой части капли. Рис. 3. Схематичная иллюстрация сталагмометрического метода Fig. 3. Schematic illustration of the stalagmometric method Рис. 4. Схематичная иллюстрация метода висячей капли: O - точка на профиле капли; Φ - угол между касательной в точке O и горизонтальной осью; s - соответствующая длина дуги от вершины капли до точки O; x - расстояние точки O от вертикальной оси; z - высота точки O от горизонтальной оси; R1 - радиус кривизны в вершине капли; R2 - радиус кривизны в точке O Figure 4. Schematic illustration of the pendant drop method: O, point on the profile of the drop; Φ, the angle between the tangent at point O and the horizontal axis; s, the respective arc length from the apex of the drop to the point O; x, the distance of the point O from the vertical axis; z, the height of the point O from the horizontal axis; R1, the radius of curvature at the drop apex; R2, the radius of curvature at the point O Для сталагмометрического метода была предложена численная процедура по определению параметра α. По заданной геометрии (см. рис. 3) был определен диаметр проволоки d и масса всей капли m. По известному коэффициенту поверхностного натяжения σ вычислялось, при каком ускорении свободного падения капля должна оторваться от проволоки: После этого последовательно запускался расчет при разных α, и с использованием метода деления отрезка пополам был определен параметр α, соответствующий физическому коэффициенту поверхностного натяжения σ. Очевидный недостаток данного метода заключается в том, что с его помощью невозможно определить коэффициент поверхностного натяжения по заданному параметру α, а лишь узнать, больше он или меньше необходимого значения. Метод висячей капли (pendant drop method) позволяет определить коэффициент поверхностного натяжения по форме капли. Силы, определяющие форму висящей капли, - поверхностное натяжение и гравитация. С одной стороны, поверхностное натяжение стремится минимизировать площадь поверхности и придать капле сферическую форму. С другой стороны, гравитация вытягивает каплю из этой сферической формы, и в результате получается типичная грушевидная форма (рис. 4). В общем случае для свободной поверхности и при наличии приложенного «избыточного давления» ΔP на границе раздела в равновесии существует баланс между гидростатическим давлением и эффектами поверхностного натяжения . Уравнение Юнга - Лапласа принимает следующий вид: (19) Исходя из геометрических построений, изображенных на рис. 4, получаем: 1/R1 = dФ/ds, 1/R2 = sin (Ф)/x, dx/ds = cos (Ф), dz/ds = sin (Ф). При использовании параметра формы капли выражение (19) примет следующий вид: (20) При известном параметре β, можно определить коэффициент поверхностного натяжения: (21) где - радиус вписанной в каплю окружности. Для определения параметра b необходимо найти максимальный диаметр капли и диаметр капли на высоте . После чего параметр b вычисляется по следующему соотношению [39]: а b c Рис. 5. Визуализация отрыва капли при численном описании сталагмометрического метода в условиях вибрационных воздействий: а - ωА = 0,05; b - ωА = 0,15; c - ωА = 0,25 Fig. 5. Visualization of drop detachment at numerical description of stalagmometric method under vibration influences: а - ωА = 0.05; b - ωА = 0.15; c - ωА = 0.25 (22) Метод висячей капли позволяет определить текущее значение коэффициента поверхностного натяжения, таким образом, идентификация параметра a происходит за одну итерацию. Верификация модели проводилась в расчетах для капли воды. Варьируемыми параметрами являлись масса капли и амплитуда колебаний, частота колебаний равнялась 6 кГц. Определяющим параметром оказалось значение амплитуды скорости ωА. Величина коэффициента поверхностного натяжения монотонно снижается с ростом амплитуды скорости. На рис. 4 представлена визуализация капли воды на подвесе вибрирующем в вертикальном направлении в одинаковые моменты времени. Повышение амплитуды скорости колебаний приводит к изменению характера движения жидкости в капле и к отрыву капли от подвеса, что равносильно снижению коэффициента поверхностного натяжения. При ωА = 0,05 (рис. 5, а) капля устойчиво удерживалась силами поверхностного натяжения, в то время как при при ωА = 0,15 (рис. 5, b) начинался отрыв, а при ωА = 0,25 (рис. 5, с) капля уже оторвалась. Результаты расчета (рис. 6) показали хорошее качественное и количественное (расхождение не превышает 10 %) совпадение с экспериментом (см. рис. 2). Верифицированная модель использовалась для исследования поведения поверхностного натяжения расплава металла при вибрационных воздействиях (рис. 7). В этой серии численных экспериментов также не учитывалась теплопроводность и связанная с ней сила Марангони. В качестве исследуемого материала была выбрана сталь марки 12Х18Н10Т с характеристиками, представленными в табл. 1. Рис. 6. Результаты расчета зависимости коэффициента поверхностного натяжения воды от амплитуды скорости в сравнении с экспериментальными данными при частоте колебаний 6 кГц Fig. 6. Calculation results of the water tension coefficient dependence on the peak velocity in comparison with the experimental data (ω = 6 kHz) Варьируемыми параметрами являлись масса капли и циклическая частота колебаний, амплитуда колебаний. Амплитуда скорости варьировалась от 0,0 м/с до 2,0 м/с с шагом 0,2 м/с. В табл. 2 приведены параметры расчета. Результаты проведенных численных экспериментов приведены на рис. 8. Эффект снижения коэффициента поверхностного натяжения вследствие увеличения скорости колебаний подводимых вибраций хорошо согласуется с известными экспериментальными данными. При пиковой скорости колебаний величиной 2,0 м/с наблюдается снижение коэффициента поверхностного натяжения на 30 %. Снижение коэффициента поверхностного натяжения должно облегчить реализацию непрерывного стекания металла с проволоки [22], что может положительно сказаться на формировании металла при проволочной наплавке, а также позволит дополнительно управлять формой наплавляемых валиков. Рис. 7. Визуализация численного моделирования метода отрыва капли для стали марки 12Х18Н10Т при ωА = 0,5 Fig. 7. Visualization of numerical simulation of the drop detachment method for 12X18H10T steel grade with ωА = 0.05 Таблица 1/Table 1 Принятые при расчете теплофизические характеристики наплавляемого материала (сталь 12Х18Н10Т) Thermal and physical characteristics of the weld material used in the simulation (steel 12Х18Н10Т) Характеристика Обозначение Размерность Значение Температура ликвидус TL °С 1515 Температура солидус TS °С 1459 Удельная теплоемкость C Дж·кг-1·K-1 710 Плотность ρ кг·м-3 7680 Теплопроводность λ Вт·м-1·K-1 26 Скрытая теплота плавления Hf Дж/кг 276 000 Динамическая вязкость μ Па·с 0,007 Коэффициент поверхностного натяжения s Н·м-1 1,615 Температурный коэффициент поверхностного натяжения Н·м-1·K-1 -0,00043 Таблица 2/Table 2 Параметры численного моделирования Numerical simulation parameters Параметр Обозначение Размерность Значение Параметр дискретизации h мм 0,1 Количество частиц n шт. 57709…81815 Шаг по времени dt с 2,5e-06 Количество шагов N шт. 100 000 Физическое время t с 0,25 Рис. 8. График зависимости коэффициента поверхностного натяжения стали 12Х18Н10Т от пиковой скорости вибрационных воздействий Fig. 8. The graph of dependence of surface tension coefficient of steel 12Х18Н10Т on the peak velocity of vibration influences Заключение В ходе работы была реализована математическая модель образования капли металлического расплава из проволоки в условиях вибрационного воздействия в формализме метода гидродинамики сглаженных частиц. Для верификации разработанной модели была проведена серия численных экспериментов по моделированию процесса отрыва капли воды от цилиндрического подвеса, который подвергался вибрационным воздействиям с разными значениями амплитуды скорости. Результаты численных экспериментов качественно и количественно совпали с результатами известных натурных экспериментов. Верифицированная математическая модель позволяет исследовать влияние вибрационных воздействий на коэффициент поверхностного натяжения для большинства жидкостей, в том числе и для жидких металлов. Более того, математическое моделирование становится единственным доступным способом исследования такого явления для жидких материалов с большой температурой. С использованием разработанной математической модели была проведена серия численных экспериментов по моделированию процесса отрыва капли расплавленного металла (сталь марки 12Х18Н10Т) от цилиндрического подвеса из того же материала. Выявлено, что при повышении амплитуды скорости вибрационных воздействий коэффициент поверхностного натяжения стали марки 12Х18Н10Т тоже снижается. При амплитуде скорости колебаний величиной 2,0 м/с наблюдается снижение коэффициента поверхностного натяжения на 30 %. Таким образом предложенная математическая модель позволяет явно моделировать эффект влияния вибрационных воздействий на поверхностное натяжение различных жидкостей и может в дальнейшем быть использована при описании процесса проволочной наплавки с применением вибрационных воздействий.

About the authors

R. P Davlyatshin

Perm National Research Polytechnic University

A. V Perminov

Perm National Research Polytechnic University

Y. V Bayandin

Perm National Research Polytechnic University; Institute of Continuous Media Mechanics UB RAS

F. R Saucedo-Zendejo

Autonomous University of Coahuila

D. N Trushnikov

Perm National Research Polytechnic University

References

Statistics

Views

Abstract - 577

PDF (Russian) - 260