MATHEMATICAL MODELING AND OPTIMAL CONTROL OF REVERSE MODE OF ELECTROLYTE COATING APPLICATION IN MULTI-ANODE BATH
Journal: VESTNIK OF ASTRAKHAN STATE TECHNICAL UNIVERSITY. SERIES: MANAGEMENT, COMPUTER SCIENCE AND INFORMATICS
Volume 2015 № 2
, 2015
Rubrics: MANAGEMENT, MODELING, AUTOMATION
Authors:
1.
Tambov State Technical University
2.
Tambov State Technical University
3.
Tambov State Technical University
Type:
Article
Pages:
from 7
to 15
Status:
Published
Received:
05.04.2015
Accepted:
28.02.2020
Published:
25.04.2015
Subject area:
VAC 05.12.2013 Системы, сети и устройства телекоммуникаций
VAC 05.13.01 Системный анализ, управление и обработка информации (по отраслям)
VAC 05.13.06 Автоматизация и управление технологическими процессами и производствами (по отраслям)
VAC 05.13.10 Управление в социальных и экономических системах
VAC 05.13.18 Математическое моделирование, численные методы и комплексы программ
VAC 05.13.19 Методы и системы защиты информации, информационная безопасность
UDK 681.5
GRNTI 20.01 Общие вопросы информатики
GRNTI 28.01 Общие вопросы кибернетики
GRNTI 49.01 Общие вопросы связи
GRNTI 50.01 Общие вопросы автоматики и вычислительной техники
GRNTI 82.01 Общие вопросы организации и управления
VAC 05.13.01 Системный анализ, управление и обработка информации (по отраслям)
VAC 05.13.06 Автоматизация и управление технологическими процессами и производствами (по отраслям)
VAC 05.13.10 Управление в социальных и экономических системах
VAC 05.13.18 Математическое моделирование, численные методы и комплексы программ
VAC 05.13.19 Методы и системы защиты информации, информационная безопасность
UDK 681.5
GRNTI 20.01 Общие вопросы информатики
GRNTI 28.01 Общие вопросы кибернетики
GRNTI 49.01 Общие вопросы связи
GRNTI 50.01 Общие вопросы автоматики и вычислительной техники
GRNTI 82.01 Общие вопросы организации и управления
Language:
Russian
Keywords:
multiband bath, reverse mode, electrolyte processes, irregularity of coating, mathematic modeling, optimal control
Abstract (English):
The existing methods of reduction in the nonuniformity of the distribution of the thickness of electrolyte coating over the surface of the subject with the use of baths with various anodes are examined. The new way of reduction in the nonuniformity of coating, which is consisted in the combination of the mode of reversing current with turning off the corresponding anodic sections in multiband bath, is proposed. The block diagram of such galvanic bath is presented. The problem of the optimal control of the described process is stated. The algorithms of the solution of the system of equations of the mathematical model on the basis of grid and iteration methods are examined. The algorithm of the search for the solution of the problem of control on the basis of the method of complete enumeration is determined. The effectiveness of the proposed method in comparison with the classical mode of reversing in the bath with monolithic anode is demonstrated by the example of electrolyte coating application in 9 anodic galvanic nickel plating bath on the flat-like object. The effectiveness of reduction in the nonuniformity of electrolyte coating proves the expediency of applying the proposed method.
The existing methods of reduction in the nonuniformity of the distribution of the thickness of electrolyte coating over the surface of the subject with the use of baths with various anodes are examined. The new way of reduction in the nonuniformity of coating, which is consisted in the combination of the mode of reversing current with turning off the corresponding anodic sections in multiband bath, is proposed. The block diagram of such galvanic bath is presented. The problem of the optimal control of the described process is stated. The algorithms of the solution of the system of equations of the mathematical model on the basis of grid and iteration methods are examined. The algorithm of the search for the solution of the problem of control on the basis of the method of complete enumeration is determined. The effectiveness of the proposed method in comparison with the classical mode of reversing in the bath with monolithic anode is demonstrated by the example of electrolyte coating application in 9 anodic galvanic nickel plating bath on the flat-like object. The effectiveness of reduction in the nonuniformity of electrolyte coating proves the expediency of applying the proposed method.
Keywords:
multiband bath, reverse mode, electrolyte processes, irregularity of coating, mathematic modeling, optimal control
multiband bath, reverse mode, electrolyte processes, irregularity of coating, mathematic modeling, optimal control
Введение Получение равномерного гальванического покрытия является одной из основных задач автоматизации гальванических процессов. Для решения данной задачи применяются различные способы, одним из которых является применение ванн со многими анодами [1]. Данные ванны функционируют в режиме подвижных анодов, подключаются к нескольким источникам питания, реализуют циклическое переключение анодов. К недостаткам перечисленных режимов функционирования относят: сложность технической реализации, высокая стоимость необходимого оборудования, снижение производительности гальванической линии в целом. Наиболее перспективным способом снижения неравномерности распределения гальванического покрытия с использованием ванн со многими анодами является реализация реверсивного режима в многоанодной ванне, что позволит устранить существующие недостатки [2]. Однако эффективность применения такого режима невозможна без разработки алгоритмов моделирования и оптимального управления для предлагаемого процесса. Целью исследований являлась разработка алгоритмов для математического моделирования и оптимального управления реверсивным режимом нанесения гальванического покрытия в многоанодной ванне для снижения неравномерности наносимого покрытия. Гальванический процесс нанесения покрытия в многоанодной ванне, протекающий в реверсивном режиме В гальваническом процессе монолитный анод разбивается на систему секций одинакового размера, размещенных в одной плоскости на равном расстоянии друг от друга, стенок ванны и напротив изделия-катода. На рис. 1 представлена система из 9-анодных секций (размерностью M N = 3 3) и катод плоской формы. Данная система секций подключается к программируемому источнику питания типа «Пульсар Про», «Flex Kraft» через устройство, позволяющее отключать требуемые анодные секции для прямого и обратного режимов реверсирования тока для снижения неравномерности распределения гальванического покрытия, вследствие того, что время протекания прямого и обратного тока скорость растворения металла покрытия на более близких к анодным секциям поверхностях катода будет больше, чем на удаленных. Рис. 1. Схема 9-анодной гальванической ванны Постановка математической задачи оптимального управления рассматриваемым процессом звучит следующим образом: найти матрицы активных анодных секций для прямого и обратного режимов реверсивного гальванического процесса, а также длительности , периодов, доставляющие минимум критерию неравномерности распределения толщины наносимого покрытия по поверхности катода Sk: , (1) при этом должны выполняться следующие ограничения: , (2) (3) где - заданная толщина гальванического покрытия; - минимальная толщина гальванического покрытия; - максимальная толщина гальванического покрытия; Т - продолжительность гальванического процесса; Тmax - максимальная продолжительность гальванического процесса; x, y, z - координаты точки в пространстве гальванической ванны; τ - произвольный момент времени; υ - количество периодов включения прямого и обратного тока. Матрица активных анодных секций для режима «прямого» включения имеет следующий вид: . (4) Для «обратного» режима: , (5) где Алгоритм решения системы уравнений математической модели и задачи оптимального управления реверсивным режимом нанесения гальванического покрытия в многоанодной ванне В [3] предложена математическая модель гальванического процесса, система уравнений которой содержит закон Фарадея и закон Ома в дифференциальной форме, а также дифференциальное уравнение Лапласа в частных производных, описывающее распределение потенциала в пространстве многоанодной гальванической ванны. При этом конфигурация гальванической ванны задается геометрической моделью рецепторного типа. Отличительной особенностью разработанной математической модели является то, что в граничные условия для дифференциального уравнения Лапласа в частных производных добавлена проверка включения соответствующей m, n-й анодной секции для «прямого» и «обратного» режимов реверсирования. Данную математическую модель предлагается решать сеточным методом. Для сеточного представления решения исходной задачи в [3] приводится конечно-разностная аппроксимация 7-ми точечным шаблоном уравнения Лапласа с соответствующими краевыми условиями. При этом полученная система алгебраических уравнений решается комбинацией методов простой итерации (по распределению потенциала) и нижней релаксации (по плотности тока) [4]. Задача поиска оптимального управления гальваническим процессом (1) с ограничениями (2)-(3) представляет собой частично целочисленную задачу нелинейного бинарного программирования. При этом данная задача принадлежит к классу NP - полных задач. Из чего следует, что все без исключения методы решения задач данного типа имеют экспоненциальную оценку зависимости трудоемкости решения от ее размерности. Трудность решения таких задач основана не столько на не качественности методов и малой производительности современных электронно-вычислительных машин, сколько на принципиальной сложности самих задач. К классическим методам решения задач целочисленного нелинейного программирования относятся комбинаторные методы, а также метод ветвей и границ [5]. Однако оптимальное решение для метода ветвей и границ может быть получено задолго до остановки алгоритма, но обнаружить это нельзя, т. к. оптимальность устанавливается только по исчерпании списка задач. Использование же конечного множества решений в комбинаторных методах и замене полного перебора сокращенным, путем отсеивания неперспективных подмножеств решений, заведомо не содержащих оптимальных решений, носит эвристический характер и не гарантирует принятия оптимального решения. Согласно приведенным недостаткам основных подходов к решению задач нелинейного программирования нами был выбран метод полного перебора для поиска оптимальной конфигурации системы анодных секций для «прямого» (4) и «обратного» (5) режимов, при этом отыскание длительностей τпрям и τобр осуществляется методом Хука - Дживса. Выбор данного метода обосновывается требованиями к минимальному количеству вычислений критерия неравномерности покрытия посредством решения уравнения математической модели для поиска данных управляющих воздействий. На рис. 2 представлена схема алгоритма поиска решения задачи оптимального управления (1) реверсивным режимом нанесения гальванического покрытия в многоанодной ванне. Рис. 2. Алгоритм поиска решения задачи оптимального управления предложенным режимом нанесения гальванического покрытия После ввода исходных данных (размерности системы анодных секций и ее конфигурации, а также геометрии изделия-катода и значений ограничений на минимальную заданную толщину покрытия (2) и длительность процесса (3)) начинаются два цикла по формированию матриц активных анодных секций для «прямого» (4) и «обратного» (5) режимов нанесения покрытия предлагаемым режимом. Число итераций в каждом цикле составляет 2MN раз, что объясняется возможностью принятия только бинарных значений для элементов матриц. При этом подбор длительностей τпрям, τобр осуществляется согласно выбранному алгоритму. После формирования матриц анодных секций и подбора длительностей реверсивного режима согласно формуле (1) рассчитывается критерий неравномерности распределения толщины гальванического покрытия для i, j-го сочетания активных анодов (первый индекс соответствует i-му сочетанию активных анодов для «прямого» включения, второй - j-му сочетанию для «обратного» включения, где i, j = 1, 2, ..., 2MN). После расчета всех значений критериев неравномерности выбираются те оптимальные управляющие воздействия, которые обеспечат его минимум. Математическое моделирование и оптимальное управление реверсивным режимом нанесения гальванического покрытия в многоанодной ванне Согласно описанным выше алгоритмам, приведем пример математического моделирования и оптимального управления для гальванического процесса нанесения никелевого покрытия для плоского стального катода Sk квадратной формы площадью 900 см2 (L1 L2 = 30 30) и 9-ти анодной гальванической ванны, функционирующей в режиме реверсирования тока (см. на рис. 1). Гальваническая ванна размерами X Y Z = 40 40 45 см заполнена по уровню 40 см стандартным электролитом никелирования Уоттса следующего состава: NiSO4 ∙ 7H2O - 280…300 г/л, NiCl2 ∙ 6H2O - 50…70 г/л, H3BO3 - 30…40 г/л. В качестве анодных секций использовались никелевые квадратные пластины размерами 10 10 см, с одинаковыми межанодными расстояниями - 2 см (рис. 3). При этом расстояние между системой анодных секций и катодом составляет 30 см. Рис. 3. Схема размещения анодных секций в ванной (сечение): X, 0, Z - сечение, перпендикулярное дну ванны и проходящее через систему секций Для моделирования описанного процесса были заданы следующие параметры [6]: - напряжение для «прямого» и «обратного» тока 4 В; - удельная проводимость электролита 0,2 (Ом ∙ см)-1; - анодный и катодный выходы по току равные 0,95. Согласно ограничениям на минимальную заданную толщину δzad = 10 мкм (2) и максимальную продолжительность процесса Тmax = 7 200 с (3), найдены матрицы активных анодных секций (4), (5), обеспечивающие оптимальный режим функционирования анодных секций для «прямого» (рис. 4, а) и «обратного» (рис. 4, б) режимов нанесения гальванического покрытия. На рис. 4, а, б единичному значению соответствуют следующие элементы матрицы активных анодных секций: для «прямого» включения - , , , , ; для «обратного» включения - (рис. 4, б). Остальные элементы находятся в отключенном состоянии. При этом найденные оптимальные длительности включения реверсивного режима составляют = 7 с, = 2 с, а количество циклов переключения υ = 800 раз. а б Рис. 4. Динамика изменения включения активных анодных секций: а - для «прямого» режима нанесения гальванического покрытия в многоанодной ванне; б - для «обратного» режима нанесения гальванического покрытия в многоанодной ванне На рис. 5, а, б представлено распределение толщины покрытия δ, по поверхности катода для найденного оптимального и классического реверсивного режимов нанесения гальванического покрытия. Эффективность применения предлагаемого нами способа модификации реверсивного режима может быть рассчитана по формуле: , (6) где R1, R2 - критерии неравномерности распределения покрытия для классического и модифицированного режимов нанесения гальванического покрытия соответственно. а б Рис. 5. Распределение толщины покрытия δ, мкм, по поверхности катода, для различных режимов реверсирования: а - для предложенного режима; б - для классического режима Согласно (1), рассчитанные критерии неравномерности распределения толщины покрытия для классического режима реверсирования с монолитным анодом и его модификации составили R1 = 3,154 и R2 = 2,428 соответственно. При этом вычисленная по формуле (6) эффективность предложенного способа по сравнению с классическим режимом реверса тока составила 23 %, что подтверждает применимость предлагаемой модификации существующего режима реверсирования тока для гальванической ванны со многими анодами. Заключение Описаны существующие режимы функционирования многоанодных гальванических ванн, а также предложена модификация режима реверсирования тока для многоанодной ванны. Рассмотрены алгоритмы решения системы уравнений математической модели, описывающей предложенный процесс, а также метод поиска оптимального управления процессом по критерию неравномерности распределения гальванического покрытия по поверхности катода. Доказана эффективность применения предлагаемой модификации реверсивного режима для многоанодной ванны по сравнению с классическим режимом реверсирования тока в гальванической ванне с монолитным анодом на примере нанесения никелевого покрытия в электролите Уоттса для плоского катода и 9-анодной ванны. Для рассмотренного примера эффективность составила 23 %. Дальнейшие исследования будут посвящены разработке программно-аппаратной части автоматизированной системы управления предложенным технологическим процессом.
1. Milovanov I. V. Modelirovanie i optimizaciya tokovyh rezhimov v processah naneseniya gal'vanopokrytiy / I. V. Milovanov // Vestn. Tambov. gos. tehn. un-ta. 2002. T. 8. № 4. S. 603-611.
2. Litovka Yu. V. Simulation and optimization of electroplating with current reversal / Yu. V. Litovka, A. V. Romanenko, A. V. Afanasev // Theoretical Foundations of Chemical Engineering. 1998. Vol. 32. No. 3. P. 266-269.
3. Konkina V. V. Matematicheskaya postanovka zadachi optimal'nogo upravleniya gal'vanicheskim processom v reversivnom rezhime dlya mnogoanodnoy vanny / V. V. Konkina, D. S. Solov'ev // Teoreticheskie i prakticheskie aspekty tehnicheskih nauk: sb. st. Mezhdunar. nauch.-prakt. konf. Ufa: Aeterna. 2014. S. 41-43.
4. Zenkevich O. Konechnye elementy i approksimaciya / O. Zenkevich, K. Morgan. M.: Mir, 1986. 318 s.
5. Bandi B. Metody optimizacii. Vvodnyy kurs / B. Bandi. M.: Radio i svyaz', 1988. 126 s.
6. Ekilik V. V. Elektrohimicheskie metody zaschity metallov: metod. posob. po speckursu / V. V. Ekilik. Rostov n/D: Izd-vo RGU, 2004. 50 s.
