• Journal of Korean Tunnelling and Underground Space Association
• /
• v.18 no.6
• /
• pp.569-579
• /
• 2016

In order to ensure the stability of tunnel portal slope, reinforcement method such as anchors, soil nails and rock bolts have been used in Korea. When selecting slope reinforcement methods in tunnel portal area such as reinforcement arrangement and length, trial and error method can be very time-consuming and it was also not easy to verify the selection of an optimum condition. In this study, using the FISH language embedded in the finite difference code FLAC3D program, the optimization technique was developed with the Differential Evolution Algorithm (DEA). After building a database on the soil nailing method in tunnel portal area, this system can be selected to an optimum arrangement plan based on the factor of safety through the FLAC3D analysis. Through the results of numerical analysis, it was confirmed that the number of analysis was decreased by about 8 times when DEA based optimization technique was used compared to the full combination (FC). In case of the design of slope reinforcement in tunnel portal area, if this built-system is used, it is expected that the selection of an optimum arrangement plan can be relatively easier.

• Tunnel and Underground Space
• /
• v.14 no.5
• /
• pp.327-338
• /
• 2004

The purpose of this research is to establish a concept of the factor of safety of tunnels which is a quantitative estimate of the stability of tunnels. Based on this concept, a numerical technique which calculates the factor of safety of tunnels was developed. To obtain the safety factor of a tunnel, the strength reduction technique in which a series of analyses are repeated with reduced ground strength until the tunnel collapses were employed. With this technique, the failure plane, as well as the factor of safety, can be obtained without prescribing the trial failure plane. Analyses were conducted with FLA $C^{2D}$(ver3.3), a geotechnical analysis program which is based on the finite difference method. From the result, the location of plastic zones, the maximum convergence and the maximum stress generated in the support system were also analyzed. The result shows that factors of safety are higher for the 1st and 2nd rock classes, and lower for the lower rock classes. Furthermore, factor of safety is higher when $K_{0}$ =0.5 compared to at in case of $K_{0}$ =2.0. Through this research, it is found that the factor of safety defined in this research can be used as a good quantitative index representing the stability of tunnels. Also, close examination of the results can help adjustment of the quantity and location of additional supports.s.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.19 no.6
• /
• pp.614-627
• /
• 2018

Urban railway lines have been constructed adjacent to residential buildings and urban areas. The expansion of transportation networks and reconstruction of residential buildings in highly populated urban areas require deep excavations in areas adjacent to urban railways. Mobilized soil stresses and changes in the groundwater level induced by deep excavations results in track irregularities in urban railways. In this study, a three-dimensional finite difference model using the commercial program FLAC3D was adopted to estimate the horizontal displacements of earth retaining structures, settlements of backfill, the stability of track irregularity and underground box structure based on the criteria of each railway organization and its relationships. In deep excavations, a change in groundwater level induces relatively very small differences for track gauge irregularities, whereas relatively large differences for longitudinal irregularities of 72.5%, twist irregularities of 83.3%, cross level irregularities of 61.9%, and alignments of 43.3% were found to be the maximum differences when the horizontal displacement of earth retaining wall and settlement of backfill were 65.1% and 21.4%, respectively, because the groundwater level (GWL) on the ground surface-mobilized tensile strength of the underground box structure exceeds the allowable value. Therefore, three-dimensional numerical analysis was performed in this study. Overall, real-time monitoring should be carried out to prevent railway accidents in advance when a deep excavation adjacent to urban railway structures is constructed.

• Journal of the Korean Geotechnical Society
• /
• v.28 no.12
• /
• pp.5-15
• /
• 2012

In this study, a series of CFD (Computational Fluid Dynamics) numerical analyses were performed in order to evaluate the thermal performance of six full-scale closed-loop vertical ground heat exchangers constructed in a test bed located in Wonju. The circulation HDPE pipe, borehole and surrounding ground formation were modeled using FLUENT, a finite-volume method (FVM) program, for analyzing the heat transfer process of the system. Two user-defined functions (UDFs) accounting for the difference in the temperatures of the circulating inflow and outflow fluid and the variation of the surrounding ground temperature with depth were adopted in the FLUENT model. The relevant thermal properties of materials measured in laboratory were used in the numerical analyses to compare the thermal efficiency of various types of the heat exchangers installed in the test bed. The simulation results provide a verification for the in-situ thermal response test (TRT) data. The CFD numerical back-analysis with the ground thermal conductivity of 4 W/mK yielded better agreement with the in-situ thermal response tests than with the ground thermal conductivity of 3 W/mK.

