• 沒有找到結果。

Analytical study of the stability of a two-ball automatic balancer

N/A
N/A
Protected

Academic year: 2021

Share "Analytical study of the stability of a two-ball automatic balancer"

Copied!
13
0
0

加載中.... (立即查看全文)

全文

(1)

Analytical study of the stability of a two-ball automatic balancer

Chung-Jen Lu



, Ming-Cheng Wang, Shih-Hsuan Huang

Department of Mechanical Engineering, National Taiwan University, No. 1, Roosevelt Road, Section 4, Taipei 10617, Taiwan, ROC

a r t i c l e

i n f o

Article history:

Received 22 August 2007 Received in revised form 23 June 2008 Accepted 26 June 2008 Keywords: Automatic balancer Stability Vibration Optical disk

a b s t r a c t

A ball-type automatic balancer, which consists of several balls moving on a circular orbit, can effectively reduce the imbalance vibrations of an optical-disk drive. Under proper working conditions, the balls can counter-balance the imbalance of a disk by positioning the appropriate angles relative to the mass center of the disk. These particular positions are referred to as the perfect balancing positions. The proper working conditions depend on the stability of the perfect balancing positions. Much work has been done to determine the stability of the equilibrium positions numerically for some specific values of the system parameters. This paper, in contrast, analytically studies the stability characteristics of the equilibrium positions under several conditions. A theoretical model of an optical-disk drive packed with a two-ball automatic balancer is constructed first. The governing equations of the theoretical model are derived using Lagrange’s equations. The closed-form formulae and the existence conditions for the equilibrium positions are presented. Experiments are conducted to verify the theoretical model. The stability of each equilibrium position is determined by Routh’s criterion. Finally, general guidelines on the stability of the equilibrium positions are proposed.

&2008 Elsevier Ltd. All rights reserved.

1. Introduction

The eccentric mass of an optical disk can cause serious vibrations at high rotating speeds. An effective way to reduce the imbalance vibrations is to install a ball-type automatic balancer to the optical-disk drive. A ball-type automatic balancer consists of several balls moving on a circular orbit (Fig. 1). In order to have a clear picture of the working principle of the automatic balancer, consider an undamped rotor system. Assuming that the total mass of the balls is much less than that of the rotor, which includes the disk and the spindle, the behavior of the system is similar to that of a Jeffcott rotor[1]. Because of imbalance, the geometric center C of the disk whirls with the same rotational speed as the disk. The configurations of the system when the rotational speed is subcritical (less than the natural frequency of the system) and supercritical (above the natural frequency) are depicted inFigs. 2a and b, respectively. Note that the center of mass G tends to stay away from the fixed point O in the subcritical speed range while it approaches O in the supercritical speed range. On the other hand, the balls, which can move freely along the orbit, may tend to stay away from O irrespective of the rotational speed. Furthermore, under proper working conditions, the balls will move to specific equilibrium positions where the imbalance vibrations are totally suppressed (Fig. 2c). These equilibrium positions are called the perfect balancing positions and their stability characteristics dictate the proper working conditions.

Contents lists available atScienceDirect

journal homepage:www.elsevier.com/locate/jnlabr/ymssp

Mechanical Systems and Signal Processing

0888-3270/$ - see front matter & 2008 Elsevier Ltd. All rights reserved. doi:10.1016/j.ymssp.2008.06.008

Corresponding author. Tel:. +886 2 3362704; fax: +882 2 23631755. E-mail address:cjlu@ntu.edu.tw (C.-J. Lu).

(2)

Many researchers have studied the properties of the ball-type automatic balancer by either numerical or perturbation methods. Alexander [2] and Cade [3] analyzed the ball-type automatic balancer, but did not describe the working mechanism of the automatic balancer. Bo¨vik and Ho¨gfors[4]examined the stability of the perfect balancing positions of a two-ball balancer by the method of multiple scales. Rajalingham et al.[5]investigated the nonlinear system that consists of an undamped rotor and a single-ball automatic balancer. The stability of the equilibrium position was determined by the eigenvalues of the corresponding linearized system. Lee and Moorham[6]numerically determined the stable regions of the perfect balancing positions by the Floquet theory. Hwang and Chung[7]studied the dynamic characteristics of a two-ball balancer with double races. Chung and Ro[8]and Kang et al.[9]determined the stable regions of the perfect balancing positions of a two-ball automatic balancer in the design parameter space. The stability was checked with the variations for a pair of design parameters. Huang et al. [10]studied the dynamic characteristics of a one-ball automatic balancer

Ball

Disk Spindle

Orbit

Fig. 1. Schematic diagram of the automatic balancer and rotating disk.

C O G G C O G C O trace of C trace of C

Fig. 2. Schematics of the configurations of the undamped disk-balancer system relative to a fixed reference frame for (a) subcritical speeds, (b) supercritical speeds, and (c) perfect balancing.

(3)

considering the effects of rolling resistance and orbit eccentricity. Chung and Jang [11] considered the effects of the flexibility of the shaft on the performance of the ball-type automatic balancer. The stability of the perfect balancing positions of a two-ball balancer was analyzed by the perturbation method. Kim et al.[12]developed a three-dimensional model and examined the dynamic behavior of a two-ball balancer. Chao et al.[13]evaluated the performance of a two-ball balancer using the method of multiple scales considering the in-plane rotational motions. Rajalingham and Bhat [14] considered the application of a two-ball balancer to a vertical flexible rotor. They pointed out that the system may settle to either a balanced state or a limit cycle. Green et al.[15]conducted a nonlinear bifurcation analysis of a two-ball balancer for eccentric rotors.

