Dynamic Behavior Analysis of a High-Rise Traction System with Tensioned Pulley Acting on Compensating Rope

: In this paper, dynamic characteristics of the symmetrical traction system with tensioned pulley acting on compensating rope is theoretically investigated. Due to the excitations from drum, the traction system will occur longitudinal and transverse vibration. In order to explore the differences between traditional traction system and tensioned traction system with different tensioned methods and seek the optimal method of vibration suppression, the damping cylinder and terminal tension acting on compensating rope between tensioned pulley and ground are placed. Caused by the change of the rope’s property, the system will produce different dynamic responses. Here, the di ﬀ erential-algebraic equations (DAEs) are derived using Hamilton principle. The transverse and longitudinal nonlinear coupling of ropes are considered. The generalized- α method is selected to solve the DAEs. Based on the response characteristics of the system, the time-frequency characteristics with different terminal damping are obtained by CWT (continuous wavelet transform) and FFT (fast Fourier transform). From results, it can be seen that tensioned pulley plays an important role in suppressing transverse and longitudinal vibration of the symmetrical traction system compared with traditional traction system, especially by adding damping cylinder. The amplitude of system decreases exponentially with the increase of the terminal damping acting on tensioned pulley. Different running speeds of tensioned traction system are discussed. The results can inform the development of relevant mitigating strategies to minimize the effects of excessive vibrations.