• Journal of the Korean Geotechnical Society
• /
• v.23 no.6
• /
• pp.97-102
• /
• 2007

Though the stability analysis of soil slopes widely employs the limit equilibrium method, the study on the jointed rock slopes must consider the direction of joint and the characteristics of Joint at the same time. This study analyzes the result of the change in the factors which show the characteristics of discontinuity and the shape factor of rock slopes, and so on, in an attempt to validate the propriety as to the interpretation of jointed rock slope stability which uses the general finite element program. First, the difference depending on the flow rules was compared, and the factor effect study was conducted. The selected independent variables included the direction of joint which displays the mechanical characteristics of discontinuity, adhesive cohesion, friction angle, the inclination and height of rock slope which reveal the shape of slope and surcharge load. And the horizontal displacement was numerically interpreted at the 1/3 point below the slope, a dependent variable, to compare the relative degree of factor effects. The findings of study on factor effects led to the validation that the result of horizontal displacement for each factor satisfied various engineering characteristics, making it possible to be applied to stability interpretation of jointed rock slope. A modelling is possible, which considers the application of the result of real geotechnical surveys & laboratory studies and the non-linear characteristics when designing the rock slope. In addition, the stress change which may result from the natural disaster, such as precipitation, and the construction, can be expressed. Furthermore, as the complicated rock condition and the ground supporting effect can be considered through FEM, it is considered to be very useful in making an engineering decision on the cut-slope, reinforcement and so on.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.1C
• /
• pp.53-62
• /
• 2008

A constitutive model, which can simulate the effect of principal stress rotation associated with direct simple shear test, is proposed in this study. The model is based on two mobilized planes. The plastic strains occur from the two mobilized planes, and depend on stress state, and they are added. The first plane is a plane of maximum shear stress, which rotates about the horizontal axis, and the second plane is a horizontal plane which is spatially fixed. The second plane is used to consider the effect of principal stress rotation on simple shear tests under different stress states. The soil skeleton behavior observed in drained simple shear tests is captured in the model. This constitutive model is incorporated into the dynamic coupled stress-flow finite difference program FLAC. The model is first calibrated with drained simple shear tests on loose Fraser River sand. The measured shear stress and volume change are partially induced by principal stress rotation and compared with model calculations. The model is verified by comparing predicted and measured settlements due to rigid footing resting on loose sands. Settlements predicted by the proposed model were very similar to measured settlements. Mohr-Coulomb model can not consider the effect of principal stress rotation and its prediction was only 20% of measured settlements.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.26 no.5D
• /
• pp.797-804
• /
• 2006

A development of regression model for asphalt concrete pavements using Falling Weight Deflectometer deflections is presented in this paper. A backcalculation program based on layered elastic theory was used to generate the synthetic modulus database, which was used to generate 95% confidence intervals of modulus in each layer. Using deflection basins of FWD data used in developing this procedure were collected from Pavement Management System in flexible pavements. Assumptions of back-calculation are that one is 3 layered flexible pavement structure and another is depth to bedrock is finite. It is found that difference of between 95% confidence intervals and modulus ranges of other papers does not exist. So, the data of 95% confidence intervals in each layer was used to develop multiple regression models. Multiple regression equations of each layer were established by SPSS, package of Statics analysis. These models were proved by regression diagnostics, which include case analysis, multi-collinearity analysis, influence diagnostics and analysis of variance. And these models have higher degree of coefficient of determination than 0.75. So this models were applied to predict modulus of domestic asphalt concrete pavement at FWD field test.

• Journal of the Korean Society of Hazard Mitigation
• /
• v.3 no.3 s.10
• /
• pp.151-163
• /
• 2003

In this study, the algorithm of groundwater flow process was established for koreanized groundwater program development dealing with the geographic and geologic conditions of the aquifer have dynamic behaviour in groundwater flow system. All the input data settings of the 3-DFM model which is developed in this study are organized in Korean, and the model contains help function for each input data. Thus, it is designed to get detailed information about each input parameter when the mouse pointer is placed on the corresponding input parameter. This model also is designed to easily specify the geologic boundary condition for each stratum or initial head data in the work sheet. In addition, this model is designed to display boxes for input parameter writing for each analysis condition so that the setting for each parameter is not so complicated as existing MODFLOW is when steady and unsteady flow analysis are performed as well as the analysis for the characteristics of each stratum. Descriptions for input data are displayed on the right side of the window while the analysis results are displayed on the left side as well as the TXT file for this results is available to see. The model developed in this study is a numerical model using finite differential method, and the applicability of the model was examined by comparing and analyzing observed and simulated groundwater heads computed by the application of real recharge amount and the estimation of parameters. The 3-DFM model is applied in this study to Sehwa-ri, and Songdang-ri area, Jeju, Korea for analysis of groundwater flow system according to pumping, and obtained the results that the observed and computed groundwater head were almost in accordance with each other showing the range of 0.03 - 0.07 error percent. It is analyzed that the groundwater flow distributed evenly from Nopen-orum and Munseogi-orum to Wolang-bong, Yongnuni-orum, and Songja-bong through the computation of equipotentials and velocity vector using the analysis result of simulation which was performed before the pumping started in the study area. These analysis results show the accordance with MODFLOW's.

