Комплексное математическое моделирование процессов плазменно-дугового проволочного напыления покрытий Текст научной статьи по специальности «Физика»

ВЕСТНИК ЮГОРСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА

2015 г. Выпуск 2 (37). С. 33-41

УДК 621.793.74

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

М. Ю. Харламов, И. В. Кривцун, В. Н. Коржик, А. И. Демьянов

Работа выполнена в рамках совместного проекта РФФИ № 14-08-90428 и НАН Украины № 06-08-14, между Институтом теоретической и прикладной механики СО РАН с Югорским государственным университетом и Институтом электросварки им. Е. О. Патона, г. Киев.

Введение

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

К одним из новых перспективных методов газотермического напыления относится плазменно-дуговое проволочное напыление [3]. В ИЭС им. Е. О. Патона НАН Украины выполняются исследования, направленные на развитие данного метода напыления. В частности, опубликована серия работ, посвященных теоретическими исследованиями и построением математических моделей процессов плазменно-дугового напыления [4-10]. Большой интерес представляет обобщение данного материала с целью получения комплексной математической модели процесса плазменно-дугового напыления. Такая модель, с соответствующим набором входных и выходных параметров будет полезна для анализа закономерностей формирования плазменно-дуговых покрытий с учетом особенностей влияния тех или иных технологических факторов и других параметров. Формулирование комплексной математической модели протекающих при плазменно-дуговом напылении процессов, а также рассмотрении особенностей математического моделирования его отдельных этапов и составляло цель данной работы.

Плазменно-дуговое проволочное напыление покрытий

Схема процесса плазменно-дугового распыления проволочных материалов, принятая при построении комплексной математической модели, представлена на рис. 1 . Дуга постоянного тока I горит между тугоплавким водоохлаждаемым катодом и токоведущей проволокой, находящейся за срезом сопла плазмотрона на расстоянии Z2 от начального сечения расчетной области (z = 0), расположенного вблизи рабочего конца катода (см. рис. 1), и при z > Z2 имеет место бестоковое инерционное движение плазмы. Распыляемая проволока круглого сечения

33

М. Ю. Харламов, И. В. Кривцун, В. Н. Коржик, А. И. Демьянов

радиусом Rw подается в плазменную дугу со скоростью vw. Формируемая плазмотроном электрическая дуга замыкается на конце проволоки, являющейся анодом дуги. Вся проволока нагревается протекающим через нее током дуги I. Подаваемый в плазмоформирующее сопло длиной Zi и радиусом Rc плазмообразующий газ с массовым расходом Gi нагревается электрической дугой и истекает во внешнюю газовую среду. Открытый участок течения дуговой плазмы может обдуваться коаксиальным потоком внешнего газа с массовым расходом G2. При этом обдувающий газ подается через кольцевой канал Ri < r < R2 под углом а к оси плазмотрона.

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

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

Рисунок 1 - Схема процесса плазменно-дугового проволочного напыления: 1 - катод; 2 - сопло;

3 - канал подачи обдувающего газа; 4 - проволока-анод; 5 - дуговая плазма; 6 - обдувающий газ;

7 - область смешения; 8 - внешняя газовая среда

Структура комплексной математической модели

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

В первом приближении комплексная модель может быть представлена семью взаимосвязанными моделями, описывающими ту или иную подсистему (этап) технологического процесса (рис. 2): 1 - модель плазменной струи, генерируемой плазмотроном с внешней проволокой-анодом; 2 - модель нагрева и плавления проволоки; 3 - модель формирования жидкой прослойки на рабочем конце проволоки; 4 - модель формирования и срыва капель расплавленного металла; 5 - модель нагрева и движения расплавленных частиц в плазменном потоке с учетом их дробления; 6 - модель динамического и теплового взаимодействия напыляемых частиц с основой и 7 - модель формирования покрытия. При описании схемы комплексной

34

Комплексное математическое моделирование процессов плазменно-дугового проволочного напыления покрытий

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

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

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

