- About this Journal ·
- Abstracting and Indexing ·
- Advance Access ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents

Advances in Mechanical Engineering

Volume 2013 (2013), Article ID 359757, 29 pages

http://dx.doi.org/10.1155/2013/359757

## Research on Dynamic Models and Performances of Shield Tunnel Boring Machine Cutterhead Driving System

^{1}Department of Information Service and Intelligent Control, Shenyang Institute of Automation, Chinese Academy of Sciences, Shenyang 110016, China^{2}University of Chinese Academy of Sciences (Graduate School of Chinese Academy of Sciences), Beijing 100039, China^{3}NHI Shenyang Heavy Machinery Group TBM Company, Shenyang 110025, China

Received 13 June 2012; Revised 2 November 2012; Accepted 2 November 2012

Academic Editor: Marco Ceccarelli

Copyright © 2013 Xianhong Li 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.

#### Abstract

A general nonlinear time-varying (NLTV) dynamic model and linear time-varying (LTV) dynamic model are presented for shield tunnel boring machine (TBM) cutterhead driving system, respectively. Different gear backlashes and mesh damped and transmission errors are considered in the NLTV dynamic model. The corresponding multiple-input and multiple-output (MIMO) state space models are also presented. Through analyzing the linear dynamic model, the optimal reducer ratio (ORR) and optimal transmission ratio (OTR) are obtained for the shield TBM cutterhead driving system, respectively. The NLTV and LTV dynamic models are numerically simulated, and the effects of physical parameters under various conditions of NLTV dynamic model are analyzed. Physical parameters such as the load torque, gear backlash and transmission error, gear mesh stiffness and damped, pinions inertia and damped, large gear inertia and damped, and motor rotor inertia and damped are investigated in detail to analyze their effects on dynamic response and performances of the shield TBM cutterhead driving system. Some preliminary approaches are proposed to improve dynamic performances of the cutterhead driving system, and dynamic models will provide a foundation for shield TBM cutterhead driving system's cutterhead fault diagnosis, motion control, and torque synchronous control.

#### 1. Introduction

The shield TBM is the large-scale underground equipment which is used to excavate tunnels and subways. The shield TBM has the advantages of high safety and reliability, low manpower, minor environmental damages, and rapid excavation speed compared with other excavation methods. Therefore, the shield TBM is widely used in underground tunnel and subway projects. The shield TBM is mainly constituted by three core systems which include the cutterhead driving system, hydraulic thrust system, and attitude positioning system. The cutterhead driving system, hydraulic thrust system, and attitude positioning system play different important roles in shield TBM, and they are shown in Figures 1, 2, and 3, respectively. The cutterhead driving system will make the shield TBM cutterhead rotate at certain speed. The hydraulic thrust system will thrust the shield TBM cutterhead along the tunnel axis at certain speed. The attitude positioning system will make the shield TBM cutterhead follow the tunnel excavation direction. The cutterhead usually has two driven modes. One is the hydraulic-driven model and the other is motor-driven mode. This paper focuses on studying motor-driven mode shield TBM cutterhead driving system.

The shield TBM cutterhead cuts rocks and soil. Different shield TBM cutterhead forms can adapt to different geology conditions. Figure 4 shows different shield TBM cutterhead forms which mainly include the spokes form (Figure 4(a), [1, 2]), panel form (Figures 4(b)-4(c)), and hybrid form (Figure 4(d), [1, 2]). Figure 5 shows shield TBM cutterhead profile and inner mechanical transmission structure. Figures 5(a)–5(c) [1, 2] and Figure 5(d) [1, 2] are the cutterhead profile and inner mechanical transmission structure, respectively. The transmission structure in Figure 5(d) shows that multiple pinions mesh a central large gear, each pinion is driven by an induction motor, and this structure is called the multiple gears transmission structure. The multiple gears transmission structure is employed to synthesize the driving motors’ electrical magnetic torque (EMT). The motor-driven shield TBM cutterhead driving system’s physical components include driving motors, couplings, reducers, bearings, main bearing and gearbox, and frequency inverters. Figure 6 [1] shows the induction motor, reducer, and coupling of the cutterhead driving system. The reducer is the torque converter which slows down the motor speed, thus the reducer’s output torque is amplified. The overall schematic of the cutterhead driving system is shown in Figure 7. The cutterhead and central large gear have identical central shaft, thus the large gear and cutterhead have the same rotation speed. Once the central large gear rotates, then the cutterhead is driven immediately.

At present, many publications [3–16] usually introduce or study the cutterhead load torque, hydraulic thrust system, cutterhead driving structure and composition, and shield TBM overall composition. In [12], the hydraulic-driven cutterhead driving system is studied for hydraulic-driven mode shield TBM. In [13–16], shield TBM cutterhead’s load torque under various geology conditions is investigated in detail. In [17, 18], the cutterhead hydraulic thrust system is studied for hydraulic-driven mode shield TBM. In [1], a nonlinear time-varying dynamic model of the shield TBM cutterhead driving system is established and simulated. The effects of parts of physical parameters on some system performance items of the shield TBM cutterhead driving system are analyzed. In [1], the gear mesh damped is not considered. In addition, the nonlinear dynamic model in [1] considers that the shield TBM cutterhead driving system has identical gear backlash and transmission error. In [2], a linear time-varying model of shield TBM cutterhead driving system is established and analyzed, and the multiple gears transmission structure is studied by automatic dynamic analysis of mechanical systems (ADAMS) software. However, the linear dynamic model in [2] is not directly simulated.

