# Dynamic Response of Parallel Hoisting System under Drive Deviation between Ropes with Time-Varying Length

Authors: Guohua Cao; Xiang Cai; Naige Wang; Weihong Peng; Jishun Li
Journal: Shock and Vibration
 (2017)

---

## 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.


---

## Body

## 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.Figure 1
Parallel hoisting system.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 [6–9] 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 Figure1 can be simplified as k parallel ropes with an attached conveyance which is restrained by rigid guidance, where k is the number of hoisting ropes, as shown in Figure 2. Length of each rope is l(t) at time t excluding the excitation at upper end. The vertically translating velocity and acceleration of ropes are v(t)=l˙(t) and a(t)=l¨(t), where the overdot denotes time differentiation, and the same representation is used as follows.Figure 2
Model of parallel hoisting system.e i ( t ) represents the longitudinal drive deviation on ith rope at upper end, which might be caused by the radius’ error or irregular outline of drum. θc and φc denote the angle displacements of conveyance about vertical axis xc and horizontal axis zc as shown in Figure 2. The following assumptions constrain the analysis. (1) Guidance devices are always keeping in touch with rigid guidance under preload; (2) the whole friction can be neglected; (3) the influence of lateral vibration of rope on the system is small enough to be neglected; (4) the ropes’ mechanical parameters remain constant.


### 2.2. Equations of Motion

The kinetic energy of thek ropes and conveyance is given by(1)T=∑i=1k∑j=1n12∫j-1lt/njlt/nρDuiDt+v2dx+∑i=1k∑j=1n12∫j-1lt/njlt/nJDθiDt2dx+12mcu˙c+v2+12Jcxθ˙c2+12Jczφ˙c2,where n is the number of elements in each rope, ρ is density of hoisting rope per unit length, J is the moment of inertia of hoisting rope about its central axis, and j denotes the jth element of the ropes. mc, Jcx, and Jcz are mass and the moment of inertia about its vertical and horizontal axis of conveyance, respectively. uc is the dynamic longitudinal displacement of conveyance. The operator D/Dt is given by D/Dt=∂/∂t+v(∂/∂x).In (1), ui and θi are given as (2)uix,t=Njx,ltqi,jt,θix,t=Njx,ltpi,jt,xj≤x≤xj+1,j=1,2,…,n,where qi,j(t)=ui,2j-1ui,2jui,2j+1T and pi,j(t)=θi,2j-1θi,2jθi,2j+1T are longitudinal and torsional displacement vectors, Nj(x,l)=Nj1Nj2Nj3 is shape function matrix depending on l(t), where Nj1=2j2-j+nx(1-4j)/l+2(nx/l)2, Nj2=4[j-j2+nx(2j-1)/l-(nx/l)2], and Nj3=1-3j+2j2+nx(3-4j)/l+2(nx/l)2.Thejth and j+1th elements of ith rope are shown in Figure 3, in which T andM are the tension and torsion between the two elements.Figure 3
Three-node element of rope.Considering the dynamic displacement and static elastic deformation, the strain energy of the hoisting system is formulated as(3)Ee=∑i=1k∑j=1n∫j-1lt/njlt/nTic+12Tidεi+Mic+12Midϕidx+8·12kcr1θc2+4·12kcr2φc2,where εi=ui,x and ϕi=θi,x denote the longitudinal and torsional strain, respectively, where the subscript x denotes partial differentiation with respect to x, and the same representation is used as follows. The static tension and torsion in ith rope can be expressed as(4)Tic=ρlt-xg+mcgk,(5)Mic=TicQicQia.The dynamic tension and torsion can be given, respectively, as(6)Tid=Qiaεi+Qibϕi,(7)Mid=Qicεi+Qidϕiin which Qia and Qid denote longitudinal and torsional stiffness coefficients of ith rope, respectively. Qib and Qic denote the coupled stiffness coefficient [18]. kc is the equivalent stiffness of guide devices; r1 is the horizontal distance between the center of conveyance and rigid guidance; r2 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(8)Eg=-∑i=1k∑j=1n∫j-1lt/njlt/nρguidx-mcguc.According to the viscoelastic properties of wire rope [19], the equivalent external nodal load on the jth element of ith rope is formulated as(9)feuj=∫j-1lt/njlt/nfucNjTdx,feθj=∫j-1lt/njlt/nfθcNjTdx,where fuc=ζu(Qiaεi,t+Qibϕi,t), fθc=ζθ(Qicεi,t+Qidϕi,t), and ζu and ζθ are longitudinal and torsional damping coefficients, respectively.The global equivalent external nodal load can be obtained as∑feuj and ∑feθj, where ∑‍ denotes the assembly operator of the Vd-FEM.The damping force of guide device on conveyance is given as(10)Feθc=-8ccr1θ˙c,Feφc=-4ccr2φ˙c,where cc is equivalent damping of guide devices.Sinceφc is very small, it can be obtained approximately by the dynamic displacement of the lower ends of ropes as(11)g1=φc-u1,2n+1-u2,2n+1d=0,where d is the distance between adjacent ropes.The geometric relationships between dynamic longitudinal displacement of conveyance andk ropes at lower ends are given as(12)g2=uc-u1,2n+1-dφck-12=0,gi=ui,2n+1-u1,2n+1-i-1dφc=0,i=3,4,…,k.The dynamic torsional displacements of lower ends of ropes can be expressed by the rotation angle of conveyance as(13)gk+i=θi,2n+1-θc=0,i=1,2,…,k.Substituting (1)–(10) into Lagrange’s equations of the first kind [20](14)ddt∂T∂Q˙i-∂T∂Qi+∂Ee+Eg∂Qi=Fei+∑j=12kλj∂gj∂Qiand using a standard assembly procedure of the FEM [12], the dynamic equations ofk ropes and conveyance without the excitations at upper ends can be obtained as (15)MQ¨+CQ˙+KQ=F+GTλ,(16)gQ,t=0which is a system of DAEs, where Q=u1,1,θ1,1,u2,1,θ2,1,…,uk,2n+1,θk,2n+1,uc,θc,φcT is the displacement vector of all the nodes on k ropes and three degrees of freedom of conveyance.M = M D M c, C=CDCc, K=KDKc, and F=FDTFcTT are the matrixes of mass, damping, stiffness, and force. MD, CD, KD, and FD are the global matrixes of ropes assembled from element matrixes, which are listed in the Appendix. And Mc=diag⁡mc,Jcx,Jcz, Cc=diag(0,8ccr12,4ccr22), Kc=diag(0,8kcr12,4kcr22), and Fc=mcg-a00T. g(Q,t) is a 2k-order vector which includes the constraint equations (11)–(13); G=∂g/∂Q, λ=λ1λ2⋯λ2kT are Lagrange multipliers.Considering the following DAEs calculation, according to (11)–(13), the 2k×2k2n+1+3 matrix G can be written as G=G0G1, where G0 is 2k×(4nk+3) matrix and G1 is 2k×2k matrix.Rewriting the displacement vector asQ=Q0TQ1TT, where Q0 is a 4nk+3 vector and Q1 is a 2k vector, hence, Q can be expressed by Q0 as(17)Q=UQ0,where U=I-G1-1G0 and I is 4nk+3 identity matrix.Differentiating (17) twice with respect tot yields(18)Q˙=UQ˙0,Q¨=UQ¨0.Substituting (17) and (18) into (15) and premultiplying it by UT, the dynamic equations of ropes and conveyance without excitations at upper ends are formulated as(19)M~Q¨0+C~Q˙0+K~Q0=F~,which is a system of ordinary differential equations (ODEs), where M~=UTMU, C~=UTCU, K~=UTKU, and F~=UTF.When longitudinal excitationsui,1=ei(t) and torsional constraint θi,1=0 are imposed on the upper ends of ropes, the dynamic equations should be converted into the form as follows:(20)M-Q-¨0+C-Q-˙0+K-Q-0=F-,where M-, C-, and K- are obtained from M~, C~, and K~ by deleting the 1 to 2k rows and the 1 to 2k columns and Q-0 is obtained from Q0 by deleting the 1 to 2k row. And F- is obtained form F~ by the operation that (21)F-a=F~a+2k-∑p=1kM~a+2k,2p-1e¨pt-∑p=1kC~a+2k,2p-1e˙pt-∑p=1kK~a+2k,2p-1ept,where a=1,2,…,4nk+3-2k.


### 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(22)gi=θi,2n+1-θc=0,i=1,2.If rope 1 was slack, the relationship between dynamic longitudinal displacement of conveyance and ropes at lower ends can be given as(23)g3=φc+2du2,2n+1-uc=0.For the convenient of dealing with the constraint condition, the sequence of elements inQ is transformed, and the expression of G0 can be expressed as (24)G01,42n+2=G02,42n+2=-1,G03,42n+3=-G03,42n+4=2d.The process of transforming DAEs to ODEs and obtaining the solution is the same as above, and it could be solved by Newmark-β method.


## 2.1. Description of Parallel Hoisting System

Parallel hoisting system depicted in Figure1 can be simplified as k parallel ropes with an attached conveyance which is restrained by rigid guidance, where k is the number of hoisting ropes, as shown in Figure 2. Length of each rope is l(t) at time t excluding the excitation at upper end. The vertically translating velocity and acceleration of ropes are v(t)=l˙(t) and a(t)=l¨(t), where the overdot denotes time differentiation, and the same representation is used as follows.Figure 2
Model of parallel hoisting system.e i ( t ) represents the longitudinal drive deviation on ith rope at upper end, which might be caused by the radius’ error or irregular outline of drum. θc and φc denote the angle displacements of conveyance about vertical axis xc and horizontal axis zc as shown in Figure 2. The following assumptions constrain the analysis. (1) Guidance devices are always keeping in touch with rigid guidance under preload; (2) the whole friction can be neglected; (3) the influence of lateral vibration of rope on the system is small enough to be neglected; (4) the ropes’ mechanical parameters remain constant.


## 2.2. Equations of Motion

The kinetic energy of thek ropes and conveyance is given by(1)T=∑i=1k∑j=1n12∫j-1lt/njlt/nρDuiDt+v2dx+∑i=1k∑j=1n12∫j-1lt/njlt/nJDθiDt2dx+12mcu˙c+v2+12Jcxθ˙c2+12Jczφ˙c2,where n is the number of elements in each rope, ρ is density of hoisting rope per unit length, J is the moment of inertia of hoisting rope about its central axis, and j denotes the jth element of the ropes. mc, Jcx, and Jcz are mass and the moment of inertia about its vertical and horizontal axis of conveyance, respectively. uc is the dynamic longitudinal displacement of conveyance. The operator D/Dt is given by D/Dt=∂/∂t+v(∂/∂x).In (1), ui and θi are given as (2)uix,t=Njx,ltqi,jt,θix,t=Njx,ltpi,jt,xj≤x≤xj+1,j=1,2,…,n,where qi,j(t)=ui,2j-1ui,2jui,2j+1T and pi,j(t)=θi,2j-1θi,2jθi,2j+1T are longitudinal and torsional displacement vectors, Nj(x,l)=Nj1Nj2Nj3 is shape function matrix depending on l(t), where Nj1=2j2-j+nx(1-4j)/l+2(nx/l)2, Nj2=4[j-j2+nx(2j-1)/l-(nx/l)2], and Nj3=1-3j+2j2+nx(3-4j)/l+2(nx/l)2.Thejth and j+1th elements of ith rope are shown in Figure 3, in which T andM are the tension and torsion between the two elements.Figure 3
Three-node element of rope.Considering the dynamic displacement and static elastic deformation, the strain energy of the hoisting system is formulated as(3)Ee=∑i=1k∑j=1n∫j-1lt/njlt/nTic+12Tidεi+Mic+12Midϕidx+8·12kcr1θc2+4·12kcr2φc2,where εi=ui,x and ϕi=θi,x denote the longitudinal and torsional strain, respectively, where the subscript x denotes partial differentiation with respect to x, and the same representation is used as follows. The static tension and torsion in ith rope can be expressed as(4)Tic=ρlt-xg+mcgk,(5)Mic=TicQicQia.The dynamic tension and torsion can be given, respectively, as(6)Tid=Qiaεi+Qibϕi,(7)Mid=Qicεi+Qidϕiin which Qia and Qid denote longitudinal and torsional stiffness coefficients of ith rope, respectively. Qib and Qic denote the coupled stiffness coefficient [18]. kc is the equivalent stiffness of guide devices; r1 is the horizontal distance between the center of conveyance and rigid guidance; r2 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(8)Eg=-∑i=1k∑j=1n∫j-1lt/njlt/nρguidx-mcguc.According to the viscoelastic properties of wire rope [19], the equivalent external nodal load on the jth element of ith rope is formulated as(9)feuj=∫j-1lt/njlt/nfucNjTdx,feθj=∫j-1lt/njlt/nfθcNjTdx,where fuc=ζu(Qiaεi,t+Qibϕi,t), fθc=ζθ(Qicεi,t+Qidϕi,t), and ζu and ζθ are longitudinal and torsional damping coefficients, respectively.The global equivalent external nodal load can be obtained as∑feuj and ∑feθj, where ∑‍ denotes the assembly operator of the Vd-FEM.The damping force of guide device on conveyance is given as(10)Feθc=-8ccr1θ˙c,Feφc=-4ccr2φ˙c,where cc is equivalent damping of guide devices.Sinceφc is very small, it can be obtained approximately by the dynamic displacement of the lower ends of ropes as(11)g1=φc-u1,2n+1-u2,2n+1d=0,where d is the distance between adjacent ropes.The geometric relationships between dynamic longitudinal displacement of conveyance andk ropes at lower ends are given as(12)g2=uc-u1,2n+1-dφck-12=0,gi=ui,2n+1-u1,2n+1-i-1dφc=0,i=3,4,…,k.The dynamic torsional displacements of lower ends of ropes can be expressed by the rotation angle of conveyance as(13)gk+i=θi,2n+1-θc=0,i=1,2,…,k.Substituting (1)–(10) into Lagrange’s equations of the first kind [20](14)ddt∂T∂Q˙i-∂T∂Qi+∂Ee+Eg∂Qi=Fei+∑j=12kλj∂gj∂Qiand using a standard assembly procedure of the FEM [12], the dynamic equations ofk ropes and conveyance without the excitations at upper ends can be obtained as (15)MQ¨+CQ˙+KQ=F+GTλ,(16)gQ,t=0which is a system of DAEs, where Q=u1,1,θ1,1,u2,1,θ2,1,…,uk,2n+1,θk,2n+1,uc,θc,φcT is the displacement vector of all the nodes on k ropes and three degrees of freedom of conveyance.M = M D M c, C=CDCc, K=KDKc, and F=FDTFcTT are the matrixes of mass, damping, stiffness, and force. MD, CD, KD, and FD are the global matrixes of ropes assembled from element matrixes, which are listed in the Appendix. And Mc=diag⁡mc,Jcx,Jcz, Cc=diag(0,8ccr12,4ccr22), Kc=diag(0,8kcr12,4kcr22), and Fc=mcg-a00T. g(Q,t) is a 2k-order vector which includes the constraint equations (11)–(13); G=∂g/∂Q, λ=λ1λ2⋯λ2kT are Lagrange multipliers.Considering the following DAEs calculation, according to (11)–(13), the 2k×2k2n+1+3 matrix G can be written as G=G0G1, where G0 is 2k×(4nk+3) matrix and G1 is 2k×2k matrix.Rewriting the displacement vector asQ=Q0TQ1TT, where Q0 is a 4nk+3 vector and Q1 is a 2k vector, hence, Q can be expressed by Q0 as(17)Q=UQ0,where U=I-G1-1G0 and I is 4nk+3 identity matrix.Differentiating (17) twice with respect tot yields(18)Q˙=UQ˙0,Q¨=UQ¨0.Substituting (17) and (18) into (15) and premultiplying it by UT, the dynamic equations of ropes and conveyance without excitations at upper ends are formulated as(19)M~Q¨0+C~Q˙0+K~Q0=F~,which is a system of ordinary differential equations (ODEs), where M~=UTMU, C~=UTCU, K~=UTKU, and F~=UTF.When longitudinal excitationsui,1=ei(t) and torsional constraint θi,1=0 are imposed on the upper ends of ropes, the dynamic equations should be converted into the form as follows:(20)M-Q-¨0+C-Q-˙0+K-Q-0=F-,where M-, C-, and K- are obtained from M~, C~, and K~ by deleting the 1 to 2k rows and the 1 to 2k columns and Q-0 is obtained from Q0 by deleting the 1 to 2k row. And F- is obtained form F~ by the operation that (21)F-a=F~a+2k-∑p=1kM~a+2k,2p-1e¨pt-∑p=1kC~a+2k,2p-1e˙pt-∑p=1kK~a+2k,2p-1ept,where a=1,2,…,4nk+3-2k.


## 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(22)gi=θi,2n+1-θc=0,i=1,2.If rope 1 was slack, the relationship between dynamic longitudinal displacement of conveyance and ropes at lower ends can be given as(23)g3=φc+2du2,2n+1-uc=0.For the convenient of dealing with the constraint condition, the sequence of elements inQ is transformed, and the expression of G0 can be expressed as (24)G01,42n+2=G02,42n+2=-1,G03,42n+3=-G03,42n+4=2d.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 Figure4.Figure 4
Prescribed movement profiles.
(a) (b) (c)Following the regulations that are mentioned above for parameter symbols, system parameters used in calculation are listed in Table1, in which R is the radius of drum.Table 1
System parameters used in calculation.
Parameter Value ρ 11.5 kg/m J 4.518 × 1 0 - 3 kg⋅m2 Q 1 a, Q2a 1.746 × 1 0 8 N Q 2 d, Q2d 1.227 × 1 0 4 N⋅m2 - Q 1 b, Q2b 8.539 × 1 0 5 N⋅m - Q 1 c, Q2c 7.969 × 1 0 5 N⋅m ζ u, ζθ 1.3 × 1 0 - 4 d 0.35 m mc 20 t Jcx 8.313 × 1 0 3 kg⋅m2 Jcz 1.155 × 1 0 4 kg⋅m2 kc 8 × 1 0 5 N/m cc 6.7 × 1 0 4 Ns/m r1 1.1 m r2 0.8 m R 2.5 mThe longitudinal drive deviationei(t) on ith rope could be denoted as (25)eit=esit+erit,i=1,2,where es1(t)=eSsin⁡n1·(lt-l0)/R and es2(t)=eSsin⁡n2·lt-l0/R+α denote swing amplitude of noncircular drum, respectively. The eS and α are expressed as eS=κ·R and α=ς·π, in which κ and ς are scale coefficients; in this paper, κ=0.001 and ς=0. The coefficients n1 and n2 denote degree of defect distribution, respectively; in this paper, n1=1 and n2=2. er1(t)=0 and er2(t)=(er/R)(l(t)-l0) denote relative deviation induced by radius’ error of drum, respectively. The er is expressed as er=ξ·R, 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 ropeer 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 er and the dynamic responses are investigated. In addition, the variable-domain finite element number n=5 is used for the convergence of the solution, and the results are shown in Figures 5~8.Figure 5
(a) Tensions, (b) torsions of the end ropes, and (c) accelerations of conveyance, calculated usingξ=0.0005.
(a) (b) (c)As shown in Figures5~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ξ=0.0015 (Figure 7) is changed more than that with the coefficient ξ=0.001 (Figure 6) and ξ=0.0005 (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.Figure 6
(a) Tensions, (b) torsions of the end ropes, and (c) accelerations of conveyance, calculated usingξ=0.001.
(a) (b) (c)Figure 7
(a) Tensions, (b) torsions of the end ropes, and (c) accelerations of conveyance, calculated usingξ=0.0015.
(a) (b) (c)Figure 8
Longitudinal, torsional, and excitation frequencies.


## 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 Figure4.Figure 4
Prescribed movement profiles.
(a) (b) (c)Following the regulations that are mentioned above for parameter symbols, system parameters used in calculation are listed in Table1, in which R is the radius of drum.Table 1
System parameters used in calculation.
Parameter Value ρ 11.5 kg/m J 4.518 × 1 0 - 3 kg⋅m2 Q 1 a, Q2a 1.746 × 1 0 8 N Q 2 d, Q2d 1.227 × 1 0 4 N⋅m2 - Q 1 b, Q2b 8.539 × 1 0 5 N⋅m - Q 1 c, Q2c 7.969 × 1 0 5 N⋅m ζ u, ζθ 1.3 × 1 0 - 4 d 0.35 m mc 20 t Jcx 8.313 × 1 0 3 kg⋅m2 Jcz 1.155 × 1 0 4 kg⋅m2 kc 8 × 1 0 5 N/m cc 6.7 × 1 0 4 Ns/m r1 1.1 m r2 0.8 m R 2.5 mThe longitudinal drive deviationei(t) on ith rope could be denoted as (25)eit=esit+erit,i=1,2,where es1(t)=eSsin⁡n1·(lt-l0)/R and es2(t)=eSsin⁡n2·lt-l0/R+α denote swing amplitude of noncircular drum, respectively. The eS and α are expressed as eS=κ·R and α=ς·π, in which κ and ς are scale coefficients; in this paper, κ=0.001 and ς=0. The coefficients n1 and n2 denote degree of defect distribution, respectively; in this paper, n1=1 and n2=2. er1(t)=0 and er2(t)=(er/R)(l(t)-l0) denote relative deviation induced by radius’ error of drum, respectively. The er is expressed as er=ξ·R, 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 ropeer 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 er and the dynamic responses are investigated. In addition, the variable-domain finite element number n=5 is used for the convergence of the solution, and the results are shown in Figures 5~8.Figure 5
(a) Tensions, (b) torsions of the end ropes, and (c) accelerations of conveyance, calculated usingξ=0.0005.
(a) (b) (c)As shown in Figures5~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ξ=0.0015 (Figure 7) is changed more than that with the coefficient ξ=0.001 (Figure 6) and ξ=0.0005 (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.Figure 6
(a) Tensions, (b) torsions of the end ropes, and (c) accelerations of conveyance, calculated usingξ=0.001.
(a) (b) (c)Figure 7
(a) Tensions, (b) torsions of the end ropes, and (c) accelerations of conveyance, calculated usingξ=0.0015.
(a) (b) (c)Figure 8
Longitudinal, torsional, and excitation frequencies.


## 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.



---

*Source: 6837697-2017-01-29.xml*