Modeling and Simulation of High Speed Spindle, Current Problems and Optimizations
Modeling and Simulation of High Speed Spindle, Current Problems and Optimizations
Modeling and Simulation of High Speed Spindle, Current Problems and Optimizations
Andrei Gheorghita 1, Mariana Rotariu 2, Marius Turnea 3, Mihai Ilea 4, Dragos Arotaritei 5
1
University Politehnica of Bucharest,
Splaiul Independentei 313, Bucharest, Romania
andu_ghe@yahoo.com, +40720272072
2
University of Medicine and Pharmacy “Grigore T. Popa”,
Str. Universitatii 16, Iasi, Romania
rotariu29@yahoo.com, +40723133560
3
University of Medicine and Pharmacy “Grigore T. Popa”,
Str. Universitatii 16, Iasi, Romania
marius.turnea@gmail.com, +40766416165
4
University of Medicine and Pharmacy “Grigore T. Popa”,
Str. Universitatii 16, Iasi, Romania
ileamihai2004@yahoo.com, +40744869328
5
University of Medicine and Pharmacy “Grigore T. Popa”,
Str. Universitatii 16, Iasi, Romania
dragos.arotaritei@gmail.com, +40742181273
1. Introduction:
Rotating machineries like lathes, grinding machines, milling machines are very important
and useful for modern manufacturing processes. The assessment of a good functionality and a
performance issue involves the deformation, vibrations and temperature due basically to unbalance of
the main spindle, friction effects in bearings, cutting forces, chatter, and cooling systems. The effects
can be a poor quality of the manufacturing product but also the wear or even destruction of machine
as effect of a long time stress reflected by unfit functionality.
Spindle bearing systems belong to rotating machinery class and the operations of cutting
material involves an analysis of dynamical model of the spindle bearing system. Many of the research
articles use the approaches that are proposed for rotor dynamic analysis but in the spindle shaft case,
the approach is more complex.
The dynamic behavior of spindle-bearing system is essential to be understood in order to
design a performant machine that uses this system. The design can be less or more complex and the
models used in simulation can vary from case to case depending on the tradeoff between accuracy of
result and complexity of model and calculus.
There are many types of bearings used for rotors in rotating machinery and one of the most
frequently type is angular contact ball bearings. The angular contact ball bearing offers some
advantages in rotating machinery worth to be mentioned: good stiffness properties, low friction, long
lifetime, and a good ration between cost and effectiveness especially at high speed. The critical
subsystem in the rotating machinery is made by ball bearings that transfer the contact force, heat due
to friction and vibration that goes thru the other elements to the cage of machine. The modeling of
bearing dynamic influences the results of dynamic analysis for shaft bearing system. The spindle
bearing system has at least two sets of bearings mounted in different positions depending of design
condition and optimization objectives.
There are many studies that analyze the theoretically and experimentally the dynamic
behavior of machines that use spindle shafts, lining, rolling bearings, and spindle-bearing systems.
There are also few review papers that deal with this subject. The purpose of this paper is to present
structured view of the main approaches in dynamic modeling of spindle-bearing systems along with
main possible optimizations.
3
Q
( )
L= c
Q (1)
0. 41 1 .39 0. 3
R r (1∓γ ) γ
Qc =99 .1 ( D r−R ) (1±γ ) 1/3 cosα ( ) D1 .8 Z−1/3
(2)
0.41 1 .39 0.3
2f (1∓γ ) γ
Qc =99.1 ( ) 2 f −1 (1±γ ) 1/3 cosα ( ) D1 .8 Z−1/3
(3)
The sign above is for the inner ring and the one below is for the outer ring. Taking into
account both the inner ring i and outer ring o, their lifetime is:
−1/e
L=( L−e −e
i +Lo ) (4)
For the contact point e=10/9 (ball bearing, angular contact), while for the roller e=9/8. For
rolls, with a contact line l, the formulas become:
4
Q
L= c
Q ( ) (5)
29/27 2/9
(1∓γ ) γ
Qc =552
(1±γ )1/4 cosα ( ) D29/27 l7/9 Z−1/4
(6)
The sign Q is a sum for all the other Z balls, while each force Q j =Q(ψ ) , 2/(Z-1):
Z 1/3
1
(
Q= ∑ Q3j
Z j=1 ) (7)
It can be noticed that the durability of ball bearings is a function that depends on the
geometrical characteristics and the material of the bearing, on the pre-gripping force and the
rotational speed of the shaft which, in turn, determines the angle of contact [5].
The angle of contact of the ball bearing towards its inner and outer running track is
constant in the absence of a centrifugal force. However, the value of this angle is greater if you apply
a pre-gripping force. Under conditions of free pre-gripping:
Pd d o−d i−2 D
cos α 0 =1− =1−
2 BD 2(f o + f i−1 )D (8)
Preloading for a contact angle is calculated iteratively, method used to assess the durability
of a bearing 7012AC, where i is the step:
Fa 1. 5
cos α 0
α (i+1 )=α (i+1 )+
ZD 2 K
−sin α
cos α ( −1 )
1. 5 0. 5
cosα 0 cos α 0
cos α (
cosα
−1 ) +1 .5 tan(cos α
−1 ) cos α 0
(9)
In the formula (9) the displacement constant is a function of curvature and is taken from
graphics or tables. Flow-through, however, should be considered as values of angles for each ball
(Fig. 1- 2).
Fig. 1 Positions of the balls in the radial plane yz (x is the bearing axis) and the forces
acting upon the bearing [1].
Fig. 2 Shifting centers of curvature. The position of the ball bearings [1]
Finally, calculating the maximum angles, which are valid for ψ=0 , for a ball bearing
type 218, one obtains:
Fig. 3 Angles
αi and αo (for =0)
The balls are made of steel, but in the case of ceramic balls, both the ball mass, material
constant and the centrifugal force and moments. As a result, it is expected that the internal and
external angles of contact with the inner and outer ring should be smaller.
LT’s lifetime, calculated approximately, once the iterations don’t improve the solution, is:
Fig. 4 LT function of pre-gripping force and the speed of rotation of the shaft, flow-through.
It is observed that minimum solutions may be obtained for pairs of values, whether f1 is
lower and f2 greater or vice versa, a so-called Pareto front convex, which through the selection of a
point, in the compliance with an objective, be it for example f1, f2, is obtained for a minimum of the
cost function.
4. Conclusions
The calculation of durability of the entire bearing involves some laborious calculations and
some rather large matrices. On the other hand, the stiff matrices of the bearing should be included in
the shaft-bearings matrix assembly, which could lead to difficult numerical resolutions of the radial
displacements in the nodes that shape the tree with finite elements.
The equations which result from the angles i and o present a lot of circular dependencies;
as a result, the convergence of solutions in the case of numerical methods depends on the initial
solution, which in this case must be pretty close to the final solution. Since the variations in relation
to the static case can be very large, the viable solution is to fragment the space of combinations of the
two angles for an i, in the range 0-2. If this fragmentation is very fine (the domain of values for
which we are trying different initial solutions), it can be explosive and as such infeasible.
Bibliografie
[1] Tedric A Harris, Rolling Bearing Analysis, 4th Edition, Wiley-Interscience, 2001
[2] Genta, G., 2005, Dynamics of Rotating Systems, Springer, NY. [10] Rao J.S., Rotor
Dynamics, Third ed., New Age, New Delhi, 1996. [11] J.S. Rao, History of Rotating
Machinery Dynamics, 2011
[3] M.I. Friswell, J.E.T. Penny, S.D. Garvey, A.W. Lees, Dynamics of Rotating Machines,
Cambridge University Press, 2010.
[4] Dana Codruța Vișan, Modelarea Dinamică a Sistemelor Rotor-Lagăre, Editura PRINTECH,
2008.
[5] D. L. Logan, A First Course in the Finite Element Method, 5th Edition, CL Engineering,
2011
[6] L. Meirovitch, Analytical Methods in Vibration. The Macmillan Company, New York, 1967.
[7] L. Meirovitch, Elements of Vibration Analysis. McGraw-Hill, Inc., New York, 1986.
[8] A. Gheorghita, D. Arotaritei, M. Turnea, G. Constantin, Modelling and simulation of high
seed spindle, current problems and optimizations – a survey ,ICMS, 2016
[9] H. D. Nelson, A Finite Rotating Shaft Element Using Timoshenko Beam Theory, ASME J.
Mech. Des., 102, 1980, pp. 793-803.
[10] Y. Cao, Y. Altintas, A General Method for the Modeling of Spindle-Bearing Systems,
Journal of Mechanical Design, Vol. 126, 2004, pp. 1089-1108
[11] Y. Cao, Modeling of High-Speed Machine-Tool Spindle Systems, PhD thesis, The
University of British Columbia, 2006.
[12] K.E. Rouch, J.S. Kao, A tapered beam finite element for rotor dynamics analysis, J. Sound
Vibration, Vol. 66, 1979, pp. 119-140.
[13] M. Turnea, D. Arotaritei, R. Filep, A. Gheorghita, M. Rotariu, C. Budacu, Interface analysis of
temperature
inside the prosthesis with effect in improving the quality of life of 98 patients, ISMH 2016.
[14] H. Cao, X. Zhang, X. Chen, The concept and progress of intelligent spindles: A review,
International Journal of Machine Tools & Manufacture, Vol. 112, 2017, pp. 21-52.
[15] D. Arotaritei, M. Turnea, R. Filep, M. Ilea, M. Rotariu, A. Gheorghita, Graphical Interface for
Analyze of Temperature of Skin at Skin-Liner Interface, The 5th IEEE International Conference on E-
Health and Bioengineering - EHB 2015