Case Study 1
Case Study 1
Research Article
Dynamic Modal Analysis of Vertical Machining
Centre Components
The paper presents a systematic procedure and details of the use of experimental and analytical modal analysis technique for
structural dynamic evaluation processes of a vertical machining centre. The main results deal with assessment of the mode shape
of the different components of the vertical machining centre. The simplified experimental modal analysis of different components
of milling machine was carried out. This model of the different machine tool’s structure is made by design software and analyzed by
finite element simulation using ABAQUS software to extract the different theoretical mode shape of the components. The model
is evaluated and corrected with experimental results by modal testing of the machine components in which the natural frequencies
and the shape of vibration modes are analyzed. The analysis resulted in determination of the direction of the maximal compliance
of a particular machine component.
Copyright © 2009 Anayet U. Patwari et al. 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.
structures, leading to lighter, stronger, and safer structure 3. Experimental Modal Analysis
with better performance. In modal analysis, a mathematical
model of a structure’s dynamic behavior is obtained. The 3.1. Measurement Hardware. A vibration measurement gen-
mathematical model consists of a set of mode shapes each erally requires several hardware components. The basic
with an associated natural frequency and modal damping. hardware elements required consist of a source of excitation,
These modal parameters provide a complete description of called an exciter (Impulse hammer), for providing a known
the structure’s dynamic behavior. Baker and his coresearchers or controlled input force to the structure, a transducer
used finite element method to analyze the instability of to convert the mechanical motion of the structure into
machining process [5]. They created a structural model of a electrical signal, a signal conditioning amplifier, and an
analysis system in which modal analysis program resides.
machine tool system using the commercial FE code, ANSYS,
The schematic diagram of hardware used performing
without any experimental tests. In this model the bed of
in a vibration test is shown in Figure 1. The different
machine tool, spindle, and tool holder, as steel blocks, are
equipments that have been used are listed as follows: Pulse
modeled. The integrity of these models is not confirmed by
Front-end (Data Acquisition), Impact Hammer, USB Don-
experimental results. In another research with the aim of
gle, Accelerometers, Impact Hammer cable, Accelerometer
analysis of chatter phenomena, the tool’s natural frequencies
cables, Pulse Front-End Power Supply, TCP/IP Cross Cable,
and the shape of their vibration modes were obtained by
and Bee’s wax.
modal testing results. In this case the variations of acoustic
emission signal during chatter are analyzed, so that it can be
used for chatter detection in machining duration [6]. Many 3.2. Test Procedures. The different milling machine compo-
researchers [7–9] tried to analyze the static and dynamic nents were identified which play a dominating role for the
analysis of the structure involved in machining system by chatter generation. The natural frequency of the different
resting using stiffness measurements and modal analysis. components was measured using modal analysis under static
Talantov and Amin have observed that chatter arising during and dynamic conditions and consequently the different
turning is a result of resonance, caused by mutual interaction mode shapes were identified. Initially excited frequencies
of the vibrations due to serrated elements of the chip were monitored during the operational mode under no-
and the natural vibrations of the system components, for load condition. It is easy to record a response in vibration
example, the spindle and the tool holder [10–12]. The chatter during machining but almost impossible to measure the
phenomena were indicated by the some of the researchers mentioned dynamic force. Therefore, the force measurement
as a resonance effect where system components played a was replaced by measurement of impulse response to the
vital role. So it is important to extract the accurate mode impact force excited by a hammer, whose tip was fitted
shape of the dominating components of machine structure by a force sensor. As the goal of these measurements was
to identify the chatter formation causes. The paper is focused to evaluate frequency transfer function, the responses at
at dynamic properties of a vertical milling machine, namely, various machine points with respect to a reference point were
at the resonance frequencies and vibration shapes of a recorded and analyzed. The reference point was selected at
vertical machining centre components. All this properties are the different location shown in Figure 2.
identified by measurements. The machine tool vibration was
excited by impulse force and a response of excited vibration (i) Knocking test. The natural frequencies of the different
was recorded. The measurement points for vibration were components were extracted from the recorded FFT
selected at the different location of spindle, tool, and diagram. One accelerometer was connected to the
collet. component; the natural frequency data from the
FFT graph was recorded by knocking the different
components using the impact hammer.
2. Simplified Modal Analysis of Milling Machine
(ii) Operational test. The dominating frequencies were
Aiming to investigate the vibration phenomena occur- identified considering high values based on the
ring occasionally at the different components of milling natural frequencies obtain from the knocking test
machine experimental and analytical modal analyses were during no-load operating condition. Accelerometers
performed. The study focused on extracting the mode shape are connected to the components and data were
of the dominating components of the milling machine recorded in Auto spectrum graph. The time exci-
in order to ensure resonance phenomena as a cause of tation, the coherence response excitation, and time
chatter. In a first step the significant eigen-frequencies with response excitation were also recorded to investigate
corresponding mode shapes were obtained by means of the quality of the signal as shown in Figure 3.
an experimental modal analysis (EMA). Subsequently, the
dynamic behavior of the machine components was simulated Experimental Modal Analysis is based on determining
using an ABAQUS FE model. The comparison of the eigen- the modal parameters by testing, unlike Analytical Modal
frequencies based on FE calculations with their experimental Analysis, where the modal parameters are derived from
counterparts proved in general quite satisfactory correla- Finite Element Models (FEMs). There are two ways of doing
tion. Experimental Modal Analysis: Classical Modal Analysis and
Advances in Acoustics and Vibration 3
(mA2 /N)
Accelerometer
2
Impact hammer 1
0
0 4 8 12 16 20
(kHz)
Knocking test
Accelerometer
Impact hammer
Knocking test
Accelerator position
Knocking position
Figure 2: Position of Accelerometer and knocking point at different components of milling machine.
(m)
0.0075 0.8
0.005
0.4
0.0025
0 0
0 2500 5000 7500 10000 12500 0 200 400 600 800 1000 1200 1400 1600
(Hz) (Hz)
Y/t chart0
(a) (b)
400
1
(m/s2 )
0 0
(N)
−1
−400
−2
0 10 20 30 40 50 60 70 0 10 20 30 40 50 60 70
(ms) (ms)
(c) (d)
0.0175 0.02
0.015 0.0175
0.015
0.0125
0.0125
0.01
0.01
0.0075
0.0075
0.005
0.005
0.0025 0.0025
0 0
0 2500 5000 7500 10000 12500 0 25 50 75 100 125 150 175 200 225
×102
(Hz) (Hz)
Figure 4: Experimental modal analysis power spectrum responses (FFT); (a) Collet and (b) Spindle casing of vertical machining centre.
Advances in Acoustics and Vibration 5
3
2
(m/s2 )
(m/s2 )
2
1
1
0 0
0 4 8 12 16 20 0 4 8 12 16 20
(kHz) (kHz)
1000 RPM 3500 RPM
(a) Inner spindle operational deflection shape with no load. Spindle mode shape: 4394 Hz, 6866 Hz, 8294 Hz
(m/s2 )
8
1
4
0 0
0 4 8 12 16 20 0 4 8 12 16 20
(kHz) (kHz)
1000 RPM 3500 RPM
(b) Inner spindle operational deflection shape with load. Spindle mode shape: 4394 Hz, 6866 Hz, 8294 Hz
Figure 5: Operational Modal analysis of inner spindle with no-load and load conditions.
Tool holder
of an FE model for dynamic simulation purposes and Table 1: Natural frequencies analysis of spindle.
to detect any possible improvements. The modal-analyzed
Different major mode shape (Hz)
natural frequencies of these models via ABAQUS software
of different components at different mode shapes are shown Spindle Inner Spindle Outer spindle
in Figures 5–7. The results of analysis of the calculated Condition 1st 2nd 3rd 1st 2nd
and experimental mode shape are mentioned in the result Experimental 4394 6866 8294 9009 11371
analysis section. Theoretical 4322 6389 8004 9520 10291
Error (%) 1.67 7.47 3.62 5.37 10.49
Model Number 1: Complete Model of Spindle Both Outer and
Inner of Eigen-Frequency Analysis. Finite element analysis of Table 2: Natural frequencies analysis of Collet.
the inner and outer spindle is performed using ABAQUS.
The calculated distortions of the each element are shown in Collet Different major mode shape (Hz)
Figure 7 at different mode shapes. The elements used in the Condition 1st 2nd 3rd 4th
FEM model for mesh generation are quadratic tetrahedral Experimental 659 2032 3955 5169.7
element. It has been observed from the calculated results that Theoretical 770.77 2318.9 3673.6 5318.2
there are five prominent mode shape of the spindle. In most Error(%) 14.50 12.37 7.66 2.79
of the mode shape the distortion in the inner spindle is more
significant than the outer spindle.
Table 3: Natural frequencies analysis of Tool Holder.
Model Number 2: Complete Model of Collet for Eigen- Tool Holder Different major mode shape (Hz)
Frequency Analysis. Finite element analysis of the collet is Condition 1st 2nd 3rd
again performed using ABAQUS. The calculated distortions Experimental 2081 8892 11352
of the each element are shown in Figure 8 at different mode Theoretical 2944 10099 10454
shapes.
Error (%) 29.31 11.95 8.59
U, Magnitude U, Magnitude
0.000e + 00 1.160e + 00
0.000e + 00 1.063e + 00
0.000e + 00 9.666e − 01
0.000e + 00 8.700e − 01
0.000e + 00 7.733e − 01
0.000e + 00 6.766e − 01
0.000e + 00 5.800e − 01
0.000e + 00 4.833e − 01
0.000e + 00 3.867e − 01
0.000e + 00 2.900e − 01
0.000e + 00 1.933e − 01
0.000e + 00 9.666e − 02
0.000e + 00 0.000e + 00
ODB: spindle.odb ABAQUS/standard version 6.6 × 1 ODB: spindle.odb ABAQUS/standard version 6.6 × 1
Wed Jun 04 12:12:56 malay peninsula Wed Jun 04 12:12:56 malay peninsula standard
Step: Frequency, frequency analysis Step: Frequency, frequency analysis
Increment 0: Base state Mode 4: Value = 7.37518E + 08 Freq = 4322.2 cycles/time
Primary var: U, magnitude Primary var: U, magnitude
Deformed var: U deformation scale factor: 1.000e + 00 Deformed var: U deformation scale factor: 3.102e + 01
(a) Base state (b) Mode: 4322 Hz
U, Magnitude U, Magnitude
1.000e + 00 1.000e + 00
9.167e − 01 9.167e − 01
8.334e − 01 8.333e − 01
7.500e − 01 7.500e − 01
6.667e − 01 6.667e − 01
5.834e − 01 5.833e − 01
5.000e − 01 5.000e − 01
4.167e − 01 4.167e − 01
3.334e − 01 3.333e − 01
2.500e − 01 2.500e − 01
1.667e − 01 1.667e − 01
8.334e − 02 8.333e − 02
0.000e + 00 0.000e + 00
ODB: spindle.odb ABAQUS/standard version 6.6 × 1 ODB: spindle.odb ABAQUS/standard version 6.6 × 1
Wed Jun 04 12:12:56 malay peninsula standard Wed Jun 04 12:12:56 malay peninsula standard
Step: Frequency, frequency analysis Step: Frequency, frequency analysis
Mode 6: Value = 1.61157E + 09 Freq = 6389.2 cycles/time Mode 7: Value = 2.52912E + 09 Freq = 8004 cycles/time
Primary var: U, magnitude Primary var: U, magnitude
Deformed var: U deformation scale factor: 3.102e + 01 Deformed var: U deformation scale factor: 3.102e + 01
(c) Mode: 6389 Hz (d) Mode: 8004 Hz
U, Magnitude U, Magnitude
1.000e + 00 1.418e + 00
9.167e − 01 1.300e + 00
8.333e − 01 1.182e + 00
7.500e − 01 1.064e + 00
6.667e − 01 9.454e − 01
5.833e − 01 8.272e − 01
5.000e − 01 7.091e − 01
4.167e − 01 5.909e − 01
3.333e − 01 4.727e − 01
2.500e − 01 3.545e − 01
1.667e − 01 2.364e − 01
8.333e − 02 1.182e − 01
0.000e + 00 0.000e + 00
ODB: spindle.odb ABAQUS/standard version 6.6 × 1 ODB: spindle.odb ABAQUS/standard version 6.6 × 1
Wed Jun 04 12:12:56 malay peninsula standard Thu Jun 05 12:19:22 malay peninsula standard
Step: Frequency, frequency analysis Step: Frequency, frequency analysis
Mode 10: Value = 3.57864E + 09 Freq = 9520.9 cycles/time Mode 12: Value = 4.18066E + 09 Freq = 10291 cycles/time
Primary var: U, magnitude Primary var: U, magnitude
Deformed var: U deformation scale factor: 3.102e + 01 Deformed var: U deformation scale factor: 3.102e + 01
(e) Mode: 9520 Hz (f) Mode: 10291 Hz
U, Magnitude U, Magnitude
0.000e + 00 1.005e + 00
0.000e + 00 9.215e − 01
0.000e + 00 8.378e − 01
0.000e + 00 7.540e − 01
0.000e + 00 6.702e − 01
0.000e + 00 5.864e − 01
0.000e + 00 5.027e − 01
0.000e + 00 4.189e − 01
0.000e + 00 3.351e − 01
0.000e + 00 2.513e − 01
0.000e + 00 1.676e − 01
0.000e + 00 8.378e − 02
0.000e + 00 0.000e + 00
ODB: assembly spindle.odb ABAQUS/standard version 6.6 × 1 ODB: collet.odb ABAQUS/standard version 6.6 × 1
Mon Jun 09 18:11:02 malay peninsula standard Thu May 22 15:18:55 malay peninsula standard
Step: Frequency, frequency analysis Step: Frequency analysis, frequency analysis
Increment 0: Base state Mode 1: Value = 2.34537E + 07 Freq = 770.77 cycles/time
Primary var: U, magnitude Primary var: U, magnitude
Deformed var: U deformation scale factor: 1.000e + 00 Deformed var: U deformation scale factor: 2.300e + 01
(a) Base state (b) Mode: 770.77 Hz
U, Magnitude U, Magnitude
1.000e + 00 1.006e + 00
9.167e − 01 9.224e − 01
8.333e − 01 8.385e − 01
7.500e − 01 7.547e − 01
6.667e − 01 6.708e − 01
5.833e − 01 5.870e − 01
5.000e − 01 5.031e − 01
4.167e − 01 4.193e − 01
3.333e − 01 3.354e − 01
2.500e − 01 2.516e − 01
1.667e − 01 1.677e − 01
8.333e − 02 8.385e − 02
0.000e + 00 0.000e + 00
ODB: collet.odb ABAQUS/standard version 6.6 × 1 ODB: collet.odb ABAQUS/standard version 6.6 × 1
Thu May 22 15:18:55 malay peninsula standard Thu May 22 15:18:55 malay peninsula standard
Step: Frequency analysis, frequency analysis Step: Frequency analysis, frequency analysis
Mode 3: Value = 2.12287E + 08 Freq = 2318.9 cycles/time Mode 4: Value = 5.32789E + 08 Freq = 3673.6 cycles/time
Primary var: U, magnitude Primary var: U, magnitude
Deformed var: U deformation scale factor: 2.300e + 01 Deformed var: U deformation scale factor: 2.300e + 01
(c) Mode: 2318.9 Hz (d) Mode: 3673.6 Hz
U, Magnitude
1.000e + 00
9.167e − 01
8.333e − 01
7.500e − 01
6.667e − 01
5.833e − 01
5.000e − 01
4.167e − 01
3.333e − 01
2.500e − 01
1.667e − 01
8.333e − 02
0.000e + 00
ODB: tool holder.odb ABAQUS/standard version 6.6 × 1 ODB: tool holder.odb ABAQUS/standard version 6.6 × 1
Wed Mar 05 13:26:05 malay peninsula Wed Mar 05 13:26:05 malay peninsula
ODB: tool holder.odb ABAQUS/standard version 6.6 × 1 ODB: tool holder.odb ABAQUS/standard version 6.6 × 1
Wed Mar 05 13:26:05 malay peninsula Wed Mar 05 13:26:05 malay peninsula
components are within the accepted ranges and the high will be up to 12000 Hz. This research work will help to
error in some of them might be referred to the boundary find out the natural frequencies of the components and
conditions specification, because it is not easy to simulate the hence predicting the chatter formation zone as resonance
realistic boundary conditions for such complicated system. phenomena.
9. Conclusion References
In this paper a finite element model is used to analyze [1] D. J. Inman, Engineering Vibration, Prentice-Hall, Englewood
the mode frequencies and shapes of different machining Cliffs, NJ, USA, 3rd edition, 2007.
components and hence compare the results with the exper- [2] T. C. Huang, Modal Analysis Modeling, Diagnostics and
imental one. This model is produced in CATIA software Control-Analytical and Experimental, ASME, New York, NY,
based on the real dimensions of vertical machining centre USA, 1991.
(model: MCFV 1060LR) machine and analysis was done [3] C. W. de Silva, Vibration and Shock Hand Book, Taylor &
by ABAQUS software. According to the model analysis, the Francis, Boca Raton, Fla, USA, 2005.
natural frequencies and vibration modes shape of the model [4] N. M. M. Maia and J. M. M. E. Silva, Theoretical and
in spindle, collet, and tool holder cases were determined Experimental Modal Analysis, John Wiley & Sons, New York,
and evaluated. The comparison between natural frequencies NY, USA, 1997.
of finite element modeling and model testing shows the [5] J. R. Baker and K. E. Rouch, “Use of finite element structural
closeness of the results. From the results, it has been models in analyzing machine tool chatter,” Finite Elements in
observed that the suitable frequency ranges for end milling Analysis and Design, vol. 38, no. 11, pp. 1029–1046, 2002.
10 Advances in Acoustics and Vibration
Rotating
Machinery
The Scientific
Engineering
Journal of Advances in
Journal of Mechanical
World Journal
Hindawi Publishing Corporation Hindawi Publishing Corporation Hindawi Publishing Corporation
Sensors
Hindawi Publishing Corporation
Engineering
Hindawi Publishing Corporation
http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014
International Journal of
Distributed Advances in
Sensor Networks
Hindawi Publishing Corporation
Civil Engineering
Hindawi Publishing Corporation
http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014
Journal of
Robotics
Hindawi Publishing Corporation
Hindawi Publishing Corporation Volume 2014 http://www.hindawi.com Volume 2014
http://www.hindawi.com
VLSI Design
Modelling &
International Journal of Simulation
International Journal of Navigation and in Engineering Advances in
Chemical Engineering Observation Acoustics and Vibration
Hindawi Publishing Corporation Hindawi Publishing Corporation Hindawi Publishing Corporation Hindawi Publishing Corporation
http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014
Hindawi Publishing Corporation
http://www.hindawi.com Volume 2014
Journal of
Control Science
and Engineering