Besides opening geometry, in situ stress and material properties, opening support also has significant effects on stress redistribution around a roadway. To investigate these effects of rock bolts on the stress redist...Besides opening geometry, in situ stress and material properties, opening support also has significant effects on stress redistribution around a roadway. To investigate these effects of rock bolts on the stress redistribution around a roadway, a series of numerical studies were carried out using the finite difference method. Since the stress changes around a roadway caused by rock bolting is small relative to the in situ stress, they cannot obviously be observed in stress contour plots. To overcome this difficulty, a new result processing methodology was developed using the contouring program Surfer. With this methodology, the effects of rock bolts on stress redistribution can obviously be analyzed. Numerical results show that in the three patterns of rock bolts installed in the roof, in the roof and the two lateral sides, and in all the four sides of the rectangular roadway, the maximum stress magnitude of the increase is 0.931 MPa, 2.46 MPa,and 6.5 MPa, respectively; the bolt number of 5 can form an integrated ground arch; the appropriate length and pre-tensioned force of the rock bolt is 2.0 m and 60 k N, respectively. What is more, the ground arch action under the function of rock bolting is able to be effectively examined. The rock bolts dramatically increase the minor principal stress around a roadway which results in significant increase in material strength. Consequently, the major principal stress that the material can carry will greatly increase.With adequate supports, an integrated ground arch which is critical for the stability of roadway will be formed around the roadway.展开更多
The purpose of this study was to investigate the effect of bolt profile on load transfer mechanism of fully grouted bolts in jointed rocks using analytical and numerical methods. Based on the analytical method with de...The purpose of this study was to investigate the effect of bolt profile on load transfer mechanism of fully grouted bolts in jointed rocks using analytical and numerical methods. Based on the analytical method with development of methods, a new model is presented. To validate the analytical model, five different profiles modeled by ANSYS software. The profile of rock bolts T3 and T4with load transfer capacity,respectively 180 and 195 kN in the jointed rocks was selected as the optimum profiles. Finally, the selected profiles were examined in Tabas Coal Mine. FLAC analysis indicates that patterns 6+7 with2 NO flexi bolt 4 m better than other patterns within the faulted zone.展开更多
Based on the safety coefficient method,which assigns rock failure criteria to calculate the rock mass unit,the safety coefficient contour of surrounding rock is plotted to judge the distribution form of the fractured ...Based on the safety coefficient method,which assigns rock failure criteria to calculate the rock mass unit,the safety coefficient contour of surrounding rock is plotted to judge the distribution form of the fractured zone in the roadway.This will provide the basis numerical simulation to calculate the surrounding rock fractured zone in a roadway.Using the single factor and multi-factor orthogonal test method,the evolution law of roadway surrounding rock displacements,plastic zone and stress distribution under different conditions is studied.It reveals the roadway surrounding rock burst evolution process,and obtains five kinds of failure modes in deep soft rock roadway.Using the fuzzy mathematics clustering analysis method,the deep soft surrounding rock failure model in Zhujixi mine can be classified and patterns recognized.Compared to the identification results and the results detected by geological radar of surrounding rock loose circle,the reliability of the results of the pattern recognition is verified and lays the foundations for the support design of deep soft rock roadways.展开更多
A complete rock failure process usually involves opening/sliding of preexisting discontinuities as well as frac- turing in intact rock bridges to form persistent failure sur- faces and subsequent motions of the genera...A complete rock failure process usually involves opening/sliding of preexisting discontinuities as well as frac- turing in intact rock bridges to form persistent failure sur- faces and subsequent motions of the generated rock blocks. The recently developed numerical manifold method （NMM） has potential for modelling such a complete failure process. However, the NMM suffers one limitation, i.e., unexpected material domain area change occurs in rotation modelling. This problem can not be easily solved because the rigid body rotation is not represented explicitly in the NMM. The discontinuous deformation analysis （DDA） is specially de- veloped for modelling discrete block systems. The rotation- induced material area change in the DDA modelling can be avoided conveniently because the rigid body rotation is represented in an explicit form. In this paper, a transition technique is proposed and implemented to convert a NMMmodelling to a DDA modelling so as to simulate a complete rock failure process entirely by means of the two methods, in which the NMM is adopted to model the early fracturing as well as the transition from continua to discontinua, while the DDA is adopted to model the subsequent motion of the generated rock blocks. Such a numerical approach also im- proves the simulation efficiency greatly as compared with a complete NMM modelling approach. The fracturing of a rock slab with pre-existing non-persistent joints located on a slope crest and the induced rockfall process are simulated. The validity of the modelling transition from the NMM to the DDA is verified and the applicability of the proposed nu- merical approach is investigated.展开更多
The goal of this review paper is to provide a summary of selected discrete element and hybrid finitediscrete element modeling techniques that have emerged in the field of rock mechanics as simulation tools for fractur...The goal of this review paper is to provide a summary of selected discrete element and hybrid finitediscrete element modeling techniques that have emerged in the field of rock mechanics as simulation tools for fracturing processes in rocks and rock masses. The fundamental principles of each computer code are illustrated with particular emphasis on the approach specifically adopted to simulate fracture nucleation and propagation and to account for the presence of rock mass discontinuities. This description is accompanied by a brief review of application studies focusing on laboratory-scale models of rock failure processes and on the simulation of damage development around underground excavations.展开更多
Anisotropy of the strength and deformation behaviors of fractured rock masses is a crucial issue for design and stability assessments of rock engineering structures, due mainly to the non-uniform and non- regular geom...Anisotropy of the strength and deformation behaviors of fractured rock masses is a crucial issue for design and stability assessments of rock engineering structures, due mainly to the non-uniform and non- regular geometries of the fracture systems. However, no adequate efforts have been made to study this issue due to the current practical impossibility of laboratory tests with samples of large volumes con- taining many fractures, and the difficulty for controlling reliable initial and boundary conditions for large-scale in situ tests. Therefore, a reliable numerical predicting approach for evaluating anisotropy of fractured rock masses is needed. The objective of this study is to systematically investigate anisotropy of strength and deformability of fractured rocks, which has not been conducted in the past, using a nu- merical modeling method. A series of realistic two-dimensional （2D） discrete fracture network （DFN） models were established based on site investigation data, which were then loaded in different directions, using the code UDEC of discrete element method （DEM）, with changing confining pressures. Numerical results show that strength envelopes and elastic deformability parameters of tested numerical models are significantly anisotropic, and vary with changing axial loading and confining pressures. The results indicate that for design and safety assessments of rock engineering projects, the directional variations of strength and deformability of the fractured rock mass concerned must be treated properly with respect to the directions of in situ stresses. Traditional practice for simply positioning axial orientation of tunnels in association with principal stress directions only may not be adequate for safety requirements. Outstanding issues of the present study and su^zestions for future study are also oresented.展开更多
This paper presents three case studies comprising failure mechanisms in phyllite mine slopes at Quadrila- tero Ferrifero, State of Minas Gerais, Brazil. Numerical modeling techniques were used in this study. Fail- ure...This paper presents three case studies comprising failure mechanisms in phyllite mine slopes at Quadrila- tero Ferrifero, State of Minas Gerais, Brazil. Numerical modeling techniques were used in this study. Fail- ure mechanisms involving discontinuities sub parallel to the main foliation are very common in these mines. Besides, failure through the rock material has also been observed due to the low strength of phyl- lites in this site. Results of this work permitted to establish unknown geotechnical parameters which have significant influence in failure processes, like the in situ stress field and the discontinuity stiffness.展开更多
In order to control asymmetric floor heave in deep rock roadways and deformation around the surrounding rock mass after excavation, in this paper we discuss the failure mechanism and coupling control countermeasures u...In order to control asymmetric floor heave in deep rock roadways and deformation around the surrounding rock mass after excavation, in this paper we discuss the failure mechanism and coupling control countermeasures using the finite difference method （FLAC^3D） combined with comparative analysis and typical engineering application at Xingcun coal mine, It is indicated by the analysis that the simple symmetric support systems used in the past led to destruction of the deep rock roadway from the key zone and resulted in the deformation of asymmetric floor heave in the roadway. Suitable rein- forced support countermeasures are proposed to reduce the deformation of the floor heave and the potential risk during mining. The application shows that the present support technology can he used to better environmental conditions. The countermeasures of asymmetric coupling support can not only effectively reduce the discrepancy deformation at the key area of the surrounding rock mass, hut also effectively control floor heave, which helps realize the integration of support and maintain the stability of the deep rock roadways at Xingcun coal mine.展开更多
Permafrost （perennially frozen ground） appears widely in the Golmud-Lhasa section of the Qinghai-Tibet railway and is characterized by high ground temperature （≥1℃） and massive ground ice. Under the scenarios of...Permafrost （perennially frozen ground） appears widely in the Golmud-Lhasa section of the Qinghai-Tibet railway and is characterized by high ground temperature （≥1℃） and massive ground ice. Under the scenarios of global warming and human activity, the permafrost under the railway will gradually thaw and the massive ground ice will slowly melt, resulting in some thaw settlement hazards, which mainly include longitudinal and lateral cracks, and slope failure. The crushed rock layer has a thermal semiconductor effect under the periodic fluctuation of natural air. It can be used to lower the temperature of the underlying permafrost along the Qinghai-Tibet railway, and mitigate the thaw settlement hazards of the subgrade. In the present paper, the daily and annual changes in the thermal characteristics of the embankment with crushed rock side slope （ECRSS） were quantitatively simulated using the numerical method to study the cooling effect of the crushed rock layer and its mitigative ability. The results showed that the ECRSS absorbed some heat in the daytime in summer, but part of it was released at night, which accounted for approximately 20% of that absorbed. Within a year, it removed more heat from the railway subgrade in winter than that absorbed in summer. It can store approximately 20% of the ＂cold＂ energy in subgrade. Therefore, ECRSS is a better measure to mitigate thaw settlement hazards to the railway.展开更多
The spacing–burden(S/B) ratio plays significant role on rock fragmentation and proper utilization of explosive energy to minimize the undesirable damage.Low S/B ratio generates fine fragments due to pressure rings co...The spacing–burden(S/B) ratio plays significant role on rock fragmentation and proper utilization of explosive energy to minimize the undesirable damage.Low S/B ratio generates fine fragments due to pressure rings coalescence of two blast holes,whereas boulder generations were observed above optimum S/B ratio.Both conditions are not acceptable because of wastage of explosive energy.Therefore,to resolve this issue,a numerical model study was conducted to optimize the S/B ratio and to envisage its effect on rock fragmentation based on utilization of explosive energy.Finite element simulation tool was used to see the extent of two blast hole influence area variation with varying S/B ratio.The better results were obtained at S/B ratio of 1:2 with optimum utilization of peak explosive energy.The performance was observed based on peak kinetic energy,peak pressure,radial and hoop stresses on centre of the two blast holes,where pressure rings coalescence.展开更多
The island coal face arises in coal mines with the purpose of preventing gas explosion or maintaining the balance between mining and tunneling. However, its particular stress conditions in the surrounding rock may inc...The island coal face arises in coal mines with the purpose of preventing gas explosion or maintaining the balance between mining and tunneling. However, its particular stress conditions in the surrounding rock may increase the difficulty of stress control in the coal face and in its mining roadways, especially when the coal seam, the roof, and the floor have rock-burst propensities, The high energy accumulated in the island coal face and in its roof and floor will intensify rock-burst propensity or even induce rock burst, which further result in great casualties and financial losses. Taking island coal face 2321 in Jinqiao coal mine as a case, we propose a method for the prediction of rock-burst-threatened areas in an island coal face with weak rock-burst propensity. Based on the anaHysis of the movement of the overlying roof and characteristics of stress distribution, this method combined numerical simulation with drilling bits to ensure the prediction accuracy. The effects of coal pillars with different widths on the mitigation of stress concentration in the coal face and on the prevention of rock burst are analyzed together with the mech- anism behind. Finally, corresponding measures against the rock burst in the island coal face are proposed.展开更多
A numerical model based on a Continuum-based Distinct Element Method(CDEM) was used to carry out a dynamic simulation of the interaction between shield and rock strata movement in longwall mining. In Northern China, t...A numerical model based on a Continuum-based Distinct Element Method(CDEM) was used to carry out a dynamic simulation of the interaction between shield and rock strata movement in longwall mining. In Northern China, the Ordos coal field geological conditions and operational characteristics were used as a case example. The CDEM was constructed on Ordos coal field shield's operation characteristics and geological conditions. Numerical modelling was carried out to investigate the effects of different mining heights on the caving process, movement characteristics, equilibrium and stability conditions of overburden as the interaction between shield and surrounding rocks. With the numerical model, the internal factors for changes in shield resistance under different mining heights was found. The quantitative relationship between mining heights and shield resistance was also obtained by the numerical simulation.展开更多
To study the law of gas dilatation energy release of rock cross-cut coal uncovering face, according to the analysis of the physical parameters distribution features of coal and rock mass in front of crosscut face,the ...To study the law of gas dilatation energy release of rock cross-cut coal uncovering face, according to the analysis of the physical parameters distribution features of coal and rock mass in front of crosscut face,the equations of elastic potential of coal and gas dilatation energy theory were set up to process a contrast calculation of the sizes of two kinds of energy. The results show that gas dilatation energy is the uppermost energy source causing outburst occurrence. Furthermore, the mathematical model of spherical flow field gas dilatation energy release was established and MATLAB software was applied to make a numerical calculation analysis on the law of gas dilatation energy release. The results indicate that the gas dilatation energy is closely related to gas parameters and its energy index does reflect the possibility of coal seam outburst.展开更多
As is known, high-level radioactive waste （HLW） is commonly heat-emitting. Heat output from HLWwilldissipate through the surrounding rocks and induce complex thermo-hydro-mechanical-chemical（THMC） processes. In hi...As is known, high-level radioactive waste （HLW） is commonly heat-emitting. Heat output from HLWwilldissipate through the surrounding rocks and induce complex thermo-hydro-mechanical-chemical（THMC） processes. In highly consolidated clayey rocks, thermal effects are particularly significantbecause of their very low permeability and water-saturated state. Thermal impact on the integrity of thegeological barriers is of most importance with regard to the long-term safety of repositories. This studyfocuses on numerical analysis of thermal effects on hydro-mechanical properties of clayey rock using acoupled thermo-mechanical multiphase flow （TH2M） model which is implemented in the finite elementprogramme OpenGeoSys （OGS）. The material properties of the numerical model are characterised by atransversal isotropic elastic model based on Hooke＇s law, a non-isothermal multiphase flow model basedon van Genuchten function and Darcy＇s law, and a transversal isotropic heat transport model based onFourier＇s law. In the numerical approaches, special attention has been paid to the thermal expansion ofthree different phases： gas, fluid and solid, which could induce changes in pore pressure and porosity.Furthermore, the strong swelling and shrinkage behaviours of clayey material are also considered in thepresent model. The model has been applied to simulate a laboratory heating experiment on claystone.The numerical model gives a satisfactory representation of the observed material behaviour in thelaboratory experiment. The comparison of the calculated results with the laboratory findings verifies thatthe simulation with the present numerical model could provide a deeper understanding of the observedeffects. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.展开更多
Large underground caverns are commonly used in variety of applications. In many cases, because of the geomechanical limitations of dimensions and requirement of high volume, several parallel caverns are used. Plastic ...Large underground caverns are commonly used in variety of applications. In many cases, because of the geomechanical limitations of dimensions and requirement of high volume, several parallel caverns are used. Plastic zone integration requires a larger rock pillar distance of theses adjacent caverns while eco- nomic and access reasons require smaller distance. In lran many underground projects are located in West and South West, Asmari formation covers a large part of these regions. The stability of underground spaces that are constructed or will be constructed in this formation has been investigated. A proper cross section based on plastic analysis and a stability criterion has been proposed for each region. Finally, in each case, allowable rock pillar between adjacent caverns with similar dimension was determined with two methods （numerical analysis and fire service law）. Results show that Fire Service Law uses a very con- servative safety factor and it was proposed to use a correction factor for allowable distance based on application of underground space.展开更多
It is important to investigate the dynamic behaviors of deep rocks near explosion cavity to reveal the mechanisms of deformations and fractures. Some improvements are carried out for Grigorian model with focuses on th...It is important to investigate the dynamic behaviors of deep rocks near explosion cavity to reveal the mechanisms of deformations and fractures. Some improvements are carried out for Grigorian model with focuses on the dilation effects and the relaxation effects of deep rocks, and the high pressure equations of states with Mie-Grüneisen form are also established. Numerical calculations of free field parameters for deep underground explosions are carried out based on the user subroutines which are compiled by means of the secondary development functions of LS-DYNA9703 D software. The histories of radial stress, radial velocity and radial displacement of rock particles are obtained, and the calculation results are compared with those of U.S. Hardhat nuclear test. It is indicated that the dynamic responses of free field for deep underground explosions are well simulated based on improved Grigorian model, and the calculation results are in good agreement with the data of U.S. Hardhat nuclear test. The peak values of particle velocities are consistent with those of test, but the waveform widths and the rising times are obviously greater than those without dilation effects. The attenuation rates of particle velocities are greater than the calculation results with classic plastic model, and they are consistent with the results of Hardhat nuclear test. The attenuation behaviors and the rising times of stress waves are well shown by introducing dilation effects and relaxation effects into the calculation model. Therefore, the defects of Grigorian model are avoided. It is also indicated that the initial stress has obvious influences on the waveforms of radial stress and the radial displacements of rock particles.展开更多
In order to better understand the mechanical properties of graded crushed rocks （GCRs） and to optimize the relevant design, a numerical test method based on the particle flow modeling technique PFC2D is developed fo...In order to better understand the mechanical properties of graded crushed rocks （GCRs） and to optimize the relevant design, a numerical test method based on the particle flow modeling technique PFC2D is developed for the California bearing ratio （CBR） test on GGRs. The effects of different testing conditions and micro-mechanical parameters used in the model on the CBR numerical results have been systematically studied. The reliability of the numerical technique is verified. The numerical results suggest that the influences of the loading rate and Poisson＇s ratio on the CBR numerical test results are not significant. As such, a loading rate of 1.0-3.0 mm/min, a piston diameter of 5 cm, a specimen height of 15 cm and a specimen diameter of 15 cm are adopted for the CBR numerical test. The numerical results reveal that the GBR values increase with the friction coefficient at the contact and shear modulus of the rocks, while the influence of Poisson＇s ratio on the GBR values is insignificant. The close agreement between the CBR numerical results and experimental results suggests that the numerical simulation of the CBR values is promising to help assess the mechanical properties of GGRs and to optimize the grading design. Be- sides, the numerical study can provide useful insights on the mesoscopic mechanism.展开更多
基金Financial support for this work provided by the National Key Scientific Apparatus Development of Special Item (No.2012YQ24012705)is deeply appreciated
文摘Besides opening geometry, in situ stress and material properties, opening support also has significant effects on stress redistribution around a roadway. To investigate these effects of rock bolts on the stress redistribution around a roadway, a series of numerical studies were carried out using the finite difference method. Since the stress changes around a roadway caused by rock bolting is small relative to the in situ stress, they cannot obviously be observed in stress contour plots. To overcome this difficulty, a new result processing methodology was developed using the contouring program Surfer. With this methodology, the effects of rock bolts on stress redistribution can obviously be analyzed. Numerical results show that in the three patterns of rock bolts installed in the roof, in the roof and the two lateral sides, and in all the four sides of the rectangular roadway, the maximum stress magnitude of the increase is 0.931 MPa, 2.46 MPa,and 6.5 MPa, respectively; the bolt number of 5 can form an integrated ground arch; the appropriate length and pre-tensioned force of the rock bolt is 2.0 m and 60 k N, respectively. What is more, the ground arch action under the function of rock bolting is able to be effectively examined. The rock bolts dramatically increase the minor principal stress around a roadway which results in significant increase in material strength. Consequently, the major principal stress that the material can carry will greatly increase.With adequate supports, an integrated ground arch which is critical for the stability of roadway will be formed around the roadway.
文摘The purpose of this study was to investigate the effect of bolt profile on load transfer mechanism of fully grouted bolts in jointed rocks using analytical and numerical methods. Based on the analytical method with development of methods, a new model is presented. To validate the analytical model, five different profiles modeled by ANSYS software. The profile of rock bolts T3 and T4with load transfer capacity,respectively 180 and 195 kN in the jointed rocks was selected as the optimum profiles. Finally, the selected profiles were examined in Tabas Coal Mine. FLAC analysis indicates that patterns 6+7 with2 NO flexi bolt 4 m better than other patterns within the faulted zone.
基金provided by the National Natural Science Foundation of China(Nos.51322401,51309222,51323004,51579239 and 51574223)the Opening Project Fund of Shandong Provincial Key Laboratory of Civil Engineering Disaster Prevention and Mitigation(No.CDPM2014KF03)+2 种基金the State Key Laboratory for GeoMechanics Opening Project Fund of Shandong Provincial Key Laboratory of Civil Engineering Disaster Prevention and MitigationDeep Underground Engineering,China University of Mining&Technology(No.SKLGDUEK1305)China Postdoctoral Science Foundation(Nos.2014M551700and 2013M531424)
文摘Based on the safety coefficient method,which assigns rock failure criteria to calculate the rock mass unit,the safety coefficient contour of surrounding rock is plotted to judge the distribution form of the fractured zone in the roadway.This will provide the basis numerical simulation to calculate the surrounding rock fractured zone in a roadway.Using the single factor and multi-factor orthogonal test method,the evolution law of roadway surrounding rock displacements,plastic zone and stress distribution under different conditions is studied.It reveals the roadway surrounding rock burst evolution process,and obtains five kinds of failure modes in deep soft rock roadway.Using the fuzzy mathematics clustering analysis method,the deep soft surrounding rock failure model in Zhujixi mine can be classified and patterns recognized.Compared to the identification results and the results detected by geological radar of surrounding rock loose circle,the reliability of the results of the pattern recognition is verified and lays the foundations for the support design of deep soft rock roadways.
基金supported by the Research Fund for the Doctoral Program of Higher Education of China (20090101120057)the Scientific Research Fund of Zhejiang Provincial Education Department (Y200909163)
文摘A complete rock failure process usually involves opening/sliding of preexisting discontinuities as well as frac- turing in intact rock bridges to form persistent failure sur- faces and subsequent motions of the generated rock blocks. The recently developed numerical manifold method （NMM） has potential for modelling such a complete failure process. However, the NMM suffers one limitation, i.e., unexpected material domain area change occurs in rotation modelling. This problem can not be easily solved because the rigid body rotation is not represented explicitly in the NMM. The discontinuous deformation analysis （DDA） is specially de- veloped for modelling discrete block systems. The rotation- induced material area change in the DDA modelling can be avoided conveniently because the rigid body rotation is represented in an explicit form. In this paper, a transition technique is proposed and implemented to convert a NMMmodelling to a DDA modelling so as to simulate a complete rock failure process entirely by means of the two methods, in which the NMM is adopted to model the early fracturing as well as the transition from continua to discontinua, while the DDA is adopted to model the subsequent motion of the generated rock blocks. Such a numerical approach also im- proves the simulation efficiency greatly as compared with a complete NMM modelling approach. The fracturing of a rock slab with pre-existing non-persistent joints located on a slope crest and the induced rockfall process are simulated. The validity of the modelling transition from the NMM to the DDA is verified and the applicability of the proposed nu- merical approach is investigated.
文摘The goal of this review paper is to provide a summary of selected discrete element and hybrid finitediscrete element modeling techniques that have emerged in the field of rock mechanics as simulation tools for fracturing processes in rocks and rock masses. The fundamental principles of each computer code are illustrated with particular emphasis on the approach specifically adopted to simulate fracture nucleation and propagation and to account for the presence of rock mass discontinuities. This description is accompanied by a brief review of application studies focusing on laboratory-scale models of rock failure processes and on the simulation of damage development around underground excavations.
文摘Anisotropy of the strength and deformation behaviors of fractured rock masses is a crucial issue for design and stability assessments of rock engineering structures, due mainly to the non-uniform and non- regular geometries of the fracture systems. However, no adequate efforts have been made to study this issue due to the current practical impossibility of laboratory tests with samples of large volumes con- taining many fractures, and the difficulty for controlling reliable initial and boundary conditions for large-scale in situ tests. Therefore, a reliable numerical predicting approach for evaluating anisotropy of fractured rock masses is needed. The objective of this study is to systematically investigate anisotropy of strength and deformability of fractured rocks, which has not been conducted in the past, using a nu- merical modeling method. A series of realistic two-dimensional （2D） discrete fracture network （DFN） models were established based on site investigation data, which were then loaded in different directions, using the code UDEC of discrete element method （DEM）, with changing confining pressures. Numerical results show that strength envelopes and elastic deformability parameters of tested numerical models are significantly anisotropic, and vary with changing axial loading and confining pressures. The results indicate that for design and safety assessments of rock engineering projects, the directional variations of strength and deformability of the fractured rock mass concerned must be treated properly with respect to the directions of in situ stresses. Traditional practice for simply positioning axial orientation of tunnels in association with principal stress directions only may not be adequate for safety requirements. Outstanding issues of the present study and su^zestions for future study are also oresented.
文摘This paper presents three case studies comprising failure mechanisms in phyllite mine slopes at Quadrila- tero Ferrifero, State of Minas Gerais, Brazil. Numerical modeling techniques were used in this study. Fail- ure mechanisms involving discontinuities sub parallel to the main foliation are very common in these mines. Besides, failure through the rock material has also been observed due to the low strength of phyl- lites in this site. Results of this work permitted to establish unknown geotechnical parameters which have significant influence in failure processes, like the in situ stress field and the discontinuity stiffness.
基金support from the National Natural Science Foundation of China (Nos. 51134005, 51374214, 41172116, and U1261212)the New Century Excellent Talents Foundation in University (No. NCET-07-0800)the Special Fund of Basic Research and Operating of China University of Mining & Technology in Beijing (No. 2009QL03)
文摘In order to control asymmetric floor heave in deep rock roadways and deformation around the surrounding rock mass after excavation, in this paper we discuss the failure mechanism and coupling control countermeasures using the finite difference method （FLAC^3D） combined with comparative analysis and typical engineering application at Xingcun coal mine, It is indicated by the analysis that the simple symmetric support systems used in the past led to destruction of the deep rock roadway from the key zone and resulted in the deformation of asymmetric floor heave in the roadway. Suitable rein- forced support countermeasures are proposed to reduce the deformation of the floor heave and the potential risk during mining. The application shows that the present support technology can he used to better environmental conditions. The countermeasures of asymmetric coupling support can not only effectively reduce the discrepancy deformation at the key area of the surrounding rock mass, hut also effectively control floor heave, which helps realize the integration of support and maintain the stability of the deep rock roadways at Xingcun coal mine.
基金supported by the National Natural Science Foundation of China(Grant Nos 40801022 and 40821001)the Chinese Academy of Sciences(CAS) Knowledge Innovation Key Directional Program(Grant Nos KZCX2-YW-Q03-04 and KZCX2-YW-311)+2 种基金CAS Western Project(Grant No.KZCX2-XB2-10)China Postdoctoral Science Foundation(No.20080430110)CAS 100-Talent Programs"Stability of Linear Engineering Foundations in Warm Permafrost Regions under a Changing Climate"and"Deformation and Stability of Roadbed in Permafrost Regions",CAS West Light Foundation for PhD G.Y.Li,Funds of the State Key Laboratory of Frozen Soils Engineering,CAS(Grant Nos SKLFSE-ZQ-02 and SKLFSE-ZY-03)
文摘Permafrost （perennially frozen ground） appears widely in the Golmud-Lhasa section of the Qinghai-Tibet railway and is characterized by high ground temperature （≥1℃） and massive ground ice. Under the scenarios of global warming and human activity, the permafrost under the railway will gradually thaw and the massive ground ice will slowly melt, resulting in some thaw settlement hazards, which mainly include longitudinal and lateral cracks, and slope failure. The crushed rock layer has a thermal semiconductor effect under the periodic fluctuation of natural air. It can be used to lower the temperature of the underlying permafrost along the Qinghai-Tibet railway, and mitigate the thaw settlement hazards of the subgrade. In the present paper, the daily and annual changes in the thermal characteristics of the embankment with crushed rock side slope （ECRSS） were quantitatively simulated using the numerical method to study the cooling effect of the crushed rock layer and its mitigative ability. The results showed that the ECRSS absorbed some heat in the daytime in summer, but part of it was released at night, which accounted for approximately 20% of that absorbed. Within a year, it removed more heat from the railway subgrade in winter than that absorbed in summer. It can store approximately 20% of the ＂cold＂ energy in subgrade. Therefore, ECRSS is a better measure to mitigate thaw settlement hazards to the railway.
文摘The spacing–burden(S/B) ratio plays significant role on rock fragmentation and proper utilization of explosive energy to minimize the undesirable damage.Low S/B ratio generates fine fragments due to pressure rings coalescence of two blast holes,whereas boulder generations were observed above optimum S/B ratio.Both conditions are not acceptable because of wastage of explosive energy.Therefore,to resolve this issue,a numerical model study was conducted to optimize the S/B ratio and to envisage its effect on rock fragmentation based on utilization of explosive energy.Finite element simulation tool was used to see the extent of two blast hole influence area variation with varying S/B ratio.The better results were obtained at S/B ratio of 1:2 with optimum utilization of peak explosive energy.The performance was observed based on peak kinetic energy,peak pressure,radial and hoop stresses on centre of the two blast holes,where pressure rings coalescence.
基金provided by the National Natural Science Foundation of China (Nos.51304208 and 51474208)
文摘The island coal face arises in coal mines with the purpose of preventing gas explosion or maintaining the balance between mining and tunneling. However, its particular stress conditions in the surrounding rock may increase the difficulty of stress control in the coal face and in its mining roadways, especially when the coal seam, the roof, and the floor have rock-burst propensities, The high energy accumulated in the island coal face and in its roof and floor will intensify rock-burst propensity or even induce rock burst, which further result in great casualties and financial losses. Taking island coal face 2321 in Jinqiao coal mine as a case, we propose a method for the prediction of rock-burst-threatened areas in an island coal face with weak rock-burst propensity. Based on the anaHysis of the movement of the overlying roof and characteristics of stress distribution, this method combined numerical simulation with drilling bits to ensure the prediction accuracy. The effects of coal pillars with different widths on the mitigation of stress concentration in the coal face and on the prevention of rock burst are analyzed together with the mech- anism behind. Finally, corresponding measures against the rock burst in the island coal face are proposed.
基金funded by the National Natural Science Foundation of China (Nos. U1261207 and 51274086)
文摘A numerical model based on a Continuum-based Distinct Element Method(CDEM) was used to carry out a dynamic simulation of the interaction between shield and rock strata movement in longwall mining. In Northern China, the Ordos coal field geological conditions and operational characteristics were used as a case example. The CDEM was constructed on Ordos coal field shield's operation characteristics and geological conditions. Numerical modelling was carried out to investigate the effects of different mining heights on the caving process, movement characteristics, equilibrium and stability conditions of overburden as the interaction between shield and surrounding rocks. With the numerical model, the internal factors for changes in shield resistance under different mining heights was found. The quantitative relationship between mining heights and shield resistance was also obtained by the numerical simulation.
文摘To study the law of gas dilatation energy release of rock cross-cut coal uncovering face, according to the analysis of the physical parameters distribution features of coal and rock mass in front of crosscut face,the equations of elastic potential of coal and gas dilatation energy theory were set up to process a contrast calculation of the sizes of two kinds of energy. The results show that gas dilatation energy is the uppermost energy source causing outburst occurrence. Furthermore, the mathematical model of spherical flow field gas dilatation energy release was established and MATLAB software was applied to make a numerical calculation analysis on the law of gas dilatation energy release. The results indicate that the gas dilatation energy is closely related to gas parameters and its energy index does reflect the possibility of coal seam outburst.
基金supported by BMWi (Bundesministerium für Wirtschaft und Energie,Berlin)
文摘As is known, high-level radioactive waste （HLW） is commonly heat-emitting. Heat output from HLWwilldissipate through the surrounding rocks and induce complex thermo-hydro-mechanical-chemical（THMC） processes. In highly consolidated clayey rocks, thermal effects are particularly significantbecause of their very low permeability and water-saturated state. Thermal impact on the integrity of thegeological barriers is of most importance with regard to the long-term safety of repositories. This studyfocuses on numerical analysis of thermal effects on hydro-mechanical properties of clayey rock using acoupled thermo-mechanical multiphase flow （TH2M） model which is implemented in the finite elementprogramme OpenGeoSys （OGS）. The material properties of the numerical model are characterised by atransversal isotropic elastic model based on Hooke＇s law, a non-isothermal multiphase flow model basedon van Genuchten function and Darcy＇s law, and a transversal isotropic heat transport model based onFourier＇s law. In the numerical approaches, special attention has been paid to the thermal expansion ofthree different phases： gas, fluid and solid, which could induce changes in pore pressure and porosity.Furthermore, the strong swelling and shrinkage behaviours of clayey material are also considered in thepresent model. The model has been applied to simulate a laboratory heating experiment on claystone.The numerical model gives a satisfactory representation of the observed material behaviour in thelaboratory experiment. The comparison of the calculated results with the laboratory findings verifies thatthe simulation with the present numerical model could provide a deeper understanding of the observedeffects. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.
文摘Large underground caverns are commonly used in variety of applications. In many cases, because of the geomechanical limitations of dimensions and requirement of high volume, several parallel caverns are used. Plastic zone integration requires a larger rock pillar distance of theses adjacent caverns while eco- nomic and access reasons require smaller distance. In lran many underground projects are located in West and South West, Asmari formation covers a large part of these regions. The stability of underground spaces that are constructed or will be constructed in this formation has been investigated. A proper cross section based on plastic analysis and a stability criterion has been proposed for each region. Finally, in each case, allowable rock pillar between adjacent caverns with similar dimension was determined with two methods （numerical analysis and fire service law）. Results show that Fire Service Law uses a very con- servative safety factor and it was proposed to use a correction factor for allowable distance based on application of underground space.
基金Project(51378498)supported by the National Natural Science Foundation of ChinaProject(BK20141066)supported the Natural Science Foundation of Jiangsu Province,China+1 种基金Project(SKLGDUEK1208)supported by State Key Laboratory for Geo Mechanics and Deep Underground Engineering(China University of Mining & Technology),ChinaProject(DPMEIKF201301)supported by State Key Laboratory of Disaster Prevention & Mitigation of Explosion & Impact(PLA University of Science and Technology),China
文摘It is important to investigate the dynamic behaviors of deep rocks near explosion cavity to reveal the mechanisms of deformations and fractures. Some improvements are carried out for Grigorian model with focuses on the dilation effects and the relaxation effects of deep rocks, and the high pressure equations of states with Mie-Grüneisen form are also established. Numerical calculations of free field parameters for deep underground explosions are carried out based on the user subroutines which are compiled by means of the secondary development functions of LS-DYNA9703 D software. The histories of radial stress, radial velocity and radial displacement of rock particles are obtained, and the calculation results are compared with those of U.S. Hardhat nuclear test. It is indicated that the dynamic responses of free field for deep underground explosions are well simulated based on improved Grigorian model, and the calculation results are in good agreement with the data of U.S. Hardhat nuclear test. The peak values of particle velocities are consistent with those of test, but the waveform widths and the rising times are obviously greater than those without dilation effects. The attenuation rates of particle velocities are greater than the calculation results with classic plastic model, and they are consistent with the results of Hardhat nuclear test. The attenuation behaviors and the rising times of stress waves are well shown by introducing dilation effects and relaxation effects into the calculation model. Therefore, the defects of Grigorian model are avoided. It is also indicated that the initial stress has obvious influences on the waveforms of radial stress and the radial displacements of rock particles.
基金supported by the Program for New Century Excellent Talents in University (NCET-08-0749)Fundamental Research Funds for the Central Universities (CHD2012JC054)
文摘In order to better understand the mechanical properties of graded crushed rocks （GCRs） and to optimize the relevant design, a numerical test method based on the particle flow modeling technique PFC2D is developed for the California bearing ratio （CBR） test on GGRs. The effects of different testing conditions and micro-mechanical parameters used in the model on the CBR numerical results have been systematically studied. The reliability of the numerical technique is verified. The numerical results suggest that the influences of the loading rate and Poisson＇s ratio on the CBR numerical test results are not significant. As such, a loading rate of 1.0-3.0 mm/min, a piston diameter of 5 cm, a specimen height of 15 cm and a specimen diameter of 15 cm are adopted for the CBR numerical test. The numerical results reveal that the GBR values increase with the friction coefficient at the contact and shear modulus of the rocks, while the influence of Poisson＇s ratio on the GBR values is insignificant. The close agreement between the CBR numerical results and experimental results suggests that the numerical simulation of the CBR values is promising to help assess the mechanical properties of GGRs and to optimize the grading design. Be- sides, the numerical study can provide useful insights on the mesoscopic mechanism.