• Geophysics and Geophysical Exploration
• /
• v.1 no.2
• /
• pp.127-135
• /
• 1998

The Born approximation is widely used for solving the complex scattering problems in electromagnetics. Approximating total internal electric field by the background field is reasonable for small material contrasts as long as scatterer is not too large and the frequency is not too high. However in many geophysical applications, moderate and high conductivity contrasts cause both real and imaginary part of internal electric field to differ greatly from background. In the extended Born approximation, which can improve the accuracy of Born approximation dramatically, the total electric field in the integral over the scattering volume is approximated by the background electric field projected to a depolarization tensor. The finite difference and elements methods are usually used in EM scattering problems with a 2D model and a 3D source, due to their capability for simulating complex subsurface conductivity distributions. The price paid for a 3D source is that many wavenumber domain solutions and their inverse Fourier transform must be computed. In these differential equation methods, all the area including homogeneous region should be discretized, which increases the number of nodes and matrix size. Therefore, the differential equation methods need a lot of computing time and large memory. In this study, EM modeling program for a 2D model and a 3D source is developed, which is based on the extended Born approximation. The solution is very fast and stable. Using the program, crosshole EM responses with a vertical magnetic dipole source are obtained and the results are compared with those of 3D integral equation solutions. The agreement between the integral equation solution and extended Born approximation is remarkable within the entire frequency range, but degrades with the increase of conductivity contrast between anomalous body and background medium. The extended Born approximation is accurate in the case conductivity contrast is lower than 1:10. Therefore, the location and conductivity of the anomalous body can be estimated effectively by the extended Born approximation although the quantitative estimate of conductivity is difficult for the case conductivity contrast is too high.

• The Journal of the Petrological Society of Korea
• /
• v.6 no.2
• /
• pp.133-149
• /
• 1997

In the metapelites around the Ogcheon granite, the metamorphic grade increases from the biotite zone through the andalusite zone to the sillimanite zone towards the intrusion contact. In the metabasites around the Boeun granite, the metamorphic grade increases from transitional zone between the greenchist and amphibolite facies through the amphibolite facies to the upper amphibolite facies towards the intrusion contact. In the Doiri area locating near the intrusion contact of the Boeun granite, sillimanite- and andalusite-bearing metapelites are found with in 500 m away from the contact. The evidence described above indicates that the Ogcheon and Boeun granites caused low-P/T type contact metamorphism to the country rocks. The P-T condition of contact metamorphism due to the intrusion of the Ogcheon granite is $540{\pm}40^{circ}C, 2.8{\pm}0.9$ kb. The temperature condition of contact metamorphism due to the intrusion of the Boeun granite is $698{\pm}28^{\circ}C$. The wide compositional range of amphibole and plagioclase in the metabasites around the Boeun granite is due to the immisibility gab of amphibole and plagioclase and unstable relict composition resulted from an incomplete metamorphic reaction. The compositional range of stable amphibole and plagioclase decreases as a metamorphic grade increases due to a close of immiscibility gab. The thermal effect of contact metamorphism due to the intrusion of the Ogcheon and Boeun granites, are calculated using the CONTACT2 program based on a two dimensional finite difference method. In order to estimate the thermal effect of an introduced pluton, a circle with 10 km diameter and a triangle with 20 km side are used for the intrusion geometries of the Ogcheon granite and the Boeun granite, respectively. The results from the field and modeling studies suggest that the intrusion temperatures of the Ogcheon granite close to $800^{\circ}C$ and the intrusion temperature of the Boeun granite is higher than $1000^{\circ}C$. However, the intrusion temperatures can be lower than the suggested temperature, if the geothermal gradient prior to the intrusion of the Ogcheon and Boeun granites was higher than the normal continental grothermal gradient.