Математическая модель задачи о профилировании инструмента при электрохимической обработке Текст научной статьи по специальности «Физика»

УПРАВЛЕНИЕ, ВЫЧИСЛИТЕЛЬНАЯ ТЕХНИКА И ИНФОРМАТИКА

УДК 621.794

А. И. Воронкова, Л. М. Котляр

МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ЗАДАЧИ О ПРОФИЛИРОВАНИИ ИНСТРУМЕНТА ПРИ ЭЛЕКТРОХИМИЧЕСКОЙ ОБРАБОТКЕ

В работе построена математическая модель плоской задачи электрохимической обработки (ЭХО) металла с учетом влияния гидродинамики течения электролита в межэлектродном канале (МЭК). Получено решение задачи о профилировании участка безкавитационного катода в областях с большим градиентом давления. Электрохимия; анод; катод; потенциал; кавитация; градиент; гидродинамика; электродинамика

Электрохимическая обработка металлов занимает важное место в современном машиностроении при производстве изделий, обработка которых традиционными механическими способами затруднена или невозможна. Одной из основных проблем в технологии ЭХО является точность анодного формообразования, повышение которой, в частности, связано с развитием математических методов расчета поверхности детали. Заметное влияние на точность формообразования может оказать переменный выход по току, высокая скорость прокачки электролита, начальная его загазованность, выделение газообразных продуктов реакции [1, 2, 3].

1. СОСТОЯНИЕ ВОПРОСА

Одной из актуальных задач теории ЭХО является изучение особенностей гидродинамики потока электролита и учета этих особенностей при проектировании формы катода. Высокая скорость течения электролита в МЭК, начальная его загазованность, выделение газообразных продуктов реакции приводят к образованию кавитационных полостей на границе катода в окрестности областей с большими градиентами давления. Каверны частично экранируют катод, что существенно влияет на форму границы обрабатываемой поверхности. Для устранения кавитационных полостей необходимо проектировать форму катода.

В данной работе решена задача о проектировании участка катода в области с большим градиентом давления при ЭХО металла катодами прямоугольной формы.

2. ПОСТАНОВКА ЗАДАЧИ

Рассматривается плоская задача формообразования анодной границы Ж при стационарном режиме ЭХО катодом ЛСОБ с линейными гра-

Контактная информация: (8552)39-66-74

ницами АС, ОБ и криволинейным участком СО (на рис. 1, в плоскости г = х + іу, показана правая половина схемы МЭК, А8 - ось симметрии). Течение электролита в МЭК направлено от точки А к точке Б. Начало координат выбрано в точке О. Ось абсцисс перпендикулярна направлению подачи катода и направлена по грани ВБ. Ширина МЭК в окрестности точки Б равна Н, где Н - характерная длина, равная ширине МЭК в случае стационарной ЭХО при плоских параллельных границах электродов и выходе по току, равном единице.

х

Б

Рис. 1. Схема МЭК

Для построения безкавитационного катода, форму криволинейного участка СО, получаемого в результате сглаживания острой кромки катода, будем строить таким образом, чтобы значение скорости течения электролита на этом участке было постоянным, равным У0. Скорость течения электролита на выходе из МЭК равна У1. Расход электролита равен q (# = У1Н).

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

что комплексные потенциалы W и WE гидродинамического и электродинамического полей связаны соотношением

dW dW dWr,

dWE dz

на анодной границе ББ dj

V = (а + Ьсоз(о)) ,

dФ е

(1)

где V - модуль скорости течения электролита на анодной границе, 0 - угол наклона вектора скорости к оси абсцисс, фЕ, ф - потенциалы электростатического и гидродинамического полей соответственно, a и Ь постоянные, зависящие от свойства электролита.

Требуется определить форму поверхности анода и криволинейного участка катода.

3. МЕТОДИКА ИССЛЕДОВАНИЯ

Для решения задачи введем вспомогательное комплексное переменное u = 2 + і'п, изменяющееся в области Ви - прямоугольнике со сторонами п/2, пт/4 (т = і|т|) (рис. 2), и будем искать функцию Z(u), конформно отображающую область Ви на область течения в плоскости г = х + іу с соответствием точек, указанном на рис. 1, 2.

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

С(и ) =1п

^у ёг(и1_> ^ 0 dW(u)

Область изменения W(u)

(2)

представлена на рис. 3. Функция dW / du мнима на вертикальных сторонах, вещественна на горизонтальных сторонах прямоугольника Ви и имеет про-стые полюса в точках А(и = іа), Б(и = п/2) и простые нули в точках С(и = пт/4), В(и = п/2 + пт/4), Б(и = 0). Согласно принципу симметрии, функцию dW / du можно аналитически продолжить на всю плоскость.

С

Б

Б

А Б

Рис. 3. Область изменения W По известным особенностям найдем [5]:

dW

= Ы-

Д (2и)#1 (и)

du Ф2 (и)Ф3 (и) ^1 (и - іа)^1 (и + іа)

Ф4(и)

-х

х— з

Ф4 (и - /а)Ф4 (и + /а)

где Фк (и), (к = 1,4) - тета-функции для пе-

риодов р и К1 [6]. Постоянная N определяется из условия, что расход жидкости в струе равен д. Определяя вычет функции Ж(и) в точке Б, получим

N = ф2(/а)ф2(/а)ф4(0).

р

Рассмотрим функцию Жуковского для течения электролита

С(и ) =1п

(

У

dZ(u)

Ш(и)

\

:1п| УУ0 I + і0(и).

На полигональных границах катода мнимая часть х(и) кусочно постоянна, а на криволинейном участке реальная часть %(и) равна нулю. На анодной границе ББ выполняется условие стационарности формообразования (1). Таким образом, для функции х(и) имеем следующие граничные условия:

1ш%(и) = -к/ 2 (и = і^),

Іш С(и) = 0 (и = к/ 2 + іГ|),

Яе с(и) = 0 (и = Х + кт/ 4),

Яе %(и) = 1п

dФ£

У

(и = X).

dф а + Ь со8 0(и)

Для определения функции dфE / dф при и = 2, отобразим область Ви (рис. 2) на верхнюю полуплоскость Вт (рис. 4) с соответствием точек, указанном на рисунках, преобразованием

ю(и) = sn(2K(2и/к -1/2), к),

К = уЛ(0, т), к =

( А2(0, т) ^2 V ^3(0, т

где т(и) - функция Якоби [6].

q

Рис. 4. Область изменения переменой m

Рассмотрим комплексный потенциал WE(u) = фЕ(и) + /уЕ(и) электростатического поля. Область изменения ^^Е(и) представлена на рис. 5.

A C D B

►

А Б

Рис. 5. Область изменения WE

Производные функций, отображающих область Вт (рис.4) на области изменения функций W(u), ^Е(и) (рис. 3, 4) с соответствием точек, указанном на рисунках, имеют вид:

= д(1 - а)

dw к(1 -ю)(ю-а)’

= V 2(1 -а)

ёю к(ю-1)^(1 + ю)(ю- а) ’

где а = ю(іа).

Используя полученные формулы, найдем

dWv

1

2(ю- a)

dW q\ (1 + ю)(1 -a) На анодной границе dWE djE

dW dj

(-1 <ю< 1).

Функцию Жуковского х(и) для течения электролита будем искать в виде

С(и ) = Хо(и) + / (u),

где Х0(и) = Г0(и) + /00(и) - функция Жуковского для течения идеальной жидкости по схеме (рис. 1), когда на границе ББ 0о = 0, а на границе СО V = ¥0. Кроме того, х0(и) имеет в области Ои те же особенности, что и функция х(и), а Ди) -аналитическая в Ои и непрерывная в Ои функция.

Следуя работе [5], построим функцию х0(и) и уравнение для определения Ди).

Функция dx0 / du вещественна при и = £,, мнима при и = /'п, и = п/2 + /'п, и = 2, + пт/4

и имеет в точке Б(и = 0) простой полюс с вычетом равным 12 . Функцию ё%0 / ёи можно продолжить на всю плоскость согласно принципу симметрии. С учетом граничных условий, эллиптическую функцию ё%0 / ёи с периодами п и пт представим в виде линейной комбинации логарифмических производных тета-функций и затем интегрированием найдем

( (и ) ^

%0(и) = 1п '

J1(u)

+ A0 u + B0

Используя граничные условия, найдем, что

Л0 = 0, Б0 = °.

Сравнивая граничные условия для функций х(и) и х0(и), получим граничные условия для неизвестной функции Ди).

1т /(и) = 0 (и = гл, и = р2 + /п), (3)

Яе / (и) = 0 (и = Х + рт/ 4), (4)

(Е V ^

Re %(u) = ln

Re %0

(5)

dj a + b cos 9(u)

(u =X).

Для определения функции f(u) отобразим предварительно область Du на полукольцо Dt (рис. 6) с помощью функции

t = exp(2i (u - pt/4)). (6)

Рис. 6. Область изменения переменной /

Учитывая граничное условие (3), функцию Ди) можно аналитически продолжить на все кольцо и представить в виде ряда Лорана с вещественными коэффициентами

/ (0 = Е сГ,

п=-¥

или, с учетом (6),

/(и) = Есп ехр(2п/(и - га/4)).

П=-¥

Используя граничное условие (4), найдем, что с0 = 0 и сп = - с-п. Окончательно для функции Ди) получим следующее представление

А(и) = 2іЕсп 8Іп(2п(и -кт/4).

п=1

При и = 2 функция Аи) имеет вид

(

А (X) = 2 Е Сп

С08(2пХ)5Й

(

+ і • $іп(2п%)ск

2

V У

+

(7)

кт п

V 2 ,

V /у

Граничное условие (5) с учетом (7) запишется в виде

2 • Е Сп С08(2пХ)5Й

1

ёф Е У0

п=1

= 1п

r0, Г0 = ЯеС0.

ёф а + Ь со8(0)

Отсюда получим бесконечную систему

уравнений для определения Сп

2

1п

V V

фЕ

У

л

л

- г

ёф а + Ь со8(0)

М = к • 5Й(к|т|п/2), п = 1, ¥ .

Безразмерные координаты точек анодной границы и границы криволинейного участка катода определяются из (2) по формуле

Z (и) 1 и'ёШ /ч,

—1— =-------[----ехртёи

Н У0 Н 0 ёи ^

при и = 2 и и = 2 + пт/4, соответственно.

4. ОБСУЖДЕНИЕ РЕЗУЛЬТАТОВ

На рис. 7 представлены результаты расчетов формы границы сглаженного участка катода и границы анода при ЭХО с постоянным выходом по току и одном отношении У0 / У = 2 для различных значений Н1 / Н, где Н1 - половина ширины МЭК в окрестности точки А.

Рис. 7. Результаты расчета анодных границ при стационарной ЭХО:

1) Н1 / Я = 3; 2) Н1 / Я = 2; 3) Н1 / Я = 1,5

Рис. 8. Результаты расчета анодных границ при стационарной ЭХО: 1-3 ЭХО безкавитационными катодами 1) Н1 / Н = 3; 2) Н1 / Н = 2; 3) Н1 / Н = 1,5;

