Исследование режимов работы серийного плазмотрона Metco F4 с использованием плазмообразующих газов аргона и азота
Получены аналитические зависимости полной мощности и напряжения плазменной дуги серийного плазмотрона MetcoF4 от тока дуги и расходов аргона и азота. Установлено, что увеличение расхода азота оказывает наибольшее влияние на увеличение напряжения дуги. Предложена расчетно-графическая методика оценки средневзвешенной температуры и скорости плазменного потока на срезе сопла серийного плазмотрона MetcoF4.
Введение
В настоящее время для нанесения защитных покрытий на крупногабаритные детали ГТД применяются установки плазменного, детонационного и высокоскоростного напыления [1–8].
При плазменном напылении между вольфрамовым катодом и медным водоохлаждаемым анодом (соплом) возбуждается электрическая дуга, нагревающая рабочий газ, в результате чего он истекает из сопла в виде плазменной струи. Порошковый материал, подаваемый в плазменную струю потоком транспортирующего газа (аргона или азота), с большой энергией ударяется о поверхность обрабатываемой детали. Благодаря высоким температурам при плазменном напылении можно наносить тугоплавкие материалы. Кроме того, регулирование температуры и скорости плазменной струи позволяет расширить класс наносимых материалов. На рис. 1 представлена принципиальная схема плазмотрона для нанесения покрытий методом плазменного напыления [9].

Рис. 1. Принципиальная схема плазмотрона
Плазменное напыление в зависимости от давления в камере напыления разделяют на атмосферное плазменное напыление (APS – Air Plasma Spray), плазменное напыление при пониженном давлении (LPPS – Low Pressure Plasma Spray) и вакуумное плазменное напыление (VPS – Vacuum Plasma Spray). Плазменными методами наносят широкий класс материалов: металлические сплавы систем Me–Cr–Al–Y, металлокерамические композиции типа MeC–NiCr/Co, покрытия на основе оксидов типа Al2O3/ZrO2–Y2O3. Оборудование для атмосферного плазменного напыления фирм Sulzer Metco и AMT используется на предприятиях авиадвигателестроения «НПЦ газотурбостроения «Салют», ПАО «ОДК-УМПО», АО «ОДК-ПМ» и других научно-производственных объединениях.
Наиболее распространенным плазмотроном для нанесения покрытий методом атмосферного плазменного напыления является плазмотрон Metco F4 производства фирмы Sulzer Metco (в настоящее время – фирма Oerlikon Metco), который представлен на рис. 2.

