**World Journal of Engineering and Technology**

Vol.07 No.02(2019), Article ID:92373,9 pages

10.4236/wjet.2019.72B001

Evaluation Method for Tunneling Stability of TBM Cutterhead

Zhaohui Xu^{1*}, Ye Er-ken∙Zha Mu-ti^{2}, Lin Xue^{1}, Fengtao Wang^{1}, Jing Chen^{3}, Junzhou Huo^{1}^{ }

^{1}School of Mechanical Engineering, Dalian University of Technology, Dalian, China

^{2}Ergis River Basin Development Project Construction Management Bureau, Urumchi, China

^{3}School of Navigation and Shipbuilding Engineering, Dalian Ocean University, Dalian, China

Received: January 14, 2019; Accepted: May 14, 2019; Published: May 17, 2019

ABSTRACT

In the process of tunneling of tunnel boring machine (TBM), different geological conditions often correspond to different working conditions, and the randomness of geological conditions also causes the order of occurrence of each working condition to be different. Under the conversion of different working conditions, this makes the vibration of different types of cutterheads different. How to choose the appropriate type of cutterhead according to different geological conditions is very important for saving engineering cost and increasing cutterhead life. In view of the above situation, this paper proposes a stability evaluation method during the TBM tunneling process to select the appropriate cutterhead type. Firstly, the corresponding relationship between geology and working conditions is established according to different geological conditions, and the input loads corresponding to geological conditions are obtained. Then, it is substituted into the dynamic model of the cutterhead system, the vibration response boundaries of each degree of freedom are obtained by solving. And the average value of the maximum boundary amplitude of each degree of freedom is taken to represent the extreme vibration of the cutterhead under the corresponding working conditions. Finally, by comparing the fluctuation of the ultimate vibration amplitude of each type of cutterhead in the process of working condition conversion, the results are as follows: when the transition between homogeneous strata and composite strata is normal and there is no large turning and deviation correction, the vibration response of the two-part cutterhead is the smallest, and the two-part cutterhead is the best choice. Otherwise, the five-part cutterhead is the best choice, while the stability of the integrated cutterhead is the worst.

**Keywords:**

Cutterhead, Working Condition Conversion, Stability, Type Selection of Cutterhead

1. Introduction

Mechanised tunneling is an efficient tunneling technology for the construction fnew underground infrastructures, in particular in urban environments [1]. Tunnel boring machine (TBM) is a factory-built tunnel construction system integrating mechanical, electrical, hydraulic and other multidisciplinary applications. It is widely used in underground projects such as water conservancy, transportation, national defense and energy. Cutterhead system is the most critical component, its engineering failure will greatly affect the progress of the project, increase time and money costs. Cutterhead system tunneling process is often different because of different geological conditions. How to carry out dynamic analysis more effectively and accurately is of great significance for studying the tunneling stability of the cutterhead system and selecting the type of cutterhead.

At present, scholars have done a lot of research on the dynamic characteristics of the key parts of TBM. Zhou et al. [2] proposed and established an equivalent coupling model through the secondary development of the finite element software ABAQUS, and carried out the corresponding dynamic analysis. Huoet al. [3] established a multi-degree-of-freedom coupling dynamic model of TBM, and studied the dynamic response of cutterhead under impact load. Li et al. [4] established a general nonlinear time-varying (NLTV) dynamic model and a linear time-varying (LTV) dynamic model for the shield tunneling machine drive system, and the effects of physical parameters under various conditions of dynamic model are analyzed. Zhang, Li et al. [5] [6] established a dynamic model of shield machine under complex geological conditions, calculated the dynamic response of the cutterhead slewing system, and explored the influence of key parameters. Cooley and Parker [7] studied the effect of the speed of the planetary carrier on the natural frequency of the system. Han [8] introduced a nonlinear spring element to simulate the interaction between the cutterhead and the rock, and established a finite element model to simulate the dynamic process of the tunneling process.

