• Journal of Korean Tunnelling and Underground Space Association
• /
• v.8 no.2
• /
• pp.115-128
• /
• 2006

The structural anisotropy and heterogeneity of rock mass, caused by discontinuities and weak zones, have a great influence on the deformation behavior of a tunnel. A tunnel construction in these complex ground conditions is very difficult. No matter how excellent a geological investigation is, local uncertainties of rock mass conditions still remain. Under these uncertain circumstances, an accurate forecast of the ground conditions ahead of the advancing tunnel face is indispensable to a safe and economic tunnel construction. This paper presents the effect of anisotropy and heterogeneity of the rock masses to be excavated by numerical analysis. The influences of distance from weak zone, the size or dimension, the different stiffness and the orientation of weak zones are analysed by 3-D finite element analysis. By analysing these numerical results, the tunnel behavior due to excavation can be well understood and the prediction of rock mass condition ahead of the tunnel face can be possible.

• Journal of the Korean GEO-environmental Society
• /
• v.4 no.3
• /
• pp.79-88
• /
• 2003

This paper presents a numerical study of the spatial behavior of a linear absorbing solute in a heterogeneous porous medium. The spatially correlated log-normal hydraulic conductivity field is generated in a given two-dimensional domain by using the geostatistical method (Turning Bands algorithm). The velocity vector field is calculated by applying the two-dimensional saturated groundwater flow equation to the Galerkin finite element method. The simulation of solute transport is carried out by using the random walk particle tracking model with CD(constant displacement) scheme in which the time interval is automatically adjusted. In this study, the spatial behavior of a solute is analyzed by the longitudinal center-of-mass displacement, longitudinal spatial spread moment and longitudinal plume skewness.

• Journal of the Korean Society for Precision Engineering
• /
• v.30 no.4
• /
• pp.397-402
• /
• 2013

In this study, we carried out the each lines of section, using GC (green silicon carbide) whetstone, the SCM415 material which separated by after and before heat treatments process, in 3+2 axis machining centers for integrated grinding after cutting end mill works, the spindle speed 8000 rpm and feed rate 150 mm/min. For the analysis of the centerline average roughness (Ra), we measured by 10 steps stages. Using Finite element analysis, we found the result of the load analysis effect of the assembly parts, when applied the 11 kg's load on both side of the ATC (Automatic tool change) arm. The result is as follows. For the centerline average roughness (Ra) in the non-heat treatment work pieces, are appeared the most favorable in the tenth section are $0.510{\mu}m$, that were shown in the near the straight line section which is the smallest deformation of curve. In addition, the bad surface roughness appears on the path is to long by changing angle, the more inclined depth of cut, because the chip discharging is not smoothly.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.22 no.1
• /
• pp.47-57
• /
• 2010

Subgrid scale stabilization method is applied to Woo and Liu(2004)'s extended Boussinesq FEM numerical model to eliminate the 2dx wiggles. In order to optimize the computational efficiency, Hessian operator is introduced and the matrix of velocity vector is combined to one matrix for solving matrix equations. The mass lumping technique is also applied to the matrix equations of auxiliary variables. The newly developed code is applied to simulate Vincent and Briggs(1989)' wave transformation experiments and the results show that the numerical solution is almost wiggle-free and it matches very well with experimental data. Due to improvement of computational efficiency and wiggle reduction, it is plausible to apply this model to a realistic problem such as harbor oscillation problems.

• Journal of the Korean Society for Railway
• /
• v.12 no.3
• /
• pp.437-442
• /
• 2009

In the case of a Linear Induction Motor (LIM), numerical analysis method like Finite Element Method (FEM) has been mainly used to analyze the travelling magnetic field problem which includes the velocity-induced electromotive force. If the problem including the velocity-induced electromotive force is analyzed by FEM using the Galerkin method, the solution can be oscillated according to the Peclet Number, which is determined by conductivity, permeability, moving velocity and size of mesh. Consequently, the accuracy of the solution can be low and the vortexes of flux can be occurred at the secondary back-iron. These vortexes of the flux occurred at the secondary back-iron does not exist physically, but it can be occurred in the analysis. In this case, the vortexes of the flux can be generally removed by using Up-Wind method which is impossible to apply a conventional S/W tool (Maxwell 2D). Therefore, in this paper, authors examined the vortexes of the flux occurred at the secondary back-iron of the LIM according to variations of the Peclet Number, and analyzed whether these vortexes of the flux affect on the dynamic force characteristics of the LIM or not.

• The Transactions of The Korean Institute of Electrical Engineers
• /
• v.63 no.10
• /
• pp.1377-1383
• /
• 2014

