# Research on Stress Analysis and Control of Surrounding Rock of Trapezoidal Roadways Based on Complex Variable Theory

Authors: Zhiqiang Wang; Chao Wu; Jianqiao Luo; Wenyu Lv; Lei Shi; Peng Wang
Journal: Shock and Vibration
 (2021)

---

## Abstract

Aiming at the problem of the serious deformation of the mining roadways in the trapezoidal section of the coal mine, the method of combining theoretical analysis and field application is used to study the surrounding rock control method of the trapezoidal roadways. The conformal mapping function of the trapezoidal roadways is calculated by the theory of complex change, and then from the analytical solution of the tangential stress distributed in the surrounding rock of trapezoidal roadways which is under the influence of the bidirectional unequal pressure, homogeneous, isotropic, and elastic rock mass is obtained. Research studies show that the roof-stress distribution of the trapezoidal roadways is uniform and the confining pressure is small, while the two sidewalls and the floor are opposite. The stress distribution of the two sidewalls and the floor varies greatly, and the stress concentration factor is large. The top corner of the trapezoidal roadways is basically not affected by stress concentration, but the stress concentration coefficient at the bottom corner is relatively large, and reinforcement measures are required in the roadway support. Based on the aforementioned research results, the multisupport scheme of “bolting with wire mesh and anchor cable + W-type steel belt + joist steel shed support + anchor cable grouting” was proposed to the surrounding rock of trapezoidal roadways with large stress caused by mining influence, thus solving the actual mining problem.


---

## Body

## 1. Introduction

With the continuous increase of coal mining depth, the difficulty of roadway support technology is rising. Due to high ground pressure and complex ground stress, the surrounding rock of some roadways has become broken and cannot bear the mining pressure well, resulting in the continuous increase of the deformation of the roadway. Conventional support methods can hardly guarantee the expected support effect, and the roadway must be repaired with reinforcement or structure optimization [1–3]. By apprehending the stress distribution of the surrounding rock of the roadway, it can effectively help us to choose the shape, size, and support measures of the roadway more reasonably.Trapezoidal roadways are widely distributed among coal mines in China. Such roadways are usually under a high confining pressure caused by complex stresses, with the uneven stress distribution. At present, approximate solutions are mostly used for nonaxisymmetric elastoplastic problems. Many scholars have conducted a lot of research on this method and obtained analytical solutions of the stress field and the plastic zone around the chamber under the condition of bidirectional unequal pressure [4–10]. References [11, 12] obtained the elastic-plastic perturbation solutions of the circular roadways under the action of bidirectional unequal pressure by using the perturbation method. Reference [13] based on the unified strength theory derived the boundary equation of the plastic zone of the roadway surrounding rock with the nonuniform stress field. Ma’s team [14–16] proposed the “butterfly shape” failure theory of circular roadways under the condition of bidirectional unequal pressure and gave the judgment criterion. At the same time, the concept of the plastic zone shape coefficient was put forward, and the calculation formula of the plastic zone morphological coefficient was derived, which can be used to distinguish the morphological characteristics of the plastic zone.The analytical method of complex function is the most superior in the plane problem of the orifice. Its characteristic is applying conformal transformation and the mapping function to transform the area of the complex orifice’s (e.g., the rectangular, trapezoidal, arch, and other noncircular curves that are often used in the roadway section)Z plane into the unit circle on the ζ plane. Namely, the orifice boundary is transformed into the unit circle boundary, which is conducive to simplifying the boundary conditions, thereby deducing the elastic analytical formula for the stress of the surrounding rock of the orifice, which has become a solution to the common method of different section shape chamber mechanics problems. Reference [17] deduced the analytic solution of surrounding rock stress of rectangular roadways based on complex variation theory and analyzed the distribution regularities of tangential stress in the surrounding rock of roadways with different lateral pressure coefficients. Reference [18] also derived the analytical solution of the surrounding rock stress of a circular roadway based on the theory of complex variation and combined it with Mohr-Coulomb criterion to derive the formula of broken range of surrounding rock plastic zone.In summary, there are kinds of researches on the nonaxisymmetric elastic-plastic problems of circular, rectangular, and arched roadways. In view of the engineering background where most coal mine roadway still uses trapezoidal roadways as mining roadway and serious deformation, this paper is based on the theory of complex variation theory to solve the stress field of surrounding rock of trapezoidal roadways under the condition of bidirectional unequal pressure, and the control technology of the surrounding rock of trapezoidal roadways is studied based on this theoretical basis.


## 2. Engineering Overview and Deformation Characteristics of Roadway Surrounding Rock

The Laogongyingzi Mine is located in Chifeng City, Inner Mongolia Autonomous Region, China. The roadways section of the coal mine working face is a trapezoid with a top length of 3 m, a bottom length of 4.6 m, and a vertical height of 3 m. The buried depth is 400 m, the average thickness of the mined coal seam is about 4.2 m, and 0～0.5 m false roofs will fall along with mining. The direct roof is dominated by mudstone and fine sandstone, with an average thickness of about 20 m, while the floor is dominated by fine sandstone. During the mining, the working face advances at a speed of 7 m/d, and the mining intensity is severe.The original support scheme of the roadway adopts the combined support of “cable between shelf and shed + steel belt.” The roof and the two sidewalls are, respectively, laid with 6 rows of I-beam with the length of 800 mm, and the spacing is 500 mm, and then the back plate is laid. After that, the two sidewalls are supported by anchor cables. The anchor cables are made of Φ22 mm × 4000 mm high-strength prestressed low-relaxation steel strands. Three anchor cables are laid in each row with a distance of 1000 mm. The first anchor cable is 300 mm away from the roof. When supporting the second anchor cable and the third anchor cable, a 150 mm × 150 mm × 10 mm pallet is used to match the I-beam.During the mining of the working face, the surrounding rock of the mining roadway is unstable and deforms greatly, and the floor needs to be cleaned and repaired many times, which seriously restricts the normal mining efficiency of the working face, as shown in Figure1. Its instability deformation characteristics are mainly as follows:(1)
The large deformation of the coal pillar sidewall shows the characteristics of sustainability and large range. Severe bulging and collapse of the coal pillar sidewall of the roadway happened; along the direction of the roadway, there is a range of more than 15 m, and the bulging amount is about 1000 mm. The bulging causes the I-beam to undergo obvious bending and deformation. A large area of the coal pillar sidewall collapsed, resulting in a large-scale cavity between the I-beam and the coal pillar sidewall, and the depth of part of the coal pillar sidewall was about 1500 mm, resulting in failure of the bolt and I-beam support. The coal mine lays wooden boards between the I-beam and the broken coal body to increase the force-bearing area in order to achieve the purpose of controlling the deformation of the roadside. However, the effect was not good, the wooden boards were broken, and the bulging and collapse continued.(2)
The degree of deformation of solid coal sidewall is lower than that of coal pillar sidewall. During roadway excavation, the amount of deformation is small. When it is affected by the mining of the working face, the roadway deforms seriously. The displacement of the bulge is 500～750 mm. Part of the I-beam is bent and deformed, and the damage range along the axial direction of the roadway is about 3 m.(3)
The roadway roof subsidence is small, the subsidence is about 200 mm, and the bending deformation of the steel strip is obvious, especially when it is affected by the mining of the working face, and part of the roof leaks. The floor heave of the roadway is obvious, and the measured maximum floor heave volume reaches more than 1000 mm, presenting a state of extrusion failure which is high in the middle and low at both ends, as well as even tensile failure and forming cracks.Figure 1
The mine pressure of the trapezoidal roadway. (a) The sidewalls of coal pillar heave for 1500 mm. The I-beam bent. (b) The roof declined 500 mm. The floor heaves seriously. The I-beam bent.
(a)(b)


