Numerical modeling for the coupled thermo-mechanical processes and spalling phenomena in sp Pillar Stability Experiment (APSE)

In this paper, the coupled thermo-mechanical （TM） processes in the AEspoe Pillar Stability Experiment （APSE） carried out by the Swedish Nuclear Fuel and Waste Management Company （SKB） were simulated using both continuum and discontinuum based numerical methods. Two-dimensional （2D） and three- dimensional （3D） finite element method （FEM） and 2D distinct element method （DEM） with particles were used. The main objective for the large scale in situ experiment is to investigate the yielding strength of crystalline rock and the formation of the excavation disturbed/damaged zone （EDZ） during excavation of two boreholes, pressurizing of one of the boreholes and heating. For the DEM simulations, the heat flow algorithm was newly introduced into the original code. The calculated stress, displacement and temperature distributions were compared with the ones obtained from in situ measurements and FEM simulations. A parametric study for initial microcracks was also performed to reproduce the spalling phenomena observed in the APSE.

Numerical study on coupled thermo-mechanical processes in sp Pillar Stability Experiment

This paper presents a study of the full three-dimensional thermo-mechanical （TM） behavior of rock pillar in,Aspo Pillar Stability Experiment （APSE） using a self-developed numerical code TM-EPCA3D. The transient thermal conduction function was descritized on space and time scales, and was solved by using cellular automaton （CA） method on space scale and finite difference method on time scale, respectively. The advantage of this approach is that no global, but local matrix is used so that it avoids the need to develop and solve large-scale linear equations and the complexity therein. A thermal conductivity versus stress function was proposed to reflect the effect of stress on thermal field. The temperature evolution and induced thermal stress in the pillar part during the heating and cooling processes were well simulated by the developed code. The factors that affect the modeling results were discussed. It is concluded that, the complex TM behavior of Aspo rock pillar is significantly influenced by the complex boundary and initial conditions.

Analysis of sp Pillar Stability Experiment: Continuous thermo-mechanical model development and calibration

The paper describes an analysis of thermo-mechanical （TM） processes appearing during the Aspo Pillar Stability Experiment （APSE）. This analysis is based on finite elements with elasticity, plasticity and dam- age mechanics models of rock behaviour and some least squares calibration techniques. The main aim is to examine the capability of continuous mechanics models to predict brittle damage behaviour of gran- ite rocks. The performed simulations use an in-house finite element software GEM and self-developed experimental continuum damage MATLAB code. The main contributions are twofold. First, it is an inverse analysis, which is used for （1） verification of an initial stress measurement by back analysis of conver- gence measurement during construction of the access tunnel and （2） identification of heat transfer rock mass properties by an inverse method based on the known heat sources and temperature measurements. Second, three different hierarchically built models are used to estimate the pillar damage zones, i.e. elas- tic model with Drucker-Prager strength criterion, elasto-plastic model with the same yield limit and a combination of elasto-plasticity with continuum damage mechanics. The damage mechanics model is also used to simulate uniaxial and triaxial compressive strength tests on the ,Aspo granite.

Stability prediction of hard rock pillar using support vector machine optimized by three metaheuristic algorithms

Hard rock pillar is one of the important structures in engineering design and excavation in underground mines.Accurate and convenient prediction of pillar stability is of great significance for underground space safety.This paper aims to develop hybrid support vector machine(SVM)models improved by three metaheuristic algorithms known as grey wolf optimizer(GWO),whale optimization algorithm(WOA)and sparrow search algorithm(SSA)for predicting the hard rock pillar stability.An integrated dataset containing 306 hard rock pillars was established to generate hybrid SVM models.Five parameters including pillar height,pillar width,ratio of pillar width to height,uniaxial compressive strength and pillar stress were set as input parameters.Two global indices,three local indices and the receiver operating characteristic(ROC)curve with the area under the ROC curve(AUC)were utilized to evaluate all hybrid models’performance.The results confirmed that the SSA-SVM model is the best prediction model with the highest values of all global indices and local indices.Nevertheless,the performance of the SSASVM model for predicting the unstable pillar(AUC:0.899)is not as good as those for stable(AUC:0.975)and failed pillars(AUC:0.990).To verify the effectiveness of the proposed models,5 field cases were investigated in a metal mine and other 5 cases were collected from several published works.The validation results indicated that the SSA-SVM model obtained a considerable accuracy,which means that the combination of SVM and metaheuristic algorithms is a feasible approach to predict the pillar stability.

Stability of coal pillar in gob-side entry driving under unstable overlying strata and its coupling support control technique

Considering the situation that it is difficult to control the stability of narrow coal pillar in gob-side entry driving under unstable overlying strata, the finite difference numerical simulation method was adopted to analyze the inner stress distribution and its evolution regularity, as well as the deformation characteristics of narrow coal pillar in gob-side entry driving, in the whole process from entry driving of last working face to the present working face mining. A new method of narrow coal pillar control based on the triune coupling support technique (TCST), which includes that high-strength prestressed thread steel bolt is used to strain the coal on the goaf side, and that short bolt to control the integrity of global displacement zone in coal pillar on the entry side, and that long grouting cable to fix anchor point to constrain the bed separation between global displacement zone and fixed zone, is thereby generated and applied to the field production. The result indicates that after entry excavating along the gob under unstable overlying strata, the supporting structure left on the gob side of narrow coal pillar is basically invalid to maintain the coal-pillar stability, and the large deformation of the pillar on the gob side is evident. Except for the significant dynamic pressure appearing in the coal mining of last working face and overlying strata stabilizing process, the stress variation inside the coal pillar in other stages are rather steady, however, the stress expansion is obvious and the coal pillar continues to deform. Once the gob-side entry driving is completed, a global displacement zone on the entry side appears in the shallow part of the pillar, whereas, a relatively steady fixed zone staying almost still in gob-side entry driving and present working face mining is found in the deep part of the pillar. The application of TCST can not only avoid the failure of pillar supporting structure, but exert the supporting capacity of the bolting structure left in the pillar of last sublevel entry, thus to jointly maintain the stability of coal pillar.

