• Journal of Korean Tunnelling and Underground Space Association
• /
• v.15 no.3
• /
• pp.345-356
• /
• 2013

Numerical simulation on rock cutting by a TBM disc cutter was carried out using SPH (Smoothed Particle Hydrodynamics) code. AUTODYN3D, a commercial software program based on finite element method, was used in this study. The three-dimensional geometry of a disc cutter and a rock specimen were modeled by Lagrange and SPH code respectively. The numerical simulation was carried out for Hwangdeung granite for 10 different cutting conditions. The results of the numerical simulation, i.e. the relation between cutter force and failure behavior, had a good agreement with those from LCM test. The cutter forces measured in the numerical simulation had 10% deviation from the LCM test results. Moreover, the optimum cutter spacing was almost identical with the experimental results. These results indicate that SPH code can be successfully used had applicability for simulation on rock cutting by a TBM disc cutter. However, further study on Lagrange-SPH coupled modelling would be necessary to reduce the computation time.

• Journal of Biosystems Engineering
• /
• v.39 no.4
• /
• pp.310-317
• /
• 2014

Purpose: This study was conducted to analyze the air flow characteristics in a plant factory with different inlet and outlet locations using computational fluid dynamics (CFD). Methods: In this study, the flow was assumed to be a steady-state, incompressible, and three-dimensional turbulent flow. A realizable k-${\varepsilon}$ turbulent model was applied to show more reasonable results than the standard model. A CFD software was used to perform the numerical simulation. For validation of the simulation model, a prototype plant factory ($5,900mm{\times}2,800mm{\times}2,400mm$) was constructed with two inlets (${\Phi}250mm$) and one outlet ($710mm{\times}290mm$), located on the top side wall. For the simulation model, the average air current speed at the inlet was $5.11m{\cdot}s^{-1}$. Five cases were simulated to predict the airflow pattern in the plant factory with different inlet and outlet locations. Results: The root mean square error of measured and simulated air current speeds was 13%. The error was attributed to the assumptions applied to mathematical modelling and to the magnitude of the air current speed measured at the inlet. However, the measured and predicted airflow distributions of the plant factory exhibited similar patterns. When the inlets were located at the center of the side wall, the average air current speed in the plant factory was increased but the spatial uniformity was lowered. In contrast, if the inlets were located on the ceiling, the average air current speed was lowered but the uniformity was improved. Conclusions: Based on the results of this study, it was concluded that the airflow pattern in the plant factory with multilayer cultivation shelves was greatly affected by the locations of the inlet and the outlet.

• Proceedings of KOSOMES biannual meeting
• /
• 2003.05a
• /
• pp.191-200
• /
• 2003

Gunsan coastal area is one of region increasing pollution problems. One of the most important factors that cause eutrophication is nutrient materials containing nitrogen and phosphorus which stem from excreation of terrestial sources. At this study, the three-dimensional numerical hydrodynamic and ecosystem model, which was developed by Institute for Resources and Environment of Japan, were applied to analyze the processes affecting the eutrophication. The residual currents, which were obtained by integrating the simulated tidal currents over 1 tidal cycle, showed the presence of a typical. Density driven currents were generated westward at surface and eastward at the bottom in Geum estuary area where the fresh waters are flowing into. The ecosystem model was calibrated with the data surveyed in the field of the study area in annual average. The simulated results of DIN were fairly good coincided with the observed values within relative error of 32.39%. correlation coefficient(r) of 0.99. In the case of DIP, the simulated results were fairly good coincided with the observed values within relative error of 24.26%, correlation coefficient (r) of 0.82. The simulations of DIN and DIP concentrations were performed using ecosystem model under the conditions of 20 ∼ 80% pollution load reductions from pollution sources. In study area, concentration of DIN and DIP were reduced to 20∼80% and under 10% in case of the 80% reduction of the input loads from fresh water respectively. But pollution loads from sediment had hardly affected DIN and DIP concentration. For the environment management of coastal areas, in case of Kunsan area, the most important pollution sources affecting eutrophication phenomenon were found to be the input loads from fresh water.

• Geophysics and Geophysical Exploration
• /
• v.10 no.1
• /
• pp.69-76
• /
• 2007

