Wave displacement of rigid bodies and particles

Kazimieras Ragulskis1 , Liutauras Ragulskis2

1Kaunas University of Technology, Kaunas, Lithuania

2Vytautas Magnus University, Kaunas, Lithuania

Journal of Vibroengineering, Vol. 15, Issue 2, 2013, p. 794-804.
Received 2 April 2013; accepted 3 June 2013; published 30 June 2013

Copyright © 2013 Vibroengineering This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Creative Commons License
Abstract.

The displacement of various particles and mechanical systems by waves finds wide use in technical devices, technological processes, takes place in non living and living nature. Here the system displaced by waves with four degrees of freedom is analyzed when one member of the system, while contacting with the working profile of the input member performing wave motion, provides motion to the output system.The obtained differential equations of motion of the system have been analyzed analytically and numerically. For the analytical investigation modification of the asymptotic method is used, which is based on the division of motion into the slow and quick motions. This method is justified when the frequencies of variation of slow motions are much smaller than the frequencies of quick motions.When the working surface of the input member moves according to the harmonic Rayleigh waves more detailed full investigations have been performed. The characteristics of transition to steady state and of stationary regimes of motions have been obtained, such as the conditions of existence and stability of the stationary regimes, curves of bifurcation, existing complicated motions.

Keywords: traveling waves, transient motions, stationary state, stability, bifurcation.

1. Introduction

The conveyance of particles and bodies by propagating waves is an important scientific and engineering problem with numerous applications. Manipulation of bioparticles and gene expression profiling using travelling wave dielectrophoresis [1-3], segregation of particles in suspensions subject to alternative current electric fields [4], transport of sand particles and oil spills in coastal waters [5, 6], powder transport by piezoelectrically excited ultrasonic waves [7, 8] are just few examples of problems involving interaction between propagating waves and transported objects.

There is a large variety of waves of various types according to their shapes and laws of propagation [9-13]. This paper is focused on the displacement of systems with lumped parameters by one-dimensional periodic waves. This is important in a number of mechanical engineering applications. For the case of Rayleigh harmonic waves the investigations are performed in more detail.

2. The model of the system

The model of the system is presented in Fig. 1.

2.1. Kinematics

The coordinates of the moving points in the frame XOY read:

(1)
A 1 u , ξ ,       A 2 u , v ,       B 1 u + u 1 , v + v 1 ,       B 2 u - R s i n φ , v - R c o s φ ,

where:

(2)
u = x + η ,       v = ξ ,       R = A 1 B 2 = A 2 B 2 .

The working profile of the input member 1 is given by:

(3)
η x , t , ξ x , t ,

and the components of velocities of their contact point A1 are:

(4)
η ' t , ξ ' t .

Fig. 1. The model of the system: 1 – the working profile of the input member, 2 – contact member of the output system, 3 – member 1 may move according to two orthogonal directions with respect to member 2, 4 – a pendulum is attached to the member 2 by a hinge

 The model of the system: 1 – the working profile of the input member, 2 – contact member of the output system, 3 – member 1 may move according to two orthogonal directions with respect to member 2,  4 – a pendulum is attached to the member 2 by a hinge

Fig. 2. Forces acting in the system

 Forces acting in the system

The components of velocities and accelerations of contact point A2 of member 2 read:

(5)
u ˙ = 1 + η ' x x ˙ + η ' t ,       v ˙ = ξ ' x x ˙ + ξ ' t ,
(6)
u ¨ = 1 + η ' x x ¨ + η x x x ˙ 2 + 2 η x t x ˙ + η t t ,       v ¨ = ξ ' x x ¨ + ξ x x x ˙ 2 + 2 ξ x t x ˙ + ξ t t .

Further the following relationships are met:

(7)
F 1 = 1 + η ' x u ¨ + ξ ' x v ¨ = a 1 x ¨ + a 2 x ˙ 2 + 2 a 3 x ˙ + a 4 ,
F 2 = - ξ ' x u ¨ + 1 + η ' x v ¨ = b 2 x ˙ 2 + 2 b 3 x ˙ + b 4 ,

where:

