Численное моделирование двумерных электростатических полей при электрохимическом полировании Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

УДК658.512.24

Л. А. Бабкина, И. Я. Шестаков, А. С. Квасов

ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ ДВУМЕРНЫХ ЭЛЕКТРОСТАТИЧЕСКИХ ПОЛЕЙ ПРИ ЭЛЕКТРОХИМИЧЕСКОМ ПОЛИРОВАНИИ

Рассматривается задача моделирования двумерного электростатического поля при электрохимическом полировании на основе биполярного электрода. Проводится сравнительный анализ результатов численного моделирования в программном комплексе COSMOS/Mи пакете математических вычислений Maple.

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

Электрохимическое полирование (ЭХП) широко применяется в производстве деталей летательных аппаратов и занимает ведущее место в доводочных операциях [1].

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

При электрохимическом полировании проблема равномерности распределения тока по поверхности деталей является первостепенной. От характера распределения тока зависит качество поверхности, в том числе ее шероховатость [3; 4].

В качестве объекта исследования рассматривается электрическое поле, возникающее при доводке поверхности детали электрохимическим полированием на основе биполярного электрода. Задачами исследования являются создание численной двумерной модели распределения электростатического поля в межэлектродном зазоре (МЭЗ) и разработка алгоритма реализации численного решения в среде интегрированных пакетов конечно-элементных расчетов.

Исследованию электрических полей МЭЗ посвящен ряд работ. Одни исследователи выявляли распределение потенциала экспериментальным путем, т. е. непосредственным измерением или моделированием, другие пытались выполнить аналитическое решение. В этих работах описаны случаи с практически неизменными границами поля [4].

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

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

тизма является переход от системы уравнений Максвелла к уравнениям второго порядка относительно напряженности электрического поля или напряженности магнитного поля [6].

(1)

(2)

(3)

(4)

Рис. 1. Принципиальная схема электрохимического полирования

Уравнения Максвелла, используемые в конечно-элементном расчете, имеют следующий вид [6]:

Ух н = з+8Б/ 5/,

Ух Е = - дБ/ 8/,

У-Б = 0,

У-Б = р ,

где Н и Е - напряженность магнитного и электрического полей соответственно; Б, Б - плотности магнитного и электрического потоков; 3 - плотность тока в проводнике; р - плотность электрического заряда; У - оператор градиента.

Составляющие уравнения имеют следующий вид:

Б = ц( Б)( Н + Нс), (5)

3 = о - Е, (6)

Б = е - Е, (7)

где ст - электропроводность; ц - магнитная проницаемость; е - диэлектрическая постоянная; Нс - коэрцитив-ность постоянного магнита.

При двумерном расчете плотность потока Б выражена как функция векторного потенциала А:

В = У- А , (8)

где А удовлетворяет условию уникальности

У- А = 0 (9)

(14)

(15)

Подставляя уравнения (8) и (5) в уравнения (1) и (2) и сокращая частную производную по времени dA/dt, получают следующее [6]:

Vx(v (B)Vx A) = J + Vx Hc , (10)

VxE = Vx(dA/dt) , (11)

где v - удельное магнитное сопротивление, V = (1/ д).

Уравнение (11) предполагает

E = dA/dt -Vcp , (12)

где j - сокращенный электрический скалярный потенциал.

Объединение этого уравнения с составляющим уравнением (1-6) приводит к следующему:

Je =-с (dA/ dt + Vcp). (13)

Полная плотность тока в уравнении (10) может состоять из токов источника J и вихревых токов J, поэтому уравнение (10) можно записать в следующем виде [6]:

Vx ( v ( B ) Vx A) + с dA dt +

+gV9 = JS +Vx HC.

Для двухмерных моделей в плоскости xy единственные ненулевые компоненты A и Vj - это z-компоненты, являющиеся функциями только от х иy, и не изменяются по направлению z. Таким образом, это уравнение принимает следующую скалярную форму:

d/dx ( v dA dx ) + d/dy (vdA/ dy ) -

-с dA dt = - JS - dHCy/ dx + dHCx/ dy.

В случае осесимметричности, если принять z за азимутальное направление, единственным ненулевым компонентом A будет азимутальный. В этом случае уравнение (15) принимает следующую скалярную форму: d/ dx (vdA/ dx + vA/x ) + d/ dy (vdA/dy ) -

-с dA dt = - JS + dHCy/ dx - dHCx/ dy.

В осесимметричном случае x и y соответствуют радиальному и осевому компонентам цилиндрических координат [6].

Для решения поставленной задачи использована программа GEOSTAR и модуль ESTAR для моделирования электрических полей в двух- и трехэлектродных ячейках. GEOSTAR выполняет функции пре- и постпроцессора для системы COSMOS/M. Геометрические возможности GEOSTAR базируются на методе смешанных граничных представлений и параметрических кубических уравнениях. Теория, используемая в ESTAR, основана на применении теории потенциальных функций к уравнениям Максвелла. Двухмерный и осесимметричный расчеты в ESTAR основаны на уравнениях (15) и (16), уравнения (1)-(14) выступают в качестве исходных данных [6].

Моделируется реальная схема электрохимического полирования цилиндрической трубки, установленной симметрично относительно двух плоских электродов (рис. 1). После цикла обработки полярность подключения электродов меняют на противоположную, либо поворачивают деталь на 180 ° относительно ее оси. Особенностью схемы является то, что обрабатываемая деталь играет роль биполярного электрода, т. е. на ее поверхности одновременно протекают анодные и катодные реакции. Ток к обрабатываемой детали подводится через электролит сразу по всей ее поверхности. Такое решение токоподво-

(16)

да исклю-чает структурные изменения металла и прижо-ги, исчезает погрешность обработки [7].

Для понимания явлений, протекающих на биполярном электроде, в работе [2] рассмотрена причина, вызывающая одновременное протекание катодной и анодной реакций на ее поверхности. Физическая природа этого явления связана с появлением поверхностных зарядов разных знаков при протекании тока через границу неоднородных сред (рис. 2).

Рис. 2. Стационарное электрическое поле в проводнике, состоящем из участков с разным удельным сопротивлением(р2 > р1 ) [2]

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

ДЕ = Е2 - Е1 — Е1 І -1

І Рі

где Е - нормальная составляющая электрического поля; р12 - удельное сопротивление среды.

Поверхностная плотность появляющихся зарядов равна

о = є 0 Ё І ^ -1

I р1

где є - диэлектрическая проницаемость.

В случае р2 >> р1 плотность зарядов

О — єо— єор2./2 .

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

Исследованию распределения тока по биполярному электроду в технической электрохимии посвящен ряд работ [2; 4; 8]. Основное внимание в этих работах уделяется изучению распределения тока по биполярному электроду.

Электростатическое поле исследовалось на модели, исходная схема которой показана на рис. 3. Здесь область, ограниченная точками (1, 2, 3, 4) соответствует сечению А-А (рис. 1) и представляет собой электропроводящую среду (в ЭХП - электролит); линии 12, 13, 14, 15, 7, 8, 5, 6 -

поверхности электродов; линии 9-10-11, 16-17-18 - поверхности обрабатываемой детали. По линиям 12-13, 14-15, 7-8, 5-6 моделируемая область имеет приложенные нагрузки в виде электрического напряжения, граничный потенциал вдоль этих линий будет равен: ф12-13 = +16 В; ф7-8 = 0 В . Для конечно-элементного разбиения расчетной области модели выбран плоский треугольный элемент.

I Із II

.1 Ні

18 5

З 5

Рис. 3. Исходная схема для расчета в программе COSMOS/M

В процессе моделирования двумерного электростатического поля при электрохимическом полировании рассмотрены случаи, когда ширина плоского электрода (L) в три и более раз больше наружного диаметра детали (D), сопротивление единичного поперечного сечения электролита (R) равно или больше поляризационного сопротивления (r).

Сопротивление поперечного сечения электролита

— = х • b, где с - удельная электропроводность (таблич-R

ные значения [7]); b - межэлектродное расстояние.

Для рассматриваемых случаев обработка ведется в нейтральном водном растворе солей при плотности тока i = 0,2...10 А/см2; для достижения необходимой плотности тока требуется напряжение U до 16 В; материал электродов должен быть устойчив при подаче анодного потенциала; материал образца - медь; D = 47 мм; L = 150 мм; b = 140 мм; X = 24,7з 1/ Ом- м для 5М NaCl [7]; анодное поляризационное сопротивление r » 0,014 Ом; ка-

тодное поляризационное сопротивление г » 0,4 Ом (по анодным и катодным поляризационным кривым [7]).

1 — х.ь — 24,73 • 0,014 — 0,346 (Я » 3 Ом);

Я

г = 2г + 2г = 0,828 Ом, т. е. Я > г.

а к 7 7

Также рассмотрены случаи, когда ширина плоского электрода меньше наружного диаметра детали Ь/Ю < 3, сопротивление единичного поперечного сечения электролита равно или больше поляризационного.

В результате расчетов установлено распределение напряженности электрического поля в исследуемой области, на основании этого получено распределение плотности тока по обрабатываемой поверхности (рис. 4, 5).

Значения потенциалов в узлах конечно-элементной сетки (рис. 5, б) по обрабатываемой поверхности (рис. 5, а) приведены в таблице.

Неравномерность распределения потенциала по поверхности обрабатываемой детали составляет » 0,78 %.

Наиболее равномерное распределение по всей поверхности обрабатываемой детали, обращенной к катоду, потенциалов и силовых линий электрического тока происходит, когда ширина плоского электрода в три и более раз больше наружного диаметра детали Ь/Ю > 3 (рис. 4, 5).

Рис. 4. Картина распределения потенциалов при ЬЮ > 3 (Е8ТЛК)

В рассматриваемых случаях, когда ширина плоского электрода в три и более раз меньше наружного диаметра детали Ь/Ю < 3, возникают краевые эффекты и происхо-

Таблица

Номер Значение

конечно-элементного узла потенциала, j

СГСЭ СИ

229 0,026 46 7,938

230 0,026 46 7,938

231 0,026 47 7,941

232 0,026 47 7,941

233 0,026 48 7,944

234 0,026 49 7,947

235 0,026 51 7,952

236 0,026 52 7,956

237 0,026 54 7,962

238 0,026 56 7,968

239 0,026 58 7,974

240 0,026 60 7,980

241 0,026 62 7,986

242 0,026 65 7,995

243 0,026 67 8,000

дит неравномерное распределение потенциалов по поверхности обрабатываемой детали.

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

В работе [8], посвященной разработке и исследованию аналитических и численных моделей электрических полей, возникающих в системах на основе биполярного электрода, получены следующие разностные схемы для вычисления распределения электрического поля:

(е0 + ег )иг + (е0 + е1 )и1 +

+(е 0 + е, )и, + (е 0 + е ь )иь

и (х, у) =

и (х, у) =

(е 0 + е г ) + (е 0 + е1 ) +

+(е 0 + е1) + (е 0 + еь )

(00 + 0г )иг + (00 + ©1 )и1 +

+(о0 + О )Ц, + (о0 + 0ь )иь (00 + 0Г) + (©0 + О,) + ’

+(о0 + О ) + (о0 + 0ь )

(17)

(18)

Пг = и (х + Кх, у);

П1 = и (х - Нх, у)

где

и, = и (х у+К); иь = и (х, у - К); и 0 = и (х, у)

е0 = е(х, у); ег = е(х + Нх, у); е, = е(х - Кх, у)

е, =е(х, у+К); еь = е(х у - К); ©0 = о(х, у)

о г = о(х + Нх, у); о, = о(х - Нх, у); о, = о(х, у + Ку) оь = о(х, у - Ку); К - размер ячейки сетки.

При этом в ячейках сетки, соответствующих электродам, устанавливается постоянное напряжение, равное

напряжению моделируемого электрода; в ячейках, граничащих с проводниками, потенциал вычисляется по формуле (18), а в остальных ячейках поле вычисляется по формуле (17). Можно задавать различные конфигурации электродов и параметры электролита, решать задачу оптимизации параметров распределения поля [8].

Исходя из этих условий, нами реализована предложенная в [8] математическая модель двумерного электростатического поля применительно к электрохимическому полированию в пакете математических вычислений Мар1е. Для реализации численного решения в Мар1е разработан алгоритм расчета, заключающийся в следующем: задается исследуемая область; исследуемая область разбивается на конечно-разностную сетку; задаются численные значения для каждого моделируемого объекта области (определяется проводящая среда - задается электропроводность среды, 1/Ом-м; потенциал электродов, В); составляется матрица размером п х т; выполняется расчет, используя разностные схемы (17) и (18).

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

Картины распределения потенциалов и силовых линий электрического тока в рассматриваемом случае, когда ширина плоского электрода в три и более раз больше наружного диаметра детали Ь/Б > 3, распределяются почти равномерно по всей поверхности обрабатываемой

с*! г.чи зт пт* хгы гча ггм г*и

а б

Рис. 5. Силовые линии электрического тока при Ь/Б > 3(Е8ТЛБ.): а - увеличенный фрагмент; б - номера узлов сетки

детали, обращенной к катоду (рис. 6; 7, а). Неравномер- геометрические параметры Ь и Б, а также физико-хими-ность распределения потенциала по поверхности обра- ческие свойства электролита.

батываемой детали составляет и 1,28 % (рис. 7, б).

■» « я

1

1 J- 1?

1 Г и

J 12 13 п-1:

>:■ I:

Г !■>

S и

п i :

ii “J-Г

J-i J-J

Рис. 6. Распределение электрических потенциалов при L/D > 3 (Maple)

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

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

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

Рис. 8. Распределение электрического потенциала по обрабатываемой поверхности: 1 - кривая распределения электрического потенциала по обрабатываемой поверхности при L/D > 3, Maple; 2 - кривая распределения электрического потенциала по обрабатываемой поверхности при L/D > 3, ESTAR

Таким образом, в работе реализована предложенная в [9] математическая модель двумерного электростатического поля применительно к электрохимическому полированию. Численное моделирование электростатических полей в МЭЗ выполнено в среде интегрированного пакета конечно-элементных расчетов COSMOS/M (модуль ESTAR) и конечно-разностный расчет в Maple. Сравнительный анализ результатов численного моделирования в Maple и COSMOS/M показал расхождение » 2 %.

I ? t > .т

I t f * .• У ^ * -г ^ .

; 1 f f .• у .f , „ „ .

f Г f /V Л J

| I Г i! У Г У .-.p.,

" м , 14 ^ '■ '* V ’■ I i ' ’ ’■ '-A V k r

\ ki 1 ■ ' * * 14 \ , S

•"ЧЪЧ j f L 1

r i- *■ -Г -#■ -

■. v-!

„-r___________'

Г .' .* _> j- i._______________________________________________________________________________

/ i У I

-i; #4 * i / -1

■ ~ я . r i - - V - t

”--v - \...........", - - ‘7--------------

...... - .—,----------------'---------------------

;v <3

1-. i ■■ i[i ■ mil i :i £TXJJ } .. Ai'l if

дткйвл ft ii w;

S, 2-1 S.'X'O LID?

S.2S 7.^-П 1. fl9u

r 27 7.0^0 1.22*

7SMZ L.-^i

*,29 7.911 L.®

I,» 7.310 (I.S87

2,3 I 7.SS7 fj.RLlI

1*31 7.^-7 U. Г2

a 6

Рис. 7. Силовые линии электрического тока при L/D > 3 (Maple): а - увеличенный фрагмент (черным цветом выделены координаты ячеек, описывающих область поверхности детали); б - изменение распределения потенциалов и плотности тока по обрабатываемой поверхности

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

Результаты численного моделирования подтверждают технологические возможности метода ЭХП с БПЭ по обработке полых цилиндрических деталей.

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

Библиографический список

1. Воробей, В. В. Теоретические основы проектирования технологических процессов ракетных двигателей. Технология производства жидкостных ракетных двигателей / В. В. Воробей, В. Е. Логинов. М. : Дрофа, 2007.

2. Мороз, И. И. Биполярный метод электрохимической обработки и некоторые его технологические возможности // И. И. Мороз, В. Ф. Орлов, Б. И. Чугунов // Электронная обработка материалов. 1982. N° 6. С. 19-23.

3. Экслер, Л. И. Классификация параметров шероховатости / Л. И. Экслер // Технологические методы повышения качества поверхности деталей машин : сб. ст. Л. : Изд-во Ленингр. ун-та, 1978. С. 140-147.

4. Вдовенко, В. Г Эффективность электрохимической обработки деталей : моногр. / В. Г Вдовенко. Красноярск : Изд-во Краснояр. гос. ун -та, 1991.

5. Monk, P. A finite element method for approximating the time-harmonic Maxwell equations / P. Monk // Numer. Math. 1992. Vol. 63. P. 243-261.

6. COSMOS Advanced Modules. Part 2 ESTAR / Low Frequency Electromagnetic Analysis. 1996.

7. Пат. Российская Федерация, 7С25F3/16. Способ электрохимического полирования / Шестаков И. Я., Бабкина Л. А. № 2229543 ; заявл. 15.07.2002 ; опубл. 27.05.2004, Бюл. № 15. Приоритет от 15.07.2002.

8. Мурашев, Д. А. Математическое моделирование электрических полей в системах с биполярным электродом : автореф. дисс. ... канд. физ.-мат. наук / Д. А. Мурашов. Саратов : Саратов. гос. техн. ун-т, 2006.

9. Шестаков, И. Я. О возможностях электрохимического полирования / И. Я. Шестаков, Л. А. Бабкина, А. Н. Жмурко // Решетневские чтения : материалы VII Всерос. научн.-практ. конф. ; Сиб. гос. аэрокосм. ун -т. Красноярск, 2003.

L. A. Babkina, I. Y а. Shestakov, A. S. Kvasov

NUMERICAL MODELING BY TWO-DIMENSIONAL ELECTROSTATIC FIELDS AT ELECTROCHEMICAL POLISHING

Modeling problem of the two-dimensional electrostatic field at electrochemical polish is considered. Comparative analysis of results of numerical modeling is executed in the programs COSMOS/M and Maple.

Keywords: numerical modeling, electrostatic field, electrochemical polish.

УДК 519.876

И. М. Митасов, А. Н. Завьялкин

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

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

Ключевые слова: метод, факторы, модель, выборка, проектирование.

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

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

Необходимость решения систем уравнений с квадратной матрицей является следствием необходимого условия минимума целевой функции метода НК. Фундаментальная интерпретация метода НК состоит в том, что оп-