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

ТЕХНОЛОГИЯ ЭЛЕКТРОХИМИЧЕСКИХ ПРОИЗВОДСТВ

УДК 682.1

Л. М. Котляр, Г. Ю. Шакирова, И. Х. Исрафилов, М. Ф. Шаехов

МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ДИСПЕРГИРОВАНИЯ МЕТАЛЛОВ В ЭЛЕКТРИЧЕСКОМ РАЗРЯДЕ С ДВУМЯ ЭЛЕКТРОДАМИ - АНОДАМИ И ЖИДКИМ КАТОДОМ

Ключевые слова: катод, анод, электролит, диспергирование, плотность тока, потенциал, комплексный потенциал.

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

Keywords: cathode, anode, electrolyte, dispersing the current density, potential, the potential of the complex.

In this paper a mathematical model of the plane problem of dispersion of one or two metal rods in an electric discharge with a liquid electrolyte has been proposed. In solving the problem, the analogy with the plane problem of the theory of jet flows of an ideal fluid, which is reduced to the definition of an analytic function in a region with one section of the unknown border was used. An approximate formula for determining the mass of dispersed metal powder, assuming axial symmetry of the rod in the special case of one electrode was derived.

В последние годы успешно развивается применение ферромагнитных порошков для изготовления деталей повышенной прочности. Композиционные порошки получают синтезом в режиме технологического горения с последующими операциями размола и классификации по фракциям. Одним из методов обработки металлических материалов являются разряды в жидкостях [1, 2], которым можно найти применение в диспергировании металлов и сплавов. Установка для синтеза ферромагнитного порошка состоит из источника постоянного тока, электролитической ванны с электролитом и токопроводящего стального электрода с различным содержанием углерода (от 0,2 до 1 ). Результаты проведённых испытаний

показали большую зависимость количества и качества (размер гранул) порошка от многих геометрических и физических параметров.

Для решения задачи отметим, прежде всего, что в процессе образования порошка зависимость выхода по току ^ от анодной плотности тока ]а идентична аналогичной зависимости ]а) для

электрохимической обработки металлов [3, 4]. Исходя из этого, будем считать, что скорость снятия металла Ут с поверхности анода на единицу массы, определяется по аналогии с законом Фарадея ¥т = ]а^е , где ^ = ]а) - выход по току, равный доле энергии затраченной на образование порошка в электрическом разряде, ]а - плотность тока, е -коэффициент, соответствующий электрохимическому эквиваленту металла и определяется экспериментально. Процесс диспергирования можно разделить на две стадии: на первой, нестационарной, происходит диспергирование на

V и анода

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

Va = v • cos 8

где 8 - угол между вектором скорости подачи анода и единичным вектором na внешней нормали

к поверхности анода. В этом случае общая схема процесса не изменяется со временем и процесс можно считать установившемся или стационарным.

Из (1) установившееся распределение плотности тока ja на стационарной границе анода определяется соотношением

Ц( ja ) j

PV с = — cos 8.

где р - плотность материала анода.

Зависимость выхода по току ]а) можно аппроксимировать уравнением гиперболы и считать, что на границе анода выполняется соотношение [5]

ру /

]а =-(а

а0 + Ъ0соъ 8)

е

В работе [6] была решена задача о диспергировании стержня, расположенного в центре межэлектродного промежутка (МЭП).

Рассмотрим двумерную модель процесса с двумя электродами, расположенными симметрично относительно оси симметрии МЭП с жидким электродом (рис. 1).

В этом случае, в силу симметрии относительно оси А1С, можно рассматривать только одну, например, правую половину течения. Эта

ь

задача отличается от рассмотренной в работе [7] отсутствием симметрии течения в МЭП. Учитывая это, рассмотрим, следуя работе [7], более простую задачу о диспергировании одного электрода, расположенного симметрично в МЭП. Введём декартовую систему координат х1, у1, связанную с анодом, который движется в направлении оси ординат.

A, A

A

Bi Bi

с

D

Рис. 1 - Область МЭП

Будем считать, пренебрегая

приэлектродными явлениями, что в МЭП существует потенциал электрического поля у1, удовлетворяющий уравнению Лапласа

Ду1 = 0

и на границах электродов выполняется условие постоянства потенциалов у1а = иа , у1к = ик . В силу уравнения (3), существует функция ф1 гармонически сопряжённая у1 и можно ввести комплексный потенциал электростатического поля Ж (х1, у1) = ф1 (х1, у1) + /у 1 (х1, у1), являющейся

аналитической функцией в области 21 = х1 + у

Введём характерные значения плотности тока ¡0 = ру/8 , длины Н = к(иа - ик )/;о (к -удельная электропроводность среды) и перейдём к безразмерным переменным

х1 у1 х = —, у = —, 2 = х + / у ,

НН

Ж - /и к

Ж = ф + /у = —1--

иа - ик

Тогда, с учётом (2), функция у

удовлетворяет в межэлектродной области уравнению Лапласа и граничным условиям на поверхности электродов

где

V fl = 1, V k = 0 b

dV а

