• Journal of Korean Tunnelling and Underground Space Association
• /
• v.12 no.1
• /
• pp.31-41
• /
• 2010

This paper presents both numerical and physical modeling approaches for the dilation and spalling of rock recognized as typical process of rock around an underground opening at depth. For physical approach, laboratory testing of rectangular beams using a synthetic rock was used to investigate the onset of dilation and spalling. The beams are axially compressed and subjected to 4-point bending to provide non-uniform compressive stresses which are similar to the maximum tangential stress distribution around circular openings. Discrete element numerical analyses using commercial code $PFC^{2D}$ (Particle Flow Code) were performed to evaluate the stress path at various locations in the beams. The findings from these approaches suggest that the onset of dilation in laboratory tests appears to be a good indicator for assessing the stress magnitudes required to initiate spalling.

• International Journal of Highway Engineering
• /
• v.16 no.2
• /
• pp.43-52
• /
• 2014

PURPOSES : To solve problems in current compaction control DCPT(Dynamic Cone Penetrometer Test), highly correlated with various testing methods, simple, and economic is being applied. However, it、s hard to utilize DCPT results due to the few numerical analyses for DCPT have been performed and the lack of data accumulation. Therefore, this study tried to verify the validation of numerical modeling for DCPT by comparing and analyzing the results of numerical analyses with field tests. METHODS: The ground elastic modulus and PR(Penetration Rate) value were estimated by using PFC(Particle Flow Code) 3D program based on the discrete element method. Those values were compared and analyzed with the result of field tests. Also, back analysis was conducted to describe ground elastic modulus of field tests. RESULTS : Relative errors of PR value between the numerical analyses and field tests were calculated to be comparatively low. Also, the relationship between elastic modulus and PR value turned out to be similar. CONCLUSIONS : Numerical modeling of DCPT is considered to be suitable for describing field tests by carrying out numerical analysis using PFC 3D program.

• Structural Engineering and Mechanics
• /
• v.66 no.4
• /
• pp.445-452
• /
• 2018

In this paper the effects of particle size and model scale of concrete have been investigated on point load index, tensile strength, and the failure processes using a PFC2D numerical modeling study. Circular and semi-circular specimens of concrete were numerically modeled using the same particle size, 0.27 mm, but with different model diameters of 75 mm, 54 mm, 25 mm, and 12.5 mm. In addition, circular and semi-circular models with the diameter of 27 mm and particle sizes of 0.27 mm, 0.47 mm, 0.67 mm, 0.87 mm, 1.07 mm, and 1.27 mm were simulated to determine whether they can match the experimental observations from point load and Brazilian tests. The numerical modeling results show that the failure patterns are influenced by the model scale and particle size, as expected. Both Is(50) and Brazilian tensile strength values increased as the model diameter and particle sizes increased. The ratio of Brazilian tensile strength to Is(50) showed a reduction as the particle size increased but did not change with the increase in the model scale.

• Computers and Concrete
• /
• v.22 no.5
• /
• pp.469-479
• /
• 2018

The Brazilian tensile strength of concrete samples is a key parameter in fracture mechanics since it may significantly change the quality of concrete materials and their mechanical behaviors. It is well known that porosity is one of the most often used physical indices to predict concrete mechanical properties. In the present work the influence of porosity shape on concrete tensile strength characteristics is studied, using a bonded particle model. Firstly numerical model was calibrated by Brazilian experimental results and uniaxial test out puts. Secondly, Brazilian models consisting various pore shapes were simulated and numerically tested at a constant speed of 0.016 mm/s. The results show that pore shape has important effects on the failure pattern. It is shown that the pore shape may play an important role in the cracks initiation and propagation during the loading process which in turn influence on the tensile strength of the concrete samples. It has also been shown that the pore size mainly affects the ratio of uniaxial compressive strength to that of the tensile one in the simulated material samples.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2009.03a
• /
• pp.58-63
• /
• 2009

Biaxial compression test was conducted on a transversely isotropic synthetic jointed rock model for the understanding of the fracture behaviors of a sedimentary or metamorphic rocks with well developed bedding or foliation in uni-direction. The joint angles employed for the model are 30, 45, and 60 degrees to the horizontal, and the synthetic rock mass was made of early strength cement. From the biaxial compression test, initiation propagation of tensile cracks at norm to the joint angle was found. The propagated tensile cracks eventually developed rock blocks, which was dislodged from the rock mass. Furthermore, the propagation process of the tensile cracks varies with joint angle: lower joint angle model shows more stable and progressive tensile crack propagation. The experiment results were validated from the simulation by using discrete element method PFC 2D. From the simulation, as has been observed from the test, a rock mass with lower joint angle produces wider damage region and rock block by tensile cracks. In addition, a rock model with lower joint angle shows a progressive tensile cracks generation around the opening from the investigation of the interacted tensile cracks.

• Smart Structures and Systems
• /
• v.22 no.4
• /
• pp.459-468
• /
• 2018

