Modelling and processing for condition monitoring and technologies analysis
Автор: Anakhin Vladimir D.
Журнал: Вестник Бурятского государственного университета. Математика, информатика @vestnik-bsu-maths
Рубрика: Математическое моделирование и обработка данных
Статья в выпуске: 4, 2017 года.
Бесплатный доступ
In this article the most general model that can be used is one in which the basic parameters we need in equations for describing advanced vibratory processing machinery are estimated, based upon obtained exact solutions for a relatively few second-order differential equations; the more important of these exact solutions are listed. They are exact in the sense that the solution is given in an expression that can be evaluated numerically to any desired degree of accuracy. The model of VSM that experiences simple harmonic motion and the governing equations sufficient to describe the behavior of the physical system adequately are used as the model for parameter estimation and to study equipment design problems.
Condition monitoring, modeling and vibratory processing, machinery analysis, optimization in cm
Короткий адрес: https://sciup.org/14835241
IDR: 14835241 | DOI: 10.18101/2304-5728-2017-4-64-72
Текст научной статьи Modelling and processing for condition monitoring and technologies analysis
In the field of the mechanical design and application of the vibration forces arising from any type of dynamic excitations for technologies mathematical concepts involved in theoretical analysis of the method are described in the papers [1, 2]. Process motivated by the acceleration of a deck undergoing harmonic motion occurs when the deck surface has both longitudinal and lateral tilts and follows a displacement pattern in directions: (1) not having a lateral tilt, and (2) defined by a lateral tilt. Interesting surveys of the entire realm of the occurrence of new phenomena are predicted, are explainable, and commercially available with particular application associated with the design and operation of sufficiently new equipment. Every phase of the above mentioned process is affected by the goals which are established, particularly with regard to that phase upon (1) that has been studied and the greatest progress has been made in treating vibration technologies. This item describes largely in mathematical terms certain of the more important features of vibration theory likely to be encountered in practice, and provides furthest information concerning features of the methods of analysis which find ready application.
The systems treated are systems with a finite number of degree-of-freedom which can be defined by a finite number of simultaneous ordinary differential equations. It should be pointed out that vibration processes are based on the velocity transfer from a vibrating deck to the flow of material on it due to the exis- tence of a vibrating force. In practice, this requires the intimate contact of the deck and particles for some period of time, during which equilibrium is approached as velocity transfer proceeds. Process synthesis represents the inventive aspect of process design. Because of to maximize efficiency, the contacting time in a particular stage may not be long enough. It is usually desirable to use vibration defined by the acceleration amplitudes of the deck expressed as a dimensionless multiple of the gravitational acceleration in a range from 6 to 10g. At such values of vibration acceleration particles will be intensively tossed upwards after a brief contact with the deck.
The plane of the X and Z axes is a plane of symmetry. The vibration force is acting along the X axis at the angle of vibration to its direction lengthwise and to the plane of the deck. The angle of lateral tilt of the deck surface (angle of ε ) can not be varied. The key difference over the design of machinery is the following: there is no inclination in direction of the acting force at the angle of ε .
Suppose that X, Y, Z fixed in the deck is a convenient set of axes with its origin at the point of material loading. The Y axis is perpendicular to the plane of the deck; the X and Z axes are oriented along and transverse. The direction of the motion of the deck and its change of position are characterized by the angles of vibration β , longitudinal and lateral tilts a and ε , respectively. The values of these angles can be varied as follows
- π /2 < a < π /2,
0 < β < π /2 - a , 0 < ε < π /2.
(1.1)
The initial conditions are when t0= t0* : x = 0, y = 0, z = 0, dx dx * * dy dy * *dz
= ≡(vx)*0, = ≡(vy)*0, dt dt 0 dt dt 0 dt
dz dt
*
≡ ( vz )*0.
The differential equations of particle motion in the X, Y, Z directions are d2x dt2
=- g sin a + A ω 2 cos β sin ω t ,
d 2 y dt 2
= ( - g cos a + A ω 2sin β sin ω t )cos ε ,
(1.2)
d 2 z dt 2
=(gcosa-Aω2 sinβsinωt) sinε, where g is the gravitational acceleration; A is the amplitude of the displacement of the deck; ω is the angular frequency of the simple harmonic motion.
By integrating once and twice, respectively, the equations are written as follows
*
dx ( t ) dx
=- g ( t - t )sin a - A ω cos β (cos ω t - cos ω t ) + , dt 00 dt 0
( t - t *)2
x ( t ) =- g 0 sin a - A cos β (sin ω t - sin ω t 0*) +
(1.3)
+Aω(t-t0*)cosβcosωt0*+(vx)*0(t-t0*), dy(t)
=-g(t-t)cosacosε-Aωsinβcosε(cosωt- dt 0
*
dy cos ω t ) + , 0 dt 0
( t - t *)2 *
y ( t ) =- g 0 cos a cos ε - A sin β cos ε (sin ω t - sin ω t *) +
2 0 (1.4)
+Aω(t-t0*)sinβcosεcosωt0*+(vy)*0(t-t0*), dz(t)
=- g ( t - t )cos a sin ε + A ω sin β sin ε (cos ω t - dt 0
*
dz cos ω t ) + , 0 dt 0
( t - t *)2
z ( t ) = g 0 cos a sin ε + A sin β sin ε (sin ω t - sin ω t 0*) -
- A ω ( t - t 0*)sin β sin ε cos ω t 0 * + ( vz )*0( t - t 0*).
The phase angles at the time of flight initiation of the particle over the vibrating deck and of its downfall on the deck, respectively, are **** δ 0 = ω t 0, δ n = ω tn .
For the period of particle flight over the deck defined by T = pT 0 , where T 0 is the time interval to complete one cycle in the deck vibration, the following expression can be written
δ n - δ 0 * = 2 π p .
(1.6)
The velocity of particle fall down on the surface of the deck ( vy ) n and the velocity of particle reflection and flight initiation ( vy )*0 are known, in accordance with Newton's theory of impact, as
( v y )* 0 R =
.
( vy ) n
Using expression (1.6), the following is obtained
2 π p
ω tn = 2 π p + ω t 0 , t = tn =+ ω t 0. ω
(1.7)
(1.8)
The Eq. (1.4) give the following
( vy ) n = - g ( tn - t 0*)cos a cos ε - A ω sin β cos ε (cos ω tn - cos ω t 0*) + ( vy )*0. (1.9)
Since cosωtn =cos(2πp+ωt0) =cosωt0, this gives
( v y ) n = - g ( t n - t 0* )cO s « cos S + ( v y ) 0 .
Using relation (1.8), the following result is obtained
, X (2np к V
( v y ) n =- g I I cos « cos S + ( v y ) 0 .
к ® J
In accordance with equality (1.7) the preceding equation can be written
. .*f. 1 1 f 2 n pg 1
(v,,)n 1 + — = ------ cos« cosS, y0
к R J к ю J so that the equation can be expressed as
* ( 2 n pg Y R к
(1.10)
( v ,,)n = ------- ------ cos « cos S .
y 0
к ю Jk 1 + R J
Eq.(1.8) makes tn -10 = —p, substituting this relation, Eq (1.4) becomes ю f 2n p 1 cos a cos s , . „ , . , . A .
-g I —— I--A sin в cos s(sin юt - sin юt0) + к ю J 2
. , 2np • n . 2np, „
+ А ю cos s ——sin в cos ю t 0 + —— ( v )0 = 0. ю ю y
Since cos ю t 0 = cos 5* 0 and sin ю tn = sin ю t 0 , dividing 2 n p and
A sin в cos s , the following equation is obtained
*
-g cos a (vy )o
(1.11)
—2------np + cos 50 +-----y=
А ю sin в А ю sin в cos s
, _ Аю2 sin в • .
where ю 0 =--------- is the coefficient of dynamic load.
g cos a
The phase angle of flight initiation of the particle over the vibrating deck can be calculated by the formula np L 2R к cos 50 =——I 1--I юо к 1 + R J or
o* n p 1 - R
(1.11a)
cos 5 =--- 0 ю 0 1 + R
The particle displacement in the X direction at the time interval to complete one cycle of the particle flight, using the appropriate Eq. (1.3), is
( t„ - t n
S X = x ( t ) = - g n ^ 0 7 sin a - A cos в (sin to t n - sin to t 0 ) + + A to ( t n - t О ) cos в cos to t О + ( v x ) 0 ( t - t * ).
' - 5* .
n 0 , this can be written го
5 * 5
Since t n = —, t 0 = — and t n - t 0 = — to го
S
X _
0.5 g
sin a ( 5 n - 5 0 * )2 + A cos в (sin 5 n - sin 5 0* ) +
+ A cos в cos 5 0* ( 5 n - 5 0* ) + ( v x M f n— 5 0 ).
to
The phase angles at time of flight initiation of the particle over the vibrating deck and of its downfall on the deck are related as follows
5 n = 2 п P + 5 0* .
Using this relation,
( 2 n p 1 2
S X = 0.5| —— I g sin a + A cos в P sin(2 n p + 5 0 ) - sin 5 0 1 +
V to ) L J
, c* , 2 np / X*
+ 2 n pA cos в cos 5 0 +---- ( vx )0
to
Since sin(2np + 50*) - sin 50* = 0,
_* g cos a
( v x )* 0
A to sin в cos s’
cos 5 =--- x----n p -
0 Ato2 sin в from the previous equation
S x =
—
g ( 2 n p ) . , . o
—|—— I sin a + 2 n pA cos в 2 1 to )
g cos a ------np A to sin в
( v y )* 0
A to 2 sin в cos s
+
2л p . .* g ( 2 np 12 ( cos a . 1 2 np (, .* (К )0
+— -( Vx )n =— —--- sin a +—- ( vx )0---- x 0 x 0
to 2 V to ) ^ tan в ) to ^ cos s tan в )
cos a . cos( a + в ) .
Letting--sin a = — ------, the preceding equation can be written tan в sin в
, _ g ( 2 n p 1 2 cos( a + в ) + 2 n p X 2 V to ) sin в to
( vx )*0
( V y )0 "
cos s tan в )
(1.12)
Eq.(1.10) gives
( V y ) 0 _ ^ ^ pg to
R
-----cos a cos s . R + 1
Similarly, vx ** becomes
VX* = ( v x ) 0 = (1 — g sin a , (1.13)
to where X is the coefficient of instantaneous friction. Substituting these terms (vx)*0 and (vy)*0 into (1.12), the following equation is obtained
„ g ( 2np Y ( cos(a + В) . 1 - X R cos a А
S = — sin a•
X 2 ( to Д 2sin в X R +1 tan в J.
The mathematical equation of particle velocity for motion in the direction of the X-axis
V = S x to = 2 n pg
X 2 n p to
n (1 R А ■ (1 1—X cosacot в---— sin a---
(2 R +1J (2 X
Simplifying this, the equation can be reduced to n pg (1 — R , „ 2 — X . А
V X =-----1----- cos a cot в --sin a l . (1.14)
to V1 + R X )
Since Eqs. (1.3) and (1.5) are of the same form, the solutions are the same except for the directions of particle motion. The solution of Eq.(1.5) is of the form in which the particle displacement in the Z direction is considered
S Z = z ( t n ) =
= |^0.5 g ( t n — t 0 )2 cos a + A sin в (sin to t n — sin to t 0 ) — A to ( t n — t * )sin в cos to t 0 J cos s +
+ ( v z ) 0 ( t n — t 0 ).
Using the relations of Eqs. (1.6) and (1.8), this equation can be expressed as follows
S z = 0.5 g
2 n p Y т л n<-* ----l cos a — 2 n pA sin в cos o 0 to ;
. 2 np . .* sin s +----- ( v z ) 0 .
to
(1.15)
It is assumed that
( v x' ) n
(v z j n = 1 — X . ( vz ) n
(1.16)
( vx ) n
After a determined interval of time At = 2np, the increment of particle veto locity toward the lateral direction is
A( Vz ) n = (Vz' ) n — (Vz )0 = (Vz ) n (1 — X) — (Vz )0, where (vz )*0 is the projection of particle velocity on the Z axis at the time of flight initiation. Using Eq. (1.5), the following is obtained
A ( v z ) n = [ g ( t - t 0 ) cos a + A to sin в (cos to t - cos to t 0 ) J sin s (1 - X ) +
* * 2n p , (1 .1 7)
+ (vz )0(1 - X) - (vz )0 =----g cos asin s(1 - X) - X( vz )0.
to
Projection of particle velocity on the Z axes is lx 2n p
)—— g cos a sin s .
vz = ( vz ) 0
to
Substitution this relation into Eq. (1.15) gives
Sz =
g ( 2np I . . . o .
—I—— I cos a - 2 n pA sin в cos 5 0 sin s +
21 to J
(2 n p )2 to T
1 - X x ^x g cosasin s=
g cos a
(1.18)
2 n pA sin в cos 5 0* sin s .
Since cos 5 0 is determined by (1.11a), Eq.(1.18) can be written
S Z = 2 g cos a
= 2 g cos a sin s
n 2 p p to 2
n 2 p p to 2
2 - X „ n2 p2 1 - R | .
---2 g cos a — --- sin s =
X to2 1 + R J
2 - X _ 1 - R
X 1 + R
(1.19)
which can be reduced to
Vz =
SZto
n pg [ 2 - X 1 - R i
—— I---I cos a sin s .
2np to ( X 1 + R J
(1.20)
Applicable mathematical expressions
The equations (1.14) and (1.20) for the particle velocity in its movement in the longitudinal and lateral directions are the basic equations in the theory of vibration separation. For the purpose of studying the effects of vibration these equations can be transformed to other forms required for particular applications.
At a specific instant of the phase angle corresponding to the particle flight initiation motivated by the maximum value of the acceleration of vibration, when cos 5 0* = 1 , from the expression (1.11a)
1 < 1 + RI p = --- to0.
n (1 - R J
In this case, if it is assumed that
X 1 - R q =--- ,
2 - X 1 + R is the parameter of separation, Eqs. (1.14) and (1.20) can be converted to the form
VX = Ato
cos в -I
sin в tan a |
I, q J
(1.21)
I 1
VZ = ^sin e sin s —1 I .
I q J
(1.22)
Dividing Eq. (1.22) by Eq.(1.21)
sin в sin s | 1 -1 | dz _ VZ ^q J (1 -q)sins
dx Vx n sin e tan a q cot в - tan a
X cos в
.
(1.23)
q
The Eq. (1.23) can be used to construct the trajectory of particle movement. This is done as illustrated in Fig.1 of [7]: resulting trajectories are drawn by broken-dotted lines through a point taken near the origin of the separation load. The directions of movement of the components are such as to take the trajectories away from the origin. The angle of particle orientation in its movement is indicated by the formula
(1 -q)sins tan^I = ———-----. (1.24)
q cot в - tan a
The XZ plane of the deck in is a plane of symmetry. Each parallel deck with its longitudinal and lateral axis in the X and Z directions, respectively, has the length l and the width b of the deck, respectively. So x = b /tan ^ , z = l tan ^ .
Separation ability is defined analytically as follows
= dx X = dq
z = b
Dz =- dz dq
, d I 1 ) = b^HI, dq ^ tan^ J
1 d
--( tan w ) .
x = 1/2
2 dq( )
Solving Eq.(1.23)
D X I
= bd । q cot в - tan a dq ( (1 - q )sin s
b cot в - tan a sin s (1 - q )2
1 d I (1 - q )sin s ) 1 . cot в - tan a D Z = ---———-----I = —sin s -----------^.
I 2 dq ( q cot в - tan a J 2 ( q cot в - tan a )
(1.25)
(1.26)
Conclusions
-
1. In the operation of modeling and processing for CM, there are several factors to be considered that are important regardless of the particular methodology. The most important consideration involves using the basic equations for evaluating the design parameters technologies.
-
2. Mathematical model method for parameter estimation is used as the theoretical basis for the description of processing for CM with the second-order differential equations whose solutions are obvious from adequate analytical models.
-
3. The method has been applied extensively in studies of the governing equations and some of the general mathematical expressions thereby obtained are given here, both as for future evaluations and as a collection of useful relations which find ready application in practice.
-
4. The solution procedure involves representation of the function behavior in the form of response curves versus parameters. The response curves for a particular system acted on by a harmonic exciting force may be derived from the obtained equations for the condition monitoring.
Список литературы Modelling and processing for condition monitoring and technologies analysis
- Anakhin V. D. Procedures for the design of vibratory equipment//Integrating Dynamics, Condition Monitoring and Control for the 21st Century/In A. Starr et al. (eds). Rotterdam: Balkema, 1999. P. 451-457.
- Anakhin V. D., Anakhin Т. V. Process consideration and concepts in data analysis. Proceeding of the Condition Monitoring. Oxford, Coxmoor publishing Company, 2005. 490 p. P. 85-91.