Abstract

The dynamic responses of parallel hoisting system with time-varying length and rigid guidance under drive deviation are investigated considering tension and torsion characteristics of the ropes. The variable-domain three-node elements of rope are employed and the corresponding differential algebraic equations (DAEs) are derived using Lagrange’s equations of the first kind. The slack situation of the rope is considered, and the dynamic equations which are systems of DAEs are transformed to ordinary differential equations (ODEs). The dynamic responses of tension, torsion, and acceleration are analyzed considering radius’ error of the drums, which indicates that the drive deviation between ropes can cause large influence on the tension difference and even cause one of the ropes to slack. However, the torsion of the corresponding rope is active. And unreasonable discordance between ropes should be controlled for the design and manufacture of drum on super deep parallel hoisting system.

1. Introduction

Parallel hoisting system, which has the advantage of large carrying capacity, small rotation of conveyance, excellent security performance, and fine economic performance, works at the manner that one conveyance is hoisted cooperatively by several cables [1]. Parallel hoisting system is widely used due to the increase of lifting heavy and large volume loads, such as elevator hoisting system (Figure 1(a)), mine hoisting system (Figure 1(b)), and parallel suspension platform of vertical shaft in construction.

The parallel hoisting system will vibrate subject to the excitation sources including ripple of the motor, turbulence of transmission device, abrasion of winding drum, and the eccentricity of head sheave. These vibrations will affect the stability of the conveyance and may also cause damage to the ropes, which will further lead to accidents.

Dynamic characteristics of single rope hoisting system with time-varying length have been studied by many researchers. Zhu and Ni investigated the linear lateral dynamics of beam and string with an arbitrarily varying length by Hamilton’s principle and Galerkin’s method [2]. Kaczmarczyk and Ostachowicz described the coupled lateral–longitudinal dynamic response of the cables in deep mine hoisting system with a depth of 2085 m and formulated a discrete mathematical model by Rayleigh–Ritz procedure [3]. Wang et al. investigated the lateral response of the moving hoisting conveyance in cable-guided hoisting system [4]. All of these researches in hoisting system with single cable are the foundations for the dynamic research of parallel hoisting system.

As the application of parallel hoisting systems is increasing, people are paying more and more attention on its research. Zhu and Ren developed a spatial discretization and substructure method to accurately calculate the dynamic responses of single and multiple suspension cables with length-variant distributed-parameter components, which simplified the handle of boundary conditions [5]. References [69] investigated the dynamics of cable-suspended parallel system, but the vibrations of cables themselves are not considered. So, it could not get accurate results. Shao et al. investigated the dynamics of a sinking winch mechanism in the framework of nonsmooth dynamics considering the unilateral property of cable and presented a numerical simulation method which is suitable for the dynamic analysis of the sinking winch mechanism [10]. Du et al. addressed the dynamic modeling of large CDPMs using a variable-domain finite element method; the effects of cable length variation and the resulting mass variation are also considered [1]. Wang et al. established the longitudinal vibration model of parallel hoisting system with tension autobalance device attached to the ends of all hoisting ropes and showed the influence on longitudinal response of different coefficients and excitations [11].

Finite element method (FEM), as a well-known and efficient numerical method, has been widely used in engineering problems [12]. It can be used to deal with the dynamic problem of continuum [13]. Stylianou and Tabarrok analyzed the dynamic characteristic of an axially moving beam [14, 15]. Wang et al. investigated three-dimensional vibrations of an underwater geometrically nonlinear cable with a weight at the lower end [16]. Moustafa et al. considered the modeling problem of the dynamics of overhead cranes with flexible cable and load hoisting or lowering during crane travel [17]. All these researches are investigated by variable-domain FEM (Vd-FEM). The Vd-FEM does not need to find trial functions to satisfy the boundary conditions while shape functions are used as interpolation functions. And the displacement boundary conditions at upper ends of cables can be satisfied accurately by modifying the global matrix properly.

Wire rope, for its complex helical structure, will produce twist motion or torsion under axial load. Costello derived the coupled stiffness coefficient of wire rope, which can be used to express the relation between axial force and torsion with strain [18]. Thus, when ropes vibrate longitudinally, the torsional vibrations will also occurred. The different upper ends excitations will cause different tension and torsion of ropes. The torsion difference will cause the conveyance rotating around the vertical axis, and tension difference in ropes also causes the conveyance rotating around the horizontal axis. Sometimes the rotation of conveyance can be uncomfortable for the passengers or even damage the guide devices. In this paper, the dynamic responses of parallel hoisting system with time-varying length and rigid guidance under drive deviation are investigated considering tension and torsion characteristics of ropes. The variable-domain three-node elements of rope are used and the corresponding differential algebraic equations (DAEs) are derived by Lagrange’s equations of the first kind, which are for working out several challenging problems including the tensions and torsions between the conveyance and cables and the conveyance motion. The dynamic responses of parallel hoisting system are discussed by considering radius’ error of drum.

2. Model for Parallel Hoisting System

2.1. Description of Parallel Hoisting System

Parallel hoisting system depicted in Figure 1 can be simplified as parallel ropes with an attached conveyance which is restrained by rigid guidance, where is the number of hoisting ropes, as shown in Figure 2. Length of each rope is at time excluding the excitation at upper end. The vertically translating velocity and acceleration of ropes are and , where the overdot denotes time differentiation, and the same representation is used as follows.

represents the longitudinal drive deviation on th rope at upper end, which might be caused by the radius’ error or irregular outline of drum. and denote the angle displacements of conveyance about vertical axis and horizontal axis as shown in Figure 2. The following assumptions constrain the analysis. Guidance devices are always keeping in touch with rigid guidance under preload; the whole friction can be neglected; the influence of lateral vibration of rope on the system is small enough to be neglected; the ropes’ mechanical parameters remain constant.

2.2. Equations of Motion

The kinetic energy of the ropes and conveyance is given bywhere is the number of elements in each rope, is density of hoisting rope per unit length, is the moment of inertia of hoisting rope about its central axis, and denotes the th element of the ropes. , , and are mass and the moment of inertia about its vertical and horizontal axis of conveyance, respectively. is the dynamic longitudinal displacement of conveyance. The operator is given by .

In (1), and are given as where and are longitudinal and torsional displacement vectors, is shape function matrix depending on , where , , and .

The th and th elements of th rope are shown in Figure 3, in which and M are the tension and torsion between the two elements.

Considering the dynamic displacement and static elastic deformation, the strain energy of the hoisting system is formulated aswhere and denote the longitudinal and torsional strain, respectively, where the subscript denotes partial differentiation with respect to , and the same representation is used as follows. The static tension and torsion in th rope can be expressed as

The dynamic tension and torsion can be given, respectively, as in which and denote longitudinal and torsional stiffness coefficients of th rope, respectively. and denote the coupled stiffness coefficient [18]. is the equivalent stiffness of guide devices; is the horizontal distance between the center of conveyance and rigid guidance; is the vertical distance between guide devices and the center of conveyance.

Setting the initial position as the place of zero gravitational potential energy, the gravitational potential energy of system can be written as

According to the viscoelastic properties of wire rope [19], the equivalent external nodal load on the th element of th rope is formulated aswhere , , and and are longitudinal and torsional damping coefficients, respectively.

The global equivalent external nodal load can be obtained as and , where denotes the assembly operator of the Vd-FEM.

The damping force of guide device on conveyance is given aswhere is equivalent damping of guide devices.

Since is very small, it can be obtained approximately by the dynamic displacement of the lower ends of ropes aswhere is the distance between adjacent ropes.

The geometric relationships between dynamic longitudinal displacement of conveyance and k ropes at lower ends are given as

The dynamic torsional displacements of lower ends of ropes can be expressed by the rotation angle of conveyance as

Substituting (1)–(10) into Lagrange’s equations of the first kind [20]and using a standard assembly procedure of the FEM [12], the dynamic equations of k ropes and conveyance without the excitations at upper ends can be obtained as which is a system of DAEs, where is the displacement vector of all the nodes on ropes and three degrees of freedom of conveyance.

, , , and are the matrixes of mass, damping, stiffness, and force. , , , and are the global matrixes of ropes assembled from element matrixes, which are listed in the Appendix. And , , , and . is a 2k-order vector which includes the constraint equations (11)–(13); , are Lagrange multipliers.

Considering the following DAEs calculation, according to (11)–(13), the matrix can be written as , where is matrix and is matrix.

Rewriting the displacement vector as , where is a vector and is a vector, hence, can be expressed by aswhere and is identity matrix.

Differentiating (17) twice with respect to t yields

Substituting (17) and (18) into (15) and premultiplying it by , the dynamic equations of ropes and conveyance without excitations at upper ends are formulated aswhich is a system of ordinary differential equations (ODEs), where , , , and .

When longitudinal excitations and torsional constraint are imposed on the upper ends of ropes, the dynamic equations should be converted into the form as follows:where , , and are obtained from , , and by deleting the 1 to rows and the 1 to columns and is obtained from by deleting the 1 to row. And is obtained form by the operation that where .

2.3. Rope Slack Condition

Now consider a conveyance hoisted by two ropes, the dynamics of the system will be more complex if one of ropes is slack due to the deviation of drums. One rope is completely slack or tight which is judged by the tension at the connection point of the conveyance. Before the rope is completely slack, the dynamic response can be obtained by the equations derived before. But, if the slack rope is completely slack, the constraint condition between ropes and conveyance will change.

When one rope is slack, it will not bear loads from the conveyance while the torsion still exists. Then, the dynamic torsional displacements of lower ends of ropes can be expressed by the rotation angle of conveyance as

If rope 1 was slack, the relationship between dynamic longitudinal displacement of conveyance and ropes at lower ends can be given as

For the convenient of dealing with the constraint condition, the sequence of elements in  Q is transformed, and the expression of can be expressed as

The process of transforming DAEs to ODEs and obtaining the solution is the same as above, and it could be solved by Newmark-β method.

3. Results and Analysis

3.1. Parameters

In mine hoisting system, travelling distance ranges is large and continuously being increased. For example, in the South African gold mining industry, the depth of Kloof Gold Mine shaft is 2085 m, and the hoisting distance of South Deep Gold Mine has reached 3000 m. Now in China, new mines with shaft depths in the region of 1000–2000 m, based on parallel hoisting system, have recently been considered, for example, mine shafts in Chihong and Linyi. In this section, two-rope parallel hoisting system is discussed. The initial and maximum lengths of ropes are reset as 30 m and 1566 m, respectively. The downward movement profiles are shown in Figure 4.

Following the regulations that are mentioned above for parameter symbols, system parameters used in calculation are listed in Table 1, in which is the radius of drum.

The longitudinal drive deviation on th rope could be denoted as where and denote swing amplitude of noncircular drum, respectively. The and are expressed as and , in which and are scale coefficients; in this paper, and . The coefficients and denote degree of defect distribution, respectively; in this paper, and . and denote relative deviation induced by radius’ error of drum, respectively. The is expressed as , in which is a scale coefficient.

3.2. Dynamic Response under Drive Deviation

The dynamic responses of the hoisting system are affected by many factors, but the influence of the radius’ error of the hoisting rope is more important because it brings large tension difference. In this paper, the following parametric study will concentrate on the influence of the radius’ error of the hoisting rope and the dynamic responses are investigated. In addition, the variable-domain finite element number is used for the convergence of the solution, and the results are shown in Figures 5~8.

As shown in Figures 5~7, the tension and torsion of both ropes fluctuate greatly at the beginning because of hoisting acceleration, and the fluctuation decreases with the hoisting velocity at constant speed, but resonances occur while the natural frequency is close to the excitation frequency, which is shown in Figure 8. It is obvious that the longitudinal and torsional frequencies about hoisting ropes system decrease as the length increases, but the rotational frequencies about conveyance keep almost unchanged. Intersection points between the natural frequencies and the excitation frequencies denote the resonance locations. It should be noted that there are multiple resonance locations during time interval 60~90 s, which validate the results of Figures 5~7. Although the torsion of ropes has approximately the same trend of vibration with tension in the beginning, there is some difference that the torsions of two ropes increase even if the tension of one rope drops gradually or the rope is slack. The reason for this phenomenon is that the weight of rope along varying length goes up and the tension of upper ends of rope rises gradually; therefore the torsions of the two ropes move up which also can be explained by (5) and (7). The accelerations of conveyance change violently in the horizontal direction during one rope alternating between the slack and tight states.

The dynamic tension difference between the ropes with the coefficient (Figure 7) is changed more than that with the coefficient (Figure 6) and (Figure 5), which indicates that the deviation between ropes caused by the winding radius’ error of drum can bring large influence to the tension difference. It is obvious that even if the deviation was only 1/1000 of the radius, the tension of one rope would dramatically drop, even close to zero, which means that there is an obvious tension imbalance between two ropes. It is worth noting that intermittent tension of one rope would happen and last for a long time, which is shown in Figure 6(a), if the deviation (for example, 1/1000 of the radius) was inappropriate. It should be emphasized that unreasonable discordance between ropes would bring tension imbalance and should be controlled during the design and manufacture of drum, especially in super deep parallel hoisting system.

4. Conclusions

This paper addresses the modeling of parallel hoisting system considering the drive deviation, and the Vd-FEM and Lagrange’s equations of the first kind are used to derive the equations of motion of the system. The dynamic model using Vd-FEM proposed for parallel hoisting system can be readily utilized to solve several challenging problems, including the slack or tight condition of ropes, the tensions and torsions between conveyance and ropes, and the rotation of the conveyance.

It should be noted that the torsions of two ropes increase even if the tension of one rope drops gradually or the rope is slack. The drive deviation (only 1/1000 of the radius) between ropes caused by the winding radius’ error of drum will bring large influence to the tension difference and even lead to slack rope. And unreasonable discordance between ropes should be controlled for the design and manufacturing of drum on super deep parallel hoisting system.

Appendix

The matrix of element mass, damping, stiffness, and force is listed as follows:whereThe matrixes of coefficients are as follows:where the derivatives of shape function about t and x are given as

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

This work is supported by the National Key Basic Research Program of China (2014CB049401), the National Natural Science Foundation of China (51475456), the Fundamental Research Funds for the Central Universities (2014YC06), and the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD).