• Journal of the Korean Geotechnical Society
• /
• v.35 no.5
• /
• pp.5-19
• /
• 2019

Considering the natural phenomenon in which steep slopes ($65^{\circ}{\sim}85^{\circ}$) consisting of rock mass remain stable for decades, slopes steeper than 1:0.5 (the standard of slope angle for blast rock) may be applied in geotechnical conditions which are similar to those above at the design and initial construction stages. In the process of analysing the stability of a good to fair continuum rock slope that can be designed as a steep slope, a general method of estimating rock mass strength properties from design practice perspective was required. Practical and genealized engineering methods of determining the properties of a rock mass are important for a good continuum rock slope that can be designed as a steep slope. The Genealized Hoek-Brown (H-B) failure criterion and GSI (Geological Strength Index), which were revised and supplemented by Hoek et al. (2002), were assessed as rock mass characterization systems fully taking into account the effects of discontinuities, and were widely utilized as a method for calculating equivalent Mohr-Coulomb shear strength (balancing the areas) according to stress changes. The concept of calculating equivalent M-C shear strength according to the change of confining stress range was proposed, and on a slope, the equivalent shear strength changes sensitively with changes in the maximum confining stress (${{\sigma}^{\prime}}_{3max}$ or normal stress), making it difficult to use it in practical design. In this study, the method of estimating the strength properties (an iso-angle division method) that can be applied universally within the maximum confining stress range for a good to fair continuum rock mass slope is proposed by applying the H-B failure criterion. In order to assess the validity and applicability of the proposed method of estimating the shear strength (A), the rock slope, which is a study object, was selected as the type of rock (igneous, metamorphic, sedimentary) on the steep slope near the existing working design site. It is compared and analyzed with the equivalent M-C shear strength (balancing the areas) proposed by Hoek. The equivalent M-C shear strength of the balancing the areas method and iso-angle division method was estimated using the RocLab program (geotechnical properties calculation software based on the H-B failure criterion (2002)) by using the basic data of the laboratory rock triaxial compression test at the existing working design site and the face mapping of discontinuities on the rock slope of study area. The calculated equivalent M-C shear strength of the balancing the areas method was interlinked to show very large or small cohesion and internal friction angles (generally, greater than $45^{\circ}$). The equivalent M-C shear strength of the iso-angle division is in-between the equivalent M-C shear properties of the balancing the areas, and the internal friction angles show a range of $30^{\circ}$ to $42^{\circ}$. We compared and analyzed the shear strength (A) of the iso-angle division method at the study area with the shear strength (B) of the existing working design site with similar or the same grade RMR each other. The application of the proposed iso-angle division method was indirectly evaluated through the results of the stability analysis (limit equilibrium analysis and finite element analysis) applied with these the strength properties. The difference between A and B of the shear strength is about 10%. LEM results (in wet condition) showed that Fs (A) = 14.08~58.22 (average 32.9) and Fs (B) = 18.39~60.04 (average 32.2), which were similar in accordance with the same rock types. As a result of FEM, displacement (A) = 0.13~0.65 mm (average 0.27 mm) and displacement (B) = 0.14~1.07 mm (average 0.37 mm). Using the GSI and Hoek-Brown failure criterion, the significant result could be identified in the application evaluation. Therefore, the strength properties of rock mass estimated by the iso-angle division method could be applied with practical shear strength.

• Tunnel and Underground Space
• /
• v.31 no.4
• /
• pp.270-288
• /
• 2021

We proposed a numerical method to simulate the hydro-mechanical behavior of rock fracture using a grain-based distinct element model (GBDEM) in the paper. As a part of DECOVALEX-2023 Task G, we verified the method via benchmarks with analytical solutions. DECOVALEX-2023 Task G aims to develop a numerical method to estimate the coupled thermo-hydro-mechanical processes within the crystalline rock fracture network. We represented the rock sample as a group of tetrahedral grains and calculated the interaction of the grains and their interfaces using 3DEC. The micro-parameters of the grains and interfaces were determined by a new methodology based on an equivalent continuum approach. In benchmark modeling, a single fracture embedded in the rock was examined for the effects of fracture inclination and roughness, the boundary stress condition and the applied pressure. The simulation results showed that the developed numerical model reasonably reproduced the fracture slip induced by boundary stress condition, the fracture opening induced by fluid injection, the stress distribution variation with fracture inclination, and the fracture roughness effect. In addition, the fracture displacements associated with the opening and slip showed good agreement with the analytical solutions. 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 experiments.