(8)
a 1 = 1 + η ' x 2 + ξ ' x 2 ,         a 2 = 1 + η ' x η x x + ξ ' x ξ x x ,
a 3 = 1 + η ' x η x t + ξ ' x ξ x t ,       a 4 = 1 + η ' x η t t + ξ ' x ξ t t ,
b 2 = 1 + η ' x ξ x x - ξ ' x η x x ,       b 3 = 1 + η ' x ξ x t - ξ ' x η x t ,
b 4 = 1 + η ' x ξ t t - ξ ' x η t t .

The angle α of the tangent at the point of the profile A1 with the axis OX:

(9)
t a n α = ξ ' x 1 + η ' x .

The velocity of slippage of point A2 with respect to point A1 in the tangential direction takes the form:

s ˙ 21 = ( v ˙ - ξ ' t ) s i n α + ( u ˙ - η ' t ) c o s α ,

or:

(10)
s ˙ 21 = a 1 x ˙ .

2.2. Differential equations of motion

The equations of dynamic equilibrium of the member 2 (Fig. 2) are:

(11)
m + m 2 u ¨ + m 2 R P u - F u + f 1 s ˙ 21 c o s α = - N s i n α + f 0 s g n s ˙ 21 c o s α ,
m + m 2 v ¨ + g + m 2 R P v - F v + f 1 s ˙ 21 s i n α = N c o s α - f 0 s g n s ˙ 21 s i n α ,

where N is the normal force acting to the member 2 at the points A1 and A2, f0 and f1 are the coefficients of dry and viscous friction:

(12)
P u = - φ ¨ c o s φ + φ ˙ 2 s i n φ , P v = φ ¨ s i n φ + φ ˙ 2 c o s φ ,
F u = H u u ˙ 1 + C u u 1 ,       F v = H v v ˙ 1 + C v v 1 ,

and g is the acceleration of gravity of the Earth.

Rearrangement of the equations (11) and taking into account the equations (2)-(9) yields:

(13)
1 + η ' x I u + ξ ' x I v + - ξ ' x I u + 1 + η ' x I v f 0 s g n s ˙ 21 + h α 1 x ˙ = 0 ,
(14)
R φ ¨ - u ¨ c o s φ + v ¨ s i n φ + g s i n φ + K φ m 2 R φ ˙ = 0 ,
(15)
u ¨ + u ¨ 1 + F u = 0 ,
(16)
v ¨ + v ¨ 1 + F v = 0 ,

where:

(17)
I u = u ¨ + μ 2 R P u - μ 1 F u m 1 , I v = v ¨ + g + μ 2 R P v - μ 1 F v m 1 ,
μ 1 = m 1 m + m 2 , μ 2 = m 2 m + m 2 , h = f 1 m + m 1 ,        ˙ = d d t ,

K φ is coefficient of viscous friction.

3. Systems excited by traveling waves

In this case the equations (3) take the form:

(18)
η = η ω t - k x ,       ξ = ξ ( ω t - k x ) ,

where η and ξ are periodic functions of their arguments.

The following notations are introduced:

(19)
ψ = ω t - k x ,       τ = ω t ,       β = k x ,       β 1 = k u 1 ,       β 2 = k v 1 ,     ' = d d τ .

In this case equations (5)-(10) by taking into account the equations (18), (19) take the following form:

(20)
k u ' = 1 - c 4 ψ ' ,       k v ' = c 5 ψ ' ,       k u = - c 4 ψ + c 6 ψ ' 2 ,       k v = c 5 ψ + c 7 ψ ' 2 ,
k F 1 = - c 1 ψ + c 2 ψ ' 2 ,         k F 2 = c 3 ψ ' 2 ,

where:

(21)
c 1 = 1 - k η ' 2 + k 2 ξ ' 2 ,       c 2 = 1 - k η ' k η - k 2 ξ ' ξ ,       c 3 = 1 - k η ' k ξ + k 2 ξ ' η ,
c 4 = 1 - k η ' ,         c 5 = k ξ ' ,       c 6 = k η ,         c 7 = k ξ ,         r = k R .

Differential equations of motion (13)-(17) by taking into account equations (18)-(21) read:

(22)
- c 1 ψ + c 2 ψ ' 2 + c 4 μ 2 r p u - μ 1 f u - c 5 k g ω 2 + μ 2 r p v - μ 1 f v +   - c 3 ψ ' 2 + c 5 μ 2 r p u - μ 1 f u + c 4 k g ω 2 + μ 2 r p v - μ 1 f v f 0 s g n s ˙ 21 + c 1 h ω 1 - ψ ' = 0 ,
(23)
r φ + ( c 4 ψ - c 6 ψ ' 2 ) c o s φ + c 5 ψ + c 7 ψ ' 2 s i n φ + h φ ω φ ' + k g ω 2 s i n φ = 0 ,
(24)
- c 4 ψ + c 6 ψ ' 2 + β 1 + f u = 0 ,
(25)
c 5 ψ + c 7 ψ ' 2 + β 2 + f v = 0 ,

where:

(26)
p u = - φ c o s φ + φ ' 2 s i n φ ,       p v = φ s i n φ + φ ' 2 c o s φ ,
f u = h u ω β ' 1 + n u ω 2 β 1 ,       f v = h v ω β ' 2 + n v ω 2 β 2 ,
h = f 1 m + m 2 ,       h u = H u m 2 ,       h v = H v m 2 ,       n u 2 = C u m 2 ,       n v 2 = C v m 2 ,       h φ = k K φ m 2 R .

3.1. Case: harmonic traveling waves:

Following equalities are obtained according to (21) and taking into account (27):

(27)
η = A c o s ω t - k x ,       ξ = B s i n ( ω t - k x ) .
(28)
c 1 = ( 1 + a s i n ψ ) 2 + b 2 c o s 2 ψ ,       c 2 = - ( 1 + a s i n ψ ) a c o s ψ + b 2 s i n ψ c o s ψ ,
c 3 = - ( 1 + a s i n ψ ) b s i n ψ - a b s i n ψ c o s ψ ,
c 4 = 1 + a s i n ψ ,       c 5 = b c o s ψ ,       c 6 = - a c o s ψ ,       c 7 = - b s i n ψ ,       a = k A ,       b = k B .

3.2. Analytical investigation for the case of slow motions when the average velocity of the output system is much smaller than the velocity of traveling waves

In this case the steady state motions and the transient ones which are near to them are sought in the following form:

(29)
ψ = ψ 0 + ψ ~ ,       φ = φ 0 + φ ~ ,       ψ 0 = ψ ' - τ + ψ - ,       φ 0 = φ ' - τ + φ - ,

where ψ'- and φ'-=mnψ'- are average velocities in the steady state regimes:

(30)
m , n = 0 , ± 1 , ± 2 , ,       m n 1 ,

ψ - and φ- are slowly varying quantities, ψ~ and φ~ are quickly varying quantities.

For the determination of quick motions linear differential equations with constant coefficients are obtained. The latter are obtained from the main equations (22)-(26) by assuming in them that the dash over separate members denotes averaging with respect to time and the wave over the members means that their average values with respect to the period of time are equal to zero. Nonlinear autonomous differential equations are obtained for the determination of slow motions by assuming in the initial equations (22)-(26) ψ and φ according to the equations (30) and performing linearization over one time period.

Differential equations describing quick motions are obtained from equations (22)-(25) when f0=0 by taking into account the equations (30):

(31)
c - 1 ψ ~ + h ω ψ ' ~ = ψ ' - 2 c 2 ψ 0 + μ 1 r φ ' - 2 c 4 ψ 0 s i n φ 0 ~ - c - 4 μ 1 f u
                  - c 5 ψ 0 k g ω 2 + μ 2 r φ ' - 2 c o s φ 0 ~ + c 1 ψ 0 ~ h ω = 0 ,
(32)
r φ ~ + k φ ω φ ' ~ - ψ ' - 2 c 6 ψ 0 c o s φ 0 + φ ' - 2 c 7 ψ 0 s i n φ 0 ~ + k g ω 2 s i n φ 0 = 0 ,
(33)
- c - 4 ψ ~ + ψ ' - 2 c 6 ( ψ 0 ) + β 1 + f u = 0 ,
(34)
ψ ' - 2 c 7 ( ψ 0 ) + β 2 + f v = 0 ,

where the upper dash means averaging with respect to the period, fu and fv are determined according to the equations (26), ci (i= 1, …,7) are determined according to the equations (28), (29).

The differential equations of slow motion respectively are only for the determination of ψ- and φ-, namely:

(35)
c - 1 ψ - + h ω ψ ' - = ψ ' - 2 c 2 ψ 0 + ψ ~ - + μ 1 r φ ' - 2 c 4 ψ 0 + ψ ~ sin φ 0 + φ ~ ¯ - μ 1 c 4 ψ 0 + ψ ~ f u -
            -   k g ω 2 c 5 ψ 0 + ψ ~ - - μ 2 r φ ' - 2 c 5 ψ 0 + ψ ~ cos φ 0 + φ ~ - + h ω c 1 ψ 0 + ψ ~ - = 0 ,
(36)
r φ - + k φ ω φ ' - - ψ ' - 2 c 6 ( ψ 0 + ψ ~ ) c o s ( φ 0 + φ ~ ) - + φ ' - 2 c 1 ( ψ 0 + ψ ~ ) s i n ( φ 0 + φ ~ ) ¯
            + k g ω 2 s i n ( φ 0 + φ ~ ) ¯ = 0 .

Further systems are analyzed starting from the simplest ones.

3.3. Case 1: system with one degree of freedom

First the problem is analyzed when the particle moves with a velocity much smaller than the velocity of the wave, that is δ-1=ε0,ε is a small quantity.

In this case the differential equations of the slow and quick motions are obtained on the basis of the equations (31) and (35) by taking into account the equations (28):

(37)
1 + 0.5 a 2 + b 2 ψ ~ + h ω ψ ' ~ = ψ ' - 2 - a c o s ψ 0 - 0.5 a 2 - b 2 s i n 2 ψ 0
              - k g ω 2 b c o s ψ 0 + h ω [ 2 a s i n ψ 0 - 0.5 ( a 2 - b 2 ) c o s 2 ψ 0 ] = 0 ,
(38)
1 + 0.5 a 2 + b 2 ψ - + h ω ψ ' - = - ψ ' - 2 a cos ψ 0 + ψ ~ - + 0.5 a 2 - b 2 s i n 2 ψ 0 + ψ ~ ¯ + k g ω 2 b cos ψ 0 + ψ ~ -
              + h ω 1 - 0.5 a 2 - b 2 + 2 a sin ψ 0 + ψ ~ ¯ - 0.5 a 2 - b 2 c o s 2 ψ 0 + ψ ~ - = 0 .

Steady state regimes of quick motion according to the equation (37) are of the following type: ψ~=R1cosψ0+S1sinψ0+R2cos2ψ0+S2sin2ψ0,

and by substituting it into the equation (38), before this performing the expansion of it into a power series with respect to ψ~ and taking into account only the linear part and after averaging it is obtained:

(39)
δ = 1 - a 2 2 d 1 S - δ 2 + k g ω 2 d 2 2 + 4 h ω 2 δ 2 + h ω 2 + 0.25 a 2 1 - d 2 2 2 δ 4 + h ω 2 4 δ 2 + h ω 2 ,

where:

(40)
d 1 = 1 + 0.5 a 2 + b 2 ,       d 2 = b a .

It is clear that:

(41)
δ = 1 - ε ,

where ε is a small quantity, on the right side of the equation (39) it is possible to assume ε=0 and:

(42)
δ | δ = 1 = 1 - a 2 2 d 1 - 1 + k g ω 2 d 2 2 + 4 h ω 2 1 + h ω 2 + 0.25 a 2 1 - d 2 2 2 1 + h ω 2 4 + h ω 2 .

In the case when the particle moves with a velocity of the excited wave, it is:

(43)
δ = 0 .

In this case:

(44)
ψ = 0 τ + ψ - .

The differential equation of slow motion according to the equation (32) is the following one:

(45)
d 1 ψ - + h ω ψ ' - = a - k g ω 2 d 2 c o s ψ - + h ω d 1 + 2 s i n ψ - - 0.5 a 1 - d 2 2 c o s 2 ψ - = Ψ ψ - .

Stationary values of ψ- are found from the following equation:

(46)
Ψ ( ψ - ) = 0 .

Stable values of ψ- are those which satisfy the following inequality:

(47)
Ψ ψ - ψ - = a k g ω 2 d 2 s i n ψ - + h ω 2 c o s ψ - + a 1 - d 2 2 s i n 2 ψ - < 0 ,

and in case of the opposite sign of the inequality (47) there are unstable values.

The equation (47) may have four or two or none real solutions with respect to ψ-. This means that the motions with the velocity of the wave exist or do not exist. When the motions with the velocity of the wave exist, which of them are stable and which are unstable is determined according to the roots of the characteristic equation.

The characteristic equation at the stationary points according to the equation (45) is the following one:

(48)
d 1 λ 2 + d 1 h ω λ - Ψ ψ - = 0 ,

where:

(49)
Ψ ψ - = a k g ω 2 s i n ψ - + h ω 2 c o s ψ - + a 1 - d 2 2 s i n 2 ψ - .

Roots of the equation (49) are:

(50)
λ 1 ; 2 = - 1 2 h ω ± 1 4 h ω 2 + d 1 2 Ψ ψ - .

Unstable points ψ- separate stable regimes, among them including stable regimes taking place at the velocity of the wave and at the velocity smaller than the velocity of the wave. Curves going through the unstable points are bifurcation curves.

Because in the stationary points of ψ-, which are determined by the equation (46), may be four, two or zero real values, then also (λ1,2)j may be j= 1, …,4, j= 1, 2, j= 0.

Solving equation (45) with initial conditions from the infinitely small vicinity of the unstable point of ψ-, which is determined by the equations (46) and (47), bifurcation lines are obtained.

Because stationary regimes may exist:

(51)
ψ ' - = δ  and  ψ ' - = 0 ,

so for the global representation the phase plane is used:

(52)
ψ ' = ψ ' ( ψ ) .

From the unstable points, which are determined by the equations (46) and (47), by integrating the differential equation (22), in which pu=pv=fu=fv=f0=0 and c1,c2,c5 are determined by the equations (28), bifurcation curves are obtained.

When b=1.5a,  hω=0.5,  kgω=2 slowly changing a=a(ετ) graphical relationships have been obtained ψ'=ψ'(a),ψ'-=ψ'-(a),ψ'max-ψ'min=Δψ',ψ'=ψ'(ψ), (Fig. 3, 4) which reflect the dynamical qualities of the system.

Fig. 3. Transition to the steady state process when a= 0.25,hω= 0.5,kgω= 2 from the unstable point ψ-= 1.08825,Ψ'ψ-= 0.748257 (Fig. 3(a)) and from the stable point ψ-= 4.5421,Ψ'ψ-= –0.794566 (Fig. 3(b)), phase trajectories of motion and attractors (Fig. 3(c))

Transition to the steady state process when a= 0.25,hω= 0.5,kgω= 2 from the unstable point ψ-= 1.08825,Ψ'ψ-= 0.748257 (Fig. 3(a)) and from the stable point ψ-= 4.5421,Ψ'ψ-= –0.794566 (Fig. 3(b)), phase trajectories of motion and attractors (Fig. 3(c))

a)

