Specifications of the test bench.

## Abstract

Ball screw feed drive system is the most widely used linear drive system in the field of industrial automation. The continuous search for efficiency puts forward higher requests to the machine tool for high speed and high acceleration, which makes the feed drive system of lightweight-designed and large-size machine tools more likely to produce vibration during high-speed and high-acceleration feed operation. Electromechanical co-simulation for ball screw feed drive dynamics is an important technique for solving vibration problems occurring in the feed motion. This chapter elaborates on this technology from three aspects: modeling and simulation of dynamic characteristics of ball screw feed drive, modeling and simulation of servo control system, and the electromechanical co-simulation of ball screw feed drive system. In this chapter, the basic theoretical models, the establishment of simulation models and the comparison between simulation and experiment results of ball screw feed drive system are comprehensively introduced to provide technical references for readers.

### Keywords

- ball screw feed drive system
- dynamic characteristics
- electromechanical co-simulation
- vibration
- lumped mass model

## 1. Introduction

Ball screw feed system is the most widely used linear drive system in the field of industrial automation [1]. In order to enhance the speed and accuracy of present systems further, current research focuses on the vibration reduction and avoidance of the feed drive. Additional damping modules or structures are integrated in the feed drive system to achieve this goal, such as semi-active damping system, set point filtering, etc. Active damping system only reacts once a vibration is present, and set point filtering can lead to path deformation [2, 3, 4]. Another way to solve this problem is to generate a smoother trajectory. For this purpose, numbers of trajectory algorithms are found out, and the frequency contents of the trajectory orders are discussed and compared [5, 6]. The vibration caused by the trajectory is difficult to analyze on hardware because of the coupling factor of variety excitation sources. All these researches need a simulation method to help the researchers or engineers study or optimize the design and parameter setting of the feed drive system [10].

Finite element model of ball screw feed drive system can predict the accurate dynamic characteristics. However, it is difficult to integrate with the simulation model of servo control system. Lumped parameter model of ball screw feed drive system can simplify the simulation model by reducing the number of degrees of freedom (DOF) of the whole system. More importantly, it can easily integrate with the simulation model of the servo control system. A reasonable simplification of the lumped parameter model is the key to accurately predict the vibration of feed drive system [7, 8, 9].

In this chapter an electromechanical co-simulation method for ball screw feed drive system was established, which can be used to study the dynamic characteristics and vibration behavior of the feed drive system. An optimized dynamic modeling and simulation method of a ball screw feed drive based on the lumped mass model was firstly presented, and the optimized calculation method of the equivalent parameters was given. Then, a model of servo control system was built up, and based on it, the electromechanical co-simulation of ball screw feed drive system was established. Finally, a simulative and experimental test is conducted based on a ball screw feed drive system test bench. The result shows that electromechanical co-simulation of ball screw feed drive system could achieve a very good predictability.

## 2. Dynamic characteristic modeling and simulation of ball screw feed system

### 2.1. Lumped mass model of ball screw feed system

A typical ball screw feed system consists of a servomotor, coupling, ball screw, work table, and base (Figure 1). The ball screw is supported by two sets of bearing, which are fixed to the base. The servomotor torque is transmitted through a coupling onto the ball screw shaft to drive the work table. The linear guideway constrains the movement of the work table in an axial direction. The base is fixed on the machine bed or placed on the ground. The transformation from the rotational movement of the screw shaft into the linear motion of the work table is realized by the ball screw system with its transmission ratio

Low-order modes are the main factors affecting the dynamic characteristics of the ball screw feed drive system of machine tools. Typically, the first axial and rotational modes of the ball screw show a dominant influence on the overall dynamics, while the relevance of higher-order modes for most technical applications is rather small [8].