• Tunnel and Underground Space
• /
• v.31 no.6
• /
• pp.593-609
• /
• 2021

We proposed a numerical method for the thermo-mechanical behavior of rock fracture using a grain-based distinct element model (GBDEM) and simulated thermally induced fracture slip. The present study is the benchmark simulation performed as part of DECOVALEX-2023 Task G, which aims to develop a numerical method to estimate the coupled thermo-hydro-mechanical processes within the crystalline rock fracture network. We represented the rock sample as an assembly of Voronoi grains and calculated the interaction of the grains (blocks) and their interfaces (contacts) using a distinct element code, 3DEC. Based on an equivalent continuum approach, the micro-parameters of grains and contacts were determined to reproduce rock as an elastic material. Then, the behavior of the fracture embedded in the rock was characterized by the contacts with Coulomb shear strength and tensile strength. In the benchmark simulation, we quantitatively examined the effects of the boundary stress and thermal stress due to heat conduction on fracture behavior, focusing on the mechanism of thermally induced fracture slip. The simulation results showed that the developed numerical model reasonably reproduced the thermal expansion and thermal stress increment, the fracture stress and displacement and the effect of boundary condition. 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 experiments.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.22 no.8
• /
• pp.1833-1841
• /
• 1997

Rectangular dielectric waveguides, the most fundamental and indispensible elements in integrated optics, have been investigated by many researchers with various approaching methods including from the relatively approximate techniques to the numerical method. In this paper, the optimum equivalent waveguide model is adopted which is determined by a perturbation feedback process for analyzing the propagation constant by means of computer simulation, we have ascertained that the propagation constant from perturbation feedback method gives the best approximate value because it coincide with more exact value than obtained by other approximating methods. The technique also provides analytical expression for the modal field profile that should be useful in the design of various integrated optical devices.

• Magazine of the Korea Concrete Institute
• /
• v.7 no.5
• /
• pp.172-178
• /
• 1995

Concrete contains numerous microcracks at initially poured. The growth and propagation of nicrockacsk are believed tc finally incur the faiure of concrete. These processings are understood as a damage. Damage IS represented as a second-order tensor and crack is treated as a con tinuum phenomenon. In this paper, damage is characterized through the effective stress concept together with the hypothesis of elastic energy equivalence, and damage evolution law and constitutive equation of a damage model are derived by using the Helmholtz frte eriergy and the dissipation potential by means of the thermodynamic principles. The constitutive equation of the model includes the effects of elasticity, anisotropic damage and plasticity of concrete. There are two effective tangent stiffness tensors in this model : one is for elastic-darnage and the other for plastic damage. For the verification of the model, finite element analysis was performed for the analysis of concrete subjec:t to uniaxial and biaxial loading and the results obtained were compared with test results.

• Tunnel and Underground Space
• /
• v.25 no.5
• /
• pp.450-461
• /
• 2015

Extensive numerical experiments have been carried out to investigate effect of block size and fracture geometry on hydraulic characteristics of fractured rock masses based on connected pipe flow in DFN systems. Using two fracture sets, a total of 72 2-D fracture configurations were generated with different combinations of fracture size distribution and deterministic fracture density. The directional block conductivity including the theoretical block conductivity, principal conductivity tensor and average block conductivity for each generated fracture network system were calculated using the 2-D equivalent pipe network method. There exist significant effects of block size, orientation, density and size of fractures in a fractured rock mass on its hydraulic behavior. We have been further verified that it is more difficult to reach the REV size for the fluid flow network with decreasing intersection angle of two fracture sets, fracture plane density and fracture size distribution.

• Tunnel and Underground Space
• /
• v.26 no.4
• /
• pp.283-292
• /
• 2016