These works have indeed furthered our understanding of the dynamic characteristics of the ball-type automatic balancer. However, it is difficult to draw general guidelines on the stability of the equilibrium positions purely based on numerical analyses. Lu[16]analytically studied the stability of the equilibrium positions of a one-ball balancer. He showed that two among the five different types of equilibrium positions are always unstable. Moreover, there is at most one stable equilibrium position at any rotational speed. In this study, we extend Lu’s work to a two-ball balancer. The stability of each equilibrium position of a rotating disk equipped with a two-ball automatic balancer is studied analytically. First, the nonlinear governing equations are derived from Lagrange’s equations with respect to a co-rotating coordinate system. Then equilibrium positions and the associated linearized equations are obtained. The closed-form formulae and the existence conditions for the equilibrium positions are presented. Experiments are conducted to validate the mathematical model. Finally, the stability of each equilibrium position is determined analytically by Routh’s criterion and the results are summarized.

2. Mathematical model and governing equations

Fig. 3depicts schematically the system considered, which consists of a rotating unbalanced disk, an automatic balancer, and the suspension system. The disk with mass mdrotates with a constant angular velocity

o

. The mass center G of the disk is located at a distance e from the geometric center C. As a result, the imbalance of the disk is mde. The automatic balancer is composed of an annular groove containing two balls and a damping fluid. The balls move freely along the groove and their movements are subject to viscous damping only. Let d denote the radius of the groove, cbthe viscous damping constant of the damping fluid, and mbthe mass of the ball. The product mbd is the amount of imbalance caused by a single ball and is referred to as the imbalance of a single ball hereafter. The suspension system is isotropic and characterized by equivalent linear springs with stiffness k and viscous dampers with damping constant c.

The xy-reference frame rotates about a fixed point O with the same speed as the disk. When the disk is at rest, its geometric center C coincides with O at which the supporting springs are undeformed. While rotating, the disk undergoes a translatory motion in this rotating frame and its position is defined by the coordinates ð^x; ^yÞ of the geometric center C. The position of the i-th ball is given by angle

b

irelative to the mass center G direction. The reason for the choice of this rotating frame is that the equations of motion expressed in this frame are autonomous.

C X y Y G e d c c k k O θ r x  = t 1 2 

(4)

The equations of motion are derived from Lagrange’s equations given by d dt

q

T

q

_qk   

q

T

q

qk þ

q

V

q

qk þ

q

R

q

_qk ¼0 (1)

where T is the kinetic energy, V the potential energy, R the Rayleigh’s dissipation function, and qk the generalized coordinates. The kinetic energy can be expressed as

T ¼1 2J

o

2 þ1 2M½_^x 2 þ _^y2þ2ð^x_^y  _^x ^yÞ

o

þ ð ^x2þ ^y2Þ

o

2  þ1 2mdð2^xe

o

2 þe2

o

2 þ2e_^y

o

Þ þ1 2mb X2 i¼1

½d2ð _

b

o

Þ2þ2dð _^y þ

o

^xÞð _

b

o

Þcos

b

iþ2dð

o

^y  _^xÞð _

b

o

Þsin

b

i (2)

where J is the moment of inertia of the equivalent rotor and M the total mass of the system. The potential energy is given by V ¼1

2kð^x

2

þ ^y2Þ (3)

Rayleigh’s dissipation function can be represented as R ¼1 2c½_^x 2 þ _^y2þ2ð^x_^y  _^x^yÞ

o

þ ð^x2þ ^y2Þ

o

2  þ1 2cbd 2X2 i¼1 _

b

2i (4)

Substitution of Eqs. (2)–(4) into Eq. (1) yields the nonlinear equations of motion as follows: M €^x  mbdP

2 i¼1

½ð _

b

o

Þ2cos

b

iþ €

b

isin

b

i mde

o

2þc _^x  2M

o

_^y þ ðk  M

o

2Þ^x  c

o

^y ¼ 0

M €^y  mbdP 2 i¼1

½ð _

b

o

Þ 2

sin

b

i €

b

icos

b

i þc _^y þ 2M

o

_^x þ ðk  M

o

2Þ ^y þ c

o

^x ¼ 0

mbd½ð€^y þ 2

o

_^x 

o

2^yÞ cos

b

i ð€^x  2

o

_^y 

o

2^xÞ sin

b

iþd €

b

i þcbd 2_

b

i¼0 i ¼ 1; 2 8 > > > > > > > < > > > > > > > : (5)

In order to simplify the subsequent analysis, the following dimensionless variables are introduced: x ¼^x d; y ¼ ^y d;

t

¼

o

nt c ¼ 2M

Bo

n; cb¼2mb

B

b

o

n;

o

n¼ ffiffiffiffiffi k M r

O

¼

o

o

n ;

m

b¼ mb M;

Z

¼ mbd mde (6) It is worth noting that

O

indicates the ratio of the rotational speed to the natural frequency of the system and

Z

