AP01_45.vp 1 Introduction In July 1995 our institute was awarded funding for a pro- ject on “Model experiments on the pressure distribution in a coal seam, or in a wide coal pillar, before and after bump initiation” in the framework of the Volkswagen Foundation. The project was led by the Technische Universität München, represented by Prof. Dr. Dr.h.c. Horst Lippmann, and the coapplicant was the of the Czech Technical University in Prague, Klokner Institute represented by Jaroslav Vacek. TUM supervised the project and supplied technicians to operate the 10MN testing machine and to help with assem- bling and disassembling the device. This took a total of 30 days, in which the testing machine was exclusively reserved for the Volkswagen Project. All other research work was carried out in Prague. Re- search work is now continuing in the framework of a Czech grant. 2 Testing devices 2.1 Loading cell Fig. 1 shows the loading cell. It consists of the lower steel tank, which is designed for the horizontal forces caused by the vertical load in araldite specimens. The loading cell is equipped with plexiglass on its sides which allows the samples to be observed during the tests. The tank is shown on Photo 2. This loading cell models (simulates) the rock mass in the vicinity of the seam. In the loading cell we placed two araldite specimens (with dimensions of 160/400/40 mm), which model a real seam. The gap between them corresponds to the width of a working gallery in a mine. We observed the mechanism and the history of coal bumps. The araldite specimen was covered with a soft dural sheet, and a force meters were placed on it in the following manner: 5 comparatively thick force meters were placed near its outer edge and another 15 thin- ner force meters were placed next to them, (see Fig. 1 and 54 © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ Acta Polytechnica Vol. 41 No. 4 – 5/2001 Stress Distribution in a Coal Seam before and after Bump Initiation J. Vacek This paper deals with to the behaviour of open rock that occurs, for example, during longwall mining in coal mines, in deep tunnel, or shaft excavation. Longwall instability leads to extrusion of rock mass into an open space. This effect is mostly referred to as a bump, or a rock burst. For bumps to occur, the rock has to possess certain particular rock burst properties leading to accumulation of energy and the potential to release this energy. Such materials may be brittle, or the bumps may arise at the interfacial zones of two parts of the rock, that have principally different material properties. The solution is based on experimental and mathematical modelling. These two methods have to allow the problem to be studied on the basis of three presumptions: • the solution must be time dependent • the solution must allow the creation of crack in the rock mass • the solution must allow an extrusion of rock into an open space (bump effect) Keywords: rock bursts, bumps, mining, rock mechanics, mathematics and physical modelling. Fig. 1: Scheme of loading cell Photo 2). In order to embed the force meters properly and to prevent them from tilting, another double steel sheet, 1 mm thick, was placed over the force meters. A 300 mm high block of duraluminium was placed over this sheet. This block simulates the handing wall and facilitates stress distribution similar to that in reality (see Photo 1). 2.2 Force meters Photo 2 shows the force meters. The force meters are 160 mm in length, 68 mm in height, and 16 or 32 mm in width. There are 4 strain gauges on each force meter – 2 on one side 30 mm from the edge of the force meter and 2 on the other side 60 mm from the edge of the force meter. These allow us to measure the deformation along its full length. The strain gauges are connected in series in order to be able to gauge each force meter separately and at the same time to increase the gauging sensitivity. The force meters indicated by numbers 1, 2, … 30 (width 16 mm), 31, 32, … 40 (width 32 mm) were calibrated within the expected range of forces, i.e., 0 to 250 (500) kN. 2.3 Experimental verification of the measuring system The stability of the measuring system in the unloaded state was verified for approximately 26 hours. The values of the readings in the unloaded force meters were measured at 5 minute intervals. The random error of the measuring system is character- ised by the floating deviation (derived from 480 values); it approximately 0.50 kN. The errors of the gauging electronic system are only small, with a value less than approximately 0.3 % of the maximum gauged range. 3 Survey of recorded bumps Table 1 surveys the forces that acted during the recorded bumps. The evaluation of their intensity is subjective. We can conclude that • the first bumps were recorded at forces 1100–1700 kN. The lack of bumps for VW 01-03 is probably because we failed to recognise them until we had gained some experi- ence. The best way to recognize a bump is by noticing the drop in force at the press or the changes in the stress registered by the force meters. © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ 55 Acta Polytechnica Vol. 41 No. 4 – 5/2001 Photo 1: Testing device. Two camcorders recorded the test. Photo 2: Loading cell. Right with sample and force meters. Results of the tests These results are presented in the form of pi /q. An exam- ple is given in Fig. 2. Here pi is the stress registered by the i-th force meter and q is the average stress inside the loading cell. We can observe the following kinds of dependence: • For bumps which occurred when the load was low (1–1.5 MN), when there were not too many cracks and the araldite was brittle, the maximum p /q ratio was greater then 4. • The bumps that occur at higher loads exhibit a lower ratio; for 3 MN the maxima are only three times higher and over 5 MN the araldite becomes plastic and the values of pi are only slightly higher than twice that of q. This is despite the fact that the empty space has increased by 30–50 mm on each side. For six groups of stronger bumps, which occur for similar loads from 2.582 MN onward (where the bumps are listed 56 © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ Acta Polytechnica Vol. 41 No. 4 – 5/2001 ModelVW 01 02 03 04 06 07 11 12 13 14 20 21 22 Test location ZUS KI TUM TUM TUM TUM TUM TUM TUM TUM TUM TUM TUM The total force (MN) acting during the bump and eveluation of its intensity 1.500 W 1.150 W 1.400 W 1.100 W 1.600 W 1.150 W 1700 W 1.900 W 2.900 W 2.600 W 2.723 M 2.200 W 2.150 W 2.700 S 2.222 W 2.900 S 2.500 W 2.760 W 2550 W 3.500 S 3.280 S 3.189 S 3.367 S 3.100 W 3.080 M 3.400 M-S 3.169 M 3.500 M 3.700 W 3.700 W 4.080 M-S 4.147 S 4.060 M 4.300 W 4.270 M-S 4.170 VS 3.970 M 3.790 W 4.400 S 4.260 S 4.200 M 4.880 S 4.900 S 4.980 S 5.365 S 4.520 M 5.000 S 4.940 S 4.960 M-S 5.170 W 4.640 M-S 4.900 W 5.300 VS 5.500 W 5.400 M 5.540 VS 5.592 S 5.540 S 5.670 S 5.480 S 5.550 W 5.570 W 5.760 M-S 6.000 S 6.200 VS 6.470 VS 6.000 S 6.445 VS 6.000 VS 6.280 M 6.320 S 6.820 S W – WEAK M – MEDIUM M-S – MEDIUM STRONG S – STRONG VS – VERY STRONG Table 1: Total number of determined bumps Fig. 2: Stress distribution of coal seam loading before (thick line) and after bump VW 03 © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ 57 Acta Polytechnica Vol. 41 No. 4 – 5/2001 Fig. 3: Mean stress distribution of coal seam loading before (thick line) and after a bump. The force is 4170– 4265, covered bump of models VW 03, 06, 07, 11, 12, 13, 14, 20, 21, 22. � [MPa] VW 21 90–100 80–90 70–80 60–70 50–60 40–50 30–40 20–30 10–20 0–10 100 90 80 70 60 50 40 30 20 10 0 Fig. 4: Loading of the coal seam during test VW 21 in axonometry Photo 3: Detail of araldite samples at force 60 kN together with average values and corresponding model num- bers), we calculated the average values of pi/q. Fig. 3 provides an example. • The maximum pi equals approximately 2q (1.8–2.2), which is the case for all loads on a loading cell bigger than 3 MN. • The maximum for 3.2 MN occurs at the fifth force meter from the gap in the middle (80 mm), for 5.4 MN at the 8th to 12th force meter (130 –190 mm) and for a load of 6.4 MN at the 12th force meter, i.e., 190 mm from the original edge of the gap in the middle. 4 Mathematical model PFC 2D(Particle Flow Code in Two Dimensions) developed by Itasca, USA was used for the numerical modelling part of the project. A physical problem concerning the movement and interaction of circular particles may be modelled directly by PFC2D . PFC2D models the movement and interaction of circular particles by the distinct element method (DEM), as described by Cundall and Strack (1979). Particles of arbitrary shape can be created by bonding two or more particles together: these groups of particles act as autonomous objects, provided that their bond strength is high. As a limiting case, each particle may be bonded to its neighbour: the resulting assembly can be regarded as a “solid” that has elastic properties and is capable of “fractur- ing” when the bonds break in a progressive manner. PFC2D contains extensive logic to facilitate the modelling of solids as close-packed assemblies of bonded particles; the solid may be homogeneous, or it may be divided into a number of dis- crete regions of blocks. 58 Acta Polytechnica Vol. 41 No. 4 – 5/2001 Photo 4: Detail of araldite samples at force 6.5 MN. At the centre extruded material can be seen. Fig. 5: Detail of the model immediately after a bump (with ball velocities).The first cracks are cleary visible. The calculation method is a time stepping, explicit scheme. Modelling with PFC2D involves the execution of many thousands of time steps. At each step, Newton’s second law (force = mass × acceleration) is integrated twice for each particle to provide updated velocities and new positions, given a set of contact forces acting on the particle. Based on these new particle positions, contact forces are derived from the relative displacements for pairs of particles: a linear or non-linear force/displacement law at contacts may be used. Fig. 5 and Fig. 6 show details of the mathematical model after a bump with ball velocities, and Fig. 7 shows the typical stress distribution along the coal mine. The stress grows (I, II) until the first bump initiation (III). This occurs at the 5th mea- surement cycle approximately, then the stress decreases in this location, but it is increasing simultaneously at the 7th measurement cycle, where the second bump initiation occurs (IV). The subsequent bump initiations can be expected in the 8th and 10th measurement cycles (V, VI). 5 Conclusion A combination of experimental and mathematical models appears very appropriate for a study of the stress distribu- tion in a coal seam before and after bump initiation. Both methods enable a time dependent study of the problem, and make it possible to study the development of cracks during bump initiation, and thes extrusion of material into an open space during a bump. They thus offer a description of the problem that is very close to reality. Acknowledgements This research and this article have been sponsored by the Grant Agency of Czech Republic, (GAČR), grant number 103/00/0530 “Experimental and mathematical study of coal seam loading before and after a bump” (“Experimentální a teoretická studie zatížení uhelné sloje před a po otřesu”). © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ 59 Acta Polytechnica Vol. 41 No. 4 – 5/2001 Fig. 6: Detail of the model some time after a bump (with ball velocities) Fig. 7: Stress distribution along the coal mine (left side) I – 5000 cycles, II – 7000 cycles, III – 9000 cycles, IV – 12000 cycles, V – 15000 cy- cles, VI – 26000 cycles References [1] Foss, M. M., Westman, E. C.: Seismic Method for in-seam coal mine ground control problems. SEG International Exposition and 64th Annual Meeting, Los Angeles, 1994, pp. 547 – 549 [2] Goodman, R. E.: Introduction to Rock Mechanics. John Wiley & sons, 1989, p. 562 [3] Torańo, J., Rodríguez, R., Cuesta, A.: Using experimen- tal measurements in elaboration and calibration of numeri- cal models in geomechanics. Computation Methods and Experimental Measurements X, Alicante, 2001, pp. 457– 476 [4] Vacek, J., Procházka, P.: Rock bumps occurence during mining. Computation Methods and Experimental Mea- surements X, Alicante, 2001, pp. 437 – 446 [5] Vacek, J., Bouška, P.: Stress distribution in coal seam before and after bump initiation. Geotechnika 2000, Glivice- -Ustroň 2000, pp. 55 – 66 [6] Vacek, J., Procházka, P.: Behaviour of brittle rock in extreme depth. 25th Conference on Our World in Concrete & Structures, Singapore, 2000, pp. 653 – 660 [7] Williams, E. M., Westman, E. C.: Stability and Stress Evalu- ation in Mines Using In-Seam Seismic Methods. 13th Confer- ence on ground control in mining, US Bureau of Mines, 1994, pp. 290 – 297 [8] Wood, D. M.: Soil Behaviour and Critical State Soil Me- chanics. Cambridge University Press, 1996, p. 462 Doc. Ing. Jaroslav Vacek, DrSc. phone: +420 2 2435 3549 fax: +420 2 2435 3855 e-mail: Vacek@klok.cvut.cz Czech Technical University in Prague Klokner Institute Šolínova 7, 166 08 Praha 6, Czech Republic Acta Polytechnica Vol. 41 No. 4 – 5/2001 60