• Tunnel and Underground Space
• /
• v.32 no.6
• /
• pp.568-585
• /
• 2022

In the present study, we proposed a numerical method for simulating thermally induced fracture slip using a grain-based distinct element model (GBDEM). As a part of DECOVALEX-2023, the thermo-mechanical loading test on a saw-cut rock fracture conducted at the Korea Institute of Civil Engineering and Building Technology was simulated. In the numerical model, the rock sample including a saw-cut fracture was represented as a group of random Voronoi polyhedra. Then, the coupled thermo-mechanical behavior of grains and their interfaces was calculated using 3DEC. The key concerns focused on the temperature evolution, thermally induced principal stress increment, and fracture normal and shear displacements under thermo-mechanical loading. The comparisons between laboratory experimental results and the numerical results revealed that the numerical model reasonably captured the heat transfer and heat loss characteristics of the rock specimen, the horizontal stress increment due to constrained displacement, and the progressive shear failure of the fracture. However, the onset of the fracture slip and the magnitudes of stress increment and fracture displacement showed discrepancies between the numerical and experimental results. We expect the numerical model to be enhanced by continuing collaboration and interaction with other research teams of DECOVALEX-2023 Task G and validated in further study.

• Tunnel and Underground Space
• /
• v.33 no.3
• /
• pp.189-207
• /
• 2023

In the present study, the thermoshearing experiment on a rough rock fracture were modeled using a three-dimensional grain-based distinct element model (GBDEM). The experiment was conducted by the Korea Institute of Construction Technology to investigate the progressive shear failure of fracture under the influence of thermal stress in a critical stress state. The numerical model employs an assembly of multiple polyhedral grains and their interfaces to represent the rock sample, and calculates the coupled thermo-mechanical behavior of the grains (blocks) and the interfaces (contacts) using 3DEC, a DEM code. The primary focus was on simulating the temperature evolution, generation of thermal stress, and shear and normal displacements of the fracture. Two fracture models, namely the mated fracture model and the unmated fracture model, were constructed based on the degree of surface matedness, and their respective behaviors were compared and analyzed. By leveraging the advantage of the DEM, the contact area between the fracture surfaces was continuously monitored during the simulation, enabling an examination of its influence on shear behavior. The numerical results demonstrated distinct differences depending on the degree of the surface matedness at the initial stage. In the mated fracture model, where the surfaces were in almost full contact, the characteristic stages of peak stress and residual stress commonly observed in shear behavior of natural rock joints were reasonably replicated, despite exhibiting discrepancies with the experimental results. The analysis of contact area variation over time confirmed that our numerical model effectively simulated the abrupt normal dilation and shear slip, stress softening phenomenon, and transition to the residual state that occur during the peak stress stage. The unmated fracture model, which closely resembled the experimental specimen, showed qualitative agreement with the experimental observations, including heat transfer characteristics, the progressive shear failure process induced by heating, and the increase in thermal stress. However, there were some mismatches between the numerical and experimental results regarding the onset of fracture slip and the magnitudes of fracture stress and displacement. This research was conducted as part of DECOVALEX-2023 Task G, and we expect the numerical model to be enhanced through continued collaboration with other research teams and validated in further studies.

• The Journal of Engineering Geology
• /
• v.15 no.3
• /
• pp.337-348
• /
• 2005

Jurassic granite from Pocheon area were tested to investigate the effect of microcracks on mechanical properties of the granite. Three oriented core specimens were used for uniaxial compressive tests and each core specimen are perpendicular to the axes'R'(rift plane),'c'(grain plane) and'H'(hardway plane), respectively Among vacious elastic constants, the variation of Poisson's ratio as function of the directions was examined. From the related chart between ratio of failure strength and Poisson's ratio, H-specimen shows the highest range in Poisson's ratio and Poisson's ratio decreases in the order of C-specimen and R-specimen. The curve pattern is nearly linear in stage $I\simIII$ but the slope increases abruptly in stage H-3. As shown in the related chart, diverging point of a curve is formed when ratio of failure strength is $0.92\sim0.96$ Stage IV -3 is out of elastic region. The behaviour of rock in the four fracturing stages was analyzed in term of the stress-volumetric strain me. From the stress increment-volumetric strain equations governing the behaviour of rock, characteristic material constants, a, n, Q, m and $\varepsilon_v^{mcf}$, were determined. Among these, inherent microcrack porosity$(a, 10^{-3})$ and compaction exponent(n) in the microcrack closure region(stage I ) show an order of $a^R(3.82)>a^G(3.38)>a^H(2.32)\;and\;n^R(3.69)>n^G(2.79)>n^H(1.99)4, respectively. Especially, critical volumetric microcrack strain($\varepsilon_v^{mcf}$) in the stage W is highest in the H-specimen, normal to the hardway plane. These results indicate a strong correlation between two major sets of microcracks and mechanical properties such as Poisson's ratio and material constants. Correlation of strength anisotropy with microcrack orientation can have important application in rock fracture studies.

• Tunnel and Underground Space
• /
• v.20 no.3
• /
• pp.217-224
• /
• 2010

Three relatively homogeneous granitic rocks were studied to investigate the relationship between their microstructural properties and fracture toughness. Fracture toughness and ultrasonic velocity were varied with the orientation of mineral's long axis and microcrack, obtained from optical microscope. The lowest fracture toughness values are obtained, when the fracture propagates parallel to weakness planes which have the orientation of mineral's long axis and microcrack, in other words, when weakness planes develop perpendicular to the direction of tensile stress agrees with that of rift plane. The fracture toughness values, measured with the short rod method, varied from 1.63 to 2.62 MPa $m^{0.5}$, and their values are related with the average grain size and average microcrack length.

