Published: 30 June 2013

Wave displacement of rigid bodies and particles

Kazimieras Ragulskis1
Liutauras Ragulskis2
1Kaunas University of Technology, Kaunas, Lithuania
2Vytautas Magnus University, Kaunas, Lithuania
Views 11
Reads 4
Downloads 1213

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.

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
A1u,ξ, A2u,v, B1u+u1,v+v1, B2u-Rsinφ,v-Rcosφ,

where:

2
u=x+η, v=ξ, R=A1B2=A2B2.

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. 1The 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. 2Forces 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+η'xx˙+η't, v˙=ξ'xx˙+ξ't,
6
u¨=1+η'xx¨+ηxxx˙2+2ηxtx˙+ηtt, v¨=ξ'xx¨+ξxxx˙2+2ξxtx˙+ξtt.

Further the following relationships are met:

7
F1=1+η'xu¨+ξ'xv¨=a1x¨+a2x˙2+2a3x˙+a4,
F2=-ξ'xu¨+1+η'xv¨=b2x˙2+2b3x˙+b4,

where:

8
a1=1+η'x2+ξ'x2, a2=1+η'xηxx+ξ'xξxx,
a3=1+η'xηxt+ξ'xξxt, a4=1+η'xηtt+ξ'xξtt,
b2=1+η'xξxx-ξ'xηxx, b3=1+η'xξxt-ξ'xηxt,
b4=1+η'xξtt-ξ'xηtt.

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

9
tanα=ξ'x1+η'x.

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

s˙21=(v˙-ξ't)sinα+(u˙-η't)cosα,

or:

10
s˙21=a1x˙.

2.2. Differential equations of motion

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

11
m+m2u¨+m2RPu-Fu+f1s˙21cosα=-Nsinα+f0sgns˙21cosα,
m+m2v¨+g+m2RPv-Fv+f1s˙21sinα=Ncosα-f0sgns˙21sinα,

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
Pu=-φ¨cosφ+φ˙2sinφ,Pv=φ¨sinφ+φ˙2cosφ,
Fu=Huu˙1+Cuu1, Fv=Hvv˙1+Cvv1,

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+η'xIu+ξ'xIv+-ξ'xIu+1+η'xIvf0sgns˙21+hα1x˙=0,
14
Rφ¨-u¨cosφ+v¨sinφ+gsinφ+Kφm2Rφ˙=0,
15
u¨+u¨1+Fu=0,
16
v¨+v¨1+Fv=0,

where:

17
Iu=u¨+μ2RPu-μ1Fum1,Iv=v¨+g+μ2RPv-μ1Fvm1,
μ1=m1m+m2,μ2=m2m+m2,h=f1m+m1, ˙=ddt,

Kφ is coefficient of viscous friction.

3. Systems excited by traveling waves

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

18
η=ηωt-kx, ξ=ξ(ωt-kx),

where η and ξ are periodic functions of their arguments.

The following notations are introduced:

19
ψ=ωt-kx, τ=ωt, β=kx, β1=ku1, β2=kv1, '=ddτ.

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

20
ku'=1-c4ψ', kv'=c5ψ', ku=-c4ψ+c6ψ'2, kv=c5ψ+c7ψ'2,
kF1=-c1ψ+c2ψ'2, kF2=c3ψ'2,

where:

21
c1=1-kη'2+k2ξ'2, c2=1-kη'kη-k2ξ'ξ, c3=1-kη'kξ+k2ξ'η,
c4=1-kη', c5=kξ', c6=kη, c7=kξ, r=kR.

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

22
-c1ψ+c2ψ'2+c4μ2rpu-μ1fu-c5kgω2+μ2rpv-μ1fv+ -c3ψ'2+c5μ2rpu-μ1fu+c4kgω2+μ2rpv-μ1fvf0sgns˙21+c1hω1-ψ'=0,
23
rφ+(c4ψ-c6ψ'2)cosφ+c5ψ+c7ψ'2sinφ+hφωφ'+kgω2sinφ=0,
24
-c4ψ+c6ψ'2+β1+fu=0,
25
c5ψ+c7ψ'2+β2+fv=0,

where:

26
pu=-φcosφ+φ'2sinφ, pv=φsinφ+φ'2cosφ,
fu=huωβ'1+nuω2β1, fv=hvωβ'2+nvω2β2,
h=f1m+m2, hu=Hum2, hv=Hvm2, nu2=Cum2, nv2=Cvm2, hφ=kKφm2R.

3.1. Case: harmonic traveling waves:

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

27
η=Acosωt-kx, ξ=Bsin(ωt-kx).
28
c1=(1+asinψ)2+b2cos2ψ, c2=-(1+asinψ)acosψ+b2sinψcosψ,
c3=-(1+asinψ)bsinψ-absinψcosψ,
c4=1+asinψ, c5=bcosψ, c6=-acosψ, c7=-bsinψ, a=kA, b=kB.

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,, mn1,

ψ- 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ωψ'~=ψ'-2c2ψ0+μ1rφ'-2c4ψ0sinφ0~-c-4μ1fu
-c5ψ0kgω2+μ2rφ'-2cosφ0~+c1ψ0~hω=0,
32
rφ~+kφωφ'~-ψ'-2c6ψ0cosφ0+φ'-2c7ψ0sinφ0~+kgω2sinφ0=0,
33
-c-4ψ~+ψ'-2c6(ψ0)+β1+fu=0,
34
ψ'-2c7(ψ0)+β2+fv=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ωψ'-=ψ'-2c2ψ0+ψ~-+μ1rφ'-2c4ψ0+ψ~sinφ0+φ~¯-μ1c4ψ0+ψ~fu-
- kgω2c5ψ0+ψ~--μ2rφ'-2c5ψ0+ψ~cosφ0+φ~-+hωc1ψ0+ψ~-=0,
36
rφ-+kφωφ'--ψ'-2c6(ψ0+ψ~)cos(φ0+φ~)-+φ'-2c1(ψ0+ψ~)sin(φ0+φ~)¯
+kgω2sin(φ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.5a2+b2ψ~+hωψ'~=ψ'-2-acosψ0-0.5a2-b2sin2ψ0
-kgω2bcosψ0+hω[2asinψ0-0.5(a2-b2)cos2ψ0]=0,
38
1+0.5a2+b2ψ-+hωψ'-=-ψ'-2acosψ0+ψ~-+0.5a2-b2sin2ψ0+ψ~¯+kgω2bcosψ0+ψ~-
+hω1-0.5a2-b2+2asinψ0+ψ~¯-0.5a2-b2cos2ψ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-a22d1S-δ2+kgω2d22+4hω2δ2+hω2+0.25a21-d222δ4+hω24δ2+hω2,

where:

40
d1=1+0.5a2+b2, d2=ba.

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-a22d1-1+kgω2d22+4hω21+hω2+0.25a21-d2221+hω24+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
d1ψ-+hωψ'-=a-kgω2d2cosψ-+hωd1+2sinψ--0.5a1-d22cos2ψ-=Ψψ-.

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

46
Ψ(ψ-)=0.

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

47
Ψψ-ψ-=akgω2d2sinψ-+hω2cosψ-+a1-d22sin2ψ-<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
d1λ2+d1hωλ-Ψψ-=0,

where:

49
Ψψ-=akgω2sinψ-+hω2cosψ-+a1-d22sin2ψ-.

Roots of the equation (49) are:

50
λ1;2=-12hω±14hω2+d12Ψψ-.

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. 3Transition 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. 4Steady 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
d1ψ+hωψ'=r0+r1cosψ0+s1sinψ0+s1φsinφ0+r2cos2ψ0
+s2sin2ψ0+r2ψφcos(ψ0+φ0),
55
φ+hφωφ'=r0*+s1φ*sinφ0+r2φ*cosψ0+φ0,

where it is denoted:

56
r0=1+0.5ω2+b2hω, r1=-a+kgω2b, s1=2hωa, s1φ=μ2rδ2,
r2=-0.5a2-b2hω, s2=-0.5a2-b2δ2, r2ψφ=-μ1rδ2,
r0*=-0.5a2-b2δ2cosψ--φ-, s1φ*=-hφgω2, r2φ*=0.5a-bδ.

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

57
ψ~=R1cosψ0+S1sinψ0+R1φcosφ0+S1φsinφ0+R2cos2ψ0+S2sin2ψ0
+R2ψφcosψ0+φ0+S2ψφsinψ0+φ0,
58
φ~=R1φ*cosφ0+S1φ*sinφ0+R2φ*cosψ0+φ0+S2φ*sinψ0+φ0,

and the slow motions are:

59
d1ψ-+hωψ'-=r0-r1ψ~sinψ0¯+s1ψ~cosψ0-+s1φφ~cosφ0-
-2r2ψ~sin2ψ0¯+2s2ψ~cos2ψ0--2r2ψφ(ψ~+φ~)sin(2δτ+ψ-+φ-)¯,
60
rφ-+hφωφ'-=r0*+s1φ*φ~cosφ0--r2φ*(ψ~+φ~)sin(2δ+ψ-+φ-)¯.

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

61
d1(ψ-+hωψ'-)=r0-0.5r1S1+0.5s1R1+0.5s1φR1φ*
+r2-S2-R2ψφsinψ--φ--S2ψφcosψ--φ-
+s2R2+R2ψφcosψ--φ--S2ψφsinψ--φ-
+r2ψφR2sinψ--φ-+S2cosψ--φ--S2ψφ-S2φ*,
62
rφ-+hφωφ'-=s0*+0.5s1φ*R1φ*+r2φ*R2sinψ--φ--S2cosψ--φ--S2ψφ-S2φ*,

or:

63
d1ψ-+hωψ'-=A1cosψ--φ-+A2sinψ--φ-+A3,
rφ-+hφωφ'-=A1*cosψ--φ-+A2*sinψ--φ-+A3*,

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

64
d1hωδ=A1cosψ--φ-+A2sinψ--φ-+A3=0,
hφωδ=A1*cosψ--φ-+A2*sinψ--φ-+A3*=0.

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

65
hhφA1*-A1cosψ--φ-+hhφA2*-A2sinψ--φ-+hhφA3*-A3=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

  • 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.
  • Cui L., Korgan H. Design and fabrication of travelling wave dielectrophoresis structures. J. Micromech. Microeng., Vol. 10, 2000, p. 72-79.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • Ragulskis K. M. Mechanisms on a Vibrating Basis (Dynamics and Stability). Kaunas, Institute of the Lithuanian Academy of Sciences, 1963, 232 p., (in Russian).
  • Ganiev R. F., Ukrainsky L. E. Non-Linear Wave Mechanics & Technologies. Moscow: R&C Dynamics, 2008, 712 p.
  • 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).
  • Ragulskis K., Spruogis B., Ragulskis M. Transformation of Rotational Motion by Inertia Couplings. Vilnius, Technika, 1999, 236 p.
  • 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.

About this article

Received
02 April 2013
Accepted
03 June 2013
Published
30 June 2013
Keywords
traveling waves
transient motions
stationary state
stability
bifurcation