High-frequency electromagnetic (EM) wave propagation associated with borehole ground-penetrating radar (GPR) is a complicated phenomenon. To improve the understanding of the governing physical processes, we employ a finite-difference time-domain solution of Maxwell's equations in cylindrical coordinates. This approach allows us to model the full EM wavefield associated with crosshole GPR surveys. Furthermore, the use of cylindrical coordinates is computationally efficient, correctly emulates the three-dimensional geometrical spreading characteristics of the wavefield, and is an effective way to discretise explicitly small-diameter boreholes. Numerical experiments show that the existence of a water-filled borehole can give rise to a strong waveguide effect which affects the transmitted waveform, and that excitation of this waveguide effect depends on the diameter of the borehole and the length of the antenna.

• Journal of Korean Society for Geospatial Information Science
• /
• v.14 no.4 s.38
• /
• pp.61-70
• /
• 2006

This paper presents photogrammetric processing to generate digital elevation models using SPOT-5 HRG stereo images and deals with the accuracy potential of HRG (High Resolution Geometry) supermode imagery for DEM generation. After bundle adjustment was preformed for sensor modelling, digital surface models were generated through the procedures of Epipolar image resampling and image matching. The DEM extracted from HRG imagery was compared along several test sections with the the refernce DEM which was obtained from the digital topographic maps of a scale of 1 to 5000. The ratio of the zone with DEM errors less than 5m to the whole zone was 53.8%, and about 2.5m RMSE was showed when assuming that the zones larger than 5m were affected by clouds, water bodies and buildings and excluding those zones from accuracy evaluation. In addition, the three-dimensional bird's eye view model and 3D building model were producted based on the DSM which was extracted from SPOT-5 HRG data.

• Structural Engineering and Mechanics
• /
• v.58 no.6
• /
• pp.1045-1075
• /
• 2016

Shallow footings are one of the most common types of foundations used to support mid-rise buildings in high risk seismic zones. Recent findings have revealed that the dynamic interaction between the soil, foundation, and the superstructure can influence the seismic response of the building during earthquakes. Accordingly, the properties of a foundation can alter the dynamic characteristics (natural frequency and damping) of the soil-foundation-structure system. In this paper the influence that shallow foundations have on the seismic response of a mid-rise moment resisting building is investigated. For this purpose, a fifteen storey moment resisting frame sitting on shallow footings with different sizes was simulated numerically using ABAQUS software. By adopting a direct calculation method, the numerical model can perform a fully nonlinear time history dynamic analysis to realistically simulate the dynamic behaviour of soil, foundation, and structure under seismic excitations. This three-dimensional numerical model accounts for the nonlinear behaviour of the soil medium and structural elements. Infinite boundary conditions were assigned to the numerical model to simulate free field boundaries, and appropriate contact elements capable of modelling sliding and separation between the foundation and soil elements are also considered. The influence of foundation size on the natural frequency of the system and structural response spectrum was also studied. The numerical results for cases of soil-foundation-structure systems with different sized foundations and fixed base conditions (excluding soil-foundation-structure interaction) in terms of lateral deformations, inter-storey drifts, rocking, and shear force distribution of the structure were then compared. Due to natural period lengthening, there was a significant reduction in the base shears when the size of the foundation was reduced. It was concluded that the size of a shallow foundation influences the dynamic characteristics and the seismic response of the building due to interaction between the soil, foundation, and structure, and therefore design engineer should carefully consider these parameters in order to ensure a safe and cost effective seismic design.

• The Journal of Engineering Geology
• /
• v.10 no.1
• /
• pp.13-25
• /
• 2000

A three-dimensional discrete fracture network model based on probabilistic characteristics of fracture geometry and transmissivity was designed to calculate the conductivity tensor and to estimate theanisotropy of conductivity. The conductivities, $K_p$, obtained from the numerical simulation of single-holepacker test corresponded well to those from the field tests. From this, it can be concluded that thefracture network model designed in this study can represent hydraulic characteristics of in-situ fractured rock mass. Block-scale conductivities, $K_b$, estimated from the modelling of steady-state flow through the REV-scale block were ranged between the arithmetic mean and harmonic mean of theconductivity estimates from packer tests. The conductivity along north-south direction was 1.4 timesgreater than that along the east-west direction. It was concluded that the anisotropy of conductivitywas insignificant. It was also found that there was a little correlation between $K_b$ and $K_p$. This would be to that the conductivities from the packer test simulation was strongly dependent on thetransmissivity and the number of fractures within the packer test intervals.