This paper mainly focuses on studying dynamic models and system performances of the shield TBM cutterhead driving system. In this paper, different gear mesh damped, gear backlashes, and transmission errors are considered in the shield TBM cutterhead driving system, thus the nonlinear and linear dynamic models are presented again. The general NLTV dynamic model of this paper considers different mesh damped, gear backlashes, and transmission errors compared with previous dynamic models in [1, 2]. Furthermore, all the effects of physical parameters of the NLTV dynamic model on dynamic response and system performances of shield TBM cutterhead driving system are analyzed in detail. In addition, the nonlinear dynamic model and linear dynamic model are simulated under different conditions, and some preliminary approaches are proposed to improve dynamic performances of the shield TBM cutterhead driving system. Cutterhead driving system’s potential problems may be recognized by studying dynamic model. The dynamic models provide theoretical foundation for design, analysis, and control of shield TBM cutterhead driving system in the future. This paper is arranged as follows. In Section 2, the preliminary and preparation is presented. In Section 3, a general nonlinear dynamic model and linear time-varying dynamic model of shield TBM cutterhead driving system are presented, and the linear dynamic models are analyzed. In Section 4, dynamic models are simulated and analyzed under conditions of different physical parameters, and the simulation results and system performances are investigated in detail. In Section 5, study contents and results are briefly reviewed and concluded, and some conclusions are made.

#### 2. General NLTV Dynamic Model of Cutterhead Driving System

##### 2.1. Preliminary and Preparation

Gear backlash and transmission error are inevitable due to manufacturing and installing process variations, mesh requirement differences, and other uncertainties. Therefore, gear mesh process is a nonlinear transmission process. Moreover, transmission error naturally appears in the driving process and usually varies with time. As mesh point diversifies, the actual backlash varies with time. Such gear backlash variation is not considered in this paper. To establish the nonlinear dynamic model, gear mesh process’s backlash of cutterhead driving system is simplified and shown in Figures 8(a)–8(b) [1]. The gear nonlinear mesh model and linear mesh model are shown in Figures 8(c) and 8(d), respectively. Figure 8(c) shows that gear backlash and transmission error are considered in the cutterhead driving system. When gear backlash and transmission error are not considered, the gear mesh process can be seen as ideal transmission. The spring and damper are usually used to depict the gear mesh process. In fact, the mesh stiffness and mesh damped are the time-varying system parameters.

Based on gear mesh dynamic in [19–24], the relative position function, elastic mesh force, and mesh torque are obtained, respectively. Where and are the transmission error and total backlash of the th gear pair, respectively, Then, the elastic mesh force and mesh torque could be expressed as Therefore, the elastic mesh force and torque of th pinion can be equivalently expressed as The system physical parameters are defined as where is the relative backlash and is a dimensionless parameter, and are elastic mesh force coefficients, and are the mesh torque coefficients, and is the gear transmission ratio. The mesh stiffness , mesh damped , and transmission error are time-varying parameters which can be expressed as the Fourier series form where is the mesh frequency, , , and are the mean mesh stiffness, mean mesh damped, and mean transmission error, respectively, and , , and are the phase angle of mesh stiffness, mesh damped, and transmission error, respectively. Nonlinear parameters are introduced to simplify elastic mesh force and torque, thus the nonlinear parameters and are defined as Finally, the elastic mesh force and torque can be written as the following general form:

##### 2.2. Nonlinear Dynamic Model