Numerical investigation into pillar failure induced by time-dependent skin degradation

This paper focuses on the instability mechanism of an isolated pillar, caused by time-dependent skin degradation and strength heterogeneity. The time-dependent skin degradation is simulated with a non-linear rheological model capable of simulating tertiary creep, whereby two different pillar failure cases are investigated. The first case is of an isolated pillar in a deep hard rock underground mine and subjected to high stresses. The results show that pillar degradation is limited to the regions near the surface or the skin until two months after ore extraction. Afterwards degradation starts to extend deeper into the pillar, eventually leaving a highly-stressed pillar core due to stress transfer from the failed skin.Rockburst potential indices show that the risk increases exponentially at the core as time goes by. It is then demonstrated that the progressive skin degradation cannot be simulated with conventional strain-softening model assuming brittle failure. The parametric study with respect to the degree of heterogeneity reveals that heterogeneity is key to the occurrence of progressive skin degradation. The second case investigated in this study is pillar failure taking place in a very long period. Such failure becomes significantly important when assessing the risk for ground subsidence caused by pillar collapse in an abandoned mine. The analysis results demonstrate that the employed non-linear rheological model can simulate gradual skin degradation taking place over several hundred years. The percentage of damage zone volume within the pillar is merely 1% after a lapse of one days and increases to 50% after one hundred years, indicating a high risk for pillar collapse in the long term. The vertical displacements within the pillar also indicate the risk of subsidence. The proposed method is suitable for evaluating the risk of ground surface subsidence above an abandoned mine.

Numerical modeling of coupled thermo-mechanical response of a rock pillar

Understanding the rock mass response to excavation and thermal loading and improving the capability of the numerical models for simulating the progressive failure process of brittle rocks are important for safety assessment and optimization design of nuclear waste repositories.The international cooperative DECOVALEX-2011 project provides a platform for development,validation and comparison of numerical models,in which the sp pillar stability experiment（APSE） was selected as the modeling target for Task B.This paper presents the modeling results of Wuhan University（WHU） team for stages 1 and 2 of Task B by using a coupled thermo-mechanical model within the framework of continuum mechanics.The rock mass response to excavation is modeled with linear elastic,elastoplastic and brittle-plastic models,while the response to heating is modeled with a coupled thermo-elastic model.The capabilities and limitations of the model for representation of the thermo-mechanical responses of the rock pillar are discussed by comparing the modeling results with experimental observations.The results may provide a helpful reference for the stability and safety assessment of the hard granite host rock in China＇s Beishan preselected area for high-level radioactive waste disposal.

Determination of the load bearing capacity of pre-stressed expandable props for ground support in underground mines

This paper aims to determine the load bearing capacity of pre-stressed expandable props with different geometries and load eccentricities for flexible support in underground mining or excavation.It is deduced that the expandable device could have much higher strength(>89 MPa)by laboratory tests,and the load bearing capacity of the expandable prop may depend on the stability of the supporting steel pipe structure.A good agreement was found between the laboratory test and numerical results in terms of the load bearing capacity and the final macro-bending failure pattern for expandable props with heights of 1.5 and 2.7 m,and the theoretical calculation for the strength of traditional steel structures is not directly suitable for the expandable props.Moreover,additional numerical simulations were performed for the expandable props with different normalized slenderness ratiosλ_(n)and loading eccentric distances e.The variation of stability coefficient of the expandable prop is in line with the Perry-Robertson equation and its correlation coefficients are fitted as a of 0.979 and b of 0.314.For estimating the load bearing capacity of the expandable props,the strength equation for traditional steel structures is improved by introducing a bending magnification factor and by modifying the normalized slenderness ratio to a converted slenderness ratio.Based on the underground field monitoring for the strength of expandable props with different heights,the empirical eccentric distances were back calculated,and a safety factor is introduced to obtain the designed strength of the expandable prop.In addition,a four-step design procedure is proposed for the expandable prop.

Thermo-mechanical coupling analysis of APSE using submodels and neural networks

The ,Aspoe Pillar Stability Experiment （APSE） is an in situ experiment for investigating the spalling mechanism under mechanical and thermal loading conditions in a crystalline rock. In this study, the thermo-mechanical behaviors in the APSE were investigated with three models： （1） a Full model with rough meshes for calculating the influence of tunnel excavation; （2） a Submodel with fine meshes for predicting the thermo-mechanical behavior in the pillar during the borehole drilling, heating, and cool- ing phases; and （3） a Thin model for modeling the effect of slot cutting for de-stressing around the pillar. In order to import the stresses calculated from the Full model to the Submodel and to define the complex thermal boundary conditions, artificial neural networks （NNs） were utilized. From this study, it was pos- sible to conclude that the stepwise approach with the application of NNs was useful for predicting the complex response of the pillar under severe thermo-mechanical loading conditions.