Establishment and simulation of dynamic model of backfilling hydraulic support with six pillars
Zhen Tian^{1} , Shuangxi Jing^{2} , Shan Gao^{3} , Jianjun Zhang^{4}
^{1, 3}School of Mechanical and Electrical Engineering, Zhoukou Normal University, Zhoukou, 466000, China
^{2}School of Mechanical and Power Engineering, Henan Polytechnic University, Jiaozuo, 454000, China
^{4}National Energy Extraction Machinery and Equipment Research and Development Center, Zhangjiakou, 075025, China
^{1, 2}Corresponding authors
Journal of Vibroengineering, Vol. 22, Issue 3, 2020, p. 486497.
https://doi.org/10.21595/jve.2019.20512
Received 9 January 2019; received in revised form 14 June 2019; accepted 5 July 2019; published 15 May 2020
JVE Conferences
A backfilling hydraulic support with six pillars used for working face roof support and goaf backfilling in coal mine is designed, and the structure and working principle of the backfilling hydraulic support are described. In order to improve the working stability of backfilling hydraulic support, the differential equations of motion and the state space model of backfilling hydraulic support are established based on Lagrange method and space coordinate system. According to the support structure and related parameters, the differential equation of motion of the system is solved by MATLAB. The influence law of disturbance frequency and amplitude on the top beam vertical vibration, roll and pitch vibration is obtained. The results show that the vertical vibration and roll vibration of the top beam are more severe in the low frequency range. And the degree of vibration gradually decreases as the disturbance frequency increases. As the disturbance amplitude increases, the vibration of the top beam is more severe. The vibration of the backfilling hydraulic support and the deformation distribution nephogram of the top beam are obtained by the finite element analysis, the validity of the dynamic model is verified by finite element simulation. The results provide the basis for the optimization design and the stability evaluation of backfilling hydraulic support.
Keywords: backfilling hydraulic support, dynamics model, disturbance, dynamic characteristics, finite element simulation.
1. Introduction
In order to realize the safe production of coal mine and the green mining of coal mine, people pay more and more attention to the backfilling technology of fully mechanized mining [1, 2]. As one of the important equipment in backfilling mining technology, backfilling hydraulic support can play the role of traditional hydraulic support in fixing support, isolating mining area and goaf. The conveying mechanism and tamping mechanism can be arranged on the backfilling hydraulic support to transport and tamp the filling materials, which can improve the stress state of overlying strata and restrain the roof subsidence [3]. The application of backfilling support in fully mechanized face is shown in Fig. 1.
The performance of backfilling the hydraulic support has an important impact on the normal operation of the relevant equipment, which is related to the safe production of the entire mine. Therefore, many experts and scholars have studied the backfilling of hydraulic supports. Zhou Yuejin analyzed the force of the top beam of the hydraulic support, the movement process of the main components of the backfilling hydraulic support was simulated, and the corresponding motion characteristic curve was obtained [4]. Xu Junming analyzed the structure of the backfilling hydraulic support, established a mechanical relationship model between the backfilling hydraulic support and the surrounding rock [5, 6]. Through the tilting and sliding experiments of hydraulic support on the test bench, Gongpeilin concluded that the mining height, inclination angle, height of center of gravity of the support have important influence on the tilting and sliding of the support [7, 8]. By analyzing the interaction between the hydraulic support and the direct jack, Qi Fangkun established the force model of the fourpillar backfilling support and calculates the working resistance of the pillar and the jack [9]. The above research mainly focuses on the analysis of the structural stress of the backfilling hydraulic support and its relationship with the surrounding rock, and there are few studies on dynamic models and dynamic characteristics of backfilling hydraulic support.
The structure and working principle of backfilling hydraulic support with six pillars are analyzed in this paper. The multifreedom vibration model and the state space model of the backfilling hydraulic support are established. The differential equations are solved by numerical integration method, and the dynamic response of the backfilling hydraulic support under different excitations is obtained.
Fig. 1. Application of backfilling support in working face
2. The structure of backfilling hydraulic support with six pillars
The backfilling hydraulic support with six pillars is a special hydraulic support that integrates the functions of support and backfilling [10, 11].^{}It is mainly composed of telescopic beam, top beam, pillar fourbar linkage mechanism, tamping institution, filling scraper conveyor, support base, etc. The structure and threedimensional model of the backfilling hydraulic support are shown in Fig. 2.
Fig. 2. The Threedimensional model of backfilling hydraulic support with six pillars
The top beam is supported by the front and center pillars and is balanced by a fourbar linkage. The rear top beam and the top beam are hinged to each other, the beam bears the pressure of the roof through the support of pillars and controls the roof subsidence before backfilling. A backfilling conveyor is suspended on the rear top beam to realize backfilling in goaf. The tamping institution is articulated with the support base, the pendulum beam jack controls the tamping mechanism to swing up and down, the tamping jack controls the tamping beam to expand and contract, and the filling material is tamped through the shovel. In order to improve the support strength of the rear top beam, the hydraulic support adds two rear pillars to the back of the support, and the support resistance of the top beam increases, which can reduce the roof subsidence and improve the backfilling effect.
3. Establishment of dynamic model
The dynamic model of the backfilling hydraulic support with six pillars is shown in Fig. 3. ${c}_{1}$, ${c}_{2}$, ${c}_{3}$, ${c}_{4}$, ${c}_{5}$, ${c}_{6}$ are damping coefficient of the backfilling hydraulic support; ${k}_{1}$, ${k}_{2}$, ${k}_{3}$, ${k}_{4}$, ${k}_{5}$, ${k}_{6}$ are stiffness coefficient of are damping coefficient of the backfilling hydraulic support; $m$, ${m}_{1}$, ${m}_{2}$, ${m}_{3}$ are the total mass of the top beam and the concentrated mass of the top beam to the pillar, kg; ${I}_{1}$, ${I}_{2}$ are pitch moment and roll moment of inertia, kg·m^{2}; $a$, $b$, $c$ and $d$ are the width of the top beam and the center of the top beam to the distance between the front, middle and rear pillars, mm; $x$ are vertical displacement of top beam centroid, mm; ${x}_{1}$is vertical displacement of concentrated mass under left front beam, mm; ${x}_{2}$ is vertical displacement of concentrated mass under right front girder, mm; ${x}_{3}$ is vertical displacement of concentrated mass under left middle front beam, mm; ${x}_{4}$ is vertical mass displacement under the right middle front beam,mm;${x}_{5}$ is vertical mass displacement under left rear beam ,mm; ${x}_{6}$ is vertical mass displacement under right rear girder, mm; $\theta $ and $\phi $ are side slope angle and top beam elevation angle of top beam, °; ${x}_{11}{x}_{66}$ are vertical displacement of joint between top beam and pillar, mm.
The backfilling hydraulic support with six pillars is a multistructure system. In order to describe the convenience of the system movement, generalized coordinates can be used to represent the independent parameters in the motion of the system. The vertical displacement of joint between top beam and pillar can be represented by generalized coordinates [12]:
Fig. 3. The dynamic model of backfilling hydraulic support with six pillars
According to the Lagrangian equation [13], the system vibration differential equation is established by using the system kinetic energy $T$, potential energy $V$and energy loss functions $R$.
The kinetic energy of the system is:
$+\frac{1}{2}{m}_{2}{{\dot{x}}^{2}}_{3}+\frac{1}{2}{m}_{2}{{\dot{x}}^{2}}_{4}+\frac{1}{2}{m}_{3}{{\dot{x}}^{2}}_{5}+\frac{1}{2}{m}_{3}{{\dot{x}}^{2}}_{6}.$
Potential energy of the system is:
$+\frac{1}{2}{k}_{4}(x\frac{1}{2}a\theta b\phi {x}_{2}{)}^{2}+\frac{1}{2}{k}_{5}(x+\frac{1}{2}a\theta +c\phi {x}_{3}{)}^{2}+\frac{1}{2}{k}_{5}(x\frac{1}{2}a\theta +c\phi {x}_{4}{)}^{2}$
$+\frac{1}{2}{k}_{6}(x+\frac{1}{2}a\theta +d\phi {x}_{5}{)}^{2}+\frac{1}{2}{k}_{6}(x\frac{1}{2}a\theta +d\phi {x}_{6}{)}^{2}.$
Energy dissipation function is:
$+\frac{1}{2}{c}_{4}(\dot{x}\frac{1}{2}a\dot{\theta}b\dot{\phi}{\dot{x}}_{2}{)}^{2}+\frac{1}{2}{c}_{5}(\dot{x}+\frac{1}{2}a\dot{\theta}+c\dot{\phi}{\dot{x}}_{3}{)}^{2}+\frac{1}{2}{c}_{5}(\dot{x}\frac{1}{2}a\dot{\theta}+c\dot{\phi}{\dot{x}}_{4}{)}^{2}$
$+\frac{1}{2}{c}_{6}(\dot{x}+\frac{1}{2}a\dot{\theta}+d\dot{\phi}{\dot{x}}_{5}{)}^{2}+\frac{1}{2}{c}_{6}(\dot{x}\frac{1}{2}a\dot{\theta}+d\dot{\phi}{\dot{x}}_{6}{)}^{2}.$
Lagrange equation expression is:
Substituting Eqs. (2)(4) into the Lagrangian equation gives the equation of motion of the system as:
Top beam roll motion differential equation:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}+\frac{1}{2}{c}_{4}a\left({x}_{1}{x}_{2}\right)+\frac{1}{2}{c}_{5}a\left({x}_{3}{x}_{4}\right)+\frac{1}{2}{c}_{6}a\left({x}_{5}{x}_{6}\right)\frac{1}{4}\left({k}_{4}+{k}_{5}+{k}_{6}\right){a}^{2}\dot{\theta}.$
Top beam pitching motion differential equation:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}({c}_{4}{b}^{2}+{c}_{5}{c}^{2}+{c}_{6}{d}^{2})\phi {c}_{4}b({\dot{x}}_{1}+{\dot{x}}_{2})+{c}_{5}c({\dot{x}}_{3}+{\dot{x}}_{4})+{c}_{6}d({\dot{x}}_{5}+{\dot{x}}_{6})$
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}+2\left({c}_{4}b{c}_{5}c{c}_{6}d\right)\dot{x}\left({c}_{4}{b}^{2}+{c}_{5}{c}^{2}+{c}_{6}{d}^{2}\right)\dot{\phi}.$
The dynamic differential equations of Eqs. (6)(8) can be expressed in matrix form [14]:
The mass matrix $M$ in the formula is:
The stiffness matrix $K$ in the formula is:
The damping matrix $C$ in the formula is:
In order to facilitate the solution of dynamic differential equations, Eq. (9) can be expressed by state variables, and the mathematical model of system state variables is:
where $x$ is displacement of structure, $\dot{x}$ is velocity of structure, $A$ is matrix of system, $B$ is input matrix, $C$ is output matrix.
4. Simulation analysis
There are 9 degrees of freedom and 14 state variables in the backfilling hydraulic support system, which are: the vertical mass displacement under the beam and its corresponding vertical velocity, the vertical displacement of the top beam centroid position and its movement speed, top beam roll angle and its angular velocity, top beam pitch angle and its angular velocity. Taking a backfilling hydraulic support as the research object, the relevant parameters are shown in Table 1.
Table 1. Related parameters of backfilling hydraulic support with six pillars
Parameter