To establish a general nonlinear dynamic model and not to lose generality, it assumes that the cutterhead is driven by motors. If the coupling mass is far less than the motor rotor mass, the coupling inertia can be ignored. In fact, the motor rotor mass cannot be larger than the coupling mass significantly, so the coupling inertia is not ignored. According to torque balance principle, the torque balance equations of the first induction motor yield For the second induction motor, the corresponding torque balance equations are obtained: Likewise, the torque balance equations of the th induction motor will yield where is the output torque of the th motor and is the input torque of the th reducer . The equations (9), (11), and (13) are substituted into (8), (10), and (12), respectively, and then the following equations will be obtained: The system parameters are defined as . Then (14) can be rewritten as the general form The input and output relationships of the th reducer can be described by the equation where is the output torque of the th reducer. Torque balance equation of the reduce- or pinion- will yield where is the input torque of pinion- and is the elastic mesh torque of pinion- or gear pair-. Substituting (18) into (17), then the torque balance equation can be written as The system parameters are defined as . Then, (19) can be rewritten as the following general form: The elastic mesh torque of the pinion- or gear pair- has been presented in Section 2.1, and the mesh torque of pinion- will be obtained as follows: The torque balance equations of the large gear can be obtained as where is total load torque of the shield TBM cutterhead. Many factors such as the cutterhead’s design parameters, physical structure, and geology conditions directly affect the load torque. The cutterhead design parameters and physical structure include the cutting depth, cutterhead opening ratio, excavation diameter, cutter tools distribution, and cutterhead forms. Geology conditions include soil properties, rocks properties, and boulders. The actual load torque will change with the geology conditions. The shield TBM cutterhead’s load torque is estimated and calculated in [13–16]. The load torque consists of the following six parts which are presented in [15]: where is the resistant torque of soil and rocks, is the friction torque of soil and rocks which chafes with the front of the cutterhead, is the friction torque of soil and rocks which chafes with the back of the cutterhead, is the friction torque of soil and rocks which chafes with bulkhead, is the stirring torque of soil and rocks by cutterhead stirring rod, and is the friction torque of cutterhead bearings and sealed chamber. Detailed derivations and descriptions of the cutterhead load torque and its parameters can be found in [15]. Therefore, calculation formulas of the load torque are directly presented here as Equations (15)–(23) constitute the NLTV dynamic model of shield TBM cutterhead driving system. Therefore, the general NLTV dynamic model is obtained as The NLTV dynamic model (25) can be combined and simplified, and then the nonlinear dynamic model is expressed as The nonlinear dynamic model can be transformed into the state space dynamic model via selecting suitable state variables, output variables, and control variables. Where is a state vector, is an output vector, and is a control vector, The MIMO state space model of the nonlinear time-varying dynamic model is obtained as where is the state matrix, is a control matrix, and is an output matrix. The matrices , , and are presented in Appendix A.

#### 3. General LTV Dynamic Model of Cutterhead Driving System

To establish general linear dynamic model, nonlinear factors in cutterhead driving system are ignored. Based on gear mesh dynamic in [19–24], the relative position, elastic mesh force, and mesh torque of linear dynamic model will be obtained, respectively, as follows: Therefore, from the previous sections, the torque balance equations of th induction motor yield where system parameters and are and , respectively. The input and output relationships of the th reducer can be described by the equation Torque balance equation of the reduce- or pinion- will yield where system parameters and are and , respectively. The mesh torque of pinion-will be obtained as follows: The torque balance equations of the large gear can be obtained as Equations (30)–(34) constitute the LTV dynamic model of the shield TBM cutterhead driving system. Therefore, the general LTV dynamic model is obtained as The linear dynamic model (35) can be combined and simplified. Therefore, the linear dynamic model is expressed as The linear dynamic model can be transformed into the state space model via selecting suitable state variables, output variables, and control variables. Where is a state vector, is an output vector, and is a control vector, The MIMO state space model of the linear time-varying dynamic model is obtained as where is the state matrix, is a control matrix, and is an output matrix. The matrices , , and are presented in Appendix B.

##### 3.1. Analysis of General LTV Dynamic Model

The general linear dynamic model of cutterhead driving system is established in the previous section. In order to analyze linear dynamic model in theory, some time-varying system parameters such as mesh stiffness, mesh damped, and load torque are replaced by mean mesh stiffness, mean mesh damped, and mean load torque. Therefore, the linear dynamic model seemed to be the linear time-invariant system. Then, the Laplace transform of (36) will be obtained as follows:

For the symbol simplicity, the system parameters are defined as . Then, (39) can be equivalently expressed as Substituting (41) into (40), then (40) can be written as It has the relationship , and then (42) can be further equivalently transformed into the following equations: Therefore, the cutterhead speed transfer function is finally obtained. The cutterhead speed transfer function shows that the dynamic performances are affected by the number of driving motors and physical parameters. To further simplify transfer function (45), it assumes that all the components of shield TBM cutterhead driving system have identical physical parameters and properties. It is stated as follows: Therefore, the cutterhead speed equation (45) could be equivalently expressed as the following: where is the denominator (characteristic polynomial) of the cutterhead speed equation (47) and and are the torque numerator and load torque numerator, respectively, Thus, the speed torque transfer function and speed load torque transfer function for the shield TBM cutterhead driving system are acquired. Where is the total torque of all the driving motors, is the cutterhead speed torque transfer function of driving motor-, is the cutterhead speed total torque transfer function, and is thef cutterhead speed load torque transfer function. The transfer function equation (49) shows that the is the poles polynomial and and are zeros polynomial. The Routh-Hurwitz stability of poles polynomial is determined by the following equation: is the set of stable Hurwitz polynomial which has negative real part. All physical parameters of the cutterhead driving system are the positive real; therefore, it has relationships , , , and . It needs to check the relationships , and then the following equation will be obtained: The item holds positive real at any time. Thus, the pole polynomial is a Hurwitz polynomial, and its characteristic roots have negative real part. The pole polynomial is a third-order system , when all the physical components of shield TBM cutterhead driving system have identical system parameters. To analyze the speed gain of the motor torque, speed torque gain, speed total torque gain, and speed load torque gain are obtained, respectively, where is the speed torque gain, is the speed total torque gain, and is the speed load torque gain. The speed gain equations show that gain , , and will be gradually reduced by increasing the number of driving motors (pinions), large gear viscous damped, motor rotor viscous damped, and pinion viscous damped. The limit of the gain , , and equals zero.

