Zablodskiy N. N., doctor of Sciences, prof.
Gritsyuk V. Yu., Ph.D., head teacher, Zablodskaya T. P. head teacher (DonSTU, Alchevsk, Ukraine)
MODELLING OF THREE DIMENSIONAL ELECTROMAGNETIC FIELD DISTRIBUTION IN THE ACTIVE PART ELECTROMECHANICAL DISINTEGRATOR
Developed three-dimensional finite element mathematical model of electromechanical disintegrator. Obtained the results of the spatial distribution of electromagnetic field in its active portion. Identified the features distribution of eddy currents and electromagnetic forces in the working conductive camera of disintegrator.
Key words: electromechanical disintegrator, three-dimensional finite element model, eddy currents, electromagnetic forces.
_ISSN 2077-1738. 36ipnuK HayKOBUx npa^ ffonffTY. 2014. № 1 (42)
_E.tektpotexhika. paaiotexhika
UDC 621. 313: 621.318.122
The problem and its relation to scientific and practical tasks. One of the ways to implement the various processes fine and superfine grinding, homogeneous mixing of liquid and solid powder materials (preparation of emulsions, suspensions, etc.), acceleration of chemical reactions is the use of a fundamentally new class of electromechanical energy converters — electromechanical disintegrator (EMD) multifactor impact [1, 2, 3]. Work area EMD is an area in which are counter running electromagnetic fields affecting workers ferromagnetic body (WFB), which in turn become magnetic dipoles and interact with the resultant field — the original source of energy. This results in a number of effects which, in addition to mechanical and thermal effects of the WFB, acting directly on the material, changing its physic and chemical properties.
With using the two-dimensional mathematical modelling, given in [4, 5], determine the complexity of the electromagnetic field distribution in the active part of EMD taking into account all the spatial components of the eddy current in the conductive working camera, is almost impossible. This problem makes it necessary to develop a three-dimensional mathematical model of EMD.
Statement of the problem. The objective of this work is to create a three-dimensional finite element model of EMD and determination based its features of electromagnetic field distribution in the active part.
Presentation of the material and its results. The general form of the nonlinear differential equation of the electromagnetic field in the partial derivatives to the vector magnetic potential A can be represented as
rot(- rot A) - y ^ - y(v x rot A) = -Jext, ^ ot
where |i = |i0 — absolute permeability; Y — specific conductivity; v — vector velocity of the electroconductive medium relative to the magnetic field source; Jext — density
extraneous current.
Algorithm for the numerical calculation of electromagnetic fields EMD, the physical properties of computational domains, as well as boundary conditions similar to those considered in [4]. Material working camera — nonmagnetic stainless steel with the conductivity — 1,1107 S/m.
To reduce the design time calculations and size of the needed amount of hardware resources, three-dimensional electromagnetic field calculation is made for a bipolar version of EMD. To create a three-dimensional finite element model was used software package Comsol Multiphysics 3.5a. On figures 1 and 2 presents the three-dimensional geometry and finite-element mesh of model, respectively.
It should be noted that the use of a PC based 4-core processor, AMD Phenom II X4 955 yielded good results for the three-
ISSN2077-1738. 36ipnuK uayKoeux npa^b flouflTY. 2014. № 1 (42)
E^EKTPQTEXHIKA Pajiqtexhika
dimensional model with a high degree of sampling grid elements. Thus, for comparison on figures 3 and 4 show the distribution curves of the normal component of the magnetic induction in the middle of the air gap in the yoke and along one pole EMD, respectively, for the three- and two-dimensional electromagnetic calculation options.
Figure 1 — Three-dimensional geometry
Figure 3 — Distribution of the normal component of the magnetic induction in the middle of the air gap (top) and in the yoke for three-dimensional calculation of EMD
Figure 2 — Finite-element mesh of model
Figure 5 shows the distribution of the *<>•«j \ \ V j
normal component of the magnetic induction 006 j \ \ \t
in the active part of the EMD in several cross 0Mr......................................1
sections. The highest value of magnetic in- 0 02° »■« A°-gth »■« »■« »■» in-
duction is observed in narrow areas of magnetic teeth — 2 T, value insignificantly de- Figure 4 — Distribution of the normal component creases as we move towards the middle of the of the magnetic induction in the middle of the air active part. gap (top) and in the yoke for two-dimensional
calculation of EMD
ISSN 2077-1738. 36ipnuK HayKOBUx npa^ ffonffTY. 2014. № 1 (42)
E.tehtpotexhika. Pa^iotexhika
Comsol Multiphysics allows defining the resultant force which acts on a conductive working camera by integrating elementary forces by volume, as well as efforts to separate spatial components, we are interested most. To deformation of the working camera EMD causes efforts normal to the surface of the camera wall. Numerical calculations showed on the pole pitch in a conductive camera distributed force action that reaches 25 newtons.
Figure 5 — Distribution of the normal component of the magnetic induction in the active part of the EMD in several cross sections
Three-dimensional graph of the distribution of the magnetic induction in the air gap EMD is shown on figure 6. Distribution of the vector magnetic potential in the active part of EMD in a three-dimensional graph shown on figure 7.
Calculation of the eddy currents in the conductive camera EMD work is shown on figure 8, distribution of instantaneous z-component of the current density (left), as well as a line current (right) corresponding to the classic "contours" of eddy currents.
A consequence of the interaction between the running magnetic fields of inductor and eddy currents flowing in the conductive camera walls (adjacent to the inductor), is the occurrence of the characteristic deflection of the camera walls within each pole pitch under the action of electromagnetic forces.
Distribution curve deflection depth along the length of the wall of the working camera EMD, built on the results of experimental measurements [5] and the nature of the distribution of electromagnetic forces obtained by solving the three-dimensional field problem by finite element method (figure 9), shown qualitative agreement and confirmed the presence of localization of existing efforts within the pole dividing inductors.
Figure 6 — Three-dimensional graph of the distribution of the magnetic induction in the air gap EMD
Figure 7 — Distribution of the vector magnetic potential in the active part of EMD in a three-dimensional graph
ISSN 2077-1738. 36ipuuK uayKoeux npaub ffouffTY. 2014. № 1 (42)
EflEKTPQTEXHIKA. Paaiqtexhika
Figure 8 — Distribution of eddy current in the working conductive camera EMD
Figure 9 — Three-dimensional distribution of electromagnetic forces in the working conductive camera EMD
ISSN 2077-1738. Збгрник наукових праць ДонДТУ. 2014. № 1 (42)
Електротехшка. Радютехшка
Given the pulsatory character of the action forces and impact load, to ensure the mechanical strength of the working camera, one can recommend its manufacturing of nonmagnetic material (stainless steel, titanium) with a minimum thickness of 4-5 mm.
Conclusions and directions for further research. Developed three-dimensional finite element mathematical model of electromechanical disintegrator. Obtained results of the
Reference List
spatial distribution of electromagnetic field in its active portion. The features of the distribution of eddy currents and electromagnetic forces in the working conductive camera disintegrator.
Further research should be aimed at a detailed study of the influence of the geometric dimensions of the active part of EMD on the power characteristics of the electromagnetic field in order to obtain needed operating properties.
1. Shinkarenko V. F. Evolutionary synthesis of new types of electromechanical energy converters technological purposes using models macroevolution / V. F. Shinkarenko, S. А. Bezsonov // Proceedings of the National technical university "KPI" - Kharkiv: NTU «KPI», 2001. - № 16. - P. 171-173.
2. Pat. 43635 Ukraine, МРК B01F 13/00. Electromechanical disintegrator / Shinkarenko V. F., Av-gustinovich А. А., Lisak V. V., Vahnovetska М. О. ; patent from the National technical university of Ukraine "KPI". - № u200902781, decl. 25.03.2009; publ. 25.08.2009, bull. № 16.
3. Filatov М. А. Assessing the possibility of using electromechanical disintegrator in coal-water fuel preparation technology / М. А. Filatov, V. Yu. Grytsyuk, V. V. Lisak // Collection of scientific works DonSTU. - Alchevsk: DonSTU. - 2012. - № 37. - P. 325-329.
4. Zablodskiy N. N. Modelling of electromagnetic field electromechanical disintegrator / N. N. Zablodskiy, М. А. Filatov, V. Yu. Grytsyuk // Collection of scientific works DonSTU. - Alchevsk: DonSTU. - 2013. - №. 39. - P. 221-226.
5. Zablodskiy N. N. Evaluation of electromagnetic forces in the working area of electromechanical disintegrator / N. N. Zablodskiy, V. Yu. Grytsyuk, S. V. Sheka // Collection of scientific works DonSTU. -Alchevsk: DonSTU. - 2013. - № 39. - P. 212-216.
Recommended for printing by the prof. Lushchik V. D.,
prof. Sivokobilenko V. F.
Paper received 27.02.14.
д.т.н. Заблодський М. М., к.т.н. Грицюк В. Ю., Заблодська Т. П. (ДонДТУ, м. Алчееськ, Украгна) МОДЕЛЮВАННЯ ТРИВИМ1РНОГО РОЗПОД1ЛУ ЕЛЕКТРОМАГН1ТНОГО ПОЛЯ В АКТИВН1Й ЧАСТИН1 ЕЛЕКТРОМЕХАН1ЧНОГО ДЕЗ1НТЕГРАТОРА
Розроблено mриеимiрну ктцеео-елементну математичну модель електромехатчного дезт-тегратора. Отримано результати простороеого розподшу електромагнтного поля е його ак-тиетй частит. Визначено особлиеостi розподшу еихроеих струмiе й електромагнтних зусиль е робочт електропроеiднiй камерi дезттегратора.
Ключовi слова: електромехатчний дезттегратор, триеимiрна ктцеео-елементна модель, еихроеi струми, електромагттт зусилля.
д.т.н. Заблодский Н. Н., к.т.н. Грицюк В. Ю., Заблодская Т. П. (ДонГТУ, г. Алчееск, Украина)
моделирование трехмерного распределения электромагнитного
ПОЛЯ В АКТИВНОЙ ЧАСТИ ЭЛЕКТРОМЕХАНИЧЕСКОГО ДЕЗИНТЕГРАТОРА
Разработана трехмерная конечно-элементная математическая модель электромеханического дезинтегратора. Получены результаты пространстеенного распределения электромагнитного поля е его актиеной части. Определены особенности распределения еихрееых токое и электромагнитных усилий ерабочей электропроеодящей камере дезинтегратора.
Ключевые слова: электромеханический дезинтегратор, трехмерная конечно-элементная модель, еихрееые токи, электромагнитные усилия.