• Smart Structures and Systems
• /
• v.30 no.1
• /
• pp.49-62
• /
• 2022

Experimental and discrete element methods were used to investigate the effects of triple joints lengths and triple joint angle on the failure behavior of rock mass under uniaxial compressive test. Concrete samples with dimension of 20 cm × 20 cm × 5 cm were prepared. Within the specimen, three imbedded joint were provided. The joint lengths were 2 cm, 4cm and 6 cm. In constant joint lengths, the angle between middle joint and other joints were 30°, 60°, 90°, 120° and 150°. Totally 15 different models were tested under compression test. The axial load rate on the model was 0.05 mm/min. Concurrent with experimental tests, the models containing triple joints, length and joint angle are similar to the experiments, were numerical by Particle flow code in two dimensions (PFC2D). Loading rate in numerical modelling was 0.05 mm/min. Tensile strength of material was 1 MPa. The results show that the failure behaviors of rock samples containing triple joints were governed by both of the angle and the length of the triple joints. The uniaxial compressive strengths (UCS) of the specimens were related to the fracture pattern and failure mechanism of the discontinuities. Furthermore, it was shown that the compressive behavior of discontinuities is related to the number of the induced tensile cracks which are increased by decreasing the joint length. Along with the damage failure of the samples, the acoustic emission (AE) activities are excited. There were only a few AE hits in the initial stage of loading, then AE hits rapidly grow before the applied stress reached its peak. In addition, every stress drop was accompanied by a large number of AE hits. Finally, the failure pattern and failure strength are similar in both methods i.e., the experimental testing and the numerical simulation methods.

• Proceedings of the Korean Society for Rock Mechanics Conference
• /
• 2000.09a
• /
• pp.105-115
• /
• 2000

As large underground projects such as radioactive waste disposal, hot water and heat storage, and geothermal energy become influential, the study, which consider all aspects of thermics, hydraulics and mechanics would be needed. Thermo Hydro-Mechanical coupling analysis is one of the most complex numerical technique because it should be implemented with the combined three governing equations to analyze the behavior of rock mass. In this study, finite element code, which is based on Biot's consolidation theory, was developed to analyze the thermo-hydro-mechanical coupling in continuum rock mass. To verify the implemented program, one-dimensional consolidation model under the isothermal and non-isothermal conditions was analyzed and was compared with the analytic solution. The parametric study on two-dimensional consolidation was also performed and the effects of several factors such as poisson's ratio and hydraulic anisotropy on rock mass behavior were investigated. In the future, this program would be revised to be used for analysis of general discontinuous media with incorporating discrete joint model.

• Tunnel and Underground Space
• /
• v.10 no.3
• /
• pp.355-365
• /
• 2000

As large underground projects such as radioactive waste disposal, hot water and heat storage, and geothermal energy become influential, the study, which consider all aspects of thermics, hydraulics and mechanics would be needed. Thermo-Hydro-Mechanical coupling analysis is one of the most complex numerical technique because it should be implemented with the combined three governing equations to analyze the behavior of rock mass. In this study, finite element code, which is based on Biot's consolidation theory, was developed to analyze the thermo-hydro-mechanical coupling in continuum rock mass. To verify the implemented program, one-dimensional consolidation model under the isothermal and non-isothermal conditions was analyzed and was compared with the analytic solution. The parametric study on two-dimensional consolidation was also performed and the effects of several factors such as poisson's ratio and hydraulic anisotropy on rock mass behavior were investigated. In the future, this program would be revised to be used for analysis of general discontinuous media with incorporating discrete joint model.

• Structural Engineering and Mechanics
• /
• v.89 no.2
• /
• pp.181-197
• /
• 2024

A numerical method is presented in this paper, for buckling analysis of thin arbitrary stiffened composite cylindrical shells under axial compression. The stiffeners can be placed inside and outside of the shell. The shell and stiffeners are operated as discrete elements, and their interactions are taking place through the compatibility conditions along their intersecting lines. The governing equations of motion are obtained based on Koiter's theory and solved by utilizing the principle of the minimum potential energy. Then, the buckling load coefficient and the critical buckling load are computed by solving characteristic equations. In this formulation, the elastic and geometric stiffness matrices of a single curved strip of the shell and stiffeners can be located anywhere within the shell element and in any direction are provided. Moreover, five stiffened composite shell specimens are made and tested under axial compression loading. The reliability of the presented method is validated by comparing its numerical results with those of commercial software, experiments, and other published numerical results. In addition, by using the ANSYS code, a 3-D finite element model that takes the exact geometric arrangement and the properties of the stiffeners and the shell into consideration is built. Finally, the effects of Poisson's ratio, shell length-to-radius ratio, shell thickness, cross-sectional area, angle, eccentricity, torsional stiffness, numbers and geometric configuration of stiffeners on the buckling of stiffened composite shells with various end conditions are computed. The results gained can be used as a meaningful benchmark for researchers to validate their analytical and numerical methods.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.14 no.4
• /
• pp.293-308
• /
• 2012

The shear behavior at the particle/surface interface such as rock joint can determine the mechanical behavior of whole structure. Therefore, a fundamental understanding of the mechanisms governing its behavior and accurately estimation of the interface strength is essential. In this paper, PFC, a numerical analysis program of discrete element method was used to investigate the effects of the surface roughness on interface strength. The surface roughness was characterized by smooth, intermediate, and rough surface, respectively. In order to investigate the effects of particle shape and crushing on particle/surface interface behavior, one ball, clump, and cluster models were created and their results were compared. The shape of particle was characterized by circle, triangle, square, and rectangle, respectively. The results showed that as the surface roughness increases, interface strength and friction angle increase and the void ratio increases. The one ball model with smooth surface shows lower interface strength and friction angle than the clump model with irregular surface. In addition, a cluster model has less interface strength and friction angle than the clump model. The failure envelope of the cluster model shows non-linear characteristic. From these findings, it is verified that the surface roughness and particle shape effect on the particle/surface interface shear behavior.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.14 no.4
• /
• pp.321-336
• /
• 2012

