Search results

The interaction between rock mass structural planes and dynamic stress levels is important to determine the stability of rock mass structures in underground geotechnical engineering. In this work, the authors aim to focus on the degradation effects of fracture geometric parameters and unloading stress paths on rock mechanical properties.

A three-dimensional Particle Flow Code (PFC3D) was used for a systematic numerical simulation of the strength failure and cracking behavior of granite specimens containing prefabricated cracks under conventional triaxial compression and triaxial unilateral unloading. The authors demonstrated the unique mechanical response of prefabricated fractured rock under two conditions. The crack initiation, propagation, and coalescence process of pre-fissured specimens were analyzed in detail.

The authors show that the prefabricated cracks and unilateral unloading conditions not only deteriorate the mechanical strength but also have significant differences in failure modes. The degrading effect of cracks on model strength increases linearly with the decrease of the dip angle. Under the condition of true triaxial unilateral unloading, the deterioration effect of peak strength of rock is very significant, and unloading plays a role in promoting the instability failure of rock after peak, making the rock earlier instability failure. Associating with the particle vector diagram and crack coalescence process, the authors find that model failure mode under unilateral loading conditions is obviously distinct from that in triaxial loading. The peak strain in the unloading direction increases sharply, resulting in a new shear slip.

This study is expected to improve the understanding of the strength failure and cracking behavior of fractured rock under unilateral unloading.

Among the various applications involving the use of microwave energy, its growing utility within the mining industry is particularly noteworthy. Conventional grinding processes are often overburdened by energy inefficiencies that are directly related to machine wear, pollution and rising project costs. In this work, we numerically investigate the effects of microwave pretreatment through a series of compression tests as a means to help mitigate these energy inefficiencies.

We investigate the effects of microwave pretreatment on various rock samples, as quantified by uniaxial compression tests. In particular, we assign sample heterogeneity based on a Gaussian statistical distribution and invoke a damage model for elemental tensile and compressive stresses based on the maximum tensile stress and the Mohr–Coulomb theories, respectively. We further couple the electromagnetic, thermal and solid displacement relations using finite element modeling.

(1) Increased power intensity during microwave pretreatment results in decreased axial compressive stress. (2) Leveraging statistics to induce variable compressive and tensile strength can greatly facilitate sample heterogeneity and prove necessary for damage modeling. (3) There exists a nonlinear trend to the reduction in smax with increasing power levels, implying an optimum energy output efficiency to create the maximum degradation-power cost relationship.

Previous research in this area has been largely limited to two-dimensional thermo-electric models. The onset of high-performance computing has allowed for the development of high-fidelity, three-dimensional models with coupled equations for electromagnetics, heat transfer and solid mechanics.

The determination of mode-I fracture toughness of brittle structural materials by means of the notched Brazilian disc configuration is studied. Advantage is taken of a recently introduced analytical solution and, also, of data provided by an experimental protocol with notched marble specimens under diametral compression using the loading device suggested by International Society for Rock Mechanics (ISRM) and also the three-dimensional digital image correlation (3D-DIC) technique.

The analytical solution highlighted the role of geometrical factors, like, for example, the width of the notch, which are usually disregarded. The data of the experimental protocol were comparatively considered with those concerning the response of the specific material under uniaxial tensile load.

This combined study provided interesting data concerning some open issues, as it is the exact crack initiation point and the level of the critical load causing crack initiation. It was definitely indicated that the crack initiation point is not a priori known (even for notched specimens) and, also, that the maximum recorded load does not correspond by default to the critical load responsible for the onset of catastrophic macroscopic fracture.

It was suggested that the load considered critical one for the determination of mode-I fracture toughness K_IC is erroneous. At a load equal to about 70% of the maximum one, a process zone is formed (zone of non-reversible phenomena) around the notch's crown, designating termination of the validity of any linear elastic solution used to determine the normalized stress intensity factors (SIFs). Moreover, at a load level equal to about 95% of the macroscopically observed fracture load, crack propagation has already begun. Therefore, the experimental procedure must be monitored with additional equipment, providing an overview of the displacement field developed during loading.

Despite the many advantages this type of equipment offers, there are still some major drawbacks. Linear cutting machine (LCM) cannot accurately simulate the true rock-cutting process as 1. it does not account for the circular path along which tunnel boring machine (TBM) disk cutters cut the tunnel face, 2. it does not accurately model the position of a disk cutter on the cutterhead, 3. it cannot perfectly replicate the rotational speed of a TBM. To enhance the knowledge of these issues and in order to mimic the real rock-cutting process, a new lab testing equipment was developed by Hyundai Engineering and Construction.

A new testing machine called rotary cutting machine (RCM) is designed to simulate the excavation process of hard-rock TBMs and includes features such as TBM cutterhead, RPM simulation, constant normal force mode and constant penetration rate mode. Two sets of tests were conducted on Hwandeung granite using different disk cutter sizes to analyze the cutting forces in various excavation modes. The results are analyzed using statistical analysis and dimensional analysis. A new model is generated using dimensional analysis, and its results are compared against the results of actual cases.

The effectiveness of the new RCM test was demonstrated in its ability to apply various modes of excavation. Initial analysis of chip size revealed that the thickness of the chips is largely dependent on the cutter spacing. Tests with varying RPM showed that an increase in RPM results in an increase in the normal force and rolling force. The cutting coefficient (CC) demonstrated a linear correlation with penetration. The optimal specific energy is achieved at an S/p ratio of around 15. However, a slightly lower S/p ratio can also be used in the design if the cutter specifications permit. A dimensional analysis was utilized to develop a new RCM model based on the results from approximately 1200 tests. The model's applicability was demonstrated through a comparison of TBM penetration data from 26 tunnel projects globally. Results indicated that the predicted penetration rates by the RCM test model were in good agreement with actual rates for the majority of cases. However, further investigation is necessary for softer rock types, which will be conducted in the future using concrete blocks.

The originality of the research lies in the development of Hyundai Engineering and Construction’s advanced full-scale laboratory rotary cutting machine (RCM), which accurately replicates the excavation process of hard-rock tunnel boring machines (TBMs). The study provides valuable insights into cutting forces, chip size, specific energy, RPM and excavation modes, enhancing understanding and decision-making in hard-rock excavation processes. The research also presents a new RCM model validated against TBM penetration data, demonstrating its practical applicability and predictive accuracy.

This study aims to propose an approach to assess the security of supply (SS) in a coal-fired electricity generation supply chain subject to public price regulation in Brazil. This study characterizes the Brazilian scenario of coal-fired electricity generation, which represents less than 3.5% of the energy sources.

Data from six mining companies that supply a coal plant were analyzed in a case study. The risks were characterized and objectively estimated through a synthetic multidimensional index. Structural changes in the earnings before interest, taxes, depreciation, amortization and exploration indicator time series of coal companies (CC) were statistically detected.

Empirical evidence demonstrates that the supply chain has a low disruption risk (SS index equal to 0.74). However, when suppliers are individually analyzed, 48.64% of all coal shows moderated disruption risk, and 2.51% is under high risk. In addition, this study finds a drop in the financial results of CC related to public regulation of coal prices. This impacts the security of coal supply.

This study discusses the influence of legal and regulatory policy risks in a coal power generation supply chain and the implications of the SS index as a management tool.

A novel SS index is presented and empirically operationalized, and its dimensions – environmental, occupational, operational, economic-financial and supply capacity – are analyzed.

This paper aims to employ a multivariate nonlinear regression analysis to establish a predictive model for the final fracture area, while accounting for the impact of individual parameters.

This analysis is based on the numerical simulation data obtained, using the hybrid finite element–discrete element (FE–DE) method. The forecasting model was compared with the numerical results and the accuracy of the model was evaluated by the root mean square (RMS) and the RMS error, the mean absolute error and the mean absolute percentage error.

The multivariate nonlinear regression model can accurately predict the nonlinear relationships between injection rate, leakoff coefficient, elastic modulus, permeability, Poisson’s ratio, pore pressure and final fracture area. The regression equations obtained from the Newton iteration of the least squares method are strong in terms of the fit to the six sensitive parameters, and the model follow essentially the same trend with the numerical simulation data, with no systematic divergence detected. Least absolutely deviation has a significantly weaker performance than the least squares method. The percentage contribution of sensitive parameters to the final fracture area is available from the simulation results and forecast model. Injection rate, leakoff coefficient, permeability, elastic modulus, pore pressure and Poisson’s ratio contribute 43.4%, −19.4%, 24.8%, −19.2%, −21.3% and 10.1% to the final fracture area, respectively, as they increased gradually. In summary, (1) the fluid injection rate has the greatest influence on the final fracture area. (2)The multivariate nonlinear regression equation was optimally obtained after 59 iterations of the least squares-based Newton method and 27 derivative evaluations, with a decidability coefficient R2 = 0.711 representing the model reliability and the regression equations fit the four parameters of leakoff coefficient, permeability, elastic modulus and pore pressure very satisfactorily. The models follow essentially the identical trend with the numerical simulation data and there is no systematic divergence. The least absolute deviation has a significantly weaker fit than the least squares method. (3)The nonlinear forecasting model of physical parameters of hydraulic fracturing established in this paper can be applied as a standard for optimizing the fracturing strategy and predicting the fracturing efficiency in situ field and numerical simulation. Its effectiveness can be trained and optimized by experimental and simulation data, and taking into account more basic data and establishing regression equations, containing more fracturing parameters will be the further research interests.

The nonlinear forecasting model of physical parameters of hydraulic fracturing established in this paper can be applied as a standard for optimizing the fracturing strategy and predicting the fracturing efficiency in situ field and numerical simulation. Its effectiveness can be trained and optimized by experimental and simulation data, and taking into account more basic data and establishing regression equations, containing more fracturing parameters will be the further research interests.

Thermal fractures initiated under cooling at the surfaces of a 2-D or 3-D structure propagate, arrest and coalesce, leading to its structural failure and material-property changes, while the same processes can happen in the rock mass between parallel hydraulic fractures filled with cold fluid, leading to enhanced fracture connectivity and permeability.

This study used a 2-D plane strain fracture model for mixed-mode thermal fractures from two parallel cooling surfaces. Fracture propagation was governed by the theory of linear elastic fracture mechanics, while the displacement and temperature fields were discretized using the adaptive finite element method. This model was validated using two numerical benchmarks with strong fracture curvature and then used to simulate the propagation and coalescence of thermal fractures in a long rock mass.

Modeling results show two regimes: (1) thermal fractures from a cooling surface propagate and arrest by following the theoretical solutions of half-plane fractures before the unfractured portion decreases to 20% rock-mass width and (2) some pairs of fractures from the opposite cooling surfaces tend to eventually coalesce. The fracture coalescence time is in a power law with rock-mass width.

These findings are relevant to both subsurface engineering and material engineering: structure failure is a key concern in the latter, while fracture coalescence can enhance the connectivity of thermal and hydraulic fractures and thus reservoir permeability in the former.

The goal of this study is to determine the values of the process parameters that should be used during the machining of ceramic tile using the abrasive water jet (AWJ) process in order to achieve the lowest possible values for surface roughness and kerf taper angle.

In the present work, ceramic tile is processed by the AWJ process and experimental data were recorded using the RSM approach based Box–Behnken design matrix. The input process factors were water jet pressure, jet traverse speed, abrasive flow rate and standoff distance, to determine the surface roughness and kerf taper angle. ANOVA was used to check the adequacy of model and significance of process parameters. Further, the elite opposition-based learning grasshopper optimization (EOBL-GOA) algorithm was implemented to identify the simultaneous optimization of multiple responses of surface roughness and kerf taper angle in AWJ.

The suggested EOBL-GOA algorithm is suitable for AWJ of ceramic tile, as evidenced by the error rate of ±2 percent between experimental and predicted solutions. The surfaces were evaluated with an SEM to assess the quality of the surface generated with the optimal settings. As compared with initial setting of the SEM image, it was noticed that the bottom cut surface was nearly smooth, with less cracks, striations and pits in the improved optimal results of the SEM image. The results of the analysis can be used to control machining parameters and increase the accuracy of AWJed components.

The findings of this study present an innovative method for assessing the characteristics of the nontraditional machining processes that are most suited for use in industrial and commercial applications.

Multi-well hydrofracturing is an important technology to create new fractures and expand existing fractures to increase reservoir permeability. The propagation morphology of the fracture network is affected by the disturbance between the fractures initiation sequences and spacings between adjacent wells. However, it remains unclear how well spacing and initiation sequences lead to fracture propagation, deflection and connection.

In this study, the thermal-hydro-mechanical coupling effect in the hydrofracturing process was considered, to establish a finite element-discrete element model of multistage hydrofracturing in a horizontal well. Using typical cases, the unstable propagation of hydraulic fractures in multiple horizontal wells was investigated under varying well spacing and initiation sequences. Combined with the shear stress shadow caused by in situ stress disturbed by fracture tip propagation, the quantitative indexes of fracture propagation such as length, volume, displacement vector, deflection and unstable propagation behavior of the hydrofracturing fracture network were analyzed.

The results show that the shear stress disturbance caused by multiple hydraulic fractures is a significant factor in multi-well hydrofracturing. Reducing the spacing between multiple wells increases the stress shadow area and aggravates the mutual disturbance and deflection between the fractures. The quantitative analysis results show that a decrease of well spacing reduces the total length of hydraulic fractures but increases the total volume of the fracture; compared with sequential and simultaneous fracturing, alternate fracturing can effectively reduce stress shadow area, alleviate fracture disturbance and generate larger fracture propagation length and volume.

The numerical models and results of the unstable propagation and stress evolution of the hydraulic fracture network under thermal-hydro-mechanical coupling obtained in this study can provide useful guidance for the evaluation and design of rock mass fracture networks in deep unconventional oil and gas reservoirs.

This paper aims to develop an efficient algorithm combining straightforward response surface functions with Monte Carlo simulation to conduct seismic reliability analysis in a systematical way.

The representative slip surfaces are identified and based on to calibrate multiple response surface functions with acceptable accuracy. The calibrated response surfaces are used to determine the yield acceleration in Newmark sliding displacement analysis. Then, the displacement-based limit state function is adopted to conduct seismic reliability analysis.

The calibrated response surface functions have fairly good accuracy in predicting the yield acceleration in Newmark sliding displacement analysis. The seismic reliability is influenced by such factors as PGA, spatial variability and threshold value. The proposed methodology serves as an effective tool for geotechnical practitioners.

The multiple sources of a seismic slope response can be effectively determined using the multiple response surface functions, which are easily implemented within geotechnical engineering.