КРИТЕРИАЛЬНЫЙ МЕТОД РАСЧЕТА РАСПРЕДЕЛЕНИЯ ТОЛЩИНЫ ПОКРЫТИЯ НА КАТОДЕ ДЛЯ ЭЛЕКТРОЛИТА ХРОМИРОВАНИЯ С НЕМОНОТОННОЙ КРИВОЙ КАТОДНОЙ ПОЛЯРИЗАЦИИ
Елизаров А.М., Литовка Ю.В.
Тамбовский государственный технический университет
392000 г.Тамбов, ул.Советская, 106, тел.(0752)722601,
Важной особенностью хромирования является немонотонность N-образной поляризационной кривой стандартного хромировочного электролита, что сказывается на равномерности распределения покрытия и осложняет расчет математической модели гальванического процесса.
Математическая модель распределения электрического поля в электролитической ячейке включает уравнение Фарадея, уравнение Ома в дифференциальной форме и уравнение Лапласа с нелинейными краевыми условиями, учитывающими поляризационные явления на электродах.
Толщина покрытия в каждой точке катода рассчитывается по формуле, полученной из закона Фарадея:
(1)
где Э - электрохимический эквивалент хрома Э = 8,98×10-8 [кг/(А с)]; r- плотность вещества (для хрома r = 7000 [кг/м3]); h- катодный выход по току; - плотность тока в точке с координатами xk, yk, zk в момент времени τ; - концентрации хромового ангидрида и серной кислоты соответственно; Tp – время нанесения покрытия.
Катодный выход по току является функцией температуры t и плотности тока ik. Выражение, аппроксимирующее экспериментальные данные, имеет вид:
(2)
Исходя из закона Ома в дифференциальной форме, плотность тока на катоде будет рассчитываться по формуле:
, (3)
где c - удельная электропроводность электролита; n – нормаль к поверхности катода; j- потенциал электрического поля электролизера.
Для нахождения распределения потенциала j в ванне используется дифференциальное уравнение Лапласа
(4)
со следующими краевыми условиями:
(5)
где Sn – площадь поверхности изолятора; Sа– площадь поверхности анода; F1(ia)– функция анодной плотности тока, учитывающая поляризацию на аноде; F2(ik)– функция катодной плотности тока, учитывающая поляризацию на катоде; UА – напряжение на аноде.
Так как в процессе хромирования используется нерастворимый свинцовый анод, то его поляризация не учитывается, поэтому F1(ia) = 0.
Функция катодной поляризации определяется кусочно-заданной N-образной функцией:
F2(ik) = (6)
Расчет распределение электрического поля в гальванической ванне осуществляется в соответствии с сеточным методом верхней релаксации, так как данный метод обеспечивает наибольшую скорость сходимости. Основная проблема при решении дифференциального уравнения Лапласа заключается в немонотонности катодной поляризационной кривой, вследствие чего в процессе решения задачи возникает необходимость выбора одной из трех ветвей N-образной кривой катодной поляризации.
Для решения данной проблемы нами предлагается метод, суть которого заключается в следующем.
Основная идея алгоритма заключается в сведении задачи поиска решения системы уравнений математической модели к задаче оптимизации.
Как видно из краевого условия на катоде (5), прикатодный потенциал принят равным нулю, следовательно, потенциал на катоде численно равен потенциалу поляризации. Тогда, согласно формулам (3, 5), итерационный процесс сойдется в том случае, когда значение плотности тока на катоде, исчисленное по закону Ома (3), будет близко с точностью ε к значению плотности тока, полученному из краевого условия на катоде (5). Таким образом, целевая функция будет иметь следующий вид:
, (7)
где m – количество узлов сетки, приходящихся на катод; Nk – номер узла сетки; - функция, обратная функции катодной поляризации (6).
Для получения очередного n+1-го приближения аргумента целевой функции используется метод итераций. При этом на каждой итерации необходимо выбрать одну из трех ветвей N-образной поляризационной кривой. С этой целью вводится n+0.5-я итерация, на которой ставится дополнительная задача оптимизации.
Найти оптимальное значение коэффициента K*, доставляющее минимум целевой функции следующего вида.
(8)
где K – некоторый коэффициент пропорциональности; - оптимальное смещение прикатодного потенциала в каждой точке поверхности катода, необходимое для выбора ветви поляризационной кривой.
Как правило, расхождение итерационного метода расчета распределения плотности тока происходит из-за резких перепадов значений потенциалов поляризации в соседних узлах сетки, приходящихся на катод, что характерно для немонотонной поляризационной кривой. Чтобы избежать таких ситуаций, оптимальное смещение прикатодного потенциала во всех узлах сетки, приходящихся на катод, будем искать, исходя из значений плотности тока в этих узлах и единого коэффициента пропорциональности K в соответствии с формулой:
(9)
где - плотность тока на катоде на n+0.5-ой итерации, определяющаяся по закону Ома в дифференциальной форме (3) после предварительного расчета по методу верхней релаксации; - прикатодный потенциал на n+0.5-итерации, определяеющийся из краевого условия на катоде (5) по формуле:
, (10)
где - вспомогательный потенциал (критерий выбора ветви поляризационной кривой), определяющийся по формуле:
(11)
Диапазон возможных значений коэффициента K выбирается таким образом, чтобы значения произведения были соизмеримы со значениями в формуле (11) и не допускали смещение потенциала в обратном направлении.
Для решения данной подзадачи можно использовать любой метод одномерной оптимизации, например, метод золотого сечения.
Далее по полученному распределению плотности тока и потенциалов на катоде происходит сравнение значений целевой функции (7) на n+1-ой и n-ой итерациях.
В случае, если значение целевой функции на очередной итерации хуже, чем на предыдущей, сканируются значения аргумента между значениями на предыдущей и на текущей итерациях пропорционально расстоянию между соответствующими значениями аргумента. Если лучшее значение найдено, то оно является очередным приближением к оптимуму. Иначе осуществляется попытка найти лучшее значение целевой функции, смещая значение аргумента относительно значения на предыдущей итерации в более электроотрицательную область. Если лучшее значение найдено, то оно является очередным приближением оптимуму, в противном случае значение на предыдущей итерации считается оптимальным.
Альтернативным критерием выхода считается достижение целевой функцией (7) значения меньшего заданной точности e.
Сравнение результатов работы данного алгоритма с экспериментальными данными показало, что среднеквадратичное относительное отклонение экспериментальных и расчетных значений толщины покрытия, рассчитанное по формуле , составило 23,03%, что соизмеримо с погрешностью измерений.