• 대한기계학회논문집A
• /
• 제20권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.

• 한국철도학회:학술대회논문집
• /
• 한국철도학회 2004년도 추계학술대회 논문집
• /
• 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
• /
• 제15권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.

• 조명전기설비학회논문지
• /
• 제19권4호
• /
• pp.108-116
• /
• 2005

1계 미분항이 포함되는 미분방정식의 수치해를 구하고자 할 때 중앙차분을 사용한 유한차분법이나 Galerkin법을 사용한 유한요소법은 그 해가 매우 불안하여 요소분할을 세밀하게 하여야만 해를 얻을 수 있다. 이러한 해의 불안 정성이 일어나는 이유는 대류항의 크기가 커질수록 후류에서의 경계조건이 해의 급격한 변화를 요구하는데 수치해가 급격한 변화에 적응하지 못하기 때문이다. 이러한 문제를 해결하기 위해 1970년대부터 upwind법이 개발되어 왔다. 본 논문은 1계 미분항이 표현되는 속도기전력이 발생하는 전자계 문제를 유한요소법을 이용하여 해석할 때 발생하는 해의 진동 문제를 해결하기 위해 Heinrich에 의해 제안된 upwind법을 적용하였다.

• Journal of Advanced Marine Engineering and Technology
• /
• 제32권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
• /
• 제26권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
• /
• 제21권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.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2007년도 정기 학술대회 논문집
• /
• 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
• /
• 제13권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.

• 대한전기학회논문지
• /
• 제36권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.