• Proceedings of the Korean Geotechical Society Conference
• /
• 2001.03a
• /
• pp.97-104
• /
• 2001

The magnitude and the orientation of in-situ stresses contribute to ground displacement and stresses in the field of underground space. This paper investigates in-situ stresses at various depth on the basis of 392 data which were determined by over-coring and hydro-fracturing test methods in the Korea peninsula. The result shows that in-situ stress distribution are more or less non-uniform through the Granite and Gneiss sub-area, and that the K-value in the Volcanic sub-area are below 1 at the deep depth. Also, the result of three dimensional numerical analyses of tunnel shows that the direction and magnitude of displacement around tunnel are much effected by the stress difference between the maximum and the minimum horizontal stress.

• Journal of Advanced Marine Engineering and Technology
• /
• v.20 no.5
• /
• pp.45-57
• /
• 1996

Flow characteristics and contamination behavior in a three-dimensional room with the desk-type obstacle have been investigated numerically by the k-${\varepsilon}$ two equation turbulence model. The room model has one supply opening on the ceiling and two exhaust openings on the side walls. Thirty six sets of calculation have been performed for twelve contamination sources of $1{\times}10^{-4}kg_c/(m^3{\cdot}s)$ strength at different inlet velocities(0.1, 1, 10m/s). This study can conclude that the source points of contaminant are located near the obstacle edge of Z-axis, at which the maximum contaminant diffusion fields are occured.

• Advances in concrete construction
• /
• v.3 no.2
• /
• pp.127-143
• /
• 2015

Steel fibre reinforced concrete (SFRC) is an anisotropic material due to the random orientation of the fibres within the cement matrix. Fibres under different inclination angles provide different strength contribution of a given crack width. For that the pull-out response of inclined fibres is of great importance to understand SFRC behaviour, particularly in the case of fibres with hooked ends, which are the most widely used. The paper focuses on the numerical modelling of the pull-out response of this kind of fibres from high-strength cementitious matrix in order to study the effects of different inclination angles of the fibres to the load-displacement pull-out curves. The pull-out of the fibres is studied by means of accurate three-dimensional finite element models, which take into account the nonlinearities that are present in the physical model, such as the nonlinear bonding between the fibre and the matrix in the early stages of the loading, the unilateral contact between the fibre and the matrix, the friction at the contact areas, the plastification of the steel fibre and the plastification and cracking of the cementitious matrix. The bonding properties of the fibre-matrix interface considered in the numerical model are based on experimental results of pull-out tests on straight fibres.

• Computers and Concrete
• /
• v.12 no.3
• /
• pp.261-283
• /
• 2013

This paper aims to contribute to the three-dimensional generalization of numerical prediction of crack propagation through the formulation of finite elements with embedded discontinuities. The analysis of crack propagation in two-dimensional problems yields lines of discontinuity that can be tracked in a relatively simple way through the sequential construction of straight line segments oriented according to the direction of failure within each finite element in the solid. In three-dimensional analysis, the construction of the discontinuity path is more complex because it requires the creation of plane surfaces within each element, which must be continuous between the elements. In the method proposed by Chaves (2003) the crack is determined by solving a problem analogous to the heat conduction problem, established from local failure orientations, based on the stress state of the mechanical problem. To minimize the computational effort, in this paper a new strategy is proposed whereby the analysis for tracking the discontinuity path is restricted to the domain formed by some elements near the crack surface that develops along the loading process. The proposed methodology is validated by performing three-dimensional analyses of basic problems of experimental fractures and comparing their results with those reported in the literature.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.29 no.6
• /
• pp.335-349
• /
• 2017

Submerged breakwaters installed under the water surface are a representative coastal structure to prevent coastal erosion, and various types of submerged breakwaters have been proposed and discussed so far. Generally, submerged breakwaters make the complex wave fields due to abrupt change in water depth at the crown of the breakwater. In this study, wave heights and mean water level formed around a breakwater are examined numerically for three-dimensional permeable submerged breakwaters. OLAFOAM, CFD open source code, is applied in the numerical analysis, and the comparisons are made with available experimental results on the permeable upright wall and the impermeable submerged breakwater to verify its applicability to the three-dimensional numerical analysis. Based on the applicability of OLAFOAM numerical code, the wave height and mean water level distribution formed around the permeable submerged breakwaters are investigated under the formation condition of salient. The numerical results show that as the gap width between breakwaters decreases, the wave height in the center of the gap increases, while it decreases behind the gap, and the installing position of the breakwater from the shoreline has little influence on the change of the wave height. Furthermore, it is found that the decrease of the mean water level near the gap between breakwaters increases with decreasing of the gap width.

