• Tunnel and Underground Space
• /
• v.32 no.6
• /
• pp.598-610
• /
• 2022

Avoiding a fault zone would be a best practice for safety in underground construction, which is only sometimes possible because of many restrictions and other field conditions. For instance, there is an ongoing conception of Korea-Japan subsea tunnels that inevitably cross a massive fault system in the Korea Strait. Therefore it was deemed necessary to find an efficient way of predicting the likely behaviour of underground structures under fault slip. This paper presents the findings from simple numerical analysis for investigating the stress induced at a normal fault with a dip of 45 degrees. We used a boundary element software that assumed constant displacement discontinuity, which allowed the displacement to be estimated separately at both the fault's hangingwall and footwall sides. The results suggested that a principal stress rotation of 45 degrees occurred at the edges of the fault during the slip, which was in agreement with the phenomenon for fault plane suggested in the body of literature. A simple numerical procedure presented in this paper could be adopted to investigate other fault-related issues associated with underground structure construction.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.6 no.2
• /
• pp.139-149
• /
• 1994

In this paper, a finite element technique is applied to the prediction of the wave resonance phenomena in harbors. The mild-slope equation is used with a partial reflection boundary condition introduced to model the energy dissipating effects on the solid boundary. For an efficient modeling of the radiation condition at infinity, a new infinite element is developed. The shape function of the infinite element is derived from the asymptotic behavior of the first kind of the Hankel's function in the analytical boundary series solutions. For the computational efficiency, the system matrices of the element are constructed by performing the relevant integrations in the infinite direction analytically. Comparisons with the results from experiments and other solution methods show that the present model gives fairly good results. Numerical experiments are also carried out to determine the proper distance to the infinite elements from the mouth of the halter, which directly affect the accuracy and efficiency of the solution.

• The Korean Journal of Petroleum Geology
• /
• v.5 no.1_2 s.6
• /
• pp.9-15
• /
• 1997

The study area is located in the Central Bransfield Basin, Antarctica. To analyze the morphology of seafloor, structure of basement, and seismic stratigraphy of the sedimentary layers, we have acquired, processed, and interpreted the multi-channel seismic data. The northwest-southeastern back-arc extension dramatically changes seafloor morphology, volcanic and fault distribution, and basin structure along the spreading ridges. The northern continental shelf shows a narrow, steep topography. In contrast, the continental shelf or slope in the south, which is connected to the Antarctic Peninsula, has a gentle gradient. Volcanic activities resulted in the formation of large volcanos and basement highs near the spreading center, and small-scale volcanic diapirs on the shelf. A very long, continuous normal fault characterizes the northern shelf, whereas several basinward synthetic faults probably detach into the master fault in the south. Four transfer faults, the northwest-southeastern deep-parallel structures, controlled the complex distributions of the volcanos, normal faults, depocenters, and possibly hydrocarbon provinces in the study area. They have also deformed the basement structure and depositional pattern. Even though the Bransfield Basin was believed to be formed in the Late Cenozoic (about 4 Ma), the hydrocarbon potential may be very high due to thick sediment accumulation, high organic contents, high heat flow resulted from the active tectonics, and adequate traps.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.3
• /
• pp.185-192
• /
• 1993

A numerical model has been used for the prediction of wave agitations in a harbor which are induced by the intrusion and transformation of incident waves. Based on linear wave theory a mild-slope equation has been used. A partial absorbing boundary condition has been used on solid boundary. Functional has been derived following Chen and Mei(l974)'s technique based on Hybrid Element Method which uses finite discretisation in the inner region and analytical solution of Helmholtz equation in the outer region. Final simultaneous equation has been solved using the Gaussian Elimination Method. Helmholtz natural period and second peak period of seiche in Donghae Harbor coincide very well with the results from numerical calculation. Computed amplification factors show good agreement, especially when the reflection coefficient on solid boundary is 0.99, with those of measurements.

• Journal of the Korean GEO-environmental Society
• /
• v.19 no.12
• /
• pp.25-32
• /
• 2018

This paper investigates the effect of rock discontinuities on a shear stress wave that is induced by earthquake or blasting and provides the result of numerical parametric studies. The numerical tests of different conditions of rock and discontinuity have been carried out after confirming that the numerical approach is valid throughout a verification analysis from which the test results were compared with a theoretical solution. In-situ stress condition was considered as a rock condition and internal friction angle and cohesive value, which are the shear strength parameters, were considered as discontinuities condition. The joint inclination angle was also taken into account as a parameter. With the various conditions of different parameters, the test results showed that a shear stress wave propagating through a mass is highly influenced by the shear strength of discontinuities and the condition of joint inclination angle as well as in-situ stress. The study results indicate that when earthquake or blasting-induced dynamic loading propagates through a jointed rock mass or a stratified soil ground the effect of in-situ stress and discontinuities including a stratum boundary should be taken into account when evaluating the dynamic effect on nearby facilities and structures.

• The Journal of Engineering Geology
• /
• v.10 no.3
• /
• pp.205-216
• /
• 2000

