Dynamic Modeling and Experiment Research on Twin Ball Screw Feed System Considering the Joint Stiffness

It is of great significance to study the dynamic characteristics of twin ball screw (TBS) feed system to improve the precision of gantry-type dual-driven computer numerical control (CNC) machine tools. In this paper, an equivalent dynamic model of the TBS feed system is established utilizing lumped mass method considering the stiffness of joints. Equivalent axial stiffness of screw-nut joints and bearing joints are both calculated by Hertz contact theory. Furthermore, a friction model is proposed because the friction force of the screw nut affects the stiffness of the joints. Then, the friction parameters are obtained by using the nonlinear system identification method. Meanwhile, a finite element model (FEM) is developed to assess the dynamic characteristics of TBS feed system under the stiffness of joints. Finally, validation experiments are conducted, and the results show that the positions of the nut and the velocities of worktable greatly affect the dynamic characteristics of the TBS feed system. Compared with the theoretical calculation, FEM and experiments indicate that the dynamic modeling proposed in this article can reach a higher accuracy.


Introduction
Twin ball screw (TBS) redundant feed system has been widely used in precision computer numerical control (CNC) equipment with high accuracy, high rigidity and high stability.Its dynamic characteristics have great effects on the control performance and machining precision and stability [1].This configuration can also improve the control bandwidth, speed and accuracy for large gantry-type machine tools, so studying the dynamic characteristics of the TBS system plays an important role in improving the positioning and machining accuracy [2][3][4][5][6].As per the platform shown in Figure 1a, it is considered beneficial to drive the worktable of a large single feed axis with two ball screws/motors which are arranged in parallel to provide a joint thrust in the same direction.This would significantly increase the rigidity of system and overall thrust for driving the stage.Although this layout has potential advantages, it does not come without additional concerns.For the TBS system, much of the performance is applied as important indexes which are always of concern by engineers and researchers, such as nature frequencies, tracking errors, positioning accuracy and non-synchronous errors.In the design phase, in order to investigate and predict the dynamic characteristics as accurately as possible for machine tools, a volume of dynamic models have been proposed in research.García-Herreros I et al. [7] established a lumped parameter model of ETEL gantry-type dual-driven worktable using the Lagrange's equation, but the simplified model ignored the joint stiffness and damping.Experiments show that the model's non-synchronous error is large.Liu Y [8] established the second-order wave equation based on the transmission chain of single screw fixed form under different static and dynamic stiffness of the model.The correctness of the model was verified through experiments, but the model of the bearing components was equivalent to rigid bodies without considering the joint stiffness of the supporting bearing.John EH et al. [9] found that the elastic interaction between the machine tool and the specimen greatly affects the stiffness of machine tool through theoretical calculation and more than 200 tests.Li YX et al. [10] have studied the reasons of non-synchronous error in the TBS system on gantry-type machine tools and proposed a synchronous control model based on dynamical theory which described the influence of mechanical coupling and biaxial friction characteristics on non-synchronous error.The experiment verifies friction force has a great effect on dynamic characteristics.Cheng Y et al. [11] developed the dynamic model of the translational and torsional pendulum motion of the TBS system using the motion superposition method and obtained the influence of the system's construction parameters such as stiffness ratio, friction damping ratio and the movement of the center of gravity of the mechanical coupling part on the non-synchronous error.Xie LM et al. [12] studied the load movement and transmission stiffness of cross beam of gantry-type machine tool and obtained the quantitative relationship between the displacement of cross beam load center of gravity and the system stiffness on the non-synchronous error of the TBS feed system.However, all the literature above didn't consider the joint stiffness of the contact surface and process parameters including the displacement, velocity, acceleration etc., which seriously affect the dynamic characteristics.
Many research suggests that approximately 60% of the total dynamic stiffness and about 90% of the total damping of the entire machine tool structure originates in the joints [13][14][15].Therefore, it is very important to study the joint stiffness of the transmission chain joints of the TBS system.There are so many theoretical research methods on studying the joint stiffness and damping, such as Hertz Contact theory, Takashi Yoshimura theory and Virtual Material theory etc. Zhang HJ et al. [16] used the lumped mass method to establish dynamic model of high speed feed system considering the joint stiffness, calculated the joint stiffness by using Hertz contact theory, and discussed the effect of different factors consisting of axial force, friction force and preload on stiffness of the transmission chain, but these conclusions were not extended to the TBS system.Chen et al. [17] developed a dynamic model of the high-speed ball screw feed system.The screw-nut joint is equivalent to a constant linear spring in the model.Zhu JM et al. [18] obtained the fixed and rolling joints stiffness of the feed system by Takashi Yoshimura and Hertz contact theory, respectively, and established a finite element model (FEM).The result shows high accuracy of the model.Apart from the above factors which affect the dynamic characteristics of the TBS system greatly, the coupling force caused by the effect of the coupling beam is also important.Yao WS et al. [19] constructed a dynamic model of the controlled servo mechanism by using the multi-variable identification method considering the mechanical coupling of twin parallel-axis linear servo mechanism.It can be concluded that synchronous control of high speed machining with mechanical coupling can be achieved by the proposed system model.In this paper, the mechanical coupling is considered, and the coupling force is decoupled as the screw-nut friction force.Meanwhile, the friction model is proposed, and the model parameters are obtained through nonlinear system identification.
Obviously, the state of the machine tools varies with the machining process.The different states will lead to different static and dynamic characteristics of the machine tools.FEM and Modal test are both popular methods to predict the static and dynamic characteristics.Altintas et al. [20] summarized the use of finite element analysis of the dynamic behavior of machine tools in a review.Brussel HV et al. [21] studied the dynamics of a three-axis machine tool by using FEM and found that it possessed position-dependent dynamics.Modal test based on experiments can be further divided into experimental modal analysis (EMA) and operational modal analysis (OMA).All these researches above have made great contributions to the understanding of the dynamic characteristics of a ball screw feed system.However, most of the researches are about the dynamic modeling of a single ball screw system, and few are about the dynamic characteristics of the TBS system.In this paper, the EMA of worktable at different positions and the OMA of worktable at different speeds are both conducted without any load.
The nature frequency is an important performance index of the TBS system, and the system bandwidth and the anti-jamming property can be improved by increasing it [22].Dynamic models are commonly used to analyze the nature frequency.In this paper, the proposed model and FEM are applied in the analysis of the 1st order nature frequency.This paper reveals that joint stiffness and friction force of screw-nut have great effects on dynamic characteristics of the TBS feed system.Therefore, a lumped dynamic model for studying the transmission stiffness of the TBS feed system is established.
Meanwhile, a friction model is proposed to study the influence of the friction force on screw-nut joint stiffness.In order to conduct experiments, a tested worktable with twin ball screws and dual linear guides are constructed.The axial stiffness of the supporting bearings and the double pre-tightening nut is tested.The vibration tests are carried out to measure the first nature frequency of the TBS system, and the relationship between friction force and the first nature frequency is obtained.Finally, contrastive analysis between the experimental results and the simulated results of three models is conducted.The results show that the proposed model agrees significantly moreso with the experimental results than the discrete model and hybrid model do.