## 3. Solution of Complex Variable Function of Stress in Surrounding Rock of Trapezoidal Roadways

### 3.1. Steps to Solve Complex Variable Function of Roadways

From the polar coordinate stress transformation formula of elastic mechanics, the complex variable function expression of the stress component in theζ plane polar coordinate is obtained according to the complex variable theory [19–21]:(1)σθ+σρ=2Φζ+Φζ¯=4ReΦζ,σθ−σρ+2iτρθ=2ζ2ρ2ω′ζ¯ωζ¯Φ′ζ+ω′ζΨζ,(2)Φζ=φ′ζ=φ′ωζ=φ′ζω′ζ,Ψζ=ψ′ζ=ψ′ωζ=ψ′ζω′ζ,where σθ, σρ, and τρθ are stress components of elastic body in the curvilinear coordinate system; φζ, ψζ, Φζ, and Ψζ are the complex potential analytic functions of complex variables ζ; and ωζ is a conformal mapping function.From elasticity theory and complex variable function theory, we obtained the following equation:(3)φζ=1+μ8πFx¯+iFy¯lnζ+Bωζ+φ0ζ,ψζ=3−μ8πFx¯−iFy¯lnζ+B′+iC′ωζ+ψ0ζ.Consideringσ=ωζ=eiθ, σ is the value of complex variable ζ at roadway boundary. Then, we get the boundary conditions expressed by φ0ζ and ψ0ζ, which can be calculated by the following Cauchy integral equations:(4)φ0ζ+12πi∫σωσω′σ¯φ0′σ¯σ−ζdσ=12πi∫σf0σ−ζdσ,(5)ψ0ζ+12πi∫σωσ¯ω′σ¯φ0′σσ−ζdσ=12πi∫σf0¯σ−ζdσ,where f0 is the known function for σ, and its expression is(6)f0=i∫Fx¯+iFy¯ds−Fx¯−iFy¯2πlnσ−1+μ8πFx¯−iFy¯×ωσω′σ¯σ−2Bωσ−B′−iC′ωσ¯,where μ is the elastic body shear modulus and Fx¯ and Fy¯ are the sum of surface force components fx and fy at orifice which are already known. B, B′, and C′ are related to the stress which is far from the orifice. The expressions are(7)B=14σx∞+σy∞,B′=12σy∞−σx∞,C′=τxy∞,(8)φ0ζ=∑k=1nαkζk,(9)ψ0ζ=∑k=1nβkζk.The following are the general steps of the complex function method to solve the orifice problems.


#### 3.1.1. Conformal Transformation

Transform the area occupied by the elastic body on theZ plane into the area on the ζ plane, and the general form of the trapezoidal roadway mapping function is [22](10)z=ωζ=Rζ+c0+c1ζ−1+c2ζ−2+⋯+cnζ−n=Rζ+∑k=0∞ckζ−k,where R is a real number, which is related to the orifice size, and ck is a complex number. The trapezoidal roadway discussed in this paper has an axis of symmetry, while ck in equation (10) must be a real number. The greater the k series, the higher the accuracy, and generally k = 3.


#### 3.1.2. Solving for Stress Components

Mapping formula (10) for different shapes of orifices can be found from books on complex functions or established by the theory of complex functions. Substitute functions ω(ζ) and equations (6)–(9) into equations (4) and (5), respectively, to get the equations of φ0ζ and ψ0ζ by applying Cauchy integral formula. Then calculate the equations of φζ and ψζ, which can be obtained by equation (3). Substitute the equations of φζ and ψζ into equation (2) to find the equations of Φζ and Ψζ. Finally, substitute them into equation (1) to obtain the stress components σθ, σρ, and τρθ in the curve coordinates.


### 3.2. Analytical Solution of Tangential Stress of Surrounding Rock around Trapezoidal Roadways

#### 3.2.1. Conformal Mapping of Trapezoidal Roadways

The mechanical model of the trapezoidal roadway is established as shown in Figure2. When the roadway has not been affected by mining stress, the stress distribution around the roadway is not affected by the leading support stress and lateral support stress, and it is subjected to bidirectional unequal compressive stresses; and the horizontal and vertical original rock stresses are evenly distributed. The length of the roadway in the z-axis direction is much larger than the length of the other two directions. The stress state is not considered. Therefore, the stress on any section of the roadway is simplified to an axisymmetric plane stress problem. The trapezoidal roadway on the Z plane is mapped to the unit circle on the ζ plane through the conformal transformation method of the complex theory, as shown in Figure 3.Figure 2
Trapezoidal roadways geometry on theZ plane.Figure 3
Unit circle of conformal mapping on theζ plane.In the analysis of the left half of the symmetrical trapezoidal roadway, to facilitate calculation, divide the left half into 3 parts. Let the polar coordinates of any point around the trapezoidal roadway be (ri, αi). Let the polar coordinate of the corresponding mapping point on the unit circle be (1, βi). Because the periphery of the roadway corresponds to the unit circle circumference on the plane, there are(11)rieiαi=Reiβi+∑k=0∞cke−ikβi,(12)ri=Reiβi−αi+∑k=0∞cke−ikβi+αi.Separate the imaginary and real parts of equation (12) to obtain (13)sinβi−αi+∑k=0∞cksinkβi+αi=0,(14)ri=Rcosβi−αi+∑k=0∞ckcoskβi+αi.According to the mathematical expression of the mapping function, now findR and ck in the mapping function.Select a rectangular coordinate system in theZ plane, and set the x-axis as the symmetry axis of the roadway. On the periphery of the roadway on the left of the x-axis, take m + 1 points in the counterclockwise direction. The rectangular coordinate value and the corresponding polar coordinate value of the i-th point are (xi, yi) and (ri, αi), respectively. According to the geometric relationship of the trapezoidal roadway in Figure 2, the following equation can be found:(15)r1=ccosα0≤α≤arctanbc,r2=x22−x12+x32+2dx32x1cosθ3−θ31−x2cosθ31arctanbc≤α≤π−arctanac,0≤θ31≤α−arctanbc,r3=ccosπ−απ−arctanac≤α≤π.Among them,(16)x1=a2+c2,x2=b2+c2,x3=x22+x12−2x1x2cosα−arctanbc,d=x22−x12−x32+2x1x3cosθ22x1cosθ2+2x2cosθ1−2x3,θ1=π2−θ+arctanbc,θ2=π2+θ−α,θ31=α−arctanbc,θ=arctan2ca−b.Leti = 1 correspond to the polar coordinate point (r1, 0). Let the polar coordinate of the corresponding point on the circumference of the plane unit circle ζ be (1, 0). Then, R1 can be obtained from equation (14) and α1 = β1 = 0, where(17a)R1=r11+∑k=0∞ck.Leti = 3 correspond to the polar coordinate point (r2, π/2). Let the polar coordinate of the corresponding point on the circumference of the plane unit circle ζ be (1, π/2). Then, R2 can be obtained from formula (14) and α2 = β2 = π/2, where(17b)R2=r21+∑k=0∞ckcosk+1/2π.Leti = 5 correspond to the polar coordinate point (r3, π). Let the polar coordinate of the corresponding point on the circumference of the plane unit circle ζ be (1, π). Then, R3 can be obtained from equation (14) and α3 = β3 = π, where(17c)R3=r31+∑k=0∞ckcosk+1π.The three formulas above give the relationship betweenR and ck. From equation (14), when k = 0, 1, 2, …, n, n takes a finite positive integer, and equation (14) cannot be accurately satisfied, but ck (k = 0, 1, 2, …, n) can be determined by optimization techniques. Meanwhile, the objective function can be taken as(18)f=∑i=1m+1ri−Ricosβi−αi+∑k=0nckcoskβi+αi2.In equation (18), ri and αi are known values, ck is an unknown value, and βi is also an unknown value, but it can be connected with ck and αi through equation (13).Whenm > n, ck∗ that minimizes f can be obtained by equation (18). ck∗ is the mapping function coefficient ck which is to be determined, and now fmin > 0.According to the actual parameters of the trapezoidal roadways in one coal mine, leta = 2.3, b = 1.5, c = 1.5, n = 3, and m = 4. In order to facilitate the calculation, take 5 special points of the trapezoidal roadway, and their polar coordinates are (c, 0), (b2+c2, arctanb/c), (a+b/2, π/2), (a2+c2, π−arctana/c), (c, π). Substitute them into equations (8) and (9); ck is calculated by computer as c0 = −0.1011, c1 = −0.1101, c2 = 0.0859, c3 = 0.0014, and R = 1.69, and the mapping function of the trapezoidal section is determined as(19)z=ωζ=1.69×ζ−0.0896−0.1101ζ+0.0859ζ2+0.0014ζ3.