• Journal of the Korean GEO-environmental Society
• /
• v.10 no.5
• /
• pp.59-67
• /
• 2009

Three-dimensional (3D) numerical analyses have been conducted to study the behaviour of a single pile to tunnelling. The numerical analysis has included soil slip at the pile-soil interface. In the numerical analyses the interaction between the tunnel and the pile constructed in weathered soil and rock has been analysed. The study includes the pile settlement, the relative shear displacement between the pile and the soil and the shear stresses at the interface and the axial force on the pile. In particular, the shear stress transfer mechanism at the pile-soil interface related to the tunnel advancement has been rigorously analysed. Due to changes in the relative shear displacement at the pile-soil interface during the tunnel advancement, the shear stress and the axial force distributions along the pile have been changed. Upward shear stress developed at most part of the pile (Z/L=0.0-0.8), while downward shear stress is mobilised near the pile tip (Z/L=0.8-1.0) resulting in tensile force on the pile, where Z is the pile location and L is the pile length. Some insights into the pile behaviour to tunnelling obtained from the numerical analyses will be reported and discussed.

• Steel and Composite Structures
• /
• v.1 no.4
• /
• pp.427-440
• /
• 2001

The growing use of unprotected or partially protected steelwork in buildings has caused a lively debate regarding the safety of this form of construction. A good deal of recent research has indicated that steel members have a substantial inherent ability to resist fire so that additional fire protection can be either reduced or eliminated completely. A performance based philosophy also extends the study into the effect of structural continuity and the performance of the whole structural totality. As part of the structural system, thermal expansion during the heating phase or contraction during the cooling phase in most beams is likely to be restrained by adjacent parts of the whole system or sub-frame assembly due to compartmentation. This has not been properly addressed before. This paper describes an experimental programme in which unprotected steel beams were tested under load while it is restrained between two columns and additional horizontal restraints with particular concern on the effect of catenary action in the beams when subjected to large deflection at very high temperature. This paper also presents a three-dimensional mathematical modelling, based on the finite element method, of the series of fire tests on the part-frame. The complete analysis starts with an evaluation of temperature distribution in the structure at various time levels. It is followed by a detail 3-D finite element analysis on its structural response as a result of the changing temperature distribution. The principal part of the analysis makes use of an existing finite element package FEAST. The effect of columns being fire-protected and the beam being axially restrained has been modelled adequately in terms of their thermal and structural responses. The consequence of the beam being restrained is that the axial force in the restrained beam starts as a compression, which increases gradually up to a point when the material has deteriorated to such a level that the beam deflects excessively. The axial compression force drops rapidly and changes into a tension force leading to a catenary action, which slows down the beam deflection from running away. Design engineers will be benefited with the consideration of the catenary action.

• Geomechanics and Engineering
• /
• v.12 no.2
• /
• pp.295-312
• /
• 2017

Cast in situ and grouted concrete helical piles with 150-200 mm diameter half cylindrical ribs have become an economical and effective choice in Shanghai, China for uplift piles in deep soft soils. Though this type of pile has been successful used in practice, the reinforcing mechanism and the contribution of the ribs to the total resistance is not clear, and there is no clear guideline for the design of such piles. To study the inclusion of ribs to the contribution of shear resistance, the shear behaviour between silty sand and concrete slabs with parallel ribs at different spacing and angles were tested in a large direct shear box ($600mm{\times}400mm{\times}200mm$). The front panels of the shear box are detachable to observe the soil deformation after the test. The tests were modelled with three-dimensional finite element method in ABAQUS. It was found that, passive zones can be developed ahead of the ribs to form undulated failure surfaces. The shear resistance and failure mode are affected by the ratio of rib spacing to rib diameter. Based on the shape and continuity of the failure zones at the interface, the failure modes at the interface can be classified as "punching", "local" or "general" shear failure respectively. With the inclusion of the ribs, the pull out resistance can increase up to 17%. The optimum rib spacing to rib diameter ratio was found to be around 7 based on the observed experimental results and the numerical modelling.