Рис. 2. Плазмотрон Metco F4 (а), электрод (б) и сопло (в)
При нанесении покрытий из различных порошковых материалов варьируются следующие параметры процесса: мощность дуги (кВт), расход плазмообразующего газа при атмосферном (нормальном) давлении (н.л/мин), расход транспортирующего газа (н.л/мин), количество подаваемого материала (г/мин), дистанция напыления (мм), скорость перемещения пятна напыления (мм/с). С увеличением мощности дуги возрастают температура и энтальпия плазменной струи. Мощность дуги определяется силой тока и напряжением, которое, в свою очередь, определяется видом и расходом плазмообразующего газа [9].
В работе проведено исследование режимов работы плазмотрона, определены зависимости мощности, подводимой к плазмотрону, от тока дуги и расхода плазмообразующих газов аргона и азота.
Работа выполнена в рамках реализации комплексных научных проблем 17.3. «Многослойные жаростойкие и теплозащитные покрытия, наноструктурные упрочняющие эрозионно- и коррозионностойкие, износостойкие, антифреттинговые покрытия для защиты деталей горячего тракта и компрессора ГТД и ГТУ» и 17.4. «Многослойные защитные покрытия и плазмохимическое оборудование для осаждения защитных и упрочняющих покрытий из газовых потоков плазмы, содержащих прекурсоры элементов синтезируемого покрытия» («Стратегические направления развития материалов и технологий их переработки на период до 2030 года») [1].
Материалы и методы
Объектом исследования в данной работе является процесс атмосферного плазменного напыления с использованием серийного плазмотрона Metco F4 и плазмообразующих газов аргона и азота. Проведено исследование режима работы плазмотрона с варьированием следующих параметров: ток дуги 300–500 А с шагом 50 А, расход аргона 15–25 н.л/мин с шагом 5 н.л/мин, расход азота 0–15 н.л/мин (0, 11 и 15 н.л/мин в связи с технологическими ограничениями самой установки плазменного напыления). На каждом режиме работы фиксировали мощность P, которую подводили к плазмотрону. В свою очередь – P=IU, где I – ток дуги, А; U – напряжение, В. Результаты исследования приведены в табл. 1.
Таблица 1
Мощность плазмотрона при различных режимах работы
Ток дуги I, А | Мощность плазмотрона Р, кВт, при расходе газа, н.л/мин | |||
азота QN | аргона QAr | |||
15 | 20 | 25 | ||
300 | 0 | 8,1 | 8,6 | 9,1 |
11 | 21,2 | 21,7 | 22,2 | |
15 | 23,9 | 24,3 | 24,5 | |
350 | 0 | 9,3 | 10 | 10,7 |
11 | 24 | 24,6 | 25,2 | |
15 | 27 | 27,7 | 28,5 | |
400 | 0 | 10,8 | 11,5 | 12,1 |
11 | 26,5 | 26,8 | 27,2 | |
15 | 29,7 | 30,5 | 31,4 | |
450 | 0 | 12,1 | 13 | 13,7 |
11 | 29,1 | 29,3 | 29,7 | |
15 | 32,2 | 32,6 | 33,4 | |
500 | 0 | 13,5 | 14,5 | 15,3 |
11 | 31,3 | 31,9 | 32,5 | |
15 | 34,6 | 35,4 | 35,9 | |
Полагая, что полная мощность плазмотрона складывается из индивидуальных вкладов мощности от использования аргона и азота, т. е.:
P=PAr+PN, (1)
из полной мощности следует вычесть мощность дуги при расходе азота 0 н.л/мин. Результат вычислений представлен в табл. 2.
Таблица 2
Мощность плазмотрона при использовании одного плазмообразующего газа – азота
Ток дугиI, А | Мощность плазмотрона Р, кВт, при расходе газа, н.л/мин | Среднее значение PN, кВт | |||
азота QN | аргона QAr | ||||
15 | 20 | 25 | |||
300 | 0 | 0 | 0 | 0 | 0 |
11 | 13,1 | 13,1 | 13,1 | 13,1 | |
15 | 15,8 | 15,7 | 15,4 | 15,633 | |
350 | 0 | 0 | 0 | 0 | 0 |
11 | 14,7 | 14,6 | 14,5 | 14,6 | |
15 | 17,7 | 17,7 | 17,8 | 17,733 | |
400 | 0 | 0 | 0 | 0 | 0 |
11 | 15,7 | 15,3 | 15,1 | 15,367 | |
15 | 18,9 | 19 | 19,3 | 19,067 | |
450 | 0 | 0 | 0 | 0 | 0 |
11 | 17 | 16,3 | 16 | 16,433 | |
15 | 20,1 | 19,6 | 19,7 | 19,8 | |
500 | 0 | 0 | 0 | 0 | 0 |
11 | 17,8 | 17,4 | 17,2 | 17,467 | |
15 | 21,1 | 20,9 | 20,6 | 20,867 | |
Целью исследования является выявление зависимостей индивидуальных мощностей PAr и PN от тока дуги и расходов газов аргона и азота. В связи с этим применены методы аппроксимации линеаризованных функций [10].
Если массив экспериментальных данных X (X1, X2, …, Xi, … Xn), Y (Y1, Y2, …, Yi, … Yn) описывается функцией вида f(g)=agb, то ее преобразуют в следующий вид:
ln(f(g))=b·ln(g)+lna, (2)
что соответствует виду y(x)=ax+b.
Затем коэффициенты в выражении (2) вычисляют по методу наименьших квадратов как:
(3)
(4)
где xi=ln(Xi), yi=ln(Yi).
Если массив экспериментальных данных X (X1, X2, …, Xi, …Xn), Y (Y1, Y2, …, Yi, … Yn) описывается функцией вида f(g)=agk, где k – константа, то ее аппроксимируют методом наименьших квадратов следующим образом:
(5)
(6)
Достоверность аппроксимации определяется коэффициентом детерминации R2, который, в свою очередь, является квадратом коэффициента корреляции rx, y:
(7)
(8)
где f(Xi) – значение функции аппроксимирующей зависимости в точке Xi.
Коэффициент корреляции rx, y характеризует тесную взаимосвязь экспериментальных параметров с уравнением регрессии. По шкале Чеддока значения 0,7≤rx, y≤0,9 соответствуют высокой корреляционной связи, при 0,9≤rx, y≤0,99 – очень высокой [10].
Результаты и обсуждение
Расчет зависимости мощности плазмотрона при использовании
плазмообразующего газа аргона
При расчете используются данные из табл. 1 при расходе азота QN=0. Очевидно, что при расходе аргона QAr=0, мощность, подводимая к плазмотрону, также будет равна PAr=0. С учетом этого предполагается, что зависимость мощности плазмотрона от расхода аргона при постоянном токе дуги будет иметь вид:

