• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.19 no.10
• /
• pp.24-30
• /
• 2018

As the demand for underground space increases, many researchers have been studying the load reduction method using high compressible materials to solve for the stability problem of the overhead load and for the increase of the earth pressure which decreases the function of the underground structure. This paper determines the optimum soft zone and the effect of the using EPS Geofoam as a load reduction material to arch structures. A finite element analysis program, ABAQUS, is used to analyze the soil-structure interaction and the behavior of buried arch structures considering different four EPS Geofoam forms to confirm the most conservative shape. The optimum cross-sectional shape was determined by comparing the results of earth pressure reduction rate in accordance with the change of span-rise ratio and span length of the arch structure. It was confirmed that the earth pressure generated in the arch structure using the optimal soft zone selected by the numerical analysis was reduced by an average of 78%. In this study, the effect of EPS Geofoam on soil pressure reduction and its applicability to underground arch structures will provide an economical and conservative way to design underground structures and will help to increase the usability of deep underground space.

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

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.19 no.3
• /
• pp.389-407
• /
• 2017

In the current work, a series of three-dimensional finite element analysis was carried out to understand the behaviour of pile when the tunnel passes through the lower part of a single pile or group piles. At the current study, the numerical analysis analysed the results regarding the ground reinforcement condition between the tunnel and pile foundation. In the numerical modelling, several key issues, such as the pile settlements, the axial pile forces, the shear stresses and the total displacements near the tunnel have been thoroughly analysed. The pile head settlements of the single pile with the maximum level of reinforcement decreased by about 16% compared to the pile without ground reinforcement. Furthermore, the maximum axial force of the single pile with the maximum level of ground reinforcement experienced a 30% reduction compared to the pile without reinforcement. It has been found that the angle of ground reinforcement in the transverse direction affects the pile behaviour more so than the length of the ground reinforcement in the longitudinal direction. On the other hand, in the case of the pile group with the reinforced pile cap, the ground displacement near the pile tip appears to be similar to the corresponding ground displacement without reinforcement. However, it was found that the pile cap near the pile head greatly restrained the pile head movement and hence the axial pile force increased by about 2.5 times near the pile top compared to the piles in other analysis conditions. The behaviour of the single pile and group piles, depending on the amount of ground reinforcement, has been extensively examined and analysed by considering the key features in great details.

• Journal of Korean Society of Environmental Engineers
• /
• v.30 no.11
• /
• pp.1146-1153
• /
• 2008

Numerical simulation was carried out to study the trichloroethylene (TCE) degradation by permeable reactive barrier (PRB), and revealed the effect of concentration of TCE, iron medium mass, and concentration of iron-reducing bacteria (IRB). Newly developed model was based on axial dispersion reactor model with chemical and biological reaction terms and was implemented using MATLAB ver R2006A for the numerical solutions of dispersion, convection, and reactions over column length and elapsed time. The reaction terms include reactions of TCE degradation by zero-valent iron (ZVI, Fe$^0$) and ferrous iron (Fe$^{2+}$). TCE concentration in the column inlet was maintained as 10 mg/L. Equation for Fe$^0$ degradation includes only TCE reaction term, while one for Fe$^{2+}$ has chemical and biological reaction terms with TCE and IRB, respectively. Two coupled equations eventually modeled the change of TCE concentration in a column. At Fe$^0$ column, TCE degradation rate was found to be more than 99% from 60 hours to 235 hours, and declined to less than 1% in 1,365 hours. At the Fe$^{2+}$ and IRB mixed column, TCE degradation rate was equilibrated at 85.3% after 210 hours and kept it constant. These results imply that the ferrous iron produced by IRB has lowered the TCE degradation efficiency than ZVI but it can have higher longevity.http://kci.go.kr/kciportal/ci/contents/ciConnReprerSearchPopup.kci#

• Journal of KIISE:Software and Applications
• /
• v.30 no.3_4
• /
• pp.212-218
• /
• 2003

