• Journal of Korea Water Resources Association
• /
• v.31 no.3
• /
• pp.347-360
• /
• 1998

Numerical simulations for fluid flow and solute transport in a fracture rock masses are performed by using a transient flow model, which is based on the three-dimensional stochastic and discrete fracture network model (DFN model) and is coupled hydraulic model with mechanical model. In the numerical simulations of the solute transport, we used to the particle following algorithm which is similar to an advective biased random walk. The purpose of this study is to predict the response of the tracer test between two deep bore holes (GPK1 and GPK2) implanted at Soultz sous Foret in France, in the context of the geothermal researches.l The data sets used are obtained from in situcirculating experiments during 1995. As the result of the transport simulation, the mean transit time for the non reactive particles is about 5 days between two bore holes.

• Geomechanics and Engineering
• /
• v.21 no.3
• /
• pp.247-258
• /
• 2020

The shear strength of unsaturated soils, a research hotspot in geotechnical engineering, has great guiding significance for geotechnical engineering design. Although kinds of calculation models for the shear strength of unsaturated soil have been put forward by predecessors, there is still need for new models to extensively consider the nonlinear variation of shear strength, particularly for the nonlinear effect of the net normal stress on the shear strength of unsaturated soil. Here, the shear strength of unsaturated soils is explored to study the nonlinear effects of net normal stress with the introduction of a general nonlinear Mohr-Coulomb (M-C) strength criterion, and the relationship between the matric suction (or suction stress) and degree of saturation (DOS) constructed by the soil-water characteristics curve (SWCC) of van Genuchten is also applied for unsaturated soil. Then, two calculation models (i.e., an envelope shell model and an effective stress model) are established for the shear strength of unsaturated soils under the nonlinear strength theory. In these two models, the curve of the shear strength of unsaturated soils versus the net normal stress exhibits a tendency to gently. Moreover, the proposed formulas have flexibility and convenience with five parameters (for the effective stress model) or six parameters (for the envelope shell model), which are from the M-C strength parameters of the saturated soil and fitting parameters of SWCC of van Genuchten. Thereafter, by comparison with the classical theory of the shear strength of unsaturated soils from some actual cases, the rationality and accuracy of the present models were verified.

• Geomechanics and Engineering
• /
• v.24 no.4
• /
• pp.307-321
• /
• 2021

In an earlier publication (Serrano et al. 2014), the theoretical basis for evaluating the shear strength in rock joints was presented and used to derive an equation that governs the relationship between tangential and normal stresses on the joint during slippage between the joint faces. In this paper, the theoretical equation is applied to two non-linear failure criteria by using non-associated flow laws, including the modified Hoek and Brown and modified Mohr-Coulomb equations. The theoretical model considers the geometric dilatancy, the instantaneous friction angle, and a parameter that considers joint surface roughness as dependent variables. This model uses a similar equation structure to the empirical law that was proposed by Barton in 1973. However, a good correlation with the empirical values and, therefore, Barton's equation is necessary to incorporate a non-associated flow law that governs breakage processes in rock masses and becomes more significant in highly fractured media, which can be induced in a rock joint. A linear law of dilatancy is used to assess the importance of the non-associated flow to obtain very close values for different roughness states, so the best results are obtained for null material dilatancy, which considers significant changes that correspond to soft rock masses or altered zones of weakness.

• Geomechanics and Engineering
• /
• v.29 no.3
• /
• pp.281-290
• /
• 2022

The direct shear test is commonly used to evaluate the shear behavior of frozen soil-structure interfaces under normal stress. However, failure criteria, such as the Mohr-Coulomb failure criterion, are needed to obtain the unconfined shear strength. Hence, the punch shear test, which is usually used to estimate the shear strength of rocks without confinement, was examined in this study to directly determine the adfreezing strength. It is measured as the shear strength of the frozen soil-structure interface under unconfined conditions. Different soils of silica sand, field sand, and field clay were prepared inside the steel and concrete ring structures. Soil and ring structures were frozen at the target temperature for more than 24 h. A punch shear test was then conducted. The test results show that the adfreezing strength increased with a decrease in the target temperature and increase in the initial water content, owing to the increase in ice content. The adfreezing strength of field clay was the smallest when compared with the other soil specimens because of the large amount of unfrozen water content. The field sand with the larger normalized roughness showed greater adfreezing strength than the silica sand with a lower normalized roughness. From the experiment and analysis, the applicability of the punch shear test was examined to measure the adfreezing strength of the frozen soil-structure interface. To find a proper sample dimension, supplementary experiments or numerical analysis will be needed in further research.

• Geomechanics and Engineering
• /
• v.36 no.2
• /
• pp.167-181
• /
• 2024

The soil-structure relative stiffness is a key factor affecting the seismic response of underground structures. It is of great significance to study the soil-structure relative stiffness for the soil-structure interaction and the seismic disaster reduction of subway stations. In this paper, the dynamic shear modulus ratio and damping ratio of an inhomogeneous soft soil site under different buried depths which were obtained by a one-dimensional equivalent linearization site response analysis were used as the input parameters in a 2D finite element model. A visco-elasto-plastic constitutive model based on the Mohr-Coulomb shear failure criterion combined with stiffness degradation was used to describe the plastic behavior of soil. The damage plasticity model was used to simulate the plastic behavior of concrete. The horizontal and vertical relative stiffness ratios of soil and structure were defined to study the influence of relative stiffness on the seismic response of subway stations in inhomogeneous soft soil. It is found that the compression damage to the middle columns of a subway station with a higher relative stiffness ratio is more serious while the tensile damage is slighter under the same earthquake motion. The relative stiffness has a significant influence on ground surface deformation, ground acceleration, and station structure deformation. However, the effect of the relative stiffness on the deformation of the bottom slab of the subway station is small. The research results can provide a reference for seismic fortification of subway stations in the soft soil area.

• Economic and Environmental Geology
• /
• v.38 no.3 s.172
• /
• pp.305-317
• /
• 2005

This study addresses the preliminary results of rock slope stability analyses including hazard assessments for slope failure conducted on the selected sections of rural road cut slope which are about 4 km long. The study area is located in the Mt. Chuntae northeast of Busan and mainly composed of Cretaceous rhyolitic ash-flow tuff', fallout tuff, rhyolitc and andesite. The volcanic rock mass in the area has a number of discontinuities that produce a potentially unstable slope, as the present cut slope is more than 70 degrees in most of the slope sections. Discontinuity geometry data were collected at selected 8 scanline sections and analyzed to estimate important discontinuity geometry parameters to perform rock slope kinematic and block theory analyses. Kinematic analysis for plane sliding has resulted in maximum safe slope angles greater than $65^{\circ}$ for most of the discontinuities. For most of the wedges, maximum safe cut slope angles greater than $45^{\circ}$ were obtained. Maximum safe slope angles greater than 80" were obtained fur most of the discontinuities in the toppling case. The block theory analysis resulted in the identification of potential key blocks (type II) in the SL4, SL5, SL6 and SL8 sections. The chance of sliding taking place through a type ll block under a combined gravitational and external loading is quite high in the investigated area. The results support in-field observations of a potentially unstable slope that could become hazardous under external forces. The results obtained through limit equilibrium slope stability analyses show how a stable slope can become an unstable slope as the water pressure acting on joints increases and how a stable slope under Barton's shear strength criterion can fail as the worst case scenario of using Mohr-Coulomb criterion.

• Journal of the Korean Geotechnical Society
• /
• v.19 no.3
• /
• pp.53-64
• /
• 2003