If all the driving motors of shield TBM cutterhead have identical electrical magnetic torque , then the transfer function of the cutterhead driving system will be obtained. Where , , and are speed torque transfer function, speed total torque transfer function, and speed load torque transfer function, respectively, when all the driving motors have identical electrical magnetic torque, Therefore, the speed torque gain, speed total torque gain, and speed load torque gain are obtained, respectively, when all the driving motors of shield TBM cutterhead have identical torque. Then, the speed torque gain, speed total torque gain, and speed load torque gain will have the relationships where ,, and are the speed torque gain, speed total torque gain, and speed load torque gain, respectively, when all the driving motors have identical electrical magnetic torque. Therefore, the transfer function and gain equations (53) and (54) under case show that speed torque gain and cutterhead speed torque transfer function are the times of the previous ones. The speed total torque transfer function and cutterhead speed load torque transfer function keep the same as the previous ones. The speed gain equation (54) show that speed gain , , and will be gradually reduced by increasing large gear viscous damped, motor rotor viscous damped, and pinion viscous damped. The gain and will be gradually reduced by increasing the number of driving motors (pinions); however, speed torque gain will be gradually increased by increasing the number of driving motors (pinions). The limit of the gain and is equal to zero; however, the limit of speed torque gain is equal to . When the number of the driving motors (pinions) is fixed, it can select suitable reducer ratio and transmission ratio to obtain the maximum speed gain. The maximum gain of speed torque gain and speed total torque gain will be obtained by solving the following optimization problem; therefore, the corresponding optimal transmission ratio and optimal reducer ratio will be obtained as follows: Finally, the optimal transmission ratio and optimal reducer ratio are obtained. Through analyzing the linear dynamic model, some conclusions could be made for shield TBM cutterhead driving system’s linear dynamic model: (i) dynamic performances of the shield TBM cutterhead driving system are affected by physical parameters such as the number of driving motors (pinions), transmission ratio, reducer ratio, large gear inertia and viscous damped, motor rotor inertia and viscous damped, pinion inertia and viscous damped, mesh stiffness, and mesh damped. (ii) System order of the cutterhead speed or cutterhead driving system is when system parameters are not equal to each other. The shield TBM cutterhead driving system degenerates into a third-order system when system parameters are equal to each other, and the cutterhead driving system can be seen as a single input system at this time. (iii) The cutterhead driving system is an open-loop stable system with negative zeros, the speed torque gain , speed total torque gain , and speed load torque gain will be gradually reduced by increasing the number of driving motors (pinions), large gear viscous damped, motor rotor viscous damped, and pinion viscous damped. In addition, the limit of the gain , , and is equal to zero. (iv) When all the driving motors’ electrical magnetic torque has relationship with , speed torque gain and cutterhead speed torque transfer function are the times of the previous ones. Speed torque gain will be gradually increased by increasing the number of driving motors (pinions); however, the gain and will be gradually reduced. The limit of the gain and is equal to zero; however, the limit of speed torque gain is equal to . (v) It can select suitable reducer ratio or transmission ratio to heighten the speed torque gain and speed total torque gain . The corresponding optimal reducer ratio and optimal transmission ratio are obtained by maximizing the gain or .

#### 4. Simulation Results and Analysis of Shield TBM Cutterhead Driving System

There are not good theoretical analysis methods for NLTV dynamic model and LTV dynamic model. To facilitate the simulation and the analysis of dynamic models of the shield TBM cutterhead driving system, some time-varying system parameters such as mesh stiffness, mesh damped, load torque, and transmission error are replaced by the mean mesh stiffness, mean mesh damped, mean load torque, and mean transmission error in the simulation. Time-varying physical parameters mean that physical parameters will change and be affected by many uncertain factors, and these uncertain factors may change with time. Therefore, from the perspective of the timeline, these physical parameters could be seen as time varying. The general NLTV and LTV dynamic models of shield TBM cutterhead driving system are simulated and studied. Effects of the NLTV dynamic model’s physical parameters on the dynamic response and performances of shield TBM cutterhead driving system are investigated in detail under different conditions. The differences of dynamic response and performances between NLTV model and LTV model are discussed under different load torque conditions. The simulation parameters are presented in Appendix C. When discussing a specific physical parameter’s effect on the dynamic response and performances, other physical parameters are all fixed in the simulation. The corresponding simulation results are shown in Figures 9–21.

##### 4.1. Physical Parameters’ Effects on Dynamic Response and Performances

The effects of all the physical parameters such as load torque, mesh stiffness and mesh damped, motor rotor inertia and viscous damped, pinion inertia and viscous damped, gear backlash, transmission error, and large gear inertia and viscous damped on the dynamic response and performances of the shield TBM cutterhead driving system are analyzed. The simulation results of the NLTV dynamic model are shown and investigated under different conditions.

