Rock Mechanics
Ehsan Taheri; Reza Mohammadpour; Mohammad Hossein Mokhtarzadeh
Abstract
In recent years, the demand for new trenchless methods has dramatically risen. Pipe jacking is a trenchless method widely used in recent years. Ground deformation is one of the significant parameters that may lead to unrepairable harm to facilities and even people. So, ground deformation analysis is ...
Read More
In recent years, the demand for new trenchless methods has dramatically risen. Pipe jacking is a trenchless method widely used in recent years. Ground deformation is one of the significant parameters that may lead to unrepairable harm to facilities and even people. So, ground deformation analysis is necessary for safety and design reasons. The present study analyzes the factors affecting ground deformation during pipe jacking. This is a descriptive-interventional study. Pipe jacking causes soil displacement in three dimensions (3-D). Therefore, 3-D numerical methods were applied for analysis. In this study, numerical simulation was performed using PLAXIS finite element numerical software, taking the case study into account. The effect of each parameter on the ground deformation pattern was studied in three directions; the uplift and their exact position were then analyzed. It should be noted that displacement analyses were performed in two areas: pipe crown and ground surface. Also, the relation of each parameter was estimated with the ground subsidence. Finally, the effect of each different factor and their sensitivity index were determined using sensitivity analysis. The highest subsidence occurs at the end of the shield due to stress relaxation. Considering the results, it was found that the relationship between the internal friction angle and subsidence is linear and direct. The relationship between the elastic modulus and subsidence is also linear but indirect. The results indicate that the most sensitive factor of ground deformation is the diameter, but the least sensitive factor is the face pressure.
Rock Mechanics
Aref Jaberi; Shokroallah Zare
Abstract
Unlike the mechanical properties of intact rock, which can be obtained on a laboratory scale, estimating the mechanical properties of the jointed rock mass is very difficult due to the presence of different joints and the complexity of the joints. Therefore, to calculate the mechanical parameters of ...
Read More
Unlike the mechanical properties of intact rock, which can be obtained on a laboratory scale, estimating the mechanical properties of the jointed rock mass is very difficult due to the presence of different joints and the complexity of the joints. Therefore, to calculate the mechanical parameters of the jointed rock mass and use the continuous media theory of the jointed rock mass, it is necessary to calculate the Representative Element Volume (REV) of the rock mass. In this study, the Discrete Element Method (DEM) and the mechanical index of strength were used to investigate the effect of persistent and non-persistent joint angles, as well as model size on the REV in x, y, and z directions. The numerical results showed that by changing the joint angles and side length, both the strength and the REV of the rock mass were affected. The maximum representative side length for the persistent joint in the x and z directions occurred at angles of 60° and 75°, respectively. The minimum strength was obtained for joints in the x and z directions at a 45° angle. Finally, the REV for persistent and non-persistent joints is calculated as 10*0.5*8m and 4*0.5*4m, respectively.
Rock Mechanics
Naeem Abbas; Li Kegang
Abstract
The study examined the influence of cohesion, friction angle, and tunnel diameter on stability within engineering and geotechnical frameworks, while considering the consequences of nearby excavations on the overall stability assessment. The results show that a higher angle of internal friction leads ...
Read More
The study examined the influence of cohesion, friction angle, and tunnel diameter on stability within engineering and geotechnical frameworks, while considering the consequences of nearby excavations on the overall stability assessment. The results show that a higher angle of internal friction leads to a decrease in soil stability number and weighting coefficient. Tunnel diameter significantly affects face support pressure, with larger diameters requiring stronger support due to increased stress. Higher friction angles help stabilize tunnel faces and mitigate diameter-related pressure effects. Stress redistribution around the tunnel is significant within 2 meters from the center, transitioning to elastic behavior elsewhere. A safety factor of 1.3 ensures tensile failure prevention in single and twin tunnels. Balanced stress distribution between tunnels with a slight difference is observed under isotropic in-situ stress. Numerical modeling enhances stress estimations and reveals changes during tunnel excavation, weakening the rock mass. Ground reaction curve analysis with support measures shows reduced tunnel convergence after implementation, suggesting support strategies like extended bolts using updated rock mass rating. The study improves tunnel design and stability assessment by comprehensively understanding stress redistribution and support strategies.
Exploitation
Behnam Alipenhani; Mehran Jalilian; Abbas Majdi; Hassan Bakhshandeh Amnieh; Mohammad Hossein Khosravi
Abstract
The paper presents the effect of the dip of joints, joint spacing, and the undercutting method on the height of the caving in block caving. The obtained results show that among the three investigated parameters, respectively, the dip of joints, undercutting method, and joint spacing have the greatest ...
Read More
The paper presents the effect of the dip of joints, joint spacing, and the undercutting method on the height of the caving in block caving. The obtained results show that among the three investigated parameters, respectively, the dip of joints, undercutting method, and joint spacing have the greatest effect on increasing the height of the caving zone. Comparing the data obtained from physical and numerical modeling shows a 97% match. Also, by increasing the joint spacing from 4 to 6 cm, 14%, from 6 to 8 cm, about 35%, and from 8 to 10, about 50%, the height of the caving zone has decreased. Regarding the dip of the joint, with the dip increasing from 30 to 45 degrees, about 3% of the caving height decreases. By increasing the dip of the joint from 45 to 60 degrees, the caving height has decreased by 42%. By increasing this value from 60 to 75 degrees, the caving height has increased by 50%. Also, changing the undercutting method from symmetric to advanced undercutting has increased the caving height by 40%. Additionally, three mathematical models have been proposed based on the shape of the caving zone in physical modeling.
Vaibhav Sharma; Andy Kwame Yeboah; Joshua Asare; Natillio Pillay; Jaspreet Singh
Abstract
The presence of any underground cavity in the soil stratum can seriously harm the structural performance of the overlying facility. These may develop because of mining, tunneling, water, and gas networks or outdated channels. In the present investigation, a circular void is considered, and its effect ...
Read More
The presence of any underground cavity in the soil stratum can seriously harm the structural performance of the overlying facility. These may develop because of mining, tunneling, water, and gas networks or outdated channels. In the present investigation, a circular void is considered, and its effect on the surface strip footing (in the form of ultimate load (UL), ultimate settlement (US), footing tilting, and footing horizontal displacement (HD)) is studied using numerical simulation. The variable parameters are load eccentricity (e), load inclination (α), and geogrid reinforcement location (u). It is observed that as the load inclination and eccentricity increases, the UL decreases. For instance, in the unreinforced soil, u/B = 0, at load inclination of α = 0°, 10°, 20°, and 30°, the UL is 249, 200, 142, and 97 kN/m, respectively. Moreover, as the geo-grid location is changed, the UL first increases when placed near the footing (u/B = 0.10), and thereafter, starts to decrease as the distance between footing and geo-grid increases. For instance, the UL is 249, 278, 267, 260, 259, and 256 kN/m when e/B = 0.0, α = 0°, and u/B varies from 0 to 0.5 with an increment of 0.1. The tilting increases as the eccentricity is increased; for example, u/B = 0.0 for α = 0°; the tilting values are 0°, 0.12°, 0.31°, and 0.61°. Moreover, as the load eccentricity increases, the HD decreases (for u/B = 0.1 and α = 10°, the HD is 4.20, 3.5, 3.00, and 2.60 mm, respectively.
Alankrit Walia; Amrit Kumar Roy
Abstract
The complex geography of the Himalayan mountain range, along with the natural circumstances that already exist and the ways in which people have influenced and intervened in the region- makes various regions of the range vulnerable to slope instability. The slope stability of the area that is the subject ...
Read More
The complex geography of the Himalayan mountain range, along with the natural circumstances that already exist and the ways in which people have influenced and intervened in the region- makes various regions of the range vulnerable to slope instability. The slope stability of the area that is the subject of this work is evaluated in Palampur, which is in the Kangra district of Himachal Pradesh. The primary objective of this work is to ensure that the slope remains stable so that the nearby three-sided residential structures and the highway remain protected. After the site visit, the geo-technical studies, which include testing in the form of bore holes and in the laboratory, are carried out. After evaluating the geo-technical technical report, the next step in the process is to begin the analysis of the slope's stability. In order to do an analytical analysis of the slope stability, the area has been subdivided into three portions, and labelled A-A, B-B, and C-C, respectively. Using the numerical modelling approaches, the mitigation design parameters for the area and the circular slip failure are computed. These calculations are based on the geo-technical characteristics of the studied area that have been specified. The factor of safety is calculated for both the natural and stable scenarios by the program. Because of this, some preventative steps and a few improvements are suggested.
B. Alipenhani; A. Majdi; H. Bakhshandeh Amnieh
Abstract
The present work aims at implementing Response Surface Methodology (RSM) in order to generate a statistical model for Minimum Required Caving Span (MRCS) and estimate both the individual and mutual effects of the rock mass parameters on rock mass cavability. The adequate required data is obtained from ...
Read More
The present work aims at implementing Response Surface Methodology (RSM) in order to generate a statistical model for Minimum Required Caving Span (MRCS) and estimate both the individual and mutual effects of the rock mass parameters on rock mass cavability. The adequate required data is obtained from the result of numerical modeling. In this work, various arrays of numerical simulations (480 models) are carried out using the UDEC software in order to study the rock mass cavability thoroughly. The effect of each individual parameter and their mutual effect on MRCS are investigated by means of ANOVA. ANOVA indicates that all the chosen parameters (depth, dip of the joint, number of joints, angle of friction of the joint surface, and joint spacing) highly affect MRCS. In other words, the results of ANOVA are in high agreement with the results of the conventional sensitivity analysis. Moreover, a combination of joint spacing and joint inclination has the highest mutual effect on MRCS, and a combination of undercut depth and joint spacing has the lowest effect on MRCS.
F. Abbasi; M.H. Khosravi; A. Jafari; Ali Bashari; B. Alipenhani
Abstract
The instability of the roof and walls of the tunnels excavated in coal mines has always attracted the attention of the miners and experts in this field. In this work, the instability of the main tunnels of the mechanized Parvadeh coal mine in Tabas, Iran, at the intersection with coal seam is studied. ...
Read More
The instability of the roof and walls of the tunnels excavated in coal mines has always attracted the attention of the miners and experts in this field. In this work, the instability of the main tunnels of the mechanized Parvadeh coal mine in Tabas, Iran, at the intersection with coal seam is studied. The main tunnels of this mine show significant horizontal displacements due to the complex ground conditions and great depth. The behavior of the rock mass surrounding the tunnel is investigated using various experimental methods, and according to the results obtained, the surrounding rock mass has squeezing conditions. In order to analyze the stability of the main tunnels, a series of 2D and 3D numerical modelings are performed using the FLAC2&3D finite difference software, and the results obtained are compared with the actual displacement values recorded in the walls of the main tunnels of the mine. The analysis results show that the tunnels under study are unstable with a steel frame support, and therefore, the use of different support systems for the stabilization is investigated. The results of modeling different types of support systems show that the use of shotcrete instead of galvanized sheet (as strut) does not have a significant effect on the reduced displacements. Also although the installation of steel sets is very effective in preventing the displacement of the walls, due to the swelling problems in the tunnel bottom and the placement of the conveyor and haulage rail, it cannot be used in practice. Finally, the use of truss bolt has yielded good results, and it can be proposed as a new support system in these tunnels. In addition, the modeling results show that in case the coal seam is higher than the tunnel foot, less displacement will occur in the tunnel walls compared to the other cases. In other words, changing the tunnel level in the future excavations can help reduce the displacements.
Seyed M. Pourhashemi; K. Ahangari; J. Hassanpour; Seyed M. Eftekhari
Abstract
Mechanized tunneling in rocks is based on fracture propagation and rock fragmentation under disc cutters. Rock chipping is an efficient kind of fragmentation process, while the grinding process may occur under special conditions. The cutter-head penetration is an appropriate parameter involved in order ...
Read More
Mechanized tunneling in rocks is based on fracture propagation and rock fragmentation under disc cutters. Rock chipping is an efficient kind of fragmentation process, while the grinding process may occur under special conditions. The cutter-head penetration is an appropriate parameter involved in order to distinguish between the chipping and grinding processes in rock cutting. In this work, the grinding and chipping processes are investigated in the Uma-Oya water conveyance tunnel in Sri Lanka. The Uma-Oya project is a water transfer, hydropower, and irrigation system in the SE part of the central highland region of Sri-Lanka. From a geological viewpoint, most parts of the tunnel route in the studied section consist of very strong and abrasive metamorphic rocks that potentially are susceptible to grinding occurrence during the boring process under disc cutters. In this work, firstly, data processing is performed in order to identify the boundary between chipping and grinding. Then the chipping and grinding processes are modeled using the practical numerical and artificial intelligent methods. In the numerical modeling stage, we try to make the modeling as realistic as possible. The results obtained from these modeling methods show that for the penetrations less than 3 mm/rev, the grinding process is dominant, and for the penetrations more than 3 mm/rev, rock chipping occurs. Also, in the numerical modeling, no significant fracture expansion is observed in the rock when the penetration is less than 3 mm/rev. Moreover, it can be seen in the numerical modeling of the chipping process that the propagated fractures come together and the chips are created.
M. Zahiri; K. Goshtasbi; J. Khademi Hamidi; K. Ahangari
Abstract
There is a direct relationship between the efficiency of mechanized excavation in hard rocks and that of disc cutters. Disc cutter wear is an important effective factor involved in the functionality of tunnel boring machines. Replacement of disc cutters is a time-consuming and costly activity that can ...
Read More
There is a direct relationship between the efficiency of mechanized excavation in hard rocks and that of disc cutters. Disc cutter wear is an important effective factor involved in the functionality of tunnel boring machines. Replacement of disc cutters is a time-consuming and costly activity that can significantly reduce the TBM utilization and advance rate, and has a major effect on the total time and cost of the tunneling projects. When these machines bore through hard rocks, the cutter wear considerably affects the excavation process. To evaluate the performance of the cutters, first, it is essential to figure out how they operate the rock cutting mechanism; secondly, it is important to identify the key factors that cause the wear. In this work, we attempt to introduce a comprehensive numerical method for estimation of disc cutter wear. The field data including the actual cutter wear more than 1000 pieces and the geological parameters along the Kani-Sib transmission tunnel in the northwest of Iran are compiled in a special database that is subjected to a statistical analysis in order to reveal the genuine wear rule. The results obtained from the numerical method indicate that with an increase in the wear of disk cutter up to 25 mm, the applied normal and rolling forces can be multiplied by 2.9 and 2.7, respectively, and by passing the critical wear, the disk cutters lose their optimal performance. This method also shows that confining pressure will increase the wear of the disc cutter. By the proposed formulation, the cutter consumption rate can be predicted with a high accuracy.
Rock Mechanics
H. Sarfaraz; M. Amini
Abstract
One of the most important instabilities of rock slopes is toppling failure. Among the types of toppling failure, block-flexural failures are more common instability which occurs in nature. In this failure, some rock blocks break because of tensile stresses, and some overturn under their weights, and ...
Read More
One of the most important instabilities of rock slopes is toppling failure. Among the types of toppling failure, block-flexural failures are more common instability which occurs in nature. In this failure, some rock blocks break because of tensile stresses, and some overturn under their weights, and next to all of them topple together. Physical and theoretical modeling of this failure is studied by Amini et al. in 2015. Due to the complexity of this failure mechanism, no appropriate numerical model has been proposed so far. In this research, first, a literature review of toppling failure is summarized. Then, using UDEC software as distinct element method (DEM), the experimental models were analyzed numerically, and Voronoi joint model was applied to simulate the failure. The results of numerical simulations are compared with the outcomes of physical models and analytical solution. The comparison illustrates that numerical modeling has good agreement with corresponding experimental tests and theoretical approach. Also, the results show that although the mechanism of block-flexural toppling failure is complicated, the numerical code is well capable to analyze of this failure.
Rock Mechanics
H. Sarfaraz; M. H. Khosravi; M. Amini
Abstract
In layered and blocky rock slopes, toppling failure is a common mode of instability that may occur in mining engineering. If this type of slope failure occurs as a consequence of another type of failure, it is referred to as the secondary toppling failure. “Slide-head-toppling” is a type ...
Read More
In layered and blocky rock slopes, toppling failure is a common mode of instability that may occur in mining engineering. If this type of slope failure occurs as a consequence of another type of failure, it is referred to as the secondary toppling failure. “Slide-head-toppling” is a type of secondary toppling failures, where the upper part of the slope is toppled as a consequence of a semi-circular sliding failure at the toe of the slope. In this research work, the slide-head-toppling failure is examined through a series of numerical modeling. Phase 2, as a software written based on the finite element method, is used in this work. Different types of slide-head-toppling failures including blocky, block-flexural, and flexural are simulated. A good agreement can be observed when the results of the numerical modeling are compared with those for the pre-existing physical modeling and analytical method.
Rock Mechanics
E. Khorasani; M. Amini; M. F. Hossaini
Abstract
Bimrocks are complex geomaterials that are defined as mixtures of rocks composed of geotechnically significant blocks within a matrix of finer texture. Bimslopes are made from bimrocks and are usually seen in weathered and shallow environments. Some characteristics of blocks affecting the strength of ...
Read More
Bimrocks are complex geomaterials that are defined as mixtures of rocks composed of geotechnically significant blocks within a matrix of finer texture. Bimslopes are made from bimrocks and are usually seen in weathered and shallow environments. Some characteristics of blocks affecting the strength of bimrocks include VBP (Volumetric Block Proportion), orientation, and arrangement, which have important roles in the stability of bimslopes. Previous studies show that bimrocks usually have a specific block size distribution, and for a bimslope with height of “H”, the size of blocks is changed from 0.05H to 0.75H. In this paper, the influence of large blocks position on bimslope stability was investigated by the physical and numerical models. The blocks that had a dimension larger than 0.5H were considered as “large blocks”. In this work, first, thirty physical models were created and tested using a titling table machine. These models have a specific block size distribution and VBP with ellipsoidal blocks. The main variable of the models is large blocks position, where three categories including lower part of bimslope, upper part of bimslope, and sporadic state are considered. Based on the results of physical trials, thirty numerical models at the laboratory scale were generated using the finite element method. After comparing the physical and numerical models, which showed a good accordance, the numerical models were developed to the natural scale. The theoretical bimslopes investigated in this work showed that the position of large blocks had a significant influence on the stability of bimslopes.
Rock Mechanics
A. Akrami; M. Hosseini; H. Sodeifi
Abstract
Hydraulic fracturing is used in the oil industry in order to increase the index of production and processing in the wells whose efficiencies have been dropped due to a long-term harvest or the rocks around the wells are of low permeability. Since the hydraulic fracturing operation is costly, it is of ...
Read More
Hydraulic fracturing is used in the oil industry in order to increase the index of production and processing in the wells whose efficiencies have been dropped due to a long-term harvest or the rocks around the wells are of low permeability. Since the hydraulic fracturing operation is costly, it is of special importance to the project managers to determine the pressure required for hydraulic fracturing and the suitable pump for this operation. The numerical modelings used in this work are aimed to investigate the fracture pressure in the carbonate rocks of Bangestan reservoir in Ahvaz, Iran, and to determine a relationship between the pressure required for fracturing and the confining pressure. In this work, unlike the other ones in this field, the developed numerical models had no initial crack or fracture, and the path of the crack and how the crack grows were studied without any pre-determination and presumption. The results obtained show that, in most cases, the crack starts from the central part of the sample, and is extended to its two ends. The crack extension direction was along the borehole axis inside the sample and perpendicular to the lateral stress. The numerical modeling results were well-consistent with the experimental ones, indicating that the pump capacity constraints in the laboratory could be overcome through numerical modelings.
M. Nikkhah
Abstract
Nowadays acoustic emission (AE) testing based on the Kaiser Effect (KE) is increasingly used to estimate the in-situ stress in laboratories. In this work, this effect is assessed on cylindrical specimens in numerical simulations of the cyclic loadings including loading, unloading, and re-loading cycles ...
Read More
Nowadays acoustic emission (AE) testing based on the Kaiser Effect (KE) is increasingly used to estimate the in-situ stress in laboratories. In this work, this effect is assessed on cylindrical specimens in numerical simulations of the cyclic loadings including loading, unloading, and re-loading cycles using a 3D code called the particle flow code (PFC) based upon the distinct element method. To achieve this objective, at first, the numerical model is calibrated using a laboratory test performed on the selected sandstone specimens. The results obtained show that PFC and the distinct element code are useful tools used to investigate the damage and KE of a brittle rock. Also the results obtained by the triaxial modeling show that a combination of triaxial loading stresses change the results of uniaxial loading. Further, KE is influenced under confining stresses so that larger confining stresses lead to greater differences between the KE stress during the uniaxial and pre-stress loadings.
