# Analysis on Force Transmission Characteristics of Two-Legged Shield Support under Impact Loading

Authors: Qing-liang Zeng; Zhao-sheng Meng; Li-rong Wan; Cheng-long Wang
Journal: Shock and Vibration
 (2018)

---

## Abstract

To study the load transfer characteristics of a two-legged shield powered support, a numerical simulation model of the support was established using the multibody dynamics software ADAMS. The model took full account of the hydraulic-elastic deformation characteristics of the support, as a series spring-damper system was used to replace the leg and the equilibrium jack. The canopy, goaf shield, lemniscate bars, and equilibrium jack are equivalent to flexible bodies. The setting force of the leg was provided by the preload of the equivalent spring, the static roof load was simulated using a slope signal, and the impact load was simulated using a step signal. Using the model, the impact and excitation effects of each hinge joint of the support were analyzed under different impact load conditions across the canopy. The results show that the location of the impact load affects the force transmissions of all hinge points of the support. Both the impact effect and the excitation effect are at a minimum when the impact force is located near the leg action line. These results are useful for the adaptive control and structural design optimization of the support.


---

## Body

## 1. Introduction

Self-advancing hydraulic support is a technically proven and well-established technique for longwall mining. It helps to hold the face, push the armored face conveyor (AFC), and provide a safe environment for all underground mining operations [1]. Its adaptability to roof loading is a key factor in the safe and efficient production of the working face [2, 3]. The type of loads applied to the shield supports in longwalls can be divided into two categories, namely, static load of the working face (the gravity load of the overlying strata) and impact load [4–6]. Among them, the impact load is the main cause of the support damage. Therefore, studying the dynamic response and force transmission of the shield support under impact loading is helpful for the strength design of the support and safe operations at the working face. Since the effect of impact loading on coal mining operations is very huge, scholars all over the world have studied the influencing factors of the impact load on the shield support and the decisions behind reasonable shield support capacity. For instance, the authors of [7–9] studied the influence of geological conditions on the rock burst tendency using a Finite Element Method. The research reported in [10–15] constructed the numerical model of stope and studied the influence of the rock burst on the support load. Pioneering work presented in [16–19] showed that data from monitoring systems (instrumentation) or shield control systems could be used to predict the occurrence time of the impact load. The supports are treated as mechanic-elastic structures in the above literatures and hydraulic-elastic deformation characteristics of the legs are neglected. Literature [20] studied the dynamic response of the support system when an impact load is applied, and the leg is regarded as a linear elastic structure. Up until now, few literatures have taken the varying characteristics of the leg stiffness into consideration.Given the disparity identified in previous studies, this paper investigates the varying characteristics of the leg stiffness and an equivalent stiffness model of a doubly telescopic leg under isentropic condition was established. Subsequently, a large cutting height shield support was taken as an example. A multibody dynamic model of the support was constructed using ADAMS. By using an equivalent stiffness model, the doubly telescopic leg is equivalent to an elastic element with variable stiffness. The force transfer coefficient variation law at the hinge point of the support was studied considering that the impact load was applied at different locations along its canopy.This paper is arranged as follows. In Section2, the derivation of the equivalent stiffness model is presented. This includes a holistic analysis of the working process of a shield support along with the mathematical formulation of a doubly telescopic leg. In Section 3, a numerical multibody dynamics model of the selected shield support is established. A series spring-damper system (based on the stiffness model described in Section 2) is used to replace the leg and the equilibrium jack. The model is then tested to verify whether the simulations satisfy the precision requirements for the subsequent analysis. In Section 4, the influence of the impact loading position on the hinge point force is analyzed. Finally, conclusions are drawn in Section 5.


## 2. Equivalent Stiffness Analysis Model of a Doubly Telescopic Leg

### 2.1. Analysis of the Working Process of a Shield Support

