Mathematical model of complex control of the vibratory transportation and technological process

The vibratory transportation and technological process is a dynamically sensitive operation which includes physically different components: vibro-exciter, elastic system, working member (absolutely rigid or of finite rigidity) and various friable loads. Interaction of these components predetermines the behavior of the friable material on the surface of the working member (WM). At the same time, existing simple models or physical experiments cannot provide sufficient precision to adequately research the mentioned complex process. Therefore, it is necessary to develop a more precise mathematical model ensuring the study and revelation of the still hidden factors influencing the vibratory process. A new generalized dynamical spatial model of the loaded vibratory technologic machine (vibro-exciter, working member, load) developed on the basis of the systemic approach is presented in the work and a system of interconnected equations of movements of the constituent masses considering dynamical, geometrical and physical parameters, is obtained. The change of parameters is reflected on the variation of dynamical characteristics of the system that allows a thorough study of the technological process with the help of mathematical modeling. Using the presented model, it is possible to find the physical parameters and their combinations, realization of which will promote the improvement of the technological process. Some results of the modeling are presented. A new design of the vibro-exciter developed on the basis of the results of modeling is presented as well.


Introduction
The vibratory transportation and technologic machines are widely used in various spheres of industry for transportation of friable materials and individual parts, measured feeding and sorting and also for carrying out various technologic operations on them [1][2][3][4][5][6][7].
The vibratory conveyors are especially effective for transportation of powdery, explosive, chemically aggressive, abrasive, heated and other friable materials, for the displacement of which, other types of conveyors are less suitable [1,8].
The advantages of these machines are: -Possibility of hermetization of the transportation and technologic processes that decrease the pollution of the environment; -Elimination of losses and dirtying of the materials to be transported; -Minimal wear of the working member; -Compatibility of the process of transportation with scattering, heating, cooling, washing, mixing and other technologic operations.
In spite of the great number of researches in the sphere of the vibratory transportation and technologic machines, many unsolved problems remain in the direction of theory, design and fabrication. For example, the reasons of generation of the working member parasitic (non-working) vibrations [9] and their influence on the technologic process are studied insufficiently; such reasons may be incorrect transfer of the exciting force, features of the elastic system (springs) or other constructional errors ( Fig. 1) [10].
In Fig. 1 is presented a two-mass system "vibro-drive -working member" and are shown spatial deviations of the working member caused by various possible errors: I -nominal (designed) position of the working member; II -position of the working member ( ) considering the errors including the eccentricities , , of displacement of the center of gravity from position to point ; possible deviation of the exciting force; -possible deflection of the elastic element; , , -turnings of the coordinate axes caused by the assembly errors of the vibratory machine and transition from position to position . 1 -basic elastic system of the vibro-machine; 2 -suspensions of the vibro-machine; -exciting force; the vibro-exciter ( -) can also have the similar changes that are not shown in the drawing. The axes and of the active and reactive parts ( and ) in this figure coincide in the initial position. The main interaction of the mentioned masses is realized along ( , ) axis that is reflected in the Eq. (7). The influence of vibrations of mass on mass and respectively on mass in other directions is small in contrast to the influence of mass on mass , because in addition to elastic forces the interaction of inertial forces takes place here.
The aim of the work is to draw up the generalized mathematical model allowing to study the influence of the parameters and characteristics of masses of the system on the vibratory technologic process [11][12][13][14][15][16]. In this regard, it is expedient to consider spatial, relative-translatory and absolute movement of the three-mass dynamic system; it is analogous to the three-mass vibrational technological system, including the following parts: vibro-exciter, working member, load to be processed (to be transported).
A dynamical model of the mentioned machine ( Fig. 1(b)) is shown in Fig. 1(a), where to each mass the corresponding coordinate systems , , are connected; -immobile (inertial) coordinate system; 1 -basic elastic system connecting a vibro-exciter to the working member, 2 -suspension of the vibratory machine, 3 -conventional elastic system connecting a friable load to the working member surface; -exciting vibratory force. Though these masses are integrated into the common system, each of them is characterized by the proper physical and mechanical properties significantly different from each other; they should be taken into account at drawing up a generalized mathematical model of their movement.
For the inclusion of the technologic load (mass ) in the common spatial system (Fig. 4) and imparting it a generalized character, we present it as a rigid body connected to the WM ( ) by the conventional elastic system 3 ( Fig. 2), describing elastic and damping properties of the friable material.
At a fixed moment of time, elastic system 3 (as well as 1 and 2 in Fig. 2 and 3) is decomposed into three components, describing elastic properties of the material in space. A distinctive feature of the elastic system 3 is a non-retaining character of its connection with the WM in dynamics. With the help of the elastic-damping elements are described characteristics between the layers of the technologic friable materials, interaction between the layers and between the lower layer and WM. In contrast to the existing models [1,2,14,16] all degrees of freedom are considered in it, i.e. it can be enclosed in the model of the general spatial system (Fig. 2, 3) and depending on the concrete problem, it can be reduced to the simpler form (plane, linear, null-dimensional). The deformation of the friable TL layer is modeled at movement by the elastic elements with coefficients of elasticity , , , , , (Fig. 3). The dissipation of energy at the deformation of the layer is considered by dampers with coefficients of resistance , , , , , (not shown in Figure). Therefore, direct contact of the TL with WM is replaced by elastic-frictional connections.

