• Geophysics and Geophysical Exploration
• /
• v.21 no.3
• /
• pp.150-161
• /
• 2018

Recently, unmanned aircraft EM (electromagnetic) survey based on ICT (Information and Communication Technology) has been widely utilized because of the efficiency in regional survey. We performed the theoretical study on the unmanned airship EM system developed by KIGAM (Korea Institute of Geoscience and Mineral resources) as part of the practical application of unmanned aircraft EM survey. Since this system has different configurations of transmitting and receiving loops compared to the conventional aircraft EM systems, a new technique is required for the appropriate interpretation of measured responses. Therefore, we proposed a method to calculate the EM field for the arbitrary shaped transmitter and verified its validity through the comparison with analytic solution for circular loop. In addition, to simulate the magnetic responses by three-dimensionally (3D) distributed anomalies, we have adapted our algorithm to 3D frequency-domain EM modeling algorithm based on the edge-FEM (finite element method). Though the analysis on magnetic field responses from a subsurface anomaly, it was found that the response decreases as the depth of the anomaly increases or the flight altitude increases. Also, it was confirmed that the response became smaller as the resistivity of the anomaly increases. However, a nonlinear trend of the out-of-phase component is shown depending on the depth of the anomaly and the used frequency, that makes it difficult to apply simple analysis based on the mapping of the magnitude of the responses and can cause the non-uniqueness problem in calculating the apparent resistivity. Thus, it is a prerequisite to analyze the appropriate frequency band and flight altitude considering the purpose of the survey and the site conditions when conducting a survey using the unmanned aircraft EM system.

• Journal of the Korean Society for Railway
• /
• v.20 no.5
• /
• pp.668-682
• /
• 2017

Alignments of railways recently constructed in Korea have been straightened due to the advent of high-speed rail, which means increasing the numbers of tunnels and bridges. Overbreak during tunnel construction may be unavoidable, and is very influential on overall stability. Over-excavation in tunneling is also one of the most important factors in construction costs. Overbreak problems around crown areas have decreased with improvements of excavation methods, but overbreak problems around bottom areas have not decreased because those areas are not very influential on tunnel stability compared with crown areas. The filling costs of 10 cm thickness of overbreak at the bottom of a tunnel are covered under construction costs by Korea Railway Authority regulations, but filling costs for more than the covered thickness are considered losses of construction cost. The filling material for overbreak bottoms of tunnels should be concrete, but concrete and mixed granular materials with fractured rock are also used for some sites. Tunnels in which granular materials with fractured rock are used may have a discontinuous section under the concrete slab track. The discontinuous section influences the propagation of waves generated from train operation. When the bottom of a tunnel is filled with only concrete material, the bottom of the tunnel can be considered as a continuous section, in which the waves generated from a train may propagate without reflection waves. However, a discontinuous section filled with mixed granular materials may reflect waves, which can cause resonance of vibration. The filled materials and vibration propagation characteristics are studied in this research. Tunnel bottom filling materials that have ratios of granular material to concrete of 5.0 %, 11.5 %, and 18.0 % are investigated. Samples were made and tested to determine their material properties. Static numerical analyses were performed using the FEM program under train operation load; test results were found to satisfy the stability requirements. However, dynamic analysis results show that some mixed ratios may generate resonance vibration from train operation at certain speeds.

• The Journal of Engineering Geology
• /
• v.29 no.2
• /
• pp.85-97
• /
• 2019