State invariant is a property that holds in every reachable state. It can be used not only in understanding and analyzing complex software systems, but it can also be used for system verifications such as checking safety, liveness, and consistency. For these reasons, there are many vital researches for deriving state invariant from finite state machine models. In previous works every reachable state is to be considered to generate state invariant. Thus it is likely to be too complex for the user to understand. This paper seeks to answer the question `how to simplify state invariant\ulcorner`. Since the complexity of state invariant is strongly dependent upon the size of states to be considered, so the smaller the set of states to be considered is, the shorter the length of state invariant is. For doing so, we let the user focus on some interested scopes rather than a whole state space in a model. Computation Tree Logic(CTL) is used to specify scopes in which he/she is interested. Given a scope in CTL, mixed reachability analysis is used to find out a set of states inside it. Obviously, a set of states calculated in this way is a subset of every reachable state. Therefore, we give a weaker, but comprehensible, state invariant.

• Journal of the Korean Association of Oral and Maxillofacial Surgeons
• /
• v.38 no.6
• /
• pp.337-342
• /
• 2012

Objectives: Implants connect the internal body to its external structure, and is mainly supported by alveolar bone. Stable osseointegration is therefore required when implants are inserted into bone to retain structural integrity. In this paper, we present an implant with a "wing" design on its area. This type of implant improved stress distribution patterns and promoted changes in bone remodeling. Materials and Methods: Finite element analysis was performed on two types of implants. One implant was designed to have wings on its cervical area, and the other was a general root form type. On each implant, tensile and compressive forces ($30N/m^2$, $35N/m^2$, $40N/m^2$, and $45N/m^2$) were loaded in the vertical direction. Stress distribution and displacement were subsequently measured. Results: The maximum stresses measured for the compressive forces of the wing-type implant were $21.5979N/m^2$, $25.1974N/m^2$, $29.7971N/m^2$, and $32.3967N/m^2$ when $30N/m^2$, $35N/m^2$, $40N/m^2$, and $45N/m^2$ were loaded, respectively. The maximum stresses measured for the root form type were $23.0442N/m^2$, $26.9950N/m^2$, $30.7257N/m^2$, and $34.5584N/m^2$ when $30N/m^2$, $35N/m^2$, $40N/m^2$, and $45N/m^2$ were loaded, respectively. Thus, the maximum stresses measured for the tensile force of the root form implant were significantly higher (about three times greater) than the wing-type implant. The displacement of each implant showed no significant difference. Modifying the design of cervical implants improves the strength of bone structure surrounding these implants. In this study, we used the wing-type cervical design to reduce both compressive and tensile distribution forces loaded onto the surrounding structures. In future studies, we will optimize implant length and placement to improve results. Conclusion: 1. Changing the cervical design of implants improves stress distribution to the surrounding bone. 2. The wing-type implant yielded better results, in terms of stress distribution, than the former root-type implant.

• Journal of the Korean Geotechnical Society
• /
• v.15 no.5
• /
• pp.217-228
• /
• 1999

The three methods of analysis (i) REV(representative elemental volume), (ii) non-REV and (iii) joint network analysis are introduced in this paper to analyze the groundwater flow in jointed rock mass and the inflow into underground excavations. The results from those methods are compared one another to reveal their characteristics by varying the number of joints and the diameter of the opening. The pre-processor, the so-called sequential analysis, is introduced to predict the equivalent hydraulic conductivity of a jointed rock mass having a number of intersecting joints. Using the finite element mesh, joint map and sequential analysis, the equivalent hydraulic conductivities are calculated for all 445 elements. The hydraulic inhomogeneity and the determination of the representative properties of jointed rock masses are discussed. In the REV analysis where the entire rock mass is homogenized through the representative properties, the inflow is increased regularly and consistently by increasing the joint density, the opening size and the conductivity contrast value. Though the non-REV analysis showed irregular variation of the inflow due to the local inhomogeneity allowed to individual elements, the inflow approached the REV results as the characteristic length increases. The joint network analysis showed the most sensitive reaction to the joint density, the opening size and the presence of the network crossing the opening. The reliability of the network analysis depends on the geometric data of individual joints. In view of the limited field data on joint geometry and possible uncertainty the REV and non-REV methods are considered more practical and rational than the joint network analysis.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.36 no.2
• /
• pp.181-188
• /
• 2012