the ratio of the imbalance of a single ball to that of the disk. Substituting Eqs. (6) into Eqs. (5) yields the dimensionless equations of motion as follows: x00

m

bP 2 i¼1 ½ð

b

0 iþ

O

Þ2cos

b

b

00 i sin

b

i þ2

B

x02

O

y0þ ð1 

O

2 Þx  2

B

O

y ¼

m

b

O

2

Z

y00

m

bP 2 i¼1 ½ð

b

0 iþ

O

Þ 2sin

b

i

b

00 icos

b

i þ2

B

y0þ2

O

x0þ ð1 

O

2 Þy þ 2

B

O

x ¼ 0

m

b½ðy00þ2

O

x0

O

2 yÞ cos

b

i ðx002

O

y0

O

2 xÞ sin

b

b

00 i þ2

B

b

m

b

b

0 i¼0; i ¼ 1; 2 8 > > > > > > > < > > > > > > > : (7)

where ð Þ0indicates differentiation with respect to

t

. Let q ¼ ½x; y;

b

1;

b

2T, Eqs. (7) can be put in the matrix form as

Mq00 þDq0 þg þ f ¼ 0 (8) where M ¼ 1 0 

m

bsin

b

1 

m

bsin

b

2 0 1

m

bcos

b

1

m

bcos

b

2 

m

bsin

b

1

m

bcos

b

1

m

b 0 

m

bsin

b

2

m

bcos

b

2 0

m

b 2 6 6 6 6 4 3 7 7 7 7 5 D ¼ 2

B

2

O

2

m

b

O

cos

b

1 2

m

b

O

cos

b

2 2

O

2

B

2

m

b

O

sin

b

1 2

m

b

O

sin

b

2 2

m

b

O

cos

b

1 2

m

b

O

sin

b

1 2

m

b

B

b 0 2

m

b

O

cos

b

2 2

m

b

O

sin

b

2 0 2

m

b

B

b 2 6 6 6 6 4 3 7 7 7 7 5

(5)

g ¼ 

m

b P2 i¼1 ð

b

0 iÞ2cos

b

i 

m

b P2 i¼1 ð

b

0 iÞ2sin

b

i 0 0 2 6 6 6 6 6 6 6 6 4 3 7 7 7 7 7 7 7 7 5 f ¼ ð1 

O

2Þx  2

B

O

y 

m

b

O

2 ðcos

b

1þcos

b

2Þ 

m

b

O

2

Z

2

B

O

x þ ð1 

O

2Þy 

m

b

O

2 ðsin

b

1þsin

b

m

b

O

2 ðx sin

b

1y cos

b

m

b

O

2 ðx sin

b

2y cos

b

2Þ 2 6 6 6 6 6 6 6 4 3 7 7 7 7 7 7 7 5 3. Equilibrium positions

The first step in analyzing a nonlinear system is to identify the equilibrium positions. The equilibrium positions, denoted by ¯q ¼ ½ ¯x; ¯y;

b

¯1;

b

¯2T, can be deduced from Eqs. (7) by suppressing the time derivatives and solving the resulting

algebraic equations ð1 

O

2Þ¯x  2

B

O

¯y 

m

b

O

2 ðcos ¯

b

1þcos ¯

b

2Þ 

m

b

O

2

Z

¼0 2

B

O

¯x þ ð1 

O

2Þ¯y 

m

b

O

2 ðsin ¯

b

1þsin ¯

b

2Þ ¼0

m

b

O

2

ð¯x sin ¯

b

i¯y cos ¯

b

iÞ ¼0; i ¼ 1; 2

8 > > > > < > > > > : (9)

To solve Eqs. (9), we introduce the polar coordinates ¯x ¼ ¯r cos ¯

y

; ¯y ¼ ¯r sin ¯

y

.

Here ¯r denotes the residual vibration of the rotating disk and ¯

y

the angle of OC relative to the x-axis (Fig. 3). Then Eqs. (9) can be rewritten as ð1 

O

2Þ¯r cos ¯

y

2

B

O

¯r sin ¯

y



m

b

O

2 ðcos ¯

b

1þcos ¯

b

2Þ 

m

b

O

2

Z

¼0 2

B

O

¯r cos ¯

y

þ ð1 

O

2 Þ¯r sin ¯

y



m

b

O

2ðsin ¯

b

1þsin ¯

b

2Þ ¼0

m

b

O

2 ¯r sinð ¯

b

i ¯

y

Þ ¼0; i ¼ 1; 2 8 > > > > < > > > > : (10a2 c)

It is worth noting that

z

bdoes not appear in the above equations and hence has no effects on the equilibrium positions. Eqs. (10c) require that sinð ¯

b

i ¯

y

Þ ¼0 or ¯r ¼ 0. The equilibrium positions associated with sinð ¯

b

i ¯

y

Þ ¼0 and ¯r ¼ 0 are called

the unbalanced positions and the perfect balancing positions, respectively. The perfect balancing positions are denoted by P0, where the subscript 0 indicates zero residual vibration. If the balls are indistinguishable, the unbalanced solutions can

be further divided into three cases: (i) ¯

b

1¼ ¯

b

2¼ ¯

y

, (ii) ¯

b

1¼ ¯

b

2¼ ¯

y

þ

p

, and (iii) ¯

b

1¼ ¯

y

and ¯

b