Dynamic Equivalent Model of the TBS System
Figure 1b shows the two-dimensional dual-driven worktable studied in this paper, the X and Y directions adopt dual-driven structure, which improves the transmission rigidity of the two-feed direction.Taking X direction as an example, the typical structure mainly consists of dual servo motors, twin ball screws, dual linear guide rails, bearings, worktable and other parts.
The transmission stiffness of the worktable depends on the stiffness of the dual flexible screw pairs and the joints, including the fixed joint and rolling joint.The fixed joint consists of worktable-slider joint and guide-bed joint which are connected by bolts.The rolling joint contains a screw-nut joint, bearing joint and guide-slider joint.Many studies [15] indicate that radial stiffness has little effect on the overall stiffness, therefore, the axial stiffness of the feed system is only considered in this paper.
Due to coupling of the worktable, the preload between screw and nut varies with coupling force and determines the screw-nut joint stiffness.The coupling force is mainly caused by friction force of the screw-nut joint, so the effect of friction between the screw and nut should not be neglected.In order to carry out research easily, the friction between the guides and sliders is not considered.In this section, a nonlinear friction model is developed, which deduces the relationship between the worktable speed and friction force, and the influence of speed on dynamic characteristics of the TBS system is further discussed.
In the equivalent dynamic model as shown in Figure 2, the screw, the screw-nut joint and the bearing joint are equivalent to spring damping elements.The worktable is equivalent to a lumped mass model.In the case of ignoring the influence of servo system stiffness, the transmission stiffness of the TBS system depends on joint stiffness.Where, C b1 and C b2 are the damping of the bearings of the twin axes.K ls1 , K rs1 and C ls1 , C rs1 are the equivalent stiffness and damping of the left and right sides of the X 1 shaft's nut, respectively.K ls2 , K rs2 and C ls2 , C rs2 are the equivalent stiffness and damping of the left and right sides of the X 2 shaft's nut, respectively.K nut1 , K nut2 and C nut1 , C nut2 are the screw-nut joint stiffness and damping of the X 1 and X 2 shaft, respectively.

