CC BY
9
13. Filimonikhin, G. Investigation of the possibility of balancing aerodynamic imbalance of the impeller of the axial fan by correction of masses [Text] / G. Filimonikhin, L. Olijnichenko // Eastern-European Journal of Enterprise Technologies. - 2015. - Vol. 5, Issue 7 (77). - P. 30-35. doi: 10.15587/1729-4061.2015.51195
14. Goncharov, V. An increase of the balancing capacity of ball or roller-type auto-balancers with reduction of time of achieving auto-balancing [Text] / V. Goncharov, G. Filimonikhin, A. Nevdakha, V. Pirogov // Eastern-European Journal of Enterprise Technologies. - 2017. - Vol. 1, Issue 7 (85). - P. 15-24. doi: 10.15587/1729-4061.2017.92834
Розглянуто задачу визначення межi переходу руху з колових на квазiпараболiчнi траектори час-тинок зернистого заповнення цилшдричног каме-ри, що обертаеться навколо горизонтальног оы. Поведтку заповнення описано за допомогою осеред-нених величин. Застосовано пластичну реологiч-ну модель зернистого середовища. На основi розра-хунку напружено-деформованого стану визначено положення поверхш переходу, як лти ковзання кар-тини руху. Выявлено залежностi межi вгд швидтсно-го, геометричних i реологiчних параметрiв системи
Ключовi слова: зернисте заповнення, обертова камера, квазтвердотшьна зона, поверхня руйну-
вання, л^я ковзання
□-□
Рассмотрена задача определения границы перехода движения с круговых на квазипараболические траектории частиц зернистого заполнения цилиндрической камеры, вращающейся вокруг горизонтальной оси. Поведение заполнения описано с помощью осредненных величин. Применена пластическая реологическая модель зернистой среды. На основе расчета напряженно-деформированного состояния определено положение поверхности перехода, как линии скольжения картины движения. Выявлены зависимости границы от скоростного, геометрических и реологических параметров системы
Ключевые слова: зернистое заполнение, вращающаяся камера, квазитвердотельная зона, поверхность разрушения, линия скольжения -□ □-
UDC 621.926.5 : 539.215
|DOI: 10.15587/1729-4061.2017.96447|
MODELING OF
fracture surface
OF THE QuASI SOLID-BODY ZONE OF MOTION OF THE GRANuLAR FILL IN A ROTATING CHAMBER
Yu. Naumenko
Doctor of Technical Sciences, Associate Professor Department of construction, road, reclamation, agricultural machines and equipment National University of Water and Environmental Engineering Soborna str., 11, Rivne, Ukraine, 33028 E-mail: [email protected]
1. Introduction
The main equipment in the large-tonnage processing of granular materials in chemical, mining, metallurgical, building materials industry, power engineering and many other sectors of industrial production is still traditional machines of the drum type [1]. The main varieties of such a wide class of machines are rotary kilns [2] and coolers [3], drying drums [4], ball [5] and tube [6] mills.
A total number of different materials, which annually undergo processing in such machines in the world, are several tens of billions of tons. A widespread use of machines of this class is due to the possibility of their large performance capacity, reliability in operation, ease of use, versatility and high efficiency.
However, a paradoxical feature of this technological equipment is the combination of the utmost simplicity of
design and behavior of the treated medium, which is extremely difficult to describe. Significant difficulties in predicting the implementation of processes in drum machines considerably limit a possible technological scope of their application.
Given that the implementation of working processes in the drum-type machines is predetermined by the character of motion mode of the treated granular material in a chamber, the task on establishing the patterns of change in the parameters of this mode appears to be rather relevant.
2. Literature review and problem statement
Motion modes of granular fill in the rotating chambers significantly affect the realization of technological processes and energy efficiency of the drum-type machines drive [7].
©
Modeling of such modes is of interest in the examination of various rotating systems [8].
A practical relevance of the task for predicting the working processes of these machines constantly attracts increased research attention to the description of behavior of the treated granular media. An enhanced complexity of this problem forces to improve traditional techniques as well as to apply new theoretical and experimental methods.
Many attempts were made to numerically solve the problem on determining the motion of granular fill in a rotating cylindrical chamber. The modeling of motion modes of the fill in a wide range of variations in the degree of filling and rotation velocity of a chamber, using the method of discrete elements, was performed in [9]. Parameters of dense two-dimensional established granular flow in a chamber were described in [10] based on the Euler modeling with the Boussinesq approximation. Article [11] presents results of numerical modeling of a three-dimensional non-established flow of active layer of granular particles under the rolling over mode, which was based on the theory of Euler biphase flow and kinetic theory of granular flow. An intermediate, between the cascade and cataract modes, granular flow at the values of filling degree of a chamber at 0.5 and the Froude number of 0.56, employing the method of discrete elements, was examined in [12]. Results of testing and verifying the source data, with the application of the method of discrete elements, for modeling the motion of granular fill in a rotating chamber, are presented in [13]. Two motion modes of granular fill with a flat and a curved free surface were explored in [14] based on the theory of deep integration that was applied for non-dense flows. In [15], using the method of discrete elements, authors obtained characteristics of the distribution of active and passive zones of a three-dimensional motion of granular fill in a chamber.
However, the initial conditions of the examined problem are a priori uncertain while the boundary conditions are of non-physical nature. This causes a significant constraint for the accuracy of numerical calculations, the results of which do not meet practical requirements.
The motion of granular loading in a rotating chamber was also investigated by various experimental methods. Dynamic characteristics and velocity fields of granular fill were determined, based on the method of the visualization of flow pictures in the cross-section of a chamber, in [16].
Limited capabilities for a visual analysis of behavior of the fill predetermined a widespread use of positron-emission methods of tomographic analysis. Article [17] examined dynamics of non-spherical particles in a rotating chamber using the method of radioactive particles. The application of such technique to establish trajectories of particles under the rolling over mode is described in [18].
The application of positron-emission observation of particles to assess the impact of variations in the speed of chamber rotation on the flow mode of spherical particles is described in [19]. Testing of a new high-precision method of positron-emission particle observation for determining the spatial and temporal kinematic characteristics of a granular flow is covered in [20]. The behavior of active and passive layers of cylindrical and spherical particles in a rotating chamber, using the method of multiple radioactive tracking of particles, was studied in [21]. [22] describes characteristics of the flow of granular loading in drum mills near free surface at big Froude numbers for the cascade and cataract motion modes that were measured by using the positron-emission particle tracking.
A comparison of results in determining the parameters of motion of granular fill in a rotating chamber, which was obtained numerically by the method of discrete elements and experimentally by the method of tracking the behavior of radioactive particles, is presented in [23].
A significant progress has been observed lately in the application of tomographic methods to explore the behavior of filling the closed chambers. However, considerable complexity in their implementation and limited measurement resolution reduce reliability and accuracy of the results obtained.
In this regard, [24] evaluated the complexity and significance of adequate determining of rheological properties of granular media for the case of solving a problem on defining the motion of filling a rotating chamber by analytical method. It was indicated that the rheological characteristics of such media significantly vary depending on the flow type.
The character of motion modes of granular fill in a rotating chamber largely depends on the position and shape of surface in the transition zone from a solid-body motion to the zone of falling and shattering. Such a surface is formed on top of the chamber and defines dynamic parameters of active movable part of the fill. In order to solve the problem on determining such a characteristic element in the flow of array of fill, analytical research methods were mostly applied.
At the same time, a traditional hypothesis on the realization of operational processes when processing granular materials in the drum machines [25] is based on the concept of a separate element of the fill in a chamber, isolated from the surrounding medium. Such element moves under the action only of gravity force G, centrifugal inertia C and reaction of the limiting surface N and T (Fig. 1).
Fig. 1. Trajectory of motion of a granular fill's element in a rotating chamber and the forces applied to it by the traditional hypothesis
According to this hypothesis, the interaction between elements is disregarded. It is considered that elements of the fill are captured by the surface of a rotating chamber until the moment when the gravity forces exceed the centrifugal forces, after which these elements start to move by parabola. Trajectories of elements' motion do not intersect. It is assumed that the fill motion mode is a two-phase «waterfall» with a solid-body zone and the toss and drop zone in the cross section of a chamber. The locus of transition of the fill's elements from the circular to the parabolic trajectories (Fig. 2) is assigned by an auxiliary circle of diameter a
g
(1)
where g is the gravitational acceleration, ra is the angular rotation velocity of a chamber.
a =
2
Fig. 2. Locus of transition points of the granular fill's elements in a rotating chamber from circular to parabolic trajectories by the traditional hypothesis
However, the real motion mode of the fill is a three phase one when the quasi solid-body and non-free falling zones in the cross-section of a rotating chamber are complemented by a zone of the shifting layer. The main process of processing the granular materials is implemented in this particular zone.
Meanwhile, the problem on defining the surface of destruction in the solid-body zone, based on the equilibrium of a separate element in the granular fill of a rotating chamber, was considered in [26]. The stressed-strained state of the medium was not examined. An equation for the transition surface was obtained in a simplified statement depending on a limited number of parameters.
[27-31] analytically examined in a static statement the sliding surfaces when considering the carrying capacity of slopes in granular materials, mainly in the application to soils.
Ref. [27] considered a problem on the analytical calculation of non-movable soil slope stability on the sliding surfaces of a simplified circular-cylindrical shape. A method of hypothetical circular sliding surfaces with regard to the soil non-uniformity was applied in [28]. In [29], based on the Mohr theory, authors examined positions of sliplines in a perfect plastic medium. The stability of slopes in granular material when reaching the boundary condition in a static statement was investigated in [30]. [31] explored the boundary stressed-strained state of slopes' stability in a granular medium.
Paper [32], based on the model of pre-boundary elasto-plastic deformation, defined the stressed-strained state of a granular medium only at slow rotation of the chamber.
Article [33] analytically in a plastic statement attempted to determine the stressed-strained state of the fill in the form of simplified picture of the motion and velocity field of a shifting layer at low rotation rate.
However, the overwhelming computational difficulties and low reliability of the hardware control restrict effectiveness of the known methods for determining the surface of mode transition in the fill motion. That is why the obtained results of numerical calculations and experiments approach the real modes of motion of the examined medium only by qualitative characteristics and external attributes. By the quantitative indicators, they are substantially divergent.
In view of the aforementioned, an analytic modeling of the required boundary for the case of non-low rotation velocity of a chamber appears to be quite an important applied problem. Of particular interest is taking account of variations
in geometric and rheological characteristics within a wide range of variations in the parameters of the system.
3. The aim and tasks of the study
The aim of present work is to create a mathematical model for the behavior of a granular fill in the transition from the circular, during the ascent, to the quasi-parabolic, when lowering, trajectories of sustainable motion in the cross-section of a cylindrical chamber that permanently revolves around a horizontal axis. This will allow us to determine dynamic parameters in the motion of a sliding layer of the fill and to predict effectiveness in the implementation of technological processing of a granular material in a rotating chamber.
To accomplish the set aim, the following tasks were to be solved:
- to perform analytical modeling of stress field in the mass of a granular fill in the cross-section of a chamber;
- to establish conditions for disequilibrium in filling the chamber;
- to define the shape and position of the surface of destruction of the quasi solid-body zone in the fill motion;
- to establish effect of the system's parameters on the boundary of transition in the motion zones.
4. Procedure of examining the stressed-strained state of the fill in a chamber
4. 1. General conceptual approach to the research technique
A granular fill in a rotating chamber was considered to be a multi-phase polydispersed system with a non-uniform dispersed phase. It was thought that the size of inhomogeneties is much less than at the distance at which the macroscopic or averaged parameters of multiphase medium of the fill substantially changed. A realization of the adopted assumption made it possible to describe a change in the averaged parameters of a multiphase system within representations and based on conventional equations and methods of mechanics of the entire medium. Granular fill of a rotating chamber was considered as a whole medium with parameters, which are averaged by volume and uniformly distributed in the space that a multi-phase medium occupies. Mathematical description of the fill motion was performed by using averaged magnitudes. A method of calculating the stressed-strained state of the fill was employed. The fill flow was thought to be slow with predominance of the forces of plastic origin over the viscous inertial forces. A plastic rheological model of granular medium of the fill was accepted.
4. 2. Applied methods for modeling the behavior of granular medium of the fill
Subsequent numerical solution of the problem on determining the limit of strength of a solid-body array, at sustainable motion of the granular fill, averaged by volume, in a rotating chamber, is to be found in a plastic statement. Experimental data revealed a relatively large value of the angle of repose in the motion along free surface at a moderate rotation velocity of the chamber. This made it possible to accept a pseudo plastic rheological model of a medium. This model is a particular case in the manifestation of viscoplastic properties of a rheologically-complex fill.
Shift resistance in the point of filling consists of the resistance caused by internal friction and adhesion and is expressed by a dependence that occurs at disequilibrium [30, 34]
|tn| = Sn tg j + k,
where on and tn are the normal and tangent components of pressure, n is a normal to the sliding surface, 0 is the angle of internal friction, k is the adhesion coefficient.
A two-dimensional stressed state in point P of a granular fill is characterized by a general scheme (Fig. 3).
Fig. 3. General scheme of a flat stressed state in the point of filling
In this point, real tension p and resulting tension p' act. Real tension p forms angle 8 with normal n and has normal on and tangent tn components. Resulting tension p' forms angle 8' with normal n and has normal on+H and tangent tn components, where H = k-tgj is the coefficient of temporary resistance to the comprehensive uniform stretching.
Slipping as a result of shift of the medium along the examined elementary area will not take place while maintaining the conditions of boundary equilibrium of the fill
|xn| <(on + H)tg j,
in this case, on>-H.
A visual representation of the stressed state in the point of granular fill is defined by a stress diagram (Fig. 4), where
s=(amax+omin)/2; t=(amax-omin)/2; omax and omin are the main normal stresses; l, || and v are the angles of direction cosines between the normal to the area and main axes. We can write
sin A
° P sin(A-kf80'
(1 -kf )n k f A-8'
1=--I——+---+ mn,
(2)
Fig. 4. Stress diagram of a flat stressed state in the point of filling
A rectangular coordinate system is used with a highlighted element with coordinates x, y and dimensions dx and dy (Fig. 5) where ox, oy and txy = tyx are three components of a pressure tensor.
Fig. 5. Schematic of a flat stressed-strained state of the fill medium
In the zone of boundary equilibrium, each point is crossed by two lines of sliding that comprise a system of two families (Fig. 6), where Q is the angle between direction omax and the x axis, Q ± e are the inclination angles of sliplines to the x axis, 1 and 2 are the slip lines.
Fig. 6. Directions of slip lines in the zone of flat boundary equilibrium of the fill
where kf = ±1, m is any integer, sinA = sind'/sinj (|A| < n/2). Two signs kf = -1 and kf = + 1 that are included in these equations correspond to the two stressed states - minimum and maximum.
Thus, in the zone of boundary equilibrium, through each point in the xy plane (Fig. 6), traces of sliding surfaces form two slip lines. The entire zone of the boundary equilibrium is covered by the two families of slip lines. However, only one
of the slip lines is active or is a break line. Position of the active family is defined by the direction along which a breach in boundary equilibrium occurs. This direction depends to a large degree on the type and initial conditions of the medium loading.
5. Results of modeling the fracture surface of a solid-body zone of the fill in a chamber
Subsequently, we consider the boundary equilibrium of an array of loose granular fill of a cylindrical chamber with radius R that rotates around a horizontal axis with constant angular velocity ra. A flat free surface of the array, before the onset of shift, is limited by the x axis and tilted to the horizontal at angle a (Fig. 7). At the array's surface, we evenly distributed at constant angle 8 to the normal reduced pressure p, caused by the action of applied centrifugal forces of inertia. Positive value 8, based from the normal to the surface of the array, is assigned in a counterclockwise direction (Fig. 8). Expression for angle 8 can be represented in the form
a.
8 = y-P,
where al is the lifting angle of the fill in a rotating chamber; b is the half of the central angle of a segment cross-section of the fill at rest, which is determined from equation
2b- sin 2b = 2nk,
where k is the degree of filling the chamber with a granular material.
Expression for angle a takes the form
Coriolis forces of inertia. Additional action of variational forces is caused by the separation of elements from the array with subsequent non-free fall. Additional action of Coriolis forces is caused by relative motion of the slip layer of the fill.
Fig. 8. Scheme of angles for the calculation of slip lines of the fill
The ratio of pressure p and specific bulk weight of the filler Y can be approximately represented by expression
P / \
— = —(CD2 - CG2). Y 2g '
(3)
Based on the solution of oblique-angled triangle AGC (Fig. 8) with angles ZGAC = n/ 2 - b, ZCGA = n/ 2 + b - a^2 and ZACG = a i/2, it is possible to write
CG = R-
cos b
b-
a = a i-ß + -.
Fig. 7. General estimated scheme of slip lines of the fill
Then relation (3) takes the form
p ra2R2 Y
2g
1 -
cos2 b
b-
Normal and tangent components of pressure along the x axis can be expressed in the form
O y + H = pcos 5, t xy = psin 5.
(4)
In a general case, stress field in the array is determined by a system of equations of flat boundary equilibrium of a granular medium
Fx-g x Y
F- g Fy Y
3Ox dt xy
3x 3y
3Oy 3tyx
= 0,
= 0,
3y 3x
(ox - oy )2 + 4txy = (ox + O y + 2k • ctg j)2 sin2 j,
The value of reduced pressure p is taken for the central part of the array - point S (Fig. 8) that approaches point D at the surface of a chamber. The magnitude of this pressure takes into account the effect of centrifugal, variational and
where Fx and Fy are the projections of mass forces. The first and the second equations of the system are differential equations of equilibrium, the third equation is the condition for boundary equilibrium.
A ■
2
Under condition of the x-axis inclination to the horizontal at angle a:
Fx=ysina, Fy=ycosa.
Then the stress field that occurs in the examined array does not depend on x and is described by the differential equations of equilibrium
"dy" dy
= y sin a,
= y cos a.
and the width of band y0 < y < y* where 0 < 8y < p is equal to p sinpsin(a-8)
b = (y.- yo ) =--:--r.
v ' y sin a sin (a-p)
In order to determine the slip lines of the examined stress field, it is necessary to apply differential equation
(5)
After integration of equations (5) with regard to (4), it is possible to obtain
I = tg (Q * ^
(7)
y is possible to express through Q in (7) by using relation
sin 2Q
y - a = (y - yo ) = bsin (a-p)—-—
sin a + sin p sin (2Q + a)
oy + H = pcos 8 +Yycos a, t „ = psin 8 +Yysin a.
(6)
Based on (6), the reduced stress that acts on the inclined lines, parallel to the x axis, is determined by magnitudes
tg8y =
psin 8 +Yysin a
pcos 8 + Yycos a' py = pcos (8-8 y) +Yycos (a-8y),
where py and 8y are the stress and its inclination angle to the normals to the slanted straight line with ordinate y. By analogy to (2), it is possible to write
sin A„
° pysin (Ay -kf8y )
(1 -kf )n k fAy-8y
Q = -—+-2—_ + mn,
where
sinAy = sin8y/sinp (|Ay| < n/2).
Then the y coordinate is expressed through dy or Q in the following way
psin (8-8y) psin 8-sin p sin (2Q + 8)
y Y sin (a-8y) Y sin a-sin p sin (2Q + a)
Straight line, along which txy=0 and which satisfies condition 8y=0, has ordinate (Fig. 7)
p sin8
a = (y0 ) =---:-.
v ' y sin a
If a < p, then the boundary equilibrium is possible in all half-plane 0 < y < ~ since inequality |8y| < p is fulfilled automatically. If a > p, then the boundary equilibrium is possible only in some plane 0 < y < y», in which inequality |8y| < p is satisfied.
The ordinate that satisfies condition 8y=p is (Fig. 7) psin (p-8)
c = (y.) = ——i-t ,
v ' Y sin (a-p)
Then it is possible to obtain after transforms
dx sin2Q± cos p
-t—= -2bsin a sin (a-p)-j. (8)
dQ [sin a + sin p sin (2Q + a)^
Equation (8) can be integrated in elementary functions
x (Q) =
sin 2Q ± cos p
= -2bsin a sin (a-p)-^ dQ+const.
[sin a + sin p sin (2Q+a)
Finally, the coordinates of required slip line AB (Fig. 7) can be represented in the parametric form
sin2Q
y = bsin (a-p)- . ,
sin a + sin p sin (2Q + a)
2bsin a sin (a-p)
x1,2 = : 2 : 2 x sin2 a- sin2 p
[± cos p(+ sin p) - sin a] cos (2Q+a) 2 [sin a+(+ sin p) sin (2Q + a)]
± cos p sin a - (+ sin p)
+ I = ^
^/sin2 a- sin2 p
x arctg
sin a- tg
Q+
sin p)
^/sin2 a- sin2 p
[± cos p(+ sin p) - sin a] cos a 2 [sin a+(+sin p) sin a]
± cos p sin a - (+ sin p)
^/sin2 a- sin2 p
x arctg
sin a-tg
2
V_/
-(+ sin p)
^/sin2 a- sin2 p
(9)
(10)
A parametric form of the equations derived allows us to define and analyze the coordinates of a slip line, depending on the variation of parameters of the examined system.
6. Discussion of results of examining the impact of parameters of a filled rotating chamber on the fracture surface of a solid-body zone of the fill
Mathematical description of the transition boundary position in the motion of granular fill in a rotating chamber, from the circular to the quasi parabolic trajectories, was performed according to a classical scheme with the use of magnitudes averaged by volume. We made an assumption about the slow motion of a loose granular medium with the predominance of forces, plastic in origin, over the viscous inertial forces. This made it possible to obtain equations that approximately establish dependences of the required limit coordinates in the parametric form on the series of parameters of the examined system:
- geometric: radius R of the chamber and the degree of its filling k;
- kinematic: angular rotation velocity of the chamber ra;
- inertial: specific bulk weight of the fill y;
- rheological: the angle of internal friction of the fill j and its lifting angle in arotating chamber al.
The advantage of the proposed approach, in comparison with the traditional hypothesis (1), is the possibility to derive such limit from parameters R, k, ra, y, j and a i. According to (9) and (10), the shape of such a surface is curvilinear non-cylindrical. Instead, according to (1), this limit is of cylindrical shape with a diameter that depends only on ra.
An analysis of (9) and (10) reveals that at increase in the values of parameters ra, y, j and al and decrease in the values of R and k, location of the transition limit moves away from the free surface of the fill and approaches surface of the chamber.
The disadvantages of the developed algorithm of numerical calculation of the limit include the need for a preliminary experimental determining of ascent angle of the fill in a rotating chamber al. For this purpose, it is
possible to use, for example, the method of visual analysis of pictures of motion in the cross section of a chamber.
In the future, it is expedient, with regard to the established characteristics of the limit of breaking the quasi solid-body zone of the fill in a rotating chamber, to investigate, based on analysis of the stressed-strained state, parameters of the shear flow of the fill near a free surface.
7. Conclusions
1. We formalized a stress field in the mass of filler in the cross-section of a rotating chamber by using a system of differential equations of flat boundary equilibrium of a granular medium. Equations are obtained that approximately establish positions of the transition limit of particles in the fill from the circular to the quasi parabolic trajectories, depending on the rotation velocity, geometric and rheological parameters of the system.
2. We received a condition for boundary equilibrium in the fill. It is demonstrated that disequilibrium is accompanied by slip layers. The patterns were discovered in the transition of quasi solid-body motion zone of the fill to the zone of nonfree fall and shear flow. It is found that such a transition is accompanied by the destruction of a solid-body zone at the border, which is a slip line of the picture of motion of plastic medium in the cross-section of a chamber.
3. We defined a position and shape of slip lines in the stress field of the fill. It is established that the breaking limit of a solid-body zone in the fill is a curvilinear non-cylindrical surface. It was found that the position of the extreme top line in the transition coincides with the surface of the chamber and is determined by ascent angle of the fill.
4. It is established that the shape and position of the transition limit depend on the magnitude of radius R, degree of filling k and angular rotation velocity of the chamber ra, specific bulk weight y and the angle of internal friction j of the fill and its ascent angle in a rotating chamber al. It was found that at increase in the values of parameters ra, j, y and al and decrease in R and k, position of the transition limit approaches a surface of the chamber.
References
1. Pershyn, V. F. Pererabotka sypuchykh materyalov v mashynakh barabannoho typa [Text] / V. F. Pershyn, V. H. Odnolko, S. V. Per-shyna. - Moscow: Mashynostroenye, 2009. - 220 p.
2. HOST 27120-86. Pechy khymycheskykh proyzvodstv s vrashchaiushchymysia barabanamy obshcheho naznachenyia. Osnovnye parametry y razmery [Text]. - Moscow: Yzd-vo standartov, 1992. - 17 p.
3. HOST 11875-86. Apparaty teploobmennye s vrashchaiushchymysia barabanamy obshcheho naznachenyia. Kholodylnyky. Osnovnye parametry y razmeruy [Text]. - Moscow: Yzd-vo standartov, 1988. - 8 p.
4. HOST 27134-86. Apparaty sushylnye s vrashchaiushchymysia barabanamy. Osnovnye parametry y razmery [Text]. - Moscow: Yzd-vo standartov, 2002. - 4 p.
5. HOST 10141-91. Melnytsy sterzhnevye y sharovye. Obshchye tekhnycheskye trebovanyia [Text]. - Moscow: Yzd-vo standartov, 1991. - 19 p.
6. HOST 12367-85. Melnytsy trubnye pomolnykh ahrehatov. Obshchye tekhnycheskye uslovyia [Text]. - Moscow: Yzd-vo standar-tov, 1986. - 13 p.
7. Naumenko, Yu. V. The antitorque moment in a partially filled horizontal cylinder [Text] / Yu. V. Naumenko // Theoretical Foundations of Chemical Engineering. - 1999. - Vol. 33, Issue 1. - P. 91-95.
8. Naumenko, Yu. V. Determination of rational rotation speeds of horizontal drum machines [Text] / Yu. V. Naumenko // Metallurgical and Mining Industry. - 2000. - Issue 5. - P. 89-92.
9. Arntz, M. M. H. D. Granular mixing and segregation in a horizontal rotating drum: A simulation study on the impact of rotational speed and fill level [Text] / M. M. H. D. Arntz, W. K. den Otter, W. J. Briels, P. J. T. Bussmann, H. H. Beeftink, R. M. Boom // AIChE Journal. - 2008. - Vol. 54, Issue 12. - P. 3133-3146. doi: 10.1002/aic.11622
10. Bonamy, D. Euler-like modelling of dense granular flows: application to a rotating drum [Text] / D. Bonamy, P.-H. Chavanis, P.-P. Cortet, F. Daviaud, B. Dubrulle, M. Renouf // The European Physical Journal B. - 2009. - Vol. 68, Issue 4. - P. 619-627. doi: 10.1140/epjb/e2009-00123-6
11. Yin, H. Numerical simulation of three-dimensional unsteady granular flows in rotary kiln [Text] / H. Yin, M. Zhang, H. Liu // Powder Technology. - 2014. - Vol. 253. - P. 138-145. doi: 10.1016/j.powtec.2013.10.044
12. Arntz, M. M. H. D. Segregation of granular particles by mass, radius, and density in a horizontal rotating drum [Text] / M. M. H. D. Arntz, H. H. Beeftink, W. K. den Otter, W. J. Briels, R. M. Boom // AIChE Journal. - 2013. - Vol. 60, Issue 1. -P. 50-59. doi: 10.1002/aic.14241
13. Marigo, M. Discrete Element Method (DEM) for Industrial Applications: Comments on Calibration and Validation for the Modelling of Cylindrical Pellets [Text] / M. Marigo, E. H. Stitt // KONA Powder and Particle Journal. - 2015. - Vol. 32, Issue 0. -P. 236-252. doi: 10.14356/kona.2015016
14. Hung, C.-Y. Granular flow regimes in rotating drums from depth-integrated theory [Text] / C.-Y. Hung, C. P. Stark, H. Capart // Physical Review E. - 2016. - Vol. 93, Issue 3. doi: 10.1103/physreve.93.030902
15. Yang, S. DEM study of granular flow characteristics in the active and passive regions of a three-dimensional rotating drum [Text] / S. Yang, A. Cahyadi, J. Wang, J. W. Chew // AIChE Journal. - 2016. - Vol. 62, Issue 11. - P. 3874-3888. doi: 10.1002/aic.15315
16. Chou, H.-T. Cross-sectional and axial flow characteristics of dry granular material in rotating drums [Text] / H.-T. Chou, C.-F. Lee // Granular Matter. - 2008. - Vol. 11, Issue 1. - P. 13-32. doi: 10.1007/s10035-008-0118-y
17. Dube, O. Dynamics of non-spherical particles in a rotating drum [Text] / O. Dube, E. Alizadeh, J. Chaouki, F. Bertrand // Chemical Engineering Science. - 2013. - Vol. 101. - P. 486-502. doi: 10.1016/j.ces.2013.07.011
18. Alizadeh, E. Characterization of Mixing and Size Segregation in a Rotating Drum by a Particle Tracking Method [Text] / E. Alizadeh, O. Dube, F. Bertrand, J. Chaouki // AIChE Journal. - 2013. - Vol. 59, Issue 6. - P. 1894-1905. doi: 10.1002/aic.13982
19. Morrison, A. J. The shape and behaviour of a granular bed in a rotating drum using Eulerian flow fields obtained from PEPT [Text] / A. J. Morrison, I. Govender, A. N. Mainza, D. J. Parker // Chemical Engineering Science. - 2016. - Vol. 152. - P. 186-198. doi: 10.1016/j.ces.2016.06.022
20. Tupper, G. B. Testing of a new dynamic Ergun equation for transport with positron emission particle tracking [Text] / G. B. Tup-per, I. Govender, D. N. De Klerk, M. C. Richter, A. N. Mainza // AIChE Journal. - 2015. - Vol. 62, Issue 3. - P. 939-946. doi: 10.1002/aic.15081
21. Rasouli, M. Investigating the dynamics of cylindrical particles in a rotating drum using multiple radioactive particle tracking [Text] / M. Rasouli, O. Dube, F. Bertrand, J. Chaouki // AIChE Journal. - 2016. - Vol. 62, Issue 8. - P. 2622-2634. doi: 10.1002/aic.15235
22. Govender, I. A positron emission particle tracking investigation of the scaling law governing free surface flows in tumbling mills [Text] / I. Govender, M. C. Richter, A. N. Mainza, D. N. De Klerk // AIChE Journal. - 2016. - Vol. 63, Issue 3. - P. 903-913. doi: 10.1002/aic.15453
23. Alizadeh, E. Comparison of DEM results and Lagrangian experimental data for the flow and mixing of granules in a rotating drum [Text] / E. Alizadeh, F. Bertrand, J. Chaouki // AIChE Journal. - 2013. - Vol. 60, Issue 1. - P. 60-75. doi: 10.1002/aic.14259
24. Govender, I. Granular flows in rotating drums: A rheological perspective [Text] / I. Govender // Minerals Engineering. - 2016. -Vol. 92. - P. 168-175. doi: 10.1016/j.mineng.2016.03.021
25. Andreev, S. E. Droblenye, yzmelchenye y hrokhochenye poleznykh yskopaemykh [Text] / S. E. Andreev, V. A. Perov, V. V. Zvere-vych. - Moscow: Nedra, 1980. - 415 p.
26. Sverdlyk, H. Y. O strukture sechenyia materyala, peresypaiushchehosia vo vrashchaiushchemsia barabane [Text] / H. Y. Sverdlyk, H. H. Hryhorev // Yzv. vuzov. Chern. Metallurhyia. - 1977. - Issue 8. - P. 169-172.
27. Tsytovych, N. A. Mekhanyka hruntov [Text] / N. A. Tsytovych. - Moscow: Vyschaya shkola, 1983. - 288 p.
28. Kharr, M. E. Osnovy teoretycheskoi mekhanyky hruntov [Text] / M. E. Kharr. - Moscow: Hosstroiyzdat, 1971. - 320 p.
29. Naday, A. Plastychnost y razrushenye tverdykh tel. Vol. 2 [Text] / A. Naday. - Moscow: Yzd-vo Myr, 1969. - 864 p.
30. Sokolovskyi, V. V. Statyka sypuchei sredy [Text] / V. V. Sokolovskyi. - Moscow: Vyschaya shkola, 1960. - 241 p.
31. Kondratev, D. S. Upruhoplastycheskoe deformyrovanye y predelnoe ravnovesye sypuchykh sred [Text] / D. S. Kondratev, P. V. Ste-tsenko, Y. V. Shyrko // Prykladnaia matematyka y mekhanyka. - 2005. - Vol. 69, Issue 1. - P. 117-134.
32. Bushmanova, O. P. Dopredelnoe plastycheskoe deformyrovanye sypuchei sredy vo vrashchaiushchemsia barabane [Text] / O. P. Bush-manova, O. A. Revuzhenko // Fyzyko-tekhnycheskye problemy razrabotky poleznykh yskopaemykh. - 2004. - Issue 6. - P. 58-67.
33. Slanevskyi, A. V. Osnovy mekhanyky sypuchei sredy vo vrashchaiushchykhsia pechakh y melnytsakh [Text] / A. V. Slanevskyi. -Sankt-Peterburg, 1996. - 238 p.
34. Henyev, H. A. Voprosy dynamyky sypuchei sredy [Text] / H. A. Henyev. - Moscow Hosstroiyzdat, 1958. - 122 p.