4-5 ЭХО прямоугольными катодами 4) Н1 / Н = 3; 5) Н1 / Н = 2; 6) Н1 / Н = 1,5

Приведенные результаты и сравнение их с результатами расчета анодных границ, полученными с учетом влияния присоединенной кавитации (рис. 8), показывают, что кавитация достаточно сильно влияет на форму обрабатываемой детали. Проектирование безкавитацион-ного катода позволяет получить поверхность детали с меньшими затратами на последующую механическую обработку. Наличие каверны на катоде может привести к неустойчивости течения в МЭК, что также влияет на качество обработки детали.

На рис. 9 представлены результаты расчетов формы границы сглаженного участка катода и границы анода при ЭХО с постоянным выходом по току и одном отношении Н1 / Н = 2 для различных значений У0 / У1.

Рис. 9. Результаты расчета анодных границ при стационарной ЭХО:

1) = 1,5; 2) ^ = 2;3) ^ = 2,5

У1

У1

У1

1

-1

0

Рис. 10. Результаты расчета анодных границ при стационарной ЭХО: 1) а = 1/3; Ь = 2/3; 2) а = 0; Ь = 1; 3) а = -1/5; Ь = 6/5

п=1

2

С

На рис. 10 представлено влияние переменного выхода по току на электрохимическое формообразование.

5. ПРИЛОЖЕНИЕ РЕЗУЛЬТАТОВ

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

ВЫВОДЫ

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

СПИСОК ЛИТЕРАТУРЫ

1. Зайцев А. Н., Житников В. П. Высокоскоростное анодное растворение в условиях нестацио-нарности электродных потенциалов. Уфа: Гилем, 2005. 220 с.

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

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

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

5. Киселев О. М., Котляр Л. М. Решение нелинейных задач теории струйных течений тяжелой жидкости. Казань: КГУ, 1978. 154 с.

6. Уитекер Э. Т., Ватсон Дж. Н. Курс современного анализа. Ч. 2. М.: Наука, 1962. 686 с.

ОБ АВТОРАХ

Котляр Леонид Михайлович, зав. каф. высш. матем. Камск. гос. инж.-экон. академии (ИНЭКА). Д-р физ.-мат. наук (Казанск. гос. ун-т, 1974), проф. Иссл. в обл. гидродинамики, газовой динамики.

Воронкова Анна Ивановна, доц. той же каф. Канд. физ.-мат. наук (Казанск. гос. ун-т, 1997). Иссл. в обл. гидродинамики.