Модели нагрева и плавления проволоки и формирования жидкой прослойки на ее конце [7, 8] являются согласованными и решаются совместно. Решение определяется по балансу расплавленной части проволоки и толщины удерживаемой на конце проволоки прослойки жидкого металла, которое устанавливается при определенном положении конца проволоки относительно оси плазменной струи. Поэтому входные параметры I для данных моделей являются общими и включают в себя скорость подачи проволоки, ее положение относительно среза сопла плазмотрона, теплофизические параметры материала проволоки, а также распределенные и интегральные характеристики плазменного потока Oj, являющиеся выходными

данными первой модели. Выходными данными O23 будут выступать температурное поле в

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

Рисунок 2 - Схема комплексного математического моделирования процесса плазменно-дугового

проволочного напыления

При моделировании процесса диспергирования проволоки (формирования и срыва капель) [9] входными данными I будут выступать выходные данные моделей первой, второй и

35

М. Ю. Харламов, И. В. Кривцун, В. Н. Коржик, А. И. Демьянов

третей стадий процесса, в частности: пространственные распределения скорости и температуры плазмы, положение расплавленного конца проволоки в струе, толщина жидкой прослойки, удерживаемой на конце проволоки, и скорость течения в ней расплава. В качестве основных выходных данных O4 будут выступать диаметр срываемых капель, их частота срыва и начальная скорость движения.

Выходные параметры четвертой стадии процесса O4, наряду с выходными параметрами O23 и O1 являются входными параметрами для модели ускорения, нагрева и дробления расплавленных частиц в плазменной струе [5]. Выходные параметры O5 данной модели - скорость и траектория движения частиц и их фрагментов, тепловое состояние частиц и фрагментов, а также распределения напыляемых частиц по поверхности основы по скорости, температуре и размерам.

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

И, наконец, при моделировании формообразования покрытия входными параметрами I7

являются выходные параметры O5 и O6, а также набор входных параметров I . Выходные параметры O7 : форма и размеры единичного слоя покрытия, параметры микроструктуры единичного слоя покрытия, распределение напряжений и деформаций и др.

Рассмотрим особенности математического описания рассмотренных этапов (стадий процесса).

Математические модели процессов плазменно-дугового напыления

Модель плазменной струи, формируемой плазмотроном с внешней проволокой-анодом. При построении математической модели использовались следующие допущения [4-6]: рассматриваемая плазменная система обладает цилиндрической симметрией, а протекающие процессы предполагаются стационарными; обдувающий газ подается осесимметричным потоком через кольцевой канал, течение этого газа в канале полагается ламинарным; плазма находится в состоянии локального термодинамического равновесия, собственное излучение плазмы -объемное; основным механизмом нагрева плазмы является джоулево тепловыделение (работой сил давления и вязкой диссипацией можно пренебречь), а перенос энергии в столбе происходит за счет теплопроводности и конвекции (естественная конвекция в расчет не принимается); течение плазмы вязкое, дозвуковое, режим течения турбулентный; внешние магнитные поля отсутствуют; плазменная компонента смеси является инертной (Аг) и не вступает в химические реакции с внешним газом.

Основу математической модели формирования плазменной струи составляет система магнитогазодинамических (МГД) уравнений в приближении турбулентного пограничного слоя для осредненных по времени значений температуры и скорости плазмы, которая имеет вид [ 4]:

-(pu) + -—(rpv ) = О, (1)

dz r dr

f du _du

p\ u---+ v —

l dz dr

1 d f _du\ d f И2ф Л,

r dr l rnr J-& lp+^0—);

pC

p

f dT dT Л

u + v

v d z dr y

1 d

r dr

f

v

rX

dT

dr

Л

a

(2)

(3)

36

Комплексное математическое моделирование процессов плазменно-дугового проволочного напыления покрытий