In this paper, the effect of non-persistent joints was determined on the behavior of concrete specimens subjected to biaxial loading through numerical modeling using particle flow code in two dimensions (PFC2D). Firstly, a numerical model was calibrated by uniaxial, Brazilian and triaxial experimental results to ensure the conformity of the simulated numerical model's response. Secondly, sixteen rectangular models with dimension of 100 mm by 100 mm were developed. Each model contains two non-persistent joints with lengths of 40 mm and 20 mm, respectively. The angularity of the larger joint changes from $30^{\circ}$ to $90^{\circ}$. In each configuration, the small joint angularity changes from $0^{\circ}$ to $90^{\circ}$ in $30^{\circ}$ increments. All of the models were under confining stress of 1 MPa. By using of the biaxial test configuration, the failure process was visually observed. Discrete element simulations demonstrated that macro shear fractures in models are because of microscopic tensile breakage of a large number of bonded discs. The failure pattern in Rock Bridge is mostly affected by joint overlapping whereas the biaxial strength is closely related to the failure pattern.

• Journal of Industrial Technology
• /
• v.30 no.A
• /
• pp.17-23
• /
• 2010

Recently, a large number of open-pit mines are planning to change their mining method to underground types because the environmental concerns and legal regulations are increased with a rise in the standard of living. The K silica mine, which is one of them and located in Kyunggi province, is planning the establishment of a vertical shaft which will be used for ore-pass channel in their new glory hole mining method. This vertical shaft will be designed to join with a horizontal gangway excavated from the ground level. In this new mining system, the excavated ore particles will be stored inside a shaft and transported out with a help of a conveyor belt. Therefore the hang-up of ore particles in a shaft, the control of gate at the bottom of a shaft, the installation of dog-leg at the gate should be investigated identically. In this study, the PFC-2D code which is one of the discrete element numerical methods has been applied to simulate the particle flow mechanism in a shaft, and the optimum mine design has been proposed to maximize the productivity and to minimize the system damage.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2009.09a
• /
• pp.591-604
• /
• 2009

This paper presents some important issues in modeling rock behaviour around an underground opening at depth which characterized as stress-induced fractural failure of rock. Unlike other conventional modeling approaches, stress-induced rock failure is highly complex process due to its own heterogeneous and discrete natures. Because of this complexity, many researchers has been struggled to mimic such processes as close as possible to reality with various approaches in both analytical, and numerical approaches for past few decades. Such approaches which are based on continuum mechanics, analytical fracture mechanics, and DEM(Discrete Element Method) were explored in this paper, and fundamental shortcomings for each approaches were illustrated here. In addition, DEM approach using $PFC^{2D}$(Particle Flow Code) was also implemented and illuminated in this paper and discuss the improvement and considerations for the future research.

• Explosives and Blasting
• /
• v.28 no.2
• /
• pp.1-16
• /
• 2010

This study focuses on the phenomenon that the blast damaged zone developed on rock slope surfaces can be affected by joint characteristics rather than by explosive power when the pre-splitting is applied to excavate a jointed rock slope. The characteristics of rock joints on a slope were investigated and categorized them into 4 cases. Also an image processing system has been used for comparing the distribution pattern of rock blocks. From this investigation, it was found that the rock blocks bigger than 2,000 mm occupied 42% in the case of single joint set and it showed the well efficiency of pre-splitting blast. In cases of 2~3 parallel joint sets and 2~3 intersecting joint sets are developed on rock surfaces, the rock blocks in the range of 1,000~2,000 mm occupied 43.6% and 35.8%, respectively, and it showed that the efficiency of pre-splitting was decreased. When more than 3 joint sets are randomly developed, however, the rock blocks in the range of 250~500 mm occupied 35% and there was no block bigger than 1,000 mm. This denotes that the blasting with pre-splitting was not effective. The numerical analysis using PFC2D showed that the blast damaged zone in a rock mass could be directly influenced by the pre-splitting. It is, therefore, required to investigate the discontinuity pattern on rock surfaces in advance, when the pre-splitting method is applied to excavate jointed rock slopes and to apply a flexible blating design with a consideration of the joint characteristics.

• Geomechanics and Engineering
• /
• v.7 no.6
• /
• pp.665-678
• /
• 2014

In discrete element modeling, 2D software has been widely used in order to gain further insights into the fundamental mechanisms with less computational time. The porosities used in 2D DEM studies should be determined with appropriate approaches based on 3D laboratory porosities. This paper summarizes the main approaches for converting porosities from 3D to 2D for DEM studies and theoretical evaluations show that none of the current approaches can be widely used in dealing with soil mechanical problems. Therefore, a parabolic equation and a criterion have been suggested for the determination of 2D porosities in this paper. Moreover, a case study has been used to validate that the 2D porosity obtained from the above suggestion to be rational with both the realistic contact force distribution in the specimen and the good agreement of the DEM simulation results of direct shear tests with the corresponding experimental data. Therefore, the parabolic equation and the criterion are suggested for the determination of 2D porosities in a wide range of polydisperse particle systems, especially in dealing with soil mechanical problems.