## DETERMINATION OF STRESS AT THE CONTOUR POINTS OF SINGLE UNDERGROUND HORIZONTAL MINE WORKINGS, SUBJECT TO UNIFORM PRESSURE, AND CALCULATION OF THEIR PERMISSIBLE LAYING DEPTH

**Authors:**Bogomolov A.N^{1}, Ushakov A.N^{1}, Bogomolova O.A^{1}**Affiliations:**- Volgograd State Technical University

**Issue:**Vol 12, No 3 (2021)**Pages:**117-129**Section:**ARTICLES**URL:**https://ered.pstu.ru/index.php/CG/article/view/1917**DOI:**https://doi.org/10.15593/2224-9826/2021.3.12- Cite item

# Abstract

The results of constructing the contour of an underground horizontal mine working of given dimensions and geometric structure based on the method proposed by the authors are presented, the stress state at the points of the contours of the working, the cross-sections of which have the shape of a vault with inclined walls and an ellipse, are given, the permissible laying depths are determined at given values of uniform pressure applied at the points of their contours. The criterion for determining the permissible laying depth is the condition that there are no points on the contours of the workings at which the normal tangential stresses exceed the tensile and compressive strength of the enclosing rock.

# Full Text

It is well known that one of the central tasks of geomechanics is the problem of assessing the stress state of rocks around single workings [1-4] in an elastic isotropic massif, moreover, the problem associated with the study of the stress state at the points of the mine working contours is of considerable interest, since the solution of this problem is associated with the problem of their strength. The geometric structure of the contours of mine workings when solving two-dimensional problems of the linear theory of elasticity [5-8] can be obtained using complex variable functions that carry out a conformal mapping of one of the canonical domains to the domains of interest to us. The mapping functions are polynomials representing the simplest and most well-studied class of functions [9-11]. Let’s use the mapping function presented in [12; 13], which looks like this (1) where - are real numbers; this function maps the interior of the unit circle to the exterior of an infinite simply connected region, the boundary of which is a family of simple closed curves. The mapping function (1) was previously used by the authors of this article in [14-16] when solving problems related to the study of stresses on the contour of mine workings and to determine the permissible laying depth of horizontal workings with a given cross-sectional configuration. Consider an underground horizontal working, in which shape of its section is determined using the mapping function (1). Let’s assume that this mine working is located at a sufficient depth H, and an internal uniform pressure of intensity p acts along the mine working contour, which makes it possible to consider the working as an underground storage, for example, of hydrocarbons of a given size and configuration. The purpose of this work is to study the stress state at the points of the contour of an underground horizontal working, built on the basis of the algorithm proposed by the authors, and to determine the permissible laying depth of mine working at various values of the internal uniform tensile pressure acting on this contour. Let’s assume that there are no points on the mine working contour at which the normal tangential stresses exceed the tensile and compressive strengths that is [15]. As an example, consider two types of mine workings of a given size: a vault with inclined walls and also a working with an elliptical cross-section. Let's proceed to the construction of the contour of the working, the cross-section of which is a vault with inclined walls of 4 m height and 5 m width. Now let’s examine a vaulted hole and associate an coordinate system with it. Let be its height, - distance from the lower base to the abscissa axis, parameters and be the values of half the lengths of the lower and upper bases of the contour (Fig. 1). In (1) let Then, separating real component from the imaginary, we obtain the equations for the hole contours in the parametric form: (2) where To determine the values of the coefficientsand parameter : 1) at i.e. ; (3) 2) at , i.e. . (4) Subtracting and adding expressions (3) and (4), we get (5) (6) 3) at , i.e. (7) where с - is a real positive number. Hence, taking into account (6) and (7), we have (8) Let the points with arguments , be located on the lower base and on the side of the contour respectively. Then: at Taking into account (2) and (8), we get (9) where (10) at Considering (2), we have (11) where (12) Solving the system of equations (9) and (11), we find (13) (14) Substituting the obtained values of and into the formulas (5), (7), (8), we obtain the values of all coefficients of the mapping function (1) and parameter The following value is assumed as a width of mine working . (15) Considering the parametric equations of curves (2) describing the contours of underground horizontal workings, it should be noted that the following equalities hold: . Thus, for any values of coefficients the contours will be symmetric about the ordinate axis, and so let’s assume In this case, the width of the working will be equal to twice the value of the function at the point of it’s minimum. Let be an argument of the function (16) at which this function has a minimum value. Then argument can be calculated by the equation (17) By construction, function (16) takes negative values for any values of the argument , moreover, since , then, according to Rolle theorem [17], we can conclude that function (16) at interval has a single extremum point, and the equation (17) has a unique real solution. Note that when calculating numerical values of and , we first assume . Substituting the obtained coefficients of the mapping function (1) or parametric equations of the mine working contour (2) into the equation (17), we obtain the real root of that equation, providing the minimum value of function (16), at that ,. Further, assuming that and substituting the resulting value c into expressions for and , let’s recalculate the values of coefficients and parameter Now, using the obtained values of coefficients, let’s construct a vault with inclined walls with specific values of its dimensions: high and wide. Let . Then the angle of inclination of the side of mine working contour to its base will be Let and . Then, using formulas (13) and (14), considering (10) and (12), we get (18) at Substituting coefficients (18) into the equation (17) results in and gives a unique solution . Then, according to (2), we have Thus, and, therefore, Substituting the obtained value into the ratios for and , according the formulas (13), (14), (5) and (6) we get (19) at Equation (17), taking into account the values of the coefficients (19), takes the form and has a solution . Substituting this value of argument into (2), we get Now and, therefore, , which clarifies the result. It is easy to see that it is possible to construct a set of homothetic contours of mine workings, for which it is necessary to multiply the coefficients of the mapping function by the homothety coefficient k<0. The new contours of the workings will have the same configuration, but with linear dimensions changed k times. Note that mapping function (1) with coefficients (19) performs a conformal mapping of the interior of the unit circle on the obtained area. To verify this fact, it is necessary to establish that the derivative of function (1), regular in the unit circle everywhere, except for the origin, at which it has a simple pole, and for does not vanish at any point on the unit circle . That means, that the equation should not have solutions inside the unit circle . Applying Rouche theorem [18], we get the condition , which holds when the values of the coefficients (19) of the mapping function (1) are substituted into it. Thus, the required underground horizontal working with a specified cross-sectional contour is constructed. The mine working contour described by parametric equations (2) with coefficients (19) is shown on Fig. 2. Following [12], notice that formula, describing the stress state at the points of the mine working contour, the configuration of which is determined using the mapping function (1) under the condition of all-round uniform pressure of a given intensity p applied to the contour, has the form (20) where (21) ; ; (22) ; , , , (23) where volumetric weight of rock; - lateral thrust coefficient; H - laying depth of mine working,the value of uniform pressure applied to the contour of the hole; moreover, according to [2], we assume that at the contour of mine working experiences compression of constant value p, and at expansion of the same intensity. Formula (20) was obtained under the condition that the laying depth H of the working is large enough. Following [4], let’s presume , (24) where the largest linear dimension of the mine working section. Finding zero values of normal tangential stresses is reduced to solving the equation (25) where , . Taking into account the results of [16], note that the extreme values of the function can be obtained from the following equations (26) where (27) Now, using the above relations, let’s consider the following problems. Problem А. Conduct a study of the stress state on the contours of underground horizontal workings with a cross-sections in the shape of a vault with inclined walls and a shape of ellipse under the pressure of uniform all-round expansion of a given intensity at a certain value of the lateral thrust coefficient of the enclosing massif. Problem B. Determine the permissible values of the laying depth of the underground horizontal workings with a cross-sections in the shape of a vault with inclined walls and the shape of an ellipse under the action of uniform all-round expansion of a given intensity at a given coefficient of lateral expansion of the enclosing massif. Consider as a granite host rock with volumetric weight t/m3 tinsel strength 1735 т/m2 and compression strength 20400 т/m2. Since the largest linear size of the considered mine working is its width, equal to 5m, then, according to (24), we set Then 625 т/m2. When solving these problems, let’s set the value of the lateral thrust coefficient to , which corresponds to the value of Poisson's ratio, typically accepted equal to [4] for rock. Solution to problem A 1. Vault with inclined walls Substituting the values of the coefficients (19) into (1), we obtain an underground horizontal working with a cross-section shown in Fig. 2 and described by parametric equations (2). Coefficients (22) are independent of the lateral thrust coefficient and are equal to (28) Substituting the values of the coefficients (19) into formulas (21), we obtain . (29) As the values of the stresses applied to the mine working contour, we enter: 0 t/m2,102 t/m2, 408 t/m2. (30) Then direct calculations by formula (25) taking into account (28), (29) give us the following values of and for which : for 0 т/ ; for 102 т/ ; for 408 т/ . The performed calculations make it possible to identify areas on which tensile stresses act at given values of uniform pressure. For 0 t/m2 we have ; For 102 t/m2 we have ; For 408 t/m2 we have . On other sections of interval compression stresses act. Diagrams of tangential normal stresses for mine working of considered shape of cross-section at given values of uniform pressure and are shown on Fig. 3. а b c Fig. 3. Diagrams of tangential normal stresses on the contour at: 0 t/m2 (а), 102 t/m2 (b), 408 t/m2 (c) 2. Mine working with elliptical section Assuming in the mapping function (1) we get (31) Function (31) carries out a conformal mapping of the interior of the unit circle on an infinite plane with an elliptical hole, and the circles corresponds to an ellipse centered at the origin and semiaxes , . By setting values of and , ellipses of any shape and size can be obtained. If , than ellipse becomes a circle. In the extreme case of the ellipse turns into axis segment of length between the dots , and the region turns into an infinite plane with a rectilinear slit. According to (22) (32) Taking into account relations (32), formulas (21) give us: ,, . (33) Substituting the values of the coefficients (33) into equation (25), we obtain the equation for finding zero values in the form (34) The formula for normal tangential stresses, taking into account the values of the coefficients (32), becomes (35) Assuming in (35), we obtain the formula for normal tangential stresses for generating a circular cross-section in the form (36) Let , , . The elliptical holes attained at these values are shown in Fig. 4. а b c Fig. 4. Elliptical shapes of the section at: (а), (b), (c) Assuming, like in previous case, 625 t/m2, and, sequentially choosing the values and , using equation (34), let’s find the areas on which tensile stresses act at given values of uniform pressure (30). We have for : at 0 t/m2 ; at 102 t/m2 ; at 408 t/m2 for : at 0 t/m2 ; at 102 t/m2 ; at 408 t/m2 for : at 0 t/m2 ; at 102 t/m2 ; at 408 t/m2 ; compressive stresses act on other segments of the interval . As an example, below are diagrams of normal tangential stresses for wine working with circular cross-section at given values of uniform pressure (Fig. 5). а b c Fig. 5. Diagrams of normal tangential stresses for mine working with circular cross-section at: 0 t/m2 (а), 102 t/m2 (b), 408 t/m2 (c) Solution to problem B 1. Vault with inclined walls Let . Then, substituting (28) and (29) into (27) and (26), we obtain an equation for the extreme values of the function , the root of which is the value . By attaching values and to it, we get , . Based on the results of [14; 15], we obtain , . Therefore, at the depth the considered horizontal working in a rock massif, folded by granite, loses strength. Linear dimensions of the working, according to (24), should not exceed , Let a uniform confining the sile pressure of the102 t/m2 magnitude act at the points of the mine working contour. Consider the expression (37) where tensile strength of the rock as a function of the mine working laying depth . By differentiating expression (37), and by obtaining for determining the extreme values of the function together with equation trigonometric equation of the ninth degree, let’s look for the minimum value of the function near the value of the argument at which the depth function reaches a minimum at . This value at 102 t/m2 is . And so, the permissible laying depth of the vault with inclined walls at 102 t/m2 equals . Now let the internal uniform tensile pressure at the points of the working contour be 408 t/m2. The result of caring out the necessary calculations, is that the permissible laying depth of the considered horizontal working is equal to . 2. Mine working with elliptical section Let . Taking into account relations (27), (32), and (33) on the interval equation (26) results in value . With and we obtain values of arguments for which function has extreme values. Further, let . Then , . Taking into account the results of [14; 15], we obtain , . Therefore, at the laying dept the considered mine working in the rock massif, folded by granite, loses its strength. According to (24), the linear dimensions of the working should not exceed , , where . Carrying out further calculations similar to those done for the previous case, we obtain: for : at 102 т/ ; at 408 т/ ; for : at 0 т/ ; at 102 т/ ; at 408 т/ ; for : at 0 т/ ; at 102 т/ ; at 408 т/ . All problems have been solved.# About the authors

