Renyuan Wu^{1} , Zhencai Zhu^{2} , Guohua Cao^{3}
^{1, 2, 3}School of Mechatronic Engineering, China University of Mining and Technology, Xuzhou 221116, China
^{1, 2, 3}Jiangsu Key Laboratory of Mine Mechanical and Electrical Equipment, China University of Mining and Technology, Xuzhou 221116, China
^{2}Corresponding author
Journal of Vibroengineering, Vol. 17, Issue 2, 2015, p. 978987.
Received 14 December 2014; received in revised form 25 January 2015; accepted 10 February 2015; published 31 March 2015
Copyright © 2015 JVE International Ltd. 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.
The behavior of ropeguided conveyances is so complicated that the ropeguided hoisting system hasn’t been understood thoroughly so far. In this paper, with userdefined functions loaded, ANSYS FLUENT 14.5 was employed to simulate the flowinduced vibration of ropeguided conveyances under different ventilation air speed. With ropeguided mine cages taken into account, results show that the ventilation affects the lateral displacement of conveyance greatly. With the increase of ventilation air speed, the maximum lateral and side displacements of ascending conveyances also increase, while those of descending conveyances don’t always increase, because the ventilation air flows downcast. With the thrust bearings equipped with the hoist rope attachment and the tail rope attachment, the rotation of conveyance about vertical axis is very small.
Keywords: flowinduced vibration, ventilation, ropeguided conveyance.
Shaft hoisting system plays a significant role in underground mining industry, which is used to transport ore, equipment and personnel, and is called “the throat” of shaft. Rope guide is characterized by short installation time, lower capital costs, lower maintenance requirements and smooth travel [1, 2]. Rope guides have been used widely in mining shafts. However, the stiffness of rope guides is very small, so the horizontal displacement of ropeguided conveyances is larger than that of fixedguided conveyances. Aerodynamic forces, residual unbalanced head/tail rope torque and Coriolis force are perceived as the main types of disturbing forces which induce the transverse displacement of ropeguided conveyances [25].
Aerodynamic force results from the asymmetrical air flow around the conveyance. Slonina and Stuhler [1] reported that at Merlebach mine, the ropeguided cage was displaced almost 600 mm towards the shaft wall when two cages were travelling side by side against the ventilation air flow, and they found out that the air flow plays an important role in generating oscillation of ropeguided conveyance. Based on the empirical formula, the magnitude of the steady aerodynamic force depends on the relative velocity between ventilation air speed and conveyance speed [25]. Aerodynamic buffeting force was measured by Hurlin [6] on a scale model and he gave an empirical model to predict the buffeting force. However, the literature about the behavior of ropeguided conveyance is quite less, and the effect of ventilation on the moving ropeguided conveyance has not been reported.
A fluidstructure interaction (FSI) method based on computational fluid dynamic (CFD) was presented to simulate the lateral oscillations of ropeguided conveyances [7]. In this study, we use the FSI method to investigate the behaviour of ropeguided conveyances when different ventilation air speeds are applied to shaft.
The detailed fluidstructure interaction (FSI) technique to simulate the lateral oscillations of ropeguided conveyances has been presented in our previous article [7].
Fig. 1 gives the mathematical model of ropeguided hoisting system. According to the Newton’s second law, the equations of horizontal motion for the conveyance in $x$direction and $y$direction can be represented as:
where $m$ represents the mass of conveyance, $x$ and $y$ represent the displacements of conveyance in $x\text{direction}$ and $y\text{direction}$ respectively,${F}_{x}$ and ${F}_{y}$ represent the disturbing forces on conveyance in $x$direction and $y$direction respectively, and $k$ represents the lateral equivalent spring stiffness for conveyance. The lateral equivalent spring stiffness of ropeguided conveyance can be found in references [4, 7].
The equations of rotation for the conveyance about the vertical axis can be represented as:
where $J$ represents the mass moment of inertia of the conveyance about the vertical axis, $\theta $ represents the angular displacements of conveyance about vertical axis, ${M}_{x}$ represents total disturbing torque exerted on the conveyance about vertical axis, and ${k}_{t}$ represents the equivalent angular spring stiffness for conveyance about vertical axis. In some Chinese mine hoisting practices, the hoist rope attachment and the tail rope attachment are equipped with the thrust bearings, which can eliminate the torsion from hoist rope and tail rope, so the total disturbing torque only takes the aerodynamic moment on the conveyance about vertical axis into consideration.
Fig. 1. Mathematical model of ropeguided hoisting system
Fig. 2. Shaft layout of mine cages
Fig. 2 gives shaft layout of mine cage, and the main parameters of shaft layouts are given in Table 1. The ventilation air speeds are 0, 2 m/s, 4 m/s, 6 m/s and 8 m/s respectively. The timespeed diagram in vertical direction of mine cage is shown in Fig. 3.
The commercial code, ANSYS FLUENT 14.5, has been employed to simulate these cases with PISO scheme. In this study, the ventilation air flow is downcast, so the upper entrance is a velocity inlet boundary condition and the lower exit is pressure outlet boundary condition. All other boundary conditions are no slip at the wall. Userdefined functions (UDFs) were written to define the vertical velocity of conveyances and to solve the equations of horizontal motion for conveyances with the Newmarkbeta method $\text{(}\gamma =$ 1/2, $\beta =$ 1/4) [8]. The movements of conveyance have been implemented using a dynamic mesh with the local cell remeshing method in ANSYS FLUENT. As the local cell remeshing method only affects the tetrahedral cell types in the mesh [9], tetrahedral meshes were generated in the domain. Fig. 4 illustrates the Compute meshes and there are about 4×10^{5} control volumes. And the shearstress transport (SST) $k\omega $ model [10] was used for the turbulent air flow. These five transient cases were resolved using a characteristic time step of 0.002 seconds, and the total time is 16.6 seconds for the hoisting conveyances. There are 8300 time steps in total and 50 iterations for each maximum time step to resolve each case. About 20 h of CPU time in a 4nodes HPC cluster were required to complete each numerical simulation.
Fig. 3. Timespeed diagram in vertical direction
Fig. 4. Compute meshes
Table 1. The main shaft parameters
Hoist ropes

4 off, 13 mm, 0.586 kg/m

Tail ropes

1 off, 26 mm, 2.33 kg/m

Guide ropes

4 off per cage, 24 mm, 3.22 kg/m

Conveyance mass

1.6 t

Hoisting distance

62.4 m

Hoisting speed

4.0 m/s

Payload

0.8 t

Method of tensioning

Weight

Tension load indicated

1000 kg per guide rope

Ventilation direction

Downcast

When two cages pass each other, pressure contours under different ventilation air speed are shown in Fig. 5. Here, ${v}_{in}$ denotes the ventilation air speed.
Fig. 5. Pressure contours under different ventilation air speed
a)${v}_{in}=$0 m/s
b)${v}_{in}=$2 m/s
c)${v}_{in}=\text{4}$ m/s
d)${v}_{in}=$6 m/s
e)${v}_{in}=$8 m/s
Fig. 6. Simulation results for ascending cage 1 in northsouth direction
a) Lateral aerodynamic force
b) Lateral acceleration
c) Lateral velocity
d) Lateral displacement
Fig. 6 shows the lateral aerodynamic force, lateral acceleration, velocity and displacement of ascending cage 1 in northsouth direction, and Fig. 7 shows the lateral aerodynamic force, lateral acceleration, velocity and displacement of descending cage 2 in northsouth direction. Fig. 8 shows the lateral aerodynamic force, lateral acceleration, velocity and displacement of ascending cage 1 in eastwest direction, and Fig. 9 shows the lateral aerodynamic force, lateral acceleration, velocity and displacement of descending cage 2 in eastwest direction. Here, ${F}_{cr}$ denotes the Coriolis force.
Fig. 7. Simulation results for descending cage 2 in northsouth direction
a) Lateral aerodynamic force
b) Lateral acceleration
c) Lateral velocity
d) Lateral displacement
Fig. 8. Simulation results for ascending cage 1 in eastwest direction
a) Side disturbing force
b) Side acceleration
c) Side velocity
d) Side displacement
Fig. 9. Simulation results for descending cage 2 in eastwest direction
a) Side disturbing force
b) Side acceleration
c) Side velocity
d) Side displacement
Fig. 10. Simulation results for ascending cage 1 around $z$axis
a) Aerodynamic torque
b) Angular acceleration
c) Angular velocity
d) Angular displacement
Fig. 11. Simulation results for descending cage 2 around $z$axis
a) Aerodynamic torque
b) Angular acceleration
c) Angular velocity
d) Angular displacement
As Fig. 5 shows, the ventilation directly affects the pressure contours around the cages. As Fig. 69 show, the ventilation has a significant effect on horizontal displacement of conveyance, because the aerodynamic force is directly related to the ventilation air speed. As Fig. 10 and Fig. 11 show, the rotation of conveyance about vertical axis is very small, because the hoist rope attachment and the tail rope attachment equipped with the thrust bearings eliminate the torsion from hoist rope and tail rope.
As Fig. 12 shows, with the increase of ventilation air speed, the maximum lateral and side displacements of ascending conveyances also increase, while those of descending conveyances do not always increase, because the ventilation air flows downcast. And the maximum lateral displacement is much larger than the maximum side displacement, because the lateral aerodynamic forces are much larger than the side aerodynamic forces.
Fig. 12. Maximum horizontal displacement
With the UDFs loaded by Fluent, this paper investigates the behavior of ropeguided conveyances under different ventilation air speed. Results show that the ventilation greatly affects the deflection of ropeguided conveyances. When designing a ropeguided hoisting system, CFD simulation is a useful tool to predict the behavior of ropeguided conveyance.
This work was supported by the Major State Basic Research Development Program of China (973 Program) (Grant No. 2014CB049404), Program for Changjiang Scholars and Innovative Research Team in University (Grant No. IRT1292) and the Project Funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD).
The authors acknowledge Dr. Geoffrey J. Krige, Dr. Malcolm E. Greenway and Rodney S. Hamilton for their counsel and articles.
We are grateful to the Advanced Analysis and Computation Center of CUMT for the award of CPU hours to accomplish this work.