The lumped mass model can reasonably reduce the number of degrees of freedom (DOF) of the simulation model while preserving the low-order modes of the system to simplify calculations. Figure 2 shows the lumped mass model of a ball screw feed drive system. The influence of the shaft on the rotational mode and axial mode of the drive system is explicitly included into the lumped mass model here. Therefore, the shaft is separated into two different branches, an axial branch and a rotational branch, while the coupling once more is realized using constrained equations. Since all components are expressed by discrete springs and dampers, the rigidity values of shaft, coupling, and bearing are combined to an overall axial

In this model the inertial component parameters are defined as the following: rotary inertia of servomotor

The equivalent rigidity parameters in the model are defined as the following: equivalent torsional rigidity

The equivalent damping parameters are defined as the following: servomotor torsional damping

The DOF parameters of the lumped mass model are defined as the following: angular rotation of the servomotor

The speed parameters of equivalent damping are defined as the following: servomotor equivalent damping speed

According to the Lagrange’s equations of the second kind, the dynamic model of the ball screw feed drive system is built up. The total kinetic energy

According to the definition of the system lumped mass, we have the independent coordinates system

The force inputs of the ball screw feed system are the servomotor torque

With

where,

The dynamic model of the ball screw feed system shown in Eq. (8) was decomposed into three subsystems: screw shaft torsional vibration system, screw shaft axial vibration system, and the table vibration system. The simulation model of the ball screw feed system can be established as Figure 3. The input of the simulation model is the motor torque

### 2.2. Equivalent parameter calculation method of ball screw feed system lumped mass model

Accurate ball screw feed system dynamic model requires a reasonable equivalent parameter calculation method of the lumped mass model. As mentioned the shaft has influence on the rotational mode and the axial mode of the drive system; the shaft is separated into two different branches, an axial branch and a rotational branch, while the coupling is realized using constrained equations. The dynamic characteristics of the feed system should be analyzed to select the appropriate equivalent parameter calculation method. The inertia of the axial system and the inertia of rotational system are not only the mass or inertia of the component itself but also the mass or inertia converted to the independent coordinate system component of the dynamic system.

The equivalent rotary inertia of the screw

With the material density

The screw equivalent mass

The axial rigidity of the ball screw feed system is related to the installation method of the screw. Here is an example of the screw-fixed-support method used on most machine tools (Figure 4). Servomotor side of the screw shaft uses fixed support to provide screw axial support, and the end of the shaft is free support. Therefore, the axial rigidity

With the material elastic modulus

The torsional rigidity of the ball screw feed system

With the nut reference rigidity

## 3. Servo control system modeling and simulation

### 3.1. Modeling of permanent magnet synchronous motor

Permanent magnet synchronous motors can be divided into two types according to the rotor type, salient pole rotor and non-salient pole rotor. The structure is shown in Figure 5; in the surface-mounted permanent magnet synchronous motor (Figure 5(a)), the magnetic circuit of the rotor is symmetrical, and the magnetic permeability and air gap permeability of the permanent magnet material are approximately the same. In the rotor two-phase coordinate system, the direct-axis inductance and the quadrature-axis inductance are equal, that is

The rotor magnetic paths of plug-in-type (Figure 5(b)) and built-in-type (Figure 5(c)) permanent magnet synchronous motors are asymmetrical, and the quadrature-axis inductance is greater than the direct-axis inductance, that is

Taking non-salient pole rotor permanent magnet synchronous motor as an example, we simplify the motor model with the following conditions: neglecting the saturation of the motor core; no eddy current and hysteresis loss; permanent magnet material has zero conductivity; three-phase windings are symmetrical; and induced potential in the winding is sinusoidal. Then, a schematic diagram of the physical model of the motor shown in Figure 6 can be obtained.

The axis of the sinusoidal magnetomotive wave generated by a flowing forward current through the phase winding is defined as the axis of the phase winding. Take axis A as the spatial reference coordinate of the ABC coordinate system. It is assumed that the positive direction of the induced electromotive force is opposite to the positive direction of the current (motor principle); take the counterclockwise direction as the positive direction of the speed and electromagnetic torque, and the positive direction of the load torque is the opposite. The physical model equation of permanent magnet synchronous motor is as follows:

In Formula (18),