Many transmission towers, high-rise buildings and bridges are constructed near steep slopes and are supported by large-diameter piles. These structures may be subjected to large lateral loads, such as violent winds and earthquakes. Widely used types of foundations for these structures are pier foundations, which have large-diameters with high stiffness. The behavior of a pier foundation subjected to lateral loads is similar to that of a short rigid pile because both elements seem to fail by rotation developing passive resistance on opposite faces above and below the rotation point, unlike the behavior of a long flexible pile. This paper describes the results of several numerical studies performed with a three-dimensional finite element method (FEM) of model tests of a laterally loaded short pile located near slopes, respectively. In this paper, the results of model tests of single piles and pile groups subjected to lateral loading, in homogeneous sand with 30$^{\circ}$ slopes and horizontal ground were analyzed by the 3-D FE analyses. The pile was assumed to be linearly elastic. The sand was assumed to have non-associative characteristics, following the MC-DP model. The failure criterion is governed by the Mohr-Coulomb equation and the plastic potential is given by the Drucker-Prager equation. The main purpose of this paper is the validation of the 3-D elasto-plastic FEM by comparisons with the experimental data.

• 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.221-242
• /
• 2021

A numerical study of the performance assesment of coupled thermo-hydro-mechanical (THM) processes in improved Korean reference disposal system (KRS⁺) for high-level radioactive waste is conducted using TOUGH2-MP/FLAC3D simulator. Decay heat from high-level radioactive waste increases the temperature of the repository, and it decreases as decay heat is reduced. The maximum temperature of the repository is below a maximum temperature criterion of 100℃. Saturation of bentonite buffer adjacent to the canister is initially reduced due to pore water evaporation induced by temperature increase. Bentonite buffer is saturated 250 years after the disposal of high-level radioactive waste by inflow of groundwater from the surrounding rock mass. Initial saturation of rock mass decreases as groundwater in rock mass is moved to bentnonite buffer by suction, but rock mass is saturated after inflow of groundwater from the far-field area. Stress changes at rock mass are compared to the Mohr-Coulomb failure criterion and the spalling strength in order to investigate the potential rock failure by thermal stress and swelling pressure. Additional simulations are conducted with the reduced spacing of deposition holes. The maximum temperature of bentonite buffer exceeds 100℃ as deposition hole spacing is smaller than 5.5 m. However, temperature of about 56.1% volume of bentonite buffer is below 90℃. The methodology of numerical modeling used in this study can be applied to the performance assessment of coupled THM processes for high-level radioactive waste repositories with various input parameters and geological conditions such as site-specific stress models and geothermal gradients.

• Economic and Environmental Geology
• /
• v.10 no.1
• /
• pp.37-48
• /
• 1977

The earth's crust is unceasingly undergoing deformations because of the forces acting upon it. The relationship between the tectonic forces and the resulting deformations are found from the states of stresses in the earth's crust induced by these forces. The study has been attempted to analyze the deformations of the overlying sedimentary layers, which are deformed by the magma intrusion along its lower boundary. The elastic model is constructed to analyze the geologic structures, by means of the theory of elasticity, and then the appropriate boundary conditions are given. The solution of the Airy stress function which satisfies the given boundary conditions is derived from the analytic method. The internal stress distributions of the deformed elastic model layer are portrayed by principal stress trajetories, and then the corresponding potential faults and joints systems are predicted from the Coulomb-Mohr failure criterion. The internal displacement distributions are shown by the calculated displacement components vectors, namely horizontal, vertical and net components. Results of the numerical calculations show the developments of some geologic structures as follows; (1) one set of shear joints and or two sets of shear joints which are oppisite directions, and one set of extension joints parallel to the ${\sigma}_1$ direction, (2) one set of high angle thrusts and normal faults, (3) symmetric fold; both limbs are dipping in opposite direction with low angle. The field work at the Wall-A San area, located near Jinju City, in southern Korea, had accomplished to compare the field structures with the predicted ones. The results of the comparison exhibits the developments of joint and fault systems satisfactorily consistent with each others. But the area does not show any type of folding, in spite of the intrusion of a granodiorite massif, this fact is one of the important features of the whole Kyungsang sedimentary basins of Mesozoic age distributed at the south-eastern parts of Korea. For this reason, it is thought that the magma intrusion had occurred with extremly low pressure. The geologic structures have been modified by the erosion and weathering throughout the geologic time, and the conditions of the sedimentary layers (width, thickness and radius of magma) are not the same as before, being intruded by the magma. To enlighten this, it is preferable to study these geologic structures with analyses of various types of rheological models.