A typical working cycle of a shield support is the period from the time it is set against the roof to the time when it is released for advancing after the shearer has made its pass [18]. During this process, the role of the supports is to control the roof from an excessive roof-to-floor convergence X, as shown in Figure 1. X is determined by the displacement of the leg h.Figure 1
Supporting schematic diagram of a shield support.Previous studies have analyzed the working process of the support according to its rate of pressure change, segmented into three sections [16, 18]. In this paper, the working process of the support, with doubly telescopic legs (Figure 2), is redefined as five segments from the point of view of its stiffness variations. As can be seen from Figure 2, the leg is equivalent to a series spring model (the calculation of its stiffness is given in Section 2.2). Y is the length of the mobile column, y1 (y2) is the length of the emulsion in the bottom stage (top stage), and D1 (D2) and d1 (d2) are the external and internal diameters of the bottom stage (top stage), respectively.Figure 2
The doubly telescopic leg and its equivalent spring model.Figure3 shows the working curve of a shield support during a whole working cycle. X0~X4 are the displacements of the roof and the support at different stages, respectively. PO′ is the active setting pressure of the support, PO is the effective pump pressure (passive setting pressure), and Pn is the yielding pressure of the relief valve. ka (kb) is the stiffness of the support at stage b (stage c). Figure 4 represents the working process of the doubly telescopic leg (before reaching Pn) by using the equivalent spring model. h1 (h2) is the displacement of the top stage (bottom stage).Figure 3
Working curve of the shield support.Figure 4
Working process of the doubly telescopic leg.
(a)
Stage a (b)
Stage b (c)
Stage cAccording to Figures3 and 4, the five segments of the working process of the support are as follows. The first segment is the active setting (stage a). The shield support is set against the roof using PO. Due to the pressure loss of the bottom valve and smaller diameter of the top stage, the active setting force PO′ formed by the leg is usually less than PO. In this portion, the roof loading P is transferred to the support gradually (Figure 4(a)). However, the load is smaller than PO′ and the leg does not retract (X0 = 0). When P is bigger than Po′ (Figure 4(b)), the top stage begins to retract with the stiffness of ka (stage b). As the roof pressure continues to rise, the support enters the third segment—the rapid increase (stage c). At this time, both stages of the leg retract, resulting in the reduction of the equivalent stiffness (Figure 4(c)). As the roof pressure continues to rise to Pn, the support enters into the overflow segment. At this point, the retracting speed of the support is determined by the performance of the relief valve. Displacement formed in this stage is irreversible; therefore, it can be referred to as the “plastic deformation segment.” The final segment (stage e) is the release segment. The roof loading applied to the hydraulic support decreased rapidly and the leg retracted quickly.The ground response curve (GRC) theory [19, 21–24] suggests that the factor that determines the suitability of the shield support is the stiffness of the support. It holds the view that if the GRC and support reaction line intersect before stage d, the selected support is considered to be reasonable and effective. This type of stiffness analysis of the support was reported previously [25], but the plastic deformation characteristic of the support was neglected in the study.


### 2.2. Mathematical Model of an Equivalent Stiffness Doubly Telescopic Leg

Assuming that the compression deformation process is an isobaric expansion under the isentropic condition, the stiffness of a single-telescopic cylinder can be determined using the following equation:(1)Kz=11/Kl+1/Kp=1hl+he/Fp+Y/EhAh,where Kz is the equivalent stiffness of the single-telescopic leg, Kl is the equivalent stiffness formed by cylinder expansion and emulsion compression, Kp is the stiffness of the piston, hl is the compression length of the emulsion, he is the displacement of the rod caused by cylinder wall expansion, Fp is the axial force applied on the piston rod, Eh is the elastic modulus of the rod, and Ah is the cross-sectional area of the rod.h e can be calculated using the deformation theory of thick-walled cylinders (the leg is a typical high pressure thick wall structure). Combining the force equilibrium equations and geometric equations of the elements, the deformation equation of the cylinder can be expressed as(2)d2udr2+1rdudr-ur2=0,where u is the radial displacement of the cylinder and r is diameter of the cylinder. Substituting the boundary condition (see (3)) into (2) and ignoring the atmospheric pressure, the analytical solution for the radial displacement of the leg can be expressed as (4).(3)σr=d=p,σr=D=p2,(4)u=Δpdm2Eg,where m=1+μD2+(1-μ)d2/D2-d2, Δp=p-p2≈p, μ is Poisson’s ratio, d is the internal diameter of the cylinder, D is the external diameter of the cylinder, p is the cylinder pressure, p2 is the atmospheric pressure, and Eg is the elastic modulus of the cylinder. Then the expansion volume of the leg ΔVg can be expressed as (5) and he can be obtained using (6).(5)ΔVg=πd2yΔpm2Eg,(6)he=ΔVgA=2yΔpmEg,where A is the internal cross-sectional area of the cylinder. hl can be obtained using the compressible fluid theory (see (7)).(7)hl=4κV0Δpπd2=κΔpy,where κ is the compressibility of the emulsion and V0 is the initial volume of the emulsion.By combining (1), (6), and (7), the formula for calculating the equivalent stiffness of ith stage cylinder is given by(8)Kzi=ΔFzihi=πEgAh4κyi/di2Eg+8yimi/di2Ah+πY,where Kzi is the equivalent stiffness of ith stage cylinder, ΔFzi is the larger part of the load experienced by stage i, and hi is the displacement of the rod of stage i. If the volume force acting on the fluid is neglected, ΔFzi=ΔFp.For a doubly telescopic cylinder, it can be regarded as a series spring system (Figure2). Then its stiffness can be calculated as [18](9)K=Kz1Kz2Kz1+Kz2,where K is the equivalent stiffness of doubly telescopic leg.