#### 3.2.2. Analytical Formula for Tangential Stress of Surrounding Rock around Trapezoidal Roadway

Now perform calculation according to the steps of the complex variable function solution method; the following equation is derived by equation (6):(20)B=14P+λP,B′=12λP−P,C′=0.From equation (11), the following equation is derived:(21)ωσ=Rσ+c0+c11σ+c21σ2+c31σ3,ωσ¯=R1σ+c0+c1σ+c2σ2+c3σ3,ω′σ=R1−c11σ2−c22σ3−c33σ4,ω′σ¯=R1−c1σ2−2c2σ3−3c3σ4.Introducing equation (21) into equation (9) gives(22)f0=−2BRσ+c0+c11σ+c21σ2+c31σ3−B′R1σ+c0+c1σ+c2σ2+c3σ3.Introducing equation (22) into equation (4) to solve with Harnack’s theorem, the integral on the right is(23)12πi∫σf0σ−ζdσ=−2BRζ−B′Rc1ζ−B′Rc2ζ2−B′Rc3ζ3.Calculate by combining equations (4) and (13); then the solutions from α1 to αi are found:(24)α1=−BR,α2=−0.2746B′R,α3=−0.0096B′R,α4=α5=⋅⋅⋅=αn=0.Introduce them into formula (7) to get(25)φ0ζ=−BRζ−0.2746B′Rζ2−0.0096B′Rζ3.Substituting equations (10) and (25) into equation (3) gives(26)φζ=BRζ+c0+c1ζ−1+c2ζ−2+⋯+cnζ−n−BRζ−0.2746B′Rζ2−0.0096B′Rζ3.From equation (2),(27)Φζ=φ′ζω′ζ=B−c1ζ2−2c2ζ−3c3ζ4−c1ζ2−2c2ζ−3c3−B′0.5492ζ5+0.0288ζ7ζ4−c1ζ2−2c2ζ−3c3.Letρ of the unit circle which has already been mapped be 1; then ζ=ρeiθ=eiθ=cosθ+isinθ. Substitute ζ into equation (27). After the calculation of the first half of equation (27), we have the following result:(28)B−c1ζ2−2c2ζ−3c3ζ4−c1ζ2−2c2ζ−3c3=−Bc1cos2θ+2c2cosθ+3c3+isin2θ+2c2sinθcos4θ−c1cos2θ−2c2cosθ−3c3+isin4θ−c1sin2θ−2c2sinθ,(29)Z1=x1,y1=c1cos2θ+2c2cosθ+3c3+isin2θ+2c2sin2θ,Z2=x2,y2=cos4θ−c1cos2θ−2c2cosθ−3c3+isin4θ−c1sin2θ−2c2sinθ.Then,(30)ReZ1Z2=x1x2+y1y2x22+y22.From equation (30),(31)x1=c1cos2θ+2c2cosθ+3c3,x2=cos4θ−c1cos2θ−2c2cosθ−3c3,y1=sin2θ+2c2sinθ,y2=sin4θ−c1sin2θ−2c2sinθ.The calculation of the last half part of equation (27) gives(32)−B′0.5492ζ5+0.0288ζ7ζ4−c1ζ2−2c2ζ−3c3=−B′0.5492cos5θ+0.0288cos7θ+i0.5492sin5θ+0.0288sin7θcos4θ−c1cos2θ−2c2cosθ−3c3+isin4θ−c1sin2θ−2c2sinθ,(33)Z3=x3,y3=0.5492cos5θ+0.0288cos7θ+i0.5492sin5θ+0.0288sin7θ.Then,(34)ReZ3Z2=x2x3+y2y3x22+y22.From equation (34),(35)x3=0.5492cos5θ+0.0288cos7θ,y3=0.5492sin5θ+0.0288sin7θ.Substitute equations (28), (30), (32), and (34) into formula (1). Because the value of σρ around the roadway is 0, we have σθ=4ReΦζ. The analytical equation for the tangential stress around the trapezoidal roadways is calculated as(36)σθ=P+λPx1x2+y1y2x22+y22+2P−λPx2x3+y2y3x22+y22.


## 3.1. Steps to Solve Complex Variable Function of Roadways

From the polar coordinate stress transformation formula of elastic mechanics, the complex variable function expression of the stress component in theζ plane polar coordinate is obtained according to the complex variable theory [19–21]:(1)σθ+σρ=2Φζ+Φζ¯=4ReΦζ,σθ−σρ+2iτρθ=2ζ2ρ2ω′ζ¯ωζ¯Φ′ζ+ω′ζΨζ,(2)Φζ=φ′ζ=φ′ωζ=φ′ζω′ζ,Ψζ=ψ′ζ=ψ′ωζ=ψ′ζω′ζ,where σθ, σρ, and τρθ are stress components of elastic body in the curvilinear coordinate system; φζ, ψζ, Φζ, and Ψζ are the complex potential analytic functions of complex variables ζ; and ωζ is a conformal mapping function.From elasticity theory and complex variable function theory, we obtained the following equation:(3)φζ=1+μ8πFx¯+iFy¯lnζ+Bωζ+φ0ζ,ψζ=3−μ8πFx¯−iFy¯lnζ+B′+iC′ωζ+ψ0ζ.Consideringσ=ωζ=eiθ, σ is the value of complex variable ζ at roadway boundary. Then, we get the boundary conditions expressed by φ0ζ and ψ0ζ, which can be calculated by the following Cauchy integral equations:(4)φ0ζ+12πi∫σωσω′σ¯φ0′σ¯σ−ζdσ=12πi∫σf0σ−ζdσ,(5)ψ0ζ+12πi∫σωσ¯ω′σ¯φ0′σσ−ζdσ=12πi∫σf0¯σ−ζdσ,where f0 is the known function for σ, and its expression is(6)f0=i∫Fx¯+iFy¯ds−Fx¯−iFy¯2πlnσ−1+μ8πFx¯−iFy¯×ωσω′σ¯σ−2Bωσ−B′−iC′ωσ¯,where μ is the elastic body shear modulus and Fx¯ and Fy¯ are the sum of surface force components fx and fy at orifice which are already known. B, B′, and C′ are related to the stress which is far from the orifice. The expressions are(7)B=14σx∞+σy∞,B′=12σy∞−σx∞,C′=τxy∞,(8)φ0ζ=∑k=1nαkζk,(9)ψ0ζ=∑k=1nβkζk.The following are the general steps of the complex function method to solve the orifice problems.