An unstable natural slope located at Icheon-Ri, Dasa-Eup, Dalseong-Kun which has been severely cracked during heavy rain fall in September, 1998, was investigated and analysed for remedial measures. Various geological data were obtained by performing a surface geological survey, a ground penetrating radar(GPR) exploration, four geotechnical borings with field and laboratory tests, and the precipitation records were also reviewed. Based on the these data, the probable sliding plane was determined as a bedding plane with dip angle of 20 to 25 degrees between sandstone and shale layers at depth of 6 to 8 meters. By performing back analyses with the limit equilibrium method, the friction angle and cohesion of the sliding plane were determined as 18$^{\circ}$ and 0, respectively. Based on these results, a series of stability analyses such as stereo-graphic projection method, limiting equilibrium method and numerical method were carried out. The remedial measures for the unstable slope were selected based on the results of the limit equilibrium analyses by considering various conditions expected at the site. Ground grouting or anchoring method with proper surface drainage systems could be employed for the long term safety of the slope.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.1
• /
• pp.169-178
• /
• 1991

A plane-strain finite element analysis of sheet metal forming is carried out by using the rigid-plastic FEM based on the membrane theory. The sheet material is assumed to possess normal anisotropy and to obey Hill's new yield criterion and its associated flow rule. A formulation of initial guess generation for the displacement field is derived by using the nonlinear elastic FEM. A method of contact treatment is proposed in which the skew boundary condition for arbitrarily shaped tools is successively used during iteration. In order to verify the validity of the developed method, plane-strain drawing with tools in analytic expression and with arbitrarily shaped tools is analyzed and compared with the published results. The comparison shows that the present method can be effectively used in the analysis of plane-strain sheet metal forming and thus provides the basis of approximate sectional analysis of panel-like sheet forming.

• Water for future
• /
• v.29 no.2
• /
• pp.143-151
• /
• 1996

The dynamic responses in the subsurface outflow, surface saturation area, soil moisture storage are established by numerical experiments with Richards equation. In addition to this, the dynamical relationship between surface saturation area and subsurface outflow, and between surface saturation area and soil moisture storage are also determined by varying the hillslope shape, soil type, and boundary condition. The simulation results indicate that the dynamical relationships between surface saturation area and subsurface outflow, and between surface saturation area and soil moisture storage are approximated by the steadystate relationships. And the dynamic responses of subsurface outflow and surface saturation area are characterized by the double peaks although the rainfall pattern is asimple pulse input. As a result of numerical simulation, the physical mechanism for the occurrence of the double peaks is explained using the concept of variable source area.

• Composites Research
• /
• v.12 no.2
• /
• pp.82-90
• /
• 1999

The mechanical properties of generally orthotropic materials are conventionally measured by performing off-axis tensile and flexure tests. However, the inevitable coupling between tension and shear in case of tensile test or bending and twisting in flexure test case induces nonuniform displacement and stress fields. Consequential stress concentration along the boundary of specimens would result in inaccurate modulus and underestimated strength. This paper proposes the variation of specimen geometry in terms of appropriate obliquity of loaded boundary. For the purpose, classical lamination theory is transformed into skewed coordinate, and characteristic equations for both of unidirectional and laminated composite specimens are formulated. Finite element analysis is employed to show the validity of the skewedness in tensile and bending test specimens.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2021.06a
• /
• pp.105-105
• /
• 2021

본 연구에서는 2차원 비선형 천수모형과 수심평균된 스칼라 이송모형을 해석하는 수치모형에 대해 기술하였다. 수치모형의 정확성을 보장함과 동시에 안정성을 높이기 위해 유한체적법, 플럭스 재구성 및 minmod 제한자를 사용하였다. 비선형 천수방정식의 이송항과 바닥 경사항은 계산된 수심의 양수 보존과 흐름의 정상 상태를 보장하기 위한 second order well-balanced positivity preserving central-upwind method를 이용하여 수치적으로 이산화되었다. 마찬가지로, 이송-확산 방정식 내 이송항은 동일한 2차 풍상차분법을 통해 수치적으로 풀이하였다. 격자점 경계면에서의 불연속으로 인한 수치진동을 방지하기 위해 이송항의 계산에 포함된 보존항의 차이로 인해 발생하는 스칼라의 수치확산을 최소화하기 위해 무차원의 비소산함수를 도입하였다. 또한, 확산항은 유한차분법을 이용하여 이산화하였다. 제안된 수치모형은 시간미분항의 계산을 위해 오일러 기법을 적용하여 계산된 수심 및 스칼라의 양수 보존여부와 함께 정지된 흐름의 정상 상태의 보존여부를 확인하였다. 제안된 수치모형의 해석 정확성을 평가하기 위해 1, 2차원 공간 내 다양한 흐름 조건에서의 해석해를 이용한 3개의 벤치마크 테스트를 수행하였다. 평균 제곱근 오차(Root Mean Squared Error, RMSE)를 산정하여 수치모형의 성능을 정량적으로 평가하였으며, 비소산함수를 적용함에 따라 스칼라의 수치확산이 감소하게 되었음을 확인하였다. 또한, 세 차례의 벤치마크 테스트 결과는 공통적으로 수치모형에 의해 계산된 결과값이 비소산함수를 고려함에 따라 해석해와 잘 일치함을 확인하였다.