It can be observed that steep slopes ($65^{\circ}$ to $80^{\circ}$) consist of rock masses were kept stable for a long time. In rock-mass slopes with similar ground condition, steeper slopes than 1 : 0.5 ($63^{\circ}$) may be applied if the discontinuities of rock-mass slope are distributed in a direction favorable to the stability of the slope. In making a decision the angle of the slope, if the preliminary rock mass conditions applicable to steep slope are quantitatively setup, they may be used as guidance in design practice. In this study, the above rock mass was defined as a good continuum rock mass and the quantitative setup criterion range was proposed using RMR, SMR and GSI classifications for the purpose of providing engineering standard for good continuum rock mass conditions. The methods of study are as follows. The stable slope at steep slopes ($65^{\circ}$ to $80^{\circ}$) for each rock type was selected as the study area, and RMR, SMR and GSI were classified to reflect the face mapping results. The results were reviewed by applying the calculated shear strength to the stable analysis of the current state of rock mass slope using the Hoek-Brown failure criterion. It is intended to verify the validity of the preliminary criterion as a rock mass condition that remains stable on a steep slope. Based on the analysis and review by the above research method, it was analyzed that a good continuum rock mass slope can be set to Basic RMR ${\geq}50$ (45 in sedimentary rock), GSI and SMR ${\geq}45$. The safety factor of the LEM is between Fs = 14.08 and 67.50 (average 32.9), and the displacement of the FEM is 0.13 to 0.64 mm (average 0.27 mm). This can be seen as a result of quantitative representation and verification of the stability of a good continuum rock mass slope that has been maintained stable for a long period of time with steep slopes ($65^{\circ}$ to $80^{\circ}$). The setup guideline for a good continuum rock mass slope will be able to establish a more detailed setup standard when the data are accumulated, and it is also a further study project. If stable even on steep slopes of 1 : 0.1 to 0.3, the upper limit of steep slopes is 1 : 0.3 with reference to the overseas design standards and report, thus giving the benefit of ensuring economic and eco-friendlyness. Also, the development of excavation technology and plantation technology and various eco-friendly slope design techniques will help overcome psychological anxiety and rapid weathering and relaxation due to steep slope construction.

• 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.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.23 no.1
• /
• pp.43-49
• /
• 2011

For the improvement of water quality in a harbor, several studies have been carried out on SEB (Seawater Exchange Breakwater) in recent years, but a problem has been shown whereby the water on the inside area far from the SEB cannot be easily exchanged. In order to solve the problem of the SEB, the Manifold channel, a new concept of the SEB, is introduced in this paper. By using the manifold channel, it is possible to exchange the water of the inside area for seawater from the outside. Here, to assess the outflow gates of the manifold channel governing flow behavior, a virtual manifold channel controlled the location, width and direction of outflow gates applied to the Jumunjin fishery port, where the SEB has been established. In addition, the desirable flow pattern of the port by utilizing the two layer current model is identified, and five general cases of the manifold channel are described in this paper. The model is verified by comparing with observation of the SEB model, and the results are in general agreement. From the results of the manifold channel, in the case of the Jumunjin fishery port, the small circulation of counter clockwise is necessary for the water exchange on the inside area, but it should be controlled by the outflow gates for other areas. Using the two layer current model, the desirable flow pattern of the port can be predicted, and the water exchange for the upper and lower layer can be examined. For the practical use of the manifold channel, further studies on the manifold channel will be necessary, and it may then be used broadly for the design of breakwater in the future.

• Journal of the Korean Society of Groundwater Environment
• /
• v.7 no.2
• /
• pp.89-96
• /
• 2000

Slug tests are the most widely used field method for quantification of hydraulic conductivity of porous media. Well recovery is affected by well casing, borehole radii, screened length, hydraulic conductivity, and specific storage of porous media. In this study, a new slug tests model was developed through finite element approximation and the validity and usefulness of the model were tested in various ways. Water level fluctuation in a well under slug test and cons-equent groundwater flow in the surrounding porous medium were appropriately coupled through estimation of well-flux using an iteration technique. Numerical accuracy of the model was verified using the Cooper et al. (1967) solution. The model has advantages in simulations for monitored slug tests, partial penetration, and inclusion of storage factor. Volume coverage of slug tests is significantly affected by storage factor. Magnitude and speed of propagation of head changes from a well increases as storage factor becomes low. It will be beneficial to use type curves of monitored head transients in the surrounding porous formation for estimation of specific storage. As the vertical component of groundwater flow is enhanced, the influence of storage factor on well recovery decreases. For a radial-vertical flow around a partially penetrated well, deviations between hydraulic estimates by various methods and data selection of recovery curve are negligible on practical purposes, whereas the deviations are somewhat significant for a radial flow.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.22 no.6
• /
• pp.615-624
• /
• 2016

In this study, two rectangular barges in close proximity were simulated to analyze the characteristics of motion responses due to hydrodynamic interactions. Using a numerical solution from DNV-GL SESAM, coupled stiffness matrix terms for these same FEM models were added to the multiple body modes in the surge direction. Potential theory was used to calculate the first order radiation and diffraction effects on the simulated barge models. In the results, the sheltering effect of the barges was not shown at 1.3 rad/s with hull separation of 20 m in transverse waves. The separation effect between the barges was more clear with longitudinal waves and a shallow water depth. However, sway forces were influenced by hull separation with transverse waves. The peaks for sway and heave motion and sway force occurred at higher frequencies as hull separation narrowed with longitudinal and transverse waves. Given a depth of 10 m, the sway motion on the lee side of a coupled barge made a significant difference in the range of 0.2-0.8 rad/s with transverse and oblique waves. Also, the peaks for sway force were situated at lower frequencies, even when incident waves changed.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.18 no.3
• /
• pp.209-214
• /
• 2017