The load torque of shield TBM cutterhead driving system is determined by outside geology conditions and cutterhead design parameters. The response results of the NLTV dynamic model in various load torque cases are shown in Figures 9(a)–9(f). The results reveal that the cutterhead driving system shows good dynamic response and performances when load torque is increased from 0 KN·m (unload case) to 500 KN·m and 1000 KN·m. The cutterhead steady-state speed is about 17.21 rpm, 14.32 rpm, and 11.33 rpm, respectively, when the load torque is 0 KN·m, 500 KN·m, and 1000 KN·m, respectively. The response results show that dynamic response and performances are deteriorated when the load torque is increased to 2500 KN·m. The cutterhead speed shows severe fluctuations and also shows a lot of nonlinear characteristics which include the irregular jump and oscillation, jump discontinuity, and harmonic phenomenon. The cutterhead speed also shows few nonlinear characteristics in unload case, but nonlinear characteristics do not obviously appear when the load torque is 500 KN·m and 1000 KN·m. The response results in Figures 9(a)–9(f) also show that the oscillation amplitude of cutterhead initial speed will be increased and the steady-state speed amplitude will be reduced when increasing the load torque. Therefore, load torque has large effects on dynamic response and performances of the shield TBM cutterhead driving system.

The gear backlash and transmission error cannot be avoided in the shield TBM cutterhead driving system. The response results of the NLTV dynamic model under various gear backlashs and transmission errors (in unload case) are shown in Figures 10(a)–10(d) and Figures 11(a)–11(f), respectively. The dynamic response results in Figures 10(a)–10(d) show that the cutterhead speed shows peak pulse and jump discontinuity behaviors when gear relative backlash is increased. The response results in Figures 10(a)–10(d) obviously reveal that cutterhead speed shows numerous nonlinear characteristics and behaviors such as irregular jump and jump discontinuity when the gear relative backlash is increased from 0.001 to 0.01, 0.05, and 0.1. However, nonlinear characteristics and behaviors are dramatically reduced when the gear relative backlash is increased from 0.01 to 0.05. Numerous nonlinear behaviors appear again when relative backlash is increased from 0.05 to 0.1. The maximum jump amplitude and minimum jump amplitude of the cutterhead speed are about 26.88 rpm and 5.09 rpm, 16.38 rpm and 7.38 rpm, and 25.56 rpm and 6.13 rpm when the relative gear backlash is 0.01, 0.05, and 0.1, respectively. So, gear relative backlash has complex and important effects on dynamic response and performances of the shield TBM cutterhead driving system.

The response results in Figures 11(a)–11(f) obviously reveal that the transmission error has complex and important effects on the dynamic response and performances of the cutterhead driving system. The results show that cutterhead speed shows numerous nonlinear behaviors such as irregular jump and oscillation, jump peak, and jump discontinuity when the transmission error is increased from 0.001 m to 0.004 m, 0.006 m, and 0.008 m. However, nonlinear behaviors and characteristics of the cutterhead speed are reduced when transmission error is increased from 0.006 m to 0.008 m. The steady-state speed ripples and fluctuations will be increased when the transmission error is increased. The dynamic response results in Figures 11(a)–11(f) also show that oscillation amplitude of cutterhead initial speed will be increased when transmission error is heightened. The oscillation amplitude of cutterhead initial speed is about −3.24 rpm, −12.20 rpm, −18.69 rpm, and −26.95 rpm, respectively, when transmission error is 0.001 m, 0.004 m, 0.006 m, and 0.008 m, respectively. Thus, transmission error has important effects on dynamic response and performances of the shield TBM cutterhead driving system.

The nonlinear behaviors and characteristics can be directly affected by gear relative backlash and transmission error. Hence, nonlinear behaviors would reflect actual working conditions of the cutterhead driving system in some degree. If relative backlash and transmission error exceed a certain level, they will cause noise, speed ripple, and speed discontinuity jump which will cause damage to cutterhead driving system in the long run. The response results indicate that the gear backlash and transmission error will tend to strengthen nonlinearity and worsen dynamic performances. The cutterhead driving system will reap benefits for reducing gear relative backlash and transmission error.

The response results of the NLTV dynamic model under various conditions of inertia and viscous damped of large gear (in unload case) are shown in Figures 12(a)–12(d) and Figures 13(a)–13(d), respectively. The response results in Figures 12(a)–12(d) show that large gear inertia has complex effects on the dynamic response and performances of the cutterhead driving system. The cutterhead speed shows few nonlinear behaviors as the large gear inertia is 56.93 kg m^{2}. The cutterhead speed also shows lots of nonlinear behaviors such as irregular jump, irregular oscillation, and jump discontinuity when the large gear inertia is increased form 56.93 kg m^{2} to 86.93 kg m^{2}, 126.93 kg m^{2}, and 156.93 kg m^{2}. The cutterhead speed shows numerous nonlinear behaviors such as irregular jump and jump discontinuity when the large gear inertia is increased from 86.93 kg m^{2} to 126.93 kg m^{2} and 156.93 kg m^{2}. Therefore, large gear inertia has important effects on dynamic response and performances of the shield TBM cutterhead driving system.