Предварительный расчет показал, что коэффициент b≈0,25, тогда с учетом того, что
методом наименьших квадратов по формулам (5) и (6) получим зависимости мощности плазмотрона Metco F4 при использовании плазмообразующего газа аргона при разных значениях тока дуги. Результаты расчетов представлены в табл. 3.
Таблица 3
Результаты расчетов зависимости мощности дуги от расхода аргона 
Ток дуги I, А | Коэффициент a | Коэффициент детерминации R2 | Коэффициент корреляции rx, y |
300 | 4,082 | 0,987 | 0,994 |
350 | 4,749 | 0,989 | 0,995 |
400 | 5,442 | 0,985 | 0,992 |
450 | 6,140 | 0,999 | 0,999 |
500 | 6,852 | 1 | 1 |
Зависимости мощности плазмотрона от расхода аргона при токах дуги 300–500 А с шагом 50 А и их уравнения регрессии представлены на рис. 3.
Рис. 3. Зависимость мощности дуги от расхода аргона при токе дуги 300 (♦); 350 (■);
400 (▲); 450 (■) и 500 А (■)
Из данных табл. 3 видно, что коэффициент a монотонно возрастает при увеличении тока дуги. Полагая, что при токе дуги 0 А мощность также составит 0 кВт, зависимость a(I) будет иметь вид f(g)=agb. Следовательно, уравнение регрессии будет найдено методом наименьших квадратов по формулам (3) и (4):
a(I)=0,012·I1,015;
R2=1;
rx, y=1.
Очевидно, что зависимость мощности плазмотрона от тока дуги и расхода аргона находится из выражения:
(9)
Зависимость (9) представлена на рис. 4.

Рис. 4. Зависимость мощности плазмотрона от тока дуги и расхода аргона
Расчет зависимости мощности плазмотрона при использовании
плазмообразующего газа азота
При расчете используются средние значения мощности PN из табл. 2 при расходах азота QN, н.л/мин: 0; 11 и 15. Очевидно, что при расходе азота QN=0, мощность, подводимая к плазмотрону, также будет равна PN=0. С учетом этого предполагается, что зависимость мощности плазмотрона от расхода азота при постоянном токе дуги будет иметь вид:

Предварительный расчет показал, что коэффициент b≈0,6, тогда с учетом того, что методом наименьших квадратов по формулам (5) и (6) получим зависимости мощности плазмотрона Metco F4 при использовании плазмообразующего газа азота при разных значениях тока дуги. Результаты расчетов представлены в табл. 4.
Таблица 4
Результаты расчетов зависимости мощности дуги от расхода азота 
Ток дуги I, А | Коэффициент a | Коэффициент детерминации R2 | Коэффициент корреляции rx, y |
300 | 3,091 | 0,997 | 0,999 |
350 | 3,481 | 0,998 | 0,999 |
400 | 3,710 | 0,982 | 0,991 |
450 | 3,899 | 1 | 1 |
500 | 4,126 | 0,998 | 0,999 |
Зависимости мощности плазмотрона от расхода азота при токах дуги 300–500 А с шагом 50 А и их уравнения регрессии представлены на рис. 5.

Рис. 5. Зависимость мощности дуги от расхода азота при токе дуги 300 (♦); 350 (■);
400 (▲); 450 (■) и 500 А (■)
Из данных табл. 4 видно, что коэффициент a монотонно возрастает при увеличении тока дуги. Полагая, что при токе дуги 0 А мощность также составит 0 кВт, зависимость a(I) будет иметь вид f(g)=agb. Следовательно, уравнение регрессии будет найдено методом наименьших квадратов по формулам (3) и (4):
a(I)=0,139·I0,547;
R2=0,986;
rx, y=0,993.
Очевидно, что зависимость мощности плазмотрона от тока дуги и расхода азота находится из выражения:
(10)
Зависимость (10) представлена на рис. 6.
Рис. 6. Зависимость мощности плазмотрона от тока дуги и расхода азота
Очевидно, что полная мощность серийного плазмотрона Metco F4 при использовании плазмообразующих газов аргона и азота вычисляется как:
(11)
Исследование зависимости полной мощности плазмотрона при использовании
плазмообразующих газов аргона и азота
Оценка адекватности математической модели проведена путем вычисления относительной погрешности εотн:

где Pэксп, Ррасч – соответственно экспериментальные значения мощности дуги из табл. 1 и расчетные значения мощности дуги, полученные по формуле (11).
Результаты оценки представлены в табл. 5.
Таблица 5
Оценка адекватности математической модели, описываемой формулой (11)
Ток дуги I, А | Относительная погрешность модели, %, при расходе газа, н.л/мин | |||
азота QN | аргона QAr | |||
15 | 20 | 25 | ||
300 | 0 | 4,721 | 3,569 | 3,639 |
11 | 1,007 | 0,636 | 0,731 | |
15 | 0,835 | 0,099 | -1,028 | |
350 | 0 | 2,960 | 3,023 | 4,168 |
11 | 2,246 | 1,894 | 2,022 | |
15 | 2,171 | 2,214 | 3,007 | |
400 | 0 | 4,309 | 3,432 | 2,955 |
11 | 2,397 | 0,614 | -0,267 | |
15 | 2,218 | 2,256 | 3,028 | |
450 | 0 | 3,744 | 3,727 | 3,405 |
11 | 3,056 | 0,753 | -0,328 | |
15 | 1,868 | 0,408 | 0,643 | |
500 | 0 | 3,988 | 3,944 | 3,744 |
11 | 2,530 | 1,334 | 0,697 | |
15 | 1,454 | 0,951 | 0,105 | |
Из данных табл. 5 видно, что максимальная относительная погрешность модели составляет не более 5%, следовательно, математическая модель достаточно точно описывает зависимость мощности дуги серийного плазмотрона Metco F4 от тока дуги в рабочем диапазоне 300–500 А и расходов аргона и азота в диапазонах 15–25 и 11–15 н.л/мин соответственно.
Известно [9, 11–16], что для наиболее эффективного использования энергии при плазменном напылении необходимо выбирать режим работы плазмотрона с максимальным напряжением. С учетом того, что P=UI, из выражения (11) получаем:
(12)
Как несложно заметить, U(I, QAr, QN) является функцией трех переменных, поэтому при ее исследовании попеременно фиксируют два параметра и рассматривается зависимость функции от третьего параметра. В качестве фиксированных параметров принимаются значения середины диапазона вариации:
н.л/мин,
н.л/мин, тогда