Most of the present potential transformers of train vehicles are of the oil-type filled with insulating oil and are susceptible to problems such as explosion due to the increase in the internal pressure during train operation and poor reliability near the end of their life cycle. As a solution to this problem, it is necessary to develop a molded dry-type potential transformer with excellent pressure-resistance performance using insulating resin. In order to localize the product, the Korea Railroad Research Institute has been developing a molded dry-type potential transformer. As part of this research, it is necessary to analyze the vibration characteristics of the developed product and to check the transformer performance in a vibration environment. In this study, a resonance test and simulated long-term life test of the developed product were conducted according to the KS R 9144 and IEC 61373 standards, respectively, which are vibration test methods for railway vehicle parts. Their natural frequencies were analyzed by comparing the results of the numerical modal analysis and resonance test, in order to confirm their adherence to the standards. Also, the performance test after the simulated long-term life test confirmed that the operation of the developed transformer was not problematic even in a long-time vibration environment.

• Geophysics and Geophysical Exploration
• /
• v.6 no.1
• /
• pp.13-22
• /
• 2003

TSP (Tunnel Seismic Profiling) survey is a technique for imaging and characterizing geological structures ahead of a tunnel face. The seismic modeling algorithm and the synthetic data could be helpful for TSP surveys. However, there is few algorithm to describe the propagation of the elastic waves around the tunnel. In this study, existing 2-dimensional seismic modeling algorithm using finite element method was modified to make a suitable algorithm for TSP modeling. Using this algorithm, TSP modeling was practiced in some models. And the synthetic data was analyzed to examine the propagation characteristics of the elastic waves. First of all, the modeling for the homogeneous tunnel model was practiced to examine the propagation characteristics of the direct waves in the vicinity of the tunnel. And the algorithm was applied to some models having reflector which is perpendicular or parallel to the excavation direction. From these, the propagation characteristics of the reflected waves were examined. Furthermore, two source-receiver arrays were used in respective models to investigate the properties of the two arrays. These modeling algorithm and synthetic data could be helpful in interpreting TSP survey data, developing inversion algorithm and designing new source-receiver arrays.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.24 no.5 s.176
• /
• pp.1075-1083
• /
• 2000

The prediction of the inelastic behavior of the structure is an essential part of reliability assessment procedure, because most of the failures are induced by the inelastic deformation, such as creep and plastic deformation. During decades, there has been much progress in understanding of the inelastic behavior of the materials and a lot of inelastic constitutive equations have been developed. These equations consist of the definition of inelastic strain and the evolution of the state variables introduced to quantify the irreversible processes occurred in the material. With respect to the definition of the inelastic strain, the inelastic constitutive models can be categorized into elastoplastic model, unified viscoplastic model and separated viscoplastic model and the different integration methods have been applied to each category. In the present investigation, the generalized integration method applicable for various types of constitutive equations is developed and implemented into ABAQUS by means of UMAT subroutine. The solution of the non-linear system of algebraic equations arising from time discretization with the generalized midpoint rule is determined using line-search technique in combination with Newton method. The strategy to control the time increment for the improvement of the accuracy of the numerical integration is proposed. Several numerical examples are considered to demonstrate the efficiency and applicability of the present method. The prediction of the inelastic behavior of the structure is an essential part of reliability assessment procedure, because most of the failures are induced by the inelastic deformation, such as creep and plastic deformation. During decades, there has been much progress in understanding of the inelastic behavior of the materials and a lot of inelastic constitutive equations have been developed. These equations consist of the definition of inelastic strain and the evolution of the state variables introduced to quantify the irreversible processes occurred in the material. With respect to the definition of the inelastic strain, the inelastic constitutive models can be categorized into elastoplastic model, unified viscoplastic model and separated viscoplastic model and the different integration methods have been applied to each category. In the present investigation, the generalized integration method applicable for various types of constitutive equations is developed and implemented into ABAQUS by means of UMAT subroutine. The solution of the non-linear system of algebraic equations arising from time discretization with the generalized midpoint rule is determined using line-search technique in combination with Newton method. The strategy to control the time increment for the improvement of the accuracy of the numerical integration is proposed. Several numerical examples are considered to demonstrate the efficiency and applicability of the present method.