جستجوی مقالات مرتبط با کلیدواژه « finite element method » در نشریات گروه « مهندسی معدن »

Determining the appropriate blasting pattern is important to prevent any damage to the tunnel perimeter in conventional tunneling by blasting operation in hard rocks. In this research work, the LS-DYNA software and numerical finite element method (FEM) are used for simulation of the blasting process in the Miyaneh-Ardabil railway tunnel. For this aim, the strong explosive model and nonlinear kinematic plastic material model are considered. Furthermore, the parameters required for the Johnson-Holmquist behavioral model are based on the Johnson-Holmquist-Ceramic material model relationships and are determined for the andesitic rock mass around studied tunnel. The model geometry is designed using AUTOCAD software and Hyper-mesh software is applied for meshing simulation. After introducing elements properties and material behavioral models and applying control and output parameters in LS-PrePost software, the modeling process is performed by LS-DYNA software. Different patterns of blastholes including 66, 23, and 19 holes, with diameters of 40 and 51 mm, and depths of 3 to 3.8 m are investigated by three-dimensional FEM. The borehole pressure caused by the ammonium nitrate-fuel oil (ANFO) detonation is considered based on the Jones-Wilkins-Lee (JWL) equation of state in the LS-DYNA software. The outer boundaries of the model are considered non-reflective to prevent the wave’s return. The results showed that LS-DYNA software can efficiently simulate the blasting process. Moreover, the post-failure rate of the blasting is reduced by more than 30% using the main charge with less explosive power and reducing the distance and diameter of contour holes.

This paper presents a comprehensive study on the stability of the deep underground closed Kolar Gold Fields mine (3.2 km deep) under varying seismic loading conditions. The study utilized the Finite Element Method (FEM)-based Midas GTS NX software tool to conduct numerical simulations of seismic loads of varying intensities under multiple conditions of water level in the mine voids. The seismic loads applied were equivalent to the intensity of maximum mining-induced seismicity experienced in the mine. The study also examined the influence of the Mysore North Fault and its effects on the surface above the mining area. A seismic hazard vulnerability map of the mining area was developed based on the results for all simulated numerical model combinations. The results inferred that for a seismic load of PGA, 0.22 g, for fault and actual water level combination, very strong shaking and moderate potential surface damage were observed at vulnerable zones with a maximum PGA of   0.196 g and Peak Ground Velocity (PGV) of 0.49 m/s. The study highlights the importance of monitoring post-mining induced seismic activities using a dedicated microseismic monitoring system with sensors placed at the most vulnerable zone locations assessed from the numerical modelling studies carried out. Remedial measures suggested include regular dewatering of mine workings based on water accumulation and backfilling of mine voids with suitable fill material. The dynamic modelling approach using Midas GTS NX was found to be a more reliable, feasible, efficient, and simple method for assessing the stability of closed mines.

Investigating the crack propagation mechanism is of paramount importance in analyzing the failure process of most materials. This process may be exposed during each kind of loading on the materials. In this work, the cracking mechanism in rock-like materials is studied using the numerical methods and compared with the experimental test results. However, the mechanism of crack growth in brittle materials such as rocks is influenced by different parameters. This research work focuses on the effect of the initial crack angles on the crack growth paths of these materials. Some cubic samples containing pre-existing cracks are tested in compression by considering different flaw orientations. The specimens are made of cement, water, and sand. Moreover, the mentioned process is numerically simulated using three different methods the finite difference method for discontinuous bodies or discrete element method, the displacement discontinuity method, and the versatile finite element method. The micro-parameters for simulation are gained by the trial-and-error procedure for the discrete element method. Eventually, the crack growth paths observed in the experiments are compared with the numerically simulated models. The results obtained show that these central cracks propagate in two ways, which are dependent on their initial angle. By increasing the initial crack angle to greater than 30° (α > 30°), the wing crack path moves further away from the initial crack, and by decreasing α to smaller than 30° (α < 30°), only the shear cracks are initiated. Therefore, the validity and accuracy of the results are manifested by comparing all the corresponding results obtained by different methods. Based on these results, it can generally be concluded that the strength of the cubic (rock material) specimens increases with increase in the crack angles with respect to the applied loading direction.

During an earthquake, the better performance of segmental tunnel lining, compared to the continuous in-cast concrete lining, is generally related to the joints between segments. In order to better understand the influence of the segment joints, their effect on the internal forces induced in tunnel lining simultaneously with the effects of the other influential parameters should be considered. In this work, the segmental joints were simulated by the representative stiffnesses and effects of these characteristics in relation to the other parameters such as the soil-liner interface behavior, number of segments in each ring and thickness of segments on the internal forces induced in structure were investigated. For this purpose, 2D numerical analyses were performed and the results obtained were discussed. Results showed that under the seismic condition, the components that had the most significant role on the internal axial forces induced in the segmental lining were rotational stiffness and axial stiffness of joints. Also the bending moments were more affected by the rotational stiffness. Generally, the radial joint stiffness had a less effect on the induced internal forces. With increase in the number of segments and their thickness, the effect of joint stiffness on the internal forces increases and the design of joints should be given more attention; however, the effects of joint stiffness and frictional behavior at the soil-liner interface on the maximum induced forces are almost independent from each other. Also in a specified joint behavior, by variation in each one of the other parameters including the soil-liner interface condition, number of segments and their thickness, the absolute magnitude of the maximum induced internal forces sometimes change significantly.

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.