Drawing up the mathematical model
For the deduction of the equations of movement, consider a spatial dynamical model of the system (Fig. 4).
The masses are considered in two positions ( Fig. 4(b)): I -ideal-immobile, determined according to the design drawings, II -dynamical displacement under the action of the exciting force determined by the direction cosines in accordance with Tables 1 and 2, where 1, 2, 3 (mass numbers), I, II (positions of masses). Direction cosines in Table 2 are presented in the form of the angles of Euler-Krylov [17], where, because of small values of the rotary displacements, products of the angles no greater than second degree are taken into account at the expansion. Table 2 shows direction cosines for mass (working member) only. For other masses, forms of expressions of the angles will be similar.
After dynamical displacement the positions of the coordinate systems will be: , , (Fig. 4(b)).  (centers of masses),  ,  , -, ,radius-vectors of the points relative to the origin of the coordinate axes of the corresponding masses, , , -, , -radius-vectors of the points relative to the origin of coordinate axes of own masses.
With the help of Fig. 2, 4 and methods of deduction of equations of motion of the material points and rigid bodies [5,17], we obtain equations of spatial movement of masses and . The following ideas were taken into account in the establishment of the differential equations of motion: -Spatial motion equations were established for the interconnected motion masses of the working body ( ) and the technological load ( ), as relative-translatory in relation to each other.
-Because is connected to mass by the potential and damping forces ( = , . , = 1, 2; = , , … , ) whose influence on the motion of the TL is insignificant, its spatial movement is described by linear differential Eq. (2); at the same time, as an electromagnetic vibro-drive realizes vibratory motion through angle with respect to mass (to the axis of WM) by means of main elastic system 1 ( Fig. 1 and 4).
-The influence of the mass of the load ( ) on the working body ( ) is indicated by coefficient , the value of which could vary between 0-1 depending on a regime of the load movement (sporadic or with working body).
-In Eqs. (1) and (3), the multiplication type 2nd order nonlinear members will be considered mostly in the high amplitude (resonance) vibrations.
Working member motion equations will be as follows: where = + , -the angle of inclination of the working member, -the angle of vibrations, , , -components of the exciting force; , , -elastic and damping forces from elastic systems 1 and 3 and weight of the TL; , , … -sum of the moments of inertia of masses relative corresponding axes, e.g.: In the equation systems Eqs. (1) and (3), not all rotational movement equations are shown (in , , , directions). Spatial vibratory movements of mass are described by the equations: where , , are principal moments of inertia; , , -components of the exciting force; , , -elastic and damping forces from elastic systems 1 (they are functions   The numerical values of stated (empirical) coefficients [3,5] correspond to the ordinary coal with 5 % humidity, or ore with maximum fragment size no greater than 100 mm and with 2-3 % humidity.
For the description of the friction force between the WM and TL of the friable type various approaches [3,[18][19][20] are used, the essence of which is that reaction of the TL on WM is proportional to the velocity and deformation of TL (Fig. 3): that will be written in the expanded form as follows: where takes the values , , ; consequently, friction forces of the TL on the surfaces of WM (e.g. of the tray-shape form) and their moments will have the form: where -coefficient of friction of the TL on the surface of WM; -distance between the friction surface and center of mass of the TL along the coordinate .
Since the rotary movement of the material is relatively small than the linear one, and only the material's longitudinal displacement is taken into consideration, only linear coordinates are provided for in the material friction forces Eqs. (5), (6).
Decomposing expression ( ) along the surfaces of the WM, we obtain: where , , -coefficients of friction along , , ; , -normal reactions on the surfaces of the WM with tray-shape form; " " -non-linear function depending on sign of velocities , , : = 1, at ( , ) < 0; = -1, at ( , ) > 0. Moments of the friction forces will be determined by putting projections Eq. (6) into expressions of Eq. (5). As it was mentioned, the TL (friable material) is unilaterally connected to WM (Fig. 3). It is periodically compressed and released depending on the action of the working member, or the pressure on the WM is increased and decreased (Fig. 2). The coefficients and , change respectively in the expression = + as well as in Eq. (3) in the process of modeling depending on ; When ≤ 0 (TL is displaced together with WM) values of and are significantly greater than when > 0 (TL loses contact partially or fully with WM) and these changes as well as the coefficients are different depending on the TL properties (dry, humid etc.) [3].
For illustration, in Fig. 5 is shown the trajectory of the friable material layer ( ) relative to the working member vibration ( ) [21].
Analysis of the oscillograms shows that the vibratory movement trajectory of the friable material layer, at various amplitudes of the working member vibration, is different: a) the material moves on the working surface with slipping without losing the contact with it; b) with increase of the amplitude the material loses the contact partially with the surface; c) with further increase of the amplitude a single touch of the material with the surface takes place during one period of the vibration followed by its jump (vibration velocity at losing the contact = ); d) with increase of the vibration amplitude (with increase of the velocity) the length of the material jump increases but since the vibration velocity is small at the moment of the material falling it cannot lose the contact until the velocity reaches the maximum again (the oscillogram of the vibration velocity is not given in the Figure and it is meant that it is shifted from the displacement by the phase 90°).
The modeling was carried out according to the parameters of the resonant vibro-feeder with 500 W electro-magnetic vibro-exciter [5,[23][24][25][26]. At that, equations of movement of mass are replaced by the equation of variation of the electromagnetic flow Ф [27]: depending on which the exciting force is determined: where -number of windings of the coil, -initial clearance between the magnet pole and anchor; -magnet penetration of the air; -area of the electromagnet core section; -total resistance; -displacement of the reduced mass of the active and reactive parts of the vibratory feeder.
The following assumptions were considered during the modeling: -Movement of the material is presented at modeling by system Eq. (3) and reaction of the TL on WM is not provided for in Eq. (1) of the working member or = 0, which we often find in the researches [1][2][3][4]7] and it does not change the qualitative characteristics of the results (in general, varies between 0-1 and its value depends on the type of friable material and vibration regime: that is the next stage of our research and the subject of our next publication).
-Only linear coordinates are provided for in the material friction forces Eqs. (5), (6), which is explained above. Fig. 6 shows some oscillograms of movement of the vibro-feeder with the load: vibrations of the working member -, displacement of the friable material -and variation of the electromagnetic flow Ф.
The mentioned oscillograms are obtained from the solutions of Eqs. (1), (3), (7) and (8) equations, when special non-working vibrations , , , , are far from the working resonance (50 Hz). The process of the vibro-feeder's entering the resonance, when amplitude of the working member and respectively angle of rise of the material displacement increase, is depicted in the picture; It should be noted in this connection that the equation of the material longitudinal displacement (first equation of (3)) does not include elastic deformation of the material (Fig. 3) because displacement of the material is not restricted in this direction (the working member groove is open in the longitudinal direction) in contrast to directions and , when the working member bottom and walls restrict the material displacement and therefore, their equations include and ; Consequently, the material vibratory displacement has the form [3], shown in Fig. 6 and with increase of the velocity, angle of its rise increases. The mathematical model was adjusted to reach the identicality with the physical experiment [3,5,15] to achieve the reliability of the modeling results: a vibratory displacement velocity of the friable material (ground coal) was measured on the real vibro-feeder. Then a simple numerical experiment was carried out by variation of the material model coefficients (Fig. 3)   The graphs on Fig. 9 are obtained from the solutions of Eqs. (1), (3), (7) and (8), when all of the non-working vibrations , , , , are far from the working resonance regime ( = 4 mm, = 50 Hz) and their influence on the process is insignificant; the graphs show the dependence of the friable material's velocity on the angle of vibration ( ) at fixed values of the angle of inclination of the working member ( ). As it is seen from the graph, the best result is reached at angles: 10° and 10°. The working resonance frequency of the WM during the modeling process is constant: = 50 Hz = const. The results of the modeling have shown that some combinations of the spatial (non-working) vibrations and basic working vibrations may have a positive influence on the transportation and technologic process of the material (increase both, transportation velocity of the friable material and the intensity of movement). Fig. 10 shows a design of the vibro-feeder with a new vibro-exciter. It is developed on the basis of the results of the modeling and it allows generation of various forms (I, II, III) of the working vibrations with the help of the variation of the angle between electromagnetic and elastic forces and thereby control the vibratory technologic process of the friable material.

Conclusions
1) A new mathematical model of spatial movement of the three-mass system -analog of the vibratory (Fig. 4) machine with technological load, is developed. A system of differential equations describes the interconnected movement of component parts of the loaded machine and allows many-sided research into vibratory transportation and the technologic process of the friable materials.
2) The presented graphs of the results show that mathematical modeling may reveal new nuances promoting improvement of both, the technologic process and the machine design.
3) It was, for example, ascertained that some spatial non-working (parasitic) vibrations in combination with basic (working) vibrations can significantly increase velocity of the material's vibratory displacement (see Figs. 7, 8).
4) A new design of the electromagnetic vibratory feeder was developed on the base of the modeling results that allows obtaining various configurations of the working vibrations and increasing velocity of the material transportation by the change of the elastic force direction. The obtained design (Fig. 10) confirms reliability of the developed mathematical model because the influence of the vertical partial (non-working) vibration is qualitatively identical for both -the modeling (Fig. 7 -) and the design (III- * ).