EngineeringTribology Paper
EngineeringTribology Paper
EngineeringTribology Paper
net/publication/243972353
CITATIONS READS
13 502
2 authors:
Some of the authors of this publication are also working on these related projects:
All content following this page was uploaded by Pranab Samanta on 03 June 2014.
The manuscript was received on 24 February 2007 and was accepted after revision for publication on 1 August 2007.
DOI: 10.1243/13506501JET282
Abstract: Survey of patents on bearings indicates the maturity of hydrodynamic and rapid devel-
opment of magnetic bearings. Active magnetic bearings are costlier compared with permanent
magnetic bearings. To understand the performance characteristics of permanent magnetic bear-
ings, an experimental setup has been developed. Experimental studies on radial permanent
magnetic bearings demonstrated the drawbacks, such as high axial thrust and low load capacity.
This has led the authors to hybridize the permanent magnet with hydrodynamic technology and
to explore the possibility of achieving the low starting torque of a permanent magnetic bearing
and the medium to high load carrying capacity of a hydrodynamic bearing in a single bearing
arrangement.
Simulation is carried out in order to reduce axial force-effect and enhance the radial force
supported by the permanent magnetic bearing. Results of simulation on permanent magnetic
bearing have been compared with that of published research papers. Finally an algorithm has
been developed to investigate the coupling of forces generated by permanent magnets and
hydrodynamic actions. Results of load sharing have been reported. The experimentally mea-
sured displacements of the shaft running at 500, 2000, and 3000 r/min have been plotted. The
effect of hydrodynamics on shaft orbit has been illustrated.
Keywords: axial force of permanent magnetic bearings, sector-shaped magnetic bearings, test
setup for horizontal journal bearings, modelling of magnetic bearings, simulation of magnetic
bearings
JET282 © IMechE 2007 Proc. IMechE Vol. 221 Part J: J. Engineering Tribology
882 H Hirani and P Samanta
Static Dynamic
Wear friction friction
coefficient coefficient coefficient
Type [1] [2] [2]
Proc. IMechE Vol. 221 Part J: J. Engineering Tribology JET282 © IMechE 2007
Hybrid journal bearings 883
Two non-contact type eddy current proximity been used to record the completion of cycle (shaft rota-
sensors (range 0–2 mm) have been used to measure tion by 360◦ ) that helps ‘time synchronous averaging’.
the film thickness of the bearing. A sleeve has been A photograph of the setup developed is shown in Fig. 3.
mounted on the shaft as a target material for the prox- A loading arrangement with leverage ratio of 1:2, as
imity sensors. The displacement sensors are placed at shown in Fig. 3, is attached to the test setup to apply
a phase angle of 90◦ so that their measurements corre- the static load on the rotor. The load arm contacts the
spond to displacement in the X - and Y -axis directions. rotor through a ball bearing that facilitates application
A fibre optic displacement sensor (range 0–12 mm) has of static load on the rotating rotor.
JET282 © IMechE 2007 Proc. IMechE Vol. 221 Part J: J. Engineering Tribology
884 H Hirani and P Samanta
Before mounting the shaft-rotor system, alignment dial indicator (accuracy 10 μm) by adjusting the posi-
between motor shaft-centre and centres of bearings tion of housing with respect to the motor shaft. On
housings was carried out. First two bearings, as shown obtaining a well-aligned position (dial indicator read-
in Fig. 4, were mounted on the lower part of housing- ings within 10 μm) of shafts, the bearing-housings
slot. Two aluminium sleeves of same size of rotor were bolted to the frame.
magnets were mounted on the shaft at the position of The developed bearing setup, shown in Fig. 3, was
the rotor magnets. These non-magnetic sleeves were used to test hydrodynamic bearings made of alu-
required to avoid magnetic disturbance between shaft minium and neodymium iron boron (NdFeB) mate-
and bearing magnet during alignment. The shaft along rials. Hydrodynamic bearings, having outer diameter
with these sleeves was inserted in bearing. The shaft (OD) of rotating part = 39.94 ± 0.02 mm and inner
centre was aligned with the motor shaft centre using diameter (ID) of stationary part = 40.0 mm, are made
of aluminium material. A photograph of an aluminium
bearing is shown in Fig. 4. Magnetic bearings, hav-
ing OD of rotating part = 39.94 ± 0.02 mm and ID of
stationary part = 40.0 ± 0.02 mm, are made of NdFeB
permanent magnets. This magnetic material has been
selected due to its high energy (280 kJ/m) and high
coercive force (915 kA/m). Aluminium bearings were
tested for 3 h operation at 115 r/min shaft speed by
supplying lubricant of 5 mPa s viscosity at 1.5 bar feed
pressure. Photographs of the worn and unworn alu-
minium bearing are shown in Fig. 5. The marked
area on this figure indicates wear that has occurred
due to lack of the hydrodynamic film. This visual
inspection is sufficient to point out drawbacks of
hydrodynamic bearing at low operating speed. Sub-
sequently, the NdFeB bearings were tested under the
same operating conditions for 3 h. The magnetic bear-
ing was expected to provide non-contact support but
the opposite was observed with full contact occurring.
Fig. 4 Hybrid bearing Figure 6 shows the relatively high wear of the magnetic
Proc. IMechE Vol. 221 Part J: J. Engineering Tribology JET282 © IMechE 2007
Hybrid journal bearings 885
JET282 © IMechE 2007 Proc. IMechE Vol. 221 Part J: J. Engineering Tribology
886 H Hirani and P Samanta
Substitution of equations (2) and (3) in equation (1) along z-axis, so that
leads to
Br2 Rr3 dRr3 Rs2 dRs2 dα dβ
dFrlasx =
Br2 ds1 ds2 R r3s2 4πμ0
dF r3s2 = or,
4πμ0 R3r3s2 2 1 1
× 3
− 3
− 3
Rr3s2 Rr3s1 Rr4s2
Br2 Rr3 dRr3 Rs2 dRs2 dα dβR r3s2
dF r3s2 = (4) × (e + Rr3 cos α + Rs2 cos β) (10)
4πμ0 R3r3s2
Br2 Rr3
dRr3 Rs2 dRs2 dα dβ
dFrlasy =
where Rr3s2 = [x02 + (e + Rr3 cos α − Rs2 cos β)2 + (Rr3 4πμ0
sin α − Rs2 sin β)2 ]1/2 ×
2 1
− 3 − 3
1
3
Similarly, considering no misalignment between Rr3s2 Rr3s1 Rr4s2
lower-arc-stator and rotor, the elemental magnetic × (Rr3 sin α − Rs2 sin β) (11)
force between faces 1 and 4, faces 2 and 4, faces 1 and
3 can be written as Br2 Rr3
dRr3 Rs2 dRs2 dα dβ
dFrlasz =
4πμ0
Br2 Rr3 dRr3 Rs2 dRs2 dα dβR r4s1 2 1 1
dF r4s1 = (5) × 3
− 3 − 3
3
4πμ0 Rr4s1 Rr3s2 Rr3s1 Rr4s2
Br2 Rr3 dRr3 Rs2 dRs2 dα dβR r4s2 2x0 (L + x0 ) (L − x0 )
dF r4s2 = − (6) × 3
+ 3
− 3
(12)
3
4πμ0 Rr4s2 Rr3s2 Rr3s1 Rr4s2
Br2 Rr3 dRr3 Rs2 dRs2 dα dβR r3s1 On integrating equations (10) to (12)
dF r3s1 = − 3
(7)
4πμ0 Rr3s1 π/2 2π R4 R2
Br2 2 1 1
Frlasx = 3
− 3 − 3
4πμ0 −π/2 0 R3 R1 Rr3s2 Rr3s1 Rr4s2
respectively, where
× (e + Rr3 cos α + Rs2 cos β)
Rr4s1 = Rr3s2 = [x02 + (e + Rr3 cos α + Rs2 cos β) 2
× Rr3 dRr3 Rs2 dRs2 dα dβ (13)
+ (Rr3 sin α − Rs2 sin β) ] 2 1/2
π/2 2π R4 R2
Br2 2 1 1
Rr4s2 = [(L − x0 ) + (e + Rr3 cos α + Rs2 cos β)
2 2 Frlasy = 3
− 3
− 3
4πμ0 −π/2 0 R3 R1 Rr3s2 Rr3s1 Rr4s2
+ (Rr3 sin α − Rs2 sin β) ] 2 1/2
× (Rr3 sin α − Rs2 sin β)
Rr3s1 = [(L + x0 )2 + (e + Rr3 cos α + Rs2 cos β)2
× Rr3 dRr3 Rs2 dRs2 dα dβ (14)
+ (Rr3 sin α − Rs2 sin β)2 ]1/2
π/2 2π R4 R2
Br2 2x0 (L + x0 )
To account for attractive forces between faces 1 and 3, Frlasz = 3
+ 3
4πμ0 −π/2 0 R3 R1 Rr3s2 Rr3s1
and faces 2 and 4, a negative sign in equations (6)
and (7) has been introduced. The total elemental force (L − x0 )
− 3
Rr3 dRr3 Rs2 dRs2 dα dβ (15)
between the lower-arc-stator and rotor is obtained by Rr4s2
adding all the elemental forces
Similarly, forces Fruasx , Fruasy , and Fruasz between upper-
arc-stator and rotor magnets can be calculated. The
dF rlas = dF r3s2 + dF r4s1 + dF r4s2 + dF r3s1 (8) total forces (Fx , Fy , and Fz ) in x-, y-, and z-directions,
respectively, will be the summation of the forces, due
to the lower-arc-stator (Frlasx , Frlasy , and Frlasz ) and the
Symmetry of the magnetic rings makes dF r3s2 = dF r4s2 ,
upper-arc-stator (Fruasx , Fruasy , and Fruasz ). The pres-
therefore equation (8) can be rewritten as
ence of the upper arc, reduces the resultant radial
force (force in y-direction) and increases the axial force
dF rlas = 2 · dF r3s2 + dF r4s2 + dF r3s1 (9) (force in z-direction).
The resultant elemental forces between the lower- 4 SIMULATION OF MAGNETIC BEARINGS
stator-arc and rotor can be divided into three compo-
nents, a radial component along y-axis, a transverse The magnitude of the force components in x-, y-,
component along x-axis and an axial component and z-directions, described in the previous section,
Proc. IMechE Vol. 221 Part J: J. Engineering Tribology JET282 © IMechE 2007
Hybrid journal bearings 887
depend on the thickness, air-gap (eccentricity), arc- 180◦ arc-stator, while the minimum axial thrust occurs
extent and length of the magnets. To perform a for the minimum possible extent of stator arc.
parametric study, MATLAB code has been developed. Ohji et al. [11] described the need for an upper arc
To ensure the validity of the developed code; results to enhance the stiffening of magnetic bearing. They
of radial force are compared with published results considered the case with R1 = 20 mm, R2 = 40 mm,
in references [10] and [11]. Qingchang et al. [10] R3 = 45 mm, R4 = 55 mm, L = 7 mm, Br1 = 1.17 T, and
considered full 360◦ rotor as well as stator. The Br2 = 0.38 T. Radial clearance between stator and rotor
input data used in published paper [10] are: inner is equal to 5 mm. Using these input parameters in
radius of rotor (R1 ) = 10 mm, outer radius of rotor the MATLAB code, values of radial force are obtained.
(R2 ) = 19 mm, inner radius of stator (R3 ) = 21 mm, These results are presented in Fig. 9. Again there is
outer radius of stator (R4 ) = 30 mm, length of bearing a close matching between the results obtained in the
(L) = 7 mm, and remanent (residual) induction (Br ) = present study and those given in reference [11].
1.25 T. Radial clearance between stator and rotor is
equal to 2 mm. The results obtained for these inputs
5 HYBRIDIZATION OF MAGNETIC AND
are plotted in Fig. 8, which shows a very close match
HYDRODYNAMIC BEARINGS
with the published results of Qingchang et al. [10].
Figure 8 provides very interesting observations that
A number of problems can be envisaged in hybridiza-
magnetic bearings can be operated with very high
tion of hydrodynamic and magnetic bearings. First
radial clearance (i.e. 2 mm) and eccentricity ratio
and foremost, the problem comes from the clearance
(eccentricity/radial-clearance) follow approximately
required for hydrodynamic action, which approaches
linear relation.
the tolerances required for the permanent magnets.
This exercise prompted the curiosity to explore the
For example, a diametral clearance of 50 μm is
effect of the extent of the stator. Rings with angular
required for a 50 mm bearing, and the minimum toler-
extent of 60◦ , 120◦ , 180◦ , 240◦ , 300◦ , and 360◦ were con-
ance required for the magnet is also 50 μm. Magnetic
sidered in the analysis. The results for radial force and
bearings can be operated with larger radial clear-
axial force are listed in Table 2. It is evident that the
ance, for example Qingchang et al. [10] used 2 mm
maximum value of the radial force is obtained for a
as radial clearance, while Ohji et al. [11] used 5 mm
clearance. Unfortunately, with such large clearance
appropriate hydrodynamic film thickness cannot be
achieved. The second problem comes from the angu-
lar extent of the bearing. The hydrodynamic bearing
requires a full 360◦ arc bearing for minimization of
oil flow, while the magnetic bearing requires angu-
lar extent of approximately 180◦ . The third problem
comes from the attitude angle of the hydrodynamic
bearing, which influences the magnitude of magnetic
force and couples magnetism with hydrodynamics.
Keeping these problems under consideration, the
bearing shown in Fig. 4, has been manufactured.
This design of bearing maintains larger radial clear-
ance between rotor-magnet and stator-magnet and
maintains low clearance between rotor-magnet and
Fig. 8 Radial magnetic force for full ring stator aluminium stator ring.
Table 2 Force components (in N) obtained in the present study for the case with
L = 15 mm, e = 45 μm, CLCE = 50 μm, Br = 1.2 T, R1 = 10 mm, R2 = 16 mm,
R3 = R2 + CLCE, R4 = 22 mm
JET282 © IMechE 2007 Proc. IMechE Vol. 221 Part J: J. Engineering Tribology
888 H Hirani and P Samanta
To analyse the manufactured bearing a code has Table 3 Simulated results of magnetic and hydrody-
been developed, which involves the following steps. namic forces for hybrid bearing
Proc. IMechE Vol. 221 Part J: J. Engineering Tribology JET282 © IMechE 2007
Hybrid journal bearings 889
No 500 74.8 75 73.4 76.4 77.2 75.36 1.03 0.99 1.06 1.01 1.02 1.022 58
Yes 500 75.7 77 76.3 76.4 74.2 75.92 1.03 1.02 1.01 1.01 1.01 1.016 58
No 910 129 122 124 123 123 124.2 0.88 0.87 0.86 0.87 0.87 0.87 81
Yes 910 120 122 122 126 122 122.4 0.86 0.87 0.87 0.86 0.85 0.862 79
No 2000 269 267 267 267 267 267.4 0.85 0.85 0.85 0.85 0.85 0.85 170
Yes 2000 268 268 268 268 268 268 0.85 0.85 0.85 0.85 0.85 0.85 171
No 3000 337 336 336 338 337 336.8 0.74 0.74 0.74 0.74 0.74 0.74 187
Yes 3000 337 337 336 337 337 336.8 0.74 0.73 0.74 0.74 0.74 0.738 186
JET282 © IMechE 2007 Proc. IMechE Vol. 221 Part J: J. Engineering Tribology
890 H Hirani and P Samanta
REFERENCES
Proc. IMechE Vol. 221 Part J: J. Engineering Tribology JET282 © IMechE 2007
Hybrid journal bearings 891
JET282 © IMechE 2007 Proc. IMechE Vol. 221 Part J: J. Engineering Tribology