Здесь T - осредненная температура плазмы; v = (pv + pv') / р, где v - осредненная радиальная скорость, р- осредненная плотность плазмы, р’ и v’ - пульсации плотности и радиальной скорости; и - осредненная аксиальная скорость плазмы; p - давление; Cp - удельная теплоемкость при постоянном давлении; а - удельная электропроводность плазмы; j - вектор плотности электрического тока; щ- объемная плотность мощности собственного излучения; П и х - полные коэффициенты динамической вязкости и теплопроводности плазмы; ро - универсальная магнитная постоянная; Hv - азимутальная составляющая магнитного поля тока

дуги, которая связана с компонентами плотности электрического тока следующими соотношениями:

Л

1 д

r dr

(rHd); j

dHL

dz

(4)

С целью более корректного описания электромагнитных характеристик дуги (не делая допущения о малости радиальной компоненты плотности электрического тока по сравнению с аксиальной) используется уравнение для напряженности магнитного поля тока дуги [6]

d 1 d 1 д \ 1 dHA

dr ra dr dz a dz

(5)

Для описания турбулентности была использована двухпараметрическая k-s модель, получившая широкое распространение в мировой практике. При этом, коэффициенты динамической вязкости и теплопроводности плазмы, используемые в приведенных выше уравнениях, имеют вид п = V+Vt; X = X + Xt, где пи х - коэффициенты молекулярной вязкости; п и Xt -коэффициенты турбулентной вязкости и теплопроводности.

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

( dm _ dm

р\ и—1 + v—1 I dz dr

1 d

r dr

ГРР12

dm1 | dr J ’

(6)

m1 + m 2 = 1.

Здесь щ (r, z) и m2 (r, z) - относительные массовые концентрации плазмообразующего и внешнего газа в плазменной смеси, щ. = pjр, где р. - приведенные плотности компонент (i

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

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

Модели нагрева и плавления проволоки-анода и формирования жидкой прослойки на ее конце. В условиях плазменно-дугового распыления тепловое состояние проволоки-анода будет определяться совокупностью следующих физических процессов [7]: конвективно-кондук-тивным теплообменом плазменного потока и окружающего газа с боковой поверхностью проволоки Qe'; обменом энергией теплового излучения между плазмой и поверхностью проволоки Qr; действием электрической дуги, вводящей тепло через анодное пятно на торце проволоки Qa; объемным джоулевым нагревом проволоки протекающим электрическим током; потерями

37

М. Ю. Харламов, И. В. Кривцун, В. Н. Коржик, А. И. Демьянов

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

Qv.

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

Y C v

/ w w w

дТ 1 д ( дТ Л д

(

дг r дг

X г----

V дг

дг,.

дТ

___м

дг,.

Л

+ J Р

(7)

"w V ^w J

Здесь Tw (r,zw) - пространственное распределение температуры в проволоке; Yw(T), Cw (Т), xw (Т) - плотность, удельная теплоемкость, коэффициент теплопроводности и удельное электрическое сопротивление материала проволоки соответственно; J - плотность электрического тока.

Граничные условия к уравнению (7) записываются следующим образом:

дТ

( дТ Л ( дТ Л

— X w I Aw л V дг J = Qc + Q - Qv ; г=Rw X w Aw ^ V Dzw j

= Q - Q - Qt; — = 0.

\Za Xdv Xdt дг

(8)

Zw =Lw

Описание граничных условий к уравнению (7) и замыкающих соотношений для определения тепловых потоков Qс, Qr, Qv и Qt приводится в [7].

Под действием дуги и обтекающего проволоку высокотемпературного плазменного потока проволока нагревается и на ее конце образуется слой расплавленного материала толщиной Lliq, некоторый объем которого будет сноситься в тонкую струю расплава обтекающим

проволоку плазменным потоком. В результате сноса части расплава с конца проволоки, условия теплового равновесия последней будут нарушены. Стремясь к равновесному состоянию, проволока займет такое пространственное положение относительно оси плазменной струи Lp = Z*, при котором объем удерживаемой на конце проволоки жидкой прослойки Vb будет

соответствовать объему расплавленного металла проволоки Vliq = nRWLHq, т. е. будет выполняться условие V, = Vig. Объем расплавленной части vu определяется по модели тепловых

процессов в проволоке [7], а выражение для объема удерживаемой на торце проволоки жидкой прослойки Vb, наряду с выражениями для оценки температуры, объема и скорости течения расплава в жидкой прослойке, удерживаемой на конце проволоки, приводятся в работе [8].

Модель формирования и срыва капель расплавленного металла. Для описания течения и распада струи расплавленного металла проволоки примем допущения об осесимметричности рассматриваемой гидродинамической системы, а также ламинарном режиме течения расплава. В этом случае можно показать, что течение струи жидкого металла с достаточно высокой точностью будет описываться в рамках системы одномерных уравнений Навье-Стокса для тонкой струи [9], записанной с учетом силы вязкого трения, действующего на расплав со стороны плазменного потока:

дv Dv 1 ф 1 д (, 2 Tv Л 2т

— + v— =---— + 3uw—— h — I-----,

дt дг Pw дг h дг V дг J hYw

дг D(vF)

— + ——- = 0,

дt дг

Ld < z <0;

(9)

(10)

где v = v(z, t) - аксиальная компонента скорости течения расплава; p - давление в расплаве;

h = h( г, t) - радиус поперечного сечения струи;

F (г, t) = nh 2( z , t) - площадь ее поперечного сечения; тp = тp (z) - напряжение трения на обтекаемой поверхности [8];

38

Комплексное математическое моделирование процессов плазменно-дугового проволочного напыления покрытий

Yw, uw - плотность и кинематическая вязкость расплавленного металла проволоки;

Ld - полная длина струи.

Величина давления в струе определяется выражением p = 2а K + p , где а - коэффициент поверхностного натяжения; K - средняя кривизна поверхности струи; pext - давление во внешней среде.

Постановка граничных и начальных условий для (9) и (10) описана в [9] и зависит от параметров удерживаемой на рабочем конце проволоки жидкой прослойки [8].

Отрыв капли в точке z = z фиксировался при выполнении условия h(z*) < h*, где h* ^ 0

. В этом случае объем отрывающейся капли определяется соотношением у * = п j h2( z, t )dz

- Ld

и предполагается выполнение следующих условий:

h (z* ) =0; Ld = z*.

(11)

Модель движения, нагрева и дробления капель расплавленного металла в плазменной струе. При моделировании поведения расплавленных частиц в плазменном потоке полагается, что координаты точки ввода жидких частиц в плазменный поток xo, zo совпадают с местоположением расплавленного конца проволоки. Начальные значения диаметра do и скорости w 0 движения жидкой частицы, определяются на основе модели струйного течения расплавленного металла проволоки и формировании капель жидкого металла в спутном высокоскоростном газовом потоке [9].

После отрыва капли ее движение в плазменной струе описывается уравнениями:

d (mw > = f, — = w. (12)

dt ’ dt

Здесь m(t), w(t) = (wx, Wy, wz) и r(t) = (x, y, z) - текущие значения массы (

m(0) = m = — 7td2p ), вектора скорости и радиус-вектора положения частицы в выбранной

V / 0 3 0 • m

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

F = 0,5 CdSp(v - w)|v - w|, (13)

где Cd - коэффициент лобового (аэродинамического) сопротивления;

S - площадь миделева сечения частицы;

v - вектор невозмущенной скорости плазмы в точке нахождения частицы; p - плотность плазмы.

При движении в плазменной струе жидкая частица деформируется. Полагалось что, деформируясь, частица принимает форму эллипсоида вращения (сфероида), геометрические размеры которого будем характеризовать безразмерным параметром y = d miijdv , представляющим собой отношение диаметра миделева сечения частицы dmid к диаметру эквивалентной по объему сферы dv. Уравнение для учета изменения формы жидкой частицы с краевыми условиями и замыкающими соотношениями приводится в [10].

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

Cd использовался интерполяционный многочлен, позволяющему вычислять Cd для эллипсоида по известным коэффициентам аэродинамического сопротивления для диска C , сферы

C, и сжатого эллипсоида с коэффициентом формы C, [10].

dSphere d0.5

39

М. Ю. Харламов, И. В. Кривцун, В. Н. Коржик, А. И. Демьянов

Тепловое состояние частиц при плазменно-дуговом напылении определялось при помощи нестационарного уравнения теплопроводности с учетом допущения о сферической форме частицы, с использованием текущего значения диаметра эквивалентной по объему сферы dV:

77 дТ

Р C -------m

г m m /-n .

dt

1A

r2 д r

f

V

2 дТ

rl у ____m

Am /-n

д r

(14)

Здесь Tm(r, t) - пространственно-временное распределение температуры в частице; %m (Т) и C (Т) - коэффициент теплопроводности и эффективная теплоемкость ее материала, определяемая как Cm = cm (Т) + W(m)S(T - Т(m)), где Cm(T) - удельная теплоемкость материала; Т^

- температура плавления; Wm) - скрытая теплота плавления; ^(х) - дельта-функция.

Краевые условия для уравнения (14) записываются с учетом конвективно-кондуктивного и радиационного тепловых потоков, а также удельных потерь тепла, связанные с поверхностным испарением материала частицы, и подробно приводятся в [10].

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

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

Выводы

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

Литература

1. Jorn Larsen-Basse. Surface engineering and the new millennium // Surface Engi-neering. -1998. - Vol. 14, No. 2. - P. 81-83.

2. Кудинов, В. В. Нанесение покрытий напылением. Теория, технология и оборудование [Текст] / В. В. Кудинов, Г. В. Бобров. - М. : Металлургия, 1992. - 432 с.

3. Технология и оборудование для плазменно-дугового напыления для восстановления ответственных деталей железнодорожного транспорта [Текст] / В. Н. Коржик, М. Ю. Харламов, С. В. Петров [и др.] // Вестник Восточноукр. национ. ун-та им. В. Даля. - 2011. - № 14. - С. 76-82.

4. Математическая модель дуговой плазмы, генерируемой плазмотроном с проволокой-анодом [Текст] / М. Ю. Харламов, И. В. Кривцун, В. Н. Коржик [и др.] // Автоматическая сварка. - 2007. - № 12. - С. 14-20.

40

Комплексное математическое моделирование процессов плазменно-дугового проволочного напыления покрытий

5. Влияние рода газа спутного потока на характеристики дуговой плазмы, создаваемой плазмотроном с проволокой-анодом [Текст] / М. Ю. Харламов, И. В. Кривцун, В. Н. Коржик [и др.] // Автоматическая сварка. - 2008. - № 6. - С. 19-24.

6. Об уточнении математической модели электрической дуги в плазмотроне с внешней токоведущей проволокой [Текст] / М. Ю. Харламов, И. В. Кривцун, В. Н. Коржик, С. В. Петров, А. И. Демьянов // Автоматическая сварка. - 2009. - № 1. - С. 53-56.

7. Нагрев и плавление проволоки-анода при плазменно-дуговом напылении [Текст] / М. Ю. Харламов, И. В. Кривцун, В. Н. Коржик, С. В. Петров // Автоматическая сварка. - 2011.

- № 5. - С. 5-11.

8. Формирование прослойки жидкого металла на торце проволоки-анода при плазменнодуговом напылении [Текст] / М. Ю. Харламов, И. В. Кривцун, В. Н. Коржик, С. В. Петров // Автоматическая сварка. - 2011. - № 12. - С. 3-8.

9. M. Yu. Kharlamov, I. V. Krivtsun and V. N. Korzhyk, Dynamic Model of the Wire Dispersion Process in Plasma-Arc Spraying, J. Therm. Spray Technol. - 2014. - 23 (3). - Р. 420-430.

10. M. Yu. Kharlamov, I. V. Krivtsun, V. N. Korzhyk, Y. V. Ryabovolyk, O. I. Demyanov, Simulation of Motion, Heating, and Breakup of Molten Metal Droplets in the Plasma Jet at Plasma-Arc Spraying, J. Therm. Spray Technol. - 2015. - 24 (4). - Р. 659-670.

11. Методы контроля температуры и скорости частиц конденсированной фазы в процессе плазменно-дугового напыления [Текст] / М. П. Бороненко [и др.] // Фундаментальные исследования. - 2013. - № 10-6. - С. 1194-1199.

12. Моделирование технологических процессов плазменного напыления покрытий наноразмерной толщины [Текст] / П. Ю. Гуляев, И. П. Гуляев // Системы управления и информационные технологии. - 2009. - № 1.1 (35). - С. 144-148.

13. Оценка скорости и температуры дисперсной фазы в струях плазменно-дугового напыления [Текст] / М. П. Бороненко [и др.] // Фундаментальные исследования. - 2014. - № 11-10.

- С. 2135-2140.

14. Экспериментальное исследование процесса плазменно-дугового проволочного напыления [Текст] / П. Ю. Гуляев, И. П. Гуляев, В. Н. Коржик [и др.] // Автоматическая сварка. -2015. - № 3-4. - С. 37-43.

41