In conclusion, although scholars have done a lot of research on the key parts of TBM, most of the dynamic analysis of TBM has not considered the influence of geological conditions, and few studies have been carried out on geological transformation and cutterhead selection. In view of the above problems, this paper presents a stability evaluation method in TBM tunneling process. Firstly, the corresponding relationship between geology and working conditions is established according to different geological conditions, and the input loads corresponding to geological conditions are obtained. Then, it is substituted into the dynamic model of the cutterhead system, the vibration response boundaries of each degree of freedom are obtained by solving. And the average value of the maximum boundary amplitude of each degree of freedom is taken to represent the extreme vibration of the cutterhead under the corresponding working conditions. Finally, by comparing the fluctuation of the ultimate vibration amplitude of each type of cutterhead in the process of working condition conversion, the stability of each type of cutterhead under different geological conditions is obtained.

2. Transformation of Geological and Working Conditions

According to the statistics of engineering geology, the location of TBM contains many strata in the process of tunneling, and the distribution characteristics of strata are uncertain. According to the research of relevant scholars, the changing trend of strata on the cutterhead tunneling interface and the distribution and combination of different geology have great influence on the fluctuation and magnitude of cutter head load. Therefore, it is necessary to deal with the stratigraphic distribution. This article uses classical working conditions to equate possible geological forms. The working conditions that may occur during the tunneling process are attributed to the maximum thrust condition, the maximum overturning condition, the turning and correcting condition and out of difficult condition. The typical working conditions are shown in Figure 1. The maximum thrust condition refers to the equivalent condition of TBM driving along a straight line under uniform geology. At this time, all the cutters are subjected to rock breaking load, and the cutterhead system bears the maximum load, and its force section is shown in Figure 1(a); Maximum overturning condition refers to the equivalent condition of TBM driving in composite strata. At this time, half of the cutter of the cutterhead contacts the rock, and the cutterhead system bears large loads and bending moments, and its force section is shown in Figure 1(b); The turning and correcting condition means that there is a certain angle between the cutterhead surface and the tunneling surface, only part of the cutter is subjected to force, and its force section is shown in Figure 1(c). When TBM is driving, it may encounter problems such as the lone stone or the edge cutter being stuck. At this time, the cutterhead torque must be increased. The cutter head is subjected to the maximum torque under this condition, and its force section is shown in Figure 1(d).

Considering the corresponding relationship between the geological distribution and the working conditions, the calculation process of the cutterhead load boundary is as follows.

1) Combining simulation and field test data to obtain the load waveform of a single cutter;

2) According to the geological survey report, the geological conditions are statistically analyzed and classified, and several typical working conditions with the largest proportion are found;

3) According to the obtained working conditions and the load waveform of the single cutter, the time history of the equivalent load of the cutterhead under working conditions with the largest proportion is fitted, and the minimum value of the cutterhead load is obtained.

4) The maximum value of the cutterhead load is generated when the load on each cutter is nominal. Therefore, the load interval of the cutterhead under

Figure 1. Transformation of geology and working conditions.

several typical working conditions with the largest proportion can be obtained and used as the load input of the dynamic model.

3. Establishment of Dynamic Model

3.1. Equivalent of Mechanical Model

The TBM cutterhead system includes cutterhead, ring gear, pinion, shield, reducer and motor. The typical cutterhead forms are integrated cutterhead (IC), two-part cutterhead (TPC), and five-part cutterhead (FPC). The cutterhead, ring gear, pinion, shield and motor are equivalent to mass points, and the pinion and red tarder are equivalent to a mass point. The equivalent mechanical model of cutterhead is established, as shown in Figure 2.

m_{δ}(δ = L, r, d, pj) represents the equivalent mass of cutterhead, ring gear, shield and pinion, respectively. k_{eqδ}(δ = x, y, L, r, z, d, dz, p) represents horizontal equivalent support stiffness of cutterhead, longitudinal equivalent support stiffness of cutterhead, radial equivalent support stiffness of ring gear, axial equivalent support stiffness of ring gear, horizontal equivalent support stiffness longitudinal equivalent support stiffness and axial equivalent support stiffness of shield, and the equivalent support stiffness of the pinion respectively. C_{eqδ}(δ = x, y, L, r, z, d, dz, p) represents horizontal damping coefficient of cutterhead, longitudinal damping coefficient of cutterhead, radial damping coefficient of ring gear, axial damping coefficient of ring gear, horizontal damping coefficient, longitudinal damping coefficient and axial damping coefficient of shield, and the damping coefficient of the pinion respectively. T_{L}, T_{mj}, F_{Xi}, F_{Y}, M_{X}, M_{Y}, F_{L} represents cutterhead torque, motor torque, cutterhead horizontal load, cutterhead longitudinal load, cutterhead overturning moment, cutterhead axial load.