### 3.1.1. Conformal Transformation

Transform the area occupied by the elastic body on theZ plane into the area on the ζ plane, and the general form of the trapezoidal roadway mapping function is [22](10)z=ωζ=Rζ+c0+c1ζ−1+c2ζ−2+⋯+cnζ−n=Rζ+∑k=0∞ckζ−k,where R is a real number, which is related to the orifice size, and ck is a complex number. The trapezoidal roadway discussed in this paper has an axis of symmetry, while ck in equation (10) must be a real number. The greater the k series, the higher the accuracy, and generally k = 3.


### 3.1.2. Solving for Stress Components

Mapping formula (10) for different shapes of orifices can be found from books on complex functions or established by the theory of complex functions. Substitute functions ω(ζ) and equations (6)–(9) into equations (4) and (5), respectively, to get the equations of φ0ζ and ψ0ζ by applying Cauchy integral formula. Then calculate the equations of φζ and ψζ, which can be obtained by equation (3). Substitute the equations of φζ and ψζ into equation (2) to find the equations of Φζ and Ψζ. Finally, substitute them into equation (1) to obtain the stress components σθ, σρ, and τρθ in the curve coordinates.


## 3.1.1. Conformal Transformation

Transform the area occupied by the elastic body on theZ plane into the area on the ζ plane, and the general form of the trapezoidal roadway mapping function is [22](10)z=ωζ=Rζ+c0+c1ζ−1+c2ζ−2+⋯+cnζ−n=Rζ+∑k=0∞ckζ−k,where R is a real number, which is related to the orifice size, and ck is a complex number. The trapezoidal roadway discussed in this paper has an axis of symmetry, while ck in equation (10) must be a real number. The greater the k series, the higher the accuracy, and generally k = 3.


## 3.1.2. Solving for Stress Components

Mapping formula (10) for different shapes of orifices can be found from books on complex functions or established by the theory of complex functions. Substitute functions ω(ζ) and equations (6)–(9) into equations (4) and (5), respectively, to get the equations of φ0ζ and ψ0ζ by applying Cauchy integral formula. Then calculate the equations of φζ and ψζ, which can be obtained by equation (3). Substitute the equations of φζ and ψζ into equation (2) to find the equations of Φζ and Ψζ. Finally, substitute them into equation (1) to obtain the stress components σθ, σρ, and τρθ in the curve coordinates.


## 3.2. Analytical Solution of Tangential Stress of Surrounding Rock around Trapezoidal Roadways

### 3.2.1. Conformal Mapping of Trapezoidal Roadways

The mechanical model of the trapezoidal roadway is established as shown in Figure2. When the roadway has not been affected by mining stress, the stress distribution around the roadway is not affected by the leading support stress and lateral support stress, and it is subjected to bidirectional unequal compressive stresses; and the horizontal and vertical original rock stresses are evenly distributed. The length of the roadway in the z-axis direction is much larger than the length of the other two directions. The stress state is not considered. Therefore, the stress on any section of the roadway is simplified to an axisymmetric plane stress problem. The trapezoidal roadway on the Z plane is mapped to the unit circle on the ζ plane through the conformal transformation method of the complex theory, as shown in Figure 3.Figure 2
Trapezoidal roadways geometry on theZ plane.Figure 3
Unit circle of conformal mapping on theζ plane.In the analysis of the left half of the symmetrical trapezoidal roadway, to facilitate calculation, divide the left half into 3 parts. Let the polar coordinates of any point around the trapezoidal roadway be (ri, αi). Let the polar coordinate of the corresponding mapping point on the unit circle be (1, βi). Because the periphery of the roadway corresponds to the unit circle circumference on the plane, there are(11)rieiαi=Reiβi+∑k=0∞cke−ikβi,(12)ri=Reiβi−αi+∑k=0∞cke−ikβi+αi.Separate the imaginary and real parts of equation (12) to obtain (13)sinβi−αi+∑k=0∞cksinkβi+αi=0,(14)ri=Rcosβi−αi+∑k=0∞ckcoskβi+αi.According to the mathematical expression of the mapping function, now findR and ck in the mapping function.Select a rectangular coordinate system in theZ plane, and set the x-axis as the symmetry axis of the roadway. On the periphery of the roadway on the left of the x-axis, take m + 1 points in the counterclockwise direction. The rectangular coordinate value and the corresponding polar coordinate value of the i-th point are (xi, yi) and (ri, αi), respectively. According to the geometric relationship of the trapezoidal roadway in Figure 2, the following equation can be found:(15)r1=ccosα0≤α≤arctanbc,r2=x22−x12+x32+2dx32x1cosθ3−θ31−x2cosθ31arctanbc≤α≤π−arctanac,0≤θ31≤α−arctanbc,r3=ccosπ−απ−arctanac≤α≤π.Among them,(16)x1=a2+c2,x2=b2+c2,x3=x22+x12−2x1x2cosα−arctanbc,d=x22−x12−x32+2x1x3cosθ22x1cosθ2+2x2cosθ1−2x3,θ1=π2−θ+arctanbc,θ2=π2+θ−α,θ31=α−arctanbc,θ=arctan2ca−b.Leti = 1 correspond to the polar coordinate point (r1, 0). Let the polar coordinate of the corresponding point on the circumference of the plane unit circle ζ be (1, 0). Then, R1 can be obtained from equation (14) and α1 = β1 = 0, where(17a)R1=r11+∑k=0∞ck.Leti = 3 correspond to the polar coordinate point (r2, π/2). Let the polar coordinate of the corresponding point on the circumference of the plane unit circle ζ be (1, π/2). Then, R2 can be obtained from formula (14) and α2 = β2 = π/2, where(17b)R2=r21+∑k=0∞ckcosk+1/2π.Leti = 5 correspond to the polar coordinate point (r3, π). Let the polar coordinate of the corresponding point on the circumference of the plane unit circle ζ be (1, π). Then, R3 can be obtained from equation (14) and α3 = β3 = π, where(17c)R3=r31+∑k=0∞ckcosk+1π.The three formulas above give the relationship betweenR and ck. From equation (14), when k = 0, 1, 2, …, n, n takes a finite positive integer, and equation (14) cannot be accurately satisfied, but ck (k = 0, 1, 2, …, n) can be determined by optimization techniques. Meanwhile, the objective function can be taken as(18)f=∑i=1m+1ri−Ricosβi−αi+∑k=0nckcoskβi+αi2.In equation (18), ri and αi are known values, ck is an unknown value, and βi is also an unknown value, but it can be connected with ck and αi through equation (13).Whenm > n, ck∗ that minimizes f can be obtained by equation (18). ck∗ is the mapping function coefficient ck which is to be determined, and now fmin > 0.According to the actual parameters of the trapezoidal roadways in one coal mine, leta = 2.3, b = 1.5, c = 1.5, n = 3, and m = 4. In order to facilitate the calculation, take 5 special points of the trapezoidal roadway, and their polar coordinates are (c, 0), (b2+c2, arctanb/c), (a+b/2, π/2), (a2+c2, π−arctana/c), (c, π). Substitute them into equations (8) and (9); ck is calculated by computer as c0 = −0.1011, c1 = −0.1101, c2 = 0.0859, c3 = 0.0014, and R = 1.69, and the mapping function of the trapezoidal section is determined as(19)z=ωζ=1.69×ζ−0.0896−0.1101ζ+0.0859ζ2+0.0014ζ3.


### 3.2.2. Analytical Formula for Tangential Stress of Surrounding Rock around Trapezoidal Roadway

Now perform calculation according to the steps of the complex variable function solution method; the following equation is derived by equation (6):(20)B=14P+λP,B′=12λP−P,C′=0.From equation (11), the following equation is derived:(21)ωσ=Rσ+c0+c11σ+c21σ2+c31σ3,ωσ¯=R1σ+c0+c1σ+c2σ2+c3σ3,ω′σ=R1−c11σ2−c22σ3−c33σ4,ω′σ¯=R1−c1σ2−2c2σ3−3c3σ4.Introducing equation (21) into equation (9) gives(22)f0=−2BRσ+c0+c11σ+c21σ2+c31σ3−B′R1σ+c0+c1σ+c2σ2+c3σ3.Introducing equation (22) into equation (4) to solve with Harnack’s theorem, the integral on the right is(23)12πi∫σf0σ−ζdσ=−2BRζ−B′Rc1ζ−B′Rc2ζ2−B′Rc3ζ3.Calculate by combining equations (4) and (13); then the solutions from α1 to αi are found:(24)α1=−BR,α2=−0.2746B′R,α3=−0.0096B′R,α4=α5=⋅⋅⋅=αn=0.Introduce them into formula (7) to get(25)φ0ζ=−BRζ−0.2746B′Rζ2−0.0096B′Rζ3.Substituting equations (10) and (25) into equation (3) gives(26)φζ=BRζ+c0+c1ζ−1+c2ζ−2+⋯+cnζ−n−BRζ−0.2746B′Rζ2−0.0096B′Rζ3.From equation (2),(27)Φζ=φ′ζω′ζ=B−c1ζ2−2c2ζ−3c3ζ4−c1ζ2−2c2ζ−3c3−B′0.5492ζ5+0.0288ζ7ζ4−c1ζ2−2c2ζ−3c3.Letρ of the unit circle which has already been mapped be 1; then ζ=ρeiθ=eiθ=cosθ+isinθ. Substitute ζ into equation (27). After the calculation of the first half of equation (27), we have the following result:(28)B−c1ζ2−2c2ζ−3c3ζ4−c1ζ2−2c2ζ−3c3=−Bc1cos2θ+2c2cosθ+3c3+isin2θ+2c2sinθcos4θ−c1cos2θ−2c2cosθ−3c3+isin4θ−c1sin2θ−2c2sinθ,(29)Z1=x1,y1=c1cos2θ+2c2cosθ+3c3+isin2θ+2c2sin2θ,Z2=x2,y2=cos4θ−c1cos2θ−2c2cosθ−3c3+isin4θ−c1sin2θ−2c2sinθ.Then,(30)ReZ1Z2=x1x2+y1y2x22+y22.From equation (30),(31)x1=c1cos2θ+2c2cosθ+3c3,x2=cos4θ−c1cos2θ−2c2cosθ−3c3,y1=sin2θ+2c2sinθ,y2=sin4θ−c1sin2θ−2c2sinθ.The calculation of the last half part of equation (27) gives(32)−B′0.5492ζ5+0.0288ζ7ζ4−c1ζ2−2c2ζ−3c3=−B′0.5492cos5θ+0.0288cos7θ+i0.5492sin5θ+0.0288sin7θcos4θ−c1cos2θ−2c2cosθ−3c3+isin4θ−c1sin2θ−2c2sinθ,(33)Z3=x3,y3=0.5492cos5θ+0.0288cos7θ+i0.5492sin5θ+0.0288sin7θ.Then,(34)ReZ3Z2=x2x3+y2y3x22+y22.From equation (34),(35)x3=0.5492cos5θ+0.0288cos7θ,y3=0.5492sin5θ+0.0288sin7θ.Substitute equations (28), (30), (32), and (34) into formula (1). Because the value of σρ around the roadway is 0, we have σθ=4ReΦζ. The analytical equation for the tangential stress around the trapezoidal roadways is calculated as(36)σθ=P+λPx1x2+y1y2x22+y22+2P−λPx2x3+y2y3x22+y22.


## 3.2.1. Conformal Mapping of Trapezoidal Roadways

The mechanical model of the trapezoidal roadway is established as shown in Figure2. When the roadway has not been affected by mining stress, the stress distribution around the roadway is not affected by the leading support stress and lateral support stress, and it is subjected to bidirectional unequal compressive stresses; and the horizontal and vertical original rock stresses are evenly distributed. The length of the roadway in the z-axis direction is much larger than the length of the other two directions. The stress state is not considered. Therefore, the stress on any section of the roadway is simplified to an axisymmetric plane stress problem. The trapezoidal roadway on the Z plane is mapped to the unit circle on the ζ plane through the conformal transformation method of the complex theory, as shown in Figure 3.Figure 2
Trapezoidal roadways geometry on theZ plane.Figure 3
Unit circle of conformal mapping on theζ plane.In the analysis of the left half of the symmetrical trapezoidal roadway, to facilitate calculation, divide the left half into 3 parts. Let the polar coordinates of any point around the trapezoidal roadway be (ri, αi). Let the polar coordinate of the corresponding mapping point on the unit circle be (1, βi). Because the periphery of the roadway corresponds to the unit circle circumference on the plane, there are(11)rieiαi=Reiβi+∑k=0∞cke−ikβi,(12)ri=Reiβi−αi+∑k=0∞cke−ikβi+αi.Separate the imaginary and real parts of equation (12) to obtain (13)sinβi−αi+∑k=0∞cksinkβi+αi=0,(14)ri=Rcosβi−αi+∑k=0∞ckcoskβi+αi.According to the mathematical expression of the mapping function, now findR and ck in the mapping function.Select a rectangular coordinate system in theZ plane, and set the x-axis as the symmetry axis of the roadway. On the periphery of the roadway on the left of the x-axis, take m + 1 points in the counterclockwise direction. The rectangular coordinate value and the corresponding polar coordinate value of the i-th point are (xi, yi) and (ri, αi), respectively. According to the geometric relationship of the trapezoidal roadway in Figure 2, the following equation can be found:(15)r1=ccosα0≤α≤arctanbc,r2=x22−x12+x32+2dx32x1cosθ3−θ31−x2cosθ31arctanbc≤α≤π−arctanac,0≤θ31≤α−arctanbc,r3=ccosπ−απ−arctanac≤α≤π.Among them,(16)x1=a2+c2,x2=b2+c2,x3=x22+x12−2x1x2cosα−arctanbc,d=x22−x12−x32+2x1x3cosθ22x1cosθ2+2x2cosθ1−2x3,θ1=π2−θ+arctanbc,θ2=π2+θ−α,θ31=α−arctanbc,θ=arctan2ca−b.Leti = 1 correspond to the polar coordinate point (r1, 0). Let the polar coordinate of the corresponding point on the circumference of the plane unit circle ζ be (1, 0). Then, R1 can be obtained from equation (14) and α1 = β1 = 0, where(17a)R1=r11+∑k=0∞ck.Leti = 3 correspond to the polar coordinate point (r2, π/2). Let the polar coordinate of the corresponding point on the circumference of the plane unit circle ζ be (1, π/2). Then, R2 can be obtained from formula (14) and α2 = β2 = π/2, where(17b)R2=r21+∑k=0∞ckcosk+1/2π.Leti = 5 correspond to the polar coordinate point (r3, π). Let the polar coordinate of the corresponding point on the circumference of the plane unit circle ζ be (1, π). Then, R3 can be obtained from equation (14) and α3 = β3 = π, where(17c)R3=r31+∑k=0∞ckcosk+1π.The three formulas above give the relationship betweenR and ck. From equation (14), when k = 0, 1, 2, …, n, n takes a finite positive integer, and equation (14) cannot be accurately satisfied, but ck (k = 0, 1, 2, …, n) can be determined by optimization techniques. Meanwhile, the objective function can be taken as(18)f=∑i=1m+1ri−Ricosβi−αi+∑k=0nckcoskβi+αi2.In equation (18), ri and αi are known values, ck is an unknown value, and βi is also an unknown value, but it can be connected with ck and αi through equation (13).Whenm > n, ck∗ that minimizes f can be obtained by equation (18). ck∗ is the mapping function coefficient ck which is to be determined, and now fmin > 0.According to the actual parameters of the trapezoidal roadways in one coal mine, leta = 2.3, b = 1.5, c = 1.5, n = 3, and m = 4. In order to facilitate the calculation, take 5 special points of the trapezoidal roadway, and their polar coordinates are (c, 0), (b2+c2, arctanb/c), (a+b/2, π/2), (a2+c2, π−arctana/c), (c, π). Substitute them into equations (8) and (9); ck is calculated by computer as c0 = −0.1011, c1 = −0.1101, c2 = 0.0859, c3 = 0.0014, and R = 1.69, and the mapping function of the trapezoidal section is determined as(19)z=ωζ=1.69×ζ−0.0896−0.1101ζ+0.0859ζ2+0.0014ζ3.


## 3.2.2. Analytical Formula for Tangential Stress of Surrounding Rock around Trapezoidal Roadway

Now perform calculation according to the steps of the complex variable function solution method; the following equation is derived by equation (6):(20)B=14P+λP,B′=12λP−P,C′=0.From equation (11), the following equation is derived:(21)ωσ=Rσ+c0+c11σ+c21σ2+c31σ3,ωσ¯=R1σ+c0+c1σ+c2σ2+c3σ3,ω′σ=R1−c11σ2−c22σ3−c33σ4,ω′σ¯=R1−c1σ2−2c2σ3−3c3σ4.Introducing equation (21) into equation (9) gives(22)f0=−2BRσ+c0+c11σ+c21σ2+c31σ3−B′R1σ+c0+c1σ+c2σ2+c3σ3.Introducing equation (22) into equation (4) to solve with Harnack’s theorem, the integral on the right is(23)12πi∫σf0σ−ζdσ=−2BRζ−B′Rc1ζ−B′Rc2ζ2−B′Rc3ζ3.Calculate by combining equations (4) and (13); then the solutions from α1 to αi are found:(24)α1=−BR,α2=−0.2746B′R,α3=−0.0096B′R,α4=α5=⋅⋅⋅=αn=0.Introduce them into formula (7) to get(25)φ0ζ=−BRζ−0.2746B′Rζ2−0.0096B′Rζ3.Substituting equations (10) and (25) into equation (3) gives(26)φζ=BRζ+c0+c1ζ−1+c2ζ−2+⋯+cnζ−n−BRζ−0.2746B′Rζ2−0.0096B′Rζ3.From equation (2),(27)Φζ=φ′ζω′ζ=B−c1ζ2−2c2ζ−3c3ζ4−c1ζ2−2c2ζ−3c3−B′0.5492ζ5+0.0288ζ7ζ4−c1ζ2−2c2ζ−3c3.Letρ of the unit circle which has already been mapped be 1; then ζ=ρeiθ=eiθ=cosθ+isinθ. Substitute ζ into equation (27). After the calculation of the first half of equation (27), we have the following result:(28)B−c1ζ2−2c2ζ−3c3ζ4−c1ζ2−2c2ζ−3c3=−Bc1cos2θ+2c2cosθ+3c3+isin2θ+2c2sinθcos4θ−c1cos2θ−2c2cosθ−3c3+isin4θ−c1sin2θ−2c2sinθ,(29)Z1=x1,y1=c1cos2θ+2c2cosθ+3c3+isin2θ+2c2sin2θ,Z2=x2,y2=cos4θ−c1cos2θ−2c2cosθ−3c3+isin4θ−c1sin2θ−2c2sinθ.Then,(30)ReZ1Z2=x1x2+y1y2x22+y22.From equation (30),(31)x1=c1cos2θ+2c2cosθ+3c3,x2=cos4θ−c1cos2θ−2c2cosθ−3c3,y1=sin2θ+2c2sinθ,y2=sin4θ−c1sin2θ−2c2sinθ.The calculation of the last half part of equation (27) gives(32)−B′0.5492ζ5+0.0288ζ7ζ4−c1ζ2−2c2ζ−3c3=−B′0.5492cos5θ+0.0288cos7θ+i0.5492sin5θ+0.0288sin7θcos4θ−c1cos2θ−2c2cosθ−3c3+isin4θ−c1sin2θ−2c2sinθ,(33)Z3=x3,y3=0.5492cos5θ+0.0288cos7θ+i0.5492sin5θ+0.0288sin7θ.Then,(34)ReZ3Z2=x2x3+y2y3x22+y22.From equation (34),(35)x3=0.5492cos5θ+0.0288cos7θ,y3=0.5492sin5θ+0.0288sin7θ.Substitute equations (28), (30), (32), and (34) into formula (1). Because the value of σρ around the roadway is 0, we have σθ=4ReΦζ. The analytical equation for the tangential stress around the trapezoidal roadways is calculated as(36)σθ=P+λPx1x2+y1y2x22+y22+2P−λPx2x3+y2y3x22+y22.


## 4. Engineering Application Examples

The actual geological conditions and roadway parameters of the coal mine are as follows: the long side of the trapezoidal roadway is 4.6 m, the short side is 3 m, the height is 3 m,P is 10 MPa, and the lateral pressure coefficient is 0.8.Substitute the parameters above into equations (36), (34), and (32) to obtain the tangential stress values of surrounding rocks around the trapezoidal roadway, as in Table 1. According to Table 1, draw the tangential stress fitting curve of surrounding rocks around the trapezoidal roadways as shown in Figure 4; and the tangential stress distribution diagram of surrounding rock around the trapezoidal roadways is shown in Figure 5.Table 1
Tangential stress value at each point around trapezoidal roadways (unit: MPa).
0°15°30°45°60°75°90°105°120°135°150°165°180°Tangential stress13.79212.67810.7079.4969.53911.13214.78519.619.70614.4510.1397.9947.372Figure 4
Fitting curve of tangential stress value of surrounding rock around trapezoidal roadways.Figure 5
Schematic diagram of tangential stress distribution in surrounding rocks of trapezoidal roadways.As illustrated in Figures4 and 5, the tangential stress distribution of surrounding rocks around trapezoidal roadways is symmetrical, and the tangential stress rangeability from the roof (0° angle) to the top angle (45° angle, 315° angle) is relatively small, while the two sidewalls (60° angle, 300° angle) transit to the bottom angle (120° angle, 240° angle) and then to the floor (180° angle), the tangential stress varies greatly, and the maximum difference of stress is 12.3 MPa.After specific analysis, as shown in Figure4, the tangential stress concentration degree of surrounding rock around trapezoidal roadway bottom angle (110° angle) > two sidewalls (90° angle) > roof (0° angle) > roof (45° angle) > floor (180° angle), where the maximum stress concentration factor of the bottom angle is 2, the second is 1.48 in the middle of the two sidewalls, and the minimum of the floor is 0.74. According to the observation from the working face of a certain mine, the roadway roof deformation is not obvious. The main deformation positions are on the two sidewalls and the floor. The maximum deformation reaches 1500 mm and 1000 mm, respectively, and the deformation patterns of the two sidewalls are similar to the tangential stress distribution state in Figure 5. The deformation of the upper part on the sidewall is small and the deformation of the lower sidewall is large, which has seriously affected the normal mining demand of the working face. Analyzing the reasons, it is believed that the tangential stress variation range of the surrounding rock around the roof is relatively gentle, so that the overall roof has a strong bearing capacity, and the stress concentration coefficient is not large, which is the reason for the fact that the deformation of the roof is not obvious; while the tangential stress of the surrounding rock from top to bottom of the two sidewalls has a big rangeability, and the surrounding rock’s self-stabilization ability is feeble. Under the influence of high-strength mining, large deformation occurs at the lower sidewall of the coal and rock mass where the stress concentration is greater. Similarly, the tangential stress from the bottom corner to the middle of the floor has a big rangeability, and the stress concentration coefficient in the middle of the floor is less than 1. The self-stable load-bearing capacity of the floor is weaker, and the large deformation of the two sidewalls reduces the load-bearing capacity, resulting in the transfer of high ground stress to the floor, which causes the floor to produce an extrusion type floor heave and tensile fracture.Based on the above complex variable theory calculation and field engineering survey, it can be seen that, for surrounding rock of trapezoidal roadways, the large deformation of surrounding rock is mainly on the two sidewalls and floor, where the tangential stress has a big rangeability. Therefore, the two sidewalls and floor are the key parts to support the trapezoidal roadway.