2¼ ¯

y

þ

p

. Including the perfect

balancing position (iv) ¯r ¼ 0, there are totally four different cases of solutions to Eqs. (10). Cases (i) and (ii) are equivalent to a single ball with mass 2mblocated at ¯

b

¼ ¯

y

and ¯

b

¼ ¯

y

þ

p

, respectively. For simplicity, the equilibrium positions associated with (i) and (ii) are referred to as PS1and PS2, respectively. In case (iii), since the two balls counter-balance each other, the resulting system is equivalent to that with no balancer. The equilibrium positions associated with case (iii) are denoted by PNO. Equilibrium positions associated with these four cases are discussed separately below.

(i) PS1( ¯

b

1¼ ¯

b

2¼ ¯

y

)

Substituting ¯

b

1¼ ¯

b

2¼ ¯

y

into Eqs. (10a) and (10b) gives, after some rearrangements,

ð

O

21Þ¯r þ 2

m

b

O

2 ¼ 

m

b

O

2 cos ¯

y

Z

2

B

O

¯r ¼ 

m

b

O

2 sin ¯

y

Z

8 > > > > < > > > > : (11a,b)

By squaring and adding we obtain the two possible equilibrium positions, denoted by P

S1. The corresponding radial

vibrations are denoted by ¯r

S1and expressed as ¯rþ S1 ¯r S1 ) ¼R 2

Z

p pffiffiffiffiD (12)

(6)

where R ¼

O

2

m

b

Z

ðp2þq2Þ and D ¼ p 2þq24

Z

2q2 (13)

in which p and q are functions of

O

as defined below

p ¼

O

21 and q ¼ 2

B

O

(14)

The corresponding ¯

y

0s are determined by ¯

y

S1¼tan1 q  ¯r S1 p  ¯r S12

m

b

O

2 (15) Note that Eq. (11b) implies that sin ¯

y

p0 and hence 

p

o ¯

y

p0.

Before proceeding further, we need to derive the conditions under which ¯r is real and non-negative. As can be seen from Eq. (12), ¯r is real if and only if D is non-negative. From Eq. (13), it follows that:

DX0 if and only if

Z

p

Z

cr (16) where

Z

cr¼ 1 2 ffiffiffiffiffiffiffiffiffiffiffiffiffiffi 1 þp 2 q2 s ¼1 2 ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi 1 þð

O

2 1Þ2 4

z

2

O

2 s (17) It is easy to see that

Z

crX1=2 and the minimum value 1/2 is reached at

O

¼1. Moreover,

Z

cris monotonically decreasing for 0o

O

o1 and monotonically increasing for 1o

O

. Assuming that ¯r is real, i.e.

Z

p

Z

cr, we proceed to derive the conditions for

non-negative ¯r. Consider the case ¯rþ

S1first. As can be seen from Eqs. (12) and (14), in the subcritical speed range,

O

o1, p is

negative, which in turn implies that ¯rþ

S1 is positive. In the supercritical speed range,

O

41, ¯r þ

S1 is non-negative if

2

Z

p þpffiffiffiffiDX0, or equivalently

Z

p1=2. Similarly, we can derive the conditions under which ¯rS1 is non-negative. In summary, ¯rþ S1X0 if

Z

p

Z

cr for 0o

O

o1

Z

p1 2 for 1p

O

( (18) ¯r S1X0 if 1 2p

Z

p

Z

cr for 0o

O

o1 (19) (ii) PS2( ¯

b

1¼ ¯

b

2¼ ¯

y

þ

p

)

Substitution of ¯

b

1¼ ¯

b

2¼ ¯

y

þ

p

into Eqs. (10a) and (10b) yields two possible equilibrium positions, denoted by PS2. The

corresponding amplitudes of radial vibrations are denoted by ¯r

S2and expressed as rþ S2 r S2 ) ¼R 2

Z

p pffiffiffiffiD (20)

where R and D are defined by Eq. (13). The corresponding ¯

y

0s are determined by ¯

y

S2¼tan1 q  ¯r S2 p  ¯r S2þ2

m

b

O

2 (21)

By a similar process as before, we can obtain the conditions under which ¯r

S2are non-negative. The results are summarized

below rþ S2Xif 0

Z

p1 2 for 0o

O

o1

Z

p

Z

cr for 1p

O