According to the equivalent mechanical model, the generalized displacement matrix of the cutter head system is:

$\begin{array}{l}\left\{\delta \right\}=\{{X}_{L},{Y}_{L},{Z}_{L},{\theta}_{Lx},{\theta}_{Ly},{\theta}_{L},{X}_{r},{Y}_{r},\\ \text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{Z}_{r},{\theta}_{x},{\theta}_{y},{\theta}_{r},{X}_{d},{Y}_{d},{Z}_{d},{H}_{pj},{V}_{pj},{\theta}_{pj},{\theta}_{mj}{\}}^{T}\end{array}$ (j = 1 - pinion number) (1)

In the matrix, X_{L}, Y_{L}, Z_{L}, θ_{Lx}, θ_{Ly}, θ_{L} represents the X-direction displacement of the cutterhead, the Y-direction displacement of the cutterhead, the axial

(a) (b) (c) (d)

Figure 2. Equivalent mechanical model of cutterhead. (a) Mechanics Model of Rotary System; (b) Torsional Mechanics Model; (c) Mechanical Model of Axis and Overturning; (d) Radial Mechanics Model.

displacement of the cutterhead, the overturning vibration displacement of the cutterhead, and the torsional vibration displacement of the cutterhead respectively; X_{r}, Y_{r}, Z_{r}, θ_{x}, θ_{y}, θ_{r} represents the X-direction displacement of the ring gear, the Y-direction displacement of the ring gear, the axial displacement of the ring gear, the overturning vibration displacement of the ring gear, and the torsional vibration displacement of the ring gear respectively; X_{d}, Y_{d}, Z_{d} represents the X-direction displacement of the shield, the Y-direction displacement of the shield, the axial displacement of the shield respectively; H_{pj}, V_{pj}, θ_{pj} represents tangential displacement of pinion, radial displacement of pinion, the torsional vibration displacement of pinion respectively; θ_{mj} represent angular displacement of torsional vibration of motor.

3.2. Establishment of Dynamics Differential Equation of Cutterhead System

According to the dynamic equivalent model of the cutterhead and the transformation relationship between geology and working conditions, the corresponding dynamic differential equations are obtained.

1) Cutterhead

$\{\begin{array}{l}{m}_{L}{\ddot{X}}_{L}+{C}_{eqx}({\dot{X}}_{L}-{\dot{X}}_{r})+{k}_{eqx}({X}_{L}-{X}_{r})={F}_{X}^{I}\\ {m}_{L}{\ddot{Y}}_{L}+{C}_{eqy}({\dot{Y}}_{L}-{\dot{Y}}_{r})+{k}_{eqy}({Y}_{L}-{Y}_{r})={F}_{Y}^{I}\\ {m}_{L}{\ddot{Z}}_{L}+{\displaystyle \sum _{i=1}^{4}[{C}_{eqLi}({\dot{Z}}_{Li}-{\dot{Z}}_{ri})+{k}_{eqLi}({Z}_{Li}-{Z}_{ri})}]={F}_{L}^{I}\\ {I}_{Lx}{\ddot{\theta}}_{Lx}+{r}_{br}[{C}_{eqL1}({\dot{Z}}_{L1}-{\dot{Z}}_{r1})-{C}_{eqL3}({\dot{Z}}_{L3}-{\dot{Z}}_{r3})+{k}_{eqL1}({Z}_{L1}-{Z}_{r1})\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}-{k}_{eqL3}({Z}_{L3}-{Z}_{r3})]={M}_{X}^{I}\\ {I}_{Ly}{\ddot{\theta}}_{Ly}+{r}_{br}[{C}_{eqL2}({\dot{Z}}_{L2}-{\dot{Z}}_{r2})-{C}_{eqL4}({\dot{Z}}_{L4}-{\dot{Z}}_{r4})+{k}_{eqL2}({Z}_{L2}-{Z}_{r2})\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}-{k}_{eqL4}({Z}_{L4}-{Z}_{r4})]={M}_{Y}^{I}\\ {I}_{L}{\ddot{\theta}}_{L}+{C}_{rLQ}({\dot{\theta}}_{L}-{\dot{\theta}}_{r})+{k}_{rLQ}({\theta}_{L}-{\theta}_{r})=-{T}_{L}^{I}\end{array}$ (2)

2) Ring gear

$\{\begin{array}{l}{m}_{r}{\ddot{X}}_{r}+{\displaystyle \sum _{j=1}^{n}({F}_{prj}+{D}_{prj})\mathrm{sin}({\phi}_{j}+\alpha )+{C}_{eqpr}({\dot{X}}_{r}-{\dot{X}}_{d})}+{C}_{eqx}({\dot{X}}_{r}-{\dot{X}}_{L})\\ +{k}_{eqr}({X}_{r}-{X}_{d})+{k}_{eqx}({X}_{r}-{X}_{L})=0\\ {m}_{r}{\ddot{Y}}_{r}-{\displaystyle \sum _{j=1}^{n}({F}_{p\text{\hspace{0.05em}}rj}+{D}_{p\text{\hspace{0.05em}}rj})\mathrm{cos}({\phi}_{j}+\alpha )+{C}_{eqr}({\dot{Y}}_{r}-{\dot{Y}}_{d})}+{C}_{eqy}({\dot{Y}}_{r}-{\dot{Y}}_{L})\\ +{k}_{eqr}({Y}_{r}-{Y}_{d})+{k}_{eqy}({Y}_{r}-{Y}_{L})=0\\ {m}_{r}\ddot{Z}{}_{r}+{\displaystyle \sum _{i=1}^{4}[{C}_{eqLi}({\dot{Z}}_{ri}-{\dot{Z}}_{Li})+{C}_{eqzi}({\dot{Z}}_{ri}-{\dot{Z}}_{d})+{k}_{eqLi}({Z}_{ri}-{Z}_{Li})}+{k}_{eqzi}({Z}_{ri}-{Z}_{d})]=0\\ {I}_{rx}{\ddot{\theta}}_{x}+{r}_{br}[{C}_{eqL1}({\dot{Z}}_{r1}-{\dot{Z}}_{L1})-{C}_{eqL3}({\dot{Z}}_{r3}-{\dot{Z}}_{L3})+{C}_{eqz1}({\dot{Z}}_{r1}-{\dot{Z}}_{d})-{C}_{eqz3}({\dot{Z}}_{r3}-{\dot{Z}}_{d})]\\ +{r}_{br}[{k}_{eqL1}({Z}_{r1}-{Z}_{L1})-{k}_{eqL3}({Z}_{r3}-{Z}_{L3})+{k}_{eqz1}({Z}_{r1}-{Z}_{d})-{k}_{eqz3}({Z}_{r3}-{Z}_{d})]=0\\ {I}_{ry}{\ddot{\theta}}_{y}+{r}_{br}[{C}_{eqL2}({\dot{Z}}_{r2}-{\dot{Z}}_{L2})-{C}_{eqL4}({\dot{Z}}_{r4}-{\dot{Z}}_{L4})+{C}_{eqz2}({\dot{Z}}_{r2}-{\dot{Z}}_{d})-{C}_{eqz4}({\dot{Z}}_{r4}-{\dot{Z}}_{d})]\\ +{r}_{br}[{k}_{eqL2}({Z}_{r2}-{Z}_{L2})-{k}_{eqL4}({Z}_{r4}-{Z}_{L4})+{k}_{eqz2}({Z}_{r2}-{Z}_{d})-{k}_{eqz4}({Z}_{r4}-{Z}_{d})]=0\\ {I}_{r}{\ddot{\theta}}_{r}-{\displaystyle \sum _{j=1}^{n}({F}_{p\text{\hspace{0.05em}}rj}+{D}_{prj}){r}_{br}+{C}_{rLQ}({\dot{\theta}}_{r}-{\dot{\theta}}_{L})+{k}_{rLQ}({\theta}_{r}-{\theta}_{L})=0}\end{array}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}$ (3)

3) Shield

$\{\begin{array}{l}{m}_{d}{\ddot{X}}_{d}+{C}_{eqr}({\dot{X}}_{d}-{\dot{X}}_{r})+{C}_{eqd}{\dot{X}}_{d}+{k}_{eqr}({X}_{d}-{X}_{r})+{k}_{eqd}{X}_{d}=0\\ {m}_{d}{\ddot{Y}}_{d}+{C}_{eqr}({\dot{Y}}_{d}-{\dot{Y}}_{r})+{C}_{eqd}{\dot{Y}}_{d}+{k}_{eqr}({Y}_{d}-{Y}_{r})+{k}_{eqd}{Y}_{d}=0\\ m{}_{d}\ddot{Z}{}_{d}+{\displaystyle \sum _{i=1}^{4}[{C}_{eqzi}({\dot{Z}}_{d}-{\dot{Z}}_{ri})}+{k}_{eqzi}({Z}_{d}-{Z}_{ri})]+{C}_{eqdz}{\dot{Z}}_{d}+{k}_{eqdz}{Z}_{d}=0\end{array}\text{\hspace{0.05em}}$ (4)

4) Pinion

$\{\begin{array}{l}{m}_{pj}{\ddot{H}}_{pj}+({F}_{p\text{\hspace{0.05em}}rj}+{D}_{p\text{\hspace{0.05em}}rj})\mathrm{cos}\alpha +{C}_{eqpj}{\dot{H}}_{pj}+{k}_{eqpj}{H}_{pj}=0\\ {m}_{pj}{\ddot{V}}_{pj}+({F}_{p\text{\hspace{0.05em}}rj}+{D}_{p\text{\hspace{0.05em}}rj})\mathrm{sin}\alpha +{C}_{eqpj}{\dot{V}}_{pj}+{k}_{eqpj}{V}_{pj}=0\\ {I}_{pj}{\ddot{\theta}}_{pj}+({F}_{p\text{\hspace{0.05em}}rj}+{D}_{p\text{\hspace{0.05em}}rj}){r}_{b\text{\hspace{0.05em}}p}+{C}_{mpQ}({\dot{\theta}}_{pj}-{\dot{\theta}}_{mj})+{k}_{mpQ}({\theta}_{pj}-{\theta}_{mj})=0\end{array}$ (5)

5) Motor

${I}_{mj}{\ddot{\theta}}_{mj}+{C}_{mpQ}({\dot{\theta}}_{mj}-{\dot{\theta}}_{pj})+{k}_{mpQ}({\theta}_{mj}-{\theta}_{pj})={T}_{mj}$ (6)

4. Result Analysis

4.1. Data Statistics

According to the different working conditions under different geological conditions, this paper establishes the relationship between geological and working conditions, obtains the load corresponding to the geology, and substitutes it into the dynamic model of the cutterhead system. The vibration response boundaries of each degree of freedom of the cutter head are obtained by solving the dynamic model. The maximum response intervals of the degree of freedom of the IC, the TPC and the FPC under different working conditions are counted. This paper takes the integrated cutter head as an example, as shown in Table 1.

4.2. Analysis of the Conversion Process of Cutterhead

Firstly, the upper boundaries of the respective degrees of freedom under the various working conditions of the cutterheadare separately summed. Then the average of the upper boundaries of the respective degrees of freedom response boundary is used to represent the response of each cutterhead under various working conditions. Finally, the condition transformation analysis is carried out to obtain the cutterhead stability and cutterhead selection under different geological conditions. In this paper, the first cut-off time is from the beginning of tunneling to the first entry into the Out of difficult condition. The comparative analysis of the conversion differences between different working conditions of different cutterheads are shown in Figure 3.

a) Working conditions conversion(Maximum thrust condition/Maximum overturning condition/Turning and correcting condition/Out of difficult condition); b) Working conditions conversion (Maximum thrust condition/Turning and correcting condition/Maximum overturning condition/Out of difficult condition); c) Working conditions conversion (Maximum overturning condition/Maximum thrust condition/Turning and correcting condition/Out of difficult condition); d) Working conditions conversion (Maximum overturning condition/Turning and correcting condition/Maximum thrust condition/Out of difficult condition).

Through engineering experience, during the normal excavation of the cutterhead, the most common working conditions of the cutterhead are the maximum thrust condition and the maximum overturning condition. Since the geological transformation process corresponds to different vibration conditions, the fluctuation of the vibration difference in different geological transformations can be

Table 1. Statistics on the maximum response boundary of the IC.

Figure 3. Working conditions conversion of different cutterheads.

reflected in the stability of different types of cutterheads during geological conditions conversion. Through the analysis of the transformation process of different working conditions and the analysis of the correspondence between geology and working conditions, we can draw the following conclusions: When the maximum thrust condition and the maximum overturning condition are normal transition (for example, the transition without turning and correcting), that is, the normal transition between the uniform rock stratum and the composite rock stratum, the vibration response of the TPC is the smallest, its working condition transition is the smoothest. In this case, the TPC is preferred. (2) When the maximum thrust condition and the maximum overturning condition are abnormal transitions, the vibration response of the FPC is the smallest, its working condition transition is the smoothest. In this case, the FPC is preferred.

Acknowledgements

This work is supported by the National Natural Science Foundation of China (Grant No. 51875076), NSFC-Liaoning United Key fund (Grant No. U1708255).

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

Cite this paper

Xu, Z.H., Er-ken∙Zha Mu-ti, Y., Xue, L., Wang, F.T., Chen, J. and Huo, J.Z. (2019) Evaluation Method for Tunneling Stability of TBM Cutterhead. World Journal of Engineering and Technology, 7, 1-9. https://doi.org/10.4236/wjet.2019.72B001

References

- 1. Freitag, S., Cao, B.T., Ninic, J., et al. (2017) Recurrent Neural Networks and Proper Orthogonal Decomposition with Interval Data for Real-Time Predictions of Mechanised Tunneling Processes. Computers & Structures. https://doi.org/10.1016/j.compstruc.2017.03.020
- 2. Zhou, Y., Kai, S. and Wu, H.G. (2016) Dynamic Evolution Mechanism of Seepage in Deep-Buried TBM Tunnel. Journal of Central South University.
- 3. Huo, J., Wu, H., Sun, W., et al. (2017) Electromechanical Coupling Dynamics of TBM Main Drive System. Nonlinear Dynamics, No. 2, 1-24.
- 4. Li, X.H., et al. (2013) Research on Dynamic Models and Performances of Shield Tunnel Boring Machine Cutterhead Driving System. Advances in Mechanical Engineering, 29. https://doi.org/10.1155/2013/359757
- 5. Cooley, C.G. and Parker, R.G. (2012) Vibration Properties of High-Speed Planetary Gears with Gyroscopic Effects. Journal of Vibration and Acoustics, 134, Article ID: 061014. https://doi.org/10.1115/1.4006646
- 6. Zhang, K., Yu, H., Liu, Z., et al. (2010) Dynamic Char-acteristic Analysis of TBM Tunneling in Mixed-Face Conditions. Simulation Modelling Practice & Theory, 18, 1019-1031. https://doi.org/10.1016/j.simpat.2010.03.005
- 7. Xian, H.L., Hai, B.Y., Ming, Z.Y., et al. (2010) Dynamic Modeling and Analysis of Shield TBM Cutterhead Driving System. Journal of Dynamic Systems Measurement and Control, 132, 1-14. https://doi.org/10.1115/1.4000818
- 8. Han, M., Cai, Z. and Qu, C. (2014) Study on the Critical Driving Speed of the Tunnel Boring Machine Cutterhead Based on the Dynamic Stability. Journal of Mechanical Engineering, 50, 10-16. https://doi.org/10.3901/JME.2014.21.010