### 2.3. Case Study

A certain type of shield support was taken as an example in this section. The main parameters of the leg are shown in Table1 (Figure 2) and its equivalent stiffness was calculated using (9). It is assumed that Po′ is 34 MPa (the active setting force is about 3856 kN). The results obtained by using (9) are shown in Figure 5. Lines OA, AB, and CD represent the stiffness of the leg at stage a, stage b, and stage c, respectively.Table 1
The main parameters of the doubly telescopic leg.
Name of the project Top stage Bottom stage Mobile column Structure lengthY (mm) / / 2747 Initial emulsion lengthyi (mm) / 2113 2197 Quality(t) 2.1 1.9 1.276 External/internal diameterDi/di (mm) 625/530 500/380 355 Setting/yielding pressurePo/Pn (MPa) 37.5/47.6 34/92.6 / Modulus of elasticity (GPa) 196 (20°C) Poisson’s ratio 0.3Figure 5
Stiffness variation curve of the doubly telescopic leg.The above settings are substituted into (10) and the calculated results are compared with the actual parameters of the support.(10)FPO˙=FPO′+∫OBKzdx=8394kN,FPn˙=FPO′+∫ODKzdx=10683kN,where FPO· is the calculated value of the setting force and FPn· is the calculated value of the yielding resistance. The error estimate is shown in (11). Given the fact that the errors are less than 5%, the calculation results of formula (9) are deemed to be of high accuracy.(11)δFPO=FPO˙-FPOFPO=1.462%,δFPn=FPn˙-FPnFPn=1.743%,where FPO is the rated setting force (8273 kN) and FPn is the rated working resistance (10500 kN).


## 2.1. Analysis of the Working Process of a Shield Support

A typical working cycle of a shield support is the period from the time it is set against the roof to the time when it is released for advancing after the shearer has made its pass [18]. During this process, the role of the supports is to control the roof from an excessive roof-to-floor convergence X, as shown in Figure 1. X is determined by the displacement of the leg h.Figure 1
Supporting schematic diagram of a shield support.Previous studies have analyzed the working process of the support according to its rate of pressure change, segmented into three sections [16, 18]. In this paper, the working process of the support, with doubly telescopic legs (Figure 2), is redefined as five segments from the point of view of its stiffness variations. As can be seen from Figure 2, the leg is equivalent to a series spring model (the calculation of its stiffness is given in Section 2.2). Y is the length of the mobile column, y1 (y2) is the length of the emulsion in the bottom stage (top stage), and D1 (D2) and d1 (d2) are the external and internal diameters of the bottom stage (top stage), respectively.Figure 2
The doubly telescopic leg and its equivalent spring model.Figure3 shows the working curve of a shield support during a whole working cycle. X0~X4 are the displacements of the roof and the support at different stages, respectively. PO′ is the active setting pressure of the support, PO is the effective pump pressure (passive setting pressure), and Pn is the yielding pressure of the relief valve. ka (kb) is the stiffness of the support at stage b (stage c). Figure 4 represents the working process of the doubly telescopic leg (before reaching Pn) by using the equivalent spring model. h1 (h2) is the displacement of the top stage (bottom stage).Figure 3
Working curve of the shield support.Figure 4
Working process of the doubly telescopic leg.
(a)
Stage a (b)
Stage b (c)
Stage cAccording to Figures3 and 4, the five segments of the working process of the support are as follows. The first segment is the active setting (stage a). The shield support is set against the roof using PO. Due to the pressure loss of the bottom valve and smaller diameter of the top stage, the active setting force PO′ formed by the leg is usually less than PO. In this portion, the roof loading P is transferred to the support gradually (Figure 4(a)). However, the load is smaller than PO′ and the leg does not retract (X0 = 0). When P is bigger than Po′ (Figure 4(b)), the top stage begins to retract with the stiffness of ka (stage b). As the roof pressure continues to rise, the support enters the third segment—the rapid increase (stage c). At this time, both stages of the leg retract, resulting in the reduction of the equivalent stiffness (Figure 4(c)). As the roof pressure continues to rise to Pn, the support enters into the overflow segment. At this point, the retracting speed of the support is determined by the performance of the relief valve. Displacement formed in this stage is irreversible; therefore, it can be referred to as the “plastic deformation segment.” The final segment (stage e) is the release segment. The roof loading applied to the hydraulic support decreased rapidly and the leg retracted quickly.The ground response curve (GRC) theory [19, 21–24] suggests that the factor that determines the suitability of the shield support is the stiffness of the support. It holds the view that if the GRC and support reaction line intersect before stage d, the selected support is considered to be reasonable and effective. This type of stiffness analysis of the support was reported previously [25], but the plastic deformation characteristic of the support was neglected in the study.