• Geomechanics and Engineering
• /
• v.8 no.1
• /
• pp.97-111
• /
• 2015

The kinematical approach of limit analysis is used to estimate the three dimensional stability analysis of rock slopes with nonlinear Hoek-Brown criterion under earthquake forces. The generalized tangential technique is introduced, which makes limit analysis apply to rock slope problem possible. This technique formulates the three dimensional stability problem as a classical nonlinear programming problem. A nonlinear programming algorithm is coded to search for the least upper bound solution. To prove the validity of the present approach, static stability factors are compared with the previous solutions, using a linear failure criterion. Three dimensional seismic and static stability factors are calculated for rock slopes. Numerical results of indicate that the factors increase with the ratio of slope width and height, and are presented for practical use in rock engineering.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1999.05a
• /
• pp.5-8
• /
• 1999

The present work is concerned with three-dimensional finite element analysis of the hollow section extrusion process using a porthole die. For economic computation, mismatching refinement, an efficient domain decomposition method with different mesh density for each subdomain, is implemented. The proposed method improves the computational efficiency significantly, especially fur the three-dimensional analysis of extrusion problems. As a numerical example, extrusion of the underframe part of a railroad vehicle are analyzed. For three-dimensional mesh generation of a complicated shape with hexahedral elements, a modified grid-based approach with the surface element layer is utilized. The analysis results are then successfully reflected on the industrial porthole die design.

• Nuclear Engineering and Technology
• /
• v.46 no.5
• /
• pp.641-654
• /
• 2014

A new computer code, named CORONA (Core Reliable Optimization and thermo-fluid Network Analysis), was developed for the core thermo-fluid analysis of a prismatic gas cooled reactor. The CORONA code is targeted for whole-core thermo-fluid analysis of a prismatic gas cooled reactor, with fast computation and reasonable accuracy. In order to achieve this target, the development of CORONA focused on (1) an efficient numerical method, (2) efficient grid generation, and (3) parallel computation. The key idea for the efficient numerical method of CORONA is to solve a three-dimensional solid heat conduction equation combined with one-dimensional fluid flow network equations. The typical difficulties in generating computational grids for a whole core analysis were overcome by using a basic unit cell concept. A fast calculation was finally achieved by a block-wise parallel computation method. The objective of the present paper is to summarize the motivation and strategy, numerical approaches, verification and validation, parallel computation, and perspective of the CORONA code.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.6 no.2
• /
• pp.297-306
• /
• 2014

The paper presents the numerical analysis of wave run-up over rubble-mound breakwaters covered by antifer units using a technique integrating Computer-Aided Design (CAD) and Computational Fluid Dynamics (CFD) software. Direct application of Navier-Stokes equations within armour blocks, is used to provide a more reliable approach to simulate wave run-up over breakwaters. A well-tested Reynolds-averaged Navier-Stokes (RANS) Volume of Fluid (VOF) code (Flow-3D) was adopted for CFD computations. The computed results were compared with experimental data to check the validity of the model. Numerical results showed that the direct three dimensional (3D) simulation method can deliver accurate results for wave run-up over rubble mound breakwaters. The results showed that the placement pattern of antifer units had a great impact on values of wave run-up so that by changing the placement pattern from regular to double pyramid can reduce the wave run-up by approximately 30%. Analysis was done to investigate the influences of surface roughness, energy dissipation in the pores of the armour layer and reduced wave run-up due to inflow into the armour and stone layer.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.35 no.1
• /
• pp.54-59
• /
• 1999

High-speed electronic digital computers have enabled engineers to employ various numerical discretization techniques for solutions of complex problems. The Finite Element Method is one of the such technique. The Finite Element Method is one of the numerical analysis based on the concepts of fundamental mathematical approximation. Three dimensional plate structures used often in partition of ship, box girder and frame are analyzed by Finite Element Method. In design of structures, the static deflections, stress concentrations and dynamic deflections must be considered. However, these problem belong to geometrically nonlinear mechanical structure analysis. The analysis of each element is independent, but coupling occurs in assembly process of elements. So, to overcome such a difficulty the shell theory which includes transformation matrix and a fictitious rotational stiffness is taken into account. Also, the Mindlin's theory which is considered the effect of shear deformation is used. The Mindlin's theory is based on assumption that the normal to the midsurface before deformation is "not necessarily normal to the midsurface after deformation", and is more powerful than Kirchoff's theory in thick plate analysis. To ensure that a small number of element can represent a relatively complex form of the type which is liable to occur in real, rather than in academic problem, eight-node quadratic isoparametric elements are used. are used.