Introduction
The traction system has been widely used in many lifting fields, such as elevators, mining hoisting cableway, belt conveyors and cranes. The lifting system consists of four major components: driving drum, lifting rope, compensating rope and conveyances. Caused by the elasticity of the ropes, vibration phenomenon will occur under the excitation of the drum, which will even occur large resonance phenomenon. The dynamical characteristic and vibration controlling method of rope drives have been studied by many researchers. Li [1] established an electric pulley cableway model in banana plantation by Lagrange equations. The dynamic characteristics were analyzed theoretically and experimentally. Feng Huirong [2] analyzed the carrying rope bouncing and the free vibration due to the instant elastic energy release at cargo's unhooking based on the theory of string vibration. The results provide some guidance for the installation design, the improvement of the technology and the safety code of practice for forestry aerial cableways. Pi [3] addresses the vibration control of beams subjected to a moving mass and modeled based on Hamiltonian principle. A method combining SMC (sliding mode control) and PPF (positive position feedback control) is proposed to suppress the transverse vibration of the beam when the moving mass is on and off the beam. Liu [4] established an axially moving belt system with input restriction and proposed a boundary control strategy reduce the vibration. Nguyen [5] proposed a control algorithm that suppresses the transverse vibrations of an axially moving belt system by the conjugate gradient method. In this paper, the proposed control method is to regulate the axial transport velocity of the membrane so as to track a desired profile according to which the vibration energy of the membrane at the end of transport decays most quickly. Different from cableway and belt conveyors, this traction system is a typical of lifting system with distributed mass and time-varying length. Therefore, its dynamic characteristics can't be expressed by an analytical solution. For this problem, many researches have been carried out. As for the mining lifting rope system, Stefan Kaczmarczyk [6] established longitudinal and transverse coupled mathematical model base on the Hamilton's principle for a catenary-vertical lifting system to analyze passage through resonances based on a multi-scale method. Response of the transverse and longitudinal motion of ropes with conveyance is shown in [7], which is observed on the fast time scale and the slow time scale. Here, the results reveal adverse dynamic behavior of the catenary largely due to the auto-parametric interactions between the in-and out-of-plane modes and finds out the resonance modal interactions can be avoided to a large extent if the winding velocity is increased to an appropriate level. The transverse response of cable-guided hoisting system with timevarying length and the influence of vessel eccentricity is studied in [8,9]. Here, the transverse response of the moving hoisting conveyance in cable-guided hoisting system is investigated in this paper. In elevator system, the coupled vibration and vibration control method are also investigated by many researchers. Bao [10][11][12] and Zhang [13] established a unilateral transverse-longitudinal coupled vibration model based on Hamilton's principle and solved the governing equations by the Galerkin method. The stability of the parameter vibration of the lifting system and natural frequency characteristics are analyzed. The results show that transient resonance will occur in elevator system when the frequency of transverse excitation of traction system is equal to or close to the natural frequency of the system. Parameter excitation has a greater impact on the longitudinal vibration of the lifting system, but a smaller impact on the transverse vibration. When the frequency of parameter excitation is close to or equal to the natural frequency of the lifting system, parameter transient will also occur in the system which will become beat vibration. In [14], the results of a study of vertical vibrations caused by torque ripple generated at the elevator drive system is presented. Based on the experimental study, the reason of the system resonance is obtained in this paper, which is consistent with the simulation results.
The traction system is composed by multiple components. And here, the rope is regarded as continuum. Therefore, the established equation is differential-algebraic mixed equation with Lagrange multipliers, which can't be solved by ordinary differential solver. For this kind of equation, many solutions are proposed. An accurate calculation with spatial discretization and substructure method is proposed in [15] to accurately calculate dynamic responses of one dimensional structure with multiple constraints. When the geometric matching condition is linear, DAEs (differential algebraic equations) can be transformed into ODEs (ordinary differential equations) by the substructure method. Martin and Olivier [16][17][18][19] proposed the generalized-α scheme for constrained mechanical systems. Here, the DAEs can be solved directly without conversion, which greatly improved the efficiency of solution.
An elevator system is a distributed parameter system, which will greatly increase the difficulty of vibration control. Vibration control in the lifting process is an essential goal. As for transverse vibration suppression method, a nonstationary sliding mode control method is proposed in [20], which can suppress the transverse vibration of ropes and conveyances to a certain degree. Wang [21] constructed a dynamic model of the flexible guided lifting bucket based on the continuum theory and the generalized Hamiltonian principle. In this paper, the adaptive boundary control method is proposed for position regulation and vibration degradation of the system, where adaptation laws and a boundary disturbance observer are formulated to handle system uncertainties. Wang [22] analyzes the traction system of the tensioned compensating rope and traditional suspension compensating rope, and finds out that the longitudinal vibration of the system can restrain to a certain extent by tensioning compensating rope. In [10,12], the vibration control method of lifting system is proposed. Here, the lifting system is regarded as a unilateral system, and the vibration can be suppressed by tensioning unilateral compensating rope.
As for the dynamic interactions between components of setting tension pulley at compensating rope ring, a comprehensive mathematical model of a high-rise elevator system, taking into account the combined lateral stiffness of the roller guides and guide rails is established in [23]. It is found that the autoparametric internal resonance takes place between the transverse in plane and out of plane modes due to the effects of the cubic nonlinear terms and excitation of the building structure. Seiji Watanabe [24] established a mathematical model regarding the rope as a single degree of freedom spring without mass. The results show the vertical motion of the compensating sheave when the elevator in emergency and normal cases. The simulation results were validated by experiments.
Previous works have mainly focused on a unilateral lifting system or a lumped parameter system. In addition, the dynamic model is not complete for the lifting of a systematic traction system with compensating pulley tensioned at the bottom of the compensating rope. In particular, there is no systematic description of the dynamic response of the traction system under the condition of a compensating rope tension. Further, previous works are focused on the control method of transverse vibration of cableway development, belt conveyors, and lifts. The methods of longitudinal vibration suppression of traction systems haven't been studied deeply. For this kind of large height vertical traction system with distributed parameters, it is impossible to realize the full status monitoring of any position of ropes and conveyances in reality. Active control needs status monitoring to obtain the feedback of the whole system's state. Therefore, active control is not applicable for this kind of large height and distributed mass longitudinal traction lifting system. Compared with active control, the passive control method is more stable and reliable, and doesn't need feedback. There is no study on the passive vibration suppression of vertical traction systems so far. Therefore, a passive control method is proposed to suppress system vibration by setting a tensioned pulley with a tensioned force and a damping cylinder at the bottom of the compensating rope. In this paper, a comprehensive coupled mathematical model of a high-rise traction system is established based on Hamilton principle. The dynamical of the traction system can be obtained by generalized-α , which can reveal the vibration displacement of the system and the position of the resonance zone. Additionally, the natural frequency of the traction system can be obtained, which presents the veering phenomenon of the conveyances, lifting and compensating ropes subsystem. By applying different tensioning methods to the tensioning pulley, the dynamic response of traction lifting system under different working conditions can be obtained. Therefore, the effect of this passive vibration suppress method is also shown and discussed. The ideal tensioning method can be obtained to suppress the vibration of the system by comparing, which can be used to develop a relevant strategy to minimize the effects of excessive vibration and improve ride quality.