дп Jo - постоянные,

= J = а0 + b0 cos 8

a0, ь0 - иисшяпимь, учитывающие зависимость выхода по току от плотности тока.

Исходя из этого, задача моделируется плоскопараллельным потенциальным течением идеальной несжимаемой жидкости. При этом, аналогом напряжённости электрического поля E является скорость фиктивного течения V и векторы E и V - ортогональны. Вдоль линии у = const выполняется равенство ду/дп = V , V = |V|. На поверхности анода V = a0 + b0 cos 6 .

На рис. 1, в плоскости z = х + iy, представлена область МЭП ограниченная линиями симметрии A1C и B1D, CD - граница катода, A1ABB1 - граница анода, A1, B1 - бесконечно удалённые точки. Граница анода разделена на две области: прямолинейные участки A1A , B1B , где

Va = 0 и участок AB на котором выполняется условие (2), аналогично симметричной задаче [7].

Функция у удовлетворяет уравнению

Лапласа Ау = 0, на границах электродов

выполняются условия:

CD

= 0,

lA.ABB,

= 1,

ду/дп| = a0 + b0 cos 6 , на линиях симметрии A1C и B1D - условия ду/дп|А C = 0 , оу/дп|в£) = 0 .

Для рассматриваемой симметричной задачи, поток фиктивного течения создаётся системой непрерывно распределённых источников и стоков: источники - вдоль линии A1C, стоки - вдоль линии B1D. На неизвестной границе AB выполняется условие V = a0 + b0 cos 6 .

Пусть в плоскости вспомогательного комплексного переменного u = § + /т области течения Dz соответствует область

Du =ju| < 1, т> 0} (рис. 2) и функция z(u) конформно отображает область Du на область Dz с соответствием точек, указанном на рис. 1, 2.

Рис. 2- Область Б„

Определим две функции: комплексный потенциал фиктивного течения Ж(и) = ф(и) + /у(и) и функцию Жуковского [5]

Х

где г (и) = 1п(к0/ V), V - модуль скорости фиктивного течения, 0 - угол наклона вектора скорости к оси х, V, = а0 - значение скорости в точке А .

Комплексный потенциал Ж

удовлетворяет граничным условиям

У(и) = |°, " =*, С, ,

[1, и = £, а 1; и = е ,0 <ст<тс

ф(и) =

0, u = £, с <£< d ф0, u = - d < £ < -с

B

Область изменения функции Ж (и) -прямоугольник БЖ = {р + 0 < ф < ф0, 0 < у < 1}.

Рис. 3 - Область

Для определения производной йЖ ¡йи комплексного потенциала отобразим область Би на верхнюю полуплоскость с соответствием точек, указанном на рис. 2, 3, преобразованием

ю(и) =

2и

1 + и2

и, используя формулу Кристоффеля-Шварца [8], найдём функцию, отображающую область Бю на область изменения функции Ж (и):

Ж (ю) = N

у V С

йю

ю2 -а2)

у а

Ф0 =

2К

К (к) = |

уК

1

йю

0 \>Ю2 - 1)

где К(к) - полный эллиптический интеграл первого рода [9]. Дифференцируя Ж (ю) по переменной ю и учитывая зависимость (5), получим:

йЖ йи

2N

^4и2 - а2

Для функции Жуковского х(и) имеем следующие граничные условия

8(и) =

л/2, и -1 <§<-а 0, и ||| < а , -л/2, и а <§< 1

Ъ0

г

и = е/О,0 <ст<л .

Функцию Жуковского х(и) будем искать в виде х(и) = Х0(и) + / (и), где функция X 0(и) = г0(и) + /8 0(и) удовлетворяет граничным условиям

80(и) = 8(и), и 1,

г0(и) = 0, и = е/О,0 <ст<л и имеет в области Би те же особенности, что и Х(и), функция /(и) - аналитическая в Би и

непрерывная в Би. Построим %0(и) методом особых точек [5]:

1