The response results in Figures 13(a)–13(d) obviously reveal that large gear viscous damped has great and complex effects on the dynamic response and performances of the cutterhead driving system. The cutterhead speed does not obviously show nonlinear behaviors and characteristics such as the irregular jump and jump discontinuity phenomenon when viscous damped is increased from 0.921 kg m^{2}/rad s^{−1} to 1.921 kg m^{2}/rad s^{−1}. However, the cutterhead speed does show numerous irregular jump and jump discontinuity behaviors when large gear viscous damped is increased from 1.921 kg m^{2}/rad s^{−1} to 4.921 kg m^{2}/rad s^{−1}, and the jump discontinuity behaviors keep a long time. But the nonlinear behaviors obviously disappear when the large gear viscous damped is increased from the 4.921 kg m^{2}/rad s^{−1} to 7.921 kg m^{2}/rad s^{−1}. The amplitude of cutterhead steady-state speed is not affected when the large gear viscous damped is increased. Therefore, response results indicate that large gear viscous damped has important and complex effects on dynamic response and system performances of the shield TBM cutterhead driving system.

The response results of the NLTV dynamic model under various motor rotor inertia and viscous damped conditions (in unload case) are shown in Figures 14(a)–14(d) and Figures 15(a)–15(d), respectively. The response results in Figures 14(a)–14(d) reveal that motor rotor inertia has large effects on the dynamic response and performances of the cutterhead driving system. The dynamic response results in Figures 14(a)–14(d) obviously do not show nonlinear behaviors and characteristics such as irregular jump, jump discontinuity, irregular oscillation, and harmonics when the motor rotor inertia is continuously increased from and to and . The transient performance items of cutterhead driving system such as rise time, delay time, and setting time are gradually increased when the motor rotor inertia is continuously increased. However, the amplitude of cutterhead steady-state speed would not be reduced when the motor rotor inertia is continuously increased. The dynamic response results show that transient performances have been affected significantly. Therefore, motor rotor inertia has large effects on dynamic response and performances of the shield TBM cutterhead driving system.

The response results in Figures 15(a)–15(d) reveal that motor rotor viscous damped has important and complex effects on the dynamic response and performances of the cutterhead driving system. The dynamic response results in Figures 15(a)–15(d) obviously do show numerous nonlinear behaviors and characteristics such as irregular jump, jump discontinuity, and irregular oscillation when the motor rotor viscous damped is increased from and to and . However, the dynamic response results also show that numerous nonlinear behaviors and characteristics such as irregular jump and jump discontinuity dramatically vanish when the motor rotor viscous damped is increased from and to and and and , respectively. The dynamic response results obviously demonstrate that the amplitude of cutterhead steady-state speed is continuously reduced when the motor rotor viscous damped is continuously increased. Therefore, motor rotor viscous damped has significant and complex effects on dynamic response and performances of the shield TBM cutterhead driving system.

The response results of the NLTV dynamic model under various conditions of pinion inertia and viscous damped (in unload case) are shown in Figures 16(a)–16(d) and Figures 17(a)–17(d), respectively. The response results in Figures 16(a)–16(d) reveal that pinion inertia has complex effects on the dynamic response and performances of the cutterhead driving system. The dynamic response results do not obviously show nonlinear behaviors and characteristics such as irregular jump, jump discontinuity, and irregular oscillation when the pinion inertia is increased from to . However, response results show that numerous nonlinear behaviors such as irregular jump and jump discontinuity dramatically appear when the pinion inertia is increased from 10.00 kg m^{2} to 20.00 kg m^{2}. But the nonlinear behaviors vanish again when the pinion inertia is increased from 20.00 kg m^{2} to 50.00 kg m^{2}. The response results show that the pinion inertia has complex affects on dynamic response and performances of the cutterhead driving system.

The response results in Figures 17(a)–17(d) reveal that pinion viscous damped has complex effects on the dynamic response and performances of the cutterhead driving system. The dynamic response results in Figures 17(a)–17(d) show numerous nonlinear behaviors and characteristics such as irregular jump, jump discontinuity, and irregular oscillation when the pinion viscous damped is increased from and to and . The dynamic response results show that numerous nonlinear behaviors and characteristics of the cutterhead speed obviously disappear when pinion viscous damped is increased from and to and and and , respectively. These nonlinear behaviors and characteristics dramatically vanish when the pinion viscous damped is further increased. The dynamic response results obviously demonstrate that the amplitude of cutterhead steady-state speed is not reduced when the pinion viscous damped is continuously increased. Therefore, pinion viscous damped has complex and significant effects on dynamic response and performances of the shield TBM cutterhead driving system.

