Research Article  Open Access
Vibration Characteristics of HeavyDuty CNC Machine ToolFoundation Systems
Abstract
Vibration characteristics of heavy CNC machine tools are directly affected by their foundations. To analyze vibrations of heavy CNC machine tools caused by internal and external loads, a system dynamics model of a rigidflexible coupled heavyduty CNC machine toolfoundation system was established based on the multibody transfer matrix method. Since joint surfaces can seriously affect the accuracy of system mechanics models, the substructure synthesis method was first used to establish a dynamic model of the joint surface. The frequency response function was then used to identify model parameters. Moreover, to improve the accuracy of parameter identification of the joint surface, a residual frequency compensation function was used to reconstruct the frequency response function. Finally, the multibody system model was implemented by combining surface elements. To verify the system dynamics model, an experimental model of the heavyduty machinefoundation system was built, taking into consideration joint surface factors, and the theoretical model was validated by comparing theoretical, simulation, and experimental results. Using the theoretical model, the influence of different forms of concrete foundations, materials, and soil properties on the vibration characteristics of heavyduty CNC machine tools was analyzed, thus providing a theoretical basis for optimizing and improving CNC machine tools.
1. Introduction
Vibrations generated by a machine’s own excitations can seriously interfere with working conditions, particularly those of sensitive equipment. More importantly, relative vibration displacement can occur between a cutting tool and a machined part, which will directly affect the accuracy and surface quality of the machined part. Owing to interactions between the heavyduty CNC machine tool and its foundation, there are significant differences in the calculation results obtained by directly using acceleration time histories as the inertial input and considering foundation. Analyzing the dynamic response of the machine tool structure alone is insufficient, and the influence of the concrete foundation on its dynamic characteristics should be considered.
To study the characteristics of the dynamic mechanicalfoundation system, it is necessary to establish an accurate system model. Previous research has been carried out on mechanical models of dynamic mechanicalfoundation systems. Ghosh [1] used the threedimensional finitedifference method to establish a dynamic mechanicalfoundation system model and studied the influence of a range of different artificial boundary conditions on the sensitivity of the static and dynamic characteristics of the system. When the artificial boundary value exceeds ten times the machine tool overall size, the sensitivity of the system significantly decreases. A threedimensional finiteelement model was proposed by Liu et al. [2] that considered a soilpower machine interaction system. Based on this model, the main factors influencing the dynamic response of the largescale machine toolfoundation were revealed. Aşık et al. [3] established an analytical solution for a bar and circular dynamic mechanical foundation model based on the variational principle and Hamiltonian principle of energy minimization. Cai et al. [4] conducted a series of analyses on the characteristics of foundationfoundation systems using cloud computing. Furthermore, Stanescu and Tabacu [5] considered the nonlinearity of soil and established a twodegreeoffreedom (DOF) machinefoundation model to study the stability of the equilibrium position.
The purpose of analyzing the characteristics of power machineryfoundation systems is to study their reliability, stability, and accuracy. However, system characteristics are often influenced by multifactor coupling. To reveal the influence of various factors, research has been conducted on the main factors affecting dynamic mechanical properties. Vivek and Ghosh [6] used the finiteelement method to study the dynamic response of adjacent interacting mechanical foundations and analyzed the influence of distance between foundation profiles on the dynamic characteristics of the system. In another study, Štimac et al. [7] established a finiteelement model of a turbogeneratorfoundation system with the objectives of avoiding resonance caused by higher harmonics of the turbine and generator and optimizing the cross section of the foundation. Haldar and Sivakumar Babu [8] used a method that arranges spatial net stiffeners in the soil foundation to increase the natural frequency of the soilfoundation system, thereby avoiding resonance due to vibrations generated by the machine’s rotations. Werner [9] deduced a vibration analysis model of the fundamental frequency of a machinefoundation system and demonstrated that the fundamental frequency of the system can be increased by increasing the contact area between the machine and foundation.
In each of the studies mentioned above, a dynamic mechanicalfoundation system model based on theoretical and experimental analyses of the dynamic characteristics of the entire superstructure was established for various types of power machinery. However, the influence of joint surfaces was not considered in the system modeling process, and these studies lack sufficient validation of the system models.
In this paper, a mechanical model of a heavyduty CNC machine toolfoundation system that considers the joint surfaces is established based on the multibody transfer matrix method. The theoretical model and simulations were validated by experiments performed in the laboratory on an equivalent scaled CNC machine toolfoundation system. Finally, the effects of different forms, materials, and reinforcement conditions of the concrete foundation, as well as different soil properties, on the dynamic characteristics of heavyduty CNC machine tools were analyzed.
2. Multibody System Theory Model
The advantage of the transfer matrix method is that the stiffness matrix depends on the DOF of each substructure. Compared to the finiteelement method, the total stiffness matrix can be effectively reduced. A multibody mechanical model was established according to the structural characteristics of the heavyduty CNC machine tool and its foundations by considering the upper and lower structural subsystems. A dynamic equation can be established to describe the entire heavyduty CNC machinefoundation system according to the state vector of the subsystem connection interface. A multibody dynamic model of the heavyduty CNC machine toolfoundation system was established, as shown in Figure 1, in which the pillar, beam, connecting beam, bed body, and foundation were considered elastic elements due to their largescale characteristics and the more rigid structures were considered rigid elements.
2.1. Joint Surface Element
Joint surfaces are responsible for 60% of the vibration problems of machine tools. For this reason, a transfer matrix at the joint surface needs to be established. State vectors of the input and output can be defined using the modal coordinate system as and [10].
According to the force balance equation, the joint surface transfer model can be obtained as follows:where is the transfer matrix of the bonding surface element. The transfer matrix of the transfer surface is a function of the state variable and excitation frequency, represented as follows:where is a diagonal matrix incorporating both stiffness and damping. The diagonal matrix is given bywhere are the equivalent stiffness in x, y, and zdirections; are the equivalent damping in x, y, and zdirections; are the equivalent stiffness about the rotational displacement around the x, y, and zdirections; and are the equivalent damping about the rotational displacement around the x, y, and zdirections.
2.2. Rigid Body Elements
During the modeling process, structural members of the heavyduty CNC machine toolfoundation system with high stiffness values are regarded as rigid components. This greatly reduces the DOF of the substructure and also reduces the dimensionality of the system’s stiffness matrix.
The rigid body space vibration transfer model can be expressed as follows:where is the rigid body vibration transfer matrix. The vibration transfer matrix is given bywhere is the vector coordinate matrix of the output point position with respect to the input point position, is the matrix of the center of mass coordinates relative to the input point position vector, and is the inertia matrix of the rigid body with respect to the input point.
2.3. Elastic Elements
An elastic body can be divided into sections, and according to the unit dynamic equation based on the principle governing system dynamics, modal coordinates are introduced. The unit input and unit output states can be substituted into the equation of motion to obtain the transfer matrix of the th unit. Combining units, the transfer model of the entire structure [11] can be defined as follows:where is the overall transfer matrix, , and is the transfer matrix of the th unit, given bywhere is the mass matrix of the th element and is the stiffness matrix of the th element.
2.4. Transfer Matrix of the HeavyDuty CNC Machine ToolFoundation System
The system, China dynamics model of the heavy machine tool, is presented in Figure 2.
State vectors , , , , , , , , , , , , , , , , , , , , , and can be defined. Elements 7, 11, and 13 are joint surface units which are input at both ends and output at one end. The state vectors can be defined by the following relations:
The transfer matrices for each component are
According to Equation (9), the following relations can be obtained:
According to Equation (10),
The relationships at the joint surfaces of elements 7, 11, and 13 can be given by
Equations (11) and (12) can be written in the matrix form to obtain the total transfer matrix of the system, given bywhere is the overall state vector of the superstructure and . Thus, the overall transfer matrix of the superstructure is
Due to their largescale characteristics, concrete foundations of heavyduty gantry CNC machine tools must be considered as elastic bodies (Figure 3). Considering the concrete foundation as an elastic component, the position of the machine tool relative to the concrete foundation can be determined and dynamic characteristics of the machine tool at different positions on the foundation can be more accurately analyzed. Therefore, the finiteelement transfer matrix of the concrete foundation can be expressed as follows:
The upper and lower structural boundary conditions are generally expressed in the following form:
Therefore, the overall transfer matrix of the upper and lower structure is
Then, Equation (17) can be simplified as follows:
3. Parameter Identification of the Joint Surface
The contact stiffness of a joint surface often contributes to the overall stiffness of a mechanical structure and, occasionally, may be the weakest part of the structure contributing to the overall stiffness. Therefore, to study the dynamic characteristics of mechanical structures, the contact stiffness of joint surfaces must be fully considered [12]. In this paper, metaltometal joints and the metaltoconcrete joints are analyzed. Due to the relatively low elasticity of concrete, as pressure acts across the large area of a concrete foundation, most of the protrusions on the concrete surface will break once a certain elastic limit is exceeded. This leads to uncertainty in the real contact area between the concrete and metal structure. Thus, an experimental identification method can be used to establish a unified parameter identification model for the metaltometal and concretetoconcrete interfaces.
3.1. Dynamic Modeling of the Joint Surface
Heavy machine tools can be considered large fixedweight structures owing to their size, and the load acts vertically on joint surfaces. As such, the influence of torsional freedom at the joint surface on the overall structure can be ignored [13]. An overall model of the substructure and interface can be established, as illustrated in Figure 4, and consists of three parts: substructure (1), substructure (2), and bolted joints.
Ignoring the mass of the bolted joints, the relationship between the frequency response function of the substructure and the displacement and external force can be expressed as follows:
The force balance condition at the joint surface of the substructure is
As defined by the frequency response function, the displacement generated by internal forces at the bolted joint surface iswhere .
According to Equations (19) and (22),assuming .
Based on Equations (19), (20), and (23), it follows that
According to Equation (24),where and .
Then, Equation (25) can be modified to obtain
Substituting Equation (24) into Equation (26), we obtain
3.2. Parameter Identification
3.2.1. Identification Theory
From Equation (27), known as the identification equation, it can be seen that, by obtaining the frequency response function of the substructure and overall structure through testing, the equivalent dynamic parameters of the joint surface can be determined, thereby avoiding the need to measure the frequency response function of the integral structural joint. Thus, Equation (27) can be modified such that
Expressing Equation (29) as a linear system of equations, we obtain
Then, modifying Equation (30), we obtain
Finally, using the least squares method to solve Equation (31), we obtain
3.2.2. Weighting
Since the value of the frequency response function varies significantly with changes in frequency, especially at the resonance frequency, identified joint surface parameters are almost completely controlled by resonance frequencies. In other words, when using the least squares method to solve the overdetermined equation, larger coefficients will have a greater influence on the solution. To make full use of the coefficient values, the effects of the coefficients can be averaged using a weight function method.
Simultaneously multiplying both sides of Equation (31) by the weighting matrix , we obtainwhere , in which is the weighting factor, , where is the parameter of matrix .
3.3. Method for Obtaining Frequency Response Function
Based on the identification method described in this paper, the frequency response function of each substructure, as well as the overall structure, must be acquired through testing and substituted into Equation (27). Then solving Equation (32), parameter values at the joint surfaces can be obtained. To identify the joint surface parameter, an experimental device was built. In this paper, the metalconcrete joint surface system is presented as an example. The test system is composed of the steel structural specimen, concrete specimen, strain bolt, force hammer, acceleration sensor, dynamic strain gauge, and signal acquisition and analysis system. Strain bolts were embedded into the concrete specimens at equidistant positions. The pretightening forces of bolts were obtained as strain values during the initial calibration. At the same time, the steel test piece was processed with a corresponding bolt hole, and a 1 mm high boss was machined around the bolt hole on the contact surface between the steel test piece and concrete to ensure full contact between the steel and concrete. The metalconcrete test system is shown in Figure 5.
(a)
(b)
3.4. Frequency Response Function Reconstruction Method
Due to modal truncation and clutter interference during the test process, some modal data are lost and the original experimental data will inevitably contain errors. To reduce these errors, a test frequency response function is reconstructed and the original experimental data are replaced with reconstructed data to improve the accuracy of parameter identification at the joint surface [14, 15].
The frequency response function of a typical mechanical joint iswhere denotes the th natural frequency and is the th modal constant. In Equation (34), the middle term, , is based on data obtained from testing, whereas the loworder term, , and highorder term, , are truncated by the test conditions and are, therefore, residual items that must be compensated.
A series of points on the frequency response function curve obtained during the experiment can be selected, denoted as , where is the frequency. The interval is the actual desired frequency range. Substituting these data points into Equation (34) yieldswhere can be visually observed and extracted from the frequency response function curve; therefore, Equation (35) is an nary linear system of equations. The value of can be obtained by solving the system of equations, and this value can be substituted into Equation (34) to obtain the reconstructed frequency response function. The reconstructed frequency response function covers the frequency range, which is richer than the original frequency response function. This eliminates random signal interference to give a smoother frequency response curve. The reconstructed frequency response curve more clearly illustrates the response characteristics of major orders and eliminates clutter interference, which is useful for subsequent calculations.
4. Dynamic Parameter Test
4.1. Test Device
Field tests cannot be used to verify the effect of differentsized foundation, foundation forms, and soils on the dynamic characteristics of the system. Moreover, the characteristics of heavyduty CNC machine tools make it impossible to carry out dynamic testing. Therefore, it is necessary to establish laboratory models of heavyduty gantry CNC machinefoundations. The laboratorybased test device can solve a number of issues: wasted capital caused by the incompatibility of heavyduty gantry CNC machinefoundation systems is reduced; the influence of joint surfaces between the structural members of the heavyduty gantry CNC machine tool, structural members, and concrete, and the concrete and soil are truly reflected; and studies relating to the influence of environmental factors on machining precision can be successfully carried out. In this paper, a test device that considers the influence of joint surfaces on a heavyduty gantry CNC machine toolfoundation system is established.
4.1.1. Structure of Gantry, Concrete, and Foundation
The gantry system is composed of connecting beams, columns, supporting beams, and beds. Columns are used to machine four sets of holes connected to the crossbeams in the zdirection. Positioning the crossbeam through each set of holes achieves positioning relative to the four column positions. Thus, dynamic analysis of the crossbeam at four different positions can be carried out. The foundation is made of C40 concrete with an overall size of 1 m × 0.5 m. The test device is versatile and can be used to test different machine tools, foundations, and soils. In this paper, several types of concrete foundation blocks were prepared, including blocks of different sizes and varying reinforcement layers. Anchor bolts were arranged at the same position on the concrete foundation to allow the upper structure to be easily replaced.
4.1.2. Soil Box for Tests
The test device was constructed with a finite layer of soil. The size of the soil box is constrained by size of the laboratory, which should be large enough to still allow suitable excitations acting on the device to be produced. In designing the size of the soil box, pressure was taken into account, and the final dimensions were 2 m × 1.5 m × 1.5 m. The inner side of the soil box was covered with a layer of rubber to prevent leakage of soil and provide a better lateral restraint. In addition, the lateral restraint prevents outward expansion of the soil and absorbs lateral boundary waves, thereby simulating a realistic soil boundary.
4.2. Test Analysis
The test system included a gantry structureconcrete foundationsoil system test stand, acceleration sensor, LMS test equipment and analysis software, and force hammer, as shown in Figure 6. For the experiments, three acceleration sensors were arranged on both sides of the base (an area of relatively high rigidity), four at the connection beams and crossbeams (relatively weak rigidity), three in the longitudinal direction of the column, and three in the ydirection of the bed. To reduce noise in the data, a force window was used for signal processing to ensure a low signaltonoise ratio even in the presence of noise. By hammering in three directions above the gantry pillar, different formations were obtained.
A dynamic analysis system was used to perform modal experiments, and a comparison between the simulation and experimental model is shown in Figure 7.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
(k)
(l)
From Figure 7, it can be seen that the simulation results align well with the experimental test results. Local modal shapes of the concrete foundation are only observed above 1600 Hz since the shape of the foundation used in experiments was monolithic; therefore, the modal frequency was relatively high. Moreover, since a scaled prototype model was used, the form matches that of the physical object; thus, the experiment provides preliminary preparation for field testing.
4.3. Comparison of Theory, Simulation, and Test Results
Considering Equation (18), the state vector of the zero elements in the matrix is removed to obtain the matrix , and the matrix represents the structural parameters and natural frequency of the system. For any eigenvalue of the system, there must be a nonzero solution; therefore, the characteristic equation of the system is
By solving Equation (36), the natural frequency of the system can be obtained. Moreover, by solving the transfer equation for each element, the state vector at each point, as well as the vibration mode of the system, can be obtained. Theoretical values were also calculated, and a comparison of the theoretical, simulation, and experimental results is presented in Table 1.

