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 rigid-flexible coupled heavy-duty CNC machine tool-foundation 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 heavy-duty machine-foundation 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 heavy-duty 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 heavy-duty 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 mechanical-foundation system, it is necessary to establish an accurate system model. Previous research has been carried out on mechanical models of dynamic mechanical-foundation systems. Ghosh [1] used the three-dimensional finite-difference method to establish a dynamic mechanical-foundation 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 three-dimensional finite-element model was proposed by Liu et al. [2] that considered a soil-power machine interaction system. Based on this model, the main factors influencing the dynamic response of the large-scale machine tool-foundation 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 foundation-foundation systems using cloud computing. Furthermore, Stanescu and Tabacu [5] considered the nonlinearity of soil and established a two-degree-of-freedom (DOF) machine-foundation model to study the stability of the equilibrium position.

The purpose of analyzing the characteristics of power machinery-foundation 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 finite-element 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 finite-element model of a turbogenerator-foundation 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 soil-foundation 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 machine-foundation 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 mechanical-foundation 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 heavy-duty CNC machine tool-foundation 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 tool-foundation 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 heavy-duty 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 finite-element method, the total stiffness matrix can be effectively reduced. A multibody mechanical model was established according to the structural characteristics of the heavy-duty CNC machine tool and its foundations by considering the upper and lower structural subsystems. A dynamic equation can be established to describe the entire heavy-duty CNC machine-foundation system according to the state vector of the subsystem connection interface. A multibody dynamic model of the heavy-duty CNC machine tool-foundation 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 large-scale 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 z-directions; are the equivalent damping in x-, y-, and z-directions; are the equivalent stiffness about the rotational displacement around the x-, y-, and z-directions; and are the equivalent damping about the rotational displacement around the x-, y-, and z-directions.

2.2. Rigid Body Elements

During the modeling process, structural members of the heavy-duty CNC machine tool-foundation 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 Heavy-Duty CNC Machine Tool-Foundation 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 large-scale characteristics, concrete foundations of heavy-duty 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 finite-element 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, metal-to-metal joints and the metal-to-concrete 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 metal-to-metal and concrete-to-concrete interfaces.

3.1. Dynamic Modeling of the Joint Surface

Heavy machine tools can be considered large fixed-weight 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 metal-concrete 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 metal-concrete test system is shown in Figure 5.

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 low-order term, , and high-order 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 n-ary 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 different-sized foundation, foundation forms, and soils on the dynamic characteristics of the system. Moreover, the characteristics of heavy-duty CNC machine tools make it impossible to carry out dynamic testing. Therefore, it is necessary to establish laboratory models of heavy-duty gantry CNC machine-foundations. The laboratory-based test device can solve a number of issues: wasted capital caused by the incompatibility of heavy-duty gantry CNC machine-foundation systems is reduced; the influence of joint surfaces between the structural members of the heavy-duty 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 heavy-duty gantry CNC machine tool-foundation 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 z-direction. 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 structure-concrete foundation-soil 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 y-direction of the bed. To reduce noise in the data, a force window was used for signal processing to ensure a low signal-to-noise 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.

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 finite-element 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 heavy-duty 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 metal-concrete 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 tool-concrete 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 heavy-duty 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 tool-foundation 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 heavy-duty CNC machine tool-concrete foundation-soil 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.


This work was supported by the Natural Science Foundation of Liaoning Province, China (grant no. 20170540431), and Beijing Natural Science Foundation (grant no. 3162003).