Transition to the steady state process when a= 0.25,hω= 0.5,kgω= 2 from the unstable point ψ-= 1.08825,Ψ'ψ-= 0.748257 (Fig. 3(a)) and from the stable point ψ-= 4.5421,Ψ'ψ-= –0.794566 (Fig. 3(b)), phase trajectories of motion and attractors (Fig. 3(c))

b)

Transition to the steady state process when a= 0.25,hω= 0.5,kgω= 2 from the unstable point ψ-= 1.08825,Ψ'ψ-= 0.748257 (Fig. 3(a)) and from the stable point ψ-= 4.5421,Ψ'ψ-= –0.794566 (Fig. 3(b)), phase trajectories of motion and attractors (Fig. 3(c))

c)

Fig. 4. Steady state motions of the system

Steady state motions of the system

3.4. Case 2: system with two degrees of freedom according to the coordinates ψ and φ

In this case in the equations of motion (22) and (23) it is assumed fu=fv=f0=0.

Here the motions are investigated, for which:

(53)
ψ ' - = φ ' - = δ .

By taking into account the equations (27) and (28), the quick and slow equations (22) and (23) become the following ones:

(54)
d 1 ψ + h ω ψ ' = r 0 + r 1 c o s ψ 0 + s 1 s i n ψ 0 + s 1 φ s i n φ 0 + r 2 c o s 2 ψ 0
              + s 2 s i n 2 ψ 0 + r 2 ψ φ c o s ( ψ 0 + φ 0 ) ,
(55)
φ + h φ ω φ ' = r 0 * + s 1 φ * s i n φ 0 + r 2 φ * cos ψ 0 + φ 0 ,

where it is denoted:

(56)
r 0 = 1 + 0.5 ω 2 + b 2 h ω ,       r 1 = - a + k g ω 2 b ,       s 1 = 2 h ω a ,         s 1 φ = μ 2 r δ 2 ,
r 2 = - 0.5 a 2 - b 2 h ω ,       s 2 = - 0.5 a 2 - b 2 δ 2 ,       r 2 ψ φ = - μ 1 r δ 2 ,
r 0 * = - 0.5 a 2 - b 2 δ 2 cos ψ - - φ - ,       s 1 φ * = - h φ g ω 2 ,       r 2 φ * = 0.5 a - b δ .

From the equations (54) and (55) the quick motions are:

(57)
ψ ~ = R 1 c o s ψ 0 + S 1 s i n ψ 0 + R 1 φ c o s φ 0 + S 1 φ s i n φ 0 + R 2 c o s 2 ψ 0 + S 2 s i n 2 ψ 0
              + R 2 ψ φ cos ψ 0 + φ 0 + S 2 ψ φ sin ψ 0 + φ 0 ,
(58)
φ ~ = R 1 φ * c o s φ 0 + S 1 φ * s i n φ 0 + R 2 φ * cos ψ 0 + φ 0 + S 2 φ * sin ψ 0 + φ 0 ,

and the slow motions are:

(59)
d 1 ψ - + h ω ψ ' - = r 0 - r 1 ψ ~ s i n ψ 0 ¯ + s 1 ψ ~ c o s ψ 0 - + s 1 φ φ ~ c o s φ 0 -
              - 2 r 2 ψ ~ s i n 2 ψ 0 ¯ + 2 s 2 ψ ~ c o s 2 ψ 0 - - 2 r 2 ψ φ ( ψ ~ + φ ~ ) s i n ( 2 δ τ + ψ - + φ - ) ¯ ,
(60)
r φ - + h φ ω φ ' - = r 0 * + s 1 φ * φ ~ c o s φ 0 - - r 2 φ * ( ψ ~ + φ ~ ) s i n ( 2 δ + ψ - + φ - ) ¯ .

By taking into account the equations (57) and (58), the equations (59) and (60) become the following ones:

(61)
d 1 ( ψ - + h ω ψ ' - ) = r 0 - 0.5 r 1 S 1 + 0.5 s 1 R 1 + 0.5 s 1 φ R 1 φ *
                + r 2 - S 2 - R 2 ψ φ sin ψ - - φ - - S 2 ψ φ cos ψ - - φ -
                + s 2 R 2 + R 2 ψ φ cos ψ - - φ - - S 2 ψ φ sin ψ - - φ -
                + r 2 ψ φ R 2 sin ψ - - φ - + S 2 cos ψ - - φ - - S 2 ψ φ - S 2 φ * ,
(62)
r φ - + h φ ω φ ' - = s 0 * + 0.5 s 1 φ * R 1 φ * + r 2 φ * R 2 sin ψ - - φ - - S 2 cos ψ - - φ - - S 2 ψ φ - S 2 φ * ,

or:

(63)
d 1 ψ - + h ω ψ ' - = A 1 cos ψ - - φ - + A 2 sin ψ - - φ - + A 3 ,
r φ - + h φ ω φ ' - = A 1 * cos ψ - - φ - + A 2 * sin ψ - - φ - + A 3 * ,

or by taking into account the equation (59) it is obtained:

(64)
d 1 h ω δ = A 1 cos ψ - - φ - + A 2 sin ψ - - φ - + A 3 = 0 ,
h φ ω δ = A 1 * cos ψ - - φ - + A 2 * sin ψ - - φ - + A 3 * = 0 .

From the equations (64) δ and ψ--φ- are found. On the right side by assuming δ=1 in the coefficients Aj,Aj* it is obtained:

(65)
h h φ A 1 * - A 1 cos ψ - - φ - + h h φ A 2 * - A 2 sin ψ - - φ - + h h φ A 3 * - A 3 = 0 .

From the equation ψ--φ- is found. One of them can be freely assigned any value, because the system of differential equations is autonomous.

Further investigation is performed in the same way as in case 1. The only difference here is in the fact that four stationary regimes at 0<δ<1 can not exist. In this case only two stationary regimes exist, one of them is stable and one of them is unstable. Also even one stable and one unstable regime may not exist.

When the output system moves at the velocity of the exciting wave, then the obtained results are of the same type as for the case 1.

4. Conclusions

The main objective of this paper is to analyze transverse and longitudinal deformations of the working profile of the input member which displaces the output member, consisting from rigid members according to one coordinate. More detailed analysis is performed when the working profile of the input member is moving according to traveling harmonic Rayleigh waves, and one member of the output system has point contact with the profile of the input member.

Approximate analytical modified asymptotic method is used for the investigation of the stationary motions of the output system, which is based on the division of motion into the slow and quick motions. The characteristics of motion have been obtained, which are valid for the regimes when the output system moves at a much lower average velocity than the velocity of the exciting wave.

Bifurcations and complicated motions and control of the system have been determined by using analytical and numerical methods.

References

  1. Chou C. F., Tegenfeldt J. O., Bakajin O., Chan S. S., Cox E. C., Darnton N., Duke T., Austin R. H. Electrodeless dielectrophoresis of single- and double-stranded DNA. Biophys. J., Vol. 83, 2002, p. 2170-2179. [Search CrossRef]
  2. Cui L., Korgan H. Design and fabrication of travelling wave dielectrophoresis structures. J. Micromech. Microeng., Vol. 10, 2000, p. 72-79. [Search CrossRef]
  3. Talary M. S., Burt J. P. H., Tame J. A., Pethig R. Electromanipulation and separation of cells using travelling electric fields. J. Phys. D Appl. Phys., Vol. 29, 1996, p. 2198-2203. [Search CrossRef]
  4. Dusiaud A. D., Khusid B., Acrivos B. A. Particle segregation in suspensions subject to high gradient ac electric fields. J. Appl. Phys., Vol. 88, 2000, p. 5463-5473. [Search CrossRef]
  5. Wael N. H., Ribberink J. S. Transport processes of uniform and mixed sands in oscillatory sheet flow. Coastal Eng., Vol. 52, No. 9, 2005, p. 745-770. [Search CrossRef]
  6. Wang S. D., Shen Y. M., Zheng Y. H. Two-dimensional numerical simulation for transport and fate of oil spills in seas. Ocean Eng., Vol. 32, 2005, p. 1556-1571. [Search CrossRef]
  7. Mracek M., Wallaschek J. A system for powder transport based on piezoelectrically excited ultrasonic progressive waves. Mater. Chem. Phys., Vol. 90, No. 2-3, 2005, p. 378-380. [Search CrossRef]
  8. Ragulskis M., Sanjuan M. A. F. Transport of particles by surface waves: a modification of the classical bouncer model. New Journal of Physics, Vol. 10, 2008, article no. 083017. [Search CrossRef]
  9. Ragulskis K. M. Mechanisms on a Vibrating Basis (Dynamics and Stability). Kaunas, Institute of the Lithuanian Academy of Sciences, 1963, 232 p., (in Russian). [Search CrossRef]
  10. Ganiev R. F., Ukrainsky L. E. Non-Linear Wave Mechanics & Technologies. Moscow: R&C Dynamics, 2008, 712 p. [Search CrossRef]
  11. Nonlinear Problems of Theory of Vibrations and Control Theory. Vibration Mechanics. For the 80-th Anniversary from the Birth of I. I. Blekhman. Edited by Beleckiy V. V. and others, St. Petersburg, Nauka, 2009, 528 p., (in Russian). [Search CrossRef]
  12. Ragulskis K., Spruogis B., Ragulskis M. Transformation of Rotational Motion by Inertia Couplings. Vilnius, Technika, 1999, 236 p. [Search CrossRef]
  13. Ragulskis M., Sakyte E. An explicit equation for the dynamics of a particle conveyed by a propagating wave. Mathematical and Computer Modeling of Dynamical Systems, Vol. 15, No. 4, 2009, p. 395-405. [Search CrossRef]