Научный журнал
Современные наукоемкие технологии
ISSN 1812-7320
"Перечень" ВАК
ИФ РИНЦ = 1,021

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

Трифонов Г.И. 1, 2 Поленов В.С. 1 Жачкин С.Ю. 2
1 Военный учебно-научный центр Военно-воздушных сил «Военно-воздушная академия имени профессора Н.Е. Жуковского и Ю.А. Гагарина»
2 ФГБОУ ВО «Воронежский государственный технический университет»
В статье рассматриваются вопросы математического описания физических и химических преобразований композитного порошкового вещества под действием транспортирующего и плазмообразующего газов, а также электрической дуги. В качестве объекта исследования выбран плазменный поток с композитными частицами порошка. Стоит отметить, что зачастую использованию плазменных технологий препятствует недостаток систематизированной информации о закономерностях и возможностях соответствующих технологических процессов. Поэтому при планировании использования технологии плазменного напыления, в частности при прогнозировании кинематических характеристик плазменного потока, была разработана математическая модель, описывающая распространение упругих волн в трёхкомпонентной среде. Причем для более точного описания эксплуатационных свойств и параметров композитных материалов введено понятие – идеализированная среда. Целесообразность использования идеализированной среды заключается в том, что она в полной мере отражает изменения эксплуатационных и технологических параметров реальных сред в заданном интервале нагрузок и температур. В итоге смоделирована уточненная математическая модель на основе распространения нестационарных волн в трёхкомпонентной среде при образовании плазменного потока напыления. Также в ходе математического моделирования получена таблица теоретическими показателями тензоров напряжений и весовых коэффициентов процесса плазменного напыления. Построены графики зависимости скорости распространения продольных волн в материале плазменного покрытия от физико-механических характеристик – модуля Юнга и коэффициента Пуассона.
математическое моделирование
плазменное напыление
упругие волны
трёхкомпонентная среда
модуль Юнга
коэффициент Пуассона
1. Пузряков А.Ф. Теоретические основы технологии плазменного напыления. М.: Изд-во МГТУ им. Н.Е. Баумана, 2008. 360 с.
2. Трифонов Г.И., Жачкин С.Ю., Краснова М.Н., Пеньков Н.А. Моделирование кинематики плазменного напыления на сложнопрофильные поверхности // Труды ГОСНИТИ. 2017. Т. 128. С. 133–139.
3. Biot M.A. Theory propagation of elastic waves in a fluid-saturated porous solid I. Low-Frequency Range. J. Acoust. Soc. America. 1956. vol. 28. no 2. P. 168–178.
4. Biot M.A. Theory propagation of elastic waves in a fluid- saturated porous solid. II. Higher Frequency Range. J. Acoust. Soc. America. 1956. vol. 28. no 2. P. 179–191.
5. Biot M.A. Mechanics of Deformation and Acoustic Propagation in Porous Media. Journal of Applied Physic. 1962. vol. 33. no 4. P. 1483–1498.
6. Поленов В.С., Масликова Т.И. О распространении нестационарных упругих волн в однородных пористых средах // Изв. РАН. МТТ. 2005. № 1. С. 104–108.
7. Ильющенко А.Ф., Шевцов А.И., Оковитый В.А., Громыко Г.Ф. Процессы формирования газотермических покрытий и их моделирование. Минск: Беларус. навука, 2011. 357 с.
8. Поленов В.С. Основные уравнения динамического деформирования трехкомпонентных сред // Современные тенденции развития естествознания и технических наук: сборник научных трудов по материалам МНПК. Белгород, 2018. С. 31–34.
9. Соснин Н.А., Ермаков С.А., Тополянский П.А. Плазменные технологии. Руководство для инженеров. СПб.: Изд-во Политехи, ун-та, 2013. 406 с.
10. Кукарских Л.А., Белых А.Г. Волновые процессы в пористых средах // V МНПК. Белгород, 2015. № 5. ч. 1. С. 13–16.
11. Симонов Ю.Н., Георгиев М.Н., Симонов М.Ю. Основы физики и механики разрушения. Пермь: Изд-во Перм. нац. исслед. политехн. ун-та, 2012. 203 с.
12. Просветов Г.И. Линейная алгебра и аналитическая геометрия: задачи и решения. М.: Альфа-Пресс, 2009. 208 с.
13. Выгодский М.Я. Справочник по высшей математике. Изд.: Астрель, 2010. 992 с.

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

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

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