### A. N Bogomolov

Volgograd State Technical University

### A. N Ushakov

Volgograd State Technical University

### O. A Bogomolova

Volgograd State Technical University

# References

- Turchaninov I.A., Iofis M.A., Kasparyan E.V. Basics of Rock Mechanics. Leningrad, Nedra, 1989, 488 p.
- Bulychev N.S. Mechanics of underground structures. Moscow, Nedra, 1989, 270 p.
- Glushko V.T., Gavelya S.P. Assesment of the tensly-deformed rock massifs. Moscow, Nedra, 1986, 221 p.
- Tsimbarevich P.M., Rock mechanics. Moscow, Ugletechizdat, 1948, 184 p.
- Muskhelishvili N.I. Some basic tasks of mathematical theory of elasticity. Moscow, Science, 1966, 708 p.
- Kolandia A.I. Mathematical methods of two-dimensional elasticity. Moscow, Science, 1973, 304 p.
- Khan H. Theory of elasticity. The basics of the linear theory of its application. Moscow, Peace, 1988, 344 p.
- Paron V.Z., Perlin P.I. Methods of mathematical theory of elasticity. Moscow, Science, 1981, 688 p.
- Filchakov P.F. Approximate conformal mapping methods. Kiev, Naukova Dumka, 1964, 536 p.
- Filchakova V.P. Сonformal mappings of regions of a special type. Kiev, Naukova Dunka, 1972, 252 p.
- Savin G.N. Distribution of stresses near the holes. Kiev, Naukova Dumka, 1968, 888 p.
- Tsvetkov V.K. Determination of the shapes of sections of mine workings with specified stresses on the contours. Physical and technical problems of mining, 1986, no. 2, pp. 24-29.
- Tsvetkov V.K. Calculation of rational parameters of mine workings. Moscow, Nedra, 1993, 251 p.
- Bogomolov A.N., Ushakov A.N. Methods of the theory of functions of a complex variable in problems of geomechanics. Volgograd, VSACU, VSSPU “Change,” 2014, 227 p.
- Bogomolov A.N., Bogomolova O.A., Ushakov A.N. Determination of the laying depth of the horizontal working based on the analysis of the stress state of the enclosing massif. Bases, foundation and soil mechanics, 2017, no. 5, pp. 2 - 9.
- Bogomolov A.N., Bogomolova O.A., Ushakov A.N. On stresses at the contour points of single underground workings of various cross-sections, subject to all-round uniform pressure. Vestnik PNIPU. Stroitel’stvo and arkhitektura, 2018, vol. 9, iss. 3, pp. 54-70. doi: 10.15593/2224-9826/2018.3.06.
- Nikolski S.M. Mathematical analysis course. Мoscow, Science, 1990, 528 p.
- Sidorov U.V., Fedoruk M.V., Schabunin M.I. Lectures on the theory of functions of a complex variable. Мoscow, Science, 1989, 480 p.

# Statistics

#### Views

**Abstract** - 816

**PDF** (Russian) - 81

### Refbacks

- There are currently no refbacks.