(-.2 „2 Л

Х0(и) = -1п

22 ч 1 - а и

+ А0и + В0 .

Сравнивая граничные условия для функций х

Х(и) = - 1п

1 + -^соб 8

и = е/О,0 <ст<л .

Учитывая граничные условия

можно аналитически продолжить на весь круг и представить в виде ряда Лорана:

с2пи

2п

Г (и) = ^ с

п=0

где с2п - действительные коэффициенты. Из

ад

условия / (1) = 0 находим с0 = с

2п

Условие (7), с учётом представления функции / (и) в виде ряда, имеет вид

с0 + с

^С2п С0Б(2ПСТ) = - 1п

п=1

1 + -^соб 8

Умножая обе части равенства (8) на соб(2пст) , интегрируя по О в пределах 0, л/2, получим бесконечную систему уравнений для определения коэффициентов с2п :

42 ( Ъ Л

— Г 1п| 1 + —соб(80(ст) + |а(ст)) С0Б(2пст)йст

л 0 I а0 )

с2п =--

п = 1, ад.

Безразмерные координаты точек границы анода АВ определяем из (4) по формуле

1 и <1Ж

2

а0 1 йи

Для решения задачи необходимо определить математические параметры с и а. Это можно сделать, задав расстояния |СО| и |А1В^ : ¿(-с) - 1{с) = |С£| , г(-а) - ¿(а) = \АХВХ\.

Массу диспергированного порошка можно определить, приближенно считая стержень осесимметричным, используя формулу для вычисления площади поверхности вращения.

V 2

т = 4ла0р^ 2 Г х

й2 йо

йо

где р - плотность металла.

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

а

0

ю

п=1

а

0

а

0

Для задачи по схеме (рис. 1) решение строится аналогично предыдущему, но оно зависит от большего числа параметров. Область вспомогательного переменного и = § + /т будет иметь вид (рис.2) и, в силу отсутствия симметрии, в задаче появятся 3 параметра (а, Ь, с), которые определяются из задания геометрических размеров на схеме (рис.1).

Некоторые результаты расчета формы криволинейного участка электродов АВ представлены на рис. 4.

.....L .J..... -4.....1.....\..... ! ! ! .....[..... L

!

_____!_____j_____!_____L_____I_____j_____ .....L...

i i i l i ; .....i

VJ ! i ! L-- .....r..... 3

: : : 4.....1.....i.....

..... .....;""'

Г " "1"" """

Рис. 4 - Результаты расчётов формы границы анода Н1 - ширина анода (Н1 = |АА|), Н2 -ширина катода (Н2 = ), Ь - ширина МЭП в окрестности точки А1 а = 0.8 , Ь = 0.8, Н1 = 0.5 ,

Н2 = 2 1. ь = 0.5; 2 Ь = 0.3; 3. Ь = 0.2

электрода не симметрична относительно его оси (см. рис.4). Для расчета линейного расхода массы используются полученные формулы (9) и подинтегральное выражение из (10).

Литература

1. Гайсин А.Ф., Абдуллин И.Ш., Фахрутдинова И.Т., Гайсин А.Ф., Гайсин Ф.М. Об особенностях получения наночастиц с использованием электрического разряда с жидкими электродами. Вестник Казанского технологического университета. 2012. Т. 15. № 22. С. 64-65.

2. Гайсин А.Ф., Абдуллин И.Ш., Гайсин Ф.М. Электроплазменная обработка стоматологических зубных мостов и коронок.Вестник Казанского технологического университета. 2012. Т. 15. № 21. С. 47-48.

3. Давыдов А.Д., Козак Е. Высокоскоростное электрохимическое формообразование. М. Наука, 1990.

4. Седыкин Ф.В. Размерная электрохимическая обработка деталей машин. М. Машиностроение, 1976.

5. Гуревич М.И. Теория струй идеальной жидкости. М. Наука, 1979.

6. Котляр Л.М., Миназетдинов Н.М. Определение формы анода с учётом свойств электролита в задачах электрохимической размерной обработки металлов. ПМТФ., 2003, Т.44 №3. С.179-184.

7. Котляр Л. М., Мельничук Е.Д., Шакирова Г.Ю., Шакиров Ю.И. Математическая модель диспергирования металла в электрическом разряде с жидким электролитом. Монография под редакцией Та<1еш7 7аЪого%г8к1, Ооко%г Ш11кр., 2012г.

8. Лаврентьев М. А., Шабат Б.В. Методы теории функций комплексного переменного. М. Наука, 1987.

9. Янке Е., Эмде Ф., Лёш Ф. Специальные функции. М. Наука, 1977.

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

© Л. М. Котляр - д-р физ.-мат. наук, профессор кафедры математики, Набережночелнинский институт Казанского (Приволжского) федерального университета, kotlyar-leonid@mail.ru; Г. Ю. Шакирова - старший преподаватель кафедры сервиса транспортных систем, Набережночелнинский институт Казанского (Приволжского) федерального университета, guzya81@mail.ru; И. Х. Исрафилов - д-р техн. наук, профессор кафедры высокоэнергетических процессов и агрегатов, Набережночелнинский институт Казанского (Приволжского) федерального университета, irmaris@yandex.ru; М. Ф. Шаехов -д-р техн. наук, профессор кафедры плазмохимических и нанотехнологий высокомолекулярных материалов, КНИТУ, shaechov@kstu.ru.

© L. M. Kotlyar - Doctor of Physical and Mathematical Sciences, Professor of the Department of Mathematics, Institute of Naberezhnye Chelny Kazan (Volga Region) Federal University, kotlyar-leonid@mail.ru; G. Yu. Shakirov - Senior Lecturer of the department of Service of transport systems, Institute of Naberezhnye Chelny Kazan (Volga Region) Federal University; guzya81@mail.ru; 1 H. Israfilov - Doctor of Technical Sciences, Professor of the Department of high-energy processes and units, Institute of Naberezhnye Chelny Kazan (Volga Region) Federal University, irmaris@yandex.ru; M. F. Shaehov - Doctor of Technical Sciences, Professor of the department plazmohimicheskim macromolecular materials and nanotechnology, KNITU, shaechov@kstu.ru.