• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.4
• /
• pp.395-405
• /
• 2004

The filling pattern and an adaptive grid refinement based on the finite element method and Eulerian mesh advancement approach have been developed to analyze incompressible transient viscous flow with free surfaces. The governing equation for flow analysis is Navier-Stokes equation including inertia and gravity effects. The mixed FE formulation and predictor-corrector method are used effectively for unsteady numerical simulation. The flow front surface and the volume inflow rate are calculated using the filling pattern technique to select an adequate pattern among four filling patterns at each triangular control volume. By adaptive grid refinement, the new flow field that renders better prediction in flow surface shape is generated and the velocity field at the flow front part is calculated more exactly. In this domain the elements in the surface region are made finer than those in the remaining regions for more efficient computation. Using the proposed numerical technique, the collapse of a water dam has been analyzed to predict flow phenomenon of fluid and the predicted front positions with respect to time have been compared with the reported experimental results.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.11
• /
• pp.1348-1358
• /
• 2004

The filling pattern and an adaptive grid refinement based on the finite element method and Eulerian mesh advancement approach have been developed to analyze incompressible transient viscous flow with free surfaces. The governing equation fur flow analysis is Navier-Stokes equation including inertia and gravity effects. The mixed FE formulation and predictor-corrector method are used effectively for unsteady numerical simulation. The flow front surface and the volume inflow rate are calculated using the filling pattern technique to select an adequate pattern among seven filling patterns at each tetrahedral control volume. By adaptive grid refinement, the new flow field that renders better prediction in flow surface shape is generated and the velocity field at the flow front part is calculated more exactly. In this domain the elements in the surface region are made finer than those in the remaining regions for more efficient computation. The collapse of a water dam and the filling of a fluidity spiral have been analyzed. The numerical results have been in good agreement with the experimental results and the efficiency of the adaptive grid refinement and filling pattern techniques have been verified.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.6
• /
• pp.703-714
• /
• 2010

This paper proposes an approach for staged finite element modeling with coupled seepage and stress analysis. The stage modeling is based on the predefined inter-relationship between the base model and the unit stage models. A unit stage constitutes a complete finite element model, of which the geometries and attributes are subject to changes from stage to stage. The seepage analysis precedes the mechanical stress analysis at every stage. Division of the wet and dry zone and the pore pressures are evaluated from the seepage analysis and used in determining input data for the stress analysis. The results of the stress analysis may also be associated with the pore water pressures. For consolidation analysis, the pore pressure and the displacement variables are mixed in a coupled matrix equation. The time marching solution produces the dissipation of excess pore pressure and variation of stresses with passage of time. For undrained analysis, the excess pore pressures are computed from the stress increment due to loading applied in the unit stage and are used in revising the hydraulic head. The solution results of a unit stage are inherited and accumulated to the subsequent stages through the relationship of the base model and the individual unit stages. Implementation of the proposed approach is outlined on the basis of the core procedures, and numerical examples are presented for demonstration of its application.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.28 no.2
• /
• pp.123-130
• /
• 2015

This paper presents a new formulation for transient simulations of microwave propagation in heterogeneous unbounded domains. In particular, perfectly-matched-layers(PMLs) are introduced to allow for wave absorption at artificial boundaries used to truncate the infinite extent of the physical domains. The development of the electromagnetic PML targets the application to engineering mechanics problems such as structural health monitoring and inverse medium problems. To formulate the PML for plane electromagnetic waves, a complex coordinate transformation is introduced to Maxwell's equations in the frequency-domain. Then the PML-endowed partial differential equations(PDEs) for transient electromagnetic waves are recovered by the application of the inverse Fourier transform to the frequency-domain equations. A mixed finite element method is employed to solve the time-domain PDEs for electric and magnetic fields in the PML-truncated domain. Numerical results are presented for plane microwaves propagating through concrete structures, and the accuracy of solutions is investigated by a series of error analyses.

• Journal of the Korean GEO-environmental Society
• /
• v.18 no.2
• /
• pp.53-58
• /
• 2017

