Research Article  Open Access
Zhongkui Zhang, Feng Gao, Yan Li, Han Zhang, "Dynamic and Stability Analysis of the Machine Hydrostatic Slide with Magnetorheological Damper", Shock and Vibration, vol. 2019, Article ID 3605830, 11 pages, 2019. https://doi.org/10.1155/2019/3605830
Dynamic and Stability Analysis of the Machine Hydrostatic Slide with Magnetorheological Damper
Abstract
Tangential dynamic behaviors of the machine hydrostatic slide with a magnetorheological (MR) fluid damper are studied, and the effect of the MR damper to control the vibration of the hydrostatic slide is discussed. The dynamic model of the hydrostatic slide with the MR damper is established, and the tangential vibration equation of linear and nonlinear is derived. The multidimensional incremental harmonic balance method (MIHBM) with discrete Fourier transform (DFT) is derived by which the nonlinear response and stability of the system are studied. The resonance response of the Duffing equation under the combined action of harmonic excitation and constant excitation is obtained. In order to investigate the vibration response of the hydrostatic slide with the MR damper in detail, the bifurcation diagram, phase diagram, and Poincaré map are given. Finally, the dynamic response of the machine hydrostatic slide with the MR damper is discussed, and it is verified that the MR damper can suppress the tangential vibration of the hydrostatic slide effectively and the constant controller can control the chaotic behavior of the system well.
1. Introduction
A hydrostatic slide is widely used in precision machines because of the characteristic of low motion error and high stiffness [1, 2]. The normal vibration characteristics of the joint surface of the hydrostatic slide have been extensively studied by scholars at home and abroad [3–5]. The dynamic characteristics of the machine table is changed instantaneously and abruptly during the startstop, acceleration to constant speed, and constant speed to deceleration because of the small friction coefficient and tangential low damping of the hydrostatic slide [6–8]. In addition, the system is forced to vibrate under the cutting force, ground vibration, residual vibration, and oil film nonlinear fluid fluctuation force which lead to poor stability of the hydrostatic slide in the tangential direction and will affect the positioning accuracy and work effectiveness. However, there are few research reports on the tangential vibration of the hydrostatic slide.
An MR damper and electrorheological (ER) damper are new intelligent vibration suppressors, which have widely been used in industrial and architectural fields in recent years. The magnitude of the force produced by the MR damper is mainly determined by the performance of internal MR fluid. The magnetorheological fluids are intelligent fluids which are very sensitive to magnetic field. When magnetic field is applied to magnetorheological fluids, the viscosity of magnetorheological fluids increases significantly and the MR fluid becomes viscoelastic solid [9]. The output damping force of the MR damper is regulated by input current. In a certain range, the greater the current, the stronger the electromagnetic field and the greater the output damping force. The working principle of the ER damper is similar to that of the MR damper. Kishore et al. [10] developed an electromagnetorheological damper that can be used to control the process dynamics of the cutting tool and reduced cutting force by 12.9%. Tarng et al. [11] used a piezoelectric inertia actuator which is mounted on the cutting tool to suppress chatter, and it was shown that chatter can be effectively suppressed. Fei et al. [12] experimentally verified the problem of the motion damper to suppress chatter in milling flexible parts. Hosseinabadi and Altintas [13] used active damping to effectively suppress the vibration of the machine tool and increase the bandwidth of the drive. Mori et al. [14, 15] proposed a passive viscoelastic damping system which can be applied to various machine tools to reduce the residual vibration. Aoyama and Inasaki [16] developed an electrorheological fluid damper which improved the stability of machine tool table systems with linear motion rolling element bearings. Chen [17] studied the dynamic characteristics of the linear feed system with the MR damper by virtual prototyping technology and proved the system has shorter positioning time, more stable speed characteristics, and better acceleration and deceleration characteristics. Bayat and Pakar [18–20] studied the nonlinear vibration response of beams using the Hamiltonian approach and variational approach, and it was demonstrated that two methods were very accurate and effective and therefore may find wide applicability in engineering and other sciences. In recent years, artificial neural networks were also developed for parameter identification and vibration prediction of nonlinear models [21, 22]. Most of the above studies improved the machining accuracy by suppressing tool chatter using dampers. However, vibration of the hydrostatic slide is also one of the important factors affecting the machining accuracy, and the influence of vibration of hydrostatic guide on processing accuracy has not been studied in the above literature.
In this paper, we studied the dynamic behaviors and stability of the machine hydrostatic slide with the magnetorheological damper. In the next section, the mathematical model of the machine hydrostatic slide with the magnetorheological damper was described and the tangential singlefreedom dynamic model equations are derived. In Section 3, the MIHBM with DFT was derived and its accuracy is verified. Section 4 analyzed the linear and nonlinear dynamic responses and stability of the machine hydrostatic slide with the magnetorheological damper. Suppression of the vibration of the hydrostatic slide by using the magnetorheological damper is studied in detail in Section 5. Finally, some conclusions are drawn briefly in Section 6.
2. Dynamic Model of Hydrostatic Slide with MR Damper
In this section, the dynamic theoretical models of the MR damper and hydrostatic slide system were established, and the singledegreeoffreedom vibration equation of the hydrostatic slide with the MR damper is derived.
2.1. Design and Characterization of MR Damper
Figure 1 shows the structure of the MR damper, which is used to control the vibration of hydrostatic slide. According to the Bingham model of the MR damper [23, 24], the damping force is equivalent to the parallel connection of viscous damping force and Coulomb damping force, i.e.,where is the damping coefficient of MR fluid, is the Coulomb damping force, and is the relative velocity, which can be deduced by the structure of the magnetorheological damper:where is the length of the piston, is the inner diameter of the cylinder, is the effective area of the piston, is the diameter of the piston damping hole, is the yield shear stress, is the viscosity coefficient, and sgn () is the symbolic function.
According to Ampere’s loop law and magnetorheological fluid properties, the relationship between yield stress and current is thatwhere is the number of coil, is the input current, and is the magnetic permeability of the magnetorheological fluid. The substitution of equations (2) and (3) into equation (1) yields
The hysteresis threedimensional model of the MR damper is obtained by parameter identification as shown in Figure 2.
Damping force increases with the increase of input current, but it does not increase linearly. When the current reaches 1.5 A, the output force tends to be saturated, which reflects the characteristics of nonlinear saturation of the MR damper. When the speed changes in different directions, the trajectory of the damping force is different, which reflects the hysteresis characteristics of the magnetorheological damper.
2.2. Dynamic Model of Hydrostatic Slide with MR Damper
The structure of the hydrostatic slide with the magnetorheological damper is shown in Figure 3. The cylinder of the MR damper is fixed. The damper rod is connected with the slide and generates damping force to control vibration of the worktable. The slide moves in the X direction, and the kinematics equation is as follows:where is the structural damping force, is the elastic force, is the friction, is the load, and is the driving force.
The structural damping force is affected by the material and internal resistance of the structure, i.e.,where is the structural damping coefficient. Assume that the friction is linear, i.e.,where is the linear friction coefficient.
In order to analyze the tangential dynamic characteristics of the hydrostatic slide with the magnetorheological damper, the system was simplified to the singledegreeoffreedom mechanical system model as shown in Figure 4. Let , and the vibration equation of the system is expressed aswhere .
3. MIHBM with DFT Method and Feasibility Analysis
3.1. MIHBM with DFT
Wang et al. [25, 26] verified that the dynamic behaviors of the NC table have chaotic characteristics. So the system of the hydrostatic slide with the magnetorheological damper exhibits nonlinear properties. The elastic force [27]: . Let , dimensionless processing, the equation is transformed towhere is the small parameter, .
MIHBM with DFT [28–32] is a new algorithm which replaces the Galerkin procedure [33, 34] with DFT. Equation (9) is a typical Duffing equation under the combined action of harmonic excitation and constant excitation. And the steadystate response of equation (9) is solved by the MIHBM with the DFT algorithm, which improves the computational efficiency.
Suppose , ω_{0} is the solution of the vibration equation, and increments are introduced as follows:
By introducing equation (10) into equation (9) and omitting higherorder small quantities, the incremental equations with and are obtained:
The excitation , response , and its increment can be expressed as the following Fourier series:where , , , , and is the number of harmonics.
According to reference [35], is discretized to a matrix and introducing the differential operator matrix as follows:where is the number of the data points in a period. So the incremental equation (11) can be expressed as
When equation (14) is multiplied by transpose matrix Γ^{T}, the following linearized incremental equation is obtained:where , , , and , and R are the Jacobian matrix, unbalance force matrix, and corrective matrix, respectively.
Adding the constraint equation by the arc length method [31, 35]:where is the arc length parameter. Figure 5 shows the flow chart of the algorithm.
3.2. Feasibility Analysis of MIHBM with DFT
In order to verify the exactness and accuracy of the method proposed in the previous section, Figure 6 shows a comparisons between the approximate analytical solutions based on the MIHBM with DFT and the numerical results based on the Runge–Kutta methods. The vibration waveform and phase diagram of primary resonance response with , , and ω = 1.0 of the hydrostatic slide with the MR damper are shown in Figure 6. NMM represents the numerical results based on the Runge–Kutta method. As demonstrated, the differences between both methods are very slight, so the MIHBM with DFT is feasible to analyze the dynamic behaviors of the hydrostatic slide with the MR damper.
(a)
(b)
4. Dynamic and Stability Analysis of Hydrostatic Slide with MR Damper
In this section, the dynamic behaviors and stability of the hydrostatic slide with the MR damper were studied by MIHBM with DFT.
4.1. Linear System Analysis
Assume the stiffness of the hydrostatic slide is much greater so that it can be regarded as a constant, namely, . Thus, the hydrostatic slide can be regarded as a linear system. Let , dimensionless processing, the equation is transformed to
The steadystate response of equation (17) is
Figure 7 shows the vibration response.
The larger the Coulomb damping force , that is, the greater the output force of the magnetorheological damper, the better the vibration suppression effect. When = 1.2, the hydrostatic slide reaches a stable state quickly. However, it is known that the force produced by magnetorheological damping tends to be saturated when the current increases to a certain value.
The force produced by magnetorheological damping has a maximum value and cannot be increased indefinitely. The amplitudefrequency curve linear system shows that the larger the damping coefficient, the smaller the amplitude of the system, as shown in Figure 8. Therefore, the viscous damping force of the MR damper can also reduce the vibrate displacement of the hydrostatic slide.
4.2. Effect of MR Damper on AmplitudeFrequency
Combined action of harmonic excitation and constant excitation is a prominent feature of vibration equation (9). In this subsection, the effects of constant excitation Coulomb force on the response of system are investigated using the amplitudefrequency curves. Let ξ = 0.015, ε = 1, and F = 0.05, for = 0.2, 0.6, 0.9, 1.1, 1.2, 2.0, respectively, and the amplitudefrequency responses of the system are shown in Figure 9.
(a)
(b)
(c)
(d)
(e)
(f)
It can be seen from Figure 9(a) that the amplitudefrequency curve of the vibration system shows complete hardeningtype nonlinearity when . With the increase in Coulomb damping force to 0.9, the amplitudefrequency curve of the vibration system not only shows hardeningtype nonlinearity but also exhibits slight softeningtype nonlinearity. However, there is no overlap region between hardeningtype nonlinearity and softeningtype nonlinearity when . The amplitudefrequency curve of the vibration system show hardeningtype nonlinearity and softeningtype nonlinearity when , and there is overlap region between hardeningtype nonlinearity and softeningtype nonlinearity. When , the region of the softeningtype nonlinearity covers the region of hardeningtype nonlinearity completely. The amplitudefrequency curve of the vibration system shows the softeningtype nonlinearity and slight hardeningtype nonlinearity.
With the increase in Coulomb damping force F to 2.0, hardeningtype nonlinearity of the amplitudefrequency curve of the vibration system disappears completely and the vibration system exhibits softeningtype nonlinearity.
4.3. Dynamic and Stability Analysis of the Hydrostatic Slide with MRD
The Floquet method [36, 37] is an effective method to determine whether the periodic solution of a nonlinear dynamic system is stable. The periodic solution of the equation is obtained by MIHBM with the DFT method, and by substituting it into equation (11), the following equation is obtained:
Introduce
Equation (20) can be rewritten as
According to the Hsu method proposed by Friedmann et al. [38], periodic T is divided into equallength subintervals and interval size is
The state transition matrix of equation (20) iswhere is the unit matrix and is the number of the polynomial expanded by the exponential matrix.
The eigenvalues and of matrix are the Floquet multipliers of the system. If the modulus of all the Floquet multipliers of the system is less than 1, the solutions are stable and unstable otherwise. Moreover, if a real Floquet multiplier passes through the unit circle along the negative real axis and the rest are located in the unit circle, there is perioddoubling bifurcation.
In this subsection, the steadystate solution of the hydrostatic slide with the magnetorheological damper system is studied by the Floquet method. Figure 10(a) shows the amplitudefrequency curve for , and it exhibits hardeningtype nonlinearity and softeningtype nonlinearity in different regions. And there is overlap region between hardeningtype nonlinearity and softeningtype nonlinearity. Each frequency value of the overlap region corresponds to five amplitude solutions, three of which are stable and two of which are unstable, and stable and unstable solutions alternate. At the lower excitation frequency (i.e., ω ∈ [0, 1.5269]), the change of the vibration system steady response is very slight until the limit point at “2” (ω = 1.5269). With the decrease in frequency ω, the amplitude A becomes unstable until the second limit point reaches “7” (ω = 1.4803). By further increasing the frequency ω, the amplitude of the stable response increases obviously and the response amplitude reaches the peak value at point “9” (ω = 1.5237). With the decrease in frequency ω, the vibration system becomes unstable until the point at “5” (ω = 0.5052). Then, the vibration system returns to stable response once again.
(a)
(b)
With the increase of the excitation frequency ω from point “1” to point “4,” the path of the stable response is 1⟶ 2 ⟶ 3 ⟶ 4. There appears the jumping phenomenon of stable solutions, and the point “2” is the jumped point. With the increase of the excitation frequency from point “4” to point “1”, the path of the stable response is 4 ⟶ 5⟶ 6 ⟶7 ⟶ 8 ⟶ 1. There appears the jumping phenomenon of stable solutions, and the points “5” and “7” are jumped points. The stable solution between points “6” and “9” is achieved by the path 4 ⟶5 ⟶ 6 ⟶ 9 ⟶ 10 ⟶ 4, and the jumped points are “5” and “9.”
Figure 10(b) presents the bifurcation diagram of the hydrostatic slide with the magnetorheological damper system with the parameter ω. At the lower excitation frequency (i.e., ω ∈ [0, 0.06]), the vibration system exhibits chaos motion. In the range of ω ∈ [0.06, 1.87], the vibration system exhibits various types of motion including periodicn and chaotic motions, as shown in Figures 11 and 12. With the increasing of excitation frequency, the vibration system exhibits chaos motion again in the range of ω ∈ [1.87, 3], as shown in Figure 13.
(a)
(b)
(a)
(b)
(a)
(b)
5. The Effect of MR Damper to Control Vibration of Hydrostatic Slide
The relationship between amplitude A and constant excitation Coulomb damping force is shown in Figure 14. For ω = 0.2, 1, 1.49, 1.6, 2.0, respectively, the curve is quite different. At lower Coulomb damping force, there appears the jumping phenomenon of stable solutions. With the increase of Coulomb damping force , the amplitude decrease monotonously for the curve ω = 0.2, 1, 1.49, respectively. The fact is that the curve decreases monotonously as long as ω ≤ 1.49. Moreover, the amplitude decreases rapidly when Coulomb damping force increases to a certain value as shown in Table 1. By further increasing Coulomb damping force , the amplitude tends to zero for the curve ω = 0.2, 1, 1.49, respectively. However, there appears the jumping phenomenon of stable solutions for the curve of ω = 1.6 and 2. It can be seen from the above analysis that the MR damper can reduce the vibration amplitude of the hydrostatic slide effectively.

Figure 15(a) presents the bifurcation diagram of the hydrostatic slide with the magnetorheological damper system with the parameter . The system exhibits chaotic motion after a very short perioddoubling bifurcation. With an increase in , the vibration system displays perioddoubling bifurcation again in the interval [0.76, 1.09]. When > 1.47, the system present perioddoubling bifurcation stably. Figures 15(b) and 15(c) present the phase diagrams and Poincaré diagrams, respectively, when = 1.7. It can be seen that the constant controller can control the chaotic behavior of the system well, which is consistent with that conclusion [39]. So the effect on the magnetorheological damper to control the vibration of the hydrostatic slide is excellent.
(a)
(b)
(c)
6. Conclusions
In this paper, the tangential dynamic characteristics and stability of the machine hydrostatic slide with the MR damper are studied by using a singlefreedom dynamic model. The nonlinear responses the vibration equation are calculated using the IHBM with DFT, and the stability of the solutions is checked using the Floquet method. Finally, the effect of the MR damper to control the vibration of the hydrostatic slide is discussed, and it is verified that the constant controller can control the chaotic behavior of the system well.(1)The MIHBM with DFT is derived which is proven to be feasible to analyze nonlinear dynamic behavior.(2)For the machine hydrostatic slide with the MR damper, the hardening and softeningtype nonlinearities are observed in amplitudefrequency curves. There is an overlap region between hardeningtype nonlinearity and softeningtype nonlinearity when F = 1.1 and each frequency value of the overlap region corresponds to five amplitude solutions. With the increase in frequency, hardeningtype nonlinearity disappears completely and the vibration system exhibits softeningtype nonlinearity. In addition, the maximum amplitude of the machine hydrostatic slide with the MR damper decreases obviously.(3)The bifurcation diagram, phase diagram, and Poincaré map for the vibration system under different excitation frequency are illustrated. With the increasing excitation frequency, the vibration system exhibits more complex dynamic behaviors including wide variety of period motion, chaotic behavior motion, and jumping discontinuous phenomena.(4)The appropriate damping force provided by using the MR damper can reduce the tangential vibration of the hydrostatic slide effectively.(5)Constant controller can control the chaotic behavior of the Duffing system well.
Overall, the researches in this paper can not only understand the dynamic response of the machine hydrostatic slide with the MR damper but also provide some reference to suppress vibration and design and optimize the vibration system.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This work was financially supported by the National Natural Science Foundation of China (NSFC) (No. 51775432), the Science and Technology Project of Weifang University of Science and Technology, China (No. 2018KJWZ08), and Facility Horticulture Laboratory of Universities in Shandong Program (No. 2018YY049).
References
 J. Hwang, C.H. Park, W. Gao, and S.W. Kim, “A threeprobe system for measuring the parallelism and straightness of a pair of rails for ultraprecision guideways,” International Journal of Machine Tools and Manufacture, vol. 47, no. 78, pp. 1053–1058, 2007. View at: Publisher Site  Google Scholar
 P. Zhang, Y. Chen, C. Zhang, J. Zha, and T. Wang, “Influence of geometric errors of guide rails and table on motion errors of hydrostatic guideways under quasistatic condition,” International Journal of Machine Tools and Manufacture, vol. 125, pp. 55–67, 2018. View at: Publisher Site  Google Scholar
 J. Hu, C. Liu, Y. Guo, and X. Liang, “Static and dynamic characteristics of a hydrostatic guideway system based on electrorheological fluid,” Journal of Vibration and Shock, vol. 35, no. 22, pp. 25–30, 2016. View at: Google Scholar
 M. Masuko and T. Nakahara, “The influences of the fluid capacitance in the oil feed line system on the transient response of hydrostatic guideways,” International Journal of Machine Tool Design and Research, vol. 14, no. 3, pp. 233–244, 1974. View at: Publisher Site  Google Scholar
 Z. Wang, W. Zhao, Y. Chen, and B. Lu, “Prediction of the effect of speed on motion errors in hydrostatic guideways,” International Journal of Machine Tools and Manufacture, vol. 64, pp. 78–84, 2013. View at: Publisher Site  Google Scholar
 J. Lee, A. H. Ghasemi, C. E. Okwudire, and J. Scruggs, “A linear feedback control framework for optimally locating passive vibration isolators with known stiffness and damping parameters,” Journal of Vibration and Acoustics, vol. 139, no. 1, Article ID 011006, 2017. View at: Publisher Site  Google Scholar
 J. Lee and C. E. Okwudire, “Reduction of vibrations of passivelyisolated ultraprecision manufacturing machines using mode coupling,” Precision Engineering, vol. 43, pp. 164–177, 2016. View at: Publisher Site  Google Scholar
 D. Yoon and C. E. Okwudire, “Active assist device for simultaneous reduction of heat and vibration in precision scanning stages,” Precision Engineering, vol. 46, pp. 193–205, 2016. View at: Publisher Site  Google Scholar
 A. Fakhar and R. Kolahchi, “Dynamic buckling of magnetorheological fluid integrated by viscopiezoGPL reinforced plates,” International Journal of Mechanical Sciences, vol. 144, pp. 788–799, 2018. View at: Publisher Site  Google Scholar
 R. Kishore, S. K. Choudhury, and K. Orra, “Online control of machine tool vibration in turning operation using electromagneto rheological damper,” Journal of Manufacturing Processes, vol. 31, pp. 187–198, 2018. View at: Publisher Site  Google Scholar
 Y. S. Tarng, J. Y. Kao, and E. C. Lee, “Chatter suppression in turning operations with a tuned vibration absorber,” Journal of Materials Processing Technology, vol. 105, no. 12, pp. 55–60, 2000. View at: Publisher Site  Google Scholar
 J. Fei, B. Lin, S. Yan et al., “Chatter mitigation using moving damper,” Journal of Sound and Vibration, vol. 410, pp. 49–63, 2017. View at: Publisher Site  Google Scholar
 A. H. H. Hosseinabadi and Y. Altintas, “Modeling and active damping of structural vibrations in machine tools,” CIRP Journal of Manufacturing Science and Technology, vol. 7, no. 3, pp. 246–257, 2014. View at: Publisher Site  Google Scholar
 K. Mori, D. Kono, I. Yamaji, and A. Matsubara, “Vibration reduction of machine tool using viscoelastic damper support,” Procedia CIRP, vol. 46, pp. 448–451, 2016. View at: Publisher Site  Google Scholar
 K. Mori, D. Kono, I. Yamaji, and A. Matsubara, “Modelling of viscoelastic damper support for reduction in low frequency residual vibration in machine tools,” Precision Engineering, vol. 50, pp. 313–319, 2017. View at: Publisher Site  Google Scholar
 T. Aoyama and I. Inasaki, “Application of electrorheological fluid dampers to machine tool elements,” CIRP Annals, vol. 46, no. 1, pp. 309–312, 1997. View at: Publisher Site  Google Scholar
 Z. H. Chen, Study on Dynamics Performance of Linear Feed System with MR Damping, Guangdong University of Technology, Guangzhou, China, 2008, in Chinese.
 M. Bayat and I. Pakar, “Nonlinear free vibration analysis of tapered beams by Hamiltonian approach,” Journal of Vibroengineering, vol. 13, no. 4, pp. 654–661, 2011. View at: Google Scholar
 I. Pakar, M. Bayat, and M. Bayat, “On the approximate analytical solution for parametrically excited nonlinear oscillators,” Journal of Vibroengineering, vol. 14, no. 1, pp. 423–429, 2012. View at: Google Scholar
 I. Pakar and M. Bayat, “Analytical study on the nonlinear vibration of EulerBernoulli beams,” Journal of Vibroengineering, vol. 14, no. 1, pp. 216–224, 2012. View at: Google Scholar
 H. Tavakoli and S. Soleimani Kutanaei, “Evaluation of effect of soil characteristics on the seismic amplification factor using the neural network and reliability concept,” Arabian Journal of Geosciences, vol. 8, no. 6, pp. 3881–3891, 2015. View at: Publisher Site  Google Scholar
 S. Rezaei, A. J. Choobbasti, and S. S. Kutanaei, “Site effect assessment using microtremor measurement, equivalent linear method, and artificial neural network (case study: Babol, Iran),” Arabian Journal of Geosciences, vol. 8, no. 3, pp. 1453–1466, 2015. View at: Publisher Site  Google Scholar
 N. M. Wereley, L. Pang, and G. M. Kamath, “Idealized hysteresis modeling of electrorheological and magnetorheological dampers,” Journal of Intelligent Material Systems and Structures, vol. 9, no. 8, pp. 642–649, 1998. View at: Publisher Site  Google Scholar
 C. Graczykowski and P. Pawłowski, “Exact physical model of magnetorheological damper,” Applied Mathematical Modelling, vol. 47, pp. 400–424, 2017. View at: Publisher Site  Google Scholar
 L. H. Wang, Research on Nonlinear Dynamic Characteristics Identifition of NC Table, Huazhong University of Science & Technology, Wuhan, China, 2009, in Chinese.
 L. H. Wang, R. S. Du, B. Wu, and S. Z. Yang, “Nonlinear dynamic characteristics of NC table,” China Mechanical Engineering, vol. 20, no. 13, pp. 1513–1519, 2009, in Chinese. View at: Google Scholar
 Q. Wu, The Study on Dynamics Characteristics and Stability of Feed System in BigInertia, Lanzhou University of Technology, Lanzhou, China, 2012, in Chinese.
 C. Duan and R. Singh, “Superharmonics in a torsional system with dry friction path subject to harmonic excitation under a mean torque,” Journal of Sound and Vibration, vol. 285, no. 45, pp. 803–834, 2005. View at: Publisher Site  Google Scholar
 C. Duan and R. Singh, “Dynamic analysis of preload nonlinearity in a mechanical oscillator,” Journal of Sound and Vibration, vol. 301, no. 3–5, pp. 963–978, 2007. View at: Publisher Site  Google Scholar
 T. C. Kim, T. E. Rook, and R. Singh, “Super and subharmonic response calculations for a torsional system with clearance nonlinearity using the harmonic balance method,” Journal of Sound and Vibration, vol. 281, no. 3–5, pp. 965–993, 2005. View at: Publisher Site  Google Scholar
 X. Kong, W. Sun, B. Wang, and B. Wen, “Dynamic and stability analysis of the linear guide with timevarying, piecewisenonlinear stiffness by multiterm incremental harmonic balance method,” Journal of Sound and Vibration, vol. 346, pp. 265–283, 2015. View at: Publisher Site  Google Scholar
 O. T. Sen, J. T. Dreyer, and R. Singh, “Envelope and order domain analyses of a nonlinear torsional system decelerating under multiple order frictional torque,” Mechanical Systems and Signal Processing, vol. 35, no. 12, pp. 324–344, 2013. View at: Publisher Site  Google Scholar
 S. Wang, L. Hua, C. Yang, X. Han, and Z. Su, “Applications of incremental harmonic balance method combined with equivalent piecewise linearization on vibrations of nonlinear stiffness systems,” Journal of Sound and Vibration, vol. 441, pp. 111–125, 2019. View at: Publisher Site  Google Scholar
 S. Zhou, G. Song, M. Sun, and Z. Ren, “Nonlinear dynamic analysis of a quarter vehicle system with external periodic excitation,” International Journal of NonLinear Mechanics, vol. 84, pp. 82–93, 2016. View at: Publisher Site  Google Scholar
 E. A. de Souza Neto and Y. T. Feng, “On the determination of the path direction for arclength methods in the presence of bifurcations and ‘snapbacks’,” Computer Methods in Applied Mechanics and Engineering, vol. 179, no. 12, pp. 81–89, 1999. View at: Publisher Site  Google Scholar
 H. Wu, Study on Nonlinear Dynamics of Tilting Pad Journal BearingRotor System, Harbin Institute of Technology, Harbin, China, 2007, in Chinese.
 Q. C. Zhang, H. L. Wang, Z. W. Zhu, F. Shen, A. D. Ren, and H. Y. Liu, Theory and Application of Bifurcation and Chaos, Tianjin University Press, Tianjin, China, 2005, in Chinese.
 P. Friedmann, C. E. Hammond, and T.H. Woo, “Efficient numerical treatment of periodic systems with application to stability problems,” International Journal for Numerical Methods in Engineering, vol. 11, pp. 1117–1136, 1997. View at: Publisher Site  Google Scholar
 W. Xu, Numerical Analysis Methods for Stochastic Dynamical System, Science Press, Beijing, China, 2013, in Chinese.
Copyright
Copyright © 2019 Zhongkui Zhang 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.