A computer program code was developed to estimate hydraulic behavior of the 2-D connected pipe network system, and implemented to evaluate the effect of joint aperture on hydraulic parameters of fractured rock masses through numerical experiments. A total of 216 stochastic 2-D DFN(discrete fracture network) blocks of $20m{\times}20m$ were prepared using two joint sets with fixed input parameters of joint orientation, frequency and size distribution. Two different cases of joint aperture variation are considered in this study. The hydraulic parameters were estimated for generated 2-D DFN blocks. The hydraulic anisotropy and the chance for equivalent continuum behavior of the DFN system were found to depend on the variability of joint aperture.

• Journal of the Korean Geotechnical Society
• /
• v.31 no.4
• /
• pp.45-55
• /
• 2015

The stability of the adjcent structures or slopes under blasting is typically evaluated using an empirical vibration attenuation curve or dynamic numerical analysis. To perform a dynamic analysis, it is necessary to determine the blast load and the damping ratio of rock mass. Various empirical methods have been proposed for the blast load. However, a study on representative values of damping ratio of a rock mass has not yet been performed. Therefore, the damping ratio was either ignored or selected without a clear basis in performing a blast analysis. Selection of the dampring ratio for the rock mass is very difficult because the vibration propagation is influenced by the layout and properties of the rock joints. Besides, the vibration induced by blasting is propagated spherically, whereas plane waves are generated by an earthquake. Since the geometrical spreading causes additional attenuation, the damping ratio should be adjusted in the case of a 2D plane strain analysis. In this study, we proposed equivalent damping ratios for use in continuum 2D plane strain analyses. To this end, we performed 2D dynamic analyses for a wide range of rock stiffness and investigated the characteristics of blast vibration propagation. Based on numerical simulations, a correlation between the attenuation equation, shear wave velocity, and equivalent damping ratio of rock mass is presented. This novel approach is the first attempt to select the damping ratio from an attenuation relationship. The proposed chart is easy to be used and can be applied in practice.

• Explosives and Blasting
• /
• v.24 no.2
• /
• pp.41-50
• /
• 2006

Excavation by explosives blasting necessarily involves noise and vibration, which is highly prone to face claims on the environmental and structural aspects from the neighbors. When the blasting carried out in the vicinity of a structure, the effect of blasting vibration on the stability of the structure should be carefully evaluated. In the conventional method of evaluation, an equation for blast vibration is obtained from test blasting which is later used to determine the amount of charge. This method, however, has limitations in use since it does not consider topography and change in ground conditions. In order to overcome the limitations, dynamic numerical analysis is recently used in continuum or discontinuous models, where the topography and the ground conditions can be exactly implemented. In the numerical analysis for tunnels and rock slopes, it is very uncommon to simulate multi-hole blasting. A single-hole blasting pressure is estimated and the equivalent overall pressure at the excavation face is used. This approach based on an ideal case usually does not consider the ground conditions. And this consequently results in errors in calculation. In this presentation of a case study, a new approach of using blast waves obtained in the test blast is proposed. The approach was carried out in order to improve the accuracy in calculating blasting pressure. The stability of a structure in the vicinity of a slope blasting was examined using the newly proposed method.

• Composites Research
• /
• v.37 no.2
• /
• pp.132-138
• /
• 2024

In this study, a molecular dynamics simulation study was performed to investigate the size-dependent electroelastic properties of single-walled boron nitride nanotubes(BNNT). To describe the elasticity and polarization of BNNT under mechanical loading, the Tersoff potential model and rigid ion approximation were adopted. For the prediction of piezoelectric constants and Young's modulus of BNNTs, piezoelectric constitutive equations based on the Maxwell's equation were used to calculate the strain-electric displacement and strain-stress relationships. It was found that the piezoelectric constants of BNNTs gradually decreases as the radius of the tubes increases showing a nonnegligible size effect. On the other hand, the elastic constants of the BNNTs showed opposites trends according to the equivalent geometrical assumption of the tubular structures. To establish the structure-property relationships, localized configurational change of the primarily bonded B-N bonded topology was investigated in detail to elucidate the BNNT curvature dependent elasticity.