( (22) ¯r S2X0 if 1 2p

Z

p

Z

cr for 1p

O

(23) (iii) PNOð ¯

b

1¼ ¯

y

; ¯

b

2¼ ¯

y

þ

p

Þ

In this case, the two balls are located at the two ends of a diameter of the circular orbit and counter-balance each other. Substituting ¯

b

1¼ ¯

y

and ¯

b

2¼ ¯

y

þ

p

into Eqs. (10a) and (10b) gives, after some rearrangements

ð

O

21Þ¯r cos ¯

y

þ2

B

O

¯r sin ¯

y

¼

m

b

O

2

Z

ð

O

21Þ¯r sin ¯

y

2

B

O

¯r cos ¯

y

¼0 8 > < > :

(7)

By squaring and adding we obtain the associated amplitude of radial vibration ¯rNOas ¯rNO¼R ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi p2þq2 q (24) where R is defined by Eq. (13). The corresponding ¯

y

NOis determined by

¯

y

NO¼tan1 q p Unlike P S1 and P  S2, PNOexists unconditionally. (iv) P0ð¯r ¼ 0Þ

In this case, the system is perfectly balanced and has no residual vibration. Substituting ¯r ¼ 0 into Eqs. (10) yields ðcos ¯

b

1þcos ¯

b

2Þ ¼  1

Z

ðsin ¯

b

1þsin ¯

b

2Þ ¼0 8 > < > :

which has the solution ¯

b

1¼cos1  1 2

Z

  and

b

¯2¼  ¯

b

1 (25)

The results indicate that, when the system is perfectly balanced, the two balls are located symmetrically on the two sides of a diameter passing the center of mass G of the disk. As can be seen from Eq. (25), ¯

b

iexists only if

Z

X(1/2). Consequently, P0 exists only if

Z

X(1/2).

The above analysis indicates that

Z

,

O

and

z

are major factors affecting the equilibrium positions of the system. The existence conditions of an equilibrium position can be seen clearly by identifying its existence regions in the

O



Z

plane. Fig. 4illustrates the existence regions of the equilibrium positions in the

O



Z

plane for

z

¼0.1. The

O



Z

plane is divided into four regions by the curve

Z

crand the line

Z

¼0.5. Note that the value of

Z

cr, as defined by Eq. (17), varies with

z

, but the qualitative characteristics of the existence regions remain the same. The region below the line

Z

¼0.5 is denoted as region I. In this region, perfect balancing cannot be achieved and only Pþ

S1, P þ

S2, and PNOmay exist. The region above

Z

¼0.5 is further divided into three regions by the two branches of

Z

cr. The region between the two branches of

Z

cris indicated as region III. The regions to the left and right of region III are called regions II and IV, respectively. In region II, only P

S1, P0, and PNOmay exist; in region III, only PNOand P0may exist; in region IV, only PS2, PNO, and P0may exist. The equilibrium positions that may exist in each region are listed in the figure. As can be seen from the figure, for a given imbalance ratio

Z

, at most four equilibrium positions may exist at a fixed rotational speed. The equilibrium positions enclosed by & are unstable in that region. The stability of the equilibrium positions will be studied, experimentally and analytically, later. Note that region III is divided into two parts by the line

O

¼1 to clearly show the stable region of P0. For the sake of clearness, the configurations of the system associated with each equilibrium position in different

O



Z

regions are also shown inTable 1.

4. Experimental study

Experimental study is performed to investigate the fundamental properties of the ABS and validate the mathematical model.Fig. 5shows a photograph of the experimental apparatus that consists of an ABS, a clear plastic disk, a DC motor, and

0 2 0 1 2 3 Ω PS1 PS1 PNo PNo PNo PNo PNo P0 P0 P0 P0 (I) (II) (III) (IV)  1 3 + -PS2 + PS2 + PS1 + cr

(8)

the suspension base. The ABS and the disk were attached firmly to the motor and then mounted on the suspension base. The suspension base was built with four L-shaped beams to realize isotropic stiffness and damping. Modal analysis was applied to determine the natural frequencies of the system. The first two natural frequencies were found to be 25 and 68 Hz. A finite element model of the experimental apparatus was constructed to identify the vibration modes of the system at these two frequencies. The results indicate that the first natural frequency, which is a repeated natural frequency, corresponds to two translational modes in the x- and y-directions, respectively. On the other hand, the second natural frequency corresponds to a torsion mode about the spindle. When the system is vibrating at low frequencies where the motion of the system is primarily composed of the two translational modes, the L-shaped beams behave like linear springs directly connected to the geometric center of the disk, which is exactly the case considered in the mathema-tical model. Therefore, the first measured natural frequency, 25 Hz, is identified as the natural frequency

o

n of the mathematical model and was used to normalize the rotational speed henceforth. The experimental procedure is briefly described below.

The rotor system was balanced prior to experiments. After that, a bolt was locked on the disk at a radial distance 25 mm to represent the eccentric mass of the disk. Bolts of different sizes were used to adjust the value of the imbalance ratio

Z

. A laser Doppler vibrometer was used to measure the vibration of one of the L-shaped beams. Video and still cameras were used in conjunction with a strobe light to record the angular positions of the balls. While conducting the experiment, the

Table 1

Configurations and existence regions of the equilibrium positions

—rotating coordinate system. —center of mass of the disk. —geometric center of the disk.

—unstable equilibrium position.

(9)

motor was powered to accelerate the disk up to desired speeds. The vibration levels of the L-shaped beam and the ball positions were recorded.

The case

Z

o0.5 was considered first. Fig. 6(a) shows the levels of steady-state vibration versus rotational speed with and without the ABS for

Z

¼0.285. The vertical dotted lines indicate the first two natural frequencies of the system. No data could be obtained in the vicinity of the natural frequencies due to the large vibration levels. At sub-critical rotational speeds,

O

o1, the two balls stick together and stay close to the eccentric mass.Fig. 6(b) is a typical result showing the positions of the balls at

O

¼0.84. In this rotational speed range, the vibration level is increased by the presence of the ABS. On the other hand, at supercritical rotational speeds,

O

41, the two balls still stick together but stay on the other side of the eccentric mass. Fig. 6(c) shows the positions of the balls at

O

¼1.56. In this case, the ABS partially counterbalances the eccentric mass and hence reduces the steady-state vibration level. Table 1indicates that, for

Z

o0.5, there exist three equilibrium positions Pþ

S1, PþS2, and PNO. ComparingFigs. 6(b) and (c)

with the configurations of different equilibrium positions shown in Table 1, one can see that the equilibrium positions observed in the experiment are Pþ

S1. These results imply that, for

Z

o0.5, P þ

S1is stable but P þ

S2 and PNOmay be unstable.

Then we examined the case

Z

40.5. A larger bolt was used to adjust the value of imbalance ratio

Z

to be 0.832. The associated steady-state vibration levels at different rotational speeds with and without the ABS are shown inFig. 7(a). At subcritical speeds, the balls stick together and stay close to the eccentric mass, as can be seen inFig. 7(b). In this case, the ABS increases the vibration levels. On the other hand, in the speed range 1.4o

O

o2, the two balls stay on the other side of the eccentric mass and open a suitable angle (Fig. 7(c)) such that the disk is almost completely counterbalanced. As a result, the vibration levels are reduced significantly. Comparison with the equilibrium configurations shown inTable 1indicates that the equilibrium positions observed in the experiment are Pþ

S1 and P0 for subcritical and supercritical speeds, respectively. These results imply that P0is stable in some supercritical speed range but unstable for subcritical speeds. The stability of the equilibrium positions is studied analytically in the next section.

5. Stability analysis

For an autonomous system, the stability of an equilibrium position can be determined by the eigenvalues of the associated linearized equations. The eigenvalues are roots of the corresponding characteristic polynomial as follows:

FðsÞ ¼ a0snþa1sn1þa2sn2þ    þan1s þ an (26)

where the coefficients ai, which are functions of the equilibrium position and system parameters, are given explicitly in the Appendix. According to Routh’s criterion, information about the stability of the equilibrium position can be obtained directly from the coefficients of the characteristic polynomial[17]. A necessary condition for the equilibrium position to be stable is that all the coefficients of the characteristic polynomial are of the same sign. As can be seen from Eqs. (A.1) and (A.2) in Appendix, both a0and a1are positive. Consequently, the equilibrium position is unstable if any of the coefficients

Motor Eccentric Mass ABS Disk Suspension Base

(10)

from a2to a8is negative. This criterion is used in the following paragraphs to analytically show that some equilibrium positions are unstable under certain conditions. The results are listed below.

(S1) The equilibrium position P

S1 is unstable.

Substitution of ¯

b

1¼ ¯

b

2¼ ¯

y

and ¯rS1 (as defined by Eq. (12)) for ¯r into Eq. (A.9) gives

a8¼ 

m

b

O

6pffiffiffiffiD

Z

! ¯r S1o0 Therefore P S1is unstable.

(S2) The equilibrium position P

S2 is unstable.

Substitution of ¯

b

1¼ ¯

b

2¼ ¯

y

þ

p

and ¯rS2(as defined by Eq. (20)) for ¯r into Eq. (A.9) gives

a8¼ 

m

b

O

6 ffiffiffiffi D p

Z

! ¯r S2o0 Therefore P S2is unstable.

(S3) The equilibrium position Pþ

S2 is unstable.

In this case, the two balls are located at the same position ¯

b

1¼ ¯

b

2¼ ¯

y

þ

p

. As a result, this particular equilibrium

position of the two-ball system is the same as that of a one-ball system with a ball of mass 2mblocated at ¯

b

¼ ¯

y

þ

p

. In addition, if the two balls are constrained to move together, then the dynamic characteristics of the two-ball system will be exactly the same as that of the one-ball system. Since this type of equilibrium position for a one-ball system

0 30 40 50 60 70 80 0 400 V ibration V elocity (mm/s) with ABS (η = 0.285) Unbalanced Disk 0 2 Ω 300 200 100 10 20 Rotational Speed (Hz) 1 3

(11)

has been shown to be unstable [16], we can conclude that, for a two-ball system, the equilibrium position Pþ S2 is

unstable.

(S4) The equilibrium position PNOis unstable.

Substituting of ¯

b

1¼ ¯

y

, ¯

b

2¼ ¯

y

þ

p

and ¯rNOfor ¯r into Eq. (A.9) gives a8¼ ¯r2NO

O

4

ðp2þq2Þo0.

Therefore PNOis unstable.

(S5) The equilibrium position P0is unstable for

O

o1.

Substitution of ¯

b

2¼  ¯

b

1and ¯r ¼ 0 into Eq. (A.8) gives a7¼4

m

b

B

b

O

O

21Þ.

Since a7o0 for

O

o1, P0is unstable for

O

o1.

In order to clearly show the above results, the unstable equilibrium positions inFig. 4are enclosed by &. However, it must be noted that the equilibrium positions that are not enclosed by & are not necessarily stable in their existence regions. To verify the above analytical results, the stable regions of the equilibrium positions in the

O



Z

plane are determined numerically by the eigenvalues of the corresponding linear systems.Fig. 8shows a typical result for

z

¼0.1 and

z

b¼0.5. The light gray and dark gray areas are stable regions of PþS1and P0, respectively. No stable equilibrium positions exist in the white area. ComparingFigs. 4 and 8, one can see that the stable region of P0is much smaller than its existence region. Specifically, P0is unstable in some areas inside its existence region; the exact size and shape of the stable region of P0depend on the values of the system parameters. As a result, to achieve perfect balancing, one needs to properly set the values of the system parameters such that P0is stable.

0 30 40 50 60 70 80 0 V ibration V elocity (mm/s) with ABS (η = 0.832) Unbalanced Disk 0 2 Ω 200 150 50 100 10 20 Rotational Speed (Hz) 1 3

(12)

6. Conclusion

Ball-type automatic balancers have been used in the Optical-disk drive industry to suppress rotational vibrations induced by an unbalanced optical disk. Proper conditions under which the imbalance vibrations can be totally suppressed are closely related to the stability of the perfect balancing positions. We have investigated the stability property of a system representative of an optical-disk drive equipped with a two-ball automatic balancer and a suspension system. Experiments were performed to verify the theoretical models. Let

Z

denote the ratio of the imbalance of a single ball to that of the disk and

O

the ratio of the rotational speed to the natural frequency of the system. The existence regions of each equilibrium position in the

Z



O

have been identified analytically. General conditions under which some of the equilibrium positions are unstable are proposed. The analytical nature of our investigation has yielded unique understanding of the stability property of the system and insights into the conditions under which perfect balancing can be achieved. The results of the analysis are summarized below.

1. There are totally six different equilibrium positions and four of them are always unstable. 2. There is at most one stable equilibrium position at a rotational speed.

3. The perfect balancing position is unstable for

O

o1.

4. In the areas

O

41 and

Z

4(1/2), the only possible stable equilibrium position is the perfect balancing position. The stable region of the perfect balancing position, whose size and shape vary with the values of system parameters, is only a subset of this area.

Acknowledgement

This work was supported by the National Science Council of R.O.C. under Grant NSC96-2221-E-002-216.

Appendix

The coefficients of the characteristic polynomial are given below: a0¼ ð1  2

m

bÞ þ 1 2

m

2 b½1  cos 2ð ¯

b

1 ¯

b

2Þ (A.1) a1¼4ð1 

m

bÞð

B

þ

B

bÞ (A.2) a2¼ ð1 

m

bÞ½2ð1 þ

O

2 þ4

BB

bÞ þ¯r

O

2ðcosð ¯

b

1 ¯

y

Þ þcosð ¯

b

2 ¯

y

ÞÞ þ4ð

B

þ

B

bÞ2þ2

m

2b

O

2 ½1  cos2ð ¯

b

1 ¯

b

2Þ (A.3) a3¼2f2½4

B

2

B

B

bð2 

m

bÞð1 þ

O

2 Þ þ

B

ð1 þ 4

B

2 bþ ð1 þ 2

m

O

2 Þ

þ¯r

O

B

ð2 

m

bÞ þ

B

bðcosð ¯

b

1 ¯

y

Þ þcosð ¯

b

2 ¯

y

ÞÞg (A.4)

0 0 1 2 3 P0 PS1+ η 1 2 3 Ω

(13)

a4¼ ð

O

21Þ2þ8

B

B

bþ2

B

ð1 þ

BB

bÞ þ4½

B

2þ3

m

bþ4

BB

bð1 þ

m

bÞ þ2

B

2b

O

2 þ2

m

b

O

4 þ3

m

2 b

O

4 ð1  cos 2ð ¯

b

1 ¯

b

2ÞÞ þ¯r2

O

4 cosð ¯

b

1 ¯

y

Þcosð ¯

b

2 ¯

y

Þ þ¯r

O

2½4

B

ð

B

þ2

B

bÞ þ ð2 

m

bÞð1 þ

O

2

Þðcosð ¯

b

1 ¯

y

Þ þcosð ¯

b

2 ¯

y

ÞÞ (A.5)

a5¼4f¯r2

B

O

4cosð ¯

b

1 ¯

y

Þcosð ¯

b

2 ¯

y

Þ þ

B

O

21Þ2þ4

BB

bð1 þ

O

þ2ð2

B

2þ3

m

B

b

O

m

bð3

B

þ

B

O

4þ¯r

O

B

þ ð1 þ 2

B

B

b

þ ð

B

þ

Bm

B

O

2ðcosð ¯

b

1 ¯

y

Þ þcosð ¯

b

2 ¯

y

ÞÞg (A.6)

a6¼24

BB

b

m

b

O

4 þ2

m

b

O

4 ½ð

O

21Þ þ 2

m

b

O

2 sin2ð ¯

b

1 ¯

b

2Þ þ4

B

2b½ð

O

2 1Þ2 þ4

B

2

O

2  þ2¯r2

O

4 ð1 þ

O

2þ2

B

2 Þcos ð ¯

b

1 ¯

y

Þcos ð ¯

b

2 ¯

y

Þ þ¯r

O

2 ½ð

O

21Þ2 þ4

B

2

O

2 þ8

BB

bð1 þ

O

2 Þ þ

m

b

O

2

ð6 þ

O

2Þðcosð ¯

b

1 ¯

y

Þ þcosð ¯

b

2 ¯

y

ÞÞ (A.7)

a7¼2

O

2f2

m

b

B

b

O

O

21Þ þ ¯r½2¯r

B

O

2ð1 þ

O

2Þcosð ¯

b

1 ¯

y

Þcosð ¯

b

2 ¯

y

Þ

þ ½

B

O

21Þ2þ

B

O

2ð4

BB

bþ3

m

b

O

2

Þðcosð ¯

b

1 ¯

y

Þ þcosð ¯

b

2 ¯

y

ÞÞg (A.8)

a8¼

O

4f¯r2½ð

O

21Þ2þ4

B

2

O

2cosð ¯

b

1 ¯

y

Þcosð ¯

b

2 ¯

y

Þ þ

m

b

O

2 ½¯rð

O

21Þðcosð ¯

b

1 ¯

y

Þ þcosð ¯

b

2 ¯

y

ÞÞ þ

m

b

O

2sin2 ð ¯

b

1 ¯

b

2Þg (A.9) References

[1] L. Meirovitch, Elements of Vibration Analysis, McGraw-Hill, New York, 1986.

[2] J.D. Alexander, An automatic dynamic balancer, in: Proceedings for the Second Southeastern Conference, 1964, pp. 415–426. [3] J.W. Cade, Self-compensating balancing in rotating mechanism, Design News (1965) 234–239.

[4] P. Bo¨vik, C. Ho¨gfors, Autobalancing of rotors, Journal of Sound and Vibration 111 (1986) 429–440.

[5] C. Rajalingham, R.B. Bhat, S. Rakheja, Automatic balancing of flexible vertical rotors using a guided ball, International Journal of Mechanical Sciences 40 (1998) 825–834.

[6] J. Lee, W.K.V. Moorham, Analytical and experimental analysis of a self-compensating dynamic balancer in a rotating mechanism, ASME Journal of Dynamic Systems, Measurement and Control 118 (1996) 468–475.

[7] C.-H. Hwang, J. Chung, Dynamic analysis of an automatic ball balancer with double races, JSME International Journal, Series C 42 (1999) 265–272. [8] J. Chung, D.S. Ro, Dynamic analysis of an automatic dynamic balancer for rotating mechanisms, Journal of Sound and Vibration 228 (1999)

1035–1056.

[9] J.R. Kang, C.P. Chao, C.L. Huang, C.K. Sung, The dynamics of a ball-type balancer system equipped with a pair of free-moving balancing masses, ASME Journal of Vibration and Acoustics 123 (2001) 456–465.

[10] W.Y. Huang, C.P. Chao, J.R. Kang, C.K. Sung, The application of ball-type balancers for radial vibration reduction of high-speed optic disk drives, Journal of Sound and Vibration 250 (2002) 415–430.

[11] J. Chung, I. Jang, Dynamic response and stability analysis of an automatic ball balancer for a flexible rotor, Journal of Sound and Vibration 259 (2003) 31–43.

[12] W. Kim, D.-J. Lee, J. Chung, Three-dimensional modelling and dynamic analysis of an automatic ball balancer in an optical disk drive, Journal of Sound and Vibration 285 (2005) 547–569.

[13] P.C.P. Chao, C.-K. Sung, C.-C. Wang, Dynamic analysis of the optical disk drives equipped with an automatic ball balancer with consideration of torsional motions, ASME Journal of Applied Mechanics 72 (2005) 826–842.

[14] C. Rajalingham, R.B. Bhat, Complete balancing of a disk mounted on a vertical cantilever shaft using a two ball automatic balancer, Journal of Sound and Vibration 290 (2006) 169–191.

[15] K. Green, A.R. Champneys, N.J. Lieven, Bifurcation analysis of an automatic dynamic balancing mechanism for eccentric rotors, Journal of Sound and Vibration 291 (2006) 861–881.

[16] C.J. Lu, Stability analysis of a single-ball automatic balancer, ASME Journal of Vibration and Acoustics 128 (2006) 122–125. [17] K. Ogata, Modern Control Engineering, Prentice-Hall, Englewood Cliffs, NJ, 2001.

數據

Fig. 1. Schematic diagram of the automatic balancer and rotating disk.
Fig. 3 depicts schematically the system considered, which consists of a rotating unbalanced disk, an automatic balancer, and the suspension system
Fig. 4 illustrates the existence regions of the equilibrium positions in the O  Z plane for z ¼ 0.1
Table 1 indicates that, for Z o0.5, there exist three equilibrium positions P þ S1 , P þ S2 , and P NO
+4

參考文獻

相關文件

Animal or vegetable fats and oils and their fractiors, boiled, oxidised, dehydrated, sulphurised, blown, polymerised by heat in vacuum or in inert gas or otherwise chemically

Milk and cream, in powder, granule or other solid form, of a fat content, by weight, exceeding 1.5%, not containing added sugar or other sweetening matter.

For 5 to be the precise limit of f(x) as x approaches 3, we must not only be able to bring the difference between f(x) and 5 below each of these three numbers; we must be able

In words, the Product Rule says that the derivative of a product of two functions is the first function times the.. derivative of the second function plus the second function

[This function is named after the electrical engineer Oliver Heaviside (1850–1925) and can be used to describe an electric current that is switched on at time t = 0.] Its graph

Isakov [Isa15] showed that the stability of this inverse problem increases as the frequency increases in the sense that the stability estimate changes from a logarithmic type to

You are given the wavelength and total energy of a light pulse and asked to find the number of photons it

refined generic skills, values education, information literacy, Language across the Curriculum (