Цель исследования: разработка математической модели на основе распространения нестационарных волн в трехкомпонентной среде при образовании плазменной струи.

Материалы и методы исследования

Вопросу по распространению упругих волн в двухкомпонентных средах посвящено множество научных трактатов и работ. Среди них стоит выделить труды заслуженного ученого М.А. Био [3–5].

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

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

Доказано, что в такой среде в общем случае распространяется три волны, скорости которых существенным образом зависят от направления распространения волновой поверхности.

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

На базе идеализированной среды запишем систему уравнений, определяющую динамическое поведение упругой, насыщенной жидкостью и газом трехкомпонентной среды в перемещениях компонент [6], выражая коэффициенты Ламе λ, μ через модуль Юнга E и коэффициент Пуассона v по формулам [8]:

trif01.wmf, trif02.wmf

trif03.wmf (1)

– полный тензор напряжений в скелете при наличии жидкости и газа в порах:

trif72a.wmf

trif72b.wmf (2)

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

trif05.wmf (3)

trif06.wmf

– уравнения движения пористой среды:

trif07.wmf

trif08.wmf (4)

trif09.wmf

Здесь ρ11, ρ22, ρ33 – эффективные плотности порошковой механической смеси для плазменного напыления, плазмы и плазмообразующего газа соответственно (кг/м3); ρ12 < 0, ρ13 < 0, ρ23 < 0 – коэффициенты динамической связи скелета; R0(2), R0(3) – коэффициенты сжимаемости компонент, заполненных жидкостью и газом; uj(α) – перемещения компонент. В скобках цифры вверху обозначают: 1 – твердая компонента, 2 – жидкость, 3 – газ. Предполагается ρji, = ρij. По повторяющимся индексам здесь и в дальнейшем проводится суммирование от 1 до 3.

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

Продифференцируем соотношения (2) и (3) по t и получим

trif73a.wmf

trif73b.wmf

trif11.wmf (5)

trif12.wmf

Запишем уравнениям (4) и соотношения (5) разрывов [9]:

trif13.wmf

trif14.wmf

trif15.wmf (6)

trif16.wmf

trif17.wmf

trif18.wmf

К соотношениям (6) применим кинематические и геометрические условия совместности первого порядка на поверхности разрыва [10]:

trif19.wmf, trif20.wmf, trif21.wmf, trif22.wmf

trif23.wmf, trif24.wmf, trif25.wmf, trif26.wmf trif27.wmf (7)

Здесь sik,η, γ, λi(1), λi(2), λi(3) – величины, характеризующие скачки первых производных напряжений, сил, действующих на плазму и плазмообразующий газ и скоростей перемещения компонент; vi – единичный вектор нормали к волновой поверхности; G – скорость распространения волновой поверхности пористой среды.

Используя условия (7), формулы (6) запишем в виде

trif28.wmf

trif29.wmf

trif30.wmf (8)

trif31.wmf

trif32.wmf

trif33.wmf

Исключая из (8) величины sij, η, γ, получим однородную систему уравнений относительно λk(1), λk(2), λk(3) [11]:

trif34.wmf

trif35.wmf (9)

trif36.wmf

Полагая, что λk(α)vk ≠ 0, на волновой поверхности, умножим (9) на vi и просуммируем по повторяющемуся индексу i, получим однородную систему трех линейных уравнений относительно wα = λk(α)vk (α = 1, 2, 3):

trif37.wmf

trif38.wmf (10)

trif39.wmf

Введем следующие обозначения:

trif40.wmf (11)

trif41.wmf

С учетом (11) систему (10) запишем в безразмерной форме:

trif42.wmf

trif43.wmf (12)

trif44.wmf, trif45.wmf, trif46.wmf

Условие существования нетривиальных решений системы (12), однородной относительно w1, w2, w3 определяет три скорости волн ускорений в насыщенной жидкостью и газом пористой среде, которые находятся из определителя третьего порядка, составленного из коэффициентов при неизвестных w1, w2, w3 системы (12):

trif47.wmf (13)

Для вычисления определителя Δ справедлива формула разложения данного определителя по элементам i-го столбца [12]:

trif48.wmf trif49.wmf (14)

где i – номер столбца; j – номер строки; aij – элемент определителя, стоящий на пересечении i-й строки и j-го столбца; Mij – минор элемента матрицы третьего порядка.