The response results of the NLTV dynamic model under various mesh stiffness and mesh damped conditions (in unload case) are shown in Figures 18(a)–18(f) and Figures 19(a)–19(f), respectively. The response results in Figures 18(a)–18(f) obviously show that mesh stiffness has important and complex effects on the dynamic response and performances of the cutterhead driving system. The cutterhead speed does not obviously show nonlinear behaviors and characteristics such as the jump discontinuity and irregular jump when the gear mesh stiffness is reduced from 10^{5} N·m/rad to 10^{4} N·m/rad. However, the cutterhead speed shows lots of nonlinear characteristics and behaviors such as irregular jump, irregular oscillation, and jump discontinuity when mesh stiffness is increased form 10^{5} N·m/rad to 10^{6} N·m/rad. Furthermore, the cutterhead speed also shows numerous nonlinear behaviors such as irregular jump, jump discontinuity, irregular oscillation, unstable oscillation, and subharmonic behavior when mesh stiffness is increased from 10^{6} N·m/rad to 10^{7} N·m/rad. The response results show that cutterhead speed contains lots of subharmonic and oscillation behaviors when mesh stiffness is increased beyond a certain range. Dynamic performances tend to be worsened by increasing mesh stiffness. Meanwhile, the oscillation amplitude of cutterhead initial speed and cutterhead steady-state speed ripple are increased, and nonlinear characteristics are motivated as well when the mesh stiffness is increased. Therefore, mesh stiffness has important and complex effects on dynamic response and performances of the shield TBM cutterhead driving system.

The response results in Figures 19(a)–19(f) obviously reveal that mesh damped has important and complex effects on the dynamic response and performances of the cutterhead driving system. The cutterhead speed does not obviously show nonlinear behaviors and characteristics such as the irregular jump, irregular oscillation, and jump discontinuity when the mesh damped is reduced from 100 N·m/ to 50 N·m/. However, the cutterhead speed shows numerous nonlinear characteristics and behaviors such as irregular jump, irregular oscillation, and jump discontinuity when mesh damped is further reduced from 50 N·m/ to 10 N·m/ and 5 N·m/. In addition, cutterhead steady-state speed contains numerous subharmonic and stable oscillation behaviors when mesh damped is reduced to 10 N·m/ and 5.0 N·m/. Meanwhile, results show that the oscillation amplitude of cutterhead initial speed and cutterhead steady-state speed ripple would be increased, and nonlinear characteristics will be motivated too when mesh damped is gradually reduced. Dynamic performances of the cutterhead driving system tend to be worsened by reducing mesh damped. Therefore, mesh damped has very important and complex effects on dynamic response and performances of the shield TBM cutterhead driving system.

The results indicate that mesh stiffness is increased considerably which tends to aggravate the dynamic performances and increase nonlinearity behaviors. The results also show that the steady-state speed ripple and nonlinear phenomenon such as complex harmonic, irregular jump, and jump discontinuity behaviors will be motivated by increasing mesh stiffness or reducing mesh damped. As a result, the cutterhead steady-state speed ripple, cutterhead harmonic speed, and nonlinear behaviors could be reduced by reducing mesh stiffness or increasing mesh damped. Consequently, the mesh stiffness and mesh damped are critical parameters which should be reasonably designed for the cutterhead driving system.

The differences between the NLTV dynamic model and LTV dynamic model of cutterhead driving system are comparatively studied, and their simulation results under various load torque conditions are shown in Figures 20(a)–20(d) (in unload case) and Figures 21(a)–21(f), respectively. The response results in Figures 20(a)–20(d) show that nonlinear model and linear model nearly have the same dynamic response in unload case. They nearly have the same cutterhead speed and transient response process. The nonlinear model shows a few nonlinear behaviors such as jump discontinuity and irregular jump, whereas linear model does not have these behaviors. The speed response of the nonlinear model shows speed oscillation in initial stage. The response results of nonlinear model and linear model under various load torque conditions are shown in Figures 21(a)–21(f), where response results indicate that nonlinear model and linear model have roughly the same response results when the load torque is 500 KN·m and 1000 KN·m, respectively. Small differences between the nonlinear model and linear model are cutterhead speed response and speed oscillation amplitude in initial stage when the load torque is fixed at 500 KN·m and 1000 KN·m, respectively. However, the response results in Figures 21(a)–21(f) show that nonlinear model and linear model have large differences when load torque is 2500 KN·m. The nonlinear model shows larger cutterhead speed fluctuation and steady-state speed ripple than those of linear model. The cutterhead speed of nonlinear model shows lots of nonlinear behaviors such as irregular jump, jump discontinuity, and subharmonics. Of course, the linear model also shows some speed ripple and jump discontinuity behaviors, but linear model’s speed fluctuation and speed ripple are smaller than those of the nonlinear model.

#### 5. Conclusion

In this paper, a general NLTV dynamic model and LTV dynamic model are presented for the shield TBM cutterhead driving system. The general NLTV dynamic model considers different transmission errors, gear backlashes, and gear mesh damped compared with the previous nonlinear dynamic model. MIMO state space forms of the NLTV and LTV dynamic model are presented as well. Through theoretical analysis of linear dynamic model, the analysis results reveal the following.(i)Dynamic response and performances of the cutterhead driving system will be affected by the number of driving motors (active pinions), transmission ratio, reducer ratio, large gear inertia and viscous damped, motor rotor inertia and viscous damped, pinion inertia and viscous damped, and mesh stiffness and mesh damped.(ii)The model order of cutterhead driving system is the when physical parameters are not equal to each other. However, the cutterhead driving system degenerates into a third-order system when physical parameters are equal to each other, and it can be seen as a single input system at this time.(iii)The cutterhead driving system is a minimum phase system (open-loop stable system with negative zeros). The speed torque gain, speed total torque gain, and speed load torque gain will be gradually reduced by increasing the number of driving motors (active pinions), large gear viscous damped, motor rotor viscous damped, and pinion viscous damped. The limits of all the speed torque gain, speed total torque gain, and speed torque gain are equal to zero.(iv)When motor electrical magnetic torque has , speed torque gain and speed torque transfer function are the times of the previous ones. Speed torque gain will be gradually increased by increasing the number of driving motors (active pinions), but speed total torque gain and speed load torque gain will be gradually reduced. The limits of speed total torque gain and speed load torque gain are equal to zero, but the limit of speed torque gain is equal to .(v)It can select suitable reducer ratio or transmission ratio to heighten the speed torque gain and speed total torque gain, and corresponding optimal reducer ratio and optimal transmission ratio are derived, respectively. The optimal TR and RR could help cutterhead driving system to acquire higher rotation speed.

