• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.3
• /
• pp.933-944
• /
• 1996

In this paper, two new techniques, the pattern filling and the refined flow field regeneration, based on the finite element method and Eulerian mesh advancement approach have been developed to analyze incompressible viscous flow with free surfaces. The gorerning equation for flow analysis is Navier-Stokes equation including inertia and gravity effects. The penalty and Newton-Raphson methods are used effectively for finite element formulation. The flow front surface and the volume inflow rate are calculated using the pattern filling technique to select an adequate pattern among five filling patterns at each quadrilateral control volume. By the refined flow field regeneration technique, the new flow field which renders better prediction in flow surface shape is generated and the velocity field at the flow front part is calculated more exactly. Using the new thchniques to be developed, the dam-breaking problem has been analyzed to predict flow phenomenon of fluid and the predicted front positions versus time have been compared with the reported experimental result.

• Proceedings of the KSR Conference
• /
• 2004.10a
• /
• pp.1089-1094
• /
• 2004

Numerical methods to estimate behaviors of jointed rock mass can be roughly divided into two method : discontinuous model and continuum model. Generally, distinct element method (DEM) is applied in discontinuous model, and finite element method (FEM) or finite difference method (FDM) is utilized in continuum model. To predict a behavior of discontinuous model by DEM, it is essential to understand characteristics of joints developed in rock mass through field tests. However, results of field tests can not provide full information about rock mass because field tests is conducted in limited area. In this paper, discontinuous analysis by UDEC and continuous analysis by FLAC is utilized to estimate a behavior of a tunnel in jointed rock mass. For including discontinuous analysis in continuous analysis, joints in rock mass is considered by reducing rock mass properties obtained by RMR and decreasing shear strength of rock mass. By comparing and revising two analysis results, analysis results similar with practical behavior of a tunnel can be induced and appropriate support system is decided.

• Journal of Mechanical Science and Technology
• /
• v.15 no.11
• /
• pp.1499-1506
• /
• 2001

Various transition elements are used in general for the effective finite element analysis of complicated mechanical structures. In this paper, a solid-to-beam transition finite element, which can b e used for connecting a C1-continuity beam element to a continuum solid element, is proposed. The shape functions of the transition finite element are derived to meet the compatibility condition, and a transition element equation is formulated by the conventional finite element procedure. In order to show the effectiveness and convergence characteristics of the proposed transition element, numerical tests are performed for various examples. As a result of this study, following conclusions are obtained. (1) The proposed transition element, which meets the compatibility of the primary variables, exhibits excellent accuracy. (2) In case of using the proposed transition element, the number of nodes in the finite element model may be considerably reduced and the model construction becomes more convenient. (3) This formulation method can be applied to the usage of higher order elements.

• Journal of the Korean Institute of Illuminating and Electrical Installation Engineers
• /
• v.19 no.4
• /
• pp.108-116
• /
• 2005

This paper presents finite element analyses solution in the travelling magnetic field problem. The travelling magnetic field problem is subject to convective-diffusion equation. Therefore, the solution derived from Galerkin-FEM with linear interpolation function may oscillate between the adjacent nodes. A simple model with Dirichlet, Neumann and Periodic boundary condition respectively, have been analyzed to investigate stabilities of solutions. It is concluded that the solution of Galerkin-FEM may oscillate according to boundary condition and element type, but that of Upwind-FFM is stable regardless boundary condition.

• Journal of Advanced Marine Engineering and Technology
• /
• v.32 no.5
• /
• pp.759-767
• /
• 2008

Bending deformation of an ionic polymer actuator(IPA) on applied low electric field across its thickness is dominated by electroosmosis of hydrated ions and self-diffusion of free water molecules. In the study by Popovic et al., two processes are assumed to occur sequentially in the way that fast electroosmosis is followed by self-diffusion and finite element formulation for the basic field equations are proposed. However the motions of hydrated ions and water molecules occur at the same time. In this study, those two processes are considered simultaneously and finite element formulation is conducted for the basic field equations governing electrochemical response of an IPA. Some numerical studies for IPA are carried out in order to show the validity of the present formulation.

• Earthquakes and Structures
• /
• v.26 no.2
• /
• pp.163-174
• /
• 2024

Dynamic equilibrium equations for finite element analysis were derived for the free field one-dimensional shear wave propagation through the horizontally layered soil deposits with the elastic half-space. We expressed Rayleigh's viscous damping consisting of mass and stiffness proportional terms. We considered two cases where damping matrices are defined in the total and relative displacement fields. Two forms of equilibrium equations are presented; one in terms of total motions and the other in terms of relative motions. To evaluate the performance of new equilibrium equations, we conducted two sets of site response analyses and directly compared them with the exact closed-form frequency domain solution. Results show that the base shear force as earthquake load represents the simpler form of equilibrium equation to be used for the finite element method. Conventional finite element procedure using base acceleration as earthquake load predicts exact solution reasonably well even in soil deposits with unrealistically high damping.

• Journal of Welding and Joining
• /
• v.21 no.6
• /
• pp.77-83
• /
• 2003

The finite element analysis was performed for the cracks existing in residual stress fields in order to investigate the effects of configuration of residual stress distribution to the fatigue crack opening behaviour. And the variation of stress distributions adjacent to the crack caused by uploading was examined. The finite element model with contact elements for the crack plane and plane stress elements for the base material and the analytical method based on the superposition principle to estimate crack opening behaviour and the stress distribution adjacent to the crack subjected to uploading were used. The results of the analysis showed that crack opening behaviors and variations of stress distribution caused by uploading were changed depending on the configuration of residual stress distribution. When the crack existed in the region of compressive residual stress and the configuration of compressive residual stress distribution were inclined, a partial crack opening just behind of a crack tip occurred during uploading. Based on the above results, it was clarified that the crack opening behaviour in the residual stress field could be predicted accurately by the finite element analysis using these analytical method and model.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2007.04a
• /
• pp.789-794
• /
• 2007

Recently, the necessity of efficient and exact method to analyze structures is increasing with the importance of the seismic analysis. But the finite element method used in many field do not give the exact solution unless the length of the element is very short enough to represent the deformation of the element. Because the amount of computer calculation increase with the increasing of the number of degree of freedoms, the finite element method for the exact dynamic analysis of structures would not be efficient. To solve these problems, spectral clement method combined spectral method using the principle of wave mechanics and finite element method for the analysis of discrete models is applied to evaluate the behavior of the spatial structures. As a result of analysis. it becomes clear that the spectral element method is faster and more exact than the finite clement method.

• Journal of Power Electronics
• /
• v.13 no.4
• /
• pp.552-557
• /
• 2013

This paper proposed an analytical model for the distributed magnetic field analysis of interior permanent magnet-type blush-less direct current motors under the stator-turn fault condition using the winding function theory. Stator-turn faults cause significant changes in electric and magnetic characteristic. Therefore, many studies on stator-turn faults have been performed by simulation of the finite element method because of its non-linear characteristic. However, this is difficult to apply to on-line fault detection systems because the processing time of the finite element method is very long. Fault-tolerant control systems require diagnostic methods that have simple processing systems and can produce accurate information. Thus analytical modeling of a stator-turn fault has been performed using the winding function theory, and the distributed magnetic characteristics have been analyzed under the fault condition. The proposed analytical model was verified using the finite element method.

• The Transactions of the Korean Institute of Electrical Engineers
• /
• v.36 no.4
• /
• pp.252-258
• /
• 1987

Time periodic finite element solutions for sinusoidally excited electromagnetic field problems in moving media are presented. Solutions by the Galerkin method contain spurious oscillations when grid Peclet number is more than one. To suppress these oscillations an upwind finite element method using two different time periodic test functions is introduced. One is multiplied to second and first-order space derivative terma and the other to the time derivative term. Test functions are obtained from trial functions by adding or subtracting quadratic bias functions with appropriate scaling factors. Phase differences are considered between trial functions and bias functions. For simple interpretations of the phase differences, complex scaling factors are used. The proposed method is developed to give nodally exact solutions for uniform grid spacing in one dimensional problems. Based on the one dimensional results, a two dimensional upwinding scheme is also derived.