Description of the System
The schematic diagram and technical solution of traditional and tensioned traction system is shown in Figure 1. As shown in Figure 1a,c, the traditional model is composed of these parts: driving drum, lifting ropes, conveyances, compensating ropes, guide rails, and guide wheels. In a traditional traction system, rotation of the driving drum drives the movement of ropes, and then the up and down movements of the conveyances. The guide wheels are set at the conveyances and limit the swing of conveyances through its matching with the guide rail. The compensating rope is in free suspension. Besides the above parts, as shown in Figure 1b, the tensioned model includes a tensioned pulley with a damping cylinder and tensioned blocks. The tensioned pulley is set at the bottom of compensating ropes, which is driven by compensating ropes. The technical solution of the tensioned traction system is shown in Figure 1d. It is supported by bearing. Further, it is guided vertically by linear slide rail. The tensioned blocks and damping cylinder are respectively placed on the upper and lower part of the tensioned pulley. Using self-gravity of tensioned pulley and tensioned blocks, the tensioned pulley can tension the compensating ropes and limit the swing of the tail of compensating rope. The system is arranged symmetrically. First of all, there are several assumptions made here to restrict the analysis.
1) Since the cable bending stiffness is quite small [25], it is neglected in this paper; 2) The mass per unit length ρi, (i = 1, 2, 3, 4) of ropes are both assumed to be uniform; 3) Due to restriction of the guide rail, rotation of the conveyance is small enough which can be neglected. Therefore, longitudinal and transverse vibration of the conveyance is considered in this paper; 4) The bending sections of the ropes in the groove arc of the driving drum and tensioned pulley are equivalent to that attached to the groove and its motion has no effect on the vibration characteristics of the system. Therefore, its vibration characteristics are ignored here. Therefore, the rope of the system is divided into four sections: left lifting rope, right lifting rope, left compensating rope and right compensating rope [23]. Here, the compensation rope and lifting rope on the same side are complementary. Therefore, equality can be obtained: 3 the length of conveyances and L is the length of guide rail. Because the tensioned pulley is equipped with guide rail, the tensioned pulley is restricted to move in the longitudinal direction only. Here, The potential energy ( ) p E t of traction system is a function of vibration displacement, which is given as follows: where the symbols E and A stand for the modulus of elasticity and cross-sectional area of ropes.  The gravitational potential energy ( ) g E t of traction system is given as follows: During the lifting progress, the cable will be affected by distributed damping force. Therefore, the virtual work ( ) c W t δ done by the non-conservative force on the system can be represented as: where, 1i c and 2i c , Except that, there still exists forces concentrated at the boundary affecting the system. Additionally, those are: the interaction forces between the conveyances and lifting(compensating) ropes and the interaction forces between the tensioned pulley and compensating ropes. Specific Here, the governing of the equation can be derived from Hamilton's principle: Substituting Equations (1)-(4) and Appendix B into the Hamilton's principle Equation (5), the governing equations of the system can be obtained as follows: , 0, 0, , 1 ,2,3,4 Meanwhile, the time boundary conditions and geometric boundary conditions of the system are applied in the above operation process, namely: The longitudinal and transverse force boundary conditions can be obtained as Equation (11).
And the geometric boundary conditions are: Furthermore, governing equations and the boundary condition which consists of geometric boundary conditions and force boundary conditions of the traction system are shown in Equations (11) and (12).
As shown in above section, the governing equation is a mixed partial and ordinary differential equations with infinite dimension in mathematical form. Many parameters are time-varying, so they can't be obtained by analytical method, only by numerical method.
In order to facilitate solving the governing equation, the lifting rope with non-homogeneous boundary condition should be translation in two parts: homogeneous ( ) Substituting Equation (13) into Equation (7), the transverse governing equation can be obtained: The Specific expressions of ( , ) i i Q x t are provided in Appendix C.
As for the transverse displacements of compensating ropes 2 2 ( , ) y x t and 3 3 ( , ) y x t , the transverse excitation mainly comes from the excitation introduced from the boundary of conveyances. Therefore, the excitation is introduced through the boundary. Further, the following transformation is introduced: Substituting Equation (16) into Equation (7), the transverse governing equation of the compensating ropes can be obtained: where, 2 2 , , , Here, a discrete model can be derived via the Galerkin method. We define the one nondimensional parameter as ( ) According to the chain derivation rule, the partial derivatives of time t and space x can be obtained as follows: i i Based on the continuous characteristics of ropes, the dynamic displacement coordinates can be approximated using the assuming modal method (AMM). The transverse and longitudinal displacement of lifting and compensating rope are assumed as follows: , Applying Green's formula yields the discretized equations as follows: (20) where, which can be written as:

M p + C p + K p = Q + N + Φ λ
( ) Φ q represents the constraint equations. u Φ and y Φ is a Jacobian matrix of the constraint equations, which can be written as: The entries of the matrices and force vector of the Equation (20) are listed in Appendix E. Rewritten as a unified matrix form: Equation (23) is a mixed-differential-algebraic equation. Because the constraint equations and their Lagrangian multipliers need to be discrete, there are difficulties with a numerical solution. Here, the generalized-α algorithm is selected to solve the mixed differential-algebraic equation.
The basic iterative form of the algorithm is as follows： Optimal algorithmic parameters m α , f α , γ and β are proposed to second-order ODEs by Chung. The numerical solution is obtained by a predictor-corrector strategy as for the algorithm. At the correction step, it needs to use Newton iterative algorithm to solve the moving unknown force and the generalized coordinates in each integral step simultaneously. The iterative equation can be transformed into the following form: (27) where, q res and Φ res denote the residual vectors of equilibrium and constraints. Further, t S denotes the augmented tangent matrix expressed in the following form: where, Through multi-step iteration, the differential-algebraic equations can be solved, meaning the dynamical responses of the traction system can be obtained.

Study Case
To obtain dynamical responses, numerical simulations were applied to bilateral traction system. The typical parameters for system are defined as following Table 1.  In the following results, the static displacement caused by gravitational potential energy is ignored. As shown in the following figures, the generation of the excitation is due to the nonroundness of the drum, so the frequency of the excitation is related to the running speed of the drum, which can be expressed as /

Comparison between Traditional Traction System and Tensioned Traction System
The model of traditional traction lifting system is shown in Figure 1a. The dynamical behaviour of traditional traction system has been studied by many scholars. In previous studies, most scholars have considered the compensating rope as a time-varying concentrated mass attached to the conveyance. Wu [26] obtained the transverse vibration characteristics of a traditional traction system. In this paper, by assuming the equivalent mass of balance rope with time-varying length lumped on conveyance, a transverse vibration model for mine friction hoist was established in accordance with the Hamilton's principle. Chang [27] established the finite element vibration model of multi-rope friction lifting system with tension balancing device. In this paper, the compensating rope is also equivalent to the mass block added on the container. The compensating rope is directly suspended and not tensioned in traditional traction system, which means its property is only inertia and no stiffness. Therefore, it can be regarded as a time-varying mass attached to the conveyance. The mass of compensating rope is equal to the product of line density and the length of the compensating rope. Therefore, this value varies with time. Using a numerical solution, the dynamic characteristics of the traditional traction system can be obtained. Here, gravitational potential energy and static displacement of the system under terminal tension is ignored. Figure 2 shows the longitudinal vibration at the conveyance of tensioned and traditional traction system. In this figure, it can be seen that longitudinal displacement of tensioned traction system is much smaller than that of traditional traction system, especially in resonance regions. Caused by the change of physical properties of compensation rope, the resonance regions are also varying between two traction system. This can be proven by Figure 3. In this figure, the longitudinal natural frequency of traditional traction system is much higher than that of tensioned traction system, which results in variation of the intersection between frequency of excitation and natural frequency. Furthermore, it affects positions of the resonance regions.  From the result, it can be seen that the longitudinal vibration of lifting ropes can be well restrained in the whole operation stage due to the terminal tension acting on the compensating rope. Especially when the rope resonates, the amplitude of longitudinal displacement of tensioned traction is much smaller than that of traditional system.    Figure 6 shows the comparison of transverse vibration at the conveyance of tensioned and traditional traction system. In this figure, it can be seen that in ascending process, the terminal tension acting on the compensating rope doesn't affect the transverse vibration at the conveyance. But in the descending progress, the application of tension can restrain the transverse vibration to some extent. Especially in the resonance region from 20s to 60 s, it can better suppress the transverse vibration. The left figure is the transverse displacement of tensioned traction system. The right figure is the transverse displacement of traditional traction system. Compared with the traditional traction system, it is found that the application of terminal tension to tension compensating rope can restrain the transverse vibration of lifting ropes, especially in the resonance region, which is consistent with the change trend at the conveyance. The suppression effect can be clearly seen in Regions A and B of Figures 7 and 8. In Region A, it can be seen that the resonance phenomenon, which occurs in traditional traction system, is eliminated with terminal tension acting on compensating rope to make it in tensioned state. Further, the overall transverse amplitude of the system is also reduced a lot.