The shear behavior at the particle/surface interface such as rock joint can determine the mechanical behavior of whole structure. Therefore, a fundamental understanding of the mechanisms governing its behavior and accurately estimation of the interface strength is essential. In this paper, PFC, a numerical analysis program of discrete element method was used to investigate the effects of the surface roughness crushing on interface strength. The surface roughness was characterized by smooth, intermediate, and rough surface, respectively. Particle shape was classified into one ball model of circular shape and 3 ball model of triangular shape. The surface shape was modelled by wall model of non-crushing surface and ball model of crushing surface. The results showed that as the bonding strength of ball model decreases, lower interface strength is induced. After the surface roughness crushing was occurred, the interface strength tended to converge and higher bonding strength induced lower surface roughness crushing. Higher friction angle was induced in wall model and higher surface roughness induced the higher friction angle. From these findings, it is verified that the surface roughness and surface roughness crushing effect on the particle/surface interface shear behavior.

• Geomechanics and Engineering
• /
• v.14 no.6
• /
• pp.545-552
• /
• 2018

Blast-induced ground vibrations by a significant amount of explosives may cause many problems for mining slope stability. Geological discontinuities have a significant influence on the transmission of dynamic pressure of detonation and according to their position relative to the slope face may have damaging or useful impacts on the slope stability. In this study, the effect of geological discontinuities was investigated by modelling a slope with geological discontinuities through applying the dynamic pressure in three-dimensional discrete element code (3DEC). The geological discontinuities in four states that generally apperceived in mine slopes are considered. Given the advantages of the pressure decay function defined by some researcher, this type of function was used to develop the pressure-time profile. The peak particle velocities (PPV) values were monitored along an axis by utilization of Fish programming language and the results were used as an indicator to measure the effects. As shown in the discontinuity-free model, PPV empirical models are reliable in rocks lacking discontinuities or tightly jointed rock masses. According to the other results, the empirical models cannot be used for the case where the rock mass contains discontinuities with any direction or dip. With regard to PPVs, when the direction of discontinuities is opposite to that of the slope face, the dynamic pressure of detonation is significantly damped toward the slope direction at the surface of discontinuities. On the other hand, when the discontinuities are horizontal, the dynamic pressure of detonation affects the rock mass to a large distance.

• Explosives and Blasting
• /
• v.23 no.3
• /
• pp.27-42
• /
• 2005

An explosion modeling technique was developed by using the spherical discrete element code, $PFC^{3D}$, which can be used to model the dynamic stress wave propagation phenomenon. The modeling technique is simply based on an idea that the explosion pressure should be applied to a $PFC^{3D}$ particle assembly not in the form of an external force (body force), but in the form of a contact force (surface force). The stress wave propagation modeling was conducted by simulating the experimental approach based on the Hopkinson's effect combined with the spatting phenomenon that had previously been developed to determine the dynamic tensile strength of Inada granite. As a result, the stress wave velocity obtained by the proposed modeling technique was 4167 m/s, which is merely $3\%$ lower than the actual wave velocity of 4300 m/s for an Inada granite.

• Geomechanics and Engineering
• /
• v.11 no.3
• /
• pp.401-420
• /
• 2016

The stability of deep coal roadways with large sections and thick top coal is a typical challenge in many coal mines in China. The innovative Universal Discrete Element Code (UDEC) trigon block is adopted to create a numerical model based on a case study at the Dongtan coal mine in China to better understand the failure mechanism and stability control mechanism of this kind of roadway. The failure process of an unsupported roadway is simulated, and the results suggest that the deformation of the roof is more serious than that of the sides and floor, especially in the center of the roof. The radial stress that is released is more intense than the tangential stress, while a large zone of relaxation appears around the roadway. The failure process begins from partial failure at roadway corners, and then propagates deeper into the roof and sides, finally resulting in large deformation in the roadway. A combined support system is proposed to support roadways based on an analysis of the simulation results. The numerical simulation and field monitoring suggest that the availability of this support method is feasible both in theory and practice, which can provide helpful references for research on the failure mechanisms and scientific support designing of engineering in deep coal mines.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2010.03a
• /
• pp.314-323
• /
• 2010

Recently, granular soils having a large particle size are frequently used as a filling material in the construction of foundation, harbor, dam, and so on. The shear behavior of this granular soil plays a key role in the stability of structures. For example, soil particle crushing occurring at the interface between structure and soil and/or within soil mass can cause the disturbance of ground characteristics and consequently induce an issues in respect of stability of structures. In order to investigate the shear behavior according to an existence and nonexistence of particle crushing, numerical analyses were conducted by using the DEM(Discrete Element Method)-based software program PFC(Particle Flow Code). Using the crushing model and non-crushing model which were created in this study, numerical analyses of ring shear test were conducted and their results were analyzed and compared. In general, landslide and slope stability are accompanied by a large displacement and consequently not only a peak strength but also a residual strength are very important in the analysis of landslide and slope stability. However the direct shear test which has been commonly used in the determination of shear strength parameters has a limitation on displacement therefore the residual strength parameters can not be obtained. The characteristics of residual shear behavior were investigated through the numerical analyses in this study.