Очевидно, скорость изменения функции оценивается по величине первой производной, тогда:

На рис. 7 показаны зависимости функций U(I), U(QAr), U(QN), а также их первые производные.
Рис. 7. Зависимости функций U(I) (а); dU(I)/dI(б); U(QAr) (в); dU(QAr)/dQAr (г); U(QN) (д); dU(QN)/dQN (е)
Видно (рис. 7, а, б), что увеличение тока дуги при сохранении значений остальных параметров снижает напряжение дуги, а повышение расхода аргона увеличивает напряжение дуги (рис. 7, в, г). Видно также (рис. 7, д, е), что повышение расхода азота увеличивает напряжение дуги. Сравнивая значения первых производных функций U(QAr) и U(QN) на рис. 7, г, е, видим, что расход азота вносит больший вклад в увеличение напряжения дуги, чем расход аргона.
Расчет температуры плазменного потока
Образование плазмы из молекул двухатомного газа представлено на рис. 8. При низких температурах до 10 К молекулы совершают поступательные движения (а), при повышении температуры до 1000 К у молекул появляется вращательное движение (б). Столкновения молекул провоцируют их диссоциацию (в), и при дальнейшем повышении температуры происходит их ионизация (г) [11].

Рис. 8. Механизм образования плазмы:
1 – положительный ион; 2 – электрон
Типичное распределение температур в плазменной струе в зависимости от расстояния от среза сопла показано в работе [12] при силе тока дуги 400 А и расходе аргона 10 н.л/мин и представлено на рис. 9.

Рис. 9. Распределение температуры в плазменной струе
Известно [12], что эффективность использования энергии плазменной дуги на нагрев плазмообразующего газа составляет величину ~30%. Среднеинтегральная объемная энтальпия определяется как [15]:
[Дж/м3],
где ΔE=0,3·P·t – энергия, передаваемая газу при нагреве в течение времени t с учетом эффективности процесса, Дж;
– объем газа в течение времени t, проходящий через срез сопла, м3.
Тогда:
(13)
(14)
Энтальпия газа, вычисленная по формулам (13) и (14), определяет индивидуальные температуры аргона TAr и азотаTN по графику [11], представленному на рис. 10.

Рис. 10. Зависимость энтальпии от температуры для аргона и азота в областях диссоциации (1) и ионизации (2)
Среднюю температуру аргон-азотной плазмы TAr–N можно приблизительно оценить из уравнения теплосодержания системы:
(15)
где CAr=523 Дж/(кг·К) – теплоемкость аргона при нормальных условиях; CN=1051 Дж/(кг·К) – теплоемкость азота при нормальных условиях; ρAr=1,784 кг/м3 – плотность аргона при нормальных условиях; ρN=1,25 кг/м3 – плотность азота при нормальных условиях.
При режиме работы плазмотрона
н.л/мин,
н.л/мин приблизительная температура аргон-азотной плазмы составляет TAr–N=7305 К.
Среднемассовая скорость потока vAr–N плазменной струи на срезе сопла может быть определена из уравнения неразрывности течения и уравнения состояния [15]:

где
, м/с – скорость истечения из сопла плазмотрона Metco F4 диаметром dс=6 мм холодного газа; Tхол=293 К – температура холодного газа.
Тогда:
(16)
Из уравнения (16) c учетом TAr–N=7305 К получаем среднемассовую скорость потока плазменной струи vAr–N=486,3 м/с.
Заключения
В результате исследования режимов работы серийного плазмотрона Metco F4 с использованием плазмообразующих газов аргона и азота получены следующие результаты.
1. Получена аналитическая зависимость полной мощности плазменной дуги от тока дуги и расходов аргона и азота. Относительная погрешность предложенной модели составляет не более 5%.
2. Получена аналитическая зависимость напряжения плазменной дуги от тока дуги и расходов аргона и азота. В результате исследования зависимости и ее частных производных установлено, что увеличение тока дуги снижает напряжение дуги, а повышение расходов аргона и азота повышает напряжение дуги, причем расход азота вносит больший вклад в повышение напряжения дуги.
3. Предложена расчетно-графическая методика оценки средневзвешенной температуры и скорости плазменного потока на срезе сопла серийного плазмотрона Metco F4.
- Каблов Е.Н. Инновационные разработки ФГУП «ВИАМ» ГНЦ РФ по реализации «Стратегических направлений развития материалов и технологий их переработки на период до 2030 года» // Авиационные материалы и технологии. 2015. №1 (34). С. 3–33. DOI: 10.18577/2071-9140-2015-0-1-3-33.
- Каблов Е.Н., Мубояджян С.А. Теплозащитные покрытия с керамическим слоем пониженной теплопроводности на основе оксида циркония для лопаток турбины высокого давления перспективных ГТД // Сб. докл. конф. «Современные достижения в области создания перспективных неметаллических композиционных материалов и покрытий для авиационной и космической техники». М.: ВИАМ, 2015. С. 3.
- Каблов Е.Н. Основные итоги и направления развития материалов для перспективной авиационной техники // 75 лет. Авиационные материалы. М.: ВИАМ, 2007. С. 20–26.
- Мубояджян С.А., Будиновский С.А., Гаямов А.М., Матвеев П.В. Высокотемпературные жаростойкие покрытия и жаростойкие слои для теплозащитных покрытий // Авиационные материалы и технологии. 2013. №1. С. 17–20.
- Будиновский С.А., Чубаров Д.А., Матвеев П.В. Современные способы нанесения теплозащитных покрытий на лопатки газотурбинных двигателей (обзор) // Авиационные материалы и технологии. 2014. №S5. С. 38–44. DOI: 10.18577/2071-9140-2014-0-s5-38-44.
- Матвеев П.В., Будиновский С.А., Мубояджян С.А., Косьмин А.А. Защитные жаростойкие покрытия для сплавов на основе интерметаллидов никеля // Авиационные материалы и технологии. 2013. №2. С. 12–15.
- Матвеев П.В., Будиновский С.А. Исследование свойств защитных жаростойких покрытий для интерметаллидных никелевых сплавов типа ВКНА для рабочих температур до 1300°С // Авиационные материалы и технологии. 2014. №3. С. 22–26. DOI: 10.18577/2071-9140-2014-0-3-22-26.
- Кашин Д.С., Стехов П.А. Защитные покрытия для жаропрочных сплавов на основе ниобия // Труды ВИАМ: электрон. науч.-технич. журн. 2015. №6. Ст. 01. URL: http://www.viam-works.ru (дата обращения 16.04.2018). DOI: 10.18577/2307-6046-2015-0-6-1-1.
- Ильющенко А.Ф., Шевцов А.И., Оковитый В.А., Громыко Г.Ф. Процессы формирования газотермических покрытий и их моделирование. Минск: Беларуская навука, 2011. 357 с.
- Кибзун А.И., Горянова Е.Р., Наумов А.В., Сиротин А.Н. Теория вероятностей и математическая статистика. Базовый курс с примерами и задачами. М.: Физматлит, 2002. 224 с.
- Хасуи А., Моригаки О. Наплавка и напыление. М.: Машиностроение, 1985. 240 с.
- Хасуй А. Техника напыления. М.: Машиностроение, 1975. 288 с.
- Кудинов В.В. Плазменные покрытия. М.: Наука, 1977. 184 с.
- Газотермическое напыление: учеб. пособие / под общ. ред. Л.Х. Балдаева. М.: Маркет ДС, 2007. 344 с.
- Кудинов В.В., Бобров Г.В. Нанесение покрытий напылением. Теория, технология и оборудование: учеб. для вузов. М.: Металлургия, 1992. 432 с.
- Сидоров А.И. Восстановление деталей машин напылением и наплавкой. М.: Машиностроение, 1987. 192 с.