### 3.2. Coordinate transformation

From the above physical model, we can see that in the ABC coordinate system, the PMSM rotor is asymmetric in the magnetic and electrical structures. The motor equation is a set of nonlinear time-varying equations related to the instantaneous position of the rotor, which makes the analysis of the dynamic characteristics of the PMSM very difficult. It is usually necessary to convert the motor equations by coordinate transformation to facilitate analysis and calculation. The coordinate system used in the vector control of the permanent magnet synchronous motor and their relationship is shown in Figure 7. In the figure, the

#### 3.2.1. Clarke transformation

Clarke transformation simplifies the voltage loop equations on the original three-phase windings into the voltage loop equations on the two-phase windings, from the three-phase stator ABC coordinate system to the two-phase stator

#### 3.2.2. Park transformation

In the physical sense, the Park transformation is equivalent to projecting the currents

#### 3.2.3. Mathematical model of permanent magnet synchronous motor in coordinate system

The mathematical model of the permanent magnet synchronous motor in the ABC coordinate system can be transformed into any two-phase coordinate system through coordinate transformation, so that it is possible to simplify the decoupling of the motor flux linkage equation and the electromagnetic torque equation. If the mathematical model of the motor is transformed into a

With the equivalent flux

Take the rotor permanent magnet flux linkage

where

The electromagnetic torque equation in the *p* is the number of pole pairs,

Using the components of the

### 3.3. Servo control modeling of ball screw feed system

In order to model the servo control of the ball screw feed system, the modeling of the three-loop cascade control architecture of the vector control and servo control system of the permanent magnet synchronous servomotor is studied, which is commonly used in the ball screw feed system.

#### 3.3.1. Modeling of permanent magnet synchronous motor vector control.

The principle of space vector pulse width modulation (SVPWM) is based on vector equivalents. The magnitude and direction of the current vector can be indirectly controlled by the timing of the six switching elements of the inverter through the three-phase winding of the permanent magnet synchronous motor, so that the winding produces a constant amplitude circular magnetic field that rotates according to a given demand, thus dragging the permanent magnet to rotate. The voltage inverter circuit is shown in Figure 8. The simulation model of SVPWM control system for permanent magnet synchronous motor is shown in Figure 9.

#### 3.3.2. Modeling of cascade control system

The SVPWM control system for permanent magnet synchronous motor is based on the three-phase current information and rotor position information fed back by the motor. AC motor is equivalent to a direct current motor by formula transformation to control the position and amplitude of the stator current.

The control system schematic is shown in Figure 10(a). The system includes a cascaded control structure with a P-position controller, a PI-velocity controller, and a PI-current controller. In the cascade control system, the servomotor feedback speed

### 3.4. Electromechanical co-simulation modeling of ball screw feed drive system

Based on the lumped mass model of ball screw feed system and the servo control system simulation model, an electromechanical co-simulation model of the ball screw feed drive system was constructed. The co-simulation schematic is shown in Figure 10; as described above (a) is the semi-closed-loop cascade control system simulation model, while (b) is the lumped mass model of ball screw feed system. The inverter outputs ABC three-phase voltage to servomotor stator, which generates rotating magnetic field and produces magnetic torque on the servomotor rotor. This magnetic torque is the output torque

The electromechanical co-simulation model of the ball screw feed drive system is shown in Figure 11. The S_Cal module on the left side generates the trajectory command for the feed drive system according to the acceleration/deceleration strategy. Under the cascade control system, which consists of position controller, velocity controller, and current controller, the servomotor drive and the ball screw accomplish the motion command accordingly.

## 4. Experimental verification of the electromechanical co-simulation model of ball screw feed drive system

The electromechanical co-simulation model in this chapter has been tested on a single-axis ball screw drive system test bench shown in Figure 12. The test bench uses an i5 CNC system and servo system of Shenyang Machine Group, which use a semi-closed-loop cascade control structure. The specifications of the test bench are listed in Table 1, which are either obtained from the manufacturers’ catalogs, approximated from prior knowledge, or calculated from computer-aided design (CAD). According to the modeling method described in Chapter 2, the lumped mass model of this ball screw feed system test bench was built up. The equivalent parameters of the lumped mass model were calculated by using the specifications in Table 1, and the other calculated lumped parameters are listed in Table 2.

Parameter of the component | Value | Parameter of the component | Value |
---|---|---|---|

Work table mass | 206 | Rotary inertia of coupling | |

Base mass | 3820 | Rotary inertia of motor | |

Coupling mass | 1.18 | Torsional rigidity of coupling | |

Motor rotor mass | 10.9 | Screw bearing rigidity | |

Screw pitch length | Nut reference rigidity | ||

Screw diameter | Nut basic dynamic load | 37.4 | |

Screw length | 1 | Ball screw length at table position | 0.35 |

Parameter of the component | Value | Parameter of the component | Value |
---|---|---|---|

Screw equivalent mass | 11.28 | Equivalent rotary inertia of screw | |

Axial rigidity of screw | Rotary rigidity of screw | ||

Axial rigidity of base | Contact rigidity of the screw nut |

Taking the servomotor torque as input and the axial acceleration of work table as output, the frequency response characteristics of the lumped parameter model of the test bench are analyzed. The bode diagram is shown in Figure 13, and simulation result shows that the work table has four-order natural frequencies, which are 26.2, 76.7, 247, and 633 Hz. Further study shows that 76.7 Hz is the main axial vibration frequency of the work table, 26.2 Hz is the main axial vibration frequency of the base, and 247 and 633 Hz are the rotational vibration frequencies.

To establish the simulation model of servo control system, the servomotor parameters are needed as shown in Table 3.

Parameter name | Value | Parameter name | Value |
---|---|---|---|

Rated power | 4.4 | Number of pole pairs | 4 |

Rated torque | 18.6 | Stator resistance per phase | 1.44 |

Rotor inertia | Inductance | ||

Rated speed | 1500 | Permanent magnetic flux | 0.21 |

In order to compare and verify the simulation results with the experimental results, the motion command parameters and the control parameters of the experimental test and simulation are set in Table 4.

Parameter name | Value | Parameter name | Value |
---|---|---|---|

Position instruction | 400 | Position loop gain | 50 |

Maximum velocity | 400 | Velocity loop gain | 10 |

Maximum acceleration | 2000 | Current loop gain | 30 |

Maximum jerk | 20,000 |

Using the same operation parameters as set in the simulation model, a feed motion experiment was conducted and the work table position was measured. The simulation results are compared to the experimental results. Figures 14 and 15 exemplarily show simulated and measured reference velocity and feedback velocity of the servomotor at the given operating conditions. The simulation result has a similar curve to the experimental result.

Figure 16 shows frequency contents of the work table acceleration signals from simulation result and the experimental result. Comparing the simulation result with the experimental result, the co-simulation model of ball screw feed drive system can predict the vibration that occurs in the feed operation. Both results show that in this case the second-order natural frequency (about 75 Hz) but not the first-order natural frequency is the main factor influencing the performance of feed drive system.

## 5. Conclusions

In this chapter an electromechanical co-simulation model of the ball screw feed drive system was constructed based on lumped mass model of ball screw feed system and the servo control system simulation model, which can be used to study the dynamic characteristics and vibration behavior of the feed drive system. Simulative and experimental tests were conducted based on a ball screw feed drive system test bench. The result shows that the co-simulation model of ball screw feed drive system can predict the vibration that occurs in the feed operation. Because of the integration of lumped parameter model into the detailed modeled cascade control simulation model, the electromechanical co-simulation of ball screw feed drive system could achieve a very good predictability for control performance and vibration behavior study of ball screw feed drive system, which may be affected by the servo controller, ball screw feed system, or the coupling between them.

## Acknowledgments

This work is supported by the major national science and technology projects “high-end CNC machine tools and basic manufacturing equipment” (2012ZX04005031).