This paper deals with eddy current loss of magnetic coupling with radial permanent magnet (PM) using analytical method such as a space harmonic method. Superposition of two kinds analysis model is used to analyze eddy current loss induced in inner PM and outer PM of magnetic coupling. When the eddy current is induced, the environmental temperature increases, and the permanent magnet(PM) characteristics are degraded because the performance of PM is greatly influenced by temperature rise. Hence, the calculation of eddy current loss becomes an important factor in the magnetic coupling. In order to analyze eddy current loss, first, on the basis of the magnetic vector potential and two-dimensional(2-D) polar-coordinate system, the magnetic field solutions of the radial magnetized PM are obtained. And we obtain the analytical solutions for the eddy current density produced by permanent magnet. Lastly, analytical solutions for eddy current loss are derived by using equivalent, electrical resistance calculated from magnet volume and analytical solution for eddy current density. This analytical results are validated by comparing with the 2-D finite element analysis (FEA).

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.3
• /
• pp.155-162
• /
• 2014

Using an isogeometric approach, a continuum-based shape design optimization method is developed for steady state power flow problems at high frequencies. In case the isogeometric method is employed to the shape design optimization, the NURBS basis functions used in CAD geometric modeling are directly utilized to embed the exact geometry into the computational framework so that the design parameterization for shape optimization is much easier than that in the finite element method and consequently provides the enhanced smoothness of design perturbations. Thus, exact geometric models can be used in both the response and the shape sensitivity analyses, where normal vector and curvature are continuous over the whole design space so that enhanced shape sensitivity can be expected. Through numerical examples, the developed isogeometric sensitivity is compared with finite difference one to provide excellent agreement. Also, it turns out that the proposed method works very well in the shape optimization problems.

• The Journal of the Acoustical Society of Korea
• /
• v.34 no.5
• /
• pp.335-342
• /
• 2015

Phononic crystals are composite materials consisting of a periodic arrangement of scattering inclusions in a host material. One of the most important properties of phononic crystals is the existence of band gaps, i.e., ranges of frequencies at which acoustic waves cannot propagate through the structure. The present study aims to investigate theoretically and experimentally the acoustic band structures in two-dimensional (2D) phononic crystals consisting of periodic square arrays of stainless steel solid cylinders with a diameter of 1 mm and a lattice constant of 1.5 mm in water. The theoretical dispersion relation that depicts the relationship between the frequency and the wave vector was calculated along the ${\Gamma}X$ direction of the first Brillouin zone using the finite element method to predict the band structures in the 2D phononic crystals. The transmission and the reflection coefficients were measured in the 2D phononic crystals with 1, 3, 5, 7, and 9 layers of stainless steel cylinders stacked in the perpendicular direction to propagation at normal incidence. The theoretical dispersion relation exhibited five band gaps at frequencies below 2 MHz, the first gap appearing around a frequency of 0.5 MHz. The location and the width of the band gaps experimentally observed in the transmission and the reflection coefficients appeared to coincide well with those determined from the theoretical dispersion relation.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.18 no.1
• /
• pp.688-695
• /
• 2017

This paper deals with characteristics analysis of a permanent magnet (PM) linear generator using analytical methods for wave energy harvesting. The wave energy is carried out from the movement of a yo-yo system. A linear generator using permanent magnets to generate a magnetic force itself does not require a separate power supply and has the advantage of simple maintenance. In addition to the use of a rare earth, a permanent magnet having a high-energy density can be miniaturized and lightweight, and can obtain high energy-conversion efficiency. We derived magnetic field solutions produced by the permanent magnet and armature reaction based on 2D polar coordinates and magnetic vector potential. Induced voltage is obtained via arbitrary sinusoidal input. In addition, electrical parameters are obtained, such as back-EMF constant, resistance, and self- and mutual-winding inductances. The space harmonic method used in this paper is confirmed by comparing it with finite element method (FEM) results. These facilitate the characterization of the PM-type linear generator and provide a basis for comparative studies, design optimization, and machine dynamic modeling.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.44 no.1
• /
• pp.80-87
• /
• 2016

In this work, we conduct a thermal pointing error analysis of the observation satellites considering seasonal and daily temperature variation with interpolated temperature based on prescribed average temperature (PAT) method. Maximum 200 degree temperature excursion is applied to the observation satellites during on-orbit operation, which cause the line of sight (LOS) to deviate from the designated pointing direction due to thermo-elastic deformation. To predict and adjust such deviation, the thermo-elastic deformation analysis with a fine structural finite element model is accomplished with interpolated thermal maps calculated from the results of on-station thermal analysis with a coarse thermal model. After verifying the interpolated temperatures by PAT with two benchmark problems, we evaluate the thermal pointing error.