The present study on nitrogen-diluted non-premixed counterflow flames with finite burner diameters experimentally investigates the important role of the outer edge flame in flame extinction. Flame stability diagrams mapping the flame extinction response of nitrogen-diluted non-premixed counterflow flames to varying global strain rates in terms of the burner diameter, burner gap, and velocity ratio are explored. There exists a critical nitrogen mole fraction beyond which the flame cannot be sustained, and also the curves of the critical nitrogen mole fraction versus the global strain rate have C-shapes in terms of burner diameter, burner gap, and velocity ratio. In flames with sufficiently high strain rates, the curves of the critical nitrogen mole fractions versus global strain rate collapse into one curve, and the flames can have the 1-D flame response of typical diffusion flames. Three flame extinction modes are identified: flame extinctions through the shrinkage of the outer edge flame with and without an oscillation of the outer edge flame prior to the extinction and flame extinction through a flame hole at the flame center. The measured flame surface temperature and a numerical evaluation of the fractional contribution of each term in the energy equation show that the radial conductive heat loss at the flame edge destabilizes the outer edge flame, and the conductive and convection heat addition to the outer edge from the trailing diffusion flame stabilizes the outer edge flame. The radial conductive heat loss at the flame edge is the dominant extinction mechanism acting through the shrinkage of the outer edge flame.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.6
• /
• pp.1299-1308
• /
• 1994

This paper is considered on the dynamic behavior and the dynamic impact coefficient on the cable-stayed bridge under the vehicle load. The method of static analysis, that is, the transfer matrix method is used to get influence values about displacements, section forces of girder and cable forces. Gotten influence values were used as basic data to analyse dynamic behavior. This paper used the transfer matrix method because it is relatively simpler than the finite element method, and calculating speed of computer is very fast and the precision of computation is high. In the process of dynamic analysis, the uncoupled equation of motion is derived from simultaneous equation of the motion of cable-stayed bridge and vehicle travelling by using mode shape, which was borne from system of undamped free vibration. The solution of the uncoupled equation of motion, that is, time history of response of deflections, velocity and acceleration on reference coordinate system, is found by Newmark-${\beta}$ method, a kind of direct integral method. After the time history of dynamic response was gotten, and it was transfered to the time history of dynamic response of cable-stayed bridge by linear transformation of coordinates. As a result of this numerical analysis, effect of dynamic behavior for cable-stayed bridge under the vehicle load has varied depending on parameter of design, that is, the ratio of span, the ratio of main span length, tower height, the flexural rigidity of longitudinal girder, the flexural rigidity of tower, and the cable stiffness, investigated. Very good agreements with the existing solution in the literature are shown for the uncracked plate as well as the cracked plate.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.22 no.1
• /
• pp.82-89
• /
• 2016

In this study, we proposed a pier of pontoon-type, "Very Large Floating Structure" (VLFS), with the length of 500m, breadth of 200 m and height of 2 m in Yeosu domestic port. Since this structure ought to endure wave loads for long periods at sea, it is essential to analyze the wave response characteristics. Direct-method is used to analyze the fluid-structure problem and the coupled motion of equation is used to obtain response results. The structural part is calculated by using finite element method (FEM) and the fluid part is analyzed by using boundary element method (BEM). Dynamic responses caused by the elastic deformation and rigid motion of structure are analyzed by numerical calculation. To investigate response characteristics of the pier in regular waves, several factors such as the wavelength, water depth, wave direction and flexural rigidity of structure are considered. As a result, wave response of pier changed at the point of $L/{\lambda}$ 1.5 and represented the torsional phenomenon according to the various incident waves. And the responses showed increasing tendency as the water depths increase at the incident point in case of $L/{\lambda}=8.0$ and peak point of vertical displacement amplitude moved from side to side as the flexural rigidity of structure changes.