To numerically study and analyze the NLTV dynamic model, the NLTV dynamic model is simulated under various physical parameters conditions. Physical parameters such as the load torque, gear backlash and transmission error, gear mesh stiffness and mesh damped, pinions inertia and damped, large gear inertia and damped, and motor rotor inertia and damped are investigated in detail to analyze their effects on the dynamic response and performances of the shield TBM cutterhead driving system. The response results of NLTV dynamic model show lots of interesting nonlinear characters and behaviors such as irregular jumps, jump discontinuity, irregular oscillations, and harmonics and subharmonics especially for cutterhead driving system under conditions of unload, heavy load, large backlash and transmission error, large mesh stiffness, small mesh damped, great pinion damped, great large gear inertia and damped, and large motor rotor damped. Through numerical studies and analysis of NLTV dynamic model, the numerical simulation results reveal the following.(i)The load torque has large effects on dynamic response and performances of cutterhead driving system. By increasing load torque, it will lead to increasing the oscillation amplitude of cutterhead initial speed and reducing the amplitude of cutterhead steady-state speed. When load torque is increased considerably, the cutterhead speed will show numerous nonlinear behaviors and characteristics.(ii)Gear relative backlash and transmission error have complex and important effects on dynamic response and performances of the cutterhead driving system. The cutterhead speed will show numerous nonlinear characteristic and behaviors such as irregular jump, and jump discontinuity when the relative backlash and transmission error are increased. Increasing transmission error, the oscillation amplitude of cutterhead initial speed will be increased as well.(iii)The large gear inertia and damped have complex and great effects on dynamic response and performances of the cutterhead driving system. When large gear inertia and damped are increased to a certain extent, the cutterhead speed will show numerous nonlinear behaviors such as irregular jump, irregular oscillation, and jump discontinuity. However, the nonlinear behaviors will disappear when damped is increased beyond a certain level.(iv)The motor rotor inertia and damped have complex and large effects on dynamic response and performances of the cutterhead driving system. The transient performance items such as rise time, delay time, and setting time will be increased when increasing the motor rotor inertia. The cutterhead speed will show numerous nonlinear behaviors when motor rotor damped is increased to a certain extent. Furthermore, cutterhead steady-state speed will be reduced by increasing the motor rotor damped. However, nonlinear behaviors will disappear when the motor rotor damped is increased beyond a certain level.(v)The pinion inertia and damped have complex and large affects on dynamic response and performances of the cutterhead driving system. The cutterhead speed will show lots of nonlinear behaviors and characteristics such as the irregular jump, irregular oscillation, and jump discontinuity when the pinion damped is increased to a certain extent. However, the nonlinear behaviors will disappear when they are increased beyond a certain level.(vi)The mesh stiffness and damped have complex and great affects on dynamic response and performances of the cutterhead driving system. The cutterhead speed will show lots of nonlinear behaviors and characteristics such as irregular jump, and jump discontinuity when the mesh stiffness is increased to a certain extent. Increasing mesh stiffness, the oscillation amplitude of cutterhead speed in initial stage will be increased as well. When the mesh stiffness is increased beyond a certain level, the cutterhead speed will show severe oscillations and harmonics. The cutterhead speed will show numerous nonlinear behaviors such as irregular jump and jump discontinuity when mesh damped is reduced to a certain extent. Reducing mesh damped, the oscillation amplitude of cutterhead speed in initial stage will be increased as well.(vii)The nonlinear model and linear model nearly have the same dynamic response in unload case and certain load torque case. The nonlinear model shows a few nonlinear behaviors such as jump discontinuity and irregular jump, whereas linear model does not show these behaviors. The speed response of the nonlinear model shows speed oscillation in initial stage. However, the nonlinear model and linear model show large differences of dynamic response when load torque is beyond a certain level.Through studying physical parameters’ effects on dynamic response and performances, several approaches are proposed to improve system performances. The approaches include the following. It can design reasonable reducer ratio and transmission ratio, pinion inertia and damped, motor rotor inertia and damped, large gear inertia and damped, and mesh stiffness and damped and also reduce gear backlash and transmission error. The NLTV and LTV dynamic models provide the basis for control and fault diagnosis of shield TBM cutterhead driving system.

#### Appendices

#### A. Matrices A, B, and C for NLTV Dynamic Model

We have the following: