Search Documents

By Eugie Kabwe

Rock bolts are components of tunnel support structures and their efficacy in reinforcement requires the appropriate estimation of the bolt pre-tension. However, cases of ineffective rock bolt application due to loosening are common which leads to the potential failure of the support structures, thus negatively affecting the tunnel stability. The effective employment of pre-tensioned bolt within support structures requires an appropriate installation tension as a result of sufficient bolt head/nut pre-load. Henceforth, this paper presents a load indicator washer (LIW) capable of withstanding assigned static loading as well as provides an appropriate pre-tension. The different LIW breadths are numerically analysed in the finite element code (ABAQUS) to estimate deformation when a designated 60-kN pre-load is applied. The numerical simulation employs a bi-linear isotropic hardening constitutive law governed by the von Mises yield criterion. The simulation also considers the nominal temperature of the nutLIW-rock bolt connection. Therefore, a long-term relaxation effect on bolt pre-tension as a result of stress relaxation is neglected. The compressive force required to achieve the LIW-predetermined deformation is the bolt pre-tension. It is drawn from the simulation results that the 1.5-, 2.0-, 3.0-, 3.5-, and 4.0-mm-thick LIWs are completely flattened when subjected to 28-, 42.5-, 65-, 76-, and 88-kN pre-load respectively.

By Eugie Kabwe

Rock bolts are components of tunnel support structures and their efficacy in reinforcement requires the appropriate estimation of the bolt pre-tension. However, cases of ineffective rock bolt application due to loosening are common which leads to the potential failure of the support structures, thus negatively affecting the tunnel stability. The effective employment of pre-tensioned bolt within support structures requires an appropriate installation tension as a result of sufficient bolt head/nut pre-load. Henceforth, this paper presents a load indicator washer (LIW) capable of withstanding assigned static loading as well as provides an appropriate pre-tension. The different LIW breadths are numerically analysed in the finite element code (ABAQUS) to estimate deformation when a designated 60-kN pre-load is applied. The numerical simulation employs a bi-linear isotropic hardening constitutive law governed by the von Mises yield criterion. The simulation also considers the nominal temperature of the nutLIW-rock bolt connection. Therefore, a long-term relaxation effect on bolt pre-tension as a result of stress relaxation is neglected. The compressive force required to achieve the LIW-predetermined deformation is the bolt pre-tension. It is drawn from the simulation results that the 1.5-, 2.0-, 3.0-, 3.5-, and 4.0-mm-thick LIWs are completely flattened when subjected to 28-, 42.5-, 65-, 76-, and 88-kN pre-load respectively.

By J. Stacho, Frankovská. J.

"Abstract The use of continuous flight auger (CFA) piles in Slovakia has increased significantly in recent decades. High-rise building of administrative centre Lakeside Park in Bratislava is based on piled-raft foundation. The technology of CFA has been used due to the advantages of this technology in comparison with driven and bored piles in local geotechnical conditions. Piles have been designed in accordance with Eurocode 7. Resistance of piles was calculated by using an analytical calculation model for bored piles according to Slovak technical standard and CFA technology has been taken into account by coefficient equal to 1.1. Settlement of piles was calculated by using a calculation model which is being used in Czech Republic and Slovakia. Two well instrumented static load tests were carried out to verify the pile design. The results of static load tests showed that the design was conservative. The paper presents analysis of numerical modelling of CFA pile and static load tests. Load-settlement curves and load distribution curves over pile shafts are compared. Software Plaxis 2011 has been used for numerical modelling. Hardening soil model allows considering over consolidation of Tertiary Silts and water pressure in deep sandy lenses. The article also presents the parametric study where the influences of pile-soil interface reduction factor, over consolidation ratio and material models of soils are analysed. The results of numerical modelling showed that the most appropriate results have been obtained when the interface reduction factor equal 0.9 and real behaviour of groundwater were considered. The conclusion of the paper confirms that mechanical behaviour of the ground can be better understood and determined by using FEM than analytical methods.1. IntroductionPiles are among the most frequently used elements of deep foundations. Piles transfer loads down to a suitable bearing stratum when the soil immediately below a construction is unsuitable. Piles are also used to reduce settlement of structure. At present, the requirements of limit settlement of structures are more stringent. Examples of these structures are commercial and administrative centres that form a complex of different tall buildings. High-rise buildings cause greater contact stresses than connected low-rise buildings. Pile foundations in such cases should be placed under the high-rise buildings to achieve an appropriate settlement.Resistance of piles can be determined from static load test results, in-situ tests and empirical experiences or calculations. Static load tests give the most reliable results, but these tests are also the most expensive. This is the reason for constant development of computational methods. Especially numerical methods give more accurate results by using more sophisticated material models e. g. hardening soil or hypoplastic model."

Jan 1, 2014

By P. C. Sui

Numerical investigation and parametric simulations of the coupled electrical conduction and mass diffusion in the GDL of a PEMFC are reported. The coupled transport is solved using a 2D, "under-the-rib" computational domain for the cathode GDL. The current density on the GDL/catalyst layer interface, which constitutes one of the boundary conditions for the GDL domain and reflects the activation overpotential in the catalyst layer and the ohmic loss in the membrane, is solved iteratively using a novel numerical algorithm. A parametric study is performed to investigate the effects on current density distribution of various operating conditions such as oxygen concentration and membrane resistance, and of design factors such as GDL geometry, anisotropic transport properties, and deformation under the land area. The results show that the current density distribution under the land area can be dominated by either electron transport or mass transport, depending on the operating regime.

Jan 1, 2005

By Baokuan Li, Fang Wang, Meilong Shan

"As the advantage of uniform bus loading, the electroslag remelting (ESR) furnaces with three - phase installation by delta or star connection are mainly used to produce the large or super large ingot. The electrical efficiency is significantly improved as the short circuiting through the molten slag. However, little is found about fiindamental research on the electromagnetic field of the ESR system with three -phase installation. A three-dimensional (3D) finite element model was developed to simulate the current density, magnetic field, electromagnetic force and Joule heating for a system of electrode, slag and ingot in the ESR system with three -phase installation.IntroductionThe metallurgical performance of the ingot is the fundament and precondition to obtain the large forging of high quality, which has a critical role in the producing of forging. The electroslag process (ESR) to product the metal materials of high homogenization and clarity become more and more popular[1].The electroslag process is a special smelting technology of producing small ingots, which originate in 1950[2].At present the three, four or six electrode is used in large scale electroslag furnace. It has been found that the deeper molten pool and the wider of the two-phase zone with only one large diameter of electrode to product large ingot effect the solidification microstructure of ingot[3]. However, the molten pool becomes shallow relatively with several electrodes producing the large ingot, which is benefit to the solidification microstructure of ingot. Three electrodes electroslag remelting furnace take the three phase equilibrium power supply to making the load and electric power network balance, avoiding the distortion phenomenon induced by single phase power supply and saving electric power. It also can uniformize the distribution of heat source."

Jan 1, 2012

The N3500 orebody is the northern continuation of the 3500 orebody at Xstrata Copper's Enterprise Mine and occurs at a depth of 1.3 to 1.8 km. At the nearby S3500 orebody, extreme deformation was observed in access drives where there was a high number of bedding planes per metre combined with either graphite or talc infill. The observed failure mechanism was largely controlled by the geometry of the inclined bedding submitted to high horizontal stresses. Due to the anisotropic strength and deformation behaviour of the shale rock mass, traditional isotropic methods of geomechanical property estimation based upon GSI are not appropriate. A numerical modelling approach, based upon the distinct element codes UDEC (2D) and 3DEC (3D), has been developed to accurately simulate the failure mechanisms and deformation observed in the initial infrastructure surrounding the N3500 orebody.

Nov 25, 2010

By K. Sasaki, T. Murakami, S. Akibayashi, Q. Wei

A numerical analysis and an experiment were done on groundwater flow past a bore hole with a cylindrical obstacle at its center. Flow around the obstacle in the bore hole was studied in detail. Two dimensional governing equations to simulate this flow were written for the hydraulic head and seepage velocity in a porous medium and for the vorticity, stream function, and velocity in the bore hole. These equations were non-dimensionalized and solved by using the alternating direction implicit and successive-over-relaxation finite difference methods. Flow visualization experiments and velocity measurements in the bore hole were done with a two dimensional model composed of a porous medium. The numerical analysis procedures were validated by a comparison of numerical and experimental results. Effects of the Reynolds number and obstacle size on flow in the bore hole are clearly shown in the velocity and streamline profiles. The quantitative relation between the seepage velocity in the porous medium and the velocity of flow around the obstacle in the bore hole also were estimated.

Jan 1, 1996

By H. Rafiei Renani, C. D. Martin

"The Cohesion Weakening Friction Strengthening (CWFS) is adopted in numerical analysis to capture the observed progressive failure of hard rock pillars. As a case study, pillars at the Elliot Lake mines with width to height ratios of 0.5 to 3 were analyzed in FLAC3D. Pillars with width to height ratios of 0.5 and 1 show a very similar strain softening behavior while with width to height ratios of 2 and especially 3 indicated prolonged strain hardening. The CWFS model successfully captured the mechanism of progressive failure observed in hard rock pillars including initiation of spalling at the corners, progress of spalling towards the centre and formation of hour-glass shaped pillar, and eventually shearing through the width of pillar. Both the CWFS model and site specific empirical formula successfully separated the stable and failed cases in the Elliot Lake mines. While both models are in reasonable agreement for width to height ratios of less than 2, the empirical formula gives increasingly lower strength values for wider pillars which may result in excessively conservative and uneconomical design.INTRODUCTIONPillars are the portion of in situ rock left between the adjacent openings to maintain the stability of excavations. Conditions of equilibrium require that the loads previously carried by the excavated rocks be transferred to the pillars resulting in elevated stresses within the pillar. Hence, prediction of the behavior of pillars under these loads is important for underground design. Determination of the strength of pillars is a classic subject in underground mining. In conventional mining methods such as room and pillar operations, economic incentives point towards higher excavation ratios i.e., leaving smaller pillars between the mined excavations, while safety and stability requirements favor wider stronger pillars. In some mining methods such as caving operations, the problem is even more challenging. On one hand, the whole operation is directly dependent on the proper continuous flow of caved ore between the drawpoints which requires optimizing the spacing between drawpoints and thereby often leaving narrower pillars. On the other hand, the pillars should be wide enough to ensure the safety of mining operation. Both the blockage of drawpoints as a result of wide pillars and uncontrolled failure of narrow pillars can jeopardize the investments and mining operation.Therefore, design of pillar poses a somewhat unique challenge to rock engineers where simple conservative approaches frequently used in other design methods are no longer acceptable. The behavior of pillar under elevated loads should be predicted accurately."

Jan 1, 2015

By Hwa-Soo Park, J. K. Yoon, Hoseok Nam

In a thin slab casting process, control of fluid flow in mold is particularly difficult due to the high casting speed and large aspect ratio of the mold and the understanding of the fluid flow motion as well as heat transfer in mold is regarded as the most important information because the melt delivery is directly related to the problems of shell growth nonuniformity, surface turbulence and mold powder entrapment. A mathematical model has been developed for the coupled analysis of fluid flow, heat transfer and solidification in funnel type mold using finite volume method based on a body fitted coordinate and the characteristics of the transport phenomena in the mold of thin slab caster were analyzed by numerical simulation. As a result of simulation, the basic flow pattern is characterized as four large recirculations and two small eddies near the narrow faces of mold, which is somewhat different from that of parallel type mold or other conventional slab caster. The heat transfer model accurately predicted the mold positions susceptible to cracking as a result of thermal stress and the predicted shell thickness shows good agreements with the measured shell thickness on breakout products. We also investigated the influences of some operation parameters such as casting speed, nozzle type, submerging depth and mold material on fluid flow, heat transfer and distribution of solidifying shell and various information for the determination of optimal operating conditions were obtained.

Jan 1, 2000

By Khamis Y. Haramy

Quality assurance and quality control of drilled shafts has become a concern due to difficulties in accurately locating construction anomalies and determining load bearing capacity of defected drilled shafts. Various non-destructive evaluation techniques have been developed to estimate the integrity of the concrete in drilled shafts. While these techniques have been widely accepted, variables and unknowns can affect the measurement results. Results are typically difficult to interpret, leading to unnecessary construction delays and possible litigation over shaft integrity. In addition, influences of surrounding ground materials, stress states under different load conditions, and crack and residual stress development during concrete curing further complicate determination of shaft performance. This paper focuses on the evaluation of load bearing capacity of defected drilled shafts under various conditions by numerical analysis and modeling.

Jan 1, 2007

By V. A. Shea-Albin

Sandstone escarpments surrounding western coal-bearing plateaus serve as critical habitat for animal species. Effects of mining subsidence on escarpment stability are not well documented; therefore, regulatory agencies may limit extraction of coal reserves underlying escarpments. Because millions of tons of this coal exist on Federal Lands, associated public revenues from Federal royalties, and valuable coal reserves, risk being lost. The U. S. Bureau of Mines is using numerical modeling techniques to analyze effects of subsidence on escarpment stability. Finite element models that incorporate bedding plane slip, gob formation, and movement along joints simulate escarpment subsidence observed in the field.

Jan 1, 1994

By Shinji Muraishi, Shinji Kumai, Yasutaka Kedo, Takashi Kambe

"Magnetic pulse forming (MPF) is one of the extreme high-speed forming processes. A metal sheet which is accelerated by electromagnetic force is collided onto a mold and fine and sharp details of the product can be obtained. The deformation mechanism of the MPFed metal cannot be understood only by the microstructure observation of the deformed metal. The deformation behavior of the metal should be investigated also by considering the local microstructure change with strain, strain rate and temperature change during deformation. Numerical analysis is an effective method to achieve this requirement. In this study, pure Al and Cu sheets were deformed by MPF. Impact velocity of the metal sheet was calculated by using ANSYS Emag-Mechanical. Based on the obtained impact velocity, the deformation behavior of the metal sheet was reproduced, and local strain, strain rate, and temperature were analyzed by using ANSYS AUTODYN. The relationship between the local microstructure change of Al and Cu sheets and the deformation parameters was examined.INTRODUCTIONMagnetic pulse forming (MPF) is one of the high-speed forming processes using electromagnetic force. Fine and sharp details of the product can be obtained, and springback after deformation is extremely reduced in comparison to the conventional forming processes(Psyk, Risch, Kinsey, Tekkaya, & Kleiner, 2011). Figure 1 shows a schematic diagram of MPF. A metal sheet is placed over the coil in the discharge circuit. The discharge current runs through the coil, and magnetic fieldis producedaround the coil. An eddy current is inducedin the metal sheet. Interaction between magnetic field and eddy current creates electromagnetic force. The metal sheet accelerated by electromagnetic forcecollides ontothe mold. Forming is completed in a few tens of micro-seconds.MPF is suitable for metals with a high electrical conductivity like Al and Cu.The MPFed metal may have characteristic microstructure because of high strain rate deformation."

Jan 1, 2018

By Juhyung Lee

Numerical analyses were performed to analyze the behavior of two different type of suction piles (mono type and group type) under various soil layer conditions (uniform silty clay layer, uniform silty sand layer, and two-layer deposit consisting of clay and sand layers) and lateral loading locations (top, middle, and bottom of the suction pile). According to the results, mono suction pile shows maximum resistance when load is applied on its middle, regardless of the soil layer conditions, while group suction piles varies significantly according to the allowance on head rotation. When head rotation is allowed, the behavior of group type piles become similar to the behavior of mono suction pile. Whereas, when rotation is prevented, resistance is maximized when load is applied on the top head. Moreover, group suction piles show higher efficiency regardless of the ground condition and geological stratum. For further advances, concrete block on the top of group suction piles are suggested as an effective alternative method to prevent partial rotation and maximize ground resistance for group suction piles.

By Ilhan Chang, Ernest Nsabimana, Kiseok Kwak, Juhyung Lee, Young-Hoon Jung

"Numerical analyses were performed to analyze the behavior of two different type of suction piles (mono type and group type) under various soil layer conditions (uniform silty clay layer, uniform silty sand layer, and two-layer deposit consisting of clay and sand layers) and lateral loading locations (top, middle, and bottom of the suction pile). According to the results, mono suction pile shows maximum resistance when load is applied on its middle, regardless of the soil layer conditions, while group suction piles varies significantly according to the allowance on head rotation. When head rotation is allowed, the behavior of group type piles become similar to the behavior of mono suction pile. Whereas, when rotation is prevented, resistance is maximized when load is applied on the top head. Moreover, group suction piles show higher efficiency regardless of the ground condition and geological stratum. For further advances, concrete block on the top of group suction piles are suggested as an effective alternative method to prevent partial rotation and maximize ground resistance for group suction piles.IntroductionWith the development of offshore technology, applications and real site installations of suction piles are increasing significantly nowadays. They are generally used as mooring anchors in various floating structures, marine platforms, and offshore wind turbines. According to its convenience in installation and removal, suction type foundations are becoming good alternatives for deep ocean support systems. The bearing capacity is similar to the bearing behavior of ordinary offshore structures (piles or drilled shafts). The most common use of suction caissons is anchors for catenary mooring lines, where the chain angle at the mudline is close to zero (Randolph, 2002; Sharma, 2004).Recently, it is reported that a huge isolated single suction pile can be replaced by a group of moderatesized suction piles (Randolph et al. 2004). A single monopile is easy to be produced and installed in site. It is expected that monopile foundation is not applicable for deep depths because it has a limitation on penetration according to its huge side skin area (i.e. friction). Moreover, stiffness requirements sometimes render thicker pile wall conditions which are impractical to be produced."

Jan 1, 2017

By Alexander B. Reeb, Richard J. Varuso, George M. Filz, Jehu B. Johnson

"T-walls, which are concrete floodwalls in an inverted “T” shape, are an integral component of the Greater New Orleans Hurricane and Storm Damage Risk Reduction System (HSDRRS). In order to resist large lateral flood loads, T-walls are typically supported on long batter piles that extend through the weak and compressible clays and peats commonly encountered in southern Louisiana. Since the piles are battered, settlement of the compressible foundation soils can result in loads acting normal to the piles and induce pile bending moments. Therefore, the design of pile-supported T-walls is a complex soil-structure interaction problem that requires detailed and careful numerical analysis to accurately estimate the stresses and bending moments in the piles. For this paper, numerical analyses were performed for three cross sections that each included a two-pile T wall system that was installed within a 15-ft tall new levee embankment. The cross sections were identical except that the elevation of the T-wall base was varied. This paper discusses the analyses, and it compares the influence of settlements and flood loading on the pile stresses and bending moments. In addition, this paper includes general recommendations for the analysis and design of pile-supported T-walls in practice.INTRODUCTION AND BACKGROUNDThe recently completed Greater New Orleans Hurricane and Storm Damage Risk Reduction System (HSDRRS) includes 350 miles of levees and floodwalls and is designed to defend against a 100-year storm surge event. The floodwalls predominantly consist of I-walls and T-walls. I-walls are comprised of a sheet pile wall curtain capped with a monolithic concrete wall above grade. T-walls are similar to I-walls, but they include a wide concrete base at grade, which forms an inverted “T” shape, and they are generally supported by long batter piles that provide additional lateral and rotational stability for resisting flood loads.During Hurricane Katrina, several I-walls failed, resulting in large breaches and devastating flooding (Duncan et al. 2008, IPET 2009). T-walls on the other hand performed quite well, including T-walls that were overtopped and eroded on the protected side, yet still did not collapse. As part of the HSDRSS, more conservative designs and height restrictions were required for I-walls (USACE 2011), and a number of I-wall sections were replaced with T-walls. However, the design of pile-supported T-walls is complex and requires careful consideration of both flood loading and settlements."

Jan 1, 2017

By Zhou Haizuo, Diao Yu, Zheng Gang

"The conventional design method for the embankment supported by DM columns usually assumes the simultaneous failure of all the columns. However, the progressive failure may occur in practice, especially considering that many researches have justified the critical load of DM columns is bending moment. In this study, a finite difference method was adopted to simulate the embankment on soft clay improved by DM columns. An increasing surcharge was applied on the embankment surface to investigate the mechanism of the progressive failure of columns beneath embankment. To capture the brittle bending failure of the DM columns, the model employed a cut-off method, in which the column was separated on failed section as long as the ultimate bending failure was reached. Furthermore, high-strength DM columns were used to enhance the global stability of embankment in this study. The numerical results showed that the ultimate bending failure was first reached at the column beneath embankment toe, and then a continuous collapse extended to the column beneath the embankment shoulder. It is found that the reinforcement could be obviously effective when the reinforced area of high-strength DM columns was extended to a certain area.INTRODUCTIONDeep mixed (DM) columns have been commonly used as a soil-improvement technique to increase stability and to reduce settlement of embankment. Research on failure modes of DM columns owing to embankment fill and surcharge has been the subject of interest for a number of investigators (Miyake 1991; Kitazume 2000; Han 2005; 2006).The potential failure modes of DM column include bending, inclination, shearing, rotating, soil flowing and a combination of these modes (Kitazume and Terashi 1991; Broms 1999). Many researchers have found that bending failure of DM columns was more critical than the shear failure. A centrifuge test was performed by Kitazume (2000) and the results showed that the columns failed in bending before shear failure. In addition, several numerical computations illustrated that design method using the shear strength overestimated the stability of embankment over DM columns (Han and Huang 2005; 2006; Navin 2006; Liu 2007; Wang 1996). In practice, an equivalent shear strength method based on bending failure could better reflect the stability of embankment (Zheng 2010; China Code JGJ79-2011)."

Jan 1, 2015

By A. Manouchehrian, M. Cai

"Rockburst, which often occurs around underground openings in highly stressed hard rocks, is a sudden and violent unstable rock failure and is associated with energy release. Rockburst is very dangerous to workers and mining equipment, and it has been a major concern for deep underground construction and mining. Before excavation, rocks are loaded under a true-triaxial condition. Upon excavation, the rocks at the excavation boundaries are loaded in the maximum principal stress direction and unloaded in the minimum principal stress direction and this provides a necessary condition for unstable rock failure to occur. Understanding rock behavior under true-triaxial unloading condition is critical for developing measures to control unstable rock failure in underground construction.In this article, an explicit FEM tool is used to simulate rock behavior under true-triaxial unloading condition, with a focus on analyzing stable and unstable failure and the released energy during failure. Firstly, a uniaxial compression test was simulated. The results indicated that the adopted numerical tool was able to simulate stable and unstable failure by varying the Loading System Stiffness (LSS). In this research, Loading System Reaction Intensity (LSRI) was proposed as an indicator to identify stable and unstable failure. Secondly, unstable rock failure under true-triaxial unloading condition was simulated. A rock specimen was loaded in three directions; the lateral pressures (s3 and s2) were kept constant and the vertical stress (s1) was increased to a pre-defined stress. The specimen was suddenly unloaded in the minimum principal stress direction (s3) and loading in the s1 direction was continued until failure occurred. Different LSS values were considered and the results showed that the LSS influenced the failure mode; a soft loading system resulted in more violent failure and released more energy than that in a stiff loading system. The influence of specimen height to width ratio on rock failure was also investigated. It was found that failure in taller specimens was more violent than that in shorter specimens. Insights gained from the simulation are useful for better understanding of the influencing factors on the released kinetic energy during unstable rock failure and to find means to control rockburst damage in engineering practice."

Jan 1, 2015

By Dao Hong Quang

This paper describes the numerical analysis of some geotechnical factors influencing the application of Longwall Top Coal Caving (LTCC) method. The study was aimed at comparing the differences in stress distribution in the chain pillars of the face ends and caving mechanism of the top coal in LTCC method between the flat and inclined thick coal seams. The outcomes from this analysis will benefit mine planners in designing the location of roadway development in the next longwall panel. In addition, results from the study will extend the understanding of caving mechanism in top coal section when applying the LTCC method (with the direction of the longwall face is down dip, retreating along strike) in inclined coal seams. The results demonstrate a clear effect of dipping factor on cavability of the top coal and on the stress distribution in the chain pillar of the face ends.

Jan 1, 2008

By N. E. Mkadmi

Backfilling increases the stability of underground stopes and may help reduce the environmental footprint from the mining operation. The barricades holding the backfill in place must be designed to prevent failure based on the total and effective stress states, which depend on fill properties, stope geometry, and filling rate. This paper presents results from a numerical study aimed at analyzing the stress states associated with the response of paste backfill following stope filling, with an emphasis on the short term behavior. The simulations illustrate how stresses evolve with backfill drainage and filling rate in a vertical stope.

May 1, 2011

By Qingwen Shi, Brijes Mishra

The abstract discusses the use of numerical simulations to investigate roof fall in underground mines. It highlights the advantages of the discrete element method (DEM) in simulating large strain problems and fracturing processes, while acknowledging its computational limitations for large-scale problems. The study presents a 3D numerical model coupling the finite difference method (FDM) and DEM to simulate laminated roof failure in a coal mine. The model represents the laminated roof with bonded particles and discontinuities, while the surrounding rock is represented by continuum zones. The research analyzes the fracturing of the laminated roof as the entry advances, demonstrating the coupled method's ability to reproduce laminated roof failure.

Mar 17, 2022