Different Tension Acting on Tensioned Pulley
Figures 9-11 shows the longitudinal response of the bilateral traction system at the ascending conveyance (Figure 9), descending conveyance ( Figure 10) and tensioned pulley ( Figure 11). Furthermore, the red ,green, blue and yellow line stand for the response with terminal tension of Terminal tension = 7×10 4 Terminal tension = 1×10 5 Terminal tension = 1.5×10 5 Terminal tension = 2×10 5 However, in Figures 9 and 10, the trend of longitudinal displacement is not horizontal. Due to the ascending and descending run of conveyances, the stiffness of lifting and compensating rope is time-varying. Here, gravity is ignored. However, the static displacement caused by the terminal tension acting on the tensioned pulley does still exist. This is the reason why the trend of the longitudinal displacement is declining. Further, it is evident that the angle of inclination is proportional to the magnitude of terminal tension. It can be concluded that the magnitude of terminal tension only affects the static displacement of the system. However, in terms of the characteristics of vibration, increasing the value of terminal tension has no effect on changes of longitudinal vibration characteristics of the bilateral traction system.

Different Damping Acting on Tensioned Pulley
In the above results, the compensating rope is just tensioned by adding terminal tension acting on the tensioned pulley. Compared with the traditional traction system, terminal tension can suppress the longitudinal vibration of the system to a certain extent. However, just as mentioned in Section 3.2.1, only increasing terminal tension can't further suppress the vibration of the bilateral traction system.
Here, in order to suppress the vibration of the bilateral traction system more deeply, the damping cylinder is installed between the tensioned pulley and ground. Further, the dynamical response of system with terminal tension and damping cylinder acting on the tensioned pulley is shown in Figures 12-15. Here, the terminal tension is set as 4