Value

Total beam mass $m$ / kg

5.0×10^{3}

Concentrated mass of top beam relative to front pillar ${m}_{1}$ / kg

1.5×10^{3}

Concentrated mass of top beam relative to center pillar ${m}_{2}$ / kg

1.7×10^{3}

Concentrated mass of rear pillar relative to rear pillar ${m}_{3}$/kg

2.1×10^{3}

Front base damping coefficient ${c}_{1}$_{}/ (N·s·m^{1})

1.0×10^{3}

Medium base damping coefficient ${c}_{2}$_{}/ (N·s·m^{1})

1.5×10^{3}

Rear base damping coefficient ${c}_{3}$/ (N·s·m^{1})

1.2×10^{3}

Front pillar damping coefficient ${c}_{4}$/ (N·s·m^{1})

400

Neutral pillar damping coefficient ${c}_{5}$_{}/ (N·s·m^{1})

600

Rear pillar damping coefficient ${c}_{6}$_{}/ (N·s·m^{1})

500

Front base stiffness coefficient ${k}_{1}$ / (N·s·m^{1})

1.0×10^{6}

Middle base stiffness coefficient ${k}_{2}$ / (N·s·m^{1})

1.2×10^{6}

Rear base stiffness coefficient ${k}_{3}$_{}/ (N·s·m^{1})

5.0×10^{5}

Front pillar stiffness coefficient ${k}_{4}$_{}/ (N·s·m^{1})

6.0×10^{5}

Neutral pillar stiffness coefficient ${k}_{5}$_{}/ (N·s·m^{1})

5.0×10^{5}

Rear pillar stiffness coefficient ${k}_{6}$_{}/ (N·s·m^{1})

7.0×10^{5}

Top beam pitching moment of inertia $I$_{}/ (N·s·m^{1})

8.2×10^{3}

Top beam roll moment of inertia $J$_{}/ (N·s·m^{1})

9.5×10^{3}

Distance between left and right pillars $a$ / mm

1.6×10^{3}

Top beam center of mass to front pillar distance $b$/ mm

3.0×10^{3}

Top beam center to center pillar distance $c$/ mm

0.5×10^{3}

Top beam center of mass to rear pillar distance $d$/ mm

2.6×10^{3}

The vibration of the coal shearer when cutting the coal and rock, operation of backfilling conveyor and tamping institution will affect the hydraulic support. These vibration sources propagate through coal and rock mass in the form of stress waves to the hydraulic support. Because of the complex nature of coal and rock mass, the disturbance of the hydraulic support has strong randomness. In the dynamic analysis of hydraulic support, the external excitation can be simplified to harmonics with different amplitudes, frequencies and phase angles. According to the research in the literature [1517], it is known that the dominant frequency of the disturbance generated by the coal shearer cutting is within the low frequency bandwidth of less than 120 rad/s, the backfilling conveyor and the compaction of the compactor also have low frequency properties. Based on the system of differential equations and state space models, the fourthorderfiveorder RungeKutta algorithm is used to solve the differential equations. The timedomain response of the vertical vibration of the top beam at 20 rad/s is shown in Fig. 4.The original static equilibrium of the hydraulic support is broken by the disturbance, the top beam has a transient violent vibration within 00.35 s, the maximum amplitude is about 17.7 mm, and then the vibration frequency and the disturbance frequency are gradually consistent under the action of hydraulic damping. The violent vibration easily causes the fatigue of the connecting pin between the top beam and the pillar. In order to reduce the vibration of the top beam, the number and position of hydraulic pillars can be optimized.
Fig. 4. Vertical vibration response of top beam
Fig. 5. Roll and pitch vibration response of top beam
The time domain response of the top beam to pitch and roll vibration is shown in Fig. 5. It can be seen from Fig. 5 that the pitch vibration of the top beam has a great fluctuation under the action of disturbance, and its maximum amplitude is about 5.5×10^{3} rad .The roll vibration is relatively small, and its maximum amplitude is about 4.3×10^{3}rad. Under the continuous action of disturbance the lateral amplitude of the top beam is slightly higher than the pitch amplitude, but the degree of vibration can be maintained in a controllable range. It will not have a destructive effect on the system.
In order to study the influence of disturbance frequency on the hydraulic support, the dynamic characteristics of the top beam are studied when the disturbance frequency varies from 20 to 80 rad/s, and the vibration response of the top beam is obtained as shown in Fig. 6. It can be seen from Fig. 6 that the top beam has a nonlinear violent vibration at the initial moment. As the disturbance frequency increases, the change in the displacement of the top beam caused by the vibration is reduced, and the amplitude eventually exhibits a periodic variation within a certain range.
The vibration of the top beam under different disturbance frequencies is obtained by the statistics of the top beam displacement, the roll angle and the pitch angle when the disturbance is changed in the range of 20 rad/s80 rad/s as shown in Fig. 7. It can be seen that the vertical vibration of the top beam decreases gradually with the increase of the disturbance frequency, but the amplitude of the decrease tends to be gentle; the variation trend of the lateral vibration of the support is similar to that of the vertical vibration, but the variation of pitching vibration with the change of the frequency is not particularly obvious.
Fig. 6. Vertical vibration response of top beam with different disturbance frequency
a) 40 rad/s
b) 60 rad/s
c) 80 rad/s
Fig. 7. The relationship between the beam vibration and disturbance frequency
a) Variation of top beam displacement with disturbance frequency
b) The variation of the top beam angle with disturbance frequency
Fig. 8. Vertical vibration response of top beam with different disturbance amplitude
a) 15×10^{3}^{}N
b) 20×10^{3}^{}N
Fig. 9. The relationship between the beam vibration and disturbance amplitude
a) Variation of top beam displacement with disturbance amplitude
b) Variation of top beam angle with disturbance amplitude
The dynamic characteristics of the system are studied when the external disturbance is 15 kN, 20 kN and 25 kN respectively. The vertical vibration response of the lower beam is obtained as shown in Fig. 8. It can be seen that as the amplitude of the disturbance increases, the vibration of the top beam is more severe.
Based on the statistics of the top beam displacement, roll angle and pitch angle, the vibration of the top beam under different disturbances changes as shown in Fig. 9. Vertical, lateral and pitch vibration of the top beam become more and more severe with the increase of the external disturbance, and the maximum amplitude of increase is more obvious than the effective value. The magnitude of external disturbance has a more significant effect on the vertical and lateral vibration of the top beam, but less on the pitching vibration.
5. Finite element analysis
According to the design drawings, the backfilling hydraulic support is modeled by Solidworks software, and then imported into ANSYS software through special interface. The model is meshed as shown in Fig. 10. When the disturbance frequency is 80 Hz, the deformation distribution of the top beam is shown in Fig. 11. The large deformation area mainly concentrates in the middle of the top beam, and presents a symmetrical variation law.
By applying external disturbances of different frequencies, the variation of the displacement of the top beam with frequency is obtained as shown in Fig. 12. It can be seen that the results obtained by the theoretical model are basically consistent with the results obtained obtained by the finite element method, which proves the correctness of the model.
Fig. 10. Mesh generation of backfilling hydraulic support with six pillars
Fig. 11. Deformation distribution of top beam
Fig. 12. Comparison of the results of two methods
6. Conclusions
A roof support and backfilling equipment for the goaf was designed. The structure of the backfilling hydraulic support with six pillars and its working principle are introduced. The important effect of backfilling hydraulic support on improving the stress state of the overburden layer above the goaf and restraining the subsidence of the roof was expounded.
Based on the Lagrangian method, a differential equation of motion and a state space model of backfilling hydraulic support were established. The fourthfifth order RungeKutta algorithm were used to solve the differential equation, and the dynamic response of the support under different disturbance frequencies and amplitudes was obtained. It is found that the vertical, lateral and pitch vibrations of the top beam fluctuate greatly at the moment of disturbance, then gradually stabilize and exhibit periodic vibration consistent with the disturbance frequency. With the increase of disturbance frequency, the vertical and roll vibration of the top beam decreases gradually, while the change of disturbance frequency has little effect on the pitch vibration of the top beam.
Through the finite element analysis of the model of the backfilling hydraulic support, the deformation distribution nephogram of the top beam is obtained. The results show that the finite element analysis data are consistent with the theoretical analysis data, which proves the correctness of the dynamic model.
Acknowledgements
The authors are grateful for the funding of this work by National Natural Science Foundation of China (No. 51674134) and Scientific and Technological Research Projects of Henan Province (182102210606).
References
 Senapati P. K., Mishra B. K. Design considerations for hydraulic backfilling with coal combustion products (CCPs) at high solids concentrations. Powder Technology, Vol. 229, Issue 6, 2012, p. 119125. [Publisher]
 Karfakis M. G., Bowman C. H., Topuz E. Characterization of coalmine refuse as backfilling material. Geotechnical and Geological Engineering, Vol. 214, Issue 2, 1996, p. 129150. [Publisher]
 Li L., Aubertin M. An improved analytical solution to estimate the stress state in subvertical backfilled stopes. Canadian Geotechnical Journal, Vol. 45, Issue 45, 2008, p. 14871496. [Publisher]
 Zhou Y. J., Zhang J. X., Nie S. J., et al. Mechanical analysis and kinematic simulation of a hydraulic support used in backfilling and the coal mining face. Journal of China University of Mining and Technology, Vol. 41, Issue 3, 2012, p. 366370. [CrossRef]
 Xu J. M., Tan F. Q., Feng J. U., et al. Principle and structure of backfill supports with six props. China Mining Magazine, Vol. 20, Issue 4, 2011, p. 101104. [CrossRef]
 Xu J. M., Zhang J., Zhou N., et al. Effect factors on equivalent mining height in fully mechanized backfilling mining with solid materials. China Coal, Vol. 37, Issue 3, 2011, p. 6668. [CrossRef]
 Gong P., Jin Z. Research on influencing factors of tilt for fullymechanized mining support with large mining height. Journal of Taiyuan University of Technology, Vol. 32, Issue 6, 2001, p. 666669. [CrossRef]
 Jiang D. H., Gong P. L., Du Z. D. Confirmation of rating working resistance of hydraulic support in residual coal mining passed through empty roadway group. Coal Mine Machinery, Vol. 36, Issue 2, 2015, p. 210212. [CrossRef]
 Qi F. K., Zhou Y. J., Wang E. Q., et al. Mechanical analysis of new backfill hydraulic support with four pillars. Mining Research and Development, Vol. 36, Issue 1, 2016, p. 9093. [CrossRef]
 Miao X. X., Zhang J. X. Key technologies of integration of coal mininggangue washingbackfilling and coal mining. Journal of China Coal Society, Vol. 39, Issue 8, 2014, p. 14241433. [CrossRef]
 Zhang Q., Zhang J., Tai Y., et al. The theoretical research on basic characteristics of backfilling hydraulic support and its application. Journal of Mining and Safety Engineering, Vol. 31, Issue 6, 2014, p. 845851. [CrossRef]
 Lu J. N., Mao J., Xie M., et al. Dynamics model of advanced powered support in heading under full support situation. Journal of China Coal Society, Vol. 40, Issue 1, 2015, p. 5057. [CrossRef]
 Fang Z. F., Deng Z. X. Research on the state space model of vehicle vibration system. China Mechanical Engineering, Vol. 16, Issue 4, 2005, p. 353354. [CrossRef]
 Yang S., Wang T. Y. Dynamic modeling and analyzing of a largescale nonlinear resonance screen. Journal of China Coal Society, Vol. 34, Issue 3, 2009, p. 405409. [CrossRef]
 Zhang Q., Zhang J. X., Ju F., et al. Backfillbody’s compression ratio design and control research insolid backfill coal mining. Journal of China Coal Society, Vol. 39, Issue 1, 2014, p. 6471. [CrossRef]
 Tiryaki B., Hekimoglu O. Z. Effects of drum vibration on performance of coal shearers. Transactions of the Institution of Mining and Metallurgy, Vol. 106, 1997, p. 9194. [CrossRef]
 Luan L. J., Fan Fu M., Gao J. S., et al. Design and research on discharge hole of porous bottom unloading conveyor. Coal Mine Machinery, Vol. 37, Issue 11, 2016, p. 57. [CrossRef]