## 2.2. Mathematical Model of an Equivalent Stiffness Doubly Telescopic Leg

Assuming that the compression deformation process is an isobaric expansion under the isentropic condition, the stiffness of a single-telescopic cylinder can be determined using the following equation:(1)Kz=11/Kl+1/Kp=1hl+he/Fp+Y/EhAh,where Kz is the equivalent stiffness of the single-telescopic leg, Kl is the equivalent stiffness formed by cylinder expansion and emulsion compression, Kp is the stiffness of the piston, hl is the compression length of the emulsion, he is the displacement of the rod caused by cylinder wall expansion, Fp is the axial force applied on the piston rod, Eh is the elastic modulus of the rod, and Ah is the cross-sectional area of the rod.h e can be calculated using the deformation theory of thick-walled cylinders (the leg is a typical high pressure thick wall structure). Combining the force equilibrium equations and geometric equations of the elements, the deformation equation of the cylinder can be expressed as(2)d2udr2+1rdudr-ur2=0,where u is the radial displacement of the cylinder and r is diameter of the cylinder. Substituting the boundary condition (see (3)) into (2) and ignoring the atmospheric pressure, the analytical solution for the radial displacement of the leg can be expressed as (4).(3)σr=d=p,σr=D=p2,(4)u=Δpdm2Eg,where m=1+μD2+(1-μ)d2/D2-d2, Δp=p-p2≈p, μ is Poisson’s ratio, d is the internal diameter of the cylinder, D is the external diameter of the cylinder, p is the cylinder pressure, p2 is the atmospheric pressure, and Eg is the elastic modulus of the cylinder. Then the expansion volume of the leg ΔVg can be expressed as (5) and he can be obtained using (6).(5)ΔVg=πd2yΔpm2Eg,(6)he=ΔVgA=2yΔpmEg,where A is the internal cross-sectional area of the cylinder. hl can be obtained using the compressible fluid theory (see (7)).(7)hl=4κV0Δpπd2=κΔpy,where κ is the compressibility of the emulsion and V0 is the initial volume of the emulsion.By combining (1), (6), and (7), the formula for calculating the equivalent stiffness of ith stage cylinder is given by(8)Kzi=ΔFzihi=πEgAh4κyi/di2Eg+8yimi/di2Ah+πY,where Kzi is the equivalent stiffness of ith stage cylinder, ΔFzi is the larger part of the load experienced by stage i, and hi is the displacement of the rod of stage i. If the volume force acting on the fluid is neglected, ΔFzi=ΔFp.For a doubly telescopic cylinder, it can be regarded as a series spring system (Figure2). Then its stiffness can be calculated as [18](9)K=Kz1Kz2Kz1+Kz2,where K is the equivalent stiffness of doubly telescopic leg.


## 2.3. Case Study

A certain type of shield support was taken as an example in this section. The main parameters of the leg are shown in Table1 (Figure 2) and its equivalent stiffness was calculated using (9). It is assumed that Po′ is 34 MPa (the active setting force is about 3856 kN). The results obtained by using (9) are shown in Figure 5. Lines OA, AB, and CD represent the stiffness of the leg at stage a, stage b, and stage c, respectively.Table 1
The main parameters of the doubly telescopic leg.
Name of the project Top stage Bottom stage Mobile column Structure lengthY (mm) / / 2747 Initial emulsion lengthyi (mm) / 2113 2197 Quality(t) 2.1 1.9 1.276 External/internal diameterDi/di (mm) 625/530 500/380 355 Setting/yielding pressurePo/Pn (MPa) 37.5/47.6 34/92.6 / Modulus of elasticity (GPa) 196 (20°C) Poisson’s ratio 0.3Figure 5
Stiffness variation curve of the doubly telescopic leg.The above settings are substituted into (10) and the calculated results are compared with the actual parameters of the support.(10)FPO˙=FPO′+∫OBKzdx=8394kN,FPn˙=FPO′+∫ODKzdx=10683kN,where FPO· is the calculated value of the setting force and FPn· is the calculated value of the yielding resistance. The error estimate is shown in (11). Given the fact that the errors are less than 5%, the calculation results of formula (9) are deemed to be of high accuracy.(11)δFPO=FPO˙-FPOFPO=1.462%,δFPn=FPn˙-FPnFPn=1.743%,where FPO is the rated setting force (8273 kN) and FPn is the rated working resistance (10500 kN).


