CC BY
38ПРИБОРОСТРОЕНИЕ, МЕТРОЛОГИЯ И ИНФОРМАЦИОННО-ИЗМЕРИТЕЛЬНЫЕ ПРИБОРЫ И СИСТЕМЫ
ИНФОРМАЦИОННО-ИЗМЕРИТЕЛЬНЫЕ И УПРАВЛЯЮЩИЕ СИСТЕМЫ
УДК 519.853, 004.942
ПОИСК СИЛЫ ТОКА И СТЕПЕНИ ЗАГРУЗКИ БАРАБАНА ДЛЯ СНИЖЕНИЯ НЕРАВНОМЕРНОСТИ ГАЛЬВАНИЧЕСКОГО ПОКРЫТИЯ МЕЛКИХ ДЕТАЛЕЙ1
Статья поступила в редакцию 21.03.2022, в окончательном варианте - 25.04.2022.
Соловьева Инна Александровна, 'Тамбовский государственный технический университет, 392000, Российская Федерация, г. Тамбов, ул. Советская, 106; 2Тамбовский государственный университет имени Г.Р. Державина, 392036, Российская Федерация, г. Тамбов, ул. Интернациональная, 33, 'аспирант, 2ассистент, ORCID: 0000-0002-1798-1859, e-mail: good.win32@yandex.ru Соловьев Денис Сергеевич, Тамбовский государственный университет имени Г.Р. Державина, 392036, Российская Федерация, г. Тамбов, ул. Интернациональная, 33,
кандидат технических наук, доцент, ORCID: 0000-0001-6613-3218, e-mail: solovjevdenis@mail.ru Литовка Юрий Владимирович, Тамбовский государственный технический университет, 392000, Российская Федерация, г. Тамбов, ул. Советская, 106,
доктор технических наук, профессор, ORCID: 0000-0002-5792-340X, e-mail: polychem@list.ru
В работе проводится анализ существующих исследований по оптимизации конструктивных параметров барабанов и режимов электролиза для снижения неравномерности гальванического покрытия мелких деталей. На основе проведенного анализа формулируется задача поиска силы тока и степени загрузки барабана для минимизации неравномерности получаемого покрытия. Для поставленной задачи предлагается математическая модель, которая учитывает глубину проникновения тока внутрь барабана и потери металла на истирание. Приводится анализ решения задачи оптимизации путем ее декомпозиции на две подзадачи, решаемые методами «золотого сечения» и пристрелки. Демонстрация решения задачи оптимизации приводится для никелирования установочного цилиндрического штифта M 5x36.
Ключевые слова: неравномерность покрытия, гальванический процесс, барабан, мелкие детали, математическая модель, оптимизация
SEARCH THE CURRENT AND THE DEGREE OF BARREL-SHAPED CAGE LOADING TO REDUCE THE NON-UNIFORMITY FOR THE PLATING COATING OF SMALL DETAILS
The article was received by the editorial board on 21.03.2022, in the final version — 25.04.2022.
Solovjeva Inna A., 1Tambov State Technical University, 106 Sovetskaya St., Tambov, 392000, Russian Federation; 2Tambov State University named after G.R. Derzhavin, 33 Internatsionalnaya St., Tambov, 392036, Russian Federation,
postgraduate student, 2Assistant, ORCID: 0000-0002-1798-1859, e-mail: good.win32@yandex.ru Solovjev Denis S., Tambov State University named after G.R. Derzhavin, 33 Internatsionalnaya St., Tambov, 392036, Russian Federation,
Cand. Sci. (Engineering), Associate Professor, ORCID: 0000-0001-6613-3218, e-mail: solovjevdenis@mail.ru
Litovka Yuri V., Tambov State Technical University, 106 Sovetskaya St., Tambov, 392000, Russian Federation,
Doct. Sci. (Engineering), Professor, ORCID: 0000-0002-5792-340X, e-mail: polychem@list.ru
The article provides an analysis for existing research on optimizing the design parameters of barrel-shaped cages and electrolysis modes to reduce the coating non-uniformity of small details. The problem of searching the current and the degree of barrel-shaped cage loading to minimize the non-uniformity of the plating coating is formulated on the analysis basis. A mathematical model that considers the current penetration depth into the load and metal losses due to abrasion is proposed to solve the problem. The analysis to the optimization problem solution is carried out by decomposition into two subproblems, which are solved by the golden section and shooting methods. Demonstration of the optimization problem solution is given for nickel plating of the locating cylindrical pin M 5x36.
Keywords: non-uniformity coating; plating process; barrel-shaped cage; small details; mathematical model; optimization
1 Исследование выполнено при финансовой поддержке РФФИ в рамках научного проекта № 20-37-90019.
Graphical annotation (Графическая аннотация)
Введение. Гальванические покрытия используются для защиты деталей от коррозии, придания им специальных свойств и улучшения внешнего вида [1]. К основным факторам, определяющим качество получаемого покрытия, относится его равномерность [2]. Как правило, гальваническое покрытие наносится путем «подвешивания» покрываемой металлической детали в растворе электролита. Однако данный подход имеет следующие недостатки: трудоемкость «подвешивания» больших объемов мелких деталей; использование подвесок может оставить следы на деталях, что приведет к ухудшению качества гальванического покрытия. Поэтому гальванические покрытия на мелкие детали наносят, как правило, в насыпном виде, загрузив их в ванну с барабаном [3]. Вращающиеся барабаны используются для избавления от трудоемких операций монтажа деталей на подвесочные приспособления, которые снижают производительность технологического процесса в целом. Условие получения покрытия на деталях с использованием вращающихся барабанов имеет существенные отличия от получения покрытия с использованием подвесок [4-6]. В процессе нанесения покрытия детали непрерывно перемешиваются, то есть периодически они могут оказываться на внешней поверхности или внутри загрузки. На поверхности деталей изменяется плотность тока из-за изменения их положения относительно анодов. Вследствие потери контакта с катодными токоподводами или массой загруженных деталей они могут оказываться неполяризо-ванными. Кроме того, происходит механическое взаимодействие деталей друг с другом. Поэтому каждая из одновременно покрываемых деталей может иметь в конце процесса случайное значение толщины, которое обусловливается влиянием совокупности факторов. Исследованию влияния факторов на неравномерность толщины получаемого таким способом покрытия посвящено множество работ. В работе [7] исследуется влияние диаметра барабана. В работе [8] исследуется влияние количества отверстий и их размеров на поверхности барабана. В работе [9] исследуется влияние температуры и состава электролита. В работе [10] исследуется влияние скорости вращения барабана и длительности гальванического процесса. Из проведенного анализа следует, что отсутствуют работы, посвященные исследованию влияния силы тока и степени загрузки барабана с использованием математических моделей процесса, которые учитывают глубину проникновения тока внутрь барабана и потери металла на истирание.
Целью статьи является поиск силы тока и степени загрузки барабана для снижения неравномерности гальванического покрытия мелких деталей.
Материалы и методы. Гальванический процесс нанесения покрытия предполагает размещение мелких деталей внутри барабана, изготавливаемого обычно из полипропилена. Барабан содержит центральные стержни, которые проводят электрический ток. Барабан с деталями медленно вращается при погружении в электролитический раствор, тем самым обеспечивая перекатывание деталей. При этом ион металла восстанавливается до металлического покрытия на поверхности деталей под действием электрического поля. В свою очередь свежий раствор электролита снаружи барабана непрерывно поступает внутрь через небольшие отверстия в его стенке, а прореагировавшие раствор и водород, образующийся в процессе нанесения покрытия, также выводятся из барабана через них. Сущность гальванического процесса нанесения покрытия на мелкие детали с использованием барабана показана на рисунке 1.
Рисунок 1 - Гальванический процесс нанесения покрытия на мелкие детали с использованием барабана
Критерием оптимизации технологического процесса выступает неравномерность толщины покрытия мелких деталей в барабане. Неравномерность толщины гальванического покрытия определяется согласно:
--Kh¡ Jlr,
(1)
где К - коэффициент разброса; к - средняя толщина покрытия; Дт - длительность процесса.
Для связи критерия оптимизации с искомой силой тока и степенью загрузки барабана составим математическую модель.
Средняя толщина покрытия определяется согласно закону Фарадея:
h = ^iflAr/[p(i + £)],
(2)
где к - электрохимический эквивалент; п - выход по току; га - анодная плотность тока; р - плотность; £ - коэффициент потерь на истирание (0.2.. .0.4 - для мягких металлов, 0.4.. .0.6 - для твердых металлов).
В качестве допущения предположим: поверхность внутри барабана представляет собой идеально гладкий цилиндр и загруженные в него детали движутся в режиме переката вокруг центра в точке С (рис. 2). Последнее допущение применимо в случае, когда барабан загружен деталями менее чем на половину.
Рисунок 2 - Схема переката деталей внутри барабана
Коэффициент разброса определяется согласно:
K = y¡a(]Df -1)(1+360®D sin 0.5a/(av))/(360®£), (3)
где а - угол к крайним точкам AB области с деталями в барабане; D - диаметр барабана; f - фактор загрузки; w - скорость вращения барабана; у - параметр распределения тока; Е - эффективность перемешивания; v - скорость падения деталей.
Угол к крайним точкам AB области с деталями в барабане определяется согласно:
a = arcsinsin2Cí п + sin2Cí п , (4)
где Ct - коэффициент загрузки барабана.
Параметр распределения тока определяется согласно:
, (5)
где 8С - площадь детали; Р - постоянная Фарадея; Ус - объем детали; г - количество электронов; К - универсальная газовая постоянная; х - удельная электропроводность электролита; ¡0 - плотность тока обмена; Т - температура электролита.
Фактор загрузки определяется согласно:
f = 0.125(«/360-0.5/пsma-Vd/Vb)/(/0 ■ а/360-1/п8т0.5а), (6)
где / - коэффициент перфорации стенок барабана; Уа, Уь - объемы катодных контактов и барабана. Эффективность перемешивания определяется согласно:
E=4(1-(2/3-ео80.5а+Ш)2 )/^(а-та)), (7)
где Ь - длина барабана.
Объем барабана определяется согласно:
Vb = 0.25^ф2 . (8)
С одной стороны, коэффициент загрузки барабана определяется согласно:
с, = М/ (РпУЬ ), (9)
с другой стороны:
С, = №, (10)
где М - масса загрузки; рт - насыпная плотность; V - объем загрузки. Количество деталей определяется согласно:
п = [М/МС ], (11)
где Мс - масса детали.
Площадь загрузки определяется согласно:
Sz = пБс. (12)
При размещении в барабане между деталями образуются зазоры, в связи с чем засыпанные детали представляют собой пористый электрод. Вследствие этого распределение тока неравномерно по всей структуре расположения деталей. Для получения наиболее равномерного покрытия необходимо учитывать факт проникновения тока внутрь загрузки. В качестве допущения примем, что детали, засыпанные в барабан, представляют собой электрод, сечением которого является окружность - эквивалентный электрод (рис. 3а). Тогда площадь деталей в активной зоне (зоне проникновения тока) определяется согласно:
Ба = кЖЬБ2 (ф-Ж , (13)
где Ж - глубина проникновения тока.
Плотность тока на глубине Ж от внешней поверхности электрода определяется согласно:
/ = 1а ехр [- Ж^З^Ж )) , (14)
где д - перенапряжение на границе «электролит - металл»; га - радиус эквивалентного электрода.
Сила тока через площадь деталей в активной зоне определяется согласно закону Ома в дифференциальной форме:
I = . (15)
Из (14) следует, что глубина проникновения тока внутрь загрузки зависит от отношения площади деталей в активной зоне к объему загрузки. Тогда глубина проникновения тока в барабанах меньшего размера будет выше, чем в барабанах большего размера при одинаковом объеме загрузки (рис. 3б). Таким образом, распределение плотности тока в барабане большего размера будет более неравномерным по сравнению с барабаном меньшего размера.
Радиус эквивалентного электрода определяется согласно:
. (16) Перенапряжение на границе «электролит - металл» определяется согласно уравнению Тафеля:
З = 2.3 RT \п(1а/1)/^). (17)
Известно, что в электролитах с высокой рассеивающей способностью получается гальваническое покрытие с наименьшей неравномерностью. Рассеивающая способность определяется согласно:
Н = (Щ )1(гХ }д81д1. (18)
На поверхности пористого электрода при Н > 1 распределение плотности тока близко к равномерному. При Н < 1 из-за концентрации плотности тока вблизи свободной поверхности пористого электрода равномерность снижается.
Рисунок 3 - Поперечное сечение барабана с эквивалентным электродом (а) и глубина проникновения тока внутрь загрузки при уменьшении размеров барабана (б)
Сформулируем задачу оптимизации. Найти силу тока I и коэффициент загрузки Ct деталей в барабане, для которых неравномерность (1) гальванического покрытия будет минимальной при уравнениях (2)-(18) и ограничениях:
Ат<ттах, (19)
h ^ hmin , (20)
Mmin < M < Mmax , (21)
i < i < i , (22)
min — a — max' V ^
где Tmax - максимальная длительность процесса (определяется производительностью гальванической линии в целом); Mmin - минимальная масса загрузки (определяется заданной производительностью барабана); Amin - минимальная толщина покрытия, несоблюдение которой приводит к браку детали; Mmax -максимальная масса загрузки (обусловлена особенностью конструкции барабана); imax - предельно допустимая плотность тока, выше которой качество покрытия снижается; imjn - минимальная плотность тока, ниже которой гальваническое покрытие не осаждается.
Проведем анализ влияния ограничений (19), (20) на критерий (1). Согласно уравнению (3), коэффициент разброса не зависит от времени. Тогда с учетом Дт = const минимальное значение неравномерности
толщины гальванического покрытия (1) будет получаться при h = hmin в (2). Следовательно, (1) примет
вид О = const/-JAt и будет иметь минимальное значение при Дт = Tmax. Коэффициент загрузки (3) не зависит от силы тока I, а зависит от степени загрузки Ct. Для расчета I по (15), согласно (13), (14), необходимо знать глубину проникновения тока внутрь загрузки W. Однако для расчета распределения плотности тока i в (14) необходимо знать распределение перенапряжения S, которое, согласно (17), также зависит от плотности тока. Поэтому для расчета распределения перенапряжения используется формула (15), которая получает первоначальное распределение плотности тока.
Таким образом, задача оптимизации разбивается на две подзадачи: внешнюю и внутреннюю. Внешняя подзадача сводится к одномерной задаче поиска минимума коэффициента разброса одним из методов одномерной оптимизации на отрезке [ Mmin/ (pmVb ), Mmax/(pmVb ) ]. По найденному минимуму коэффициента разброса во внутренней подзадаче определяется распределение перенапряжения одним из методов решения краевых задач и осуществляется поиск значения глубины проникновения тока внутрь загрузки, исходя из распределения плотности тока и перенапряжения.
Экспериментальная часть. Решение задачи оптимизации осуществим на примере поиска силы тока и коэффициента загрузки барабана для никелирования установочного цилиндрического штифта M 5 x 36. В таблице 1 приводятся значения параметров, используемых в расчете уравнений (2)-(18).
В таблице 2 приводятся значения ограничений (19)-(22), используемых при решении задачи оптимизации неравномерности гальванического покрытия.
Для решения внешней подзадачи использовался метод «золотого сечения» [11]. Для решения внутренней подзадачи использовался метод пристрелки [12]. Поиск решения задачи оптимизации осуществлялся в разработанном программном обеспечении [13].
Таблица 1 - Значения параметров
Параметр Значение
Барабан L, см 60
D, см 22,5
ю, об/ч 480
fo 0,22
Vd, см3 800
v, см/ч 31500
Деталь Sc, см2 6,04
Vc, см3 0,71
Mc, г 5,53
Pm, г/см3 5
Гальванический процесс ¿о, А/см2 0,001
0,4
P, г/см3 8,9
k, г/(А-ч) 1,095
Н 0,97
Z 2
T, К 313
X, 1/(Ом-см) 0,5
F, Кл/моль 96484,56
R, Дж/(моль-К) 8,31
Таблица 2 - Значения ограничений
Параметр Значение
Барабан Mmm, г 10000
Mm„, г 40000
Гальванический процесс Amin, МКМ 8
^max? ч 0,5
¿min, А/СМ2 0,015
¿max, А/СМ2 0,030
Результаты и их обсуждение. На рисунке 4 показаны зависимости плотности тока и рассеивающей способности электролита от глубины проникновения тока внутрь загрузки. На глубине Ж > 0,93 см от внешней поверхности 8а плотность тока г отсутствует, а рассеивающая способность электролита снижается (Н < 1), тем самым увеличивая неравномерность покрытия.
0 ОД 0,4 0,6 0,S 1 0 0,2 0,4 0,6 0,6 1
Рисунок 4 - Зависимость плотности тока (слева) и рассеивающей способности электролита (справа) от глубины проникновения тока внутрь загрузки
На рисунке 5 показана область допустимых решений (располагается между прямыми С,(МШш) и С,(Мшшх), /ш1п и /шах) и оптимальное решение (С,*; /*) для задачи оптимизации.
LA Ct(M™¿J Ct(bUJ
0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4
Рисунок 5 - Область допустимых решений и оптимальное решение задачи оптимизации
В таблице 3 приводятся найденные значения параметров в результате решения задачи оптимизации неравномерности гальванического покрытия.
Таблица 3 - Найденные значения параметров
Па раметр Значение
Барабан Ct 0,21
Sz, см2 19270
M, г 250б2
Гальванический процесс ia, А/см2 0,002
Sa, см2 1379
/, А 27,б4
а, мкм 0,508
Таким образом на n = [M/Mc] = [25G62/5,53] = 4532 деталях установочного цилиндрического штифта
M 5 x Зб средняя толщина никелевого покрытия составит h = Km + G = 8 + °.5°В = 8.5°8 мкм.
Заключение. Проблема получения гальванического покрытия с заданной толщиной и минимальной неравномерностью возникает при управлении технологическим процессом. Для гальванического процесса нанесения покрытий на мелкие детали с использованием барабанов имеется ряд специфических особенностей, связанных с необходимостью улучшения качества наносимых покрытий в условиях ужесточающихся требований к ним. Решение проблемы усложняется противоречивостью экспериментальных исследований о воздействии параметров процесса и барабана на равномерность покрытия. Выходом из сложившейся ситуации является оптимизация параметров гальванического процесса и барабана с использованием разработанной математической модели, которая учитывает глубину проникновения тока внутрь барабана и потери металла на истирание. Улучшение качества наносимых покрытий повлечет за собой уменьшение брака и перерасхода металла, что в итоге снизит себестоимость обработки деталей.
Библиографический список
1. Гамбург, Ю. Д. Теория и практика электроосаждения металлов / Ю. Д. Гамбург, Дж. Зангари. - Москва : БИНОМ. Лаборатория знаний, 2G15. - 438 с.
2. Помогаев, В. М. Прогнозирование равномерности осаждения гальванических покрытий / В. М. Помогаев, А. В. Волкович, И. В. Петроченкова, А. Е. Шувакин // Известия высших учебных заведений. Серия: Химия и химическая технология. - 2GG7. - Т. 5G, № 3. - С. 1G3-1G6.
3. Singleton, R. Fundamentals of barrel plating / R. Singleton, E. Singleton // Metal Finish. - 2G11. - Vol. 1G9 (8). -P. 14-22.
4. Purcar, M. 3D electrode shape change simulation in electroplating / M. Purcar, C. Munteanu, V. Topa // Revue Roumaine des Sciences Techniques. Serie Electrotechnique et Energetique. - 2G13. - Vol. 58. - P. 252-262.
5. Лютов, А. Г. Моделирование процесса нанесения гальванических покрытий с учетом геометрических конфигураций электродов / А. Г. Лютов, А. Р. Ишкулова // Вестник Уфимского государственного авиационного технического университета. - 2G15. - Т. 19, № 4 (7G). - С. 45-48.
6. Соловьев, Д. С. Применение решения обратной задачи математического моделирования гальванического процесса для оптимизации неравномерности толщины покрытия / Д. С. Соловьев, И. А. Соловьева, Ю. В. Литовка // Прикаспийский журнал: управление и высокие технологии. - 2G2G. - № 1 (49). - С. 131-143.
7. Караваев, В. И. Оптимизация гальванических процессов в ваннах с барабанами с целью повышения равномерности покрытия / В. И. Караваев, Ю. В. Литовка // Журнал прикладной химии. - 2004. - Т. 77, № 9. - С. 1481-1486.
8. Han, K. H. Barrel plating process specification for undercoating with copper cyanate / K. H. Han, S. B. Lee, I. K. Hong // Journal of Industrial and Engineering Chemistry. - 2012. - Vol. 18 (3). - P. 888-897.
9. Hong, I. K. Optimization of barrel plating process for electroless Ni-P plating / I. K. Hong, H. Kim, S. B. Lee // Journal of Industrial and Engineering Chemistry. - 2014. - Vol. 20 (5). - P. 3767-3774.
10. Yanagishita, T. Fabrication of Porous Si Particles by Barrel Anode Etching / T. Yanagishita, S. Ueno, K. Taniguchi, H. Masuda // Chemistry Letters. - 2016. - Vol. 45 (7). - P. 708-710.
11. Измаилов, А. Ф. Численные методы оптимизации / А. Ф. Измаилов, М. В. Солодов. - Москва : Физматлит, 2008. - 320 c.
12. Самарский, А. А. Численные методы математической физики / А. А. Самарский, А. В. Гулин. - Москва : Научный мир, 2000. - 358 c.
13. Соловьева, И. А. Оптимизация гальванических процессов для деталей, покрываемых насыпью в барабанах / И. А. Соловьева, Д. С. Соловьев, Ю. В. Литовка // Свидетельство об официальной регистрации программы для ЭВМ № 2022613030, дата регистрации 01.03.2022.
References
1. Gamburg, Yu. D., Zangari, G. Teoriya ipraktika elektroosazhdeniya metallov [Theory and practice of metal elec-trodeposition]. Moscow, BINOM. Laboratoriya znaniy Publ., 2015. 438 p.
2. Pomogaev, V. M., Volkovich, A. V., Petrochenkova, I. V., Shuvakin, A. E. Prognozirovanie ravnomernosti osa-zhdeniya galvanicheskih pokrytiy [Uniformity of electroplate deposition prediction]. Izvestiya vysshikh uchebnykh zavedeniy. Seriya: Khimiya i khimicheskaya tekhnologiya [News of Higher Educational Institutions. Series: Chemistry and chemical technology], 2007, vol. 50, no. 3, pp. 103-106.
3. Singleton, R., Singleton E. Fundamentals of barrel plating. Metal Finish, 2011, vol. 109 (8), pp. 14-22.
4. Purcar, M., Munteanu, C., Topa, V. 3D electrode shape change simulation in electroplating. Revue Roumaine des Sciences Techniques. Serie Electrotechnique etEnergetique, 2013, vol. 58, pp. 252-262.
5. Lutov, A. G., Ishkulova, A. R. Modelirovanie protsessa naneseniya galvanicheskikh pokrytiy s uchetom geometrich-eskikh konfiguratsiy elektrodov [Modeling of plating process in view of geometrical configurations of electrodes]. Vestnik Ufimskogo gosudarstvennogo aviatsionnogo tekhnicheskogo universiteta [Bulletin of the Ufa State Aviation Technical University], 2015, vol. 19, no. 4 (70), pp. 45-48.
6. Solovyev, D. S., Solovyeva, I. A., Litovka, Yu. V. Primenenie resheniya obratnoy zadachi matematicheskogo mo-delirovaniya galvanicheskogo protsessa dlya optimizatsii neravnomernosti tolshchiny pokrytiya [Inverse problem solution of the mathematical modeling for galvanic process to optimize the non-uniformity of the coating thickness]. Prikaspiyskiy zhur-nal: upravlenie i vysokie tekhnologii [Caspian Journal: Control and High Technologies], 2020, no. 1 (49), pp. 131-143.
7. Karavaev, V. I., Litovka, Yu. V. Optimizatsiya galvanicheskikh protsessov v vannakh s barabanami s tselyu po-vysheniya ravnomernosti pokrytiya [Optimization of electroplating processes in drum baths to obtain a more uniform coating]. Zhurnalprikladnoy khimii [Russian Journal of Applied Chemistry], 2004, vol. 77, no. 9, pp. 1481-1486.
8. Han, K. H., Lee, S. B., Hong, I. K. Barrel plating process specification for undercoating with copper cyanate. Journal of Industrial and Engineering Chemistry, 2012, vol. 18 (3), pp. 888-897.
9. Hong, I. K., Kim, H., Lee, S. B. Optimization of barrel plating process for electroless Ni-P plating. Journal of Industrial and Engineering Chemistry, 2014, vol. 20 (5), pp. 3767-3774.
10. Yanagishita, T., Ueno, S., Taniguchi, K., Masuda, H. Fabrication of Porous Si Particles by Barrel Anode Etching. Chemistry Letters, 2016, vol. 45 (7), pp. 708-710.
11. Izmailov, A. F., Solodov, M. V. Chislennye metody optimizatsii [Numerical optimization methods]. Moscow, Fiz-matlit Publ., 2008. 320 p.
12. Samarsky, A. A., Gulin, A. V. Chislennye metody matematicheskoy fiziki [Numerical methods of mathematical physics]. Moscow, Nauchnyy mir Publ., 2000. 358 p.
13. Solovyeva, I. A., Solovyev, D. S., Litovka, Yu. V. Optimizatsiya galvanicheskikh protsessov dlya detaley, pokryvaemykh nasypyu v barabanakh [Optimization of electroplating processes for parts coated in bulk barrels]. Svidetelstvo ob ofitsialnoy registratsii programm dlya EVM [The Certificate on official registration of the computer program], no. 2022613030, registration date 01.03.2022.