Раскрывая определитель (13) по формуле (14) и сделав замену trif50.wmf, получим кубическое уравнение относительно z:

trif51.wmf,

trif52.wmf

trif53.wmf (15)

trif54.wmf

Коэффициенты σ11… σ33 находятся по формулам (11). Решение кубического уравнения (15) находим по формулам Кардана [13]. Для решения уравнения (15) разделим его на k и введем новую переменную trif55.wmf, тогда после преобразования получим:

trif56.wmf trif57.wmf trif58.wmf (16)

Вычислим дискриминант trif59.wmf. Если D < 0, то рассматриваемое уравнение обладает тремя действительными корнями, которые записываются как комплексные величины. Если D > 0, то уравнение имеет несколько различных решений, а именно одно действительное и два мнимых. Если D = 0, то имеем три совпавших нулевых корня или имеем три действительных решения, два из которых совпадают.

Таким образом, в трехкомпонентной пористой среде распространяются три продольные волны в зависимости от дискриминанта кубического уравнения. Если связь между компонентами жидкость – газ и упругость – газ в среде отсутствует γ13 = 0, γ23 = 0, σ23 = 0, то кубическое уравнение (15) сводится к квадратному уравнению:

trif60.wmf trif61.wmf, trif62.wmf, trif63.wmf (17)

Уравнение (17) совпадает с уравнением работы [6].

Положим в системе (9) λi(α)vi (α = 1, 2, 3).Тогда получим G = Gt:

trif64.wmf trif65.wmf

trif66.wmf (18)

trif67.wmf

trif68.wmf, trif69.wmf

Условием существования ненулевых решений системы (18) является ее определитель, составленный из коэффициентов при λi(α) (α = 1,2,3) и который должен быть равен нулю:

trif70.wmf (19)

Раскрывая уравнение (19) по формуле (14), получим выражение для нахождения скорости поперечной волны, распространяющейся в трехкомпонентной среде:

trif71.wmf (20)

Если связь между плазмой и плазмообразующим газом, композитным материалом и плазмообразующим газом отсутствует, γ23 = 0, γ13 = 0, то из (20) следует

trif72.wmf (21)

Формула (21) совпадает с формулой, полученной в работе [6].

Результаты исследования и их обсуждение

Согласно внесенным данным в идеализированную среду, а также результатам расчетов уравнений (11), получаем сводную таблицу значений.

Таблица значений тензоров напряжений и весовых коэффициентов

№ п/п

1

2

3

4

5

σ11

1

1

1

1

1

σ12

6,5E-13

6,14E-13

5,86E-13

5,63E-13

5,46E-13

σ22

1,15E-13

1E-13

8,75E-14

7,68E-14

6,75E-14

σ13

1,08E-12

1,06E-12

1,05E-12

1,01E-12

9,71E-13

σ23

1,08E-12

1,06E-12

1,05E-12

1,01E-12

9,71E-13

σ33

1,91E-13

1,73E-13

1,58E-13

1,37E-13

1,2E-13

γ34

1,049211

1,044258

1,039306

1,094341

1,061377

γ12

–0,01

–0,015

–0,02

–0,015

–0,012

γ22

0,020493

0,01647

0,030446

0,030423

0,017399

γ13

–0,04

–0,03

–0,02

–0,08

–0,05

γ23

–0,01

–0,001

–0,01

–0,015

–0,005

γ33

0,050296

0,031272

0,030248

0,095236

0,055224

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

trif1.wmf

Рис. 1. Зависимость коэффициента Пуассона от скорости распространения продольной волны

trif2.wmf

Рис. 2. Зависимость модуля Юнга от скорости распространения продольной волны

Выводы

Разработана математическая модель на основе распространения нестационарных волн в трехкомпонентной среде при образовании плазменной струи.

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


Библиографическая ссылка

Трифонов Г.И., Поленов В.С., Жачкин С.Ю. МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ПРОЦЕССОВ ПЛАЗМЕННОГО НАПЫЛЕНИЯ ПОКРЫТИЙ ТРЕХКОМПОНЕНТНЫХ СРЕД // Современные наукоемкие технологии. – 2018. – № 10. – С. 131-136;
URL: http://top-technologies.ru/ru/article/view?id=37208 (дата обращения: 20.06.2021).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.074