## 3. Numerical Multibody Dynamics Model of the Shield Support

### 3.1. Preparation of the Simulation Model

The selected prototype of the support is shown in Figure6(a). Its operating height ranges from 3.8 m to 8.2 m and its working resistance reaches up to 21000 kN. The canopy, goaf shield, front bar, and rear bar are modelled as flexibles using HyperMesh (Figure 6(b)). The data are transferred to ADAMS in MNF format. The base is fixed on the ground as a rigid body; the legs are replicated using the equivalent stiffness model described in Section 2.3. Given the fact that the equilibrium jack is single extension with asymmetric tension and compression characteristics, it is modelled as a constant stiffness spring.Figure 6
Model of the shield support.
(a)
Prototype model of the support (b)
Grid model of the key parts of the support


### 3.2. The Simulation Experiment of the Support under Normal Working Conditions

The aim of the simulation test is to preliminarily verify whether the simulation model could satisfy the precision requirements for the subsequent tests. In order to simulate the stress state of the support more realistically, a simulated roof is placed on the canopy. The simulated roof can rotate freely in the horizontal plane and it connects with the roof at impact (Figure7). Applying a vertical downward static load of 20000 kN to the simulated roof, the response of the leg system is shown in Figure 8. As shown in the figure, the contact force increases to 18504 kN in 0.7 s (part of the load is bore by the equilibrium jack) and then remains stable. The load on the single leg is approximately half of the contact force (8421 kN). As mentioned previously, the elastic deformation process of the support can be divided into three segments (the lifting segment and the releasing segment are not included). Stage a is the active setting. In this portion, the roof load applied on the support is smaller than the active setting force (3856 kN), so essentially the support does not lower (due to the extra gravity of the simulated roof, the support lowers about 0.6 mm). In stage b, the roof load increases gradually beyond the setting force of the leg and the displacement of the roof decreased linearly. When the working resistance of the leg reaches 8421 kN, the rapid increase segment begins. As can be seen from Figure 8, the slope of the displacement and force change curves are clearly altered. These changes can be attributed to the change in the stiffness of the support. The overall displacement of the leg is h′ = 7071.18-6995.45 = 75.73 mm and the calculated result from Section 2.3 is h = 76 mm. Since h≈h′, the model is of high feasibility and reliability.Figure 7
Numerical model with the simulated roof.Figure 8
Dynamic response curves of the leg system.In order to facilitate the placement of the impact load, another model is tested. The same load is applied on the canopy above the leg (socket) without the simulated roof (Figure9). The results are the same. Therefore, the static pressure of the simulated roof can be replicated by applying the concentrated force on the socket.Figure 9
Numerical model with the concentrated force on the socket.


## 3.1. Preparation of the Simulation Model

The selected prototype of the support is shown in Figure6(a). Its operating height ranges from 3.8 m to 8.2 m and its working resistance reaches up to 21000 kN. The canopy, goaf shield, front bar, and rear bar are modelled as flexibles using HyperMesh (Figure 6(b)). The data are transferred to ADAMS in MNF format. The base is fixed on the ground as a rigid body; the legs are replicated using the equivalent stiffness model described in Section 2.3. Given the fact that the equilibrium jack is single extension with asymmetric tension and compression characteristics, it is modelled as a constant stiffness spring.Figure 6
Model of the shield support.
(a)
Prototype model of the support (b)
Grid model of the key parts of the support


## 3.2. The Simulation Experiment of the Support under Normal Working Conditions

The aim of the simulation test is to preliminarily verify whether the simulation model could satisfy the precision requirements for the subsequent tests. In order to simulate the stress state of the support more realistically, a simulated roof is placed on the canopy. The simulated roof can rotate freely in the horizontal plane and it connects with the roof at impact (Figure7). Applying a vertical downward static load of 20000 kN to the simulated roof, the response of the leg system is shown in Figure 8. As shown in the figure, the contact force increases to 18504 kN in 0.7 s (part of the load is bore by the equilibrium jack) and then remains stable. The load on the single leg is approximately half of the contact force (8421 kN). As mentioned previously, the elastic deformation process of the support can be divided into three segments (the lifting segment and the releasing segment are not included). Stage a is the active setting. In this portion, the roof load applied on the support is smaller than the active setting force (3856 kN), so essentially the support does not lower (due to the extra gravity of the simulated roof, the support lowers about 0.6 mm). In stage b, the roof load increases gradually beyond the setting force of the leg and the displacement of the roof decreased linearly. When the working resistance of the leg reaches 8421 kN, the rapid increase segment begins. As can be seen from Figure 8, the slope of the displacement and force change curves are clearly altered. These changes can be attributed to the change in the stiffness of the support. The overall displacement of the leg is h′ = 7071.18-6995.45 = 75.73 mm and the calculated result from Section 2.3 is h = 76 mm. Since h≈h′, the model is of high feasibility and reliability.Figure 7
Numerical model with the simulated roof.Figure 8
Dynamic response curves of the leg system.In order to facilitate the placement of the impact load, another model is tested. The same load is applied on the canopy above the leg (socket) without the simulated roof (Figure9). The results are the same. Therefore, the static pressure of the simulated roof can be replicated by applying the concentrated force on the socket.Figure 9
Numerical model with the concentrated force on the socket.


## 4. Influence of Impact Loading Position on the Hinge Point Force

### 4.1. Analysis of the Position of the Impact Load

With the continuous advance of the working face, breaks in the immediate roof and the main roof form the periodic weighting. This, combined with the random rock bursts, forms the impact loading that acts on the support. Ideally, the canopy of the support is in full contact with the roof. Thus, the resultant force from the impact load is located near the leg action line. The force is transferred to the goaf shield and the base through the legs and the equilibrium. However, in reality, the underground conditions are complex and changeable and the shearer usually cuts the roof unevenly. Thus, contact between the canopy and the roof is usually sporadic and irregular. In this case, when impact loading occurs, poor contact conditions can lead to unexpected failure of the support. Therefore, it is of significant importance to analyze the force transfer characteristics of the hinge joints of the hydraulic support when the impact load is acting on different parts of the canopy [20].The width direction of the canopy is defined as theX axis and the length direction as the Y axis. The impact point is taken every 1000 mm (ΔL) along the length direction and every 400 mm (ΔW) along the width direction (Figure 10). A static load of 18500 kN is applied on the socket and an impact load of 1500 kN is applied at each impact point. In order to evaluate the results, scale transmissions are done using (12), namely, the impact coefficient and excitation coefficient.(12)Ix,yi=Fimax-FioFI,Lx,yi=Fimax-FisFis×100%,where I(x,y)i is the impact coefficient of the ith hinge point, while the impact load acts at (x,y); L(x,y)i is the excitation coefficient of the ith hinge point, while the impact load acts at (x,y) (if not specifically pointed out, Li is the average value of the five points in the same width direction); Fimax is the peak force response of the ith hinge point caused by the impact loading; Fio is the force on the ith hinge point caused by the static load; Fis is the stable force response of the ith hinge point caused by the impact loading; and FI is the impact force (1500 kN in this study).Figure 10
Location of the impact load.


### 4.2. Influence on the Hinge Point Force of the Equilibrium Jack

The impact coefficient(Ia) and the excitation coefficient (La) at the hinge point of the equilibrium jack are shown in Figure 11 as the position of the impact load changes. Due to the fact that the selected support has only one equilibrium jack, the variation of the equilibrium jack is strictly symmetrical, while the impact loading varies along the width direction of the support. When the position of impact load changes along the width of the canopy, the change of Ia (2%) is far less than that caused by the change in the length direction (201%). When the impact load moves backward, Ia decreases gradually and then reverses (except the leg dominant area). This means that the equilibrium jack gradually changes from a tension to a compression state to adapt to the change of the external load position. Besides, the support has a trend of back tilting (see Figure 12). In the leg dominant area, the impact coefficient is almost unchanged. At this time, the impact load is mainly borne by the leg; the equilibrium jack does not need to balance the additional torque caused by the impact load.Figure 11
The impact coefficient and excitation coefficient of the equilibrium jack.Figure 12
Rotation angle of the canopy.With a decrease in equilibrium jack force,La increases gradually. La reaches a maximum when the impact force acts on the leg action line (60%). It can also be seen that the jack force is close to its rated working resistance when the impact force acts at the front end of the canopy. When the impact force acts at the back end of the canopy, the rod side of the jack reaches the setting pressure of the relief valve (yield area). The results show that the self-stability and carrying capacity of the support are very poor when the load acts on both ends of the support under this single region loading condition.


### 4.3. Influence on the Hinge Point Force of the Lemniscate Bars

The impact coefficients of the lemniscate bars at the hinge point are shown in Figures13(a) and 13(b) as the position of the impact load changes. When the impact load moves backward, the force of the front bar (Ic) is basically unchanged in the beginning, while the force of the rear bar (Id) decreased gradually. When the impact load passes through the leg action line, Ic and Id increase rapidly. Along the width direction, as the load moves rightward (towards the monitoring point), Ic increases gradually, while Id decreases initially but then increases slowly.Figure 13
The impact coefficients of the lemniscate bars.
(a)
The front bar (b)
The rear barSince the steady-state forces of the lemniscate bars have opposite trends, a new coefficient called “steady-state force impact coefficient (ISi)” is introduced. ISi can be calculated using (13). The steady-state force impact coefficients of the bars are shown in Figures 14(a) and 14(b). When the impact load moves backward, the steady-state force of the front (ISc) and rear (ISd) bar hinge joints increases gradually. Along the width direction, as the load moves rightward (toward the monitoring point), ISc increases slowly, while ISd decreases.(13)ISi=Fis-FioFio.The excitation coefficients of the lemniscate bars are shown in Figure 15 (average value across the width direction). As can be seen, the excitation force (25%~30%) of the lemniscate bars caused by impact loading acting on the front end of the canopy is much larger than that of the force acting on the rear end (1%~3%).Figure 14
The steady-state force impact coefficients of the lemniscate bars.
(a)
The front bar (b)
The rear barFigure 15
The excitation coefficients of the lemniscate bars.


### 4.4. Influence on the Hinge Point Force between the Canopy and the Goaf  Shield

The impact coefficient(Ib) and the excitation coefficient (Lb) of the hinge point force between the canopy and goaf shield are shown in Figure 16. When the impact load moves backward and rightward, Ib initially decreases but then rises with the leg action line as the boundary (similar to Id, but the impact effect is more severe). The excitation effect (Lb) at both ends of the canopy is about 10% but is only 0.4%~3% near the leg action line. Figure 16 indicates that the impact loading has little effect on Lb.Figure 16
The impact coefficient and the excitation coefficient of the hinge point force between the canopy and goaf shield.


## 4.1. Analysis of the Position of the Impact Load

With the continuous advance of the working face, breaks in the immediate roof and the main roof form the periodic weighting. This, combined with the random rock bursts, forms the impact loading that acts on the support. Ideally, the canopy of the support is in full contact with the roof. Thus, the resultant force from the impact load is located near the leg action line. The force is transferred to the goaf shield and the base through the legs and the equilibrium. However, in reality, the underground conditions are complex and changeable and the shearer usually cuts the roof unevenly. Thus, contact between the canopy and the roof is usually sporadic and irregular. In this case, when impact loading occurs, poor contact conditions can lead to unexpected failure of the support. Therefore, it is of significant importance to analyze the force transfer characteristics of the hinge joints of the hydraulic support when the impact load is acting on different parts of the canopy [20].The width direction of the canopy is defined as theX axis and the length direction as the Y axis. The impact point is taken every 1000 mm (ΔL) along the length direction and every 400 mm (ΔW) along the width direction (Figure 10). A static load of 18500 kN is applied on the socket and an impact load of 1500 kN is applied at each impact point. In order to evaluate the results, scale transmissions are done using (12), namely, the impact coefficient and excitation coefficient.(12)Ix,yi=Fimax-FioFI,Lx,yi=Fimax-FisFis×100%,where I(x,y)i is the impact coefficient of the ith hinge point, while the impact load acts at (x,y); L(x,y)i is the excitation coefficient of the ith hinge point, while the impact load acts at (x,y) (if not specifically pointed out, Li is the average value of the five points in the same width direction); Fimax is the peak force response of the ith hinge point caused by the impact loading; Fio is the force on the ith hinge point caused by the static load; Fis is the stable force response of the ith hinge point caused by the impact loading; and FI is the impact force (1500 kN in this study).Figure 10
Location of the impact load.


## 4.2. Influence on the Hinge Point Force of the Equilibrium Jack

The impact coefficient(Ia) and the excitation coefficient (La) at the hinge point of the equilibrium jack are shown in Figure 11 as the position of the impact load changes. Due to the fact that the selected support has only one equilibrium jack, the variation of the equilibrium jack is strictly symmetrical, while the impact loading varies along the width direction of the support. When the position of impact load changes along the width of the canopy, the change of Ia (2%) is far less than that caused by the change in the length direction (201%). When the impact load moves backward, Ia decreases gradually and then reverses (except the leg dominant area). This means that the equilibrium jack gradually changes from a tension to a compression state to adapt to the change of the external load position. Besides, the support has a trend of back tilting (see Figure 12). In the leg dominant area, the impact coefficient is almost unchanged. At this time, the impact load is mainly borne by the leg; the equilibrium jack does not need to balance the additional torque caused by the impact load.Figure 11
The impact coefficient and excitation coefficient of the equilibrium jack.Figure 12
Rotation angle of the canopy.With a decrease in equilibrium jack force,La increases gradually. La reaches a maximum when the impact force acts on the leg action line (60%). It can also be seen that the jack force is close to its rated working resistance when the impact force acts at the front end of the canopy. When the impact force acts at the back end of the canopy, the rod side of the jack reaches the setting pressure of the relief valve (yield area). The results show that the self-stability and carrying capacity of the support are very poor when the load acts on both ends of the support under this single region loading condition.


## 4.3. Influence on the Hinge Point Force of the Lemniscate Bars

The impact coefficients of the lemniscate bars at the hinge point are shown in Figures13(a) and 13(b) as the position of the impact load changes. When the impact load moves backward, the force of the front bar (Ic) is basically unchanged in the beginning, while the force of the rear bar (Id) decreased gradually. When the impact load passes through the leg action line, Ic and Id increase rapidly. Along the width direction, as the load moves rightward (towards the monitoring point), Ic increases gradually, while Id decreases initially but then increases slowly.Figure 13
The impact coefficients of the lemniscate bars.
(a)
The front bar (b)
The rear barSince the steady-state forces of the lemniscate bars have opposite trends, a new coefficient called “steady-state force impact coefficient (ISi)” is introduced. ISi can be calculated using (13). The steady-state force impact coefficients of the bars are shown in Figures 14(a) and 14(b). When the impact load moves backward, the steady-state force of the front (ISc) and rear (ISd) bar hinge joints increases gradually. Along the width direction, as the load moves rightward (toward the monitoring point), ISc increases slowly, while ISd decreases.(13)ISi=Fis-FioFio.The excitation coefficients of the lemniscate bars are shown in Figure 15 (average value across the width direction). As can be seen, the excitation force (25%~30%) of the lemniscate bars caused by impact loading acting on the front end of the canopy is much larger than that of the force acting on the rear end (1%~3%).Figure 14
The steady-state force impact coefficients of the lemniscate bars.
(a)
The front bar (b)
The rear barFigure 15
The excitation coefficients of the lemniscate bars.


## 4.4. Influence on the Hinge Point Force between the Canopy and the Goaf  Shield

The impact coefficient(Ib) and the excitation coefficient (Lb) of the hinge point force between the canopy and goaf shield are shown in Figure 16. When the impact load moves backward and rightward, Ib initially decreases but then rises with the leg action line as the boundary (similar to Id, but the impact effect is more severe). The excitation effect (Lb) at both ends of the canopy is about 10% but is only 0.4%~3% near the leg action line. Figure 16 indicates that the impact loading has little effect on Lb.Figure 16
The impact coefficient and the excitation coefficient of the hinge point force between the canopy and goaf shield.


## 5. Conclusions

( 1 ) From the point of view of stiffness, the working process of a shield support is analyzed and is redefined in five stages, namely, active setting, passive setting, rapid increase, overflow, and unloading. The rationality of this definition is discussed by means of constructing the equivalent stiffness model of the doubly telescopic leg and the numerical simulation of the support loading.( 2 ) The position of the impact load has different influence on the hinge point force of the shield support. In general, due to the large length-width ratio (270%), the impact effect and the excitation effect are more obvious when the impact load moves along the length direction.( 3 ) The impact effect, Ia and Ib, is most significant when the position of the impact load changes (400%). The excitation effect, Lc and Ld, is most significant when the position of the impact load changes (25%~30). The additional moment produced by the impact load is minimised when it is acting on the leg action line. At this point, both the impact effect and the excitation effect reach their lowest (0.1%~3%), which means that the optimal equilibrium region of the support is located near the leg action line.( 4 ) There is a positive correlation between the working resistance of the equilibrium jack and the position of the impact load in the length direction. When the impact load moves backward, the working resistance of the equilibrium jack decreases gradually and is eventually reversed. The equilibrium capacity of the equilibrium jack is poor when the impact load is applied behind the leg action line. The bearing capacity of the support can be improved by increasing the front-to-rear ratio of the canopy within a reasonable range.



---

*Source: 3854684-2018-03-15.xml*