• Computers and Concrete
• /
• v.17 no.5
• /
• pp.639-653
• /
• 2016

The present article discusses the effect of the ratio of bridge surface to total shear surface, number of bridge areas and normal stress on the failure behavior of the planar non-persistent open joints. Totally, 38 models were prepared using plaster and dimensions of $15cm{\times}15cm{\times}15cm$. The bridge area occupied $45cm^2$, $90cm^2$ and $135cm^2$ out of the shear surface. The number of rock bridges increase in fixed area. Two similar samples were prepared on every variation in the rock bridges and tested for direct shear strength under two high and low normal loads. The results indicated that the failure pattern and the failure mechanism is mostly influenced by the ratio of bridge surface to total shear surface and normal stress so that the tensile failure mode change to shear failure mode by increasing in the value of introduced parameters. Furthermore, the shear strength and shear stiffness are closely related to the ratio of bridge surface to total shear surface, number of bridge areas and normal stress.

• Tunnel and Underground Space
• /
• v.20 no.2
• /
• pp.125-130
• /
• 2010

Model blast experiments of mortar blocks were performed to investigate the effect of the guide hole on crack growth. The mortar block specimens have a blast charge hole and 8 guide holes. Two of circular guide hole, notched guide hole, diamond shaped guide hole and diamond shaped guide holder are installed around 110 mm, 165 mm and 220 mm apart from the charge hole for each specimen. From the blast experiments, it was revealed that all the guide hole used in this study were effective for controlling the crack growth at the fracture control.

• Geomechanics and Engineering
• /
• v.19 no.5
• /
• pp.421-434
• /
• 2019

It is of great significance to study the mechanical properties and failure mechanism of the defected rock for geological engineering. The defected sandstone modeling with power-law distribution of pre-cracks was built in this paper by Particle Flow Code software. Then the mechanical properties of sandstone and the corresponding failure process were meticulously analyzed by changing the power-law index (PLI) and the number of pre-cracks (NPC). The results show that (1) With the increase of the PLI, the proportion of prefabricated long cracks gradually decreases. (2) When the NPC is the same, the uniaxial compressive strength (UCS) of sandstone increases with the PLI; while when the PLI is the same, the UCS decreases with the NPC. (3) The damage model of rock strength is established based on the Mori-Tanaka method, which can be used to better describe the strength evolution of damaged rock. (4) The failure mode of the specimen is closely related to the total length of the pre-crack. As the total length of the pre-crack increases, the failure intensity of the specimen gradually becomes weaker. In addition, for the specimens with the total pre-crack length between 0.2-0.55 m, significant lateral expansion occurred during their failure process. (5) For the specimens with smaller PLI in the pre-peak loading process, the concentration of the force field inside is more serious than that of the specimens with larger PLI.

• Tunnel and Underground Space
• /
• v.19 no.2
• /
• pp.146-157
• /
• 2009

An experimental study of accelerated weathering on mudstone sample specimens from Haman, Gyeongnam was performed to investigate the variations of physico-mechanical properties of deteriorated rocks due to freeze-thaw weathering. Each complete cycle of freeze and thaw lasted 24 hours, comprising 2 hours of saturating in vacuum chamber, 8 hours of freezing at $-16{\pm}1^{\circ}C$ and 14 hours of thawing at room temperature. Total of 55 cycles of freeze-thaw were completed with measuring the index properties as well as geometries of microfractures. The measured specific gravity and P-wave velocity found to decrease with increasing freeze-thaw cycles. On the other hand, absorption ratio and effective porosity were continuously increased with increasing freeze-thaw cycles. It was found that the index properties of deteriorated sample specimen depend on its initial properties and flaws in rock. The size and density of the traces of the microfracture on slab specimen exhibited abrupt changes after 30 cycles of freeze-thaw weathering. The results obtained in this study show that the box fractal dimension($D_B$) given in this paper has the strong capability of quantifying the size and density of the microfracture.

• Smart Structures and Systems
• /
• v.22 no.6
• /
• pp.663-673
• /
• 2018

In this study, the effect of particle size on the cracks propagation and coalescence or cracking pattern of the edge notched disc specimens are investigated. Firstly, calibration of PFC3D was performed using Brazilian experimental test output. Then micro parameters were used to build edge notched disc specimen. The horizontal wall of the assembly is let to move downward with a standard low speed of 0.016 m/s. The numerical results show that the tensile cracks are dominant failure pattern for the modeled discs. These tensile cracks initiate from the pre-existing notch tip and propagate parallel to the loading direction then interact with the upper boundary of the modeled specimen. As the size of the balls (ball diameter) decrease the number of tensile cracks increase. The tensile fracture toughness of the samples also decreases as the particle size increases. Understanding the crack propagation and crack coalescence phenomena in brittle materials such as concretes and rocks is of paramount importance in the stability analyses for engineering structures such as rock slopes, underground structures and tunneling.

• Tunnel and Underground Space
• /
• v.19 no.6
• /
• pp.558-566
• /
• 2009

Time-dependent behavior is a basic mechanical property of rocks. Predicting the failure time of rock structures by analyzing the time-dependent characteristic is important and problematic. It is tried to predict the failure time of tunnel, slope & laboratory creep test specimen from measured displacement(or strain) and rate with relationship suggested by Voight($\ddot{\Omega}=A\dot{\Omega}^\alpha$, where $\Omega$ is a measurable quantity such as strain & displacement and A & $\alpha$ are constants). A & $\alpha$ are estimated through applying the nonlinear least square method to the single and double integrated Voight's equations and utilized to predict the failure time. Predicted failure time is in accordance with real one except minor error. Linear inverse rate method applied to creep strain and rate yields a poor linear correlation of data and precision of predicted failure time is not better than methods using strain and rate.