10
× N at all times. The dynamical response of the conveyances and tensioned pulley are shown in Figure 12a,b,c-1, Figure 13a,b,c-1, Figure 14a,b,c-1 and Figure 15a,b,c-1. In order to see the influence of bottom damping on the time and frequency domain of vibration characteristics, the continuous wavelet transform and Fourier transform are applied to the dynamic response, which are shown in Figure  12a,b,c-2, Figure 13a,b,c-2, Figure 14a,b,c-2, Figure 15a,b,c-2 and Figure 12a,b,c-3, Figure 13a,b,c-3, Figure 14a,b,c-3 and Figure 15a,b,c-3, respectively. From Figure 12a,b,c-2, Figure 13a,b,c-2, Figure  14a,b,c-2 and Figure 15a,b,c-2, it can be seen that the maximum energy in the running process, which can correspond to the dynamic response which is shown in Figure 12a,b,c-1, Figure 13a,b,c-1, Figure  14a,b,c-1 and Figure 15a,b,c-1, and can obtain the corresponding frequency. By compare the corresponding frequency of Figure 12a,b,c-2, Figure 13a,b,c-2, Figure 14a,b,c-2 Figure 15a,b,c-2 and Figure 12a,b,c-3, Figure 13a,b,c-3, Figure 14a,b,c-3 and Figure 15a,b,c-3, the amplitude corresponding to the frequency of the maximum energy point can be obtained. By comparing Figure 12a,b,c-1, Figure 13a,b,c-1, Figure 14a,b,c-1 and Figure 15a,b,c-1, Figure 12a,b,c-2, Figure 13a,b,c-2, Figure  14a,b,c-2 Figure 15a,b,c-2, Figure 12a,b,c-3, Figure 13a,b,c-3, Figure 14a,b,c-3 and Figure 15a,b,c-3, the same location and frequency of each figure can be marked one by one correspondingly. Here, we mark the region with the largest energy point and mark it uniformly with Zones A, C and D. There is a special point marked in the graph after fast Fourier transform (FFT) transformation, which is set as Point B. This is the excitation frequency of the system at constant speed. Figure 12 shows time-frequency characteristics of longitudinal response of the traction system at the ascending conveyance. It can be seen that energy concentration points are mainly concentrated in Zone A and Zone C. Zone A is the increase of vibration amplitude caused by the sudden change of running acceleration because the rope is long, and the vibration frequency is about 0.6-0.7 Hz. The increase of amplitude and frequency of Zone C is due to the system's resonance, and the vibration frequency is about 1-1.1 Hz. Under the same frequency excitation, the vibration frequency in the resonance region is about twice of that in the non-resonance region.
Comparing Figure 12a-c, it can be seen that vibration amplitude of Region A is basically completely dissipated due to the increase of damping. Although the energy of Region C also dissipates, the dissipating speed of that is not as fast as that of Region A, so that the maximum energy region of system changes from Region A to Region C with the increase of damping. It can be seen that the suppressing effect of terminal damping on resonance suppression is not ideal. However, the suppressing effect on the non-resonance region is very strong. In the process of increasing damping, the excitation frequency which is shown as Point B is gradually prominent, so the proportion of its frequency's amplitude gradually increases.  Figure 13 shows time-frequency characteristics of longitudinal response of the bilateral traction system at the descending conveyance. Energy concentration points are mainly concentrated in Zone A. The increase of amplitude and frequency of Zone A is due to the resonance of the system. Zone C is the non-resonance region. By comparing the time-domain signal with the frequency-domain signal after the FFT transformation, the suppressing effect of terminal damping on Region C is better than that of Region A. Just as the same phenomenon of Figure 12, in the process of increasing damping, the excitation frequency which is shown as Point B is gradually prominent, and the proportion of its frequency's amplitude gradually increases. From Figures 12 and 13, it can be concluded that terminal damping has a vibration suppression effect on both resonance and non-resonance regions of dynamical response of ascending and descending conveyance in general. However, the suppression effect of terminal damping on resonance region is not as good as that on non-resonance region. Further, the excitation frequency which is shown as Point B is gradually prominent with the increase of terminal damping. Therefore, as shown in Figure 12a,b,c-3 and Figure 13a,b,c-3 by Fourier transform, when the terminal damping is large enough, the amplitude of the resonance frequency and excitation frequency component will be prominent, and the amplitude of other vibration frequency components will be dissipated, even eliminated. The excitation frequency will also be dissipated, but it will never disappear. Figures 14 and 15 show time-frequency characteristics of longitudinal response and rotation of the bilateral traction system at the tensioned pulley. In these figures, there are two resonances occurred in the process, which are also the energy concentration point, mainly because they are affected by the longitudinal vibration of both ascending and descending side comprehensively. Comparing Regions A, C and D, the energy is gradually decreasing due to the increase of terminal damping. In Figure 14, it can be seen that the main frequency affecting the longitudinal vibration of the tensioned pulley is the excitation frequency shown as Point B. Just as with the same conclusion in Figures 12 and 13, the resonance frequency and excitation frequency component will be the major components with the increase of terminal damping.  Figure 16 shows the variation of largest amplitude of longitudinal vibration bilateral tensioned traction system along with time. It is evident that the maximum amplitude of the longitudinal vibration of the system decreases with the increase of damping. Especially in resonance C, the resonance amplitude of the system is greatly reduced due to the damping effect, even eliminate resonance.  Among them, 0 represents the head of the lifting rope, 1 represents the end of compensating rope, 0.5 represents the position where the conveyance sits, and its position is relative. From the following results, almost all dangerous sections are basically concentrated at the head of the lifting rope, the end of the compensating rope and the position of the conveyance. The most dangerous section is at the end of compensating rope, and the secondary hazard section is at the head of the lifting rope. The third place is the conveyance. Caused by damping directly applied to the tensioner where it is closest to the end of compensation rope, it can reduce its proportion as the dangerous section to a certain extent, and the proportion transfer from the end of the compensation rope to the conveyance. However, the end of the compensating rope takes the largest proportion as the dangerous section which won't change with the variation of terminal damping.  Figures 18 and 19 shows the transverse displacement of the ascending and descending conveyances. It is evident that adding damping cylinder at the connection between the tensioned pulley and the ground won't affect the transverse vibration of the bilateral traction system.  In order to obtain the quantitative relationship between the increase of damping and the decrease of vibration amplitude of the system, the response characteristics of tensioned traction systems with different damping are calculated in Figure 20. Here, the damping acting on the tensioned pulley is set from 0 to  In the above results, it can be seen that increasing damping between tensioned pulley and ground both can suppress the longitudinal vibration of the system. When damping of tensioned pulley is in the interval of 0 − 5 2 10 × N/s, the sum of squares and standard deviation of the vibration's amplitude of the tensioned system is decreased rapidly. When damping of tensioned pulley is in the interval of 5 2 10 × − 5 8 10 × N/s, the decrease rate of sum of squares and standard deviation of vibration's amplitude is obviously reduced. When the damping is very large, with the increase of damping, the variance and square sum of the system basically have no change. It can be seen that there is not a linear relationship between the increase of damping of damping cylinder and the decrease of vibration amplitude. The results show that the amplitude of system decreases exponentially with the increase of the terminal damping acting on tensioned pulley.