Both the theoretical and simulation values are larger than the experimental values (Table 1), most likely due to the influence of the artificial boundary conditions of the soil. This effect can be seen by comparing the modal frequency values obtained from theory, simulation, and experiment, which are similar. A maximum error of 3.8% was obtained; therefore, the finiteelement model and theoretical model were validated.
5. Influence of Foundation Form on Overall Dynamic Characteristics
To study the influence of the concrete foundation on the dynamic characteristics of a number of different heavyduty CNC machine tools, a mechanics model based on the multibody transfer matrix system was used. First, according to the boundary conditions, the system transfer equation was solved and the boundary state vector at time was obtained. Then, the transfer equation of the component was used to obtain the state vector at each connection point of the system at time . Incrementing by each time, the process is repeated until reaching the required time . Finally, by changing the contour dimensions of the foundation, several properties were studied including the concrete foundation material, reinforcement conditions, and soil properties. Moreover, the displacement of the machine tool nose due to the impact load was determined. Thus, the effects of various influencing factors on machining accuracy were uncovered.
5.1. Influence of Foundation Contour Dimensions on Dynamic Response
As the contour dimensions of the foundation increase under the same impact, the displacement response at the tool nose becomes smaller, as shown in Figure 8, which suggests that contours of foundations significantly affect the dynamic response of heavy CNC machine tools. Increasing the contour size of the foundation also increases the mass of the foundation, and the joint surface between the foundation and soil becomes larger. This means the rigidity of the foundation boundary conditions increases, and the displacement of the tool nose decreases.
5.2. Influence of Concrete and Reinforcement
As the size of the concrete foundation increases, the maximum value of the dynamic response at the tool tip becomes smaller, as shown in Figure 9, and the response time is also reduced. This is because when the type of concrete changes, the modulus of elasticity of the concrete foundation also changes. If the elasticity increases, the increased rigidity of the concrete foundation also increases the stiffness at the metalconcrete joint surface. Changes in the type of concrete will therefore have a greater influence on the dynamic characteristics than changes in reinforcement. However, due to the relatively large amount of concrete in the foundation, the costs associated with the concrete are much larger than those associated with reinforcement. Therefore, the concrete model and reinforcement layers should be optimized during the design process to both improve the dynamic characteristics of the heavy machine toolconcrete foundation and reduce costs.
5.3. Influence of Soil Properties
As the soil becomes more rigid, the maximum dynamic response displacement at the tool tip decreases, as shown in Figure 10. However, the number of oscillations increases. This suggests that a harder soil can have an amplifying effect on the impact force, and soft ground can provide more shock absorption. For heavyduty CNC machine tools, accuracy is important; therefore, the amplifying effects of harder soil should be avoided. Thus, selecting a suitable soil type is necessary since soil properties clearly play a key role in machining accuracy.
6. Conclusions
(1)Taking into account the influence of joint surfaces and soil properties at various locations within the machine toolfoundation system, a flexible coupled system made of steel was established based on the multibody transfer matrix method. Contour sizes of the foundation, as well as different concrete and soil properties, various mechanical models, and different forms of reinforcement were analyzed.(2)In this paper, the substructure synthesis method was used to establish an overall dynamic model, and the frequency response function was applied to identify the parameters. The residual frequency compensation function was used to reconstruct the frequency response function, obtained via testing, thereby overcoming the errors associated with incomplete test data and clutter in the signal during the operation process to effectively improve the identification accuracy.(3)To perform laboratory tests based on different machine tool structures, concrete foundations, and soil properties, a scaled model of the heavyduty CNC machine toolconcrete foundationsoil system was built. Tests were performed, and the theoretical model and simulations were validated against the experimental results.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by the Natural Science Foundation of Liaoning Province, China (grant no. 20170540431), and Beijing Natural Science Foundation (grant no. 3162003).
References
 P. Ghosh, “FLAC based numerical studies on dynamic interference of two nearby embedded machine foundations,” Geotechnical and Geological Engineering, vol. 30, no. 5, pp. 1161–1181, 2012. View at: Publisher Site  Google Scholar
 J. B. Liu, Z. Y. Wang, K. F. Zhang et al., “3D finite element analysis of large dynamic machine foundation considering soilstructure interaction,” Engineering Mechanics, vol. 19, no. 3, pp. 3435, 2002. View at: Google Scholar
 M. Z. Aşık and C. V. G. Vallabhan, “A simplified model for the analysis of machine foundations on a nonsaturated, elastic and linear soil layer,” Computers & Structures, vol. 79, no. 31, pp. 2717–2726, 2001. View at: Publisher Site  Google Scholar
 L. Cai, Y. Tian, and Z. Liu, “Application of cloud computing to simulation of a heavyduty machine tool,” International Journal of Advanced Manufacturing Technology, vol. 84, no. 1–4, pp. 291–303, 2015. View at: Publisher Site  Google Scholar
 N. D. Stanescu and S. Tabacu, “Stability analysis for NeoHookean machinetools foundations,” Wseas Transactions on Applied and Theoretical Mechanics, vol. 5, no. 1, pp. 33–44, 2010. View at: Google Scholar
 P Vivek and P. Ghosh, Dynamic Interaction of Two Nearby Machine Foundations on Homogeneous Soil, American Society of Civil Engineers, Reston, VA, USA, 2012.
 G Štimac, S Braut, and R Žigulić, “Structural optimization of turbine generator foundation with frequency constraint,” Strojarstvo, vol. 53, pp. 389–398, 2011. View at: Google Scholar
 S. Haldar and G. L. Sivakumar Babu, “Improvement of machine foundations using reinforcement,” Proceedings of the Institution of Civil EngineersGround Improvement, vol. 162, no. 4, pp. 199–204, 2009. View at: Publisher Site  Google Scholar
 U. Werner, “Derivation of a plane vibration model for electrical machines on soft machine foundations,” Forschung im Ingenieurwesen, vol. 74, no. 4, pp. 185–205, 2010. View at: Publisher Site  Google Scholar
 X. T. Rui, B. He, B. Rong et al., “Discrete time transfer matrix method for dynamics of a multirigidflexiblebody system moving in plane,” Proceedings of the Institution of Mechanical Engineers, Part K: Journal of Multibody Dynamics, vol. 223, no. 1, pp. 23–42, 2009. View at: Publisher Site  Google Scholar
 X. T. Rui, Y. U. HaiLong, H. E. Bin et al., “Finite element transfer matrix method of multibody system for naval gun vibration analysis,” Acta Armamentarii, vol. 28, no. 9, pp. 1036–1040, 2007. View at: Google Scholar
 S. Jiang, “Dynamic characteristic parameters of linear guideway joint with ball screw,” Journal of Mechanical Engineering, vol. 46, no. 1, pp. 92–99, 2010. View at: Publisher Site  Google Scholar
 L. Li, L. Cai, T. Guo et al., “Identification of characteristic parameters of joints by substructure synthesis method,” Journal of Vibration Measurement & Diagnosis, vol. 31, no. 4, pp. 439–444, 2011. View at: Google Scholar
 Z. Liu, X. Song, Y. Zhao et al., “Stiffness identification of spindletoolholder joint based on finite difference technique and residual compensation theory,” Advances in Mechanical Engineering, vol. 5, Article ID 753631, 2013. View at: Google Scholar
 Y. Zhao, X. Song, L. Cai et al., “Surface fractal topographybased contact stiffness determination of spindletoolholder joint,” Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science, vol. 230, no. 4, pp. 1989–1996, 2015. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 Yang Tian 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.