Dynamic Modeling of the TBS System
Figure 2 shows the model of the TBS feed system which can be formulated based on Lagrange's energy method [23][24][25].First, the kinetic energy of the system can be computed as follows: The potential energy of the feed system is given by Let L = T − V, and The general coordinates and general forces are defined as where, F a1 , F a2 are the axial forces of the twin ball screw, respectively.T 1 , T 2 are the output torque of the dual servo motor, respectively.The dynamic equation of the system can be expressed matrix form [26].
where, M, C and K are, respectively, the mass matrix, damping matrix, and stiffness matrix of the system.For simplifying analysis, the case is considered without damping.
Table 1 lists the parameters and related values used in the model.In this paper, the first order mode is mainly studied.It's obvious that the first and second modes are from axial vibration, and the third mode is from torsional vibration.Therefore, the torsional stiffness of the TBS system is neglected.Eigenvectors can be studied to clarify the modeling.Many other parameters are obtained in the next section.

Transmission Stiffness Modeling of the TBS System
In this section, the calculation of the stiffness of the screw-nut joints and bearing joints and contact parameters identification are discussed.

Calculating the Stiffness of the Screw-Nut Joints
The ball screw pair is mainly composed of ball screw, nut and other parts.In order to eliminate the axial gap and improve the axial contact stiffness of the ball screw, the double-nut pretension structure is adopted in this paper.As shown in Figure 3, the double nut is composed of nut A and nut B which are connected by an elastic spacer.The elastic spacer can provide preload by elastic deformation.
Assuming that the contact between the ball and the raceway obeys Hertz contact theory [27], therefore, the screw-nut joint stiffness can be calculated by Hertz contact theory.The ball screw used in this paper is HIWIN FDW20-5B2, and its basic parameters are shown in Table 1.Without considering the influence of the axial force of the screw, we can get Equation (7).
The synthetic curvature of the contact point between the ball and the raceway is: The axial deformation of the screw nut and the axial stiffness of the joints can be obtained from the geometric relation and Hooke's law:

Calculating the Axial Stiffness of Screw-Nut Joints Caused by Friction Force
In the process of the worktable movement, the screw-nut joint stiffness varies with the displacement and velocity of the worktable due to coulomb friction and axial force influence between the screw and nut.As shown in Figure 4, F a1 and F a2 are the axial forces of the axis 1 and axis 2, respectively.When the worktable is in uniform motion, the axial force of the screw depends on the friction force of the screw nut.According to the second law of Newton, the Equation ( 10) is obtained as follows: ( When a = 0, the axial force of the screw is equal to the friction force.Furthermore, the friction model is established and the relationship between the friction force and the output torque is [29]: where, T out and P T , P T = 0.90 are the output torque and efficiency of the motor, respectively.In the uniform motion of the worktable, there is Take axis 1 as the main axis, and the axis 2 as the slave axis.When the velocity of worktable is in the range of 0-50 mm/s, the output torque of the dual motors is synchronously measured.The relationship between the friction force of the screw nut and the velocity of worktable can be obtained by the nonlinear identification method.As shown in Figure 4, the nonlinear fitting curves of the friction force and the velocity are given.Table 2 shows the friction model parameters of two axis obtained by nonlinear identification.Considering the friction force of the screw nut, the Equation ( 12) can be obtained: Similarly, the axial deformation of the nut A and the joint stiffness can be obtained: According to the coordination condition of elastic deformation, the deformation relation of nut A and nut B is given by Equation ( 14): Therefore, the contact deformation of nut B: As it can be seen from the above Equation ( 15) that when δ AT ≥ 2δ A , the contact deformation of the nut B is 0, which causes a sharp change of the screw-nut joint stiffness.Similarly, the equivalent axial contact stiffness of the nut B is: Then the equivalent axial stiffness of screw nut of the axis 1 depends on the larger value of K nutAT and K nutBT , expressed as:

Calculating the Axial Stiffness of the Bearing Joints
Figure 5 shows the bearing assembly used in this paper.The ball screw is mounted with one end fixed and one end free, so the fixed end with a pair of angular contact ball bearings can bear axial force.The deep groove ball bearings with only supporting functionality are used at the other end of the screw.Therefore, the axial stiffness can be ignored.F p and α are the preload and contact angle of the angular contact bearing, respectively.Table 3 shows the basic parameters of the angular contact bearings used in this paper.The normal force and axial deformation of a single ball can be obtained by Hertz contact theory.
Therefore, the axial stiffness K 1 of a single bearing joint is: Since the contact between the two angular contact ball bearings is in series, the stiffness of bearing joints K b1 is: where, zis the number of balls, δ b1 , δ b2 are the contact deformation of the ball with the inner and outer rings, respectively.

Calculating the Transmission Stiffness of the Overall TBS System
As shown in Figure 2, the equivalent axial stiffness of the TBS system depends on the axial stiffness of the screw, the stiffness of the screw-nut joints and the bearing joints.On the left side of the nut, the angular contact bearings and the screw are in series.Thus, the equivalent axial stiffness of the left subsystem of the lead screw is: where, x is the length from the nut to the left side of the screw, and L is the working length of the screw.
On the right side of the nut, the axial stiffness of the screw can be neglected.Furthermore, the left subsystem and the nut are in series, and axis 1 and axis 2 are in parallel.The axial stiffness of axis 1, K a1 and X direction of the TBS worktable, K Ax can be expressed as: Then, the equivalent axial stiffness K Ax in the Y direction can be obtained in the same way.The stiffness of the X and Y direction is in series, and then the overall transmission stiffness of the two-dimensional TBS worktable is as follows [30]:

FE Analysis of the TBS Worktable
ANSYS finite element analysis software is applied to perform modal analysis of the TBS worktable.At first, the 3D geometry model of the TBS worktable is established in the 3D software Solidworks, and then imports it into ANSYS.In order to calculate conveniently, the model is simplified properly, and the features of the screw thread, the small holes, fillet etc., which have few influences on the dynamic characteristics, are neglected.Motors and other components that have few effects on calculation results are ignored.The solid model is modeled by Solid 186 unit.An automatic meshing method is used in ANSYS to divide the solid model.And the Combine 214 spring elements are used to simulate the stiffness of the screw-nut joints, the bearing joints and the guide-slider joints, ignoring the stiffness and damping of the fixed joints.This section mainly studies the finite element analysis of the TBS worktable considering the dynamic stiffness of the joints.

FEM of the Screw-Nut Joints
As shown in Figure 6a, the equivalent stiffness of screw-nut joints is taken as springs.Four springs are arranged along the circumferential direction of the screw.The stiffness of each spring is equal to 1/4 of the stiffness of screw-nut joints.In order to reduce the number of mesh and improve the calculation efficiency, the structure of screw and nut is simplified, so the FEM is obtained.The equivalent axial stiffness of screw-nut joints, K nut = 1.28 × 10 8 N/m can be deduced by the third section.

FEM of the Bearing Joints
As shown in Figure 6b, similarly, the equivalent stiffness of bearing joints is seen as springs.Four springs are arranged along the circumferential direction of the bearing.The stiffness of each spring is also equal to 1/4 of the stiffness of bearing joints.The equivalent axial stiffness, K b = 4.2 × 10 7 N/m can be similarly drawn from the former section, the FEM of the bearing can be obtained by simplifying the ball.

FEM of the Guide-Slider Joints
As shown in Figure 6c, the material of the guide is steel 55, the Young modulus, E = 206 GPa.The density is 7800 kg/m 3 , and the Poisson's ratio is 0.3.The equivalent stiffness of the guide-slider joints is taken as springs, and the damping effect is ignored.The dynamic model and FEM are established, and the equivalent normal stiffness and tangential stiffness of the guide-slider joints are K n = 3.2 × 10 7 N/m and K τ = 1.2 × 10 7 N/m, respectively.The two parameters can be calculated by the third section.

FEM of the TBS Worktable
The FEM of the entire TBS worktable is shown in Figure 6d.In the ANSYS, the boundary conditions are set.The bottom surface of the worktable is fixed.The joint stiffness among all components is set, which is obtained by the above value.The method of structured grid combined with unstructured grid is used to divide the mesh size to 2 mm.

Finite Element Analysis of the TBS Worktable
By solving and analyzing the FEM, the first order nature frequencies of the worktable can be obtained at different positions.As shown in Table 4, the nature frequencies obtained by the FE analysis are larger than that obtained by experiments because the stiffness of the fixed joints is ignored in FE analysis.The overall stiffness is larger than the actual; so are the nature frequencies.However, these will be improved in the following study.Block Lanczos method is used for modal solution in ANSYS when the worktable is located at the middle position of the stroke.Meanwhile, the vibration modes of two-dimensional TBS worktable are obtained.

Parameters
Given Value M w : worktable mass 24.43 kg M n1 , M n2 : nut 1 and nut 2 mass 0.575 kg p: the lead of the ball screw 5 × 10 −3 m J m1 , J m2 : inertial moment of the motor 1 and 2 1.76 × 10 −4 kgm 2 J b1 , J b2 : inertial moment of the ball screw 1 and 2 7.85 × 10 −5 kgm 2 J w : inertial moment of the worktable 0.505 kgm 2 X n1 : axial displacement of the nut 1 state variable X n2 : axial displacement of the nut 2 state variable X w : displacement of the worktable state variable θ n1 : rotating angle of the screw 1 state variable θ n2 : rotating angle of the screw 2 state variable θ w : twist angle of the worktable state variable Figure 7 shows the first four modes of the TBS worktable considering the stiffness of joints.The Figure 7a,b show the first and second low order translational vibration modes of the worktable in the Y and X direction, respectively.The Figure 7c shows the third vibration mode is high order yawing vibration of the worktable in the Z direction.The Figure 7d shows the fourth vibration mode is the pitching vibration of the worktable around Y axis.

Dynamic Experiments
The experimental study is carried out on the basis of theoretical calculation and FE analysis.According to the transmission stiffness model, the transmission stiffness of the TBS feed system is related to the position of the worktable.The modal tests of the worktable are carried out, because it is of great significance to determine the dynamic characteristics of the worktable at different positions for practical machining process.
As shown in Figure 1a, the worktable is bolted to the bed.The direction of X and Y are respectively driven by "dual servo motors + twin ball screws" to realize the linkage of two directions, and different contour trajectory.The parameters of the ball screw are shown in Table 1.Reference to the relevant information [31], the initial preload of the screw nut is (0.07~0.1)P d = 1.115KN.A pair of angular contact ball bearings are mounted on the left end of the ball screw (model: FAG B7202CT P4S.UL), supporting and bearing the axial force and the worktable of the coupling force.The axial preload of the bearing,F p = 1 KN is applied to the bearing lock nut by means of a torque wrench.The relevant parameters are shown in Table 2.
The Econ dynamic analyzer is applied to test the dynamic characteristics of the worktable at different positions and velocities.As shown in Figure 8, four single axes accelerate sensors (model: BK4507B) are arranged at four corners with large vibration amplitude of the worktable.The impulse force hammer (model: DYTRAN 5800B5) is utilized by single point excitation multi-point vibration pickup method to conduct the experiments.Econ MI7008 acquisition card is employed to realize high-precision data acquisition.Data processing is carried out by Modal Genius software.The frequency response curve and mode shape are obtained by data fitting.

The Dynamic Characteristics of the Worktable at Different Positions
Because the transmission stiffness of the TBS system is related to the nut position, its dynamic characteristics will change when the nut position varies in the X and Y direction of the worktable.Figure 9 shows the relationship between the transmission stiffness of the worktable and the position of the nuts.The abscissa indicates the ratio of the nut position to the screw stroke, and the longitudinal coordinate is the overall transmission stiffness of the worktable.It can be seen from the Figure 9 that the overall transmission stiffness of the TBS system decreases with the increase of the nut position in the two directions.The transmission stiffness is maximum when the nut is located at the left end of the ball screw.Similarly, the nut at the right end of the transmission system shows the minimum transmission stiffness.So, the overall stiffness of the system is poor, which is not conducive to the processing operation.Therefore, the case that the nut is at the right end of the ball screw should be avoided in the actual processing, especially for high precision machine tools.In order to make the experiment convenient and without losing generality, only the position variation of the nut in the X direction is taken into account in the experiment, and the nut is located at the middle point of the Y direction.Figure 10 shows the frequency response curves of the TBS worktable by experiments when the nut is located at different positions in the X direction.It can be seen from the curves that the first order nature frequency increases with the increase of the overall stiffness of the worktable.From Table 5 and Figure 11, the first nature frequencies of the TBS worktable are very close when the nut is located at different positions because the screw stroke is too short.The theoretical calculation is closer to the experimental data, and the errors are basically within 10%, which verifies the correctness of the theoretical model.However, the FE analysis has larger errors than the experimental results do.The maximum error is up to 19.2% because the worktable structure is simplified properly in the process of establishing the FEM, such as damping, chamfer and small pore structure.These factors lead to big errors, but they will be improved in the future work.

The Dynamic Characteristics of the Worktable at Different Velocities
The stiffness of the screw-nut joints and the whole TBS system are changed due to the different friction characteristics of the screw nut at different velocities.For the convenience of research, the X direction of the system is taken as the research object, and the middle point of the screw stroke is selected as the reference position.In the process of the experiments, the master-slave synchronization control strategy is adopted in the TBS system.The axis 1 and axis 2 are regarded as the main axis and slave axis, respectively.The output velocity of the master serves as the velocity reference for the slaves.
The OMA reflects the dynamic characteristics of the machine tools under different working states.In the process of nonlinear friction force identification, the dynamic characteristics of nuts at the middle point of the ball screw at different velocities are studied, without external excitation, and the excitation is produced by the servo motor to drive the ball screw, which is called self-excited.The experimental device is shown in the last section.
The equivalent stiffness of the screw-nut joints, K nutAT and K nutBT at different velocities can be calculated by Equations ( 10)- (17).As shown in Figure 12, it can be seen that the stiffness of the screw-nut joints decreases with the increase of the velocity of the worktable, and the stiffness of the screw nut B reduces faster.When the velocity is 35 m/min, the stiffness of the screw nut A remains stable, but the stiffness of the screw nut B decreases sharply to 0. This is because the friction force of the screw nut B is equal to the driving force applied on the nut B at this velocity, which means that the velocity is the maximum critical velocity of the TBS system.Meanwhile, when the velocity exceeds 35 m/min, the nut B will fail and does not play the role in providing preload and eliminating clearance.
The equivalent stiffness of the bearing joints, K b under different velocities can be calculated by Equations ( 18)- (20).As shown in Figure 11, when the velocity is in the range of 0-30 m/min, the equivalent stiffness of the bearing joints decreases slowly.When the velocity is 30 m/min, the equivalent stiffness of the bearing joints is minimum.This is because the preload is equal to the axial force of the bearing.When the velocity exceeds 30 m/min, the stiffness of bearing joints increases rapidly because the axial force of the bearing is greater than the preload.Figure 13 shows the frequency response curves at different velocities.It can be seen from the curves that when the velocity is in the range of 0-30 m/min, the first order nature frequency of the TBS system decreases because the overall stiffness of the TBS system is minimum at this time.Therefore, this velocity should be avoided in the actual machining process.Furthermore, when the velocity is over 30 m/min, the nature frequency of the system increases gradually due to the increase of the overall stiffness of the TBS system.By comparing Figures 10 and 12, it is found that the dynamic nature frequency is higher than the static nature frequency.

Conclusions
This article considers the stiffness of the typical joints of the TBS feed system, takes the self-developed TBS two-dimension worktable as the research object, and establishes the equivalent transmission stiffness model of the TBS worktable.What's more, the friction force model of the screw nut is proposed.Finally, the dynamic characteristics of the two-dimensional TBS worktable at different positions and at different velocities are studied, and the effectiveness of the model is verified by the experiments.The main conclusions of this article are as follows: 1.
The lumped mass method and FE method were employed to establish the dynamic model of the TBS feed system considering the stiffness of screw-nut and bearing joints, and to build the stiffness model of the transmission chain.Furthermore, the calculation formulas of the static overall stiffness of the two-dimensional TBS system were deduced.It also gives the mapping relation between the overall stiffness and the position of the nut in the X and Y direction.

2.
The equivalent stiffness model of the joints is established by using Hertz contact theory, and the stiffness of the screw-nut joints, bearing joints and guide-slider joints were calculated by it.Then, the calculation results of the joints' stiffness were applied to the FEM.Compared with theoretical calculation and FE analysis, the maximum error is up to 19.2% because the stiffness of the fixed joints was neglected.

3.
The dynamic characteristics of the worktable at different positions were studied, and the maximum error of the nature frequency of the system was 11.1%.Thus, the correctness of the theoretical model is verified by experiments.The overall stiffness of the TBS system decreases with the increase of the nut positions of X and Y direction.In actual machining process, the furthermost right side of the worktable should be avoided.

4.
The dynamic characteristics of the worktable at different velocities are studied.Due to the friction force variations with the velocities, a friction model was proposed, and the parameters of the model were obtained through nonlinear system identification.The experiment showed that the maximum critical speed of the system is 35 m/min.When the velocity exceeds it, K nutBT = 0, the bearing will fail.When the velocity is 30 m/min, the stiffness of the bearing is minimum and the same as the overall TBS system.By comparing Figures 10 and 13, it shows the dynamic nature frequency is larger than the static nature frequency.
The main implication of this paper was developing the transmission chain stiffness model and researching the EMA and OMA of the TBS feed system considering the stiffness of the joints, which will guide the dynamic design of the CNC machine tools.In further study, the stiffness of the fixed joints such as worktable-slider joints connected by bolts will be considered.The load and cutting force will also be considered.Furthermore, the dynamic characteristics at different accelerations will be studied.

Figure 1 .
Figure 1.(a) Dual-drive grinding wheel machine tools; (b) The 3D model of the TBS feed worktable.

Figure 2 .
Figure 2. Equivalent dynamic model of the TBS feed system of the X direction.

d b cos λ 2 sin 5 β cos 5 λ 1
where, P p is the axial preload of ball screw, P A and P B are the normal forces between the ball and the nut A and B of the screw, respectively.N is the number of working balls and, N = i•πd/d b cos λ. i is the total number of load bearing ring of the single nut.d and d b are the nominal diameter and ball diameter of the ball screw, respectively.And β is the contact angle between the raceway and the ball.λ is the helix angle of the ball screw.K i and a i are Hertz coefficients, which can be obtained from the literature[28].δ A1 and δ A2 are the contact deformation of the nut A. Similarly, δ B1 and δ B2 are the contact deformation of the nut B. µ 1 , µ 2 and E 1 , E 2 are the Poisson's ratio and elastic modulus of ball and screw, respectively.Where, µ 1 = µ 2 = 0.3, E 1 = E 2 = 210 GPa.f is the tightness of the ball screw.

Figure 3 .
Figure 3. Cross-section view of a gasket-type double-nut ball screw joints.

Figure 4 .
Figure 4. Friction force and motor velocity curve of the TBS system.(a) Friction force and motor velocity curve of axis 1; (b) Friction force and motor velocity curve of axis 2.

Figure 5 .
Figure 5. Assembly structure of the support bearing units.

Figure 6 .
Figure 6.FEM of the TBS worktable.(a) FEM of the ball screw and nut; (b) FEM of the angular contact bearing; (c) FEM model of the linear guide; (d) FEM model of the whole TBS system.

Figure 7 .
Figure 7.The typical mode of the TBS system.(a) The first mode: linear vibration; (b) The second mode: linear vibration; (c) The third mode: yawing vibration; (d) The fourth mode: pitching vibration.

Figure 8 .
Figure 8. Dynamics measurement of the TBS feed system.

Figure 9 .
Figure 9.The transmission Stiffness of the dual-screw system.

Figure 10 .
Figure 10.Frequency response of the TBS system under the different nut position.

Figure 11 .
Figure 11.The 1st nature frequency of the TBS system at different nut position.

Figure 12 .
Figure 12.The equivalent stiffness of the joints.

Figure 13 .
Figure 13.Frequency response of the TBS system under the different velocity.

Table 1 .
Basic parameters of the screw-nut joints.

Table 2 .
The coefficients of the friction model.

Table 3 .
Basic parameters of the angular contact bearing.

Table 4 .
Parameters and related values used in modeling the TBS system.

Table 5 .
Comparison of theoretical calculation, finite element analysis and experiment.