CC BY
13
UDC 004.942
Doi: 10.31772/2587-6066-2020-21-3-424-432
For citation: Bocharova O. A., Murygin A. V, Bocharov A. N., Zaitsev R. V. Simulation of the induction soldering process of waveguide paths from aluminum alloys. Siberian Journal of Science and Technology. 2020, Vol. 21, No. 3, P. 424-432. Doi: 10.31772/2587-6066-2020-21-3-424-432
Для цитирования: Моделирование процесса индукционной пайки волноводных трактов из алюминиевых сплавов / О. А. Бочарова, А. В. Мурыгин, А. Н. Бочаров, Р. В. Зайцев // Сибирский журнал науки и технологий. 2020. Т. 21, № 3. С. 424-432. Doi: 10.31772/2587-6066-2020-21-3-424-432
SIMULATION OF THE INDUCTION SOLDERING PROCESS OF WAVEGUIDE PATHS
FROM ALUMINUM ALLOYS
O. A. Bocharova*, A. V. Murygin, A. N. Bocharov, R. V. Zaitsev
Reshetnev Siberian State University of Science and Technology 31, Krasnoyarskii rabochii prospekt, Krasnoyarsk, 660037, Russian Federation E-mail: shyx_89@mail.ru
A system of waveguide paths is a complex structure of various elements with various geometries. Induction soldering based on the induction heating method is one of the promising methods for waveguides fabricating. Induction soldering of waveguide paths has a number of technological features: the melting temperature of the base material AD31 (695-663 0C) slightly differs from the melting temperature of St. AK12 solder (577-580 0C) at an average induction heating rate of 20-25 0C / sec; a wide variety of standard sizes of waveguide paths elements complicates the development and subsequent reproduction of technological parameters of the induction soldering process; zones of maximum heating of waveguide paths elements do not coincide with zones of soldering. Therefore, to solve the problems of controlling the waveguides soldering process, it is necessary to simulate this process. The paper deals with the problem of simulating the process of heating a waveguide during induction soldering. Requirements for the process model have been formed. The model is built on the basis of the differential heat conduction equation. The formed model requirements take into account the geometric parameters of waveguides, the physical parameters of materials, the initial and boundary conditions, as well as the uneven distribution of eddy current density in the waveguide. It is proposed to use the finite difference method for the numerical solution of the heat conduction equation. The process of calculating the temperature at the grid nodes is shown. The authors propose a two-stage solution. At the first stage, at an intermediate time step, the temperature at the grid nodes along the X axis is calculated. At the second stage, the temperature at the grid nodes along the Y axis is calculated. The numerical solution of the difference equations along the X and Y axes is carried out by the sweep method. An algorithm for the numerical solution of the heat conduction equation has been developed.
Keywords: waveguide path, induction soldering, model of the waveguide heating process, differential heat conduction equation, finite difference method.
МОДЕЛИРОВАНИЕ ПРОЦЕССА ИНДУКЦИОННОЙ ПАЙКИ ВОЛНОВОДНЫХ ТРАКТОВ
ИЗ АЛЮМИНИЕВЫХ СПЛАВОВ
О. А. Бочарова*, А. В. Мурыгин, А. Н. Бочаров, Р. В. Зайцев
Сибирский государственный университет науки и технологий имени академика М. Ф. Решетнева Российская Федерация, 660037, г. Красноярск, просп. им. газ. «Красноярский рабочий», 31
E-mail: shyx_89@mail.ru
Система волноводных трактов представляет собой сложную конструкцию из различных элементов с разнообразной геометрией. Одним из перспективных способов изготовления волноводов является индукционная пайка, основанная на методе индукционного нагрева. Индукционная пайка волноводных трактов обладает рядом технологических особенностей: относительно небольшая разница температуры плавления основного материла АД31 (695-663 0С) и припоя Св. АК12 (577-580 0С) при средней скорости индукционного нагрева 20-25 0С/сек; большое разнообразие типоразмеров элементов волноводных трактов представляет сложность при отработке и последующем воспроизведении технологических параметров процесса индукционной пайки; зоны максимального нагрева элементов волноводных трактов не совпадают с зонами пайки. Поэтому для решения задач управления процессом пайки волноводов необходимо провести моделирование данного процесса.
В статье рассмотрена задача моделирования процесса нагрева волновода при индукционной пайке. Сформированы требования к модели. Модель строится на основе дифференциального уравнения теплопроводности. Сформированные требования к модели учитывают геометрические параметры волноводов, физические параметры материалов, начальные и граничные условия, а также неравномерное распределение плотности вихревого тока в волноводе. Предлагается для численного решения уравнения теплопроводности использовать метод конечных разностей. Показан процесс расчета температуры в узлах сетки. Решение осуществляется в два этапа. На первом этапе на промежуточном временном шаге проводится расчет температуры в узлах сетки по оси X, на втором этапе вычисляется температура в узлах сетки по оси У. Численное решение разностных уравнений по оси X и У осуществляется методом прогонки. Разработан алгоритм численного решения уравнения теплопроводности.
Ключевые слова: волноводный тракт, индукционная пайка, модель процесса нагрева волновода, дифференциальное уравнение теплопроводности, метод конечных разностей.
Introduction. Systems of waveguide paths are widely used in spacecraft to ensure their functioning and transmission of electromagnetic energy in them [1].
A system of waveguide paths is a complex structure of various elements with various geometries (fig. 1).
Various methods are used to connect the structural elements of waveguide paths: argon-arc welding, soldering in salt baths, soldering using a laser source [2-5].
Induction soldering is one of the promising methods for waveguide paths fabricating [1; 3].
Induction soldering provides contactless, fast and local heating of the soldering point. It reduces the oxidation of current-carrying surfaces and waveguide warping. The equipment providing induction heating is rather compact and easily controlled, which allows automated controlling the soldering process [6].
However, induction soldering of waveguide paths has a number of features [7]:
1. The melting temperature of the base material AD31 (695-663 0C) slightly differs from the melting temperature of St. AK12 solder (577-580 0C) at an average induction heating rate of 20-25 0C / sec.
2. A wide variety of standard sizes of waveguide paths elements complicates the development and subsequent reproduction of technological parameters of the induction soldering process.
3. Zones of maximum heating of the waveguide paths elements do not coincide with zones of soldering.
Therefore, it is necessary to simulate the process of heating the waveguide before setting the problem of controlling the process of waveguide paths induction soldering.
Fig. 1. An example of a section of a waveguide distribution system: a - a direct element; b - a curved element with a variable radius of curvature; c - a curved element with a constant radius of curvature; d - a coupling; e - a flange; f- a flexible section; g - an intermediate support
Рис. 1. Пример участка волноводно-распределительной системы: а - прямой элемент; б - криволинейный элемент с переменным радиусом кривизны; в - криволинейный элемент с постоянным радиусом кривизны; г - соединительная муфта; д - фланец; е - гибкая секция; ж - промежуточная опора
а б
Fig. 2. The researched system "inductor - waveguide": 1 - a waveguide pipe; 2 - a flange; 3 - an inductor
Рис. 2. Исследуемая система «индуктор - волновод»: 1 - волноводная труба; 2 - фланец; 3 - индуктор
Construction of a waveguide heating model. The
connection of the straight section of the waveguide tube and the flange was chosen as a typical object for considering the process of induction heating (fig. 2).
The inductor (without considering the power part of the induction equipment) is a ring made of copper (fig. 2). The inductor is a source of alternating magnetic field, which forms eddy currents flowing in the waveguide. The calculation of a simplified two-dimensional waveguide-inductor circuit showed that the eddy currents in the waveguide are unevenly distributed [8-10].
Based on the above, it is possible to construct a mathematical model of heating the waveguide path during induction soldering, taking into account the following assumptions:
- the spatial configuration of the system allows us to restrict ourselves to a simplified two-dimensional model;
- a general model can be used to simulate a pipe-flange system, since the waveguide elements can be represented by plates in a simplified form.
For a generalized model of thermal processes occurring during induction soldering of waveguide paths, we take the differential equation of heat conduction [11]. Since the heating of the waveguide is due to eddy currents flowing in it, a permanent heat source is added to the formula:
dT ( x, y, t) _
= a
dt
(d 2T ( x, y, t ) + d 2T ( x, y, t )Л
dx 2
5y2
where a is the coefficient of thermal diffusivity of the material; T is temperature; t is time; x, y are Cartesian coordinates; q(x, y, t) is a permanent source of heat.
Equation (1) describes a variety of options for the waveguide heating process flow. Therefore, in order to simulate the thermal processes occurring during induction soldering, it is necessary to add single-valued conditions to equation (1). Single-valued conditions may contain geometric, physical, initial and boundary conditions. These conditions are as follows:
- a simplified geometric two-dimensional model of the waveguide is adopted, shown in fig. 3;
- eddy currents are supposed to flow over the surface of the waveguide [12];
- it is taken into account that the eddy current density is distributed unevenly due to the shape of the waveguide and the location of the inductor relative to the waveguide [12];
- thermophysical characteristics of the waveguide material are constant;
- the initial temperature of the waveguide is the same along its entire length:
T(x,y) = Ti = const;
- heat exchange with the environment (Newton's boundary condition) is carried out by the inner and outer surface of the pipe:
ÔT ( x, y, t ).
dx
- = b(T (x, y, t) - Ten),
+ q(x, y, t), (1)
where b is the heat exchange coefficient; Ten is the temperature of the environment; 0 < x < lp (inner side); 0 < x < lp - hf (outer side);
Fig. 3. A simplified geometric model of the waveguide: lp - length of a pipe; hp - width of a pipe; hf - width of a flange; lf - length of a flange
Рис. 3. Упрощенная геометрическая модель волновода: 1Т - длина трубы; hT - ширина трубы; hf - ширина фланца; lf - длина фланца
Fig. 4. A difference mesh of the solution domain Рис. 4. Разностная сетка области решения
- heat exchange with the environment is carried out by the upper and side faces of the flange. For the upper face:
dT ( x, y, t)
where hp < y < hp
dy + lf-
■ = b(T(x, y, t) - Ten),
For the side face:
a dT (x, y, t)
dx
where lp - hf < x < lp;
= b(T(x, y, t) - Ten),
- heating power remains constant q(x, y, t) = const;
- heat exchange between the pipe and the flange is assumed to be zero.
For the numerical solution of equation (1), we apply the finite difference method [13-15]. Let us consider the method applied to a waveguide tube.
To approximate equation (1) by the finite difference method, we represent the pipe as a set of nodes with coordinates (fig. 4): xt = i-hx, y{ = j hy, tk = kz, where hx, hy is the step along the coordinate grid x, y; t is the time step.
Let us sample equation (1) in two stages. At the first stage, at an intermediate time step t / 2, we carry out sampling along the X axis; at the second stage, we carry out sampling along the Y axis. We obtain equation (1) in the difference form [15]:
2 _ Tk
Tij
t 2
Tk+1 _ Tk + 2 T.j T.j
k+1 k+1 k+1 ^ T 2 _ 2Tk+ 2 + Tk + 2 T+1, j 2Tij + T-1, j
+q
k+1
= a
f Tk+1 _ 2Tk+1 + Tk+1 ^ Ti, j+i _ 21i,j + Ti,j-i
+q
k+1
(2)
(3)
where Tk = T(^, y-, tk ) ; q^ = q(^, y-, tk ).
Solving first equation (2), and then equation (3), we determine the temperature field at a whole time step.
Approximation of the boundary conditions for a pipe is as follows:
Tk Tk
a^-1-^ = b(Tkj - Ten), if 0 < i < Nx, j = 0,
h
y
Tk Tk
a^j^-u_ = b(Tkj _ Ten). if 0 < < Nx. j = Ny.
h
y
k+1 Tk + 2 Ti, j , 2 —---+ a
f k+iA
Tk + 2 Ti. j
t
V У
T 2 Tk+1 Tk , 1
T+1. j , Ti_1.j , Ti,j , k+2
= a—+ a—+ + qi 2
h,2 h,2 T ^.
Then
1 2q
t h2
V 2 x y
Tk+1=aTk+2 + aTk+2 +1 Tk + qk+1
T i.j ~ h2 T i+1.j h2 1 i_1.j T 1 i.j qi.j '
Hence we get
rpk+2 , ^k+i k+| k+2
a1T1, j 2 = brTr +w + CiTi_w + di + % j 2.
(4)
where
aoTok+2 = Vj2 + d0.
„k+
k+1
aNTN. / = CNTN/ + dN •
(5)
(6)
„k+
k+1
To. j2 = ад-.j2 + e.
0T i.j
(7)
where P0 = —, Q0 = — are determined by the initial
ao a0
and boundary conditions of the simulation.
We substitute relation (7) into (4) for i = 0. Continuing
k+1
the process of sequential substitution Tt -2 can be ex-
k+1
pressed through Ti+12 :
k+1 k+1 Tk+2 = ptk 2 + q Ti. j PiTi+1. j ei ,
(8)
where P,- h ^ are new coefficients obtained during the substitution process. Using expressions (8) and (4), we can obtain formulas for calculating p, h Qi:
Pi = —b-, (9)
Qi =
ai _ cip
d + c,Q _1
ai _ сгрг_1
(10)
When calculating PNx and QN we get that PNx = 0,
= 0. Therefore, TN will be equal to QN . Having
calculated TN in this way, one can start the backward-sweep process to obtain TN_i, TN-2, T2, Ti, T0 [9].
k+1 k+i Ti_+2 = Pi_\Ti_\j + Qi_1 •
(11)
The sweep method [14] is suitable for solving difference equations (2) and (3). Let us consider this method using the example of solving equation (2). We transform the equation (2) [14]:
Based on the above, we can formulate a solution algorithm:
1. We calculate P0 and Q0, the coefficients a0, b0, J0 are determined using the initial and boundary conditions.
2. We calculate P,- and Qi for i = 1, 2, ..., Nx, using expressions (9) and (10).
3. We suppose TN = QN.
4. We substitute TN into formula (11) and determine
TN—1, TN^ T2, T1, T0.
The algorithm for solving equation (3) is similar. The coefficients included in equation (3) will take the following form:
1 , a , a 1 k+,2
a =—, b = ^ , d; =— T 2
ht
h2
'j
1 , 2q ^ 1 k «i =T + bi = ci = ^ =TTk-.
2 hx hx 2
For the boundary points 0 and Nx, we write equation (4) in the following form:
The sweep algorithm begins by writing equation (5) in the form:
The sweep will be carried out on index j.
The Elcut software package was used to calculate the eddy current density distribution in the waveguide pipe and flange. The results confirm the previously accepted assumptions when setting the simulation problem. Fig. 5 and 6 show that the eddy current density increases sharply on the corner surfaces of the pipe and the waveguide flange. In addition, in fig. 5 the eddy current density decreases sharply in the thickness of the pipe and flange walls. Eddy current flow areas are well localized (fig. 6).
Fig. 7 shows the results of heating the waveguide according to the proposed model.
To check the simulation results, the temperature change at control points 1 and 2 was taken (fig. 7). The control points coincide with the points of temperature measurement during the process of heating the waveguide using induction equipment [16; 17].
As it can be seen from fig. 8 and 9, the simulation results quite accurately describe the process of the waveguide elements heating. Consequently, this model can be used to simulate an automated control system for soldering waveguide paths and analyzing control algorithms.
а b
Fig. 5. The distribution of pipe eddy current density for the size 35x15 mm: a - in the large side of the pipe; b - in the thickness of the pipe
Рис. 5. Распределение плотности вихревого тока трубы для типоразмера 35x15 мм: а - по большой стороне трубы; б - по толщине трубы
Fig. 6. The distribution of eddy current density in the "inductor - waveguide" system for the size 35x15 mm Рис. 6. Распределение плотности вихревого тока в системе «индуктор - волновод» для типоразмера 35x15 мм
Fig. 7. Temperature fields of the waveguide during induction heating Рис. 7. Температурные поля волновода при индукционном нагреве
Fig. 8. Heating of the 35 x 15 waveguide pipe: 1 - calculated values; 2 - measured data
Рис. 8. Нагрев трубы волновода типоразмера 35 x 15: 1 - расчетные значения; 2 - измеренные данные
Fig. 9. Heating of the 35 x 15 waveguide flange: 1 - calculated values; 2 - measured data
Рис. 9. Нагрев фланца волновода типоразмера 35 x 15: 1 - расчетные значения; 2 - измеренные данные
Conclusion. 1. The process of induction soldering of waveguide paths has a number of features, which complicates the selection of process parameters for various types of waveguide paths and process control and suggests the need to simulate the heating of waveguides based on the differential equation of heat conduction.
2. To simulate the process of the pipe and the waveguide flange heating, it is most expedient to use the finite difference method.
Acknowledgments. The reported study was funded by Russian Foundation for Basic Research, Government of Krasnoyarsk Territory, Krasnoyarsk Regional Fund of
Science, within the framework of the research project: "Mathematical and physical simulation of processes occurring during induction soldering of pipelines in protective environments", project № 18-48-242006.
Благодарности. Исследование выполнено при финансовой поддержке Российского фонда фундаментальных исследований, Правительства Красноярского края, Красноярского краевого фонда науки в рамках научного проекта: «Математическое и физическое моделирование процессов, происходящих при индукционной пайке трубопроводов в защитных средах», проект № 18-48-242006.
References
1. Zlobin S. K., Mikhnev M. M., Laptenok V. D., Bo-charov A. N., Dolgopolov B. B. [Features of production of waveguide-distribution paths of antenna-feeder devices of space vehicles]. Vestnik SibGAU. 2013, No 6, P. 196-201 (In Russ.).
2. Brovko A. V [Problems of automatic welding of radar waveguides]. Izvestiya vuzov: Mashinostroenie. 2013, No. 1, P. 50-54 (In Russ.).
3. Bushminsky I. P. Izgotovlenie ehlementov kon-struktsii SVCh. Volnovody i volnovodnye ustroistva [Manufacturing of elements of microwave structures. Waveguides and waveguide devices]. Moscow, Vysshaya shkola Publ., 1974, P. 304.
4. Full in-house production facilities. Available at: http ://www. advancedmicrowave. com/our-facilities (accessed: 10.05.2020).
5. Pamin S. et al. Joining of aluminum waveguides using pulsed laser radiation. Microwave Conference (APMC), 2015 Asia-Pacific. - IEEE, 2015, vol. 3, P. 1-3.
6. Rapoport E., Pleshivtseva Y. Optimal Control of Induction Heating Processes. CRC Press, NY, 2007, 349 p.
7. Zlobin S. K. [Features of soldering elements of waveguide-distribution paths from aluminum alloys with the use of an induction heating source]. Materialy XVI Mezhdunar. nauch. konf. "Reshetnevskie chteniya" [Materials XVI Intern. Scientific. Conf "Reshetnev reading"]. Krasnoyarsk, 2012, Vol. 1, P. 16-17 (In Russ.).
8. Bocharova O. A., Tynchenko V S., Bocharov A. N., Oreshenko T. G., Murygin A. V, Panfilov I. A. Induction heating simulation of the waveguide assembly elements. Journal of Physics: Conference Series. 2019, Vol. 1353, P. 012040.
9. Patidar B., Hussain M. M., Sanjoy Das, Mukherjee D, Tiwari A. P. Simulation and Experimental Validation of Induction Heating of MS Tube for Elevated Temperature.
NDT Application Excerpt from the Proceedings of the COMSOL Conference in Pune, 2015, 6 p.
10. Rhein S., Tilman U., Knut G. Optimal control of induction heating processes using FEM software. Conference: 2015 European Control Conference (ECC), 2015, P. 515-520.
11. Lykov A. V. Teoriya teploprovodnosti [Theory of thermal conductivity]. Moscow, Vysshaya shkola Publ., 1967, 599 p.
12. Babat G. I. Induktsionnyi nagrev metallov i ego promyshlennoe primenenie [Induction heating of metals
and its industrial application]. Moscow - Leningrad, Energy Publ., 1965, 552 p.
13. Paskonov V M., Polezhaev V I., Chudov L. A. Chislennoe modelirovanie protsessov teplo-massoobmena [Numerical modeling of heat and mass transfer processes]. Moscow, Nauka Publ., 1984, 288 p.
14. Patankar S. V., Kalabin E. V., Yankov G. G. Chislennoe reshenie zadach teploprovodnosti i konvektiv-nogo teploobmena pri techenii v kanalakh [Numerical solution of problems of thermal conductivity and convec-tive heat transfer during flow in channels]. Moscow, Mehi Publ., 2003, 312 p.
15. Samara A. A. Teoriya raznostnykh skhem [Theory of difference schemes]. Moscow, Nauka Publ., 1977, 388 p.
16. Zlobin S. K., Mikhnev M. M., Laptenok V. D., Seregin Yu. N., Bocharov A. N., Tynchenko V. S., Dubets Yu. P., Dolgopolov B. B. [Automated equipment and technology for soldering waveguide paths of spacecraft]. Vestnik SibGAU. 2014, No. 4 (56), P. 219-229 (In Russ.).
17. Murygin A. V, Tynchenko V. S., Laptenok V. D., Emilova O. A., Seregin Y. N. Modeling of thermal processes in waveguide tracts induction soldering. IOP Conference Series: Materials Science and Engineering. 2017, Vol. 173(1), P. 012026.
Библиографические ссылки
1. Особенности производства волноводно-распре-делительных трактов антенно-фидерных устройств космических аппаратов / С. К. Злобин, М. М. Михнев,
B. Д. Лаптёнок и др. // Вестник СибГАУ 2013. № 6.
C. 196-201.
2. Бровко А. В. Проблемы автоматической сварки волноводов радиолокационных станций // Известия вузов: Машиностроение. 2013. № 1. С. 50-54.
3. Бушминский И. П. Изготовление элементов конструкций СВЧ. Волноводы и волноводные устройства. Учебное пособие для вузов. М. : Высшая школа, 1974. 304 с.
4. Full in-house production facilities [Электронный ресурс]. URL: http://www.advancedmicrowave.com/our-facilities (дата обращения: 10.05.2020).
5. Pamin S. et al. Joining of aluminum waveguides using pulsed laser radiation // Microwave Conference (APMC), 2015 Asia-Pacific. - IEEE, 2015. Vol. 3. P. 1-3.
6. Rapoport, E., Pleshivtseva Y. Optimal Control of Induction Heating Processes. CRC Press. N. Y., 2007.
7. Особенности пайки элементов волноводно-распределительных трактов из алюминиевых сплавов с применением источника индукционного нагрева / С. К. Злобин, М. М. Михнев, В. Д. Лаптенок и др. // Решетневские чтения : материалы XVI междунар. на-учн. конф. : в 2 ч. Красноярск, 2012 . Ч. 1. С. 16-17.
8. Induction heating simulation of the waveguide assembly elements / O. A. Bocharova, V. S. Tynchenko, A. N. Bocharov et al. // Journal of Physics: Conference Series. IOP Publishing, 2019. Vol. 1353. P. 012040
9. Patidar, B. Simulation and Experimental Validation of Induction Heating of MS Tube for Elevated Temperature / B. Patidar, M. M. Hussain, Sanjoy Das et al. // NDT
Application Excerpt from the Proceedings of the COM-SOL Conference in Pune. 2015.
10. Rhein S., Tilman U., Knut G. Optimal control of induction heating processes using FEM software // Conference: 2015 European Control Conference (ECC). 2015, P. 515-520.
11. Лыков А. В. Теория теплопроводности. М. : Высшая школа, 1967. 599 с.
12. Бабат Г. И. Индукционный нагрев металлов и его промышленное применение. М. - Л. : Энергия, 1965. 552 с.
13. Пасконов В. М., Полежаев В. И., Чудов Л. А. Численное моделирование процессов тепло-массо-обмена. М. : Наука, 1984. 288 с.
14. Патанкар С. В. Численное решение задач теплопроводности и конвективного теплообмена при те-
чении в каналах : пер. с англ. Е. В. Калабина ; под ред. Г. Г. Янькова. М. : Изд-во МЭИ, 2003. 312 с.
15. Самарский А. А. Теория разностных схем. М. : Наука, 1977. 388 с.
16. Автоматизированное оборудование и технология для пайки волноводных трактов космических аппаратов / С. К. Злобин, М. М. Михнев, В. Д. Лапте-нок и др. // Вестник СибГАУ 2014. № 4(56). С. 219-229.
17. Modeling of thermal processes in waveguide tracts induction soldering / A. V Murygin, V. S. Tynchenko, V. D. Laptenok et al. // IOP Conference Series: Materials Science and Engineering 173(1) 012026. 2017.
©Bocharova O. A., Muiygin A. V, Bocharov А. N., Zaitsev R. V., 2020
Bocharova Olesya Andreevna - senior lecturer of the department of Information and control systems; Reshetnev Siberian State University of Science and Technology. E-mail: shyx_89@mail.ru.
Murygin Aleksandr Vladimirovich - Dr. Sc., professor, head of the department of Information and control systems; Reshetnev Siberian State University of Science and Technology. E-mail: avm514@mail.ru.
Bocharov Aleksei Nikolaevich - Cand. Sc., docent of the department of Information and control systems; Reshetnev Siberian State University of Science and Technology. E-mail: sibalexbo@gmail.com.
Zaitsev Roman Viktorovich - a post-graduate student; Reshetnev Siberian State University of Science and Technology. E-mail: spoon27@yandex.ru.
Бочарова Олеся Андреевна - старший преподаватель кафедры информационно-управляющих систем; Сибирский государственный университет науки и технологий имени академика М. Ф. Решетнева. E-mail: shyx_89@mail.ru.
Мурыгин Александр Владимирович - доктор технических наук, профессор, заведующий кафедрой информационно-управляющих систем; Сибирский государственный университет науки и технологий имени академика М. Ф. Решетнева. E-mail: avm514@mail.ru.
Бочаров Алексей Николаевич - кандидат технических наук, доцент, доцент кафедры информационно-управляющих систем; Сибирский государственный университет науки и технологий имени академика М. Ф. Решетнева. E-mail: sibalexbo@gmail.com.
Зайцев Роман Викторович - аспирант; Сибирский государственный университет науки и технологий имени академика М. Ф. Решетнева. E-mail: spoon27@yandex.ru.