## 5. Surrounding Rock Control Technology of Trapezoidal Roadway

### 5.1. Support Plan and Key Parameters

Based on the above-mentioned trapezoidal roadway surrounding rock tangential stress distribution characteristics and field engineering roadway surrounding rock deformation characteristics, propose a support plan: bolting with wire mesh and anchor cable + W-type steel belt + joist steel shed support + anchor cable grouting. The construction steps are as follows.First of all, after the excavation of the roadway, the bolting with wire mesh and anchor cable + W-type steel belt should be timely supported. Then observe the change rate of stress and displacement of the roadway surrounding rock, and carry out the secondary support of the joist steel shed support with high strength and large stiffness when the deformation rate of roadway is slow. According to the on-site monitoring results of the coal mine, the support position of the joist steel shed generally lags behind the tunneling face by 20～40 m. For the geological conditions of the trapezoidal roadway, the coal seam and the roof and the floor rock strength are low, and the two sidewalls of the roadway will produce a large amount of deformation, which cannot be stable by using the joist steel bracket only. At the same time, the stress of the surrounding rock at the roof is imposed on the top beam, and the top beam is transmitted to the floor through the shed’s legs. If the floor is not supported, the floor heave phenomenon will inevitably occur. Therefore, three times of anchor cable grouting support should be carried out at the end. The support time is selected after the secondary support when the two sidewalls of the roadway are further deformed and compressed backplane, that is, shortly after being affected by the mining stress but before the influence of the advanced support stress of the working face. The secondary cracks in the roadway are relatively developed at this time, the grouting effect is good, and the damage and destruction of the joist steel bracket and the floor heave caused by the large deformation of the two sidewalls under the influence of mining are avoided.The initial and timely support layout is shown in Figure6.Figure 6
Initial and timely support layout. (a) Roadway section support diagram. (b) Roof support plan. (c) Sidewalls support plan.
(a)(b)(c)Combined with the aforementioned theoretical analysis results of the tangential stress distribution of the trapezoidal roadways, it is believed that the stress distribution of the roof is stable, and the stress concentration is only slightly higher in the middle of the roof. Therefore, one anchor cable is installed in each row in the middle of the roof. The tangential stress distributions of the two sidewalls are extremely uneven, and the stress concentration of the lower sidewall is large, so the lower sidewall is the most deformed part. When the timely support is operating, do not install the anchor cable, allow the sidewalls to deform properly so that it will release part of the deformation energy first, which could increase the development degree of secondary cracks, and create conditions for the third-time anchor cable grouting support. Only one anchor cable is installed in each row of the upper sidewall to increase the strength of the upper sidewall support. Rock bolts should be laid first, and anchor cables should be laid uniformly when carrying out joist steel shed support. Support parameters are mainly selected for the existing support equipment in the mine.The secondary support arrangement is shown in Figure7.Figure 7
Joist steel shed support diagram.Six rows of joist steel are laid on the roof as well as the two sidewalls at a distance of 500 mm, and then the backplate is laid for anchor cable support. The anchor cable adopts Φ22 mm × 4000 mm high-strength prestressed and low-relaxation steel strands. The installation position is shown in Figure6; the roof and upper anchor cables use 300 mm × 300 mm × 14 mm pallets and 800 mm I-beam to provide a reaction force for the support, which is equivalent to adding some nodes at the support and effectively enhancing the bearing capacity of the metal support.The third-time anchor cable grouting support arrangement is shown in Figure8.Figure 8
Third-time anchor cable grouting support arrangement plan.The deep hole grouting adopts a 5 m grouting anchor cable. First of all, the grouting anchor cable is driven into the surrounding rock, and then deep grouting is carried out. Generally, the grouting pressure is controlled at 3～4 MPa, and the grouting diffusion radius is 1.5～2.0 m. Drilling at an angle of 15° can make the anchor cable hit the rock layer of the floor, and the spreading range of the grout can spread to the lower side and floor where the stress concentration is large. By improving the mechanical parameters of the surrounding rock in the grouting reinforcement area and enhancing the integrity of the loose coal body, the surrounding rock of the two sidewalls and the floor is transformed from the load-applying body of the joist steel shed support into the joint bearing body.In conclusion, the joint support plan of the surrounding rock of trapezoidal roadway is obtained finally as shown in Figure9.Figure 9
Combined support scheme for trapezoidal roadway. (a) Sectional drawing of joint support scheme. (b) Two-sidewall plan of joint support plan.
(a)(b)


### 5.2. Effect of Application

A rail roadway in this mine is chosen for the experiment. The rail roadway was influenced twice by the mining activities: one occurred in 105(7) working face and the other occurred in 105(8) working face. The above supporting scheme was applied in the actual construction of test roadway. Using cross distribution method to observe the deformation of roadway surrounding rock and the observation station layout is shown in Figure10.Figure 10
Observation station layout.The monitoring data of No. I observation station are shown in Figures11 and 12. The monitoring data indicate that after mining influence of 105(7) working face and the goaf is gradually stable, the maximum displacement of two sidewalls of roadway reaches 210 mm, the maximum displacement of the roof and floor reaches 140 mm, and the floor heave is up to 120 mm. During the advancing process of 105(8) working face, two sidewalls of roadway amplification are 50 mm, and floor heave amplification is 20 mm, which all mean the deformation of roadway is small. The support effect is shown in Figure 13. The surface of the roadway is flat, and the pillars of the roadway’s advanced support section do not bend or break under the condition of wood support. The support scheme is proved to meet the requirements of high-intensity mining activities and will ensure the safe and efficient mining in the whole working face.Figure 11
Observation of roadway displacement influenced by mining at upper section working face(105 (7) working face).Figure 12
Observation of roadway displacement influenced by mining in this section working face(105 (8) working face).Figure 13
Optimization supporting effect of trapezoidal roadway. (a) Middle section. (b) Advanced support section.
(a)(b)


## 5.1. Support Plan and Key Parameters

Based on the above-mentioned trapezoidal roadway surrounding rock tangential stress distribution characteristics and field engineering roadway surrounding rock deformation characteristics, propose a support plan: bolting with wire mesh and anchor cable + W-type steel belt + joist steel shed support + anchor cable grouting. The construction steps are as follows.First of all, after the excavation of the roadway, the bolting with wire mesh and anchor cable + W-type steel belt should be timely supported. Then observe the change rate of stress and displacement of the roadway surrounding rock, and carry out the secondary support of the joist steel shed support with high strength and large stiffness when the deformation rate of roadway is slow. According to the on-site monitoring results of the coal mine, the support position of the joist steel shed generally lags behind the tunneling face by 20～40 m. For the geological conditions of the trapezoidal roadway, the coal seam and the roof and the floor rock strength are low, and the two sidewalls of the roadway will produce a large amount of deformation, which cannot be stable by using the joist steel bracket only. At the same time, the stress of the surrounding rock at the roof is imposed on the top beam, and the top beam is transmitted to the floor through the shed’s legs. If the floor is not supported, the floor heave phenomenon will inevitably occur. Therefore, three times of anchor cable grouting support should be carried out at the end. The support time is selected after the secondary support when the two sidewalls of the roadway are further deformed and compressed backplane, that is, shortly after being affected by the mining stress but before the influence of the advanced support stress of the working face. The secondary cracks in the roadway are relatively developed at this time, the grouting effect is good, and the damage and destruction of the joist steel bracket and the floor heave caused by the large deformation of the two sidewalls under the influence of mining are avoided.The initial and timely support layout is shown in Figure6.Figure 6
Initial and timely support layout. (a) Roadway section support diagram. (b) Roof support plan. (c) Sidewalls support plan.
(a)(b)(c)Combined with the aforementioned theoretical analysis results of the tangential stress distribution of the trapezoidal roadways, it is believed that the stress distribution of the roof is stable, and the stress concentration is only slightly higher in the middle of the roof. Therefore, one anchor cable is installed in each row in the middle of the roof. The tangential stress distributions of the two sidewalls are extremely uneven, and the stress concentration of the lower sidewall is large, so the lower sidewall is the most deformed part. When the timely support is operating, do not install the anchor cable, allow the sidewalls to deform properly so that it will release part of the deformation energy first, which could increase the development degree of secondary cracks, and create conditions for the third-time anchor cable grouting support. Only one anchor cable is installed in each row of the upper sidewall to increase the strength of the upper sidewall support. Rock bolts should be laid first, and anchor cables should be laid uniformly when carrying out joist steel shed support. Support parameters are mainly selected for the existing support equipment in the mine.The secondary support arrangement is shown in Figure7.Figure 7
Joist steel shed support diagram.Six rows of joist steel are laid on the roof as well as the two sidewalls at a distance of 500 mm, and then the backplate is laid for anchor cable support. The anchor cable adopts Φ22 mm × 4000 mm high-strength prestressed and low-relaxation steel strands. The installation position is shown in Figure6; the roof and upper anchor cables use 300 mm × 300 mm × 14 mm pallets and 800 mm I-beam to provide a reaction force for the support, which is equivalent to adding some nodes at the support and effectively enhancing the bearing capacity of the metal support.The third-time anchor cable grouting support arrangement is shown in Figure8.Figure 8
Third-time anchor cable grouting support arrangement plan.The deep hole grouting adopts a 5 m grouting anchor cable. First of all, the grouting anchor cable is driven into the surrounding rock, and then deep grouting is carried out. Generally, the grouting pressure is controlled at 3～4 MPa, and the grouting diffusion radius is 1.5～2.0 m. Drilling at an angle of 15° can make the anchor cable hit the rock layer of the floor, and the spreading range of the grout can spread to the lower side and floor where the stress concentration is large. By improving the mechanical parameters of the surrounding rock in the grouting reinforcement area and enhancing the integrity of the loose coal body, the surrounding rock of the two sidewalls and the floor is transformed from the load-applying body of the joist steel shed support into the joint bearing body.In conclusion, the joint support plan of the surrounding rock of trapezoidal roadway is obtained finally as shown in Figure9.Figure 9
Combined support scheme for trapezoidal roadway. (a) Sectional drawing of joint support scheme. (b) Two-sidewall plan of joint support plan.
(a)(b)


## 5.2. Effect of Application

A rail roadway in this mine is chosen for the experiment. The rail roadway was influenced twice by the mining activities: one occurred in 105(7) working face and the other occurred in 105(8) working face. The above supporting scheme was applied in the actual construction of test roadway. Using cross distribution method to observe the deformation of roadway surrounding rock and the observation station layout is shown in Figure10.Figure 10
Observation station layout.The monitoring data of No. I observation station are shown in Figures11 and 12. The monitoring data indicate that after mining influence of 105(7) working face and the goaf is gradually stable, the maximum displacement of two sidewalls of roadway reaches 210 mm, the maximum displacement of the roof and floor reaches 140 mm, and the floor heave is up to 120 mm. During the advancing process of 105(8) working face, two sidewalls of roadway amplification are 50 mm, and floor heave amplification is 20 mm, which all mean the deformation of roadway is small. The support effect is shown in Figure 13. The surface of the roadway is flat, and the pillars of the roadway’s advanced support section do not bend or break under the condition of wood support. The support scheme is proved to meet the requirements of high-intensity mining activities and will ensure the safe and efficient mining in the whole working face.Figure 11
Observation of roadway displacement influenced by mining at upper section working face(105 (7) working face).Figure 12
Observation of roadway displacement influenced by mining in this section working face(105 (8) working face).Figure 13
Optimization supporting effect of trapezoidal roadway. (a) Middle section. (b) Advanced support section.
(a)(b)


## 6. Conclusion

This paper uses complex variable theory and elastic theory to study the tangential stress of surrounding rock of trapezoidal roadways and draws the following conclusions:(1)
The analytical solution is obtained by the complex change theory, the stress value obtained by calculation is only related to the ground stress, the lateral pressure coefficient, and the size of the roadway section, and it has nothing to do with the mechanical parameters of the surrounding rock of the roadway. Therefore, this value is not the actual stress value of the engineering roadway’s surrounding rock. The analytical solution is consistent only in the regularities of stress distribution, which can reflect the stress state of a roadway with a certain shape section and can provide a theoretical basis for the selection of the shape and size of the roadway section and the support scheme.(2)
The tangential stress fitting curve of surrounding rock around the roadway shows that the top corner of the trapezoidal roadways is not affected by stress concentration, and the maximum stress concentration degree at the bottom corner is 2. The overall stress distribution of the roof is uniform, while the stress distribution of the two sidewalls and the floor varies greatly, and the maximum stress difference is 12.3 MPa. In addition, the stress concentration of the lower part of the sidewall and bottom corner is large, so reinforcement measures should be taken in the design of roadway support.(3)
The results of the study on the surrounding rock control of the trapezoidal roadways show that the design of the trapezoidal roadway support should take the two sidewalls and floor near the bottom corner as the main support focus. Based on the design of bolting with wire mesh and anchor cable + W-type steel belt + joist steel shed support, this paper proposes three-time grouting reinforcement for the middle of the two sidewalls, the lower sidewall, and the floor. The field application results show that the support scheme can meet the high-intensive mining requirements of the working face, and there is no need for repair during the service period.



---

*Source: 6381785-2021-08-04.xml*