In this paper, we propose a new method to debond between mixed finite elements for porous media in ABAQUS (2014). ABAQUS just provides debonding algorithm for the u-p model using cohesive elements in standard version. However, this approach has a drawback that it is hard to simulate complex debonding problems like element separation, rigid body motion, and contact between separated elements in standard version. ABAQUS-explicit can resolve these complex problems, but cohesive elements for the u-p model cannot be applied. We introduce a new algorithm for debonding for porous media instead of using cohesive elements. In this method, subroutines VUMAT to apply constitutive models and VDISP to separate elements in ABAQUS are used to simulate debonding problems. In addition, a simple 2-D example is demonstrated in the ABAQUS-explicit solver.

• Composites Research
• /
• v.12 no.3
• /
• pp.66-74
• /
• 1999

A modified mixed mode bending test was performed to investigate the mixed mode delamination for carbon/epoxy composite material. Various mixed mode ratios could be produced by changing the applied load position on the loading lever and the bending load position on the specimen. The modified mixed mode bending test was analyzed to obtain strain energy release rates using beam theory, compliance method and finite element method, This results were in good agreement with the experimental result, which confirmed the validity of this test.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.9
• /
• pp.1327-1334
• /
• 2003

Flow characteristics of the fluid laden with many particles in the two-dimensional channel are investigated using the Navier-Stokes equations coupled with the equation of motion of particles by direct numerical simulation. A four-step fractional step method with Crank-Nicolson scheme and ALE technique is used for P2P1 mixed finite element method. The motion and distribution of particles in the fluid is virtually described as a result of direct numerical simulation and the increase of viscosity is compared with theoretical equations. The effect of channel height on the relative viscosity and the tubular pinch effect are discussed.

• Structural Engineering and Mechanics
• /
• v.36 no.5
• /
• pp.575-598
• /
• 2010

The shape and size of the plastic zone around the crack tip are analyzed under pure mode I, pure mode II and mixed mode (I+II) loading for small scale yielding and for both plane stress and plane strain conditions. A new analytical formulation is presented to determine the radius of the plastic zone in a non-dimensional form. In particular, the effect of T-stress on the plastic zone around the crack tip is studied. The results of this investigation indicate that the stress field with a T-stress always yields a larger plastic zone than the field without a T-stress. It is found that under predominantly mode I loading, the effect of a negative T-stress on the size of the plastic zone is more dramatic than a positive T-stress. However, when mode II portion of loading is dominating the effect of both positive and negative T-stresses on the size of the plastic zone is almost equal. For validating the analytical results, several finite element analyses were performed. It is shown that the results obtained by the proposed analytical formulation are in very good agreements with those obtained from the finite element analyses.

• Journal of the Korean GEO-environmental Society
• /
• v.11 no.7
• /
• pp.5-14
• /
• 2010

The paper processed settlement analysis using Finite Elements Method(FEM). Because Stress Distribution Ratio has to be decreased, for settlement analysis of soft clay deposit improved by sand compaction piles(SCP). Back analysis was processed comparing the measured settlements of laboratory model tests and finite element analysis where the SCP treated area was assumed as mixed ground with clay deposit rather than being a composite ground. The paper proposes a methodology which employs a compression index($C_c$) for settlement analysis of soft clay deposit improved by sand compaction piles from the back analysis. This approach is applied to a field measurement case(A revetment founded on the SCP improved clay deposit with the replacement ratio of 45%).

• Proceedings of the KIEE Conference
• /
• 1998.07a
• /
• pp.70-72
• /
• 1998

This paper discusses the development of a new optimization algorithm for DC motor design. In principle, the new algorithm utilizes a mixed method that consists of genetic algorithms in conjunction with direct search method. The genetic algorithms are used for locating the global optimum region while the direct search method is used to achieve objective function convergence. In order to validate the effectiveness, the new algorithm has been applied to an actual DC motor. Field and torque characteristics of the DC motor are computed using finite element method and the principle of virtual work, respectively.