Discussion on Different Running Speed of Traction System
During the operation of bilateral tensioned system, different running speed will produce different dynamical response. Figures 21-24 shows the longitudinal response of ascending and descending conveyances and tensioned pulley under different running speed. The black, red, green line stand for the system of running speed is 12 m/s, 13.66 m/s, 15 m/s. From the perspective of dynamic response, it can be seen that the resonance positions of the system are also different with the different running speed. However, with the increase of the system speed, it is found that the overall amplitude of the longitudinal response of conveyances and tensioned pulley is decreasing, especially in the resonance region, which even will eliminate resonance phenomenon of the system. Therefore, running speed of the system is not positively related to the vibration amplitude of the system, even in a negative proportion. This is a strange phenomenon. The faster the speed is, the smaller the amplitude is. The slower the speed is, the more violent the vibration is, and the larger the amplitude is. The reason can be obtained through the frequency chart, as shown in Figure 25. This is because the running speed will affect the excitation frequency, and the change of excitation frequency will have a huge impact on the vibration response of the system. When the excitation frequency is close to the natural frequency of the system, there will be resonance in the system, and the vibration amplitude will also increase significantly. Therefore, the change of system operation speed is positively related to the excitation frequency, which affects the relationship between the excitation frequency and the natural frequency, and then affects the response change of the system. Therefore, it is not that the smaller the speed of the system, the smaller the vibration amplitude of the system.

Conclusions
In this paper, the complete mathematical model of bilateral traction system is established based on Hamilton's principle, which takes into account the transverse and longitudinal coupling effect of ropes. The governing equations and boundary conditions are obtained. From the results, it is evident that when the compensating rope is in tensioned state, the amplitudes of longitudinal and transverse vibration of ropes and conveyances will be much smaller than that of traditional traction system when the compensating rope is in free suspended state. Further, and considering that the length of lifting rope and compensating rope is complementary, the veering phenomena appears in the frequency of each order. Two tensioned methods are discussed. From the results, it is found that the terminal tension acting on compensating rope can only change the static displacement of ropes, which won't change the dynamical displacement. Besides, setting damping cylinder between tensioned pulley and ground can suppress the longitudinal vibration and rotation vibration of the tensioned bilateral traction system. Especially at the resonance region, the increase of the damping's value can even eliminate longitudinal resonance region. When the terminal damping is large enough, the amplitude of the resonance frequency and excitation frequency component will be prominent, which can be proven by continuous wavelet transform (CWT) and FFT. The amplitude of system decreases exponentially with the increase of the terminal damping acting on the tensioned pulley, which is a non-linear relationship. Except that, different running speed of tensioned traction system is discussed. It can be seen that the magnitude of vibration amplitude of tensioned traction system is not necessarily proportional to the magnitude of the velocity, but rather, inversely proportional. The prediction can then inform the development of relevant mitigating strategies to minimize the effects of excessive vibrations and improve ride quality.

Appendix D
The trial functions [28] can be set as follows: